# Patent

## Title
System for determining torque and tire forces using integrated sensing system

## Abstract

A vehicle (10) includes a control system (18) that is used to control a vehicle system. The control system determines an axle torque, and longitudinal forces at each tire in response to the axle torque. Lateral forces at each tire are determined in response to the longitudinal forces. The control system of the vehicle is determined in response to the longitudinal and lateral forces.

## Background
TECHNICAL FIELD
The present invention relates generally to a control apparatus for controlling a system of an automotive vehicle in response to sensed dynamic behavior, and more specifically, to a method and apparatus for determining tire torque and the tire forces between the tires and the road surface and controlling a vehicle system based on the tire torque and tire forces.
BACKGROUND
Various automotive vehicles have recently begun including vehicle dynamic control systems. Such vehicle dynamic control systems include yaw stability control systems, roll stability control systems, integrated vehicle dynamic control systems, etc. The ongoing goal of vehicle controls is to achieve a coordinated system level vehicle performances for ride, handling, safety and fuel economy.
With current advances in mechatronics, vehicle controls have increased opportunities for achieving performances, which were previously reserved for spacecraft and aircraft. For example, gyro sensors, previously only used in aircraft, have now been incorporated in various vehicle controls, and the anti-lock brake systems invented for airplanes are now standard automotive control systems. Current sensor technology generates ever-increasing opportunities for vehicle control. A typical vehicle control system senses 3-dimensional dynamic vehicle motions. For example, during yaw stability and roll stability controls, the control task involves three-dimensional motions along the vehicle roll, pitch, and yaw directions and along the vehicle longitudinal, lateral and vertical directions.
The coupling between different motion directions may not be as strong as in an aircraft or a spacecraft, however, they cannot be neglected in most maneuvers that involve vehicle rolling over or yawing. For example, the excessive steering of a vehicle will lead to excessive yaw and lateral motion, which further cause large rolling motion towards the outside of the turning. If the driver brakes the vehicle during the excessive steering, then the vehicle will also experience roll and pitch motions together with lateral and longitudinal accelerations. Hence, a successful vehicle dynamics control involves an accurate determination of the vehicle roll, pitch and relative yaw attitude (same as the so-called sideslip angle).
With the aforementioned vehicle attitude determination needs, a new vehicle sensing technology which contains an inertial measurement unit (IMU) and all the other sensors used in vehicle controls is desirable. This sensing system is called an Integrated Sensing System (short to ISS). IMUs have been used in inertial navigation system (INS) for aircraft and satellite for decades. Typically an INS system determines the attitude of a flight vehicle through the sensor signals from the IMU sensors. The IMU sensor set includes three gyros and three linear accelerometers. The INS contains an IMU and a processor unit to compute the navigation solutions necessary for navigation, attitude reference and various other data communication sources. At the same token, the ISS will also be used (but not limited) to vehicle attitude determination.
With the use of IMU sensor cluster and the other standard sensors equipped with a vehicle, accurate estimation of the vehicle operating states is possible. One of the important states are the forces and torques applied to the wheels, including tire longitudinal and lateral forces, the driving torques and braking torques applied to the wheels. Those torques and forces may be used to determine the intention of the driver, the road surface condition and to facilitate the vehicle dynamics controls like yaw stability control and roll stability control. A four-wheel model is used in the computation rather than a bicycle model in prior known systems. This results in a more accurate determination of the forces and torques. Existing computations of the tire lateral forces using a bicycle model are conducted along the vehicle body frame; hence normal loading will contaminate the computation. The bicycle model also cannot differentiate vehicle yaw motion due to longitudinal force deviation between the left and right sides and is inaccurate during brake intervention.
It would therefore be desirable to accurately estimate the tire lateral and longitudinal forces applied to the wheels, and the braking and driving torques applied to each wheel.
SUMMARY
The present invention provides a method and system to dynamically predict the tire torque and the forces between the tire and road surface by using the available sensors typically equipped in an integrated vehicle dynamics control system. The computation unit for this will be described below as the Force and Torque Estimation (FATE) unit. The output of the FATE unit may be used in dynamic control systems. For example, the rollover event happens when one side of the wheels experience extremely high lateral tire forces. The accurate determination of the lateral tire forces in the FATE unit can be used to assess the potential for a rollover, which will help cut unnecessary activations that could occur if only the lateral acceleration is used.
In one aspect of the invention, a method of controlling a vehicle having tires includes determining an axle torque, determining longitudinal forces at each tire in response to the axle torque, determining lateral forces at each tire in response to the longitudinal forces, and controlling a vehicle system in response to the longitudinal forces and the lateral forces.
In a further aspect of the invention, a method of controlling a vehicle includes determining a road frame longitudinal acceleration, determining a road frame lateral acceleration of the vehicle, determining a road frame yaw rate, determining an axle torque, determining a steered wheel angle, and determining tire longitudinal accelerations in response to the axle torque. The method further includes determining lateral and longitudinal tire forces in response to the road frame longitudinal accelerations, road frame lateral acceleration, the wheel speed derivatives, the drive torque from powertrain and the steered wheel angle. A vehicle system is controlled in response to the computed tire longitudinal and lateral forces.
One advantage of the invention is that no additional sensors in addition to the dynamic control system sensors are required for the system. Therefore, the tire torque and tire forces may be determined in a cost effective manner.
Other advantages and features of the present invention will become apparent when viewed in light of the detailed description of the preferred embodiment when taken in conjunction with the attached drawings and appended claims.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagrammatic view of a vehicle with variable vectors and coordinate frames according to the present invention.
FIG. 2 is a block diagram of a stability system according to the present invention.
FIG. 3 is a front view of an automotive vehicle illustrating various angles according to the present invention.
FIG. 4 is a side view of an automotive vehicle illustrating various variables thereon.
FIG. 5 is a top view of an automotive vehicle having the various forces acting thereon.
FIG. 6 is a side view of a wheel of an automotive vehicle having the wheel speed and torques acting thereon.
FIG. 7 is a logic flow diagram of a method for controlling a vehicle dynamic system in accordance with another embodiment of the present invention.



DESCRIPTION OF THE PREFERRED EMBODIMENT
In the following figures the same reference numerals will be used to identify the same components. The present invention is preferably used in conjunction with vehicle control systems, which include, but are not limited to a yaw stability control system, a roll stability control system, an integrated vehicle dynamics control system, or a total vehicle control system for achieving fuel economy and safety and other vehicle system level performances. The system is also described with respect to an inertial measurement unit (IMU) or integrated sensing system (ISS). However, discrete sensors may also be used to implement the present invention.
Referring to FIG. 1, an automotive vehicle 10 with a safety system of the present invention is illustrated with the various forces and moments thereon during a rollover condition. Vehicle 10 has front right (FR) and front left (FL) wheel/tires 12a and 12b and rear right (RR) wheel/tires 13a and rear left (RL) wheel/tires 13, respectively. The vehicle 10 may also have a number of different types of front steering systems 14a and rear steering systems 14b, including having each of the front and rear wheels configured with a respective controllable actuator, the front and rear wheels having a conventional type system in which both of the front wheels are controlled together and both of the rear wheels are controlled together, a system having conventional front steering and independently controllable rear steering for each of the wheels, or vice versa. Generally, the vehicle has a weight represented as Mg at the center of gravity of the vehicle, where g=9.8 m/s2 and M is the total mass of the vehicle.
As mentioned above, the system may also be used with safety systems including active/semi-active suspension systems, anti-roll bar, or airbags or other safety devices deployed or activated upon sensing predetermined dynamic conditions of the vehicle.
The sensing system 16 is coupled to a control system 18. The sensing system 16 may comprise many different sensors including the sensor set typically found in a roll stability control or a rollover control system (including lateral accelerometer, yaw rate sensor, steering angle sensor and wheel speed sensor which are equipped for a traditional yaw stability control system) together with a roll rate sensor, or in addition a pitch rate sensor and a vertical accelerometer. The various sensors will be further described below. The sensors may also be used by the control system in various determinations such as to determine a lifting event, determine normal forces on the wheels, determine a height and position of a mass, etc. The wheel speed sensors 20 are mounted at each corner of the vehicle and generate signals corresponding to the rotational speed of each wheel. The rest of the sensors of sensing system 16 may be mounted directly on the center of gravity of the vehicle body, along the directions x, y and z shown in FIG. 1. As those skilled in the art will recognize, the frame from b1,b2 and b3 is called a body frame 22, whose origin is located at the center of gravity of the car body, with the b1, corresponding to the x axis pointing forward, b2 corresponding to the y axis pointing off the driving side (to the left), and the b3 corresponding to the z axis pointing upward. The angular rates of the car body are denoted about their respective axes as ωx for the roll rate, ωy for the pitch rate and ωz for the yaw rate. Calculations may take place in an inertial frame 24 that may be derived from the body frame 22 as described below.
The roll, pitch and yaw angular rate sensors and the lateral, longitudinal, and vertical accelerometers may be mounted on the vehicle car body along the body frame directions b1,b2 and b3 which are the x-y-z axes of the sprung mass of the vehicle. The sensors may be in an IMU or ISS or discrete sensors. The sensors may be off the center of gravity but numerically translated to the center of gravity.
The longitudinal acceleration sensor is mounted on the car body located at the center of gravity, with its sensing direction along b1 axis, whose output is denoted as ax. The lateral acceleration sensor is mounted on the car body located at the center of gravity, with its sensing direction along b2 axis, whose output is denoted as ay.
The other frame used in the following discussion includes the road frame, as depicted in FIG. 1. The road frame system r1r2r3 is fixed on the driven road surface, where the r3 axis is along the average road normal direction computed from the normal directions of the four-tire/road contact patches.
In the following discussion, the Euler angles of the body frame b1b2b3 with respect to the road frame r1r2r3 are denoted as θxbr and θybr, which are also called the relative Euler angles (i.e., relative roll and relative pitch angles, respectively).
Referring now to FIG. 2, control system 18 is illustrated in further detail having a controller 26 used for receiving information from a number of sensors which may include a yaw rate sensor 28, a speed sensor 20, a lateral acceleration sensor 32, a vertical accelerometer sensor 33, a roll angular rate sensor 34, a longitudinal acceleration sensor 36, a pitch rate sensor 37, steering angle (of the wheels or actuator) position sensor 38, a suspension position (height) sensor 40. It should be noted that various combinations and sub-combinations of the sensors may be used.
The sensor cluster 16 may be within a housing 43, and includes a roll rate sensor 34 generating a roll rate signal, a pitch rate sensor 37, generating a pitch rate signal, a yaw rate sensor 38 generating a yaw rate signal, a longitudinal acceleration sensor 36 generating a longitudinal acceleration signal, a lateral acceleration sensor 32 generating a lateral acceleration signal, and a vertical acceleration sensor 33 generating a vertical acceleration sensor 56 generating a vertical acceleration signal.
If discrete or individual sensors are used lateral acceleration, roll orientation and speed may be obtained using a global positioning system (GPS). Based upon inputs from the sensors, controller 26 may control a safety device 44. Depending on the desired sensitivity of the system and various other factors, not all the sensors may be used in a commercial embodiment. Safety device 44 may control an airbag 45 or a steering actuator 46a-46d at one or more of the wheels 12a, 12b, 13a, 13b of the vehicle. Also, other vehicle components such as a suspension control 48 may be used to adjust the suspension to prevent rollover.
Roll angular rate sensor 34 and pitch rate sensor 37 may sense the roll condition or lifting of the vehicle based on sensing the height of one or more points on the vehicle relative to the road surface. Sensors that may be used to achieve this include but are not limited to a radar-based proximity sensor, a laser-based proximity sensor and a sonar-based proximity sensor. The roll rate sensor 34 may also use a combination of sensors such as proximity sensors to make a roll rate determination.
Roll rate sensor 34 and pitch rate sensor 37 may also sense the roll condition or lifting based on sensing the linear or rotational relative displacement or displacement velocity of one or more of the suspension chassis components. This may be in addition to or in combination with suspension position sensor 40. The suspension position sensor 40, roll rate sensor 34 and/or the pitch rate sensor 37 may include a linear height or travel sensor, a rotary height or travel sensor, a wheel speed sensor used to look for a change in velocity, a steering wheel position sensor, a steering wheel velocity sensor and a driver heading command input from an electronic component that may include steer by wire using a hand wheel or joy stick.
The yaw rate sensor 28, the roll rate sensor 34, the lateral acceleration sensor 32, and the longitudinal acceleration sensor 36 may be used together to determine that the wheel has lifted. Such sensors may be used to determine wheel lift or estimate normal loading associated with wheel lift.
The roll condition such as the relative roll angle of the vehicle may also be established by one or more of the following translational or rotational positions, velocities or accelerations of the vehicle including a roll gyro, the roll rate sensor 34, the yaw rate sensor 28, the lateral acceleration sensor 32, the vertical acceleration sensor 33, a vehicle longitudinal acceleration sensor 36, lateral or vertical speed sensor including a wheel-based speed sensor 20, a radar-based speed sensor, a sonar-based speed sensor, a laser-based speed sensor or an optical-based speed sensor.
Controller 26 may include sensing algorithms including but not limited to reference attitude and reference directional velocity determinations, global/relative attitude determination, directional velocity determination, sensor plausibility check, sensor signal conditioning, sensor error compensation, road profile, surface condition determination, and abnormal state monitoring.
The controller 26 includes various control units controlling the aforementioned sensing algorithms. These units may include: a reference signal unit 70 (reference signal generator (RSG)), which includes an attitude reference computation and a velocity reference computation, a road profile unit 72 (road profile determination unit (RPD)), an attitude unit or relative attitude determination unit 74 (RAD), a global attitude unit 76 (global attitude determination unit (GAD) and a directional unit 78 (directional velocity determination unit (DVD)), a sensor plausibility unit 80 (sensor plausibility check unit (SPC)), an abnormal state unit 82 (abnormal state monitoring unit (ASM)), a sensor signal compensating unit 84 (SSC), an estimation unit 86 (force and torque estimation unit (FATE)), a car body to fixed reference frame unit 88 (body to reference unit (B2R)), a normal loading unit 90 (normal loading determination unit (NLD)), and a vehicle parameter unit 92 (vehicle parameter determination unit (VPD)). Signals generated from any one of the aforementioned units are referred to prediction of vehicle operation states signals. The present invention lies primarily in the normal loading detection unit (NLD) 90.
Safety device 44 may control the position of the front right wheel actuator 46A, the front left wheel actuator 46B, the rear left wheel actuator 46C, and the right rear wheel actuator 46D. Although as described above, two or more of the actuators may be simultaneously controlled. For example, in a rack-and-pinion system, the two wheels coupled thereto are simultaneously controlled.
Safety device 44 may be coupled to a brake controller 60. Brake controller 60 controls the amount of brake torque at a front right brake 62a, front left brake 62b, rear left brake 62c and a rear right brake 62d. Other safety systems such as an antilock brake system 64, a yaw stability control system 66 and a traction control system 68 may also benefit from the knowledge of the roll gradient, roll rate parameter, roll acceleration coefficient, additional mass and position of the mass. This information may impact the control strategy such as modifying the brake force.
Speed sensor 20 may be one of a variety of speed sensors known to those skilled in the art. For example, a suitable speed sensor may include a sensor at every wheel that is averaged by controller 26. The controller may translate the wheel speeds into the speed of the vehicle. Yaw rate, steering angle, wheel speed and possibly a slip angle estimate at each wheel may be translated back to the speed of the vehicle at the center of gravity. Various other algorithms are known to those skilled in the art. Speed may also be obtained from a transmission sensor. For example, if speed is determined while speeding up or braking around a corner, the lowest or highest wheel speed may not be used because of its error. Also, a transmission sensor may be used to determine vehicle speed. The wheel speed sensors are mounted at wheel locations and are denoted as w1,w2,w3,w4 for left-front, right-front, left-rear and right-rear wheels, respectively.
The roll condition of a vehicle can be characterized by the relative roll angle between the vehicle body and the wheel axle and the wheel departure angle (between the wheel axle and the average road surface). Both the relative roll angle and the wheel departure angle may be calculated in relative roll angle estimation module (RAD 74) by using the roll rate and lateral acceleration sensor signals. If both the relative roll angle and the wheel departure angles are large enough, the vehicle may be in either single wheel lifting or double wheel lifting. On the other hand, if the magnitude of both angles is small enough, the wheels are likely all grounded. In case that both of them are not small and the double wheel lifting condition is detected or determined, the sum of those two angles will be used by the feedback control module to compute the desired actuation command for achieving rollover control performance.
The roll condition of a vehicle can be characterized by rolling radius-based wheel departure roll angle, which captures the angle between the wheel axle and the average road surface through the dynamic rolling radii of the left and right wheels when both of the wheels are grounded. Since the computation of the rolling radius is related to the wheel speed and the linear velocity of the wheel, such rolling-radius based wheel departure angle will assume abnormal values when there are large wheel slips. This happens when a wheel is lifted and there is torque applied to the wheel. Therefore, if this rolling radius-based wheel departure angle is increasing rapidly, the vehicle might have lifted wheels. Small magnitude of this angle indicates the wheels are all grounded.
The roll condition of the vehicle can be seen indirectly from the wheel longitudinal slip. If during a normal braking or driving torque the wheels at one side of the vehicle experience increased magnitude of slip, then the wheels of that side are losing longitudinal road torque. This implies that the wheels are either driven on a low mu surface or lifted up. The low mu surface condition and wheel-lifted-up condition can be further differentiated based on the chassis roll angle computation, i.e., in low mu surface, the chassis roll angle is usually very small. Hence, an accurate determination of chassis roll is desired.
The roll condition of the vehicle can be characterized by the normal loading sustained at each wheel. Theoretically, when a normal loading at a wheel decreases to zero, the wheel is no longer contacting the road surface. In this case a potential rollover is underway. Large magnitude of this loading indicates that the wheel is grounded. Normal loading is a function of the calculated chassis roll and pitch angles. Hence, an accurate determination of chassis roll and pitch angles is desired.
The roll condition can be identified by checking the actual road torques applied to the wheels and the road torques, which are needed to sustain the wheels when they are grounded. The actual road torques can be obtained through torque balancing for each wheel using wheel acceleration, driving torque and braking torque. If the wheel is contacting the road surface, the calculated actual road torques must match or be larger than the torques determined from the nonlinear torques calculated from the normal loading and the longitudinal slip at each wheel.
The roll condition of a vehicle can be characterized by the chassis roll angle itself, i.e., the relative roll angle θxr between the vehicle body and the wheel axle. If this chassis roll angle is increasing rapidly, the vehicle might be on the edge of wheel lifting or rollover. Small magnitude of this angle indicates the wheels are not lifted or are all grounded. Hence, an accurate determination of the chassis roll angle is beneficial for determining if the vehicle is in non-rollover events.
The roll condition of a vehicle can also be characterized by the roll angle between the wheel axle and the average road surface, this is called wheel departure angle. If the roll angle is increasing rapidly, the vehicle has lifted wheel or wheels and aggressive control action needs to be taken in order to prevent the vehicle from rolling over. Small magnitude of this angle indicates the wheels are not lifted.
The center of gravity C is also illustrated with nominal mass M. A roll axis is also illustrated at a distance D from the center of gravity. The symbol ay is the lateral acceleration.
The safety device may include control function/priority logic. As illustrated, the logic resides within safety device 44, but may be part of controller 26 and/or brake controller 60.
The FATE module 86 may be coupled to an engine and/or engine controller 94 to receive a signal corresponding to the engine torque generated by the engine. As will be described below, from this value, various other torques including the transmission axle and tire torque may be determined.
Referring now to FIG. 3, the relationship of the various angles of the vehicle 10 relative to the road surface 11 is illustrated. In the following, a reference road bank angle θbank is shown relative to the vehicle 10 on a road surface. The vehicle has a vehicle body 10a and wheel axle 10b. The wheel departure angle θwda is the angle between the wheel axle and the road. The relative roll angle θxr is the angle between the wheel axle 10b and the body 10a. The global roll angle θx is the angle between the horizontal plane (e.g., at sea level) and the vehicle body 10a. 
Another angle of importance is the linear bank angle. The linear bank angle is a bank angle that is calculated more frequently (perhaps in every loop) by subtracting the afore-mentioned relative roll angle from the calculated global roll angle. If all things were slowly changing without drift, errors or the like, the linear bank angle and reference road bank angle terms would be equivalent.
During an event causing the vehicle to roll, the vehicle body is subject to a roll moment due to the coupling of the lateral tire force and the lateral acceleration applied to the center of gravity of vehicle body. This roll moment causes suspension height variation, which in turn results in a vehicle relative roll angle (also called chassis roll angle or suspension roll angle). The relative roll angle is an important variable that is used as an input to the activation criteria and to construct the feedback brake pressure command, since it captures the relative roll between the vehicle body and the axle. The sum of such a chassis roll angle and the roll angle between wheel axle and the road surface (called wheel departure angle) provides the roll angle between the vehicle body and the average road surface, which is one of the important variables feeding back to the roll stability control module.
Referring now to FIG. 4, an automotive vehicle 10 is illustrated with various parameters illustrated thereon. The side view of automotive vehicle 10 is illustrated. A front suspension 82 and a rear suspension 82r is illustrated. The suspensions are coupled to the body at a respective suspension point 84f, 84r. The distance from the suspension point 84f to the center of the wheel is labeled zsh. The distance from the center of gravity CG to the front suspension is labeled as lf. The distance from the CG to the rear suspension point 84r is labeled as lr. The vertical distance between the center of gravity and the suspension point are labeled respectively as hf and hr A portion of the body axis b3 and the road axis r3 are illustrated. The angle therebetween is the relative pitch angle θyr. The rolling radius of the tire is labeled as zw.
Referring now to FIG. 5, a top view of vehicle 10. Lateral and longitudinal velocities of the center of gravity are denoted as vx and vy, a yaw angular rate is denoted as ωz, a front wheel steering angle is denoted as δs, lateral acceleration is represented by ay, and longitudinal acceleration is represented by ax.
Also illustrated is the front track width tf and rear track width tr. That is, the track widths are half the track widths from the center line or center of gravity of the vehicle to the center of the tire.
Referring now to FIG. 6, a wheel and tire assembly 12b is illustrated. The wheel speed is illustrated as ωi. The braking torque on the wheel is τbi. The drive torque of the wheel is illustrated as τdi. It should be noted that the subscript i represents an integer of which there are four typically in a vehicle.
Referring now to FIG. 7, the method of the present invention begins in start block 110. In step 102, signals from the various sensors are received. In step 104, the plausibility of each of the sensors is determined. In step 106, if the sensor signals are not plausible, step 108 is executed in which sensor fault processing logic is performed. In step 110, sensor signal compensation is performed. Sensor signal compensation may include compensating the sensor signals for various offsets.
In the following all the deviations are based on a four-wheel vehicle model (in contrast to a single track bicycle model). The four-wheel vehicle is driven on an average road surface by averaging the four tire contact patches. A yaw road frame r1r2r3 is used here to facilitate the computation. Such a yaw road frame has its longitudinal axle r1 and lateral axle r2 lie on the average road surface and its vertical axle r3 perpendicular to the average road surface, and it yaws and moves with vehicle. The attitude difference between the vehicle body fixed frame and the yaw road frame is calculated as the relative attitude from the Relative Attitude Determination (RAD) unit in ISS.
The rotational motion of the ith wheel in a four-wheel vehicle model shown in FIG. 6 can be described by the following equation

Iwi{dot over (ω)}i=τdi−τbi−FxiRi  (1)

where Fxi is the longitudinal force applied to the ith wheel from the road, which is called the longitudinal road force; Ri is the rolling radius of the ith wheel; τdi is the driving torque applied to the wheel and τbi is the braking torque applied to the wheel; ωi is the wheel rotational angular velocity measured through an ABS wheel speed sensor; Iwi is the moment of inertia of the wheel tire assembly around the wheel rotational axis. i=1, 2, 3, 4 correspond to the front left, the front right, the rear left and the rear right wheel respectively.

Notice that, the braking torque for the ith wheel can be directly calculated from the wheel pressure

τbi=τsipi  (2)

where pi is the brake pressure and τsi is the so-called specific torque for the ith brake, which depends on the following for disk brake systems:

the brake efficiency;
the friction coefficient between the brake lining and the disk;
the radius from the center of the brake lining to the center of the wheel;
the area of the brake piston; and
the number of pistons per wheel brake.
The drive torque τdi at the ith wheel is related to the engine torque τe. In step 112 the engine torque is determined. The following description presumes a front wheel drive vehicle with an open differential. Of course, these equations may be easily modified for a rear drive condition. In this case the following is true




  
    
      
        
          
            τ
            
              d
              ⁢
              
                
              
              ⁢
              1
            
          
          =
          
            
              τ
              
                d
                ⁢
                
                  
                
                ⁢
                2
              
            
            =
            
              τ
              fa
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            ω
            fa
          
          =
          
            
              
                w
                1
              
              
                2
                ⁢
                
                  R
                  1
                
              
            
            +
            
              
                w
                2
              
              
                2
                ⁢
                
                  R
                  2
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




where ωa denotes the rotational, speed of the drive axle. Through inertia relationships, the axle torque at the front axle may be determined from the engine torque in step 114. The engine torque may be determined from the torque mapping based on the engine speed and loading. Such mapping is typically done in engine development and is stored in a memory in the engine controller.

Ie is the engine rotational inertia, ωe the engine rotational speed, τtci the torque at the input of the torque converter,

Ie{dot over (ω)}e=τe−τtci  (4)

The τtco torque output of the torque converter is related to the torque input τtci 

τtco=ατtci  (5)

ωt is the transmission rotation speed, It the transmission inertia, rt the transmission ratio, then the output torque τt of the transmission satisfies

It{dot over (ω)}t=τtco−rtτt  (6)

ωd is the rotational speed of the drive axle, Id the rotational inertia of the drive axle,

Id{dot over (ω)}d=τt−rdaτd  (7)

Based on the above equations, the axle torque τa can be computed as in the following




  
    
      
        
          τ
          fa
        
        =
        
          
            
              
                c
                1
              
              ⁡
              
                (
                
                  
                    τ
                    e
                  
                  -
                  
                    
                      I
                      e
                    
                    ⁢
                    
                      ω
                      .
                    
                  
                
                )
              
            
            ⁢
            
              λ
              ⁡
              
                (
                
                  
                    
                      c
                      2
                    
                    ⁢
                    
                      ω
                      fa
                    
                  
                  
                    ω
                    e
                  
                
                )
              
            
          
          -
          
            
              c
              3
            
            ⁢
            
              I
              t
            
            ⁢
            
              
                ω
                .
              
              fa
            
          
          -
          
            
              c
              4
            
            ⁢
            
              I
              d
            
            ⁢
            
              ω
              fa
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where c1,c2,c3,c4 are 4 coefficients which are functions of the involved gear ratios, λ(●) is a nonlinear function.

Similar consideration may be conducted for a rear wheel drive, and four wheel drive vehicle. In the following discussion only the front wheel drive is considered.
Using the calculated axle torque for the front axle, the wheel motions for the four wheel rotational motions satisfying the following (subscript 1 and 2 are front wheels, 3 and 4 represent rear wheels):

Iwf{dot over (ω)}1=τfa−τb1−Fx1R1 

Iwf{dot over (ω)}2=τfa−τb2−Fx2R2 

Iwr{dot over (ω)}3=−τb1−Fx3R3 

Iwr{dot over (ω)}4=−τb2−Fx4R4  9)

For the four-wheel vehicle model shown in FIG. 5, the IMU sensor cluster is mounted on the vehicle body, i.e., along the body-fixed xyz frame, while the tire forces are along the yaw or road frame. Hence the body-fixed accelerations are converted to the road frame in step 116 by the following

ax-road=axs cos(θyr)+ays sin(θxr)sin(θyr)+azs sin(θxr)cos(θyr)

ay-road=ays cos(θxr)cos(θyr)+azs sin(θxr)  (10)

where axs is the longitudinal sensor output, ays is the lateral sensor output, θyr is a calculated relative pitch angle, θxr is a calculated relative roll angle.

The yaw rate signal around the axis perpendicular to the average road surface may be computed as in the following

ωz-road=−ωxs sin(θyr)+ωys sin(θxr)cos(θyr)+ωzs cos(θxr)cos(θyr),  (11)

where ωzs is the output of the yaw rate sensor, ωxs is the output of the roll rate sensor, ωys is the output of the pitch rate sensor.

The longitudinal road forces are also the forces to move the whole vehicle. Hence based on the four-wheel vehicle model in FIGS. 4 and 5, the following equations describe the total vehicle's longitudinal motion along the yaw road axle r1 

Mtax-road=(Fx1+Fx2)cos(δs)+Fx3+Fx4−(Fy1Fy2)sin(δs)  (12)

where Mt is the total vehicle mass; axs is the vehicle longitudinal acceleration measured by a sensor; δs is the steered wheel angle, which can be calculated through the driver's steering wheel (hand wheel) angle and the known steering gear ratio; θy is the vehicle body's pitch angle; Fy1, Fy2, Fy3 and Fy4 are the lateral forces applied to the front-left, front-right, rear-left and rear right wheels (tires) from the road, Fx1, Fx2, Fx3 and Fx4 and the longitudinal forces applied to the wheels/tires.

The lateral motion of the four-wheel vehicle model should obey the following equation of motion along the yaw road axle r2 and r3 

Iz{dot over (ω)}z-road=Mlong+Mlat 

Mtay-road=(Fx1+Fx2)sin(δs)+(Fy1+Fy2)cos(δs)+Fy3Fy4  (13)

where Mlong and Mlat are the yawing moments generated from the longitudinal forces and lateral forces respectively and they are along the yaw road axle r3 

Mlong=Fx1l cos(γ+δs)−Fx2l cos(γ−δs)+(Fx3−Fx4)tr 

Mlat=Fy1l sin(γ+δs)+Fy2l sin(γ−δs)−(Fy3+Fy4)lr  (14)

where

l=√{square root over (tf2+lf2)},γ=a tan(tf/lf).  (15)

Notice that for a front wheel steering vehicle, the rear wheel is kept in the same direction during travel. Therefore, the resultant effect of the lateral force on the two wheels of the rear axle Fy3 and Fy4 are the same as the effect of the their sum. At the same token, for the rear wheel steering case, the resultant effect of the lateral force on the two wheels of the rear axle Fy1 and Fy2 are the same as the effect of the their sum. Hence in front wheel steering case the two unknowns Fy3 and Fy4 are combined as one unknown variable

Fyr=Fy3+Fy4  (16)

and in the rear wheel steering case, the two unknowns Fy1, and Fy2 are combined as one unknown variable

Fyf=Fy1+Fy2  (17)

Therefore in the front wheel steering case, the following seven wheel forces are determined from seven equations: Fx1,Fx2,Fx3,Fx4,Fy1,Fy2 and Fyr; and in rear wheel driving case the following seven wheel forces are needed to be determined from seven equations: Fx1,Fx2,Fx3,Fx4,Fyf,Fy3 and Fy4. In the following discussion, the front wheel steering case is set forth. The rear wheel steering case may be similarly obtained.
Define a longitudinal force vector as




  
    
      
        
          F
          x
        
        =
        
          [
          
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
            
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
            
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    3
                  
                
              
            
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    4
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        18
        )
      
    
  




and a lateral force vector





  
    
      
        
          F
          y
        
        =
        
          [
          
            
              
                
                  F
                  
                    y
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
            
            
              
                
                  F
                  
                    y
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
            
            
              
                
                  F
                  
                    y
                    ⁢
                    
                      
                    
                    ⁢
                    r
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        19
        )
      
    
  




From Equation (9), the longitudinal forces for each wheel Fx may be determined in step 118 as in the following




  
    
      
        
          F
          x
        
        =
        
          [
          
            
              
                
                  
                    1
                    
                      R
                      1
                    
                  
                  ⁢
                  
                    (
                    
                      
                        τ
                        af
                      
                      -
                      
                        τ
                        
                          b
                          ⁢
                          
                            
                          
                          ⁢
                          1
                        
                      
                      -
                      
                        
                          I
                          wf
                        
                        ⁢
                        
                          
                            
                              w
                              .
                            
                            1
                          
                          
                            R
                            1
                          
                        
                      
                    
                    )
                  
                
              
            
            
              
                
                  
                    1
                    
                      R
                      2
                    
                  
                  ⁢
                  
                    (
                    
                      
                        τ
                        af
                      
                      -
                      
                        τ
                        
                          b
                          ⁢
                          
                            
                          
                          ⁢
                          2
                        
                      
                      -
                      
                        
                          I
                          wf
                        
                        ⁢
                        
                          
                            
                              w
                              .
                            
                            2
                          
                          
                            R
                            2
                          
                        
                      
                    
                    )
                  
                
              
            
            
              
                
                  
                    1
                    
                      R
                      3
                    
                  
                  ⁢
                  
                    (
                    
                      
                        -
                        
                          τ
                          
                            b
                            ⁢
                            
                              
                            
                            ⁢
                            3
                          
                        
                      
                      -
                      
                        
                          I
                          wr
                        
                        ⁢
                        
                          
                            
                              w
                              .
                            
                            3
                          
                          
                            R
                            3
                          
                        
                      
                    
                    )
                  
                
              
            
            
              
                
                  
                    1
                    
                      R
                      4
                    
                  
                  ⁢
                  
                    (
                    
                      
                        -
                        
                          τ
                          
                            b
                            ⁢
                            
                              
                            
                            ⁢
                            4
                          
                        
                      
                      -
                      
                        
                          I
                          wr
                        
                        ⁢
                        
                          
                            
                              w
                              .
                            
                            4
                          
                          
                            R
                            4
                          
                        
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        20
        )
      
    
  




From Equations (12) and (13), the matrix form of the following computation of Fy from Fx is determined




  
    
      
        
          F
          x
        
        =
        
          
            
              C
              
                -
                1
              
            
            ⁡
            
              (
              
                δ
                s
              
              )
            
          
          ⁢
          
            (
            
              
                M
                ⁡
                
                  [
                  
                    
                      
                        
                          a
                          
                            x
                            -
                            road
                          
                        
                      
                    
                    
                      
                        
                          
                            ω
                            .
                          
                          
                            z
                            -
                            road
                          
                        
                      
                    
                    
                      
                        
                          a
                          
                            y
                            -
                            road
                          
                        
                      
                    
                  
                  ]
                
              
              +
              
                
                  D
                  ⁡
                  
                    (
                    
                      δ
                      s
                    
                    )
                  
                
                ⁢
                
                  F
                  x
                
              
            
            )
          
        
      
    
    
      
        (
        21
        )
      
    
  




where the matrices C(δs),M and D(δs) are defined as in the following





  
    
      
        
          
            C
            ⁡
            
              (
              δ
              )
            
          
          =
          
            [
            
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  0
                
              
              
                
                  
                    l
                    ⁢
                    
                      
                    
                    ⁢
                    
                      sin
                      ⁡
                      
                        (
                        
                          
                            δ
                            s
                          
                          +
                          γ
                        
                        )
                      
                    
                  
                
                
                  
                    
                      -
                      l
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      sin
                      ⁡
                      
                        (
                        
                          
                            δ
                            s
                          
                          -
                          γ
                        
                        )
                      
                    
                  
                
                
                  
                    -
                    
                      l
                      r
                    
                  
                
              
              
                
                  
                    cos
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  
                    cos
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  1
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          M
          =
          
            [
            
              
                
                  
                    -
                    
                      M
                      t
                    
                  
                
                
                  0
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    I
                    z
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  0
                
                
                  
                    M
                    t
                  
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          
            D
            ⁡
            
              (
              δ
              )
            
          
          =
          
            [
            
              
                
                  
                    cos
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  
                    cos
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
                
                  1
                
                
                  1
                
              
              
                
                  
                    
                      -
                      l
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      cos
                      ⁡
                      
                        (
                        
                          
                            δ
                            s
                          
                          +
                          γ
                        
                        )
                      
                    
                  
                
                
                  
                    l
                    ⁢
                    
                      
                    
                    ⁢
                    
                      cos
                      ⁡
                      
                        (
                        
                          
                            δ
                            s
                          
                          -
                          γ
                        
                        )
                      
                    
                  
                
                
                  
                    -
                    
                      t
                      r
                    
                  
                
                
                  
                    t
                    r
                  
                
              
              
                
                  
                    -
                    
                      sin
                      ⁡
                      
                        (
                        
                          δ
                          s
                        
                        )
                      
                    
                  
                
                
                  
                    -
                    
                      sin
                      ⁡
                      
                        (
                        
                          δ
                          s
                        
                        )
                      
                    
                  
                
                
                  0
                
                
                  0
                
              
            
            ]
          
        
      
    
    
      
        (
        22
        )
      
    
  




For the magnitude of the steering angle |δs|≧δ with a small non-zero value threshold δ, (21) may be determined in step 120 and may be further expressed as in the following




  
    
      
        
          
            F
            
              y
              ⁢
              
                
              
              ⁢
              r
            
          
          =
          
            
              
                M
                t
              
              ⁡
              
                (
                
                  
                    a
                    
                      y
                      -
                      road
                    
                  
                  +
                  
                    
                      a
                      
                        x
                        -
                        road
                      
                    
                    
                      tan
                      ⁡
                      
                        (
                        
                          δ
                          s
                        
                        )
                      
                    
                  
                
                )
              
            
            -
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                +
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
              
                sin
                ⁡
                
                  (
                  
                    δ
                    2
                  
                  )
                
              
            
            -
            
              
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    3
                  
                
                +
                
                  F
                  
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    4
                  
                
              
              
                tan
                ⁡
                
                  (
                  
                    δ
                    s
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            F
            
              y
              ⁢
              
                
              
              ⁢
              1
            
          
          =
          
            
              -
              
                
                  
                    l
                    r
                  
                  ⁢
                  
                    F
                    
                      y
                      ⁢
                      
                        
                      
                      ⁢
                      r
                    
                  
                
                
                  
                    l
                    f
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
              
            
            +
            
              
                
                  
                    -
                    
                      M
                      t
                    
                  
                  ⁢
                  
                    a
                    
                      x
                      -
                      road
                    
                  
                  ⁢
                  l
                  ⁢
                  
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        γ
                        -
                        
                          δ
                          s
                        
                      
                      )
                    
                  
                
                -
                
                  
                    I
                    z
                  
                  ⁢
                  
                    
                      ω
                      .
                    
                    
                      z
                      -
                      road
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
              
              l
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            F
            
              y
              ⁢
              
                
              
              ⁢
              2
            
          
          =
          
            
              
                
                  l
                  r
                
                ⁢
                
                  F
                  
                    y
                    ⁢
                    
                      
                    
                    ⁢
                    r
                  
                
              
              
                
                  l
                  f
                
                ⁢
                
                  sin
                  ⁡
                  
                    (
                    
                      δ
                      s
                    
                    )
                  
                
              
            
            +
            
              
                
                  
                    M
                    t
                  
                  ⁢
                  
                    a
                    
                      x
                      -
                      road
                    
                  
                  ⁢
                  l
                  ⁢
                  
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        γ
                        +
                        
                          δ
                          s
                        
                      
                      )
                    
                  
                
                +
                
                  
                    I
                    z
                  
                  ⁢
                  
                    
                      ω
                      .
                    
                    
                      z
                      -
                      road
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        δ
                        s
                      
                      )
                    
                  
                
              
              l
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




When the steering angle |δs|<δ threshold, the above computation may no longer be valid. In this case, since Fy1 and Fy2 are along the same direction, their effects to the vehicle dynamics may not be distinguishable. In this case only the sum of Fy1 and F2, and the sum of Fy3 and Fy4 may be computed from the following relationship

Iz{dot over (ω)}z-road=Mlong+Mlat 

Mtay-road=(Fy1+Fy2)+(Fy3+Fy4)

Mlong=Fx1l cos(γ)−Fx2l cos(γ)+(Fx3−Fx4)tr 

Mlat=(Fy1+Fy2)l sin(γ)−(Fy3+Fy4)lr  (24)

That is, the following computation may be used




  
    
      
        
          
            
              F
              
                y
                ⁢
                
                  
                
                ⁢
                1
              
            
            +
            
              F
              
                y
                ⁢
                
                  
                
                ⁢
                2
              
            
          
          =
          
            
              
                [
                
                  
                    
                      
                        
                          
                            I
                            z
                          
                          ⁢
                          
                            
                              ω
                              .
                            
                            
                              z
                              -
                              road
                            
                          
                        
                        -
                        
                          
                            F
                            
                              x
                              ⁢
                              
                                
                              
                              ⁢
                              1
                            
                          
                          ⁢
                          l
                          ⁢
                          
                            
                          
                          ⁢
                          
                            cos
                            ⁡
                            
                              (
                              γ
                              )
                            
                          
                        
                        -
                      
                    
                  
                  
                    
                      
                        
                          
                            F
                            
                              x
                              ⁢
                              
                                
                              
                              ⁢
                              2
                            
                          
                          ⁢
                          l
                          ⁢
                          
                            
                          
                          ⁢
                          
                            cos
                            ⁡
                            
                              (
                              γ
                              )
                            
                          
                        
                        +
                        
                          
                            (
                            
                              
                                F
                                
                                  x
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  3
                                
                              
                              -
                              
                                F
                                
                                  x
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  4
                                
                              
                            
                            )
                          
                          ⁢
                          
                            t
                            r
                          
                        
                      
                    
                  
                
                ]
              
              +
              
                
                  M
                  t
                
                ⁢
                
                  a
                  
                    y
                    -
                    road
                  
                
                ⁢
                
                  l
                  r
                
              
            
            
              
                l
                r
              
              +
              
                l
                ⁢
                
                  
                
                ⁢
                
                  sin
                  ⁡
                  
                    (
                    γ
                    )
                  
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              F
              
                y
                ⁢
                
                  
                
                ⁢
                3
              
            
            +
            
              F
              
                y
                ⁢
                
                  
                
                ⁢
                4
              
            
          
          =
          
            
              
                -
                
                  [
                  
                    
                      
                        
                          
                            
                              I
                              z
                            
                            ⁢
                            
                              
                                ω
                                .
                              
                              
                                z
                                -
                                road
                              
                            
                          
                          -
                          
                            
                              F
                              
                                x
                                ⁢
                                
                                  
                                
                                ⁢
                                1
                              
                            
                            ⁢
                            l
                            ⁢
                            
                              
                            
                            ⁢
                            
                              cos
                              ⁡
                              
                                (
                                γ
                                )
                              
                            
                          
                          -
                        
                      
                    
                    
                      
                        
                          
                            
                              F
                              
                                x
                                ⁢
                                
                                  
                                
                                ⁢
                                2
                              
                            
                            ⁢
                            l
                            ⁢
                            
                              
                            
                            ⁢
                            
                              cos
                              ⁡
                              
                                (
                                γ
                                )
                              
                            
                          
                          +
                          
                            
                              (
                              
                                
                                  F
                                  
                                    x
                                    ⁢
                                    
                                      
                                    
                                    ⁢
                                    3
                                  
                                
                                -
                                
                                  F
                                  
                                    x
                                    ⁢
                                    
                                      
                                    
                                    ⁢
                                    4
                                  
                                
                              
                              )
                            
                            ⁢
                            
                              t
                              r
                            
                          
                        
                      
                    
                  
                  ]
                
              
              +
              
                
                  M
                  t
                
                ⁢
                
                  a
                  
                    y
                    -
                    road
                  
                
                ⁢
                l
                ⁢
                
                  
                
                ⁢
                
                  sin
                  ⁡
                  
                    (
                    γ
                    )
                  
                
              
            
            
              
                l
                r
              
              +
              
                l
                ⁢
                
                  
                
                ⁢
                
                  sin
                  ⁡
                  
                    (
                    γ
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




In step 122, various systems of the vehicle may be controlled in response to the axle torque, the longitudinal forces and the lateral forces. Such systems may include traction control systems, yaw control systems, rollover control systems, suspension systems, engine systems, various vehicle safety systems including deployable safety systems such as airbags and rollover bars. Various portions of the vehicle suspension may also be controlled such as an anti-roll bar.
While particular embodiments of the invention have been shown and described, numerous variations and alternate embodiments will occur to those skilled in the art. Accordingly, it is intended that the invention be limited only in terms of the appended claims.

## Claims

1. A method of controlling a vehicle having tires comprising:
determining an axle torque;
determining longitudinal forces at each tire in response to the axle torque;
determining a steering wheel angle;
determining a vehicle velocity;
determining a plurality of lateral forces;
when the steering wheel angle is at or exceeds a predetermined steering wheel angle threshold value, said determining of said plurality of lateral forces comprises determining lateral forces for each individual front tire and a lateral tire force for a rear tire in response to the longitudinal forces and the steering wheel angle;
when the steering wheel angle is below the predetermined steering wheel angle threshold value, said determining of said plurality of lateral forces comprises determining a sum of lateral forces for the front tires and a sum of lateral forces for the rear tires in response to the longitudinal forces and the steering wheel angle; and
controlling lateral motion of a vehicle system in response to the plurality of lateral forces.

2. A method as recited in claim 1 wherein determining said plurality of lateral forces comprises determining said plurality of lateral forces in response to the longitudinal forces and sensor signals used in a vehicle stability control system.

3. A method as recited in claim 1 wherein determining said plurality of lateral forces comprises determining said plurality of lateral forces in response to the longitudinal forces and a road frame longitudinal acceleration.

4. A method as recited in claim 3 further comprising determining the road frame longitudinal acceleration in response to a sensor lateral acceleration, a sensor vertical acceleration, a sensor longitudinal acceleration, a relative pitch angle and a relative roll angle.

5. A method as recited in claim 1 wherein determining said plurality of lateral forces comprises determining said plurality of lateral forces in response to the longitudinal force, a road frame longitudinal acceleration, a road frame lateral acceleration, a road frame yaw rate, a steered wheel angle.

6. A method as recited in claim 5 further comprising determining the road frame yaw rate in response to a sensor roll rate, a sensor yaw rate, a sensor pitch rate, a relative pitch angle and a relative roll angle.

7. A method as recited in claim 5 further comprising determining the road frame lateral acceleration in response to a sensor lateral acceleration, a sensor vertical acceleration, a relative pitch angle, and a relative roll angle.

8. A method of controlling a vehicle as recited in claim 1 wherein the step of controlling lateral motion of a vehicle is related to a vehicle safety system.

9. A method of controlling a vehicle as recited in claim 1 wherein the step of controlling lateral motion of a vehicle is related to a vehicle suspension.

10. A method of controlling a vehicle as recited in claim 1 wherein the step of controlling lateral motion of a vehicle is related to a dynamic control system.

11. A method of controlling a vehicle as recited in claim 10 wherein the dynamic control system comprises a roll stability control system.

12. A method of controlling a vehicle having a rear axle, said method comprising:
determining a road frame longitudinal acceleration;
determining a road frame lateral acceleration of the vehicle;
determining road frame yaw rate;
determining an axle torque in response to wheel speeds at the drive axle and an engine speed;
determining a steered wheel angle;
determining tire longitudinal forces in response to axle torque;
determining a plurality of tire lateral forces;
when the steered wheel angle is at or exceeds a predetermined steered wheel angle threshold value, said determining of said plurality of tire lateral forces comprises determining front individual tire lateral forces and a sum of tire lateral forces at the rear axle in response to the tire longitudinal forces, the road frame longitudinal acceleration, the road frame lateral acceleration and the steered wheel angle;
when the steered wheel angle is below the predetermined steered wheel angle threshold value, said determining of said plurality of tire lateral forces comprises determining a sum of tire lateral forces at the front axle and a sum of tire lateral forces at the rear axle in response to the tire longitudinal forces, the road frame longitudinal acceleration and the road frame lateral acceleration; and
controlling lateral motion of a vehicle system in response to the tire plurality of tire lateral forces.

13. A method as recited in claim 12 further comprising determining the road frame longitudinal acceleration in response to a sensor lateral acceleration, a sensor vertical acceleration, a sensor longitudinal acceleration, a relative pitch angle and a relative roll angle.

14. A method as recited in claim 12 further comprising determining the road frame yaw rate in response to a sensor roll rate, a sensor yaw rate, a sensor pitch rate, a relative pitch angle and a relative roll angle.

15. A method as recited in claim 12 further comprising determining the road frame lateral acceleration in response to a sensor lateral acceleration, a sensor vertical acceleration, a relative pitch angle and a relative roll angle.

16. A method of controlling a vehicle as recited in claim 12 wherein controlling lateral motion of a vehicle comprises controlling at least one of a vehicle safety system, a vehicle suspension, a dynamic control system, a roll stability control system, a yaw stability control system and a vehicle suspension.