Elasticity Solutions for Sandwich Arches considering Permeation Effect of Adhesive

Huo, Ruili; Liu, Yichen; Wu, Peng; Fang, Hai; Liu, Weiqing; Zhou, Ding

doi:https://doi.org/10.1155/2020/7358930

Advances in Polymer Technology

On this page

Abstract Introduction Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Applications of Fiber Reinforced Polymer Composites in Civil Engineering

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 7358930 | https://doi.org/10.1155/2020/7358930

Elasticity Solutions for Sandwich Arches considering Permeation Effect of Adhesive

Ruili Huo,¹Yichen Liu,¹Peng Wu,¹Hai Fang,¹Weiqing Liu,¹and Ding Zhou¹

Guest Editor: Guangming Chen

Received21 Oct 2019

Accepted07 Mar 2020

Published28 Mar 2020

Abstract

In this work, analytical solution of simply supported sandwich arches considering permeation effect of adhesives is presented. The permeation layer is described by the functionally graded material, exponentially graded in the radial direction. The stresses and deformations of each layer are based on the two-dimensional (2D) elasticity theory in the polar coordinate. The governing equations of the arch are solved by the layer-wise method, which turns the differential equations with variable coefficients into constant coefficients. The solution can be obtained efficiently by means of the recursive matrix method, especially for the arch with many layers. The present solution agrees well with the finite element solution with a fine mesh, while the finite element method is time consuming in mesh division and calculation. The one-dimensional (1D) solution based on the Euler–Bernoulli theory is close to the present one; however, the error increases as the arch becomes thick. The effect of permeation layer thickness on the stresses is studied. It is indicated that the stress distributions tend to be smooth along the radial direction as the permeation layer thickness increases.

1. Introduction

Owing to various advantages such as corrosion resistance, antifatigue, and high specific strength and stiffness, composite structures are increasingly used in different branches of engineering. In addition, their mechanical behavior can be optimized by ingeniously choosing the material type and composite pattern. A typical application is the sandwich structures [1–5], usually composed of two stiff face layers and a relatively soft core layer, bonded together by the adhesive. It is known that the interfacial damage exists widely in sandwich structures due to the discontinuity of material property at the interface. In practice, the adhesive will permeate into adjacent macrovoid materials such as glass fiber, balsa wood, and aluminum honeycomb, which leads to a refined sandwich model, as shown in Figure 1. The permeation layer is the mixture of the face (or core) layer and the adhesive. This makes the material property smoothly vary along the radial direction, like of functionally graded material (FGM). Such a problem deserves to be deeply investigated. For sandwich plates considering permeation effect, analytical solution has been proposed by Huo et al. [6]. In addition to sandwich plates, sandwich arches have also received a lot of attention due to their advantage of artistic appearance and excellent load-carrying capacity. The present study aims at extending their work to investigate the sandwich arch with permeation effect.

A review of the literature indicates that many analytical models have been proposed for the mechanical behaviors of sandwich structures. Galuppi and Royer-Carfagni [7] studied the bending behavior of the sandwich arch formed by two elastic facial layers and a thin polymeric interlayer. In their report, the deformation of the arch was based on the Euler–Bernoulli assumption, which was suitable for thin arches with moderate curvature. The stress distribution in a sandwich plate with the FG face was studied by Raissi et al. [8] by the use of the layer-wise method based on the first-order shear deformation theory. The vibration of sandwich microplates with nanocomposite face layers under magnetic and electric fields was studied by Mohammadimehr et al. [9] by virtue of the sinusoidal shear deformation theory. Livani et al. [10] investigated the supersonic flutter of doubly curved sandwich panels with variable thickness subjected to thermal load. In their study, the governing equations were based on a new higher-order shear deformation theory, in which the vertical and in-plane displacement components were assumed as quadratic and cubic functions, respectively. Bouchafa et al. [11] presented a new refined hyperbolic shear deformation theory with only four unknown functions to analyze the thermoelastic bending of FG sandwich plates. Bending, vibration, and buckling analysis of FG sandwich plates were analyzed by Nguyen et al. [12] according to the refined higher-order shear deformation theory. By using a higher-order sandwich theory, Rahmani et al. [13] tackled the free vibration of the composite sandwich shell with a flexible core.

In addition to the aforementioned analytical models, some numerical methods are efficient for the analysis of sandwich structures. By means of the differential quadrature finite element method, Liu et al. [14] investigated bending of the sandwich shell consisting of two homogenous face sheets and a FG core. Natarajan and Manickam [15] presented a QUAD-8 shear flexible element method based on higher-order shear deformation theory to study the bending and free flexural vibration of FG sandwich plates. Dynamic behavior of corrugated-core sandwich plates was analyzed by Peng et al. [16] via a mesh-free Galerkin method, in which the domain of the plate was discretized by scattered nodes, and no elements were required. By virtue of the isogeometric method, Thai et al. [17] studied the free vibration and buckling analysis of sandwich plates. This method has the advantage to exactly represent domains and achieve approximation with arbitrarily high smoothness. Mantari et al. [18] developed a layer-wise finite element formulation-based trigonometric layer-wise shear deformation theory for the bending of thick sandwich plates.

The analytical solution based on the elasticity theory is important because it can be used as a benchmark to assess other simplified solutions [19]. Based on the elasticity theory, Pagano [20] proposed an analytical solution for the bending of simply supported composite and sandwich plates. The elasticity solution of the cylindrical sandwich shell subjected to radial and axial loads was presented by Kardomateas [21]. Zenkour [22] studied the bending problem of cross-ply laminated plates by the use of the elasticity equations linked with the state space method. Kashtalyan and Menshykova [23] developed elasticity solution for sandwich plates with a FG core.

However, in the above investigations, the permeation effect of the adhesive between layers is neglected. In this paper, 2D elasticity solution is presented for a refined sandwich arch model considering the permeation effect of the adhesive. In the present work, the permeation layer is described as FGM, exponentially graded along the radial direction. The general solution of stresses and deformations with unknown coefficients is obtained by means of the layer-wise method. The coefficients are determined via the recursive matrix method, which is efficient for arches with many layers. The solution obtained is compared with the finite element solution and 1D solution based on the Euler–Bernoulli theory, respectively.

2. Analytical Model

2.1. Basic Equations

Without loss of generality, we consider a sandwich arch with internal radius R₁, external radius R₂, thickness H, and angle , composed of p FG layers with each thickness h_i, as shown in Figure 2. The adjacent layers are perfectly bonded, i.e., the deformations are continuous along the radial direction. A polar coordinate system r − θ is established. The arch is simply supported at two edges and is subjected to radial load acting on the external surface. The permeation layer is the material consisting in the mixture of face (or core) layer material and the adhesive material. Here, we assume its elastic modulus to be exponentially graded along the radial direction, i.e., , in which represents the elastic modulus at the internal surface of the i-th layer and denotes the gradient index. Poisson’s ratio is constant in each layer. The homogeneous layers, i.e., face, core, and adhesive layers, can be included in the model by letting k_i = 0 only.

According to the 2D elasticity theory in the polar coordinate [24], the constitutive relations are given byin which and are the circumferential and radial normal stresses, respectively; is the shear stress; and are the circumferential and radial normal strains, respectively; and is the shear strain. The geometrical equations arein which and are the deformation components in the circumferential and radial directions, respectively. The equilibrium equations are given by

The simply supported boundary conditions of the arch are

The load condition on the external and internal surfaces of the arch can be expressed by

2.2. General Solution for a Sublayer

By combining equations (1)–(3), two partial differential equations with respect to and are derived:

Here, the above equations are solved by the layer-wise method [25]. According to this method, each arch layer is divided into sublayers each with very small thickness so that the variable r in equation (6) can be taken as a constant. The radius, the middle plane of the j-th (j = 1,2, …, pλ) sublayer, is denoted by . The distances from the arch center O to the internal and external surfaces of the j-th sublayer are, respectively, represented by and . Then, we expand the deformation components of the j-th sublayer into Fourier series:where . Substituting equation (7) into equation (6), we obtain two ordinary differential equations:where [j/λ] means the integer part of j/λ. The solution of equation (8) isin which are unknown coefficients, and are the four roots of the following quartic equation:and can be expressed by as follows:

Substituting equation (9) into equation (7), the general solutions of the deformation components with unknown coefficients are obtained:

Substitution of equation (12) into equations (1) and (2) gives the general solutions of the stress components

2.3. Recursive Matrix Method

The out-of-plane variables, i.e., , , , and , of each sublayer can be rearranged into the vector form

By substituting equations (12) and (13) into equation (14), , , , and can be further rearranged into the matrix form as follows:in which

By replacing r in the above equation with and , respectively, one has

By eliminating in the two equations of equation (17), we obtain

The out-of-plane variables are continuous between the adjacent sublayers, i.e.,

By combining equation (18) with equation (19), the relation between the internal and external surfaces of the arch is obtained:

We definein which , , , and are the 2 × 2 submatrixes. Equation (20) can be written as

Then, equation (22) is further decomposed into two matrix equations as follows:

The load function in equation (5) can be expanded into Fourier series

Substituting the load condition into equation (23), one has

By reusing equations (18) and (19), is obtained:

Finally, the coefficients are determined:

Substituting the coefficients back into equations (12) and (13), the solution of stress and deformation components is obtained.

It can be observed that, as the layer number increases, only the computation effort in equation (23) slightly increases. Therefore, the present method is highly efficient for arches with many layers. Moreover, other boundary conditions can also be applied in the present model. For example, the clamped condition can be equivalent to the simply supported one acted by the unknown horizontal reaction, which can be further determined by the zero displacement condition at the clamped end [26].

3. Example and Discussion

In the following examples, the infinite series in the solution are truncated into the finite term, i.e., m = 1,2, …, N.

3.1. Convergence Study

Consider a simply supported sandwich arch with R₁ = 1000 mm, R₂ = 1240 mm, and β = π/4. The face, core, and adhesive layers are made of carbon fiber, polyurethane foam, and epoxy, respectively. The thickness, elastic modulus, and Poisson’s ratio of each arch layer are listed in Table 1.

Since the governing equations of equation (6) are solved by the layer-wise method, the effect of the divided number on the solution accuracy should be analyzed first. We assume the arch is subjected to radial load q(θ) = sin(π/β) N/mm. For sinusoidal load, one series term is sufficient, i.e., N = 1. Table 2 gives the present solution of stresses and deformations in the face layer and permeation layer with different divided number λ, respectively. It can be found that the solution tends to be convergent with the increase of λ. Moreover, the solution in the face layer converges faster than that in the permeation layer. The solution is accurate with at least three significant figures when λ = 10. Therefore, the divided number is taken as λ = 10 for the following calculations.

Then, the convergence property associated with the series terms for the solution is studied. The arch is acted by uniform load q(θ) = 1 N/mm. The present solution of stresses and deformations with different series number N is given in Table 3. It can be found that the present solution is rapidly convergent and has high convergence precision with four significant digits when N = 9. Thus, the number of series terms is fixed at N = 9 unless stated.

3.2. Comparison Study

We still consider the simply supported sandwich arch in Section 3.1; however, the acting load is q(θ) = 1 N/mm and β is the variable. The present solution is compared with the finite element (FE) solution obtained by ANSYS. In the FE modeling, the PLANE-182 element is used to simulate all the layers. Due to the symmetry of the arch, only half part of the arch is modeled, as shown in Figure 3. The FE mesh is created by dividing the arch length into 50 elements, while the thicknesses of face, adhesive, core, and permeation layers is divided into 4, 1, 20, and 10 elements, respectively. Figure 4 shows the comparison of stresses and deformations between the present solution and FE solution with different arch angle β. It can be found that the FE solution is in good agreement with the present one; however, the FE solution with the fine mesh is time consuming, especially in mesh division and calculation. Besides, the stresses and deformations, in absolute values, increase with the increase of β.

(a)

(b)

(c)

(d)

By letting = 0 in the basic equations, the present solution can be used to predict stresses and deformations for the fully homogeneous layered arch. Consider a sandwich arch composed two face layers sandwiching a relatively soft thin interlayer, which was studied by Galuppi and Royer-Carfagni [7] based on the 1D Euler–Bernoulli theory. The geometric and material parameters are fixed at q(θ) = 0.5 N/mm², E₁ = E₃ = 70000 MPa, E₂ = 2.6 MPa, = = = 0.3, R₁ = 1000 mm, β = π/4, h₁ = h₃, and h₂ = 2 mm. We define three variables: = at θ = 0.5β, r = ; = at θ = 0, r = + h₁; and = at θ = 0.5β, r = . The comparison of , , and between the present, FE, and 1D solutions with different layer thicknesses is shown in Table 4. It can be found that the FE solution agrees well with the present solution for all cases. The 1D solution is close to the present one for the thin arch, while the error of 1D solution increases when the arch becomes thick. The errors of and in 1D solution reach 10.1% and 18.5%, respectively.

3.3. Parametric Study

Consider a simply supported sandwich arch with five layers with R₁ = 2000 mm, R₂ = 2100 mm, and β = π/8 acted by distributed load q(θ) = sin(π/β) N/mm. Here, only the permeation effect between the adhesive and core layers is considered. The thickness and the material property of each layer are given in Table 5. The distributions of , , and along the radial direction with different permeation thicknesses h₃ = 0, 10, 15, and 20 mm are shown in Figure 5, in which h₃ = 0 means the special case of the sandwich arch without the permeation layer. It can be found that (i) the absolute value of at the internal and external surfaces in the face layers decreases with the increase of h₃; (ii) the absolute peak value of in the face layers decreases with the increase of h₃, while that in the core layer increases with the increase of h₃; (iii) the distributions of and tend to be smooth with the increase of h₃; and (iv) the distribution of is less affected by h₃.

(a)

(b)

(c)

(d)

(e)

4. Conclusions

Based on the elasticity theory in the polar coordinate, a refined model of the simply supported sandwich arch considering the permeation effect of the adhesive is presented. By means of the recursive matrix method, the solution can be efficiently obtained for the arch with many layers. Conclusions are summarized as follows:(1)The present solution is rapidly convergent with the increase of series term and sublayer number.(2)The finite element solution with the fine mesh is in agreement with the present solution, while the finite element method is time consuming in mesh division and calculation.(3)The 1D solution, based on the Euler–Bernoulli theory, is close to the present solution for thin arches; however, the error increases with the increase of arch thickness.(4)The permeation thickness has considerable effect on the stress and deformation distributions. The distributions of circumferential normal stress and the shear stress tend to be smooth in the radial direction as the permeation thickness increases.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (Grant nos. 51608264, 51978336, and 51778285) and the Natural Science Foundation for distinguished Young Scholars of Jiangsu Province (Grant no. BK20190034).

References

J. B. Yan, J. R. Liew, M. H. Zhang, and K. M. A. Sohel, “Experimental and analytical study on ultimate strength behavior of steel-concrete-steel sandwich composite beam structures,” Materials and Structures, vol. 48, no. 5, pp. 1523–1544, 2015.
View at: Publisher Site | Google Scholar
S. Pappadà, R. Rametta, A. Passaro, L. Lanzilotto, and A. Maffezzoli, “Processing, mechanical properties, and interfacial bonding of a thermoplastic core-foam/composite-skin sandwich panel,” Advances in Polymer Technology, vol. 29, no. 3, pp. 137–145, 2010.
View at: Publisher Site | Google Scholar
C.-S. Chen, F.-H. Liu, and W.-R. Chen, “Dynamic characteristics of functionally graded material sandwich plates in thermal environments,” Mechanics of Advanced Materials and Structures, vol. 24, no. 2, pp. 157–167, 2016.
View at: Publisher Site | Google Scholar
S. Cheng, P. Qiao, F. Chen, W. Fan, and Z. Zhu, “Free vibration analysis of fiber-reinforced polymer honeycomb sandwich beams with a refined sandwich beam theory,” Journal of Sandwich Structures & Materials, vol. 18, no. 2, pp. 242–260, 2016.
View at: Publisher Site | Google Scholar
V. Tahouneh, “Vibration and mode shape analysis of sandwich panel with MWCNTs FG-reinforcement core,” Steel and Composite Structures, vol. 25, no. 3, pp. 347–360, 2017.
View at: Google Scholar
R. Huo, W. Liu, P. Wu, and D. Zhou, “Analytical solutions for sandwich plates considering permeation effect by 3-D elasticity theory, steel,” Steel and Composite Structures, vol. 25, no. 2, pp. 127–139, 2017.
View at: Google Scholar
L. Galuppi and G. Royer-Carfagni, “Analytical approach à la Newmark for curved laminated glass,” Composites Part B: Engineering, vol. 76, pp. 65–78, 2015.
View at: Publisher Site | Google Scholar
H. Raissi, M. Shishehsaz, and S. Moradi, “Stress distribution in a five-layer sandwich plate with FG face sheets using layerwise method,” Mechanics of Advanced Materials and Structures, vol. 26, no. 14, pp. 1234–1244, 2019.
View at: Publisher Site | Google Scholar
M. Mohammadimehr, H. B. Zarei, A. Parakandeh, and A. G. Arani, “Vibration analysis of double-bonded sandwich microplates with nanocomposite facesheets reinforced by symmetric and un-symmetric distributions of nanotubes under multi physical fields,” Structural Engineering and Mechanics, vol. 64, no. 3, pp. 361–379, 2017.
View at: Google Scholar
M. Livani, K. MalekzadehFard, and S. Shokrollahi, “Higher order flutter analysis of doubly curved sandwich panels with variable thickness under aerothermoelastic loading,” Structural Engineering and Mechanics, vol. 60, no. 1, pp. 1–19, 2016.
View at: Publisher Site | Google Scholar
A. Bouchafa, M. B. Bouiadjra, M. S. A. Houari, and A. Tounsi, “Thermal stresses and deflections of functionally graded sandwich plates using a new refined hyperbolic shear deformation theory,” Steel and Composite Structures, vol. 18, no. 6, pp. 1493–1515, 2015.
View at: Publisher Site | Google Scholar
K. T. Nguyen, T. H. Thai, and T. P. Vo, “A refined higher-order shear deformation theory for bending, vibration and buckling analysis of functionally graded sandwich plates,” Steel and Composite Structures, vol. 18, no. 1, pp. 91–120, 2015.
View at: Publisher Site | Google Scholar
O. Rahmani, S. M. R. Khalili, and K. Malekzadeh, “Free vibration response of composite sandwich cylindrical shell with flexible core,” Composite Structures, vol. 92, no. 5, pp. 1269–1281, 2010.
View at: Publisher Site | Google Scholar
B. Liu, A. J. M. Ferreira, Y. F. Xing, and A. M. A. Neves, “Analysis of functionally graded sandwich and laminated shells using a layerwise theory and a differential quadrature finite element method,” Composite Structures, vol. 136, pp. 546–553, 2016.
View at: Publisher Site | Google Scholar
S. Natarajan and G. Manickam, “Bending and vibration of functionally graded material sandwich plates using an accurate theory,” Finite Elements in Analysis and Design, vol. 57, pp. 32–42, 2012.
View at: Publisher Site | Google Scholar
L.-X. Peng, S.-T. Yan, G.-K. Mo, and X. Zhang, “Free vibration analysis of corrugated-core sandwich plates using a meshfree Galerkin method based on the first-order shear deformation theory,” International Journal of Mechanical Sciences, vol. 78, pp. 8–18, 2014.
View at: Publisher Site | Google Scholar
C. H. Thai, A. J. M. Ferreira, S. P. A. Bordas, T. Rabczuk, and H. Nguyen-Xuan, “Isogeometric analysis of laminated composite and sandwich plates using a new inverse trigonometric shear deformation theory,” European Journal of Mechanics—A/Solids, vol. 43, pp. 89–108, 2014.
View at: Publisher Site | Google Scholar
J. L. Mantari, A. S. Oktem, and C. Guedes Soares, “A new trigonometric layerwise shear deformation theory for the finite element analysis of laminated composite and sandwich plates,” Computers & Structures, vol. 94-95, pp. 45–53, 2012.
View at: Publisher Site | Google Scholar
W. Q. Chen and K. Y. Lee, “Three-dimensional exact analysis of angle-ply laminates in cylindrical bending with interfacial damage via state-space method,” Composite Structures, vol. 64, no. 3-4, pp. 275–283, 2004.
View at: Publisher Site | Google Scholar
N. J. Pagano, “Exact solutions for composite laminates in cylindrical bending,” Journal of Composite Materials, vol. 3, no. 3, pp. 398–411, 1969.
View at: Publisher Site | Google Scholar
G. A. Kardomateas, “Elasticity solutions for sandwich orthotropic cylindrical shells under external/internal pressure or axial force,” AIAA Journal, vol. 39, no. 4, pp. 713–719, 2001.
View at: Publisher Site | Google Scholar
A. M. Zenkour, “Three-dimensional elasticity solution for uniformly loaded cross-ply laminates and sandwich plates,” Journal of Sandwich Structures & Materials, vol. 9, no. 3, pp. 213–238, 2007.
View at: Publisher Site | Google Scholar
M. Kashtalyan and M. Menshykova, “Three-dimensional elasticity solution for sandwich panels with a functionally graded core,” Composite Structures, vol. 87, no. 1, pp. 36–43, 2009.
View at: Publisher Site | Google Scholar
S. P. Timoshenko and J. C. Goodier, Theory of Elasticity, McGraw-Hill, New York, NY, USA, 1970.
W. Q. Chen, Y. F. Wang, J. B. Cai, and G. R. Ye, “Three-dimensional analysis of cross-ply laminated cylindrical panels with weak interfaces,” International Journal of Solids and Structures, vol. 41, no. 9-10, pp. 2429–2446, 2004.
View at: Publisher Site | Google Scholar
Y.-p. Xu, D. Zhou, and Y. K. Cheung, “Elasticity solution of clamped-simply supported beams with variable thickness,” Applied Mathematics and Mechanics, vol. 29, no. 3, pp. 279–290, 2008.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Ruili Huo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

308

Downloads

658

Citations