Published: 21 September 2023

Buckling of mono-axially compressed rectangular grids

Francesca Pancella1
Yuri De Santis2
Angelo Luongo3
1, 2, 3Department of Civil, Construction-Architectural and Environmental Engineering, University of L’Aquila, L’Aquila, Italy
Corresponding Author:
Francesca Pancella
Views 26
Reads 7
Downloads 77

Abstract

A theoretical analysis of the buckling of mono-axially compressed rectangular grid is carried out. The grid is composed of two orthogonal orders of continuous beams, simply supported at the ends. The critical load is determined by two different models: (i) the Kirchhoff’s homogeneous orthotropic plate, and (ii) a system of beams, solved by the Rayleigh-Ritz method. Parametric studies are carried out and numerical examples developed. The analytical solutions are validated by comparison with numerical Finite Element analyses.

1. Introduction

Grid structures are structural systems extensively used in many engineering applications for their many advantages: (a) they are efficient in transferring concentrated loads, (b) participate as a whole to the load carrying action, (c) are easy to inspect or repair because are open structures.

In the past, grid structures consisting of ribs and skins have been largely studied. Ref. [1] grouped analytical models for grid structures into two categories: exact models, which consider the grid geometry simulating ribs and skin individually, and equivalent models, which smear the grid as a homogenized plate. Grid structures can have different configuration: with or without skins (at the inner and outer surface of the shell), made of composite materials, multiple geometries, and shapes. In recent years, many studies have been carried out on the buckling behavior of reinforced cylindrical shells with cross stiffeners [2]. However, rectangular beam grids are the most common in practice and their buckling load is an important parameter to know.

In this work, the buckling behavior of rectangular beam grids is studied analytically. The critical load of mono-axially compressed rectangular beam grids is derived by using two different models: (i) the Kirchhoff's orthotropic plate, and (ii) a planar system of beams, both axially loaded in-plane. The two models are intrinsically different: the plate is meant as an equivalent model of homogenized continuum, in which the beams “are smeared” on the plate; the frame, instead, preserves the discrete nature of the constituent members. To solve the plate, the classic separation of the variables is adopted; to solve the frame, it is conjectured that the envelope of the buckled beams is a smooth function, so that a variational procedure can be applied.

Indeed, the Rayleigh-Ritz energy method had already been adopted long time ago [3], but it seems to have been forgotten nowadays. Here, the method of Ref. [3] is reproposed with slight changes. On the other hand, several orthotropic plate theories are available in literature [4]-[7]. The Kirchhoff's orthotropic model is chosen here for its simplicity.

The paper is organized as follows. In Section 2 the problem is posed. In Section 3 the Kirchhoff’s orthotropic theory is used for an equivalent plate, and the bifurcation condition is derived. In Section 4 the Rayleigh-Ritz energy method is worked out for the beam assembly, and the relevant critical load is computed. In Section 5 some parametric studies are carried out and the analytical solutions are validated by comparison with numerical Finite Element results. Finally, in Section 6 some conclusions are drawn.

2. Problem position

A rectangular l×b beam grid with rectangular Δx×Δy mesh is considered (see Fig. 1). The grid is composed of two orders of continuous beams, simply supported at the ends. The x-direction and y-direction beams are prestressed by forces Px and Py, respectively, positive when they generate compression. Therefore, in general, the grid is biaxially prestressed. p and q are the numbers of beams in x-direction and in y-direction, respectively. It is consequently l/Δx=q+1 and b/Δy=p+1. Each order of beams consists of equal beams of uniform section. The flexural and torsional stiffness of the single beam in i-direction are denoted by EIi and GJi, respectively.

Fig. 1Beam grid

Beam grid

3. Equivalent orthotropic plate model

A rectangular Kirchhoff’s orthotropic plate, simply supported on the whole contour is considered. The plate is monoaxially prestressed by forces px uniformly distributed and normal to the edges. Denoting by wx,y the out of plane displacement of the point x,y, the indefinite equilibrium equation is written as:

1
Dxw,xxxx+2Hw,xxyy+Dyw,yyyy+pxw,xx=0,

where Di i=x,y and H are the flexural and the torsional stiffnesses of the plate, respectively. They are evaluated by enforcing energy equivalences between the plate and the grid, and turn out to be equal to:

2
Dx=EIxΔy, Dy=EIyΔx, H=12GJxΔy+GJyΔx.

A linear combination with unknown coefficients of products of sinusoidal functions is used as trial functions, namely:

3
wx,y=n=1Nm=1MAnmsinnπlxsinmπby,

where Anm are unknown coefficients and n, m are integer numbers, even unknowns. Each term of the series describes a deformation in which the plate bends in n half-waves (of length l/n) in the x-direction and m half-waves (of length b/m) in the y-direction. Denoting by αl/b the aspect ratio and substituting Eq. (3) in Eq. (1), 2 critical loads are found, each associated with the pair n,m, i.e:

4
pxcn,m=Dyπ2b2DxDyn2α2+2HDym2+m4α2n2.

Since we are interested in the smallest of these loads, we need to minimize pxcn,m with respect to the two variables. Accordingly, m=1 must be taken, while n must be determined by trial and error. By letting δDx/Dy, γH/Dy and pxc=μxcπ2Dx/b2, the nondimensional critical load μxc is determined as:

5
μxc=n2α2+2γδ+α2δ n2.

It is observed that, when α is an integer, the minimum critical load occurs at n=α, and it is equal to μxc,min=1+2γδ+1δ. The same limit value is found for any α sufficiently large. When the plate is isotropic (i.e., when δ=γ=1), then, the well-known value μxc,min=4 is recovered.

4. Beam assembly model

The grid is considered as an assembly of beams undergoing flexure and torsion. It is conjectured that, in the buckled configuration, the beams lie on a smooth surface w=w(x,y), which is still given by Eq. (3), so that, by using the Raileigh-Ritz model, the problem can be solved in closed-form.

The Total Potential Energy (TPE) of the grid is defined as sum of flexural (f), torsional (t) and geometric (g) energy of the two orders of beams, Π=Πxf+Πyf+Πxyt+Πxg+Πyg.

By accounting for θx=w,y, θy=w,x, the single energy terms read:

6
Πxf=i=1p120lEIxw,xx2xi,yi-dx, Πyf=j=1q120bEIyw,yy2xj-,yjdy,
Πxyt=i=1p120lGJxw,yx2xi,yi-dx+j=1q120bGJyw,xy2xj-,yjdy,
Πxg=i=1p120lPxw,x2xi,yi-dx, Πyg=j=1q120bPyw,y2xj-,yjdy.

Assume, for simplicity's sake, that EIx=EIy=EI and GJx=GJy=GJ. By imposing stationarity of the TPE, the equilibrium equation follows, from which the buckling critical load is derived as:

7
Pxc=π2b2l2m2n2blGJbj=1qcos2j n Δxlπ+li=1pcos2i m Δybπm2l ρj=1qsin2j n Δxlπ+n2bi=1psin2i m Δybπ+EIm4l3j=1qsin2j n Δxlπ+n4b3i=1psin2i m Δybπm2lρj=1qsin2j n Δxlπ+n2bi=1psin2i m Δybπ,

where, ρPy/Px is the ratio between the magnitudes of the two loads and n, m are the values which minimize Pxc. In the special case of ρ=0, by introducing the dimensionless parameters ΔξΔx/l, ΔηΔy/b and by letting αl/b the aspect ratio of the grid, β=GJ/EI the torsional-to-flexural stiffness ratio, and Pxc=μxcπ2EI/b2, the critical nondimensional load μxc is determinated:

8
μxc=m2β1αj=1qcos2jn Δξπ+i=1pcos2imΔηπi=1psin2i m Δη π+m4n2αj=1qsin2jnΔξπ+n21α2i=1psin2i mΔη πi=1psin2i m Δη π.

It should be noticed that when m=p+1, Eq. (8) becomes indeterminate, because the flexural and prestress energies are zero (indeed, the beams are located at the nodal lines of the deformed configuration).

5. Numerical results

A parametric study is carried out by using the beam assembly model (Eq. (8)). Since μxc=μxc(ξ,η,α,β,n,m), as a first analysis, ξ, η are fixed and μxc (minimized with respect to n,m) plotted vs the aspect ratio α for different stiffness ratios β (Fig. 2(a)). The classic signature curves, typical of homogeneous and isotropic plates, is obtained. However, differently from those, it is observed that the critical load slowly decreases with α. This result can be interpreted by making use of the equivalent plate model. Indeed, since the nondimensional steps ξ, η are kept constants, also the number of beams parallel to the x and y-direction are constant, so that, while the aspect ratio increases, the equivalent stiffnesses per unit length of the plate decrease, entailing a decreasing of the critical load.

As a second analysis, only η is kept constant, while the Δξ=Δη/α is taken, thus obtaining the plot in Fig. 2(b). It is seen that the critical load now tends to an asymptotic value for large α. Indeed, the case analyzed is representative of a family of grids in which the width b is kept constant together with the number of beams in the x-direction, while the length l is increased, together with the number of beams in the y-direction, to keep their dimensional step Δx constant. Again, by resorting to the equivalent plate model, since the plate stiffnesses per unit length are equal for any member of the family, the tendency to an asymptotic value is consistently observed.

Fig. 2Frame model nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β=0,…,1: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α

Frame model nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β=0,…,1: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α

a)

Frame model nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β=0,…,1: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α

b)

The analytical results are compared with a Finite Element analysis for validation. The FE model is made of an assembly of Euler-Bernoulli beams. Constraints at the ends do not permit z-displacements for all the beams, and x- and y-displacement for beams in the y and x-directions, respectively. The loads are applied as concentrated forces at the ends of the beams, directed as the beam longitudinal axes (see Fig. 1). One of the (analytical) curves in Figs. 2(a) and 2(b) is compared in Figs. 3(a) and 3(b) with FE results, and a very good agreement is found. In the same figures, the results provided by the Kirchhoff’s orthotropic plate theory (Eq. (5)) are also reported; they turn out to be close to those relevant to the assembly beam model.

Fig. 3Nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β= 0.25: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α. Rayleigh-Ritz energy method (solid black line), Kirchhoff's orthotropic plate theory (dashed cyan line), FE analysis (dots)

Nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β= 0.25: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α. Rayleigh-Ritz energy method (solid black line), Kirchhoff's orthotropic plate theory (dashed cyan line), FE analysis (dots)

a)

Nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the aspect ratio α, for β= 0.25: a) Δξ=Δη=0.09; b) Δη=0.01, Δξ= Δη/α. Rayleigh-Ritz energy method (solid black line), Kirchhoff's orthotropic plate theory (dashed cyan line), FE analysis (dots)

b)

As a third analysis, a study is performed, in which the dimension b, l of the grid are kept constant, while the number of beams is increased. Fig. 4 compares the two analytical model with FE results. It refers to a square dominion (b= l) with an equal number of beams (p=q). It is seen that when the number of beams is sufficiently large (e.g., q=9, corresponding to Δη=0.1) the error is reasonable (about 6 %), and tends to zero with q. However, when q is smaller, the grid model gives a result better than the homogenized plate.

Fig. 4Comparison of results: a) nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the number of beams q=p, for β=0.25; b) percentage error (PE) from FE results. Rayleigh-Ritz energy method (blue), Kirchhoff's orthotropic plate theory (magenta), FE analysis (dots)

Comparison of results: a) nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the number of beams q=p, for β=0.25; b) percentage error (PE) from FE results. Rayleigh-Ritz energy method (blue), Kirchhoff's orthotropic plate theory (magenta), FE analysis (dots)

a)

Comparison of results: a) nondimensional critical load of the rectangular grid, subject to mono-axial compression vs the number of beams q=p, for β=0.25; b) percentage error (PE) from FE results. Rayleigh-Ritz energy method (blue), Kirchhoff's orthotropic plate theory (magenta), FE analysis (dots)

b)

6. Conclusions

A buckling analysis of a simply supported rectangular grid, subject to mono-axial compression, was carried out. Two models were used: (a) a homogenized Kirchhoff’s orthotropic plate, and (b) an assembly of beams tackled via an energy method. Analytical results were compared with Finite Element analyses. The following conclusions were drawn.

1) The Rayleigh-Ritz energy method supplies results in good agreement with FE models for all the cases analyzed.

2) The Kirchhoff's orthotropic plate model is also in good agreement with FE results only when the number of beams is sufficiently large.

The results should be important to the research, since based on analytical solutions that, unlike the numerical approaches, allow to study the influence of the parameters without redoing the model. The main contribution of this study is the interpretation of the grid mechanical behavior. Indeed, it has been proven here, that the surface enveloping the buckled grid is close to that describing the buckling of a homogeneous plate, while giving better results.

References

  • H.-J. Chen and S. W. Tsai, “Analysis and optimum design of composite grid structures,” Journal of Composite Materials, Vol. 30, No. 4, pp. 503–534, Mar. 1996, https://doi.org/10.1177/002199839603000405
  • M. Zarei, G. H. Rahimi, M. Hemmatnezhad, and F. Pellicano, “On the buckling load estimation of grid-stiffened composite conical shells using vibration correlation technique,” European Journal of Mechanics – A/Solids, Vol. 96, p. 104667, Nov. 2022, https://doi.org/10.1016/j.euromechsol.2022.104667
  • H. L. Cox and H. E. Smith, “The buckling of grids of stringers and ribs,” Proceedings of the London Mathematical Society, Vol. s2-48, No. 1, pp. 1–26, 1945, https://doi.org/10.1112/plms/s2-48.1.1
  • G. Z. Harris, “Buckling and Postbuckling of Orthotropic Plates,” AIAA Journal, Vol. 14, No. 11, pp. 1505–1506, Nov. 1976, https://doi.org/10.2514/3.61487
  • G. S. Johnston, “Buckling of orthotropic plates due to biaxial in-plane loads taking rotational restraints into account,” Fibre Science and Technology, Vol. 12, No. 6, pp. 435–443, Nov. 1979, https://doi.org/10.1016/0015-0568(79)90017-4
  • L. P. Kollar and I. A. Veres, “Buckling of rectangular orthotropic plates subjected to biaxial normal forces,” Journal of Composite Materials, Vol. 35, No. 7, pp. 625–635, 2001.
  • I. Hwang and J. S. Lee, “Buckling of orthotropic plates under various inplane loads,” KSCE Journal of Civil Engineering, Vol. 10, No. 5, pp. 349–356, Sep. 2006, https://doi.org/10.1007/bf02830088

About this article

Received
27 May 2023
Accepted
24 June 2023
Published
21 September 2023
SUBJECTS
Mathematical models in engineering
Keywords
beam grid systems
buckling analysis
Rayleigh-Ritz method
Kirchhoff’s orthotropic plate
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.