模拟管道中任意裂纹扩展的基于CB壳单元的扩展有限元方法的研究.pdfVIP

模拟管道中任意裂纹扩展的基于CB 壳 单元的扩展有限元方法研究 * * * *1) 许丹丹 ，王恒 ，成斌斌 ，庄茁 * （清华大学航天航空学院，北京 100084） 摘要：本文将扩充形函数引入到基于连续体的 （continuum-based, CB ）壳单元中，以模拟曲面壳体中沿任意方向 扩展的裂纹。扩充形函数可以代表独立于网格的间断。计算三维裂纹应力强度因子时，应用“等效域积分”方法， 采用最大能量释放率准则作为裂纹扩展判据。该方法能够获得管道中独立于单元网格的任意裂纹扩展路径。由于 裂纹与网格独立，使得在裂纹扩展过程中无需对网格进行重新划分。文中给出了管道中不同断裂模式的数值算例。 关键词：CB 壳单元；扩展有限元；三维任意裂纹扩展；能量释放率。 1 引言 由于管道在工程中的广泛应用以及因裂纹非稳态扩展而引发的灾难性后果，管道中的断裂问题已引起 高度重视，人们已开展了大量相关研究工作。数值模拟在研究管道裂纹止裂方面发挥了重要作用。许多运 [1] [2] 用有限元思想的方法已被用于处理裂纹扩展问题，如节点力释放 和移动有限单元 等。近年来，扩充形 函数的概念很大的丰富了管道断裂的数值研究。此概念使用独立于网格的能够描述间断和奇异性的扩充形 函数，从而能对管道裂纹扩展进行更加有效而精确的数值模拟。扩展有限元方法首先由Belytschko 等人于 [3,4] [5,6] [7] [8] 1999 年提出 ，并在裂纹扩展 ，剪切带演化 和多相流 等不同的间断问题上得到令人满意的结果。使 用该方法能够获得管道中独立于单元网格的任意裂纹扩展路径。由于裂纹与网格独立，使得在裂纹扩展过 程中无需对网格进行重新划分。 本文将扩充形函数的概念引入到基于连续体的壳单元中，以模拟管道中沿任意方向扩展的裂纹。结构 中使用8 节点CB 壳单元，不仅可以用于模拟曲面壳体，还可以有效缓解自锁现象。在数值计算中，CB 壳 单元可以考虑壳体厚度的变化，这一点在研究厚壳问题时不容忽视。CB 壳单元有限元格式由三维实体单 元推导得出，扩充形函数也同样建立于六面体单元，通过使用类似于标准 CB 壳单元的方法将扩充自由度 引入壳单元中。计算三维非平面裂纹应力强度因子时，本文应用 “等效域积分”方法，以计算壳体中的断 裂参数。文中采用最大能量释放率准则作为裂纹扩展判据。 在第2 部分中介绍了扩展有限元方法的基本思想，并在第3 部分介绍建立于CB 壳单元的扩展有限元 格式。第4 部分介绍了此方法中的完全的拉格朗日格式。断裂应力强度因子的计算及壳体中的裂纹扩展准 则将在第5 部分中加以阐述。第6 部分则给出了运用这种算法计算的不同断裂形式下的数值算例。 1Email: zhuangz@tsinghua.edu.cn 2 扩展有限元基本思想 扩展有限元方法采用扩充形函数的办法来表示计算域内的间断。将一些具有不连续性质的形函数引入 [3] 到位移场的构造中。经过离散，域内试函数表达式为 ： u N (X)u   (X)a i  I iI  J iJ 1 ()