# Mathematical Identification of Influential Parameters on the Elastic Buckling of Variable Geometry Plate

Authors: Mirko Djelosevic; Jovan Tepic; Ilija Tanackov; Milan Kostelac
Journal: The Scientific World Journal
 (2013)

---

## Abstract

The problem of elastic stability of plates with square, rectangular, and circular holes as well as slotted holes was discussed. The existence of the hole reduces the deformation energy of the plate and it affects the redistribution of stress flow in comparison to a uniform plate which causes a change of the external operation of compressive forces. The distribution of compressive force is defined as the approximate model of plane state of stress. The significant parameters of elastic stability compared to the uniform plate, including the dominant role of the shape, size, and orientation of the hole were identified. Comparative analysis of the shape of the hole was carried out on the data from the literature, which are based on different approaches and methods. Qualitative and quantitative accordance of the results has been found out and it verifies exposed methodology as applicable in the study of the phenomenon of elastic stability. Sensitivity factor is defined that is proportional to the reciprocal value of the buckling coefficient and it is a measure of sensitivity of plate to the existence of the hole. Mechanism of loss of stability is interpreted through the absorption of the external operation, induced by the shape of the hole.


---

## Body

## 1. Introduction

Thin-walled plates are the basic components of a large number of modern supporting structures of open and closed type or profiled and box girder that are used in the construction of bridges, airplanes, boats, and other authorized facilities. Plate shaped elements of no uniform geometry are often incorporated in these structures which implies a variation in thickness and existence of various shapes holes. Requirements for using the plates with perforated holes and stepwise variable thickness result from functional, assembly, service, control, and other conditions, and often they are the result of weight minimization in order to achieve economic effect in the process of optimum design of structures. Mathematical identification of phenomenon of plate stability and general stability of girder and columns as their components is based on the application of the methods of Kantorovich and Ritz. The first method is based on solving partial differential equations of the fourth order (1), and it is derived from the equilibrium conditions of thin plate of constant geometry. The second method is based on the principle of minimum total potential energy, which provides a stable state of plate balance. Theoretical basis for both methods were systematically presented in [1–3]. The issue of stability of thin-walled elements occupies an important place in the process of supporting structures design [4].For plates of constant thickness affecting parameters include geometric sizes, supporting conditions, and load conditions [5–7], while the stability of the element with a hole in addition to these characteristics is affected by the shape, size, and position of the hole [8, 9].Most researches in the area of stability refer to the idealized case of a simply supported plate. In contrast to this is the plate with restraint edges along which the displacements and rotations are disabled. These two extreme and diametrically opposite cases practically very rarely occur in pure form, but they are most often occur in combination with other basic supporting conditions. In the case of real girders there is a case of elastic supported plate meaning that along the individual edges are beams of certain bending and torsional rigidity. Analysis of constant geometry plate shows that the supporting conditions greatly affect the stability [7]. In addition, implementation of staying is the way to increase stability [10–13].Researches of plates with holes and/or stepwise variable thickness as components of girder are of special interests [14–24]. It is important to emphasize that very important methods are those having universal application in the study of the stability of plates and shells. It refers to the elements of the support that may have more holes, usually arranged equidistant (perforation) as well as sudden or gradual changes in thickness.Rational design of structures requires the use of these structures which includes the identification of parameters influencing stability. In the last ten years, development of computer techniques based on the use of finite elements method (FEM) has enabled researchers to make a significant contribution to this phenomenon [25]. Researches [18, 22] show that the existence of the hole leads to a reduction of buckling coefficient for a certain value, which depends only on the shape, size, and position of the hole. To this the fact that the distance between the holes in perforated plates is also a parameter which in interaction with others affects the value of the critical buckling stress should be added. Due to the complexity of the mathematical model between the buckling coefficients, that is, minimum value of the critical buckling stress and characteristic parameters of the hole, researchers have mostly kept on semianalytic, approximate, and numerical methods and expressions obtained experimentally.


## 2. Analysis and Problem Formulation

A significant number of researchers have examined the phenomenon of linear or elastic buckling of plates indicating the importance of such an analysis in the design of structures. Studies of plate elements of girder or columns with perforated holes are of particular interest. These studies are realized using the finite element method, finite strip method [26–28], and experimental tests. Due to the geometric complexity and stress complexity the study of such problems mathematically is aimed to the partial analysis of individual components.Analysis of the elastic buckling of plates a number of researchers has relied on the application of von Karman's theory based on the fourth-order partial differential equation [1, 2, 29]:
(1)∇4w=-1D(Nx∂2w∂x2+Ny∂2w∂y2+2Nxy∂2w∂x∂y).The function of the deflection depends on the type of load and boundary conditions and it is presented as a linear combination of trigonometric and hyperbolic functions. As analyzed plates have ideally planar geometry without transversal pressure and initial curvature, integration constants that define the contour conditions remain indefinite.In this case, the deflection function with indefinite coefficients can be used to determine the coefficient of buckling. Under the influence of compression load plate is in the position of stable equilibrium until a critical intensity when it suddenly loses stability and fracture occurs. This situation is a consequence of the idealization of the initial state of plate geometry. The application of this method is limited to a plate of constant thickness.Numerous studies of linear buckling plates, especially those that are based on mathematical models, refer only to the problem of uniform plates. Study subject are plates of constant or variable thickness in one or both directions under the influence of combined plate loads. Analytical solution for pressure plate of variable thickness composed of two segments of the same length is given in [5]. The energy method based on the Rayleigh-Ritz principle is applied to analyze the plate [29, 30]. More complex cases are where the plates of gradual variable thickness in two normal directions are considered by the finite difference method [6].Thin-walled elements are potentially critical places of the structure due to the interactive effect between planar state of stress and elastic buckling phenomenon. The existence of holes in such elements followed by a stress concentration generally reduces their strength to buckling. Elastic buckling is the first step in the analysis of stability of structures. Measures taken to increase the critical stress of elastic buckling effect on favorable on behavior in the zone of inelastic buckling and on mechanism of boundary conditions.The subject of this paper simply supported uniaxial loaded rectangular plate with square, rectangular, and circular holes as well as slotted holes. Elastic buckling of plates with so-defined geometry and load conditions is modeled in terms of energy method (detailed information is given in Section3). The mathematical model is based on the principle of minimum total energy of plate that provides a state of stable equilibrium. Model of the elastic buckling of plate with a hole is tightly coupled with the model of a planar state of the considered element. The aim of the authors is to verify and complement the existing results regarding the elastic buckling of plates with various shapes holes. The analysis based on the physical and mathematical model enables exact identification of mechanism of elastic stability.


## 3. The Mathematical Model

The approach, in this paper, based on an analytical approach, requires forming a mathematical model according to defined physical conditions of supporting and load. Beside the model of elastic stability, it is necessary to form a model of plane state of stress of plate with rectangular hole. In this paper, the research is directed to analyze plates with holes of different shapes in the domain of elastic material behavior. Supporting structures are formed by different geometry plates interconnected into a functional entirely.The model formed on the basis of the plates system provides a real picture of the behavior of the whole cross section. Approach to solving was given in [31], where the problem of the static load capacity of box girder partially exposed to horizontal loads is analyzed. Partial analysis of stability of girder includes independent examination of certain plates [32], which greatly simplifies the mathematical model, which will be presented in this chapter.


### 3.1. Core of the Buckling Problem

The mathematical model is based on the application of energy method using total potential energy of plate. The most general application of the energy method refers to the case of transversely loaded plates and/or plates with initial curvature simultaneously exposed to the effects of plane forces. Transversal and plane forces represent known external load whose effect is manifested through external workEW, which leads to a deformation of plate. As a result of external load we have occurrence of deformation energy EI, manifested through the process of bending and torsion which leads to the induction of an additional plane forces due to the existence of connections of plate supports, if any. Transversal loading affects the bending of plate. Plane forces affect increase or reducing of the plate deflection in dependence on whether they are pressed or tightening. The existence of initial curvature is equivalent to the existence of a fictive transversal load [1].The mathematical model is formed based on the following assumptions:(i)
considered plates have ideal geometry (without significant warping);(ii)
behavior of isotropic plates is in the area of elastic stability;(iii)
the effect of supporting to inducing forces in a plate’s plane is ignored;(iv)
the impact of the external transversal loading is eliminated.Generally, the energy balance of plate, which is located in the state of stable equilibrium is(2)EW=EI.Work of external forces (EW) is defined by the following formula:
(3)EW=12∫0a∫0b[Nx(∂w∂x)2+Nx(∂w∂x)20000000000000+2Nxy∂w∂x∂w∂y(∂w∂x)2]dxdy-∫0a∫0bwqdxdy.The internal deformation energy (EI) is
(4)EI=D2∫0a∫0b{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)200000000000-2(1-υ)[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dxdy.A necessary and sufficient condition for plate in order to be in stable equilibrium position is that the total energyE has a minimum. Total energy of plate whose components are given by (3) and (4) and in pursuance of [1, 2] is
(5)E=EW+EIE=12∫0a∫0b[Nx(∂w∂x)2+Nx(∂w∂x)2000000000000+2Nxy∂w∂x∂w∂y(∂w∂x)2]dxdyEm+D2∫0a∫0b{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)2-2(1-υ)0000000000000×[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dxdyEm-∫0a∫0bwqdxdy.First term refers to the deformation energy or to work performed by forces acting in the plane of the plate. The induction of these forces is a consequence of boundary conditions or restrictions of plate; their effect has passive character, and it manifests itself by reducing deflection, and thus reducing the deformation energy plate, too. An example is limited supporting plate, that is, plate supported on fixed supports to prevent extension caused by the deflection of plate. The occurrence of passive plane forces is always connected to the phenomenon of deformation or internal energy and it is characteristics of geometry, material, and boundary conditions of plate. Conversely, if the plate is exposed to external planar load then their effect is active because initiate strain and the corresponding energy is classified as external work.


### 3.2. Approximate Model of Planar State of Stress of Plate with Rectangular Hole

To determine the work carried out by plane compression forces (3) it is necessary to know their distributions across the whole plate surface. If a plate is of constant thickness, with no holes, then the distribution of these forces is uniform and this case of load is the easiest case for studying stability. In case where there is a hole on the plate and/or stepwise variable thickness, disregarding the edges of plates exposed to constant planar load, there is a redistribution of these forces in the highly nonuniform functions. This phenomenon is explained by the existence of discontinuous changes through which the load cannot be transferred, but it is diverted to other areas of plates that are “bottlenecks” in the transfer of load. These are the most usually areas above and below the hole in whose immediate vicinity there is a concentration of the load, while in front of and back of the plate density of load stream lines decreases (Figure 1).Figure 1
Redistribution of stress flow of plate with a hole.Model of planar state of plate with rectangular hole (Figure2) is formed under the following conditions:(i)
distribution of load at place where elements of plate are separated is linearized;(ii)
the effect of the transversal forces along the connecting line among the elements of plate is ignored.Figure 2
Model of planar state of stress of plate with rectangular hole.Justification for these assumptions is reflected in the fact that linearized force functionNx has the same mean value as real distribution that the difference of the external work is a small value in comparison to the work of real load. External shearing forces do not affect the plate and the induced deformation energy due to the fact that shear stress is small value in comparison to the normal stress.The components of the reaction forceN1′(y) and N2′(y) defined by (6) and (7) are determined from the equilibrium condition of plate element “1” plates:
(6)N1′(y)=bb1b-b22b-b1-b2N=f1N,0≤y≤b1(7)N2′(y)=bb2b-b12b-b1-b2N=f2N,(b-b2)≤y≤b.Elasticity of plate affects the induction of additional reactions of element “1” which manifest as moments (8) and linear variable force (9):
(8)M1=s2e2+s4e1s2s3-s1s4N=g1N,M2=(e1s2+s1s2s2e2+s4e1s2s3-s1s4)N=g2N,M3=[t1r3-r2r3e1s2-s2e2+s4e1s2s3-s1s4(s1s2r2r3+r1r3)]N=g3N,M4=[t1r6-r8r6e1s2-s2e2+s4e1s2s3-s1s4(s1s2r8r6+r2r6)]N=g4N.The coefficientsr1,…,8, s1,…,4, e1,…,2, and t1,…,2 are defined in Appendix A.Distribution of linear variable forces, caused by the influence of the moments of elastic clamp, is given by(9)N1′′(y)=(2yb1-1)2M1b1=f1N,0≤x≤b1,N1′′(y)=(2(b-y)b2-1)2M2b2,(b-b2)≤x≤b.The distribution of reaction forces (6), (7), and (9) along x=a1 is stepwise variable, which means that at part b1<y<(b1+h) does not act load. Developing into a single trigonometric series and by their adding coming up are to unique function (12) which determines the conditions at the edge x=a1 of segment “1” of plate is obtained:
(10)N1(x=a1,y)=∑n=1∞[(-2bb1sinnπb12b)4f1Nnπsinnπb12bsinnπb12b-8g1Nn2π2b1200000000×cosnπ2(nπcosnπb12b00000000000000000-2bb1sinnπb12b)8g1Nn2π2b12]000000×sinnπyb=N∑n=1∞Pnsinnπyb,(11)N2(x=a1,y)=∑n=1∞[(nπcosnπb22b-2bb2sinnπb22b)4f2Nnπsinnπ(2b-b2)2bsinnπb22b0000000+8g2Nn2π2b22cosnπ(2b-b2)2b0000000×(nπcosnπb22b-2bb2sinnπb22b)nπ(2b-b2)2b]0000000×sinnπyb=N∑n=1∞Qnsinnπyb,(12)N(x=a1,y)=N1(x=a1,y)+N2(x=a1,y)=N∑n=1∞(Pn+Qn)sinnπyb.Defining a planar state of stress of plateAiry stress function is applied. This function has the following form:
(13)∂4U∂x4+∂4U∂x2∂y2+∂4U∂y4=0.A function that satisfies (13) is presented by the trigonometric series:
(14)U(x,y)=∑n=1∞[Ansinhnπxb+Bn(nπxb)sinhnπxb+Cncoshnπxb+Dn(nπxb)coshnπxb]×sinnπyb.Planar forcesNx and Nxy are defined through the component stresses σx and σxy formulated by
(15)Nx=σxδ=1δ∂2U∂y2Nx=π2δb2∑n=1∞{n2[(2coshnπxb+nπxbsinhnπxb)Ansinhnπxb000000000000000+Bn(2coshnπxb+nπxbsinhnπxb)000000000000000+Cncoshnπxb000000000000000+Dn(2sinhnπxb+nπxbcoshnπxb)]}000000000000×sinnπyb,Nxy=σxyδ=-1δ∂2U∂x∂yNxy=-π2δb2∑n=1∞{n2[(sinhnπxb+nπxbcoshnπxb)Ancoshnπxb000000000000000000+Bn(sinhnπxb+nπxbcoshnπxb)000000000000000000+Cnsinhnπxb000000000000000000+Dn(coshnπxb+nπxbsinhnπxb)]}00000000000000×cosnπyb.Integration constantsAn, Bn, Cn, and Dn are determined from the boundary conditions as follows:
(16)(1)Nx=N=2Nπ∑n=1∞1n(1-cosnπ)sinnπyb=const,0000000000000000000000000000000Nxy=0,x=0,(2)Nx=N∑n=1∞(Pn+Qn)sinnπyb,Nxy=0,x=a1.By substitution of (16) in (15) required coefficients defined by the coefficients (17) are obtained;
(17)An=An′N,Bn=Bn′N,Cn=Cn′N,Dn=Dn′N,
where
(18)αn=nπa1b.CoefficientsAn′, Bn′, Cn′, and Dn′ are given in Appendix B.The coefficients given by (17)-(18) are derived for element “1,” and for applying them to element “2” is necessary to correct factor αn by ratio a1/a2, and then its value is αn=nπa2/b. If the hole is not centric along the x-axis, that is, if the center of hole is nearer to one of the sides x=0 or x=b, then the elements of plates “1” and “2” are different, which is manifested through coefficient αn. Especially, if the hole is centric, then we have αn=nπa/2b.


### 3.3. Determining the Buckling Coefficient and Critical Stress

Analysis of elastic buckling of perforated plate by energy method is presented on the example of simply supported plate of ideal planar geometry, and without effect of transversal load, noting that this does not reduce the generality of the methodology. Deflection is assumed by trigonometric series (7), which satisfies the boundary conditions of a simply supported plate:
(19)w(x,y)=∑m=1∞∑n=1∞Amnsinmπxasinnπyb.The work carried out by uniaxial variable load (Nx) affecting the plate with a rectangular hole is
(20)EW=12∫ANx(∂w∂x)2dA=12∫0a1∫0bNx(∂w∂x)2dxdy22222+12∫a2a∫0bNx(∂w∂x)2dxdy22222+12∫a1a2∫0b1Nx(∂w∂x)2dxdy22222+12∫a1a2∫b2bNx(∂w∂x)2dxdy.The total external work is determined by the components that correspond to elements of plates “1–4,” according to Figure2:
(21)EW,1=12∫0a1∫0bNx(∂w∂x)2dxdyEW,1=12π4m2n2a2b2δAmn2EW,1m×(AnImn+BnJmn+CnKmn+DnLmn)EW,1=12π4m2n2a2b2δAmn2EW,1m×(An′Imn+Bn′Jmn+Cn′Kmn+Dn′Lmn)NEW,1=εmnAmn2N,EW,2=12∫a-a2a∫0bNx(∂w∂x)2dxdyEW,2=12π4m2n2a2b2δAmn2EW,2m×(EnVmn+FnWmn+GnRmn+HnTmn)EW,2=12π4m2n2a2b2δAmn2EW,2m×(En′Vmn+Fn′Wmn+Gn′Rmn+Hn′Tmn)NEW,2=ρmnAmn2N.Integral forms given by coefficientsImn, Jmn, Kmn, and Lmn which correspond to element “1” or by coefficients Vmn, Wmn, Rmn, and Tmn belonging to element “2” are formulated in Appendix C:
(22)EW,3=12∫a1a2∫b-b2bN1(y)(∂w∂x)2dxdyEW,3=14Amn2QnNm2bπna2(c+amπcos2mπpasinmπca)EW,3m×[13(cos3nπ-cosnπ(b-b2)b)00000000-(cosnπ-cosnπ(b-b2)b)(cos3nπ-cosnπ(b-b2)b)]EW,3=φmnAmn2N,EW,4=12∫a1a2∫0b1N2(y)(∂w∂x)2dxdyEW,4=14Amn2PnNm2bπna2(c+amπcos2mπpasinmπca)EW,4m×[13(cos3nπb1b-1)-(cosnπb1b-1)]EW,4=ωmnAmn2N.Adequate internal or deformation energy of plate consisted of bending and torsion of the components is given by(23)EI=D2∫A{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)200000000000-2(1-υ)[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dA.Flexion component of the internal energy is(24)EI,f=18DAmn2π4(m4a4+n4b4+2υm2n2a2b2)EI,fm×[ab-(c-amπcos2mπpasinmπca)00000000×(h-bnπcos2nπqbsinnπhb)]EI,f=λmnAmn2.Torsional component of the internal energy is:(25)EI,t=14D(1-υ)Amn2π4m2n2a2b2EI,tm×[ab-(c+amπcos2mπpasinmπca)0000000×(h+bnπcos2nπqbsinnπhb)]EI,t=μmnAmn2.The plate is in a state of stable equilibrium if the following condition is satisfied:(26)∂ET∂Amn=0,resp.∂(EI,f+EI,t+EW)∂Amn=0,
where the total potential energy of the plate is given by
(27)ET=EI,f+EI,t+EW.
Substitution of (20), (24), and (25) in (26) leads to (28) which defines the critical stress of elastic buckling:
(28)σcr=Ncrδ=λmn+μmnεmn+ρmn+φmn+ωmn.The critical buckling force can be written in compact form:(29)σcr=kπ2Db2δ,
where k is the coefficient of elastic buckling and it can be presented as
(30)khole=λmn+μmnεmn+ρmn+φmn+ωmnkholem×b2δπ2D.Ratio of buckling coefficient of uniform plate and plate with a hole is defined by sensitivity factort which is a measure of plate sensitivity to discontinuous changes in geometry. This factor is proportional to the reciprocal value of the coefficient of elastic buckling khole, and it is important because it represents the ability of plate with hole and a particular form and orientation, to keep its behavior in the domain of elastic stability:
(31)t=kunholekhole=4khole.Characteristic value of correction factort is 1, corresponding to a uniform plate (Figure 3(a)). If the plate is with hole this value is adjusted by coefficient (1/khole) and sensitivity factor t can be larger, smaller, or possibly equal to 1, depending on the values of influential parameters.Uniaxial stressed rectangular plate: (a) uniform geometry; (b) with a rectangular hole.
(a)
(b)


## 3.1. Core of the Buckling Problem

The mathematical model is based on the application of energy method using total potential energy of plate. The most general application of the energy method refers to the case of transversely loaded plates and/or plates with initial curvature simultaneously exposed to the effects of plane forces. Transversal and plane forces represent known external load whose effect is manifested through external workEW, which leads to a deformation of plate. As a result of external load we have occurrence of deformation energy EI, manifested through the process of bending and torsion which leads to the induction of an additional plane forces due to the existence of connections of plate supports, if any. Transversal loading affects the bending of plate. Plane forces affect increase or reducing of the plate deflection in dependence on whether they are pressed or tightening. The existence of initial curvature is equivalent to the existence of a fictive transversal load [1].The mathematical model is formed based on the following assumptions:(i)
considered plates have ideal geometry (without significant warping);(ii)
behavior of isotropic plates is in the area of elastic stability;(iii)
the effect of supporting to inducing forces in a plate’s plane is ignored;(iv)
the impact of the external transversal loading is eliminated.Generally, the energy balance of plate, which is located in the state of stable equilibrium is(2)EW=EI.Work of external forces (EW) is defined by the following formula:
(3)EW=12∫0a∫0b[Nx(∂w∂x)2+Nx(∂w∂x)20000000000000+2Nxy∂w∂x∂w∂y(∂w∂x)2]dxdy-∫0a∫0bwqdxdy.The internal deformation energy (EI) is
(4)EI=D2∫0a∫0b{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)200000000000-2(1-υ)[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dxdy.A necessary and sufficient condition for plate in order to be in stable equilibrium position is that the total energyE has a minimum. Total energy of plate whose components are given by (3) and (4) and in pursuance of [1, 2] is
(5)E=EW+EIE=12∫0a∫0b[Nx(∂w∂x)2+Nx(∂w∂x)2000000000000+2Nxy∂w∂x∂w∂y(∂w∂x)2]dxdyEm+D2∫0a∫0b{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)2-2(1-υ)0000000000000×[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dxdyEm-∫0a∫0bwqdxdy.First term refers to the deformation energy or to work performed by forces acting in the plane of the plate. The induction of these forces is a consequence of boundary conditions or restrictions of plate; their effect has passive character, and it manifests itself by reducing deflection, and thus reducing the deformation energy plate, too. An example is limited supporting plate, that is, plate supported on fixed supports to prevent extension caused by the deflection of plate. The occurrence of passive plane forces is always connected to the phenomenon of deformation or internal energy and it is characteristics of geometry, material, and boundary conditions of plate. Conversely, if the plate is exposed to external planar load then their effect is active because initiate strain and the corresponding energy is classified as external work.


## 3.2. Approximate Model of Planar State of Stress of Plate with Rectangular Hole

To determine the work carried out by plane compression forces (3) it is necessary to know their distributions across the whole plate surface. If a plate is of constant thickness, with no holes, then the distribution of these forces is uniform and this case of load is the easiest case for studying stability. In case where there is a hole on the plate and/or stepwise variable thickness, disregarding the edges of plates exposed to constant planar load, there is a redistribution of these forces in the highly nonuniform functions. This phenomenon is explained by the existence of discontinuous changes through which the load cannot be transferred, but it is diverted to other areas of plates that are “bottlenecks” in the transfer of load. These are the most usually areas above and below the hole in whose immediate vicinity there is a concentration of the load, while in front of and back of the plate density of load stream lines decreases (Figure 1).Figure 1
Redistribution of stress flow of plate with a hole.Model of planar state of plate with rectangular hole (Figure2) is formed under the following conditions:(i)
distribution of load at place where elements of plate are separated is linearized;(ii)
the effect of the transversal forces along the connecting line among the elements of plate is ignored.Figure 2
Model of planar state of stress of plate with rectangular hole.Justification for these assumptions is reflected in the fact that linearized force functionNx has the same mean value as real distribution that the difference of the external work is a small value in comparison to the work of real load. External shearing forces do not affect the plate and the induced deformation energy due to the fact that shear stress is small value in comparison to the normal stress.The components of the reaction forceN1′(y) and N2′(y) defined by (6) and (7) are determined from the equilibrium condition of plate element “1” plates:
(6)N1′(y)=bb1b-b22b-b1-b2N=f1N,0≤y≤b1(7)N2′(y)=bb2b-b12b-b1-b2N=f2N,(b-b2)≤y≤b.Elasticity of plate affects the induction of additional reactions of element “1” which manifest as moments (8) and linear variable force (9):
(8)M1=s2e2+s4e1s2s3-s1s4N=g1N,M2=(e1s2+s1s2s2e2+s4e1s2s3-s1s4)N=g2N,M3=[t1r3-r2r3e1s2-s2e2+s4e1s2s3-s1s4(s1s2r2r3+r1r3)]N=g3N,M4=[t1r6-r8r6e1s2-s2e2+s4e1s2s3-s1s4(s1s2r8r6+r2r6)]N=g4N.The coefficientsr1,…,8, s1,…,4, e1,…,2, and t1,…,2 are defined in Appendix A.Distribution of linear variable forces, caused by the influence of the moments of elastic clamp, is given by(9)N1′′(y)=(2yb1-1)2M1b1=f1N,0≤x≤b1,N1′′(y)=(2(b-y)b2-1)2M2b2,(b-b2)≤x≤b.The distribution of reaction forces (6), (7), and (9) along x=a1 is stepwise variable, which means that at part b1<y<(b1+h) does not act load. Developing into a single trigonometric series and by their adding coming up are to unique function (12) which determines the conditions at the edge x=a1 of segment “1” of plate is obtained:
(10)N1(x=a1,y)=∑n=1∞[(-2bb1sinnπb12b)4f1Nnπsinnπb12bsinnπb12b-8g1Nn2π2b1200000000×cosnπ2(nπcosnπb12b00000000000000000-2bb1sinnπb12b)8g1Nn2π2b12]000000×sinnπyb=N∑n=1∞Pnsinnπyb,(11)N2(x=a1,y)=∑n=1∞[(nπcosnπb22b-2bb2sinnπb22b)4f2Nnπsinnπ(2b-b2)2bsinnπb22b0000000+8g2Nn2π2b22cosnπ(2b-b2)2b0000000×(nπcosnπb22b-2bb2sinnπb22b)nπ(2b-b2)2b]0000000×sinnπyb=N∑n=1∞Qnsinnπyb,(12)N(x=a1,y)=N1(x=a1,y)+N2(x=a1,y)=N∑n=1∞(Pn+Qn)sinnπyb.Defining a planar state of stress of plateAiry stress function is applied. This function has the following form:
(13)∂4U∂x4+∂4U∂x2∂y2+∂4U∂y4=0.A function that satisfies (13) is presented by the trigonometric series:
(14)U(x,y)=∑n=1∞[Ansinhnπxb+Bn(nπxb)sinhnπxb+Cncoshnπxb+Dn(nπxb)coshnπxb]×sinnπyb.Planar forcesNx and Nxy are defined through the component stresses σx and σxy formulated by
(15)Nx=σxδ=1δ∂2U∂y2Nx=π2δb2∑n=1∞{n2[(2coshnπxb+nπxbsinhnπxb)Ansinhnπxb000000000000000+Bn(2coshnπxb+nπxbsinhnπxb)000000000000000+Cncoshnπxb000000000000000+Dn(2sinhnπxb+nπxbcoshnπxb)]}000000000000×sinnπyb,Nxy=σxyδ=-1δ∂2U∂x∂yNxy=-π2δb2∑n=1∞{n2[(sinhnπxb+nπxbcoshnπxb)Ancoshnπxb000000000000000000+Bn(sinhnπxb+nπxbcoshnπxb)000000000000000000+Cnsinhnπxb000000000000000000+Dn(coshnπxb+nπxbsinhnπxb)]}00000000000000×cosnπyb.Integration constantsAn, Bn, Cn, and Dn are determined from the boundary conditions as follows:
(16)(1)Nx=N=2Nπ∑n=1∞1n(1-cosnπ)sinnπyb=const,0000000000000000000000000000000Nxy=0,x=0,(2)Nx=N∑n=1∞(Pn+Qn)sinnπyb,Nxy=0,x=a1.By substitution of (16) in (15) required coefficients defined by the coefficients (17) are obtained;
(17)An=An′N,Bn=Bn′N,Cn=Cn′N,Dn=Dn′N,
where
(18)αn=nπa1b.CoefficientsAn′, Bn′, Cn′, and Dn′ are given in Appendix B.The coefficients given by (17)-(18) are derived for element “1,” and for applying them to element “2” is necessary to correct factor αn by ratio a1/a2, and then its value is αn=nπa2/b. If the hole is not centric along the x-axis, that is, if the center of hole is nearer to one of the sides x=0 or x=b, then the elements of plates “1” and “2” are different, which is manifested through coefficient αn. Especially, if the hole is centric, then we have αn=nπa/2b.


## 3.3. Determining the Buckling Coefficient and Critical Stress

Analysis of elastic buckling of perforated plate by energy method is presented on the example of simply supported plate of ideal planar geometry, and without effect of transversal load, noting that this does not reduce the generality of the methodology. Deflection is assumed by trigonometric series (7), which satisfies the boundary conditions of a simply supported plate:
(19)w(x,y)=∑m=1∞∑n=1∞Amnsinmπxasinnπyb.The work carried out by uniaxial variable load (Nx) affecting the plate with a rectangular hole is
(20)EW=12∫ANx(∂w∂x)2dA=12∫0a1∫0bNx(∂w∂x)2dxdy22222+12∫a2a∫0bNx(∂w∂x)2dxdy22222+12∫a1a2∫0b1Nx(∂w∂x)2dxdy22222+12∫a1a2∫b2bNx(∂w∂x)2dxdy.The total external work is determined by the components that correspond to elements of plates “1–4,” according to Figure2:
(21)EW,1=12∫0a1∫0bNx(∂w∂x)2dxdyEW,1=12π4m2n2a2b2δAmn2EW,1m×(AnImn+BnJmn+CnKmn+DnLmn)EW,1=12π4m2n2a2b2δAmn2EW,1m×(An′Imn+Bn′Jmn+Cn′Kmn+Dn′Lmn)NEW,1=εmnAmn2N,EW,2=12∫a-a2a∫0bNx(∂w∂x)2dxdyEW,2=12π4m2n2a2b2δAmn2EW,2m×(EnVmn+FnWmn+GnRmn+HnTmn)EW,2=12π4m2n2a2b2δAmn2EW,2m×(En′Vmn+Fn′Wmn+Gn′Rmn+Hn′Tmn)NEW,2=ρmnAmn2N.Integral forms given by coefficientsImn, Jmn, Kmn, and Lmn which correspond to element “1” or by coefficients Vmn, Wmn, Rmn, and Tmn belonging to element “2” are formulated in Appendix C:
(22)EW,3=12∫a1a2∫b-b2bN1(y)(∂w∂x)2dxdyEW,3=14Amn2QnNm2bπna2(c+amπcos2mπpasinmπca)EW,3m×[13(cos3nπ-cosnπ(b-b2)b)00000000-(cosnπ-cosnπ(b-b2)b)(cos3nπ-cosnπ(b-b2)b)]EW,3=φmnAmn2N,EW,4=12∫a1a2∫0b1N2(y)(∂w∂x)2dxdyEW,4=14Amn2PnNm2bπna2(c+amπcos2mπpasinmπca)EW,4m×[13(cos3nπb1b-1)-(cosnπb1b-1)]EW,4=ωmnAmn2N.Adequate internal or deformation energy of plate consisted of bending and torsion of the components is given by(23)EI=D2∫A{[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2](∂2w∂x2+∂2w∂y2)200000000000-2(1-υ)[∂2w∂x2∂2w∂y2-(∂2w∂x∂y)2]}dA.Flexion component of the internal energy is(24)EI,f=18DAmn2π4(m4a4+n4b4+2υm2n2a2b2)EI,fm×[ab-(c-amπcos2mπpasinmπca)00000000×(h-bnπcos2nπqbsinnπhb)]EI,f=λmnAmn2.Torsional component of the internal energy is:(25)EI,t=14D(1-υ)Amn2π4m2n2a2b2EI,tm×[ab-(c+amπcos2mπpasinmπca)0000000×(h+bnπcos2nπqbsinnπhb)]EI,t=μmnAmn2.The plate is in a state of stable equilibrium if the following condition is satisfied:(26)∂ET∂Amn=0,resp.∂(EI,f+EI,t+EW)∂Amn=0,
where the total potential energy of the plate is given by
(27)ET=EI,f+EI,t+EW.
Substitution of (20), (24), and (25) in (26) leads to (28) which defines the critical stress of elastic buckling:
(28)σcr=Ncrδ=λmn+μmnεmn+ρmn+φmn+ωmn.The critical buckling force can be written in compact form:(29)σcr=kπ2Db2δ,
where k is the coefficient of elastic buckling and it can be presented as
(30)khole=λmn+μmnεmn+ρmn+φmn+ωmnkholem×b2δπ2D.Ratio of buckling coefficient of uniform plate and plate with a hole is defined by sensitivity factort which is a measure of plate sensitivity to discontinuous changes in geometry. This factor is proportional to the reciprocal value of the coefficient of elastic buckling khole, and it is important because it represents the ability of plate with hole and a particular form and orientation, to keep its behavior in the domain of elastic stability:
(31)t=kunholekhole=4khole.Characteristic value of correction factort is 1, corresponding to a uniform plate (Figure 3(a)). If the plate is with hole this value is adjusted by coefficient (1/khole) and sensitivity factor t can be larger, smaller, or possibly equal to 1, depending on the values of influential parameters.Uniaxial stressed rectangular plate: (a) uniform geometry; (b) with a rectangular hole.
(a)
(b)


## 4. Comparative Analysis and Verification of Results

By analyzing the mathematical model of elastic buckling it can be concluded that the work of external load has a crucial role in the stability of plate. This conclusion comes to the fore when the plate with elongated holes is considered. Specifically, then the change of deformation energy of plate is minor while the existence of the hole initiate redistribution of stresses and significant change in the external energy. The consequence of this phenomenon is the correction of buckling coefficient or sensitivity factor, reducing or increasing their value. In this paper, an analysis of the influential parameters for some of the characteristic forms will be considered.


### 4.1. The Case of Rectangular Hole

Figure4 shows the comparative analysis of the buckling coefficient by presented method with data for the square plate whose dimensions are axa, with centric rectangular hole dimensions of cxh, where c=0,25a represents width of the hole. Commonly, the buckling coefficient k is expressed as a function of dimensionless value (h/a), where h is the height of the hole (Figure 3(b)). By analyzing the diagram (Figure 4) can be concluded that the data [9, 16, 22, 24] have the same distribution trend while the values differ.Figure 4
Comparative analysis of buckling coefficient of square plate with a rectangular hole.The main difference between presented procedure and literature is reflected in the shape of the distribution until the valuesy/a=0.25, where the function k(h/a) has a concave shape. According to [22] this form is linear, while other researches show that the function is convex in the whole domain. On the domain h/a>0.25 function of buckling coefficient asymptotic tends to the distribution according to [16], while its values are the closest to the data corresponding to [22]. The research [18] by analyzing the vertical slotted hole, which is closest in shape to a rectangular hole, shows the existence of the concave side in the field to h/a=0.25, as presented in the discussion on the slotted hole.A special case of the previous considerations is shown in Figure5, where the square plate with a centric square hole by dimensions of dxd is analyzed. Comparative analysis of this procedure is given in studies [14, 22]. By increasing dimensions of the hole in comparison with plate dimensions the buckling coefficient decreases permanently. The function k(d/a), the presented method, agrees well with the distribution corresponding to [22], both in terms of trend and quantity. According to [14] in the area above the d/a>0.4 the buckling coefficient k first stagnates and then tends to grow slowly.Figure 5
Comparative analysis of buckling coefficient of square plate with a square hole.In the special case when there is no rectangular hole we have the case of geometric uniform simply supported and uniaxial stressed plate (Figure3(a)) with buckling coefficient (32), corresponding to the data from [1, 2]:
(32)k=(mba+amb)2.The minimum value of the critical stress is obtained for(m,n)=1:
(33)σcr(m,n=1)=4π2Db2,
where k=4 is buckling coefficient of uniform square plate.


### 4.2. A Rectangular Plate with a Circular Hole

Determination of the coefficient of elastic buckling of a rectangular plate with a circular hole (Figure6) by previously presented procedure is not possible in a direct manner. The reason for this relates to the complexity of integrant due to circular contours and interdependent coordinates x and y. Therefore, it is necessary to use approximate solving approach that is based on the results of the previous analysis. Namely, the essence of this procedure is that the circular hole with diameter d is divided in s rectangular gaps with dimensions hi and ci whose positions are defined by pi and qi. The number of rectangular gaps has to be large enough in order to follow the contour of circular hole in a best way.Figure 6
A rectangular plate with a circular hole.In order to be applicable for a circular hole, it is necessary to make a correction of model of planar state of stress developed for rectangular hole, according to Figure7. Element “1” dimensions of a1xb is adjusted to the length of a1′ to which the smaller influence of moments corresponds (distribution tends to be uniform form).Figure 7
Distribution of the circular hole in rectangular gaps.Defined geometric values of rectangular gaps (34) are used in the expression for the deformation energy (24) and (25) instead of values p, c, q, and h:
(34)pi=p-d2+2i-12c,i=1,2,…,s,ci=di,i=1,2,…,s,qi=q=const,hi=d2-[(2i-1)c]2,i=1,2,…,s.The function of the buckling coefficientk(d/a) has a tendency to decline in value over the area of consideration (Figure 8). Progressive decline of value is in the interval d/a = (0.1–0.4) and values over 0.9. The obtained data are qualitatively and quantitatively consistent with [14], which justifies the assumptions in the process of model formation and its implementation to other similar contour shapes.Figure 8
Comparative analysis of buckling coefficient of square plate with a circular hole.


### 4.3. Rectangular Plate with a Slotted Hole

Research of plate with rectangular hole has shown that the orientation of the hole to the plate position has a great influence on the stability. On the other hand, analysis of plate with hole indicates an increase in work performed by an external forceNx causing a decrease of the buckling coefficient. Slotted hole represents a combination of the previous two cases.Analysis of buckling coefficientk indicates the dependence of the hole orientation. Figure 9(a) corresponds to a horizontal hole position in relation to the direction of the force Nx. Buckling coefficient according to [18], defined by FEM method, is linearly decreasing. The distribution of k(d/a) obtained by the presented methodology has the same trend as the maximum deviation of 5% for d/a=0.25 (Figure 10).A rectangular plate with a slotted hole (a) horizontal (b) vertical orientation.
(a)
(b)Figure 10
Comparative analysis of buckling coefficient of square plate with a slotted hole (horizontal position).However, if the hole is rotated by an angle of 90° (Figure9(b)) while the other features unchanged, the behavior of the plate is very different. The function k(d/a) can be divided into two parts in the domain of examination. The first part refers to the interval [0,0.2], where k is a concave function and it decreases regressively with increasing d/a. In the domain [0.2,0.5] function has convex form that affects primarily the stagnation and then to progressively increase. Functions of buckling coefficients for the case with slotted hole and rectangular hole are analogous. The value of k is larger for the rectangular hole in relation to the slotted hole for the same value of dimensionless argument as a result of lower external work (Figure 11).Figure 11
Comparative analysis of buckling coefficient of square plate with a slotted hole (vertical position).


## 4.1. The Case of Rectangular Hole

Figure4 shows the comparative analysis of the buckling coefficient by presented method with data for the square plate whose dimensions are axa, with centric rectangular hole dimensions of cxh, where c=0,25a represents width of the hole. Commonly, the buckling coefficient k is expressed as a function of dimensionless value (h/a), where h is the height of the hole (Figure 3(b)). By analyzing the diagram (Figure 4) can be concluded that the data [9, 16, 22, 24] have the same distribution trend while the values differ.Figure 4
Comparative analysis of buckling coefficient of square plate with a rectangular hole.The main difference between presented procedure and literature is reflected in the shape of the distribution until the valuesy/a=0.25, where the function k(h/a) has a concave shape. According to [22] this form is linear, while other researches show that the function is convex in the whole domain. On the domain h/a>0.25 function of buckling coefficient asymptotic tends to the distribution according to [16], while its values are the closest to the data corresponding to [22]. The research [18] by analyzing the vertical slotted hole, which is closest in shape to a rectangular hole, shows the existence of the concave side in the field to h/a=0.25, as presented in the discussion on the slotted hole.A special case of the previous considerations is shown in Figure5, where the square plate with a centric square hole by dimensions of dxd is analyzed. Comparative analysis of this procedure is given in studies [14, 22]. By increasing dimensions of the hole in comparison with plate dimensions the buckling coefficient decreases permanently. The function k(d/a), the presented method, agrees well with the distribution corresponding to [22], both in terms of trend and quantity. According to [14] in the area above the d/a>0.4 the buckling coefficient k first stagnates and then tends to grow slowly.Figure 5
Comparative analysis of buckling coefficient of square plate with a square hole.In the special case when there is no rectangular hole we have the case of geometric uniform simply supported and uniaxial stressed plate (Figure3(a)) with buckling coefficient (32), corresponding to the data from [1, 2]:
(32)k=(mba+amb)2.The minimum value of the critical stress is obtained for(m,n)=1:
(33)σcr(m,n=1)=4π2Db2,
where k=4 is buckling coefficient of uniform square plate.


## 4.2. A Rectangular Plate with a Circular Hole

Determination of the coefficient of elastic buckling of a rectangular plate with a circular hole (Figure6) by previously presented procedure is not possible in a direct manner. The reason for this relates to the complexity of integrant due to circular contours and interdependent coordinates x and y. Therefore, it is necessary to use approximate solving approach that is based on the results of the previous analysis. Namely, the essence of this procedure is that the circular hole with diameter d is divided in s rectangular gaps with dimensions hi and ci whose positions are defined by pi and qi. The number of rectangular gaps has to be large enough in order to follow the contour of circular hole in a best way.Figure 6
A rectangular plate with a circular hole.In order to be applicable for a circular hole, it is necessary to make a correction of model of planar state of stress developed for rectangular hole, according to Figure7. Element “1” dimensions of a1xb is adjusted to the length of a1′ to which the smaller influence of moments corresponds (distribution tends to be uniform form).Figure 7
Distribution of the circular hole in rectangular gaps.Defined geometric values of rectangular gaps (34) are used in the expression for the deformation energy (24) and (25) instead of values p, c, q, and h:
(34)pi=p-d2+2i-12c,i=1,2,…,s,ci=di,i=1,2,…,s,qi=q=const,hi=d2-[(2i-1)c]2,i=1,2,…,s.The function of the buckling coefficientk(d/a) has a tendency to decline in value over the area of consideration (Figure 8). Progressive decline of value is in the interval d/a = (0.1–0.4) and values over 0.9. The obtained data are qualitatively and quantitatively consistent with [14], which justifies the assumptions in the process of model formation and its implementation to other similar contour shapes.Figure 8
Comparative analysis of buckling coefficient of square plate with a circular hole.


## 4.3. Rectangular Plate with a Slotted Hole

Research of plate with rectangular hole has shown that the orientation of the hole to the plate position has a great influence on the stability. On the other hand, analysis of plate with hole indicates an increase in work performed by an external forceNx causing a decrease of the buckling coefficient. Slotted hole represents a combination of the previous two cases.Analysis of buckling coefficientk indicates the dependence of the hole orientation. Figure 9(a) corresponds to a horizontal hole position in relation to the direction of the force Nx. Buckling coefficient according to [18], defined by FEM method, is linearly decreasing. The distribution of k(d/a) obtained by the presented methodology has the same trend as the maximum deviation of 5% for d/a=0.25 (Figure 10).A rectangular plate with a slotted hole (a) horizontal (b) vertical orientation.
(a)
(b)Figure 10
Comparative analysis of buckling coefficient of square plate with a slotted hole (horizontal position).However, if the hole is rotated by an angle of 90° (Figure9(b)) while the other features unchanged, the behavior of the plate is very different. The function k(d/a) can be divided into two parts in the domain of examination. The first part refers to the interval [0,0.2], where k is a concave function and it decreases regressively with increasing d/a. In the domain [0.2,0.5] function has convex form that affects primarily the stagnation and then to progressively increase. Functions of buckling coefficients for the case with slotted hole and rectangular hole are analogous. The value of k is larger for the rectangular hole in relation to the slotted hole for the same value of dimensionless argument as a result of lower external work (Figure 11).Figure 11
Comparative analysis of buckling coefficient of square plate with a slotted hole (vertical position).


## 5. The Influence of the Hole on the External Work and Elastic Stability of Plate

Previously it was concluded that the work of external forces dominantly influence the elastic stability of plate. Deformation work of plate exclusively depends on the hole size, while the size, position, and orientation of the hole have secondary impact. Deformation work is a linear function of plate rigidity. Distributions of planar compressive forces cause induction of external work, so their identification is crucial in analysis of stability. Model of planar state of stress was developed for the rectangular hole; it was implemented with correction to circular and slotted form. The analysis carried out this section clearly shows that increasing of size of the one-dimensional hole (square and circular) leads to reduction in elastic stability. In the case of two-dimensional holes (rectangular shape and slotted) it should distinguish horizontal and vertical orientation. The first case is analogous to the previous analysis, while the second variant affects the elastic stability. It firstly decreases and then increases dominantly and even exceed the value of 4, which is characteristic for the plate without hole. Rectangular hole whose greater dimension is normal to the direction of load action is exposed to considerable stress concentration, which is negative in terms of static, but it is perfect condition for stability. Considering that the work of the external compression forces directly affects the buckling coefficient, it is necessary to introduce the coefficient of external workλ, which presents the absorbed energy caused by the shape of the hole in relation to the uniform plate.The coefficient of external work is defined as the ratio of external energy of plate with hole and uniform geometry, through relation(35)λ=EW,holeEW,unhole.Value of coefficientλ for slotted hole in the horizontal position is greater than 1 in the interval from 0 to 0.5 (Figure 12). This value decreases for larger domains. Coefficient λ for the other shapes of hole has a tendency to fall, since their shape affects the reduction of work of compressive forces in relation to the uniform plate. The minimum value of λ is for the vertical orientation of slotted hole.Figure 12
Change of coefficientλ in relation to the parameter (h/a) for the considered hole.There is certain similarity between the diagrams shown in Figures12 and 13. Slotted hole in horizontal position has the most sensitivity. This form redistributes the load above and below the hole over nominal value of N, while in front of and behind the hole circular configuration prevents turbulence of fluid flow and distribution near the value Nx is formed. Because the value of Ew is greater than the value corresponding to the uniform plate, where the pressure distribution along the plate corresponds to the nominal value of Nx. The situation of plate with rectangular hole is opposite to the above mentioned. As a result, we have maximum sensitivity for horizontal orientation slotted shape, while the same shape rotated by 90° has a minimum. Other variations are between these two cases.Figure 13
Sensitivity factort as a function of the (h/a) for the holes of the considered forms.The diagram given in Figure13 is a basic guideline for selection of the hole shape in terms of the plate stability. Vertical position of slotted hole (Figure 9(b)), in this respect represent of optimal shape, which is important in the construction process of thin-walled supporting structures.


## 6. Conclusions

In this paper the problem of elastic stability of plate with a hole using a mathematical model developed on the basis of the energy method has been addressed. The analysis includes consideration of four types of holes: rectangular, square, circular, and slotted (the combination of rectangular and circular). The existence of the hole reduces the deformation energy of plate, depending on the form which most commonly affects the decrease in external work due to compressive forces. The dominant factor on the buckling coefficientk is the work of external forces. That is the reason why coefficient λ which manifests absorption ability of plate with the hole in relation to the uniform plate is defined. Through the exposed methodology analyzed mechanism of elastic stability of plates using energy balance in a state of stable equilibrium is analyzed. The mathematical identification of influential parameters on stability which include dimensions and type of plate material, shape, size, position, and orientation of the hole was carried out.Verification of results of the presented method has been verified by studies [9, 14, 16, 18, 22, 24]. Sensitivity factor t defines criteria for selection of the best form of hole in terms of stability. Results of the application of this study are important for the optimization of the structures, with special emphasis on thin-walled erforated elements of high bay warehouses in order to increase the economy of distribution of supply chains. The developed model (Section 3) is applicable to the methodology [31] and it can be implemented in order to further studying of the phenomenon of buckling of thin-walled structures.



---

*Source: 268673-2013-11-24.xml*