Methodology for Comprehensive Comparison of Energy Harvesting Shock Absorber Systems

<jats:p>In recent years, there has been a lot of work related to Energy Harvesting Shock Absorbers (EHSA). These devices harvest energy from the movement of the vehicle’s shock absorbers caused by road roughness, braking, acceleration and turning. There are different technologies that can be used in these systems, but it is not clear which would be the best option if you want to replace a conventional shock absorber with an EHSA. This article presents a methodology to compare the performance of different EHSA technologies that can replace a shock absorber with a given damping coefficient. The methodology allows to include different analysis options, including different types of driving cycles, computer vehicle models, input signals and road types. The article tests the methodology in selecting the optimal EHSA technology for a particular shock absorber and vehicle, optimizing at the same time energy recovery. In addition, a study of parameters in each type of technology is included to analyze its influence on the final objective. In the example analyzed, the EHSA technology with a rack and pinion system turned out to be the best. The proposed methodology can be extrapolated to other case studies and design objectives.</jats:p>

section_content: Given the paradigm of the world's increasing energy needs and environmental problems due to pollution, reducing energy consumption, or improving the efficiency of existing or new systems is a must. One of the sectors with the most opportunities for improvement is transport, which since 1980 has been the largest consumer of fuel, and also the largest generator of greenhouse gases. 

For light vehicles, regulations are stricter, and this has led to a trend to bring cars with higher fuel efficiency and lower emissions to the market [1]. This applies to electric, hybrid and internal combustion vehicles. Electric and hybrid vehicles are more efficient than internal combustion vehicles, with energy losses from the drive system of 15 to 20%, compared to 64 to 75% for the gasoline engine. In general, electric vehicles are 60-73% efficient, depending on the driving cycle [2]. 

The performance of cars can be improved by recovering the energy they waste. There are different technologies such as: Kinetic Energy Recovery Systems (KERS), Energy Harvesting Shock Absorber Systems (EHSA) and Waste Heat Recovery Units. 

The goal is to increase the effectiveness and life of electric and hybrid vehicle batteries to reduce high costs. Energy recovery in vehicles can improve their performance and cover different electrical power needs within the vehicle, from feeding back the vehicle battery to covering the needs of specific sensors without the need for additional wiring. 

Energies 2020, 13, 6110 3 of 25 64 W for the front shock absorber and 40 W for the rear one, the work of Lafarge et al. [17] with a 4 DOF model, recovering 40 W average. Models with 6 DOF that capture vertical displacement and pitch mode of both the sprung mass and the seat with an average recovered energy of 31. 63-55. 39 W have been reported by Nakano et al. [18]. Finally, full vehicle models, with 9 DOF and 24. 75 W recovery for 60 km/h are reported by Mi et al. [19] and 480 W for 90 km/h with 7 DOF by Pham et al. [20]. The work reported by Abdelkareem et al. [21] showed a model of 7 DOF with energy recovery of 51 W for a speed of 36 km/h, 100 W for a speed of 72 km/h, and 151 W with a speed of 108 km/h.
section_name: Introduction

section_content: In the simulation of the previous models, the researchers use different techniques to validate their developments. In some cases, simple signals are used. Sometimes, regulated driving cycles are selected. In other cases, signals recorded directly from vehicles are applied. The driving cycles are usually applied under standardized laboratory conditions to new vehicles before they enter a market. They consist of a sequence of events of braking and acceleration, attached to typical real road conditions. These tests are performed to determine the fuel consumption and CO 2 emissions of internal combustion engine vehicles [22, 23]. 

To evaluate the performance of energy recovery systems, sinusoidal functions [24, 25] can be used to random signals of variable amplitude to simulate standardized driving events from different environments and road types. Several driving cycles have been used in the energy recovery study of EHSA systems that provide better input information than harmonic inputs can. 

The use of drive cycles to evaluate EHSA performance began in 2016, using the New European Driving Cycle (NEDC) in a 3-degree of freedom model [26]. They used 10 different cycles, 5 of which were city cycles. The road roughness profile was included in a half vehicle model with 4 DOF [27] in the work of M. Abdelkareem et al. [21]. They conclude that steady-state driving cycles (composed of periods of constant acceleration and deceleration, NEDC, EUDC (Extra Urban Driving Cycle) at high speeds can recover more energy than transient state speed cycles (composed of a wide range of speeds and accelerations Worldwide Harmonised Light Vehicles Test Procedure (WLTP), FTP (Federal Test Procedure), HWFET (Highway Fuel Economy Test). They also add to the pitch movement a speed delay function between the rear and front wheels using a 7 DOF model. 

M. Abdelkareem et al. [28] test 10 driving cycles (HDUDD, NYCC, LA92, US06, SC06, CUEDC, IM240, J1015, RTS95, CSC) that were applied with random road profiles at different speeds and under two driving condition (aggressive and non-aggressive). They also include the loaded and unloaded car variant in a 6-degree of freedom half truck model. Within the 6 degrees of freedom, pitching is included. 

On the other hand, road profile modeling has been used extensively to simulate the roughness of road profiles and to obtain an estimate of the recoverable energy in them [29] [30] [31]. Road roughness has a fundamental influence on the energy that can be recovered in EHSAs. In theoretical roads without any roughness, the energy recovery would be limited to the pitching or turning movement of the vehicle. The amount of energy that is possible to recover is increased in roads with higher roughness as shown in the work performed by Lincoln Bowen and Yuxin Zhang [32, 33] where they compare different types of roads (A, B, C, and D) at different speeds and obtain the power. In the same way, mechanical efficiency is reduced by increasing the input frequency as shown by Zhongjie Li [9]. From these works, it can be concluded that the use of standardized driving cycles is a good approximation to estimate the energy recovered, which can be improved by including the vertical dynamics of the vehicle and the turns. 

To consider the EHSA energy recovery calculation of a complete vehicle it is important to ensure that pitch and roll capabilities are included in the vehicle model. This makes it interesting to include the lateral dynamics as well as the vertical and longitudinal dynamics of the vehicle. Published works on the subject sometimes include a trajectory profile delay function on each of the wheels, so that it causes rotations in the vehicle. If these delays are included in driving cycle tests (such as the WLTP), the simulation of recoverable energy would be closer to reality. 

In summary, after reviewing the articles published on the design, simulations, development and testing of EHSA systems, it can be concluded that there is a great heterogeneity of studies. There are different technologies applicable to EHSAs. In addition, the computer models used can include 2 DOFs or reach up to 9 DOFs. The models can be simulated in different types of roads with different roughness, at different speeds or different driving cycles, incorporating in the model the vertical and lateral dynamics. At the test level, they can be applied from basic sinusoidal signals to signals acquired directly from the vehicle instrumentation. This diversity in the design, simulation, testing and validation makes it almost impossible to have a clear idea of the performance of one EHSA system compared to the others. However, this comparison is essential if we are to take a leap forward and try to make this technology viable for the manufacture of functional prototypes. 

The motivation for this article is the difficulty that currently exists to contextualize and put into a common frame of reference all the data available in the relevant literature. From the point of view of an EHSA system designer it is desirable to be able to compare and know which technology is the best, and under which conditions it is so. It is also important to know which parameters are the most sensitive and most influential in the design. For this reason, this paper proposes a methodology that allows the comparison of results obtained through different technologies and EHSA designs. With this methodology, the comparative criteria of different EHSA prototypes for energy recovery are standardized, in a structured and agile way, based on different driving cycle conditions, different types of roads (both urban and highway), and with different driving parameters and vehicle characteristics. In this way, different EHSA systems can be compared with each other on the basis of a standardized working environment. This methodology for the analysis and comparison of EHSA systems aims at facilitating the design of an EHSA system to choose between the different available technologies and their design parameters, comparing them with each other under consistent and comparable validation conditions. 

The novelty of this article, therefore, aims to respond to this challenge by presenting an inclusive methodology and a unique integrative tool, which can take into account within it every one of the variables that have been discussed (vehicle models, types of EHSA, driving cycles, road profiles, etc. ) that the designer can choose to make the comparison. For this purpose, a new global methodology is proposed and its functionality is presented through a case study in which three of the technologies currently used in the design of rotational electromagnetic EHSA systems are compared. 

With the idea of being more inclusive in terms of the degrees of freedom that can be modelled, we will use a complete vehicle model that allows for the inclusion of movements such as pitching and rolling, movements that are normally experienced in real conditions. 

The article is distributed as follows. In Section 2, the methodology used in the work is shown, the simulation of the mathematical model is used with the incorporation of the WLTP driving cycle, going deeper into the modeling of turnings., In Section 3, the methodology applied is shown in the comparison of the different EHSA systems. Section 4 focuses on showing the results and discussions obtained from the three systems analyzed, and the parameters of greater effect selected, tested in the complete vehicle model with 7 DOF. Section 5 concludes from the results obtained and, finally, presents recommendations for future work.
section_name: Use of Driving Cycles in the Analysis of EHSA Systems

section_content: The methodology goes beyond establishing the same conditions for the modeling of all EHSA systems. In this way, four different stages were established (Figure 1 ). 

Input. In this first stage, the selected standardized driving cycles can be combined with the selected road roughness. The modeling includes the vehicle model and parameters such as: spring stiffness of the EHSA system, masses, and wheel stiffness. Finally, the user can select the preferred EHSA technology to be used. parameters of the EHSA system that have the greatest impact on energy harvesting. Once the inputs are available, the experiments or the vehicle model for the simulations can be performed. 

Analysis. At this stage, all the necessary iterations are carried out to obtain reliable, consistent, and congruent results. 

Validation. At this later step, it is possible to compare, under the same design and validation conditions, the different EHSA technologies selected for the design.
section_name: Materials and Methods

section_content: In this case, a full 7 DOF car model was chosen, which will be able to include in the energy recovery calculation vertical, pitching and rolling vehicle movements. The mathematical model is described in Figure 2. 

For demonstrative purposes, mass and rigidity parameters are simplified as follows: = = =, = = =. Experiments/Simulation. To reach this stage it is important to identify the most effective parameters of the EHSA system that have the greatest impact on energy harvesting. Once the inputs are available, the experiments or the vehicle model for the simulations can be performed. 

Analysis. At this stage, all the necessary iterations are carried out to obtain reliable, consistent, and congruent results. 

Validation. At this later step, it is possible to compare, under the same design and validation conditions, the different EHSA technologies selected for the design.
section_name: Suspension and Vehicle Models

section_content: In this case, a full 7 DOF car model was chosen, which will be able to include in the energy recovery calculation vertical, pitching and rolling vehicle movements. The mathematical model is described in Figure 2. 

∑ = on the sprung mass (pitch) For demonstrative purposes, mass and rigidity parameters are simplified as follows:

Equations for the complete vehicle are deducted from Newton's second law: ma = F at sprung mass

Iθ = T on the sprung mass (pitch)

I∅ = T in the sprung mass (roll)

To exemplify the car model, the parameters of a compact car was included (Renault Twitzy 45, shown in Table 1 ).
section_name: Suspension and Vehicle Models

section_content: As an example, in this work, three of the most used EHSA systems in the literature were selected: ball screw, rack and pinion, and cable transmission [14]. 

The EHSA study is based on the methodology proposed by L. Bowen et al. [32] for the analysis of these EHSA systems. Figure 3 shows the DC motor model used in the three selected EHSA systems.
section_name: Selection of EHSA Technology

section_content: As an example, in this work, three of the most used EHSA systems in the literature were selected: ball screw, rack and pinion, and cable transmission [14]. 

The EHSA study is based on the methodology proposed by L. Bowen et al. [32] for the analysis of these EHSA systems. Figure 3 shows the DC motor model used in the three selected EHSA systems. 

(mm) 1110 (mm) 1227
section_name: Selection of EHSA Technology

section_content: As an example, in this work, three of the most used EHSA systems in the literature were selected: ball screw, rack and pinion, and cable transmission [14]. 

The EHSA study is based on the methodology proposed by L. Bowen et al. [32] for the analysis of these EHSA systems. Figure 3 shows the DC motor model used in the three selected EHSA systems. From the model of the DC motor the relation between the current of the conductor i of the circuit, the constant of the torque k i and the moment is established by

On the other hand, applying Newton's second law

where τ g is the generator input torque; J g is the motor inertia; α g is the spindle rotation. 

From the lows Kirchhorff's voltage laws:

where U em f is the counter electromotive force that depends directly on the speed change and k e represents the counter electromotive voltage constant

Energies 2020, 13, 6110 8 of 25

By simplifying the above equations and applying Laplace's transform. 

Sorting and clearing up the torque in Laplace's domain:

Grouping in equivalent factors, C g and K g are expressed as:
section_name: Selection of EHSA Technology

section_content: The Ball-Screw Energy Harvesting Shock Absorber (BS-EHSA) system, which is defined as a linear actuator that transforms the linear movement into angular, is characterized by handling pitch angles from 1 to 10 degrees, friction coefficients from 0. 01 to 0. 003 that accounts for efficiencies above 90% [34]. 

The ball screw system is formed when the end of the screw is connected to the electric generator, while the nut is connected to the shock absorber. These systems are often used because they represent a simple solution related to the number of components connected and they are ease of assembly and robust. In Figure 4, the generator is coupled to the spindle, where τ g is the generator torque and α is the angular displacement. Some disadvantages of this technology is that some mechanical failures might occur, such as the possibility of reaching critical speed rotation, the wear of balls and tracks and reduction of system life due to fatigue in components depending on the maximum load and speed determined by the manufacturer. 

Considering the efficiency, the inertia, and the ratio of the dynamic equations of the ball screw in the transformation from linear to angular motion we have:

The relationship between the relative displacement z and angular displacement αs of Ball Screw and lead: Some disadvantages of this technology is that some mechanical failures might occur, such as the possibility of reaching critical speed rotation, the wear of balls and tracks and reduction of system life due to fatigue in components depending on the maximum load and speed determined by the manufacturer. 

Considering the efficiency, the inertia, and the ratio of the dynamic equations of the ball screw in the transformation from linear to angular motion we have:

α s = α g (13) The relationship between the relative displacement z and angular displacement α s of Ball Screw and lead:

where z represents the axial displacement of a nut on the spindle and l is the lead. Derivating and clearing τ g. τ g = J g. . 

where z represents the axial displacement of a nut on the spindle and l is the lead. 16 )

and where F a represents the input force of the path or excitation.
section_name: Ball Screw Model

section_content: The rack and pinion drive system (Figure 5 ) is very robust mechanically, although it has the disadvantage of gear backlash and therefore limited wear and tear and fatigue durability. By using the same generator for all EHSA systems, the formulation obtained in the Equations ( 15 )-(18) ball screw system can be used. 

Making a balance of power on the pinion shaft,

The balance on the gearbox shaft. 

= + (20)

From the relationship of linear and rotational movement between rack and pinion, By using the same generator for all EHSA systems, the formulation obtained in the Equations ( 15 )-(18) ball screw system can be used. 

Making a balance of power on the pinion shaft,

Energies 2020, 13, 6110 10 of 25

The balance on the gearbox shaft. 

From the relationship of linear and rotational movement between rack and pinion,

. 

In Equation ( 22 ), r p represents the pinion radius. The ratio of the angular displacement of the gearbox between the input and output. 

Substituting in Equation ( 19 ) and grouping the terms M eqRP, C eqRP, K eqRP.
section_name: Rack and Pinion System

section_content: The third technology for the comparison is the CD-EHSA (Cable Dynamics Energy Harvesting Shock Absorber [5], an innovative system that has the advantage of having a good ratio between rigidity and weight, good durability as shown in Figure 6. 

Energies 2020, 13, x FOR PEER REVIEW 11 of 27
section_name: Cable Driving Transmission System

section_content: The third technology for the comparison is the CD-EHSA (Cable Dynamics Energy Harvesting Shock Absorber [5], an innovative system that has the advantage of having a good ratio between rigidity and weight, good durability as shown in Figure 6. The equations for this model are shown below:

Based on the balance of forces condition in the cable,

Balancing energy on the pulley in the drive pulley, The equations for this model are shown below: Based on the balance of forces condition in the cable,

Balancing energy on the pulley in the drive pulley,

. 

. 

Making a balance of energy in the multiplier,

If we express the torque from the mass-spring-damper form,

The relationship between the input and output angular position of the gearbox. 

Clearing for F a, and grouping the equivalent coefficients become

(40)
section_name: Cable Driving Transmission System

section_content: Once the EHSA technologies to be compared were chosen, the driving cycle for which the comparison will be made is selected. In this case, the Worldwide Harmonised Light Vehicles Test Procedure (WLTP) cycle was chosen, as it has become mandatory for car manufacturers since 2017, which is supposed to be closer to reality [35]. 

The cycle is divided into different classes. In this study, we will focus on class 3, which is further subdivided into two subclasses according to its maximum speed. Class 3a with a maximum speed of 120 km/h and class 3b with a maximum speed greater than 120 km/h. The cycle is divided into 4 driving categories, low, medium, high, and super high as shown in Figure 7. 

* = (( + ) + )
section_name: Driving Cycle

section_content: Once the EHSA technologies to be compared were chosen, the driving cycle for which the comparison will be made is selected. In this case, the Worldwide Harmonised Light Vehicles Test Procedure (WLTP) cycle was chosen, as it has become mandatory for car manufacturers since 2017, which is supposed to be closer to reality [35]. 

The cycle is divided into different classes. In this study, we will focus on class 3, which is further subdivided into two subclasses according to its maximum speed. Class 3a with a maximum speed of 120 km/h and class 3b with a maximum speed greater than 120 km/h. The cycle is divided into 4 driving categories, low, medium, high, and super high as shown in Figure 7.
section_name: Driving Cycle

section_content: There are currently two main methods for measuring road roughness. The first one obtains the power distribution over different frequencies at which the signal is formed. It can be used to create
section_name: Road Roughness

section_content: There are currently two main methods for measuring road roughness. The first one obtains the power distribution over different frequencies at which the signal is formed. It can be used to create an artificial road profile and to characterize the test road segment. The second one is the International Roughness Index (IRI). 

To simulate the vibration generated by the road roughness, the road profile is obtained through the methodology established by ISO 8608 2016 [36]. According to this methodology, road roughness can be classified from a very good condition (A) to a very bad condition (H). In our case, a medium type C classification road was selected. A random road profile is generated through a unit input signal of white noise filtered using a low pass filter as follows:

where G q = 2. 56 × 10 -6 m/cycle is a mean value according to the class C road roughness; u corresponds to the vehicle speed; n 0 is the reference spatial equal to the value of 0. 1 m -1 ; and f 0 is the minimum cut-off frequency with a value of 0. 0628 Hz. Figure 8 shows the class of road type C obtained applying Equation (44). Equation ( 44 ) is a transfer function that can be used to find the road roughness in the time domain (Equation ( 45 )) [37, 38] :

. 

where z(t) is the road roughness and w(t) is a white Gaussian noise. The obtained road profile is shown in Figure 9 for a variable speed according to the WLTP cycle.
section_name: Road Roughness

section_content: where = 2. 56 10 m/cycle is a mean value according to the class C road roughness; u corresponds to the vehicle speed; is the reference spatial equal to the value of 0. 1 m -1 ; and is the minimum cut-off frequency with a value of 0. 0628 Hz. 

Figure 8 shows the class of road type C obtained applying Equation (44). Equation ( 44 ) is a transfer function that can be used to find the road roughness in the time domain (Equation ( 45 )) [37, 38] :

where () is the road roughness and w(t) is a white Gaussian noise. The obtained road profile is shown in Figure 9 for a variable speed according to the WLTP cycle. This profile is applied to the 4 wheels of the vehicle model, with a delay between the front wheels of 0. 1 s and 0. 2 s between the rear wheels. In this way, pitch and roll modes of the vehicle model are excited.
section_name: 𝐺(𝑠) = 𝑠 + 2𝜋𝑓

section_content: Once the vehicle model is set up (9 DOF), the EHSA technologies are selected (BS-EHSA, RP-EHSA and CD-EHSA), and road profile (type C) and testing cycles(WLTP driving cycle) are chosen, it is time to apply the proposed methodology (Figure 10 ) to select the best design option. Regarding the design of the EHSA, different objectives can be selected, and the methodology can be adapted to them. In this case, the design is aimed at selecting the best EHSA technology to replace a real light vehicle suspension [32], where a fixed damping factor of 1500 Ns/m was set up. This objective is at the same time realistic and useful if one intends to introduce EHSA technology in real vehicles. Once the design can cope with the damping objective, it can be further optimised for maximum energy recovery. 

Based on this design objective, the three EHSA models are studied and geometrically adapted to achieve the damping objective. 

The methodology will be applied according to the following flow chart This profile is applied to the 4 wheels of the vehicle model, with a delay between the front wheels of 0. 1 s and 0. 2 s between the rear wheels. In this way, pitch and roll modes of the vehicle model are excited.
section_name: Applied Methodology

section_content: Once the vehicle model is set up (9 DOF), the EHSA technologies are selected (BS-EHSA, RP-EHSA and CD-EHSA), and road profile (type C) and testing cycles(WLTP driving cycle) are chosen, it is time to apply the proposed methodology (Figure 10 ) to select the best design option. Regarding the design of the EHSA, different objectives can be selected, and the methodology can be adapted to them. In this case, the design is aimed at selecting the best EHSA technology to replace a real light vehicle suspension [32], where a fixed damping factor of 1500 Ns/m was set up. This objective is at the same time realistic and useful if one intends to introduce EHSA technology in real vehicles. Once the design can cope with the damping objective, it can be further optimised for maximum energy recovery. 

the same time realistic and useful if one intends to introduce EHSA technology in real vehicles. Once the design can cope with the damping objective, it can be further optimised for maximum energy recovery. 

Based on this design objective, the three EHSA models are studied and geometrically adapted to achieve the damping objective. 

The methodology will be applied according to the following flow chart Tables 2 and 3 show the characteristics of the generator and vehicle parameters to be used in the three EHSA models. According to our methodology, it is necessary to have the same testing conditions for our EHSA systems to be compared. We selected a motor whose parameters were used in the work reported by Lincoln [14], obtained from a commercial motor (Maxon Motor 218010) [39]. Based on this design objective, the three EHSA models are studied and geometrically adapted to achieve the damping objective. 

The methodology will be applied according to the following flow chart. Tables 2 and 3 show the characteristics of the generator and vehicle parameters to be used in the three EHSA models. According to our methodology, it is necessary to have the same testing conditions for our EHSA systems to be compared. We selected a motor whose parameters were used in the work reported by Lincoln [14], obtained from a commercial motor (Maxon Motor 218010) [39]. 

Table 2. The Generator parameters used in the three EHSA models.
section_name: Applied Methodology

section_content: Table 3. The commons vehicle parameters used in the three EHSA models.
section_name: Generator Parameters

section_name: Vehicle Parameters

section_content: Modelling conditions for each EHSA system are shown in Figure 1. 

Based on Equation ( 17 ), the parameters that intervene in the calculation of the damping are identified. To determine which of those parameters have greater influence within the system, a design of experiments (DOE) was carried out. The results of this study show the relative influence of the screw lead, the speed ratio of the gearbox and the external resistance (Figure 11 ):

Based on these results, we can modify some of these parameters to obtain the damping objective. The damping value of 1500 Ns/m can be obtained in different parameter combinations (Figure 12 ). In this case, resistance of 1 Ohm and 8 mm lead screw were selected due to the market availability. 

Modelling conditions for each EHSA system are shown in Figure 1. Based on Equation ( 17 ), the parameters that intervene in the calculation of the damping are identified. To determine which of those parameters have greater influence within the system, a design of experiments (DOE) was carried out. The results of this study show the relative influence of the screw lead, the speed ratio of the gearbox and the external resistance (Figure 11 ): Based on these results, we can modify some of these parameters to obtain the damping objective. The damping value of 1500 Ns/m can be obtained in different parameter combinations (Figure 12 ). In this case, resistance of 1 Ohm and 8 mm lead screw were selected due to the market availability. Table 4 shows the parameters used for the ball screw model simulation according to the results obtained by DOE to reach the damping target value of 1500 Ns/m
section_name: Case I. BS-EHSA

section_content: The same process was followed (Figure 1 ) with the rack and pinion system. A DOE study showed the most influencing parameters (Figure 13 ). Table 4 shows the parameters used for the ball screw model simulation according to the results obtained by DOE to reach the damping target value of 1500 Ns/m.
section_name: Case II. RP-EHSA

section_content: The same process was followed (Figure 1 ) with the rack and pinion system. A DOE study showed the most influencing parameters (Figure 13 ).
section_name: Case II. RP-EHSA

section_content: The same process was followed (Figure 1 ) with the rack and pinion system. A DOE study showed the most influencing parameters (Figure 13 ). In this case, the combinations of = 0. 01 m and = 6 Ohm were chosen. Table 5 shows the parameter values used for the rack and pinion in the simulation according to the results obtained by DOE to reach the damping target value. In this case, the combinations of r p = 0. 01 m and R e = 6 Ohm were chosen. 

Table 5 shows the parameter values used for the rack and pinion in the simulation according to the results obtained by DOE to reach the damping target value.
section_name: Case II. RP-EHSA

section_content: In this case, the methodology is applied to the CD-EHSA; following the same structure as it was shown in Figure 1. 

Based on Equation ( 27 ), the damping of the cable transmission system depends on several parameters. To analyse their influence, a DOE study is also carried out on them (Figure 15 ). Figure 16 shows the values of the effect of the external resistance and the radius of the drive pulley concerning the damping. Figure 16 shows the values of the effect of the external resistance and the radius of the drive pulley concerning the damping. 

In this case, due to construction availability, the selected combination of an external resistance of 22 ohms and a pulley radius of 0. 008 m was preferred. 

Table 6 shows the parameters of the cable transmission system used in the simulation according to the results obtained by DOE to reach the damping target value of 1500 Ns/m. Figure 16 shows the values of the effect of the external resistance and the radius of the drive pulley concerning the damping. In this case, due to construction availability, the selected combination of an external resistance of 22 ohms and a pulley radius of 0. 008 m was preferred. 

Table 6 shows the parameters of the cable transmission system used in the simulation according to the results obtained by DOE to reach the damping target value of 1500 Ns/m.
section_name: Case III. CD-EHSA

section_content: From the results obtained above, it was possible to limit the universe of combinations to achieve the damping objective value (1500). The three EHSA technologies were tested in a WLTP driving cycle, combined with a type of road C. The parameter selection in each ESHA model focused on those values that were functional from the practical point of view, seeking to have the least number of elements and to reduce the weight of the system. The selection of the parameters was based on DOE studies as in Figures 11, 13 and 15. In the first case, the BS-EHSA system in Figure 12 presented a wide range of possible combinations. However, to make the system functional, the parameters were reduced to integer values where a balance between efficiency and feasible ball screw lead was chosen. 

In the second case, the RP-EHSA system was studied (Figure 14 ). The system was very sensitive to changes in the radius of the pinion and the external resistance. However, to reach the damping target value, it was necessary to incorporate an external gearbox. 

In the third case, the CD-EHSA system showed many possible solutions (see Figure 16 ). The solution was narrowed by looking for the smallest pulley diameter to have the best speed ratio combined with a multiplier gearbox with a minimum speed ratio. 

Comparing the results of the road profile (Figure 9 ) vs. the displacement of the suspended mass (Figure 17 ) is evident the damping in the sprung mass. 

In the third case, the CD-EHSA system showed many possible solutions (see Figure 16 ). The solution was narrowed by looking for the smallest pulley diameter to have the best speed ratio combined with a multiplier gearbox with a minimum speed ratio. 

Comparing the results of the road profile (Figure 9 ) vs. the displacement of the suspended mass (Figure 17 ) is evident the damping in the sprung mass. Figure 19 shows the roll angle derived from the delays (Figure 20 ) at each wheel corner. The sequence of road input values presented to each shock absorber in the simulation is shown in Figure 20. Figure 19 shows the roll angle derived from the delays (Figure 20 ) at each wheel corner. The sequence of road input values presented to each shock absorber in the simulation is shown in Figure 20. From the results obtained for the BS-EHSA, a total average harvest of = 89. 58 Watts was obtained. 

The values obtained for the EHSA rack and pinion system are shown in Figure 22. The sequence of road input values presented to each shock absorber in the simulation is shown in Figure 20. 

In these conditions, energy was recovered for each of the EHSA systems from the WLTP drive cycle on each wheel. Figures 21 22 23 show the values of each shock absorber without considering the full-stop times in the driving cycle. We start by comparing the four wheels shown the energy recovered by the BS-EHSA Figure 21. In these conditions, energy was recovered for each of the EHSA systems from the WLTP drive cycle on each wheel. Figures 21 22 23 show the values of each shock absorber without considering the full-stop times in the driving cycle. We start by comparing the four wheels shown the energy recovered by the BS-EHSA Figure 21. From the results obtained for the BS-EHSA, a total average harvest of = 89. 58 Watts was obtained. 

The values obtained for the EHSA rack and pinion system are shown in Figure 22. From the results obtained for the RP-EHSA, a total average harvest of = 283. 38 W was obtained, with the highest value of energy recovered. This system uses a multiplier with a speed ratio of 1:10, which supports a higher value than the ball screw system. 

The values obtained for the last CD-EHSA system are shown in Figure 23. From the results obtained for the RP-EHSA, a total average harvest of = 283. 38 W was obtained, with the highest value of energy recovered. This system uses a multiplier with a speed ratio of 1:10, which supports a higher value than the ball screw system. 

The values obtained for the last CD-EHSA system are shown in Figure 23. From the results obtained for the CD-EHSA (Figure 23 ), a total average energy harvesting of = 132. 63 Watts was obtained, being the second-best. This system uses a multiplier with a speed ratio of 1:12 (Table 6 ), which supports a higher value than the BS-EHSA system but lower than the RP-EHSA systems. From the results obtained for the BS-EHSA, a total average harvest of P Tave = 89. 58 Watts was obtained. 

The values obtained for the EHSA rack and pinion system are shown in Figure 22. 

From the results obtained for the RP-EHSA, a total average harvest of P Tave = 283. 38 W was obtained, with the highest value of energy recovered. This system uses a multiplier with a speed ratio of 1:10, which supports a higher value than the ball screw system. 

The values obtained for the last CD-EHSA system are shown in Figure 23. 

From the results obtained for the CD-EHSA (Figure 23 ), a total average energy harvesting of P Tave = 132. 63 Watts was obtained, being the second-best. This system uses a multiplier with a speed ratio of 1:12 (Table 6 ), which supports a higher value than the BS-EHSA system but lower than the RP-EHSA systems.
section_name: Results and Discussion

section_content: After reviewing the relevant literature, it is clear that there are different technologies used in EHSA. It is also clear that each author designs, develops, simulates and validates these systems in a different way than those with similar systems or technologies. This paper presented a methodology that can compare the performance of different EHSA technologies against a common design objective. Specifically, a given damping coefficient was set as a design objective, corresponding to that of a real vehicle, with the challenge of achieving a functional EHSA shock absorber with each of the chosen EHSA technologies. It is shown that the methodology is general and applicable to any type of EHSA. In this case, it was applied to two of the most used EHSAs in the literature (Ball Screw and Rack and Pinion) and a little-known one based on cable dynamics, which validates its versatility. 

Therefore, three different EHSA technologies with different characteristics were compared to obtain a given damping factor. The same methodology was applied to the three EHSA technologies (with different design parameters), simulating a given drive cycle, a light-duty vehicle, and a specific road roughness profile producing normalized inputs. 

The selection of the EHSA parameters was based on those with the greatest effect on damping, which were simulated with the mathematical modelling of the EHSA system and with a design of experiments. 

Once the damping objective was achieved for each EHSA technology, the energy recovered was calculated in each system. Given these results, the highest average instantaneous energy recovered is the RP-EHSA system with the value of P Tave = 283. 38 Watts, followed by the CD-EHSA system with P Tave = 132. 63 Watts, and finally, the BS-ESHA system with P Tave = 89. 58 Watts. These results assure that each of the technologies achieves the objective damping factor. In this regard, the BS-ESHA was able to reach the damping factor with no gearbox. On the other hand, the RP-EHSA and the CD-EHSA systems require a gearbox to hit that damping value. The BS-EHSA system is the simplest design and does not add extra weight to the car that will, eventually, affect fuel consumption and represents higher costs. In this regard, the presented methodology has proven helpful to compare three different EHSA technologies for a given design objective with a fixed damping ratio. 

Therefore, using this methodology, it is possible to compare different EHSA technologies under the same testing conditions regarding vehicle model, road profile, driving cycle and so on. Design objectives can vary from one case to another, but they can be reliably compared. 

In future studies, different objective parameters should be considered, balancing the energy recovered and the weight and cost of the total system.
section_name: Conclusions

section_content: Funding: The authors thank the Global Nebrija-Santander Chair of Energy Recovery in Surface Transport for their financial support. This work was also supported by the Comunidad de Madrid [grant SEGVAUTO 4. 0-CM-P2018EEMT-4362 ] and the Agencia Estatal de Investigación [grant RETOS 2018-RTI2018-095923-B-C22 ].

section_name: Conflicts of