The present invention generally relates to vehicle dynamics control. More particularly, the invention relates to determining longitudinal vehicle velocity in a vehicle dynamics control system.
Accurate knowledge of longitudinal vehicle velocity can be used to improve many vehicle dynamics control systems. For example, the Anti-lock Braking System (ABS) controller can maximize the braking force and minimize the stopping distance. However, the vehicle's longitudinal velocity is usually not available in production vehicles. Wheel speeds are used instead in many situations. But wheel speeds do not represent the vehicle's longitudinal velocity when excessive tire slip exists.
Several methods have been proposed to determine longitudinal vehicle velocity. Certain methods use wheel speeds and longitudinal acceleration to determine longitudinal vehicle velocity and, therefore, need compensation for sensor bias, changes of rolling radii of wheels, and longitudinal slip. It has also been proposed to use the global positioning satellite system (GPS) in vehicle velocity calculations. However, this method suffers from low bandwidth because common GPS receivers have less than 10 Hz update rate and it is too low for vehicle dynamics control applications.
A vehicle includes a longitudinal accelerometer and a GPS. A method for determining longitudinal acceleration of the vehicle includes measuring a signal from the longitudinal accelerometer. Further, the method includes determining a sensitivity and bias of the longitudinal accelerometer using a measure of vehicle longitudinal velocity provided by the GPS, a measure of vehicle lateral velocity, and a measure of vehicle yaw rate. Finally, the method includes applying the sensitivity and bias to the measured signal from the longitudinal accelerometer. Vehicle longitudinal velocity may be estimated by integrating the longitudinal acceleration.
Furthermore, individual wheel speeds are determined by calculating longitudinal and lateral velocity at the vehicle center of gravity (CG) using OPS velocity and vehicle yaw rate; calculating individual wheel longitudinal and lateral velocities from the longitudinal and lateral velocity at the vehicle CG and wheel distribution relative to the vehicle CG; and calculating individual wheel speeds using the individual wheel longitudinal and lateral velocities and steering angles. Vehicle longitudinal velocity may be estimated by integrating the longitudinal acceleration or from individual wheel speeds.
The invention may take physical form in certain parts and arrangement of parts, the preferred embodiment of which will be described in detail and illustrated in the drawings incorporated hereinafter, wherein:
Referring to
In an embodiment of the invention, a method 30 for correcting longitudinal accelerometer and determining longitudinal vehicle velocity using GPS is shown in
If GPS velocity is available (54) and if ABS and TCS are not active (60), the following condition shown in Equation 1 is checked (64) to avoid situations where excessive slip exists at the wheels.
|ax|<ax,th (1)
wherein, ax is the longitudinal acceleration of the vehicle and ax,th is a vehicle dependent threshold for longitudinal acceleration.
If the absolute value of ax is less than ax,th (68), then the longitudinal velocities at all wheels are calculated using the GPS velocity. The longitudinal and lateral velocities at the center of gravity are initially calculated (56) from the GPS velocity by taking into account the GPS antenna position and yaw rate. Once the longitudinal and lateral velocities at the center of gravity are calculated (56), the longitudinal and lateral velocities at each wheel are calculated (84) from the longitudinal and lateral velocities at the center of gravity using Equation 2:
vx,fl=vx−rTf/2, vy,fl=vy+ra
vx,fr=vx+rTf/2, vy,fr=vy+ra
vx,rl=vx−rTr/2, vy,rl=vy−rb
vx,rr=vx+rTr/2, vy,rr=vy−rb (2)
The wheel speeds at each wheel are calculated (84) from the longitudinal and lateral velocities using Equation 3:
vws,fl=vx,fl cos δf+vy,fl sin δf
vws,fr=vx,fr cos δf+vy,fr sin δf
vws,rl=vx,rl
vws,rr=vx,rr (3)
wherein, vws,fl, vws,fr, vws,rl, and vws,rr represent the speeds of each respective wheel, 22-28.
The rolling radii of the wheels are calculated (86) from the following relationship:
vws,fl=rw,flww,fl
vws,fr=rw,frww,fr
vws,rl=rw,rlww,rl
vws,rr=rw,rrww,rr (4)
wherein rw,fl, rw,fr, rw,rl, and rw,rr are the rolling radii of each respective wheel, 22-28, and wherein ww,flww,fr, ww,rl, and ww,rr are the rotational speeds of each respective wheel, 22-28.
The rolling radii of the wheels are then updated (88) using either a statistical estimation or updating method such as, but not limited to moving average method or a recursive least square method as is known in the art. Statistical estimation methods are disclosed in “Applied Optimal Estimation” by Arthur Gelb, 1974, The MIT Press, Massachusetts Institute of Technology, Cambridge, Mass. 02142, which is herein incorporated by reference.
If the condition the absolute value of ax is not less than ax,th (66), a determination is made if ax is positive, i.e. greater than zero (69). If ax is positive (72) for a vehicle which has an un-driven wheel (78), the rolling radii of each un-driven wheel is updated (82) using calculations (80), (90) similar to the calculations (84), (86), using Equations 2, 3, and 4. If ax is equal to or less than zero, the rolling radii are not updated.
|r|<rth (5)
After storing the sensor information in the buffer (98), a determination is made whether GPS velocity data is available (112). Bias and sensitivity errors of the longitudinal accelerometer are estimated (118) as shown in more detail in
With reference now to
ax={dot over (v)}x−rvy (6)
wherein ax represents the acceleration of the vehicle along the vehicle's longitudinal axis (i.e. the longitudinal acceleration of the vehicle).
For the longitudinal accelerometer 18, the acceleration of the vehicle along the accelerometer's longitudinal axis, ax,sensor, can be modeled as shown in Equation 7:
wherein θv and θr represent the pitch angles from the vehicle suspension and road grade respectively as shown in
In addition, a determination of an effect of the sensitivity and bias on the acceleration measurement is made using the following Equations 8-9:
ax,m=saxax,sensor+bax (8)
wherein ax,m is the measured longitudinal acceleration adjusted for sensor sensitivity and offset, wherein sax and bax are the sensitivity and bias of the longitudinal accelerometer, respectively.
Equation 7 is substituted into Equation 8, and then, since both θr, and θv are small angles Equations 7 and 8 become Equation 9:
ax,m=sax{dot over (v)}x−saxrvy+saxgθr+saxgθvbax (9)
Since the GPS system provides both horizontal and vertical velocities of the vehicle, road grade, θr, is estimated (124) as shown in Equation 10:
wherein vV,GPS and vH,GPS are horizontal and vertical velocities of the vehicle based on GPS values.
Pitch angle, as shown in
wherein kθ is the pitch stiffness of the vehicle.
Substituting Equation 6 into Equation 11, Equation 11 then becomes Equation 12:
Substituting Equation 12 into Equation 9, the measured longitudinal acceleration, ax,m, can be modeled as Equation 13:
wherein {dot over (v)}x is the acceleration of the longitudinal accelerometer, and sax,e is the effective sensitivity of the longitudinal accelerometer, which includes the pitch angle effect on the accelerometer due to gravity and is defined in Equation 14:
Since every term in Equation 13 except the sensitivity and bias of the longitudinal accelerometer, sax,e and bax, is measured or estimated, sax,e and bax can be estimated (130) through use of an estimation method such as a least square method using information stored in the buffer 19.
To estimate sax,e and bax, ux,delta(t) is defined in Equation 15:
wherein ux,m(t) is defined in Equation 16:
ux,m(t)=vx(t)+∫Or(−rvy+gθr)dt (16)
The problem is then formulated as in the following matrix form shown in Equation 17:
Equation 17 can be rewritten as Equation 18:
zux=Hux{circumflex over (x)}us (18)
wherein,
From Equation (18), sax,e and bax are estimated using a least square method as shown in Equation 20:
{circumflex over (x)}ux=(HuxTHus)−1HTzux (20)
Four longitudinal velocities of the vehicle 12 are initially calculated (134) from the associated four wheel speeds using Equations 2-4 and the rolling radii are updated (82) or (88). A determination (136) is made to determine if a brake is on. When the vehicle is braking, the maximum of the four velocities represents the longitudinal speed of the vehicle (160). When the vehicle 12 is not braking (138), a determination is made whether ax exceeds the threshold ax,th, or if the TCS is on (142). If ax exceeds the threshold, ax,th, as described in Equation (1), or if the TCS is on (146), a determination is made whether an un-driven wheel exists (148). An average of the velocities from the un-driven wheels is calculated (152) and represents a vehicle longitudinal speed for a vehicle which has un-driven wheels. For a vehicle that does not have un-driven wheels (162), a minimum of the four driven wheel velocities is calculated (158) to represent the vehicle longitudinal speed for a vehicle. If either conditions ax>ax,th or TCS on are not met, an average of the velocities from all four wheel speeds is calculated (154) and represents the longitudinal speed of the vehicle.
The wheel speeds are calculated using the rolling radii, which are updated in accordance with (82) or (88), as described with reference to
|ax,m−ax,ws|<Δax,th (21)
wherein ax,ws is the longitudinal acceleration of the vehicle estimated using the wheel speeds, and Δas,th is a vehicle dependent threshold. When condition shown in Equation 21 is not met (148), the longitudinal velocity estimate (38) from the longitudinal accelerometer (38) as shown in
The invention has been described with specific reference to the exemplary embodiments and modifications thereto. Further modifications and alterations may occur to others upon reading and understanding the specification. It is intended to include all such modifications and alterations insofar as they come within the scope of the invention.
Number | Name | Date | Kind |
---|---|---|---|
5862511 | Croyle et al. | Jan 1999 | A |
6029111 | Croyle | Feb 2000 | A |
6308134 | Croyle et al. | Oct 2001 | B1 |
6671587 | Hrovat et al. | Dec 2003 | B2 |
6671595 | Lu et al. | Dec 2003 | B2 |
6847875 | Lu et al. | Jan 2005 | B2 |
6915193 | Lu et al. | Jul 2005 | B2 |
7096116 | Tanaka et al. | Aug 2006 | B2 |
7426437 | Breed et al. | Sep 2008 | B2 |
20030130775 | Lu et al. | Jul 2003 | A1 |
20030130778 | Busse | Jul 2003 | A1 |
20030149512 | Hrovat et al. | Aug 2003 | A1 |
20030191568 | Breed | Oct 2003 | A1 |
20040172183 | Lu et al. | Sep 2004 | A1 |
20040204798 | Imada et al. | Oct 2004 | A1 |
20050143916 | Kim et al. | Jun 2005 | A1 |
20070027583 | Tamir et al. | Feb 2007 | A1 |
20080086251 | Lu et al. | Apr 2008 | A1 |
20080177495 | Gold | Jul 2008 | A1 |
20080294342 | Hoshizaki et al. | Nov 2008 | A1 |
20080306687 | Ryu et al. | Dec 2008 | A1 |
20090037058 | Senneff et al. | Feb 2009 | A1 |
20100017128 | Zeng | Jan 2010 | A1 |
20100047744 | Miura | Feb 2010 | A1 |
20100060443 | Wakao et al. | Mar 2010 | A1 |
20100312473 | Hoshizaki | Dec 2010 | A1 |
20100332266 | Tamir et al. | Dec 2010 | A1 |
20110066376 | Takaoka | Mar 2011 | A1 |
Number | Date | Country |
---|---|---|
10247994 | Apr 2004 | DE |
102004003877 | Aug 2004 | DE |
10200404805 | Sep 2004 | DE |
2003-017751 | Jan 2003 | JP |
2003-067324 | Feb 2003 | JP |
2006-335883 | Dec 2006 | JP |
2008-212757 | Aug 2008 | JP |
2009-211997 | Sep 2009 | JP |
2003-97845 | Dec 2003 | KR |
PCTIL2004000610 | Jul 2004 | WO |
Number | Date | Country | |
---|---|---|---|
20080306687 A1 | Dec 2008 | US |