A new method for optimizing the parameters of torsional vibration dampers

Torsional vibration dampers (TVDs) are essential components for reducing the torsional vibration of a vehicle power transmission system (VPTS). This paper presents a new parameter optimization method for designing TVDs. The method combines the modal inertia and energy methods by adjusting the modal inertia method using an adjustment factor that is optimized using the energy method. The optimization of TVD parameters seeks to minimize the maximum torsional elastic potential energy of the rear axle near the resonance speed, so that the design variables can be optimized by a manual search process. The proposed method is employed to optimize the parameters of single-stage, two-stage parallel, and two-stage series TVDs coupled to a model VPTS. The damping effects of TVDs optimized by the modal inertia method, the energy method, and the proposed method were compared and analyzed, and the calculation efficiencies of the methods were evaluated. Results show that the proposed method provides better damping effects than the modal inertia method, and also provides far better computing efficiency than the energy method.


Introduction
Torsional vibration is a significant type of vibration that generally occurs in vehicle power transmission systems (VPTSs).Severe torsional vibration can greatly affect not only the driving comfort of a vehicle, but also the durability of vehicle components.In addition, the noise, vibration, and harshness (NVH) performance requirements of vehicles are subject to increasingly severe constraints, and commonly employed clutch torsional vibration dampers (TVDs), which are relatively inexpensive, may not conform to modern damping requirements.Moreover, dual mass flywheel type TVDs, which provide better NVH performance, are too expensive for use in economical passenger vehicles.As a result, equipping VPTSs with economical and suitably high-performance TVDs has become increasingly urgent.
From a design perspective, TVDs provide a better damping effect with increasing moment of inertia (abbreviated henceforth as inertia) under optimally tuned conditions [1], and under constant inertia conditions, the damping effect can be improved by dividing a single-stage TVD into a multi-stage TVD comprising some number of stages [2].However, providing an ever-increasing inertia requires increasing the mass of the components, which runs counter to the requirements for designing lightweight vehicles.In addition, increasing the number of TVD stages generally increases the costs, and each TVD stage requires additional space, which introduces additional design and installation challenges.From this perspective, improving current design optimization methods provides a necessary avenue for designing relatively low cost TVDs that meet both design specifications and damping performance requirements.
Current TVD designs have widely employed rubber components due to their low cost and convenience.At present, comprehensive investigations have been conducted to optimize and match the kinetic parameters of rubber TVDs [2][3][4][5][6][7].However, these investigations have mainly focused on suppressing the torsional vibrations of crankshafts.The operational principles of TVDs are similar to those of dynamic vibration absorbers (DVAs) and tuned mass dampers (TMDs), where a matched TVD facilitates anti-resonance in the VPTS, and greatly reduces vibration at the nature frequency.Nonetheless, intense vibration may occur in the TVD.
Research regarding the vibration suppression of undamped and damped systems with a single degree of freedom (SDOF) have been well-established, and evaluation criteria have primarily included the suppression of their vibration amplitudes ( ) [8], nonlinear responses [9], energy over their entire frequency range ( ) [10], transient responses [11] and their accelerations [12], and the suppressed frequency bandwidth and stability of the main system [13].However, studies regarding the vibration suppression of systems with multiple degrees of freedom (MDOF) started relatively late.Lagrange's interpolation algorithm has been utilized to calculate the system transfer function [14], and optimization was conducted using the damped least square feasible direction method.However, only dynamic systems with no more than 8 degrees of freedom (DOF) could be effectively studied due to the poor computing power at that time and the limited calculation accuracy of the Lagrange interpolation algorithm.More recently, single or multiple vibration modes of a main system were suppressed by rationally optimizing the parameters of one or more DVAs according to a minimum power flow criterion [15].Ozer and Royston [16] employed invariant point theory developed by Den Hartog for undamped MDOF systems to establish equations for the optimal parameters of DVAs through Sherman-Morrison matrix inversion theory.Thus, a semi-analytical optimal design method was developed.Ozer and Royston [17] again utilized matrix inversion theory to further investigate suppressing the vibrations of a damped MDOF system by optimizing the parameters of DVAs based on the weighted and minimized responses of SDOF or more DOF.Vakakis and Paipetis [18] utilized the polynomial method to minimize the transmissibility between two DOF to optimize the parameters of DVAs in an undamped MDOF system.Kitis et al. [19] proposed a method for optimizing the parameters of dual DVAs by minimizing the transmissibility between two DOF for a damped MDOF system.Cunniff [20] employed the Newton iteration method to search the maximum vibration responses of a main system, and the optimal parameters of DVAs were obtained by optimization using the quasi-Newton Method.Marano et al. [21] achieved multi-objective optimization by minimizing the maximum standard deviation of the acceleration of a SDOF and the failure probability of TMDs in damped MDOF systems.Lavan and Daniel [22] proposed a fast-converging and efficient method to minimize the total mass of multiple TMDs, and the allowable acceleration of each DOF was employed as a constraint for the suppression of multiple vibration modes in three dimensions for irregularly shaped buildings (an MDOF system).
Regarding the parameter optimization of TVDs for systems with MDOF, we note that, while minimizing the vibration amplitude peak of a SDOF can certainly reduce the vibrational energy of a main MDOF system, the optimum damping effect of TVDs on the other DOF cannot be assured.In addition, minimizing the sum of the vibration amplitude peaks of all DOF will certainly suppress the vibration of every SDOF; however, phase differences between each SDOF will prevent each SDOF from simultaneously attaining a maximum response.Therefore, this method cannot guarantee the full absorption of the vibrational energy of the entire system over any single period of oscillation.
Currently, the modal inertia and energy methods are most widely employed for the parameter optimization of TVDs when the vibration control of the overall VPTS or several DOF is taken as the evaluation criterion.According to Den Hartog's optimally tuned condition, an MDOF system in the modal inertia method is equivalent to a modal system based on the mode shape employed to match its TVD, and an analytical solution can be obtained efficiently without calculating vibration responses [2,23].However, the influence of damping effects between DOF is not taken into account, which is not appropriate for a damped MDOF system such as a VPTS.Therefore, the TVD parameters obtained by the modal inertia method does not represent an optimal solution when damping effects are taken into account.In contrast, the energy method considers the influence of damping by means of calculating responses.However, while this provides an optimal solution when damping effects are taken into account, it does so with a much lower computational efficiency.As a result, obtaining an optimal solution requires considerable time, particularly a multi-stage TVD resonating on a main system with a large number of DOF.
In this paper, we combine the advantages of the modal inertia and energy methods and mitigate their disadvantages by adjusting the modal inertia method using an adjustment factor that is optimized using the energy method.Our work is based on a modeled VPTS of a front-engine, rear-drive vehicle.Through the proposed approach, a reasonable matching between a TVD and the VPTS are realized, and the torsional vibration amplitude and maximum torsional elastic potential energy (TEPE) of the VPTS are reduced.

MDOF torsional vibration model of a VPTS
Due to the complicated structure of a VPTS, the lumped mass method was utilized to simplify the torsional vibration model, which means that a forced vibration model was established by taking the engine output torque as the excitation source.In accordance with the principle of simplification and equivalent calculation method [3], the inertia of each SDOF, the torsional stiffness, and damping coefficient beginning with the intermediate shaft of the gearbox were rendered equal to the engine crankshaft, depending on the transmission ratio of the gear wheels.Accordingly, a 39-DOF torsional vibration model was established, as shown in Fig. 1.The definitions of the inertia terms for the different components of the VPTS represented by each SDOF are listed in Table 1.

Fig. 1. VPTS torsional vibration model
The torsional vibration equation of the VPTS is: where θ, , and represent angular displacement, angular velocity and angular acceleration 39×1 column vectors, respectively.represents the column vector of the excitation torque, i.e., = [0,0,0,0, , 0, , 0, , 0, , 0, ⋯ ,0] × , and , , and represent the 39×39 inertia matrix, the torsional damping matrix, and the torsional stiffness matrix, respectively, which are expressed as follows: The output torque of each cylinder can be presented according to a Fourier series expansion to determine each harmonic excitation torque.Considering that the firing order is 1-3-4-2, the output torque values of the four cylinders can be derived as follows: Here, represents the average output torque of a single cylinder, represents the amplitude of the th order torque, represents the excitation frequency, is time, and represents the phase angle of the th order torque of the first cylinder ( = 0.5, 1, 1.5, 2, …).
The harmonic superposition method was employed to solve for the vibration amplitude of each SDOF.As shown in Fig. 2, the vibration amplitudes at the input ends of the gearbox, driveshaft, and rear axle all peak at an engine speed of about 1500 rpm, which represents a resonance phenomenon in the VPTS.This phenomenon can be reduced by a TVD.Prior to optimization, a natural characteristics analysis was firstly conducted for the torsional vibration model.Here, an undamped free vibration model is generally utilized to simplify the calculation: The characteristic equation was solved to obtain the third order natural frequency (50.58 Hz) corresponding to the resonance rotation speed at about 1500 rpm.The corresponding mode shape is shown in Fig. 3, which indicates that the torsional vibration of the complete engine crankshaft (DOF -) is inconspicuous, the clutch (DOF and ) undergoes large deformation due to its low torsional stiffness, which is much less than that of the other shaft segments, the vibration enlarges continuously as the transmission power passes through the input shaft (DOF -), intermediate shaft (DOF and ), and output shaft(DOF -) of the gearbox, the complete driveshaft (DOF -) and main reducer (DOF and ) vibrate increasingly violently, and the vibration decreases rapidly as the transmission power passes through the differential mechanism (SDOF ), half-axles (DOF and ), and reaches the tires (SDOF ).In an MDOF system, the SDOF having the largest amplitude is equipped with a TVD to absorb the greatest amount of vibrational energy.According to the vibration characteristics of the model VPTS and actual situations, a TVD can be optionally equipped at the back end (SDOF ) of the driveshaft.The mechanical model is shown in Fig. 4, where , , and represent the equivalent inertia, torsional stiffness, and torsional damping coefficient of the crankshaft for the TVD, respectively.Then, the torsional vibration equation is established as follows:

Modal inertia method
The modal inertia method regards the MDOF system as a SDOF system in the vibration mode according to the equivalent vibrational energy [23].By assuming a harmonic solution = Θ sin( ), where is the vibration response and Θ is the vibration amplitude of the th DOF.Hence, the velocity of the th DOF is = Θ cos( ), and the kinetic energy of the MDOF system is = ∑ .Assuming an installation of the TVD on the th DOF, and the motion of the th DOF can be replaced by SDOF , such that the kinetic energy is = , assuming that = , the equivalent inertia (modal inertia) is obtained as follows: The relative vibration amplitude of each SDOF in the vibration mode for an undamped system can be represented by the mode shape ⋯ ⋯ , which yields: A mechanical model of subjected to unit harmonic excitation coupled to a TVD is shown in Fig. 5, where the modal stiffness ( ) can be calculated by the expression = ⁄ , and represents the natural frequency of the main system.The optimal torsional stiffness ( ) and damping coefficient ( ) of a single-stage TVD can be solved from analytic solutions obtained in accordance with Den Hartog's optimally tuned formula [2,3]: where: = ⁄ represents the natural frequency of the TVD, = ⁄ represents the ratio of the inertia of the TVD to ; and = 2 represents the critical damping coefficient.

Energy method
In this method, the torsional stiffness and torsional damping coefficient of the TVD are taken as design variables, and as the constraint condition, the torsional vibration torque of each shaft section in the torsional vibration system should be less than the corresponding allowable torque.Because the torsional vibration of the rear axle dominates the overall NVH performance, the objective function seeks to minimize the maximum TEPE of the rear axle (beginning with component = 33), which is equipped with a TVD.The optimization model is as follows: Here, is the maximum TEPE of the rear axle when = , is the maximum TEPE of the rear axle near the resonance speed.we also note that the corresponding frequency near the resonance speed was divided into -1 equal parts, i.e., , , … , , which represents an even number of points; and is the time period of the torsional vibration system when = .The crankshaft completes two revolutions during a single excitation period, and the 0.5th order excitation has the largest time period in multiple frequency excitation; therefore = 4 /(0.5 ).Finally, ( , ) represents the inequality constraints.The responses of various-order excitation torque are superimposed to obtain the response of each SDOF based on the linear superposition principle: where represents the column vector of the vibration amplitudes subjected to the th order excitation torque of the torsional vibration system.The corresponding steady-state solution of vibration responses is given by the imaginary part of Eq. ( 15).

The proposed method
To consider the damping effect between various shaft sections, the MDOF system can be equivalent to a SDOF system, as shown in Fig. 6, in which is the modal damping coefficient.The value of can be obtained in accordance with the equal dissipated energy for the MDOF system and the equivalent SDOF system in the vibration mode.The dissipated energy of the MDOF system in the vibration mode can be obtained as follows: The dissipated energy of the SDOF system is given by: By equating and , we obtain the expression for : Fig. 6.Equivalent damped SDOF system Fig. 7. TVD and an equivalent damped SDOF As for a damped MDOF system subject to a harmonic torque, we assume a harmonic solution = Θ sin( + ), where is the phase angle of the th DOF.Hence, the velocity of the th DOF is = Θ cos( + ), and the kinetic energy of the MDOF system is = ∑ .Assuming that the TVD is installed on the th DOF, and the motion of the th DOF can be replaced by a vibrational SDOF.On the basis of energy-equivalence principle described in section 3.1, the inertia of the SDOF can be obtained as follows: When = , discrepancies develop between the actual relative vibration amplitude relationships of various DOF and the mode shape, and the damping effect may lead to a decrease in the relative vibration.By assuming that the th DOF has the largest vibration amplitudes in undamped and damped MDOF systems, we obtain: The viscous damping effect induces phase differences between various DOF, and we note from Eq. ( 19) that varies with respect to during a period of oscillation (i.e., = 2 / ).Consequently, the TVD optimized by the modal inertia method is inaccurate from the perspective of equivalent energy.
In Fig. 7, we assume that there is a set of optimal parameters of , , and can provide the TVD with best damping effect for MDOF system.Assume that the relationship of and is given by: where is defined as adjustment factor.To ensure that and damping ratio ( (2 ) ⁄ ) do not change, the and are also adjusted proportionately with respect to , thus: For an undamped MDOF system, the TVD with parameters optimized by the modal inertia method splits the original energy peak of the main system into two smaller and equal peaks.
However, for a damped MDOF system, the modal inertia method cannot provide two smaller energy peaks of equal height.In Fig. 7, may be considered as being adjusted properly to provide two smaller energy peaks of the same height.A TVD is employed to absorb excess system energy, and the modal inertia is adjusted to absorb more energy.From this perspective, the appropriate adjustment factor can be optimized by minimizing the maximum TEPE of the rear axle near the resonance speed.
For the optimization of TVD parameters in Fig. 7, a numerical optimization method is generally employed to minimize the maximal vibration amplitude of , thus to ensure that all peak values are of equivalent heights.Through calculation, the expression for the parameters of the TVD in terms of can be obtained: The substitution of Eq. ( 23) into Eq.( 14) yields the following TEPE optimization model for : . .( ) ≤ 0. (24)

Optimization of single-stage TVD parameters
Based on the VPTS torsional vibration model presented in Section 2, the value of calculated by Eq. ( 12) is 0.0058 Kg•m², and the value of calculated by Eq. ( 18) is 0.1159 N•m•s/rad.By assuming that the value of is 0.3.The TVD parameters were calculated using Eq. ( 23) and the maximum TEPE near the resonance speed of the rear axle was calculated by Eq. ( 24), and the results are listed in Table 2 for different values of .The maximum TEPE decreases continuously while increases from 1 to 1.9, and increases when = 2.0.Thus, we preliminarily conclude that the optimal is close to 1.9.The maximum TEPE was calculated for values 1.89 and 1.91, and the maximum TEPE was found to obtain a minimum when = 1.9.
The energy method given by Eq. ( 14) provides the following optimized TVD parameters: = 121.4572N•m/rad, = 0.1840 N•m•s/rad, and the maximum TEPE is 0.0235 J.The optimized parameters obtained using the modal inertia method, the energy method, and the proposed method were substituted into Eq.( 10), and the calculated vibration amplitudes at the input end of the rear axle (main reducer) are shown in Fig. 8.The figure indicates that the amplitude of the damped system employing TVD parameters optimized by the proposed method is further reduced at approximately 1500 rpm and increases slightly between 1000-1300 rpm compared with the amplitude when employing the TVD parameters optimized by the modal inertia method.On the other hand, the amplitude increases slightly between 1200-1500 rpm and decreases slightly between 1500-1700 rpm and is essentially equivalent at all other speeds compared with the amplitude when employing the TVD parameters optimized by the energy method.The proposed method obtains a better optimal solution than the modal inertia method, and within a much shorter time using a manual search of the adjustment factor than the energy method.From comparison tests, the proposed method is found to require less than 2 min of the computing time, but the energy method consumes 30 min.As such, the computing efficiency increased greatly while simultaneously ensuring a near optimal damping effect.

Optimization of parameters for two-stage parallel and series TVDs
In addition to single-stage TVDs, the proposed method is also applicable for the optimization of parameters for multi-stage TVDs.For enhancing the damping effect, a single-stage TVD may be converted into a two-stage parallel or series TVD, whose mechanical models are shown in Fig. 9.In contrast, the mechanical model of the matched MDOF system is shown in Fig. 10, whose torsional vibration equation is: In Fig. 9, a two-stage TVD will induce two new resonance frequencies near the original resonance frequency of .The presence of three close peaks tends to trap the optimization results in local attractive basins, which results in a less than optimal damping effect.A comparisons of the damping effects of the two-stage TVDs with parameters optimized using the local optimization solver and the global optimization solver are shown in Fig. 11 when = 1, where the excited frequency ratio represents ⁄ , and the vibration amplitude ratio represents the ratio of the vibration amplitude of to its static deflection (1 ⁄ ).Fig. 11 shows that only two peaks are obtained due to the local optimization of the two-stage parallel TVD, which represents a much lower damping effect than when using the global optimization solver, and the damping effect of the two-stage series TVD with parameters optimized using the local optimization solver is slightly weaker than that when employing the global optimization solver.It must be noted that, the results of two-stage TVD parameter local optimization are strongly correlated to the initial values selected, and thus a global optimization solver must be employed.
The optimization strategy is as follows: the total inertia of the two-stage TVDs remains unchanged and equal to the inertia of a single-stage TVD for optimization each time the value of is changed.The optimal parameters of the two-stage parallel and series TVDs and the maximum TEPE obtained based on Eq. ( 25) for different values of are listed in Tables 3 and 4, respectively.It must be noted that, the torsional damping coefficient of rubber is generally between 0.05-0.4N•m•s/rad, which should be the range of damping coefficient for optimization.
For the two-stage parallel TVD, the maximum TEPE decreases continuously with increasing from 1 to 2.1, and then increases for = 2.2.Thus, we preliminarily conclude that the optimal value of is close to 2.1.The maximum TEPE was calculated for ranging from 2.09 to 2.06, and the maximum TEPE obtained a minimum value when = 2.07.
For the two-stage series TVD, the maximum TEPE decreases continuously with increasing from 1 to 1.9, and then increases for = 2.0.Thus, we preliminarily conclude that the optimal value of is close to 1.9.The maximum TEPE was calculated for ranging from 1.89 to 1.87, and the maximum TEPE obtained a minimum value when = 1.88.
The damping effects of the TVD with parameters optimized by the modal inertia method and the proposed method were compared, and the vibration amplitudes at the input end of the rear axle are shown in Fig. 12.The figure indicates that the amplitude near 1500 rpm has decreased considerably, but a small increase is observed in the range 1000-1250 rpm compared with the damping effect of the TVD with parameters optimized by the modal inertia method.The improvement in the damping effect of the two-stage series TVD is particularly obvious for a very small increase at engine speeds below the resonance speed.
The following optimization model is employed for optimizing the parameters of a two-stage TVD using the energy method: (26) The problem of local convergence also exists when using the energy method.The tendency for the optimization results to become trapped in local attractive basins can observed in Fig. 13, where the damping effects of the two-stage TVDs with parameters optimized using the energy method with local optimization may be weaker than the damping effects obtained when the parameters are optimized with global optimization.Thus, a global optimization solver must again be employed, and the optimization results are listed in Table 5.The damping effects of the two-stage parallel and series TVDs optimized by the modal inertia method, the energy method with global optimization, and the proposed method were compared, and the vibration amplitudes at the input end of the rear axle are shown in Fig. 14.We can observe that the energy method with global optimization provided a slightly better damping effect than the proposed method in the 1300-1500 rpm range, although the proposed method generally produced a better damping effect for engine speeds greater than 1500 rpm.
The energy method with global optimization has very low computing efficiency, particularly for the large number of DOF of the main system, such that the method requires about 4 h computation time when the 2nd and 4th order torques which are only considered as the main causes for torsional vibration.The multiple design variables must be optimized by the energy method using a complex model like that shown in Fig. 1, for which the vibration responses must be solved.In contrast, the proposed method translates the multiple design variables into a single optimal design variable obtained by global optimization on the basis of the simplified model in Fig. 9, which provides low complexity and high computing efficiency, and requires about only 10 min.In addition, it is far more efficient to apply global optimization to the model in Fig. 9 than in the VPTS model shown in Fig. 1.

Practical applicability of the proposed damping technique
The proposed method can be improved at the condition that the natural frequency of the vibration system closes to the boundaries of excitation frequency domain.Fig. 15 shows a damped SDOF model consists of , , and , and subjected to unit harmonic excitation coupled to a TVD.If closes to the lower bound of , the height of the right peak can be adjusted lower than that of the left peak, the reason for this is that the probability of reaches the left peak frequency is largely less than that of the right peak.Through comparative calculation, the value of could be selected as 1.8.Then the graph of vibration amplitude ratio of against ⁄ is shown in Fig. 16(a).In contrast, if closes to the upper bound of , the height of the left peak can be adjusted lower than that of the right peak, and could be selected as 0.65.Then the vibration amplitude graph is shown in Fig. 16(b).The modal inertia method and the energy method cannot adjust the heights of the two peaks in these special conditions, and the proposed method after improvement can provide different peak heights as required by adjusting β.
Dr. Xiaodong Tan handled the work of writing the full manuscript.Dr. Lin Hua worked out a technology scheme of using torsional vibration damper to solve the problem of torsional vibration of the vehicle, and revised the manuscript at last.Dr. Chihua Lu put forward the main innovation point, and implemented parameters optimization by MATLAB programming.Dr. Can Yang established the torsional vibration model.Dr. Yongliang Wang plotted the plane graphics, and organized the optimization data, and created the tables.Dr. Sheng Wang translated the Chinese manuscript into English.

Conclusions
By considering the inherent torsional vibration characteristics of the entire VPTS as the starting point, we combined the modal inertia method with the energy method to propose a new method for optimizing the parameters of a TVD.These three methods were applied for optimizing the parameters of single-stage, two-stage parallel, and two-stage series TVDs for the model VPTS shown in Fig. 1.The proposed method provides peak values of equivalent heights (excluding the special conditions described in the research extension) in the energy curve by adjusting , and the local convergence is not existed when the TVD parameters listed in Tables 3 and 4 are in the condition of global convergence solutions.In addition, the proposed method has been accepted by the cooperative enterprise, and the production and performance test of TVDs are now being considered.Due to limited space, these contents will be discussed in the following research.
The main conclusions are as follows: 1) The vibration amplitude at the input end of the rear axle for TVD parameters optimized using the proposed method may increase slightly for engine speeds less than the original resonance speed relative to that obtained for TVD parameters optimized using the modal inertia method.Moreover, the vibration amplitude growth rate may rise gradually along with increasing .
2) The method proposed in this paper combines the advantages of high computing efficiency provided by the modal inertia method and the excellent damping effect of the TVD using parameters optimized by the energy method with global optimization.On the other hand, the disadvantages of insufficient damping effect provided by the TVD using parameters optimized by the modal inertia method and the poor computing efficiency of the energy method with global optimization are both avoided.These merits make the proposed method very practical and effective in engineering applications, which can be extended to the optimal design of DVAs and TMDs coupled to an MDOF system.
3) The merits of the proposed method are increasingly apparent as the number of DOF of the system increases, particularly when employed for the parameter optimization of multi-stage TVDs coupled to the system.However, the advantages are not particularly obvious when a single-stage TVD is coupled to a system with a smaller number of DOF.

Fig. 4 . 5 .
Fig. 4. Mechanical matching model of the TVD Fig. 5. TVD and modal inertia 3. Description of modal inertia and energy methods

Fig. 8 .
Fig. 8. Vibration amplitudes at the input end of the rear axle for the undamped VPTS and damped VPTS employing TVD parameters optimized by the various methods

9 . 10 .Fig. 11 .
Fig. 11.Comparison of damping effects for the two-stage TVDs optimized by the local and global optimization solvers NEW METHOD FOR OPTIMIZING THE PARAMETERS OF TORSIONAL VIBRATION DAMPERS.XIAODONG TAN, LIN HUA, CHIHUA LU, CAN YANG, YONGLIANG WANG, SHENG WANG

Fig. 12 .Fig. 13 .
Fig. 12.Comparison of damping effects for the modal inertia and proposed methods a) Two-stage parallel TVD b) Two-stage series TVD Fig. 14.Comparison of damping effects by different methods

Table 1 .
Definitions of the different DOF Symbol Definition Inertia of the torsion damper at the free end of the engine crankshaft , , , and Inertia values of the stepped shafts , , , , and Inertia values of the main journal , , , and Inertia values of the 4 cylinder crank-connecting rod mechanisms Inertia values of the flywheel and driving part of the clutch , , , and Inertia values of the driven part of the clutch and input shaft of the gearbox and Inertia values of the intermediate shaft of the gearbox , , , and Inertia values of the output shaft of the gearbox Inertia of the front universal joint and internal spline Inertia of the front half shaft Inertia of the intermediate spline Inertia of the intermediate universal joint , , Inertia values of the rear half shaft Inertia of the rear universal joint and Inertia values of the main reducer driving and driven bevel gear Inertia of the differential mechanism and Inertia values of the two half-axles Equivalent inertia of wheels and the entire vehicle ( = 1, 2, …, 38) Torsional stiffness between various DOF Torsional stiffness of tires ( = 1, 2, …, 38) Torsional damping coefficient between various DOF Torsional damping coefficient of tires , , , and External damping coefficients of the 4 cylinder piston crank-connecting rod mechanisms

Table 2 .
Optimized single-stage TVD parameters obtaining the smallest maximum TEPE value

Table 3 .
Optimized results of the two-stage parallel TVD shown in Fig.9(a)

Table 4 .
Optimized results of the two-stage series TVD shown in Fig.9(b)

Table 5 .
TVD parameter results optimized by the energy method with global optimization