Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel

Wang, Qiang; Zhao, Youqun; Lin, Fen; Fu, Hongxun; Du, Xianbin

doi:10.21595/jve.2016.16952

Journal of Vibroengineering

Browse Journal

Submit article

Published: 31 December 2016

Check for updates

Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel

Qiang Wang¹

Youqun Zhao²

Fen Lin³

Hongxun Fu⁴

Xianbin Du⁵

^{1, 2, 3, 4, 5}College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China

Corresponding Author:

Youqun Zhao

Cite the article Download PDF

Downloads 1591

WoS Core Citations 11

CrossRef Citations 11

Abstract

This paper presents the vibration characteristics and its key influencing factors of a new mechanical elastic wheel (MEW). The MEW was modeled as a ring on elastic foundations (REF) with distributed spring stiffness in the radial and tangential directions. The general forced solutions of inextensible vibration were derived by the use of a modal expansion technique and Arnoldi method, and the accuracy of the solutions had been validated by FEM simulation and modal test under free suspension and various loading situations. The natural frequencies and mode shapes of the rotating MEW could be obtained under free suspension and various loading situations. Moreover, the effects of various rotational speed, loading and different number of hinges on natural frequencies were investigated. Finally, the effect of different number of hinges on the damping ratio of the MEW radial modes was also analyzed. The analysis results reflect the objective law of the actual vibration characteristics of the MEW, and provide a reference for the MEW structure optimization and the vibration characteristics of the whole vehicle.

1. Introduction

To change the existing pneumatic tire puncture damage and high speed blowout etc., developing run-flat and anti-puncture tire to guarantee high performance and security has become a consensus of the world’s major tire manufacturers. Therefore, the researchers have recently focused their attention on non-pneumatic tires with different structures. Based on that, the new mechanical elastic wheel (MEW) as a non-pneumatic safety tire structure is proposed. The MEW can be realized the basic function of traditional pneumatic tire, in additional, the problems such as stinging, puncturing and blasting damage have been avoided. Thus, the MEW is greatly satisfied with the requirements of safety performance for the special vehicles. The tire is the main component of vibration reduction, and its dynamic characteristics directly affect the driving smoothness and riding comfortableness of vehicle to a certain extent. So that the good dynamic characteristics are also one of the key factors for tire design. The structural vibration parameters of tire can also predict the interaction between tire and road excitation, suspension system and so on, and it can provide guidance for the optimization design of tire structure and suspension system.

The free and forced vibration of stationary or rotating rings has been widely studied for various ring models and for various boundary conditions. I. F. Kozhevnikov [1-2] investigated the vibration of a tyre using the model of a wheel with a reinforced tyre under various loading and various rotating situations. The non-linear vibration of a rotating ring was investigated by W. B. Bickford [3] and W. Kim [4]. S. C. Huang [5-6] derived the general equations of motion that governed both transverse and circumferential motions of rotating rings and expanded the solution to forced response. S. J. Kim [7] and Y. T. Wei [8] investigated the forced response for rotating tires under various loading using the REF model. C. R. Dohrmann [9] and S. C. Huang [10] investigated the dynamic analysis of a tire-wheel-suspension assembly using an inextensible circular ring on a foundation connected to the wheel. Meanwhile, D.H. Guan et al. [11] and G. Q. Zhao et al. [12] studied the natural frequencies and modes of tires using numerical techniques and experimental modal analysis. However, researchers have rarely focused their attention on the vibration characteristics of non-pneumatic tire. Paul F. Joseph et al. [13-14] investigated the radial stiffness of non-pneumatic tire using the curved beam model, and analyzed the influence of flexible ring material properties and the number of spokes on the radial stiffness. Doo-Man Kim et al. [15] investigated the vibration characteristics of non-pneumatic tire under various loading and rotating situations using the finite element method.

The mechanical properties, trafficability and traction ability of MEW were systematic investigated by our research group, the results indicated that the proposed wheel had good trafficability, traction ability and small rolling resistance [16-18]. To further investigate the vibration characteristics of the MEW, the MEW is modeled as a REF-model with distributed spring stiffness and damping coefficients in the radial and tangential directions in this paper. The general forced solutions of inextensible vibration are derived by the use of a modal expansion technique and Arnoldi method under free suspension and various loading situations. And the present analytical solution is compared with the results of finite element simulation and modal test. Moreover, the effects of various rotational speed, loading and different number of hinges on natural frequencies are presented.

2. Structure of mechanical elastic wheel

MEW is mainly consisted of elastic wheel (rubber tread, elastic bead ring, clasps), hub, pins, hinges and other accessories as shown in Fig. 1.

Fig. 1Structure of MEW

In the process of the operated MEW, the vertical loading and torque are transferred from the axle to hub, and it is transmitted to the elastic wheel through hinges. The state of hinges is changed from equilibrium to preload. The obvious flexure deformation appeared in the tangent parts between elastic wheel and the ground under the action of vertical loading. The deformed tendency of the radial shrink appears in the upper parts with respect to the free conditions, it is similar to the elliptic type. Then the generating pulling force overcomes the static friction force of the wheel in ground contact to impelling the wheel forward. Based on the design of wheel structure, the hinges only load on pulling force rather than pressure. The hub is suspended in the elastic wheel depending on the tensile strength of hinges. The hinges that closer to the ground are slightly curved, the hub slips a distance relative to the free condition downward under the action of vertical loading. The elastic wheel endures most of the excitation from the road when the wheel is moving, displays elastic deformation instantaneously, and instantaneous bending of hinges has been relieved accordingly. Meanwhile, it is shown that the transmitting type of the MEW structure has become not only high efficiency and excellent adhesion ability, but also favorable buffer damping and passing ability. Therefore, the buffer damping performance of MEW is different from the ordinary pneumatic tire [16].

3. Theoretical modeling

3.1. Model description

Fig. 2 shows an unconstrained ring rotating at a constant angular speed $Ω$ about the $Y$ -axis where the $X Y Z$ -coordinate system is a space-fixed inertial frame. The MEW is modeled by a rotating ring on an elastic foundation with distributed spring stiffness $k_{r}$ and $k_{v}$ in the radial and tangential directions. The location of an element of MEW is described using the cylindrical coordinates ( $R$ , $ϕ$ ) in the non-rotating coordinate system, or ( $R$ , $θ$ ) in the rotating coordinate system, where $R$ is the radius of undeformed centroidal line of the ring. The ring thickness, width, applied radial force, tangential force and moment are denoted by $h$ , $b$ , $q_{r}$ , $q_{v}$ and $q_{β}$ , respectively. Only planar deformations of the ring are considered.

Fig. 2Schematic of MEW ring model and elastic wheel section structure

3.2. The equations of motion

It is assumed that the transverse shear deformation and the warping of the cross-section due to torsion negligibly small. According to the Lagrange strain theory, the flexural strain of the ring can be written as:

1

ε_{θ} = \frac{\bar{v'} + \bar{w}}{R} + \frac{1}{2} {(\frac{\bar{v'} + \bar{w}}{R})}^{2} + \frac{1}{2} {(\frac{\bar{v} - \bar{w'}}{R})}^{2},

where $\bar{v}$ and $\bar{w}$ are circumferential and radial displacements of the ring, respectively, and the prime denotes the partial derivative with respect to $θ$ . Note that the nonlinear term has to be included in the strains in order to compute the initial stress work correctly [8].

Inserting the expressions $\bar{v} = v + z β$ , $\bar{w} = w$ , and $β = (v - w' / R)$ into Eq. (1) yields:

2

ε_{θ} = \frac{v' + w}{R} + \frac{1}{2} {(\frac{v' + w}{R})}^{2} + \frac{1}{2} {(\frac{v - w'}{R})}^{2} + \frac{z}{R^{2}} (v^{'} - w^{''}),

where $v$ and $w$ represent the mid-plane displacements.

Based on the Hamilton principle, all the energies can be expressed as follows:

3

U = b \int_{0}^{2 π} \int_{- (h / 2)}^{(h / 2)} [\frac{1}{2} σ_{θ} ε_{θ} + σ_{θ}^{0} ε_{θ}] R d z d θ + \int_{0}^{2 π} [\frac{1}{2} k_{r} w^{2} + \frac{1}{2} k_{v} v^{2}] R d θ,

4

T = \int_{0}^{2 π} \frac{1}{2} ρ A R [{(\dot{w} - Ω v)}^{2} + {(\dot{v} + Ω (w + R))}^{2}] d θ,

5

W = \int_{0}^{2 π} (q_{r} w + q_{v} v + q_{β} \frac{v' - w}{R}) R d θ,

where $U$ , $T$ , and $W$ are potential, kinetic, and external force energy, respectively, $σ_{θ}^{0}$ , $A$ , and $ρ$ represent the initial stress, the ring section area and density.

Hamilton principle can be expressed as:

6

δ \int_{t_{1}}^{t_{2}} [U - T - W] d t = 0,

where $δ$ is the vibrational symbol.

Substituting Eqs. (2)-(5) into Eq. (6), and rearranging the terms yields the final equations of motion:

7

\begin{array}{l} - \frac{E I}{R^{4}} (v^{(6)} + 2 v^{(4)} + v^{''}) + σ_{θ}^{0} A \frac{1}{R^{2}} (v^{(4)} + 2 v^{''} + v) \end{array},

+ ρ A (\ddot{v} - {\ddot{v}}^{''} - 4 Ω {\dot{v}}^{'} + Ω^{2} (v^{''} - v)) - k_{r} v^{''} + k_{v} v = q_{r}^{'} + q_{v} + \frac{1}{R} q_{β}^{''} + \frac{1}{R} q_{β},

where the inextensible condition, i.e., $v' + w = 0$ , is used. $E$ is effective modulus of the belt and $I$ the moment of inertia in Eq. (7), which means $E I$ is the effective bending stiffness of the ring.

The natural frequencies of the non-rotating wheel are used to determine the model parameters and validate the model. With $Ω$ , $q_{r}$ , $q_{v}$ and $q_{β}$ being set to zero, Eq. (7) can be written as:

8

- \frac{E I}{R^{4}} (v^{(6)} + 2 v^{(4)} + v^{''}) - k_{r} v^{''} + k_{v} v + ρ A (\ddot{v} - {\ddot{v}}^{''}) = 0 .

Assume the free vibration mode in the sinusoidal series as follows:

9

v (θ, t) = \sum_{n = 0}^{\infty} A_{n} \sin (n θ + ω_{n} t) .

Inserting Eq. (9) into Eq. (8) yields the following natural frequency expression:

10

f_{n} = \frac{ω_{n}}{2 π} = \frac{1}{2 π} {\{\frac{1}{ρ A (1 + n^{2})} [\frac{E I}{R^{4}} (n^{6} - 2 n^{4} + n^{2}) + k_{r} n^{2} + k_{v}]\}}^{\frac{1}{2}} .

The undamped vibration frequencies of the rotating MEW can be determined by setting $q_{r}$ , $q_{v}$ and $q_{β}$ equal to zero and assuming a solution of the form:

11

v (ϕ, t) = B_{0} s i n ω_{0} t + (B_{1} s i n ϕ + {\bar{B}}_{1} c o s ϕ) s i n ω_{1} t + \sum_{n = 2}^{\infty} B_{n} \sin (n ϕ + ω_{n} t),

where:

12

ϕ = θ + Ω t .

The assumed form of the solution used in the summation of Eq. (11) reflects the fact that the mode shapes for $n >$ 1 are complex.

For $n >$ 1, substituting Eq. (11) into Eq. (7) yields the following vibration frequency expression:

13

{\bar{f}}_{n} = |\frac{2 Ω n}{n^{2} + 1} \pm \sqrt{\frac{Ω^{2} n^{2} (n^{2} - 1)^{2}}{(n^{2} + 1)^{2}} + {f_{n}}^{2}}| .

The vibration frequencies of the rotating MEW in the reference frame given by Eq. (13) are consistent with previous results in which the equations of motion are formulated in a rotating rather than a fixed reference frame.

In order to obtain the dynamic response of MEW in ground contact with a flat frictionless surface, the MEW is modeled as a ring on stiffness foundation with viscous damping coefficients $c_{r}$ and $c_{v}$ in the radial and tangential directions. The equation of motion Eq. (7) can be modified as:

14

\begin{array}{l} - \frac{E I}{R^{4}} (v^{(6)} + 2 v^{(4)} + v^{''}) + \frac{σ_{θ}^{0} A}{R^{2}} (v^{(4)} + 2 v^{''} + v) + ρ A (\ddot{v} - {\ddot{v}}^{''} - 4 Ω {\dot{v}}^{'} + Ω^{2} (v^{''} - v)) \end{array}

- k_{r} v^{''} + k_{v} v - c_{r} {\dot{v}}^{''} + c_{v} \dot{v} = q_{r}^{'} + q_{v} + \frac{q_{β}^{''}}{R} + \frac{q_{β}}{R} .

As an approximation, the functions $v (ϕ, t)$ and ${q'}_{r} + q_{v}$ are expanded as Fourier series:

15

v (ϕ, t) = \frac{α_{0}}{2} + \sum_{n = 0}^{N} [α_{n} (t) c o s n ϕ + β_{n} (t) s i n n ϕ],

16

{q'}_{r} + q_{v} = \frac{g_{0}}{2} + \sum_{n = 0}^{N} [g_{n} (t) c o s n ϕ + h_{n} (t) s i n n ϕ],

where:

17

g_{n} = \frac{1}{π} \int_{0}^{2 π} (q_{v} c o s n ϕ + n q_{r} s i n n ϕ) d ϕ,

18

h_{n} = \frac{1}{π} \int_{0}^{2 π} (q_{v} s i n n ϕ - n q_{r} c o s n ϕ) d ϕ .

Substituting Eq. (15), (16) and $q_{β} = 0$ into Eq. (14) yields a set of linear second-order ordinary differential equations:

19

n = 0 : ρ A {\ddot{α}}_{0} + c_{v} {\dot{α}}_{0} + k_{v} α_{0} = g_{0},

20

n \neq 0 : [\begin{array}{l} m_{n} \\ m_{n} \end{array}] \{\begin{array}{l} {\ddot{α}}_{n} \\ {\ddot{β}}_{n} \end{array}\} + [\begin{matrix} c_{n} & γ_{n} \\ - γ_{n} & c_{n} \end{matrix}] \{\begin{array}{l} {\dot{α}}_{n} \\ {\dot{β}}_{n} \end{array}\} + [\begin{matrix} k_{n} & τ_{n} \\ - τ_{n} & k_{n} \end{matrix}] \{\begin{array}{l} α_{n} \\ β_{n} \end{array}\} = \{\begin{array}{l} g_{n} \\ h_{n} \end{array}\},

where:

21

m_{n} = ρ A (1 + n^{2}),

22

c_{n} = n^{2} c_{r} + c_{v},

23

γ_{n} = - 4 ρ A n Ω,

24

k_{n} = \frac{E I n^{2} (n^{2} - 1)^{2}}{R^{4}} + k_{r} n^{2} + k_{v} - ρ A (n^{2} + 1) Ω^{2},

25

τ_{n} = n Ω c_{n} .

The assumption of a frictionless surface implies that the contact forces acting on the MEW are in the vertical direction alone. Thus, the radial and tangential forces are given by:

26

q_{r} (ϕ) = \bar{F} (ϕ) s i n ϕ,

27

q_{v} (ϕ) = \bar{F} (ϕ) c o s ϕ,

where $\bar{F} (ϕ)$ is the contact force.

Consistent with the use of Eq. (15) one can verify that the ring model of the MEW only contacts the flat surface at a finite number of angular locations. Thus, the contact force can be expressed as:

28

\bar{F} (ϕ) = \sum_{i = 1}^{m} {\bar{F}}_{c i} η (ϕ - ϕ_{c i}),

where $m$ is the number of contact locations, $η$ denotes the delta function, and $ϕ_{c i}$ is the $i$ th angular location where contact occurs. The number of contact locations $m$ depends on the ring and foundation stiffness [9].

For the problems considered in this study, contact occurs at four locations. Force equilibrium in the vertical direction and contact at the angular locations $ϕ_{c i}$ imply that:

29

\sum_{i = 1}^{m} {\bar{F}}_{c i} = \frac{F}{R b},

30

r (ϕ_{c i}) = r (ϕ_{c, i + 1}), (i = 1, \dots, m - 1),

where $r (ϕ)$ denotes the vertical position of the ring at the angular location $ϕ$ due to concentrated forces at the angular locations $ϕ_{c i}$ . Using the inextensibility assumption $w = - v'$ , one obtains:

31

r (ϕ) = (R - v^{'} (ϕ)) s i n ϕ + v (ϕ) c o s ϕ .

And at the point of contact:

32

r^{'} (ϕ_{c i}) = 0, (i = 1, \dots, m) .

Making use of Eqs. (17, 18) and Eqs. (26)-(29) in the absence of vibration yields:

33

α_{0} = \frac{g_{0}}{k_{v}},

34

α_{n} = \frac{(k_{n} g_{n} - n Ω c_{n} h_{n})}{(k_{n}^{2} + n^{2} Ω^{2} c_{n}^{2})},

35

β_{n} = \frac{(k_{n} h_{n} + n Ω c_{n} g_{n})}{(k_{n}^{2} + n^{2} Ω^{2} c_{n}^{2})},

where:

36

g_{n} = \sum_{i = 1}^{m} \frac{{\bar{F}}_{c i} (c o s n ϕ_{c i} c o s ϕ_{c i} + n s i n n ϕ_{c i} s i n ϕ_{c i})}{π},

37

h_{n} = \sum_{i = 1}^{m} \frac{{\bar{F}}_{c i} (s i n n ϕ_{c i} c o s ϕ_{c i} - n c o s n ϕ_{c i} s i n ϕ_{c i})}{π} .

Substituting Eq. (15) into Eq. (31), one obtains:

38

r (ϕ) = R s i n ϕ + \frac{α_{0}}{2} c o s ϕ

+ \sum_{n = 1}^{N} [α_{n} (n s i n n ϕ s i n ϕ + c o s n ϕ c o s ϕ) - β_{n} (n c o s n ϕ s i n ϕ - s i n n ϕ c o s ϕ)] .

Differentiating Eq. (38) with respect to $ϕ$ yields:

39

r^{'} (ϕ) = R c o s ϕ - \frac{α_{0}}{2} s i n ϕ + \sum_{n = 1}^{N} (n^{2} - 1) (α_{n} c o s n ϕ s i n ϕ + β_{n} s i n n ϕ s i n ϕ),

40

r'' (ϕ) = - R s i n ϕ - \frac{α_{0}}{2} c o s ϕ

+ \sum_{n = 1}^{N} (n^{2} - 1) [α_{n} (c o s n ϕ c o s ϕ - n s i n n ϕ s i n ϕ) + β_{n} (s i n n ϕ c o s ϕ + n c o s n ϕ s i n ϕ)] .

According to the Newton’s method, the solution of all cases can be obtained. The Eq. (33)-(39) and constraint Eq. (30) and (32) can be expressed solely in terms of the unknown ${\bar{F}}_{c i}$ and $ϕ_{c i}$ . The exception is for cases in which there is only one contact location and the MEW is stationary ( $Ω =$ 0) and undamped. The solution is given simply by ${\bar{F}}_{c 1} = F / (R b)$ and $ϕ_{c 1} = 0$ .

Assume the perturbation variables ${\bar{α}}_{n}$ , ${\bar{β}}_{n}$ , $\bar{F}$ and ${\bar{ϕ}}_{c i}$ of the static equilibrium as follows:

41

{\bar{α}}_{n} = α_{n} - {\tilde{α}}_{n},

42

{\bar{β}}_{n} = β_{n} - {\tilde{β}}_{n},

43

{\overset{̿}{F}}_{c i} = {\bar{F}}_{c i} - {\tilde{\bar{F}}}_{c i}

44

{\bar{ϕ}}_{c i} = ϕ_{c i} - {\tilde{ϕ}}_{c i} .

Linearizing the Eq. (19) and Eq. (20) about the static equilibrium yields:

45

n = 0 : \frac{ρ A {\ddot{\bar{α}}}_{0} + c_{v} {\dot{\bar{α}}}_{0} + k_{v} {\bar{α}}_{0}}{2} = \frac{(d_{01}^{T} {\overset{̿}{F}}_{c} + d_{03}^{T} {\tilde{f}}_{c} {\bar{ϕ}}_{c})}{2 π},

46

n \neq 0 : [\begin{array}{l} m_{n} \\ m_{n} \end{array}] \{\begin{array}{l} {\ddot{\bar{α}}}_{n} \\ {\ddot{\bar{β}}}_{n} \end{array}\} + [\begin{matrix} c_{n} & γ_{n} \\ - γ_{n} & c_{n} \end{matrix}] \{\begin{matrix} {\dot{\bar{α}}}_{n} \\ {\dot{\bar{β}}}_{n} \end{matrix}\} + \{\begin{matrix} k_{n} & τ_{n} \\ - τ_{n} & k_{n} \end{matrix}\} \{\begin{matrix} {\bar{α}}_{n} \\ {\bar{β}}_{n} \end{matrix}\} = \frac{1}{π} [\begin{matrix} d_{n 1}^{T} {\overset{̿}{F}}_{c} + d_{n 3}^{T} {\tilde{f}}_{c} {\bar{ϕ}}_{c} \\ d_{n 2}^{T} {\overset{̿}{F}}_{c} + d_{n 4}^{T} {\tilde{f}}_{c} {\bar{ϕ}}_{c} \end{matrix}],

where:

47

{\overset{̿}{F}}_{c i} = {[{\overset{̿}{F}}_{c i}, \dots, {\overset{̿}{F}}_{c m}]}^{T},

48

{\bar{ϕ}}_{c} = [{\bar{ϕ}}_{c 1}, \dots, {\bar{ϕ}}_{c m}]^{T},

49

{\tilde{f}}_{c} = d i a g [{\tilde{\bar{F}}}_{c i}, \dots, {\tilde{\bar{F}}}_{c m}],

50

d_{n 1} (i) = c o s n {\tilde{ϕ}}_{c i} c o s {\tilde{ϕ}}_{c i} + n s i n n {\tilde{ϕ}}_{c i} s i n {\tilde{ϕ}}_{c i},

51

d_{n 2} (i) = s i n n {\tilde{ϕ}}_{c i} c o s {\tilde{ϕ}}_{c i} - n c o s n {\tilde{ϕ}}_{c i} s i n {\tilde{ϕ}}_{c i},

52

d_{n 3} (i) = (n^{2} - 1) c o s n {\tilde{ϕ}}_{c i} s i n {\tilde{ϕ}}_{c i},

53

d_{n 4} (i) = (n^{2} - 1) s i n n {\tilde{ϕ}}_{c i} s i n {\tilde{ϕ}}_{c i}

Linearizing the contact boundary conditions $r (ϕ_{c i}) = r (ϕ_{c, i + 1})$ and $r' (ϕ_{c i}) = 0$ about the static equilibrium yields:

54

H_{1}^{T} Q = 0,

55

Φ {\bar{ϕ}}_{c} + H_{2}^{T} Q = 0,

where:

56

H_{1} = [d_{01} / 2 d_{11} d_{12} \dots d_{M 1} d_{M 2}]^{T},

57

H_{2} = [d_{03} / 2 d_{13} d_{14} \dots d_{M 3} d_{M 4}]^{T},

58

Q = [{\bar{α}}_{0} {\bar{α}}_{1} {\bar{β}}_{1} \dots {\bar{α}}_{N} {\bar{β}}_{N}]^{T},

59

Φ = d i a g [r^{''} ({\tilde{ϕ}}_{c 1}), \dots, r^{''} ({\tilde{ϕ}}_{c m})] .

Eq. (45), (46) and Eq. (54)-(55) can be expressed equivalently as:

60

[\begin{matrix} {\hat{m}}_{n} & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{matrix}] [\begin{matrix} \ddot{Q} \\ {\ddot{\bar{\bar{F}}}}_{c} \\ {\ddot{\bar{ϕ}}}_{c} \end{matrix}] + [\begin{array}{l} \hat{λ} & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{array}] [\begin{matrix} \dot{Q} \\ {\dot{\bar{\bar{F}}}}_{c} \\ {\dot{\bar{ϕ}}}_{c} \end{matrix}] + [\begin{matrix} \hat{ψ} & - H_{1} / π & - H_{2} {\tilde{f}}_{c} / π \\ H_{1}^{T} & 0 & 0 \\ H_{2}^{T} & 0 & Φ \end{matrix}] [\begin{matrix} Q \\ {\overset{̿}{F}}_{c} \\ {\bar{ϕ}}_{c} \end{matrix}] = [\begin{matrix} 0 \\ 0 \\ 0 \end{matrix}],

where:

61

{\hat{m}}_{n} (n, n) = ρ A (n^{2} + 1),

62

{\hat{m}}_{1} (1, 1) = \frac{ρ A}{2}, \hat{λ} (1, 1) = \frac{c_{v}}{2}, \hat{ψ} (1, 1) = \frac{k_{v}}{2},

63

\hat{λ} (2 n, 2 n) = c_{n}, \hat{λ} (2 n, 2 n + 1) = γ_{n},

\hat{λ} (2 n + 1, 2 n) = - γ_{n}, \hat{λ} (2 n + 1, 2 n + 1) = c_{n},

64

\hat{ψ} (2 n, 2 n) = k_{n}, \hat{ψ} (2 n, 2 n + 1) = τ_{n},

\hat{ψ} (2 n + 1, 2 n) = - τ_{n}, \hat{ψ} (2 n + 1, 2 n + 1) = k_{n} .

The mode shapes and frequencies of the system in ground contact can be determined by transforming Eq. (60) to solving the associated eigenvalue problem.

3.3. Model parameters of MEW

The geometrical and structural model parameters are obtained directly from the MEW design, as shown in Table 1.

Table 1Model geometrical and structural parameters of MEW

Parameters type	Value
$ρ$	2.56×10³kg/m³
$h$	0.05 m
$b$	0.317 m
$A$	0.01585 m²
$R$	0.425 m

Through analyzing the finite element simulation and experimental mode shape corresponding to each natural frequency, the physical parameters of theoretical model, including $k_{r}$ , $k_{v}$ and $E I$ are determined.

4. Modal analysis and model validation

To validate the accuracy of the analytical model of the MEW, the finite element modal analysis and experimental modal analysis of the wheel are presented, respectively. And some physical parameters of the analytical model are identified by the results of the modal comparative analysis.

4.1. Finite element modal analysis

The MEW model is appropriately simplified without considering the influence of tread pattern and the irregular interface shape of hub on the calculation precision. The hinges can be simplified into 12 groups of three bar linkage. The three-dimensional geometric model of MEW is imported into ABAQUS, and the model is divided by C3D10M elements. The elastic bead ring is also modeled by solid elements, the cord ply comprises cord-rubber composite which displays the mechanics anisotropy and non-linearity, and it is modeled using the rebar layer. The generated finite element model is shown in Fig. 3. In the simulation, the stiffness calculations for the rebar elements use the same integration nodes as the calculation for the underlying rubber shell elements.

Fig. 3FEM model of the MEW

According to the results of finite element simulation, the in-plane radial modal parameters of MEW in first six order are obtained under the number of circumferential hinges is 12 sets, the simulation results are shown in Fig. 4.

4.2. Experimental modal analysis

The test wheel is a mechanical elastic wheel, it is mounted on the Suspended platform and the multifunctional testbed separately. And the test and analysis apparatus are listed in Table 2.

Table 2Apparatus for test

Apparatus	Type
Multichannel data acquisition and processing system	SCM205
Modal test analysis software	LMS Test.Lab 13A
Force hammer	086C03
Acceleration sensor	333B30
Computer workstation	Dell

Modal test system for the free suspension MEW, the MEW is under the free suspension state that is shown in Fig. 5(a). According to the hammering method, the radial vibration excitation test is executed. The wheel response in the radial direction is measured with the acceleration sensor, and at some locations, the response in all three translational directions is found. Set the frequency band of 512 Hz and the frequency resolution of 0.83 Hz, and all data are logged on an LMS Test.Lab system, connects to a computer workstation. Since the natural frequency of the suspension system is very low (less than 1 Hz), and is far below the tire first natural frequency (about 30 Hz), the influence of the suspension system on modal parameters can be ignored.

Fig. 4Simulation natural frequencies and mode shapes of the MEW

a) 30.31 Hz (1st mode)

b) 90.83 Hz (2nd mode)

c) 107.46 Hz (3rd mode)

d) 134.27 Hz (4th mode)

e) 158.64 Hz (5th mode)

f) 207.36 Hz (6th mode)

Fig. 5Modal test of the MEW

a) Suspension state

b) Loaded state

According to the parameter identification method of frequency response function matrix and PolyMAX, the vibration mode and natural frequencies of MEW are obtained under the number of circumferential hinges is 12 sets, as shown in Fig. 6.

4.3. Parameters identification and model validation

Through the contrastive analysis, the simulation modal parameters and experimental modal parameters are in good agreement. The error of mode shape is probably due to the damping of the elastic wheel and the influence of the environment. To obtain more accurate physical parameters of the analytical model, the physical parameters of Eq. (10) are fitted by the nonlinear least squares optimization method based on the natural frequencies of the wheel modal test. Combining with the natural frequencies of $n =$ 1, 4, 6, then solving for the model physical parameters lead to $k_{r} =$ 2.229×10⁷ N/m², $k_{v} =$ –1.898×10⁷ N/m², and $E I =$ 1218.22 N·m². Substituting the obtained physical parameters into Eq. (10) yields all natural frequencies, and which are substituted into Eq. (9) to obtain the corresponding mode shapes, the calculation results are shown in Fig. 7.

Fig. 6Experimental natural frequencies and mode shapes of the MEW

a) 32.14 Hz (1st mode)

b) 88.12 Hz (2nd mode)

c) 105.47 Hz (3rd mode)

d) 131.73 Hz (4th mode)

e) 153.65 Hz (5th mode)

f) 202.59 Hz (6th mode)

Fig. 7Theoretical natural frequencies and mode shapes of the MEW

a) 32.26 Hz (1st mode)

b) 94.24 Hz (2nd mode)

c) 112.16 Hz (3rd mode)

d) 132.05 Hz (4th mode)

e) 160.81 Hz (5th mode)

f) 203.13 Hz (6th mode)

The finite element simulation, modal test results and analytical solution are in excellent agreement for the first six order natural frequency and mode shape in Fig. 4, Fig. 6 and Fig. 7. The maximum error between finite element method results and analytical solution is 4.19 %, and the maximum error between modal test results and analytical solution is only 6.49 %. The good agreement among the predictions, simulation and experiments, to a certain extent, validates the analytical model, and it indicates that the REF model of wheel will be used as a theoretical basis to study the mechanical dynamic characteristics.

5. Influencing factors analysis

Through the vibration modal analysis of the free situations and different working conditions, it can be seen that the influence of changes of the working conditions on the mode shapes of MEW are less. Thus, according to the single variable method, the influence of different working conditions on the vibration characteristics of MEW is studied by the change of natural frequency.

5.1. Influence of loading on vibration

The vibration characteristics of MEW under different vertical loads are analyzed by the method of experimental modal and theoretical prediction. Fig. 5(b) gives the loaded state of wheel, and the vertical load during the test is real-timely monitored by means of a static strain sensor. The vibration frequencies of the MEW in ground contact under a certain load ( $F =$ 3000 N) are shown in Table 3.

Table 3Natural frequencies of the loaded MEW (F= 3000 N): experimental vs. REF model

Rank	$n =$ 1	$n =$ 2	$n =$ 3	$n =$ 4	$n =$ 5	$n =$ 6
$f_{n}$ : REF model	35.39	101.13	119.87	142.08	171.82	219.86
$f_{n}$ : Experimental	33.26	96.68	116.59	136.92	164.27	211.45

The theoretical prediction value and experimental results are in good agreement for the natural frequencies of the loaded MEW. According to the damping ratios of the experimental $n =$ 1 and $n =$ 2 modes, and making using of equations $c_{v} = 2 ζ_{0} \sqrt{k_{v} ρ h}$ and $c_{r} = 2 ζ_{1} \sqrt{2 (k_{r} + k_{v}) ρ h} - c_{v}$ yields the damping coefficient $c_{r}$ and $c_{v}$ .

Effect of vertical load on natural frequencies of the MEW radial vibration under the number of circumferential hinges is 12 sets as shown in Fig. 8.

The above analysis shows that the influence of changes of the vertical load on the MEW natural frequency is less. However, the modal frequency increases with the increasing of the vertical load. When the vertical load is greater, the corresponding vibration frequency is also greater when compare with the natural frequencies of the free suspension MEW. When the vertical load is larger than 4000 N, the frequency of each mode has a mutation. Because the contact area of the MEW produces a partial ‘warp’ phenomenon when the load exceeds a certain value.

Fig. 8Effect of vertical load on vibration frequency of radial mode of MEW

Fig. 9Effect of rotating on vibration frequency of the damped system in ground contact

5.2. Influence of rotating on vibration

The vibration modal of MEW under different rotating speed are analyzed by the finite element simulation and the linearized equations of motion Eq. (60) under a certain load, respectively. The analysis indicates that the finite element method results and analytical solution are in excellent agreement. It can be seen that the influence of rotating on the mode shapes of MEW are less, and Fig. 9 illustrates the influence of rotating on natural frequencies in the condition that the number of hinges remain unchanged.

The analysis results show the radial natural frequency decreases with the increasing of rotational speed in the damped system. The exception is for the first-order form in which the rotational speed has no effect on the first-order frequency. However, the frequencies of modes (2nd)-(6th) all have a noticeable dependence on $Ω$ , and the low speed has little influence on the natural frequency of MEW.

5.3. Influence of the hinges changing on vibration

According to the different matching vehicle models, the number of MEW’s hinges can be adjusted appropriately to cope with the requirements for ride comfort performance. The MEW style can be divided into 6, 9, 12 and 18 sets of hinges. Hence the influence of the number of hinges on the dynamic characteristics of the MEW is studied. Fig. 10 illustrates the influence of different number of the hinges on natural frequency of radial vibration under the free suspension state. The radial natural frequency of each mode has become increased with number of the hinges. Effect of different number of the hinges on frequencies of radial vibration under a certain load ( $F =$ 5100 N) as shown in Fig. 11.

Fig. 10Effect of number of the hinges on natural frequency under the free suspension state

Fig. 11Effect of number of the hinges on natural frequency under a certain load

Compare with the free suspension MEW, more number of the hinges, higher the corresponding natural frequency. The natural frequency of higher mode shape is sensitive to the changing number of the hinges.

Table 4 shows the damping ratio of the MEW radial modes of free suspension under different number of hinges. The damping ratio decreases with the increasing of the number of hinges, which makes the MEW vibration amplitudes of the corresponding natural frequencies increase, especially for high frequencies. Therefore, the increase in the number of hinges helps reduce the rolling resistance of MEW, and improves the life of MEW. However, for the same vehicle model, the ride comfort performance is reduced.

Table 4Damping ratio of MEW radial modes under different number of hinges

Damping ratio (%)		Number of hinges
Damping ratio (%)		6	9	12	18
Mode order	1	6.85	6.08	4.91	3.92
	2	5.94	5.22	4.36	3.43
	3	5.13	4.37	3.59	2.71
	4	4.56	3.83	3.27	2.39
	5	4.24	3.51	2.82	2.12
	6	4.05	3.33	2.49	1.96

6. Conclusions

1) This paper develops combined numerical analysis and test method to deal with vibration characteristics and its key influencing factors of a new mechanical elastic wheel. The MEW is modeled as a ring on elastic foundations with distributed spring stiffness in the radial and tangential directions. The vibration frequencies and mode shapes of the MEW can be analyzed by using the matrix perturbation theory and the Arnoldi method under free suspension and various loading situations. And the results of the analytical solution are validated through finite element simulation and modal test.

2) The analytical solutions are compared with the finite element simulation and the modal test under the condition of free suspension and various loading situations, it can be seen that the developed theory produces quite accurate results. With such an approach, the influence of different working conditions on vibration frequencies can be analyzed. It is shown that the modal frequency increases with the increasing of vertical loading, the corresponding vibration frequency is also greater when compare with the natural frequencies of the free suspension MEW. And the radial vibration frequency decreases with the increase of rotational speed in the damped system. Compare with the free suspension MEW, more number of the hinges, higher the corresponding natural frequency. Moreover, the damping ratio decreases with the increasing of the number of hinges, which makes the MEW vibration amplitudes of the corresponding natural frequencies increase. The analysis results reflect the objective law of the actual vibration characteristics of the MEW, and provide a reference for the MEW structure optimization and the vibration characteristics of the whole vehicle.

References

Kozhevnikov I. F. The vibrations of a free and loaded tyre. Journal of Applied Mathematics and Mechanics, Vol. 70, Issue 2, 2006, p. 250-256.

Publisher
Kozhevnikov I. F. Vibrations of a rolling tyre. Journal of Sound and Vibration, Vol. 331, Issue 7, 2012, p. 1669-1685.

Publisher
Huang S. C., Soedel W. Effects of coriolis acceleration on the free and forced in-plane vibrations of rotating rings on elastic foundation. Journal of Sound and Vibration, Vol. 115, Issue 2, 1987, p. 253-274.

Publisher
Huang S. C., Su C. K. In-plane dynamics of tires on the road based on an experimentally verified rolling ring model. Vehicle System Dynamics, Vol. 21, Issue 1, 1992, p. 247-267.

Publisher
Bickford W. B., Reddy E. S. On the in-plane vibrations of rotating rings. Journal of Sound and Vibration, Vol. 101, Issue 1, 1985, p. 13-22.

Publisher
Kim W., Chung J. Free non-linear vibration of a rotating thin ring with the in-plane and out-of-plane motions. Journal of Sound and Vibration, Vol. 258, Issue 1, 2002, p. 167-178.

Publisher
Kim S. J., Savkoor A. R. Contact problem of in-plane rolling of tires on a flat road. Vehicle System Dynamics, Vol. 27, Issue 27, 1997, p. 189-206.

Publisher
Wei Y. T., Nasdala L., Rothert H. Analysis of forced transient response for rotating tires using REF models. Journal of Sound and Vibration, Vol. 320, Issue 1, 2009, p. 145-162.

Publisher
Dohrmann C. R. Dynamics of a tire-wheel-suspension assembly. Journal of Sound and Vibration, Vol. 210, Issue 5, 1998, p. 627-642.

Publisher
Huang S. C., Hsu B. S. An approach to the dynamic analysis of rotating tire-wheel-suspension units. Journal of Sound and Vibration, Vol. 156, Issue 3, 1992, p. 505-519.

Publisher
Yam L. H., Guan D. H., Zhang A. Q. Three-dimensional mode shapes of a tire using experimental modal analysis. Experimental Mechanics, Vol. 40, Issue 4, 2000, p. 369-375.

Publisher
Guan Yanjin, Cheng Gang, Zhao Guoqun, et al. Investigation of the vibration characteristics of radial tires using experimental and numerical techniques. Journal of Reinforced Plastics and Composites, Vol. 30, Issue 24, 2011, p. 2035-2050.

Publisher
Gasmi Amir, Joseph Paul F., Rhyne Timothy B., et al. Closed-form solution of a shear deformable extensional ring in contact between two rigid surfaces. International Journal of Solids and Structures, Vol. 48, Issue 5, 2011, p. 843-853.

Publisher
Gasmi Amir, Joseph Paul F., Rhyne Timothy B., et al. Development of a two-dimensional model of a compliant non-pneumatic tire. International Journal of Solids and Structures, Vol. 49, Issue 13, 2012, p. 1723-1740.

Publisher
Lee Chihun, Ju Jaehyung, Kim Doo-Man Vibration analysis of non-pneumatic tires with hexagonal lattice spokes. Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Chicago, USA, 2012.

Publisher
Wang Wei, Zhao Youqun, Wang Jian, et al. Structure analysis and ride comfort of vehicle on new mechanical elastic tire. Proceedings of the FISITA World Automotive Congress, Springer-Verlag, Berlin, Heidelberg, Vol. 10, 2013, p. 199-209.

Publisher
Zang Liguo, Zhao Youqun, Li Bo, et al. Influence of mechanical elastic wheel configuration for traction ability. Journal of Harbin Engineering University, Vol. 35, Issue 11, 2014, p. 1-6, (in Chinese).

Search CrossRef
Zhao Youqun, Zang Liguo, Chen Yueqiao, et al. Non-pneumatic mechanical elastic wheel natural dynamic characteristics and influencing factors. Journal of Central South University, Vol. 22, Issue 5, 2015, p. 1707-1715.

Publisher
Kint P., Sas P., Desmet W. Development and validation of a three-dimensional ring-based structural tyre Model. Journal of Sound and Vibration, Vol. 326, Issues 3-5, 2009, p. 852-869.

Publisher
Kint P., Sas P., Desmet W. Measurement and analysis of rolling tire vibrations. Optics and Lasers in Engineering, Vol. 47, Issues 3-4, 2009, p. 443-453.

Publisher
Sugiyama Hiroyuki, Suda Yoshihiro Modeling nonlinear flexible tire belt in the study of in-plane tire dynamics. Proceeding of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Nevada, USA, 2007.

Publisher
Kim Byoung Sam, Chi Chang Heon, Lee Tae Keun A study on radial directional natural frequency and damping ratio in a vehicle tire. Applied Acoustics, Vol. 65, Issue 5, 2007, p. 538-556.

Publisher

Cited by

Vibration analysis of radial tire using the 3D rotating hyperelastic composite REF based on ANCF

Bo Fan | Zhongmin Wang

(2024)

Nonlinear Dynamic Characteristics of Rotating Thick Ring on the Elastic Foundation Based on Plane Stress Theory

Bo Fan | Zhongmin Wang | Qingbo Wang

(2024)

Steady-state rolling resistance prediction model of non-pneumatic tyres considering tread temperature: theory and experiment

Shuo Guo | Youqun Zhao | Fen Lin | Tao Xu | Danyang Li

(2024)

Effect of Different Road Obstacles on the Structural Behavior of a Honeycomb Nonpneumatic Tire

Afshin Ashofteh | Amirhossein Shahdadi

(2023)

Nonlinear forced transient response of rotating ring on the elastic foundation by using adaptive ANCF curved beam element

Bo Fan | Zhongmin Wang | Qingbo Wang

(2022)

Radial stiffness analysis of circular ring with uncertain boundary conditions: Mechanical elastic wheel as an example

(2022)

Research on non-pneumatic tire with gradient anti-tetrachiral structures

Taoyu Wu | Mingxuan Li | Xiaolei Zhu | Xiaofeng Lu

(2020)

Multi-objective optimization for ride comfort of hydro-pneumatic suspension vehicles with mechanical elastic wheel

Youqun Zhao | Han Xu | Yaoji Deng | Qiuwei Wang

(2019)

Transient Dynamic Characteristics of a Non-Pneumatic Mechanical Elastic Wheel Rolling Over a Ditch

You-Qun Zhao | Yao-Ji Deng | Fen Lin | Ming-Min Zhu | Zhen Xiao

(2018)

Simulation of steady-state rolling non-pneumatic mechanical elastic wheel using finite element method

(2018)

Thermal modal analysis of novel non-pneumatic mechanical elastic wheel based on FEM and EMA

(2018)

About this article

Received

04 March 2016

Accepted

18 August 2016

Published

31 December 2016

SUBJECTS

Vibration in transportation engineering

DOI

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

Keywords

mechanical elastic wheel

vibration characteristics

REF model

finite element analysis

modal test

influencing factors analysis

Acknowledgements

This research work is partially supported by the Major Exploration Project of the General Armaments Department of China (No. NHA13002), the National Natural Science Foundation of China (No. 11072106), the Fundamental Research Funds for the Central Universities (No. NP 2016412) and Funding of Jiangsu Innovation Program for Graduate Education and the Fundamental Research Funds for the Central Universities (No. KYLX_0241).

Author Contributions

The research ideas of this paper and solution of the problem were proposed by Qiang Wang. Professor Youqun Zhao contributed to the guarantor of integrity of entire study, and gave a guide to the solution of the research problem and the writing style. Fen Lin provided some help in the writing style. Hongxun Fu helped perform the partial data processing. Xianbin Du provided some help in the process of finite element modeling.

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

2024 04 04

Modal and vibration response characteristics of rotating support platform

Qian Zhao

Research article

2023 08 01

Effect of excitation vibration on mechanical property and stress corrosion resistance of cast steel

Ruiying Shao, Hongjun Wang, Kunyan Lu, Juan Song

Research article

2016 09 30

Vibration performance analysis of vehicle with the non-pneumatic new mechanical elastic wheel in the impulse input experiment

Wang Wei, Zhu Kai, Bei Shaoyi, Zhang Lanchun, Wang Yongzhi

Research article

2016 09 30

Numerical analysis and optimization of wheel vibrations and radiation noises of the high-speed train

Xue-yu Lu, Zheng-qiao Lu, Li-zong Lin

Q. Wang, Y. Zhao, F. Lin, H. Fu, and X. Du, “Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel,” Journal of Vibroengineering, Vol. 18, No. 8, pp. 5337–5352, Dec. 2016, https://doi.org/10.21595/jve.2016.16952

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2016.16952
UR  - https://doi.org/10.21595/jve.2016.16952
TI  - Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel
T2  - Journal of Vibroengineering
AU  - Zhao, Youqun
AU  - Wang, Qiang
AU  - Lin, Fen
AU  - Fu, Hongxun
AU  - Du, Xianbin
PY  - 2016
DA  - 2016/12/31
PB  - JVE International Ltd.
SP  - 5337-5352
IS  - 8
VL  - 18
SN  - 1392-8716
ER  - 

Copy Ris

Copied to clipboard!

@article{Zhao_2016,
	doi = {10.21595/jve.2016.16952},
	url = {https://doi.org/10.21595/jve.2016.16952},
	year = 2016,
	month = {dec},
	publisher = {{JVE} International Ltd.},
	volume = {18},
	number = {8},
	pages = {5337--5352},
	author = {Youqun Zhao and Qiang Wang and Fen Lin and Hongxun Fu and Xianbin Du},
	title = {Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]Y. Zhao, Q. Wang, F. Lin, H. Fu, and X. Du, “Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel,” Journal of Vibroengineering, vol. 18, no. 8, pp. 5337–5352, Dec. 2016, doi: 10.21595/jve.2016.16952.

Copy IEEE

Copied to clipboard!

Zhao, Youqun, Qiang Wang, Fen Lin, Hongxun Fu, and Xianbin Du. “Research on Vibration Characteristics and Its Key Influencing Factors of New Mechanical Elastic Wheel.” Journal of Vibroengineering 18, no. 8 (December 31, 2016): 5337–52. https://doi.org/10.21595/jve.2016.16952.

Copy Chicago

Copied to clipboard!

Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel

Abstract

1. Introduction

2. Structure of mechanical elastic wheel

3. Theoretical modeling

3.1. Model description

3.2. The equations of motion

3.3. Model parameters of MEW

4. Modal analysis and model validation

4.1. Finite element modal analysis

4.2. Experimental modal analysis

4.3. Parameters identification and model validation

5. Influencing factors analysis

5.1. Influence of loading on vibration

5.2. Influence of rotating on vibration

5.3. Influence of the hinges changing on vibration

6. Conclusions

References

Cited by

About this article

Related Articles