接触问题求解方法刘怀举的博客.docVIP

接触问题求解方法 (2010-10-02 00:25:00)转载▼ 标签： 杂谈 不管在接触边界之间是否有间隙存在，接触作用的出现对结构受载荷之后的接触状态和应力分布都有直接的影响，有些结构正是由于接触作用，使不连续的部分共同工作，从而提高了整个结构的承载力和刚度；而正是由于接触的存在使得有些结构出现局部高应力，很容易使材料屈服或发生裂缝，如果再受到循环荷载的影响，还可能产生疲劳失效。因此对于接触问题的研究具有重要的工程实际意义。 赫兹接触理论及后来其它学者发展的弹性接触理论称为经典接触力学，它们都是通过封闭的解析解来解决接触问题的，但其应用范围有限，因而接触力学的进展主要与消除这些限制有关。接触问题属于数学上的混合边值问题，Boussinesq积分方程是其主导方程。按所用数学方法的不同，可将接触问题的理论解法大致分为经典解法和非经典解法。以经典的数学工具如积分变换法和复变函数法求解接触问题的方法称为经典解法，以有限元法、边界元法等求解接触问题的方法称为非经典法。以传统有限元法等数值方法为基础的非经典解法主要有罚函数法、拉格朗日乘子法、增广拉格朗日乘子法、摄动拉格朗日乘子法及数学规划法。 Trefftz有限元模型利用辅助网线位移场或面力场，在一种杂交意义上将单元域内位移场关联起来。单元域内位移场精确满足控制微分方程，它可表达为微分方程的特解、适当截断的Trefftz完备解系与待定参数乘积的和的形式。利用定义在每个单元边界上的独立的网线位移场，单元间的连续性就在一种近似意义上得到满足。在单元一级上消去内部待定参数后即可得到标准的力-位移关系式（即单元刚度方程）。 变分泛函是Trefftz有限元法的核心，它在单元公式推导中起着至关重要的作用。 由于Trefftz有限元法继承了传统有限元法和传统边界元法的优点：(1)公式中只含有对单元的边界积分，这样就可以生成任意多边形单元甚至曲边单元。因此，Trefftz有限元法可被认为是一种特殊形式的边界型求解方法，单元的边界类似于一种特殊形式的边界单元，其刚度矩阵对称、计算简单，而不像传统边界元那样要进行复杂的计算（如计算边界奇异积分的复杂积分规则，处理非齐次方程的特殊积分，处理间断问题的双节点技术等等）。(2)Trefftz有限元在处理各种几何、荷载奇异性或局部效应（如角点、孔洞、裂纹、夹杂、集中载荷等）问题有较高的效率和许多优点。因为只需将插值函数取成既满足控制方程，又满足给定的边界条件即可，而这种特殊函数可容易地从许多文献中找到。这样，在遇到角点、孔洞、裂纹、夹杂以及由各种载荷引起的应力集中或各种奇异性时，不需局部网格细分，只需在插值函数上作调整，就能达到所希望的精度。 （以上From弹性接触问题的杂交Trefftz有限元解法 王克用 天津大学博士论文 2006） 用有限元法求解不可逆接触问题，与求解同等规模的非接触问题相比，计算机内存及计算时间大大增加，为克服此缺点，先后提出了各种求解方法，包括子结构法、广义子结构法及混合有限元法。这些方法都是根据接触力与位移的条件建立有限元方程求解，对接触边界条件的处理都是基于一般力学的方法。 对于接触边界条件处理方面最优成果的工作是将接触问题考虑为有约束的最小值问题，对这种最小值问题的一种解法就是罚函数法，将接触区域的非嵌入条件作为罚项引入到接触系统的总势能中，将约束变分问题转化为罚优化问题求解。 Wriggers和Zavarise，Heegaard和Curnier，Simo和Laursen由变分不等式出发建立了带有库伦摩擦/大滑动/非线性本构关系的接触问题的罚有限元公式，即增广拉格朗日法，并可由变分不等式讨论其解的存在性和唯一性。增广拉格朗日法既改善收敛性，同时还改善了罚函数法的精度，成为工程界流行的求解摩擦接触问题的算法。这种方法的精度强烈依赖于所用罚因子的数量，但是各作者均采用相应的方法使得适当的罚因子大小能精确满足各种条件。目前一些商业软件如Ansys，Abaqus等都可以选用此方法。 接触约束算法就是通过对接触边界的适当处理，将约束优化问题转化为无约束优化问题求解。根据无约束优化方法的不同，一般可分为罚函数法、拉格朗日乘子法、增广拉格朗日乘子法和直接约束法。 罚函数法实际上是将接触非线性问题转化为材料非线性问题。根据处理方法不同又可分为障碍函数法和惩罚函数法。障碍函数法假设接触面之间充满某虚拟物质，在未接触时其刚度趋于零，不影响物体的自由运动，在接触区域其刚度变得足够大，能阻止接触物体的相互嵌入。罚函数法对接触约束条件的处理是通过在势能泛函中增加一个惩罚势能，从而将接触问题等价为无约束优化问题。其 原理是一旦接触区域发生穿透，罚函数便夸大这种误差的影响，从而使系统的求解无法正常进行。罚函数法不增加系统求解规模，但由于人为