The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ...The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.展开更多
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 the...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.展开更多
In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in th...In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in the central area,and the radial point interpolation method(RPIM)based on the global weakness is utilized to construct the meshfree shape function.The Governing equations in each background gird are solved by Gaussian integration.In the extended area where the points are sparsely distributed,to avoid the instability of the meshfree method,finite element method(FEM)with regular grids is used to solve the governing equation.Finally,the meshfree and finite element governing equations are coupled by the continuity of the field at the interfaces,and the direct solution technique is used to realize the 3D MT forward modeling.Numerical experiments of several typical electrical models are used to verify the effectiveness of the method.展开更多
A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using ...A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.展开更多
A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(F...A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.展开更多
In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis...In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis functions for interpolation.All these radial basis functions include shape parameters.The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory.The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM.The RPIM is studied based on the global Galerkin weak form performed using two integration technics:classical Gaussian integration and the strain smoothing integration scheme.The numerical performance of this method is tested on their behavior on curve fitting,and on three elastic mechanical problems with regular or irregular nodes distributions.A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system.All resulting RPIM methods perform very well in term of numerical computation.The Smoothed Radial Point Interpolation Method(SRPIM)shows a higher accuracy,especially in a situation of distorted node scheme.展开更多
基金supported by the Key Program of the National Natural Science Foundation of China (Grand No. 51138001)the China-German Cooperation Project (Grand No. GZ566)+1 种基金the Innovative Research Groups Funded by the National Natural Science Foundation of China (Grand No. 51121005)the Special Funds for the Basic Scientific Research Expenses for the Central University (Grant No. DUT13LK16)
文摘The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.
基金Project(2010CB732103)supported by the National Basic Research Program of ChinaProject(51179092)supported by the National Natural Science Foundation of ChinaProject(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.
基金Supported by Project of National Natural Science Foundation of China(No.42074120).
文摘In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in the central area,and the radial point interpolation method(RPIM)based on the global weakness is utilized to construct the meshfree shape function.The Governing equations in each background gird are solved by Gaussian integration.In the extended area where the points are sparsely distributed,to avoid the instability of the meshfree method,finite element method(FEM)with regular grids is used to solve the governing equation.Finally,the meshfree and finite element governing equations are coupled by the continuity of the field at the interfaces,and the direct solution technique is used to realize the 3D MT forward modeling.Numerical experiments of several typical electrical models are used to verify the effectiveness of the method.
基金supported by the Natural Science Foundation of China (10972075,10802029)the Hunan Provincial Natural Science Foundation of China(10JJ3036)Scientific Research Fund of Hunan Provincial Education Depart ment of China(08C230)~~
基金supported by the National 973 Scientific and Technological Innovation Project (No. 2004CB719402)National Natural Science Foundation of China (No. 10672055)+3 种基金Key Project of NSFC (No. 60635020)Natural Science Foundation for Out standing Youth of China (No. 50625519)Hunan Provincial Natural Science Foundation of China (No. 07JJ6002)Scientific Research Fund of Hunan Provincial Education Department of China (No. 08C230)
文摘A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.
基金co-supported by the National Key R&D Program of China(Nos.2018YFF01012401-05)the National Natural Science Foundation of China(No.51975243)+2 种基金Jilin Provincial Department of Education(No.JJKH20180084KJ),Chinathe Fundamental Research Funds for the Central Universities and Jilin Provincial Department of Science&Technology Fund Project,China(Nos.20170101043JC and 20180520072JH)Graduate Innovation Fund of Jilin University,China(No.101832018C184).
文摘A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.
文摘In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis functions for interpolation.All these radial basis functions include shape parameters.The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory.The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM.The RPIM is studied based on the global Galerkin weak form performed using two integration technics:classical Gaussian integration and the strain smoothing integration scheme.The numerical performance of this method is tested on their behavior on curve fitting,and on three elastic mechanical problems with regular or irregular nodes distributions.A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system.All resulting RPIM methods perform very well in term of numerical computation.The Smoothed Radial Point Interpolation Method(SRPIM)shows a higher accuracy,especially in a situation of distorted node scheme.
