三维近不可压弹性问题的数值方法及其快速求解算法研究-一般力学与力学基础专业论文.docxVIP

万方数据 万方数据 湘潭大学 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所 取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任 何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡 献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的 法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意 学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文 被查阅和借阅。本人授权湘潭大学可以将本学位论文的全部或部分内容编 入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇 编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 III 湘潭大学硕士学位论文 湘潭大学硕士学位论文 摘 要 有限元(FEM)方法是求解三维弹性力学问题的一类重要的数值方法。在实际 计算时，诸如橡胶、塑料等材料呈现出近不可压缩（即泊松比 ν ??0.5 ）的性质， 利用通常的有限元（如线性元）进行求解会出现所谓的体积闭锁现象，往往需要 采用一些特殊的方法。本文，首先基于 ANSYS 平台系统研究了六面体网格剖分 下高阶单元法、减缩积分法及基于 u / p 格式的混合高阶元法对求解混合边界条件 的三维近不可压缩问题的有效性和鲁棒性（robustness）。数值结果表明：这三种 协调有限元法均能有效地克服三维弹性材料的体积闭锁现象，其中混合高阶元法 最为精确，计算所得位移值均随网格尺寸变小而稳定地收敛于理论解。但混合元 方法得到的总体刚度矩阵通常为一半正定矩阵，且计算规模比位移法大一倍，在 选取快速求解器将带来不便。我们希望对近不可压缩问题进行离散后得到的是一 个对称、正定矩阵，这样便于选取更为有效的求解器。本文第二部分针对混合边 界条件的三维近不可压缩弹性问题，利用基于能量泛函极小的罚函数有限元方法 克服体积闭锁现象，详细推导了相应的计算格式，分析了该方法实施成功的条件， 并通过数值实验验证了该方法对解决体积闭锁现象的有效性和鲁棒性。在三维有 限元分析中，剖分网格的质量将对计算精度和求解效率产生很大影响，实际计算 时若能采用各向同性网格，则对问题的分析将具有更好的收敛性。本文，最后针 对罚函数有限元分析中形成的大型的、稀疏的和高度病态的正定方程组，设计了 几种预处理共轭梯度（PCG）法，包括基于块对角逆预条件子的 PCG 法（即 M1 -PCG 和 M 2 -PCG）和基于整体矩阵的 RS-GAMG-PCG 法，并对这些方法在 求解悬臂梁问题和 Cook 膜问题罚函数二次元方程的计算效率进行了数值测试与 分析，结果表明，对近不可压缩弹性问题，若能利用容易获知的部分几何与分析 信息（如方程类型，节点自由度信息），再结合经典 AMG 法中的网格粗化技术 及插值算子构造方法，可设计具有更好计算效率的 AMG 法及相应的 PCG 法， 将大大提高其有限元分析整体效率。 关键词：三维弹性问题；有限元方法；Locking 现象；罚函数有限元；病态矩阵； 代数多重网格法；预条件子 I Abstract The finite element (FEM) method is an important class of numerical methods for solving the three-dimensional elastic problems. In practical calculation，such as rubber and plastic, which show nearly incompressibility material property, i.e., Poissons ratio ν is close to 0.5. The volume locking phenomenon will happen when the commonly used finite elements (such as a linear element) are applied to the solution of this nearly incompressibility problems. Thus, we need to use some special methods. In this paper, first of all, based on the ANSYS software, we make some systematic studies for the efficiency and robustness of high