Efficient Analytical Model for Time-Dependent Behavior of Layered Functionally Graded Plates with Viscoelastic Interlayers
Layered functionally graded (FG) plate plays an important role in engineering constructions. In this work, an efficient analytical solution is proposed to investigate the time-dependent characters of layered FG plate with viscoelastic interlayers, in which each FG layer is described by the elasticity theory, and the interlayer viscoelasticity is simulated by the generalized Maxwell model. The constitutive equations in the interlayer are simplified, and then, the analytical solutions of stresses and displacements for the layered FG plate are solved by virtue of the efficient recursive matrix method. Some examples are analyzed to investigate the influences of geometric and material characteristics on the long-term behaviors for the FG plate.
The layered plates are commonly used in various engineering [1–4], by right of their outstanding merits, such as high strength, lightweight, and corrosion resistance. Enough bonding stiffness between adjacent members is the guarantee of mechanical performance of layered plates. In many cases, the connection between adjacent layers is not rigid, and thus, the interfacial slip often happens [5, 6]. Besides, the bond behavior in layered plates exhibits viscoelastic property, due to the use of polymer adhesive [7, 8]. As a result, the mechanical performance of the whole layered plates is time-dependent under sustained loads .
Although the layered plates exhibit good mechanical performance, it is necessary to reduce immediately the modulus difference between adjacent layers, which is a major cause of large interlaminar stress . This requirement has led to the development of a novel material, called functionally graded (FG) materials, which can be designed to have continuous mechanical properties [11–13]. Due to its outstanding behaviors, FG materials have been applied in various areas, such as piezoelectric fiber-reinforced composite , heat-resistant materials , and nonuniform porous materials . The concept of FG materials has also been introduced into the layered plates as face and/or transition layers, which can optimize the stress distribution and deformation [17, 18].
A number of works have been proposed to study the mechanical responses of layered FG plates. Gunes et al.  presented a backpropagation artificial neural network method for investigating the three-dimensional free vibration response of an adhesively bonded wide and narrow FG plates. Demirbas and Apalak  performed a thermoelastic investigation of FG circular plates with adhesive bond and gave the thermal stress and strain distributions with different gradient indexes. Based on the equivalent single-layer higher-order theory, Tornabene et al.  analyzed the dynamic response of anisotropic doubly curved shells with arbitrary geometry and variable thickness. According to the four-variable refined plate theory, an analytical solution was deduced by Li et al.  for the thermomechanical bending analysis of FG sandwich plate. Wang et al.  derived an analytical solution for investigating the heat transfer behavior in FG Sandwich plates. Tornabene et al. [24, 25] used the higher-order theory to study the dynamic behavior of anisotropic doubly curved shells. An extended cohesive damage model was proposed by Ghimire and Chen  for the analysis of geometrical ratio influences on the failure mechanisms of FG sandwiches, and they found that the layered core can provide significant improvements in loading capacity. On the basis of the four-variable plate theory, Trinh el al.  presented a Levy solution with state space concept for static, dynamic, and buckling analyses for sandwich FG plates. Khorshidi and Karimi  established an analytical model based on a modified shear deformation theory for flutter investigation of sandwich plates with FG face sheets under thermal condition. A full layerwise method was employed by Nikbakht et al.  to obtain the yielding initiation of FG sandwich plates subjected to bi-sinusoidal distributed loading under general boundary conditions. Based on the first-order shear deformation theory, Tornabene  investigated the dynamic behavior of moderately thick functionally graded conical, cylindrical shells, and annular plates.
In the above literature, most works in regard to layered plates were based on the assumption of perfectly bonding or static slip, while the time-dependent characters resulting from viscoelastic bonding interlayers were neglected. Furthermore, the solutions based on the simplified plate theories have considerable error for thick plates.
This work proposes an efficient analytical model for investigating the time-dependent bending behavior of layered FG plates with considering viscoelastic bonding interlayer. In the analytical model, the mechanical behavior of each FG layer is expressed by the three-dimensional elasticity theory, which renounces the shear deformation assumption. Thus, it is highly accurate for structures with any thickness. The interlayer viscoelasticity is described by the generalized Maxwell model. The constitutive equations in the interlayer are simplified, and then, the analytical solutions of stresses and displacements for the layered FG plate are solved by virtue of the efficient recursive matrix method. Additionally, some examples provided by the present solution are conducted to investigate the influences of geometric and material characteristics on the time-dependent behaviors of the FG plate.
2. Efficient Analytical Model
In Figure 1, a layered FG plate formed by p FG layers bonded by viscoelastic bonding interlayer is considered, which is simply supported and bears to a sustained load q(x,y) on its top. The three-dimensional size of the plate is a × b × H, and the thicknesses for FG layer and bonding interlayer are, respectively, hi and Δh, in which the label i means the layer index. The elastic modulus of each FG layer, denoted by , varies according to the exponential law along the thickness direction, and is expressed bywhere is the graded index and . The generalized Maxwell model, as shown in Figure 2, is employed to model the interlayer viscoelasticity, with time-dependent shear modulus given by
2.1. Governing Equations for a FG Layer
According to the elasticity theory, the constitutive equations for the ith (i = 1, 2, …, p) FG layer arewhere , , , , , and are stresses, , , , , , and are strains and and are lame parameters expressed as follows:
The strain-displacement relations are given bywhere , , and represent the displacements in x, y, and z direction, respectively. The stress components of each FG layer should meet the following static equilibrium equations:
The simply supported boundary conditions can be expressed by
Four steps are taken to obtain the differential equation involving displacement components. Firstly, by substituting (1) and (5) into (3) and eliminating the strain components, one has
By substituting (8) into (6) and eliminating the stress components, the governing equations only involving the displacement components are obtained:
The present structural problem is displacement-based, and the stresses can be determined by (8) after the displacements are solved. The above governing equations are actually partial differential equations which can hard to be solved directly. Secondly, for simply supported boundary conditions, the displacements of each FG layer can be expanded in Fourier series as follows:where and . Thirdly, by applying the Fourier series expansion of (10) to (9), the partial differential equations become ordinary differential equations as follows:
At last, a fourth-order differential equation of is derived out by eliminating and in (11):
The general solution of iswhere are the undetermined time-dependent coefficients, which can be obtained according to the boundary and continuity conditions in the subsequent sections; the details of are defined in Appendix A. By substituting into (13) and (11), and are obtained:where the details of are listed in Appendix A. By substituting of (13), (14), and (8) into (9), the stress components involving undetermined coefficients can be written as
2.2. Governing Equations for a Viscoelastic Interlayer
According to the Boltzmann superposition principle, the constitutive equations for the ith (i = 1, 2, …, p-1) interlayer are given by
The above convolution integral equations means that the stress of viscoelastic interlayer depends on the total strain history, which leads to the heavy calculation and time-consuming. In the present work, they are simplified as
In comparison with the exact solution, this simplification leads to small error, and the present solution is always on the side of safety . Considering that Δh is far less than hi, the shear strains can be assumed to be constant in z direction and given by
The shear stress continuity relationships between adjacent layers are
2.3. Recursive Matrix Method
The stresses on the top and bottom surfaces of the plate are
Since the stresses and displacements are in series form, q(x,y) should also be expanded as
By substitution of (13)–(15) into (8), the general solution of stresses and displacements can be transferred into the matrix form aswherewhere the nonzero coefficients in matrix can be found in Appendix B. By combining (17)–(19), the continuity relationships between the adjacent FG layers are rewritten in the matrix form aswhere
By substituting and into the z-coordinate in (24), respectively, we obtain
Elimination of in (26) yields
By reusing (24) and (27), from i = 1 to i = p, one obtains
Four 3 × 3, submatrices are used to define the matrix multiplication in (28) as follows:
Thus, (28) can be written as
By the decomposition of (30), two submatrix equations are given as
By solving (31), the displacement boundary value on the bottom surface of the structure can be expressed as
Similar to the relationship in (28), for any FG layer is obtained from (24) and (27) as follows:
The time-dependent coefficients of the ith FG layer are further obtained as
Finally, the substitution of (34) into (13)–(15) yields the solution of stress and displacement components for each FG layer.
3. Numerical Results and Discussion
In the following, the FG sandwich plate (p = 3) is taken as an example. The variables with the superscript r, e.g., , are defined as the result of stress or displacement at the point of , , and . The variables with two vertical lines denote their absolute values, e.g., .
3.1. Validation of the Present Solution
The convergence property of the present solution is assessed first. The series are truncated into a finite number M for actual calculations here. The parameters in the FG plate are taken as q(x,y) = 1 N/mm2, a = 1000 mm, b = 800 mm, ∆h = 0.5 mm, h1 = h3 = 30 mm, h2 = 40 mm, = = 80 GPa, = 40 GPa, k1 = k3 = 0.05, k2 = 0, and = = = = 0.3. The material of polyvinyl butyral (PVB) is chosen for the viscoelastic interlayer and its viscoelastic parameters of which are given in Table 1. Table 2 lists the present results with different series terms, respectively. It is found that the present results tend to be constant as series terms increase, and the convergence accuracy is four significant digits.
The present solution is compared with the Kirchhoff–Love (KL) solution , respectively. The parameters in the FG plate are fixed at q(x,y) = sin(πx/a)sin(πy/b) N/mm2, a = b = 1000 mm, ∆h = 0.2 mm, k2 = 0, = = 80 GPa, = 40 GPa, = = = = 0.3, t = 104 s, = 0.5 MPa, = 500 MPa, = 1 s, h1:h2:h3 = 1 : 2:1, and k3 = −k1. The above parameters mean that the geometric and material properties of the present structure are symmetric about the midplane in the direction of thickness. Figure 3 represents the relative error between KL solution and the present one for a different length-to-height ratio a/H and graded index k1, respectively. A good agreement is found for KL solution in thin plate case with small |k1|, while the relative error of KL results increases gradually as the plate become thick or k1 decreases. The errors for , , and are 15.4%, 6.6%, and 15.9%, respectively, as a/H = 8 and k1 = −0.09. The imprecision of KL solutions mainly results from the neglect of the transverse shear deformation and the deviation of neutral plane of the FG layer.
3.2. Parameter Research
In this section, some examples provided by the present solution are conducted to investigate the influences of geometric and material characteristics on the time-dependent behaviors of the FG plate. Some of the parameters are fixed at q(x,y) = sin(πx/a)sin(πy/b) N/mm2, a = b = 1000 mm, k1:k3 = −1, k2 = 0, = , = 25Gpa, and = = = = 0.3; the viscoelastic constants are given in Table 1, while the others are variable.
Figure 4 plots the stress and displacement distribution in z direction when t = 1 day, 1 year, and 10 years and the corresponding results in perfectly bonded (PB) case. The constants are fixed at h1 = h3 = 20 mm, h2 = 40 mm, ∆h = 0.2 mm, and k1 = −0.05. It can be obtained from Figure 4 that , , , , and considerably increase with t, while changes slightly with time. and show zig-zag distributions, and they, in the facial layers, are obviously curve distribution, which is different from isotropic material. also show a zig-zag distribution but keeps straight in each FG layer. gives a multipeak distribution. Compared with the results of PB case, the maximum values of , , , , and increase by 200.1%, 26.2%, 200.2%, 604.7%, and 1221%, respectively, at t = 10 years.
The effect of the graded index on the stress and displacement of the FG plate, as well as the elastic modulus distributions along the thickness direction, are illustrated in Figure 5. The parameters are taken as h1 = h3 = 10 mm, h2 = 30 mm, ∆h = 0.2 mm, and t = 1 day, and the average modulus of FG facial layer, i.e., (i = 1, 3), is fixed at 80 GPa. From the results shown in Figures 5(b) and 5(c), the decline of k1 reduces the maximum value of and but enlarges that of . As k1 goes down from 0.2 to -0.2, and fall by 10.0% and 15.7%, respectively, while increases by 156.7%. In the other word, the stress and displacement distribution can be optimized by adjusting the graded index.
Figure 6 shows the influences of t, k1, and Δh on the maximum value of the stress and displacement. The parameters are fixed at the same as those in Figure 5, except t, k1 and Δh are variable. It can be seen from Figure 6 that, as k1 = −0.1, , , and increase monotonously and keeps invariant as t or Δh increases. Due to the increase of t and Δh resulting in a reduced shear modulus of the interlayer, and as t⟶∞, approaches to the fixed value . Similarly, and increase monotonously and tend to definite values with the increase of k1 or t as Δh = 0.2 mm, which results from that the difference in elastic modulus between the adjacent layer surfaces grows gradually as k1 increases. For a given t, decreases first and then increases as k1 grows, which is owing to the fact that the location of transfers from the outer surface of the facial layer to the inner surface in the early stage, while in the medium or long term, the location of transfers from the facial layer to the core layer as k1 increases. Similarly, for a small certain k1, increases monotonously with t, while for a large k1, firstly decreases and then increases as t goes on, which is due to the transfer of the location from the facial layer to the core layer.
The elasticity theory with the recursive matrix method was proposed to analyze the time-dependent behavior of the layered functionally graded plates with viscoelastic interlayer. The following conclusions can be provided:(1)In the thin plate case, the KL solution with small graded index approaches the present one. Since the transverse shear deformation is neglected and the neutral plane gradually deviates from the middle plane, the relative errors increase gradually as length-to-height ratio or graded index increases.(2)As time goes on, the influence of the interlayer shear modulus degeneration for the adhesive bonding case is obvious, which leads to the reduced interface shear stress and the increased interfacial slip as well as the greater deflection.(3)By adjusting the graded index, the stress and displacement distribution of the case with fixed material consumption can be optimized and the location of the maximum value of the normal stress transfers between the facial layer and the core layer. The maximum value of deflection and shear stress decreases with the reduction of the graded index and the interlayer thickness.
A. Details of Coefficients in General Solutions
In (13) and (14), the details of (j = 1–6) and (l = 1–4) are given as follows:where
B. Details of Elements in Coefficient Matrix
in (22) is with the following nonzero elements:
The raw/processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study.
Conflicts of Interest
The authors declare no conflicts of interest with respect to the research, authorship, and/or publication of this article.
This research was financially supported by National Natural Science Foundation of China (no. 52108220) and Natural Science Foundation of Jiangsu Province (no. BK20190668).
F. Moleiro, E. Carrera, A. J. M. Ferreira, and J. N. Reddy, “Hygro-thermo-mechanical modelling and analysis of multilayered plates with embedded functionally graded material layers,” Composite Structures, vol. 233, Article ID 111442, 2020.View at: Publisher Site | Google Scholar
V. Gulizzi, I. Benedetti, and A. Milazzo, “A high-resolution layer-wise discontinuous Galerkin formulation for multilayered composite plates,” Composite Structures, vol. 242, Article ID 112137, 2020.View at: Publisher Site | Google Scholar
L. V. Tran and S. E. Kim, “Static and free vibration analyses of multilayered plates by a higher-order shear and normal deformation theory and isogeometric analysis,” Thin-Walled Structures, vol. 130, pp. 622–640, 2018.View at: Publisher Site | Google Scholar
I. Benedetti and A. Milazzo, “Advanced models for smart multilayered plates based on Reissner Mixed Variational Theorem,” Composites Part B: Engineering, vol. 119, pp. 215–229, 2017.View at: Publisher Site | Google Scholar
L. Galuppi and G. Royer-Carfagni, “Laminated beams with viscoelastic interlayer,” International Journal of Solids and Structures, vol. 49, no. 18, pp. 2637–2645, 2012.View at: Publisher Site | Google Scholar
P. Foraboschi, “Analytical model for laminated-glass plate,” Composites Part B: Engineering, vol. 43, no. 5, pp. 2094–2106, 2012.View at: Publisher Site | Google Scholar
A. Tounsi, M. S. A. Houari, S. Benyoucef, and E. A. A. Bedia, “A refined trigonometric shear deformation theory for thermoelastic bending of functionally graded sandwich plates,” Aerospace Science and Technology, vol. 24, no. 1, pp. 209–220, 2013.View at: Publisher Site | Google Scholar
R. Xu and Y. F. Wu, “Two-dimensional analytical solutions of simply supported composite beams with interlayer slips,” International Journal of Solids and Structures, vol. 44, no. 1, pp. 165–175, 2007.View at: Publisher Site | Google Scholar
A. Zemanová, J. Zeman, and M. Šejnoha, “Comparison of viscoelastic finite element models for laminated glass beams,” International Journal of Mechanical Sciences, vol. 131-132, pp. 380–395, 2017.View at: Publisher Site | Google Scholar
M. K. Apalak and R. Gunes, “Thermal residual stresses in an adhesively-bonded functionally graded single-lap joint,” Journal of Adhesion Science and Technology, vol. 20, no. 12, pp. 1295–1320, 2006.View at: Publisher Site | Google Scholar
R. M. J. Groh and P. M. Weaver, “On displacement-based and mixed-variational equivalent single layer theories for modelling highly heterogeneous laminated beams,” International Journal of Solids and Structures, vol. 59, pp. 147–170, 2015.View at: Publisher Site | Google Scholar
K. Susanto, “Vibration analysis of piezoelectric laminated slightly curved beams using distributed transfer function method,” International Journal of Solids and Structures, vol. 46, no. 6, pp. 1564–1573, 2009.View at: Publisher Site | Google Scholar
N. T. Khiem, H. T. Tran, and D. Nam, “Modal analysis of cracked continuous Timoshenko beam made of functionally graded material,” Mechanics Based Design of Structures and Machines, vol. 48, no. 4, pp. 459–479, 2020.View at: Publisher Site | Google Scholar
A. M. Zenkour and R. A. Alghanmi, “Bending of exponentially graded plates integrated with piezoelectric fiber-reinforced composite actuators resting on elastic foundations,” European Journal of Mechanics - A: Solids, vol. 75, pp. 461–471, 2019.View at: Publisher Site | Google Scholar
S. K. Jalali, M. H. Naei, and A. Poorsolhjouy, “Thermal stability analysis of circular functionally graded sandwich plates of variable thickness using pseudo-spectral method,” Materials & Design, vol. 31, no. 10, pp. 4755–4763, 2010.View at: Publisher Site | Google Scholar
Y. Zhang, G. Jin, M. Chen, T. Ye, C. Yang, and Y. Yin, “Free vibration and damping analysis of porous functionally graded sandwich plates with a viscoelastic core,” Composite Structures, vol. 244, Article ID 112298, 2020.View at: Publisher Site | Google Scholar
M. H. G. Babaki and M. Shakouri, “Free and Forced Vibration of sandwich Plates with Electrorheological Core and Functionally Graded Face Layers,” Mechanics Based Design Structures and Machines, vol. 49, no. 5, pp. 689–706, 2019.View at: Publisher Site | Google Scholar
M. Dorduncu, “Stress analysis of sandwich plates with functionally graded cores using peridynamic differential operator and refined zigzag theory,” Thin-Walled Structures, vol. 146, Article ID 106468, 2020.View at: Publisher Site | Google Scholar
R. Gunes, M. A. Kemal, and I. Ozkes, “Free vibration analysis of adhesively bonded single lap joints with wide and narrow functionally graded plates,” Composite Structures, vol. 92, no. 1, pp. 1–17, 2010.View at: Publisher Site | Google Scholar
M. D. Demirbas and M. K. Apalak, “Investigation of the thermo-elastic response of adhesively bonded two-dimensional functionally graded circular plates based on theory of elasticity,” Iranian Journal of Science and Technology, Transactions of Mechanical Engineering, vol. 42, no. 4, pp. 415–433, 2018.View at: Publisher Site | Google Scholar
F. Tornabene, M. Viscoti, and M. Dinitri, “Equivalent single layer higher order theory based on a weak formulation for the dynamic analysis of anisotropic doubly-curved shells with arbitrary geometry and variable thickness,” Thin-Walled Structures, vol. 174, Article ID 109119, 2022.View at: Publisher Site | Google Scholar
D. D. Li, Z. B. Deng, G. P. Chen, H. Xiao, and L. Zhu, “Thermomechanical bending analysis of sandwich plates with both functionally graded face sheets and functionally graded core,” Composite Structures, vol. 169, pp. 29–41, 2017.View at: Publisher Site | Google Scholar
X. Wang, Z. Wang, T. Zeng, S. Cheng, and F. Yang, “Exact analytical solution for steady-state heat transfer in functionally graded sandwich slabs with convective-radiative boundary conditions,” Composite Structures, vol. 192, pp. 379–386, 2018.View at: Publisher Site | Google Scholar
F. Tornabene, M. Viscoti, M. Dinitri, and J. N. Reddy, “Higher order theories for the vibration study of doubly-curved anisotropic shells with a variable thickness and isogeometric mapped geometry,” Composite Structures, vol. 267, Article ID 113829, 2021.View at: Publisher Site | Google Scholar
F. Tornabene, M. Viscoti, and R. Dimitri, “Generalized higher order layerwise theory for the dynamic study of anisotropic doubly -curved shells with a mapped geometry,” Engineering Analysis with Boundary Elements, vol. 134, pp. 147–183, 2022.View at: Publisher Site | Google Scholar
S. Ghimire and J. Chen, “An extended cohesive damage model study of geometrical ratio effects on failure mechanisms of functionally graded sandwiches with multi-layered cores,” Composite Structures, vol. 224, Article ID 110999, 2019.View at: Publisher Site | Google Scholar
L. C. Trinh, T. P. Vo, H. T. Thai, T. K. Nguyen, and P. Keerthan, “State-space Levy solution for size-dependent static, free vibration and buckling behaviours of functionally graded sandwich plates,” Composites Part B: Engineering, vol. 149, pp. 144–164, 2018.View at: Publisher Site | Google Scholar
K. Khorshidi and M. Karimi, “Flutter analysis of sandwich plates with functionally graded face sheets in thermal environment,” Aerospace Science and Technology, vol. 95, Article ID 105461, 2019.View at: Publisher Site | Google Scholar
S. Nikbakht, S. J. Salami, and M. Shakeri, “A 3D full layer-wise method for yield achievement in Functionally Graded Sandwich Plates with general boundary conditions,” European Journal of Mechanics - A: Solids, vol. 75, pp. 330–347, 2019.View at: Publisher Site | Google Scholar
F. Tornabene, “Free vibration analysis of functionally graded conical, cylindrical shell and annular plate structures with a four-parameter power-law distribution,” Computer Methods in Applied Mechanics and Engineering, vol. 198, pp. 37–40, 2009.View at: Publisher Site | Google Scholar