Generality of physical or mechanical problems are modeled by partial differential equations (PDEs). Moreover, a high number of numerical methods and techniques have been developed for the approximation of solution of the PDEs in much research that focused mainly on the improvement of accuracy and efficiency of them. In recent years, attention has turned on the development of meshless methods, especially for the numerical solution of partial differential equations. The meshless local Petrov-Galerkin (MLPG) approach  -  has become very attractive, as a promising method for solving problems. This method based on a local weak form and Moving Least Squares (MLS) approximation     . The main advantage of this method over the widely used ﬁnite element methods (FEM) is that it does not need any mesh, either for the interpolation of the solution variables or for the integration of the weak forms   . The MLPG method and its variations have been applied for elastodynamic and elastostatic problems by several authors. For example, it is applied to free and forced vibration analysis for solids by Gu Y. T. et al.  they found that the parameter which decides the size of the sub-domain needs to be chosen carefully. Long S. Y. et al., have applied MLPG5 method used the Heaviside function as the test function for elastic dynamic problems  . They found a good agreement compared with the results obtained by (FEM). This method also applied by Ping Xia et al. in elastic dynamic analysis of moderately thick plate using meshless LRPIM  . They have used the Newmark method for solving the dynamic problem and have studied the eﬀects of the size of the quadrature subdomain and the inﬂuence domain on the dynamic properties. They found if appropriate sizes are selected, good results and stability can be obtained. A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis by De-An Hu et al.  , they have used the standard implicit Newmark’s time integration scheme for solving the global dynamic system equations obtained by assembling all local discrete equations. Recent results founded by Moussaoui and al concerning the effects of support domain for elastostatic problem by MLPG  .
The MLPG formulation is proposed in this paper to extend the MLPG method to dynamic analysis and for solving the problem of a thin elastodynamic homogeneous rectangular plate  . The Laplace transform  is applied to eliminate the time variable, then, the obtained equations by the local formulation becomes in function with coefficient of Laplace transform. The Stehfest inversion method is applied to obtain the time-dependent solutions  . The result presented for different values of and with regular distribution of nodes . After, the results are presented with fixed values of and for different time-step. We found large domains of , and time-step by using Laplace transform method. The integral equations have a very simple non-singular form. Moreover, both the contour and domain integrations can be easily carried out in rectangular sub-domains.
This paper is organized as follows: Section 2 introduces the least squares approximation (MLS) for combination of shape function. In Section 3, the Basic equations of elastodynamics and their Laplace transforms are proposed. In Section 4, the MLPG formula including the local weak form in Laplace-transformed domain are developed, using the weighted residual method locally from the dynamic partial differential equation. The numerical results and discussions for 2D problem example are given in Section 5. Finally, the paper ends with the conclusion.
2. Moving Least Square (MLS) Approximation
We Consider a sub-domain the neighbourhood of point , which is located within the problem domain . To approximate a function in , a finite set of monomial basis functions , is considered in the space coordinates in two-dimension is given by:
The approximation function of a field variable is defined in a subdomain by:
where is the monomial basis function of the spatial coordinates for two dimensional problem, and m is the number of the monomial basis functions. is a function of point and time and is a vector of coefficients given by:
The coefficient is obtained at any point by minimizing a weigh- ted discrete L2 norm
where n is the number of nodes in the support domain of for which the weight function , and is the nodal parameter of at .
The stationarity of with respect to gives:
which leads to the following linear relations:
is the vector that collects the nodal displacements for all nodes in the support domain:
is called the weighted moment matrix defined by:
and the matrix B is defined by:
Solving from Equation (6) as:
and substituting the above Equation (10) back into Equation (2) we obtain:
where is the matrix of MLS shape functions corresponding n nodes in the support domain of the point and can be written as:
The shape function for the ith node is defined by:
The following quartic spline function is used, it has the following form of
In which is the distance from node to point , and is the size of the influence domain for the weight function.
3. Basic Equations of Elastodynamics
The governing equations for a linear two-dimensional elastodynamic problem on a domain Ω bounded by a boundary Γ are:
where the mass density, c is the damping coefficient, is the acceleration, the velocity, the stress tensor, the body force tensor, and denotes
The initial and boundary conditions are given as follows:
where and are the prescribed displacement and traction on the boundary and respectively, and denote the initial displacement and initial velocity, is the unit outward normal to the boundary . and are complementary subsets of .
The Laplace-transform  of a function is deﬁned by:
Then the Laplace-transform of the basic Equation (15) gives the general expression of the partial differential equation:
4. The MLPG Weak Formulation in Laplace-Transformed Domain
The MLPG (meshless local Petrov-Galerkin) method constructs the weak form over local subdomain such as , which is a small region taken for each node inside the global domain. The local subdomains overlap, and cover the whole global domain . In the present paper, the local subdomains are taken to be of a quadrature shape. The local weak form   of the governing Equation (25) can be written as:
where is a test function:
and applying the Gauss divergence theorem we can write:
is the boundary of the local subdomain, which is consisted of three parts, i.e.
(See Figure 1)
is the local boundary that is totally inside the global domain,
is the part of the local boundary, which lies on the global boundary with prescribed tractions, i.e.
is the part of the local boundary that lies on the global boundary with prescribed displacements, i.e.
Figure 1. The support domain ΩS and integration domain Ωq for node I.
The local weak form in Equation (29) is leading to the following local integral equation:
The strains can be obtained using the approximated displacements:
Considering Equation (11) with Laplace transform:
The constitutive equation gives the relationship between the stress and the strain:
The traction vectors at a boundary point are approximated by numbers of nodal values as:
where the matrix related to the normal vector of by:
The matrix represented by the gradients of the shape functions as:
The stress-strain matrix D for plane stress is deﬁned by:
In which E is the Young’s modulus and is the Poisson’s ratio.
Obeying the boundary conditions at those nodal points on the global boundary, where displacements are prescribed, and making use of the approximation formulae Equation (11), we obtains the discretized form of the displacement boundary conditions:
Substituting equations. (36) and (37) into Equation (32) we obtain the discretized Local integral equations:
where is a matrix of weight functions given by:
is a matrix that collects the derivatives of the weight:
The cubic spline functions are used as the test functions for the local weak form:
In which is the distance from node to point , and is the size of the influence domain for the weight function.
Collecting the discretized local integral equations together with the discretized boundary conditions for displacements, we get the complete system of algebraic equations for computation of nodal displacements, which are the Laplace transforms of ﬁctitious parameters .
The time dependent values of the transformed variables can be obtained by an inverse Laplace transform. There are many inversion methods available for the Laplace transformation. In the present analysis, the Stehfest algorithm  is used. If
is the Laplace-transform of
, an approximate value
for a speciﬁc time
The selected number N = 10 with a single precision arithmetic is optimal to receive accurate results. It means that for each time t, it is needed to solve N boundary value problems for the corresponding Laplace parameters force with . If M denotes the number of the time instants in which we are interested to know g(t), the number of the Laplace transform solutions is then M × N.
5. Numerical Results and Discussions
In this section, numerical results will be presented to illustrate the implementation and effectiveness of the proposed method. We present a numerical study for elastodynamic 2-D problem of a rectangular homogeneous isotropic plate  by using MLPG method, subjected to a dynamic force at the right end (Figure 2). A plane stress problem is considered, and a unit thickness is used. The dimensions of plate are length and height . The external excitation force , where (simple harmonic load with ) is a function of time, the damping coefficient is fixed for all numerical computation. A total number uniformly distributed nodes is used, as shown in Figure 3, to represent the problem domain.
The dimension of the quadrature domain is set to , where is the size of the quadrature domain, is a distance to the ﬁrst nearest neighbouring point from node i. The dimension of the inﬂuence domain is set to , where is the size of the support domain.
The isotropic rectangular plate analyzed with the following material properties: in Steel. Some important parameters on the performance of the method have been investigated.
Figure 4 displays the variations of displacements as a function of time at point B under the harmonic load for different values of size of support domain , where the time step is . It can be seen
Figure 2. Rectangular homogeneous isotropic plate subjected to a dynamic force at the right end of the plate.
Figure 3.Conﬁguration and nodal arrangement for the plate.
Figure 4. Displacements uy at the middle point B at the free end of the plate excited by the time-step load for different values of where and .
from this figure that the size of support domain inﬂuences on the results if , and has a small effect on the results if ). When , the results obtained by the present method is very good compared with the other authors  that have used the Newmark method. In the following analysis, is employed.
The time variation of displacements are given in Figure 5 with different values of size of the quadrature domain , where . It is found that the size of the quadrature domain inﬂuences seriously the results if and has a small effect on the results if .
Figure 5. The time variation of displacements uy for different values of where and .
Figure 6. The time variation of displacements in the y direction at the point B with different time steps using the fixed values .
However, if the size of the quadrature domain is too large, the result obtained for displacements by MLPG is great. We found good results with comparing with the results obtained by Long S. Y et al.  , they have used .
Many time steps are used in computation to check the stability of the presented MLPG formulation. The displacement variations as a function of time and results for diﬀerent time steps are plotted in Figure 6. It can be found that when the time step is less than 0.02 s, perfect results have been obtained using the Laplace transform comparing with the results obtained by other authors   that have used the Newmark method. It also can be found that when a time step is larger than 0.02 s, the results are not convergent and not accurate.
The present method MLPG that uses the cubic spline test function is used to analyze elastodynamic problem. The equation formulation based on MLPG method in Laplace transform and time domain with MLS approximation has been successfully implemented to solve elastodynamic problems in isotropic solids, subjected to a dynamic force at the right end of the plate. We found that the amplitude of the vibration decreases with time because the effects of damping and the harmonic load, the response should converge to the static deformation. We found that when the time step is less than 0.02 s, perfect results have been obtained by using the Laplace transform and when a time step is larger than 0.02 s the results are not accurate. We found that the size of the quadrature domain has a small effect on the results if , and the size of the support domain influences on the results if . The sizing parameters and , which decides the size of the subdomain needs to be chosen carefully, especially, in the dynamic analysis.