Numerical evaluation of lightweight ultra high strength steel sandwich for energy absorption

Legislation regarding greenhouse gas emissions forces automotive manufacturers to bring forth new and innovative materials and structures for weight reduction of the body-in-white. The present work evaluates a lightweight ultra high strength steel sandwich concept, with perforated cores, for energy absorption applications. Hat-profile geometries, subjected to crushing, are studied numerically to evaluate specific energy absorption for the sandwich concept and solid hat-profiles of equivalent weight. Precise discretization of the perforated core requires large computational power. In the present work, this is addressed by homogenization, replacing the perforated core with a homogeneous material with equivalent mechanical properties. Input data for the equivalent material is obtained by analyzing a representative volume element, subjected to in-plane loading and out-of-plane bending/twisting using periodic boundary conditions. The homogenized sandwich reduces the number of finite elements and thereby computational time with approximately 95%, while maintaining accuracy with respect to force–displacement response and energy absorption. It is found that specific energy absorption is increased with 8–17%, when comparing solid and sandwich hat profiles of equivalent weight, and that a weight saving of at least 6% is possible for equivalent performance.


Introduction
The automotive industry faces long-term goals regarding reduction of greenhouse gas (GHG) emissions. Effort is therefore put into finding new design solutions and material selections, reducing the life-cycle ecological footprint through energy efficient production and reduced weight of e.g. body-in-white parts (BIW). At the same time, innovative design and new materials have to fulfill the demand regarding final product performance, such as crashworthiness.
Steel, a common material in many sectors, due to its high strength to cost ratio, has been utilized for manufacturing lightweight components by various methods of heat treatment. According to Kawajiri et al. [1], replacing mild steel with high strength steel has been one of the most effective material substitutions with regard to lowering GHG emissions. In the 1970's, press hardening was invented by Luleå university of technology and NJA (now SSAB) and industrialized by Plannja HardTech (now Gestamp Hardtech). In the press hardening process, a blank is austenitized and quenched, forming a mainly martensitic microstructure. The result is an ultra high strength steel (UHSS) with high shape accuracy, suitable for intrusion protection in vehicles. Due to the superior mechanical properties of UHSS, gauge thicknesses can be reduced, thus lowering the weight of the BIW. Karbasian and Tekkaya [2] describe the process steps of the methods, providing an insight into the method of press hardening. By controlling the temperature gradients of the tools, tailored microstructures are attainable for the hot formed components. Such components are desirable since the varying mechanical properties can be utilized for passenger safety by optimizing structural integrity and energy absorption. Extensive research has been carried out on UHSS for predicting the microstructure evolution during forming, fracture and fatigue. Models for predicting microstructure evolution has been carried out by several authors, e.g. Oldenburg and Lindkvist [3], and an overview of modeling of press hardening has been given by Oldenburg [4] and Merklein et al. Merklein et al. [5]. Research on modeling fracture in UHSS has been studied by several authors. Östlund et al. [6] suggest a strain based failure criteria, dependent on stress triaxiality and independent of element size, for shell elements. Good agreement was found between numerical and experimental tensile test results. The fracture toughness of low alloyed boron steel with varying microstructures was studied by Golling et al. [7], providing a better understanding of fracture initiation and allowing for manufacturing of components with improved crashworthiness. In Golling et al. [8] fracture toughness was studied under varying strain rates for different grades of high strength steels. It was found that the crack propagation resistance increased with the strain rate.
An additional method for reducing vehicle weight is incorporation of sandwich structures into the BIW, see e.g. the work by Hara and Özgen [9] and Balakrishnan et al. [10]. A typical sandwich consists of stiff face plates separated by a lightweight core. The ideal core has the lowest possible weight with sufficient stiffness to withstand transverse and shear loads, maintaining the initial distance between the face plates. A wide range of cores exist, such as foams, geometrical patterns, and solid cores.
Foam has been studied by several authors, and it was shown by Gibson [11] that the mechanical properties of foams are strongly influenced by the bulk material. Thus, the selection of bulk material is of importance when the foam is selected. A review of steel foam was done by Smith et al. [12], including manufacturing processes and structural applications. Additional work regarding steel foam was conducted by Park and Nutt [13], Smith et al. [14], Szyniszewski et al. [15] and Szyniszewski et al. [16]. Zhang et al. [17] studied aluminum foam experimentally, identifying parameters for constitutive modeling. Sulong et al. [18] studied advanced pore morphology foam elements in compression, using experiments and numerical methods. Xiao et al. [19] studied aluminum foam filled bumper beams for energy absorption. Micromechanical models were adopted by Marsavina et al. [20] for predicting mechanical properties of foam. Özer et al. [21] studied syntactic foams in crash boxes for energy absorption.
Honeycomb cores subjected to crushing were studied by Aktay et al. [22], where methods of modeling honeycomb were validated by experiments, with good agreement. Nayak et al. [23] utilized honeycomb sandwich panels and illustrated their use for blast load mitigation. Additional work on honeycomb was conducted by Sun et al. [24]. The study contributed to increase the knowledge of crashworthiness and collapse modes of honeycomb cores. Honeycomb patterns were also studied by Wu et al. [25], where various honeycomb geometries were investigated at low impact velocities. The core consisted of aluminum whereas carbon reinforced polymer (CFRP) was used as face plates, which was compared to a panel solely based on CFRP. The increased performance of the sandwich, with respect to crashworthiness, was evident. Hammarberg et al. [26] manufactured a UHSS sandwich panel with a bidirectionally corrugated core, intended for stiffness applications. Homogenization of the core was carried out by analyzing a representative volume element under periodic boundary conditions, reducing computational time with maintained accuracy.
In the work by Mohr and Wierzbicki [27], a sandwich based on a perforated, shear-rigid, core was studied numerically for crash box columns. The perforated core was modeled as a continuum, using a honeycomb material model. Detailed modelling, with precise discretization of the perforated core, was not included to check the accuracy of the honeycomb model. The box columns were loaded axially and specific energy absorption was studied for a set of sandwich cores. No increases in specific energy was observed for the numerical models. Furthermore, it was discussed by the authors that skin fracture is likely to occur and must therefore be taken into account when predicting sandwich performance.
The present work suggests a lightweight UHSS sandwich concept with a perforated core, similar to what was studied by Mohr and Wierzbicki [27], for energy absorption applications. Thus, combining the benefits of UHSS and sandwich structures. The concept is investigated numerically by subjecting sandwich hat-profiles to a quasi-static crushing load using the finite element method (FEM). The obtained energy absorption is compared to solid hat-profiles of equivalent weight. Accuracy of the predicted sandwich performance is improved by including an incremental stress state dependent damage model. Failure strain is prescribed dependent on stress triaxiality, allowing for study of skin fracture. Precise discretization of the core, using solid finite elements, causes long computational time. To address this, the core is homogenized using an equivalent material formulation, greatly reducing the number of finite elements and thereby computational time. Input data for the equivalent material is obtained by characterizing a representative volume element (RVE) for in-plane loading and out-of-plane bending/twisting, using periodic boundary conditions. From the RVE analysis equivalent properties for the perforated sections of the core are determined The objectives of the paper is to investigate if the sandwich concept can be utilized for increasing specific energy absorption of hat-profiles, and to study if homogenization can be utilized for reducing computational time while maintaining accuracy. The novelty of the work lies in the study of a UHSS sandwich concept, intended for energy absorption, using refined finite element models with the inclusion of skin fracture, as well as homogenization of the perforated core using a stress-resultant based equivalent representation.
Following the introduction is a section describing the finite element models of the study, where the detailed models are presented as well as the damage model and the homogenization procedure. This is followed by a discussion of the results as well as conclusions drawn.

Finite element modeling
In the present section, the finite element models used for studying the energy absorption of the sandwich hatprofile are introduced in Fig. 1. The geometry and material of the sandwich hat-profiles are presented in Sect. 2.1. A parameter study is also conducted and presented in the section, investigating how the sandwich constituents contribute to the energy absorption. This is followed by a presentation of the detailed finite element models used for analyzing energy absorption, see Sect. 2.2. In the detailed models, core and skins are represented by solid-and shell elements respectively. The performance of sandwiches is limited by their failure modes, such as skin fracture, shear failure of the core and delamination. In the detailed FE models, skin fracture is studied by including an incremental stress-state dependent damage model, calibrated for plane stress, described in Sect. 2.3. Inclusion of a damage model for the skins is necessary since sandwich structures are generally thicker than the component it is replacing, and thus more likely to fail under fracture. By including the damage model it is possible to approximately predict the amount of energy the sandwich can absorb before skin failure, at which point the simulations are terminated. Core failure and delamination are excluded from the FE analysis. The absorbed energy is compared to the performance of a solid hat-profile of equivalent weight and intrusion. The results indicate if the presented sandwich concept is a realistic contender for lightweight energy absorption applications. A large amount of solid finite elements are required for precise discretization of the perforated core in the detailed models, causing long computational times. To address this issue, the perforated sections of the core are homogenized, using a stress-resultant based constitutive model, with equivalent mechanical properties to that of the perforated section of the core. Equivalent material data is obtained by subjecting a representative volume element to periodic boundary conditions from which desired data is extracted, presented in the final part of this section. All simulations are run using the multi-physics solver LS-DYNA R11.0 [28]. Explicit-and implicit time integration are used for the hat-profiles and RVE analysis, respectively.

Geometries and materials
The middle surface of the hat-profile geometry is presented in Fig. 2. Due to symmetry a quarter model is utilized. Dimensions and planes of symmetry are presented in Fig. 2a, b. Thicknesses of core and skins are varied in the study, see Table 1 and Fig. 2c. The core of the sandwich consists of perforated sections, see Fig. 2d, which are only placed over the plane sections and not the radii, forming a partly perforated core. For the perforated sections of the core, the relative core density, * , is defined as the ratio between the densities of the core, core and the bulk material, bulk , according to Equation (1) is rewritten in terms of the hole diameter, D, the spacings d A and d C , and the angle, (see Fig. 2d) where In accordance with Mohr and Wierzbicki [27], the limiting case is when the spacing between the holes approaches zero ( d A = d C = 0 and = 60 • ), corresponding to a relative core density, * , of 9.3% for the perforated sections of the core.
The UHSS sandwich concept is evaluated by subjecting a hat-profile geometry to crushing. A frontal view is presented in a and a perspective view is presented in b. The barrier on top crushes down on the hat profile, while the flange edges are clamped and the support below is fully fixed Two different hole diameters are used for the sandwich: 3 mm and 1.5 mm. The size of the hole diameters are chosen to be relatively small compared to the overall structure of the hat-profile, while still not requiring tremendous amounts of finite elements for geometrical resolution. Three parameters are varied to investigate their influence on the energy absorption and fracture of the sandwich with a hole diameter of 3 mm. The parameters studied are the hole distribution, d A , core/skin distribution with maintained initial thickness, and reduced total thickness of the sandwich while maintaining initial distribution between core and skins. The last parameter is also investigated for a hole diameter of 1.5 mm. The angle, , is kept constant at 45 • and d C is determined from Eq. (3). A summary of the tested parameters is presented in Table 1, where it should be noted that the mass presented corresponds to a quarter of the mass of the hat profile according to the symmetry conditions used in the numerical models. The entire sandwich is based on the fully hardened UHSS, 22MnB5. Young's modulus and density is 200 GPa and 7800 kg∕m 3 respectively and Poison's ratio is 0.3.

Detailed finite element models
Hat-profiles, subjected to crushing, are used to numerically study the energy absorption of the UHSS sandwich. A quarter finite element model of a hat profile is presented in Fig. 3. Planes of symmetry are indicated in Fig. 2a and To ensure quasi-static conditions, i.e. inertia forces and kinetic energy may be neglected, sandwich #1 of Table 1 is analyzed by subjecting it to the three velocity curves of Fig. 7a. By comparing the quotient between kineticand internal energy, as well as the force-displacement response, quasi-static conditions can be ensured.
The sandwich core of the hat-profile is discretized using fully integrated 8-node solid elements, whereas plane stress shell elements are used for the skins. During crushing the hat-profile is subjected to bending. To resolve the bending mode and ensure mesh independent results a convergence check is conducted regarding the required number of solid elements through the thickness. This is done by studying a sandwich consisting of solid finite elements. Brick elements with equal side lengths are utilized, starting with one element through the thickness. For each subsequent mesh refinement, the number of elements through the thickness is doubled. Evaluation of mesh convergence is done by comparing internal energies and force-displacement response, computational time is also considered.

Modeling of fracture
To account for skin fracture, a generalized incremental stress-state dependent damage model (GISSMO) is utilized [28]. The fracture model allows for an incremental formulation of the damage accumulation, D , according to In Eq. (4), D corresponds to the damage value in the interval 0 ≤ D ≤ 1 , where elements are deleted if D = 1 . is a damage parameter, p is the incremental plastic strain, f is the equivalent plastic failure strain of Fig. 5. The data of the figure is obtained by calibration against tensile specimens according to [29] using LS-OPT. The equivalent plastic failure strain can be given as a function of both stress triaxiality, , and the Lode parameter, L, determined as where ij corresponds to the Cauchy stress tensor, J 2 to the second deviatoric stress invariant, to the Lode angle (counterclockwise angle from the first principal stress axis in the deviatoric plane) and J 3 to the third deviatoric stress invariant. To account for instability and localization, the instability measure, F, is defined in an incremental formulation according to where F is in the range of 0 ≤ F ≤ 1 . p,loc is the equivalent plastic instability strain of Fig. 5, and similarly to f it can be given as a function of triaxiality and the Lode parameter. Once F has reached unity, the current value of D is stored, and from this point coupled to the flow stress.
Calibration of the fracture model can be performed by using the force-displacement response from the five specimens presented in [29]. The curves for plastic failure strain, f , and the plastic instability strain, p,loc , were determined by inverse modeling using LS-OPT [30].

Homogenized finite element model
Using solid elements for detailed discretization of the sandwich core results in long computational times. To address this issue, the perforated regions of the core (see Fig. 3b) are homogenized using a stress-resultant based constitutive model, implemented for shell elements. Modeling of the solid regions of the core is done by the piece-wise linear plasticity model presented in Sect. 2.2. The stress-resultant model treats in-plane behavior separately from bending, typically necessary for perforated materials where differences in the elastic constants, such as Young's modulus, may differ greatly for in-plane loading and bending. The in-plane, plane stress elastic tensor, C p , is defined as where out of plane shear is accounted for by Q 55 and Q 66 . The constituents of the matrix are given by  (8) and (10), as defined in laminate theory [31], in the following manner where N and M are the membrane forces and bending moments, whereas and are the in-plane strains and curvatures. In the present context B = 0 and A and D are defined as where t corresponds to the thickness of the laminate. The elastic constants of Eqs. (8) and (10), as well as the yield curves for the resultants, are determined by subjecting the representative volume element (RVE) in Fig. 6, to the corresponding stress states. In-plane properties are obtained using periodic boundary conditions (PBC). PBCs are necessary for periodicity to be maintained and to not overpredict stiffness when characterizing a single RVE, see e.g. Hammarberg et al. [26], Suquet [32], Xia et al. [33], Omairey et al. [34]. A short description of PBC is given in Sect. 2.5.

Periodic boundary conditions
The local strain, ik , is rewritten as where ̄i k and ik (u * ) are the average and fluctuating parts of the strain, and u * is a periodic displacement field. The where i and j relate to the degree of freedom (DOF) and the face connected to the node, and i, j = (1, 2, 3) for a three dimensional RVE.
In the present work, LS-DYNA is used to apply PBC on the RVE. The constant, c j i , in Eq. (20) corresponds to the displacement of dummy nodes, generated by the software, from which displacements and forces can be obtained. From the dummy nodes, the homogenized quantities for stress and strain are obtained. Initially, the deformation gradient, F ij , is determined from the displacements of the dummy nodes as (16) x j + ij (no summation over indices) .
ij and x i are the Kronecker delta and the dimension of the RVE respectively. Forces, f ij , are extracted from the dummy nodes from which the macroscopic first Piola-Kirchhoff stress tensor, P ij , is determined. The Green-Lagrange strain tensor, E, is obtained from the deformation gradient according to The Cauchy stress tensor is obtained as where J is the determinant of the deformation gradient, i.e. the Jacobian. From Eqs. (22) and (23), the homogenized constitutive relationship for various loading conditions can be determined.

Results and discussion
Several finite element simulations were utilized to investigate the energy absorption of the ultra high strength sandwich concept. Initially, a quasi-static check and a mesh convergence study were conducted. In the subsequent simulations, energy absorption of the sandwich and its constituents were studied for several sandwich layouts (see Table 1).

Quasi-static check and mesh convergence
To ensure quasi-static conditions were met, sandwich #1 of Table 1, was subjected to the three velocity curves of Fig. 7a. Comparing the quotient between kinetic and internal energy, see Fig. 7b, and the force-displacement curves of Fig. 7c, where the accumulated internal energy is presented in the legend, indicated that velocity curve #01 ensured quasi-static conditions. The difference between velocity curve #01 and #02, in terms of accumulated internal energy, was ≈ 1.8% and since the force-displacement The results obtained from the mesh convergence study are presented in Fig. 8a, where internal energy together with CPU time is presented for each level of mesh refinement. The corresponding force-displacement response is presented in Fig. 8b. It is noted that the model with two solid elements through the thickness produces adequate results with respect to the internal energy as only small differences are observed for the subsequent mesh refinements. When comparing the force-displacement responses it is clear that four elements through the thickness is necessary to get a converged response over the entire range of displacement. Thus four solid elements through the thickness was chosen, providing a good trade off between accuracy and CPU time.

Energy absorption of sandwiches
Eleven sandwiches (Table 1) were evaluated numerically for energy absorption in order answer the question if the suggested UHSS sandwich concept is a contender for increasing specific energy absorption, as compared to solid hat profiles of equivalent weight. In the models, Fig. 7 The response obtained for the three tool velocity curves in a are investigated with respect to the kinetic energy in b and the force response in c to ensure negliable dynamic effects. The results are presented for Sandwich 1 in Table 1 Fig. 8 For the mesh convergence study four levels of mesh refinements are used: 1, 2, 4 and 8 fully integrated cubic solid elements through the thickness. In a internal energy and computational time are are presented versus number of elements. For each mesh size the force-displacement response is presented in b skin fracture was included by a damage model calibrated for shell elements under plane stress conditions. Results are presented until the point of sandwich skin fracture, comparing both sandwiches and equivalent models for a given intrusion. To illustrate the deformation mode, the initial and final state are presented for sandwich #1 and #3 in Fig. 9.
In Fig. 10a, sandwiches #1-3, with hole diameters of 3 mm, are represented by the solid lines. Each sandwich was compared to a solid hat profile of equivalent weight, represented by the dashed lines in the figure. For each sandwich, the results are presented up until the point of skin fracture, at which point the curves are cut off. For sandwich #1, the accumulated damage just before fracture is presented in Fig. 11a and after fracture in Fig. 11b. In the figure, the concentration of damage occurs near the top rows of holes, corresponding to the concentration of effective plastic strain, see Fig. 12a.
To evaluate the impact of the hole diameter on the sandwich response, a diameter of 1.5 mm was also evaluated, corresponding to sandwich #8-10. The force-displacement response is presented in Fig. 10b, cut off at the point of skin fracture. The accumulated damage is presented before and after skin fracture in Fig. 11c, d for sandwich #8. As for a hole diameter of 3 mm, the concentration of damage is located near the top row of the holes with some additional concentrations at the top radius. However, the main part of the fracture for sandwich #8 still occurs in a similar location as for sandwich #1.
From the results presented in Fig. 10, it was noted that by utilizing the sandwich concept, the specific energy absorption was increased by approximately 8-17%, as compared to the solid hat profiles. A main difference observed was the impact of the hole diameter on the sandwich ductility. For all the sandwiches, the smaller hole diameter allowed for a greater displacement before fracture, but the effect is particularly clear for sandwich #1  . 9 The initial and final states are presented for sandwich #1 and #3 to illustrate the loading case and deformation mode used for evaluating the energy absorption of the sandwich concept and #8 where the difference in allowable displacement is approximately 50% in favor of the latter. This also emphasises the importance of including fracture models for evaluations of sandwich concepts, as the point of fracture is greatly affected by the layout of the sandwich. Thus, it would be beneficial to also include a fracture model calibrated for solid elements, making it possible to predict core fracture as well. In this work, the fracture model had been calibrated for plane stress and no data was available for solid elements which is why only skin fracture was studied.
To evaluate the contribution by the sandwich constituents, core and skins, the hole distribution was varied while the thickness of the layers were held constant. For sandwich #4 and #5 the hole spacing was increased and decreased by 25% respectively. The force-displacement response is presented in Fig. 13a, with sandwich #1 as a reference. Additionally, for sandwich #6 and #7, the hole distributions and total thicknesses were held constant while the layer thicknesses were varied. The force-displacement response is presented in Fig. 13b. The intention with comparing sandwiches #1 and #4-7 was to evaluate how the layout of the sandwich affected the energy absorption. From the results it was observed that adding mass to the sandwich increased its energy absorption, but more importantly, it was noted that the placement of the mass mattered. The further away from the neutral layer that mass was added the greater was the effect on specific energy absorption. When sandwich #4 and #5 were compared, the difference in mass was found to be 7.4% whereas the difference in SEA was only 5.9%. A similar comparison for sandwich #6 and #7 indicated an increase of SEA of 30% whereas the weight had only been increased by 8%. A comparison between sandwich #5 and #6 was also done. Sandwich #5 with a mass of 18.9 g exhibited a lower SEA compared to sandwich #6 with a mass of 18.6 g. The greater SEA was attributed to the fact that sandwich #6 had a larger portion of its mass further away from the neutral layer, where the bending stress is zero, due to its core thickness. The distribution of the internal energy of the sandwich constituents is presented in Fig. 14. It was noted that the major difference between the two sandwiches was the contribution by the lower skin of sandwich #6, which had a significantly higher energy absorption compared to the other layers. A summary of the performance of the sandwiches in terms of SEA is presented in Table 2.
The weight saving potential of the sandwich was also evaluated. This was done by comparing the weights of a solid and a sandwich hat-profile with similar energy absorption. The obtained results are presented in Fig. 15, from which it was found that the weight could be lowered by approximately 7%.
An effort was made to reduce computational time by replacing the perforated sections of the core with an equivalent material. The equivalent material was modeled using a constitutive routine based on a stress-resultant formulation. Two sets of elastic constants were required, corresponding to in-plane loading and bending respectively, presented in Table 3. Plasticity was modeled by specifying uncoupled yield curves for the resultants, see Fig. 16a, b. For the elastic constants differences were observed for both Young's modulus and Poisson's ratio. This indicated the necessity for the constitutive routine to make a distinction between in-plane loading and bending, motivating the chosen material model. From the results presented in  Fig. 16c, it was found that the homogenized model captured the overall force-displacement response quite well while reducing computational time with approximately 95%. The energy absorption of the homogenized model corresponds well to the detailed model, both absorbing 112.4 J. However, the homogenized model was not able to predict fracture corresponding to what was observed in the detailed models. In the detailed models fracture consistently occurred due to the holes of the core, whereas the homogenized predicts failure at a much later stage occurring over the top radius in the skin.
The results of the present study should be compared to the results obtained by Mohr and Wierzbicki [27], where a similar sandwich concept was utilized for crash boxes.
According to Mohr and Wierzbicki [27] the specific energy absorption remained constant regardless of the distribution of core and skins, and redistribution of mass away from the neutral layer did contribute to increasing the SEA. A major difference between [27] and the present study is the usage of detailed FE models. In the work done by Mohr and Wierzbicki [27] homogenized models, based on a honeycomb model available within LS-DYNA [28], was utilized for representing the core, where as the present work used detailed FE models for analyzing the sandwich concept.
Debonding is an important failure mode for laminates and skins, which was not studied in the present work but would be necessary for better predictions of the sandwich performance. Fig. 11 The accumulated damage parameter before and after fracture is presented for sandwich 1 in a, b and for sandwich 8 in c, d Fig. 12 Effective plastic strain before fracture for sandwich 1 in a and sandwich 8 in b Fig. 13 Comparing forcedisplacement response and energy absorption when varying hole distribution and core and skin thickness. The total thickness and hole diameter are kept constant at 1.848 mm and 3 mm respectively The present work highlights the importance of accurate representation of the perforated core, with respect to the distribution of mass around the neutral layers. Alternatively, a constitutive model accounting for the phenomena observed in the present study should be utilized. A constitutive model, developed for perforated sheets, should also be considered as it could drastically reduce computational time.

Conclusions
The present study investigated the possibility of utilizing ultra-high strength steel in sandwich structures for energy absorption applications. It was found that the suggested sandwich concept provides a means for increasing specific energy absorption of components, and a weight saving potential for equal performance. The models used can be a tool for lightweight component development for energy absorption applications. Distribution of mass was shown to be of importance, since displacing the mass away from the neutral layers increased specific energy absorption. Inclusion of fracture models are required for correctly predicting sandwich performance, since skin fracture occurred earlier in the sandwich hat-profiles as compared Fig. 15 To illustrate the weight saving potential of the sandwich concept, a solid hat is compared to a sandwich with equivalent energy absorption  to the solid hat-profiles. Furthermore, constitutive models must distinguish between in-plane loading and bending due to the differences in elastic constants for each loading cases. Homogenization of the perforated core can be utilized for reduction of computational time with maintained accuracy. A limit of the present study, is that a perfect bond was assumed between the layers of the sandwiches and delamination was not studied. In summary, methods for reducing weight, by introducing light and strong materials and structures, are of importance for meeting demands regarding reduction of emissions from the automotive industry. The numerical models of the suggested UHSS sandwich concept provides a way of increasing specific energy absorption. This indicates that a weight saving potential exists by utilizing sandwich structures for energy absorption applications.
In summary, methods for reducing weight, by introducing light and strong materials and structures, are of importance for meeting demands regarding reduction of emissions from the automotive industry. The numerical models of the suggested UHSS sandwich concept provides a way of increasing specific energy absorption. This indicates that a weight saving potential exists by utilizing sandwich structures for energy absorption applications.
Acknowledgements The authors would like to thank Mikael Schill at DYNAmore Nordic for the calibrated GISSMO data used in the present study. Economic support is supplied through the Swedish lightweight innovation programme-LIGHTer, which is gratefully acknowledged.

Compliance with ethical standards
Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.