Design and Structural Simulations of a Custom Li-Po Accumulator for Low Range, Lightweight, Single-Seater, Open Cockpit, and Open-Wheeled Racecar

: Electric, hybrid, and fuel cell vehicles are the future of the automobile industry, and power source design is one of the most crucial steps in designing these vehicles. This paper aims to design and structurally simulate a custom accumulator—which powers an electric vehicle, for a lightweight, single-seater formula-style racecar. The work is dependent on the model-based design and CAD model approach. Mathematical modeling on SCILAB is used to model equations to get the characteristics of the accumulator, such as the energy, capacity, current, voltage, state of charge, and discharge rates. The output of this model gives the conﬁguration of the battery pack as several cells in series and parallel to adequately power the tractive system. An accumulator container is designed to safeguard the cells from external impacts and vibrational loads, which otherwise can lead to safety hazards. Following this, the Finite Element Analysis (FEA) performed on the accumulator resulted in maximum peak deformation of 0.56 mm, ensuring the safety check against various external loads. Further, the ﬁner stability of the battery pack was virtually validated after performing the vibrational analysis, resulting in a deformation of 3.5493 mm at a 1760.8 Hz frequency.


Introduction
With climate change becoming a severe issue, the world is embracing renewable energy sources and accepting the different sustainable transportation modes.Growing concerns over carbon emission have lately sparked worldwide interest in more efficient and cleaner transportation systems.To mitigate pollution, the world is shifting to electric vehicles as better alternatives to fuel-based automobiles.To support this movement, the Government of India has pledged to make India an all-electric vehicle nation by 2030 [1].This project is aligned with the global agenda of manufacturing zero-emission vehicles hence enforcing stricter Environmental Protection Acts worldwide.This paper discusses the design and structural simulations of a custom accumulator for use in a lightweight, single-seater formula-style racecar, which will encourage motorsports teams and companies to shift from the traditionally used internal combustion engines to an EV powertrain for use in racing events [2].The high manufacturing cost, increased weight, increased volume, and need for proper thermal management of Li-ion battery cells are currently the main

Methodology
The accumulator or the battery pack, which powers the electric vehicle, is a custombuilt battery pack that includes all the cells, busbars, maintenance plugs, Battery Management System (BMS), Accumulator Isolation Relay (AIRs), fuses, and cooling system to manage the battery and the power distribution for the motor controller.This paper discusses the design of a safe, high-performance accumulator that could provide the maximum power of 52 kW to the tractive system and store sufficient energy to allow the car to travel up to 22 km in one charge.As the battery pack is intended to be used in motorsport applications, performance and safety are paramount.The research demonstrated uses several tools.Drive cycle data with speed at closely spaced time steps are acquired from a lap time simulation software (Optimum Lap) for the vehicle model as its primary input.Mathematical modeling on SCILAB is then used to model the various equations governing the resistive forces, speed, torque, and power delivered by the propulsion system, which in our case is the electric motor, and then ultimately the characteristics of the accumulator such as the energy, capacity, current, voltage, state of charge and discharge rates.The output of this model gives the configuration of the battery pack giving us the number of cells in series and parallel to power the tractive system adequately.As the accumulator would be mounted on the vehicle, it needs to withstand the static and vibrational loads encountered in regular use.Keeping this in mind, the accumulator is designed on CAD software (Fusion 360) and simulated structurally (on Simscale) to handle these loads using finite element analysis.Since the battery packs operate at voltages over 200 V, which can be dangerous, the containers also protect the individuals who work around them from electrical shocks.The materials and components selected for the design considered safety the main feature to be satisfied, and then a lightweight and reliable container was designed that can efficiently run the vehicle.Figure 1 shows the flow chart of the methodology followed.The accumulator is also designed to be easily accessible and serviceable to facilitate easy maintainability.Further, the finer stability of the battery pack was virtually validated after performing the vibrational analysis using ANSYS.The car under consideration is a singleseater, lightweight, formula-style car that can compete in amateur racing competitions such as Formula Society of Automotive Engineers (FSAE) events held worldwide [12].

Cell Selection
In an electric vehicle, cell selection dictates the design of the accumulator [8].Hence, choosing the right type of cell for the design to start the design is a crucial step.This decision, taken very early in the design process, affects the tractive system design and the efficiency of the entire vehicle.The most common type of cells used in automobile applications is lithium-based cells, which have lithium metal oxides as positive electrodes and lithium salts as electrolytes for the transfer of electrons.As the accumulator is designed for lightweight vehicles, energy density by weight and volume become important factors.Cells with specific energy densities are preferred as they would result in a smaller and lighter battery pack.As the vehicle is an open cockpit, thermal management of the cells becomes difficult, especially in tropical climates.Cells that can withstand higher temperatures are referable as some cell chemistries tend to decrease charge holding capacities at higher temperatures.This cell selection method is not easy to produce, as cells come in several different chemistries and different types of packaging, manufactured by different manufacturers, each with its advantages and disadvantages.Various cell chemistries are available today; each cell chemistry was studied and analyzed for our application.Table 1 shows the comparison of some of the shortlisted cell chemistries.Based on the above comparison, lithium phosphate cell chemistry was selected as it had the highest energy density and was economical.It can also provide a high operating current (1 C-5 C) and has good performance even at temperatures near the top of its operating limits.They also have good overcharge capabilities and long life cycles, ideal for use in small-range vehicles.They also have a high tolerance to abuse due to their relatively moderate energy density.Internationally, there are many manufacturers of cells, including A123 systems, Kokam, EIG, K2 Energy, Thundersky, and Melasta [9].Based on carefully examining the data sheets received from each vendor, pricing, availability, and packaging interpretations, A123 systems were selected.A123 systems are a worldwide leader in lithium-ion energy storage.They provide cells and cell modules for use in automotive applications [10].Cells are mainly available in three types of enclosures or geometries-cylindrical cells, prismatic cells, and pouch cells.Cylindrical cells are enclosed in a cylindrical outer shell, which allows them to maintain their shape against the

Cell Selection
In an electric vehicle, cell selection dictates the design of the accumulator [8].Hence, choosing the right type of cell for the design to start the design is a crucial step.This decision, taken very early in the design process, affects the tractive system design and the efficiency of the entire vehicle.The most common type of cells used in automobile applications is lithium-based cells, which have lithium metal oxides as positive electrodes and lithium salts as electrolytes for the transfer of electrons.As the accumulator is designed for lightweight vehicles, energy density by weight and volume become important factors.Cells with specific energy densities are preferred as they would result in a smaller and lighter battery pack.As the vehicle is an open cockpit, thermal management of the cells becomes difficult, especially in tropical climates.Cells that can withstand higher temperatures are referable as some cell chemistries tend to decrease charge holding capacities at higher temperatures.This cell selection method is not easy to produce, as cells come in several different chemistries and different types of packaging, manufactured by different manufacturers, each with its advantages and disadvantages.Various cell chemistries are available today; each cell chemistry was studied and analyzed for our application.Table 1 shows the comparison of some of the shortlisted cell chemistries.Based on the above comparison, lithium phosphate cell chemistry was selected as it had the highest energy density and was economical.It can also provide a high operating current (1-5 C) and has good performance even at temperatures near the top of its operating limits.They also have good overcharge capabilities and long life cycles, ideal for use in small-range vehicles.They also have a high tolerance to abuse due to their relatively moderate energy density.Internationally, there are many manufacturers of cells, including A123 systems, Kokam, EIG, K2 Energy, Thundersky, and Melasta [9].Based on carefully examining the data sheets received from each vendor, pricing, availability, and packaging interpretations, A123 systems were selected.A123 systems are a worldwide leader in lithium-ion energy storage.They provide cells and cell modules for use in automotive applications [10].Cells are mainly available in three types of enclosures or geometries-cylindrical cells, prismatic cells, and pouch cells.Cylindrical cells are enclosed in a cylindrical outer shell, which allows them to maintain their shape against the expansion during charging.Prismatic cells also have a metal external enclosure.This makes the cells heavier and reduces the energy density per unit weight of the cell.Considering the various types of cell geometry, we chose the A123 pouch cell with 3.2 V and 20 Ah capacity.Using a pouch cell reduces the weight of the accumulator considerably, hence aiding the vehicle performance.The electrodes in pouch cells lay parallel to each other in a single plane, forming a flat rectangular shape with the terminals at the same edge of the cell.Lithium-ion pouch cells need to be under pressure to operate at peak performance [10].However, the electrical, mechanical, and thermal integration of cells into packs and packs into electric vehicles is paramount to ensure long and safe operation [11,12].

Lap Time Simulations
Optimum Lap is a lap-time simulation software that uses a point mass model of the car and simulates its run around a user-defined track.The car's specifications, like weight, torque, rpm, suspension characteristics, etc., can vary as per the user input [13].The track considered for the simulation is used in many motorsport races and Formula Bharat events at Kari Motor Speedway, Coimbatore, India.For the car's inputs, a typical light formula-style racecar used in FSAE events was taken as a reference with a weight of 300 Kg, powered by an Emrax 188 motor with a peak power of 52 kW and peak torque of 90 nm [14].Figure 2 shows the variation of speed on the simulated track.To calculate the motor torque and speed, the final reduction ratio of the transmission setup had to be fixed.The track was carefully mapped, and simulations were run for varying final drive ratios (FDR) for the car to achieve the fastest possible times.An FDR of 5 was chosen for further discussion as it provides good results in both full lap and acceleration runs and would also be easy to package mechanically.
Considering the various types of cell geometry, we chose the A123 pouch cell w 3.2 V and 20 Ah capacity.Using a pouch cell reduces the weight of the accumulator co siderably, hence aiding the vehicle performance.The electrodes in pouch cells lay par lel to each other in a single plane, forming a flat rectangular shape with the terminals the same edge of the cell.Lithium-ion pouch cells need to be under pressure to opera at peak performance [10].However, the electrical, mechanical, and thermal integrati of cells into packs and packs into electric vehicles is paramount to ensure long and sa operation [11,12].

Lap Time Simulations
Optimum Lap is a lap-time simulation software that uses a point mass model of t car and simulates its run around a user-defined track.The car's specifications, li weight, torque, rpm, suspension characteristics, etc., can vary as per the user input [1 The track considered for the simulation is used in many motorsport races and Formu Bharat events at Kari Motor Speedway, Coimbatore, India.For the car's inputs, a typi light formula-style racecar used in FSAE events was taken as a reference with a weig of 300 Kgs, powered by an Emrax 188 motor with a peak power of 52 kW and pe torque of 90 nm [14].Figure 2 shows the variation of speed on the simulated track.calculate the motor torque and speed, the final reduction ratio of the transmission set had to be fixed.The track was carefully mapped, and simulations were run for varyi final drive ratios (FDR) for the car to achieve the fastest possible times.An FDR of 5 w chosen for further discussion as it provides good results in both full lap and accelerati runs and would also be easy to package mechanically.

SCILAB Modeling and Simulations-Powertrain Sizing
Models can be classified according to the different modeling approaches used [1 The major categories are mathematical models, electrochemical models, and electri equivalent circuit networks [16].For our application, the mathematical modeling a proach is used using SCILAB.SCILAB uses a numerically oriented programming la guage to simulate mathematical models.The powertrain sizing model, modeled usi SCILAB, gives the total number of cells, the number of cells in series and parallel, a battery pack volume based on various input parameters like drive-cycle data, wheel

SCILAB Modeling and Simulations-Powertrain Sizing
Models can be classified according to the different modeling approaches used [15].The major categories are mathematical models, electrochemical models, and electrical equivalent circuit networks [16].For our application, the mathematical modeling approach is used using SCILAB.SCILAB uses a numerically oriented programming language to simulate mathematical models.The powertrain sizing model, modeled using SCILAB, gives the total number of cells, the number of cells in series and parallel, and battery pack volume based on various input parameters like drive-cycle data, wheel radius, and all the forces acting under total tractive effort load.Apart from the drive cycle data given to the model, these parameters are shown below in Table 2.

Transmission Parameters
Gear ratio 5 Transmission efficiency 0.85 The powertrain-sizing model consists of five main system blocks-Chassis, Transmission, Motor, Motor Controller, and Battery in Figure 3.The input parameters like drive-cycle data, wheel radius, and all the forces acting under total tractive effort given to the Chassis system block result in wheel speed and wheel torque.The formulae used for this calculation are mentioned below in Equations ( 1) and (2).
where, The output of the transmission system, i.e., motor speed and motor torque, is used to calculate motor power in the motor system block.The battery power obtained as an output from the motor controller system block is used to calculate the total number of cells, cells in series and parallel in the last battery system block.The range of output values of each system block can be seen in the above figure.Figure 4 shows the graphs obtained as an output from these system blocks.The output of the transmission system, i.e., motor speed and motor torque, is used to calculate motor power in the motor system block.The battery power obtained as an output from the motor controller system block is used to calculate the total number of cells, cells in series and parallel in the last battery system block.The range of output values of each system block can be seen in the above figure.Figure 4 shows the graphs obtained as an output from these system blocks.The output of the transmission system, i.e., motor speed and motor torque, is used to calculate motor power in the motor system block.The battery power obtained as an output from the motor controller system block is used to calculate the total number of cells, cells in series and parallel in the last battery system block.The range of output values of each system block can be seen in the above figure.Figure 4 shows the graphs obtained as an output from these system blocks.At any instant, the maximum resistive force will be the summation of peak values from each of the above graphs.Here, it is evident that Aerodynamic peak force = 200 N, Acceleration peak force = 1000 N, Gradient peak force = 0 N, Rolling peak force = 44 N. Hence the maximum instantaneous resistive load that our motor should be able to handle is calculated as the sum of the four forces (peak) which is equal to 1245 N (approximately).The rolling force here is constant because rolling resistance depends on the product of weight of the object and coefficient of rolling resistance.In our case, we consider the weight of the car to be constant throughout and coefficient is a constant.
Figure 5 represents the speed and torque for wheels and motor, respectively, and Figure 6 depicts the required motor power.These serve as the minimum requirements that must be sought while sizing our components.We can surmise from the above graphs and relations that any deviation in the values obtained from resistive forces directly impacts the motor power calculation, which will eventually affect the battery sizing, as we will see from the relations further on.The battery system block shown in Figure 7 is divided further into various sub-systems.Each sub-system contains its calculation related to the battery.The various calculations are completed using the following Equations ( 3)- (9).
where, BP = Battery power in kW MC e f f iciency = Motor controller e f f iciency Energy per km = ( BP)/(D × 3.6 × 10 6 ) where, D = Drive cycle range in km Total energy = Energy per km × D (5) where, BC rate = Battery Discharge rate BC = Battery nominal capacity in Ah where, SOC = Battery state o f charge where, where, The battery power, voltage, current, c-rate, and SOC are calculated under this system to ensure that the cell configuration achieved from them successfully powers the vehicle to complete the distance in a single charge.Figure 7 shows the battery system block modeled in Scilab software.
Figure 8 shows the battery power and battery current graphs obtained from these sub-system blocks, which helps determine the battery pack capacity using the relations mentioned above.
Another critical application of a battery model is to estimate battery states, such as state of charge and state of health, that are key battery operating states yet not directly measurable, using model-based estimation methods [17,18].
The simulated battery discharge C rates and state of charge curve are depicted in Figure 9.A discharge rate curve helps us understand the usage of our cells.It helps us determine how much current over the continuous and peak rating of the individual cell are we drawing from the accumulator.The results show that the continuous discharge rate is 3 C while peak discharge is at 5.3 C.However, the datasheet of 20 Ah A123 cells suggests that the usable continuous discharge current can go up to 200 A, which is 10 C, and the usable peak discharge current can go up to 600 A, which is 30 C. This can be seen in the image below (Table 3).The battery power, voltage, current, c-rate, and SOC are calculated under this system to ensure that the cell configuration achieved from them successfully powers the vehicle to complete the distance in a single charge.Figure 7 shows the battery system block modeled in Scilab software.Figure 8 shows the battery power and battery current graphs obtained from these sub-system blocks, which helps determine the battery pack capacity using the relations mentioned above.The battery power, voltage, current, c-rate, and SOC are calculated under this system to ensure that the cell configuration achieved from them successfully powers the vehicle to complete the distance in a single charge.Figure 7 shows the battery system block modeled in Scilab software.Figure 8 shows the battery power and battery current graphs obtained from these sub-system blocks, which helps determine the battery pack capacity using the relations mentioned above.The battery power, voltage, current, c-rate, and SOC are calculated under this system to ensure that the cell configuration achieved from them successfully powers the vehicle to complete the distance in a single charge.Figure 7 shows the battery system block modeled in Scilab software.Figure 8 shows the battery power and battery current graphs obtained from these sub-system blocks, which helps determine the battery pack capacity using the relations mentioned above.Another critical application of a battery model is to estimate battery states, such as state of charge and state of health, that are key battery operating states yet not directly measurable, using model-based estimation methods [17,18].
The simulated battery discharge C rates and state of charge curve are depicted in Figure 9.A discharge rate curve helps us understand the usage of our cells.It helps us determine how much current over the continuous and peak rating of the individual cell are we drawing from the accumulator.The results show that the continuous discharge rate is 3 C while peak discharge is at 5.3 C.However, the datasheet of 20 Ah A123 cells suggests that the usable continuous discharge current can go up to 200 A, which is 10 C, and the usable peak discharge current can go up to 600 A, which is 30 C. This can be seen in the image below (Table 3).The battery pack sizing model gives us the required total energy, which, when divided by the required motor operating voltage, gives the required pack capacity in Ah.This is then divided by the chosen cell capacity to obtain the number of cells in parallel.The battery pack sizing model gives us the required total energy, which, when divided by the required motor operating voltage, gives the required pack capacity in Ah.This is then divided by the chosen cell capacity to obtain the number of cells in parallel.Several cells in series are found by dividing the pack voltage known by a single shortlisted cell voltage.
Hence, our usage is well within safe operating limits and bolsters our belief that despite the undulating high and low discharge rates, the battery pack charge does not fall below 20% at the end of the 2500 s endurance instance.The motor power obtained from the output graph of Scilab's transmission system model is in the range between nominal power and the peak power of the chosen motor.Hence, we can conclude that the battery pack would provide sufficient power to the motor during its operation, even during the 22 km of the car's endurance.
The number of cells obtained in series is 90, and the number of cells obtained in parallel is 1.Hence, a 90S1P was selected as the best cell configuration to get the best performance out of the accumulator for dynamic events.

CAD Modeling
The cells and cell configuration obtained from the previous steps must be arranged and packaged in an accumulator box for safe use.As there are no guides for the mechanical design of an accumulator, the rulebook from Formula Student Germany was used as a guide for the design [19].This document provides a set of rules that must be followed to ensure the safe use of the accumulator under racing conditions.
The rules state that one accumulator segment cannot exceed a maximum voltage of 120 V DC, the maximum energy of 6 MJ, and a mass of 12 kg.To facilitate this, the accumulator was divided into six smaller "modules," containing 15 in 3 columns of 5 cells each, as shown in Figure 10.The cells were connected in series.The modules were designed to enclose the cells to provide structural stability.The individual modules were connected electrically in series to achieve the required 90S1P configuration.The six modules were placed inside an external accumulator box and the maintenance plugs and other necessary components.The accumulator box must protect the cells under lateral, longitudinal, and vertical loading.The accumulator must also be electrically insulated.As the design uses pouch cells, constant pressure is applied to the face of the cells for maximum efficiency.The design must also provide easy assembly, accessibility, and ease of maintenance.
Several iterations of the accumulator box and modules were made, each an improvement over the other.The goal was to maximize the structural strength of the accumulator without increasing the weight too much.
Neoprene rubber foam, a heat retardant, and insulating material were used to maintain pressure on each cell face while providing the necessary electrical insulation.To maintain the series connection, busbars and clamps were designed according to the cell orientation.Positive locking through double tab washers was used to safeguard against the motor vibrations acting on the module.Considering the active cooling system, the module container was designed with slots for fans to vent the air in and out of the system using.Handles were designed for easy access to the module containers.A 2.5 mm thick aluminum 7075-T6 material was chosen for the box, considering its strength to weight ratio [20].
A maintenance plug for 3/0 AWG gauge wire was used with a latching connector and a plug terminal on the other end to connect the modules in series.The lug terminal sides of the maintenance plug are connected to cell tabs, and the two latching connectors are used to connect and disconnect two adjacent modules, making it safe to disconnect for maintenance.An accumulator container was designed around these modules.For ease of assembly and disassembly, one side opening is provided apart from the top opening, as seen in Figure 11.The thickness of the base plate is 3.5 mm, while the thickness of the vertical walls is 2.5 mm.The mass of the whole accumulator system was 75.3 kg considering all the PCBs for monitoring cell temperature and voltage, maintenance plugs, fasteners, and brackets for mounting the accumulator container, and it was acceptable.Figure 12 clearly shows this assembled accumulator configuration.A maintenance plug for 3/0 AWG gauge wire was used with a latching connector and a plug terminal on the other end to connect the modules in series.The lug terminal sides of the maintenance plug are connected to cell tabs, and the two latching connectors are used to connect and disconnect two adjacent modules, making it safe to disconnect for maintenance.An accumulator container was designed around these modules.For ease of assembly and disassembly, one side opening is provided apart from the top opening, as seen in Figure 11.The thickness of the base plate is 3.5 mm, while the thickness of the vertical walls is 2.5 mm.The mass of the whole accumulator system was 75.3 kg considering all the PCBs for monitoring cell temperature and voltage, maintenance plugs, fasteners, and brackets for mounting the accumulator container, and it was acceptable.Figure 12 clearly shows this assembled accumulator configuration.

Finite Element Analysis
FEA and finite element modeling (FEM) are numeric techniques that use computational power to approximate the solution to stress and strain estimation equations.They are used when the shapes or forces are too complicated to be solved by analytical methods [21,22].Simscale, an online cloud computing Computer-Aided Engineering (CAE)

Finite Element Analysis
FEA and finite element modeling (FEM) are numeric techniques that use computational power to approximate the solution to stress and strain estimation equations.They are used when the shapes or forces are too complicated to be solved by analytical methods [21,22].Simscale, an online cloud computing Computer-Aided Engineering (CAE) software, is used for structural simulations to provide access to higher computational power.Following the rulebook, the container, the mounting of the accumulator container to the chassis, and the mounting of the cells to the container must be designed to withstand the following accelerations with mounting points as fixed support: For analysis of ductile materials with complex geometries, Von Mises equivalent stress is used.Von Mises stress, σv, is given below in Equation ( 10): 10) For the design to be safe, the maximum value of the σv should be less than the material's yield point.

Vibrational Analysis
One typical approach to characterize the observed dynamic behavior of a mechanical system is performing modal analysis [23].Modal analysis or structural vibrational analysis is the study of the dynamic properties of a body under natural vibrational fre-

Vibrational Analysis
One typical approach to characterize the observed dynamic behavior of a mechanical system is performing modal analysis [23].Modal analysis or structural vibrational analysis is the study of the dynamic properties of a body under natural vibrational frequency conditions.This would simulate the displacement or deformation of the part when subjected to vibration and loads during its lifetime [24].This technique can explain a structure in terms of its natural characteristics.The natural characteristics are natural frequency, damping, and mode shapes [25].The modal model of the system is formulated mathematically, and the modal data is the characteristic information of the system.Modal analysis is mostly dispensed to seek out the system's natural frequency and for different mode shapes at different eigenvalues [26].The modal analysis relies on expressing the response of a linear time-invariant dynamic system as a linear simple harmonic motion set of natural vibrations.The damping factor is defined by the solver depending upon the material properties and the part's dimensions.Examples of strategies for modal analysis range from phase-retrieval combined with direct mode projection, correlating the modes with spectral or temporal degrees of freedom, combining principal component analysis after adapting the detection system with a training data set, to performing a coordinate transformation that converts the beam into a more convenient basis [27].The fundamental of modal analysis using measured frequency response function data is curve fitting the data using a predefined mathematical model of the measured structure.This model assumes the number of Degrees of Freedom (DoFs) of the structure, its damping type, and possibly the number of vibration modes within the measured frequency range.These assumptions dictate the mathematical expression of each Frequency Response Function (FRF) curve from measurement and a curve-fitting process trying to derive all modal parameters in a mathematical formula of an FRF using measurement data [28].The analysis is carried out with a flow structure of Static Structural → Modal → Random Vibrations → Lifespan estimation.
Modal analysis is carried out for frequencies ranging from 20 to 2000 Hz.This includes all the possible frequencies generated by the motor and vibrations.

Static Structural Analysis
The analysis is carried out on ANSYS Mechanical.The CAD model is imported into the software, and certain prerequisites are determined and set.The given CAD model meshes in ANSYS Meshing, where the body is divided into several elements suitable for mathematical analysis.The material for the accumulator is defined to be Aluminum Alloy 7075-T6 grade.Remote points were defined on the module walls on which the pouch cells are mounted.The mass of the modules is defined as point masses with remote attachments to said remote points.The weight of each pouch cell was taken to be 496 g [29], and the weight of one module is defined to be 7.44 kg.The accumulator's mounting holes are defined as the fixed points in static structural analysis to calculate the total deformation due to the weight of the pouch cells.The total deformation in the static condition is shown in Figure 15.

Modal Analysis
The results of the static structural analysis are linked to the modal frequency analysis to correlate the results.For frequencies between 20 and 2000 Hz [30], the number of modes must be analyzed.Each mode denotes a particular natural frequency, modal damping factor, and displacement pattern, and mode shape [23].Due to these natural modes, the vibrations are determined by the properties of the vibration source and mode shapes of the system.The mode v/s frequency graph is shown in Figure 16.cells are mounted.The mass of the modules is defined as point masses with remote attachments to said remote points.The weight of each pouch cell was taken to be 496 g [29], and the weight of one module is defined to be 7.44 kg.The accumulator's mounting holes are defined as the fixed points in static structural analysis to calculate the total deformation due to the weight of the pouch cells.The total deformation in the static condition is shown in Figure 15.

Modal Analysis
The results of the static structural analysis are linked to the modal frequency analysis to correlate the results.For frequencies between 20 and 2000 Hz [30], the number of modes must be analyzed.Each mode denotes a particular natural frequency, modal damping factor, and displacement pattern, and mode shape [23].Due to these natural modes, the vibrations are determined by the properties of the vibration source and mode shapes of the system.The mode v/s frequency graph is shown in Figure 16.The total deformation is calculated for modes 1-650, ranging from 17.895 to 1767.3 Hz.The minimum to maximum total deformation experienced by the accumulator ranges from 0.4890 to 3.5493 mm for 17.85-1760.8Hz, respectively.The deformation results can be seen in Figure 17.The total deformation is calculated for modes 1-650, ranging from 17.895 to 1767.3 Hz.The minimum to maximum total deformation experienced by the accumulator ranges from 0.4890 to 3.5493 mm for 17.85-1760.8Hz, respectively.The deformation results can be seen in Figure 17.

Random Vibrations
The mode shapes are extracted from the total deformation results as the solution to modal analysis.The solution of modal analysis is used as an input to the setup of random vibrations.The Segalman-Reese method is used for the dynamic analysis.The displacement response PSD (Power Spectral Density) can be simplified for a single base excitation [30].A log-log plot between Actual Spectral Density (ASD) and frequency is used for performing random vibration analysis.The extracted modes are used as precur-

Random Vibrations
The mode shapes are extracted from the total deformation results as the solution to modal analysis.The solution of modal analysis is used as an input to the setup of random vibrations.The Segalman-Reese method is used for the dynamic analysis.The displacement response PSD (Power Spectral Density) can be simplified for a single base excitation [30].A log-log plot between Actual Spectral Density (ASD) and frequency is used for performing random vibration analysis.The extracted modes are used as precursors to the setup of random vibrations.The power spectral density acceleration graph is defined in the setup.The ASD (G2/Hz) V/s frequency (Hz) graph is shown in Figure 18, and the input values are shown in the table.The software will hence interpolate the values for G2/Hz for the frequencies for each mode.The input data determines the directional deformation values for scale factors sigma 1, sigma 2, and sigma 3. Maximum directional deformation was observed in the X direction, with negligible deformation in the Y and Z directions.The equivalent stresses under Von Mises stress are calculated for the scale factors.Table 4 shows scale factor values, respective probability, and maximum directional deformation.The input data determines the directional deformation values for scale factors sigma 1, sigma 2, and sigma 3. Maximum directional deformation was observed in the X direction, with negligible deformation in the Y and Z directions.The equivalent stresses under Von Mises stress are calculated for the scale factors.Table 4 shows scale factor values, respective probability, and maximum directional deformation.As the probability of the sigma 3 scale factor is the highest, it is hence considered for total directional deformation and equivalent Von Mises stress.The results of directional deformation of scale factors sigma 3 are shown in Figure 19.
where, The analysis resulted in a total lifespan of 1.047 × 1026 h, which is exceptionally high.The significant damage throughout this time is observed at the mounting holes of the accumulator while the rest of the accumulator sustained low to no damage at all.

Discussion
From most of the results, we comprehend that there is still much scope for further improving the results.Our state of charge curve shows that only 20% of the charge is left in the cell at the end of the drive cycle run.Implementing a regeneration concept and

Lifespan
Steinberg's Formula (Equation (11)) is used to calculate life via the fatigue tool in random vibrations in ANSYS.Steinberg formulation utilizes all three stress occurrences (1σ, 2σ, 3σ) and their rate of occurrence along with Miner's rule to compute fatigue damage of the system.
where, n 1σ = actual number o f cycles at or below 1σ level n 2σ = actual number o f cycles at or below 2σ level n 3σ = actual number o f cycles at or below 3σ level n 1σ , n 2σ and n 3σ = an Allowable number o f cycles ( f rom f atigue curve) at 1σ, 2σ, and 3σ stress levels The analysis resulted in a total lifespan of 1.047 × 1026 h, which is exceptionally high.The significant damage throughout this time is observed at the mounting holes of the accumulator while the rest of the accumulator sustained low to no damage at all.

Discussion
From most of the results, we comprehend that there is still much scope for further improving the results.Our state of charge curve shows that only 20% of the charge is left in the cell at the end of the drive cycle run.Implementing a regeneration concept and conjugating its model with our existing mathematical model can further help us bring down the battery pack capacity by 1.5-2 Ah, depending upon the level of regeneration we want to achieve.A detailed study of the charging and discharging characteristics of the cell can help develop its charging technology, and the independent SOC study of each cell can help increase the lifespan of these energy storage devices.Physical testing of the pouch cell on cell testing rigs can help us establish a Battery Management System (BMS).A BMS enables and assists in keeping the charge, voltage, and current balanced and uniform among all the energy device cells.This concept can help downsize our accumulator by heightening the SOC limits towards the lower ends.Our mathematical model can further incorporate a heat generation model so that heat rejection within the energy storage device can be simulated and analyzed to improvise the accumulator's material selection, orientation, and configuration.Software like IPG carmaker and Carsim can further model the intricacies and integration of the energy storage device with the rest of the car and estimate the differences in lap times and overall efficiency.Generative design can optimize the topology of the energy system containers making them more reliable and cost-effective for alternative applications.

Result
The final accumulator model was successfully designed against all the predetermined objectives of the project.With a final drive reduction of 5, the total lap time was determined as 1 min 45 s per lap.The resistive force value of 1224 N helped estimate the motor power required in the range of 20-30 kW (continuous), which helps in sizing the powertrain and justifies the selection of the Emrax 188 motor.The battery sizing models provide a 90S1P cell configuration for the selected LFP pouch cells; this configuration helps us safely use the accumulator for the desired period without running out of power with approximately 20% charge left at the end of the run.The accumulator module experiences a peak deformation of 0.56 mm along the longitudinal axis.The accumulator brackets generate peak Von Mises stress of 150 MPa and experience a peak deformation of 0.18 mm.The accumulator container experiences a peak deformation of 3 mm.All the stresses generated in the container are below the ultimate tensile stress for the material, with a safety factor ranging from 0.1 to 0.5.For the vibrational analysis, the total deformation experienced by the accumulator is 3.5493 mm at 1760.8 Hz.The accumulator design is hence safe to use in electric Formula Student racecars.

Conclusions
The methodology described in the paper can be used to design the accumulators for passenger cars and public transport.Even though each application is different and will follow other safety and design requirements, the underlying design methodology used in this paper can be applied.The driving conditions and the driver psychology while driving different segment vehicles are very diverse; hence by interpreting these attributes carefully, designers can tune their electric motors with suitable battery performance maps to get close to the optimal driving experience.This idea can be extended to all X-EVs and other electric mobility platforms by modifying them as needed.

Funding:
The authors extend their appreciation to King Saud University for funding this work through Researchers Supporting Project number (RSP-2021/164), King Saud University, Riyadh, Saudi Arabia.

Figure 2 .
Figure 2. The output of lap-time simulation from optimum lap shows the speed at each instant the track.Note that due to the lack of long straight lines the top speed of the vehicle is around kmph.

Figure 2 .
Figure 2. The output of lap-time simulation from optimum lap shows the speed at each instant on the track.Note that due to the lack of long straight lines the top speed of the vehicle is around 77 kmph.

Figure 5 .
Figure 5. (a) Motor speed.(b) Motor torque for the entire run.

Figure 6 .
Figure 6.Motor power throughout the entire run.

Figure 6 .
Figure 6.Motor power throughout the entire run.

Figure 6 .
Figure 6.Motor power throughout the entire run.

Figure 5 .
Figure 5. (a) Motor speed.(b) Motor torque for the entire run.

Figure 6 .
Figure 6.Motor power throughout the entire run.

Figure 9 .
Figure 9. (a) Battery discharge C-rate and (b) SOC for the entire run.

Figure 10 .
Figure 10.Cell connection using cell tabs.(a) Cell connection in series.Double tab washers are used to make sure that the bolts do not loosen due to vibration; (b) connections of cell tabs are made by bending the cell tabs over a standard bus bar.This eliminates the need to punch holes in the cell tabs, which might damage the cell tabs.

Figure 10 .Figure 11 .
Figure 10.Cell connection using cell tabs.(a) Cell connection in series.Double tab washers are used to make sure that the bolts do not loosen due to vibration; (b) connections of cell tabs are made by bending the cell tabs over a standard bus bar.This eliminates the need to punch holes in the cell tabs, which might damage the cell tabs.Energies 2021, 14, x FOR PEER REVIEW 13 of 22

Figure 11 .
Figure 11.(a) Accumulator container without the modules.(b) Assembled accumulator container with the modules.

Figure 12 .Figure 13 .
Figure 12.(a) Longitudinal deformation.(b) Lateral deformation.(c) Vertical deformation.The base of the modules is considered as fixed points and load is applied to the module walls.

Figure 12 .
Figure 12.(a) Longitudinal deformation.(b) Lateral deformation.(c) Vertical deformation.The base of the modules is considered as fixed points and load is applied to the module walls.

•
40 g in the longitudinal direction (forward/aft), • 40 g in the lateral direction (left/right), • 20 g in the vertical direction (up/down), • Gussets attaching the accumulator to the chassis should withstand the load of 20 kN in all directions.

Figure 12 .Figure 13 .
Figure 12.(a) Longitudinal deformation.(b) Lateral deformation.(c) Vertical deformation.The base of the modules is considered as fixed points and load is applied to the module walls.

Figure 13 . 22 Figure 14 .
Figure 13.(a) Equivalent Von Mises stress for the bracket (MPa).(b) Total deformation (mm) for the bracket.The base of the brackets is taken as fixed support, and the forces are applied to the walls of the brackets.Energies 2021, 14, x FOR PEER REVIEW 15 of 22

Figure 15 .
Figure 15.Static structural simulation due to cell weight under gravity results in negligible deformation of 2.47 × 10 −5 mm on the modules.

Figure 15 .
Figure 15.Static structural simulation due to cell weight under gravity results in negligible deformation of 2.47 × 10 −5 mm on the modules.Energies 2021, 14, x FOR PEER REVIEW 17 of 22

Figure 16 .
Figure 16.Mode (X-axis) v/s frequency.Each mode has a respective frequency for analysis.

Figure 16 .
Figure 16.Mode (X-axis) v/s frequency.Each mode has a respective frequency for analysis.

Figure 16 .Figure 17 .
Figure 16.Mode (X-axis) v/s frequency.Each mode has a respective frequency for analysis.The total deformation is calculated for modes 1-650, ranging from 17.895 to 1767.3 Hz.The minimum to maximum total deformation experienced by the accumulator ranges from 0.4890 to 3.5493 mm for 17.85-1760.8Hz, respectively.The deformation results can be seen in Figure17.

Figure 17 .
Figure 17.(a) Minimum deformation for modal analysis at 17.895 Hz frequency is 0.489 mm.(b) Maximum deformation for modal analysis 1760.8 Hz frequency is 3.5 mm.

Energies 2021 ,
14,  x FOR PEER REVIEW 18 of 22 sors to the setup of random vibrations.The power spectral density acceleration graph is defined in the setup.The ASD (G2/Hz) V/s frequency (Hz) graph is shown in Figure18, and the input values are shown in the table.The software will hence interpolate the values for G2/Hz for the frequencies for each mode.

Energies 2021 , 22 Figure 19 .
Figure 19.Directional deformation (mm) for 3 Sigma scale factor.The maximum deformation is 0.45 mm towards the top of the accumulator module.2.6.4.Lifespan Steinberg's Formula (Equation (11)) is used to calculate life via the fatigue tool in random vibrations in ANSYS.Steinberg formulation utilizes all three stress occurrences (1σ, 2σ, 3σ) and their rate of occurrence along with Miner's rule to compute fatigue damage of the system.

Figure 19 .
Figure 19.Directional deformation (mm) for 3 Sigma scale factor.The maximum deformation is 0.45 mm towards the top of the accumulator module.

Table 2 .
Model input: the inputs used to define the characteristics of the car.

Table 3 .
Datasheet of A123 cells considered for the design.

Table 3 .
Datasheet of A123 cells considered for the design.

Table 4 .
Scale factor values, probability of occurrence, and maximum directional deformation result for random vibration analysis.

Table 4 .
Scale factor values, probability of occurrence, and maximum directional deformation result for random vibration analysis.