Roland's homepage

My random knot in the Web

Corrugation against buckling

Composite sandwich products loaded in bending tend to fail by buckling of the laminate under compression. The author’s intuition is that corrugating the surface under compression should help. So the question is; does it help, and how much. That is what will be investigated in this article.

Introduction

Target of the comparison are two sandwiches, 34 mm thick in total. Both are 1.5 m long, 48 cm wide. The laminates on both sides are 2 mm thick. The core is PMI foam. In one of the sandwiches the bottom laminate is corrugated in a sinusoidal shape with an initial wavelength of 60 mm and amplitude of 5 mm. The laminate is mostly unidirectional carbon fibers with one 200 g/m² woven layer on the outside. The sides of the sandwich are closed with 2 mm thick quasi-isotropic laminates. The cross-sections of both are shown below.

Both sandwiches are clamped at one end and loaded with a identical load (via a rigid body constraint) at the other end.

The behavior under load will be simulated with the CalculiX finite element analysis software, with the geometry being generated by the preprocessor with the help of a Python script to interpolate a sinusoid as a spline.

The meshes are shown below.

mesh-flat mesh-corr

The properties of the laminate are as follows:

Generated by lamprop version 2022.06.22
laminate: Mflat
thickness: 1.98 mm, density: 1.51 g/cm³
fiber volume fraction: 56%, fiber weight fraction: 66.6%
laminate weight: 3004 g/m², resin consumption: 1004 g/m²
num weight angle   vf fiber
    [g/m²]   [°]  [%]
1    100     0   56 T700SC
2    100    90   56 T700SC
3    300     0   56 T700SC
4    300     0   56 T700SC
5    300     0   56 T700SC
6    300     0   56 T700SC
7    300     0   56 T700SC
8    300     0   56 T700SC
In-plane engineering properties:
E_x  = 125348 MPa, E_y  = 16102 MPa, E_z  = 13097 MPa
G_xy  = 3780 MPa, G_xz  = 3166 MPa, G_yz  = 3888 MPa
ν_xy = 0.21477
α_x = 2.682e-07 K⁻¹, α_y = 1.815e-05 K⁻¹
Engineering properties derived from 3D stiffness matrix:
E_x = 126081 MPa, E_y = 19069 MPa, E_z = 13783 MPa
G_xy = 3780 MPa, G_xz = 3825 MPa, G_yz = 4639 MPa
ν_xy = 0.186, ν_xz = 0.306, ν_yz = 0.396
** Material data for CalculiX / Abaqus (SI units):
*MATERIAL,NAME=Mbot
*ELASTIC,TYPE=ORTHO
1.291e+11,6.027e+09,2.179e+10,6.048e+09,6.438e+09,1.583e+10,3.78e+09,3.825e+09,
4.639e+09,293
*DENSITY
1514

The breaking strain of the fibers used is at least 1.7%. That means that the tensile strength in the 0° direction is 2140 MPa. For these kinds of laminates, we use the rule of thumb that the compression strength is around 60% of the tensile strength, or -1280 MPa.

Orienting laminate properties

The top laminate is aligned with the global coordinate system used in the FEA. So it needs no special attention.

The planes of the side laminates are at 90° to the XY plane. So a *ORIENTATION card is needed in the solver input to align the material properties with the orientation of the laminate.

The corrugated side is more complicated. The elements on that side are curved and their normals point in different directions. Creating orientations for those elements by hand, while possible, would be cumbersome and error-prone. So the author wrote a program for that: auto-orient.

This analyzes the sets of elements it is invoked on, and generates suitable orientations, element sets and *SOLID SECTION cards. These can simply be included in the solver job file.

Corrugated version 1 vs flat version 1

disp-flat disp-corr

The sandwich with the corrugated bottom deflects more (-115 mm) than the flat-bottomed sandwich (-99 mm). That is to be expected; its cross-section has a lower second area moment because most of the bottom laminate is closer to the top laminate than the flat version.

stress-flat stress-corr

Results for flat1:

displacement -98.95 mm
compression stress -406 MPa
compression stress factor 3.17  
tensile stress 406 MPa
tensile stress factor 5.28  
buckling factor 2.26  

Results for corrugated1:

displacement -114.79 mm
compression stress -434 MPa
compression stress factor 2.96  
tensile stress 429 MPa
tensile stress factor 5.00  
buckling factor 2.50  

The maximum stresses are higher in the sandwich with the corrugated bottom than those in the regular sandwich. This is also due to the lower stiffness. The ratio between the compression strength and the occurring stress (the compression factor) is therefore higher for the flat panel than for the corrugated panel.

The buckling factor is slightly higher then for the flat panel. But in both cases the buckling factor is smaller than the compression stress factor. So both will fail in buckling.

Corrugated version 2: updated bottom laminate

In the bottom laminate one layer of unidirectional fibers was replaced by a layer of ±45° fiber of equivalent weight.

This has changed the laminate properties from those shown in above to:

Generated by lamprop version 2022.06.22
laminate: Mbot
thickness: 1.98 mm, density: 1.51 g/cm³
fiber volume fraction: 56%, fiber weight fraction: 66.6%
laminate weight: 3004 g/m², resin consumption: 1004 g/m²
num weight angle   vf fiber
    [g/m²]   [°]  [%]
1    100     0   56 T700SC
2    100    90   56 T700SC
3    300     0   56 T700SC
4    300     0   56 T700SC
5    150    45   56 T700SC
6    150   -45   56 T700SC
7    300     0   56 T700SC
8    300     0   56 T700SC
9    300     0   56 T700SC
In-plane engineering properties:
E_x  = 109107 MPa, E_y  = 19647 MPa, E_z  = 13097 MPa
G_xy  = 8358 MPa, G_xz  = 3250 MPa, G_yz  = 3804 MPa
ν_xy = 0.39130
α_x = 4.135e-08 K⁻¹, α_y = 1.377e-05 K⁻¹
Engineering properties derived from 3D stiffness matrix:
E_x = 110279 MPa, E_y = 22902 MPa, E_z = 14148 MPa
G_xy = 8418 MPa, G_xz = 3893 MPa, G_yz = 4571 MPa
ν_xy = 0.349, ν_xz = 0.243, ν_yz = 0.377
** Material data for CalculiX / Abaqus (SI units):
*MATERIAL,NAME=Mcorr
*ELASTIC,TYPE=ORTHO
1.155e+11,1.066e+10,2.609e+10,6.081e+09,6.405e+09,1.583e+10,8.418e+09,3.893e+09,
4.571e+09,293
*DENSITY
1514

The Young’s modulus in the length direction has been reduced significantly. But the shear modulus has been more than doubled.

Displacement and stress images:

disp-corr2 stress-corr2

Displacement and stress values:

displacement -121.54 mm
compression stress -404 MPa
compression stress factor 2.78  
tensile stress 429 MPa
tensile stress factor 5.00  
buckling factor 2.87  

Since the longitudinal stiffness of the bottom laminate has been reduced, it is not suprising that the deflection has been increased.

The buckling factor has been increased from 2.5 to 2.87 with just the one small change to the laminate.

The compression stress factor is now slightly lower than the buckling factor. So this sandwich will likely fail in pure compression rather than buckling.

This bottom laminate will be used for the corrugated laminates from now on.

Flat version 2: changed bottom laminate

To discern how much improvement in the previous alternative was due to the material properties, a sandwich with a flat bottom laminate with a layer of ±45° fiber was also calculated.

Results overview:

displacement -104.54 mm
compression stress -385 MPa
compression stress factor 2.93  
tensile stress 405 MPa
tensile stress factor 5.29  
buckling factor 2.12  

The displacement has increased as expected because of the lower longitudinal stiffness of the bottom laminate.

It’s buckling factor is only 2.12. So the corrugation and ±45° layer really have a strong synergetic effect. But the ±45° layer is not really useful for the flat bottom laminate.

Corrugated version 3: increased amplitude

For the next test, the amplitude of the corrugation is changed from 5 mm to 10 mm. The position of the top laminate is adjusted so the total thickness does not exceed 34 mm. The sides of the sandwich are closed with a quasi-isotropic laminate, as before. The cross-section is shown below.

Deflection and stress are shown below.

disp-corr3 stress-corr3

Deflection and stress values:

displacement -140.36 mm
compression stress -482 MPa
compression stress factor 2.33  
tensile stress 441 MPa
tensile stress factor 4.85  
buckling factor 4.51  

The deflection is again higher because more of the bottom laminate is closer to the top laminate, leading to a reduced stiffness.

The buckling strength of this sandwich is twice that of the original sandwich, with the same overall thickness.

However in this case the compression stress limit would be reached with a load factor of 2.3, so long before buckling occurs.

Corrugated version 4: thicker sandwich with 5 mm amplitude

In this variant, the total thickness is increased to 36.5 mm so that the middle of the bottom laminate is 34 mm from the top. The amplitude of the corrugation is the original 5 mm.

Displacement and stress images:

disp-corr4 stress-corr4

Displacement and stress values:

displacement -103.92 mm
compression stress -378 MPa
compression stress factor 2.98  
tensile stress 403 MPa
tensile stress factor 5.32  
buckling factor 3.14  

Deflection and stresses are close to the “flat2” variant, which has the same laminate build-up.

The buckling factor is slightly larger than the compression stress factor. So in this case the laminate would likely also fail in compression just before it buckles.

Corrugated version 5: corrugation with 120 mm wavelength

This is a variant on version 2, where the wavelength of the corrugation has been changed from 60 mm to 120 mm.

The amplitude of the corrugation and the total height have not been changed compared to version 1 and 2.

The maximum deflection is -122 mm and the maximum stresses are +424 MPa and -400 MPa. This is similar to corrugated version 2.

However, the buckling factor is only 2.04, which is worse than the flat version 1! So the wavelength of the corrugation has a huge influence, and shorter is better.

So the amplitude will be increased to 13 mm and the total thickness to 40.5 mm.

Displacement and stress images:

disp-corr5 stress-corr5

Displacement and stress results:

displacement -100.26 mm
compression stress -410 MPa
compression stress factor 2.74  
tensile stress 382 MPa
tensile stress factor 5.62  
buckling factor 2.97  

Since buckling factor is higher than the compression stress factor, this sandwich will likely fail in compression. So the larger the wavelength, the larger the amplitude needs to be to generate the desired improvement in buckling strength.

Corrugated version 6: corrugation with 30 mm wavelength and amplitude 3 mm

The question here is what happens is we reduce the wavelength and the amplitude.

The cross-section of this version is shown above. From a standpoint of manufacturability this is not ideal. Getting the fibers from the bottom layers to not bridge over the corrugations will probably require tooling. E.g. heavy metal rods placed in the corrugations in the middle and rolled outwards from valley to valley when putting down each layer.

Displacement and stress images:

disp-corr6 stress-corr6

Displacement and stress results:

displacement -114.93 mm
compression stress -399 MPa
compression stress factor 2.82  
tensile stress 431 MPa
tensile stress factor 4.97  
buckling factor 2.95  

Here as well the laminate will probably fail in compression before buckling occurs.

With a smaller wavelength, a smaller amplitude is sufficient to reach a high enough buckling factor.

Conclusion

A corrugated laminate on the compression side of a sandwich can significantly improve the buckling strength. It is possible to make the calculated buckling strength at least equal to the pure compression strength of the laminate.

A sandwich with the corrugated bottom is less stiff than the sandwich with a flat bottom of the same total thickness. This is to be expected:

  • Most of the bottom laminate is closer to the top laminate, which means a reduced second area moment and so a reduced geometric stiffness.
  • The corrugated side needs some ±45° fibers to increase its shear stiffness. This decreases the longitudinal material stiffness.

For the sandwich in question (34 mm thick, 2 mm skins), several different corrugations were proven suitable in this way;

  • 60 mm wavelength with an amplitude of 5 mm.
  • 30 mm wavelength with an amplitude of 3 mm.

There is clearly a correlation between the wavelength and amplitude; larger wavelengths require larger amplitudes. Approximately, the amplitude should be 10% of the wavelength.

From a standpoint of manufacturing, longer wavelengths are going to be easier than shorter ones. However, from the standpoint of x-ray image quality, lower amplitudes are probably going to be preferred.


For comments, please send me an e-mail.


Related articles


←  Approximating elliptical arcs in CalculiX Graphics Meshing a circle with hex elements in CalculiX  →