Advances in Materials Science and Engineering

Advances in Materials Science and Engineering / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 2812748 | 16 pages | https://doi.org/10.1155/2019/2812748

A Coupling Electromechanical Inhomogeneous Cell-Based Smoothed Finite Element Method for Dynamic Analysis of Functionally Graded Piezoelectric Beams

Academic Editor: Ayman E. Hamada
Received16 Apr 2019
Accepted07 Jul 2019
Published06 Aug 2019

Abstract

To accurately simulate the continuous property change of functionally graded piezoelectric materials (FGPMs) and overcome the overstiffness of the finite element method (FEM), we present an electromechanical inhomogeneous cell-based smoothed FEM (ISFEM) of FGPMs. Firstly, ISFEM formulations were derived to calculate the transient response of FGPMs, and then, a modified Wilson-θ method was deduced to solve the integration of the FGPM system. The true parameters at the Gaussian integration point in FGPMs were adopted directly to replace the homogenization parameters in an element. ISFEM provides a close-to-exact stiffness of the continuous system, which could automatically and more easily generate for complicated domains and thus significantly decrease numerical errors. The accuracy and trustworthiness of ISFEM were verified as higher than the standard FEM by several numerical examples.

1. Introduction

Because of their outstanding electromechanical properties, easy fabricability, and preparation flexibility, piezoelectric materials are extensively applied as sensors and actuators to monitor and modulate the response of structures [1, 2]. Piezoelectric actuators and sensors are innovative for microscopic electromechanical systems and intelligent material systems, particularly in aerospace and medical fields [3]. Conventional piezoelectric sensors and actuators comprise multiple layers of various piezoelectric materials [47]. Moreover, piezoelectric layers with uniform material properties are limited by large bending displacement, stress concentration, creeping at high temperature, and failure from interfacial unbounding. All these phenomena are induced by mechanical or electric loading at layer interfaces [8].

To overcome the above limitations, Zhu et al. introduced and fabricated functionally graded piezoelectric material (FGPM) sensors and actuators [9, 10]. FGPMs change nonstop in one or more directions without generating internal stress concentration despite the production of large displacements. Takagi et al. fabricated FGPM bimorph actuators by using a mixed system of lead-zirconate-titanate (PZT) and Pt [11]. Nowadays, FGPMs are widely used intelligent materials for sensors and actuators in microstructural engineering. Many efforts have been made to analyze the behaviors and static/dynamic responses of FGPMs (e.g., shells, beams, and plates), such as the wave propagation study of FGPM plates based on the laminate theory [12]. Moreover, exact 3D analysis of FGPM rectangular plates was conducted by using a state-space approach [13] and to investigate the natural frequencies and mode shapes after being poled perpendicular to the middle plane [14]. The above method was also applied to explore the free vibration of rectangular FGPM plates [15]. The semianalytical finite element method (FEM) was used to investigate the static response of anisotropic and linear functionally graded magneto-electro-elastic plates [16]. Lezgy-Nazargah et al. [17] carried out static and dynamic analyses on piezoelectric beams by using a refined sinus model. The results have been found in good agreement with the reference solutions for various electrical and mechanical constrained conditions. Meanwhile, the 3D exact state-space solution [18] and Peano series solution [19] were developed for the cylindrical bending vibration of the FGPM laminates, respectively. Layerwise FEM was adopted to investigate the displacement and stress responses of an FGPM bimorph actuator [20]. Qiu et al. inhibited the vibration of a smart flexible clamped plate by using piezoelectric ceramic patch sensors and actuators [21]. The Timoshenko beam theory was used to analyze the static and dynamic responses of FGPM actuators to thermo-electro-mechanical loading [22]. The first-order shear deformation theory was used to study the static bending, free vibration, and dynamic responses of FGPM plates under electromechanical loading [23]. The free and induced vibrations of FGPM beams under thermo-electro-mechanical loading were characterized using the 3-order shear deformation beam theory [24]. A high-order theory for FGPM shells was proposed based on the generalized Hamilton’s principle [8]. Although the FEM (h-version) is adequate for low-frequency vibration analysis, it is not well suited to the vibration analysis of medium- or high-frequency regimes [25]. The spectral finite element method (SFEM) [26, 27] and the weak form quadrature method (QEM) [28, 29] are developed for the dynamic analysis of FGPM beams and structures.

Though FEM is the most widely used and effective numerical approach in practical issues in research and engineering (including mechanics of vibration), it is not necessarily fully perfect or cannot be further improved. For example, the probable overestimation of stiffness in solid structures may lead to locking behavior and inaccurate stress-solving [30]. By adding strain smoothing into FEM[31], Liu et al. established a series of cell-based [3236], node-based [37, 38], edge-based [3941], or face-based [42, 43] smoothed FEMs (S-FEMs) and their combinations [4447]. These S-FEMs with different properties can be used to get desired solutions for a variety of benchmarks and practical mechanic issues [4850]. The strain-smoothing operations can reduce or alleviate the overstiffness of standard FEM, significantly improving the accuracy of both primal and dual quantities [51]. Moreover, owing to absence of parametric mapping, the shape function derivatives and S-FEM models established in elasticity are not required to be insensitive to mesh distortion [52]. S-FEMs have been successfully extended to analyze the dynamic control of piezoelectric sensors and actuators, topological optimization of linear piezoelectric micromotors, statics, frequency, or defects of smart materials [5363]. Zheng et al. [64] utilized the cell-based smoothed finite element method with the asymptotic homogenization method to analyze the dynamic issues on micromechanics of piezoelectric composite materials. Zhou et al. [65, 66] deduced the linear and nonlinear cell-based smoothed finite element method of functionally graded magneto-electro-elastic (MEE) structures and further examined the transient responses of MEE sensors or energy harvest structures considering the damping factors. However, there is little literature reported concerning the dynamic response of FGPMs using the electromechanical inhomogeneous cell-based smoothed finite element method. Because of versatility, S-FEMs become convenient and efficient numerical approaches to address different physical issues.

Given the continuous change of the gradient of material properties along the thickness x3 direction and with cell-based gradient smoothing, we deduced the basic formula of ISFEM and a modified Wilson-θ method to solve the integral solution of the FGPM system. The displacements and potentials of FGPM cantilever beams under sine wave load, cosine wave load, step wave load, and triangular wave load were analyzed in comparison with FEM.

2. Basic Equations for Piezoelectric Materials

2.1. Geometry and Coordinate System

Each beam has a rectangular uniform cross section and is made of Ni layers either completely or partially composed of FGPM beams. The Cartesian coordinate system (x1, x2, x3) and geometric parameters are illustrated in Figure 1.

2.2. Constitutive Equations

At the k-th layer, 3D linear constitutive equations are polarized along its global coordinates as follows:where and are the stress tensor and infinitesimal strain tensor, respectively, and are electric field and electric displacement vector components, respectively; , , and are the piezoelectric, elastic, and dielectric material constants, respectively. Different from homogeneous piezoelectric materials, the three constants are dependent on coordinate . We assume that the material properties along the thickness direction are arbitrarily distributed as follows:where is an arbitrary function and , , and are values at the plane x3 = 0.

In an FGPM beam, equations (1) and (2) are reduced towhere

2.3. Weak Formulation

The principle of virtual work for a piezoelectric medium of volume Ω and regular boundary surface Γ can be written aswhere Fs, Fv, u, and φ are the vectors of surface force, mechanical body force, node displacements, and node electrical potentials, respectively; , , and are the electrical body charge, surface charge, and mass density, respectively; and is the virtual quantity.

3. Electromechanical ISFEM

The solving domain Ω is discretized into np elements, which contain Nn nodes; the approximation displacement and the approximation electrical potential for the FGPM problem can be expressed aswhere and are the ISFEM displacement shape function and electrical potential shape function, respectively.

Four-node element is divided into four smoothing subdomains. Field nodes, edge smoothing nodes, center smoothing nodes and edge Gaussian points, the outer normal vector distribution, and the shape function values are shown in Figure 2.

At any point in the smoothing subdomains , the smoothed strain and the smoothed electric field arewhere and are the strain and electric field in FEM, respectively, and is the constant function:where

Substituting equation (10) into equations (8) and (9), then we havewhere is the boundary of and and are the outer normal vector matrixes of the smoothing domain boundary

Equations (12) and (13) can be rewritten aswhere is the number of smoothing elements.

At the Gaussian point , equations (16) and (17) are where and are the Gaussian point and the length of the smoothing boundary, respectively, and is the total number of boundaries for each smoothing subdomain. As the shape function is linearly changed along each side of the smoothing subdomain, one Gauss point is sufficient for accurate boundary integration [30].

The essential distinction between ISFEM and FEM is that FEM needs to construct the shape function matrix of the element, while ISFEM only needs to use the shape function at the Gaussian point of the smoothing element boundary and does not require to involve the shape function derivatives. The above can reduce the continuity requirement of the shape function, and therefore, the accuracy and convergence of the method are improved.

The dynamic model of the FGPM electromechanical system can be derived from the Hamilton principle in the following form:wherewhere , (i = 1, 2, 3, 4) is the mass of the i-th smoothing element corresponding to node i, T is the smoothing element thickness, and is the density of Gaussian integration point of the i-th smoothing subdomain:where .

The application of the inhomogeneous smoothing element is to calculate stiffness matrix of the element. The parameters of four smoothing subdomains (i = 1, 2, 3, 4) are various in the elements, so the actual parameters at the Gaussian integration point are taken directly in order to reflect the changes of material property in each element.

4. Modified Wilson-θ Method

The modified Wilson-θ method is an important scheme and an implicit integral way to solve the dynamic system equations [63]. If θ > 1.37, the solution is unconditionally stable. The detailed procedures are showed as follows:

4.1. Initial Calculation

(1)Formulate generalized stiffness matrix , mass matrix , and damping matrix (2)Calculate initial values of , , (3)Select the time step Δt and the integral constant θ (θ = 1.4)(4)Formulate an effective generalized stiffness matrix :

4.2. For Each Time Step

(1)Calculate the payload at time t + θΔt:(2)Calculate the generalized displacement at time t + θΔt:(3)Calculate the generalized acceleration, generalized speed, and generalized displacement at time t + Δt:

5. Numerical Examples

Four numerical examples were conducted under sine wave load, cosine wave load, step wave load and triangular wave load, respectively. FGPM cantilever beams of the same dimensions (length L = 40 mm, width h = 5 mm and thickness b = 1 mm) were subjected to forced vibration (Figure 3). The material constants are shown in Table 1. And initial conditions were and at t = 0 moment. The FGPM beams were made of PZT-4 or PZT-5H on basis of exponentially graded piezoelectric materials with the following material properties:where and n is the gradient parameter.


Materialc11 (GPa)c12 (GPa)c13 (GPa)c33 (GPa)c44 (GPa) (kg/m³)e31 (C/m2)e33 (C/m2)e15 (C/m2)χ11/χ0χ33/χ0χ0 (C/N·m2)

PZT-4139.077.8074.30115.025.607600−5.2015.1012.701476.01301.08.854 × 10−12
PZT-5H126.079.6084.10114.023.307500−6.5017.4423.301468.91698.38.854 × 10−12

5.1. Sine Wave Load

The load F applied to the free end and the load waveform is demonstrated in Figure 4. A convergence investigation with respect to meshes was first carried out. Four smoothing subcells were used for electromechanical ISFEM with ∆t = 1 × 10−3 s. The variations of displacement u3 and electrical potential φ at the loading point of the PZT-4-based FGPM beam combined with respect to time are shown in Figure 5. The results at n = −5 with the element number of 480, 800, 1200, or 1680 were compared with the reference solution [67]. The variations of u3 and φ at the loading point combined with respect to time at n = −1, 0, 1, and 5 in comparison with the reference solution are shown in Figure 6 [67]. Figure 7 illustrates the total energy norm Err versus the mesh density at t = 0.002 s and t = 0.01 s. The simulation results are well consistent among different numbers of meshes, which demonstrate the high convergence of ISFEM.

Figure 8 shows a comparison of calculation time between ISFEM and FEM at Intel ® Xeon ® CPU E3-1220 v3 @ 3.10 GHz, 16GB RAM. The bandwidth of the system matrices for the FEM and ISFEM is identical. However, in the ISFEM, only the values of shape functions (not the derivatives) at the quadrature points are needed and the requirement of traditional coordinate transform procedure is not necessary to perform the numerical integration. Therefore, the ISFEM generally needs less computational cost than the FEM for handling the dynamic analysis problems.

The variations of u3 and φ at the loading point combined with t = 0.0001 s, 0.0004 s, 0.0025 s, and 0.004 s in the PZT-4-based FGPM cantilever beam are shown in Table 2. The 80 × 10 meshes of ISFEM at n = −5, −1, 0, 1, and 5 are shown in Figure 9, and FEM is considered 160 × 20 elements. Clearly, the results of ISFEM with 80 × 10 elements are the same as the calculated results of FEM using 160 × 20 elements, suggesting ISFEM has higher accuracy.


Time (s)Field variablesn = −5n = −1n = 0n = 1n = 5
ISFEMFEMISFEMFEMISFEMFEMISFEMFEMISFEMFEM

0.0001u3 × 10−10 m12.85212.6582.7602.7511.7071.7021.0061.0040.07560.0754
φ × 10−2 V11.84811.9120.4420.4270.2120.2070.09610.09440.001120.00110

0.0004u3 × 10−9 m28.51728.3674.7354.7152.8552.8421.7381.7310.1870.187
φ × 10−1 V5.8515.7470.5970.5820.3490.3390.1930.1880.08050.782

0.0025u3 × 10−8 m13.34513.1911.5381.5280.9170.9130.5640.5610.08820.0873
φ (V)2.4472.3910.1980.1910.1120.1090.06210.06050.003860.00372

0.004U3 × 10−9 m77.79977.68210.15410.1466.1656.1233.7273.7250.5140.516
φ (V)1.4341.4060.1290.1250.07560.07400.04140.04010.002240.00221

The variations of u3 and φ at the loading point combined with time of the PZT-5H-based FGPM cantilever beam are listed in Table 3. Clearly, when n changes from −5 to 5, the maximum u3 and φ decrease, which is consistent with the PZT-4-based FGPM cantilever beam. Furthermore, the results of ISFEM with 80 × 10 elements are the same as the calculated results of FEM using 160 × 20 elements, suggesting ISFEM has higher accuracy.


Time (s)Field variablesn = −5n = −1n = 0n = 1n = 5
ISFEMFEMISFEMFEMISFEMFEMISFEMFEMISFEMFEM

0.0001u3 × 10−10 m13.69113.4242.9512.9421.8261.8211.0761.0740.08020.0800
φ × 10−2 V12.09812.4930.4130.4030.1900.1940.08250.08430.0006240.000649

0.0004u3 × 10−9 m31.39431.1965.5835.5613.3813.3672.0502.0420.2070.206
φ × 10−1 V5.9016.0570.5990.6070.3490.3500.1920.1900.07360.0735

0.0025u3 × 10−8 m16.12216.0062.0172.0031.1181.1170.7410.7360.1070.106
φ (V)2.5222.5840.2190.2190.1220.1230.06910.06890.003920.00393

0.004u3 × 10−9 m76.01975.05210.33910.2646.2646.2613.7953.7680.5010.496
φ (V)1.4021.4370.1130.1130.06490.06530.03510.03500.001790.00178

5.2. Cosine Wave Load

The cosine load F applied to the free end and load waveform is shown in Figure 10. The variations of u3 and φ at the loading point combined with time of PZT-4- and PZT-5H-based FGPM cantilever beams are listed in Tables 4 and 5, respectively. The calculated results of ISFEM with 80 × 10 elements are the same as those of FEM using 160 × 20 elements, implying that ISFEM possesses higher accuracy.


Time (s)Field variablesn = −5n = −1n = 0n = 1n = 5
ISFEMFEMISFEMFEMISFEMFEMISFEMFEMISFEMFEM

0.0001u3 × 10−9 m20.42720.1594.3874.3822.7132.7101.6001.5990.1200.120
φ × 10−2 V188.326189.8217.0296.8963.3693.2691.5281.50440.01780.0175

0.0004u3 × 10−8 m21.20821.2342.7382.7431.6101.6111.0061.0040.1410.142
φ × 10−1 V27.88927.9493.3153.2160.1930.1871.1131.0980.06480.0639

0.0025u3 × 10−9 m−33.487−32.845−12.058−11.952−4.445−4.345−4.436−4.392−0.225−0.221
φ (V)−0.138−0.137−0.132−0.128−0.0491−0.0479−0.0486−0.0470−0.00086−0.00085

0.004u3 × 10−9 m−11.111−11.501−0.984−0.969−5.549−5.616−0.351−0.339−0.0981−0.0913
φ (V)−1.408−1.365−0.0176−0.0174−0.0653−0.06350.002340.002290.0005960.000579


Time (s)Field variablesn = −5n = −1n = 0n = 1n = 5
ISFEMFEMISFEMFEMISFEMFEMISFEMFEMISFEMFEM

0.0001u3 × 10−9 m2.9022.9014.6914.6852.9022.9011.7101.7100.1280.127
φ × 10−2 V192.287198.9676.5616.7423.0173.0921.3121.3420.009920.0101
0.0004u3 × 10−8 m23.94123.9483.4863.4932.0692.0641.2811.2840.1600.159
φ × 10−1 V28.17829.0943.58413.6132.1082.1211.2031.2010.06210.0616
0.0025u3 × 10−9 m117.958122.357−4.831−4.723−7.090−6.942−1.777−1.7150.7950.8136
φ (V)0.508630.5017−0.0407−0.0393−0.0688−0.0652−0.0151−0.01480.003380.00342
0.004u3 × 10−9 m−183.26−179.45−27.847−27.063−9.636−9.387−10.233−9.998−1.287−1.263
φ (V)−2.273−2.176−0.284−0.277−0.0931−0.0885−0.0951−0.0949−0.00466−0.00451

5.3. Step Wave Load

The step load F applied to the free end and load waveform is indicated in Figure 11. The variations of u3 and φ at the loading point combined with time of PZT-4- and PZT-5H-based FGPM cantilever beams are illustrated in Figures 12 and 13, respectively. It is clearly shown that ISFEM possesses higher accuracy than FEM for the calculated results of ISFEM using 80 × 10 elements are the same as FEM using 160 × 20 elements.

5.4. Triangular Wave Load

The triangular load F applied to the free end and load waveform is shown in Figure 14. The variations of u3 and φ at the loading point combined with time of PZT-4- and PZT-5H-based FGPM cantilever beams are shown in Figures 15 and 16, respectively. It shows that the solutions of CS-FEM with less elements are the same as the solutions of FEM using more elements.

6. Conclusions

An electromechanical ISFEM was proposed given the continuous changes of the gradient of material properties along the thickness x3 direction and with cell-based gradient smoothing. The modified Wilson-θ method was deduced to solve the integral solution of the FGPM system. The displacements and potentials of cantilever beams combining with sine load, cosine load, step load, and triangular load were analyzed by ISFEM in comparison with FEM.(1)ISFEM is correct and effective in solving the dynamic response of FGPM structures(2)ISFEM can reduce the systematic stiffness of FEM and provides calculations closer to the true values(3)ISFEM is more efficient than FEM and takes less computation time at the same accuracy

This study indicates a possibility to select suitable grading controlled by the power law index according to the application.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Liming Zhou and Bin Cai performed the simulations. Bin Cai contributed to the writing of the manuscript.

Acknowledgments

Special thanks are due to Professor Guirong Liu for the S-FEM Source Code in http://www.ase.uc.edu/∼liugr/software.html. This work was financially supported by the National Key R&D Program of China (grant no. 2018YFF01012401-06), Jilin Provincial Department of Science and Technology Fund Project (grant no. 20170101043JC), Jilin Provincial Department of Education (grant nos. JJKH20180084KJ and JJKH20170788KJ), Fundamental Research Funds for the Central Universities, Science and Technological Planning Project of Ministry of Housing and Urban–Rural Development of the People's Republic of China (2017-K9-047) and Graduate Innovation Fund of Jilin University (grant no. 101832018C184).

References

  1. K. M. Liew, X. Q. He, and S. Kitipornchai, “Finite element method for the feedback control of FGM shells in the frequency domain via piezoelectric sensors and actuators,” Computer Methods in Applied Mechanics and Engineering, vol. 193, no. 3–5, pp. 257–273, 2004. View at: Publisher Site | Google Scholar
  2. G. Song, V. Sethi, and H.-N. Li, “Vibration control of civil structures using piezoceramic smart materials: a review,” Engineering Structures, vol. 28, no. 11, pp. 1513–1524, 2006. View at: Publisher Site | Google Scholar
  3. B. Legrand, J.-P. Salvetat, B. Walter, M. Faucher, D. Théron, and J.-P. Aimé, “Multi-MHz micro-electro-mechanical sensors for atomic force microscopy,” Ultramicroscopy, vol. 175, pp. 46–57, 2017. View at: Publisher Site | Google Scholar
  4. F. Pablo, I. Bruant, and O. Polit, “Use of classical plate finite elements for the analysis of electroactive composite plates. Numerical validations,” Journal of Intelligent Material Systems and Structures, vol. 20, no. 15, pp. 1861–1873, 2009. View at: Publisher Site | Google Scholar
  5. M. D’Ottavio and O. Polit, “Sensitivity analysis of thickness assumptions for piezoelectric plate models,” Journal of Intelligent Material Systems and Structures, vol. 20, no. 15, pp. 1815–1834, 2009. View at: Publisher Site | Google Scholar
  6. P. Vidal, M. D’Ottavio, M. Ben Thaïer, and O. Polit, “An efficient finite shell element for the static response of piezoelectric laminates,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 7, pp. 671–690, 2011. View at: Publisher Site | Google Scholar
  7. S. B. Beheshti-Aval, M. Lezgy-Nazargah, P. Vidal, and O. Polit, “A refined sinus finite element model for the analysis of piezoelectric-laminated beams,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 3, pp. 203–219, 2011. View at: Publisher Site | Google Scholar
  8. X.-H. Wu, C. Chen, Y.-P. Shen, and X.-G. Tian, “A high order theory for functionally graded piezoelectric shells,” International Journal of Solids and Structures, vol. 39, no. 20, pp. 5325–5344, 2002. View at: Publisher Site | Google Scholar
  9. X. H. Zhu and Z. Y. Meng, “Operational principle, fabrication and displacement characteristic of a functionally gradient piezoelectric ceramic actuator,” Sensors and Actuators A. Physical, vol. 48, no. 3, pp. 169–176, 1995. View at: Publisher Site | Google Scholar
  10. C. C. M. Wu, M. Kahn, and W. Moy, “Piezoelectric ceramics with functional gradients: a new application in material design,” Journal of the American Ceramic Society, vol. 79, no. 3, pp. 809–812, 2005. View at: Publisher Site | Google Scholar
  11. K. Takagi, J.-F. Li, S. Yokoyama, R. Watanabe, A. Almajid, and M. Taya, “Design and fabrication of functionally graded PZT/Pt piezoelectric bimorph actuator,” Science and Technology of Advanced Materials, vol. 3, no. 2, pp. 217–224, 2002. View at: Publisher Site | Google Scholar
  12. G. R. Liu and J. Tani, “Surface waves in functionally gradient piezoelectric plates,” Journal of Vibration and Acoustics, vol. 116, no. 4, pp. 440–448, 1994. View at: Publisher Site | Google Scholar
  13. Z. Zhong and E. T. Shang, “Three-dimensional exact analysis of a simply supported functionally gradient piezoelectric plate,” International Journal of Solids and Structures, vol. 40, no. 20, pp. 5335–5352, 2003. View at: Publisher Site | Google Scholar
  14. C. Piotr, “Three-dimensional natural vibration analysis and energy considerations for a piezoelectric rectangular plate,” Journal of Sound and Vibration, vol. 283, no. 3–5, pp. 1093–1113, 2005. View at: Publisher Site | Google Scholar
  15. W. Q. Chen and H. J. Ding, “On free vibration of a functionally graded piezoelectric rectangular plate,” Acta Mechanica, vol. 153, no. 3-4, pp. 207–216, 2002. View at: Publisher Site | Google Scholar
  16. R. K. Bhangale and N. Ganesan, “Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates,” International Journal of Solids and Structures, vol. 43, no. 10, pp. 3230–3253, 2006. View at: Publisher Site | Google Scholar
  17. M. Lezgy-Nazargah, P. Vidal, and O. Polit, “An efficient finite element model for static and dynamic analyses of functionally graded piezoelectric beams,” Composite Structures, vol. 104, pp. 71–84, 2013. View at: Publisher Site | Google Scholar
  18. M. Lezgy-Nazargah, “A three-dimensional exact state-space solution for cylindrical bending of continuously non-homogenous piezoelectric laminated plates with arbitrary gradient composition,” Archive of Mechanics, vol. 67, pp. 25–51, 2015. View at: Google Scholar
  19. M. Lezgy-Nazargah, “A three-dimensional Peano series solution for the vibration of functionally graded piezoelectric laminates in cylindrical bending,” Scientia Iranica, vol. 23, no. 3, pp. 788–801, 2016. View at: Publisher Site | Google Scholar
  20. H.-J. Lee, “Layerwise laminate analysis of functionally graded piezoelectric bimorph beams,” Journal of Intelligent Material Systems and Structures, vol. 16, no. 4, pp. 365–371, 2005. View at: Publisher Site | Google Scholar
  21. Z.-C. Qiu, X.-M. Zhang, H.-X. Wu, and H.-H. Zhang, “Optimal placement and active vibration control for piezoelectric smart flexible cantilever plate,” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 521–543, 2007. View at: Publisher Site | Google Scholar
  22. J. Yang and H. J. Xiang, “Thermo-electro-mechanical characteristics of functionally graded piezoelectric actuators,” Smart Materials and Structures, vol. 16, no. 3, pp. 784–797, 2007. View at: Publisher Site | Google Scholar
  23. B. Behjat, M. Salehi, A. Armin, M. Sadighi, and M. Abbasi, “Static and dynamic analysis of functionally graded piezoelectric plates under mechanical and electrical loading,” Scientia Iranica, vol. 18, no. 4, pp. 986–994, 2011. View at: Publisher Site | Google Scholar
  24. A. Doroushi, M. R. Eslami, and A. Komeili, “Vibration analysis and transient response of an FGPM beam under thermo-electro-mechanical loads using higher-order shear deformation theory,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 3, pp. 231–243, 2011. View at: Publisher Site | Google Scholar
  25. G. W. Wei, Y. B. Zhao, and Y. Xiang, “A novel approach for the analysis of high-frequency vibrations,” Journal of Sound and Vibration, vol. 257, no. 2, pp. 207–246, 2002. View at: Publisher Site | Google Scholar
  26. P. Kudela, M. Krawczuk, and W. Ostachowicz, “Wave propagation modelling in 1D structures using spectral finite elements,” Journal of Sound and Vibration, vol. 300, no. 1-2, pp. 88–100, 2007. View at: Publisher Site | Google Scholar
  27. Y. Kim, S. Ha, and F.-K. Chang, “Time-domain spectral element method for built-in piezoelectric-actuator-induced lamb wave propagation analysis,” AIAA Journal, vol. 46, no. 3, pp. 591–600, 2008. View at: Publisher Site | Google Scholar
  28. H. Zhong and Z. Yue, “Analysis of thin plates by the weak form quadrature element method,” Science China Physics, Mechanics and Astronomy, vol. 55, no. 5, pp. 861–871, 2012. View at: Publisher Site | Google Scholar
  29. X. Wang, Z. Yuan, and C. Jin, “Weak form quadrature element method and its applications in science and engineering: a state-of-the-art review,” Applied Mechanics Reviews, vol. 69, no. 3, Article ID 030801, 2017. View at: Publisher Site | Google Scholar
  30. G. R. Liu and T. Nguyen-Thoi, Smoothed Finite Element Methods, CRC Press, Taylor and Francis Group, Boca Raton, FL, USA, 2010.
  31. J.-S. Chen, C.-T. Wu, S. Yoon, and Y. You, “A stabilized conforming nodal integration for Galerkin mesh-free methods,” International Journal for Numerical Methods in Engineering, vol. 50, no. 2, pp. 435–466, 2001. View at: Publisher Site | Google Scholar
  32. K. Y. Dai and G. R. Liu, “Free and forced vibration analysis using the smoothed finite element method (SFEM),” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 803–820, 2007. View at: Publisher Site | Google Scholar
  33. S. P. A. Bordas and S. Natarajan, “On the approximation in the smoothed finite element method (SFEM),” International Journal for Numerical Methods in Engineering, vol. 81, no. 5, pp. 660–670, 2010. View at: Publisher Site | Google Scholar
  34. C. V. Le, H. Nguyen-Xuan, H. Askes, S. P. A. Bordas, T. Rabczuk, and H. Nguyen-Vinh, “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: Publisher Site | Google Scholar
  35. G. R. Liu, W. Zeng, and H. Nguyen-Xuan, “Generalized stochastic cell-based smoothed finite element method (GS_CS-FEM) for solid mechanics,” Finite Elements in Analysis and Design, vol. 63, pp. 51–61, 2013. View at: Publisher Site | Google Scholar
  36. K. Nguyen-Quang, H. Dang-Trung, V. Ho-Huu, H. Luong-Van, and T. Nguyen-Thoi, “Analysis and control of FGM plates integrated with piezoelectric sensors and actuators using cell-based smoothed discrete shear gap method (CS-DSG3),” Composite Structures, vol. 165, pp. 115–129, 2017. View at: Publisher Site | Google Scholar
  37. Y. H. Bie, X. Y. Cui, and Z. C. Li, “A coupling approach of state-based peridynamics with node-based smoothed finite element method,” Computer Methods in Applied Mechanics and Engineering, vol. 331, pp. 675–700, 2018. View at: Publisher Site | Google Scholar
  38. G. Liu, M. Chen, and M. Li, “Lower bound of vibration modes using the node-based smoothed finite element method (NS-FEM),” International Journal of Computational Methods, vol. 14, no. 4, Article ID 1750036, 2016. View at: Publisher Site | Google Scholar
  39. Z. C. He, G. Y. Li, Z. H. Zhong et al., “An ES-FEM for accurate analysis of 3D mid-frequency acoustics using tetrahedron mesh,” Computers & Structures, vol. 106-107, pp. 125–134, 2012. View at: Publisher Site | Google Scholar
  40. X. Y. Cui, G. Wang, and G. Y. Li, “A nodal integration axisymmetric thin shell model using linear interpolation,” Applied Mathematical Modelling, vol. 40, no. 4, pp. 2720–2742, 2016. View at: Publisher Site | Google Scholar
  41. X. Y. Cui, X. B. Hu, and Y. Zeng, “A copula-based perturbation stochastic method for fiber-reinforced composite structures with correlations,” Computer Methods in Applied Mechanics and Engineering, vol. 322, pp. 351–372, 2017. View at: Publisher Site | Google Scholar
  42. T. Nguyen-Thoi, G. R. Liu, K. Y. Lam, and G. Y. Zhang, “A face-based smoothed finite element method (FS-FEM) for 3D linear and geometrically non-linear solid mechanics problems using 4-node tetrahedral elements,” International Journal for Numerical Methods in Engineering, vol. 78, no. 3, pp. 324–353, 2009. View at: Publisher Site | Google Scholar
  43. Z. C. He, G. Y. Li, Z. H. Zhong et al., “An edge-based smoothed tetrahedron finite element method (ES-T-FEM) for 3D static and dynamic problems,” Computational Mechanics, vol. 52, no. 1, pp. 221–236, 2013. View at: Publisher Site | Google Scholar
  44. E. Li, Z. C. He, X. Xu, and G. R. Liu, “Hybrid smoothed finite element method for acoustic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 283, pp. 664–688, 2015. View at: Publisher Site | Google Scholar
  45. C. Jiang, Z.-Q. Zhang, G. R. Liu, X. Han, and W. Zeng, “An edge-based/node-based selective smoothed finite element method using tetrahedrons for cardiovascular tissues,” Engineering Analysis with Boundary Elements, vol. 59, pp. 62–77, 2015. View at: Publisher Site | Google Scholar
  46. X. B. Hu, X. Y. Cui, Z. M. Liang, and G. Y. Li, “The performance prediction and optimization of the fiber-reinforced composite structure with uncertain parameters,” Composite Structures, vol. 164, pp. 207–218, 2017. View at: Publisher Site | Google Scholar
  47. X. Cui, S. Li, H. Feng, and G. Li, “A triangular prism solid and shell interactive mapping element for electromagnetic sheet metal forming process,” Journal of Computational Physics, vol. 336, pp. 192–211, 2017. View at: Publisher Site | Google Scholar
  48. G. R. Liu, H. Nguyen-Xuan, and T. Nguyen-Thoi, “A theoretical study on the smoothed FEM (S-FEM) models: properties, accuracy and convergence rates,” International Journal for Numerical Methods in Engineering, vol. 84, no. 10, pp. 1222–1256, 2010. View at: Publisher Site | Google Scholar
  49. E. Li, Z. C. He, G. Wang, and G. R. Liu, “An efficient algorithm to analyze wave propagation in fluid/solid and solid/fluid phononic crystals,” Computer Methods in Applied Mechanics and Engineering, vol. 333, pp. 421–442, 2018. View at: Publisher Site | Google Scholar
  50. W. Zeng, G. R. Liu, D. Li, and X. W. Dong, “A smoothing technique based beta finite element method (βFEM) for crystal plasticity modeling,” Computers & Structures, vol. 162, pp. 48–67, 2016. View at: Publisher Site | Google Scholar
  51. H. Nguyen-Xuan, G. R. Liu, T. Nguyen-Thoi, and C. Nguyen-Tran, “An edge-based smoothed finite element method for analysis of two-dimensional piezoelectric structures,” Smart Materials and Structures, vol. 18, no. 6, Article ID 065015, 2009. View at: Publisher Site | Google Scholar
  52. H. Nguyen-Van, N. Mai-Duy, and T. Tran-Cong, “A smoothed four-node piezoelectric element for analysis of two-dimensional smart structures,” Computer Modeling in Engineering and Sciences, vol. 23, no. 3, pp. 209–222, 2008. View at: Publisher Site | Google Scholar
  53. P. Phung-Van, T. Nguyen-Thoi, T. Le-Dinh, and H. Nguyen-Xuan, “Static and free vibration analyses and dynamic control of composite plates integrated with piezoelectric sensors and actuators by the cell-based smoothed discrete shear gap method (CS-FEM-DSG3),” Smart Materials and Structures, vol. 22, no. 9, Article ID 095026, 2013. View at: Publisher Site | Google Scholar
  54. Z. C. He, E. Li, G. R. Liu, G. Y. Li, and A. G. Cheng, “A mass-redistributed finite element method (MR-FEM) for acoustic problems using triangular mesh,” Journal of Computational Physics, vol. 323, pp. 149–170, 2016. View at: Publisher Site | Google Scholar
  55. W. Zuo and K. Saitou, “Multi-material topology optimization using ordered SIMP interpolation,” Structural and Multidisciplinary Optimization, vol. 55, no. 2, pp. 477–491, 2017. View at: Publisher Site | Google Scholar
  56. E. Li, Z. C. He, J. Y. Hu, and X. Y. Long, “Volumetric locking issue with uncertainty in the design of locally resonant acoustic metamaterials,” Computer Methods in Applied Mechanics and Engineering, vol. 324, pp. 128–148, 2017. View at: Publisher Site | Google Scholar
  57. L. Chen, Y. W. Zhang, G. R. Liu, H. Nguyen-Xuan, and Z. Q. Zhang, “A stabilized finite element method for certified solution with bounds in static and frequency analyses of piezoelectric structures,” Computer Methods in Applied Mechanics and Engineering, vol. 241–244, pp. 65–81, 2012. View at: Publisher Site | Google Scholar
  58. L. Zhou, M. Li, Z. Ma et al., “Steady-state characteristics of the coupled magneto-electro-thermo-elastic multi-physical system based on cell-based smoothed finite element method,” Composite Structures, vol. 219, pp. 111–128, 2019. View at: Publisher Site | Google Scholar
  59. L. Zhou, S. Ren, G. Meng, X. Li, and F. Cheng, “A multi-physics node-based smoothed radial point interpolation method for transient responses of magneto-electro-elastic structures,” Engineering Analysis with Boundary Elements, vol. 101, pp. 371–384, 2019. View at: Publisher Site | Google Scholar
  60. H. Nguyen-Van, N. Mai-Duy, and T. Tran-Cong, “A node-based element for analysis of planar piezoelectric structures,” CMES: Computer Modeling in Engineering & Sciences, vol. 36, no. 1, pp. 65–96, 2008. View at: Google Scholar
  61. E. Li, Z. C. He, Y. Jiang, and B. Li, “3D mass-redistributed finite element method in structural-acoustic interaction problems,” Acta Mechanica, vol. 227, no. 3, pp. 857–879, 2016. View at: Publisher Site | Google Scholar
  62. L. Zhou, M. Li, G. Meng, and H. Zhao, “An effective cell-based smoothed finite element model for the transient responses of magneto-electro-elastic structures,” Journal of Intelligent Material Systems and Structures, vol. 29, no. 14, pp. 3006–3022, 2018. View at: Publisher Site | Google Scholar
  63. L. Zhou, M. Li, H. Zhao, and W. Tian, “Cell-based smoothed finite element method for the intensity factors of piezoelectric bimaterials with interfacial crack,” International Journal of Computational Methods, vol. 16, no. 7, Article ID 1850107, 2019. View at: Google Scholar
  64. J. Zheng, Z. Duan, and L. Zhou, “A coupling electromechanical cell-based smoothed finite element method based on micromechanics for dynamic characteristics of piezoelectric composite materials,” Advances in Materials Science and Engineering, vol. 2019, Article ID 4913784, 16 pages, 2019. View at: Publisher Site | Google Scholar
  65. L. Zhou, S. Ren, C. Liu, and Z. Ma, “A valid inhomogeneous cell-based smoothed finite element model for the transient characteristics of functionally graded magneto-electro-elastic structures,” Composite Structures, vol. 208, pp. 298–313, 2019. View at: Publisher Site | Google Scholar
  66. L. Zhou, M. Li, B. Chen, F. Li, and X. Li, “An inhomogeneous cell-based smoothed finite element method for the nonlinear transient response of functionally graded magneto-electro-elastic structures with damping factors,” Journal of Intelligent Material Systems and Structures, vol. 30, no. 3, pp. 416–437, 2019. View at: Publisher Site | Google Scholar
  67. R. X. Yao and Z. F. Shi, “Steady-state forced vibration of functionally graded piezoelectric beams,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 8, pp. 769–779, 2011. View at: Publisher Site | Google Scholar

Copyright © 2019 Bin Cai and Liming Zhou. 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.


More related articles

521 Views | 211 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.