基于蜂窝数值流形元的静弹性力学问题求解.docVIP

基于蜂窝数值流形元的静弹性力学问题求解 [关键词]数值流形方法；蜂窝单元；静弹性力学问题 [摘 要]数值流形方法是一种求解固体力学问题的新颖高效数值方法，而n边形单元（n4）则具有网格划分灵活，能高效进行材料分析等优点。本文利用数值流形方法中数学覆盖系统独立于物理域的优势，基于正六边形数学单元（蜂窝单元）求解二维静弹性力学问题。对典型问题的分析结果表明：在同等情况下，采用正六边形单元的求解精度与正四边形单元的相当，但明显高于正三角形单元上的精度。 Investigation of Elastostatic Problems by the Numerical Manifold Method on Honeycomb Elements Keywords: numerical manifold method; honeycomb elements; elastostatic problems Abstract: The numerical manifold method (NMM) is widely applied to solve both continuous and discontinuous problems in solid mechanics with high efficiency. At the same time, the n-sided elements (n4) are also very attractive due to their greater flexibility in meshing and higher efficiency in materials modeling, compared with the frequently used triangular and quadrilateral elements. To take advantage of the special feature that the mathematical cover system is independent of the physical domain in the NMM, in the present paper, the NMM is applied to solve 2-D elastostatic problems on regular hexagonal (honeycomb) elements. Our results show that the accuracy on hexagonal elements is very close to that on regular quadrilateral elements but much higher than that on regular triangular elements. 数值流形方法是由美籍华人石根华博士在上世纪90年代提出的一种新颖高效的数值方法[1, 2]，该方法能统一求解连续和不连续固体力学问题。自提出至今，该方法在理论和应用层面均取得了较大的发展。理论研究的代表性工作包括：Chen等[3]通过使用高阶多项式覆盖函数发展了高阶数值流形方法；Cheng等[4]在数值流形方法中引入Wilson 非协调元并推导了相关公式；Li等[5]基于单位分解的思想和有限覆盖近似理论提出了无网格数值流形方法；Cheng和Zhang[6]基于四面体和六面体两种基本单元将二维数值流形方法的基本方程推广至三维情形利用数值流形方法模拟了材料中过程数值流形方法以数值流形为核心，在不连续变形分析的块体系统不连续运动学理论基础上，融入有限单元法的连续变形分析能力，而形成一种新的统一计算方法。基本思想是基于有限覆盖技术，将占局部区域的、又相互重叠的有限覆盖组成一个数学覆盖系统，从而去覆盖整个求解域，在各个物理覆盖（数学覆盖和物理域的交集）上构造覆盖函数组成总体函数去逼近求解域的场函数，这也是和有限单元法思想整体区域离散为多个局部区域再集合上的根本不同点。 数值流形方法包括数学覆盖，物理覆盖及流形单元三个基本。 数学覆盖是由用户定义形状任意几何分片。一系列数学覆盖组成数学覆盖系统，该系统不需要与物理域的外边界一致，也不依赖于物理域中诸如夹杂、孔洞和裂纹等不连续的内部结构，只需足够大以覆盖整个物理域。在系统的每个数学覆盖上均定义一个权函数，其满足 (1) (2) 其中：表示自身连续的函数空间。对所有包含点的数学覆盖，定义在它们上的权函数应满足单位分解条件，即