Issue 
Mechanics & Industry
Volume 24, 2023
High fidelity models for control and optimization



Article Number  43  
Number of page(s)  7  
DOI  https://doi.org/10.1051/meca/2023040  
Published online  18 December 2023 
Research Article
On the response of reentrant auxetic structure under dynamic crushing
^{1}
LaMé, Univ. Orléans, Univ. Tours,
INSA CVL 63 av de Lattre de Tassigny,
18020
Bourges, France
^{2}
CEDREM, ECOPARC Domaine de Villemorant
41210
NeungsurBeuvron, France
^{3}
LAMIH UMR CNRS 8201, Université Polytechnique HautsdeFrance Campus MontHouy,
59313
Valenciennes Cedex 9, France
^{*} email: adeline.petit@etu.univorleans.fr
Received:
14
February
2023
Accepted:
13
November
2023
The present work focuses on the dynamic crushing response of 2D reentrant auxetic honeycomb. The aim of this work is to propose analytic formulations to predict the crushing strength of the considered structure by including design parameters (specific geometrical ratios and/or material properties). The model for dynamic impact is based on shock waves propagation analogy in a rigid, perfectly plastic, locking (RPPL) material model. The formulations depend on the geometric and the material characteristics of the auxetic but also on the impact velocity. Finite Element simulations of the impact of a plate on an auxetic structure were carried out using the RADIOSS^{®} explicit solver. The impact velocities varies from 35 m/s to 200 m/s. The numerical simulations presented in this study show good accordance between analytical and Finite Element results.
Key words: Auxetic structure / impact / finite element model / explicit simulation / analytical solution
© A. Petit et al., Published by EDP Sciences, 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Structured materials, known for their multifunctional application [1], are remarkable for their particular stressstrain curves: after an elastic part where the curve is quasilinear, a plateau stress interval is reached where the stress is almost constant and finally the curves rise abruptly at densification strain [2]. During the past few decades [3,4], a new class of structured materials have emerged called ≪ auxetic ≫ (from Greek ≪ auxetos ≫: that may be increased). This kind of structure is characterized by a negative Poisson’s ratio which gives it, compared to conventional, cellular structure a higher shear modulus, better indentation resistance and energy absorption ([5–7]).
The present work focuses on one of these improved properties: the energy absorption capacity of auxetic structures (AS). Among the large range of structures having a negative Poisson’s ratio, the reentrant honeycomb is one of the most investigated geometry. Experimental samples are relatively easy to develop using 3D printing technologies [8]. In addition, this geometry can be implemented in FE models without significant difficulties. That is why it has been chosen for this study.
Inplane auxetic crushing test is the benchmark analysis to determine quasistatic and dynamic crushing strength of AS. This strategy is adopted by Ruan et al. [9], Hu et al. [10,11], and Hou et al. [12] to develop their models. In 1993, Reid et al. [13], on the basis that wood behave like cellular structure, proposed a shock model for wood using the rigid, perfectlyplastic, locking (RPPL) simplification. This model of the dynamic crushing of cellular structure is validated in 2003 by Ruan et al. [9] thanks to numerical inplane dynamic crushing of honeycombs.
Based on energy conservation principle and shockwave analogy, a theoretical (analytical) formulation of stress will be established. Finite Element simulations have confirmed the analytical model.
The objective is to propose a general methodology to assess the strength of an AS loaded by the impact of a plate.
2 Characteristics of the auxetic structure under consideration
2.1 Geometry of the structure
The complete auxetic structure, which will be investigated in compression along the Yaxis in this article, is an auxetic “reentrant honeycomb” represented in Figure 1.
The structure is composed of cells, as the one shown in Figure 2 with the notations used for the geometrical data. Note that S_{0} is the total projected area on the X–Z plane, as shown in Figure 1.
A unitcell geometry is characterised by its neutral axis (black dotted line in Fig. 2). The horizontal edge is denoted H and the bevel edge is denoted L. The angle formed by edges H and L is denoted θ. The cellwall thickness is denoted t_{thick} and b is the outofplane thickness.
The edges dimensions h and ℓ are defined as follow :
The initial cell height is defined by:
At the initial time (t = 0), the total height in the Ydirection for n cell rows (layer) is:
Two characteristic ratios are defined:
Fig. 1 View of the complete auxetic structure (“AS”). 
Fig. 2 Reentrant auxetic cell configuration. 
2.2 Relative density and densification strain
The relative density of a cellular structure is defined as the ratio of the cellular structure mass density ρ_{c} to the mass density ρ_{s} of the material of which it is constituted. For the reentrant honeycomb, it is defined as follow :
where V_{c} and V_{s} are the cellular structure and material volumes, respectively. The relative density can be written in terms of the nondimensional ratio τ and η :
The strain will be defined by the ratio of the height variation D_{0} − D to the initial height D_{0}, D being the height at time t:
At densification the strain is ε_{DS} and the height is D_{s}, while the travelled distance of the upper wall is ε_{DS}D_{0}.
At this characteristic strain, the cells are all collapsed so that all cellwalls are in contact and the compressive stress is almost the compressive strength of the material itself. So, the height is D_{s} = (4 n + 1) t_{thick} and the densification strain ε_{DS} is therefore expressed as:
3 Analytic prediction of stress
3.1 Stress enhancement through the shock front
Now, an impact on the AS is considered. It corresponds to impact the top surface (y = D_{0}, Fig. 1) with a plate of mass M, while the other face (y = 0) is fixed. The velocity of the plate is υ_{0} when it impacts the AS at time t_{0} = 0, and υ for t > 0.
When the AS is dynamically crushed, the compaction starts from the impacted face and propagates like a shockwave through the structure, involving a strong discontinuity (or “front”) between the undeformed and deformed regions. The main steps of the shockwave analogy applied to cellular material have been presented by [13]. As a similar approach is adopted in this study, these steps are briefly recalled in Figure 3.
The AS is idealized with the RPPL model. Therefore, the stress is raised instantaneously to the static plateau stress σ_{stat} at the initial time. The σ_{stat} expression can be obtained by considering the plastic collapse of the cells network [14] :
with σ_{y} the yield stress of the material.
Figures 3c and d show the deformed states at times t and t + dt. Behind the front, the stress is raised to the dynamic value σ and the structure is compressed up to densification strain ε_{DS}. The structure is supposed to keep the same area S_{0} and to deform only in the uniaxial direction.
In Figure 3, the shock wave propagates from proximal to distal end with the wave speed υ_{f}. From kinematic analysis, this velocity in lagrangian coordinates is expressed as:
where u is the displacement of the plate and λ is the deformed length of the crushed AS. As the current velocity is υ = du/dt, the front velocity relative to the impactor is :
The conservation of mass gives :
where ρ_{DS} is the AS density when the structure is crushed up to densification.
During an infinetely small duration dt, the element of length dD (hashed in Fig. 3c) passes from an undeformed to a crushed state (hashed in Fig. 3d) and its lengths is dλ. The mass included within this volume hashed in Figure 3d, S_{0} dλ, is expressed in terms of the cellular density ρ_{c} :
During this time interval, the volume element S_{0} dλ considered above is submitted to the stress σ_{stat} on its right face and to the stress σ on its left face. The velocity of the element increases from 0 to υ during this duration. Thus the conservation of momentum of this element reads:
From equation (16) it can be deduced that:
This formulation is the same as the one proposed by Reid [13] to assess the stress inside the crushed part (shaded part in Fig. 3). After the impact, the crushing velocity υ and the reaction forces of the AS vary constantly with time due to the dynamic coupling between the plate and the AS at the interface. It is proposed to evaluate this evolution with respect to time by using the energy conservation approach.
Fig. 3 Schematics and parameters defining a shock propagation into cellular bar after a rigid mass impact. 
3.2 Analytical determination of displacement, velocity and stress
The conservation of energy of the whole system during the period t to t + dt (Figs. 3c and 3d) reads:
The left hand side of equation (18) includes a part of the structure which is set in motion by the impactor. This part shares a fraction of the system kinetic energy with the impactor. A fraction of the initial kinetic energy is transferred to the structure in the form of internal (strain) energy.
Equation (18) is expanded with the help of equations (13), (12), (14) for dλ, dD and ρ_{DS} respectively. In addition the second order term ρ_{c} S_{0} du dυ/ɛ_{DS} is neglected in the left hand side of equation (18). Finally, the nonlinear differential equation describing the plate decelleration is:
Where the constants B and C are defined by:
Equation (19) can be written in term of the displacement u only by using the second and first time derivatives of u:
Note that equation (22) holds only if u, du/ dt are all both non constant functions of time.
Solving equation (22) with the initial conditions du/ dt = υ_{0} and u = 0 at t = 0 gives :
The time derivative of of u(t) gives the impactor velocity:
Substituting equation (24) into equation (17) gives the formulation of dynamic stress with respect to time and impact conditions.
When u(t) and υ(t) are known the stress σ is directly calculated by equation (17). In addition, the following physical quantities of interest can also be determined :
– Stop time for the rigid plate. According to equation (24), if the time t equals
the velocity of the plate is theoretically zero.
– Time for densification. At the particular time value t_{DS} the AS begins to densified and the corresponding strain is ε_{DS} given by equation (10). By equating the crushed length to ε_{DS}D_{0}, the time for densification is the the appropriate (physical) time root t_{DS} of the following equation:
Data for the AS and plates in the FE model (lengths are in millimetre).
Dimensions of the plates.
Material properties for the FE model.
4 Finite element model
A finite element model was developed to simulate the impact of the plate on the AS represented in Figure 1 and to evaluate the validity of the previous formulations (Sect. 3.2). The commercial explicit nonlinear FE software RADIOSS^{®} (Altair HyperWorks) [15] is used to implement this model. The model is composed of a top falling plate, an auxetic structure supported by a fixed plate at y = 0. The auxetic is made of six cells in height Ydirection and five cells in horizontal Xdirection. The characteristics dimensions for the AS and the plates are given in Tables 1 and 2.
An elasticperfectly plastic aluminium is used for the auxetic material. A perfectlyelastic material is used for the two plates. The material properties are given in Tables 3. In the simulation, the crushing plate is supposed elastic with an elastic Young modulus equal to that of steel. The mass density is adjusted so that the crushing plate (with its imposed dimensions) has a mass equal to 4 g. This gives the plate a density of ρs = 1.11 × 10^{4} kg/m^{3}.
In this simulation, Poisson’s ratio is not a critical parameter, as the behaviour of the plates is assumed to be elastic. As its value can vary between 0.28 and 0.34, we have chosen 0.33 to represent an average value for the steel. Poisson’s ratio also has no influence on the perfectly plastic elastic behaviour of the auxetic structure material. Thus, this material property is not a critical parameter in the simulation. For the sake of simplicity, we have set the same value of 0.33 for plates and aluminium.
The mesh size sensitivity has been analyzed with a coarse mesh and four gradually refined meshes (from 2 to 5 elements in the thickness, corresponding to sizes from t_{thick}/2 to t_{thick}/5). Variations of the impactor velocity results υ for these meshes are presented in Figure 4. It is found that the order of the deviation is of 1% between the one and five elements in thickness meshes. If the element sizes are equal or lower than t_{thick}/2, the results differ of less than 0.5%. Therefore, a mesh with two elements in the thickness (0.15 mm size) was retained.
The mesh of the three parts is composed of underintegrated 8node hexaedral solid elements with physical hourglass stabilization. The plates are meshed with elements of 1 mm length. The auxetic structure is meshed with three elements in Zdirection. The fixed plate is constraint in all directions, the top plate is able to move only in Ytranslation direction and, to avoid outofplane displacement, the auxetic is fixed in Ztranslation direction. A penalty method general contact is applied to model all contact interactions at the plate/auxetic and auxetic/auxetic interfaces. The contact friction is not taken into account.
Fig. 4 Finite element impactor velocity results υ for several mesh sizes (from one to five elements in the thickness t_{thick}) at initial velocity υ_{0} = 200 ms^{−1}. 
5 Results and dicussion
We consider the crushing of the AS by a falling plate of initial velocity υ_{0} in Ydirection. The simulations are done for a 4 g mass plate with initial velocity υ_{0}: 35, 50 and 200 ms^{−1}. The dynamic response of the AS obtained through the FE simulation is expressed in terms of stress σ (left scale of the graph) and the instantaneous plate velocity υ (right scale of the graph) υs. time are plotted in Figures 5b, 6b and 7b, respectively. Figures 5a, 6a and 7a give respectively the AS deformation at selected times.
The plate is supposed always in contact with the AS. The stress σ_{FE} is calculated by the ratio of the contact force (crushing plate – AS) to the initial transverse (XZ plane) area S_{0}. The dynamic stress σ is calculated with equation (17) by using the instantaneous velocity equation (24). The constant plateau stress σ_{stat}, equation (11), is also plotted. The times t_{stop} and t_{DS} are calculated with equation (25) and equation (26), respectively.
The calculated stress, equation (17), leads to accurate estimation of the AS average stress values for the different impact conditions.
In Figures 5b, 6b and 7b an initial peak stress can be observed for a short duration which causes an important absorption of energy. Hence, during this same time lapse, the FE velocity decreases sharply. This initial peak is not taken into account in the theory developed in Section 3.2. That is why FE and analytic velocity differ from each other at the beginning.
At 35 ms^{−1}, the FE simulation shows that the plate is completely halted 1.85 ms after impact. The time t_{stop}, calculated with equation (25), is 2.16ms, which indicates that the method presented in Section 3.2 overestimates the stop time for low initial velocity. After the period of initial peak stress, values of plate FE velocity and analytical velocity, determined with equation (24), are very close. However, once the velocity is lower than about 20 m s^{−1} , the two values come apart from each other.
It can be seen, in Figure 5a, that the deformed region differs from that of the shockwave model used to established the formulation of σ. In fact, the deformed region observed in the FE results corresponds to deformation modes, which depend upon υ_{0} and the AS topology and geometry, (see for example Ruan et al. [9] for constant crushing velocity).
At 50 m s^{−1} and 200 m s^{−1} , the AS was not able to stop the plate before densification. When densification begins the plate velocity is not zero. The plate will continue to crush the structure while the stress increases up to the peak and rebound of the plate will occurs. The prediction of t_{DS} with equation (26), is reasonably accurate regarding FE results. For both, the evaluation of analytic plate velocity are very close to the FE velocity.
The formulations established in Section 3.2 seem to be more accurate regarding FE results for medium ( 50 m s^{−1} ) and high initial impact velocity (200 ms^{−1}).
The alloy properties chosen in this article are close to those of the AL3024T351 alloy. Khan et al. [16] have shown that if the strain rate does not exceed 1500 s^{−1} , the behaviour of the material is not very sensitive to strain rate. As we have verified that in the simulations of the dynamic compression of the structure, the dε/dt rate remains below 100 s^{−1}, we have chosen not to take the strain rate into account in the behaviour of the material making up the AS.
6 Conclusion
In this paper, closedform formulas have been proposed for predicting the crushing strength of architectured auxetic lattices loaded by the impact of a plate. The approach is based on energy conservation principle and shockwaves propagation analogy in a rigid, perfectly plastic, locking material model. It enables the analysis of periodic collapse of the structure. We have extended formulation in order to include more design parameters (specific geometrical ratios and/or material properties). Our analysis enables to theoretically predict the dynamic crushing strength. The formulation depends on the geometric and the material characteristics of the auxetic but also on the impact velocity.
The results of analytical calculations are consistent with the results of finiteelement simulations. It is then possible to adapt the geometrical and material parameters in order to optimize the energy absorption ability for facing specific load cases. Indepth studies are in progress to understand the influence of the deformation modes on energy absorption.
Fig. 5 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 35 ms^{−1}. (a) AS deformation at t = 0.77 ms. (b) Stress and plate velocity υs. time. 
Fig. 6 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 50 ms^{−1}. (a) AS deformation at t = 0.45 ms. (b) Stress and plate velocity υs. time. 
Fig. 7 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 200 ms^{−1}. (a) AS deformation at t = 0.17 ms. (b) Stress and plate velocity υs. time. 
Acknowledgement
The authors acknowledge the French Region CentreVal de Loire for having provided the funds for the study under the project AxSur.
References
 L.J. Gibson, M.F. Ashby, Cellular Solids: Structure and Properties, 2nd. edn., Cambridge University Press, Cambridge, 1997 [Google Scholar]
 M.F. Ashby, A.G. Evans, N.A. Fleck, L.J. Gibson, J.W. Hutchinson, H.N.G. Wadley, Metal Foams: A Design Guide, ButterworthHeinemann, Burlington, 2000 [Google Scholar]
 R. Lakes, Foam structures with a negative Poisson’s ratio, Science 235, 1038–1040 (1987) [CrossRef] [PubMed] [Google Scholar]
 K.E. Evans, M.A. Nkansah, I.J. Hutchinson, S.C. Rogers, Molecular network design, Nature 353, 124 (1991) [CrossRef] [Google Scholar]
 J. Simpson, Z. Kazanci, Crushing investigation of crash boxes filled with honeycomb and reentrant (auxetic) lattices, ThinWalled Struct. 150, 106676 (2020) [CrossRef] [Google Scholar]
 N. Abdullah, M. Sani, M. Salwani, N. Husain, A review on crashworthiness studies of crash box structure, ThinWalled Struct. 153, 106795 (2020) [CrossRef] [Google Scholar]
 V. Carneiro, J. Meireles, H. Puga, Auxetic materials – a review, Mater. Sci. (Poland) 31, 561–571 (2013) [CrossRef] [Google Scholar]
 T. Fíla, P. Zlámal, O. Jiroušek, J. Falta, P. Koudelka, D. Kytýř, T. Doktor, J. Valach, Impact testing of polymerfilled auxetics using split Hopkinson pressure bar: impact testing of polymerfilled auxetics, Adv. Eng. Mater. 19, 1700076 (2017) [CrossRef] [Google Scholar]
 D. Ruan, G. Lu, B. Wang, T.X. Yu, Inplane dynamic crushing of honeycombs – a finite element study, Int. J. Impact Eng. 28, 161–182 (2003) [CrossRef] [Google Scholar]
 L.L. Hu, T.X. Yu, Dynamic crushing strength of hexagonal honeycombs, Int. J. Impact Eng. 37, 467–474 (2010) [CrossRef] [Google Scholar]
 L.L. Hu, T.X. Yu, Mechanical behavior of hexagonal honeycombs under lowvelocity impact – theory and simulations, Int. J. Solids Struct. 50, 3152–3165 (2013) [CrossRef] [Google Scholar]
 X. Hou, Z. Deng, K. Zhang, Dynamic crushing strength analysis of auxetic honeycombs, Acta Mechanica Solida Sinica 29, 490–501 (2016) [CrossRef] [Google Scholar]
 S.R. Reid, C. Peng, Dynamic uniaxial crushing of wood, Int. J. Impact Eng. 19, 531–570 (1997) [CrossRef] [Google Scholar]
 L.J. Gibson, M.F. Ashby, G.S. Schajer, C.I. Robertson, The mechanics of twodimensional cellular materials, Proc. Roy. Soc. Lond. Ser. A Math. Phys. Sci. 382, 25–42 (1982) [Google Scholar]
 Altair^{®} radioss^{®}. URL https://www.altair.com/radioss/ [Google Scholar]
 A.S. Khan, H. Liu, Variable strain rate sensitivity in an aluminum alloy: response and constitutive modeling, Int. J. Plasticity 36, 1–14 (2012) [CrossRef] [Google Scholar]
Cite this article as: A. Petit, D. Beaujouan, R. Delille, F. Chaari, A. Langlet, On the response of reentrant auxetic structure under dynamic crushing, Mechanics & Industry 24, 43 (2023)
All Tables
All Figures
Fig. 1 View of the complete auxetic structure (“AS”). 

In the text 
Fig. 2 Reentrant auxetic cell configuration. 

In the text 
Fig. 3 Schematics and parameters defining a shock propagation into cellular bar after a rigid mass impact. 

In the text 
Fig. 4 Finite element impactor velocity results υ for several mesh sizes (from one to five elements in the thickness t_{thick}) at initial velocity υ_{0} = 200 ms^{−1}. 

In the text 
Fig. 5 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 35 ms^{−1}. (a) AS deformation at t = 0.77 ms. (b) Stress and plate velocity υs. time. 

In the text 
Fig. 6 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 50 ms^{−1}. (a) AS deformation at t = 0.45 ms. (b) Stress and plate velocity υs. time. 

In the text 
Fig. 7 AS crushing simulations results for impact conditions: M = 4 g and υ_{0} = 200 ms^{−1}. (a) AS deformation at t = 0.17 ms. (b) Stress and plate velocity υs. time. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.