(57)【要約】【課題】予定された経路、特に信号交差点を通過しようとする自動運転車両が他の車両などとの位置関係を把握し安全を判断して走行する方法を検討するには、自車の走行占有領域の検討、及び他車両等との接触を回避して走行する制御則検討に供する道具が必要である。【解決手段】  軌跡設計走行制御検証方法を実行する装置は、行程域の地図を作成する地図作成手段、目標軌跡を幾何学的に且つ動力学的に辿る自車両外形の移動画を描く描画手段、行程域の他車両等の目標軌跡を目的的に作成できその軌跡を幾何学的に辿る他車両外形の移動画を描く描画手段、緯度経度座標上で、自車両と他車両などとの相対位置変化を運動力学的、視覚的に確認する確認手段を含む。【選択図】  図1
2020119060.tif 000002
【特許請求の範囲】
【請求項１】
運行経路上の、或る場所の或る車線を決めてその或る場所に入り出るまでの行程において、検討対象の自動運転車両が移動する際の占有空間が刻々と描けて、自車両に備える周辺認識装置によって周辺交通事象を認識して、その行程域を通過する過程の制御則を見出す検討に使用する自動運転車両の軌跡設計走行制御検証方法。
【請求項２】
請求項１に記載の自動運転車両の軌跡設計走行制御検証方法において、道路線形の変化点座標を取得して、緯度経度のXY座標へ変換して、各変化点間を補完し、コーステーブルを作成する自動運転車両の軌跡設計走行制御検証方法。
【請求項３】
請求項１に記載の自動運転車両の軌跡設計走行制御検証方法において、横方向のつり合い式と回転のつり合い式を備える車両運動モデル式を用いるシミュレーションにより、極低速時から車速上昇時の車両旋回姿勢変化を算出して走行占有領域の検討と通常の車両移動解析を可能にする自動運転車両の軌跡設計走行制御検証方法。
【請求項４】
請求項１に記載の自動運転車両の軌跡設計走行制御検証方法において、自車の計画車速、他車の計画車速を設定し、その車速が信号現示時間に対応するようにして、信号時間の余裕検討を可能にする自動運転車両の軌跡設計走行制御検証方法。
【請求項５】
請求項１に記載の自動運転車両の軌跡設計走行制御検証方法において、右折ギャップを時間軸と距離軸で表現して、他車との接触回避距離と接触回避時間を算出して、自車速度の増減検討ができるようにする自動運転車両の軌跡設計走行制御検証方法。
【請求項６】
請求項１に記載の自動運転車両の軌跡設計走行制御検証方法において、自車の重心点を原点とする極座標で他車位置を捉え、自車の巾方向、前後方向の他車干渉余裕代を設けて自車の速度増減検討ができる自動運転車両の軌跡設計走行制御検証方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自動運転車両が狭隘路や交通混雑を通り抜ける際の走行占有領域の検討、及び他車両等との接触を回避して走行する制御則検討に供する自動運転車両の軌跡設計走行制御検証方法に関する。
【背景技術】
【０００２】
交通流シミュレータの報告がある。事前に予測することが難しい「イベントの開催、道路工事等のための交通管制／規制の影響、ETC等の導入効果など」を予測・評価するため、或いは「信号の変化や車両の合流、分岐、歩行者の横断といった道路上で起こる事象を、１台ごとの車両に反映して、数千台の車両からなる渋滞などを予測するためのものである。
【０００３】
交通渋滞のような交通現象のみでなく交通事故などのシミュレーションも行える予防安全シミュレータの報告がある。道路設計のために大型連結車両等の走行軌跡空間を確認するための建築設計のためのCADがある。
【０００４】
自動運転車両の開発が進行している。予定された経路上、特に交差点で生じる事象に対応して安全かつ円滑に走行する方法を理論的視覚的に検討する計算解析ツールの提供が望まれる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８−１２９８７４号公報
【特許文献２】特開２０１１−８６１１０号公報
【特許文献３】特開２００８−２０３９１３号公報
【非特許文献】
【０００６】
【非特許文献１】自動車技術会交通事故予測シミュレーション検証マニュアル
【非特許文献２】建築設計製図CAD事例
【発明の概要】
【発明が解決しようとする課題】
【０００７】
予定された経路、特に信号交差点を通過しようとする自動運転車両が他の車両などとの位置関係を把握し安全を判断して走行する方法を検討するには、その自動運転車両が交差点等で遭遇する事象を洗い出す手法としての道具と、洗い出された事象に対処する装備と制御についてシミュレーションを実施して評価・確認する道具が必要である。
【０００８】
上述の先行技術文献が参考になるも、特許文献１はシミュレーション対象通路のレーンを決めて、自動車、自転車または歩行者の刻々と移動する先の位置を計算して、複数種別の複数の移動体が移動する位置の時系列変化をシミュレーションする装置であり、自動運転車両の走行占有領域の検討や他車両などとの接触等を回避する制御の在り方の検討に供するものではない。
【０００９】
特許文献２は、運転支援装置に対するドライバの反応を考慮して交通流をシミュレーションするもので、自動運転車両の走行占有領域の検討や他車両などとの接触等を回避する制御の在り方の検討に供するものではない。
【００１０】
特許論文３は、仮想ドライバの知覚ミス、認識ミス、判断ミスによる交通事故率をシミュレーションする装置であって、自動運転車両の走行占有領域の検討や他車両などとの接触等を回避する制御の在り方の検討に供するものではない。
【００１１】
非特許文献１は、交通事故予測シミュレータが備えるべき仕様についてのマニュアルであり、自動運転車両の走行占有領域の検討や他車両などとの接触等を回避する制御の在り方の検討に供するものではない。
【００１２】
非特許文献２は、大型車両や大型連結車両が通り抜けるために必要な道路空間の検討等に用いるCADであり、自動運転車両の走行占有領域の検討の参考にはなるが、他車両などとの接触等を回避する制御の在り方の検討に供するものではない。
【課題を解決するための手段】
【００１３】
本発明は、旅程（トリップスケジュール）上の、或る場所の或る車線を決めてその或る場所に入り出るまでの行程域において、検討対象の自動運転車両が移動する際の占有空間が刻々と描けて、自車両に備えるLiDAR等の周辺認識装置によって周辺交通事象を認識して、その行程域を通過する過程の制御則を見出す検討に使用する自動運転車両の軌跡設計走行制御検証方法を提供する。
【００１４】
即ち、本発明に係る軌跡設計走行制御検証方法は、運行経路上の、或る場所の或る車線を決めてその或る場所に入り出るまでの行程において、検討対象の自動運転車両が移動する際の占有空間が刻々と描けて、自車両に備える周辺認識装置によって周辺交通事象を認識して、その行程域を通過する過程の制御則を見出す検討に使用する構成である。
【００１５】
換言すれば本発明は、検討対象となる自動運転車両が所定の走行区間を走行する際の軌跡をコンピュータの画面に表示することで、走行軌跡の設計や走行制御の検証に資する軌跡設計走行制御検証方法であって、前記コンピュータの画面への表示は周辺認識装置によって認識した周辺交通事象を反映するとともに自動運転車両の移動につれた占有面積が表示される構成である。
【００１６】
軌跡設計走行制御検証方法を実行する装置は、行程域の地図を作成する地図作成手段、目標軌跡を幾何学的に且つ動力学的に辿る自車両外形の移動画を描く描画手段、行程域の他車両等の目標軌跡を目的的に作成できその軌跡を幾何学的に辿る他車両外形の移動画を描く描画手段、緯度経度座標上で、自車両と他車両などとの相対位置変化を運動力学的、視覚的に確認する確認手段を含む。
【００１７】
自車両の重心位置を原点とする極座標上で、自車両と他車両などとの相対位置変化を運動力学的、視覚的に確認することで、接触回避のための他車に対する自車の制御要件の検討を行う。
【００１８】
前記したように、本発明に係る軌跡設計走行制御検証方法は、軌跡設計の検証及び走行制御の検証のうち一方のみの場合を含み、また検証には設計した軌跡により車両が狭い空間を通過できるか或いは設計した軌跡（クロソイド曲線）によりスムーズな走行が可能かの検証が含まれ、走行制御の検証には設計した軌跡に対応した車両速度の検証が含まれる。
【発明の効果】
【００１９】
また本発明に係る「自動運転車両の軌跡設計制御検証方法」によれば、自車両の重心位置を原点とする極座標上で、自車両と他車両などとの相対位置関係変化の解析検討ができるので、自車両の周辺認識のため装備するカメラやレーザーなどに備えるべき仕様の設計検討に直接的に適用することができ、自動運転車両開発のための有効な開発ツールになる。
【図面の簡単な説明】
【００２０】
【図１】当計算解析ツールの概念図である。
【図２】当計算解析ツールの開発工程図である。
【図３】目標軌跡の作成要領である。
【図４】運動計算シムリンク図である。
【図５】車両旋回姿勢の説明図である。
【図６】車両旋回占有領域説明図である。
【図７】右折ギャップの捉え方の説明図である。
【図８】右折ギャップにおける自車速度制御の説明図である。
【図９】自車両重心位置を原点とする極座標上での干渉検討要領である。
【発明を実施するための形態】
【００２１】
本発明は、自動運転車両が狭隘路などを通り抜ける際の走行占有領域の検討、及び青信号や右折信号の合間に他の車両等との接触を回避して走行する制御方法を検討するための計算解析ツールに関するものである。自動運転車両の開発の現場の求めに直に答える計算解析ツールを提供するものであるので、現実の経路での運行計画に即して解析しての現実解が求められる。故に、現地現物に即しての有効な結果が得られるものでなくてはならない。
【００２２】
図１に、当計算解析ツールの概念図を示す。事例として、新宿５丁目東交差点の場合で示す。本発明が目的の一つとする“（１）走行占有領域の検討”の概念を図中（１）、（２）、（３）、（４）に示す。（１）は縁石に乗り上げない後輪軌跡、（２）は隣接車線に張り出さない軌跡、（３）は円滑操舵軌跡である。（２）の対策によって（４）の他車追突防止ができる。
【００２３】
本発明の目的の二つ目である“青信号や右折信号の合間に他の車両等との接触を回避して走行する制御方法を検討”の概念を図中（２）の（５）に示す。右折車の走行軌跡と直進車の走行軌跡が交差する点から、直進車の車頭までの距離時間を右折ギャップと称する。そのギャップを確保する制御則はどうあるべきか。そのギャップは、曲がった先の状況（３）の(６）も予測しての行動が求められる。
【００２４】
図１の右に自車の重心点を原点とする極座標を示す。この極座標で捉えた他車軌跡を、自車に搭載するジャイロ、加速度計から算出する自車の平面軌跡に重畳させると図１の左側部に示すように、他車座標を計算でき、右折ギャップの計算ができる。また、極座標上で自車の前後長、横巾に対しする接触余裕を確保する自車速度調整ができる。
【００２５】
図２に当計算解析ツールの開発工程図を示す。工程（１）で、対象交差点等の現場検証をする。現地に赴いて解析に必要な現地データを取得する。工程（２）で対象交差点等の地図を作成する。Google Earthを利用するか。国土地理院の地図を入手するかなどをする。現地に赴いて現場計測もする。工程（３）で目標軌跡を作成する。直線部は線形補間し、曲線部はクロソイド関数による曲線補間をすることによって自動運転に適した軌跡を設計する。工程（４）で自車・他車運動シムリンクを作成する。自車の移動軌跡、他社の移動軌跡、自車と他車の相対移動の解析ができる。自車軌跡は車両モデルを用いるも、他車についてはモデルを省略して点移動としても良い。工程（５）で、車両挙動と衝突回避の解析をする。自車移動領域・空間の成立確認、右折ギャップの余裕時間を求めて制御に用いる。
【００２６】
対象交差点の現場解析の例としては、運行ルートの始点から終点に至る自車両の任務（Mission）を分析して、解析因子とその条件を定める。具体的には、走行ルートを辿り、曲率、勾配、舗装非舗装、車線数…等の道路仕様、交通標識と設置場所、路肩利用や交差点マナーなどの利用特性、利用車両の分布等々を調べ車頭間隔や車速分布や信号現示時間などの実測値を取り上げる。
【００２７】
図３は、目標軌跡の作成要領である。（１）にて、始点、終点を含む経路線形の変化点の緯度経度を取得する。（２）にて緯度経度をＸＹ座標に変換する。この変換は、「平面１４年国土交通省告示第９号」で定義されている平面直角座標系（http://www.gsi.go.jp/sokuchikijun/jpc.html）を使用する。全国１９の平面直角座標系で構成されており、各座標系の原点がそれぞれ定められている。例えば、関東地方に福島県を含むIX(9)系の原点は、東経139度50分0秒、北緯36度0分0秒と定められている。（３）にて、各変化点間を直線ないし曲線補間し、それを（４）にて地図に重ねる乃至グラフにして補間状態を確認し、（５）のコーステーブルにする。このテーブルには、時間ではなく距離での同期をとるためのインデックス（ＩＤ）を付して、Ｘ座標、Ｙ座標、方位角、曲率、経路長、計画速度、停止線位置が記録され制御に供される。
【００２８】
図４は、運動計算シムリンクである。（Ａ）の部分が自車解析、（Ｂ）の部分が他車解析、（Ｃ）の部分が自車他車相対解析に関する。
【００２９】
(Ａ)の自車解析について説明する。図３の（５）のコーステーブルから計画車速（詳細（13）で説明）、目標軌跡曲率が（1）にて式（１）により前輪実舵角が算出される。
【００３０】
【数１】
2020119060.tif 000003
【００３１】
式（１）によって算出された実舵角δに修正実舵角δＦＢ（詳細（11）で説明）が加算され（2）の車両運動モデル式に計画車速と共に代入される。車両運動モデル式（詳細後述）は、車両の横運動のつり合い式と回転運動のつり合い式による連立微分方程式である。車両重心点の横速度を前後速度で除算した値である車両横すべり角（β）と、車両重心点まわりの回転角速度ヨーレイト（γ）が算出される。このヨーレイトを積分するとヨー角（φ）になる。
【００３２】
このヨー角（φ）と車両横すべり角（β）を（3）にて合算して、その余弦成分と正弦成分を得る。（4）において余弦成分に車速を乗じて積分すると前後移動距離（X1）が算出され、（5）において正弦成分に車速を乗じて積分すると横移動距離(Y1)が算出される。X1とY1の二乗和の平方根（6）から（7）の距離インデックスＩＤが得られる。このインデックスＩＤを、図３（５）のコーステーブルのＩＤに照らして、座標、方位角、曲率、経路長、計画車速、停止線位置を摘出して制御入力とする。
【００３３】
前後移動距離（X1）、横移動距離(Y1)、車体横すべり角（β）及びヨー角（φ）から、ＸＹ座標における自車位置と自車姿勢が（15）の様に描かれる。これを演算時間毎に差し替えるとアニメーション（動画）になる。
【００３４】
目標軌跡曲率に車速を（8）で乗じると道路線形から生じるヨーレイトになる。この道路線形ヨーレイトを車両ヨーレイトから差引き（9）にて積分すると道路線形に対する車両軌跡接線との角度偏差（e3）が得られる。（10）にて、横移動距離（Y1）から目標軌跡座標のY値を差引くと横偏差（e2）が得られる。
【００３５】
得られた横偏差（e2）と角度偏差（e3）を（11）の修正実舵角の式に代入すると修正実舵角（δＦＢ）が得られる。
【００３６】
【数２】
2020119060.tif 000004
【００３７】
式（２）のｅ３Ｖｔｃは制御サイクル時間ｔｃの間に生じる角度偏差得e3による横偏差である。それに、横遍差e2を加えた値に曲率（ρ＝ｖ/γ）を乗じて、修正操舵相当の曲率とし、その修正操舵相当曲率にｌ（エル）/ｌｒを乗じて前輪実舵角とする。ここにｌは、ホィールベース、ｌｒは、重心から後軸までの距離である。δFBは、ＰＩＤコントローラ（12）にて、積分ゲインkI、比例ゲインkp、微分ゲインｋDによって合わせ込まれ偏差e2,e3を最小化する。画して、車両モデルの軌跡は、目標軌跡に整合する。
【００３８】
目標軌跡には停止線位置が仕込まれている。現示信号、対向車線などの状況から通過可否判定がされ、通過可能との判定であれば計画車速で通過し、通過不可能との判定であれば停止線で停止する。その制御を（13）の現示信号時間の適正化検討も含めて速度制御部で実施する。
【００３９】
図４の（13）に信号１と信号２の信号現示例を示す。この場合、信号２が自車線方向、信号１が交差車線方向である。公共バスとしての車内転倒事故防止の観点から抑えられるべき減速度があり、貨物トラックとして荷痛み防止の観点から押さえられるべき減速度がある。その減速度を“Gain Block”に仕込み、その減速度で停止するための停止距離に対応する距離インデックスを“Constant Block”に仕込んで、信号現示と“AND Block”に取り込んで、信号に同期しての停止線停止ができる。
【００４０】
図４の(Ｂ)の他車解析について説明する。（Ａ）の自車には車両の運動を表現するモデル式を備えたけれども、他車については、目標軌跡を計画計画車速で移動して信号の切り替わりに対応すればよい。ここでは、一台のみを示すが、解析目的に応じて複数設定もする。他車の座標はXn、Ynとして出力される。ここに、サフィックスのｎは他車台数により付記される。
【００４１】
図４の(C)の他車解析について説明する。（A）の自車解析から自車座標X1,Y1、車体横すべり角β、ヨー角φが、（B）の他車解析から他車座標Xn,Ynが（C）へ入力されて、XY座標上にプロットされ、計算周期毎に移動する。自車の重心点に原点を置き自車前後方向をｘ軸、横方向をｙ軸とする移動座標（ｘｙ座標）が描かれる。
【００４２】
自車位置（X1,Y1）と他車位置（Xn,Yn）の相対位置関係は式（３）、式（４）、式（５）になる。
【００４３】
【数３】
2020119060.tif 000005
【００４４】
【数４】
2020119060.tif 000006
【００４５】
【数５】
2020119060.tif 000007
【００４６】
かくして、自車に対する他車の位置は、自車他車間相対距離（ｌ(エル）と相対角度（θ）とする極座標(x,y)で表現され、自車に搭載されるカメラ、レーザー等の周辺認識装置との対応がとれる。
【００４７】
図５は車両旋回姿勢の説明図である。（Ａ）に車両運動モデル図、（Ｂ）に極低速時の旋回姿勢図、（Ｃ）に車速上昇時の旋回姿勢図を示す。（Ａ）について説明する。ここでは、前1軸、後1軸の前後2軸車の場合について示す。他に、前1軸、後2軸の3軸車、前2軸、後1軸の3軸車、前2軸、後2軸の4軸車も備える。
【００４８】
車両の前後軸（x軸）上に左右輪を合わせ置いた前1輪・後1輪の自転車モデルとする。速度一定と仮定すると、このモデル横方向のつり合い式（6）、回転のつり合い式（7）、補助方程式(8),(9),(10),(11)で表現され、補助方程式を式（6）、式（7）に代入すると式(12)、式(13)になる。
【００４９】
【数６】
2020119060.tif 000008
【００５０】
【数７】
2020119060.tif 000009
【００５１】
【数８】
2020119060.tif 000010
【００５２】
【数９】
2020119060.tif 000011
【００５３】
【数１０】
2020119060.tif 000012
【００５４】
【数１１】
2020119060.tif 000013
【００５５】
【数１２】
2020119060.tif 000014
【００５６】
【数１３】
2020119060.tif 000015
【００５７】
ここに、記号の説明は次の通りである。
Ｍ：車両質量、Ｉ：車両慣性モーメント、ｌ：ホィールベース、ｌｆ：重心点から前軸までの距離、ｌｒ：重心点から前軸までの距離、Ｖ：重心点速度、Ｖｆ：前軸速度、Ｖｒ：後軸速度、δ：前輪実舵角、β：車体横すべり角、βf：前輪横すべり角、βr：後輪横すべり角、ψ’（＝ω）:ヨーレイト、ＣＦ：前軸のコーナリングフォース、ＣＲ：後軸のコーナリングフォース、Ｋｆ：前軸コーナリング係数。Ｋｒ：後軸コーナリング係数
【００５８】
式（12）、式（13）は、タイヤ横すべりがある図６の（ｃ）の状態を表現する式である。この式を次の様に展開すると、図５の（Ｂ）の極低速（≒速度ゼロ）の状態を表現する式になる。
【００５９】
即ち、式（１２）は、式（１２‐１）から式（１２−４）に展開され、式（１３）は、式（１３‐１）から式（１３−６）に展開されて、図５の（Ｂ）を表現する式になる。
式（１２）を展開する。
【００６０】
【数１２−１】
2020119060.tif 000016
【００６１】
【数１２−２】
2020119060.tif 000017
【００６２】
【数１２−３】
2020119060.tif 000018
【００６３】
【数１２−４】
2020119060.tif 000019
【００６４】
【数１３−１】
2020119060.tif 000020
【００６５】
【数１３−２】
2020119060.tif 000021
【００６６】
【数１３−３】
2020119060.tif 000022
【００６７】
【数１３−４】
2020119060.tif 000023
【００６８】
【数１３−５】
2020119060.tif 000024
【００６９】
【数１３−６】
2020119060.tif 000025
【００７０】
図６は、車両旋回占有領域説明図である。式（１3−６）で旋回すると、車両後部に対して、車両前部は外側に張り出して旋回する。刻々と変化する車両位置を動画として示すことで、車両旋回占有領域を視覚的に確認できる。
【００７１】
図７は、右折ギャップの捉え方の説明図である。自車が右折する軌跡と対向車線を直進してくる他車との交点ＰGと直進他車との距離を右折ギャップと称する。自車からＰu点までの距離ｌ１と角度αＵ及び自車からＰＧ点までの距離ｌ２と角度αＧとすると、右折ギャップｌＧは、式（14）になる。
【００７２】
【数１４】
2020119060.tif 000026
【００７３】
ｌ１の微分値即ち距離変化速度、ｌ2の微分即ち距離変化速度、αＧの微分値即ち角度変化速度、αＵの微分値即ち角度変化速度から、ｌＧの距離変化速度が分かる。即ち、ｌ2の微分値は即ち自車速度であり、自明であるので、ｌ１、αＧ、αＵをライダー等で把握することにより、右折ギャップの変化速度、即ち、直進他車の車速が分かる。
【００７４】
図８は、右折ギャップにおける自車速度の説明図である。式(14)に対応するシムリンクである。(1)にてｌ１を微分してαＵの余弦を乗じた値から、（2）のｌ２を微分してαＧの余弦を乗じた自車速度（3）を引くと（4）の他車速度ＶｐＯになり、これを積分して（5）の元のｌＧから差引いた距離（6）は、（7）の様に時間変化する線になる。一方、自車速度（3）を（8）の元のｌ２から差引いた距離（9）は、（10）の２本の点線の様に時間変化する。２本の点線の上側は自車が減速する場合であり、上側は増速する場合である。画して、（11）が他車に対する自車の接触回避距離であり、（12）が他車に対する自車の接触回避時間である。この回避距離、回避時間は自車の速度制御で拡大することが出来る。
【００７５】
図９は、自車両重心位置を原点とする極座標上での干渉検討要領である。自車の重心を原点とする極座標である。（1）の半径rの右手前方に相手車両が存在する場合を例示している。式(14)は相手の横距離の式である。（1）のβは自車の車体横すべり角β、(2)のθは相手の存在方向角、（3）のγはヨーレイトで、それを成分しているからヨー角である（1）+（2）+（3）余弦に車速を乗じ積分した値は相手方向の自車の速度ベクトルである。それをｒから差引いた距離に方向角θの正弦を乗じ自車に対する相手の横位置を求め、これが接触余裕代を含む巾Bの２分の１以上であれば、接触回避できる。
【００７６】
式(15)は相手の前後距離の式である。式(15)の（5）の正弦部分を（7）の余弦に、（6）の巾Bの２分の一を前後長さLfに置き換えた式である。相手の存在位置がθ＜±π/２の場合は自車速度を下げる方向が接触回避方向であり、θ＞±π/２の場合は自車速度を上げる方向が接触回避方向である。
【図１】
2020119060.tif 000027
【図２】
2020119060.tif 000028
【図３】
2020119060.tif 000029
【図４】
2020119060.tif 000030
【図５】
2020119060.tif 000031
【図６】
2020119060.tif 000032
【図７】
2020119060.tif 000033
【図８】
2020119060.tif 000034
【図９】
2020119060.tif 000035