In a method for determining the vehicle longitudinal velocity in a vehicle, a wheel longitudinal force is determined in an estimation equation as a function of the wheel speeds and the vehicle longitudinal velocity is ascertained as a function of the wheel longitudinal force.

Patent
   8473174
Priority
Aug 31 2009
Filed
Aug 17 2010
Issued
Jun 25 2013
Expiry
Aug 17 2030
Assg.orig
Entity
Large
2
8
all paid
10. A regulating or control unit configured to ascertain wheel speed using a sensor, determine a wheel longitudinal force as a function of the wheel speed, and ascertain longitudinal velocity as a function of the wheel longitudinal force;
wherein the wheel longitudinal force (Fw,x)is ascertained as a function of the wheel contact force (Fw,z)a roadway coefficient of frictionroad),and a longitudinal slip coefficient of frictionw) which is a function of a wheel longitudinal slip(k):

Fw,x=max{μroadw{κ}}·Fw,z;
wherein, for equations describing wheel dynamic response

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}

ŷ=h{{circumflex over (x)}}
with

{circumflex over (x)}=[{circumflex over (ω)}w{circumflex over (μ)}road]T

u=[MpMbFw,zυx]T
an observer equation system of the form

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}−k{{circumflex over (x)},u}({circumflex over (y)}31 y),

ŷ=h{{circumflex over (x)}}
is established, the amplification vector (k{{circumflex over (x)},u}) containing two amplification factors k1, k2), which are one of constant or dependent on a state variable.
11. A driver assistance system in a vehicle having a regulating or control unit, the regulating or control unit configured to ascertain wheel speed using a sensor, determine a wheel longitudinal force as a function of the wheel speed, and ascertain longitudinal velocity as a function of the wheel longitudinal force;
wherein the wheel longitudinal force (Fw,x) is ascertained as a function of the wheel contact force (Fw,z),a roadway coefficient of frictionroad), and a longitudinal slip coefficient of frictionw)which is a function of a wheel longitudinal slip(k):

Fw,x=max{μroadw{κ}}·Fw,z;
wherein, for equations describing wheel dynamic response

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}

ŷ=h{{circumflex over (x)}}
with

{circumflex over (x)}=[{circumflex over (ω)}w{circumflex over (μ)}road]T

u=[MpMbFw,zυx]T

ŷ={circumflex over (ω)}w,
an observer equation system of the form

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}−k{{circumflex over (x)},u}({circumflex over (y)}−y),

ŷ=h{{circumflex over (x)}}
is established, the amplification vector (k{{circumflex over (x)},u}) containing two amplification factors (k1, k2), which are one of constant or dependent on a state variable.
1. A method for determining the vehicle longitudinal velocity in a vehicle, comprising:
ascertaining wheel speed using a sensor;
determining a wheel longitudinal force as a function of the wheel speed; and
ascertaining vehicle longitudinal velocity as a function of the wheel longitudinal force;
wherein the wheel longitudinal force (Fw,x) is ascertained as a function of the wheel contact force (Fw,z), a roadway coefficient of frictionroad), and a longitudinal slip coefficient of frictionw) which is a function of a wheel longitudinal slip (K):

Fw,x=max{μroadw{κ}}·Fw,z;
wherein, for equations describing wheel dynamic response

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}

ŷ=h{{circumflex over (x)}}
with

{circumflex over (x)}=[{circumflex over (ω)}w{circumflex over (μ)}road]T

u=[MpMbFw,zvx]T

ŷ={circumflex over (ω)}w,
an observer equation system of the form

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u }−k{{circumflex over (x)},u }({circumflex over (y)}−y),

ŷ=h{{circumflex over (x)}}
is established, the amplification vector (k{{circumflex over (x)},u}) containing two amplification factors k1, k2, which are one of constant or dependent on a state variable.
12. An electronic stability program system having a regulating or control unit configured to ascertain wheel speed using a sensor, determine a wheel longitudinal force as a function of the wheel speed, and ascertain longitudinal velocity as a function of the wheel longitudinal force, and having a sensor system including a yaw rate sensor, a lateral acceleration sensor, and one wheel speed probe per wheel;
wherein the wheel longitudinal force (Fw,x) is ascertained as a function of the wheel contact force (Fw,z), a roadway coefficient of frictionroad) and a longitudinal slip coefficient of frictionw), which is a function of a wheel longitudinal slip (k):

Fw,x=max{μroadw{κ}}·Fw,z;
wherein, for equations describing wheel dynamic response

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}

ŷ=h{{circumflex over (x)}}
with

{circumflex over (x)}=[{circumflex over (ω)}w{circumflex over (μ)}road]T

u=[MpMbFw,zvx]T

ŷ={circumflex over (ω)}w,
an observer equation system of the form

{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}−k{{circumflex over (x)},u}({circumflex over (y)}−y),

ŷ=h{{circumflex over (x)}}
is established, the amplification vector (k{{circumflex over (x)},u}) containing two amplification factors (k1, k2),which are one of constant or dependent on a state variable.
2. The method as recited in claim 1, wherein the wheel longitudinal slip (K) is a function of the wheel speeds (ωw), the tire radius (rw), and the vehicle longitudinal velocity (vx):
κ = { ω w r w - v x ω w r w for ω w r w > v x ω w r w - v x v x for ω w r w < v x ,
for driving and braking, respectively.
3. The method as recited in claim 1, wherein the longitudinal slip coefficient of frictionw) increases linearly with the wheel longitudinal slip (K) and remains at a constant value upon reaching a limiting value.
4. The method as recited in claim 1, wherein a wheel rotational movement is determined as a function of a wheel moment of inertia (Jw), a wheel drive torque (Mp), a wheel braking torque (Mb), and the wheel longitudinal force (Fw,x) via a differential equation:
ω w = 1 J w ( M p - M b - F w , x r w ) .
5. The method as recited in claim 1, wherein the roadway coefficient of frictionroad) is determined as a function of a wheel drive torque(Mp), a wheel braking torque (Mb), and the wheel contact force (Fw,z) via a differential equation:

μroadroadg{Mp,Mb,Fw,z}.
6. The method as recited in claim 1, wherein at least one amplification factor (k1, k2) of the amplification vector (k{{circumflex over (x)},u})

k{{circumflex over (x)},u}=[k1k2]T
is determined as a function of one of a the relative wheel velocity ({circumflex over (v)}rel)or a variable correlated therewith according to

k1=k1{{circumflex over (v)}rel}

k2 =k2{{circumflex over (v)}rel}
the relative wheel velocity({circumflex over (v)}rel) being determined as a function of a calculated wheel speed ({circumflex over (ω)}w),the wheel radius (rw), and a calculated vehicle longitudinal velocity ({circumflex over (v)}x) according to

{circumflex over (v)}rel={circumflex over (ω)}wrw−{circumflex over (v)}x.
7. The method as recited in claim 6, wherein the first amplification factor (k1) of the amplification vector (k{{circumflex over (x)},u}) is set to an at least approximately constant, positive value.
8. The method as recited in claim 7, wherein the second amplification factor (k2) of the amplification vector (k{{circumflex over (x)},u}) is set to an at least approximately constant, positive value in the case of a braked wheel and is set to an at least approximately constant, negative value in the case of a driven wheel, the second amplification factor (k2) being an at least approximately linear function of the relative wheel velocity ({circumflex over (v)}rel) in the transition from positive to negative values.
9. The method as recited in claim 1, wherein the calculated vehicle longitudinal velocity ({circumflex over (v)}x)is ascertained as a function of the calculated wheel longitudinal forces ({circumflex over (F)}w,x)on the vehicle wheels from the movement equation
v . ^ x = 1 m F ^ w , x .

The present invention relates to a method for determining the vehicle longitudinal velocity in a vehicle.

A method for determining a vehicle reference velocity is described in German Patent No. DE 102 59 272 A1, in which various calculation paths are followed on the basis of acceleration values ascertained by sensors, in order to obtain preliminary individual values for the vehicle speed. The individual values for the vehicle velocity are weighted, the desired vehicle reference velocity being calculated as the mean value of the weighted individual velocities.

Ascertaining weighted variables for the individual wheel velocity values on the basis of a sensor system of an electronic stability program (ESP), which includes a lateral acceleration sensor, a yaw rate sensor, and four wheel speed probes, and calculating a mean value, which represents a support variable, as a function of the weighted wheel velocity values, is described in German Patent No. DE 199 36 710 A1.

Both in German Patent No. DE 102 59 272 A1 and also in German Patent No. DE 199 36 710 A1, the velocity variable to be ascertained is calculated on the basis of the free rolling wheel velocities, which are obtained from mathematical models, which only have a restricted validity, however. In the case of more complex models, via which the vehicle velocity may be ascertained with higher quality, a nonlinear model or a nonlinear algorithm is typically used, which may only be analyzed with increased computing effort, however.

An object of the present invention is to determine the vehicle longitudinal velocity in a vehicle with high quality and reduced computing effort.

Using an example method according to the present invention, the vehicle longitudinal velocity in a vehicle may be determined, preferably solely on the basis of an ESP (electronic stability program) sensor system, which includes one rotational speed probe per vehicle wheel, a yaw rate sensor, and an acceleration sensor. The example method requires only a comparatively low computing effort and nonetheless allows the ascertainment of the vehicle longitudinal speed with comparatively high quality.

The foundation of the example method is the determination of wheel longitudinal forces individual to the wheel in an estimation equation as a function of the wheel speeds and subsequently ascertaining the vehicle longitudinal velocity as a function of the particular wheel longitudinal force. The wheel longitudinal force is a function of various variables, in particular the roadway coefficient of friction, a longitudinal slip coefficient of friction, which is caused by the longitudinal slip on the tires, and the wheel contact force. The wheel contact force is known per se; it is provided either as a stationary value, which is calculated from the vehicle mass, or is determined from estimation equations, which are stored in a regulating or control unit in the vehicle. The longitudinal slip coefficient of friction is ascertained as a function of the wheel speed and the vehicle longitudinal velocity; the roadway coefficient of friction, which is typically not detected by technical measurement, is calculated from a differential equation system as a function of wheel torques and forces. Furthermore, the wheel speeds are expediently also computed in addition to the measurement, the difference between measured and estimated variables being evaluated in an observer equation system. The driving and braking torques on the wheel may also be computed or determined from measured variables, for example, the engine drive torque or the brake pressure in the brake circuit.

The equations which describe the wheel dynamic response are expediently based on an observation system, in particular a nonlinear Luenberger observer, the observer equation in the observer part including an amplification vector having scalar amplification factors, which are either predefined as constant variables or are ascertained as a function of state variables of the vehicle. In a preferred embodiment, at least one of the amplification factors is a function of the relative wheel velocity, which is ascertained as the difference of the calculated wheel speed multiplied by the wheel radius and calculated vehicle velocity. For example, a first amplification factor is set to an at least approximately constant positive value, while in contrast the second amplification factor in the amplification vector is set to different values in the case of a braked wheel and in the case of a driven wheel. In the case of a braked wheel, the second amplification factor expediently assumes an approximately constant positive value; in the case of a driven wheel, it expediently assumes an at least approximately constant negative value. In the transition between positive and negative values, the amplification factor runs at least approximately linearly as a function of the relative wheel velocity.

These equations adequately describe the wheel dynamic response. The wheel longitudinal force, on which the calculation of the vehicle longitudinal velocity is based, may be determined from the equations of the wheel dynamic response. In this case, a differential equation is expediently established for the vehicle longitudinal velocity, which is a function of the sum of the individual wheel longitudinal forces per vehicle wheel, divided by the vehicle mass.

Overall, in this way a differential-algebraic equation system is obtained, which includes an observer, the equation system being able to be solved with high precision and comparatively little computing time. Accordingly, in the case of an implementation in a regulating or control unit in the vehicle, the vehicle longitudinal velocity may be ascertained rapidly and with high quality, fundamentally solely based on the ESP sensor system. The regulating or control unit in which the method is executed is expediently part of a driver assistance system, for example, an electronic stability program or an antilock braking system, or is assigned to such a driver assistance system.

FIG. 1 shows a diagram having a simplified view of the longitudinal slip coefficient of friction as a function of the wheel longitudinal slip.

FIG. 2 shows a diagram having a first amplification factor, which is used in an observer, as a function of the relative wheel velocity.

FIG. 3 shows a diagram having a second amplification factor, which is part of the observer, also as a function of the relative wheel velocity.

FIG. 4 shows a block diagram of the estimator of the vehicle longitudinal velocity.

The following exemplary embodiment is based on a vehicle having four wheel speed probes, for example, from a standard ESP sensor system, which additionally has a yaw rate sensor and a lateral acceleration sensor, so that the wheel speeds at each vehicle wheel, the yaw rate, and the lateral acceleration in the vehicle may be determined by sensors. A requirement for the performance of the example method is the information of the wheel speed probes, the information of the yaw rate sensor and of the lateral acceleration sensor optionally being able to be used to improve the precision. In addition, the drive and braking torques of individual wheels are also to be available as estimated values, the brake pressures in the wheel brake cylinders optionally being able to be measured individually in the brake system of the vehicle, whereby the estimation of the wheel torques is improved.

The relationships which describe the wheel dynamic response will be explained first:

Longitudinal force Fw,x which is effectively utilized in the contact area between the tire and the roadway is a function of various factors. For the estimation of longitudinal velocity vx, it may approximately be represented by
Fw,x=max{μroadw{κ}}·Fw,z.

In this case, μroad is the coefficient of friction caused by the roadway, μw{κ} is the coefficient of friction caused by the longitudinal slip, and Fw,z is the contact force of the wheel.

As shown in FIG. 1, longitudinal slip coefficient of friction μw initially increases linearly and with greater slope with wheel longitudinal slip (κ) and remains at an at least approximately constant value upon reaching a limiting value.

Longitudinal slip (κ) is given by

κ = { ω w r w - v x ω w r w for ω w r w > v x ω w r w - v x v x for ω w r w < v x
for driving and braking, respectively, where ωw designates the wheel speed and rw designates the wheel radius. The wheel movement may therefore be described as a function of wheel moment of inertia Jw, wheel drive torque Mp, wheel braking torque Mb, and wheel longitudinal force Fw,x by the following differential equation:

ω w = 1 J w ( M p - M b - F w , x r w ) .

The above-described equations form the foundation for the description of the wheel dynamic response, but are completed by the observation of roadway coefficient of friction μroad, which reflects the traction of the current roadway. This variable may not be detected by technical measurement in mass production vehicles, but it has a great influence on the wheel dynamic response and therefore also on the vehicle dynamic response. For these reasons, the unknown behavior of coefficient of friction μroad having a dynamic response as a function of the system excitation is modeled according to the following differential equation:
μroadroadg{Mp,Mb,Fw,z},
where g{Mp, Mb, Fw,z} designates a standardized algebraic function, which contains the wheel torques and forces.

Wheel speed ωw is used as output y of the model:
y=ωw.

The mathematical system description is introduced for the further observations in functional form:
{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u},
ŷ=h{{circumflex over (x)}},
with
{circumflex over (x)}=[{circumflex over (ω)}w{circumflex over (μ)}road]T, u=[MpMbFw,zvx]T={circumflex over (ω)}w.

A nonlinear Luenberger observer may be established for this system as follows:
{dot over ({circumflex over (x)})}=f{{circumflex over (x)},u}−k{{circumflex over (x)},u}({circumflex over (y)}−y),
ŷ=h{{circumflex over (x)}}.

Therein, k{{circumflex over (x)},u} designates an amplification vector of the observer, which has both amplification factors k1, k2:
k{{circumflex over (x)},u}=[k1k2]T,
where amplification factors k1, k2 are ascertained as a function of relative wheel velocity {circumflex over (υ)}rel according to
k1=k1{{circumflex over (v)}rel}
k2=k2{{circumflex over (v)}rel}.

Relative wheel velocity {circumflex over (v)}rel is determined as a function of calculated wheel speed {circumflex over (ω)}w, wheel radius rw, and calculated vehicle longitudinal velocity {circumflex over (v)}x according to
{circumflex over (v)}rel={circumflex over (ω)}wrw−{circumflex over (v)}x

As may be inferred from FIG. 2, first amplification factor k1 is set to an at least approximately constant, positive value, a linear or nonlinear dependency on relative wheel velocity {circumflex over (v)}rel optionally also coming into consideration.

The curve of second amplification factor k2 is shown in FIG. 3. Second amplification factor k2 is set to an at least approximately constant, positive value in the case of a braked wheel (left half of diagram) and is set to an at least approximately constant, negative value in the case of a driven wheel (right half of the diagram), second amplification factor being an at least approximately linear function, but with greater slope, of relative wheel velocity {circumflex over (v)}rel in the transition from the positive to the negative value.

Alternatively, the dependency of the amplification factors may also be related to another variable related to relative velocity {circumflex over (v)}rel, for example, longitudinal slip κ.

With the aid of the present model and the estimator, the wheel dynamic response and the unknown roadway coefficient of friction may be ascertained. Current estimated tire force {circumflex over (F)}w,z in the longitudinal direction may thus be calculated at any time.

Fundamentally, other methods may also be used for calculating amplification vector k. Comparable results may be achieved using an unscented Kalman filter (UKF) and using a linear Kalman filter (KF). However, the above-described Luenberger observer requires significantly less calculation time and fewer memory resources.

Estimated tire forces {circumflex over (F)}w,z of all wheels are used to ascertain longitudinal velocity {circumflex over (v)}z of the vehicle according to the movement equation

v . ^ x = 1 m F ^ w , x .

If this movement equation is combined with the equations describing the wheel dynamic response, the complete sequence for estimating the longitudinal velocity results, see FIG. 4. This concept may be used in the case of front-wheel drive, rear-wheel drive, and all-wheel drive.

As shown in FIG. 4, drive torque Mp, braking torque Mb, wheel speed ωw, and wheel contact force Fw,z for each wheel are entered in the calculation block as input variables. In addition, a previously ascertained value for vehicle longitudinal velocity {circumflex over (v)}x is supplied to the calculation block.

In the calculation block, the wheel dynamic response is calculated according to the above-described wheel equations, which are carried out for each vehicle wheel. Calculated wheel longitudinal force Fw,x for each vehicle wheel, which is used in a further calculation block for ascertaining vehicle longitudinal velocity {circumflex over (v)}x, is available at the output of the calculation block. Since this is a recursive algorithm, output-side vehicle longitudinal velocity {circumflex over (v)}x is supplied to the first calculation block as an input variable.

Antonov, Sergey, Kroeger, Patrick

Patent Priority Assignee Title
10246063, Jul 14 2016 GOODRICH CORPORATION Wheel reference balance estimator
9701289, Mar 25 2013 CONTINENTAL TEVES AG & CO OHG Method for determining a vehicle reference speed and vehicle controller having such a method
Patent Priority Assignee Title
6015192, Jul 18 1996 Nissan Motor Co., Ltd. System for estimating vehicle body speed and road surface friction coefficient
6611781, Aug 06 1999 Robert Bosch GmbH Method and device for determining a speed value
6659570, Nov 20 2000 Toyota Jidosha Kabushiki Kaisha Vehicle motion control device and method
20030163226,
DE10259272,
DE19936710,
FR2871889,
WO2007048947,
///
Executed onAssignorAssigneeConveyanceFrameReelDoc
Aug 17 2010Robert Bosch GmbH(assignment on the face of the patent)
Jan 09 2012ANTONOV, SERGEYRobert Bosch GmbHASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0278670513 pdf
Jan 17 2012KROEGER, PATRICKRobert Bosch GmbHASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0278670513 pdf
Date Maintenance Fee Events
Dec 19 2016M1551: Payment of Maintenance Fee, 4th Year, Large Entity.
Dec 18 2020M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
Dec 10 2024M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
Jun 25 20164 years fee payment window open
Dec 25 20166 months grace period start (w surcharge)
Jun 25 2017patent expiry (for year 4)
Jun 25 20192 years to revive unintentionally abandoned end. (for year 4)
Jun 25 20208 years fee payment window open
Dec 25 20206 months grace period start (w surcharge)
Jun 25 2021patent expiry (for year 8)
Jun 25 20232 years to revive unintentionally abandoned end. (for year 8)
Jun 25 202412 years fee payment window open
Dec 25 20246 months grace period start (w surcharge)
Jun 25 2025patent expiry (for year 12)
Jun 25 20272 years to revive unintentionally abandoned end. (for year 12)