Vibration of MEE Composite Conical Shell Surrounded by Nonlinear Elastic Foundation Considering the Effect of Geometrical Nonlinearity

1. Mohammadrezazadeh

Faculty of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, 1991943344, Iran

1.     Introduction

Structural mechanics is one of the interesting topics considered in numerous researches [1-11].  In addition, several studies considered nano-scale structures [12-15]. Furthermore, an example of a study on the micro-scale is the study that was done by Ouakad and ŻUR [16]. The focus of this study is on the vibration of conical shell structures on an elastic foundation composed of magneto-electro-elastic (MEE) material. Magneto-electro-elastic composite materials are smart materials that have the ability to convert magnetic, electrical, and mechanical energies to each other [17]. These materials can be used for actuators, sensors, and also vibration control purposes [17]. There are several works in literature that take attention to the study of structures containing magneto-electro-elastic materials. Bhangale and Ganesan [18] studied the vibration of functionally graded magneto-electro-elastic cylindrical shells via a finite element model. Annigeri et al. [19] investigated the vibration behavior of a magneto-electro-elastic cylindrical shell subjected to simply supported boundary conditions using a series solution and finite element model. Tsai and Wu [20] presented three-dimensional free vibration responses of simply supported doubly curved shells composed of functionally graded magneto-electro-elastic material with open-circuit surface conditions via an asymptotic method. Kumaravel et al. [21] adopted the finite element method to study the vibration and buckling of magneto-electro-elastic cylinders under clamped-clamped boundary conditions. Lang and Xuewu [22] accomplished a study about the vibration and buckling of cylindrical shells from functionally graded magneto-electro-thermo-elastic material by means of higher-order shear deformation theory. Razavi and Shooshtari [17] carried out an investigation about the free vibration of simply supported thin magneto-electro-elastic doubly-curved shells surrounded by a foundation with the help of Donnell theory. Shooshtari and Razavi [23] carried out the large amplitude vibration study of simply-supported magneto-electro-elastic curved panels on the basis of the Donnell shell theory and the Galerkin approach. Shooshtari and Razavi [24] hired the Galerkin method and Lindstedt-Poincare perturbation approach to investigate linear and nonlinear vibration of magneto-electro-elastic laminated doubly-curved thin shell resting on an elastic foundation. Mohammadimehr et al. [25] utilized first-order shear deformation theory, energy method, and Hamilton principle to research the free vibration of magneto-electro-elastic composite cylindrical panels reinforced by carbon nanotubes considering closed and open circuit boundary conditions. Vinyas and Harursampath [26] investigated nonlinear vibration of higher order shear deformable Carbon nanotube reinforced magneto-electro-elastic doubly curved shells on the basis of the Donnell shell theory and von-Karman nonlinear approach. Rostami and Mohammadimehr [27] dealt with the vibration control of a sandwich rotating cylindrical shell containing a nanocomposite face sheet and porous core and integrated with functionally graded magneto-electro-elastic layers applying differential quadrature method and first-order shear deformation theory of shells. Ye et al. [28] studied transient dynamic and free vibration of composite magneto-electro-elastic cylindrical shells via the scaled boundary finite element method. Rostami et al. [29] researched the vibration control of the sandwich rotating cylindrical shell which had functionally graded core and functionally graded magneto-electro-elastic layers on the basis of the first-order shear deformation theory of shells and differential quadrature method. Ni et al. [30] employed a Hamiltonian system to study the vibration of porous magneto-electro-thermo-elastic functionally graded cylindrical shells under magneto-electro-thermal loadings. Zhou et al. [31] obtained the time-dependent responses of MEE structures in a hygrothermal environment using the cell-based smoothed finite element technique in conjunction with the modified Newmark procedure. Farajpour et al. [32] studied large amplitude vibration of magneto-electromechanical mass nanosensors. The nanomechanical sensor utilized vibrating MEE nanoplates with numerous locations in order to trap nanoparticles [32]. It is obvious that references [17] to [32] are not about the study of conical shells.

Conical shells are from engineering structures that have a lot of applications including piping, pressure vessels, and ship structures [33]. An investigation of the free vibration of magneto-electro-elastic conical shells under clamped-free, simply supported and clamped-clamped boundary conditions based on the finite element method was done by Srikantamurthy et al. [34]. In addition, Srikantamurthy and  Annigeri [35] studied the free vibration of multiphase MagnetoElectro-Elastic uniform thickness conical shells under Clamped-Free boundary conditions. References [34] and [35] considered the vibration of MEE conical shells, but they did not study the vibration of MEE conical shells on elastic foundations.

Vibration characteristics of conical structures on linear elastic foundations were studied by many researchers [36-50]. It is necessary to mention that the behavior of the foundation is generally nonlinear [51]. The author found that there are limited papers in the literature about the vibration of conical structures on nonlinear elastic foundations. Zhu et al. [52] accomplished the smart control of large amplitude vibration of porous piezoelectric sandwich conical panels surrounded by nonlinear elastic foundation via first-order shear deformation theory, von Kármán nonlinear approach, and harmonic balance method. The considered conical panel contained a viscoelastic core as well as two porous piezoelectric layers [52]. Molla-Alipour et al. [53] studied free vibration of the bidirectional functionally graded cylindrical and conical shells as well as annular plates on nonlinear elastic foundations. References [52] and [53] investigated the vibration of conical structures on a nonlinear elastic foundation, but the conical structures are not from MEE material.

To the best of the author's knowledge, there is not any study in the literature which considers nonlinear vibration of simply supported MEE composite conical shells surrounded by the nonlinear elastic foundation. This study considers the nonlinear vibration responses of simply supported MEE composite conical shells on nonlinear elastic foundations subjected to electric or magnetic potential. The nonlinearity of the system is due to the geometric nonlinearity modeled via the von Karman approach as well as foundation nonlinearity. The effects of shear deformation as well as rotary inertia are taken into account in the process of shell modeling.

Unfortunately, it is not possible to obtain an exact solution for a great number of mathematical problems. Numerical methods are effective methods that are extremely used for the solution of mathematical problems. For example, the discrete singular convolution method which is a numerical method is used in several researches such as [54-56]. Regardless of the advantages of the numerical methods, analytical methods are used in order to have parametric studies, consider the physics of problems and validate the numerical results [57-58]. In reference [59], a semi-analytical method was employed in order to solve equations. Because of the advantages of approximate analytical methods, in this study, the Lagrange method as an approximate analytical method is handled to convert the partial differential equations of the conical shell to an ordinary differential equation. Lindstedt-Poincare method and modal analysis are employed to obtain the nonlinear response of the system. Several results from published literature are compared with the results of this study for validation purposes. In the next step, the effects of several parameters including nonlinear and linear constants of elastic foundation, electric and magnetic potentials, thickness, and length on the fundamental linear frequency as well as nonlinear vibration responses are illustrated.

2.     Modeling of the Shell

2.1. Strain Relations

Figure 1 shows the coordinate system and schematic of the considered MEE composite conical shell. In addition, a schematic of the conical shell on a nonlinear elastic foundation is depicted in Figure 2.

The longitudinal, circumferential, and normal directions of the coordinate system are shown with , and , respectively. The variable  depicts the length, is the thickness,  implies a semi-vertex angle and and  are respectively the radii of small and large edges on the middle surface of the shell. The radius of each point on the middle surface of the shell is denoted by .

Equation (1) indicates the relations of strains ( ) with membrane strains
( ) as well as curvatures ( ) [61]:

Considering the von Karman approach for geometrical nonlinearity and shear deformation leads to equations (2) and (3) for membrane strains and curvatures, respectively [61-62]. In addition, shear strains in the middle surface of the conical shell are extracted through equation (4) [61]:

Fig. 1. Image of the schematic of the MEE composite conical shell and coordinate system [60]

Fig. 2. Image of schematic of the conical shell on the nonlinear elastic foundation [39]

It should be mentioned that  , and denote the displacements of an arbitrary point on the middle surface of the conical shell located at point ( ) [61]. In addition, and refer to the total angular rotations of the line normal to the middle surface about  and  axes, respectively [61].

2.2. Stress-strain Relations

The coupled constitutive relations for the MEE shell can be expressed as following [63, 64 quoted from 65-68]:

while , , and  refer to stress, strain, electric field, and magnetic field vectors, respectively [64]. Besides,  and  denote electric displacement and magnetic induction vectors, respectively [63-64]. In addition,  , and imply matrices of elastic, piezoelectric, and piezomagnetic coefficients, respectively [63-64]. It should be mentioned that ,  and  are respectively the matrices of dielectric, magnetoelectric, and magnetic permeability coefficients [64].

If the electric vector and magnetic intensity are respectively written as gradients of the scalar electric ( ) and magnetic ( ) potentials, vector equations of Maxwell in the quasi-static approximation are satisfied [63]:

In this paper, because of the thin nature of the shell, the components of in-plane electric and magnetic fields are ignored
( ) [69 quoted from 70].

Equations (9) and (10) show Gauss’ laws for electrostatics and magnetostatics which can be used to obtain the electromagnetic state of the shell [17 quoted from 71]:

Considering magneto-electric boundary conditions ( , , , ) [17] and the previous assumptions that , ,  and  lead to the following relations for and :

while  and  refer to electric and magnetic potentials, respectively [17].

3.     Lagrange Method

In order to extract the nonlinear ordinary differential equation of the system, for the considered problem, the Lagrange method equation [61] can be used:

Whereas . It should be mentioned that and  refer to kinetic and strain energies [61] while denotes energy due to a nonlinear foundation. Required relations of these energies can be obtained through equations (13) to (15):

It is valuable to emphasize that equations (13) and (14) are written based on reference [61] while equation (15) is derived on the basis of references [37, 51]. In equation (14), , and denote mass moments of inertia terms [62]. It should be noted that in equation (15), and demonstrate linear and nonlinear constants of the Winkler elastic foundation while indicates the shear stiffness constant of the elastic foundation [51]. In order to derive the nonlinear ordinary differential equation of the vibration, the following relations which satisfy geometrical boundary conditions of the simply supported conical shells    [61, 72] can be used:

The displacement and rotation relations of equation (16) are written based on the reference [61] while and  demonstrate half-wave numbers of displacements in the length and circumference of the shell, respectively [61], and  refers to the vibration mode shape. Substituting equation (16) into equations (13) to (15) and then using equation (12) lead to extraction of the nonlinear ordinary differential equation of the system:

while  and  respectively demonstrate mass and stiffness matrices and  is the vector of nonlinear coefficients caused by geometric nonlinearity and nonlinearity of elastic foundation.

4.     Lindstedt-Poincare Method and Modal Analysis

In order to use the Lindstedt-Poincare method, a new independent variable is defined [73]; thus the relation is acquired. In addition, the variables , and should be considered as following [73]:

while  denotes the fundamental linear frequency ( ) of the conical shell. Using mentioned relations, one can rewrite equation (17) as shown in equation (20):

By performing some mathematical operations and setting the coefficients of , and equal to zero, equations (21) to (23) are acquired which are respectively coefficients of , and :

In order to obtain the response of equation (21), applying modal analysis relation  [61] results in equation (24). It is important to mention that and denote modal matrix and modal coordinates vector, respectively [61].

The row of equation (24) which includes  term of matrix  and its response is as following:

Substituting modal analysis relations and  [61] into equation (22) leads to the acquisition of equation (26):

Substituting equation (25) into the row of equation (26) which contains component of the matrix eventuates to and particular solution of . Putting modal analysis relations ,  and  [61] into equation (23) and doing some simplifications lead to the obtaining of equation (27):

Neglecting the terms that contain other frequencies than  from the row of equation (27) which contains term of matrix and utilizing equation (25) and also and  lead to the following equation:

In order to derive a finite response for equation (28), it is necessary to put the coefficient of equal to zero                       ( ) which leads to  and also a particular solution of . Finally, the term of  which is only related to fundamental linear frequency ( ) can be obtained as:

In addition, putting the results obtained for ,  and  into equation (18) leads to obtain of the following relation for nonlinear frequency ratio ( ):

Equation (30) indicates that the nonlinear frequency ratio has a direct relation with nonlinear parameter . Furthermore,  has a direct relation with the square of the amplitude parameter ( ) while its relation with the square of the fundamental linear frequency is inverse.

5.     Results and Discussions

The purpose of this paper is the investigation of the nonlinear vibration of MEE composite conical shells on a nonlinear elastic foundation under electric or magnetic potential. In order to validate this study, Table 1 compares the frequency parameter ( ) results of this study with literature (references [74] and [75]) for an isotropic simply supported conical shell with constants of and while and  denote Young modulus and Poisson ratio, respectively. According to this table, one can conclude that in most cases there is good agreement between the results of this study and the literature.

Table 2 shows the frequency results (Hz) obtained for an isotropic cylindrical shell on a linear elastic foundation with constants of ,  ,  , ,  ,  and  and compares the results with the results of reference [76]. It is noteworthy to mention that   and   refer to the radius of the cylindrical shell on the middle surface and the mass density, respectively. The results of Table 2 are for different values of the constants of the foundation and half-wave number of displacement in the circumference of the shell ( ). This table demonstrates very good agreement between the results of this paper and the reference [76] which could be evidence for the validity of the present research.

After validation studies, it is time to investigate the vibration responses of the MEE conical shell on a nonlinear elastic foundation. It is important to mention that BaTiO₃-CoFe₂O₄ composite material is chosen as MEE material. The constants used in this study are extracted based on the constants introduced in reference [64] for piezoelectric BaTiO₃ material and magnetostrictive CoFe₂O₄ material unless the values of the densities of the materials.

The values of densities are  for CoFe₂O₄ and for BaTiO₃ [77]. Because of the importance of obtaining of the constants of BaTiO₃-CoFe₂O₄ composite material, some explanations are required. References [64] and [77] are provided the constants of BaTiO₃ and CoFe₂O₄ materials, separately. In order to obtain the constants of BaTiO₃-CoFe₂O₄ composite material, in this paper, the average of the constants of BaTiO₃ and CoFe₂O₄ materials is used. The constants of reference [64] are for a situation that normal strain ( ) and normal stress ( ) are considered. In the present paper, and ; so it is necessary to do some operations which are based on reference [62].

The equation of MEE material can be written as  [64 quoted from 65-68] which can be easily converted to . Therefore, it is necessary to obtain matrices ,  as well as and eliminate rows and columns that correspond to and . Doing this simple operation leads to obtain of the constants of the MEE composite material as shown in Table 3:

Table 1. Comparison of frequency parameter results of this study and literature for isotropic conical shells
with different semi-vertex angle values

Table 2. Comparison of the frequencies (Hz) of the present study with literature for isotropic cylindrical shells
surrounded by elastic foundation

Table 3. The values of the constants of the composite
MEE material

In addition, the average of the densities is derived as .

It is necessary to mention that the results of this study are acquired for MEE composite conical shell with , ,  and  on nonlinear elastic foundation with constants ,  and . All constants used in this paper are as mentioned unless other values are emphasized. In addition, all results for MEE composite conical shell are acquired for the mode shape .

Table 4 presents the values of fundamental linear frequency as well as nonlinear parameter
( ) for different values of the constants of the elastic foundation for shells subjected to electric or magnetic potential.

This table indicates that the value of fundamental linear frequency doesn’t have any change with the change of the value of . This is because of the fact that is a nonlinear parameter; so it has no effect on the fundamental linear frequency. On the other hand, according to this table, one can conclude that as one of the linear constants (  or ) of the elastic foundation increases, the fundamental linear frequency gets higher values; because the increase of or increases the linear stiffness of the system which leads to greater fundamental linear frequency. Also, Table 4 indicates that increases with an increase of while an increase of the linear constants of the elastic foundation ( or ) causes a decrease of . This is because of the fact that  is a nonlinear parameter and its increase leads to the increase of the nonlinearity of the system which appears in . On the other hand, when  or , which are linear parameters, increases, the rate of the nonlinearity of the system decreases which leads to the decrease of .

Figures 3 (a) and (b) depict the effect of the nonlinear constant of the elastic foundation on the curves of nonlinear frequency ratio versus amplitude parameter in the presence of electric and magnetic potentials, respectively. These figures illustrate that in the presence of electric or magnetic potential, the increase of the nonlinear constant of the elastic foundation causes the increase of the nonlinear frequency ratio. This is because of the fact that as shown in Table 4, the increase of  leads to the increase of the nonlinear parameter ( ) which according to equation (30) causes a greater nonlinear frequency ratio.

Table 4. The influence of the constants of the elastic foundation on the vibration characteristics of the considered
MEE composite conical shell

Figure 3. The influence of the nonlinear constant of the foundation on the  diagrams in the presence of,
a: electric potential, b: magnetic potential

Figures 4 (a) and (b) show diagrams of nonlinear frequency ratio versus amplitude parameters for different values of  in the presence of electric and magnetic potentials, respectively. These figures demonstrate that the change of has a low impact on the nonlinear frequency ratio versus amplitude parameter diagrams. However, more precision in Figures 4 (a) and (b) reveals that, for the constant amplitude parameter, the increase of leads to the decrease of the nonlinear frequency ratio. This is due to the reason that the increase of leads to greater linear stiffness of the system and so as depicted in Table 4, greater fundamental linear frequency.

In addition, according to Table 4,  gets smaller values with an increase of . As shown in equation (30), the mentioned facts lead to the decrease of the nonlinear frequency ratio of the considered MEE composite conical shell.

Figures 5 (a) and (b) display diagrams of nonlinear frequency ratio against amplitude parameters in the presence of electric and magnetic potentials, respectively, while  getting different values. Looking at these figures leads to the conclusion that for a constant value of , the increase of  leads to a decrease in the value of the nonlinear frequency ratio.

This conclusion can be explained as following: As shown in Table 4, the increase of causes smaller values for  and greater values for fundamental linear frequency. It is apparent from equation (30) that the decrease of and the increase of result in a smaller nonlinear frequency ratio.

In addition to the noted results, Figures 3 to 5 indicate that the greater the amplitude parameter, the greater the nonlinear frequency ratio; because the increase of the amplitude parameter as shown in equation (30), makes the nonlinear effects more apparent.

Figure 4. Diagrams of  for different values of the foundation constant  when, a:  , b: , exists

Figure 5. The curves of nonlinear frequency ratio against amplitude parameter for different values of in the condition that,
 a:  , b:

Figure 6 represents a diagram of the nonlinear frequency ratio versus the electric potential for different values of amplitude parameter while . It may be helpful to mention that [-90000, -9000, -900, -90, -9, 0, 9, 90, 900, 9000, 90000] V.   Figure 6 shows that as the electric potential moves from -90000 to 90000, the nonlinear frequency ratio gets higher values. The reason for this behavior can be explained as follows: the movement of electric potential from negative values to positive ones leads to a decrease of the linear part of stress which strengthens the nonlinearity of the system; this is apparent from equations (5) and (11) as well as Table 3. Besides, in order to illustrate the reason for this behavior, Table 5 could be helpful.  Table 5 shows the values of fundamental linear frequency and nonlinear parameters for different values of electric potential while the magnetic potential is considered to be zero. This table indicates that the nonlinear parameter increases and in most cases, the fundamental linear frequency decreases with the movement of electric potential from -90000 V to 90000 V which according to equation (30), justifies the increase of the nonlinear frequency ratio in the considered situations. Another result of Figure 6 is that for a constant value of  , as the value of the amplitude parameter increases, the nonlinear frequency ratio gets higher values.

Figure 6. The curves of nonlinear frequency ratio against the electric potential for different amplitude parameter values

Table 5. The impression of the electric potential on the values of fundamental linear frequency and nonlinear parameter ( )

Figure 7 depicts the curves of nonlinear frequency ratio versus the magnetic potential for various amplitude parameter values while and [-90000, -9000, -900, -90, -9, 0, 9, 90, 900, 9000, 90000] A.

Figure 7 shows that for constant amplitude parameter, the increase of the magnetic potential from -90000 A to 90000 A results in the decrease of the nonlinear frequency ratio. This is because of the fact that the movement of magnetic potential from -90000 A to 90000 A leads to an increase of the linear part of stress which decreases the effect of nonlinearity; according to equations (5) and (11) and also Table 3, this result is obvious. Furthermore, the reason for this behavior can be illustrated with the help of Table 6 as following: Table 6 demonstrates the values of fundamental linear frequency and nonlinear parameters for different magnetic potential values whereas . This table shows that with the movement of the magnetic potential from -90000 A to 90000 A, the fundamental linear frequency gets higher values because of the increase of the linear stiffness while the nonlinear parameter becomes smaller. It is apparent from equation (30) that the mentioned factors cause a decrease in the nonlinear frequency ratio. In addition, one can conclude from Figure 7 that for the constant value of magnetic potential, the greater the amplitude parameter, the higher the nonlinear frequency ratio.

Figure 7. The curves of  versus magnetic potential for different values of amplitude parameter

Table 6. The values of vibration characteristics for various values of magnetic potential

Figure 8 (a) demonstrates the curves of nonlinear frequency ratio versus amplitude parameter in the presence of positive electric potential, in the condition without electric and magnetic potentials as well as in the presence of positive magnetic potential. This figure implies that for constant , the presence of positive magnetic potential leads to the smallest value of the nonlinear frequency ratio. Besides, a comparison of the responses of the system in the presence of positive electric potential and the system without electric and magnetic potentials reveals that for the constant amplitude parameter, the value of the nonlinear frequency ratio is slightly greater for the condition that . The reason for the mentioned results of Figure 8 (a) can be illustrated as follows: in the presence of positive electric potential, the linear part of the stress is smaller in comparison with the situation without electric and magnetic potentials. On the other hand, in the presence of positive magnetic potential, the linear part of stress is higher in comparison with the situation without electric and magnetic potentials. It should be mentioned that greater linear stress leads to a decrease in the nonlinearity of the system. As can be seen, the effect of the magnetic potential on the diagram is more apparent; this is because of the fact that according to Table 3, the value of  is much greater than the absolute value of . Figure 8 (b) shows the nonlinear frequency ratio against amplitude parameter diagrams for the shells under the negative value of electric potential, without electric and magnetic potentials, and negative value of magnetic potential. It can be concluded from this figure that for a constant value of the amplitude parameter, the value of  is the highest in the presence of negative magnetic potential and is the smallest in the presence of negative electric potential. One can explain the reason for the mentioned behaviors of the diagrams of Figure 8 (b) as following: when the system is subjected to negative electric potential, its linear stress is higher than at . Conversely, when the system is under negative magnetic potential, its linear stress is smaller in comparison with the system without electric and magnetic potentials. It is apparent that the greater the linear stress, the weaker the effect of the nonlinearity. Because of this fact that the value of is much greater than the absolute value of , the effect of magnetic potential on the diagram is more apparent. In addition, one can justify the mentioned consequences of Figures 8 (a) and (b) using equation (30) as well as the values of  and  shown in Tables 5 and 6 for conditions considered in Figures 8 (a) and (b).

Table 7 represents the values of vibration characteristics including fundamental linear frequency as well as nonlinear parameters for different thickness values of the MEE composite conical shell. This table demonstrates that as increases, acquires smaller values; this is because of this fact that the increase of leads to an increase in the mass of the shell which leads to a decrease in the fundamental linear frequency. Another outcome of Table 7 is that the value of  decreases as the value of increases. Figures 9 (a) and (b) show the influence of the thickness of conical shells on the curves of in the situations with electric and magnetic potentials, respectively. One can figure out from these figures that the value of the nonlinear frequency ratio becomes greater with an increase in the thickness because of the decrease in the fundamental linear frequency as depicted in Table 7.

Figure 8. The influence of the electric and magnetic potentials on the nonlinear frequency ratio- amplitude parameter diagrams

Table 7. The values of fundamental linear frequency and nonlinear parameter for MEE composite conical shells
with different thickness values

Figure 9. Diagrams of nonlinear frequency ratio versus amplitude parameter for different thickness values for shells under, a: electric potential, b: magnetic potential

Table 8 represents the influence of the length on the fundamental linear frequency and nonlinear parameter ( ) of the conical shells under electric or magnetic potential. One can confirm from this table that the increase of the length leads to a decrease in the fundamental linear frequency and . The reason for the decrease of fundamental linear frequency is the increase of mass due to the increase in length. The effects of the length of the MEE composite conical shell on the nonlinear frequency ratio against amplitude parameter diagrams are shown in Figures 10 (a) and (b). It is apparent that Figures 10 (a) and (b) are for the shells subjected to electric and magnetic potentials, respectively. One can infer from Figures 10 (a) and (b) that for a determined amplitude parameter, the increase of the length makes smaller values for the nonlinear frequency ratio which is due to the decrease of the nonlinear parameter as depicted in Table 8.

Another result of Figures 8 to 10 is that the nonlinear frequency ratio becomes higher with an increase in the amplitude parameter.

Table 8. The effect of the length of the MEE composite conical shell on the vibration characteristics

Figure 10. The influence of the length of the MEE composite conical shell on the curves in the presence of,
 a: electric potential, b: magnetic potential

6.     Conclusions

This paper deals with the nonlinear vibration of an MEE composite conical shell surrounded by a nonlinear elastic foundation subjected to electric or magnetic potential. The relations of strains are extracted considering the effect of shear deformation with the help of the von Karman nonlinear approach.  Stress, electric displacement, and magnetic induction vectors are derived using coupled relations of MEE material.  Applying quasi-static approximation of Maxwell's vector equations, Gauss’ laws for electrostatics and magnetostatics, and considering the thin nature of the MEE composite conical shell lead to the extraction of electric and magnetic fields. The nonlinear ordinary differential equation of the system is extracted via the Lagrange technique while the effect of rotary inertia is considered in the extraction of kinetic energy. Lindstedt-Poincare method and modal analysis are employed to obtain the nonlinear responses of the MEE composite conical shell. The results of the literature are compared with this study's results to investigate the accuracy of the results of this research. The effects of several parameters including the nonlinear and linear constants of elastic foundation, the presence of electric or magnetic potential, thickness and length on the fundamental linear frequency, nonlinear parameter, and the curves of nonlinear frequency ratio versus amplitude parameter are investigated which can be classified as mentioned below:

1. Fundamental linear frequency does not have any change as the value of the nonlinear constant of the elastic foundation changes.
2. The value of fundamental linear frequency increases with an increase of the linear constants of the elastic foundation or the magnetic field and decreases with an increase in the thickness or the length.
3. Nonlinear parameter gets greater values with an increase of the nonlinear constant of the elastic foundation or the electric potential and decreases with an increase of the linear constants of elastic foundation, magnetic potential, thickness, or length.
4. For constant amplitude parameters, as the value of the nonlinear constant of the elastic foundation or thickness increases, the nonlinear frequency ratio becomes greater.
5. For a determined amplitude parameter, the value of the nonlinear frequency ratio becomes smaller with an increase of the linear constants of elastic foundation or the length.
6. The increase in the amplitude parameter leads to an increase of the nonlinear frequency ratio.

References

[1]   Duc, N. D., Anh, V. T. T. and Cong, P. H., 2014. Nonlinear axisymmetric response of FGM shallow spherical shells on elastic foundations under uniform external pressure and temperature. European Journal of Mechanics-A/Solids, 45, pp.80-89.

[2]   Duc, N. D., 2014. Nonlinear static and dynamic stability of functionally graded plates and shells. Vietnam: Vietnam Natl Univ Press. Google Scholar.

[3]   Duc, N. D., Cong, P. H., Tuan, N. D., Tran, P. and Van Thanh, N., 2017. Thermal and mechanical stability of functionally graded carbon nanotubes (FG CNT)-reinforced composite truncated conical shells surrounded by the elastic foundations. Thin-Walled Structures, 115, pp.300-310.

[4]   Duc, N. D., Seung-Eock, K. and Chan, D. Q., 2018. Thermal buckling analysis of FGM sandwich truncated conical shells reinforced by FGM stiffeners resting on elastic foundations using FSDT. Journal of Thermal Stresses, 41(3), pp.331-365.

[5]   Do, Q. C., Pham, D. N., Vu, D. Q., Vu, T. T. A. and Nguyen, D. D., 2019. Nonlinear buckling and post-buckling of functionally graded CNTs reinforced composite truncated conical shells subjected to axial load. Steel and Composite Structures, 31.

[6]   Quan, T. Q., Kim, S. E. and Duc, N. D., 2019. Nonlinear dynamic response and vibration of shear deformable piezoelectric functionally graded truncated conical panel in thermal environments. European Journal of Mechanics-A/Solids, 77, 103795.

[7]   Chan, D. Q., Long, V. D. and Duc, N. D., 2019. Nonlinear buckling and postbuckling of FGM shear-deformable truncated conical shells reinforced by FGM stiffeners. Mechanics of Composite Materials, 54(6), pp.745-764.

[8]   Manh, D. T., Anh, V. T. T., Nguyen, P. D. and Duc, N. D., 2020. Nonlinear post-buckling of CNTs reinforced sandwich-structured composite annular spherical shells. International Journal of Structural Stability and Dynamics, 20(02), 2050018.

[9]   Van Thanh, N., Khoa, N. D. and Duc, N. D., 2020. Nonlinear dynamic analysis of piezoelectric functionally graded porous truncated conical panel in thermal environments. Thin-Walled Structures, 154, 106837.

[10] Duc, N. D., Quan, T. Q. and Cong, P. H., 2021. Nonlinear vibration of auxetic plates and shells.

[11] Zhou, L., Nie, B., Ren, S., Żur, K. K. and Kim, J., 2021. On the hygro-thermo-electro-mechanical coupling effect on static and dynamic responses of piezoelectric beams. Composite Structures, 259, 113248.

[12] Jankowski, P., 2020. Effect of kerr foundation and in-plane forces on free vibration of fgm nanobeams with diverse distribution of porosity. Acta Mechanica et Automatica, 14(3), pp.135-143.

[13] Jankowski, P., Żur, K. K., Kim, J., Lim, C. W. and Reddy, J. N., 2021. On the piezoelectric effect on stability of symmetric FGM porous nanobeams. Composite Structures, 267, 113880.

[14] Żur, K. K., Farajpour, A., Lim, C. W. and Jankowski, P., 2021. On the nonlinear dynamics of porous composite nanobeams connected with fullerenes. Composite Structures, 274, 114356.

[15] Ghobadi, A., Golestanian, H., Beni, Y. T. and Żur, K. K., 2021. On the size-dependent nonlinear thermo-electro-mechanical free vibration analysis of functionally graded flexoelectric nano-plate. Communications in Nonlinear Science and Numerical Simulation, 95, 105585.

[16] Ouakad, H. M. and ŻUR, K. K., 2022. On the snap-through buckling analysis of electrostatic shallow arch micro-actuator via meshless Galerkin decomposition technique. Engineering Analysis with Boundary Elements, 134, pp.388-397.

[17] Razavi, S. and Shooshtari, A., 2014. Free vibration analysis of a magneto-electro-elastic doubly-curved shell resting on a Pasternak-type elastic foundation. Smart Materials and Structures, 23(10), 105003.

[18] Bhangale, R. K. and Ganesan, N., 2005. Free vibration studies of simply supported non-homogeneous functionally graded magneto-electro-elastic finite cylindrical shells. Journal of Sound and Vibration, 288(1-2), pp.412-422.

[19] Annigeri, A. R., Ganesan, N. and Swarnamani, S., 2006.  Free vibrations of simply supported layered and multiphase magneto-electro-elastic cylindrical shells.  Smart Materials and Structures, 15(2), 459.

[20] Tsai, Y. H. and Wu, C. P., 2008. Dynamic responses of functionally graded magneto-electro-elastic shells with open-circuit surface conditions. International journal of engineering science, 46(9), pp.843-857.

[21] Kumaravel, A., Ganesan, N. and Sethuraman, R., 2010. Buckling and vibration analysis of layered and multiphase magneto‐electro‐elastic cylinders subjected to uniform thermal loading. Multidiscipline Modeling in Materials and Structures.

[22] Lang, Z and Xuewu, L., 2013. Buckling and vibration analysis of functionally graded magneto-electro-thermo-elastic circular cylindrical shells.  Applied Mathematical Modelling, 37(4), pp.2279-2292.

[23] Shooshtari, A. and Razavi, S., 2016.  Large amplitude free vibration of magnetoelectroelastic curved panels. Scientia Iranica, 23(6), pp.2606-2615.

[24] Shooshtari, A., and Razavi, S., 2015. Linear and nonlinear free vibration of a multilayered magneto-electro-elastic doubly-curved shell on elastic foundation. Composites Part B: Engineering, 78, pp. 95-108.

[25] Mohammadimehr, M., Okhravi, S. V. and Akhavan Alavi, S. M., 2018. Free vibration analysis of magneto-electro-elastic cylindrical composite panel reinforced by various distributions of CNTs with considering open and closed circuits boundary conditions based on FSDT.  Journal of Vibration and Control, 24(8), pp.1551-1569.

[26] Vinyas, M and Harursampath, D, 2020. Nonlinear vibrations of magneto-electro-elastic doubly curved shells reinforced with carbon nanotubes. Composite Structures, 253, 112749.

[27] Rostami, R. and Mohammadimehr, M., 2020. Vibration control of rotating sandwich cylindrical shell-reinforced nanocomposite face sheet and porous core integrated with functionally graded magneto-electro-elastic layers. Engineering with Computers, pp.1-14.

[28] Ye, W., Liu, J., Zang, Q. and Lin, G., 2020. Magneto-electro-elastic semi-analytical models for free vibration and transient dynamic responses of composite cylindrical shell structures. Mechanics of Materials, 148, 103495.

[29] Rostami, R., Irani Rahaghi, M. and Mohammadimehr, M., 2021. Vibration control of the rotating sandwich cylindrical shell considering functionally graded core and functionally graded magneto-electro-elastic layers by using differential quadrature method. Journal of Sandwich Structures & Materials, 23(1), pp. 132-173.

[30] Ni, Y., Zhu, S., Sun, J., Tong, Z., Zhou, Z., Xu, X. and Lim, C. W., 2021. An accurate model for free vibration of porous magneto-electro-thermo-elastic functionally graded cylindrical shells subjected to multi-field coupled loadings. Journal of Intelligent Material Systems and Structures, 1045389X20986894.

[31] Zhou, L., Li, M., Tang, J., Li, F. and Żur, K. K., 2021. Evaluation of performance of magneto-electro-elastic sensor subjected to thermal-moisture coupled load via CS-FEM. Thin-Walled Structures, 169, 108370.

[32] Farajpour, A., Żur, K. K., Kim, J. and Reddy, J. N., 2021. Nonlinear frequency behaviour of magneto-electromechanical mass nanosensors using vibrating MEE nanoplates with multiple nanoparticles. Composite Structures, 260, 113458.

[33] Sofiyev, A. H., Kuruoglu, N. and Halilov, H. M., 2010. The vibration and stability of non-homogeneous orthotropic conical shells with clamped edges subjected to uniform external pressures. Applied mathematical modelling, 34(7), pp.1807-1822.

[34] Srikantamurthy, J. S., Annigeri, A. R. and Raghavendra, B. V., 2020. Free Vibration Analysis of Layered Magneto-Electro-Elastic Truncated Conical Shells. Journal of Applied Science and Engineering, 23(4), pp.669-676.

[35] Srikantamurthy, J. S. and Annigeri, A. R., 2021. Free Vibration Analysis of Multiphase Magneto-Electro-Elastic Composite Conical Shells. Indian Journal of Science and Technology, 14(19), pp. 1525-1533.

[36] Zarouni, E., Rad, M. J. and Tohidi, H., 2014. Free vibration analysis of fiber reinforced composite conical shells resting on Pasternak-type elastic foundation using Ritz and Galerkin methods. International Journal of Mechanics and Materials in Design, 10(4), pp.421-438.

[37] Sofiyev, A. H., 2014. The combined influences of heterogeneity and elastic foundations on the nonlinear vibration of orthotropic truncated conical shells. Composites Part B: Engineering, 61, pp.324-339.

[38] Asanjarani, A., Satouri, S., Alizadeh, A. and Kargarnovin, M. H., 2015. Free vibration analysis of 2D-FGM truncated conical shell resting on Winkler–Pasternak foundations based on FSDT. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 229(5), pp.818-839.

[39] Wu, S., Qu, Y. and Hua, H., 2015. Free vibration of laminated orthotropic conical shell on Pasternak foundation by a domain decomposition method. Journal of Composite Materials, 49(1), pp.35-52.

[40] Deniz, A., Zerin, Z. and Karaca, Z., 2016. Winkler-Pasternak foundation effect on the frequency parameter of FGM truncated conical shells in the framework of shear deformation theory. Composites Part B: Engineering, 104, pp.57-70.

[41] Nguyen Dinh, D. and Nguyen, P. D., 2017. The dynamic response and vibration of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) truncated conical shells resting on elastic foundations. Materials, 10(10), 1194.

[42] Nguyen, P. D., Quang, V. D., Anh, V. T. T. and Duc, N. D., 2019. Nonlinear vibration of carbon nanotube reinforced composite truncated conical shells in thermal environment. International Journal of Structural Stability and Dynamics, 19(12), 1950158.

[43] Chan, D. Q., Anh, V. T. T. and Duc, N. D., 2019. Vibration and nonlinear dynamic response of eccentrically stiffened functionally graded composite truncated conical shells surrounded by an elastic medium in thermal environments. Acta Mechanica, 230(1), pp.157-178.

[44] Sofiyev, A. H., Mammadov, Z., Dimitri, R. and Tornabene, F., 2020. Vibration analysis of shear deformable carbon nanotubes‐based functionally graded conical shells resting on elastic foundations. Mathematical Methods in the Applied Sciences.

[45] Eyvazian, A., Musharavati, F., Tarlochan, F., Pasharavesh, A., Rajak, D. K., Husain, M. B. and Tran, T. N., 2020. Free vibration of FG-GPLRC conical panel on elastic foundation. Structural Engineering and Mechanics, 75(1), pp.1-18.

[46] Mohammadrezazadeh, S. and Jafari, A. A., 2021. Nonlinear vibration suppression of laminated composite conical shells on elastic foundations with magnetostrictive layers. Composite Structures, 258, 113323.

[47] Fu, T., Wu, X., Xiao, Z., Chen, Z. and Li, B., 2021. Analysis of vibration characteristics of FGM sandwich joined conical–conical shells surrounded by elastic foundations. Thin-Walled Structures, 165, 107979.

[48] Jamalabadi, M. Y. A., Borji, P., Habibi, M. and Pelalak, R., 2021. Nonlinear vibration analysis of functionally graded GPL-RC conical panels resting on elastic medium. Thin-Walled Structures, 160, 107370.

[49] Amirabadi, H., Farhatnia, F. and Civalek, Ӧ., 2021. Frequency response of rotating two-directional functionally graded GPL-reinforced conical shells on elastic foundation. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 43(7), pp.1-24.

[50] Mohammadrezazadeh, S. and Jafari, A. A., 2021. The study of the nonlinear vibration of SS and CC laminated composite conical shells on elastic foundations through an approximate method. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 09544062211021439.

[51] Sofiyev, A. H. and Kuruoglu, N., 2017. Combined effects of transverse shear stresses and nonlinear elastic foundations on the nonlinear dynamic response of heterogeneous orthotropic cylindrical shells. Composite Structures, 166, pp.153-162.

[52] Zhu, C., Fang, X., Liu, J. and Nie, G., 2020. Smart control of large amplitude vibration of porous piezoelectric conical sandwich panels resting on nonlinear elastic foundation. Composite Structures, 246, 112384.

[53] Molla-Alipour, M., Shariyat, M. and Shaban, M., 2020. Free Vibration Analysis of Bidirectional Functionally Graded Conical/Cylindrical Shells and Annular Plates on Nonlinear Elastic Foundations, Based on a Unified Differential Transform Analytical Formulation. Journal of Solid Mechanics, 12(2), pp.385-400.

[54] Civalek, Ö. and Kiracioglu, O., 2010. Free vibration analysis of Timoshenko beams by DSC method. International Journal for Numerical Methods in Biomedical Engineering, 26(12), pp.1890-1898.

[55] Civalek, Ö. M. E. R., 2014. Geometrically nonlinear dynamic and static analysis of shallow spherical shell resting on two-parameters elastic foundations. International Journal of Pressure Vessels and Piping, 113, pp.1-9.

[56] Mercan, K., Demir, Ç. and Civalek, Ö., 2016. Vibration analysis of FG cylindrical shells with power-law index using discrete singular convolution technique. Curved and Layered Structures, 3(1).

[57] Pirbodaghi, T., Ahmadian, M. T. and Fesanghary, M., 2009. On the homotopy analysis method for non-linear vibration of beams. Mechanics Research Communications, 36(2), pp.143-148.

[58] Motallebi, A. A., Poorjamshidian, M. and Sheikhi, J., 2014. Vibration analysis of a nonlinear beam under axial force by homotopy analysis method. Journal of Solid Mechanics, 6(3), pp.289-298.

[59] Dastjerdi, S., Akgöz, B. and Civalek, Ö., 2020. On the effect of viscoelasticity on behavior of gyroscopes. International Journal of Engineering Science, 149, 103236.

[60] Lam, K. Y. and Hua, L., 1997. Vibration analysis of a rotating truncated circular conical shell. International Journal of Solids and Structures, 34(17), pp.2183-2197.

[61] Rao, S. S., 2007. Vibration of Continuous Systems. JOHN WILEY & SONS, INC.

[62] Reddy, J. N., 2004. Mechanics of laminated composite plates and shells: theory and analysis. Second Edition, CRC Press.

[63] Li, Y. and Zhang, J., 2013. Free vibration analysis of magnetoelectroelastic plate resting on a Pasternak foundation. Smart materials and structures, 23(2), 025002.

[64] Pan, E. and Heyliger, P. R., 2002. Free vibrations of simply supported and multilayered magneto-electro-elastic plates. Journal of Sound and Vibration, 252(3), pp.429-442.

[65] Harshe, G., Dougherty, J. P. and Newnham, R. E., 1993. Theoretical modelling of multilayer magnetoelectric composites. International journal of applied electromagnetics in materials, 4(2), pp.145-145.

[66] Nan, C. W., 1994.  Magnetoelectric effect in composites of piezoelectric and piezomagnetic phases. Physical Review B, 50(9), 6082.

[67] Benveniste, Y., 1995. Magnetoelectric effect in fibrous composites with piezoelectric and piezomagnetic phases. Physical Review B, 51(22), 16424.

[68] Pan, E., 2001. Exact solution for simply supported and multilayered magneto-electro-elastic plates. Journal of Applied Mechanics-TRANSACTIONS OF THE ASME, 68(4), pp.608-618.

[69] Sladek, J., Sladek, V., Krahulec, S. and Pan, E., 2013. Analyses of functionally graded plates with a magnetoelectroelastic layer. Smart materials and structures, 22(3), 035003.

[70] Liu, M. F. and Chang, T. P., 2010. Closed form expression for the vibration problem of a transversely isotropic magneto-electro-elastic plate. Journal of Applied Mechanics- TRANSACTIONS OF THE ASME, 77(2).

[71] Milazzo, A., 2014. Large deflection of magneto-electro-elastic laminated plates. Applied Mathematical Modelling, 38(5-6), pp.1737-1752.

[72] Qatu, M. S., 2004. Vibration of laminated shells and plates. Elsevier.

[73] Nayfeh, A. H. and Mook, D. T., 1995. Nonlinear Oscillations. John Wiley & Sons, Inc.

[74] Irie, T., Yamada, G. and Tanaka, K., 1984. Natural frequencies of truncated conical shells. Journal of Sound and Vibration, 92(3), pp.447-453.

[75] Li, F. M., Kishimoto, K. and Huang, W. H., 2009. The calculations of natural frequencies and forced vibration responses of conical shell using the Rayleigh–Ritz method. Mechanics Research Communications, 36(5), pp.595-602.

[76] Shah, A. G., Mahmood, T., Naeem, M. N. and Arshad, S. H., 2011. Vibration characteristics of fluid-filled cylindrical shells based on elastic foundations. Acta mechanica, 216(1), pp.17-28.

[77] Ramirez, F., Heyliger, P. R. and Pan, E., 2006. Discrete layer solution to free vibrations of functionally graded magneto-electro-elastic plates. Mechanics of Advanced Materials and Structures, 13(3), pp.249-266.