二维弹性力学问题的光滑无网格伽辽金法

马文涛

doi:10.6052/0459-1879-18-135

二维弹性力学问题的光滑无网格伽辽金法

doi: 10.6052/0459-1879-18-135

马文涛^2), ,

宁夏大学数学统计学院, 银川 750021

基金项目: 1)国家自然科学基金项目(51769026)和宁夏自然科学基金重点项目(NZ17002)资助.

详细信息

作者简介:
2)马文涛, 教授, 主要研究方向: 计算力学及其工程应用.E-mail:wt-ma2002@163.com

通讯作者:
马文涛

中图分类号: O241,O343;
计量
- 文章访问数: 905
- HTML全文浏览量: 96
- PDF下载量: 242
- 被引次数: 0
出版历程
- 收稿日期: 2018-04-25
- 刊出日期: 2018-09-18

A SMOOTHED MESHFREE GALERKIN METHOD FOR 2D ELASTICITY PROBLEM

Ma Wentao^{2)
, ,}

School of Mathematics and Statistics, Ningxia University, Yinchuan 750021, China

摘要

摘要: 计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM) 的深入发展. 为了提高EFGM 的计算速度, 本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法. 该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数; 基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变, 根据广义光滑伽辽金弱形式建立系统离散方程. 两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成. 为了提高方法的精度, 由Richardson外推法确定两层光滑网格上的最优光滑应变. 几个数值算例验证了该方法的精度和计算效率. 数值结果表明, 随着光滑积分网格数目的增加, 光滑无网格伽辽金法的计算精度逐步接近EFGM 的, 但计算效率要远远高于EFGM的. 另外, 光滑无网格伽辽金法的边界条件可以像有限元那样直接施加. 从计算精度和效率综合考虑, 光滑无网格伽辽金法比EFGM具有更好的数值表现, 具有十分广阔的发展空间.
- 无网格法 /
- 光滑应变 /
- 两层嵌套光滑三角形网格 /
- 计算效率
Abstract: Despite clear general progress with element-free Galerkin method (EFGM), its low computational efficiency becomes a technical issue in the simulation of realistic problems. To improve the efficiency of EFGM, a smoothed meshfree Galerkin method is presented for the 2D elasticity problem. In the method presented, displacement fields are constructed using the moving least square (MLS) approximation and strains are smoothed over two-level nesting smoothed triangular cells based on the generalized gradient smoothing technique. Then, the generalized smoothed Galerkin (GS-Galerkin) weak form is used to create the discretized system equations. Each two-level nesting smoothed triangular cells include the triangular background cell itself and four equal-area triangular sub-cells, respectively. According to the Richardson extrapolation method, an optimal combination of the two-level smoothed strains can be obtained. Since the present method uses the linear interpolation on the boundary of problem domain, the boundary conditions including the essential and natural boundary conditions can directly impose as that in FEM. Several examples, including the cantilever beam, infinite plate with a circle hole, infinite plate with a central crack and the twin-arched tunnel, are investigated to demonstrate the accuracy and efficiency of the present method. The numerical results show that with more smoothing sub-cells by using in the smoothed meshfree Galerkin method, higher numerical accuracy can be obtained. In addition, the present method is higher efficient than EFGM. As a consequence, the smoothed meshfree Galerkin method with two-level nesting smoothed triangular cells significantly outperforms the EFGM and is very successful and attractive numerical method for solving the elasticity problems.
- meshfree method /
- strain smoothing /
- two-level nesting smoothed triangular cells /
- computational efficiency

HTML全文

参考文献(1)

[1]	Belytschko T, Lu YY, Gu L.Element-free Galerkin methods. International Journal for Numerical Methods in Engineering, 1994, 37: 229-256
[2]	Belytschko T, Kronggauz Y, Organ D, et al.Meshless methods: An overview and recent development. Computational Methods in Applied Mechanics and Engineering, 1996, 139: 3-47
[3]	Liu WK, Jun S, Zhang YF.Reproducing kernel particle methods. International Journal for Numerical Methods in Fluids, 1995, 20: 1081-1106
[4]	Atluri SN, Shen SP.The Meshless Local Petrov-Galerkin (MLPG) Method. Tech Science, 2002
[5]	Nguyen VP, Rabczuk S, Bordas S, et al.Meshless Methods: A review and computer implementation aspects. Mathematics and Computers in Simulation, 2008, 79: 763-813
[6]	Rabczuk T, Bordas S, Zi G.A three-dimensional meshfree method for continuous multiplecrack initiation, nucleation and propagation in statics and dynamics. Computational Mechanics, 2007, 40: 473-495
[7]	马文涛, 许艳, 马海龙. 修正的内部基扩充无网格法求解多裂纹应力强度因子. 工程力学, 2015, 32(10): 18-24
[7]	(Ma Wentao, Xu Yan, Ma Hailong.Solving stress intensity factors of multiple cracks by using a modified intrinsic basis enriched meshless method. Engineering Mechanics, 2015, 32(10): 18-24 (in Chinese))
[8]	Atluri SN, Zhu T.A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics. Computational Mechanics, 1998, 22: 117-127
[9]	Wang JG, Liu GR.A point interpolation meshless method based on radial basis function.International Journal for Numerical Methods in Engineering, 2002, 54: 1623-1648
[10]	杨建军, 郑健龙. 无网格局部强弱法求解不规则域问题. 力学学报, 2017, 49(3): 659-666
[10]	(Yang Jianjun, Zheng Jianlong.Meshless local strong-weak (MLSW) method for irregular domain problems. Journal of Theoretical and Applied Mechanics, 2017, 49(3): 659-666 (in Chinese))
[11]	Duan QL, Li XK, Zhang HW, et al.Second-order accurate derivatives and integration schemes for meshfree methods. International Journal for Numerical Methods in Engineering, 2012, 92: 399-424
[12]	Beissl S, Belytschko T, Nodal integration of the element-free Galerkin method. Computer Method in Applied Mechanics and Engineering, 1996, 139(19): 49-64
[13]	Chen JS, Wu CT, Yoon S, et al.A stabilized conforming nodal integration for Galerkin meshfree methods. International Journal for Numerical Methods in Engineering, 2001, 50: 435-466
[14]	Chen JS, Yoon J, Wu CT.Nonlinear version of stabilized conforming nodal integration for Galerkin meshfree methods. International Journal for Numerical Methods in Engineering, 2002, 53: 2587-2615
[15]	Liu GR.A generalized gradient smoothing technique and smoothed bilinear form for Galerkin formulation of a wide class of computational methods. International Journal of Computational Methods, 2008, 5(2): 199-236
[16]	Liu GR, Zhang GY.A normed G space and weaked weak (W) formulation of a cell-based smoothed point interpolation method. International Journal of Computational Methods, 2009, 6(1): 147-179
[17]	Liu GR, Nguyen TT, Smoothed Finite Element Method. Boca Raton: CRC Press, 2010
[18]	Liu GR, Zhang GY.Smoothed Point Interpolation Methods: G Space Theory and Weakened Weak Forms. New Jersey: World Scientific , 2013
[19]	Bordas PAS, Rabczuk T, Hung NX, et al.Strain smoothing in FEM and XFEM. Computers and Structures, 2010, 88: 1419-1443
[20]	Chen L, Rabczuk T, Bordas PAS, et al. Extended finite element method with edge-based strain smoothing (ESm-XFEM) for linear elastic crack growth. Computational Methods in Applied Mechanics and Engineering, 2012, 209-212: 250-265
[21]	Liu GR, Smoothed finite element methods (S-FEM): An overview and recent developments. Arch Computat. Methods Eng, 2018, 25: 397-435
[22]	Cui XY, Li GY.Smoothed Galerkin methods using cell-wise strain smoothing technique. Engineering Analysis with Boundary Elements, 2012, 36: 825-835
[23]	杜超凡, 章定国. 光滑节点插值法: 计算固有频率下界值的新方法. 力学学报, 2015, 47(5): 839-847
[23]	(Du Chaofan, Zhang Dingguo.Node-based smoothed point interpolation method: A new method for computing lower bound of natural frequency. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(5): 839-847
[24]	Liu F, Yu CY, Yang YT.An edge-based smoothed numerical manifold method and its application to static, free and forced vibration analyses. Engineering Analysis with Boundary Elements, 2018, 86: 19-30
[25]	Duan QL, Wang BB, Gao X, et al.Quadratically consistent nodal integration for second order meshfree Galerkin methods. Computational Mechanics, 2014, 54: 353-368
[26]	Wang BB, Duan QL, Shao YL, et al.An efficient nodal integration with quadratic exactness for three-dimensional meshfree Galerkin methods. Engineering Analysis with Boundary Elements, 2016, 70: 99-113
[27]	邵玉龙, 段庆林, 李锡夔等. 功能梯度材料的二阶一致无网格. 工程力学, 2017, 34(3): 15-21
[27]	(Shao Yulong, Duan Qinglin, Li Xi kui, et al. Quadratically consistent meshfree method for functionally graded materials. Engineering Mechanics, 2017 34(3): 15-21 (in Chinese))
[28]	Wang DD, Wu JC.An efficient nesting sub-domain smoothing integration algorithm with quadratically exactness for Galerkin meshfree methods. Computer Methods in Applied Mechanics and Engineering, 2016, 298: 485-519
[29]	Brezinski C, Zaglia MR.Extrapolation Methods: Theory and Practice. North-Holland, 1991
[30]	Timoshenko SP, Goodier JN.Theory of Elasticity, 3rd edition. New York: McGraw-Hill, 1987
[31]	Williams ML.On the stress distributions at the base of a stationary crack. Journal of Applied Mechanics, 1957, 24: 109-113
[32]	徐荣桥. 结构分析的有限元法与MATLAB程序设计. 北京: 人民交通出版社, 2006
[32]	(Xu Rongqiao, Finite Element Method and Matlab Programing Design for Structural Analysis. Beijing: China Communications Press, 2006 (in Chinese))