Influence of hydraulic engine mounts on engine shake based on full vehicle model
Rong Guo^{1} , Jun Gao^{2} , Xiaokang Wei^{3}
^{1, 2, 3}Clean Energy Automotive Engineering Center, Tongji University, Shanghai, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 10, 2016, p. 376-381.
Received 8 November 2016; accepted 9 November 2016; published 8 December 2016
JVE Conferences
This work proposes a full vehicle model including hydraulic engine mounts (HEMs) to better describe the characteristics of the vehicle engine shake performance. The model consists of 14 degree of freedoms (DOFs), namely 6 of the powertrain, 1 of the fluid within the inertia track of the HEM, 3 of the car body and 4 of the unsprung mass. Simulation based on the model is performed to demonstrate how the powertrain mounting system including the HEM influences the frequency response of vehicle subjected to harmonic excitations from road.
Keywords: engine shake, full vehicle model, hydraulic engine mount.
1. Introduction
Present automotive industry is witnessing a neck to neck fight to produce highly developed models with better performance. One of the performance requirements is accurate engine shake models to give better evaluation of the influence on the ride comfort caused by powertrain systems. Some popular models for vehicle ride analysis have been implemented by many researchers. Wen-Bin Shangguan et al. [1] established a vehicle model with 13 Degree of Freedoms (DOFs). Based on the model, the vibration at the seat track is calculated with input from road profile. Feng Wang [2] compared and analyzed the vibration characteristics based on a 6 DOF powertrain model and 13 DOF model with the powertrain system. The simulation results based on the latter model coincided with the experiment results better.
However, the models they set up are based on the assumption that the powertrain are on the elastic support connected to the vehicle body, which means, only rubber mounts are taken into account. It is known that hydraulic engine mounts (HEMs) are widely used in cars, therefore it is necessary to study the vibration response of the vehicle body based on the engine shake models with the mounting system including HEMs.
Vilsecker J. [3] considered the powertrain and the body as two subsystems, with the HEM between them. Tao C et al. [4] established a 3 DOF model with powertrain system, and then demonstrated that HEMs could better cancel the vehicle vibration. Fukazawa M. et al. [5] developed a CAE method to predict the vibration transfer function of the HEM on a vehicle with sufficient precision. The models including HEMs and vehicle body are simple and not adequate for more accurate analysis.
This work focuses on providing more accurate engine shake model by replacing a rubber mount in the 13 DOF vehicle model with a HEM. Simulation based on the model is performed to evaluate the effect of the HEM on the engine shake performance by comparing the results obtained from the known 13 DOF model.
2. Full vehicle model with powertrain for engine shake analysis
2.1. Mechanical model of the HEM
The linear mechanical model of HEM for the vertical vibration analysis under low frequency and large amplitude excitations is established, as shown in Fig. 1. The parameter ${m}_{i}$, ${b}_{i}$, ${A}_{i}$, ${X}_{i}$ denotes the effective fluid mass in the inertia track, the damping coefficient for the linear fluid motion, the area of the inertia track, and the linear displacement of the fluid in the inertia track. ${A}_{p}$, ${k}_{r}$, ${b}_{r}$ represents equivalent area of the upper chamber, stiffness and damping properties of the main spring. The compliance of upper chamber and lower chamber is defined as ${C}_{1}$, ${C}_{2}$. The transmitted force to the force is derived for the mechanical model to ensure this system is mathematically equivalent to the original system equations. Observing the transmitted force at each pin, the force equations are determined:
${F}_{{T}_{3}}={k}_{2}\frac{{A}_{i}}{{A}_{i}}{X}_{i},{F}_{{T}_{4}}={b}_{i}{\dot{X}}_{i}.$
Fig. 1. Low frequency, large amplitude linear model
Fig. 2. Full vehicle model including HEMs
2.2. Full vehicle model for engine shake including hydraulic engine mount
A 14 DOFs full vehicle model shown in Fig. 2 for the analysis of engine shake is proposed in this work. The model includes 6 DOFs of the powertrain, 3 DOFs of car body, 4 DOFs of the four unsprung mass, and 1 DOF of the fluid within the inertia track of the HEM.
In static equilibrium status, the following coordinate systems are defined. The Vehicle Coordinate System (VCS) is denoted with ${O}_{o}$-${X}_{o}{Y}_{o}{Z}_{o}$; the Powertrain Coordinate System (PCS) is denoted ${O}_{p}$-${X}_{p}{Y}_{p}{Z}_{p}$; the Car Body Coordinate system (BCS) is denoted with ${O}_{b}$-${X}_{b}{Y}_{b}{Z}_{b}$. The origins of ${O}_{p}$-${X}_{p}{Y}_{p}{Z}_{p}$ and ${O}_{b}$-${X}_{b}{Y}_{b}{Z}_{b}$ are located at the center of gravity (CoG) of powertrain and car body, respectively. The coordinate systems for the unsprung mass (UCS)are denoted as ${O}_{ui}$-${X}_{ui}{Y}_{ui}{Z}_{ui}$ ($i=$1,2, 3, 4) with its origin located at the CoG of each unsprung mass. The ${X}_{o}$- and ${Y}_{o}$- axis in the VCS are parallel to the horizontal plane, and ${Z}_{o}$- axis is normal to the horizontal plane and the positive direction is upward. The positive ${X}_{o}$- axis points to the rear of the vehicle. The ${X}_{p}$-, ${X}_{b}$-, ${X}_{i}$ and ${X}_{uj}$- axes are parallel to ${X}_{o}$- axis; the ${Y}_{p}$-, ${Y}_{b}$-, ${Y}_{i}$- and ${Y}_{uj}$- axes are parallel to the ${Y}_{o}$- axis; the ${Z}_{p}$-, ${Z}_{b}$-, ${Z}_{i}$ and ${Z}_{uj}$- axes are parallel to the ${Z}_{o}$- axis, respectively. In static equilibrium, each dynamic coordinate system coincides with its fixed coordinate system. A right-hand Mount Coordinate System (MCS) denoted as ${o}_{i}$-${u}_{i}{v}_{i}{w}_{i}$_{}($i=$1, 2, 3, 4) is built for each mount, with its origin ${o}_{i}$ at the center of the mount. The stiffness and damping of the rubber mounts in ${u}_{i}$-, ${v}_{i}$- and ${w}_{i}$- directions are ${k}_{ui}$, ${k}_{vi}$, ${k}_{wi}$ and ${c}_{ui}$, ${c}_{vi}$, ${c}_{wi}$, while ${k}_{s1}$, ${k}_{s2}$, ${k}_{s3}$, ${k}_{s4}$ and ${c}_{s1}$, ${c}_{s2}$, ${c}_{s3}$, ${c}_{s4}$ are the stiffness and damping of the suspensions. ${k}_{t1}$, ${k}_{t2}$, k_{t3}. k_{t4} denotes the stiffness of the tires.
The displacement excitation from ground to wheels is defined as ${\mathbf{q}}_{h}={\left({z}_{h1},{z}_{h2},{z}_{h3},{z}_{h4}\right)}^{T}$, and the displacements of CoGs for powertrain ${m}_{p}$, car body ${m}_{b}$ and unsprung mass ${m}_{ui}$ are defined as ${\mathbf{q}}_{p}={\left({x}_{p},{y}_{p},{z}_{p},{\theta}_{xp},{\theta}_{yp},{\theta}_{zp}\right)}^{T}\text{,}$${\mathbf{q}}_{b}\mathbf{}={\left({z}_{b},{\theta}_{xb},{\theta}_{xb}\right)}^{T}$ and ${\mathbf{q}}_{u}={\left({z}_{u1},{z}_{u2},{z}_{u3},{z}_{u4}\right)}^{T}\text{,}$ respectively, where $T$ denotes the transpose of a vector. The displacement vector is defined as $\mathbf{q}={\left[{\mathbf{q}}_{p}^{T},{\mathbf{q}}_{b}^{T},{\mathbf{q}}_{u}^{T},{z}_{i}\right]}^{T}$. The component in ${\mathbf{q}}_{h}$ denotes the excitation to left front (LF) wheel, right front (RF) wheel, left rear (LR) wheel and right rear (RR) wheel, respectively. Displacement of the unsprung mass $i$ is defined as ${z}_{ui}$. The component in ${\mathbf{q}}_{b}$ represents the vertical displacement in $Z$ direction and angular displacement around $X$ and $Y$ directions; ${\theta}_{xp}$, ${\theta}_{yp}$, ${\theta}_{zp}$ are the angular displacement around the $X$, $Y$, $Z$ directions respectively in the PCS. The coordinate system for the fluid in the inertia track of the HEM is defined as ${O}_{i}$-${X}_{i}{Y}_{i}{Z}_{i}$ with its origin located at the COG of the fluid mass.
2.2.1. The kinetic energy
The kinetic energy of the whole system consists of translational and rotational energy of vehicle body, as well as the translational energy of the wheels. The kinetic energy of the powertrain can be written as:
$+\frac{1}{2}\cdot {I}_{zz}{\dot{\theta}}_{zp}^{2}-\frac{1}{2}\cdot {I}_{xy}{{\dot{\theta}}_{x}}_{p}{\dot{\theta}}_{yp}-\frac{1}{2}\cdot {I}_{yz}{\dot{\theta}}_{yp}{\dot{\theta}}_{zp}-\frac{1}{2}\cdot {I}_{zx}{{\dot{\theta}}_{z}}_{p}{{\dot{\theta}}_{x}}_{p},$
where the moments of inertia (${I}_{xx}$, ${I}_{yy}$, ${I}_{zz}$) and inertia products (${I}_{xy}$, ${I}_{yz}$, ${I}_{xz}$) of a powertrain are defined in the PCS.
Similarly, the kinetic energy of the body can be expressed as follows:
where ${J}_{x}$, ${J}_{y}$ denote moments of inertia of the body. The kinetic energy of the unsprung mass become:
Assuming that there are $N$ HEMs in the mounting system, then the kinetic energy of the fluid in the inertia track ${m}_{i}$ is:
where ${z}_{i}$ represents the displacement of the fluid. The total kinetic energy can be calculated:
Eq. (6) can be expressed in the quadratic form:
2.2.2. The potential energy
The potential energy of the mounting system is:
The displacement of ${m}_{i}$ in ${O}_{i}$-${X}_{i}{Y}_{i}{Z}_{i}$ are defined as ${x}_{mi}$, ${y}_{mi}$, ${z}_{mi}$, then the potential energy of the HEM become:
$+1/2\cdot {k}_{2}{\left(\frac{{A}_{i}}{{A}_{p}{z}_{i}}\right)}^{2}=\frac{1}{2}\cdot {\mathbf{q}}^{T}{\mathbf{K}}_{2}^{\text{'}}\mathbf{q}.$
Similarly, the potential energy of the entire suspension system can be obtained:
The total potential energy of the tire system becomes:
As a result, the total potential energy of the 14DOF system proposed can be derived:
Also, the stiffness matrix of the 14DOF system can be derived:
2.2.3. The dissipation energy
The dissipation energy of the rubber mounts can be written as:
The dissipation energy of the HEM can also be derived:
The dissipation energy of the suspension system is:
Therefore, the total dissipation energy of the 14DOF system proposed can be described as:
The stiffness matrix of the 14DOF system becomes:
2.2.4. Response of the seat track
The transfer function between the inputs and outputs can be derived:
where $\mathbf{H}\left(w\right)$ is a (14×4) matrix, each column of which represents the response of the 14 DOF to the excitation at corresponding tires. Therefore, the transmissibility of vehicle body is $\left|\mathbf{H}\left(w\right)\right|$. Then the acceleration response magnitude becomes:
Then acceleration response magnitude at the seat track becomes:
where ${\mathbf{H}}_{a}^{7}\left(w\right)$, ${\mathbf{H}}_{a}^{8}\left(w\right)$, ${\mathbf{H}}_{a}^{9}\left(w\right)$ denote the 7th, 8th, and 9th row of ${\mathbf{H}}_{a}\left(w\right)$, respectively. ${\mathbf{r}}_{xst}$, ${\mathbf{r}}_{yst}$ represents the $X$ axis and $Y$ axis coordinates of seat track in the VCS.
3. Simulation and analysis
The acceleration response magnitudes of the vehicle body under different inputs are calculated, as demonstrated in Fig. 3.
Fig. 3. Acceleration response magnitude of the vehicle body under different inputs: a) RF wheel input; b) LF wheel input; c) RR wheel input; d) LR wheel input
Fig. 3(a)-(b) illustrate that the magnitude of bounce resonance in the full vehicle model without HEMs is much higher than that in the 14 DOF model. After the engine rubber mount is replaced by HEM, the magnitude peak value of the acceleration response to the RF wheel and LF wheel input decrease by approximately 24.81 % and 33.34 % respectively. The HEM with high resistance is responsible for the restraint of the resonance, while the damping coefficients of rubber mounts are too small to mitigate the resonance. The resonant frequencies switch from 7.83 Hz to 10.19 Hz and from 7.81 Hz to 10.29 Hz respectively.
Fig. 3(c)-(d) show that the HEM has little influence on the response magnitude when the input is from RF wheel or RR wheel, for the reason that the curves are so close to each other. Therefore, the following analysis neglects evaluating the acceleration response to LR wheel and RR wheel inputs. Based on the above analysis, we derive the acceleration responses of seat track under right front and LR wheel inputs, which are shown in Fig. 4.
In Fig. 4(a), the bounce resonance at 7.6 Hz, dominates the curve when the powertrain is mounted without HEMs. However, when the mounting system features the HEM, the bounce resonance is severely restrained by up to 68.78 %. The resonance frequency switches to approximately 10 Hz. Compared with the green curve, the blue curve fluctuates more smoothly with a maximum value of 81.64. When the HEM is adopted, the maximum value in Fig. 4(b) is 3.97 times as much as that in Fig. 4(a). It can be seen that the resonance of the green curve at 7.74 Hz is limited by 19.18 % due to the damping of the HEM.
Fig. 4. Acceleration response magnitude of the seat track under different inputs: a) RF wheel input; b) LF wheel input
4. Conclusions
In this work, a full vehicle model with 14 DOFs including the fluid in the inertia track of the HEM is proposed to calculate the engine shake at vehicle body and seat track. It is concluded that the mounting system with the HEM has greater restriction on the body and seat track shake subjected to front wheel inputs, compared with the mounting system simply with rubber mounts. Furthermore, the response magnitude due to LF wheel input is much higher than that due to RF wheel input for the reason that the HEM is at the right hand side of the powertrain.
References
- Shangguan Wen-Bin, et al. Design method of automotive powertrain mounting system based on vibration and noise limitations of vehicle level. Mechanical Systems and Signal Processing, Vol. 76, 2016, p. 677-695. [Search CrossRef]
- Wang Feng, Jin Y. J., Zhang J. W. Vibration simulation and optimization of a powertrain mounting system based on a full vehicle model. Journal of Vibration and Shock Vol. 27, Issue 4, 2008, p. 134-138. [Search CrossRef]
- Vilsecker J., Mack W., Falkner A., et al. On the complexity of hydromount-models in full vehicle simulation. Proceedings in Applied Mathematics and Mechanics, Vol. 11, Issue 1, 2011, p. 341-342. [Search CrossRef]
- Tao C., Zhu H., Xu P., et al. Frequency-dependent hydraulic engine mount with five-parameters fractional derivative model in vehicle model. SAE Technical Papers, 2015. [Search CrossRef]
- Masahiro Fukazawa, Murao T., Unigame S. The method to predict the vibration transfer function of hydraulic engine mount on a vehicle. Aquichan, Vol. 4, Issue 4, 2009, p. 70-72. [Search CrossRef]
- Ms Trupti Bhoskar, et al. Genetic algorithm and its applications to mechanical engineering: a review. Materials Today: Proceedings, Vol. 2, Issue 4, 2015, p. 2624-2630. [Search CrossRef]