[go: up one dir, main page]

JP2962025B2 - Body slip angle estimation device - Google Patents

Body slip angle estimation device

Info

Publication number
JP2962025B2
JP2962025B2 JP4020592A JP2059292A JP2962025B2 JP 2962025 B2 JP2962025 B2 JP 2962025B2 JP 4020592 A JP4020592 A JP 4020592A JP 2059292 A JP2059292 A JP 2059292A JP 2962025 B2 JP2962025 B2 JP 2962025B2
Authority
JP
Japan
Prior art keywords
vehicle
slip angle
vehicle body
yaw rate
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP4020592A
Other languages
Japanese (ja)
Other versions
JPH05185948A (en
Inventor
浩 稲垣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to JP4020592A priority Critical patent/JP2962025B2/en
Publication of JPH05185948A publication Critical patent/JPH05185948A/en
Application granted granted Critical
Publication of JP2962025B2 publication Critical patent/JP2962025B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Steering Control In Accordance With Driving Conditions (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は車両重心点における車体
横すべり角を推定する装置に関するものであり、特にそ
れの推定精度を向上させる技術に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an apparatus for estimating a side slip angle of a vehicle at the center of gravity of a vehicle, and more particularly to a technique for improving the accuracy of the estimation.

【0002】[0002]

【従来の技術】車両の操縦性および安定性を車両の実際
の運動状態量との関係において適正に制御したいという
要望があり、その運動状態量の一つに車両重心点におけ
る車体横すべり角が存在する。この車体横すべり角は一
般に、それ自身を直接取得するのは比較的困難であるた
め、それに関連するパラメータを検出してそれから間接
に取得することが行われている。
2. Description of the Related Art There is a demand for appropriately controlling the maneuverability and stability of a vehicle in relation to the actual amount of motion of the vehicle, and one of the amounts of motion is the vehicle body slip angle at the center of gravity of the vehicle. I do. Generally, it is relatively difficult to directly obtain the vehicle body side slip angle itself, and therefore, a parameter related thereto is detected and then obtained indirectly.

【0003】車体横すべり角を間接に取得する技術の一
例が特開昭62−105776号公報に記載されてい
る。これは、車両の横すべり運動に関して、 β′=Gy /V−γ ただし、 β′:車体横すべり角βの時間微分値 Gy :車両重心点における横加速度 V:車体の走行速度である車速 γ:ヨーレート なる運動方程式が成立するという事実に基づき、その式
の右辺の値を積分することによって車体横すべり角βを
推定する技術である。
An example of a technique for indirectly obtaining a vehicle body slip angle is described in Japanese Patent Application Laid-Open No. Sho 62-105776. This is related to the side slip motion of the vehicle. Β ′ = G y / V−γ where β ′: Time derivative of the vehicle side slip angle β G y : Lateral acceleration at the center of gravity of the vehicle V: Vehicle speed γ which is the running speed of the vehicle : Yaw rate Based on the fact that the equation of motion holds, the vehicle body side slip angle β is estimated by integrating the values on the right side of the equation.

【0004】[0004]

【発明が解決しようとする課題】この従来の車体横すべ
り角推定技術(以下、単に従来技術という)は、具体的
には、横加速度センサからの出力信号に基づく横加速度
の検出値を車速センサからの出力信号に基づく車速Vの
検出値で割った値から、ヨーレートセンサからの出力信
号に基づくヨーレートγの検出値を差し引いた値を時間
に関して積分することによって車体横すべり角βを逐次
推定するものである。
This conventional vehicle body slip angle estimating technique (hereinafter simply referred to as "prior art") specifically, detects a lateral acceleration detection value based on an output signal from a lateral acceleration sensor from a vehicle speed sensor. The vehicle body slip angle β is successively estimated by integrating the value obtained by subtracting the detected value of the yaw rate γ based on the output signal from the yaw rate sensor from the value divided by the detected value of the vehicle speed V based on the output signal of the vehicle with respect to time. is there.

【0005】しかし、車両に搭載されるセンサからの出
力信号には普通に、温度変化等に基づいてドリフトが発
生する。そのため、この従来技術を採用する場合には、
そのドリフトに基づいて検出値に発生する誤差も積分さ
れて車体横すべり角βが推定されてしまい、その推定精
度を十分には高めることができないという問題がある。
[0005] However, an output signal from a sensor mounted on the vehicle usually drifts due to a temperature change or the like. Therefore, when adopting this conventional technology,
An error generated in the detected value based on the drift is also integrated to estimate the vehicle body slip angle β, and there is a problem that the estimation accuracy cannot be sufficiently improved.

【0006】本発明はこの問題を解決することを課題と
して為されたものである。
[0006] The present invention has been made to solve this problem.

【0007】[0007]

【課題を解決するための手段】この課題を解決するため
に本発明の要旨は、車両重心点における車体横すべり角
を推定する装置を、(a) 車体に発生する横加速度を検出
する横加速度センサと、(b) 車体の走行速度である車速
を検出する車速センサと、(c) 車体に発生するヨーレー
トを検出するヨーレートセンサと、(d) それらに接続さ
れて車両重心点における車体横すべり角βを逐次推定す
る車体横すべり角推定手段であって、車両重心点の回り
に発生するヨーイングモーメントYMを車両のヨー慣性
モーメントIZ で割った変数であるYM/IZ の現在値
を取得するとともに、車両の横すべり運動とヨーイング
運動とから成る車両の平面運動に関して想定される、車
体横すべり角βとヨーレートγとをそれぞれ状態変数と
し、車両重心点における横加速度Gy を車速Vで割った
変数であるGy /Vと、YM/IZ から車体横すべり角
βの前回推定値を差し引いた変数とをそれぞれ入力変数
とし、ヨーレートγを出力変数とする可観測な線形定係
数システムにおいてオブザーバを用いて、横加速度セン
サ,車速センサおよびヨーレートセンサの各々からの出
力信号に基づく横加速度Gy ,車速Vおよびヨーレート
γの現在値と、取得したYM/IZ の現在値と、車体横
すべり角βの前回推定値とに基づいて車体横すべり角β
の今回値を推定するものとを含むものとしたことにあ
る。
SUMMARY OF THE INVENTION In order to solve this problem, the gist of the present invention is to provide a device for estimating a vehicle body slip angle at the center of gravity of a vehicle, comprising: (a) a lateral acceleration sensor for detecting a lateral acceleration generated in the vehicle body; (B) a vehicle speed sensor that detects a vehicle speed that is the traveling speed of the vehicle body, (c) a yaw rate sensor that detects a yaw rate generated on the vehicle body, and (d) a vehicle body slip angle β at the vehicle center of gravity connected to them. a vehicle slip angle estimating means for sequentially estimating acquires the current value of a variable by dividing the yawing moment YM generated around the vehicle center of gravity in the yaw inertia moment I Z of the vehicle YM / I Z, The vehicle body side slip angle β and the yaw rate γ assumed for the plane movement of the vehicle including the vehicle side slip movement and the yawing movement are set as state variables, respectively. The rate G y and the G y / V is a variable divided by the vehicle speed V, the from YM / I Z respectively input variables and variables obtained by subtracting the previous estimated value of the vehicle slip angle beta, soluble to output variable yaw rate γ Using an observer in the observed linear constant coefficient system, the current values of the lateral acceleration G y , the vehicle speed V and the yaw rate γ based on the output signals from the lateral acceleration sensor, the vehicle speed sensor and the yaw rate sensor, and the obtained YM / I Z Of the vehicle body slip angle β based on the current value of
And that for estimating the current value.

【0008】なお、ここにおける「車体横すべり角推定
手段」はYM/IZ の現在値を取得する方式に種々のも
のを採用することができる。例えば、ヨー慣性モーメン
トIZ は不変であるとの前提の下にヨーイングモーメン
トYMの現在値を取得することによってYM/IZ の現
在値を取得する方式であって、車両の各車輪の車輪すべ
り角α,各車輪のコーナリングパワーC(定数),各車
輪に発生する回転力T等に基づいてヨーイングモーメン
トYMの現在値を取得するものを採用することもできる
のである。
[0008] Here, the definitive "vehicle slip angle estimating means" may employ various things modes of obtaining the current value of YM / I Z. For example, a method of obtaining the current value of YM / I Z by the yaw inertia moment I Z to obtain the current value of the yawing moment YM under the assumption that it is invariant, for each wheel of the vehicle wheel slip It is also possible to adopt a device that acquires the current value of the yawing moment YM based on the angle α, the cornering power C (constant) of each wheel, the rotational force T generated on each wheel, and the like.

【0009】また、「車体横すべり角推定手段」は、横
加速度センサが車体に車両重心点と一致するように取り
付けられる場合には、その横加速度センサにより検出さ
れた横加速度をそのまま車両重心点における横加速度G
y として用いるものとすることができるが、横加速度セ
ンサが車両重心点と一致しない場合には、その横加速度
センサにより検出される横加速度を、その横加速度セン
サと車両重心点との相対位置関係とヨーレートγの時間
微分値とによって補正することによって車両重心点にお
ける横加速度Gy を取得するものとすることができる。
When the lateral acceleration sensor is mounted on the vehicle body so as to coincide with the center of gravity of the vehicle, the "lateral body slip angle estimating means" uses the lateral acceleration detected by the lateral acceleration sensor at the center of gravity of the vehicle as it is. Lateral acceleration G
If the lateral acceleration sensor does not match the center of gravity of the vehicle, the lateral acceleration detected by the lateral acceleration sensor is used as the relative positional relationship between the lateral acceleration sensor and the center of gravity of the vehicle. it can be made to obtain the lateral acceleration G y of the vehicle center of gravity by correcting by the time differential value of the yaw rate γ and.

【0010】[0010]

【作用】本発明に係る車体横すべり角推定装置において
は、例えば上述のようにして車両重心点の回りに発生す
るヨーイングモーメントYMを車両のヨー慣性モーメン
トIZ で割った変数であるYM/IZ の現在値が取得さ
れる。
In the vehicle body slip angle estimating apparatus according to the present invention, for example, YM / I Z is a variable obtained by dividing the yawing moment YM generated around the vehicle center of gravity by the yaw inertia moment I Z of the vehicle as described above. Is obtained.

【0011】車両を急に加,減速したり、ステアリング
ホイールを素早く操作したりする場合の過渡的な車両運
動を除けば、車速Vの変化も車両のローリング運動も無
視することができる。このように車両がローリング運動
なしで水平面内を一定速度で走行する運動すなわち車両
の平面運動にのみ着目することができる場合には、車両
の横すべり運動とヨーイング運動との二つのみを考えれ
ばよいことになり、この場合には一般に、次の二つの運
動方程式が成立する。
Except for the transient vehicle motion when the vehicle is suddenly accelerated or decelerated or when the steering wheel is quickly operated, the change in the vehicle speed V and the rolling motion of the vehicle can be ignored. As described above, when it is possible to pay attention only to the motion in which the vehicle travels in the horizontal plane at a constant speed without the rolling motion, that is, only the plane motion of the vehicle, it is only necessary to consider only two of the sideslip motion and the yawing motion of the vehicle. In this case, the following two equations of motion generally hold.

【0012】 mV(β′+γ)=CF ・・・(1) なる横すべり運動に関する式(1) と、 IZ γ′=YM ・・・(2) なるヨーイング運動に関する式(2) との二つが成立する
のである。ただし、 m:車両の慣性質量 V:車速 β:車両重心点における車体横すべり角(左方向が正、
右方向が負) β′:車体横すべり角βの時間微分値である車体横すべ
り角微分 γ:ヨーレート(反時計方向が正、時計方向が負) γ′:ヨーレートγの時間微分値であるヨーレート微分 CF:車両の前輪にそれのすべり角αf と車両のエンジ
ンおよびブレーキの作動とに基づいて発生する前輪コー
ナリングフォースFf と、後輪にそれのすべり角αr
エンジンおよびブレーキの作動とに基づいて発生する後
輪コーナリングフォースFr との和 IZ :車両のヨー慣性モーメント YM:前輪コーナリングフォースFf と後輪コーナリン
グフォースFr とに基づいて車両重心点の回りに発生す
るヨーイングモーメント
MV (β ′ + γ) = CF (1) Equation (1) relating to the side slipping motion, and I Z γ ′ = YM (2) Equation (2) relating to the yawing motion: One is established. Where: m: inertial mass of the vehicle V: vehicle speed β: body slip angle at the center of gravity of the vehicle (left is positive,
The right direction is negative) β ': The vehicle side slip angle differential which is the time differential value of the vehicle side slip angle β γ: The yaw rate (positive in the counterclockwise direction and negative in the clockwise direction) γ': The yaw rate differential which is the time differential value of the yaw rate γ CF: a front wheel cornering force F f generated based on the operation of its slip angle alpha f and the vehicle engine and braking to the front wheels of the vehicle, to the operation of the engine and brake and its slip angle alpha r on the rear wheel sum I Z of the wheel cornering force F r after generated based: yaw moment of inertia of the vehicle YM: yawing moment generated around the vehicle center of gravity based on the front wheel cornering force F f and the rear wheel cornering force F r

【0013】本発明に係る車体横すべり角推定装置にお
いては、コーナリングフォースCFは直接には検出され
ないが、前記式(1) におけるV(β′+γ)が車両重心
点における横加速度Gy に一致するとの事実、すなわ
ち、 V(β′+γ)=Gy ・・・(3) なる式(3) で表される関係を利用し、間接に検出され
る。
In the vehicle body slip angle estimating apparatus according to the present invention, the cornering force CF is not directly detected, but when V (β '+ γ) in the above equation (1) matches the lateral acceleration G y at the vehicle center of gravity. That is, V (β ′ + γ) = G y (3) The relation is indirectly detected by using the relationship expressed by the equation (3).

【0014】したがって、式(1) は β′=Gy /V−γ ・・・(1) ′ なる式(1) ′に変形でき、また、式(2) は γ′=YM/IZ ・・・(2) ′ なる式(2) ′にそれぞれ変形できる。Therefore, equation (1) can be transformed into equation (1) 'of β' = G y / V-γ (1) ', and equation (2) can be expressed as γ' = YM / I Z .. (2) ′.

【0015】これら式(1) ′および(2) ′における物理
量のうち、横加速度Gy ,車速V,ヨーレートγおよび
YM/IZ は検出可能であるから、本発明に係る車体横
すべり角検出装置においては、それら式(1) ′および
(2) ′と検出可能な物理量とから車体横すべり角βが時
々刻々推定される。それら式(1) ′および(2) ′で記述
される車両の平面運動に関して、検出不能な車体横すべ
り角βを状態変数として含む可観測な線形定係数システ
ム(線形時不変システムともいわれる)が想定され、さ
らに、それに対してオブザーバが想定されることによっ
て車体横すべり角βが推定されるのである。
Since the lateral acceleration G y , the vehicle speed V, the yaw rate γ, and the YM / I Z can be detected among the physical quantities in the equations (1) ′ and (2) ′, the apparatus for detecting the vehicle body slip angle according to the present invention is provided. In formulas (1) ′ and
(2) The vehicle body slip angle β is momentarily estimated from 'and the detectable physical quantity. With respect to the plane motion of the vehicle described by equations (1) 'and (2)', an observable linear constant coefficient system (also called a linear time-invariant system) is assumed, which contains the undetectable body slip angle β as a state variable. Further, the vehicle body slip angle β is estimated by assuming an observer therewith.

【0016】車両の平面運動に関して線形定係数システ
ムを想定するために、例えば次のような入力変数,出力
変数および状態変数を導入することが考えられる。 第1の入力変数u1 (t) :Gy /V 第2の入力変数u2 (t) :YM/IZ 出力変数y(t) :γ 第1の状態変数x1 (t) :β 第2の状態変数x2 (t) :γ しかし、それら入力変数,出力変数および状態変数で
は、可観測な線形定係数システムを構成することができ
ない。以下、この理由を説明する。
In order to assume a linear constant coefficient system for the plane motion of the vehicle, it is conceivable to introduce, for example, the following input variables, output variables and state variables. First input variable u 1 (t): G y / V Second input variable u 2 (t): YM / I Z Output variable y (t): γ First state variable x 1 (t): β Second state variable x 2 (t): γ However, these input variables, output variables and state variables cannot form an observable linear constant coefficient system. Hereinafter, the reason will be described.

【0017】このような入力変数,出力変数および状態
変数を導入する場合には、次のような行列形式で表され
る状態方程式(4) が誘導される。
When such input variables, output variables and state variables are introduced, a state equation (4) expressed in the following matrix form is derived.

【0018】[0018]

【数1】 (Equation 1)

【0019】さらに、次のような行列形式で表される出
力方程式(5) も誘導される。
Further, an output equation (5) expressed in the following matrix form is derived.

【0020】[0020]

【数2】 (Equation 2)

【0021】一方、線形定係数システムが可観測である
ためには、可観測行列U0 のランクが状態変数の次元n
(今回は2)に一致することが必要である。可観測行列
0 は、後述のシステム行列Aと出力行列Cとを用いて
表され、次元nの値が2である場合には、
On the other hand, in order for the linear constant coefficient system to be observable, the rank of the observable matrix U 0 must be equal to the dimension n of the state variable.
(This time, 2) must match. The observable matrix U 0 is represented using a system matrix A and an output matrix C, which will be described later, and when the value of the dimension n is 2,

【0022】[0022]

【数3】 (Equation 3)

【0023】で表される。上述の状態方程式(4) および
出力方程式(5) においては、
## EQU2 ## In the above state equation (4) and output equation (5),

【0024】[0024]

【数4】 (Equation 4)

【0025】C=[0 1] であるから、結局、Since C = [0 1], after all,

【0026】[0026]

【数5】 (Equation 5)

【0027】となり、これのランクは1であって、次元
nである2ではないから、今回の線形定係数システムは
可観測ではない。
Since the rank is 1 and is not 2 which is the dimension n, the present linear constant coefficient system is not observable.

【0028】そこで、本発明に係る車体横すべり角推定
装置においては、車体横すべり角βは時間の経過につれ
て比較的緩やかに変化する特性を持っていて、車体横す
べり角βの今回推定値β(i) と前回推定値β(i-1) との
差がほぼ0とみなし得るという事実を利用し、新たに、
次のような入力変数,出力変数および状態変数を導入す
る。 第1の入力変数u1 (t) :Gy /V 第2の入力変数u2 (t) :YM/IZ −β(i-1) 出力変数y(t) :γ 第1の状態変数x1 (t) :β(i) 第2の状態変数x2 (t) :γ したがって、次のような状態方程式(6) が誘導される。
Therefore, in the vehicle body slip angle estimating apparatus according to the present invention, the vehicle body slip angle β has a characteristic that changes relatively slowly with time, and the current estimated value β (i) of the vehicle body slip angle β Using the fact that the difference between the estimated value and the previous estimated value β (i-1) can be regarded as almost 0,
Introduce the following input, output, and state variables. First input variable u 1 (t): G y / V Second input variable u 2 (t): YM / I Z(i-1) Output variable y (t): γ First state variable x 1 (t): β (i) Second state variable x 2 (t): γ Accordingly, the following state equation (6) is derived.

【0029】[0029]

【数6】 (Equation 6)

【0030】さらに、次のような出力方程式(7) も誘導
される。
Further, the following output equation (7) is derived.

【0031】[0031]

【数7】 (Equation 7)

【0032】式(6) においては、システム行列AがIn equation (6), the system matrix A is

【0033】[0033]

【数8】 (Equation 8)

【0034】で表されるため、今回の可観測行列U0 Thus, the observable matrix U 0 of this time is

【0035】[0035]

【数9】 (Equation 9)

【0036】で表されることになる。この可観測行列U
0 のランクは2であって、状態変数の次元nである2に
等しいため、それら式(6) および(7) で記述される線形
定係数システムは可観測である。
## EQU2 ## This observable matrix U
Since the rank of 0 is 2 and is equal to 2 which is the dimension n of the state variable, the linear constant coefficient systems described by the equations (6) and (7) are observable.

【0037】以上要するに、本発明に係る車体横すべり
角推定装置においては、車両の横すべり運動とヨーイン
グ運動とから成る車両の平面運動に関して、車体横すべ
り角βとヨーレートγとをそれぞれ状態変数とし、車両
重心点における横加速度Gy を車速Vで割った変数であ
るGy /Vと、YM/IZ から車体横すべり角βの前回
推定値を差し引いた変数とをそれぞれ入力変数とし、ヨ
ーレートγを出力変数とする可観測な線形定係数システ
ムが想定されているのであり、その線形定係数システム
においてオブザーバが用いられて、横加速度Gy ,車速
Vおよびヨーレートγの現在値と、取得したYM/IZ
の現在値と、車体横すべり角βの前回推定値β(i-1)
に基づいて車体横すべり角βの今回値β(i) が推定され
るのである。
In short, in the vehicle body slip angle estimating apparatus according to the present invention, the vehicle body slip angle β and the yaw rate γ are set as state variables for the vehicle plane movement including the vehicle skid movement and the yawing movement, and the vehicle center of gravity is set. the lateral acceleration G y is the G y / V is a variable divided by the vehicle speed V, the from YM / I Z respectively input variables and variables obtained by subtracting the previous estimated value of the vehicle body slip angle β at the point, the output variable yaw rate γ An observable linear constant coefficient system is assumed, and an observer is used in the linear constant coefficient system to obtain the current values of the lateral acceleration G y , the vehicle speed V and the yaw rate γ, and the obtained YM / I Z
And the current value of is the previous estimate of the vehicle body slip angle β β (i-1) between the present value of the vehicle slip angle beta, based on the beta (i) is estimated.

【0038】オブザーバは、横加速度Gy およびヨーレ
ートγを直接積分することなく車体横すべり角βの今回
値β(i) を推定するものである。そのため、前述のよう
に、それら横加速度Gy 等にドリフトが発生しても、そ
れに基づく誤差が積分されずに済み、車体横すべり角β
の推定精度がそのドリフトによって低下せずに済む。
The observer estimates the current value β (i) of the vehicle body slip angle β without directly integrating the lateral acceleration G y and the yaw rate γ. Therefore, as described above, even if a drift occurs in the lateral acceleration Gy or the like, an error based on the drift is not integrated, and the vehicle body side slip angle β
Does not need to be reduced due to the drift.

【0039】ところで、各種センサによる検出値は普通
ノイズを含んでいるため、このような検出値を積分する
と車体横すべり角βを精度よく推定することができな
い。これに対し、本発明に係る車体横すべり角推定装置
においては、オブザーバを用いることにより、そのよう
な検出値を積分することなく車体横すべり角βが推定さ
れるから、推定された車体横すべり角βが検出値のノイ
ズの影響をそれほど受けずに済む。
By the way, since values detected by various sensors usually include noise, it is impossible to accurately estimate the vehicle body slip angle β by integrating such detected values. On the other hand, in the vehicle body slip angle estimating device according to the present invention, the use of the observer allows the vehicle body slip angle β to be estimated without integrating such a detected value. The detection value is less affected by noise.

【0040】[0040]

【発明の効果】このように、本発明によれば、各種セン
サによる検出値を積分することなく車体横すべり角βを
推定することができるから、推定された車体横すべり角
βが検出値のノイズ,ドリフト等の影響をそれほど受け
ずに済むという効果が得られる。
As described above, according to the present invention, the vehicle body slip angle β can be estimated without integrating values detected by various sensors. The effect is obtained that the influence of drift or the like is not so much affected.

【0041】[0041]

【実施例】以下、本発明の一実施例である車体横すべり
角推定装置を含む後輪舵角制御装置を図面に基づいて詳
細に説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS A rear wheel steering angle control device including a vehicle body slip angle estimation device according to an embodiment of the present invention will be described in detail with reference to the drawings.

【0042】この後輪舵角制御装置は、エンジンの出力
をオートマチックトランスミッションを経て左右前輪と
左右後輪とにそれぞれ伝達する4輪駆動式自動車である
車両に設けられており、図2に示すように、ステアリン
グホイール10の操舵角に応じて左右の前輪12の舵角
δf を変化させる前輪操舵機構16と、駆動源としての
ステップモータ20により左右の後輪22の舵角δr
変化させる後輪操舵機構26とを備えている。
This rear wheel steering angle control device is provided in a vehicle which is a four-wheel drive type vehicle that transmits the output of the engine to the left and right front wheels and the left and right rear wheels via an automatic transmission, as shown in FIG. to a front wheel steering mechanism 16 for changing the steering angle [delta] f of the left and right front wheels 12 in accordance with the steering angle of the steering wheel 10, to vary the steering angle [delta] r of the left and right rear wheels 22 by a step motor 20 as a driving source And a rear wheel steering mechanism 26.

【0043】この後輪舵角制御装置は、さらに、同図に
示すように、前輪舵角δf を検出する前輪舵角センサ3
0,後輪舵角δr を検出する後輪舵角センサ32,車体
の走行速度である車速Vを検出する車速センサ34,ヨ
ーレートγを検出するヨーレートセンサ36,および車
両重心点における横加速度Gy を検出する横加速度セン
サ40を備えている。それらセンサ30,32,34,
36,40はコントローラ50の入力部に接続され、そ
れの出力部には前記ステップモータ20が接続されてい
る。
[0043] The rear wheel steering angle control apparatus further includes, as shown in the figure, the front wheel steering angle sensor 3 for detecting a front wheel steering angle [delta] f
0, wheel steering angle sensor 32 after detecting the rear wheel steering angle [delta] r, the vehicle speed sensor 34 for detecting the vehicle speed V is a vehicle speed, lateral acceleration in the yaw rate sensor 36, and the vehicle center of gravity detecting a yaw rate gamma G A lateral acceleration sensor 40 for detecting y is provided. These sensors 30, 32, 34,
Reference numerals 36 and 40 are connected to an input section of a controller 50, and the output section thereof is connected to the step motor 20.

【0044】前輪舵角センサ30および後輪舵角センサ
32はそれぞれ、反時計方向の舵角δf ,δr を正、時
計方向の舵角δf ,δr を負として検出するものであ
る。ヨーレートセンサ36は、反時計方向のヨーレート
γを正、時計方向のヨーレートγを負として検出するも
のである。横加速度センサ40は、車両の重心点に配置
されていて、車両重心点における横加速度Gy を左方向
を負、右方向を正として検出するものである。
[0044] and detects each front wheel steering angle sensor 30 and the rear wheel steering angle sensor 32, the counterclockwise steering angle [delta] f, the [delta] r positive, clockwise steering angle [delta] f, as a negative the [delta] r . The yaw rate sensor 36 detects the counterclockwise yaw rate γ as positive and the clockwise yaw rate γ as negative. The lateral acceleration sensor 40, located on the center-of-gravity point of the vehicle, and detects the lateral acceleration G y of the vehicle center of gravity negative leftward, rightward as positive.

【0045】コントローラ50は、図示しないCPU,
ROM,RAMおよびバスを含むコンピュータを主体と
して構成されており、そのROMに、図1のフローチャ
ートで表される車体横すべり角推定ルーチンおよび図3
のフローチャートで表される後輪舵角制御ルーチンを始
めとする各種プログラム,図4,図5および図6のグラ
フで表される目標後輪舵角決定のためのマップを始めと
する各種マップ等が予めストアされている。
The controller 50 includes a CPU (not shown),
The computer mainly includes a ROM, a RAM, and a bus. The ROM includes a vehicle body slip angle estimation routine shown in the flowchart of FIG.
Various programs including a rear wheel steering angle control routine represented by the flowchart of FIG. 4, various maps including a map for determining a target rear wheel steering angle represented by the graphs of FIGS. 4, 5, and 6, and the like. Is stored in advance.

【0046】図1の車体横すべり角推定ルーチンは、次
の式(8) に従って車両重心点の回りに発生するヨーイン
グモーメントYMを推定することが可能とされている。 YM=af f −ar r ・・・(8) ただし、 af :前輪車軸と車両重心点との距離(定数) ar :後輪車軸と車両重心点との距離(定数) Ff :左右前輪12に発生する合計コーナリングフォー
ス(以下、単に前輪コーナリングフォースという) Fr :左右後輪22に発生する合計コーナリングフォー
ス(以下、単に後輪コーナリングフォースという) なお、前輪コーナリングフォースFf は、 −αf f (1−hf f ) で定義され、一方、後輪コーナリングフォースFr は、 −αr r (1−hr r ) で定義される。ただし、 Cf :左右前輪12に係る合計コーナリングパワー(以
下、単に前輪コーナリングパワーという。定数) Cr :左右後輪22に係る合計コーナリングパワー(以
下、単に後輪コーナリングパワーという。定数) αf :左右前輪12の車輪すべり角 αr :左右後輪22の車輪すべり角 Tf :エンジンおよび前輪ブレーキの作動に基づいて左
右前輪12に発生する合計回転力(以下、単に前輪回転
力という) Tr :エンジンおよび後輪ブレーキの作動に基づいて左
右後輪22に発生する合計回転力(以下、単に後輪回転
力という) hf :前輪回転力Tf による前輪コーナリングパワーC
f の減少係数(定数) hr :後輪回転力Tr による後輪コーナリングパワーC
r の減少係数(定数)
The vehicle body slip angle estimation routine shown in FIG. 1 is capable of estimating the yawing moment YM generated around the vehicle center of gravity according to the following equation (8). YM = a f F f -a r F r ··· (8) However, a f: front wheel axle and the distance between the vehicle center of gravity (constant) a r: rear wheel axle and the distance between the vehicle center of gravity (constant) F f : Total cornering force generated on left and right front wheels 12 (hereinafter simply referred to as front wheel cornering force) F r : Total cornering force generated on right and left rear wheels 22 (hereinafter simply referred to as rear wheel cornering force) Note that front wheel cornering force F f is defined by −α f C f (1−h f T f ), while the rear wheel cornering force F r is defined by −α r C r (1−h r T r ). Here, C f : the total cornering power of the left and right front wheels 12 (hereinafter simply referred to as front wheel cornering power; constant) C r : the total cornering power of the left and right rear wheels 22 (hereinafter simply referred to as rear wheel cornering power. Constant) α f : Wheel slip angle of the left and right front wheels 12 α r : Wheel slip angle of the left and right rear wheels 22 T f : Total torque generated in the left and right front wheels 12 based on the operation of the engine and the front wheel brake (hereinafter simply referred to as front wheel torque) T r : total rotational force generated on the left and right rear wheels 22 based on the operation of the engine and the rear wheel brake (hereinafter, simply referred to as rear wheel rotational force) hf : front wheel cornering power C by front wheel rotational force Tf
reduction factor f (constant) h r: rear wheel rotational force T r rear wheel cornering power by C
Reduction factor of r (constant)

【0047】また、前輪すべり角αf は、 β+af γ/V−δf で定義され、一方、後輪すべり角αr は、 β−ar γ/V−δr で定義される。The front wheel slip angle α f is defined by β + a f γ / V−δ f , while the rear wheel slip angle α r is defined by β−a r γ / V−δ r .

【0048】また、減少係数hf は、 δf /(αf f ) で定義され、一方、減少係数hr は、 δr /(αr r ) で定義される。The reduction coefficient h f is defined by δ f / (α f C f ), while the reduction coefficient h r is defined by δ r / (α r C r ).

【0049】本実施例においては、車両の平面運動が線
形領域にあるとして前輪コーナリングパワーCf も後輪
コーナリングパワーCr も減少係数hf ,hr もそれぞ
れ定数とされている。しかし、前輪,後輪すべり角
αf ,αr はそれぞれ、上述の各定義式を用いた演算に
よって取得される。具体的には、前輪すべり角αf につ
いては、それの定義式、すなわち、 β+af γ/V−δf なる式において、βが後述の前回推定値β(i-1) とさ
れ、γ,Vおよびδf がそれぞれ各種センサによる検出
値とされることによって取得され、一方、後輪すべり角
αr については、それの定義式、すなわち、 β−ar γ/V−δr なる式において、βが後述の前回推定値β(i-1) とさ
れ、γ,Vおよびδr がそれぞれ各種センサによる検出
値とされることによって取得される。また、前輪,後輪
回転力Tf ,Tr もそれぞれ、次のようにして演算によ
って取得される。
In the present embodiment, there is a respective constant wheel cornering power C f also the rear wheel cornering power C r also decreases the coefficient h f, h r as the plane motion of the vehicle is in the linear region. However, the front wheel and rear wheel slip angles α f , α r are obtained by the calculations using the above-described respective definition expressions. Specifically, for the front wheel slip angle α f , β is defined as a previously estimated value β (i−1) described later in a definition equation thereof, that is, an equation of β + a f γ / V−δ f , and γ, V and [delta] f is obtained by being a value detected by various sensors, respectively, whereas, for the rear wheel slip angle alpha r, its definition formula, i.e., the β-a r γ / V- δ r becomes formula , beta is the previous estimate of below β (i-1), γ , V and [delta] r is obtained by being a value detected by various sensors, respectively. The front wheel and rear wheel rotational forces T f and Tr are also obtained by calculation as follows.

【0050】すなわち、前輪回転力Tf は、車両のオー
トマチックトランスミッション(以下、単にトランスミ
ッションという)から前輪12に伝達される前輪駆動力
から車両の前輪ブレーキによって前輪12に作用させら
れる制動力を差し引いたものとして取得され、前輪駆動
力は、エンジンのスロットルバルブの開度(以下、単に
スロットル開度という)θ,エンジンの回転数Ne ,ト
ランスミッションのギヤ比R,そのトランスミッション
からの出力が前輪12に配分される比率等に基づいて取
得され、一方、前輪制動力は、前輪ブレーキのホイール
シリンダ圧等に基づいて取得される。後輪回転力Tr
それと同様にして取得される。具体的には、トランスミ
ッションから後輪22に伝達される後輪駆動力から車両
の後輪ブレーキによって後輪22に作用させられる後輪
制動力を差し引いたものとして取得され、後輪駆動力
は、スロットル開度θ,エンジンの回転数Ne ,トラン
スミッションのギヤ比R,そのトランスミッションから
の出力が後輪22に配分される比率等に基づいて取得さ
れ、一方、後輪制動力は、後輪ブレーキのホイールシリ
ンダ圧等に基づいて取得される。
That is, the front wheel rotational force Tf is obtained by subtracting the braking force applied to the front wheels 12 by the front wheel brake of the vehicle from the front wheel driving force transmitted from the automatic transmission (hereinafter simply referred to as transmission) of the vehicle to the front wheels 12. The driving force of the front wheels is obtained as follows: the opening degree of the throttle valve of the engine (hereinafter simply referred to as the throttle opening degree) θ, the engine speed Ne , the transmission gear ratio R, and the output from the transmission to the front wheels 12. The front wheel braking force is obtained based on the allocated ratio and the like, while the front wheel braking force is obtained based on the wheel cylinder pressure of the front wheel brake and the like. The rear wheel rotational force Tr is obtained in the same manner. Specifically, the rear wheel driving force transmitted from the transmission to the rear wheels 22 is obtained by subtracting the rear wheel braking force applied to the rear wheels 22 by the rear wheel brake of the vehicle. The throttle opening θ, the engine speed N e , the transmission gear ratio R, and the output from the transmission are distributed based on the ratio of the rear wheels 22 and the like. Is obtained based on the wheel cylinder pressure of the vehicle.

【0051】したがって、前輪コーナリングフォースF
f も後輪コーナリングフォースFr も確定し、前述の、
車両のヨーイング運動に関する式(2) は、結局、 γ′=(af f −ar r )/IZ ・・・(9) なる式(9) に変形されることになる。この式においてヨ
ー慣性モーメントIZ は定数扱いされるから、結局、
(af f −ar r )/IZ すなわちヨーレート微分
γ′の現在値は検出可能な物理量であるということがで
きる。
Therefore, the front wheel cornering force F
f and the rear wheel cornering force Fr are determined,
Expression for the yawing motion of the vehicle (2) is, after all, γ '= (a f F f -a r F r) / I becomes Z · · · (9) made to be deformed into the equation (9). In this equation, since the yaw moment of inertia I Z is treated as a constant, after all,
(A f F f -a r F r) / it current value of Z i.e. the yaw rate differentiation gamma 'can be said to be detectable physical quantity.

【0052】さらに、この車体横すべり角推定ルーチン
は、車体横すべり角βを逐次推定し、具体的には、車両
の横すべり運動とヨーイング運動とから成る車両の平面
運動に関して想定される可観測な線形定係数システムに
おいてオブザーバを想定して、横加速度Gy ,車速Vお
よびヨーレートγの現在値と、ヨーイングモーメントY
Mを車両のヨー慣性モーメントIZ で割ることによって
取得したYM/IZ の現在値(=(af f −a
r r )/IZ )と、車体横すべり角βの前回推定値β
(i-1) とに基づいて車体横すべり角βの今回値β(i)
推定することも可能とされている。
Further, the vehicle body slip angle estimation routine successively estimates the vehicle body slip angle β, and specifically, an observable linear constant assumed for a plane motion of the vehicle including a vehicle slip motion and a yawing motion. Assuming an observer in the coefficient system, the current values of the lateral acceleration G y , the vehicle speed V and the yaw rate γ, and the yawing moment Y
The current value of YM / I Z obtained by dividing M by the yaw moment of inertia I Z of the vehicle (= (a f F f −a
r F r ) / I Z ) and the previous estimated value β of the vehicle body slip angle β
It is also possible to estimate the current value β (i) of the vehicle body side slip angle β based on (i-1) .

【0053】その可観測な線形定係数システムは、具体
的には、次のような入力変数,出力変数および状態変数
を導入されており、 第1の入力変数u1 (t) :Gy /V 第2の入力変数u2 (t) :(af f −ar r )/I
Z −β(i-1) 出力変数y(t) :γ 第1の状態変数x1 (t) :β(i) 第2の状態変数x2 (t) :γ
In the observable linear constant coefficient system, specifically, the following input variables, output variables and state variables are introduced, and a first input variable u 1 (t): G y / V second input variable u 2 (t): (a f F f -a r F r) / I
Z- β (i-1) Output variable y (t): γ First state variable x 1 (t): β (i) Second state variable x 2 (t): γ

【0054】[0054]

【数10】 (Equation 10)

【0055】なる状態方程式(10)と、The state equation (10) becomes:

【0056】[0056]

【数11】 [Equation 11]

【0057】なる出力方程式(11)とで記述される。This is described by the following output equation (11).

【0058】この線形定係数システムは、上述の説明か
ら明らかなように、2次元の状態変数を持っているが、
その一方であるγについては検出可能であり、他方であ
るβ(i) のみが検出不能であるから、今回は、二つの状
態変数すべてを推定するのではなく車体横すべり角βの
みを推定し、オブザーバの次元を減らして最小次元オブ
ザーバを構成することとする。
As is clear from the above description, this linear constant coefficient system has a two-dimensional state variable.
On the other hand, one γ can be detected, and only the other β (i) cannot be detected.This time, instead of estimating all two state variables, only the vehicle body slip angle β is estimated, The dimension of the observer is reduced to form the minimum dimension observer.

【0059】ところで、線形定係数システムは一般に、 X(t) ′=AX(t) +BU(t) なる状態方程式と、 Y(t) =CX(t) なる出力方程式とで記述される。ただし、 X(t) :n次元の状態ベクトル Y(t) :r次元の出力ベクトル U(t) :m次元の入力ベクトル A:n×n次元のシステム行列 B:n×m次元の入力行列 C:r×n次元の出力行列Incidentally, the linear constant coefficient system is generally described by a state equation of X (t) '= AX (t) + BU (t) and an output equation of Y (t) = CX (t). X (t): n-dimensional state vector Y (t): r-dimensional output vector U (t): m-dimensional input vector A: n × n-dimensional system matrix B: n × m-dimensional input matrix C: r × n-dimensional output matrix

【0060】そして、このシステムについては一般に、 ω(t) ′=A′ω(t) +KY(t) +B′U(t) X(t) =Dω(t) +HY(t) なる式で表される最小次元オブザーバを構成することが
できる。ただし、 ω(t) :X(t) から変換され、かつX(t) より次元が低
い変数 ω(t) ′:ω(t) の時間微分値 A′:(n−r)×(n−r)次元の行列 K:(n−r)×r次元の行列 B′:(n−r)×m次元の行列 D:n×(n−r)次元の行列 H:n×r次元の行列 以下、最小次元オブザーバの一般的な構成法について説
明するが、これは、文献『システム制御理論入門 (実
教出版株式会社 平成3年1月20日発行)』に詳細に
記載されているため、ここでは簡単に説明する。
Then, this system is generally expressed by the following equation: ω (t) ′ = A′ω (t) + KY (t) + B′U (t) X (t) = Dω (t) + HY (t) Can be constructed. Where ω (t): a variable converted from X (t) and lower in dimension than X (t) ω (t) ′: time derivative of ω (t) A ′: (n−r) × (n -R) dimensional matrix K: (nr) × r dimensional matrix B ′: (nr) × m dimensional matrix D: nx (nr) dimensional matrix H: nxr dimensional matrix Matrix Hereinafter, a general construction method of the minimum-dimensional observer will be described. This is described in detail in the document "Introduction to System Control Theory (Jikkyo Shuppan Co., Ltd., published on January 20, 1991)". Here, a brief description will be given.

【0061】各行列内の各要素を決定するためにまず、
n×n次元の行列Sを決定する。この行列Sは
To determine each element in each matrix, first
An n × n-dimensional matrix S is determined. This matrix S is

【0062】[0062]

【数12】 (Equation 12)

【0063】で表され、ここにおいてCは前述の出力行
列であるが、Wは行列Sの行列式を0にしない適当な
(n−r)×n次元の行列である。次に、
Where C is the output matrix described above, and W is an appropriate (n−r) × n-dimensional matrix that does not set the determinant of the matrix S to zero. next,

【0064】[0064]

【数13】 (Equation 13)

【0065】における行列A11,A12,A21およびA22
の各々を決定し、さらに、
The matrices A 11 , A 12 , A 21 and A 22 in
Is determined, and further,

【0066】[0066]

【数14】 [Equation 14]

【0067】における行列B1 およびB2 の各々を決定
する。続いて、(n−r)×r次元の行列Lを設計パラ
メータとして導入し、さらに、行列A′を、 A′=A22−LA12 なる式で定義する。そして、この行列A′の固有値がオ
ブザーバの極(任意の値)となるように行列Lの各要素
を決定する。その後、行列K,B′,DおよびH内の各
要素を K=A′L+A21−LA11 B′=−LB1 +B2
Determine each of the matrices B 1 and B 2 in Subsequently, an (n−r) × r-dimensional matrix L is introduced as a design parameter, and a matrix A ′ is defined by the following equation: A ′ = A 22 −LA 12 . Then, each element of the matrix L is determined such that the eigenvalue of the matrix A 'becomes the pole (arbitrary value) of the observer. Thereafter, the matrix K, B ', each element in the D and H K = A'L + A 21 -LA 11 B' = - LB 1 + B 2

【0068】[0068]

【数15】 (Equation 15)

【0069】[0069]

【数16】 (Equation 16)

【0070】なる式が成立するように決定する。The following equation is established.

【0071】以上、最小次元オブザーバの一般的な構成
法を簡単に説明したが、本実施例における構成法を具体
的に説明する。前述のように、本実施例においては、
The general configuration of the minimum-dimensional observer has been briefly described above, but the configuration in the present embodiment will be specifically described. As described above, in this embodiment,

【0072】[0072]

【数17】 [Equation 17]

【0073】[0073]

【数18】 (Equation 18)

【0074】C=[0 1] である。このとき、 W=[1 0] とすれば、行列Sの行列式が0でないという要件を満た
すことができる。また、今回は行列Lは1個の要素から
成る行列になるため、行列A′の固有値は−Lで表さ
れ、今回はそれが今回の最小次元オブザーバの極である
−pに一致すると仮定すれば、 L=p となる。なお、このオブザーバがその成立条件を満たす
ためには、行列A′が安定行列であることが必要である
ため、オブザーバの極−pの値は負とされ、結局、その
pの値が正とされることになる。
C = [0 1]. At this time, if W = [1 0], the requirement that the determinant of the matrix S is not 0 can be satisfied. In this case, since the matrix L is a matrix composed of one element, the eigenvalue of the matrix A 'is represented by -L, and this time, it is assumed that it coincides with -p which is the pole of the minimum dimension observer this time. Then, L = p. In order for the observer to satisfy the condition, the matrix A 'needs to be a stable matrix. Therefore, the value of the pole -p of the observer is negative, and the value of p is eventually positive. Will be done.

【0075】したがって、前記行列A′,K,B′,D
およびHは今回は、以下のようになる。 A′=−p K=−p2 −1 B′=[1 −p]
Therefore, the matrices A ', K, B', D
And H this time are as follows. A ′ = − p K = −p 2 −1 B ′ = [1−p]

【0076】[0076]

【数19】 [Equation 19]

【0077】[0077]

【数20】 (Equation 20)

【0078】今回は前記ωは1次元であって、第1の状
態変数x1 (t) である車体横すべり角βが変換されたも
のとなるから、今回の最小次元オブザーバの式は、
In this case, ω is one-dimensional, and the vehicle state slip angle β, which is the first state variable x 1 (t), is converted.

【0079】[0079]

【数21】 (Equation 21)

【0080】なる式(12)と、Equation (12)

【0081】[0081]

【数22】 (Equation 22)

【0082】すなわち、 β(i) =ω+pγ ・・・(13) なる式(13)との二つとなる。したがって、まず、先の式
(12)に従い、既知または検出可能な物理量であるp,
γ,Gy /V,(af f −ar r )/IZ およびβ
(i-1) からωを推定し、その後、式(13)に従い、そのω
と、既知または検出可能な物理量であるpおよびγから
β(i) 、すなわち、車体横すべり角βの今回値を推定す
る。ただし、式(12)は連続時間システムを記述するた
め、これを離散時間化することによってωを逐次推定す
ることにする。
That is, β (i) = ω + pγ (13) Therefore, first, the previous expression
According to (12), a known or detectable physical quantity p,
γ, G y / V, ( a f F f -a r F r) / I Z and β
ω is estimated from (i-1) , and then, according to equation (13),
And β (i) , that is, the current value of the vehicle body slip angle β, is estimated from the known or detectable physical quantities p and γ. However, since equation (12) describes a continuous-time system, ω is successively estimated by making this a discrete-time system.

【0083】ところで、前述のように、線形定係数シス
テムは一般に、 X(t) ′=EX(t) +FU(t) ・・・(14) なる状態方程式(14)で表され、これは一般に、U(t) の
サンプリング周期をdtで表すことにすると、 X(t) (i) =E′X(t) (i-1) +F′U(t) (i) ・・・(15) なる式(15)に変換されるという規則が既に知られてい
る。なお、この規則についても前記文献『システム制御
理論入門』に記載されているため、詳細な説明は省略す
る。ただし、 X(t) :状態ベクトル X(t) ′:状態ベクトルの時間微分値 U(t) :入力ベクトル E:システム行列 F:入力行列 X(t) (i) :X(t) の今回推定値 X(t) (i-1) :X(t) の前回推定値 U(t) (i) :U(t) の今回検出値 E′:eEdt F′:(eEdt −I)E-1F I:単位行列
As described above, the linear constant coefficient system is generally expressed by a state equation (14) as follows: X (t) '= EX (t) + FU (t) (14) , U (t) is represented by dt. X (t) (i) = E'X (t) (i-1) + F'U (t) (i) (15) The rule that the expression is converted into Expression (15) is already known. Since this rule is also described in the above-mentioned document "Introduction to System Control Theory", detailed description is omitted. X (t): state vector X (t) ': time derivative of state vector U (t): input vector E: system matrix F: input matrix X (t) (i) : this time of X (t) Estimated value X (t) (i-1) : Previous estimated value of X (t) U (t) (i) : Current detected value of U (t) E ': e Edt F': (e Edt -I) E -1 FI: unit matrix

【0084】この規則に従って式(12)を変換するのであ
るが、まず、式(14)の構成を式(12)に揃えると、 X(t) ′=EX(t) +F1 1(t)+F2 2(t)+F3 3(t) ・・・(16) なる式(16)となり、その結果、この式における各要素と
式(12)における各要素との対応関係は、 E=−p X(t) =ω(t) F1 =−p2 −1 F2 =1 F3 =−p U1(t)=γ U2(t)=Gy /V U3(t)=(af f −ar r )/IZ −β(i-1) となる。また、 F1 ′=(eEdt −I)E-11 2 ′=(eEdt −I)E-12 3 ′=(eEdt −I)E-13 であり、また、 E-1=1/(−p) であるから、結局、式(12)は、 ω(i) =e1 ω(i-1) +e2 γ(i) +e3 y (i) /V(i) +e4 (((af f −ar r )/IZ (i) −β(i-1) ) ・・・(17) なる式(17)に変換される。
Equation (12) is converted in accordance with this rule. First, when the construction of equation (14) is adjusted to equation (12), X (t) '= EX (t) + F 1 U 1 (t ) + F 2 U 2 (t) + F 3 U 3 (t) (16) As a result, the relationship between each element in this equation and each element in equation (12) is E = -p X (t) = ω (t) F 1 = -p 2 -1 F 2 = 1 F 3 = -p U 1 (t) = γ U 2 (t) = G y / V U 3 ( t) = become (a f F f -a r F r) / I Z -β (i-1). F 1 ′ = (e Edt −I) E −1 F 1 F 2 ′ = (e Edt −I) E −1 F 2 F 3 ′ = (e Edt −I) E −1 F 3 , In addition, since E −1 = 1 / (− p), the equation (12) eventually becomes: ω (i) = e 1 ω (i-1) + e 2 γ (i) + e 3 G y (i) is converted to / V (i) + e 4 (((a f F f -a r F r) / I Z) (i) -β (i-1)) ··· (17) becomes equation (17) .

【0085】ただし、 ω(i) :ωの今回推定値 ω(i-1) :ωの前回推定値 γ(i) :ヨーレートγの今回サンプリング値 Gy (i) :横加速度Gy の今回サンプリング値 V(i) :車速Vの今回サンプリング値 ((af f −ar r )/IZ (i) :((af f
−ar r )/IZ ) の今回演算値
Where ω (i) : the current estimated value of ω ω (i-1) : the previous estimated value of ω γ (i) : the current sampling value of the yaw rate γ G y (i) : the current value of the lateral acceleration G y sampled values V (i): this sampling value of the vehicle speed V ((a f F f -a r F r) / I Z) (i): ((a f F f
-A r F r ) / I z )

【0086】式(17)においてe1 ,e2 ,e3 およびe
4 はいずれも定数であって、具体的には、 e1 :e-pdt2 :(e-pdt−1)(p2 +1)/p e3 :−(e-pdt−1)/p e4 :e-pdt−1 となる。
In the equation (17), e 1 , e 2 , e 3 and e
4 is a both a constant, specifically, e 1: e -pdt e 2 : (e -pdt -1) (p 2 +1) / p e 3 :-( e -pdt -1) / p e 4 : e −pdt −1

【0087】したがって、式(17)にω(i-1) ,γ(i)
y (i) ,V(i) および((af f −ar r )/I
Z (i) の各値を代入すればω(i) の値が取得でき、こ
れとγ(i) の値とを、前記式(13)に相当する式であっ
て、 β(i) =ω+pγ なる式に代入すればβ(i) の値が取得できる。なお、式
(17)では、V(i) =0であるとω(i) の演算が不能とな
るため、このときはω(i) を0とするとともに、β(i)
も0とすることにする。
Therefore, the equation (17)(i-1), Γ(i),
Gy (i), V(i)And ((afF f-ArFr) / I
Z)(i)Substituting each value of ω(i)Can be obtained.
And γ(i)Is a formula corresponding to the formula (13).
And β(i)= Ω + pγ, then β(i)Can be obtained. Note that the expression
In (17), V(i)= 0 if = 0(i)Operation becomes impossible
Therefore, at this time ω(i)To 0 and β(i)
Is also set to 0.

【0088】以上詳記した事情を背景とし、図1の車体
横すべり角推定ルーチンは、サンプリング周期dtが経
過するごとに実行される。本ルーチンの各回の実行時に
は、まず、ステップS1(以下、単にS1で表す。他の
ステップについても同じとする)において、本ルーチン
の今回の実行がコントローラ50の電源が投入された後
の初回であるか否かが判定される。今回はそうであるか
ら判定がYESとなり、S2においてωの値が0、βの
値がpγとされる。pは予定された定数であり、γは現
在のヨーレートである。本ステップは、ωおよびβのそ
れぞれの初期値を設定するためのステップなのである。
Against the background described above, the vehicle body slip angle estimation routine of FIG. 1 is executed every time the sampling period dt elapses. In each execution of this routine, first, in step S1 (hereinafter simply referred to as S1; the same applies to other steps), this execution of this routine is performed for the first time after the power of the controller 50 is turned on. It is determined whether there is. Since this is the case this time, the determination is YES, and the value of ω is set to 0 and the value of β is set to pγ in S2. p is a predetermined constant and γ is the current yaw rate. This step is a step for setting respective initial values of ω and β.

【0089】続いて、S3において、車速Vが0である
か否かが判定される。今回はそうではないと仮定すれ
ば、判定がNOとなり、S4において、まず、前輪,後
輪コーナリングフォースFf ,Fr の現在値が前記式に
従って演算される。本ステップにおいては、さらに、前
記式(17)に相当する式(図1に示す式であって、以下、
演算式という)に、ω,γ,Gy ,Vおよびβのそれぞ
れの現在値(ただし、今回は、ωおよびβについては上
述の初期値に等しい)とFf およびFr のそれぞれの今
回演算値とが代入されることによってωの現在値が更新
される。なお、この式におけるe1 ,e2 ,e3
4 ,af ,ar およびIZ はそれぞれ定数である。そ
の後、S5において、ωの現在値とpγの現在値との和
が演算されることによってβの現在値が更新され、S6
において、そのβの現在値が、車体横すべり角βの今回
推定値としてRAMにストアされる。以上で本ルーチン
の一回の実行が終了する。
Subsequently, in S3, it is determined whether or not the vehicle speed V is zero. Assuming that it is not the case this time, a negative decision (NO) is obtained in S4, first, the front wheel, rear wheel cornering force F f, the current value of F r is calculated according to the equation. In this step, an expression corresponding to the expression (17) (the expression shown in FIG.
), Current values of ω, γ, G y , V and β (however, this time, ω and β are equal to the above initial values) and current calculations of F f and Fr respectively The current value of ω is updated by substituting the value. Note that e 1 , e 2 , e 3 ,
e 4 , a f , ar and I Z are constants, respectively. Thereafter, in S5, the current value of β is updated by calculating the sum of the current value of ω and the current value of pγ, and in S6
, The current value of β is stored in the RAM as the current estimated value of the vehicle side slip angle β. This completes one execution of this routine.

【0090】これに対して、車速Vが0である場合に
は、S3の判定がYESとなり、S7において、ωの現
在値もβの現在値も0とされ、その後、S6において、
そのβの現在値すなわち0が車体横すべり角βの今回推
定値としてRAMにストアされる。
On the other hand, if the vehicle speed V is 0, the determination in S3 is YES, and in S7, the current value of ω and the current value of β are also set to 0, and thereafter, in S6,
The current value of β, that is, 0, is stored in the RAM as the current estimated value of the vehicle body slip angle β.

【0091】なお、本ルーチンの2回目以後の各々の実
行時には、S1の判定がNOとなってS2の実行がスキ
ップされるが、車速Vが0でなくS3の判定がNOとな
ってS4が実行される場合には、このS4において、本
ルーチンの前回の実行において取得されたωおよびβの
それぞれの値、すなわち、ωおよびβのそれぞれの現在
値と、γ,Gy ,V,Ff およびFr のそれぞれの現在
値とが演算式に代入されることによってωの現在値が更
新され、さらに、S5において、そのωの現在値と、p
とヨーレートγの現在値との積との和が演算されること
によってβの現在値が更新されることになる。
In each of the second and subsequent executions of this routine, the determination in S1 is NO and the execution of S2 is skipped. However, the vehicle speed V is not 0 and the determination in S3 is NO and S4 is executed. If executed, in S4, the respective values of ω and β obtained in the previous execution of this routine, that is, the current values of ω and β, and γ, G y , V, F f and the current value of omega by each and the current value is assigned to the arithmetic expression of F r is updated, further, in S5, the current value of the omega, p
The current value of β is updated by calculating the sum of the product of the current value and the current value of the yaw rate γ.

【0092】一方、図3の後輪舵角制御ルーチンは、 δr =Kδf +λγ なる式に従って後輪舵角δr の目標値を決定し、それが
実現されるようにステップモータ20を制御するもので
ある。この式においてKは舵角ゲインであって、図4の
マップに従って車速Vに応じて決定される。また、λは
ヨーレートゲインであって、図5のグラフに従って車速
Vに応じて決定される第1部分ゲインλ1 と図6のグラ
フに従って車体横すべり角βの絶対値に応じて決定され
る第2部分ゲインλ2 との和として決定される。なお、
図4〜図6のグラフは、車両が定常状態にあるか過渡状
態にあるかを問わず車体横すべり角βの実際値ができる
限り0に近づくようにするための特性を持っている。
[0092] On the other hand, wheel steering angle control routine after 3 determines a target value of the rear wheel steering angle [delta] r in accordance with δ r = Kδ f + λγ made wherein control step motor 20 as it is implemented Is what you do. In this equation, K is a steering angle gain, which is determined according to the vehicle speed V according to the map of FIG. Λ is a yaw rate gain, and the second partial gain λ 1 is determined according to the vehicle speed V according to the graph of FIG. 5 and the second partial gain λ 1 is determined according to the absolute value of the vehicle body slip angle β according to the graph of FIG. It is determined as the sum of the partial gains lambda 2. In addition,
4 to 6 have characteristics to make the actual value of the vehicle body slip angle β approach 0 as much as possible regardless of whether the vehicle is in a steady state or a transient state.

【0093】図3の後輪舵角制御ルーチンも繰り返し実
行される。なお、図1の車体横すべり角推定ルーチンは
図3の後輪舵角制御ルーチンの一回の実行において一回
だけ実行されるようになっていて、結局、それらルーチ
ンは互いに同期して実行されるようになっている。
The rear wheel steering angle control routine of FIG. 3 is also repeatedly executed. Note that the vehicle body slip angle estimation routine of FIG. 1 is executed only once in one execution of the rear wheel steering angle control routine of FIG. 3, and eventually these routines are executed in synchronization with each other. It has become.

【0094】図3の後輪舵角制御ルーチンの各回の実行
時には、まず、S101において、各種センサから前輪
舵角δf ,ヨーレートγ,横加速度Gy および車速Vが
取り込まれ、続いて、S102において、車速Vに対応
する舵角ゲインKの値が図4のマップを用いて、車速V
に対応する第1部分ゲインλ1 の値が図5のマップを用
いてそれぞれ決定され、S103において、図1の車体
横すべり角推定ルーチンが実行されて車体横すべり角β
の現在値が推定される。その後、S104において、そ
の車体横すべり角βの現在値の絶対値に対応する第2部
分ゲインλ2 の値が図6のマップを用いて決定され、さ
らに、第1部分ゲインλ1 とその第2部分ゲインλ2
の和がヨーレートゲインλの今回値に決定される。続い
て、S105において、舵角ゲインKと前輪舵角δf
の積とヨーレートゲインλとヨーレートγとの積との和
が後輪舵角δr の今回の目標値に決定され、S106に
おいて、その今回の目標値が実現されるようにステップ
モータ20に対して駆動信号が出力される。以上で本ル
ーチンの一回の実行が終了する。
[0094] The wheel steering angle at each time of execution of the control routine after 3, first, in S101, the front wheel steering angle [delta] f from the various sensors, yaw rate gamma, lateral acceleration G y and the vehicle speed V is taken, followed by, S102 , The value of the steering angle gain K corresponding to the vehicle speed V is calculated using the map shown in FIG.
The first portion gain lambda 1 value corresponding to the determined respectively by using a map of Figure 5, in S103, the vehicle body slip angle is executing the vehicle body slip angle estimation routine of FIG. 1 beta
Is estimated. Thereafter, in S104, the value of the second partial gain λ 2 corresponding to the absolute value of the current value of the vehicle body slip angle β is determined using the map of FIG. 6, and further, the first partial gain λ 1 and its second the sum of the partial gains lambda 2 is determined to the current value of the yaw rate gain lambda. Subsequently, in S105, the sum of the product of the product and the yaw rate gain λ and the yaw rate γ and the steering angle gain K and the front wheel steering angle [delta] f is determined as the current target value of the rear wheel steering angle [delta] r, in S106 A drive signal is output to the step motor 20 so that the current target value is realized. This completes one execution of this routine.

【0095】以上の説明から明らかなように、本実施例
においては、各種センサによる検出値が積分されること
なくオブザーバにより車体横すべり角βが推定されるか
ら、その推定値が検出値のノイズ,ドリフト等の影響を
それほど受けずに済み、推定値の信頼性すなわち精度を
十分に高めることができるという効果が得られる。
As is clear from the above description, in this embodiment, the vehicle body slip angle β is estimated by the observer without integrating the detection values from the various sensors, and the estimated value is used as the noise of the detection value. The effect that the influence of the drift or the like is not so large can be obtained, and the reliability of the estimated value, that is, the accuracy can be sufficiently improved.

【0096】さらに、本実施例においては、車体横すべ
り角βがオブザーバにより推定され、しかも、このオブ
ザーバにおいては車体横すべり角βの推定誤差の減衰速
度(収束速度)が設計パラメータとしての行列Lすなわ
ちオブザーバの極(−p)の値の如何によって任意に設
定できる。そのため、推定誤差の減衰速度を比較的自由
かつ簡単に望ましい値に近づけることができるという効
果も得られる。
Further, in this embodiment, the vehicle body slip angle β is estimated by the observer, and in this observer, the decay speed (convergence speed) of the estimation error of the vehicle body slip angle β is determined by the matrix L as a design parameter, that is, the observer. Can be arbitrarily set depending on the value of the pole (−p) of Therefore, an effect is obtained that the decay speed of the estimation error can be relatively freely and easily approximated to a desired value.

【0097】以上の説明から明らかなように、本実施例
においては、コントローラ50のうち図1の車体横すべ
り角推定ルーチンを実行する部分が、車速センサ34,
ヨーレートセンサ36および横加速度センサ40と共同
して本発明に係る車体横すべり角推定装置を構成してお
り、このうち、コントローラ50のうち図1の車体横す
べり角推定ルーチンを実行する部分が、本発明における
「車体横すべり角推定手段」の一態様を構成しているの
である。
As is apparent from the above description, in this embodiment, the part of the controller 50 that executes the vehicle body slip angle estimation routine of FIG.
The yaw rate sensor 36 and the lateral acceleration sensor 40 cooperate with each other to constitute a vehicle body slip angle estimating apparatus according to the present invention, and the part of the controller 50 that executes the vehicle body slip angle estimating routine of FIG. Of the vehicle body slip angle estimating means.

【0098】なお付言すれば、本実施例においては、車
両の平面運動が線形領域にあるとして前輪,後輪コーナ
リングパワーCf ,Cr がそれぞれ定数とされていた
が、例えば、車速V,車輪すべり角αf ,αr 等の車両
運動状態量を変数とする関数として定義したり、前輪コ
ーナリングパワーCf と前輪すべり角αf との関係と、
後輪コーナリングパワーCr と後輪すべり角αr との関
係とをそれぞれ予め取得してマップ等としてROMにス
トアしておくこともでき、このようにすれば、車両の平
面運動が線形領域にあるか非線形領域にあるかを問わ
ず、前輪コーナリングフォースFf も後輪コーナリング
フォースFr も確定することができることになる。
It should be noted that, in this embodiment, the front wheel and rear wheel cornering powers C f and C r are assumed to be constants assuming that the plane motion of the vehicle is in the linear region. It is defined as a function using vehicle motion state quantities such as slip angles α f and α r as variables, and the relationship between the front wheel cornering power C f and the front wheel slip angle α f is defined as
The relationship between the rear wheel cornering power C r and the rear wheel slip angle α r can be obtained in advance and stored in a ROM or the like as a map or the like. regardless of whether a certain or non-linear region, also will be able to determine the rear wheel cornering force F r also the front wheel cornering force F f.

【0099】さらに付言すれば、本実施例は、車両の運
動状態を制御する車両運動制御装置の一例である後輪舵
角制御装置において用いられる車体横すべり角βを推定
することを目的とする車体横すべり角推定装置であった
が、他の車両運動制御装置、すなわち、例えば、車体の
進行方向と向きとができる限り互いに一致するようにエ
ンジンの駆動トルクを前輪12と後輪22とにそれぞれ
配分する比率を制御する駆動トルク配分制御装置等にお
いて用いられる車体横すべり角βを推定することを目的
とする車体横すべり角推定装置として本発明を実施する
こともできる。
[0099] In addition, the present embodiment aims at estimating the vehicle side slip angle β used in a rear wheel steering angle control device which is an example of a vehicle motion control device for controlling the motion state of the vehicle. Although it is a side slip angle estimating apparatus, another vehicle motion control apparatus, that is, for example, distributes the driving torque of the engine to the front wheels 12 and the rear wheels 22 so that the traveling direction and the direction of the vehicle body match each other as much as possible. The present invention can also be implemented as a vehicle body slip angle estimating apparatus for the purpose of estimating a vehicle body slip angle β used in a drive torque distribution control device or the like for controlling a ratio to be controlled.

【0100】以上、本発明の一実施例を図面に基づいて
詳細に説明したが、この他にも、特許請求の範囲を逸脱
することなく、当業者の知識に基づいて種々の変形,改
良を施した態様で本発明を実施することができる。
While the preferred embodiment of the present invention has been described in detail with reference to the drawings, various modifications and improvements may be made based on the knowledge of those skilled in the art without departing from the scope of the claims. The present invention can be carried out in the mode in which it is performed.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の一実施例である車体横すべり角推定装
置が用いる車体横すべり角推定ルーチンを示すフローチ
ャートである。
FIG. 1 is a flowchart illustrating a vehicle body slip angle estimation routine used by a vehicle body slip angle estimation apparatus according to an embodiment of the present invention.

【図2】その車体横すべり角推定装置を含む後輪舵角制
御装置を示すシステム図である。
FIG. 2 is a system diagram showing a rear wheel steering angle control device including the vehicle body slip angle estimation device.

【図3】その後輪舵角制御装置が用いる後輪舵角制御ル
ーチンを示すフローチャートである。
FIG. 3 is a flowchart showing a rear wheel steering angle control routine used by a rear wheel steering angle control device.

【図4】図3の後輪舵角制御ルーチンにおいて用いられ
る車速Vと舵角ゲインKとの関係を示すグラフである。
FIG. 4 is a graph showing a relationship between a vehicle speed V and a steering angle gain K used in a rear wheel steering angle control routine of FIG.

【図5】図3の後輪舵角制御ルーチンにおいて用いられ
る車速Vと第1部分ゲインλ1 との関係を示すグラフで
ある。
FIG. 5 is a graph showing a relationship between a vehicle speed V and a first partial gain λ 1 used in a rear wheel steering angle control routine of FIG. 3;

【図6】図3の後輪舵角制御ルーチンにおいて用いられ
る車体横すべり角βの絶対値と第2部分ゲインλ2 との
関係を示すグラフである。
6 is a graph showing a relationship between an absolute value of a vehicle body side slip angle β and a second partial gain λ 2 used in the rear wheel steering angle control routine of FIG. 3;

【符号の説明】[Explanation of symbols]

34 車速センサ 36 ヨーレートセンサ 40 横加速度センサ 50 コントローラ 34 vehicle speed sensor 36 yaw rate sensor 40 lateral acceleration sensor 50 controller

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.6,DB名) B62D 6/00 B62D 101:00 B62D 111:00 B62D 137:00 ──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int. Cl. 6 , DB name) B62D 6/00 B62D 101: 00 B62D 111: 00 B62D 137: 00

Claims (1)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 車体に発生する横加速度を検出する横加
速度センサと、 車体の走行速度である車速を検出する車速センサと、 車体に発生するヨーレートを検出するヨーレートセンサ
と、 それらに接続されて車両重心点における車体横すべり角
βを逐次推定する車体横すべり角推定手段であって、車
両重心点の回りに発生するヨーイングモーメントYMを
車両のヨー慣性モーメントIZ で割った変数であるYM
/IZ の現在値を取得するとともに、車両の横すべり運
動とヨーイング運動とから成る車両の平面運動に関して
想定される、車体横すべり角βとヨーレートγとをそれ
ぞれ状態変数とし、車両重心点における横加速度Gy
車速Vで割った変数であるGy /Vと前記YM/IZ
ら車体横すべり角βの前回推定値を差し引いた変数とを
それぞれ入力変数とし、ヨーレートγを出力変数とする
可観測な線形定係数システムにおいてオブザーバを用い
て、前記横加速度センサ,車速センサおよびヨーレート
センサの各々からの出力信号に基づく横加速度Gy ,車
速Vおよびヨーレートγの現在値と、前記取得したYM
/IZ の現在値と、車体横すべり角βの前回推定値とに
基づいて車体横すべり角βの今回値を推定するものとを
含むことを特徴とする車体横すべり角推定装置。
1. A lateral acceleration sensor for detecting a lateral acceleration generated in a vehicle body, a vehicle speed sensor for detecting a vehicle speed which is a running speed of the vehicle body, a yaw rate sensor for detecting a yaw rate generated in the vehicle body, and connected thereto. A vehicle slip angle estimating means for sequentially estimating a vehicle slip angle β at a vehicle center of gravity, wherein YM is a variable obtained by dividing a yawing moment YM generated around the vehicle center of gravity by a yaw inertia moment I Z of the vehicle.
/ I Z is obtained, and the vehicle body slip angle β and the yaw rate γ, which are assumed for the plane motion of the vehicle including the vehicle slip motion and the yawing motion, are set as state variables, respectively, and the lateral acceleration at the center of gravity of the vehicle is obtained. the G y and each input variable and a variable obtained by subtracting the previous estimated value of the G y / V is a variable divided by the vehicle speed V the YM / I Z vehicle body slip angle from the beta, observable to output variable yaw rate γ The linear acceleration coefficient system uses an observer to calculate the current values of the lateral acceleration G y , the vehicle speed V and the yaw rate γ based on the output signals from the lateral acceleration sensor, the vehicle speed sensor and the yaw rate sensor, and the obtained YM
A vehicle slip angle estimating apparatus characterized by including a device that estimates a current value of the vehicle body slip angle β based on a current value of / I Z and a previous estimated value of the vehicle body slip angle β.
JP4020592A 1992-01-09 1992-01-09 Body slip angle estimation device Expired - Lifetime JP2962025B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP4020592A JP2962025B2 (en) 1992-01-09 1992-01-09 Body slip angle estimation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP4020592A JP2962025B2 (en) 1992-01-09 1992-01-09 Body slip angle estimation device

Publications (2)

Publication Number Publication Date
JPH05185948A JPH05185948A (en) 1993-07-27
JP2962025B2 true JP2962025B2 (en) 1999-10-12

Family

ID=12031524

Family Applications (1)

Application Number Title Priority Date Filing Date
JP4020592A Expired - Lifetime JP2962025B2 (en) 1992-01-09 1992-01-09 Body slip angle estimation device

Country Status (1)

Country Link
JP (1) JP2962025B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014148262A (en) * 2013-02-01 2014-08-21 Toyota Motor Corp Controller for vehicle

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7885750B2 (en) * 2006-08-30 2011-02-08 Ford Global Technologies Integrated control system for stability control of yaw, roll and lateral motion of a driving vehicle using an integrated sensing system to determine a sideslip angle

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014148262A (en) * 2013-02-01 2014-08-21 Toyota Motor Corp Controller for vehicle

Also Published As

Publication number Publication date
JPH05185948A (en) 1993-07-27

Similar Documents

Publication Publication Date Title
JP3060923B2 (en) Vehicle state estimation device
US7845218B2 (en) Tire state estimator and tire state estimation method
US4898431A (en) Brake controlling system
US5275475A (en) Method for controlling vehicle dynamics
EP0983919B1 (en) A method for detecting a bank angle experienced by a moving vehicle
EP1627790B1 (en) Estimating method for road friction coefficient and vehicle slip angle estimating method
JP2932589B2 (en) Vehicle motion control device
US6615124B1 (en) Vehicular dynamic controlling apparatus and method
JP4926715B2 (en) Method and apparatus for assisting a vehicle operator in stabilizing a vehicle
US6745112B2 (en) Method of estimating quantities that represent state of vehicle
US6922617B2 (en) Vehicle control device and vehicle control method
JP2002145037A (en) Road surface friction coefficient estimation device for vehicles
JP2003306092A (en) Method for estimating vehicle state quantity
CN114375269B (en) Apparatus and method for estimating road friction coefficient
JP5211995B2 (en) Vehicle deceleration control apparatus and method
JPH08202403A (en) Vehicle state quantity estimation device
JP3191708B2 (en) Vehicle skid state quantity detection device
JP2962025B2 (en) Body slip angle estimation device
JP3271956B2 (en) Road surface friction coefficient estimation device for vehicles
JPH06221968A (en) Road surface friction coefficient detection device
US20080167777A1 (en) Method for Controlling the Steering Orientation of a Vehicle
JP6428497B2 (en) Vehicle control device
JP2962026B2 (en) Body slip angle estimation device
JPH07246948A (en) Compensator of quantity of vehicle motion state
JP3669174B2 (en) Vehicle motion control device

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070806

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080806

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080806

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090806

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100806

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100806

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110806

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110806

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120806

Year of fee payment: 13

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120806

Year of fee payment: 13