Shock and Vibration
Volume 2018 (2018), Article ID 5141060, 17 pages
https://doi.org/10.1155/2018/5141060

Research Article

An Inhomogeneous Cell-Based Smoothed Finite Element Method for Free Vibration Calculation of Functionally Graded Magnetoelectroelastic Structures

Yan Cai, Guangwei Meng, and Liming Zhou

School of Mechanical Science and Engineering, Jilin University, Changchun 130025, China

Correspondence should be addressed to Liming Zhou; lmzhou@jlu.edu.cn

Received 27 October 2017; Revised 23 December 2017; Accepted 15 January 2018; Published 7 March 2018

Academic Editor: Marco Alfano

Copyright © 2018 Yan Cai 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.

Abstract.

To overcome the overstiffness and imprecise magnetoelectroelastic coupling effects of finite element method (FEM), we present an inhomogeneous cell-based smoothed FEM (ICS-FEM) of functionally graded magnetoelectroelastic (FGMEE) structures. Then the ICS-FEM formulations for free vibration calculation of FGMEE structures were deduced. In FGMEE structures, the true parameters at the Gaussian integration point were adopted directly to replace the homogenization in an element. The ICS-FEM provides a continuous system with a close-to-exact stiffness, which could be automatically and more easily generated for complicated domains, thus significantly decreasing the numerical error. To verify the accuracy and trustworthiness of ICS-FEM, we investigated several numerical examples and found that ICS-FEM simulated more accurately than the standard FEM. Also the effects of various equivalent stiffness matrices and the gradient function on the inherent frequency of FGMEE beams were studied.

1. Introduction

Functionally graded magnetoelectroelastic (FGMEE) materials are generally multiphase composites with continuously varying mechanical properties. FGMEE materials can convert magnetic, electric, and mechanical energy from one type into another and have received wide attention recently due to their electroelastic, magnetoelastic, and electromagnetic coupling effects [1, 2]. Therefore, FGMEE materials have been adopted in various smart structures, such as magnetic field probes, smart vibration sensors, optoelectronic devices, and medical ultrasonic transducers [3, 4]. The smart FGMEE structures are commonly fabricated in the beam pattern. However, to better apply FGMEE beams, researchers must analyse the statics and free vibration, and to predict the coupled response of FGMEE beams for practical applications, they should accurately calculate the properties of free vibrations.

Several computational techniques were proposed to investigate the electroelastic, magnetoelastic, and electromagnetic coupling effects of smart structures, such as finite element method (FEM), mesh-free method, and scaled boundary FEM [5–10]. Bhangale and Ganesan analyzed the static behaviors of linear anisotropic FGMEE plates using semianalytical FEM and investigated the free vibration of FGMEE plates and cylindrical shells [11, 12]. A layerwise partial mixed FEM was proposed to model MEE plates [13]. Phoenix et al. analysed the static and dynamic behaviors of coupled MEE plates using FEM with the Reissner mixed variational theorem [14]. Buchanan used FEM to study the free vibrations of infinite magnetoelectroelastic cylinders [15]. However, these FEMs overestimated the stiffness of solid structures and were limited by low accuracy. Therefore, Sladek et al. proposed a mesh-free method to more accurately study the static behavior of a circular FGMEE plate [16], but this method reduced the computational effectiveness. Recently, Liu et al. solved the deformations of a nonuniform MEE plate using scaled boundary FEM [17]. By incorporating the nonlocal theory into scaled boundary FEM, Ke and Wang more accurately and effectively studied the free vibrations of MEE beams [18]. However, the effectiveness of scaled boundary FEM is still low and should be improved.

FEM, as a powerful computational tool to investigate MEE coupling behaviors, yet overestimates the stiffness of FGMEE structures, which may result in locking behavior and inaccurate eigenvalue solutions [19]. To overcome these limitations, a series of cell-based smoothed FEM (CS-FEM) [20–26] and node-based, edge-based, or face-based smoothed FEMs [27–32] were proposed. In recent years, many CS-FEM-based formulations were proposed. Moreover, CS-FEM does not require the shape function derivatives or high generosity of program and is insensitive to mesh distortion because of the absence of isoparametric mapping.

CS-FEM has been successfully extended into dynamical control of piezoelectric sensors and actuators, topological optimization of linear piezoelectric micromotor, and analysis of static behaviors, frequency, and defects of piezoelectric structures [33–43]. Due to its versatility, CS-FEM becomes a simple and effective numerical tool to solve numerous electric and mechanical physical problems. However, the application of CS-FEM to investigate MEE properties is still a challenge.

In this work, free vibrations of FGMEE structures were studied. Inhomogeneous CS-FEM (ICS-FEM) for FGMEE materials was formulated by incorporating gradient smoothing into the standard FEM for the multi-physics field of FGMEE. Then the equations of free vibration computation were deduced under the multi-physics coupling field for FGMEE materials. Finally, FGMEE beams in the functional gradient exponential form or power law form were calculated under different boundary conditions. ICS-FEM outperformed FEM when compared with the reference solution.

This paper is organized as follows: Section 2 introduces the basic formulations for FGMEE materials. Section 3 briefly describes the properties of the FGMEE materials. Section 4 contains the detailed formulation of ICS-FEM. In Section 5, two numerical examples and the model of typical MEMS-based FGMEE energy harvester are investigated in detail. Final conclusions from the numerical results are drawn in Section 6.

2. Basic Formulations

The material properties of a functionally graded material (FGM) plate vary continuously, which is an advantage over the discontinuity across adjoining layers in a laminated plate. The wide range of engineering applications of FGM has attracted many scientists to investigate the behaviors of FGM.

Considering the transverse isotropy of the FGMEE medium [9, 44] and for the plane stress problem, we set stress components = = = 0, electric displacement component = 0, and magnetic induction component = 0. The geometric parameters and the chosen Cartesian coordinate system are illustrated in Figure 1.

Figure 1: FGMEE beams: Cartesian coordinate system and geometric parameters.

The basic formulations for MEE materials include equilibrium equations, geometric equations, and constitutive equations. The equilibrium equations are

where , and denote stress components; and are electric displacement components; and are magnetic induction components.

The geometric equations are

where , , and denote strain components; and are displacement components; and are electric field components; is electrical potential; and are magnetic field components; is magnetic potential.

The constitutive equations are

where , , and are the elastic, dielectric, and magnetic permeability coefficients, respectively; , , and are piezo-electric, piezomagnetic, and magnetoelectric coefficients, respectively. For FGMEE materials, we have

where is an arbitrary function; and , , , , and are the values on the plane = 0.

3. FGMEE Materials

An FGMEE structure is characterized by the high heterogeneity of material properties with a distribution prescribing the volume fractions of constituent phases. For particular analysis, it is functional to idealize them as continua with smooth gradual variation of material properties in the spatial coordinates. Hence, the proper micromechanical model should be able to characterize the material property distribution of a system in accurate sense.

Previous literatures focus on two types of gradation methods widely applied to solve many problems. Among various methods for composites, some are also used for FGMEE materials, including the exponential and Voigt rule of mixture scheme.

For FGMEE materials with exponential variation in the thickness direction (-direction), (4) can be rewritten as

where is the exponential factor governing the degree of -direction gradient, is the thickness, the superscript 0 indicates the -independent coefficients, and = 0 in homogeneous MEE materials.

The volume fraction of an FGMEE structure across the thickness direction is assumed as a simple power law type as follows:

where and is the power law index. The bottom surface of the material () is whereas the top surface () is . The total volume of the constituents should be

Based on (6) and (7), the effective material property is defined as follows:

where “MC” is general notation for material property. With (3), the effective coefficients can be written as

where “eff” stands for effective properties corresponding to a specific value of . Material coefficients of piezoelectric BaTiO₃, magnetostrictive CoFe₂O₄, and MEE BatiO₃-CoFe₂O₄ are given in Table 1. Figure 2 depicts the through-the-thickness distribution of the volume fraction changing with different values of . For = 1.0, the variation of effective material property is linear.

Table 1: Magnetoelectroelastic coefficients of material properties [6].

Material constants
CoFe₂O₄
BatiO₃-CoFe₂O₄
BatiO₃

10⁹ N/m²
286
200
166

10⁹ N/m²
173
110
77

10⁹ N/m²
170
110
78

10⁹ N/m²
269.5
190
162

10⁹ N/m²
45.3
45
43

C/m²
0
−3.5
−4.4

C/m²
0
11
18.6

C/m²
0
0
11.6

10⁻⁹ C/Vm
0.08
0.9
11.2

10⁻⁹ C/Vm
0.093
7.5
12.6

10⁻⁴ Ns²/C²
−5.9
−1.5
0.05

10⁻⁴ Ns²/C²
1.57
0.75
0.1

N/Am
580
200
0

N/Am
700
260
0

N/Am
560
180
0

10⁻¹² Ns/VC
0
6.0
0

10⁻¹² Ns/VC
0
2500
0

kgm⁻³
5730
5730
5730

Figure 2: Variation of the volume fraction function versus the nondimensional thickness with varying .

4. ICS-FEM

The solution domain is discretized into elements containing nodes, the approximation displacement , the approximation electrical potential , and the approximation magnetic potential . For FGMEE materials, we have

where u, , and are the vectors of node displacement, node electrical potential, and node magnetic potential, respectively; , , and are displacement shape, electrical potential shape, and magnetic potential shape functions of ICS-FEM, respectively. , , and were expressed in similar shape functions. Four-node element divided into four smoothing subdomains [27], field nodes, edge smoothing nodes, center smoothing nodes, edge Gaussian point, outer normal vector distribution, and shape function values are shown in Figure 3.

Figure 3: Smoothing subdomains and the values of shape functions.

At any point in the smoothing subdomain , the smoothed strain , smoothed electric field , and smoothed magnetic field are

where , , and are the strain, electric field, and magnetic field in FEM, respectively; is the constant function:

where

Substituting (12) into (11), we get

where is the boundary of ; , , and are the outer normal vector matrices of the boundary:

Eqs. (14) can be rewritten as

where is the number of smoothing elements

At the Gaussian point of the smoothing boundary , (17) are rewritten as

where is the length of the smoothing boundary; is the total number of boundaries for each smoothing subdomain.

As for the essential difference, FEM has to derive the shape function matrix of the element, but ICS-FEM avoids this step and simply uses the shape function at , which reduces the requirement for continuity of the shape function and improves the accuracy and convergence by the use of gradient smoothing.

The thermodynamic potential of a 2D FGMEE problem is given as

where S, E, and H are independent variables of strain, electric field, and magnetic field, respectively.

By applying (1) into (19), we get the variational expression of MEE plane:

By minimizing (20) for nodal variables of shape functions for strain-displacement, electric field–electric potential, and magnetic field–magnetic potential, we get the ICS-FEM equations for MEE plane:

where is the eigenvalues.

Different elemental stiffness matrices used for FGMEE beams are expressed as follows:

where ; is the mass of smoothing element ; is the smoothing element thickness; is density of Gaussian integration point in smoothing subdomain ; , , , , , and are the matrices of elastic constant, dielectric coefficient, magnetic permeability, piezomagnetic coefficient, piezomagnetic coefficient, and magnetoelectric coefficient, respectively.

The inhomogeneous smoothing element was adopted to calculate its stiffness matrix. Because the parameters of four smoothing subdomains differed in element , the actual parameters at the Gaussian integration point were taken directly to simulate the changes of material property in each element.

By eliminating the terms of electric and magnetic potentials using a condensation technique, we get the equivalent stiffness matrix :

where

The eigenvectors corresponding to and are given as

To study the effect of magnetoelectric constant on system frequencies, we derived by neglecting the magnetoelectric coupling effect.

By making , we get the reduced cell-based finite element equations:

The reduced stiffness matrix is

To evaluate the effect of PE phase on beam frequency, we derived the stiffness matrix by setting the magnetic potential = 0:

To study the magnetic effect of PM phase on system frequency, we obtained by plugging electric potential to zero in (26):

5. Results and Discussion

5.1. C-F Beam

The free vibrations on FGMEE beams were calculated by changing the exponential factor (Figure 4). The material properties of FGMEE beams were governed by the -direction exponential variation. The following geometrical parameters were considered: length = 0.3 m and width = 0.02 m with the assumption of plane stress. Boundary conditions were = = = = 0 at the clamped end. Table 2 gives the various structural frequencies in the study.

Table 2: List of f used in the study.

Structural frequency
Matrix used to compute the structural frequency

Figure 4: Geometry of an FGMEE beam and the coordinates.

Firstly, the convergence of ICS-FEM was verified by using BaTiO₃-CoFe₂O₄ FGMEE beams, with properties listed in Table 1. The natural frequencies of these beams were calculated using ICS-FEM with different meshes (30 × 2, 60 × 4, 120 × 8, 150 × 10) (Figure 5). The simulation results with different meshes agree well, which prove the good convergence of ICS-FEM. The first- to third-order modes of clamp-free BaTiO₃–CoFe₂O₄ FGMEE beams using different elements with exponential factor = 1.0 and 5.0 were calculated by ICS-FEM, and the results were summarized in Figures 6 and 7, respectively. It was found that the first- to third-order modes of BaTiO₃–CoFe₂O₄ FGMEE beams in the same gradient distribution were basically not affected by equivalent stiffness matrix or mesh number. the first- to third-order modes of BaTiO₃–CoFe₂O₄ FGMEE beams were basically not different between = 1.0 and 5.0.

(a)

(b)

(c)

(d)

(e)

Figure 5: Structural frequencies (a) ; (b) ; (c) ; (d) ; (e) of clamp-free BaTiO₃–CoFe₂O₄ FGMEE beams.

Figure 6: The first- to third-order modes of clamp-free BaTiO₃–CoFe₂O₄ FGMEE beams using different elements with exponential factor = 1.0 by ICS-FEM.

Figure 7: The first- to third-order modes of clamp-free BaTiO₃–CoFe₂O₄ FGMEE beams using different elements with exponential factor = 5.0 by ICS-FEM.

Secondly, the free vibration frequencies of CoFe₂O₄ FGMEE beams were studied by both ICS-FEM and FEM using extremely irregular elements, with domain discretization shown in Figure 8. The frequencies of clamp-free CoFe₂O₄ FGMEE beams with different values of exponential factor are shown in Figure 9; the first eleven natural frequencies calculated by ICS-FEM are smaller than those calculated by FEM. The validity of ICS-FEM is verified by the agreements between the calculations and the reference solutions. The shape of quadrilateral element in FEM cannot be severely distorted but was eliminated in ICS-FEM. ICS-FEM abstains from calculating the derivative of the shape functions of an element, and the area integral of the solution domain is converted to the boundary integral. The stiffness of FGMEE structures is improved because ICS-FEM does not require continuity of the shape function. The ICS-FEM provides a continuous system with a close-to-exact stiffness, which could be automatically and more easily generated for complicated domains, thus significantly decreasing the numerical error. The free vibration of CoFe₂O₄ FGMEE beam, a pure CoFe₂O₄ material without piezoelectric or magnetoelectric material coefficients, influences structural frequency because the magnetic effect is marginally higher compared with . coincides with since piezoelectric phase is absent in CoFe₂O₄. Similarly, coincides with of the CoFe₂O₄ FGMEE beam as the magnetoelectric effect is absent in pure CoFe₂O₄ FGMEE beam. The natural frequencies of CoFe₂O₄ FGMEE beams increase with the rise of the exponential factor.

Figure 8: Domain discretization using four-node extremely irregular elements.

(a)

(b)

(c)

(d)

(e)

Figure 9: Structural frequencies of clamp-free CoFe₂O₄ FGMEE beams with exponential factor = (a) 0, (b) 0.5, (c) 1.0, (d) 2.0, and (e) 5.0.

Finally, the comparison of calculation time between ICS-FEM and FEM at Intel (R) Xeon (R) CPU E3-1220 v3 @ 3.10 GHz, 16 G RAM is shown in Figure 10, with the element number of 60, 240, 960, and 1500. As showed in Figure 10, the time required to solve algebraic equations by ICS-FEM is similar to that of FEM. Because the stiffness construction of ICS-FEM is based on smoothing cells inside each element, no coupling occurs between nodal degrees-of-freedom that are the distance of up to two elements. In other words, the bandwidth of ICS-FEM stiffness matrix is the same as that of FEM. Nevertheless, ICS-FEM is more effective in terms of generalized displacement (including displacement, electrical potential and magnetic potential) and computational efficiency (computation time for the same accuracy).

Figure 10: Comparison of computational efficiency.

5.2. S-S Beam

The free vibrations on BaTiO₃–CoFe₂O₄ FGMEE beams were studied by changing the gradient function form in Figure 4. The geometrical parameters were the same as the C-F beam. The simply supported boundary conditions were used: at (, ) and at (, ).

Firstly, the free vibration frequencies for BaTiO₃–CoFe₂O₄ FGMEE beams were calculated by both ICS-FEM with 60 × 4 meshes and FEM with 240 × 16 meshes (Figure 11). Results show the first eleven natural frequencies calculated by ICS-FEM are closer to the reference solutions than those calculated by FEM, indicating ICS-FEM is more efficient than FEM due to the reduced number of meshes. The differences in natural frequencies , , and are marginal, so the magnetic effect only slightly impacts the natural frequencies of FGMEE beams. The natural frequencies increase with the rise of the exponential factor.

(a)

(b)

(c)

(d)

(e)

Figure 11: Structural frequencies for simply supported BaTiO₃-CoFe₂O₄ FGMEE beams with exponential factor = (a) 0, (b) 0.5, (c) 1.0, (d) 2.0, and (e) 5.0.

Secondly, as for BaTiO₃–CoFe₂O₄ FGMEE beams with power law type (, 1.0, 5.0, 10.0, 20.0), the bottom surface () is BaTiO₃–CoFe₂O₄ whereas the top surface () is CoFe₂O₄. The free vibration frequencies for FGMEE beams calculated by ICS-FEM with 60 × 4 meshes and FEM with 240 × 16 meshes are shown in Figure 12. Results show the first eleven natural frequencies calculated by ICS-FEM are closer to the reference solutions than those calculated by FEM. Also the differences in , , and are marginal, so the magnetic effect does not largely impact the natural frequencies of FGMEE beams. Meanwhile, the natural frequency of BaTiO₃–CoFe₂O₄ FGMEE beams is between those of BaTiO₃–CoFe₂O₄ MEE beams and CoFe₂O₄ MEE beams.

(a)

(b)

(c)

(d)

(e)

Figure 12: Structural frequencies of simply supported BaTiO₃-CoFe₂O₄ FGMEE beams with power law index = (a) 0, (b) 1.0, (c) 5.0, (d) 10.0, and (e) 20.0.

5.3. Typical MEMS-Based FGMEE Energy Harvester

The model of FGMEE energy harvester developed by ICS-FEM is shown in Figure 13. The free vibrations on the FGMEE energy harvester were studied by changing the exponential factor. The geometrical parameters were  mm,  mm, and its structure was fabricated with BaTiO₃ FGMEE.

(a)

(b)

Figure 13: Typical MEMS-based energy harvester fabricated with BaTiO₃ FGMEE. (a) Model of the energy harvester; (b) simplified model of ICS-FEM.

The free vibration frequencies for the FGMEE energy harvester calculated by ICS-FEM with 60 × 4 meshes and FEM with 240 × 16 meshes are shown in Figure 14. The first eleven natural frequencies calculated by ICS-FEM are closer to the reference solutions than those calculated by FEM, indicating ICS-FEM is more efficient than FEM owing to the reduced number of meshes. The ICS-FEM does not take the derivative of the shape functions of the element and can be much easily generated automatically for complicated domains, thus significantly decreasing the numerical errors. The natural frequencies and agree well with each other since the piezomagnetic phase is absent from the BaTiO₃ FGMEE energy harvester. Moreover, coincides with of the BaTiO₃ FGMEE energy harvester as the magnetoelectric effect is absent in pure BaTiO₃ materials. The and are very close, so the piezoelectric effect only slightly affects the natural frequencies of the energy harvester, which increase with the rise of the exponential factor.

(a)

(b)

(c)

(d)

(e)

Figure 14: Structural frequencies for BaTiO₃ FGMEE energy harvester with exponential factor = (a) 0, (b) 0.5, (c) 1.0, (d) 2.0, and (e) 5.0.

The Wilson-θ method and the equivalent stiffness matrix [] were employed to solve the dynamic response of the FGMEE energy harvester. The parameters were set as time step = 0.005 s, ; without damping; sine-wave transient load with a time period of 2 s; 4 cycles of loading (Figure 15). The dynamic behaviors for the harvester calculated by ICS-FEM with 60 × 4 meshes and FEM with 240 × 16 meshes are shown in Figures 16 and 17, respectively. The temporal variations of displacement and electric potential calculated by ICS-FEM are closer to the reference solutions than those by FEM, which validate the accuracy of ICS-FEM. The temporal variations of and Φ of the FGMEE energy harvester decrease with the increase of the exponential factor when the material properties of the harvester are in exponential distribution.

Figure 15: Sine-wave load at point of the BaTiO₃ FGMEE energy harvester.

Figure 16: Temporal variations of displacement at point of BaTiO₃ FGMEE energy harvester with different values of exponential factor calculated by ICS-FEM and FEM.

Figure 17: Temporal variations of electric potential at point of BaTiO₃ FGMEE energy harvester with different values of exponential factor calculated by ICS-FEM and FEM.

6. Conclusions

The free vibrations on FGMEE structures were studied. Firstly, ICS-FEM for FGMEE materials was formulated by incorporating gradient smoothing into the FEM-based computation for the FGMEE multi-physics field. Then the equations of free vibration computation were deduced for the multi-physics coupling field of FGMEE materials. Finally, the FGMEE beams were calculated with functional gradient exponential form or power law form under different boundary conditions.

(i) ICS-FEM reduced the systematic stiffness of the finite element, which improved the computational accuracy compared with FEM under the same element number. ICS-FEM was more efficient than FEM in terms of computation time for the same accuracy.

(ii) Due to the material property changes in each smoothing element, the true parameters at the Gaussian integration point were adopted directly. ICS-FEM avoided the derivative of the shape functions but only transformed the area integral to the boundary integral in the solution domain, which omitted the requirement of continuity of the shape function.

(iii) The magnetic effect slightly influenced the natural frequencies of FGMEE beams, which increased with the increase of the exponential factor when the material properties of FGMEE beams were under exponential distribution. The natural frequency of FGMEE beams lied in between those of the FGMEE beams using the bottom surface of materials and the FGMEE beams using the upper surface, when the material properties of FGMEE beams were under the power law distribution.

(iv) The natural frequencies and general displacements of the FGMEE energy harvester developed by ICS-FEM were more accurate compared with FEM, owing to the reduced number of meshes.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

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 Natural Science Foundation of China (Grant no. 11502092), Jilin Provincial Department of Science and Technology Fund Project (Grant nos. 20160520064JH, 20170101043JC), and the Fundamental Research Funds for the Central Universities.

References

W. Eerenstein, N. D. Mathur, and J. F. Scott, “Multiferroic and magnetoelectric materials,” Nature, vol. 442, no. 7104, pp. 759–765, 2006.

J. F. Scott, “Applications of magnetoelectrics,” Journal of Materials Chemistry, vol. 22, no. 11, pp. 4567–4574, 2012.

A. Zhang and B. Wang, “Effects of crack surface electrostatic tractions on the fracture behaviour of magnetoelectric composite materials,” Mechanics of Materials, vol. 102, pp. 15–25, 2016.

J. Zhai, Z. Xing, S. Dong, J. Li, and D. Viehland, “Detection of pico-Tesla magnetic fields using magneto-electric sensors at room temperature,” Applied Physics Letters, vol. 88, no. 6, Article ID 062510, 2006.

D. J. Huang, H. J. Ding, and W. Q. Chen, “Analytical solution for functionally graded magneto-electro-elastic plane beams,” International Journal of Engineering Science, vol. 45, no. 2-8, pp. 467–485, 2007.

A. Milazzo, C. Orlando, and A. Alaimo, “An analytical solution for the magneto-electro-elastic bimorph beam forced vibrations problem,” Smart Materials and Structures, vol. 18, no. 8, Article ID 085012, 2009.

W. Q. Chen, K. Y. Lee, and H. J. Ding, “On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates,” Journal of Sound and Vibration, vol. 279, no. 1-2, pp. 237–251, 2005.

C.-P. Wu and Y.-C. Lu, “A modified Pagano method for the 3D dynamic responses of functionally graded magneto-electro-elastic plates,” Composite Structures, vol. 90, no. 3, pp. 363–372, 2009.

A. R. Annigeri, N. Ganesan, and S. Swarnamani, “Free vibration behaviour of multiphase and layered magneto-electro-elastic beam,” Journal of Sound and Vibration, vol. 299, no. 1-2, pp. 44–63, 2007.

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.

R. K. Bhangale and N. Ganesan, “Free vibration studies of simply supported non-homogeneous functionally graded magneto-electro-elastic finite cylindrical shells,” Journal of Sound and Vibration, vol. 288, no. 1-2, pp. 412–422, 2005.

R. K. Bhangale and N. Ganesan, “Free vibration of simply supported functionally graded and layered magneto-electro-elastic plates by finite element method,” Journal of Sound and Vibration, vol. 294, no. 4, pp. 1016–1038, 2006.

R. G. Lage, C. M. M. Soares, C. A. M. Soares, and J. N. Reddy, “Layerwise partial mixed finite element analysis of magneto-electro-elastic plates,” Computers & Structures, vol. 82, pp. 1293–1301, 2004.

S. S. Phoenix, S. K. Satsangi, and B. N. Singh, “Layer-wise modelling of magneto-electro-elastic plates,” Journal of Sound and Vibration, vol. 324, no. 3-5, pp. 798–815, 2009.

G. R. Buchanan, “Free vibration of an infinite magneto-electro-elastic cylinder,” Journal of Sound and Vibration, vol. 268, no. 2, pp. 413–426, 2003.

J. Sladek, V. Sladek, S. Krahulec, C. S. Chen, and D. L. Young, “Analyses of circular magnetoelectroelastic plates with functionally graded material properties,” Mechanics of Advanced Materials and Structures, vol. 22, no. 6, pp. 479–489, 2015.

J. Liu, P. Zhang, G. Lin, W. Wang, and S. Lu, “High order solutions for the magneto-electro-elastic plate with non-uniform materials,” International Journal of Mechanical Sciences, vol. 115-116, pp. 532–551, 2016.

L.-L. Ke and Y.-S. Wang, “Free vibration of size-dependent magneto-electro-elastic nanobeams based on the nonlocal theory,” Physica E: Low-Dimensional Systems and Nanostructures, vol. 63, pp. 52–61, 2014.

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.

S. P. 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.

G. R. Liu and G. Y. Zhang, “A normed G space and weakened weak (W-2) formulation of a cell-based smoothed point interpolation method,” International Journal of Computational Methods, vol. 6, no. 1, pp. 147–179, 2009.

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.

H. Feng, X. Cui, and G. Li, “A stable nodal integration method for static and quasi-static electromagnetic field computation,” Journal of Computational Physics, vol. 336, pp. 580–594, 2017.

S. Li, X. Cui, H. Feng, and G. Wang, “An electromagnetic forming analysis modelling using nodal integration axisymmetric thin shell,” Journal of Materials Processing Technology, vol. 244, pp. 62–72, 2017.

S. Natarajan, A. J. M. Ferreira, S. P. A. Bordas, E. Carrera, and M. Cinefra, “Analysis of composite plates by a unified formulation-cell based smoothed finite element method and field consistent elements,” Composite Structures, vol. 105, pp. 75–81, 2013.

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.

T. Nguyen-Thoi, H. C. Vu-Do, T. Rabczuk, and H. Nguyen-Xuan, “A node-based smoothed finite element method (NS-FEM) for upper bound solution to visco-elastoplastic analyses of solids using triangular and tetrahedral meshes,” Computer Methods Applied Mechanics and Engineering, vol. 199, no. 45-48, pp. 3005–3027, 2010.

T. Nguyen-Thoi, G. R. Liu, H. Nguyen-Xuan, and C. Nguyen-Tran, “Adaptive analysis using the node-based smoothed finite element method (NS-FEM),” International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 2, pp. 198–218, 2011.

T. N. Tran, G. R. Liu, H. Nguyen-Xuan, and T. Nguyen-Thoi, “An edge-based smoothed finite element method for primal-dual shakedown analysis of structures,” International Journal for Numerical Methods in Engineering, vol. 82, no. 7, pp. 917–938, 2010.

H. Nguyen-Xuan, L. V. Tran, T. Nguyen-Thoi, and H. C. Vu-Do, “Analysis of functionally graded plates using an edge-based smoothed finite element method,” Composite Structures, vol. 93, no. 11, pp. 3019–3039, 2011.

W. Li, Y. Chai, M. Lei, and G. R. Liu, “Analysis of coupled structural-acoustic problems based on the smoothed finite element method (S-FEM),” Engineering Analysis with Boundary Elements, vol. 42, pp. 84–91, 2014.

T. Nguyen-Thoi, G. R. Liu, H. C. Vu-Do, and H. Nguyen-Xuan, “A face-based smoothed finite element method (FS-FEM) for visco-elastoplastic analyses of 3D solids using tetrahedral mesh,” Computer Methods Applied Mechanics and Engineering, vol. 198, no. 41-44, pp. 3479–3498, 2009.

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.

W. Zeng, G. R. Liu, D. Li, and X. W. Dong, “A smoothing technique based beta finite element method (beta FEM) for crystal plasticity modeling,” Computers & Structures, vol. 162, pp. 48–67, 2016.

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.

Y. Onishi, R. Iida, and K. Amaya, “F-bar aided edge-based smoothed finite element method using tetrahedral elements for finite deformation analysis of nearly incompressible solids,” International Journal for Numerical Methods in Engineering, vol. 109, no. 11, pp. 1582–1606, 2017.

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.

L. M. Zhou, G. W. Meng, F. Li, and H. Wang, “Cell-based smoothed finite element method-virtual crack closure technique for a piezoelectric material of crack,” Mathematical Problems in Engineering, vol. 2015, Article ID 371083, 10 pages, 2015.

C. V. Le, “Estimation of bearing capacity factors of cohesive-frictional soil using the cell-based smoothed finite element method,” Computers & Geosciences, vol. 83, pp. 178–183, 2017.

L. M. Zhou, G. W. Meng, F. Li, and S. Gu, “A cell-based smoothed XFEM for fracture in piezoelectric materials,” Advances in Materials Science and Engineering, vol. 2016, Article ID 4125307, 14 pages, 2016.

A. L. Pramod, S. Natarajan, A. J. Ferreira, E. Carrera, and M. Cinefra, “Static and free vibration analysis of cross-ply laminated plates using the Reissner-mixed variational theorem and the cell based smoothed finite element method,” European Journal of Mechanics-A/Solids, vol. 62, pp. 14–21, 2017.

F. Wu, L. Y. Yao, M. Hu, and Z. C. He, “A stochastic perturbation edge-based smoothed finite element method for the analysis of uncertain structural-acoustics problems with random variables,” Engineering Analysis with Boundary Elements, vol. 80, pp. 116–126, 2017.

J. L. Mantari, E. M. Bonilla, and C. Guedes Soares, “A new tangential-exponential higher order shear deformation theory for advanced composite plates,” Composites Part B: Engineering, vol. 60, pp. 319–328, 2014.

A. Daga, N. Ganesan, and K. Shankar, “Behaviour of magneto-electro-elastic sensors under transient mechanical loading,” Sensors and Actuators A: Physical, vol. 150, no. 1, pp. 46–55, 2009.