Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation

Liao, Xin; Sun, Xiaojuan; Wang, Hu

doi:10.21595/jve.2021.21921

Journal of Vibroengineering

Browse Journal

Submit article

Published: 04 May 2021

Check for updates

Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation

Xin Liao¹

Xiaojuan Sun²

Hu Wang³

¹School of Mechanical Engineering, Shijiazhuang Tiedao University, Shijiazhuang, China

¹State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures, Shijiazhuang Tiedao University, Shijiazhuang, China

²School of Mechanical Engineering, Taiyuan University of Science and Technology, Taiyuan, China

³14th Institute of China Electronic Technology Corporation, Nanjing, China

Corresponding Author:

Xiaojuan Sun

Cite the article Download PDF

Abstract

In order to improve the ride comfort of the cab of construction machinery and reduce the harm of vibration to the driver, the lumped parameter model of passive hydraulic damping rubber mount (PHDRM) considering the nonlinearity of multi-inertial tracks is first established. Then the overall FE (finite element) fluid-solid coupling model of the vibration isolator is proposed by using ABAQUS combined simulation mode in this paper. Meanwhile, Combined with bidirectional strong coupling algorithm, the parameters of the rubber main spring and inertial track are identified respectively. Also, the dynamic characteristic experiments of the PHDRM samples are carried out. The experiment results show that the analytical solution and FE modeling method of the HDRM system are reliable and accurate under larger amplitude excitation. On the basis of this, the effects of the rubber main spring and inertial track of PHDRM on low frequency dynamic characteristics under larger amplitude excitation are respectively analyzed and discussed. It provides practical reference value for nonlinear modeling method research and dynaimc performance analysis of PHDRM with complex structure further.

1. Introduction

In the recent years, as modern construction and transportation equipment, the construction machinery, such as loader, earthmover, etc., has had a profound impact on various fields of national economic construction and developed rapidly. The idea of “people-oriented” has been paid attention more and more, especially the cab will be subjected to different levels of excitations from the road surface during the operation of construction machinery [1-3]. If the driver has been always in the forced vibration working environment for a long time, it will cause high fatigue [4, 5] and lead to frequent accidents, and the driver's health will be seriously damaged. However, the traditional research mainly focuses on the vibration isolation of engines and vehicle mounted equipment [6-8], and the driver's physical feelings are ignored. Therefore, the research and design of cab suspension system has been highly valued by all walks of life year by year. Improving the bad working environment of drivers is an important development direction of cab suspension system of construction machinery.

Related studies have shown that the human body is a complex resonance system, especially sensitive to low frequency vibration [9-13]. When the vibration and shock caused by various excitations of road surface and engine are transmitted to the cab through the frame, it will cause complex vibration coupling characteristics. In order to improve the vibration of the cab, many studies are devoted to the design of semi-active vibration isolators. However, for large heavy-duty construction machinery, the manufacturing cost of vehicles has been greatly increased. The passive hydraulic damping rubber mount (PHDRM) has remarkable nonlinear characteristics, and its stiffness and damping characteristics are obviously superior to those of the traditional rubber or spring vibration isolators [14, 15]. It is a kind of vibration absorber with excellent performance, which is more used in construction machinery cab suspension system. At present, there is still a lack of research on the nonlinear modeling method of PHDRM, and its complex flow-solid coupling mechanical properties are the key research contents for reducing cab vibration. According to the literature [16], the experimental research and the parameter identification are respectively carried out for two kinds of hydraulic damping suspension with different shapes and sizes of orifice, and the dynamic characteristics of orificing disk hydraulic damping suspension are analyzed. Compared with the experimental results, the correctness of the mechanism of orificing disk is verified. A lumped parameter model for hydraulic damping suspension with multi-inertial channels is established in the literature [17] in order to analyze dynamic characteristics. The dynamic stiffness and the calculation method of its hysteresis angle are deduced, and the formula of peak frequency of hysteresis angle is also given. The literature [18] has studied the relationship between design parameters of a hydraulic suspension system and dynamic stiffness, hysteresis angle. And by parameter identification technique, the nonlinear lumped parameter model is established and the influences of fluid track on the dynamic characteristics of the system are analyzed. The results show that the simulation datum are consistent with the experimental results, and the parameter identification method is effective. However, from another perspective, the lumped parameter modeling method is used to establish the hydraulic damping suspension model, which cannot calculate the distribution of liquid flow in fluid track and pressures in different chamber of the vibration isolator. But the flow-solid coupling FE simulation analysis method can effectively avoid the above problems, and is widely used in engineering and has achieved good results. Many scholars have adopted the FE method in the structural design of hydraulic damping rubber mount. A nonlinear finite element model of rubber material of hydraulic damping suspension system is established in the literature [19]. The dynamic equivalent elastic modulus of elastic material is determined by iterative model updating program. The nonlinear FE model of hydraulic damping suspension is established by using INTESIM strong liquid- solid FE method in literature [20]. The static and dynamic stiffness characteristics of the vibration isolator are analyzed, and the validity of the simulation results is verified. Also, the literature [21] describes a two-dimensional fluid-solid coupling FE model for hydraulic damping rubber mount and static analysis and structural optimization have been carried out. By using nonlinear FE method, the performance parameters of PHDRM structure are identified in the literature [22] and the superplastic constitutive relation is used to accurately describe the incompressible properties of rubber materials, and then the dynamic behaviors of two different hydraulic damping suspensions are analyzed. However, in the existing studies, fewer reports have been found on the research and calculation of the overall FE modeling method for PHDRM. Therefore, it is necessary to establish an accurate simulation model of hydraulic damping suspension system for construction machinery cab, and the low frequency dynamic characteristics are analyzed according to the main parameters of the vibration isolator.

In this paper, the low frequency dynamic characteristics of passive hydraulic damping rubber mount under larger amplitude excitation are examined. A lumped parameter model of passive hydraulic damping rubber mount considering the nonlinearity of multi-inertial tracks is proposed in the following section. Then, the overall FE fluid-solid coupling model of the vibration isolator is developed. Meanwhile, based on the established FSI (Fluid-structure-Interaction) models of rubber main spring and inertial track, the main parameters are respectively identified. The identified results are applied to the nonlinear lumped parameter model with multi-inertial tracks above. Moreover, the low frequency larger amplitude excitation experiments of the PHDRM samples are investigated, while validating good agreement with the above simulation results. Finally, the influences of the rubber main spring and inertial track of PHDRM on low frequency dynamic characteristics under larger amplitude excitation are analyzed and discussed in detail.

2. Modeling of cab hydraulic damping suspension system

Due to strong vibration or shaking produced during operation for the cab of construction machinery, it is very difficult to improve the bearing capacity of hydraulic damping rubber mount by redesigning the structural dimensions and weights of the cab, especially for medium and heavy construction machinery. Fig. 1 shows the 3D model diagram of the inertial track-decoupled disk hydraulic damping rubber mount, which consists of rubber main spring, inertial track, rubber bottom membrane, decoupled disk and outer metal skeleton, etc. When hydraulic damping rubber mount is working, the top of the vibration isolator is connected to the cab, and the bottom of the vibration isolator is connected to the frame, and the external vibration makes the liquid reciprocate quickly in the inertial track by the rubber spring. The secondary suspension system of construction machinery is widely considered to be mainly low frequency vibration isolation. In addition, the structures of vibration isolator internal components such as rubber spring and inertial channel are optimized, which greatly improves the bearing capacity of the vibration isolator and has the low frequency and large damping characteristics.

Fig. 1The 3D model of passive hydraulic damping rubber mount (after simplied)

2.1. Modeling of lumped parameters of hydraulic damping suspension system

In order to analyze and compare the low frequency dynamic characteristics of hydraulic damping rubber mount, a lumped parameter model considering the nonlinearity of multiple inertial tracks is established, as shown in Fig. 2. The main parameters in the model include the dynamic stiffness $K_{r}$ of the rubber main spring, and equivalent viscous damping coefficient $B_{r}$ , which has very small influence on the dynamic characteristics of the vibration isolator by experiments and can be negligible in the low frequency range, and also the equivalent piston area $A_{e}$ , the volume stiffness of the upper liquid chamber $K_{1} = {1 / C}_{1}$ , the volume stiffness of the lower liquid chamber $K_{2} = {1 / C}_{2}$ , the inertia coefficient $I_{i}$ and the flow damping coefficient $B_{i}$ of the liquid flowing through the inertial track. The system variables in the lumped parameter model are respectively the average pressure of the upper liquid chamber and the lower liquid chamber $P_{1} (t)$ and $P_{2} (t)$ , the liquid rate of flow through the inertial track $Q_{i} (t)$ , the displacement excitation at the input end $x (t)$ and the support reaction at the output end $F (t)$ .

Fig. 2Lumped parameter model of hydraulic damping rubber mount considering the nonlinearity of multi-inertial tracks

When the external excitation frequency is low and the amplitude is large, the decoupling disk of hydraulic damping rubber mount is in the unstarted state in most of the time, thus it is approximate that the liquid flows only through the inertial track between the upper and lower liquid chambers. Considering that the inertia force and damping force of liquid flowing through inertial track are nonlinear, so the nonlinear momentum equation of liquid flow in inertial track can be derived as follows:

1

P_{1} (t) - P_{2} (t) = I_{i} {\dot{Q}}_{i} (t) + (B_{i 1} + B_{i 2} |Q_{i} (t)|) Q_{i} (t) .

Due to the incompressibility of the liquid, the relationship between the motion velocity of the rubber main spring equivalent piston and the flow velocity of the liquid in the inertial track can be expressed as:

2

Q_{i} (t) = A_{e} {\dot{x}}_{d} = S_{i} {\dot{x}}_{i},

where, $x_{i}$ is the fluid displacement in the inertial track.

The continuity equation of liquid flow can be given as:

3

{\dot{P}}_{2} (t) = K_{2} Q_{i} (t) .

By substituting Eq. (2) and Eq. (3) into Eq. (1), the nonlinear momentum equation of liquid flow in inertial track is then given by:

4

P_{1} (t) - P_{2} (t) = I_{i} A_{e} {\ddot{x}}_{d} + B_{i 1} A_{e} {\dot{x}}_{d} + B_{i 2} A_{e}^{2} |{\dot{x}}_{d}| {\dot{x}}_{d},

where, $|{\dot{x}}_{d}| {\dot{x}}_{d}$ is quadratic function of vibration velocity of rubber main spring equivalent piston, which can be linearized by methods in the literature [23]:

5

\begin{array}{l} |{\dot{x}}_{d}| {\dot{x}}_{d} = (X_{d} ω)^{2} |c o s ω t| c o s ω t \approx \{\frac{ω}{π} \int_{0}^{2 π / ω} |c o s ω t| c o s ω t \cdot c o s ω t \cdot d t\} (X_{d} ω)^{2} c o s ω t \\ = \frac{8}{3 π} X_{d} ω {\dot{x}}_{d}, \end{array}

where, $X_{d}$ is the displacement amplitude of rubber main spring equivalent piston.

Thus, Eq. (4) is rewritten as:

6

P_{1} (t) - P_{2} (t) = I_{i} A_{e} {\ddot{x}}_{d} + B_{i 1} A_{e} {\dot{x}}_{d} + B_{i 2} A_{e}^{2} \frac{8}{3 π} X_{d} ω {\dot{x}}_{d} .

The force equilibrium equation of the suspension system is establihsed as follows:

7

\{\begin{array}{l} F (t) = K_{r} x + B_{r} \dot{x} + A_{e} P_{1} (t), \\ A_{e} P_{1} (t) = K_{b} (x - x_{d}), \end{array}

where, $K_{b}$ is the equivalent linear stiffness of the upper liquid chamber, $K_{b} \approx K_{1}$ .

The state equation of the system is established by nonlinear momentum equation and continuity equation of system:

8

\dot{X} = X A + B_{i},

where:

A = [\begin{matrix} 0 & 0 & - K_{b} & \dots & - K_{b} \\ 0 & 0 & K_{2} & \dots & K_{2} \\ {1 / I}_{i} & - 1 / I_{i} & - B_{i} / I_{i} & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋱ & ⋮ \\ 1 / I_{i} & - 1 / I_{i} & 0 & \dots & - B_{i} / I_{i} \end{matrix}],

X^{T} = [\begin{array}{l} P_{1} & P_{2} & Q_{i 1} & Q_{i 2} & \dots & Q_{i N} \end{array}], {B_{i}}^{T} = [\begin{array}{l} A_{e} {\dot{x}}_{d} + C_{1} {\dot{P}}_{1} & 0 & 0 & 0 & \dots & 0 \end{array}],

B_{i} = B_{i 1} + B_{i 2} |Q_{i} (t)| .

Letting ${Q_{i}}_{1} (t) + {Q_{i}}_{2} (t) + \dots + {Q_{i}}_{N} (t) = Q_{i}$ , according to the above equations, by the Las transformation and substituting Eq. (8) into Eq. (7), the dynamic reaction force of the hydraulic damping rubber mount is obtained as follows:

9

F (s) = K_{r} X (s) + s B_{r} X (s) + A_{e}^{2} K_{b} \{\frac{s^{2} I_{i} + s (B_{i 1} + B_{i 2} A_{e} \frac{8}{3 π} X_{d} ω) + N K_{2}}{s^{2} I_{i} + s (B_{i 1} + B_{i 2} A_{e} \frac{8}{3 π} X_{d} ω) + N K_{2} + N \frac{K_{b}}{A_{e}^{2}}}\} X (s) .

Since the rubber bottom membrane is thin, $K_{2} \approx 0$ , thus the complex stiffness of the isolation system is given by:

10

K = \{K_{r} + A_{e}^{2} K_{b} [\frac{ω^{4} I_{i}^{2} + ω^{2} ({B_{i}}^{2} - \frac{N I_{i} K_{b}}{A_{e}^{2}})}{(- ω^{2} I_{i} + \frac{N K_{b}}{A_{e}^{2}})^{2} + {B_{i}}^{2} ω^{2}}]\} + j ω \{B_{r} + \frac{N B_{i} K_{b}^{2}}{(- ω^{2} I_{i} + \frac{N K_{b}}{A_{e}^{2}})^{2} + {B_{i}}^{2} ω^{2}}\} .

When $N = 1$ , Eq. (7) can be considered as the theoretical formula of lumped parameter model with the single inertial track. The above lumped parameters are required to be identified further by the FEM method in the next Part. Then the identified results can be substituted into Eq. 10, and the low frequency dynamic characteristics of the hydraulic damping rubber mount can be obtained. The analytical results are shown in Fig. 15-16.

2.2. FE Modeling of hydraulic damping suspension system

The liquid in the passive hydraulic damping rubber mount will reciprocate between the upper chamber and lower chamber, and the interaction of fluid and solid exsits in the chamber of the vibration isolator. Therefore, based on fluid-solid coupling method, the overall finite element model of the hydraulic damping rubber mount is established by combined simulation mode in ABAQUS. Due to the small deformation of rubber main spring metal skeleton, it usually can be neglected. Because the rubber main spring is the typical superplastic material, the deformation characteristics within 80 % of rubber can be better described by using the Mooney-Rivilin material model in ABAQUS, which also has good stability. The strain energy function of the Mooney-Rivilin material model can be expressed as follows:

11

U = \sum_{i + j = 1}^{N} C_{i j} (I_{1} - 3)^{i} (I_{2} - 3)^{j} + \sum_{i = 1}^{N} \frac{1}{D_{1}} (\sqrt{I_{3}} - 1)^{2 i} .

When $N =$ 1, the strain energy function of second order Mooney-Rivilin material model can be given by:

12

U = C_{10} (I_{1} - 3) + C_{01} (I_{2} - 3) + \frac{1}{D_{1}} (J_{e l} - 1)^{2},

where, $C_{10}$ , $C_{01}$ and $D_{1}$ are respectively the temperature-related material parameters, $I_{1}$ and $I_{2}$ are first and second order strain invariant, respectively. $J_{e l}$ is the elastic volume ratio, and assuming that the rubber material is isotropic and incompressible, the influence of elastic volume ratio can be neglected. The parameters of rubber material are shown in Table 1. Other structures are set to be steel materials.

Table 1The parameters of rubber material

Rubber density	Superelasticity
$ρ$ (Kg/m³)	C₁₀	C₀₁	D₁
0.97×10^-9	0.44	0.11	0.001

In the simulation environment, the assembly position of the fluid model and the solid model in the same coordinate system should be guaranteed to match each other. Therefore, the fluid model of simplified vibration isolator is shown in Fig. 2, which is divided into three parts: the fluid model in the upper chamber, the fluid model in the inertial track, and the fluid model in the lower chamber. Table 2 shows the material parameters of the fluid finite element model. Then, the overall FE modeling of hydraulic damping rubber mount is carried out by ABAQUS. First, the FE models of solid and fluid are respectively established, which can be discretized using different element and mesh sizes. Secondly, the flow-solid coupling interface is defined as follows: (1) the contact surfaces between the fluid in the upper chamber and rubber main spring, metal skeleton and track cover, (2) the contact surfaces between the fluid in the lower chamber and rubber bottom membrane, track cover, (3) the contact surfaces between the fluid in the inertial track and track cover. Each pair of flow-solid coupling surfaces are respectively selected in the solid model and the fluid model. In order to make the FE mesh size of the fluid model and the solid model basically consistent and improve the calculation efficiency, each pair of contact surfaces in the fluid and the solid model should be created with the same area during the pretreatment process. When the location of the fluid and solid nodes on the flow-solid coupling surface can not be coincided, the displacement of the fluid node is obtained by the displacement interpolation of the solid node, and the stress of the solid node is obtained by the pressure interpolation of the corresponding fluid node. The information on the structural materials and FE unit nodes of each part of hydraulic damping rubber mount are shown in Table 3.

Table 2The parameters of fuild material

Density of Silicone oil / T·(mm³)^-1	Dynamic viscosity / MPa·s
1.1×10^-9	5.6×10^-8

Table 3Information on the structural materials and FE unit nodes of hydraulic damping rubber mount

Structure	Material	Unit length	Unit type	Number of units	Number of nodes
Rubber main spring	Rubber	6	Tetrahedron	22364	11327
Inertial track	Steel 45	8	Tetrahedron	13791	3865
Rubber bottom membrane	Rubber	6	Hexahedron	3662	4201
Liquid	Silicone oil	4	Hexahedron	83694	15712
Decoupled disk	Steel 45	8	Tetrahedron	4938	4010
Track cover	Steel 45	8	Tetrahedron	2941	1085

Because of the incompressibility of the fluid in the hydraulic damping rubber mount, thus the combined simulation mode in ABAQUS should be created to carry out the fluid-solid coupling calculation. When modeling, the dynamic analysis module and fluid dynamics analysis module are combined with each other together. The solid FE model is applied in the ABAQUS/Explicit Model and the fluid FE model is applied in the ABAQUS/CFD Model. By discretizing the dynamic equations of fluid and solid and introducing coupling boundary conditions, the nonlinear FE equation of fluid-solid coupling system can be yielded:

13

F (X) = [\begin{array}{l} F_{f} [X_{f}, d_{s} (X_{s})] \\ F_{s} [X_{s}, τ_{f} (X_{f})] \end{array}] = 0,

where, $F_{f}$ and $F_{s}$ are respectively the FE equations corresponding to the equations of motion of fluid and solid, $X = {(X_{f}, X_{s})}^{T}$ , $X_{f}$ is the vector consisting of the velocity vector of nodes in the fluid FE model and the pressure vectors of the elements, $X_{s}$ is the vector composed of displacement vectors of nodes in the solid FE model, $d_{s}$ and $τ_{f}$ are respectively the displacements of the structures and the stresses of the fluid on the coupling boundary.

Fig. 3Schematic diagram of fluid model

3. Parameter identification of hydraulic damping suspension system

3.1. Parameter identification of rubber main spring

The FE mesh of rubber main spring is built by using ten node quadratic tetrahedron elements C3D10H, and the displacement load of 5 mm is carried out. The displacement cloud map of the rubber main spring is obtained, as shown in Fig. 4. The total reaction force can be obtained by extracting support reaction of all the nodes subjected to the displacement load at the top of the vibration isolator. By using the total reaction force and the displacement as the horizontal and vertical coordinates respectively, the force-displacement curve of the rubber main spring is obtained as shown in Fig. 5. Thus, the static stiffness of the rubber main spring in the vertical direction is calculated as 236.71 N/mm, and its dynamic stiffness can be estimated by the following formula:

14

K_{r} = f ∙ k_{s t a t i c} .

Normally, when $f =$ 1.4, the dynamic stiffness is calculated as 331.39 N/mm.

Fig. 4Displacement cloud map of rubber main spring

Fig. 5Force-displacement curve of rubber main spring

When the rubber main spring produces displacement under preload, the fluid model in the upper chamber is added into the rubber spring solid model. By using the combined simulation mode, a pressure “head” is set for the fluid and the fluid will flow down the bottom surface under pressure. So the bottom pressure of the fluid is set to 0 MPa, which makes the fluid flows down. Also, experiments show that when the excitation displacement exceeds 3 mm, the equivalent piston area of the rubber main spring is basically unchanged and independent on the excitation frequency. Fig. 6 shows the displacement cloud map of the fluid model in the fluid-solid coupled model, from which the values of the displacement of fluid downward $x_{d}$ , the area of fluid bottom $A_{d}$ , and the upper end displacement of solid $x$ can be obtained respectively. According to Eq. (15):

15

A_{e} = \frac{A_{d} \cdot x_{d}}{x},

the equivalent piston area of rubber main spring is obtained $A_{e} =$ 4035 mm². Fig. 7 shows that the equivalent piston area of rubber main spring after stable motion is about 4008.48 mm². The volume stiffness of upper liquid chamber is written as:

16

K_{1} = Δ P / Δ V,

where, $Δ P$ is the variable quantity of fluid pressure for inner chamber of rubber main spring, $Δ V$ is the variable quantity of volume in the inner chamber. By loading displacement, the FSI model of rubber main spring can be solved to obtain the volume stiffness of the upper liquid chamber $K_{1} =$ 4.21×10¹⁰ N/m⁵.

Fig. 6Displacement cloud map of the fluid model (FSI)

Fig. 7Equivalent piston area curve over time

3.2. Parameter identification of inertial track

The fluid FE model in the inertial track is meshed with tetrahedron element FC3D4. The inlet cross section of the inertial track is set to the velocity load surface, and the outlet cross section pressure is set to 0 MPa, and the other surfaces are set to the fluid wall. The velocities of fluid inflow are respectively given as 682 mm/s, 1364 mm/s, 2046 mm/s, 2728 mm/s, 3410 mm/s, 4093 mm/s, 4775 mm/s. Fig. 8 shows the pressure distribution cloud map and velocity cloud map of inertial track at inlet fluid velocity of 4775 mm/s by simulation and analysis. The relationship curve between pressure difference and cross section flow at each flow rate is shown in Fig. 9. The momentum equation is given by:

17

Y = U β,

where:

Y = [\begin{matrix} Δ p^{1} \\ Δ p^{2} \\ \dots \\ Δ p^{n} \end{matrix}], U = [\begin{matrix} \dot{Q_{1}} & Q_{1} & |Q_{1}| Q_{1} \\ \dot{Q_{2}} & Q_{2} & |Q_{2}| Q_{2} \\ \dots & \dots & \dots \\ \dot{Q_{n}} & Q_{n} & |Q_{n}| Q_{n} \end{matrix}], β = [\begin{matrix} I_{1} \\ B_{i 1} \\ B_{i 2} \end{matrix}] .

$Δ p^{n}$ and $Q_{n}$ represent the pressure difference and flow rate between the inlet and outlet of the inertial track respectively, and by MATLAB calculation the flow damping coefficient of inertial track can be obtained as $B_{i} =$ 1.1×10⁸N⋅s/m⁵. The inertia coefficient of liquid flowing through inertial track is given as $I_{i}$ = 2.58×10⁶kg/m⁴, the cross section area for inertial track is given as $S_{i} =$ 80 mm², and the length of inertial track is given as $L =$ 102.6 mm.

Fig. 8Simulation of inertial track at inlet fluid velocity of 4775 mm/s

a) Pressure distribution diagram

b) Flow rate diagram

Fig. 9Relationship curve of flow and pressure difference in inertial track

4. Analysis of dynamic performance of hydraulic damping suspension system

4.1. Dynamic characteristic analysis of flow-solid coupling

The main dynamic characteristics of passive hydraulic damping rubber mount include dynamic stiffness and hysteretic angle, which can be obtained by dynamic reaction-displacement curve of hydraulic damping suspension system during one vibration period.

The equation of dynamic stiffness, hysteretic angle and damping coefficient are respectively as follows:

18

K_{d} = \frac{F}{D},

19

b = \frac{K_{d} s i n (ϕ)}{ω_{0}},

20

ϕ = s i n^{- 1} \frac{S}{π F D},

where, $F$ are $D$ are the force and displacement amplitudes, respectively, and $S$ is the area of hysteresis loop.

For the simulation of dynamic characteristics, the vertical displacement load is first applied on the top surface of the rubber main spring:

21

x (t) = x_{0} + A s i n (2 π f t),

where, $x_{0}$ is the preload displacement for hydraulic damping suspension system, relating to the static working load of the vibration isolator, $A$ is the dynamic displacement excitation amplitude, $f$ is the excitation frequency. By using ABAQUS multi-field coupling solution platform, the collaborative analysis task can be created based on Dell workstation with ZQ E5-2665 2.40 GHz 8-core processor. Then, the curve of the displacement of the active end of rubber main spring and dynamic reaction force is obtained, as shown in Fig. 10(a). It can be seen that there is hysteresis in the relationship between dynamic reaction force and displacement, which indicates that the suspension system has a larger damping hysteretic angle. Compared with Fig. 10(b), it is found that the hysteresis produced at 10 Hz is slightly larger than that at 20 Hz, and the damping hysteresis angle of the system is smaller at 20 Hz excitation frequencies.

Fig. 10The hysteretic curves of rubber main spring and dynamic reaction force for hydraulic damping suspension system

a) $x_{0} = -$ 0.05 mm, $A =$ 0.2 mm, $f =$ 10 Hz

b) $x_{0} = -$ 0.05 mm, $A =$ 0.2 mm, $f =$ 20 Hz

The motion states and deformations of the fluid at any time during the whole excitation process can be observed by the visual results of nonlinear FE calculation. Fig. 11 shows the pressure distribution cloud diagram of the fluid. It can be seen that there is obvious pressure gradient in inertial track at different times in one cycle, and the pressure distributions of the upper and lower liquid chambers are basically uniform. When $t =$ 0.15 s, the fluid pressure of the upper and lower liquid chambers changes obviously. The liquid flows between the upper and lower chambers through the inertial track, and the pressures at the inlet and outlet of the inertial track changes with the reciprocating motion of the fluid. The maximum flow velocity is 1472 mm/s, as shown in Fig. 12. Fig. 13 shows the load velocity curve applied to the top surface of the rubber main spring. It indicates that the maximum velocity is 26.84 mm/s. Thus, the ratio of the two velocities is about 54.84, which is not different from the ratio of the equivalent piston area of the rubber main spring to the cross section area of the inertial track. The correctness of the FE model can be verified.

Fig. 11The pressure distribution diagrams of Fluid Model at different Times (A= 0.2 mm, f= 10 Hz)

a) $t =$ 0.11 s

b) $t =$ 0.14 s

c) $t =$ 0.15 s

d) $t =$ 0.16 s

e) $t =$ 0.17 s

f) $t =$ 0.18 s

g) $t =$ 0.19 s

h) $t =$ 0.20 s

Fig. 12The velocity of flow distribution in inertial track (t= 0.15 s)

Fig. 13Load velocity curve under displacement excitation

4.2. Experiment verification

In order to further verify the effectiveness of lumped parameter modeling and overall FEM modeling, the dynamic performance experiments of the samples of passive hydraulic damping rubber mount were carried out. By using the MTS elastomer test platform, the working end of the test platform is connected with the connecting end of the bottom plate of the cab, and the vibration isolator is fixed on the test platform, as shown in Fig. 14. According to the total mass of the cab, the steady-state harmonic displacement excitation with the amplitude of $A =$ 0.2 mm is first employed, and the preload displacement is 4 mm, and the excitation frequency range is 0-50 Hz. The displacement signal at the working end of the test platform and the reaction signal of the base of the vibration isolator are respectively recorded. According to Eqs. (18-20), the experimental datum cane be processed. Finally, the experimental results are shown in Fig. 15. It indicates that under the low frequency excitation, the dynamic characteristic curves of the hydraulic damping rubber mount by lumped parameter modeling and overall FEM modeling have good consistency with the experimental results, respectively, which illustrates above modeling methods are credible. Also, it can be found that there are still some deviations between the analytical solution of the dynamic characteristics and the nonlinear FE solution.

When the construction machinery is subjected by the road excitation with low frequency and large amplitude, the secondary suspension system of its cab will also be excited by the larger amplitude. Therefore, the larger amplitude excitation experiment of the vibration isolator is conducted. Fig. 16 shows the dynamic characteristics of hydraulic damping rubber mount under larger amplitude excitation. It illustrates the above two modeling methods still have reliability.

Fig. 14Mechanical performance test platform of vibration isolator

a) Schematic illustration of the vibration isolation clamping

b) Test site

Fig. 15Dynamic Characteristics of hydraulic damping rubber mount when A= 0.2 mm

a) Dynamic stiffness

b) Damping coefficient

Fig. 16Dynamic characteristics of hydraulic damping rubber mount when A= 2.0 mm

a) Dynamic stiffness

b) Damping coefficient

4.3. Effect of rubber main spring on low frequency dynamic characteristics under larger amplitude excitation

To deeply analyze the effects of rubber spring on the low frequency dynamic characteristics on the vibration isolator under larger amplitude excitation, the simulation analysis with setting frequency 5-50 Hz, 5 Hz per interval, $A =$ 2.0 mm are carried out based on above fluid-solid coupling model. And the dynamic stiffness characteristics and damping coefficient at each excitation frequency are obtained, as shown in Fig. 17-18. The peak values of dynamic stiffness of hydraulic damping suspension have reached the maximum near 30 Hz, and then have decreased gradually, and also have increased with the increase of dynamic stiffness and equivalent piston area of rubber main spring. At the excitation frequency of 22 Hz, the damping coefficient of the vibration isolator reaches the peak, and then tends to decrease gently. Because of the poor working and transportation environment of construction machinery, hydraulic damping rubber mount can produce better vibration isolation effect with the inhibition of strong vibration or shaking of cab under larger amplitude excitation.

Fig. 17Effect of different dynamic stiffness of rubber main springs on dynamic characteristics of hydraulic damping rubber mount

Fig. 18Effect of different equivalent piston areas on dynamic characteristics of hydraulic damping rubber mount

a) Dynamic stiffness

b) Damping coefficient

4.4. Effect of inertial track on low frequency dynamic characteristics under larger amplitude excitation

In order to analyze the influence of the physical parameters of the inertial track on the low frequency larger amplitude dynamic characteristics of the vibration isolator, the four most important factors, namely, the cross section area of inertial track, the length of inertial track, the number of inertial track, and the damping coefficient of inertial track flow, are respectively changed on the basis of the original model. Multiple groups of the dynamic stiffness and damping coefficient curves at different excitation frequencies are obtained by the above modeling method, as shown in Figs. 19-22, respectively. It can be seen that the peak values of dynamic stiffness of the hydraulic damping rubber mount have become larger with the increase of the cross section area and the number of the inertial track, while the peak variations of damping coefficient have shown the opposite result, and the peak frequencies of the damping coefficients have been gradually increasing with the increase of the number of inertial tracks. Therefore, in the structural design of the vibration isolator, the influence of the number of inertial tracks on the low frequency dynamic characteristics and the manufacturing cost should be considered synthetically. When the length of inertial track has increased gradually, the peak values of dynamic stiffness of the vibration isolator have first increased, and then decreased, and its peak frequencies has been decreasing. The results illustrates that the length of inertial track in a certain range will change the low frequency dynamic characteristics of hydraulic damping rubber mount under larger amplitude excitation. If the length value is too large, the stiffness will be small and the vibration isolation performance will become poor although its peak frequencies will be reduced. Since the sensitive frequency range of human body is 0.5-80 Hz according to ISO-2631, the driver's ride comfort should be considered fully for the vibration isolation performance of the cab hydraulic suspension system. If the length of the inertial track is slightly small, it would cause the peak frequency to be too large. Low frequency vibration isolation performance will also be reduced, especially under larger amplitude excitation. Meanwhile, the peak values of damping coefficient of the vibration isolation have gradually increased with the increase of length, and the peak frequencies have decreased continuously.

Fig. 19Effect of different cross section areas of inertial tracks on dynamic characteristics of hydraulic damping rubber mount

a) Dynamic stiffness

b) Damping coefficient

Fig. 20Effect of different number of inertial tracks on dynamic characteristics of hydraulic damping rubber mount

a) Dynamic stiffness

b) Damping coefficient

Fig. 21Effect of different length of inertial tracks on dynamic characteristics of hydraulic damping rubber mount

a) Dynamic stiffness

b) Damping coefficient

Fig. 22Effect of different damping coefficients of inertial tracks on dynamic characteristics of hydraulic damping rubber mount

a) Dynamic stiffness

b) Damping coefficient

Moreover, It can be found from Fig. 22 that when the damping coefficient of the inertial track has been changed, the overall variation amount of the dynamic stiffness and damping coefficient curve of the vibraiton isolator under the low frequency larger amplitude is not very large. And the peak values of the dynamic stiffness, the peak values of the damping coefficients and their peak frequencies have been all decreased with the increase of the damping coefficient of the inertial track. Therefore, when the liquid viscosity of hydraulic damping rubber mount is selected, if the viscosity is smaller, it can not produce obvious damping effect when flowing through the inertial track, which affects the stiffness of hydraulic damping suspension system and reduces the vibration isolation performance. But if the viscosity is larger, However, if the viscosity is larger, it will also hinder the movement of fluid, or block the inertial track, which can not play a good vibration isolation effect.

5. Conclusions

In this work, the modeling method and nonlinear dynamic characteristic analysis for PHDRM of construction machinery cab are studied. Considering the nonlinearity of multi-inertial tracks, the nonlinear lumped parameter model of PHDRM is developed to predict the suspension system's low frequency dynamic characteristics under larger amplitude excitation. Meanwhile, an overall FE model of PHDRM with fluid-solid coupling is established, and the parameters of the rubber main spring and inertial track are identified respectively by using ABAQUS combined simulation mode. In addition, the experiment results of the PHDRM samples are in good agreement with the simulation results under low frequency larger amplitude excitation. Furthermore, the effects of the main parameters of PHDRM on low frequency dynamic characteristics under larger amplitude excitation are also respectively addressed. The analysis and discussion provide the analysis basis for the design and evaluation of the secondary suspension system of the construction machinery cab.

References

Kordestani A., Rakheja S., Marcotte P., et al. Analysis of ride vibration environment of soil compactors. SAE International Journal of Commercial Vehicles, Vol. 3, Issue 1, 2010, p. 259-272.

Publisher
Yu Man, Zhou Chen Yu, Wei Lang, et al. Takagi-sugeno fuzzy-model-based automobile electrohydraulic active suspension H-infinity control. China Journal of Highway and Transport, Vol. 31, Issue 8, 2018, p. 205-217.

Search CrossRef
Sun X. J., Zhang J. Optimization of low-frequency vibration isolation for cab ride comfort of construction machineries. Journal of Mechanical Engineering, Vol. 21, 2012, p. 52-60.

Search CrossRef
Burstrom L., Aminoff A., Bjor B., et al. Musculoskeletal symptoms and exposure to whole-body vibration among open-pit mine workers in the Arctic. International Journal of Occupational Medicine and Environmental Health, Vol. 30, Issue 4, 2017, p. 553-564.

Publisher
Liu C., Qiu Y., Griffin M. J., et al. Dynamic forces over the interface between a seated human body and a rigid seat during vertical whole-body vibration. Journal of Biomechanics, Vol. 61, 2017, p. 176-182.

Publisher
Zeng X., Liette J., Noll S., et al. Analysis of motor vibration isolation system with focus on mount resonances for application to electric vehicles. SAE International Journal of Alternative Powertrains, Vol. 4, Issue 2, 2015, p. 370-377.

Publisher
Zhang Hui Jie, Hao Hui Rong, Guo Zhi Ping Analysis and optimization design of six degrees of freedom engine mount system. Journal of Vibration Engineering, Vol. 32, Issue 5, 2019, p. 801-810.

Search CrossRef
Wang Min, Yao Guo Feng, et al. A novel design of semi-active hydraulic mount with wide-band tunable notch frequency. Journal of Sound and Vibration, Vol. 333, Issue 8, 2014, p. 2196-2211.

Publisher
Deprez K., Moshou D., et al. Improvement of vibrational comfort on agricultural vehicles by passive and semi-active cabin suspensions. Computers and Electronics in Agriculture, Vol. 49, Issue 3, 2005, p. 431-440.

Publisher
Morales A. L., Nieto A. J., et al. A semi-active vehicle suspension based on pneumatic springs and magnetorheological dampers. Journal of Vibration and Control, Vol. 24, Issue 4, 2018, p. 808-821.

Publisher
Bovenzi M., Rui F., Negro C., et al. An epidemiological study of low back pain in professional drivers. Journal of Sound and Vibration, Vol. 298, Issue 3, 2006, p. 514-539.

Publisher
Maciejewski I., Krzyzynski T., et al. Shaping the vibro-isolation properties of horizontal seat suspension. Journal of Low Frequency Noise Vibration and Active Control, Vol. 36, Issue 3, 2017, p. 203-213.

Publisher
Ning Dong Hong, Sun Shuai Shuai, Du Hai Ping, et al. Integrated active and semi-active control for seat suspension of a heavy duty vehicle. Journal of intelligent material systems and structures, Vol. 29, Issue 1, 2018, p. 91-100.

Publisher
De T. J., Deprez K., Anthonis J., et al. Conceptual cab suspension system for a self-propelled agricultural machine – part 2: operator comfort optimization. Biosystems Engineering, Vol. 90, Issue 3, 2005, p. 271-278.

Publisher
Xu Pei Jun, Bill B., Ken T., et al. Optimal mounting design for cab vibration isolation. International Journal of Vehicle Design, Vol. 57, Issue 2, 2011, p. 292-304.

Publisher
Fan Rang Lin, Lv Zhen Hua Nonlinear dynamic characteristics and parameters identification method for hydraulic engine mount. Journal of Mechanical Engineering, Vol. 43, Issue 7, 2007, p. 145-151.

Publisher
Shangguan Wen Bin, Song Zhi Shun, Zhang Yun Qing, et al. Experimental study and simulation analysis of hydraulic engine mounts with multiple inertia tracks. Journal of Vibration Engineering, Vol. 18, Issue 3, 2005, p. 318-323.

Search CrossRef
Yang C., Yin Z., Shangguan W., et al. A study on the dynamic performance for hydraulically damped rubber bushings with multiple inertia tracks and orifices: parameter identification and modeling. Shock and Vibration, Vol. 2016, 2016, p. 3695950.

Publisher
Jahani K., Dehnad M. Identifying the frequency dependent material property of a hydraulic engine mount through an iterative procedure using 3D finite element modeling. Journal of Mechanical Science and Technology, Vol. 28, Issue 6, 2014, p. 2041-2047.

Publisher
Cao Zheng Lin, Li Jun, Guo Kong Hui, et al. Simulations of dynamic characteristics of hydraulic engine mounts by FEM with strong fluid-structure coupling. Journal of Vibration and Shock, Vol. 31, Issue 10, 2012, p. 4-8.

Search CrossRef
Sun X. J., Zhang J. Low-Frequency shaking of the construction machine cab and design of vibration isolators for protecting. Applied Mechanics and Materials, Vol. 141, 2011, p. 49-53.

Publisher
Christopherson J., Jazar G. N. Dynamic behavior comparison of passive hydraulic engine mounts. Part 2: finite element analysis. Journal of Sound and Vibration, Vol. 290, Issue 3, 2006, p. 1071-1090.

Search CrossRef
Christopherson J., Nakhaie Jazar G. Dynamic behavior comparision of passive hydraulic engine mounts. Part 2: Finite element analysis. Journal of Sound and Vibration, Vol. 290, 2006, p. 1071-1090.

Publisher

Cited by

Time domain response analysis of multi-flow channel hydraulic mount

(2024)

Measurement and optimization of nonlinear damping systems for agricultural engineering vehicle cab

Xin Zhang | Yuanyou Liu | Zhanlong Li | Zengliang Xiao

(2024)

Structural Improvement in Crane Boom Rest

CA Thulasi Dharan | U Sudalai Vadivoo | SA Mohammed Hussain Basith | Pathmarajan R | P Ganesh

(2024)

Structural Improvement in Crane Boom Rest

CA Thulasi Dharan | U Sudalai Vadivoo | SA Mohammed Hussain Basith | Pathmarajan R | P Ganesh

(2024)

Effect of bell plate structure on high- and low-frequency characteristics of hydraulic mount

Zhihong Lin | Yuedong Huang

(2024)

Genetic algorithm-based shape optimization of rubber mount for tractor-cabin mounting system

Min-woo Kang | Hyun-woo Han | Jung-ho Park | Seung-je Cho | Young-jun Park

(2024)

Advances in Engineering Research and Application

Canh Chi Huan | Dang Viet Ha | Le Anh Vu | Tran Van Thoan | Le Van Quynh

(2023)

Advances in Engineering Research and Application

Le Van Quynh | Le Anh Vu | Bui Van Cuong | Hoang Anh Tan | Le Xuan Long

(2022)

Modeling the Dynamics of a Gyroscopic Rigid Rotor with Linear and Nonlinear Damping and Nonlinear Stiffness of the Elastic Support

Zharilkassin Iskakov | Kuatbay Bissembayev | Nutpulla Jamalov | Azizbek Abduraimov

(2021)

About this article

Received

16 February 2021

Accepted

09 April 2021

Published

04 May 2021

SUBJECTS

Mechanical vibrations and applications

DOI

https://doi.org/10.21595/jve.2021.21921

Keywords

mechanical engineering

hydraulic damping rubber mount

flow-solid coupling combined simulation

nonlinear dynamic characteristics

rubber main spring

inertial track

Acknowledgements

This work has been funded by The National Natural Science Foundation of China (No. 11902207), the Natural Science Foundation of Hebei Province (A2020210018), Education and Teaching Reform Project of Shijiazhuang Tiedao University (Y2020-15).

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2023 12 05

Time domain response analysis of multi-flow channel hydraulic mount

Zhihong Lin, Chunrong Wang, Yunxiao Chen, Mingzhong Wu, Feijie Zheng, Jian Wei

Research article

2023 11 27

Analysis of the influence of the interaction between decoupled membrane channels and inertial channels on the mount' characteristics

Zhihong Lin, Yunxiao Chen, Mingzhong Wu, Feijie Zheng

Research article

2021 12 02

Effect of mounting and fluid parameters on dynamic damping characteristics of a hydraulic damper

Wenlin Wang, Rongyong Li, Shan Zhu, Yongming Wu

Research article

2021 06 28

The experimental study of elastic-hysteresis properties of rubber elements of sleeve-pin couplings

Dmitry Nasonov, Vladimir Ilichev, Vladimir Raevsky

X. Liao, X. Sun, and H. Wang, “Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation,” Journal of Vibroengineering, Vol. 23, No. 3, pp. 542–558, May 2021, https://doi.org/10.21595/jve.2021.21921

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2021.21921
UR  - https://doi.org/10.21595/jve.2021.21921
TI  - Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation
T2  - Journal of Vibroengineering
AU  - Sun, Xiaojuan
AU  - Liao, Xin
AU  - Wang, Hu
PY  - 2021
DA  - 2021/05/04
PB  - JVE International Ltd.
SP  - 542-558
IS  - 3
VL  - 23
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

@article{Sun_2021,
	doi = {10.21595/jve.2021.21921},
	url = {https://doi.org/10.21595/jve.2021.21921},
	year = 2021,
	month = {may},
	publisher = {{JVE} International Ltd.},
	volume = {23},
	number = {3},
	pages = {542--558},
	author = {Xiaojuan Sun and Xin Liao and Hu Wang},
	title = {Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]X. Sun, X. Liao, and H. Wang, “Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation,” Journal of Vibroengineering, vol. 23, no. 3, pp. 542–558, May 2021, doi: 10.21595/jve.2021.21921.

Copy IEEE

Copied to clipboard!

Sun, Xiaojuan, Xin Liao, and Hu Wang. “Modeling and Dynamic Analysis of Hydraulic Damping Rubber Mount for Cab under Larger Amplitude Excitation.” Journal of Vibroengineering 23, no. 3 (May 4, 2021): 542–58. https://doi.org/10.21595/jve.2021.21921.

Copy Chicago

Copied to clipboard!

Modeling and dynamic analysis of hydraulic damping rubber mount for cab under larger amplitude excitation

Abstract

1. Introduction

2. Modeling of cab hydraulic damping suspension system

2.1. Modeling of lumped parameters of hydraulic damping suspension system

2.2. FE Modeling of hydraulic damping suspension system

3. Parameter identification of hydraulic damping suspension system

3.1. Parameter identification of rubber main spring

3.2. Parameter identification of inertial track

4. Analysis of dynamic performance of hydraulic damping suspension system

4.1. Dynamic characteristic analysis of flow-solid coupling

4.2. Experiment verification

4.3. Effect of rubber main spring on low frequency dynamic characteristics under larger amplitude excitation

4.4. Effect of inertial track on low frequency dynamic characteristics under larger amplitude excitation

5. Conclusions

References

Cited by

About this article

Related Articles