Stress Distribution Inside a Lithium-Ion Battery Cell during Fast Charging and Its Effect on Degradation of Separator

The automotive industry is rapidly transitioning to electric vehicles (EVs) in response to the global efforts to reduce greenhouse gas emissions. Lithium-ion battery (LIB) has emerged as the main tool for energy storage in electric vehicles. A widespread adoption of EVs, however, requires a fast-charging technology that can significantly reduce charging time while avoiding any unsafe conditions including short circuits due to failure of the separator in an LIB cell. Therefore, it is necessary to understand the mechanical stresses during fast charging and their long-term effect on the integrity of the separator. This paper presents a novel hybrid model for the prediction of the stress distribution in the separator of a pouch cell under various charging speeds, ambient temperatures, and pack assembly conditions, such as compressive pressures. The proposed hybrid model consists of three sub-models, namely, an electrochemical cell model, a lumped-parameter model, and a solid mechanics model. A robust parameter identification scheme is implemented to determine the model parameters using the experimental data. The separator within the test setup will experience maximum von Mises stress of 74 MPa during 4C charging, i.e., when the charge current in A is four times as high as the capacity of the battery cell in Ah. To assess the evolution of the damage in the separator under the estimated stress during fast charging, creep and fatigue tests are conducted on the separator. Their results indicate a progressive accumulation of damage in the separator, further emphasizing the importance of understanding and mitigating mechanical degradation in separator materials.


Introduction
Due to the volatility of oil markets and strict mandates to limit greenhouse gas emissions, the automotive industry has been swiftly switching towards producing electrical vehicles as well as hybrid electric vehicles (HEVs). As the market share of EVs increases, there will be a greater demand for high efficiency energy storage systems such as lithiumion batteries (LIBs). This is partly due to the high specific energy and energy density of LIBs, which can be two times those of other cell technologies such as nickel-metal hydride [1].
An LIB cell is composed of four main components, i.e., anode, cathode, separator, and electrolyte [2]. In an LIB cell, a critical component is the separator, which is a porous membrane positioned between the electrodes. The essential function of the separator is to prevent direct contact between the electrodes or an internal short while allowing ionic flows through the electrolyte [3]. The separator, therefore being considered the most crucial component for safety, should meet several requirements in regard to its thermal and mechanical properties. First, the separator should be chemically stable under the reductive and oxidative conditions encountered during battery operations [4]. Second, the separator should exhibit low thermal shrinkage, e.g., in an amount of less than 5% after being soaked at 90 • C for 60 min [5]. This requirement exists to avoid any exposure of electrode edges due to shrinkage of the separator at high temperatures, thereby preventing internal shorts and thermal runaways. Moreover, the separator should melt above an elevated temperature to avoid a thermal runaway when the cell is exposed to high temperatures. A shutdown temperature of 130 • C or lower is usually required. Lastly, the separator should be mechanically robust to withstand the mechanical stresses during cell and pack assembly, operations, and any other abuse conditions [6]. Both tensile and puncture strengths of the separator should be high enough to avoid internal shorts under such conditions [7]. A gradual and up to 50% decrease in ductility of the separator was observed after the battery cell underwent 400, 800, and 1600 charge-discharge cycles involving 4C-rate charging, i.e., charging at a rate with current that is four times as high as the battery capacity [8]. A scanning electron micrograph (SEM) analysis of the cycled separators revealed accumulated damage in the separator in the form of pore closure along the thickness direction and fiber fracture and cracks along the machine direction. The degradation was attributed to the damage caused by mechanical creep which can accumulate over the cycles. Given its integral role in battery safety and functionality, an accurate characterization and modeling of the stress distribution in the separator is necessary for addressing the safety concerns under fast-charging conditions.
During battery operation, LIBs are subject to electrochemical processes, thermal energy, and mechanical stresses that may lead to safety concerns [9]. Both the electrochemical process and heat generation induce mechanical loads in the battery components. The former process induces a volumetric change of the active material due to the intercalation processes of lithium-ions, while the latter causes thermal expansion of cell components [10]. Both mechanical loads can squeeze the separator [11], which can lead to shrinkage of its pores, reduced ion transfer, and increase in internal resistance [12]. However, no practical method exists for in situ measurement of the stress inside the battery. Multiphysics modeling based on a finite element method is often adopted instead to predict the stress distribution. Xiao et al. [13] proposed a multiphysics and multiscale model for prediction of the stress in a polymeric separator by coupling a 1-dimensional (1-D) macroscopic battery model and a meso-scale battery sub-model. Shi et al. [14] built a multiphysics cell model for a pouch cell to predict the stress in the separator. A 1-D thermo-electrochemical cell model calculates lithium concentration and generates heat while a macroscale 2-D model of a pouch cell determines its stress distributions. However, such cell models by themselves lack the capability to predict the stress inside a cell according to its external conditions imposed by the battery pack. A battery cell will be subject to varying compressive pressure and ambient thermal conditions according to both design and use parameters of the battery pack. Swelling of a battery cell during charging and discharging and varying compressive pressure when the cell is constrained inside a battery pack are often modeled by a lumpedparameter modeling approach. Oh et al. [15] developed a thermal swelling model to predict the swelling in the battery due to temperature variation. In their other works, Oh et al. [16] used a phenomenological force and swelling model to predict the dynamic swelling in the battery. Xu et al. [17] developed a five Kelvin-Voigt elements model to characterize the viscoelastic behavior of the separator. This model enables the prediction of creep for two specific commercial separators. In their additional research, Xu et al. [18] formulated an anisotropic homogeneous model for predicting the mechanical behavior of an LIB under bending, indentation, and compression loads. Moreover, the numerical model took into account the state of charge dependence on the mechanical behavior of the lithium-ion battery (LIB). Liu et al. [19] created a multiphysics model that captures the mechanical behavior of cylindrical cells. The model also addresses short-circuit and thermal runaway behaviors. Duan et al. [20] developed a 2-D coupled multiphysics model to identify various internal short circuit modes. The model is capable of describing the entire evolution process of the battery starting from the initial deformation, triggering of internal short circuit, evolution of internal short circuit, up to the eventual thermal runaway. Yuan et al. [21] developed a generalized short-circuit criterion for two common failure modes, using volumetric strain and equivalent strain of the separator as key parameters. In this paper, a computationally efficient hybrid model is developed by linking the cell model in [14] with a lumped-parameter model to account for the pack parameters in addition, thereby allowing optimal design of both cell and pack parameters for safety. The proposed hybrid model consists of three main components: (1) a 1-D electrochemical model, (2) a lumped-parameter model (LPM), and (3) a 2-D solid mechanics model. The 1-D electrochemical model predicts the heat generated in the battery cell and lithium-ion concentration which are input to the LPM. The LPM represents the mechanical and thermal behaviors of a battery cell in a battery pack array. The predicted lithium-ion concentration and temperature are input to the solid mechanics model for prediction of the stress in the separator. The LPM parameters are determined via non-linear regression based on the temperature, mechanical load, and displacement data collected from a fast charging and discharging test on an LIB cell. The hybrid model, thus developed, allows calculating the stress in the separator during fast charging. New insights into damage accumulation in the separator under fast cyclic charging are obtained by performing creep and fatigue tests on the separator under the predicted stress.
This paper is organized as follows: Section 2 describes the hybrid model proposed in this study. Then, Section 3 provides detailed explanation on the material properties of the tested LIB cell and the testing setup. Section 4 presents the experimental results along with the parameter identification process. In Section 5, a simulation study is conducted using the developed model to check the effect of charging rate, pack assembly condition, and ambient temperature on the stress induced in the separator. Section 6 focuses on presenting the results of creep and fatigue tests conducted under the predicted stress during fast charging. Finally, the conclusion is presented in Section 7.

Hybrid Model
During a charge and discharge cycle, an LIB cell is subject to electrochemical, thermal, and mechanical phenomena occurring simultaneously. Each phenomenon is described by a set of equations and coupled with others. The framework of the proposed hybrid model is divided into three components. The first component is a 1-D model for predicting transport of lithium-ions, electrochemical reactions, and generated heat given charge and discharge rates. The 1-D model is composed of three domains of anode, separator, and cathode. The second model component is a lumped-parameter model (LPM) that will predict the load, displacement, and temperature in the battery cell and its fixture. This model is composed of three parts: a thermal model for heat transfer, a spring model for the mechanical forces, and a linking equation which links the thermal model to the spring model. The LPM takes both predicted heat and state of charge (SOC) from the 1-D electrochemical model as inputs. Additional input variables of the LPM include ambient temperature and pre-torque on the fixture. The third component is a 2-D solid mechanics model which will calculate the stress distribution in the separator given (1) cell core temperature and boundary pressure predictions from the LPM and (2) lithium concentration in the electrodes as predicted by the 1-D battery model. Figure 1 shows all three components with the input/output variables. In what follows, a detailed description of each model component is provided.

One-Dimensional (1-D) Battery Model
The 1-D battery model, shown in Figure 2, accounts for the electrochemical reaction and heat generation in a battery cell. An LIB cell Li-ion is composed of three components placed in the electrolyte, namely, cathode, anode, and separator.

One-Dimensional (1-D) Battery Model
The 1-D battery model, shown in Figure 2, accounts for the electrochemical reaction and heat generation in a battery cell. An LIB cell Li-ion is composed of three components placed in the electrolyte, namely, cathode, anode, and separator.

One-Dimensional (1-D) Battery Model
The 1-D battery model, shown in Figure 2, accounts for the electrochemical and heat generation in a battery cell. An LIB cell Li-ion is composed of three com placed in the electrolyte, namely, cathode, anode, and separator.
is the active material volume fraction, is the electrolyte volume frac rosity), and denotes the binder and filler volume fraction. The superscr notes the electrode under consideration ( for positive electrode and for nega trode). In what follows, the governing equations for the diffusion, kinetics, an conservation [22,23] are summarized for the solid and electrolyte phases.

Solid Phase
The following equations apply to the solid phase of the anode and cathode. T material in the electrodes is assumed to be composed of spherical particles. Th second law describes the transport in the solid electrode particles and predicts th of concentration field with time in spherical coordinates: = where is the particle radius coordinate, is the lithium concentration in the s ticle phase, and is the diffusion coefficient of lithium-ion in solid electrode p The equation assumes a constant diffusion coefficient of lithium-ions in the soli According to the symmetry condition, the first boundary condition yields − since there is no flux at the center of the particle. The second boundary conditio the pore wall flux across the interface and the rate of diffusion of lithium-ions solid/liquid interface, yielding − = , where is the particle radius a the pore wall flux of lithium-ions defined through the Butler-Volmer equation (E (6)). The conservation of charge in the solid phase is governed by Ohm's law: The volume fraction of the electrodes is defined by: where ε i s is the active material volume fraction, ε i l is the electrolyte volume fraction (porosity), and ε i f iller denotes the binder and filler volume fraction. The superscript i denotes the electrode under consideration (p for positive electrode and n for negative electrode). In what follows, the governing equations for the diffusion, kinetics, and charge conservation [22,23] are summarized for the solid and electrolyte phases.

Solid Phase
The following equations apply to the solid phase of the anode and cathode. The active material in the electrodes is assumed to be composed of spherical particles. The Fick's second law describes the transport in the solid electrode particles and predicts the change of concentration field with time in spherical coordinates: where r is the particle radius coordinate, c i s is the lithium concentration in the solid particle phase, and = J i , where r p is the particle radius and J i is the pore wall flux of lithium-ions defined through the Butler-Volmer equation (Equation (6)).
The conservation of charge in the solid phase is governed by Ohm's law: where τ i e f f is the effective conductivity determined by τ i e f f = ε i s 1.5 τ i s , φ i s is the potential in the solid phase, a v is the particle specific surface area, and F r denotes the Faraday's constant. The first boundary condition applies to the interface between the cur-

Liquid Phase
The following equations apply to the electrolyte in the electrodes and separator. To account for lithium species transport in the liquid, the concentrated solution theory is adopted [22,24,25]. The material balance for a salt in liquid phase is given by: where where κ e f f ,i is the effective conductivity, ϕ l is the potential in the liquid phase, R is the ideal gas constant, and T is the temperature. At the two ends of the cell, there is no current flow in the liquid phase, i.e.,

Charge Transfer
The pore wall flux at the interface between the electrode and electrolyte is determined by the Butler-Volmer equation which describes the relationship between the electrical current, electrode potential, and local concentration at the interface. Specifically, the rate of transfer of lithium-ions from the electrolyte to an electrode is controlled by deviation of the surface potential from its equilibrium value as follows: where i ex is the exchange current density, η i is the overpotential of the intercalation reaction, and α a and α c are the charge transfer coefficients for oxidation and reduction, respectively, which are assumed to be 0.5. Both the forward rate of the anodic process and the backward rate of the cathodic process are given in Equation (6). The exchange current density is given by: where k i is the reaction rate constant, c i l is the concentration of lithium-ions in the electrolyte, c i s,sur f is the concentration of lithium-ions on the surface of the solid electrode, and c i s,max − c i s,sur f is the available space for lithium intercalation. The overpotential is given by: where φ i s is the potential of the solid electrode, φ i l is the potential of the electrolyte, and U i is the open circuit potential.

Heat Generation
Two main mechanisms of heat generation are considered in this study, namely, reversible entropy heat and irreversible joule heat. The corresponding equation for heat generation rate, Q gen , is as follows: where (φ s − φ l − U i )J represents the joule heat while T ∂U i ∂T J represents the entropy heat.

Lumped-Parameter Model
In a lumped-parameter model, a battery cell is represented based on equivalent material properties such as thermal conductivity and specific heat [15,26,27]. In addition, a uniform heat generation is assumed throughout the battery cell [28][29][30].

Thermal Model
A schematic diagram of the pouch cell sandwiched between two fixture plates is shown in Figure 3. The lumped-parameter model assumes a 1-D thermal behavior of the pouch cell with a conduction boundary condition via contacts with the plates. Moreover, a convective boundary condition exists on both the top and bottom of the plates.
where is the potential of the solid electrode, is the potential of the electrolyte, an is the open circuit potential.

Heat Generation
Two main mechanisms of heat generation are considered in this study, namely, re versible entropy heat and irreversible joule heat. The corresponding equation for heat gen eration rate, , is as follows: where − − ) represents the joule heat while represents the entrop heat.

Lumped-Parameter Model
In a lumped-parameter model, a battery cell is represented based on equivalent ma terial properties such as thermal conductivity and specific heat [15,26,27]. In addition, uniform heat generation is assumed throughout the battery cell [28][29][30].

Thermal Model
A schematic diagram of the pouch cell sandwiched between two fixture plates i shown in Figure 3. The lumped-parameter model assumes a 1-D thermal behavior of th pouch cell with a conduction boundary condition via contacts with the plates. Moreove a convective boundary condition exists on both the top and bottom of the plates. The transient one-dimensional conductive heat transfer equation is as follows: where is the cell heat capacity, ) is the heat generated in the battery, i the core temperature of the battery, and , ) is the surface temperature of the batter while Κ, , and denote the thermal conductivity, area, and length, respectively. Th convective heat transfer at the top boundary of the plate can be expressed as: Moreover, a two-state representation of the 1-D thermal model is given by: The transient one-dimensional conductive heat transfer equation is as follows: where C cell is the cell heat capacity, Q gen (t) is the heat generated in the battery, T C is the core temperature of the battery, and T S,1 (t) is the surface temperature of the battery while K, A, and L denote the thermal conductivity, area, and length, respectively. The convective heat transfer at the top boundary of the plate can be expressed as: Moreover, a two-state representation of the 1-D thermal model is given by: where T S,2 (t) is the surface temperature of the plate, T ∞ is the ambient temperature, and h air is the air convection coefficient.

Spring Model
A spring model was developed to represent the mechanical behavior of a battery cell in the fixture. Figure 4 illustrates the behavior of the cell at constrained conditions. Two equivalent stiffness elements are adopted, i.e., the equivalent stiffness s Total , representing the cell and the plate, and the equivalent stiffness s Rod , representing the rod. These equivalent stiffness values refer to the relationship between the force applied on the fixture and the displacement of the fixture along the vertical axis. A linear model is adopted to predict the dependence of the force on the displacement (Branch A): where s Total is the linear equivalent stiffnesses, and δ c,% Total represents the constrained displacement of the battery and plates. The superscript c represents the constrained displacement conditions. There are four rods connected in parallel at the four edges of the plate. The total force in Branch B can be described by: where ∆ bolt is the displacement in the rod due to the tightening of the bolt and δ f ,% b is the free unconstrained displacement of the battery. The total measured displacement is given by: where δ c 0 is the initial displacement.

Spring Model
A spring model was developed to represent the mechanical behavior in the fixture. Figure 4 illustrates the behavior of the cell at constrained c equivalent stiffness elements are adopted, i.e., the equivalent stiffness the cell and the plate, and the equivalent stiffness , representing the ro alent stiffness values refer to the relationship between the force applied on the displacement of the fixture along the vertical axis. A linear model is ad the dependence of the force on the displacement (Branch A): is the linear equivalent stiffnesses, and ,% represents displacement of the battery and plates. The superscript represents the c placement conditions. There are four rods connected in parallel at the fo plate. The total force in Branch B can be described by: where ∆ is the displacement in the rod due to the tightening of the b the free unconstrained displacement of the battery. The total measured given by: is the initial displacement. At equilibrium, the plate and the battery cell will be compressed whi will be stretched. Therefore, the plates and the battery cell will be under co and the rods will be subject to tensile force. At mechanical equilibrium, t will be equal: At equilibrium, the plate and the battery cell will be compressed while the four rods will be stretched. Therefore, the plates and the battery cell will be under compressive force and the rods will be subject to tensile force. At mechanical equilibrium, these two forces will be equal: Batteries 2023, 9, 502 8 of 24

Coupling between Thermal and Mechanical Models
In the proposed LPM, the total swelling of the battery cell, ∆L Total , consists of (1) swelling due to thermal expansion, ∆L Thermal , and (2) swelling due to Li-ion intercalation process, ∆L Intercalation . Assuming uniform heat generation across the cell, the temperature distribution in the battery will follow a parabolic shape (See Figure 3). The thermal swelling normal to the battery cell is given by: (18) where ∆L T is the thermal swelling (due to the cell heating/cooling), L Cell is the thickness of the cell, α is the lumped coefficient of thermal expansion, and T is the temperature. The subscript S denotes the surface initial state. Integrating Equation (18) based on the parabolic temperature distribution across the cell gives: where T C (t), T S (t), and T Re f denote cell-core temperature, cell-surface temperature, and reference temperature, respectively. Swelling due to the Li-ion intercalation during charging and discharging is represented by a linear equation: where ξ is the coefficient of expansion according to the state-of-charge (SOC) of the cell. An average SOC over thickness was used in this model. The thermal model in Section 2.2.1 and the spring model in Section 2.2.2 are linked via the free swelling term as follows:

Solid Mechanics Model
A 2-D solid mechanics model is adopted to predict the stress in a battery cell. The model was developed in COMSOL Multiphysics 5.6. Schematic representation of the mechanics model for a pouch cell is shown in Figure 5. The part modeled is highlighted using a dotted box. The negative current collector, anode electrode, separator, cathode electrode, and positive current collector are shown in different colors. Two types of boundary conditions exist in the model, namely, internal and external boundary conditions. During battery production, the separator is tightly wrapped around the electrodes after electrodes and current collectors are bonded together. Therefore, a set of bonded electrode and current collector is considered as a continuous medium. In contrast, a contact-type boundary condition, allowing a slip at the interface, was assigned between cathode and separator and between anode and separator, respectively. External boundary conditions include a symmetric boundary condition on the left and bottom side of the model and a free boundary on the right side. Three main inputs were used in the solid mechanic model; a distributed normal load is applied at the top according to the predicted force on the battery and fixture plates by the LPM, i.e., F Branch A (see Figure 4). Additionally, the core temperature from the LPM is assumed to be uniform throughout the solid mechanics model since the cell is very thin [14], and an average lithium concentration was used for the electrodes. The constitutive equations of the model include those based on mechanical strains , thermal strain [14], and intercalation-induced strain [13,31] to predict the stress in the battery components as follows: where is the Young's modulus, is the Poisson's ratio, is the stress tensor, and is the Kronecker delta, is the thermal expansion coefficient, ∆ is the change in temperature from a reference value, ∆ is the change in the diffusion species concentration from a reference (stress-free) value, and Ω is the partial molar volume [31,32]. The partial molar volume can be calculated using [31]:

Remark 2. Equation (22) only considers elastic deformation and does not consider any plastic
deformation. This is a safe assumption since the amount of deformation in this study is very small and does not exceed the elastic limit [8,14,[33][34][35]. It should be also noted that the strain induced by intercalation, denoted by , is exclusively applicable to the positive and negative electrodes. Additionally, Equation (22) considers a constant Young's modulus in the mechanical strains for the electrodes. The Young's modulus was considered independent of the lithium concentration. In [36], the authors showed that elastic stiffening increases the elastic modulus, which increases the stress. However, to maintain simplicity in this work, no concentration-dependent material properties were considered.  Without an extra cathode on top of the positive current collector, the cathode below will form an asymmetric beam which tends to deflect upward away from the separator due to intercalation-induced stresses.

Battery Cell
The constitutive equations of the model include those based on mechanical strains ε me ij , thermal strain ε T ij [14], and intercalation-induced strain ε D ij [13,31] to predict the stress in the battery components as follows: where E v is the Young's modulus, ν v is the Poisson's ratio, σ ij is the stress tensor, and δ ij is the Kronecker delta, α v ij is the thermal expansion coefficient, ∆T is the change in temperature from a reference value, ∆c is the change in the diffusion species concentration from a reference (stress-free) value, and Ω i is the partial molar volume [31,32]. The partial molar volume can be calculated using [31]: (22) only considers elastic deformation and does not consider any plastic deformation. This is a safe assumption since the amount of deformation in this study is very small and does not exceed the elastic limit [8,14,[33][34][35]. It should be also noted that the strain induced by intercalation, denoted by ε D ij is exclusively applicable to the positive and negative electrodes. Additionally, Equation (22) considers a constant Young's modulus in the mechanical strains for the electrodes. The Young's modulus was considered independent of the lithium concentration. In [36], the authors showed that elastic stiffening increases the elastic modulus, which increases the stress. However, to maintain simplicity in this work, no concentration-dependent material properties were considered. More details about the PE separator are as follows. The PE separator (12 µm thick) is fabricated through a co-extrusion blowing process. The separator porosity is 43% with a crystallinity ranging between 60% and 80%. The tensile strength in the machine direction (MD) is 235 MPa while that in the transversal direction (TD) is 225 MPa. The PE separator properties are summarized in Table 1. In a battery cell, the pores in the electrodes and separator are filled with electrolyte. Therefore, the effective thermal expansion coefficient and Young's modulus of the electrodes and separator were calculated considering their porosities. The effective thermal expansions of the electrodes and separator were estimated by the rule of mixture as follows:

Battery Cell
where α i and V i are the thermal expansion coefficient and the volume fraction of the ith phase, respectively. The effective Young's modulus can be calculated in a similar manner. The material properties for the 1-D battery cell model and the 2-D solid mechanics model are listed in Table 2. The use of '*' next to a value indicates that it is an effective material property.

Cell Fixture and Instrumentation
A fixture, shown in Figure 6, was used to emulate a battery module assembly. The fixture was built by holding two steel plates (each with a size of 125 mm × 99 mm × 12.7 mm) with four threaded rods (9.5 mm in diameter). A thermocouple was attached on the top of the battery cell to measure the cell surface temperature while another thermocouple was installed on the top of the upper plate. An additional thermocouple was used to monitor the ambient temperature. Four compression load cells (Omega LC9150-500-1K) were placed between the upper plate and hex nuts. A high precision displacement sensor (Keyence LK-G10) was placed 10 mm above the upper plate to measure the expansion of the battery cell. The whole setup was placed inside a thermal chamber whose temperature was set to 20 • C. After the battery cell was charged at a constant current-constant voltage (CC-CV) to 3.6 volts at 0.5 C, with a 5% current cutoff (0.0875 A), it was allowed to rest for 1 h for thermal equilibrium. Then, a pre-torque in the amount of 1 Nm was applied on each hex nut to compress the battery cell. top of the battery cell to measure the cell surface temperature while another thermocouple was installed on the top of the upper plate. An additional thermocouple was used to monitor the ambient temperature. Four compression load cells (Omega LC9150-500-1K) were placed between the upper plate and hex nuts. A high precision displacement sensor (Keyence LK-G10) was placed 10 mm above the upper plate to measure the expansion of the battery cell. The whole setup was placed inside a thermal chamber whose temperature was set to 20 °C. After the battery cell was charged at a constant current -constant voltage (CC-CV) to 3.6 volts at 0.5 C, with a 5% current cutoff (0.0875 A), it was allowed to rest for 1 h for thermal equilibrium. Then, a pre-torque in the amount of 1 Nm was applied on each hex nut to compress the battery cell.

Identification of LPM Parameters
This section describes how the parameters of the lumped-parameter model (LPM) are identified.

Experiments for Measurement of Stiffness
Both force and displacement were measured under several torques to find the equivalent spring constant of battery cell and plates. The fixture shown in Figure 6 was used for the experiment. The battery was initially charged to 3.6 volts and placed between the plates. No thermocouples were used in this experiment. Four compression load cells (Omega LC9150-500-1K) were placed between the upper plate and hex nuts. The torque on each hex nut was increased by an increment of 0.2 N·m from 0 to 1 N·m and by an increment of 0.1 N·m from 1 to 1.2 N·m. Force and displacement measurements were recorded at each increment. The equivalent stiffness was fitted on the force and the displacement experimental data are shown in Figure 7a. The spring constant for the battery was found to be = 9.0892 × 10 N m .

Identification of LPM Parameters
This section describes how the parameters of the lumped-parameter model (LPM) are identified.

Experiments for Measurement of Stiffness
Both force and displacement were measured under several torques to find the equivalent spring constant of battery cell and plates. The fixture shown in Figure 6 was used for the experiment. The battery was initially charged to 3.6 volts and placed between the plates. No thermocouples were used in this experiment. Four compression load cells (Omega LC9150-500-1K) were placed between the upper plate and hex nuts. The torque on each hex nut was increased by an increment of 0.2 N·m from 0 to 1 N·m and by an increment of 0.1 N·m from 1 to 1.2 N·m. Force and displacement measurements were recorded at each increment. The equivalent stiffness was fitted on the force and the displacement experimental data are shown in Figure 7a. The spring constant for the battery was found to be s Total = 9.0892 × 10 6 N·m −1 .
Another experiment was carried out to find the stiffness coefficient of the threaded rod. A rod was placed vertically with a lock-nut at the bottom. An Omega load cell was placed between the lock-nut and a hex nut. As the applied torque on the hex nut increased from 0 to 1.6 N·m, the force and elongation of the rod were measured. The rod spring constant was fitted on the force and the elongation data are shown in Figure 7b. The turn angle was used to calculate the elongation in the rod using the pitch. Note this spring constant corresponds to the portion of the rode between the two nuts. The Young's modulus can then be found using: Using the Young's Modulus, we can then calculate the spring constant of the rod for its length in the fixture. The spring constant for the rod was found to be S Rod = 1.0311 × 10 7 N·m −1 .
(Omega LC9150-500-1K) were placed between the upper plate and hex nuts. The torque on each hex nut was increased by an increment of 0.2 N·m from 0 to 1 N·m and by an increment of 0.1 N·m from 1 to 1.2 N·m. Force and displacement measurements were recorded at each increment. The equivalent stiffness was fitted on the force and the displacement experimental data are shown in Figure 7a. The spring constant for the battery was found to be = 9.0892 × 10 N · m .

Charging and Discharging Test
The force, displacement, and temperatures during charging and discharging were measured after the battery cell was installed in the fixture described in Section 3.2. The collected data were used for tuning the LPM parameters as described in Section 4.3. The test was performed under ambient temperature of 20 • C.
The battery was initially charged to 3.6 volts with a constant current of 1.75 A (0.5C), and then the voltage was held at 3.6 volts until the current dropped to 5% of its initial value. Figure 8a shows the voltage and current at the start of discharge from 3.6 volts. The data consist of three parts based on constant current-constant voltage (CC-CV) protocols. The first part of the curve shows CC-CV discharging at 0.5C from 3.6 to 3 volts (0-2178 s). The second part is a period of CC-CV fast charging at 4C from 3 to 4.2 volts (2179-3407 s), and the third part corresponds to CC-CV discharging at 0.5C from 4.2 to 3 volts (3408-9470 s). Figure 8 also shows variations of load and displacement overlayed on the voltage curve. Starting with the displacement figure (see Figure 8b), the discharge between 3.6 volts and 3 volts decreases the displacement (height) of the top surface of the fixture due to shrinkage of the battery cell. During the period of fast charging at 4C, a sharp increase in the displacement due to rapid expansion of the battery cell can be observed. During the voltage-hold period after fast charging, a small increase in the displacement is followed by a decrease. This is due to an initial increase in the temperature during charging followed by a decrease, causing cell shrinkage, during the hold period. The maximum displacement reached during fast charging was 45.8 µm. During discharging, the displacement decreases sharply between 4.2 volts and 3.7 volts, followed by a small decrease between 3.7 and 3.6 volts and another sharp decrease between 3.6 and 3 volts. The load follows the same trend of the displacement curve. After tightening the bolts with a 1 N·m torque, the initial load applied on the battery was 1314 N at the beginning of the test. The load decreases to 983 N during the first discharging, then increases significantly during charging to reach 2333 N at maximum followed by a small decrease to 2213 N during the hold period. During the second discharging period, the load decreases to 1070 N.

LPM Parameters
Some parameters of the LPM were obtained from the literature while others identified by fitting the model to the experimental data. The thermal conductivity o cell for the LPM, for example, can be determined by combining the material prop reported in the literature. The thermal conductivity of a battery cell, composed of m layers of different conductivities, will vary depending on the direction of heat flow [3 The through-plane thermal conductivity, in the direction of the heat flow consider the LPM, can be expressed as [37]: where Κ is the thermal conductivity of the material layer and is the thickness o material layer. The equivalent thermal conductivity of the cell, thus determined, is f to be 0.633 . The thermal conductivity and thermal expansion coefficient of the are 167 , and 25.2 × 10 −6 , respectively, from the literature. The experimental data used for tunning of the remaining LPM parameters in those collected during CC discharging from 3.6 to 3 volts, CC fast charging from 3 volts, and CC discharging from 4.2 to 3 volts. The responses and predictors of the are summarized in Table 3. Four LPM parameters were estimated based on the ex mental data: (1) the heat capacity of the cell, , (2) the heat transfer coefficient o ℎ , (3) the thermal expansion coefficient of the cell, , and (4) the SOC coefficie expansion due to lithiation, . The parameters were estimated such that the su squared errors in predicting the temperatures, , and , , the force, , and the placement, ,% , is minimized as follows:

LPM Parameters
Some parameters of the LPM were obtained from the literature while others were identified by fitting the model to the experimental data. The thermal conductivity of the cell for the LPM, for example, can be determined by combining the material properties reported in the literature. The thermal conductivity of a battery cell, composed of many layers of different conductivities, will vary depending on the direction of heat flow [30,37]. The through-plane thermal conductivity, in the direction of the heat flow considered in the LPM, can be expressed as [37]: where K v is the thermal conductivity of the material layer and t v is the thickness of the material layer. The equivalent thermal conductivity of the cell, thus determined, is found to be 0.633 W mK . The thermal conductivity and thermal expansion coefficient of the plate are 167 W mK , and 25.2 × 10 −6 1 K , respectively, from the literature. The experimental data used for tunning of the remaining LPM parameters include those collected during CC discharging from 3.6 to 3 volts, CC fast charging from 3 to 4.2 volts, and CC discharging from 4.2 to 3 volts. The responses and predictors of the LPM are summarized in Table 3. Four LPM parameters were estimated based on the experimental data: (1) the heat capacity of the cell, C Cell , (2) the heat transfer coefficient of air, h air , (3) the thermal expansion coefficient of the cell, α Cell , and (4) the SOC coefficient for expansion due to lithiation, ξ. The parameters were estimated such that the sum of squared errors in predicting the temperatures, T s,1 and T s,2 , the force, F, and the displacement, δ EXP,% , is minimized as follows:

Ambient temperature T Amb
A non-linear regression was performed using the least-squares curve fitting in MAT-LAB. The parameter values, thus estimated, are listed in Table 4. The estimated parameter values are in the same order of magnitude as those in the literature, e.g., in [15]. The results of model-fitting are shown in Figure 10. The experimental battery surface temperature, the plate surface temperature, and the fixture displacement are in good agreement with the LPM predictions. However, a relatively large error in prediction of force can be observed at the end of charging and during discharging. In order to evaluate the goodness of model-fitting, the root mean square error (RMSE) is calculated (see Table 5): where y actual i is the experimental value, y predicted i is the value predicted in the model, and N is the number of data points. Additionally, the RMSE divided by the range and average are also listed. It can be seen that, for each output, its RMSE is below 20% of its range indicating a good model-fit.  In order to evaluate the goodness of model-fitting, the root mean square error (RMSE) is calculated (see Table 5): where is the experimental value, is the value predicted in the model, and N is the number of data points. Additionally, the RMSE divided by the range and average are also listed. It can be seen that, for each output, its RMSE is below 20% of its range indicating a good model-fit.

Simulation Study
As mentioned in the introduction, a significant degradation in the separator's integrity can be observed after the battery cell undergoes cyclic 4C charging [8]. Therefore, the Figure 10. Comparison of the simulated output and the experimental data after model fitting.

Simulation Study
As mentioned in the introduction, a significant degradation in the separator's integrity can be observed after the battery cell undergoes cyclic 4C charging [8]. Therefore, the stress in the separator under 4C charging is estimated using the developed hybrid model. Moreover, the effect of charging rate, pre-torque, and ambient temperature on the stress in the separator is investigated.

Estimated Stress in the Separator
The stress in the separator under 4C charging condition of the cycling test in [8] was calculated using the hybrid model. Specifically, the core temperature, boundary pressure, and lithium concentration predicted using the electrochemical model and the LPM model were used as inputs for the solid mechanics model (See Figure 1). Figure 11a shows a homogeneous distribution of the von Mises stress in the middle section of the separator at the end of 4C charging, while Figure 11b presents the evolution of the maximum von Mises stress over time during charging. The von Mises stress is calculated as follows [38]: where σ 1 , σ 2 , and σ 3 are the principal stress components. The von Mises stress serves as a measure of the overall stress experienced by the material at a specific location, accounting for both normal and shear stresses.
where , , and are the principal stress components. The von Mises stress serves as a measure of the overall stress experienced by the material at a specific location, accounting for both normal and shear stresses.
The von Mises stress showed a homogeneous distribution within the central part of the separator during the 4C charging. The maximum von Mises stress is found to be 74.1 MPa at the end of fast charging. Creep and fatigue behavior of the separator under the estimated stress will be presented in Section 6.

Effect of Charging Rate
The hybrid model allows examination of the effect of charging rate on the separator. Evolution of the von Mises stress in the separator over time was analyzed at four different charging rates of 1C, 2C, 4C, and 5C with other conditions, such as the pre-torque and ambient temperature, fixed at those in Section 5.1. The CC-CV charging protocol was employed for each charging rate. Figure 12 shows the temperature variations during charging for different charging rates. When the charging rate is 1C, the cell temperature slightly increases from 20 °C to the maximum value of 20.25 °C at the end of the constant current (CC) period. As the charging rate increases, the maximum cell temperature increases, e.g., 24.6 °C at 5C. It can be seen that a higher charging rate results in an elevated cell temperature, which can potentially compromise battery safety and lead to accelerated degradation of cell components including the separator. Therefore, a careful consideration of the charging rate is essential for optimal battery performance and longevity. The von Mises stress showed a homogeneous distribution within the central part of the separator during the 4C charging. The maximum von Mises stress is found to be 74.1 MPa at the end of fast charging. Creep and fatigue behavior of the separator under the estimated stress will be presented in Section 6.

Effect of Charging Rate
The hybrid model allows examination of the effect of charging rate on the separator. Evolution of the von Mises stress in the separator over time was analyzed at four different charging rates of 1C, 2C, 4C, and 5C with other conditions, such as the pre-torque and ambient temperature, fixed at those in Section 5.1. The CC-CV charging protocol was employed for each charging rate. Figure 12 shows the temperature variations during charging for different charging rates. When the charging rate is 1C, the cell temperature slightly increases from 20 • C to the maximum value of 20.25 • C at the end of the constant current (CC) period. As the charging rate increases, the maximum cell temperature increases, e.g., 24.6 • C at 5C. It can be seen that a higher charging rate results in an elevated cell temperature, which can potentially compromise battery safety and lead to accelerated degradation of cell components including the separator. Therefore, a careful consideration of the charging rate is essential for optimal battery performance and longevity. Evolution of the von Mises stress in the separator over time at the four different charging rates is shown in Figure 13. Regardless of the charging rate, the von Mises stress Evolution of the von Mises stress in the separator over time at the four different charging rates is shown in Figure 13. Regardless of the charging rate, the von Mises stress exhibits a linear increase during the constant charge (CC) period followed by a slower increase during the constant voltage (CV) period. A faster increase in the von Mises stress is observed with an increase in the charging rate. It should be noted that the proposed hybrid model is capable of predicting the evolution of stress for any charging protocols. The hybrid model, therefore, opens the possibility of optimizing the charging schedule for achieving multiple objectives such as prolonging both electrode health and separator integrity. Evolution of the von Mises stress in the separator over time at the four different charging rates is shown in Figure 13. Regardless of the charging rate, the von Mises stress exhibits a linear increase during the constant charge (CC) period followed by a slower increase during the constant voltage (CV) period. A faster increase in the von Mises stress is observed with an increase in the charging rate. It should be noted that the proposed hybrid model is capable of predicting the evolution of stress for any charging protocols. The hybrid model, therefore, opens the possibility of optimizing the charging schedule for achieving multiple objectives such as prolonging both electrode health and separator integrity.

Effect of Pre-Torque
The effect of pre-torque on separator stress is examined by applying various pretorque values to the hybrid model. Different pre-torques ranging from 0.5 to 12 N·m were simulated. Note, in the experiment in Section 3, a pre-torque of 1 N·m was applied to the fixture. Figure 14 shows only a minor change in the von Mises stress with increasing pre-

Effect of Pre-Torque
The effect of pre-torque on separator stress is examined by applying various pretorque values to the hybrid model. Different pre-torques ranging from 0.5 to 12 N·m were simulated. Note, in the experiment in Section 3, a pre-torque of 1 N·m was applied to the fixture. Figure 14 shows only a minor change in the von Mises stress with increasing pre-torque. The reason why the pre-torque plays a minor role in the separator stress is as follows. The stress in the separator can be attributed to two main mechanisms. First, during charging, the anode expands in all directions due to lithiation and a temperature increase. Expansion of anode will exert a tensile load on the separator as the anode pulls the separator in the horizontal direction, as seen in Figure 5. Second, the pressure due to the pre-torque and structural constraint imposed by the fixture will compress the separator in the vertical direction. The stress due to pulling is much larger than the one due to the pressure by at least two orders in this particular configuration. Therefore, the change in pressure due to varying pre-torque will have only a minor effect on the von Mises (overall) stress.

Effect of Ambient Temperature
The effect of ambient temperature on separator stress was investigated by applying six different ambient temperatures ranging from 0 to 50 • C to the hybrid model. Figure 15 shows the relationship between the von Mises stress in the separator at the end of charging and the ambient temperature. The results show a decrease in the von Mises stress with increasing ambient temperature.
Expansion of anode will exert a tensile load on the separator as the anode pulls the separator in the horizontal direction, as seen in Figure 5. Second, the pressure due to the pretorque and structural constraint imposed by the fixture will compress the separator in the vertical direction. The stress due to pulling is much larger than the one due to the pressure by at least two orders in this particular configuration. Therefore, the change in pressure due to varying pre-torque will have only a minor effect on the von Mises (overall) stress.

Effect of Ambient Temperature
The effect of ambient temperature on separator stress was investigated by applying six different ambient temperatures ranging from 0 to 50 °C to the hybrid model. Figure 15 shows the relationship between the von Mises stress in the separator at the end of charging and the ambient temperature. The results show a decrease in the von Mises stress with increasing ambient temperature.

Effect of Ambient Temperature
The effect of ambient temperature on separator stress was investigated by applying six different ambient temperatures ranging from 0 to 50 °C to the hybrid model. Figure 15 shows the relationship between the von Mises stress in the separator at the end of charging and the ambient temperature. The results show a decrease in the von Mises stress with increasing ambient temperature. The cell temperature will increase when the ambient temperature rises. This will lead to expansion of battery components in both horizontal and vertical directions (see Figure 5). Since the thermal expansion coefficient of the separator is significantly larger than those of the other battery components, the separator will experience a pronounced length expansion in the horizontal direction, which leads to less pulling by the anode and lower stresses.

Creep and Fatigue Testing
As previously discussed, noticeable degradation of the separators' mechanical properties (by as much as 50%) was observed when the LIB, made of the same separator as the one in this study, was subject to repeated fast (4C) charging [8]. In Section 5, the maximum von Mises stress during 4C charging was found to be 74 MPa while only a small temperature rise was observed. Therefore, creep and fatigue tests were carried out on a pristine separator at room temperature to examine its behavior under the estimated stress during 4C charging.
Strip-shaped uniaxial tensile specimens were prepared according to ASTM D882. The samples were cut from a new separator roll using a razor. The longer axis the sample was aligned in the machine direction was with a length of 100 mm and a width of 10 mm. The gage length was chosen to be 50 mm. The creep and fatigue experiments were conducted using an ElectroPuls E1000. Foamy jaws were used to hold the specimen from each end.
Creep testing was conducted at five target stresses ranging from 50 to 90 MPa, which are lower than the yield stress of the separator (ranging from 100 to 140 MPa in the machine direction depending on the strain rate, according to the authors' previous study [8]). Note the maximum von Mises stress during 4C charging was estimated at 74.1 MPa in Section 5.1. In each creep testing, the loading increased linearly from 0 MPa to the target stress at a rate of 0.1 MPa/s. Once a target stress is reached, the loading remained at the same stress level until the displacement reaches 60 mm (except for the target stress of 50 MPa since the maximum duration of the creep test was 6 days). The creep test results are presented in Figure 16.  The fatigue test was designed to mimic the stress evolution in the separator during cyclic fast charging for 68 cycles. Each cycle consists of (1) an increasing load along a linear ramp (0 to 70 MPa in 13 min) representing a 4C charging duration, (2) a stress hold representing a CV period for 10 min, and (3) a linear decreasing ramp (70 to 0 MPa in 90 min) representing a discharging duration. Figure 17 shows the variations of load and displacement over time during the first cycle. The displacement curve shows a rapid increase during the charging period, followed by a gradual increase in the displacement during the holing period and the beginning part of discharging. This is followed by decrease in the displacement for the remaining period of discharging. Even after the load returns to 0, a non-zero elongation of 1.3 mm remains, indicating a permanent damage in the separator. The fatigue test was designed to mimic the stress evolution in the separator during cyclic fast charging for 68 cycles. Each cycle consists of (1) an increasing load along a linear ramp (0 to 70 MPa in 13 min) representing a 4C charging duration, (2) a stress hold representing a CV period for 10 min, and (3) a linear decreasing ramp (70 to 0 MPa in 90 min) representing a discharging duration. Figure 17 shows the variations of load and displacement over time during the first cycle. The displacement curve shows a rapid increase during the charging period, followed by a gradual increase in the displacement during the holing period and the beginning part of discharging. This is followed by decrease in the displacement for the remaining period of discharging. Even after the load returns to 0, a non-zero elongation of 1.3 mm remains, indicating a permanent damage in the separator.
representing a discharging duration. Figure 17 shows the variations of load and displacement over time during the first cycle. The displacement curve shows a rapid increase during the charging period, followed by a gradual increase in the displacement during the holing period and the beginning part of discharging. This is followed by decrease in the displacement for the remaining period of discharging. Even after the load returns to 0, a non-zero elongation of 1.3 mm remains, indicating a permanent damage in the separator. The results of fatigue testing for 68 cycles are shown in Figure 18. The displacement over time is shown in Figure 18a while the maximum displacement of each cycle is plotted in Figure 18b. In Figure 18b, the displacement increases rapidly during the first 10 cycles and then shows a linear increase for the remaining cycles. It should be noted that the fatigue test demonstrates an extreme case of separator degradation during cycling. Once the The results of fatigue testing for 68 cycles are shown in Figure 18. The displacement over time is shown in Figure 18a while the maximum displacement of each cycle is plotted in Figure 18b. In Figure 18b, the displacement increases rapidly during the first 10 cycles and then shows a linear increase for the remaining cycles. It should be noted that the fatigue test demonstrates an extreme case of separator degradation during cycling. Once the separator elongates permanently after the first few cycles, it will be subject to stresses much lower than 70 MPa in the proceeding cycles. Nevertheless, the test results indicate that the load on the separator during 4C charging is high enough to cause damage to the separator. Moreover, the damage can be accumulated in the separator via repeated loading, explaining the separator degradation after cyclic fast charging in the previous study [8]. separator elongates permanently after the first few cycles, it will be subject to stresses much lower than 70 MPa in the proceeding cycles. Nevertheless, the test results indicate that the load on the separator during 4C charging is high enough to cause damage to the separator. Moreover, the damage can be accumulated in the separator via repeated loading, explaining the separator degradation after cyclic fast charging in the previous study [8].

Conclusions
Combining three sub-models enabled predicting the stress in the separator under various use conditions including those for battery pack assembly and thermal management. The parameters of a lumped-parameter model were determined based on the experimental data showing good model-fit, while the material properties in other sub-mod-

Conclusions
Combining three sub-models enabled predicting the stress in the separator under various use conditions including those for battery pack assembly and thermal management. The parameters of a lumped-parameter model were determined based on the experimental data showing good model-fit, while the material properties in other sub-models were obtained from the literature. The developed model allowed estimating the maximum von Mises stress as 74 MPa during the fast-charging condition, under which a significant separator degradation was observed after 1600 cycles in the previous study. In order to investigate the damage to the separator under this loading condition, creep and fatigue tests were conducted on a pristine separator. The high creep sensitivity and creep deformation in the test results indicate that such loading can produce and accumulate significant damage in the separator. However, care should be taken to avoid concluding that the mechanical stress is the only driving force behind separator degradation. The high temperature, which was not considered due to the low temperature rise of the setup in this study, may play a role by causing thermal degradation, which is a major mechanism for degradation of polymers in general. While the current state of the model allows for an efficient way in predicting the stresses in the cell across various charge rates, torque levels, and ambient temperature conditions, further investigations should be conducted to explore a more complex solid mechanics model which includes plasticity and damage in the constitutive equations. Additionally, the model can be further improved by incorporating concentrationdependent material properties to improve the mechanical predictions.  Data Availability Statement: Data will be available on request.

Acknowledgments:
The authors would like to thank Tom Brackett for conducting the Li-ion batteries charging and discharging test.

Conflicts of Interest:
The authors declare no conflict of interest.