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.展开更多
A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal p...A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal points;one class is that the source term is time dependent, and the other class is that the source term is time and space dependent. Some numerical experiments are presented and 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.展开更多
Laminated composite materials are widely implemented in several engineering constructions. For its relative light weight, these materials are suitable for aerospace, military, marine, and automotive structural applica...Laminated composite materials are widely implemented in several engineering constructions. For its relative light weight, these materials are suitable for aerospace, military, marine, and automotive structural applications. To obtain safe and economical structures, the modelling analysis accuracy is highly relevant. Since meshless methods in the recent years achieved a remarkable progress in computational mechanics, the present work uses one of the most flexible and stable interpolation meshless technique available in the literature—the Radial Point Interpolation Method(RPIM).Here, a 2 D approach is considered to numerically analyse composite laminated beams. Both the meshless formulation and the equilibrium equations ruling the studied physical phenomenon are presented with detail. Several benchmark beam examples are studied and the results are compared with exact solutions available in the literature and the results obtained from a commercial finite element software. The results show the efficiency and accuracy of the proposed numeric technique.展开更多
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.展开更多
Combining the radial point interpolation method (RPIM), the dual reciprocity method (DRM) and the hybrid boundary node method (HBNM), a dual reciprocity hybrid radial boundary node method (DHRBNM) is proposed ...Combining the radial point interpolation method (RPIM), the dual reciprocity method (DRM) and the hybrid boundary node method (HBNM), a dual reciprocity hybrid radial boundary node method (DHRBNM) is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square (MLS) in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it's computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.展开更多
基金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.
文摘A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal points;one class is that the source term is time dependent, and the other class is that the source term is time and space dependent. Some numerical experiments are presented and 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 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.
基金the funding provided by Ministério da Ciência,Tecnologia e Ensino Superior--Fundaca para a Ciência e a Tecnologia(Portugal)(Grants SFRH/BPD/75072/2010,SFRH/BPD/111020/2015)the inter-institutional projects from LAETA(Grant UID/EMS/50022/2013)+1 种基金the funding of Project NORTE-010145-FEDER-000022-SciTech-Science and Technology for Competitive and Sustainable Industriescofinanced by Programa Operacional Regional do Norte(Grant NORTE2020),through Fundo Europeu de Desenvolvimento Regional(FEDER)
文摘Laminated composite materials are widely implemented in several engineering constructions. For its relative light weight, these materials are suitable for aerospace, military, marine, and automotive structural applications. To obtain safe and economical structures, the modelling analysis accuracy is highly relevant. Since meshless methods in the recent years achieved a remarkable progress in computational mechanics, the present work uses one of the most flexible and stable interpolation meshless technique available in the literature—the Radial Point Interpolation Method(RPIM).Here, a 2 D approach is considered to numerically analyse composite laminated beams. Both the meshless formulation and the equilibrium equations ruling the studied physical phenomenon are presented with detail. Several benchmark beam examples are studied and the results are compared with exact solutions available in the literature and the results obtained from a commercial finite element software. The results show the efficiency and accuracy of the proposed numeric technique.
基金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.
基金Project supported by the National Basic Research Program of China (No. 2010CB732006)the CAS/SAFEA International Partnership Program for Creative Research Teams (No. KZCX2-YW-T12)the National Natural Science Foundation of China (No. 11002154)
文摘Combining the radial point interpolation method (RPIM), the dual reciprocity method (DRM) and the hybrid boundary node method (HBNM), a dual reciprocity hybrid radial boundary node method (DHRBNM) is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square (MLS) in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it's computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.