The single track model

Similar documents
Single-track models of an A-double heavy vehicle combination

Tutorial 1 - Drive fundamentals and DC motor characteristics

Detection of Critical Driving Situations using Phase Plane Method for Vehicle Lateral Dynamics Control by Rear Wheel Steering

Derivation of a Six Degrees-of-Freedom Ground-Vehicle Model for Automotive Applications

SOLUTION di x = y2 dm. rdv. m = a 2 bdx. = 2 3 rpab2. I x = 1 2 rp L0. b 4 a1 - x2 a 2 b. = 4 15 rpab4. Thus, I x = 2 5 mb2. Ans.

Simple Car Dynamics. Outline. Claude Lacoursière HPC2N/VRlab, Umeå Universitet, Sweden, May 18, 2005

UNIVERSITY OF BOLTON SCHOOL OF ENGINEERING

MECH 3140 Final Project

Chapter 10 Single Track Models

PHYSICS (B) v 2 r. v r

Problems. 66 km/h B km/h 30 A. v A. 1.5 ft

Low Complexity MPC Schemes for Integrated Vehicle Dynamics Control Problems

Proposal of Step Climbing of Wheeled Robot Using Slip Ratio Control

The University of Melbourne Engineering Mechanics

INTI INTERNATIONAL UNIVERSITY FOUNDATION IN SCIENCE (CFSI) PHY1203: GENERAL PHYSICS 1 FINAL EXAMINATION: SEPTEMBER 2012 SESSION

Step Climbing Control of Wheeled Robot Based on Slip Ratio Taking Account of Work Load Shift by Anti-Dive Force of Suspensions and Accerelation

CHAPTER 8: ROTATIONAL OF RIGID BODY PHYSICS. 1. Define Torque

SECTION A. 8 kn/m. C 3 m 3m

Simulation of the Stick-Slip Friction between Steering Shafts Using ADAMS/PRE

Jackknife stability of a tractor semi-trailer combination

Simulation of an articulated tractor-implement-trailer model under the influence of lateral disturbances

Control of Mobile Robots Prof. Luca Bascetta

Road Vehicle Dynamics

Modeling and Validation of a Complex Vehicle Dynamics Model for Real-time Applications

Team-Exercises for DGC 100 Modelica Course

Chapter 1: Kinematics of Particles

Hierarchical steering control for a front wheel drive automated car

16.07 Dynamics Final Exam

Estimation of Tire-Road Friction by Tire Rotational Vibration Model

1. Linear Motion. Table of Contents. 1.1 Linear Motion: Velocity Time Graphs (Multi Stage) 1.2 Linear Motion: Velocity Time Graphs (Up and Down)

1 MR SAMPLE EXAM 3 FALL 2013

CONTROL DESIGN FOR AN OVERACTUATED WHEELED MOBILE ROBOT. Jeroen Ploeg John P.M. Vissers Henk Nijmeijer

Time-Optimal Automobile Test Drives with Gear Shifts

The basic principle to be used in mechanical systems to derive a mathematical model is Newton s law,

On my honor, I have neither given nor received unauthorized aid on this examination.

Plane Motion of Rigid Bodies: Forces and Accelerations

ENGINEERING COUNCIL DYNAMICS OF MECHANICAL SYSTEMS D225 TUTORIAL 3 CENTRIPETAL FORCE

Phys101 Second Major-173 Zero Version Coordinator: Dr. M. Al-Kuhaili Thursday, August 02, 2018 Page: 1. = 159 kw

Mechatronics. MANE 4490 Fall 2002 Assignment # 1

5. Plane Kinetics of Rigid Bodies

2. Kinetic friction - The force that acts against an object s motion. - Occurs once static friction has been overcome and object is moving

Kinematics, Dynamics, and Vibrations FE Review Session. Dr. David Herrin March 27, 2012

Q2. A machine carries a 4.0 kg package from an initial position of d ˆ. = (2.0 m)j at t = 0 to a final position of d ˆ ˆ

Electric Vehicle Performance Power and Efficiency

Problem Set x Classical Mechanics, Fall 2016 Massachusetts Institute of Technology. 1. Moment of Inertia: Disc and Washer

Final Exam April 30, 2013

Design a SSV. Small solar vehicle. Case SSV part 1

Chapter 6 Dynamics I: Motion Along a Line

Phys101 Third Major-161 Zero Version Coordinator: Dr. Ayman S. El-Said Monday, December 19, 2016 Page: 1

SOLUTION 8 7. To hold lever: a+ M O = 0; F B (0.15) - 5 = 0; F B = N. Require = N N B = N 0.3. Lever,

Physics 107 HOMEWORK ASSIGNMENT #9b

We reserve the right to make changes in the course of technical development MAN Nutzfahrzeuge Aktiengesellschaft

Final Examination Thursday May Please initial the statement below to show that you have read it

is acting on a body of mass m = 3.0 kg and changes its velocity from an initial

UNIVERSITY OF SASKATCHEWAN GE MECHANICS III FINAL EXAM APRIL 18, 2011 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS

Investigation of Steering Feedback Control Strategies for Steer-by-Wire Concept

Drive train. Steering System. Figure 1 Vehicle modeled by subsystems

[1] In your answer, you should use appropriate technical terms, spelled correctly.... A satellite. Fig. 4.1

Statement: This paper will also be published during the 2017 AUTOREG conference.

Mechanics II. Which of the following relations among the forces W, k, N, and F must be true?

Lecture 6 Physics 106 Spring 2006

State-Estimator Design for the KTH Research Concept Vehicle

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS

Appendix W. Dynamic Models. W.2 4 Complex Mechanical Systems. Translational and Rotational Systems W.2.1

Supplementary Problems

Chapter 10 Practice Test

Bench Test of Minimum Time Autonomous Driving for Electric Vehicle Based on Optimization of Velocity Profile Considering Energy Constraint

Vehicle Planar Motion Stability Study for Tyres Working in Extremely Nonlinear Region

Terramechanics V MARYLAND U N I V E R S I T Y O F. Terramechanics V. ENAE 788X - Planetary Surface Robotics

PROBLEM Copyright McGraw-Hill Education. Permission required for reproduction or display. SOLUTION. ω = 29.6 rad/s. ω = = 36 3.

MCE 366 System Dynamics, Spring Problem Set 2. Solutions to Set 2

String tyre model for evaluating steering agility performance using tyre cornering force and lateral static characteristics

Improving EV Lateral Dynamics Control Using Infinity Norm Approach with Closed-form Solution

7.6 Journal Bearings

4) Vector = and vector = What is vector = +? A) B) C) D) E)

A. unchanged increased B. unchanged unchanged C. increased increased D. increased unchanged

Modeling of Vehicle Dynamics using Matrix-Vector Oriented Calculation in Matlab.

Note on Posted Slides. Motion Is Relative

Dept. of EEE, KUET, Sessional on EE 3202: Expt. # 1 2k15 Batch

On my honor, I have neither given nor received unauthorized aid on this examination.

Selection Calculations For Motorized Actuators

Name Student ID Score Last First. I = 2mR 2 /5 around the sphere s center of mass?

MODELLING AND CONTROL OF A VEHICLE WITH SINGLE-WHEEL CHASSIS ACTUATORS. Ralf Orend

Supplementary Information Microfluidic quadrupole and floating concentration gradient Mohammad A. Qasaimeh, Thomas Gervais, and David Juncker

EQUATIONS OF MOTION: GENERAL PLANE MOTION (Section 17.5) Today s Objectives: Students will be able to analyze the planar kinetics of a rigid body

AP Physics 1 Summer Assignment (2014)

Physics 111. Lecture 23 (Walker: 10.6, 11.1) Conservation of Energy in Rotation Torque March 30, Kinetic Energy of Rolling Object

Vehicle Dynamics Control for Rollover Mitigation

(a) During the first part of the motion, the displacement is x 1 = 40 km and the time interval is t 1 (30 km / h) (80 km) 40 km/h. t. (2.

Quiz Number 4 PHYSICS April 17, 2009

Last Time: Start Rotational Motion (now thru mid Nov) Basics: Angular Speed, Angular Acceleration

Physics 106 Common Exam 2: March 5, 2004

First Name: Last Name: Section: 1. March 26, 2008 Physics 207 EXAM 2

(35+70) 35 g (m 1+m 2)a=m1g a = 35 a= =3.27 g 105

A study on wheel force measurement using strain gauge equipped wheels

Closed-form Method to Evaluate Bike Braking Performance

PLANAR KINETIC EQUATIONS OF MOTION: TRANSLATION

Rigid Body Kinetics :: Virtual Work

Chapter 6. Force and motion II

Chapter 12: Rotation of Rigid Bodies. Center of Mass Moment of Inertia Torque Angular Momentum Rolling Statics

Transcription:

The single track model Dr. M. Gerdts Uniersität Bayreuth, SS 2003 Contents 1 Single track model 1 1.1 Geometry.................................... 1 1.2 Computation of slip angles........................... 2 1.3 Longitudinal tyre forces............................ 2 1.4 Lateral tyre forces................................ 4 1.5 Air resistance.................................. 4 1.6 Equations of motion I.............................. 4 1.7 Equations of motion II............................. 5 1.8 Constraints................................... 5 1.9 Parameters................................... 5 i

1 1 Single track model 1.1 Geometry l l h e SP α F Lx δ α ψ F Ly F s F u α h h F uh F sh Notation:, h elocity front/rear wheel δ steering angle elocity α, α h slip angle front/rear wheel β side slip angle ψ yaw angle F s, F sh lateral tyre forces F u, F uh longitudinal tyre forces l, l h distance from center of graity to front/rear wheel distance from center of graity to drag mount point e SP F Lx, F Ly air resistance in longitudinal and lateral direction m mass of car The steering angle δ is related to the steering wheel angle δ w by δ w = i L δ.

2 1 SINGLE TRACK MODEL 1.2 Computation of slip angles The slip angles are gien by ( ) l ψ sin β α = δ arctan, cos β ( ) lh ψ + sin β α h = arctan. cos β Explanation: Since the car body does not expand or shrink, the elocity components in the longitudinal direction of the car body hae to be equal: cos β = h cos α h = cos(δ α ). In the lateral direction the difference between the elocities is gien by the yaw angle elocity: h sin α h = l h ψ + sin β, sin(δ α ) = l ψ sin β. Combining these four equations yield the aboe formulas for the respectie slip angles. 1.3 Longitudinal tyre forces The car has rear wheel drie. The drier controls the braking force F B 0, the gear i {1, 2, 3, 4, 5} and the accelerator pedal position φ. The latter will result in the torque M wheel (φ, i) = i g (i) i t M mot (φ, i), where M mot (φ, i) = f 1 (φ) f 2 (w mot (i)) + (1 f 1 (φ))f 3 (w mot (i)) denotes the motor torque and w mot (i) = i g(i) i t R 1 1 S denotes the rotary frequency of the motor depending on the gear i and the longitudinal slip S. For conenience, the slip is neglected, e.g. S = 0. If the slip is not neglected, it is defined by 1, if R ϕ, R ϕ S = R ϕ 1, otherwise

1.3 Longitudinal tyre forces 3 ϕ denotes the rotary frequency of the wheel and is gien by the differential equation I R ϕ = F uh R. The functions f 1, f 2 and f 3 are gien by f 1 (φ) = 1 exp( 3φ), f 2 (w mot ) = 37.8 + 1.54 w mot 0.0019 wmot, 2 f 3 (w mot ) = 34.9 0.04775 w mot. braking force: The braking force is distributed on the front and rear wheels by the formulas F B = 2 3 F B, F Bh = 1 3 F B such that F B + F Bh = F B holds. rolling resistance: The rolling resistance force at the front and rear wheel, respectiely, is gien by F R = f R () F z, F Rh = f R () F zh where f R () = f R0 + f R1 100 + f R4 is the friction coefficient and F z = m l h g l + l h, F zh = m l g l + l h ( ) 4 ( in [km/h]), 100 denote the static tyre loads at the front and rear wheel, respectiely. longitudinal force front wheel: longitudinal force rear wheel: F u = F B F R. F uh = M wheel(φ, i) R F Bh F Rh

4 1 SINGLE TRACK MODEL 1.4 Lateral tyre forces The lateral tyre forces are functions of the respectie slip angles (and the tyre loads, which are constant in our model). A simple model is the AT-model (arcustangens-model): F s = c AT 1 arctan(c AT 2 α ), F sh = c AT 1 arctan(c AT 2 α h ). A famous model is the magic formula of Pacejka: F s = D sin (C arctan (B α E (B α arctan(b α )))), F sh = D h sin (C h arctan (B h α h E h (B h α h arctan(b h α h )))). The slope of F s at α = 0 is gien by B C D and similar for F sh. 1.5 Air resistance F Lx = 1 2 c w ρ A 2, F Ly = 1 2 c y ρ A 2 R Notation: c w ρ A c y R air drag coefficient air density effectie flow surface elocity lateral air drag coefficient lateral air elocity 1.6 Equations of motion I ẋ = cos(ψ β), ẏ = sin(ψ β), = 1 m [(F uh F Lx ) cos β + F u cos(δ + β) (F sh F Ly ) sin β F s sin(δ + β)], β = w z 1 m [(F uh F Lx ) sin β + F u sin(δ + β) + (F sh F Ly ) cos β + F s cos(δ + β)], ψ = w z, ẇ z = 1 I zz [F s l cos δ F sh l h F Ly e SP + F u l sin δ]

1.7 Equations of motion II 5 1.7 Equations of motion II ẍ = 1 m [(F uh F Lx + F u cos δ F s sin δ) cos ψ (F sh F Ly + F u sin δ + F s cos δ) sin ψ] = 1 m [(F uh F Lx ) cos ψ + F u cos(δ + ψ) F s sin(δ + ψ) (F sh F Ly ) sin ψ] ÿ = 1 m [(F uh F Lx + F u cos δ F s sin δ) sin ψ + (F sh F Ly + F u sin δ + F s cos δ) cos ψ] = 1 m [(F uh F Lx ) sin ψ + F u sin(δ + ψ) + F s cos(δ + ψ) + (F sh F Ly ) cos ψ] ψ = 1 I zz [F s l cos δ F sh l h F Ly e SP + F u l sin δ] side slip angle: absolute elocity: (ẏ ) β = ψ arctan ẋ = ẋ 2 + ẏ 2 1.8 Constraints The steering angle is restricted by δ 0.5263 [rad]. The steering angle elocity is restricted by δ 0.5 [rad/s]. The braking force F B is restricted by 0 F B 15000 [N]. The accelerator pedal position φ is restricted by 0 φ 1. 1.9 Parameters Car:

6 1 SINGLE TRACK MODEL m 1239 [kg] car mass g 9.81 [m/s 2 ] acceleration due to graity l 1.19016 [m] distance from center of graity to front wheel l h 1.37484 [m] distance from center of graity to rear wheel e SP 0.5 [m] distance from center of graity to drag mount point R 0.302 [m] wheel radius I zz 1752 [kgm 2 ] moment of inertia i L 21.1 steering wheel transmission I R 1.5 moment of inertia of wheel Drag: c w 0.3 air drag coefficient ρ 1.249512 [N/m 2 ] air density A 1.4378946874 [m 2 ] effectie flow surface c y 0.3 lateral air drag coefficient Gear shift: i g (1) 3.91 first gear i g (2) 2.002 second gear i g (3) 1.33 third gear i g (4) 1.0 fourth gear i g (5) 0.805 fifth gear i t 3.91 motor torque transmission Tyre: AT-model c AT 1 3317.99 tyre coefficient AT-model c AT 1h 2872.29 tyre coefficient AT-model c AT 2 30.0 tyre coefficient AT-model c AT 2h 30.0 tyre coefficient AT-model Tyre: Pacejka-model

1.9 Parameters 7 B 10.96 tyre coefficient Pacejka-model (stiffness factor) C 1.3 tyre coefficient Pacejka-model (shape factor) D 4560.40 tyre coefficient Pacejka-model (peak alue) E 0.5 tyre coefficient Pacejka-model (curature factor) B h 12.67 tyre coefficient Pacejka-model (stiffness factor) C h 1.3 tyre coefficient Pacejka-model (shape factor) D h 3947.81 tyre coefficient Pacejka-model (peak alue) E h 0.5 tyre coefficient Pacejka-model (curature factor) Rolling resistance: f R0 0.009 coefficient f R1 0.002 coefficient f R4 0.0003 coefficient