A Study on the Dynamic Performance for Hydraulically Damped Rubber Bushings with Multiple Inertia Tracks and Orifices : Parameter Identification and Modeling

Hydraulically damped rubber bushings (HDBs) are important for vehicle noise, vibration, and harshness (NVH) performance as they are able to decay the vehicle’s oscillation induced by engine and road.The dynamic stiffness and loss angle of anHDB are crucial and it is significant to investigate the relations between the design parameters with the dynamic stiffness and loss angle. Therefore, the force-deflection relation of the HDB is measured statically and the dynamic stiffness and loss angle are measured dynamically and the test data are analyzed with a view to examine how the measurement results are influenced by the design parameters (the number of the fluid tracks). Compared with the results predicted by a nonlinear lumped parameter model whose parameters are extracted by a parameter identification technique, using the model, the effect of the main rubber and the fluid track on the dynamic stiffness and the loss angle is investigated. A unified analytical model of HDB is also developed with the purpose of predicting the static and dynamic characteristics, and the predictions are shown to be well correlated with the measurement data. The good correlation suggests the validity of the model and the parameter identification implementation.


Introduction
Motivation.Rubber bushings are commonly used in engine subsystems, vehicle body, and vehicle suspension to damp their oscillations excited by dynamic loads, thus improving driving safety, ride comfort, and handling performances [1][2][3][4][5][6][7][8][9].Compared with conventional rubber bushings, an hydraulically damped rubber bushing (HDB) can provide a high viscous damping coefficient in certain frequency range, while exhibiting amplitude dependent properties and frequency dependent properties [5][6][7][8][9], which is one of the advantages of an HDB.In designing the frequency characteristics of an HDB, designers take care of the value and peak frequency of the loss angle of the HDB mostly, which are dependent upon the performance of the main rubber spring, the fluid property, the shapes of two chambers, the fluid tracks, and so forth.Some effective methods and models have to be invited to calculate the influences of those factors on the dynamic characteristics of an HDB during the design stage.The dynamic performances of an HDB are often characterized with two terms, dynamic stiffness and loss angle, which are defined as the amplitude ratio of the load response to the displacement response and the phase difference between the load and displacement responses in frequency domain, respectively [4][5][6][7][8][9][10].
The lumped parameter (LP) model [10] is often employed to estimate the dynamic performances of an HDB [1,2,[4][5][6][7][8][9].Each lumped parameter represents a certain physical meaning, and the dynamic characteristics of the HDB can be described as a function of lumped parameters explicitly [9,10].However, the traditional method of parameter identification is a high-cost and time-consuming task for some lumped parameters, such as the HDB's effective piston area and the compliances of two fluid chambers.To obtain these design Shock and Vibration variables of an HDB, a new parameter identification method is proposed based on the characteristic frequency points of the dynamic performances of HDB in this paper.
During the design stage of an HDB, it is not easy to adjust the maximum frequency of loss angle to the same frequency of the most disruptive excitations with one fluid element due to geometry limitations of the fluid track [11].Besides, the hydraulic damping mechanism is relatively sensitive to the operating environment while multiple damped dynamic vibration absorbers have good robustness against it [12].However, the influences of fluid track number on the dynamic characteristics of an HDB with multiple inertia tracks and orifices as well as its practical application issues are not fully investigated.
Therefore, the motivations of this paper are to develop a superior parameter identification method for the HDB and to reveal the influences of the geometry size and the number of fluid tracks on the dynamic characteristics of the HDB by test and simulation approaches.
Literature Review.The first HDB was introduced in the mid-1980s on MacPherson-type front suspension systems [11].Similar designs are already used as HEM and are adapted for application with cylindrical sleeve-type bushing joints simply [11,13].There are also many patents [14][15][16][17][18][19][20] that claim the dynamic characteristics features of some specific HDBs in frequency domain; however, no analytical analysis or even measured dynamic properties are provided [4][5][6][7][8].Further, very limited articles have reported the dynamic properties characterization and modeling issues of HDB.
Lu and Ari-Gur [1,2] present a linear LP model for an HDB with two inertia tracks and derived the natural frequency analytically.Then, the model is applied to the simulation of an HDB with multiple inertia tracks.However, their model is not validated experimentally.Sauer and Guy [3] present a fluid system model for an HDB with a bypass track and one inertia track in parallel, which serves as a relief valve securing the HDB and the adjacent structure from undesirable shocks, but the detailed calculation conclusions are not carried out.Shangguan and Xu [9] propose a nonlinear lumped parameter state-space model (by introducing a nonlinear track resistance   ) for an HDB with one inertia track, and the model is verified by experiments, while the analytical study is conducted based on the linearized model.Sevensson and Hakansson [21] and Pan et al. [22] provide an empirical mathematic model including a nonlinear spring with several fluid elements and elastoplastic elements coupled in parallel, but the calculation results show that the accuracy of the model is not satisfactory for low amplitude excitation.Arzanpour and Golnaraghi [23,24] present a linear LP model of an HEM employing the same approach as that for an HDB with one inertia track.Chai et al. [4][5][6][7][8] initiate a laboratory device with alternative internal configurations, and the dynamic characteristics for different configurations are investigated comparatively, as well as the steady state and transient responses of the device.However, the influence of the number of the fluid tracks on the dynamic characteristics is not discussed.The parameters of LP models of HDB with multiple inertia tracks are identified by methods of curve fitting and Nyquist diagrams [4,8].However, the phenomenon of the excitation amplitude-invariant fixed points on the dynamic performance curves of HDBs with different internal configurations is not investigated.In summary, no uniform LP models have been developed to analyze the influences of the number of long fluid tracks and orifices on the dynamic characteristics of HDBs.A thorough understanding of the mechanism of HDBs is needed.

Scope and Objectives.
As a passive vibration isolation component in the vehicle suspension and steering wheel assemblies, HDBs play an important role in dealing with vehicle system vibration issues, such as brake judder and shimmy [3,11].The configuration of a typical HDB and the runner plate are illustrated in Figures 1 and 2, respectively.Its geometric and material parameters are listed in Table 1.
The elastomeric rubber spring, enclosed by inner and outer metal sleeves, has two functions: to support the static load of the vehicle and to provide partial static stiffness and somewhat damping of the suspension system.The typical passive hydromechanical bushing usually consists of two almost same fluid chambers connected by one or more inertia tracks and/or orifices.Typically, the fluid chambers filled with an antifreeze and water mixture which has appropriate dynamic and thermal properties required for isolation and control [25].In the event of a dynamically oscillatory excitation, the chamber pressures change in radical direction and the volume compensation will take place between two fluid chambers.Thus, the fluid flows back and forth through the fluid tracks as an oscillating mass to provide hydraulic damping [4-9, 11, 13].
Based on the above analysis, HDB's engineering designs are conceptually identical to those of HEM, but their structures, working mechanisms, and dynamic properties are essentially different [4][5][6][7][8].In particular, the lower chamber of the HEM is very thin, and the pressure in the lower chamber is usually smaller than that in the upper chamber.On the contrary, the two chamber pressures of the HDB are relatively high.Therefore, the LP models developed for HEM cannot be directly applied to HDB [4][5][6][7][8].
The disadvantage of the HDB with inertia track is that the flow velocity of the liquid column in the inertia track becomes very slow due to the viscosity of the fluid once the frequency excitation reaches its resonant frequency; that is, it usually can be assumed that the flow "shut off" [25][26][27][28], ( ( (1)  which is undesirable.However, the HDB with large orifice can provide sufficient damping in a wide frequency range, while the damping in the low frequency range is deteriorated [11].
The fluid tracks plates of HDBs usually have irregular geometry and vary from sample to sample significantly.In this paper, the HDBs with different fluid track plates are manufactured and studied, which can provide insights into the features needed for scientific verification.
The objectives of this paper include the following:

Design of HDB with Multiple Configurations and Experimental Study
2.1.Types and Experimental Procedure.A group of HDB types with varied configurations are manufactured to investigate the influences of the inertia tracks and orifices on the HDB performances.The basic features of the 10 tested HDBs are summarized in Table 2.For each type (HB1-HB5), two identical samples are prepared.Within these HBDs, the type of HB6 is special because it denotes the main rubber of the HDB.
Because the fluid of HB1 and HB3 does not provide hydraulic damping, the LP models are not considered in this study.Different LP models for those HDBs (HB1, HB2, HB4, and HB5) are illustrated in Figure 3.   is the complex stiffness of the main rubber and   is the damping ratio of the main rubber. 1 and  2 are the volumetric compliances of the two fluid chambers induced by the rubber spring and  1 and  2 are the volumetric stiffness of the two fluid chambers, respectively. 1 and  2 represent the effective pumping areas of the upper and the lower chambers, assuming that the behavior of the rubber spring can be approximated by Q i (t) Outer sleeve

Inner sleeve
Outer sleeve x r (t) Outer sleeve

Inner sleeve
Outer sleeve x r (t) x o (t) Outer sleeve

Inner sleeve
Outer sleeve x r (t) x o1 (t) x o2 (t) that of an artificial cylinder piston. 1 () and  2 () denote the pressures in the upper and the lower chambers.The inertia tracks and orifices are characterized by the fluid inertial coefficients   and   and fluid resistances   and   , respectively.The static preload and a sinusoidal displacement excitation,   () =   sin , are applied to the inner metal sleeve, while the outer metal sleeve is fixed.  () is the dynamic force which is transmitted to the outer metal sleeve by the hydraulic and main rubber paths [4][5][6][7][8].
The static and dynamic tests for the above-mentioned HDB with multiple internal configurations are conducted using elastomer test machine (MTS 831).Detailed descriptions of the test method and data processing technique for obtaining the dynamic properties of an HDB refer to [24].In addition, it is assumed that the relative displacement between the inner and outer metal sleeves varies in radial direction.Therefore, only radial damping is considered in this study.

Experimental Results.
Figure 4 shows the tests results of static stiffness of the HDBs with different internal configurations, that is, HB1, HB2, HB4, and HB5.It can be seen that the force-displacement curves for these 4 samples are linear and have similar gradients under small amplitude excitations (0∼ 3 mm).However, they tend to differ with each other and be nonlinear for large displacement excitations (>3 mm).When the displacement excitation is greater than 4 mm, the outer metal sleeve contacts with the stopper member; thus, the force-displacement curves exhibit significant nonlinearity.
In general, the static characteristics of these samples correlated well with each other.This indicates that the vulcanization process conditions for these samples are stable and consistent.
The results of the dynamic tests of the main rubber and the HDBs with multiple internal configurations under sinusoidal excitations are shown in Figure 5.

The Dynamic Performances of the Main Rubber Spring.
As shown in Figures 5(a) and 5(b), a sinusoidal excitation is applied to the inner sleeve of the main rubber (HB6) with a preload of 387 N. The changes of the dynamic stiffness and loss angle of HB6 with excitation frequency are very small compared to those of HB1, HB2, HB3, HB4, and HB5.Therefore, the dynamic stiffness and the loss angle of the HB6 can be reasonably treated as two constants, for example, taking the average of the test data (  is about 438340 N/m and   is about 102.06 Ns/m), in the subsequent analysis.

The Dynamic Performances of HDBs with Multiple Configurations.
Figure 5 shows that the dynamic performances of the HDBs with multiple configurations are strongly excitation frequency and amplitude dependent.
As the excitation frequency increases, the dynamic stiffness of both HB3 and HB6 can be considered as constant.The difference between them can be attributed to the incompressibility of the liquid, which causes a stiff.HDB.The dynamic stiffness of HDBs with multiple configurations tends to be a horizontal line under high frequency band, respectively, so they can be considered as amplitude invariant in this scenario.
Under the 0.8 mm excitation, the frequencies of the peak for the loss angles of HB1, HB2, HB4, and HB5 are 12, 20, 52, and 88 Hz, respectively.As the frequency increases, the loss angles of HB1, HB2, HB4, and HB5 gradually approach a horizontal line which overlaps the curves of HB3 and HB6.This phenomenon indicates that, in higher frequency range, the response of the fluid in tracks is roughly attenuated and the damping of the HDB is mainly determined by its main rubber.The bandwidths of the loss angle for the HDBs with different internal configurations are different.With the increase of the number of inertia tracks, the HDB can provide damping over a wider frequency band.While increasing the number of orifices, the HDB can provide larger damping in a higher and wider frequency band.
The high frequency dynamic performances of the HDBs with multiple internal configurations can be examined by observing Figures 5(c) and 5(d).As the number of orifices increases, the notch frequency    min shifts to a higher value and the level of the dynamic stiffness   decreases simultaneously.Moreover, it is shown that the corresponding loss angles are all 90 ∘ which means there is a resonance at    min .
It can be found that the loss angle shifts abruptly from about 0 ∘ to 180 ∘ .This means the resonance is dominated not only by the fluid flow resistance but also by the fluid inertial effect, suggesting that HDB with multiple configurations cannot be modelled by a single degree of freedom (SDOF) system under high frequency small amplitude excitation [25,27].Since the fluid motion is highly turbulent at higher frequencies, its influence is not modelled here.

The Fixed Points on Dynamic Performances of HDBs with Multiple Configurations.
It should be noted that many excitation amplitude-invariant fixed points   ( = 1 ∼ 4) and   ( = 1 ∼ 6) exist on the dynamic test results curves.With the number of inertia tracks increases, the dynamic stiffness (   ) of the fixed point is decreased.As the number of orifices increases,    decreases rapidly.These fixed points contain important information of an HDB, and a new parameter identification method can be developed.

The Nonlinear LP Model for an HDB with One Inertia Track
This section presents the modeling of the HDB with one inertial track illustrated in Figure 3(a).The following assumptions are made [29]: (i) the influence of gravity is ignored and the fluid in chamber is incompressible; (ii) the inertia and damping of the fluid in lower and upper chambers are neglected and the fluid pressure in each chamber is uniform; (iii) the relative velocity of fluid is constant along the inertia channel and the properties of inlet and outlet cross sections are uniform; (iv) the dynamic viscosity and density properties of the liquid are constant, and the cross sections along inertia track have the same shape.
With the above assumptions, the continuity equations can be written as (1) According to Bernoulli equation for nonstationary flow [30], the momentum equation of the fluid column mass in the inertia track is yielded as follows:

Shock and Vibration
where   is the major loss coefficient and the calculation details of   can be found in Wang et al. [29] and  is the total local loss factor including contraction loss coefficient of inlet  in , expansion loss coefficient of outlet  out , and the track bend local loss factor   determined by the curvature radius and diameter of fluid track according to empirical formulas [29]. in and  out are the cross-sectional areas of inlet and outlet, respectively.
The dynamic force transmitted to the outer sleeve can be expressed as where   () is the transmitted force of the rubber path and  ℎ () is the transmitted force of the hydraulic path.

Analytical Analysis of Dynamic Performance of an HDB
with One Inertia Track.Substituting (1) into (2a), (2b), and (2c), the equation of the liquid column flowing in the inertia track can be obtained Then, the resonance frequency of liquid column flowing in the inertia track can be derived as follows: Equation (5) indicates that the value of the resonant frequency   closely relates not only to the bulk stiffness values of two chambers, but also to ,   ,   .
Since | ẋ  | ẋ  is a quadratic function [31] of the liquid column velocity, the exact solution of (4) cannot be achieved.An alternative is to linearize the nonlinear term as | ẋ  | ẋ  ≈ 8   ẋ  /3 utilizing the describing function analysis [31].Then, (4) can be rewritten as where the damping ratio of the system is Substituting   =     ,   =    (− 1 ) , and   = F  (−) into ( 6), it gives Therefore, the transfer function of the fluid flowing in an inertia track can be written as where Using ( 1)-( 3) and ( 9), the sinusoidal complex stiffness can be obtained Thus, the complex stiffness can be rewritten as where  *  is the dynamic stiffness of the rubber spring and  *  represents dynamic characteristic induced by the subsystem of chambers and liquid in the inertia track;   represents the storage dynamic stiffness and   denotes the loss dynamic stiffness.According to (12) and ( 13), the dynamic performances of HDB with one inertia track can be investigated analytically.

When the Ratio of Excitation Frequency and Resonance
Frequency of Liquid Column Tends to 0. When the excitation frequency tends to zero, that is,  → 0,   and  *  become When  2 ≈  1 , then As shown in (14),  ( → 0) is positive, which indicates that the displacement excitation   and the corresponding response   are in-phase, and the liquid pumped by the equivalent piston area of the rubber spring is pumped into Shock and Vibration the other chamber in-phase.Thus, there is almost no pressure fluctuation in both chambers.Equation (15) means that the dynamic performance of the HDB is mainly contributed by the rubber spring.
In practice,  2 ≈  1 ,  *  is not equal to zero, and the value of  ( → 0) is slightly bigger than that of  ( → 0) .As the test results shown, when the excitation amplitude is 1.00 mm,  (=0) = 485.47N/mm,  (=0) = 434.14N/mm, and the ratio of  (=0) and  (=0) is 1.12; that is, they are almost the same.Therefore,   can be obtained by testing  *  of an HDB with the liquid under a very low frequency excitation.

When the Ratio of Excitation Frequency and Resonance
Frequency of Liquid Column Tends to 1. Using (9), when the excitation frequency approaches the resonance frequency (  ) of the inertia track subsystem, that is,  → 1,   and  *  become Substituting (17a) into (11) leads to Using (18), the loss angle of the HDB with one inertia track can be obtained as As shown in (10), (19), many parameters such as ,   , ,   , and   can be selected to adjust the loss angle.The smaller the excitation displacement amplitude   , the larger the   , and  (=1) became larger too.This is the reason why HDB has amplitude-variant dynamic performances.This conclusion matches well the experimental results shown in Figure 5(b).

When the Ratio of Excitation Frequency and Resonance
Frequency of Liquid Column Tends to ∞.Using ( 9) and ( 13), when the high excitation frequency is applied on an HDB with one inertia track, that is,  → ∞,   and  *  become As shown in (20), the displacement response of liquid column is drastically weakened; this effect is sometimes referred to as flow "shut-off" [27] under high excitation frequency.In this situation, the liquid pumped by the piston of the rubber spring is conformed all by two chambers due to their elasticity of bulk, which results in pressure fluctuations.Substituting ( 20) into (11), it can be obtained that As shown in (21), the storage dynamic stiffness under high excitation frequency is not varied, and  ( → ∞) is constant considering that the influence of the damping of the main rubber is negligible [28], which means that these parameters   ,  1 ,  2 ,  1 , and  2 can be regarded as constants.This analysis result agrees well with the test results shown in Figure 5(a).

Parameter Identification Method Based on the Fixed Points.
Comparing ( 18) with ( 21), then It can be seen that the storage dynamic stiffness tends to be a constant in higher frequency range and the values of   ( → 1) and   ( → ∞) are identical.This indicates that the storage dynamic stiffness curves for different amplitude excitations will intersect at   and have the value of   +  2  2 2 +  1  2 1 .As shown in Figure 6(a), all the tested curves of storage dynamic stiffness of HB1 pass through the fixed point  5 , and all the curves of storage dynamic stiffness of HB3 pass through the fixed point  5 approximately.All the curves almost approach a horizontal line in higher frequency range.Thus, the analytical conclusion is verified by test results, which implies the validity of the nonlinear LP model.
Set the excitation frequency at  5 as  5 and the corresponding storage dynamic stiffness as   (=1) .Based on the above analytical analysis, the frequency  5 of the fixed point ( 5 ) and   of the liquid column have the following relationship: Frequency (Hz)  Since the curves near the fixed point  5 are almost vertical to the frequency axis as shown in Figure 6(a), small experiment error of the storage dynamic stiffness will not cause a notable change in  5 .Therefore, the experimental data of  5 can be employed to identify   .Supposing  1 and  2 are proportional, the values of  1 and  2 can then be obtained according to (23).
Besides, the curves of storage dynamic stiffness are almost perpendicular to the longitudinal axis when  → ∞, so a higher variation of excitation frequency will only lead to a nearly smaller fluctuation of storage dynamic stiffness in high frequency range.On the contrary, the slopes of these curves around the resonance frequency (when  = 1) are almost vertical to the frequency axis; thus, a smaller fluctuation of excitation frequency will result in a significant variation of storage dynamic stiffness.So, the parameter identification results for  1 and  2 based on   ( → ∞) are more precise and credible.
According to the above method and Table 1, these parameters identifications are conducted and the results are summarized as follows:   = 4.86 × 10 5 (N/m),  2 = 2.31 × 10 11 (N/m 5 ),  1 = 2.08 × 10 11 (N/m 5 ),  1 = 1.91 × 10 −3 (m 2 ), and  2 = 1.72 × 10 −3 (m 2 ).Then, the dynamic performances of HB1 are calculated and compared with the experiment data, as shown in Figure 7.It is shown that the calculation results match well the experiment data, which Outer sleeve Inner sleeve Outer sleeve further validates the analytical analysis presented in this study.

The LP Models for HDB with Multiple Inertia Tracks and Orifices
As shown in Figure 5, the dynamic performances of HB2, HB4, and HB5 are very similar to that of HB1.This suggests that the LP model of HDB with one inertia track is also applicable for those with multiple inertial tracks.Thus, it will be adopted in the following analyses.

The LP Model for an HDB with Nonidentical Inertia Tracks.
An LP model of HDB with three nonidentical inertial tracks is illustrated in Figure 8.   ,   ( = 1, 2, 3) denote the lengths and the cross-sectional areas of the inertia tracks and  represents the inertial track number;  1 and  2 ( = 1, 2, 3) are the linear and nonlinear damping coefficients of the fluid motion in the inertia tracks, respectively.The inertia   and the linear and nonlinear resistances  1 and  2 ( = 1, 2, 3) of the fluid in the inertia track can be obtained by The fluid volume flow rate in three inertia tracks can be denoted as   ( = 1, 2, 3), respectively.The state equation of the LP model can be written as follows with the state vector defined as where In general, the nonlinear resistances  2 ( = 1, 2, 3) are negligible; that is,  2 = 0 ( = 1, 2, 3).Since   ( = 1, 2, 3) is finite, there will be some end effects due to contraction loss of inlet and expansion loss of outlet, and the tracks are curvilinear.These additional effects (including turbulence) are expected to add more damping to the HDB.
To compensate these effects, an empirical coefficient  is introduced to better approximate the actual system and to enhance the capillary tube resistance  1 ( = 1, 2, 3) as follows [32]: Outer sleeve

Inner sleeve
Outer sleeve x r (t) where  is the dynamic viscosity of the HDB fluid and A then becomes a constant matrix A  as Substituting Q  = ( 1 ,  2 ,  3 ) and  2 = 0 ( = 1, 2, 3) into (25), the motion equation of the lumped fluid mass can be written as where M, C, K, and F() are mass matrix, damping matrix, stiffness matrix, and force vector, respectively, and can be expressed as ] The characteristic equation of ( 29) is According to (31), the resonance frequency   can be obtained as This indicates that there is only one peak frequency for the loss angle, and the system has only one nonzero resonant frequency.

The LP Model for an HDB with Multiple Identical Inertia Tracks and Identical
Orifices.An LP model for an HDB with two groups of identical inertia tracks and identical orifices is shown in Figure 9.The inertia track group consists of  1 identical inertia tracks with the length of   and crosssectional area of   for each inertia track.Each inertia track is characterized by the fluid inertances   , linear fluid resistances  1 , and nonlinear resistances  2 ( = 1, 2, . . .,  1 ) of the lumped fluid mass.  can be thought of as the resistance to change in flow rate Q to the pressure gradient Δ applied on the fluid element.  of these elements is fundamentally the energy dissipation of fluid flow   over the pressure gradient Δ of the element.
Based on the steady flow assumption, the orifice-type models can also be used to describe the damping of orificetype tracks for   concentrated over a short distance.Because the orifice-type track has an aspect ratio of 2  /  ≈ 1, it can be taken as a short-tube orifice for dynamic performance prediction of HDB.The discharge coefficient   of such a short-tube orifice has the following empirical formulation [32] based on the Reynolds number (Re): Using the square-root relation [32,33] and taking into account the impact of the fluid jet, the fluid inertances   and the linearized resistance  1 ( = 1, 2, . . .,  2 ) of an orifice can be expressed as where  and  are empirical coefficients for adjusting HDB's inertances and damping in nonideal cases.  denotes the flow through the orifice [32].
Assuming the state variable of the LP model is X and is defined as X  = ( 1 ,  2 ,  3 ,  4 ) = ( 1 ,  2 ,   ,   ), the state equation then can be obtained as where ] The dynamic force transmitted to the outer sleeve by the rubber and hydraulic paths is Accordingly, the transfer function can be derived from (34), (35), and (37) in the Laplace domain as follows: where Using the same technique presented in Section 4.1, the peak frequency of the loss angle   can be explicitly expressed as The comparisons between the tested and calculated dynamic performances of HDBs with multiple identical inertia tracks and identical orifices under excitation amplitude of 0.8 mm and lower frequency band (1∼100 Hz) are presented in Figure 10.In general, the calculated results agree with the experiment data reasonably.Obvious errors for HB4 and HB5 can be seen around frequencies of 35 Hz and 55 Hz, respectively.This mainly attributes to the assumption that the fluid-jet length of the short-tube orifice can be regarded as the length of inertia track, which underestimates the impact of liquid turbulence in the orifice.However, the calculated dynamic performances of HB4 and HB5 can predict the dynamic performances of the HDBs in the whole frequency range.Therefore, the proposed model for an HDB with multiple identical inertia tracks and identical orifices is validated.

Model Verification and Discussion
Based on the geometric and material parameters of the HDB listed in Table 1, the peak frequency of the loss angle   of the HDB with multiple identical inertia tracks and identical orifices can be calculated according to (40).Table 3 shows the calculated results and the corresponding experiment data.It can be seen that they match well each other.Thus, the calculation method is validated.

Influences of Key Parameters.
Based on a lot of parametric analyses of the HDB with one inertia track, the same dynamic performances as that of the HDB with two identical inertia tracks can only be achieved by adjusting the length of the fluid track.HB1 l i = 0.4l i HB1 l i = 0.8l i HB1 l i = 4.0l i HB1 l i = 0.4l i HB1 l i = 0.8l i HB1 l i = 4.0l i As shown in Figure 11, when 0.4  <   < 0.6  , the curves of loss angle for HB1 approach that of HB2, and when 0.2  <   < 0.6  , the curves of dynamic stiffness for HB1 show consistency with that of HB2.This indicates that the same dynamic performances of an HDB with two identical inertia tracks cannot be achieved simultaneously by only changing the track length of an HDB with one inertia track.
As well known, the HDBs have irregular geometry and are constructed with alternate materials from elastomers to metals.It will be a very challenging work to reach an ideal dynamic performance by changing those parameters, considering the huge numbers of experiments.Varying the number of inertia tracks can be the best choice.
In general, an HDB with  1 ≥ 2 is more versatile for tuning since there are now many   and   to change, but as the number of inertia tracks increases, the maximum dynamic stiffness, the maximum loss angle, and the peak frequency of the loss angle also increase as shown in Figure 12, which is good for motion control of, say, the automotive suspension bounce mode.

Prediction of Dynamic
Performances of the HDBs.The influences of the number of inertia tracks on the dynamic performances of an HDB can be analytically investigated with the LP model proposed in Section 4.2.Assume that an HDB consists of  1 ( 2 = 0, 1 ≤  1 ≤ 7) identical inertia tracks.The influences of the number of inertia tracks on the dynamic performances of the HDB are shown in Figure 12.It is expected that an increase in the number of inertia tracks will introduce more fluid oscillating force contributing to   ().
Comparing the results for the inertia track cases of  1 = 2 to  1 = 7 shown in Figure 12, it shows that the maximum dynamic stiffness increases from approximately 2200 N/mm to 3200 N/mm.Additionally, the maximum loss angle increases from approximately 50 ∘ to 87 ∘ and the peak frequency of the loss angle increases from approximately 13.5 Hz to 52 Hz.Thus, the sensitivity of the HDB's performance increases near the maximum loss angle as well, as the increase of the number of inertia tracks.Further, the number of inertia tracks can be used to tune the device over a wider range of frequencies via increasing the peak frequency of the loss angle.
The influences of the track type on the performances of an HDB can be numerically analyzed with the proposed LP mode in Section 4.2 too.Assume that an HDB consists of  1 identical inertia tracks and  2 identical orifice tracks.The dynamic performances of these HDBs are shown in Figure 13.When comparing the results for track cases of  1 = 1 ( 2 = 2) and  1 = 2 ( 2 = 2) to that for the track case of  1 = 3 ( 2 = 2), it can be seen that the influences of the number of inertia tracks are not significant since the existence of orifices.
Comparing the results for the inertia track cases of  1 = 2 to  1 = 7 shown in Figure 12, it is shown that the maximum dynamic stiffness increases from approximately 2200 N/mm to 3200 N/mm.The maximum loss angle increases from approximately 50 ∘ to 87 ∘ and the peak frequency of the loss angle increases from approximately 13.5 Hz to 52 Hz.The HDB's performances near the maximum loss angle become more sensitive as the number of inertia tracks increases.Thus, increasing the number of inertia tracks can improve the device performance over a wider frequency range.
Observing the results for the track case of  1 = 0 ( 2 = 4) and  1 = 0 ( 2 = 6) shown in Figure 13, it can be seen that the maximum dynamic stiffness decreases from approximately 2300 N/mm to 2220 N/mm, but the dynamic stiffness increases from approximately 1050 N/mm to 1250 N/mm at a very low frequency, for example, 1 Hz.In addition, the maximum loss angle decreases from approximately 55 ∘ to 54 ∘ and the peak frequency of the loss angle increases from approximately 154 Hz to 196 Hz.In a word, as the number of orifices increases, the maximum dynamic stiffness and loss angle decrease slightly and the peak frequency of the loss angle increases significantly.Moreover, the HDB with more orifices can provide large damping in a wide frequency range.Therefore, the influences of the number of orifices on the performances of the HDBs are remarkable.

Conclusions
Hydraulically damped bushing (HDB) is crucial for the handling and NVH performances of vehicle.In this study, the dynamic stiffness and the loss angle of the HDBs with multiple internal configurations, which are the essential characteristics of the HDBs, are tested and the influences of

Figure 1 :
Figure 1: The cross-sectional view of HDB.

Figure 2 :
Figure 2: The fluid track plate of hydraulic bushing.
(a) design a controlled experiment to test the static and dynamic performances of HDB with multiple inertia tracks and orifices by change of the fluid track plate and investigate the excitation amplitude-invariant fixed points on the dynamic characteristics of HDBs with different fluid track plates; (b) present a new parameter identification technique and develop a nonlinear unified LP model for the HDB with multiple inertia tracks and orifices to analyze the influences of the number of different tracks on the dynamic characteristics of the HDB and estimate the peak frequency of the loss angle; (c) validate the unified analytical model by comparing simulation results with test data and predict the influences of the number of fluid tracks on the dynamic performances of the HDB.

Figure 3 :
Figure 3: Fluid system model for HDBs with multiple configurations.(a) One inertia track; (b) two inertia tracks; (c) one inertia and an orifice; (d) two orifices.

Figure 4 :
Figure 4: The tests results static characteristics of the HDBs with different internal configurations.

2 Figure 5 :
Figure 5: The dynamic test results of HDBs with different internal configurations.(a) Dynamic stiffness under low frequency band, large amplitude excitations; (b) loss angle under low frequency bands, large amplitude excitations; (c) dynamic stiffness under high frequency band, small amplitude excitations; (d) loss angle under high frequency band, small amplitude excitations.

Figure 6 :FrequencyFigure 7 :
Figure 6: The test results of storage and loss dynamic stiffness of the hydraulic bushing.(a) Storage stiffness.(b) Loss stiffness.

Figure 8 :
Figure 8: Fluid system model for an HDB with three nonidentical inertia tracks.

Figure 9 :
Figure 9: Fluid system model of hydraulic bushing with multiple identical inertia tracks and identical orifices.

5. 1 .
Model Verification.As shown in Figure9, when  2 = 0 and  1 = 1, the LP model is the same as the fluid system model of an HDB with one inertia track shown in Figure3(a); when  2 = 0 and  1 = 2, the LP model depicts the fluid system model of an HDB with two inertia tracks shown in Figure3(b); when  2 = 1 and  1 = 1, the LP model represents the fluid system model of an HDB with one inertia track and one orifice element shown in Figure 3(c); when  2 = 2 and  1 = 0, the LP model denotes the fluid system model of an HDB with two orifices shown in Figure 3(d).

Figure 12 :
Figure 12: Calculated dynamic performances of HDB with inertia tracks increasing.(a) Dynamic stiffness.(b) Loss angle.

Table 1 :
Geometric and material parameters of a typical HDB.Parameter Value Length of inertia track   /m 0.19 Cross-sectional area of the track   /  /m 2 9.83 × 10 −6 Wet perimeter of the inertia track   /m 1.27 × 10 −2 Hydraulic diameter of the inertia track   /m 3.087 × 10 −3

Table 2 :
The configurations of HDB samples.

Table 3 :
Comparison of the liquid column resonance frequency.