In this simulation conducted with the LS-DYNA solver, the focus was on a volume representation of an 18650 battery, a widely used cylindrical lithium-ion battery cell, which was then employed to simulate the behavior of the battery under compression mechanical loading. This comprehensive simulation aimed to provide a comparable finite element model with the component test, contributing valuable data for further computing optimization and research into a new method of defining an 18650 battery for large-scale and resource-consuming simulation.
3.1. Model Definition
The computational model intricately replicates the experimental test scenario, featuring a rigid component exerting pressure on the battery, with a rigid wall serving as fixed boundary to emulate the lower support plate of the press. The battery geometry is represented using solid elements within a singular part, characterized by material properties derived from empirical data obtained during the experimental phase. This singular part design aims to encapsulate the nuanced and intricate behavior of the battery in a manner that is both straightforward and comparable. Notably, key parameters, such as the stress–strain curve and compressive force were crucial for accurately capturing the dynamic response of the battery. Given the inherent complexity of the actual system, parameters including density and Young’s modulus underwent optimization through the simulation process to align the model with the intricate characteristics observed in the real-world scenario.
The simulation model comprises a rigid barrier, defined as a shell part consisting of eight quadrilateral elements, each assigned an LS Dyna rigid material property. This rigid barrier is in contact with a battery, which is situated on an LS Dyna rigid wall surface specifically defined as a flat infinite plane (
Figure 8a).
The force reading function of the rigid wall is activated to monitor and record the forces exerted during the simulation. The battery, a solid part, is represented by 2400 hexahedral elements, each with dimensions of approximately 2.5 × 1.7 × 1.7 mm. These elements collectively form a detailed and comprehensive representation of the battery geometry. This is a homogenization model of the hole battery without the caps. The boundary conditions are established by the rigid wall surface, and the force reading function allows for the real-time tracking of the forces applied to the battery during the simulation. This setup, seen in
Figure 8a, enables a thorough examination of the structural behavior and response of the battery under the influence of the rigid barrier, providing important insights into the mechanical aspects of the system.
The same simulation set up is done for the Jelly Roll, but with a diameter offset to match the real dimensions. To add more credibility to the model the case from [
21] is added in a battery model composed of the jelly roll volume and a case (
Figure 8b) tied together with an LS Dyna contact TIED_SURFACE_SURFACE_OFFSET. The thickness of the case is doubled compared to [
21], in order to mimic the usage of the cases in Tesla batterie to 0.64 mm. LS Dyna MAT24 MAT_PIECEWISE_LINEAR_PLASTICITY was used to obtain robust results.
3.2. Material Defining and Parameters’ Calibration
Utilizing LS-DYNA’s MAT24 material model facilitates a comprehensive analysis of structural integrity and deformation behavior in the 18650 lithium-ion battery under compression. MAT24 in LS-DYNA is a versatile material model that combines linear elasticity with viscoelasticity. It is suitable for simulating the time-dependent response of materials, making it valuable for analyzing the dynamic behavior of structures. The material properties in MAT24 are defined in a piecewise manner, meaning that different regions of the stress–strain curve and Prony series can be specified to accurately represent different phases of material behavior. The Prony series allows MAT24 to capture the time-dependent response of the material, essential for simulating dynamic loading conditions. MAT24 uses the Young modulus for the elastic area of deformation but also incorporates a stress–strain curve, allowing for a detailed representation of the material’s plastic deformation characteristics. This curve is defined in a piecewise manner, providing flexibility to capture different phases of material behavior. The material also includes a damping coefficient that accounts for energy dissipation within the material and density.
MAT24 lacks the capability to define anisotropy. However, this limitation has not been a central concern in the present study, as cylindrical lithium-ion batteries primarily exhibit deformation in only two directions. The focus of this paper is specifically on the most critical direction of the battery, during high-impact vehicle crashes, where lateral forces dominate. Axial deformation, while important, occurs less frequently in these high-impact events and is thus outside the scope of this research.
The parameters associated with this material were defined based on the raw testing data and the numerical fitting of the load curve and deformation of the battery. LS SYNA MAT24, a plasticity material model, serves as a foundational element in our simulation framework. This model is adept at simulating the plastic deformation of metals and provides an isotropic, kinematic hardening plasticity model. A notable feature employed in this study is MAT_PIECEWISE_LINEAR_PLASTICITY, which allows for the definition of a piecewise linear stress–strain curve, enhancing the model’s precision in mimicking nonlinear material responses. It is crucial to highlight that, for this homogenized battery model, a stress–strain curve derived from hardware tests was utilized as base guideline. Notably, the piecewise definition was scaled with the method of curve fitting as for the rest of the parameters, the resulting values shown in
Table 2. Additional material parameters include a Young’s modulus of 400 MPa scaled down from the average 515 MPa measured by curve fitting in the elastic area, a Poisson ratio of 0.4, and a density of 2 × 10
−6 (kg/mm
3). For the Jelly Roll, the same method is used but with a difference in E Modulus: the value uses is 200 MPa compared with the average measured 242 MPa. The case material, as also found in [
22], uses an E Modulus of 200 GPa. The plastic domain is defined by a stress–strain curve starting from the yield point of 360 MPa.
It is crucial to acknowledge the inherent challenges associated with the homogenization of a complex system like a Li-ion battery. While the stress–strain curve parameters derived from hardware tests, with modifications in stress and strain scaling, were incorporated into our simulation, it is essential to recognize that these parameters may not possess direct and easily interpretable physical meanings in the context of the homogenized model. The intricacies of a Li-ion battery’s internal structure and composition make it inherently challenging to directly correlate individual material properties from hardware tests with the behavior of the homogenized model. Therefore, while these parameters provide a basis for simulation, their direct interpretation may be limited in the broader context of the complex and multifaceted mechanical interactions within the homogenized Li-ion battery model. The study aims to navigate these challenges and extract valuable insights into the overall compression behavior of Li-ion batteries through careful consideration of the homogenized model and its associated material parameters.
3.3. Discussion
In comparing the experimentally derived load–displacement curve from the uniaxial compression test on the 18650 Li-ion battery with a simulated counterpart, a comprehensive evaluation can be conducted across the three distinctive zones that characterize the mechanical response of the battery. The validation regards the three zones of the force curve described in chapter 2 (
Figure 3 and
Figure 5). In the elastic domain Zone 1, where the battery exhibits linear elasticity, the comparison is focused on the level of force applied and the corresponding displacement. The simulation ideally replicates the reversible deformation observed in this zone, showcasing the battery’s ability to withstand external stress without undergoing irreversible changes. Moving to Zone 2, which marks the transition from elastic to plastic deformation, the comparison is centered around the critical juncture where plastic deformation occurs. The simulation captures the onset of plastic deformation without a significant increase in load, leading to a discernible load plateau. Attention is given to ensuring that the simulated structural compromise aligns with the experimental observations in terms of force level at the displacement of critical points, and the absence of a proportional increase in load. In the final densification zone, Zone 3, the comparison is critical in evaluating the simulated battery’s response to compression-induced compaction. The focus here is on the substantial increase in load as the battery undergoes pronounced densification. The simulation mirrors the experimental findings in terms of force magnitude, displacement characteristics, and the rapid rise in applied load.
Additionally, the rate at which the force increases the stiffness of the battery in each zone is carefully compared between the simulated and experimental curves. Given that we are homogenizing a complex model, a pragmatic approach to analysis was taken. We consider the curve profile to be piecewise, as we are primarily interested in key stiffness points within the compression phenomenon. These points include the beginning of the plastic zone, the onset of densification, and the intensity of densification. Discrepancies in the stiffness application could indicate an error deviation in the simulated material properties or structural behavior. This systematic evaluation based on the criteria of force levels, displacement characteristics, and stiffness across the three zones is pivotal for validating the simulation model. It ensures reliability in informing design, safety, and performance optimization considerations across various applications, providing a robust understanding of the battery’s mechanical behavior.
In the examination of the battery volume deformation, specific focus is directed towards the compression behavior seen in the alignment of the force curves seen in
Figure 9 and the analytical values from
Table 3.
In the actual response curve in the elastic region, designated as Zone 1, the termination of elastic deformation occurs at 0.6 mm press displacement after the initial contact with the battery, registering a force of 2.1 kN and a stiffness of 3.56 kN/mm. Observing
Figure 9, it is discerned that the MAT24 battery volume force curve demonstrates a commendable alignment with the actual curve, but with a more linear behavior. However, comparative analysis, as delineated in
Table 3, also reveals an analytical view of the elastic region with no error noted in the values.
Within the plastic deformation zone, denoted as Zone 2, specific metrics characterize the deformation behavior. The displacement in this region extends to 2.7 mm, resulting in a maximum end force of 2.67 kN and a stiffness of 0.98 kN/mm.
Figure 9 visually indicates a favorable alignment between the force curve and the MAT24 battery volume model. In a quantitative assessment, the MAT24 battery volume model exhibits a peak end force of 2.64 kN with a 1.1% error, and a stiffness of 0.97 kN/mm, reflecting a small 1% discrepancy.
In the densification zone, denoted as Zone 3, a controlled displacement limit of 3.5 mm was imposed to facilitate the calculation of densification rates. Key parameters considered in this context include the initiation point and the subsequent rate of densification. In the actual model, densification commences at 2.7 mm, exhibiting a stiffness of 1.57 kN/mm. The MAT24 battery volume model indicates a stiffness of 1.57 kN/mm, reflecting no analytical deviation. The alignment between the curves during the densification phase is evident in
Figure 9.
In the examination of battery Jelly Roll deformation, as in the battery volume model, specific focus is directed towards the compression behavior observed in the alignment of force curves (
Figure 10) and analytical values (
Table 4).
In the elastic region (Zone 1), termination of elastic deformation occurs at 0.5 mm displacement, with a force of 0.6 kN and a stiffness of 1.2 kN/mm in the actual model. Notably, the MAT24 Jelly Roll volume model demonstrates commendable alignment with the actual curve. However, comparative analysis reveals that the model registers a force of 0.6 kN and a stiffness of 1.2 kN/mm, resulting in a no difference in force and stiffness.
Transitioning to the plastic deformation zone (Zone 2), displacement extends to 2.5 mm, resulting in a maximum end force of 0.99 kN and a stiffness of 0.39 kN/mm.
Figure 10 visually indicates favorable alignment between the force curve and the MAT24 Jelly Roll volume model. However, quantitatively, the model displays values of a peak end force of 1.16 kN, a 17% error, and a stiffness of 0.46 kN/mm, a 17.9% discrepancy. A potential enhancement for future versions of the model would be to introduce more piecewise points within the densification zone. This would allow for a finer resolution of the material’s response during this critical phase, leading to better representation of the compaction process and increased accuracy.
In the densification zone (Zone 3), with a displacement of 3.5 mm, densification commences at 2.5 mm in the actual model, at a peak force of 4, with a stiffness of 1.14 kN/mm. The MAT24 Jelly Roll volume model has a peak force of 3.98, with a 0.5% error, accompanied by a stiffness of 1.13 kN/mm, reflecting 0.8% deviation. This good alignment of the curves in Zone 2 and 3 can also be seen in
Figure 10.
When analyzing the Jelly Roll volume alongside the battery case model, particularly in Zone 1, which denotes the elastic region, distinct areas of elastic deformation become evident. The battery case exhibits higher stiffness initially, followed by a combination of Jelly Roll elasticity and case plasticity, as seen in
Figure 11. This conclusion is not seen in the analytical results presented in
Table 5 because, in cases of high deformation, the point of the force and displacement at the beginning of the plastic area is in focus.
The comparison between various modeling approaches for the battery can be elucidated through the examination of
Figure 11 and
Figure 12.
Figure 12 suggests that the elasticity of the battery volume model aligns more closely with the actual curve in the elastic zone. This can be attributed to the simplified definition of material properties within a homogeneous volume model compared to the complexities involved in modeling materials for heterogeneous models. The deviations observed in the plastic region of the model are also influenced by the interaction with the battery case, though the case tends to reduce these discrepancies to 8.6% and 9.1%, respectively.
3.4. Model Sensibility Test
To ensure the model’s validity and applicability across different configurations, a sensitivity analysis was performed on a battery model composed of a jelly roll volume and its protective case. The analysis involved scaling the plastic strain and case thickness to observe how these modifications affect the densification process and the load curve.
When the plastic strain was scaled down by 10% and 20%, the densification occurred at lower deformations, as expected. Specifically, the initial point of densification shifted from 2.7 mm (original) to 2.5 mm and 2.3 mm, respectively. This indicates that reducing plastic strain sensitivity decreases the deformation required for densification to begin, altering the material response during high compression loads (
Figure 13).
For the case thickness, scaling it down to half, which is representative of a thinner, commercially available Li-ion cell casing, resulted in a nearly halved load curve, reflecting the lower force resistance of the thinner case. Conversely, doubling the case thickness led to an exponential increase in force resistance, with the load curve scaling nearly threefold. This notable effect can be seen in the elastic region, where the increased thickness led to the elasticity of the case having a higher weight in the elasticity of the battery, as shown in
Figure 14. These sensitivity tests highlight how both plastic strain and casing thickness directly influence the battery’s overall mechanical behavior under compression.