Research Article | Open Access

Volume 2016 |Article ID 5803457 | https://doi.org/10.1155/2016/5803457

Hongfen Gao, Gaofeng Wei, "Complex Variable Meshless Manifold Method for Elastic Dynamic Problems", Mathematical Problems in Engineering, vol. 2016, Article ID 5803457, 8 pages, 2016. https://doi.org/10.1155/2016/5803457

# Complex Variable Meshless Manifold Method for Elastic Dynamic Problems

Accepted28 Feb 2016
Published21 Mar 2016

#### 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) . 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 . 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 . Fedeliński studied the dynamic stress intensity factors by using the boundary element method .

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 , the finite point method (FPM) [14, 15], reproducing kernel particle method (RKPM) , multiscale reproducing kernel particle method , the hp cloud method , the partition of unity finite element method [19, 20], generalized finite element method , meshless local Petrov-Galerkin method , least-square collocation meshless method , meshless weighted least-square method , complex variable meshless method (CVMM) [25, 26], and complex variable meshless manifold method (CVMMM) .

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.

At the moment . The free end of the beam is subjected to a suddenly applied uniform distributed load , as shown in Figure 2.

Initial condition for the structure is , the analytical solution of the beam can be obtained by using the mode superposition method, and the axial displacement is

The axial stress iswhere is frequency, which is , is length of plate, is the cross-sectional area of plate, , and is the exciting force acting on the right side of the plate, which is .

In this example, the complex variable meshless manifold method is used to discrete the spatial domain and together with the Newmark time integration method is used to simulate the dynamics. We obtain the change of the displacement and stress of the point B at the free end versus time and the change of the displacement and stress of the central point C versus time. The numerical results are compared with the analytical solutions.

The distributions of regular nodes of are adopted, as shown in Figure 3. The linear basis function and Gauss weighted function are used in analysis, and the backward difference method of Newmark method is used in dynamic simulation; the time step is .

The relationships between the displacement of points B and C versus time are shown in Figures 4 and 5, the relationships between the stress of points B and C are given in Figures 6 and 7, respectively, and the numerical results are compared with the analytical solutions. From Figures 47, the numerical solutions are in good agreement with the analytical solutions. The results show that the CVMMM and the discrete forms are correct and the method is feasible. The figures show that the Newmark method does not have the ability to dissipate the high frequency, so the stresses have a dramatic oscillation phenomenon.

#### 5. Conclusion

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 results illustrate that the CVMMM and the discrete forms are correct and the method is feasible. Compared with the traditional meshless manifold method, the CVMMM has higher accuracy in the same distribution of nodes.

Because finite covering technology of manifold analysis is used, the CVMMM can handle the domain of continuous and discontinuous problems and can overcome the mapping process of long crack in finite element method or meshless method. The accuracy of numerical solution is improved.

The CVMMM can be applied to the problems such as crack propagation or physical nonlinearity and geometric nonlinearity and can be further used to solve practical engineering problems.

#### Competing Interests

The authors declare that they have no competing interests.

#### Acknowledgments

The work is supported by the National Natural Science Foundation of China through Grant 11271234 and the Science Program of Shandong Provincial Education Department through Grant J13LJ51.

1. K. Arakawa, T. Mada, and K. Takahashi, “Correlations among dynamic stress intensity factor, crack velocity and acceleration in brittle fracture,” International Journal of Fracture, vol. 105, no. 4, pp. 311–320, 2000. View at: Publisher Site | Google Scholar
2. J.-B. Leblond and P. Mouro, “Crack propagation from a pre-existing flaw at a notch root. I. Introduction and general form of the stress intensity factors at the initial crack tip,” International Journal of Fracture, vol. 104, no. 3, pp. 211–224, 2000. View at: Google Scholar
3. M. L. Wilkins, “Numerical analysis of dynamic crack problems,” in Mechanics of Fracture, G. C. Sih, Ed., vol. 4, P. Noordhoff, Leyden, Netherlands, 1977. View at: Google Scholar
4. M. H. Aliabadi, “Boundary element formulations in fracture mechanics,” Applied Mechanics Reviews, vol. 50, no. 2, pp. 83–96, 1997. View at: Publisher Site | Google Scholar
5. J. Sládek and V. Sládek, “Dynamic stress intensity factors studied by boundary integro-differential equations,” International Journal for Numerical Methods in Engineering, vol. 23, no. 5, pp. 919–928, 1986. View at: Publisher Site | Google Scholar
6. P. Fedeliński, “Boundary element method in dynamic analysis of structures with cracks,” Engineering Analysis with Boundary Elements, vol. 28, no. 9, pp. 1135–1147, 2004. View at: Publisher Site | Google Scholar
7. O. C. Zienkiewicz, The Finite Element Method, McGraw-Hill, 3rd edition, 1977.
8. C. A. Brebbia, The Boundary Element Method for Engineers, Pentech Press, London, UK, 1978. View at: MathSciNet
9. L. B. Lucy, “A numerical approach to the testing of the fission hypothesis,” The Astronomical Journal, vol. 8, no. 12, pp. 1013–1024, 1977. View at: Publisher Site | Google Scholar
10. R. A. Gingold and J. J. Moraghan, “Smoothed particle hydrodynamics—theory and application to non-spherical stars,” Monthly Notices of the Royal Astronomical Society, vol. 181, pp. 375–389, 1977. View at: Publisher Site | Google Scholar
11. T. Belytschko, Y. Y. Lu, and L. Gu, “Element-free Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994. View at: Publisher Site | Google Scholar | MathSciNet
12. H.-J. Chung and T. Belytschko, “An error estimate in the EFG method,” Computational Mechanics, vol. 21, no. 2, pp. 91–100, 1998. View at: Publisher Site | Google Scholar | MathSciNet
13. E. Oñate, S. Idelsohn, O. C. Zienkiewicz, and R. L. Taylor, “A finite point method in computational mechanics. Applications to convective transport and fluid flow,” International Journal for Numerical Methods in Engineering, vol. 39, no. 22, pp. 3839–3866, 1996. View at: Publisher Site | Google Scholar
14. E. Onate and S. Idelsohn, “A mesh-free finite point method for advective-diffusive transport and fluid flow problems,” Computational Mechanics, vol. 21, no. 4-5, pp. 283–292, 1998. View at: Publisher Site | Google Scholar | MathSciNet
15. H. Fatahi, J. Saberi-Nadjafi, and E. Shivanian, “A new spectral meshless radial point interpolation (SMRPI) method for the two-dimensional Fredholm integral equations on general domains with error analysis,” Journal of Computational and Applied Mathematics, vol. 294, pp. 196–209, 2016. View at: Publisher Site | Google Scholar | MathSciNet
16. W. K. Liu, S. Jun, and Y. F. Zhang, “Reproducing kernel particle methods,” International Journal for Numerical Methods in Fluids, vol. 20, no. 8-9, pp. 1081–1106, 1995. View at: Publisher Site | Google Scholar | MathSciNet
17. W. K. Liu, S. Jun, D. T. Sihling, Y. Chen, and W. Hao, “Multiresolution reproducing kernel particle methods for computational fluid dynamics,” International Journal for Numerical Methods in Fluid, vol. 24, no. 12, pp. 1391–1415, 1997. View at: Publisher Site | Google Scholar
18. C. A. Duarte and J. T. Oden, “H-p clouds—an h-p meshless method,” Numerical Methods for Partial Differential Equations, vol. 12, no. 6, pp. 673–705, 1996. View at: Publisher Site | Google Scholar
19. J. M. Melenk and I. Babuška, “The partition of unity finite element method: basic theory and applications,” Computer Methods in Applied Mechanics and Engineering, vol. 139, no. 1–4, pp. 289–314, 1996. View at: Publisher Site | Google Scholar | MathSciNet
20. Y. C. Cai, L. Han, L. G. Tian, and L. Y. Zhang, “Meshless method based on Shepard function and partition of unity for two-dimensional crack problems,” Engineering Analysis with Boundary Elements, vol. 65, pp. 126–135, 2016. View at: Publisher Site | Google Scholar | MathSciNet
21. J. T. Oden, A. Duarte, and O. C. Zienkiewicz, “A new cloud-based hp finite element method,” International Journal for Numerical Methods in Engineering, vol. 50, pp. 160–170, 1998. View at: Google Scholar
22. S. N. Atluri and T. Zhu, “A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics,” Computational Mechanics, vol. 22, no. 2, pp. 117–127, 1998. View at: Publisher Site | Google Scholar | MathSciNet
23. X. Zhang, X.-H. Liu, K.-Z. Song, and M.-W. Lu, “Least-squares collocation meshless method,” International Journal for Numerical Methods in Engineering, vol. 51, no. 9, pp. 1089–1100, 2001. View at: Publisher Site | Google Scholar | MathSciNet
24. K. M. Liew, J. Ren, and S. Kitipornchai, “Analysis of the pseudoelastic behavior of a SMA beam by the element-free Galerkin method,” Engineering Analysis with Boundary Elements, vol. 28, no. 5, pp. 497–507, 2004.
25. Y. M. Cheng and J. H. Li, “A meshless method with complex variables for elasticity,” Acta Physica Sinica, vol. 54, no. 10, pp. 4463–4471, 2005. View at: Google Scholar | MathSciNet
26. Y. J. Deng, C. Liu, M. J. Peng, and Y. M. Cheng, “The interpolating complex variable element-free galerkin method for temperature field problems,” International Journal of Applied Mechanics, vol. 7, no. 2, Article ID 1550017, 2015. View at: Publisher Site | Google Scholar
27. H. F. Gao and Y. M. Cheng, “A complex variable meshless manifold method for elasticity problems,” Acta Mechanica Sinica, vol. 41, no. 3, pp. 1–9, 2009. View at: Google Scholar
28. K. M. Liew, C. Feng, Y. Cheng, and S. Kitipornchai, “Complex variable moving least-squares method: a meshless approximation technique,” International Journal for Numerical Methods in Engineering, vol. 70, no. 1, pp. 46–70, 2007. View at: Publisher Site | Google Scholar | MathSciNet
29. F. X. Sun, J. F. Wang, Y. M. Cheng, and A. X. Huang, “Error estimates for the interpolating moving least-squares method in n-dimensional space,” Applied Numerical Mathematics, vol. 98, pp. 79–105, 2015. View at: Publisher Site | Google Scholar | MathSciNet
30. H. F. Gao and Y. M. Cheng, “A complex variable meshless manifold method for fracture problems,” International Journal of Computational Methods, vol. 7, no. 1, pp. 55–81, 2010. View at: Publisher Site | Google Scholar | MathSciNet
31. H. F. Gao and G. F. Wei, “Stress intensity factor for interface cracks in bimaterials using complex variable meshless manifold method,” Mathematical Problems in Engineering, vol. 2014, Article ID 353472, 8 pages, 2014. View at: Publisher Site | Google Scholar