A linear elastic model for the in-plane behavior of a honeycomb structure
G.J. Oosterbaan (ID: 0543487)
MT 08.28
Bachelor Thesis Report
Supervisor: dr.ir. J.A.W. van Dommelen
Eindhoven University of Technology Department of Mechanical Engineering Section Mechanics of Materials
Eindhoven (Netherlands), July 2008 Abstract
The in-plane behavior of regular hexagonal honeycomb structures will be discussed in this report. In chapter 2 a homogenization procedure is used to determine the plane stress com- pliance matrix which gives a linear relation between strain and stress. For the determination of the compliance matrix basic beam bending equations are used. These basic beam bending equations use the assumptions of Bernoulli/Euler/Kirchhoff. The compliance matrix is only valid in the linear elastic regime. The linear elastic regime is given by the initial yield surface shown in section 3.1. Initial yielding will occur if the first material point reaches the yield stress. If the total cross section of a cell wall reaches the yield stress plastic collapse will occur. The plastic collapse yield surface is given in section 3.2. The initial yield surface has the shape of a parallelogram the plastic collapse yield surface has an ellipse shape. These shapes are due the bending dominated character for uniaxial loading and the tensile character for equi-biaxial loading. For validation purposes and visualization of the Von-Misses stress distribution there are done some numerical simulations in chapter 4. From these numerical experiments it can be concluded that for high slenderness ratios the compliance matrix is accurate. Furthermore it can be seen that the locations with the highest Von-Misses stresses are near the cell wall intersections.
Preface/Acknowledgements
The demand for a material model of a hexagonal honeycomb structure is initialized by the University Racing Eindhoven Team (URE). The URE uses honeycomb sandwich panels as a body for their race car. Because the URE has switched to a carbon fiber body in 2008 there was no demand from their side anymore. Because there was no more demand from the URE this report has a more scientific and theoretical approach.
Furthermore I would like to thank my coaches Varvara Kouznetsova and Hans van Dom- melen from the Mechanics of Materials Department of the TU/e for their inspiration and patience with this research.
Contents
1 Introduction 7
2 Linear elastic model for in-plane loading 9 2.1 Homogenization procedure ...... 9 2.1.1 Constitutive relations ...... 9 2.1.2 Applying constitutive relations to unit cell ...... 11
3 Yielding under biaxial loading 15 3.1 Initial Yielding ...... 16 3.2 Full plasticity collapse ...... 17
4 Numerical simulations 21 4.1 Modeling in Marc-Mentat ...... 21 4.2 Results ...... 27 4.2.1 Subproblem 1: loading in 11-direction ...... 27 4.2.2 Subproblem 2: loading in 22-direction ...... 28 4.2.3 Subproblem 3: shear loading ...... 28 4.3 Von Mises Stress distribution ...... 30
5 Conclusion 35
A Constitutive relations 37 A.1 Derivation of energy density ...... 37 A.2 Derivation of strain tensor ...... 38 A.3 Derivation of stress tensor ...... 38
B Point symmetry and displacement relations 41 B.1 Subproblem 1 ...... 41 B.2 Subproblem 2 ...... 41 B.3 Subproblem 3 ...... 42
C Displacements as a function of the stress σ 45 C.1 Subproblem 1 ...... 46 C.2 Subproblem 2 ...... 47 C.3 Subproblem 3 ...... 49
3 4 List of Symbols
Scalars symbol description dimension c distance neutral bending line of cell wall [m] d size of a cell [m] fi force in ii-direction [N] h cell wall height [m] I second moment of area [m4] L size RVE [m] l cell wall length [m] M bending moment [Nm] Mp moment of full plasticity [Nm] ds(k) unit cell surface area belonging to point (k) [m2] t cell wall thickness [m] (k) ui displacement point (k) in ii-direction [m] W energy density [J/m2] V unit cell volume [m3] εij macroscopic strain element in ij-direction [-] ν Poisson’s ratio honeycomb [-] νs Poisson’s ratio solid material [-] σij macroscopic stress element in ij-direction [Pa] σa axial stress in cell wall [Pa] σb,max maximum bending stress in cell wall [Pa] σys yield stress material [Pa] ∗ σiy uniaxial initial macroscopic yield stress [Pa] ∗ σpl uniaxial full plasticity macroscopic yield stress [Pa]
Columns symbol description dimension ~e orthonormal base vector column [-] f˜ force column [N] u˜ displacement column [m] ˜ε strain column [-] σ˜ stress column [Pa] ˜ 5 Vectors symbol description dimension f~(k) force vector of point (k) [N] ~n(k) normal vector of point (k) [m] ~p(k) stress vector of point (k) [Pa] ~u(k) displacement vector of point (k) [m] ~x(k) coordinate vector from the origin O to point (k) [m] ~0 zero vector [-]
Tensors symbol description dimension ε strain tensor [-] σ stress tensor [Pa]
Matrices symbol description dimension C elasticity matrix [Pa] R rotation matrix [-] S compliance matrix [1/Pa]
Mathematical Operators symbol description Pm k=1 sum over the elements k=1,2,3,...,m · inner product × crossproduct : double inner product aT transpose of a ˜ ˜
6 Chapter 1
Introduction
In this report a linear elastic in-plane honeycomb material model will be given. The model will be for a regular hexagonal honeycomb structure. In figure 1.1.a a beehive honeycomb is shown and in figure 1.1.b a man made aluminium honeycomb is shown. Because nature has a tendency for mechanically efficient materials the honeycomb structure can be found in a lot of natural cellular materials as can be seen in figure 1.2. The man made honeycomb structures are commonly used as the core material of honeycomb sandwich panels or as kinetic energy absorbers. In structural applications the material should remain in the elastic regime, whereas for energy absorbtion plastic deformation is needed. Honeycomb structures are not just used in sandwich panels or as energy absorbers but have many other applications such as air directionalization, thermal panels, acoustic panels, light diffusion and radio frequency shielding. The main advantage of honeycomb panels is its high stiffness and bending strength combined with a low weight. The in-plane behavior of Honeycomb structures has to be studied to know the contribution of the honeycomb core to the total panel strength. Also cellular solids are often approximated by this linear elastic honeycomb model. This report will focus on the linear elastic behavior for structural applications. The main parameters of interest are:
(1) a relationship between the applied stress on the honeycomb and the resulting macroscopic strains; (2) the maximal stress that can be applied on the honeycomb before yielding occurs; (3) the maximal stress that can be applied on the honeycomb before total collapse occurs.
Chapter 2 will be devoted to (1). A homogenization technique is used to derive the compli- ance matrix S which gives a relation between the macroscopic strain and macroscopic stress of the honeycomb continuum. Also the macroscopic Poisson’s ratio will follow from this ho- mogenization procedure. In chapter 3 the maximal stress that can be applied before yielding occurs and before the plastic collapse of the honeycomb structure occurs are considered. This results in two yield surfaces: one for initial yielding and one for plastic collapse of the honey- comb. In chapter 4 some numerical analyses are done with Marc-Mentat for validation of the analytical model and visualization of the deformations and the Von Mises stress distribution. In the last chapter a conclusion will be drawn and some recommendations for further research will be given.
7 (a) Beehive honeycomb. (b) Aluminium hexagonal honeycomb.
Figure 1.1: Regular honeycomb structures.
Figure 1.2: Natural cellular materials: (a) cork (b) balsa (c) sponge (d) trabecular bone (e) coral (f) cuttlefish bone (g) iris leaf (h) plant stalk.
8 Chapter 2
Linear elastic model for in-plane loading
2.1 Homogenization procedure
In this chapter the linear elastic model for the honeycomb core will be set up. The model for the honeycomb core is largely based on the models used by Onck [1] and Gibson and Ashby [2]. The honeycomb can be seen as a continuum with properties that are obtained by homogenization over a Representative Volume Element (RVE). This homogenization is valid when the micro structural cell size d is much smaller as the length scale L of the RVE and the micro structural geometry is randomly distributed or if the microstructural geometry is equal and periodically distributed. In the case of the regular honeycomb the structure is equal and periodically distributed. First of all some basic linear elastic constitutive relations will be set up. Thereafter these constitutive relations are applied to the honeycomb to derive the compliance matrix S and the elasticity matrix C. The starting relations are force and moment equilibrium and energy conservation. The honeycomb material is characterized by Young’s modulus Es, the yield stress σY , the cell wall thickness t, the cell wall height h and the length l of the cell wall. It will be assumed that a uniform stress (static boundary conditions) or a uniform strain (kinematic boundary conditions) will be applied. The procedure for static boundary conditions is visualized in figure 2.1. The honeycomb structure used is a regular hexagonal cell structure. Because this cell structure is equal and periodically distributed the smallest RVE is in size equal to a unit cell. The geometric structure and the used unit cell are shown in figure 2.2.
2.1.1 Constitutive relations For the unit cell, force equilibrium should be satisfied,
m m m X X X f~(k) = ~p(k)ds(k) = σ · ~n(k)ds(k) = ~0 (2.1) k=1 k=1 k=1 where f~(k) is the force on a cell wall (k), m is the number of cell walls in the RVE. The stress vector ~p(k) is acting on surface ds(k) and can be calculated by taking the inner product between the uniform stress tensor σ and the normal vector ~n(k) on surface ds(k).
9 r e2 r r e1 e3 s ji e ij
Figure 2.1: Homogenization with static boundary conditions, source [1].
A moment equilibrium should be satisfied, m m X X ~x(k) × f~(k) = ~x(k) × (σ · ~n(k)ds(k)) = ~0 (2.2) k=1 k=1 The vector ~x(k) is a Eulerian coordinate vector from the origin O to the boundary point of cell wall (k) on which the force f~(k) is acting. The conservation of energy principle dictates that the internal work has to be the same as the external work. Because uniform stress is assumed the energy density is equal to the work done by the applied forces averaged over the RVE, m 1 1 X 1 W = σ : ε = f~(k) · ~u(k) (2.3) 2 V 2 k=1 where the first term represents the internal stored energy W and the second term the external
10 r r (2) 60° ds f f (1) l r r n(2) n(1) (2) (1)
c 120°
3l l O
c' (3) (4) r r r n(3) n(4) e2 r r r (3) (4) e2 f f r e1 r e1 (a) Regular hexagonal cell structure. (b) Unit cell.
Figure 2.2: Unit cell structure.
work We averaged over the volume V of the RVE. In the first term ε represents the uniform strain tensor. The derivation of the energy density function can be found in appendix A.1. Inserting the relation f~(k) = σ·~n(k)ds(k) in equation (2.3) and making use of the symmetry of the stress tensor σ yields a relation for the strain tensor ε,
m 1 X 1 ε = ~u(k)~n(k) + ~n(k)~u(k) ds(k) (2.4) V 2 k=1 the derivation of the strain tensor can be found in appendix A.2. The stress tensor σ is derived by substituting the displacement vector ~u(k) = ε · ~x(k) in the energy density equation (2.3),
m 1 X 1 σ = ~x(k)f~(k) + f~(k)~x(k) (2.5) V 2 k=1
For the derivation of equation (2.5) see appendix A.3.
2.1.2 Applying constitutive relations to unit cell The center point O can be seen as a point of symmetry, see appendix B for a visualization. This means that the forces at beam (1) and beam (2) are related to beam (3) and beam (4) respectively. Thus, it suffices to take only the upper half part of the unit cell with beams (1) and (2) into account.
11 √ Applying equation (2.4) and filling in for the normal vectors ~n(1) = 3 1 ~e, ~n(2) = 2 2 √ √ ˜ − 3 1 ~e and for the surface area ds = 3l, see figure 2.2.b, leads to the following strain 2 2 ˜ tensor ε which is dependent on the displacement vectors ~u(1) and ~u(2), √ √ 2 3 u(1) − u(2) 1 u(1) + u(2) + 3 u(1) − u(2) T 3lh 1 1 3lh 1 1 2 2 ε = ~e √ ~e ˜ 1 (1) (2) (1) (2) 2 (1) (2) ˜ 3lh u1 + u1 + 3 u2 − u2 3lh u2 + u2 (2.6)
Now only a relation for the displacement vectors ~u(k) as a function of the stress tensor σ remains to be determined. This relation can be found by treating the cell walls as cantilever beams which are clamped in their intersection points and thus there are no hinges in the intersection points. The principle of superposition allows to handle the problem (total loading) as three different subproblems (load cases) and sum these afterwards. The different load cases are one for simple tension in 1-direction, one for simple tension in 2-direction, and one for shear in 12-direction. The results for the displacement vectors ~u(1) and ~u(2) as a function of the stress tensor σ are given below, for the derivation see appendix C.
Displacements subproblem 1 (tension in 11-direction)
2 4 9l σ11 3l σ11 ! + 3 (1) T √16Esht 16√Esht ~u = ~e 2 4 (2.7) 3 3l σ11 3 3l σ11 ˜ − 3 16Esht 16Esht
2 4 9l σ11 3l σ11 ! − − 3 (2) T √16Esht 16√Esht ~u = ~e 2 4 (2.8) 3 3l σ11 3 3l σ11 ˜ − 3 16Esht 16Esht
Displacements subproblem 2 (tension in 22-direction)
2 4 3l σ22 3l σ22 ! − 3 (1) T √16Esht 16√Esht ~u = ~e 2 4 (2.9) 9 3l σ22 3 3l σ22 ˜ + 3 16Esht 16Esht
2 4 3l σ22 3l σ22 ! − + 3 (2) T √16Esht 16√Ehst ~u = ~e 2 4 (2.10) 9 3l σ22 3 3l σ22 ˜ + 3 16Esht 16Esht
Displacements subproblem 3 (shear in 12-direction) √ √ 2 4 3 3l σ12 3 3l σ12 ! − 3 (1) T 8Esht 4Esht ~u = ~e 2 4 (2.11) 3l σ12 3l σ12 ˜ + 3 8Esht 2Esht
√ √ 2 4 3 3l σ12 3 3l σ12 ! − 3 (2) T 8Esht 4Esht ~u = ~e 2 4 (2.12) 3l σ12 3l σ12 ˜ − − 3 8Esht 2Esht By substitution of the expressions above in the strain tensor equation (2.6) the compliance matrix can be determined. The compliance matrix gives a relation for the strain column
12 T T ε = (ε11 ε22 ε12) as a function of the stress column σ = (σ11 σ22 σ12) . The compliance ˜matrix is given by ˜
√ l 2 l 2 3 + t 1 − t 0 3l 2 2 S = 1 − l 3 + l 0 (2.13) 4E ht t t s l 2 0 0 2 + 2 t