0滚滑接触下钢轨热力耦合分析
- 格式:pdf
- 大小:788.75 KB
- 文档页数:7
第27卷第8期 Vol.27 No.8 工 程 力 学 2010年 8 月 Aug. 2010 ENGINEERING MECHANICS 199
——————————————— 收稿日期:2009-02-19;修改日期:2009-05-11 基金项目:国家自然科学基金项目(50705078, 50821063);国家重点基础研究发展计划项目(2007CB714702);霍英东教育基金会青年教师基金项目(114028);四川省青年科技基金项目(08ZQ026-021) 作者简介:李 伟(1985―),男,湖北监利人,硕士生,从事轮轨滚动接触疲劳研究(E-mail: 1022liwei@163.com); 温泽峰(1976―),男(壮族),广西南宁人,副研究员,博士,从事轮轨接触力学研究(E-mail: zfwen@home.swjtu.edu.cn); 吴 磊(1981―),男,贵州毕节人,博士生,主要从事轮轨热力耦合研究(E-mail: leiwulanxia@)163.com); *金学松(1956―),男,江苏扬州人,教授,博士,博导,从事轮轨关系研究(E-mail: xsjin@home.swjtu.edu.cn).
文章编号:1000-4750(2010)08-0199-06 滚滑接触下钢轨热力耦合分析
李 伟,温泽峰,吴 磊,*金学松 (西南交通大学牵引动力国家重点实验室,成都 610031)
摘 要:采用热弹塑性有限元法,对滚滑接触状态下的钢轨进行热力耦合分析,计算钢轨接触区附近的温升和应力/应变,并分析蠕滑率对它们的影响。有限元模型中,考虑材料参数随温度的变化,通过在钢轨表面上移动边界条件来模拟车轮的移动。分析结果表明:热影响区主要在很薄的钢轨表层;钢轨表层温升和热应变均随着蠕滑率的增大而增大。热载荷在钢轨表层产生的周向和轴向残余应力表现为拉应力,其在蠕滑率小时有减小机械载荷产生的残余压应力作用;在蠕滑率较大时热响应对残余应力和残余应变影响明显,且当蠕滑率增大到一定值时,考虑热影响时热力耦合作用下的周向残余应力和轴向残余应力都表现为拉应力,而不考虑热影响时二者表现为压应力,随着蠕滑率的继续增大,考虑热影响时的残余应力和残余应变随之增大。 关键词:轮轨滚滑接触;蠕滑率;热力耦合;温升;残余应力/应变;有限元法 中图分类号:TH117.1; U213.4 文献标识码:A
THERMO-MECHANICAL COUPLING ANALYSIS OF RAIL IN ROLLING-SLIDING CONTACT
LI Wei , WEN Ze-feng , WU Lei , *JIN Xue-song (State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China)
Abstract: A thermo-mechanical coupling analysis of the rail under a rolling-sliding state is conducted by the thermo-elasto-plastic finite element method. The temperature rise and the stress/strain near the wheel-rail contact patch of the rail are calculated. The effect of the wheel-rail creepage is also investigated. The temperature-dependent material properties are taken into consideration in the finite element model. The movement of boundary conditions of the rail surface is used to simulate the movement of wheels. The results indicate that the thermally affected zone exists mainly in a very thin layer of the material near the rail contact surface. Both the temperature rise and thermal strain increase with increasing the creepage. The circumferential and axial residual stresses generated by the thermal load in the surface layer of the rail appear to be tensile. The thermal load caused by small creepages can reduce the residual compressive stresses generated by the mechanical load. When the creepage becomes larger, however, the thermal response has a significant influence on the residual stresses and residual strains. As the creepage reaches a certain value, the circumferential and axial residual stresses in the thermo-mechanical case become tensile stresses, while in the mechanical case they appear to be compressive stresses. With further increasing the creepage, the residual stresses and residual strains in the thermo-mechanical case increase as well. 200 工 程 力 学 Key words: wheel-rail rolling-sliding contact; creepage; thermo-mechanical coupling; temperature rise; residual stress/strain; finite element method
热损伤是高速重载车轮和钢轨主要失效形式之一[1―2]。研究轮轨在各种运行状况下的温度场和
应力/应变场分布状况与特点以及各种因素的影响是揭示轮轨热损伤机理和研究轮轨损伤减缓措施的基础。很多研究者对轮轨摩擦热问题作了许多的研究,大多数只涉及摩擦温升的分析[3―12],而对摩
擦热应力的分析相对较少。文献[13]建立平面应变有限元模型对轮轨摩擦热弹性应力进行了计算,得到各应力分布情况,模型中没有考虑机械载荷(轮重和切向力)的作用。文献[14]利用有限元软件ABAQUS计算了无限体平面的温度、应力和塑性应变的分布情况,模型中没有考虑温度对材料参数的影响。文献[15―17]利用有限元法,建立了轮轨摩擦耦合热弹性有限元分析模型,分析了轮轨摩擦温升、应变和应力分布状况及特点,但都没有考虑塑性影响和温度变化对材料参数的影响。文献[18]采用弹塑性有限元法,计算了多步非稳态载荷下的钢轨滚动接触应力和变形,给出了钢轨接触表面的残余应力和残余应变情况,但忽略了温度的影响。本文在文献[15―17]的基础上,基于热力完全耦合方法,利用有限元软件ABAQUS建立了轮轨接触热力耦合有限元分析模型,模型中考虑轮轨自由表面与环境间的热对流影响和温度对材料参数的影响。分析了蠕滑率对温升、残余应力和残余应变的影响,并与不考虑热影响的情况进行了对比分析。
1 数值计算模型
1.1 轮轨接触模型 轮轨经一段时间磨损和变形后,轮轨接触斑变得细长,可将轮轨接触视为二维问题[19]。轮轨接触斑上的法向压力分布由Hertz接触理论给出[20]:
202
()1xpxp
a=−
(1)
式中:a表示椭圆接触斑宽度之半;p0为接触斑上的最大接触压力。 考虑轮轨接触斑处于滚滑状态,则接触斑分为滑动区和粘着区。在粘着区,接触质点对具有相同的滑动速度。在滑动区,接触质点对发生相对滑动,切向力正比于法向压力。Carter提出了粘滑区切向力分布的计算表达式[19]:
(), ()()(),ssqxaxcdqxqxqxcdxa
′−−+⎧
=⎨
′′′+−+
⎩
≤≤≤≤ (2)
式中: ()()qxfpx′= (3)
1/22
0()1sxdc
qxfp
ac
⎧⎫−⎪⎪⎛⎞
′′=−−
⎨⎬⎜⎟
⎝⎠⎪⎪⎩⎭
(4)
式(4)中:f为摩擦系数;c为粘着区宽度之半,且ds =a−c。图1是式(2)的示意图。
图1 切向力分布 Fig.1 Tangential force distribution
文献[20]指出,当两个接触物体可用半空间近似且具有相同的弹性常数时,切向力对接触斑的形状和大小以及法向压力分布没有影响。本文假设车轮和钢轨具有相同的弹性常数,且接触斑大小远小于轮轨的几何特征尺寸和接触斑附近的曲率半径,因此本文假设的轮轨接触模型中忽略了切向力对法向力的影响。文献[21]分析认为温升不大时,热变形小,热对法向力和切向力分布的影响可以忽略,因而本文计算时忽略了温升对轮轨接触法向力和切向力的影响。 这里还假定车轮在平直的直线轨道上运行,即轮轨作用力大小不随时间变化的影响,关于动态载荷作用下的轮轨接触数值模拟可参考文献[18],本文主要研究轮轨蠕滑率对轮轨接触热响应的影响,动态载荷作用下的轮轨接触热响应分析将另文 讨论。 式(2)中,粘着区的大小由下式确定:
()qx′()qx
()qx′′a cs
d
滑动区 粘着区