Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Unsupervised Anomaly Detection Approach for Time-Series in Multi-Domains Using Deep Reconstruction Error
Next Article in Special Issue
Construction of Analytic Solution to Axisymmetric Flow and Heat Transfer on a Moving Cylinder
Previous Article in Journal
Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model
Previous Article in Special Issue
Optimal Perturbation Iteration Method for Solving Fractional Model of Damped Burgers’ Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Buckling of Corrugated Ring under Uniform External Pressure

by
Igor I. Andrianov
1,*,
Igor V. Andrianov
2,
Alexander A. Diskovsky
3 and
Eduard V. Ryzhkov
3
1
Rhein Energie AG, Parkgürtel 24, D-50823 Köln, Germany
2
Institut für Allgemeine Mechanik, RWTH Aachen University, Templergraben 64, D-52056 Aachen, Germany
3
Department of Economic and Information Security, Dnipro State University of Foreign Affairs, Gagarina Ave 26, UA-49005 Dnipro, Ukraine
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(8), 1250; https://doi.org/10.3390/sym12081250
Submission received: 16 May 2020 / Revised: 22 June 2020 / Accepted: 1 July 2020 / Published: 29 July 2020
(This article belongs to the Special Issue Asymptotic Methods in the Mechanics and Nonlinear Dynamics)

Abstract

:
Stability analysis of a corrugated ring subjected to uniform external pressure is under consideration. Two main approaches to solving this problem are analyzed. The equivalent bending stiffness approach is often used in engineering practice. It is based on some plausible assumptions about the behavior of a structure. Its advantage is the simplicity of the obtained relations; the disadvantage is the difficulty in estimating the area of applicability. In this paper, we developed an asymptotic homogenization method for calculating the critical pressure for a corrugated ring, which made it possible to mathematically substantiate and refine the equivalent bending stiffness approach. To evaluate the results obtained using the equivalent stiffness approach and asymptotic homogenization method, the imperfection method is used. The influence of the corrugation parameters on buckling pressure is analyzed.

1. Introduction

Corrugation of smooth shells is widely used to change their stiffness. In particular, for cylindrical shells under external pressure, the corrugation of the generatrix is used [1,2,3,4,5,6,7]; corrugation of the directrix significantly increases the stability of such shells under axial compression [8]. Therefore, the study of the stability of corrugated shells under external pressure is of practical importance.
The methods currently used to calculate corrugated plates and shells can be subdivided into three groups: computational models based on finite element analysis (FEA), the equivalent bending stiffness approach, based on some physical assumptions, and the asymptotic homogenization method.
FEA in principle makes it possible to calculate any corrugated shell [6,7,8,9,10,11,12,13,14,15], but, for the early stages of design, simple analytic solutions are very useful to engineers. Analytical solutions can be also used as benchmark examples for numerical algorithms.
Before the computational revolution, the main method for calculating corrugated shells was equivalent stiffnesses modeling [16,17]. A common approach is to replace a corrugated shell with an anisotropic one that has equivalent stiffness properties. This method has been used in the analysis of corrugated shells since the 1930s [17]. The main disadvantage of this approach is the difficulty in estimating accuracy and the area of applicability [18,19,20,21,22,23,24,25,26].
In other papers [27,28,29,30,31], the asymptotic homogenization method was used. The original equations of the corrugated shell are projected on the basic surface, equally spaced from the corrugation crests. This allows the establishment of a relation between the homogenized and real components of the stress–strain state. In the zeroth approximation, the slowly changing components of the solution are determined. The subsequent approximations are rapidly oscillating periodic functions (correctors) and they are found in the solution to the cell problem. This approach allows the solution of a number of optimization problems and the consideration of structures with functionally graded corrugations [29,30,31].
The buckling problem for longitudinally corrugated cylindrical shells under external pressure was investigated by Semenyuk and Babich [9]. They analyzed simply supported cylindrical shells, of which cross sections are described in polar coordinates by Fourier series R ( φ ) = R 0 + Σ k = 1 R k cos ( k φ ) . The solution was represented as the product of the sine in the longitudinal direction on the Fourier series in the circumferential direction. Numerical results show the decrease in critical pressure for a corrugated shell in comparison with a cylindrical shell of the same length and radius R0.
The nonlinear elastic deformation of a thin flexible ring under the external hydrostatic pressure is described in [15].
In our work, we study the stability of the corrugated ring subjected to external pressure. Obtained results can be used for long cylindrical shells, when the number of waves in the circumferential direction is two.
The paper is organized as follows. First, we employ the basic relations in Section 2. In Section 3, the equivalent bending stiffness approach is carried out. Section 4 deals with homogenization of the basic relations. Section 5 presents an imperfection method. Finally, Section 6 presents concluding remarks.

2. A Statement of the Problem

The radius of the corrugated ring in the polar coordinate system is shown in Figure 1.
r = R ( 1 + h g ( n φ ) )
where R is the radius of the basic circular ring, equally spaced from the crests of the corrugation, h = H R ,   H is the corrugation depth (Figure 2), g ( n φ ) is the periodic function with period 2 π / n (pitch of corrugation), n is the number of corrugations, and 0 φ 2 π , | h g | 1 .
In this case, arc ( φ 1 , φ 2 ) length s and radius of curvature ρ of the ring are:
s = φ 1 φ 2 A d φ ;   ρ = A 3 r 2 + 2 r φ 2 r r φ φ
where A = r 2 + r φ 2 , ( · ) φ = d ( · ) d φ .
We suppose
h 1 ; 2 π n 1
Further, we choose the corrugation profile most convenient for analysis (Figure 1):
r = R ( 1 + h c o s ( n φ ) )
We assume that, for the ring, the hypothesis of flat cross-section is satisfied and the magnitude of the bending moment is proportional to the change in the curvature of the middle surface of the ring, in the section under consideration.
We deal with the general instability of the corrugated ring, i.e., the transition to the form with non-circular basic ring. Local buckling modes are absent. The following assumptions are also valid: spatial forms of equilibrium are not considered; during deformation the load remains directed normal to the deformed ring, and its intensity does not change.
The problem is solved in a linear formulation, using the static criterion of bifurcational stability loss. As shown in a number of studies [32,33,34,35], if the initial problem is conservative and the load is applied quasistatically, a similar approach to studying the stability of elastic systems is justified.

3. Equivalent Bending Stiffness Approach

Let us replace the corrugated ring with a circular one, radius R, with the equivalent bending stiffness. To determine this stiffness, following [16], we single out one corrugation wave, and due to the large number of such waves, the curvature of the basic circle can be neglected (Figure 2). Effective bending stiffness, D, is stiffness that ensures equality of deformations of the original curvilinear and approximating rectilinear beams. As a result, one obtains
D = E I k
where k = 1 l s = A 1 ,   = 1 2 π R 0 2 π ( ) d φ ; E is the Young’s modulus of the ring material; I is the moment of inertia of the ring cross-section; l s = 0 2 π / n A d φ   and   l = 2 π R / n are the lengths of the wave of the corrugation and its projection.
To evaluate the intensity of external pressure p 0 acting on an equivalent basic ring, we select arcs d s and d l (Figure 3). Projecting the load shows that the normal pressure on the basic ring is the same as the original one.
Note that the projection also gives a tangential load in the basic circular ring, self-balancing along one pitch of corrugation. Additional tangential load is not considered in the equivalent bending stiffness approach. Above, we deal with pressure which is permanently normal to the arc of the ring. If the external pressure is caused by radial forces directed permanently to the center (see [32], problem 134), then intensity p must be replaced to p k .
Buckling pressure for the circular ring with bending stiffness (5) is [34]
p b = k
where p b = p b 1 p b 0 ,     p b 1 = 3 k E I R 3 is the buckling pressure of the circular ring with bending stiffness (5); p b 0 = 3 E I R 3   is   the buckling pressure of the circular ring with radius, R, and bending stiffness E I .
Corrugation profile (4) has two important characteristics, namely, n and h . The buckling pressure (6) depends on the parameter n h (Table 1), which can be called the corrugation parameter. For h = 0 , one has k = 1 ; by increasing n h , the critical pressure decreases (Figure 4).
Comparing of the corrugated and circular rings with the same perimeter shows that corrugation increases the critical pressure
p b p = k 2
where p b p = P b P b p ,   P b p = 3 E I ( k 1 R ) 3 is the buckling pressure of the circular ring of the same perimeter.
Radius R 1 of the circular ring, with the perimeter equal to the perimeter of the corrugated one, can be calculated as follows R 1 = 0 2 π A d φ / 2 π .
Using Formulas (6) and (7), one easily determines, for corrugation profile (4), that the buckling load of the corrugated ring is 31% higher than for the circular ring of the same perimeter. On the other hand, the corrugated ring restricts 23.6% less area. Comparison using (6) with the basic circular ring (Figure 1) shows a decrease in the critical pressure for the corrugated one by 12.6% and a decrease in the restricted area by 0.125%. In this case, the perimeter of the corrugated ring is 14.5% larger than the perimeter of the basic ring.
The questions remain unclear:
(a)
the applicability of the equivalent bending stiffness approach to the stability problems;
(b)
estimation of the error caused by neglecting the curvature of the basic circular ring in determining equivalent bending stiffness;
(c)
estimation of the accuracy of the buckling pressure obtained with the equivalent bending stiffness approach and the possibility of its refinement.

4. Asymptotic Homogenization Method

Suppose that under the influence of an external pressure, the initial equilibrium state of the corrugated ring transits to a new equilibrium state, which has two perpendicular axes of symmetry. We write the equation describing this state and find the minimum value of pressure at which such a state becomes stable.
The symmetry axes of the new equilibrium form are taken as the x, y coordinate axes and consider the part of the ring lying above the horizontal axis of symmetry (Figure 5).
The bending moment M in the cross-section B is
M = M ˜ + N ˜ | A D ¯ | p | A B ¯ | 2 2
where N ˜ = p ( r w ˜ ) .
Expression (8) can be written as
M = M ˜ p ( | A B ¯ | 2 2 | A O ¯ | | A D ¯ | )
From triangle ABO (Figure 5), one obtains
| A B ¯ | 2 2 | A O ¯ | | A D ¯ | = 1 2 ( | O B ¯ | 2 | O A ¯ | 2 ) = 1 2 ( ( r w ) 2 ( r ˜ w ˜ ) 2 )
Substituting expression (10) into Equation (9), after linearization with respect to w and w ˜ , we obtain
M = M ˜ 1 2 p ( r 2 r ˜ 2 + 2 r ˜ w ˜ 2 r w )
We express the bending moment (11) through the projections of displacements on the basic circular ring (Figure 6)
M = E I 1 A ( 1 A ( r u + r φ w ρ ( r w r φ u A ) φ ) ) φ
Expression (12) is reduced to the form
M = E I 1 A ( 1 A 2 ( r φ ε r ϑ ) ) φ
where
ε = u φ + w ;   ϑ = w φ u ,
ε ,   ϑ are the circumferential deformation and the angle of rotation of the tangent to the basic circular ring [34].
We suppose inextensibility conditions for basic ring
u φ + w = 0
Taking into account condition (15), expression (13) can be rewritten as follows
M = E I 1 A ( r ( w φ u ) A 2 ) φ .
Substituting the expression for moment (16) into Equation (11), one obtains
( r ( w φ u ) A 2 ) φ + p ¯ A r w = A 2 ( p ¯ ( r 2 r ˜ 2 + 2 r ˜ w ˜ ) 2 M ¯ )
where p ¯ = p E I ;   M ¯ = M ˜ E I .
For eigenvalue problem (15), (17) must be supplemented with symmetry conditions (Figure 5)
w φ ( 0 ) = w φ ( π 2 ) = 0
So, we obtained an ODE with rapidly oscillating coefficients. This kind of problem can be efficiently solved by the homogenization method [36,37,38] based on the two-scale asymptotic expansions.
Let us introduce a new variable, ξ = n φ , which is assumed to be the independent on φ , hence
d d φ = φ + n ξ
The projections of the displacements u ,   w are represented as series
u = k = 0 n k u k ( φ , ξ ) ;   w = k = 0 n k w k ( φ , ξ ) ,   corr
where u k ,   w k are the ξ-periodic functions with period 2 π .
Substituting expansions (20) into Equations (15) and (17), after splitting in powers of n k , we obtain
ξ ( a ( ξ ) w 0 ξ ) = 0
u 0 ξ = 0
ξ ( a ( ξ ) w 1 ξ ) + ξ ( a ( ξ ) ( w 0 φ u 0 ) ) = 0
ξ ( a ( ξ ) w 2 ξ ) + ξ ( a ( ξ ) ( w 1 φ u 1 ) ) + a ( ξ ) 2 w 1 φ ξ + a ( ξ ) φ ( w 0 φ u 0 ) + p ¯ A r w 0 = A 2 ( p ¯ ( r 2 r ˜ 2 + 2 r ˜ w ˜ ) 2 M ¯ )
u 0 φ = w 0
where a ( ξ ) = r A 2 .
Using Equations (21) and (22), one obtains
w 0 = w 0 ( φ ) ; u 0 = u 0 ( φ )
Then Equation (23) can be rewritten as follows
ξ ( a ( ξ ) w 1 ξ ) = a ( ξ ) ξ ( w 0 φ u 0 )
Integrating Equation (27) over ξ , we obtain
w 1 ξ = w 0 φ + u 0 + C ( φ ) a ( ξ )   .
The integration “constant” C ( φ ) is determined from the periodicity condition for the function w 1 with respect to ξ :
C ( φ ) = a ^ ( w 0 φ u 0 ) , a ^ = ( 1 2 π 0 2 π a 1 d ξ ) 1
Excluding derivative w 1 ξ from Equation (24), one obtains
ξ ( a ( ξ ) w 2 ξ ) + ξ ( a ( ξ ) ( w 1 φ u 1 ) ) + a ^ φ ( w 0 φ u 0 ) + p ¯ A r w 0 = A 2 ( p ¯ ( r 2 r ˜ 2 + 2 r ˜ w ˜ ) 2 M ¯ )
We apply to Equation (30) the averaging operator, 0 2 π ( ) d ξ . The first two terms vanish due to the periodicity in ξ , and finally we have
a ^ d d φ ( d w 0 d φ u 0 ) + p ¯ s ˙ w 0 = p ¯ 2 ( s ¨ s r ˜ ( r ˜ 2 w ˜ ) ) s M ¯
where s ˙ = 1 2 π 0 2 π A r d ξ   ; s ¨ = 1 2 π 0 2 π A r 2 d ξ   ;   s = 1 2 π 0 2 π A d ξ   .
Using expression (25), Equation (31) is reduced to
d 2 w 0 d φ 2 + β 2 w 0 = α
where β 2 = 1 + p ¯ k ^ ;   k ^ = a ^ ˙ s ˙ ;   α = p ¯ 2 a ^ ( s ¨ s r ˜   ( r ˜ 2 w ˜ ) ) s M ¯ a ^ .
The general solution of the ODE (32) is
w 0 = C 1 sin β φ + C 2 cos β φ + α β 2
where C 1 ,   C 2 are the integration constants.
For function w 0 , one has
d w 0 d φ = 0   at   φ = 0 ,   π 2
Substituting expression (33) into Equation (34), we find
C 1 = 0 ; β = 2 m , m = 1 , 2 ,
Assuming m = 1 , one finds from expression (35) the buckling pressure of the corrugated ring
p ¯ b 2 = 3 k ^
Thus, the rationale for the use of the equivalent bending stiffness approach to the problem of the stability of a corrugated ring is obtained.
Figure 7 shows how the value of the critical pressure of the corrugated ring p ¯ b 2 decreases compared to a value for the circular ring of radius R with an increase in the depth of the corrugation h .
Thus, an increase in the depth of the corrugation decreases the critical pressure for the corrugated ring.
Figure 8 shows a comparison of buckling pressures p ¯ b 2   and   p ¯ b 2 e for corrugated and circular rings having the same volume of material. It is obvious that it is more optimal to increase the buckling pressure by increasing the cross section of a circular ring than by corrugating it.
Asymptotic homogenization method gives less of a value of a buckling load than equivalent bending stiffness approach. This difference increases with the increase of the depth of the corrugation h (Figure 9).
The question of estimating the accuracy of the results obtained using the equivalent bending stiffness approach and asymptotic homogenization method remains open.

5. Using the Imperfection Method

“The imperfection method is characterized by the question: What is the value of the load for which the static displacements of the imperfect system become excessive or even infinite?” [33]. This approach was widely used by Timoshenko [34]. Critical re-examination of the imperfection method was described by Gubanova and Panovko [35].
Let us write the ODE of the deformed midline of the ring in the framework of the flat cross-section hypothesis
1 ρ 1 ρ ¯ = M E I
where 1 ρ ¯ = ( r w ) 2 + 2 ( r w ) φ 2 ( r w ) ( r w ) φ φ A ¯ 3 ;   A ¯ = ( r w ) 2 + ( r w ) φ 2 .
Substituting the expansion A ¯ in powers of w R and expression (11) for the bending moment M to the equilibrium Equation (37), after ignoring the second-order terms, one obtains
r w φ φ + 3 A ( r r φ + 2 r φ 3 r 4 3 r φ A r φ r φ φ ) w φ + 1 A ( A r φ φ 2 r A + 3 r 2 + 6 r φ 2 3 r r φ φ + p r A 4 ) w = ( p ¯ ( r 2 r ˜ 2 + 2 r ˜   w ˜ ) 2 M ¯ ) A 3 2
The solutions of Equation (38) must satisfy the symmetry conditions (18).
Numerical solutions of the boundary value problem (38), (18) were obtained on the basis of the following iterative scheme: at the first iteration, we suppose w ˜ = 0 ;   M ¯ = 0 and calculate w ˜ = w ( 0 ) ;   M ¯ = 1 ρ ( 0 ) 1 ρ ¯ ( 0 ) . Then, new values of w ˜ and M ¯ are substituted into Equation (39) and used in the second iteration. Subsequent iterations are carried out similarly. This algorithm was used for the following corrugation profile
r = 1 + h c o s ( 2 π n φ )
Symmetry conditions are
w φ ( 0 ) = w φ ( 1 4 ) = 0
Comparison of the results of calculation of the deflection at the point φ = 0 for n = 16 , h = 0.01 , showed a satisfactory convergence of the used iterative scheme (Table 2).
To estimate the critical pressure of the corrugated ring, the dependences of the deflection at the point were φ = 0 on the external pressure for corrugation profile (39) with parameters n = 16, h = 0.01 (Figure 10a); n = 8, h = 0.02 (Figure 10b).
The pressure at which a sharp increase in the deflection occurs gives estimation of buckling load. The same figures show that the buckling pressure values found using the equivalent bending stiffness approach p ¯ b 1 , asymptotic homogenization method p ¯ b 2 , and buckling pressure for circular ring p ¯ b 0 = 3 .

6. Conclusions

One can increase the critical pressure of a circular ring of a given perimeter using corrugation. However, a corrugated ring restricts less area, which is of practical importance in the design of cylindrical tanks. Comparison with the basic circular ring shows a decrease in critical pressure. Thus, long cylindrical circumferentially-corrugated shells have less rigidity at external pressure than circular shells with basic radius.
It is well known [39,40,41] that initial imperfections lead to the decrease of critical pressure for shells. However, in this case, one deals with initial imperfections whose change in circumferential direction coincides with a buckling form of the ideal shell. A circular ring buckles along two waves; therefore, many-waves corrugation increases its stability.
We solved our problem under the assumption that the pre-buckling state is momentless. For circular cylindrical shells under axisymmetric transverse pressure, this assumption is correct [42]. In our case, the correctors to the homogenized solution are rapidly variable, therefore, their projection onto the slow changing buckling form is insignificant.
The critical pressure of the corrugated ring depends on the product of two parameters of the corrugation: number of corrugations n and their depth h . Rings with the same values of parameter nh, have the same critical pressure.
The asymptotic homogenization method gave the possibility of mathematically justifying the equivalent bending stiffness approach and corrected the values of reduced bending rigidity and buckling pressure.
In our paper, corrugation profile has a cosine-shaped configuration. In engineering practice, a lot of other configuration, such as trapezoidal, rectangular, and zigzag, profiles are used. The main conclusions that we obtained based on the analysis of the simplest case remain unchanged. Namely: the corrugation reduces the critical pressure, the decrease depends mainly on the corrugation parameter nh, with increasing nh, and the critical pressure decreases.
Asymptotic homogenization method is applicable for the above profiles according to the scheme developed in this paper. Only the solutions to the cell problem will change, and therefore the values of equivalent stiffnesses will be different.
Replacing the original corrugated ring with a circular one of the same perimeter leads to a symmetrization of the problem. It was shown in [43] that an increase in symmetry in eigenvalue problems for ODEs with periodic coefficients leads to an increase in eigenvalues. Therefore, the found-above values of buckling loads are overestimated.
A few words concerning possible generalizations: experimental studies [7] indicate the need to take into account, in many cases, inelastic general instability and local buckling modes. When deriving nonlinear equations by taking into account inelastic deformation, the results obtained in [44,45] can be used. The models proposed in [25,26] could be used for improving the original model, accounting more exactly for the coupling between tangential and bending loads.

Author Contributions

Software, I.I.A. and E.V.R.; methodology, I.V.A.; formal analysis, A.A.D. All authors discussed all stages of the research and preparing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The research and publication of this article were funded by the authors themselves. The authors thank the anonymous referees whose valuable comments and suggestions favored improvement of the paper.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Data Availability

All the data required for the publication have been included in the research article. These are in the form of tables and plots in the research paper.

Abbreviations

The following abbreviations are used in this manuscript:
ODEOrdinary Differential Equation
FEAfinite element analysis

References

  1. Dayyani, I.; Shaw, A.D.; Flores, E.I.S.; Friswell, M.I. The mechanics of composite corrugated structures: A review with applications in morphing aircraft. Comp. Struct. 2015, 133, 358–380. [Google Scholar] [CrossRef] [Green Version]
  2. Adriaenssens, S.; Dejong, M. Form finding of corrugated shell structures for seismic design and validation using non-linear pushover analysis. Eng. Struct. 2019, 181, 362–373. [Google Scholar]
  3. Malek, S.; Williams, C.J.K. The equilibrium of corrugated plates and shells. Nexus Netw. J. 2017, 19, 619–627. [Google Scholar] [CrossRef] [Green Version]
  4. Manko, Z.; Beben, D. Research on steel shell of a road bridge made of corrugated plates during backfilling. J. Bridge Eng. 2005, 10, 592–603. [Google Scholar] [CrossRef]
  5. Wang, K.; Zhou, M.; Hassanein, M.F.; Zhong, J.; Ding, H.; An, L. Study on elastic global shear buckling of curved girders with corrugated steel webs: Theoretical analysis and FE modelling. Appl. Sci. 2018, 8, 2457. [Google Scholar] [CrossRef] [Green Version]
  6. Semenyuk, N.P. Stability of corrugated arches under external pressure. Int. Appl. Mech. 2013, 49, 211–219. [Google Scholar] [CrossRef]
  7. Ross, C.; Humphries, M. The buckling of corrugated circular cylinders under uniform external pressure. Thin Walled Struct. 1993, 17, 259–271. [Google Scholar] [CrossRef]
  8. Semenyuk, N.P.; Zhukova, N.B.; Ostapchuk, V.V. Stability of corrugated composite noncircular cylindrical shells under external pressure. Int. Appl. Mech. 2007, 43, 1380–1389. [Google Scholar] [CrossRef]
  9. Semenyuk, N.P.; Babich, I.Y. Stability of longitudinally corrugated cylindrical shells under uniform surface pressure. Int. Appl. Mech. 2007, 43, 1236–1247. [Google Scholar] [CrossRef]
  10. Semenyuk, N.P.; Zhukova, N.B.; Neskhodovskaya, N.A. Stability of orthotropic corrugated cylindrical shells under axial compression. Mech. Comp. Mater. 2002, 38, 243–250. [Google Scholar] [CrossRef]
  11. Norman, A.D.; Seffen, K.A.; Guest, S.D. Multistable corrugated shells. Proc. R. Soc. A 2008, 464, 1653–1672. [Google Scholar] [CrossRef]
  12. Wang, X.J.; Wang, Z.M.; Wang, N. Finite element analysis for the stiffness and the buckling of corrugated tubes in heat exchanger. Adv. Mater. Res. 2012, 468–471, 1675–1680. [Google Scholar] [CrossRef]
  13. Wennberg, D.; Wennhage, P.; Stichel, S. Orthotropic models of corrugated sheets in finite element analysis. Int. Sci. Res. Netw. Mech. Eng. 2011, 2011, 979532. [Google Scholar] [CrossRef] [Green Version]
  14. Sowiński, K. Buckling of shells with special shapes with corrugated middle surfaces–FEM study. Eng. Struct. 2019, 179, 310–320. [Google Scholar] [CrossRef]
  15. Dolgikh, D.V.; Kiselev, V.V. Corrugation of a flexible ring under external hydrostatic compression. J. Appl. Math. Mech. 2010, 74, 204–213. [Google Scholar] [CrossRef]
  16. Andreeva, L.E. Elastic Elements of Instruments; Israel Program for Scientific Translations: Jerusalem, Israel, 1966. [Google Scholar]
  17. Seydel, E.B. Schubknickversuche Mit Wellblechtafeln; Jahrbuch d. Deutsch; Versuchsanstallt für Luftfahrt. E.V.: Munich/Berlin, Germany, 1931; Volume 4, pp. 233–235. [Google Scholar]
  18. Briassoulis, D. Equivalent orthotropic properties of corrugated sheets. Comput. Struct. 1986, 23, 129–138. [Google Scholar] [CrossRef]
  19. Donnell, L.H. The flexibility of corrugated pipes under longitudinal forces and bending. Trans. ASME 1932, 54, 69–75. [Google Scholar]
  20. Ye, Z.; Berdichevsky, V.L.; Yu, W. An equivalent classical plate model of corrugated structures. Int. J. Solids Struct. 2014, 51, 2073–2083. [Google Scholar] [CrossRef] [Green Version]
  21. Xia, Y.; Friswell, M.I.; Flores, E.I.S. Equivalent models of corrugated panels. Int. J. Solids Struct. 2012, 49, 1453–1462. [Google Scholar] [CrossRef] [Green Version]
  22. Nguyen-Minh, N.; Tran-Van, N.; Bui-Xuan, T.; Nguyen-Thoi, T. Static analysis of corrugated panels using homogenization models and a cell-based smoothed Mindlin plate element. Front. Struct. Civ. Eng. 2019, 13, 251–272. [Google Scholar] [CrossRef]
  23. Kolpakov, A.G.; Rakin, S.I. Calculation of the effective stiffness of the corrugated plate by solving the problem on the plate cross-section. J. Appl. Mech. Tech. Phys. 2016, 57, 757–767. [Google Scholar] [CrossRef]
  24. Kolpakov, A.A.; Kolpakov, A.G. Discussion of the effective stiffnesses in: Ye, Berdichevsky, and Yu [Int. J. Solids Struct. 51 (2014) 2073–2083]. Int. J. Solids Struct. 2019, 174–175, 145–146. [Google Scholar] [CrossRef]
  25. Gimena, L.; Gimena, F.; Gonzaga, P. Structural analysis of a curved beam element defined in global coordinates. Eng. Struct. 2008, 30, 3355–3364. [Google Scholar] [CrossRef]
  26. Mercuri, V.; Balduzzi, G.; Asprone, D.; Auricchio, F. Structural analysis of non-prismatic beams: Critical issues, accurate stress recovery, and analytical definition of the Finite Element (FE) stiffness matrix. Eng. Struct. 2020, 213, 110252. [Google Scholar] [CrossRef]
  27. Andrianov, I.V.; Diskovsky, A.A.; Kholod, E.G. Homogenization method in the theory of corrugated plates. Tech. Mech. 1998, 18, 123–133. [Google Scholar]
  28. Syerko, E.; Diskovsky, A.A.; Andrianov, I.V.; Comas-Cardona, S.; Binetruy, C. Corrugated beams mechanical behavior modeling by the homo–genization method. Int. J. Solids Struct. 2013, 50, 928–936. [Google Scholar] [CrossRef] [Green Version]
  29. Andrianov, I.V.; Awrejcewicz, J.; Diskovsky, A.A. Design optimization of FGM beam in stability problem. Eng. Comput. 2019, 36, 248–270. [Google Scholar] [CrossRef]
  30. Andrianov, I.I.; Awrejcewicz, J.; Diskovsky, A.A. The optimal design of a functionally graded corrugated cylindrical shell under axisymmetric loading. Int. J. Nonlinear Sci. Numer. Simul. 2019, 20, 387–398. [Google Scholar] [CrossRef]
  31. Andrianov, I.V.; Diskovsky, A.A.; Syerko, E. Optimal design of a circular diaphragm using the homogenization approach. Math. Mech. Solids 2017, 22, 283–303. [Google Scholar] [CrossRef]
  32. Feodosiev, V.I. Advanced Stress and Stability Analysis, Worked Examples; Springer: Berlin, Germany, 2005. [Google Scholar]
  33. Bensoussan, A.; Lions, J.-L.; Papanicolaou, G. Asymptotic Analysis for Periodic Structures; North-Holland: Amsterdam, The Netherlands, 1978. [Google Scholar]
  34. Andrianov, I.V.; Awrejcewicz, J.; Manevitch, L.I. Asymptotical Mechanics of Thin-Walled Structures: A Handbook; Springer: Berlin, Germany, 2004. [Google Scholar]
  35. Manevich, L.I.; Andrianov, I.V.; Oshmyan, V.G. Mechanics of Periodically Heterogeneous Structures; Springer: Berlin, Germany, 2002. [Google Scholar]
  36. Ziegler, H. Principles of Structural Stability, 2nd ed.; Springer: Basel, Switzerland, 1977. [Google Scholar]
  37. Timoshenko, S.P.; Gere, J.M. Theory of Elastic Stability, 2nd ed.; McGraw-Hill: New York, NY, USA, 1961. [Google Scholar]
  38. Panovko, Y.a.G.; Gubanova, I.I. Stability and Oscillations of Elastic Systems Paradoxes, Fallacies, and New Concepts; Consultants Bureau: New York, NY, USA, 1965. [Google Scholar]
  39. Grigolyuk, E.I.; Kabanov, V.V. Stability of Shells; Nauka: Moscow, Russia, 1978. (In Russian) [Google Scholar]
  40. Vol’mir, A.S. Stability of Deformable Systems; Foreign Technology Division, Air Force Systems Command; Wright-Patterson Air Force Base: Dayton, OH, USA, 1967. [Google Scholar]
  41. Thielemann, W.; Esslinger, M. Beul- und Nachbeulverhalten isotropic Zylinder unter Äussendruck. Stahlbau 1967, 36, 161–174. [Google Scholar]
  42. Myachenkov, V.I. The stability of cylindrical shells under axisymmetric transverse pressure. Sov. Appl. Mech. 1970, 6, 19–23. [Google Scholar] [CrossRef]
  43. Zevin, A.A. Estimates of eigenvalue of self-adjont boundary-value problems with periodic coefficients. Ukr. Math. J. 1998, 50, 719–722. [Google Scholar] [CrossRef]
  44. Dell’Isola, F.; Romano, A. On the derivation of thermomechanical balance equations for continuous systems with a nonmaterial interface. Int. J. Eng. Sci. 1987, 25, 1459–1468. [Google Scholar] [CrossRef] [Green Version]
  45. Dell’Isola, F. Linear growth of a liquid droplet divided from its vapour by a “Soap Buble”–like fluid Interface. Int. J. Eng. Sci. 1989, 27, 1053–1067. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Corrugated ring with profile (4), h = 0.05 ;   n = 16 .
Figure 1. Corrugated ring with profile (4), h = 0.05 ;   n = 16 .
Symmetry 12 01250 g001
Figure 2. One pitch of corrugation.
Figure 2. One pitch of corrugation.
Symmetry 12 01250 g002
Figure 3. Projection of external pressure to the basic ring: p ,   p 0 are the intensities of external pressure acting on the corrugated and basic circular ring; θ is the angle between the tangents to the corrugated and the basic circlular ring; d s   and   d l are the arcs of the corrugated and basic circular ring, d l = d s   c o s θ ; P = p d s = p d l c o s θ ;   P n = P c o s θ = p d l ;   p 0 = p .
Figure 3. Projection of external pressure to the basic ring: p ,   p 0 are the intensities of external pressure acting on the corrugated and basic circular ring; θ is the angle between the tangents to the corrugated and the basic circlular ring; d s   and   d l are the arcs of the corrugated and basic circular ring, d l = d s   c o s θ ; P = p d s = p d l c o s θ ;   P n = P c o s θ = p d l ;   p 0 = p .
Symmetry 12 01250 g003
Figure 4. Ratio of buckling loads of corrugated (4) and basic circular rings p b 0 on parameter n h   ( h = 0.05 ; n = 0 100 ) .
Figure 4. Ratio of buckling loads of corrugated (4) and basic circular rings p b 0 on parameter n h   ( h = 0.05 ; n = 0 100 ) .
Symmetry 12 01250 g004
Figure 5. The part of the ring lying above the horizontal axis of symmetry: numbers 1 and 2 correspond to the middle lines of the basic circular ring before and after buckling; N ˜ ,   M ˜ ,   w ˜ ,   r ˜ are the internal force, moment, deflection, and radius at the points A and C.
Figure 5. The part of the ring lying above the horizontal axis of symmetry: numbers 1 and 2 correspond to the middle lines of the basic circular ring before and after buckling; N ˜ ,   M ˜ ,   w ˜ ,   r ˜ are the internal force, moment, deflection, and radius at the points A and C.
Symmetry 12 01250 g005
Figure 6. Projections of displacements on the basic circular ring. Here u ,   w are the tangential and normal displacements of the corrugated ring; u ,   w are the projections of these displacements on the basic circular ring; u = r u + r φ w A ; w = r w r φ u A .
Figure 6. Projections of displacements on the basic circular ring. Here u ,   w are the tangential and normal displacements of the corrugated ring; u ,   w are the projections of these displacements on the basic circular ring; u = r u + r φ w A ; w = r w r φ u A .
Symmetry 12 01250 g006
Figure 7. Decrease of the buckling pressure of a corrugated ring with an increase in the depth of the corrugation.
Figure 7. Decrease of the buckling pressure of a corrugated ring with an increase in the depth of the corrugation.
Symmetry 12 01250 g007
Figure 8. Comparison of buckling pressures p ¯ b 2   and   p ¯ b 2 e for corrugated and circular rings having the same volume of material.
Figure 8. Comparison of buckling pressures p ¯ b 2   and   p ¯ b 2 e for corrugated and circular rings having the same volume of material.
Symmetry 12 01250 g008
Figure 9. Comparison of buckling pressure obtained using the asymptotic homogenization method (36) and equivalent bending stiffness approach (7).
Figure 9. Comparison of buckling pressure obtained using the asymptotic homogenization method (36) and equivalent bending stiffness approach (7).
Symmetry 12 01250 g009
Figure 10. Diagrams of loading a corrugated ring with a corrugation (39): (a) n = 16 ,   h = 0.01 ; (b) n = 8 ,   h = 0.02 ;   p ¯ b 0 = 3 ,   p ¯ b 1 = 2.98 ,   p ¯ b 2 = 2.94 are buckling pressures for the smooth ring and for the corrugated ring obtained by the equivalent bending stiffness approach and asymptotic homogenization method.
Figure 10. Diagrams of loading a corrugated ring with a corrugation (39): (a) n = 16 ,   h = 0.01 ; (b) n = 8 ,   h = 0.02 ;   p ¯ b 0 = 3 ,   p ¯ b 1 = 2.98 ,   p ¯ b 2 = 2.94 are buckling pressures for the smooth ring and for the corrugated ring obtained by the equivalent bending stiffness approach and asymptotic homogenization method.
Symmetry 12 01250 g010
Table 1. Buckling loads (6) for the circular ring with the equivalent bending stiffness for constant value of nh.
Table 1. Buckling loads (6) for the circular ring with the equivalent bending stiffness for constant value of nh.
nh = 1.6n = 160;
h = 0.01
n = 80;
h = 0.02
n = 40;
h = 0.04
n = 20;
h = 0.08
n = 16;
h = 0.1
pb0.67670.67670.67670.67660.6764
Table 2. Comparison of the results of calculation of the deflection at the point φ = 0 for n = 16 ,   h = 0.01 ; here Δ i = w i 1 w 1 w i 1 100 ,   i = 1 4 .
Table 2. Comparison of the results of calculation of the deflection at the point φ = 0 for n = 16 ,   h = 0.01 ; here Δ i = w i 1 w 1 w i 1 100 ,   i = 1 4 .
Number of Iteration/Pressurep = −1p = −1.5p = −2
3Δ3 = 0.8%Δ3 = 1.09%Δ3 = 12%
4Δ4 = 0.01%Δ4 = 1.06%Δ4 = 3.18%

Share and Cite

MDPI and ACS Style

Andrianov, I.I.; Andrianov, I.V.; Diskovsky, A.A.; Ryzhkov, E.V. Buckling of Corrugated Ring under Uniform External Pressure. Symmetry 2020, 12, 1250. https://doi.org/10.3390/sym12081250

AMA Style

Andrianov II, Andrianov IV, Diskovsky AA, Ryzhkov EV. Buckling of Corrugated Ring under Uniform External Pressure. Symmetry. 2020; 12(8):1250. https://doi.org/10.3390/sym12081250

Chicago/Turabian Style

Andrianov, Igor I., Igor V. Andrianov, Alexander A. Diskovsky, and Eduard V. Ryzhkov. 2020. "Buckling of Corrugated Ring under Uniform External Pressure" Symmetry 12, no. 8: 1250. https://doi.org/10.3390/sym12081250

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