General Formulas and Program Design for Manifold Method Based on Independent Covers Ⅱ: General Program Design

SU Hai-dong, YANG Zhen, XIE Zhi-qiang, QI Yong-feng, GONG Ya-qi

Journal of Changjiang River Scientific Research Institute ›› 2025, Vol. 42 ›› Issue (4) : 202-210.

PDF(6477 KB)
PDF(6477 KB)
Journal of Changjiang River Scientific Research Institute ›› 2025, Vol. 42 ›› Issue (4) : 202-210. DOI: 10.11988/ckyyb.20240113
Numerical Manifold Method Based on Independent Covers

General Formulas and Program Design for Manifold Method Based on Independent Covers Ⅱ: General Program Design

Author information +
History +

Abstract

Based on the general calculation formula of the manifold method based on independent covers presented in the previous article, we provide the flowchart of the calculation program. First, we summarize the integration methods for various geometric shapes (such as partitions, stripes, and boundary faces) that may appear in one- to three-dimensional spaces. On this basis, we develop integration programs according to simplex geometric elements of points, lines, faces, and bodies. This approach ensures the universality for any mesh shape. Next, we propose a programming strategy that separates the integration module from the integrand function module. The arbitrary combination of these two modules endows the program with extensibility and the potential to achieve universality in solving partial differential equations. Moreover, the universality of series is realized through the determination of series formulas, corresponding coordinates, coordinate transformation matrices, and series matrices. In addition, all calculation parameters can be input via formulas using user subroutines, thus achieving universality of input parameters. Ultimately, with relatively less program code, we can conduct one- to three- dimensional steady-state and transient analyses of the differential equations of motion in elasticity, conduction equations, and wave equations, including one to three types of boundary conditions.

Key words

partial differential equations / series solutions / mesh division / exact geometry / independent covers / numerical manifold method

Cite this article

Download Citations
SU Hai-dong , YANG Zhen , XIE Zhi-qiang , et al . General Formulas and Program Design for Manifold Method Based on Independent Covers Ⅱ: General Program Design[J]. Journal of Changjiang River Scientific Research Institute. 2025, 42(4): 202-210 https://doi.org/10.11988/ckyyb.20240113

References

[1]
苏海东. 独立覆盖流形法的通用计算公式和通用程序设计:(一)通用计算公式[J]. 长江科学院院报, 2025, 42(4):193-201,210.
(SU Hai-dong. General Formulas and Program Design for Manifold Method Based on Independent Covers Ⅰ. General Formulas[J]. Journal of Changjiang River Scientific Research Institute, 2025, 42(4):193-201,210.(in Chinese )
[2]
SHI G H. Simplex Integration for Manifold Method, FEM, DDA and Analytical Analysis[C]//Proceedings of the First International Forum on Discontinuous Deformation Analysis and Simulations of Discontinuous Media. Berkeley, Califonia, USA, 1996:205-262.
[3]
LIN S, XIE Z. A New Recursive Formula for Integration of Polynomial over Simplex[J]. Applied Mathematics and Computation, 2020, 376: 125140.
[4]
ZIENKIEWICZ O C, TAYLOR R L. The Finite Element Method[M].5th Edition. Oxford, UK: Butterworth-Heinemann, 2000.
[5]
苏海东, 付志, 颉志强. 基于任意形状网格和精确几何边界的数值计算[J]. 长江科学院院报, 2020, 37(7):167-174.
Abstract
有限元网格形状要尽可能规则,网格之间必须通过结点连接,这些要求给复杂形状求解域的数值计算带来很大的前处理工作负担,而且实际的曲线边界一般要离散成有限单元能够描述的形式,难以模拟CAD模型的精确几何。针对这些问题,基于独立覆盖流形法提出任意形状且任意连接的覆盖网格,在CAE分析中模拟CAD模型的精确几何边界及其边界条件:将求解域划分为可包含曲线边的任意形状的块体网格,可以采用单纯形解析积分和数值积分2种方式进行块体积分;仅需在积分过程中考虑块体之间的窄条形(包括曲线条)的覆盖重叠区域,而不必在计算模型中生成这些条形;通过边界条实现本质边界条件的严格施加,包括曲线上的边界条件;给出2个数值算例验证了方法的有效性。任意形状的覆盖网格将为实现基于精确几何模型的数值计算及其完全自动化的前处理开辟新的路径。
(SU Hai-dong, FU Zhi, XIE Zhi-qiang. Numerical Computations Based on Cover Meshes with Arbitrary Shapes and on Exact Geometric Boundaries[J]. Journal of Yangtze River Scientific Research Institute, 2020, 37(7): 167-174.) (in Chinese)
Finite element meshes should keep regular shape as much as possible, and ensure correct connections through nodes. These requirements pose a great burden to the pre-processing procedure of numerical computations for solving domains with complex shapes. On the other hand, curve boundaries in practical situations are usually discretized into shapes which finite element meshes can describe, resulting in an imprecise simulation of exact geometry defined in CAD. In view of this, cover meshes with arbitrary shapes and arbitrary connections are implemented using Manifold Method based on independent covers. Exact geometric boundaries of CAD models and boundary conditions are simulated in CAE analyses. The solving domain is divided into block meshes with arbitrary shapes which can contain curve boundaries. And two approaches, including analytical integration method with simplexes and numerical integration method, can be used for the block integration. The thin strips for cover overlapping are considered only in the integration process, but are not necessarily involved in the generation of computation models. Essential boundary conditions are strictly applied through boundary strips, including the boundary conditions on curves. Moreover, two numerical examples are given to illustrate the validity of the method. Cover meshes with arbitrary shapes bring about a new path for numerical computations based on exact geometric models and automatic pre-processing procedures.
[6]
SU H D, QI Y F, GONG Y Q, et al. Preliminary Research of Numerical Manifold Method Based on Partly Overlapping Rectangular Covers[C]//DDA Commission of International Society for Rock Mechanics. Proceedings of the 11th International Conference on Analysis of Discontinuous Deformation (ICADD11). Fukuoka, Japan, August 27-29, 2013, London: Taylor &Francis Group, 2013: 341-347.
[7]
祁勇峰, 苏海东, 龚亚琦. 基于新型流形法的三维应力强度因子求解[J]. 土木与环境工程学报(中英文), 2021, 43(5):58-65.
(QI Yong-feng, SU Hai-dong, GONG Ya-qi. Computing 3D Stress Intensity Factors Based on New Manifold Method[J]. Journal of Civil and Environmental Engineering, 2021, 43(5):58-65.) (in Chinese)
[8]
苏海东, 宋文硕, 龚亚琦, 等. 独立覆盖流形法的通用计算公式和通用程序设计:(三)算例验证[J]. 长江科学院院报, 2025, 42(4):211-218.
(SU Hai-dong, SONG Wen-shuo, GONG Ya-qi, et al. General Formulas and Program Design for Manifold Method Based on Independent Covers Ⅲ: Example Verification[J]. Journal of Changjiang River Scientific Research Institute, 2025, 42(4):211-218.) (in Chinese)
[9]
SMITH I M, GRIFFITHS D V. 有限元方法编程[M]. 王崧, 周坚鑫, 王来, 等, . 北京: 电子工业出版社, 2003.
(SMITH I M, GRIFFITHS D V. Programming the Finite Element Method[M]. Translated by WANGSong, ZHOUJian-xin, WANGLai, et al. Beijing: Publishing House of Electronics Industry, 2003.) (in Chinese)
PDF(6477 KB)

Accesses

Citation

Detail

Sections
Recommended

/