JournalofMechanicalStrength
2022ꎬ44(1):162 ̄168
DOI:10 16579/j.issn.1001 9669 2022 01 022
∗20200916收到初稿ꎬ20201011收到修改稿ꎮ国家自然科学基金项目(51975379)ꎬ上海市自然科学基金项目(19ZR1435500)资助ꎮ∗∗宋有硕ꎬ男ꎬ1987年生ꎬ江苏南通人ꎬ上海理工大学机械工程学院讲师ꎬ博士ꎬ研究方向为机械结构的设计与分析∗∗∗马㊀跃ꎬ男ꎬ1996年生ꎬ江苏溧阳人ꎬ上海理工大学机械工程学院硕士研究生ꎬ研究方向为机械结构的设计与分析ꎮ
恶劣工况下越野车用等速万向节接触应力和保持架强度分析∗
CONTACTSTRESSANDCAGESTRENGTHANALYSISOFCONSTANTVELOCITYJOINTFOROFF ̄ROADVEHICLEUNDERBADWORKING
CONDITIONS
宋有硕∗∗㊀马㊀跃∗
∗∗
(上海理工大学机械工程学院ꎬ上海200093)
SONGYouShuo㊀
MAYue
(CollegeofMechanicalEngineeringꎬUniversityofShanghaiforScienceandTechnologyꎬ
Shanghai200093ꎬChina)
摘要㊀军用越野车球笼式等速万向节受到的载荷比传统汽车大ꎬ且工况恶劣ꎬ有可能发生接触磨损ꎬ以及大轴间角度下保持架的断裂ꎬ因此对恶劣工况下接触应力和保持架强度的研究非常重要ꎮ以某典型的军用越野车球笼等速万向节为对象ꎬ首先对钢球与内外滚道之间的接触应力进行理论解析和仿真ꎬ将仿真和解
析结果进行对比ꎬ验证仿真模型的正确性ꎻ其次ꎬ根据万向节关键部位之间的载荷传递关系和物理实验ꎬ分析保持架断裂的原因和过程ꎻ最后ꎬ进行极限轴间角下静扭实验ꎬ并对保持架断面进行硬度测试ꎬ通过硬度 ̄强度转换关系得到沿厚度的强度分布ꎬ结果表明保持架在正常工况下强度满足要求ꎬ但是恶劣工况下有可能会破碎ꎮ可以为越野车用等速万向节设计及制造提供理论参考ꎮ
关键词㊀万向节㊀保持架㊀接触分析㊀应力分析㊀强度中图分类号㊀U463㊀216+.3
Abstract㊀Theloadofballcagetypeconstantvelocityjointofmilitaryoff ̄roadvehicleislargerthanthatoftraditional
vehicleꎬandtheworkingconditionisbadꎬwhichmaycausecontactwearꎬaswellasthefractureofthecageundertheanglebetweenlargeaxes.Thereforeꎬitisveryimportanttostudythecontactstressandthestrengthofthecageunderthebadworkingcondition.Inthispaperꎬatypicalmilitaryoff ̄roadvehicleballcageconstantvelocityjointistakenastheobject.Firstlyꎬ
thecontactstressbetweenthesteelballandtheinsideandoutsideracewayistheoreticallyanalyzedandsimulated.Thesimulationandanalyticalresultsarecomparedtoverifythecorrectnessofthesimulationmodel.Secondlyꎬaccordingtotheloadtransferrelationshipbetweenthekeypartsoftheuniversaljointandthephysicalexperimentꎬthereasonandprocessofcagefractureareanalyzed.Finallyꎬthestatictorsiontestundertheultimateaxialangleiscarriedoutꎬandthehardnesstestiscarriedoutonthesectionofthecage.Thestrengthdistributionalongthethicknessisobtainedthroughtherelationshipbetweenhardnessand
strengthconversion.Theresultsshowthatthestrengthofthecagemeetstherequirementsundernormalworkingconditionsꎬbutmaybebrokenunderbadworkingc
onditions.Thisresearchcanprovidetheoreticalreferenceforthedesignandmanufactureofconstantvelocityjointforoff ̄roadvehicle.
Keywords㊀UniversaljointꎻCageꎻContactanalysisꎻStressanalysisꎻStrengthCorrespondingauthor:SONGYouShuoꎬE ̄mail:wssys2006@163.comTheprojectsupportedbytheNationalNaturalScienceFoundationofChina(No.51975379)ꎬandtheNaturalScienceFoundationofShanghai(No.19ZR1435500).
Manuscriptreceived20200916ꎬinrevisedform20201011.
㊀㊀前言
等速万向节是车辆传动系统中的重要组成部分ꎮ现今对普通乘用车等速万向节的设计研究已经比较成
熟ꎬ但是对于特种越野车万向节的设计及制造仍然高度依赖国外技术ꎬ需要参考德国GKN和日本NTN
等国外企业提供的设计图纸和工艺规范ꎮ军用越野车球笼等速万向节受到的载荷比传统汽车大ꎬ工况恶劣ꎬ有
㊀第44卷第1期宋有硕等:恶劣工况下越野车用等速万向节接触应力和保持架强度分析163㊀㊀
可能发生接触磨损及大轴间角度下保持架的断裂ꎬ因此对恶劣工况下接触应力和保持架强度研究非常重要ꎮ
但是现今对特种车辆球笼万向节的研究仍缺乏系统性ꎬ大部分研究集中于普通乘用车用万向节ꎬ但是对于普通乘用车而言ꎬ行驶工况并不恶劣ꎬ不太可能发生万向节接触应力过大或保持架的破裂现象ꎮ现今对恶劣工况下军用特种车万向节接触应力和保持架强度的研究较少ꎬ对于万向节接触应力和保持架在大轴间角㊁大扭矩工况下的仿真模型有待进一步修正ꎮ
卢曦等对万向节的装配以及球面间的配合进行了分析[1]ꎮ侯文亮等对保持架进行了仿真研究ꎬ分析了保持架强度分布规律[2]ꎮ陈静等对万向节内部的接触应力进行了分析ꎬ表明内滚道比外滚道更容易受损[3 ̄4]ꎮ缪莹赟等对钢球传递运动中产生的摩擦力进行了研究ꎬ求解了各种复杂工况下的滚动轴承的力学特性[5]ꎮ李玉璇和王立新通过Hertz接触理论计算得出万向节沟道的接触应力[6]ꎮRyuI等研究了万向节的准静态断裂过程ꎬ并分析了各个部件之间的相互作用力[7]199 ̄205ꎮ石宝枢等研究了汽车球笼式等速万向节的结构主参数的合理范围[8]ꎮHanJN和HeinC等研究涉及了数值模拟
中的网格划分技术及塑性变形的应力求解[9 ̄10]ꎮ在接触应力方面ꎬ除了Hertz接触理论ꎬ很多学者还提出了非线性弹簧阻尼接触应力模型ꎬ论述了如何确定阻尼系数[11 ̄12]ꎮ
本文以某典型的军用越野车球笼等速万向节为对象ꎬ首先对钢球与内外滚道之间的接触应力进行理论解析和仿真ꎬ将仿真和理论解析结果进行对比ꎬ验证仿真模型的正确性ꎻ其次ꎬ根据万向节关键部位之间的载荷传递关系和物理实验ꎬ分析保持架断裂的原因和过程ꎻ最后ꎬ进行极限轴间角下静扭实验ꎬ并对保持架断面进行硬度测试ꎬ通过硬度 ̄强度转换关系得到沿厚度的强度分布ꎬ为越野车用球笼等速万向节的设计提供理论参考ꎮ
1㊀万向节结构
以某军用越野车用RF125型球笼等速万向节为对象ꎬ万向节的结构如图1所示ꎬ由钟形壳1㊁钢球2㊁保持架3和星形套4组成ꎮ钟形壳(又称外套)和星形套(又称内套)各有六个凹槽和六个球面ꎬ凹槽和球面把钢球固定住使得钟形壳和星形套同步旋转ꎮ保持架上的穿孔部分是保持架窗口(3a)ꎬ其固定钢球并将钢球限制在中心内ꎬ笼柱(3b)位于笼窗之间ꎬ它们的球面与钟形壳和星形套相互作用
图1㊀球笼式等速万向节结构示意图
Fig.1㊀Schematicdiagramofball ̄cagetypeCVjoint
万向节各个零件的材料及机械特性如表1所示ꎮ
表1㊀各零部件的材料及特性
Tab.1㊀Materialsandcharacteristicsofparts
零件
Component材料名称
Material
弹性模量
Young s
modulus/GPa
密度
Density/
(kg/m3)
泊松比
Poisson s
ratio
钟形壳
Outerrace55#钢55#steel20178500 3
星形套
Innerrace20CrMnTi20678500 29
钢球
BallGCr1520678500 29
保持架
Cage20CrMnTi20678500 292㊀接触应力解析模型
钢球与内㊁外套椭圆式沟道为点接触ꎬ传递扭矩时6个钢球同时承载ꎬ作用在钢球与沟道接触点的法向压力为
Q=M
ZRPsinβ(1)式中ꎬM为传递力矩ꎻRP为钢球分布圆半径ꎻβ为接触压力角ꎻZ为钢球数ꎮ
假设钢球与沟道接触表面处于弹性应力状态ꎬ而且接触区域的尺寸比沟道接触点曲率半径小的多ꎬ引用赫兹理论来计算沟道表面接触应力ꎮ钢球与内外滚道几何模型如图2所示ꎮ
接触点到内沟道沟底沿X方向的距离
h=2fsin2β
1+2fcosβ
2(2)㊀㊀内沟道 ̄钢球接触当量主曲率半径
Rx=d(R+h)
2(R+h)+dꎬRy=
2f-1d(3)式中ꎬf为滚道瞬时曲率系数ꎬ其值为0 525ꎮ
则钢球与内外套滚道的最大接触应力Pmax为
Pmax=32QEᶄ2
36πε2KðR2
æ
è
ç
ö
ø
÷
1/3
(4)式中
㊀164㊀机㊀㊀械㊀㊀强㊀㊀度2022年
图2㊀圆弧滚道几何模型Fig.2㊀Arcracewaygeometricmodel
Eᶄ=
21-ν2aEa
168汽车网+
1-ν2bEb
为当量弹性模量ꎻε=1 003+
0 5968
Ry/Rx
为第二类椭圆积分ꎻK=1 0339Ry/Rx()0 6360
接触椭圆参数ꎻðR=Rx+Ry为主曲率半径ꎻνa㊁νb分别为两接触体的泊松比ꎮ
此型号万向节的相关参数如下:钢球直径d=23 813mmꎬ接触角β=45ʎꎬ内外沟道接触点沟曲率
系数fi=fe=0 525ꎬ节圆半径RP=39 6mmꎬ内外沟道底径分别为Ri=27 644mmꎬRe=51 557mmꎬ最大传递扭矩为Mmax=5782Nmꎬ材料参数E=2 06ˑ1011N/
m2ꎬ泊松比μ=0 3ꎮ将上述参数带入上述理
论解析公式可得接触应力的解析解ꎬ求得零轴间角下钢球与万向节外沟道的接触应力为4465MPaꎬ与内沟道的接触应力为4386MPaꎮ
3㊀仿真模型
3 1㊀仿真模型的建立
图3所示为仿真模型的边界和载荷条件ꎬ在钟形壳上施加固定约束ꎮ将驱动扭矩施加在星形套花键上ꎬ同时约束花键除切向以外的转动ꎮ保持架和钢球不施加任何约束ꎬ由各个零件之间的接触关系自动达到平衡ꎮ
根据万向节各个部件之间的运动关系ꎬ选择接触方式为面面接触ꎮ取接触面为带摩擦接触ꎬ在二硫化钼润滑条件下设置摩擦因数为0
1ꎮ
图3㊀边界和载荷条件Fig.3㊀Boundaryandloadconditions
为了使接触区域网格足够细ꎬ计算结果更加准确ꎬ先运用整体网格划分ꎬ尺寸为3mmꎬ然后用局部网格细化的方法对接触体进行网格细化ꎬ网格划分采用六面体单元ꎬ并尽量使单元规则细化ꎮ根据赫兹接触应力ꎬ在接触点中心会产生接触椭圆ꎬ因此在内外套与钢球接触处切割出圆弧形状ꎬ设置其网格尺寸为0 1
mmꎮ由于整体的计算量会非常大ꎬ因此取六分之一模型研究ꎬ扭矩相应取六分之一ꎮ图4所示为沟道 ̄钢球接触处的网格细化模型
图4㊀沟道 ̄钢球接触部位网格细化
Fig.4㊀Themeshrefinementofthegroove ̄ballcontact
3 2㊀仿真模型结果与验证
加载2500Nm和5782Nm扭矩后钢球与钟形壳最大接触主应力的运算结果如图5所示ꎮ从图5中可以看出钟形壳接触带呈椭圆分布ꎬ随着扭矩的提高ꎬ接触部位的最大主应力也随着变大ꎬ最大值约为
4125MPaꎮ
图5㊀钢球与钟形壳最大接触主应力的剖面视图
Fig.5㊀Maximumcontactprincipalstressbetweenballandouterrace
当承受转矩为5782Nm时ꎬ零度轴间角下ꎬ钢球
㊀第44卷第1期宋有硕等:恶劣工况下越野车用等速万向节接触应力和保持架强度分析165㊀㊀
与星形套㊁钟形壳的整体接触应力如图6所示ꎮ万向节内钢球与星形套滚道最大接触应力约为4658MPaꎬ这一数值明显大于钟形壳滚道的应力ꎬ因此ꎬ在实际工作过程中ꎬ更容易磨损的是星形套内滚道ꎮ这与球笼式万向节在实际工作过程中ꎬ内滚道更容易产生坑洼㊁点蚀等现象一致
图6㊀零度轴间角下钢球与星形套㊁钟形壳的整体
接触应力
Fig.6㊀Theoverallcontactstressnebulogramofsteelballꎬstar
sleeveandbellshellunderzerointeraxialangle
零度轴间角下接触应力的仿真值和理论解析值对比结果如表2所示ꎮ
表2㊀接触应力理论值与仿真值对比
Tab.2㊀Comparisonoftheoreticalandsimulationvaluesof
contactstress
沟道
Channel理论值Theoreticalvalue/MPa仿真值Simulationval
ue/MPa相对误差Relativeerror/%
外沟道Outerchannel446541258内沟道Innerchannel
4386
4658
通过表2发现ꎬ内㊁外沟道上的最大接触应力相对误差分别为6%与8%ꎬ此相对误差产生的原因可能是理论解析假设条件理想化以及仿真方法本身的误差性ꎬ此误差大小在合理误差范围内ꎬ可以证明本文所建立的仿真模型的准确性和可靠性ꎮ
4㊀恶劣工况下保持架强度分析
4 1㊀大轴间角及扭矩下保持架应力分布
保持架在大轴间角时很容易损伤ꎬ为了分析保持架在大轴间角下的损伤方式和原因ꎬ利用上述建立的仿真模型ꎬ进行大轴间角及极限扭矩下保持架应力分布的求解(轴间角为40ʎꎬ扭矩Mmax=5782Nm)ꎬ求解结果如图7所示ꎮ图7显示了一个破损的保持架样
图7㊀保持架等效应力分布与试验损伤样品
Fig.7㊀Equivalentstressdistributionofretainerandtestdamagedsamples
4 2㊀万向节关键部位之间的相互作用力
为了弄清保持架损伤的原因和过程ꎬ首先需要对万向节关键部位之间的相互作用力进行分析ꎬ进而求解出保持架的受力状态ꎮ钢球通过几个平行的路径传递扭矩ꎮ六个钢球在小轴间角下具有相似的传递载荷ꎬ但在大轴间角下的载荷分布有显著不同ꎮ因此需
要分析最大轴间角度下各关键部位之间的相互作用力ꎬ以下将分析钢球与钟形壳ꎬ以及钢球与保持架窗口的相互作用力ꎮ
4 2 1㊀钢球与钟形壳的相互作用力万向节大部分驱动扭矩通过钢球传递ꎮ每个钢球
从星形套接收驱动扭矩并将其传递到钟形壳ꎮ将轴间角固定在40ʎꎬ轴逆时针旋转一圈ꎬ驱动扭矩保持在
5782Nm不变ꎬ图8所示为万向节的等效应力ꎮ建立圆柱坐标系(见图8)ꎬ定义该坐标系中的θ
逆时针方向为正ꎬ钢球与钟形壳端面最远距离处为基准零点(图8中的2a处)ꎮ在θ=30ʎ位置时星形套的压痕最明显ꎬ在θ=260ʎ位置时钟形壳的压痕最大ꎬ而在其他角度处的应力值相对较小ꎮ由此可知ꎬ万向节在大轴间角下(40ʎ)ꎬ钢球和星形套㊁钟形壳不同位置的接触应力值具有显著的波动性ꎬ在θ=30ʎ处星形套的应力值最大ꎬθ=260ʎ处钟形壳的应力最大
图8㊀万向节在40ʎ轴间角下的等效应力
Fig.8㊀Equivalentstressofuniversaljointat40degreeAnglebetweenaxes
图9所示为钟形壳滚道与钢球之间相互作用的等效应力ꎮ其中ꎬ切向载荷Fθ是传递驱动力矩的驱动载荷ꎬ轴向载荷Fz是平行于驱动轴的反作用载荷ꎬ径向载荷Fr是沿径向的反作用载荷ꎮ从图9中可以更加
㊀166㊀机㊀㊀械㊀㊀强㊀㊀度2022年㊀
清晰的看出ꎬ钟形壳滚道上的应力峰值分别在θ=30ʎ和θ=260ʎ处ꎬ但是在θ=120ʎ和θ=180ʎ之间几乎不受力ꎮ
4 2 2㊀保持架窗口与钢球的相互作用力图10所示为保持架窗口与钢球之间相互作用的等效应力ꎬ保持架应力在θ=120ʎ处达到最大值ꎬ并且只有在θ=70ʎ和θ=190ʎ之间才有明显的应力值ꎬ该范围大致与图9所示钢球在钟形壳滚道上空载区域一致ꎮ由此可知ꎬ保持架受到载荷的根本原因在于钢球传递载荷的不平衡
图9㊀钟形壳滚道与钢球之间相互作用的等效应力Fig.9㊀Equivalentstressoftheinteractionbetweenouterrace
racewayandsteelball
4 3㊀大轴间角下保持架的凸出
图11所示为最大轴间角下万向节各个部件的位置与旋转半径ꎮ在θ=180ʎ附近ꎬ星形套与保持架向钟形壳外端面凸出ꎬ此时保持架窗口受到来自钢球㊁星形套和钟形壳共同作用下的剪切荷载ꎬ这种凸出间隙是造成保持架断裂的重要原因ꎮ
根据几何关系可得保持架的凸出间隙计算公式如下
[7]199 ̄205
g=Ror+Rir()sinarctancosθtanα()-βor+βir2éëêêùû
úú(5)
式中ꎬg是间隙ꎬRor是星形套外球面的半径ꎬRir是钟形壳内球面的半径ꎬβor是钟形壳中心到端面的角度ꎬ
βir是星形套中心到端面的角度ꎬα为轴间角ꎮ根据计算得到ꎬ保持架的凸出间隙从θ=150ʎ到
θ=210ʎ位置处变为正(如果为负ꎬ表示没凸出间隙
)ꎬ
图10㊀保持架窗口与钢球之间相互作用的等效应力Fig.10㊀Theequivalentstressoftheinteractionbetweenthe
cagewindowandthe
ball
图11㊀星形套与保持架向钟形壳外端面的凸出间隙Fig.11㊀Protrudingclearanceofinnerraceandcagetowardthe
outerendofouterrace
㊀㊀
这与图10所示的保持架应力仿真结果吻合ꎬ即保持架应力在θ=120ʎ时达到峰值ꎬ然后从θ=150
ʎ至θ=180ʎ附近急剧降低到接近0ꎬ整个圆周承受载荷极不均匀ꎮ
因此ꎬ在大轴间角下保持架的凸出间隙过大时ꎬ万向节保持架在恶劣工况下有可能因受力不均而断裂ꎮ因此ꎬ在设计阶段降低保持架的凸出间隙对于提高万向节的强度具有重要意义ꎬ然而过分限制保持架的凸出间隙又会导致万向节的轴间角范围受限ꎬ需要根据实际情况综合考虑合理的取值ꎮ
5㊀极限轴间角下静扭实验
进行极限轴间角下静扭实验ꎬ实验中保持轴间角在40ʎ不变ꎬ固定钟形壳ꎬ然后在驱动轴上施加逐渐变大的扭矩直到万向节破坏ꎬ试验结果如表3所示ꎮ