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

作　　者：彭翀 袁会娜 张丙印 张琰 

机构地区：[1]State Key Laboratory of Hydroscience and Engineering (Department of Hydraulic Engineering,Tsinghua University

出　　处：《Journal of Central South University》2013年第11期3229-3246,共18页中南大学学报（英文版）

基　　金：Project(2010CB732103)supported by the National Basic Research Program of China;Project(51179092)supported by the National Natural Science Foundation of China;Project(2012-KY-02)supported by the State Key Laboratory of Hydroscience and Engineering,China

摘　　要：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.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.

关 键 词：meshfree method radial point interpolation method shape parameter influence domain size nodal distribution 

分 类 号：O34[理学—固体力学]