Next Article in Journal
Research on Adaptive Control of Grinding Force for Carbide Indexable Inserts Grinding Process Based on Spindle Motor Power
Previous Article in Journal
An Alpha/Beta Radiation Mapping Method Using Simultaneous Localization and Mapping for Nuclear Power Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Numerical Method for Theoretical Tire Model Simulation

State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130025, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(9), 801; https://doi.org/10.3390/machines10090801
Submission received: 17 August 2022 / Revised: 31 August 2022 / Accepted: 8 September 2022 / Published: 11 September 2022
(This article belongs to the Section Vehicle Engineering)

Abstract

:
Theoretical tire models are often used in tire dynamics analysis and tire design. In the past, scholars have carried out a lot of research on theoretical model modeling; however, little progress has been made on its solution. This paper focuses on the numerical solution of the theoretical model. New force and moment calculation matrix equations are constructed, and different iterative methods are compared. The results show that the modified Richardson iteration method proposed in this paper has the best convergence-stability in the steady and unsteady state calculation, which mathematically solves the problem of nonconvergence of discrete theoretical models in the published reference. A novel discrete method for solving the total deformation of tires is established based on the Euler method. The unsteady characteristics of tire models are only related to the path frequency without changing its parameters, so the unsteady state ability of the tire model can be judged based on this condition. It shows that the method in the references have significant differences at different speeds with the same path frequency under turn slip or load variations input, but the method proposed in this paper has good results.

1. Introduction

For tire and vehicle dynamics simulations and analyses, different types of tire models have been developed and categorized as empirical and theoretical models [1,2]. The empirical model, or combined theoretical and empirical model, is used for vehicle dynamics, and its model parameters are obtained by fitting measurements [2], including the MF/PAC2002 [2,3,4,5], TMeasy [6,7,8], UniTire [9,10], Hankook-Tire [11], TameTire [12,13], MF-Swift [14,15], FTire [16,17], CDTire [18,19] and RMOD-K [20,21] models, etc.
However, compared with the theoretical tire model, the empirical model is not suitable for the mechanical analysis of tire dynamics and tire design due to having too many model parameters and being dependent on experimental data. Tire design engineers should find the optimal design parameters that can satisfy the engineering requirements simultaneously and, for more efficient product development, a reliable theoretical model is essential because it enables systematic parameter study and design optimization. The theoretical models have advantages over the finite-element models in terms of computational efficiency, and they can be incorporated into the parameter study and optimization framework more easily [22]. Furthermore, combined with the static simulation data of the finite-element, using the theoretical model to achieve the steady and unsteady simulation should be a good prospect for future tire virtual development.
Fromm and Julien [3,23] proposed the brush tire model, which laid the foundation for tread modeling and tire theoretical analysis. The beam [1,3,24,25,26,27,28,29,30,31,32,33,34,35,36] and string models [3] are two classical belt/carcass modeling methods. Finally, the latest beam model includes belt tension, flexural rigidity, shear rigidity [1,22,25,26] and lateral shear force distribution [1], and the relationship between the parameters of the beam model and the design parameters of belts is established [1,22].
H. Sakai [28] pointed out that the previous practice of assuming the shape of the contact-patch as rectangular, the pressure distribution as lateral uniform and the longitudinal parabolic distribution does not reflect the reality; therefore, he studied the influence of different load, inflation pressure, camber, side slip angle and lateral force on the contact length in detail, and established a semi-empirical model of contact length with different rib along the width. Through the study of the distribution of contact patch-pressure on each rib, the empirical formula of contact patch pressure on each rib with load was established [29]. Guo Konghui established the general function of pressure distribution based on a rectangular contact patch in which the convexity, uniformity and fore-aft shift of pressure-distribution were considered [37,38,39]. Patrick Gruber, Robin S. Sharp and Andrew D studied normal and shear forces in the contact based on the finite-element model and established a 2D semi-empirical contact patch, considering longitudinal grooves as a function of vertical load and camber angle [40]. Ch. Oertel established the geometric relationship between the contact patch and vertical load, camber angle and tread profile curvature [41]; however, the contact patch length calculated by this method will be larger than the actual one, and the contact patch area caused by camber will be larger. Nakajima, Y. and Hidano, S. proposed two-dimensional contact patch considering that shear force, contact shape and contact pressure distribution are changed in the fore–aft and lateral direction caused by the camber, fore–aft and lateral external force [1,42], and the deficiency is that the contact patch shape is simply expressed as a rectangle or trapezoid.
Analytical [3,39,43,44,45,46,47,48,49,50] and numerical methods are usually used to solve theoretical models. The advantage of the former is high efficiency, but it is only suitable for simple theoretical models. With the increase of the complexity of the theoretical model, the numerical method becomes more and more important. TreadSim [3] is a typical representative for the steady-state numerical solution; however, the iteration process turned out to be unstable for certain combinations of vertical load, belt yaw and belt bend stiffness [51], and the method proposed in [51] does not actually solve the problem essentially by iteration. For unsteady-state numerical solutions, there are two different (but essentially the same) methods that are both based on the Lagrange description. An element that leaves the contact patch at one side is redefined in the model so that it will enter the contact patch at the opposite side. After a new element has entered the contact patch, the elements are also renumbered so that element number one is always the first element in the contact patch, and when the traveling distance in an increment time is larger than the distance between adjacent tread elements, a linear-interpolation method is used [52,53,54,55]. However, this method has some problems under the condition of turn slip or load variations inputs, which will be discussed in this paper. There are many methods for the iterative solution of matrix equations, including the traditional [56,57] and the latest [58,59], and different methods are suited for different application scenarios. For dynamic equations of tire models which usually consider mass, Newton’s method and the Newmark-beta method are always used [60,61]. However, in this paper, we do not intend to focus on the development of new iterative methods, but to reconstruct the matrix equations that can meet the application requirements of different iterative methods and compare several common iterative methods to select the most suitable for the discrete theoretical model and make appropriate optimizations.
In this paper, a novel matrix equation for the force and moment calculation is constructed and different iterative methods are compared under steady-state constant side slip angle and unsteady-state conditions, including the step side slip, sine side slip, load variations at constant side slip, sine side slip angle at constant turn slip and sine longitudinal slip at constant side slip inputs. It shows that only the modified Richardson method proposed in this paper can complete all of the simulations without any convergence. A novel, discrete solution method for tire total deformation is proposed based on the Euler method, and time derivatives are replaced by spatial gradients. Furthermore, based on the condition of the unsteady characteristics of tire models only relating to the path frequency without changing their parameters, the discrete solution method proposed is compared with reference [52,53,54,55] under different unsteady-state conditions, including the step side slip angle, sine side slip angle, load variations at constant side slip angle, step turn slip simulation, sine longitudinal combined with a constant side slip, sine turn combined with a constant side slip and load variations at constant longitudinal and side slip inputs. It shows that the simulation results have little difference by using the proposed method, but the reference method has a significant difference under turn slip or load variations inputs at different speed or, more precisely, the ratio between the distance of the contact patch center traveling in space within the time interval and distance between adjacent tread elements in the contact patch. However, this paper only studies the numerical solution method of the theoretical model in detail, and the results of the comparison between the model and the experiment show that a refined contact patch model considering camber and the influence of lateral force and longitudinal force on the contact patch model demonstrate the need for more in-depth research on the dynamic friction model, considering that the influences of contact pressure and slip velocity are added to the theoretical model.
The sections of this paper are arranged as follows. In Section 2, the tire total deformation discrete solution method based on the Euler method will be given. In Section 3, the modeling method of belt/carcass deformation, contact patch, stress direction in the sliding region and anisotropic of tread stiffness is described. In Section 4, according to the adhesion and sliding region, the matrix equations for the calculation of force and torque are established, and the advantages and disadvantages of different iterative methods are compared. In Section 5, the theoretical model is compared with the experiment to verify the rationality of the theoretical model parameters, and based on the principle that the unsteady state only depends on the path frequency, the discrete solution method proposed in this paper is compared with reference [52,53] unsteady-state conditions. Finally, the discussion and conclusions will be given in Section 6.

2. Tire Total Deformation Discrete Solution Method Based on the Euler Method

The total deformation of tires caused by friction consists of tread and carcass. When the tire deformation is described, two coordinate systems are used; one is the road coordinate system ( O , X , Y , Z ) and the other is the contact patch coordinate system ( C , x , y , z ) , named as I S O W . When the tire slips relative to the road surface, the friction forces makes the tire deform ( u , v ) relative to the contact patch coordinates ( x , y ) , as shown in Figure 1. P is the contact point at undeformed state and p is the vector of P relative to the ( O , X , Y , Z ) system. C represents the contact patch center and c is the vector C relative to the ( O , X , Y , Z ) system. ( x , y ) represents the contact point vector at undeformed state, ( u , v ) represents the total tire deformation vector and q is the vector of P relative to C . α represents the sideslip angle and ψ represents the yaw angle of the x -axis with respect to the X -axis.
The position of the tread location relative to ( O , X , Y , Z ) can be expressed as the equation:
O P = O C + C P   or   p = c + q
Herein, q = ( x + u ) e x + ( y + v ) e y . So, there is the slip velocity of the tread relative to the road:
V g = p ˙ = c ˙ + q ˙ = V c + ( x ˙ + u ˙ ) e x + ( y ˙ + v ˙ ) e y + ψ ˙ { ( x + u ) e y ( y + v ) e x }
where V c is the contact patch center velocity and e x and e y are the ( C , x , y , z ) base vector. When there is no relative slip between the tread and the road surface, V g = 0 , then
u ˙ = V s x + ψ ˙ ( y + v )
v ˙ = V c y ψ ˙ ( x + u )
Herein, x ˙ = d x d t = Ω R e = V r , V s x = V c x + x ˙ , y ˙ = 0 , expression of Equations (3) and (4) by the Euler method
u ˙ = u x d x d t + u y d y d t + u S d S d t = V s x + ψ ˙ ( y + v )
v ˙ = v x d x d t + v y d y d t + v S d S d t = V c y ψ ˙ ( x + u )
Then, there are:
u x + u S V c V r = S x + φ ( y + v )
u x + u S V c V r = S x + φ ( y + v )
Herein, d S d t = V c , S x = V s x V r = κ κ + 1 , S y = V c y V r = 1 κ + 1 V c y V x = 1 κ + 1 tan α , φ = ψ ˙ V r , and κ is practical slip, α is side slip angle, S x is the longitudinal theoretical slip, S y is the lateral theoretical slip and φ is the turn theoretical slip. Equations (7) and (8) are only related to the coordinates of the contact patch ( x , y ) and the space position of the contact patch center. The expression of Equations (7) and (8) in a discrete method:
u { S k + 1 } ( x i + 1 , y j + 1 ) u { S k + 1 } ( x i , y j + 1 ) x i + 1 x i + u { S k + 1 } ( x i + 1 , y j + 1 ) u { S k } ( x i + 1 , y j + 1 ) S k + 1 S k V c V r = S x + φ [ y j + 1 + v { S k + 1 } ( x i + 1 , y j + 1 ) ]
v { S k + 1 } ( x i + 1 , y j + 1 ) v { S k + 1 } ( x i , y j + 1 ) x i + 1 x i + v { S k + 1 } ( x i + 1 , y j + 1 ) v { S k } ( x i + 1 , y j + 1 ) S k + 1 S k V c V r = S y φ [ x j + 1 + u { S k + 1 } ( x i + 1 , y j + 1 ) ]
According to Equations (9) and (10), there are tire total deformation:
u { S k + 1 } ( x i + 1 , y j + 1 ) = p u 1 [ p u 2 + φ 1 1 d s V c V r 1 d x p u 3 ]
v { S k + 1 } ( x i + 1 , y j + 1 ) = p v 1 [ p v 2 φ 1 1 d s V c V r 1 d x p v 3 ]
Herein:
d x = ( x i + 1 x i ) represents the distance between adjacent tread elements in the contact patch.
d s = S k + 1 S k represents the distance of the contact patch center traveling in space within the time interval d t .
p u 1 = p v 1 = 1 1 + ( φ 1 d s V c V r 1 d x ) 2 1 1 d s V c V r 1 d x
p u 2 = S x + φ y j + 1 u { S k + 1 } ( x i , y j + 1 ) 1 d x + u { S k } ( x i + 1 , y j + 1 ) 1 d s V c V r
p u 3 = S y φ x i + 1 v { S k + 1 } ( x i , y j + 1 ) 1 d x + v { S k } ( x i + 1 , y j + 1 ) 1 d s V c V r
p v 2 = S y φ x i + 1 v { S k + 1 } ( x i , y i + 1 ) 1 d x + v { S k } ( x i + 1 , y j + 1 ) 1 d x V c V r
p v 3 = S x + φ y j + 1 u { S k + 1 } ( x i , y j + 1 ) 1 d x + u { S k } ( x i + 1 , y j + 1 ) 1 d s V c V r
If V c V r is set to zero, it can be transformed into a steady state.

3. Model Description

3.1. Belt/Carcass Deformation [23,54]

Assuming that longitudinal translation deformation only occurs under longitudinal force [54]:
u c = F x K c x 0
where K c x 0 is the longitudinal translational stiffness of carcass. Under lateral force, the carcass is expressed as the deformation of an infinite beam under concentrated force. The beam has in-plane flexural rigidity E I z , where I z is the moment of the inertia of the area, E is the Young’s modulus of the belt in the circumferential direction, k s is the lateral fundamental spring rate per unit length in the circumferential direction and the tension T 0 is uniformly distributed in the belt width direction [23]. M represents bending moment and Q represents shear force, as shown in Figure 2. The left diagram shows the equilibrium of forces of the beam section and the right diagram shows model of the beam on the elastic foundation.
The equilibrium of forces acting on the element of length d x along the belt width results in the following equation:
w ( x ) d x + Q + d Q + T 0 ( y x + 2 y x 2 d x ) Q T 0 y x k s y d x = 0
and   for   the   Euler   Bernoulli   beam :   Q = d M d x
M = E I z d 2 y d x 2
Furthermore, according to (14)–(16), the differential equation of the lateral deformation of the beam can finally be expressed as:
E I z 4 y x 4 T 0 2 y x 2 + k s y = w ( x )
According to boundary conditions, (1) x , y ( x ) 0 ; (2) x = 0 , y ˙ ( x ) = 0 ; and (3) F y = 2 k s 0 y d x , there is:
y ( x ) = δ F y 4 k s e λ 1 x ( cos λ 2 x + λ 1 λ 2 sin λ 2 x )
Herein:
λ 1 = k s 4 E I z 4 ( 1 + T 0 4 E I z k s ) ,   λ 2 = k s 4 E I z 4 ( 1 T 0 4 E I z k s )
δ = λ 1 2 + λ 2 2 λ 1 ,   T 0 = ξ 4 E I z k s
Torsional deformation of carcass under aligning moment [54].
y ( x ) = M z N θ x
N θ is the carcass torsion stiffness.
Therefore, the lateral deformation of the carcass under lateral force and an aligning moment can be expressed as
v c = δ F y 4 k s e λ 1 x ( cos λ 2 x + λ 1 λ 2 sin λ 2 x ) + M z N θ x

3.2. Contact Patch

Reference [41] presents a method for calculating the shape of a contact patch based on vertical force. Due to the fact that this method usually leads to a large contact patch area following camber input, this paper does not consider camber.
G ( x i , y j ) = R l R 2 x i 2 R 2 ( y j R y ) n y
R l = R d
F z = p F z 1 d + p F z 2 d 2
where R is free radius, R l is the loaded radius, d is the vertical deflection, R y is the lateral radius and n y is the lateral direction curvature exponent.
However, directly according to this method, the obtained contact length will be larger. A correction method is proposed as follows:
G lim = G 1 d 2 + G 2 d
Then, the contact patch can be obtained when G ( x i , y j ) < G lim .
The pressure distribution of ground imprinting is expressed as follows [41,54]:
η ( x i h j ) = A ( 1 ( x i h j ) 2 n ) ( 1 + λ ( x i h j ) 2 n ) ( 1 B x i h j )
F z i j ( x i , y j ) = F 0 i j η ( x i h j ) ( 1 + a F z 1 ( y j b i ) 2 ( a F z 1 + 1 ) ( y j b i ) 6 )
assuming that F 0 i j is a constant, which can be obtained according to F z = F z i j ( x i , y j ) .
h j is the contact length at position y j , b i is the contact width at position x i and a F z 1 is the convexity factor along the contact width.
A = ( 2 n + 1 ) ( 4 n + 1 ) 2 n ( 4 n + 1 + λ ) ,   B = 3 ( 2 n + 3 ) ( 4 n + 3 ) ( 4 n + 1 + λ ) ( 2 n + 1 ) ( 4 n + 1 ) ( 4 n + 3 + 3 λ ) Δ h j
where λ is convexity factor of contact pressure along contact length, Δ h j is the contact pressure center-offset factor and n is the contact pressure-uniformity factor.
When considering the tread pattern longitudinal groove, the F z i j ( x i , y j ) of the corresponding position is set to zero. The contact patch example is shown as in Figure 3.

3.3. Direction of Stress in the Sliding Regions

As shown in Figure 4, B i , P i , T i are the upper (which is attached to the belt) and lower end of the bristle and the bristle deformation, respectively, at previous times. B i + 1 , P i + 1 , T i + 1 correspond to the current time. d G is the bristle slip displacement vector relative to the road.
d G = V g d t = ( V B + V T ) d t = d B + d T = d B + T i + 1 T i
then
d G T i + 1 = d B T i = P i B i + 1
For sliding regions, the sliding direction V g determines the direction of deformation, so the vector d G is in the same direction as T i + 1 and T i + 1 is on vector P i B i + 1 ; that is to say P i B i + 1 determines the direction of T i + 1 .

3.4. Stiffness of the Tread Element

The pattern block structure usually leads to the anisotropy of tread element stiffness. It is assumed that the stiffness satisfies the elliptic equation, as shown in Figure 5:
( k t cos θ k t x ) 2 + ( k t sin θ k t y ) 2 = 1
then:
k t = 1 ( cos θ k t x ) 2 + ( sin θ k t y ) 2
where θ is the deformation direction angle of the tread element.

4. Iterative Method for Force and Moment Calculation

4.1. Deformation and Stress Calculation of Tread Element

When the tread element is in the adhesion region, the deformation is calculated as follows
u T = u u c
v T = v v c
Herein, u and v are obtained from Section 2, Equations (11) and (12); u c and v c are obtained from Section 3.1, Equations (13) and (17). The deformation direction angle of the tread element is expressed as
tan θ = v T u T
The longitudinal and lateral stresses of tread element are expressed as
q x a = k t u T ,   q y a = k t v T ,   q a = q x a 2 + q y a 2
If q a μ q z
q x = q x a ,   q y = q y a
else
q x = q x s = q x a q a μ q z = μ q z cos θ
q y = q y s = q y a q a μ q z = μ q z sin θ
u T = q x k t ,   v T = q y k t ;   u = u T + u c ,   v = v T + v c
Above are the stress and deformation of tread element in the sliding region.

4.2. Calculation of Force and Moment

4.2.1. Fx Calculation

F x = q x d A = k t u T d A + q x s d A = k t u d A F x k t K c x 0 d A + q x s d A
then:
( 1 + p F x ) F x = F x r
where:
p F x = k t K c x 0 d A ,   F x r a = k t u d A ,   F x r s = q x s d A ,   F x r = F x r a + F x r s

4.2.2. Fy Calculation

F y = q y d A = k t v T d A + q y s d A = k t v d A F y k t η d A M z k t ς d A + q y s d A
then:
( 1 + p F y ) F y = F y r + M z p M T F
where:
η = δ 4 k s e λ 1 x [ cos ( λ 2 x ) + λ 1 λ 2 sin ( λ 2 x ) ] ,   ς = x N θ
p F y = k t η d A ,   p M T F = k t ς d A ,   F y r a = k t v d A ,   F y r s = q y s d A
F y r = F y r a + F y r s

4.2.3. Mz Calculation

M z = q y ( x + u ) d A q x ( y + v ) d A = k t v T ( x + u ) + q y s ( x + u ) k t u T ( y + v ) q x s ( y + v )
then:
( 1 + p M z ) M z = M z r + F y p F T M
where:
p M z = k t ς ( x + u ) d A ,   p F T M = k t η ( x + u ) d A
M z r x a = k t u T ( y + v ) d A ,   M z r x s = q x s ( y + v ) d A
M z r y a = k t v ( x + u ) d A ,   M z r y s = q y s ( x + u ) ,
M z r = M z r x a + M z r x s + M z r y a + M z r y s
According to Equations (33) and (35), there are
( 1 + p F y p F T M p M T F 1 + p M z ) F y = F y r + p M T F 1 + p M z M z r
( 1 + p M z p F T M p M T F 1 + p F y ) M z = M z r + p F T M 1 + p F y F y r
Make:
A 11 = 1 + p F x ,   A 22 = 1 + p F y p F T M p M T F 1 + p M z ,   A 33 = 1 + p M z p F T M p M T F 1 + p F y
B 1 = F x r ,   B 2 = F y r + p M T F 1 + p M z M z r ,   B 3 = M z r + p F T M 1 + p F y F y r
There are
A 11 0 0 0 A 22 0 0 0 A 33 F x F y M z = B 1 B 2 B 3
Make:
A = A 11 0 0 0 A 22 0 0 0 A 33 ,   B = B 1 B 2 B 3 ,   X = F x F y M z
Finally, there are
A X = B
A reasonable iterative method can be constructed according to this matrix equation.

4.3. Iterative Strategy

4.3.1. Matrix Splitting Method

Make:
A = M N
then:
X m + 1 = M 1 N X m + M 1 B
where N i j M i j = I R , I R < 1 and M i j = A i j 1 I R , N i j = I R 1 I R A i j .

4.3.2. Steepest-Descent Method

r m = B A X m p m = r m α m = r m T p m p m T A p m X m + 1 = X m + α m p m

4.3.3. Richrdson Method

X m + 1 = X m + ω ( B A X m )
where ω = 2 λ 1 + λ n , λ 1 and λ n are the maximum and minimum eigenvalues of A , respectively.

4.3.4. Extrapolation Acceleration Method

According to matrixes M and N from Section 4.3.1, the extrapolation acceleration method can be made:
X m + 1 = ( 1 ω ) X m + ω ( M 1 N X m + M 1 B )
where ω = 2 2 ( λ 1 + λ n ) , λ 1 and λ n are the maximum and minimum eigenvalues of M N , respectively.

4.3.5. Method in References [54,55]

This method obtains the solution directly from the matrix equation:
X m + 1 = A 1 B
This method is equivalent to Section 4.3.1 where I R is equal to zero.

4.3.6. Method in References [3,51]

This method is equivalent to that of the matrix M whose coefficients are equal to one, and the method is the same as in Section 4.3.1.

4.3.7. Iterative Error

E R R = B A X m + 1
If E R R T E R R T o l , the iterative reaches convergence, where T o l is the iterative error tolerance.

4.3.8. Comparison of Iterative Methods

A comparison of the iterative methods at a steady state under a constant side slip angle = three deg is shown in Figure 6, and method 4 and 5’s iterative processes coincide completely.
The unsteady comparison results of different iterative methods are shown in Table 1.
Figure 6 shows that method 6 has the fastest iteration speed at a steady state, followed by method 3. It can be seen from Table 1 that although method 3 has the slowest iteration speed in unsteady state, it has the most stable convergence. As such, method 3 is selected as the iterative strategy of the discrete theoretical model. However, it should be noted that in the case of unsteady large slip, due to the fact that the coefficients of the A matrix parameters tend to one but the B matrix changes, the solution is sometimes not convergent, so it is necessary to improve the iterative method. Through research, it is found that the following equation can improve convergence:
ω = 2 λ 1 + λ n + p λ λ 1 + λ n
The adjustment of p λ can better improve the convergence, but it will reduce a certain convergence speed. The greater p λ is, the better the convergence will be, but the iteration speed will be reduced, and this parameter can be reasonably selected for specific situations. This paper recomments p λ = 20 , and the simulation in Table 1 uses this parameter.

5. Simulation of the Discrete Theoretical Model

5.1. Model Parameters and Verification

The model parameters are listed in Table 2, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the comparison between the model and the test and the test information is shown in Table 3.
The contact pressure is performed on the TekScan. The other three conditions are performed on the MTS CT Flat-trac machine (MTS System Corporation). For the sideslip angle step condition, the test procedure is to turn the steer angle first, then load and finally apply the belt speed of the MTS to 10kph within one second. It should be noted that the sign of steer angle in the table is the opposite to that of the slip angle which is determined by the MTS Flat-trac machine.
It should be noted that the response of normalized lateral force instead of lateral force with travel distance is used, which does not affect the assessment of transient behavior.
As shown in Figure 11, the aligning moment is quite different from the experimental data, and a reasonable reason is whether the local camber of the carcass caused by the lateral force is not considered, which can lead to greater aligning moment if taken into consideration.
The unsteady characteristics of tire models are only related to the path frequency without changing its parameters, so the unsteady state-ability of the tire model can be judged based on this condition and changing the ratio between the distance of the contact patch center traveling in space within the time interval and the distance between adjacent tread elements in the contact patch by changing speed. The path frequency is defined as:
ω s = 2 π f V
where ω s [ m 1 ] represents the path frequency, f [ s 1 ] represents the time frequency and V [ m / s ] represents the contact patch center velocity.
Next, the discrete method proposed for tire total-deformation calculation in this paper (Section 2) will be compared with the references [52,53], and the comparison results of two methods are given. The comparison error calculation method is as follows:
e r r = ( y p y r ) 2 y p 2 × 100 %
where y p represents the proposed method in this paper and y r represents the proposed method in reference [52,53].

5.2. Pure Side Slip Simulation

5.2.1. Step Side Slip Angle

The simulation results of step side slip angle with the proposed method and the reference method are as shown in Figure 12.

5.2.2. Sine Side Slip Angle

The simulation results of sine side slip angle with the proposed method and the reference method are as shown in Figure 13.

5.2.3. Load Variations at Constant Side Slip Angle

The simulation results of load variations at constant side slip angle with the proposed method and the reference method are as shown in Figure 14.

5.2.4. Discussion

The above simulations have the same path frequency, and the simulation results for different speeds should also be the same, in theory. The simulation results of the proposed method are slightly different, but significant differences occur in reference methods under the load variations condition at different speeds. The simulation errors of the two methods at 3m/s are listed in Table 4.

5.3. Step Turn Slip Simulation

The simulation results of step turn slip with the proposed method and the reference method are as shown in Figure 15.
The method proposed in the references is obviously not suitable for the turn slip simulation because the steady-state simulation results at different speeds are quite different. The reason is that the method proposed in the reference is based on the principle of linear interpolation in the leading edge, while the tread deformation is nonlinear during turn slip and the method fails when d s is much larger than d x . At the same time, it can be seen that the aligning moment responds quickly in the initial position, which is explained in Section 5.5.2. The simulation errors of the two methods at 3m/s are listed in Table 5.

5.4. Combined Slip Simulation

5.4.1. Sine Longitudinal Combined with a Constant Side Slip

The simulation results of sine longitudinal combined with a constant side slip with the proposed method and the reference method are as shown in Figure 16.

5.4.2. Sine Turn Combined with a Constant Side Slip

The simulation results of sine turn combined with a constant side slip with the proposed method and the reference method are as shown in Figure 17.

5.4.3. Load Variations at Constant Longitudinal and Side Slip

The simulation results of load variations at constant longitudinal and side slip with the proposed method and the reference method are as shown in Figure 18.

5.4.4. Discussion

In the case of load variations or turn slip input the reference method has obvious problems, and the conclusions are consistent with Section 5.2 and Section 5.3. The simulation errors of the two methods at 3m/s are listed in Table 6.

5.5. Deformation Analysis of Tread Element

5.5.1. Deformation of Tread Element under Step Side Slip Angle

In the side slip angle step condition, with the increase of the travel distance, the deformation of the tread element approximately changes shape from rectangular to trapezoidal and finally to triangular, as shown in Figure 19.

5.5.2. Deformation of the Tread Element under Step Turn Slip

In the turn slip step condition, the lateral deformation is first symmetrically distributed relative to the longitudinal center of contact patch, which is the reason why the slope of lateral force with the travel distance is zero at the starting position and also explains why aligning moment responds so quickly, as shown in Figure 20.
With the increase of the travel distance, the intersection point of the lateral deformation curve with the contact patch longitudinal axis gradually shifts from the center to the intersection point of the steady state. In this process, the lateral force increases continuously, and the increment gradient of the aligning moment caused by the lateral deformation decreases, which leads to a peak of the aligning moment which then decreases and finally tends to the steady state value.
At the same time, the lateral deformation is approximate to the nonlinear change of the parabolic shape, which also explains why the reference method will have problems in the turn slip input simulation.

5.5.3. Deformation of the Tread Element under Load Variations at a Constant Side Slip Angle

Figure 21 shows the that the lateral deformation of the proposed method is more regular and stable than reference. It explains why under load variations conditions at the same path frequency, when the traveling distance in an increment time is larger than the distance between adjacent tread elements, the lateral force and aligning moment of the reference method will be different.

6. Discussion and Conclusions

6.1. Discussion

Figure 7 and Figure 8 show that the contact patch parameters are reasonable, and the simulation results are consistent with the experiments. Figure 9, Figure 10 and Figure 11a show that the belt/carcass parameters are reasonable, and simulation and test results are consistent as well. However, Figure 11b shows that the aligning moment is quite different from the experimental data, and a reasonable reason is that the local camber of the carcass caused by the lateral force is not considered, which can lead to greater aligning moment if considered.
Figure 6 shows that different iterative methods converge to the same value and that the iterative method in reference [54,55] has the fastest iteration speed, followed by the modified Richardson method proposed in this paper. Table 1 shows that the modified Richardson has the best convergence and then method in reference [54,55] and the worst is the method in reference [3,51].
Figure 12, Figure 13 and Figure 16 show that whether the discrete method proposed or the references [52,53,54,55] are considered, the simulation results at different speeds with the same path frequency are only slightly different under side slip angle or longitudinal slip input. Figure 14, Figure 15, Figure 17 and Figure 18 show that for the discrete proposed method the simulation results still have slight differences, but the results of the method in the references [52,53,54,55] have significant differences at different speeds with the same path frequency, and the greater the speed difference, or more precisely the greater the ratio between the distance of the contact patch center traveling in space within the time interval and distance between adjacent tread elements in the contact patch, the greater the simulation difference. The calculation errors in Table 4, Table 5 and Table 6 also confirm the above analysis. Figure 19 shows that the tread deformation varies nearly linearly at the leading edge of the contact patch under step side slip angle, but Figure 20 shows that the tread deformation varies nonlinearly in an approximate parabolic under step turn slip. Figure 21 shows that the tread deformation is more stable by using the proposed method under the load variations input. Tread deformation explains why the methods in the reference [52,53,54,55] have problems under the turn slip or load variations inputs, since the method performs a linear interpolation process at the leading edge which is not suitable for nonlinear deformation and deformation instability. The above conclusions may be able to explain the reasons for the reference [3] “at each time step in which the wheel is rolled further over a distance equal to the interval between two successive tread elements” simulation settings.

6.2. Conclusions

In this paper, novel numerical methods for discrete theoretical tire model steady- and unsteady-state simulations are proposed. Some conclusions are as follows.
The comparison between the theoretical model and the experiment shows that the accuracy of the aligning moment needs further improvement, which may be affected by tire local carcass camber caused by lateral force, but the model parameters are reasonable.
New basic matrix equations for force and moment calculation are constructed for the convenient application of different iterative methods. Different iterative methods are studied and compared. The results show that the new modified Richardson iteration method proposed in this paper has the best convergence stability in the unsteady calculation, which mathematically solves the problem of nonconvergence of discrete theoretical models in the published references [3,51,55].
A novel discrete method for solving the total deformation of tires is established based on the Euler method. The unsteady characteristics of tire models are only related to the path frequency without changing its parameters, so the unsteady-state ability of the tire model can be judged based on this condition. It shows that the methods in the references [52,53,54,55] have significant differences at different speeds with the same path frequency, and the greater the speed difference, or more precisely the greater the ratio between the distance of the contact patch center traveling in space within the time interval and distance between adjacent tread elements in the contact patch, the greater the simulation difference under the turn slip or load variations inputs. However, the method proposed in this paper has good results.
This paper only studies the numerical solution method of the theoretical model in detail. In the future, the following work needs to be carried out. The establishment of a refined contact patch model considering camber and the influence of lateral force and longitudinal force on the contact patch model and of a dynamic friction model considering the influence of contact pressure and slip velocity is added to the theoretical model. System verification of theoretical models with experiments. Furthermore, it is hoped that theoretical models can be used in tire virtual development, since the computational efficiency of the finite element model is low and virtual data simulation, especially for turn slip input because the existing test machines cannot carry out this condition. Combined with the static simulation data of the finite element, using the theoretical model to achieve the steady and unsteady simulations should be a good prospect for the future of tire virtual development.

Author Contributions

Conceptualization, Q.L. and D.L.; methodology, Q.L. and D.L.; software, Q.L.; validation, Q.L., Y.M. and D.X.; formal analysis, Q.L., Y.M. and D.X.; investigation, Q.L.; resources, D.L.; data curation, Q.L., Y.M. and D.X.; writing—original draft preparation, Q.L.; writing—review and editing, D.L., Y.M. and D.X.; visualization, Q.L.; supervision, D.L.; project administration, D.L.; funding acquisition, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Strategic Priority Research Program of Chinese Academy of Sciences, grant number XDC06060100.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

SymbolDescription
κ Practical longitudinal slip
α Side slip angle
S x Theoretical longitudinal slip
S y Theoretical lateral slip
φ Theoretical turn slip
d x Distance between adjacent tread elements
d s Distance of the contact patch center traveling in space within the time interval d t
K c x 0 Longitudinal translational stiffness of carcass
E I z In-plane flexural rigidity of belt
k s the lateral fundamental spring rate per unit length in the circumferential direction
N θ Carcass torsion stiffness
R Tire free radius
R l Tire loaded radius
d Tire vertical deflection
R y Tire lateral radius
n y Lateral direction curvature exponent
G 1 Contact patch shape correction quadratic coefficient
G 2 Contact patch shape correction primary coefficient
λ Convexity factor of contact pressure along contact length
n Contact pressure uniformity factor
Δ Contact pressure center offset factor
μ Friction coefficient
T o l Iterative error tolerance
k t x Tread element longitudinal stiffness
k t y Tread element lateral stiffness
k t Tread element stiffness
p F z 1 Convexity factor of contact pressure along the contact width
p λ Iterative correction coefficient
F y s s Steady-state lateral force

References

  1. Nakajima, Y.; Hidano, S. Theoretical Tire Model Considering Two-Dimensional Contact Patch for Force and Moment. Tire Sci. Technol. 2022, 50, 27–60. [Google Scholar] [CrossRef]
  2. Kuiper, E.V.O.J.; Van Oosten, J.J.M. The PAC2002 advanced handling tire model. Veh. Syst. Dyn. 2007, 45, 153–167. [Google Scholar] [CrossRef]
  3. Pacejka, H.B. Tire and Vehicle Dynamics, 3rd ed.; Elsevier: Oxford, UK, 2012. [Google Scholar]
  4. Pacejka, H.B.; Bakker, E. The magic formula tyre model. Veh. Syst. Dyn. 1992, 21, 1–18. [Google Scholar] [CrossRef]
  5. Pacejka, H.B.; Besselink, I.J.M. Magic formula tyre model with transient properties. Veh. Syst. Dyn. 1997, 27, 234–249. [Google Scholar] [CrossRef]
  6. Hirschberg, W.; Rill, G.; Weinfurter, H. Tire model tmeasy. Veh. Syst. Dyn. 2007, 45, 101–119. [Google Scholar] [CrossRef]
  7. Rill, G. TMeasy—A Handling Tire Model based on a three-dimensional slip approach. In Proceedings of the XXIII International Symposium on Dynamic of Vehicles on Roads and on Tracks (IAVSD 2013), Quingdao, China, 19–23 August 2013; pp. 19–23. [Google Scholar]
  8. Rill, G. Road Vehicle Dynamics: Fundamentals and Modeling, 1st ed.; Taylor & Francis: Boca Raton, FL, USA, 2011. [Google Scholar]
  9. Guo, K.; Lu, D.; Chen, S.K.; Lin, W.C.; Lu, X.P. The UniTire model: A nonlinear and non-steady-state tyre model for vehicle dynamics simulation. Veh. Syst. Dyn. 2005, 43 (Suppl. S1), 341–358. [Google Scholar] [CrossRef]
  10. Guo, K.; Lu, D. UniTire: Unified tire model for vehicle dynamic simulation. Veh. Syst. Dyn. 2007, 45, 79–99. [Google Scholar] [CrossRef]
  11. Gim, G.; Choi, Y.; Kim, S. A semiphysical tyre model for vehicle dynamics analysis of handling and braking. Veh. Syst. Dyn. 2005, 43 (Suppl. S1), 267–280. [Google Scholar] [CrossRef]
  12. Février, P.; Fandard, G. Thermal and mechanical tyre modelling for handling simulation. ATZ Worldw. 2008, 110, 26–31. [Google Scholar] [CrossRef]
  13. Pearson, M.; Blanco-Hague, O.; Pawlowski, R. TameTire: Introduction to the model. Tire Sci. Technol. 2016, 44, 102–119. [Google Scholar] [CrossRef]
  14. Jansen, S.T.; Verhoeff, L.; Cremers, R.; Schmeitz, A.J.; Besselink, I.J. MF-Swift simulation study using benchmark data. Veh. Syst. Dyn. 2005, 43 (Suppl. S1), 92–101. [Google Scholar] [CrossRef]
  15. Schmeitz, A.J.C.; Besselink, I.J.M.; Jansen, S.T.H. Tno mf-swift. Veh. Syst. Dyn. 2007, 45, 121–137. [Google Scholar] [CrossRef]
  16. Gipser, M. FTire-the tire simulation model for all applications related to vehicle dynamics. Veh. Syst. Dyn. 2007, 45, 139–151. [Google Scholar] [CrossRef]
  17. Gipser, M. FTire and puzzling tyre physics: Teacher, not student. Veh. Syst. Dyn. 2016, 54, 448–462. [Google Scholar] [CrossRef]
  18. Gallrein, A.; De Cuyper, J.; Dehandschutter, W.; Bäcker, M. Parameter identification for LMS CDTire. Veh. Syst. Dyn. 2005, 43 (Suppl. S1), 444–456. [Google Scholar] [CrossRef]
  19. Gallrein, A.; Bäcker, M. CDTire: A tire model for comfort and durability applications. Veh. Syst. Dyn. 2007, 45, 69–77. [Google Scholar] [CrossRef]
  20. Oertel, C.; Fandre, A. Ride comfort simulations and steps towards life time calculations: RMOD-K tyre model and ADAMS. In Proceedings of the International ADAMS Users’ Conference, Berlin, Germany, 17–19 November 1999. [Google Scholar]
  21. Oertel, C.; Fandre, A. Tire model RMOD-K 7 and misuse load cases. SAE Tech. Pap. 2009. [Google Scholar] [CrossRef]
  22. Gil, G.; Lee, S. In-Plane Bending and Shear Deformation of Belt Contributions on Tire Cornering Stiffness Characteristics. Tire Sci. Technol. 2021, 49, 276–297. [Google Scholar] [CrossRef]
  23. Clark, S.K. (Ed.) Mechanics of Pneumatic Tires, 1st ed.; US Government Printing Office: Washington, DC, USA, 1971. [Google Scholar]
  24. Nakajima, Y. Advanced Tire Mechanics, 1st ed.; Springer: Singapore, 2019. [Google Scholar]
  25. Sarkisov, P. Physical Understanding of Tire Transient Handling Behavior, 1st ed.; Cuvillier Verlag: Göttingen, Germany, 2019. [Google Scholar]
  26. Sarkisov, P.; Prokop, G.; Kubenz, J.; Popov, S. Physical understanding of transient generation of tire lateral force and aligning torque. Tire Sci. Technol. 2019, 47, 308–333. [Google Scholar] [CrossRef]
  27. Sakai, H. Tire Engineering; Guranpuri-Shuppan: Tokyo, Japan, 1987. (In Japanese) [Google Scholar]
  28. Sakai, H. Theoretical and experimental studies on the dynamic properties of tyres Part 1: Review of theories of rubber friction. Int. J. Veh. Des. 1981, 2, 78–110. [Google Scholar]
  29. Sakai, H. Theoretical and experimental studies on the dynamic properties of tyres. Part 2: Experimental investigation of rubber friction and deformation of a tyre. Int. J. Veh. Des. 1981, 2, 182–226. [Google Scholar]
  30. Sakai, H. Theoretical and experimental studies on the dynamic properties of tyres: Part 3: Calculation of the six components of force and moment of a tyre. Int. J. Veh. Des. 1981, 2, 335–372. [Google Scholar]
  31. Sakai, H. Theoretical and experimental studies on the dynamic properties of tyres. Part 4: Investigations of the influences of running conditions by calculation and experiment. Int. J. Veh. Des. 1982, 3, 333–375. [Google Scholar]
  32. Gwanghun, G.; Parviz, E.N. An Analytical Model of Pneumatic Tyres for Vehicle Dynamic Simulations. Part 1: Pure slips. Int. J. Veh. Des. 1990, 11, 589–618. [Google Scholar]
  33. Gim, G.; Kim, S. A semi-physical tire model for practical cornering simulations of vehicles. In IPC-9 Paper 971283; Bali, Indonesia, November 1997. [Google Scholar]
  34. Gim, G.; Nikravesh, P.E. A three dimensional tire model for steady-state simulations of vehicles. SAE Trans. 1993, 102, 150–159. [Google Scholar]
  35. Kabe, K.; Miyashita, N. A new analytical tire model for cornering simulation. Part I: Cornering power and self-aligning torque power. Tire Sci. Technol. 2006, 34, 84–99. [Google Scholar] [CrossRef]
  36. Miyashita, N.; Kabe, K. A new analytical tire model for cornering simulation. Part II: Cornering force and self-aligning torque. Tire Sci. Technol. 2006, 34, 100–118. [Google Scholar] [CrossRef]
  37. Guo, K.H. A Unified Tire Model for Braking Driving and Steering Simulation (No. 891198); SAE Technical Paper; SAE: Warrendale, PA, USA, 1989. [Google Scholar]
  38. Lu, D. Modeling of Tire Cornering Properties under Transient Vertical Load and Its Effect on Vehicle Handling. Ph.D. Thesis, Jilin University, Changchun, China, 2003. [Google Scholar]
  39. Guo, K.; Liu, Q. A generalized theoretical model of tire cornering properties in steady state condition. SAE Trans. 1997, 106, 469–475. [Google Scholar]
  40. Gruber, P.; Sharp, R.S.; Crocombe, A.D. Normal and shear forces in the contact patch of a braked racing tyre. Part 1: Results from a finite-element model. Veh. Syst. Dyn. 2012, 50, 323–337. [Google Scholar] [CrossRef]
  41. Oertel, C. RMOD-K Formula Documentation; Brandenburg University of Applied Sciences: Brandenburg an der Havel, Germany, 2011. [Google Scholar]
  42. Nakajima, Y.; Hidano, S. Theoretical Tire Model for Wear Progress of Tires with Tread Pattern Considering Two-Dimensional Contact Patch. Tire Sci. Technol. 2021. [Google Scholar] [CrossRef]
  43. Sui, J. Study on Longitudinal and Side Slip Characteristics of Tire Under Steady State Conditions. Ph.D. Thesis, Jilin University, Changchun, China, 1995. [Google Scholar]
  44. Liu, Q. Analysis, Modeling and Simulation of Unsteady Sideslip Characteristics of Tire. Ph.D. Thesis, Jilin University, Changchun, China, 1996. [Google Scholar]
  45. Guo, K.; Liu, Q. Modelling and simulation of non-steady state cornering properties and identification of structure parameters of tyres. Veh. Syst. Dyn. 1997, 27, 80–93. [Google Scholar] [CrossRef]
  46. Romano, L. Advanced Brush Tyre Modelling, 1st ed.; Springer: Cham, Switzerland, 2022. [Google Scholar]
  47. Romano, L.; Timpone, F.; Bruzelius, F.; Jacobson, B. Analytical results in transient brush tyre models: Theory for large camber angles and classic solutions with limited friction. Meccanica 2022, 57, 165–191. [Google Scholar] [CrossRef]
  48. Romano, L.; Bruzelius, F.; Hjort, M.; Jacobson, B. Development and analysis of the two-regime transient tyre model for combined slip. Veh. Syst. Dyn. 2022, 1–35. [Google Scholar] [CrossRef]
  49. Romano, L.; Sakhnevych, A.; Strano, S.; Timpone, F. A novel brush-model with flexible carcass for transient interactions. Meccanica 2019, 54, 1663–1679. [Google Scholar] [CrossRef]
  50. Romano, L.; Bruzelius, F.; Jacobson, B. Unsteady-state brush theory. Veh. Syst. Dyn. 2021, 59, 1643–1671. [Google Scholar] [CrossRef]
  51. Hoogh, J.D. Implementing Inflation Pressure and Velocity Effects into the Magic Formula Tyre Model. Master’s Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2005. [Google Scholar]
  52. Zegelaar, P.W.A. The Dynamic Response of Tyres to Brake Torque Variations and Road Unevennesses. Ph.D. Thesis, Delft University of Technology, Delft, The Netherlands, 1998. [Google Scholar]
  53. Maurice, J.P. Short Wavelength and Dynamic Tyre Behaviour under Lateral and Combined Slip Conditions. Ph.D. Thesis, Delft University of Technology, Delft, The Netherlands, 2000. [Google Scholar]
  54. Xu, N.; Yang, Y.; Guo, K. A Discrete Tire Model for Cornering Properties Considering Rubber Friction. Automot. Innov. 2020, 3, 133–146. [Google Scholar] [CrossRef]
  55. Bai, F. Research on Tire Model Based on the Analysis of the Contact Patch. Ph.D. Thesis, Jilin University, Changchun, China, 2016. [Google Scholar]
  56. Available online: https://math.ecnu.edu.cn/~jypan/Teaching/MatrixIter/lect03_Iter_s.pdf (accessed on 1 August 2022).
  57. Greenbaum, A.; Chartier, T.P. Numerical Methods: Design, Analysis, and Computer Implementation of Algorithms; Princeton University Press: Princeton, NJ, USA, 2012. [Google Scholar]
  58. Zha, Y.; Yang, J.; Zeng, J.; Tso, C.H.M.; Zeng, W.; Shi, L. Review of numerical solution of Richardson–Richards equation for variably saturated flow in soils. Wiley Interdiscip. Rev. Water 2019, 6, e1364. [Google Scholar] [CrossRef]
  59. Singh, J.; Ganbari, B.; Kumar, D.; Baleanu, D. Analysis of fractional model of guava for biological pest control with memory effect. J. Adv. Res. 2021, 32, 99–108. [Google Scholar] [CrossRef]
  60. Yamashita, H.; Jayakumar, P.; Sugiyama, H. Physics-based flexible tire model integrated with lugre tire friction for transient braking and cornering analysis. J. Comput. Nonlinear Dyn. 2016. [Google Scholar] [CrossRef]
  61. Huang, S.; Wu, H.; Guo, K.; Lu, D.; Lu, L. An in-plane dynamic tire model with real-time simulation capability. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022. [Google Scholar] [CrossRef]
Figure 1. Top view of tire contact area and its deformations ( u , v ) with respect to ( C , x , y , z ) .
Figure 1. Top view of tire contact area and its deformations ( u , v ) with respect to ( C , x , y , z ) .
Machines 10 00801 g001
Figure 2. Model of the beam on the elastic foundation. (a) The equilibrium of forces for the beam section; (b) Belt/carcass deformation applied F y .
Figure 2. Model of the beam on the elastic foundation. (a) The equilibrium of forces for the beam section; (b) Belt/carcass deformation applied F y .
Machines 10 00801 g002
Figure 3. The contact patch model results.
Figure 3. The contact patch model results.
Machines 10 00801 g003
Figure 4. Tread element deformation and its slip displacement vector relative to the road.
Figure 4. Tread element deformation and its slip displacement vector relative to the road.
Machines 10 00801 g004
Figure 5. Tread element stiffness ellipse.
Figure 5. Tread element stiffness ellipse.
Machines 10 00801 g005
Figure 6. Curves of lateral force and aligning moment with iteration times (a) Fy responses with iteration times; (b) Mz responses with iteration times.
Figure 6. Curves of lateral force and aligning moment with iteration times (a) Fy responses with iteration times; (b) Mz responses with iteration times.
Machines 10 00801 g006
Figure 7. Contact patch shape comparison between model and test.
Figure 7. Contact patch shape comparison between model and test.
Machines 10 00801 g007
Figure 8. Contact pressure comparison between model and test. (a) Contact pressure along x center line; (b) Contact pressure along y center line.
Figure 8. Contact pressure comparison between model and test. (a) Contact pressure along x center line; (b) Contact pressure along y center line.
Machines 10 00801 g008
Figure 9. Parking comparison between model and test.
Figure 9. Parking comparison between model and test.
Machines 10 00801 g009
Figure 10. Step side slip angle comparison between model and test. (a) Side slip angle = 1 deg; (b) Side slip angle = 4 deg.
Figure 10. Step side slip angle comparison between model and test. (a) Side slip angle = 1 deg; (b) Side slip angle = 4 deg.
Machines 10 00801 g010
Figure 11. Pure side slip comparison between model and test at 0.05 Hz at V = 10 km/h. (a) Side force Fy; (b) Aligning moment Mz.
Figure 11. Pure side slip comparison between model and test at 0.05 Hz at V = 10 km/h. (a) Side force Fy; (b) Aligning moment Mz.
Machines 10 00801 g011
Figure 12. Response of lateral force and aligning moment under side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Figure 12. Response of lateral force and aligning moment under side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Machines 10 00801 g012
Figure 13. Response of lateral force and aligning moment under sine side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Figure 13. Response of lateral force and aligning moment under sine side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Machines 10 00801 g013aMachines 10 00801 g013b
Figure 14. Response of lateral force and aligning moment under load variations at constant side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Figure 14. Response of lateral force and aligning moment under load variations at constant side slip angle input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Machines 10 00801 g014
Figure 15. Response of the lateral force and aligning moment under step turn slip input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Figure 15. Response of the lateral force and aligning moment under step turn slip input. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Machines 10 00801 g015aMachines 10 00801 g015b
Figure 16. Response of lateral force and aligning moment under sine longitudinal combined with a constant side slip. (a) Response of longitudinal force using the proposed method. (b) Response of longitudinal force using the reference method. (c) Response of lateral force using the proposed method. (d) Response of lateral force using the reference method. (e) Response of aligning moment using the proposed method. (f) Response of aligning moment using the reference method.
Figure 16. Response of lateral force and aligning moment under sine longitudinal combined with a constant side slip. (a) Response of longitudinal force using the proposed method. (b) Response of longitudinal force using the reference method. (c) Response of lateral force using the proposed method. (d) Response of lateral force using the reference method. (e) Response of aligning moment using the proposed method. (f) Response of aligning moment using the reference method.
Machines 10 00801 g016aMachines 10 00801 g016b
Figure 17. Response of lateral force and aligning moment under sine turn combined with a constant side slip. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Figure 17. Response of lateral force and aligning moment under sine turn combined with a constant side slip. (a) Response of lateral force using the proposed method. (b) Response of lateral force using the reference method. (c) Response of aligning moment using the proposed method. (d) Response of aligning moment using the reference method.
Machines 10 00801 g017
Figure 18. Responses of longitudinal force and lateral force aligning moments under load variations at constant longitudinal and side slip input. (a) Response of longitudinal force using the proposed method. (b) Response of longitudinal force using the reference method. (c) Response of lateral force using the proposed method. (d) Response of lateral force using the reference method. (e) Response of aligning moment using the proposed method. (f) Response of aligning moment using the reference method.
Figure 18. Responses of longitudinal force and lateral force aligning moments under load variations at constant longitudinal and side slip input. (a) Response of longitudinal force using the proposed method. (b) Response of longitudinal force using the reference method. (c) Response of lateral force using the proposed method. (d) Response of lateral force using the reference method. (e) Response of aligning moment using the proposed method. (f) Response of aligning moment using the reference method.
Machines 10 00801 g018
Figure 19. Lateral deformation of the tread element at different traveling distances under the step side slip angle input. (a) Response of lateral force and aligning moment. (b) Deformation of the tread element at the contact patch center line under different travel distance. (c) Deformation of the tread element at A position. (d) Deformation of the tread element at B position. (e) Deformation of the tread element at C position. (f) Deformation of the tread element at D position. (g) Deformation of the tread element at E position. (h) Deformation of the tread element at F position.
Figure 19. Lateral deformation of the tread element at different traveling distances under the step side slip angle input. (a) Response of lateral force and aligning moment. (b) Deformation of the tread element at the contact patch center line under different travel distance. (c) Deformation of the tread element at A position. (d) Deformation of the tread element at B position. (e) Deformation of the tread element at C position. (f) Deformation of the tread element at D position. (g) Deformation of the tread element at E position. (h) Deformation of the tread element at F position.
Machines 10 00801 g019aMachines 10 00801 g019b
Figure 20. Lateral deformation of the tread element at different traveling distances under the step turn slip input. (a) Response of lateral force and aligning moment. (b) Deformation of the tread element at the contact patch center line under different travel distance. (c) Deformation of the tread element at A position. (d) Deformation of the tread element at B position. (e) Deformation of the tread element at C position. (f) Deformation of the tread element at D position. (g) Deformation of the tread element at F position. (h) Deformation of the tread element at E position.
Figure 20. Lateral deformation of the tread element at different traveling distances under the step turn slip input. (a) Response of lateral force and aligning moment. (b) Deformation of the tread element at the contact patch center line under different travel distance. (c) Deformation of the tread element at A position. (d) Deformation of the tread element at B position. (e) Deformation of the tread element at C position. (f) Deformation of the tread element at D position. (g) Deformation of the tread element at F position. (h) Deformation of the tread element at E position.
Machines 10 00801 g020
Figure 21. Lateral deformation of tread element at different traveling distance under load variations at a constant side slip angle. (a) Response of vertical force and lateral force for proposed method. (b) Response of vertical force and lateral force for the reference method. (c) Tread deformation at different positions for the proposed method. (d) Tread deformation at different positions for the reference method.
Figure 21. Lateral deformation of tread element at different traveling distance under load variations at a constant side slip angle. (a) Response of vertical force and lateral force for proposed method. (b) Response of vertical force and lateral force for the reference method. (c) Tread deformation at different positions for the proposed method. (d) Tread deformation at different positions for the reference method.
Machines 10 00801 g021
Table 1. The iterative CPU times of different methods at different conditions.
Table 1. The iterative CPU times of different methods at different conditions.
Iterative MethodIterative CPU Time [s] and the Simulation Time Is Set to 3 s
Step Side Slip Sine Side Slip Load Variations at Constant Side Slip Sine Side Slip Angle at Constant Turn SlipSine Longitudinal Slip at Constant Side Slip
Method 135.1194.1069.40no convergenceno convergence
Method 233.9771.4659.99no convergenceno convergence
Method 355.34136.21141.71154.96280.46
Method 433.6165.0759.67no convergenceno convergence
Method 532.5066.2660.14no convergenceno convergence
Method 6no convergenceno convergenceno convergenceno convergenceno convergence
Table 2. The discrete tire model parameters.
Table 2. The discrete tire model parameters.
Parameters NameParameters ValueParameters NameParameters Value
k t x [ N / m 3 ] 1.3377 × 108 K c x 0 [ N / m ] 4.3735 × 105
k t y [ N / m 3 ] 1.0332 × 108 E I z [ N m 2 ] 1 × 103
n [ ] 2 k s [ N / m 2 ] 1.25 × 105
λ [ ] 0 ξ [ ] 0.2
Δ [ ] 0 N θ [ N m / r a d ] 1.2994 × 104
a F z 1 [ ] −0.05 μ [ ] 1.11
p F z 1 [ N / m ] 2.01 × 105 R 0 [ m ] 0.3465
p λ [ ] 20 n y [ ] 5.4
G 1 [ ] 3.64 R y [ m ] 0.145
G 2 [ ] −0.74Tol10
d x [ m m ] 2 d y [ m m ] 2
Table 3. The test information.
Table 3. The test information.
Test No.Test ConditionsInflation PressureRoad SpeedSteer AngleSteer CycleVertical Load
1Contact Pressure250 kPa000 s5415 N
2Parking Maneuver250 kPa015~−15°32 s (Slope)5386 N
3Sideslip Angle Step250 kPa0 → 10 kph−1°-5414 N
250 kPa0 → 10 kph−4°-5407 N
4Sine Sideslip Angle250 kPa10 kph−15~15°20 s (0.05 Hz)5406 N
Table 4. The simulation errors at pure side slip simulation.
Table 4. The simulation errors at pure side slip simulation.
Simulation ConditionsCompare ItemsError [%]
Step Side Slip Angle F y 0.28
M z 1.12
Sine Side Slip Angle F y 1.37
M z 2.18
Load Variations at Constant Side Slip Angle F y 5.86
M z 7.59
Table 5. The simulation errors at the step turn slip simulation.
Table 5. The simulation errors at the step turn slip simulation.
Simulation ConditionsCompare ItemsError [%]
Step Turn Slip Simulation F y 141.40
M z 142.20
Table 6. The simulation errors at combined slip simulation.
Table 6. The simulation errors at combined slip simulation.
Simulation ConditionsCompare ItemsError [%]
Step Turn Slip Simulation F x 2.54
F y 1.23
M z 3.67
Sine Turn Combined with a Constant Side Slip F y 61.01
M z 6.43
Load Variations at Constant Longitudinal and Side Slip F x 3.95
F y 8.49
M z 8.86
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Q.; Lu, D.; Ma, Y.; Xia, D. A Novel Numerical Method for Theoretical Tire Model Simulation. Machines 2022, 10, 801. https://doi.org/10.3390/machines10090801

AMA Style

Liu Q, Lu D, Ma Y, Xia D. A Novel Numerical Method for Theoretical Tire Model Simulation. Machines. 2022; 10(9):801. https://doi.org/10.3390/machines10090801

Chicago/Turabian Style

Liu, Qianjin, Dang Lu, Yao Ma, and Danhua Xia. 2022. "A Novel Numerical Method for Theoretical Tire Model Simulation" Machines 10, no. 9: 801. https://doi.org/10.3390/machines10090801

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop