A Novel E-DVA Module Synthesis Featuring of Synergy between Driving and Vibration Attenuation

To attenuate the negative effects brought by heavy unsprung mass of the decentralized driving electric vehicle, a novel e-DVA module featuring of synergy between driving and vibration attenuation is proposed in this paper. It presents the advantages of compact structure and low cost. Structure design proves the feasibility of the e-DVA module. Kinematic analysis of the slidercrank mechanism is carried out to conclude the transmission ratio ripple under road excitation. After parameter matching and optimization of the e-DVA module based on the H 2 /H ∞ norm criterions, vertical dynamics analyses in both frequency and time domains are conducted theoretically to prove the performance improvements on the ride comfort and handling stability under the constraint of DVA deflection bound.


Introduction
Regarding fossil fuel saving and environmental protection, some research institutions and corporations have proposed a variety of electric vehicles, the majority of which are based on the conventional internal combustion engine platform [1].Although it is still in the stage of concept promotion featuring of individual driving wheels, the decentralized driving electric vehicle has shown some advantages.Eliminating the components including differential in the simplified chassis, the space within vehicle body is preserved.In addition, as driving torque on each wheel can be achieved and controlled precisely, it provides new possibilities for the dynamics control strategies [2].However, fixing the motor in wheel or close to wheel increases the unsprung mass, which deteriorates the performances of ride comfort, handling stability, and safety.As a result, its promotion is hindered [3].
To reduce the unsprung mass, considerable investigations have been done by researchers around the world.Except lightweight designs on material and structure [4,5], positioning the motor on the vehicle body is a common idea to transfer it as a sprung mass [6], as depicted in Figure 1.However, with the introduction of universal joints, the transmission system efficiency is reduced, and the stiffness of the elongated shaft is comparatively small.
In recent years, a lot of studies on active suspension and semiactive suspension are carried out by researchers [7,8].In spite of the improvements on ride comfort and handling stability, high cost actuators and complicated controllers contribute to the failure to be widely promoted.Consequently, designs based on passive suspension are more practicable approaches for vehicles in production at present.
Our research group has been investigating techniques to attenuate the negative effects brought by heavy unsprung mass.Taking the gearbox as a trailing arm, Figure 2(a) illustrates an integral of a driving system and chassis in a top view [9].Taking the vehicle body as a reference, as the motor swings far away from the wheel under excitation, the vibration amplitude of the motor is much smaller than the wheel.Consequently, the negative influences of unsprung mass are reduced comparing with the architecture of hub motor direct driving, and it improves compactness as the shaft and universal joints are eliminated.Furthermore, with the introduction of the reducer between the motor and wheel, it fully utilizes the high power density of high speed motor and reduces the cost of motor, especially in the state of low speed and high torque.But the attenuation performance depends on the structure parameter.Study results show that the longer the distance between axles of wheel and motor, the lower the equivalent unsprung mass [10].And Figure 2(b) shows a prototype application in a decentralized driving electric vehicle.Dynamic vibration absorber (DVA) has been widely applied to attenuate the resonance vibration of primary system in mechanical systems [11].But in vehicle suspensions, both the introduction of an auxiliary mass and the space restriction near the wheel limit its application [12,13].
Taking the driving motor itself as a DVA, a novel electric wheel concept named e-DVA module throughout the paper based on slider-crank mechanism is proposed in Section 2, featuring of the synergy between electric driving and dynamic vibration absorbing.In Section 3, kinematics analyses of the slider-crank mechanism composed by two gear-pole composite mechanism are carried out, and the wheel speed ripple is concluded.In Section 4, under a random excitation, the parameters of the e-DVA module are matched and optimized based on the  2 norms of body acceleration, suspension deflection, and tire dynamic load under the  ∞ norm constraint of DVA motion bound.At last, comparing with the in-wheel motored electric wheel, the vertical dynamics analyses are carried out in both frequency and time domains to verify the improvements on ride comfort and handling stability.

E-DVA Module Architecture
With the advantages of integration and adaptation, modular design is very convenient for mass production and recycling, which shortens the period of investigation and manufacturing [14].Michelin and Siemens corporations, respectively, proposed the electric wheel modules integrating the hub motor, active suspension, and brake system in wheels [15,16].However, both wheel modules are of complicated structures and high cost for direct driving of hub motor and introducing active suspension.A novel e-DVA module architecture is proposed in this section, and a simplified schematic diagram based on slidercrank mechanism is demonstrated in Figure 3.It is mainly composed of the tire, rim, bearing, hub, braking disc, caliper, gearbox, and motor.The output shaft of the motor  is connected with gear  by spline, and the driving shaft is connected with gear  by spline.Gear  is engaged with gear , while gear  is engaged with gear , and gear  shares the same shaft with gear .The two ends of link  are connected with gear  and gear , respectively, while the two ends of link  are connected with gear  and gear .The upper end of the spring is fixed with the motor , and the lower end of the spring is fixed with the knuckle.The motor , the link , the link , the spring, and damper construct a slider-crank mechanism, in which the motor  acts as the slider.Taking the knuckle as a reference, under the excitation of uneven road, the motor  would move up and down along the common axle of the spring and damper, while both link  and link  conduct planer motions, guaranteeing the engagements of the gears in motion.
Detailed structure design is demonstrated in Figure 4. To avoid interference with other components, the motor is constrained moving vertically along the diameter of the rim to fully utilize the space in wheel.For compactness, a novel bearing is adopted, whose inner ring is regarded as the knuckle, while the outer ring is fixed with the rim.So that the components including transmission system can be arranged    in wheel, which is originally occupied by the conventional hub bearing and knuckle in sequence.A guide mechanism is introduced to restrict the motor motion.
In the e-DVA module designed for a low speed electric vehicle, a 15-inch tire is selected.According to the maximum speed requirement of 70 km/h, the power system parameters are designed and matched as shown in Table 1.
In conclusion, the outstanding feature of the e-DVA system proposed in this paper is the synergy between electric driving and dynamic vibration absorbing.Serving as a DVA, the motor attenuates the resonance vibration of vehicle body without any auxiliary mass.Performing as the links in the slider-crank mechanism, the gearbox decouples the rigid connection between the motor and the wheel and decreases the torque requirements for the motor.

Kinematics Analyses
In Figure 3, as gear  is engaged with gear  and link  performs the planer motion when the e-DVA module is under excitation, an equivalent planetary gear mechanism composed of gear , gear , and link  is analyzed in kinematics in this section.Supposing that the rotational speeds of gear , gear , and link  are   ,   , and   , respectively, and the teeth of gear  and gear  are   and   , then the rotational speed relationship shown below is satisfied: In a similar way, for the planetary gear mechanism composed of gear , gear , and link , supposing that the Gear  and gear  share the same shaft by splines, which means   =   .
After computation, the speed of gear  is obtained as follows: It is concluded from the equation above that, as the motor moves along the common axle of spring and damper with respect to knuckle, the transmission ratio is not stable under uneven road excitation.The speed ripple of the wheel can be depicted as follows, which is related to the rotational speeds of link  and link : The influence of wheel speed ripple would be analyzed in Section 4 combining with a specific example.

Vertical Dynamics Analyses
4.1.Equivalent Mass of Gear System.In this paper, the gear system is introduced to reduce the torque requirement of the motor, and the motor cost is reduced meanwhile.To establish the model precisely, the gear system should be taken into consideration.The influence of the gear system mass on the vertical dynamics is supposed to be analyzed by distributing the gear mass to the motor and wheel, respectively.According to the parallel axis theorem, the first stage reducer is transferred to the motor mass when the moments of inertial are kept equivalent: in which  * 3 is the equivalent DVA mass.In the same way, the second stage reducer is transferred to the wheel mass: in which  * 1 is the equivalent wheel mass.It is concluded that, with the introduction of the gear system, the equivalent DVA mass and wheel mass are increased a little, while the unsprung mass is decreased significantly without in-wheel motor.

Model Establishment.
According to the e-DVA module architecture, a quarter vehicle model of three-degree-offreedom (DOF) is established in Figure 5, where  * 1 is the equivalent wheel mass (the unsprung mass),  2 is the sprung mass, and  * 3 is the equivalent DVA mass;  1 is the tire stiffness,  2 is the suspension stiffness, and  3 is the DVA stiffness;  2 and  3 represent the suspension damping and DVA damping, respectively, and the tire damping is small enough to be neglected; ,  1 ,  2 , and  3 are the vertical displacements of the road, wheel, vehicle body, and DVA, respectively.
According to Newton's second law, the dynamic differential equation is established, as shown in where the mass matrix is as follows: damping matrix is as follows: stiffness matrix is as follows: The system generalized force is the road excitation K 1 , in which Performing Laplace transformation on (7), the transfer functions of the output displacements with road vertical displacement  are concluded, as depicted in As a result, to reflect the ride comfort, the transfer function of body acceleration z 2 with the input vertical speed q is shown in The transfer function of the wheel dynamic load   with the input vertical speed q is depicted as follows, which demonstrates the handling stability: The suspension deflection  2 −  1 and DVA deflection  3 −  1 should be in their motion bounds, and the transfer functions are illustrated similarly in Regarding the parameter optimization of a single freedom DVA with damper, extensive researches have been conducted [17].The researches aimed at determining dynamic vibration system parameters including the tuning frequency and damping ratio to attenuate the main system vibration.However, in the vehicle suspension design, besides the ride comfort, objectives including handling stability and motion interference should also be considered.Furthermore, regarding an actual system, the driving motor parameters are previously determined by the power requirement of the vehicle.This paper aims at improving the ride comfort and handling stability through optimizing the architecture of the power system without influencing the power performance.Restricted by the motor structure characteristic and manufacturing, the power performance is probably influenced if the motor mass is varied.As a result, although it is feasible to include the motor mass in the optimization, we take the motor mass as a hard constraint.In conclusion, the optimized parameters involve the DVA stiffness and damping under multiobjectives.
As is previously acknowledged, the performances of ride comfort and handling stability are reflected by the vehicle body acceleration and tire dynamic load, while suspension deflection should be in the permitted range; otherwise the system dynamics would become extremely nonlinear.As for the e-DVA module in this paper, in addition to the factors considered above, the deflection of the DVA should also be restricted to avoid interference with other components.
However, for a determined primary system and motor mass, adjusting the stiffness and damping of the DVA system has little influence on the peak values at the first-order resonance frequency.Alternately, the objective is to reduce the total vibration energy of the system over all frequencies when the primary system is subjected to random excitation.In this optimization criterion, the RMS values of the frequency response curves are optimized to be minimums.If the input road excitation is independent, the RMS values can be represented by the  2 norms in the physical sense.Consequently,  2 norms of the body acceleration, suspension deflection, and tire dynamic load are taken as the objectives of ride comfort and handling stability.
According to the variation ranges of the above three items, specific weighting coefficients are evaluated to formulate the criterion function of the optimization objective, as depicted in ( 16) in detail.One has In a similar way, the minimal of the DVA deflection amplitude at the resonance frequency, represented by  ∞ norm, is evaluated as the objective of motion interference. is set as the travel bound, as depicted in       ( 3 − 1 )∼ q    ∞ < .
In summary, considering the performance improvements and interference, the objective can be depicted as follows: (18) Parameter Matching of a Specific E-DVA Module.Supposing a vehicle drives on  class road with a speed of 20 m/s, the power spectrum density (PSD) of road is depicted in In the previous equation,  0 is the referenced space frequency;  0 = 0.1 m −1 ; ( 0 ) is the uneven coefficient;  is the vehicle speed.The frequency range of integration is from 0.1 Hz to 100 Hz.
According to ( 5)-( 6), considering the influence of gear system, the equivalent DVA mass and wheel mass are computed, and simulation parameters of a specific example are listed in Table 2.According to the  2 / ∞ norm criterions, for the e-DVA system with different stiffness and damping matches, Figures 6-8 show the RMS maps of body acceleration, suspension deflection, and tire dynamic load, as well as peak value map of DVA deflection in Figure 9.
The weighting coefficients are determined by the relative importance of the body acceleration, suspension deflection, tire dynamic load, and DVA deflection.However, in the optimization process, it is concluded that the suspension deflection is hardly influenced by the DVA stiffness and damping as its perturbation range is only 0.4 mm, and it locates in a suitable range, as depicted in Figure 7.As a result, the suspension deflection is neglected.In conclusion, as the ride comfort is preferred for passengers, the weighted performances of body acceleration, tire dynamic load, and DVA deflection are set as 6 : 5 : 4 in the variation ranges of both DVA stiffness and damping, while the weighting   coefficient of suspension deflection is set as 0, as shown in the following: From Figures 7-9, the variation ranges of suspension performances are concluded as follows: Consequently, the weighting matrices  and  are determined: The suboptimal solution of ( 16) can be concluded, in which the DVA stiffness  3 is 42000 N/m and damping  3 is 1200 N⋅s/m.From Figures 10-13, it is shown that, in the human sensitive frequency range of 4-8 Hz, the power spectral density of body acceleration in the e-DVA module decreases significantly comparing with the in-wheel motored system.In the whole frequency range, the second peak is shifted and damped, as well as in the PSD diagrams of suspension deflection and tire dynamic load.It can be explained as follows.As in-wheel motored system is a typical twodegree-of-freedom system, there are two peak values in the frequency domain response, as depicted from Figures 10-12, in which, the first peak value corresponds to the amplitudes of suspension performances at the body natural frequency, while the second peak value corresponds to the amplitudes at the wheel natural frequency.In this paper, transferring the motor mass as the DVA mass from the wheel, the ratio of sprung mass to the unsprung mass is increased, and the wheel natural frequency is reduced.As a result, the second peak value is shifted to the lower frequency.On the other hand, as the motor is transferred from the wheel, the negative effect of heavy unsprung mass is decreased.Further, the DVA owns the function of attenuating the vibration amplitude at the resonance frequency.Consequently, the amplitude is damped at the second resonance frequency.

Simulation in Time Domain
. Also under the excitation of  class road at the speed of 20 m/s, the body accelerations, suspension deflections, and tire dynamic loads are compared between the e-DVA module and the in-wheel motored system in time domain, as depicted from Figures 14-16.And Figure 17 shows DVA deflection response.
And the RMS values are computed according to the simulation above, as depicted in Table 3.
The suspension performances are evaluated in both frequency and time domains.The second peak value is very important, which presents significant reduction of the suspension performance amplitude in the range near resonance frequency.In time domain responses, although the attenuation effect is especially significant in the range near the second resonance frequency, there is no significant difference as the excitation is random including many frequencies.RMS values are utilized to compare the performances numerically, as indicated in Table 3.It is concluded that both the suspension performances of ride comfort and handling stability are improved, as the body acceleration, tire dynamic load, and suspension deflection are reduced 17.88%, 15.82%, and 2.70%, respectively.

Influence Analysis of Speed Ripple.
At last, as the transmission ratio is not stable as concluded in Section 2, the influence of speed ripple on wheel speed is analyzed combining with a specific e-DVA module with parameters above.Supposing that the rotational speed of the motor is 5000 r/min, then the wheel speed is 3538 ∘ /s without any road excitation.In-wheel motored system e-DVA module  When the motor performs vertical sinusoidal motion with the amplitude of 6.5 mm and angular frequency of 2 under vibration, the speed ripple is 12.34 ∘ /s as concluded in ADAMS software.Comparing with wheel speed of 3538 ∘ /s, the speed ripple is much smaller, which has little influence for the common driving and proves the feasibility of application.

Conclusions
(1) To suppress the negative effects brought by the heavy unsprung mass of the in-wheel motored system, by taking the motor as a dynamic vibration absorber, a novel e-DVA module is proposed in this paper, featuring of the synergy between driving and vibration attenuation.Serving as the crank, the reducer not only transmits the power but also decouples the rigid connection between the motor and the unsprung mass.The e-DVA module presents the advantages of compact structure and low cost.
(2) The stiffness and damping of the e-DVA system are matched and optimized based on weighted  2 / ∞ norms of suspension performances.After theory deduction and simulation analysis, comparing with the in-wheel motored system, significant improvements on ride comfort and handling stability are achieved through vertical dynamic analysis in both frequency and time domains.There is little influence of speed ripple on the wheel speed analyzed in ADAMS software, which proves the feasibility of application.

Figure 1 :
Figure 1: Schematic diagram of positioning the motor on the vehicle body.

Figure 2 :
Figure 2: Integral of electric drive and single trailing arm suspension.

1 Figure 5 :
Figure 5: A quarter vehicle model of three DOF.

1 )
ff n e s s (N •m − 1 ) D am pi ng (N •m •s −

Figure 9 :
Figure 9: Peak value map of DVA deflection.

Figures 10 -
And Figure13shows the power spectral density of DVA deflection.From, it is shown that, in the human sensitive frequency range of 4-8 Hz, the power spectral density of body acceleration in the e-DVA module decreases significantly comparing with the in-wheel motored system.In the whole frequency range, the second peak is shifted and damped, as well as in the PSD diagrams of suspension deflection and tire dynamic load.It can be explained as follows.As in-wheel motored system is a typical twodegree-of-freedom system, there are two peak values in the frequency domain response, as depicted from Figures10-12, in which, the first peak value corresponds to the amplitudes of suspension performances at the body natural frequency, while the second peak value corresponds to the amplitudes at the wheel natural frequency.In this paper, transferring the motor mass as the DVA mass from the wheel, the ratio of sprung mass to the unsprung mass is increased, and the wheel natural frequency is reduced.As a result, the second peak value is shifted to the lower frequency.On the other hand, as the motor is transferred from the wheel, the negative effect of heavy unsprung mass is decreased.Further, the DVA owns the function of attenuating the vibration amplitude at the resonance frequency.Consequently, the amplitude is damped at the second resonance frequency.

Figure 16 :
Figure 16: Tire dynamic load variation with time.

Table 1 :
Power system parameters.
rotational speeds of gear , gear , and link  are   ,   , and   , respectively, and the teeth of gear  and gear  are   and   , then the following equation illustrates the speed relationship: