接触问题有限元分析方法综述.pdfVIP

维普资讯 第 21卷第 3期 水 利 水 电科 技 进 展 2001年 6月 接触问题有限元分析方法综述 孙林松 ，王德信 ，谢能刚 (1．河海大学土木工程学院．江苏 南京 210098；2．华东冶金学院机槭工程系．安徽 马鞍山 243000) 摘要：从直接遮代法、接触约束法和数学规划法等方面综述接触问题有限元分析的基本方法．直接迭 代法是一种 “试验一误差”方法，概念清楚，实施方便，但计算工作量较大，而且不能保证迭代一定收敛． 接触约束法主要利用罚函数方法或 Iatgia~e乘子法将接触问题转化为无约束问题求解 数学规划法 利用接触问题的互补条件、非穿透条件等，将其归结为二次规划(线性~-41)问题求解，这是一种非迭 代类解法，收敛平稳、迅速，计算工作量较小．文末简单介绍接触问题研究的一些新趋势． 关键词：接触问题；有限元；迭代法；接触约束法；数学规划法 中图分类号：03433 文献标识码：A 文章编号：1006-7647f2001)03-00l8一O3 接触问题广泛存在于机械工程、土术工程等领 有限元混合法在弹性接触问题的求解中也得到 域，如齿轮的啮合、坝体的接缝等．这类 问题的特点 较广泛的应用．它以结点位移和接触力为未知量，并 是具有单边约束和未知接触区域，接烛区域的确定 采用有限元形函数插值，将接触区域的位移约束条 依赖于加载方式、荷载水平、接触面性质等因素，属 件和接触力约束条件均反映到刚度矩阵中去，构成 于边界待定问题 ． 有限元混合法控制方程[3】 接触问题 的研究很早就引起 了人们的重视 ．早 )={) 在 1882年，H．Herz就比较系统地研究了弹性体的接 触 问题 ，并提 出经典的Herz接触理论 ．随着数值解 其中，置，Ⅱ，，分别为通常有限元位移法中的剐度矩 法的兴起和发展，出现了许多求解接触 问题的非经 阵、未知结点位移 向量和结点荷载 向量；．，，工，c分 典方法，有限单元法作为解决复杂工程 问题的最有 别为接触力约束矩阵、位移约束矩阵和接触力向量． 效的数值方法，也成为求解接触问题的一种主要方 由于求解式 (1)中非对称方程组将多耗机时、 法．以有限元为基础的接触问题数值解法，主要可分 多 占内存 ，不少学者在对称化算法上进行 了研究L4j． 为直接迭代法、接触约束算法和数学规划法等． 对弹塑性接触问题，在求解过程中接触非线性和 材料非线性都需要迭代求解 ．通常是利用系统刚度矩 1 直接迭代法 阵的变化来反映材料非线性的影响，在每次塑性修正 迭代法是解决非线性问题的常用方法，它在接 迭代过程中都要结合对接触状态的判断进行接触迭 触问题的研究中也首先得到了应用．在用有限元位 代计算．拟弹性叠加双重迭代法 模拟弹性叠加原理 移法求解接触问题时，首先假设初始接触状态形成 建立有限元方程，再利用内外循环迭代求出方程的近 系统刚度矩阵，求得位移和接触力后，根据接触条件 似解，在整个计算过程 中不改变接触体刚度阵，仅增 不断修改接触状态，重新形成刚度矩阵求解，反复迭 加平衡力修正项，在一定程度上提高了计算效率 代直至收敛…．在上述方法中，每次迭代都要重新形 在用迭代法求解接触问题，尤其是有摩擦接触问 成刚度矩阵，求解控制方程，而实际上接触问题的非 题时．为了保证收敛到正确的结果，荷载增量往往要受 线性主要反映在接触边界上，因此，通常采用静力凝 到一定的限制 ．对弹塑性接触问题，荷载增量更是受 聚技术，使得每次迭代只是对接触点进行 ，大大提高