`ISRN Mechanical EngineeringVolume 2011 (2011), Article ID 979532, 9 pageshttp://dx.doi.org/10.5402/2011/979532`
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.

Abstract

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, 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: For an arbitrary orthotropic, thin plate, the flexural rigidities may be expressed as 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] Rearranging (4) gives Utilising (3) to express the Young's modulus as a function of the bending rigidity, we get

The extensional rigidities and moduli of an orthotropic plate (here denoted and to distinguish from the flexural expressions) follow the relations [2]

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: 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]

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 Briassoulis' equations are modified equations from, for example, Easley [5].

A third variant of the flexural rigidities is given by Samanta and Mukhopadhyay [1] 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: while Briassoulis has derived different versions of and (see below) 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.

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 by  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] 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 The critical stress is then calculated as

where is the cross-section area of the half-corrugation.

The analytical local buckling stress is given by [8] 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 This gives a critical stress of

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)

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

This gives us a critical stress ofThe moment of inertia for the half corrugation beam is

Furthermore,  GPa, and  mm. This gives an analytical Euler buckling load, according to (18), of

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

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: 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 𝜌 : Density 𝐼 : Moment of inertia 𝐴 : Cross-section area 𝑃 c r : Buckling load 𝑢 , 𝑣 , 𝑤 : Displacements in 𝑥 , 𝑦 , and 𝑧 , respectively.

Acknowledgment

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.

References

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.
2. D. Briassoulis, “Equivalent orthotropic properties of corrugated sheets,” Computers and Structures, vol. 23, no. 2, pp. 129–138, 1986.
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.