20-Node Rational Elements for 3D Anisotropic Elastic Problems
-
摘要: 常规有限元方法的插值函数通常仅仅从数学层面上考虑单元的几何性质,忽视了与物理问题相关的物性参数,因此可能降低数值分析的效果.理性有限元的构造方法与常规有限元法不同,其插值函数使用的是控制微分方程解析解的线性组合,求解过程是在物理域内直接列式,对单元的应变场和应力场同时进行插值,并在单元级别考虑分片实验的要求并直接进行修正,最终形成与问题物性参数紧密相关的单元刚度阵.该方法避免了传统方法对物理问题和数学问题的割裂,可显著提高数值分析的稳定性和精度.利用空间各向异性问题的基本解,从最小势能原理出发,构造出两种满足分片实验要求的二十节点理性块体单元.数值算例表明,所给出的理性单元不仅具有较高的求解精度,而且具有良好的数值稳定性.Abstract: For the conventional finite element method, only the geometry and node locations of elements were considered in the interpolation functions, while the physical parameters which reflect the key features of the physical problems were ignored, so its numerical performance may be not satisfying in some cases. The construction of the rational finite element method was different from that of the conventional finite element method. The linear combinations of the fundamental solutions to the problem’s controlling differential equations were used as the interpolation functions, so the stress and strain fields were interpolated directly in the physical domain at the same time. The transfer matrix was modified at the element level to pass the patch test, and the resulting element stiffness matrix was related closely to the physical parameters of the problem. The rational finite element avoids the separation between the mathematical and physical aspects of a problem, so the stability and accuracy of numerical analysis could be improved significantly. Two kinds of 20-node rational brick elements based on the principle of minimum potential energy and satisfying the requirements of the patch test, were constructed according to the fundamental solutions to general 3D anisotropic problems. Numerical examples show that the rational elements give numerical results with not only high accuracy, but also good numerical stability.
-
Key words:
- general anisotropy /
- rational finite element /
- 20-node brick element
-
[1] Zienkiewicz O C, Taylor R L.The Finite Element Method for Solid and Structural Mechanics [M]. Burlington: Butterworth-Heinemann, 2005: 24-28. [2] Ergatoudis I, Irons B M, Zienkiewicz O C. Curved, isoparametric,“quadrilateral” elements for finite element analysis[J].International Journal of Solids and Structures,1968,4(1): 31-42. [3] Pian T H H, Chen D. On the suppression of zero energy deformation modes[J].International Journal for Numerical Methods in Engineering,1983,19(12): 1741-1752. [4] Punch E F, Atluri S N. Development and testing of stable, invariant, isoparametric curvilinear 2-and 3-D hybrid-stress elements[J].Computer Methods in Applied Mechanics and Engineering,1984,47(3): 331-356. [5] Chen W, Cheung Y K. Three-dimensional 8-node and 20-node refined hybrid isoparametric elements[J].International Journal for Numerical Methods in Engineering,1992,35(9): 1871-1889. [6] 钟万勰, 纪峥. 理性有限元[J]. 计算结构力学及其应用, 1996,13(1): 1-8.(ZHONG Wan-xie, JI Zheng. Rational finite element[J].Computational Structural Mechanics and Applications, 1996,13(1): 1-8.(in Chinese)) [7] 钟万勰, 纪峥. 平面理性元的收敛性证明[J]. 力学学报, 1997,29(6): 676-685.(ZHONG Wan-xie, JI Zheng. The convergence proof of the plane rational finite element[J].Acta Mechanica Sinica,1997,29(6): 676-685.(in Chinese)) [8] 纪峥, 钟万勰. 平面理性四节点及五节点四边形有限元[J]. 计算力学学报, 1997,14(1): 19-27.(JI Zheng, ZHONG Wan-xie. Rational plane quadrilateral four and five nodes finite elements[J].Chinese Journal of Computational Mechanics,1997,14(1): 19-27.(in Chinese)) [9] 王永富, 钟万勰. 空间理性八节点块体元[J]. 应用力学学报, 2003,20(3): 131-135.(WANG Yong-fu, ZHONG Wan-xie. A rational hexahedron 8-node finite element[J].Chinese Journal of Applied Mechanics,2003,20(3): 131-135.(in Chinese)) [10] 王莹, 孙秀珍, 刘迎曦, 于驰, 张军, 于申. OSAHS患者与正常人上呼吸道流场特性比较[J]. 大连理工大学学报, 2009,49(4): 476-481.(WANG Ying, SUN Xiu-zhen, LIU Ying-xi, YU Chi, ZHANG Jun, YU Shen. Comparisons of flow characteristics of upper airway between patients with OSAHS and normal adults[J].Journal of Dalian University of Technology,2009,49(4): 476-481.(in Chinese)) [11] 刘钧玉, 林皋, 杜建国. 基于 SBFEM 的多裂纹问题断裂分析[J]. 大连理工大学学报, 2008,48(3): 392-397.(LIU Jun-yu, LIN Gao, DU Jian-guo . Fracture analysis of multiple crack problems based on SBFEM[J].Journal of Dalian University of Technology ,2008,48(3): 392-397.(in Chinese)) [12] Ozaki S, Hikida K, Hashiguchi K. Elastoplastic formulation for friction with orthotropic anisotropy and rotational hardening[J].International Journal of Solids and Structures,2012,49(3): 648-657. [13] 王华丽. 理性有限元插值函数的构造[J]. 教育教学论坛, 2010(23): 180-182, 185.(WANG Hua-li. The interpolation function construction of rational finite element[J].Education Teaching Forum,2010(23): 180-182, 185.(in Chinese)) [14] Macneal R H, Harder R L. A proposed standard set of problems to test finite element accuracy[J].Finite Elements in Analysis and Design,1985,1(1): 3-20.
点击查看大图
计量
- 文章访问数: 1140
- HTML全文浏览量: 93
- PDF下载量: 984
- 被引次数: 0