弹塑性接触分析（精编版）

1、题1： 表面光滑的刚性圆柱体与弹性平面的接触问题。有以下假设：接触体材料均匀连续，各向同性， 在接触区内只产生服从虎克定律的弹性变形，接触区相比接触体表面很小且在其附近的表面是光滑的，压力垂直于物体接触面，接触面上的摩擦力忽略不计。各参数为：计算区域宽度为 L=0.128mm ，圆柱体半径 R=0.5mm ，弹性模量 E=210GPa，泊松比，平面应变问题， P=50N/m, =0.31) 用有限元法求弹性平面应力分布；2) 用有限元法求的弹性平面表面接触压力分布曲线，并与Hertz理论解作对比。解：1、使用有限元方法求解（1） 建立有限元模型图 1 有限元模型如图 1 有限元模型， 刚性圆弧

2、半径为0.5mm ，AB 边长为 0.128mm 。可变形体采用PLANE42单元，如图 2 设置为处理平面应变问题。材料参数为： 弹性模量 E=210000MPa ，泊松比0.3。图 2 PLANE42的单元设置（2） 接触对设置按照图 3 所示的各图完成接触对的设置；在接触对的设置过程中，将圆弧线定义为刚体， 同时在坐标原点y 方向上 0.1mm 处定义刚体的控制节点，利用此节点施加刚体的边界条件； 选择图 1 所示的 AB 边作为可变形体的接触区域；最后使用翻转法线方向的命令，保证两接触对的法线方向相对。最后进行模型检测，看间隙是否过大，在接触单元Options中选择cnof/icont

3、 中选闭合Gap。接触算法采用软件默认的设置，不定义摩擦系数。9 /9图 3 设置接触对（3） 施加边界条件如图 4 所示施加边界条件。约束可变形平面底边的所有自由度，约束刚体控制点x 方向的位移，并在刚体控制点上施加负y 方向 50N 的压力。图 4 施加边界条件（4） 计算结果进行求解，获得的两接触对的接触压力如图5 所示，最大接触压力值为2685MPa，位于加载的中心。 可变形体内部的Mises 应力分布如图6 所示，最大 Mises 应力值为1665MPa， 位于接触区域以下。 （中间变细， 弹性模量变成平面应变模量。contact presure2853.von mises 1742

4、）图 5 接触压力 / 赫兹接触应力对应（单位：MPa）应力是内部材料抵抗变形而产生能量的反映，压力是作为力传递的外载在表面材料上的分布图 6 Mises 应力图（单位：MPa）2、使用解析法求解（选用 mm 、MPa 计算结果与选用m、 Pa计算的结果不同） 根据 Hertz 接触理论，接触半宽的计算公式为：a4PR E其中： a 为接触半宽；P 为外载荷；R 为刚性圆弧的半径；E 是可变形体的等效弹性模量。EE / (12 )接触区域的接触压力分布为：第一主应力沿深度方向的分布：p( x)2P ( a2a 2x2 )1/2x第二主应力沿深度方向的分布：p0 ( a2a2z2 )(a 2z2

5、 )1/ 22zzp0a(a 2z2 )1/ 2主剪应力沿深度方向的分布：1p0a zz2 ( a 2z2 )1/ 2 计算得到的最大接触压力为2700MPa，压力分布如图7 所示： 由 Ansys 计算得到的接触压力分布如图8 所示，最大接触压力为2682MPa（ 2720MPa），误差为 0.67%。300025002000150010005000-0.015-0.01-0.00500.0050.010.015图 7 由 Hertz 公式计算得到的接触压力分布General Postproc>>path operations>>define path>>

6、map on paths>>plot path item.图 8 由 Ansys 计算得到的接触压力分布（峰值数据和云图数据对不上,调整好数据点数就可以调出峰值）题 2：结合自己的专业与研究方向，选择一三维接触问题进行建模、求解及计算分析。要求有问题的提出、建模、求解、分析及结论各部分的详细论述。解：行星滚柱丝杠副基本结构如图9所示，丝杠 A是具有螺旋升角的多头螺纹；螺母B具有内螺纹，牙型为三角形，牙型半角与滚柱牙型半角相同；C是滚柱，单线螺纹，为了增大接触 面积， 其螺旋面通常加工出具有较大接触半径的圆形轮廓，滚柱与螺母的螺旋角相同，以保证滚柱在螺母内滚动时无相对轴向位移。同时为

7、了消除丝杠螺旋升角对滚柱产生的倾斜力矩，在滚柱两端加工有直齿D，与内齿圈 E啮合，以确保滚柱轴线平行于丝杠轴线而正常滚动。F为滚柱保持架， 使滚柱沿圆周均匀分布；滚柱保持架则由弹簧挡圈G定位； H为密封圈， I是定位销， J是润滑油孔。图9 行星滚柱丝杠结构行星滚柱丝杠的接触分析是对其后续进行摩擦、润滑和磨损分析的基础。所以有必要对其接触状态进行分析，故本文使用ABAQUS软件进行行星滚柱丝杠的单螺纹啮合计算。为了节约计算时间只进行单螺纹啮合计算，对模型进行简化。 创建的行星滚柱丝杠的单螺纹啮合几何模型如图10所示。由于该行星滚柱丝杠具有5个滚柱，所以取螺母以及丝杠的1/5进行分析。为了减少网

8、格的划分个数去掉多余的螺纹部分。对行星滚柱丝杠的网格划分结果如图11所示。采用 C3D8R（8节点线性三维积分缩减单元）以及 C3D4（4节点线性四面体单元）进行网格的划分。如图12在螺纹啮合的部分进行网格的细化，共划分了403702个单元。图10 行星滚柱丝杠单螺纹啮合几何模型图11 行星滚柱丝杠单螺纹啮合有限元模型图 12 接触区域的网格细化如图 13，在可能发生接触的区域创建接触对，同时为了方便施加载荷，在螺母以及丝杠轴线的方向上创建两个参考点，同时建立点对面的耦合关系。如图14约束螺母、 滚柱、丝杠的x方向以及 z方向的位移，在A点施加载荷，约束B点所有的自由度。图13 接触对以及自由

9、度耦合的创建图14 边界调节与施加外载荷有限元的计算结果如图1518所示。行星滚柱丝杠螺母螺纹接触部位的位置如图15（ a） 所示，接触椭圆的中心位于滚柱轴线和螺母轴线所形成的平面内。螺母接触部位的局部放大 如图 15（ b）所示， 应力剖面如图 15（ c）所示。 滚柱与螺母接触部位的应力分布如图16所示， 滚柱与丝杠接触部位的应力分布如图17所示，可以看出滚柱与螺母接触椭圆的中心以及滚柱的圆心和滚柱与丝杠接触部位的中心并不像大多数文献假设的位于同一条直线上。滚柱与丝杠接触部位的中心偏离了中心线。丝杠接触部位的应力分布如图18所示， 模型中各个零件的装配位置关系是按照实际给出的，但是丝杠螺纹的最初接触位置并没有位于螺纹节线附近，而处于螺纹的顶端，这种受力情况不利于行星滚柱丝杠使用