Table of Contents
ISRN Mechanical Engineering
Volume 2011, Article ID 979532, 9 pages
Research Article

Orthotropic Models of Corrugated Sheets in Finite Element Analysis

Department of Aeronautical and Vehicle Engineering, The Royal Institute of Technology (KTH), Teknikringen 8, 100 44 Stockholm, Sweden

Received 28 January 2011; Accepted 21 March 2011

Academic Editors: F. Findik and A. Tounsi

Copyright © 2011 David Wennberg et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


To reduce computational effort of finite element (FE) calculations a corrugated sheet is replaced with an orthotropic plate. Analytical expressions for the mechanical properties are studied and compared to finite Element calculations in extension, free vibration, and buckling. Good similarity is shown in the stiffened and transverse direction of the corrugated sheet; however, the orthotropic models do not give an accurate twisting behavior. The stiffened direction of the corrugated sheet best matches the analytical expressions. Keeping in mind the presented limitation, the orthotropic model presented herein can be used to drastically reduce the number of elements needed when modelling corrugated sheet with finite elements.

1. Introduction

Modeling corrugated sheets using the finite element method (FEM) usually requires an extensive amount of elements. To reduce the required amount of elements, 2-dimensional, orthotropic models may be used to represent the properties of a corrugated sheet. Samanta and Mukhopadhyay [1], Briassoulis [2], and Liew et al. [3] have derived analytical expressions for the extensional and flexural rigidities of a corrugated sheet. These expressions are here benchmarked against each other with finite element analysis (FEA) in extension, modal analysis, checking both frequency and mode shape, and buckling analysis. Furthermore, a 3-dimensional corrugated sheet model is used as a reference.

The specific construction studied in this paper is the structural floor of a rail vehicle. It is designed of corrugated sheet metal and a metal framework as the load carrying structure. To reduce computational effort, and enable fast parameter studies of the entire rail car, the original FE model has to be reduced.

Another goal of this exercise is to find a more efficient substitution to the existing structure (efficient here meaning lighter with the same mechanical properties). A possible competitor to the corrugated sheet and frame construction may be a sandwich panel structure. The expressions tested and formulated in this paper may, as well as simplifying FE modeling, also be used as design parameters in a sandwich-panel material and geometry selection process.

In this paper, parameters without subindices are the regular, isotropic parameters for the material in question. Sub-indexed parameters refer to the orthotropic model, and indices are with respect to directions illustrated in Figure 1.

Figure 1: Cross-section of one repeated corrugation of a larger corrugated sheet.

2. Method

The corrugated sheet was studied in free vibration, extension, and buckling. FE calculations were performed with the software package HyperWorks, where HyperMesh was used as a preprocessor, RADIOSS as the solver, and HyperView as the postprocessor. To create an orthotropic equivalent of the corrugated sheet in HyperMesh, the PSHELL element property was used and the MID1 (membrane) and MID2 (bending) material options in the PSHELL properties were activated [4]. Material properties for MID1 and MID2 were defined by the MAT2 constitutive material matrix presented below,𝜎1𝜎2𝜏12=𝐺11𝐺12𝐺13𝐺21𝐺22𝐺23𝐺31𝐺32𝐺33𝜖1𝜖2𝛾12.(1) For each orthotropic model, two such constitutive matrices are thus calculated, one for bending and one for the membrane properties of the corrugated sheet.

The comparison of free vibration properties was performed using four FE models of a large, unsupported sheet. One model was a thorough, 3-dimensional representation of the corrugated sheet, and the other three models were of orthotropic plates with equivalent flexural rigidities calculated according to Samanta and Mukhopadhyay [1], Briassoulis [2], and Liew et al. [3], see Figure 2. These rigidities were used to calculate the input for the MID2 material properties defined above. Both the natural frequencies and the mode shapes were compared between the different models.

Figure 2: Models used in modal analysis and extensional test. Geometric parameters from the corrugated sheet are used in analytical expressions to calculate equivalent orthotropic properties for the 2D orthotropic plate.

Extensional rigidities were calculated according to Samanta and Mukhopadhyay [1] and Briassoulis [2]. These values were benchmarked in extension, again with a fully modeled corrugated sheet as a reference.

For the comparison of buckling behavior, four models were created, cf. Figure 3. Model B:1 is a supported corrugated sheet which was loaded along two edges. Buckling loads for this model were calculated by means of finite element analysis (FEA). Model B:2 is an equivalent orthotropic plate of the same size as B:1, loaded and constrained in the same way. Model B:3, which illustrates a representative part of a corrugation, was also modeled with FE and compared to analytical Euler buckling loads and local buckling loads calculated from Model B:4. Material properties for the orthotropic plate model were chosen as the material properties that best matched the fully corrugated model from the previous two tests in free vibration and extension.

Figure 3: Models used in buckling analysis. In B:1 it is the section 𝐿𝑥 by 𝐿𝑦 which is studied in buckling, not the entire structure. Geometric parameters from B:1 are used in analytical expressions to calculate equivalent orthotropic properties for B:2. B:3 is a representative part of B:1, this model is compared to the analytical result of Euler buckling found in B:4.

An overview of the analyses performed can be found in Table 1.

Table 1: Summary of the different analyses performed, cf. Figures 2 and 3.
2.1. Orthotropic Plates

Assuming we have large corrugated sheets, that is, the corrugation size is small in comparison to the overall size of the sheet, we may be able to model the sheet as a thin orthotropic plate with the stress-strain relationship as follows:𝜎𝑥𝜎𝑦𝜏𝑥𝑦=11𝜈𝑥𝑦𝜈𝑦𝑥𝐸𝑥𝜈𝑦𝑥𝐸𝑥0𝜈𝑥𝑦𝐸𝑦𝐸𝑦0001𝜈𝑥𝑦𝜈𝑦𝑥𝐺𝑥𝑦𝜖𝑥𝜖𝑦𝛾𝑥𝑦.(2) For an arbitrary orthotropic, thin plate, the flexural rigidities may be expressed as𝐷𝑥=𝐸𝑥𝑡3121𝜈𝑥𝑦𝜈𝑦𝑥,𝐷𝑦=𝐸𝑦𝑡3121𝜈𝑥𝑦𝜈𝑦𝑥,𝐷𝑥𝑦=𝐺𝑥𝑦𝑡36.(3) Furthermore, for a corrugated sheet, the Poisson's ratio 𝜈𝑦𝑥, is the same as that of the isotropic material [2], that is, 𝜈𝑦𝑥=𝜈. To find 𝜈𝑥𝑦, we utilise the relation [3]𝜈𝑥𝑦𝜈𝑦𝑥=𝐸𝑥𝐸𝑦=𝐷𝑥𝐷𝑦.(4) Rearranging (4) gives𝜈𝑥𝑦=𝐷𝑥𝐷𝑦𝜈𝑦𝑥=𝐷𝑥𝐷𝑦𝜈.(5) Utilising (3) to express the Young's modulus as a function of the bending rigidity, we get𝐸𝑥=121𝜈𝑥𝑦𝜈𝑦𝑥𝑡3𝐷𝑥,𝐸𝑦=121𝜈𝑥𝑦𝜈𝑦𝑥𝑡3𝐷𝑦,𝐺𝑥𝑦=6𝐷𝑥𝑦𝑡3.(6)

The extensional rigidities and moduli of an orthotropic plate (here denoted 𝐵 and 𝐸𝑒 to distinguish from the flexural expressions) follow the relations [2]𝐵𝑥=𝐸𝑒𝑥𝐵𝑡,𝑦=𝐸𝑒𝑦𝐵𝑡,𝑥𝑦=𝐺𝑒𝑥𝑦𝑡.(7)

3. Corrugated Model

The model used in this study is 12.164 m long (𝑦-dir) and 2.336 m wide (𝑥-dir), which is a significant part of a rail vehicle's floor area, see Figure 4 for definition of directions. The model weight is 475.1 kg. Remaining properties can be found in Table 2.

Table 2: Material data and parameters for the studied corrugated sheet. Parameters can be found in Figure 1.
Figure 4: Directions in the corrugated sheet.

The density used in the orthotropic models was estimated as (𝑙/𝑐)𝜌, where 𝜌 is the density of steel. The thickness of the orthotropic plate is equal to that of the corrugated sheet, 𝑡.

4. Flexural Properties

To derive the expressions for the flexural moduli, Liew et al. [3] used the following expressions for the flexural rigidities:𝐷𝑥=𝐸𝑡3121𝜈2𝑐𝑙𝐷,(8)𝑦=𝐸𝑡3121𝜈2+𝐸𝑡𝑐𝐷𝛼,(9)𝑥𝑦=𝐸𝑡312(1+𝜈),(10) where 𝐸 is the Young's modulus of the material and 𝜈 is Poisson's ratio; the other parameters, all but 𝛼, are geometric variables described in Figure 1. The parameter 𝛼 is given by [3]𝑓𝛼=33tan𝜃+𝑓21𝑏+3tan2𝜃𝑐3𝑓𝑏+tan𝜃3𝑐(2𝑓+𝑏tan𝜃)tan𝜃2𝑓𝑏+tan𝜃2+(2𝑓+𝑏tan𝜃)2𝑓𝑐𝑏.tan𝜃(11)

Equations (8)–(10) were derived by Briassoulis [2]. Liew et al. have, however, made changes to the calculation of 𝐷𝑦. Briassoulis original expression for 𝐷𝑦 is given by𝐷𝑦=𝐸𝑡3121𝜈2+𝐸𝑡𝑓22.(12) Briassoulis' equations are modified equations from, for example, Easley [5].

A third variant of the flexural rigidities is given by Samanta and Mukhopadhyay [1]𝐷𝑥=𝑐𝑙𝐸𝑡3𝐷12,(13)𝑦=𝐸𝐼𝑦𝑐𝐷,(14)𝑥𝑦=𝑙𝑐𝐸𝑡3.6(1+𝜈)(15) These equations are similar to the ones used by Briassoulis, and, furthermore, Samanta and Mukhopadhyay and Briassoulis have referenced the same original authors [57].

The three orthotropic models (Liew's, Samanta's, and Briassoulis') were compared to a 3D model of the corrugated plate in free vibration.

4.1. Flexural Analysis

Bending rigidities calculated with the three methods presented in Section 4 can be found in Table 3. The elasticity moduli corresponding to the bending rigidities found in Table 3 are presented in Table 4. These values were used to calculate the constitutive relation in (1), which was used to describe the orthotropic models' bending characteristics for free vibration analysis.

Table 3: Comparison of flexural rigidities calculated with methods described in Section 1.
Table 4: Orthotropic moduli calculated with (6) and bending rigidities found in Table 3.

The first six eigen modes of the reference model, that is, the corrugated model, are depicted in Figure 5.

Figure 5: Modal analysis of corrugated sheet, first six eigen modes.

Results from the modal analysis of the three orthotropic models in comparison to the reference model are found in Table 5.

Table 5: Comparison of the first six eigen modes of the orthotropic models and the original model. The different types refer to the type of eigen mode defined in the original model, for example, type b is the first bending mode around the 𝑥-axis, cf. Figure 5(b). The reduction of frequency compared to the original model is represented by Δ.

5. Extensional Properties

For extensional rigidities, both Samanta and Mukhopadhyay [1] and Briassoulis [2] have derived similar equations. Samanta and Mukhopadhyay have presented the following expression:𝐵𝑥=𝐸(𝑡/𝑓)2𝑡6[],𝐵(2/𝑐)((𝑐𝑙cos𝜃)/(1cos𝜃))+(4𝑓/3𝑐sin𝜃)𝑦𝑙=𝐸𝑡𝑐,𝐵𝑥𝑦=𝑐𝑙𝐸𝑡2(1+𝜈)(16) while Briassoulis has derived different versions of 𝐵𝑥 and 𝐵𝑥𝑦 (see below)𝐵𝑥=𝐸𝑡1+(𝑓/𝑡)261𝜈2𝑙2/𝑐2,𝐵(𝑙/2𝜋𝑐)sin(2𝜋𝑙/𝑐)𝑦𝑙=𝐸𝑡𝑐,𝐵𝑥𝑦=𝐸𝑡.2(1+𝜈)(17) An extensional benchmark calculation between a fully modeled corrugated metal sheet and orthotropic equivalents with values according to Briassoulis and Samanta and Mukhopadhyay was performed with FEA. A tensional load of 10 kN was applied along each side of the sheets, separately, in two calculations. The models, the orthotropic and corrugated, are oriented as shown in Figure 4, that is, 𝑦 is along the corrugations and 𝑥 is across. The loading condition is depicted in Figure 6.

Figure 6: Loading and boundary conditions for two extensional benchmark tests of corrugated and orthotropic plate. Dashed lines are simply supported boundary.
5.1. Extensional Loading

Extensional moduli for the methods presented by Samanta and Mukhopadhyay [1] and Briassoulis [2], cf. Section 5, are presented in Table 6.

Table 6: Extensional moduli comparison between methods found in [1, 2].

The results from the FE extensional benchmark comparison between the method presented by Briassoulis, Samanta, and the corrugated reference can be found in Table 7.

Table 7: Strain of corrugated sheet and orthotropic models due to 10 kN loading in 𝑥 and 𝑦 directions, cf. Figures 6 and 4.

6. Buckling

A suitable replacement for the corrugated sheet must also have a sufficient buckling strength. Therefore buckling characteristics of the corrugated sheet and the orthotropic model are studied and compared. Buckling analysis was performed on a part of the corrugated sheet with properties according to Table 2 and dimensions 2271.5(𝐿𝑥) by 746.6(𝐿𝑦) mm. This is equivalent to a part of the corrugated sheet between two support beams in the floor of the studied rail vehicle, cf. Figure 3. The plate was subjected to compression in the corrugated direction, that is, the 𝑦-direction. Long edges were clamped to simulate a continuous plate attached to the structural beams, and short edges were simply supported, cf. Figure 7.

Figure 7: Illustration of corrugated section subjected to buckling analysis, see also Figure 3. Dashed lines are clamped and dotted lines are simply supported.

An orthotropic model was also studied. Here values were chosen as the ones that showed best coherence to the fully corrugated model in the previous two tests, that is, flexural properties according to Samanta and Mukhopadhyay and extensional properties according to Briassoulis.

To derive analytical expressions for the buckling stress, half a corrugation was studied analytically and by means of FE analysis. Modeled with appropriate boundary conditions, a half corrugation will buckle at approximately the same buckling-stress as the entire sheet. Geometry of a half corrugation beam is illustrated in Figure 8. Boundary conditions applied during FE analysis can be found in Table 8. Euler buckling and local buckling was calculated analytically.

Table 8: Boundary conditions for half-corrugation buckling analysis.
Figure 8: Geometry of half a corrugation of the corrugated sheet presented in Figure 1.

The analytical buckling load, for Euler buckling, is given by [8]𝑃𝐸cr=𝜋2𝐸𝐼𝑙2𝑒,(18) where 𝑙𝑒 is the effective length of the column and is defined as half the length of the column for a case with both edges clamped, cf. Megson, [8, page 258]. 𝐼, which in this case may be written 𝐼𝑥, is the second moment of area of the cross section defined as𝐼𝑥=𝐶𝑧2𝑑𝐴.(19) The critical stress is then calculated as𝜎cr,half=𝑃𝐸cr𝐴𝑥𝑧,half,(20)

where 𝐴𝑥𝑧,half is the cross-section area of the half-corrugation.

The analytical local buckling stress is given by [8]𝜎cr,local=𝑘𝜋2𝐸121𝜈2𝑡𝑙𝑏2𝑡=3.6𝐸𝑙𝑏2,(21) where 𝑘 is the buckling coefficient and 𝑙𝑏 is the shorter side length of the plate, in this case the loaded side lengths 13.5 mm and 27.5 mm. For the half corrugation it will be the web that is most critical.

6.1. Buckling Results

The first buckling mode of the small corrugated sheet section situated between structural beams, cf. Figure 7, is shown in Figure 9. Buckling occurred at a total load of about 6.21 MN. The total cross sectional area of this sheet is 𝐴=4823mm2.(22) This gives a critical stress of𝜎cr=𝑃cr𝐴=6.211064823=1288MPa.(23)

Figure 9: First buckling mode of corrugated sheet, critical load 6.21 MN.

The first buckling mode of the orthotropic plate model corresponded to a load of 7.16 MN. Since this model was created with the same sheet thickness as the sheet metal used in the fully corrugated model, that is, 1.5 mm, we get a significantly higher buckling stress (the total cross-section area is a factor 𝑙/𝑐 smaller)𝜎cr=𝑃cr𝐴=7.161061.52271.5=2101MPa.(24)

Corrected with the factor 𝑙/𝑐, we get a comparable buckling stress of 1484 MPa.

The first critical buckling load of the half corrugation model presented in Section 6 was calculated to 111.8 kN; the buckling mode is illustrated in Figure 10.

Figure 10: First buckling mode of half corrugation.

The area of this cross sections is𝐴𝑥𝑧,half=82mm2.(25)

This gives us a critical stress of𝜎cr,half=𝑃cr,half𝐴𝑥𝑧,half=111.810382=1363MPa.(26)The moment of inertia for the half corrugation beam is𝐼𝑥=8487mm4.(27)

Furthermore, 𝐸=200 GPa, and 𝑙=746.6 mm. This gives an analytical Euler buckling load, according to (18), of𝑃𝐸cr,half=𝜋22.01058487(0.5746.6)2=120kN.(28)

The finite element analysis gave a first local buckling mode, for the half corrugation model, at a compressive load of 173 kN, which equals a critical stress of 2110 MPa. For the section of the corrugated sheet between beams, cf. Figure 7, local buckling occurred at a force of 9.60 MN; this equals a critical stress of 1990 MPa. The analytical local buckling stress for the web of the half corrugation, cf. (21), was calculated to𝜎cr,local𝑡=3.6𝐸𝑙𝑏2=2140MPa.(29)

The orthotropic model does not have the same local buckling mode due to the change of geometry.

A summary of the critical stresses for buckling of the corrugated sheet, orthotropic model, the half corrugation, and the analytical calculations can be found in Table 9.

Table 9: Critical buckling stresses calculated for the corrugated sheet section situated between structural beams as described in Section 6 and Figure 7, as well as the critical stresses for a half corrugation beam, and analytically calculated references.

7. Discussion

The orthotropic models suggested by Samanta and Mukhopadhyay [1] and Liew et al. [3] match the results from the modal analysis on the reference model with regard to order or modal shape. The frequencies, all except the first, also match the reference model, especially the ones calculated with the model suggested by Samanta and Mukhopadhyay, which differ with an average of 3% from the reference model (Note: first mode not considered), cf. Table 5.

For the first eigen mode, the method presented by Briassoulis and Liew et al. gives better results than the method by Samanta and Mukhopadhyay. This mode is largely influenced by the twisting rigidity 𝐷𝑥𝑦. One may conclude from this that Briassoulis expression for 𝐷𝑥𝑦 is better then the expression for 𝐷𝑥𝑦 presented by Samanta and Mukhopadhyay.

The modal order has changed for the model presented by Briassoulis [2] in comparison to the other models, cf. mode frequencies for type d and e in Table 5.

The frequency of the eigen mode type b, is similar for all orthotropic models because they all have similar values for 𝐸𝑥, see Table 4. The difference in bending rigidities between Samanta and Briassoulis/Liew lies in the simplification: 𝜈20 made by Samanta, cf. (13) and (8).

Mode d is coupled to the bending rigidity 𝐷𝑦, Young's modulus 𝐸𝑦, see Figure 5(d). Here, Liew et al. [3] seem to have made improvements on Briassoulis' equation. However, the best result is given by Samanta and Mukhopadhyay's expressions [1].

The extensional moduli in 𝑦-direction, that is, 𝐸𝑒𝑦, derived by Briasssoulis and presented in Table 6, correspond well to the extensional property of the corrugated sheet in this direction, cf. the extensional benchmark results in Table 7.

The expressions for 𝐵𝑥 given by Briassoulis and Samanta seem to under estimate the extensional rigidity. However, Briassoulis' expressions are the better alternative, at least for this specific case.

The critical buckling stresses of the half corrugation model and the analytical calculation match each other rather well, cf. Table 9. The difference in the first critical buckling stress may be explained by the fact that the half corrugation's neutral axis actually has an inclination to the 𝑥𝑦-plane. The boundary conditions applied, however, reduce the effect of twisting of the cross-section.

The orthotropic model's buckling stress is very close to that of the analytical. The corrugated model, however, has a significantly lower first buckling stress. This may be explained by two factors: the bending rigidities used in the orthotropic model showed slightly stiffer characteristics in the free vibration comparison compared to the corrugated model for the two natural modes that best fit the buckling mode, that is, mode type b and d in Figure 5, cf. Table 5, results under Samanta compared to Original. Furthermore the corrugated model cross-section shows the same effect as that of the half corrugation model, and may twist slightly to reduce the bending rigidity.

The local buckling stresses are better matched, especially for the half corrugation model and the analytical beam reference.

The buckling analysis also showed that a safety margin towards buckling of the corrugated sheets has been used in this specific construction. The first buckling mode occurred at a stress of 1288 MPa, well above the yield stress of the most common steels used for this type of construction.

8. Conclusions

Depending on application, the orthotropic model may be a good substitution to reduce number of elements needed in an FE model. Computational time can be reduced with an orthotropic model since the number of elements needed can be lowered below the smallest amount needed to fully model each corrugation (the entire plate may be modeled as one element; however, this would put somewhat extreme restrictions on the type of calculations, and boundary conditions that can be performed and used with accurate results).

If accurate buckling stresses are sought after, a better orthotropic model may be able to approximate the plate thickness of the orthotropic model as 𝑡(𝑙/𝑐) instead of approximating the density as 𝜌(𝑙/𝑐). However, this has implications on all bending and extensional moduli presented herein since these are derived using 𝑡 as the orthotropic sheet thickness. This would affect the calculations made in (6) through (7). In this study, the thickness 𝑡 of the orthotropic sheet was set to 1.5 mm in these equations.

A mix of bending rigidities from the different models could be used, that is, 𝐷𝑥 and 𝐷𝑦 from Samanta and Mukhopadhyay [1] and 𝐷𝑥𝑦 from Briassoulis [2] and for the extensional rigidities the expression given by Briassoulis alone.

List of Symbols

𝑡:Geometric parameter of corrugation, see Figure 1
𝑐:Geometric parameter of corrugation, see Figure 1
𝑙:Geometric parameter of corrugation, see Figure 1
𝑑:Geometric parameter of corrugation, see Figure 1
𝑓:Geometric parameter of corrugation, see Figure 1
𝜃:Geometric parameter of corrugation, see Figure 1
𝑏:Geometric parameter of corrugation, see Figure 1
𝑥,𝑦,𝑧:Directional coordinates, cf. Figure 1
𝜎:Stress component
𝜖:Strain component
𝜈:Poisson's ratio
𝐸:Young's modulus for bending
𝐺:Shear modulus for bending
𝐷:Bending rigidity
𝐸𝑒:Young's modulus for extension and compression
𝐺𝑒:Shear modulus for extension and compression
𝐵:Extensional rigidity
𝛼:Offset contribution to the moment of inertia of half a corrugation
𝐼:Moment of inertia
𝐴:Cross-section area
𝑃cr:Buckling load
𝑢,𝑣,𝑤:Displacements in 𝑥, 𝑦, and 𝑧, respectively.


This work is part of the project: “A Light Weight Car body for High-Speed Trains”, a Ph.D. project within “Multi-functional body-panels” under the Centre for Eco2 Vehicle Design at the Royal Institute of Technology in Stockholm, Sweden, with funding from Vinnova, Bombardier Transportation, SAAB automobile, and A2Zound.


  1. A. Samanta and M. Mukhopadhyay, “Finite element static and dynamic analyses of folded plates,” Engineering Structures, vol. 21, no. 3, pp. 277–287, 1999. View at Google Scholar
  2. D. Briassoulis, “Equivalent orthotropic properties of corrugated sheets,” Computers and Structures, vol. 23, no. 2, pp. 129–138, 1986. View at Google Scholar
  3. K. Liew, L. Peng, and S. Kitipornchai, “Buckling analysis of corrugated plates using a mesh-free Galerkin method based on the first-order shear deformation theory,” Computational Mechanics, vol. 38, no. 1, pp. 61–75, 2006. View at Publisher · View at Google Scholar
  4. HyperWorks 10, Online Help and Documentation,
  5. T. J. Easley, “Buckling formulas for corrugated metal shear diaphragms,” Journal of the Structural Division, vol. 101, pp. 1403–1417, 1975. View at Google Scholar
  6. T. J. Easley and E. D. McFarland, “Buckling of light gage corrugated metal shear diaphragms,” Journal of the Structural Division, vol. 95, no. 7, pp. 1497–1516, 1969. View at Google Scholar
  7. E. D. McFarland, An investigation of the static stability of corrugated rectangular loaded in pure shear, Ph.D. thesis, University of Kansas, Lawrence, Kan, USA, 1967.
  8. T. H. G. Megson, Aircraft Structures for Engineering Students, Butterworth-Heinemann, Boston, Mass, USA, 4th edition, 2007.