Development of the Flexible Ring on an Elastic Continuous Foundation Tire Model for Planar Vibration of the Heavy Load Radial Tire

Featured Application

Tire model is suitable for vehicle dynamics simulation and tire design.

Abstract

This paper investigates the planar vibration characteristic of heavy load radial tires with a large flat ratio. A proposed tire model with a flexible ring on an elastic continuous foundation is investigated utilizing kinematic modeling and experimental modal analysis. Planar coupling deformation of the radial and tangential direction is considered to enrich the kinematic characteristic of the flexible belt and the continuous sidewall; a flexible ring on an elastic continuous foundation tire model is proposed to investigate the coupling vibration characteristic between the flexible belt and the continuous sidewall. In-extensibility assumption is utilized to simplify the proposed tire model and the planar vibration modal features of the heavy load radial tire are discussed. The variation of the inflation pressure on the radial and tangential stiffness of the sidewall spring model is enriched into the flexible ring on an elastic continuous foundation tire model to extend the modal prediction of the tires with a different inflation pressure. Taking the relative error between the experimental and analytical modal resonance frequency of the tested tire with a different inflation pressure as the object value, structural parameters of the proposed tire model are identified by a backward genetic algorithm. Experimental and theoretical results show that: the planar coupling vibration characteristic of the heavy load radial tire can be predicted precisely with the flexible ring on an elastic continuous foundation tire model; meanwhile, considering the linear variations of the radial and tangential sidewall stiffness due to the inflation pressure, the proposed tire model can be extended to analyze the vibration characteristic of the heavy load radial tire with a different inflation pressure.

1. Introduction

The challenges to engineer a vehicle with the smooth-riding [1], pleasant-sounding [2], and safe-maneuverability [3] characteristics are becoming increasingly difficult and emergently needed. Moreover, the assessments of the vehicle quality and comfort are greatly influenced by their expectations and perceptions of the vehicle Noise Vibration and Harness (NVH) [4], especially for the heavy load vehicle [5] at higher speed. Regarding the vehicle dynamics, the dynamic behavior of the tires is of utmost importance as the tires generating most of the external forces and moments acting on a vehicle for acceleration, steering, and braking conditions. Being the only element of the vehicle contacting with the road, the tire is a key factor regarding the force and vibration generation and transmission. Meanwhile, the accurate prediction of the tire dynamics and the ground-tire interaction [6] have always been of interest. Different tire models with the different model complexity levels have been developed and researched for different applications.

Simplified tire models with an analytical expression or purely empirically/semi-empirical formulae, such as the magic model [7], brush model [8], and LuGre model [9], etc., are used to investigate the longitudinal and lateral dynamics of the tire-road interactions with the computationally effective feature, while the structural elasticity and the flexible deformation of the tire belt and tire sidewall are ignored, leading to the existing deviation under certain manoeuvres. At the other extreme, taking the advantage the powerful computation and the composite characteristic [10], detailed tire models, such as the finite element tire model [6], absolute nodal coordinate formulation (ANCF) finite element tire model [11], and spectral element tire model [12], present the detail material feature and geometrical feature at the expense of much computation resource. Additionally, the detailed tire/road contacting feature can be investigated with the detailed tire and road model, while the large computation resource, the extensive finite element analysis (FEA)-oriented pre-processing, and the difficultly-obtained material parameters limit the extended application for vehicle dynamics.

With less computation resources required than the detailed tire model, and a higher precision than the simplified tire model, the mechanical structure-based tire model [13] is proposed and much research has been devoted to focus on the kinematic modeling and numerical solving. Being the typical tire model respective of the mechanical structure-based tire model, the flexible belt tire model simplifies the tire as the flexible belt on an elastic foundation (referred to as the tire sidewall) and can range from one dimension to three dimensions. The circumference vibration of the flexible belts can be featured as the one-dimensional tire model [14] (referred to as the radial deformation in the polar coordinate system), or the two-dimensional tire model [15] (referred to as the radial and tangential deformation in the polar coordinate system), while the cross-section vibration of the belts [16] is added to form the three-dimensional tire model (referred to as the radial, tangential, and lateral deformation [17] in polar coordinates). The stiffness dimension of the sidewall spring model is consistent with the belt dimension ranging from one dimension [18] to three dimensions [19]. The radial deformation resulting from the inflation pressure and bending stiffness of the flexible belt is validated to be the main influence factor within the low/medial frequency band.

The coupling vibration characteristic between the tire belt and the tire sidewall was revealed as the main factor influencing the planar characteristic of the heavy load radial tire with a large flat ratio. Figure 1 illustrates the coupling vibration between the flexible belt and the continuous sidewall with the same direction. The vibration characteristic of heavy load radial tire was theoretically and experimentally investigated in the prior paper with the flexible beam on an elastic continuous sidewall beam tire model [20], in which the one-dimensional vibration characteristic, due to the radial deformation of the flexible belt, is considered, while the radial and tangential deformation of the flexible belt and the continuous sidewall is ignored. Considering the coupling feature between the radial and tangential deformations of the flexible belt, the flexible ring on an elastic continuous sidewall ring tire model is proposed and developed from the flexible beam on an elastic continuous sidewall beam tire model.

In the present paper, the radial and tangential deformation, respectively, of the flexible belt ring and continuous sidewall ring are investigated to enrich the two-dimensional tire model for the heavy load radial tire with a large flat ratio. Additionally, with the stiffness variation of the sidewall spring model to the inflation pressure, the flexible ring on the elastic continuous foundation tire model is extended to the modal prediction of the tires with a different inflation pressure.

2. Theoretical Derivation of the Two Dimension Tire Model

The theoretical modeling of the heavy load radial tire with a large flat ratio was investigated with the flexible beam on an elastic continuous sidewall beam tire model [20] and planar experimental modal analysis was researched. To conclude, the prior research is introduced briefly in the following section.

2.1. Flexible Beam on an Elastic Continuous Sidewall Beam Tire Model

The schematic diagram of the flexible beam on an elastic continuous sidewall beam tire model is shown in Figure 2. The bending deformation of the flexile belt and the pre-tension axle force resulting from the inflation pressure was considered. Additionally, the sectional stiffness of the sidewall radial spring and the bending feature of the continuous sidewall beam were added to enrich the one-dimensional tire model. $u_r$ and $u_{sr}$ are, respectively, the horizontal displacement of the flexible belt beam and the continuous sidewall beam at time $t$ at point $\theta$; $EI$ and $E_s{I}_s$ are, respectively, the bending stiffness of the flexible belt beam and the continuous sidewall beam; $\rho A$ is the density per line of the belt beam and $m_s$ is the density per rad of the continuous sidewall beam; $k_{r1}$ and $k_{r2}$ are the radial stiffness of the sectional sidewall spring model; $R$ and $R_s$ are, respectively, the belt radii and sidewall ring radii; $F_N$ is the pre-tension force due to the inflating pressure and $F_N=PRb$, where $P$ is inflation pressure; and $b$ is belt width.

Equation (1) illustrates the coupling kinematic equations of the flexible beam on an elastic continuous sidewall beam tire model: A. Equation (1a) shows the bending equation of flexible belt beam; B. Equation (1b) implies the bending equation of continuous sidewall beam; C. Equation (1c) illustrates the vibration equation of wheel.

$$\{\begin{array}{ll}\frac{EI}{R^4}\frac{{\partial }^4{u}_r}{\partial {\theta }^4}+k_{r1}\cdot (u_r-u_{sr\theta })+\rho A\frac{{\partial }^2{u}_r}{\partial t^2}-\frac{Pb}{R}\frac{{\partial }^2{u}_r}{\partial {\theta }^2}=0& (\mathrm{a})\\ \frac{E_s{I}_s}{R_s{}^4}\frac{{\partial }^2{u}_{sr}}{\partial {\theta }^2}+m_s\frac{{\partial }^2{u}_{sr}}{\partial t^2}=k_{r1}\cdot (u_r-u_{sr})-k_{r2}(u_{sr}-u_{wr}\cos \theta )& (\mathrm{b})\\ m_w{\ddot{u}}_{wr}={\displaystyle {\int }_0^{2\cdot \pi }[k_{r2}\cdot (u_{sr}-u_{wr}\cos \theta )}\cos \theta ]\mathrm{d}\theta & (\mathrm{c})\end{array}$$

(1)

The radial deformation of the flexible belt and the continuous sidewall is considered in the flexible beam on an elastic continuous sidewall tire model, while the coupling deformation characteristic of the continuum mediums, such as the tire belt and the tire sidewall between the radial and tangential deformations exists. Thus, the flexible ring on an elastic continuous foundation tire model is proposed and developed with the radial and tangential deformation of flexible belt ring and the continuous sidewall ring.

2.2. Flexible Ring on an Elastic Continuous Sidewall Ring Tire Model

The schematic diagram of the proposed two dimension tire model for the heavy load radial tire with a large flat ratio is illustrated in Figure 3. Meanwhile, the three-dimensional coordinate system, named as OXYZ, is shown, in which XOY shows the planar circumferential vibration characteristic of the heavy load radial tire, and YOZ and XOZ illustrate the cross-sectional and out-of-plane vibration characteristic of the heavy load radial tire.

The main focus of the paper is the circumferential vibration of the heavy load radial tire, so the theoretical modeling and experimental modal analysis of the heavy load radial tire with a large flat ratio is limited within XOY.

With the coupling feature between the radial and tangential deformation of the flexible belt ring and the continuous sidewall ring, the two dimension tire model with a flexible ring on an elastic continuous foundation is developed from the flexible beam on an elastic continuous sidewall beam tire model.

Some important assumptions are highlighted as:

(1)

The cross-sections of the flexible belt ring and the continuous sidewall ring keep plane feature and is perpendicular to the middle section during the tire deformation;

(2)

The temperature and the inflation pressure keep stable during the tire deformation; and

(3)

Considering the decoupling feature between the out of-plane vibration and the planar circumferential vibration, the out of-plane vibration of the tire is ignored.

The flexible ring on continuous sidewall ring tire model consists of four parts:

• The radial and tangential stiffness of the sectional sidewall spring model;
• The continuous sidewall ring (marked in blue in Figure 3);
• The flexible belt ring (marked in green in Figure 3); and
• The wheel with the concentrating mass.

The radial stiffness of the tire sidewall is modeled with the radial spring $k_{sr1}$ (connecting the flexible belt ring and the continuous sidewall ring) and $k_{sr2}$ (connecting the continuous sidewall ring and the wheel). Additionally, the tangential stiffness of the tire sidewall is featured with the tangential spring $k_{s\theta 1}$ (connecting the flexible belt ring and continuous sidewall ring) and $k_{s\theta 2}$ (connecting the continuous sidewall ring and the tire wheel).

A polar coordinate system with the origin point located at the ring center is adapted to the theoretical model. Two polar coordinate systems, named as $XOY$ and $X_1{O}_1{Y}_1$, are defined and respectively refer to the polar coordinate system with a non-rotating (fixed) and a rotating coordinate origin frame in Figure 4. $P$,$P_s$ are the arbitrary points of, respectively, the flexible belt and the continuous sidewall. $u_{\theta }$, $u_r$ respectively refer to the tangential and radial deformation of $P$ at rotating coordinate $X_1{O}_1{Y}_1$, while $u_{s\theta }$, $u_{sr}$ respectively refer to the tangential and radial deformation of $P_s$ at rotating coordinate $X_1{O}_1{Y}_1$.

The interval angles between $P$-axis and $X$-axis, $P$-axis and $X_1$-axis are respectively $\phi$ and $\theta$. $\mathsf{\Omega }$ is the angular velocity of the rotating ring with respect to the fixed center $O$/$O_1$”:

$$\phi =\theta +\mathsf{\Omega }t$$

(2)

The transformation between the coordinate $(x,y)$ of the rotating coordinate frame and coordinate $(x^{*},y^{*})$ of the non-rotating (fixed) coordinate frame is referred as:

$$\{\begin{array}{c}x=x^{*}\cos (\mathsf{\Omega }t)-y^{*}\sin (\mathsf{\Omega }t)\\ y=x^{*}\sin (\mathsf{\Omega }t)+y^{*}\cos (\mathsf{\Omega }t)\end{array}$$

(3)

2.3. Governing Equations of Hamilton Variation Principle

The kinematic equation of the rotating flexible ring on an elastic continuous foundation tire model can be obtained by applying the Hamilton principle [21]:

$$\delta {\displaystyle \underset{t_1}{\overset{t_2}{\int }}(U-T-W)dt}=0$$

(4)

2.3.1. Mechanical Energy Analysis of the Tire/Wheel Assembly

The kinetic energy of the flexible belt ring, the tire wheel and the continuous sidewall ring can be referred to as $T_c$, $T_w$, and $T_s$ in Equation (5):

$$\{\begin{array}{ll}{T}_c={\displaystyle \underset{0}{\overset{2\pi }{\int }}\frac{1}{2}}\rho AR{\left|\overrightarrow{R}\right|}^{2}d\theta & (\mathrm{a})\\ T_w=\frac{1}{2}{m}_w({\dot{x}}^2+{\dot{y}}^2)+\frac{1}{2}{I}_w{(\mathsf{\Omega }+\dot{\gamma })}^2& (\mathrm{b})\\ T_s={\displaystyle \underset{0}{\overset{2\pi }{\int }}\frac{1}{2}}{m}_s{R}_s{\left|{\overrightarrow{R}}_s\right|}^{2}d\theta & (\mathrm{c})\end{array}$$

(5)

where $\rho A$ is the line density of the flexible belt ring; $m_s$ is the rad density of the continuous sidewall ring; $m_w$ is the wheel mass; $I_w$ is the wheel moment; $R$ and $R_s$ are, respectively, the belt radii and sidewall ring radii; and $\gamma$ is the rotating angle of the tire wheel.

The whole kinetic energy of the tire/wheel assembly is calculated as Equation (6):

$$T=T_c+T_w+T_s$$

(6)

The potential energy of the tire/wheel assembly, including strain energy of flexible belt ring and continuous sidewall ring and the spring energy of the sectional sidewall spring respectively refer to $U_c$, $U_s$, and $U_{sp}$.

The pre-tension strain of flexible belt ring resulting from inflation pressure is taken into consideration and the strain energy of flexible belt ring is shown in Equation (7) [22]:

$$U_c=b{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{-\frac{h}{2}}{\overset{\frac{h}{2}}{\int }}(\frac{1}{2}{\sigma }_r{\epsilon }_r+\frac{1}{2}{\sigma }_{\theta }{\epsilon }_{\theta }+{\sigma }_{\theta }^0{\epsilon }_{\theta })}}Rdzd\theta$$

(7)

where the first part of Equation (7) is the radial strain energy of the flexible belt ring and ${\epsilon }_r=\frac{\partial u_r}{\partial R}=0$; the second part of Equation (7) is the circumferential strain energy of the flexible belt ring and the nonlinear strain ${\epsilon }_{\theta }$:${\epsilon }_{\theta }=\frac{(u_{\theta }{}^{\prime }+u_r)}{R}+\frac{1}{2}{(\frac{u_{\theta }{}^{\prime }+u_r}{R})}^2+\frac{1}{2}{(\frac{u_r{}^{\prime }-u_{\theta }}{R})}^2$ [23]; and the third part of Equation (7) is the initial strain energy of flexible belt ring and ${\sigma }_{\theta }^0$ is the initial strain resulting from the inflation pressure and the centrifugal force caused by rotation and calculated as Equation (8) in Figure 5:

$${\sigma }_{\theta }^{0}A=T_0=({\displaystyle \underset{0}{\overset{\pi }{\int }}(\rho A{\mathsf{\Omega }}^{2}R+p_{0}b)Rd\theta })\sin \theta$$

(8)

Thus, ${\sigma }_{\theta }^{0}A=\rho A{\mathsf{\Omega }}^2{R}^2+p_{0}bR$.

The strain energy of continuous sidewall ring $U_s$ is shown in Equation (9):

$$U_s={\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{-\frac{h}{2}}{\overset{\frac{h}{2}}{\int }}(\frac{1}{2}{\sigma }_{sr}{\epsilon }_{sr}+\frac{1}{2}{\sigma }_{s\theta }{\epsilon }_{s\theta }+{\sigma }_{s\theta }^0{\epsilon }_{s\theta })}}{R}_{s}dzd\theta$$

(9)

where the first part of Equation (9) is the radial strain of continuous sidewall ring; the second part of Equation (9) is the circumferential strain energy of continuous sidewall ring and the nonlinear strain ${\epsilon }_{s\theta }$:${\epsilon }_{s\theta }=\frac{(u_{s\theta }{}^{\prime }+u_{sr})}{R_s}+\frac{1}{2}{(\frac{u_{s\theta }{}^{\prime }+u_{sr}}{R_s})}^2+\frac{1}{2}{(\frac{u_{sr}{}^{\prime }-u_{s\theta }}{R_s})}^2$; the third part of Equation (9) is the initial strain energy of continuous sidewall ring; and ${\sigma }_{s\theta }^0$ is the initial strain resulting from the centrifugal force caused by rotation, ${\sigma }_{s\theta }^0{A}_s=m_s{\mathsf{\Omega }}^2{R}_s{}^2$.

The spring energy of the radial and tangential sectional sidewall stiffness:

$$U_{sp}={\displaystyle \underset{0}{\overset{2\pi }{\int }}\begin{array}{l}(\frac{1}{2}{k}_{sr1}{(u_r-u_{sr})}^2+\frac{1}{2}{k}_{s\theta 1}{(u_{\theta }-u_{s\theta })}^2+\frac{1}{2}{k}_{sr2}{(u_{sr}-x^{*}\cos \theta -z^{*}\sin \theta )}^2\\ +\frac{1}{2}{k}_{s\theta 2}{(u_{s\theta }+x^{*}\sin \theta -z^{*}\cos \theta -r\gamma )}^2)\end{array}}Rd\theta$$

(10)

The potential energy of tire/wheel assembly is summed up as Equation (11):

$$U=U_c+U_s+U_{sp}$$

(11)

2.3.2. Work Done by the Applied Forces

The work done by the applied forces consists of two parts, including $W_1$ due to the inflation pressure, and $W_2$ related with $q_{\theta }$ and $q_r$ as illustrated in Figure 4.

$$W=W_1+W_2$$

(12)

$W_1$ can refer to the product of the middle section area of the flexible belt and the ring width, in which the middle section area of flexible belt is calculated as:

$$\begin{array}{l}\tilde{A}=\frac{1}{2}{\displaystyle \oint (-y_{b}d{x}_b+x_{b}d{y}_b)}\\ =\frac{1}{2}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left\{\begin{array}{l}\left(R\sin \phi +u_{\theta }\cos \phi +u_r\sin \phi +y_w\right)[(R+u_{\theta }{}^{\prime }+u_r)\sin \phi -(u_r{}^{\prime }-u_{\theta })\cos \phi ]\\ +\left(R\cos \phi -u_{\theta }\sin \phi +u_r\cos \phi +x_w\right)[(R+u_{\theta }{}^{\prime }+u_r)\cos \phi +(u_r{}^{\prime }-u_{\theta })\sin \phi ]\end{array}\right\}}d\phi \end{array}$$

(13)

So, $\tilde{A}=\frac{1}{2}{\displaystyle \underset{0}{\overset{2\pi }{\int }}(R^2+u_r{}^2+u_{\theta }{}^2+2u_{r}R+u_{\theta }{}^{\prime }{u}_r-u_{\theta }{u}_r{}^{\prime })}d\phi$ and $W_2$ is obtained as:

$$W_1=\frac{p_{0}Rb}{2}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left[2u_r+\frac{1}{r}(u_r{}^2+u_{\theta }{}^2++u_{\theta }{}^{\prime }{u}_r-u_{\theta }{u}_r{}^{\prime })\right]}d\phi$$

(14)

where, $\frac{u_r}{r}\ll 1$, $\frac{u_r{}^{\prime }}{r}\ll 1$, $\frac{u_{\theta }{}^{\prime }}{r}\ll 1$ and Equation (14) can be simplified as:

$$W_1={\displaystyle \underset{0}{\overset{2\pi }{\int }}{p}_{0}rb{u}_r}d\phi$$

(15)

$W_2$, related with $q_{\theta }$ and $q_r$, is derived as:

$$W_2=b{\displaystyle \underset{0}{\overset{2\pi }{\int }}(q_r{u}_r+q_{\theta }{u}_{\theta })Rdzd\theta }$$

(16)

2.3.3. Kinematic Dynamics of Flexible Ring on an Elastic Continuous Sidewall Ring Tire Model

Equation (5) (kinetic energy of the tire/wheel assembly), Equation (11) (potential energy of the tire/wheel assembly) and Equation (12) are substituted into Equation (4) and the independent variables, including $u_r$, $u_{\theta }$, $u_{sr}$, $u_{s\theta }$, $x^{*}$, $y^{*}$, and $\gamma$, are manipulated. With the means of the variation principle, the planar coupling kinematic equations between the flexible belt and the continuous sidewall are obtained as:

$$\{\begin{array}{ll}\frac{EI}{R^4}(u_r{}^{(4)}-u_{\theta }{}^{(3)})+\frac{EA}{R^2}(u_r+u_{\theta }{}^{\prime })+\frac{{\sigma }_{\theta }^{0}A}{R^2}(u_{\theta }{}^{\prime }-u_r{}^{″})+k_{sr1}(u_r-u_{sr})-\frac{p_{0}b}{R}(u_r+u_{\theta }{}^{\prime })& \\ +\rho A({\ddot{u}}_r-2\mathsf{\Omega }{\dot{u}}_{\theta }-{\mathsf{\Omega }}^2{u}_r)=q_r& (\mathrm{a})\\ \frac{EI}{R^4}(u_r{}^{(3)}-u_{\theta }{}^{″})-\frac{EA}{R^2}(u_r{}^{\prime }+u_{\theta }{}^{″})+\frac{{\sigma }_{\theta }^{0}A}{R^2}(u_{\theta }-u_r{}^{\prime })+k_{s\theta 1}(u_{\theta }-u_{s\theta })+\frac{p_{0}b}{R}(u_r{}^{\prime }-u_{\theta })& \\ +\rho A({\ddot{u}}_{\theta }+2\mathsf{\Omega }{\dot{u}}_r-{\mathsf{\Omega }}^2{u}_{\theta })=q_v& \left(\mathrm{b}\right)\\ \frac{E_s{I}_s}{R_s{}^4}(u_{sr}{}^{(4)}-u_{s\theta }{}^{(3)})+\frac{E_s{A}_s}{R_s{}^2}(u_{sr}+u_{s\theta }{}^{\prime })+\frac{{\sigma }_{s\theta }^0{A}_s}{R_s{}^2}(u_{s\theta }{}^{\prime }-u_{sr}{}^{″})+m_s({\ddot{u}}_{sr}-2\mathsf{\Omega }{\dot{u}}_{s\theta }-{\mathsf{\Omega }}^2{u}_{sr})& \\ +k_{sr2}(u_{sr}-y^{*}\sin \theta -x^{*}\cos \theta )-k_{sr1}(u_r-u_{sr})=0& \left(\mathrm{c}\right)\\ \frac{E_s{I}_s}{R_s{}^4}(u_{sr}{}^{(3)}-u_{s\theta }{}^{″})-\frac{E_s{A}_s}{R_s{}^2}(u_{sr}{}^{\prime }+u_{s\theta }{}^{″})+\frac{{\sigma }_{s\theta }^0{A}_s}{R_s{}^2}(u_{s\theta }-u_{sr}{}^{\prime })+m_s({\ddot{u}}_{s\theta }+2\mathsf{\Omega }{\dot{u}}_{s\theta }-{\mathsf{\Omega }}^2{u}_{s\theta })& \\ +k_{s\theta 2}(u_{s\theta }-y^{*}\cos \theta +x^{*}\sin \theta -\gamma r)-k_{s\theta 1}(u_{\theta }-u_{s\theta })=0& \left(\mathrm{d}\right)\\ m_w({\ddot{x}}^{*}-2\mathsf{\Omega }{\dot{y}}^{*}-{\mathsf{\Omega }}^2{x}^{*})+\pi R_s(k_{sr2}+k_{s\theta 2})x^{*}-R_s{\displaystyle \underset{0}{\overset{2\pi }{\int }}(k_{sr2}{u}_{sr}\cos \theta -k_{s\theta 2}{u}_{s\theta }\sin \theta )d\theta }=F_{wx}& \left(\mathrm{e}\right)\\ m_w({\ddot{y}}^{*}+2\mathsf{\Omega }{\dot{x}}^{*}-{\mathsf{\Omega }}^2{y}^{*})+\pi R_s(k_{sr2}+k_{s\theta 2})y^{*}-R_s{\displaystyle \underset{0}{\overset{2\pi }{\int }}(k_{sr2}{u}_{sr}\sin \theta +k_{s\theta 2}{u}_{s\theta }\cos \theta )d\theta }=F_{wy}& \left(\mathrm{f}\right)\\ I_w\ddot{\gamma }+2\pi k_{s\theta 1}{R}_s{}^3\gamma -R_s{}^2{\displaystyle \underset{0}{\overset{2\pi }{\int }}{k}_{s\theta 2}{u}_{s\theta }d\theta }=T_w& \left(\mathrm{g}\right)\end{array}$$

(17)

where Equation (17a) gives the radial vibration equation of the flexible belt ring; Equation (17b) shows the tangential vibration equation of the flexible belt ring; Equation (17c) gives the radial vibration equation of the continuous sidewall ring; Equation (17d) presents the tangential vibration equation of the continuous sidewall ring; Equation (17e,f,g) respectively illustrate the horizontal, vertical, and rotating movement of the wheel.

3. Planar Modal Characteristic of Heavy Load Radial Tire

The planar modal experiment of the heavy load radial tire with a large flat ratio was implemented to investigate the coupling feature between the flexible belt ring and continuous sidewall ring.

3.1. System Implementation of the Experimental Modal Analysis

Figure 6b shows the schematic diagram of the planar modal experiment of the heavy load radial tire and an experimental hammer test is implemented, the principal features of which are shown in Figure 6a. The system implementation of the experimental modal analysis consists of three main parts, including the excitation part, response part and the data analysis part.

Vibration excitation is realized using a B&K hammer sensor (Naerum, Denmark), with the input force of the hammer sensor being measured by a force transducer. The PCB accelerometer sensors (PCB Piezotronics, Inc., Depew, NY, USA) are used to measure the resulting acceleration of the tire belt, the tire sidewall and the tire wheel illustrated in Figure 6c, in which the acceleration sensors are attached to the tire with a thin layer of mounting wax. Figure 6d presents the hardware system of the experimental modal analysis, in which by measuring the excitation force with the hammer sensor and the vibration response with the acceleration sensor, the input excitation force and output response acceleration are gathered using high-speed data acquisition equipment, namely a DEWE-43 (Austria).

The experimental transfer function is calculated using spectrum analysis by combining the input excitation force and the output response acceleration. Two additional methods were implemented to further improve the data quality for the measured transfer function:

(1)

Linear averaging was used to calculate the transfer function for every excitation point by averaging 10 repetitions of the transfer function measurement; and

(2)

The coherence function was checked to evaluate the quality of the measured transfer function to assess the measurement process for every excitation.

The frequency transfer functions for a heavy-load radial tire with a large aspect ratio are shown in Figure 7 to illustrate the vibration characteristics, in which, Figure 7a illustrates the three-dimensional planar transfer function between the excitation at point A and the response. Figure 7b presents the coherence value calculated from the input force and output acceleration response as a means of checking the planar transfer function; furthermore the following conclusion can be made on the basis of Figure 7:

(1)

The three-directional transfer function and coherence values for the excitation force and the response acceleration compared in Figure 7a show that the coherence value of the radial excitation force and the radial acceleration response have a difference of 1 for frequencies below 13 Hz. There is a similar difference of 1 between the coherence value for the radial excitation force and the tangential acceleration response for frequencies below 50 Hz. Beyond the above frequency band, the coherence value of the radial excitation force and the lateral acceleration response fluctuate around the frequency of 75 Hz; and

(2)

The larger amplitude of the radial transfer function than the tangential and lateral transfer functions implies that planar radial deformation of the flexible belt is the main vibration source for heavy-load radial tires with a large aspect ratio.

In addition to this, the sectional planar transfer function within 300 Hz can be divided into two parts: 0–160 Hz and 160–280 Hz. The sectional transfer function results for the coupling vibration between the flexible belt and the continuous sidewall, moreover, indicates that coupling vibration between the flexible belt and the continuous sidewall generally needs to be considered when undertaking vibration analysis of heavy-load radial tires with a large aspect ratio.

(3)

The structure-cavity coupling vibration peaks at 120 Hz (the first-order structure-cavity resonance frequency) and 240 Hz (the second-order structure-cavity resonance frequency) as shown in Figure 7a. For comparison, the first-order resonance frequency for tires with a small aspect ratio, including bus and truck tires, is typically within 230–250 Hz.

Meantime, the experimental transfer functions of the heavy load radial tire with a different inflation pressure are investigated in Figure 8.

3.2. Modal Resonance Parameters Identification

The least square complex exponential method is utilized for the modal parameters identification, including the modal resonance frequency, modal damping, and modal shape.

Figure 7 and Figure 8 present the frequency domain vibration characteristic, meanwhile the impulse is needed. Utilizing the inverse FFT method, the impulse response function $h(t)$ is given as:

$$h(t)=\mathrm{Re}\left({\displaystyle \sum _{r=1}^{2N}{A}_r{e}^{s_{r}t}}\right)$$

(18)

where $\mathrm{Re}$ is the real part.

The impulse response function of the finite time form $t_k=k\Delta t$ transforming from the impulse response function is given as:

$$h(k\Delta t)={\displaystyle \sum _{r=1}^{2N}{A}_r{e}^{s_{r}k\Delta t}}={\displaystyle \sum _{r=1}^{2N}{A}_r{V}_r^k}$$

(19)

where $V_r=e^{s_r\Delta t}$,$k=0,1,\cdots ,L$.

Impulse response functions of the different time can be concluded with the finite time form:

$$\{\begin{array}{c}{h}_0={\displaystyle \sum _{r=1}^{2N}}{A}_r{V}_r^0=A_1+A_2+\cdots +A_{2N}\\ h_1={\displaystyle \sum _{r=1}^{2N}}{A}_r{V}_r^1=A_1{V}_1+A_2{V}_2+\cdots +A_{2N}{V}_{2N}\\ h_2={\displaystyle \sum _{r=1}^{2N}}{A}_r{V}_r^2=A_1{V}_1^2+A_2{V}_2^2+\cdots +A_{2N}{V}_{2N}^2\\ ⋮\\ h_L={\displaystyle \sum _{r=1}^{2N}}{A}_r{V}_r^L=A_1{V}_1^L+A_2{V}_2^L+\cdots +A_{2N}{V}_{2N}^L\end{array}$$

(20)

where ${\beta }_{2N}=1$.

$${\displaystyle \sum _{r=1}^{2N}}{\beta }_k{V}_{}^k={\displaystyle \prod _{r=1}^N}(V-V_k)(V-V_k^{\ast })=0$$

(21)

Multiplying Equation (21) by ${\beta }_k$, the calculated equations are summed up:

$${\displaystyle \sum _{k=1}^{2N}}{\beta }_k{h}_k={\displaystyle \sum _{k=1}^{2N}}{\beta }_k{\displaystyle \sum _{r=1}^{2N}}{A}_r{V}_r^k={\displaystyle \sum _{r=1}^{2N}}{A}_r{\displaystyle \sum _{k=1}^{2N}}{\beta }_k{V}_r^k$$

(22)

Substituting ${\displaystyle \sum _{k=0}^{2N}}{\beta }_k{V}_r^k=0$ and ${\beta }_{2N}=1$, Equation (22) is transformed as:

$${\displaystyle \sum _{k=1}^{2N-1}}{\beta }_k{h}_k=-h_{2N}$$

(23)

Constructing Equation (23) with the time domain result to solve the autoregressive coefficient ${\beta }_k$, the equation is set with the same time interval $\Delta t$:

$$\{\begin{array}{c}{\displaystyle \sum _{k=1}^{2N-1}}{\beta }_k{h}_k={\beta }_0{h}_0+{\beta }_1{h}_1+\cdots +{\beta }_{2N-1}{h}_{2N-1}=-h_{2N}\\ {\displaystyle \sum _{k=1}^{2N-1}}{\beta }_k{h}_{k+1}={\beta }_0{h}_1+{\beta }_1{h}_2+\cdots +{\beta }_{2N-1}{h}_{2N}=-h_{2N+1}\\ \\ {\displaystyle \sum _{k=1}^{2N-1}}{\beta }_k{h}_{k+M-1}={\beta }_0{h}_{M-1}+{\beta }_1{h}_M+\cdots +{\beta }_{2N-1}{h}_{L-1}=-h_L\end{array}$$

(24)

in which, $M=L-2N$.

Transforming Equation (24) in the matrix form, the equation set is given as Equation (25):

$$\left[\begin{array}{ccccc}{h}_0& h_1& h_2& \cdots & h_{2N-1}\\ h_1& h_2& h_3& \cdots & h_{2N}\\ ⋮& ⋮& ⋮& & ⋮\\ h_{M-1}& h_M& h_{M+1}& \cdots & h_{L-1}\end{array}\right]\left[\begin{array}{c}{\beta }_0\\ {\beta }_1\\ ⋮\\ {\beta }_{2N-1}\end{array}\right]=\left[\begin{array}{c}{h}_{2N}\\ h_{2N+1}\\ ⋮\\ h_L\end{array}\right]$$

(25)

Utilizing the least square method, the autoregressive coefficient ${\beta }_k$ is gained:

$$\left\{\beta \right\}={({[h]}^T[h])}^{-1}([h]\{h^{\prime }\})$$

(26)

$${\displaystyle \sum _{r=1}^{2N}}{\beta }_k{V}_{}^k={\beta }_0+{\beta }_1{V}_{}^{}+{\beta }_2{V}_{}^2+\cdots +{\beta }_{2N-1}{V}_{}^{2N-1}+V_{}^{2N}=0$$

(27)

In which $V_r$ can be given by solving the polynomial equation formed with the autoregressive coefficient ${\beta }_k$. The modal parameters, including modal resonance frequency ${\omega }_r$ and modal damping ${\xi }_r$ are calculated in Equation (28) and modal shape coefficient $A_r$ is obtained in Equation (29):

$$\{\begin{array}{c}{R}_r=\ln V_r=s_r\Delta t\\ {\omega }_r=\left|R_r\right|/\Delta t\\ {\xi }_r=\sqrt{\frac{1}{1+{(\mathrm{Im}(R_r)/\mathrm{Re}(R_r))}^2}}\end{array}$$

(28)

$$\left[\begin{array}{ccccc}1& 1& 1& \cdots & 1\\ V_1& V_2& V_3& \cdots & V_{2N}\\ ⋮& ⋮& ⋮& & ⋮\\ V_1^L& V_2^L& V_3^L& \cdots & V_{2N}^L\end{array}\right]\left[\begin{array}{c}{A}_1\\ A_2\\ ⋮\\ A_{2N}\end{array}\right]=\left[\begin{array}{c}{h}_0\\ h_1\\ ⋮\\ h_L\end{array}\right]$$

(29)

The radial and tangential transfer functions are taken as the object transfer function; moreover, utilizing the least square complex exponential method, the modal resonance frequency and modal damping of the heavy load radial tire is identified and the identified transfer function can be calculated. Figure 9 shows the compared results between the experimental and identified transfer functions of the heavy load radial tire with a different inflation pressure. The coupling modal characteristic between the flexible belt ring and the continuous sidewall ring is identified and presented in Figure 10, including modal resonance frequency (Figure 10a), modal damping (Figure 10b), and modal shape (Figure 10c).

The compared transfer functions illustrated in Figure 9 show that the discrepancy between the identified transfer functions and the experimental transfer functions of the heavy load radial tire with a different inflation pressure is small, which validates that the least square complex exponential method can satisfy the identification precision as the sixteen orders resonance frequencies existing within 300 Hz.

To conclude:

(1)

The sixteen orders modal parameters of the heavy load radial tire are identified within 300 Hz;

(2)

The two typical features of the heavy load radial tire resulting from the large flat ratio appear as:

• The sectional tendency of the modal resonance frequency (Figure 10c,d); and
• The same/opposite direction (Figure 10c,d) between the flexible belt and continuous sidewall.

(3)

The modal shapes of the heavy load radial tire is harmonic, while the vibration direction of flexible belt and continuous sidewall is different in the low and medium frequency band. The modal shape of the heavy load radial tire within the low frequency is the coupling vibration of flexible belt and continuous sidewall with the same direction, while the modal shape within the medium frequency band is the coupling vibration of flexible belt and continuous sidewall with the opposite direction.

4. Analytical Coupling Modal of the Flexible Belt and the Continuous Sidewall

4.1. Simplification Forms of the Kinematic Equation

The in-extensibility assumption [23] (the circumferential length of the middle surface of the flexible ring, such as the tire belt and the tire sidewall, is constant during tire deformation) is utilized to simplify the dynamic equation of flexible belt:

$$u_{\theta }{}^{\prime }+u_r=0$$

(30)

Substituting Equation (30) into Equation (17a,b), then differentiating Equation (17a) with respect to $\theta$ and adding it to Equation (17b), the kinematic equation of the flexible belt ring with the in-extensibility belt ring assumption is transformed as:

$$\begin{array}{l}-\frac{EI}{R^4}(u_{\theta }{}^{(6)}+2u_{\theta }{}^{(4)}+u_{\theta }{}^{″})+\frac{{\sigma }_{\theta }^{0}A}{R^2}(u_{\theta }{}^{(4)}+2u_{\theta }{}^{″}+u_{\theta })-k_{sr1}(u_{\theta }{}^{″}+u_{sr}{}^{\prime })+k_{s\theta 1}(u_{\theta }-u_{s\theta })-\frac{p_{0}b}{R}(u_{\theta }+u_{\theta }{}^{″})\\ +\rho A({\ddot{u}}_{\theta }-{{\ddot{u}}^{″}}_{\theta }-4\mathsf{\Omega }{{\dot{u}}^{\prime }}_{\theta }+{\mathsf{\Omega }}^2(u_{\theta }{}^{″}-u_{\theta }))={q^{\prime }}_r+q_{\theta }\end{array}$$

(31)

Moreover, the in-extensibility assumption is also used to simplify the dynamic equations of the continuous sidewall ring (the circumferential length of the middle surface of the flexible belt ring and the continuous sidewall ring is constant during the deformation).

The vibration characteristic of the continuous sidewall ring is same or opposite to the flexible belt ring and the in-extensibility assumption can be used to simplify the dynamic equation of the continuous sidewall ring with the similar form as:

$$u_{s\theta }{}^{\prime }+u_{sr}=0$$

(32)

The radial and tangential equations of the continuous sidewall ring are transformed as:

$$\{\begin{array}{ll}-\frac{EI}{R^4}(u_{\theta }{}^{(6)}+2u_{\theta }{}^{(4)}+u_{\theta }{}^{″})+\frac{{\sigma }_{\theta }^{0}A}{R^2}(u_{\theta }{}^{(4)}+2u_{\theta }{}^{″}+u_{\theta })-k_{sr1}(u_{\theta }{}^{″}-{u^{″}}_{s\theta })+k_{s\theta 1}(u_{\theta }-u_{s\theta })& \\ -\frac{p_{0}b}{R}(u_{\theta }+u_{\theta }{}^{″})+\rho A({\ddot{u}}_{\theta }-{{\ddot{u}}^{″}}_{\theta }-4\mathsf{\Omega }{{\dot{u}}^{\prime }}_{\theta }+{\mathsf{\Omega }}^2(u_{\theta }{}^{″}-u_{\theta }))={q^{\prime }}_r+q_{\theta }& \left(\mathrm{a}\right)\\ -\frac{E_s{I}_s}{R_s{}^4}(u_{s\theta }{}^{(6)}+2u_{s\theta }{}^{(4)}+u_{s\theta }{}^{″})+\frac{{\sigma }_{s\theta }^0{A}_s}{R_s{}^2}(u_{s\theta }{}^{(4)}+2u_{s\theta }{}^{″}+u_{s\theta })-k_{sr2}{u}_{s\theta }{}^{″}+k_{s\theta 2}(u_{s\theta }-\gamma r)& \\ +(k_{sr2}+k_{s\theta 2})(x^{*}\sin \theta -y^{*}\cos \theta )-k_{sr1}(-u_{\theta }{}^{″}+u_{s\theta }{}^{″})-k_{s\theta 1}(u_{\theta }-u_{s\theta })& \\ +m_s[{\ddot{u}}_{s\theta }-{\ddot{u}}_{s\theta }{}^{″}-4\mathsf{\Omega }{{\dot{u}}^{\prime }}_{s\theta }+{\mathsf{\Omega }}^2(u_{s\theta }{}^{″}-u_{s\theta })]=0& \left(\mathrm{b}\right)\\ m_w({\ddot{x}}^{*}-2\mathsf{\Omega }{\dot{y}}^{*}-{\mathsf{\Omega }}^2{x}^{*})+\pi R_s(k_{sr2}+k_{s\theta 2})x^{*}+R_s{\displaystyle \underset{0}{\overset{2\pi }{\int }}(k_{sr2}{\dot{u}}_{s\theta }\cos \theta k_{s\theta 2}{u}_{s\theta }\sin \theta )d\theta }=F_{wx}& \left(\mathrm{c}\right)\\ m_w({\ddot{y}}^{*}+2\mathsf{\Omega }{\dot{x}}^{*}-{\mathsf{\Omega }}^2{y}^{*})+\pi R_s(k_{sr2}+k_{s\theta 2})y^{*}-R_s{\displaystyle \underset{0}{\overset{2\pi }{\int }}(-k_{sr2}{\dot{u}}_{s\theta }\sin \theta +k_{s\theta 2}{u}_{s\theta }\cos \theta )d\theta }=F_{wy}& \left(\mathrm{d}\right)\\ I_w\ddot{\gamma }+2\pi k_{s\theta 1}{R}_s{}^3\gamma -R_s{}^2{\displaystyle \underset{0}{\overset{2\pi }{\int }}{k}_{s\theta 2}{u}_{s\theta }d\theta }=T_w& \left(\mathrm{e}\right)\end{array}$$

(33)

The harmonic vibration characteristic of the heavy load radial tire is influenced by the material and geometrical parameters of the tire, moreover is independent with the wheel motion. The wheel is assumed to be fixed and non-rotating, so:

$$x^{*}=0,y^{*}=0\mathrm{and}\gamma =0$$

(34)

The kinematic equations of the proposed flexible ring on an elastic continuous foundation tire model are transformed as:

$$\{\begin{array}{ll}-\frac{EI}{R^4}(u_{\theta }{}^{(6)}+2u_{\theta }{}^{(4)}+u_{\theta }{}^{″})+\frac{{\sigma }_{\theta }^{0}A}{R^2}(u_{\theta }{}^{(4)}+2u_{\theta }{}^{″}+u_{\theta })-k_{sr1}(u_{\theta }{}^{″}-{u^{″}}_{s\theta })& \\ +k_{s\theta 1}(u_{\theta }-u_{s\theta })-\frac{p_{0}b}{R}(u_{\theta }+u_{\theta }{}^{″})+\rho A({\ddot{u}}_{\theta }-{{\ddot{u}}^{″}}_{\theta })={q^{\prime }}_w+q_v& \left(\mathrm{a}\right)\\ -\frac{E_s{I}_s}{R_s{}^4}(u_{s\theta }{}^{(6)}+2u_{s\theta }{}^{(4)}+u_{s\theta }{}^{″})-k_{sr2}{u}_{s\theta }{}^{″}+k_{s\theta 2}{u}_{s\theta }-k_{sr1}(-u_{\theta }{}^{″}+u_{s\theta }{}^{″})& \\ -k_{s\theta 1}(u_{\theta }-u_{s\theta })+m_s({\ddot{u}}_{s\theta }-{\ddot{u}}_{s\theta }{}^{″})=0& \left(\mathrm{b}\right)\end{array}$$

(35)

4.2. Modal Solution of the Planar Coupling Vibration between the Flexible Belt and the Continuous Sidewall

The rotating feature and the movement of wheel is taken into the completing kinematic equations (Equation (17)) and the simplified forms of the proposed tire model with the In-extensibility assumption are derived as Equation (35). Planar coupling vibration characteristic between the flexible belt ring and the continuous sidewall ring would be further investigated and solved with modal expansion method. The free vibration mode of the flexible belt and the continuous sidewall is assumed in the sinusoidal series illustrated in Equation (36), in which, Equation (36a,b) respectively refer to the circumference vibration modal of the flexible belt ring and the continuous sidewall ring:

$$\{\begin{array}{ll}{u}_{\theta }(\theta ,t)={\displaystyle \sum _1^{\propto }{V}_{\theta ,n}\cos (n(\theta -\phi ))\sin {\omega }_{n}t}& (\mathrm{a})\\ u_{s\theta }(\theta ,t)={\displaystyle \sum _1^{\propto }{V}_s{}_{\theta ,n}\cos (n(\theta -\phi ))\sin {\omega }_{n}t}& (\mathrm{b})\end{array}$$

(36)

Substituting Equation (36) into Equation (35), the planar dynamics of the heavy load radial tire with a large flat ratio are given as:

$$\left[\begin{array}{cc}-\rho A(1+n^2){\omega }^2+K_{11}& K_{12}\\ K_{21}& -m_s(1+n^2){\omega }^2+K_{22}\end{array}\right]\left[\begin{array}{c}{V}_{\theta ,n}\\ V_{s\theta ,n}\end{array}\right]=0$$

(37)

where $K_{11}=\frac{EI}{R^4}(n^6-2n^4+n^2)+\frac{p_{0}b}{R}(n^4-n^2)+k_{sr1}{n}^2+k_{s\theta 1}$, $K_{12}=K_{21}=-k_{sr1}{n}^2-k_{s\theta 1}$, $K_{22}=\frac{E_s{I}_s}{R_s{}^4}(n^6-2n^4+n^2)+k_{sr1}{n}^2+k_{s\theta 1}+k_{sr2}{n}^2+k_{s\theta 2}$.

The kinematic equation (Equation (37)) is reorganized as:

$$\rho A{m}_s{(1+n^2)}^2{\omega }^4-\left[\rho A(1+n^2)K_{22}+m_s(1+n^2)K_{11}\right]{\omega }^2+K_{11}{K}_{22}-K_{12}{K}_{21}=0$$

(38)

The roots of Equation (38) are obtained as Equation (39). Meanwhile, the deformation relationship between the flexible belt ring and the continuous sidewall ring refers to Equation (40):

$${\omega }_n{}^2=\frac{\left[\rho A{K}_{22}+m_s{K}_{11}\right]\pm \sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}}{2\rho A{m}_s(1+n^2)}$$

(39)

$$\frac{V_{\theta ,n}}{V_{s\theta ,n}}=\frac{K_{12}}{\rho A(1+n^2){\omega }^2-K_{11}}=\frac{m_s(1+n^2){\omega }^2-K_{22}}{K_{21}}$$

(40)

Analytical modal resonance frequency of the flexible belt ring and the continuous sidewall ring is derived as:

$$f_n{}_{1,n2}=\raisebox{1ex}{${\omega }_n{}_{1,n2}$}\!\left/ \!\raisebox{-1ex}{$2\pi $}\right.$$

(41)

The deformation relationship between the flexible belt ring and the continuous sidewall ring is concluded as:

$$f=f_n{}_1=\raisebox{1ex}{$\sqrt{\frac{\left[\rho A{K}_{22}+m_s{K}_{11}\right]-\sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}}{2\rho A{m}_s(1+n^2)}}$}\!\left/ \!\raisebox{-1ex}{$2\pi $}\right.,$$

$$\frac{V_{\theta ,n}}{V_{s\theta ,n}}=\frac{2m_s{K}_{12}}{\left[\rho A{K}_{22}-m_s{K}_{11}\right]+\sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}}$$

(42)

$$f=f_n{}_2=\raisebox{1ex}{$\sqrt{\frac{\left[\rho A{K}_{22}+m_s{K}_{11}\right]+\sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}}{2\rho A{m}_s(1+n^2)}}$}\!\left/ \!\raisebox{-1ex}{$2\pi $}\right.,$$

$$\frac{V_{\theta ,n}}{V_{s\theta ,n}}=\frac{2m_s{K}_{12}}{\left[\rho A{K}_{22}-m_s{K}_{11}\right]-\sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}}$$

(43)

The coupling relationship of the flexible belt and the continuous sidewall is discussed by comparing the two roots of Equation (29).

The Equation (32) is the identical equation and substituted into Equations (42) and (43) respectively.

$$\sqrt{{\left[\rho A{K}_{22}+m_s{K}_{11}\right]}^2-4\rho A{m}_s(K_{11}{K}_{22}-K_{12}{K}_{21})}>\left|\rho A{K}_{22}-m_s{K}_{11}\right|$$

(44)

Substituting Equation (44) into Equations (42) and (43), the conclusion is made that:

(1)

$f=f_n{}_1$, substituting Equation (44) into Equation (42) and $\frac{V_{\theta ,n}}{V_{s\theta ,n}}>0$;

(2)

$f=f_{n2}$, substituting Equation (44) into Equation (43) and $\frac{V_{\theta ,n}}{V_{s\theta ,n}}<0$;

(3)

The positive ratio of $V_{\theta ,n}$ and $V_{s\theta ,n}$ with $f=f_n{}_1$ implies the same vibration direction of flexible belt ring and continuous sidewall ring, which is consistent with Figure 10c,d and the conclusion given by the flexible beam on an elastic continuous sidewall beam tire model [23]; and

(4)

The negative ratio of $V_{\theta ,n}$ and $V_{s\theta ,n}$ with $f=f_n{}_2$ implies the opposite vibration direction of flexible belt ring and the continuous sidewall ring, which is consistent with Figure 10c,d, and the conclusion given by the flexible beam on an elastic continuous sidewall beam tire model [23].

5. Structural Parameters Identification

With the disassembly of the heavy load radial tire (Figure 11), some geometrical parameters can be measured directly; meanwhile, some material parameters of the proposed tire model cannot be measured directly and the intelligent algorithm [24] is utilized to identify the unknown material parameters. The parameters of the proposed tire model are listed in

Table 1. Geometrical and structural parameters of the GL073A tire.

	Parameters	Symbol	Unit	Value
Geometrical parameters	Belt width	$b$	m	0.35
	Inflation pressure	$P$	N/m2	8 × 105
	Tire radii	$R$	m	0.65
Structural parameters	Density per rad of sidewall	$m_s$	kg/rad	10
	Density per line of belt	$\rho A$	kg/m	19.64
	Radial spring connecting the sidewall and belt	$k_{sr1}$	N/m	unknown
	Radial stiffness connecting the sidewall and wheel	$k_{sr2}$	N/m	unknown
	Tangential springconnecting the sidewall and belt	$k_{s\theta 1}$	N/m	unknown
	Tangential stiffness connecting the sidewall and wheel	$k_{s\theta 2}$	N/m	unknown
	Bending stiffness of flexible belt ring	$EI$	N/m	unknown
	Bending stiffness of continuous sidewall ring	$E_s{I}_s$	N/m	Unknown

. Taking the error between the experimental modal resonance frequency (Figure 10a modal resonance frequency of heavy load radial tire with normal inflation pressure) and analytical modal resonance frequency (Equation (41)) as the object value, a genetic algorithm (GA) is utilized to identify the unknown structural parameters illustrated in Table 1.

The optimization procedure is initialized with a population size of 20, 1000 generations, a generation gap of 0.9, a crossover rate of 0.7, and a mutation rate of 0.1. The major components of the GA are the encoding scheme, fitness evaluation, parent selection, crossover, and mutation operators, illustrated in Figure 12.

Figure 13 presents the structural parameters of the optimization procedure. The identified result utilizing the flexible ring on an elastic continuous foundation tire model is compared with that using the flexible beam on an elastic continuous sidewall beam tire model and the compared results are shown in

Table 2. Identified parameters of the GL073A tire.

Symbol	Unit	Value Identified by Flexible Ring on an Elastic Continuous Sidewall Ring Tire Model	Value Identified by Flexible Beam on an Elastic Continuous Sidewall Beam Tire Model
$k_{sr1}$	N/m	7.111 × 106	6.967 × 106
$k_{sr2}$	N/m	3.957 × 106	3.325 × 106
$k_{s\theta 1}$	N/m	1.087 × 106	╳
$k_{s\theta 2}$	N/m	9.003 × 106	╳
$EI$	N/m	48.670	23.109
$E_s{I}_s$	N/m	38.479	48.753

.

The identification results of Figure 13 and Table 2 imply that:

(1)

The object value of the parameters identification procedure (Figure 13a) is convergent to the steady point with the optimization procedure proceeding to 200 generations;

(2)

The structural parameters, including: bending stiffness $EI$,$E_s{I}_s$, and sidewall stiffness $k_{sr1}$, $k_{sr2}$, $k_{s\theta 1}$ and $k_{s\theta 2}$, converge to the optimum solution as the optimization procedure proceeding 400 generations;

(3)

The steady object value of identification procedure utilizing the flexible ring on an elastic continuous sidewall ring tire model is compared with that utilizing the flexible beam on an elastic continuous sidewall beam model in Figure 13a. The compared identification procedure implies that the model prediction error using flexible ring on an elastic continuous sidewall ring tire model is less than that utilizing the flexible beam on an elastic continuous sidewall beam model resulting from taking the coupling feature between the radial and tangential deformation of flexible belt ring and continuous sidewall ring into consideration;

(4)

Table 2 illustrates the identified structural parameters respectively using flexible ring on an elastic continuous sidewall ring tire model and flexible beam on an elastic continuous sidewall beam tire model. The structural parameter bending stiffness of the continuous sidewall ring $E_s{I}_s$, identified using the flexible ring on an elastic continuous sidewall ring tire model, is less and the structural parameters bending stiffness of flexible belt ring $EI$, and sidewall spring stiffness $k_{sr1}$, $k_{sr2}$ are larger than that using the flexible beam on an elastic continuous sidewall beam tire model;

(5)

Figure 14 shows the compared modal resonance frequency of experimental modal resonance frequency (Figure 14b) and the analytical modal resonance frequency, respectively by the flexible ring on an elastic continuous sidewall ring tire model and flexible beam on an elastic continuous sidewall beam tire model; the compared result shows that the error between the analytical modal frequency with normal inflation pressure of the flexible ring on an elastic continuous sidewall ring tire model, and the experimental modal frequency is limited within 3% and less than that of the analytical modal frequency of the flexible beam on an elastic continuous sidewall beam tire model.

6. Discussion

6.1. Variations of Inflation Pressure on the Planar Characteristic

Influence of the inflation pressure on the planar characteristic of the heavy load radial tire consists of two sections:

(1)

The pre-tension strain ${\sigma }_{\theta }^0$ of flexible belt ring resulting from the inflation pressure; and

(2)

The variation of the sidewall pre-tension force caused from inflation pressure.

Figure 15 shows the geometrical characteristic of the sidewall arc in which $R_s$ is the radius of the sidewall uniform arc; $2\theta$ is the interval angle; $l_s$ is the length of sidewall arc; $L_0$ is the initial vertical distance between the belt and rim; $O_1$ and $O_2$ are, respectively, the centers of the two-piece sidewall.

The variation of sidewall radial spring stiffness was discussed [25] and the analytical radial stiffness of sidewall was derived as:

$$k_{sr}=p_0\frac{\cos \theta +\theta \sin \theta }{\sin \theta -\theta \cos \theta }$$

(45)

where $k_{sr}$ is the radial stiffness of the sidewall spring.

Radial stiffness of the sidewall is sensitive with the inflation pressure and the geometrical feature. Related work has been done by other researchers and the tangential stiffness of the sidewall spring [25] is derived analytically as:

$$k_{s\theta }=\frac{p_0{l}_s\cos \theta }{b\theta L_0(1-L_0/4\theta )}$$

(46)

The intersection angle $\theta$ between the tire belt and the sidewall arc is assumed invariable as stretching and compressing deformation of sidewall and the radial stiffness $k_{sr}$ and the tangential stiffness $k_{s\theta }$ of the sidewall are proportional to the inflation pressure $p_0$.

Thus, radial and tangential stiffness of the tire sidewall with a different inflation pressure (illustrated in Equations (45) and (46)) refers to:

$$\{\begin{array}{cc}{k}_{sr1}=K_{sr1}+(p-p_0)*AA& (\mathrm{a})\\ k_{sr2}=K_{sr2}+(p-p_0)*BB& (\mathrm{b})\\ k_{s\theta 1}=K_{s\theta 1}+(p-p_0)*CC& (\mathrm{c})\\ k_{s\theta 2}=K_{s\theta 2}+(p-p_0)*DD& (\mathrm{d})\end{array}$$

(47)

where $K_{sr1}$, $K_{sr2}$ are, respectively, the radial stiffness of the tire sidewall with normal inflation pressure condition; $K_{s\theta 1}$, $K_{s\theta 2}$ are, respectively, the tangential stiffness of the tire sidewall with normal inflation pressure condition; $p_0$ is the normal inflation pressure; $AA$, $BB$ are, respectively, the proportional coefficients of sidewall radial stiffness; and $CC$, $DD$ are, respectively, the proportional coefficients of sidewall tangential stiffness.

6.2. Stiffness Coefficients Identification Using a GA

The proportional coefficients of sidewall stiffness sensitive with the inflation pressure are identified utilizing a GA with the modal prediction error of different inflation pressure as the object value. The optimization procedure is initialized with the same parameters list in Section 5 and the identification results of the identification procedure is presented in Figure 16.

Figure 16a presents the compared object value of flexible ring on an elastic continuous foundation tire model and flexible beam on an elastic continuous sidewall beam tire model. Additionally, Figure 16b shows the stiffness coefficients of the sidewall spring model of the identification procedure and the Figure 16c illustrates the radial and tangential stiffness of the sidewall spring model with a different inflation pressure.

Table 3. Identified stiffness coefficients sensitive to inflation pressure.

Symbol	Unit	Value	Symbol	Unit	Value
$AA$	/	−6.371	$CC$	/	−1.537
$BB$	/	−1.451	$DD$	/	−1.244

lists the identified stiffness coefficients of the sidewall spring model, which is sensitive with the inflation pressure.

Figure 17 presents the compared modal resonance frequency (inflation pressure: 0.7, 0.6, 0.5, 0.4, and 0.3 MPa) between the experimental modal resonance frequency, analytical modal resonance frequency, respectively utilizing the flexible ring on an elastic continuous foundation tire model and flexible beam on an elastic continuous sidewall beam tire model.

The identification results (Figure 16 and Figure 17 and Table 3) imply that:

(1)

The object value of the identification procedure (Figure 16a) is convergent to the steady point with the optimization procedure proceeding over 65 generations;

(2)

The steady object value of the identification procedure utilizing the proposed tire model is compared with that utilizing the flexible beam on an elastic continuous sidewall beam model in Figure 16a, and the compared identification procedure implies that the model prediction error using the proposed tire model is less than that utilizing the flexible beam on an elastic continuous sidewall beam model resulting from taking the coupling feature between the radial and tangential deformation of flexible belt ring and continuous sidewall ring into consideration;

(3)

The errors between the analytical modal frequency with a different inflation pressure (inflation pressure: 0.7, 0.6, 0.5, and 0.4 MPa illustrated in Figure 17) utilizing the flexible ring on an elastic continuous foundation tire model and experimental modal frequency are limited within 3%; meanwhile, the error of analytical modal resonance frequency with inflation pressure of 0.3 MPa is up to 6%; and

(4)

The predicted errors utilizing flexible ring on an elastic continuous foundation tire model are less than that utilizing the flexible beam on an elastic continuous beam tire model because of the coupling feature between the radial and tangential deformation of the flexible belt ring and continuous sidewall.

7. Conclusions

Utilizing the theoretical modeling and experimental modal analysis, flexible ring on an elastic continuous foundation tire model is proposed as the two-dimensional tire model to investigate the planar vibration characteristic of the heavy load radial tires with a large flat ratio. The planar coupling feature between the radial and tangential deformations of the flexible belt ring and the continuous sidewall ring is investigated. Additionally, the planar coupling vibration characteristic between the flexible belt ring and the continuous sidewall ring is modeled with a flexible ring on an elastic continuous foundation tire model. Furthermore, modal prediction error of this tire with the normal inflation pressure utilizing the parameters identified can be limited within 3%, and the simplified forms of the flexible ring on an elastic continuous foundation tire model are derived. Additionally, the stiffness variation of the sidewall spring model with respect to the inflation pressure is considered to extend the flexible ring on an elastic continuous foundation tire model. Lastly, structural parameters and the stiffness coefficients, which are sensitive with the inflation pressure, are identified utilizing a genetic algorithm (GA).

To conclude:

(1)

Considering the coupling feature of the flexible belt and the continuous sidewall, a flexible ring on an elastic continuous foundation tire models the character of the planar vibration of the heavy load radial tire with a larger flat ratio;

(2)

The pre-tension force variations of the flexible belt and the stiffness variations of the two-dimensional sidewall spring model with respect to the inflation pressure are considered and enriched into the flexible ring on an elastic continuous foundation tire model; and

(3)

The flexible ring on an elastic continuous foundation tire model with the stiffness variation of the sidewall spring model with inflation pressure can predict the modal frequency of the tires with a different inflation pressure, and the prediction errors of the planar modal frequency of the tires with a different inflation pressure are limited within 3% (inflation pressure: 0.7, 0.6, 0.5, and 0.4 MPa) and 6% (0.3 MPa).