滚轮接触应力的有限元法分析
- 格式:docx
- 大小:38.59 KB
- 文档页数:3
滚轮接触应力的有限元法分析
杨阳
【摘 要】滚轮应用领域广泛,它承载力大,转动速度慢,滚轮与轨道接触变形大,接触状态复杂.通过有限元计算完成滚轮的应力应变计算,取得应力应变分布数据.对于接触而导致的有限元计算结果不收敛问题,采用赫兹公式计算验证,力学性能满足要求.
【期刊名称】《煤矿机电》
【年(卷),期】2017(000)003
【总页数】3页(P76-78)
【关键词】滚轮;接触应力;有限元分析(FEA)
【作 者】杨阳
【作者单位】煤炭科学技术研究院有限公司,北京100013;煤炭资源高效开采与洁净利用国家重点实验室,北京100013
【正文语种】中 文
【中图分类】TB301.1
滚轮是一种重要的机械支撑传动部件,应用广泛。它是在低速重载下工作的回转支承[1],与轨道的接触状态复杂。由于试验台体积和质量都很大,需要计算滚轮和滚圈轨道间的比压,保证不被压溃。对滚轮的力学能力进行精确计算和数值模拟是设计的重点和难点。
回转支承的失效形式主要有两种[2]:一是滚道损坏,二是断齿。而滚道损坏占的比例达98%以上,评判材料是否产生塑性屈服失效的一个重要准则是Mises屈服准则。
本项目通过有限元软件计算方式完成滚轮的应力应变及分布情况计算,计算取得应力分布等数据,对于接触而导致的有限元计算结果不收敛问题,采用赫兹公式计算验证力学性能是否满足要求。
滚轮材料为40CrNi,调质处理后,屈服极限σs=785×106 Pa,密度7.85 g/cm3,质量0.94 t,弹性模量206 GPa,泊松比0.3。滚轮直径1 m,导轨直径13.6 m,承载设备质量512 t,滚轮接触面直径1 000 mm,接触面宽200 mm。将建立好的模型部件转化为Parasolid格式,导入有限元软件对模型进行力学分析计算。滚轮模型如图1所示。
为了得到便于计算的动态网格特性,缩减计算时间,边界定义过程中把滚轮分成两部分:滚轮主体与接触部分的两个区域,如图2所示。
网格设定的参数将决定网格的大小、单元形状、中间节点位置等结果,将影响分析结果准确性和过程的经济性[3]。
有限元软件中网格划分方法分为自由网格划分和映射网格划分[4]。对于滚轮主体部分采用自由网格划分形式,网格尺寸为2 mm,四面体结构网格;接触部分网格尺寸大小设定为1 mm,并采用映射网格划分方式,以提高计算精度,保证计算结果具有更好的收敛性[5],滚轮网格划分如图3所示。节点数量为500 310个,单元数量为307 750个。
1)计算参数设定。滚轮工作情况下,承受工作载荷为32 kN。
2)约束与载荷设定。根据实际受力情况,约束与载荷分别如图4、图5所示,对滚轮内环采用fixed support约束,对滚轮外接触部分施加force载荷力320 kN。
3)有限元计算。如图6、图7所示,最大应力为123 MPa,在滚轮与轨道接触点,最大变形量为0.059 4 mm,应力及应变成对称规律分布。
由弹性力学可知,当轴线方向平行的两个圆柱体相接触并受压力作用时,接触面将发生弹性形变,由线接触变为面接触,其接触面为一段矩形接触面,在接触面上产生接触应力,并且最大接触应力位于接触区中线上,对于接触而导致的计算结果不收敛问题,采用赫兹公式计算验证,其计算理论值符合赫兹公式的要求:
式中:P——使两物体接触的法向力; ρ1,ρ2——曲率半径; E1,E2——两接触物体的弹性模量; μ1,μ2——两接触物体的泊松比; b——接触长度。
中正号用于外接触,负号用于内接触,实际是外侧接触,取正号。根据上述条件,公式简化为:
式中,ρ为综合曲率半径。
σmax=0.418=79.4 MPa
最后计算安全系数,通过有限元计算835/123≈6.79;赫兹公式计算:835/79.4≈10.5。
本项目通过有限元软件计算完成滚轮的应力应变计算及分布情况,取得应力分布等数据。利用材料力学计算对滚轮进行载荷分析,验证有限元计算结果,使误差结果在合理范围内,安全系数在要求的范围内。滚轮能满足力学性能,保证了设备的安全使用,为大型试验设备、矿山机械、翻车机、大型研磨机、起重机的支撑滚轮的设计提供了新的方法和理论基础。
【相关文献】
[1] 侯宁.影响回转支承承载能力的四个参数[J].建筑机械,2002(1):21-22
[2] 徐立民,陈卓.回转支承[M].合肥:安徽科学技术出版社,1988.
[3] 王国强.实用工程数值模拟技术及其在ANSYS上的实践[M].西安:西北工业大学出版社,1999:1-76.
[4] 易日.使用ANSYS6.1进行结构力学分析[M].北京:北京大学出版社,2002:101.
[5] Xie G, Ramaeker J A H.Graded mesh generation and transformation[J]. Finite Elements
in Analysis and Design,1994,17:41-55.