Research on vibration characteristics and its key influencing factors of new mechanical elastic wheel
Qiang Wang^{1} , Youqun Zhao^{2} , Fen Lin^{3} , Hongxun Fu^{4} , Xianbin Du^{5}
^{1, 2, 3, 4, 5}College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 18, Issue 8, 2016, p. 53375352.
https://doi.org/10.21595/jve.2016.16952
Received 4 March 2016; received in revised form 14 July 2016; accepted 18 August 2016; published 31 December 2016
JVE Conferences
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.
Keywords: mechanical elastic wheel, vibration characteristics, REF model, finite element analysis, modal test, influencing factors analysis.
1. Introduction
To change the existing pneumatic tire puncture damage and high speed blowout etc., developing runflat and antipuncture 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 nonpneumatic tires with different structures. Based on that, the new mechanical elastic wheel (MEW) as a nonpneumatic 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 [12] investigated the vibration of a tyre using the model of a wheel with a reinforced tyre under various loading and various rotating situations. The nonlinear vibration of a rotating ring was investigated by W. B. Bickford [3] and W. Kim [4]. S. C. Huang [56] 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 tirewheelsuspension 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 nonpneumatic tire. Paul F. Joseph et al. [1314] investigated the radial stiffness of nonpneumatic tire using the curved beam model, and analyzed the influence of flexible ring material properties and the number of spokes on the radial stiffness. DooMan Kim et al. [15] investigated the vibration characteristics of nonpneumatic 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 [1618]. To further investigate the vibration characteristics of the MEW, the MEW is modeled as a REFmodel 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. 1. Structure 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 $\mathrm{\Omega}$ about the $Y$axis where the $XYZ$coordinate system is a spacefixed 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$, $\varphi $) in the nonrotating coordinate system, or ($R$, $\theta $) 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}_{\beta}$, respectively. Only planar deformations of the ring are considered.
Fig. 2. Schematic 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 crosssection due to torsion negligibly small. According to the Lagrange strain theory, the flexural strain of the ring can be written as:
where $\overline{v}$ and $\overline{w}$ are circumferential and radial displacements of the ring, respectively, and the prime denotes the partial derivative with respect to $\theta $. 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 $\overline{v}=v+z\beta $, $\overline{w}=w$, and $\beta =(vw\mathrm{\text{'}}/R)$ into Eq. (1) yields:
where $v$ and $w$ represent the midplane displacements.
Based on the Hamilton principle, all the energies can be expressed as follows:
where $U$, $T$, and $W$ are potential, kinetic, and external force energy, respectively, ${\sigma}_{\theta}^{0}$, $A$, and $\rho $ represent the initial stress, the ring section area and density.
Hamilton principle can be expressed as:
where $\delta $ is the vibrational symbol.
Substituting Eqs. (2)(5) into Eq. (6), and rearranging the terms yields the final equations of motion:
$+\rho A\left(\ddot{v}{\ddot{v}}^{\text{'}\text{'}}4\mathrm{\Omega}{\dot{v}}^{\text{'}}+{\mathrm{\Omega}}^{2}\left({v}^{\text{'}\text{'}}v\right)\right){k}_{r}{v}^{\text{'}\text{'}}+{k}_{v}v={q}_{r}^{\text{'}}+{q}_{v}+\frac{1}{R}{q}_{\beta}^{\text{'}\text{'}}+\frac{1}{R}{q}_{\beta},$
where the inextensible condition, i.e., $v\mathrm{\text{'}}+w=0$, is used. $E$ is effective modulus of the belt and $I$ the moment of inertia in Eq. (7), which means $EI$ is the effective bending stiffness of the ring.
The natural frequencies of the nonrotating wheel are used to determine the model parameters and validate the model. With $\mathrm{\Omega}$, ${q}_{r}$, ${q}_{v}$ and ${q}_{\beta}$ being set to zero, Eq. (7) can be written as:
Assume the free vibration mode in the sinusoidal series as follows:
Inserting Eq. (9) into Eq. (8) yields the following natural frequency expression:
The undamped vibration frequencies of the rotating MEW can be determined by setting ${q}_{r}$, ${q}_{v}$ and ${q}_{\beta}$ equal to zero and assuming a solution of the form:
where:
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:
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:
${k}_{r}{v}^{\text{'}\text{'}}+{k}_{v}v{c}_{r}{\dot{v}}^{\text{'}\text{'}}+{c}_{v}\dot{v}={q}_{r}^{\text{'}}+{q}_{v}+\frac{{q}_{\beta}^{\text{'}\text{'}}}{R}+\frac{{q}_{\beta}}{R}.$
As an approximation, the functions $v(\varphi ,t)$ and ${q\mathrm{\text{'}}}_{r}+{q}_{v}$ are expanded as Fourier series:
where:
Substituting Eq. (15), (16) and ${q}_{\beta}=0$ into Eq. (14) yields a set of linear secondorder ordinary differential equations:
where:
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:
where $\stackrel{}{F}\left(\varphi \right)$ 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:
where $m$ is the number of contact locations, $\eta $ denotes the delta function, and ${\varphi}_{ci}$ 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 ${\varphi}_{ci}$ imply that:
where $r\left(\varphi \right)$ denotes the vertical position of the ring at the angular location $\varphi $ due to concentrated forces at the angular locations ${\varphi}_{ci}$. Using the inextensibility assumption $w=v\mathrm{\text{'}}$, one obtains:
And at the point of contact:
Making use of Eqs. (17, 18) and Eqs. (26)(29) in the absence of vibration yields:
where:
Substituting Eq. (15) into Eq. (31), one obtains:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\sum}_{n=1}^{N}\left[{\alpha}_{n}\right(n\mathrm{s}\mathrm{i}\mathrm{n}n\varphi \mathrm{s}\mathrm{i}\mathrm{n}\varphi +\mathrm{c}\mathrm{o}\mathrm{s}n\varphi \mathrm{c}\mathrm{o}\mathrm{s}\varphi ){\beta}_{n}(n\mathrm{c}\mathrm{o}\mathrm{s}n\varphi \mathrm{s}\mathrm{i}\mathrm{n}\varphi \mathrm{s}\mathrm{i}\mathrm{n}n\varphi \mathrm{c}\mathrm{o}\mathrm{s}\varphi \left)\right].$
Differentiating Eq. (38) with respect to $\varphi $ yields:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\sum}_{n=1}^{N}({n}^{2}1)\left[{\alpha}_{n}\right(\mathrm{c}\mathrm{o}\mathrm{s}n\varphi \mathrm{c}\mathrm{o}\mathrm{s}\varphi n\mathrm{s}\mathrm{i}\mathrm{n}n\varphi \mathrm{s}\mathrm{i}\mathrm{n}\varphi )+{\beta}_{n}(\mathrm{s}\mathrm{i}\mathrm{n}n\varphi \mathrm{c}\mathrm{o}\mathrm{s}\varphi +n\mathrm{c}\mathrm{o}\mathrm{s}n\varphi \mathrm{s}\mathrm{i}\mathrm{n}\varphi \left)\right].$
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 ${\stackrel{}{F}}_{ci}$ and ${\varphi}_{ci}$. The exception is for cases in which there is only one contact location and the MEW is stationary ($\mathrm{\Omega}=$ 0) and undamped. The solution is given simply by ${\stackrel{}{F}}_{c1}=F/\left(Rb\right)$ and ${\varphi}_{c1}=0$.
Assume the perturbation variables ${\stackrel{}{\alpha}}_{n}$, ${\stackrel{}{\beta}}_{n}$, $\stackrel{}{F}$ and ${\stackrel{}{\varphi}}_{ci}$ of the static equilibrium as follows:
Linearizing the Eq. (19) and Eq. (20) about the static equilibrium yields:
where:
Linearizing the contact boundary conditions $r\left({\varphi}_{ci}\right)=r\left({\varphi}_{c,i+1}\right)$ and $r\mathrm{\text{'}}\left({\varphi}_{ci}\right)=0$ about the static equilibrium yields:
where:
Eq. (45), (46) and Eq. (54)(55) can be expressed equivalently as:
where:
$\widehat{\lambda}\left(2n+1,2n\right)={\gamma}_{n},\widehat{\lambda}\left(2n+1,2n+1\right)={c}_{n},$
$\widehat{\psi}\left(2n+1,2n\right)={\tau}_{n},\widehat{\psi}\left(2n+1,2n+1\right)={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 1. Model geometrical and structural parameters of MEW
Parameters type

Value

$\rho $

2.56×10^{3}^{}kg/m^{3}

$h$

0.05 m

$b$

0.317 m

$A$

0.01585 m^{2 }

$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 $EI$ 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 threedimensional 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 cordrubber composite which displays the mechanics anisotropy and nonlinearity, 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. 3. FEM model of the MEW
According to the results of finite element simulation, the inplane 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 2. Apparatus 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. 4. Simulation 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. 5. Modal 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^{7} N/m^{2}, ${k}_{v}=$ –1.898×10^{7} N/m^{2}, and $EI=$ 1218.22 N·m^{2}. 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. 6. Experimental 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. 7. Theoretical 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 realtimely 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 3. Natural 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{\zeta}_{0}\sqrt{{k}_{v}\rho h}$ and ${c}_{r}=2{\zeta}_{1}\sqrt{2({k}_{r}+{k}_{v})\rho 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. 8. Effect of vertical load on vibration frequency of radial mode of MEW
Fig. 9. Effect 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 firstorder form in which the rotational speed has no effect on the firstorder frequency. However, the frequencies of modes (2nd)(6th) all have a noticeable dependence on $\mathrm{\Omega}$, 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. 10. Effect of number of the hinges on natural frequency under the free suspension state
Fig. 11. Effect 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 4. Damping ratio of MEW radial modes under different number of hinges
Damping ratio (%)

Number of hinges


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.
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).
References
 Kozhevnikov I. F. The vibrations of a free and loaded tyre. Journal of Applied Mathematics and Mechanics, Vol. 70, Issue 2, 2006, p. 250256. [Publisher]
 Kozhevnikov I. F. Vibrations of a rolling tyre. Journal of Sound and Vibration, Vol. 331, Issue 7, 2012, p. 16691685. [Publisher]
 Huang S. C., Soedel W. Effects of coriolis acceleration on the free and forced inplane vibrations of rotating rings on elastic foundation. Journal of Sound and Vibration, Vol. 115, Issue 2, 1987, p. 253274. [Publisher]
 Huang S. C., Su C. K. Inplane dynamics of tires on the road based on an experimentally verified rolling ring model. Vehicle System Dynamics, Vol. 21, Issue 1, 1992, p. 247267. [Publisher]
 Bickford W. B., Reddy E. S. On the inplane vibrations of rotating rings. Journal of Sound and Vibration, Vol. 101, Issue 1, 1985, p. 1322. [Publisher]
 Kim W., Chung J. Free nonlinear vibration of a rotating thin ring with the inplane and outofplane motions. Journal of Sound and Vibration, Vol. 258, Issue 1, 2002, p. 167178. [Publisher]
 Kim S. J., Savkoor A. R. Contact problem of inplane rolling of tires on a flat road. Vehicle System Dynamics, Vol. 27, Issue 27, 1997, p. 189206. [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. 145162. [Publisher]
 Dohrmann C. R. Dynamics of a tirewheelsuspension assembly. Journal of Sound and Vibration, Vol. 210, Issue 5, 1998, p. 627642. [Publisher]
 Huang S. C., Hsu B. S. An approach to the dynamic analysis of rotating tirewheelsuspension units. Journal of Sound and Vibration, Vol. 156, Issue 3, 1992, p. 505519. [Publisher]
 Yam L. H., Guan D. H., Zhang A. Q. Threedimensional mode shapes of a tire using experimental modal analysis. Experimental Mechanics, Vol. 40, Issue 4, 2000, p. 369375. [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. 20352050. [Publisher]
 Gasmi Amir, Joseph Paul F., Rhyne Timothy B., et al. Closedform 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. 843853. [Search CrossRef]
 Gasmi Amir, Joseph Paul F., Rhyne Timothy B., et al. Development of a twodimensional model of a compliant nonpneumatic tire. International Journal of Solids and Structures, Vol. 49, Issue 13, 2012, p. 17231740. [Search CrossRef]
 Lee Chihun, Ju Jaehyung, Kim DooMan Vibration analysis of nonpneumatic tires with hexagonal lattice spokes. Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Chicago, USA, 2012. [Search CrossRef]
 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, SpringerVerlag, Berlin, Heidelberg, Vol. 10, 2013, p. 199209. [Search CrossRef]
 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. 16, (in Chinese). [Search CrossRef]
 Zhao Youqun, Zang Liguo, Chen Yueqiao, et al. Nonpneumatic mechanical elastic wheel natural dynamic characteristics and influencing factors. Journal of Central South University, Vol. 22, Issue 5, 2015, p. 17071715. [Publisher]
 Kint P., Sas P., Desmet W. Development and validation of a threedimensional ringbased structural tyre Model. Journal of Sound and Vibration, Vol. 326, Issues 35, 2009, p. 852869. [Search CrossRef]
 Kint P., Sas P., Desmet W. Measurement and analysis of rolling tire vibrations. Optics and Lasers in Engineering, Vol. 47, Issues 34, 2009, p. 443453. [Search CrossRef]
 Sugiyama Hiroyuki, Suda Yoshihiro Modeling nonlinear flexible tire belt in the study of inplane tire dynamics. Proceeding of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Nevada, USA, 2007. [Search CrossRef]
 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. 538556. [Search CrossRef]
Cited By
Mechanics of Advanced Materials and Structures
Taoyu Wu, Mingxuan Li, Xiaolei Zhu, Xiaofeng Lu

2020

Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering
Youqun Zhao, Han Xu, Yaoji Deng, Qiuwei Wang

2019

International Journal of Automotive Technology
YouQun Zhao, YaoJi Deng, Fen Lin, MingMin Zhu, Zhen Xiao

2018

AIP Advances
Youqun Zhao, Mingmin Zhu, Fen Lin, Zhen Xiao, Haiqing Li, Yaoji Deng

2018

Simulation Modelling Practice and Theory
Yaoji Deng, Youqun Zhao, Fen Lin, Zhen Xiao, Mingmin Zhu, Haiqing Li

2018
