CN105574249A - 蜗杆状砂轮磨齿的磨削力模型 - Google Patents

蜗杆状砂轮磨齿的磨削力模型 Download PDF

Info

Publication number
CN105574249A
CN105574249A CN201510927729.8A CN201510927729A CN105574249A CN 105574249 A CN105574249 A CN 105574249A CN 201510927729 A CN201510927729 A CN 201510927729A CN 105574249 A CN105574249 A CN 105574249A
Authority
CN
China
Prior art keywords
grinding
rightarrow
model
gear
emery wheel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510927729.8A
Other languages
English (en)
Other versions
CN105574249B (zh
Inventor
石照耀
于渤
林家春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201510927729.8A priority Critical patent/CN105574249B/zh
Publication of CN105574249A publication Critical patent/CN105574249A/zh
Application granted granted Critical
Publication of CN105574249B publication Critical patent/CN105574249B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/09Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Finish Polishing, Edge Sharpening, And Grinding By Specific Grinding Devices (AREA)
  • Grinding-Machine Dressing And Accessory Apparatuses (AREA)

Abstract

本发明涉及蜗杆状砂轮磨齿的磨削力模型,属于齿轮机械制造领域。本发明基于G.Wener磨削力模型,提出了一种新的蜗杆状砂轮磨齿的磨削力模型。首先,根据蜗杆状砂轮磨齿过程中,接触点处珩磨轮和工件的几何形状和速度关系,建立一个近似模型;其次,以近似模型为基础建立切屑模型,计算出磨削过程中速度、磨削厚度和当量直径;最后,以近似模型为基础,建立间齿珩磨磨削力模型。这种磨削力模型适用于蜗杆砂轮磨齿和间齿珩齿加工方法。这种方法适用于计算采用蜗杆状砂轮的磨齿加工中的磨削力,包括蜗杆砂轮磨齿和间齿珩齿加工方法。这种方法将磨齿过程近似成外圆磨过程,从而简化了计算过程,填补了蜗杆状砂轮磨齿的磨削力模型的空白。

Description

蜗杆状砂轮磨齿的磨削力模型
技术领域
本发明涉及一种蜗杆状砂轮磨齿的磨削力模型,属于齿轮机械制造领域。
背景技术
磨削机理的研究可以有效的解决磨削过程中的关键问题。磨削机理的主要研究内容包括磨削参数、磨削力、磨削温度、砂轮磨损、砂轮修整、磨削液和磨削加工零件的表面完整性等。根据磨削参数和砂轮与工件的特征,可以建立磨削力模型。磨削功率则是建立在磨削力和磨削参数基础之上,功率消耗中的一部分能力转化为热量,便可建立磨削温度模型,磨削温度模型可以预测磨削烧伤现象,避免由于温度过高而产生对工件有害的磨削烧伤。通过磨削力的监测,可以有效的了解砂轮磨损情况,预测砂轮修整周期。对磨削过程建模,也可以得到表面粗糙度的经验公式。20世纪70年代G.Wener研究了材料对磨削力的影响,建立了磨削力经验模型,用于计算平面磨和外圆磨过程的磨削力。G.Wener磨削力模型为后来各种磨削过程中磨削力的研究奠定了基础。
齿轮磨削是一种复杂的磨削过程,其砂轮和工件并不是简单的圆柱或平面,同时,在磨削加工中,砂轮和工件上接触点的运动关系是实时变化的,因此很难建立其磨削力模型。BogdanW.Kruszynski等人针对Niles方法进行研究,Niles方法是Niles公司提出的一种齿轮展成磨加工方法,他们根据Niles方法中砂轮形状以及砂轮和工件的位置关系,建立切屑模型,从而建立简单的力学模型。E.Brinksmeier等人利用近似模型研究珩齿加工过程,根据实际加工过程中产生的速度建立近似模型,再针对近似模型建立磨削力模型。蜗杆状砂轮的形状更为复杂,为了更好的研究蜗杆状砂轮磨齿加工方法,有必要对蜗杆状砂轮磨齿过程的磨削力进行研究。
发明内容
为了更好的研究蜗杆状砂轮磨齿加工方法,有必要对蜗杆状砂轮磨齿过程的机理进行研究。
本发明采用的技术方案为蜗杆状砂轮磨齿的磨削力模型,首先,根据蜗杆状砂轮磨齿过程中,接触点处珩磨轮和工件的几何形状和速度关系,建立一个近似模型;其次,以近似模型为基础建立切屑模型,计算出磨削过程中速度、磨削厚度和当量直径;最后,以近似模型为基础,建立间齿珩磨磨削力模型。这种磨削力模型适用于蜗杆砂轮磨齿和间齿珩齿加工方法。
具体步骤如下:
S1.建立磨齿过程近似模型
磨齿过程中,被加工齿轮的理论齿面与砂轮表面为空间交错轴啮合,接触情况为点接触,坐标系如图1,其中Sc(Oc-Xc,Yc,Zc)和Sg(Og-Xg,Yg,Zg)分别是砂轮和被加工齿轮所在的静止坐标系,Sc和Sg分别为砂轮和被加工齿轮的坐标系名称,Oc和Og分别为砂轮和被加工齿轮的对称中心,Zc和Zg分别与砂轮和被加工齿轮的轴线重合,XcOcYc平面为砂轮齿宽中点的轴截面所在平面,XgOgYg平面为被加工齿轮齿宽中点的轴截面所在平面,l为砂轮沿着被加工齿轮轴线的移动距离。S1(O1-X1,Y1,Z1)和S2(O2-X2,Y2,Z2)分别是与砂轮和被加工齿轮固连的运动坐标系,S1和S2分别为砂轮和被加工齿轮的坐标系名称,O1和O2分别与Oc和Og重合,Z1和Z2分别与砂轮和被加工齿轮的轴线重合,为砂轮与被加工齿轮的转动角度,ro1和ro2分别为砂轮与被加工齿轮的工作节圆半径,λ是砂轮与被加工齿轮的轴交角,T是砂轮与被加工齿轮的节点。
根据啮合原理,在接触点处,被加工齿轮的齿面与砂轮表面的相对速度沿着法向的分量为零,被加工齿轮的速度vg在切平面上的投影vgt和珩磨轮的速度vc在切平面上的投影vct共线。vgt和vct是产生切削作用的速度,分别对应着磨削原理中工件速度和砂轮速度。过vgt和vct所在直线做一个平面Σ,该平面同时经过接触点法向,平面Σ分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg和Γc,这两条曲线在接触点出的曲率半径分别为r和r。将接触点处的运动类比外圆磨削,分别以曲率半径r和r为工件和砂轮的半径,建立近似模型,如图2所示。
在实际接触中,被加工齿轮与珩磨轮并非点接触,由于存在磨削余量,实际接触过程为面接触。为更符合实际情况,将类比被加工齿轮与砂轮的两个圆柱设计成中鼓形状,中鼓形状为圆弧,如图3所示。vgt和vct所在直线垂直做一个平面Σ1,该平面同时经过接触点法向,平面Σ1分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg1和Γc1,这两条曲线在接触点出的曲率半径分别为rgΓ1和rcΓ1
沿着轴向,将上述模型进行无限分割,每两个分割界面之间的部分,能够被看做外圆磨削过程,因此,利用G.Wener提出的磨削力模型,建立每一个分割部分的外圆磨削力模型,对整个磨削宽度范围内的磨削力进行积分,便能够获得当前接触点处的磨削力。G.Wener的磨削力F′n的模型为
F′n=K′n[c1]γ[vw/vs]2ε-1[ap]ε[de]1-ε
其中,vw和vs分别为被加工齿轮和砂轮的速度,ap为磨削厚度,de为当量直径,K′为单位切削面积的磨削力,c1是与磨刃密度有关的系数,γ和ε是与磨刃在砂轮圆周表面的分布有关的指数。
S2.计算磨齿过程中近似外圆磨过程的vw、vs、ap和de
(1)计算vw和vs
根据S1中所述,被加工齿轮的速度在切平面上的投影(即vw)和珩磨轮的速度在切平面上的投影(即vw)是产生切削的速度,因此需要根据间齿啮合模型推导这两个速度分量。
加工过程中砂轮转速为ωc,被加工齿轮转速为ωg,因此砂轮在Sc坐标系中的转速和被加工齿轮在Sg坐标系中的转速
ω → c c = ( 0 , 0 , ω c ) T
ω → g g = ( 0 , 0 , ω g ) T
Sc坐标系中,珩磨轮在接触点处的线速度
v → c c = ω → c c × r → 1 c = ( - y 1 c ω c , x 1 c ω c , 0 ) T
在Sg坐标系中,被加工齿轮在接触点的线速度
v → g g = ω → g g × r → 2 g = ( - y 2 g ω g , x 2 g ω g , 0 ) T
将珩磨轮在接触点的速度由Sc坐标系转化到Sg坐标系中,
v → c g = M g c v → c c = ( y 1 c ω c , - x 1 c ω c c o s λ , x 1 c ω c s i n λ )
那么接触点处的相对速度
v → c g g = v → c g - v → g g = y 1 c ω c + y 2 g ω g - x 1 c ω c c o s λ - x 2 g ω g x 1 c ω c sin λ
在Sg坐标系中,被加工齿面的法向量为n2g。在齿面上一点处的切向量有无数条,利用齿面方程对齿面参数θ2求导,得到一个切矢量都垂直的方向向量
ρ → l 2 g = n → 2 g × ρ → θ 2 g
也是齿面的一条切向量,张成一个接触点处的切平面。相互垂直,并且都是单位向量,在接触点处构成了一个新的直角坐标系St
在坐标系St中,由于沿着齿面法向的速度并不起到切削的作用,因此只研究速度在切平面内的投影。被加工齿轮的速度和珩磨轮的速度方向上的投影分别为
v → c ρ θ 2 = v → c g · ρ → θ 2 g v → c ρ l = v → c g · ρ → l 2 g
v → g ρ θ 2 = v → g g · ρ → θ 2 g v → g ρ l = v → g g · ρ → l 2 g
因此,珩磨轮的速度vc在切平面上的投影的夹角为 ψ c g = arccos | v → c ρ l | | v → c t | . 被加工齿轮的速度vg在切平面上的投影 | v → g t | = ( v → g ρ θ 2 ) 2 + ( v → g ρ l ) 2 , v → g t 的夹角为 ψ g g = arccos | v → g ρ l | | v → g t | .
(2)计算ap
两个鼓形圆柱相交部分,除去已经被磨削的部分,即为切屑近似模型,事实上这个切屑模型为相交部分的四分之一,如图4所示。
将近似模型沿轴向分割成无限小切片,那么每一层切片为一个外圆磨模型。每个小外圆磨模型中,磨削深度最大的位置是在过两个轴线的平面上。不同切削宽度位置的磨削深度分布在过两个轴线的平面与切屑近似模型相交的截面上,某磨削宽度位置bi处的磨削深度为api,几何关系如图所示。
图4中,bi为线段CcCg到线段OgΓ1OcΓ1的距离,因为线段AcBc远小于直径,这个公式简化为
b i 2 ≈ 2 r c Γ 1 A c B c
A c B c ≈ b i 2 / 2 r c Γ 1
同理线段AgBg
A g B g ≈ b i 2 / 2 r g Γ 1
综上所述,磨削厚度api
a p i = a p - A c B c - A g B g = a p - b i 2 ( 1 / 2 r c Γ 1 - 1 / 2 r g Γ 1 )
上述公式中,rcΓ1和rgΓ1是未知的,ap为最大磨削厚度。
根据S1节中分析可知,rcΓ1和rgΓ1是分别为珩磨轮和被加工齿轮的齿面在接触点处于相对速度垂直的方向上的曲率半径。
(3)计算de
在磨削力模型当中,与砂轮和工件几何尺寸相关的参数为当量直径de,无论是哪种类型的磨削过程,都有一个对应的当量直径。根据外圆磨当量直径公式通过公式看出与当量直径相关的砂轮和工件的几何尺寸为各自的半径,ds和dw分别对应着S1中建立的近似模型中接触点处的曲率半径r和r。曲率半径为曲率的倒数,因此需要求得曲率半径r和r分别对应的法曲率k和k
当近似模型中,当量直径de
d e = 2 r c Γ r g Γ r c Γ + r g Γ
实际接触时,接触区域很小,在接触区域内所有位置的当量直径变化很小,因此把接触点的当量直径近似看做接触区域的当量直径。
S3.计算磨削力
磨齿过程被类比成S1中提出的磨削近似模型,在近似模型当中,每一个垂直于轴的截面都看做外圆磨,那么根据G.Wener的模型有
F′n=K′n[c1]γ[vgt/vct]2ε-1[api]ε[de]1-ε
在同一接触点,切屑磨削宽度上不同位置的vgt、vct和de是相同的,沿着磨削宽度不同位置的磨削深度是变化的,那么接触位置法向磨削力的模型表示为,
F n = ∫ o b a F n ′ db i
将单位宽度上磨削力模型带入上述公式,并进行积分,并得到
Fn=Kn[c1]γ[vgt/vct]2ε-1[de]1-ε{[ap]ε+1-[ap-kpba 2]ε+1}
其中 k p = 1 / 2 r c Γ 1 - 1 / 2 r g Γ 1 , K n = K n ′ 2 k p ( ϵ + 1 ) .
所述的磨削力模型适用于蜗杆砂轮磨齿和间齿珩齿加工方法。
附图说明
图1空间交错轴啮合坐标系。
图2近似外圆磨模型。
图3中鼓形状近似模型。
图4切屑模型。
具体实施方式
以下结合具体加工实例对本发明进行说明:
采用蜗杆砂轮磨齿技术加工斜齿圆柱齿轮,被加工齿轮参数:法向压力角αn2=27.5°,法向模数mn2=3,齿数Z2=27,螺旋角β2=18°,右旋。蜗杆状砂轮参数:法向压力角αn1=27.5°,法向模数mn1=3,头数Z1=4,螺旋升角γ1=8°,右旋。轴交角λ=100°,中心距a=85.696。
S1.建立磨齿过程近似模型
磨齿过程中,被加工齿轮的理论齿面与砂轮表面为空间交错轴啮合,接触情况为点接触,坐标系如图1。根据啮合原理,在接触点处,齿面与砂轮表面的相对速度沿着法向的分量为零,被加工齿轮的速度vg在切平面上的投影vgt和珩磨轮的速度vc在切平面上的投影vct共线。vgt和vct是产生切削作用的速度,分别对应着磨削原理中工件速度和砂轮速度。过vgt和vct所在直线做一个平面Σ,该平面同时经过接触点法向,平面Σ分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg和Γc,这两条曲线在接触点出的曲率半径分别为r和r。将接触点处的运动类比外圆磨削,分别以曲率半径r和r为工件和砂轮的半径,建立近似模型如图2。
在实际接触中,被加工齿轮与珩磨轮并非点接触,由于存在磨削余量,实际接触过程为面接触。为了更符合实际情况,将类比被加工齿轮与砂轮的两个圆柱设计成中鼓形状,中鼓形状为圆弧,如图3所示。vgt和vct所在直线垂直做一个平面Σ1,该平面同时经过接触点法向,平面Σ1分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg1和Γc1,这两条曲线在接触点出的曲率半径分别为rgΓ1和rcΓ1
沿着轴向,将上述模型进行无限分割,每两个分割界面之间的部分,可以被看做外圆磨削过程,因此,利用G.Wener提出的磨削力模型,建立每一个分割部分的外圆磨削力模型,对整个磨削宽度范围内的磨削力进行积分,便可以获得当前接触点处的磨削力。G.Wener的磨削力模型为
F′n=K′n[c1]γ[vw/vs]2ε-1[ap]ε[de]1-ε
其中,vw和vs分别为工件和砂轮的速度,ap为磨削厚度,de为当量直径,K′为单位切削面积的磨削力,c1是与磨刃密度有关的系数,γ和ε是与磨刃在砂轮圆周表面的分布有关的指数。
S2.计算磨齿过程中近似外圆磨过程的vw、vs、ap和de
(1)计算vw和vs
根据S1中所述,被加工齿轮的速度在切平面上的投影(即vw)和珩磨轮的速度在切平面上的投影(即vw)是产生切削的速度,因此需要根据间齿啮合模型推导这两个速度分量。
加工过程中砂轮转速为ωc,被加工齿轮转速为ωg,因此
ω → c c = ( 0 , 0 , ω c ) T
ω → g g = ( 0 , 0 , ω g ) T
珩磨轮在接触点处的线速度
v → c c = ω → c c × r → 1 c = ( - y 1 c ω c , x 1 c ω c , 0 ) T
被加工齿轮在接触点的线速度
v → g g = ω → g g × r → 2 g = ( - y 2 g ω g , x 2 g ω g , 0 ) T
将珩磨轮在接触点的速度由Sc坐标系转化到Sg坐标系中,
v → c g = M g c v → c c = ( y 1 c ω c , - x 1 c ω c c o s λ , x 1 c ω c s i n λ )
那么接触点处的相对速度 v → 12 g
v → c g g = v → c g - v → g g = y 1 c ω c + y 2 g ω g - x 1 c ω c c o s λ - x 2 g ω g x 1 c ω c sin λ
在Sg坐标系中,被加工齿面的法向量为n2g。在齿面上一点处的切向量有无数条,利用齿面方程对齿面参数θ2求导,可以得到一个切矢量都垂直的方向向量
ρ → l 2 g = n → 2 g × ρ → θ 2 g
也是齿面的一条切向量,张成一个接触点处的切平面。相互垂直,并且都是单位向量,在接触点处构成了一个新的直角坐标系St
在坐标系St中,由于沿着齿面法向的速度并不起到切削的作用,因此只研究速度在切平面内的投影。被加工齿轮的速度和珩磨轮的速度方向上的投影分别为
v → c ρ θ 2 = v → c g · ρ → θ 2 g v → c ρ l = v → c g · ρ → l 2 g
v → g ρ θ 2 = v → g g · ρ → θ 2 g v → g ρ l = v → g g · ρ → l 2 g
因此,珩磨轮的速度vc在切平面上的投影的夹角为 ψ c g = arccos | v → c ρ l | | v → c t | . 被加工齿轮的速度vg在切平面上的投影 | v → g t | = ( v → g ρ θ 2 ) 2 + ( v → g ρ l ) 2 , v → g t 的夹角为 ψ g g = arccos | v → g ρ l | | v → g t | .
(2)计算ap
两个鼓形圆柱相交部分,除去已经被磨削的部分,即为切屑近似模型,事实上这个切屑模型为相交部分的四分之一,如图4所示。
将近似模型沿着轴向分割成无限小切片,那么每一层切片可看做一个外圆磨模型。每个小模型中,磨削深度最大的位置是在过两个轴线的平面上。不同切削宽度位置的磨削深度分布在过两个轴线的平面与切屑近似模型相交的截面上,某磨削宽度位置bi处的磨削深度为api,几何关系如图所示。
图中,因为AcBc远小于直径,这个公式可以简化为
b i 2 ≈ 2 r c Γ 1 A c B c
A c B c ≈ b i 2 / 2 r c Γ 1
同理,
A g B g ≈ b i 2 / 2 r g Γ 1
综上所述
a p i = a p - A c B c - A g B g = a p - b i 2 ( 1 / 2 r c Γ 1 - 1 / 2 r g Γ 1 )
上述公式中,rcΓ1和rgΓ1是未知的。
根据S1节中分析可知,rcΓ1和rgΓ1是分别为珩磨轮和被加工齿轮的齿面在接触点处于相对速度垂直的方向上的曲率半径。
(3)计算de
在磨削力模型当中,与砂轮和工件几何尺寸相关的参数为当量直径de,无论是哪种类型的磨削过程,都有一个对应的当量直径。根据外圆磨当量直径公式通过公式,可以看出与当量直径相关的砂轮和工件的几何尺寸为各自的半径,ds和dw分别对应着S1中建立的近似模型中接触点处的曲率半径r和r。曲率半径为曲率的倒数,因此需要求得曲率半径r和r分别对应的法曲率k和k
当近似模型中,当量直径为
d e = 2 r c Γ r g Γ r c Γ + r g Γ
实际接触时,接触区域很小,可以认为在接触区域内所有位置的当量直径变化很小,因此可以把接触点的当量直径近似看做接触区域的当量直径。
S3.计算磨削力
磨齿过程被类比成S1中提出的磨削近似模型,在近似模型当中,每一个垂直于轴的截面都可以看做外圆磨,那么根据G.Wener的模型有
F′n=K′n[c1]γ[vgt/vct]2ε-1[api]ε[de]1-ε
在同一接触点,切屑磨削宽度上不同位置的vgt、vct和de可认为是相同的,沿着磨削宽度不同位置的磨削深度是变化的,那么接触位置法向磨削力的模型可以表示为,
F n = ∫ o b a F n ′ db i
将单位宽度上磨削力模型带入上述公式,并进行积分,可以得到
Fn=Kn[c1]γ[vgt/vct]2ε-1[de]1-ε{[ap]ε+1-[ap-kpba 2]ε+1}
其中 k p = 1 / 2 r c Γ 1 - 1 / 2 r g Γ 1 , K n = K n ′ 2 k p ( ϵ + 1 ) .

Claims (2)

1.蜗杆状砂轮磨齿的磨削力模型,其特征在于:首先,根据蜗杆状砂轮磨齿过程中,接触点处珩磨轮和工件的几何形状和速度关系,建立一个近似模型;其次,以近似模型为基础建立切屑模型,计算出磨削过程中速度、磨削厚度和当量直径;最后,以近似模型为基础,建立间齿珩磨磨削力模型;这种磨削力模型适用于蜗杆砂轮磨齿和间齿珩齿加工方法;
具体步骤如下:
S1.建立磨齿过程近似模型
磨齿过程中,被加工齿轮的理论齿面与砂轮表面为空间交错轴啮合,接触情况为点接触,坐标系如图1,其中Sc(Oc-Xc,Yc,Zc)和Sg(Og-Xg,Yg,Zg)分别是砂轮和被加工齿轮所在的静止坐标系,Sc和Sg分别为砂轮和被加工齿轮的坐标系名称,Oc和Og分别为砂轮和被加工齿轮的对称中心,Zc和Zg分别与砂轮和被加工齿轮的轴线重合,XcOcYc平面为砂轮齿宽中点的轴截面所在平面,XgOgYg平面为被加工齿轮齿宽中点的轴截面所在平面,l为砂轮沿着被加工齿轮轴线的移动距离;S1(O1-X1,Y1,Z1)和S2(O2-X2,Y2,Z2)分别是与砂轮和被加工齿轮固连的运动坐标系,S1和S2分别为砂轮和被加工齿轮的坐标系名称,O1和O2分别与Oc和Og重合,Z1和Z2分别与砂轮和被加工齿轮的轴线重合,为砂轮与被加工齿轮的转动角度,ro1和ro2分别为砂轮与被加工齿轮的工作节圆半径,λ是砂轮与被加工齿轮的轴交角,T是砂轮与被加工齿轮的节点;
根据啮合原理,在接触点处,被加工齿轮的齿面与砂轮表面的相对速度沿着法向的分量为零,被加工齿轮的速度vg在切平面上的投影vgt和珩磨轮的速度vc在切平面上的投影vct共线;vgt和vct是产生切削作用的速度,分别对应着磨削原理中工件速度和砂轮速度;过vgt和vct所在直线做一个平面Σ,该平面同时经过接触点法向,平面Σ分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg和Γc,这两条曲线在接触点出的曲率半径分别为r和r;将接触点处的运动类比外圆磨削,分别以曲率半径r和r为工件和砂轮的半径,建立近似模型;
在实际接触中,被加工齿轮与珩磨轮并非点接触,由于存在磨削余量,实际接触过程为面接触;为更符合实际情况,将类比被加工齿轮与砂轮的两个圆柱设计成中鼓形状,中鼓形状为圆弧;vgt和vct所在直线垂直做一个平面Σ1,该平面同时经过接触点法向,平面Σ1分别于被加工齿轮的理论齿面与珩磨轮表面相交,得到两条交线Γg1和Γc1,这两条曲线在接触点出的曲率半径分别为rgΓ1和rcΓ1
沿着轴向,将上述模型进行无限分割,每两个分割界面之间的部分,能够被看做外圆磨削过程,因此,利用G.Wener提出的磨削力模型,建立每一个分割部分的外圆磨削力模型,对整个磨削宽度范围内的磨削力进行积分,便能够获得当前接触点处的磨削力;G.Wener的磨削力F′n的模型为
F′n=K′n[c1]γ[vw/vs]2ε-1[ap]ε[de]1-ε
其中,vw和vs分别为被加工齿轮和砂轮的速度,ap为磨削厚度,de为当量直径,K′为单位切削面积的磨削力,c1是与磨刃密度有关的系数,γ和ε是与磨刃在砂轮圆周表面的分布有关的指数;
S2.计算磨齿过程中近似外圆磨过程的vw、vs、ap和de
(1)计算vw和vs
根据S1中所述,被加工齿轮的速度在切平面上的投影(即vw)和珩磨轮的速度在切平面上的投影(即vw)是产生切削的速度,因此需要根据间齿啮合模型推导这两个速度分量;
加工过程中砂轮转速为ωc,被加工齿轮转速为ωg,因此砂轮在Sc坐标系中的转速和被加工齿轮在Sg坐标系中的转速
ω → c c = ( 0 , 0 , ω c ) T
ω → g g = ( 0 , 0 , ω g ) T
Sc坐标系中,珩磨轮在接触点处的线速度
v → c c = ω → c c × r → 1 c = ( - y 1 c ω c , x 1 c ω c , 0 ) T
在Sg坐标系中,被加工齿轮在接触点的线速度
v → g g = ω → g g × r → 2 g = ( - y 2 g ω g , x 2 g ω g , 0 ) T
将珩磨轮在接触点的速度由Sc坐标系转化到Sg坐标系中,
v → c g = M g c v → c c = ( y 1 c ω c , - x 1 c ω c c o s λ , x 1 c ω c s i n λ )
那么接触点处的相对速度
v → c g g = v → c g - v → g g = y 1 c ω c + y 2 g ω g - x 1 c ω c c o s λ - x 2 g ω g x 1 c ω c sin λ
在Sg坐标系中,被加工齿面的法向量为n2g;在齿面上一点处的切向量有无数条,利用齿面方程对齿面参数θ2求导,得到一个切矢量都垂直的方向向量
ρ → l 2 g = n → 2 g × ρ → θ 2 g
也是齿面的一条切向量,张成一个接触点处的切平面;相互垂直,并且都是单位向量,在接触点处构成了一个新的直角坐标系St
在坐标系St中,由于沿着齿面法向的速度并不起到切削的作用,因此只研究速度在切平面内的投影;被加工齿轮的速度和珩磨轮的速度方向上的投影分别为
v → c ρ θ 2 = v → c g · ρ → θ 2 g v → c ρ l = v → c c g · ρ → l 2 g
v → g ρ θ 2 = v → g g · ρ → θ 2 g v → g ρ l = v → c g g · ρ → l 2 g
因此,珩磨轮的速度vc在切平面上的投影 的夹角为 ψ c g = arccos | v → c ρ l | | v → c t | ; 被加工齿轮的速度vg在切平面上的投影 | v → g t | = ( v → g ρ θ 2 ) 2 + ( v → g ρ l ) 2 , 的夹角为 ψ g g = arccos | v → g ρ l | | v → g t | ;
(2)计算ap
两个鼓形圆柱相交部分,除去已经被磨削的部分,即为切屑近似模型,事实上这个切屑模型为相交部分的四分之一;
将近似模型沿轴向分割成无限小切片,那么每一层切片为一个外圆磨模型;每个小外圆磨模型中,磨削深度最大的位置是在过两个轴线的平面上;不同切削宽度位置的磨削深度分布在过两个轴线的平面与切屑近似模型相交的截面上,某磨削宽度位置bi处的磨削深度为api
bi 2=AcBc(2rcΓ1-AcBc),bi为线段CcCg到线段OgΓ1OcΓ1的距离,因为线段AcBc远小于直径,这个公式简化为
b i 2 ≈ 2 r cΓ 1 A c B c
A c B c ≈ b i 2 / 2 r cΓ 1
同理线段AgBg
A g B g ≈ b i 2 / 2 r gΓ 1
综上所述,磨削厚度api
a pi = a p - A c B c - A g B g = a p - b i 2 ( 1 / 2 r cΓ 1 - 1 / 2 r gΓ 1 )
上述公式中,rcΓ1和rgΓ1是未知的,ap为最大磨削厚度;
根据S1节中分析可知,rcΓ1和rgΓ1是分别为珩磨轮和被加工齿轮的齿面在接触点处于相对速度垂直的方向上的曲率半径;
(3)计算de
在磨削力模型当中,与砂轮和工件几何尺寸相关的参数为当量直径de,无论是哪种类型的磨削过程,都有一个对应的当量直径;根据外圆磨当量直径公式通过公式看出与当量直径相关的砂轮和工件的几何尺寸为各自的半径,ds和dw分别对应着S1中建立的近似模型中接触点处的曲率半径r和r;曲率半径为曲率的倒数,因此需要求得曲率半径r和r分别对应的法曲率k和k
当近似模型中,当量直径de
d e = 2 r c Γ r g Γ r c Γ + r g Γ
实际接触时,接触区域很小,在接触区域内所有位置的当量直径变化很小,因此把接触点的当量直径近似看做接触区域的当量直径;
S3.计算磨削力
磨齿过程被类比成S1中提出的磨削近似模型,在近似模型当中,每一个垂直于轴的截面都看做外圆磨,那么根据G.Wener的模型有
F′n=K′n[c1]γ[vgt/vct]2ε-1[api]ε[de]1-ε
在同一接触点,切屑磨削宽度上不同位置的vgt、vct和de是相同的,沿着磨削宽度不同位置的磨削深度是变化的,那么接触位置法向磨削力的模型表示为,
F n = ∫ o b a F n ′ db i
将单位宽度上磨削力模型带入上述公式,并进行积分,并得到
F n = K n [ c 1 ] γ [ v g t / v c t ] 2 ϵ - 1 [ d e ] 1 - ϵ { [ a p ] ϵ + 1 - [ a p - k p b a 2 ] ϵ + 1 }
其中kp=1/2rcΓ1-1/2rgΓ1
2.根据权利要求1所述的蜗杆状砂轮磨齿的磨削力模型,其特征在于:所述的磨削力模型适用于蜗杆砂轮磨齿和间齿珩齿加工方法。
CN201510927729.8A 2015-12-14 2015-12-14 蜗杆状砂轮磨齿的磨削力模型 Active CN105574249B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510927729.8A CN105574249B (zh) 2015-12-14 2015-12-14 蜗杆状砂轮磨齿的磨削力模型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510927729.8A CN105574249B (zh) 2015-12-14 2015-12-14 蜗杆状砂轮磨齿的磨削力模型

Publications (2)

Publication Number Publication Date
CN105574249A true CN105574249A (zh) 2016-05-11
CN105574249B CN105574249B (zh) 2018-08-28

Family

ID=55884379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510927729.8A Active CN105574249B (zh) 2015-12-14 2015-12-14 蜗杆状砂轮磨齿的磨削力模型

Country Status (1)

Country Link
CN (1) CN105574249B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383942A (zh) * 2016-09-12 2017-02-08 上海汽车变速器有限公司 蜗杆砂轮磨削斜齿轮的自然扭曲微观修形优化方法
CN106735612A (zh) * 2016-12-14 2017-05-31 北京工业大学 一种改善珩齿加工的方法
CN106971078A (zh) * 2017-04-11 2017-07-21 重庆大学 螺杆磨削过程中考虑运动参数的磨削力精确预测方法
CN107045556A (zh) * 2016-10-14 2017-08-15 台州市德嘉机电科技有限公司 螺旋曲面磨削过程中的磨削力预测方法
CN107133417A (zh) * 2017-05-25 2017-09-05 扬州大学 圆弧圆柱蜗杆的齿形检测方法
CN107153736A (zh) * 2017-05-11 2017-09-12 东北大学 一种修正的考虑鼓向修形的齿轮副啮合特性分析方法
CN108897955A (zh) * 2018-07-01 2018-11-27 北京工业大学 一种考虑螺纹升角和磨粒重叠系数的高速内螺纹磨削力计算方法
CN109858168A (zh) * 2019-02-14 2019-06-07 湘潭大学 一种基于离散元仿真的砂轮表面磨粒排布优化计算方法
CN114918494A (zh) * 2022-05-10 2022-08-19 重庆大学 齿轮展成磨削瞬时等效模型及其设计方法
CN115081143A (zh) * 2022-06-30 2022-09-20 合肥工业大学 一种蜗杆砂轮磨齿加工磨削力的仿真计算方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101106512B1 (ko) * 2007-08-17 2012-01-20 미츠비시 쥬고교 가부시키가이샤 배럴 웜형 공구의 드레싱 방법, 드레싱 장치 및 내접 기어 연삭기
CN103753395A (zh) * 2013-12-26 2014-04-30 西安交通大学 一种磨削力传感器及其测量方法
CN103995979A (zh) * 2014-06-05 2014-08-20 上海理工大学 非球面平行磨削加工的磨削力计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101106512B1 (ko) * 2007-08-17 2012-01-20 미츠비시 쥬고교 가부시키가이샤 배럴 웜형 공구의 드레싱 방법, 드레싱 장치 및 내접 기어 연삭기
CN103753395A (zh) * 2013-12-26 2014-04-30 西安交通大学 一种磨削力传感器及其测量方法
CN103995979A (zh) * 2014-06-05 2014-08-20 上海理工大学 非球面平行磨削加工的磨削力计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张祥雷 等: "基于多颗磨粒随机分布的虚拟砂轮建模及磨削力预测", 《航空学报》 *
谢桂芝 等: "工程陶瓷高速深磨磨削力模型的研究", 《机械工程学报》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383942B (zh) * 2016-09-12 2019-06-07 上海汽车变速器有限公司 蜗杆砂轮磨削斜齿轮的自然扭曲微观修形优化方法
CN106383942A (zh) * 2016-09-12 2017-02-08 上海汽车变速器有限公司 蜗杆砂轮磨削斜齿轮的自然扭曲微观修形优化方法
CN107045556B (zh) * 2016-10-14 2020-08-11 台州市德嘉机电科技有限公司 螺旋曲面磨削过程中的磨削力预测方法
CN107045556A (zh) * 2016-10-14 2017-08-15 台州市德嘉机电科技有限公司 螺旋曲面磨削过程中的磨削力预测方法
CN106735612A (zh) * 2016-12-14 2017-05-31 北京工业大学 一种改善珩齿加工的方法
CN106971078A (zh) * 2017-04-11 2017-07-21 重庆大学 螺杆磨削过程中考虑运动参数的磨削力精确预测方法
CN106971078B (zh) * 2017-04-11 2020-01-14 重庆大学 螺杆磨削过程中考虑运动参数的磨削力精确预测方法
CN107153736A (zh) * 2017-05-11 2017-09-12 东北大学 一种修正的考虑鼓向修形的齿轮副啮合特性分析方法
CN107153736B (zh) * 2017-05-11 2019-07-19 东北大学 一种修正的考虑鼓向修形的齿轮副啮合特性分析方法
CN107133417A (zh) * 2017-05-25 2017-09-05 扬州大学 圆弧圆柱蜗杆的齿形检测方法
CN107133417B (zh) * 2017-05-25 2020-04-10 扬州大学 圆弧圆柱蜗杆的齿形检测方法
CN108897955B (zh) * 2018-07-01 2022-09-09 北京工业大学 一种考虑螺纹升角和磨粒重叠系数的高速内螺纹磨削力计算方法
CN108897955A (zh) * 2018-07-01 2018-11-27 北京工业大学 一种考虑螺纹升角和磨粒重叠系数的高速内螺纹磨削力计算方法
CN109858168A (zh) * 2019-02-14 2019-06-07 湘潭大学 一种基于离散元仿真的砂轮表面磨粒排布优化计算方法
CN109858168B (zh) * 2019-02-14 2023-02-03 湘潭大学 一种基于离散元仿真的砂轮表面磨粒排布优化计算方法
CN114918494A (zh) * 2022-05-10 2022-08-19 重庆大学 齿轮展成磨削瞬时等效模型及其设计方法
CN114918494B (zh) * 2022-05-10 2023-12-26 重庆大学 齿轮展成磨削瞬时等效模型的设计方法
CN115081143A (zh) * 2022-06-30 2022-09-20 合肥工业大学 一种蜗杆砂轮磨齿加工磨削力的仿真计算方法及系统
CN115081143B (zh) * 2022-06-30 2024-02-13 合肥工业大学 一种蜗杆砂轮磨齿加工磨削力的仿真计算方法及系统

Also Published As

Publication number Publication date
CN105574249B (zh) 2018-08-28

Similar Documents

Publication Publication Date Title
CN105574249A (zh) 蜗杆状砂轮磨齿的磨削力模型
CN103592891B (zh) 运动学约束的复杂曲面五轴数控加工刀矢光顺方法
CN102004818B (zh) 自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法
CN103777568B (zh) 一种基于刃磨过程的整体式立铣刀容屑槽建模方法
CN103646141B (zh) 一种平底螺旋立铣刀正交车铣轴类零件的切削力建模方法
CN104635619A (zh) 基于刀具矢量插补的五轴数控加工方法
CN104759702B (zh) 圆柱齿轮的拓扑修形方法
Han et al. Cold orbital forging of gear rack
CN105817717A (zh) 用于通过对角线展成法对工件的齿轮制造机械加工的方法和装置
CN104331604B (zh) 基于元胞自动机的瓷砖抛光能耗建模的方法
Wang et al. Modeling of surface topography based on cutting vibration in ball-end milling of thin-walled parts
Kozlov et al. Modeling a cylindrical surface machined by a non-circular face tool
Lin et al. Additive manufacturing with secondary processing of curve-face gears
KR20150092319A (ko) 가변 리드 스크류를 구비한 로터 가공 방법
Olejnik et al. Calculation methodology for geometrical characteristics of the forming tool for rib cold rolling
Radzevich et al. Application priority mathematical model of operating parameters in advanced manufacturing technology
CN105728858A (zh) 一种大中型圆柱齿轮垂直交叉主轴式滚磨光整加工方法
Guo et al. Force prediction model for five-axis flat end milling of free-form surface based on analytical CWE
Han et al. Investigation on formation mechanism of surface texture and modeling of surface roughness with internal gear power honing
Pham et al. A practical approach for simulation and manufacturing of a ball-end mill using a 5-axis CNC grinding machine
Jiang et al. Prediction and simulation of cutting force in hypoid gear machining using forming method
Rababah A practical and optimal approach to CNC programming for five-axis grinding of the end-mill flutes
Ivanov et al. Recovery tread wheel pairs of machining
Yuan et al. The surface roughness modeling on turn-milling process and analysis of influencing factors
CN103065053B (zh) 热轧型钢飞剪传动曲轴旋转半径确定方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant