JP2002173044A - 路面摩擦状態推定装置 - Google Patents
路面摩擦状態推定装置Info
- Publication number
- JP2002173044A JP2002173044A JP2000370709A JP2000370709A JP2002173044A JP 2002173044 A JP2002173044 A JP 2002173044A JP 2000370709 A JP2000370709 A JP 2000370709A JP 2000370709 A JP2000370709 A JP 2000370709A JP 2002173044 A JP2002173044 A JP 2002173044A
- Authority
- JP
- Japan
- Prior art keywords
- road surface
- vehicle
- surface friction
- friction state
- steering amount
- 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.)
- Pending
Links
Landscapes
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Regulating Braking Force (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
(57)【要約】
【課題】 走行中の路面の摩擦係数をリアルタイムで推
定する。 【解決手段】 操向される車輪に作用するSAT(セル
フアライニングトルク)を反力算出部70にて算出す
る。また、走行される車輪の横すべり角を操向量算出部
80にて算出する。SATを横すべり角の一次関数とし
て表すために、算出されたSATと横すべり角からオン
ライン最小二乗法を用いて、その関数の勾配、切片を、
関数決定部90により求める。この勾配と切片より、路
面の摩擦係数を摩擦状態算出部100にて算出する。
定する。 【解決手段】 操向される車輪に作用するSAT(セル
フアライニングトルク)を反力算出部70にて算出す
る。また、走行される車輪の横すべり角を操向量算出部
80にて算出する。SATを横すべり角の一次関数とし
て表すために、算出されたSATと横すべり角からオン
ライン最小二乗法を用いて、その関数の勾配、切片を、
関数決定部90により求める。この勾配と切片より、路
面の摩擦係数を摩擦状態算出部100にて算出する。
Description
【0001】
【発明の属する技術分野】本発明は、操向される車輪を
有する車両が、走行中の路面の摩擦状態を推定する装置
に関する。
有する車両が、走行中の路面の摩擦状態を推定する装置
に関する。
【0002】
【従来の技術】車輪を有し、路面上を走行する車両の運
動は、最終的には、車輪と路面の摩擦力により実現され
ている。言い換えれば、車輪と路面間に発生する摩擦力
を超えた車両の挙動は実現できない。すなわち、制動距
離、旋回半径と速度などの限界値は、車輪と路面間の摩
擦力により決定される。したがって、車輪路面間の摩擦
力を決定する大きな因子が路面の摩擦係数である。この
摩擦係数は、路面ごとに異なることに加え、路面の状況
によっても異なる。前者は、アスファルト路面か、コン
クリート路面か、砂利道か、などであり、後者は、乾い
ているか、濡れているか、雪が積もっているか、などで
ある。このような様々な状況下で、路面の摩擦係数は大
きく変化するが、車両の挙動を安定させるためには、現
在走行中の路面の状態、または摩擦係数が取得できると
好ましい。路面の状態を、「滑りやすい」「滑りにく
い」などと運転者に報知することもできる。また、アン
チロックブレーキ、トラクション制御、車両安定性制御
などのシステムにおいて、現在の路面状況を用いた制御
を行うことも可能となる。
動は、最終的には、車輪と路面の摩擦力により実現され
ている。言い換えれば、車輪と路面間に発生する摩擦力
を超えた車両の挙動は実現できない。すなわち、制動距
離、旋回半径と速度などの限界値は、車輪と路面間の摩
擦力により決定される。したがって、車輪路面間の摩擦
力を決定する大きな因子が路面の摩擦係数である。この
摩擦係数は、路面ごとに異なることに加え、路面の状況
によっても異なる。前者は、アスファルト路面か、コン
クリート路面か、砂利道か、などであり、後者は、乾い
ているか、濡れているか、雪が積もっているか、などで
ある。このような様々な状況下で、路面の摩擦係数は大
きく変化するが、車両の挙動を安定させるためには、現
在走行中の路面の状態、または摩擦係数が取得できると
好ましい。路面の状態を、「滑りやすい」「滑りにく
い」などと運転者に報知することもできる。また、アン
チロックブレーキ、トラクション制御、車両安定性制御
などのシステムにおいて、現在の路面状況を用いた制御
を行うことも可能となる。
【0003】現在走行中の路面の摩擦係数を測定する装
置が特開平11−287749号公報に記載されてい
る。この装置は、ステアリングホイールを切り増しした
ときに、この回転角(操舵角)の変化量に対する操舵ト
ルクの変化量を求め、この値から路面の摩擦係数を算出
している。
置が特開平11−287749号公報に記載されてい
る。この装置は、ステアリングホイールを切り増しした
ときに、この回転角(操舵角)の変化量に対する操舵ト
ルクの変化量を求め、この値から路面の摩擦係数を算出
している。
【0004】
【発明が解決しようとする課題】前述の公報の装置は、
ノイズの影響を受けやすいという欠点がある。これは、
摩擦係数を算出するための物理量として、操舵角の変化
量と操舵トルクの変化量という、ノイズを増幅させる傾
向のある物理量を用いているためである。ある物理量の
変化量を求めるということは、ある物理量を微分すると
いうことであり、微分演算は、ノイズを増幅させてしま
う。よって、前記公報の装置により算出された摩擦係数
はノイズを多く含むものとなるという問題があった。
ノイズの影響を受けやすいという欠点がある。これは、
摩擦係数を算出するための物理量として、操舵角の変化
量と操舵トルクの変化量という、ノイズを増幅させる傾
向のある物理量を用いているためである。ある物理量の
変化量を求めるということは、ある物理量を微分すると
いうことであり、微分演算は、ノイズを増幅させてしま
う。よって、前記公報の装置により算出された摩擦係数
はノイズを多く含むものとなるという問題があった。
【0005】また、前記公報の装置は、切り増し時にの
み摩擦係数の算出を行うものである。したがって、これ
以上は切り増しを行えないという限界時において摩擦係
数を算出できない。最も欲しい摩擦係数は、限界時にお
ける摩擦係数であり、この装置では、この値を求めるこ
とができない。
み摩擦係数の算出を行うものである。したがって、これ
以上は切り増しを行えないという限界時において摩擦係
数を算出できない。最も欲しい摩擦係数は、限界時にお
ける摩擦係数であり、この装置では、この値を求めるこ
とができない。
【0006】さらにまた、車輪路面間の摩擦係数は、操
舵角に対してヒステリシスを有し、切り増し時と切り戻
し時とで異なる特性を有する。前記公報の装置では、切
り増し時にのみ摩擦係数の算出を行うこととし、ヒステ
リシスの影響を排除しているが、切り戻し時や一定舵角
を保持しているときなどに、摩擦係数の算出を行えない
という問題があった。
舵角に対してヒステリシスを有し、切り増し時と切り戻
し時とで異なる特性を有する。前記公報の装置では、切
り増し時にのみ摩擦係数の算出を行うこととし、ヒステ
リシスの影響を排除しているが、切り戻し時や一定舵角
を保持しているときなどに、摩擦係数の算出を行えない
という問題があった。
【0007】本発明は、前述の問題点を解決するために
なされたものであり、ノイズの影響を受けにくく、精度
が高く、また推定頻度の高い路面摩擦係数を推定する装
置を提供することを目的とする。
なされたものであり、ノイズの影響を受けにくく、精度
が高く、また推定頻度の高い路面摩擦係数を推定する装
置を提供することを目的とする。
【0008】
【課題を解決するための手段】前述の課題を解決するた
めに、本発明にかかる路面摩擦状態推定装置は、車両の
操向される車輪が路面より受ける反力を、車輪の操向量
の一次関数として、この関数の勾配と切片を求め、これ
らから路面の摩擦状態を算出するものである。勾配と切
片を分離して求めることにより、路面状態によって変化
する、反力と操向量のヒステリシスの影響を受けずに、
勾配を求めることができる。
めに、本発明にかかる路面摩擦状態推定装置は、車両の
操向される車輪が路面より受ける反力を、車輪の操向量
の一次関数として、この関数の勾配と切片を求め、これ
らから路面の摩擦状態を算出するものである。勾配と切
片を分離して求めることにより、路面状態によって変化
する、反力と操向量のヒステリシスの影響を受けずに、
勾配を求めることができる。
【0009】また、勾配と切片の算出は、取得した反力
と操向量に基づきオンライン最小二乗法を適用して行わ
れる。関数の決定に当たって取得した物理量の微分量を
排除したことにより、ノイズの影響を低減させることが
できる。
と操向量に基づきオンライン最小二乗法を適用して行わ
れる。関数の決定に当たって取得した物理量の微分量を
排除したことにより、ノイズの影響を低減させることが
できる。
【0010】また、操向量は、車両の運動に関する特性
を、モデル化してあらかじめ記憶しておき、この特性を
適用して算出するようにできる。すなわち、車両の運動
状態を表す所定の物理量を検出して、この物理量を前記
モデル化された特性に適用して車両の運動状態を解析す
る。操向量はこの解析結果に含まれる。
を、モデル化してあらかじめ記憶しておき、この特性を
適用して算出するようにできる。すなわち、車両の運動
状態を表す所定の物理量を検出して、この物理量を前記
モデル化された特性に適用して車両の運動状態を解析す
る。操向量はこの解析結果に含まれる。
【0011】前述の操向量は車輪の横すべり角とすると
ができる。
ができる。
【0012】
【発明の実施の形態】以下、本発明の実施の形態(以下
実施形態という)を、図面に従って説明する。図1は、
本実施形態にかかる操舵装置20の概略構成図である。
操舵装置20は、ラックアンドピニオン機構を有するパ
ワーステアリング装置である。運転者が車両の旋回など
を行うために操作するステアリングホイール22は、ス
テアリングシャフト24を介して、前記ラックアンドピ
ニオン機構のピニオン26に結合している。ピニオン2
6は、ラックロッド28上に設けられたラックとかみ合
い、ここでステアリングホイール22の回転運動がラッ
クロッド28の直線運動に変換される。このラックロッ
ド28の動きがタイロッドを介してナックルに伝達され
操向車輪30,32の向き、すなわち車輪の操向角が変
更される。
実施形態という)を、図面に従って説明する。図1は、
本実施形態にかかる操舵装置20の概略構成図である。
操舵装置20は、ラックアンドピニオン機構を有するパ
ワーステアリング装置である。運転者が車両の旋回など
を行うために操作するステアリングホイール22は、ス
テアリングシャフト24を介して、前記ラックアンドピ
ニオン機構のピニオン26に結合している。ピニオン2
6は、ラックロッド28上に設けられたラックとかみ合
い、ここでステアリングホイール22の回転運動がラッ
クロッド28の直線運動に変換される。このラックロッ
ド28の動きがタイロッドを介してナックルに伝達され
操向車輪30,32の向き、すなわち車輪の操向角が変
更される。
【0013】ステアリングシャフト24上には、このシ
ャフトに加わるトルクを検出するトルクセンサ40と、
シャフトの回転角度である操舵角を検出する操舵角セン
サ42が設けられている。これらのセンサの出力は、制
御部50に送られる。
ャフトに加わるトルクを検出するトルクセンサ40と、
シャフトの回転角度である操舵角を検出する操舵角セン
サ42が設けられている。これらのセンサの出力は、制
御部50に送られる。
【0014】制御部50は、CPU(中央処理装置)5
2を中心とするマイクロプロセッサとして構成されてお
り、処理プログラムを記憶したROM(読み出し専用メ
モリ)54と、一時的にデータを記憶するRAM(ラン
ダムアクセスメモリ)56と、入出力ポート(図示せ
ず)とを備える。制御部50には、トルクセンサ40か
らのトルク信号や、操舵角センサ42からの操舵角信
号、さらには、車体に設けられた車速センサ58からの
車速信号が入力される。
2を中心とするマイクロプロセッサとして構成されてお
り、処理プログラムを記憶したROM(読み出し専用メ
モリ)54と、一時的にデータを記憶するRAM(ラン
ダムアクセスメモリ)56と、入出力ポート(図示せ
ず)とを備える。制御部50には、トルクセンサ40か
らのトルク信号や、操舵角センサ42からの操舵角信
号、さらには、車体に設けられた車速センサ58からの
車速信号が入力される。
【0015】図2は、制御部50が、路面の摩擦状態、
例えば摩擦係数の算出を行う装置として機能する場合に
かかる制御ブロック図である。反力算出部70は、トル
クセンサ40の出力、すなわち操舵トルクTpから、ラ
ックアンドピニオン機構のギア比、車輪操向機構のジオ
メトリなどを考慮して、セルフアライニングトルク(以
下、SATと記す)を算出する。
例えば摩擦係数の算出を行う装置として機能する場合に
かかる制御ブロック図である。反力算出部70は、トル
クセンサ40の出力、すなわち操舵トルクTpから、ラ
ックアンドピニオン機構のギア比、車輪操向機構のジオ
メトリなどを考慮して、セルフアライニングトルク(以
下、SATと記す)を算出する。
【0016】一方、操向量算出部80は、車両の運動方
程式から、操向量として横すべり角αfを算出する。車
両の運動特性は、運動方程式から導かれた状態方程式と
して特性記憶部82に記憶されている。この記憶された
特性に従って、解析部84にて、車両の運動状態を示す
物理量、例えば車両の車速u、操舵角θsに基づき現在
の運動状態が解析される。この状態に、横すべり角αf
が含まれる。なお、横すべり角αfは、車輪を上から見
たときの車輪中心面と車輪の進行方向のなす角として定
義される。
程式から、操向量として横すべり角αfを算出する。車
両の運動特性は、運動方程式から導かれた状態方程式と
して特性記憶部82に記憶されている。この記憶された
特性に従って、解析部84にて、車両の運動状態を示す
物理量、例えば車両の車速u、操舵角θsに基づき現在
の運動状態が解析される。この状態に、横すべり角αf
が含まれる。なお、横すべり角αfは、車輪を上から見
たときの車輪中心面と車輪の進行方向のなす角として定
義される。
【0017】関数決定部90において、求められたSA
Tと横すべり角αfの関係を決定する。すなわち、SA
Tを横すべり角の一次関数とみなし、この関数の勾配と
切片を求める。この勾配と切片に基づき摩擦状態算出部
100にて、路面の摩擦状態として摩擦係数を算出す
る。
Tと横すべり角αfの関係を決定する。すなわち、SA
Tを横すべり角の一次関数とみなし、この関数の勾配と
切片を求める。この勾配と切片に基づき摩擦状態算出部
100にて、路面の摩擦状態として摩擦係数を算出す
る。
【0018】次に、個々の算出部、決定部の演算につい
て詳細に説明する。前述したように、反力算出部70で
は、トルクセンサ40からの出力信号と、操舵機構の幾
何学的条件からSATを算出する。なお、SATは、タ
イロッドの軸力を検出するセンサを設け、このセンサ出
力に基づき算出することも可能である。
て詳細に説明する。前述したように、反力算出部70で
は、トルクセンサ40からの出力信号と、操舵機構の幾
何学的条件からSATを算出する。なお、SATは、タ
イロッドの軸力を検出するセンサを設け、このセンサ出
力に基づき算出することも可能である。
【0019】操向量算出部80は、特性記憶部82に記
憶された車両の運動状態モデルに基づき、操向量として
横すべり角αfを算出する。運動状態のモデルは、次式
(1),(2)の状態方程式として記憶されている。
憶された車両の運動状態モデルに基づき、操向量として
横すべり角αfを算出する。運動状態のモデルは、次式
(1),(2)の状態方程式として記憶されている。
【0020】
【数1】
【0021】解析部84において、式(1),(2)
は、離散化された状態で取り扱われる。また、式(1)
において、入力は車速uと操舵角θsであり、横速度ν
およびヨーレートrが、出力される。式(2)におい
て、横速度ν、ヨーレートrに基づき前輪横すべり角α
fが算出される。
は、離散化された状態で取り扱われる。また、式(1)
において、入力は車速uと操舵角θsであり、横速度ν
およびヨーレートrが、出力される。式(2)におい
て、横速度ν、ヨーレートrに基づき前輪横すべり角α
fが算出される。
【0022】関数決定部90では、前述のように求めら
れたSATと横すべり角αfから、次式(3)〜(7)
で記述されるオンライン最小二乗法によって動作点A周
りの前輪横すべり角に対するSATの関数関係(勾配
k、切片T)を演算する。
れたSATと横すべり角αfから、次式(3)〜(7)
で記述されるオンライン最小二乗法によって動作点A周
りの前輪横すべり角に対するSATの関数関係(勾配
k、切片T)を演算する。
【0023】
【数2】
【0024】θの第1要素は、関数の勾配、第2要素は
切片である。SATは、図3に示す破線のように、横す
べり角αfに対して、単調に増加し、その後単調に減少
する特性を示すが、動作点Aにおける接線として近似す
る。この接線の勾配kと切片Tが前述のオンライン最小
二乗法により求められる。
切片である。SATは、図3に示す破線のように、横す
べり角αfに対して、単調に増加し、その後単調に減少
する特性を示すが、動作点Aにおける接線として近似す
る。この接線の勾配kと切片Tが前述のオンライン最小
二乗法により求められる。
【0025】路面摩擦状態算出部100においては、前
記の勾配k、切片Tより次式(8)に基づき摩擦係数μ
を算出する。
記の勾配k、切片Tより次式(8)に基づき摩擦係数μ
を算出する。
【0026】
【数3】
【0027】式(8)により、摩擦係数μを算出するの
に代えて、例えば、図4に示すようなマップを用いて路
面の摩擦状態を得ることも可能である。図4では、勾配
kと切片Tが共に小さな値となる場合は滑りやすい路面
(低μ路)、勾配kと切片Tのいずれかが大きな値とな
る場合は滑りにくい路面(高μ路)、中間を中μ路と3
段階に分けている。
に代えて、例えば、図4に示すようなマップを用いて路
面の摩擦状態を得ることも可能である。図4では、勾配
kと切片Tが共に小さな値となる場合は滑りやすい路面
(低μ路)、勾配kと切片Tのいずれかが大きな値とな
る場合は滑りにくい路面(高μ路)、中間を中μ路と3
段階に分けている。
【0028】以上のように求められた路面の摩擦状態
は、これに基づき運転者に警告を与えたり、車両挙動制
御の因子として用いたりすることができる。
は、これに基づき運転者に警告を与えたり、車両挙動制
御の因子として用いたりすることができる。
【0029】次に、本実施形態の算出方法の特徴につい
て説明する。前述したように、本実施形態は、オンライ
ンで演算を行っているので、リアルタイムで路面摩擦係
数を取得することができる。このため、路面状態が急変
した場合においても、すばやく摩擦係数を算出すること
ができる。
て説明する。前述したように、本実施形態は、オンライ
ンで演算を行っているので、リアルタイムで路面摩擦係
数を取得することができる。このため、路面状態が急変
した場合においても、すばやく摩擦係数を算出すること
ができる。
【0030】また、路面車輪間の摩擦係数は発生してい
る摩擦力により変化し、よって車両の限界状態の摩擦係
数は、最大摩擦力発生時、またはこれに近い状態のとき
算出するのが望ましい。従来技術として挙げた前述の公
報は、切り増し時のみ摩擦係数の算出を行っているの
で、最大摩擦力発生時およびその近傍に算出が行われる
とは限らず、限界の摩擦係数を求めることが必ずしもで
きない。一方、本実施形態によれば、横すべり角が極大
値に達したところで、摩擦係数を算出することにより、
より限界に近い摩擦係数を得ることができる。
る摩擦力により変化し、よって車両の限界状態の摩擦係
数は、最大摩擦力発生時、またはこれに近い状態のとき
算出するのが望ましい。従来技術として挙げた前述の公
報は、切り増し時のみ摩擦係数の算出を行っているの
で、最大摩擦力発生時およびその近傍に算出が行われる
とは限らず、限界の摩擦係数を求めることが必ずしもで
きない。一方、本実施形態によれば、横すべり角が極大
値に達したところで、摩擦係数を算出することにより、
より限界に近い摩擦係数を得ることができる。
【0031】また、本実施形態の算出手法は、切り増
し、切り戻しに関わらず、摩擦係数の算出を行える。横
すべり角とSATの間にはヒステリシスが存在し、ステ
アリングホイールを切り増すときと、切り戻す時では、
異なる軌跡を描く。図5には、SATの変化を1次近似
したヒステリシス特性が示されている。切り増し時と切
り戻し時において、勾配は変わらず、切片のみが変化し
ている。前述のように、本実施形態においては、式
(3)〜(7)により、勾配と切片が分離されるので、
切り増し時、切り戻し時双方にて精度よく摩擦係数の算
出を行うことができる。
し、切り戻しに関わらず、摩擦係数の算出を行える。横
すべり角とSATの間にはヒステリシスが存在し、ステ
アリングホイールを切り増すときと、切り戻す時では、
異なる軌跡を描く。図5には、SATの変化を1次近似
したヒステリシス特性が示されている。切り増し時と切
り戻し時において、勾配は変わらず、切片のみが変化し
ている。前述のように、本実施形態においては、式
(3)〜(7)により、勾配と切片が分離されるので、
切り増し時、切り戻し時双方にて精度よく摩擦係数の算
出を行うことができる。
【0032】さらに、前述のヒステリシス特性は、路面
の摩擦係数によって、変化する。これは、このヒステリ
シスの発生が、主に車輪のタイヤ、特にそのトレッド面
のねじれによるものであり、摩擦係数の違いによってね
じれの程度が変化するために生じるものである。よっ
て、摩擦係数が高いほどヒステリシスが大きくなる、つ
まり、切り増し、切り戻しの切片の差が大きくなる。こ
のように、路面の摩擦係数によってヒステリシス特性が
変化するため、ヒステリシス特性を一定のものとして取
り扱うことはできず、本実施形態の手法が有効である。
図6には、路面の摩擦係数μの高低によるヒステリシス
特性の違いが示されている。低摩擦係数の路面において
は、ヒステリシスがほとんどなくなっていることが分か
る。
の摩擦係数によって、変化する。これは、このヒステリ
シスの発生が、主に車輪のタイヤ、特にそのトレッド面
のねじれによるものであり、摩擦係数の違いによってね
じれの程度が変化するために生じるものである。よっ
て、摩擦係数が高いほどヒステリシスが大きくなる、つ
まり、切り増し、切り戻しの切片の差が大きくなる。こ
のように、路面の摩擦係数によってヒステリシス特性が
変化するため、ヒステリシス特性を一定のものとして取
り扱うことはできず、本実施形態の手法が有効である。
図6には、路面の摩擦係数μの高低によるヒステリシス
特性の違いが示されている。低摩擦係数の路面において
は、ヒステリシスがほとんどなくなっていることが分か
る。
【0033】図7は、高μ路において、速度40km/
hで走行中に、±40degの操舵を連続して行ったと
きの操舵角と推定された前輪横すべり角および前輪横す
べり角の推定を行うモデルの検証のためこのモデルよっ
て導出されるヨーレートと、ヨーレートの実測値の比較
を示したものである。ヨーレートのモデル出力と実測値
はよく一致しており、このことから前輪横すべり角が正
確に推定できていることがわかる。
hで走行中に、±40degの操舵を連続して行ったと
きの操舵角と推定された前輪横すべり角および前輪横す
べり角の推定を行うモデルの検証のためこのモデルよっ
て導出されるヨーレートと、ヨーレートの実測値の比較
を示したものである。ヨーレートのモデル出力と実測値
はよく一致しており、このことから前輪横すべり角が正
確に推定できていることがわかる。
【0034】図8は、推定された前輪横すべり角と、S
AT実測値からオンライン最小二乗法を適用して推定さ
れたSATの勾配kと切片T、更にこれらから推定され
た摩擦係数μの推定値を示している。摩擦係数μ推定値
は振動的ではあるが、後述する低μ路に比して高い値を
示しており、確かに高μ路を走行していることが分か
る。また、前輪横すべり角との対応から、前輪横すべり
角が大きな限界に近い領域で、摩擦係数μの推定値の変
動幅は小さくなっており、推定の精度が向上しているこ
とが分かる。さらに、切片Tに着目すると、操舵に応じ
て矩形は的に変化しており、ヒステリシスに伴う特性変
化が切片Tの値に現れていることが分かる。
AT実測値からオンライン最小二乗法を適用して推定さ
れたSATの勾配kと切片T、更にこれらから推定され
た摩擦係数μの推定値を示している。摩擦係数μ推定値
は振動的ではあるが、後述する低μ路に比して高い値を
示しており、確かに高μ路を走行していることが分か
る。また、前輪横すべり角との対応から、前輪横すべり
角が大きな限界に近い領域で、摩擦係数μの推定値の変
動幅は小さくなっており、推定の精度が向上しているこ
とが分かる。さらに、切片Tに着目すると、操舵に応じ
て矩形は的に変化しており、ヒステリシスに伴う特性変
化が切片Tの値に現れていることが分かる。
【0035】図9は、低μ路において、図7の実験と同
様に速度40km/hで走行中に±15degの操舵を
連続して行ったときの操舵角と推定された前輪横すべり
角および前輪横すべり角の推定を行うモデルの検証のた
めこのモデルよって導出されるヨーレートと、ヨーレー
トの実測値の比較を示したものである。この実験におい
ても、ヨーレートのモデル出力と実測値はよく一致して
おり、前輪横すべり角が正確に推定できていることがわ
かる。また、低μ路においても高μ路と同じパラメータ
を用いたモデルが実測値と一致しているということは、
車輪路面間の応力特性には、路面摩擦係数μの影響が現
れていないといえる。すなわち、応力特性を表すパラメ
ータであるコーナリングスティフネスは、低μ路走行時
においても変化していないことが分かる。したがって、
このような定常走行に近い走行条件では、コーナリング
スティフネス変化はほとんど見られず、コーナリングス
ティフネス変化が車両バネ上特性に及ぼす影響から低μ
路走行を判別することは不可能であることが分かる。こ
れに対し、SATには、このような定常走行に近い走行
条件においても路面μに依存した変化が現れることか
ら、SATの勾配と切片から路面摩擦状態を推定する本
実施形態では、低μ路走行を判別することが可能であ
る。
様に速度40km/hで走行中に±15degの操舵を
連続して行ったときの操舵角と推定された前輪横すべり
角および前輪横すべり角の推定を行うモデルの検証のた
めこのモデルよって導出されるヨーレートと、ヨーレー
トの実測値の比較を示したものである。この実験におい
ても、ヨーレートのモデル出力と実測値はよく一致して
おり、前輪横すべり角が正確に推定できていることがわ
かる。また、低μ路においても高μ路と同じパラメータ
を用いたモデルが実測値と一致しているということは、
車輪路面間の応力特性には、路面摩擦係数μの影響が現
れていないといえる。すなわち、応力特性を表すパラメ
ータであるコーナリングスティフネスは、低μ路走行時
においても変化していないことが分かる。したがって、
このような定常走行に近い走行条件では、コーナリング
スティフネス変化はほとんど見られず、コーナリングス
ティフネス変化が車両バネ上特性に及ぼす影響から低μ
路走行を判別することは不可能であることが分かる。こ
れに対し、SATには、このような定常走行に近い走行
条件においても路面μに依存した変化が現れることか
ら、SATの勾配と切片から路面摩擦状態を推定する本
実施形態では、低μ路走行を判別することが可能であ
る。
【0036】図10は、推定された前輪横すべり角と、
SAT実測値からオンライン最小二乗法を適用して推定
されたSATの勾配kと切片T、更にこれらから推定さ
れた摩擦係数μの推定値を示している。摩擦係数μ推定
値は振動的ではあるが、前述の項μ路に比して低い値を
示しており、確かに低μ路を走行していることが分か
る。また、図8の場合と同様、前輪横すべり角が大きな
限界に近い領域で、摩擦係数μの推定値の変動幅は小さ
くなっており、推定の精度が向上していることが分か
る。さらに、切片Tに着目すると、図8と同様、操舵に
応じて矩形波的に変化しているが、図8の切片Tの波形
に比較すると振幅が小さいことが分かる。これは路面摩
擦係数μに応じてヒステリシス特性が異なったことに起
因しているが、本手法では、このようなヒステリシス特
性の違いにかかわらず、SATの勾配kおよび摩擦係数
μの推定を適切に行える。
SAT実測値からオンライン最小二乗法を適用して推定
されたSATの勾配kと切片T、更にこれらから推定さ
れた摩擦係数μの推定値を示している。摩擦係数μ推定
値は振動的ではあるが、前述の項μ路に比して低い値を
示しており、確かに低μ路を走行していることが分か
る。また、図8の場合と同様、前輪横すべり角が大きな
限界に近い領域で、摩擦係数μの推定値の変動幅は小さ
くなっており、推定の精度が向上していることが分か
る。さらに、切片Tに着目すると、図8と同様、操舵に
応じて矩形波的に変化しているが、図8の切片Tの波形
に比較すると振幅が小さいことが分かる。これは路面摩
擦係数μに応じてヒステリシス特性が異なったことに起
因しているが、本手法では、このようなヒステリシス特
性の違いにかかわらず、SATの勾配kおよび摩擦係数
μの推定を適切に行える。
【図1】 本実施形態の装置の概略構成を示す図であ
る。
る。
【図2】 本実施形態の演算処理にかかるブロック図で
ある。
ある。
【図3】 SATを横すべり角の一次関数として例示し
た図である。
た図である。
【図4】 路面摩擦係数を算出するためのマップの例を
示す図である。
示す図である。
【図5】 SATのヒステリシス特性を例示した図であ
る。
る。
【図6】 摩擦係数の異なる路面におけるSATのヒス
テリシス特性の相違を示す図である。
テリシス特性の相違を示す図である。
【図7】 摩擦係数の高い路面において、左右に操舵を
連続し行った場合の、操舵角、横すべり角、実測値とモ
デル解析とによるヨーレートを示す図である。
連続し行った場合の、操舵角、横すべり角、実測値とモ
デル解析とによるヨーレートを示す図である。
【図8】 図7の条件で、SATを示す一次関数の勾配
と切片、およびこれらから導かれた路面摩擦係数を示す
図である。
と切片、およびこれらから導かれた路面摩擦係数を示す
図である。
【図9】 摩擦係数の低い路面において、左右に操舵を
連続し行った場合の、操舵角、横すべり角、実測値とモ
デル解析とによるヨーレートを示す図である。
連続し行った場合の、操舵角、横すべり角、実測値とモ
デル解析とによるヨーレートを示す図である。
【図10】 図9の条件で、SATを示す一次関数の勾
配と切片、およびこれらから導かれた路面摩擦係数を示
す図である。
配と切片、およびこれらから導かれた路面摩擦係数を示
す図である。
20 操舵装置、22 ステアリングホイール、24
ステアリングシャフト、26 ピニオン、28 ラック
ロッド、30,32 操向車輪、40 トルクセンサ、
42 操舵角センサ、50 制御部、58 車速セン
サ。
ステアリングシャフト、26 ピニオン、28 ラック
ロッド、30,32 操向車輪、40 トルクセンサ、
42 操舵角センサ、50 制御部、58 車速セン
サ。
フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) B62D 113:00 B62D 113:00 119:00 119:00 Fターム(参考) 2F051 AA01 BA00 3D032 CC05 DA03 DA15 DA23 DD02 EA01 EB04 EB16 EB21 GG01 3D046 BB23 HH08 HH21 HH22 HH36 HH46
Claims (4)
- 【請求項1】 操向される車輪を有する車両が走行中の
路面の摩擦状態を推定する路面摩擦状態推定装置であっ
て、 前記操向される車輪が路面より受ける反力を取得する反
力取得手段と、 前記操向される車輪の操向量を取得する操向量取得手段
と、 前記反力を前記操向量の一次関数として、前記取得した
反力と操向量に基づき、前記関数の勾配と切片を算出す
る関数決定手段と、 前記関数の勾配と切片から、路面の摩擦状態を算出する
路面摩擦状態算出手段と、を有する路面摩擦状態推定装
置。 - 【請求項2】 請求項1に記載の路面摩擦状態推定装置
であって、 前記関数決定手段は、前記取得した反力と操向量に基づ
きオンライン最小二乗法を用いて前記関数の勾配と切片
を算出するものである、路面摩擦状態推定装置。 - 【請求項3】 請求項1または2に記載の路面摩擦状態
推定装置であって、 前記操向量取得手段は、 当該車両の運動に関する特性をあらかじめ記憶している
車両特性記憶手段と、 当該車両の状態を表す所定の物理量を検出する車両状態
検出手段と、 前記車両の運動に関する特性に従って、前記車両の運動
状態を表す所定の物理量に基づき、前記車両の運動を解
析し、前記操向量を算出する操向量算出手段と、を有す
る、路面摩擦状態推定装置。 - 【請求項4】 請求項3に記載の路面摩擦状態推定装置
であって、前記操向量は横すべり角である、路面摩擦状
態推定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000370709A JP2002173044A (ja) | 2000-12-05 | 2000-12-05 | 路面摩擦状態推定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000370709A JP2002173044A (ja) | 2000-12-05 | 2000-12-05 | 路面摩擦状態推定装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2002173044A true JP2002173044A (ja) | 2002-06-18 |
Family
ID=18840547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000370709A Pending JP2002173044A (ja) | 2000-12-05 | 2000-12-05 | 路面摩擦状態推定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2002173044A (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009006940A (ja) * | 2007-06-29 | 2009-01-15 | Mitsubishi Electric Corp | 車両用操舵装置 |
KR100892480B1 (ko) | 2007-11-27 | 2009-04-10 | 현대자동차주식회사 | 조향 복원토크 추정 시스템 |
JP2009243938A (ja) * | 2008-03-28 | 2009-10-22 | Denso Corp | 荷重センサ |
US20140114547A1 (en) | 2012-04-18 | 2014-04-24 | Eaton Corporation | Method and apparatus for real time estimation of road surface friction coefficient |
-
2000
- 2000-12-05 JP JP2000370709A patent/JP2002173044A/ja active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009006940A (ja) * | 2007-06-29 | 2009-01-15 | Mitsubishi Electric Corp | 車両用操舵装置 |
JP4583413B2 (ja) * | 2007-06-29 | 2010-11-17 | 三菱電機株式会社 | 車両用操舵装置 |
KR100892480B1 (ko) | 2007-11-27 | 2009-04-10 | 현대자동차주식회사 | 조향 복원토크 추정 시스템 |
JP2009243938A (ja) * | 2008-03-28 | 2009-10-22 | Denso Corp | 荷重センサ |
US7926826B2 (en) | 2008-03-28 | 2011-04-19 | Denso Corporation | Load sensor device and steering apparatus for vehicle having the same |
JP4687739B2 (ja) * | 2008-03-28 | 2011-05-25 | 株式会社デンソー | 荷重センサ |
US20140114547A1 (en) | 2012-04-18 | 2014-04-24 | Eaton Corporation | Method and apparatus for real time estimation of road surface friction coefficient |
JP2014518533A (ja) * | 2012-04-18 | 2014-07-31 | イートン コーポレーション | 路面摩擦係数のリアルタイム推定方法及び装置 |
US8914213B2 (en) | 2012-04-18 | 2014-12-16 | Eaton Corporation | Method and apparatus for real time estimation of road surface friction coefficient |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4926715B2 (ja) | 車両を安定化させる際に車両操作者を支援するための方法及び装置 | |
US6508102B1 (en) | Near real-time friction estimation for pre-emptive vehicle control | |
JP4202872B2 (ja) | 車両用操舵装置 | |
JP4042277B2 (ja) | 車体横すべり角推定装置 | |
JP3539722B2 (ja) | 車両の路面摩擦係数推定装置 | |
JP3426000B2 (ja) | カーブ走行における車両安定性の制御方法 | |
JP3618274B2 (ja) | 車両用センサ異常検出装置 | |
JP4459561B2 (ja) | コーナリングパワー制御装置およびコーナリングパワー制御方法 | |
US6862512B2 (en) | Method and system for controlling the performance of a motor vehicle | |
US6745112B2 (en) | Method of estimating quantities that represent state of vehicle | |
JP2002012160A (ja) | 車両の路面摩擦係数推定装置 | |
US20050216155A1 (en) | Steering device | |
JP2003127888A (ja) | 物理量推定装置、路面摩擦状態推定装置、操舵角中立点推定装置、及び空気圧低下推定装置 | |
JP3970095B2 (ja) | 操舵装置 | |
US6233513B1 (en) | Method and system for computing a vehicle body slip angle in a vehicle movement control | |
JP2003306092A (ja) | 車両状態量の推定方法 | |
EP1760451A1 (en) | Method and system for road surface friction coefficient estimation | |
Tseng et al. | Technical challenges in the development of vehicle stability control system | |
JP3236003B1 (ja) | 車両の路面摩擦係数推定装置 | |
JP2021160540A (ja) | 路面μ推定装置 | |
JP4568996B2 (ja) | 操舵装置と操舵装置の抗力算出装置 | |
JP5206490B2 (ja) | 車両接地面摩擦状態推定装置及びその方法 | |
JP2003081119A (ja) | 自動車の電動パワーステアリング装置 | |
JP4363118B2 (ja) | 路面状態判定方法および路面状態判定装置 | |
JP2002173044A (ja) | 路面摩擦状態推定装置 |