Home Dynamic response of functionally graded plate under harmonic load with variable gradient parameters
Article Open Access

Dynamic response of functionally graded plate under harmonic load with variable gradient parameters

  • Jianhui Tian EMAIL logo , Hongrui Zhang , Jinjuan Sun , Jialiang Wu and Guangchu Hu
Published/Copyright: July 15, 2022
Download Article (PDF)
Become an author with De Gruyter Brill
Submit Manuscript Author Information
Science and Engineering of Composite Materials
From the journal Science and Engineering of Composite Materials Volume 29 Issue 1

Abstract

This article used the strip element method to study the dynamics problems of the functionally graded plate with variable gradient parameters under the harmonic loads. The dynamic model of the functionally graded plate is established by using the strip element method, the rationality and accuracy of the theoretical results are verified by finite element method, and the displacement response under different gradient parameters is also calculated. The results show that under the different gradient parameters, the displacement varies harmonically with time, and with the increase of gradient parameters, the fluctuation period of displacement with time increases continuously, and the displacement peak also gradually increases. The displacement along the thickness direction also shows the harmonic form. Through comparison, it is found that the gradient parameters have a greater impact on the dynamic response for the functionally graded plate; with the increase in the gradient parameters, the displacement response also increases, but the displacement response trend slows down.

Keywords: functionally graded plate; strip element method; harmonic load; dynamic response

1 Introduction

Due to the high specific strength and specific stiffness, and easy design, functionally graded materials (FGMs) can work under some extreme environmental conditions and are widely used in aerospace, biomedicine, petrochemical, heavy industry, and other fields [1,2,3]. With the development of FGMs, the research on the dynamic characteristics of FGMs becomes more and more important. In the process of dynamics research, at present, in the field of traditional FGMs research, linear parameters are generally used, which results in certain limitations in practical applications. Therefore, the dynamics of FGMs with variable gradient parameters is extremely critical and important.

Compared with other composite materials, FGMs have better mechanical properties. The key researches mainly focus on vibration, shock, and dynamic response under different external conditions and different boundary conditions by using different methods, such as theoretical solutions, numerical calculations, or experiments. Shariyat et al. [4] used finite element method and iterative algorithm to study the dynamic stress, displacement distribution, hygrothermal elastic wave propagation, and reflection response of a functionally graded hollow sphere under thermomechanical shock. Hao et al. [5] studied the nonlinear forced vibration and natural frequency of a double-bent flat shell with a rectangular sandwich FGM and verified its effectiveness using numerical simulation methods. Gupta and Talha [6] focused on the structural characteristics of FGMs plates and shells under different boundaries and environmental conditions under thermo-electromechanical loads. Bakhtiari and Kheradpisheh [7] studied the transient response of an inflatable multilayer hollow functionally graded cylinder with interlayer adhesion defects under load. This structure is widely used in aerospace structures. Bozyigit et al. [8] applied the theory of univariate shear deformation to the analysis of the free vibration and harmonious response of the multi-layer frame model considering the basic flexibility. Aris and Ahmadi [9] studied the nonlinear vibration and resonance analysis of truncated gradient tapered shells under harmonic excitation. Parandvar and Farid [10] established a nonlinear finite element model of the dynamic response of a FGM plate under thermal, static, and harmonic loads and studied the effect of initial conditions and static pressure on the dynamic response of the system. Najarzadeh et al. [11,12] studied the free vibration of thin plates under arbitrary load using the boundary element method. Parida et al. [13] used the finite element method to perform dynamic analysis on the simply supported beam structure and studied the effect of temperature as a function of dynamic parameters on the mechanical properties. Guo et al. [14] analyzed the modes of simply supported plates with uniform thickness and stepped thickness by using a dynamic shape function. Using a dynamic stiffness method (DSM), Kumar and Jana [15] studied the free vibration characteristics of rectangular FGM thin plates with S-FGM and E-FGM characteristics along the thickness direction, and they studied the free vibration characteristics of functionally gradient rectangular plate using DSM [16].

For functionally graded plates with variable gradient parameters, previous research work has involved the study of the thermal conduction of functionally graded plates with variable gradient parameters under heat source load [17], but the research on the dynamic performance of variable gradient parameters has not been carried out, so based on the strip element method, fixed boundary conditions are applied to study the dynamic response of functionally graded plates with variable gradient parameters. In addition, Karamia et al. [18] studied the influence of various boundary conditions and found that for all boundaries. Dhital et al. [19] used the properties of two materials to smoothly transition from one material to another for reducing thermal stress, residual stress, and stress concentration factor. These studies are very helpful for the development of the research on functionally graded plates with variable gradient parameters. And Liu et al. [20] verified the correctness of the strip element method. However, the study of functionally graded plates with variable gradient parameters under harmonic load has not been mentioned. In this article, the variable gradient parameter model will be established on the basis of the strip element method to calculate the displacement response of the functionally gradient plate under different gradient parameters and verify its rationality and accuracy, this work provides a new way for the variable gradient parameter model.

2 Establishment of dynamic model

2.1 Dynamic theory model

When studying the dynamic performance of a functionally graded plate, as shown in Figure 1, the domain is defined by D g = (x, z), in which x ∈ (−∞, +∞), z ∈ (0, H), and D p ∈ D g , D p are composed of boundaries B 1, B 2, B 3, and B 4. Since the functionally graded plate is isotropic in the x–y plane and the y-direction is infinite, simplifying the problem to the two-dimension relation with x and z, there are only two displacement components u and v in x-direction and z-direction, respectively.

Figure 1 Functionally graded plate discrete model.
Figure 1

Functionally graded plate discrete model.

Without considering the internal damping of the functionally graded plate, it is assumed that the dynamic control equation of the functionally graded plate in matrix form as

(1) M U ̈ − L T σ = Q ,

where M is the mass matrix, U is the displacement vector, L T is the differential operator matrix, σ is the stress vector, and Q is the load vector.

2.2 System equation

The relationship between strain and displacement in matrix form is

(2) ε = L U ,

where ε is the vector of the strain, defined as

(3) ε = ε x x ε z z γ x z .

The vector form of the displacement is

(4) U = u v ,

where u and v are the x-direction and z-direction displacements, respectively.

The differential operator matrix is

(5) L = ∂ ∂ x 0 0 ∂ ∂ z ∂ ∂ z ∂ ∂ x .

Simplify the differential operator matrix to

(6) L = L x ∂ ∂ x + L z ∂ ∂ z ,

where L x , L z is a constant matrix

(7) L x = 1 0 0 0 0 1 ,

(8) L z = 0 0 0 1 1 0 .

The relationship between stress and strain in matrix form is

(9) σ = c ε ,

where the vector form of stress is

(10) σ = σ x x σ z z τ x z .

The stiffness coefficient matrix is

(11) c = c 11 c 12 c 13 c 21 c 22 c 23 c 31 c 32 c 33 .

Transforming the equation (1) could obtain the motion balance equation in the form of displacement

(12) M U ̈ − L T c L U = Q ,

where

(13) L T c L = D x x ∂ 2 ∂ x 2 + 2 D x z ∂ 2 ∂ x ∂ z + D z z ∂ 2 ∂ z 2 ,

and D xx , D xz , and D zz could be calculated

(14) D x x = L x T c L x = c 11 c 13 c 13 c 33 ,

(15) D x z = 1 2 [ L x T c L z + L z T c L z ] = 1 2 c 13 c 33 + c 12 c 33 + c 12 2 c 23 ,

(16) D z z = L z T c L z = c 33 c 23 c 23 c 22 .

Assuming boundary stress form as

(17) R = R x R z ,

when x is a constant, the stress acting on the x plane can be expressed as

(18) R x = L x T c L U = L x T c L x ∂ U ∂ x + L x T c L z ∂ U ∂ z = D x x ∂ U ∂ x + D x z ′ ∂ U ∂ z ,

when z is a constant, the stress acting on the z-plane can be expressed as

(19) R z = L z T c L U = L z T c L x ∂ U ∂ x + L z T c L z ∂ U ∂ z = D z x ′ ∂ U ∂ x + D z z ∂ U ∂ z ,

where

(20) D x z ′ = L x T c L z = c 13 c 12 c 33 c 23 ,

(21) D z x ′ = L z T c L x = c 13 c 33 c 12 c 23 = [ D x z ′ ] T .

3 Dynamic theory of variable gradient parameter strip element method

3.1 Establishment of the theoretical model of variable gradient parameters

In the functionally graded plate, the volume fraction of ZrO2 is assumed to be

(22) V C = 1 − a 1 − z H + b 1 − z H r p ,

where a, b, r, and p is the gradient parameters and z is the coordinate thickness.

In this calculation model, gradient parameters are a = 1, b = 0, and gradient parameter p is a variable parameter from 0 to positive infinity. The volume fraction of ceramics in the model changes with the gradient value. When the gradient parameter p approaches infinite, the volume fraction of ceramics is close to zero, and when the gradient parameter p approaches zero, the volume fraction of ceramics is close to 1. Density, Poisson’s ratio and the elastic modulus of the functionally graded plate all depend on the volume fraction of ceramics. Different gradient parameters correspond to the volume fraction of ceramics in the direction of thickness as shown in Figure 2.

Figure 2 Volume fraction of ceramics under different gradient parameters.
Figure 2

Volume fraction of ceramics under different gradient parameters.

The functionally graded plate is composed of ZrO2 and Ti–6Al–4V, and its material properties are shown in Table 1.

Table 1

Material properties of ZrO2 and Ti–6Al–4V

Material Density Elastic modulus Poisson ratio
ρ (kg m−3) E(Pa)
ZrO2 5,331 1.1 × 1011 0.25
Ti–6Al–4V 4,420 3.5 × 1011 0.32

As the volume fraction in the functionally graded plate changes, the density of the model also changes accordingly. The density in the model is assumed to be

(23) ρ ( z ) = ρ M ( 1 − V C ) + ρ C V C ,

where the ρ C is the density of ZrO2, ρ M is the density of Ti–6Al–4V.

Similarly, the elastic modulus of the functionally graded plate changes as

(24) E ( z ) = E M ( 1 − V C ) + E C V C ,

where E C is the elastic modulus of ZrO2 and E M is the elastic modulus of Ti–6Al–4V.

The variation of density corresponding to different gradient parameters in the thickness direction is shown in Figure 3. The change law of elastic modulus is consistent with the change law of the mass density.

Figure 3 Density of functionally graded plates varies with gradient parameters.
Figure 3

Density of functionally graded plates varies with gradient parameters.

3.2 Variable gradient parameter strip element method theory

Divide the domain D g into N strip elements, assuming that the displacement field of the element is

(25) U ( x , z ) = N ( z ) V ( x ) exp ( − i ω t ) ,

where V(x) is the displacement vector, N(z) is the shape function matrix, t is the time, and ω is the angular frequency.

(26) N ( z ) = ( 1 − 3 z + 2 z 2 ) I 4 ( z − z 2 ) I ( − z + 2 z 2 ) I ,

(27) V = V U T V M T V L T T ,

where the I is the identity matrix, which V U T , V M T , V L T are the displacements of the upper, middle, and lower nodal lines of the elements.

Assume the form of the applied harmonic external load is

(28) Q = Q ¯ exp ( − i ω t ) ,

where Q ¯ is the amplitude vector.

According to the principle of virtual work, we could obtain

(29) δ V T ⋅ Q = δ V T ⋅ R + ∫ 0 h δ U T ( M U ̈ − L T c L U − Q ) d z .

Solve the approximate differential set of a set of elements

(30) Q ¯ = − A 2 ∂ 2 V ∂ z 2 + A 1 ∂ V ∂ z + A 0 V − ω 2 M V ,

where

(31) A 0 = 1 3 h 7 D x x − 8 D x x D x x − 8 D x x 16 D x x − 8 D x x D x x − 8 D x x 7 D x x ,

(32) A 1 = 1 3 3 ( D x z − D x z ′ ) − 4 D x z D x z 4 D x z 0 − 4 D x z − D x z 4 D x z − 3 ( D x z − D x z ′ ) ,

(33) A 2 = h 30 4 D z z 2 D z z − D z z 2 D z z 16 D z z 2 D z z − D z z 2 D z z 4 D z z ,

(34) M = ∫ 0 h ρ N T ( z ) N ( z ) d x = ρ h 30 4 I 2 I − I 2 I 16 I 2 I − I 2 I 4 I ,

in the mass matrix, I is a 2 × 2 identity matrix.

The obtained element equations are extended to the domain D g to obtain the approximate differential equations

(35) Q ¯ g = − A 2 g ∂ 2 V g ∂ x 2 + A 1 g ∂ V g ∂ x + A 0 g V g − ω 2 M g V g .

The subscript g represents the matrix or vector in the whole domain, and the matrix A i g ( i = 0 , 1 , 2 ) , M g and vector Q g , V g are formed by combining the corresponding matrices and vectors of all elements. A ig and M g are M × M matrices, where M = 4N + 2 and N is the number of elements.

Equation (35) represents a set of second-order differential equations with constant coefficients that can be accurately solved. Assuming the form of displacement and load is

(36) V g = d g exp ( i k x ) ,

(37) Q ¯ g = P g exp ( i k x ) ,

where k is the wave number.

Substituting equations (36) and (37) into equation (35) could obtain

(38) P g = [ k 2 A 2 g + i k A 1 g + A 0 g − ω 2 M g ] d g ,

where P g is the amplitude vector of the external force acting on the boundary node line.

(39) P g = 0 P .

For the given P g , the characteristic equation could be obtained:

(40) 0 = [ k 2 A 2 g + i k A 1 g + A 0 g − ω 2 M g ] d g − P g .

For a given ω, equation (39) can be transformed into a standard characteristic equation:

(41) 0 = 0 I ω 2 M g − A 0 g − i A 1 g − k I 0 0 A 2 g d g k d g − P g .

Solving equation (41) can get 2m(m = 6N − 2(N − 1) = 4N + 2) eigenvalues k j (1…2m), if the upper part of the jth eigenvectors corresponding to d g is represented by ϕ j , then there is

(42) φ j T = φ 1 j φ 2 j φ 3 j .

According to the mode superposition method, the displacement of the x-direction can be written as the superposition of all eigenvectors

(43) V g = ∑ j = 1 2 m C j φ j exp ( i k j x ) = G ( x ) C ,

where C j is the undetermined coefficient, and equation (43) is a basic solution system in the D g domain. The constant vector C can be determined from the boundary B 1, B 2 to obtain a special solution of the problem domain D p . There are m nodes and 2m boundary conditions on the boundary B 1, B 2, respectively, resulting in 2m constants C j .

According to equation (43), the constant vector C can be expressed as the displacement vector of the node on the boundary B 1, B 2

(44) C = G b − 1 V b g = G B 1 G B 2 − 1 V b g ,

where

(45) V b g = V b g B 1 V b g B 2 = u 1 v 1 u 2 v 2 ⋅ ⋅ ⋅ u m v m ,

(46) G B 1 = φ 11 Z 1 B 1 φ 12 Z 2 B 1 ⋯ φ 1 L Z B 4 B 1 φ 21 Z 1 B 1 φ 22 Z 2 B 1 ⋯ φ 2 L Z B 4 B 1 ⋮ ⋮ ⋱ ⋮ φ m 1 Z 1 B 1 φ m 2 Z 2 B 1 ⋯ φ m L Z B 4 B 1 ; ,

(47) Z j B 1 = exp ( i k j z B 1 ) ,

where B 1 is the left boundary, and B 2 is the right boundary. In equation (45), u i and v i represent the displacements in the z-direction and the x-direction at the point i on the boundaries B 1 and B 2, respectively. In equation (47), z B 1 is the coordinate in the z direction on the left boundary B 1. The matrix G B 2 in equation (44) is the same as matrix G B 1 but Z j B 1 is replaced by Z j B 2 = exp ( i k j z j B 2 ) , where z B2 is the coordinate in the z-direction on the boundary B 2.

3.3 Application of boundary conditions

Using equation (18), at any point of x is a constant, the stress vector in one element is

(48) R x = D x z ' ∂ U ∂ z + D x x ∂ U ∂ x ,

where

(49) D x z ' = c 13 c 23 c 12 c 23 ,

and

(50) R = R x | z = 0 R x | z = h / 2 R x | z = h T = R 1 V + R 2 ∂ V ∂ x ,

where

(51) R 1 = 1 h − 3 D z x ' 4 D z x ' − D z x ' − D z x ' 0 D z x ' D z x ' − 4 D z x ' 3 D z x ' ,

(52) R 2 = D x x 0 0 0 D x x 0 0 0 D x x .

Assemble matrix R 1, R 2 into the overall unit matrix and get the stress vector of the nodal line as

(53) R g = R 1 g V g + R 2 g ∂ V g ∂ x .

Simultaneous equations (43)–(43) could obtain

(54) R b g = K V b g ,

where

(55) R b g = R b g L R b g R .

R b t L and R b t R respectively represent the outer vectors of the left and right boundaries.

(56) K = R 1 g 0 0 R 1 g + R 2 g ∂ G U ∂ x G b − 1 R 2 g ∂ G L ∂ x G b − 1 ,

where K is the stiffness matrix. Equation (54) gives the relationship between stress and displacement at the left and right boundaries. For a given R bg , V bg could be calculated by equation (54), and V g can be calculated from equations (43) and (44). The global displacement in the time domain can be solved by equation (25).

3.4 Application of the dynamic strip method under simple harmonic load

The thickness of the functionally graded plate along the z direction is H = 90, and it is divided into 10-layer elements. Within the boundary formed by B 1, B 2, B 3, and B 4, the length of the functionally graded plate along the x-direction is x = H = 90. Harmonic load action time t = 2 s. The surface materials of the first layer and the tenth layer of this material are ZrO2 and Ti–6Al–4V, respectively, and the properties of the intermediate material change continuously along the thickness direction.

The amplitude of the applied harmonic load Q ¯ = 4.9 × 10 4 , by the load equation (45), and the loading form of the harmonic load varying with time are shown in Figure 4.

Figure 4 Harmonic load forms.
Figure 4

Harmonic load forms.

4 Results and discussion

Before the output of the theoretical results, the rationality and accuracy of the theoretical results need to be verified. For the theoretical in this article, further finite element simulation is needed to verify its accuracy and correctness. The study of functionally graded plates with variable gradient parameters shows that when the gradient parameters tend to infinity, the volume fraction of metal materials in functionally graded plates tends to 1. Therefore, the displacement response of pure metal material under harmonic loading is analyzed by the finite element method, and the results are in agreement with the strip element method (Figure 5).

Figure 5 Verification of strip element method results and finite element results.
Figure 5

Verification of strip element method results and finite element results.

According to the displacement form of the strip element method, the law of the overall displacement of the functionally graded plate with time at z = 0 in the lower surface, gradient parameters p = 0.5, p = 1, p = 2, p = 5, and p = 10 are calculated as shown in Figure 6 under the harmonic load. It can be obtained that the overall displacement of the functionally gradient plate under the action of harmonic load presents sinusoidal fluctuation with time. Under the different gradient parameters, the displacement is in the form of simple harmonic motion with time, and with the increase of gradient parameters, the fluctuation period of displacement with time increases continuously, and the displacement peak also gradually increases, because the metal composition of the functionally gradient plate increases with the increase of gradient parameters, and the plasticity of metal is greater than that of ceramic.

Figure 6 Overall displacement changes with time under different gradient parameters at z = 0.
Figure 6

Overall displacement changes with time under different gradient parameters at z = 0.

Figure 7 shows the variation law of displacement with gradient parameters at t = 0.2 s and t = 1.3 s. It can be found that the displacement of the functionally graded plate increases with the gradient parameters, but the increasing trend first increases and then decreases.

Figure 7 Overall displacement changes with the gradient parameter.
Figure 7

Overall displacement changes with the gradient parameter.

Figure 8(a)–(c) respectively shows at the gradient parameters p 1 = 0.5, p 2 = 1, p 3 = 2, p 4 = 5, p 5 = 10 of the functionally graded plate as the displacement changes along the thickness at x = 0 mm, x = 45 mm, and x = 90 mm. It can be obtained that the displacements along the thickness direction at different positions of x are in harmonic form within the boundary. The peak displacement appears between x = 20 and x = 80. The displacement peak value is different under different gradient parameters. As the gradient parameter increases, the displacement peak value gradually increases and the displacement phase also increases.

Figure 8 Displacement changes under different gradient parameters along the thickness with different x positions. (a) x = 0, (b) x = 45, and (c) x = 90.
Figure 8

Displacement changes under different gradient parameters along the thickness with different x positions. (a) x = 0, (b) x = 45, and (c) x = 90.

Figure 9(a)–(c) respectively shows the displacement of the upper surface, middle surface, and lower surface of the functionally graded plate with gradient parameters at x = 10 mm, x = 45 mm, and x = 90 mm. It can be obtained that the displacement of the upper surface, the middle surface, and the lower surface is related to the gradient parameters and increases with the increase in the gradient parameters. With the increase in the gradient parameters, the volume fraction of metal in the functionally graded plate increases, and the plastic deformation capacity of metals is higher than that of ceramics. Therefore, as the gradient parameters increase, the displacement increases gradually.

Figure 9 The upper, middle, and lower surface displacement changes with gradient parameters with different x positions. (a) x = 0, (b) x = 45, and (c) x = 90.
Figure 9

The upper, middle, and lower surface displacement changes with gradient parameters with different x positions. (a) x = 0, (b) x = 45, and (c) x = 90.

Figure 10(a) and (b), respectively, shows the variation of displacement along the x-direction on the upper, middle, and lower surfaces of the functionally graded plate at p = 0.5 and p = 2. It can be seen that the displacement of the functionally graded plate has a faster downward trend along the x-direction and gradually tends to zero on the boundary. Among them, the displacement when the displacement p = 0.5 is smaller than that when p = 2, and the displacement is maximum on the middle surface.

Figure 10 The upper, middle, and lower surface displacement changes along the x-direction with different parameters. (a) p = 0.5 and (b) p = 2.
Figure 10

The upper, middle, and lower surface displacement changes along the x-direction with different parameters. (a) p = 0.5 and (b) p = 2.

5 Conclusion

The strip element method is first used to study the dynamics of functionally graded plates with variable gradient parameters under harmonic loads, and useful conclusions are obtained:

  1. Under the different gradient parameters, the displacement varies harmonically with time, and with the increase of gradient parameters, the fluctuation period of displacement with time increases continuously, and the displacement peak also gradually increases.

  2. Under the fixed boundary conditions, the displacement of the functionally graded plate along the thickness direction shows the harmonic form.

  3. By comparing the displacement at different positions for gradient parameters, as the gradient parameters increase, the peak displacement of the functionally graded plate increases, but the trend of displacement changes decreases.

  4. The displacement perpendicular to the load direction gradually tends to zero along the x-axis, and the displacement peak appears on the middle surface of the functionally graded plate.

Acknowledgments

The research is supported by the National Natural Science Foundation of China (No. 11302159), the Principal Foundation of Xi’an Technological University (No. XGPY200213), Basic Research Plan of Natural Science in Shannxi Province (No. 2021JQ-650), and Scientific Research Program Funded by Shaanxi Provincial Education Department (No. 19JK0412).

  1. Conflict of interest: The authors state no conflict of interest.

References

[1] Fujita K, Furukawa H, Yamanaka M, Niino M. Solar excited laser using FGM concept for space energy network. Mater Sci Forum. 2003;423–425:807–12.10.4028/www.scientific.net/MSF.423-425.807Search in Google Scholar

[2] Obata Y, Takeuchi K, Kawazoe M, Kanayama K. Design of functionally graded wood-based board for floor heating system with higher energy efficiency. Mater Sci Forum. 2003;423–425:819–24.10.4028/www.scientific.net/MSF.423-425.819Search in Google Scholar

[3] Uemura S. The activities of FGM on new application. Mater Sci Forum. 2003;423–425:1–10.10.4028/www.scientific.net/MSF.423-425.1Search in Google Scholar

[4] Shariyat M, Jahanshahi S, Rahimi H. Nonlinear hermitian generalized hygrothermoelastic stress and wave propagation analyses of thick FGM spheres exhibiting temperature, moisture, and strain-rate material dependencies. Composite Struct. 2019;229:111364.10.1016/j.compstruct.2019.111364Search in Google Scholar

[5] Hao Y, Li Z, Zhang W, Li S, Yao M. Vibration of functionally graded sandwich doubly curved shells using improved shear deformation theory. Sci China Technol Sci. 2018;61(6):791–808.10.1007/s11431-016-9097-7Search in Google Scholar

[6] Gupta A, Talha M. Recent development in modeling and analysis of functionally graded materials and structures. Prog Aerosp Sci. 2015;79:1–14.10.1016/j.paerosci.2015.07.001Search in Google Scholar

[7] Bakhtiari M, Kheradpisheh M. Dynamics behavior and imperfection sensitivity of a Fluid-filled multilayered FGM cylindrical structure. Int J Mech Sci. 2020;176:105425.10.1016/j.ijmecsci.2020.105425Search in Google Scholar

[8] Bozyigit B, Yesilce Y, Wahab MA. Single variable shear deformation theory for free vibration and harmonic response of frames on flexible foundation. Eng Struct. 2020;208:110268.10.1016/j.engstruct.2020.110268Search in Google Scholar

[9] Aris H, Ahmadi H. Nonlinear vibration analysis of FGM truncated conical shells subjected to harmonic excitation in thermal environment. Mech Res Commun. 2020;104:103499.10.1016/j.mechrescom.2020.103499Search in Google Scholar

[10] Parandvar H, Farid M. Large amplitude vibration of FGM plates in thermal environment subjected to simultaneously static pressure and harmonic force using multimodal FEM. Composite Struct. 2016;141:163–71.10.1016/j.compstruct.2016.01.044Search in Google Scholar

[11] Najarzadeh L, Movahedian B, Azhari M. Free vibration and buckling analysis of thin plates subjected to high gradients stresses using the combination of finite strip and boundary element methods. Thin-Walled Struct. 2018;123:36–47.10.1016/j.tws.2017.11.015Search in Google Scholar

[12] Najarzadeh L, Movahedian B, Azhari M. Stability analysis of the thin plates with arbitrary shapes subjected to non-uniform stress fields using boundary element and radial integration methods. Eng Anal Bound Elem. 2018;87:111–21.10.1016/j.enganabound.2017.11.010Search in Google Scholar

[13] Parida SP, Jena PC, Dash. RR. FGM beam analysis in dynamical and thermal surroundings using finite element method. Mater Today. 2019;18:3676–82.10.1016/j.matpr.2019.07.301Search in Google Scholar

[14] Guo SJ, Keane AJ, Moshrefi-Torbati M. Vibration analysis of stepped thickness plates. J Sound Vib. 1997;204(4):645–57.10.1006/jsvi.1997.0955Search in Google Scholar

[15] Kumar S, Jana P. Application of dynamic stiffness method for accurate free vibration analysis of sigmoid and exponential functionally graded rectangular plates. Int J Mech Sci. 2019;163:105105.10.1016/j.ijmecsci.2019.105105Search in Google Scholar

[16] Kumar S, Ranjan V, Jana P. Free vibration analysis of thin functionally graded rectangular plates using the dynamic stiffness method. Composite Struct. 2018;197:39–53.10.1016/j.compstruct.2018.04.085Search in Google Scholar

[17] Tian J, Zhang H, Sun J. Study on heat conduction of functionally graded plate with the variable gradient parameters under the H(t) heat source. J Therm Stresses. 2020;43(6):724–38.10.1080/01495739.2020.1742623Search in Google Scholar

[18] Karamia B, Janghorban M, Rabczuk T. Dynamics of two-dimensional functionally graded tapered timoshenko nanobeam in thermal environment using nonlocal strain gradient theory. Compos Part B. 2020;182:107622.10.1016/j.compositesb.2019.107622Search in Google Scholar

[19] Dhital S, Rokaya A, Kaizer MR, Zhang Y, Kim J. Accurate and efficient thermal stress analyses of functionally graded solids using incompatible graded finite elements. Composite Struct. 2019;222:110909.10.1016/j.compstruct.2019.110909Search in Google Scholar PubMed PubMed Central

[20] Liu G, Xi Z, Horie Y. Elastic waves in an-isotropic laminates. Appl Mech Rev. 2001;56(2):B23.10.1115/1.1553437Search in Google Scholar
Received: 2021-10-06
Revised: 2022-06-12
Accepted: 2022-06-23
Published Online: 2022-07-15

© 2022 Jianhui Tian et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Experimental investigations of a novel pressure microfoam preparation device for dust removal
  3. Influence of hydrothermal aging on the mechanical performance of foam core sandwich panels subjected to low-velocity impact
  4. Experimental study on surface wrapping strengthening of EPS particles and its concrete performance
  5. Modification of mechanical properties of Shanghai clayey soil with expanded polystyrene
  6. A new EPS beads strengthening technology and its influences on axial compressive properties of concrete
  7. A novel superabsorbent material based on soybean straw: Synthesis and characterization
  8. Use of line laser scanning thermography for the defect detection and evaluation of composite material
  9. Research on back analysis of meso-parameters of hydraulic cemented sand and gravel based on Box-Behnken design response surface
  10. Hot deformation behavior and microstructure of a 0.5 wt% graphene nanoplatelet reinforced aluminum composite
  11. Analysis of electromagnetic characteristics of the proposed composite four-rail electromagnetic launcher
  12. Preparation and characterization of a graphene hybridizing polyurethane damping composite
  13. Effects of layup parameters and interference value on the performance of CFRP–metal interference fit joints
  14. Vibration and noise reduction of pipelines using shape memory alloy
  15. Finite element analysis of behavior and ultimate strength of composite column
  16. Dynamic response of functionally graded plate under harmonic load with variable gradient parameters
  17. Deformation behavior of rubber composite based on FEA and experimental verification
  18. Effects of Z-pin on moisture absorption property and damage mode under flexural load for carbon fiber composite
  19. Design and testing of a smart rubber stave for marine water-lubricated bearings
  20. Study of carbon nano-modifier of fly ash in cement concrete mixtures of civil engineering
  21. Analysis of multiple impact tests’ damage to three-dimensional four-directional braided composites
  22. Theoretical analysis of aluminum honeycomb sandwich panel supported by reinforced concrete wall under low-speed impact load
  23. Effects of local fiber discontinuity on the fatigue strength parameter at the fiber inclusion corner in fiber-reinforced composites
  24. Experimental investigation on compressive properties of three-dimensional five-directional braided composites in hygrothermal environment
  25. Failure process of steel–polypropylene hybrid fiber-reinforced concrete based on numerical simulations
  26. A simple method for measuring the monofilament diameter of continuous filament yarn with high bending stiffness via synthetic laser imaging
  27. Span length effect on flexural properties of composite laminate reinforced with a plain weave carbon fiber fabric in a polymer matrix
  28. Mechanical properties improving and microstructure characterization of inorganic artificial stone binder
  29. Effect of thermal treatment process on the structure of C/SiO2 composite aerogels
  30. Mechanical and corrosion resistance analysis of laser cladding layer
  31. Wear and corrosion mechanisms of Ni–WC coatings modified with different Y2O3 by laser cladding on AISI 4145H steel
  32. Damage and failure analysis of composite stiffened panels under low-velocity impact and compression after impact with damp-heat aging
  33. In-situ CT characterization of 2D woven SiCf/SiC composite loading under compression
  34. Effect of the manufacturing process on the equivalency qualification of glass fiber reinforced polymer
  35. Study of concrete properties based on crushed stone sand mixture and fiber of fly ash of thermal power plants
  36. Establishment of wear mechanism distribution diagram of ZTAp-reinforced iron matrix composites
  37. Calculation method of elastic modulus for carbon fiber-reinforced plastics considering inhomogeneous interphase
  38. An experimental study on the failure and enhancement mechanism of bolt-strengthening GFRP T-joint subjected to tensile loading
  39. The viability of cell that encapsulated in calcium alginate hydrogel beads
  40. Discussion of ceramic bar reinforced TWIP steel composite structure
  41. A theoretical framework underlying an accelerated testing method and its application to composites under constant strain rates and fatigue loading
  42. Theoretical analysis of interfacial design and thermal conductivity in graphite flakes/Al composites with various interfacial coatings
  43. Multiscale heat conduction and fractal oxidation behaviors of needle-punched carbon/carbon composites
  44. Numerical simulation of composite grid sandwich structure under low-velocity impact
  45. Wear properties of Al/TiO2 composites fabricated via combined compo-casting and APB process
  46. Review Articles
  47. Application of melanin as biological functional material in composite film field
  48. Review on research progress of cemented sand and gravel dam
  49. Communication
  50. Fabrications and microstructure analysis of cobalt-based coatings by an easy-coating and sintering process
  51. Letter to the Editor
  52. Investigation on mechanical and conductive behaviors of nano-graphite-based concrete
https://doi.org/10.1515/secm-2022-0019
Keywords for this article
functionally graded plate; strip element method; harmonic load; dynamic response
Creative Commons
BY 4.0

Articles in the same Issue

  1. Regular Articles
  2. Experimental investigations of a novel pressure microfoam preparation device for dust removal
  3. Influence of hydrothermal aging on the mechanical performance of foam core sandwich panels subjected to low-velocity impact
  4. Experimental study on surface wrapping strengthening of EPS particles and its concrete performance
  5. Modification of mechanical properties of Shanghai clayey soil with expanded polystyrene
  6. A new EPS beads strengthening technology and its influences on axial compressive properties of concrete
  7. A novel superabsorbent material based on soybean straw: Synthesis and characterization
  8. Use of line laser scanning thermography for the defect detection and evaluation of composite material
  9. Research on back analysis of meso-parameters of hydraulic cemented sand and gravel based on Box-Behnken design response surface
  10. Hot deformation behavior and microstructure of a 0.5 wt% graphene nanoplatelet reinforced aluminum composite
  11. Analysis of electromagnetic characteristics of the proposed composite four-rail electromagnetic launcher
  12. Preparation and characterization of a graphene hybridizing polyurethane damping composite
  13. Effects of layup parameters and interference value on the performance of CFRP–metal interference fit joints
  14. Vibration and noise reduction of pipelines using shape memory alloy
  15. Finite element analysis of behavior and ultimate strength of composite column
  16. Dynamic response of functionally graded plate under harmonic load with variable gradient parameters
  17. Deformation behavior of rubber composite based on FEA and experimental verification
  18. Effects of Z-pin on moisture absorption property and damage mode under flexural load for carbon fiber composite
  19. Design and testing of a smart rubber stave for marine water-lubricated bearings
  20. Study of carbon nano-modifier of fly ash in cement concrete mixtures of civil engineering
  21. Analysis of multiple impact tests’ damage to three-dimensional four-directional braided composites
  22. Theoretical analysis of aluminum honeycomb sandwich panel supported by reinforced concrete wall under low-speed impact load
  23. Effects of local fiber discontinuity on the fatigue strength parameter at the fiber inclusion corner in fiber-reinforced composites
  24. Experimental investigation on compressive properties of three-dimensional five-directional braided composites in hygrothermal environment
  25. Failure process of steel–polypropylene hybrid fiber-reinforced concrete based on numerical simulations
  26. A simple method for measuring the monofilament diameter of continuous filament yarn with high bending stiffness via synthetic laser imaging
  27. Span length effect on flexural properties of composite laminate reinforced with a plain weave carbon fiber fabric in a polymer matrix
  28. Mechanical properties improving and microstructure characterization of inorganic artificial stone binder
  29. Effect of thermal treatment process on the structure of C/SiO2 composite aerogels
  30. Mechanical and corrosion resistance analysis of laser cladding layer
  31. Wear and corrosion mechanisms of Ni–WC coatings modified with different Y2O3 by laser cladding on AISI 4145H steel
  32. Damage and failure analysis of composite stiffened panels under low-velocity impact and compression after impact with damp-heat aging
  33. In-situ CT characterization of 2D woven SiCf/SiC composite loading under compression
  34. Effect of the manufacturing process on the equivalency qualification of glass fiber reinforced polymer
  35. Study of concrete properties based on crushed stone sand mixture and fiber of fly ash of thermal power plants
  36. Establishment of wear mechanism distribution diagram of ZTAp-reinforced iron matrix composites
  37. Calculation method of elastic modulus for carbon fiber-reinforced plastics considering inhomogeneous interphase
  38. An experimental study on the failure and enhancement mechanism of bolt-strengthening GFRP T-joint subjected to tensile loading
  39. The viability of cell that encapsulated in calcium alginate hydrogel beads
  40. Discussion of ceramic bar reinforced TWIP steel composite structure
  41. A theoretical framework underlying an accelerated testing method and its application to composites under constant strain rates and fatigue loading
  42. Theoretical analysis of interfacial design and thermal conductivity in graphite flakes/Al composites with various interfacial coatings
  43. Multiscale heat conduction and fractal oxidation behaviors of needle-punched carbon/carbon composites
  44. Numerical simulation of composite grid sandwich structure under low-velocity impact
  45. Wear properties of Al/TiO2 composites fabricated via combined compo-casting and APB process
  46. Review Articles
  47. Application of melanin as biological functional material in composite film field
  48. Review on research progress of cemented sand and gravel dam
  49. Communication
  50. Fabrications and microstructure analysis of cobalt-based coatings by an easy-coating and sintering process
  51. Letter to the Editor
  52. Investigation on mechanical and conductive behaviors of nano-graphite-based concrete
Sign up now to receive a 20% welcome discount
Institutional Access
How does access work?
Have an idea on how to improve our website?
Please write us.
© 2025 De Gruyter Brill
Downloaded on 28.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/secm-2022-0019/html
Scroll to top button