Method & system for improving the driving behavior of a vehicle

0Associated
Cases 
0Associated
Defendants 
0Accused
Products 
59Forward
Citations 
0
Petitions 
1
Assignment
First Claim
1. A driving situation detection system for observing and analyzing driving situations of a vehicle, at least including the situations of straight travel and cornering, wherein the system is capable of observing and analyzing driving situations with almost steadystate courses of the vehicle at a steering angle velocity of lower than, or equal to, 120°/s, wherein the system is designed as a state machine capable of detecting at least one of the states which are members of the group consisting of straight travel with partial braking, cornering with partial braking, and cornering with a change of load;
 and wherein the system is capable of activating a controller depending on the states which reflect a course deviation.
1 Assignment
0 Petitions
Accused Products
Abstract
A method for a Sensitive Electronic Stability Program (SESP) presents a general approach for the correction of maneuvers of turning into a bend at low speed. It integrates existing methods as well as subsequent extensions. SESP supplements the standard active yaw control (AYC) function. This allows the SESP to use variables and mechanisms of AYC, on the one hand. On the other hand, AYC continues operating unimpeded in the background and will intervene as usual when SESP cannot stabilize the vehicle appropriately. When the standard AYC intervenes, SESP control operations are forbidden, or running SESP control operations are stopped. This stop can take place either abruptly or (which is more comfortable) by way of a moderate decrease of the correcting variables.
67 Citations
View as Search Results
TRAILER SWAY MITIGATION USING MEASURED DISTANCE BETWEEN A TRAILER AND A TOW VEHICLE  
Patent #
US 20110022282A1
Filed 07/24/2009

Current Assignee
Robert Bosch LLC, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch LLC, Robert Bosch GmbH

HOLISTIC CONTROL FOR STABILIZING VEHICLETRAILER SWAYING  
Patent #
US 20110029210A1
Filed 07/30/2009

Current Assignee
Robert Bosch LLC, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch LLC, Robert Bosch GmbH

VEHICLE CONDITION ESTIMATING DEVICE  
Patent #
US 20110257876A1
Filed 01/13/2009

Current Assignee
Ibaraki Toyota Jidosha Kabushiki Kaisha

Sponsoring Entity
Ibaraki Toyota Jidosha Kabushiki Kaisha

Driving Skill Improvement Device and Driving Skill Improvement Method  
Patent #
US 20100055649A1
Filed 09/01/2009

Current Assignee
Hitachi Ltd.

Sponsoring Entity
Hitachi Ltd.

APPARATUS AND METHOD FOR DETECTING DECREASE IN TIRE AIR PRESSURE AND PROGRAM FOR DETECTING DECREASE IN TIRE AIR PRESSURE  
Patent #
US 20100126264A1
Filed 09/29/2009

Current Assignee
Sumitomo Rubber Industries Limited

Sponsoring Entity
Sumitomo Rubber Industries Limited

Method and system for coordinating a vehicle stability control system with a suspension damper control subsystem  
Patent #
US 7,706,941 B2
Filed 07/20/2006

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations Incorporated

Vehicle attitude control device based on stability factor  
Patent #
US 7,647,150 B2
Filed 07/08/2005

Current Assignee
Advics Company Limited

Sponsoring Entity
Advics Company Limited

Suspension control calibration integrity  
Patent #
US 7,684,911 B2
Filed 06/20/2005

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations Incorporated

Vehicle state determination integrity  
Patent #
US 7,747,359 B2
Filed 06/20/2005

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations LLC

METHOD AND ARRANGEMENT FOR MEASURING AND ESTIMATING A BRAKE FACTOR IN A VEHICLE BRAKE SYSTEM  
Patent #
US 20100211279A1
Filed 10/05/2007

Current Assignee
Volvo Lastvagnar AB

Sponsoring Entity
Volvo Lastvagnar AB

APPARATUS AND METHOD OF CONTROLLING VEHICLE WHEEL SPIN  
Patent #
US 20090132123A1
Filed 12/28/2007

Current Assignee
Hyundai Motor Company

Sponsoring Entity
Hyundai Motor Company

THREEWHEEL VEHICLE ELECTRONIC STABILITY SYSTEM  
Patent #
US 20090152940A1
Filed 05/05/2006

Current Assignee
Bombardier Recreational Products Incorporated

Sponsoring Entity
Bombardier Recreational Products Incorporated

Method of adjusting a motor vehicle electronic stability program  
Patent #
US 7,542,838 B2
Filed 08/04/2008

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

System and method for displaying a driving profile  
Patent #
US 7,561,054 B2
Filed 06/09/2006

Current Assignee
GREENROAD DRIVING TECHNOLOGIES LTD.

Sponsoring Entity
GREENROAD DRIVING TECHNOLOGIES LTD.

Closedloop control for trailer sway mitigation  
Patent #
US 20080036296A1
Filed 08/11/2006

Current Assignee
Robert Bosch Corporation, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch Corporation, Robert Bosch GmbH

Method and System for Coordinating a Vehicle Stability Control System with a Suspension Damper Control SubSystem  
Patent #
US 20080082234A1
Filed 07/20/2006

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations LLC

VEHICLE SYSTEMS CONTROL FOR IMPROVING STABILITY  
Patent #
US 20080183353A1
Filed 01/24/2008

Current Assignee
Honda Motor Company

Sponsoring Entity


METHOD OF ADJUSTING A MOTOR VEHICLE ELECTRONIC STABILITY PROGRAM  
Patent #
US 20080288139A1
Filed 08/04/2008

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

System and method for displaying a driving profile  
Patent #
US 20070001831A1
Filed 06/09/2006

Current Assignee
GREENROAD DRIVING TECHNOLOGIES LTD.

Sponsoring Entity
GREENROAD DRIVING TECHNOLOGIES LTD.

Driver load measuring method, device, and program for a vehicle accompanied by attitude changes, and storage medium for storing the program  
Patent #
US 7,266,444 B2
Filed 09/19/2005

Current Assignee
Honda Motor Company

Sponsoring Entity
Honda Motor Company

Driver load measuring method, device, and program for a vehicle accompanied by attitude changes, and storage medium for storing the program  
Patent #
US 20060069498A1
Filed 09/19/2005

Current Assignee
Honda Motor Company

Sponsoring Entity
Honda Motor Company

Vehicle attitude control device  
Patent #
US 20060169522A1
Filed 07/08/2005

Current Assignee
Advics Company Limited

Sponsoring Entity
Advics Company Limited

Electronic stability system on a threewheeled vehicle  
Patent #
US 20060180372A1
Filed 08/18/2004

Current Assignee
Bombardier Recreational Products Incorporated

Sponsoring Entity
Bombardier Recreational Products Incorporated

Vehicle state determination integrity  
Patent #
US 20060287777A1
Filed 06/20/2005

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations LLC

Suspension control calibration integrity  
Patent #
US 20060287790A1
Filed 06/20/2005

Current Assignee
GM Global Technology Operations LLC

Sponsoring Entity
GM Global Technology Operations LLC

Method and arrangement for measuring and estimating a brake factor in a vehicle brake system  
Patent #
US 8,108,117 B2
Filed 10/05/2007

Current Assignee
Volvo Lastvagnar AB

Sponsoring Entity
Volvo Lastvagnar AB

VEHICLE BEHAVIOR CONTROL DEVICE  
Patent #
US 20120173111A1
Filed 09/09/2010

Current Assignee
Advics Company Limited

Sponsoring Entity
Advics Company Limited

Vehicle systems control for improving stability  
Patent #
US 8,229,642 B2
Filed 01/24/2008

Current Assignee
Honda Motor Company

Sponsoring Entity
Honda Motor Company

BRAKE FLUID PRESSURE CONTROL APPARATUS FOR VEHICLE  
Patent #
US 20120212043A1
Filed 02/21/2012

Current Assignee
Autoliv Nissin Brake Systems Japan Co. Ltd.

Sponsoring Entity
Autoliv Nissin Brake Systems Japan Co. Ltd.

Apparatus and method for detecting decrease in tire air pressure and program for detecting decrease in tire air pressure  
Patent #
US 8,261,606 B2
Filed 09/29/2009

Current Assignee
Sumitomo Rubber Industries Limited

Sponsoring Entity
Sumitomo Rubber Industries Limited

METHOD FOR SETTING A LIMITING VALUE OF A VEHICLE STATE VARIABLE IN THE EVENT OF AN ACCIDENT  
Patent #
US 20120265408A1
Filed 09/07/2010

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

Holistic control for stabilizing vehicletrailer swaying  
Patent #
US 8,326,504 B2
Filed 07/30/2009

Current Assignee
Robert Bosch LLC, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch LLC, Robert Bosch GmbH

TURNING CONTROL DEVICE FOR VEHICLE  
Patent #
US 20120330526A1
Filed 03/04/2011

Current Assignee
Honda Motor Company

Sponsoring Entity
Honda Motor Company

Apparatus and method of controlling vehicle wheel spin  
Patent #
US 8,340,865 B2
Filed 12/28/2007

Current Assignee
Hyundai Motor Company

Sponsoring Entity
Hyundai Motor Company

LOOK AHEAD VEHICLE SUSPENSION SYSTEM  
Patent #
US 20130035825A1
Filed 10/08/2012

Current Assignee
Enpulz LLC

Sponsoring Entity
Enpulz LLC

INTELLIGENT VEHICLE SENSOR DEVICE  
Patent #
US 20130332030A1
Filed 03/01/2012

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Vehicle yaw rate correction  
Patent #
US 8,694,224 B2
Filed 02/28/2013

Current Assignee
Magna Electronics Incorporated

Sponsoring Entity
Magna Electronics Incorporated

Closedloop control for trailer sway mitigation  
Patent #
US 8,740,317 B2
Filed 08/11/2006

Current Assignee
Robert Bosch Corporation, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch Corporation, Robert Bosch GmbH

Look ahead vehicle suspension system  
Patent #
US 8,818,629 B2
Filed 10/08/2012

Current Assignee
Enpulz LLC

Sponsoring Entity
Enpulz LLC

SLIPCONTROLLED BRAKING SYSTEM FOR ELECTRICALLY DRIVEN MOTOR VEHICLES  
Patent #
US 20140257664A1
Filed 09/27/2012

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Trailer sway mitigation using measured distance between a trailer and a tow vehicle  
Patent #
US 8,838,353 B2
Filed 07/24/2009

Current Assignee
Robert Bosch LLC, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch LLC, Robert Bosch GmbH

Device for estimating state quantity of skid motion of vehicle  
Patent #
US 8,855,885 B2
Filed 03/26/2010

Current Assignee
Honda Motor Company

Sponsoring Entity
Honda Motor Company

TRACTION SYSTEM AND A METHOD FOR CONTROLLING SAID TRACTION SYSTEM  
Patent #
US 20150081190A1
Filed 05/08/2012

Current Assignee
Volvo Lastvagnar AB

Sponsoring Entity
Volvo Lastvagnar AB

Vehicle behavior control device  
Patent #
US 9,031,759 B2
Filed 09/09/2010

Current Assignee
Advics Company Limited

Sponsoring Entity
Advics Company Limited

Towing Vehicle Controller Providing Brake Control to a Towed Vehicle and Method  
Patent #
US 20150210257A1
Filed 04/09/2015

Current Assignee
Bendix Commercial Vehicle Systems LLC

Sponsoring Entity
Bendix Commercial Vehicle Systems LLC

Vehicle condition estimating device  
Patent #
US 9,096,232 B2
Filed 01/13/2009

Current Assignee
Ibaraki Toyota Jidosha Kabushiki Kaisha

Sponsoring Entity
Ibaraki Toyota Jidosha Kabushiki Kaisha

Brake fluid pressure control apparatus for vehicle  
Patent #
US 9,260,096 B2
Filed 02/21/2012

Current Assignee
Autoliv Nissin Brake Systems Japan Co. Ltd.

Sponsoring Entity
Autoliv Nissin Brake Systems Japan Co. Ltd.

Method for controlling braking torque in a twowheeled vehicle when traveling in an inclined position  
Patent #
US 9,302,660 B2
Filed 07/18/2011

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

User interface for driver performance application  
Patent #
US 9,342,983 B1
Filed 10/23/2012

Current Assignee
GREENROAD DRIVING TECHNOLOGIES LTD.

Sponsoring Entity
GREENROAD DRIVING TECHNOLOGIES LTD.

VEHICLE CONTROL APPARATUS, VEHICLE CONTROL SYSTEM, AND VEHICLE CONTROL METHOD  
Patent #
US 20160251014A1
Filed 01/31/2016

Current Assignee
Fujitsu Ten Limited

Sponsoring Entity
Fujitsu Ten Limited

Towing vehicle controller providing brake control to a towed vehicle and method  
Patent #
US 9,592,800 B2
Filed 04/09/2015

Current Assignee
Bendix Commercial Vehicle Systems LLC

Sponsoring Entity
Bendix Commercial Vehicle Systems LLC

Traction system and a method for controlling said traction system  
Patent #
US 9,597,950 B2
Filed 05/08/2012

Current Assignee
Volvo Lastvagnar AB

Sponsoring Entity
Volvo Lastvagnar AB

Method for setting a limiting value of a vehicle state variable in the event of an accident  
Patent #
US 9,604,610 B2
Filed 09/07/2010

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

Slipcontrolled braking system for electrically driven motor vehicles  
Patent #
US 9,744,862 B2
Filed 09/27/2012

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Nonlinear compensation controller for active steering system in a vehicle  
Patent #
US 9,878,738 B2
Filed 03/15/2013

Current Assignee
Robert Bosch LLC, Robert Bosch GmbH

Sponsoring Entity
Robert Bosch LLC, Robert Bosch GmbH

Steering and braking control system  
Patent #
US 9,925,988 B1
Filed 02/24/2017

Current Assignee
Ford Global Technologies LLC

Sponsoring Entity
Ford Global Technologies LLC

Detecting activation of an operator control element in a motor vehicle  
Patent #
US 10,095,566 B2
Filed 04/08/2016

Current Assignee
Audi AG

Sponsoring Entity
Audi AG

Vehicle brake hydraulic pressure control apparatus and road surface friction coefficient estimating device  
Patent #
US 10,239,507 B2
Filed 11/28/2016

Current Assignee
Autoliv Nissin Brake Systems Japan Co. Ltd.

Sponsoring Entity
Autoliv Nissin Brake Systems Japan Co. Ltd.

Brake fluid pressure control apparatus for vehicle  
Patent #
US 10,328,916 B2
Filed 01/06/2016

Current Assignee
Autoliv Nissin Brake Systems Japan Co. Ltd.

Sponsoring Entity
Autoliv Nissin Brake Systems Japan Co. Ltd.

Method for modifying a driving stability control of a vehicle  
Patent #
US 20050004738A1
Filed 06/28/2002

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Method for regulating the driving stability of a vehicle  
Patent #
US 6,957,873 B2
Filed 10/30/2001

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Method of influencing the torque delivered by a vehicle drive of a vehicle  
Patent #
US 6,611,747 B1
Filed 11/05/2001

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

Method of improving ABS control behaviour on cornering  
Patent #
US 6,289,272 B1
Filed 09/17/1999

Current Assignee
Continental Teves AG Company HG

Sponsoring Entity
Continental Teves AG Company HG

Method for improving the controllabilty of motor vehicles  
Patent #
US 5,668,724 A
Filed 09/11/1995

Current Assignee
Robert Bosch GmbH

Sponsoring Entity
Robert Bosch GmbH

System for integrated driving stability control  
Patent #
US 5,694,321 A
Filed 06/07/1995

Current Assignee
ITT Automotive Europe GmbH

Sponsoring Entity
ITT Automotive Europe GmbH

Vehicle brake system including cornering characteristic control  
Patent #
US 5,172,961 A
Filed 07/05/1991

Current Assignee
Nissan Motor Co. Ltd.

Sponsoring Entity
Nissan Motor Co. Ltd.

Steering and brake controlling system  
Patent #
US 4,998,593 A
Filed 03/31/1989

Current Assignee
Aisin Seiki Co. Ltd.

Sponsoring Entity
Aisin Seiki Co. Ltd.

8 Claims
 1. A driving situation detection system for observing and analyzing driving situations of a vehicle, at least including the situations of straight travel and cornering, wherein the system is capable of observing and analyzing driving situations with almost steadystate courses of the vehicle at a steering angle velocity of lower than, or equal to, 120°/s, wherein the system is designed as a state machine capable of detecting at least one of the states which are members of the group consisting of straight travel with partial braking, cornering with partial braking, and cornering with a change of load;
 and wherein the system is capable of activating a controller depending on the states which reflect a course deviation.
 View Dependent Claims (2)
 3. A device for controlling the yaw performance of a twoaxle, fourwheel vehicle which is equipped with a hydraulic brake system that includes a brake for each wheel, and with at least one sensor system or calculation models for establishing values which are members of the group consisting of wheel speeds, lateral acceleration, yaw velocity, steering angle, the steering angle velocity;
 comprising at least two electronic yaw torque controllers, the first controller realizing an additional yaw torque MG_{SESP }to improve the track performance of the vehicle depending on a comparison between desired and actual yaw rates below an ESP intervention threshold and the second controller (10) putting this torque into practice to stabilize the driving state of the vehicle above the ESP intervention threshold by a corresponding brake application.
 4. A comfort control system for a vehicle with an ESP controller for controlling an additional yaw torque, with a situation detection system in which different driving situations are distinguished and taken into consideration in the control, the control being active or inactive depending on control thresholds, wherein an additional driving situation detection system is provided in which almost steadystate shapes of the course of the vehicle with a steering angle velocity of lower than, or equal to, 120°/s are observed and analyzed, and in that depending on expected or actually determined comfort control deviations in the steadystate range indicative of a request by the driver and of the vehicle state, a modification of a reference value representing the desired course of the vehicle is activated when the control by the ESP controller is inactive.
 6. A method of improving the control performance of a control system (ESP) wherein the rotational behavior of the vehicle is measured and evaluated as an input variable of the control system, wherein logical linking of the yaw rate with a desired course of the vehicle as intended by the driver is used to produce control variables for changing the vehicle course, and the control variables are set inactive or active depending on control thresholds in order to change the vehicle course, the method comprising the following steps:
 changing the control thresholds depending on driving situations, observing and analyzing deviations from the desired vehicle course with an almost steadystate shape of the course of the vehicle at a steering angle velocity of lower than, or equal to, 120°/s, and modifying a reference value representing the desired course of the vehicle in the event of expected or actually determined deviations in the steadystate range, when the control thresholds are below the control thresholds set higher steering angle velocities.
 View Dependent Claims (7)
 8. Method as claimed in claim 7, wherein the modelbased vehicle reference speed is determined from at least one member of the group of quantities consisting of yaw velocity, steering wheel angle, lateral acceleration, steering wheel angle velocity, yaw acceleration, a quantity representative of any of the foregoing quantities, respectively.
1 Specification
This patent application is a ContinuationInPart of the International Patent Application PCT/EP2004/052689, international filing date Oct. 28, 2004.
BACKGROUND OF THE INVENTION1. General Structure of Driving Stability Control (DSC)
The term driving stability control (DSC) covers four principles of influencing the driving behavior of a vehicle by means of predeterminable pressures in individual wheel brakes and by interfering with the engine management of the driving engine. These include brake slip control (ABS) which is to prevent locking of individual wheels during a brake operation, traction slip control (TCS), which prevents the spinning of the driven wheels; the electronic brake force distribution system (EBD), which controls the ratio of the brake forces between the front and the rear axle of the vehicle; and a yaw torque control system (YTC), which ensures stable driving conditions during travel in a curve.
Consequently, a vehicle is defined in this connection as a motor vehicle with four wheels, which is equipped with a hydraulic brake system. In a hydraulic brake system, brake pressure can be built up by the driver by means of a pedalactuated master cylinder. Each wheel has a brake, with which one inlet valve and one outlet valve each is associated. The wheel brakes communicate with the master cylinder by way of the inlet valves, while the outlet valves lead to a pressureless tank or to a lowpressure accumulator. Finally, there also is an auxiliary pressure source, which is able to build up a pressure in the wheel brakes regardless of the position of the brake pedal. The inlet and outlet valves can be electromagnetically actuated for pressure regulation in the wheel brakes.
To detect states in the dynamics of the vehicle movement, there are four speed sensors, one per wheel, one yaw rate meter, one lateral accelerometer, and at least one pressure sensor for the brake pressure generated by the brake pedal. The pressure sensor may be replaced with a pedal travel or pedal force meter if the auxiliary pressure source is arranged such that a brake pressure built up by the driver is not distinguishable from that of the auxiliary pressure source. In addition, it is possible to poll information about the condition of the transmission, e.g. about gearshift control, etc.
A fallback solution is advantageously put into practice in light of such a large number of sensors. This means that, in the case of failure of part of the sensor system, only the component of the control system that depends on that part is switched off. If, for example, the yaw rate meter fails, no yaw torque control can be performed, but the ABS, TCS and EBD systems continue to function. The driving stability control can consequently be limited to these other three functions.
In a driving stability control, the driving behavior of a vehicle is influenced such that the driver will be better able to master the vehicle in critical situations, or critical situations will be avoided to begin with. A critical situation is defined herein as an unstable driving condition in which, in the extreme case, the vehicle does not follow the driver's instructions. The function of the driving stability control is consequently to impart to the vehicle the behavior desired by the driver in such situations within the physical limits.
While the longitudinal slip of the tires on the road surface is mainly of significance for the brake slip control system, the traction slip control system and the electronic brake force distribution system, the yaw torque control system (YTC) also involves additional variables, e.g., the yaw rate {dot over (Ψ)}.
Various vehicle reference models may be used for yaw torque control. The calculation is simplest on the basis of a singletrack model, i.e., the front wheels and the rear wheels are integrated in this model into one wheel each, which is located on the longitudinal axis of the vehicle. The calculations become considerably more complicated if they are based on a twotrack model. However, since lateral displacements of the center of gravity (rolling movements) can also be taken into account in the twotrack model, the results are more accurate.
The system equations<maths id="MATHUS00001" num="1"><math overflow="scroll"><mtable><mtr><mtd><mrow><mover><mi>β</mi><mo>.</mo></mover><mo>=</mo><mrow><mrow><msub><mi>c</mi><mn>11</mn></msub><mo></mo><mfrac><mi>β</mi><mi>v</mi></mfrac></mrow><mo></mo><mover><mi>Ψ</mi><mo>.</mo></mover><mo>+</mo><mrow><msub><mi>c</mi><mn>12</mn></msub><mo></mo><mfrac><mover><mi>Ψ</mi><mo>.</mo></mover><msup><mi>v</mi><mn>2</mn></msup></mfrac></mrow><mo>+</mo><mrow><msub><mi>c</mi><mn>13</mn></msub><mo></mo><mfrac><mi>δ</mi><mi>v</mi></mfrac></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1.1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>Ψ</mi><mi>¨</mi></mover><mo>=</mo><mrow><mrow><msub><mi>c</mi><mn>21</mn></msub><mo></mo><mi>β</mi></mrow><mo>+</mo><mrow><msub><mi>c</mi><mn>22</mn></msub><mo></mo><mfrac><mover><mi>Ψ</mi><mo>.</mo></mover><mi>v</mi></mfrac></mrow><mo>+</mo><mrow><msub><mi>c</mi><mn>23</mn></msub><mo></mo><mi>δ</mi></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1.2</mn></mrow></mtd></mtr></mtable></math></maths>can be written in the phase space diagram for a singletrack model.
The sideslip angle β and the yaw rate {dot over (Ψ)} represent the phase variables of the system. The input variable acting on the vehicle is the steering angle δ, as a result of which the vehicle receives the yaw rate {dot over (Ψ)} as an output variable. The model coefficients c<sub>ii </sub>are formed as follows:<maths id="MATHUS00002" num="2"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>c</mi><mn>11</mn></msub><mo>=</mo><mrow><mo></mo><mfrac><mrow><msub><mi>c</mi><mi>h</mi></msub><mo>+</mo><msub><mi>c</mi><mi>v</mi></msub></mrow><mi>m</mi></mfrac></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><msub><mi>c</mi><mn>12</mn></msub><mo>=</mo><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow><mi>m</mi></mfrac></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><msub><mi>c</mi><mn>13</mn></msub><mo>=</mo><mfrac><msub><mi>c</mi><mi>v</mi></msub><mi>m</mi></mfrac></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><msub><mi>c</mi><mn>21</mn></msub><mo>=</mo><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow><mi>Θ</mi></mfrac></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><msub><mi>c</mi><mn>22</mn></msub><mo>=</mo><mrow><mo></mo><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msubsup><mi>l</mi><mi>h</mi><mn>2</mn></msubsup></mrow><mo>+</mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msubsup><mi>l</mi><mi>v</mi><mn>2</mn></msubsup></mrow></mrow><mi>Θ</mi></mfrac></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><msub><mi>c</mi><mn>23</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mi>Θ</mi></mfrac></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1.3</mn></mrow></mtd></mtr></mtable></math></maths>
c<sub>h </sub>and c<sub>v </sub>are the resulting rigidities from the elasticity of the tire, wheel suspension and steering on the rear axle and the front axle, respectively. l<sub>h </sub>and l<sub>v </sub>are the distances of the rear axle and the front axle, respectively, from the center of gravity of the vehicle. Θ is the moment of inertia of the vehicle, i.e., the moment of inertia of the vehicle around its vertical axis.
Longitudinal forces and displacements of the center of gravity are taken into account in this model. This approximation is also valid only for low angular velocities. Consequently, the accuracy of this model decreases with decreasing curve radii and increasing velocities. However, the amount of calculations is manageable. Further explanations of this singletrack model can be found in the book Fahrwerktechnik: Fahrverhalten [Chassis Engineering: Driving Behavior] by Adam Zomotor, Vogel Buchverlag, Wurzburg, 1987.
A twotrack model, whose accuracy is superior to that of a singletrack model, is proposed for a vehicle in DE40 30 704 A1. The yaw rate {dot over (Ψ)} and the sideslip angle β form the phase variables in this case as well. However, when a twotrack model is used, it is necessary to consider the fact that an enormous calculation capacity is needed to make it possible to perform a control intervention in a relatively short time.
The methods and the control systems are used to create an additional torque by targeted interventions at the individual brakes of a vehicle, which torque leads by way of the actually measured yaw variation per time unit (actual yaw rate) of a vehicle to the yaw variation per time unit (desired yaw rate) which is influenced by the driver. Hence, a method and a control system will particularly intervene in a supporting manner into the steering performance of the vehicle when due to certain conditions (e.g. excessive speed, slippery roadway) the curved track actually covered by the vehicle does not correspond to the curved track desired by the driver without additional torque. In principle, methods and control systems of this type to improve driving stability have already been described comprehensively and, therefore, shall not be explained again in detail. In methods and control systems of this type, input parameters which result from the curved track desired by the driver (e.g. steering wheel angle, driving speed) are always sent to a vehicle model circuit which, on the basis of a known singletrack model or another driving model, determines a nominal yaw rate (ψ<sub>Desired</sub>) from these input parameters and from parameters being characteristic of the driving performance of the vehicle, yet also from quantities (coefficient of friction of the roadway) predetermined by the characteristics of the ambience, which yaw rate is compared to the measured actual yaw rate (ψ<sub>actual</sub>). The difference of the yaw angles (Δψ<sub>Diff</sub>) is converted by means of a yaw torque controller into an additional yaw torque M<sub>G </sub>which forms the input parameter of a distribution logic unit.
The distribution logic unit itself defines the brake pressure to be applied to the individual brakes, optionally in dependence on the braking request of the driver demanding a defined brake pressure at the wheel brakes. This brake pressure shall produce an additional torque at the vehicle, in addition to the optionally desired brake effect, which supports the driving behavior of the vehicle in the direction of complying with the steering request of the driver.
If due to external conditions during driving or conditions being due to the performance of the driver, variations of the vehicledynamic driving performance (e.g. changes in the coefficient of friction) will occur, e.g. a change in the engine torque, e.g. due to release of the accelerator or pushdown of the accelerator, or due to braking, the driving behavior of the vehicle will change because, among others, there will be a change of the axle load and, thus, of forces which is induced by the interaction of several influences such as tire influences, kinematical influences and elastokinematical influences.
For example, when driving through a curve and release of the accelerator occurs, the driving forces F<sub>a </sub>are active at the drive wheels prior to the release of the accelerator. Due to the lateral deformation of the tire contact area, the longitudinal driving force F<sub>A</sub>=2×F<sub>a </sub>will act in dependence on the lateral forces slightly outside the wheel center plane. An understeering yaw torque ({dot over (Ψ)}<sub>understeer</sub>) is produced due to the longitudinal driving force F<sub>A </sub>that acts asymmetrically in relation to the vehicle longitudinal axis.
After the release of the accelerator, the engine (and other resistances) slows down the vehicle, the (longitudinal) driving forces become negative. In addition, the deceleration produces an inertia force mx in the point of gravity SP, with the result that the axle load increases at the front wheels and decreases by the same amount at the rear wheels. Thus, the distribution of the transmittable lateral forces changes. The lateral force change (lateral force on the front axle rises slightly and decreases greatly at the rear axle) produces an oversteering yaw torque ({dot over (Ψ)}<sub>oversteer</sub>), the slip angles at the rear axle will increase and the vehicle turns into the curved track. When the engine torque is changed from driving force into brake force, the reversal of theses torques will induce a change of the driving behavior of the vehicle from an oversteering to an understeering driving behavior.
There are methods fulfilling only partial aspects of the mentioned requirement.
Thus, the function ‘ABSplus’ is known for partial brake operations in the curve. This function achieves stabilization of the vehicle by pressure reduction on the inside wheels in a turn. However, ABSplus detects the vehicle performance exclusively from the measured wheel speeds.
When the driver slows down in a curve to such an extent that ABS control is triggered, already the ABS function itself is often capable of counteracting the tendency of turning into a bend. The reason for this is that a greater vertical force and, thus, a higher potential of longitudinal force prevails on the exterior curve side than on the interior curve side. ABS safeguards the optimum utilization of the longitudinal force potential. The resulting unbalance of forces will then bring about a stabilizing yaw torque.
It is disadvantageous that these methods do not have a controller of their own but share the yaw torque controllers with the standard ESP. Their effect evolves because they take influence on the parameters of the YTC controllers (e.g. decrease of the control thresholds).
The previous methods suffer from the following disadvantages:
 1. Each method acts only in one or in a few defined driving situations and is limited to one defined intervention strategy only.
 2. Each method includes nonoptimal partial solutions; e.g. the comfort is optimal in one method because the hydraulic pump provides its full output for pressure buildup. In the other method, the reference yaw rate is not used continuously.
 3. To cover a rising number of driving situations by a simultaneous activation of an increasing rate of single methods will quickly reach limits because the ranges of influence (driving situations) of the methods can overlap each other what is not desired, or they may leave gaps, and the intervention strategies contradict each other in part.
In view of the above, an object of the invention is to provide a method and a control system for improving the driving behavior of a vehicle which influences the driving behavior of a vehicle in such a manner that it is adapted early and comfortably to the desired driving behavior depending on the driving situation. Another objective is to enhance the sensitivity of response of the control.
SUMMARY OF THE INVENTIONThis object is achieved in that driving situations with almost steadystate courses of the vehicle at a steering angle velocity of lower than, or equal to, 120°/s are observed and analyzed in the driving situation detection system, and the driving situation detection system is designed as a state machine which detects at least one of the states such as straight travel with partial braking, cornering with partial braking, and/or cornering with a change of load, and activates a controller depending on the states which reflect a course deviation.
One embodiment of the driving situation detection system is connected to a reference signal model in which a YTC reference variable with an offset is corrected, with the driving situation detection system determining the time at which the offset correction is started and terminated and the offset value is superposed on the reference variable at the starting point.
The device for controlling the yaw performance of a twoaxle, fourwheel vehicle is equipped with a hydraulic brake system that includes a brake for each wheel, and with at least one sensor system or models for gathering or calculating, respectively, the wheel speeds, the lateral acceleration, the yaw velocity, the steering angle, and/or the steering angle speed. It is also equipped with at least two electronic yaw torque controllers, with the first controller realizing an additional yaw torque MG<sub>SESP </sub>to improve the track performance of the vehicle depending on a comparison between nominal and actual yaw rates below the ESP intervention threshold and the second controller putting this torque into practice to stabilize the driving state of the vehicle above the ESP intervention threshold by a corresponding brake application.
In a first situation detection system, different driving situations such as straight travel, cornering, and the like, are distinguished and taken into consideration in the control. This control is active or inactive depending on control thresholds. Another driving situation detection system is provided in which almost steadystate shapes of the course of the vehicle with a steering angle speed of lower than, or equal to, 120°/s are observed and analyzed, and depending on expected or actually determined deviations in the steadystate range indicative of the driver's request and the vehicle state, and activates a modification of a reference value representing the desired course of the vehicle when the control by the ESP controller is inactive.
An embodiment of the driving stability control further determines the desired course of the vehicle in the ESP control and evaluates deviations from the desired course, and it determines during a desired course of driving if there is a tendency of a subsequent deviation from course due to a partial braking situation or load change situation.
The rotational performance of the vehicle is measured and evaluated as an input variable of the control system, wherein logical linking of the yaw rate with the course of the vehicle desired by the driver is used to produce control variables for changing the vehicle course, and the control variables are set inactive or active to change the vehicle course depending on control thresholds. The control thresholds are changed depending on driving situations, and with an almost steadystate shape of the course of the vehicle at a steering angle velocity of lower than, or equal to, 120°/s deviations from the desired vehicle course are observed and analyzed. With expected or actually determined deviations in the steadystate range, a reference value representing the desired course of the vehicle is modified when the control thresholds are below the ESP control thresholds.
An embodiment of the method provides that the modelbased vehicle reference speed is calculated from the sensor signals of the driving stability control DSC.
Another detail of the method includes that the modelbased vehicle reference speed is determined from the yaw velocity, the steering wheel angle, and/or the lateral acceleration and/or their derivatives or substitute signals, respectively.
The method SESP (Sensitive Electronic Stability Program) is applicable to slow maneuvers of turning into a bend. Such a driving situation exists when all subsequently mentioned conditions are satisfied:
 1. The driver wishes steadystate driving, that, means his steering dynamics is at a relatively low rate. This may concern both straight travel and cornering.
 2. The vehicle follows the driver's request at least in the beginning.
 3. A deviation from the course is slowly developing during steadystate driving in the sense that the vehicle leaves the desired straight travel or turns to the inside of a bend in the desired cornering maneuver (‘oversteering’).
 4. There is a measurable cause advanced in time for the deviation from the course.
Due to the restriction to oversteering situations, in particular such causes are at issue in this respect which result in an asymmetrical distribution of the wheel longitudinal forces (e.g. inhomogeneous roadway, asymmetrical brake wear) during straight travel, while they induce a shift of the wheel lateral force potential to the front axle (e.g. braking maneuver or change of load) during cornering.
Meanwhile, there is a demand in the automotive industry that the ESP controller (ESP (YTC)=Electronic Stability Program) assists the driver in governing such slow maneuvers of turning into a bend. As this occurs, it is desired that the correcting interventions shall be unnoticed by the driver to the greatest extent possible.
The standard AYC(YTC)module (AYC=Active Yaw Control) of the ESP controller satisfies this demand only insufficiently because it intervenes only in the event of noticeable course deviations for reasons of robustness, and these interventions are often felt distinctly by the driver, e.g. by feedback of the pressure increases through the brake pedal and through the steering wheel.
Therefore, the AYCmodule must be extended to comply with the desired purpose.
The new method SESP presents a general approach for the correction of maneuvers of turning into a bend in the sense of the requirement described hereinabove. It offers a structure which is appropriate to integrate existing methods as well as subsequent extensions. FIG. 35 shows the structure of SESP.
SESP is designed to be a supplement to the standard AYC function. This allows the SESP to use variables and mechanisms of AYC, on the one hand. On the other hand, AYC continues operating unimpeded in the background and will intervene as usual when SESP cannot stabilize the vehicle appropriately. When the standard AYC intervenes, SESP control operations are forbidden, or running SESP control operations are stopped. This stop can take place either abruptly or (what is more comfortable) by way of a moderate decrease of the correcting variables, FIG. 35.
The single elements of SESP will be described briefly in the following.
The driving situation detection system uses information about the driver's request (e.g. steering wheel angle, engine torque demand, brake pressure) and the current vehicle state (e.g. lateral acceleration, estimated curve radius, speed) in order to decide whether there is a potential turning situation. The detection is configured as a state machine. FIG. 30 shows as an example the possible states and permitted state transitions in the current SESP realization. The use of a state machine safeguards the unambiguousness of the detected driving situation. Making a distinction between several states is necessary in order to adapt the SESP controller optimally to the current driving situation. It is constantly checked in a parallel fashion (e.g. by way of the AYC state, ABS state, transversal inclination of the roadway, vehicle speed) whether SESP interventions are permitted. The state ‘SESP inactive’ is instantaneously assigned in the event of a negative result. Some causes of turning (e.g. change of load) naturally have an effect on the vehicle performance which is only limited in time. In these cases, the respective state in the SESP driving situation detection system is left again after a defined time. Another argument for the temporal limitation of such states is that the driver is frequently able to master slow turning maneuvers after a certain ‘time of adaptation’ by himself.
As has been mentioned before, the detected state of the driving situation detection system is used to activate the SESP controllers and to adapt their parameters optimally to the current driving situation (e.g. the control thresholds). In order to be able to trigger the sensitive SESP controllers as early as possible, the SESP driving situation detection system not only assesses the vehicle reaction but already the causes before the reaction. In addition to influencing the parameters of the controller, another objective of the driving situation detection system is to indicate to the reference signal formation unit when the offset for the offset correction of the reference variables shall be stored. This is likewise done already whenever a possible cause for a turning situation is detected.
SESP uses the principle of an offset correction for the reference signal formation. The driver's request is represented by the SESP reference yaw rate. Compared to the AYC reference yaw rate, the SESP reference yaw rate has an offset which is just so high that the SESP yaw rate deviation is zero at the commencement of the load change or the braking intervention.
Due to this offset, possible uncompensated offsets of the yaw rate sensor are compensated, on the one hand. This is necessary to prevent control errors of the sensitive SESP yaw rate controller. On the other hand, only the deviations from the current vehicle performance are regulated by the offset compensation.
FIG. 31 shows the formation of the SESP reference yaw rate for the case of a positive yaw rate sensor offset.
The number and type of the SESP controllers can be made dependent on the case of application. The regulating variable of the controller is the yaw rate, and an additional yaw torque is the correcting variable. However, it is also easily possible to consider e.g. the regulating variable, namely, sideslip angle velocity, or the correcting variable, namely, additional engine torque. Further possible correcting variables are (if corresponding actuators are provided): additional steering angle, change of the spring/damper characteristics, change of the engine torque distribution between front and rear axles. The purpose of the yaw torque controller is to calculate the additional yaw torque, which is necessary for the vehicle stabilization, from the deviation of the actual yaw rate from the SESP reference yaw rate.
The controller configured as a pure pcontroller includes a threshold for the SESP yaw rate deviation as an activation criterion apart from the existence of a permitted driving situation. Depending on the driving situation, this threshold is maximally as high as in the case of the standard AYC yaw rate controller, see FIG. 33. Another threshold is used for the situation ‘load change in the curve’.
In the SESP module ‘arbitration’, all equal correcting variables (e.g. additional yaw torques) of the SESP controllers are combined according to defined priority rules. The selection of these priority rules can again be made dependent on the case of application. Maximum, addition and suppressing contradicting demands are wellknown examples for such rules. An additional objective of the arbitration is to coordinate the SESP correcting variables with the demands of other controllers (in particular of standardAYC) and to ensure ‘smooth’ transitions if requested.
The correcting variables are put into practice in the last step. In doing so, it has to be taken into account that the driver must not feel the SESP interventions, if possible. This also involves that he/she does not receive any information about SESP interventions from the combination instrument. Standard mechanisms of the ESP controller are coused in this step. An interface conveys the strategy desired by SESP to these mechanisms.
The application of the additional yaw torque in the current SESP realization will be explained as an example in the following, cf. FIG. 32. It is initially checked whether the driver brakes using a minimum pressure. In the affirmative, the additional yaw torque is realized by way of pressure reduction at the inside rear wheel in a turn. This type of intervention is comfortable because it manages, on the one hand, with highly reduced valve and pump activities, as compared to an active pressure buildup, and avoids a reaction through the steering system by intervention at the rear axle, on the other hand. In addition, the increase of the lateral force potential of the inside rear wheel in a turn supports the transverse stability of the vehicle. Pressure is additionally reduced at the inside front wheel in a turn, if required. When the brake pressure of the driver is below the initially mentioned minimum pressure, the additional yaw torque is realized by way of pressure increase at the outside rear wheel in a turn and, if required, at the outside front wheel in a turn in addition. To maintain the lateral force potential of the wheel concerned, the pressure is limited in an upward direction to the level of the maximum longitudinal force between wheel and road (by employing an already existing ‘slip monitor’). To prevent this pressure level from being exceeded in any case, the ABS controller is moved to a sensitive mode with a view to achieving double safety. Simultaneously, the EBD controller is disabled to prevent undesirable pressure limitations. To ensure the comfort of the SESP pressure increase mode, the wheel pressure is adjusted by means of the already existing lownoise method ‘EVR (Electrical Reversing Valve) control’.
To minimize the risk of a destabilization of the rear axle due to active pressure buildup when a low coefficient of friction of the roadway prevails, it is suggested to distribute the more pressure from the curveoutside rear wheel to the curveoutside front wheel, the lower the estimated coefficient of friction of the roadway is. Further, it can be favorable to choose a gradual transition instead of the ‘binary’ decision between pressure reduction mode and pressure increase mode in such a manner that initially the potential for pressure reduction is utilized. The part of the additional yaw torque which thereafter is possibly not yet taken into account will then be realized by a pressure increase in the manner described above within the limits of the available longitudinal force potential. This part would be left out of account in the case of a purely ‘binary’ decision, FIG. 32.
If the information is available in the ESP controller that the vehicle has a trailer, it is suggested to use this information for SESP to the effect that the vehicle longitudinal deceleration, which is produced by SESP interventions (such as pressure buildup), is limited to a level that does not destabilize the trailer. In the simplest case, SESP pressure buildup can e.g. be prevented when a trailer is detected.
In order to comply with correcting variables, it is to be noted that if several actuators are provided for one correcting variable, it is necessary to proportionally share the correcting variable among the actuators. E.g. comfort or efficacy is applicable as distribution criteria. The criterion reaction time can be considered e.g. according to the following principle. To begin with, the correcting variable is submitted completely to the quickest responding actuator. The component which cannot be realized by this actuator is then submitted to the second quickest actuator, etc. A minimum total reaction time is reached with this principle.
It is to be understood in addition that in order to comply with the additional yaw torque, it is initially checked whether the driver is braking with a minimum pressure of e.g. 20 bar. In the affirmative, the additional yaw torque is put into practice by a comfortable pressure reduction at the inside rear wheel in a turn (and, optionally, at the inside front wheel in a turn).
When the brake pressure of the driver is below the initially mentioned minimum pressure, the total additional yaw torque is realized by way of pressure buildup at the outside rear wheel in a turn. In order to maintain the lateral force potential of the wheel concerned, the pressure is limited in an upward direction to the level of the maximum longitudinal force between wheel and road (by employing an already existing ‘wheel slip monitor’). To prevent this pressure level from being exceeded in any case, the ABS controller is set to a sensitive mode with a view to achieving double safety. Simultaneously, activation of the EBD controller is prevented to rule out undesirable pressure limitations.
To ensure the comfort of the SESP also in the pressure increase mode, the wheel pressure in this mode is adjusted by means of the already existing lownoise method ‘EVR (Electrical Reversing Valve) control’.
Referring to the example of a left turn, FIG. 33 shows on which wheels SESP performs pressure reduction (−) and pressure increase (+).
Advantages of the New Method SESP
The new method SESP has the following advantages over existing methods:
 1. The structure of SESP allows combining existing methods which, when looked at individually, cover only a small range of effect, in order to achieve an overall method with a larger range of action.
 2. The structure of SESP is open to extensions. It is possible to add new situations, regulating variables and correcting variables, as desired. Thus, it has the potential of dealing with turning maneuvers within the limits of the available sensor system and actuator system.
 3. The clearcut isolation of the partial modules of controller and correcting variable realization from the driving situation detection system allows using the same working mechanisms for various driving situations. This fact ensures a constant control quality and facilitates the application and further development of SESP.
 4. Compared to the modelbased approach of the ‘sideslip angle control’, SESP within the limits of its range of action is not dependent on additional sensor equipment (such as individual wheel force sensors or wheel braking pressure sensors). The conventional ESP sensor equipment is sufficient.
 5. The consistent evaluation of the temporally preceding causes of turning maneuvers and the only selective controller activation allows the SESP to employ relatively sensitive control thresholds. This enables SESP in contrast to standard AYC to detect the slow turning maneuvers, on the one hand, and to correct them early using relatively little correcting energy, on the other hand.
 6. Using own SESP controllers allows the standard AYC controller to continue running in the background and constantly checking whether the stabilizing effect of SESP is sufficient. If this is not the case, AYC will intervene as usual.
 7. The continuous use of the described offsetcorrected reference variables reduces the disturbing influence of signal errors.
 8. The proposed strategy for realizing the additional yaw torque offers great comfort both in the pressure reduction mode and the pressure increase mode. It is additionally designed consistently to maintain vehicle stability.
The significant improvement is achieved by the clearcut separation between driving situation detection, control, and the realization of correcting variables. SESP integrates in this manner existing methods, optimizes them and is open to subsequent extensions. This fact enables SESP to significantly extend the range of action of the ESP controller without using additional sensor equipment or actuator equipment.
BRIEF DESCRIPTION OF THE DRAWINGSIn the drawings,
FIG. 1 is a block diagram relating to the overall structure of a system for driving stability control;
FIG. 2 is a block diagram relating to the structure of the yaw torque controllers;
FIG. 3 is a flow chart relating to the detection of a driving situation such as a cornering maneuver;
FIGS. 4 and 5 show in each case a flow chart relating to the determination of the coefficient of friction of the roadway, and FIG. 5 should be inserted into FIG. 4;
FIGS. 6 and 8 show block diagrams relating to a combined method to determine the current values of the sideslip angle velocity and of the sideslip angle in different types of representation;
FIG. 7 is a block diagram for the direct determination of the sideslip angle velocity from kinematical considerations as a part of the combined method of FIG. 6;
FIG. 9 is a control circuit for driving stability control with a drivingspeedresponsive change of the computing model for the vehicle;
FIGS. 10 and 11 show diagrams exhibiting the dependency of the slip angle difference of a vehicle on the sideslip angle and the speed vector of the individual wheels;
FIGS. 12 to 15 show a block diagram of a control circuit for controlling the driving stability wherein the variables being compared with each other in the comparator represent derivatives of the yaw velocity;
FIG. 16 is a control circuit for determining the driving stability wherein the pressure gradient and/or the valve on/off time of the vehicle brakes is employed as a regulating variable;
FIG. 17 is a block diagram for describing the controller in order to calculate the additional yaw torque;
FIG. 18 is a block diagram for describing a lowpass filter;
FIG. 19 is a flow chart for calculating a corrected nominal yaw velocity;
FIG. 20 is a block diagram for calculating a corrected additional yaw torque;
FIG. 21 is a schematic view of a motor vehicle;
FIG. 22 is a block diagram for describing the distribution logic unit;
FIG. 23 is a schematic view of a motor vehicle and the acting forces when the steering wheel is turned;
FIG. 24 is a diagram for describing the coefficients of the lateral force and longitudinal force depending on the wheel slip;
FIG. 25A, B is a schematic illustration of motor vehicles for describing the understeering and oversteering performance;
FIG. 26 is a flow chart with a decision logic unit within the distribution logic unit;
FIG. 27 is a block diagram for calculating onoff times for inlet and outlet valves;
FIG. 28 is a diagram for describing time intervals within a calculation loop;
FIG. 29 is a basic block diagram for determining the wheel brake pressure;
FIG. 30 is a basic block diagram of the driving situation detection system;
FIG. 31a, b shows the formation of the SESP reference yaw rate for the case of a positive yaw rate sensor offset;
FIG. 32a, b illustrates the application of additional yaw torque in the current SESP realization;
FIG. 33 shows on which wheels SESP performs pressure reduction (−) and pressure increase (+), referring to the example of a left turn;
FIG. 34 illustrates how the risk of SESP control errors is reduced;
FIG. 35 shows the structure of SESP.
DETAILED DESCRIPTION OF THE DRAWINGSA general description of the course of a driving stability control (DSC) is shown in FIG. 1.
The vehicle 1 forms the socalled controlled system.
The variables given by the driver, namely, the driver brake pressure p<sub>brake </sub>or p<sub>TMC</sub>, the accelerator pedal position, gearshift information, and the steering angle δ, act on the vehicle 1. The variables resulting from this, namely, the motor moment M<sub>motor</sub>, the lateral acceleration a<sub>lat</sub>, the yaw rate {dot over (Ψ)}, the wheel speeds and hydraulic signals, such as wheel brake pressures p<sub>TMC</sub>, are measured on the vehicle, and information about gearshift is determined. To evaluate this data, the DSC system has five electronic controllers 7, 8, 9, 10 and 20, which are associated with the antilock system ABS, the traction slip control system TCS, the electronic brake force distribution system EBD, and the yaw torque control system YTC, respectively, and the sensitive yaw torque control SESP. The electronic controllers for ABS 7, TCS 8, EBD 9, and YTC 10 may basically correspond to the state of the art.
The wheel speeds are sent to the controllers for the antilock system 7, the traction slip control system 8, the electronic brake force distribution system 9, and the sensitive yaw torque control 20. The controller 8 of the traction slip control system additionally receives data on the actual engine torque, the motor moment M<sub>motor</sub>. This information is also sent to the controller 10 for the yaw torque control YTC and the controller 20 for the sensitive yaw torque control SESP. The SESP controller additionally receives gearshift information. In addition, the controllers receive the data on the lateral acceleration a<sub>lat </sub>and the yaw rate {dot over (Ψ)}<sub>meas </sub>of the vehicle from the sensors. Since a vehicle reference velocity v<sub>Ref</sub>, on the basis of which an excess brake slip of one of the wheels can be determined, is determined in the controller 7 of the ABS by way of the individual wheel speeds of the vehicle wheels, such a reference velocity does not need to be calculated in the YTC controller 10, but it is taken over from the ABS controller 7. Where the vehicle reference speed is calculated or whether a separate calculation is performed for the yaw torque control makes only a slight difference for the process of the yaw torque control YTC. This also applies, e.g., to the longitudinal acceleration a<sub>long </sub>of the vehicle. Accordingly, the value for this could also be determined in the ABS controller 7, and sent to the YTC controller 10. This applies to the determination of the coefficient of friction μ of the road surface with restrictions only, because a more accurate coefficient of friction determination than is determined for the antilock system is desirable for yaw torque control and for the sensitive yaw torque control 20.
All five electronic controllers of the DSC, i.e., the controllers for YTC 10, ABS 7, TCS 8 and EBD 9, develop brake pressure set values P<sub>SESP</sub>, P<sub>YTC</sub>, P<sub>ABS</sub>, P<sub>TCS</sub>, P<sub>EBD </sub>for the individual wheels simultaneously and independently from one another based on their own control strategies.
In addition, preset values M<sub>TCS </sub>and M<sub>AdjustM </sub>for the engine torque are calculated in parallel by the TCS controller 8 and the YTC controller 10 and SESP controller.
The pressure preset values P<sub>YTC </sub>and P<sub>SESP </sub>of the YTC controller 10 and the SESP controller 20 for the individual wheel brake pressures are determined as follows:
The YTC controller 10 and the SESP controller 20 first calculate an additional yaw torque MG<sub>YTC </sub>or MG<sub>SESP</sub>, which leads to stabilization and improvement of the driving condition within a curve if it is generated by a corresponding brake actuation. These MG<sub>YTC </sub>and MG<sub>SESP </sub>are sent to a distribution logic unit 2 and 21, respectively, which could also be represented as part of the YTC controller 10 or SESP controller 20, respectively. In addition, the possible desire of the driver to decelerate the vehicle, which is recognized from the driver brake pressure (=P<sub>TMC</sub>), is also sent to distribution logic unit 2 and 20. The distribution logic unit 2 and 21 calculate from the preset yaw torques MG<sub>YTC </sub>and MG<sub>SESP </sub>and from the desired driver brake pressure, yaw torque control brake pressures P<sub>YTC </sub>and P<sub>SESP </sub>for the wheel brakes, which may differ from each other very greatly for the individual wheels. These yaw torque control brake pressures P<sub>YTC </sub>and P<sub>SESP </sub>are sent to a priority circuit 3 for the wheel brake pressures, for function optimization along with the pressure preset values calculated by the other controllers 7, 8 and 9 for ABS, TCS and EBD. This priority circuit 3 determines desired wheel pressures p<sub>Desired </sub>for optimal driving stability and optimal driving behavior, taking into account the driver's desire.
These desired pressures may either correspond to the pressure preset values of one of these four controllers, or represent a superimposition. The procedure followed in the case of the engine torque is similar to the procedure with the wheel brake pressures. While ABS and EBD act only on the wheel brakes, intervention with the engine torque is also provided in the case of SESP, YTC and TCS. The preset values M<sub>SESP</sub>, M<sub>AdjustM </sub>and M<sub>TCS </sub>calculated separately for the engine torque in the ESP controller 20, the YTC controller 10 and in the TCS controller 8 are again evaluated in a priority circuit 4 and superimposed to a desired torque. However, this desired torque M<sub>Desired </sub>may also just as well correspond only to the calculated preset value of one of the three controllers.
Driving stability control by intervention and/or improvement of the driving behavior by way of the brake and engine intervention can now be performed based on the calculated desired preset values for the wheel brake pressure P<sub>Desired </sub>and for the engine torque M<sub>Desired</sub>. Hydraulic signals or values, which reflect the actual wheel brake pressure, are also sent for this purpose to the pressure control unit 5. From this, the pressure control unit 5 generates valve signals, which are sent to the control valves of the individual wheel brakes in the vehicle 1. The engine management 6 controls the drive motor of the vehicle according to M<sub>Desired</sub>, as a result of which a changed motor moment is again generated. This will then again lead to new input variables for the five electronic controllers 7, 8, 9, 10 and 20 of the DSC system.
2. Structure of the Yaw Torque Controller (YTC) including the Sensitive Yaw Torque Controller (SESP)
FIG. 2 shows in a block diagram how the additional yaw torques M<sub>G </sub>and MG<sub>SESP </sub>are determined within the YTC controller 10 and the sensitive yaw torque controller 20 for the distribution logic units 2 and 21. The steering angle δ, the vehicle reference velocity V<sub>Ref </sub>from the ABS controller 7, the measured lateral acceleration a<sub>lat</sub>, and the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>are entered for this purpose as input variables. The vehicle reference velocity V<sub>Ref </sub>passes through a filter 17, which sets a constant value above zero at low velocities, so that the denominator of a fraction will not become equal to zero during the further calculations. The filtered value of V<sub>Ref </sub>is sent to an SESP situation detection system 22 and an activation logic unit 25, while the unfiltered value is sent to an activation logic unit 11, which recognizes the standstill of the vehicle.
This direct determination of the vehicle reference velocity V<sub>Ref </sub>by the activation logic unit 11 may also be omitted if standstill of the vehicle is assumed when the filtered vehicle reference velocity V<sub>RefFil </sub>assumes its constant minimum.
A vehicle reference model 12, which calculates a preset value for a change in the yaw rate Δ{dot over (Ψ)} on the basis of the steering angle δ, the filtered vehicle reference velocity V<sub>RefFil </sub>as well as the measured yaw rate {dot over (Ψ)}<sub>Meas</sub>, is stored in the YTC controller.
To keep the preset values within the physically possible range, the coefficient of friction μ of the road surface,
which is calculated as an estimated value {circumflex over (μ)} in a coefficientoffriction and situation detection system 13, is also needed for these calculations. If the coefficient of friction determined within the framework of the antilock control has sufficient accuracy, this coefficient of friction may be used as well, or the coefficient of friction calculated in the YTC controller 10 may be taken over in the ABS controller 7.
The coefficientoffriction and situation detection system 13 uses for its calculations the filtered reference velocity V<sub>RefFil</sub>, the measurement vehicle lateral acceleration a<sub>lat</sub>, the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>and the steering angle δ. The situation detection system 13 distinguishes different cases, such as straight travel, travel in curves, reverse travel and standstill of the vehicle. Standstill of the vehicle is assumed when the filtered vehicle reference velocity V<sub>RefFil </sub>assumes its constant minimum. Instead of the unfiltered vehicle reference velocity, this information may also be sent to the activation logic unit 11 to recognize standstill of the vehicle. The fact that at a given steering angle δ, the orientation of the measurement yaw rate {dot over (Ψ)} is opposite that during forward travel is utilized to recognize reverse travel. The measured yaw rate {dot over (Ψ)}<sub>Meas </sub> is compared for this purpose with the desired yaw rate {dot over (Ψ)}<sub>desired </sub>preset by the vehicle reference model 12. If the signs are always opposite, and this also applies to the time derivatives of the two curves, the vehicle is traveling in reverse, because {dot over (Ψ)}<sub>desired </sub>is always calculated for forward travel, since the usual speed sensors do not detect information on the direction of rotation of the wheels.
The output signals of the coefficientoffriction and situation detection system 13 are sent to the SESP situation detection system 22 which determines on the basis of the detected driving situations of the YTC controller 10 and further information about the driver's request (δ, M<sub>Motor</sub>, P<sub>TMC</sub>, gearshift information) and the present driving state (v<sub>wheel</sub>, a<sub>lat</sub>, v<sub>ref</sub>, {dot over (Ψ)}<sub>Meas</sub>, 13) whether there is a potential deviation from the course desired by the driver (rotary condition of the vehicle about its vertical axis, in particular turning situation). At least the steering angle δ, the motor moment demand M<sub>Motor</sub>, and the driver brake pressure P<sub>TMC </sub>are included in the detection of the driver's request, while the vehicle state is determined from at least {dot over (Ψ)}<sub>Meas</sub>, v<sub>wheel</sub>, a<sub>lat</sub>, v<sub>ref</sub>. The separate SESP situation detection system 22 safeguards the unambiguousness of the detected driving situation. Making a distinction between at least the states such as cornering with a load change 31, cornering with partial braking 32, straight travel with partial braking 33, and straight travel with load change 34 is necessary in order to adapt the program 23 for realizing the SESP yaw torque control law optimally to the current driving situation. Simultaneously, it is constantly checked by way of conditions such as the YTC state, ABS state, transversal roadway inclination, and vehicle speed whether SESP interventions are permitted. The state 35 (FIG. 30) ‘SESP inactive’ is instantaneously assigned in the event of a negative result. Some causes of turning of the vehicle, such as change of load, naturally have an effect on the vehicle performance which is only limited in time. In these cases, the respective state in the SESP driving situation detection system is left again after a defined time. The temporal limitation of such states permits the driver to correct slow turning maneuvers after a certain ‘time of adaptation’ by himself.
The state detected in the driving situation detection system 22 triggers the activation of the SESP control law unit 23 and induces an optimal adaptation of the parameters (e.g. the control thresholds) to the current driving situation. The driving situation detection system 22 evaluates to this end any causes preceding the vehicle reaction such as P<sub>TMC</sub>, gearshift information, motor moment information so that the sensitive SESP control law unit 23 is activated at an early time. Apart from influencing the SESP controller parameters, the driving situation detection system 22 has the additional objective of indicating to the reference signal formation unit 24 at what moment the offset for the offset correction of the reference variables shall be stored. This is also done already upon detection of any potential cause for a turning situation of the vehicle 1.
Finally, a kinematic velocity of the sideslip angle determination, referred to as kinematic determination 14 for short, is performed in the YTC controller 10 on the basis of the filtered vehicle reference velocity V<sub>RefFil</sub>, the measured vehicle lateral acceleration a<sub>lat</sub>, and the measured yaw rate {dot over (Ψ)}<sub>Meas</sub>.
To cut off peaks in the case of great variations in the sideslip angles, the calculated value of velocity of the sideslip angle passes through a firstorder lowpass filter 15, which sends an estimated value <img id="customcharacter00001" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> for the velocity of the sideslip angle to the activation logic unit 11 and optionally to the SESP activation logic unit 25, to a program 16 for putting the yaw torque control law into practice, and to the SESP reference signal formation unit 24. Program 16 and SESP reference signal formation unit 24 also use the preset values for changing Δ{dot over (Ψ)} for the yaw rate, which is the difference of the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>and the desired yaw rate {dot over (Ψ)}<sub>desired </sub>calculated on the basis of the vehicle reference model 12. The additional yaw torque M<sub>G </sub>for the vehicle, which is to be conveyed by way of the brake pressures, is calculated from this.
The program 16 and the SESP controller 20 operate permanently to keep ready current control variables all of the time. However, whether these controlling torques are transmitted to the distribution logic units 2, 21 shown in FIGS. 1 and 2 depends on the activation logic units 11, 25.
The activation logic unit 11 of the YTC controller 10 receives not only the value of the unfiltered vehicle reference velocity v<sub>Ref </sub>and, as was described, the velocities of the sideslip angles <img id="customcharacter00002" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, but also the amount of the deviation Δ{dot over (Ψ)} of the desired yaw rate {dot over (Ψ)}<sub>desired </sub>from the measured yaw rate {dot over (Ψ)}<sub>Meas</sub>, and information from the situation detection system 13 during reverse travel.
If the vehicle is traveling in reverse, the transmission of M<sub>G </sub>is interrupted. This also happens when standstill of the vehicle is recognized or when neither the estimated velocity of the sideslip angle <img id="customcharacter00003" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> nor the preset value for the change in the yaw rate, Δ{dot over (Ψ)}, reaches a value that would make control necessary. The logic circuit for calculating the engine controlling torque M<sub>AdjustM </sub>is not shown.
The activation logic 25 of the SESP controller 20 not only receives the value of the SESP reference signal formation unit 24 and the activation logic unit 11, but also the deviation <img id="customcharacter00004" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> and/or Δ{dot over (Ψ)}<sub>SESP</sub> of the desired yaw rate {dot over (Ψ)}<sub>Desired </sub>from the measured yaw rate {dot over (Ψ)}<sub>Meas</sub>, but also receives information from the SESP situation detection system 22 when one of the SESP states such as a cornering maneuver with load change, a cornering maneuver with partial braking, straight travel with partial braking, and straight travel with load change prevails or is expected.
If the vehicle has not assumed any one of the determined states, or if states of the YTC controller 10 or the ABS controller 7, the TCS controller 8 are superposed on these states while defined conditions are complied with, the transfer of MG<sub>SESP </sub>is interrupted. The same applies if a standstill of the vehicle is detected or if neither the estimated sideslip angle velocity <img id="customcharacter00005" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, nor the specification for the yaw rate change Δ{dot over (Ψ)}<sub>SESP </sub>reach a value which necessitates a control operation. The logic circuit for calculating the engine controlling torque M<sub>SESPadjustM </sub>is not shown.
2.1 CoefficientofFriction and Situation Detection System of the YTC Controller 10
The logic processes taking place in the coefficientoffriction and situation detection system 13 are shown in the form of flow charts in FIGS. 3, 4 and 5.
FIG. 3 pertains to the situation detection. Eight different driving situations can be distinguished with the process shown:
<0> Standstill of the vehicle
<1> Constant straight travel
<2> Accelerated straight travel
<3> Decelerated straight travel
<6> Reverse travel
<7> Constant travel in curve
<8> Accelerated travel in curve
<9> Decelerated travel in curve.
Logic branches are represented as blocks in the flow chart.
Based on a given situation 51 to be determined, it is first determined in block 52 whether or not the vehicle is at a standstill. If the filtered vehicle reference velocity V<sub>RefFil </sub>assumes its minimum v<sub>min</sub>, standstill of the vehicle, i.e., situation <0>, is assumed. If V<sub>RefFil </sub>is grater than v<sub>min</sub>, the result of the preceding run of situation detection is polled in block 53.
If the situation last determined was recognized as reverse travel, i.e., situation <6>, reverse travel continues to be present, because standstill of the vehicle did not occur in the meantime, because situation <0> would have otherwise been recognized in block 52 in the meantime.
If the preceding run of the situation detection recognized a situation other than <6>, the value of the lateral acceleration a<sub>lat </sub>is polled in block 54. If this is lower than a defined threshold value a<sub>trans min</sub>, it is assumed that the vehicle is traveling straight, i.e., that one of the situations <1> through <3> prevails.
This is also true when the value of the measured lateral acceleration a<sub>lat </sub>is above the threshold value a<sub>latmin</sub>, but it is recognized in block 55 in the next step that the value of the steering angle δ is lower than a threshold value δ<sub>min</sub>. In this case, the measured lateral acceleration a<sub>lat </sub>is an error of measurement that results from the fact that lateral accelerometers are usually securely mounted in the transverse axis of the vehicle and are tilted with the vehicle in the case of an inclination of the road surface, so that a lateral acceleration that does not actually occur is indicated.
Consequently, if the vehicle is traveling straight, the value of the longitudinal acceleration a<sub>long </sub>is examined in block 59. If this value is lower than a threshold value a<sub>longmin</sub>, constant straight travel is assumed. However, if the value of the longitudinal acceleration a<sub>long </sub>is greater than this threshold value, block 60 makes a distinction between positive and negative longitudinal acceleration. If the value of a<sub>long </sub>is above the threshold value a<sub>longmin</sub>, the vehicle is in an accelerated straight travel, i.e., in situation <2>. If the value of a<sub>long </sub>is below the threshold value a<sub>longmin</sub>, this means nothing else but negative longitudinal acceleration, e.g., decelerated straight travel, namely, situation <3>.
If none of the situations <0> through <3> occurs and a steering angle value δ that is greater than the threshold value δ<sub>min </sub>is recognized in block 55, a polling is performed in block 56 to determine whether the vehicle is currently traveling in reverse. The detection of reverse travel is necessary only at this point, because the yaw rate {dot over (Ψ)} hardly differs from zero during straight travel anyway, and no control intervention is therefore performed. Reverse travel must be ruled out with certainty only when travel in a curve is recognized, in which the yaw torque control itself becomes active. This is not possible based solely on the signals of the wheel speed sensors, because such sensors only transmit the value of the speed, without making it possible to infer the direction of travel from it.
As was described above, situation <6> is determined by comparing the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>with the desired yaw rate {dot over (Ψ)}<sub>desired </sub>determined in the vehicle reference model 12. If the signs are opposite, and if this is also true of the time derivatives of the two variables, namely, the yaw accelerations {dot over (Ψ)}<sub>Meas </sub>and {dot over (Ψ)}<sub>desired</sub>, the vehicle is in a curve, traveling in reverse. The signs of the yaw accelerations are therefore compared, to rule out that the opposite signs of the yaw rates originate not only from a phase shift, which is due to the timedelayed calculation of the desired values.
If the conditions for reverse travel are not satisfied, there is travel in a curve in the forward direction. Whether or not this travel in the curve takes place at constant velocity is investigated in block 57. As was done before in blocks 59 and 60 in the case of straight travel, the value of the longitudinal acceleration a<sub>long </sub>is first examined in block 57. If it is lower than the threshold value a<sub>longmin</sub>, there is constant travel in a curve, i.e., situation <7>. In the case of longitudinal acceleration a<sub>long </sub>whose value is greater than the threshold value a<sub>longmin</sub>, it is further examined in block 58 whether the longitudinal acceleration a<sub>long </sub>is positive or negative. The vehicle is in an accelerated travel in a curve, i.e., in situation <8>, in the case of positive longitudinal acceleration a<sub>long</sub>, while a decelerated travel in a curve, corresponding to situation <9>, is recognized in the case of negative longitudinal acceleration a<sub>long</sub>.
The longitudinal acceleration a<sub>long </sub>can be recognized in different ways. It can be determined, e.g., from the reference velocity v<sub>Ref </sub>provided by the ABS controller 7, in which case it should be borne in mind that such a reference velocity v<sub>Ref </sub>may deviate from the actual vehicle velocity during an ABS intervention. Consequently, a correction of v<sub>Ref </sub>is justified in an ABS case. However, the longitudinal acceleration a<sub>long </sub>can also be taken over under certain circumstances directly from the ABS controller if such calculation is performed there.
The situation detection according to FIG. 3 is continually run through again, and the situation last determined remains stored and available in block 53.
A possible process for determining the coefficient of friction of the road surface is shown in FIGS. 4 and 5. According to this process, the coefficient of friction is determined only when the yaw torque controller enters the control. However, since no estimated coefficient of friction is available at the time of entry into the control, the coefficient of friction μ=1 is set at the beginning of the control.
If the yaw torque control system responds on the basis of an instantaneous driving situation, it can be assumed that the vehicle is at least in the vicinity of the borderline range to unstable driving situations. The instantaneous coefficient of friction of the road surface can therefore be inferred from an examination of the current measured variables of the vehicle. The coefficient of friction then determined at the time of entry into the control will subsequently serve as the basis for limiting the desired yaw rate {dot over (Ψ)}<sub>desired </sub>and consequently also for the control deviation for the yaw rate Δ{dot over (Ψ)}, which is transmitted to the YTC control law unit 16. The coefficient of friction is determined for the first time at the time of entry into the control, associated with a subsequent phase of updating for the limitation of the desired yaw rate to physically meaningful values. Based on the originally preset coefficient of friction μ=1, a maximum coefficient of friction {circumflex over (μ)} is determined at the time of entry into the control, and the calculation of the additional yaw torque M<sub>G </sub>will then be based on this value.
An internal coefficient of friction {circumflex over (μ)}<sub>int </sub>is first calculated for this purpose from the measured lateral acceleration a<sub>lat </sub>and a calculated value for the longitudinal acceleration a<sub>long</sub>, which value corresponds to the instantaneous coefficient of friction if complete utilization of the frictional connection is assumed. However, since it must be assumed that the maximum frictional connection is not yet reached at the time of entry into the control, a higher coefficient of friction {circumflex over (μ)} is associated with the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>by means of a table, a characteristic curve or a constant factor. This coefficient of friction {circumflex over (μ)} is then sent to the control system. It is thus possible to calculate with a desired yaw rate {dot over (Ψ)}<sub>desired </sub>adjusted to the coefficient of friction of the road surface in the next step of the calculation and to improve the control. While the yaw torque control system acts on the vehicle, the estimated coefficient of friction {circumflex over (μ)} must be further updated, because a change in the coefficient of friction might take place during the control. If the control system is not activated based on the adjustment of the coefficient of friction in the vehicle reference model due to the resulting changed control deviation of the yaw rate Δ{dot over (Ψ)}, the coefficient of friction {circumflex over (μ)} is further updated in Tμ<sub>End </sub>number of steps. If the yaw torque control system is not activated even during this phase of updating, the estimated coefficient of friction {circumflex over (μ)} is reset to 1.
The adjustment or updating of the estimated coefficient of friction {circumflex over (μ)} may also be omitted in certain situations. Such situations are, e.g., straight travel, travel in reverse or standstill of the vehicle, i.e., situations <0> through <4>. These are situations in which no yaw torque control is performed anyway, so that an estimation of the coefficient of friction may be omitted if the time derivative of the coefficient of friction {circumflex over (μ)}, i.e., <img id="customcharacter00006" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00001.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, is negative and the value of the time derivative of the steering angle δ, i.e. {dot over (δ)}, exceeds a predetermined threshold. It can be assumed in the latter case that a change in the lateral acceleration a<sub>lat </sub>is based on a change in the steering angle δ, rather than on a change in the coefficient of friction.
It is generally true of the coefficient of friction calculated in this manner that it is a mean coefficient of friction for all four wheels of the vehicle. The coefficient of friction cannot be determined in this manner for the individual wheels.
The process of the coefficientoffriction determination will now be explained on the basis of FIG. 4. The behavior of the vehicle is affected by the prevailing coefficient of friction of the road surface according to field 61 in each driving situation. To determine the corresponding coefficient of friction of the road surface, the measured lateral acceleration a<sub>lat </sub>is first filtered according to step 62, i.e., either the measured values are smoothed, or the curve passes through a lowpass filter, so that no extreme peaks appear. Step 63 comprises the situation detection according to FIG. 3. The driving situation recognized is later significant for the phase of updating in step 74. Polling is performed in block 64 to determine whether a control intervention is necessary. Such a calculation is first based on the initial coefficient of friction μ=1. If control is considered to be necessary, polling is performed in block 65 to determine whether this was also the condition at the end of the preceding run of the coefficientoffriction determination. If an entry into control is involved here, control was not recognized before, so that an internal coefficient of friction {circumflex over (μ)}<sub>int </sub>is determined for the first time in step 67. It is calculated from the following equation:<maths id="MATHUS00003" num="3"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mover><mi>μ</mi><mo>^</mo></mover><mi>int</mi></msub><mo>=</mo><mfrac><msqrt><mrow><msubsup><mi>a</mi><mi>trans</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mi>long</mi><mn>2</mn></msubsup></mrow></msqrt><mi>g</mi></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.1</mn></mrow></mtd></mtr></mtable></math></maths>
Here, g is the gravitation constant g=9.81 m/sec<sup>2</sup>.
The parameter reg<sub>old </sub>for step 65 is set at 1 in step 68. In addition, the counting parameter Tμ is set at 1, corresponding to the fact that the first determination of the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>has been performed. An estimated coefficient of friction {circumflex over (μ)} is associated with the calculated internal coefficient of friction {circumflex over (μ)}<sub>int </sub>in step 69. This is done under the assumption that the existing acceleration components are not based on a complete utilization of the frictional connection, either. The estimated coefficient of friction {circumflex over (μ)} is consequently usually between the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>thus determined and 1. The determination of the coefficient of friction is thus concluded.
Consequently, assuming an unchanged driving situation, reg<sub>old</sub>=1 is decided in block 65 during the next run of this coefficientoffriction determination. A value {circumflex over (μ)}<sub>int</sub>, which replaces {circumflex over (μ)}<sub>int</sub>, determined in the preceding run, is later determined here as well. The parameters determined in field 68 are not updated, because the updating of {circumflex over (μ)}<sub>int </sub>was performed during a control. reg<sub>old </sub>had been set at 1 already in the run before that, and it remains unchanged. The number Tμ of runs performed continues to be 1, because counting is continued only if no control takes place. As was described above, an estimated coefficient of friction {circumflex over (μ)} is also associated with the updated value of {circumflex over (μ)}<sub>int </sub>by means of a table, a nonlinear relationship, or a constant factor.
If it is determined in one run in block 64 that control is not necessary, a polling is then performed in block 71 to determine whether the parameter reg<sub>old </sub>for the control was last set at 0 or 1. If it was set at 1 in the last run, the number Tμ of runs is polled in block 72. This Tμ equals 1 if control was performed in the last run. If control was performed only in the run before last, Tμ=2, etc. If Tμ has not yet reached a certain TμEnd in step 72, it is increased by 1 in step 73, and a repeated updating of the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>is performed in step 74. If the number TμEnd is then reached in one of the next runs without control having taken place, the parameter reg<sub>old </sub>is again reset to 0 for the control. The estimated coefficient of friction {circumflex over (μ)} is equated with the initial coefficient of friction μ=1. The phase of updating for the coefficient of friction {circumflex over (μ)} is thus terminated.
If it is then again recognized in the next run in block 64 that no control is necessary, the initial coefficient of friction μ=1 is retained in field 76 in block 71 with reg<sub>old</sub>=0. A coefficientoffriction determination is again performed only if the necessity of a control intervention is recognized in block 64.
The criteria for updating the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>after step 74 are shown in FIG. 5. Based on the instruction in field 77 that the internal coefficient of friction {circumflex over (μ)}<sub>int </sub>is to be updated, the time derivatives of the estimated coefficients of friction {circumflex over (μ)} or {circumflex over (μ)}<sub>int </sub>formed before, as well as of the steering angle δ are formed in step 78.
When it is then recognized in block 79 that the vehicle is neither at a standstill nor is it traveling straight, i.e., that one of the situations <6> through <9> occurs, the results from step 78 are evaluated in step 80. A coefficient of friction determination is performed, as was explained above, only if a decreasing coefficient of friction cannot be attributed to a steering maneuver. No updating of the coefficient of friction is performed if the vehicle is traveling straight, forward or in reverse, or if it is at a standstill, or if a reduction in the estimated coefficient of friction {circumflex over (μ)} can be attributed to a steering maneuver.
2.1.1 SESP Situation Detection System
The SESP situation detection system 22 is designed as a state machine. The state machine detects the state (the driving situations) of the vehicle by way of the quantity of the attribute values (input signals). State diagrams indicate at which event the vehicle passes over from a defined state into a neighboring state. A subsequent state depends on the initial state and on the event that occurred. FIG. 30 shows as an example the possible states, namely, cornering maneuver with load change 31, cornering maneuver with partial braking 32, straight travel with partial braking 33, and straight travel with load change 34, and the permitted state transitions in the SESP situation detection system 22. The use of a state machine safeguards the unambiguousness of the detected driving situation. Making a distinction between several states is necessary in order to adapt the programs 23 for realizing the SESP yaw torque control law optimally to the current driving situation. All driving situations are determined in the SESP situation detection system 22 which are relevant for the calculation of the SESP driving situations states 31 to 34. To this end, the situation detection system 22 determines the following driving situations which can then be used as input information in the state machine 22.1.
SteadyState Straight Ahead Driving Maneuver
The SESP situation detection system 22 detects if the driver wants to drive straight ahead and uses the driving situation <1>, <2> and <3> detected by the coefficientoffriction and situation detection system 13 and the steering angle δ, i.e. the driver's steering input. The result is memorized in flag Sesp_straight_ahead.
Flag Sesp_straight_ahead is set to True if all of the following conditions are fulfilled:
 i. The driving situation from 13 is <1>, <2> or <3>
 ii. The amount of the steering angle δ<threshold value k1
 iii. The amount of the steering angle velocity {dot over (δ)}<threshold value k2
Otherwise flag Sesp<sub>—straight</sub>_ahead is set to False.
SteadyState Curve Driving
The SESP situation detection system 22 detects if the driver wants to drive a steady curve, using the driving situation <7>, <8> and <9> detected by the coefficientoffriction and situation detection system 13, the lateral acceleration, the driver's steering input δ and the roadway radius estimated from a<sub>lat </sub>and v<sub>ref </sub>in the coefficientoffriction and situation detection system 13 or in the YTC controller 10, respectively. The result is stored in the flag.
Flag Sesp_steady_curve is set to True if all of the following conditions are fulfilled for a predetermined duration (counter):
 iv. The driving situation from 13 is <7>, <8> or <9>
 v. The estimated curve radius (from 10, 13)>threshold value k3
 vi. The amount of the steering angle velocity {dot over (δ)}<threshold value k4
If any one of the conditions is not fulfilled, flag Sesp_steady_curve is set to False; the counter is reset to 0.
Application of the Brakes Possible
The SESP situation detection system 22 detects if the driver is braking in a specific manner, involving the risk that ‘brake pull (deviation from the desired track)’ of the vehicle occurs because the driver is braking, and the information whether the driver is braking in a specific manner is determined using the master cylinder pressure P<sub>TMC </sub>and its gradient. The result is stored in flag Sesp_brake_pull_possible.
Flag Sesp_brake_pull_possible is set to True if all of the following conditions are fulfilled:
 vii. The YTC controller 10 detects braking by the driver
 viii. The driver brake pressure P<sub>TMC</sub>>threshold value k5
Otherwise flag Sesp<sub>—brake</sub>_pull_possible is set to False.
Turning Possible because the Driver is Braking
The SESP situation detection system 22 detects if there is an oversteering tendency (turning tendency) in curves because the driver is braking in a specific manner, which can cause turning of the vehicle into the bend, and the information whether the driver is braking in this specific manner is determined using the master cylinder pressure P<sub>TMC </sub>and its gradient. The result is stored in flag Sesp_brake_ov_possible.
Flag Sesp_brake_ov_possible is set to True if all of the following conditions are fulfilled:
 ix. The YTC controller 10 detects braking by the driver
 x. The gradient of P<sub>driver</sub>>threshold value k6; the condition x. is maintained within a predetermined time window (period) because after detecting a turning maneuver by the driver, it is still possible during a defined period that the vehicle turns into the bend.
 xi. The driver brake pressure P<sub>TMC</sub>>threshold value k7
Otherwise flag Sesp<sub>—brake</sub>_ov_possible is set to False.
Oversteering due to Engine Drag Torque Possible
The SESP situation detection system 22 detects, if an engine drag torque or engine brake torque exists, respectively, which can cause turning of the vehicle into a bend. The result is stored in flag Sesp_drag_ov_possible.
Flag Sesp_drag_ov_possible is set to True if the following conditions are fulfilled:
 xii. M<sub>motor</sub><threshold value k8
 and simultaneously one of the following conditions:
 xiii. Sum of the wheel speeds vwheel of the nondriven axle−sum of the wheel speeds vwheel of the driven axle<threshold value k9; the signal of the speed difference of the wheel speeds is filtered in a lowpass filter of first order. The speed difference of the wheels indicates a load change;
 or
 xiv. gradient of the engine torque M<sub>motor</sub><negative threshold value k10; the condition xiv. must be satisfied within a predetermined time window (period) in retrospect, if the condition xii. is currently determined;
 or
 xv. a gear shift takes place or has taken place, respectively.
Otherwise flag Sesp_drag_ov_possible is set to False after a predetermined period. The flag is reset with delay because the detected load changes will still have consequences on the driving performance of the vehicle during a defined period after their detection, and may cause turning of the vehicle into a bend.
Detection of a Turning Tendency
The SESP situation detection system 22 detects if the vehicle tends to turning into a bend during straight driving or cornering, with the yaw rate and its acceleration being used. The result is stored in flag Sesp_oversteer_tendency.
Flag Sesp_oversteer_tendency is set to True, if all of the following conditions are fulfilled:
 xvi. {umlaut over (ψ)}<sub>meas </sub>(=derivative of {dot over (Ψ)}<sub>Meas</sub>)>threshold value k11
 xvii. Sign of {umlaut over (ψ)}<sub>Meas</sub>=the sign of {dot over (ψ)}<sub>Meas </sub>
Otherwise flag Sesp_oversteer_tendency is set to False.
These conditions are sufficient for the SESP because SESP is limited to situations in which the driver desires a yaw rate that remains almost constant. SESP reference values are not used in this case because they are not always available.
ABS active at Rear Axle
The SESP situation detection system 22 detects, if at least one wheel of the rear axle is ABScontrolled. The result is stored in flag Sesp_abs_active_at_ra.
Flag Sesp_abs_active_at_ra is set to True, if one of the following conditions is fulfilled:
 xviii. Right rearwheel is ABScontrolled
 or
 xix. Left rear wheel is ABScontrolled
Otherwise flag Sesp_abs_active_at_ra is set to False.
SESP Driving State
The SESP situation detection system 22 calculates the SESP driving situation state SESP_DRIVE_STATE. This is done in the state machine 22.1 (FIG. 30) using the results of the individual abovementioned SESP detections of the driving situations and of the SESP deactivation control system as input. SESP_DRIVE_STATE is used in order to adapt the SESP control to the driving situation.
The following state transitions are fixed in connection with FIG. 30 for SESP_DRIVE_STATE:<tables id="TABLEUS00001" num="1"><table frame="none" colsep="0" rowsep="0"><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><thead><row><entry/><entry/></row><row><entry/><entry namest="OFFSET" nameend="1" align="center" rowsep="1"/></row></thead><tbody valign="top"><row><entry/><entry>Transition from state 35 to state 33.</entry></row><row><entry/><entry>Sesp_straight_braked:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_straight_ahead == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_brake_pull_possible == True</entry></row><row><entry/><entry> AND</entry><entry>Abs_cycle == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_f OR bidden == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 35 to state 32</entry></row><row><entry/><entry>Sesp_curve_braked:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_steady_curve == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_brake_ov_possible == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_oversteer_tendency == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_abs_active_at_ra == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_f OR bidden == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 35 to state 31</entry></row><row><entry/><entry>Sesp_curve_drag_tq:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_steady_curve( ) == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_drag_ov_possible( ) == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_oversteer_tendency( ) == True</entry></row><row><entry/><entry> AND</entry><entry>Ayc_driver_braking( ) == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_f OR bidden( ) == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 35 to state 34</entry></row><row><entry/><entry>Transition from state 33 to state 35</entry></row><row><entry/><entry>Sesp_straight_braked</entry></row><row><entry/><entry>Sesp_drive_idle:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_straight_ahead == False</entry></row><row><entry/><entry> OR</entry><entry>Abs_cycle == True</entry></row><row><entry/><entry> OR</entry><entry>Sesp_f OR bidden == True</entry></row><row><entry/><entry> OR</entry><entry>Sesp_in_cycle == False</entry></row><row><entry/><entry> AND</entry><entry>(Sesp_brake_pull_possible == False)</entry></row><row><entry/><entry/><entry> )</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 32 to state 31</entry></row><row><entry/><entry>Sesp_curve_braked</entry></row><row><entry/><entry>Sesp_curve_drag_tq:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_in_cycle == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_steady_curve == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_drag_ov_possible == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_brake_ov_possible == False</entry></row><row><entry/><entry> AND</entry><entry>Ayc_driver_braking == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_f OR bidden == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 32 to state 35</entry></row><row><entry/><entry>Sesp_drive_idle:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_steady_curve == False</entry></row><row><entry/><entry> OR</entry><entry>Sesp_brake_ov_possible == False</entry></row><row><entry/><entry> OR</entry><entry>Sesp_abs_active_at_ra == True</entry></row><row><entry/><entry> OR</entry><entry>Sesp_f OR bidden == True</entry></row><row><entry/><entry> OR</entry><entry>Sesp_in_cycle == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_oversteer_tendency == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 31 to state 32</entry></row><row><entry/><entry>Sesp_curve_drag_tq</entry></row><row><entry/><entry>Sesp_curve_braked:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_in_cycle == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_steady_curve == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_brake_ov_possible == True</entry></row><row><entry/><entry> AND</entry><entry>Sesp_abs_active_at_ra == False</entry></row><row><entry/><entry> AND</entry><entry>Sesp_f OR bidden == False</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="189PT" align="left"/><tbody valign="top"><row><entry/><entry>Transition from state 31 to state 35</entry></row><row><entry/><entry>Sesp_drive_idle:</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="3"><colspec colname="OFFSET" colwidth="28PT" align="left"/><colspec colname="1" colwidth="35PT" align="left"/><colspec colname="2" colwidth="154PT" align="left"/><tbody valign="top"><row><entry/><entry/><entry>Sesp_steady_curve( ) == False</entry></row><row><entry/><entry> OR</entry><entry>Sesp_drag_ov_possible == False</entry></row><row><entry/><entry> OR</entry><entry>Sesp_f OR bidden == True</entry></row><row><entry/><entry> OR</entry><entry>(Sesp_in_cycle == False</entry></row><row><entry/><entry> AND</entry><entry>(Ayc_driver_braking == True</entry></row><row><entry/><entry> OR</entry><entry>Sesp_oversteer_tendency == False))</entry></row><row><entry/><entry namest="OFFSET" nameend="2" align="center" rowsep="1"/></row></tbody></tgroup></table></tables>
Outside SESP control, a transition always passes over the state Sesp_drive_idle 35. Within SESP control, the transitions pass e.g. between the two states Sesp_curve_braked 32 and Sesp_curve_drag_tq 31 (or between the states 34 and 33), yet do not pass the state Sesp_drive_idle 35, if the control shall be continued in the new state. This is to avoid a reset of the SESP reference yaw rate 24.
If the conditions for the states Sesp_curve_braked 32 and Sesp_drag_tq 31 are fulfilled at the same time, the state with the partial braking, e.g. Sesp_curve_braked 32, always has the higher priority compared to the state with a load change.
2.1.1.1 Common Inputs for all SESP Controllers 20, 23
The calculation of signals is performed by the YTC controller 10 and can be commonly used by the SESP controllers as input.
Minimum Friction:
Estimates the road friction, using the longitudinal acceleration of the vehicle 1 and the friction signal that is calculated by the coefficientoffriction and situation detection system 13. The result is stored in SESP_MY_MIN.
In the majority of situations being relevant for SESP control, the vehicle does not fully use the available road friction. Therefore, SESP_MY_MIN represents only that part of the road friction which is normally utilized by vehicle 1.
2.2. Determination of {dot over (β)} and <img id="customcharacter00007" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>
The prevailing sideslip angle β as well as its time derivative, the velocity of the sideslip angle {dot over (β)}, are an indicator of the stability of a driving condition. The determination of these values will be explained below.
2.2.1 Kinematic Determination of {dot over (β)}
The kinematic determination of {dot over (β)}, 14, is nothing else but the determination of the velocity of the sideslip angle {dot over (β)}, separated from any vehicle models, from measured variables or from variables calculated on the basis of measured values, according to purely physical considerations:
The acceleration a<sub>trans </sub>of the center of gravity of the vehicle at right angles to its longitudinal axis in the plane of movement is measured. The center of gravity of the vehicle moves with the velocity vector <u style="single">v</u> relative to an inertial system:<maths id="MATHUS00004" num="4"><math overflow="scroll"><mtable><mtr><mtd><mrow><munder><mi>v</mi><mi>_</mi></munder><mo>=</mo><mrow><mi>v</mi><mo></mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mi>cos</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>sin</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.2</mn></mrow></mtd></mtr></mtable></math></maths>
The yaw angle is designated by Ψ and the sideslip angle by β. The acceleration vector <u style="single">a</u> is obtained as a derivative over time t as:<maths id="MATHUS00005" num="5"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><munder><mi>a</mi><mi>_</mi></munder><mo>=</mo><mrow><mfrac><mo>ⅆ</mo><mrow><mo>ⅆ</mo><mi>t</mi></mrow></mfrac><mo></mo><munder><mi>v</mi><mi>_</mi></munder></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mrow><mrow><mover><mi>v</mi><mo>.</mo></mover><mo></mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mi>cos</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>sin</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>+</mo><mrow><mrow><mi>v</mi><mo></mo><mrow><mo>(</mo><mrow><mover><mi>Ψ</mi><mo>.</mo></mover><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mo></mo><mrow><mi>sin</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>cos</mi><mo></mo><mrow><mo>(</mo><mrow><mi>Ψ</mi><mo>+</mo><mi>β</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.3</mn></mrow></mtd></mtr></mtable></math></maths>
The acceleration sensor measures the projection of the acceleration vector to the transverse axis of the vehicle:<maths id="MATHUS00006" num="6"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>a</mi><mi>trans</mi></msub><mo>=</mo><mrow><msup><munder><mi>a</mi><mi>_</mi></munder><mi>r</mi></msup><mo></mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mrow><mo></mo><mi>sin</mi></mrow><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Ψ</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>cos</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Ψ</mi></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.4</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mi>trans</mi></msub><mo>=</mo><mrow><mrow><mover><mi>v</mi><mo>.</mo></mover><mo></mo><mi>sin</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>β</mi></mrow><mo>+</mo><mrow><mrow><mi>v</mi><mo>·</mo><mrow><mo>(</mo><mrow><mover><mi>Ψ</mi><mo>.</mo></mover><mo>+</mo><mover><mi>β</mi><mo>.</mo></mover></mrow><mo>)</mo></mrow></mrow><mo></mo><mi>cos</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>β</mi></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.5</mn></mrow></mtd></mtr></mtable></math></maths>
After linearization of the trigonometric functions (sin β=β; cos β=1), the equation can be rewritten as<maths id="MATHUS00007" num="7"><math overflow="scroll"><mtable><mtr><mtd><mrow><mover><mi>β</mi><mo>.</mo></mover><mo>=</mo><mrow><mfrac><msub><mi>a</mi><mi>trans</mi></msub><mi>v</mi></mfrac><mo></mo><mover><mi>Ψ</mi><mo>.</mo></mover><mo></mo><mrow><mfrac><mover><mi>v</mi><mo>.</mo></mover><mi>v</mi></mfrac><mo></mo><mi>β</mi></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.6</mn></mrow></mtd></mtr></mtable></math></maths>
The velocity of the sideslip angle {dot over (β)} can now be calculated corresponding to the above differential equation. Besides the lateral acceleration a<sub>lat</sub>, the yaw rate Ψ, the scalar velocity of the vehicle v and its time derivative &vdot; are included as measured variables. To determine β, {dot over (β)} from the previous calculation can be numerically integrated, and &vdot;=0 is assumed for the first determination of {dot over (β)}. A simplification is obtained if the last term is generally ignored, so that no β needs to be determined.
The proposed procedure offers the advantage that the velocity of the sideslip angle {dot over (β)} is directly derived from the sensor signals, and thus it can also be determined in the nonlinear range of the transverse dynamics. The disadvantages are the sensitivity of the procedure to measurement noise and the cumulative integration of errors of measurement, as a result of which the determination of the sideslip angle may become highly inaccurate.
These disadvantages are circumvented by the combination with a modelsupported procedure. FIG. 6, which can be inserted in place of the block 18 drawn in broken line in FIG. 2, shows such a combination of the kinematic determination with the observer modelsupported determination of the velocity of the sideslip angle {dot over (β)}. The steering angle δ, which is indicated by an arrow drawn in broken line, is also included as an additional input variable in such a modelsupported procedure. The mutual influences and correction of the combined methods of determination of the velocity of the sideslip angle <img id="customcharacter00008" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> also make it possible to calculate the sideslip angle β itself with less error, so that it can then also be made available to the control as {circumflex over (β)}. This is also indicated by an arrow drawn in broken line.
2.2.2. Combination of the Kinematic Determination of {dot over (β)} with an Observer Vehicle Model
The area 18 bordered in broken lines in FIG. 2 can also be replaced with the representation according to FIG. 6. It will thus become possible to determine not only the existing velocity of the sideslip angle {dot over (β)}, but also the prevailing sideslip angle β. Contrary to a purely kinematic calculation of the velocity of the sideslip angle <img id="customcharacter00009" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, an observer vehicle model 84 is used here to determine the driving condition, in addition to the kinematic determination of {dot over (β)}, i.e. 83. Just like the vehicle reference model 12 for determining the yaw rate, the observer vehicle model 84 receives the steering angle δ as the input variable. The filtered vehicle reference velocity v<sub>RefFil </sub>is included as a parameter. The measurable output variables, namely, the lateral acceleration a<sub>lat </sub>and the yaw rate {dot over (Ψ)}<sub>Meas</sub>, are needed for the kinematic determination of {circumflex over (β)}, 83, but not for the observer vehicle model 84, which creates these variables, in principle, itself. Another term Y, which is identical in the simplest case to the additional yaw torque calculated by the yaw torque control law unit, represents the changes in the vehicle behavior, which are caused by a control intervention. Y is hence used to expose the observer's simulated vehicle to the same conditions as the real vehicle.
Besides a velocity of the sideslip angle {dot over (β)}<sub>Obs</sub>, the observer vehicle model also gives a value for the yaw acceleration {umlaut over (Ψ)}<sub>Obs</sub>. The variable for the velocity of the sideslip angle <img id="customcharacter00010" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, which originates from the kinematic determination of {dot over (β)}, is multiplied by a weighting factor k after passing through the lowpass filter, while the variable for the velocity of the sideslip angle {dot over (β)}<sub>Obs </sub>Y after addition with a correction factor from the measured yaw angle velocity, is multiplied by a factor h determining the degree of the correction, i.e. a weighting factor (1−k). The value of k is always between 0 and 1. We would have k=1 without the observer vehicle model. After adding the two velocities of the sideslip angles, the sum is integrated into an estimated sideslip angle {circumflex over (β)}. Besides the kinematic velocity of the sideslip angle <img id="customcharacter00011" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>, this is also made available to the control. In addition, the sideslip angle {circumflex over (β)} is transmitted to both the kinematic determination of β, 83, and the observer vehicle model 84. A similar correcting variable is the yaw acceleration {umlaut over (Ψ)}Obs calculated by the observer vehicle model 84.
This is first integrated to a yaw rate and returns to the observer vehicle model 84, on the one hand, and is subtracted from the measured yaw rate {dot over (Ψ)}meas, on the other hand. This difference is multiplied by a factor h<sub>2 </sub>that determines the value of the next control steps in the correction of the observer vehicle model 84 and is provided with the dimension 1/s. The yaw rate multiplied by this factor h<sub>2 </sub>has consequently the same dimension as the yaw acceleration {umlaut over (Ψ)}, so that the two variables can be added up and form a returning correcting variable for the yaw rate after further integration. In the course of a yaw torque control, the term Y assumes values different from zero, corresponding to the additional yaw torque M<sub>G </sub>applied. By being divided by the moment of inertia in yaw 0 of the vehicle, the term Y also acquires the dimension of yaw acceleration and is added to the sum of the yaw accelerations, so that the integrated correction variable also takes into account the control effects or influences.
If an observer vehicle model 84 according to FIG. 6 prevails, which makes possible a more reliable determination of the sideslip angle β than would be possible with a purely kinematic determination of the velocity of the sideslip angle <img id="customcharacter00012" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> and integration, the sideslip angle {circumflex over (β)} thus determined can also be transmitted to the actual yaw torque controller 10.
The kinematic determination of {dot over (β)}, which takes place in combination with an observer vehicle model, is shown in FIG. 7. As is apparent already from FIG. 6, the lateral acceleration a<sub>lat </sub>and the yaw rate {dot over (Ψ)}<sub>Meas </sub>are included in the calculation 91 according to Equation F 2.6 as measured output variables.
The filtered vehicle reference velocity V<sub>RefFil </sub>is differentiated in field 93 to provide the vehicle reference velocity &vdot;<sub>Ref</sub>, which is divided in field 94 by the filtered vehicle reference velocity V<sub>RefFil</sub>, which leads to a factor fβ after nonlinear multiplication 95. This nonlinear multiplication 95 leads to the factor fβ being set to equal zero at low quotients of &vdot;<sub>Ref </sub>and v<sub>RefFil </sub>so that this factor, which precedes the sideslip angle {circumflex over (β)}, can be ignored. The sideslip angle {circumflex over (β)} is taken into account in the kinematic determination of {dot over (β)} only when the vehicle acceleration &vdot;<sub>Ref </sub>reaches a significant value. The {circumflex over (β)} used here is the combined {circumflex over (β)}, which is used both as a variable for the control and for feedback according to FIG. 6. After calculation 91, the value determined for the velocity of the sideslip angle passes through a lowpass filter 92, as was described above, and it yields the estimated velocity of the sideslip angle <img id="customcharacter00013" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>.
The filtered vehicle reference velocity v<sub>RefFil </sub>is differentiated in field 93 to provide the vehicle reference velocity &vdot;<sub>Ref</sub>, which is divided in field 94 by the filtered vehicle reference velocity v<sub>RefFil</sub>, which leads to a factor fβ after nonlinear multiplication 95. This nonlinear multiplication 95 leads to the factor fβ being set to equal zero at low quotients of &vdot;<sub>Ref </sub>and v<sub>RefFil </sub>so that this factor, which precedes the sideslip angle {circumflex over (β)}, can be ignored. The sideslip angle {circumflex over (β)} is taken into account in the kinematic determination of {dot over (β)} only when the vehicle acceleration &vdot;<sub>Ref </sub>reaches a significant value. The {circumflex over (β)} used herein is the combined {circumflex over (β)}, which is used both as a variable for the control and for feedback according to FIG. 6. After calculation 91, the value determined for the velocity of the sideslip angle passes through a lowpass filter 92, as was described above, and it yields the estimated velocity of the sideslip angle <img id="customcharacter00014" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>.
FIG. 8 shows how the observer vehicle model 84 from FIG. 6 operates. A matrix representation was selected, in which “→” are scalar and “=>” multidimensional formations.
The matrix representation is based on Equations F 1.1 through F 1.3 The phase variables β and {dot over (Ψ)} are combined into a phase vector <u style="single">x</u>(t), so that the following set of equations is obtained:<FORM><u style="single">{dot over (x</u>)}(t)=<u style="single">A</u>(v(t))<u style="single">x</u>(t)+<u style="single">B</u>(v(t))<u style="single">u</u>(t) F 2.7 </FORM>with the system matrix <u style="single">A</u>(v(t)), the input matrix B(v(t)), the phase vector <u style="single">x</u>(t) and the input vector <u style="single">u</u>(t):<maths id="MATHUS00008" num="8"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><munder><mi>A</mi><mi>_</mi></munder><mo></mo><mrow><mo>(</mo><mrow><mi>v</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mo></mo><mfrac><mrow><msub><mi>c</mi><mi>h</mi></msub><mo>+</mo><msub><mi>c</mi><mi>v</mi></msub></mrow><mrow><mi>mv</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd><mtd><mrow><mrow><mo></mo><mn>1</mn></mrow><mo>+</mo><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow><mrow><msup><mi>mv</mi><mn>2</mn></msup><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow><mi>Θ</mi></mfrac></mtd><mtd><mrow><mo></mo><mfrac><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo></mo><msubsup><mi>l</mi><mi>h</mi><mn>2</mn></msubsup></mrow><mo>+</mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msubsup><mi>l</mi><mi>v</mi><mn>2</mn></msubsup></mrow></mrow><mrow><mi>Θ</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mrow><mi>v</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow></mfrac></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><mrow><munder><mi>B</mi><mi>_</mi></munder><mo></mo><mrow><mo>(</mo><mrow><mi>v</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mfrac><msub><mi>c</mi><mi>v</mi></msub><mrow><mi>mv</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>c</mi><mi>v</mi></msub><mo></mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mi>Θ</mi></mfrac></mtd><mtd><mfrac><mn>1</mn><mi>Θ</mi></mfrac></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>;</mo></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><mrow><munder><mi>x</mi><mi>_</mi></munder><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mi>β</mi><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>Ψ</mi><mo>.</mo></mover><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>;</mo></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><munder><mi>u</mi><mi>_</mi></munder><mo></mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mi>δ</mi></mtd></mtr><mtr><mtd><mi>Y</mi></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.8</mn></mrow></mtd></mtr></mtable></math></maths>
The input vector <u style="single">u</u>(t) contains as the input variables the steering angle δ and the term Y, which is the additional yaw torque generated by the yaw torque control system.
Instead of weighting factors, a weighting matrix <u style="single">K</u><sub>1 </sub>and a weighting vector <u style="single">k</u><sub><u style="single">2</u></sub><u style="single"></u> are used for the weighted addition of the variables determined.<maths id="MATHUS00009" num="9"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><msub><munder><mi>K</mi><mi>_</mi></munder><mn>1</mn></msub><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><mn>1</mn><mo></mo><mi>k</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>;</mo></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><munder><msub><mi>k</mi><mn>2</mn></msub><mi>_</mi></munder><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mi>k</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>;</mo></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mn>0</mn><mo>≤</mo><mi>k</mi><mo>≤</mo><mn>1</mn></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.9</mn></mrow></mtd></mtr></mtable></math></maths>
To eliminate the process variables, two vectors, <u style="single">c</u>β and c{dot over (Ψ)}, are introduced, with cancel one component of phase vector x(t) each:<maths id="MATHUS00010" num="10"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><munder><msub><mi>c</mi><mi>β</mi></msub><mi>_</mi></munder><mo>=</mo><mrow><mo>[</mo><mrow><mn>1</mn><mo>,</mo><mn>0</mn></mrow><mo>]</mo></mrow></mrow><mo>;</mo></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><mi>c</mi><mo></mo><mover><mi>Ψ</mi><mo>.</mo></mover></mrow><mo>=</mo><mrow><mo>[</mo><mrow><mn>0</mn><mo>,</mo><mn>1</mn></mrow><mo>]</mo></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.10</mn></mrow></mtd></mtr></mtable></math></maths>
The dynamics of the observer vehicle model, i.e., the value of the correction steps, is determined by a vector <u style="single">h</u>, whose first component, h<sub>1</sub>, is dimensionless, and whose second component, h<sub>2</sub>, has the dimension (1/s):<maths id="MATHUS00011" num="11"><math overflow="scroll"><mtable><mtr><mtd><mrow><munder><mi>h</mi><mi>_</mi></munder><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><msub><mi>h</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>h</mi><mn>2</mn></msub></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.11</mn></mrow></mtd></mtr></mtable></math></maths>
Based on the vehicle model in the description of the phase space (F 1.1 and F 1.2), the structure described below is then obtained for determining the sideslip angle β by means of an observer according to FIG. 8.
The vehicle 101 is shown in FIG. 8 only to distinguish between input variables and output variables. It is not a part of the combined procedure for determining the velocity of the sideslip angle <img id="customcharacter00015" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>.
The system equations according to F 2.7 are formed in the adder 104. To do so, the system matrix <u style="single">A</u> is multiplied by the phase vector <u style="single">x</u>, and the input matrix <u style="single">B</u> is multiplied by the input variables δ and y, i.e., with the input vector <u style="single">u</u>.
The current vehicle reference velocity v<sub>RefFil </sub>is included as the only variable parameter in both the system matrix <u style="single">A</u> and the input matrix <u style="single">B</u>. The time derivative &xdot; of the phase vector <u style="single">x</u>, formed in the adder 104 by addition, is now multiplied by the weighting matrix<maths id="MATHUS00012" num="12"><math overflow="scroll"><munder><msub><mi>K</mi><mn>1</mn></msub><mi>_</mi></munder></math></maths>according to F 2.9 and is sent to another adder 105.
Simultaneously to these processes, a velocity of the sideslip angle <img id="customcharacter00016" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is estimated in the direct procedure 103. The filtered vehicle reference velocity v<sub>RefFil</sub>, as well as its time derivative &vdot;<sub>Ref</sub>, determined in the differentiator 102 (identified by 93 in FIG. 7), the measured lateral acceleration a<sub>lat</sub>, as well as the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>according to Equation F 2.6 are used for this. The last term of the equation is ignored in the first step, because no value of the sideslip angle β is available as yet. After the velocity of the sideslip angle is determined, it still passes through the lowpass filter 92, as was shown in FIG. 7, after which the resulting estimated velocity of the sideslip angle <img id="customcharacter00017" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is made available for the further calculation. This <img id="customcharacter00018" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> corresponds to the <img id="customcharacter00019" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> which is output from the shaded field in FIG. 2. The scalar <img id="customcharacter00020" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is multiplied by the weighting factor <u style="single">k</u><sub><u style="single">2</u></sub><u style="single"></u>, so that a vector is obtained from this, whose first component has the dimension of an angular velocity, and whose second component equals zero. This vector is also sent to the adder 105. The vector resulting from the sum of the time derivative &xdot; of the phase vector <u style="single">x</u> formed according to Equation F 2.7 and of the vector obtained from the multiplication with <u style="single">k</u><sub><u style="single">2</u></sub><u style="single"></u> is integrated in the integrator 106 into the phase vector <u style="single">x</u>. One of the components β and {dot over (Ψ)} is eliminated from the phase vector by scalar multiplication with {dot over (Ψ)} and is further processed. While the estimated <img id="customcharacter00021" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is sent to the YTC control law unit 16, on the one hand, and to the direct process 103, on the other hand, the calculated {dot over (Ψ)} is used within the combined process only as a state variable within the observer and for determining the error of estimation. The difference between the yaw rate {dot over (Ψ)} determined from the observer vehicle model and the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>is formed for this purpose in the adder 107. This difference is multiplied by a vector <u style="single">h</u>, whose first component is dimensionless and sets the value of the correction steps for the velocity of the sideslip angle {dot over (β)}, and whose second component has the dimension s<sup>−1 </sup>and determines the value of the control steps during the correction of the yaw rate {dot over (Ψ)}.
The sideslip angle <img id="customcharacter00022" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is also returned as a correcting variable; specifically, it is fed back into the direct procedure of the kinematic determination of {dot over (β)} according to FIG. 7, so that the last term of Equation F 2.6 can also be assigned a value in the subsequent control step.
A substantially more accurate determination of the sideslip angle <img id="customcharacter00023" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is possible due to the mutual correction of the two calculation procedures, i.e., the calculation on the basis of a vehicle model and the calculation on the basis of kinematic considerations, so that this sideslip angle can also be sent as a controlled variable to the YTC control law unit 16.
2.3 Vehicle Reference Models
The vehicle reference model will be explained below on the basis of FIGS. 9 through 15.
FIG. 9 shows an even more simplified version of the control circuit according to FIG. 1 and FIG. 2 for controlling the driving stability of a vehicle. The controllers 7 through 9 in FIG. 1, the associated priority circuit 3 and the motor management 6 are omitted, and the distribution logic unit 2 is shown combined with the pressure control unit 5. An additional yaw torque M<sub>G </sub>around the vertical axis is calculated and set within the control circuit, so that the curve path desired by the driver is maintained. The additional yaw torque M<sub>G </sub>is generated by specific braking processes on the individual wheels, and the course of the braking processes and the selection of the wheels to be braked are set by the distribution logic unit 2. The desired direction of travel is set by the driver by selecting a corresponding angular position of the steering wheel. The steering wheel is coupled with the steered wheels in a fixed transmission ratio (steering ratio). A defined steering angle δ of the wheels is thus set.
2.3.1. Dynamic SingleTrack Model
A socalled vehicle reference model 12 (FIG. 2)=302 (FIG. 9), which is supplied with input data (velocity v, represented by v<sub>Ref</sub>, steering angle δ), is provided in the YTC controller 10. The size of the change in the yaw angle (yaw rate {dot over (Ψ)}<sub>Desired</sub>) per time unit is calculated in the vehicle reference model 302 on the basis of the input data. The desired value of the yaw rate {dot over (Ψ)}<sub>Desired </sub>is compared with the measured actual value of the yaw rate {dot over (Ψ)}<sub>Meas </sub>in a downstream comparison unit 303. The comparison unit 303 sends as an output value an output variable Δ{dot over (Ψ)}, which corresponds to the difference between {dot over (Ψ)}<sub>Desired </sub>and {dot over (Ψ)}<sub>Meas</sub>. The difference value thus determined is sent to a control law unit 16 for controlling the yaw torque. On the basis of Δ{dot over (Ψ)}, the control law unit calculates an additional yaw torque M<sub>G</sub>, which is sent to the distribution logic unit 2. Based on the additional yaw torque M<sub>G </sub>and possibly the driver's desire to build up pressure in the brakes, p<sub>Brake</sub>, the distribution logic unit 2 sets output variables. These may be brake pressure values or valve switching times.
Optimal mode of operation of the vehicle reference model 302 is also important in the range of low velocities. To ensure this, the vehicle reference model 302 may also be provided with a steadystate circular travel model 306, in addition to the abovedescribed linear dynamic singletrack model 311.
For the steadystate circular travel:<maths id="MATHUS00013" num="13"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mi>desired</mi></msub><mo>=</mo><mrow><mrow><mi>δ</mi><mo>*</mo><mfrac><mi>v</mi><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac><mo>*</mo><mfrac><mi>l</mi><mrow><mn>1</mn><mo>+</mo><mfrac><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>v</mi><mi>ch</mi><mn>2</mn></msubsup></mfrac></mrow></mfrac></mrow><mo>+</mo><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mi>korr</mi></msub></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.12</mn></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mi>β</mi><mo>=</mo><mrow><mrow><mi>δ</mi><mo>*</mo><mfrac><mrow><mn>1</mn><mo></mo><mrow><mfrac><msub><mi>l</mi><mi>v</mi></msub><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac><mo>*</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi>m</mi><mo>*</mo><msup><mi>v</mi><mn>2</mn></msup></mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo>*</mo><mrow><mo>(</mo><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mrow><mi>h</mi><mo>)</mo></mrow></msub></mrow></mrow></mrow></mfrac></mrow><mo>)</mo></mrow></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>v</mi><mi>ch</mi><mn>2</mn></msubsup></mfrac></mrow></mfrac></mrow><mo>+</mo><msub><mi>β</mi><mi>korr</mi></msub></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mi>in</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>which</mi></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.13</mn></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>v</mi><mi>ch</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mi>v</mi></msub><mo>*</mo><msub><mi>c</mi><mi>h</mi></msub><mo>*</mo><msup><mrow><mo>(</mo><mrow><msub><mi>l</mi><mi>h</mi></msub><mo>+</mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>m</mi><mo>*</mo><mrow><mo>(</mo><mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo>*</mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><msub><mi>c</mi><mi>v</mi></msub><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow><mo>)</mo></mrow></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.14</mn></mrow></mtd></mtr></mtable></math></maths>
Here, v=front; h=rear; m=weight; l=distance between the axle and the center of gravity; {dot over (Ψ)}<sub>korr</sub>, β{dot over (Ψ)}korr=correction terms for, {dot over (Ψ)}, β respectively.
The system equations F 1.1 and F 1.2 are valid for the linear dynamic singletrack model.
The switching over between the calculation models 306 and 311 is performed automatically by a changeover switch (not shown in the drawing) in the vehicle reference model 302 as a function of the velocity of the vehicle. A hysteresis of a few km/h is provided for switchover processes from one model to the other. Below the switching threshold, the desired yaw rate {dot over (Ψ)}<sub>Desired </sub>is calculated according to the model of steadystate circular travel. If the velocity, increasing from a lower value, exceeds the threshold that applies to this direction, the calculation of the desired value of the yaw rate {dot over (Ψ)}<sub>Desired </sub>is performed by means of the dynamic singletrack model 311. The dynamic processes that are particularly important for control at higher velocities are thus incorporated in the model.
The desired values calculated by the circular travel model, such as {dot over (Ψ)}<sub>Desired </sub>and β, are used as the starting values for the singletrack model when switching over from the circular travel model 306 to the singletrack model 311. As a result, transient effects during switchover are avoided. Further calculation is performed by means of the singletrack model 311 until the velocity drops below the velocity threshold, which is lower for decreasing velocity. To minimize transient effects here as well, the correction factors {dot over (Ψ)}<sub>korr </sub>and β<sub>korr </sub>necessary for the circular travel model are calculated with the values for {dot over (Ψ)}<sub>Desired </sub>and β, which were calculated before in the singletrack model, as well as with the velocity v<sub>Ref </sub>and the steering angle δ as the input variables.
The correction values are as follows:<maths id="MATHUS00014" num="14"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mi>korr</mi></msub><mo>=</mo><mrow><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mi>desired</mi></msub><mo></mo><mrow><mi>δ</mi><mo>*</mo><mfrac><mi>v</mi><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac><mo>*</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mfrac><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>v</mi><mi>ch</mi><mn>2</mn></msubsup></mfrac></mrow></mfrac></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.15</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>β</mi><mi>korr</mi></msub><mo>=</mo><mrow><mi>β</mi><mo></mo><mrow><mi>δ</mi><mo>*</mo><mfrac><mrow><mn>1</mn><mo></mo><mrow><mfrac><msub><mi>l</mi><mi>v</mi></msub><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac><mo>*</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi>m</mi><mo>*</mo><msup><mi>v</mi><mn>2</mn></msup></mrow><mrow><msub><mi>c</mi><mi>h</mi></msub><mo>*</mo><mrow><mo>(</mo><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow></mrow></mfrac></mrow><mo>)</mo></mrow></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>v</mi><mi>ch</mi><mn>2</mn></msubsup></mfrac></mrow></mfrac></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.16</mn></mrow></mtd></mtr></mtable></math></maths>
The effect of these correction factors decreases exponentially over time according to the equation:<FORM>korr(n+1)=korr(n)*λ F 2.17 </FORM>in which λ may assume values between 0 and less than 1. The calculation runs are counted with n and n+1.
Sudden changes are avoided as a result, because the two calculation methods yield different results in the steadystate case. Thus, the changeover between calculation models offers the possibility of determining the desired values for the control system at a rather high accuracy to velocities of v=0 km/h.
It was explained in connection with FIG. 9 that different models can be considered for use as vehicle calculation models. The steadystate circular travel may be a preferred model. The yaw rate {dot over (Ψ)}<sub>Desired </sub>can be calculated according to this model from the above formula. If such a vehicle calculation model is to be represented, it is possible to send the measured values λ and v<sub>Ref </sub>to a calculation circuit and to subsequently poll the desired value of the yaw rate {dot over (Ψ)}<sub>Desired </sub>as an output value.
2.3.3. Simplified Model
An extremely simple model for determining a desired yaw rate will be described below. It shall be an alternative to the abovedescribed combination model. It is characterized in that an acceptable result is obtained with a small amount of calculations.
The desired yaw rate {dot over (Ψ)}<sub>Desired </sub>is calculated according to this model as follows:<maths id="MATHUS00015" num="15"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mi>desired</mi></msub><mo>=</mo><mfrac><mrow><mi>δ</mi><mo>*</mo><mi>v</mi></mrow><mi>l</mi></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.18</mn></mrow></mtd></mtr></mtable></math></maths>
This equation is obtained from F 2.12, with equations F 2.14 and F 2.15 if the rigidities c<sub>r </sub>and c<sub>l </sub>are assumed to be very high.
This hypothesis is based on the following considerations.
In the vehicle reference model described above, the desired yaw rate {dot over (Ψ)}<sub>Desired </sub>is calculated either by means of a dynamic vehicle model (e.g., a singletrack model) or by a static model (called steadystate circular travel value) and is compared with the measured yaw rate {dot over (Ψ)}<sub>Meas</sub>. However, the preset value (and consequently also the control intervention) depend directly on the quality of the vehicle model in each of these hypotheses. Since these are linear equivalent models, the model markedly differs in some cases from the actual behavior of the vehicle.
If the real behavior of the vehicle additionally changes due to, e.g., load or wear of individual components, the model describes the vehicle only insufficiently. Consequently, adaptation of the model should be performed by means of a continuous parameter estimation, in connection with which the following problems arise:
An excitation must be present for the estimation, i.e., the driver should sufficiently excite the vehicle by means of a steering instruction in the linear range (<0.4 g). This hardly applies to normal driving.
Furthermore, it is not possible to directly estimate all parameters of the linear singletrack model. Thus, certain parameters should be preselected as fixed parameters.
Consequently, control on the basis of model hypotheses can always offer a satisfactory solution only regarding the model preset values. It may therefore be sufficient in many cases to proceed according to a simpler control principle.
One important goal of driving stability control is to coordinate the driving behavior such that the response of the vehicle to steering, braking and gas pedal inputs of the driver is always predictable and readily controllable. Consequently, understeering and oversteering operating conditions of the vehicle must be recognized and corrected to neutral behavior by a corresponding braking or engine management intervention.
The idea of simplified control principle is that a direct indicator of the understeering/oversteering behavior is used as a controlled variable. According to a definition of the steering behavior of a motor vehicle, the mean slip angles of the front axle and rear axle (α<sub>v</sub>, α<sub>H</sub>) are compared for this purpose. In the case of greater slip angles of the front axle, the vehicle thus exhibits an understeering behavior, and, in the opposite case, an oversteering behavior. According to the definition, neutral behavior is present if the slip angles front and rear are equal. Thus,<FORM>>0: understeering </FORM><FORM>α<sub>v</sub>−α<sub>h</sub>=0: neutral </FORM><FORM><0: oversteering applies. F 2.19 </FORM>
Based on the difference of the slip angles, it is consequently possible to directly determine the instantaneous driving condition of the vehicle. If the singletrack vehicle model (FIG. 10) is used as a hypothesis, the slip angles can be derived from this as a function of the steering angle δ, the sideslip angle β, the yaw rate {dot over (Ψ)} and the velocity of the vehicle v, as follows:<maths id="MATHUS00016" num="16"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>a</mi><mi>v</mi></msub><mo>=</mo><mrow><mi>δβ</mi><mo></mo><mfrac><mrow><mover><mi>Ψ</mi><mo>.</mo></mover><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mi>v</mi></mfrac></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.20</mn><mo></mo><mi>a</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mi>h</mi></msub><mo>=</mo><mrow><mi>β</mi><mo>+</mo><mfrac><mrow><mover><mi>Ψ</mi><mo>.</mo></mover><mo>*</mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mi>v</mi></mfrac></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.20</mn><mo></mo><mi>b</mi></mrow></mtd></mtr></mtable></math></maths>
Since the sideslip angle cannot be directly measured or calculated in a simple manner, an explicit calculation of the individual slip angles must be performed. However, if their difference is formed, it is possible to calculate this variable on the basis of the existing measured variables (steering angle, yaw rate), of the vehicle reference velocity v<sub>Ref </sub>known from the ABS controller and from the constant wheel base 1.<maths id="MATHUS00017" num="17"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>a</mi><mi>v</mi></msub><mo></mo><msub><mi>a</mi><mi>h</mi></msub></mrow><mo>=</mo><mrow><mi>δ</mi><mo></mo><mfrac><mrow><mi>l</mi><mo>*</mo><mover><mi>Ψ</mi><mo>.</mo></mover></mrow><mi>v</mi></mfrac></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.21</mn></mrow></mtd></mtr></mtable></math></maths>
Thus, a variable that can be used as an indicator of understeering/oversteering is available.
If the known relationship between the instantaneous curve radius R of the curve path of the center of gravity of the vehicle and the difference of the slip angles is also considered<maths id="MATHUS00018" num="18"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>δ</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>a</mi><mi>v</mi></msub><mo></mo><msub><mi>a</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.22</mn></mrow></mtd></mtr></mtable></math></maths>
It can be recognized that if a neutral state of the vehicle (F 2.19) is assumed<FORM>a<sub>v</sub>−a<sub>h</sub>=0 F 2.23 </FORM>the curve radius R can be determined only by the steering angle, namely,<maths id="MATHUS00019" num="19"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mi>δ</mi></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.24</mn></mrow></mtd></mtr></mtable></math></maths>
A control that directly uses the calculated slip angle difference as the controlled variable is therefore possible. The instruction for this control is to keep the value of this controlled variable as small as possible in order thus to achieve an approximately neutral behavior. It may be meaningful to assume this tolerance threshold to be asymmetric, so that the tolerance can be selected to be smaller in the direction of oversteering behavior.
The desired yaw rate {dot over (Ψ)}<sub>Desired </sub>can be calculated according to these considerations (F 2.18). This yaw rate {dot over (Ψ)}<sub>Desired </sub>is then compared with {dot over (Ψ)}<sub>Meas </sub>and is used as the basis of the control according to FIG. 1.
2.3.5 Desired Value Limitation in the YTC Controller
Controlling the driving behavior of the vehicle makes sense only as long as the adhesion of the wheels of the vehicle on the road surface permits the calculated additional torque to act on the vehicle.
It is undesirable, e.g., for the control to always force the vehicle to the curve path predetermined by the steering angle δ when the steering wheel was turned in excessively or too rapidly in relation to the existing velocity of the vehicle.
{dot over (Ψ)}<sub>Desired </sub>should therefore be prevented from always being selected as the preset value under all circumstances, according to the vehicle reference model selected, because if the reference model alone is followed, it may happen under unfortunate circumstances that if the steering wheel angle is accidentally set at an excessively high value, and the velocity is also high at the same time, the actual yaw rate {dot over (Ψ)} will be changed so much, due to the fact that {dot over (Ψ)}<sub>Desired </sub>is also too high in this case, that the vehicle will rotate around its own axis in the extreme case, while its center of gravity is moving in an essentially straight line at the same time. This condition is even much more unfavorable for the driver than the condition in which the vehicle is unable to obey the driver's desire due to the poor friction conditions and pushes out in a strongly understeering manner. This is true because the vehicle will at most only travel straight in this case, without also rotating around its own axis. To avoid these consequences, which are disadvantageous in special cases, calculation algorithms, which make it possible to set the maximum yaw rate {dot over (Ψ)}<sub>Desired max </sub>valid for the velocity just measured by way of the coefficient of friction {circumflex over (μ)}, are additionally provided in the vehicle reference model. {circumflex over (μ)} is determined in the coefficientoffriction detection system 13. The calculation algorithms are based on the theory of steadystate circular travel, for which {dot over (Ψ)}=a<sub>trans</sub>/v
(F 2.18).
The maximum allowable lateral acceleration a<sub>latlim </sub>can be determined essentially as a function of the coefficient of friction, the velocity v, the longitudinal acceleration a<sub>long</sub>, and possibly other parameters. Thus, F 2.25<FORM>a<sub>q lim</sub>=f(μ, v, a<sub>long</sub>, . . . ) </FORM>
The maximum yaw rate can be calculated as:<maths id="MATHUS00020" num="20"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mover><mi>Ψ</mi><mo>.</mo></mover><mrow><mi>Desired</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>max</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>a</mi><mrow><mi>q</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>lim</mi></mrow></msub><mi>v</mi></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>2.26</mn></mrow></mtd></mtr></mtable></math></maths>
It is therefore possible to set a limit value for the yaw rate, which does not take the driver's wish directly into account any longer, but it contributes to preventing the vehicle from additionally rotating around its vertical axis when it swings out.
Details of the suitable determination of μ will be described under 2.1.
Provisions can also be made to permit a control intervention only under certain prevailing conditions. One possibility for this may be, e.g., for the activation logic unit 11 in FIG. 2 to not transmit any current M<sub>G </sub>to the distribution logic unit 2 when an excessively large sideslip angle {circumflex over (β)} is determined, which can happen depending on the just occurring velocity.
2.4 Control Law Unit of the YTC Controller
The program structure of the control law unit 16 of the yaw torque controller 10 will be described below. From four input variables, the program calculates the additional yaw torque M<sub>G </sub>around the vertical axis of the vehicle that is necessary to obtain a stable vehicle behavior especially during travel in a curve. The yaw torque M<sub>G </sub>calculated is the basis for the calculations of the pressures to be applied to the wheel brakes.
The following input variables are available for the control law unit (see FIG. 17):
At input 500: Δ{dot over (Ψ)}
at input 501: Δ{umlaut over (Ψ)}
at input 502: <img id="customcharacter00024" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>
at input 503: {circumflex over (β)}
If the slip angle difference is used as a basis, Δ{dot over (Ψ)} is present at the input 500 and Δ{umlaut over (Ψ)} is present at the input 501.
Input 503 is facultative. It is available especially when a socalled observer vehicle model 84 is provided in the overall calculation system.
The value at input 500 is obtained as the difference between the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>and the desired yaw rate {dot over (Ψ)}<sub>Desired </sub>calculated by means of a vehicle reference model 12.
The value at input 501 is obtained either as a change in the variable at input 500 over time from one calculation loop to the next, divided by the loop time T<sub>0</sub>, or as a difference between the time derivative of the measured yaw rate and the time derivative of the calculated desired yaw rate.
A calculation loop is defined as a calculation run through the DSC driving stability controller according to FIG. 1. Due to its structure, such a loop requires a certain amount of real time, the loop time T<sub>0</sub>. This must be kept sufficiently short for an effective control.
The values at the inputs 500 and 501, namely, Δ{dot over (Ψ)} and Δ{umlaut over (Ψ)}, are first sent to a respective lowpass filter 510 and 511.
The two lowpass filters are, in principle, of the same design, and have the structure shown in FIG. 18.
The input variable 520 of the lowpass filter according to FIG. 18 is designated by u, and the output variable 521 is designated by y. The output variable 521 is sent to a register 522 and is available as a previous value y(k−1) at the time of the next calculation. The output value 521 for the calculation loop can then be calculated according to the formula<FORM>y(k)=λ*y(k−1)+(1−λ)*u*k<sub>p</sub> F 2.27 </FORM>in which λ may assume values between 0 and 1. λ describes the quality of the lowpass filter. The recursion function is eliminated at the limit value λ=0: The previous values y(k−1) are of no significance for the calculation of the new output value 521. The more closely λ approaches the value of 1, the stronger will be the effect of the previous values, so that the current input value 520 becomes established as an output value 521 only slowly.
k<sub>p </sub>is a linear rating factor.
The lowpass filtration just described is performed for both input values 500 and 501, and it leads to filtered values 515, 516.
An identical lowpass filtration 512 is performed for the input variable 502, namely, <img id="customcharacter00025" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>. The filtered value 517 is sent, just as the unfiltered value 503, to nonlinear filters 523, 524. These filters 523, 524 have the task of setting the output value to 0 for low input values and of transmitting an input value reduced by the limit value for input values that are above a certain limit value. The limitation is performed in the negative and positive ranges alike. The limit values {dot over (β)}th and βth may be fixed values implemented in the program, but they may also be variables that depend on other parameters, e.g., the coefficient of friction between the tires and the road surface. The limit values are calculated separately as a linear function of the coefficient of friction in this case.
All four variables, namely, 515, 516, 517 and 518, are weighted with a linear factor each in a next step 530, 531, 532 and 533, respectively.
These factors are implemented as fixed values in the calculation system. They can be calculated, in terms of their order of magnitude, from corresponding vehicle models, but they need, in general, a fine adjustment by driving tests. A corresponding set of linear factors is thus set for each vehicle or for each model of vehicle. The input variables 500, 501, 502, 503 thus weighted are added up, and (addition member 540) the additional yaw torque M<sub>G </sub>is obtained, which is used as the basis for the further calculation process of the program.
However, it was found in practice that modifications of the calculated yaw torque are still necessary.
Two approaches can be made to this end:
1. The input variables, especially Δ{dot over (Ψ)}, are modified.
2. The calculated yaw torque M<sub>G </sub>is subjected to filtration.
Attempts are made with both approaches to perform the control not only in consideration of the yaw rate, but also in consideration of the sideslip angle.
2.4.1 Modification of the Input Variables
As explained, a desired value is calculated for the yaw rate by means of a vehicle reference model. Since the vehicle reference model cannot completely agree with the actual conditions, it is usually necessary to correct the result of the model calculation once again. The values which are provided by a yaw rate sensor, as well as a steering angle sensor, are essentially evaluated in the reference model. Correction of the calculated desired yaw rate can be performed by additionally taking into account the values provided by a lateral acceleration sensor.
The evaluation may be performed in various manners. One way is proposed below, according to which the measured lateral acceleration is first converted into a velocity of the sideslip angle <img id="customcharacter00026" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>. A correction of the desired value for the yaw rate is performed with this value.
The calculation of <img id="customcharacter00027" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is performed, e.g., by the kinematic determination of {dot over (β)} 14, 15 (FIG. 2).
The procedure is carried out according to the scheme shown in FIG. 19. The estimated value of the velocity of the sideslip angle <img id="customcharacter00028" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is compared with a first threshold value th<sub>1 </sub>(block 400), if desired, after a lowpass filtration. The meaning of this comparison will appear only after a correction of the desired value of yaw rate {dot over (Ψ)}<sub>Desired</sub>, and it is therefore explained in greater detail below.
If <img id="customcharacter00029" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>>th<sub>1</sub>, the value of <img id="customcharacter00030" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is compared with a second threshold value th<sub>2 </sub>(block 401), and the second threshold value is higher than the first threshold value th<sub>1</sub>. If this threshold value is also exceeded, integration 402 of the velocity of the sideslip angle <img id="customcharacter00031" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> over time is first performed. To do so, the velocity of the sideslip angle <img id="customcharacter00032" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is multiplied by the loop time T<sub>0 </sub>and added to the previous integration result Intg<sub>i1</sub>. The integration steps are counted with n, so that the number n is increased by 1 after the integration (step 403). The integration time is thus represented by the number n of integration steps performed. The integration result Intg<sub>n </sub>(<img id="customcharacter00033" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/>) is compared with a threshold value βs (block 404). The amount of the threshold value represents a maximum allowable deviation from a sideslip angle that is theoretically to be maintained. The threshold value βs is in the order of magnitude of approximately 5°.
If this threshold value is exceeded, the desired yaw rate {dot over (Ψ)}<sub>Desired </sub>is newly evaluated by an additive constant S (step 405), which depends on the instantaneous velocity of the sideslip angle <img id="customcharacter00034" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> and the number n of integration steps. This means that the desired yaw rate is further reduced with each new loop in which the threshold value βs is exceeded. The additive constant S is either added or subtracted, depending on the sign of {dot over (Ψ)}<sub>Desired</sub>, so that the value of the desired yaw rate is reduced at any rate. If Intg<sub>n </sub>does not reach the threshold value βs, {dot over (Ψ)} is not limited (step 407).
The estimated velocity of the sideslip angle is checked again in a repeated loop to determine whether its value is lower than the threshold th<sub>1</sub>. If so, this is interpreted as meaning that the vehicle has again stabilized. The consequence of this is that n in step 406 is again set at 0 and that the further calculation in step 407 is based on a desired yaw rate that is not corrected, i.e., it is identical to the value obtained as the result of the vehicle reference model. In addition, the start value Intg<sub>n1 </sub>of the integration is set to equal zero.
If the value of a velocity of the sideslip angle exceeds th<sub>1</sub>, but not th<sub>2</sub>, the old value Intg<sub>n </sub>remains unchanged, i.e., the integration is omitted for one loop. The previous limitation is preserved. Should the threshold value th<sub>2 </sub>be exceeded again, the integration is continued.
2.4.2 Correction of M<sub>G </sub>
Another possibility is to manipulate the yaw torque M<sub>G</sub>, which is calculated by the control law unit 16. To do so, the difference between the previous value M<sub>1</sub>(k−1) and the current value M<sub>1</sub>(k) is formed. The subscript 1 indicates that these values are the direct result of the yaw torque controller, i.e., they were not yet calculated on the basis of the next correction. This difference is related to the loop time T<sub>0 </sub>and yields ΔM<sub>1</sub>. A correction gradient, which is obtained from <img id="customcharacter00035" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> multiplied by a correction factor, is added to this gradient ΔM<sub>1</sub>. The gradient thus corrected is multiplied by the loop time T<sub>0 </sub>and is added to the yaw torque M<sub>1</sub>(k−1) of the preceding calculation. This leads to the current moment M<sub>G</sub>(k), which is used as the basis for the further calculation.
This calculation is performed by a logic unit as is shown in FIG. 20. The calculated moments, which are obtained from the “control law unit 16” subprogram, are sent into a shift register 420. The current value M<sub>1</sub>(k) stands in the first place 421 of the shift register 420; the previous value M<sub>1</sub>(k−1) stands in the second place 422 of the shift register 420. As soon as a new value M<sub>1 </sub>is available, the value is shifted from register 421 into register 422, and the value in register 421 is replaced with the new value. The values in the registers 421 and 422 are sent to a calculation logic unit 430, which calculates a ΔM according to the following formula:<FORM>ΔM=M<sub>1</sub>(k)−M<sub>1</sub>(k−1)+a*{dot over (β)}*T<sub>0</sub> F 2.28 </FORM>
In addition, the estimated velocity of the sideslip angle <img id="customcharacter00036" he="3.56mm" wi="1.44mm" file="US20060273657A120061207P00002.TIF" alt="custom character" imgcontent="character" imgformat="tif"/> is sent to the calculation logic unit 430 from the kinematic determination of {dot over (β)}. Furthermore, a value for a correction factor a, with which the velocity of the sideslip angle is converted into a change in moment, is set in a memory. The new moment M(k) is calculated according to the formula<FORM>M(k)=M(k−1)+ΔM F 2.29 </FORM>
The current value of the corrected moment is stored in register 431, and the value from the previous calculation is stored in register 432. The value in register 431 is used as the basis for the further calculation.
2.5 SESP Yaw Rate Control Law Unit
The structure of the SESP system is schematically shown in FIG. 35.
The program of the SESP yaw rate control law unit 23 compares the measured yaw rate with a specific SESP reference yaw rate. When a deviation exists, the controller demands an additional yaw torque MG<sub>SESP </sub>for the correction of the deviation.
2.5.1 SESP Reference Signal Formation Unit
The SESP reference signal formation unit 24 calculates the SESP reference yaw rate which is the yaw rate of the vehicle intended by the driver. The result is stored in sesp_psip_ref.
The SESP yaw rate controller is more sensitive than the YTC yaw rate controller 10. For this reason, it requires a special mechanism to avoid erroneous interventions due to sensor and model errors. The SESP reference yaw rate contains an offset SESP_DPSIP_STORED compensating these errors for a limited time.
If a turning situation is not suspected, the offset SESP_DPSIP_STORED represents the deviation between the measured or estimated yaw rate {dot over (ψ)}Meas and the reference yaw rate {dot over (ψ)}<sub>Desired </sub>calculated in the vehicle model 12.
However, as soon as there is a suspected turning situation, the offset is kept constant. A turning situation is suspected when the following condition is fulfilled:
 xx. SESP state 31 to 34 prevails (unlike 35)
 or
 xxi. a possible previous cause i to xviiii for the turning maneuver has been detected.
Now the SESP reference yaw rate {dot over (ψ)}<sub>SESP(ref) </sub>is calculated. The SESP reference yaw rate represents the driver's request. For the calculation, the offset is added to the reference yaw rate of the YTC (=AYC) reference yaw rate:<FORM>{dot over (ψ)}YTC<sub>ref</sub>+Offset={dot over (ψ)}<sub>SESP(ref) </sub></FORM>
The offset is the difference between {dot over (ψ)}<sub>Desired</sub>−{dot over (ψ)}<sub>Meas </sub>in the event of e.g. a load change.
Compared to the YTC (=AYC) reference yaw rate, the SESP reference yaw rate comprises an offset which is just so great that the SESP yaw rate deviation at the beginning of the load change or the initial phase of braking is zero.
This offset is used to compensate possible uncompensated offsets of the yaw rate sensor, on the one hand. This is necessary to avoid control errors of the sensitive SESP yaw rate controller. On the other hand, the offset compensation only allows adjusting the deviations from the current vehicle performance.
FIG. 31 shows the formation of the SESP reference yaw rate for the case of a positive yaw rate sensor offset.
SESP reference signal formation unit 24 calculates the yaw rate deviation Δ{dot over (ψ)}<sub>SESP </sub>being used as an input for the SESP yaw rate controller. Further, FIG. 2 shows that this input Δ{dot over (ψ)}<sub>SESP </sub>is also made available to the SESP activation logic unit 25.
The SESP yaw rate deviation is calculated as the difference between the measured yaw rate {dot over (ψ)}<sub>Meas </sub>and the SESP reference yaw rate.
Further refinement is envisaged when producing the offset correction for defined applications.
A state machine is provided for this purpose which makes a distinction between the following situations:
Sesp_refcomp_idle: SESP reference variables are not required
Sesp_refcomp_straight: SESP reference variables can be determined, they are valid for straight travel
Sesp_refcomp_curve: SESP reference variables can be determined, they are valid for cornering maneuvers
Sesp_refcomp_uncertain: SESP reference variables are required, but cannot be exactly determined due to the dynamic driving situation.
The (additional) state machine is so designed that a changeover to a more appropriate state may be performed even during a current reference offset correction. If required, a new determination of the offset is triggered in this case.
The SESP reference variables are now produced depending on the current situation as will be described in the following:
Sesp_refcomp_idle: is not produced (as a substitute SESPreference=actual quantity)
Sesp_refcomp_straight: SESP reference=AYC reference+offset, with the offset being determined as described and e.g. illustrated in FIG. 31
Sesp_refcomp_curve: if understeering prevails upon entry into this situation: SESP reference=AYCreference; otherwise: SESP reference=AYC reference+offset.
The risk of SESP control errors is reduced due to this distinction of cases. This can be illustrated in FIG. 34 by an extension, and
Sesp_refcomp_uncertain: SESP reference=AYC reference (=−best guess”).
applies.
Due to the inaccurately defined SESP reference, the activation thresholds of the SESP controllers are slightly raised in this case.
2.5.2 SESP Activation Logic Unit
The SESP yaw rate control law unit 23 includes time delays in order to avoid very short activations or interruptions of the control. The delay for the controller input is stored as sesp_dpsip_in_delay. The delay for the controller output is stored as sesp_dpsip_out_delay.
The short initial delay sesp_dpsip_out_delay=Sesp_dpsip_in_delay_tab[0] is chosen if there is a tendency to a greater oversteering; that means if all of the following conditions are satisfied:
 a. Δ{dot over (ψ)}<sub>SESP</sub>>threshold value k20; the threshold value k20 is situationresponsive, e.g. responsive to the vehicle speed (states 31, 32 more sensitive than 33, 34), of the transversal inclination of the roadway and the like,
 b. {dot over (ψ)}<sub>SES Pref</sub>>{dot over (ψ)}<sub>Mess</sub>, wherein the signs of {dot over (ψ)}<sub>SES Pref</sub> and {dot over (ψ)}<sub>Meas</sub> must be identical; the condition indicates that there is a driving situation with an oversteering behavior; is only applicable to the states 31 or 32; the entry time is responsive to the situation, e.g. to {umlaut over (ψ)}<sub>Meas</sub>.
The oversteering situation is only used as a precondition for the controller activation. Consequently, the first activity of the SESP yaw rate controller will always be an oversteering intervention. But the following activities can also be understeering interventions in order to prevent understeering reactions of the vehicle. The oversteering flag for the standard GTC(AYC) reference yaw rate can also be used as a precondition for the controller activation to ensure that no SESP oversteering interventions can be started during an actual understeering situation. This option is shown in FIG. 34 as variant (2).
 c. The state machine 22.1 adopts a state unlike 35
It applies for the conditions a. to c. that they must be fulfilled for a defined duration.
The deactivation depends on the following conditions:
 d. {dot over (ψ)}<sub>SES Pref</sub><threshold value k21; threshold value k21 is situationresponsive, it is lower than the entry threshold k20;
 or
 e. state machine 22.1 adopts state 35; as this occurs, the state ‘SESP inactive’ must prevail for a defined period before the deactivation takes place. The exit delay is responsive to the situation, it is zero if state 35 prevails, otherwise it depends on {umlaut over (ψ)}<sub>Meas. </sub>2.5.3 SESP Control Law Unit
The SESP control law unit 23 calculates the additional yaw torque which must be demanded from the SESP yaw rate controller in order to correct the SESP yaw rate deviation. The result is stored in the value sesp_dpsip_req_yaw_tq.
((Initially, the SESP yaw rate deviation Δ{dot over (ψ)}<sub>SESP</sub>))* is calculated by the exit threshold k21 and limited downwards to the minimum value zero. The quantity calculated there is the input value for a Pcontroller.
The amplification factor depends on the estimated coefficient of friction of the road. As the coefficient of friction of the coefficientoffriction estimation unit 13 is not active for the SESPrelevant driving situations, a substitute value is used which is determined from the measured a<sub>lat </sub>and depending on vref from the vehicle longitudinal acceleration. This coefficient of friction represents the currently utilized coefficient of friction. The amplification factor of the Pcontroller depends on this coefficient of friction. The additional yaw torque MG<sub>SESP </sub>prevails at the output of the control law unit 23.
2.5.4 SESP Distribution Logic Unit
The SESP distribution logic unit distributes the additional yaw torque request MG<sub>SESP </sub>from SESP to the wheel brake pressures. There are different variants depending on the brake pressure of the driver. The present yaw torque distribution state of the SESP is stored in SESP_TQDIS_STATE.
SESP influences the brake pressure of at least one wheel. The following variants are defined for this intervention:
 Pressure reduction mode: The wheel brake pressure at the inside rear wheel in a turn and, if required, also additionally at the inside front wheel in a turn is reduced up to an envisaged pressure bottom limit ayc_pdec_sec_pressure_limit. If the potential by pressure reduction is not sufficient to realize MG<sub>SESP</sub>, pressure increase is not performed at this time for comfort reasons. Pressure increase is, however, also possible according to a strategy in conformity with comfort aspects.
 Pressure increase mode: The wheel brake pressure at the outside rear wheel in a turn is increased until the upper limit of pressure is reached which is defined by the maximum possible force in a longitudinal direction. The upper limit of pressure is ensured by the wheel slip controller 26 with the brake slip monitoring function in a combination with suppression of the EBD function. The pressure is increased for comfort reasons by way of the function ETR (=changeover valve/separating valve control).
 Distribution of the additional yaw torque from SESP to the wheels is arranged for defined applications.
 In the case of SESP pressure increase, it is now possible to build up pressure also on the outside front wheel in a turn.
 How much percent of the pressure of the rear wheel is built up at the front wheel in addition can be adjusted by way of a characteristic curve depending on the estimated minimum coefficient of friction SESP_MY_MIN.
 The following dependency has proven suitable for a brake system with hydraulic brake force boosting: 50% on ice, 100% on medium coefficients of friction, 0% on high coefficients of friction.
 A value of 0% on high coefficients of friction is appropriate because the driver considers pressure increase at the front axle as uncomfortable in this case. On the other hand, a sole pressure increase at the rear axle is almost sufficient herein for stabilization because more longitudinal force can be built up on high coefficients of friction.
 The possibility of pressure increase at the front axle is illustrated in FIG. 32a.
 When SESP pressure increase is terminated, the pressure at the rear wheel is no longer decreased abruptly but slowly. This fact increases driving comfort.
 Slow pressure reduction is also employed if SESP must be discontinued due to AYC control beginning.
 However, pressure is decreased abruptly as soon as the sign of the AYC additional yaw torque changes compared to the previous SESP intervention. The option of an additional pressure reduction on the front axle, which has been mentioned before regarding the SESP pressure reduction case (page 17) and in FIG. 32, is preferred in particular for vehicles equipped with a brake system with vacuum brake booster and additional hydraulic brake force boosting, e.g. by way of the hydraulic return pump of the brake system.
The decision whether pressure reduction or pressure increase shall be executed depends on the following conditions:
Pressure Reduction
 f. YTC controller 10 has detected braking by driver P<sub>TMC </sub>and
 g. P<sub>TMC</sub>>threshold value k30
In the absence of these conditions, pressure increase takes place.
If the pressure increase mode is active, MG<sub>SESP </sub>is reduced depending on vref. It is thereby taken into account that the additional yaw torque MG<sub>SESP </sub>can be realized more effectively during pressure increase than during pressure reduction.
The decision whether the brake intervenes on the vehicle side which is on the left or the right in a turn is taken according to the same rules as in the YTC yaw torque distribution logic unit 2. Of course, it is also possible to include the SESP distribution logic unit 21 in the process of reaching a decision.
Individual wheel pressure demands of SESP can be taken from the distribution logic unit 21 in consideration of the boundaries of the wheel slip controller 26. In addition, it is indicated to the other controllers 7, 8, 9 by way of a flag that the SESP control is active.
This is favorable because the YTC controller 10 adds the offset to the input signal Δ{dot over (ψ)} of the activation logic unit 11 in order to make it possible for the SESP to stabilize the vehicle to a certain degree also above its control thresholds. This possibility allows SESP control in a transition range in borderline situations, with almost unstable driving situations, however, only if SESP has become active before. The offset is variable and can be reduced until zero.
3. Distribution Logic Unit of the YTC Controller
3.1 Additional Yaw Torque by the Application of Brake Forces
To achieve stable travel of the vehicle even in a curve, it is first necessary to determine the steering angle. The steering angle represents the curved path of the vehicle desired by the driver. In the case of stable, steadystate travel in a curve, the vehicle shall travel through the curve at an approximately constant sideslip angle and constant yaw rate. Deviations from this sideslip angle or from this yaw rate must be compensated by the driver by steering in the opposite direction. However, this is not always possible when the driver travels through the curve at the limit velocity for the curve. It is necessary in such situations to specifically brake the vehicle and to apply additional moments around the vertical axis to the vehicle, which are to bring about an adjustment of the actual yaw rate to the desired yaw rate. Control algorithms which describe these relationships have been described before, so that they need not be explained in greater detail herein.
However, there remains the problem that an additional yaw torque M<sub>G </sub>calculated by the control algorithm must be put into practice in an appropriate manner by specifically applying brake forces.
In the chase of hydraulic brakes, the task is therefore practically to set a brake pressure for every individual wheel brake. The moment to be obtained around the vertical axis shall be obtained with the lowest possible pressures in the individual brakes. It is therefore proposed that a coefficient be determined for each wheel and that the brake pressures be calculated from the vehicle yaw torque to be generated and the actual weighted coefficient.
As explained above, it is favorable, especially in vehicle brake systems operating hydraulically, to determine the coefficients such that the brake pressure for the individual wheel brakes can be directly determined. The weighting of the coefficients is performed by dividing every individual coefficient by the sum of the squares of all coefficients.
Each coefficient determines the relationship between the wheel brake pressure and the individual wheel brake forces thus generated as a percentage of the yaw torque of the vehicle. Parameters which change during the travel of a vehicle are included as variables in the determination of the individual coefficients. They are, in particular,
 the steering angle δ
 the coefficient of friction μ between the tire and the road surface,
 the vehicle mass m, and
 the axle load distribution N<sub>z</sub>.
Variables which are included in the calculation of the coefficients and are vehiclespecific or brakespecific are, e.g., the following, for a disc brake system:
 the area A of the brake pistons,
 the number n of pistons per wheel brake,
 the coefficient of friction μ<sub>R </sub>between the disc and the brake lining,
 the ratio s of the effective friction radius to the dynamic tire radius, and
 the efficiency η of the brake.
The method of calculation proposed has the advantage that the corresponding brake pressures can be calculated very rapidly from a predetermined additional yaw torque. Should the abovedescribed parameters change during travel, this is taken into account by way of a change in the coefficients in the calculation of the brake pressure.
While some influencing variables are used linearly in the calculation of the coefficients, especially the dependency of the coefficients on the steering angle δ is nonlinear.
However, it was found that a linearized estimation of the dependency between the individual coefficients and the steering angle yields sufficiently good results.
FIG. 21 schematically shows a vehicle during straight travel with four wheels 601, 602, 603, 604. A wheel brake 605, 606, 607, 608 is associated with each of the wheels. These can be actuated independently from one another, and brake forces are generated by the wheel braking moments exerted by the wheel brakes on the contact surfaces of the tires on the road surface. For example, a braking force F, which in turn generates a moment M (positive in the example) around the vertical axis, is generated on wheel 601 when the wheel brake 605 is actuated.
Such moments around the vertical axis of the vehicle can be used specifically to keep a vehicle stable on a path desired by the driver.
Furthermore, sensors are provided in the vehicle. They include wheel sensors, which detect the angular velocity of the wheels 601, 602, 603, 604. In addition, the steering wheel angle is detected with a steering sensor 612. In addition, a sensor 613 for the yaw rate is provided.
A yaw torque, which, when applied, is able to make the yaw rate of the vehicle as well as its sideslip angle agree with the driver's desire, can be calculated with these sensors, which detect the driver's desire, on the one hand, and the behavior of the vehicle, on the other hand. The wheel brakes 605, 606, 607, 608 are actuated independently from one another for this purpose, using a control device which is part of a complex program for controlling the driving stability.
The general situation is shown in FIG. 22. A program module, which calculates the yaw torque M<sub>G</sub>, is designated by reference numeral 16. FIG. 22 shows a control device, which calculates the pressures p<sub>xx </sub>that are to be introduced into the individual wheel brakes 605, 606, 607, 608. The pressure values 622, 623, 624, 625 determined can be subjected to further evaluation and can be converted into corresponding control signals for the wheel brakes 605, 606, 607, 608.
The control device itself consists of two parts, namely, a first part 630, in which coefficients c<sub>xx </sub>for the individual wheels are calculated. The coefficients c<sub>xx </sub>establish a linear relationship between the pressure in the wheel brake and the proportionate yaw torque, which is brought about by the brake force on the corresponding wheel. The individual pressure values p<sub>xx </sub>622, 623, 624, 625 are calculated in the second part 631 by weighting the individual coefficients and taking into account the yaw torque M<sub>G </sub>to be applied.
The pressure values as well as the coefficients are designated with subscripts:
V: front h: rear
l: left r: right
X: either v/l or h/r
The first calculation part 630 takes into account the steering angle, which is made available to the calculation process by way of an evaluation 632 of the steering sensor 612. To calculate the coefficient, the coefficient of friction p, which is derived from the wheel rotation behavior in an evaluation unit 633 (cf. Section 2.1), is taken into account. The wheel rotation behavior is in turn determined by a signal of the wheel sensors at the individual wheels. The mass of the vehicle as well as the load distribution N<sub>Z</sub>, which are determined in an evaluation unit 634, in which the behavior of the vehicle is analyzed in different situations, are included as well. The first program part 630 has access to a memory 635, which contains the abovementioned vehiclespecific and wheel brakespecific values.
A coefficient c<sub>xx </sub>is calculated from the abovementioned values for each wheel; the values 640, 641, 642, 643 may be calculated simultaneously or consecutively. The calculation is performed according to a function implemented in the program. The known relationships between the brake pressure and the brake force are taken into account in this function. The relationship is usually linear. Only the steering angle δ must be taken into account separately. How the steering angle can be taken into account in a suitable manner will be described below.
The pressure values for the individual wheel brakes are determined in the second calculation step 631 either simultaneously or consecutively from the individual coefficients 640, 641, 642, 643 according to the following formula:<maths id="MATHUS00021" num="21"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>p</mi><mrow><mi>x</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1</mn></mrow></msub><mo>=</mo><mrow><mfrac><msub><mi>c</mi><mi>xl</mi></msub><mrow><msubsup><mi>c</mi><mrow><mi>v</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1</mn></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mrow><mi>v</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>r</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>hl</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>hr</mi><mn>2</mn></msubsup></mrow></mfrac><mo>*</mo><msub><mi>M</mi><mi>G</mi></msub></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>3.1</mn><mo></mo><mi>a</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>p</mi><mi>xr</mi></msub><mo>=</mo><mrow><mfrac><msub><mi>c</mi><mi>xr</mi></msub><mrow><msubsup><mi>c</mi><mrow><mi>v</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>1</mn></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>vr</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mrow><mi>h</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>l</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>hr</mi><mn>2</mn></msubsup></mrow></mfrac><mo>*</mo><msub><mi>M</mi><mi>G</mi></msub></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>3.1</mn><mo></mo><mi>b</mi></mrow></mtd></mtr></mtable></math></maths>
Calculating the individual pressures according to this formula offers the advantage that only relatively low pressures must be introduced into the wheel brakes to reach the calculated braking moment. Furthermore, the brake pressure control is able to respond very sensitively and rapidly to changes especially in the steering angle and in the coefficients of friction.
The steering angle δ is taken into account in the calculation of the coefficients as follows: FIG. 23 shows for this a schematic representation of a vehicle, in which the front wheels 601 and 602 are shown turned in. The distance between the front wheels is designated by S, and the distance between the center of gravity 610 and the front axle is designated by l<sub>v</sub>.
The wheel planes 650, 651 form steering angles 652, 653 with the longitudinal axis of the vehicle. The steering angles δ 652, 653 are assumed to be equal for simplicity's sake. The effective lever arm h<sub>l </sub>or h<sub>r </sub>relative to the brake force F, which acts in the wheel plane 650, 651, is calculated from approximation considerations for small steering angles as follows.<maths id="MATHUS00022" num="22"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>h</mi><mi>r</mi></msub><mo>=</mo><mrow><mfrac><mi>s</mi><mn>2</mn></mfrac><mo>+</mo><mrow><mi>δ</mi><mo>*</mo><munder><mi>l</mi><mi>v</mi></munder></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>3.2</mn><mo></mo><mi>a</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>h</mi><mi>l</mi></msub><mo>=</mo><mrow><mfrac><mi>s</mi><mn>2</mn></mfrac><mo></mo><mrow><mi>δ</mi><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>3.2</mn><mo></mo><mi>b</mi></mrow></mtd></mtr></mtable></math></maths>
Since the “small steering angle” approximation is not always satisfied, it was found to be favorable to calculate with the following formula, if applicable.<maths id="MATHUS00023" num="23"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>h</mi><mi>r</mi></msub><mo>=</mo><mrow><mfrac><mi>s</mi><mn>2</mn></mfrac><mo>+</mo><mrow><mi>δ</mi><mo>*</mo><msqrt><mrow><mfrac><msup><mi>s</mi><mn>2</mn></msup><mn>4</mn></mfrac><mo>+</mo><msubsup><mi>l</mi><mi>v</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mrow><mn>3.3</mn><mo>.</mo><mi>a</mi></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>h</mi><mi>l</mi></msub><mo>=</mo><mrow><mfrac><mi>s</mi><mn>2</mn></mfrac><mo></mo><mrow><mi>δ</mi><mo>*</mo><msqrt><mrow><mfrac><msup><mi>s</mi><mn>2</mn></msup><mn>4</mn></mfrac><mo>+</mo><msubsup><mi>l</mi><mi>v</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>3.3</mn><mo></mo><mi>b</mi></mrow></mtd></mtr></mtable></math></maths>
Should the calculated lever arms become smaller than zero, they are set equal to zero.
The wheel coefficients c<sub>xx </sub>can be calculated as follows:<FORM>c<sub>xx</sub>=c<sub>hydxx</sub>*h<sub>l,r</sub>. F 3.4 </FORM>and all parameters except for the steering angle δ are taken into account in c<sub>hydxx</sub>.
The coefficients can thus be represented as the product of two terms, in which one term determines the effective lever arm, and the other term is independent from the steering angle.
3.2 Additional Yaw Torque by Reducing Lateral Forces
One method of applying brake forces acting on one side is to actuate the wheel brakes such that the wheels will be braked with different intensity. One procedure that brings this about was described in the preceding section.
This procedure reaches a limit when driving stability control is to be performed during pedal braking, i.e., when a certain brake pressure has already been set in the wheel brakes because of braking by the driver. The abovedescribed procedure can be applied, in principle, to this case as well. Instead of absolute pressures, changes in the brake pressures already set are determined.
However, the following problems arise. If a very high pressure has already been introduced into a wheel brake, so that very high brake forces are reached, an increase in the brake pressure would not necessarily lead to an increase in the brake force, because the limit of adhesion between the tire and the road surface has been reached. The linear relationship between the brake pressure and the brake force, which was assumed in the abovementioned model, is no longer present in this case.
The limit of the brake force on one side of the vehicle, which is not to be exceeded, can be compensated in terms of a yaw torque control by reducing the braking force on the other side of the vehicle.
However, this has the disadvantage that the deceleration of the vehicle is also reduced with the reduction in the brake force. This is not always acceptable, because the vehicle is to be stopped over the shortest possible distance when a braking process is initiated by the driver. Therefore, an excessive reduction in the actual deceleration of the vehicle compared with the driver's desire cannot generally be accepted. The following approach is taken to solve this problem.
The wheel brakes of at least one wheel are actuated such that the longitudinal slip 2 of the wheel is set such that it is greater than the longitudinal slip at which the maximum frictional connection is reached. This procedure is based on the fact that the brake force transmitted, i.e., the longitudinal force on the tire, reaches its maximum at a longitudinal slip of approximately 20% (0%=freely rolling wheel; 100%=locked wheel), and the brake force that can be transmitted decreases only slightly at values above 20%, so that there is no appreciable loss during the deceleration of the vehicle at wheel slips between 20% and 100%.
However, if the lateral force that can be transmitted, i.e., the force that acts at right angles to the wheel plane, is also taken into account at the same time, a strong dependence on wheel slip is seen, which is manifested in that the lateral force that can be transmitted greatly decreases with increasing slip. In the slip range above 50%, the wheel exhibits a behavior similar to that of a locked wheel, i.e., hardly any lateral forces are applied.
Controlled skidding of the vehicle can be provoked by judiciously selecting the wheels on which a high longitudinal slip is set, and the change in the yaw angle brought about by the skidding shall correspond to the desired change. Since the longitudinal forces are essentially preserved in this procedure, but the lateral forces are markedly reduced, the yaw rate can be controlled without excessively reducing the deceleration of the vehicle.
The wheel that is driven, at least briefly, with an increased longitudinal slip is selected according to the following rules. Let us examine travel in a curve to the right, which is intended by the driver. Corresponding “mirrorimage” rules apply to traveling in a curve to the left. The case may occur in which the vehicle will not turn into the curve as sharply as expected. In other words, the vehicle is understeered. The rear wheel that is the inner wheel in the curve is operated with increased slip values in this case. However, if the vehicle turns too sharply into the curve—this case is called oversteering—the front wheel that is the other wheel in the curve is operated with high slip values.
In addition, the pressure can be prevented from decreasing on one front wheel. This is done according to the following rules. In a driving situation in which the vehicle exhibits understeering behavior, the brake pressure is prevented from decreasing on the front wheel that is the outer wheel in the curve. The pressure is prevented from decreasing on the front wheel that is the inner wheel in the curve in a situation in which the vehicle exhibits oversteering behavior.
The actual control of the brake pressure may be performed as follows. As was explained before, the brake pressure in the individual wheel brakes is determined individually as a function of the yaw torque to be reached and the weighted wheel coefficients.
A factor which is independent from the brake slip can be introduced in the calculation of the coefficients; this factor is adjusted such that the abovedescribed desired brake slip will become established. The reduction in pressure on a wheel can be limited by setting a lower threshold for the corresponding coefficient.
The procedure implemented in the control program of the brake system will be explained in greater detail below.
Based on weighted coefficients, the control program calculates the brake pressure that must be produced in every individual wheel brake. The calculation becomes more problematic when the vehicle is braked, especially when it is being decelerated while utilizing the limit of frictional connection between the tire and the road surface. It is quite possible in such cases that an antilock control will first begin before a superimposed driving stability control becomes necessary.
The basic considerations for an unbraked vehicle cannot be taken over in such cases, because, e.g., the corresponding brake force does not increase linearly upon the increase in pressure in a wheel brake, since the limit of frictional connection has been reached. An increase in the pressure in this wheel brake would not consequently produce any additional brake force and consequently any additional moment.
Even though the same effect of generating an additional yaw torque can be produced by reducing the wheel brake pressure of the other wheel of the axle, this would cause, on the whole, a reduction in the braking force, which in turn conflicts with the requirement that the vehicle is to be stopped over the shortest possible distance.
The behavior of vehicle wheels shown in FIG. 24 is therefore utilized. This diagram shows slip values λ between 0% and 100% on the X axis, where 0% indicates a freely rolling wheel and 100% a locked wheel. The Y axis shows the frictional force and lateral force values μ<sub>B </sub>and μ<sub>s</sub>, respectively, in the range of values from 0 to 1. The solid lines show the dependence of the coefficient of friction on slip for different slip angles α. It is seen, especially in the case of small slip angles, that the curve has a maximum in the slip range of λ=20%. The coefficient of friction slightly decreases toward 100%. The maximum coefficient of friction equals approximately 0.98 for a slip angle of 2°, while it is still 0.93 at λ=100%. However, an examination of the values of the lateral force shows an extreme reduction over the slip range, especially for great slip angles. The value of the lateral force for a slip value of 0% is 0.85 at a slip angle of 10°, to drop to 0.17 for slip values of almost 100%.
Thus, it can be determined from the curves in FIG. 24 that relatively strong brake forces, but weak lateral forces can be transmitted at slip values in the range of 40% to 80%.
This behavior of the wheel can be utilized to specifically reduce the lateral force of a given wheel of the vehicle. The wheel is selected according to the following scheme, which will be explained in greater detail on the basis of FIGS. 25a and 25b.
FIGS. 25a and b show a schematic representation of a vehicle in a right curve. Corresponding to the radius of the curve and the velocity of the vehicle, the vehicle must turn around its vertical axis, i.e., there must be a defined clockwise yaw rate.
As was explained above, the vehicle has a yaw angle sensor. If the measured yaw rate {dot over (Ψ)}<sub>Meas </sub>deviates from the {dot over (Ψ)}<sub>Desired </sub>to be reached, an additional moment M<sub>G </sub>around the vertical axis of the vehicle must be applied.
If the measured yaw rate deviates from the yaw rate to be reached to such an extent that the vehicle does not turn sufficiently, a socalled understeering behavior is present. An additional moment, which is counted as negative in this situation, must be applied. It shall cause the vehicle to turn into the curve. This could be achieved in this case by increasing the brake pressure in the righthand wheels of the vehicle.
However, if the vehicle is already being braked by the driver, it may be possible that these wheels already transmit maximum brake force. If this is determined by an electronic evaluation unit, the pressure in the right rear wheel brake is increased such that the wheel runs at slip values in the range of 40% to 80%. Wheel 604 is therefore marked with a “λ.” As was explained above, this leads to a considerable reduction in the lateral force. Consequently, only weak lateral forces are built up on the right rear wheel, as a consequence of which the vehicle swings out with its tail to the left, i.e., a clockwise turning begins. The minimization of the lateral force is maintained until the actual yaw rate {dot over (Ψ)}<sub>Meas </sub>corresponds to the desired {dot over (Ψ)}<sub>Desired </sub>of the vehicle.
FIG. 25b shows the situation of an oversteering vehicle. The vehicle turns around the vertical axis faster than it would correspond to a calculated desired yaw rate. It is proposed that the lateral force on the front left wheel 601 be reduced in this case. This is also done by introducing slip values between 40% and 80% on this wheel. Wheel 601 is therefore marked with a “λ.”
A subprogram that brings about a further reduction in pressure on the front wheel 601 (that is, the outer wheel in the curve for the case of understeering (FIG. 25a)) or on the front wheel 602 (that is, the inner wheel in the curve for the case of oversteering (FIG. 25b)) can be inserted in the control program for both cases. These wheels are marked with “p<sub>min</sub>.” The corresponding actuations are laterally reversed for travel in a curve to the left.
The pressure in the individual wheels can be controlled by determining a coefficient, which describes the relationship between the change in pressure and the calculated additional yaw torque M<sub>G</sub>, for every individual wheel.
These coefficients are a function of parameters that describe the vehicle or the wheel brakes, and of variables which change during travel. These are especially the steering angle δ and the coefficient of friction μ for the road/tire pairing (cf. Section 3.1). A dependence on the longitudinal slip of the corresponding wheel is now additionally introduced for the abovementioned control. The pressure on individual wheels can be prevented from decreasing by defining lower limits for the coefficients, replacing the calculated value of the coefficients with the minimum if the actual value drops below the minimum.
A corresponding algorithm is shown in FIG. 26. The additional yaw torque M<sub>G </sub>is first calculated (program 640). The corresponding changes in the brake force and in the brake pressure are calculated from this moment for the individual wheels (program part 641). The brake pressures determined are compared with thresholds p<sub>th</sub>, which are determined, among other things, by the road/tire coefficientoffriction pairing (block 642). The thresholds p<sub>th </sub>determine whether a further increase in the wheel brake pressure with a simultaneous increase in brake forces is possible. If the pressures to be introduced remain below these limit values, the control is performed according to the procedure mentioned in Section 3.1. If the calculated brake pressures are above these threshold values, the pressures are calculated according to the scheme 644 described above.
4. Priority Circuit
The pressures to be introduced into the wheel brakes are calculated from the additional yaw torque M<sub>G </sub>by means of a distribution logic unit (Section 3).
Based on these pressure values, control signals for inlet and outlet valves are sent by a subordinate pressure control circuit. The actual wheel brake pressures are harmonized with the calculated ones in this subordinate pressure control circuit.
If control signals of other controllers (ABS7, TCS8, EBD9) are to be included as well (Section 1), it is also necessary first to convert their control signals into pressure values by means of a hydraulic model of the wheel brakes stored in the computer.
The pressure requirements of the YTC controller 10 are then related to the pressure requirements of the ABS controller and other controllers. This is done in a priority circuit, which decides what requirements are to be prioritized, and whether averaged pressures are to be sent to the pressure control unit 5 for the wheel brakes. The pressure control unit 5 in turn converts the pressures into valve switching times.
Instead of desired pressures, desired changes in pressure may also be sent to the priority circuit (cf. Section 7).
In this case, the priority circuit 3 sends the changes in pressure Δp to its output according to the rule that the requirement to reduce the pressure on one of the wheels is preferentially satisfied, and the requirement to maintain the pressure in one wheel brake has priority over the requirements to increase the pressure. Thus, the individual requirements on the priority circuit are processed according to the rule that when there is a requirement to reduce the pressure, requirements to maintain the pressure or to increase pressure are ignored. In the same manner, no pressure is increased when maintenance of pressure is required.
5. Priority Circuit with Direct Comparison of Valve Switching Times
Another method can also be used as an alternative to this.
The distribution logic unit calculates valve switching times directly, like the other controllers as well, rather than pressures, from the additional M<sub>G</sub>. The valve switching times of the YTC can thus be compared with the required valve switching times of the ABS. Unlike before, different valve switching times rather than different pressure requirements are then evaluated in the priority circuit.
To obtain valve switching times, the distribution logic unit first calculates changes in pressure to be set for each wheel brake.
Switching times for actuating the individual wheel brakes are calculated from the changes in pressure by means of a downstream, nonlinear control element.
This nonlinear control element may be, e.g., a counter.
This counter converts the preset changes in pressure into cycle counts. To do so, the loop time T<sub>0 </sub>is divided into approximately 3 to 10 switching intervals (cycles). The maximum number of cycles per loop time is a fixed quantity, which is determined according to the quality of control to be reached.
How long a valve within a loop time is to be actuated is determined by the calculated cycle count.
Since there are, in general, two valves per wheel brake, with one valve (inlet valve) regulating the feed of the pressure medium to the wheel brake, and the other valve (outlet valve) regulating the release of the pressure medium from the wheel brake, a total of eight signals are to be generated.
Theses cycle counts are sent to the priority circuit, which receives the cycle counts of other controllers in additional channels.
The priority circuit decides which controller is to be given preference, i.e., which cycle count is taken over for the actual valve control.
The response of the vehicle to the brake forces generated by the actuation of the wheel brakes is a changed yaw rate. This is detected by the YTC controller 10, which will again determine a new additional yaw torque.
Consequently, brake pressures are not calculated or set at any point of the control circuit. Therefore, the control algorithms need no information on the wheel brake, and, in particular, no information on the relationship between the volume received by the wheel brakes and the resulting brake pressures.
One possibility of calculating the cycle times is explained on the basis of FIG. 27.
Brake pressures, which are to be built up in the individual wheel brakes, are calculated from the additional yaw torque M<sub>G </sub>by way of the distribution logic unit 700. How this is done can be found described in Sections 3.1. and 3.2. As a result of the calculation within the distribution logic unit, there are four pressure values p<sub>1 </sub>through p<sub>4 </sub>for a fourwheel vehicle. These variables must be converted into switching times for the valves, which control the feed of pressure medium (pressure buildup) and the release of the pressure medium (pressure reduction) and from the wheel brakes. As was mentioned above, the switching times for the valves are calculated from the change in the preset pressure value rather than from the absolute values of the preset pressure value. Each value p<sub>n </sub>(n=1 through 4) is therefore sent to a shift register 701. The current value is written to the first register place 702. The previous value from the first register place 702 is received in the second register place 703, so that the pressure requirement from the preceding calculation loop is written there. This value is designated by p<sub>n</sub>*.
The current pressure requirement is read from the first register place 702 in the next step 705. If this value is 0 or lower than a minimum, the program branches into a loop 706, with which it shall be ensured that so much pressure medium is removed from the wheel brake that the pressure becoming established becomes zero. To do so, the inlet valve is closed and the outlet valve is opened for at least one loop time T<sub>0</sub>.
If the current required pressure value is above this minimum, the difference of the two register values 702 and 703 is formed. This is done in the subtractor 707. The calculated change in pressure Δp may be either greater or less than 0. If it is greater than 0, the pressure must be increased in the corresponding wheel brake. If it is less than 0, the pressure must be reduced in the corresponding wheel brake. In the case of a pressure buildup, the program runs through the righthand decision path 710. Taking the pressure difference to be set and the pressure requirement or, if corresponding signals are present, based on the actual pressure in the wheel brake, an opening time Δt<sub>in </sub>is calculated for the inlet valve. The opening time Δt<sub>out </sub>of the outlet valve is set to zero. Conversely (decision path 711), the opening time Δt<sub>in </sub>of the inlet valve is set to zero if a reduction in pressure is required, while the opening time Δt<sub>out </sub>of the outlet valve is calculated from the required pressure difference and the actual pressure in the wheel brake or the required pressure, which is written in the first register place 702.
As a rule, there is a linear relationship between the opening time Δt and the intended change in pressure Δp.
As was explained, the calculation is performed with cycle counts rather than with the opening times. This is explained in greater detail in the diagram in FIG. 28. The abovedescribed calculations are performed at constant time intervals (loop time T<sub>0</sub>), and the control signals for the valves of the wheel brakes in the next loop are set as the result of a calculation. One loop time T<sub>0 </sub>is approximately 3 msec.
Depending on how fine the control is to operate, each loop time T<sub>0 </sub>is divided into N time intervals.
The diagram in FIG. 28 shows a division into 6 steps. The switching times for the valves are no longer issued as time variables, but as the number of cycles within one loop, during which the valve is to be opened. As can be determined from FIG. 28, an opening time of 1.5 msec is obtained, e.g., for n=3.
Should the required opening time be longer than the loop time, n is set at the corresponding maximum value N (to 6 in the example shown).
This calculation is performed for each wheel brake, i.e., four times for a fourwheel vehicle. The calculations may be performed simultaneously or consecutively. As a result, eight values are available; four values for inlet valves and four values for outlet valves. These values are sent to a modified priority circuit 720. The switching time requirement, likewise expressed in cycle times, of an ABS controller and additional controllers are sent to this priority circuit 720 as well.
This actuation is performed such that a change in the pressure in the wheel brakes is obtained. The pressure forces and consequently the moments exerted on the vehicle will thus change. Thus, a change is obtained in the variables which describe the driving dynamics of the vehicle. These are directly or indirectly detected by sensors and are in turn sent to the calculation.
This again leads to a changed moment requirement, which, as was described above, is converted into new control signals for the valves.
The calculation of the pressure differences to be set is based on the pressure requirements from the preceding calculation loop. However, these do not have to have been actually set, so that the actual pressures in the wheel brakes differ from the corresponding calculated pressure requirements. It is therefore necessary to adjust the actual pressure in the wheel brake to the pressure requirements in certain situations. This can be done in the simplest manner when the pressure requirement is zero, i.e., the distribution logic unit 700 requires a value that corresponds to the pressure zero in a wheel brake. The difference from the preceding value is not formed, and the control signals are not derived from this in such a case, but it is branched off in step 705 into the loop 706 for calculating the switching times, and this loop is to ensure that a pressure value of zero is indeed set. This is done by setting the switching time Δt<sub>out </sub>for the outlet valve to at least the loop time T<sub>0</sub>.
It may also become necessary to send corresponding information to the priority circuit 720, so that this time requirement, which is to lead to zero pressure in a wheel brake, will not be superimposed by preset values of the other controllers. In addition, it can be determined in this information that the reduction in pressure shall take place over several loop times, so that it is ensured that a complete pressure reduction will indeed take place.
6. Wheel Brake Pressure Detection
The DSC pressure controller described up to Section 4 provides brake pressure values for the wheel brakes as a result. These preset values must be put into practice. One method is to measure the pressures in the wheel brakes and to compare them with the preset values. A pressure controller that operates according to the usual laws adjusts the wheel brake pressure to the predetermined desired value. This procedure requires one pressure sensor per wheel brake, i.e., four pressure sensors for a fourwheel vehicle.
Attempts will be made, in general, even for cost reasons to make do with as few sensors as possible. In addition, each sensor represents another potential source of disturbance. The failure of one sensor may lead to the necessity of switching off the entire control system.
It is therefore proposed that an evaluation system be provided, which derives a pressure variable that corresponds to the pressure in the wheel brakes on the basis of data available from the already existing sensors. The following concept is proposed for doing so.
As was explained above, the pressure in each wheel brake is controlled by two valves. The inlet valve controls the feed of the pressure medium, while the outlet valve controls the release of the pressure medium.
The signals sent by a pressure controller are therefore control times which indicate how long a valve shall be opened or closed. One loop time is divided into a fixed number of time intervals (cycles). The control times can thus be represented as a cycle count, which indicates over how many time intervals a valve shall be opened or closed.
The basic consideration is that these control signals shall be sent not only to the wheel brakes, but as calculated variables also to a vehicle model. The real vehicle responds to the brake pressures introduced, and a certain velocity v of the center of gravity and wheel speeds ω<sub>i </sub>of the individual wheels will become established. The velocity of the vehicle is not directly measured, but it is also derived from the speeds ω<sub>i </sub>of the individual wheels in special calculation steps. They are therefore called the reference velocity v<sub>Ref</sub>.
Corresponding values can also be simulated within one vehicle model.
A correcting variable for the pressure in the individual wheel brakes can be determined from a comparison of the actual values of ω<sub>i</sub>, v<sub>Ref </sub>with the calculated values of ω<sub>i </sub>and v<sub>Ref </sub>or on the basis of the values of ω<sub>i </sub>and v<sub>Ref </sub>estimated on the basis of the vehicle model, and a pressure calculated by way of a hydraulic model can be modified by means of the correcting variable, so that a better estimate of the wheel brake pressures can be given.
The general structure just described is explained in greater detail in FIG. 29.
A pressure control unit, which has number 5 in FIG. 1, is designated by 800. The pressure control unit calculates control times for the valves of the wheel brakes from a first value 801, which characterizes the pressure to be set, and from a second value 802, which marks an existing, estimated or measured pressure in the wheel brake. The control times are represented as an output variable 803 here. The vehicle is designated by 810. This is to illustrate that the vehicle responds to forces which are caused by the pressure set in the wheel brakes. The speeds ω<sub>i </sub>of the individual wheels will change now as well.
Wheel sensors, which detect the speeds of the wheels, so that the ω<sub>i </sub>values are immediately available, shall also belong to the vehicle 810.
An evaluation unit ω<sub>i </sub>also belongs to the vehicle 810; this evaluation unit usually represents a partial area of an ABS controller, which calculates a socalled reference velocity v<sub>Ref</sub>, which is to correspond to the actual velocity of the vehicle, form the wheel speeds ω<sub>i </sub>of the individual wheels under certain boundary conditions.
A slip λ<sub>i </sub>can be calculated for each wheel from the individual wheel speeds and the vehicle reference velocity.
The values ω<sub>i</sub>, v<sub>Ref </sub>are available as output values 811. The slip λ<sub>i </sub>is available as the value 812.
The calculation model used is designated as a whole by 820. It contains three sub models, namely,
a hydraulic model 821,
a vehicle model 822, and
a tire model 823.
In two approximation formulas, the hydraulic model 821 describes the relationship between the brake pressure p and the volume V enclosed in the wheel brake and the change ΔV in volume when the inlet or outlet valve is opened for a certain time.<FORM>p=a*V+b*V<sup>2</sup> F 6.1 </FORM><FORM>ΔV=±c*t<sub>in/out</sub>*√{square root over (Δp)} F 6.2 </FORM>
The parameters a, b and c are variables which describe the brake system and are stored as values in corresponding memories. P describes the current pressure in the wheel brake. V describes the current volume enclosed in the wheel brake.
Δp is measured either across the inlet valve or across the outlet valve; the difference between a pressure source and p is determined in the case of measurement across the inlet valve, while the difference between p and the pressure in a tank, which is usually 1 bar and therefore can be ignored, is determined in the case of measurement across the outlet valve.
If it is assumed that the pressure in the wheel brakes and the enclosed volume can be set to zero at the beginning of a control, the change in volume and, hence, the change in pressure in the individual wheel brakes can be reconstructed by monitoring the valve opening times.
At any rate, it is clear that the formulas shown can describe the actual conditions only very approximately, so that a corresponding correction is necessary. In model 822, the vehicle is described, in general, by a rigid body, which stands on a plane in four contact points (tire contact points).
The body can move in parallel to the plane, i.e., in the x and y directions, and rotate around its center of gravity, with the axis of rotation being at right angles to the plane of movement.
The forces acting on the body are brake forces in the contact surface of the tires and air resistance forces.
The wheel loads F<sub>z,v </sub>and F<sub>z,h </sub>which are directed perpendicular to the plane can be calculated based on these considerations as follows:<maths id="MATHUS00024" num="24"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><msub><mi>F</mi><mrow><mi>z</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>=</mo><mfrac><mrow><mrow><mi>m</mi><mo>*</mo><mi>g</mi><mo>*</mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo>+</mo><mrow><mi>h</mi><mo>*</mo><mrow><mo>(</mo><mrow><mrow><mo></mo><msub><mi>F</mi><mrow><mi>x</mi><mo>,</mo><mi>v</mi></mrow></msub></mrow><mo></mo><msub><mi>F</mi><mrow><mi>x</mi><mo>,</mo><mi>h</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow></mrow><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mrow><mi>m</mi><mo>*</mo><mi>g</mi><mo>*</mo><msub><mi>l</mi><mi>h</mi></msub></mrow><mo></mo><mrow><mi>h</mi><mo>*</mo><mi>m</mi><mo>*</mo><msub><mover><mi>v</mi><mo>.</mo></mover><mi>ref</mi></msub></mrow></mrow><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>6.3</mn><mo></mo><mi>a</mi></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><msub><mi>F</mi><mrow><mi>z</mi><mo>,</mo><mi>h</mi></mrow></msub><mo>=</mo><mfrac><mrow><mrow><mi>m</mi><mo>*</mo><mi>g</mi><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mo>+</mo><mrow><mi>h</mi><mo>*</mo><mrow><mo>(</mo><mrow><mrow><mo></mo><msub><mi>F</mi><mrow><mi>x</mi><mo>,</mo><mi>v</mi></mrow></msub></mrow><mo></mo><msub><mi>F</mi><mrow><mi>x</mi><mo>,</mo><mi>h</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow></mrow><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mrow><mi>m</mi><mo>*</mo><mi>g</mi><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mo></mo><mrow><mi>h</mi><mo>*</mo><mi>m</mi><mo>*</mo><msub><mover><mi>v</mi><mo>.</mo></mover><mi>ref</mi></msub></mrow></mrow><mrow><msub><mi>l</mi><mi>v</mi></msub><mo>+</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mfrac></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>6.3</mn><mo></mo><mi>b</mi></mrow></mtd></mtr></mtable></math></maths>
Such a model is usually sufficient for performing the desired pressure correction. The model can, of course, be refined, if necessary. For the further calculation, the model provides essentially the loads F<sub>x </sub>of the tire contact surfaces as a function of the deceleration of the center of gravity. The wheel is considered to be a rotatable disc, which has a certain moment of inertia.<maths id="MATHUS00025" num="25"><math overflow="scroll"><mtable><mtr><mtd><mrow><mover><mi>ω</mi><mo>.</mo></mover><mo>=</mo><mfrac><mrow><mrow><msub><mi>R</mi><mi>rad</mi></msub><mo>*</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mo></mo><msub><mi>M</mi><mi>Br</mi></msub></mrow><mi>Θ</mi></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>6</mn><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>.4</mi></mrow></mtd></mtr></mtable></math></maths>
The decelerating torques acting on the wheel are determined linearly from the wheel brake pressure.<FORM>M<sub>Br</sub>=C<sub>Br</sub>*p F 6.5 </FORM>
It is assumed in the tire model that the utilization of the frictional connection, f, namely, the ratio of the braking force to the wheel load, changes linearly with the slip of the wheel.<FORM>F<sub>x</sub>˜λ*F<sub>z</sub> F 6.6 </FORM>
The equations given make it possible to calculate the wheel speed of each wheel and the reference velocity of the vehicle model.
These values can be compared with the actual values 811. This is done at the reference point 830.
Taking a correction factor k into account, an additional volume can be determined from the difference between the measured and estimated speeds of each wheel.
This additional pressure medium volume ΔV is added to the calculated desired volume to obtain the new desired volume, from which a wheel brake pressure, which corresponds to the actual wheel brake pressure relatively accurately, can be derived according to formula F 6.1.
The accuracy of the estimation depends, of course, on the correction factor k, which may have to be determined by experiments in advance.
This factor differs from one vehicle to the next, and it also depends, among other things, on how well the vehicle model describes the actual conditions.
The additional volume may also include a tolerance volume, with which the fact that the volume throughput through the valves is not proportional to the switching times is taken into account. The opening cross section of the valve increases or decreases only slowly during the opening and closing of a valve, so that only a reduced volume will flow during the time intervals in which the actual opening cross section still increases toward or decreases from the full opening crosssection.
7. Substitution of a Yaw Rate Meter
The yaw rate is a particularly distinctive variable for the abovementioned control, because it is used as a controlled variable, whose deviation Δ{dot over (Ψ)} is to be minimized. However, as will be described below, other controlled variables may be advantageously used as well. The following designations will be used in this section for simplification:
Δ{dot over (Ψ)}=g<sub>I </sub>as the measured actual value of the yaw rate,
Δ{dot over (Ψ)}<sub>Meas</sub>=ġ<sub>I </sub>as the measured actual value of the yaw acceleration,
d/dtΔ{umlaut over (Ψ)}<sub>Meas</sub>=&guml;<sub>I </sub>as the measured actual value of the change in yaw acceleration (yaw angle pressure).
This also applies analogously to the desired values according to FIG. 9, which are always marked with the subscript “s.”
The measured yaw rate in FIG. 12 is usually determined by means of a yaw rate sensor 321, which issues the output signal g<sub>I</sub>. However, such known yaw rate sensors with direct issuance of the yaw rate are of a rather complicated design and therefore very expensive. This is also true of the downstream comparison unit and the controller belonging to the control circuit. It is therefore desirable in this respect to seek a way out and to offer simpler sensor systems and a controller of a simpler design.
FIG. 13 shows the sketch of the mode of operation of a novel sensor 321, which has a first lateral accelerometer 322 and a second lateral accelerometer 323. The two acceleration meters 322, 323 are arranged on the longitudinal axis of the vehicle above the front axle and the rear axle, respectively. The lateral accelerometers may be arranged, in principle, at any point outside the center of gravity SP, in which case a corresponding conversion is performed.
FIG. 15 indicates the rectangular outline 324 of a vehicle with its tires 325 and sensors. Based on this arrangement, the front lateral accelerometer 322 measures the lateral acceleration a<sub>qv </sub>at the level of the front axle 326, and the rear lateral accelerometer 323 measures the lateral acceleration a<sub>qh </sub>at the level of the rear axle 327.
The two lateral accelerometers are able to furnish a variable that depends on the yaw rate. It can be shown from mathematical deductions that the yaw acceleration and the lateral acceleration a<sub>trans </sub>of the center of gravity SP can be determined from the measurement results of the lateral accelerometers as follows:<maths id="MATHUS00026" num="26"><math overflow="scroll"><mtable><mtr><mtd><mrow><mover><mi>Ψ</mi><mi>¨</mi></mover><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mi>qh</mi></msub><mo></mo><msub><mi>a</mi><mi>qv</mi></msub></mrow><mrow><msub><mi>l</mi><mi>h</mi></msub><mo>+</mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mi>lat</mi></msub><mo>≈</mo><mrow><mrow><mrow><mo></mo><msub><mover><mi>v</mi><mo>.</mo></mover><mi>ref</mi></msub></mrow><mo>*</mo><mi>β</mi></mrow><mo></mo><mfrac><mrow><mrow><msub><mi>a</mi><mi>qh</mi></msub><mo>*</mo><msub><mi>l</mi><mi>v</mi></msub></mrow><mo>+</mo><mrow><msub><mi>a</mi><mi>qv</mi></msub><mo>*</mo><msub><mi>l</mi><mi>h</mi></msub></mrow></mrow><mrow><msub><mi>l</mi><mi>h</mi></msub><mo>+</mo><msub><mi>l</mi><mi>v</mi></msub></mrow></mfrac></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.2</mn></mrow></mtd></mtr></mtable></math></maths>
As is apparent from FIG. 13, l<sub>v</sub>, l<sub>h </sub>are the distances between the respective lateral accelerometers 322, 323, on the one hand, and the center of gravity SP, on the other hand, while v is the velocity of the vehicle, and β is the sideslip angle. The yaw acceleration ġ<sub>I </sub>can thus be determined from the lateral accelerations and the distances of the acceleration meters 322, 323.
It is therefore proposed that the yaw rate acceleration ġ<sub>I </sub>be used instead of the yaw rate proposed in the previous sections, or it is also possible to perform a linear weighting of the individual input values for the comparison unit, similarly to the priorart condition control. The yaw rate g and the sideslip angle β can be calculated from the yaw angle pressure ġ<sub>I </sub>and the velocity of the sideslip angle {dot over (β)} by means of a bandlimited integration or a firstorder, scaled, lowpass filter in order to obtain variables whose dimension corresponds to the output variables of the vehicle reference model 302 (Section 2.3.1.) from sensor 321.
For the bandlimited integration:<maths id="MATHUS00027" num="27"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>G</mi><mo></mo><mrow><mo>(</mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></mrow><mo>=</mo><mfrac><mrow><msup><mi>λ</mi><mn>2</mn></msup><mo>*</mo><mfrac><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></mfrac><mo>*</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo></mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow><mo>*</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo></mo><mrow><mi>λ</mi><mo>*</mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup></mrow></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.3</mn></mrow></mtd></mtr></mtable></math></maths>applies, while the following dependence is obtained by using a lowpass filter:<maths id="MATHUS00028" num="28"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>G</mi><mo></mo><mrow><mo>(</mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></mrow><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mn>1</mn></msub><mo>*</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo></mo><mi>λ</mi></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo></mo><mrow><mi>λ</mi><mo>*</mo><msup><mi>z</mi><mrow><mo></mo><mn>1</mn></mrow></msup></mrow></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.4</mn></mrow></mtd></mtr></mtable></math></maths>
The velocity of the sideslip angle is obtained after evaluating the equation<FORM>a<sub>lat</sub>=v+({dot over (Ψ)}+{dot over (β)}) F 7.5 </FORM>
Thus, it is seen that even though a priorart yaw rate meter can be replaced by using two lateral accelerometers, the measures just described must be taken to transform the yaw acceleration into the yaw rate. However, the measures just described must be taken to transform the yaw acceleration into the yaw rate. After forming Δg and Δġ, the control law unit 16 from FIG. 1 can follow unchanged. The moment M<sub>G </sub>thus calculated is additionally converted in the control law unit 16 into a change in moment M by a derivation with respect to time.
However, it is more expedient under certain circumstances to pass over to a nonlinear control according to FIG. 17, in which the yaw acceleration ġ is sent to the comparison unit 303 both as an actual value and as a desired value as a result from the vehicle reference model 302. To do so, corresponding derivatives must be formed within the vehicle reference model.
As a consequence, the deviation of the yaw acceleration Δġ, rather than the yaw rate difference Δg, is present at the output of the comparison unit 303 and is sent as an input variable to the control law unit 16. Furthermore, as is apparent from FIG. 15, the velocity of the sideslip angle {dot over (β)} can be additionally sent to the yaw torque control law unit 16 for the more accurate determination of the change in the moment.
As was mentioned in connection with FIG. 14, it is possible to abandon an additional yaw torque M<sub>G </sub>as an output signal of the control law unit 16, and to use the change in moment &Mdot;, as the output signal, instead. The change in moment, &Mdot;, i.e., the derivative of the additional yaw torque M<sub>G</sub>, is converted into individual changes in pressure in a modified distribution logic unit. This means that the changes in pressure are distributed among the individual wheel brakes such that the desired additional yaw torque M<sub>G </sub>is obtained, on the whole. Details of this will be described below in connection with FIG. 16.
It should be borne in mind that at the same time, there may be a certain pressure distribution in the wheel brakes due to the driver actuating the brake. It is more favorable in this case to determine the moment M<sub>G </sub>by integrating the change in moment &Mdot;, whereafter the pressure differences that must be brought about with respect to the pressure occurring in every individual wheel brake can be directly determined from the moment M<sub>G</sub>. The abovedescribed advantageous variant, in which the derivatives of the controlled variables used in Sections 1 through 3 are used, may also be combined with the distribution logic unit according to Section 3. Two control principles are available here; one of them yields an additional yaw torque M<sub>G</sub>, and the other a change in the additional yaw torque &Mdot; as a preset value. Switching over between the principles may be provided for. Switching over to the respectively other control principle must be performed especially when the other calculation of additional controlled variables (sideslip angle, etc.) according to one principle cannot be performed with sufficient accuracy (cf., e.g., Section 2.2.2.). It should also be noted that Δ&guml; can also be sent as a correcting variable to the control law unit 16 according to FIG. 15, in addition to Δġ.
Besides adapting amplifiers k1, k2, k3, two threshold value switches S2, S3 are shown in the control law unit 16 according to FIG. 15; these threshold value switches are to improve the control behavior within the control law unit 16 and to optimally adapt the influence of the introduced variables to the ideal control behavior as a function of the velocity. The amplifiers k1 through k3 have a comparable task. The individual values are then added in an adder and sent as an output signal to the YTC controller 10. General explanations with respect to the control law unit which apply here can be found in Section 2.4.
How the pressure preset values at the output of the controllers 7, 8, 9 are linked with the pressure preset value of a distribution logic unit 2 in a priority circuit 3 was shown in connection with FIG. 1. The use of pressure preset values always requires a corresponding prior conversion in the devices that issue these preset values. The effort involved in the exchange of information between the program modules of the control circuit can be simplified by the measures described below.
The control circuit for controlling the driving stability according to FIGS. 9 and 14 is shown in an even more simplified form in FIG. 16; the designations introduced there are maintained.
The YTC controller 10 according to FIG. 1 is modified here inasmuch as the change &Mdot; in the additional yaw torque M<sub>G</sub>, which is sent to the distribution logic unit 2 together with the pressure distribution on the brakes desired by the driver (desire to brake), is present at the output. FIG. 12 is referred to for the calculation of &Mdot;.
The distribution logic unit 2 has a logic block 340 and a pressure gradient circuit 341. The essential task of the logic block 340 is to ensure that despite the intervention of the driving stability control, the vehicle as a whole is not braked more strongly than is desired by the driver by presetting a pressure signal at the input of the distribution logic unit 2. This is to prevent instabilities from being additionally introduced by the driving stability control system. Consequently, when a brake pressure is provided on a wheel based on the driver's desire to brake, and, on the other hand, a pressure buildup on one or two wheels is required by way of the DSC controller and a reduction in pressure on the opposite wheels is required in order to reach the additional yaw torque, there may be mutually contradictory requirements with respect to the individual wheels, namely, a pressure buildup with a simultaneous reduction in pressure. Regarding other wheels, it may be required to increase the pressure not only based on the driver's desire to brake, but at the same time also based on the stability control. The logic block ensures that the brake pressure is first reduced in the corresponding wheels, after which an increase in brake pressure beyond the driver's desire up to a certain limit value can take place. It is thus ensured that the average brake force will not become greater, considering all wheels and taking the additional torque brought about by the DSC control into account, than that desired by the driver.
As was explained in Section 3.2., a specific increase in the longitudinal slip λ on one wheel can be used to reduce the lateral forces, while the brake force is preserved in the longitudinal direction. Consequently, a yaw torque can thus be generated without the deceleration of the vehicle decreasing.
The changes in pressure ΔP<sub>xx </sub>on the individual wheels xx are calculated in the pressure gradient circuit 341 of the distribution logic unit 2 on the basis of predetermined constants c<sub>xx </sub>and the change in moment &Mdot;, and the difference between the brake pressure desired by the driver, P<sub>Brake</sub>, and the brake pressure actually measured, P<sub>xxist</sub>, is also included in the calculation. Thus, the following equation applies<maths id="MATHUS00029" num="29"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><mi>Δ</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><msub><mi>P</mi><mi>xx</mi></msub></mrow><mo>=</mo><mrow><mrow><mo>±</mo><mfrac><mrow><msub><mi>d</mi><mi>xx</mi></msub><mo>*</mo><mover><mi>M</mi><mo>.</mo></mover></mrow><mrow><mo>∑</mo><msubsup><mi>c</mi><mi>xx</mi><mn>2</mn></msubsup></mrow></mfrac></mrow><mo>+</mo><mrow><msub><mi>g</mi><mi>l</mi></msub><mo>*</mo><mrow><mo>(</mo><mrow><msub><mi>P</mi><mi>Brake</mi></msub><mo></mo><msub><mi>P</mi><mi>xxist</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mi>and</mi><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mi>xx</mi><mo>∈</mo><mrow><mo>[</mo><mrow><mi>vr</mi><mo>,</mo><mi>vl</mi><mo>,</mo><mi>hr</mi><mo>,</mo><mi>hl</mi></mrow><mo>]</mo></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mrow><mi>and</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><msub><mi>g</mi><mi>l</mi></msub></mrow><mo>=</mo><mrow><mi>proportionally</mi><mo></mo><mrow><mstyle><mtext> </mtext></mstyle><mo></mo><mstyle><mtext> </mtext></mstyle></mrow><mo></mo><mrow><mi>factor</mi><mo>.</mo></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.6</mn></mrow></mtd></mtr></mtable></math></maths>
The actual brake pressure P<sub>xxist </sub>is determined either by a pressure gauge at the corresponding wheel, or it is calculated by way of a brake model, which follows the changes in pressure specified for the wheel and is therefore an image of the pressure occurring on the wheel (Section 6).
The pressure requirements calculated are sent to a priority circuit 3 and they are evaluated there (See section 4, above).
The above description presupposes that pressure gradients were directly processed in the priority circuit. However, this is not necessary. It is also possible to process valve switching times Δt in the priority circuit 3 (Section 5). However, a valve switching time circuit 343 must be inserted in this case between the distribution logic unit 2 and the priority circuit 3, and valve switching times Δt will be sent by the other controllers 7, 8, 9 as well. The priority circuit now processes the valve switching times Δt entered according to a corresponding scheme, as was described in Section 4 for the brake pressures. The output variables of the priority circuit are valve switching times. The required changes in pressure Δt<sub>xx </sub>of the individual wheels xx are converted into valve switching times Δp according to the equation<FORM>S<sub>xx</sub>=Krp<sub>xxist</sub>·Δp<sub>xx</sub> F 7.7 </FORM>
Here, Kr<sub>xx </sub>is a gain factor that depends on the actual pressure of the individual wheels and is calculated during pressure buildup according to the following rule:<maths id="MATHUS00030" num="30"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>Kr</mi><mi>xx</mi></msub><mo></mo><mrow><mo>(</mo><msub><mi>P</mi><mi>xxist</mi></msub><mo>)</mo></mrow></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>Dv</mi><mi>up</mi></msub><mo>*</mo><msub><mi>T</mi><mn>0</mn></msub><mo>*</mo><msqrt><mrow><msup><mi>a</mi><mn>2</mn></msup><mo>+</mo><mrow><mn>4</mn><mo>*</mo></mrow></mrow></msqrt><mo></mo><mi>b</mi><mo>*</mo><msub><mi>P</mi><mi>xxist</mi></msub><mo>*</mo><msqrt><mrow><mn>160</mn><mo></mo><msub><mi>P</mi><mi>xxist</mi></msub></mrow></msqrt></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.8</mn></mrow></mtd></mtr></mtable></math></maths>while<maths id="MATHUS00031" num="31"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>Kr</mi><mi>xx</mi></msub><mo></mo><mrow><mo>(</mo><msub><mi>P</mi><mi>xixst</mi></msub><mo>)</mo></mrow></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>Dv</mi><mi>down</mi></msub><mo>*</mo><msub><mi>T</mi><mn>0</mn></msub><mo>*</mo><msqrt><mrow><msup><mi>a</mi><mn>2</mn></msup><mo>+</mo><mrow><mn>4</mn><mo>*</mo></mrow></mrow></msqrt><mo></mo><mi>b</mi><mo>*</mo><msub><mi>P</mi><mi>xxist</mi></msub><mo>*</mo><msqrt><msub><mi>P</mi><mi>xxist</mi></msub></msqrt></mrow></mfrac></mrow></mtd><mtd><mrow><mi>F</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mn>7.9</mn></mrow></mtd></mtr></mtable></math></maths>applies to a reduction in pressure. Here, xx is again a subscript indicating the position of the individual wheels.