JP3317196B2 - 荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置 - Google Patents
荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置Info
- Publication number
- JP3317196B2 JP3317196B2 JP18160897A JP18160897A JP3317196B2 JP 3317196 B2 JP3317196 B2 JP 3317196B2 JP 18160897 A JP18160897 A JP 18160897A JP 18160897 A JP18160897 A JP 18160897A JP 3317196 B2 JP3317196 B2 JP 3317196B2
- Authority
- JP
- Japan
- Prior art keywords
- wheel
- acceleration
- deceleration
- road surface
- vehicle body
- 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 - Fee Related
Links
Landscapes
- Regulating Braking Force (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、荷重配分推定装
置、車体加減速度演算装置、及び路面状態推定装置に係
り、より詳しくは、車輪にかかる荷重と共振レベルとが
正の相関関係にあることを利用して各車輪の荷重配分を
推定する荷重配分推定装置、荷重移動による共振レベル
の相違に基づいて車体の加減速度を推定する車体加減速
度推定装置、及び車体加減速度と車輪加減速度との間の
関係により路面状態を推定する路面状態推定装置に関す
る。
置、車体加減速度演算装置、及び路面状態推定装置に係
り、より詳しくは、車輪にかかる荷重と共振レベルとが
正の相関関係にあることを利用して各車輪の荷重配分を
推定する荷重配分推定装置、荷重移動による共振レベル
の相違に基づいて車体の加減速度を推定する車体加減速
度推定装置、及び車体加減速度と車輪加減速度との間の
関係により路面状態を推定する路面状態推定装置に関す
る。
【０００２】
【従来の技術】従来より、ハンドリング、ブレーキ制
動、加速などの運転操作に伴うロール、ノーズダイブ、
スクォート量を少なくし、最適なアライメント状態を保
持することで車両の操縦安定性を向上させる姿勢制御技
術が提案されている。
動、加速などの運転操作に伴うロール、ノーズダイブ、
スクォート量を少なくし、最適なアライメント状態を保
持することで車両の操縦安定性を向上させる姿勢制御技
術が提案されている。
【０００３】この姿勢制御技術のうちアンチロール制御
では、旋回時、車両の左右横方向に働く加減速度をＧセ
ンサ（加速度センサ）で検知して制御量を決定する。旋
回時は遠心力によって重心の移動が発生し、外輪側は目
標車高値より車高が下がり、内輪側は上がるので、アン
チロール制御では、外輪側のサスペンションのシリンダ
圧を上昇させると共に内輪側のシリンダ圧を下げ、遠心
力に起因して車体をロールさせようとする力を打ち消
し、横方向の車両姿勢をほぼフラットにする。なお、旋
回時初期の制御応答遅れを補償するために、操舵角と車
速との関係からロール開始を予測してロール予見制御を
行う。
では、旋回時、車両の左右横方向に働く加減速度をＧセ
ンサ（加速度センサ）で検知して制御量を決定する。旋
回時は遠心力によって重心の移動が発生し、外輪側は目
標車高値より車高が下がり、内輪側は上がるので、アン
チロール制御では、外輪側のサスペンションのシリンダ
圧を上昇させると共に内輪側のシリンダ圧を下げ、遠心
力に起因して車体をロールさせようとする力を打ち消
し、横方向の車両姿勢をほぼフラットにする。なお、旋
回時初期の制御応答遅れを補償するために、操舵角と車
速との関係からロール開始を予測してロール予見制御を
行う。
【０００４】また、アンチダイブ、アンチスクォート制
御では、減速時或いは加速時に、車両に作用する前後方
向の加減速度をＧセンサで検知し、検知した値に応じて
前後のサスペンションのシリンダ圧を上昇又は下降させ
ることにより前後方向の車両姿勢をほぼフラットにす
る。
御では、減速時或いは加速時に、車両に作用する前後方
向の加減速度をＧセンサで検知し、検知した値に応じて
前後のサスペンションのシリンダ圧を上昇又は下降させ
ることにより前後方向の車両姿勢をほぼフラットにす
る。
【０００５】このように従来の姿勢制御技術では、Ｇセ
ンサにより検出した横方向の加減速度又は前後方向の加
減速度に基づいて、サスペンションの制御量を決定して
いる。
ンサにより検出した横方向の加減速度又は前後方向の加
減速度に基づいて、サスペンションの制御量を決定して
いる。
【０００６】また、車両旋回時では、タイヤのコーナリ
ング特性の限界を越えてスピンが起こらないように、各
車輪にかかる荷重を制御する必要がある。
ング特性の限界を越えてスピンが起こらないように、各
車輪にかかる荷重を制御する必要がある。
【０００７】
【発明が解決しようとする課題】しかしながら、上記従
来技術では、以下のような問題がある。
来技術では、以下のような問題がある。
【０００８】すなわち、車両姿勢の制御やアンチスピン
制御などのような車両安定化制御を行う場合、各車輪の
荷重配分の推定が重要であるが、この荷重配分は、車体
特性は勿論のこと、乗員や荷物の荷重分布、サスペンシ
ョンの状態量などにより異なってくる。このため、Ｇセ
ンサなどにより検出された加速度などを用いた間接的な
推定では、推定が複雑になり過ぎると共に、推定値が不
正確となるおそれがあるので、簡単に良好な安定化制御
ができないという問題がある。また、路面状態によって
も制御を変える必要がある。
制御などのような車両安定化制御を行う場合、各車輪の
荷重配分の推定が重要であるが、この荷重配分は、車体
特性は勿論のこと、乗員や荷物の荷重分布、サスペンシ
ョンの状態量などにより異なってくる。このため、Ｇセ
ンサなどにより検出された加速度などを用いた間接的な
推定では、推定が複雑になり過ぎると共に、推定値が不
正確となるおそれがあるので、簡単に良好な安定化制御
ができないという問題がある。また、路面状態によって
も制御を変える必要がある。
【０００９】本発明は、上記事実に鑑みて成されたもの
で、各車輪の荷重配分を簡単かつ正確に推定できる荷重
配分推定装置、車体加減速度を正確に演算できる車体加
減速度演算装置、走行中の路面の状態を的確に推定でき
る路面状態推定装置を提供することを目的とする。
で、各車輪の荷重配分を簡単かつ正確に推定できる荷重
配分推定装置、車体加減速度を正確に演算できる車体加
減速度演算装置、走行中の路面の状態を的確に推定でき
る路面状態推定装置を提供することを目的とする。
【００１０】
（請求項１の発明）上記目的を達成するために請求項１
の発明は、車輪速度が微小振動しているときの各車輪の
荷重配分を推定する荷重配分推定装置において、車体と
車輪と路面とから構成される振動系の共振周波数での車
輪速度の微小振幅に基づく共振レベルを、各車輪毎に検
出する検出手段と、車輪にかかる荷重と該車輪の共振レ
ベルとが正の相関関係にあることに基づいて、前記検出
手段により検出された各車輪の共振レベルにより、各車
輪の荷重配分を推定する荷重配分推定手段と、を有する
ことを特徴とする。 （請求項１の発明の原理）重量Ｗの車体を備えた車両が
車体速度ωv で走行している時の車輪での振動現象、す
なわち車体と車輪と路面とによって構成される振動系の
振動現象を、車輪回転軸で等価的にモデル化した図２に
示すモデルを参照して考察する。
の発明は、車輪速度が微小振動しているときの各車輪の
荷重配分を推定する荷重配分推定装置において、車体と
車輪と路面とから構成される振動系の共振周波数での車
輪速度の微小振幅に基づく共振レベルを、各車輪毎に検
出する検出手段と、車輪にかかる荷重と該車輪の共振レ
ベルとが正の相関関係にあることに基づいて、前記検出
手段により検出された各車輪の共振レベルにより、各車
輪の荷重配分を推定する荷重配分推定手段と、を有する
ことを特徴とする。 （請求項１の発明の原理）重量Ｗの車体を備えた車両が
車体速度ωv で走行している時の車輪での振動現象、す
なわち車体と車輪と路面とによって構成される振動系の
振動現象を、車輪回転軸で等価的にモデル化した図２に
示すモデルを参照して考察する。
【００１１】図２のモデルにおいて、ブレーキ力は、路
面と接するタイヤのトレッド１１５の表面を介して路面
に作用するが、このブレーキ力は実際には路面からの反
作用（制動力）として車体に作用するため、車体重量の
回転軸換算の等価モデル１１７はタイヤのトレッドと路
面との間の摩擦要素１１６（路面μ）を介して車輪１１
３と反対側に連結したものとなる。これは、シャシーダ
イナモ装置のように、車輪下の大きな慣性、すなわち車
輪と反対側の質量で車体の重量を模擬することができる
ことと同様である。
面と接するタイヤのトレッド１１５の表面を介して路面
に作用するが、このブレーキ力は実際には路面からの反
作用（制動力）として車体に作用するため、車体重量の
回転軸換算の等価モデル１１７はタイヤのトレッドと路
面との間の摩擦要素１１６（路面μ）を介して車輪１１
３と反対側に連結したものとなる。これは、シャシーダ
イナモ装置のように、車輪下の大きな慣性、すなわち車
輪と反対側の質量で車体の重量を模擬することができる
ことと同様である。
【００１２】図２でタイヤリムを含んだ車輪１１３の慣
性をＪw 、リムとトレッド１１５との間のばね要素１１
４のばね定数をＫ、車輪半径をＲ、トレッド１１５の慣
性をＪt 、トレッド１１５と路面との間の摩擦要素１１
６の摩擦係数をμ、車体の重量の回転軸換算の等価モデ
ル１１７の慣性をＪV とすると、ホイールシリンダ圧に
より生じるブレーキトルクＴb ’から車輪速ωw までの
伝達特性は、車輪運動の方程式より、
性をＪw 、リムとトレッド１１５との間のばね要素１１
４のばね定数をＫ、車輪半径をＲ、トレッド１１５の慣
性をＪt 、トレッド１１５と路面との間の摩擦要素１１
６の摩擦係数をμ、車体の重量の回転軸換算の等価モデ
ル１１７の慣性をＪV とすると、ホイールシリンダ圧に
より生じるブレーキトルクＴb ’から車輪速ωw までの
伝達特性は、車輪運動の方程式より、
【００１３】
【数１】
【００１４】となる。なお、ｓはラプラス変換の演算子
である。タイヤが路面にグリップしている時は、トレッ
ド１１５と車体等価モデル１１７とが直結されていると
考えると、車体等価モデル１１７とトレッド１１５との
和の慣性と、車輪１１３の慣性とが共振する。すなわ
ち、この振動系は、車輪と車体と路面とから構成された
車輪共振系とみなすことができる。このときの車輪共振
系の共振周波数ω∞は、(1) 式の伝達特性において、 ω∞＝√｛（Ｊw ＋Ｊt ＋Ｊv ）Ｋ／Ｊw （Ｊt ＋Ｊv ）｝／２π (2) となる。この状態は、路面μ特性（スリップ速度に対す
るμの変化特性）において、ピークμに達する前のスリ
ップ速度の領域に対応する。
である。タイヤが路面にグリップしている時は、トレッ
ド１１５と車体等価モデル１１７とが直結されていると
考えると、車体等価モデル１１７とトレッド１１５との
和の慣性と、車輪１１３の慣性とが共振する。すなわ
ち、この振動系は、車輪と車体と路面とから構成された
車輪共振系とみなすことができる。このときの車輪共振
系の共振周波数ω∞は、(1) 式の伝達特性において、 ω∞＝√｛（Ｊw ＋Ｊt ＋Ｊv ）Ｋ／Ｊw （Ｊt ＋Ｊv ）｝／２π (2) となる。この状態は、路面μ特性（スリップ速度に対す
るμの変化特性）において、ピークμに達する前のスリ
ップ速度の領域に対応する。
【００１５】逆に、タイヤの摩擦係数μがピークμに近
づく場合には、タイヤ表面の摩擦係数μがスリップ率に
対して変化し難くなり、トレッド１１５の慣性の振動に
伴う成分は車体等価モデル１１７に影響しなくなる。つ
まり等価的にトレッド１１５と車体等価モデル１１７と
が分離され、トレッド１１５と車輪１１３とが共振を起
こすことになる。このときの車輪共振系は、車輪と路面
とから構成されているとみなすことができ、その共振周
波数ω∞’は、(2) 式において、車体等価慣性Ｊv を０
とおいたものと等しくなる。すなわち、 ω∞' ＝√｛（Ｊw ＋Ｊt ）Ｋ／Ｊw Ｊt ）｝／２π (3) となる。この状態は、路面μ特性においてピークμ近傍
のスリップ速度領域に対応する。
づく場合には、タイヤ表面の摩擦係数μがスリップ率に
対して変化し難くなり、トレッド１１５の慣性の振動に
伴う成分は車体等価モデル１１７に影響しなくなる。つ
まり等価的にトレッド１１５と車体等価モデル１１７と
が分離され、トレッド１１５と車輪１１３とが共振を起
こすことになる。このときの車輪共振系は、車輪と路面
とから構成されているとみなすことができ、その共振周
波数ω∞’は、(2) 式において、車体等価慣性Ｊv を０
とおいたものと等しくなる。すなわち、 ω∞' ＝√｛（Ｊw ＋Ｊt ）Ｋ／Ｊw Ｊt ）｝／２π (3) となる。この状態は、路面μ特性においてピークμ近傍
のスリップ速度領域に対応する。
【００１６】(2) と(3) 式とを比較し、車体等価慣性Ｊ
v が車輪慣性Ｊw 、トレッド慣性Ｊ t より大きいと仮定
すると、(3) 式の場合の車輪共振系の共振周波数ω∞’
は(2) 式のω∞よりも高周波数側にシフトすることにな
る。
v が車輪慣性Ｊw 、トレッド慣性Ｊ t より大きいと仮定
すると、(3) 式の場合の車輪共振系の共振周波数ω∞’
は(2) 式のω∞よりも高周波数側にシフトすることにな
る。
【００１７】ところで、上記車輪共振系を構成する車輪
が、内力又は外力の作用によって微小振動するとき、そ
の車輪速度は、平均的な車輪速度の回りに該共振系の共
振周波数をピークとする振動で微小振動する。そこで、
タイヤが路面にグリップしているときの共振周波数ω∞
での車輪速度の振動強度（振幅）に基づいて、該共振系
の共振レベルを表すことができる。なお、この共振レベ
ルとして、車輪速度の微小振幅値そのものを用いても良
い。
が、内力又は外力の作用によって微小振動するとき、そ
の車輪速度は、平均的な車輪速度の回りに該共振系の共
振周波数をピークとする振動で微小振動する。そこで、
タイヤが路面にグリップしているときの共振周波数ω∞
での車輪速度の振動強度（振幅）に基づいて、該共振系
の共振レベルを表すことができる。なお、この共振レベ
ルとして、車輪速度の微小振幅値そのものを用いても良
い。
【００１８】ここで、図１４に示すように、タイヤが路
面にグリップしている状態からピークμに近づくと、上
述したように共振周波数は、ω∞からより高周波数側の
ω∞’に上昇する。このとき、タイヤが路面にグリップ
した状態で最大値となっていたω∞での共振レベルは、
共振周波数の上昇と共に減少する。従って、この共振レ
ベルに基づいて、タイヤと路面との間の摩擦状態を判定
することができる。
面にグリップしている状態からピークμに近づくと、上
述したように共振周波数は、ω∞からより高周波数側の
ω∞’に上昇する。このとき、タイヤが路面にグリップ
した状態で最大値となっていたω∞での共振レベルは、
共振周波数の上昇と共に減少する。従って、この共振レ
ベルに基づいて、タイヤと路面との間の摩擦状態を判定
することができる。
【００１９】また、この共振レベルは、タイヤが路面に
強くグリップしているほど大きな値をとることから、車
輪にかかる荷重（輪荷重）と正の相関があることが予想
できる。実際、評価試験等により、車両旋回時におい
て、荷重が増加する旋回外輪のタイヤの共振レベルが増
加し、荷重が減少する旋回内輪の共振レベルが減少する
結果が得られている。
強くグリップしているほど大きな値をとることから、車
輪にかかる荷重（輪荷重）と正の相関があることが予想
できる。実際、評価試験等により、車両旋回時におい
て、荷重が増加する旋回外輪のタイヤの共振レベルが増
加し、荷重が減少する旋回内輪の共振レベルが減少する
結果が得られている。
【００２０】このように車輪にかかる荷重と該車輪の共
振レベルとは正の相関関係にあるので、旋回時のみなら
ず、車両が直線上で加減速度運動（加速度運動又は減速
度運動）を行うことにより前後方向に輪荷重が移動する
際にも、タイヤの共振レベルは、輪荷重が移動した方向
の車輪では増大し、その反対方向の車輪では減少する。
例えば、ブレーキ制動時には、前輪側に荷重が移動する
ため、前輪の共振レベルが増大し、後輪の共振レベルが
減少する。一方、アクセルペダルを踏んだ加速時には、
後輪側に荷重が移動するため、後輪の共振レベルが増大
し、前輪の共振レベルが減少する。
振レベルとは正の相関関係にあるので、旋回時のみなら
ず、車両が直線上で加減速度運動（加速度運動又は減速
度運動）を行うことにより前後方向に輪荷重が移動する
際にも、タイヤの共振レベルは、輪荷重が移動した方向
の車輪では増大し、その反対方向の車輪では減少する。
例えば、ブレーキ制動時には、前輪側に荷重が移動する
ため、前輪の共振レベルが増大し、後輪の共振レベルが
減少する。一方、アクセルペダルを踏んだ加速時には、
後輪側に荷重が移動するため、後輪の共振レベルが増大
し、前輪の共振レベルが減少する。
【００２１】そこで、請求項１の発明では、検出手段
が、各車輪の共振レベルＧi （i=1,2,3,4)を検出する。
そして、荷重配分推定手段が、検出された各車輪の共振
レベルにより、各車輪の荷重配分を推定する。
が、各車輪の共振レベルＧi （i=1,2,3,4)を検出する。
そして、荷重配分推定手段が、検出された各車輪の共振
レベルにより、各車輪の荷重配分を推定する。
【００２２】例えば、共振レベルに対し輪荷重が正の相
関関数Ｆで変化する場合、第ｉ輪の荷重Ｗi を、 Ｗi ＝ Ｆ（Ｇi ） (4) により求める。
関関数Ｆで変化する場合、第ｉ輪の荷重Ｗi を、 Ｗi ＝ Ｆ（Ｇi ） (4) により求める。
【００２３】また、もっと簡単に次式のような略比例関
係により求めても良い。 Ｗi ≒ ｋ・Ｇi (5) ただし、ｋは比例定数である。
係により求めても良い。 Ｗi ≒ ｋ・Ｇi (5) ただし、ｋは比例定数である。
【００２４】さらに、各車輪の荷重がわかれば、前後方
向及び横方向の移動荷重や荷重移動の配分比率などもた
だちに求めることができる。
向及び横方向の移動荷重や荷重移動の配分比率などもた
だちに求めることができる。
【００２５】このように本発明では、各車輪の荷重配分
が実際に路面に接地している車輪の共振レベルにより直
接的に求められるので、Ｇセンサなどにより検出された
加速度に基づいて推定された荷重配分よりも正確な値が
得られる。また、これにより、車両安定化制御をより正
確かつ安定に行うことができる。例えば、コーナリング
時のように横方向の加速度で走行する場合の前後方向の
荷重移動は、車体特性、乗員や荷物の荷重分布、サスペ
ンションの状態量などにより異なってくるが、このよう
な場合でも、本発明では、これら諸量により推定するこ
となく各車輪の荷重配分を直接的に求めることができ
る。 （請求項２の発明）また、請求項２の発明は、車輪の車
輪速度が微小振動しているときの車体の加減速度を演算
する車体加減速度演算装置において、車体と車輪と路面
とから構成される振動系の共振周波数での車輪速度の微
小振幅に基づく共振レベルを、前輪と後輪について各々
検出する第１の検出手段、及び右輪と左輪について各々
検出する第２の検出手段の少なくともいずれかの検出手
段と、前記第１の検出手段により検出された前輪の共振
レベルと後輪の共振レベルとに基づいて車体の前後方向
の加減速度を演算する第１の演算手段、及び前記第２の
検出手段により検出された右輪の共振レベルと左輪の共
振レベルとに基づいて車体の横方向の加減速度を演算す
る第２の演算手段の少なくともいずれかの演算手段と、
を有することを特徴とする。 （請求項２の発明の原理）図１に、輪荷重の前後方向の
配分を求めるため、車両が加減速度運動を行っていると
きのハーフカーモデルを示す。なお、図１の各物理パラ
メータは、以下の通りとなる。
が実際に路面に接地している車輪の共振レベルにより直
接的に求められるので、Ｇセンサなどにより検出された
加速度に基づいて推定された荷重配分よりも正確な値が
得られる。また、これにより、車両安定化制御をより正
確かつ安定に行うことができる。例えば、コーナリング
時のように横方向の加速度で走行する場合の前後方向の
荷重移動は、車体特性、乗員や荷物の荷重分布、サスペ
ンションの状態量などにより異なってくるが、このよう
な場合でも、本発明では、これら諸量により推定するこ
となく各車輪の荷重配分を直接的に求めることができ
る。 （請求項２の発明）また、請求項２の発明は、車輪の車
輪速度が微小振動しているときの車体の加減速度を演算
する車体加減速度演算装置において、車体と車輪と路面
とから構成される振動系の共振周波数での車輪速度の微
小振幅に基づく共振レベルを、前輪と後輪について各々
検出する第１の検出手段、及び右輪と左輪について各々
検出する第２の検出手段の少なくともいずれかの検出手
段と、前記第１の検出手段により検出された前輪の共振
レベルと後輪の共振レベルとに基づいて車体の前後方向
の加減速度を演算する第１の演算手段、及び前記第２の
検出手段により検出された右輪の共振レベルと左輪の共
振レベルとに基づいて車体の横方向の加減速度を演算す
る第２の演算手段の少なくともいずれかの演算手段と、
を有することを特徴とする。 （請求項２の発明の原理）図１に、輪荷重の前後方向の
配分を求めるため、車両が加減速度運動を行っていると
きのハーフカーモデルを示す。なお、図１の各物理パラ
メータは、以下の通りとなる。
【００２６】ｖ ： 車体速度 Ｇ ： 車両の重心 Ｍ ： 車両質量 ｇ ： 重力加速度 ｈ ： 重心Ｇの路面からの高さ Ｒr,f ： 車輪の有効半径（ｆは前輪、ｒは後輪に関
する） Ｗr,f ： 輪荷重（ｆは前輪、ｒは後輪に関する） ωr,f ： 車輪速度（ｆは前輪、ｒは後輪に関する） μr,f ： 車輪と路面との間の摩擦係数（ｆは前輪、
ｒは後輪に関する） Ｐr,f ： ブレーキ力（ｆは前輪、ｒは後輪に関す
る） Ｋr,f ： ブレーキ力のトルク変換係数（ｆは前輪、
ｒは後輪に関する） Ｌ ： 前輪中心から後輪中心までの間隔 Ｌ1 ： 前輪中心から重心Ｇの水平位置までの間隔 Ｌ2 ： 後輪中心から重心Ｇの水平位置までの間隔 図１のモデルの車両が、ｄｖ／ｄｔ（ｔは時間）の加速
度（減速度）で走行している場合、該車両には、加速度
（減速度）方向とは反対方向に慣性力が作用する。この
慣性力の作用は、重心Ｇに対し、 Ｆ＝−Ｍｄｖ／ｄｔ (4−２) の力の作用と等価である（図１の力Ｆの方向はブレーキ
制動している減速度時のもの）。
する） Ｗr,f ： 輪荷重（ｆは前輪、ｒは後輪に関する） ωr,f ： 車輪速度（ｆは前輪、ｒは後輪に関する） μr,f ： 車輪と路面との間の摩擦係数（ｆは前輪、
ｒは後輪に関する） Ｐr,f ： ブレーキ力（ｆは前輪、ｒは後輪に関す
る） Ｋr,f ： ブレーキ力のトルク変換係数（ｆは前輪、
ｒは後輪に関する） Ｌ ： 前輪中心から後輪中心までの間隔 Ｌ1 ： 前輪中心から重心Ｇの水平位置までの間隔 Ｌ2 ： 後輪中心から重心Ｇの水平位置までの間隔 図１のモデルの車両が、ｄｖ／ｄｔ（ｔは時間）の加速
度（減速度）で走行している場合、該車両には、加速度
（減速度）方向とは反対方向に慣性力が作用する。この
慣性力の作用は、重心Ｇに対し、 Ｆ＝−Ｍｄｖ／ｄｔ (4−２) の力の作用と等価である（図１の力Ｆの方向はブレーキ
制動している減速度時のもの）。
【００２７】ここで、この慣性力Ｆによる前輪の荷重の
増分をΔＷf 、後輪の荷重の増分をΔＷr とすると、慣
性力Ｆは水平方向の成分しかなく、よって車両全体で
は、鉛直方向の力の増分が無いことから、 ΔＷf ＋ΔＷr ＝０ (5−２) が成立する。
増分をΔＷf 、後輪の荷重の増分をΔＷr とすると、慣
性力Ｆは水平方向の成分しかなく、よって車両全体で
は、鉛直方向の力の増分が無いことから、 ΔＷf ＋ΔＷr ＝０ (5−２) が成立する。
【００２８】また、路面からの高さｈの重心Ｇに作用す
る慣性力Ｆによって、重心Ｇから路面まで鉛直方向に垂
ろした点Ｒの回りにはＮ＝Ｆ×ｈのトルクが作用する。
このトルクＮは、前輪荷重の増分ΔＷf 及び後輪荷重の
増分ΔＷr による点Ｒの回りのトルクと等しくなるはず
である。よって、 Ｆｈ ＝ ΔＷf Ｌ1 −ΔＷr Ｌ2 (6−２) が成立する。
る慣性力Ｆによって、重心Ｇから路面まで鉛直方向に垂
ろした点Ｒの回りにはＮ＝Ｆ×ｈのトルクが作用する。
このトルクＮは、前輪荷重の増分ΔＷf 及び後輪荷重の
増分ΔＷr による点Ｒの回りのトルクと等しくなるはず
である。よって、 Ｆｈ ＝ ΔＷf Ｌ1 −ΔＷr Ｌ2 (6−２) が成立する。
【００２９】(4−２) 〜(6−２) 式より、前輪及び後輪
の荷重増分は、
の荷重増分は、
【００３０】
【数２】
【００３１】となる。ただし、・は時間に関する１階微
分を示す。(7) 式より、前輪の荷重Ｗf 及び後輪の荷重
Ｗr を重力加速度ｇを用いて表すと、次式のようにな
る。
分を示す。(7) 式より、前輪の荷重Ｗf 及び後輪の荷重
Ｗr を重力加速度ｇを用いて表すと、次式のようにな
る。
【００３２】
【数３】
【００３３】となる。ただし、Ｗf0，Ｗr0は、静止時の
輪荷重（ｆは前輪、ｒは後輪に関する）であり、Ｗは車
体の荷重となる。
輪荷重（ｆは前輪、ｒは後輪に関する）であり、Ｗは車
体の荷重となる。
【００３４】ここで、加速度ｄｖ／ｄｔを、重力加速度
ｇに対する比α（以下、「加減速度」という）で表す
と、
ｇに対する比α（以下、「加減速度」という）で表す
と、
【００３５】
【数４】
【００３６】となる。このαを、輪荷重で表すと、(8)
〜(11)式より、
〜(11)式より、
【００３７】
【数５】
【００３８】が得られる。さらに、(12)式において、Ｗ
f0 ≒Ｗr0を用いると、
f0 ≒Ｗr0を用いると、
【００３９】
【数６】
【００４０】と近似できる。今、仮に、共振レベルが輪
荷重と比例すると仮定すると、(5) 式の比例関係を用い
て、(13)式より、
荷重と比例すると仮定すると、(5) 式の比例関係を用い
て、(13)式より、
【００４１】
【数７】
【００４２】が得られる。ただし、Ｇf は前輪の共振レ
ベル、Ｇr は後輪の共振レベルである。
ベル、Ｇr は後輪の共振レベルである。
【００４３】そこで、請求項２の発明では、第１の検出
手段が前輪の共振レベルＧf 及び後輪の共振レベルＧr
を検出し、第１の演算手段が、 (14) に基づいて前後方
向の車体加減速度αを演算する。
手段が前輪の共振レベルＧf 及び後輪の共振レベルＧr
を検出し、第１の演算手段が、 (14) に基づいて前後方
向の車体加減速度αを演算する。
【００４４】また、上記の加減速度の導出は、左右の車
輪についても同様に適用することができ、この場合、横
方向の車体加減速度α’は、次式により表される。
輪についても同様に適用することができ、この場合、横
方向の車体加減速度α’は、次式により表される。
【００４５】
【数８】
【００４６】ただし、(15)式のＧR は右輪の共振レベ
ル、ＧL は左輪の共振レベルである。また、Ｌ’は、右
輪と左輪の車輪間隔である。
ル、ＧL は左輪の共振レベルである。また、Ｌ’は、右
輪と左輪の車輪間隔である。
【００４７】そこで、請求項２の発明では、第２の検出
手段が右輪の共振レベルＧR 及び左輪の共振レベルＧL
を検出し、第２の演算手段が、 (15) に基づいて横方向
の車体加減速度α’を演算する。
手段が右輪の共振レベルＧR 及び左輪の共振レベルＧL
を検出し、第２の演算手段が、 (15) に基づいて横方向
の車体加減速度α’を演算する。
【００４８】なお、(14)、(15) 式は、１つの演算例で
あって、本発明はこれに限定されるものではなく、他の
演算方法、例えば(4) 式より演算された輪荷重に基づい
て、車体加減速度を演算する方法も適用可能である。
あって、本発明はこれに限定されるものではなく、他の
演算方法、例えば(4) 式より演算された輪荷重に基づい
て、車体加減速度を演算する方法も適用可能である。
【００４９】このように本発明では、Ｇセンサなどを用
いることなく、車体の加減速度を精度良く推定し、種々
の車両安定化制御に有効に用いることができる。 （請求項３の発明）また、請求項３の発明は、車輪の車
輪速度が微小振動している車体が走行する路面の状態を
推定する路面状態推定装置において、車体と車輪と路面
とから構成される振動系の共振周波数での車輪速度の微
小振幅に基づく共振レベルを、前輪と後輪について各々
検出する検出手段と、前記検出手段により検出された前
輪の共振レベルと後輪の共振レベルとに基づいて車体の
加減速度を演算する演算手段と、車輪の加減速度を検出
する車輪加減速度検出手段と、前記演算手段により演算
された車体の加減速度及び前記車輪加減速度検出手段に
より検出された車輪の加減速度の間の関係と、予めデー
タとして備えられている車体の加減速度及び車輪の加減
速度の間の路面状態毎の関係と、の比較に基づいて路面
状態を推定する路面推定手段と、を有することを特徴と
する。 （請求項３の発明の原理）路面から反力として各車輪に
作用することにより車体を動かす制動力は、駆動トルク
とブレーキトルクとによる作用が同じ場合でも、高い摩
擦係数μの路面ほど大きくなる。このため、車輪の加減
速度αw に対する車体の加減速度αの変化特性は、図３
に示すように、路面状態（高μ路、中μ路、低μ路）に
よって異なってくる。すなわち、図３に示すように、低
μ路、中μ路、高μ路の各場合において、車輪の加減速
度αw の同じ値に対し、高いμの路面ほど、車体の加減
速度の絶対値が大きい。また、車輪の加減速度αw の増
大に対し、高いμの路面ほど車体の加減速度の変化率が
大きい。
いることなく、車体の加減速度を精度良く推定し、種々
の車両安定化制御に有効に用いることができる。 （請求項３の発明）また、請求項３の発明は、車輪の車
輪速度が微小振動している車体が走行する路面の状態を
推定する路面状態推定装置において、車体と車輪と路面
とから構成される振動系の共振周波数での車輪速度の微
小振幅に基づく共振レベルを、前輪と後輪について各々
検出する検出手段と、前記検出手段により検出された前
輪の共振レベルと後輪の共振レベルとに基づいて車体の
加減速度を演算する演算手段と、車輪の加減速度を検出
する車輪加減速度検出手段と、前記演算手段により演算
された車体の加減速度及び前記車輪加減速度検出手段に
より検出された車輪の加減速度の間の関係と、予めデー
タとして備えられている車体の加減速度及び車輪の加減
速度の間の路面状態毎の関係と、の比較に基づいて路面
状態を推定する路面推定手段と、を有することを特徴と
する。 （請求項３の発明の原理）路面から反力として各車輪に
作用することにより車体を動かす制動力は、駆動トルク
とブレーキトルクとによる作用が同じ場合でも、高い摩
擦係数μの路面ほど大きくなる。このため、車輪の加減
速度αw に対する車体の加減速度αの変化特性は、図３
に示すように、路面状態（高μ路、中μ路、低μ路）に
よって異なってくる。すなわち、図３に示すように、低
μ路、中μ路、高μ路の各場合において、車輪の加減速
度αw の同じ値に対し、高いμの路面ほど、車体の加減
速度の絶対値が大きい。また、車輪の加減速度αw の増
大に対し、高いμの路面ほど車体の加減速度の変化率が
大きい。
【００５０】図３のような特性を逆に利用し、本発明で
は、演算手段により車体の加減速度αを演算すると共
に、車輪加減速度演算手段により車輪の加減速度αw を
随時演算し、これらの加減速度のデータの間の関係と、
予めデータとして備えられている車体の加減速度及び車
輪の加減速度の間の路面状態毎の関係と、の比較に基づ
いて路面状態を推定する。例えば、観測されたαとαw
との関係が、図３のどの路面状態の関係に最も近いかを
判定することにより、現在走行中の路面の状態を推定す
ることができる。
は、演算手段により車体の加減速度αを演算すると共
に、車輪加減速度演算手段により車輪の加減速度αw を
随時演算し、これらの加減速度のデータの間の関係と、
予めデータとして備えられている車体の加減速度及び車
輪の加減速度の間の路面状態毎の関係と、の比較に基づ
いて路面状態を推定する。例えば、観測されたαとαw
との関係が、図３のどの路面状態の関係に最も近いかを
判定することにより、現在走行中の路面の状態を推定す
ることができる。
【００５１】このように本発明では、路面状態を推定す
ることができるので、推定された路面状態に応じて制御
パラメータを変更する処理等を行えば、種々の車両安定
化制御をより高いレベルで実行することができる。
ることができるので、推定された路面状態に応じて制御
パラメータを変更する処理等を行えば、種々の車両安定
化制御をより高いレベルで実行することができる。
【００５２】なお、本発明に係る共振レベルを、次のよ
うな微小ゲインで表すようにしても良い。ただし、この
微小ゲインを演算する際には、車輪と車体と路面とから
なる振動系の共振周波数ω∞（(3) 式) でブレーキ圧Ｐ
b を微小励振することが前提となる。
うな微小ゲインで表すようにしても良い。ただし、この
微小ゲインを演算する際には、車輪と車体と路面とから
なる振動系の共振周波数ω∞（(3) 式) でブレーキ圧Ｐ
b を微小励振することが前提となる。
【００５３】まず、ブレーキ圧Ｐb が微小励振される
と、車輪速度ωw も平均的な車輪速度の回りに共振周波
数ω∞で微小振動する。ここで、このときのブレーキ圧
Ｐb の共振周波数ω∞の微小振幅をＰv 、車輪速度の共
振周波数ω∞の微小振幅をωwvとした場合、微小ゲイン
Ｇd を、 Ｇd ＝ωwv／Ｐv (16) により演算する。なお、この微小ゲインＧd を、ブレー
キ圧Ｐb に対する車輪速ωw の比（ωw ／Ｐb ）の共振
周波数ω∞の振動成分とみなし、 Ｇd ＝（（ωw ／Ｐb ）｜ｓ＝ｊω∞） (17) と表すこともできる。
と、車輪速度ωw も平均的な車輪速度の回りに共振周波
数ω∞で微小振動する。ここで、このときのブレーキ圧
Ｐb の共振周波数ω∞の微小振幅をＰv 、車輪速度の共
振周波数ω∞の微小振幅をωwvとした場合、微小ゲイン
Ｇd を、 Ｇd ＝ωwv／Ｐv (16) により演算する。なお、この微小ゲインＧd を、ブレー
キ圧Ｐb に対する車輪速ωw の比（ωw ／Ｐb ）の共振
周波数ω∞の振動成分とみなし、 Ｇd ＝（（ωw ／Ｐb ）｜ｓ＝ｊω∞） (17) と表すこともできる。
【００５４】外力による車輪速度の微小振動がきわめて
小さい場合などでも、ブレーキ圧の励振による微小ゲイ
ンは、安定した共振レベルを示すことができ、上記各発
明における推定精度を向上させることができる。また、
ブレーキ圧微小振幅で正規化しているので、ブレーキ圧
が変動する条件でも正確な値が得られる。
小さい場合などでも、ブレーキ圧の励振による微小ゲイ
ンは、安定した共振レベルを示すことができ、上記各発
明における推定精度を向上させることができる。また、
ブレーキ圧微小振幅で正規化しているので、ブレーキ圧
が変動する条件でも正確な値が得られる。
【００５５】
【発明の実施の形態】以下、本発明の実施の形態を、図
面に基づいて詳細に説明する。 （第１の実施の形態）図４には、本発明の第１の実施の
形態に係る荷重配分推定装置の構成ブロック図が示され
ている。同図に示すように、本荷重配分推定装置１０が
適用される車両の右前輪（ＦＲ）、左前輪（ＦＬ）、右
後輪（ＲＲ）、左後輪（ＲＬ）には、各車輪の車輪速度
を各々検出する車輪速センサ２０ａ、２０ｂ、２０ｃ、
２０ｄが備えられている。
面に基づいて詳細に説明する。 （第１の実施の形態）図４には、本発明の第１の実施の
形態に係る荷重配分推定装置の構成ブロック図が示され
ている。同図に示すように、本荷重配分推定装置１０が
適用される車両の右前輪（ＦＲ）、左前輪（ＦＬ）、右
後輪（ＲＲ）、左後輪（ＲＬ）には、各車輪の車輪速度
を各々検出する車輪速センサ２０ａ、２０ｂ、２０ｃ、
２０ｄが備えられている。
【００５６】そして、本実施の形態の荷重配分推定装置
１０には、車輪速センサにより検出された各車輪速度ω
1 、ω2 、ω3 、ω4 から、各車輪の微小ゲインＧd1、
Ｇd2、Ｇd3、Ｇd4を各々演算する微小ゲイン演算部２２
ａ、２２ｂ、２２ｃ、２２ｄが備えられており、これら
の微小ゲイン演算部は、演算された各車輪の微小ゲイン
に基づいて各車輪の荷重配分を推定する荷重配分推定部
２６に接続されている。また、前輪用の微小ゲイン演算
部２２ａ、２２ｂと荷重配分推定部２６との間には、微
小ゲインの演算値を車体速度に応じて遅延させて出力す
る遅延部２４ａ、２４ｂが介在されている。
１０には、車輪速センサにより検出された各車輪速度ω
1 、ω2 、ω3 、ω4 から、各車輪の微小ゲインＧd1、
Ｇd2、Ｇd3、Ｇd4を各々演算する微小ゲイン演算部２２
ａ、２２ｂ、２２ｃ、２２ｄが備えられており、これら
の微小ゲイン演算部は、演算された各車輪の微小ゲイン
に基づいて各車輪の荷重配分を推定する荷重配分推定部
２６に接続されている。また、前輪用の微小ゲイン演算
部２２ａ、２２ｂと荷重配分推定部２６との間には、微
小ゲインの演算値を車体速度に応じて遅延させて出力す
る遅延部２４ａ、２４ｂが介在されている。
【００５７】さらに、遅延部２４ａ、２４ｂの出力端及
び微小ゲイン演算部２２ｃ、２２ｄの出力端には、遅延
された前輪の微小ゲインＧd1、Ｇd2と後輪の微小ゲイン
Ｇd3、Ｇd4とに基づいて(14)式、(15)式から車体の前後
方向の加減速度及び横方向の加減速度を演算出力する車
体加減速度演算部３０が接続されている。
び微小ゲイン演算部２２ｃ、２２ｄの出力端には、遅延
された前輪の微小ゲインＧd1、Ｇd2と後輪の微小ゲイン
Ｇd3、Ｇd4とに基づいて(14)式、(15)式から車体の前後
方向の加減速度及び横方向の加減速度を演算出力する車
体加減速度演算部３０が接続されている。
【００５８】また、荷重配分推定装置１０には、車輪速
センサにより検出された各車輪速度ω1 、ω2 、ω3 、
ω4 から各車輪の加減速度を演算する車輪加減速度演算
部２８と、演算された各車輪の加減速度の時間的な変化
率を演算する加減速度変化率演算部３２と、演算された
車体の加減速度及び車輪の加減速度の間の関係が、どの
路面状態での関係に最も近いかを判断することにより現
在走行中の路面の状態を推定する路面状態推定部３４
と、が備えられている。なお、車体の加減速度及び車輪
の加減速度の間の関係は、路面状態だけでなく加減速度
変化率によってもそれぞれ異なるため、本実施の形態の
路面状態推定部３４では、演算された加減速度変化率に
応じて判断基準を変更している。
センサにより検出された各車輪速度ω1 、ω2 、ω3 、
ω4 から各車輪の加減速度を演算する車輪加減速度演算
部２８と、演算された各車輪の加減速度の時間的な変化
率を演算する加減速度変化率演算部３２と、演算された
車体の加減速度及び車輪の加減速度の間の関係が、どの
路面状態での関係に最も近いかを判断することにより現
在走行中の路面の状態を推定する路面状態推定部３４
と、が備えられている。なお、車体の加減速度及び車輪
の加減速度の間の関係は、路面状態だけでなく加減速度
変化率によってもそれぞれ異なるため、本実施の形態の
路面状態推定部３４では、演算された加減速度変化率に
応じて判断基準を変更している。
【００５９】このうち荷重配分推定部２６は、(5) 式で
示したように共振レベル（微小ゲイン）と輪荷重とが略
比例関係にあると仮定した場合、例えば、図５のように
構成することが可能となる。
示したように共振レベル（微小ゲイン）と輪荷重とが略
比例関係にあると仮定した場合、例えば、図５のように
構成することが可能となる。
【００６０】図５に示すように、本実施の形態の荷重配
分推定部２６は、各車輪の微小ゲインＧd1，Ｇd2，
Ｇd3，Ｇd4を総和することにより、車体重量Ｗに対応す
る共振レベルＧw を演算する総和器３６と、総和Ｇw に
１／４を乗じる乗算器３８とを備えている。
分推定部２６は、各車輪の微小ゲインＧd1，Ｇd2，
Ｇd3，Ｇd4を総和することにより、車体重量Ｗに対応す
る共振レベルＧw を演算する総和器３６と、総和Ｇw に
１／４を乗じる乗算器３８とを備えている。
【００６１】また、荷重配分推定部２６は、各車輪の微
小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4のいずれかと、乗算器
３８の出力Ｇw ／４との差の絶対値（以下、単に「差
分」という）｜Ｇd1−Ｇw ／４｜，｜Ｇd2−Ｇw ／４
｜，｜Ｇd3−Ｇw ／４｜，｜Ｇd4−Ｇw ／４｜を各々演
算する差分器６２ａ、６２ｂ、６２ｃ、６２ｄを備えて
いる。また、これらの差分器の出力端には、演算された
各差分値に基づいて車体の前後の荷重移動比率及び左右
の荷重移動比率を演算出力する荷重移動比率演算器６５
が接続されている。
小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4のいずれかと、乗算器
３８の出力Ｇw ／４との差の絶対値（以下、単に「差
分」という）｜Ｇd1−Ｇw ／４｜，｜Ｇd2−Ｇw ／４
｜，｜Ｇd3−Ｇw ／４｜，｜Ｇd4−Ｇw ／４｜を各々演
算する差分器６２ａ、６２ｂ、６２ｃ、６２ｄを備えて
いる。また、これらの差分器の出力端には、演算された
各差分値に基づいて車体の前後の荷重移動比率及び左右
の荷重移動比率を演算出力する荷重移動比率演算器６５
が接続されている。
【００６２】さらに、荷重配分推定部２６は、各車輪の
微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に対し比例定数ｋを
各々乗算する乗算器６４ａ，６４ｂ，６４ｃ，６４ｄを
備えている。
微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に対し比例定数ｋを
各々乗算する乗算器６４ａ，６４ｂ，６４ｃ，６４ｄを
備えている。
【００６３】次に、車体加減速度演算部３０は、(5) 式
で示したように共振レベル（微小ゲイン）と輪荷重とが
略比例関係にあると仮定した場合、例えば、図６のよう
に構成することが可能となる。
で示したように共振レベル（微小ゲイン）と輪荷重とが
略比例関係にあると仮定した場合、例えば、図６のよう
に構成することが可能となる。
【００６４】図６に示すように、車体加減速度演算部３
０は、各車輪の微小ゲインＧd1，Ｇ d2，Ｇd3，Ｇd4のう
ち、前輪の微小ゲインであるＧd1とＧd2との平均値Ｇf
を演算する前輪平均部６６、後輪の微小ゲインであるＧ
d3とＧd4との平均値Ｇr を演算する後輪平均部６７、右
輪の微小ゲインであるＧd1とＧd3との平均値ＧR を演算
する右輪平均部６８、及び左輪の微小ゲインであるＧd2
とＧd4との平均値ＧLを演算する左輪平均部６９を備え
ている。これらの平均部６６，６７，６８，６９は、後
段の処理で割り算があるため、本実施の形態では、平均
化操作の際に１／２を乗じないで、単なる加算器として
構成している。
０は、各車輪の微小ゲインＧd1，Ｇ d2，Ｇd3，Ｇd4のう
ち、前輪の微小ゲインであるＧd1とＧd2との平均値Ｇf
を演算する前輪平均部６６、後輪の微小ゲインであるＧ
d3とＧd4との平均値Ｇr を演算する後輪平均部６７、右
輪の微小ゲインであるＧd1とＧd3との平均値ＧR を演算
する右輪平均部６８、及び左輪の微小ゲインであるＧd2
とＧd4との平均値ＧLを演算する左輪平均部６９を備え
ている。これらの平均部６６，６７，６８，６９は、後
段の処理で割り算があるため、本実施の形態では、平均
化操作の際に１／２を乗じないで、単なる加算器として
構成している。
【００６５】また、前輪平均部６６及び後輪平均部６７
には、前輪平均値Ｇf と後輪平均値Ｇr とを加算する加
算器７２及び前輪平均値Ｇf から後輪平均値Ｇr を減算
する減算器７３が接続されている。また、加算器７２及
び減算器７３の出力端には、減算器７３の演算結果（Ｇ
f −Ｇr ）を加算器７２の演算結果（Ｇf ＋Ｇr ）で除
算する除算器７４が接続されており、さらに除算器７４
の出力端には、除算器７４の演算結果に定数（−Ｌ／２
ｈ）を乗じる定数器７５（図１のハーフカーモデル参
照）が接続されている。なお、この定数器７５の演算結
果が車体の前後方向の加減速度となる。
には、前輪平均値Ｇf と後輪平均値Ｇr とを加算する加
算器７２及び前輪平均値Ｇf から後輪平均値Ｇr を減算
する減算器７３が接続されている。また、加算器７２及
び減算器７３の出力端には、減算器７３の演算結果（Ｇ
f −Ｇr ）を加算器７２の演算結果（Ｇf ＋Ｇr ）で除
算する除算器７４が接続されており、さらに除算器７４
の出力端には、除算器７４の演算結果に定数（−Ｌ／２
ｈ）を乗じる定数器７５（図１のハーフカーモデル参
照）が接続されている。なお、この定数器７５の演算結
果が車体の前後方向の加減速度となる。
【００６６】また、右輪平均部６８及び左輪平均部６９
には、右輪平均値ＧR と左輪平均値ＧL とを加算する加
算器７６及び右輪平均値ＧR から左輪平均値ＧL を減算
する減算器７７が接続されている。また、加算器７６及
び減算器７７の出力端には、減算器７７の演算結果（Ｇ
R −ＧL ）を加算器７６の演算結果（ＧR ＋ＧL ）で除
算する除算器７８が接続されている。さらに除算器７８
の出力端には、除算器７８の演算結果に定数（−Ｌ’／
２ｈ）を乗じる定数器７９（図１のハーフカーモデル参
照）が接続されている。なお、この定数器７９の演算結
果が車体の横方向の加減速度となる。
には、右輪平均値ＧR と左輪平均値ＧL とを加算する加
算器７６及び右輪平均値ＧR から左輪平均値ＧL を減算
する減算器７７が接続されている。また、加算器７６及
び減算器７７の出力端には、減算器７７の演算結果（Ｇ
R −ＧL ）を加算器７６の演算結果（ＧR ＋ＧL ）で除
算する除算器７８が接続されている。さらに除算器７８
の出力端には、除算器７８の演算結果に定数（−Ｌ’／
２ｈ）を乗じる定数器７９（図１のハーフカーモデル参
照）が接続されている。なお、この定数器７９の演算結
果が車体の横方向の加減速度となる。
【００６７】次に、路面状態推定部３４の詳細な構成を
図７を用いて説明する。図７に示すように、本実施の形
態の路面状態推定部３４は、車体の加減速度αの時系列
データと、車輪の加減速度αw の時系列データとの間に
成り立つ関係を最小自乗法に基づいて演算する最小自乗
法演算部８０を備えている。なお、車体の加減速度αと
して、図４の車体加減速度演算部３０により演算された
前後方向の加減速度を用いることができ、また、車輪の
加減速度αw として、車輪加減速度演算部２８により演
算された各車輪の加減速度の平均値を用いることが可能
である。
図７を用いて説明する。図７に示すように、本実施の形
態の路面状態推定部３４は、車体の加減速度αの時系列
データと、車輪の加減速度αw の時系列データとの間に
成り立つ関係を最小自乗法に基づいて演算する最小自乗
法演算部８０を備えている。なお、車体の加減速度αと
して、図４の車体加減速度演算部３０により演算された
前後方向の加減速度を用いることができ、また、車輪の
加減速度αw として、車輪加減速度演算部２８により演
算された各車輪の加減速度の平均値を用いることが可能
である。
【００６８】この最小自乗法演算部８０は、例えば、図
３に示すように車体の加減速度と車輪の加減速度の次元
からなる平面上に各時系列データをプロットしたとき、
各データとの距離の自乗が最小となるような直線を求
め、この直線の勾配ａ及び切片ｂを出力する。この勾配
ａは、αの変化に対するαw の変化の割合（切片の次元
はαw となる）、或いはαw の変化に対するαの変化の
割合（切片の次元はαとなる）のいずれかである。
３に示すように車体の加減速度と車輪の加減速度の次元
からなる平面上に各時系列データをプロットしたとき、
各データとの距離の自乗が最小となるような直線を求
め、この直線の勾配ａ及び切片ｂを出力する。この勾配
ａは、αの変化に対するαw の変化の割合（切片の次元
はαw となる）、或いはαw の変化に対するαの変化の
割合（切片の次元はαとなる）のいずれかである。
【００６９】なお、最小自乗法演算部８０では、得られ
た加減速度のデータに基づいて直線の勾配及び切片を推
定する際に、発進、停止のようにαw が大きく変化する
ところでは、短期間のうちに演算したデータに基づいて
推定し、定常走行に近くαwの変化が小さいところで
は、より多くのデータを集めて推定する。これにより、
αw の大きさに依らず、ほぼ一定の精度で直線勾配等を
推定することが可能となる。
た加減速度のデータに基づいて直線の勾配及び切片を推
定する際に、発進、停止のようにαw が大きく変化する
ところでは、短期間のうちに演算したデータに基づいて
推定し、定常走行に近くαwの変化が小さいところで
は、より多くのデータを集めて推定する。これにより、
αw の大きさに依らず、ほぼ一定の精度で直線勾配等を
推定することが可能となる。
【００７０】また、路面状態推定部３４は、各路面状態
（本実施の形態では、例えば、低μ路、中μ路、高μ路
の３種類とする）毎に、車体の加減速度αと車輪の加減
速度αw と車輪の加減速度の変化率との間に成り立つ関
係を予め記憶した高μ路テーブル８２ａ、中μ路テーブ
ル８２ｂ、及び低μ路テーブル８２ｃを各々備えた変換
部８４ａ、８４ｂ、８４ｃを有している。なお、このよ
うに車輪の加減速度の変化率毎の関係を各テーブルに表
したのは、ドライバの踏力の大きさでブレーキトルクが
変化し、よって、その反力として路面から作用する摩擦
力の係数μが同じ種類のμ路でも変化するからである。
（本実施の形態では、例えば、低μ路、中μ路、高μ路
の３種類とする）毎に、車体の加減速度αと車輪の加減
速度αw と車輪の加減速度の変化率との間に成り立つ関
係を予め記憶した高μ路テーブル８２ａ、中μ路テーブ
ル８２ｂ、及び低μ路テーブル８２ｃを各々備えた変換
部８４ａ、８４ｂ、８４ｃを有している。なお、このよ
うに車輪の加減速度の変化率毎の関係を各テーブルに表
したのは、ドライバの踏力の大きさでブレーキトルクが
変化し、よって、その反力として路面から作用する摩擦
力の係数μが同じ種類のμ路でも変化するからである。
【００７１】これらの変換部８４ａ、８４ｂ、８４ｃ
は、各テーブルに基づいて、入力された車輪の加減速度
の変化率に応じた、車体加減速度と車輪加減速度との関
係（以下、「α−αw 特性」という）を求め、入力され
たα、αw のいずれかの演算値に対応する「α−αw 特
性」での勾配及び切片をそれぞれ変換出力する。
は、各テーブルに基づいて、入力された車輪の加減速度
の変化率に応じた、車体加減速度と車輪加減速度との関
係（以下、「α−αw 特性」という）を求め、入力され
たα、αw のいずれかの演算値に対応する「α−αw 特
性」での勾配及び切片をそれぞれ変換出力する。
【００７２】また、これらの変換部８４ａ、８４ｂ、８
４ｃの後段には、各変換部の出力と最小自乗法演算部８
０の出力とを照合する照合演算部８５が接続されてい
る。この照合演算部８５は、最小自乗法演算部８０によ
り演算された勾配及び切片と、変換部８４ａ、８４ｂ、
８４ｃにより出力された勾配及び切片との間の距離を各
々算出する。
４ｃの後段には、各変換部の出力と最小自乗法演算部８
０の出力とを照合する照合演算部８５が接続されてい
る。この照合演算部８５は、最小自乗法演算部８０によ
り演算された勾配及び切片と、変換部８４ａ、８４ｂ、
８４ｃにより出力された勾配及び切片との間の距離を各
々算出する。
【００７３】さらに、照合演算部８５の後段には、算出
された各距離の中から最小値を選択する最小値選択部８
６が接続されている。この最小値演算部８６は、距離の
最小値を与えた変換部に対応する路面状態の情報コード
を出力する。
された各距離の中から最小値を選択する最小値選択部８
６が接続されている。この最小値演算部８６は、距離の
最小値を与えた変換部に対応する路面状態の情報コード
を出力する。
【００７４】ところで、本実施の形態では、平均ブレー
キ圧Ｐm の回りに車体と車輪と路面とから構成される振
動系の共振周波数ω∞（(3) 式）のブレーキ圧微小振幅
Ｐvを印加する微小励振手段を備えており、該手段によ
りブレーキ圧を微小励振したときの微小ゲインを共振レ
ベルとして演算する。
キ圧Ｐm の回りに車体と車輪と路面とから構成される振
動系の共振周波数ω∞（(3) 式）のブレーキ圧微小振幅
Ｐvを印加する微小励振手段を備えており、該手段によ
りブレーキ圧を微小励振したときの微小ゲインを共振レ
ベルとして演算する。
【００７５】ここで、微小励振手段について説明する。
まず、平均ブレーキ圧指令及び微小励振指令を実際の車
輪への制動トルクに変換する部分（バルブ制御系）は、
図１０に示すように、マスタシリンダ４８、制御バルブ
５２、ホイールシリンダ５６、リザーバー５８及びオイ
ルポンプ６０を備えている。
まず、平均ブレーキ圧指令及び微小励振指令を実際の車
輪への制動トルクに変換する部分（バルブ制御系）は、
図１０に示すように、マスタシリンダ４８、制御バルブ
５２、ホイールシリンダ５６、リザーバー５８及びオイ
ルポンプ６０を備えている。
【００７６】このうちブレーキペダル４６は、ブレーキ
ペダル４６の踏力に応じて増圧するマスタシリンダ４８
を介して制御バルブ５２の増圧バルブ５０へ接続されて
いる。また、制御バルブ５２は、減圧バルブ５４を介し
て低圧源としてのリザーバー５８へ接続されている。さ
らに、制御バルブ５２には、該制御バルブによって供給
されたブレーキ圧をブレーキディスクに加えるためのホ
イールシリンダ５６が接続されている。この制御バルブ
５２は、入力されたバルブ動作指令に基づいて増圧バル
ブ５０及び減圧バルブ５４の開閉を制御する。
ペダル４６の踏力に応じて増圧するマスタシリンダ４８
を介して制御バルブ５２の増圧バルブ５０へ接続されて
いる。また、制御バルブ５２は、減圧バルブ５４を介し
て低圧源としてのリザーバー５８へ接続されている。さ
らに、制御バルブ５２には、該制御バルブによって供給
されたブレーキ圧をブレーキディスクに加えるためのホ
イールシリンダ５６が接続されている。この制御バルブ
５２は、入力されたバルブ動作指令に基づいて増圧バル
ブ５０及び減圧バルブ５４の開閉を制御する。
【００７７】なお、この制御バルブ５２が増圧バルブ５
０のみを開くように制御されると、ホイールシリンダ５
６の油圧（ホイールシリンダ圧）は、ドライバがブレー
キペダル４６を踏み込むことによって得られる圧力に比
例したマスタシリンダ４８の油圧（マスタシリンダ圧）
まで上昇する。逆に減圧バルブ５４のみを開くように制
御されると、ホイールシリンダ圧は、ほぼ大気圧のリザ
ーバ５８の圧力（リザーバ圧）まで減少する。また、両
方のバルブを閉じるように制御されると、ホイールシリ
ンダ圧は保持される。
０のみを開くように制御されると、ホイールシリンダ５
６の油圧（ホイールシリンダ圧）は、ドライバがブレー
キペダル４６を踏み込むことによって得られる圧力に比
例したマスタシリンダ４８の油圧（マスタシリンダ圧）
まで上昇する。逆に減圧バルブ５４のみを開くように制
御されると、ホイールシリンダ圧は、ほぼ大気圧のリザ
ーバ５８の圧力（リザーバ圧）まで減少する。また、両
方のバルブを閉じるように制御されると、ホイールシリ
ンダ圧は保持される。
【００７８】ホイールシリンダ５６によりブレーキディ
スクに加えられるブレーキ力（ホイールシリンダ圧に相
当）は、マスタシリンダ４８の高油圧が供給される増圧
時間、リザーバー５８の低油圧が供給される減圧時間、
及び供給油圧が保持される保持時間の比率と、圧力セン
サ等により検出されたマスタシリンダ圧及びリザーバー
圧とから求められる。
スクに加えられるブレーキ力（ホイールシリンダ圧に相
当）は、マスタシリンダ４８の高油圧が供給される増圧
時間、リザーバー５８の低油圧が供給される減圧時間、
及び供給油圧が保持される保持時間の比率と、圧力セン
サ等により検出されたマスタシリンダ圧及びリザーバー
圧とから求められる。
【００７９】従って、制御バルブ５２の増減圧時間をマ
スタシリンダ圧に応じて制御することにより、所望のブ
レーキトルクを実現することができる。そして、ブレー
キ圧の微小励振は、平均ブレーキ力を実現する制御バル
ブ５２の増減圧制御と同時に共振周波数に対応した周期
で増圧減圧制御を行うことにより可能となる。
スタシリンダ圧に応じて制御することにより、所望のブ
レーキトルクを実現することができる。そして、ブレー
キ圧の微小励振は、平均ブレーキ力を実現する制御バル
ブ５２の増減圧制御と同時に共振周波数に対応した周期
で増圧減圧制御を行うことにより可能となる。
【００８０】具体的な制御の内容として、図１１に示す
ように、微小励振の周期（例えば２４［ｍｓ］）の半周
期Ｔ／２毎に増圧と減圧のそれぞれのモードを切り替
え、バルブへの増減圧指令は、モード切り替えの瞬間か
ら増圧時間ｔi 、減圧時間ｔrのそれぞれの時間分だけ
増圧・減圧指令を出力し、残りの時間は、保持指令を出
力する。平均ブレーキ力は、マスタシリンダ圧に応じた
増圧時間ｔi と減圧時間ｔr との比によって定まると共
に、共振周波数に対応した半周期Ｔ／２毎の増圧・減圧
モードの切り替えによって、平均ブレーキ力の回りに微
小振動が印加される。
ように、微小励振の周期（例えば２４［ｍｓ］）の半周
期Ｔ／２毎に増圧と減圧のそれぞれのモードを切り替
え、バルブへの増減圧指令は、モード切り替えの瞬間か
ら増圧時間ｔi 、減圧時間ｔrのそれぞれの時間分だけ
増圧・減圧指令を出力し、残りの時間は、保持指令を出
力する。平均ブレーキ力は、マスタシリンダ圧に応じた
増圧時間ｔi と減圧時間ｔr との比によって定まると共
に、共振周波数に対応した半周期Ｔ／２毎の増圧・減圧
モードの切り替えによって、平均ブレーキ力の回りに微
小振動が印加される。
【００８１】ここで、この微小励振手段によりブレーキ
圧が微小励振されたときの共振レベルを表す微小ゲイン
を演算する微小ゲイン演算部２２の詳細な構成を図８を
用いて説明する。図８に示すように、微小ゲイン演算部
２２は、平均ブレーキ圧の回りに共振周波数ω∞でブレ
ーキ圧を微小励振したときの、車輪速度ωi （ｉは車輪
番号、ｉ＝1,2,3,4 ）の共振周波数ω∞の微小振幅（車
輪速微小振幅ωwv）を検出する車輪速微小振幅検出部４
０と、共振周波数ω∞のブレーキ圧の微小振幅Ｐv を検
出するブレーキ圧微小振幅検出部４２と、検出された車
輪速微小振幅ω wvをブレーキ圧微小振幅Ｐv で除算する
ことにより微小ゲインＧd を出力する除算器４４と、か
ら構成される。ここで、車輪速微小振幅検出部４０は、
共振周波数ω∞の振動成分を抽出するフィルタ処理を行
う図９のような演算部として実現できる。例えば、この
振動系の共振周波数ω∞が４０［Ｈｚ］程度であるの
で、制御性を考慮して１周期を２４［ｍｓ］、約４１．
７［Ｈｚ］に取り、この周波数を中心周波数とする帯域
通過フィルタを設ける。このフィルタにより、車輪速度
信号ωi から約４１．７［Ｈｚ］近傍の周波数成分のみ
が抽出される。さらに、このフィルタ出力を全波整流器
により全波整流、直流平滑化し、この直流平滑化信号か
ら低域通過フィルタによって低域振動成分のみを通過さ
せることにより、車輪速微小振幅を出力する。また、周
期の整数倍、例えば１周期の２４［ｍｓ］、２周期の４
８［ｍｓ］の時系列データを連続的に取り込み、４１．
７［Ｈｚ］の単位正弦波、単位余弦波との相関を求める
ことによっても実現できる。
圧が微小励振されたときの共振レベルを表す微小ゲイン
を演算する微小ゲイン演算部２２の詳細な構成を図８を
用いて説明する。図８に示すように、微小ゲイン演算部
２２は、平均ブレーキ圧の回りに共振周波数ω∞でブレ
ーキ圧を微小励振したときの、車輪速度ωi （ｉは車輪
番号、ｉ＝1,2,3,4 ）の共振周波数ω∞の微小振幅（車
輪速微小振幅ωwv）を検出する車輪速微小振幅検出部４
０と、共振周波数ω∞のブレーキ圧の微小振幅Ｐv を検
出するブレーキ圧微小振幅検出部４２と、検出された車
輪速微小振幅ω wvをブレーキ圧微小振幅Ｐv で除算する
ことにより微小ゲインＧd を出力する除算器４４と、か
ら構成される。ここで、車輪速微小振幅検出部４０は、
共振周波数ω∞の振動成分を抽出するフィルタ処理を行
う図９のような演算部として実現できる。例えば、この
振動系の共振周波数ω∞が４０［Ｈｚ］程度であるの
で、制御性を考慮して１周期を２４［ｍｓ］、約４１．
７［Ｈｚ］に取り、この周波数を中心周波数とする帯域
通過フィルタを設ける。このフィルタにより、車輪速度
信号ωi から約４１．７［Ｈｚ］近傍の周波数成分のみ
が抽出される。さらに、このフィルタ出力を全波整流器
により全波整流、直流平滑化し、この直流平滑化信号か
ら低域通過フィルタによって低域振動成分のみを通過さ
せることにより、車輪速微小振幅を出力する。また、周
期の整数倍、例えば１周期の２４［ｍｓ］、２周期の４
８［ｍｓ］の時系列データを連続的に取り込み、４１．
７［Ｈｚ］の単位正弦波、単位余弦波との相関を求める
ことによっても実現できる。
【００８２】なお、ブレーキ圧微小振幅Ｐv は、マスタ
シリンダ圧、図１１に示したバルブの増圧時間ｔi の長
さ、及び減圧時間ｔr の長さによって所定の関係で定ま
るので、図８のブレーキ圧微小振幅検出部４２は、マス
タシリンダ圧、増圧時間ｔi及び減圧時間ｔr からブレ
ーキ圧微小振幅Ｐv を出力するテーブルとして構成する
ことができる。
シリンダ圧、図１１に示したバルブの増圧時間ｔi の長
さ、及び減圧時間ｔr の長さによって所定の関係で定ま
るので、図８のブレーキ圧微小振幅検出部４２は、マス
タシリンダ圧、増圧時間ｔi及び減圧時間ｔr からブレ
ーキ圧微小振幅Ｐv を出力するテーブルとして構成する
ことができる。
【００８３】次に、本実施の形態の作用を説明する。所
定条件の成立時に、微小励振手段が共振周波数ω∞でブ
レーキ圧を微小励振する。ここで、所定条件とは、車体
が加減速度運動をするときの条件であり、例えば、ブレ
ーキペダルが踏み込まれたとき、操舵角が所定角度を超
えたとき、アクセル開度が所定値を超えたときなどが挙
げられる。ここで、図１０のブレーキ圧制御系におい
て、アクセル開度が所定値を超えたときなどのようにブ
レーキペダルが踏み込まれていない状態でも、踏力とは
無関係に増圧バルブ５０を介して制御バルブ５２へ高油
圧を加える手段を設けることにより、ブレーキ圧の微小
励振が可能となる。
定条件の成立時に、微小励振手段が共振周波数ω∞でブ
レーキ圧を微小励振する。ここで、所定条件とは、車体
が加減速度運動をするときの条件であり、例えば、ブレ
ーキペダルが踏み込まれたとき、操舵角が所定角度を超
えたとき、アクセル開度が所定値を超えたときなどが挙
げられる。ここで、図１０のブレーキ圧制御系におい
て、アクセル開度が所定値を超えたときなどのようにブ
レーキペダルが踏み込まれていない状態でも、踏力とは
無関係に増圧バルブ５０を介して制御バルブ５２へ高油
圧を加える手段を設けることにより、ブレーキ圧の微小
励振が可能となる。
【００８４】ブレーキ圧が微小励振されると、図４の微
小ゲイン演算部２２ａ，２２ｂ，２２ｃ，２２ｄが、各
車輪の車輪速度ω1 ，ω2 ，ω3 ，ω4 から微小ゲイン
Ｇd1，Ｇd2，Ｇd3，Ｇd4を演算し、荷重配分推定部２６
へ出力する。ここで、右前輪及び左前輪の微小ゲイン
は、遅延部２４ａ，２４ｂにより、車体速度Ｖに応じた
遅延時間Ｔだけ遅延されて荷重配分推定部２６へ伝達さ
れる。なお、前後の荷重配分を求める場合、この遅延時
間Ｔを、例えば、次式により演算することができる。
小ゲイン演算部２２ａ，２２ｂ，２２ｃ，２２ｄが、各
車輪の車輪速度ω1 ，ω2 ，ω3 ，ω4 から微小ゲイン
Ｇd1，Ｇd2，Ｇd3，Ｇd4を演算し、荷重配分推定部２６
へ出力する。ここで、右前輪及び左前輪の微小ゲイン
は、遅延部２４ａ，２４ｂにより、車体速度Ｖに応じた
遅延時間Ｔだけ遅延されて荷重配分推定部２６へ伝達さ
れる。なお、前後の荷重配分を求める場合、この遅延時
間Ｔを、例えば、次式により演算することができる。
【００８５】Ｔ ＝ Ｖ／Ｌ ここで、Ｌは図１に示すように前後車輪間の間隔であ
る。すなわち、この時間Ｔの遅延によって、前後の車輪
で同じ路面位置（同じμ値）での微小ゲインの比較が可
能となる。
る。すなわち、この時間Ｔの遅延によって、前後の車輪
で同じ路面位置（同じμ値）での微小ゲインの比較が可
能となる。
【００８６】次に、荷重配分推定部２６において、総和
器３６が４輪の微小ゲインをすべて加算することにより
総和Ｇw （＝Ｇd1＋Ｇd2＋Ｇd3＋Ｇd4）を出力し、乗算
器３８が総和Ｇw に１／４を乗じた値Ｇw ／４を算出す
る。なお、微小ゲインが輪荷重に比例すると仮定した場
合、総和Ｇw は車両重量Ｗに対応する共振レベルを表
し、よって、この乗算器３８の出力Ｇw ／４は、１つの
車輪の平均的な荷重に対応する共振レベルを表すことと
なる。
器３６が４輪の微小ゲインをすべて加算することにより
総和Ｇw （＝Ｇd1＋Ｇd2＋Ｇd3＋Ｇd4）を出力し、乗算
器３８が総和Ｇw に１／４を乗じた値Ｇw ／４を算出す
る。なお、微小ゲインが輪荷重に比例すると仮定した場
合、総和Ｇw は車両重量Ｗに対応する共振レベルを表
し、よって、この乗算器３８の出力Ｇw ／４は、１つの
車輪の平均的な荷重に対応する共振レベルを表すことと
なる。
【００８７】すなわち、本実施の形態では、各車輪の静
止荷重がほぼ均等の場合を想定している。ただし、各車
輪の静止荷重が均等でない場合でも、総和Ｇw を各車輪
の静止荷重の比に応じて振り分けることで対応できる。
止荷重がほぼ均等の場合を想定している。ただし、各車
輪の静止荷重が均等でない場合でも、総和Ｇw を各車輪
の静止荷重の比に応じて振り分けることで対応できる。
【００８８】そして、各差分器６２ａ、６２ｂ、６２
ｃ、６２ｄは、各微小ゲインと乗算器３８の出力Ｇw ／
４との差分｜Ｇd1−Ｇw ／４｜，｜Ｇd2−Ｇw ／４｜，
｜Ｇd3−Ｇw ／４｜，｜Ｇd4−Ｇw ／４｜を各々算出す
る。これらの差分は、荷重移動による共振レベルの増加
分或いは減少分を表している。
ｃ、６２ｄは、各微小ゲインと乗算器３８の出力Ｇw ／
４との差分｜Ｇd1−Ｇw ／４｜，｜Ｇd2−Ｇw ／４｜，
｜Ｇd3−Ｇw ／４｜，｜Ｇd4−Ｇw ／４｜を各々算出す
る。これらの差分は、荷重移動による共振レベルの増加
分或いは減少分を表している。
【００８９】各差分器により差分が算出されると、荷重
移動比率演算器６５は、前後の荷重移動配分比率Ｐ1 を
次式に従って演算する。
移動比率演算器６５は、前後の荷重移動配分比率Ｐ1 を
次式に従って演算する。
【００９０】
【数９】
【００９１】
【数１０】
【００９２】すなわち、右輪の差分の和と、左輪の差分
の和との比率である。ところで、車両がコーナリングす
るときの荷重移動の前後配分が、例えば６０：４０のと
き、各輪の荷重は、 前輪外輪＝Ｗ／４＋０．６ΔＷ 後輪外輪＝Ｗ／４＋０．４ΔＷ 前輪内輪＝Ｗ／４−０．６ΔＷ 後輪内輪＝Ｗ／４−０．４ΔＷ となる。ただし、Ｗは車両重量、ΔＷは、コーナリング
中に車両にロールを発生させる遠心力によるモーメント
Ｎに起因する左右の荷重移動である（ΔＷ＝Ｎ／Ｌ’，
Ｌ’は左右輪間の間隔）。よって、微小ゲインが輪荷重
に比例すると仮定した場合、(18)式により、前後配分の
比率（６０：４０）が直ちに求められることがわかる。
の和との比率である。ところで、車両がコーナリングす
るときの荷重移動の前後配分が、例えば６０：４０のと
き、各輪の荷重は、 前輪外輪＝Ｗ／４＋０．６ΔＷ 後輪外輪＝Ｗ／４＋０．４ΔＷ 前輪内輪＝Ｗ／４−０．６ΔＷ 後輪内輪＝Ｗ／４−０．４ΔＷ となる。ただし、Ｗは車両重量、ΔＷは、コーナリング
中に車両にロールを発生させる遠心力によるモーメント
Ｎに起因する左右の荷重移動である（ΔＷ＝Ｎ／Ｌ’，
Ｌ’は左右輪間の間隔）。よって、微小ゲインが輪荷重
に比例すると仮定した場合、(18)式により、前後配分の
比率（６０：４０）が直ちに求められることがわかる。
【００９３】このコーナリング時の前後の荷重配分は、
前後方向の加速度及び横加速度だけからは求めることが
できず、車体特性、コーナリング時のサスペンションの
状態量などにより定まるが、本実施の形態では、微小ゲ
インに基づいて直接的に荷重配分を求めるので、該荷重
配分のより正確な値を得ることができる。また、このこ
とは、加速時、ブレーキ制動時における前後及び左右の
荷重移動の場合においても成立する。
前後方向の加速度及び横加速度だけからは求めることが
できず、車体特性、コーナリング時のサスペンションの
状態量などにより定まるが、本実施の形態では、微小ゲ
インに基づいて直接的に荷重配分を求めるので、該荷重
配分のより正確な値を得ることができる。また、このこ
とは、加速時、ブレーキ制動時における前後及び左右の
荷重移動の場合においても成立する。
【００９４】さらに、本実施の形態では、乗算器６４
ａ、６４ｂ、６４ｃ、６４ｄが、(5)式に基づき各車輪
の微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に比例定数ｋを乗
算することにより、それぞれ右前輪の荷重Ｗ1 、左前輪
の荷重Ｗ2 、右後輪の荷重Ｗ3、左後輪の荷重Ｗ4 を直
ちに算出する。
ａ、６４ｂ、６４ｃ、６４ｄが、(5)式に基づき各車輪
の微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に比例定数ｋを乗
算することにより、それぞれ右前輪の荷重Ｗ1 、左前輪
の荷重Ｗ2 、右後輪の荷重Ｗ3、左後輪の荷重Ｗ4 を直
ちに算出する。
【００９５】これらの輪荷重は、加速度等から間接的に
推定されるのではなく、路面に実際に接地している車輪
の共振特性を表す微小ゲインにより演算されるので、本
実施の形態によれば、より正確な輪荷重値を得ることが
できる。
推定されるのではなく、路面に実際に接地している車輪
の共振特性を表す微小ゲインにより演算されるので、本
実施の形態によれば、より正確な輪荷重値を得ることが
できる。
【００９６】なお、微小ゲインの値は、荷重のみならず
路面μの値に依存するが、遅延部２４ａ、２４ｂによ
り、前輪の微小ゲインを車体速度に応じて遅延させるこ
とにより、ほぼ同じ路面位置（同じ路面μ）での前輪と
後輪の微小ゲインを用いて荷重配分や輪荷重を算出する
ので、車輪間の路面状態の相違による演算結果の誤差を
少なくすることができる。
路面μの値に依存するが、遅延部２４ａ、２４ｂによ
り、前輪の微小ゲインを車体速度に応じて遅延させるこ
とにより、ほぼ同じ路面位置（同じ路面μ）での前輪と
後輪の微小ゲインを用いて荷重配分や輪荷重を算出する
ので、車輪間の路面状態の相違による演算結果の誤差を
少なくすることができる。
【００９７】以上のように本実施の形態に係る荷重配分
推定部２６では、荷重配分や輪荷重を精度良く推定する
ことができるので、これらのパラメータを用いた制御を
より正確かつ安定に行うことができる。例えば、アクテ
ィブサスペンションにより荷重配分を所定範囲になるよ
うに制御する車両姿勢制御、荷重が多い車輪に大きい制
動力を加えることによる車両安定制御、及び自動変速機
の変速点の変更制御などに有効に適用することができ
る。さらに、演算した車輪の荷重を用いてタイヤのコー
ナリング特性が限界を超えないように制御することによ
り、スピンを防止するアンチスピン制御にも応用するこ
とができる。
推定部２６では、荷重配分や輪荷重を精度良く推定する
ことができるので、これらのパラメータを用いた制御を
より正確かつ安定に行うことができる。例えば、アクテ
ィブサスペンションにより荷重配分を所定範囲になるよ
うに制御する車両姿勢制御、荷重が多い車輪に大きい制
動力を加えることによる車両安定制御、及び自動変速機
の変速点の変更制御などに有効に適用することができ
る。さらに、演算した車輪の荷重を用いてタイヤのコー
ナリング特性が限界を超えないように制御することによ
り、スピンを防止するアンチスピン制御にも応用するこ
とができる。
【００９８】また、本実施の形態に係る車体加速度演算
部３０は、図６に示す構成によって、演算された各車輪
の微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に対し、(14)式及
び(15)式に相当する演算を実行することにより、前後方
向の車体加減速度及び横方向の車体加減速度を算出す
る。この車体加減速度は、次に述べる路面状態推定部３
４で用いられる以外に、車両姿勢制御などに有効に用い
ることができる。
部３０は、図６に示す構成によって、演算された各車輪
の微小ゲインＧd1，Ｇd2，Ｇd3，Ｇd4に対し、(14)式及
び(15)式に相当する演算を実行することにより、前後方
向の車体加減速度及び横方向の車体加減速度を算出す
る。この車体加減速度は、次に述べる路面状態推定部３
４で用いられる以外に、車両姿勢制御などに有効に用い
ることができる。
【００９９】さらに、本実施の形態に係る路面状態推定
部３４では、図７の最小自乗法演算部８０が、車体加減
速度演算部３０により演算された車体加減速度αの時系
列データと、車輪加減速度演算部２８により演算された
車輪加減速度αw の時系列データとの間に成り立つ関係
を、誤差の自乗和が最小となるような直線で近似し、該
直線の勾配ａ0 及び切片ｂ0 を求める。
部３４では、図７の最小自乗法演算部８０が、車体加減
速度演算部３０により演算された車体加減速度αの時系
列データと、車輪加減速度演算部２８により演算された
車輪加減速度αw の時系列データとの間に成り立つ関係
を、誤差の自乗和が最小となるような直線で近似し、該
直線の勾配ａ0 及び切片ｂ0 を求める。
【０１００】ここで、実際の路面におけるα−αw 特性
のグラフを図１２に示す。同図の横軸は重力加速度ｇに
対する車体減速度の比、縦軸は重力加速度ｇに対する車
輪減速度の比を表している。また、この図は、車体減速
度比０の状態からブレーキ制動を開始し、車体減速度比
が約１．６３のところでブレーキオフし、その後車体減
速度比が０になるまでのα−αw 特性の時間的変化を、
矢印に沿って高μ路（０．８）、中μ路（０．４）、低
μ路（０．２）毎に示したものである。
のグラフを図１２に示す。同図の横軸は重力加速度ｇに
対する車体減速度の比、縦軸は重力加速度ｇに対する車
輪減速度の比を表している。また、この図は、車体減速
度比０の状態からブレーキ制動を開始し、車体減速度比
が約１．６３のところでブレーキオフし、その後車体減
速度比が０になるまでのα−αw 特性の時間的変化を、
矢印に沿って高μ路（０．８）、中μ路（０．４）、低
μ路（０．２）毎に示したものである。
【０１０１】図１２に示すように、路面状態が異なる
と、α−αw 特性のグラフの勾配及び切片の値が異なっ
ており、これらの値から路面状態が識別できることがわ
かる。また、ブレーキオフ以降では、同じμ路でも、勾
配等が異なっており、ブレーキトルクの大きさ（本実施
の形態では、車輪加減速度の変化率等で表す）に応じて
判断基準を変更する必要があることがわかる。
と、α−αw 特性のグラフの勾配及び切片の値が異なっ
ており、これらの値から路面状態が識別できることがわ
かる。また、ブレーキオフ以降では、同じμ路でも、勾
配等が異なっており、ブレーキトルクの大きさ（本実施
の形態では、車輪加減速度の変化率等で表す）に応じて
判断基準を変更する必要があることがわかる。
【０１０２】そこで、図７の変換部８４ａ、８４ｂ、８
４ｃは、高μ路テーブル８２ａ，中μ路テーブル８２
ｂ，低μ路テーブル８２ｃから、加減速度変化率演算部
３２により演算された車輪加減速度の変化率に対応する
α−αw 特性を抽出する。そして、抽出されたα−αw
特性から、入力された車体加減速度の演算値（α1 とす
る）及び車輪加減速度の演算値（αw 1 とする）のいず
れかに対応する勾配及び切片を求める。
４ｃは、高μ路テーブル８２ａ，中μ路テーブル８２
ｂ，低μ路テーブル８２ｃから、加減速度変化率演算部
３２により演算された車輪加減速度の変化率に対応する
α−αw 特性を抽出する。そして、抽出されたα−αw
特性から、入力された車体加減速度の演算値（α1 とす
る）及び車輪加減速度の演算値（αw 1 とする）のいず
れかに対応する勾配及び切片を求める。
【０１０３】例えば、入力されたα1 （又はαw 1 ）の
前後の２つの値と、該値に対応する２つの車輪加減速度
の値（又は車体加減速度）をα−αw 特性から求め、こ
れら２つの座標点を結ぶ直線の勾配及び切片を求める。
なお、最小自乗法演算部８０が複数の加減速度のデータ
を用いているので、これらのデータに対応する複数の点
を、抽出されたα−αw 特性上で求め、該複数点に基づ
いて勾配及び切片を求めても良い。
前後の２つの値と、該値に対応する２つの車輪加減速度
の値（又は車体加減速度）をα−αw 特性から求め、こ
れら２つの座標点を結ぶ直線の勾配及び切片を求める。
なお、最小自乗法演算部８０が複数の加減速度のデータ
を用いているので、これらのデータに対応する複数の点
を、抽出されたα−αw 特性上で求め、該複数点に基づ
いて勾配及び切片を求めても良い。
【０１０４】また、処理速度を向上させるために、各テ
ーブル８２ａ、８２ｂ、８２ｃを、車輪の加減速度の変
化率及び車体加減速度（又は車輪加減速度）に対して、
直接、勾配及び切片のデータを与えるテーブルとして構
成しても良く、この場合には、各変換部は、入力された
車輪の加減速度の変化率及び車体加減速度（又は車輪加
減速度）に基づき直ちに勾配及び切片に変換することが
できる。
ーブル８２ａ、８２ｂ、８２ｃを、車輪の加減速度の変
化率及び車体加減速度（又は車輪加減速度）に対して、
直接、勾配及び切片のデータを与えるテーブルとして構
成しても良く、この場合には、各変換部は、入力された
車輪の加減速度の変化率及び車体加減速度（又は車輪加
減速度）に基づき直ちに勾配及び切片に変換することが
できる。
【０１０５】このようにして変換部８４ａ，８４ｂ，８
４ｃが、勾配ａ及び切片ｂのデータを変換出力すると、
照合演算部８５では、これらのデータと最小自乗法演算
部８０の出力（ａ0 ，ｂ0 ）とをそれぞれ照合すること
により、各変換部に対し次式のような距離値Ｄ1 ，
Ｄ2 ，Ｄ3 を各々算出する。ただし、変換部８４ａ，８
４ｂ，８４ｃの出力をそれぞれ（ａ1 ，ｂ1 ），
（ａ2 ，ｂ2 ），（ａ3 ，ｂ3 ）とする。
４ｃが、勾配ａ及び切片ｂのデータを変換出力すると、
照合演算部８５では、これらのデータと最小自乗法演算
部８０の出力（ａ0 ，ｂ0 ）とをそれぞれ照合すること
により、各変換部に対し次式のような距離値Ｄ1 ，
Ｄ2 ，Ｄ3 を各々算出する。ただし、変換部８４ａ，８
４ｂ，８４ｃの出力をそれぞれ（ａ1 ，ｂ1 ），
（ａ2 ，ｂ2 ），（ａ3 ，ｂ3 ）とする。
【０１０６】 Ｄ1 ＝｛（ａ0 −ａ1 ）2 ＋（ｂ0 −ｂ1 ）2 ｝ Ｄ2 ＝｛（ａ0 −ａ2 ）2 ＋（ｂ0 −ｂ2 ）2 ｝ (20) Ｄ3 ＝｛（ａ0 −ａ3 ）2 ＋（ｂ0 −ｂ3 ）2 ｝ 次に、最小値選択部８６では、(20)式により演算された
各距離値Ｄ1 ，Ｄ2 ，Ｄ3 の中から最小値を選択し、こ
の最小値を与えた変換部に対応する路面状態の情報コー
ドを出力する。例えば、Ｄ2 が最小値の場合、変換部８
４ｂに対応する路面状態である”中μ路”の情報コード
を出力する。
各距離値Ｄ1 ，Ｄ2 ，Ｄ3 の中から最小値を選択し、こ
の最小値を与えた変換部に対応する路面状態の情報コー
ドを出力する。例えば、Ｄ2 が最小値の場合、変換部８
４ｂに対応する路面状態である”中μ路”の情報コード
を出力する。
【０１０７】なお、(20)式の距離値Ｄ1 〜Ｄ3 は、実際
に検出されたαとαw との間の関係と、各路面状態毎に
モデル化されたαとαw との間の関係とが、どれだけ離
れているかを示すパラメータであり、従って、この距離
値が小さいということは、該距離値を与えたモデルに近
い路面状態であるということを示している。
に検出されたαとαw との間の関係と、各路面状態毎に
モデル化されたαとαw との間の関係とが、どれだけ離
れているかを示すパラメータであり、従って、この距離
値が小さいということは、該距離値を与えたモデルに近
い路面状態であるということを示している。
【０１０８】以上のように本実施の形態に係る路面状態
推定部３４は、現在走行中の路面状態を、正確に判断す
ることができる。また、既に述べた荷重配分推定部２６
による荷重配分等の推定結果や車体加減速度の演算結果
を用いた車両安定化制御において、路面状態推定部３４
により推定された路面状態に応じて制御量や目標値等を
変更することにより、路面状態に係わらず常に安定な制
御を行うことが可能となる。 （第２の実施の形態）次に、第２の実施の形態を説明す
る。なお、第２の実施の形態の構成は、第１の実施の形
態の構成と同様であるので、同一の符号を付して詳細な
説明を省略し、異なる構成部分についてのみ説明する。
推定部３４は、現在走行中の路面状態を、正確に判断す
ることができる。また、既に述べた荷重配分推定部２６
による荷重配分等の推定結果や車体加減速度の演算結果
を用いた車両安定化制御において、路面状態推定部３４
により推定された路面状態に応じて制御量や目標値等を
変更することにより、路面状態に係わらず常に安定な制
御を行うことが可能となる。 （第２の実施の形態）次に、第２の実施の形態を説明す
る。なお、第２の実施の形態の構成は、第１の実施の形
態の構成と同様であるので、同一の符号を付して詳細な
説明を省略し、異なる構成部分についてのみ説明する。
【０１０９】第１の実施の形態では、微小ゲインと輪荷
重とが略比例関係にあると仮定したが、より厳密には、
正の相関関係の関数でモデル化することができる。そこ
で、第２の実施の形態では、図１３に示すように、演算
された微小ゲインＧd を輪荷重に変換する変換部９０を
各車輪毎に設ける。
重とが略比例関係にあると仮定したが、より厳密には、
正の相関関係の関数でモデル化することができる。そこ
で、第２の実施の形態では、図１３に示すように、演算
された微小ゲインＧd を輪荷重に変換する変換部９０を
各車輪毎に設ける。
【０１１０】図１３の変換部９０は、少なくとも微小ゲ
インＧd と輪荷重との間の正の相関関係を示す関数９２
を記憶したテーブル９１を備えている。なお、微小ゲイ
ンＧ d と輪荷重以外の１又は２以上の他のパラメータに
応じて関数９２が変化する場合、このテーブル９１を、
当該他のパラメータ毎に微小ゲインＧd と輪荷重との間
の関数９２を記憶する３次元以上のテーブルとして構成
することができる。
インＧd と輪荷重との間の正の相関関係を示す関数９２
を記憶したテーブル９１を備えている。なお、微小ゲイ
ンＧ d と輪荷重以外の１又は２以上の他のパラメータに
応じて関数９２が変化する場合、このテーブル９１を、
当該他のパラメータ毎に微小ゲインＧd と輪荷重との間
の関数９２を記憶する３次元以上のテーブルとして構成
することができる。
【０１１１】なお、各車輪毎の変換部９０は、図４の微
小ゲイン演算部２２ａ，２２ｂ，２２ｃ，２２ｄと、荷
重配分推定部２６（及び車体加減速度演算部３０）との
間にそれぞれ介在される。ただし、前輪の変換部に関し
ては遅延部２４ａ、２４ｂの前段及び後段のいずれでも
良い。
小ゲイン演算部２２ａ，２２ｂ，２２ｃ，２２ｄと、荷
重配分推定部２６（及び車体加減速度演算部３０）との
間にそれぞれ介在される。ただし、前輪の変換部に関し
ては遅延部２４ａ、２４ｂの前段及び後段のいずれでも
良い。
【０１１２】次に、第２の実施の形態の作用を説明す
る。第ｉ輪（i=1,2,3,4 ）の変換部９０に、第ｉ輪の微
小ゲインＧdiが入力されると、変換部９０は、テーブル
９１により示された関数９２により、Ｇdiを第ｉ輪の荷
重Ｗi に変換して出力する。そして、後段に接続されて
いる荷重配分推定部２６及び車体加減速度演算部３０で
は、微小ゲインＧdiを輪荷重Ｗi に置き換えた演算を行
う。この輪荷重Ｗi は、微小ゲインと輪荷重とが略比例
関係にあると仮定した時よりも精度の良い値として求め
られているので、第２の実施の形態では、第１の実施の
形態より正確に荷重配分や車体加減速度を演算すること
ができる。
る。第ｉ輪（i=1,2,3,4 ）の変換部９０に、第ｉ輪の微
小ゲインＧdiが入力されると、変換部９０は、テーブル
９１により示された関数９２により、Ｇdiを第ｉ輪の荷
重Ｗi に変換して出力する。そして、後段に接続されて
いる荷重配分推定部２６及び車体加減速度演算部３０で
は、微小ゲインＧdiを輪荷重Ｗi に置き換えた演算を行
う。この輪荷重Ｗi は、微小ゲインと輪荷重とが略比例
関係にあると仮定した時よりも精度の良い値として求め
られているので、第２の実施の形態では、第１の実施の
形態より正確に荷重配分や車体加減速度を演算すること
ができる。
【０１１３】また、第２の実施の形態では、変換部９０
により直ちに輪荷重が求められるので、図５の乗算器６
４ａ，６４ｂ，６４ｃ，６４ｄを省略することができ
る。
により直ちに輪荷重が求められるので、図５の乗算器６
４ａ，６４ｂ，６４ｃ，６４ｄを省略することができ
る。
【０１１４】なお、他のパラメータ（操舵角、路面勾配
等）を考慮する場合、変換部９０は、当該他のパラメー
タが入力されると、該パラメータを１つの次元として有
するテーブル９１より、この入力値に対応する微小ゲイ
ンと輪荷重との関係を抽出し、抽出された関係に基づい
て、入力された微小ゲインＧdiを第ｉ輪の荷重Ｗi に変
換して出力する。この場合、さらに正確に荷重配分や車
体加減速度を演算することができる。
等）を考慮する場合、変換部９０は、当該他のパラメー
タが入力されると、該パラメータを１つの次元として有
するテーブル９１より、この入力値に対応する微小ゲイ
ンと輪荷重との関係を抽出し、抽出された関係に基づい
て、入力された微小ゲインＧdiを第ｉ輪の荷重Ｗi に変
換して出力する。この場合、さらに正確に荷重配分や車
体加減速度を演算することができる。
【０１１５】以上が本発明の各実施の形態であるが、本
発明は、その要旨を逸脱しない範囲において種々に変更
可能である。
発明は、その要旨を逸脱しない範囲において種々に変更
可能である。
【０１１６】例えば、第１の実施の形態では、共振レベ
ルとして微小ゲインを用いたが、すべての車輪のブレー
キ圧微小振幅が同じとなる条件では、ブレーキ圧微小振
幅値で正規化する必要はなく、車輪速微小振幅そのもの
を共振レベルとして用いることができる。この場合、図
４の構成において、微小ゲイン演算部２２ａ、２２ｂ、
２２ｃ、２２ｄを、図８の車輪速微小振幅検出部４０に
置き換えることができる。
ルとして微小ゲインを用いたが、すべての車輪のブレー
キ圧微小振幅が同じとなる条件では、ブレーキ圧微小振
幅値で正規化する必要はなく、車輪速微小振幅そのもの
を共振レベルとして用いることができる。この場合、図
４の構成において、微小ゲイン演算部２２ａ、２２ｂ、
２２ｃ、２２ｄを、図８の車輪速微小振幅検出部４０に
置き換えることができる。
【０１１７】また、荷重配分推定部２６や車体加減速度
演算部３４では、各車輪の微小ゲインを各車輪の車輪速
微小振幅に置き換えることにより、荷重配分や車体加減
速度を演算することができる。荷重配分推定部２６や車
体加減速度演算部３４では、共振レベルの比が問題とな
るので、ブレーキ圧微小振幅によらず正確な演算値が得
られる。
演算部３４では、各車輪の微小ゲインを各車輪の車輪速
微小振幅に置き換えることにより、荷重配分や車体加減
速度を演算することができる。荷重配分推定部２６や車
体加減速度演算部３４では、共振レベルの比が問題とな
るので、ブレーキ圧微小振幅によらず正確な演算値が得
られる。
【０１１８】なお、第２の実施の形態において、図１３
の変換部への入力を、車輪速微小振幅とした場合、テー
ブル９２を少なくとも車輪速微小振幅と輪荷重とブレー
キ圧微小振幅との関係を示すテーブルにすると共に、他
のパラメータとして、少なくともブレーキ圧微小振幅を
用いることが必要になる。
の変換部への入力を、車輪速微小振幅とした場合、テー
ブル９２を少なくとも車輪速微小振幅と輪荷重とブレー
キ圧微小振幅との関係を示すテーブルにすると共に、他
のパラメータとして、少なくともブレーキ圧微小振幅を
用いることが必要になる。
【０１１９】また、第１の実施の形態では、微小励振手
段により共振周波数ω∞でブレーキ圧を微小励振した
が、路面等から共振周波数成分を含む自然外乱を受ける
ときは、微小励振手段を用いずに各車輪の共振レベルを
検出し、該共振レベルに基づいて荷重配分、車体加減速
度、及び路面状態を推定することができる。
段により共振周波数ω∞でブレーキ圧を微小励振した
が、路面等から共振周波数成分を含む自然外乱を受ける
ときは、微小励振手段を用いずに各車輪の共振レベルを
検出し、該共振レベルに基づいて荷重配分、車体加減速
度、及び路面状態を推定することができる。
【０１２０】この場合、図４の構成において、微小ゲイ
ン演算部２２ａ、２２ｂ、２２ｃ、２２ｄを、図８の車
輪速微小振幅検出部４０に置き換えることができる。す
なわち、共振レベルを車輪速微小振幅で表す。また、荷
重配分推定部２６や車体加減速度演算部３４では、各車
輪の微小ゲインを各車輪の車輪速微小振幅に置き換える
ことにより、荷重配分や車体加減速度を演算することが
できる。荷重配分推定部２６や車体加減速度演算部３４
では、共振レベルの比が問題となるので、外乱振幅によ
らず正確な演算値が得られる。
ン演算部２２ａ、２２ｂ、２２ｃ、２２ｄを、図８の車
輪速微小振幅検出部４０に置き換えることができる。す
なわち、共振レベルを車輪速微小振幅で表す。また、荷
重配分推定部２６や車体加減速度演算部３４では、各車
輪の微小ゲインを各車輪の車輪速微小振幅に置き換える
ことにより、荷重配分や車体加減速度を演算することが
できる。荷重配分推定部２６や車体加減速度演算部３４
では、共振レベルの比が問題となるので、外乱振幅によ
らず正確な演算値が得られる。
【０１２１】なお、第２の実施の形態で、図１３の変換
部への入力を、外乱による車輪速微小振幅とした場合、
外乱振幅の変化によって車輪速微小振幅が変化すること
を正規化するため、テーブル９２を少なくとも車輪速微
小振幅と輪荷重と外乱振幅との関係を示すテーブルにす
ると共に、他のパラメータとして、少なくとも共振周波
数成分での外乱振幅を用いることが必要になる。
部への入力を、外乱による車輪速微小振幅とした場合、
外乱振幅の変化によって車輪速微小振幅が変化すること
を正規化するため、テーブル９２を少なくとも車輪速微
小振幅と輪荷重と外乱振幅との関係を示すテーブルにす
ると共に、他のパラメータとして、少なくとも共振周波
数成分での外乱振幅を用いることが必要になる。
【０１２２】また、路面状態推定部３４において、高μ
路、中μ路、低μ路の３種類の路面状態を識別させた
が、本発明は、これに限定されるものでなく、２種類の
路面状態や、さらにきめ細かく４種類以上の路面状態を
識別させることも可能である。
路、中μ路、低μ路の３種類の路面状態を識別させた
が、本発明は、これに限定されるものでなく、２種類の
路面状態や、さらにきめ細かく４種類以上の路面状態を
識別させることも可能である。
【０１２３】
【発明の効果】以上説明したように、請求項１の発明に
よれば、各車輪の荷重配分を、実際に路面に接地してい
る車輪の共振レベルに基づいて直接的に求めるようにし
たので、荷重配分の推定を簡単かつ正確に行うことがで
きる、という効果が得られる。
よれば、各車輪の荷重配分を、実際に路面に接地してい
る車輪の共振レベルに基づいて直接的に求めるようにし
たので、荷重配分の推定を簡単かつ正確に行うことがで
きる、という効果が得られる。
【０１２４】また、請求項２の発明によれば、各車輪の
共振レベルに基づいて正確に加減速度を演算することが
できる、という効果が得られる。
共振レベルに基づいて正確に加減速度を演算することが
できる、という効果が得られる。
【０１２５】さらに、請求項３の発明によれば、各車輪
の共振レベルに基づく車体加減速度と車輪加減速度との
間の関係を用いて路面状態を推定するようにしたので、
路面状態の推定精度が向上する、という効果が得られ
る。
の共振レベルに基づく車体加減速度と車輪加減速度との
間の関係を用いて路面状態を推定するようにしたので、
路面状態の推定精度が向上する、という効果が得られ
る。
【図１】本発明の原理、特に、車体の加減速度運動によ
る前後の荷重移動の配分を説明するためのハーフカーモ
デルを示す図である。
る前後の荷重移動の配分を説明するためのハーフカーモ
デルを示す図である。
【図２】車輪と車体と路面とから構成される振動系の等
価モデルを示す図である。
価モデルを示す図である。
【図３】本発明の路面状態推定原理を説明するための、
車体の加減速度αと車輪の加減速度αw との間の関係を
路面状態毎に示した図である。
車体の加減速度αと車輪の加減速度αw との間の関係を
路面状態毎に示した図である。
【図４】本発明の第１の実施の形態に係る荷重配分推定
装置の構成を示すブロック図である。
装置の構成を示すブロック図である。
【図５】第１の実施の形態に係る荷重配分推定装置の荷
重配分推定部の詳細な構成を示すブロック図である。
重配分推定部の詳細な構成を示すブロック図である。
【図６】第１の実施の形態に係る荷重配分推定装置の車
体加減速度演算部の詳細な構成を示すブロック図であ
る。
体加減速度演算部の詳細な構成を示すブロック図であ
る。
【図７】第１の実施の形態に係る荷重配分推定装置の路
面状態推定部の詳細な構成を示すブロック図である。
面状態推定部の詳細な構成を示すブロック図である。
【図８】第１の実施の形態に係る荷重配分推定装置の微
小ゲイン演算部の詳細な構成を示すブロック図である。
小ゲイン演算部の詳細な構成を示すブロック図である。
【図９】第１の実施の形態に係る微小ゲイン演算部の車
輪速微小振幅演算部の詳細な構成を示すブロック図であ
る。
輪速微小振幅演算部の詳細な構成を示すブロック図であ
る。
【図１０】ブレーキ部のハードウェア構成を示すブロッ
ク図である。
ク図である。
【図１１】ブレーキ圧の微小励振と平均ブレーキ力の制
御を同時に行う場合の制御バルブへの指令を示す図であ
る。
御を同時に行う場合の制御バルブへの指令を示す図であ
る。
【図１２】実際の路面でブレーキ制動し、その後ブレー
キオフしたときの車体減速度比と車輪減速度比との間の
関係を路面状態毎に示したグラフである。
キオフしたときの車体減速度比と車輪減速度比との間の
関係を路面状態毎に示したグラフである。
【図１３】本発明の第２の実施の形態に係る変換部の構
成を示すブロック図である。
成を示すブロック図である。
【図１４】共振周波数の上昇とタイヤグリップ時の共振
周波数成分のゲイン減少を示す線図である。
周波数成分のゲイン減少を示す線図である。
２２ 微小ゲイン演算部 ２４ 遅延部 ２６ 荷重配分推定部 ２８ 車輪加減速度演算部 ３０ 車体加減速度演算部 ３２ 加減速度変化率演算部 ３４ 路面状態推定部 ４０ 車輪速微小振幅演算部 ４２ ブレーキ圧微小振幅演算部 ８０ 最小自乗法演算部
───────────────────────────────────────────────────── フロントページの続き (72)発明者 山口 裕之 愛知県愛知郡長久手町大字長湫字横道41 番地の１株式会社豊田中央研究所内 (72)発明者 菅井 賢 愛知県愛知郡長久手町大字長湫字横道41 番地の１株式会社豊田中央研究所内 (56)参考文献 特開 平５−229320（ＪＰ，Ａ) 特開 平５−133831（ＪＰ，Ａ) 特開 平７−159308（ＪＰ，Ａ) 特開 平６−213801（ＪＰ，Ａ) 特開 平６−258196（ＪＰ，Ａ) 特開 平５−330322（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) G01M 17/007 G01L 1/10
Claims (3)
- 【請求項１】 車輪速度が微小振動しているときの各車
輪の荷重配分を推定する荷重配分推定装置であって、 車体と車輪と路面とから構成される振動系の共振周波数
での車輪速度の微小振幅に基づく共振レベルを、各車輪
毎に検出する検出手段と、 車輪にかかる荷重と該車輪の共振レベルとが正の相関関
係にあることに基づいて、前記検出手段により検出され
た各車輪の共振レベルにより、各車輪の荷重配分を推定
する荷重配分推定手段と、 を有することを特徴とする荷重配分推定装置。 - 【請求項２】 車輪の車輪速度が微小振動しているとき
の車体の加減速度を演算する車体加減速度演算装置であ
って、 車体と車輪と路面とから構成される振動系の共振周波数
での車輪速度の微小振幅に基づく共振レベルを、前輪と
後輪について各々検出する第１の検出手段、及び右輪と
左輪について各々検出する第２の検出手段の少なくとも
いずれかの検出手段と、 前記第１の検出手段により検出された前輪の共振レベル
と後輪の共振レベルとに基づいて車体の前後方向の加減
速度を演算する第１の演算手段、及び前記第２の検出手
段により検出された右輪の共振レベルと左輪の共振レベ
ルとに基づいて車体の横方向の加減速度を演算する第２
の演算手段の少なくともいずれかの演算手段と、 を有することを特徴とする車体加減速度演算装置。 - 【請求項３】 車輪の車輪速度が微小振動している車体
が走行する路面の状態を推定する路面状態推定装置であ
って、 車体と車輪と路面とから構成される振動系の共振周波数
での車輪速度の微小振幅に基づく共振レベルを、前輪と
後輪について各々検出する検出手段と、 前記検出手段により検出された前輪の共振レベルと後輪
の共振レベルとに基づいて車体の加減速度を演算する演
算手段と、 車輪の加減速度を検出する車輪加減速度検出手段と、 前記演算手段により演算された車体の加減速度及び前記
車輪加減速度検出手段により検出された車輪の加減速度
の間の関係と、予めデータとして備えられている車体の
加減速度及び車輪の加減速度の間の路面状態毎の関係
と、の比較に基づいて路面状態を推定する路面推定手段
と、 を有することを特徴とする路面状態推定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP18160897A JP3317196B2 (ja) | 1997-07-07 | 1997-07-07 | 荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP18160897A JP3317196B2 (ja) | 1997-07-07 | 1997-07-07 | 荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1123425A JPH1123425A (ja) | 1999-01-29 |
| JP3317196B2 true JP3317196B2 (ja) | 2002-08-26 |
Family
ID=16103790
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP18160897A Expired - Fee Related JP3317196B2 (ja) | 1997-07-07 | 1997-07-07 | 荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3317196B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4296705B2 (ja) | 2000-10-25 | 2009-07-15 | 株式会社デンソー | 車両用ブレーキ装置 |
| JP2006084424A (ja) | 2004-09-17 | 2006-03-30 | Sumitomo Rubber Ind Ltd | 荷重配分判定方法 |
| JP2019113373A (ja) * | 2017-12-22 | 2019-07-11 | 住友ゴム工業株式会社 | 輪荷重推定装置 |
| JP7463839B2 (ja) * | 2020-05-20 | 2024-04-09 | 株式会社デンソー | 判定装置 |
| CN113954850B (zh) * | 2021-11-26 | 2023-04-07 | 上海交通大学 | 基于边缘计算和软测量的车轮轮心载荷估计方法和*** |
-
1997
- 1997-07-07 JP JP18160897A patent/JP3317196B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH1123425A (ja) | 1999-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6324461B1 (en) | Road surface condition estimating apparatus and variation reduction processing apparatus | |
| EP0891904B1 (en) | Wheel condition estimating apparatus | |
| EP2853438B1 (en) | Traction control device, and traction control method | |
| KR100760486B1 (ko) | 타이어 그립도 추정 장치 및 방법, 그리고 주행 상태 제어방법 | |
| US5402345A (en) | Method for estimating speed of a braked vehicle | |
| EP0699568A2 (en) | Anti-lock brake controlling apparatus | |
| JP2002120709A (ja) | 路面状態推定装置 | |
| JP5839104B2 (ja) | 車両の旋回特性推定装置 | |
| Boisvert et al. | Estimators of wheel slip for electric vehicles using torque and encoder measurements | |
| JP3344951B2 (ja) | 物理量推定装置及びａｂｓ制御装置 | |
| US7059687B2 (en) | Braking force distribution control device | |
| JP3317196B2 (ja) | 荷重配分推定装置、車体加減速度演算装置及び路面状態推定装置 | |
| JP3619388B2 (ja) | 車両の重心高さの推定演算装置 | |
| JP3331939B2 (ja) | 車両状態量推定装置 | |
| JPH1178843A (ja) | 車輪状態推定装置 | |
| JP3435625B2 (ja) | 路面状態演算装置 | |
| JP3454090B2 (ja) | アンチロックブレーキ制御装置、トルク勾配推定装置及び制動トルク勾配推定装置 | |
| JP3424535B2 (ja) | 路面状態推定装置 | |
| JP3454086B2 (ja) | 車輪挙動量サーボ制御装置及び限界判定装置 | |
| JP3640126B2 (ja) | スリップ率サーボ制御装置 | |
| JP2001354020A (ja) | サスペンション制御装置 | |
| JP2003237560A (ja) | 最大路面摩擦係数推定装置 | |
| JP3702649B2 (ja) | 制動力推定装置及びブレーキ圧推定装置 | |
| JP4674959B2 (ja) | 路面状態推定装置 | |
| JPH0952572A (ja) | アンチロックブレーキ制御装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |