(57)【特許請求の範囲】
【請求項１】
モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、
前記モータのトルクＴを測定するモータトルク測定手段と、
駆動輪の回転速度ωおよび回転加速度
【数１】
0004936552.tif 000031
（ωｄｏｔ）を算出する車両モデル演算手段と、
前記モータトルク測定手段が測定したトルクＴと前記車両モデル演算手段が算出した回転速度ωおよび回転加速度ωｄｏｔを用いて、スリップ率に関する常微分方程式（Ａ）を計算することにより推定スリップ率
【数２】
0004936552.tif 000032
（λｈａｔ）を算出するスリップ率演算手段と
を備えたことを特徴とするスリップ率推定装置。
【数３】
0004936552.tif 000033
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【請求項２】
モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、
前記モータのトルクＴを測定するモータトルク測定手段と、
前記自動車の車体加速度ａxを測定する加速度測定手段と、
駆動輪の回転速度ωおよび回転加速度ωｄｏｔを算出する車両モデル演算手段と、
前記自動車の駆動力Ｆdの推定値を算出する駆動力演算手段と、
前記駆動力演算手段が算出した推定駆動力
【数４】
0004936552.tif 000034
（Ｆdｈａｔ）と前記加速度測定手段が測定した前記車体加速度ａxとから走行抵抗Ｆdrの推定値を算出する走行抵抗演算手段と、
前記モータトルク測定手段が測定したトルクＴと、前記車両モデル演算手段が算出した回転速度ωおよび回転加速度ωｄｏｔと、前記走行抵抗演算手段が算出した推定走行抵抗
【数５】
0004936552.tif 000035
（Ｆdrｈａｔ）を用いて、スリップ率に関する常微分方程式（Ｂ）を計算することにより推定スリップ率λｈａｔを算出するスリップ率演算手段と
を備えたことを特徴とするスリップ率推定装置。
【数６】
0004936552.tif 000036
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【請求項３】
モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、
予め指定された駆動輪の回転速度から前記モータに対するトルク指令を演算する手段と、
前記モータのトルクを当該トルク指令に基づき制御する手段と、
請求項１又は２に記載のスリップ率推定装置を備え、当該スリップ率推定装置が算出する推定スリップ率が所望の値をとるようにモータの当該トルクを制御することを特徴とするスリップ率制御装置。
【請求項４】
モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、
予め指定された駆動輪の回転速度から前記モータに対するトルク指令を演算する手段と、
前記モータのトルクを前記トルク指令に基づき制御する手段と、
入力される駆動輪の回転速度から目標トルクＴを算出する比例積分制御手段と、
前記目標トルクＴから前記駆動輪の回転速度ω及び回転加速度ωｄｏｔを算出する車両モデル演算手段と、
前記比例積分制御手段が算出した目標トルクＴ、前記車両モデル演算手段が算出した駆動輪の回転速度ω及び回転加速度ωｄｏｔを用いてスリップ率に関する常微分方程式（Ａ）を計算して推定スリップ率λｈａｔを算出するスリップ率演算手段と、
前記スリップ率推定装置が算出した推定スリップ率から車体速度を算出し、当該車体速度から目標スリップ率に対する駆動輪の目標回転速度を算出する車輪速度演算手段とを備え、
前記比例積分制御手段は、前記スリップ率演算手段が算出した推定スリップ率に応じて当該制御ゲインを変化させ、前記トルク指令を演算する手段は、当該比例積分制御手段から算出された前記目標トルクＴに基づいて演算することを特徴とするスリップ率制御装置。
【数７】
0004936552.tif 000037
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【請求項５】
モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、
前記モータに対するトルク指令を演算する手段と、
前記モータのトルクＴを当該トルク指令に基づき制御する手段と、
前記自動車の車体加速度ａxを測定する加速度測定手段と、
入力される駆動輪の回転速度から目標トルクＴを算出する比例積分制御手段と、
前記目標トルクＴから前記駆動輪の回転速度ω及び回転加速度ωｄｏｔを算出する車両モデル演算手段と、
前記自動車の駆動力Ｆdの推定値を算出する駆動力演算手段と、
前記駆動力演算手段が算出した推定駆動力Ｆdｈａｔと前記加速度測定手段が測定した前記車体加速度ａxとから走行抵抗Ｆdrの推定値を算出する走行抵抗演算手段と、
前記比例積分制御手段が算出したトルクＴと、前記車両モデル演算手段が算出した駆動輪の回転速度ωおよび回転加速度ωｄｏｔと、前記走行抵抗演算手段が算出した推定走行抵抗Ｆdrｈａｔを用いて、スリップ率に関する常微分方程式（Ｂ）を計算することにより推定スリップ率λｈａｔを算出するスリップ率演算手段と
前記スリップ率演算手段が算出した推定スリップ率λｈａｔから車体速度を算出し、当該車体速度から目標スリップ率λ*に対する駆動輪の目標回転速度を算出する車輪速度演算手段とを備え、
前記比例積分制御手段は、前記スリップ率演算手段が算出した推定スリップ率に応じて当該制御ゲインを変化させ、前記トルク指令を演算する手段は、当該比例積分制御手段から算出された前記トルクＴに基づいて演算することを特徴とするスリップ率制御装置。
【数８】
0004936552.tif 000038
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【請求項６】
モータのトルクで駆動輪を駆動する自動車において、請求項１又は２に記載のスリップ率推定装置を備え、当該スリップ率推定装置が算出する推定スリップ率が所望の値をとるようにモータの前記トルクを制御することを特徴とする自動車。
【請求項７】
モータのトルクで駆動輪を駆動する自動車において、モータの当該トルクを制御する請求項３乃至５のいずれかに記載のスリップ率制御装置を備えたことを特徴とする自動車。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スリップ率推定装置及びスリップ率制御装置に関し、より詳細には、スリップ率の時間変化を考慮した車体速度検出不要な減速時におけるスリップ率推定装置及びスリップ率制御装置に関する。
【背景技術】
【０００２】
現在、モータを駆動力として利用した電気自動車に注目が集まっている。これは、燃費・環境問題の他にも、モータを駆動力とすることにより得られる利点がある。その利点として以下の３つが挙げられる。
・モータがエンジンに比べて小さいため、各車輪に分散配置が可能
・モータがエンジンに対して、トルク応答が数百倍速い
・発生トルクを正確に把握できる
【０００３】
電気自動車は環境問題の視点から大きな注目を集めているが、以上のようなモータ特性により車両制御の点からも非常に有利であるといえる。上記の優位点を利用した、トラクションコントロールの研究が発表されている（非特許文献１、２、３参照）。
【０００４】
トラクション制御を行う際、車輪速度と車体速度によって定義されるスリップ率は非常に重要な値であるが、車輪速度に比べて車体速度は測定が困難な場合が多い。通常は非駆動輪の車輪速度、加速度センサの積分、光学センサ等から測定する。しかし、非駆動輪を用いた測定においては、機械ブレーキを使って減速した場合、制動力が四輪全てに働くため、非駆動輪が存在せず車体速度は測定できない。また四輪駆動車であればもともと非駆動輪が存在しないので、この場合も車体速度を測定することができない。これらを解消するために非駆動の第五輪を取付けるのは一般車では現実的でない。加速度センサの積分から求める場合、加速度センサの値にオフセットがあると、オフセットを積分し続けるので推定値が発散してしまう。また、光学センサは信頼できるが、高価であり実用向きでない。
【０００５】
従って、検出困難な車体速度を用いることなく正確なスリップ率を得るスリップ率推定法は極めて実用的であるといえる。非特許文献１ではスリップ率を制御しトラクションコントロールを行っているが、検出困難な車体速度を必要としている。非特許文献２では車体速度情報を用いていないが、空転検出を行うのみでスリップ率の過渡的変化は考慮していない。これに対して発明者らのグループは非特許文献３で、駆動時（加速時）における車体速度検出不要なスリップ率推定法と制御法を提案した。
【０００６】
【非特許文献１】佐渡秀夫、坂井真一郎、内田利之、堀洋一:「路面状態推定とスリップ率制御を用いた電気自動車の高性能トラクションコントロール」、平成１０ 年電気学会産業応用部門全国大会講演論文集、１９９８年、３、ｐｐ３２１−３２４
【非特許文献２】坂井真一郎、佐渡秀夫、堀洋一:「電気自動車における車体速度情報不要の新しいタイヤ空転検出法」、電学論、２０００年、Ｄ．１２０、２、ｐｐ２８１−２８７
【非特許文献３】藤井淳、藤本博志:「車体速度検出不要なスリップ率推定に基づく電気自動車のトラクション制御に関する研究」、横浜国立大学学位論文、２００７年
【非特許文献４】小竹元基、永井正夫:「操安性向上を目指した超小型電気自動車の車輪駆動トルク制御」、自動車技術会論文集、２００３年、３４、４、ｐｐ１６９−１７４
【非特許文献５】H.B.Pacejka、 and E.Bakker、 “The Magic FormulaTyre Model”、 Tyre models for vehicle dynamic analysis:proceedings of the 1st International Colloquium onTyre Models for Vehicle Dynamics Analysis、 held inDelft、 The Netherlands、 Oct 21-22 (1991)
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、駆動時と制動時ではスリップ率の定義が異なるため、制動時において非特許文献３の手法を直接利用することはできない。
【０００８】
本発明は、このような課題に鑑みてなされたもので、その目的とするところは、車体速度の検出を必要としない制動時のスリップ率推定装置及びスリップ率制御装置を提供することにある。
【課題を解決するための手段】
【０００９】
このような目的を達成するために、請求項１に記載の発明は、モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、前記モータのトルクＴを測定するモータトルク測定手段と、駆動輪の回転速度ωおよび回転加速度
【００１０】
【数１】
0004936552.tif 000002
【００１１】
（ωｄｏｔ）を算出する車両モデル演算手段と、前記モータトルク測定手段が測定したトルクＴと前記車両モデル演算手段が算出した回転速度ωおよび回転加速度ωｄｏｔを用いて、スリップ率に関する常微分方程式（Ａ）を計算することにより推定スリップ率
【００１２】
【数２】
0004936552.tif 000003
【００１３】
（λｈａｔ）を算出するスリップ率演算手段とを備えたことを特徴とする。
【００１４】
【数３】
0004936552.tif 000004
【００１５】
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【００１６】
請求項２に記載の発明は、モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、前記モータのトルクＴを測定するモータトルク測定手段と、前記自動車の車体加速度ａxを測定する加速度測定手段と、駆動輪の回転速度ωおよび回転加速度ωｄｏｔを算出する車両モデル演算手段と、前記自動車の駆動力Ｆdの推定値を算出する駆動力演算手段と、前記駆動力演算手段が算出した推定駆動力
【００１７】
【数４】
0004936552.tif 000005
【００１８】
（Ｆdｈａｔ）と前記加速度測定手段が測定した前記車体加速度ａxとから走行抵抗Ｆdrの推定値を算出する走行抵抗演算手段と、前記トルク測定手段が測定したトルクＴと、前記車両モデル演算手段が算出した回転速度ωおよび回転加速度ωｄｏｔと、前記走行抵抗演算手段が算出した推定走行抵抗
【００１９】
【数５】
0004936552.tif 000006
【００２０】
（Ｆdrｈａｔ）を用いて、スリップ率に関する常微分方程式（Ｂ）を計算することにより推定スリップ率λｈａｔを算出するスリップ率演算手段とを備えたことを特徴とする。
【００２１】
【数６】
0004936552.tif 000007
【００２２】
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【００２３】
請求項３に記載の発明は、モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、予め指定された駆動輪の回転速度から前記モータに対するトルク指令を演算する手段と、前記モータのトルクを当該トルク指令に基づき制御する手段と、請求項１又は２に記載のスリップ率推定装置を備え、当該スリップ率推定装置が算出する推定スリップ率が所望の値をとるようにモータの当該トルクを制御することを特徴とする。
【００２４】
請求項４に記載の発明は、モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、予め指定された駆動輪の回転速度から前記モータに対するトルク指令を演算する手段と、前記モータのトルクを前記トルク指令に基づき制御する手段と、入力される駆動輪の回転速度から目標トルクＴを算出する比例積分制御手段と、前記目標トルクＴから前記駆動輪の回転速度ω及び回転加速度ωｄｏｔを算出する車両モデル演算手段と、前記比例積分制御手段が算出した目標トルクＴ、前記車両モデル演算手段が算出した駆動輪の回転速度ω及び回転加速度ωｄｏｔを用いてスリップ率に関する常微分方程式（Ａ）を計算して推定スリップ率λｈａｔを算出するスリップ率演算手段と、前記スリップ率演算手段が算出した推定スリップ率から車体速度を算出し、当該車体速度から目標スリップ率に対する駆動輪の目標回転速度を算出する車輪速度演算手段とを備え、前記比例積分制御手段は、前記スリップ率演算手段が算出した推定スリップ率に応じて当該制御ゲインを変化させ、前記トルク指令を演算する手段は、当該比例積分制御手段から算出された前記目標トルクＴに基づいて演算することを特徴とする。
【００２５】
【数７】
0004936552.tif 000008
【００２６】
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【００２７】
請求項５に記載の発明は、モータのトルクで駆動輪を駆動する自動車に用いるスリップ率測定装置において、前記モータに対するトルク指令を演算する手段と、前記モータのトルクＴを前記トルク指令に基づき制御する手段と、前記自動車の車体加速度ａxを測定する加速度測定手段と、入力される駆動輪の回転速度から目標トルクＴを算出する比例積分制御手段と、前記目標トルクＴから前記駆動輪の回転速度ω及び回転加速度ωｄｏｔを算出する車両モデル演算手段と、前記自動車の駆動力Ｆdの推定値を算出する駆動力演算手段と、前記駆動力演算手段が算出した推定駆動力Ｆdｈａｔと前記加速度測定手段が測定した前記車体加速度ａxとから走行抵抗Ｆdrの推定値を算出する走行抵抗演算手段と、前記比例積分制御手段が算出したトルクＴと、前記車両モデル演算手段が算出した駆動輪の回転速度ωおよび回転加速度ωｄｏｔと、前記走行抵抗演算手段が算出した推定走行抵抗Ｆdrｈａｔを用いて、スリップ率に関する常微分方程式（Ｂ）を計算することにより推定スリップ率λｈａｔを算出するスリップ率演算手段と前記スリップ率演算手段が算出した推定スリップ率λｈａｔから車体速度を算出し、当該車体速度から目標スリップ率λ*に対する駆動輪の目標回転速度を算出する車輪速度演算手段とを備え、前記比例積分制御手段は、前記スリップ率演算手段が算出した推定スリップ率に応じて当該制御ゲインを変化させ、前記トルク指令を演算する手段は、当該比例積分制御手段から算出された前記トルクＴに基づいて演算することを特徴とする。
【００２８】
【数８】
0004936552.tif 000009
【００２９】
（ｒ：駆動輪のタイヤ半径、Ｍ：車両重量、Ｊω：駆動輪回転部慣性モーメント）
【００３０】
請求項６に記載の発明は、モータのトルクで駆動輪を駆動する自動車において、請求項１又は２に記載のスリップ率推定装置を備え、当該スリップ率推定装置が算出する推定スリップ率が所望の値をとるようにモータの前記トルクを制御することを特徴とする。
【００３１】
請求項７に記載の発明は、モータのトルクで駆動輪を駆動する自動車において、モータの当該トルクを制御する請求項３乃至５のいずれかに記載のスリップ率制御装置を備えたことを特徴とする。
【発明の効果】
【００３２】
本発明によれば、車体速度の検出なしに制動時のスリップ率の推定が可能になる。さらに推定されたスリップ率を用いて駆動輪を制御することにより、駆動輪のスリップ率を最適スリップ率に制御することができれば制動距離を短くすることができる。また、車輪が路面に粘着するので、非特許文献４のようなヨーモーメント制御の駆動力が確実に路面に伝わるため、２次元方向の運動も安定化され、車両がスピンする危険性を低減することができる。
【発明を実施するための最良の形態】
【００３３】
電気自動車は、駆動輪に接続されたモータの駆動力が駆動輪に伝達されて走行する。電源から出力された電流は、インバータを介してモータに供給される。モータは、制御手段としての電子制御ユニット（以下「ＥＣＵ」と称す）を介して電気的に接続されており、モータの出力は、ＥＣＵからの指令に基づいて制御される。ＥＣＵは、ＣＰＵ、ＲＯＭ、ＲＡＭ、入出力ポート、および記憶装置等を含むものである。またＥＣＵには、モータの発生トルクを測定するトルク測定器、モータに装備された位置センサ、車体に生じる加速度を測定する加速度センサが電気的にそれぞれ接続されている。
【００３４】
以下に、ECUに実装されるスリップ率制御装置について説明する。
【００３５】
本発明の構成は以下の通りである。まず１節で車両の運動方程式を示し、それに基づいて車両モデルのブロック図を構成する。そしてスリップ率推定法の比較対象（従来技術）を〈２・１〉節で、本発明の一実施形態を〈２・２・１〉節で示し、〈２・２・２〉節で提案法の別の実施形態を示す。また３節でシミュレーション・実験から本発明の手法の有効性を示し、４節でスリップ率を制御するシミュレーションを行う。
【００３６】
本明細書においては、便宜上、Ｘ（Ｘは任意の文字）の時間微分を表す
【００３７】
【数９】
0004936552.tif 000010
【００３８】
の代わりにＸｄｏｔ、Ｘの推定を表す
【００３９】
【数１０】
0004936552.tif 000011
【００４０】
の代わりにＸｈａｔの表記も用いる。
【００４１】
以下、図面を参照しながら本発明の実施形態について詳細に説明する。
１．車両の運動方程式
図１に、車両の各部に働く力を模式的に示す。モータの時定数が非常に小さいと仮定するとき、後輪二輪駆動の場合、車両に働く力は図１のようになり、車両の運動方程式は以下の式で表現できる。
【００４２】
【数１１】
0004936552.tif 000012
【００４３】
【数１２】
0004936552.tif 000013
【００４４】
【数１３】
0004936552.tif 000014
【００４５】
それぞれＪωi：車輪回転部慣性モーメント、ωi：モータの回転速度、Ｔi：モータトルク、ｒ：タイヤ半径、Ｆi：駆動力、Ｆdr：走行抵抗、Ｍ：車体重量、Ｖ：車体速度、Ｖωi：車輪速度である。ただし、ｉ＝ｒ、ｌであり、左右のそれぞれのモータの状態を表す。本発明では機械ブレーキを使用せず、モータの逆トルクのみで減速することを考える。そのため、式（１）には機械ブレーキによるトルクの影響が含まれていない。もしブレーキトルクが検出あるいは推定できれば、本発明で提案するスリップ率の推定法は、機械ブレーキを用いたものに容易に拡張できる。
【００４６】
ここでスリップ率λは次式で表される。
【００４７】
【数９】
0004936552.tif 000015
【００４８】
式（４）の分母はＶω、Ｖの大小関係により変わる。駆動時にはＶω＞Ｖなのでｍａｘ（Ｖω、Ｖ、ε）＝Ｖωであり、制動時にはＶω＜Ｖなのでｍａｘ（Ｖω、Ｖ、ε）＝Ｖである。またε（＜＜１）はＶωとＶが共に零の場合の零割を防ぐための小さな定数である。
【００４９】
図２に、タイヤと路面間の摩擦係数μとスリップ率λ両者の関係を表すμ−λ曲線を示す。この曲線は路面状況によって大きく異なるが、殆どの場合、駆動時（制動時）では特定のスリップ率において摩擦係数が最大値（最小値）をとり、スリップ率がそれより大きく（小さく）なると減少（増加）する。この曲線を表す方法として、今回は実験により得られたデータに近い曲線を描くように選ばれた方程式であるＭａｇｉｃ  Ｆｏｒｍｕｌａを用いる（非特許文献５参照）。
【００５０】
【数１５】
0004936552.tif 000016
【００５１】
この式はシミュレーションにおける車両モデルブロックにのみ用い、スリップ率推定には用いない。
【００５２】
路面とタイヤ間の摩擦力は摩擦係数μに垂直抗力Ｎをかけて得られる。
【００５３】
【数１６】
0004936552.tif 000017
【００５４】
式（１）〜（６）より図３のような左右それぞれの車輪を考慮した左右二輪駆動の車両モデルを得ることができる。
【００５５】
２．スリップ率推定法
〈２・１〉駆動力オブザーバによる推定法（従来技術）
〈２・１・１〉推定式の導出
ここでは式（２）第３項の走行抵抗を無視できるとし、Ｖについて解き、積分すれば以下のように車体速度を推定することができる。
【００５６】
【数１７】
0004936552.tif 000018
【００５７】
ここで式（１）より、モータにとって駆動力Ｆiは外乱とみなせるので、図４のような外乱オブザーバにより駆動力を推定する。これを駆動力オブザーバという（非特許文献１〜３参照）。
【００５８】
式（７）および式（４）を用いてスリップ率を求める。この推定法をＤＦＯ（Ｄｒｉｖｉｎｇ  Ｆｏｒｃｅ  Ｏｂｓｅｒｖｅｒ）と呼ぶ。
〈２・１・２〉推定誤差
実際には真値の車体速度と推定した車体速度には誤差がある。ここで車体速度の真値と推定値の誤差をｅv、スリップ率の真値と推定値の誤差をｅλiとすると以下のようになる。
【００５９】
【数１８】
0004936552.tif 000019
【００６０】
【数１９】
0004936552.tif 000020
【００６１】
推定した車体速度には走行抵抗が含まれていないため、推定する時間がたつにつれて積分によって推定した車体速度の誤差は拡大していく。それに伴い、スリップ率の推定にも誤差が生じる。
〈２・２〉スリップ率の時間変化を考慮した推定法（本発明）
〈２・２・１〉推定式の導出
式（４）をＶについて解き、この両辺を時間で微分し、式（１）、（２）、（３）を代入してＶω、Ｆiを消去し、λｄｏｔについて解くと次式が得られる。
【００６２】
【数２０】
0004936552.tif 000021
【００６３】
また、スリップ率推定装置は走行抵抗が小さい（Ｆdr≒０）と仮定し以下のように構成する。
【００６４】
【数２１】
0004936552.tif 000022
【００６５】
この推定法をＳＲＥ（Ｓｌｉｐ  Ｒａｔｉｏ  Ｅｓｔｉｍａｔｏｒ）と呼ぶ。
【００６６】
図１２に、本発明の一実施形態に係るスリップ率推定装置を示す。このスリップ率推定装置は、モータのトルクで駆動輪を駆動する自動車に実装され、駆動輪の回転速度ωおよび回転加速度ωｄｏｔを算出する車両モデル演算装置と、モータトルク測定手段が測定したトルクＴと車両モデル演算装置が算出した回転速度ωおよび回転加速度ωｄｏｔを用いて、スリップ率λｈａｔを算出するスリップ率演算装置（ＳＲＥ）からなる。
【００６７】
しかしこの推定法は走行抵抗を十分小さいものとし無視しているため、走行抵抗が無視できない場合正確な推定を行うことができない。そこで次節で走行抵抗を推定し、走行抵抗を考慮したスリップ率推定法について説明する。
【００６８】
〈２・２・２〉走行抵抗推定
式（２）より、車体にとって走行抵抗Ｆdrは外乱とみなせるので、駆動力Ｆiと加速度ａxから走行抵抗Ｆdrを推定する。駆動力Ｆiは図４の駆動力オブザーバにより推定したＦiｈａｔを用いる。その結果、図５のような外乱オブザーバを２段組んだようなブロック図となる（非特許文献３参照）。また、図５のＤＦＯ内のＬＰＦと加速度にかかるＬＰＦは同一のものである。
【００６９】
走行抵抗演算手段によって走行抵抗が推定できれば、式（１１）は走行抵抗を考慮した以下のような式に書き換えることができる。
【００７０】
【数２２】
0004936552.tif 000023
【００７１】
これをＤＲＥ−ＳＲＥ（Ｄｒｉｖｉｎｇ  Ｒｅｓｉｓｔａｎｃｅ  Ｅｓｔｉｍａｔｏｒ−ＳｌｉｐＲａｔｉｏ  Ｅｓｔｉｍａｔｏｒ）と呼ぶ。
【００７２】
このＤＲＥ−ＳＲＥに基づくスリップ率推定装置は、モータのトルクで駆動輪を駆動する、車体加速度ａxを測定する加速度測定手段を備えた自動車に実装され、駆動輪の回転速度ωおよび回転加速度ωｄｏｔを算出する車両モデル演算手段、自動車の駆動力Ｆdの推定値を算出する駆動力演算手段、走行抵抗演算手段、スリップ率演算装置を備えている。走行抵抗演算手段は、駆動力演算手段が算出した推定駆動力Ｆdｈａｔと前記加速度測定手段が測定した前記車体加速度ａxとから走行抵抗Ｆdrの推定値を算出する。そして、スリップ率演算手段が、モータトルク測定手段が測定したトルクＴと車両モデル演算手段が算出した回転速度ωおよび回転加速度ωｄｏｔと、走行抵抗演算手段が算出した推定走行抵抗Ｆdrｈａｔを用いて、スリップ率λｈａｔを算出する。
【００７３】
〈２・２・３〉推定誤差
ここで式（１０）、（１２）より推定誤差は、以下のようになる。またＦdrｈａｔの推定値はＤＦＯの時定数に従って真値に収束する。よってＦdr≒Ｆdrｈａｔとすると、
【００７４】
【数２３】
0004936552.tif 000024
【００７５】
となり、これを式（１）、（２）等からＶωiｄｏｔ、Ｖωi、Ｖｄｏｔを用いて書き換えると以下のようになる。
【００７６】
【数２４】
0004936552.tif 000025
【００７７】
よって、以下のような条件を満たせば誤差は収束する。
【００７８】
【数２５】
0004936552.tif 000026
【００７９】
制動時の場合スリップ率の値がとる範囲が−１＜λ＜０であるため、式（１５）より式（１４）は車輪の減速度が車体の減速度の少なくとも２倍よりも大きいときに誤差が収束することを意味する。それは即ち車輪が路面上を空転しているときであり、このときλiもλiｈａｔも−１に近づけば、式（１５）から必ず収束する。そして一度収束すれば走行抵抗Ｆdrに起因する推定誤差がない限り、正確なスリップ率を推定し続けることができる。
【００８０】
３．スリップ率推定
〈３・１〉シミュレーションによる検証
〈３・１・１〉シミュレーション条件
前章までに説明した推定法に関して、収束性の理論検討のためＳＲＥとＤＲＥ−ＳＲＥは初期誤差がない場合と、初期誤差がある場合とでスリップ率推定のシミュレーションを行う。実際に電気自動車に適用する場合、発明者らが以前に提案した加速時の推定法（非特許文献３参照）と本発明の減速時の推定法を切り替えて使うことになる。その切り替えの時に誤差が生じる可能性があるため、初期誤差がある場合の検討は実用上も意味があるものである。但し、本願明細書では簡単化のため減速時のスリップ率のみを推定し、加速時のスリップ率は推定しない。またＤＦＯによる推定法も、実用上は加速時から式（７）によって車体速度を推定するので、今回のシミュレーションは加速時に生じる車体速度の誤差を考慮し、初期誤差をもたせた検討を行っている。
【００８１】
シミュレーションは左右の車輪速度、トルク指令、スリップ率が等しいものとして左右の区別なく求めている。プラントのパラメータは実機より得られた値であるＪω＝１．２６ｋｇｍ2、Ｍ＝７００ｋｇ、ｒ＝０．３０２ｍとする。トルク指令値は−１１０Ｎｍを与える。路面状況はシミュレーション開始から１．７５ｓｅｃまで低μ路（μmax＝０．２）とし、それ以降は乾燥路（μmax＝１．０）とした。走行抵抗は実験における推定値の最大値が２００Ｎであると得られたので、それを一定値として与えた。ＳＲＥとＤＲＥ−ＳＲＥは初期誤差がない場合には車体速度、車輪速度が共に４．００ｍ／ｓｅｃから推定を開始し、初期誤差がある場合には車体速度が４ｍ／ｓｅｃ、車輪速度が３．９６ｍ／ｓｅｃとなったところから推定を開始する。式（７）における車体速度の積分の初期値は、実機で４．００ｍ／ｓｅｃまで急加速したときに約１割程度の誤差が確認されたため、車体速度が加速時に生じる車体速度の誤差を考慮して４．４０ｍ／ｓｅｃを与えた。
【００８２】
〈３・１・２〉ＤＦＯによる推定法（従来技術）
図６（ｂ）に、ＤＦＯに基づくスリップ率の真値と推定値を示す。スリップ率の絶対値が大きいところでは車輪速度が車体速度に比べて小さくなりλ≒−Ｖ／Ｖとなるので推定値は−１に収束する。しかしＶｈａｔの初期誤差によって、スリップ率の絶対値が小さなところでは大きな誤差が確認できる。積分する時間が長くなるにつれ図６（ａ）の誤差は大きくなるので推定法としては問題がある。そして車体速度が零になったとき走行抵抗の影響によって推定した車体速度は零にはならないという問題もある。
【００８３】
〈３・１・３〉スリップ率の時間変化を考慮したスリップ率推定法（本発明）
図６（ｃ）に初期誤差無しのＳＲＥに基づくスリップ率の真値と推定値を示し、図６（ｄ）に初期誤差有りのＳＲＥに基づくスリップ率の真値と推定値を示す。このように走行抵抗を考慮していない場合、低μ路において空転しているとき、すなわちスリップ率の絶対値が大きなところでは推定値が真値に漸近することが確認できるが、高μ路に突入しタイヤが粘着すると式（１４）を満たさなくなり、推定誤差が拡大することが確認できる。
【００８４】
図６（ｅ）に初期誤差無しのＤＲＥ−ＳＲＥに基づくスリップ率の真値と推定値を示し、図６（ｆ）に初期誤差無しのＤＲＥ−ＳＲＥに基づくスリップ率の真値と推定値を示す。このように走行抵抗を考慮した場合、初期誤差があると、高μ路突入までに収束しきっておらず、その影響で誤差が収束しない。しかしＤＦＯとＳＲＥよりは大きく改善されている。また初期誤差がなければたとえ高μ路に突入し、式（１４）の安定条件を満たさなくなったとしても、正確な推定が行える。
〈３・２〉オフライン実験による検証
【００８５】
〈３・２・１〉実験車両
今回の実験には本研究室で製作した電気自動車ＦＰＥＶ２−Ｋａｎｏｎを用いた。ＦＰＥＶ２−Ｋａｎｏｎには東洋電機製アウターロータ型インホイールモータを後輪２輪に装着している。本モータはダイレクトドライブ方式であり、減速ギヤによるバックラッシュの影響がない。従って反力情報がギヤで失われることなくモータ側にもどるため、本研究室で提案している各種推定法を行うにあたり非常に有効であると考えられる。表１にモータのスペックを示す。
【００８６】
【表１】
0004936552.tif 000027
【００８７】
車両制御のためのコントローラはｄＳＰＡＣＥ社のＡＵＴＯＢＯＸ−ＤＳ１１０３を採用した。本コントローラは耐振動性に優れ、８Ｖ〜６０Ｖと非常に幅広い電圧で起動し、インパネモニターによりリアルタイムでの各測定値、推定値のグラフの表示やパラメータチューニングが可能である。
【００８８】
〈３・２・２〉実験結果
今回提案した手法に関して、オフライン実験による検証を行った。低μ路はアルミの板を地面に敷き、そこに洗剤と水をまくことによって実現を図っている。モータが搭載されている後輪がアルミの板に乗るまで加速をし、その後左右どちらかの車輪速度が零になるまで−１８０Ｎｍ一定のトルク指令を与え、車輪速度が一度零になった後はトルク指令を与えない。車体速度はＣＯＲＲＳＹＳ−ＤＡＴＲＯＮ社製ＣＯＲＲＥＶＩＴ  Ｓ−４００の光学式センサを用いることにより測定する。ただし測定した車体速度は検証のみに使用し、推定には使用しない。
【００８９】
〈３・１〉節では左右の車輪速度、トルク、スリップ率が等しいとしたが実際には左右に差があるため、左右それぞれのスリップ率を推定し、検証する。シミュレーション同様、減速時のスリップ率のみ推定し、加速時のスリップ率は推定しない。図７（ｂ）、図７（ｃ）、図７（ｇ）の推定はトルク指令が負になった瞬間から推定を行っているため、車輪速度が車体速度よりもまだ大きいところから推定を開始しており初期誤差を持っている。図７（ｄ）、図７（ｈ）の推定はスリップ率の真値が零になった瞬間から推定を行った。しかし図７（ｄ）、７（ｈ）の推定は検証用であり、実際にはこのようなことはできない。また今回センサ等の不具合により走行抵抗の推定が行えなかったため、オフライン実験ではＳＲＥでの推定のみを行い、ＤＲＥ−ＳＲＥでの推定は行わないものとする。
【００９０】
図７（ｂ）と図７（ｆ）は、式（４）と図９に基づいてスリップ率推定および車輪速度制御を行っているため、車体速度の推定値の誤差拡大と共にスリップ率の推定値の誤差も拡大してしまっている。図９のように式（７）に基づき算出された車体速度は、時間がたつにつれ誤差が拡大してしまうため、ＤＦＯはスリップ率推定には不向きである。
【００９１】
図７（ｃ）と図７（ｇ）は、トルク指令が負になった瞬間から推定するという方法をとっているため、真値に対して初期誤差がある状態から推定を行っている。その結果スリップ率の絶対値が大きなところでは誤差が収束し、粘着すると誤差を持つ。
【００９２】
図７（ｄ）と図７（ｈ）は、初期誤差がないため空転しているときは正確な推定が行えている。粘着すると走行抵抗に起因する推定誤差が残ってしまっているが、ＤＦＯによる推定法より推定結果がよくなっている。
【００９３】
今回走行抵抗の推定が行えなかったが、走行抵抗推定を行いＤＲＥ−ＳＲＥで推定すれば、ＳＲＥで確認された推定誤差は収束し、さらに良好な結果が得られると考えられる。
【００９４】
４．スリップ率制御
〈４・１〉車輪速度制御
図６（ｅ）のようなスリップ率推定ができるときのスリップ率制御のシミュレーションを行う。図２に示すように、スリップ率が−０．２付近で摩擦係数が最小値をとるため、このとき最大の制動力が得られる。この最適スリップ率になるように車輪速度を制御すれば、制御無しで急制動し、タイヤがロックしてしまったときよりも制動距離が縮む。
【００９５】
非特許文献３で行われている車輪速度制御によるスリップ率制御を行う。車輪速度が検出できるため、スリップ率の推定は車体速度の推定と等価である。スリップ率推定値を基に式（４）より式（１６）のように車体速度の推定値を求めることができ、式（１７）より最適スリップ率に対する車輪角速度指令がわかる。これにより回転速度制御を用いることでスリップ率制御を行う。
【００９６】
【数２６】
0004936552.tif 000028
【００９７】
【数２７】
0004936552.tif 000029
【００９８】
ω*は車輪速度の指令値であり、λ*はスリップ率の指令値である。この場合、λ*＝−０．２とすることになる。図１０に、車輪速度制御ブロック図を示す。車輪速度制御器はＰＩ制御器を用い、プラントは車輪の慣性モーメントのみを考慮した以下の式とし、極配置法によって制御ゲインを決定する。
【００９９】
【数２８】
0004936552.tif 000030
【０１００】
速度制御系の極は５０ｒａｄ／ｓｅｃとした。
【０１０１】
〈４・２〉シミュレーションによる検証
シミュレーション開始と共に制御をかけ、スリップ率を−０．２に収束させる。路面状態は常に低μ路（μmax＝０．２）する。また、車体速度と車輪速度の初期値は８ｍ／ｓｅｃ（＝２８．８ｋｍ／ｈ）とする。その他の条件は〈３・１・１〉節と同じものを用いる。
【０１０２】
シミュレーション結果より、推定したスリップ率を用いて最適スリップ率に制御できていることがわかる。制御無しの場合と制御有りの場合の車体速度を比較すると制御有りのほうが速く減速していることが確認できる。図１１（ｃ）、１１（ｄ）より制動距離をそれぞれ求めると、制御有りが約２５．３ｍであり制御無しが２６．７ｍであり制御有りのほうが１．４ｍ程度制動距離が短くなっている。但し、低μ路の最適スリップ率における摩擦係数と、スリップ率−１における摩擦係数にはそれほど大きな差がないため制動距離があまり縮まらないと考えられる。しかし制動距離があまり変わらなかったとしてもタイヤがロックしているか粘着しているかは大きな違いがある。今回のシミュレーション結果では確認できないが、タイヤが滑り状態の場合車体が大きくスピンする可能性がある。低μ路では制動距離にそれほど大きな違いが表れなかったが、高μ路の最適スリップ率における摩擦係数と、スリップ率−１における摩擦係数には低μ路の場合よりも大きな差があるため、高μ路における急制動になればこの違いは顕著になってくると予想される。
【０１０３】
本実施形態においては、車輪速度制御の比例積分制御を採用したが、これは一例であって、トルク制御に基づく制御などの他の制御方法を用いることもできる。本発明で重要なことは、車体速度を測定せずに制動時のスリップ率の推定を高精度に行うことができることにある。そのため、本発明によって推定されたスリップ率を用いてスリップ率の制御を行うスリップ率制御装置は、その制御方法に関わりなく、より高精度なスリップ率制御を行うことができる。
【０１０４】
５．まとめ
本願明細書を作成するにあたり、スリップ率推定には不可欠であった車体速度を用いないスリップ率推定法を提案し、シミュレーションとオフライン実験を行った。駆動力オブザーバによる推定法は推定時間が長ければ長いほど走行抵抗を考慮してないことによって推定誤差が広がってしまうが、提案した推定法は初期誤差がなければ走行抵抗を考慮しなくてもそれほど大きな誤差拡大は確認されず推定が行えていることが確認できた。また、原理的には初期誤差があっても正確なスリップ率の推定を可能にする走行抵抗を考慮したスリップ率推定方法を実現できた。
【図面の簡単な説明】
【０１０５】
【図１】車両の各部に働く力を模式的に示す図である。
【図２】タイヤと路面間の摩擦係数μとスリップ率λ両者の関係を表すμ−λ曲線を示す図である。
【図３】左右２輪駆動の車両モデルブロック図である。
【図４】駆動力オブザーバブロック図である。
【図５】走行抵抗推定器ブロック図である。
【図６】（ａ）は車体速度の真値と推定値を示す図であり、（ｂ）はＤＦＯに基づくスリップ率の真値と推定値を示す図であり、（ｃ）は初期誤差無しのＳＲＥに基づくスリップ率の真値と推定値を示す図であり、（ｄ）は初期誤差有りのＳＲＥに基づくスリップ率の真値と推定値を示す図であり、（ｅ）は初期誤差無しのＤＲＥ−ＳＲＥに基づくスリップ率の真値と推定値を示す図であり、（ｅ）は初期誤差無しのＤＲＥ−ＳＲＥに基づくスリップ率の真値と推定値を示す図である。
【図７】（ａ）、（ｅ）は左車輪、右車輪のオフライン実験における車輪速度と車体速度を示す図であり、（ｂ）、（ｆ）は左車輪、右車輪のＤＦＯに基づくオフライン実験におけるスリップ率の真値と推定値を示す図であり、（ｃ）、（ｇ）は左車輪、右車輪の初期誤差有りのＳＲＥに基づくオフライン実験におけるスリップ率の真値と推定値を示す図であり、（ｄ）、（ｈ）は左車輪、右車輪の初期誤差無しのＳＲＥに基づくオフライン実験におけるスリップ率の真値と推定値を示す図である。
【図８】左右輪のトルク指令値を示す図である。
【図９】駆動力積分による車体速度の推定値を示す図である。
【図１０】車輪速度制御ブロック図である。
【図１１】（ａ）は車輪速度制御を行った場合のスリップ率の真値と推定値を示す図であり、（ｂ）は車輪速度制御を行った場合のトルクを示す図であり、（ｃ）は車輪速度制御有りの車輪速度と車体速度を示す図であり、（ｄ）は車輪速度制御無しの車輪速度と車体速度を示す図である。
【図１２】本発明の一実施形態に係るスリップ率推定装置を示すブロック図である。
【図１】
0004936552.tif 000039
【図２】
0004936552.tif 000040
【図３】
0004936552.tif 000041
【図４】
0004936552.tif 000042
【図５】
0004936552.tif 000043
【図６】
0004936552.tif 000044
【図７】
0004936552.tif 000045
【図８】
0004936552.tif 000046
【図９】
0004936552.tif 000047
【図１０】
0004936552.tif 000048
【図１１】
0004936552.tif 000049
【図１２】
0004936552.tif 000050