李佳龙, 李钢, 余丁浩. 多边形比例边界有限元非线性高效分析方法[J]. 工程力学, 2020, 37(9): 8-17. DOI: 10.6052/j.issn.1000-4750.2019.10.0634
引用本文: 李佳龙, 李钢, 余丁浩. 多边形比例边界有限元非线性高效分析方法[J]. 工程力学, 2020, 37(9): 8-17. DOI: 10.6052/j.issn.1000-4750.2019.10.0634
LI Jia-long, LI Gang, YU Ding-hao. POLYGON SCALED BOUNDARY FINITE ELEMENT NONLINEAR EFFICIENT ANALYSIS METHOD[J]. Engineering Mechanics, 2020, 37(9): 8-17. DOI: 10.6052/j.issn.1000-4750.2019.10.0634
Citation: LI Jia-long, LI Gang, YU Ding-hao. POLYGON SCALED BOUNDARY FINITE ELEMENT NONLINEAR EFFICIENT ANALYSIS METHOD[J]. Engineering Mechanics, 2020, 37(9): 8-17. DOI: 10.6052/j.issn.1000-4750.2019.10.0634

多边形比例边界有限元非线性高效分析方法

POLYGON SCALED BOUNDARY FINITE ELEMENT NONLINEAR EFFICIENT ANALYSIS METHOD

  • 摘要: 比例边界有限元法作为一种高精度的半解析数值求解方法,特别适合于求解无限域与应力奇异性等问题,多边形比例边界单元在模拟裂纹扩展过程、处理局部网格重剖分等方面相较于有限单元法具有明显优势。目前,比例边界有限元法更多关注的是线弹性问题的求解,而非线性比例边界单元的研究则处于起步阶段。该文将高效的隔离非线性有限元法用于比例边界单元的非线性分析,提出了一种高效的隔离非线性比例边界有限元法。该方法认为每个边界线单元覆盖的区域为相互独立的扇形子单元,其形函数以及应变-位移矩阵可通过半解析的弹性解获得;每个扇形区的非线性应变场通过设置非线性应变插值点来表达,引入非线性本构关系即可实现多边形比例边界单元高效非线性分析。多边形比例边界单元的刚度通过集成每个扇形子单元的刚度获取,扇形子单元的刚度可采用高斯积分方案进行求解,其精度保持不变。由于引入了较多的非线性应变插值点,舒尔补矩阵维数较大,该文采用Woodbury近似法对隔离非线性比例边界单元的控制方程进行求解。该方法对大规模非线性问题的计算具有较高的计算效率,数值算例验证了算法的正确性以及高效性,将该方法进行推广,对实际工程分析具有重要意义。

     

    Abstract: The scaled boundary finite element method (SBFEM) is a high-precision semi-analytical numerical solution method, which is especially suitable for solving problems such as unbounded media and stress singularity, and the advantage of the polygon boundary element is more obvious to simulate the crack growth process and local mesh re-segmentation problems than the finite element method. At present, the scaled boundary finite element method is more concerned with the solution of the linear elasticity problem, while the research of the nonlinear scaled boundary element is in its infancy. An efficient inelasticity-separated scaled boundary finite element method (IS-SBFEM) is proposed based on the basic theory of the SBFEM and the inelasticity-separated theory. The proposed method considers that the sector sub-element of each boundary line covered domain is independent, and its shape function and strain-displacement matrix can be obtained by the semi-analytical elastic solution. Moreover, the nonlinearity strain field of each sector sub-element is establish by introducing nonlinear strain interpolation points, and the nonlinear constitutive relations can be introduced to achieve efficient nonlinear analysis of polygon scaled boundary element. The stiffness matrix of the polygon scaled boundary element can be obtained by assembling the stiffness of each sector sub-element, the numerical integration of each sector domain can be obtained by using Gaussian integration scheme, and its accuracy remains unchanged. Because more nonlinear strain interpolation points are introduced, the dimension of the Schur complement matrix is larger. The Woodbury approximation approach is used to solve the governing equations of the inelasticity-separated scaled boundary element. This method has efficiency advantages for the calculation of large-scale nonlinear problems. Numerical examples are adopted to verify the correctness and efficiency of the algorithm. Popularize this method that is important to practical engineering analysis.

     

/

返回文章
返回