Continental ContiTech马牌人字齿同步带有限元法计算同步带弹性接触受力和弹性接触有限元理论
Continental ContiTech马牌人字齿同步带在啮合传动时与带齿的接触状态非常复杂，受到带齿面弹性变形，几何非线性和材料非线性等多种因素的影响，使得计算带齿与轮齿的接触应力变得相当困难。传统做法是用以两平行圆柱对压的赫兹理论计算接触强度，赫兹理论做了许多的假设，仅假设相接触的两物体只产生弹性变形，应力应变关系符合胡克定律这一条假设，就与同步带复杂的材料特性非常不符合。因此传统赫兹理论不足以精确描述问题，计算带齿接触强度时的精度和可靠度都很低。
有限元法是分析Continental ContiTech马牌人字齿同步带接触问题的最有效的方法之一，特别是对于复杂的接触问题，若用解析法求得精确结果十分困难，甚至不可能实现，只能采用数值近似的模拟方法。有限元法具有快速、精确和可靠性高等优势，此使用弹性接触有限元法分析弹性带齿的受力是接近真实情况并得到较精确解的最有效途径之一。人字齿同步带传动过程中带齿与带轮齿槽的啮合接触，可以看作是两个弹性物体相互接触。由前面对带齿的受力分析可以看出，带齿受到紧边拉力、松边拉力、载荷、齿顶的摩擦力等几种力的综合作用。带在拉力作用下由紧边开始进入啮合接触，带齿与轮齿相互接触后，接触面积的大小和接触处的应力都会随着两物体间合力的变化而变化，这些变化不仅与接触面上合力的大小有关，而且与带齿和带轮齿各自材料的性质有关。
2.3.1 弹性接触有限元理论
使用弹性有限元的基本理论，可以有效解决复杂的同步带的接触问题。图 2-9
表示 A 和 B 两个弹性物体相互接触，分别对物体 A、B 划分网格，形成了由节点组
成的能够计算的有限元网格，两物体初始接触面上的节点组成了多对接触点对。采用三维八节点接触面单元如图 2-10 所示，假设接触点对之间由无数微小的弹簧建立连接关系，接触面之间靠节点处的作用力联系起来，单元形心放置在原点处，z 向
表示接触压力，x、y 表示摩擦切向力。
图 2-9 两物体的弹性接触状态 图 2-10 三维八节点接触单元
物体 A 和 B 的接触状态有三种，分别是：分离、粘结接触和滑动接触[37]。对于这三种情况，接触点对的位移和力的条件各不相同，实际的接触状态在这三种接触情况中转化，导致了接触问题高度非线性特点。 PA 和 PB 表示两物体受到的外载荷， uA 和 uB 表示两物体的位移，可以得到以下方程：
|
|
|
|
ïì{K A }{u A }= {P A }-{R A} |
(2-38) |
|
|
|
|
í |
|
|
|
|
|
ï{K B }{u B }= {P B }-{RB } |
|
|
|
|
|
î |
|
|
式中 |
{K A }，{K B }——弹性体 A、B 的刚度矩阵； |
|
||
|
|
u A |
， u B |
——弹性体 A、B 的节点位移向量； |
|
|
|
{ } |
{ } |
|
|
{P A }，{P B }——A、B 的整体载荷向量；
{R A }，{R B }——A、B 的接触力向量。
两物体处于不同接触状态下时的边界条件不一样，设 unA 、 unB 分别为接触点对的法向位移，utA 、utB 为切向位移，RnA 、RnB 分别为接触点对的法向接触力，RtA 、RtB
为切向接触力，可得到不同接触状态的边界条件如下[38] ：
分离接触点对：
RnA = RnB = RtA = RtB = 0
粘结接触点对：
RnA + RnB = 0 ， RtA + RtB = 0 ， unA = unB ， utA = utB
滑动接触点对：
RnA + RnB = 0 ， RtA = -mRtB ， unA = unB ， utA = utB