1. Introduction

Figure 2. Thermal runaway mechanism of lithium-ion battery for EVs [2].

Figure 2. Thermal runaway mechanism of lithium-ion battery for EVs [2].

Figure 3. Diagram of side collision test with movable barrier according to ECE Regulation R95 [5].

Figure 3. Diagram of side collision test with movable barrier according to ECE Regulation R95 [5].

2. Materials and Methods

2.1. Methodology

We evaluated the safety of different electric vehicle configurations, by means of by correlating the simulation results with the thresholds set out in ECE Regulation R95. The reference car model was a combustion engine car (Geo Metro 3-Door Gen II), modeled with finite elements, including a mobile barrier. The side impact was simulated according to Directive 97/26/EC.

Figure 4. (a) Geo Metro 3-Door Gen II model; (b) Geo Metro 3 model modeled with finite elements. Model used as a reference.

Figure 4. (a) Geo Metro 3-Door Gen II model; (b) Geo Metro 3 model modeled with finite elements. Model used as a reference.

This model was converted to an electric vehicle by introducing two different traction battery configurations; subsequently, some absorbers were applied. A real battery of a first-generation Nissan Leaf was taken as a reference (48 modules 24 kWh), containing 48 modules with 4 cells each. The materials applied in this real battery were via measurements from some samples from the battery.

Figure 5. (a) Vehicle floor and area under rear bench of the reference vehicle; (b) adapted electric vehicle floor and area under rear bench to create room for the traction battery.

Figure 5. (a) Vehicle floor and area under rear bench of the reference vehicle; (b) adapted electric vehicle floor and area under rear bench to create room for the traction battery.

We used a battery with 28 modules with 4 pouch cells each. Some cells with higher energy density, such as prismatic and cylindrical cells, are safer [8]. However, pouch cells were selected because they are commonly used in most electric vehicles. The least favorable scenario was studied. In total, four different electric vehicle configurations were analyzed; each of them is explained bellow.

2.1.1. Electric Vehicle Configuration 1

A high-voltage battery was designed based on maximizing the use of the available space on the floor and under the rear bench and minimizing the increase in floor height under occupants’ feet. The separation between the modules and the housing wall was 50 mm, which is the recommended dimension to avoid damage to cells in case of impact [9]. The high-voltage battery housing was modeled in two parts: the bottom and top housing, reproducing a real battery scheme. The modules were characterized according to the technical data of the Nissan Leaf battery [10]. The battery dimensions were 303x223x55 mm. The inner volume of the modules was modeled with these measures as a solid. The battery housing was represented with a 0,5 mm thickness and the same overall dimensions. Moreover, a 245x135x30 mm space was reserved on one side of the modules to represent the space occupied by the battery management system (BMS). However, the BMS was not implemented in the finite element model. This did not interfere with the results obtained from the simulation because the BMS is located on the side opposite from the impact.

2.1.2. Electric Vehicle Configuration 2

In this second configuration, the results of a study on the propagation of thermal runaway were considered [11]. The authors of this study stated that the possibility of the propagation of thermal runaway is higher in the vertical than in the horizontal arrangement. In addition, the importance of the spaces between packages and their cooling to ensure safe operation was highlighted, which reduce the likelihood of a loss of heat transfer control in the battery assembly. In this configuration, protection against thermal runaway is improved because the modules on the floor of the vehicle are reorganized, and the overall width of the assembly is reduced, leaving more space for possible lateral shock absorbers. In this configuration, 90 mm of space exists at both battery set sides, which is larger than the 50 mm in configuration 1.

The method of modeling the battery housing was the same as that in configuration 1, with the dimensions adapted to the new battery layout. The modules, their respective housings, and the BMS had the same dimensions as in configuration 1.

Figure 6. High-voltage battery layout in EV (a) configuration 1; (b) configuration 2.

Figure 6. High-voltage battery layout in EV (a) configuration 1; (b) configuration 2.

2.1.3. Electric Vehicle Configuration 3

With both configurations 3 and 4, we aimed to improve the battery protection in the event of a collision by implementing shock-absorbing elements on both sides of the housing, adapting their dimensions to the battery housing shape and the available space.

Based on prior results [12], in which three alternatives for different shock absorbers were compared, the material chosen for these new components was aluminum foam (Table 1).

Configuration 3 corresponds to configuration 1 with the addition of absorbers. In configuration 1, less space was available on the sides of the battery, as the placement was lower on the vehicle floor and was therefore wider. In this case, the maximum absorber width was 10.9 cm.

Figure 7. (a) Configuration 3: configuration 1 with shock absorbers; (b) configuration 4: configuration 2 with shock absorbers.

Figure 7. (a) Configuration 3: configuration 1 with shock absorbers; (b) configuration 4: configuration 2 with shock absorbers.

2.1.4. Electric Vehicle Configuration 4

Configuration 4 was configuration 2 to which shock absorbers were added. In configuration 2, the battery distribution was higher above the vehicle floor and allowed for a narrower housing. Therefore, the shock absorbers had a maximum width of 14 cm (3.1 cm wider than in configuration 3).

Figure 8. Comparison of four electric vehicle configurations analyzed.

Figure 8. Comparison of four electric vehicle configurations analyzed.

2.2. Characterization of Battery Materials

For the finite element modeling of a traction battery to conduct simulations with various computer programs, we needed characterize the battery materials . To this end, Vickers microhardness tests were carried out on a sample of both the upper and lower cover of the original battery casing and the module casings to determine their mechanical properties for their introduction into the simulation program. On the other hand, to characterize the material inside the modules, we determined the composition of the four cells of the bag that composed the modules.

The Vickers test is based on an optical measuring system. A diamond indenter shaped like a square pyramid is used to perform the microhardness test. The indenter is used to penetrate the material studied, with a known load and for a period that is usually between 10 and 15 seconds. Subsequently, the length of the diagonal of the indentation left by the diamond indenter is measured with a microscope to calculate the area of the inclined surface of the indentation. With these data, the Vickers hardness can be calculated as the quotient between the load applied during the test and the indentation area.

We collaborated with the Materials Department of the University of Zaragoza to complete this test. First, specimens with the necessary dimensions for the tests were obtained.

Figure 9. Explanation of Vickers microhardness test.

Figure 9. Explanation of Vickers microhardness test.

Figure 10. Test specimen obtained from the (a) upper and (b) lower casing of the high-voltage battery, and from the (c) casing of a high-voltage battery module.

Figure 10. Test specimen obtained from the (a) upper and (b) lower casing of the high-voltage battery, and from the (c) casing of a high-voltage battery module.

The conducted tests are shown below.

Figure 11. Views of Vickers test on specimen from upper casing of high-voltage battery.

Figure 11. Views of Vickers test on specimen from upper casing of high-voltage battery.

Figure 12. Views of Vickers test on high-voltage battery module housing specimen.

Figure 12. Views of Vickers test on high-voltage battery module housing specimen.

The results obtained in each of the tests (two per specimen) are shown in Table 2.

Based on the results of the microhardness tests on both the upper and lower casing covers, a hardness of 110 kg/mm² was obtained. According to the interpretation in the conversion tables [13], a mechanical strength of 53 psi was deduced. Applying the conversion factor (6.9 N/mm² per psi unit), we obtained a mechanical strength of 365.7 N/mm². The yield strength can be also approximated as three times the Vickers hardness, i.e., 330 MPa. With these results, we compared the various materials that can be used to manufacture the casing, among which the low-carbon steels DC01, DC02, and DC03 were compared [14]. We chose hardened DC01, which has a yield strength of 200-380 N/mm² and a mechanical strength of 290-430 N/mm². The properties used to define the steel chosen are listed in Table 3 and Table 4.

The same procedure was applied for the individual module housings, producing a Vickers hardness of 71 kg/mm², a yield strength of 213 MPa (three times the Vickers hardness), and a mechanical strength of 227.7 MPa. Several aluminum alloys and low carbon steels were compared, and the aluminum alloy 6063-T6 was chosen, which had a yield strength of 215 MPa and a mechanical strength of 240 MPa. The properties used to define the chosen Aluminum 6063-T6 are listed in Table 5 and Table 6.

The content of the modules was defined by calculating the equivalent properties from the components of the four pouch cells inside them. These modules had a graphite anode and a LiMn₂O₄ cathode with LiNiO₂ [10]. Both separator foils and cell envelopes were neglected. LiMn₂O₄ and LiNiO₂ were assumed to have the same characteristics, and the properties of the anode and cathode were weighted with a coefficient of 0.5 each. Based on previously collected data [16], the properties of the equivalent module material were obtained, and they are summarized in Table 7.

2.3. Occupant Injury Analysis

This vehicle model was imported into MADYMO software (MAthematical DYnamic MOdel). MADYMO is software that allows the analysis of the injuries suffered by the occupants in a collision and to determine their origin. It is used to determine the behavior of occupant restraint systems (seat belts, airbags, and head restraints).

In this study, MADYMO was used to analyze the movement and stresses experienced by a right front occupant in a vehicle involved in a collision on the right side, with reference to the running direction. For this purpose, the vehicle was fitted with a seat on which the dummy was located. The dummy movement and the stresses it experienced because of the lateral collision were analyzed. The seat belt was also modeled.

2.3.1. Dummy Model

We conducted a parametric evaluation of the effect of a side impact on the kinematics of a dummy. For this purpose, a EuroSID 2 was used for each of the five vehicle configurations analyzed. The test conditions were those established in the ECE Regulation R95, in which the ES-2 dummy (EuroSID 2) is applied. The models used by MADYMO are shown below (Figure 13). They were calibrated and validated using numerous impact tests on components and complete dummies.

2.3.2. Vehicle Model

For the analysis, the car model was imported (for the five configurations), and the seat was fitted. In addition, the stresses experienced by the right front occupant, with reference to the running direction, were analyzed.

Figure 14. Image of the vehicle used, and details of the seat fitted for the dummy.

Figure 14. Image of the vehicle used, and details of the seat fitted for the dummy.

2.3.3. Lateral Crash Test with Deformable Barrier

The side crash test with a deformable barrier was analyzed according to the UN ECE R95 standard. This test simulates a side crash between two cars at a given relative angle, as might occur at a road junction, for example. To perform this test, a vehicle to be tested, a movable deformable barrier, and dummies are required. A deformable barrier with the characteristics of the barrier used in the ECE Regulation R95 was added to the above model.

In FMVSS 214, the ES-2 re dummy is used; in UN R95, the ES-2 dummy is used; and in EuroNCAP, the WorldSID dummy is used. For this reason, the ES-2 dummy was applied following ECE Regulation R95.

Figure 15 shows the whole model that was simulated including the corresponding dummy and the seat belt.

3. Results and Discussion

Table 8. Comparison of instant of maximum deformation of car model.

Figure 17. Average velocity (m/s) measured on vehicle ground during impact for each configuration.

Figure 17. Average velocity (m/s) measured on vehicle ground during impact for each configuration.

Figure 18. Average acceleration (m/s²) measured at vehicle floor during impact for each configuration.

Figure 18. Average acceleration (m/s²) measured at vehicle floor during impact for each configuration.

Table 10. Comparison of protection elements.

3.2. Intrusion into Passenger Compartment of the Vehicle

Intrusion into the passenger compartment was measured with four different parameters, measured at the instant of maximum deformation in each of the cases analyzed, to determine the effect of the battery on the structural behavior of the vehicle during a side impact.

3.2.1. Intrusion at Side Sill

The first two parameters considered were the maximum distance and the average distance that different nodes penetrate at side sill height (Figure 19). These distances were calculated as the difference between the vehicle width before and after impact at different nodes.

We observed that both the battery and the absorbers stiffened the vehicle structure, as the intrusions were reduced. By including the battery, the floor was stiffened, and the side sill intrusions were considerably reduced. The intrusions were smaller in the case of configuration 1 as the battery was wider than in the case of configuration 2; therefore, a larger zone of the vehicle floor was stiffened.

By implementing shock absorbers, the average penetration into the passenger compartment at the side sill was further reduced. However, in the second configuration with absorbers, the maximum intrusion obtained between the nodes was almost the same as that without absorbers, whereas the average intrusion considerably decreased. This means that one of the nodes had a large deviation from the line of nodes; however, generally, the right side of the vehicle reduced the intrusion of its surface.

Table 12. Comparison of protection elements.

Table 12. Comparison of protection elements.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
dTmax(mm)	333,55	224,85	253,02	205,55	253,48
dTmax(mm)	252,60	128,862	160,31	93,87	100,67

3.2.2. Intrusion at B Pillar

For this element, two different parameters were calculated to quantify the intrusion into the passenger compartment at the B pillar of the vehicle at the instant of maximum deformation.

The first, 𝑑_𝐻, is defined as the distance that a B-pillar node penetrates inside the passenger compartment with respect to the two nodes of the A and C pillars. The distances from this node to the imaginary line formed by two nodes on the A and C pillars of the vehicle at the same z-value, both before and after the impact (Figure 20), were calculated and summed to calculate the value of 𝑑_𝐻.

The second parameters, 𝑑_𝑉, is the distance penetrating the center of the B pillar in relation to its anchorages. It was calculated as the sum of the distances before and after the impact (Figure 21) from a node in the middle of the pillar to the imaginary line between two nodes each located at an anchorage of the pillar.

When analyzing the intrusions on the B- pillar of the vehicle by including the battery assembly, a reduction in d_H was achieved. After implementing the absorbers, this distance was further reduced. d_V slightly increased when the battery pack was fitted, as the floor did not deform as much as the rest of the side, and this intrusion was measured as a function of the distance from a node in the center of the pillar to the line formed by its anchorages. However, this distance decreased with the implementation of the absorbers.

Table 13. Intrusion at B pillar for the different configurations.

Table 13. Intrusion at B pillar for the different configurations.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
dH(mm)	351,48	222,91	238,96	160,19	120,95
dV(mm)	101,32	103,71	108,85	85,30	81,24

3.3. Results for Different Electric Vehicle Configurations

The movement of the dummy and the stresses to which it was subjected were analyzed for the different vehicle configurations. The injury thresholds indicated by the regulations and the private programs that analyze vehicle safety were taken as the reference values. The test configuration used was ECE Regulation R95.

Below is a sequence of frames showing the movement experienced by the passenger of a vehicle subjected to an acceleration similar to that experienced in a lateral collision according to ECE Regulation R95. The passenger is located on the right-hand side of the vehicle related to the running direction.

3.3.1. Passenger Kinematics in a Combustion Engine Vehicle

Figure 22. Kinematics of passenger in combustion vehicle subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

Figure 22. Kinematics of passenger in combustion vehicle subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

3.3.2. Passenger Kinematics in an Electric Vehicle: Configuration 1

Figure 23. Kinematics of passenger for electric vehicle configuration 1 subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

Figure 23. Kinematics of passenger for electric vehicle configuration 1 subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

3.3.3. Passenger Kinematics in Electric Vehicle: Configuration 2

Figure 24. Kinematics of passenger for electric vehicle configuration 2 subjected to lateral collision according to ECE Regulation R95, using the ES-2 dummy.

Figure 24. Kinematics of passenger for electric vehicle configuration 2 subjected to lateral collision according to ECE Regulation R95, using the ES-2 dummy.

3.3.4. Passenger Kinematics in Electric Vehicle: Configuration 1 with Energy Absorbers

Figure 25. Kinematics of passenger for electric vehicle configuration 1 with energy absorbers, subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

Figure 25. Kinematics of passenger for electric vehicle configuration 1 with energy absorbers, subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

3.3.5. Passenger Kinematics in Electric Vehicle: Configuration 2 with Energy Absorbers

Figure 26. Kinematics of passenger for electric vehicle configuration 2 with energy absorbers, subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

Figure 26. Kinematics of passenger for electric vehicle configuration 2 with energy absorbers, subjected to a lateral collision according to ECE Regulation R95, using the ES-2 dummy.

A comparison of the movement showed that with a combustion vehicle, the deformation of the lower part of the vehicle is increased, and the occupant moves more to the side opposite to the impact side. For configurations 1 and 2 with energy absorbers, the occupant moves more toward the side opposite to the impact side than in the cases of the configurations without absorbers. This occurs because the traction battery with absorbers stiffens the area and, when shock absorbers are not included, has a small clearance that allows it to deform and not directly receive the impact. With absorbers, no space exists before the traction battery plus absorbers assembly is impacted.

Figure 27. Passenger kinematic comparison for different vehicle configurations.

Figure 27. Passenger kinematic comparison for different vehicle configurations.

3.3. Influence of Battery Layout on Risk of Injury

The influence of battery layout and the presence or absence of energy absorbers on the risk of injury to the passenger is explained below.

3.2.1. Analysis of Head Injuries

According to ECE Regulation R95, the head performance criterion (CCC) applies when contact occurs with the head in a crash. The CCC is the maximum value of the following formula (1):

$$({\mathrm{t}}_2-t_1){\left(\frac{1}{t_2-t_1}{\int }_{t1}^{t2}adt\right)}^{2.5}$$

(1)

where a is the resulting acceleration at the center of gravity of the head, in meters per second squared, divided by 9,81, measured as a function of time and filtered with a channel frequency class of 1 000 Hz; t1 and t2 are any two moments between the initial and final contact.

This formula corresponds to the head injury criterion (HIC) (2):

$$\mathrm{HIC}=\underset{T0\le t_1\le t_2\le \mathit{TE}}{\max }{\left[\frac{1}{t_2-t_1}{\int }_{t1}^{t2}R\left(t\right)dt\right]}^{2.5}\left(t_2-t_1\right)$$

(2)

where t1 is the start time, t2 is the end time, and R(t) is a(t), the acceleration curve experienced by the struck head, where (t) is measured in seconds and a is measured in g's.

A value of 1000 is cited as the injury threshold in aviation regulations. First, the HIC calculation was limited to a 36 ms interval (HIC₃₆) with an injury threshold of 1000. Subsequently, it was revised, limiting the maximum time interval to 15 ms and reducing the injury threshold to 700 (referred to as HIC₁₅). An HIC₁₅ value of 700 represents a 5% risk of serious injury or an abbreviated injury scale (AIS) of 4.

The CCC should be less or equal to 1000 when no head contact occurs. In that case, the CCC is not measured or calculated. Instead, "no head contact" is indicated.

Table 14 shows the HIC36 value for each of the configurations analyzed.

The values are well above the injury threshold (1000), which is why the curtain airbag and the side airbag have been included in vehicles.

To explain the difference among the HIC in the different configurations, the acceleration experienced by the head of the dummy in each configuration was analyzed.

Figure 28 shows how the maximum acceleration value experienced by the occupant of a combustion vehicle is reduced when a traction battery is added. Moreover, this maximum deceleration is reached later. The lowest acceleration value is reached with configuration 1. However, by adding crash absorbers and reducing the space between the running board (where the deformable moving barrier is impacted) and the traction battery plus absorbers, the maximum value is reached earlier. However, when that area is stiffened, the maximum head acceleration value is higher.

3.2.2. Analysis of Thorax Injuries

According to ECE Regulation R95, the maximum chest deformation is the maximum deformation value for any rib as determined by chest displacement transducers, filtered with a channel frequency of 180 Hz.

The viscosity criterion is applied. The maximum viscosity is the maximum value for the viscosity criterion at any rib, calculated from the instantaneous product of the relative compression of the thorax with respect to the semithorax and the compression speed derived by compression differentiation, filtered with a channel frequency of 180 Hz. The normalized width of the half-thoracic cage is considered as 140 mm, where D (in meters) is the rib deformation:

$$VC=max\left(\frac{D}{\mathrm{0,14}}\bullet \frac{dD}{dt}\right)$$

(3)

The performance criteria for the thorax are as follows:

(a) In the case of the rib deformation criterion, less than or equal to 42 mm.

(b) In the case of the soft tissue criterion (viscosity criterion (VC), less than or equal to 1,0 m/s.

During a transitional period of two years from the date specified in the Regulation, the viscosity criterion value should not be a determining value for passing the approval test but should be recorded in the test report and registered by the approval authority. After the transitional period has elapsed, a viscosity criterion value of 1,0 m/s shall be applied as a pass.

The following chest injury criteria values were obtained in the upper, middle, and lower rib area with MADYMO software:

Table 15. Rib compression and viscosity criterion (VC) for each of analyzed configurations.

Table 15. Rib compression and viscosity criterion (VC) for each of analyzed configurations.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
Upper Rib Compr. (mm)	59,47	52,24	55,08	75,90	67,26
Mid Rib Compr. (mm)	69,51	46,13	48,39	48,82	42,25
Lower Rib Comp. (mm)	68,70	56,66	64,15	38,37	30,18
Upper Rib VC (m/s)	1,57	1,38	1,33	3,71	2,61
Mid Rib VC (m/s)	1,94	1,15	1,34	1,46	1,03
Lower Rib VC (m/s)	2,14	1,78	2,43	1,06	0,57

VC value in the upper part of the ribs is the lowest for configuration 2 and the highest value for configuration 1 with shock absorbers. In the middle part of the ribs, the VC value is the lowest for configuration 2 with shock absorbers and the highest in of the combustion vehicle. In the lower part of the ribs, the VC value is the lowest for configuration 2 with shock elements and the highest for configuration 2.

Figure 29. VC at upper ribs for each of the configurations analyzed.

Figure 29. VC at upper ribs for each of the configurations analyzed.

Figure 30. VC at middle ribs for each of the configurations analyzed.

Figure 30. VC at middle ribs for each of the configurations analyzed.

Figure 31. VC at lower ribs for each of the configurations analyzed.

Figure 31. VC at lower ribs for each of the configurations analyzed.

3.2.3. Analysis of Abdomen Injuries

According to ECE Regulation R95, the maximum force on the abdomen is the maximum value of the sum of the three forces measured by transducers mounted 39 mm below the surface of the impacted side, with a CFC of 600 Hz.

The abdominal performance criterion according to ECE Regulation R95 states that the maximum force on the abdomen must be less than or equal to 2,5 kN of internal force (equivalent to an external force of 4,5 kN).

Table 16. APF value (maximum force in the abdomen) for each of the configurations analyzed.

Table 16. APF value (maximum force in the abdomen) for each of the configurations analyzed.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
APF (N)	9.168,0	1.277,0	1.539,9	1.707,0	1.558,9

We observed that the inclusion of the traction battery considerably reduces the value of the maximum force in the abdomen (APF), reaching values that are below the injury threshold value (2.5 kN). The minimum value is reached with configuration 1.

Figure 32. APF value for each of the configurations analyzed.

Figure 32. APF value for each of the configurations analyzed.

3.2.4. Analysis of Pelvis Injuries

According to ECE Regulation R95, the maximum force on the pubic symphysis is the maximum force measured with a load cell on the pubic symphysis of the pelvis, filtered with a channel frequency of 600 Hz.

The criterion for the pelvic behavior according to ECE Regulation R95 states the maximum force on the pubic symphysis must be less than or equal to 6 kN.

Table 17. Maximum value of maximum force on pubic symphysis (PSPF) for each of the configurations analyzed.

Table 17. Maximum value of maximum force on pubic symphysis (PSPF) for each of the configurations analyzed.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
PSPF (N)	5792,11	4.609,2	4.864,3	3.323,9	2.642,9

The traction battery, regardless of the configuration, reduces the PSPF value. However, the PSPF is below the injury threshold (6 kN) for all cases. The lowest PSPF value is achieved with configuration 2 with shock absorbers.

Figure 33. PSPF value for each of the configurations analyzed.

Figure 33. PSPF value for each of the configurations analyzed.

The modifying parameters of the EuroNCAP assessment protocol, backplate force, T12 vertebra force, and T12 vertebra moment were also obtained.

Table 18. Value of the modifier parameters for each of the configurations analyzed.

Table 18. Value of the modifier parameters for each of the configurations analyzed.

	Without Battery	Config.1	Config.2	Config.1 with Shock-Absorbing Elements	Config.2 with Shock-Absorbing Elements
Fy force in backplate	3.648,2	797,77	1.194,1	3722,3	954,81
Fy force in T12	10.493,0	3.161,2	5.071,2	7.036,6	6.576,8
Mx torque in T12	630,57	305,85	310,55	379,45	335,50

Figure 34. F_y in backplate for each of the configurations analyzed.

Figure 34. F_y in backplate for each of the configurations analyzed.

We observed that the inclusion of the traction battery reduces the value of the force on the backplate, reaching the lowest value in the case of configuration 1.

Figure 35. F_y at the T₁₂ vertebra for each of the configurations analyzed.

Figure 35. F_y at the T₁₂ vertebra for each of the configurations analyzed.

We observed that the inclusion of the traction battery reduces the F_y of the T12 vertebra, reaching the lowest value for configuration 1, and reduces the Mx of the T12 vertebra, reaching the lowest value for configuration 1.

Figure 36. M_x at the T₁₂ vertebra for each of the configurations analyzed.

Figure 36. M_x at the T₁₂ vertebra for each of the configurations analyzed.

4. Conclusions

6. Acknowledgments

References

1. Bakker J., Sachs C., Otte D., Justen R., Hannawald L., Friesen F. Analysis of fuel cell vehicles equipped with compressed hydrogen storage systems from a road accident safety perspective. SAE Int. J. Passeng. Cars Mech. Syst. 4(2011) pp. 332-343.
2. Xuning F., Minggao O., Xiang L., Languang L., Yong X., Xiangming H. Thermal Runaway mechanism of lithium ion battery for electric vehicles: A review. Energy storage materials, 2018, Volume 10, pp. 246-267.
3. Dirección General de Tráfico (DGT). Siniestralidad mortal en vías interurbanas 2022 (cómputo de personas fallecidas a 24 horas): Datos provisionales. 2023. Observatorio Nacional de Seguridad Vial.
4. NHTSA, U.S. Department of Transportation, National Highway Traffic Administration.Laboratory Test Procedure for FMVSS No. 214, Dynamic Side Impact Protection. 2012.
5. BOE. Reglamento n°95 de Naciones Unidas – Prescripciones uniformes sobre la homologación de los vehículos en lo realtivo a la protección de sus ocupantes en caso de colisión lateral (2021/1861). 2021. Available online: https://www.boe.es/buscar/doc.php?id=DOUE-L-2021-81479 (accessed on 3 April 2023).
6. Bahouth G.T., Murakhovskiy D., Digges K.H., Rist H., Wiik R. Opportunities for reducing far-side casualties. 24^th International Technical Conference on the Enhanced Safety of Vehicles (ESV). 2015. Gothenburg, Sweden, National Highway Traffic Safety Administration.
7. Gabler, H.C., et al. Far-side impact injury risk for belted occupants in Australia and the United States. 2005. Proceedings of the 19^th International Conference on Enhanced Safety of Vehicles, 2005, Washington DC, USA.
8. Trumph. El reto: una fabricación rentable y fiable en cuanto a los procesos de baterías de iones litio de alto rendimiento para la electromovilidad. 2023. Available online: https://www.trumpf.com/es_MX/soluciones/sectores/automocion/electromovilidad/celdas-de-bateria-y-modulos/ (accessed on 12 April 2023).
9. Marklines. Nissan Leaf Teardown: Lithium-ion battery pack structure, 2023. Available online: https://www.marklines.com/en/report_all/rep1786_201811 (accessed on 12 April 2023).
10. Maluf, N. Qnovo, 2023. Available online: https://www.qnovo.com/blogs/inside-the-battery-of-a-nissan-leaf (accessed on 12 April 2023).
11. Wang Z., Mao N., Jiang F. Study on the effect of spacing on thermal runaway propagation for lithium-ion battery. Springer 2019. Environmental Science, Engineering Journal of Thermal Analysis and Calorimetry.
12. Setiawan R., Salim M.R. Crashworthiness Design for an Electric City Car against Side Pole Impact. Journal of Engineering and Technological Sciences, 2017, Vol. 49.
13. Air Force and Naval Air Systems Command. Technical manual. Engineering Series for Aircraft Repair.1967. Aerospace Metals - General Data and Usage Factors.
14. Knight Group. 2017. Acero dulce y acero de carbono. Available online: https://www.knight-group.co.uk/wp-content/uploads/PM%20Files/Spanish/Materials/PMQR78%20Spanish%20Knight%20Group%20Carbon%20and%20Mild%20Steels.pdf. (accessed on 12 April 2023).
15. ASM International. Volume 2: Properties and Selection: Nonferrous Alloys and Special-Purpose Materials. 1400 Crystal Drive, Suite 430, Arlington, VA 22209: s.n., 1990.
16. Yue Q., Louis G.H., Christine J. Lithium Concentration Dependent Elastic Properties of Battery Electrode Materials from First Principles Calculations. 2014. Journal of the electrochemical Society, 161 (11) F3010-F3018.
17. Simcenter MADYMO Model Manual. Siemens Version 2212 Model Manual, December 2022.
18. Wang Z., Mao N., Jiang F. Study on the effect of spacing on thermal runaway propagation for lithium-ion battery. Springer 2019. Environmental Science, Engineering Journal of Thermal Analysis and Calorimetry.
19. Brooke L. SAE International. 2021. Available online: https://www.mobilityengineeringtech.com/component/content/article/ae/pub/features/articles/45150?r=40733 (accessed on 9 May 2023).