地下工程中的非线性接触算法研究及数值实现.pdfVIP

第33 卷 第 12 期 岩石力学与工程学报 Vol.33 No.12 2014 年12 月 Chinese Journal of Rock Mechanics and Engineering Dec .，2014 地下工程中的非线性接触算法研究及数值实现 1 1，2 1 周墨臻 ，钱晓翔 ，张丙印 (1. 清华大学 水沙科学与水利水电工程国家重点实验室，北京 100084；2. 北京市轨道交通设计研究院有限公司，北京 100089) 摘要：基于罚函数法的点对面形式因其具有显著的简便性而在工程领域内应用广泛。针对罚因子可能引起矩阵病 态的缺陷，采用局部坐标系下的相对位移作为与接触相关的广义自由度，大幅提高了罚因子的取值上限，且无需 使用专门的矩阵预处理子即可用于求解大规模接触问题。针对点对面存在投影奇异性及结果振荡的风险，发展采 用径向基点插值法的三维曲面光滑方法，其优势在于网格适应性极强且不损失曲面精度。将双向链接边表作为三 维接触面的数据结构，有效地降低了接触邻接搜索的计算代价。基于以上改进，开发大型非线性接触分析有限元 计算程序，数值算例的结果表明，该计算程序具有较强地描述位移不连续现象的能力，能有效求解接触非线性方 程，可用于模拟大规模地下工程中的复杂多体接触问题。 关键词：地下工程；接触；有限元；罚函数法；点对面 中图分类号：TU 91 文献标识码：A 文章编号：1000 –6915(2014)12 –2390 –06 ALGORITHM AND NUMERICAL ANALYSIS METHOD FOR NONLINEAR CONTACT PROBLEMS IN UNDERGROUND ENGINEERING 1 1，2 1 ZHOU Mozhen ，QIAN Xiaoxiang ，ZHANG Bingyin (1. State Key Laboratory of Hydroscience and Engineering ，Tsinghua University，Beijing 100084，China；2. Beijing Rail and Transit Design and Research Institute Company Limited，Beijing 100089，China) Abstract ：The node-to-segment approach using penalty method was widely applied in engineering field due to its great simplicity. To deal with the problem in which a large penalty factor may cause ill-condition，the degree of freedom relevant to contact is described by relative displacement in local coordinate system，which enlarges significantly the upper-limit of penalty factor and provides an efficient solution for large-scale contact problem without using any special preconditions. A general 3D contact smoothing method based on radial point int