Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 5803457, 8 pages

http://dx.doi.org/10.1155/2016/5803457

## Complex Variable Meshless Manifold Method for Elastic Dynamic Problems

^{1}College of Mechanical and Electrical Engineering, Shandong Management University, Jinan 250357, China^{2}School of Mechanical and Automotive Engineering, Qilu University of Technology, Jinan 250353, China

Received 7 January 2016; Accepted 28 February 2016

Academic Editor: Elena Benvenuti

Copyright © 2016 Hongfen Gao and Gaofeng Wei. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Combining the finite covering technical and complex variable moving least square, the complex variable meshless manifold method can handle the discontinuous problem effectively. In this paper, the complex variable meshless method is applied to solve the problem of elastic dynamics, the complex variable meshless manifold method for dynamics is established, and the corresponding formula is derived. The numerical example shows that the numerical solutions are in good agreement with the analytical solution. The CVMMM for elastic dynamics and the discrete forms are correct and feasible. Compared with the traditional meshless manifold method, the CVMMM has higher accuracy in the same distribution of nodes.

#### 1. Introduction

The transient analysis is still an important part of the elastic mechanics in the elastic dynamic simulation. Because of the complexity of the problem, the numerical calculation has not been improved [1, 2]. Wilkins first studied dynamics of embedded wafer crack, elliptical crack, and semi elliptical surface crack with finite difference method (FDM) [3]. Because the finite difference method cannot describe the dynamic singular stress field, the method is not widely used. The finite element method (FEM) can simulate the dynamics problem, but the data preparation and calculation time are more lager.

In recent decades, the boundary element method (BEM) has been developed. It has the advantages of reducing dimension, saving computation time, and solving high accuracy and can deal with the problem of stress concentration and infinite and semi-infinite domains [4]. Therefore, the BEM has shown a greater potential than FEM and has been widely used in the field of elastic mechanics. Currently, the transient analysis of the elastic dynamic problem with BEM is mainly divided into time difference and spatial domain boundary element method, time domain method, and integral transformation method. J. Sládek and V. Sládek used Laplace transform method to study the buried wafer crack of the infinite elastic body under dynamic stress intensity factor of impact load [5]. Fedeliński studied the dynamic stress intensity factors by using the boundary element method [6].

Both the finite element method and the boundary element method are based on the element method. In the process of solving the problem of elastic dynamics, the accuracy of the solution can be influenced by the shape of the element when the element is extremely deformed [7, 8]. In recent years, with the research of the meshless methods, the meshless method is gradually taking a unique advantage in the study of dynamic problems.

There are many meshless methods, such as smoothed particle hydrodynamics method (SPH) [9, 10], the element-free Galerkin (EFG) method [11, 12], the collocation point method [13], the finite point method (FPM) [14, 15], reproducing kernel particle method (RKPM) [16], multiscale reproducing kernel particle method [17], the hp cloud method [18], the partition of unity finite element method [19, 20], generalized finite element method [21], meshless local Petrov-Galerkin method [22], least-square collocation meshless method [23], meshless weighted least-square method [24], complex variable meshless method (CVMM) [25, 26], and complex variable meshless manifold method (CVMMM) [27].

The CVMM is established by the complex variable moving least-squares (CVMLS) approximation [28, 29]. In the CVMLS approximation, the number of the unknown coefficients in the trial function is less than that of the MLS approximation, and the CVMLS approximation just needs fewer nodes in the domain than the MLS approximation. Then, combining CVMM method with manifold method, complex variable meshless manifold method (CVMMM) was presented and the method has been applied to the two-dimensional elastic problems, the fracture problems, bimaterial problems [30, 31], and so forth. This method has higher computational efficiency and precision than the traditional meshless method. In this paper, the dynamic problem will be researched through the CVMMM.

#### 2. CVMMM for Elastic Dynamic Problems

##### 2.1. Fundamental Equations of Elastic Dynamics

Setting the solving domain as , and the boundary as , the fundamental equation of the two-dimensional linear elastic dynamics is as follows.

(1) Equilibrium equation iswhere is the differential operator matrix, is the stress of any point in the domain , is the body force of any point in the domain , is mass density, and is camping coefficient; and represent the two-order derivative and first-order derivative of the displacement to time, respectively, and can be expressed as and represent inertia force and damping force, respectively.

(2) Geometrical equations arewhere and are the strain and displacement of any point in the domain , respectively, and can be given as

(3) Physical equation iswhere is elastic matrix.

(4) Boundary conditions arewhere is the known displacement on the displacement boundary , and is the known surface force on the surface force boundary . , , and is the boundary of the domain . Considerwhere and are the direction cosine of the outer normal for the boundary , respectively.

(5) Initial condition iswhere and are initial displacement and initial velocity of any node in the domain .

##### 2.2. Integral Weak Form of Elastic Dynamics

The equivalent integral weak form of elastic dynamics is

Substituting geometric equation (6) and physical equation (8) into (12), we have

##### 2.3. CVMMM for Elastic Dynamics

The displacement of any node at any time can be approximated by using node displacement in the influence domain. For CVMMM, the displacement of any node at any time within the domain can be expressed aswhere

By (14) and (16), we can obtainwhereand is the displacement of node at time ; that is,

The velocity and acceleration of any node at any time within the domain can be expressed aswhere and are the velocity and acceleration of node at time , respectively:

From (17), of (13) can be expressed aswhere and are the derivative of the shape function of CVMMM to and , respectively.

Substituting (17), (19), (20), and (23) into (13), we can obtain

Due to arbitrary of variation for node displacement, the system of equation can be obtained as follows:where , , and are the node acceleration vector, velocity vector, and displacement vector of the system, respectively; , , , and are the mass matrix, damping matrix, stiffness matrix, and load matrix of the system, respectively, and can be expressed as

##### 2.4. Processing of Boundary Conditions

In this paper, penalty function method is used to deal with displacement boundary condition. A new item is added to (13) to compulsorily satisfy the displacement boundary condition. The modified integral weak form can be written aswhere is the penalty factor, and the value is generally

When the direction of (or ) has a displacement constraint, the value of (or ) is correspondingly equal to 1; otherwise it is 0.

Substituting (17) into the last item of (28), we have

The discrete form iswhere

From (28) and discrete forms (26) and (31), the final equation of system can be obtained as

Let

Then (33) can be simplified as

If the damping is neglected, the equation of system can be simplified as

If the right side of (36) is equal to zero, the expression is a free vibration equation of the system.

#### 3. Time Integration Scheme

In the dynamic analysis, due to the inertia force and damping force appearing in the equilibrium equation, the mass matrix and damping matrix are introduced to the final solving equation, which makes the system discrete equation not algebraic equations, but ordinary differential equations.

For the two order ordinary differential equations, the solutions can be generally obtained by using the method of solving ordinary differential equations. However, in dynamic analysis of the meshless manifold method, the above method is not efficient because the order of the matrix is very high. The common method is the direct integral method and the mode superposition method in which the two methods are essentially equivalent. From the computational efficiency, the direct integral method is used in this paper.

There are two direct integration methods which are central difference method and Newmark method. The central difference method is a stable conditional algorithm and the time step must be less than a certain critical value which is decided by the solving problem; otherwise the algorithm will be unstable. Therefore, the central difference method is suitable for solving the problem of wave propagation, and not suitable for the structural dynamics. Because the low frequency component of the dynamic response of the structure is the main, a larger time step may allow to using from the calculation accuracy, and not necessary to limit the time step to smaller. As a result, a nonconditional stable implicit algorithm for the structural dynamics, which is mainly determined by the accuracy requirements, is usually adopted. The Newmark method is used in this paper.

The time domain is divided into time intervals, . The specific calculation is that a time interval is moving forward step by step; calculation of the next time is based on the previous time step. If the solutions of the moment of are obtained, the next computation is to solve the solution of the moment of .

Newmark integral method is essentially an extension of linear acceleration method. In the time step of , the recursive relation of displacement and velocity iswhere parameters of and are related to the stability and accuracy of the results. Different parameters are chosen to correspond to different integral methods. Linear acceleration method is as follows: Averaged acceleration method is as follows: Uniform acceleration method is as follows: Galerkin method is as follows: Backward difference method is as follows:

From (37), we have

Let

Equation (44) can be written as

Substituting (46) into (38), we have

In the Newmark method, the displacement solution at time of is obtained by solving (35); that is,

Substituting (46) and (47) into (48), we have

After being calculated, and can be obtained by (47) and (46), respectively.

#### 4. Numerical Example

The two-dimensional rectangular plate’s width is 2 m and length is 4 m; the free end is subjected to a suddenly applied uniform distributed load, as shown in Figure 1. Modulus of elasticity is , Poisson ratio is , and mass density is , without consideration of the influence of damping.