This contribution shows an analysis of vibration measurement on large floor-mounted traction batteries of Battery Electric Vehicles (BEV). The focus lies on the requirements for a realistic replication of the mechanical environments in a testing laboratory. Especially the analysis on global bending transfer functions and local corner bending coherence indicate that neither a fully stiff fixation of the battery nor a completely independent movement on the four corners yields a realistic and conservative test scenario. The contribution will further show what implication these findings have on future vibration & shock testing equipment for large traction batteries. Additionally, it will cover an outlook on how vibration behavior of highly integrated approaches (cell2car) changes the mechanical loads on the cells.
Shock and Vibration Testing is a crucial task for evaluating performance characteristics and ensuring safety functions and mechanical environments in nearly any mobile application. This is particularly true for Lithium-Ion based Rechargeable Energy Storage Systems (RESS), which are treated under the UN Dangerous Goods Convention and need to undergo at least a risk mitigation scheme for transportation environments described in UN38.3 [
Hence, the safety critical battery components are typically tested on various levels from cell, over module to battery pack/traction battery. Tests on vehicle level are then finally undertaken for high-level verification and validation.
As current BEV do not show a higher risk of catching fire than current Internal Combustions Cars (ICE) [
This leads to the important aspects of the location of a component in a vehicle and the interaction between vehicle and battery pack. Typically, performance and safety testing are not only done on the full product level (in this case a complete real life production car) but taken down to lower component level in the V-model [
The common approach of industry standards like ISO 16750-3 reflects this by presenting some generalized local environments for components attached in the respective region. Different tests are described for components attached to the vehicle body or to unsprung masses like the wheel suspension, for example. For large RESS battery packs or traction batteries weighting several 100 kg, this approach does not work out well anymore [
Large, heavy floor batteries with significant interaction between battery and vehicle are treated in Category 3 of ISO 19453-6 regarding the shock and vibration testing requirements. Different from the other categories of less interaction, no generalized test level for the mechanical environment is given for Category 3. Instead, the standard indicates the need for “vehicle specific profiles and testing time” or “vehicle specific time signals” for testing. This contradicts the goals of a rapid V-model design process because no trustworthy design requirements and test levels are given on the general normative level.
Here further research is needed to derive partially generalized test levels and requirements depending on the stiffness/mass interaction between battery and vehicle frame.
Further need for research arises with of a four-point-excitation of a generalized large floor battery as in [
Another interesting aspect for vibration behavior assessment is the change of mechanical properties of electrical and thermal aging. Not only the volume of a Lithium-Ion battery changes over temperature, state of charge and cycles but also key figures in the elastic and plastic behavior change, as described in [
With the given state of the art, several questions arise to fill existing gaps and make further progress in efficiency and quality in the development of lithium-ion battery packs under mechanical vibration & shock environments, such as:
This question concerns the relevant excitation frequency range of up to 500 Hz and in special cases up to 2000 Hz. Resonances with an amplification of more than Factor 2 are typically considered relevant, as the amplifications in acceleration cause an equivalent amplification in forces and stresses based on the modal mass fraction of the bending mode.
In a common road vehicle, the interaction between vehicle and ground causes a power input through tire and suspensions at the four corner points of the floor section between the wheels. If this floor section housing large batteries in modern BEVs shows considerable large deflections between the four corners and the mass center, local corner bending occurs and causes bending stresses.
A higher integration level will bring more coupling and eventually higher loads on the subcomponents like cells or their direct housings. Measurements on the bending behavior with various added loads for an integral structure as well as for decoupled battery packs in a differential approach should give insight into answering that question.
The following chapter shortly describes the background and methods in vehicle dynamics as well as signal analysis used in this paper.
Vehicle models with component movement or deformation typically range up to 200 Hz [
A larger stiff Nickel Manganese Cobalt Oxide Pouch Cell showed resonance frequencies of around 200 and 300 Hz, which is higher than typical road-induced vibrations of sprung masses in a passenger car. Internal vibration excitation may still originate from the electric motor but on a much smaller level than on an ICE.
Therefore, the analysis focus of this investigation is set to 0.5–2000 Hz with a particular focus of up to 200 Hz as long higher frequencies power spectrums remain at least one decade lower.
Rigid body movements of ground vehicles are widely described in literature, as, i.e., [
In [
A description of general elastic deformation types of a simplified large traction battery in the floor segment of a vehicle body is described in
It is easily understandable that these deformations and movements cause local stresses and deformations in the battery structure and cells. The deformations are caused by inertial forces as well as incoherent forces acting on the corners.
Local and global resonances will amplify the deformation and induced stresses depending on the damping of the corresponding mode shape. In order to describe the first relevant global and local corner bending mode shapes, the structure is described with translational movement vectors at five points, consisting of four corner and one point in the middle, see
The model can be further reduced to vertical movements at the given points, especially for the small bending deflections and local deflections to be expected. All points have a given distance in the horizontal plane (X,Y) from each other. The distance in X, Y components make up the distance values of DX… and DY… All points can make a vertical movement over time. So, at a given time step t each point will have always the same X and Y distance from each other (this simplification is only acceptable for small movements). But each point may have a different
Rotational movements at the points would be of further interest but are more difficult to measure compared to the small piezo-electric accelerometers for uniaxial measurements. Measurements with 6dof (3 translational, 3 rotational) would be possible to perform with MEMS. As this seems to be difficult to achieve together with translational movements at high frequencies and was not the main scope presently, it is left for future research.
By introducing a given distance between the points, a bending effect is indirectly described by the comparison of a completely stiff rigid body movement of all points and the actual movement at the points measured.
This leads us to the task of analysis methods in order to quantify local bending influence at the four corner points as well as global bending mode around X and Y with the midpoint near the assumed maximum of the mode shape. The further analysis will use methods in time and frequency domain, both based on the modal domain models of local corner and global bending. Further information on the different domains of dynamic analysis can be found in [
A quick method to describe the bending influence inducing stress in the traction battery for an early analysis is to look at the deflection difference at a given corner between its rigid body movement and the actual measured deflection from the vertical position (see
In a simple battery model with vertical movements at the four corners and at one point in the middle (see Global bending can be detected for first simple modes but not distinguished clearly by looking at the averaged movements of the four corner points and the point in the middle; Global torsion in its first simple modes and local corner bending can be detected but not distinguished from each other by looking at the difference in movement of one of the four corner points and the rigid body movements based on the other corner points.
As this contribution targets a first holistic view on mechanical battery loads, the simplification with five respectively four points is deemed acceptable, as it gives enough information for the intended purpose. For a more detailed view on local stresses in the structure of a fully detailed design more points of movements have to be considered.
With the model simplified as described above, one can calculate a rigid body movement as shown in the following example of Point 3 from the measured vertical movements. From three out of the four corner points the assumed rigid body movement of the one left out is calculated:
For
where
The flex movement of first torsion modes and local corner bending in time domain difference is of particular interest for later statistical analysis:
The same needs to be done for all other corner points with their opposite corner coming in negative as in the above equation. Basically, this approach detects modes with local or global torsion parts. The simple global bending modes are not detected by this approach but will be by the approach in the following sub-chapter. Higher Torsion modes with local minima and maxima between the four corner points cannot be distinguished from each other.
Remark: This procedure may overestimate local corner bending by putting all non-rigid flexure of the body into the difference at one point. Opposite points will yield the same values of
When looking for a simple indication of global bending in time domain, especially the difference function between the mean of the vertical movements of all four corner points
Remark: This procedure highly depends on putting the Point 5 in the global bending maxima of the mode, which may be unlikely for higher modes and more complex mode shapes. The approach was used as a simple and quick qualitative check if there is any significant global bending in relation to the actual measured values to be detected in the measurement files. For this it was rendered sufficient.
Different statistical key figures can be used for comparison on the differential movements in time domain. The statistical key figures deemed most relevant for the task given are: Root Mean Square (RMS)-to describe the equivalent DC-mean and (if squared) a general “power” content of the signal. Kurtosis-method to describe the statistical distribution of a signal in one number instead of looking on the “pointiness” of a probability or distribution chart over classes of the analyzed signal. In our case, a high kurtosis would describe a signal with more and higher extreme signal peaks as well as more measurement data around zero compared to a gaussian signal of the same RMS level. A Gaussian distribution has Kurtosis 3. Skewness-describes the difference between the positive and the negative side of a signal similar to Kurtosis in one number. A high skewness indicates one-sided contacts in a vibration movement. In our case that could be a passenger door contacting to the frame in one direction of the vertical movement.
Further descriptions on the use of these statistical measures in vibration analysis can be found in [
Apart from the simple and quick time domain analysis methods described above, transfer functions in frequency domain between different points (real measurement and calculated) were also used.
When looking for information on frequency domain information like resonance frequencies of the assumed bending modes or amplitude and phase, the frequency domain transfer function is calculated for any circular frequency
The transfer function to describe the flexible global bending behaviour of the Midpoint 5 is denominated by capital letters and calculated in the frequency domain by the Fast Fourier Transforms (FFT) of the calculated time domain rigid
The transfer function then needs some averaging or smoothing, especially for a high number of data inputs at a high sampling rate resulting in good frequency resolution over a broad frequency range. The following uses an arithmetic average over 10 time-segments of each measurement file if not mentioned otherwise. A suitable windowing and amplitude compensation also need to be considered. For the later often used measurements on rough cobble stone a typical measurement includes data of ca. 60 s at 15.000 Hz sampling rate giving something around 900.000 data points for each sensor. For 10 equal, non-overlapping time segments, this would render 90.000 data points of which not all may be analyzed due to FFT restrictions. Further information on transfer function generation can be found in the chapter “experimental modal analysis” in [
The transfer functions between the mean of all corner movements and the measured midpoint were analyzed for various mass loading configurations of the vehicles. As modal stiffness is assumed not to change, the added mass should lower the resonance frequency of a mode, as long as the structure the mass is added to is coupled to the resonating structure. The results indicate to what degree vehicle structure and battery pack are mechanically coupled in the later result section. This idea is later used in the analysis by repeating the same road segment measurements on the same vehicle with different amounts of added mass in form of mass loading into the passenger compartment of the vehicle. Obviously, the mass loading needs to be high enough with regard to the mass fraction already resonating.
As most mechanical influence is related to the mechanical power input, the Power Spectral Density analysis method is a very helpful tool. It is not only based on the measured acceleration levels but it enables a comparison of the power input over frequency. Basically, the acceleration input over a small frequency range (i.e., like 1 Hz wide) is squared and divided by the width of frequency range. As this contribution covers a first look at requirements for shaker testing of large floor batteries and the data is still relatively small with a few dozen measurements on two different BEVs, no PSD curves are presented here. For a proper test specification much more information on typical usage load blocks and possible fatigue damage mechanisms have to be analyzed and considered. This will be covered in future publications as research progresses.
With a larger set of measurements for more different vehicles, events and loading configurations, future research will also include an analysis based on the Fatigue Damage Spectrum (FDS). It benchmarks the occurred vibration peaks against a fictive S-N-curve to compare the damage potentially induced. This analysis is related to a Shock Response Spectrum (SRS), which is used to benchmark the measured signal regarding the maximum excitation of a fictive vibration system, see also [
In order to answer the question of how independent the four corners or the center of the batteries move, a frequency-dependent measure of coherence between two signals is needed. Here, the ordinary coherence function is used. It is calculated by the squared cross-spectrum over the product of the two auto-spectra of the two signals for each frequency point. The coherence tells the linear relationship over frequency for the two signals. Further information can be found in the chapter frequency response estimation of experimental modal analysis in [
The current contribution focusses on an analysis of the needs on how to enable a realistic environment for large floor tractions batteries in a lab regarding its real-world excitations. This analysis is based on the first results of a smaller pre-test measurement campaign, with more planned in the future. These preliminary results indicate significant changes on the approach of component vibration testing in comparison to the state of the art.
The data for analysis was measured on serial production cars in real world environments on different road surfaces under varied loading configurations. Each measurement is defined by vehicle, road surface, loading configuration and repetition. The road surfaces include but are not limited to: pot hole–30 kph, cobble stone–30 kph, gully–30 kph, city drive HH-various speeds, country road–100 kph, motorway–130 kph, country roads-various speeds, rough cobble stone–30 kph. Furthermore, the noise floor was measured without engine running (for ICE) and without movement (both BEVs and ICE).
The vehicles used for this analysis are shown in
Vehicle | Compact BEV |
Compact BEV |
Large ICE |
---|---|---|---|
Empty weight | 1320 kg | 1810 kg | 2250 kg |
Max gross weight | 1670 kg | 2270 kg | 3000 kg |
Added mass loading in addition to driver, 80 kg | 87.5 kg |
200 kg | 200 kg |
Battery energy | 33 kWh | 62 kWh | n.a. |
The table shows empty and max gross weight of the vehicles as well as information on the additional loading in the measurement campaign. Beside the mass loading of the driver and minor other masses, the additional masses loaded are given in the table. For VW ID3 and VW T5, the loaded configuration included 200 kg of extra mass, the BMW i3 was measured in the test campaign in two different extra loading configurations. All vehicles were also measured in empty configuration (including driver, etc., with around 80 kg).
The large ICE was used as a general example to analyze how a larger mass rather stiffly coupled to the vehicle floor will change vibration behavior of the overall structure. This will give a first insight in vibrational behavior of highly coupled battery and vehicle interaction as in a cell2car approach.
The measurement and data acquisition equipment used is described in
Both MEMS and IEPE/ICP sensors were used to measure the acceleration at the same location with the sensors being placed as closely together as possible. This comparison allows a verification of the measurement procedure and equipment used at this early pre-measurement campaign. The results of the comparison are not part of the scope of this contribution. Generally, apart from the minor deviations or particular disadvantages in some aspects mentioned in the following, both yielded the same results for the scope of this analysis. Of course, the general aspect of giving good data quality with both sensor types was checked on a shaker setup prior to the measurement campaign. As a pre-test verification, the parallel measurement was also carried out in the actual vehicle measurements. Mixing, combining or comparing results from both measurement types were avoided for the scope of this contribution.
The MEMS accelerometers were mainly used to get high data quality especially for low frequency range below 10 Hz. The sensors on the given evaluation boards used here in particular show a decay of voltage output for higher frequencies at constant acceleration. Therefore, a linear compensation function between 20 Hz (Amplification 1) and 600 Hz (Amplification 6) was used, which compensated the non-linear frequency behavior of up to roughly 500 Hz quite well with max 6% error. The frequency range below 20 Hz did not need any compensation on the MEMS and yields more precise values than the piezo accelerometers. As most analysis was done either on low-frequency time domain signals or on frequency domain transfer functions between two signals, the uncompensated MEMS results can be used equivalently to the piezo accelerometer data between 10–500 Hz.
For the piezo accelerometers a frequency linearity of up to 5% is specified from 1 Hz on. But in most cases the measured accelerations on ICP are several decades below the maximum measurement range, so that the ICP output voltages are a little noise-prone and need a lot of signal conditioning. Therefore, the higher output voltage per acceleration of the MEMS used here is assumed to yield slightly better results.
The MEMS and IEPE accelerometers were measured with two separate but identical USB data acquisition boxes. As no combination of signals from both measurement setups was used, the general problem of real-time synchronization between two sampling systems is not relevant here.
When looking at signals with particular high kurtosis it was striking that several ICP sensors showed a kurtosis of around 100, while the MEMS sensors in the same places showed values between 3 and 15, i.e., in a pot hole event with the same loaded vehicle. From a closer view on the signals, it is assumed that the screwed connectors of the ICP sensors could have caused short contact losses. On the MEMS sensors the cables are soldered on the printed circuit board directly. However, as any rattling on the vehicle generates significant high frequency input, the MEMS sensors with a decaying frequency sensitivity above 500 Hz will also measure less non-linear contact input prone to cause high kurtosis values on the ICP sensors that can measure much higher frequencies.
Further signal checks were performed based on the comparison of RMS
Most of the analysis is done with National Instruments DIADEM 2020. This includes digital filtering, frequency domain compensation, Fast Fourier Transformation (FFT), Power Spectral Densities (PSD), transfer functions with amplitude and unwrapped phase, coherence functions, Cross Spectral Densities (CSD), channel arithmetic for time domain difference functions, statistical distributions (random probe) and statistical key figures like Root Mean Square (RMS), min, max, kurtosis and skewness.
Obviously, many findings can be derived of measurements from the first 66 measurements on at least eight channels with the three different cars in 2–3 loading conditions, over 4–6 different events/road conditions with one or two independent measurement systems (MEMS and ICP) amounting to. The main focus of the analysis in this contribution lies on the three questions in Does relevant global bending occur? Does relevant local corner bending occur? What changes will a more integrated cell placement like cell2car bring in comparison with a more differential decoupled pack approach?
From these questions an adequate testing solution in concept, control and hardware can be developed.
Remarks on the general signal quality and measured vehicle vibration behavior: Both measurement procedures with MEMS and ICP sensors lead to very useful, comparable, consistent data even though the MEMS data needed a frequency band compensation to achieve frequency linearity up to 500 Hz. This can be circumvented by using other or modified versions of the MEMS measurement equipment used. Most driving conditions-including cobble stone and rough cobble stone-show a kurtosis around 3 (gaussian) on all vertical floor segment corners and axle sensor, both loaded and empty. The transient gully and pot hole events show a very high kurtosis >>20 for vertical axle and an elevated kurtosis for vertical floor segment corners, mostly between 5–20. In cases of high and elevated kurtosis on the signals: If the acceleration data of time domain difference functions between the different point on the floor segment or battery or difference functions between calculated mean, calculated rigid body moment equivalents and the measured counterparts is used for analysis (i.e., for calculating bending deflection) it needs to be taken account that large differences in the kurtosis of the signals subtracted may lead to a significantly higher kurtosis of the difference function, due to a mean offset. Sometimes the value may reach twice the kurtosis of one of original signals. It is not used for analysis.
Some data on the statistical key figures of two different measurements can be found in
The following syntax is used for the generalized points on the chassis and pack:
Be: Bearing between suspension and chassis (input into vehicle body); Ax: Axle of suspension (unsprung masses); RR: right rear; RF: right front; LF: Left front; LR: Left rear; X, Y, Z: directions in vehicle coordinate system; AC: Alternating values without stationary mean.
Channel | Minimum | Maximum | RMS | Skewness | Kurtosis |
---|---|---|---|---|---|
Be_RF_X AC | –0.296 | 0.204 | 0.0603 | –0.320 | 2.96 |
Be_RF_Z AC | –0.526 | 0.752 | 0.0955 | 0.644 | 12.43 |
Ax_RF_Z AC | –7.12 | 8.28 | 0.829 | 0.959 | 44.68 |
RR_Z AC | –1.25 | 0.602 | 0.0811 | 0.413 | 12.67 |
RF_Z AC | –0.543 | 0.725 | 0.0922 | 0.865 | 14.09 |
LF_Z AC | –0.383 | 0.227 | 0.0539 | –1.06 | 8.27 |
LR_Z AC | –0.365 | 0.235 | 0.0550 | –0.893 | 7.79 |
C_Z AC | –0.299 | 0.270 | 0.0465 | –0.391 | 7.63 |
RR diff of mean | –0.978 | 0.457 | 0.0552 | ||
RF diff of mean | –0.333 | 0.542 | 0.0642 | ||
LF diff of mean | –0.496 | 0.339 | 0.0568 | ||
LR diff of mean | –0.534 | 0.343 | 0.0647 | ||
C diff of mean | –0.206 | 0.394 | 0.0172 | ||
RR flex (diff local rigid) | –1.28 | 0.517 | 0.0476 | ||
RF flex (diff local rigid) | –0.517 | 1.28 | 0.0476 | ||
LF flex (diff local rigid) | –1.28 | 0.517 | 0.0476 | ||
LR flex (diff local rigid) | –0.517 | 1.28 | 0.0476 |
Note: unit: [g], approx. 9.81 [m/s²].
Channel | Minimum | Maximum | RMS | Skewness | Kurtosis |
---|---|---|---|---|---|
Be_RF_X AC | –0.628 | 0.414 | –0.128 | 0.177 | 3.30 |
Be_RF_Z AC | –2.61 | 0.912 | 0.0164 | –0.0254 | 3.11 |
Ax_RF_Z AC | –9.02 | 8.22 | –7.64E-03 | –0.0654 | 4.25 |
RR_Z AC | –1.61 | 0.807 | –2.15E-03 | –0.0780 | 3.45 |
RF_Z AC | –0.699 | 0.770 | 1.02E-03 | –0.0332 | 3.01 |
LF_Z AC | –0.701 | 0.711 | 4.03E-03 | –3.04E-03 | 2.97 |
LR_Z AC | –0.722 | 0.655 | –2.03E-03 | –0.0538 | 2.93 |
C_Z AC | –0.485 | 0.451 | 2.56E-03 | –0.0923 | 2.91 |
RR diff of mean | –1.16 | 0.609 | –2.36E-03 | ||
RF diff of mean | –0.563 | 0.823 | 8.03E-04 | ||
LF diff of mean | –0.587 | 0.604 | 3.81E-03 | ||
LR diff of mean | –0.604 | 0.613 | –2.25E-03 | ||
C diff of mean | –0.234 | 0.502 | 2.34E-03 | ||
RR flex (diff local rigid) | –1.65 | 0.696 | 2.90E-03 | ||
RF flex (diff local rigid) | –0.696 | 1.65 | –2.90E-03 | ||
LF flex (diff local rigid) | –1.65 | 0.696 | 2.90E-03 | ||
LR flex (diff local rigid) | –0.696 | 1.65 | –2.90E–03 |
Note: unit: [g], approx. 9.81 [m/s²].
The quick pre-analysis data check based on the time-domain difference functions clearly shows that relevant bending or torsion effects occur:
In both measurement examples a significant bending can be assumed from comparing the min, max and RMS of the center difference between the movement mean of all four corner points and the measurement point in the middle denominated as “C diff of mean”, amounting to more than 50% of the corresponding values of the simple measurement values of “C_Z AC”, the midpoint of the battery.
Also, both measurement examples indicate a local corner bending or first torsional modes from looking at the local difference functions between calculated rigid body movements and actual measured movements, denominated as “XX flex (diff local rigid)”, also more than 50% of the original measurement files.
For a more detailed analysis than the simple time-domain difference functions for a quick data check, the frequency domain approach is used. The global bending with the highest deflection to be expected in the middle of the floor segment is characterized by the stiff/mass ratio, which defines the resonance frequency for a single degree-of-freedom vibration (SDOF) system. Obviously, a battery model may need higher mode numbers for an adequate description, but the following general relation remains the same, also for more SDOF systems: As the stiffness of the vehicle body and chassis does not change during use, an added mass loading will change the relevant resonance frequency of the structure. This approach was used to analyze the stiffness behavior of the floor segment of an ICE as well as the battery pack in the floor region of BEVs, see
Both the VW ID3 and BMW i3 show significant global bending modes. The VW ID3 has significant global bending resonances at around 60 and 180 Hz with amplifications of 10 and 3. The BMW i3 has significant global bending at around 40 and 110 Hz with amplifications of 10 and 3.
Both BEV show a totally decoupled behavior, because the added mass does not lower the resonance frequency of the battery’s center. The battery pack in both the BMW i3 and VW ID3 is a very stiff block that is mounted into a mostly stiff battery compartment with the load carrying vehicle frame around. The battery pack seems to be well decoupled from the vehicle frame (body & chassis) carrying the added loading, both parts being rather stiff. The mechanical interaction of the battery with the load carrying vehicle structure is quite low.
Global damping properties cannot be derived really accurately due to only one center measurement point in the middle with no verification that it was mounted in the global mode shape maximum. But an amplification of around Factor 10 seen in the measurement plots indicates a damping ratio of at least 5%.
Regarding artificially generated vibration environments in a lab it is noted that highly coherent excitations in amplitude and phase at the four corners (i.e., all move up and down synchronously) yield higher amplification ratios in bending resonances. Therefore, it is of importance that the coherence of the four excitation points over frequency on a lab test stand reflects the real-world excitation coherence. Completely interlocked excitation of the four corner excitation points (
The influence of corner excitation coherence changes when looking at the local corner bending.
In this case a highly incoherent excitation in amplitude and phase at corners (
Together with the opposite effect in global bending this leads to the necessity of a close resemblance of the corner excitation coherence in the lab to the measured signals. An example comparison between the large ICE and the compact BMW i3 with stiff decoupled battery pack is shown in
The mid-level coherence of the plots reflects the strong interaction of the corners with load carrying vehicle frame (in contrast to the global center bending).
As an indicator for the trend towards higher integration levels (cell2car) a large ICE (VW T5) was measured with and without added mass in the middle of the floor segment, too. For the large ICE, the added mass in the middle of the floor segment significantly lowers two resonance frequencies in the transfer function from calculated rigid center to measured center, see
These two extremes are shown, because current developments of battery packs seem to move towards more integrated, load-carrying battery structure and away from a stiff, decoupled battery “block” pack as in the BMW i3 and VW ID3.
This early analysis of a pre-test measurement campaign is based only on measurements of two BEV and one ICE. Therefore, the results cannot be generalized regarding the given quantitative information. So far, the effects found relevant for future testing requirements can be only generalized on a qualitative scale. Additionally, a quantifiable analysis of a broader data base will use more detailed methods like FDS and Shock Response Spectrum (SRS). Therefore, no PSD plots for potential random vibration testing environments are given in this contribution as more test results are needed to derive a proper test specification regarding vehicle life and fatigue damage effects.
Generally global resonance bending is important because of multiplying stresses in the structures, especially for flexible heavy structures like large battery packs. The measurements show that all vehicles have significant global bending of amplification Factors 10 or more in the frequency range of up to 200 Hz. Even though the measurement setup with only a few points does not allow a high precision, the damping ratio of global bending modes can be assumed to be ca. 5% from amplification up to Factor 10. High modal mass fractions are likely for main global bending modes, indicating a general potential of high stresses.
The comparison of a very stiff battery pack more or less decoupled from a rather stiff frame as in VW ID3 and especially in BMW i3 with the example of highly coupled loading in a common ICE transport vehicle indicates what can be expected in a more integrated structure with a stronger coupling. Generally, the interaction of large battery masses with the vehicle structure is likely to change vibration and shock environments on the battery pack. The more integrated a cell2car approach gets and the less levels of supporting structures between cell and vehicle it has, the more bending and torsion will cause stresses on cell level.
The data shows that the forced local bending at corners is mid-level correlated with 20%–50% coherence in lower frequency band. Slightly higher coherence values can be found for stiffer structure (edge frequency between 50 Hz for lower stiffness, 80 Hz for higher stiffness). Generally, forced corner bending was found up to at least 500 Hz.
From the presented data and the resulting analysis some requirements for a component testing approach of large floor-mounted traction battery packs can be derived:
The vibration testing needs an excitation and setup that allows conservative equivalents or replications of real mechanical environments to cause global bending in resonances. Furthermore, an individual excitation of the four corners of a battery traction pack is important to replicate the real world’s coherence. However, the right mechanical impedance of test setup and interfaces is crucial: A very stiff fixture with identical excitation on all four corners potentially overtests global bending and undertests local corner bending. A completely independent movement of all four corner points potentially overtests local corner bending and undertests global bending. A matched coherence between lab and real-world excitations would solve this problem. This may be done mechanically or by control methods.
The need for further research is especially seen in creating a broader data base from more measurement campaigns with more data in higher spatial and frequency resolution, deriving simple models to describe the interaction between vehicle and battery pack so that vibration levels design requirements and for test specifications can be derived during the development process of the vehicle and not only when everything is finished, developing a suitable test infrastructure with multi-point excitation and interfaces that can replicate realistic interaction with vehicle and excitation.