Next Article in Journal
Leveraging the Advantages of Additive Manufacturing to Produce Advanced Hybrid Composite Structures for Marine Energy Systems
Next Article in Special Issue
Process Prediction for Compound Screws by Using Virtual Measurement and Recognizable Performance Evaluation
Previous Article in Journal
Index of Body Inflammation for Maxillofacial Surgery Purpose-to Make the Soluble Urokinase-Type Plasminogen Activator Receptor Serum Level Independent on Patient Age
Previous Article in Special Issue
Detecting Underground Geospatial Features for Protection of Rights in 3D Space: Korean Cases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Buckling of Rectangular Composite Pipes under Torsion

1
Department of Mechanical Engineering, Kanagawa University, Rokkakubashi, Kanagawa-ku, Yokohama 221-8686, Japan
2
East Operation Group, Mitsubishi Electric TOKKI Systems Corporation, Kamimachiya, Kamakura 247-0065, Japan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(3), 1342; https://doi.org/10.3390/app11031342
Submission received: 9 December 2020 / Revised: 24 January 2021 / Accepted: 30 January 2021 / Published: 2 February 2021
(This article belongs to the Special Issue Selected Papers from IMETI 2020)

Abstract

:

Featured Application

Thin-walled rectangular composite pipes, because of their light weight, high stiffness, and strength, can be applied to ultra-lightweight airplanes, satellites, spacecraft, automobiles, and civil structures.

Abstract

The numerical buckling load of rectangular composite pipes under torsional load was derived by using the energy method. The authors found no available simple design method or chart for the buckling loads of rectangular composite pipes, which are often used airplanes, spacecraft, and other lightweight structures, through their involvement in a Mars exploration airplane project. Thus, numerical results were obtained for length-to-width ratios (l/b) from 1 to 20, width-to-height ratios (h/b) from 1 to 6, and [0/90] layer ratios (r) from 0 to 1, which means [(0/90)r,(±45)1-r]s. The layups were assumed to be symmetric, and tension-bending, torsion-bending, and tension-shear coupling stiffnesses were ignored. To establish a simple design method, a closed-form polynomial equation for the buckling load factor was derived by minimizing the weighted residuals of the safe and non-safe side errors, which were obtained by comparing the derived numerical results with the polynomial equations. As a result, the errors of the polynomial equation for the buckling load factor were 4.95% for the non-safe side and 12.4% for the safe side. The errors are sufficiently good for preliminary design use and for parametric design studies and optimization.

1. Introduction

The Mars exploration airplane project [1,2] is a joint project of JAXA (Japan Aerospace Exploration Agency) and Japanese universities. It is intended to be the first extraterrestrial atmospheric flight exploration. Aerial exploration can avoid complex terrain where a rover cannot explore and can also explore where a satellite cannot observe, like the sides of a valley. An image of the Mars exploration airplane in flight is shown in Figure 1. To build the airplane, the authors designed a wing to be tested at high altitude, where the density of the air is close to the density of the Martian atmosphere. A rectangular pipe, which can easily align the attack angle of the wing, was used. In addition, because the density of Martian atmosphere is 1/100th that of the Earth, the airplane must be ultra-lightweight; thus, the rectangular pipe was designed using thin carbon fiber reinforced plastics (CFRP). Thin structures usually fail not for lack of material strength or delamination but through buckling load because the buckling load of thin structure lower than material failure load. The buckling load can be increased by careful design of the layup sequence of the composite layers; however, no available simple design equation or chart were found in spite of there being many studies on composite structures [3,4,5,6,7,8,9,10]. The buckling load under shear force has been studied; however, those studies were not for composite rectangular tubes, but rather for long anisotropic plates [11,12]. An approximate expression for the torsional buckling load was derived by Omidvari and Hematiyan [13]; however, the solution is for an isotropic tube. Banks and Rhodes investigated postbuckling behavior of composite box sections [14]. Loughlan analyzed buckling of composite stiffened box sections under compression and bending by finite strip method [15]. Vo and Lee conducted numerical analysis for flexural–torsional buckling of thin-walled composite box beams by finite element method [16]. Moreover, while a comprehensive design handbook for the buckling of composite structures [17] has been published, it does not contain a method for the torsional buckling load of composite rectangular tubes. Given this state of affairs, it was decided that a numerical analysis of the torsional buckling load of cantilever rectangular orthotropic pipes should be conducted by using the energy method including a consideration of layup anisotropy. In addition, to easily use this analysis for design, a simple closed-form polynomial equation for the buckling load factor was derived by minimizing the weighted residuals of the safe and non-safe side errors of the derived numerical results. Although a numerical analysis like a finite element analysis (FEA) can determine the buckling load of anisotropic structures, a simple closed-form design equation is especially useful for the initial design and optimum design of the whole system that require a huge number of iterations for analysis of not only the structure but also its aerodynamics, guidance, and control.

2. Numerical Calculation

2.1. Energy Method

To use a rectangular thin-walled pipe as a wing spar, a buckling analysis using the energy method was conducted on the cantilevered rectangular pipe under torque T shown in Figure 2. The rectangular pipe consisted of two vertical plates and two horizontal plates. The vertical and horizontal plates were assumed to be simply supported. The out-plane displacements v and w are described as follows:
w = m = 1 n = 1 B m n sin m π x l sin n π y b
v = m = 1 n = 1 H m n sin m π x l sin n π z h
Note that the study discusses about local wall buckling of rectangular pipe, not global buckling of long column. Thus, warping is not included. As in reference [18], the plates’ edge angles were assumed to be constant, as shown in Figure 3.
In this case, the out-plane displacements v and w are described using Φmn as follows:
w = m = 1 n = 1 b h Φ m n sin m π x l sin n π y b
v = m = 1 n = 1 h b Φ m n sin m π x l sin n π z h
where n = 1, 2,…; m = 1, 2,…; and
B m n b = H m n h = Φ m n b h
Φmn is a common coefficient to determine amplitudes of w and v and Equation (5) is the condition of constant the plates’ edge angles. This assumption will be validated later. The total potential energy Π can be expressed in terms of the strain energy U and work of the external force W as
Π = U W
The contributions of a single horizontal plate (in x-y plane) to U and W are written as
U p = V 1 2 ( σ x ε x + σ y ε y + τ x y γ x y ) d V
W p = t τ c r w x w y d x d y
For vertical plates, y replaced z in Equations (7) and (8) and related equations appearing hereafter. For horizontal plates (in x-y plane), by considering the bending strain due to the out-plane displacement w, the strains can be described as follows:
ε x = z 2 w x 2 ,   ε y = z 2 w y 2 ,   γ x y = 2 z 2 w x y
Considering Hooke’s Law for anisotropic layers
[ σ x σ y τ x y ] = [ E x x E x y E x s E x y E y y E y s E x s E y s E s s ] [ ε x ε y γ x y ]
and substituting Equation (10) into Equation (7), using Equation (9), and integrating Equation (7), Up becomes
U p = 1 2 { D x x ( 2 w x 2 ) 2 + 2 ( 2 D s s + D x y ) ( 2 w x 2 ) ( 2 w y 2 ) + 4 D x s ( 3 w x 3 ) ( w y ) + 4 D x s ( 3 w y 3 ) ( w x ) + D y y ( 2 w y 2 ) 2 } d x d y
For a symmetric layup, Dxx, Dxy, Dyy, and Dss are
D x x = E x t 3 12 ( 1 ν x ν y ) ,   D y y = E y t 3 12 ( 1 ν x ν y ) ,   D x y = ν x D y y ,   D s s = G x y t 3 12 ,   D x s = D y s = 0
Generally, Dxs and Dss are non-zero even in a symmetric layup; when the number of layups is large enough, the following assumption can be used.
D x s = D y s = 0
For vertical plates, “y” is changed to “z” in Equations (7)–(13). The total strain energy U and the work W of the four plates is written as
U = 2 × 1 2 0 b 0 l { D x x ( 2 w x 2 ) 2 + 2 ( 2 D s s + D x y ) ( 2 w x 2 ) ( 2 w y 2 ) + 4 D x s ( 3 w x 3 ) ( w y ) + 4 D y s ( 3 w y 3 ) ( w x ) + D y y ( 2 w y 2 ) 2 } d x d y + 2 × 1 2 0 h 0 l { D x x ( 2 v x 2 ) 2 + 2 ( 2 D s s + D x z ) ( 2 v x 2 ) ( 2 v z 2 ) + 4 D x s ( 3 v x 3 ) ( v z ) + 4 D z s ( 3 v z 3 ) ( v x ) + D z z ( 2 v z 2 ) 2 } d x d z W = 2 t 0 b 0 l τ c r w x w y d x d y + 2 t 0 h 0 l τ c r w x w y d x d z
Substituting Equations (3) and (4) into Equation (14) yields the total strain energy U and the work W of the external torsional stress τcr:
U = π 4 4 b 2 h l 3 m = 1 n = 1 Φ m n 2 { b 4 D x x m 4 + 2 b 2 ( 2 D s s + D x y ) m 2 n 2 + D y y n 4 } + π 4 4 b h 2 l 3 m = 1 n = 1 Φ m n 2 { h 4 D x x m 4 + 2 h 2 ( 2 D s s + D x z ) m 2 n 2 + D z z n 4 }
W = 8 t τ c r b h m = 1 n = 1 i = 1 j = 1 m n i j Φ m n Φ i j ( i 2 m 2 ) ( j 2 n 2 ) + 8 t τ c r h b m = 1 n = 1 i j m n i j Φ m n Φ i j ( i 2 m 2 ) ( j 2 n 2 )
where τcr = T/(2bht), i = 1, 2,…; and j = 1, 2,…. The buckling torsional stress is obtained by minimizing the total potential energy.
Φ m n ( U W ) = 0
The following equation is obtained from Equation (17).
π 4 32 τ c r t b h [ ( b l ) 3 { D x x m 4 + 2 ( 2 D s s + D x y ) ( l b ) 2 m 2 n 2 + D y y ( l b ) 4 n 4 } + ( h l ) 3 { D x x m 4 + 2 ( 2 D s s + D x z ) ( l h ) 2 m 2 n 2 + D z z ( l h ) 4 n 4 } ] Φ m n m n l b l h { ( l b ) 2 + ( l h ) 2 } i j i j Φ i j ( i 2 m 2 ) ( j 2 n 2 ) = 0
where summation indicates only those values of i and j for which m + i and n + j are odd. Equation (18) can be rewritten as following matrix form for m + n is even is
[ K 11 τ 0 4 9 β 1 0 0 0 K 13 τ 4 5 β 1 0 0 4 9 β 1 4 5 β 1 K 22 τ 4 5 β 1 20 63 β 1 0 0 4 5 β 1 K 31 τ 0 0 0 20 63 β 1 0 K 15 τ ] [ Φ 11 Φ 13 Φ 22 Φ 31 Φ 15 ] = 0 , K m n = π 4 32 t b h [ ( b l ) 3 { D x x m 4 + 2 ( 2 D s s + D x y ) ( l b ) 2 m 2 n 2 + D y y ( l b ) 4 n 4 } + ( h l ) 3 { D x x m 4 + 2 ( 2 D s s + D x z ) ( l h ) 2 m 2 n 2 + D z z ( l h ) 4 n 4 } ] β 1 = b l h l { ( l b ) 2 + ( l h ) 2 }
and for m + n is odd is
[ K 12 τ 4 9 β 1 0 4 5 β 1 0 4 9 β 1 K 21 τ 8 45 β 1 0 4 5 β 0 8 45 β 1 K 14 τ 8 7 β 1 0 4 5 β 1 0 8 7 β 1 K 23 τ 36 25 β 1 0 4 5 β 1 0 36 25 β 1 K 32 τ ] [ Φ 12 Φ 21 Φ 14 Φ 23 Φ 32 ] = 0 , K m n = π 4 32 t b h [ ( b l ) 3 { D x x m 4 + 2 ( 2 D s s + D x y ) ( l b ) 2 m 2 n 2 + D y y ( l b ) 4 n 4 } + ( h l ) 3 { D x x m 4 + 2 ( 2 D s s + D x z ) ( l h ) 2 m 2 n 2 + D z z ( l h ) 4 n 4 } ] β 1 = b l h l { ( l b ) 2 + ( l h ) 2 }
To obtain a non-zero solution of Φmn, the determinant of the coefficient matrix for Φmn should be zero. The coefficient matrix separates into one for which m + n is odd and one for which m + n is even. It is well known that the thin isotropic plate under shear load, in which m + n is even has a lower buckling load than the one in which m + n is odd except for long plate. Numerical calculations, however, both the case that m + n is even and m + n is odd, were conducted because orthotropic stiffness may lead to lower buckling load for the case that m + n is odd.
The normalized torsional buckling load factor ks is defined as
τ c r = k s π 2 D e t l 2 = k s π 2 t l 2 D x x × D y y + D x y + 2 D s s 2

2.2. Case of a Constant Edge Angle

To solve Equations (19) and (20), the terms of m and n should be limited to a finite number, mmax and nmax. To check the accuracy relative to the number of terms, a numerical calculation was conducted on a quasi-isotropic layup of T300 CFRP. The material properties of the quasi-isotropic layup are listed in Table 1. mmax and nmax = 10, 20,…, 100, and b = h = 100 mm. The results are shown in Figure 4.
The error of mmax = nmax = 20 is within 1% of mmax = nmax = 100 for l/b ≤ 6, and the buckling stress τcr is constant for l/b > 6 because of effect of boundary constraints at x = 0 and x = l vanishing; thus, mmax = nmax = 20 for l/b ≤ 6 and buckling stress τcr (l/b = 6) for l/b > 6 were used in the subsequent numerical calculations.
To check the assumption of a constant edge angle, an FE analysis was conducted using NASTRAN (version2018.1). The dimensions of the pipe were b = h = 100 mm, t = 1 mm, and l = 400 mm. A square mesh with a size of 10 was used. A rigid body element (RBE2) connected the end of the plates to the center of the section and torsional load (torque) was applied to the center of the section. Figure 5 shows the buckling mode and the changes in the angles of the pipe from 90°.

2.3. Case of a Variable Edge Angle

Supposing that the edge angle is variable, substituting Equations (1) and (2) into Equation (14) yields the total strain energy U and the work W of the external torsional stress τcr as follows:
U = π 4 b 2 4 h l 3 m = 1 n = 1 Β m n 2 { D x x m 4 + 2 ( l b ) 2 ( 2 D s s + D x y ) m 2 n 2 + ( l b ) 4 D y y n 4 } + π 4 h 2 4 b l 3 m = 1 n = 1 H m n 2 { h 4 D x x m 4 + 2 ( l h ) 2 ( 2 D s s + D x z ) m 2 n 2 + ( l h ) 4 D z z n 4 }
W = 8 t τ c r ( h b ) 2 m = 1 n = 1 i j m n i j Β m n Β i j ( i 2 m 2 ) ( j 2 n 2 ) + 8 t τ c r ( b h ) 2 m = 1 n = 1 i = 1 j = 1 m n i j H m n H i j ( i 2 m 2 ) ( j 2 n 2 )
The buckling torsional stress is obtained by minimizing the total potential energy.
B m n ( U W ) = 0 , H m n ( U W ) = 0
The following equation is obtained from Equation (17).
π 4 32 τ c r t b 2 { ( b l ) 3 [ D x x m 4 + 2 ( 2 D s s + D x y ) ( l b ) 2 m 2 n 2 + D y y ( l b ) 4 n 4 ] } H m n m n ( b h ) 2 i = 1 j = 1 i j H i j ( i 2 m 2 ) ( j 2 n 2 ) = 0
π 4 32 τ c r t h 2 { ( h l ) 3 [ D x x m 4 + 2 ( 2 D s s + D x z ) ( l h ) 2 m 2 n 2 + D z z ( l h ) 4 n 4 ] } Β m n m n ( h b ) 2 i = 1 j = 1 i j Β i j ( i 2 m 2 ) ( j 2 n 2 ) = 0
The normalized torsional buckling load factor ks was obtained by setting the determinant of the coefficient matrix for Βmn and Hmn to zero. To check the accuracy relative to the number of terms of n and m, a similar numerical calculation as shown in Figure 4 was conducted, and the same conclusion was obtained for the terms of n and m.
The numerical results for Equations (18) and (26), and FEM for h/b = 0.7 are shown in Figure 6. The results for Equations (18) and (26) for h/b = 1.0 are also plotted for comparison. In the case of l/b = 1, the FEM result of h/b = 0.7 is close to the constant-angle case; however, in the region l/b ≥ 2, the FEM result of h/b = 0.7 is close to the variable-angle case. There are no apparent differences between the results of Equation (18) for h/b = 1.0 and those of Equation (26) for h/b = 1.0 because of the section is rectangular and thus the edge angle is constant. These results show that the assumption of a constant edge angle is a stiff constraint for the region h/b < 1.0. Thus, the assumption of a variable edge angle, Equation (26), was employed in the subsequent analysis because of its more accurate and safe-side estimation. Note that both FE analysis and the method of minimizing the total potential energy are essentially the same because they are based on energy principle. If the difference appears, the reason is convergence due to mesh size of FE model or number of terms of displacement function of Equations (1) and (2). The convergence is discussed in Section 2.2 and Figure 4 and Figure 6 shows the present results gives conservative result comparing with FE results, thus the accuracy of the solution is reliable.

3. Closed-Form Polynomial Equation

To make quick and convenient estimates for torsional buckling of a composite tube, a simple closed form polynomial equation was derived by fitting the numerical solution of the orthotropic [(0/90)r,(±45)1-r]s layups. Here, r is the ratio of the 0/90 layer number to the total layer number. Because the [(0/90)r,(±45)1-r]s layup is often used, it is thus chosen. The numerical analysis described above is hard to use for design purposes and would entail a huge calculation cost if it were to be used for optimum design of the whole system. Thus, it is considered that a simple approximate design equation would still be useful. Table 2 shows the material properties of unidirectional pre-preg tape, and Table 3 shows the layup sequence used in the numerical analysis. The total ply number was 32, and the binding stiffness was systematically controlled by varying r.
Assuming that the number of layups was large enough, the bending stiffness was calculated as Eij t3/12 = Dij (i, j = x, y, s). Advanced stiffness properties like through thickness properties discussed in Refs. [19,20,21] are not used because of to make quick and convenient estimates. Figure 7, Figure 8 and Figure 9 show the buckling load factor ks for l/b, h/b, and r. The ranges of l/b, h/b, and r are 1 ≤ l/b ≤ 6, 0.7 ≤ h/b ≤ 1, and 0 ≤ r ≤ 1.
Figure 7 and Figure 8 can be fitted by a second-order polynomial, while Figure 9 can be fitted by a first-order polynomial; thus, the following polynomial was employed.
k s = ( H 0 ( h b ) + H 1 ( h b ) ( l b ) + H 2 ( h b ) ( l b ) 2 ) ( R 1 + R 2 r ) H 0 = H 00 + H 01 ( h b ) ,   H 1 = H 10 + H 11 ( h b ) ,   H 2 = H 20 + H 21 ( h b )
The coefficients H00, H01,…, H21, R1 and R2 were determined by minimizing the following weighted residual and adjusting the weight w to make the non-safe side residual ks+er less than 5%.
k s + e r = k s k s k s   ( f o r   k s k s k s 0   ) , k s e r = | k s k s k s |   ( f o r   k s k s k s < 0   ) max { k s + e r , w k s e r } < 0.05
The minimization was conducted on a Microsoft Excel solver with the GRG Nonlinear option. The following coefficients shown in Table 4 were obtained with w = 0.4; ks+er was 4.77%, and ks-er was 11.9%.
As shown in Figure 4, the buckling stress τcr does not change in the region l/b > 6. Thus, the buckling load factor for the torsional buckling load of the orthotropic rectangular pipe can be described as follows:
τ c r = { k s π 2 D e t l 2   ,   for   l / b 6 k s ( l / b = 6 ) π 2 D e ( 6 b ) 2 t   ,   for   l / b > 6
Figure 10 compares the buckling stress determined by the polynomial equation and the numerical result for mmax = nmax = 100 for the maximum ks-er error case (r = 0.5, b = l = 100 mm, t = 1 mm). The error is also plotted. The error in the region l/b > 6 is around −11%, which gives a safe and conservative side for the design. It is natural that constant gap appears in the region l/b > 6 because ks does not change in the region l/b > 6 and the constant ks’ at l/b = 6 is used. Note that the polynomial equation can be used within the ranges of 1 ≤ l/b ≤ 6, 0.7 ≤ h/b ≤ 1, and 0 ≤ r ≤ 1.

4. Conclusions

Numerical solutions for torsional buckling of an orthotropic rectangular tube were derived. Solutions for different constraint conditions on the edge angle, i.e., with a constant or variable edge angle, were compared with FE analytical results. It was found that a solution under the condition of a variable edge angle is safe and conservative. A simple closed-form polynomial equation of the buckling load factor was also derived for this condition. The equation gives a 4.77% error for the non-safe side and 11.9% error for the safe side in spite of its simple form. In addition, it is well-known that good agreement exists between the theoretical buckling load and experimental buckling load of plates and rectangular pipes, thus the equation can be easily applied to the initial design and optimum design of the whole system.

Author Contributions

Conceptualization, A.T.; methodology, A.T. and R.M.; software, R.M.; validation, A.T.; formal analysis, R.M. and A.T.; writing—original draft preparation, A.T.; writing—review and editing, R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fujita, K.; Nagai, H.; Oyama, A. A Parametric Study of Mars Airplane Concept for Science Mission on Mars. Trans. JSASS Aerosp. Technol. Jpn. 2016, 14, Pk_83–Pk_88. [Google Scholar] [CrossRef]
  2. Nagai, K.; Oyama, A.; Mars Airplane, W.G. Mission Scenario of Mars Exploration by Airplane. In Proceedings of the 2013 Asia-Pacific International Symposium on Aerospace Technology, Takamatsu, Japan, 20–22 November 2013. [Google Scholar]
  3. Airoldi, A.; Bettini, P.; Boiocchi, M.; Sala, G. Composite Elements for Biomimetic Aerospace Structures with Progressive Shape Variation Capabilities. Adv. Technol. Innov. 2016, 1, 13–15. [Google Scholar]
  4. Chang, C.Y. Segmented Compression Molding for Composite Manufacture. Proc. Eng. Technol. Innov. 2017, 5, 41–44. [Google Scholar]
  5. Sahu, V.; Gayathri, V. Strength Studies of Dadri Fly Ash Modified with Lime Sludge–A Composite Material. Int. J. Eng. Technol. Innov. 2014, 4, 161–169. [Google Scholar]
  6. Liu, C.M.; Chiang, M.S.; Chuang, W.C. Lean Transformation for Composite-Material Bonding Processes. Int. J. Eng. Technol. Innov. 2012, 2, 48–62. [Google Scholar]
  7. Takano, A. Buckling Experiment on Anisotropic Long and Short Cylinders. Adv. Technol. Innov. 2016, 1, 25–27. [Google Scholar]
  8. Ning, A.X.; Pellegrino, S. Experiments on Imperfection Insensitive Axially Loaded Cylindrical Shells. Int. J. Solids Struct. 2017, 115–116, 73–86. [Google Scholar] [CrossRef]
  9. Wagner, H.N.R.; Petersen, E.; Khakimova, R.; Hühne, C. Buckling Analysis of an Imperfection-Insensitive Hybrid Composite Cylinder under axial Compression–Numerical Simulation, Destructive and Nondestructive Experimental Testing. Compos. Struct. 2019, 225, 111152. [Google Scholar] [CrossRef]
  10. Podvornyi, A.V.; Semenyuk, N.P.; Trach, V.M. Stability of Inhomogeneous Cylindrical Shells under Distributed External Pressure in a Three-Dimensional Statement. Int. Appl. Mech. 2018, 53, 623–638. [Google Scholar] [CrossRef]
  11. Weaver, P. On Optimisation of Long Anisotropic Flat Plates Subject to Shear Buckling Loads. In Proceedings of the 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference, Palm Springs, CA, USA, 19–22 April 2004. [Google Scholar]
  12. Raju, G.; Wu, Z.; Weaver, P.M. Buckling and Postbuckling of Variable Angle Tow Composite Plates under In-Plane Shear Loading. Int. J. Solids Struct. 2015, 58, 270–287. [Google Scholar] [CrossRef] [Green Version]
  13. Omidvari, A.; Hematiyan, M.R. Approximate Closed-form Formulae for Buckling Analysis of Rectangular Tubes under Torsion. Ije Trans. B: Appl. 2015, 28, 1226–1232. [Google Scholar]
  14. Banks, W.M.; Rhodes, J. The Postbuckling Behaviour of Composite Box Sections. In Composite Structures; Marshall, I.H., Ed.; Springer: Berlin/Heidelberg, Germany, 1981; pp. 402–414. [Google Scholar]
  15. Loughlan, J. The buckling of composite stiffened box sections subjected to compression and bending. Compos. Struct. 1996, 35, 101–116. [Google Scholar] [CrossRef]
  16. Vo, T.P.; Lee, J. Flexural–torsional buckling of thin-walled composite box beams. Thin-Walled Struct. 2007, 45, 790–798. [Google Scholar] [CrossRef] [Green Version]
  17. Column Research Committee of Japan (Ed.) Handbook of Structural Stability; Corona Publishing Company: Tokyo, Japan, 1971; Part 4.
  18. Furusu, K. Studies on the Buckling of the Box Beam Composed of Thin Plates. Ph.D. Thesis, Seikei University, Tokyo, Japan, 2016. (In Japanese). [Google Scholar]
  19. Jones, R.M. Mechanics of Composite Materials; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  20. Ng, Y.C. Deriving composite lamina properties from laminate properties using classical lamination theory and failure criteria. J. Compos. Mater. 2005, 39, 1295–1306. [Google Scholar] [CrossRef]
  21. Kalkan, A.; Mecitoğlu, Z. A Method Based on Classical Lamination Theory to Calculate Stiffness Properties of Closed Composite Sections. J. Aeronaut. Space Technol. 2017, 10, 31–44. [Google Scholar]
Figure 1. Flight image of the Mars exploration airplane.
Figure 1. Flight image of the Mars exploration airplane.
Applsci 11 01342 g001
Figure 2. A cantilevered rectangular pipe under torsion.
Figure 2. A cantilevered rectangular pipe under torsion.
Applsci 11 01342 g002
Figure 3. Assumption of constant plate edge angles.
Figure 3. Assumption of constant plate edge angles.
Applsci 11 01342 g003
Figure 4. Convergence of buckling load for terms of m and n.
Figure 4. Convergence of buckling load for terms of m and n.
Applsci 11 01342 g004
Figure 5. FE analysis for torsional buckling.
Figure 5. FE analysis for torsional buckling.
Applsci 11 01342 g005
Figure 6. FE analysis of torsional buckling.
Figure 6. FE analysis of torsional buckling.
Applsci 11 01342 g006
Figure 7. ks versus l/b.
Figure 7. ks versus l/b.
Applsci 11 01342 g007
Figure 8. ks versus h/b.
Figure 8. ks versus h/b.
Applsci 11 01342 g008
Figure 9. ks versus r.
Figure 9. ks versus r.
Applsci 11 01342 g009
Figure 10. Comparison of polynomial and numerical results.
Figure 10. Comparison of polynomial and numerical results.
Applsci 11 01342 g010
Table 1. Material properties of quasi-isotropic layup of T300 CFRP.
Table 1. Material properties of quasi-isotropic layup of T300 CFRP.
ItemValue
Young’s modulus Ex = Ey [MPa]56,400
Shear modulus Gxy [MPa]21,500
Poisson’s ratio vx = vy [-]0.313
Table 2. Material properties of unidirectional pre-preg tape.
Table 2. Material properties of unidirectional pre-preg tape.
ItemValue
EL [MPa]147,000
ET [MPa]9800
GLT [MPa]5096
vL [-]0.32
Table 3. Layup sequence.
Table 3. Layup sequence.
Number of LayersRatio of 0/90 LayerBending Stiffness [N-mm2]
0/90±45rDxxDxyDyyDss
03203845299638453158
4280.1254187265441872816
8240.254529231345292474
12200.3754870197148702133
16160.55212163052121791
20120.6255553128855531450
2480.755895946.458951108
2840.8756237604.76237766.3
32016578263.16578424.7
Table 4. Coefficients.
Table 4. Coefficients.
H001.1137H010.22213
H10−3.3798H116.0957
H2011.382H21−7.4961
R11.0348R20.26278
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Takano, A.; Mizukami, R.; Kitamura, R. Buckling of Rectangular Composite Pipes under Torsion. Appl. Sci. 2021, 11, 1342. https://doi.org/10.3390/app11031342

AMA Style

Takano A, Mizukami R, Kitamura R. Buckling of Rectangular Composite Pipes under Torsion. Applied Sciences. 2021; 11(3):1342. https://doi.org/10.3390/app11031342

Chicago/Turabian Style

Takano, Atsushi, Ryo Mizukami, and Ryuta Kitamura. 2021. "Buckling of Rectangular Composite Pipes under Torsion" Applied Sciences 11, no. 3: 1342. https://doi.org/10.3390/app11031342

APA Style

Takano, A., Mizukami, R., & Kitamura, R. (2021). Buckling of Rectangular Composite Pipes under Torsion. Applied Sciences, 11(3), 1342. https://doi.org/10.3390/app11031342

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