The local Petrov-Galerkin methods (MLPG) have attracted much attention due to their great flexibility in dealing with numerical model in elasticity problems. It is derived from the local weak form (WF) of the equilibrium equations and by inducting the moving last square approach for trial and test functions in (WF) is discussed over local sub-domain. In this paper, we studied the effect of the configuration parameters of the size of the support or quadrature domain, and the effect of the size of the cells with nodes distribution number on the accuracy of the methods. It also presents a comparison of the results for the Shear stress, the deflections and the error in energy.

Recently Meshless formulations are becoming popular due to their higher adaptivity and lower cost for preparing input data in the numerical analysis. A variety of meshless methods has been proposed so far (Belytschko et al., 1994; Atluri and Shen, 2002; Liu, 2003; Atluri, 2004) [1-6]. Many of them are derived from a weak-form formulation on global domain [

The meshless local Petrov-Galerkin (MLPG) method originated by Atluri and Zhu [

The method is a fundamental base for the derivation of many meshless formulations, since trial and test functions are chosen from different functional spaces.

MLPG does not need a global mesh for either function approximation or integration. The procedure is quite similar to numerical methods based on the strong-form formulation, such as the finite difference method (FDM). However, because in the MLPG implementation, moving least squares (MLS) approximation is employed for constructing shape functions, special treatments are needed to enforce the essential boundary conditions [4,7].

The aims of this paper are to study the effect on accuracy and convergence of MLPG methods of different size parameters: and associated to support and quadrature domains respectively. The support domain is denoted be equal to influence domain. For fixed values of: and, the effect of cells numbers with nodes distribution number, on energy errors is also studied and some of our results are presented.

In this work, the MLPG method will be developed for solving the problem of a thin elastic homogenous plate. The discretization and numerical implementation are presented in Section 2 numerical example for 2D problem are given in Section 3. Then paper ends with discussions and conclusions.

Let us consider a two-dimensional problem of solid mechanics in domain bounded by whose strongform of governing equation and the essential boundary conditions are given by:

on (2)

on (3)

where in, is the stress vector and the body force vector.

On the natural boundaries is the prescribed traction, denoted the vector of unit outward normal at a point.

the displacement components in the plan and on the essential boundaries.

In the local Petrov-Galerkin approaches [

where is the local domain of quadrature for node I and is the part of the essential boundary that intersect with the quadrature domain. is the weight or test function , [

is the penalty factor, Here we use the same penalty factor for all the displacement constraint equations (essential boundary conditions) [

Generally, in meshfree methods, the representation of field nodes in the domain will be associated to other repartitions of problem domain: influence domain for nodes interpolation, is the support domain for accuracy. For each node is the weight function domain, and is the quadrature domain for local integration.

Using the divergence theorem [

where

: The internal boundary of the quadrature domain

: The part of the natural boundary that intersects with the quadrature domain

: The part of the essential boundary that intersects with the quadrature domain When the quadrature domain is located entirely within the global domain on and no boundary conditions are specified then.

Unlike the Galerkin method, the Petrov-Galerkin method chooses the trial and test functions from different spaces. The weight function is purposely selected in such a way that it vanishes on. We can then change the expression of Equation (5):

Witch is the local Petrov-Galerkin weak form. Here we require [3,11] and the simplified Petrov-Galerkin form is:

Precedent equations are used to establish the discrete equations for all the nodes whose quadrature domain falls entirely within the problem domain (Equation (7)) and to establish the discrete equations for all the boundary nodes or the nodes whose quadrature domain intersects with the problem boundary “Equation (6)”.

To approximate the distribution of the function in the support domain over a number of nodes. We shall have the approximant of [

where I denote the set of the nodes in the support domain of point.

the MLS shape function for node I that is created using nodes in the support domain of point. The discrete system in Equation (6) is given in matrix form:

where is a matrix that collects the derivatives of the weight functions in Equation (6), and

is the matrix of weight function. The stress vector defined by:

where is the symmetric elasticity tensor of the material

Substituting the differential operator

and Equation (8) into Equation (10) we obtain:

where and by using

the tractions of a point x can be written as:

Substituting Equations (8), (11) and (12) into Equation (6), we obtain the discrete systems of linear equations for the node I.

That can assembled in matrix form:

where nodal stiffness matrix

And nodal force vector with contributions from body forces applied in the problem domain, tractions applied on the natural boundary, as well as the penalty force terms.

Two independent linear equations can be obtained for each node in the entire problem domain and assembled all these equations to obtain the final global system equations:

To solve the precedent system, the standard Gauss quadrature formula is applied with 16 Gauss points [3,14] for evaluation of boundary and domain integrals in Equations (15) and (16)

In this section, numerical results are presented for Cantilever rectangular plate in

The problem data:

The height of the beam and the length of the beam:

The thickness of the plat: and Loading (integration of the distributed traction):

Young’s modulus: and Poisson’s ratio:

The standard Gaussian quadrature formula is applied with 16 Gauss points, and for MLS approximation linear polynomial basis functions are applied, the cubic spline function is used as the test function for the local Petrov-

Galerkin weak-form. In our numerical calculations we consider many regular distributions of nodes: 55 or 175. To calculate the error energy a background cells is required, then we have varying the number of cell. To obtain the distribution of the deflection and stress through the plates, size of quadrature domain and support domain are varied. Nodal configuration for a cantilever plate with 55 nodes (

The accuracy is clear for the greater value of field nodes distribution. It is also shown in this figure, on the cross-section the meshless MLPG agree well with those from analytical solution (dashed lines).