特開 2019-166904 車両状態推定装置、制御装置、サスペンション制御装置、サスペンション装置、ステアリング制御装置、及びステアリング装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】特開2019-166904(P2019-166904A)

(43)【公開日】2019年10月3日

(54)【発明の名称】車両状態推定装置、制御装置、サスペンション制御装置、サスペンション装置、ステアリング制御装置、及びステアリング装置

(51)【国際特許分類】

   B60G 17/018 20060101AFI20190906BHJP

   B60G 17/015 20060101ALI20190906BHJP

   B60W 40/10 20120101ALI20190906BHJP

【ＦＩ】

   B60G17/018

   B60G17/015 A

   B60W40/10

【審査請求】未請求

【請求項の数】19

【出願形態】ＯＬ

【全頁数】38

(21)【出願番号】特願2018-54785(P2018-54785)

(22)【出願日】2018年3月22日

(71)【出願人】

【識別番号】000146010

【氏名又は名称】株式会社ショーワ

(74)【代理人】

【識別番号】110000338

【氏名又は名称】特許業務法人ＨＡＲＡＫＥＮＺＯ ＷＯＲＬＤ ＰＡＴＥＮＴ ＆ ＴＲＡＤＥＭＡＲＫ

(72)【発明者】

【氏名】水口 祐一郎

(72)【発明者】

【氏名】エラ サンディープ

(72)【発明者】

【氏名】伊澤 正樹

【テーマコード（参考）】

3D241

3D301

【Ｆターム（参考）】

3D241BA16

3D241BA51

3D241CC17

3D241CC18

3D241CE09

3D241DA56Z

3D241DA57Z

3D241DB02Z

3D241DB08Z

3D241DB10Z

3D241DB12Z

3D241DB14Z

3D241DB16Z

3D241DB22Z

3D241DB35Z

3D241DB44Z

3D301AA15

3D301DA38

3D301EA21

3D301EA22

3D301EA35

3D301EB13

3D301EB43

3D301EC06

3D301EC07

3D301EC53

(57)【要約】

【課題】車両の状態を好適に推定する。
【解決手段】車両状態推定部１２００は、車両状態に関する状態量についての演算を行う主演算部１２１０を備えており、主演算部１２１０への入力値には、車輪の接地荷重変動が含まれる。当該車輪の接地荷重変動は、タイヤ接地荷重変動算出部１２２０中の車輪有効半径算出部が算出したタイヤ有効半径Ｒ_ｅを入力とし接地荷重を出力とするマップを参照して算出される。
【選択図】図５

【特許請求の範囲】

【請求項1】

車両状態を推定する車載用の車両状態推定装置であって、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
ことを特徴とする車両状態推定装置。

【請求項2】

前記接地荷重変動算出部は、
前記車両の車輪角速度及び車輪速度の推定値の一方または両方が０である場合に、
定常値である車輪有効半径を算出することを特徴とする、請求項１に記載の車両状態推定装置。

【請求項3】

前記接地荷重変動算出部は、スリップ比の推定値を更に参照して車輪の接地荷重変動を算出することを特徴とする請求項１又は２に記載の車両状態推定装置。

【請求項4】

前記接地荷重変動算出部が算出した前記接地荷重変動を参照して参照用車輪有効半径を算出する参照用車輪有効半径演算部をさらに備え、
前記接地荷重変動算出部は、前記参照用車輪有効半径演算部が算出した参照用車輪有効半径を更に参照して車輪の接地荷重変動を算出する
ことを特徴とする請求項３に記載の車両状態推定装置。

【請求項5】

各車輪の参照用車輪有効半径をさらに推定することを特徴とする請求項４に記載の車両状態推定装置。

【請求項6】

前記１又は複数の入力値には、操作入力が含まれることを特徴とする請求項１から５の何れか１項に記載の車両状態推定装置。

【請求項7】

車輪の接地荷重変動を少なくとも参照して、車輪の路面変位を算出する路面変位算出部を更に備え、前記１又は複数の入力値には、車輪の路面変位が含まれていることを特徴とする請求項１から６の何れか１項に記載の車両状態推定装置。

【請求項8】

前記１又は複数の入力値には、車輪に関する仮想的なバネ及びダンパによる力が含まれることを特徴とする請求項１から７の何れか１項に記載の車両状態推定装置。

【請求項9】

前記演算部への前記１又は複数の入力値には、
前後Ｇセンサが検出したバネ上前後方向の加速度、
横Ｇセンサが検出したバネ上横方向の加速度、及び、
ヨーレートセンサが検出したヨーレート
が含まれることを特徴とする請求項１から８の何れか１項に記載の車両状態推定装置。

【請求項10】

車両の前後方向のバネ上速度、車両の横方向のバネ上速度、車両のバネ上鉛直速度、車両のロールレート、車両のピッチレート、車両のヨーレート、及び、車輪のバネ下鉛直速度を推定することを特徴とする請求項１から９の何れか１項に記載の車両状態推定装置。

【請求項11】

車両のロール角、車両のピッチ角、車両のヨー角、各車輪のサスストローク変位、各車輪のバネ下鉛直変位、実舵角、実舵角速度、各車輪のスリップ比、各車輪のスリップ角、各車輪のタイヤ前後力、各車輪の横力、各車輪の接地荷重変動、各車輪のタイヤ有効半径、及び各車輪での路面変位の少なくとも何れかを推定することを特徴とする請求項１０に記載の車両状態推定装置。

【請求項12】

前記出力値は、バネ上鉛直速度、ロールレート、及びピッチレートを含む請求項１から１１の何れか１項に記載の車両状態推定装置。

【請求項13】

前記入力値は、更に、制御量演算部の出力値を含み、
前記出力値は、更に、
車両の前後方向のバネ上速度、車両の横方向のバネ上速度、車両のヨーレート、車両のロール角、車両のピッチ角、車両のヨー角、各車輪のサスストローク変位、各車輪のバネ下鉛直変位、各車輪のバネ下鉛直速度、実舵角、実舵角速度の少なくとも何れかを含んでいる
ことを特徴とする請求項１２に記載の車両状態推定装置。

【請求項14】

前記接地荷重変動算出部は、
車輪角速度、
バネ上前後速度の推定値、
ピッチレートの推定値、及び、
ヨーレートの推定値
を参照して、車輪の接地荷重変動を計算する
ことを特徴とする請求項１から１３の何れか１項に記載の車両状態推定装置。

【請求項15】

車両状態を推定する車両状態推定部と、
規範車両モデルに関する演算を行う規範車両モデル演算部と、
前記車両状態推定部の出力値である推定出力から、前記規範車両モデル演算部の出力値である規範出力を減算する減算部と、
前記減算部による減算結果を積分する積分部と、
前記車両状態推定部の演算対象である推定状態量を増幅する第１の増幅部と、
前記積分部による積分結果を増幅する第２の増幅部と、
前記規範車両モデル演算部の演算対象である状態量を増幅する第３の増幅部と、
前記第１の増幅部による増幅結果、前記第２の増幅部による増幅結果、及び前記第３の増幅部による増幅結果を加算する加算部と
を備えており、
前記車両状態推定部は、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
ことを特徴とする制御装置。

【請求項16】

サスペンションの減衰力を制御するサスペンション制御装置であって、
車両状態を推定する車両状態推定部を備え、
前記車両状態推定部は、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
ことを特徴とするサスペンション制御装置。

【請求項17】

サスペンションと前記サスペンションの減衰力を制御するサスペンション制御部とを備えたサスペンション装置であって、
前記サスペンション制御部は、車両状態を推定する車両状態推定部を備え、
前記車両状態推定部は、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
いることを特徴とするサスペンション装置。

【請求項18】

運転者が操舵操作する操舵部材に対して印加するアシストトルク又は反力トルクを制御するステアリング制御装置であって、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
ことを特徴とするステアリング制御装置。

【請求項19】

運転者が操舵操作する操舵部材と、前記操舵部材に対して印加するアシストトルク又は反力トルクを制御するステアリング制御部とを備えたステアリング装置であって、
前記ステアリング制御部は、
車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、
１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、
車両モデルを用いた演算を行う演算部と、
を備え、
前記接地荷重変動算出部は、
前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、
前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する
ことを特徴とするステアリング装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、車両の状態を推定する車両状態推定装置、制御装置、サスペンション制御装置、サスペンション装置、ステアリング制御装置、及びステアリング装置に関する。

【背景技術】

【0002】

車両の乗り心地を向上させ、また、高い操縦安定性を実現するために、車両の状態を推定し、推定した状態に基づいて車両を制御する技術が知られている。また、このような技術では、車両の状態を推定するために車両モデルが用いられる。

【0003】

例えば、特許文献１には、車両の車高に関する動的なモデルに従ってあらかじめ設定された最適フィードバックゲインに基づき、車両調整部材の制御量を定める技術が開示されている。

【0004】

また、特許文献２には、車両のモデルに基づき推定ヨーレート及び規範ヨーレートを取得し、これらのヨーレートに基づいてステアリング特性を制御する技術が開示されている。また、その他の技術として例えば、特許文献３〜５が知られている。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】昭６１−１７８２１２号公報（１９８６年８月９日）

【特許文献2】特開２００４−１８９１１７号公報（２００４年７月８日公開）

【特許文献3】特開２０１６−１０７８６２号公報（２０１６年６月２０日公開）

【特許文献4】特開２０１６−１６８８８７号公報（２０１６年９月２３日公開）

【特許文献5】特開２０１４−８８８５号公報（２０１４年１月２０日公開）

【発明の概要】

【発明が解決しようとする課題】

【0006】

本発明は、車両の状態を好適に推定することのできる車両状態推定装置及び制御装置を実現することを目的とする。また、上記車両状態推定装置及び上記制御装置の推定結果を利用することで、高い乗り心地及び高い操縦安定性を実現することができるサスペンション制御装置、サスペンション装置、ステアリング制御装置、及びステアリング装置を実現することを目的とする。

【課題を解決するための手段】

【0007】

かかる目的のもと、本発明に係る車両状態推定装置は、車両状態を推定する車載用の車両状態推定装置であって、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【0008】

また、かかる目的のもと、本発明に係る制御装置は、車両状態を推定する車両状態推定部と、規範車両モデルに関する演算を行う規範車両モデル演算部と、前記車両状態推定部の出力値である推定出力から、前記規範車両モデル演算部の出力値である規範出力を減算する減算部と、前記減算部による減算結果を積分する積分部と、前記車両状態推定部の演算対象である推定状態量を増幅する第１の増幅部と、前記積分部による積分結果を増幅する第２の増幅部と、前記規範車両モデル演算部の演算対象である状態量を増幅する第３の増幅部と、前記第１の増幅部による増幅結果、前記第２の増幅部による増幅結果、及び前記第３の増幅部による増幅結果を加算する加算部と、を備えており、前記車両状態推定部は、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【0009】

また、本発明に係るサスペンション制御装置は、サスペンションの減衰力を制御するサスペンション制御装置であって、車両状態を推定する車両状態推定部を備え、前記車両状態推定部は、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【0010】

また、本発明に係るサスペンション装置は、サスペンションと前記サスペンションの減衰力を制御するサスペンション制御部とを備えたサスペンション装置であって、前記サスペンション制御部は、車両状態を推定する車両状態推定部を備え、前記車両状態推定部は、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【0011】

また、本発明に係るステアリング制御装置は、運転者が操舵操作する操舵部材に対して印加するアシストトルク又は反力トルクを制御するステアリング制御装置であって、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【0012】

また、本発明に係るステアリング装置は、運転者が操舵操作する操舵部材と、前記操舵部材に対して印加するアシストトルク又は反力トルクを制御するステアリング制御部とを備えたステアリング装置であって、前記ステアリング制御部は、車両における車輪の接地荷重の変動分である車輪の接地荷重変動を算出する接地荷重変動算出部と、１又は複数の入力値を参照して車輪有効半径を算出する車輪有効半径算出部と、車両モデルを用いた演算を行う演算部と、を備え、前記接地荷重変動算出部は、前記車輪有効半径算出部が算出した前記車輪有効半径を入力とし接地荷重を出力とするマップを参照して、前記車輪の接地荷重変動を算出し、前記演算部は、前記接地荷重変動算出部が算出した車輪の接地荷重変動を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。

【発明の効果】

【0013】

本発明に係る車両状態推定装置によれば、車両の状態を好適に推定することができる。また、本発明に係るサスペンション制御装置、サスペンション装置、ステアリング制御装置、及びステアリング装置によれば、高い乗り心地及び高い操縦安定性を実現することができる。

【図面の簡単な説明】

【0014】

【図1】本発明の実施形態１に係る車両の概略構成例を示す図である。

【図2】本発明の実施形態１に係る懸架装置における油圧緩衝装置の概略構成例を示す概略断面図である。

【図3】本発明の実施形態１に係るＥＣＵの概略構成例を示すブロック図である。

【図4】本発明の実施形態１に係る規範車両モデル演算部の構成例を示すブロック図である。

【図5】本発明の実施形態１に係る車両状態推定部の構成例を示すブロック図である。

【図6】本発明の実施形態１に係る車両状態推定部が参照するタイヤ剛性に関するマップである。

【図7】本発明の実施形態２に係るＥＣＵの概略構成例を示すブロック図である。

【図8】本発明の実施形態２に係る車両状態推定部の構成例を示すブロック図である。

【図9】本発明の実施形態２に係る車両状態推定部が備える推定状態量再構築部の構成例を示すブロック図である。

【発明を実施するための形態】

【0015】

〔実施形態１〕
以下、本発明の実施形態１について、詳細に説明する。

【0016】

（車両９００の構成）
図１は、本実施形態に係る車両９００の概略構成例を示す図である。図１に示すように、車両９００は、懸架装置（サスペンション）１００、車体２００、車輪３００、タイヤ３１０、操舵部材４１０、ステアリングシャフト４２０、トルクセンサ４３０、舵角センサ４４０、トルク印加部４６０、ラックピニオン機構４７０、ラック軸４８０、エンジン５００、ＥＣＵ（Electronic Control Unit、制御装置、ステアリング制御装置、サスペンション制御装置、サスペンション制御部）６００、発電装置７００及びバッテリ８００を備えている。ここで、懸架装置１００、及びＥＣＵ６００は、本実施形態に係るサスペンション装置を構成する。また、操舵部材４１０、ステアリングシャフト４２０、トルクセンサ４３０、舵角センサ４４０、トルク印加部４６０、ラックピニオン機構４７０、ラック軸４８０、及びＥＣＵ６００は、ステアリング装置を構成する。なお、車両９００としては、ガソリン車、ハイブリッド電気自動車（ＨＥＶ車）、電気自動車（ＥＶ車）等を挙げることができる。

【0017】

タイヤ３１０が装着された車輪３００は、懸架装置１００によって車体２００に懸架されている。車両９００は、四輪車であるため、懸架装置１００、車輪３００及びタイヤ３１０については、それぞれ４つ設けられている。

【0018】

なお、左側の前輪、右側の前輪、左側の後輪及び右側の後輪のタイヤ及び車輪をそれぞれ、タイヤ３１０Ａ及び車輪３００Ａ、タイヤ３１０Ｂ及び車輪３００Ｂ、タイヤ３１０Ｃ及び車輪３００Ｃ、並びに、タイヤ３１０Ｄ及び車輪３００Ｄとも称する。以下、同様に、左側の前輪、右側の前輪、左側の後輪及び右側の後輪にそれぞれ付随した構成を、符号「Ａ」「Ｂ」「Ｃ」及び「Ｄ」を付して表現することがある。

【0019】

懸架装置１００は、油圧緩衝装置、アッパーアーム及びロアーアームを備えている。また、油圧緩衝装置は、一例として、当該油圧緩衝装置が発生させる減衰力を調整する電磁弁であるソレノイドバルブを備えている。ただし、これは本実施形態を限定するものではなく、油圧緩衝装置は、減衰力を調整する電磁弁として、ソレノイドバルブ以外の電磁弁を用いてもよい。例えば、上記電磁弁として、電磁流体（磁性流体）を利用した電磁弁を備える構成としてもよい。

【0020】

エンジン５００には、発電装置７００が付設されており、発電装置７００によって生成された電力がバッテリ８００に蓄積される。

【0021】

運転者の操作する操舵部材４１０は、ステアリングシャフト４２０の一端に対してトルク伝達可能に接続されており、ステアリングシャフト４２０の他端は、ラックピニオン機構４７０に接続されている。

【0022】

ラックピニオン機構４７０は、ステアリングシャフト４２０の軸周りの回転を、ラック軸４８０の軸方向に沿った変位に変換するための機構である。ラック軸４８０が軸方向に変位すると、タイロッド及びナックルアームを介して車輪３００Ａ及び車輪３００Ｂが転舵される。

【0023】

トルクセンサ４３０は、ステアリングシャフト４２０に印加される操舵トルク、換言すれば、操舵部材４１０に印加される操舵トルクを検出し、検出結果を示すトルクセンサ信号をＥＣＵ６００に提供する。より具体的には、トルクセンサ４３０は、ステアリングシャフト４２０に内設されたトーションバーの捩れを検出し、検出結果をトルクセンサ信号として出力する。なお、トルクセンサ４３０として磁歪式トルクセンサを用いてもよい。

【0024】

舵角センサ４４０は、操舵部材４１０の舵角を検出し、検出結果をＥＣＵ６００に提供する。

【0025】

トルク印加部４６０は、ＥＣＵ６００から供給されるステアリング制御量に応じたアシストトルク又は反力トルクを、ステアリングシャフト４２０に印加する。トルク印加部４６０は、ステアリング制御量に応じたアシストトルク又は反力トルクを発生させるモータと、当該モータが発生させたトルクをステアリングシャフト４２０に伝達するトルク伝達機構とを備えている。

【0026】

なお、本明細書における「制御量」の具体例として、電流値、デューティー比、減衰率、減衰比等が挙げられる。

【0027】

なお、上述の説明において「トルク伝達可能に接続」とは、一方の部材の回転に伴い他方の部材の回転が生じるように接続されていることを指し、例えば、以下のいずれかの場合を含む。

【0028】

・一方の部材と他方の部材とが一体的に成形されている場合。

【0029】

・一方の部材に対して他方の部材が直接的又は間接的に固定されている場合。

【0030】

・一方の部材と他方の部材とが継手部材等を介して連動するよう接続されている場合。

【0031】

また、上記の例では、操舵部材４１０からラック軸４８０までが常時機械的に接続されたステアリング装置を例に挙げたが、これは本実施形態を限定するものではない。本実施形態に係るステアリング装置は、例えばステア・バイ・ワイヤ方式のステアリング装置であってもよい。ステア・バイ・ワイヤ方式のステアリング装置に対しても本明細書において以下に説明する事項を適用することができる。

【0032】

ＥＣＵ６００は、車両９００が備える各種の電子機器を統括制御する。より具体的には、ＥＣＵ６００は、トルク印加部４６０に供給するステアリング制御量を調整することにより、ステアリングシャフト４２０に印加するアシストトルク又は反力トルクの大きさを制御する。

【0033】

また、ＥＣＵ６００は、サスペンション制御量を供給することによって懸架装置１００を制御する。より具体的には、ＥＣＵ６００は、懸架装置１００に含まれる油圧緩衝装置が備えるソレノイドバルブに対して、サスペンション制御量を供給することによって当該ソレノイドバルブの開閉を制御する。この制御を可能とするために、ＥＣＵ６００からソレノイドバルブへ駆動電力を供給する電力線が配されている。

【0034】

また、車両９００は、車輪３００毎に設けられ各車輪３００の車輪速（車輪の車輪角速度ω）を検出する車輪速センサ３２０を備えている。また、車両９００は、車両９００の横方向の加速度を検出する横Ｇセンサ３３０、車両９００の前後方向の加速度を検出する前後Ｇセンサ３４０、車両９００のヨーレートを検出するヨーレートセンサ３５０、エンジン５００が発生させるトルクを検出するエンジントルクセンサ５１０、エンジン５００の回転数を検出するエンジン回転数センサ５２０、及びブレーキ装置が有するブレーキ液に印加される圧力を検出するブレーキ圧センサ５３０を備える構成としてもよい。これらの各種センサによる検出結果は、ＥＣＵ６００に供給される。

【0035】

なお、図示は省略するが、車両９００は、ブレーキ時の車輪ロックを防ぐためのシステムであるＡＢＳ（Antilock Brake System）、加速時等における車輪の空転を抑制するＴＣＳ（Traction Control System）、及び、旋回時のヨーモーメント制御やブレーキアシスト機能等のための自動ブレーキ機能を備えた車両挙動安定化制御システムであるＥＳＣ（Electronic Stability Control）制御可能なブレーキ装置を備えている。

【0036】

ここで、ＡＢＳ、ＴＣＳ、及びＥＳＣは、推定した車体速に応じて定まる車輪速と、車輪速センサ３２０によって検出された車輪速とを比較し、これら２つの車輪速の値が、所定の値以上相違している場合にスリップ状態であると判定する。ＡＢＳ、ＴＣＳ、及びＥＳＣは、このような処理を通じて、車両９００の走行状態に応じて最適なブレーキ制御やトラクション制御を行うことにより、車両９００の挙動の安定化を図るものである。

【0037】

また、上述した各種のセンサによる検出結果のＥＣＵ６００への供給、及び、ＥＣＵ６００から各部への制御信号の伝達は、ＣＡＮ（Controller Area Network）３７０を介して行われる。

【0038】

（懸架装置１００）
図２は、本実施形態に係る懸架装置１００における油圧緩衝装置の概略構成例を示す概略断面図である。図２に示すように、懸架装置１００は、シリンダ１０１と、シリンダ１０１内に摺動可能に設けられたピストン１０２と、ピストン１０２に固定されたピストンロッド１０３とを備えている。シリンダ１０１は、ピストン１０２によって上室１０１ａと下室１０１ｂとに仕切られており、上室１０１ａ及び下室１０１ｂは作動油によって満たされている。なお、図２において図示は省略しているが、懸架装置１００には、バーストを防止するようにガス室が設けられている。

【0039】

また、図２に示すように、懸架装置１００は、上室１０１ａと下室１０１ｂとを連通させる連通路１０４を備えており、当該連通路１０４上には、懸架装置１００の減衰力を調整するソレノイドバルブ１０５が設けられている。

【0040】

ソレノイドバルブ１０５は、ソレノイド１０５ａと、ソレノイド１０５ａによって駆動され、連通路１０４の流路断面積を変更するバルブ１０５ｂとを備えている。

【0041】

ソレノイド１０５ａはＥＣＵ６００から供給されるサスペンション制御量に応じてバルブ１０５ｂを出し入れし、それにより連通路１０４の流路断面積が変更され、懸架装置１００の減衰力が変更される。

【0042】

なお、懸架装置１００として、アクティブサスペンションやエアサスペンションを用いてもよい。

【0043】

（ＥＣＵ６００）
以下では、参照する図面を替えて、ＥＣＵ６００について具体的に説明する。図３は、ＥＣＵ６００の概略構成例を示す図である。

【0044】

図３に示すように、ＥＣＵ６００（制御装置）は、制御量演算部１０００、及び車両状態推定部（車両状態推定装置）１２００を備えている。なお、図３に示す車両各部１３００は、制御量演算部１０００による演算結果を参照して制御される車両９００の各部、及び車両９００の状態量を取得するための各種のセンサを表している。制御対象の車両９００の各部の一例として、懸架装置１００やトルク印加部４６０が挙げられ、各種センサの一例として、各車両の角速度センサ、横Ｇセンサ３３０、前後Ｇセンサ３４０、及びヨーレートセンサ３５０が挙げられる。

【0045】

（制御量演算部）
制御量演算部１０００は、図３に示すように、規範車両モデル演算部１１００、減算部１０１２、積分部１０１４、第１の増幅部１０２１、第２の増幅部１０２２、第３の増幅部１０２３、及び加算部１０２４を備えている。

【0046】

規範車両モデル演算部１１００は、入力値に対して規範用車両モデルを用いた演算を行い、演算結果である規範出力を減算部１０１２に供給する。また、規範車両モデル演算部１１００は、演算対象である種々の状態量を、規範状態量として第３の増幅部１０２３に供給する。規範車両モデル演算部１１００が出力する規範出力は、車両制御における目標値としての意味を有する。ここで、規範出力は、上記演算対象である種々の状態量の少なくとも一部を構成する。

【0047】

規範車両モデル演算部１１００への入力の一例として、図３に示すように、車輪の車輪角速度ω_fl〜ω_rr、操作入力、及び、後述の車両状態推定部１２００において推定する各車輪の接地荷重変動dF_z0fl〜dF_z0rrから算出する各車輪の路面変位z_0fl〜z_0rrが挙げられる。ここで、操作入力には、操舵部材４１０の操舵角が含まれる。ここで、各車輪の接地荷重変動dF_z0fl〜dF_z0rrとは、各車輪の接地荷重の変動分のことを指す。

【0048】

なお、本明細書において、添え字「fl」、「fr」、「rl」、「rr」はそれぞれ左前、右前、左後、右後の車輪に関するものであることを明示するための添え字である。上記のように、これらの添え字をまとめて、例えば「ω_fl〜ω_rr」のように表示することもある。いずれのパラメータについても同様である。

【0049】

また、規範車両モデル演算部１１００が減算部１０１２に供給する規範出力の例として、車体２００の前後方向のバネ上速度u、横方向のバネ上速度v、バネ上鉛直速度w、ロールレートp、ピッチレートq、ヨーレートr、ロール角phi、ピッチ角theta、ヨー角psi、各車輪のサスストローク変位SusSt_fl〜SusSt_rr、各車輪のバネ下鉛直変位z_1flm〜z_1rrm、各車輪のバネ下鉛直速度w_1flm〜w_1rrm、実舵角δ、及び実舵角速度dδの少なくとも何れか又はそれらの組み合わせによって表現することができる物理量が挙げられる。また、規範車両モデル演算部１１００が第３の増幅部１０２３に供給する規範状態量の例として、車体２００の前後方向のバネ上速度u、横方向のバネ上速度v、バネ上鉛直速度w、ロールレートp、ピッチレートq、ヨーレートr、ロール角phi、ピッチ角theta、ヨー角psi、各車輪のサスストローク変位SusSt_fl〜SusSt_rr、各車輪のバネ下鉛直変位z_1flm〜z_1rrm、各車輪のバネ下鉛直速度w_1flm〜w_1rrm、実舵角δ、及び実舵角速度dδの少なくとも何れかが挙げられる。なお、規範車両モデル演算部１１００の具体的な構成は後述する。

【0050】

減算部１０１２は、後述する車両状態推定部１２００からの推定出力を取得し、取得した推定出力から、規範車両モデル演算部１１００が出力する規範出力を減算し、減算結果を積分部１０１４に供給する。

【0051】

車両状態推定部１２００が推定し、減算部１０１２に供給する推定出力の一例として、車体２００の前後方向のバネ上速度、横方向のバネ上速度、バネ上鉛直速度、ロールレート、ピッチレート、ロール角、ピッチ角、ヨー角、各車輪のサスストローク変位、各車輪のバネ下鉛直変位、各車輪のバネ下鉛直速度、実舵角、及び実舵角速度の少なくとも何れか又はそれらの組み合わせによって表現することができる物理量が挙げられる。また、車両状態推定部１２００が推定し、第１の増幅部１０２１に供給する推定状態量の一例として、車体２００の前後方向のバネ上速度u、横方向のバネ上速度v、バネ上鉛直速度w、ロールレートp、ピッチレートq、ヨーレートr、ロール角phi、ピッチ角theta、ヨー角psi、各車輪のサスストローク変位SusSt_fl〜SusSt_rr、各車輪のバネ下鉛直変位z_1flm〜z_1rrm、各車輪のバネ下鉛直速度w_1flm〜w_1rrm、実舵角δ、及び実舵角速度dδの少なくとも何れかが挙げられる。

【0052】

なお、本発明の車両状態推定装置により、推定可能な物理量は、車両状態推定部１２００における演算対象の「各状態量」、「各状態量の何れかの組み合わせで表現することができる物理量」、「タイヤモデル演算部１２４０で計算されるタイヤ前後力及びタイヤ横力」、「スリップ算出部１２３０で計算されるスリップ比及びスリップ角」、「タイヤ接地荷重変動算出部１２２０の出力を利用し得られるタイヤ接地荷重」、「タイヤ接地荷重変動算出部１２２０で算出するタイヤ有効半径」、「参照用タイヤ有効半径演算部１２７０で計算される参照用タイヤ有効半径」、及び「路面変位算出部１２８０で算出される路面変位」である。

【0053】

積分部１０１４は、減算部１０１２による減算結果を積分する。積分した結果は、第２の増幅部１０２２に供給される。

【0054】

第１の増幅部１０２１は、車両状態推定部１２００から供給される推定状態量を、増幅係数K1を用いて増幅し、増幅した結果を加算部１０２４に供給する。

【0055】

第２の増幅部１０２２は、積分部１０１４による積分結果を、増幅係数K2を用いて増幅し、増幅した結果を加算部１０２４に供給する。

【0056】

第３の増幅部１０２３は、規範車両モデル演算部１１００から供給される規範状態量を、増幅係数K3を用いて増幅し、増幅結果を加算部１０２４に供給する。

【0057】

加算部１０２４は、第１の増幅部１０２１による増幅結果と、第２の増幅部１０２２による増幅結果と、第３の増幅部１０２３による増幅結果とを加算し、加算した結果を車両状態推定部１２００、及び車両各部１３００に供給する。加算部１０２４による加算結果は、制御量演算部１０００による演算結果を表している。

【0058】

制御量演算部１０００は、車両状態推定部１２００の出力値である推定出力から、規範車両モデル演算部１１００の出力値である規範出力を減算する減算部１０１２と、減算部１０１２による減算結果を積分する積分部１０１４と、車両状態推定装置１２００による演算対象である推定状態量を増幅する第１の増幅部１０２１と、積分部１０１４による積分結果を増幅する第２の増幅部１０２２と、規範車両モデル演算部１１００の演算対象である規範状態量を増幅する第３の増幅部１０２３と、第１の増幅部１０２１による増幅結果、第２の増幅部１０２２による増幅結果、及び第３の増幅部１０２３による増幅結果を加算する加算部１０２４を備えているので、車両９００は、偏差なく規範モデル特性に追従することができる。

【0059】

また、制御量演算部１０００が、積分部１０１４を備えていることにより、車両９００は、偏差なく規範モデル特性に追従することができる。

【0060】

（規範車両モデル演算部）
続いて、図４を参照して規範車両モデル演算部１１００の構成について具体的に説明する。図４は、規範車両モデル演算部１１００の構成例を示すブロック図である。図４に示すように、規範車両モデル演算部１１００は、主演算部１１１０、規範車両モデルタイヤ接地荷重算出部１１２０、スリップ算出部１１３０、タイヤモデル演算部１１４０、操縦安定性・乗心地制御部１１５０、及び、参照用タイヤ有効半径演算部１１６０を備えている。

【0061】

（主演算部）
主演算部１１１０は、１又は複数の入力値を参照し、車両状態に関する状態量についての線形演算を行うことによって１又は複数の出力値を算出する。主演算部１１１０のことを単に演算部と呼称することもある。

【0062】

主演算部１１１０は、図４に示すように、第１の入力行列演算部１１１１、第２の入力行列演算部１１１２、第３の入力行列演算部１１１３、第４の入力行列演算部１１１８、加算部１１１４、積分部１１１５、システム行列演算部１１１６、観測行列演算部１１１７を備えている。ここで、第１の入力行列演算部１１１１、第２の入力行列演算部１１１２、第３の入力行列演算部１１１３、及び、第４の入力行列演算部１１１８を第１の演算部とも呼称する。

【0063】

路面入力に対する入力行列B0に関する演算を行う第１の入力行列演算部１１１１には、一例として路面変位（鉛直方向の変位）z_0fl、z_0fr、z_0rl、z_0rrが入力される。ここで、路面変位z_0fl、z_0fr、z_0rl、z_0rrとしては、後述する車両状態推定部１２００にて算出したものを用いる。

【0064】

第１の入力行列演算部１１１１は、入力された路面変位z_0fl〜z_0rrに対して、入力行列B0を演算し、演算した結果を加算部１１１４に供給する。

【0065】

操作量に対する入力行列B1に関する演算を行う第２の入力行列演算部１１１２は、一例として、操舵部材４１０の操舵角に対して入力行列B1を演算し、演算した結果を加算部１１１４に供給する。

【0066】

タイヤ前後/横力に対する入力行列B2に関する演算を行う第３の入力行列演算部１１１３は、後述するタイヤモデル演算部１１４０から供給される各車輪のタイヤ前後力F_x0fl〜F_x0rr、及び各車輪のタイヤ横力F_y0fl〜F_y0rrに対して入力行列B2を演算し、演算した結果を加算部１１１４に供給する。

【0067】

制御結果を反映した、追加サスペンション力、アシストトルク、及び反力トルクに対する入力行列B3に関する演算を行う第４の入力行列演算部１１１８は、後述する操縦安定性・乗心地制御部１１５０の出力に対して入力行列B3を演算し、演算した結果を加算部１１１４に供給する。

【0068】

加算部１１１４は、第１の入力行列演算部１１１１、第２の入力行列演算部１１１２、第３の入力行列演算部１１１３、第４の入力行列演算部１１１８、及び後述するシステム行列演算部１１１６からの出力を加算し、加算結果を積分部１１１５に供給する。

【0069】

積分部１１１５は、加算部１１１４から供給される加算結果を積分する。積分部１１１５による積分結果は、前述した第３の増幅部１０２３、システム行列演算部１１１６、及び、観測行列演算部１１１７に供給される。また、積分部１１１５による積分結果のうち、z_1flm〜z_1rrmは、規範車両モデルタイヤ接地荷重算出部１１２０に供給される。

【0070】

システム行列演算部（第２の演算部）１１１６は、積分部１１１５による積分結果に対して、システム行列Aを演算し、演算した結果を加算部１１１４に供給する。

【0071】

観測行列演算部（第３の演算部）１１１７は、積分部１１１５による積分結果に対して、観測行列Cを演算し、演算した結果を規範出力として、前述した減算部１０１２に供給する。また、観測行列Cを演算した結果はスリップ算出部１１３０にも供給される。

【0072】

なお、主演算部１１１０が備える各部における演算は、線形演算として実行される。したがって、上記の構成を有する主演算部１１１０によれば、１又は複数の入力値を参照した車両状態に関する状態量についての線形演算を好適に行うことができる。

【0073】

また、主演算部１１１０への入力は上記の例に限られるものではなく、例えば、
・操舵トルク
・各車輪の車輪角速度
・各車輪の実舵角
・各車輪の駆動トルク
の少なくとも何れかを主演算部１１１０に入力する構成とし、主演算部１１１０がこれらの入力値に対する線形演算を実行する構成としてもよい。この場合、例えば、主演算部１１１０が、各システム行列A、入力行列B、及び観測行列Cによって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、上記の入力を参照して各車両モデルを切り替える構成とすることができる。

【0074】

また、車両９００が、積載量検知手段を備える構成とし、主演算部１１１０には、当該積載量検知手段による検出値を入力する構成としてもよい。この場合、例えば、主演算部１１１０が、各積載量に応じたシステム行列A、入力行列B、及び観測行列Cによって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、積載量検知手段による検出値に応じて各車両モデルを切り替える構成とすることができる。なお、上記積載量検知手段は、センサで積載量を検知する構成であってもよいし、センサレスで積載量を検知する構成であってもよい。

【0075】

主演算部１１１０への入力は、
・ヨーレート
・前後Ｇ
・横Ｇ
・ブレーキ圧
・ＥＳＣフラグ、ＴＣＳフラグ、ＡＢＳフラグ
・エンジントルク
・エンジン回転数
の少なくとも何れかをさらに含む構成としてもよい。この場合、例えば、主演算部１１１０が、各システム行列A、入力行列B、及び観測行列Cによって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、上記の入力を参照して各車両モデルを切り替える構成とすることができる。

【0076】

（主演算部による演算対象の状態量の例）
主演算部１１１０による演算対象の状態量の一例を状態量ベクトルxとして表現すれば以下の通りである。

【0077】

以下、x方向は、車両９００の進行方向（前後方向）を示し、z方向は鉛直方向を示し、y方向はx方向及びz方向の双方に垂直な方向（横方向）のことを示す。

【0078】

【数1】

【0079】

ここで、
u、v、wは、車体２００のバネ上速度のx、y、z方向成分であり、
p、q、rは、車体２００のバネ上角速度のx軸、y軸、z軸方向成分、すなわち、ロールレート、ピッチレート、及び、ヨーレートである。また、

【0080】

【数2】

【0081】

はそれぞれオイラー角の３成分（それぞれ、phi、theta、psiとも表記する。また、phiはロール角、thetaはピッチ角、psiはヨー角を表す。）であり、
SusSt_fl〜SusSt_rrは、各車輪のサスストローク変位であり、上記は、バネ上と同じ動きをするボディ座標系で観測した状態量である。

【0082】

また、z_1flm〜z_1rrmは、各車輪のバネ下鉛直変位であり、w_1flm〜w_1rrmは、各車輪のバネ下鉛直速度である。ただし、z_1flm〜z_1rrm、w_1flm〜w_1rrmは、ｘ方向及びｙ方向の並進運動、並びにｚ軸方向の回転運動（ヨー）のみ、バネ上と同じ動きをする座標系で観測した状態量である。

【0083】

また、δは、実舵角であり、
dδは、実舵角速度である。

【0084】

なお、実舵角δ、及び実舵角速度dδは、各車輪３００それぞれについて個別に設定する構成としてもよいが、本明細書中では、一例として前輪のみについて設定されるものとし、かつ、タイヤの切戻りは考慮していない。

【0085】

（主演算部出力の例）
主演算部１１１０が出力する規範出力の種類は、観測行列Cをどのように選ぶかによって決定される。一例として、主演算部１１１０が出力する規範出力を特定状態量ベクトルyとして表現すれば以下のように、車体２００の前後方向のバネ上速度u、横方向のバネ上速度v、バネ上鉛直速度w、ロールレートp、ピッチレートq、ヨーレートr、ロール角phi、ピッチ角theta、ヨー角psi、各車輪のサスストローク変位SusSt_fl〜SusSt_rr、各車輪のバネ下鉛直変位z_1flm〜z_1rrm、各車輪のバネ下鉛直速度w_1flm〜w_1rrm、実舵角δ、及び実舵角速度dδを含んでいる。

【0086】

【数3】

【0087】

また、主演算部１１１０が出力する規範出力は、上述した状態量ベクトルxに含まれる状態量の何れか又はそれらの組み合わせによって表現することができる物理量である。

【0088】

（主演算部による演算対象の運動方程式の一例）
主演算部１１１０による演算対象の運動方程式の一例を示せば以下の通りである。また、各物理量の上に付されたドット「・」は時間微分を表す。

【0089】

・バネ上並進及び回転運動に関する下記運動方程式

【0090】

【数4】

【0091】

・オイラー角に関する下記運動方程式

【0092】

【数5】

【0093】

・バネ下鉛直運動に関する下記運動方程式

【0094】

【数6】

【0095】

・実舵に関する下記運動方程式（ただし、タイヤ切戻りを考慮せず）

【0096】

【数7】

【0097】

・各車輪のサスペンション力

【0098】

【数8】

【0099】

・サスストローク変位

【0100】

【数9】

【0101】

上記の運動方程式及び関係式において、
mは車両のバネ上質量（すなわち車体２００の質量）であり、
F_x、F_y、F_zは、車両のバネ上（すなわち車体２００）に作用するx、y、z方向の力であり、
M_x、M_y、M_zは、車両のバネ上に作用するx、y、z軸に関するモーメントであり、
I_x、I_y、I_zは、車両のバネ上のx、y、z軸に関する慣性モーメントであり、
I_zxは、ｙ軸の慣性乗積である。

【0102】

k_2f、k_2rは、前輪及び後輪に関するスプリングのバネ定数であり、
c_2f、c_2rは、前輪及び後輪に関するダンパの減衰係数であり、
F_contfl、F_contfr、F_contrl、F_contrrは、制御の結果、追加されるサスペンション力である。

【0103】

また、F_zfl、F_zfr、F_zrl、F_zrrは各車輪のサスペンション力であり、
z_fl〜z_rrは、各車輪におけるバネ上鉛直変位であり、
w_fl〜w_rrは、各車輪におけるバネ上鉛直速度であり、
z_1fl〜z_1rrは、各車輪におけるバネ下鉛直変位であり、
w_1fl〜w_1rrは、各車輪におけるバネ下鉛直速度であり、
上記は、バネ上と同じ動きをするボディ座標系で観測した状態量である。

【0104】

m1はバネ下質量、k₁はタイヤ上下変位に対する剛性である。また、F_zflm〜F_zrrmは、各輪のバネ下にかかるサスペンション反力であり、x方向及びy方向の並進運動、並びにz軸方向の回転運動（ヨー）のみ、バネ上と同じ動きをする座標系で、観測した物理量である。

【0105】

また、αは操舵角であり、
I_sは、キングピン軸周りの車輪慣性モーメントであり、
c_sは、キングピン等価粘性摩擦係数であり、
k_sは、キングピン軸周りの等価弾性係数である。

【0106】

M_contは、制御の結果、追加されるアシストトルクである。

【0107】

上記に挙げた運動方程式以外にも、主演算部１１１０は車輪回転運動に関する運動方程式を演算対象としてもよい。また、上記に挙げた運動方程式に現れる物理量は、互いに結びつけられる複数の関係式（例えば各座標変換等）が存在しており、各運動方程式は、これらの関係式と共に解かれる。

【0108】

（運動方程式の線形化と主演算部への実装）
上述した運動方程式は一般に非線形であり、以下のように表現できる。

【0109】

【数10】

【0110】

ここで、xは状態量を示すベクトルであり、f(x)、g(x)は、xの関数であり、ベクトルとして表現される。

【0111】

上記非線形の運動方程式をテイラー展開し、ヤコビ行列に各状態量の初期値を代入することにより以下に示す行列Aが得られ、同様の方法で以下に示す行列B、Cが得られる。

【0112】

その結果、線形化された運動方程式が、以下のように状態空間にて表現される。

【0113】

【数11】

【0114】

ここで、行列Aは、上述したシステム行列Aに対応し、行列Bは、上述した入力行列B0、B1、B2、B3に対応し、行列Cは、上述した観測行列Cに対応する。

【0115】

以上の説明から、図４に示した主演算部１１１０は、対象とする運動方程式を線形的に演算する構成であることが分かる。

【0116】

（規範車両モデル演算部のその他の構成）
続いて、規範車両モデル演算部１１００が備える構成のうち、主演算部１１１０以外の構成について説明する。

【0117】

（規範車両モデルタイヤ接地荷重算出部）
規範車両モデルタイヤ接地荷重算出部１１２０は、積分部１１１５による積分結果として、算出される状態量の一部である各車輪のバネ下鉛直変位z_1flm〜z_1rrm、後述する車両状態推定部１２００にて算出する各車輪の路面変位z_0fl〜z_0rr、及び、接地荷重の定常値F_z0constfl〜F_z0constrrを用いて、規範車両モデルタイヤ接地荷重F_z0fl〜F_z0rrを以下の式を用いて算出する。接地荷重の定常値F_z0constfl〜F_z0constrrは、例えば、フラット路面での静止状態における接地荷重値を予め、例えば出荷前に、計測することにより求められる。

【0118】

F_z0fl = -k₁(z_1flm-z_0fl) + F_z0constfl
F_z0rl = -k₁(z_1rlm-z_0rl) + F_z0constrl
F_z0fr = -k₁(z_1frm-z_0fr) + F_z0constfr
F_z0rr = -k₁(z_1rrm-z_0rr) + F_z0constrr
算出した規範車両モデルタイヤ接地荷重F_z0fl〜F_z0rrは、タイヤモデル演算部１１４０及び参照用タイヤ有効半径演算部１１６０に入力される。参照用タイヤ有効半径演算部１１６０は、規範車両モデルタイヤ接地荷重算出部１１２０にて算出する規範車両モデルタイヤ接地荷重F_z0fl〜F_z0rrを用い、例えば後述する参照用タイヤ有効半径に係る第二の式を用いて、参照用タイヤ有効半径R_efl〜R_errを算出する。

【0119】

（スリップ算出部）
スリップ算出部１１３０は、観測行列演算部１１１７による演算結果、参照用タイヤ有効半径演算部１１６０による演算結果、及び車輪速センサ３２０が検出した各車輪の車輪角速度ω_fl〜ω_rrを参照する。そして、各車輪のスリップ比s_fl〜s_rrを算出し、観測行列演算部１１１７による演算結果として各車輪のスリップ角β_fl〜β_rrを算出する。さらに、算出した結果を積算部１１３１及びタイヤモデル演算部１１４０に供給する。

【0120】

（タイヤモデル演算部）
タイヤモデル演算部１１４０は、主演算部１１１０による演算結果の少なくとも一部を直接的または間接的に参照した非線形演算を行う。図４に示す例では、タイヤモデル演算部１１４０は、観測行列演算部１１１７による演算により得られる各車輪のスリップ比s_fl〜s_rr、各車輪のスリップ角β_fl〜β_rr、及び規範車両モデルタイヤ接地荷重算出部１１２０が算出した規範車両モデルタイヤ接地荷重F_z0fl〜F_z0rrを参照して非線形演算を行う。すなわち、図４に示す例では、タイヤモデル演算部１１４０は、主演算部１１１０による演算結果の少なくとも一部を間接的に参照した非線形演算を行う。

【0121】

より具体的には、タイヤモデル演算部１１４０は、各車輪のスリップ比s_fl〜s_rr、各車輪のスリップ角β_fl〜β_rr、及び各車輪の規範車両モデルタイヤ接地荷重F_z0fl〜F_z0rrを参照する。そして、タイヤモデルに関する演算式を用いて、各車輪のタイヤ前後力F_x0fl〜F_x0rr、及び各車輪のタイヤ横力F_y0fl〜F_y0rrを算出する。タイヤモデル演算部１１４０による具体的な演算式は、本実施形態を限定するものではないが、例えば、左前輪を代表すると、一般に知られている近似式

【0122】

【数12】

【0123】

を用いることができる。ここで、第１式におけるF_Px0flは、直進時の左前輪のタイヤ前後力を表している。各変数は、タイヤの特性やF_z0flに依存する値である。第２式におけるF_Py0flは、タイヤ前後力を伴わない際のタイヤ横力を表している。

【0124】

ここで、本実施形態に係る規範車両モデル演算部１１００において、主演算部１１１０は線形演算を行い、タイヤモデル演算部１１４０は主演算部１１１０による演算結果の少なくとも一部を直接的または間接的に参照した非線形演算を行う。このように線形演算部と非線形演算部とを分離させた構成を採用することにより、車両モデルを用いた状態量の演算を好適に行うことができる。

【0125】

また、タイヤモデル演算部１１４０は、タイヤモデルに基づく非線形演算を行うので、非線形演算を線形演算から好適に分離することができる。

【0126】

また、上述のように、第３の入力行列演算部１１１３は、タイヤモデル演算部１１４０による非線形演算結果を入力として取り込むので、主演算部１１１０による線形演算に非線形演算結果を好適に取り込むことができる。したがって、主演算部１１１０は線形演算を行いつつ精度の高い演算を行うことができる。

【0127】

（操縦安定性・乗心地制御部１１５０）
操縦安定性・乗心地制御部１１５０は、規範車両モデル各部を制御するための制御量を決定し、各部へ供給するために、観測行列演算部１１１７が出力する規範出力に対して作用する。操縦安定性・乗心地制御部１１５０による出力は、第４の入力行列演算部１１１８に供給され、入力行列B3が演算される。

【0128】

一例として、操縦安定性・乗心地制御部１１５０は、スカイフック制御、ロール姿勢制御、ピッチ姿勢制御及びバネ下制御、実舵角制御処理及び、制御量選択処理を行う。

【0129】

ここで、スカイフック制御とは、路面の凹凸を乗り越える際の規範車両モデルの動揺を抑制し、乗り心地を高める乗り心地制御（制振制御）のことである。

【0130】

スカイフック制御では、一例として、バネ上鉛直速度、４輪のストローク速度、ピッチレート、及びロールレートを参照して、スカイフック目標制御量を決定し、その結果を制御量選択処理の対象とする。

【0131】

ロール姿勢制御では、転舵時ロールレート、及び操舵角を参照して、ロール目標制御量を算出し、その結果を制御量選択処理の対象とする。

【0132】

ピッチ姿勢制御では、加減速時ピッチレートを参照してピッチ制御を行い、ピッチ目標制御量を算出し、その結果を制御量選択処理の対象とする。

【0133】

バネ下制御では、各車輪のバネ下鉛直速度を参照して、バネ下制振制御目標制御量を決定し、決定結果を制御量選択処理の対象とする。実舵角制御では、実舵角を参照して、目標制御量を算出し、その結果を制御量選択処理の対象とする。

【0134】

制御量選択処理では、スカイフック目標制御量、ロール目標制御量、ピッチ目標制御量、バネ下制振制御目標制御量、及び実舵角目標制御量のうち、最も大きい値を有する目標制御量を選択し出力してもよい。

【0135】

（車両状態推定部）
続いて、参照する図面を替えて、車両状態推定部１２００について具体的に説明する。車両状態推定部１２００は、入力値に対して推定用車両モデルを用いた演算を行い、演算結果である推定出力を減算部１０１２に供給する。また、車両状態推定部１２００は、演算対象である種々の状態量を、第１の増幅部１０２１に供給する。

【0136】

図５は、車両状態推定部１２００の構成例を示すブロック図である。図５に示すように、車両状態推定部１２００は、主演算部１２１０、タイヤ接地荷重変動算出部（接地荷重変動算出部）１２２０、スリップ算出部１２３０、タイヤモデル演算部１２４０、仮想バネ・ダンパ力算出部１２６０、参照用タイヤ有効半径演算部（参照用車輪有効半径演算部）１２７０、路面変位算出部１２８０、積分部１２３１、積分部１２７１、及び加算部１２７５を備えている。

【0137】

（主演算部）
推定用車両モデルを用いた演算を行う主演算部１２１０は、１又は複数の入力値を参照し、車両状態に関する状態量についての線形演算を行うことによって１又は複数の出力値を算出する。主演算部１２１０のことを単に演算部と呼称することもある。

【0138】

主演算部１２１０は、図５に示すように、第１の入力行列演算部１２１１、第２の入力行列演算部１２１２、第３の入力行列演算部１２１３、第４の入力行列演算部１２１８、第５の入力行列演算部１２１９、加算部１２１４、積分部１２１５、システム行列演算部１２１６、観測行列演算部１２１７を備えている。ここで、第１の入力行列演算部１２１１、第２の入力行列演算部１２１２、第３の入力行列演算部１２１３、第４の入力行列演算部１２１８、及び、第５の入力行列演算部１２１９を第１の演算部とも呼称する。

【0139】

接地荷重変動入力に対する入力行列B00'に関する演算を行う第１の入力行列演算部１２１１には、タイヤ接地荷重変動算出部１２２０により得られる接地荷重変動dF_z0fl〜dF_z0rrが入力される。

【0140】

第１の入力行列演算部１２１１は、入力された接地荷重変動dF_z0fl〜dF_z0rrに対して、入力行列B00'を演算し、演算した結果を加算部１２１４に供給する。

【0141】

操作量に対する入力行列B1'に関する演算を行う第２の入力行列演算部１２１２は、一例として、操舵部材４１０の操舵角に対して入力行列B1'を演算し、演算した結果を加算部１２１４に供給する。なお、第２の入力行列演算部１２１２が演算する入力行列B1'は、第２の入力行列演算部１１１２が演算する入力行列B1と同じであってもよいし、異なっていてもよい。

【0142】

タイヤ前後/横力に対する入力行列B2'に関する演算を行う第３の入力行列演算部１２１３は、後述するタイヤモデル演算部１２４０から供給される各車輪のタイヤ前後力F_x0fl〜F_x0rr、及び各車輪のタイヤ横力F_y0fl〜F_y0rrに対して入力行列B2'を演算する。そして、演算した結果を加算部１２１４に供給する。なお、第３の入力行列演算部１２１３が演算する入力行列B2'は、第３の入力行列演算部１１１３が演算する入力行列B2と同じであってもよいし、異なっていてもよい。

【0143】

また、制御量演算部１０００の出力に対する入力行列B4'に関する演算を行う第４の入力行列演算部１２１８には、制御量演算部１０００の出力が入力される。第４の入力行列演算部１２１８は、制御量演算部１０００の出力に対して、入力行列B4'を演算し、演算した結果を加算部１２１４に供給する。

【0144】

仮想バネ・ダンパ力に対する入力行列B5'に関する演算を行う第５の入力行列演算部１２１９は、後述する仮想バネ・ダンパ力算出部１２６０の出力に対して入力行列B5'を演算し、演算した結果を加算部１２１４に供給する。

【0145】

加算部１２１４は、第１の入力行列演算部１２１１、第２の入力行列演算部１２１２、第３の入力行列演算部１２１３、第４の入力行列演算部１２１８、第５の入力行列演算部１２１９、及び後述するシステム行列演算部１２１６からの出力を加算する。そして、加算結果を積分部１２１５に供給する。

【0146】

積分部１２１５は、加算部１２１４から供給される加算結果を積分する。積分部１２１５による積分結果は、推定状態量として出力されると共に、システム行列演算部１２１６、観測行列演算部１２１７及び、仮想バネ・ダンパ力算出部１２６０に供給される。

【0147】

システム行列演算部（第２の演算部）１２１６は、積分部１２１５による積分結果に対して、システム行列A'を演算し、演算した結果を加算部１２１４に供給する。

【0148】

観測行列演算部（第３の演算部）１２１７は、積分部１２１５による積分結果に対して、観測行列C'を演算し、演算した結果を推定出力として、前述した減算部１０１２に供給する。また、観測行列C'を演算した結果はスリップ算出部１２３０にも供給される。

【0149】

なお、主演算部１２１０が備える各部における演算は、線形演算として実行される。したがって、上記の構成を有する主演算部１２１０によれば、１又は複数の入力値を参照した車両状態に関する状態量についての線形演算を好適に行うことができる。

【0150】

また、主演算部１１１０と同様に、主演算部１２１０への入力は上記の例に限られるものではなく、例えば、
・操舵トルク
・各車輪の車輪角速度
・各車輪の実舵角
・各車輪の駆動トルク
の少なくとも何れかを主演算部１２１０に入力する構成とし、主演算部１２１０がこれらの入力値に対する線形演算を実行する構成としてもよい。この場合、例えば、主演算部１２１０が、各システム行列A'、入力行列B'、及び観測行列C'によって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、上記の入力を参照して各車両モデルを切り替える構成とすることができる。

【0151】

また、車両９００が、積載量検知手段を備える構成とし、主演算部１２１０には、当該積載量検知手段による検出値を入力する構成としてもよい。この場合、例えば、主演算部１２１０が、各積載量に応じたシステム行列A'、入力行列B'、及び観測行列C'によって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、積載量検知手段による検出値に応じて各車両モデルを切り替える構成とすることができる。なお、上記積載量検知手段は、センサで積載量を検知する構成であってもよいし、センサレスで積載量を検知する構成であってもよい。

【0152】

主演算部１２１０への入力は、
・ヨーレート
・前後Ｇ
・横Ｇ
・ブレーキ圧
・ＥＳＣフラグ、ＴＣＳフラグ、ＡＢＳフラグ
・エンジントルク
・エンジン回転数
の少なくとも何れかをさらに含む構成としてもよい。この場合、例えば、主演算部１２１０が、各システム行列A'、入力行列B'、及び観測行列C'によって表現される各車両モデルを切り替える車両モデル切替部を備える構成とし、当該車両モデル切替部が、上記の入力を参照して各車両モデルを切り替える構成とすることができる。

【0153】

（主演算部による演算対象の状態量の例）
主演算部１２１０による演算対象の状態量は、主演算部１１１０による演算対象の状態量と同様であるのでここでは詳細な説明を省略する。なお、主演算部１２１０が出力する規範出力は、主演算部１１１０と同様に、上述した状態量ベクトルxに含まれる状態量の何れか又はそれらの組み合わせによって表現することができる物理量である。

【0154】

（運動方程式の線形化と主演算部への実装）
運動方程式の線形化と主演算部１２１０への実装は、主演算部１１１０への実装において説明した事項と同様である箇所については、詳細な説明を省略する。車両状態推定部のバネした鉛直運動に関する運動方程式については、下式とする。なお、主演算部１１１０への実装において説明した、線形化された運動方程式における行列A、Cは、主演算部１２１０における行列A'、C'に対応する。また、線形化された運動方程式における行列Bは、主演算部１２１０における行列B00'、B1'、B2'、B4'、B5'に対応する。

【0155】

【数13】

【0156】

以上の説明から、図５に示した主演算部１２１０は、対象とする運動方程式を線形的に演算する構成であることが分かる。

【0157】

（車両状態推定部のその他の構成）
続いて、車両状態推定部１２００が備える構成のうち、主演算部１２１０以外の構成について説明する。

【0158】

（スリップ算出部）
スリップ算出部１２３０は、観測行列演算部１２１７による演算結果、参照用タイヤ有効半径演算部１２７０による演算結果、及び、車輪速センサ３２０が検出した各車輪の車輪角速度ω_fl〜ω_rrを参照する。そして、各車輪のスリップ比s_fl〜s_rrを算出し、観測行列演算部１２１７による演算結果として各車輪のスリップ角β_fl〜β_rrを算出し、算出した結果を積分部１２３１、及びタイヤモデル演算部１２４０に供給する。積算部１２３１は、スリップ算出部１２３０から供給される算出結果を積分し、積分後のスリップ比ｓを算出する。積算部１２３１により算出された積分後のスリップ比ｓは、タイヤ接地荷重変動算出部１２２０に供給される。

【0159】

（タイヤ接地荷重変動算出部）
タイヤ接地荷重変動算出部１２２０は、積分部１２１５により演算された前後方向のバネ上速度uの推定値、ピッチレートｑの推定値、及びヨーレートｒの推定値、車輪速センサ３２０が検出した各車輪の車輪角速度ω_fl〜ω_rrを参照する。そして、各車輪のタイヤ有効半径（車輪有効半径）R_efl〜R_errを算出する。各車輪のタイヤ有効半径R_efl〜R_errは、各車輪の車輪速度V_fl〜V_rrを各車輪の車輪角速度ω_fl〜ω_rrで除することにより求まる。

【0160】

【数14】

【0161】

また、各車輪の車輪速度V_fl〜V_rrは、車体２００の前後方向のバネ上速度u、ピッチレートq、ヨーレートr、及び各車輪のスリップ比s_fl〜s_rrを用いて、以下のように表される。

【0162】

【数15】

【0163】

ここで、h₀は、路面から車体２００の重心までの距離を表しており、tr_fは、車体２００のフロントのトレッド幅に０．５を乗じたものを表しており、tr_rは、車体２００のリアのトレッド幅に０．５を乗じたものを表している。このように、タイヤ接地荷重変動算出部１２２０は、車輪角速度ω、バネ上前後速度u、ピッチレートq及びヨーレートrを参照して各車輪のタイヤ有効半径R_eを算出する。このように、タイヤ接地荷重変動算出部１２２０は、１又は複数の入力値を参照してタイヤ有効半径R_eを算出する車輪有効半径算出部を兼ねている。

【0164】

さらに、タイヤ接地荷重変動算出部１２２０は、fl、rl、fr及びrrのそれぞれについて、入力された車輪角速度ω及び算出した車輪速度Vの少なくとも一方がゼロであるか否かを判定する。入力されたω及び算出したVの一方または両方がゼロでない場合には、タイヤ接地荷重変動算出部１２２０は、算出したＲ_ｅを真とする。fl、rl、fr及びrrのいずれかにおける入力されたω及び算出したＶの一方または両方がゼロである場合には、タイヤ接地荷重変動算出部１２２０は、対応するfl、rl、fr及びrrのR_econstを真とする。R_econstは、定常値であるR_eであり、例えば、接地荷重F_z0constとなるときのR_eの値である。

【0165】

タイヤ接地荷重変動算出部１２２０は、真と判定したR_eを参照し、下記式（接地荷重変動に係る第一の式）

【0166】

【数16】

【0167】

によって、接地荷重変動dF_z0fl〜dF_z0rrを算出する。

【0168】

ここで、f(R_e)は、タイヤ有効半径R_eとの相関関係から決まる接地荷重F_z0の値であり、このようなR_eとF_z0との相関関係は、例えば、図６に示すマップで表される。このようなマップを参照することにより、タイヤ有効半径R_efl、R_efr、R_erl、R_errに応じた接地荷重F_z0が求められる。

【0169】

算出した接地荷重変動dF_z0fl〜dF_z0rrは、第１の入力行列演算部１２１１及び後述する路面変位算出部１２８０に供給される。また、算出した接地荷重変動dF_z0fl〜dF_z0rrは、各車輪の接地荷重定常値F_z0constfl〜F_z0constrrと加算されたうえで、各車輪の接地荷重F_z0fl〜F_z0rrとしてタイヤモデル演算部１２４０及び後述する参照用タイヤ有効半径演算部（参照用車輪有効半径演算部）１２７０に供給される。

【0170】

このように、タイヤ接地荷重変動算出部１２２０は、車輪角速度ω、バネ上前後速度u、ピッチレートq及び、ヨーレートr、及びスリップ比sを参照する。そして、車輪の接地荷重変動を計算し、フィードバック的に第１の入力行列演算部１２１１に入力する。

【0171】

なお、本実施形態において、上述した車体２００の前後方向のバネ上速度u、ピッチレートq、ヨーレートr、及び各車輪のスリップ比s_fl〜s_rrは、何れもセンサ値ではない。このように本実施形態において、ＥＣＵ６００の各部において算出した値を用いている。このため、追加のセンサを要することなく、状態量を好適に演算することができる。

【0172】

本実施形態において、参照用タイヤ有効半径演算部１２７０が算出する参照用タイヤ有効半径は、後述するようにスリップ算出部１２３０、タイヤモデル演算部１２４０等を介して結果的にタイヤ接地荷重変動算出部１２２０に参照される。このように、参照用タイヤ有効半径演算部１２７０が算出する参照用タイヤ有効半径は、タイヤ接地荷重変動算出部１２２０にフィードバック的に入力される。よって、状態量を精度よく演算することができる。

【0173】

（タイヤモデル演算部）
タイヤモデル演算部１２４０は、主演算部１２１０による演算結果の少なくとも一部を直接的または間接的に参照した非線形演算を行う。図５に示す例では、タイヤモデル演算部１２４０は、観測行列演算部１２１７による演算により得られる各車輪のスリップ比s_fl〜s_rr、各車輪のスリップ角β_fl〜β_rr、及びタイヤ接地荷重変動算出部１２２０が演算する接地荷重変動dF_z0fl〜dF_z0rrと定常値F_z0constfl〜F_z0constrrとを加算して得られる各車輪の接地荷重F_z0fl〜F_z0rrを参照して非線形演算を行う。すなわち、図５に示す例では、タイヤモデル演算部１２４０は、主演算部１２１０による演算結果の少なくとも一部を間接的に参照した非線形演算を行う。

【0174】

より具体的には、タイヤモデル演算部１２４０は、各車輪のスリップ比s_fl〜s_rr、各車輪のスリップ角β_fl〜β_rr、及び各車輪の接地荷重F_z0fl〜F_z0rrを参照する。そして、タイヤモデルに関する演算式を用いて、各車輪のタイヤ前後力F_x0fl〜F_x0rr、及び各車輪のタイヤ横力F_y0fl〜F_y0rrを算出する。タイヤモデル演算部１２４０による具体的な演算式は、本実施形態を限定するものではないが、例えば、タイヤモデル演算部１１４０が用いる数式と同様の数式を用いることができる。

【0175】

なお、主演算部１２１０による演算結果の少なくとも一部を直接的または間接的に参照した非線形演算を行うタイヤモデル演算部１２４０は、各車輪のタイヤ前後力F_x0fl〜F_x0rr、及び各車輪のタイヤ横力F_y0fl〜F_y0rrを算出するタイヤ力推定装置としても捉えることができる。

【0176】

（仮想バネ・ダンパ力算出部）
仮想バネ・ダンパ力算出部１２６０は、積分部１２１５により演算した結果の一部である各車輪のバネ下鉛直変位z_1flm〜z_1rrm、及び各車輪のバネ下鉛直速度w_1flm〜w_1rrm、を参照する。そして、各車輪に関する仮想バネ及び仮想ダンパによる力F_s1fl〜F_s1rrを以下の式により算出し、算出結果を第５の入力行列演算部１２１９に供給する。

【0177】

【数17】

【0178】

ここでk_s1f〜k_s1rは、各車輪のバネ下変位に乗ぜられる係数であり、c_s1f〜c_s1rは、各車輪のバネ下速度に乗ぜられる係数である。ここで、この仮想バネ・ダンパ力を考慮するために、バネ下鉛直運動に関する運動方程式としては、

【0179】

【数18】

【0180】

が用いられる。

【0181】

（参照用タイヤ有効半径演算部）
参照用タイヤ有効半径演算部１２７０は、入力された定常値加算後の接地荷重F_z0fl〜F_z0rrを参照する。そして、各車輪の参照用タイヤ有効半径R_efl、R_efr、R_erl、R_errを以下の式（参照用タイヤ有効半径に係る第二の式）によって演算する。さらに、演算した結果をスリップ算出部１２３０に供給する。

【0182】

【数19】

【0183】

ここで、B_Reff、D_Reff、F_Reff、F_z0nomは、実験結果等から得られるフィッティング係数である。また、k₁は、タイヤの上下変位に対する剛性であり、定数である。また、前述した接地荷重変動に係る第一の式中におけるf(R_e)マップは、上記参照用タイヤ有効半径に係る第二の式で表現されるF_z0とR_eとの関係と実質的に等価である。

【0184】

積算部１２３１は、スリップ算出部１２３０から供給される算出結果を積分し、積分後のスリップ比s_fl〜s_rrを算出する。積算部１２３１により算出された積分後のスリップ比s_fl〜s_rrは、タイヤ接地荷重変動算出部１２２０に供給される。

【0185】

参照用タイヤ有効半径R_efl〜R_errは、スリップ算出部１２３０に供給される。

【0186】

（路面変位算出部）
路面変位算出部１２８０は、主演算部１２１０が演算した各車輪のバネ下変位z_1flm〜z_1rrm、及び、タイヤ接地荷重変動算出部１２２０が算出した各車輪の接地荷重変動dF_z0fl〜dF_z0rrを用いて、

【0187】

【数20】

【0188】

により、各車輪の路面変位z_0fl〜z_0rrを算出する。算出した路面変位z_0fl〜z_0rrは、規範車両モデル演算部１１００に供給される。

【0189】

ここで、本実施形態に係る車両状態推定部１２００において、主演算部１２１０は線形演算を行い、タイヤモデル演算部１２４０は主演算部１２１０による演算結果の少なくとも一部を直接的または間接的に参照した非線形演算を行う。このように線形演算部と非線形演算部とを分離させた構成を採用することにより、車両モデルを用いた状態量の演算を好適に行うことができる。

【0190】

また、タイヤモデル演算部１２４０は、タイヤモデルに基づく非線形演算を行うので、非線形演算を線形演算から好適に分離することができる。

【0191】

また、上述のように、第３の入力行列演算部１２１３は、タイヤモデル演算部１２４０による非線形演算結果を入力として取り込むので、主演算部１２１０による線形演算に非線形演算結果を好適に取り込むことができる。したがって、主演算部１２１０は線形演算を行いつつ精度の高い演算を行うことができる。

【0192】

また、上述のように、本実施形態に係る車両状態推定部１２００の主演算部１２１０には、各車輪のタイヤ接地荷重変動算出部１２２０で算出した接地荷重変動dF_z0fl、dF_z0fr、dF_z0rl、dF_z0rrが入力される。このような構成とすることによって、追加センサを備える必要なく、状態量をより精度よく演算することができる。

【0193】

このように、本実施形態に係る車両状態推定部１２００は、標準的に装備された車載センサ情報のみを利用し、路面入力起因及び操舵入力起因により生じる車両挙動を区別せず同一の車両挙動推定部(車両モデル)を利用する、これにより、車両の状態を好適に推定することができる。特に、従来の技術では不可能であった、路面入力及び操舵入力が同時に生じる際の車両挙動の推定を実現することができる。

【0194】

以上の説明から明らかなように、本実施形態に係る車両状態推定部１２００は、車両における車輪の接地荷重F_z0の変動分である車輪の接地荷重変動dF_z0を算出するタイヤ接地荷重変動算出部１２２０と、１又は複数の入力値を参照してタイヤ有効半径R_eを算出する車輪有効半径算出部と、車両モデルを用いた演算を行う主演算部１２１０とを備えている。そして、タイヤ接地荷重変動算出部１２２０は、車輪有効半径算出部が算出したタイヤ有効半径R_eを入力とし接地荷重F_z0を出力とするマップを参照して、車輪の接地荷重変動dF_z0を算出する。また、主演算部１２１０は、タイヤ接地荷重変動算出部１２２０が算出した車輪の接地荷重変動dF_z0を含む１又は複数の入力値を参照し、車両状態に関する状態量についての演算を行うことによって１又は複数の出力値を算出する。よって、車両状態推定部１２００に入力される各車輪の接地荷重変動を好適に算出することができる。

【0195】

前述した実施形態では、タイヤ接地荷重変動算出部１２２０は、車輪有効半径算出部を兼ねているが、当該車輪有効半径算出部は、タイヤ接地荷重変動算出部１２２０から独立した構成であってもよい。たとえば、上記車輪有効半径算出部は、タイヤ接地荷重変動算出部１２２０から独立して配置され、u、q、r及びωが入力され、かつ算出したタイヤ有効半径R_eをタイヤ接地荷重変動算出部１２２０に供給するように構成されていてもよい。このような構成は、タイヤ接地荷重変動算出部１２２０の負荷を軽減する観点からより効果的である。

【0196】

また、上述のように、タイヤ接地荷重変動算出部１２２０は、車輪角速度、バネ上前後速度の推定値、ヨーレートの推定値、及びピッチレートの推定値を参照して、車輪の接地荷重変動を好適に算出することができる。

【0197】

また、上述のように、タイヤ接地荷重変動算出部１２２０は、車両の車輪速度の推定値V及び車輪角速度ωの一方または両方が０である場合に、定常値であるタイヤ有効半径R_econstを算出してもよい。この構成によれば、ホイルスピン（V_fl〜V_rrのいずれかがゼロ）の場合、あるいはホイルロック（ω_fl〜ω_rrのいずれかがゼロ）の場合においても接地荷重変動を参照する車両状態の推定制御を継続して実施することが可能である。よって、車輪の接地荷重変動dF_z0をより好適に算出することができる。

【0198】

上記の車輪速度Vの推定値及び車輪角速度ωの判定は、車輪速度V、車輪角速度ω、車輪速度Vの算出で参照されるバネ上速度u、ピッチレートq、ヨーレートr、からなる群から選ばれる一以上の数値が、算出されるタイヤ有効半径R_eとして許容される所定の範囲内に含まれているか否か、の判定に置き換えることが可能である。あるいは、上記の車輪速度Vの推定値及び車輪角速度ωの判定は、算出されたタイヤ有効半径R_eとして許容される所定の範囲内に含まれているか否か、の判定に置き換えることが可能である。このような許容される所定の範囲に基づく判定では、範囲外の場合に、前述したように定常値を採用してもよいし、あるいは範囲外の期間、接地荷重変動を参照する車両状態の制御を停止してもよい。

【0199】

また、上述のように、タイヤ接地荷重変動算出部１２２０は、スリップ比の推定値を更に参照して車輪の接地荷重変動を好適に算出することができる。さらに、車両状態推定部１２００は、参照用タイヤ有効半径R_eを算出する参照用タイヤ有効半径演算部１２７０を更に備え、スリップ算出部１２３０は、参照用タイヤ有効半径演算部１２７０が算出する参照用タイヤ有効半径R_eを、スリップ比の推定にさらに参照する。よって、タイヤ接地荷重変動算出部１２２０は、車輪の接地荷重変動をより一層好適に算出することができる。

【0200】

なお、参照用タイヤ有効半径演算部１２７０は、独自に参照用タイヤ有効半径R_eを算出するが、タイヤ接地荷重変動算出部１２２０で算出したタイヤ有効半径R_eをそのまま参照用タイヤ有効半径R_eとしてもよい。あるいは、参照用タイヤ有効半径演算部１２７０は、タイヤ接地荷重変動算出部１２２０で算出したタイヤ有効半径R_eと、加算部１２７５から供給される接地荷重F_z0とを参照して参照用タイヤ有効半径R_eを算出してもよい。

【0201】

また、前述したように、車両状態推定部１２００における１又は複数の入力値には、操作入力が含まれてもよい。このような構成によれば、車両の状態を好適に推定して高い乗り心地及び高い操縦安定性を実現する観点からより一層効果的である。

【0202】

また、本実施形態に係る車両状態推定部１２００は、車輪の接地荷重変動を少なくとも参照して、車輪の路面変位を算出する路面変位算出部を備えているので、車輪の路面変位を好適に算出することができる。

【0203】

また、上述のように、本実施形態に係る車両状態推定部１２００は、仮想バネ・ダンパ力算出部１２６０を備えており、仮想バネ・ダンパ力算出部１２６０は、仮想バネ及び仮想ダンパによる力（単に仮想的なバネの力ともいう）F_s1fl〜F_s1rrを算出し、算出結果を第５の入力行列演算部１２１９に供給する。

【0204】

発明者の知見によれば、仮想バネ及び仮想ダンパによる力を考慮しない場合、車両状態推定部１２００による状態量の推定において、推定結果が発散しやすいという傾向がある。上記のように仮想的なバネ・ダンパ力F_s1fl〜F_s1rrを算出し、算出結果を車両状態推定部１２００の主演算部１２１０に供給することによって、このような発散を抑制することができる。

【0205】

また、車両状態推定部１２００は、車両９００の前後方向のバネ上速度u、横方向のバネ上速度y、バネ上鉛直速度w、ロールレートp、ピッチレートq、ヨーレートr、及び、車輪のバネ下鉛直速度w_1flm〜w_1rrmを推定する。このような構成は、車両９００を偏差なく規範モデル特性に追従させるようにその車両状態を好適に制御する観点からより一層効果的である。

【0206】

また、車両状態推定部１２００は、車両９００のロール角、ピッチ角、ヨー角、各車輪のサスストローク変位、各車輪のバネ下鉛直変位、実舵角、実舵角速度、各車輪のスリップ比、各車輪のスリップ角、各車輪のタイヤ前後力、各車輪の横力、各車輪の接地荷重変動、各車輪のタイヤ有効半径、各車輪の参照用タイヤ有効半径、及び各車輪での路面変位の少なくとも何れかを推定する。このような構成も、車両９００を偏差なく規範モデル特性に追従させるようにその車両状態を好適に制御する観点からより一層効果的である。

【0207】

また、上記推定出力（出力値）は、バネ上鉛直速度、ロールレート、及びピッチレートを含む。このような構成も、車両９００を偏差なく規範モデル特性に追従させるようにその車両状態を好適に制御する観点からより一層効果的である。

【0208】

また、上記入力値は、更に、制御量演算部１０００の出力値を含み、当該出力値は、更に、車両９００の前後方向のバネ上速度u、車両の横方向のバネ上速度y、車両のヨーレートr、車両のロール角phi、車両のピッチ角theta、車両のヨー角psi、各車輪のサスストローク変位SusSt_fl〜SusSt_rr、各車輪のバネ下鉛直変位z_1flm〜z_1rrm、各車輪のバネ下鉛直速度w_1flm〜w_1rrm、実舵角δ、実舵角速度dδの少なくとも何れかを含んでいる。このような構成も、車両９００を偏差なく規範モデル特性に追従させるようにその車両状態を好適に制御する観点からより一層効果的である。

【0209】

また、ＥＣＵ６００は、前述した車両状態推定部１２００、規範車両モデル演算部１１００、減算部１０１２、積分部１０１４、第１の増幅部１０２１、第２の増幅部１０２２、第３の増幅部１０２３、ならびに、これらの増幅部による増幅結果を加算する加算部１０２４、を備える。そして、車両状態推定部１２００による推定結果を利用して車両状態の制御を行う。このため、高い乗り心地及び高い操縦安定性を実現することができる。

【0210】

また、本実施形態に係るサスペンション制御装置、及びサスペンション装置は、車両状態推定部１２００を備え、車両状態推定部１２００による推定結果を利用してサスペンションの制御を行う。このため、高い乗り心地及び高い操縦安定性を実現することができる。

【0211】

また、本実施形態に係るステアリング制御装置、及びステアリング装置は、車両状態推定部１２００を備え、車両状態推定部１２００による推定結果を利用してステアリングの制御を行う。このため、高い乗り心地及び高い操縦安定性を実現することができる。

【0212】

〔実施形態２〕
以下、本発明の実施形態２について、詳細に説明する。なお、既に説明した部材と同様の部材については同じ符号を付し、その説明を省略する。

【0213】

（車両９００の構成）
本実施形態に係る車両９００は、実施形態１と同様に、車輪速センサ３２０に加え、車両９００の横方向の加速度を検出する横Ｇセンサ３３０、車両９００の前後方向の加速度を検出する前後Ｇセンサ３４０、及び、車両９００のヨーレートを検出するヨーレートセンサ３５０を備えている。

【0214】

これらの各種センサによる検出結果は、ＥＣＵ６００ａに供給される。

【0215】

（ＥＣＵ６００ａ）
本実施形態に係る車両９００は、実施形態１において説明したＥＣＵ６００に代えて、ＥＣＵ６００ａを備えている。

【0216】

以下では、図７〜図９を参照して、本実施形態に係るＥＣＵ６００ａについて具体的に説明する。図７は、ＥＣＵ６００ａの概略構成例を示す図である。

【0217】

図７に示すように、ＥＣＵ６００ａは、ＥＣＵ６００と同様に、規範車両モデル演算部１１００を備えている。また、ＥＣＵ６００ａは、ＥＣＵ６００が備える車両状態推定部１２００に代えて、車両状態推定部１２００ａを備えている。

【0218】

また、ＥＣＵ６００ａには、車両各部１３００から、

【0219】

【数21】

【0220】

が入力される。ここで、ωは、実施形態１と同様に、車輪速センサ３２０が検出した車輪の角速度を表している。ドット付のu_sは、前後Ｇセンサ３４０が検出した車両のバネ上前後方向加速度を表している。ドット付のv_sは、横Ｇセンサ３３０が検出した車両のバネ上横方向加速度を表している。また、ＥＣＵ６００ａには実施形態１と同様に操作入力も入力される。なお、本実施形態において、下付き添え字「s」が付された状態量は、当該状態量がセンサによって検出されたものであることを示している。

【0221】

図７に示したＥＣＵ６００ａのその他の構成は、実施形態１と同様であるため、ここでは説明を省略する。

【0222】

（車両状態推定部）
図８を参照して本実施形態に係る車両状態推定部１２００ａについて説明する。図８は、車両状態推定部１２００ａの構成例を示すブロック図である。

【0223】

図８に示すように、本実施形態に係る車両状態推定部１２００ａは、実施形態１に係る車両状態推定部１２００が備える各構成に加え、フロント路面変位算出部１２８５、リア路面変位算出部１２９０、第１の入力行列演算部１２１１ａ、第２の入力行列演算部１２１１ｂ、推定状態量再構築部１２１０ａ、ゲイン適用部１２１１ｃ、減算部１２１１ｄ、及び第６の入力行列演算部１２１１ｅを備えている。なお、車両状態推定部１２００ａは、実施形態１に係る車両状態推定部１２００が備える路面変位算出部１２８０を備えていない。

【0224】

図８に示すように、車両状態推定部１２００ａには、実施形態１に係る車両状態推定部１２００への各種の入力に加え、
・前後Ｇセンサ３４０が検出した車両のバネ上前後方向加速度（ドット付のu_s）
・横Ｇセンサ３３０が検出した車両のバネ上横方向加速度（ドット付のv_s）
・ヨーレートセンサ３５０が検出した車両のヨーレートr_s
も入力される。また、推定状態量再構築部１２１０ａには、初期速度u_iniも入力される。

【0225】

（フロント路面変位算出部）
フロント路面変位算出部１２８５は、主演算部１２１０が演算した前輪のバネ下鉛直変位z_1flm、z_1frm、及び、タイヤ接地荷重変動算出部１２２０が算出した前輪の接地荷重変動dF_z0fl、dF_z0frを用いて、前輪の路面変位z_0fl、z_0frを算出する。フロント路面変位算出部１２８５による路面変位z_0fl、z_0frの具体的な算出処理は、路面変位算出部１２８０と同様であるのでここでは説明を省略する。

【0226】

（リア路面変位算出部）
リア路面変位算出部１２９０は、前輪の接地荷重変動、及び前輪のバネ下変位に基づき、後輪の路面変位を算出する。より具体的には、前輪の接地荷重変動dF_z0fl、dF_z0fr、及び前輪のバネ下鉛直変位z_1flm、z_1frmに基づき、以下の式を用いて、後輪の路面変位z_0rl、z_0rrを算出する。

【0227】

【数22】

【0228】

ここで、k_１は、タイヤ上下変位に対する剛性であり、定数である。

【0229】

また、リア路面変位算出部１２９０は、前輪の接地荷重変動dF_z0fl、dF_z0fr、及び前輪のバネ下鉛直変位z_1flm、z_1frmに基づき、算出した路面変位z_0rl、z_0rrを、前輪の路面変位に比べて、
t=WH/V_ave
の時間的遅れを伴わせて、第２の入力行列演算部１２１１ｂに入力する。ここで、WHは、車両９００のホイールベースの長さを表しており、V_aveは、例えば、車輪速センサ３２０による検出結果の四輪の平均値を利用できる。ここで、V_aveについては、推定のバネ上前後速度を利用してもよいし、ＧＰＳセンサの利用などの別手段で取得したものを利用してもよい。

【0230】

（第２の入力行列演算部）
リア路面変位に対する入力行列B0_r'に関する演算を行う第２の入力行列演算部１２１１ｂは、上述したリア路面変位算出部１２９０の出力に対して入力行列B0_r'を演算し、演算した結果を加算部１２１４に供給する。

【0231】

減算部１２１１ｄは、推定状態量再構築部１２１０ａの出力である、バネ上前後方向加速度の推定値、バネ上横方向加速度の推定値及びヨーレートの推定値から、前後Ｇセンサ３４０が検出した車両のバネ上前後方向加速度（ドット付のu_s）、横Ｇセンサ３３０が検出した車両のバネ上横方向加速度（ドット付のv_s）、及びヨーレートセンサ３５０が検出した車両のヨーレートr_sを減算して得られる結果を、ゲイン適用部１２１１ｃに供給する。

【0232】

（ゲイン適用部）
ゲイン適用部１２１１ｃは、減算部１２１１ｄの出力に対してオブザーバゲインを乗じることによって得られた結果を第６の入力行列演算部１２１１ｅに入力する。第６の入力行列演算部１２１１ｅは、ゲイン適用部１２１１ｃの出力に対して入力行列B6'を演算し、演算した結果を加算部１２１４に入力する。

【0233】

減算部１２１１ｄ及びゲイン適用部１２１１ｃが存在することによって、実施形態１において説明したバネ上並進及び回転運動に関する運動方程式のうち、F_x、F_y、M_zに関する運動方程式は、本実施形態に係る車両状態推定部１２００ａでは、以下のように変更される。

【0234】

【数23】

【0235】

ここで、ドット付のuは、車両のバネ上前後方向加速度の推定値を示しており、ドット付のu_sは、上述の通り、前後Ｇセンサ３４０が検出した車両のバネ上前後方向加速度である。ドット付のuと、ドット付のu_sとの差分は、上述した減算部１２１１ｄによって算出される。vについても同様である。また、rは、ヨーレートの推定値を示しており、r_sは、ヨーレートセンサ３５０が検出したヨーレートを示している。rとr_sとの差分は、上述した減算部１２１１ｄによって算出される。

【0236】

また、L₁、L₂、L₃は、ゲイン適用部１２１１ｃによって乗ぜられる各オブザーバゲインを表している。

【0237】

なお。L₁、L₂、L₃としては、それぞれ一定値を用いる構成としてもよいが、車体の速度に応じて変化する値を用いる構成としてもよい。

【0238】

（推定状態量再構築部）
図９は、推定状態量再構築部１２１０ａの構成例を示すブロック図である。推定状態量再構築部１２１０ａは、図９に示すように、加算部１２２１、及び推定状態量再計算部１２２２を備えている。推定状態量再構築部１２１０ａには、再構築前の状態量が入力される。再構築前の状態量のうち、車体２００の前後方向のバネ上速度uは、加算部１２２１にて、初期速度u_iniと加算されたうえで、推定状態量再計算部１２２２に入力される。また、推定状態量再計算部１２２２には以下が入力される。

【0239】

・車体２００の横方向のバネ上速度v
・車体２００の上下方向のバネ上速度w
・オイラー角の３成分（phi、theta、psi）
・前後Ｇセンサ３４０が検出した車両のバネ上前後方向加速度（ドット付のu_s）
・横Ｇセンサ３３０が検出した車両のバネ上横方向加速度（ドット付のv_s）
そして、推定状態量再計算部１２２２は、以下の式に基づき、再計算後の前後方向のバネ上速度u_new、再計算後の横方向のバネ上速度v_new、車両のバネ上前後方向加速度の推定値（ドット付のu）、及び車両のバネ上横方向加速度の推定値（ドット付のv）を算出する。

【0240】

【数24】

【0241】

ここで、下付き添え字「g」が付された状態量は、地上に固定されたグローバル座標系での状態量を示しており、下付き添え字「g」が付されていない状態量は、バネ上と同じ動きをするボディ座標系での状態量を示している。また、車体２００の上下方向のバネ上の加速度であるドット付のwとしては、推定された状態量を元に計算する。

【0242】

また、^GR_Bは、ボディ座標系からグローバル座標系への座標変換行列を示しており、また、^BR_Gは、グローバル座標系からボディ座標系への座標変換行列を示しており、それぞれ、以下の式で与えられる。

【0243】

【数25】

【0244】

推定状態量再計算部１２２２によって再計算された前後方向のバネ上速度u_new及び横方向のバネ上速度v_newは、図９に示すように、他の状態量と共に、再構築後の状態量として推定状態量再構築部１２１０ａから出力され、システム行列演算部１２１６及び減算部１２１１ｃに入力される。また、推定状態量再計算部１２２２によって計算された、ドット付きのｕ、ドット付きのｖは、状態量ｒとともに、減算部１２１１ｄに入力される。

【0245】

〔ソフトウェアによる実現例〕
ＥＣＵ６００、６００ａの制御ブロック（制御量演算部１０００、規範車両モデル演算部１１００、車両状態推定部１２００、１２００ａ）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ＣＰＵ（Central Processing Unit）を用いてソフトウェアによって実現してもよい。

【0246】

後者の場合、ＥＣＵ６００、６００ａは、各機能を実現するソフトウェアであるプログラムの命令を実行するＣＰＵ、上記プログラム及び各種データがコンピュータ（またはＣＰＵ）で読み取り可能に記録されたＲＯＭ（Read Only Memory）または記憶装置（これらを「記録媒体」と称する）、上記プログラムを展開するＲＡＭ（Random Access Memory）などを備えている。そして、コンピュータ（またはＣＰＵ）が上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。

【0247】

本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。

【符号の説明】

【0248】

１００ 懸架装置（サスペンション）
２００ 車体
６００、６００ａ ＥＣＵ（制御装置、サスペンション制御装置、サスペンション制御部、ステアリング制御装置）
１０００ 制御量演算部
１０１２ 減算部、
１０１４ 積分部
１０２１ 第１の増幅部
１０２２ 第２の増幅部
１０２３ 第３の増幅部
１０２４ 加算部
１１００ 規範車両モデル演算部
１１１０ 主演算部（演算部）
１１１１ 第１の入力行列演算部
１１１２ 第２の入力行列演算部
１１１３ 第３の入力行列演算部
１１１４ 加算部
１１１５ 積分部
１１１６ システム行列演算部
１１１７ 観測行列演算部
１１４０ タイヤモデル演算部
１２００、１２００ａ 車両状態推定部（車両状態推定装置）
１２１０ 主演算部（演算部、線形演算部）
１２１１ 第１の入力行列演算部（第１の演算部）
１２１２ 第２の入力行列演算部（第１の演算部）
１２１３ 第３の入力行列演算部（第１の演算部）
１２１４ 加算部
１２１５ 積分部
１２１６ システム行列演算部（第２の演算部）
１２１７ 観測行列演算部（第３の演算部）
１２４０ タイヤモデル演算部（非線形演算部）

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】