基于任意网格的二维自适应分析及其优化

宋文硕, 苏海东, 颉志强

raybet体育在线 院报 ›› 2025, Vol. 42 ›› Issue (12) : 160-169.

PDF(2408 KB)
PDF(2408 KB)
raybet体育在线 院报 ›› 2025, Vol. 42 ›› Issue (12) : 160-169. DOI: 10.11988/ckyyb.20241090
水工结构与材料

基于任意网格的二维自适应分析及其优化

作者信息 +

Two-dimensional Adaptive Analysis Based on Arbitrary Meshes and Its Optimization

Author information +
文章历史 +

摘要

独立覆盖流形法采取“分区级数解”的思想,将计算域划分为若干分区,在各个分区内部采用多项式级数等完备级数直接逼近真实物理场函数;各分区之间采用窄条带连接,通过单位分解函数将各分区级数连接成整体近似函数,在收敛意义上具有C1连续性;计算网格具有任意形状、任意连接和任意加密的特性。基于计算网格的任意性,采用“一分为二”的网格分裂算法作为网格划分和细分的策略;基于收敛意义上的整体C1连续性,采用物理场导数的连续程度作为误差控制的指标。以上组成了基于任意网格的自适应分析策略。通过数值试验验证了这一策略的可行性。在此基础上提出优化方案:采用绝对误差指标控制近似函数导数的连续性,以简化误差判断;采用网格预划分的方式优化凹角局部区域的网格分布。最后,通过方孔、重力坝等算例进行验证,表明优化方案可以大幅减少网格数量,从而节省算力。

Abstract

[Objective] This study aims to optimize the two-dimensional adaptive analysis strategy of the independent cover-based manifold method, focusing on addressing its deficiencies in error control and mesh distribution, thereby significantly enhancing computational efficiency and engineering practicality. [Methods] Based on the arbitrarily shaped and connected cover meshes of the independent cover-based manifold method, a “split-one-into-two” mesh splitting algorithm was employed for arbitrary refinement, and the degree of continuity of physical field derivatives was adopted as the error control indicator, forming an adaptive analysis strategy. An optimization scheme was proposed. 1) Adopting an absolute error indicator to replace the relative error indicator: the original relative error indicator tended to cause over-refinement in regions of minor stress and was overly sensitive in concave corner singularity regions. Using the absolute error indicator not only simplified the error judgment logic but also permitted larger error thresholds to be set near singular points such as concave corners, thereby effectively avoiding over-refinement. 2) Introducing a local mesh pre-partitioning and short strip elimination strategy: to address the issue of excessively high mesh density and irregular distribution in concave corner regions, a local pre-partitioning strategy was proposed, which pre-set the initial mesh in these regions by inwardly offsetting and reversely extending the edges of the concave corner. Simultaneously, an adjacent point merging algorithm was introduced during the mesh splitting process, which avoided the generation of extremely short connection strips and improved the conditioning of the system equations. [Results] Verification through two typical hydraulic structure examples, the square-hole and the gravity-dam model, demonstrated that the optimized scheme achieved a breakthrough improvement in computational efficiency. For the square-hole example, the original adaptive strategy generated 310 covers, corresponding to 6 520 degrees of freedom (DOFs). Under the same accuracy objective, the optimized scheme required only 59 covers and 933 DOFs. This represented a reduction of approximately 81% in the number of meshes and approximately 86% in DOFs. For the gravity-dam example, the original strategy generated 228 covers and 4 810 DOFs, whereas the optimized scheme required only 106 covers and 2 354 DOFs, achieving significant results of over 53% reduction in the number of meshes and 51% reduction in DOFs. The most notable achievement of the optimized scheme was in the effective suppression of mesh over-refinement near concave corner singularity regions. The calculation results demonstrated that the new strategy could generate more reasonable meshes, while ensuring computational accuracy, it substantially reduced the computational scale, and greatly enhanced the computational efficiency. [Conclusion] The proposed optimization strategy significantly enhances the efficiency of adaptive analysis while maintaining high accuracy. Through absolute error control and local mesh pre-partitioning, it effectively solves the problems of mesh over-refinement and unreasonable distribution near concave corner singularities, laying a foundation for subsequent three-dimensional adaptive analysis and engineering applications. Future research includes: criteria for selecting error thresholds and the highest order of cover series; further automating the local mesh pre-partitioning process to enable it to handle more complex geometries, ultimately achieving the goal of efficient and fully automatic simulation analysis for hydraulic structures.

关键词

网格自动划分 / 自适应分析 / 误差分析 / 独立覆盖 / 数值流形方法

Key words

automatic mesh generation / adaptive analysis / error analysis / independent covers / numerical manifold method

引用本文

导出引用
宋文硕, 苏海东, 颉志强. 基于任意网格的二维自适应分析及其优化[J]. raybet体育在线 院报. 2025, 42(12): 160-169 https://doi.org/10.11988/ckyyb.20241090
SONG Wen-shuo, SU Hai-dong, XIE Zhi-qiang. Two-dimensional Adaptive Analysis Based on Arbitrary Meshes and Its Optimization[J]. Journal of Changjiang River Scientific Research Institute. 2025, 42(12): 160-169 https://doi.org/10.11988/ckyyb.20241090
中图分类号: O302 (力学中的数学方法)    TV3 (水工结构)   

参考文献

[1]
SU H, LIN S, XIE Z, et al. Fundamentals and Progress of the Manifold Method Based on Independent Covers[J]. Science China Technological Sciences, 2024, 67(4): 966-991.
[2]
郎美玲. 针对非线性椭圆型方程的自适应有限体积元法研究[D]. 济南: 济南大学, 2022.
(LANG Mei-ling. Adaptive Finite Volume Element Method for Nonlinear Elliptic Equations[D]. Jinan: University of Jinan, 2022. (in Chinese))
[3]
BRENNER S C, SCOTT L R. The Mathematical Theory of Finite Element Methods[M]. New York: Springer New York, 2002.
[4]
王勖成, 邵敏. 有限单元法基本原理和数值方法[M]. 2版. 北京: 清华大学出版社, 1997: 35-36.
(WANG Xu-cheng, SHAO Min. Basic Principle and Numerical Method of Finite Element Method[M]. 2nd Ed. Beijing: Tsinghua University Press, 1997: 35-36. (in Chinese))
[5]
ZIENKIEWICZ O C, TAYLOR R L, ZHU J Z. The Finite Element Method: Its Basis and Fundamentals[M]. Oxford,Tokyo: Butterworth-Heinemann, 2013: 43-79.
[6]
BELYTSCHKO T. 连续体和结构的非线性有限元[M]. 庄茁,译. 北京: 清华大学出版社, 2002: 4-11.
(BELYTSCHKO T. Nonlinear Finite Elements for Continua and Structures[M]. Translated by ZHUANG Zhuo. Beijing: Tsinghua University Press, 2002: 4-11. (in Chinese))
[7]
BABUŠKA I. Error Estimates for Adaptive Finite Element Computations[J]. SIAM Journal on Numerical Analysis, 1978, 15(4): 736-754.
[8]
BABUŠKA I, MILLER A. A Feedback Finite Element Method with a Posteriori Error Estimation: Part I. the Finite Element Method and Some Basic Properties of the a Posteriori Error Estimator[J]. Computer Methods in Applied Mechanics and Engineering, 1987, 61(1): 1-40.
[9]
BABUŠKA I, RHEINBOLDT W C. A Posteriori Error Estimates for the Finite Element Method[J]. International Journal for Numerical Methods in Engineering, 1978, 12(10): 1597-1615.
[10]
许雯. 二维椭圆偏微分方程基于高精度技术的自适应有限元方法[D]. 湘潭: 湘潭大学, 2021.
(XU Wen. High Accuracy Techniques Based Adaptive Finite Element Method for Two-dimensional Elliptic PDEs[D]. Xiangtan: Xiangtan University, 2021. (in Chinese))
[11]
WESSELING P. Principles of Computational Fluid Dynamics[M]. Heidelberg: Springer, 2001.
[12]
SHI G H. Manifold Method of Material Analysis[C]// Proceedings of the Ninth Army Conference on Applied Mathematics and Computing. Minneapolis, Minnesota, USA, 1991:51-76.
[13]
裴觉民. 数值流形方法与非连续变形分析[J]. 岩石力学与工程学报, 1997, 16(3):279-292.
(PEI Jue-min. Numerical Manifold Method and Discontinuous Deformation Analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 1997, 16(3):279-292. (in Chinese))
[14]
祁勇峰, 苏海东, 崔建华. 部分重叠覆盖的数值流形方法初步研究[J]. raybet体育在线 院报, 2013, 30(1): 65-70.
(QI Yong-feng, SU Hai-dong, CUI Jian-hua. Preliminary Study on Numerical Manifold Method with Partially Overlapping Covers[J]. Journal of Yangtze River Scientific Research Institute, 2013, 30(1): 65-70. (in Chinese))
Finite element meshes based on completely overlapping covers are currently used as mathematical meshes in Numerical Manifold Method (NMM). The decline of calculation accuracy at some key parts is caused by the mismatch between physical meshes and mathematical meshes. In view of this, NMM with partially overlapping covers is presented for the first time. The method  adopts an analysis mode mainly concerning independent covers, and the continuity between the independent covers is kept by a small partially overlapping region. Thus NMM with completely overlapping covers is extended to the method of covers in a general meaning. Rectangular covers are preliminarily researched. The forms and formulas of partially overlapping covers are derived conveniently from completely overlapping covers by using freedom constraint among nodes of a single rectangular mathematical mesh. Finally, the validity and accuracy of the new NMM are demonstrated by an example.    
[15]
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.
[16]
苏海东, 宋文硕, 龚亚琦, 等. 独立覆盖流形法的通用计算公式和通用程序设计: (三)算例验证[J]. raybet体育在线 院报, 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))
[17]
苏海东, 韩陆超, 颉志强. 精确几何薄曲梁曲壳分析的分区级数解[J]. raybet体育在线 院报, 2022, 39(9):144-151.
摘要
薄梁板壳的数值计算涉及关于挠度的4阶微分方程,其困难在于构造C<sub>1</sub>连续的近似函数;同时,由于薄曲梁和曲壳控制方程的复杂性,通常用直梁或平板单元近似地模拟曲梁或曲壳,容易产生几何误差进而带来力学分析上的误差。前期研究采用独立覆盖流形法实现了基于厚梁板壳假设的精确几何曲梁和曲壳分析,本文在此基础上讨论了这种新型流形法的分区级数解的C<sub>1</sub>连续性,完成了基于Euler-Bernoulli梁理论和Kirchhoff-Love板壳理论的精确几何薄曲梁和曲壳分析,并解决了几何公式推导复杂的问题。详细给出了薄曲梁的计算公式,简述了薄曲壳的计算过程,将前期文献中的算例在薄梁板壳假设下重新计算,验证了方法的有效性,相比厚梁板壳假设可节省约30%的自由度。研究成果同时展示了应用独立覆盖流形法求解4阶微分方程的潜力。
(SU Hai-dong, HAN Lu-chao, XIE Zhi-qiang. Analysis of Thin Curved Beam and Curved Shell with Exact Geometry Using Piecewise-defined Series Solutions[J]. Journal of Yangtze River Scientific Research Institute, 2022, 39(9): 144-151. (in Chinese))
The numerical calculation of thin beam, plate and shell involves the fourth-order differential equation about deflection whose difficulty lies in constructing approximation functions with <i>C</i><sub>1</sub> continuity. In the meantime, due to the complexity of the governing equation, the thin curved beam and curved shell are usually simulated approximately by using straight beam or flat plate elements, which is prone to generate geometric errors and then brings errors in mechanical analysis. In our previous study, manifold method based on independent covers is used to analyze curved beam and shell with exact geometry based on the assumption of thick beam and shell. On this basis, the <i>C</i><sub>1</sub> continuity of the piecewise-defined series solutions of the new manifold method is discussed. The thin curved beam and shell with exact geometry is analyzed based on Euler-Bernoulli beam theory and Kirchhoff-Love shell theory, and the complexity of derivation of geometric formula is overcome. The calculation formula of thin curved beam is given in detail, and the process of thin curved shell is briefly described. The examples in previous study are recalculated under the assumption of thin beam, plate and shell, which verifies the effectiveness of the proposed method. Compared with the assumption of thick beam, plate and shell, the method saves about 30% of the degree of freedom. Meanwhile, the research demonstrates the potential of solving the fourth-order differential equations by applying manifold method based on independent covers.
[18]
苏海东, 颉志强, 龚亚琦, 等. 基于独立覆盖的流形法的收敛性及覆盖网格特性[J]. raybet体育在线 院报, 2016, 33(2):131-136.
摘要
针对前期提出的基于部分重叠覆盖的数值流形方法,将其内涵范围缩小,仅研究其中的一种情况——基于独立覆盖的数值流形方法。从完备性和协调性2个方面讨论该方法的收敛性,特别强调其收敛性是基于各个独立覆盖的逼近而建立起来的,独立覆盖之间条形连接区域的尺寸要取小,并由此推断及用实例说明,覆盖网格可以具备“3个任意”的优良特性——任意形状、任意连接以及由此而来的可任意加密的能力,从而有望使数值计算的前处理工作大为简化。
(SU Hai-dong, XIE Zhi-qiang, GONG Ya-qi, et al. Characteristics of Convergence and Cover Mesh in Numerical Manifold Method Based on Independent Covers[J]. Journal of Yangtze River Scientific Research Institute, 2016, 33(2): 131-136. (in Chinese))
The scope of a numerical manifold method (NMM) based on partially overlapping covers is narrowed to a special case based on independent covers. Convergence of the new method is discussed from two aspects: completeness and coordination. The convergence of the method is due to the convergence of each independent cover. Results show that the size of the strips between independent covers should be small. Moreover, the cover meshes have three excellent features: arbitrary shape, arbitrary connection, and arbitrary refinement. Finally, some illustrations are given to verify these “arbitrary” features, and the method can be used to greatly simplify the pre-processing of numerical analysis.
[19]
吴顺唐. 数学分析[M]. 南京: 南京大学出版社,2022:29-46.
(WU Shun-tang. Mathematical Analysis[M]. Nanjing: Nanjing University Press, 2022: 29-46. (in Chinese))
[20]
苏海东, 杨震, 颉志强, 等. 独立覆盖流形法的通用计算公式和通用程序设计:(二)通用程序设计[J]. raybet体育在线 院报, 2025, 42(4): 202-210.
摘要
在独立覆盖流形法通用计算公式的基础上,给出了计算程序的整体流程。对一维至三维各种几何形体(包括分区、条带和边界面)的积分方式进行总结,基于点、线、面、体的单纯形几何元素开发积分程序,实现网格形状的通用性。提出将积分模块与被积函数模块分开考虑的编程思路,然后再将两者任意组合,使程序具备了扩展性,有望实现偏微分方程求解的通用性。通过级数公式和相应的各种坐标以及坐标转换矩阵、级数矩阵的确定,实现了级数的通用性。所有计算参数都可以通过用户子程序输入公式,实现输入参数的通用性。最终可用较少的程序代码,实现弹性力学运动微分方程、传导方程、波动方程的一维至三维稳态和瞬态分析(含一类至三类边界条件)。
(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. (in Chinese))
[21]
苏海东, 龚亚琦, 颉志强, 等. 基于矩形独立覆盖初步实现结构静力分析的自动计算[J]. raybet体育在线 院报, 2016, 33(2): 144-150.
摘要
采用前期研究的基于矩形独立覆盖的新型数值流形法,提出结构线弹性静力分析的自动计算方法,包括自动的前处理、自适应分析等。根据独立覆盖的特点提出几个后验误差指标:独立覆盖之间条形连接区域的应变连续性指标;边界应力指标和独立覆盖的高阶误差指标。利用新方法的h型网格加密及p型升阶的方便性,选择一种路径尝试h-p型的混合自适应,其中,对于矩形独立覆盖采用简单的二分法实现覆盖加密。通过几个二维算例验证了新方法实现自动计算的可行性,只需人工输入结构外形、材料参数和边界条件,其它工作完全交由计算机完成,最终得到满足一定精度的计算结果。
(SU Hai-dong, GONG Ya-qi, XIE Zhi-qiang, et al. Preliminary Implementation of Automatic Computation for Static Analysis of Structures Using NMM Based on Independent Rectangular Covers[J]. Journal of Yangtze River Scientific Research Institute, 2016, 33(2): 144-150. (in Chinese))
By means of numerical manifold method (NMM) based on independent rectangular covers proposed in previous study, we present an automatic computation method for static analysis of linear-elastic structures, including automatic pre-processing, self-adaptive analyses and so on. According to the characteristics of independent covers, we give 3 indexes for posterior error such as index of strain continuity in strip area between two covers, stress index on boundary surfaces and high-order error in an independent cover. By using convenient h-version mesh refinement and p-version order increasing in the new method, we implement h-p version self-adaptivity in a selected way to realize h-version refinement of rectangular covers by using simple bisection method. Some 2D numerical examples are given to illustrate the feasibility of automatic computation, in which all the procedures are automatically accomplished by the computer, except for necessary manual input of structural outlines, material parameters, and boundary conditions. Finally, we obtain calculated data with certain precision.
[22]
苏海东, 付志, 颉志强. 基于任意网格划分的二维自动计算[J]. raybet体育在线 院报, 2020, 37(7): 160-166.
摘要
针对有限元法网格剖分和加密不方便、难以实现精确几何建模以及人工操作量大等问题,采用前期提出的独立覆盖流形法,利用其覆盖网格所具有的任意形状、任意连接和任意加密的特性,基于“凸剖分”的思路提出二维求解域的一种任意网格划分方法。在此基础上,结合前期研究的误差估计和h-p型混合自适应分析手段,尝试二维结构线弹性静力分析的自动计算,包括求解域的自动细分、多项式级数的自动升阶等过程。通过重力坝和带圆孔平板的2个算例验证了方法的可行性,其中第2个算例演示了从CAD的几何信息和计算参数输入到基于精确几何的CAE自动建模、自适应分析、成果自动输出的全过程,初步实现了CAE自动计算以及CAD与CAE的融合。
(SU Hai-dong, FU Zhi, XIE Zhi-qiang. Automatic Two-dimensional Computation Based on Arbitrary Mesh Division[J]. Journal of Yangtze River Scientific Research Institute, 2020, 37(7): 160-166. (in Chinese))
Finite Element Method (FEM) is inconvenient in mesh division and subdivision, difficult in precise modeling of exact geometry, and costs large amount of labor operations. In view of this, we propose an approach of arbitrary mesh subdivision in the 2D solving domain based on convex decomposition idea using Manifold Method based on independent covers presented previously, in which cover meshes are of arbitrary shape, arbitrary connection and arbitrary subdivision. On this basis, with the help of error estimation and h-p version self-adaptive technology in previous studies, we attempt to implement the automatic static analysis of 2D linear-elastic structure, including the automatic subdivision of the solving domain, and the automatic elevation of polynomial orders. Two numerical examples, one of which is a gravity dam and the other is a plate with a small circular hole, are given to illustrate the validity of the present method. Especially in the latter, the whole procedure is exhibited, involving the input of geometry information and computational parameters in CAD, automatic CAE modeling with exact geometry, automatic self-adaptive analysis, as well as the automatic output of computational results. Hence, the automatic CAE computation and CAD/CAE integration are realized preliminarily.
[23]
刘亚军. 基于独立覆盖流形法的流体计算研究[D]. 武汉: raybet体育在线 , 2020.
(LIU Ya-jun. Research on Fluid Calculation Using Numerical Manifold Method Based on Independent Covers[D]. Wuhan: Changjiang River Scientiffic Research Institute, 2020. (in Chinese))
[24]
CHUNG T J. 流体动力学的有限元分析[M]. 张二骏, 龚崇准, 王木兰, 等, 译. 北京: 水利水电出版社, 1980: 203-213.
(CHUNG T J. Finite Element Analysis in Fluid Dynamics[M]. Translated by ZHANG Er-jun, GONG Chong-zhun, WANG Mu-lan, et al. Beijing: Water Resources and Electric Power Press, 1980: 203-213. (in Chinese))

编辑: 刘运飞
PDF(2408 KB)

Accesses

Citation

Detail

段落导航
相关文章

/

Baidu
map