A new meshless method based on MLPG for elastic dynamic problems

The basic concept and numerical implementation of a new local Petrov-Galerkin method for solving a dynamic problem are presented in this paper. It uses a radial basis function (RBF) coupled with a polynomial basis function as a trial function, and uses the Heaviside function as a test function of the weighted residual method. The shape function has the Kronecker Delta properties for the trial-function-interpolation, and so no additional treatment to impose essential boundary conditions. The method does not involve any domain and singular integrals to generate the global effective stiffness matrix except for the mass and damping matrice; it only involves a regular boundary integral. It possesses a great flexibility in dealing with the numerical model of the elastic dynamic problem under various boundary conditions with arbitrary shapes. The Newmark family of methods is adopted in computation. The numerical results also show that using a multiquadrics (MQ) function with the polynomial basis function as the interpolation function can give quite accurate numerical results. The aQ and aS are investigated which are parameters of the radii of the sub-domain and influence domain, respectively.