Abstract

An efficient analytical method is proposed so as to study the long-term performance of laminated tubes with neighboring layers glued by viscoelastic interlayers bearing inhomogeneous long-duration loads. In the analysis, the formulations for describing the mechanical performance of each tube layer are based on the elasticity theory, while the viscoelastic behavior of the glue interlayer is modeled by the generalized Maxwell model. Making use of the quasielastic approximation and the recursive matrix method, the analytical solution for the long-term performance of the laminated tube is solved efficiently. The present analytical solution is verified by comparing it with other solutions. The influences of the interlayer thickness, the material property, as well as the load distribution on the long-term performance of the laminated tube are discussed.

1. Introduction

In the past decades, with the progress of structure research, the performance of layered structure has been greatly improved. The high bearing capacity, lightweight, and corrosion resistance characteristics make the layered structures have more possibilities in the structural engineering, for example, structural insulating panels, glued layered timbers, laminated tubes, and many more [13]. Besides the layered beams and plates, the applications about laminated tubes are increasingly applied in infrastructures [47]. Generally, a laminated tube is composed of two or more subcomponents of different materials, and the adjacent members are usually connected discontinuously by studs, nails, and dowels or continuously glued by adhesives [8]. Due to the characteristics of the constituent materials, practical problems exist for the laminated structure under different external conditions [912]. One of them is that the interfacial slip is inevitable to occur because the kinds of connections can hardly be rigid. In addition, the adhesives to bond layers are usually viscoelastic materials, which lead to long-term performance for layered systems. In practical engineering, concrete-filled fiber-reinforced polymer tubes are employed in structural applications including marine constructions and bridges. However, weak adhesion will reduce the rigidity of the pipes, resulting in the interfacial slip [13]. Additionally, the fluid-conveying pipes, which are made up of two reinforced composite layers with a viscoelastic interlayer, are widely applied in the aircraft and mechanical systems. When the pipes are subjected to sustained pressure, the viscoelasticity of the interlayer can affect the long-term performance of the pipe [14]. These phenomena deserve further study and discussion.

A lot of works about the mechanical performance of laminated tubes have been reported in the literature. Xia et al. [15] proposed three-dimensional elasticity solutions to investigate the stress and deformation fields in laminated tubes under internal pressures. Zheng et al. [16] proposed a theoretical model to predict the short-term burst pressure of polyethylene tubes with winding steel wires under different temperatures. The 3-D vibration performance of a periodic composite tube was analyzed by Shen et al. [17] by the use of the transfer matrix method. Based on the three-dimensional anisotropic elasticity, Bakaiyan et al. [18] gave the elasticity solutions for thermal stress and displacement fields of the laminated tubes under temperature environment, and they found that the shear stress and strain in the tube strongly depend on the layer stacking sequence. Hashemian and Mohareb [19] developed a finite difference model to study the buckling of sandwich tubes under external pressure. By using the finite element (FE) modeling, Das and Baishya [20] analyzed the failure modes of layer tubes with socket joints. Exact solutions were proposed by Xia et al. [21] for filament wound fiber reinforced laminated tubes under pure bending load. By employing the FE method, Xu et al. [22] studied the collapse performance of sandwich tubes under external pressures. Zhu et al. [23] analyzed a new GFRP laminated tube with a stiffened core under burial condition in a theoretical way, and they also found a way to improve the strength-to-weight ratio of the tubes by optimizing the geometric parameters. Fernandez-Valdes et al. [24] presented the analysis of friction and initial imperfection of the sandwich tubes by virtue of the FE method. Closed-form analytical solutions of sandwich tubes containing functionally graded interlayers were presented by Sburlati and Kashtalyan [25] based on the elasticity theory, and the results indicated the graded interlayers can remarkably reduce the hoop stress in the tube.

In most of the existing literature for laminated tubes, the neighboring layers are considered to be perfectly connected or have static stiffness, while the long-term performance of viscoelastic adhesives to bond adhesive layers is neglected. Moreover, in most cases, only the homogeneous loads are considered for the laminated tube. In reality, laminated tubes commonly bear inhomogeneous radial loads; for instance, the underground tubes are subjected to inhomogeneous soil pressures along the circumference of the tube [2628].

In the proposed work, an analytical model is presented for studying laminated tubes with infinite length, with neighboring layers glued by viscoelastic interlayers and subjected to inhomogeneous radial loads along the circumferential direction. The formulations for describing the mechanical performance of each tube layer are based on the elasticity theory, while the viscoelastic behavior of the glue interlayer is modeled by the generalized Maxwell model. By means of the quasielastic approximation as well as the recursive matrix method, the solutions for the long-term behavior of the laminated tube are achieved analytically. The present solution based on the elasticity theory has advantages over other solutions. Compared with the Ritz solution, the present solution has higher accuracy because the Ritz method has inherent approximation, and the accuracy of its results depends greatly on the trial functions assumed [29]. The Navier method is often used to solve the mechanical problems of plates and shells, but it has a restrictive condition for boundary [30]. In addition, compared with the Euler–Bernoulli theory [31], the present solution is more accurate since the elasticity theory introduces no assumption of transverse shear deformation.

2. Analytical Model for the Laminated Tube

The reference structure presented in Figure 1 is a laminated tube with inside radius Rin, outside radius Rout, thickness H and infinite length, composed of elastic tube layers glued by viscoelastic interlayers. The thickness of the i-th layer and the interlayer are, respectively, and , where the subscript i means the layer index, and is thin enough in comparison with . The tube bears inhomogeneous radial loads and , respectively, acting on the inner and outer surfaces.

2.1. Equations for an Elastic Tube Layer

The referential tube can be seen as the two-dimensional plan strain problem, because the material properties, geometric shape, and external loads remain unchanged in the longitudinal direction. Based on the elasticity theory [32], the constitutive equations for the i-th (i = 1, 2, …, ) tube layer can be given by the following equation:in which, , , and are the radial, circumferential, and shear stresses, respectively, while , , and are the corresponding strains. The strain-displacement relationships for the i-th tube layer arewhere and are the radial and circumferential displacements, respectively. The stress components keep the equilibrium equations as follows:

Substituting equation (2) into equation (1), one has

Then, by combining equations (3) and (4) and eliminating the stress components, we have

In order to solve the partial differential equations, the displacement components are beforehand expanded into Fourier series, i.e.,

Substitution of equation (6) into equation (5) yields

According to the orthogonality of the Fourier series, a set of ordinary differential equations of , , and is obtained, as follows:

By solving the group of differential equations and substituting the results into equation (4), one obtains the general solutions of displacements and stresses.whereand , , , , , and are unknown coefficients.

2.2. Equations for a Viscoelastic Interlayer

It is verified that the generalized Maxwell model (Figure 2) is appropriate for the viscoelastic bonding interlayer [31, 33], with the modulus given by the following equation:where the mark represents the variable belonging to the interlayer; , , , and denote Poisson’s ratio, the long-term moduli, the relaxation moduli, and the relaxation time. By the use of the quasielastic approximation, the constitutive equations of the i-th (i = 1, 2, …, ) interlayer are given by [34] the following equation:

The strain-displacement relationships for the thin interlayer arewhere and . In addition, the equilibrium relations of the interlayer and the neighboring tube layers are

By combining equations (11)–(13) and eliminating the stresses and strains of the interlayer, one obtains

2.3. Coefficient Determination by Recursive Matrix Method

In this section, the coefficients in the general solutions of equation (9) are determined by the recursive matrix method. The loads acting on the tube can be expressed by the following equation:in which loads functions are expanded aswhere

The out-of-plane displacements and stresses in equation (12) are rearranged by matrix form.in which the functions , , , , , and (n = 1,2) are also written into matrix form aswhere

By assigning r in the above equation as and , respectively, we have

By eliminating and in equation (22), one obtains

Equation (15) can also be rearranged by matrix form:where

Combining equation (23) and equation (24), one obtains

The products of the matrixes in equation (26) are then defined by 2 × 2 partitioned matrixes:

By inserting equation (27) into equation (26), the equations become

By decomposing the above matrix equations, one has

Substitution of equation (19) into equations (16) and (17) yields

The solutions of the second and fourth equations of equation (29) are

Therefore, and are determined. By reusing equations (23) and (24), and are obtained:

At last, the coefficients are determined as

What needs to be pointed out is that by the use of the recursive matrix method, the coefficients are determined by calculating a series of small matrixes instead of solving a large matrix problem; therefore, the present analytical model is very efficient.

3. Examples and Discussion

In order to investigate the effect of inhomogeneous degree of the load distribution, an elliptical load is employed and applied on the external surface of the tube, which is expressed bywhere k = b/a; a and b are half lengths of semiminor axis and major axis in the ellipse, respectively, as shown in Figure 3. The internal load is fixed at  = 0.2 N/mm2. Some representative stresses and displacements are defined in advance:

Besides, the symbol denotes the absolute value, e.g., , and the subscript max means the maximum value, e.g., .

3.1. Verification of the Present Solution

For the solution of stresses and displacements expressed in series form, the convergence characteristics should be checked firstly. For this purpose, a three-layer tube glued by polyvinyl butyral (PVB) is considered as an example, with parameters fixed at k = 2,  =   = 100 GPa,  = 50 GPa,  =   =   =   = 0.3, Rin = 900 mm, h1 = h3 = 50 mm, h2 = 100 mm, and  = 0.3 mm, and the viscoelastic constants of PVB are listed in Table 1.

Table 2 gives the values of the stresses and displacements for different numbers of terms N = 5, 10, 15, 20, respectively. It can be observed that the stress and displacement components are convergent rapidly with the precision of three significant figures when N = 10.

Then, to validate the present analytical solution, the FE solution from the software ANSYS is used to make a comparison. Due to the symmetry of the structure, only 1/4 part of the tube is modeled, and the FE mesh size is fine enough, as illustrated in Figure 4. Here, both tube layers and interlayers in the tube are simulated by the PLANE-183 element. Figure 5 displays the comparison of , , and along the radial direction as well as along the circumferential direction. It can be obtained that the FE solution matches the present one well (see Figures 5(a)5(d). However, it should be mentioned that the calculation of the FE solution is time-consuming due to fine mesh division.

3.2. Parametric Analysis

In this section, a sandwich tube with PVB interlayer under radial loads is considered with Rin = 900 mm, h1 = h3 = 50 mm, h2 = 100 mm,  =   =   =   = 0.3, while the other parameters are variable.

Figure 6 presents the variations of , , and with t; meanwhile,  =   = 100 GPa,  = 50 GPa, k = 2. It can be found from Figure 6 that and increase and remain unchanged eventually as t increases (see Figures 6(a) and 6(c)), while decreases gradually with t (see Figure 6(b)). This is mainly resulted from that is degenerated with t and approaches the fixed value as t. Moreover, the thicker interlayer leads to the increase of and (see Figures 6(a) and 6(c)), and by contrast, decreases with the interlayer thickness (see Figure 6(b)). The main reason for this phenomenon is that the shear stiffness between layers, i.e., , decreases as increases.

Figure 7 illustrates the effects of k and t on , , and , with  = 0.3 mm,  =   = 100 GPa, and  = 50 GPa. In the mass, , , and increase with t and remain unchanged eventually, while for k = 1, , and only have small and positive values at any time (see Figures 7(a)7(c)). For any of t, , , and increase with k significantly (see Figures 7(a)7(c)). In other words, the inhomogeneous degree of the load distribution has a large effect on the stresses and displacements in the tube.

4. Conclusions

According to the present work, some conclusions can be summarized:(1)The present solution has high accuracy and matches well with the FE solution. However, the FE method is time-consuming due to fine mesh division.(2)For the long-term behavior of the laminated tube, the circumferential stresses as well as the radial deformation increase with the time, while the shear stress in the interlayer decreases with the time. All the stresses and deformations tend to remain unchanged eventually because the material property in the interlayer degenerates into a constant value in the long term.(3)The long-term behavior of the laminated tube is largely affected by the inhomogeneous degree of the load distribution. When the load acting on the tube surface is homogeneous, the values of stresses and displacements are very small, while they significantly increase with the inhomogeneous degree of load distribution.

Data Availability

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 that they have no conflicts of interest with respect to the research, authorship, and/or publication of this article.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 52108220).