具体的螺旋桨轴和衬套过盈联结弹性力学计算结果如表1所示。
2.2 过盈联结的有限元模拟计算 应用ANSYS 有限元软件强大的后置处理功能,求解过盈联接后的接触应力分布与变形状况。
接触问题中产生接触的两物体须满足边界不穿透约束条件,在接触边界施加不穿透约束的方法主要有拉格朗日乘子法、罚函数法、增广拉格朗日乘子法和基于求解器的直接约束法。前两种方法处理时都具有局限性,直接约束法处理接触问题是追踪物体的运动轨迹。一旦探测出发生接触,便将接触所需的运动约束和节点力作为边界条件直接施加在产生接触的节点上,对接触的描述精度高,具有普遍适应性,因此本文有限元计算均采用基于求解器的直接约束法求解约束应力问题。为了和过盈联结联接组件的工况接近, 在有限元分析中须合理确定联接组件的相互作用关系、约束条件[6~8]。
在ANSYS软件中接触类型可分为两种:刚体—柔体接触和柔体—柔体接触。在刚体—柔体接触中,两接触面中有一个被当作刚体,与柔体相比它有大得多的刚度;柔体—柔体接触类型则假定两接触体均为变形体,变形量的大小主要由材料的弹性模量和结构刚性而定,且要求两接触体都剖分网格,这正好满足了研究需要,由此可确定选择柔体—柔体接触类型。
目前ANSYS软件可支持的接触方式有三种:点—点、点—面和面—面接触。1) 点—点接触方式,主要用于模拟点—点的接触行为,也可以模拟面—面接触行为,不过使用这种接触方式需要预先知道确切的接触位置,且只适于模拟接触面间有较小相对滑动的情况,计算量较大。2) 点—面接触方式,通过一组节点来定义接触面,生成多个接触单元,也可以模拟面—面接触问题,使用这类接触单元的优点在于:不需要预先知道确切的接触位置,接触面间的网格也可以不一致。3) 面—面接触方式,可弥补点—面接触方式生成单元太多的不足,计算量相对较小,同时面—面接触单元还可提供更好的接触结果,如法向压力。因此本文采用面—面接触方式模拟衬套过盈联结的接触行为。
基于整体接触法的思路,过盈联接组件的结构接触分析可按以下步骤进行:
1) 定义接触分析单元Contact174,定义过盈联接组件的材料属性并指定单元为轴对称;
2) 建立过盈联接组件的1/4作为分析模型;
3) 采用扫略方式分网;
4) 使用接触向导在过盈联接组件上定义面—面接触单元,其中将衬套内表面(目标面) 和轴外表面(接触面) 确定为一对接触副;
5) 添加位移约束等载荷,对模型进行求解;
6) 计算结果分析,着重分析过盈联接组件的接触压力和应力分布状况。
下面给出轴和衬套长度为L = 60 mm,轴半径为b= 50mm。过盈量δ=0.02mm, 衬套外半径c=70mm的有限元数值计算具体结果云图,图2、3、4分别表示衬套、轴的综合应力分布云图和接触面的接触应力分布云图。
由应力分布云图知,衬套的综合应力远大于轴的综合应力,螺旋桨轴和衬套的综合应力均出现在其两端边缘,即两端是危险截面。
为了对螺旋桨轴和衬套过盈联结进行比较全面的研究,本文分别对不同模式进行了详细的弹性力学和有限元计算分析,为节约篇幅,这里仅给出有代表性的结果,见表1所列。
表1 弹性力学和有限元计算主要结果
Table 1 The main results of elasticity force and finite element calculation method 模式 过盈量δ/mm 衬套外半径
c/mm 衬套的最大
综合应力σemax/MPa 螺旋桨轴的最大综合
应力σemax/MPa 螺旋桨轴和衬套间
接触应力σH/MPa 弹性力学法 有限元法 弹性力学法 有限元法 弹性力学法 有限元法 模式1 0.02 70 74.7 75.2 19.4 19.9 19.6 20.7 模式2 0.035 70 130.4 131 34.1 34.7 35.3 36.3 模式3 0.05 70 186.6 188 50.1 50 50.6 51.7 模式4 0.05 100 181 182.5 77.1 77.6 77.4 79 模式5 0.05 150 179.6 181.2 91.4 91.8 91.8 93.3
上表是给出五种模式的主要尺寸,然后分别采用弹性力学方法和有限元方法对螺旋桨轴和衬套过盈应力进行分析计算。值得注意的是,基于弹性力学理论计算的结果,是假定过盈联结面上的压力均匀分布,没有考虑应力沿轴向的变化情况,也没有考虑轴和衬套两端边缘应力的变化,而工程实际中轴和衬套过盈联结的边缘应力问题必须加以考虑,另外,在传统方法中过盈联接的计算、选用、校核,是通过第三强度理论的当量弯矩法,对强度的校核结果较为局限,无法深入地了解材料的应变的部分和范围,尤其对于应力接近或超出材料屈服极限情况下的联结难以做出准确的判断,因此,该方法应力计算结果是不符合实际且不可靠的。ANSYS 有限元软件有强大的后置处理功能,考虑了过盈配合边缘效应,对过盈联接后应力的分布、变形的范围进行非常直观的分析,其该方法更接近实际情况,由模式1、2、3计算结果得出,过盈量不同,轴和衬套的综合应力也不相同,过盈量越大,螺旋桨轴和衬套的综合应力也越大;由模式3、4、5得出,衬套外半径不同,轴和衬套的综合应力也不相同,衬套外半径越大,轴和衬套的综合应力越小;衬套的最大综合应力远大于轴的最大综合应力。 3 结论 1) 从表1的结果对比知,有限元法求解的应力值略大于弹性力学方法,主要是因为有限元法考虑了过盈联结的边缘效应,其该方法更接近实际情况。
2) 衬套的综合应力远大于轴的综合应力。
3) 螺旋桨轴和衬套过盈联结过盈量越大,螺旋桨轴和衬套的综合应力都越大;衬套外半径越大,螺旋桨轴和衬套的综合应力都越大;
4) 螺旋桨轴和衬套的综合应力均出现在其两端边缘,即两端是危险截面。
参考文献:
[1] 张松,艾兴. 高速主轴过盈联结特性研究[J]. 制造技术与机床,2003,(10).
[2] 章巧芳,贾虹,. 轮轴过盈联接有限元分析[J]. 机械强度,2006,(4).
[3] 严国平. 船用艉轴机械密封热-结构特性数值分析:[博士学位论文],湖北:武汉理工大学: 2008.
[4] HILLS D A, NOWELL D,SACKFIELD,Mechanics of ElasticContact[M].Butterworth Heinemann Lt,1993.
[5] JOHNSON K L.ContactMechanics[M].Combridge Cambridge University Press,1985.
[6] 许小强,赵洪伦. 过盈联结应力的接触非线性有限元分析[J]. 机械设计与研究,2000,(1).
[7] 王勖成,邵敏. 有限单元法基本原理和数值方法(第二版). 北京:清华大学出版社,1997.3,64~68
[8] Peter Kohnke,PH.D.ANSYS,Inc.Theory Manual.001369.Twelfth Edition. SAS.IP,Inc. 1994.