High resolution simulations of energy absorption in dynamically loaded cellular structures

Cellular materials have potential application as absorbers of energy generated by high velocity impact. CTH, a Sandia National Laboratories Code which allows very severe strains to be simulated, has been used to perform very high resolution simulations showing the dynamic crushing of a series of two-dimensional, stainless steel metal structures with varying architectures. The structures are positioned to provide a cushion between a solid stainless steel ﬂyer plate with velocities ranging from 300 to 900 m/s, and an initially stationary stainless steel target. Each of the alternative architectures under consideration was formed by an array of identical cells each of which had a constant volume and a constant density. The resolution of the simulations was maximised by choosing a conﬁguration in which one-dimensional conditions persisted for the full period over which the specimen densiﬁed, a condition which is most readily met by impacting high density specimens at high velocity. It was found that the total plastic ﬂow and, therefore, the irreversible energy dissipated in the fully densiﬁed energy absorbing cell, increase (a) as the structure becomes more rodlike and less platelike and (b) as the impact velocity increases. Sequential CTH images of the deformation processes show that the ﬂow of the cell material may be broadly divided into macroscopic ﬂow perpendicular to the compression direction and jetting-type processes (microki-Communicated netic ﬂow) which tend to predominate in rod and rodlike conﬁgurations and also tend to play an increasing role at increased strain rates. A very simple analysis of a conﬁguration in which a solid ﬂyer impacts a solid target provides a baseline against which to compare and explain features seen in the simulations. The work provides a basis for the devel-opment of energy absorbing structures for application in the 200–1000 m/s impact regime.


Introduction
A common requirement in the design of many structures and vehicles is the ability to ameliorate the effects of projectile impact or collisions. In particular, scenarios in which a high velocity projectile impacts a target may require stresses to be minimised in either the projectile or the target or, sometimes, both. This protection is often accomplished by incorporating a "cushion" of porous or cellular materials between the projectile and the target. This energy absorbing component can act to dissipate the excess energy via irreversible dynamic deformation processes.
Porous or cellular materials have demonstrated superior energy absorption under shock or impact loading when compared to monolithic materials [1,2], while also offering the benefit of improved strength to weight ratios. A variety of cellular structures have been investigated with the aim of improving understanding of compressive response and thus maximising energy absorption over a range of representative impact conditions. These include metal [3][4][5] and polymeric [6][7][8][9] foams, composites [10][11][12] and honeycombs [13][14][15]. One notable feature of this class of materials is the enhance-ment of crushing strength observed under dynamic loading due to inertial effects, as seen by Reid and Peng [16] in wood, Tan et al. [17] in aluminium foams and Xue and Hutchinson [18] and Wu and Jiang [19] in metallic honeycombs. This phenomenon has also been the subject of a numerical study by Liu et al. [20], the conclusions of which identify a critical velocity above which inertial effects become significant.
Recent investigations have expanded the range of structures to include lattice truss structures such as octet [21], Kagomé [22] or pyramidal [23] designs. For a review of the various core topologies and manufacturing methods employed in the creation of periodic cellular structures the reader is referred to Wadley [24] and Lefebvre et al. [25]. Many of these studies fielded the cellular material in a sandwich panel configuration whereby the deflection of the face sheets can be used as a metric for measuring energy absorption.
Advances in manufacturing technology, in particular the flexibility of design afforded by additive manufacturing, have led to a significant increase in the range of cellular geometries available, for example the gyroid structure as described by Yan et al. [26], or microlattice designs tested by Mines et al. [27]. Such complex structures can be realised due to the ability of additive processes to build components as an aggregate of thin layers of material, using a three-dimensional (3D) Computer Aided Design (CAD) model as a reference. As a result of these developments, there is a growing requirement to compare the relative merits of different structures and identify the most appropriate for a given loading scenario.
As described above, previous research on the response of cellular structures, both additively and conventionally manufactured, has improved our understanding of energy absorbing processes in this important class of materials and has led to refinement of cellular material design. In particular in many of the past studies, computer simulation has been used both to highlight the important deformation phenomena and to optimise material design. Clearly, an important prerequisite of such simulations should be to determine how well the chosen computational technique matches experiments performed in the appropriate loading regime. It is especially important to determine the validity of the computer code during the initial phase of the compression, in which the plastic wave generated in the lattice by the impact propagates through the full thickness of the sample. This is because, if the code cannot model the early phase of the deformation, it is very unlikely that it will model the later phases correctly. However in many previous studies the response of the system is affected by stress (usually release) waves, propagating from the lateral outer boundaries of the sample before the initial plastic wave traverses the thickness of the specimen. For example, consider a disc-shaped cellular metal sample of diameter 200 mm and thickness 40 mm impacted by an impactor travelling at 0.2 mm µs −1 . Typically in a metal release waves propagate at 5 mm µs −1 . Therefore, following impact, release waves will travel the 100 mm from the outer boundary to the centre of the system in 20 µs. During this 20 µs period the flyer will travel 4 mm or about 10 % of the thickness of the specimen. It follows, therefore, that for most of the period of compression, the response of the sample, including the region near the centre of the sample, depends on the position of the outer perimeter. Therefore, to allow meaningful comparison between experiment and simulation the whole sample must be included in the model. Further if the detailed response of the individual struts or other features of the cells is needed, then, in general, a threedimensional (3D) code must be used. In summary, in the two or three-dimensional configurations employed in much of the previous work on this topic, the need to model the whole experiment plus the requirement to model the structure of the individual cells in 3D would require unrealistically high computational resources.

Advantages of one-dimensional (1D) geometry
Significant computing resources can be saved in the simulation of cellular structures if a "pseudo 1D" modelling scheme similar to that employed in simulations by Borg and Vogler [28] is adopted. Consider a cylindrical sample made up of a regular cellular structure in which the diameter of the sample is much greater than its thickness. It can be shown that, following impact, the time to reach full densification of the whole specimen is less than the time for release waves to reach the axis of the sample if the diameter of the specimen, the density of the sample and the impact velocity are high and the thickness of the specimen is low. Under such conditions, the compression of the cells at the centre of the disc may, on a macroscopic or continuum scale, be considered 1D. At the scale of individual cells such a system is, of course, not 1D but 3D. However, cells near the centre of the disc will experience the same stress field as each other and therefore every cell will respond in the same way. In this situation, modelling a single cell or column of cells provides the response of the whole of the region over which 1D conditions apply. Note that since, in general, the deformation of an individual cell is three dimensional, rather than one dimensional, the configuration should, perhaps, be more correctly termed "Pseudo 1D". However, for convenience we will refer to this modelling concept as "1D" during the remainder of this paper. Clearly significant computing resources are saved by confining modelling to a single cell or small group of cells.
In an earlier paper, we described a study in which lattices consisting of an array of intersecting stainless steel rods were manufactured using an additive manufacturing technique known as Selective Laser Melting (SLM) [29]. The structure consisted of intersecting 316L stainless steel rods. The sam-ples, which were accelerated using a 100 mm gas gun, were discs of diameter ∼100 mm and thickness 6 mm. The cells which made up the lattice measured 1 mm × 1 mm × 1 mm and the average density of the lattice was 64 % of solid. Impact velocities ranged from 0.3 to 0.7 mm µs −1 . A rough estimate of the time for which 1D conditions persist in this geometry is as follows. The time for elastic release waves from the perimeter to reach the centre of the system was ∼10 µs. At an impact velocity of 0.5 mm µs −1 the time for full densification of the lattice was of the order of 6 µs. Therefore a useful amount of compression was achieved before the effect of release waves from the outer of the specimen was felt at the centre of the system. Note that the densities of the sample, and the impact velocity regime in our work are somewhat higher than in many previous studies of the dynamic response of cellular materials. In this regime the behaviour of the material tends to be dominated by shock formation and propagation and as such relates to severe blast and impact scenarios.
The experiments were simulated using the Sandia National Laboratories Code CTH. CTH is a multi-dimensional, multimaterial Eulerian-based hydro-code developed by Sandia National Laboratories [30]. The Eulerian meshing scheme enables it to cope with the simulation of large deformations and material distortions. The incorporation of the laws of conservation of mass, momentum and energy, together with a suite of material models which allow description of most states of matter normally encountered in shock physics, make the code appropriate for calculating the response of materials loaded under shock and high strain rate conditions. In effect the code is optimised for simulating situations, encountered at high impact energies, in which inertial (or material acceleration) effects dominate. The effects of material strength on material response may also be modelled using strength models such as Steinberg-Guinan [31] and . CTH has been successfully applied to the study of granular ceramic materials impacted in a velocity regime similar to that discussed in this paper [28]. The CTH code contrasts with implicit finite element codes in which the main aim is to simulate the response of materials in the low strain rate regime where the response of material is dominated by strength rather than by inertial effects.
A summary of the findings of our earlier study is as follows. Both "full" and "cut down" three-dimensional structural CTH calculations were run in which the stereo lithography (stl) files used to manufacture the samples were imported into CTH. In these simulations the detailed, dimensioned, structure of the whole cellular sample was modelled but in the "cut down" simulations only a 2 cell × 2 cell × 6 cell column of cells was modelled thereby much reducing the computational requirements. As expected, it was found that, since the conditions in the sample were close to 1D, there was little difference in the predictions of the full and cut down simulations. It was found that both sets of 3D simulations gave a reasonable match to all of the features observed in the experiments including the oscillations which were superimposed on the first and second velocity plateaus at the back surface of the solid steel target. In particular, the simulations provided insights into the processes at play in the shocked lattice. As the compression wave reaches a cavity, a release wave followed by a compression wave is generated. The array of cavities in the structure thereby generates stress oscillations which propagate in all directions. As these waves reach the back surface of the target, they cause the velocity oscillations observed in the experiment. In summary, our study gave us confidence in the capability of CTH to predict material response in this regime and gave some important insights into the mechanisms at play. However, because the structure was rather complex it was not ideally suitable as a basis for a systematic study of the effects of cell design on energy absorption. It was felt, therefore, that our proposed study of the effect of cell shape on energy absorption needed to be based on a structure that was much simpler than that used in the earlier study.
In this paper we describe the use of the previously validated CTH models to perform an extended study of a simplified lattice structure. In the study described here, the simulations are confined to a single 6 mm × 6 mm twodimensional cell allowing the distortions of the structure to be imaged at higher resolution than in the earlier study. The aim is to provide a clearer understanding of the energy absorbing processes occurring in metal-based cellular samples during impact than was provided by the relatively complex structures in the earlier study and, thereby, to support the optimisation of cellular structures for application as energy absorbers.

Modelling set-up
In the study described here simulations were confined to a single cell. Further significant savings on computer resources were made by adopting a two-dimensional rather than a three-dimensional configuration. Figure 1 shows an idealised two-dimensional configuration, chosen as the basis for the current study. The impactor is directed vertically upwards. The energy absorbing elements (EAE) and the impactor and target are considered to consist of infinitely long bars extending in and out of the paper. The outlines shown in Fig. 1 represent sections through these bars. For example the grey rectangle illustrates the orientation of the "rod" variant (also shown on the extreme right of Fig. 2). The flyer and the target are rectangles of dimensions 6 mm × 18 mm, again effectively extending an infinite distance in and out of the paper. In order to provide an improved understanding of the effects of lattice structure on energy absorption, it was necessary The five structures assessed using CTH. In all cases the 18 mm impactor approaches the cell from below. The 18 mm target block is positioned above, and in contact, with the EAE to evaluate a systematic range of potential strut orientations. Therefore, simulations of the five alternative EAE structures depicted in Fig. 2 were performed. In all cases the solid occupies 40 % of the total volume of the cell. The flyer impacts from below and the target, or acceptor, is above the EAE. For convenience the EAE on the extreme right ( Fig. 2e) is labelled "rod". However, since this is a 2D setup the "rod" element actually consists of a plate extending into the paper. The structure on the extreme left ( Fig. 2a) also consists of a plate extending into the paper but, it is rotated by 90 • relative to the structure on the extreme right. The structure in the centre, labelled "cross" consists of two plates at 90 • to each other. The inclusion of the rodlike, cross and platelike structures provides a simple transition from rod to plate structure. Since the shape and mass of the unit cell is held constant the strut thickness, d, depends solely on the angle α. It can be shown that the porosity, φ, defined as the fraction of the cell occupied by solid material, is given by: for options a and e for options b and d for option c where α, d and L are as defined in Fig. 2.
It is envisaged that, theoretically, each of the unit cells depicted in Fig. 2 could be repeated horizontally and vertically to provide a macroscopic structure. However, note that, with the exception of the "cross", the macroscopic structure that would be formed by repeating the unit cells studied here would not form practicable, interconnected, material Structural parameters for the five EAE designs are shown in Table 1.
The steel components were modelled using a Mie-Gruneisen equation of state and the Steinberg-Guinan (SG) strength model, the same approach used in the previously reported study. The SG model is implemented as follows. The yield stress, Y , is given by: where Y 0 is the yield stress in the reference state, p is the pressure, η the density normalized by the reference solid density, and T is the temperature. ε PL , is the equivalent plastic strain, and β, n, g, h are material-dependent parameters. The first term on the right of this equation represents work hardening, and the second term pressure hardening and thermal softening. The final, exponential, term produces a rapid decay in flow stress as the melt temperature is approached. The material parameters used in the current study are shown in Table 2. They were the same as those which have been used in the earlier investigation apart from the yield strength. The CTH default value for Y 0 (0.34 GPa) was used in the current calculations whereas in the previous study an elevated value of 0.64 GPa, based on the measured strength of SLM steel was used. The mesh size in the study reported here was 0.02 mm as compared with 0.05 mm in the earlier work. In this work, the number of cells through the thickness of each strut varied between 48 and 120 depending on which EAE was being investigated; this compares with ∼8 cells through each strut in the earlier study [29]. Key objectives of the study to be described are to obtain insights into the mechanisms at play during the interaction between the flyer and the cushioned target and to estimate from simulations which of the chosen variants maximises the irreversible energy generated in the EAE and thereby minimises the energy transferred to the target. Therefore CTH was used to generate time sequences showing materials, plastic flow, internal and kinetic energy, pressure, temperature and flow stress. The total kinetic and total internal energy of the flyer, EAE and target were also output from the code. As mentioned in Sect. 1, the regime of prime interest in this work is the period during which the porous structure (the EAE) compresses due to the transmission of the first plastic wave. However, the transmission of the wave into the solid target is also of great interest because of its potential to evaluate the accuracy of the simulations by monitoring the transmitted wave in future experiments.

Simulations
Simulations provide a picture of the mechanisms at play during the deformation process following impact. For example Fig. 3 shows a sequence generated by a 0.3 mm µs −1 impact on a rodlike EAE. The left-most frame shows a material plot at time zero. On the distance scale displayed on the extreme left, the solid stainless steel flyer extends from 2 to 20 mm; at time zero it has a velocity of 0.3 mm µs −1 directed ver-tically upwards. The EAE extends from 20 to 26 mm and the solid steel target extends from 26 to 44 mm. In the rest of the frames, the colours represent velocities with magnitudes indicated by the scale at the extreme right. The times of the frames in microseconds from impact are shown at the top of each frame. From this sequence we judge that full densification of the EAE occurs at ∼14 µs. The estimated trajectories of the waves propagating through the system are shown by the superimposed black lines. The solid black line which extends from the 2.8 to the 14.0 µs frame shows the estimated trajectory of the front face of the flyer. The dashed lines are a schematic representation of the wavelets generated as the flyer interacts with the EAE. Reference to the velocity scale shows that, by the time full densification of the EAE occurs, the velocity of the flyer has reduced from 0.3 to ∼0.24 mm µs −1 .
As expected the different EAEs respond differently. Figure 4 shows internal energy, plastic strain and temperature plots for the five EAE variants just after full densification for a simulation at 0.3 mm µs −1 . It is clear that the plastic strain and, therefore, internal energy absorbed by the EAE tends to increase as the configuration is varied from plate to rod but that there is little difference between the rodlike and rod variants.
The energies in each computational cell of each of the three regions (flyer, EAE, target) can be summed to give the total energy for a 1 cm length (measured into the paper) of each region or component as a function of time. As an illustration, Fig. 5 shows plots of internal, kinetic and total (internal plus kinetic) energy for a system with a cross EAE impacted at 0.3 mm µs −1 . The kinetic energy in the impactor reduces to a minimum and then rises again. Meanwhile the internal energy exhibits a "double peak" profile which is also seen in the target. The kinetic energy in the target exhibits a single peak over the time frame of the graph. The general shapes of the profiles, which were found to be similar for all of the EAEs, will be discussed further in Sect. 5. Figure 6a shows the total internal energy delivered to the different EAEs by a 0.3 mm µs −1 impact. It is seen that the internal energy in all of the EAE initially ramps up and then reaches a plateau. The onset of the plateau corresponds roughly to the time at which the EAE reach full densification. The greatest energy absorption is achieved with the rod and rodlike EAEs and the least energy is absorbed by the plate and platelike EAEs. The cross structure lies between these two. Corresponding results for a 0.9 mm µs −1 impact are shown in Fig. 6b. It can be seen that, again, the energy rises roughly linearly to a plateau. Again the onset of the plateau corresponds to full densification (see also the sequence in Fig. 14). There are two significant differences between the 0.3 and 0.9 mm µs −1 results. First, the internal energies dissipated in the 0.9 mm µs −1 simulation are at least twice the corresponding energies in the 0.3 mm µs −1 simula-  tions. Second, in the higher velocity simulation the energy generated in the cross structure is similar to that generated in the rod and rodlike structures whereas at 0.3 mm µs −1 the cross is intermediate between the rod and plate structures.
One of the objectives of this study was to determine the factors affecting the energy transferred to the target at late time. This is difficult to assess thoroughly using a Eulerian code for two main reasons. First, as explained in the discus-sion, the waves generated by impact will reflect back and forth in the longitudinal direction for a significantly longer time than it was possible to run the simulations in this study. Second, it will be noted later that the model used in the simulations does not allow the components to separate when tensile forces develop between them. It follows that the simulations become invalid before energy exchange between the target and the other components in the system ceases. However, despite this difficulty, we believe a useful guide to the  Figure 7a shows the total (kinetic plus internal) energies in the target following a 0.3 mm µs −1 impact onto a target protected by the five EAEs under investigation. Note that with the rod EAE the energy in the target increases at ∼3 µs after impact. We assume that this occurs as the elastic wave generated by the first impact arrives at the target. By contrast the other EAEs do not transmit a significant elastic wave. In each case the energy rises to a peak and the timing of the peak broadly corresponds to the moment that the release waves reflected from the back face of the flyer and the front face of the target meet. It is clear from Fig. 7a that the peak energy delivered to the target is less for the rod and rodlike EAE than for the plate and platelike EAEs. Figure 7b shows the target total energies for a series of simulations of 0.9 mm µs −1 impacts. In these examples the simulations were run as far as the peak energy but no further. Again the position of the peak corresponds to the time at which the stress at the impactor EAE interface is expected to go negative. In this high velocity example, it can be seen that the choice of EAE makes only a very small difference to the energy transferred to the target.

Outline
A series of simulations have been performed with the aim of ranking a series of Energy Absorbing Elements (EAEs) or "cushions" in terms of their effectiveness in minimising the energy transmitted from an impactor with velocities in the range 0.3-0.9 mm µs −1 to an initially stationary target. The results have been presented in the preceding sections. In this section, the aim is to support an understanding of the results by analysing the way the internal energy generated by the impacts is distributed.

Shock analysis
Before examining the detailed response of our chosen system (as depicted in Figs. 1, 2) it is illustrative to consider the energy distribution in a very simple system from which the EAE has been removed. Figure 8a shows a schematic distance vs. time plot for a 0.3 mm µs −1 system. For simplicity it has been assumed that energy dissipation is zero and, therefore, that the system remains isentropic. In this situation, all energy transfers are reversible. Therefore the energy is in the form of either kinetic energy or compression energy; effectively the material is treated as elastic. In this situation the wave velocity remains constant (at 4.8 mm µs −1 say) and the particle velocity in between the shocks generated by the impact in the flyer and target is half the impactor velocity. On wave reflection from the free surfaces of the flyer and the target, the surface velocity is exactly twice the in-situ particle velocity. It is also convenient to assume that when the release waves from the free surfaces meet the interface between the flyer and the target, the target and the flyer do not separate. This last assumption, albeit nonphysical, allows a direct comparison with the computer simulations, which were run in a mode which does not allow separation of the components. This aspect of the modelling will be discussed later. The assumptions listed above allow the shock parameters in the flyer and target (e.g., density, pressure, particle velocity and internal energy) to be estimated using basic shock theory. Knowledge of the particle velocity distribution in the impactor and target allows the kinetic energies (KE) of the impactor and target to be derived. Assuming the density of the steel is 7.9 g cm −3 , and that the flyer measures 0.6 cm wide × 1.8 cm long and extends 1 cm into the page, the volume of the flyer is 1.08 cm 3 and its mass, m, is 8.532 g. The initial kinetic energy of the flyer, E 300 , is given by where V is the flyer velocity. Table 3 lists the initial kinetic energy of the flyer and the specific initial kinetic energy of the flyer for impact velocities ranging from 0.3 to 0.9 mm µs −1 .
The average specific energy of the system is calculated by dividing the initial kinetic energy by the total mass of the system. Knowing the position of the shock as a function of time allows the KE distribution to be deduced. Retaining the 0.3 mm µs −1 "experiment" as an example, the illustration in Fig. 8a shows that at 0 and 15 µs the whole impactor (or flyer) has a velocity of 0.3 mm µs −1 and at 7.5 and 22.5 µs the whole flyer is stationary. The energy distributions corresponding to Fig. 8a are shown in Fig. 9. It can be seen that the KE of the flyer oscillates between a maximum of 383 J at 0, 15 µs etc. and a minimum (of zero) at 7.5, 22.5 µs etc. It can also be seen that (since kinetic energy varies as velocity 2 ) at 3.75 µs (the wave transit time) the kinetic energy of the whole system (flyer plus target) is 192 J (half the total KE at 0 µs). Therefore at 7.5 µs (and 15 µs, etc.) half the total energy of the system is in kinetic energy and half is in internal energy. Of course, in our simple elastic example, there is no increase of entropy, therefore, at this stage, all of the internal energy is in the form of mechanical compression energy. The specific internal energy in the initial shock generated by the impact is 11.25 J g −1 . Note that in both the impactor and the target, the internal energy exhibits a double peak form. The pressure and specific volume in the shock generated by the impact can be calculated if the Hugoniot of the steel is known. Assuming a linear shock velocity (U s ) vs. particle velocity (u p ) relationship and using the constants from Table 2 gives the Hugoniot: From the laws of conservation for shocks the pressure, p, and volume, v, behind the initial shock are given by: and It can be shown using a pressure versus particle velocity construction that the pressure and specific volume in a shock generated by impacts at a range of velocities are as listed in columns 2 to 4 of Table 4 as are the corresponding normal strains given by: Note from Fig. 8a that from 7.5 to 15 µs the interface between the two components is predicted to be in tension. However it is assumed in the illustration that the components will remain in contact. This assumption allows us to compare the predictions of the simple model with the predictions of CTH. Since in reality the components would probably separate rather than go into tension the predictions are not valid after the time the release wave from the front and back surfaces of the system meet (in our example beyond ∼22 µs after impact time). However the results of simulations beyond this time have been presented to permit a more complete comparison between analytical and simulated results. Note that there are models in CTH to treat tensile failure. We plan to investigate the effect of applying these models in future simulations; it is possible that these will affect the late time behaviour of the simulations.

Plasticity analysis
Although the p, v states presented in Table 4 assume that the impactor and target were strengthless, it is useful to separate the internal energy that would have been generated if the material had finite strength into its component parts. Remember that the following analysis applies to a simple metal-on-metal impact with no EAE present.
The red line in the schematic stress vs. strain plot in Fig. 10 depicts the pressure vs. strain response of a material subjected (at a low strain rate) to uniaxial strain. It is assumed that the pressure vs. strain relationship may be simplified as; where K is the bulk modulus, and ε n is given by:  The dashed black line shows the stress strain trajectory corresponding to release from X . The total energy, E, generated when the material is strained to ε n is given by (12): However E is made up of the following three components. The elastic energy, E e , is defined as the reversible energy associated without change of shape of the material with change of volume. The hydrodynamic energy, E h , is also reversible and is the energy associated with change of volume without change of shape. Finally the plastic energy, E p , is the irreversible energy associated with change of shape without change of volume. For the very simple case of uniaxial strain of an EPP material assumed in our example, it is straightforward to estimate the component energies. In fact, as in (13)- (15), E e , E h and E p are given by the grey, blue and yellow areas in Fig. 10. (Note that L is the longitudinal modulus): Estimated values of the three components of internal energy associated with simple steel-on-steel uniaxial loading are presented in Table 4. Note that the E e and E p both depend on the flow stress Y 0 . In the material model used in CTH the flow stress can vary from Y 0 to Y max . This is accounted for in Table 4 by estimating the minimum and maximum values of E e and E p . The expected effect of adding an energy absorbing element to the basic configuration shown in Fig. 8a is illustrated schematically in Fig. 8b. The leading edge of the flyer interacts with the leading edge of the EAE, and high pressure wavelets are driven back into the flyer. These wavelets reflect from the back face of the flyer leading to a gradual reduction in its velocity. Simultaneously, we would expect that wavelets driven through the structure of the EAE gradually impart movement and energy to the target.

Analysis of simulations
The wave propagation histories in the simplified model and the energy distributions derived from them furnish a baseline against which to compare the CTH simulations and provide a qualitative explanation of some of the features observed in the energy vs. time profile generated by the simulations. For example Fig. 11a, b shows computed energies for a system without an EAE for comparison with the idealised configuration depicted in Fig. 8a. It is seen that the periodic time variation in the internal energy profile generated by the simulations corresponds qualitatively to the oscillations between the back face of the flyer and the front face of the target. However, unlike the simple elastic model, the CTH simulation predicts that the amplitude of the oscillatory waves decreases with time. This is because, in the code, the passage of the compression waves converts some of the initial kinetic energy to irreversible energy by shock propagation and plastic flow.
In the example in Fig. 12 the time variation of the computed energies in an 0.3 mm µs −1 impactor following interaction with a rodlike EAE is shown for comparison with the predictions of the simplified configuration in Fig. 8a. It is seen that the shape features obtained from the simulated energy vs. time plots are similar to those implied by the simple analysis. In particular, the comparison provides a clear As with the simple uniaxial strain configuration discussed in the previous section, the total internal energy, E, generated in the simulations of the cushioned systems may be regarded as the sum of three components. Remember the compression energy, E c , is the energy associated with changes in volume without change of shape and the elastic energy, E e , is the energy associated with shear distortion prior to yielding. Compressional and elastic energy are reversible and will be available to do work on the system as the forces on the system are relaxed. The plastic energy, E p , is the energy associated with plastic flow following yielding. Plastic energy is not reversible but instead is converted to heat and not, therefore, returned to the system at late time. It follows that maximising the plastic energy generated in the EAE minimises the energy reflected back into the impactor and transmitted into the target. Ideally, the plastic energy (or plastic work) generated in each of our EAE variants during densification would be output from the CTH code. Although CTH does not appear to have the capability to output plastic energy it can provide plastic strain. For example Fig. 13 shows the plastic strain averaged over all the cells in rod EAEs generated during loading by projectiles with velocities ranging from 0.3 to 0.9 mm µs −1 . In each case there is a roughly linear increase in average plastic strain followed by a plateau. The onset of the plateau corresponds to full densification of the EAE. The simulations clearly show that, as impact velocity increases, the amount of plastic flow corresponding to full densification increases. Since the stress at which this flow takes place would be expected to increase with pressure, the increase in plastic work with velocity would be even more marked than the trend shown in Fig. 13. An alternative approach is illustrated by the data in Table 5 which lists the specific internal energies generated in the EAEs at 0.3 and 0.9 mm µs −1 . Also listed in Table 5 are the elastic and compression energies estimated for shock waves generated by steel-on-steel impact without an EAE present. We assume that these latter values represent upper limits on the elastic and compression energy generated in an EAE. It can be seen that, at both impact velocities and for all EAE variants, the elastic component is small compared with plastic energy. For the plate configuration the estimated compression energy is an appreciable fraction of the total internal energy which is not surprising since the plate configuration is geometrically similar to a simple steel-on-steel impact. However it is clear that in the other EAE configurations the total internal energy is dominated by the plastic component of the internal energy. It follows that, in most of the simulations run in support of this investigation, the total internal energy may be taken as an approximate measure of the plastic energy. This contention is supported by the observation in Fig. 4 that the patterns showing the distribution of plastic strain appear very similar to the corresponding internal energy distributions. A possible exception to this trend is the plate EAE where, as might be intuitively expected, most of the internal energy is in the form of hydrodynamic compression energy.

Conclusions from energy analysis
It is clear, and intuitively obvious, that the densification of the rod and rodlike EAEs requires more plastic flow than the platelike and plate options. Consequently, the former absorb more irreversible energy and so transmit less energy to the target.
The amount of plastic flow associated with full densification also depends on the impact velocity. As impact velocity increases the energy analysis indicates that the amount of plastic flow and therefore the irreversible energy generated increases. Additional simulations were run to explore this phenomenon.
The sequence in Fig. 14 maps internal energy density for a rodlike cell impacted at 0.9 mm µs −1 . Full densification of the cell occurs at ∼4.8 µs, after which time there is very little change in the internal energy distribution in the vicinity of the EAE. Note that intense flow occurs at the edges of the cell as mushrooming of the angled plates produces jets which eventually impinge on the end-face of the target. Note that the small mesh size (0.02 mm) employed in this study has allowed the fine detail of the jets and their interaction with the cell walls to be resolved. The highest specific energies seen in the graphics are generated as a result of this jet impact. It can be seen that, following impact, plastic waves propagate into the flyer and the target.
Comparing Fig. 6a with Fig. 6b shows that increasing the impact velocity increases the energy absorption in the EAE and the relative amplitudes of the energy profiles in the target, shown in Fig. 7, are consistent with the amplitudes of the internal energy profiles in the corresponding Table 5 CTH simulations: peak volume average specific internal energies generated in the EAE (J g −1 ) Impact velocity (mm µs −1 )  EAEs. Of particular interest here is the observation that, at 0.3 mm µs −1 , the cross structure lies about halfway between the platelike and rodlike structures, but at 0.9 mm µs −1 there is little difference between the cross and the rod and rodlike structures. The differing effectiveness of the cross structure as the impact velocity changes may be explained by comparing the deformation mechanisms of the cross and rodlike structures at 0.3 and 0.9 mm µs −1 . The simulated response of the cross and rodlike structures at 0.3 and 0.9 mm µs −1 are shown for comparison in the material plots in Fig. 15a-d. The frames in Fig. 15a shows that, according to the simulations, at 0.3 mm µs −1 the cross structure closes by axial movement of the parts of the cell driven by the impactor. By contrast, as shown in Fig. 15b, at 0.3 mm µs −1 the rodlike cell closes predominantly by lateral flow, a process accompanied by significantly more plastic deformation than with the cross variant. Therefore, at 0.3 mm µs −1 , more energy is absorbed by the rodlike than by the cross variant. At 0.9 mm µs −1 different mechanisms come into play. It is seen from Fig. 15c, d that both the cross and the rodlike variants deform by a concentration of lateral flow in the region of the cell adjacent to the impactor. This produces a mushroom feature which occu-pies the whole width of the cell and intense lateral flow and jetting predominate. Plots showing plastic strain in cross and rodlike EAEs impacted at 0.9 mm µs −1 are shown in Fig. 16.
As indicated, in the cross structure concentrations of flow develop between the impactor and the diagonal struts and between the diagonal struts themselves. Particularly prominent is a mushroom feature located near the original centre of the cross, which we believe is formed by penetration of the cross by a jet of impactor material. These concentrations of deformation, or microkinetic flows, which are not observed when the cross variant is impacted at 0.3 mm µs −1 are the probable reason for the much-increased energy absorption as velocity increases. In summary the similarity of the mechanisms in the two EAEs results in the energy generated in the EAEs by 0.9 mm µs −1 impacts being similar. The change in deformation behaviour at higher velocities to one dominated by intense plastic flow parallels that observed by Nesterenko in studies of granular materials [33]. Nesterenko proposed that the compression of granular porous material at high rates consisted of distortions of the grains themselves combined with smaller scale deformations resulting in features such as jets and localised melting at grain boundaries. Separating the collapse process into two phases enabled the concept of microkinetic energy to be introduced, defined as the energy dissipated via intense plastic flow during pore collapse. In our simulations, it appears that the occurrence of the plastic flow significantly increases the energy absorbing capability of a particular structure, representing the microkinetic energy component of the total dissipated energy. Figure 17 shows the results of loading a rod structure at velocities increasing from 0.3 to 0.9 mm µs −1 . The simulations clearly show that as the impact velocity is increased beyond 0.5 mm µs −1 there is a transition from a mechanism dominated by roughly uniform lateral flow to one domi-  [28]. These workers did not observe microjetting in their simulations. Our work is consistent with their findings in that microjetting was observed at impact velocities of 0.7 and 0.9 mm µs −1 but not at 0.3 and 0.5 mm µs −1 .
Note that intense flow concentrations have been observed both in the microjets and in the material with which the jets interact. Typically, in the simulations depicted in Figs. 14, 15, 16 and 17, flow concentrations have length scales extending down to the 100 µm regime. Although we have not yet made a systematic study of the effect of mesh size on jet formation and interaction, we have observed that in relatively coarse mesh simulations the jets appear much less well defined with a tendency to appear shorter and broader. There-fore we assume that coarsening the meshes from their current value of 20 µm would result in significant loss of the detail associated with the microkinetic flows.
Although the simulations have clearly demonstrated that increasing the impact velocity increases the energy absorbed by the EAE, they have also shown that at the highest velocities, the existence of a "single cell" cushion makes negligible difference to the energy transmitted to the target (see Fig. 7b). The reason for this is that the irreversible energy absorbed by the EAE is only a small fraction of the total energy delivered to the system by the impactor. From Table 3 the kinetic energies of the impactor are 384 and 3455 J at 0.3 and 0.9 mm µs −1 respectively. Multiplying the specific energies listed in Table 4 by the EAE mass (1.137 g) gives the maximum plastic energy generated in the EAEs as 143 J and 367 J at 0.3 and 0.9 mm µs −1 impacts respectively. The fraction of the input energy dissipated in the EAE is, therefore 37 and 10 % respectively. We note that the effectiveness of an energy absorbing cushion of the type investigated here could be increased by increasing its axial length. This could be achieved by stacking several EAE cells together along the impact direction.

Summary and conclusions
Porous materials have potential as energy absorbers in shock and impact scenarios. Computer simulation is needed to assess the relative merits of different cellular structures and materials in different loading environments. We have noted that the most useful comparisons between simulation and experiment are those in which nominally one-dimensional, or uniaxial strain, conditions persist, at least until maximum (or full) densification of the porous material is reached. Adopting this configuration has allowed very fine mesh simulation (48-120 meshes through each strut). We note that 1D conditions are most readily achieved in high density (low porosity) high velocity systems. For this reason we chose to study the response of materials of 40 % solid density in the velocity range 0.3 to 0.9 mm µs −1 .
A series of two-dimensional cellular structures designated "Energy Absorbing Elements" (EAEs) have been studied. The main variable was the angle of the flat plates which form the structure. As the angle changed the thickness of the plates was modified to maintain the average density of the structure at 40 % of solid.
The Sandia National Laboratories code, CTH, was used to simulate a system in which an 18 mm steel flyer impacts a target consisting of a 6 mm EAE in contact with an 18 mm steel target. The very fine mesh scheme (0.02 mm) employed in the code allowed the mechanisms at play during the compaction of the different EAEs to be assessed in various ways. For example, time sequences showing the distribution of particle velocity, pressure, internal energy and plastic strain were plotted. Simulations run at the relatively low velocity of 0.3 mm µs −1 showed that all EAE variants were compressed to solid density in about the same time. However, significantly more plastic flow and internal energy was generated by rod and rodlike EAEs than in the plate and platelike EAEs. It was shown that at the low impact velocity, for most cell types the cells densify by simple broadening of the elements of the EAE. By contrast, as impact velocity is increased to 0.5 mm µs −1 and beyond, jetting processes play an increasingly important role. In effect, jetting increases the amount of plastic flow associated with the full densification of the porous material. The increased plastic flow is accompanied by an increase in the irreversible energy absorbed during the densification process. Quantitative comparison of the effectiveness of the different EAEs can be provided by using CTH to calculate the energy distribution in the sys-tem. For example the pictorial evidence indicated that more energy was generated in the rod and rodlike EAE than in the platelike versions, implying that the former variant would provide a more effective cushion than the latter. This conclusion was supported by CTH calculations in which the energies generated in each cell of each of the three components of the system were summed. For example, simulations at 0.3 mm µs −1 showed that, relative to a platelike EAE, the rodlike EAE led to a lower gain in total energy of the target.
A very simple analysis of a configuration in which an 18 mm steel flyer impacts an 18 mm steel target with no cushion included has been presented. The model provides a prediction of the time variation of the kinetic and internal energies within the flyer and the target, thereby forming a baseline against which to compare simulations of the cushioned systems. The approach provides a qualitative explanation of key features of the simulated energy distribution in the cushioned system and allowed us to estimate the fraction of internal energy dissipated as irreversible energy in the cushioned systems. It was concluded that, apart from the plate variant, the great majority of the internal energy generated in the EAE is irreversible.
By applying the principles presented above when considering the merits of alternative energy absorbing structures, it should be possible to assess which is most suited to a specific impact regime. In conjunction with the close control over cellular material design which is possible using SLM, this approach offers the opportunity to manufacture and test cellular configurations which maximise generation of microkinetic energy.
An important future aim will be to extend the single cell study to multiple cells. In particular simulations will be performed to determine whether the energy absorption scales as cells are stacked in the direction of application of the load.
Further future investigations will consider how to verify the simulations described here against experimental data. One method of accessing this data is dynamic x-ray imaging of the impacted lattice, which will supply line of sight mass at different positions and times. See, for example, a recent study of lattice compaction by Eakins and Chapman [34]. Corresponding deformation profiles can then be extracted from 3D simulations for comparison with experiment.