(57)【要約】【課題】車体横すべり角を考慮してヨー運動及び横運動を制御することで、車両応答性やドライバのフィーリング等を向上させる。【解決手段】車両挙動の制御方法は、車速を取得する工程と、ステアリングの操舵に対応する車輪の実舵角を取得する工程と、取得された車速及び実舵角に基づき、実舵角に対する車体横すべり角の伝達関数を規定する、（１）実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、（２）実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値に応じて、車両のヨー運動及び横運動を制御する工程と、を有する。【選択図】図４
2018171949.tif 000002
【特許請求の範囲】
【請求項１】
車両挙動の制御方法であって、
車速を取得する工程と、
ステアリングの操舵に対応する車輪の実舵角を取得する工程と、
取得された上記車速及び上記実舵角に基づき、実舵角に対する車体横すべり角の伝達関数を規定する、（１）実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、（２）実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値に応じて、車両のヨー運動及び横運動を制御する工程と、
を有する、ことを特徴とする車両挙動の制御方法。
【請求項２】
上記制御する工程では、上記第１パラメータと上記第２パラメータとを積算した値を用いて、上記車速及び上記実舵角に応じて車体横すべり角を調整することで、上記ヨー運動と上記横運動との関係を制御する、請求項１に記載の車両挙動の制御方法。
【請求項３】
上記制御する工程では、上記ヨー運動と上記横運動との相対的な差が小さくなるように、上記第１パラメータと上記第２パラメータとを積算した値を用いて上記車体横すべり角の過渡特性を調整する、請求項２に記載の車両挙動の制御方法。
【請求項４】
上記制御する工程では、上記第１パラメータと上記第２パラメータとを積算した値が小さくなるように上記ヨー運動及び上記横運動を制御する、請求項１乃至３のいずれか一項に記載の車両挙動の制御方法。
【請求項５】
上記第１パラメータと上記第２パラメータとを積算した値は、正の値の範囲内において車速に応じて大きくなる、請求項１乃至４のいずれか一項に記載の車両挙動の制御方法。
【請求項６】
上記第１パラメータは次式により表される、請求項１乃至５のいずれか一項に記載の車両挙動の制御方法。
2018171949.tif 000017
ここで、上式において、
2018171949.tif 000018
2018171949.tif 000019
「ｍ」は車両の質量であり、「Ｖ」は車速であり、「ｌ」はホイールベースであり、「ｌf」は車両重心点と前輪の車軸との距離であり、「ｌr」は車両重心点と後輪の車軸との距離であり、「Ｋf」は前輪１輪あたりのタイヤコーナリングパワー、「Ｋr」は後輪１輪あたりのタイヤコーナリングパワーであり、「δ」は車輪の実舵角であり、「β」は車体スリップ角である。
【請求項７】
上記第２パラメータは次式により表される、請求項１乃至６のいずれか一項に記載の車両挙動の制御方法。
2018171949.tif 000020
ここで、上式において、「Ｔβ」は第２パラメータであり、「ｍ」は車両の質量であり、「Ｖ」は車速であり、「Ｉ」は車両のヨー慣性モーメントであり、「ｌ」はホイールベースであり、「ｌf」は車両重心点と前輪の車軸との距離であり、「ｌr」は車両重心点と後輪の車軸との距離であり、「Ｋf」は前輪１輪あたりのタイヤコーナリングパワー、「Ｋr」は後輪１輪あたりのタイヤコーナリングパワーであり、「β」は車体スリップ角である。
【請求項８】
ドライビングシミュレータを用いた車両挙動のシミュレーション方法であって、
上記ドライビングシミュレータに適用する車速を取得する工程と、
上記ドライビングシミュレータに適用する車輪の実舵角を取得する工程と、
取得された上記車速及び上記実舵角に基づき、実舵角に対する車体横すべり角の伝達関数を規定する、（１）実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、（２）実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値に応じて、車両のヨー運動及び横運動を模擬した運動を上記ドライビングシミュレータにより被験者に与える工程と、
を有する、ことを特徴とする車両挙動のシミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、車両挙動の制御方法及び車両挙動のシミュレーション方法に関する。
【背景技術】
【０００２】
本発明に関連する技術が、例えば特許文献１及び２などに開示されている。特許文献１には、車輪の横力及びタイヤ横すべり角を利用して、横方向の車両運動を正確に記述するようにした車両挙動再現システムが開示されている。特許文献２には、ドライバはヨー運動と横運動とのバランスを感じながらステアリング操作を行っているため、このバランスが変化すると違和感を覚えるという知見に基づき、運転者のステアリング操作に対して所定の応答遅れを適用して、後輪を転舵すると共に旋回外輪に制動トルクを発生させるようにした旋回挙動制御装置が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第４１４０７２０号公報
【特許文献２】特開２００７−１７６３１４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、ドライバは、操舵に起因するヨー運動及び横運動の大きさの変化や前後関係を前庭感覚や視覚などで認識し、車両のコーナリング時における応答性の良し悪しを感じていると考えられる。しかしながら、従来技術では、車両のコーナリング時におけるドライバにとっての車両応答性の評価を、ヨー運動及び横運動を用いて説明することが困難であった、つまり定量的に評価できなかった。
【０００５】
ここで、或る車速Ｖにおいて、横加速度とヨーレートｒと車体横すべり角βは、次のような関係を有している。なお、下式では、横加速度を、横変位ｙを２回微分したものとして、ｙの上に２つのドットを付して表記している。
2018171949.tif 000003
【０００６】
上記の関係式から、車体横すべり角βが、横運動とヨー運動との間に介在して両者の応答の違いに影響を与えているものと考えられる。車両のスピンやドリフト等を除けば、車体横すべり角βは微小量であり、「ｃｏｓβ」は１の概数となる。そのため、上式によれば、車体横すべり角βの時間変化が、主として横運動とヨー運動との相対的な差を招いているものと考えられる。また、車体横すべり角βの時間変化によって、高速になるほど、操舵に対するヨーレートの応答遅れに対して横加速度の応答遅れが相対的に大きくなるという知見もある。
【０００７】
このようなことから、車体横すべり角の特性（特に過渡特性）が変われば、ヨー運動と横運動との関係が変わり、こうしてヨー運動と横運動との関係が変わることで、車両応答性に関するドライバのフィーリングが変化するものと考えられる。従来技術においては、このような車体横すべり角の過渡特性に着眼してヨー運動と横運動に関する制御を行う試みはなされていなかった。
【０００８】
本発明は、上述した従来技術の問題点を解決するためになされたものであり、車体横すべり角を考慮してヨー運動及び横運動を制御することで、車両応答性やドライバのフィーリング等を向上させることができる車両挙動の制御方法及び車両挙動のシミュレーション方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記の目的を達成するために、本発明は、車両挙動の制御方法であって、車速を取得する工程と、ステアリングの操舵に対応する車輪の実舵角を取得する工程と、取得された車速及び実舵角に基づき、実舵角に対する車体横すべり角の伝達関数を規定する、（１）実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、（２）実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値に応じて、車両のヨー運動及び横運動を制御する工程と、を有する、ことを特徴とする。
このように構成された本発明によれば、実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値を用いることで、ヨー運動及び横運動を適切に制御できるようになる。すなわち、車速及び実舵角に適した所望のヨー運動及び横運動を実現できるようになる。これにより、本発明によれば、操舵により旋回するときの車両挙動に対するフィーリングを向上させることができる、つまりドライバが操舵したときに感じる操縦安定性（車両応答性など）を向上させることができる。
【００１０】
本発明において、好ましくは、制御する工程では、第１パラメータと第２パラメータとを積算した値を用いて、車速及び実舵角に応じて車体横すべり角を調整することで、ヨー運動と横運動との関係を制御する。
このように構成された本発明によれば、第１パラメータと第２パラメータとを積算した値を用いて、車速及び実舵角に応じて車体横すべり角を調整することで、ヨー運動及び横運動をより適切に制御できるようになる。
【００１１】
本発明において、好ましくは、制御する工程では、ヨー運動と横運動との相対的な差が小さくなるように、第１パラメータと第２パラメータとを積算した値を用いて車体横すべり角の過渡特性を調整する。
このように構成された本発明によれば、ヨー運動と横運動との相対的な差が小さくなるので、操舵時（詳しくは舵角が変化しているとき）において車両の向きに対する車両が進んでいる方向の変化を小さくすることができる。したがって、本発明によれば、車両の操縦安定性が向上して、旋回時の修正操舵を低減することができる。その結果、操縦安定性が更に向上して、ドライビングフィーリングをより一層向上させることができる。
【００１２】
本発明において、好ましくは、制御する工程では、第１パラメータと第２パラメータとを積算した値が小さくなるようにヨー運動及び横運動を制御する。
このように構成された本発明によっても、ヨー運動と横運動との相対的な差が小さくなるので、車両の操縦安定性が向上して、旋回時の修正操舵を低減することができる。
【００１３】
本発明において、好ましくは、第１パラメータと第２パラメータとを積算した値は、正の値の範囲内において車速に応じて大きくなる。
【００１４】
本発明において、好ましくは、第１パラメータは次式により表される。
2018171949.tif 000004
ここで、上式において、
2018171949.tif 000005
2018171949.tif 000006
「ｍ」は車両の質量であり、「Ｖ」は車速であり、「ｌ」はホイールベースであり、「ｌf」は車両重心点と前輪の車軸との距離であり、「ｌr」は車両重心点と後輪の車軸との距離であり、「Ｋf」は前輪１輪あたりのタイヤコーナリングパワー、「Ｋr」は後輪１輪あたりのタイヤコーナリングパワーであり、「δ」は車輪の実舵角であり、「β」は車体スリップ角である。
【００１５】
本発明において、好ましくは、第２パラメータは次式により表される。
2018171949.tif 000007
ここで、上式において、「Ｔβ」は第２パラメータであり、「ｍ」は車両の質量であり、「Ｖ」は車速であり、「Ｉ」は車両のヨー慣性モーメントであり、「ｌ」はホイールベースであり、「ｌf」は車両重心点と前輪の車軸との距離であり、「ｌr」は車両重心点と後輪の車軸との距離であり、「Ｋf」は前輪１輪あたりのタイヤコーナリングパワー、「Ｋr」は後輪１輪あたりのタイヤコーナリングパワーであり、「β」は車体スリップ角である。
【００１６】
他の観点では、上記の目的を達成するために、本発明は、ドライビングシミュレータを用いた車両挙動のシミュレーション方法であって、ドライビングシミュレータに適用する車速を取得する工程と、ドライビングシミュレータに適用する車輪の実舵角を取得する工程と、取得された車速及び実舵角に基づき、実舵角に対する車体横すべり角の伝達関数を規定する、（１）実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、（２）実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値に応じて、車両のヨー運動及び横運動を模擬した運動をドライビングシミュレータにより被験者に与える工程と、を有する、ことを特徴とする。
このように構成された本発明によれば、実舵角により生じる車体横すべり角のゲインを示す第１パラメータと、実舵角の変化に起因する車体横すべり角の過渡的変化を示す第２パラメータと、を積算した値を可変パラメータとして適用することで、車体横すべり角の過渡特性を変更したシミュレーションを適切に行うことができる。
【発明の効果】
【００１７】
本発明の車両挙動の制御方法及び車両挙動のシミュレーション方法によれば、車体横すべり角を考慮してヨー運動及び横運動を制御することで、車両応答性やドライバのフィーリング等を向上させることができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の実施形態によるドライビングシミュレータシステムの概略構成図である。
【図２】車速に対するＧβ及びＴβのそれぞれの変化を示す図である。
【図３】車速に対するＧβ・Ｔβの変化を示す図である。
【図４】本発明の実施形態において適用する水平運動の概念モデルを示すブロック図である。
【図５】本発明の実施形態によるドライビングシミュレータシステムのコントローラが実行する処理を示すフローチャートである。
【図６】本実験において適用した各種応答パラメータのステップ応答を示す。
【図７】本実験において適用した実験モードについての説明図である。
【図８】レーンチェンジ中の実舵角と横変位との関係を人の伝達関数として表したブロック図である。
【図９】本実験によりドライバから仕様１〜３について得られたドライバパラメータτLの算出結果を示す図である。
【図１０】本実験によりドライバから仕様１〜３について得られたドライバパラメータｈの算出結果を示す図である。
【図１１】本実験によりドライバから仕様１〜３について得られたドライバパラメータτhの算出結果を示す図である。
【図１２】官能評価での仕様１〜３に対するドライバのコメントを示す図である。
【図１３】官能評価においてドライバが乗りやすいと評価した仕様の順番を示す図である。
【図１４】本発明による車両挙動の制御方法を適用した車両の概略構成図である。
【発明を実施するための形態】
【００１９】
以下、添付図面を参照して、本発明の実施形態による車両挙動の制御方法及び車両挙動のシミュレーション方法について説明する。
【００２０】
１．ドライビングシミュレータシステム
まず、図１を参照して、本発明の実施形態によるドライビングシミュレータシステムについて説明する。図１は、本発明の実施形態による車両挙動の制御方法及び車両挙動のシミュレーション方法を適用したドライビングシミュレータシステムの概略構成図である。
【００２１】
図１に示すように、ドライビングシミュレータシステム１は、主に、車両挙動を模擬した運動を被験者に与えるよう構成されたドライビングシミュレータ２と、ドライビングシミュレータ２を制御するコントローラ１０と、を有する。具体的には、ドライビングシミュレータ２は、被験者によって操作されるステアリング４と、被験者が着座するシート６と、被験者に提示すべき視覚刺激が表示される表示部（スクリーン）８と、を備える。
【００２２】
また、ドライビングシミュレータ２では、シート６が、ヨー運動用アクチュエータ１２によって駆動されることで、ターンテーブル（不図示）などを介して矢印Ａ１で示すように動くようになっている。これにより、疑似的な車両のヨー運動を被験者に与えるようになっている。加えて、シート６が、横運動用アクチュエータ１４によって駆動されることで、リニアレール（不図示）などを介して矢印Ａ２で示すように動くようになっている。これにより、疑似的な車両の横運動を被験者に与えるようになっている。
【００２３】
これらのヨー運動用アクチュエータ１２及び横運動用アクチュエータ１４は、コントローラ１０によって制御される。具体的には、コントローラ１０は、少なくともステアリング４の操作に対応する操舵角や車速（例えば被験者による図示しないアクセルペダルの操作によって設定される）などに応じた疑似的なヨー運動及び横運動を被験者に与えるように、ヨー運動用アクチュエータ１２及び横運動用アクチュエータ１４を制御する。また、コントローラ１０は、表示部８に対する表示制御も行う。この場合、コントローラ１０は、上記のヨー運動及び横運動に応じて変化する画像を表示部８に表示させるようにする。
【００２４】
なお、上記したヨー運動は、ヨー角度やヨーレートなどを含む運動であり、上記した横運動は、横変位や横加速度などを含む運動である。
【００２５】
２．シミュレーション方法
次に、上記したドライビングシミュレータシステム１を利用した、本実施形態によるシミュレーション方法について説明する。
【００２６】
２−１．伝達関数
最初に、本実施形態によるシミュレーション方法において適用する車両運動モデルの伝達関数について説明する。
【００２７】
車両の重心点を原点とし車両に固定した座標系において、水平面内の車両運動を記述する一般的な方程式をラプラス変換すると、以下の式（１）及び式（２）が得られる。式（１）は、水平面内の車両運動において横運動に関するつり合いの式をラプラス変換したものであり、式（２）は、水平面内の車両運動においてヨー運動に関するつり合いの式をラプラス変換したものである。
【００２８】
2018171949.tif 000008
2018171949.tif 000009
【００２９】
なお、本明細書では、各種パラメータを指し示す記号を以下のように定義するものとする。
ｍ：車両の質量
Ｉ：車両のヨー慣性モーメント
ｌ：ホイールベース
ｌf：車両重心点と前輪の車軸との距離
ｌr：車両重心点と後輪の車軸との距離
Ｋf：前輪１輪あたりのタイヤコーナリングパワー
Ｋr：後輪１輪あたりのタイヤコーナリングパワー
δ：前輪の実舵角
ｙ：横変位
ｒ：ヨーレート
β：車体横すべり角
ｓ：ラプラス演算子
【００３０】
式（１）及び式（２）において、「β（ｓ）」は車体横すべり角βのラプラス変換を示し、「ｒ（ｓ）」はヨーレートｒのラプラス変換を示し、「δ（ｓ）」は前輪の実舵角δのラプラス変換を示している。
【００３１】
さらに、式（１）及び式（２）に示すβ（ｓ）とｒ（ｓ）に関する代数方程式（連立方程式）を機械的に解くと、以下の式（３）及び式（４）が得られる。
【００３２】
2018171949.tif 000010
2018171949.tif 000011
【００３３】
式（３）は、実舵角に対する車体横すべり角の伝達関数を示し、式（４）は、実舵角に対するヨーレートの伝達関数を示している。式（３）及び式（４）中の各種パラメータの定義は以下の通りである。
【００３４】
2018171949.tif 000012
2018171949.tif 000013
【００３５】
２−２．車速に応じた車体横すべり角の特性
上述した式（３）より、操舵に対する車体横すべり角の特性は、Ｇβ、Ｔβ、ωn、ζの応答パラメータに応じたものとなる。これら応答パラメータの中で車体横すべり角に固有のパラメータは、Ｇβ及びＴβとなる。以下では、Ｇβ及びＴβにおける車速に応じた変化について説明する。
【００３６】
なお、Ｇβは、実舵角に対する車体横すべり角（実舵角により生じる横すべり角）のゲイン（ほぼ定常値）に相当する、換言すると実舵角の変化が車体横すべり角に及ぼす影響の大きさを示す係数に相当する。また、Ｔβは、実舵角の変化に起因する車体横すべり角の過渡変化に相当する。Ｇβは、本発明における「第１パラメータ」に相当し、Ｔβは、本発明における「第２パラメータ」に相当する。
【００３７】
図２は、或る車両諸元を仮定したときの、車速に応じたＧβ及びＴβを計算した結果を示す。図２は、横軸に車速Ｖ［ｋｍ／ｈ］を示し、縦軸にＧβ及びＴβを示している。具体的には、実線で表したグラフＧ１１はＧβを示し、破線で表したグラフＧ１２はＴβを示している。
【００３８】
図２に示すように、Ｇβは、車速が高くなるにつれて減少する。具体的には、Ｇβは、所定車速未満では正の値であるが、所定車速以上では負の値となる。つまり、Ｇβは、所定車速において正負が入れ替わる。基本的には、Ｇβは、所定車速未満の領域では（つまりＧβが正である場合）、コーナリング中などにおいて車輪の方向（車体横すべり角に相当）が車両進行方向に対して外側を向き、所定車速以上の領域では（つまりＧβが負である場合）、コーナリング中などにおいて車輪の方向（車体横すべり角に相当）が車両進行方向に対して内側を向くという傾向を示すものとなる。
他方で、Ｔβは、所定車速未満の範囲においては、正の値となり、車速が高くなるにつれて増加する。一方で、所定車速以上の範囲においては、Ｔβは、負の値となり（つまり所定車速を境に正負が入れ替わる）、且つ車速が高くなるにつれて負の値の範囲内において増加する。例えば、Ｔβは、実舵角が変化しているときは、車輪の方向（車体横すべり角に相当）が少し外側に向き、実舵角が変化しなくなると車輪の方向（車体横すべり角に相当）が一定になるという傾向を示すものとなる。
【００３９】
このようにＧβにおいて正負が変わる所定車速とＴβにおいて正負が変わる所定車速とは、等しい車速となり、この車速は以下の式（１２）で表される。式（１２）は、上記した式（５）及び式（６）において「Ｖ2」が適用された部分から得られる。
【００４０】
2018171949.tif 000014
【００４１】
式（１２）に示す車速において、Ｔβは不連続点を有する。つまり、Ｔβは、当該車速において正負が変わり、値も大きく変化する。そのため、Ｔβ単体では、制御上において扱いづらいパラメータであると言える。一方で、Ｇβは、車速に比例して変化するものなので、Ｔβよりも扱いやすいパラメータであると言える。ここで、ＧβとＴβとの関係に着目すると、同じ車速で見たときに、Ｇβが減少するとＴβは増加し、Ｇβが増加すればＴβは減少するという関係を有する。したがって、ＧβとＴβの積（Ｇβ・Ｔβ）は、必ず正となる。なお、ＧβとＴβの積を取る演算は、式（３）に示した実舵角に対する車体横すべり角の伝達関数において実際に行われるものである。
【００４２】
図３は、車速に対するＧβ・Ｔβの変化を示す。図３は、横軸に車速Ｖ［ｋｍ／ｈ］を示し、縦軸にＧβ・Ｔβを示している。図３に示すように、Ｇβ・Ｔβは車速が高くなるにつれて増加し（但し、車速が高くなると、Ｇβ・Ｔβの増加度合い（変化率）が小さくなる）、また、Ｇβ・Ｔβは不連続点を有しないことがわかる。そのため、制御上において、Ｔβ単体ではコントロールしづらいが、Ｇβ・Ｔβはコントロールしやすいと言える。したがって、本実施形態では、シミュレーションにおいて制御に用いるパラメータとして、ＧβとＴβの積を採用することとした。
【００４３】
ここで、Ｇβ・Ｔβは、車体横すべり角の伝達関数（式（３）参照）において微分要素係数である。このＧβ・Ｔβは、実舵角の変化に起因する車体横すべり角を指し示すものとなる。よって、Ｇβ・Ｔβは、車速に応じて増加することから、高速になるほど、車体横すべり角の過渡的変化が大きくなり、ヨー運動と横運動の相対的な差が広がるものと考えられる。
【００４４】
２−３．車両運動モデル
次に、図４を参照して、本実施形態によるシミュレーション方法において適用する、具体的な車両運動モデルについて説明する。図４は、本実施形態において適用する水平運動の概念モデルを示す。このモデルは、上記した応答パラメータにより、ヨーと車体横すべり角特性を独立して調整するようにしたモデルである。
【００４５】
図４に示すモデルには、前輪の実舵角δと車速Ｖが入力される。次に、この実舵角δから、式（３）に基づき、車体横すべり角βが算出される。この場合、本実施形態では、式（３）中のＧβ・Ｔβを、図３に示したような関係を考慮して設定するようにする。また、このような車体横すべり角βの算出と並行して、実舵角δから、式（４）に基づき、ヨーレートｒが算出される。そして、車体横すべり角βと、当該車体横すべり角βを微分した値と、ヨーレートｒと、車速Ｖとから、横加速度が算出され、この横加速度を積分（２回積分）することで横変位ｙが算出される。さらに、ヨーレートｒを積分（１回積分）することでヨー角度が算出される。
【００４６】
本実施形態では、このようなモデルより算出された車体横すべり角β、横加速度、横変位ｙ、ヨー角度及びヨーレートｒのパラメータに基づき、コントローラ１０がヨー運動用アクチュエータ１２及び横運動用アクチュエータ１４を制御する。具体的には、コントローラ１０は、ヨー角度及びヨーレートｒに基づきヨー運動用アクチュエータ１２を制御すると共に、横変位ｙ及び横加速度に基づき横運動用アクチュエータ１４を制御することで、これらのパラメータに応じたヨー運動及び横運動を模擬した運動を被験者に与えるようにする。加えて、コントローラ１０は、このヨー運動及び横運動に応じた画像を表示させるように表示部８を制御する。
【００４７】
ここで、従来の一般的な車両運動モデルでは、モデルを積分器型に表しているが、そのような積分器型に表したモデルでは、Ｇβ・Ｔβを独立して制御することはできない。これに対して、本実施形態では、図４に示したように、モデルを伝達関数型に表している。こうすることで、Ｇβ・Ｔβを独立して制御することができるのである。
【００４８】
２−４．処理フロー
次に、図５を参照して、本実施形態によるシミュレーション方法を実現すべく、ドライビングシミュレータシステム１のコントローラ１０が行う処理フローについて説明する。
【００４９】
まず、ステップＳ１では、コントローラ１０は、事前に設定された各種パラメータをメモリなどから取得する。例えば、コントローラ１０は、車両の質量ｍ、車両のヨー慣性モーメントＩ、ホイールベースｌ、車両重心点と前輪の車軸との距離ｌf、車両重心点と後輪の車軸との距離ｌr、前輪１輪あたりのタイヤコーナリングパワーＫf、後輪１輪あたりのタイヤコーナリングパワーＫrなどを取得する。そして、コントローラ１０は、車速に応じて、Ｇr、Ｔr、ωn、ζ、Ｇβを算出する。また、コントローラ１０は、車速に応じて事前に定められたＧβ・Ｔβ（つまり車速とＧβ・Ｔβとの関係）も取得する。
【００５０】
次いで、ステップＳ２では、コントローラ１０は、ステアリング４の操舵角に対応する前輪の実舵角δとアクセルペダルの操作に対応する車速Ｖ（１つの例では車速Ｖとして一定値が適用される）とを取得する。
【００５１】
次いで、ステップＳ３では、コントローラ１０は、ステップＳ１で取得した各種パラメータ及びステップＳ２で取得した実舵角δ並びに車速Ｖに基づき、図４に示した車両運動モデルを用いて、車体横すべり角β、横加速度、横変位ｙ、ヨー角度及びヨーレートｒを算出する。この場合、コントローラ１０は、ステップＳ１で取得した車速とＧβ・Ｔβとの関係を参照して、ステップＳ２で取得した車速Ｖに対応するＧβ・Ｔβを得て、このＧβ・Ｔβを適用して車体横すべり角βを算出する。
【００５２】
次いで、ステップＳ４では、コントローラ１０は、ステップＳ３で算出した車体横すべり角β、横加速度、横変位ｙ、ヨー角度及びヨーレートｒのパラメータに基づき、ヨー運動用アクチュエータ１２及び横運動用アクチュエータ１４を制御する。具体的には、コントローラ１０は、ヨー角度及びヨーレートｒに基づきヨー運動用アクチュエータ１２を制御すると共に、横変位ｙ及び横加速度に基づき横運動用アクチュエータ１４を制御する。
なお、このステップＳ４の処理が終了したら、コントローラ１０は、ステップＳ２に戻り、ステップＳ２以降の処理を再度行う。この場合、コントローラ１０は、最初のステップＳ１の処理において取得した各種パラメータの値を用いて、ステップＳ２以降の処理を行う。
但し、他の例では、ステップＳ４の終了後、コントローラ１０は、ステップＳ１に戻り、ステップＳ１以降の処理を再度行ってもよい。この例では、コントローラ１０は、ステップＳ１において、前輪１輪あたりのタイヤコーナリングパワーＫf及び後輪１輪あたりのタイヤコーナリングパワーＫrのそれぞれの現在の値を再び取得して、これらのタイヤコーナリングパワーＫf、Ｋrの値を用いてステップＳ２以降の処理を行う。こうすることで、より精度の高いシミュレーションを実現することができる。
【００５３】
３．実験内容
次に、本実施形態によるシミュレーション方法をドライビングシミュレータシステム１に適用して行った、具体的な実験の内容について説明する。
【００５４】
３−１．実験仕様
図６は、本実験において適用した応答パラメータに関する応答を示す。具体的には、図６は、本実験において適用した前輪の実舵角δ、ヨーレートｒ及び車体横すべり角βのそれぞれの時間変化を示している。
【００５５】
上述したように、本実施形態ではＧβ・Ｔβに着眼している。このＧβ、Ｔβ以外の応答パラメータは、或る車両諸元を仮定して、車両速度を６０［ｋｍ／ｈ］に固定して算出した（式（７）〜式（１１）参照）。図６は、こうして算出した応答パラメータを車両運動モデル（図４参照）に適用した場合の、実舵角δ、ヨーレートｒ及び車体横すべり角βのそれぞれのステップ応答を示している。本実験では、このようなステップ応答を有する実舵角δ、ヨーレートｒ及び車体横すべり角βを用いた。
【００５６】
特に、車体横すべり角βについては、グラフＧ２１、Ｇ２２、Ｇ２３に示すような３つの値を用いて実験を行った。具体的には、Ｔβを調整して３つのＧβ・Ｔβを用意し、これらＧβ・Ｔβのそれぞれより得られた３つの車体横すべり角βを用いた。グラフＧ２１の車体横すべり角βに適用したＧβ・Ｔβの値を第１の値とし、グラフＧ２２の車体横すべり角βに適用したＧβ・Ｔβの値を第２の値とし、グラフＧ２３の車体横すべり角βに適用したＧβ・Ｔβの値を第３の値とすると、これらは、第２の値を基準とすると、「第１の値＝第２の値／２」及び「第３の値＝第２の値×２」という関係を有する。グラフＧ２１、Ｇ２２、Ｇ２３に示すように、Ｇβ・Ｔβの値が小さくなると、車体横すべり角βの時間変化（応答）の傾きが若干緩やかになる。以下では、第１の値をＧβ・Ｔβに適用して行った実験を「仕様１」と呼び、第２の値をＧβ・Ｔβに適用して行った実験を「仕様２」と呼び、第３の値をＧβ・Ｔβに適用して行った実験を「仕様３」と呼ぶ。
【００５７】
３−２．実験モード
図７は、本実験で適用する実験モードについての説明図である。図７において、符号ＰＹ１、Ｃ１は、車両のレーンチェンジ前の道路を規定するパイロン及び中央線をそれぞれ示し、符号ＰＹ２、Ｃ２は、車両のレーンチェンジ後の道路を規定するパイロン及び中央線をそれぞれ示している。これら道路の幅は３．０［ｍ］であり、これら道路の中央線Ｃ１、Ｃ２の間隔は２．５［ｍ］であり、これら道路間の距離、つまりレーンチェンジ区間は３５［ｍ］であるものとする。
【００５８】
図７に示すように、本実験は、車両をレーンチェンジさせるような状況において行われた。具体的には、ドライビングシミュレータシステム１の表示部８にパイロンＰＹ１、ＰＹ２や中央線Ｃ１、Ｃ２などに対応する画像が表示され、被験者は、表示部８に表示された画像を見ながら、ステアリング４の操作のみで車両のレーンチェンジタスクを行った。この場合、被験者がアクセルペダルを踏むと、車速を一定速度６０［ｋｍ／ｈ］に設定した。より詳しくは、被験者は、車両が直進後にレーンチェンジ区間に入った場所において、ダッシュボードに設置されたランプが点灯した場合にのみ、レーンチェンジを行うこととした。
【００５９】
また、実験は、被験者５人（ドライバＡ〜Ｅ）に対して行われた。具体的には、一人の被験者につき、仕様１〜３のそれぞれを９回実験した。更に、被験者が実験仕様をわからないようにするために、仕様をランダムに変更して実験した。
【００６０】
３−３．評価方法
従来の研究において、計測されたドライバの操舵行動と車両運動から、ドライバモデルのパラメータを算出し、運動性能の評価手法が提案されている。特に、当該研究では、そのような評価手法を応答パラメータωn、ζの検討に適用し、官能評価と相関するという結果が得られている。したがって、本実験でも同様の手法を利用した。
【００６１】
レーンチェンジ中の実舵角と横変位との関係は、以下の式（１３）のように人の伝達関数として表せる。また、この式（１３）を図示化すると、図８のようになる。図８は、ドライバと車両のシステムに関する閉ループのモデルを示す。
【００６２】
2018171949.tif 000015
【００６３】
式（１３）及び図８において、「ｙOL」は目標レーンチェンジ幅であり、「ｈ」、「τh」、「τL」はそれぞれドライバパラメータである。これらのドライバパラメータｈ、τh、τLは、実験において計測した車両軌跡及び実舵角の時刻歴データから算出される。具体的には、以下の式（１４）で表される評価関数Ｊを最小にするように、ドライバパラメータｈ、τh、τLを求めている。式（１４）中のσ*及びｙ*は、それぞれ、計測された実舵角及び横変位である。
【００６４】
2018171949.tif 000016
【００６５】
こうして算出されたドライバパラメータを用いて、ドライバ（被験者）の操舵行動を評価する。この場合、ｈはゲインに相当し、τhは予見時間に相当し、τLは全ての遅れを一次遅れで代表した時定数に相当する。この中でもτLは、性能評価に相関する量として考えられている。従来の研究では、τLが大きいほど、官能評価が良くなるという傾向が観測されている。本実験では、仕様１〜３のそれぞれにおける５試行のデータから、算出精度の良い試行を選んで、結果として用いた。
【００６６】
３−４．実験結果
図９、図１０、図１１は、上述した実験により、ドライバ（被験者）Ａ〜Ｅから仕様１〜３について得られたドライバパラメータτL、ｈ、τhの算出結果をそれぞれ示している。
【００６７】
他方で、上記の計測実験と同様のレーンチェンジと２レーンを用いて、ドライバに自由に走行を行わせて、官能評価を実行した。車速は６０［ｋｍ／ｈ］の一定に制御され、ドライバはステアリング操作のみで評価を行った。計測実験と同様に、仕様を伝えずに評価を行った。図１２は、仕様１〜３に対するドライバ（被験者）のコメントを示し、図１３は、ドライバが乗りやすいと評価した仕様の順番を示している。
【００６８】
ドライバパラメータの算出を行った結果、仕様１、仕様２、仕様３の順でドライバパラメータτLが大きくなることがわかった（図９参照）。また、この後に行った官能評価でも、仕様１が最も乗りやすいと評価した被験者が多いことがわかった（図１３参照）。更に、各仕様に対する被験者のコメントでも、被験者は仕様ごとに違いを感じていることがわかった（図１２参照）。被験者の中には、ヨー運動特性は変化していないにも関わらず（図６参照）、ヨー運動の変化についてコメントしている被験者がいた。
【００６９】
４．まとめ
以上の結果より、Ｇβ・Ｔβの値に応じて、ドライバパラメータ及び官能評価において性能評価が変わることがわかった。より詳しくは、Ｇβ・Ｔβの値が小さいほど（特に仕様１）、ドライバパラメータ及び官能評価の両方において、性能評価が良いことがわかった。したがって、本実施形態によれば、Ｇβ・Ｔβを可変パラメータとして適用することで、車体横すべり角の過渡特性を変更したシミュレーションを適切に行うことができることが見出された。
【００７０】
また、「２−２．車速に応じた車体横すべり角の特性」のセクションで述べたように、車速に応じてＧβ・Ｔβの値が増加する。加えて、「発明が解決しようとする課題」のセクションで述べたように、車体横すべり角はヨー運動と横運動に介在している。これらから、高速になるにつれて運転操作が難しくなることの要因として、ヨー運動と横運動との相対的な差（特に過渡的な差）が広がると考えられる。
【００７１】
５．適用例
次に、上記したシミュレーション方法の内容を実際の車両に適用した例について述べる。
【００７２】
図１４は、本発明による車両挙動の制御方法を適用した車両の概略構成図である。図１４に示すように、車両は、主に、ドライバによって操作されるステアリング２０と、前輪２３（２３ａ、２３ｂ）と、この前輪２３を転舵させるための前輪操舵機構２１と、後輪２４（２４ａ、２４ｂ）と、この後輪２４を転舵させるための後輪操舵機構２２と、前輪２３及び後輪２４を別個に転舵させるべく、前輪操舵機構２１及び後輪操舵機構２２を制御するコントローラ２５と、を有する。このように、車両は、四輪操舵（４ＷＳ（4 Wheel Steering））を実施可能に構成されている。例えば、コントローラ２５は、前輪２３ａ、２３ｂ、後輪２４ａ、２４ｂのそれぞれを舵角θF1、θF2、θR1、θR1に設定するように、前輪操舵機構２１及び後輪操舵機構２２の各々を制御する。
【００７３】
このような四輪操舵が可能なシステムでは、ヨー運動と横運動との関係を制御することができる。具体的には、横運動は、車両に働く横方向の力の和に起因する運動に相当し、ヨー運動は、車両において前側部分に働く横方向の力と後ろ側部分に働く横方向の力との差に起因する運動に相当するので、四輪操舵が可能なシステムによれば、前輪２３及び後輪２４のそれぞれの舵角を別個に変えることで、車両に働く横方向の力の差を変化させることができ、その結果、ヨー運動と横運動との関係を制御することができるのである。
【００７４】
したがって、本実施形態では、コントローラ２５は、現在の車速及び実舵角に基づき、上記したＧβ・Ｔβに応じて車両のヨー運動及び横運動を調整すべく、前輪操舵機構２１及び後輪操舵機構２２を制御する。つまり、コントローラ２５は、ヨー運動と横運動との関係が、現在の車速及び実舵角に応じた所望の関係となるように、前輪操舵機構２１及び後輪操舵機構２２を制御する。例えば、ヨー運動と横運動との所望の関係を実現するためのＧβ・Ｔβの値を、車速に対応付けたマップとして作成しておくと共に、このようなＧβ・Ｔβの値を実現するための前輪操舵機構２１及び後輪操舵機構２２の制御値（前輪２３及び後輪２４のそれぞれの実舵角に対応する）を規定したマップを作成しておき、コントローラ２５は、そのようなマップを参照して、前輪操舵機構２１及び後輪操舵機構２２に対する制御を行う。
【００７５】
１つの例では、上述したようなＧβ・Ｔβの値が小さいほど性能評価が良いという実験結果を受けて、コントローラ２５は、Ｇβ・Ｔβが小さくなるように、例えばＧβ・Ｔβが所定値（車速に応じて変化させるものとする）以下となるように、前輪操舵機構２１及び後輪操舵機構２２を制御する。
他の例では、コントローラ２５は、ヨー運動と横運動との過渡的な差が小さくなるように、車速に応じたＧβ・Ｔβに基づき、前輪操舵機構２１及び後輪操舵機構２２を制御する。この場合、ヨー運動と横運動との過渡的な差を小さくするためのＧβ・Ｔβは比較的小さな値となる。よって、当該制御は、上記したＧβ・Ｔβが小さくなるようにする制御と類似したものとなる。
【００７６】
以上述べた本実施形態によれば、Ｇβ・Ｔβを用いて車体横すべり角を調整して、ヨー運動と横運動との関係を適切に制御することができる。具体的には、Ｇβ・Ｔβを用いて車体横すべり角の過渡特性を調整して、ヨー運動と横運動との相対的な差を制御することができる。これにより、操舵により旋回するときの車両挙動に対するフィーリングを向上させることができる、つまりドライバが操舵したときに感じる操縦安定性（車両応答性など）を向上させることができる。その結果、本実施形態によれば、旋回時の修正操舵を低減することができ、操縦安定性が更に向上して、ドライビングフィーリングをより一層向上させることができる。
【００７７】
ここで、本実施形態による作用効果をより詳しく説明する。通常、ドライバは、車両の向きに対して車両が進んでいる方向を何らかの方法で認識しようとする。熟練したドライバは、経験により、車速が高いと、車両の向きに対して車両が外側に進んでいることを認識し、車速が低いと、車両の向きに対して車両が内側に進んでいることを認識する傾向にある。このとき、舵角変化が生じたときに、車両の向きに対する車両が進んでいる方向の変化が小さければ、遠方での自車の位置を予測しやすくなり、車両を意のままに操ることができるようになる。
上記した本実施形態によれば、Ｇβ・Ｔβを用いてヨー運動と横運動との関係を制御するので、操舵時（詳しくは舵角が変化しているとき）において車両の向きに対する車両が進んでいる方向の変化を小さくすることができる。したがって、本実施形態によれば、車両の操縦安定性が向上して、旋回時の修正操舵を低減することができるのである。
【００７８】
なお、上記した適用例では、本発明による車両挙動の制御方法を、四輪操舵が可能な車両に適用していたが、本発明の適用はこれに限定はされない。これ以外にも、本発明は、例えば、各輪の制駆動力を独立して制御できるよう構成された車両にも適用可能である。要は、各輪に働く力（横力だけだなく縦力も含む）の分布を制御できるように構成された種々の車両に、本発明を適用可能である。
【符号の説明】
【００７９】
１  ドライビングシミュレータシステム
２  ドライビングシミュレータ
４、２０  ステアリング
６  シート
８  表示部
１０、２５  コントローラ
１２  ヨー運動用アクチュエータ
１４  横運動用アクチュエータ
２１  前輪操舵機構
２２  後輪操舵機構
２３  前輪
２４  後輪
【図１】
2018171949.tif 000021
【図２】
2018171949.tif 000022
【図３】
2018171949.tif 000023
【図４】
2018171949.tif 000024
【図５】
2018171949.tif 000025
【図６】
2018171949.tif 000026
【図７】
2018171949.tif 000027
【図８】
2018171949.tif 000028
【図９】
2018171949.tif 000029
【図１０】
2018171949.tif 000030
【図１１】
2018171949.tif 000031
【図１２】
2018171949.tif 000032
【図１３】
2018171949.tif 000033
【図１４】
2018171949.tif 000034