JPH11139293A - Device for estimating condition of road surface - Google Patents
Device for estimating condition of road surfaceInfo
- Publication number
- JPH11139293A JPH11139293A JP9326921A JP32692197A JPH11139293A JP H11139293 A JPH11139293 A JP H11139293A JP 9326921 A JP9326921 A JP 9326921A JP 32692197 A JP32692197 A JP 32692197A JP H11139293 A JPH11139293 A JP H11139293A
- Authority
- JP
- Japan
- Prior art keywords
- road surface
- speed
- physical quantity
- braking force
- wheel
- 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.)
- Granted
Links
Landscapes
- Regulating Braking Force (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、路面状態推定装置
に係り、より詳しくは、車両が走行する路面の路面状態
を推定する路面状態推定装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a road surface condition estimating device, and more particularly, to a road surface condition estimating device for estimating a road surface condition of a road surface on which a vehicle travels.
【0002】[0002]
【従来の技術】従来より、車輪と路面との間の摩擦係数
μがピーク値を超えて車輪がロック状態に移行する直前
に、車輪に作用するブレーキトルクを低下させることに
よって、車輪のロックを防止しピークμ値に追従制御す
るアンチロックブレーキ制御装置が提案されている。2. Description of the Related Art Locking of a wheel has been conventionally performed by reducing brake torque acting on the wheel immediately before the wheel enters a locked state after the coefficient of friction μ between the wheel and the road surface exceeds a peak value. There has been proposed an anti-lock brake control device for preventing and following the peak μ value.
【0003】ところで、車両がある速度で走行している
時、ブレーキをかけていくと車輪と路面との間にスリッ
プが生じるが、車輪と路面との間の摩擦係数μは、
(1)式で表されるスリップ率sに対し、図20のよう
に変化することが知られている。[0003] When a vehicle is running at a certain speed, a slip occurs between the wheels and the road surface when the brakes are applied, and the friction coefficient μ between the wheels and the road surface is as follows.
It is known that the slip ratio s represented by the equation (1) changes as shown in FIG.
【0004】[0004]
【数1】s = (V−Vw )/V・・・(1) ただし、Vは車体速度(角速度換算)、Vw は車輪速度
であり、よって、(V−Vw )はスリップ速度ΔVとな
る。S = (V−V w ) / V (1) where V is the vehicle speed (in terms of angular speed), V w is the wheel speed, and (V−V w ) is the slip speed. ΔV.
【0005】図20に示すように、このμ−s特性で
は、あるスリップ率(図20のA2領域)で摩擦係数μ
がピーク値をとるようになる。As shown in FIG. 20, in the μ-s characteristic, at a certain slip ratio (A2 region in FIG. 20), the friction coefficient μ
Takes a peak value.
【0006】そこで、特開平1−249559号公報に
は、車体速度の近似値、及び検出された車輪速度から
(1)式よりスリップ率を演算し、演算したスリップ率
が、予め設定してある基準スリップ率(ピークμを与え
るスリップ率)に略一致するように、ブレーキ力を制御
することにより、ピークμに追従するアンチロックブレ
ーキ制御装置が提案されている。Therefore, in Japanese Patent Application Laid-Open No. 1-249559, the slip ratio is calculated from the approximate value of the vehicle body speed and the detected wheel speed by the equation (1), and the calculated slip ratio is set in advance. An anti-lock brake control device that follows a peak μ by controlling a braking force so as to substantially match a reference slip ratio (a slip ratio giving a peak μ) has been proposed.
【0007】[0007]
【発明が解決しようとする課題】しかしながら、車輪が
走行する路面状態によってピークμとなるスリップ率が
異なる。このため、上記公報記載の発明のように固定さ
れた基準スリップ率の追従制御を行った場合、路面状態
によって制動距離が大きくなり過ぎたり、或いはピーク
μを超えてブレーキ制動されることによりタイヤロック
が発生するおそれがある。However, the slip ratio at which the peak μ is obtained differs depending on the condition of the road surface on which the wheels travel. For this reason, when the follow-up control of the fixed reference slip ratio is performed as in the invention described in the above-mentioned publication, the braking distance becomes too large depending on the road surface condition, or the tire is locked due to brake braking exceeding the peak μ. May occur.
【0008】この対策として路面状態を推定演算し、演
算された路面状態に応じて基準スリップ率を変化させる
必要があるが、路面状態を精度よく推定する技術がな
い。As a countermeasure, it is necessary to estimate and calculate the road surface condition and change the reference slip ratio according to the calculated road surface condition. However, there is no technology for accurately estimating the road surface condition.
【0009】本発明は、上記事実に鑑みて成されたもの
で、路面状態を精度よく推定することの可能な路面状態
推定装置を提供することを目的とする。The present invention has been made in view of the above-described circumstances, and has as its object to provide a road surface condition estimating apparatus capable of accurately estimating a road surface condition.
【0010】[0010]
【課題を解決するための手段】上記目的達成のため本発
明は、車両の車輪の制動状態を表す第1の物理量を検出
する第1の検出手段と、前記車両の速度を検出する第2
の検出手段と、前記車輪のすべり易さを表しかつ前記車
輪が走行する路面状態、前記第1の物理量、及び車両の
速度に応じて定まる第2の物理量を検出する第3の検出
手段と、前記第1の検出手段により検出された第1の物
理量及び前記第2の検出手段により検出された車速に基
づいて、予め定められた複数の路面状態毎に第2の物理
量を演算する演算手段と、前記車輪の制動状態を表す第
1の物理量に対する前記第3の検出手段により検出され
た第2の物理量の軌跡と、前記第1の物理量に対する前
記演算手段により複数の路面状態毎に演算された第2の
物理量の軌跡の各々と、を比較する比較手段と、前記比
較手段の比較結果に基づいて、前記演算された軌跡の内
の前記検出された軌跡に最も近い軌跡に対応する路面状
態を、前記車輪が走行している実際の路面状態として推
定する推定手段と、を備えている。In order to achieve the above object, the present invention provides a first detecting means for detecting a first physical quantity representing a braking state of a vehicle wheel, and a second detecting means for detecting a speed of the vehicle.
Detection means, and a third detection means for detecting the road surface state on which the wheels travel and the second physical quantity determined according to the speed of the vehicle, indicating the ease of slipping of the wheels, Calculating means for calculating a second physical quantity for each of a plurality of predetermined road surface conditions based on the first physical quantity detected by the first detecting means and the vehicle speed detected by the second detecting means; A trajectory of a second physical quantity detected by the third detecting means with respect to a first physical quantity representing a braking state of the wheel, and a trajectory calculated for each of a plurality of road surface states by the calculating means with respect to the first physical quantity. Comparing means for comparing each of the trajectories of the second physical quantity, and a road surface state corresponding to a trajectory closest to the detected trajectory among the calculated trajectories based on a comparison result of the comparing means. , Said wheels Includes an estimating means for estimating an actual road surface condition that line, the.
【0011】第1の検出手段は、車両の車輪の制動状態
を表す第1の物理量を検出する。ここで、第1の物理量
は、車輪に対する制動力やスリップ速度がある。The first detecting means detects a first physical quantity representing a braking state of a vehicle wheel. Here, the first physical quantity includes a braking force and a slip speed on the wheel.
【0012】第2の検出手段は、車両の速度(車速)を
検出する。第3の検出手段は、車輪のすべり易さを表す
第2の物理量を検出する。この第2の物理量は、車輪が
走行する路面状態、第1の物理量、及び車両の速度に応
じて定まる。具体的には、第3の物理量は、路面と車輪
との間の摩擦係数又は車輪に対する制動トルクのすべり
速度に対する勾配である。即ち、これらの勾配が高い
と、車輪の周速が大きく、すべりにくい。一方、これら
の勾配が低いと、車輪の周速が小さく、すべり易い。よ
って、これらの勾配は、車輪のすべり易さを表す。The second detecting means detects the speed (vehicle speed) of the vehicle. The third detecting means detects a second physical quantity representing the ease of slipping of the wheel. The second physical quantity is determined according to the road surface condition on which the wheels travel, the first physical quantity, and the speed of the vehicle. Specifically, the third physical quantity is a coefficient of friction between the road surface and the wheels or a gradient of the braking torque for the wheels with respect to the slip speed. That is, when these gradients are high, the peripheral speed of the wheels is high, and slip is difficult. On the other hand, when these gradients are low, the peripheral speed of the wheels is low, and the wheels tend to slip. Thus, these slopes represent the ease with which the wheels slide.
【0013】演算手段は、第1の検出手段により検出さ
れた第1の物理量及び第2の検出手段により検出された
車速に基づいて、予め定められた複数の路面状態毎に第
2の物理量を演算する。即ち、演算手段は、第1の物理
量と、第2の物理量を演算するための第3の物理量と、
の車速に応じて定まる関係を複数の路面状態毎に予め求
めておき、検出された第1の物理量及び検出された車速
に基づいて、複数の路面状態毎に第3の物理量を演算
し、演算した第3の物理量から複数の路面状態毎に第2
の物理量を演算してもよい。また、演算手段は、第1の
物理量と第2の物理量との車速に応じて定まる関係を複
数の路面状態毎に予め求めておき、検出された第1の物
理量及び車速に基づいて、複数の路面状態毎に第2の物
理量を演算してもよい。The calculating means calculates a second physical quantity for each of a plurality of predetermined road surface conditions based on the first physical quantity detected by the first detecting means and the vehicle speed detected by the second detecting means. Calculate. That is, the calculating means includes a first physical quantity, a third physical quantity for calculating the second physical quantity,
Is determined in advance for each of the plurality of road surface states, and a third physical quantity is calculated for each of the plurality of road surface states based on the detected first physical quantity and the detected vehicle speed. From the third physical quantity obtained, the second
May be calculated. Further, the calculating means obtains in advance a relationship determined according to the vehicle speed between the first physical quantity and the second physical quantity for each of a plurality of road surface conditions, and based on the detected first physical quantity and the detected vehicle speed, obtains a plurality of relations. The second physical quantity may be calculated for each road surface condition.
【0014】比較手段は、車輪の制動状態を表す第1の
物理量に対する第3の検出手段により検出された第2の
物理量の軌跡と、第1の物理量に対する演算手段により
複数の路面状態毎に演算された第2の物理量の軌跡の各
々と、を比較する。推定手段は、比較手段の比較結果に
基づいて、上記演算された軌跡の内の上記検出された軌
跡に最も近い軌跡に対応する路面状態を、車輪が走行し
ている実際の路面状態として推定する。The comparing means calculates the trajectory of the second physical quantity detected by the third detecting means with respect to the first physical quantity representing the braking state of the wheel, and calculates the trajectory of the first physical quantity for each of a plurality of road surface conditions by the calculating means. With each of the trajectories of the second physical quantity. The estimating means estimates a road surface state corresponding to a trajectory closest to the detected trajectory in the calculated trajectories as an actual road surface state in which wheels are traveling, based on a comparison result of the comparing means. .
【0015】このように、第2の物理量を検出し、検出
された第2の物理量の軌跡と、複数の路面状態各々に対
応して演算された第2の物理量の軌跡の各々と、を比較
するので、比較結果に基づいて、複数の路面状態各々に
対応して演算された軌跡の内の検出された軌跡に最も近
い軌跡を選択すれば、選択した軌跡に対応する路面状態
は車輪が走行している実際の路面状態に対応し、実際の
路面状態を精度よく推定することができる。In this manner, the second physical quantity is detected, and the trajectory of the detected second physical quantity is compared with each of the trajectories of the second physical quantity calculated corresponding to each of the plurality of road surface conditions. Therefore, based on the comparison result, if a trajectory closest to the detected trajectory is selected from among the trajectories calculated corresponding to each of the plurality of road surface conditions, the road surface condition corresponding to the selected trajectory is determined by the traveling of the wheels. In response to the actual road surface state, the actual road surface state can be accurately estimated.
【0016】ところで、上記第2の物理量は、種々の原
因により誤差が含まれる。そこで、第1の検出手段は第
1の物理量を、第2の検出手段は車速を、第3の検出手
段は第2の物理量を、車輪に対して制動が開始されたと
きから所定時間毎に検出し、演算手段は、複数の路面状
態各々毎に第2の物理量を車輪に対して制動が開始され
たときから第1の物理量が所定量変化する毎に演算す
る。Incidentally, the second physical quantity includes an error due to various causes. Therefore, the first detecting means calculates the first physical quantity, the second detecting means calculates the vehicle speed, and the third detecting means calculates the second physical quantity at predetermined time intervals from when the braking is started on the wheels. The detecting means calculates the second physical quantity for each of the plurality of road surface conditions each time the first physical quantity changes by a predetermined amount from the time when braking of the wheel is started.
【0017】そして、比較手段は、車輪に対して制動が
開始されたときから第1の物理量が所定量変化する毎
に、該所定時間毎に検出された第2の物理量と複数の路
面状態各々毎及び第1の物理量が所定量変化する毎に演
算された第2の物理量各々とを比較するようにしても
い。Each time the first physical quantity changes by a predetermined amount from the time when braking of the wheel is started, the comparing means compares the second physical quantity detected at each of the predetermined time and the plurality of road surface conditions with each other. Each time and each time the first physical quantity changes by a predetermined amount, the second physical quantity calculated may be compared.
【0018】このように、車輪に対して制動が開始され
たときから第1の物理量が所定量変化する毎に、所定時
間毎に検出された第2の物理量と、複数の路面状態各々
毎及び第1の物理量が所定量変化する毎に演算された物
理量各々と、を比較するので、誤差の影響を少なくする
ことができ、より精度よく路面状態を推定することがで
きる。As described above, every time the first physical quantity changes by a predetermined amount from the time when the braking is started to the wheels, the second physical quantity detected at a predetermined time interval and the plurality of road surface conditions and Since each physical quantity calculated every time the first physical quantity changes by a predetermined amount is compared, the influence of an error can be reduced, and the road surface state can be estimated more accurately.
【0019】[0019]
【発明の実施の形態】以下、本発明の実施の形態を図面
を参照して詳細に説明する。Embodiments of the present invention will be described below in detail with reference to the drawings.
【0020】図1に示すように、本実施の形態に係る路
面状態推定装置は、車輪速度を検出する車輪速センサ2
8、車輪に作用する制動力を検出する制動力検出部3
4、及び路面μ勾配として後述する微小ゲインを演算す
る微小ゲイン演算部36、及び路面状態及びすべり速度
を演算する演算回路10を備えている。As shown in FIG. 1, a road surface state estimating apparatus according to the present embodiment comprises a wheel speed sensor 2 for detecting a wheel speed.
8. Braking force detector 3 for detecting the braking force acting on the wheels
4 and a minute gain calculating unit 36 for calculating a minute gain, which will be described later, as a road surface μ gradient, and an arithmetic circuit 10 for calculating a road surface state and a slip speed.
【0021】演算回路10は、予め定められた路面状
態、即ち、Dry、Snow、Iceに対応する路面μ
勾配の基準値を、制動力検出部34で検出された制動力
が所定量変化する毎に作成・記憶する基準値作成記憶回
路12D、12S、12Iと、前記制動力が所定量変化
する毎に微小ゲインを記憶する微小ゲイン記憶回路12
Jと、を備えている。基準値作成記憶回路12D、12
S、12I及び微小ゲイン記憶回路12Jには、基準値
作成記憶回路12D、12S、12Iに記憶された路面
μ勾配の基準値の時系列データと、実際に検出され、微
小ゲイン記憶回路12Jに記憶された路面μ勾配の基準
値の時系列データと、を照合して、路面状態を演算する
時系列データ軌跡照合部14が接続されている。時系列
データ軌跡照合部14の出力側には、すべり速度選択部
16が接続されている。The arithmetic circuit 10 determines a predetermined road surface condition, that is, a road surface μ corresponding to Dry, Snow and Ice.
A reference value creation storage circuit 12D, 12S, 12I for creating and storing a reference value of the gradient each time the braking force detected by the braking force detection unit 34 changes by a predetermined amount, and each time the braking force changes by a predetermined amount. Small gain storage circuit 12 for storing small gain
J. Reference value creation storage circuits 12D, 12
S, 12I and the minute gain storage circuit 12J store the time series data of the reference value of the road μ gradient stored in the reference value creation storage circuits 12D, 12S, 12I, and the time series data actually detected and stored in the small gain storage circuit 12J. A time-series data trajectory matching unit 14 for comparing the obtained time-series data of the reference value of the road surface μ gradient and calculating the road surface state is connected. The slip speed selecting unit 16 is connected to the output side of the time-series data trajectory matching unit 14.
【0022】基準値作成記憶回路12D、12S、12
Iは各々同一の構成であるので、基準値作成記憶回路1
2Dのみを説明し、他の説明を省略する。基準値作成記
憶回路12Dは、車輪速センサ28に接続された車速演
算部18を備えている。車速演算部18には、すべり速
度を演算するすべり速度演算部20が接続されている。
すべり速度演算部20には、制動力検出部34が接続さ
れている。すべり速度演算部20は、路面μ勾配の基準
値を作成するμ勾配基準値作成部22に接続されてい
る。μ勾配基準値作成部22には、車速演算部18及び
制動力検出部34が接続されている。μ勾配基準値作成
部22は、μ勾配基準値作成部22により作成された路
面μ勾配の基準値の時系列データを記憶する基準値時系
列データ記憶部24が接続されている。基準値時系列デ
ータ記憶部24は、時系列データ軌跡照合部14に接続
されている。Reference value creation storage circuits 12D, 12S, 12
Since I has the same configuration, the reference value creation storage circuit 1
Only 2D will be described, and other description will be omitted. The reference value creation storage circuit 12D includes a vehicle speed calculation unit 18 connected to the wheel speed sensor 28. A slip speed calculator 20 for calculating the slip speed is connected to the vehicle speed calculator 18.
A braking force detector 34 is connected to the slip speed calculator 20. The slip speed calculation unit 20 is connected to a μ gradient reference value creation unit 22 that creates a reference value of a road surface μ gradient. The vehicle speed calculator 18 and the braking force detector 34 are connected to the μ gradient reference value generator 22. The μ slope reference value creation unit 22 is connected to a reference value time series data storage unit 24 that stores time series data of the road surface μ slope reference value created by the μ slope reference value creation unit 22. The reference value time-series data storage unit 24 is connected to the time-series data track collation unit 14.
【0023】微小ゲイン記憶回路12Jは、検出された
微小ゲインの時系列データを記憶する微小ゲイン時系列
データ記憶部25を備えている。微小ゲイン時系列デー
タ記憶部25には、制動力検出部34及び微小ゲイン演
算部36が接続されている。微小ゲイン時系列データ記
憶部25は、時系列データ軌跡照合部14に接続されて
いる。The minute gain storage circuit 12J includes a minute gain time series data storage unit 25 for storing the time series data of the detected minute gain. The braking force detecting unit 34 and the minute gain calculating unit 36 are connected to the minute gain time-series data storage unit 25. The small-gain time-series data storage unit 25 is connected to the time-series data trajectory matching unit 14.
【0024】すべり速度選択部16には、基準値作成記
憶回路12D、12S、12Iの各すべり速度演算部2
0が接続されている。The slip speed selector 16 includes the respective slip speed calculators 2D, 12S, and 12I of the reference value creation storage circuits 12D, 12S, and 12I.
0 is connected.
【0025】制動力検出部34は、路面から車輪に対し
摩擦力として作用する制動力を、車輪の力学的モデルに
従って以下のように推定する。The braking force detector 34 estimates a braking force acting as a frictional force from the road surface to the wheels as follows according to a dynamic model of the wheels.
【0026】すなわち、車輪には、車輪に対し車輪の回
転方向と反対方向に作用するブレーキトルクTB と、車
輪に対し摩擦力として車輪の回転方向に作用する制動力
FによるタイヤトルクTf と、が作用する。ブレーキト
ルクTB は、車輪のブレーキディスクに対し車輪の回転
を妨げるように作用するブレーキ力に由来するものであ
り、制動力F及びタイヤトルクTf は、車輪と路面との
間の摩擦係数をμB 、車輪半径をr、車輪荷重をWとし
たとき、(2)式、(3)式によって表される。[0026] That is, the wheel, and the brake torque T B acting in the direction opposite to the rotational direction of the wheel to the wheel, and tire torque T f by the braking force F acting in the rotational direction of the wheel as a frictional force to the wheel , Works. Brake torque T B is derived from the braking force which acts to prevent rotation of the wheel relative to the wheel of the brake disc, the braking force F and the tire torque T f is the friction coefficient between the wheel and the road surface When μ B , the wheel radius is r, and the wheel load is W, it is expressed by the equations (2) and (3).
【0027】[0027]
【数2】F= μB W・・・(2)[Equation 2] F = μ B W (2)
【0028】[0028]
【数3】Tf = F×r = μB Wr・・・(3) 従って、車輪の運動方程式は、T f = F × r = μ B Wr (3) Therefore, the equation of motion of the wheel is
【0029】[0029]
【数4】 となる。ただし、Iは車輪の慣性モーメント、ωは車輪
の回転速度(車輪速度)である。(Equation 4) Becomes Here, I is the moment of inertia of the wheel, and ω is the rotational speed (wheel speed) of the wheel.
【0030】車輪加速度(dω/dt)を検知し、ブレ
ーキディスクに加えられるホイールシリンダ圧に基づい
てブレーキトルクTB を求めれば、(4) 式に基づいて制
動力Fを推定することができる。具体的には、アクセル
開度などから求めた車輪の駆動トルクと、外乱としての
制動力Fが車輪に作用する(4) 式と等価な力学モデルを
オブザーバとして構成する。このオブザーバでは、(4)
式を2階積分することにより得られる回転位置と実際に
検出された回転位置との偏差を0に一致させるように制
御周期毎に等価モデルの外乱及び回転速度を修正し、修
正された外乱を制動力として推定する。[0030] detecting the wheel acceleration (d [omega / dt), by obtaining a brake torque T B based on the wheel cylinder pressure applied to the brake disc, it is possible to estimate the braking force F based on the equation (4). Specifically, a dynamic model equivalent to the equation (4) in which the drive torque of the wheel determined from the accelerator opening and the like and the braking force F acting as a disturbance acts on the wheel is configured as an observer. In this observer, (4)
The disturbance and rotation speed of the equivalent model are corrected for each control cycle so that the deviation between the rotation position obtained by the second-order integration of the equation and the rotation position actually detected is equal to 0, and the corrected disturbance is corrected. Estimate as braking force.
【0031】次に、微小ゲイン演算部36を説明する。
まず、路面μ勾配としての微小ゲインを説明する。Next, the small gain calculating section 36 will be described.
First, a small gain as a road surface μ gradient will be described.
【0032】重量Wの車体を備えた車両が車体速度ωv
で走行している時の車輪での振動現象、すなわち車体と
車輪と路面とによって構成される振動系の振動現象を、
車輪回転軸で等価的にモデル化した図3に示すモデルを
参照して説明する。A vehicle having a body weight of W has a vehicle speed ω v
The vibration phenomenon at the wheels when traveling in, that is, the vibration phenomenon of the vibration system composed of the vehicle body, wheels and the road surface,
This will be described with reference to a model shown in FIG. 3 which is equivalently modeled by a wheel rotation axis.
【0033】図3のモデルにおいて、ブレーキ力は、路
面と接するタイヤのトレッド115の表面を介して路面
に作用する。しかし、このブレーキ力は実際には路面か
らの反作用(制動力)として車体に作用する。このた
め、車体重量の回転軸換算の等価モデル117は、タイ
ヤのトレッドと路面との間の摩擦要素116(路面μ)
を介して車輪113と反対側に連結したものとなる。こ
れは、シャシーダイナモ装置のように、車輪下の大きな
慣性、すなわち車輪と反対側の質量で車体の重量を模擬
することができることと同様である。In the model of FIG. 3, the braking force acts on the road surface via the surface of the tread 115 of the tire in contact with the road surface. However, this braking force actually acts on the vehicle body as a reaction (braking force) from the road surface. For this reason, the equivalent model 117 in terms of the rotation axis in terms of the vehicle body weight is a friction element 116 (road surface μ) between the tread of the tire and the road surface.
And connected to the opposite side of the wheel 113 via the. This is the same as the large inertia under the wheels, that is, the weight of the vehicle body can be simulated by the mass on the side opposite to the wheels, as in the chassis dynamo device.
【0034】図3でタイヤリムを含んだ車輪113の慣
性をJw 、リムとトレッド115との間のばね要素11
4のばね定数をK、車輪半径をR、トレッド115の慣
性をJt 、トレッド115と路面との間の摩擦要素11
6の摩擦係数をμ、車体の重量の回転軸換算の等価モデ
ル117の慣性をJV とすると、ホイールシリンダ圧に
より生じるブレーキトルクTb ’から車輪速ωw までの
伝達特性は、車輪運動の方程式より、In FIG. 3, the inertia of the wheel 113 including the tire rim is represented by J w , and the spring element 11 between the rim and the tread 115 is represented by J w .
4, the spring constant is K, the wheel radius is R, the inertia of the tread 115 is J t , and the friction element 11 between the tread 115 and the road surface.
The friction coefficient of the 6 mu, the inertia of the equivalent model 117 of the body weight of the rotary shaft conversion and J V, the transfer characteristic from the braking torque T b 'caused by the wheel cylinder pressure to the wheel speed omega w is the wheel motion From the equation,
【0035】[0035]
【数5】 (Equation 5)
【0036】となる。なお、sはラプラス変換の演算子
である。ここで、タイヤが路面にグリップしている時
は、トレッド115と車体等価モデル117とが直結さ
れていると考える。この場合、車体等価モデル117と
トレッド115との和の慣性と、車輪113の慣性とが
共振する。即ち、この振動系は、車輪と車体と路面とか
ら構成された車輪共振系とみなすことができる。このと
きの車輪共振系の共振周波数ω∞は、(5))式の伝達特性
において、## EQU1 ## Note that s is an operator of Laplace transform. Here, when the tire is gripping the road surface, it is considered that the tread 115 and the vehicle body equivalent model 117 are directly connected. In this case, the inertia of the sum of the vehicle equivalent model 117 and the tread 115 and the inertia of the wheels 113 resonate. That is, this vibration system can be considered as a wheel resonance system including the wheels, the vehicle body, and the road surface. The resonance frequency ω∞ of the wheel resonance system at this time is given by the transfer characteristic of equation (5).
【0037】[0037]
【数6】 (Equation 6)
【0038】となる。ここで、図20において(6) 式が
成立する摩擦状態は、ピークμに達する前の領域A1に
対応する。## EQU4 ## Here, the friction state in which the expression (6) is established in FIG. 20 corresponds to the region A1 before reaching the peak μ.
【0039】逆に、タイヤの摩擦係数μがピークμに近
づく場合には、タイヤ表面の摩擦係数μがスリップ率に
対して変化し難くなる。即ち、トレッド115の慣性の
振動に伴う成分は車体等価モデル117に影響しなくな
る。つまり等価的にトレッド115と車体等価モデル1
17とが分離され、トレッド115と車輪113とが共
振を起こすことになる。このときの車輪共振系は、車輪
と路面とから構成されているとみなすことができる。そ
の共振周波数ω∞’は、(6) 式において、車体等価慣性
Jv を0とおいたものと等しくなる。すなわち、Conversely, when the friction coefficient μ of the tire approaches the peak μ, the friction coefficient μ of the tire surface hardly changes with respect to the slip ratio. That is, the component accompanying the inertial vibration of the tread 115 does not affect the vehicle equivalent model 117. That is, the tread 115 is equivalent to the vehicle equivalent model 1
17 is separated, and the tread 115 and the wheel 113 resonate. At this time, the wheel resonance system can be regarded as being composed of wheels and a road surface. Its resonance frequency Omega∞ 'is equal to that placed in (6), the vehicle body equivalent inertia J v 0. That is,
【0040】[0040]
【数7】 (Equation 7)
【0041】となる。この状態は、図20では、ピーク
μ近傍の領域A2に対応する。なお、ピークμを越えて
ブレーキ制動されると、領域A3に瞬時に移行し、タイ
ヤがロックされる。## EQU1 ## This state corresponds to the area A2 near the peak μ in FIG. When the brake is applied beyond the peak μ, the operation instantaneously shifts to the area A3, and the tire is locked.
【0042】車体等価慣性Jv が車輪慣性Jw 、トレッ
ド慣性Jt より大きいと仮定する。この場合、(7) 式の
場合の車輪共振系の共振周波数ω∞’は(6) 式のω∞よ
りも高周波数側にシフトすることになる。[0042] Suppose the body equivalent inertia J v is the wheel inertia J w, greater than the tread inertia J t. In this case, the resonance frequency ω∞ ′ of the wheel resonance system in equation (7) shifts to a higher frequency side than ω∞ in equation (6).
【0043】ここで、車輪と車体と路面とからなる振動
系の共振周波数ω∞((7) 式) でブレーキ力を微小励振
すると(ここでは、ブレーキ圧Pb を微小励振するとす
る)、車輪速度ωw も平均的な車輪速度の回りに共振周
波数ω∞で微小振動する。ここで、このときのブレーキ
圧Pb の共振周波数ω∞の微小振幅をPv 、車輪速度の
共振周波数ω∞の微小振幅をωwvとした場合、微小ゲイ
ンGd を[0043] Here, when the micro exciting the braking force at the wheel and the vehicle body and the road surface and the vibration system of the resonance frequency ω∞ consisting ((7)) (in this case, a finely excite the brake pressure P b), the wheels The speed ω w also minutely vibrates around the average wheel speed at the resonance frequency ω∞. Here, small amplitude P v of the resonance frequency ω∞ brake pressure P b in this case, if the small amplitude of the resonance frequency ω∞ wheel speed to omega wv, the fine gain G d
【0044】[0044]
【数8】Gd =ωwv/Pv ・・・(8) とする。なお、この微小ゲインGd を、ブレーキ圧Pb
に対する車輪速ωw の比(ωw /Pb )の共振周波数ω
∞の振動成分とみなし、[Number 8 and G d = ω wv / P v ··· (8). It should be noted that this small gain G d is used as the brake pressure P b
The ratio of the wheel speed ω w with respect to (ω w / P b) of the resonance frequency ω
Considering the vibration component of ∞,
【0045】[0045]
【数9】 Gd =((ωw /Pb )|s=jω∞)・・・(9) と表すこともできる。G d = ((ω w / P b ) | s = jω∞) (9)
【0046】この微小ゲインGd は、(9) 式に示すよう
に(ωw /Pb )の共振周波数ω∞の振動成分であるの
で、摩擦状態がピークμ近傍の領域に至ったとき、共振
周波数がω∞’にシフトするため急激に減少する。すな
わち、微小ゲインGd は、路面μ特性を規定する物理量
であるといえる。Since the small gain G d is a vibration component of the resonance frequency ω∞ of (ω w / P b ) as shown in the equation (9), when the friction state reaches a region near the peak μ, Since the resonance frequency shifts to ω∞ ′, it sharply decreases. That is, it can be said that the minute gain Gd is a physical quantity that defines the road surface μ characteristic.
【0047】そして、微小ゲイン演算部36は、図2に
示すように、振動系の共振周波数ω∞((6) 式)でブレ
ーキ圧を微小励振したときの、車輪速度Vw の共振周波
数ω∞の微小振幅(車輪速微小振幅ωwv)を検出する車
輪速微小振幅検出部40と、共振周波数ω∞のブレーキ
圧の微小振幅Pv を検出するブレーキ圧微小振幅検出部
42と、検出された車輪速微小振幅ωwvをブレーキ圧微
小振幅Pv で除算することにより微小ゲインGd を出力
する除算器44と、から構成される。[0047] Then, fine gain calculation unit 36, as shown in FIG. 2, when the brake pressure and micro excited in the vibration system of the resonance frequency ω∞ ((6) formula), the resonance frequency of the wheel speed V w omega The wheel speed minute amplitude detecting unit 40 for detecting the minute amplitude of ∞ (wheel speed minute amplitude ω wv ), and the brake pressure minute amplitude detecting unit 42 for detecting the minute amplitude P v of the brake pressure at the resonance frequency ω∞ are detected. was composed of wheel speed small amplitude omega wv a divider 44 which outputs the micro-gain G d is divided by the braking pressure differential small amplitude P v, from.
【0048】ここで、車輪速微小振幅検出部40は、共
振周波数ω∞の振動成分を抽出するフィルタ処理を行う
図4のような演算部として実現できる。例えば、この振
動系の共振周波数ω∞が40[Hz]程度であるので、
制御性を考慮して1周期を24[ms]、約41.7
[Hz]に取り、この周波数を中心周波数とする帯域通
過フィルタ75を設ける。このフィルタにより、車輪速
度信号ωi から約41.7[Hz]近傍の周波数成分の
みが抽出される。さらに、このフィルタ出力を全波整流
器76により全波整流、直流平滑化し、この直流平滑化
信号から低域通過フィルタ77によって低域振動成分の
みを通過させることにより、車輪速微小振幅ωwvを出力
する。Here, the wheel speed minute amplitude detecting section 40 can be realized as an arithmetic section as shown in FIG. 4 which performs a filtering process for extracting a vibration component of the resonance frequency ω∞. For example, since the resonance frequency ω∞ of this vibration system is about 40 [Hz],
Considering the controllability, one cycle is 24 [ms], about 41.7.
[Hz], and a band-pass filter 75 having this frequency as a center frequency is provided. This filter only the frequency component of about 41.7 [Hz] from near the wheel speed signal omega i is extracted. Further, the filter output is full-wave rectified by a full-wave rectifier 76 and DC-smoothed, and only a low-frequency vibration component is passed from the DC-smoothed signal by a low-pass filter 77 to output a wheel speed minute amplitude ωwv . I do.
【0049】なお、周期の整数倍、例えば1周期の24
[ms]、2周期の48[ms]の時系列データを連続
的に取り込み、41.7[Hz]の単位正弦波、単位余
弦波との相関を求めることによっても車輪速微小振幅検
出部40を実現できる。Incidentally, an integral multiple of the cycle, for example, 24 of one cycle
[Ms], the wheel speed minute amplitude detection unit 40 can also be obtained by continuously taking in time series data of 48 [ms] in two cycles and obtaining a correlation with a unit sine wave and a unit cosine wave of 41.7 [Hz]. Can be realized.
【0050】ここで、平均ブレーキ圧Pm の回りに共振
周波数のブレーキ圧微小振幅Pv を印加する微小励振手
段について説明する。まず、平均ブレーキ圧指令及び微
小励振指令を実際の車輪への制動トルクに変換する部分
(バルブ制御系)は、図5に示すように、マスタシリン
ダ48、制御バルブ52、ホイールシリンダ56、リザ
ーバー58及びオイルポンプ60を備えている。[0050] Here, a description will be given small excitation means for applying a brake pressure differential small amplitude P v of the resonance frequency around the mean braking pressure P m. First, the part (valve control system) that converts the average brake pressure command and the minute excitation command into actual braking torque for the wheels is, as shown in FIG. 5, a master cylinder 48, a control valve 52, a wheel cylinder 56, and a reservoir 58. And an oil pump 60.
【0051】ブレーキペダル46は、ブレーキペダル4
6の踏力に応じて増圧するマスタシリンダ48を介して
制御バルブ52の増圧バルブ50へ接続されている。ま
た、制御バルブ52は、減圧バルブ54を介して低圧源
としてのリザーバー58へ接続されている。さらに、制
御バルブ52には、該制御バルブによって供給されたブ
レーキ圧をブレーキディスクに加えるためのホイールシ
リンダ56が接続されている。この制御バルブ52は、
入力されたバルブ動作指令に基づいて増圧バルブ50及
び減圧バルブ54の開閉を制御する。The brake pedal 46 includes the brake pedal 4
6 is connected to a pressure increasing valve 50 of a control valve 52 via a master cylinder 48 which increases the pressure in accordance with the pedaling force. The control valve 52 is connected to a reservoir 58 as a low pressure source via a pressure reducing valve 54. Further, a wheel cylinder 56 for applying the brake pressure supplied by the control valve to the brake disc is connected to the control valve 52. This control valve 52 is
The opening and closing of the pressure increasing valve 50 and the pressure reducing valve 54 are controlled based on the input valve operation command.
【0052】なお、この制御バルブ52が増圧バルブ5
0のみを開くように制御されると、ホイールシリンダ5
6の油圧(ホイールシリンダ圧)は、ドライバがブレー
キペダル46を踏み込むことによって得られる圧力に比
例したマスタシリンダ48の油圧(マスタシリンダ圧)
まで上昇する。逆に減圧バルブ54のみを開くように制
御されると、ホイールシリンダ圧は、ほぼ大気圧のリザ
ーバ58の圧力(リザーバ圧)まで減少する。また、両
方のバルブを閉じるように制御されると、ホイールシリ
ンダ圧は保持される。The control valve 52 is connected to the pressure increasing valve 5
0 is controlled to open only the wheel cylinder 5
The hydraulic pressure (wheel cylinder pressure) of the master cylinder 48 (master cylinder pressure) is proportional to the pressure obtained when the driver depresses the brake pedal 46.
To rise. Conversely, when the pressure is controlled so as to open only the pressure reducing valve 54, the wheel cylinder pressure decreases to the pressure of the reservoir 58 (reservoir pressure), which is approximately atmospheric pressure. Further, when both valves are controlled to be closed, the wheel cylinder pressure is maintained.
【0053】ホイールシリンダ56によりブレーキディ
スクに加えられるブレーキ力(ホイールシリンダ圧に相
当)は、マスタシリンダ48の高油圧が供給される増圧
時間、リザーバー58の低油圧が供給される減圧時間、
及び供給油圧が保持される保持時間の比率と、圧力セン
サ等により検出されたマスタシリンダ圧及びリザーバー
圧とから求められる。The braking force (corresponding to the wheel cylinder pressure) applied to the brake disk by the wheel cylinder 56 includes a pressure increasing time during which the high hydraulic pressure of the master cylinder 48 is supplied, a pressure reducing time during which the low hydraulic pressure of the reservoir 58 is supplied,
And the ratio of the holding time during which the supply oil pressure is held, and the master cylinder pressure and the reservoir pressure detected by a pressure sensor or the like.
【0054】従って、制御バルブ52の増減圧時間をマ
スタシリンダ圧に応じて制御することにより、所望のブ
レーキトルクを実現することができる。そして、ブレー
キ圧の微小励振は、平均ブレーキ力を実現する制御バル
ブ52の増減圧制御と同時に共振周波数に対応した周期
で増圧減圧制御を行うことにより可能となる。Therefore, a desired brake torque can be realized by controlling the increasing / decreasing time of the control valve 52 according to the master cylinder pressure. Further, the minute excitation of the brake pressure can be performed by performing the pressure increase / decrease control at a cycle corresponding to the resonance frequency simultaneously with the pressure increase / decrease control of the control valve 52 for realizing the average braking force.
【0055】具体的な制御の内容として、図6に示すよ
うに、微小励振の周期(例えば24[ms])の半周期
T/2毎に増圧と減圧のそれぞれのモードを切り替え、
バルブへの増減圧指令は、モード切り替えの瞬間から増
圧時間ti 、減圧時間tr のそれぞれの時間分だけ増圧
・減圧指令を出力し、残りの時間は、保持指令を出力す
る。平均ブレーキ力は、マスタシリンダ圧に応じた増圧
時間ti と減圧時間tr との比によって定まると共に、
共振周波数に対応した半周期T/2毎の増圧・減圧モー
ドの切り替えによって、平均ブレーキ力の回りに微小振
動が印加される。As specific contents of the control, as shown in FIG. 6, the mode of pressure increase and the mode of pressure decrease are switched every half cycle T / 2 of the period of minute excitation (for example, 24 [ms]).
Decrease pressure command to the valve, the moment between the pressure boosting t i of mode switching, and outputs the respective just time period the pressure increase, pressure reduction command decompression time t r, the rest of the time, and outputs the held command. The average braking force, as well as determined by the ratio of the pressure increasing time t i corresponding to the master cylinder pressure and the decompression time t r,
By switching between the pressure increasing and reducing modes every half cycle T / 2 corresponding to the resonance frequency, a minute vibration is applied around the average braking force.
【0056】なお、ブレーキ圧微小振幅Pv は、マスタ
シリンダ圧、図6に示したバルブの増圧時間ti の長
さ、及び減圧時間tr の長さによって所定の関係で定ま
るので、図2のブレーキ圧微小振幅検出部42は、マス
タシリンダ圧、増圧時間ti 及び減圧時間tr からブレ
ーキ圧微小振幅Pv を出力するテーブルとして構成する
ことができる。[0056] Incidentally, the brake pressure differential small amplitude P v is the master cylinder pressure, the length of the pressure increasing time t i of the valve shown in FIG. 6, and so determined in a predetermined relationship with the length of the decompression time t r, Figure second brake pressure differential small amplitude detector 42 may be a master cylinder pressure, the pressure increasing time t i and the pressure reducing time t r constituting a table for outputting the brake pressure differential small amplitude P v.
【0057】次に、本実施の形態の作用を説明する。最
初に、基準値作成記憶回路12D、12S、12Iによ
る路面μ勾配の基準値の作成方法を説明する。Next, the operation of the present embodiment will be described. First, a method of creating the reference value of the road μ gradient by the reference value creation storage circuits 12D, 12S, and 12I will be described.
【0058】路面μ勾配Gd は、すべり速度ΔVに対す
る路面と車輪との間の摩擦係数μの勾配であり、制動
力、路面状態及び車速毎に、(10)式により表され
る。The road surface μ gradient G d is the gradient of the friction coefficient μ between the road surface and the wheels with respect to the slip speed ΔV, and is expressed by equation (10) for each braking force, road surface condition, and vehicle speed.
【0059】[0059]
【数10】 (Equation 10)
【0060】よって、路面μ勾配Gd は、制動力Pc、
路面状態、車速V、すべり速度ΔV、及び摩擦係数μに
より特定することができる。Accordingly, the road surface μ gradient G d is determined by the braking force Pc,
It can be specified by the road surface condition, the vehicle speed V, the slip speed ΔV, and the friction coefficient μ.
【0061】ところで、すべり速度ΔVと制動力Pと
は、各路面状態に応じて、図7〜図9に示すように、車
速に応じて、固有の関係を有することが分かっている。It is known that the slip speed ΔV and the braking force P have a specific relationship according to the vehicle speed as shown in FIGS. 7 to 9 according to each road surface condition.
【0062】そこで、基準値作成記憶回路12D、12
S、12Iの各すべり速度演算部20は、車速に応じ
て、図7〜図9に示す制動力Pとすべり速度ΔVとの関
係式を記憶している。従って、制動力Pc及び車速Vが
入力されるとすべり速度演算部20は、上記関係式、制
動力Pc及び車速Vに基づいて、制動力Pc及び車速V
に対応するすべり速度ΔVを演算し、演算により求めた
すべり速度ΔVを、μ勾配基準値作成部22に出力す
る。Therefore, the reference value creation storage circuits 12D and 12D
Each of the slip speed calculation units 20 of S and 12I stores a relational expression between the braking force P and the slip speed ΔV shown in FIGS. 7 to 9 according to the vehicle speed. Accordingly, when the braking force Pc and the vehicle speed V are input, the slip speed calculation unit 20 calculates the braking force Pc and the vehicle speed V based on the relational expression, the braking force Pc and the vehicle speed V.
Is calculated, and the slip speed ΔV obtained by the calculation is output to the μ gradient reference value creation unit 22.
【0063】なお、すべり速度演算部20に入力される
車速Vは、車速演算部18により演算したものである。
具体的には、次のように求める。車速演算部18には、
上記すべり速度演算部20により演算されたすべり速度
ΔVと、車輪速センサ28からの車輪速Vωと、が入力
される。すべり速度ΔVは、車速Vから車輪速Vωを減
算した値(V−Vω)である。よって、車速演算部18
は、入力されたすべり速度ΔVと入力された車輪速Vω
とを加算して、車速Vを求め、すべり速度演算部20に
出力する。The vehicle speed V input to the slip speed calculating section 20 is calculated by the vehicle speed calculating section 18.
Specifically, it is obtained as follows. The vehicle speed calculation unit 18
The slip speed ΔV calculated by the slip speed calculator 20 and the wheel speed Vω from the wheel speed sensor 28 are input. The slip speed ΔV is a value obtained by subtracting the wheel speed Vω from the vehicle speed V (V−Vω). Therefore, the vehicle speed calculation unit 18
Is the input slip speed ΔV and the input wheel speed Vω
Are added to obtain the vehicle speed V, which is output to the slip speed calculator 20.
【0064】また、路面μと制動力Pとは、各路面状態
に応じて、図10〜図12に示すように、車速に応じ
て、固有の関係を有することが分かっている。It is known that the road surface μ and the braking force P have a specific relationship according to the vehicle speed as shown in FIGS. 10 to 12 according to each road surface condition.
【0065】そこで、基準値作成記憶回路12D、12
S、12Iの各μ勾配基準値作成部22は、車速に応じ
て、図10〜図12に示す制動力Pと路面μとの関係式
を記憶している。従って、制動力Pc及び車速Vが入力
されるとμ勾配基準値作成部22は、上記関係式、制動
力Pc及び車速Vに基づいて、制動力Pc及び車速Vに
対応する路面μを演算する。Therefore, the reference value creating and storing circuits 12D, 12D
Each μ gradient reference value creation unit 22 of S and 12I stores a relational expression between the braking force P and the road surface μ shown in FIGS. 10 to 12 according to the vehicle speed. Therefore, when the braking force Pc and the vehicle speed V are input, the μ gradient reference value creation unit 22 calculates the road surface μ corresponding to the braking force Pc and the vehicle speed V based on the above relational expression, the braking force Pc and the vehicle speed V. .
【0066】このように、μ勾配基準値作成部22に
は、すべり速度演算部20から、制動力Pc及び車速V
に対応するすべり速度ΔV(2ポイント、ΔV1 、ΔV
2 )が入力され、制動力Pc及び車速Vに対応する路面
μ(2ポイント、μ1 、μ2 )を演算している。そこ
で、μ勾配基準値作成部22は、これらのすべり速度Δ
Vと路面μとを、(10)式に代入することにより、路
面μ勾配Gd の基準値を演算する。そして、車輪速度セ
ンサ28、制動力検出部34からは、制動開始時から所
定時間経過まで所定サンプリング時間毎にサンプリング
された車輪速Vω、制動力Pc が入力され、μ勾配基準
値作成部22は、制動開始時から所定量制動力が経過す
るまで所定量制動力が変化する毎に路面μ勾配Gd の基
準値を演算し、各路面μ勾配Gd の基準値の時系列デー
タを、基準値時系列データ記憶部24に出力する。そし
て、基準値時系列データ記憶部24は、路面μ勾配Gd
の基準値の時系列データを記憶する。As described above, the μ gradient reference value creation unit 22 sends the braking force Pc and the vehicle speed V
Velocity ΔV (2 points, ΔV 1 , ΔV
2 ) is input and the road surface μ (2 points, μ 1 , μ 2 ) corresponding to the braking force Pc and the vehicle speed V is calculated. Therefore, the μ gradient reference value creation unit 22 calculates the slip speed Δ
The V and the road surface mu, by substituting the equation (10) calculates the reference value of the road surface mu gradient G d. Then, from the wheel speed sensor 28 and the braking force detector 34, the wheel speed Vω and the braking force Pc sampled at predetermined sampling times from the start of braking to the elapse of a predetermined time are input, and the μ gradient reference value generator 22 until a predetermined amount a braking force has elapsed from the start of braking calculates the reference value of the road surface μ gradient G d each time a predetermined amount the braking force is changed, the time-series data of the reference value for each road surface μ gradient G d, The reference value is output to the time-series data storage unit 24. The reference value time-series data storage unit 24 stores the road surface μ gradient G d
The time series data of the reference value is stored.
【0067】また、微小ゲイン時系列データ記憶部25
は、μ勾配基準値作成部22から基準値時系列データ記
憶部24への基準値の時系列データの出力タイミングに
同期して、微小ゲインの時系列データが記憶される。即
ち、微小ゲイン時系列データ記憶部25へは、制動開始
時から所定時間経過まで所定サンプリング時間毎にサン
プリングされた微小ゲインGd T、制動力Pc が入力さ
れ、微小ゲイン時系列データ記憶部25は、制動開始時
から前出の所定量制動力が経過するまで所定量制動力が
変化する毎に、微小ゲインGd Tを記憶する。Further, the minute gain time series data storage unit 25
The time series data of the minute gain is stored in synchronization with the output timing of the time series data of the reference value from the μ gradient reference value creation unit 22 to the reference value time series data storage unit 24. That is, the minute gain G d T and the braking force P c sampled at predetermined sampling times from the start of braking to the elapse of a predetermined time from the start of braking are input to the minute gain time series data storage unit 25. The reference numeral 25 stores a small gain G d T every time the predetermined amount of braking force changes from the start of braking until the aforementioned predetermined amount of braking force has elapsed.
【0068】次に、制動開始時に、時系列データ軌跡照
合部14が実行する照合処理ルーチンを図13を参照し
て説明する。Next, a comparison processing routine executed by the time-series data trajectory comparison unit 14 at the start of braking will be described with reference to FIG.
【0069】図13のステップ82で、制動開始時から
制動力が所定量変化した時刻を表す変数kを初期化し、
ステップ84で、変数kを1インクリメントする。In step 82 of FIG. 13, a variable k representing the time at which the braking force has changed by a predetermined amount from the start of braking is initialized,
In step 84, the variable k is incremented by one.
【0070】ステップ86で、変数kにより表される時
刻(k・τ,τ:サンプリング周期)において微小ゲイ
ン演算部36により演算された路面μ勾配Gd T(k)
(微小ゲイン)を取り込む。At step 86, the road surface μ gradient G d T (k) calculated by the minute gain calculator 36 at the time (k · τ, τ: sampling period) represented by the variable k.
(Small gain).
【0071】ステップ88で、路面状態を識別する変数
sを初期化し、ステップ90で、変数sで、変数kを1
インクリメントする。At step 88, a variable s for identifying the road surface condition is initialized, and at step 90, the variable k is set to 1 by the variable s.
Increment.
【0072】ステップ92で、路面状態sに対応する基
準値時系列データ24に記憶された路面μ勾配Gd の時
系列データの内、時刻kにおける基準値Gd s(k)を
取り込む。[0072] In step 92, the time of the series data of the road surface μ gradient G d stored in the reference value time series data 24 corresponding to the road surface condition s, captures the reference value G d s at time k (k).
【0073】ステップ94で、照合値zs を(11)式
より演算する。In step 94, the collation value z s is calculated from equation (11).
【0074】[0074]
【数11】 zs =zs +(Gd T(k)−Gd s(k))2 ・・・(11) ステップ96で、変数sが予め定められた路面状態の総
数s0 (本実施の形態では、3)以上か否かを判断す
る。変数s≧総数s0 でない場合には、時刻kにおける
照合値zs を演算していない路面状態があるので、ステ
ップ90に戻って、以上の処理(ステップ90〜ステッ
プ96)を実行する。変数s≧総数s0 の場合には、全
ての路面状態の時刻kにおける照合値zs を演算したの
で、ステップ98で、変数kが、制動開始から所定時間
経過した時刻を表すn以上か否かを判断する。変数k≧
nでない場合には、制動開始から所定時間経過していな
いので、ステップ84に戻って以上の処理(ステップ8
4〜ステップ98)を実行する。変数k≧nの場合に
は、制動開始から所定時間経過したので、ステップ10
0で、照合値zs (z1 〜z3 )の内の最小の照合値z
s に対応する路面状態が、現在走行している路面状態で
あると推定する。Equation 11] z s = z s + (G d T (k) -G d s (k)) 2 ··· (11) at step 96, the total number s 0 of the road surface state variable s is predetermined ( In the present embodiment, it is determined whether 3) or more. If not the variable s ≧ total number s 0, since there is a road surface condition that is not computed verification value z s at time k, returns to step 90 to execute the above processing (step 90 step 96). When the variable s ≧ the total number s 0 , the collation values z s at all times k of the road surface state have been calculated. Judge. Variable k ≧
If it is not n, since the predetermined time has not elapsed since the start of braking, the flow returns to step 84 and the above processing (step 8)
4 to step 98) are executed. If the variable k ≧ n, since a predetermined time has elapsed since the start of braking, step 10
0, the minimum collation value z of the collation values z s (z 1 to z 3 )
It is estimated that the road surface state corresponding to s is the current road surface state.
【0075】即ち、zs は(12)式により表される。That is, z s is represented by equation (12).
【0076】[0076]
【数12】 (Equation 12)
【0077】このように、zs を、制動開始から所定時
間経過までの時系列データにより演算しているのは、図
14に示すように、微小ゲイン演算部36により演算さ
れた路面μ勾配Gd T(k)は、現在走行している路面
状態の路面μ勾配Gd s(k)(図14では、Snow
の路面状態である)と必ずしも一致するとは限らず、現
在走行している路面状態の路面μ勾配Gd sの基準値ま
わりにふらつく。よって、制動開始から所定量制動力が
経過するまでの時系列データを用いて精度よく路面状態
を推定するためである。As described above, z s is calculated based on the time-series data from the start of braking to the lapse of a predetermined time because the road surface μ gradient G calculated by the minute gain calculator 36 is calculated as shown in FIG. d T (k) is the currently running to the road surface condition and the road surface μ gradient G d s (k) (FIG. 14, Snow
Does not always coincide with the road surface state), and fluctuates around the reference value of the road surface μ gradient G ds of the road surface state currently running. Therefore, the road surface condition is accurately estimated using the time-series data from when the braking starts to when the predetermined amount of braking force has elapsed.
【0078】そして、各路面状態に対応して予め記憶し
ている摩擦係数のピーク値μMAX の内から、推定した路
面状態に対応する摩擦係数のピーク値μMAX を、図示し
ないABS制御部等に出力し、推定した路面状態を表す
データをすべり速度選択部16に出力する。Then, from the friction coefficient peak value μ MAX stored in advance corresponding to each road surface condition, the estimated friction coefficient peak value μ MAX corresponding to the road surface condition is converted to an ABS control unit (not shown) or the like. And outputs data representing the estimated road surface condition to the slip speed selection unit 16.
【0079】時系列データ軌跡照合部14からピーク値
μMAX を入力したABS制御部は、入力したピーク値μ
MAX を与えるスリップ率を基準スリップ率とし、ピーク
μに追従するようにブレーキ力を制御する。The ABS control unit which receives the peak value μ MAX from the time-series data trajectory collation unit 14 inputs the peak value μ MAX.
The slip ratio giving MAX is set as a reference slip ratio, and the braking force is controlled so as to follow the peak μ.
【0080】時系列データ軌跡照合部14から路面状態
を表すデータを入力したすべり速度選択部16は、基準
値作成記憶回路12D、12S、12Iの各すべり速度
演算部20から出力されたすべり速度の内、入力した路
面状態を表すデータに基づいて、推定された路面状態に
対応するすべり速度ΔVを出力する。The slip speed selecting unit 16 which has received the data representing the road surface state from the time-series data trajectory collating unit 14 calculates the slip speed output from the slip speed calculating unit 20 of the reference value creation storage circuits 12D, 12S, 12I. The slip speed ΔV corresponding to the estimated road surface state is output based on the input data representing the road surface state.
【0081】以上のように本実施の形態では、現在走行
中の路面の路面状態やピークμ値を正確に求めることが
できるので、VSC、ABS、TRC等の車両の安定化
制御(制御ゲインの変更、制御目標値の設定))やドラ
イバへの路面状態の警告、車体横すべり角、ヨーレート
等の各車両状態の推定することが可能となる。As described above, according to the present embodiment, the road surface condition and the peak μ value of the road surface currently running can be accurately obtained, and therefore, the vehicle stabilization control (such as control gain of VSC, ABS, TRC, etc.) Change, setting of a control target value), warning of a road surface condition to a driver, and estimation of various vehicle conditions such as a vehicle body slip angle and a yaw rate.
【0082】以上説明した実施の形態では、図7〜図9
に示す制動力Pとすべり速度ΔVとの関係式、図10〜
図12に示す制動力Pと路面μとの関係式を記憶してい
るが、本発明はこれに限定されず、制動力Pc と路面μ
勾配とは図15に示す関係を有するので、図15に示す
関係式を、車速Vに対応して記憶し、車速V、制動力P
c 、及び制動力Pc と路面μ勾配との関係式から、車速
V及び制動力Pc に対応する路面μ勾配を演算するよう
にしてもよい。In the embodiment described above, FIGS.
The relational expression between the braking force P and the slip velocity ΔV shown in FIG.
Although the relational expression between the braking force P and the road surface μ shown in FIG. 12 is stored, the present invention is not limited to this, and the braking force Pc and the road surface μ are stored.
Since the gradient has the relationship shown in FIG. 15, the relational expression shown in FIG. 15 is stored in correspondence with the vehicle speed V, and the vehicle speed V and the braking force P are stored.
The road surface μ gradient corresponding to the vehicle speed V and the braking force Pc may be calculated from c and the relational expression between the braking force Pc and the road surface μ gradient.
【0083】また、前述した実施の形態では、路面μ勾
配Gd T(k)と基準値Gd s(k)との二乗法の最小
値により路面状態を推定しているが、本発明はこれに限
定されず、制動開始から所定時間経過までの路面μ勾配
Gd T(k)及び基準値Gds(k)の時系列データの
相関関係を演算して、路面状態を推定するようにしても
い。[0083] Further, in the embodiment described above, although estimates the road surface condition by a minimum value of the squares of the road surface μ gradient G d T (k) with a reference value G d s (k), the present invention is not limited thereto, and calculates the correlation between the time-series data of the road from the start of braking until a predetermined time elapses μ gradient G d T (k) and the reference value G d s (k), to estimate the road surface condition Anyway.
【0084】更に、前述した例では、図7〜図9に示す
制動力Pとすべり速度ΔVとの関係式、図10〜図12
に示す制動力Pと路面μとの関係式、又は、図15に示
す制動力Pと路面μ勾配Gd との関係式を記憶するよう
にしているが、本発明はこれに限定されず、すべり速度
ΔVと、すべり速度ΔV、路面μ、路面μ勾配Gd との
関係の各々も、制動力Pと、すべり速度ΔV、路面μ、
路面μ勾配Gd との関係の各々と同様な関係となるの
で、すべり速度ΔVと、すべり速度ΔV、路面μ、路面
μ勾配Gd との関係式を記憶し、同様に処理するように
してもよい。Further, in the above-described example, the relational expression between the braking force P and the slip speed ΔV shown in FIGS.
Relationship between the braking force P and the road surface μ shown in, or has provided to store the relationship between the braking force P and the road surface μ gradient G d shown in FIG. 15, the present invention is not limited thereto, and the sliding velocity [Delta] V, sliding velocity [Delta] V, the road surface mu, also each of the relationship between the road surface mu gradient G d, braking force and P, sliding velocity [Delta] V, the road surface mu,
Since the same as the respective relations of the relationship between the road surface mu gradient G d, and the sliding velocity [Delta] V, sliding velocity [Delta] V, the road surface mu, stores a relational expression between the road surface mu gradient G d, so as to similarly treated Is also good.
【0085】また、前述した例では、路面μ勾配Gd と
して、微小ゲインを求めているが、本発明はこれに限定
されず、微小ゲインが制動トルク勾配と等価な物理量で
あるので、微小ゲイン演算部36に代えて制動トルク勾
配演算部を備え、路面μ勾配Gd に代えて、制動トルク
勾配を求めて、同様に処理するようにしてもよい。以
下、微小ゲインが、制動トルク勾配と等価な物理量であ
ることを説明する。[0085] In the example described above, as the road surface μ gradient G d, but seeking small gains, the present invention is not limited thereto, since small gain is braking torque gradient equivalent physical quantity, small gain instead of the calculation unit 36 includes a braking torque gradient calculation unit, instead of the road surface μ gradient G d, seeking braking torque gradient may be treated similarly. Hereinafter, it will be described that the small gain is a physical quantity equivalent to the braking torque gradient.
【0086】図20、図16の特性では、スリップ速度
ΔVと、車輪−路面間の摩擦係数μとの間には、あるス
リップ速度で摩擦係数μがピークをとる関数関係が成立
する。According to the characteristics shown in FIGS. 20 and 16, a functional relationship is established between the slip speed ΔV and the friction coefficient μ between the wheel and the road surface, in which the friction coefficient μ has a peak at a certain slip speed.
【0087】ところで、ブレーキ圧を微小励振すると、
車輪速度が微小励振するので、スリップ速度もあるスリ
ップ速度の回りで微小振動する。ここで、図16の特性
を有する路面において、あるスリップ速度の回りで微小
振動したときの摩擦係数μのスリップ速度ΔVに対する
変化を考える。By the way, when the brake pressure is slightly excited,
Since the wheel speed slightly excites, the slip speed also slightly vibrates around a certain slip speed. Here, consider a change in the friction coefficient μ with respect to the slip speed ΔV when the vehicle slightly vibrates around a certain slip speed on a road surface having the characteristics shown in FIG.
【0088】このとき、路面の摩擦係数μは、At this time, the friction coefficient μ of the road surface is
【0089】[0089]
【数13】μ = μ0 +αRΔω・・・(13) と近似できる。すなわち、微小振動によるスリップ速度
の変化が小さいため、傾きαRの直線で近似できる。Μ = μ 0 + αRΔω (13) That is, since the change in the slip speed due to the minute vibration is small, it can be approximated by a straight line having the slope αR.
【0090】ここで、タイヤと路面間の摩擦係数μによ
り生じる制動トルクTb =μWRに(13)式を代入す
ると、Here, when the equation (13) is substituted into the braking torque T b = μWR generated by the friction coefficient μ between the tire and the road surface,
【0091】[0091]
【数14】 Tb = μWR = μ0 WR+αR2 ΔωW・・・(14) となる。ここで、Wは輪荷重である。(14)式の両辺
をΔωで1階微分すると、The [number 14] T b = μWR = μ 0 WR + αR 2 ΔωW ··· (14). Here, W is a wheel load. First-order differentiation of both sides of equation (14) with Δω gives:
【0092】[0092]
【数15】 (Equation 15)
【0093】を得る。よって、(15)式により、制動
トルク勾配(dTb /Δω)が、αR2 Wに等しいこと
が示された。Is obtained. Therefore, Equation (15) shows that the braking torque gradient (dT b / Δω) is equal to αR 2 W.
【0094】一方、ブレーキトルクTb ’がブレーキ圧
Pb と比例関係にあることから、微小ゲインGd は、ブ
レーキトルクTb ’に対する車輪速度ωw の比(ωw /
Tb’)の共振周波数ω∞の振動成分と比例関係にあ
る。従って、(5) 式の伝達特性により、微小ゲインGd
は次式によって表される。On the other hand, since the brake torque T b ′ is proportional to the brake pressure P b , the small gain G d is determined by the ratio of the wheel speed ω w to the brake torque T b ′ (ω w /
T b ′) is proportional to the vibration component of the resonance frequency ω∞. Therefore, the small gain G d is obtained by the transfer characteristic of equation (5).
Is represented by the following equation.
【0095】[0095]
【数16】 (Equation 16)
【0096】但し、However,
【0097】[0097]
【数17】 JA =Jt +Jv +Jw 、JB =Jt +Jv ・・・(17)J A = J t + J v + J w , J B = J t + J v (17)
【0098】[0098]
【数18】 (Equation 18)
【0099】一般に、(18)式において、In general, in equation (18),
【0100】[0100]
【数19】 |A| = 0.012 << |B| = 0.1・・・(19) となることから、(15)、(16)式より、| A | = 0.012 << | B | = 0.1 (19) From the equations (15) and (16),
【0101】[0101]
【数20】 (Equation 20)
【0102】を得る。すなわち、スリップ速度ΔVに対
する制動トルクTb の勾配(制動トルク勾配)は微小ゲ
インGd に比例する。Is obtained. That is, the gradient (the braking torque gradient) of the braking torque T b with respect to the slip speed ΔV is proportional to the micro-gain G d.
【0103】よって、制動トルク勾配を求め、求めた制
動トルク勾配に基づいて、上記と同様に処理すればよ
い。Therefore, the braking torque gradient may be determined, and the same processing as described above may be performed based on the determined braking torque gradient.
【0104】次に、制動トルク勾配を求める方法を説明
する。各車輪の車輪運動及び車体運動は、(21)式、
(22)式の運動方程式によって記述される。Next, a method for obtaining the braking torque gradient will be described. The wheel motion and the body motion of each wheel are expressed by equation (21),
It is described by the equation of motion of equation (22).
【0105】[0105]
【数21】 (Equation 21)
【0106】[0106]
【数22】 ただし、Fi ’は、第i輪に発生した制動力、Tbiは踏
力に対応して第i輪に加えられたブレーキトルク、Mは
車両質量、Rc は車輪の有効半径、Jは車輪慣性、vは
車体速度である(図8参照)。なお、・は時間に関する
微分を示す。(21)式、(22)式において、Fi ’
はスリップ速度(v/Rc −ωi )の関数として示され
ている。(Equation 22) Here, F i ′ is the braking force generated on the i-th wheel, T bi is the braking torque applied to the i-th wheel corresponding to the pedaling force, M is the vehicle mass, R c is the effective radius of the wheel, and J is the wheel The inertia, v, is the vehicle speed (see FIG. 8). In addition, * shows differentiation with respect to time. In equations (21) and (22), F i ′
Is shown as a function of the slip speed (v / R c -ω i ).
【0107】ここで、車体速度を等価的な車体の角速度
ωv で表す((23)式)と共に、制動トルクR
c Fi ’をスリップ速度の1次関数(傾きki 、y切片
Ti )として記述する((24)式)。[0107] Here, the vehicle speed represented by an equivalent body of the angular velocity omega v with ((23)), the braking torque R
c F i ′ is described as a linear function (slope k i , y-intercept T i ) of the slip speed (Equation (24)).
【0108】[0108]
【数23】v = Rc ωv ・・・(23)[Number 23] v = R c ω v ··· ( 23)
【0109】[0109]
【数24】 Rc Fi ’(ωv −ωi )=ki ×(ωv −ωi )+Ti ・・・(24) さらに、(23)式、(24)式を、(21)式、(2
2)式へ代入し、車輪速度ωi 及び車体速度ωv をサン
プル時間τ毎に離散化された時系列データωi[k] 、ω
v [k] (kはサンプル時間τを単位とするサンプル時
刻、k=1,2,.....)として表すと、(25)式、(2
6)式を得る。Equation 24] R c F i '(ω v -ω i) = k i × (ω v -ω i) + T i ··· (24) Further, equation (23), the equation (24), (21 ) Formula, (2
2), the wheel speed ω i and the vehicle speed ω v are discretized for each sample time τ, and the time series data ω i [k], ω
v [k] (k is the sample time in units of the sample time τ, k = 1, 2,...)
6) Obtain the equation.
【0110】[0110]
【数25】 (Equation 25)
【0111】[0111]
【数26】 (Equation 26)
【0112】ここで、(25)式、(26)式を連立
し、車体の等価角速度ωv を消去すると、Here, when equations (25) and (26) are combined to eliminate the equivalent angular velocity ω v of the vehicle body,
【0113】[0113]
【数27】 [Equation 27]
【0114】を得る。ところで、スリップ速度3rad/s
という条件下でRc Mg/4(gは重力加速度)の最大
制動トルクの発生を仮定すると、Is obtained. By the way, slip speed 3rad / s
Assuming the generation of the maximum braking torque of R c Mg / 4 (g is the gravitational acceleration) under the condition
【0115】[0115]
【数28】 [Equation 28]
【0116】を得る。ここで、具体的な定数として、τ
=0.01 (sec)、Rc =0.3 (m) 、M=1000(kg)を考慮す
ると、Is obtained. Here, as a specific constant, τ
= 0.01 (sec), R c = 0.3 (m) and M = 1000 (kg)
【0117】[0117]
【数29】 (Equation 29)
【0118】となり、(27)式は次式のように近似す
ることができる。Equation (27) can be approximated as follows.
【0119】[0119]
【数30】 [Equation 30]
【0120】ただし、However,
【0121】[0121]
【数31】 (Equation 31)
【0122】である。(31)式中のTb について、ブ
レーキ油圧(マスタシリンダ圧又はホイールシリンダ
圧)が計測できる場合、次のようにおくことができる。Is as follows. (31) for T b in the formula, the brake hydraulic pressure (master cylinder pressure or wheel cylinder pressure) If can be measured, can be placed as follows.
【0123】[0123]
【数32】Tb =kb ・Pc (k)・・・(32) 但し、Pc はブレーキ油圧、kb はブレーキ油圧からブ
レーキトルクへの変換係数である。この場合、(30)
式、(31)式は、(33)式、(34)式に変換され
る。T b = k b · P c (k) (32) where P c is a brake oil pressure, and k b is a conversion coefficient from the brake oil pressure to the brake torque. In this case, (30)
Equations (31) are converted to Equations (33) and (34).
【0124】[0124]
【数33】 [Equation 33]
【0125】[0125]
【数34】 (Equation 34)
【0126】但し、However,
【0127】[0127]
【数35】 (Equation 35)
【0128】ここで、現サンプリング時点kにおけるブ
レーキ油圧Pc (k)と、1サンプリング前の時点k−
1におけるブレーキ油圧Pc (k−1)とは、ほとんど
等しいとする(Pc (k)=Pc (k−1))。Here, the brake oil pressure P c (k) at the current sampling time k and the time k−
It is assumed that the brake hydraulic pressure P c (k-1) at 1 is almost equal (P c (k) = P c (k-1)).
【0129】また、(35)式中のμ勾配kj は、前回
の推定値を使用する。このように整理することにより、
(33)式は未知係数ki 、fi に関し、線形の形で記
述することが可能となり、(33)式にオンラインのパ
ラメータ同定手法を適用することにより、スリップ速度
に対する制動トルク勾配ki を推定することができる。As the μ gradient k j in the equation (35), the previous estimated value is used. By organizing in this way,
The equation (33) can describe the unknown coefficients k i and f i in a linear form. By applying the online parameter identification method to the equation (33), the braking torque gradient k i with respect to the slip speed can be calculated. Can be estimated.
【0130】すなわち、以下のステップ1及びステップ
2を繰り返すことにより、検出された車輪速度の時系列
データωi [k] から制動トルク勾配の時系列データを推
定することができる(最小自乗推定法)。That is, by repeating the following steps 1 and 2, it is possible to estimate the braking torque gradient time series data from the detected wheel speed time series data ω i [k] (least squares estimation method). ).
【0131】ステップ1:Step 1:
【0132】[0132]
【数36】φi [k] T ・θi =yi [k] ・・・(36) 但し、Φ i [k] T · θ i = y i [k] (36) where
【0133】[0133]
【数37】 (37)
【0134】[0134]
【数38】 (38)
【0135】[0135]
【数39】 yi [k] =−ωi [k] + 2ωi [k−1]−ωi [k−2]−Pci(k)・・・(39) とおく。なお、(37)式の行列φi [k] の第1要素
は、1サンプル時間での車輪速度の変化に関する物理量
であり、(39)式は、1サンプル時間の車輪速度の変
化の1サンプル時間での変化に関する物理量である。Y i [k] = − ω i [k] + 2ω i [k−1] −ω i [k−2] −P ci (k) (39) The first element of the matrix φ i [k] in the equation (37) is a physical quantity related to the change in the wheel speed in one sample time, and the equation (39) is one sample in the change in the wheel speed in one sample time. It is a physical quantity related to change with time.
【0136】ステップ2:Step 2:
【0137】[0137]
【数40】 (Equation 40)
【0138】[0138]
【数41】 [Equation 41]
【0139】[0139]
【数42】 (Equation 42)
【0140】という漸化式からθi の推定値を演算し、
θi の推定値の行列の第一要素を推定された制動トルク
の勾配として抽出する。ただし、λは過去のデータを取
り除く度合いを示す忘却係数(例えばλ=0.98)で
あり、”T ”は行列の転置を示す。The estimated value of θ i is calculated from the recurrence formula
The first element of the matrix of the estimated values of θ i is extracted as the gradient of the estimated braking torque. Here, λ is a forgetting factor (for example, λ = 0.98) indicating a degree of removing past data, and “ T ” indicates transposition of a matrix.
【0141】なお、(40)式の左辺は、車輪速度の変
化に関する物理量の履歴及び車輪速度の変化の変化に関
する物理量の履歴を表す物理量である。The left side of the equation (40) is a physical quantity representing the history of the physical quantity related to the change in the wheel speed and the history of the physical quantity related to the change in the change in the wheel speed.
【0142】図17〜図19に示すように、Dry、S
now、Ice の各路面において、(31)式から求
めた制動トルク勾配に変換係数を乗算して求めた路面μ
勾配Gd の時系列データGd T1 は、路面μ勾配Gd の
基準値の時系列データGdsに精度よく一致し、より実走
行路面の軌跡に近づいており、路面状態の推定精度が向
上する。なお、図17〜図19には、(29)式から求
めた制動トルク勾配に変換係数を乗算して求めた路面μ
勾配Gd の時系列データGd T2 を示している。時系列
データGd T1 は、時系列データGd T2 より、時系列
データGdsに精度よく一致している。よって、時系列デ
ータGd T1 は、時系列データGd T2より改善されて
いる。As shown in FIGS. 17 to 19, Dry, S
On each of the road surfaces now and Ice, the road surface μ obtained by multiplying the braking torque gradient calculated from the equation (31) by a conversion coefficient is calculated.
Time series data G d T 1 gradient G d is precisely matched to the series data G ds at the reference value of the road surface μ gradient G d, is approaching more the curve of the real road surface, the estimation accuracy of the road surface condition improves. 17 to 19 show the road surface μ obtained by multiplying the braking torque gradient obtained from the equation (29) by a conversion coefficient.
The time series data G d T 2 of the gradient G d is shown. The time series data G d T 1 matches the time series data G ds more accurately than the time series data G d T 2 . Therefore, the time series data G d T 1 is improved from the time series data G d T 2 .
【0143】[0143]
【発明の効果】以上説明したように本発明は、第2の物
理量を検出し、検出された第2の物理量の軌跡と、複数
の路面状態各々に対応して演算された第2の物理量の軌
跡の各々と、を比較するので、比較結果に基づいて、複
数の路面状態各々に対応して演算された軌跡の内の検出
された軌跡に最も近い軌跡を選択すれば、選択した軌跡
に対応する路面状態は車輪が走行している実際の路面状
態に対応し、実際の路面状態を精度よく推定することが
できる、という効果を有する。As described above, according to the present invention, the second physical quantity is detected, the trajectory of the detected second physical quantity, and the second physical quantity of the second physical quantity calculated corresponding to each of a plurality of road surface conditions are detected. Since each of the trajectories is compared with each other, if the trajectory closest to the detected trajectory is selected from the trajectories calculated for each of the plurality of road surface conditions based on the comparison result, the trajectory corresponds to the selected trajectory. This road surface state corresponds to the actual road surface state where the wheels are traveling, and has the effect that the actual road surface state can be accurately estimated.
【図1】本実施の形態に係る路面状態推定装置のブロッ
ク図である。FIG. 1 is a block diagram of a road surface state estimation device according to the present embodiment.
【図2】微小ゲイン演算部のブロック図である。FIG. 2 is a block diagram of a minute gain calculation unit.
【図3】車輪と車体と路面とから構成される振動系の等
価モデルを示す図である。FIG. 3 is a diagram showing an equivalent model of a vibration system composed of wheels, a vehicle body, and a road surface.
【図4】車輪速微小振幅検出部のブロック図である。FIG. 4 is a block diagram of a wheel speed minute amplitude detection unit.
【図5】ブレーキ圧微小振幅検出部のブロック図であ
る。FIG. 5 is a block diagram of a brake pressure minute amplitude detection unit.
【図6】ブレーキ圧の微小励振と平均ブレーキ力の制御
を同時に行う場合の制御バルブへの指令を示す図であ
る。FIG. 6 is a diagram showing a command to a control valve in a case where micro excitation of a brake pressure and control of an average braking force are simultaneously performed.
【図7】Dry路面状態における所定速度のときの制動
力とすべり速度との関係を示す線図である。FIG. 7 is a diagram showing a relationship between a braking force and a slip speed at a predetermined speed in a Dry road surface state.
【図8】Snow路面状態における所定速度のときの制
動力とすべり速度との関係を示す線図である。FIG. 8 is a diagram showing a relationship between a braking force and a slip speed at a predetermined speed in a Snow road surface state.
【図9】Ice路面状態における所定速度のときの制動
力とすべり速度との関係を示す線図である。FIG. 9 is a diagram illustrating a relationship between a braking force and a slip speed at a predetermined speed in an Ice road surface state.
【図10】Dry路面状態における所定速度のときの制
動力と路面μとの関係を示す線図である。FIG. 10 is a diagram showing a relationship between a braking force and a road surface μ at a predetermined speed in a Dry road surface state.
【図11】Snow路面状態における所定速度のときの
制動力と路面μとの関係を示す線図である。FIG. 11 is a diagram showing a relationship between a braking force and a road surface μ at a predetermined speed in a Snow road surface state.
【図12】Ice路面状態における所定速度のときの制
動力と路面μとの関係を示す線図である。FIG. 12 is a diagram showing a relationship between a braking force and a road surface μ at a predetermined speed in an Ice road surface state.
【図13】時系列データ軌跡照合部が、制動開始時に実
行する照合処理ルーチンを示したフローチャートであ
る。FIG. 13 is a flowchart illustrating a collation processing routine executed by the time-series data trajectory collation unit at the start of braking.
【図14】時系列データ軌跡照合部の照合処理(軌跡)
の様子を示した図である。FIG. 14 is a collation process (trajectory) of the time-series data trajectory collation unit;
FIG.
【図15】制動力と路面μ勾配との関係を示す線図であ
る。FIG. 15 is a diagram showing a relationship between a braking force and a road surface μ gradient.
【図16】スリップ速度に対する摩擦係数μの変化特性
を示すと共に、微小ゲインが制動トルク勾配と等価であ
ることを説明するため、微小振動の中心の回りのμの変
化が直線で近似できることを示す図である。FIG. 16 shows a change characteristic of the friction coefficient μ with respect to the slip speed, and shows that the change of μ around the center of the minute vibration can be approximated by a straight line in order to explain that the minute gain is equivalent to the braking torque gradient. FIG.
【図17】Dry路面状態における所定速度のときの制
動力と路面μ勾配との関係を示す線図である。FIG. 17 is a diagram showing a relationship between a braking force and a road surface μ gradient at a predetermined speed in a Dry road surface state.
【図18】Snow路面状態における所定速度のときの
制動力と路面μ勾配との関係を示す線図である。FIG. 18 is a diagram showing a relationship between a braking force and a road μ gradient at a predetermined speed in a Snow road surface state.
【図19】Ice路面状態における所定速度のときの制
動力と路面μ勾配との関係を示す線図である。FIG. 19 is a diagram showing a relationship between a braking force and a road μ gradient at a predetermined speed in an Ice road surface state.
【図20】タイヤと路面との間の摩擦係数μのスリップ
率sに対する特性を示す線図である。FIG. 20 is a diagram showing characteristics of a friction coefficient μ between a tire and a road surface with respect to a slip ratio s.
10 演算回路 12D、12S、12I 基準値作成記憶回路 14 時系列データ軌跡照合部 28 車輪速センサ 34 制動力検出部 36 微小ゲイン演算部 DESCRIPTION OF SYMBOLS 10 Operation circuit 12D, 12S, 12I Reference value creation storage circuit 14 Time-series data trajectory comparison unit 28 Wheel speed sensor 34 Braking force detection unit 36 Micro gain calculation unit
───────────────────────────────────────────────────── フロントページの続き (72)発明者 小野 英一 愛知県愛知郡長久手町大字長湫字横道41番 地の1株式会社豊田中央研究所内 (72)発明者 梅野 孝治 愛知県愛知郡長久手町大字長湫字横道41番 地の1株式会社豊田中央研究所内 (72)発明者 菅井 賢 愛知県愛知郡長久手町大字長湫字横道41番 地の1株式会社豊田中央研究所内 ──────────────────────────────────────────────────続 き Continued on the front page (72) Inventor Eiichi Ono 41-Cho, Yokomichi, Nagakute-cho, Aichi-gun, Aichi Prefecture Inside Toyota Central R & D Laboratories Co., Ltd. (41) Inventor Satoshi Sugai, Toyota Chuo R & D Center, 41, Nagakute-cho, Aichi-gun, Aichi-gun, Japan
Claims (1)
量を検出する第1の検出手段と、 前記車両の速度を検出する第2の検出手段と、 前記車輪のすべり易さを表しかつ前記車輪が走行する路
面状態、前記第1の物理量、及び車両の速度に応じて定
まる第2の物理量を検出する第3の検出手段と、 前記第1の検出手段により検出された第1の物理量及び
前記第2の検出手段により検出された車速に基づいて、
予め定められた複数の路面状態毎に第2の物理量を演算
する演算手段と、 前記車輪の制動状態を表す第1の物理量に対する前記第
3の検出手段により検出された第2の物理量の軌跡と、
前記第1の物理量に対する前記演算手段により複数の路
面状態毎に演算された第2の物理量の軌跡の各々と、を
比較する比較手段と、 前記比較手段の比較結果に基づいて、前記演算された軌
跡の内の前記検出された軌跡に最も近い軌跡に対応する
路面状態を、前記車輪が走行している実際の路面状態と
して推定する推定手段と、 を備えた路面状態推定装置。1. A first detecting means for detecting a first physical quantity representing a braking state of a wheel of a vehicle; a second detecting means for detecting a speed of the vehicle; Third detection means for detecting a road surface state on which the wheel travels, the first physical quantity, and a second physical quantity determined according to the speed of the vehicle; a first physical quantity detected by the first detection means And based on the vehicle speed detected by the second detection means,
Calculating means for calculating a second physical quantity for each of a plurality of predetermined road surface conditions; and a locus of a second physical quantity detected by the third detecting means with respect to the first physical quantity representing the braking state of the wheel. ,
Comparing means for comparing each of the trajectories of the second physical quantity calculated for each of the plurality of road surface states by the calculating means with respect to the first physical quantity, and the calculation based on the comparison result of the comparing means A road surface state estimating device comprising: estimating means for estimating a road surface state corresponding to a locus closest to the detected locus in the locus as an actual road surface state on which the wheels are traveling.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP32692197A JP3424535B2 (en) | 1997-11-11 | 1997-11-11 | Road surface condition estimation device |
EP98111728A EP0887241A3 (en) | 1997-06-27 | 1998-06-25 | Road surface condition estimating apparatus and variation reduction processing apparatus |
US09/106,342 US6324461B1 (en) | 1997-06-27 | 1998-06-29 | Road surface condition estimating apparatus and variation reduction processing apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP32692197A JP3424535B2 (en) | 1997-11-11 | 1997-11-11 | Road surface condition estimation device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH11139293A true JPH11139293A (en) | 1999-05-25 |
JP3424535B2 JP3424535B2 (en) | 2003-07-07 |
Family
ID=18193255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP32692197A Expired - Fee Related JP3424535B2 (en) | 1997-06-27 | 1997-11-11 | Road surface condition estimation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3424535B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1215096A3 (en) * | 2000-12-14 | 2003-06-04 | Sumitomo Rubber Industries, Ltd. | Apparatus and method for identifying tires and apparatus and method for evaluating road surface conditions |
EP1479580A1 (en) * | 2003-05-22 | 2004-11-24 | Fuji Jukogyo Kabushiki Kaisha | Device for estimating friction coefficient on road surface of vehicle |
US7483808B2 (en) * | 2007-06-29 | 2009-01-27 | Caterpillar Inc. | System and method for measuring machine rolling resistance |
CN113299060A (en) * | 2021-04-29 | 2021-08-24 | 山西省交通科技研发有限公司 | Vehicle information coefficient acquisition system based on measured data |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62255284A (en) * | 1986-04-30 | 1987-11-07 | Mazda Motor Corp | Running road state discriminator |
JPH01249559A (en) * | 1988-03-31 | 1989-10-04 | Nissan Motor Co Ltd | Antiskid control device |
JPH04224447A (en) * | 1990-04-02 | 1992-08-13 | Volkswagen Ag <Vw> | Method of monitoring frictional coupling between tire of wheel driven and road |
JPH07112659A (en) * | 1993-10-18 | 1995-05-02 | Nippondenso Co Ltd | Road surface friction coefficient estimating device |
JPH07186928A (en) * | 1993-12-24 | 1995-07-25 | Aisin Seiki Co Ltd | Estimation device for maximum coefficient of friction of vehicle traveling road surface |
-
1997
- 1997-11-11 JP JP32692197A patent/JP3424535B2/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62255284A (en) * | 1986-04-30 | 1987-11-07 | Mazda Motor Corp | Running road state discriminator |
JPH01249559A (en) * | 1988-03-31 | 1989-10-04 | Nissan Motor Co Ltd | Antiskid control device |
JPH04224447A (en) * | 1990-04-02 | 1992-08-13 | Volkswagen Ag <Vw> | Method of monitoring frictional coupling between tire of wheel driven and road |
JPH07112659A (en) * | 1993-10-18 | 1995-05-02 | Nippondenso Co Ltd | Road surface friction coefficient estimating device |
JPH07186928A (en) * | 1993-12-24 | 1995-07-25 | Aisin Seiki Co Ltd | Estimation device for maximum coefficient of friction of vehicle traveling road surface |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1215096A3 (en) * | 2000-12-14 | 2003-06-04 | Sumitomo Rubber Industries, Ltd. | Apparatus and method for identifying tires and apparatus and method for evaluating road surface conditions |
EP1479580A1 (en) * | 2003-05-22 | 2004-11-24 | Fuji Jukogyo Kabushiki Kaisha | Device for estimating friction coefficient on road surface of vehicle |
US7058541B2 (en) | 2003-05-22 | 2006-06-06 | Fuji Jukogyo Kabushiki Kaisha | Device for a vehicle for estimating friction coefficient of a road surface |
US7483808B2 (en) * | 2007-06-29 | 2009-01-27 | Caterpillar Inc. | System and method for measuring machine rolling resistance |
CN113299060A (en) * | 2021-04-29 | 2021-08-24 | 山西省交通科技研发有限公司 | Vehicle information coefficient acquisition system based on measured data |
CN113299060B (en) * | 2021-04-29 | 2022-11-08 | 山西省交通科技研发有限公司 | Vehicle information coefficient acquisition system based on measured data |
Also Published As
Publication number | Publication date |
---|---|
JP3424535B2 (en) | 2003-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6324461B1 (en) | Road surface condition estimating apparatus and variation reduction processing apparatus | |
EP0825080B1 (en) | Friction torque-gradient estimating apparatus, anti-lock brake controller and control start determination method | |
JP4331881B2 (en) | Road surface condition estimation device | |
US8825286B2 (en) | Method and device for determining a center of gravity of a motor vehicle | |
US9988057B2 (en) | Method and device for determining the mass of a motor vehicle, and a motor vehicle with a device of this type | |
KR100227559B1 (en) | Anti-lock brake controlling apparatus | |
JP5339121B2 (en) | Slip rate estimation device and method, and slip rate control device and method | |
Boisvert et al. | Estimators of wheel slip for electric vehicles using torque and encoder measurements | |
JP3435625B2 (en) | Road surface condition calculation device | |
JP3331939B2 (en) | Vehicle state quantity estimation device | |
JP3319407B2 (en) | Road surface condition estimation device | |
JPH11139293A (en) | Device for estimating condition of road surface | |
JP3454090B2 (en) | Anti-lock brake control device, torque gradient estimating device, and braking torque gradient estimating device | |
JP3436124B2 (en) | Fluctuation reduction processing device and road surface condition estimation device | |
JP3792756B2 (en) | Anti-lock brake control device | |
JPH1191539A (en) | Friction state arithmetic unit and braking force control device | |
JP4027585B2 (en) | Physical quantity estimation device and road surface friction state estimation device | |
JPH1111292A (en) | Wheel behavior amount servo control device and limit judgment device | |
JP3702649B2 (en) | Braking force estimation device and brake pressure estimation device | |
JP3769410B2 (en) | Physical quantity estimation device and ABS control device | |
JPH10264803A (en) | Braking pressure estimation device, anti-lock brake control device, and braking pressure control device | |
JP3640126B2 (en) | Slip rate servo controller | |
JP3644218B2 (en) | Anti-lock brake control device | |
JP2003237560A (en) | Maximum road surface friction coefficient-estimating device | |
JP3405189B2 (en) | Anti-lock brake control device, torque gradient estimating device, and braking torque gradient estimating device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |