Advanced Materials and Technologies for Structural Performance Improvement
View this Special IssueResearch Article  Open Access
A CellBased Smoothed XFEM for Fracture in Piezoelectric Materials
Abstract
This paper presents a cellbased smoothed extended finite element method (CSXFEM) to analyze fractures in piezoelectric materials. The method, which combines the cellbased smoothed finite element method (CSFEM) and the extended finite element method (XFEM), shows advantages of both methods. The crack tip enrichment functions are specially derived to represent the characteristics of the displacement field and electric field around the crack tip in piezoelectric materials. With the help of the smoothing technique, integrating the singular derivatives of the crack tip enrichment functions is avoided by transforming interior integration into boundary integration. This is a significant advantage over XFEM. Numerical examples are presented to highlight the accuracy of the proposed CSXFEM with the analytical solutions and the XFEM results.
1. Introduction
Because of their inherent coupling of electric and mechanical behaviors, piezoelectric materials have been widely used in sensors, actuators, signal transmitters and surface acoustic wave devices, aerospace panels, and civil structures. In those applications, piezoelectric materials may experience high mechanical stresses and electric field concentrations. As a result they may fail due to dielectric breakdown or fractures. These materials are usually inhomogeneous and brittle, with low ultimate tensile strength and fracture toughness. Therefore, defects such as cracks and voids should be detected to ensure reliability and durability of the piezoelectric structures. Numerical simulation of fractures in piezoelectric ceramics was conducted in [1], using finite element method (FEM) [2, 3], boundary element method (BEM) [4, 5], meshless method [6], and extended finite element method (XFEM) [7–9].
The theoretical fundamentals of piezoelectric fracture mechanics for cracks were presented in [10, 11]. The analytical work to investigate the fracture mechanics of piezoelectric structures was based on Lekhnitskii and Stroh formalism [12]. An elliptic hole with a major axis perpendicular to the polarization direction inside piezoelectric ceramic and the field variables around the cavity were studied in [13]. An overview and critical discussion about the present state in the field of piezoelectric fracture mechanics were given in [14]. A survey on using numerical methods of crack analyses in piezoelectric medium along with FEM to solve fracture parameters is presented in [15]. Recently, XFEM was applied to analyze the 2D crack problems and fully coupled piezoelectric effect of piezoelectric ceramics [16]. The integral and integral were used to solve the stress intensity factors and electric displacements intensity factors [17, 18]. The newly developed crack tip enrichment functions of XFEM were found suitable for cracks in piezoelectric materials [19]. An extension of XFEM for dynamic fracture in piezoelectric materials was presented in [20].
In recent years, the smoothed FEM has been well adopted to solve fracture mechanics problem [21–27]. Based on generalized gradient smoothing technique a lot of novel and powerful numerical methods have been developed [28–36]. A cellbased smoothed finite element method (CSFEM) has been developed with the smoothing domains constructed based on cell of the elements. In the method line integration was used along the boundaries of the smoothing cells instead of area integration. Moreover, CSFEM does not need mapping and derivatives. It also showed that the results are less sensitive to distorted elements.
In terms of the advantages, smoothing technique was incorporated into the extended FEM [37–39]. An edgebased smoothed XFEM was developed to combine the advantages of the edgebased smoothed FEM and the XFEM [40]. A nodebased smoothed XFEM was applied to linear elastic fracture mechanics [41].
Recently, extended finite element method, collocation boundary element, cellbased smoothed finite element method, and so forth were used to solve the problem of fracture in piezoelectric materials. The extended finite element method focuses on the definition of new enrichment functions suitable for cracks in piezoelectric structures and generalized domain integrals are used for the determination of crack tip parameters [42]. The collocation boundary element with subdomain technique is developed, whereby the fundamental solutions are computed by a fast numerical algorithm applying Fourier series [43]. The cellbased smoothed finite element method and VCCT have been used to simulate the fracture mechanics of piezoelectric materials. A piezoelectric element tailored for VCCT was used to study the crack of piezoelectric materials. CSFEM and VCCT were introduced into fracture mechanics of piezoelectric materials and CSFEMVCCT for piezoelectric material with cracks was put forward [44]. In this paper cellbased smoothed XFEM (CSXFEM) is extended to simulate flaws in piezoelectric structures. CSXFEM combining characteristics of extended finite element method and smoothed finite element method can improve the accuracy of XFEM. Nodal enrichment can model crack propagation without remeshing. CSXFEM has advantages that the crack tip element does not need fine division; the shape function is simple and no derivatives of shape functions are needed.
This paper is outlined as follows. In Section 2, the governing equations of piezoelectric materials are introduced. Section 3 focuses on the formulation of cellbased smoothed extended finite element method. Section 4 presents the electromechanical integral for 2D crack analysis. In Section 5, numerical examples with the assumption of impermeable crack face boundary conditions are presented to demonstrate the accuracy and efficiency of CSXFEM. Section 6 is the conclusion.
2. Governing Equations
The electroelastic response of a piezoelectric body of volume and regular boundary surface is governed by the mechanical and electrostatic equilibrium equations:where is mechanical body force, is electric body charge, is the symmetric Cauchy stress tensor, and is electric displacement vector components.
The constitutive equations for a twodimensional piezoelectric material in the  plane can be expressed in terms of the stress and the electric field:where , , and are the strain tensor, the electric displacement vector, and the electric field vector, respectively; , , and denote elastic stiffness at constant electric field, piezoelectric constants, and dielectric permittivity at constant strain, respectively.
The strains tensor are related to displacements byand the electric field vector is related to electric potential byThe piezoelectric body could be subjected to the following essential and natural boundary conditions (Figure 1).
Essential Boundary Conditions. Consider
Natural Boundary Conditions. Considerwhere , , , , and are mechanical displacement, electric potential, surface force components, surface charge, and outward unit normal vector components, respectively. The crack faces and are considered tractionfree. The cracks are assumed to be electrically impermeable. Nevertheless, extension to limited permeable cracks is possible.
The twodimensional matrix form of the mechanical and electrical constitutive equations can be given by [15]where are the elastic compliance constants, are piezoelectric constants, and are the dielectric constants.
3. CellBased Smoothed Extended Finite Element Method
In CSXFEM, the approximation of displacement and electric potential field in a piezoelectric material are given bywhere , , and are shape functions of the nodal displacement, while are the nodal degrees of freedom associated with node , and and are additional nodal degrees of freedom corresponding to the Heaviside function and the neartip functions, , respectively. , , and are shape functions of the nodal electric potential, while are the nodal degrees of freedom associated with node , and and are additional nodal degrees of freedom corresponding to the Heaviside function and the neartip functions, , respectively.
Nodes in set have supports split by crack and nodes in set which belong to the smoothing domains contain a crack tip. These nodes are enriched with the Heaviside and asymptotic branch function fields depicted with squares and circles, respectively. The support domain of is associated with nodes of CSFEM, shown in Figure 2.
As illustrated in Figure 2, fournode quadrilateral cells are used for cellbased strain smoothing operation. The meshing characteristics of CSXFEM are consistent with XFEM. The complex structure can adopt the fine mesh. For CSFEM, the number of subcells (SC) would affect the performance of the results. In the case that the solution of SFEM (SC = 1) is overestimated to the exact solution, there exists one optimal value SC > 1 (normally SC = 4) which gives the best results as compared to the exact ones. In the case that the solution of SFEM (SC = 1) is underestimated to the exact solution, it is suggested that we should use SC = 2 to obtain the solution. This solution will be stable and have the smallest displacement and energy norms. In practical calculation, we can use SC = 4 for all problems. The results will always be better than that of standard FEM. And in many cases (not all) this solution (SC = 1) is closest to the exact solution [21]. So we can use four subcells for each quadrilateral.
In CSXFEM, Heaviside enriched degrees of freedom are added to nodes in whose support domain is split by the crack and tip enriched degrees of freedom are added to nodes in set whose support domain contains the crack tip. In order to keep the convergence rate as high as possible, a socalled geometric enrichment which is independent from the discretization is used [45–47]:where is a point on the crack surface; see Figure 3.
For piezoelectric problems, it is advisable to use the regular enrichment functions stemming from the isotropic elasticity. It should be mentioned that similar results have been obtained independently with alternative enrichment functions for cracks in confined plasticity problems [42]. The neartip enrichment consists of functions which incorporate the radial and angular behaviors of the twodimensional asymptotic crack tip displacement field [37, 48]:where and are polar coordinates in the local crack tip coordinate system; see Figure 4.
Employing the strain smoothing operation, the smoothed strain in the domain from the displacement approximation in (8) can be written in the following matrix form:where is the smoothed strain gradient matrix for the standard CSFEM part; and correspond to the enriched parts of the smoothed strain gradient matrix associated with the Heaviside and branch functions, respectively.
When employing the electric field smoothing operation, the smoothed electric field in the domain from the displacement approximation in (9) can be written in the following matrix form:where is the smoothed electric field gradient matrix for the standard CSFEM part; and correspond to the enriched parts of the smoothed electric field gradient matrix associated with the Heaviside and branch functions, respectively. These matrixes can be written as follows:wherewhere is the number of segments of the boundary , is the number of Gauss points used in each segment, is the corresponding Gauss weights, and are the outward unit normal components to each segment on the smoothing domain boundary, and is the th Gaussian point on the th segment of the boundary .
The standard discrete system of equations is obtained:where
4. Electromechanical Integral
According to research by Rice [49] and Eshelby [50] on pure mechanical applications,Equation (18) describes the material force when an electromechanically loaded domain is virtually displaced by the vector . The term in brackets is called the piezoelectric energy momentum tensor:analogous to the energy momentum tensor, which was introduced by Eshelby [50]. If the integration path contains neither defects nor source terms and , the force vanishes and . Now we consider a path enclosing the tip of a crack; see Figure 5. The electromechanical integral vector is the material force associated with the crack tip singularity. It is defined as the limit when is shrunk towards
For the numerical analysis of the integral, an equivalent domain integral is easier to handle than a line integral. The condition for the transformation into a domain integral applying the Gaussian integral theorem, a closed integration path with an outward pointing normal vector is given, as can be seen in Figure 5. Then (21) can be written as follows:Now a weighting function is introduced, which should be continuous and satisfy the conditionThe calculation of the integral is carried out in a ring of elements surrounding the crack tip. The elements within the ring move as a rigid body. is a constant in these elements; so the derivative of with respect to is zero. For the elements outside the ring, is zero, and again the derivative of is zero. For the elements belonging to the ring, the vector is 1. If the weighting function is introduced in (22), the contribution of the integral along disappears and The transformation into an equivalent domain integral leads to
If we use the expression of for (20), we get the 2dimensional integral as a domain integral: If the material is homogenous, no volume forces or charges and no crack face tractions or charges are applied, and the crack face normal points to the directions; that is,Equation (26) can be reduced asThe component of the electromechanical integral vector has a physical meaning of the energy release rate .
The mechanical stresses and the electric fluxes behave singular as , whereas the electric potential and the mechanical displacements show a parabolic shape ~. The angular functions , , , and depend only on material constants. The coefficients ,, and are the wellknown mechanical stress intensity factors, which are complemented by the new forth “electric intensity factor” , that characterizes the electric field singularity. Mutual interdependence between mechanical and electrical crack tip parameters can be given by
If the limit to an infinitesimal crack growth is considered taking place inside the neartip solution, its relationship with the intensity factors can be found in [6]. Consider The generalized Irwinmatrix depends on the elastic, piezoelectric, and dielectric material constants and the relative orientation of the crack with respect to material’s axes and polarization vectors. For a special case when the crack lies perpendicular to the polarization, the crack tip field and the coefficients of were determined in [10, 11]. Then is reduced to the following expression, where, , , , , and are material constants:The energy release rate consists of the mechanical terms (for each opening mode) and the electric contribution.
5. Numerical Examples
5.1. Electromechanical Griffith Crack (Uniaxial Load)
In order to test the accuracy of the CSXFEM for crack analysis, the methods were applied to a crack in a plane subjected to normal uniaxial tension when MPa and electrical flux C/m^{2}. Figure 6 shows the polarization direction as . The distance of the central crack along the direction is and side length of the plate is 10 m. At the crack faces, impermeable electric boundary conditions are prescribed. Piezoelectric materials PZT4, P7, and PZTH5 were adopted for numerical simulation. The material parameters are shown in Table 1.

The exact analytical solution for the crack in the infinite plane under farfield loads , , and was given by Pak [11]. ConsiderIn order to verify the reliability of CSXFEM, we set crack length as m and use two grid models ((I) 2288 elements, (II) 1521 elements) before the fracture starts, as shown in Figure 7. There are five types of smoothing elements ((i) split smoothing element: there is one Gauss point on each boundary segment for split smoothing element; (ii) splitblending smoothing element: one Gauss point on each boundary segment is sufficient; (iii) tip smoothing element: five Gauss points on a segment of smoothing element are sufficient; (iv) tipblending smoothing element: five Gauss points on each boundary segment are sufficient; (v) standard smoothing element: one Gauss point on each boundary segment is sufficient) being used for numerical integration as mentioned in [41]. In the simulation every 4node calculation grid adopts four smoothing elements. The results are compared with those of the XFEM and the theoretical solutions.
The results of stress intensity factors and electric displacements factors produced by CSXFEM and XFEM for three loading cases are listed in Table 2. It can be seen that the results of CSXFEM are closer to the analytical solution than those of XFEM when using the same mesh. This confirms that combining cellbased smoothing technique with the XFEM can improve accuracy.

The normalized mechanical and electrical intensity factors of PZT4 and P7 with different length of cracks and model (I) grids are listed in Table 3. It is obvious that the CSXFEM can produce more accurate results than XFEM when using the same number of nodes, which indicates that the smoothing technique adopted in this work improve the calculation of normalized mechanical and electrical intensity factors for fracture in piezoelectric materials.

5.2. Electromechanical Griffith Crack (Shear Load)
In the second example of straight electromechanical crack in the plane, the normal mechanical load is replaced by a shear load. The farfield loadings are , MPa, and C/m^{2}; see Figure 8.
From Table 4, we can observe that the result of the mechanical and electrical intensity factors of PZTH5 has a high precision in two models. The accuracy of CSXFEM is higher than that of XFEM. The results prove that the CSXFEM can decrease stiffness of the system and improve solution accuracy. Also model (I) of PZTH5 was used under the CSXFEM and XFEM using gauss integral calculation efficiency. CSXFEM takes 46.247 seconds, while XFEM takes 48.252 second with the following CPU setup: Intel Core i53470 3.20 GHz, RAM: 8 G. The efficiency of CSXFEM has been improved but is not obvious.

The normalized mechanical and electrical intensity factors of PZT4 and P7 with different length of cracks when using model (I) grids are listed in Table 5. It is obvious that the CSXFEM can produce more accurate results than XFEM using the same number of nodes.

5.3. Piezoelectric Model with a Hole with Cracks
A piezoelectric model, with a center circular hole mm and horizontal cracks on the left and right, is subjected to unidirectional uniform tensile MPa and electric displacement C/m^{2} at infinity. The geometric model is expressed as in Figure 9, the length of side mm, is the direction of polarization, and material is PZTH5. The situation of meshing when mm was given by Figure 10.
The stress intensity factor and electric displacement intensity factor on the left side of the crack tip obtained by CSXFEM and XFEM in different crack lengths are shown in Table 6. The maximum error of calculation in CSXFEM is 2.9%, and the maximum error of calculation in XFEM is 4.1%. The results verified the accuracy of CSXFEM.

6. Conclusions
With the growing applications of piezoelectric structures in innovative technical areas, problems of strength and reliability become important and have to be carefully investigated. In order to quantitatively assess fracture and fatigue, sophisticated analysis of cracks’ electromechanical properties is needed. The fracture mechanics approach for cracklike defects in piezoelectric materials reveals coupled electrical and mechanical field singularities. Effective numerical methods are needed to evaluate fracture behavior of cracks in arbitrary piezoelectric structures subjected to combined electromechanical loading.
In this paper CSXFEM with crack tip enrichment functions was presented for electromechanical crack analyses. Two examples are used to verify the accuracy of CSXFEM. Through the numerical simulation conclusions can be drawn as follows:(1)The cellbased smoothing technique is extended into XFEM to combine the advantages of XFEM and CSFEM. Thus, CSXFEM has high accuracy and convergence rate. The present method also simplifies the integration of discontinuous approximation by transforming interior integration into boundary integration. More importantly, no derivatives of shape functions are needed to compute the stiffness matrix.(2)The examples show that CSXFEM can perform better in accuracy and convergence rate compared with XFEM in terms of stress intensity factors and electric displacements factors.(3)CSXFEM is superior to XFEM in regard to the calculation of the stiffness matrix and singularity in the integrand. It also avoids the mapping process, which will increase the complexity of the calculation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was financially supported by the National Natural Science Foundation of China (Grant no. 51305157), the National Key Scientific Instrument and Equipment Development Projects, China (Grant no. 2012YQ030075), and Jilin Provincial Department of Science and Technology Fund Project (Grant no. 20130305006GX).
References
 M. Li, J. X. Yuan, D. Guan, and W. M. Chen, “Application of piezoelectric fiber composite actuator to aircraft wing for aerodynamic performance improvement,” Science China Technological Sciences, vol. 54, no. 2, pp. 395–402, 2011. View at: Publisher Site  Google Scholar
 M. Kuna, “Finite element analyses of crack problems in piezoelectric structures,” Computational Materials Science, vol. 13, no. 1–3, pp. 67–80, 1998. View at: Publisher Site  Google Scholar
 J. Meckerle, “Smart materials and structures—a finite element approach—an addendum: a bibliography (1997–2002),” Modelling and Simulation in Materials Science and Engineering, vol. 11, no. 5, pp. 707–744, 2003. View at: Google Scholar
 F. GarcíaSánchez, Ch. Zhang, and A. Sáez, “2D transient dynamic analysis of cracked piezoelectric solids by a timedomain BEM,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 3340, pp. 3108–3121, 2008. View at: Publisher Site  Google Scholar
 Ch. Zhang, F. GarcíaSánchez, and A. Sáez, “Timedomain BEM analysis of cracked piezoelectric solids under impact loading,” in Computational Mechanics: Proceedings of “International Symposium on Computational Mechanics” July 30–August 1, 2007, Beijing, China, pp. 206–218, Springer, Berlin, Germany, 2009. View at: Publisher Site  Google Scholar
 J. Sladek, V. Sladek, C. Zhang, P. Solek, and L. Starek, “Fracture analyses in continuously nonhomogeneous piezoelectric solids by the MLPG,” Computer Modeling in Engineering and Sciences, vol. 19, no. 3, pp. 247–262, 2007. View at: Google Scholar  MathSciNet
 S. Nanthakumar, T. Lahmer, X. Zhuang, G. Zi, and T. Rabczuk, “Detection of material interfaces using a regularized level set method in piezoelectric structures,” Inverse Problems in Science and Engineering, pp. 1–24, 2015. View at: Publisher Site  Google Scholar
 S. S. Nanthakumar, T. Lahmer, and T. Rabczuk, “Detection of flaws in piezoelectric structures using extended FEM,” International Journal for Numerical Methods in Engineering, vol. 96, no. 6, pp. 373–389, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Nanthakumar, T. Lahmer, and T. Rabczuk, “Detection of multiple flaws in piezoelectric structures using XFEM and level sets,” Computer Methods in Applied Mechanics and Engineering, vol. 275, pp. 98–112, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Suo, C.M. Kuo, D. M. Barnett, and J. R. Willis, “Fracture mechanics for piezoelectric ceramics,” Journal of the Mechanics and Physics of Solids, vol. 40, no. 4, pp. 739–765, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. E. Pak, “Linear electroelastic fracture mechanics of piezoelectric materials,” International Journal of Fracture, vol. 54, no. 1, pp. 79–100, 1992. View at: Publisher Site  Google Scholar
 X.L. Xu and R. K. N. D. Rajapakse, “Analytical solution for an arbitrarily oriented void/crack and fracture of piezoceramics,” Acta Materialia, vol. 47, no. 6, pp. 1735–1747, 1999. View at: Publisher Site  Google Scholar
 H. Sosa, “Plane problems in piezoelectric media with defects,” International Journal of Solids and Structures, vol. 28, no. 4, pp. 491–505, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Kuna, “Fracture mechanics of piezoelectric materials—where are we right now?” Engineering Fracture Mechanics, vol. 77, no. 2, pp. 309–326, 2010. View at: Publisher Site  Google Scholar
 M. Kuna, “Finite element analyses of cracks in piezoelectric structures: a survey,” Archive of Applied Mechanics, vol. 76, no. 1112, pp. 725–745, 2006. View at: Publisher Site  Google Scholar
 R. R. Bhargava and K. Sharma, “XFEM studies on an inclined crack in a 2D finite piezoelectric media,” in Future Communication, Computing, Control and Management, vol. 141 of Lecture Notes in Electrical Engineering, pp. 285–290, Springer, Berlin, Germany, 2012. View at: Publisher Site  Google Scholar
 M. Abendroth, U. Groh, M. Kuna, and A. Ricoeur, “Finite elementcomputation of the electromechanical Jintegral for 2D and 3D crack analysis,” International Journal of Fracture, vol. 114, no. 4, pp. 359–378, 2002. View at: Publisher Site  Google Scholar
 Y. Motola and L. BanksSills, “Mintegral for calculating intensity factors of cracked piezoelectric materials using the exact boundary conditions,” Journal of Applied Mechanics, vol. 76, no. 1, Article ID 011004, 9 pages, 2009. View at: Publisher Site  Google Scholar
 R. R. Bhargava and K. Sharma, “A study of finite size effects on cracked 2D piezoelectric media using extended finite element method,” Computational Materials Science, vol. 50, no. 6, pp. 1834–1845, 2011. View at: Publisher Site  Google Scholar
 H. NguyenVinh, I. Bakar, M. A. Msekh et al., “Extended finite element method for dynamic fracture of piezoelectric materials,” Engineering Fracture Mechanics, vol. 92, pp. 19–31, 2012. View at: Publisher Site  Google Scholar
 G. R. Liu, T. T. Nguyen, K. Y. Dai, and K. Y. Lam, “Theoretical aspects of the smoothed finite element method (SFEM),” International Journal for Numerical Methods in Engineering, vol. 71, no. 8, pp. 902–930, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 T. NguyenThoi, T. BuiXuan, P. PhungVan, H. NguyenXuan, and P. NgoThanh, “Static, free vibration and buckling analyses of stiffened plates by CSFEMDSG3 using triangular elements,” Computers & Structures, vol. 125, pp. 100–113, 2013. View at: Publisher Site  Google Scholar
 C. V. Le, H. NguyenXuan, H. Askes et al., “A cell‐based smoothed finite element method for kinematic limit analysis,” International Journal for Numerical Methods in Engineering, vol. 83, no. 12, pp. 1651–1674, 2010. View at: Google Scholar
 G. R. Liu, T. NguyenThoi, and K. Y. Lam, “An edgebased smoothed finite element method (ESFEM) for static, free and forced vibration analyses of solids,” Journal of Sound and Vibration, vol. 320, no. 45, pp. 1100–1130, 2009. View at: Publisher Site  Google Scholar
 N. NguyenThanh, T. Rabczuk, H. NguyenXuan, and S. P. A. Bordas, “An alternative alpha finite element method (AαFEM) for free and forced structural vibration using triangular meshes,” Journal of Computational and Applied Mathematics, vol. 233, no. 9, pp. 2112–2135, 2010. View at: Publisher Site  Google Scholar
 T. NguyenThoi, G. R. Liu, K. Y. Lam, and G. Y. Zhang, “A facebased smoothed finite element method (FSFEM) for 3D linear and geometrically nonlinear solid mechanics problems using 4node tetrahedral elements,” International Journal for Numerical Methods in Engineering, vol. 78, no. 3, pp. 324–353, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 H. NguyenXuan, T. Rabczuk, N. NguyenThanh, T. NguyenThoi, and S. Bordas, “A nodebased smoothed finite element method with stabilized discrete shear gap technique for analysis of ReissnerMindlin plates,” Computational Mechanics, vol. 46, no. 5, pp. 679–701, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. NguyenXuan, L. V. Tran, T. NguyenThoi, and H. C. VuDo, “Analysis of functionally graded plates using an edgebased smoothed finite element method,” Composite Structures, vol. 93, no. 11, pp. 3019–3039, 2011. View at: Publisher Site  Google Scholar
 G. R. Liu, H. NguyenXuan, T. NguyenThoi, and X. Xu, “A novel Galerkinlike weakform and a superconvergent alpha finite element method (SαFEM) for mechanics problems using triangular meshes,” Journal of Computational Physics, vol. 228, no. 11, pp. 4055–4087, 2009. View at: Publisher Site  Google Scholar
 H. NguyenXuan, L. V. Tran, C. H. Thai, and T. NguyenThoi, “Analysis of functionally graded plates by an efficient finite element method with nodebased strain smoothing,” ThinWalled Structures, vol. 54, pp. 1–18, 2012. View at: Publisher Site  Google Scholar
 H. NguyenXuan, T. Rabczuk, T. NguyenThoi, T. N. Tran, and N. NguyenThanh, “Computation of limit and shakedown loads using a nodebased smoothed finite element method,” International Journal for Numerical Methods in Engineering, vol. 90, no. 3, pp. 287–310, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. NguyenXuan, G. R. Liu, N. Nourbakhshnia, and L. Chen, “A novel singular ESFEM for crack growth simulation,” Engineering Fracture Mechanics, vol. 84, pp. 41–66, 2012. View at: Publisher Site  Google Scholar
 N. Nourbakhshnia and G. R. Liu, “A quasi‐static crack growth simulation based on the singular ES‐FEM,” International Journal for Numerical Methods in Engineering, vol. 88, no. 5, pp. 473–492, 2011. View at: Publisher Site  Google Scholar
 W. Zeng, G. R. Liu, Y. Kitamura, and H. NguyenXuan, “A threedimensional ESFEM for fracture mechanics problems in elastic solids,” Engineering Fracture Mechanics, vol. 114, pp. 127–150, 2013. View at: Publisher Site  Google Scholar
 P. PhungVan, T. NguyenThoi, T. LeDinh, and H. NguyenXuan, “Static and free vibration analyses and dynamic control of composite plates integrated with piezoelectric sensors and actuators by the cellbased smoothed discrete shear gap method (CSFEMDSG3),” Smart Materials and Structures, vol. 22, no. 9, Article ID 095026, 2013. View at: Publisher Site  Google Scholar
 P. PhungVan, L. De Lorenzis, C. H. Thai, M. AbdelWahab, and H. NguyenXuan, “Analysis of laminated composite plates integrated with piezoelectric sensors and actuators using higherorder shear deformation theory and isogeometric finite elements,” Computational Materials Science, vol. 96, pp. 495–505, 2015. View at: Publisher Site  Google Scholar
 S. P. A. Bordas, T. Rabczuk, N.X. Hung et al., “Strain smoothing in FEM and XFEM,” Computers & Structures, vol. 88, no. 2324, pp. 1419–1443, 2010. View at: Publisher Site  Google Scholar
 X. Zhao, S. P. A. Bordas, and J. Qu, “A hybrid smoothed extended finite element/level set method for modeling equilibrium shapes of nanoinhomogeneities,” Computational Mechanics, vol. 52, no. 6, pp. 1417–1428, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 L. Chen, T. Rabczuk, S. P. A. Bordas, G. R. Liu, K. Y. Zeng, and P. Kerfriden, “Extended finite element method with edgebased strain smoothing (ESmXFEM) for linear elastic crack growth,” Computer Methods in Applied Mechanics and Engineering, vol. 209–212, pp. 250–265, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Jiang, T. E. Tay, L. Chen, and X. S. Sun, “An edgebased smoothed XFEM for fracture in composite materials,” International Journal of Fracture, vol. 179, no. 12, pp. 179–199, 2013. View at: Publisher Site  Google Scholar
 N. VuBac, H. NguyenXuan, L. Chen et al., “A nodebased smoothed extended finite element method (NSXFEM) for fracture analysis,” Computer Modeling in Engineering and Sciences, vol. 73, no. 4, pp. 331–355, 2011. View at: Google Scholar  MathSciNet
 E. Bechet, M. Scherzer, and M. Kuna, “Application of the XFEM to the fracture of piezoelectric materials,” International Journal for Numerical Methods in Engineering, vol. 77, pp. 1535–1599, 2009. View at: Google Scholar
 U. Groh and M. Kuna, “Efficient boundary element analysis of cracks in 2D piezoelectric structures,” International Journal of Solids and Structures, vol. 42, no. 8, pp. 2399–2416, 2005. View at: Publisher Site  Google Scholar
 L. M. Zhou, G. W. Meng, F. Li, and H. Wang, “Cellbased smoothed finite element methodvirtual crack closure technique for a piezoelectric material of crack,” Mathematical Problems in Engineering, vol. 2015, Article ID 371083, 10 pages, 2015. View at: Publisher Site  Google Scholar
 P. Laborde, J. Pommier, Y. Renard, and M. Salaün, “Highorder extended finite element method for cracked domains,” International Journal for Numerical Methods in Engineering, vol. 64, no. 3, pp. 354–381, 2005. View at: Publisher Site  Google Scholar
 N. Moës, J. Dolbow, and T. Belytschko, “A finite element method for crack growth without remeshing,” International Journal for Numerical Methods in Engineering, vol. 46, no. 1, pp. 131–150, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. Rabczuk and T. Belytschko, “Cracking particles: a simplified meshfree method for arbitrary evolving cracks,” International Journal for Numerical Methods in Engineering, vol. 61, no. 13, pp. 2316–2343, 2004. View at: Publisher Site  Google Scholar
 T. Rabczuk, G. Zi, S. Bordas, and H. NguyenXuan, “A geometrically nonlinear threedimensional cohesive crack method for reinforced concrete structures,” Engineering Fracture Mechanics, vol. 75, no. 16, pp. 4740–4758, 2008. View at: Publisher Site  Google Scholar
 J. R. Rice, “Mathematical analysis in the mechanics of fracture,” in Fracture: An Advanced Treatise, H. Liebowitz, Ed., vol. 2, pp. 191–311, Academic Press, New York, NY, USA, 1968. View at: Google Scholar
 J. D. Eshelby, “Energy relations and the energymomentum tensor in continuum mechanics,” in Inelastic Behavior of Solids, M. F. Kanninen, W. F. Adler, A. R. Rosenfield, and R. I. Jaffee, Eds., pp. 77–115, McGrawHill, New York, NY, USA, 1970. View at: Google Scholar
Copyright
Copyright © 2016 Li Ming Zhou et al. 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.