奥利康准双曲面齿轮的理论齿面推导及仿真

20 机械传动 2009年

文章编号:1004-2539(2009)02-0020-02

奥利康准双曲面齿轮的理论齿面推导及仿真

(河南科技大学机电工程学院, 河南洛阳 471003)

聂少武 邓效忠 李天兴

摘要 根据奥利康制准双曲面齿轮的加工方法和切齿加工原理,分析了刀具、摇台及工件的相对位置和相对运动关系,建立了切齿加工坐标系。运用空间啮合原理,推导出理论齿面方程。将理论齿面进行数据离散后建立方程组,数值分析求解得到齿面离散点坐标。利用三维软件进行了理论齿面的三维仿真,从而为齿面接触分析(TCA)、齿面误差测量以及有限元分析等方面的研究提供了理论基础。

关键词 奥利康制准双曲面齿轮 加工坐标系 理论齿面 三维仿真

0 引言

曲齿锥齿轮由于传动平稳、承载能力高、重合度大等优点[1],现在已成为汽车的关键传动部件。目前,曲齿锥齿轮主要分为 格里森 制和 奥利康 制,简称 格 制和 奥 制。 格 制锥齿轮在国内应用比较普遍,研究得也相当充分,而 奥 制锥齿轮在国内应用极少,对其研究还处于起步阶段,相关文献资料也较少。国内现在引进的轿车后桥齿轮有一些采用的是 奥 制,很有必要对 奥 制锥齿轮进行理论研究。

基于此,本文根据 奥 制准双曲面齿轮的加工原理和加工方法,建立了切齿加工数学模型,推导了理论齿面方程,并进行了图形仿真,这将为齿面接触分析(TCA)、齿面误差测量以及有限元分析等方面的研究提供理论基础。

奥利康铣齿刀盘是整体刀盘,刀盘上装z0组刀齿,每组刀齿有3个刀片,分为粗刀、内切刀、外切刀。图1为刀具坐标系。图中, 为刀倾角,P0为刀

齿节点,T为产形轮分度平面,T0为刀齿节平面,T1为过刀齿节点与刀盘轴线垂直的平面。G为T1面与刀刃的交点,P为刀刃上任一点,u为刀具参数, 01为T1面内刀齿方向角,hx1为由变位引起的偏移距,Se0(Xe0,Ye0,Ze0)为T1面内与刀盘固连的坐标系,Se1(Xe1,Ye1,Ze1)为T1面内辅助坐标系,初始位置,Se0与Se1重合。 i(i=1,2)为刀盘绕自身轴线当前转角。Se2(Xe2,Ye2,Ze2)为T0面辅助坐标系,Se(Xe,Ye,Ze)为分度平面内与刀具固连的坐标系,!k为刀具齿形角(外刀为负,内刀为正),P0G为刀具半径修正量。

1 加工原理和加工方法

奥 制准双曲面齿轮为等高齿摆线准双曲面齿轮,即在齿长方向齿线为延伸外摆线,沿分锥母线齿高

相等。其加工原理采用假想平面齿轮原理,连续分度双面法加工,要求刀盘转过一组齿,产形轮也转过一个齿。刀刃绕刀盘轴线自转,又相对于产形轮轴线公转,在平面产形轮上形成延伸外摆线产形面。被加工齿轮的齿面由产形轮的产形面展成而成。

奥 制准双曲面齿轮的加工方法主要有刀倾全展

成法和刀倾半展成法。刀倾全展成法通过刀轴倾斜来修正齿面接触区,可以同时对大小轮齿面修正,也可以只修正小轮齿面。刀倾半展成法加工时大轮无展成运动,小轮采用与大轮近似的圆锥产形轮对偶加工[2]28-29。本文研究的内容均采用刀倾全展成法。

2.2 产形面坐标系

加工奥利康制准双曲面齿轮时,刀盘一方面绕自身轴线自转,一方面绕产形轮轴线公转,刀盘相对产形轮的运动可视为产形轮不动,刀盘上一动圆在与产形轮固连的一定圆上做纯滚动,刀刃上一点在产形轮上形成的轨迹即为延伸外摆线,刀刃上所有点的轨迹面形成产形轮的产形面。图2为延伸外摆线形成过程。

图3为加工右旋齿轮时产形面的形成坐标系,产形轮左旋,刀盘中心位于第一象限,假定产形轮不动,刀盘绕刀轴逆时针转动,同时还绕产形轮轴线逆时针T

[3]

图1 刀具坐标系

2 坐标系的建立

2.

第33卷 第2期 奥利康准双曲面齿轮的理论齿面推导及仿真 21

中。Op0(Op)为摇台轴线与分度平面T的交点,M点为切齿参考点,坐标系Sp0(Xp0,Yp0,Zp0)为产形轮辅助坐标系,Sp(Xp,Yp,Zp)为与产形轮固连的坐标系,Oe为刀盘轴心在分度平面T上的投影,Se(Xe,Ye,Ze)为与刀盘固连的坐标系。初始位置,坐标系Sp0与Sp重合。Se与Sp0固连,转动过程中相对位置保持不变。 p2为刀盘绕产形轮轴线转动的角度。S2为径向刀位, q2为刀位极角,∀2为一已知角,其计算可参见文献[2]222-223

轮坯逆时针转动(从小端看)。其余与右旋齿轮情况类似。

。

图5 右旋齿轮加工标系

图2 延伸外摆线的形成

图4为加工左旋齿轮时产形面的形成坐标系,加工 奥 制左旋齿轮时,产形轮右旋,刀盘中心位于第四象限,采用左旋刀盘,刀盘绕刀轴轴顺时针转动,同时还绕产形轮轴线顺时针转动。各个坐标系之间的关系

与加工右旋齿轮时的情况类似。

图6 左旋齿轮加工坐标系

3 齿面理论数据的产生

3.1 齿面方程的推导

将刀具实体矢量从刀具坐标系经过一系列坐标变换转化到与被加工齿轮固连的坐标系中即可得到理论

[4-齿面方程。

6]

图3 左旋产形轮形成坐标系图4 右旋产形轮形成坐标系

2.3 切齿加工坐标系

图5为右旋齿轮加工坐标系,摇台逆时针转动,轮坯顺时针转动(从小端看)。Sm(Xm,Ym,Zm)是机床坐标系,与机床固连,XmOYm位于产形轮的分度平面内,Sp(Xp,Yp,Zp)为产形轮坐标系,与摇台固连,产成过程中,随摇台一起转动,角#g2为Sp相对Sm的当前转角,过轮坯轴线与分度平面垂直的平面为机床水平面,Sn(Xn,Yn,Zn)、坐标系为辅助坐标系,Sn与Sm平行,Xa与轮坯轴线重合,S2(X2,Y2,Z2)坐标系与轮坯固连,展成过程中,绕齿轮轴线转动,角#2为轮坯的当前转角。初始位置,S2与Sa坐标系重合。O2为设计交叉点,O2!为节锥顶点,A点为切齿加工公垂线在轮坯轴线上的垂足, 2为安装角。

,在矩阵变换中,有以下关系式成立

∃ z0pipi

ip0===z

0iip∃i#izp

mpi===

pigizi

#gi= gi∀%i

式中,i=1.2,分别表示小轮和大轮;ip0为产形轮与刀盘速比;∃pi为产形轮角速度;∃0i为刀盘角速度;∃i为轮坯角速度;zp为产形轮齿数;zi为轮坯齿数;z0为刀组数; gi为摇台转角;%i为产形轮参考点偏离角,有正负之分,符号根据齿轮旋向和偏置情况确定。

(1)由图1知,刀刃上任意一点P在Se0坐标系中矢量方程为

r0+dr0cos 01+usin!kcos 01

re0=

-dr0sin 01-usin!ksin 01ucos!k

1

式中,dr0=P0G,r0=Oe1P0,r0为刀盘名义半径。

#

产形面齿面方程

rpi(ui, i)=M 1re0

其中,M 1为Se0 Sp的转换矩阵。

M 1=Mpp0Mp0eMee2Me2e1Me1e0产形面单位法矢为

Ni

npi=

|Ni|

rpi rpi

其中 Ni=

ii

(2)求得产形面方程后,采用求共轭齿面的方法求解齿面方程和单位法矢。

ri(ui, i, gi)=M 1rpini=L inpi

其中,M 1为Sp Si的转换矩阵,L i由n i变换而得。

M 1=MiaManMnmMmp

由于矩阵繁多、结果复杂,限于篇幅有限,这里就不给出具体结果。

(3)在产形轮与轮坯的啮合过程中,在Sm坐标系下有啮合方程成立

F(ui, i, gi)=0由啮合方程可求得 gi= gi(ui, i),代入齿面方程可得ri(ui, i)。

3.2 划分齿面网格

过齿轮轴截面做旋转投影面,如图7所示,对齿面进行网格规划。齿面上与M!对应的点满足以下关系式

xi(ui, i, gi)=XL

22

y2i(ui, i, gi)+zi(ui, i, gi)=RL

[7]

pi

nmivm=

数,表2为刀具相关参数(见文献[2]337)。

表1 奥 制准双曲面齿轮几何参数和加工参数名 称轴夹角/(∃)偏置距/mm参考点法向模数/mm

全齿高/mm齿顶高/mm参考点分度圆半径/mm参考点螺旋角/(∃)

节锥角/(∃)齿数旋向齿宽/mm刀倾角/(∃)刀位/mm高变位移距/mm产形轮齿数产形轮参考点螺旋角/(∃)

分锥顶距/mm轮坯安装角/(∃)垂直轮位/mm轮位修正量/mm床位修正量/mm机床型号

1.94597516.42713617.7481.589333-1.290249

S174.6725.9407250.0105517.5888410左旋34296.17177-1.333942

41.82131938.523873

-1.99984971.30395-2.162-3.29724-4.197658

S17

小 轮

90203.3382867.50

285.73522137.1344471.9867941右旋30296.171351.333752大 轮

图7 齿面网格规划

将此关系式与啮合方程联立组成非线性方程组,利用牛顿#拉斐森算法,通过VC++和MATALB混合编程求解可得到齿面上规划点的三维坐标及单位法矢。

图8 齿面仿真表2 刀具参数

名称刀盘型号T0面刀齿方向角/(∃)T1面刀齿方向角/(∃)内刀半径修正量/mm

小轮1FS13-88L114.24814514.239818-0.737419

大轮2FS13-88R214.24814514.2398180.665472

4 图形仿真

以一对奥利康准双面齿轮为例,根据上述推导出的理论齿面及计算结果对齿轮进行了三维图形仿真,1页)

序号5

A0x

A0y

A1x

A1y

是否满足要求不合要求(存在分支)不合要求(无曲柄)是否满足要求满足要求满足要求满足要求满足要求不合要求(存在分支)不合要求(无曲柄)

[2] 罗佑新.平面四杆机构近似导引综合的优化#混沌方法[J].机械

设计,2003,20(11):34-36.

[3] 张龙庭,董湘怀.非线性方程组的混沌求解方法及在机构学中的

应用[J].机床与液压,2003(1):165-169.

[4] 罗佑新,黎大志.求一般6-SPS并联机器人机构的全部位置正解

的混沌迭代方法[J].工程设计学报,2003(2):70-74.

[5] 谢进,陈永.混沌理论和平面连杆机构综合[J].机械科学与技术,

2000,19(4):524-526.

[6] 谢进,陈永.计算运动学中牛顿迭代法的混沌现象[J].机械传动,

2000,24(1):4-6.

[7] 谢进,陈永.基于混沌的刚体导引问题Burmester点的求解方法

[J].中国机械工程,2002,13(7):608-710.

[8] LiuWB,ChenGR.Dynamicalanalysisofachaoticsystemwithtwodou

ble-scollchaoticattractors[J].J.Bifurc.&Chaos,2004,14(3):971-998.

[9] 陈志坚,孙克辉,张泰山.Liu混沌系统的非线性反馈同步控制

[J].物理学报,2005,54(6):2580-2583.

[10] WolfA,SwiftJB,SwinneyHLetal.DeterminingLyapunovexponents

fromatimeseries[J].PhysicaD:NonlinearPhenomena,1985,16(3):285-317.

[11] HendrikRichter.ThegeneralizedHenonmaps:examplesforhigher-di

mensionalchaos[J].InternationalJournalofBifurcationandChaos,2002,12:1371-1381.

[12] 刘安心,杨廷力.机械系统运动学设计[M].北京:中国石化出版

社,1999:60-67.

收稿日期:20081113 收修改稿日期:20081208

项目资助:湖南省 十一五 重点建设学科(机械设计及理论)(湘教通

2006180)、国家自然科学基金资助(No:50845038)、湖南省自然科学基金(07JJ3093)、湖南科技厅计划项目(2007FJ3030)。

作者简介:车晓毅(1966-),女,湖南常德人,副教授

-1.08806-1.36136-1.15230-0.28868

6序号123456

3.41455B0x3.414553.845715.127383.845715.127383.84571

-2.00586

B0y-2.00586-2.860821.86785-2.860821.86785-2.86082

1.89649B1x1.896491.663895.779981.663895.779981.66389

3.82482B1y3.824825.039981.125245.039981.125245.03998

5 结论

本文首次提出了基于新的Liu混沌系统产生混沌序列求解机构综合问题的方法,对平面铰链四杆机构刚体导引综合进行了研究。研究结果表明,该方法既不要构造初始机构,同时又可以得到多组机构方案,为机构综合的多方案优选打下了基础,对平面及空间机构的综合研究具有应用价值。本文提出的方法同样适用机构近似综合。

参

考

文

献

[1] 罗佑新,郭惠昕,肖伟跃.机构综合的牛顿优化迭代混沌方法研究

[J].机械传动,2005,29(6):11-13.

(上接第22页)

名称

外刀半径修正量/mm内刀齿形角/(∃)外刀齿形角/(∃)

小轮10.75582820.90527-19.08958

大轮2-0.6821116.302609-23.697467

参考文献

[1] 曾韬.螺旋锥齿轮设计与加工[M].哈尔滨:哈尔滨工业大学出

版社,1989.

[2] 董学朱.摆线齿锥齿轮及准双曲面齿轮设计和制造[M].北京:机

械工业出版社,2002.

[3] 北京齿轮厂.螺旋锥齿轮[M].北京:科学出版社,1974:12-15.[4] 吴序堂.齿轮啮合原理[M].北京:机械工业出版社,1982.[5] 李巍,李剑锋,王青云.克林贝格摆线齿锥齿轮齿面方程及图形仿

真[J].机械科学与技术,2006,25(12):1462-1466.

[6] LitvinFL,FuentesA.GearGeometryandAppliedTheory[M].2nded.

UnitedKingdom:CambridgeUniversityPress,2004.

[7] 李天兴,邓效忠,魏冰阳.基于一维测头准双曲面齿轮齿面偏差的

测量[J].中国机械工程,2007,18(8):958-961.收稿日期:20080602

作者简介:聂少武(1983-),男,河南洛阳人,硕士研究生

在图8中,图8a和图8b分别为小轮和大轮三维模型,图8c和图8d分别为大轮和小轮啮合模型。

5 结论

依据奥利康制准双曲面齿轮的加工方法和切齿原理,建立切齿加工坐标系,能够推导出理论齿面并绘出了三维模型。对所求得的齿面进行仿真,验证了数学模型和理论推导的正确性。在此研究的基础上可以进行齿面接触分析、齿面误差测量和有限元分析等方面的研究。

英文摘要 2009年

ABSTRACTS&KEYWORDS

JOURNALOFMECHANICALTRANSMISSIONVol.33.No.2,2009

ResearchonaNovelTransmissionMechanismbasedonthePrinci pleofSpaceCurveMesh%%%%%%%%%%%%%%%%

%%%%%%%%%ChenYangzhi,HuQiang,LuoLiang(1)Abstract Amicro-elastic-mesh-wheel-transmission(MEMWT)mechanismbasedontheprincipleofspacecurvemeshisintroduced.Firstofall,thebasicprincipleandcharacteristicsofthenoveltransmis sionmechanismareelaborated.Then,theresearchprogressonthisnovelmechanismissummarized,includingthetheoryofspacecurvemesh,thedesignandcalculatingsystemofspacecurvemeshtransmis siononthebasisofMATLAB,kinematicsperformancetest,transmis sionforceanalysisanddynamicsperformancetest.Finally,problems,whichshouldbesolved,aresimplyintroduced.

Keywords: Meshtransmission Microtransmissionmechanism Micro-elastic-mesh-wheel Spacecurvemesh

ResearchofSimultaneouscontrolModelforHobblingThreeOrderEllipseHelicalGear%%%%%%%%%%%%%%%%%%

%%%%%%%HuChibing,TianFangyong,CaoLiangwei(5)Abstract Onthebasisofthedevelopedsystem CNCgear-hobblingsystem thatcoverswithCNCmachiningofnon-circularverticalgear,themeshingtheoryofhelicalgearsandtoolrack-skewmethods,thesimultaneouscontrolstrategyofhobblingthreeorderellipsehelicalgearisproposed,andthegraphofCNCsimultaneouscontrolmodelisde signed.Hobblingexperimentthatthree-orderellipsehelicalgearcanbegainedprovesit.Whentwogearsengaged,theycanturnflexiblyandthebacklashiswelldistributed.

Keywords: CNCThreeorderelliptic Helicalgear Hobbling Si multaneouscontrol

ResearchonEnhancingtheContactFatigueStrengthofCastIronMachineElementswithAlloyCoatingandAlloying%%%%%

%%%%%%%%%%%%%%%%%%%%HuangJi(9)Abstract Theworkingperformanceofalloycoatingandafteritsalloy ingofHT250graycastironmachineelementssurfacewasstudied.Uti lizingtheelectronicmicroscope(EM)andOpticalmicroscope(OM)ananalysisuponsurfaceappearance,structureofsub-surfaceandlongi tudinalsection,surfaceroughness,compositionofalloycoatingandmi crohardness,bindingstrengthbetweenalloycoatingandmatrixmetal.Etc,wascarriedout.Testingresultofcontactingfatigueshowsthatthecoatingcouldenhancenotablythecontactingfatiguelife.

Keywords: Alloycoating Alloying Strengthofcontactingfatiguestrength Greycastiron Machineelements

AnalyticalMethodforComputingandAdjustingContactRatioofNon-circularGearPairs%%%%%%%%%%%%%%%%%%%%%%%%%WangLei,ZhangRui,SuZhijian(14)theofis

scribed.Firstly,theanalyticalmodelsoftoothprofilesandaddendum

curvearedevelopedgeometrically.Thentherelationshipbetweencon tactratioandpolarangleisestablished.Furthermorethefactorseffect ingoncontactratioareanalyzed.Withtheconcernoftheundercuttingandthepointingoftheteeth,theconceptofboundarycontactratioisintroduced.Finallyasystematicmethodtocontrolcontactratioaccu ratelyisproposed.

Keywords: Non-circulargear Contactratio Addendumcurve Working-pressure-anglefunction

StudyontheControlSystemofGearTransmissionbasedonVirtu alPrototyping%%%%%%%%%%%%%%%%%%%%

%%%%FangZifan,ShuGang,HeKongde,ZhaoXinze(17)Abstract AmixedtypeofgeartransmissionwasdynamicsmodeledbyPro/EandADAMS.Basedonthedynamicsmodel,theoutputdis placementofthePIDtwo-loopcontrolsystemforthegeartransmissionwasmodeledbythecontroltoolkitofADAMS/View.Thevirtualproto typeofthecontrolsystemforgeartransmissionisassembled.Thesimu lationresultsofseveralclassicalinputsignalsindicatedthatthevirtualprototypeofthecontrolsystemforgeartransmissionisvalidity,andthecontrolsystemcanmeettheengineeringcontrolaim.

Keywords: Geartransmission Controlsystem Two-loopPID Virtualprototyping Simulation

TheoreticalToothSurfaceDeductionandSimulationofOrionHy poidGears%%%NieShaowu,DengXiaozhong,LiTianxing(20)Abstract AccordingtothecuttingmethodandcuttingprincipleofOri onhypoidgears,therelativepositionsandmovementsamongcutterhead,cradleandworkpiecewereanalysis,andthecuttingcoordinatesystemswereestablished.Thetheoreticaltoothsurfaceequationswerederivedonthebasisofspaceengagementtheory.Theequationswerees tablishedafterdissertingthetheoreticaltoothsurface,andthecoordi natesofdiscretepointsweresolvedbythenumericanalysis.Onthemeantimethesimulationoftheoreticaltoothsurfacewascompleted,whichwillprovidethetheoreticalbasisforTCA,themeasurementoftoothsurfacedeviationandthefiniteelementanalysis.

Keywords: Orionhypoidgears Cuttingcoordinatesystem Theo reticaltoothsurface 3Dsimulation

PIDControlModelfortheContinuousConveyorofDisc-tubeAs semblybasedontheBPNeuralNetwork%%%%%%%%%

%%%%%%%%%ShiShulin,LuanLijun,LiXiaozhi(23)Abstract UsingBPneuralnetwork,Conveyorsystemnetworkmodelforthecontinuousconveyorofdisc-tubeassemblyandaspiralconvey orisestablished.ByoptimizationtrainingwithPID,andadjustingtheoutputparametersKP,KI,KD.Atanytime,theoptimumvalueofoutputisTheresultofis