曙批㈨ - journal of electronic science and technology.pdfVIP

第32卷第4期 北京工业大学学报 V01．32No．4 0F UNIVERSITY0FTECHNOLOGY 2006年4月 JOURNALBEUING Apr．2006 结构工程中接触问题的数值计算方法 韩 青1， 张毅刚1， 赵凯红2 摘要：接触问题涉及到应力集中、边界非线性、材料或几何非线性等问题．随着计算机的发展，有限元与边界 元理论和数值计算方法都取得了很大进展．本文对建筑工程中接触问题的数值计算方法进行了综述，分析了各 种方法的优缺点及适用范围，并对建筑工程中接触问题进一步的研究方向提出建议． 关键词：迭代法；边界元法；有限元法；接触问题；接触约束法；数学规划法 311．4 中图分类号：Tu 文献标识码：A 接触问题涉及到应力集中、边界非线性以及材料或几何非线性等问题．物体之间的接触可能对称也 可能非对称，可能弹性接触也可能弹塑性接触，可能有摩擦接触也可能无摩擦接触，可能是一维接触也可 能是多维接触，接触状态也分黏接、分离、滑动等情况．用解析的方法分析接触问题只局限于简单的问题， 而在工程中遇到的实际接触问题，其结构形状、边界形状及接触状态都非常复杂，通常无法求得解析解． 随着计算机技术的发展，出现了有限元、边界元等数值计算方法，使得进一步分析工程中遇到的复杂问题 成为可能[1·4|． 在结构分析中，有限元法是分析接触问题有效通用的方法，以有限元为基础的接触问题数值解法，又 可分为直接迭代法、接触约束算法和数学规划法等． 1直接迭代法 迭代法是解决非线性问题的常用方法，在接触问题中首先得到了应用．该方法首先对小载荷下接触 区域作初始接触状态估计，根据接触条件对工作载荷下的接触状态不断地迭代、修改直至收敛．在计算 中，每次迭代都要重新形成刚度矩阵，求解控制方程，工作量较大，因此，采用静力凝聚技术，使得每次迭代 只是对接触点进行，提高了求解效率．另外，还有虚力法[5]，用沿边界的虚拟等效压力来模拟接触状态，这 样在每次迭代中并不重新形成刚度矩阵，所做的只是回代工作． 有限元混合法在弹性接触问题的求解中也得到较广泛的应用，它以结点位移和接触力为未知量，采用 有限元形函数插值，将接触区域的位移约束条件和接触力约束条件均反映到刚度矩阵中，构成有限元混合 法控制方程 曙批㈨) (1) 式中：K为刚度矩阵；“为未知结点位移向量；，为结点荷载向量；J为接触力约束矩阵；L为位移约束矩 阵；c为接触力向量． 对弹塑性接触问题，在求解过程中接触非线性和材料非线性都需要迭代求解．通常利用系统刚度矩 阵的变化来反映材料非线性的影响，在每次塑性修正迭代过程中要结合对接触状态的判断进行接触迭代 计算． 在用迭代法求解接触问题，尤其是有摩擦接触问题时，为了保证收敛得到正确的结果，荷载增量要受 收稿日期：2004—11—11． 作者简介：韩青(1965一)，女，天津人，副教授 万方数据 322 北京工业大学学报 到一定的限制‘6|．对弹塑性接触问题，荷载增量不允许在一个增量步中出现2种非线性的限制【5|．用迭 代法计算非对称接触问题时，将耗费较多机时，占用较多内存空间，故该方法较适用于对称接触分析． 2接触约束算法 接触问题可描述为求区域内位移场u，使得系统的势能玎(u)在接触边界条件的约束下达到最小， 即 min Ⅱ(【厂)=告【，[T1Ku一【，[T1F s．￡． g0 (2) 接触约束算法是通过对接触边界约束条件的适当处理，将式(2)的约束优化问题转化为无约束优化问 题求解，可分为罚函数方法和Lagrange乘子法等． 1)罚函数方法 罚函数方法是将接触非线性问题转化为材料非线性问题．根据处理方法的不同又可分为障碍函数法 和惩