Next Article in Journal
Bending Force of Hot Rolled Strip Based on Improved Whale Optimization Algorithm and Twinning Support Vector Machine
Next Article in Special Issue
Investigation on the Electrochemical Corrosion Behavior of TP2 Copper and Influence of BTA in Organic Acid Environment
Previous Article in Journal / Special Issue
Microstructure and Properties of Laser Surface Remelting AlCoCrFeNi2.1 High-Entropy Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Collapse of Externally Pressurized Steel–Composite Hybrid Cylinders: Analytical Solution and Experimental Verification

1
School of Mechanical Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China
2
School of Naval Architecture & Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China
3
China Ship Scientific Research Center, Wuxi 214082, China
*
Authors to whom correspondence should be addressed.
Metals 2022, 12(10), 1591; https://doi.org/10.3390/met12101591
Submission received: 20 July 2022 / Revised: 17 September 2022 / Accepted: 20 September 2022 / Published: 24 September 2022

Abstract

:
To evaluate the collapse pressure of the steel–composite hybrid cylinders under external pressure without excessive computational cost, an analytical formula was derived in this study. The rationality of the derived formula was verified by the comparison with experimental and numerical results. The experimental results indicate that samples are manufactured and tested with good quality. The derived formula considered material failure and could reasonably predict the collapse pressure of the steel–composite hybrid cylinders with a maximum difference of 3.1%. Moreover, the effects of the wrap angle, thickness, and length on the collapse pressure of the hybrid cylinders were theoretically analyzed. The loading capacity of the hybrid cylinders was maximized under a wrap angle of ±55° for the composite layer. These findings are mainly because the hoop stress is twice the value of axial stress for a cylinder under uniform pressure.

1. Introduction

Steel cylindrical shells have a robust design, can be produced using mature techniques, and are widely used in ships and marine structures, such as pressure vessels and pipelines used in autonomous underwater vehicles, human-occupied vehicles, and subsea engineering applications [1,2,3,4,5]. However, steel cylinders are prone to instability or collapse under external pressure [6]. The instability of steel cylinders has been investigated analytically and experimentally [3,7,8,9,10]. Moreover, steel cylinders are sensitive to imperfections, which can considerably reduce their collapse pressure. Cylinders with imperfections thus have poor structural efficiency.
Therefore, steel–composite hybrid cylinders with improved structural efficiency have been proposed. These hybrid cylinders have attracted considerable research attention because they possess the advantages of steel and composite cylinders [11,12,13,14,15]. The inclusion of composite layers can considerably improve the structural efficiency of a steel cylinder. These layers have high specific strength and stiffness [16,17], flexible manufacturing requirements, excellent energy absorption ability [18], and high corrosion resistance [19,20]. Many studies have numerically and experimentally investigated hybrid cylinders subjected to axial and internal pressure [21,22]. Teng and Hu [21] experimentally demonstrated that external wrapping with composite layers enhanced the ductility of steel cylinders under axial compression. Subsequently, Vakili [23] used experimental and numerical approaches to demonstrate that local wrapping with composite layers substantially increases the loading capacity of an elephant-foot cylinder under axial and internal pressure. However, few studies have experimentally evaluated the collapse properties of steel–composite hybrid cylinders under external pressure.
The collapse pressure of steel and composite cylinders is typically evaluated using an analytical approach. Numerous studies have developed analytical approaches for quickly predicting the buckling or collapse pressure of steel-only cylinders. Some formulas, such as the Venstel–Krauthaer [24], Ross [25], CCS2013 (submersible specifications of the China Classification Society) [26], and ABS2012 (submersible specifications of the American Bureau of Shipping) [27] formulas, have been proposed for evaluating the loading capacity of cylinders under uniform external pressure.
For composite-only cylinders, several analytical formulas have been developed for the rapid prediction of the buckling pressure [28,29,30]. For example, Rosenow [31] developed an analytical solution based on classical lamination theory (CLT) to predict the buckling pressure of composite cylinders and reported that the filament wound pipe should be wound at 54.75° for optimal biaxial pressure loading. Imran and Shi [32,33] used linear eigenvalues to optimize the design of composite cylinders subjected to external uniform pressure. They later [34] developed an analytical model based on first-order shear deformation theory for evaluating the buckling of composite cylinders. Given multiple boundary conditions, Lopatin [35,36] predicted the critical buckling pressure of composite cylinders with various closed ends by using an analytical method. Bai [37,38,39] proposed an analytical buckling model for reinforced thermoplastic pipes that can describe cases with various load types, such as bending, tension, and external pressure. Moreover, Messager [40] proposed a linear Sanders-type buckling model of cross-ply cylinders in which winding-induced geometrical imperfections were considered.
However, no study has derived an analytical formula for predicting the collapse pressure of steel–composite hybrid cylinders under external pressure. Thus, a reasonable formula without excessive computations should be developed for predicting the collapse pressure of such cylinders.
In this study, an analytical approach and experimental approach were adopted for evaluating the collapse pressure of steel–composite hybrid cylinders under external pressure. The remainder of this paper is organized as follows: Section 2 details an analytical formula designed to evaluate the collapse pressure of steel–composite hybrid cylinders under external pressure. Section 3 describes a comparison of analytical, experimental, and numerical results for steel-only and steel–composite hybrid cylinders. This section also provides a theoretical explanation of the effects of the wrap angle, thickness, and length on the collapse pressure of such cylinders. Finally, Section 4 presents the conclusions of this study and provides suggestions for future research. The results indicate that the proposed formula can reasonably predict the collapse pressure of steel–composite hybrid cylinders. The maximum difference between the predictions and the experimental results was 3.1%. Moreover, the loading capacity of the hybrid cylinders was maximized under a wrap angle of ±55° for the composite layer.

2. Materials and Methods

An analytical formula was designed to evaluate the collapse pressure of steel–composite hybrid cylinders under external pressure. For comparison, these pressures were also calculated using the NASA SP-8700 [28] and ASME Code 2007 (RD-1172) [29] formulas on the basis of a modified stiffness matrix.

2.1. Problem Definition

A steel cylinder strengthened with composite is investigated in this study. The steel–composite hybrid cylinder comprises an outer composite layer and an inner steel layer (Figure 1). The grade of the inner steel layer is O6Cr19Ni10 according to GB/T 3280-2015 [41]. This grade is equivalent to 304 according to ASTM A959-16 [42]. The carbon-fiber-reinforced polymer (CFRP) prepreg (CFS24T300) was adopted in the outer composite layer.
The two ends of the hybrid cylinder are enclosed by heavy bungs, which impose rigid boundary conditions on the cylinder’s ends [43,44]. The outer surfaces of the hybrid cylinder and the bungs are subjected to a uniform external hydrostatic pressure P0. This loading causes the radial and axial compression of the cylinder (Figure 2).
The steel–composite hybrid cylinder has a length of L, a total thickness of h, and an inner radius of R. Moreover, the wrap angle of its composite layer is θ (Figure 3). The thicknesses of the steel and composite layers are ts and tc, respectively. The outer radius of the cylinder is Rs. The thickness of a single composite layer is tl. The steel layer can be transformed into n equivalent isotropic single layers with a thickness of tl each. Therefore, if the total number of layers is N, the numbers of steel and composite single layers are n and N − n, respectively.

2.2. Mechanical Properties

According to Hooke’s law, the stress and strain in a single layer (Figure 4) have the following relationship:
ε 1 = σ 1 E 1 ,   ε 2 = ν 12 ε 1 = ν 12 σ 1 E 1 ,   γ 12 = τ 12 G 12 ,
ε 2 = σ 2 E 2 ,   ε 1 = ν 21 ε 2 = ν 21 σ 2 E 2 ,   γ 12 = τ 12 G 12 ,
where E1 and E2 represent the Young’s moduli of the fiber and matrix orientations, respectively.
By substituting Equations (1) and (2) into Equation (3), the following matrix is obtained:
[ ε 1 ε 2 γ 12 ] = S [ σ 1 σ 2 τ 12 ] = [ S 11 S 12 0 S 21 S 22 0 0 0 S 66 ] [ σ 1 σ 2 τ 12 ] ,
where S 11 = 1 E 1 , S 22 = 1 E 2 , S 12 = S 21 = ν 21 E 2 = ν 12 E 1 , and S 66 = 1 G 12 . In isotropic materials, S 11 = S 22 = 1 E , S 12 = S 21 = ν E , and S 66 = 2 ( 1 + ν ) E .
[ σ 1 σ 2 τ 12 ] = Q [ ε 1 ε 2 γ 12 ] = [ Q 11 Q 12 0 Q 21 Q 22 0 0 0 Q 66 ] [ ε 1 ε 2 γ 12 ] ,
where Q = S 1 represents the reduced stiffness matrix. On the basis of Equation (4), the stresses in layer k can be expressed as follows (Figure 5):
{ σ } k = [ Q ¯ ] k { ε } k ,   [ σ x σ y τ x y ] = Q ¯ [ ε x ε y γ x y ] = [ Q ¯ 11 Q ¯ 12 Q ¯ 16 Q ¯ 21 Q ¯ 22 Q ¯ 26 Q ¯ 16 Q ¯ 26 Q ¯ 66 ] { [ ε x 0 ε y 0 γ x y 0 ] + t l [ q x q y q x y ] } ,
where tl and q represent the thickness of a single layer and the curvature, respectively.
Q ¯ 11 = cos 4 θ k Q 11 + sin 4 θ k Q 22 + 2 cos 2 θ k sin 2 θ k ( Q 12 + 2 Q 66 ) ; Q ¯ 22 = sin 4 θ k Q 11 + cos 4 θ k Q 22 + 2 cos 2 θ k sin 2 θ k ( Q 12 + 2 Q 66 ) ; Q ¯ 66 = ( Q 11 + Q 22 2 Q 12 ) cos 2 θ k sin 2 θ k + ( cos 2 θ k sin 2 θ k ) 2 Q 66 ; Q ¯ 12 = Q ¯ 21 = ( cos 4 θ k + sin 4 θ k ) Q 12 + ( Q 11 + Q 22 4 Q 66 ) cos 2 θ k sin 2 θ k ; Q ¯ 16 = Q ¯ 61 = [ cos 2 θ k ( Q 11 Q 12 2 Q 66 ) + sin 2 θ k ( Q 12 Q 22 + 2 Q 66 ) ] cos θ k sin θ k ; Q ¯ 26 = Q ¯ 62 = [ cos 2 θ k ( Q 12 Q 22 + 2 Q 66 ) + sin 2 θ k ( Q 11 Q 12 2 Q 66 ) ] cos θ k sin θ k ;
where θk represents the fiber orientation of the composite layer (Figure 5).
For isotropic materials, the following equations are obtained:
[ σ x σ y τ x y ] = [ σ 1 σ 2 τ 12 ] ,
[ σ x σ y τ x y ] = [ Q ¯ 11 Q ¯ 12 Q ¯ 16 Q ¯ 21 Q ¯ 22 Q ¯ 26 Q ¯ 16 Q ¯ 26 Q ¯ 66 ] { [ ε x 0 ε y 0 γ x y 0 ] + t l [ q x q y q x y ] } = [ Q 11 Q 12 0 Q 21 Q 22 0 0 0 Q 66 ] { [ ε x 0 ε y 0 γ x y 0 ] + t l [ q x q y q x y ] } ,
The equilibrium conditions applied to the force system in the composite and steel layers are presented in Figure 5 and Figure 6.
The following equilibrium equations are obtained on the basis of CLT:
{ N x N y N x y } = h / 2 h / 2 { σ x σ y σ x y } d z = k = 1 n z k 1 z n { σ x s σ y s σ z y s } k d z + k = n + 1 N z k 1 z N { σ x c σ y c σ z y c } k d z ,
{ M x M y M x y } = t 2 t / 2 { σ x σ y σ x y } z d z = k = 1 n z k 1 z n { σ x s σ y s σ z y s } k z d z + k = n + 1 N z k 1 z N { σ x c σ y c σ z y c } k z d z ,
where −s and −c represent the steel and composite layers, respectively; n = ts/tl; ts is the thickness of the steel layer (Figure 2); h is the total thickness of the hybrid cylinder; N = n + tc/tl; and tc is the thickness of the composite layer (Figure 2).
By substituting Equations (5) and (7) into Equations (8) and (9), the following equations are obtained:
N x = k = 1 n z k 1 z n ( Q 11 ) k ( ε x 0 + z q x ) + ( Q 12 ) k ( ε y 0 + z q y ) d z + k = n + 1 N z k 1 z N ( Q ¯ 11 ) k ( ε x 0 + z q x ) + ( Q ¯ 12 ) k ( ε y 0 + z q y ) + ( Q ¯ 16 ) ( γ x y 0 + z q x y ) d z ,
N y = k = 1 n z k 1 z n ( Q 21 ) k ( ε x 0 + z q x ) + ( Q 22 ) k ( ε y 0 + z q y ) d z + k = n + 1 N z k 1 z N ( Q ¯ 21 ) k ( ε x 0 + z q x ) + ( Q ¯ 22 ) k ( ε y 0 + z q y ) + ( Q ¯ 26 ) ( γ x y 0 + z q x y ) d z ,
N x y = k = 1 n z k 1 z n ( Q 66 ) k ( γ x y 0 + z q x y ) d z + k = n + 1 N z k 1 z N ( Q ¯ 16 ) k ( ε x 0 + z q x ) + ( Q ¯ 26 ) k ( ε y 0 + z q y ) + ( Q ¯ 66 ) ( γ x y 0 + z q x y ) d z ,
By merging the resultant equation, the following equations are obtained:
N x = ( k = 1 n z k 1 z n ( Q 11 ) k d z + k = n + 1 N z k 1 z N ( Q ¯ 11 ) k d z ) ( ε x 0 ) + ( k = 1 n z k 1 z n ( Q 12 ) k d z + k = n + 1 N z k 1 z N ( Q ¯ 12 ) k d z ) ( ε y 0 ) + ( k = n + 1 N z k 1 z N ( Q ¯ 16 ) k d z ) ( γ x y 0 ) + ( k = 1 n z k 1 z n ( Q 11 ) k z d z + k = n + 1 N z k 1 z N ( Q ¯ 11 ) k z d z ) ( q x ) + ( k = 1 n z k 1 z n ( Q 12 ) k z d z + k = n + 1 N z k 1 z N ( Q ¯ 12 ) k z d z ) ( q y ) + ( k = n + 1 N z k 1 z N ( Q ¯ 16 ) k z d z ) ( q x y ) = A 11 ( ε x 0 ) + A 12 ( ε y 0 ) + A 16 ( γ x y 0 ) + B 11 ( q x ) + B 12 ( q y ) + B 16 ( q x y ) ,
N y = ( k = 1 n z k 1 z n ( Q 21 ) k d z + k = n + 1 N z k 1 z N ( Q ¯ 21 ) k d z ) ( ε x 0 ) + ( k = 1 n z k 1 z n ( Q 22 ) k d z + k = n + 1 N z k 1 z N ( Q ¯ 22 ) k d z ) ( ε y 0 ) + ( k = n + 1 N z k 1 z N ( Q ¯ 26 ) k d z ) ( γ x y 0 ) + ( k = 1 n z k 1 z n ( Q 21 ) k z d z + k = n + 1 N z k 1 z N ( Q ¯ 21 ) k z d z ) ( q x ) + ( k = 1 n z k 1 z n ( Q 22 ) k z d z + k = n + 1 N z k 1 z N ( Q ¯ 22 ) k z d z ) ( q y ) + ( k = n + 1 N z k 1 z N ( Q ¯ 26 ) k z d z ) ( q x y ) = A 21 ( ε x 0 ) + A 22 ( ε y 0 ) + A 26 ( γ x y 0 ) + B 12 ( q x ) + B 22 ( q y ) + B 26 ( q x y ) ,
N x y = k = n + 1 N z k 1 z N ( Q ¯ 16 ) k d z ( ε x 0 ) + k = n + 1 N z k 1 z N ( Q ¯ 26 ) k d z ( ε y 0 ) + ( k = 1 n z k 1 z n ( Q 66 ) k d z + k = n + 1 N z k 1 z N ( Q ¯ 66 ) k d z ) ( γ x y 0 ) + k = n + 1 N z k 1 z N ( Q ¯ 16 ) k z d z ( q x ) + k = n + 1 N z k 1 z N ( Q ¯ 66 ) k z d z ( q y ) + ( k = 1 n z k 1 z n ( Q 66 ) k z d z + k = n + 1 N z k 1 z N ( Q ¯ 66 ) k z d z ) ( q x y ) ,
where Aij and Bij represent the extensional stiffness and bending extensional stiffness of the hybrid structure, respectively. By substituting Equations (5) and (7) into Equations (8) and (9) and merging the resultant equations, the following expression is obtained:
[ M x M y M x y ] = [ B 11 B 12 B 16 D 11 D 12 D 16 B 12 B 22 B 26 D 12 D 22 D 26 B 16 B 26 B 66 D 16 D 26 D 66 ] [ ε x 0 ε y 0 ε x y 0 q x q y q x y ] ,
where Dij is the bending stiffness of the hybrid structure and is calculated using Equations (5), (7)–(9), and (16).

2.3. Analytical Formulas

2.3.1. Designed Formula

For simultaneously considering the elastic buckling mode and material failure mode, an analytical formula was designed for calculating the collapse pressure of steel–composite hybrid cylinders under external pressure (Figure 6). The Merchant–Rankine formula can be used to determine the collapse load of a structure with a rigid support with satisfactory accuracy [45]. The collapse pressure (Pc) formula proposed in this study was derived from the measured collapse pressure of hybrid cylinder models. The Merchant–Rankine equation was marginally modified by introducing a numerical factor b. This equation comprises the elastic buckling pressure (Pe) and material failure pressure (Pf) for steel–composite hybrid cylinders. The modified Merchant–Rankine formula is expressed as follows:
1 P c = 1 P e + b P f ,
where P c , P e , and P f represent the collapse pressure, elastic buckling pressure, and material pressure of the hybrid cylinder, respectively, and b = 0.1.
The steel–composite hybrid cylinder deforms under external pressure, and this deformation progresses until it collapses. During deformation, internal forces act on the steel and composite layers. The steel and composite materials begin to fail under the load at which the stresses match the tensile, compressive, transverse, or shear strength of the material, and the steel–composite hybrid cylinder is assumed to collapse under the failure pressure. The force resultant (Equation (18)) was determined on the basis of the constitutive equation. The stresses under external pressure can be calculated using Equations (5) and (18).
{ ε } = [ A i j ] 1 { N } ,
where Nx = −PR/2, Ny = −PR, and Nxy = 0 are the membrane forces, and P represents the applied pressure.
For each composite single layer, the total stress is separated into the fiber-directional stress ( σ 1 c ), in-plane transversal stress ( σ 2 c ), and in-plane shear stress ( τ c ). For each steel single layer, the total stress is separated into the in-plane stress ( σ 1 s ) and shear stress ( τ s ). The ratios of these stresses to the tensile, compressive, and shear strengths can be calculated to determine the maximum failure index H. Therefore, the material failure pressure ( P f ) of the hybrid cylinder can be obtained by dividing the applied pressure (P) by the maximum failure index (H) as follows:
P f = P H ,   H = M a x { | σ 1 c | X T , | σ 1 c | X C , | σ 2 c | Y T , | σ 2 c | Y C , | τ c | S c , | σ 1 c | T s , | σ 1 c | C s , | σ 1 c | S s } ,
where X T , X C , Y T , Y C , and S c represent the fiber tensile, fiber compressive, matrix tensile, matrix compressive, and shear strengths of the composite layer, respectively. Moreover, T s , C s , and S s represent the tensile, compressive, and shear strengths of the steel layer, respectively.
The proposed analytical model for elastic buckling ( P e ) does not consider transverse shear effects. The presented linear buckling analysis is based on the principles of the Messager [40] model. Equilibrium is achieved between the internal and external forces acting in the cylindrical coordinate system on a shell element. The following three partial differential equations are satisfied:
N x x + N x y y = 0 ,
N x y x + N y y = 0 ,
2 M x x 2 + 2 2 M x y x y + 2 M y y 2 N y R + N x 0 2 ω x 2 + N y 0 2 ω y 2 = 0 ,
where N x 0 and N y 0 represent the membrane forces.
According to thin shell theory, the strain–displacement relationship [40,46] is expressed as follows:
[ ε x 0 ε y 0 γ x y 0 ] = [ x 0 0 0 y 1 R y x 0 ] [ u v w ] ,   [ q x q x q x y ] = [ 0 0 2 x 2 0 R y 2 y 2 0 R x 2 2 x y ] [ u v w ] ,
Moreover, the displacement approximation functions indicating the shape of the instability mode of the cylindrical shell can be expressed as presented in Equation (24). The cylindrical shell subjected to uniform external pressure is simply supported at edges x = 0 and x = L.
{ u = U cos α x cos β y v = V sin α x sin β y ω = W sin α x cos β y ,
where α = mπ/L and β = n/R; m and n are the numbers of longitudinal and circumferential half waves, respectively; R is the inner radius of the cylindrical shell; and U, V, and W are unknown constants that represent the maximum displacements in the x-direction, y-direction, and z-direction, respectively.
By substituting Equations (13)–(16), (23), and (24) into Equations (20)–(22) and integrating the governing equation, an eigenvalue problem with a simple form (Equation (25)) is obtained. The lowest Pe value is the elastic buckling pressure of the hybrid structure.
P e = P e m F C ,   ( [ K ] + P e m [ L ] ) [ U V W ] = [ 0 0 0 ] ,
where K 11 = A 11 α 2 + A 66 β 2 , K 12 = K 21 = ( A 12 + A 66 ) α β , K 13 = K 31 = A 12 R α + B 11 α 3 + ( B 12 + 2 B 66 ) α β 2 , K 22 = A 22 β 2 + A 66 α 2 , K 23 = K 32 = ( B 12 + 2 B 66 ) α 2 β + A 22 R β + B 22 β 3 , K 33 = D 11 α 4 + 2 ( D 12 + 2 D 66 ) α 2 β 2 + D 22 β 4 + A 22 R 2 + 2 B 12 R α 2 + 2 B 22 R β 2 , L 11 = L 12 = L 13 = L 21 = L 22 = L 23 = L 31 = L 32 = 0 , and L 33 = R 2 α 2 R β 2 . The parameters Aij, Bij, and Dij are computed using Equations (13)–(16), and FC is the safety factor.

2.3.2. Modified NASA SP-8700 Formula

The NASA SP-8700 formula [28], which is based on CLT, can be used to determine the buckling pressure of a composite cylinder. By using the modified stiffness matrix (Section 2.2), this formula (Equation (26)) can be used to calculate the collapse pressure of a steel–composite hybrid cylinder.
P N A S A m o = M i n ( P a m n ) ,   P a m n = R F N A S A [ n 2 + 1 2 ( m π R L ) 2 ] det [ C 11 C 12 C 13 C 21 C 22 C 23 C 31 C 32 C 33 ] det [ C 11 C 12 C 21 C 22 ] ,
where C 11 = A 11 ( m π L ) 2 + A 66 ( n R ) 2 , C 22 = A 22 ( n R ) 2 + A 66 ( m π L ) 2 , C 33 = D 11 ( m π L ) 4 + 2 ( D 12 + 2 D 66 ) ( m π L ) 2 ( n R ) 2 + D 22 ( n R ) 4 + A 22 R 2 + 2 B 22 R ( n R ) 2 + 2 B 12 R ( m π L ) 2 , C 13 = C 31 = A 12 R ( m π L ) + B 11 ( m π L ) 3 + ( B 12 + 2 B 66 ) ( m π L ) ( n R ) 2 , C 12 = C 21 = ( A 12 + A 66 ) ( m π L ) ( n R ) , C 23 = C 32 = ( B 12 + 2 B 66 ) ( m π L ) 2 ( n R ) + A 22 R ( n R ) + B 22 ( n R ) 3 .
The parameters R and L represent the inner radius and length of the cylinder, respectively, and m and n are integers necessary to produce the correct minimum pressure. The term FNASA represents the safety factor (1~5). The C matrices are computed using Equations (13)–(16).

2.3.3. Modified ASME Code 2007 (RD-1172) Formula

The ASME Code 2007 formula [29] (Equation (27)) describes the buckling pressure of a composite cylinder and is used to estimate the allowable pressure. By using the modified stiffness matrix of the hybrid cylinder, Equation (27) can be used to calculate the collapse pressure as follows:
P A S M E m o = ( K D ) 0.8531 γ E h f 3 / 4 E a t 1 / 4 t 5 / 2 ( 1 v x v y ) 3 / 4 L ( R ) 3 / 2 F A S M E ,
where v x = ( A B D 1 ) 5 , 4 ( A B D ) 4 , 4 , v y = ( A B D 1 ) 5 , 4 ( A B D 1 ) 5 , 5 , E a f = 12 t 3 ( A B D 1 ) 4 , 4 , E h f = 12 t 3 ( A B D 1 ) 5 , 5 , and E a t = A 11 A 22 A 12 2 A 22 t . The matrices A, B, and D are computed using Equations (13)–(16), KD is the knockdown factor (0.84), and γ is the reduction factor. If Zp ≤ 100, γ = 1 0.001 Z p . If Zp > 100, γ = 0.9. Furthermore, t is the wall thickness, and FASME is the safety factor. The parameter Zp is calculated as follows:
Z p = E h f 3 / 2 E a t 1 / 2 E a f 2 ( 1 v x v y ) 1 / 2 L 2 ( R t ) ,

3. Results and Discussion

This section presents a comparison of the analytical and experimental results obtained for steel-only and steel–composite hybrid cylinders. The elastic buckling pressures of the hybrid cylinders were calculated using analytical formulas and compared with numerical results. Moreover, the effects of the wrap angle, thickness, and length on the collapse pressure of hybrid cylinders were theoretically analyzed.

3.1. Verification of the Analytical Model

3.1.1. Steel-Only Cylinders

Steel–composite hybrid cylinders comprise an outer composite layer and an inner steel layer. The steel layer can be transformed into n equivalent isotropic single layers (Section 2). To verify the rationality of this transformation, the calculated collapse pressure of steel-only cylinders (tc = 0) was compared with the experimental results.
Table 1 presents the collapse pressure of steel-only cylinders obtained using the proposed formula (PC), using the CSS2013 formula (PCCS) [26], using the ABS2012 formula (PABS) [27], and in the experiments (Ptest) of Zhu [47]. The safety factor (FC) was 2.8. The nominal radius and thickness of the steel-only cylinders were 50 and 1 mm, respectively, in the experiments. In Table 1, L/R indicates the ratio of the cylinder length to the cylinder radius.
Table 2 presents the collapse pressure of steel-only cylinders obtained using the derived formula (PC) (Equation (17)), using the Venstel–Krauthaer formula ( P V - K   (Equation (29)) [24], using the Ross formula ( P R o s s ) (Equation (30)) [25], and in the experiments (Ptest) of Zhang [6]. The safety factor was 2.8. The nominal radius, length, and thickness of the steel-only cylinders (SL) were 35, 130, and 1 mm, respectively, in the experiments.
P V - K   = 0.92 E t l ( t r ) 1.5 ,
P R o s s = 2.6 E ( t 2 r ) 2.5 l 2 r 0.45 ( t 2 r ) 0.5 ,
The constitutive material used in the calculations was the same as that adopted in the corresponding experiments [6,47]. Table 1 and Table 2 reveal that the proposed formula determined the collapse pressure of the steel-only cylinders with acceptable accuracy in both cases. The errors were 0.020–0.069 and 0.027–0.031 for verification examples I and II, respectively. The minimum errors were 0.020 and 0.027, which are lower than those of the other analytical formulas.

3.1.2. Steel–Composite Hybrid Cylinders

To verify the rationality of the analytical model for the collapse pressure of a steel–composite hybrid cylinder, the [55/−55]4 and [90/90/0/90/90/0/90/90] orientations were chosen as examples of test validation. Five steel–composite hybrid cylinders comprising outer carbon-fiber-reinforced polymer (CFRP) layers and an inner steel layer were fabricated. The flowchart of sample fabrication is depicted in Figure 7a. The hybrid cylinders were enclosed by O-rings acting as heavy steel bungs. The thickness of the steel bungs was 30 mm, and the nominal thickness (ts) and outer radius (Rs) of the steel layer and cylinders were 1.5 and 79.5 mm, respectively. The Young’s modulus and Poisson’s ratio of the steel were 190 GPa and 0.3, respectively. These heavy steel bungs imposed relatively rigid boundary conditions on the cylinder ends. The hand lay-up method was adopted to fabricate the steel–composite hybrid structure. To achieve excellent bonding, the outer surface of the steel cylinder was abrasively blasted and then cleaned using acetone. The steel cylinder was then wrapped using eight CFRP prepreg (CFS24T300) layers (Jiangsu Boshi Carbon Fiber Technology). The nominal thickness (tl) of the one-ply prepreg was 0.15 mm. The mechanical properties of CFRP were provided by the manufacturer and were as follows: E11 = 115 GPa, E22= 7.72 GPa, G12 = 3.72 GPa, Nu12 = 0.33, XT = 1400.09 MPa, XC = 580.06 MPa, YT = 44.36 MPa, YC = 133.03 MPa, and SC = 45.04 MPa. Three hybrid cylinders exhibited wrap sequences of [55/−55]4, where 90° represents the circumferential direction of the cylinder, and were denoted as CYH1, CYH2, and CYH3. Two hybrid cylinders had wrap sequences of [90/90/0/90/90/0/90/90] and were denoted as CYL1 and CYL2. The lengths of the CYH and CYL cylinders were 320 and 280 mm, respectively. Each specimen was varnished using polyurea coating GDJN001 to prevent water absorption [48]. The polyurea-coated films had a total thickness of 0.25 mm. The polyurea layers were applied at 5 h intervals under a temperature of <30 °C.
The testing chamber was an internal cylindrical pressure vessel (Figure 7c) located at Jiangsu University of Science and Technology, Zhenjiang, China. The inner diameter, total inner height, and maximum pressure of the testing chamber were 200 mm, 1000 mm, and 8 MPa, respectively. Each sample was submerged in water without any constraints. A manual water pump was used to apply pressure. The destruction of each sample was accompanied by a loud noise and a decrease in pressure. After the flanges were removed using a heat gun (Figure 7d), the collapse modes of the five samples were observed (Figure 8). The collapse modes appear as local dents. Such local dents are associated with the initial geometric imperfections and plasticity of the steel layer. These results indicate that the collapse mode of such a hybrid structure is similar to the typical characteristic of the shell of revolution under uniform external pressure. As the composite layer has poor ductility, fiber breakage generally occurred due to dent deformation of the steel layer. The fiber breakage of CYH occurs near its center, as depicted in Figure 8a. These failure modes were mainly caused by large deformation in the center. Because of such deformation, the crack extended to the edge of the local dent first in the wrapped direction and then in the opposite direction.
The history of pressure applied to the cylinders, obtained from the hydrostatic tests, is depicted in Figure 9. Each applied pressure monotonically increased until it reached its maximum value, which is the collapse pressure of the corresponding sample. Each pressure fluctuation resembles a staggered slope because of the manual discontinuous operation of the water pump. The staggered horizontal and upward small segments indicate pressure holding and application processes, respectively. Once the maximum value is reached, the pressure on each cylinder suddenly decreases, suggesting the collapse of the cylinders. This decrease occurs because the volume of the steel cylinder is reduced due to deformation. That reveals the plasticity deformation of the steel layer occurred.
The experiments were highly repeatable. The collapse pressures of CYH1, CYH2, and CYH3 were 2.920, 2.819, and 2.820 MPa, respectively. The average collapse pressure and the maximum difference between the collapse pressure were 2.853 MPa and 3.46%, respectively. The collapse pressures of CYL1 and CYL2 were 3.232 and 2.964 MPa, with the average of and difference between these values being 3.098 MPa and 8.29%, respectively. On the basis of the modified stiffness matrix, the NASA SP-8700, ASME Code 2007 (RD-1172), and proposed formulas were used to calculate the collapse pressure of the hybrid cylinders with FC = 2.8 (Table 3). The NASA SP-8700 and ASME Code 2007 formulas exhibited errors of 25.5–36.9% and 9.7–14.2%, respectively. Although the error of the ASME Code 2007 formula was smaller than that of the NASA SP-8700 formula, both errors were unacceptably large. The proposed formula accounts for the interaction between the elastic buckling mode and the material failure mode. The average differences between the predictions of the proposed formula and the experimental results were 1.7% and 3.1% for CYH and CYL, respectively. Therefore, the proposed formula predicts collapse pressure for hybrid cylinders more accurately than the NASA SP-8700 and ASME Code 2007 formulas.
Nonlinear RIKS analysis was performed using the ABAQUS/Standard software to reveal the mechanism of collapse behavior of the samples. The finite element model of the sample was established in accordance with the sample’s real geometric shape, which was determined from scanning data. The multilayer shell-type approach was adopted to avoid interface problems associated with hybrid cylinders. This modeling approach has been successfully used to examine the collapse properties of steel–composite cylinders under axial or internal pressure [12,13,15]. As depicted in Figure 10, first, shell elements (S4R) present on the outer surface of the steel layer are assigned an outward offset section. The continuum shell elements (SC8R) are then generated using a mesh offset in the ABAQUS environment. The thickness of the inward offset is the nominal thickness of the CFRP layers (tc = (N − n) × tl). The elements of the steel and CFRP layers share the same nodes and are connected with each other (Figure 10). Two perfect surfaces are drawn and connected to the two ends of each cylinder to simulate the two heavy bungs. The elements of the cylinder and bungs share the same nodes; that is, they are connected. A unit of external pressure is applied to the external surface of each finite-element model. To prevent rigid body displacement, three-point-constraint boundary conditions were used in the analysis. These conditions do not overconstrain the problem because the pressure is equally applied. CCS2018 [49] recommends using similar boundary conditions for evaluating the buckling of various revolution shells under uniform external pressure, and these conditions have been widely used to evaluate the buckling of pressure shells [5].
For the case of the cylinder CYL1, the numerical results correlate satisfactorily with the experimental data. The numerically obtained equilibrium curve is presented in Figure 11. The curve shows an initial linear increase before a gradual decrease. The maximum von Mises equivalent stress of the steel layer at the maximum pressure point was equal to the material yield point, which suggests that nonlinear elastic–plastic buckling of the steel layer occurred under uniform external pressure. The collapse pressure of the cylinder CYL1, as determined using nonlinear RIKS analysis, is 3.024 MPa, 6.64% lower than that in the experimental data. The difference between the predictions of the derived formula and the numerical result was 5.32% for CYL1. Therefore, the derived formula could reasonably predict the collapse pressure of the hybrid cylinder. Moreover, the postcollapse mode at the end of the pressure curves in Figure 11 formed a local dent. These findings are identical to the experimental observations displayed in Figure 8b. Such a local dent is associated with the initial geometric imperfections and plasticity of the steel layer.
For further comparison, a linear eigenvalue analysis was performed using the ABAQUS/Standard software to examine the buckling pressure of the steel–composite hybrid cylinders, and the obtained numerical results were compared with the analytical results. The elastic buckling pressure of steel–composite hybrid cylinders were analytically calculated using the modified NASA SP-8700, ASME Code 2007, and proposed formulas for FC = 1. The length (L) and radius (Rs) of the cylinders were 320 and 80 mm, respectively. The ratio of the composite layer to the steel layer was 0.2, 0.8, or 1.6. The wrap sequences of the hybrid cylinders were [15/−15]4, [30/−30]4, [45/−45]4, [50/−50]4, [55/−55]4, [60/−60]4, [65/−65]4, [75/−75]4, and [85/−85]4. Elastic buckling modes obtained using ABAQUS are listed in Table 4. The analytical and numerical results for the buckling pressure of these hybrid cylinders are listed in Table 5, which also presents the wrap angles for these cylinders. The ratios of the buckling pressure calculated with the NASA SP-8007, ASME Code 2007, and proposed formulas to the numerical pressure values were 1.380–1.881, 0.359–0.514, and 0.745–1.030, respectively. For the proposed formula, the average ratio between the calculated and numerical pressure was 0.9, with the standard deviation being 0.09 and the standard error being 0.029 (Equation (31) [50]). Therefore, the proposed formula can accurately evaluate the buckling pressure of hybrid cylinders and is superior to the other two formulas in this task.
S E 2 = ( P C P a b a q u s 1 ) 2 N ( N 1 ) ,

3.2. Effects of the Wrap Angle, Thickness, and Length on the Collapse Pressure

The modified NASA SP-8700 ( P N A S A m o ), modified ASME Code 2007 ( P A S M E m o ), and proposed formula (PC) were used to determine the effects of the length, wrap angle, and number of layers on the collapse pressure of hybrid cylinders with FC = 2.8. The length (L)-to-radius (Rs) ratios of the cylinders were 2, 3, 4, 6, 8, and 10. The ratios between the thicknesses of the composite (tc) and steel (ts) layers were 0.2, 0.4, 0.8, 1.2, and 1.6. The wrap angle (θ) of the cross-ply composite layer of the hybrid cylinders was 15°–85°. The analytical results for the collapse pressure of these cylinders are listed in Table 6, Table 7, Table 8, Table 9 and Table 10. The error of the ASME formula, which was calculated using the expression ( P A S M E m o P t e s t ) / P t e s t , was lower than that of the NASA formula for L/RS values of 2, 3, and 4. The collapse pressures calculated with the ASME and proposed formulas were similar for medium-length and short cylinders.
The calculation results were obtained using the proposed formula for the effects of the wrap angle (θ) on the collapse pressure of the steel–composite hybrid cylinders with tc/ts values of 0.2, 0.8, and 1.6 (Figure 12). When tc/ts = 0.2, the collapse pressure of the hybrid cylinders was similar for θ values of 15°–85° (Figure 12a). For a tc/ts value of 0.8 or 1.6, as θ increased, the collapse pressure initially increased and then stabilized (Figure 12a,b). Therefore, the increase in θ can improve the collapse pressure, which can be increased by 113.09% (tc/ts = 1.6, L/RS = 2). Notably, the inflection points of the trends occur at θ = 55°–65°. The collapse pressure was maximized when L/RS = 2 and θ = 55°. The aforementioned results suggest that a wrap angle of ±55° can optimize the loading capacity of steel–composite hybrid cylinders. This angle is the classical wrap angle for a composite pressure shell and is widely used for composite cylinders subjected to external pressure [51,52]. These findings are mainly because the hoop stress is twice the axial stress of a cylinder subjected to external pressure.
For each length, the collapse pressure monotonically increased with an increase in tc/ts (ts = 1.5 mm, Figure 13). Notably, the rate of this increase reduced as L/RS increased and was affected by the wrap angle (θ). That is because the loading capacity of the longer cylinders considerably decreased. The thickness increasing plays a weak part in improving the loading capacity of the longer cylinders. For an L/RS value of 2, the maximum increases in the collapse pressure with the tc/ts value were 167.86% and 262.93% when θ was 15° and 85°, respectively. The collapse pressure also decreased monotonically as L/RS increased (Figure 14). This decrease was initially rapid but then became more gradual. These findings are mainly because the collapse resistance is considerably higher for shorter cylinders. The length plays a prominent part in improving the loading capacity of the shorter cylinders. The aforementioned trends are reasonably consistent with those obtained for steel-only cylinders [47].

4. Conclusions

To evaluate the collapse pressure of the steel–composite hybrid cylinders without excessive computational cost, a derived analytical formula was presented in the present work. The rationality of the derived formula was verified by the comparison with experimental and numerical results. The CYH and CYL hybrid cylinder comprises an outer composite layer and an inner steel layer. The wrap sequences of the composite layer were [55/−55]4 and [90/90/0/90/90/0/90/90] for CLH and CLY, respectively. Moreover, the effects of the wrap angle, thickness, and length on the collapse pressure of the hybrid cylinders were theoretically analyzed. The following primary conclusions were drawn:
(1)
The derived formula can determine the collapse pressure of steel-only cylinders with acceptable accuracy. The errors of this formula with respect to the results obtained in two verification experiments were 0.020–0.069 and 0.027–0.031. The minimum errors of the proposed formula in these experiments were 0.020 and 0.027, which were lower than the corresponding errors of other analytical formulas.
(2)
The experimental results obtained for the steel–composite hybrid cylinders were repeatable. The maximum difference between the experimental collapse pressure was 8.29%. These findings indicate that samples are manufactured and tested with good quality. The average difference between the collapse pressure calculated using the proposed formula and the experimental results was 1.7% and 3.1% for the CYH and CYL cylinders, respectively. The derived formula considered material failure and could reasonably predict the collapse pressure of the steel–composite hybrid cylinders.
(3)
An increase in wrap angle caused an increase in the collapse pressure (by up to 113.09%). Notably, the inflection points of the trends occur at a wrap angle ranging from 55° to 65°. The maximum collapse pressure was observed. The aforementioned results suggest that the loading capacity of steel–composite hybrid cylinders can be maximized under a wrap angle of ±55°. These findings are mainly because the hoop stress is twice the value of axial stress for cylinders under uniform pressure.
(4)
For each cylinder length, the collapse pressure monotonically increased as the thickness ratios of the composite layer to the steel layer increased. However, the rate of this increase was reduced as the length-to-radius ratio increased and was considerably affected by the wrap angle. That is because the loading capacity of the longer cylinders considerably decreases. The thickness increasing plays a weak part in improving the loading capacity of the longer cylinders. Moreover, as the length-to-radius ratio increased, the collapse pressure initially decreased rapidly and then decreased gradually. These findings are mainly because the collapse resistance is considerably higher for shorter cylinders. The length plays a prominent part in improving the loading capacity of the shorter cylinders.
Future studies should examine the effects of local damage propagation on the collapse properties of steel–composite hybrid cylinders by conducting hydrostatic tests and developing corresponding analytical models for these cylinders.

Author Contributions

Conceptualization, W.T. and J.Z.; methodology, J.Z.; software, X.Z.; validation, X.Z., Y.L., and M.Z.; investigation, M.Z.; data curation, X.Z. and Y.L.; writing—original draft preparation, X.Z.; writing—review and editing, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 52071160 and 52171258); Six Talent Climax Foundation of Jiangsu, China (grant number KTHY-068); and Graduate Research and Innovation Projects in Jiangsu Province, China (KYCX21_3498).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Magnucki, K.; Jasion, P.; Rodak, M. Strength and Buckling of an Untypical Dished Head of a Cylindrical Pressure Vessel. Int. J. Press. Vessel. Pip. 2018, 161, 17–21. [Google Scholar] [CrossRef]
  2. Magnucki, K.; Lewinski, J.; Cichy, R. Strength and Buckling Problems of Dished Heads of Pressure Vessels-Contemporary Look. J. Press. Vessel Technol. Trans. ASME 2018, 140, 041201. [Google Scholar] [CrossRef]
  3. Blachut, J. Buckling of Cylinders with Imperfect Length. In Proceedings of the ASME 2013 Pressure Vessels and Piping Conference, Paris, France, 14–18 July 2013; pp. 1–9. [Google Scholar]
  4. Błachut, J. Experimental Perspective on the Buckling of Pressure Vessel Components. Appl. Mech. Rev. 2014, 66, 30–33. [Google Scholar] [CrossRef]
  5. Blachut, J.; Magnucki, K. Strength, Stability, and Optimization of Pressure Vessels: Review of Selected Problems. Appl. Mech. Rev. 2008, 61, 0608011–06080133. [Google Scholar] [CrossRef]
  6. Zhang, J.; Zhu, Z.; Wang, F.; Zhao, X.; Zhu, Y. Buckling Behaviour of Double-Layer and Single-Layer Stainless Steel Cylinders under External Pressure. Thin-Walled Struct. 2021, 161, 107485. [Google Scholar] [CrossRef]
  7. Khamlichi, A.; Bezzazi, M.; Limam, A. Buckling of Elastic Cylindrical Shells Considering the Effect of Localized Axisymmetric Imperfections. Thin-Walled Struct. 2004, 42, 1035–1047. [Google Scholar] [CrossRef]
  8. Fatemi, S.M.; Showkati, H.; Maali, M. Experiments on Imperfect Cylindrical Shells under Uniform External Pressure. Thin-Walled Struct. 2013, 65, 14–25. [Google Scholar] [CrossRef]
  9. Ifayefunmi, O.; Fadzullah, S.H.S.M. Buckling Behaviour of Imperfect Axially Compressed Cylinder with an Axial Crack. Int. J. Automot. Mech. Eng. 2017, 14, 3837–3848. [Google Scholar] [CrossRef]
  10. Błachut, J. Buckling of Axially Compressed Cylinders with Imperfect Length. Comput. Struct. 2010, 88, 365–374. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Liu, Z.; Xin, J.; Wang, Y.; Zhang, C.; Zhang, Y. The Attenuation Mechanism of CFRP Repaired Corroded Marine Pipelines Based on Experiments and FEM. Thin-Walled Struct. 2021, 169, 108469. [Google Scholar] [CrossRef]
  12. Vukelic, G.; Vizentin, G.; Bakhtiari, R. Failure Analysis of a Steel Pressure Vessel with a Composite Wrap Repair Proposal. Int. J. Press. Vessel. Pip. 2021, 193, 104476. [Google Scholar] [CrossRef]
  13. Kr Singh, D.; Villamayor, A.; Hazra, A. Numerical and Experimental Analysis of Loctite Adhesive Composite Wrapping on En 10028 Steel Pipe. Mater. Today Proc. 2020, 44, 4158–4165. [Google Scholar] [CrossRef]
  14. Maali, M.; Kılıç, M.; Yaman, Z.; Ağcakoca, E.; Aydın, A.C. Buckling and Post-Buckling Behavior of Various Dented Cylindrical Shells Using CFRP Strips Subjected to Uniform External Pressure: Comparison of Theoretical and Experimental Data. Thin-Walled Struct. 2019, 137, 29–39. [Google Scholar] [CrossRef]
  15. Draidi, Z.; Bui, T.T.; Limam, A.; Tran, H.V.; Bennani, A. Buckling Behavior of Metallic Cylindrical Shell Structures Strengthened with CFRP Composite. Adv. Civ. Eng. 2018, 2018, 4231631. [Google Scholar] [CrossRef]
  16. Carvelli, V.; Panzeri, N.; Poggi, C. Buckling Strength of GFRP Under-Water Vehicles. Compos. Part B Eng. 2001, 32, 89–101. [Google Scholar] [CrossRef]
  17. Özbek, Ö. Axial and Lateral Buckling Analysis of Kevlar/Epoxy Fiber-Reinforced Composite Laminates Incorporating Silica Nanoparticles. Polym. Compos. 2021, 42, 1109–1122. [Google Scholar] [CrossRef]
  18. Ross, C.T.F. A Conceptual Design of an Underwater Vehicle. Ocean Eng. 2006, 33, 2087–2104. [Google Scholar] [CrossRef]
  19. Davies, P.; Choqueuse, D.; Bigourdan, B.; Chauchot, P. Composite Cylinders for Deep Sea Applications: An Overview. J. Press. Vessel Technol. Trans. ASME 2016, 138, 060904. [Google Scholar] [CrossRef]
  20. Wei, R.; Pan, G.; Jiang, J.; Shen, K.; Lyu, D. An Efficient Approach for Stacking Sequence Optimization of Symmetrical Laminated Composite Cylindrical Shells Based on a Genetic Algorithm. Thin-Walled Struct. 2019, 142, 160–170. [Google Scholar] [CrossRef]
  21. Teng, J.G.; Hu, Y.M. Behaviour of FRP-Jacketed Circular Steel Tubes and Cylindrical Shells under Axial Compression. Constr. Build. Mater. 2007, 21, 827–838. [Google Scholar] [CrossRef]
  22. Choqueuse, D.; Bigourdan, B.; Deuff, A.; Douchin, B.; Quetel, L. Hydrostatic Compression Behaviour of Steel-Composite Hybrid Tubes. In Proceedings of the 17th International Conference on Composite Materials, Edinburgh, UK, 27–31 July 2009. [Google Scholar]
  23. Vakili, M.; Showkati, H. Experimental and Numerical Investigation of Elephant Foot Buckling and Retrofitting of Cylindrical Shells by FRP. J. Compos. Constr. 2016, 20, 04015087. [Google Scholar] [CrossRef]
  24. Ventsel, E.S.; Krauthammer, T.; Carrera, E. Thin Plates and Shells: Theory: Analysis, and Applications; CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
  25. Ross, C. A Proposed Design Chart to Predict the Inelastic Buckling Pressures for Conical Shells Under Uniform External Pressure. Mar. Technol. 2007, 44, 77–81. [Google Scholar] [CrossRef]
  26. Society, China Classification, Specification for the Entry and Construction of Diving Systems and Submersible; People Traffic Press: Beijing, China, 2018.
  27. American Bureau of Shipping. Underwater Vehicles, Systems and Hyperbaric Facilities; American Bureau of Shipping: Houston, TX, USA, 2012. [Google Scholar]
  28. Buckling of Thin-Walled Circular Cylinders(NASA/SP-8007); National Aeronautics and Space Administration: Moffett Field, CA, USA, 2020.
  29. Section: X Fiber-Reinforced Plastic Pressure Vessel; ASME Boiler and Pressure Vessel Code; The American Society of Mechanical Engineers: New York, NY, USA, 2017.
  30. Liu, J.; Yu, B.; Zhou, Y.; Zhang, Y.; Duan, M. The Buckling of Spherical-Cylindrical Composite Shells by External Pressure. Compos. Struct. 2021, 265, 113773. [Google Scholar] [CrossRef]
  31. Rosenow, M.W.K. Wind Angle Effects in Glass Fibre-Reinforced Polyester Filament Wound Pipes. Composites 1984, 15, 144–152. [Google Scholar] [CrossRef]
  32. Imran, M.; Shi, D.; Tong, L.; Elahi, A.; Waqas, H.M.; Uddin, M. Multi-Objective Design Optimization of Composite Submerged Cylindrical Pressure Hull for Minimum Buoyancy Factor and Maximum Buckling Load Capacity. Def. Technol. 2021, 17, 1190–1206. [Google Scholar] [CrossRef]
  33. Imran, M.; Shi, D.; Tong, L.; Waqas, H.M. Design Optimization of Composite Submerged Cylindrical Pressure Hull Using Genetic Algorithm and Finite Element Analysis. Ocean Eng. 2019, 190, 106443. [Google Scholar] [CrossRef]
  34. Imran, M.; Shi, D.; Tong, L.; Elahi, A.; Uddin, M. On the Elastic Buckling of Cross-Ply Composite Closed Cylindrical Shell under Hydrostatic Pressure. Ocean Eng. 2021, 227, 108633. [Google Scholar] [CrossRef]
  35. Lopatin, A.V.; Morozov, E.V. Buckling of Composite Cylindrical Shells with Rigid End Disks under Hydrostatic Pressure. Compos. Struct. 2017, 173, 136–143. [Google Scholar] [CrossRef]
  36. Lopatin, A.V.; Morozov, E.V. Buckling of the Composite Sandwich Cylindrical Shell with Clamped Ends under Uniform External Pressure. Compos. Struct. 2015, 122, 209–216. [Google Scholar] [CrossRef]
  37. Bai, Y.; Xu, W.; Cheng, P.; Wang, N.; Ruan, W. Behaviour of Reinforced Thermoplastic Pipe (RTP) under Combined External Pressure and Tension. Ships Offshore Struct. 2014, 9, 464–474. [Google Scholar] [CrossRef]
  38. Bai, Y.; Ruan, W.; Cheng, P.; Yu, B.; Xu, W. Buckling of Reinforced Thermoplastic Pipe (RTP) under Combined Bending and Tension. Ships Offshore Struct. 2014, 9, 525–539. [Google Scholar] [CrossRef]
  39. Bai, Y.; Yuan, S.; Tang, J.; Qiao, H.; Cheng, P.; Cao, Y.; Wang, N. Behaviour of Reinforced Thermoplastic Pipe under Combined Bending and External Pressure. Ships Offshore Struct. 2015, 10, 575–586. [Google Scholar] [CrossRef]
  40. Messager, T. Buckling of Imperfect Laminated Cylinders under Hydrostatic Pressure. Compos. Struct. 2001, 53, 301–307. [Google Scholar] [CrossRef]
  41. GB/T 3280-2015; Cold Rolled Stainless Steel Plate, Sheet and Strip. Chinese Standard Instutute: Beijing, China, 2016.
  42. ASTM A959-16; Standard Guide for Specifying Harmonized Standard Grade Composition. American Society of Testing Materials: West Conshohoken, PA, USA, 2016.
  43. Błachut, J. On Elastic-Plastic Buckling of Cones. Thin-Walled Struct. 2011, 49, 45–52. [Google Scholar] [CrossRef]
  44. Błachut, J. Buckling of Externally Pressurized Steel Toriconical Shells. Int. J. Press. Vessel. Pip. 2016, 144, 25–34. [Google Scholar] [CrossRef]
  45. Jaspart, J.P. Extending of the Merchant-Rankine Formula for the Assessment of the Ultimate Load of Frames with Semi-Rigid Joints. J. Constr. Steel Res. 1988, 11, 283–312. [Google Scholar] [CrossRef]
  46. Shen, K.C.; Pan, G. Buckling and Strain Response of Filament Winding Composite Cylindrical Shell Subjected to Hydrostatic Pressure: Numerical Solution and Experiment. Compos. Struct. 2021, 276, 114534. [Google Scholar] [CrossRef]
  47. Zhu, Y.; Dai, Y.; Ma, Q.; Tang, W. Buckling of Externally Pressurized Cylindrical Shell: A Comparison of Theoretical and Experimental Data. Thin-Walled Struct. 2018, 129, 309–316. [Google Scholar] [CrossRef]
  48. Ross, C.T.F.; Little, A.P.F.; Haidar, Y.; Waheeb, A.A. Buckling of Carbon/Glass Composite Tubes under Uniform External Hydrostatic Pressure. Strain 2011, 47, 156–174. [Google Scholar] [CrossRef]
  49. Rules for Classification of Diving Systems and Submersibles; China Classification Society: Beijing, China, 2018; ISBN 0086010581128.
  50. Cho, Y.S.; Oh, D.H.; Paik, J.K. An Empirical Formula for Predicting the Collapse Strength of Composite Cylindrical-Shell Structures under External Pressure Loads. Ocean Eng. 2019, 172, 191–198. [Google Scholar] [CrossRef]
  51. Carroll, M.; Ellyin, F.; Kujawski, D.; Chiu, A.S. The Rate-Dependent Behaviour of ±55° Filament-Wound Glass-Fibre/Epoxy Tubes under Biaxial Loading. Compos. Sci. Technol. 1995, 55, 391–403. [Google Scholar] [CrossRef]
  52. Liu, W.; Soden, P.D.; Kaddour, A.S. Design of End Plugs and Specimen Reinforcement for Testing ± 55° Glass/Epoxy Composite Tubes under Biaxial Compression. Comput. Struct. 2005, 83, 976–988. [Google Scholar] [CrossRef]
Figure 1. Geometry of a steel–composite cylindrical shell.
Figure 1. Geometry of a steel–composite cylindrical shell.
Metals 12 01591 g001
Figure 2. Cylindrical shell under hydrostatic pressure.
Figure 2. Cylindrical shell under hydrostatic pressure.
Metals 12 01591 g002
Figure 3. Stacking sequence of the steel–composite cylindrical shell.
Figure 3. Stacking sequence of the steel–composite cylindrical shell.
Metals 12 01591 g003
Figure 4. Stress of a two-dimensional single layer: (a) fiber-directional stress; (b) in-plane transversal stress; (c) in-plane shear stress.
Figure 4. Stress of a two-dimensional single layer: (a) fiber-directional stress; (b) in-plane transversal stress; (c) in-plane shear stress.
Metals 12 01591 g004
Figure 5. Layer orientation (θk) and internal forces of the composite layer of a steel–composite cylinder.
Figure 5. Layer orientation (θk) and internal forces of the composite layer of a steel–composite cylinder.
Metals 12 01591 g005
Figure 6. Internal force of a single steel layer (Δts = tl) of a steel–composite cylinder.
Figure 6. Internal force of a single steel layer (Δts = tl) of a steel–composite cylinder.
Metals 12 01591 g006
Figure 7. Experimental flow: (a) flowchart of cylinder fabrication; (b) samples; (c) sample testing; (d) flanges removed using heat gun.
Figure 7. Experimental flow: (a) flowchart of cylinder fabrication; (b) samples; (c) sample testing; (d) flanges removed using heat gun.
Metals 12 01591 g007
Figure 8. Collapse modes of the hybrid cylinders with removed rigid bungs: (a) CYH (outer composite [±55]4) and (b) CYL (outer composite [90/90/0/90/90/0/90/90]).
Figure 8. Collapse modes of the hybrid cylinders with removed rigid bungs: (a) CYH (outer composite [±55]4) and (b) CYL (outer composite [90/90/0/90/90/0/90/90]).
Metals 12 01591 g008
Figure 9. Pressure–time curves obtained through the hydrostatic testing of the hybrid cylinders.
Figure 9. Pressure–time curves obtained through the hydrostatic testing of the hybrid cylinders.
Metals 12 01591 g009
Figure 10. Finite-element model of the hybrid cylinders.
Figure 10. Finite-element model of the hybrid cylinders.
Metals 12 01591 g010
Figure 11. Applied pressure versus collapse point displacement for the cylinder CYL1.
Figure 11. Applied pressure versus collapse point displacement for the cylinder CYL1.
Metals 12 01591 g011
Figure 12. Effect of the wrap angle (θ) of the composite layer on the collapse pressure of the hybrid cylinders: (a) tc/ts = 0.2, (b) tc/ts = 0.8, and (c) tc/ts = 1.6.
Figure 12. Effect of the wrap angle (θ) of the composite layer on the collapse pressure of the hybrid cylinders: (a) tc/ts = 0.2, (b) tc/ts = 0.8, and (c) tc/ts = 1.6.
Metals 12 01591 g012aMetals 12 01591 g012b
Figure 13. Effect of the thickness (tc) of the composite layer on the collapse pressure of the hybrid cylinders: (a) θ = 15° and (b) θ = 85°.
Figure 13. Effect of the thickness (tc) of the composite layer on the collapse pressure of the hybrid cylinders: (a) θ = 15° and (b) θ = 85°.
Metals 12 01591 g013
Figure 14. Effect of the length-to-radius ratio (L/Rs,) of hybrid cylinders on their collapse pressure: (a) θ = 15° and (b) θ = 55°.
Figure 14. Effect of the length-to-radius ratio (L/Rs,) of hybrid cylinders on their collapse pressure: (a) θ = 15° and (b) θ = 55°.
Metals 12 01591 g014
Table 1. Collapse pressure (MPa) of steel-only cylinders (hybrid cylinder with tc = 0) obtained using analytical formulas and in experiment I.
Table 1. Collapse pressure (MPa) of steel-only cylinders (hybrid cylinder with tc = 0) obtained using analytical formulas and in experiment I.
Case I [47] P C C S P A B S P t e s t [47] P N A S A m o P A S M E m o P C Difference
CCS-TestABS-TestC-Test
L/R = 1.04.5003.6004.8005.0294.2284.8980.0630.2500.020
L/R = 1.53.5202.4503.2703.4462.6703.3860.0760.2510.035
L/R = 2.02.9001.8702.9702.8052.0032.7650.0240.3700.069
L/R = ratio of the cylinder length to the inner radius; PCCS and PABS were calculated using the CCS2013 [26] and ABS2012 [27] formulas, respectively; CCS-test = |PCCSPtest|/Ptest; ABS-test = |PABSPtest|/Ptest; C-test = |PCPtest|/Ptest.
Table 2. Collapse pressure (MPa) of the steel-only cylinders (hybrid cylinder with tc = 0) obtained using analytical formulas and in experiment II.
Table 2. Collapse pressure (MPa) of the steel-only cylinders (hybrid cylinder with tc = 0) obtained using analytical formulas and in experiment II.
Case II [6] P V - K P R o s s P t e s t [6] P N A S A m o P A S M E m o P C Difference
V-K-TestRoss-TestC-Test
SL14.8024.9494.6214.8333.8504.7470.0390.0710.027
SL24.9585.1134.6834.9163.9124.8270.0590.0920.031
PV-K and PRoss were calculated using the Venstel–Krauthaer [24] and Ross [25] formulas, respectively; V-K-test = |PV-KPtest|/Ptest; Ross-test = |PRossPtest|/Ptest; C-test = |PCPtest|/Ptest.
Table 3. Experimental and analytical collapse pressure (MPa) of the steel–composite hybrid cylinders (tc-nominal/ts-nominal = 0.8).
Table 3. Experimental and analytical collapse pressure (MPa) of the steel–composite hybrid cylinders (tc-nominal/ts-nominal = 0.8).
Sampleθ (°) P t e s t P t e s t a v P N A S A m o P A S M E m o P C Difference
NASA-TestASME-TestC-Test
CYH1[±55]42.9202.8533.5812.4482.8050.2550.1420.017
CYH22.819
CYH32.820
CYL1[90/90/
0/90/90
/0/90/90]
3.2323.0984.2412.7983.1940.3690.0970.031
CYL22.964
θ = wrap angle of the outer composite layer of a hybrid cylinder.
Table 4. Elastic buckling modes obtained using ABAQUS.
Table 4. Elastic buckling modes obtained using ABAQUS.
tc/tsθ (°)
153045505560657585
0.2n = 4
Metals 12 01591 i001
n = 4
Metals 12 01591 i002
n = 4
Metals 12 01591 i003
n = 4
Metals 12 01591 i004
n = 4
Metals 12 01591 i005
n = 4
Metals 12 01591 i006
n = 4
Metals 12 01591 i007
n = 4
Metals 12 01591 i008
n = 4
Metals 12 01591 i009
0.8n = 4
Metals 12 01591 i010
n = 4
Metals 12 01591 i011
n = 4
Metals 12 01591 i012
n = 4
Metals 12 01591 i013
n = 4
Metals 12 01591 i014
n = 4
Metals 12 01591 i015
n = 4
Metals 12 01591 i016
n = 4
Metals 12 01591 i017
n = 4
Metals 12 01591 i018
1.6n = 4
Metals 12 01591 i019
n = 4
Metals 12 01591 i020
n = 3
Metals 12 01591 i021
n = 3
Metals 12 01591 i022
n = 3
Metals 12 01591 i023
n = 3
Metals 12 01591 i024
n = 3
Metals 12 01591 i025
n = 3
Metals 12 01591 i026
n = 3
Metals 12 01591 i027
n = circumferential waves.
Table 5. Numerical and analytical buckling pressure (MPa).
Table 5. Numerical and analytical buckling pressure (MPa).
tc/tsθ (°) P N A S A m o P A S M E m o P C P a b a q u s P N A S A m o / P a b a q u s P A S M E m o / P a b a q u s P C / P a b a q u s
0.2154.2661.6402.6272.8651.4890.3820.917
304.4071.6582.8232.9181.5110.3790.968
454.6411.7453.0153.0391.5270.3830.992
504.7301.7913.0473.0911.5310.3860.986
554.8201.8423.0593.1451.5330.3910.973
604.9051.8963.0563.2001.5330.3950.955
654.9831.9483.0423.2521.5320.3990.935
755.1042.0363.0073.3431.5270.4060.900
855.1672.0852.9853.3971.5210.4090.879
0.8156.1062.1753.8053.7981.6080.3821.002
307.3332.3464.4654.3371.6910.3611.030
459.3443.0185.4225.3971.7310.3731.005
509.6723.3245.5465.8131.6640.3810.954
559.9573.6385.6116.2231.6000.3900.902
6010.1923.9385.6366.6081.5420.3970.853
6510.3744.2075.6416.9511.4920.4030.812
7510.6004.6175.6407.4711.4190.4120.755
8510.6894.8325.6407.7441.3800.4160.728
1.61511.6393.9666.4416.3081.8450.4191.021
3015.8724.5437.9628.4401.8810.3590.943
4520.6276.4559.79611.8661.7380.3630.826
5022.0397.33110.20712.4051.7770.3940.823
5523.2258.18510.42912.8821.8030.4240.810
6024.1618.96210.53513.6291.7730.4380.773
6524.8509.63010.59713.6311.8230.4710.777
7525.60310.60110.68714.1231.8130.5000.757
8525.69311.08910.71514.3751.7870.5140.745
Table 6. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.2.
Table 6. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.2.
L R s θ (°) P N A S A m o P A S M E m o P C Difference L R s P N A S A m o P A S M E m o P C Difference
NASA-CASME-CNASA-CASME-C
2153.0732.1872.1200.4500.03261.0290.7290.8940.1510.185
303.1892.2112.2680.4060.025 1.0600.7370.9340.1350.211
453.3122.3272.3920.3850.027 1.1100.7760.9840.1290.211
503.3482.3882.4070.3910.008 1.1300.7960.9980.1320.202
553.3812.4572.4090.4040.020 1.1490.8191.0110.1370.190
603.4122.5282.4010.4210.053 1.1680.8421.0210.1440.175
653.4392.5972.3880.4400.088 1.1850.8661.0290.1520.159
753.4772.7142.3570.4750.151 1.2120.9051.0400.1660.130
853.4962.7812.3390.4950.189 1.2260.9271.0450.1740.113
3152.0611.4581.5840.3020.07980.8070.5470.7220.1180.242
302.1391.4741.6810.2720.123 0.8370.5530.7570.1070.269
452.2621.5521.7920.2630.134 0.8940.5820.8100.1040.281
502.2871.5921.8050.2670.118 0.9170.5970.8280.1070.279
552.3121.6381.8120.2760.096 0.9400.6140.8450.1120.274
602.3341.6851.8120.2880.070 0.9640.6320.8610.1190.266
652.3541.7311.8090.3010.043 0.9860.6320.8750.1260.278
752.3851.8091.7990.3260.006 1.0160.6790.8920.1390.240
852.4001.8531.7910.3400.034 1.0200.6950.8910.1440.220
4151.5241.0941.2460.2230.122100.6450.4370.5900.0940.258
301.5741.1061.3110.2000.157 0.6620.4420.6100.0840.275
451.6571.1641.3900.1920.163 0.6870.4650.6370.0800.269
501.6891.1941.4110.1970.154 0.6970.4780.6450.0810.259
551.7211.2281.4280.2060.140 0.7070.4910.6520.0840.246
601.7521.2641.4410.2160.123 0.7160.5050.6580.0880.232
651.7801.2991.4490.2280.104 0.7250.5190.6630.0930.217
751.8231.3571.4590.2490.070 0.7380.5430.6700.1010.190
851.8451.3901.4630.2610.050 0.7450.5560.6740.1050.175
Table 7. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.4.
Table 7. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.4.
L R s θ (°) P N A S A m o P A S M E m o P C Difference L R s P N A S A m o P A S M E m o P C Difference
NASA-CASME-CNASA-CASME-C
2153.4562.4532.5750.3420.04761.1330.7841.0180.1120.230
303.7422.4562.8180.3280.129 1.2160.8091.0990.1060.264
454.0392.7313.0770.3130.112 1.3550.9101.2260.1050.258
504.1362.8823.1380.3180.081 1.4070.9611.2340.1400.222
554.2233.0453.1790.3280.042 1.4581.0151.2700.1480.201
604.2993.2093.2060.3410.001 1.5071.0701.3030.1570.179
654.3593.3633.2230.3530.044 1.5501.1211.3290.1660.157
754.4343.6133.2410.3680.115 1.6151.2041.3670.1810.119
854.4643.7513.2490.3740.154 1.6471.2501.3850.1890.097
3152.2851.5671.8630.2260.15980.8800.5880.7920.1110.258
302.4921.6182.0460.2180.209 0.9600.6070.8680.1060.301
452.7201.8212.2470.2110.190 1.1100.6831.0000.1110.317
502.7901.9222.2970.2150.163 1.1700.7211.0480.1170.312
552.8552.0302.3360.2220.131 1.1880.7611.0600.1200.282
602.9122.1392.3660.2310.096 1.2030.8021.0690.1250.250
652.9612.2422.3890.2400.061 1.2170.8411.0760.1310.219
753.0292.4092.4210.2510.005 1.2360.9031.0860.1390.168
853.0612.5002.4360.2560.026 1.2460.9381.0900.1430.140
4151.6831.1761.4420.1670.185100.7100.4700.6510.0900.278
301.8191.2131.5690.1590.227 0.7550.4850.6970.0830.303
452.0461.3661.7660.1580.227 0.8270.5460.7640.0820.285
502.1311.4411.8310.1640.213 0.8540.5760.7870.0850.267
552.2151.5231.8890.1720.194 0.8800.6090.8080.0890.246
602.2921.6041.9400.1820.173 0.9050.6420.8270.0940.224
652.3611.6821.9820.1910.152 0.9270.6730.8430.0990.202
752.4491.8072.0350.2030.112 0.9600.7230.8670.1080.166
852.4621.8752.0410.2060.081 0.9770.7500.8790.1120.146
Table 8. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.8.
Table 8. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 0.8.
L R s θ (°) P N A S A m o P A S M E m o P C Difference L R s P N A S A m o P A S M E m o P C Difference
NASA-CASME-CNASA-CASME-C
2154.6973.1133.2060.4650.02961.4450.9061.1900.2150.239
305.4843.3353.7040.4800.100 1.7151.0431.3990.2250.255
456.5064.0244.3260.5040.070 2.1711.3411.7340.2520.226
506.8204.4324.4730.5250.009 2.3381.4771.8420.2700.198
557.0864.8504.5680.551 0.062 2.4991.6171.9350.2920.164
607.2025.2504.5840.5710.145 2.6451.7502.0120.3150.130
657.2265.6094.5600.5840.230 2.7701.8702.0730.3360.098
757.2126.1564.5120.5990.364 2.8132.0522.0830.3500.015
857.1826.4434.4840.6020.437 2.8122.1482.0770.3530.034
3153.0112.0162.0800.4470.03181.1220.7160.9620.1670.255
303.5522.1062.4220.4670.130 1.3710.7821.1620.1800.327
454.2452.6832.8430.4930.057 1.6691.0061.3980.1940.280
504.4742.9552.9510.5160.001 1.7381.1081.4470.2000.234
554.6793.2343.0260.5460.068 1.8001.2131.4880.2100.185
604.8503.5003.0750.5770.138 1.8561.3131.5200.2210.137
654.9843.7393.1060.6050.204 1.9031.4021.5460.2310.093
755.1464.1043.1370.6410.308 1.9671.5391.5800.2450.026
855.2044.2953.1460.6540.365 1.9971.6111.5960.2510.009
4152.1811.4501.6470.3240.119100.8940.5730.7890.1330.274
302.6191.5641.9490.3440.197 1.0420.6260.9170.1370.317
453.3372.0122.4050.3880.163 1.2850.8051.1180.1490.280
503.4542.2162.4700.3990.103 1.3730.8861.1850.1580.252
553.5562.4252.5130.4150.035 1.4580.9701.2460.1700.221
603.6402.6252.5400.4330.034 1.5351.0501.2980.1830.191
653.7052.8052.5560.4500.097 1.6021.1221.3410.1940.163
753.7863.0782.5730.4710.196 1.6991.2311.4020.2120.122
853.8183.2212.5800.4800.249 1.7461.2891.4320.2190.100
Table 9. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 1.2.
Table 9. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 1.2.
L R s θ (°) P N A S A m o P A S M E m o P C Difference L R s P N A S A m o P A S M E m o P C Difference
NASA-CASME-CNASA-CASME-C
2156.6274.0624.3020.5400.05661.9361.2321.6720.1580.263
308.3324.6115.1830.6070.110 2.5191.4162.1280.1840.335
4510.3175.9836.2390.6540.041 3.5181.9942.8770.2230.307
5010.6346.7346.3870.6650.054 3.8252.2453.0870.2390.273
5510.8397.4816.4420.6830.161 3.8942.4943.1270.2450.203
6010.9388.1736.4370.6990.270 3.9432.7243.1490.2520.135
6510.9538.7776.4040.7100.371 3.9782.9263.1620.2580.075
7510.8389.2776.3260.7130.466 4.0153.2223.1760.2640.015
8510.7109.6676.2760.7070.540 4.0273.3743.1820.2660.060
3154.1542.6693.1030.3390.14081.5250.9241.3560.1240.318
305.2072.9833.7740.3800.210 2.0031.0621.7480.1460.392
456.7663.9894.7360.4290.158 2.4931.4962.1530.1580.305
507.2664.4904.9970.4540.101 2.6521.6842.2750.1660.260
557.5474.9875.1160.4750.025 2.7971.8702.3780.1760.214
607.6045.4485.1160.4860.065 2.9222.0432.4620.1870.170
657.6215.8515.1000.4940.147 3.0252.1942.5290.1960.132
757.5946.4455.0640.5000.273 3.1662.4172.6200.2080.078
857.5576.7475.0430.4990.338 3.2292.5302.6620.2130.050
4152.9821.9612.3980.2430.182101.1730.7391.0710.0960.310
303.9272.1413.0530.2860.299 1.4950.8501.3480.1090.370
454.9802.9923.7860.3150.210 2.0321.1971.8000.1290.335
505.2513.3673.9530.3280.148 2.2241.3471.9530.1390.310
555.4813.7404.0750.3450.082 2.4061.4962.0890.1510.284
605.6654.0864.1590.3620.017 2.5691.6352.2060.1650.259
655.8044.3884.2160.3760.041 2.7091.7552.3040.1760.238
755.9624.8344.2820.3920.129 2.9091.9332.4410.1910.208
856.0145.0614.3060.3970.175 3.0012.0242.5050.1980.192
Table 10. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 1.6.
Table 10. Collapse pressure (MPa) of hybrid cylinders obtained using analytical formulas when tc/ts = 1.6.
L R s θ (°) P N A S A m o P A S M E m o P C Difference L R s P N A S A m o P A S M E m o P C Difference
NASA-CASME-CNASA-CASME-C
2159.3685.4075.6790.6500.04862.6481.6962.2370.1840.242
3012.5156.3687.0170.7830.092 3.6931.9422.9990.2310.353
4515.3328.7028.4160.8220.034 5.1312.8694.0240.2750.287
5015.9209.7758.6650.8370.128 5.3243.2584.1590.2800.217
5516.28210.9138.7530.8600.247 5.4763.6384.2470.2890.144
6016.43211.9498.7420.8800.367 5.5903.9834.3030.2990.074
6516.41412.8418.6910.8890.478 5.6744.2804.3400.3070.014
7516.10214.1348.5730.8780.649 5.7714.7124.3890.3150.073
8515.77914.7868.4870.8590.742 5.8074.9294.4120.3160.117
3155.6843.5724.0770.3940.12482.1291.2251.8550.1480.339
307.6124.1595.1560.4770.193 2.7671.4562.3580.1730.383
4510.1825.7386.5870.5460.129 3.6812.1523.0750.1970.300
5010.5246.5166.7750.5530.038 3.9782.4443.2900.2090.257
5510.7467.2756.8550.5680.061 4.2452.7283.4680.2240.213
6010.8657.9666.8690.5820.160 4.4762.9873.6110.2400.173
6510.9058.5606.8570.5900.248 4.6663.2103.7250.2530.138
7510.8569.4236.8180.5920.382 4.9233.5343.8810.2690.090
8510.7799.8576.7930.5870.451 5.0363.6963.9520.2740.065
4154.1572.6443.2270.2880.181101.5700.9801.4160.1090.307
305.6693.0294.1840.3550.276 2.1471.1651.8930.1340.384
457.3674.3035.2810.3950.185 3.1231.7212.6750.1670.357
507.8714.8875.5670.4140.122 3.4681.9552.9330.1820.334
558.2955.4575.7680.4380.054 3.7922.1833.1590.2000.309
608.6295.9755.9020.4620.012 4.0822.3903.3500.2190.287
658.8756.4205.9950.4800.071 4.2122.5683.4300.2280.251
759.1447.0676.1010.4990.158 4.1022.8273.3520.2240.157
859.1767.3936.1190.5000.208 4.0452.9573.3150.2200.108
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zuo, X.; Tang, W.; Zhang, J.; Li, Y.; Zhan, M. Collapse of Externally Pressurized Steel–Composite Hybrid Cylinders: Analytical Solution and Experimental Verification. Metals 2022, 12, 1591. https://doi.org/10.3390/met12101591

AMA Style

Zuo X, Tang W, Zhang J, Li Y, Zhan M. Collapse of Externally Pressurized Steel–Composite Hybrid Cylinders: Analytical Solution and Experimental Verification. Metals. 2022; 12(10):1591. https://doi.org/10.3390/met12101591

Chicago/Turabian Style

Zuo, Xinlong, Wenxian Tang, Jian Zhang, Yongsheng Li, and Ming Zhan. 2022. "Collapse of Externally Pressurized Steel–Composite Hybrid Cylinders: Analytical Solution and Experimental Verification" Metals 12, no. 10: 1591. https://doi.org/10.3390/met12101591

APA Style

Zuo, X., Tang, W., Zhang, J., Li, Y., & Zhan, M. (2022). Collapse of Externally Pressurized Steel–Composite Hybrid Cylinders: Analytical Solution and Experimental Verification. Metals, 12(10), 1591. https://doi.org/10.3390/met12101591

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop