接触问题分析

《接触问题分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西安交通大学研究生考试卷专业班级硕2040学号3112045018姓名吴甲中成绩考核课程高等计算力学考核日期2012.11考核方式考试基于ABAQUS/Standard接触问题分析及实例摘要接触问题是许多工程实践中的常见问题，其实际结构系统往往由几个非永久性连接在一起的部分组成。这些参与接触之间的部分间会有沿接触面法向的相互作用（如接触压力）和沿接触面切向的相互作用（如摩擦作用）。在有限元分析中，接触条件是一类特殊的不连续约束，它允许力从模型的一部分传递到另一部分。因为只有当两个表面发生接触时才会有约束产生，当两个面分开时，就不存在约束作用，所以这种约束作用是不连续

2、的。本文将通过分析ABAQUS/Standard对接触问题的求解模式，来探讨有限元软件在求解接触问题时的内涵，并通过分析一个冲压金属板的实例来展现更为详尽的过程。一.ABAQUS/Standard中接触问题在ABAQUS/Standard中，接触问题或是基于表面（surface）或是基于接触单元（contactelement）。因此，首先必须在Interaction模块中各模型上创建可能发生接触的表面并判断哪一对表面可能具有接触约束，即接触对，随后定义控制各接触面之间相互作用的本构模型，这些接触面相互作用的定义包括诸如摩擦行为等。这里，接触问题属于边界非线性问题，边

3、界条件不再是定解条件，而是待求结果；两接触体间接触面积与压力随外载荷的变化而变，并与接触体的刚性有关。这是该问题的特点，也是困难所在。二.接触面间的相互作用1.接触面间的法向作用两个接触面分开的距离称为间隙（clearance），当间隙变为零时，表明两个表面形成接触关系（并不意味着接触约束的形成）。令P为两个接触面之间的接触作用力，当P为零或负值时，接触面分开即接触约束被移开。当P为正值时，表明接触约束形成（如图1）。由于接触条件从开（间隙值为正）到闭（间隙值为零）时接触压力可能剧烈变化导致这一过程存在剧烈非线性，因而在Standard模块中需要更多的增量步迭代以求

4、接触约束变化的过程收敛。图1接触约束形成条件2.常见的接触面间切向作用—摩擦模型当接触表面接触约束形成时，除了法向的接触压力外，还有阻止表面之间相对滑动的切向摩擦力。在实际的有限元软件分析过程中，关于相对滑动尺度的表征是很难确定的，从而极大增加了计算成本，小滑动的模型计算成本是很小的。在ABAQUS中，小滑移的定义遵循一个基本原则：相对滑移量不超过一个单元典型尺度，即可以近似应用“小滑动”来减小计算的复杂度。切向的摩擦作用通常用库伦摩擦（Coulombfriction）来描述，即应用摩擦因子和接触压力P的乘积来表示临界剪切摩擦力：当接触面之间的剪切作用等于时，接触面

5、间的相对滑移才会发生。如图2所示，在剪切作用达到之前的状态称为粘结状态，此时其间剪切作用小于，相对滑移量为零。图2弹性滑移摩擦模型然而,在这种理想化的摩擦模型中，粘结和滑移状态的不连续可能会导致ABAUQS/Standard中的收敛问题。因此ABAQUS在大多数默认情况下使用一个允许“弹性滑动”（弹性滑动认为在粘结的接触面之间允许发生小量的相对运动）的罚摩擦公式，并且ABAQUS会自动地根据所使用网格单元的尺度自动的选择罚刚度（图2），因此这里的弹性滑移量是单元尺度的小部分。三.ABAQUS的接触算法1.接触模拟的单元选择接触算法的关键在于选择合适的单元类型，而单元

6、类型的选择又与从面节点（slavesurface）的作用力有关。在3D单元中，二阶单元的等效结点载荷可能容易引起混淆，因为对于常值接触压力，它们甚至连符号都不相同。另外，对于图3所示的二阶单元的接触模型，算法将难以区分等效节点载荷和表示常值接触压力还是在单元面上的实际变化。而一阶单元却不存在上述问题，因此对于接触问题，应当尽量选择一阶单元进行网格划分。图3表面积为A的二维二阶单元的节点等效载荷1.接触算法的计算过程在ABAQUS/Standard中，接触算法的具体流程如下图所示：图4ABAQUS/Standard接触算法流程图这里接触算法是基于Newton-Raph

7、son迭代的算法所建立的。在每个增量步开始时，ABAQUS/Standard首先检查所有接触相互作用的状态，以建立从属节点是开放还是闭合立，若为闭合，还需进一步确定是处于滑动还是粘结。对于闭合的节点，ABAQUS/Standard施加一个约束，而对从闭合到开放状态的节点则撤销约束，然后进行迭代计算。任何节点若在迭代计算后其间隙变为非正值，则其接触状态从开放变为闭合；任何节点若在迭代计算后其接触压力变为负值，则其接触状态从闭合变为开放。在一次迭代结束后，ABAQUS/Standard会修正接触约束以反映接触状态的变化情况，并试图进行下一次迭代，直到完成迭代并不改变