Wheel Dynamics

Download as pdf or txt
Download as pdf or txt
You are on page 1of 8

Proceedings of the XII International Symposium on Dynamic Problems of Mechanics (DINAME 2007), P.S.Varoto and M.A.

Trindade (editors), ABCM, Ilhabela, SP, Brazil, February 26 - March 2, 2007

Wheel Dynamics
Georg Rill1
1

University of Applied Sciences Regensburg, Galgenberstr. 30, 93053 Regensburg, Germany

Abstract: The dynamics of the wheel rotation is mainly inuenced by the driving or the braking torque and the longitudinal tire force. Today different tire models are available. Structural tire models are very complex. Here, the dynamics of the wheel rotation is dominated by the complex tire model. Within handling models, the steady state tire forces and torques are generated as functions of the longitudinal and lateral slip. Depending on the slip denition the dynamics of a wheel depends now on the vehicle velocity or the angular velocity of the wheel. Here, the wheel dynamics will become more and more stiff if the vehicle slows down or the wheel is close to a locking situation. This causes problems in drive away and braking to stand still maneuvers. In this paper a simple quarter car model is used to study the wheel dynamics. Different model approaches and the performance of explicit and implicit ode solvers are investigated. Keywords: Vehicle Dynamics, Wheel, Tire, ode-Solver

INTRODUCTION
Today, numerical simulation of system dynamics is a standard in the design of cars and trucks. Usually the vehicles are modeled by Multibody Systems (Rauh, 2003 as well as Weinfurter et al., 2004). A typical model for a passenger car consists of several subsystems, Fig. 1. Complex vehicle models are used to investigate the ride comfort and the handling

Driver Drive train Chassis Engine Suspension Wheel/Tire

Steering System

Figure 1 Vehicle modeled by subsystems

performance as well as to generate load data for life time prediction analysis. Additionally to these off-line simulation tasks real-time applications, typically used for the design of vehicle control systems and the test of electronic control units, become more and more popular. In order to achieve real time performance either the model complexity has to be reduced (Pankiewicz, 2003) or sophisticated modeling techniques and appropriate ode-solvers (Rill, 2006c) have to be used. Among many other problems (Rill and Chucholowski, 2005 as well as Rill 2006a) the wheel dynamics becomes critical in particular in drive away and braking to stand still maneuvers. The dynamics of the wheel rotation is mainly inuenced by the driving or the braking torque and the longitudinal tire force. Usually the driving torque, if present, is transmitted via the driving shaft to the wheel. The torsional stiffness of the driving shaft provides a simple torque-based interface between the wheel rotation and a separate drive train model. This interface is moderately stiff and can therefore be mastered by standard ode-solvers. The dynamics of a braked wheel is quite delicate, because a wheel may lock in an instant. Here, either a soft braking torque model or an implicit ode solver is needed.

SIMPLE VEHICLE MODEL Equations of Motion


In order to focus on the dynamics of the wheel motion a quite simple vehicle model will be used, Fig. 2. The linear momentum for the chassis and the angular momentum for the wheel result in mv = Fx , (1)

Wheel Dynamics

v m Fx r T

Figure 2 Simple chassis/wheel model

= T r Fx ,

(2)

where m is the corresponding chassis mass, and r denote the inertia and the radius of the wheel. Within this simple model the wheel suspension and the tire deection was not taken into account. Finally, T describes the driving or braking torque.

Longitudinal Tire Force


Today, different tire models are available (Lugner and Pl ochl, 2006). Structural tire models are very complex. They are computer time consuming and they need a lot a data. Usually, they are used for stochastic vehicle vibrations occurring during rough road rides and causing strength-relevant component loads. Here, the dynamics of the wheel rotation is dominated by the complex tire model. Handling models are characterized by a useful compromise between user-friendliness, model-complexity and efciency in computation time on the one hand, and precision in representation on the other hand (Hirschberg et al., 2002). Here, the steady state tire forces and torques are generated as functions of the longitudinal and lateral slip. The longitudinal force Fx can be described as a function of the longitudinal slip sx , Fig. 3. The longitudinal
4000

Longitudinal force Fx [N]

3000 2000 1000 0 -1000 -2000 -3000 -4000 -40 -30 -20 -10 0 10 20 30 40

Longitudinal slip sx [%]

Figure 3 Longitudinal force characteristics of a typical passenger car tire

slip is dened by rv , r || where r serves in this simple approach also for the dynamic rolling radius. sx = (3)

Linearized Equations of Motion


In normal driving situations the longitudinal slip is very small and the wheel is close to the rolling condition v . r Hence, disturbances of a normal driving situation with v = v0 and = v0 /r may be described by v = v0 + and For small disturbances v v0 and r sx = = v

(4)

(5) (6)

v0 + . r v0 the longitudinal slip can be simplied to v0 + r v0 |v0 + r | v r |v0 | v .

(7)

G. Rill

But, for small slip values the steady state force characteristics can be approximated by a linear function Fx dF0 sx , where dF0 describes the initial inclination of the longitudinal tire characteristics Fx (sx ). Assuming a constant driving velocity v0 = const the equations of motion can be written as v = and 1 dF0 (r m v0 v) (9) (8)

= T r dF0 (r v0

v) .

(10)

This set of linear differential equations can now be arranged in matrix form dF0 1 dF0 1 v v v0 m v0 m = dF0 r2 dF0 r2 r r v0 v0 x x A where the disturbances in the vehicle velocity variables.

(11)

0 + , T

v and the circumferential velocity of the wheel r

were used as state

Eigenvalues
The solution of the homogenous state equation x = Ax is given by x(t ) = x0 e t , where the eigen-values and eigenvectors x0 are dened by (A I ) x0 = 0 . (12) Non-trivial solutions x0 = 0 are possible if det |A I | = 0 or dF0 1 v0 m dF0 r2 dF0 1 v0 v0 m dF0 r2 v0 =0 (13)

will hold. The resulting quadratic equation has the solution 1 = 0 and 2 = dF0 v0 r2 1 + m . (14)

The second eigenvalue and hence the dynamic of the system depends on the driving velocity. For small driving velocities the system becomes very stiff. Hence, numerical integration algorithms based on explicit formulas must reduce the step size with a dropping driving velocity in order to maintain the stability.

NUMERICAL INTEGRATION MatLab-Function


The equations of motion (1) and (2) and a simple nonlinear longitudinal force characteristics are coded in the following MatLab function
function xp=f(t,x) % % simple vehicle model including wheel dynamics % % data global mass Theta r dF0 Fxmax Torque % % states v=x(1); % vehicle velocity o=x(2); % wheel ang. vel. % % longitudinal slip (nonlinear) sx = (r*o-v)/abs(r*o); %

Wheel Dynamics % % longitudinal force (simple fx = dF0*sx; if fx > +Fxmax, fx=+Fxmax; if fx < -Fxmax, fx=-Fxmax; % % equations of motion vp = fx/mass; % op = (Torque-r*fx)/Theta; % % % state derivatives xp = [ vp; op ];

nonlinear function) end end

linear momentum angular momentum

as a system of rst order differential equations.

Vehicle Data
Typical data for a passenger car are provided in Tab. 1. Chassis mass (quarter car) Inertia of wheel Wheel radius Initial incl. long. force char. Maximum long. force Driving torque m r dF0 Fxmax T = = = = = = 400 kg 1.2 kgm2 0.3 m 100 000 N / 3200 N 100 Nm

Table 1 Typical data for passenger car

Explicit Euler Integration


The explicit Euler integration x(t + h) = x(t ) + h f (t , x) is numerically stable as long as the integration step size h and the eigenvalues of the dynamic system will satisfy | 1 + h | < 1 . Applied to Eq. (14) one gets h or |v0 | h dF0 2 r2 1 + m . (18) dF0 v0 r2 1 + m 2 (17) (16) (15)

The step size of a stable explicit Euler solution depends on the vehicle data and the driving speed. Hence, in real-time applications where a constant step size must be used driving manoeuver with v0 0 cannot be handled by an explicit Euler algorithm. Here, for an integration step size of h = 0.5 ms the explicit Euler algorithm is stable as long as |v0 | > 1.9375 m/s will hold. The Integration was started at t = 0 with v(t = 0) = v0 = 2 ms and (t = 0) = v0 /r. This describes a vehicle which is rolling backwards. The constant driving torque of T = 100 Nm will then accelerate the vehicle forward. The results are plotted in Fig. 4. As predicted the explicit Euler algorithm becomes numerically unstable if the absolute value of the velocity drops below a critical value. Here, the numerical instability produces non-physical oscillations in the circumferential wheel velocity at r 1.8 m/s. On the basis of the linearized equations the critical velocity was predicted in Eq. 19 with a slightly larger value. (19)

G. Rill

-1.6 -1.7 Velocity [m/s] -1.8 -1.9 -2.0 -2.1 -2.2 -2.3 -2.4 -2.5 0 0.2 0.4 0.6 0.8 time [s] 1.0

r v

Figure 4 Results of the explicit Euler algorithm

Comparison to Implicit Formula


A comparison of the MatLab solvers ode23 (explicit Runge-Kutta (2,3) pair of Bogacki and Shampine) with ode23s (a modied implicit Rosenbrock formula of order 2) is shown in Fig. 5. For both formulas the defaults error values of
100 Step Size [s] 10-1 10-2 10-3 10-4 10-5 -2.5 explicit (ode23) implicit (ode23s)

-2.0

-1.5

-1.0 -0.5 Vehicle velocity [m/s]

Figure 5 Step sizes for explicit and implicit solver

RelTol = 103 and AbsTol = 106 were applied. In order to maintain stability the explicit Runge-Kutta algorithm has to reduce the step size nearly proportional to the absolut value of the vehicle velocity. Here, the implicit algorithm is approximately 50 times faster than the explicit one. Even at small absolute velocities it runs with considerable large step sizes.

Modied Slip
For a locked wheel or at stand still when = 0 will hold, the longitudinal slip in Eq. (3) is no longer dened. A small modication in the slip denition rv , (20) sx = r || + vnum where a small but nite velocity vnum > 0 was added to denominator avoids this problem. Now, the explicit Euler algorithm is stable if h r2 1 |v0 | + vnum > dF0 + . (21) 2 m

Wheel Dynamics

will hold. Hence, by setting vnum = 2 m/s a stable transition from negative to positive velocities is possible here, Fig. 6. The simple explicit Euler algorithm can master drive away and stand still maneuvers. But, the modied slip denition
1.0 Velocity [m/s] 0.5 0 -0.5 -1.0 -1.5 -2.0 -2.5 0 0.5 1.0 1.5 2.0 2.5 3.0 time [s] 3.5 4.0

0.1 Longitudinal slip [-] 0.08 0.06


physical

r v

0.04 0.02
modified

0 0

2 Time [s]

Figure 6 Stable solution with explicit Euler algorithm for vnum = 2 m/s

slows down the wheel dynamics especially at low wheel angular velocities and may hence produce results with poor accuracy. Figure 6 shows that the longitudinal slip based on physics will tend to innity at r 0 whereas the modied slip stay constant thus generating a constant longitudinal force which correspond to the constant driving torque. However, in normal driving situations, where r || vnum holds, the difference between the physically based slip and the modied slip will hardly be noticeable. In general, an implicit algorithm should be used to integrate the wheel rotation. Then, even the dynamics of a braked wheel can be mastered without problems (Rill, 2006a).

DYNAMIC TIRE FORCES Standard Approach


According to van der Jagt (2000) the tire forces show a dynamic behavior which can be approximated by a rst order differential equation xD = Fx FxD , Tx F (22) where the time constant Tx is derived from the relaxation length Tx =
x x,

r ||

(23)

Using Eqs. (23), (8) and (3) the differential equation (22) reads as xD = dF0 r v FxD F r || r ||
x

(24)

which, nally results in xD = dF0 (r v) r || FxD . F


x x

(25)

For = 0 this equation is still dened; it simplies to xD = dF0 v . F


x

(26)

Hence, the longitudinal force FxD is increased or decreased according to the magnitude and sign of vehicle velocity v.

Wheel/Tire Dynamics
The dynamic tire force FxD instead of the steady state value Fx must be used in the equations of motion (1) and (2) now. According to Eqs. (4) to (8) the equations of motion can be linearized. One gets m v = FxD , = T r FxD , xD = dF0 (r F
x

(27) v) |v0 |
x

FxD ,

G. Rill

where the term r || in the denominator of Eq. (23) was replaced by its linearized value |v0 |. Arranged in matrix form they read as 1 0 v 0 v 0 m 2 r (28) + T , r = 0 0 r dF dF0 v0 0 xD F FxD 0
x x x

Eigenvalues
The eigenvalues of the corresponding system matrix A are now given by v0 i 2 x dF0
x

1 = 0

and 2/3 =

1 r2 + m

v0 2 x

(29)

For v0 = 0 the eigenvalues 2/3 are purely imaginary. Hence, at stand still the wheel and the longitudinal tire force will perform undamped motions. Figure 7 shows the eigenvalues calculated with the vehicle data given in Tab. 1 and a relaxation length of x = 0.7 m which is a typical value for a normal passenger car tire. The stability region of the explicit Euler algorithm is a circle with radius R = 1/h touching the origin from the left. Again, to achieve stable solutions with the explicit Euler algorithm a very small step size must be used.
Im() 100 v0=0 Stability region of explicit Euler algorithm v0=100 2 40 20 40 1/h 20 1 -200 -150 -100 -50 0 Re()

v0=100

3 v0=0

Figure 7 Eigenvalues of vehicle with dynamic tire force and stability region of the explicit Euler algorithm

Simulation Results
For a step size of h = 0.0005 s the solution becomes unstable if the absolute circumferential velocity drops below 1 m/s, Fig. 8. Within the Euler solution the dynamic tire force FxD was limited to a maximum value of Fxmax = 3200 N . However, the dynamics of the vehicle with dynamic tire forces can be handled quite well with any kind of implicit algorithm or higher order explicit formula, Fig. 8. This simple dynamic tire model covers the whole velocity range. An enhanced dynamic tire model can be found in Rill, 2006b. Here, braking to stand still and drive away manoeuvres can be simulated in good accuracy because here, no changes in the wheel/tire dynamics at low velocities by modications in the slip denition are necessary. At stand still the dynamic tire model automatically changes to a stick slip model.

Wheel Dynamics
Euler explicit: h=0.005
2 Velocity [m/s] 1 0 -1 -2 0 3000 2000 1000 0 -1000 -2000 -3000 0 1
r v

ode23 and ode23s


2 Velocity [m/s] 1 0 -1 -2 4 0 350 340 1 3 2 Time [s] 4

10

Step sizes [s]

10
r v

-1

implicit: ode23s

10

-2

3 2 Time [s]

explicit: ode23
10
-3

Fx [N]

Fx [N]

330 320 310


10

-4

3 2 Time [s]

300 0

3 2 Time [s]

10

-5

-2.5

-2

-1.5

-1 -0.5 Velocity [m/s]

0.5

Figure 8 Simulation results and step sizes of explicit and implicit algorithms

CONCLUSION
In this paper a simple quarter car model was used to analyze the wheel/tire dynamics. Using steady state tire models the wheel/tire dynamics strongly depends on the velocity of the vehicle. A slip modication which slows down the wheel/tire dynamics at low velocities is needed to overcome the singularity at stand still. Dynamic tire models operate in the whole velocity range. It was shown further that the explicit Euler algorithm is not suitable for solving the wheel/tire dynamics especially at low velocities. Here, implicit algorithms perform best. In real time applications where a constant step size is needed the implicit Euler algorithm has proofed its efciency, Rill, 2006b. If the overall vehicle model cannot be integrated by an implicit solver a co-simulation of the overall vehicle with the subsystem wheel/tire is recommended.

REFERENCES
Hirschberg, W, Rill, G., Weinfurter, H., User-Appropriate Tyre-Modeling for Vehicle Dynamics in Standard and Limit Situations, Vehicle System Dynamics, Vol. 38, No. 2, pp. 103-125. Lisse: Swets & Zeitlinger, 2002 van der Jagt, P., The Road to Virtual Vehicle Prototyping; new CAE-models for accelerated vehicle dynamics development, PhD-Thesis, Tech. Univ. Eindhoven, 2000 Lugner, P.; Pacejka, H. and Pl ochl, M., Recent advances in tyre models and testing procedures, Vehicle System Dynamics, Vol. 43, No. 67, pp. 413-436, 2005 Pankiewicz, E. and Rulka, W., From Off-Line to Real Time Simulations by Model Reduction and Modular Vehicle Modeling, in: Proceedings of the 19th Biennial Conference on Mechanical Vibration and Noise Chicago, Illinois, 2003 Rauh, J., Virtual Development of Ride and Handling Characteristics for Advanced Passenger Cars, Vehicle System Dynamics, Vol. 40, Nos. 1-3, pp. 135-155, 2003 Rill, G., A Modied Implicit Euler Algorithm for Solving Vehicle Dynamics Equation, Multibody System Dynamics, Vol. 15, No. 1, pp 1-24, 2006 Rill, G., First Order Tire Dynamics, Proceedings of the III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering, Lisbon, Portugal, 2006 Rill, G., Vehicle Modeling by Subsystems To appear in: Journal of Brazilian Society of Mechanical Sciences and Engineering - RBCM, 2006 Rill, G. and Chucholowski, C., Modeling Concepts for Modern Steering Systems, ECCOMAS Multibody Dynamics, Madrid, Spain, 2005 Weinfurter, H., Hirschberg, W., Hipp, E., Entwicklung einer St orgr oenkompensation f ur Nutzfahrzeuge mittels Steerby-Wire durch Simulation, in: Berechnung und Simulation im Fahrzeugbau, VDI-Bericht 1846, S. 923-941, VDI Verlag D usseldorf, 2004

RESPONSIBILITY NOTICE
The author is the only responsible for the printed material included in this paper.

You might also like