Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics

Abstract

Recently, the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems, however, the accuracy of this method depends on many factors and their influences are not fully investigated yet. In this work, three main factors, i.e., the shape parameters, the influence domain size, and the nodal distribution, on the accuracy of the radial point interpolation method (RPIM) are systematically studied and conclusive results are obtained. First, the effect of shape parameters (R, q) of the multi-quadric basis function on the accuracy of RPIM is examined via global search. A new interpolation error index, closely related to the accuracy of RPIM, is proposed. The distribution of various error indexes on the R-q plane shows that shape parameters q ∈ [1.2, 1.8] and R ∈ [0, 1.5] can give good results for general 3-D analysis. This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics. Second, through numerical experiments, an average of 30–40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics. Third, it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation. Nodal distributions with better uniformity give better results. Furthermore, how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.

Keywords

meshfree method / radial point interpolation method / shape parameter / influence domain size / nodal distribution

References

[1] BelytschkoT, LuY Y, GuL. Element-free Galerkin methods [J]. International Journal for Numerical Methods in Engineering, 1994, 37(2): 229-256

[2] LiuW K, JunS, ZhangY F. Reproducing kernel particle methods [J]. International Journal for Numerical Methods in Fluids, 1995, 20(8/9): 1081-1106

[3] AtluriS N, ZhuT. A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics [J]. Computational Mechanics, 1998, 22(2): 117-127

[4] LiuG R, GuY T. A point interpolation method for two-dimensional solids [J]. International Journal for Numerical Methods in Engineering, 2001, 50(4): 937-951

[5] ChenWGribelM, SchweitzerM A. New RBF collocation methods and kernel RBF with applications [M]. Meshfree Methods for Partial Differential Equations. Vol. 1, 2000

[6] ChenW, FuZ J. Boundary particle method for inverse Cauchy problems of inhomogeneous Helmholtz equations [J]. Journal of Marine Science and Technology, 2009, 17(3): 157-163

[7] XieH, NogamiT, WangJ G. A radial boundary node method for two-dimensional elastic analysis [J]. Engineering Analysis with Boundary Elements, 2003, 27: 853-862

[8] WangJ G, LiuG R. A point interpolation meshless method based on radial basis functions [J]. International Journal for Numerical Methods in Engineering, 2002, 54(11): 1623-1648

[9] LiuG R, GuY T. A matrix triangularization algorithm for the polynomial point interpolation method [J]. Computer Methods in Applied Mechanics and Engineering, 2003, 192(19): 2269-2295

[10] LiuG R, ZhangG Y, GuY T, WangY Y. A meshfree radial point interpolation method (RPIM) for three-dimensional solids [J]. Computational Mechanics, 2005, 36(6): 421-430

[11] WangJ G, LiuG R, WuY G. A point interpolation method for simulating dissipation process of consolidation [J]. Computer Methods in Applied Mechanics and Engineering, 2001, 190(45): 907-5922

[12] WangJ G, LiuG R, LinP. Numerical analysis of Biot’s consolidation process by radial point interpolation method [J]. International Journal of Solids and Structures, 2002, 39(6): 1557-1573

[13] WangJ G, ZhangB, NogamiT. Wave-induced seabed response analysis by radial point interpolation meshless method [J]. Ocean Engineering, 2004, 31(1): 21-42

[14] WangJ G, LiuG R. On the optimal shape parameters of radial basis functions used for 2-D meshless methods [J]. Computer Methods in Applied Mechanics and Engineering, 2002, 191(23/24): 2611-2630

[15] LiuG R, GuY TAn introduction to meshfree methods and their programming [M], 2005DordrechtSpringer

[16] BuhmannM DRadial basis functions: Theory and implementations [M], 2003CambridgeCambridge University Press57-73

[17] WendlandHScattered data approximation[M], 2005CambridgeCambridge University Press101-125

[18] HardyR L. Multiquadric equations of topography and other irregular surfaces [J]. J Geophys Res, 1971, 76(8): 1905-1915

[19] HardyR L. Theory and applications of multiquadric-biharmonic method: 20 Years of discovery 1968–1988 [J]. Computers & Mathematics with Applications, 1990, 19(88): 163-208

[20] FrankeR. Scatter data interpolation: Test of some methods [J]. Mathematics of Computation, 1982, 38(157): 181-200

[21] RippaS. An algorithm for selecting a good value for the parameter c in radial basis function interpolation [J]. Advances in Computational Mathematics, 1999, 11(2): 193-210

[22] FornbergB, ZuevJ. The Runge phenomenon and spatially variable shape parameters in RBF interpolation [J]. Computers & Mathematics with Applications, 2007, 54(3): 379-398

[23] SarraS A, SturgillD. A random variable shape parameter strategy for radial basis function approximation methods [J]. Engineering Analysis with Boundary Elements, 2009, 33(11): 1239-1245

[24] XiaoJ R, McCarthyM A. A local Heaviside weighted meshless method for two-dimensional solids using radial basis functions [J]. Computational Mechanics, 2003, 31(3): 301-315

[25] LiL, ZhuJ, ZhangS. A hybrid radial boundary node method based on radial basis point interpolation [J]. Engineering Analysis with Boundary Elements, 2009, 33(11): 1273-1283

[26] XiaoJ R, GamaB A, GillespieJ W, KansaE. Meshless solutions of 2D contact problems by subdomain variational inequality and MLPG method with radial basis functions [J]. Engineering Analysis with Boundary Elements, 2005, 29(2): 95-106

[27] DinisL M J S, JorgeR M N, BelinhaJ. Analysis of 3D solids using the natural neighbour radial point interpolation method [J]. Computer Methods in Applied Mechanics and Engineering, 2007, 196(13/16): 2009-2028

[28] QianW-changElastic mechanics [M], 1956BeijingScience Press

[29] TimoshenkoS P, GoodierJ NTheory of elasticity [M], 1970New YorkMcGraw-Hill

[30] LiuG R, GuY T. Boundary meshfree methods based on the boundary point interpolation methods [J]. Engineering Analysis with Boundary Elements, 2004, 28(5): 475-487

[31] DuQ, FaberV, GunzburgerM. Centroidal Voronoi tessellations: Applications and algorithms [J]. SIAM Review, 1999, 41(4): 637-676

[32] LiuY, JieY-xin. Simulated-annealing-based algorithm for generating point sets for meshfree methods [J]. Journal of Tsinghua University: Science and Technology, 2008, 48(6): 959-962