Next Article in Journal
Global Hypothesis Test to Compare the Predictive Values of Diagnostic Tests Subject to a Case-Control Design
Previous Article in Journal
On the Omega Distribution: Some Properties and Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Buckling of Tapered Heavy Columns with Constant Volume

1
Department of Civil and Environmental Engineering, Wonkwang University, 460 Iksan-daero, Iksan-si, Jeollabuk-do 54538, Korea
2
Department of Civil Engineering, University of Seoul, 163 Seoulsiripdae-ro, Dongdaemun-gu, Seoul 02504, Korea
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(6), 657; https://doi.org/10.3390/math9060657
Submission received: 26 February 2021 / Revised: 12 March 2021 / Accepted: 15 March 2021 / Published: 19 March 2021

Abstract

:
This paper studies the buckling of standing columns under self-weight and tip load. An emphasis is placed on linearly tapered columns with regular polygons cross-section whose volume is constant. Five end conditions for columns are considered. The differential equation governing the buckling shapes of the column is derived based on the equilibrium equations of the buckled column elements. The governing equation is numerically integrated using the direct integration method, and the eigenvalue is obtained using the determinant search method. The accuracy of the method is verified against the existing solutions for particular cases. The effects of side number, taper ratio, self-weight, and end condition on the buckling load and mode shape are investigated. The contribution of self-weight acting alone to the buckling response is also explored. For a given column volume, especially, the buckling length and its stress distribution of the columns with different geometries and end conditions are estimated.

1. Introduction

Columns are elements of the structures in various engineering fields that are subjected to external compressive loads. Long and slender columns have been erected for highways, bridges, offshore facilities, plant structures, etc. In the design of slender mega-columns, self-weight effects are important and must be included in buckling analysis. Such columns are also referred to as heavy columns [1,2]. Tapered members behave differently than uniform members because their variable cross-sections create effective coupling between internal forces and efficient stress distributions [3]. Based on their space utilization, esthetics, safety, optimization, and economic benefits, tapered members are commonly used in engineering practice. Because a tapered member is controlled by its cross-sectional shape and column volume, which affect structural behaviors, various shapes of cross-section are frequently used in practical engineering. Over the past few decade, many efforts have been made to improve structural analyses, including column analysis based on the topics described above.
A short literature review of these topics is provided below. Wang and Drachman [1] investigated the self-weight buckling of a cantilever heavy column with an end load based on a second-order differential equation in terms of the arc length of the buckled column. Interestingly, they applied an inverted cantilever column, which is a column hanging from its fixed end that is subjected to an upward end load. Greenhill [4] studied the maximum stable column lengths (i.e., buckling lengths) of heavy columns such as mast poles. As indicated by the title of the paper, column buckling length was compared to the maximum height at which trees considered as cantilever columns could grow. Since then, small amounts of impactful research have been performed on the buckling analysis of heavy columns: Grishcoff [5] used the infinite series to study the buckling loads of cantilever columns by combining the effects of self-weight and axial loading; Wang and Ang [6] derived buckling load equations for a heavy column subjected to an axial compressive load and restrained by internal supports. Chai and Wang [7] determined the minimum critical buckling load of self-weighted heavy columns under various end conditions using the differential transformation technique. Duan and Wang [2] derived the exact buckling loads of heavy columns under various end conditions in terms of generalized hypergeometric functions. Lee and Lee [8] studied the buckling of a prismatic heavy column under various end conditions, where the buckling length of the column was calculated by considering only its self-weight (without any axial compressive load). Regarding the optimization of heavy columns, tall columns with variable cross-sections and constant volumes were investigated by Keller and Niordson [9], Atanackovic and Glavardanov [10], and Sadiku [11].
For tapered beam/column analysis, various taper functions [3,12,13] along the column axis, including linear, parabolic, sinusoidal, and exponential functions, have been considered. The effects of various cross-sectional shapes [14,15], including rectangular, circular, elliptical, and regular polygons, on the optimization of column buckling have been examined. Additionally, the initial imperfection affecting column behavior was discussed in the open literature [16,17]. The stability of standing heavy column with the intermediate supports, i.e., laterally braced column, was discussed by Wang [18].
Despite the considerable works discussed above, no buckling solutions have been presented in the open literature with a focus on tapered heavy columns and self-weight with regular polygon cross-sections and constant volumes. This study focused on the buckling loads and buckling self-weights of columns under various end conditions. Based on the small deflection theory, a differential equation is derived from the equilibrium equations of the buckled column elements. A direct integral method is developed for integrating the governing equation and the determinant search method is adopted for determining eigenvalues. The predicted results for the buckling load and buckling self-weight are compared to reference values. Numerical results for the buckling load, buckling length, and buckling stress with corresponding mode shapes are presented.

2. Mathematical Formulation

Figure 1a presents an ideal and linear elastic column of span length l placed in a Cartesian coordinate x , y system originating at the toe end t . The toe end t ( x = 0 ) is either hinged or clamped and the head end h ( x = l ) is either free, hinged, or clamped. Therefore, five end condition combinations are possible: “hinged-hinged (H-H)”, “hinged-clamped (H-C)”, “clamped-free (C-F)”, “clamped-hinged (C-H)”, and “clamped-clamped (C-C)”, where the former end represents the toe end and the latter end represents the head end. Columns with H-F end condition were not considered in this study because they are unstable in the structural mechanism from an engineering point of view.
The target columns are linearly tapered with cross-sectional shapes of k   3 -sided regular polygons with circumradii r measured from the centroid to a vertex at any coordinate x . At the toe end t , r is represented as r t . At the head end h , r is represented as r h . The column volume V is always constant. The cross-sectional area and second moment of the plane area at x are denoted as A and I , respectively. In the buckling analysis in this study, self-weight effects were included. Such effects are a major concern in the analysis of heavy columns. The internal self-weight intensity, which is the downward self-weight per unit of axial length induced by column mass and gravity, is represented as F w   = γ A , where γ is the weight density of the column material. The column is subjected to an external compressive load P at the head end and its own self-weight W   = γ V . When P increases and reaches the buckling load B , the column with a buckling length l buckles and forms the buckled-mode shape represented by the solid curve. After column buckling, the internal forces of the axial force N , shear force Q , and bending moment M are applied to the buckled column at x .
To express the taper function of r at x mathematically, the taper ratio n , which is defined as the ratio of the head radius r h to the toe radius r t is introduced.
n = r h r t .
The linear taper function, which is one of the most practical functions in field engineering, of r is expressed in terms of x as follows:
r = r t F 1 ,   F 1 = n 1 x l + 1 ,
where n 1 = n 1 .
By using r in Equation (2), the variable functions of A and I for the k -sided regular polygon at x can be obtained as follows [19]:
A = c 1 r 2 = c 1 r t 2 F 1 2 ,
I = c 2 r 4 = c 2 r t 4 F 1 4 ,
where c 1 and c 2 are
c 1 = k sin π k cos π k ,
c 2 = k 12 sin π k cos 3 π k 3 + tan 2 π k ,
where k   3 is the integer side number and k = for the circular cross-section.
The column volume V is determined as
V = 0 l A d x = c 1 r t 2 0 l F 1 2 d x = c 1 c 3 r t 2 l ,
where c 3 is
c 3 = 1 3 n 2 + n + 1 .
Note that the length l in Equation (7) is the buckling length of the column subjected to an external buckling load B and self-weight W   = γ V .
Based on Equation (7), the circumradius r t can be obtained in terms of V as
r t = V c 1 c 3 l .
By using Equations (3) and (4) with Equation (9), A and I can be obtained in terms of V as
A = V c 3 l F 1 2 ,
I = c 2 V 2 c 1 2 c 3 2 l 2 F 1 4 .
Figure 1b presents a free-body diagram with an infinitesimal length d x for a buckled column element, which is in an equilibrium state based on the internal forces ( N , Q , M ) and self-weight F w = γ A . By setting F x = 0 , F y = 0 , and M = 0 , the three equilibrium equations can be established as
d N d x + F w = 0 ,
d Q d x = 0 ,
d M d x N d y d x Q = 0 .
Differentiating Equation (14) yields the second derivative d 2 M / d x 2 as
d 2 M d x 2 d N d x d y d x N d 2 y d x 2 d Q d x = 0 .
Substituting Equations (12) and (13) into Equation (15) yields
d 2 M d x 2 + F w d y d x N d 2 y d x 2 = 0 .
The self-weight intensity F w at x caused by the γ value of the column material is given by
F w = γ A = γ V c 3 l F 1 2 .
Considering B and F w in Equation (17), the axial force N at x is obtained as
N = B + γ V 0 x F w d x = B + γ V 1 F 2 c 3 ,   F 2 = n 1 2 3 x 3 l 3 + n 1 x 2 l 2 + x l ,
where the term γ V is equal to the total column weight W .
The bending moment M is given by the relationship between load and deformation based on the small deflection theory [19] as
M = E I d 2 y d x 2 .
Differentiating Equation (19) twice yields
d 2 M d x 2 = E d 2 I d x 2 d 2 y d x 2 2 E d I d x d 3 y d x 3 E I d 4 y d x 4 .
Substituting Equations (17), (18), and (20) into Equation (16) yields
d 4 y d x 4 = 2 I d I d x d 3 y d x 3 1 I d 2 I d x 2 d 2 y d x 2 1 E I B + γ V 1 F 2 c 3 d 2 y d x 2 + γ V c 3 l F 1 2 E I d y d x .
From Equation (11), the first and second derivatives of I are determined, respectively:
d I d x = 4 n 1 c 2 V 2 c 1 2 c 3 2 l 3 F 1 3 ,
d 2 I d x 2 = 12 n 1 2 c 2 V 2 c 1 2 c 3 2 l 4 F 1 2 .
Substituting Equation (22a,b) into Equation (21) yields
d 4 y d x 4 = 8 n 1 l F 1 d 3 y d x 3 12 n 1 2 l 2 F 1 2 d 2 y d x 2 c 1 2 c 3 2 c 2 l 2 E V 2 B + γ V 1 F 2 c 3 1 F 1 4 d 2 y d x 2 + c 1 2 c 3 c 2 γ l E V 1 F 1 2 d y d x .
To facilitate numerical analysis and obtain the most general results for this class of problems, the following system parameters are cast into non-dimensional forms:
ξ = x l ,
η = y l ,
β = B l 4 E V 2 ,
λ = γ l 4 E V ,
where ξ , η are non-dimensional Cartesian coordinates, β is the buckling load parameter, and λ is the self-weight parameter.
By using Equations (24)–(27), Equation (23) in dimensional units can be transformed into the non-dimensional differential Equation (28), which governs the buckled shape of the heavy column as
d 4 η d ξ 4 = 8 n 1 f 1 d 3 η d ξ 3 12 n 1 2 f 1 2 d 2 η d ξ 2 c 1 2 c 3 2 c 2 β + λ 1 f 2 c 3 1 f 1 4 d 2 η d ξ 2 + c 1 2 c 3 c 2 λ f 1 2 d η d ξ ,
where f 1 = n 1 ξ + 1 and f 2 = n 1 2 / 3 ξ 3 + n 1 ξ 2 + ξ . The eigenvalues of β , λ in Equation (28) are conjugated with each other. This means that for a given λ value, the eigenvalue β is unique, and vice versa.
Now, consider the boundary conditions in Equation (28). At the top free end ( x = l ), M in Equation (19) and Q in Equation (14) are both equal to zero. Therefore, the non-dimensional boundary conditions of the head free end ξ = 1 are obtained as follows:
d 2 η d ξ 2 = 0 ,   d 3 η d ξ 3 + c 1 2 c 3 2 n 4 c 2 d η d ξ β = 0 .
For the toe and head hinged ends ( x = 0 and x = l ),   y and M are both zero and the non-dimensional boundary conditions at ξ = 0 and ξ = 1 are obtained:
η = 0 ,   d 2 η d ξ 2 = 0 .
For the toe and head clamped ends ( x = 0 and x = l ),   y and d y / d x are both zero and the non-dimensional boundary conditions at ξ = 0 and ξ = 1 are obtained:
η = 0 ,   d η d ξ = 0 .
By using the differential equation in Equation (28) subjected to the selected boundary conditions in Equations (29)–(31), the conjugate eigenvalues of β , λ can be computed using appropriate numerical solution methods for a given set of column parameters for the end conditions ( k and n ).
It is possible for a column to buckle under its self-weight W , even if no external load P = 0 is applied. The buckling self-weight parameter Γ for P = 0 was introduced using Equation (27) and can be formulated as
Γ = γ L 4 E V ,
where L is the self-weight buckling length for which the column buckles under self-weight alone. Setting β = 0 and using Γ instead of λ in Equation (28) yields the following equation:
d 4 η d ξ 4 = 8 n 1 f 1 d 3 η d ξ 3 12 n 1 2 f 1 2 d 2 η d ξ 2 c 1 2 c 3 2 c 2 Γ 1 f 2 c 3 1 f 1 4 d 2 η d ξ 2 + c 1 2 c 3 c 2 Γ f 1 2 d η d ξ ,
where the buckling self-weight parameter Γ is the eigenvalue in the differential equation of Equation (33).
After calculating the conjugate eigenvalues β , λ from Equation (28) for a given set of E , V , and γ , the buckling length l is calculated using Equation (26) or Equation (27), and the buckling stress σ at x is obtained as
l = E V 2 B β 4 ,
σ = N A = c 3 E V l 3 β f 1 2 + 1 f 2 c 3 λ f 1 2 ,
where A and N in Equation (35) are given by Equations (10) and (18), respectively. In particular, the self-weight buckling length L and self-weight buckling stress σ caused only by the self-weight W with P = 0 are obtained using Equations (36) and (37), respectively.
L = E V γ Γ 4 ,
σ = c 3 E V l 3 1 f 2 c 3 Γ f 1 2 .

3. Solution Methods

Based on the mathematical formulations above, two FORTRAN computer programs were written to solve the conjugate eigenvalues of β , λ in Equation (28) and the eigenvalue Γ in Equation (33). The input column parameters are the end conditions, as well as the side number k and taper ratio n for Equations (28) and (33), respectively. To calculate the mode shape ξ , η , Equations (28) and (33), which are boundary problems subjected to the end conditions selected from Equations (29)–(31), are integrated numerically using a direct integration method such as the Runge–Kutta method [20]. The eigenvalues β , λ and Γ are calculated using the determinant search method [15,19] enhanced by the Regula–Falsi method [20]. By using these solution methods, the eigenvalues β , λ and Γ with their corresponding mode shapes ξ , η for various end conditions can be calculated. Five different end conditions are considered, as discussed in Section 2. These types of solution methods for boundary and eigenvalue problems such as Equations (28) and (33) have been described in detail in [15,21], and interested readers should refer to these previous studies.
Before executing the numerical methods described above, it is important to choose a suitable step size Δ ξ when applying the Runge–Kutta scheme, which is computed using the following equation for a given number of dividing elements n e for the unit buckling length:
Δ ξ = 1 n e .
Convergence analysis considering the buckling load parameter β was performed to obtain a suitable n e   = 1 / Δ ξ and the results of the H-H column with a circular cross-section ( k = ) with λ = 1 and n = 0.5 are presented in Figure 2. The solution β with n e = 10 converges at a ratio of 0.9996   ( = 0.26865 / 0.26876 ) to the solution β with n e = 200 , meaning that β with n e = 10 is sufficiently converged. Additionally, the solution β with n e = 20 agrees well with that with n e = 200 (within four significant figures). All computations in this study were conducted on a PC with a GPU. The solutions for β with n e = 20 were computed within one-third of a second.

4. Results and Discussion

Numerical experiments on the effects of column parameters on the conjugate eigenvalues of β , λ in Equation (28) and the eigenvalue Γ in Equation (33) with their corresponding mode shapes ξ , η were performed. For validation purposes, the buckling loads B and buckling self-weight parameters Γ in this study and various references [2,6,8,22] are compared in Table 1 and Table 2, respectively. First, the B values for a concrete column with V = 15 m3, E = 20 GPa, and λ = 0 (i.e., without self-weight, with varying end conditions, a side number k , and taper ratio n ) are compared. The results of this study and those presented by Riley [22] are in good agreement (0.3% error). Second, the buckling self-weight parameters c 1 2 c 3 2 Γ / c 2   ( = γ A L 3 / E I ) for n = 1 (i.e., uniform column) in this study and previous studies [2,6,8] with various end conditions are compared. Note that the parameters of c 1 2 c 3 2 Γ / c 2 for the buckling self-weight parameter Γ have also been adopted in the literature [2,6]. If n = 1 , then the parameters are the same, regardless of k . The results of this study and the references [6] agree well, and the results of this study and the references [2,8] are the same to within five significant figures. Thus, the analytical theories and numerical methods developed in this study are validated when considering all of the column parameters, including the end conditions, k , and n .
Table 3 shows the effects of the side number k on the buckling load parameter β with a conjugate eigenvalue of λ = 1 and n = 0.5 for each end condition. As k increases, β decreases. One can see that an equilateral triangle k = 3 column is the strongest column with the largest β value for a given column volume. This is because the area is the same regardless of k in the same volume, but the circumradius r and the second moment of the plane area I depend on k and are greater from k = 3 to k = (see the ratio of I k / I k = 3 in the last column of the table). The β value of the equilateral triangle column is 1.464   = 0.3934 / 0.2688 times larger than the circular column k = for the H-H condition. The value of β depends heavily on the end conditions, as indicated by C-C column maximum and C-F column minimum. For the circular cross-section, the β value of the C-C column is 19.79   = 2.0759 / 0.1049 times larger than that of the C-F column.
Figure 3 presents β versus n curves for a conjugate eigenvalue of λ = 1 and circular cross-section. Columns subjected to an external load p are in the stability domain under the β versus n curves (i.e., p < b ), meaning they are not buckled. As n increases, β increases, reaches a peak coordinate β , n marked with ▲, and then decreases. At the peak point of β , n of each curve, the taper ratio n is optimized, implying that the column with the optimized n has the maximum β . For example, for the C-H column, the column achieves the maximum β = 1.2814 with an optimized n = 0.8501 . One can see that the β values with n = 1 , excluding the C-F column, are nearly identical to the β values with the optimized n . Figure 3 also highlights the stability region of n . For the C-F column, the columns with 0 < n < 0.7383 are stable (i.e., not buckled), unless p < b . In contrast, the columns with n > 0.7383 are unconditionally unstable (i.e., buckled), even if p = 0 , implying that the columns are buckled by the self-weight parameter λ = 1 . For the C-C and C-H columns, the lower limit of stability for n is n = 0 (see marks of ) and the upper limit of n does not appear until n = 1 . For the H-C and H-H columns, the lower limits of stability for n are n = 0.0949 and n = 0.1426 (see marks of ), respectively, and upper limits of n do not appear until n = 1.
Figure 4 presents a graphical chart of the conjugated eigenvalues of the buckling load parameter and self-weight parameter β , λ for a circular cross-section with n = 0.5 . In the governing differential equation, namely, Equation (28), there are two conjugated eigenvalues of β , λ that are unique. As λ increases, β decreases. β is the largest at λ = 0 when excluding the self-weight effect and the effect of λ on β is significant. For example, β with λ = 2 is 25.5% smaller than β with λ = 0 ( 2.3035 / 1.8353 = 1.255 ; see marks of ●). Eventually, β becomes zero at λ = 8.6443 (i.e., the buckling self-weight parameter Γ = 8.6443 ). Therefore, the column with λ = 8.6443 buckles under the column self-weight alone, without any external load. In this figure, values of Γ marked by are presented for a given set of column parameters.
Figure 5 presents the buckled mode shapes ξ , η for each end condition with a circular cross-section, λ = 1 , and n = 0.5 . In this figure, the buckling load parameters β shown in Table 3 and the positions ξ of the maximum deflection for each end condition are also presented. Note that the coordinate η of the deflection represents relative deflection, rather than absolute deflection. The buckling length parameter λ = 1 for each end condition is the same, but the value of β heavily depends on the end condition. The location of the maximum deflection depends on the end condition. The location of the maximum deflection of a column like a utility pole may be controlled by guywires to prevent unexpected buckling stemming from undesirable column imperfections.
Table 4 shows the effects of the side number k on the buckling self-weight parameter Γ with n = 0.5 for each end condition. As k increases, Γ decreases. An equilateral triangle k = 3 column is the strongest column with the largest Γ value for a given column volume. The Γ value of the triangle column is 1.209   = 2.1405 / 1.7701 times larger than that of the circular column k = for the H-H condition. The value of Γ depends heavily on the end conditions, as indicated by the C-C column maximum and C-F column minimum. For the circular cross-section, the Γ value of the C-C column is 5.257   = 8.6443 / 1.6443 times greater than that of the C-F column. Therefore, selecting proper end conditions is one of the most important design criteria for heavy column design, as discussed previously regarding Table 3.
Figure 6 presents Γ versus n curves for the circular cross-section, where the values of Γ with n = 0.5 listed in Figure 4 are also represented as marks. Columns with the self-weight parameter λ are in the stability domain under the Γ versus n curves (i.e., λ < Γ ) and are not buckled by self-weight. As n increases, Γ increases, reaches a peak at the coordinates Γ , n marked with ▲, and then decreases. At the peak point of Γ , n on each curve, the taper ratio n is optimized to avoid buckling under self-weight, implying that the column with the optimized n has the maximum Γ . For example, for a H-C column, the column achieves the maximum Γ = 2.7164 with an optimized n = 0.5863 .
Figure 7 presents the buckling stresses σ in dimensional units for columns subjected to (a) self-weight without an external load ( P = 0 ) and (b) an external buckling load of B = 5 MN, where the buckling column length l , stress σ t at the column toe ξ = 0 , and stress σ h at the pile head ξ = 1 are presented. The column parameters considered are a circular cross-section, n = 0.5 , V = 10 m3, E = 20 GPa, and γ = 23 kN/m3 for a concrete material. In the case of (a) self-weight, σ decreases along the column axis and σ is maximized as σ m a x = σ t at ξ = 0 , which is the expected behavior. For the buckling column length L , the C-C column is the longest and the C-F column is the shortest, which is the expected behavior. Considering the ultimate stress of σ u = 40 MPa for the concrete material, σ t values between 0.825 and 1.249 MPa are relatively small compared to σ u , meaning heavy column ruptures are caused by buckling, rather than fracturing. In the case of (b), the external load of B = 5 MN, σ increases along the column axis, where σ is minimized as σ m i n = σ t and maximized as σ m a x = σ h because the column is subjected to an external load and the column area decreases (i.e., n = 0.5 ). Additionally, the buckling length l of the C-C column is the largest and that of the C-F column is the smallest. Even when an external load is applied to the column, the column ruptures as a result of buckling, rather than fracturing, just as in the case of self-weight buckling.
Table 5 and Table 6 summarizes the tallest non-buckling column lengths of L , l provided in Figure 7. These tables also include numerical results for a steel heavy column with a square cross-section, E = 210 GPa, and γ = 77 kN/m3, with the other parameters kept constant. The buckling behavior of steel columns is similar to that of concrete columns. It is noteworthy that the self-eight buckling length L (see Equation (36)) and buckling length l (see Equation (34)) of the steel column do not increase significantly beyond those of the concrete column, despite the Young’s moduli of E = 20 GPa for the concrete column and E = 210 GPa for the steel column. Note that under the same column parameters given above, if the length of a particular column is shorter than the tallest length L or l shown in Table 5 and Table 6, the column is safe from column buckling. For example, the H-H column with a specific column length of 10 m ( < L = 62.64   m ) will not be buckling. The corresponding circumradii of the column are r t = 0.752 m and r h = 0.376 m ( n = 0.5 , V = 10   m 3 and column length = 10 m), which are practical in real engineering systems. The column stress σ t = γ V / A t at the toe end is computed as σ t = 0.129 MPa ( < σ a l l = 2   MPa ,   approximately ) , and therefore this column is safe from self-weight buckling.

5. Concluding Remarks

This paper presents the buckling of heavy column included its own self-weight. The column is linearly tapered, the cross-section is a regular polygon, and the column volume is held in constant. Five end conditions of the column are considered. Using equilibrium equations of the buckled column element based on the small deflection beam theory, the fourth-order ordinary differential equation governing the buckled mode shape of such column is derived. For integrating the differential equation, the direct integration method such as Runge–Kutta method is used and for calculating eigenvalues, i.e., buckling load and self-weigh buckling length, is applied as the solution methods. Predicted buckling loads and self-weight buckling lengths agree well with those of references. Numerical results of the buckling load, self-weight buckling length, buckled mode shape, and buckling stress are presented herein and are extensively discussed. The results of this study are expected to be utilized in the design of heavy columns including the self-weight effect.

Author Contributions

Both authors made equal contributions to this paper. Both authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (Grant Number: NRF-2020R1C1C1005374).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, C.Y.; Drachman, B. Stability of a heavy column with an end load. J. Appl. Mech. ASME 1981, 48, 668–669. [Google Scholar] [CrossRef]
  2. Duan, W.H.; Wang, C.M. Exact solution for buckling of columns including self-weight. J. Eng. Mech. ASCE 2008, 134, 116–119. [Google Scholar] [CrossRef]
  3. Lee, J.K.; Lee, B.K. Elastica of non-prismatic and non-linear elastic cantilever beams under combined loadings. Appl. Sci. 2019, 9, 877. [Google Scholar] [CrossRef] [Green Version]
  4. Greenhill, A.G. Determination of the greatest height consistent with stability that a vertical pole or mast can be made, and the greatest height to which a tree of given proportions can grow. J. Phys. Ther. 1882, 1, 337–338. [Google Scholar]
  5. Grishcoff, N. Theory of elastic stability. Bull. Acad. Sci. Russia 1930, 101. [Google Scholar]
  6. Wang, C.M.; Ang, A.K. Buckling capacities of braced heavy columns under an axial load. Comput. Struct. 1988, 28, 563–571. [Google Scholar] [CrossRef]
  7. Chai, Y.H.; Wang, C.M. An application of differential transformation to stability analysis of heavy columns. Int. J. Struct. Stab. Dyn. 2006, 6, 317–332. [Google Scholar] [CrossRef]
  8. Lee, J.K.; Lee, B.K. Buckling lengths of heavy column with various end conditions. Eng. Solid. Mech. 2020, 8, 163–168. [Google Scholar] [CrossRef]
  9. Keller, J.B.; Niordson, F.I. The tallest column. J. Math. Mech. 1966, 165, 433–446. [Google Scholar] [CrossRef]
  10. Atanackovic, T.M.; Glavardanov, V.B. Optimal shape of a heavy compressed column. Struct. Multidiscip. Optim. 2004, 28, 388–396. [Google Scholar] [CrossRef]
  11. Sadiku, S. Buckling load optimization for heavy elastic column: A perturbation approach. Struct. Multidiscip. Optim. 2008, 35, 447–452. [Google Scholar] [CrossRef]
  12. Smith, W.G. Analytical solutions for tapered column buckling. Comput. Struct. 1988, 28, 677–681. [Google Scholar] [CrossRef]
  13. Eroglu, E.; Tufekci, E. Some closed-form solutions for buckling of straight beams with varying cross-section by Variational Iteration Method with Generalized Lagrange Multipliers. Int. J. Eng. Appl. Sci. 2018, 10, 159–175. [Google Scholar] [CrossRef]
  14. Kotodziej, J.A.; Fraska, A. Elastic torsion of bars processing regular polygon in cross-section. Comput. Struct. 2005, 84, 78–91. [Google Scholar] [CrossRef]
  15. Lee, J.K. A unified model for analyzing free vibration and buckling of end-bearing piles. Ocean Eng. 2018, 152, 17–25. [Google Scholar] [CrossRef]
  16. Klasson, A.; Crocetti, R.; Hansson, E.F. Slender steel columns: How they are affected by imperfections and bracing stiffness. Structures 2016, 8, 35–43. [Google Scholar] [CrossRef]
  17. Harvey, P.S.; Cain, T.M.N. Buckling of elastic column with initial imperfection and eccentricity. Structures 2020, 23, 660–664. [Google Scholar] [CrossRef]
  18. Wang, C.Y. Stability of a braced heavy standing column with tip load. Mech. Res. Commun. 2010, 37, 210–213. [Google Scholar] [CrossRef]
  19. Gere, J.M. Mechanics of Materials; Brooks/Cole-Thomson Learning: Belmont, CA, USA, 2004. [Google Scholar]
  20. Burden, R.L.; Faires, D.J.; Burden, A.M. Numerical Analysis; Cengage Learning: Boston, MA, USA, 2016. [Google Scholar]
  21. Lee, J.K.; Park, S.H.; Kim, Y. Transverse free vibration of axially loaded tapered friction piles in heterogeneous soil. Soil Dyn. Earthq. Eng. 2019, 117, 116–121. [Google Scholar] [CrossRef]
  22. Riley, C.E. Elastic Buckling Loads of Slender Column with Variable Cross-Section by Newmark Method. Master’s Thesis, Colorado State University, Fort Collins, CO, USA, 1990. [Google Scholar]
Figure 1. (a) Schematic diagram of a buckled heavy column with a k -sided regular polygon cross-section and (b) forces imposed on a buckled element.
Figure 1. (a) Schematic diagram of a buckled heavy column with a k -sided regular polygon cross-section and (b) forces imposed on a buckled element.
Mathematics 09 00657 g001
Figure 2. Convergence analysis for a suitable n e   = 1 / Δ ξ in the Runge–Kutta scheme.
Figure 2. Convergence analysis for a suitable n e   = 1 / Δ ξ in the Runge–Kutta scheme.
Mathematics 09 00657 g002
Figure 3. Curves of β versus n .
Figure 3. Curves of β versus n .
Mathematics 09 00657 g003
Figure 4. Chart of conjugate eigenvalues of β , λ .
Figure 4. Chart of conjugate eigenvalues of β , λ .
Mathematics 09 00657 g004
Figure 5. Example buckling mode shapes for λ = 1 , k = , and n = 0.5 .
Figure 5. Example buckling mode shapes for λ = 1 , k = , and n = 0.5 .
Mathematics 09 00657 g005
Figure 6. Curves of Γ versus n .
Figure 6. Curves of Γ versus n .
Mathematics 09 00657 g006
Figure 7. Examples of buckling stress σ for a concrete column with a circular cross-section, n = 0.5 , V = 10 m3, E = 20 GPa, and γ = 23 kN/m3: (a) self-weight buckling ( P = 0 ) and (b) external load buckling ( B = 5 MN).
Figure 7. Examples of buckling stress σ for a concrete column with a circular cross-section, n = 0.5 , V = 10 m3, E = 20 GPa, and γ = 23 kN/m3: (a) self-weight buckling ( P = 0 ) and (b) external load buckling ( B = 5 MN).
Mathematics 09 00657 g007
Table 1. Comparisons of buckling load B * for λ = 0 .
Table 1. Comparisons of buckling load B * for λ = 0 .
End   Condition ,   k   and   n Buckling   Load   B in MN
This StudyRiley [22]% Error
H-H, k = 3 , n = 0.4 49.9549.950.
H-C, k = 4 , n = 0.5 109.88109.720.15
C-F, k = 5 , n = 0.6 22.0722.010.27
C-H, k = 6 , n = 0.7 132.39132.330.05
C-C, k = , n = 0.8 270.17270.170.
* l = 15 m, V = 15   m3, and E = 20 GPa for the concrete column.
Table 2. Comparisons of buckling self-weight parameter Γ in terms of c 1 2 c 3 2 Γ / c 2 for n = 1 *.
Table 2. Comparisons of buckling self-weight parameter Γ in terms of c 1 2 c 3 2 Γ / c 2 for n = 1 *.
End Condition c 1 2 c 3 2 Γ / c 2
H-HH-CC-FC-HC-C
This study18.568730.00947.837352.500774.6286
Duan and Wang [2]18.5687-7.837352.500774.6286
Wang and Ang [6]18.58-7.8453.9178.96
Lee and Lee [8]18.568730.00947.837352.500774.6286
* If n = 1 , c 1 2 c 3 2 Γ / c 2 ( = γ A L 3 / E I ) are identical regardless of k .
Table 3. Effects of side number k on β with λ = 1 and n = 0.5 .
Table 3. Effects of side number k on β with λ = 1 and n = 0.5 .
k Buckling   Load   Parameter   β Ratio of
H-HH-CC-FC-HC-C I k / I k = 3
3 (triangle)0.39341.01230.15781.28142.55951.0
4 (square)0.29701.81930.11701.08962.18500.8660
5 (pentagon)0.27890.78330.10921.05382.11520.8410
6 (hexagon)0.27340.77220.10691.04282.09360.8333
(circular)0.26880.76300.10491.03372.07590.8270
Table 4. Effect of side number k on Γ for n = 0.5 .
Table 4. Effect of side number k on Γ for n = 0.5 .
k Buckling   Self - Weight   Parameter   Γ
H-HH-CC-FC-HC-C
3 (triangle)2.14053.24971.98838.014410.453
4 (square)1.85372.81431.72196.94079.0523
5 (pentagon)1.80022.73311.67726.74038.7911
6 (hexagon)1.78372.70801.65696.67878.7106
(circular)1.77012.68741.64436.62788.6443
Table 5. Tallest buckling length L and buckling stresses σ t of heavy columns * without compressive load.
Table 5. Tallest buckling length L and buckling stresses σ t of heavy columns * without compressive load.
L and   σ t H-HH-CC-FC-HC-C
(a) Circular ( k = ) concrete column
L   (m)62.6469.5261.4987.1393.11
σ t   (MPa)0.8400.9330.8251.1691.249
(b) Square ( k = 4 ) steel column
L   (m)84.3293.6082.78117.3125.3
σ t   (MPa)3.7874.2043.7185.2695.630
* Buckling due to self-weight W .
Table 6. Tallest buckling length * l and buckling stresses σ h of heavy columns * with compressive load.
Table 6. Tallest buckling length * l and buckling stresses σ h of heavy columns * with compressive load.
l and   σ h H-HH-CC-FC-HC-C
(a) Circular concrete column
l   (m)21.5325.8717.2625.9430.77
σ h   (MPa)25.1130.1820.1430.2635.90
(b) Square steel column
l   (m)39.4146.9731.6647.3455.94
σ h   (MPa)45.9854.8036.9355.2465.26
* Buckling due to external compressive load P = 5 MN including self-weight W .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, B.K.; Lee, J.K. Buckling of Tapered Heavy Columns with Constant Volume. Mathematics 2021, 9, 657. https://doi.org/10.3390/math9060657

AMA Style

Lee BK, Lee JK. Buckling of Tapered Heavy Columns with Constant Volume. Mathematics. 2021; 9(6):657. https://doi.org/10.3390/math9060657

Chicago/Turabian Style

Lee, Byoung Koo, and Joon Kyu Lee. 2021. "Buckling of Tapered Heavy Columns with Constant Volume" Mathematics 9, no. 6: 657. https://doi.org/10.3390/math9060657

APA Style

Lee, B. K., & Lee, J. K. (2021). Buckling of Tapered Heavy Columns with Constant Volume. Mathematics, 9(6), 657. https://doi.org/10.3390/math9060657

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