JP3763476B2 - 車両及び運転者の挙動解析システム - Google Patents

車両及び運転者の挙動解析システム Download PDF

Info

Publication number
JP3763476B2
JP3763476B2 JP2003153184A JP2003153184A JP3763476B2 JP 3763476 B2 JP3763476 B2 JP 3763476B2 JP 2003153184 A JP2003153184 A JP 2003153184A JP 2003153184 A JP2003153184 A JP 2003153184A JP 3763476 B2 JP3763476 B2 JP 3763476B2
Authority
JP
Japan
Prior art keywords
vehicle
model
data
driver
assessment
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.)
Expired - Fee Related
Application number
JP2003153184A
Other languages
English (en)
Other versions
JP2004352116A (ja
Inventor
ヒーマー マーカス
隆徳 松永
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2003153184A priority Critical patent/JP3763476B2/ja
Priority to US10/680,219 priority patent/US7239945B2/en
Priority to DE10353650A priority patent/DE10353650A1/de
Publication of JP2004352116A publication Critical patent/JP2004352116A/ja
Application granted granted Critical
Publication of JP3763476B2 publication Critical patent/JP3763476B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/101Side slip angle of tyre
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/02Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to ambient conditions
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/08Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to drivers or passengers
    • B60W40/09Driving style or behaviour
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/105Speed
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/12Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to parameters of the vehicle itself, e.g. tyre models
    • B60W40/13Load or weight

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
  • Time Recorders, Dirve Recorders, Access Control (AREA)

Description

【0001】
【発明の属する技術分野】
この発明は、一般的には事故再現システムに関し、特に、危機的な運転状況を解析および査定するシステムに関するものである。
【0002】
【従来の技術】
今日、最新の車両において、乗員の安全性を高めるために多大な努力がなされている。エアバッグシステムおよび車体構造の改良に加えて、制御システムの果たす役割は大きくなってきている。ブレーキの電子制御システム、すなわち車両安定制御(VDC:vehicle dynamics control)は、最新の車両において搭乗者の安全性を著しく高めた。なお、関連の先行技術文献として、非特許文献1、2、3、4及び5がある。
【0003】
VDCシステムの普及に伴い、事故の数秒前の事象を再現することが難しくなるという弊害がある。路面のブレーキ痕が見当たらないことは非常に多い。ブレーキ痕は事故再現の専門家が調査する際の重要な手段であるため、犯罪であるか否かの問題に明確に答えることができないのである。
【0004】
さらに、保険金詐欺の数は増加傾向にある。今日では、ドイツの保険金がらみの全事象のうち10%は人為的操作によるものであり、ドイツの年間財政損失は10億ドルを超えると推定される。
【0005】
事故再現のための既存のシステム、すなわち、いわゆるイベントデータレコーダ(EDR:event data recorder)は、専用に搭載されたセンサによって得られる測定情報のみを用いる(例えば、特許文献1参照)か、または車両のエアバッグシステムを用いている。
【0006】
別のシステムは、一組の微分方程式を解くことによって車両の動きを再構築する(例えば、特許文献2参照)。
【0007】
【特許文献1】
独国特許発明第19509711号明細書(第1頁−第3頁)
【特許文献2】
独国特許発明第4132981号明細書(第1頁−第4頁、図1−図3)
【非特許文献1】
M. Burckhardt. 「Radschlupf-Regelsysteme」 Vogel Fachbuchgruppe: Fahrwerktechnik. Vogel Buchverlag, Wuerzburg, 1st edition, 1993
【非特許文献2】
S. L. Miller et al. 「Calculating longitudinal and lateral wheel slip and tire parameters using GPS velocity」 Proceedings of the American Control Conference, Society of Automotive Engineers (SAE), Arlington, VA, USA, 2001
【非特許文献3】
U. Kiencke and L. Nielsen. 「Automotive Control Systems」 Springer-Verlag, Berlin, Heidelberg, New York, 2000
【非特許文献4】
M. Mitschke.「Dynamik der Kraftfahrzeuge」 Band A: Antrieb und Bremsen. Springer-Verlag, Berlin, Heidelberg, New York, 2nd edition, 1988
【非特許文献5】
M. Mitschke. 「Dynamik der Kraftfahrzeuge」 Band C: Fahrverhalten. Springer-Verlag, Berlin, Heidelberg, New York, 2nd edition, 1990
【0008】
【発明が解決しようとする課題】
上述したような従来のシステムでは、事故等の危険な運転状況を定量的に解析し、評価するために必要なすべての情報を記録することができないという問題点があった。
【0009】
この発明は、上述のような課題を解決するためになされたもので、第1の目的は、事故等の危険な運転状況を完全に解析するために必要な情報を記録することができる車両及び運転者の挙動解析システムを得るものである。
【0010】
また、第2の目的は、主に事故再現のために使用するが、それに加えて、交通規則違反の証明、または必ずしも事故に至らなかった不安定な運転状況の解析の手段としても使用することができる車両及び運転者の挙動解析システムを得るものである。
【0011】
【課題を解決するための手段】
この発明に係る車両及び運転者の挙動解析システムは、車載電気機器の動作情報、車載データネットワークを介して得られるデータ、及び車両を構成する電気装置から得られる付加情報を内部で演算可能な形式に変換を行うデータ前処理手段と、前記データ前処理手段からの前処理済データに基づいて車両搭載センサでは測定できない変数を推定する車両モデル/推定手段と、前記データ前処理手段からの前処理済データ、及び前記車両モデル/推定手段からの推定された変数を共通の時刻管理情報に統合して事故のシナリオを再現するモデルベース再現手段と、前記モデルベース再現手段により再現された事故のシナリオに基づいて車両、及び運転者の挙動を査定する査定手段とを設け、前記車両モデル/推定手段は、前記データ前処理手段からの前処理済データに基づいて、車両に搭載されるセンサでは測定できない変数として、車両の重心速度を求める車速推定部と、前記データ前処理手段からの前処理済データに基づいて、車両の動力伝達機構を数式として構築した駆動モデルを用いて駆動トルク求める駆動モデル部と、前記前処理済データ、前記車両の重心速度、及び前記駆動トルクに基づいて、車両に搭載されるセンサでは測定できない変数として、車両の重量を推定する車重推定部と、前記前処理済データ、及び前記車両の重心速度に基づいて、車輪を数式として構築した車輪モデルを用い、車両に搭載されるセンサでは測定できない変数として、タイヤの横すべり角と車輪に作用する力を計算する車輪モデル部と、前記前処理済データ、前記車両の重心速度、前記車両の重量、及び前記車輪に作用する力に基づいて、走行時における車両の状態を数式として構築した状態空間モデルを用い、車両に搭載されるセンサでは測定できない変数として、車体の横すべり角を導出する状態空間モデル部とを含むものである。
【0012】
【発明の実施の形態】
実施の形態1.
この発明の実施の形態1に係る車両及び運転者の挙動解析システムについて図面を参照しながら説明する。図1は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの構成を示す図である。なお、各図中、同一符号は同一又は相当部分を示す。
【0013】
図1において、車両モデル/推定装置(車両モデル/推定手段)30および査定装置(査定手段)70は、本システム10の中核を構成する。この車両モデル/推定装置30において実行されるアルゴリズムによって、事故再現のために重要な測定不可能な変数を求めることができる。開示される本システム10では、微分方程式の解が得られるだけでなく、さらに、推定手法を用いて状態および車両の重量等のパラメータも得られる。たとえば、車両の重量が推定されると、これによりその車両が過荷重であるか否かを判定することができる。さらに、摩擦係数を求めて、車両の周辺環境についての情報を導出することができる。またさらに、車両の対地速度および車体の横すべり角により、危機的な運転状況における車両状態についての情報が得られる。
【0014】
車両及び運転者の挙動解析システム10は、以下の2つの理由により、常に測定誤差(雑音、オフセット、ドリフト等)を受ける測定可能なセンサ情報のみを用いるシステムよりも、より正確なデータを得ることができる。
【0015】
第1に、本システム10は、既存の冗長性も用いる。たとえば、車輪速と加速度信号との組み合わせを用いることにより、前後方向の車両の移動量の測定を改善することができる。
【0016】
第2に、システム10内で使用される車両モデルは、(モデル自体に隠された事前の知識により)さらにフィルタのように動作する。これは、フィルタリング効果により、より高い精度を達成できることを意味する。
【0017】
これらの理由に関して、システム10により、事故のシナリオのより正確な再現が保証される。しかしながら、車両及び運転者の挙動解析システム10は、単に車両の挙動再現という分野における現状の「技術状態」を改善するだけではない。車両モデル/推定装置30およびモデルベース再現装置(モデルベース再現手段)50からのデータに基づいて、危機的な運転状況における車両および運転者の挙動を、さらに査定することができる。この新たな過程は査定装置70において実行され、この過程は、事故の数秒前を明らかにすることを助け、既存のEDRシステムを改良する。これは、既存のEDRシステムが車両の移動を描写することに限られているからである。
【0018】
危機的な運転状況における車両および運転者の挙動の解析に対し、重要な作用がいくつかある。
【0019】
まず、車両ダイナミクスおよびその軌跡を正確に描写する必要がある。それゆえ、事故再現のために重要ないくつかの測定不可能な変数を求めなければならない。危機的な運転状況における運転者の挙動について推論することもまた必要である。運転者が事故前に適切な行動をとらないことは非常によくある。したがって、運転者の入力を分析すること、およびそれにより得られる車両ダイナミクスにより運転者の挙動を査定することは、本システム10の基本的な要件である。
【0020】
最後に、車両の周辺環境は、道路交通における車両および運転者の挙動に著しく影響を及ぼす。したがって、事故前後の車両の周囲の現場の状態について推論する必要がある。
【0021】
つぎに、この実施の形態1に係る車両及び運転者の挙動解析システムの詳細な動作について図面を参照しながら説明する。
【0022】
本システム10は、車両において入手可能なデータを取得する。前照灯(ヘッドライト)、尾灯(テールライト)、ブレーキ、フォッグランプ、後退灯(back light)、方向指示器、ハザードランプ、ハイビームの状態が、電気機器2によって提供される。しかしながら、ワイヤハーネス12から利用可能であるのはライトの情報だけではない。さらに、チェックランプ(check lamp)やワイパー、クラクションの状態が記録され得る。緊急車両の場合には、非常灯と救急サイレン(ambulance horn)の状態が重要であるかもしれない。
【0023】
信号は、通常、マイクロプロセッサに印加できる最大電圧よりも高い電圧を有するため、ワイヤハーネス12から信号を受信した後に、データ前処理装置(データ前処理手段)20において事前に前処理する必要がある。
【0024】
今日の最新車両には、1つまたは複数の車載データネットワーク4が装備されている。多くの場合、そのような通信ネットワークは、CAN(Controller Area Network:コントローラエリアネットワーク)に代表される。
【0025】
本システム10は、利用可能なセンサ情報をCANバス14からデータ前処理装置20に読み込む。しかしながら、本システム10は、他のいかなる車載データネットワークからのデータも読み取ることが可能であり、CANバス14に限られない。また、本システム10は、車両のエンジン制御装置(ECU:Engine Control Unit)から直接データを取得することができる。
【0026】
最後に、本システム10は、例えばエアバッグシステムやABSシステムのような車両を構成する電気装置から、付加情報6が個々の導線16を介して取得する。この情報の全てが、データ前処理装置20に入力され、ここで信号がCPU内で処理されるために準備される。
【0027】
前処理後に、データは、共通の導線25を介して車両モデル/推定装置30へと伝達される。この車両モデル/推定装置30において、前処理済データは、測定することができない状態を推定するためと、事故再現理論で関係する車両に作用する力および他の変数を計算するために用いられる。さらに、導線25上の前処理済データは、車両のセンサ情報の既存の冗長性を用いて補強される。
【0028】
車両モデル/推定装置30において車両モデルと推定の計算を実行した後、データは内部双方向バス35に書き込まれ、メモリ40に保存される。このメモリ40としては、たとえばFLASH−ROMまたはEEPROMなどのような不揮発性メモリを用いても、あるいは揮発性ランダムアクセスメモリ(RAM)を用いてもよい。後者の場合、メモリモジュールはデータ損失からシステムを守るように、バッテリで保護される。事故が起こると、車両電源が遮断されるため、車両及び運転者の挙動解析システム10は、記録されたデータを失ってしまう可能性があるためである。
【0029】
事故の発生後、データはインタフェース45を介してモデルベース再現装置50へ伝達され、そこで事故後の車両位置から後戻りするためのバックトラッキングアルゴリズムが実行される。さらに、車両モデル/推定装置30からの全情報が、共通の時刻管理情報に従い、順序立てられた形で設定される。
【0030】
さらに、車両のワイヤハーネス12から収集された情報は、運転者が与える信号を解析するために用いられる。これら信号は、データ前処理装置20から複数の導線27を介してメモリ40へ直接流れる。メモリ40に情報を適切に保存する動作は、マイクロプロセッサ(図1には図示せず)によって管理される。
【0031】
事故についての人間の付加的な知識(人による情報)8が、ユーザインタフェース18によってモデルベース再現装置50に入力される。事故再現の専門家にとって興味深い、人による付加的な知識とは、たとえば、警官または重要な目撃者により提供される現場気象状況、道路状態(凍結、泥砂の散乱など)、カーブ半径、または道路傾斜のような車両の周囲の環境である。集まった情報に基づいて、事故のシナリオがモデルベース再現装置50において再現され、内部バス65を介して査定装置70へ流れる。
【0032】
この査定装置70において、事故のシナリオについての全情報が、車両の周囲の環境条件に左右される運転者および車両の挙動を査定するために処理される。
【0033】
図2は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムのデータ前処理装置の構成を示す図である。
【0034】
電気機器2からの「複数の静的状態」についての情報(ライト、ワイパー、クラクション等)は、ワイヤハーネス12を介してデータ前処理装置20に伝達される。ワイヤハーネス12からの信号レベルは、通常、車両のバッテリ電源に等しい。この信号レベルは、マイクロプロセッサ回路には高すぎるため、一般的な導線215を介してCPU(図2には図示せず)に転送するために、第1の信号変換部210において変換する必要がある。
【0035】
CANバス14、任意の他の車載データネットワーク、またはECUから得られた情報は、信号復号器220において復号化されなければならない。この信号復号器220は、供給されたデータの供給元と優先順位を得るために、CAN識別子を解読する。
【0036】
信号復号器220での解読後、信号は導線225を介してデジタルフィルタ240に伝達され、このデジタルフィルタ240は、高周波ノイズを除去するために解読したデータをローパスフィルタリングする。このデジタルフィルタ240は、必要な処理速度に応じてデジタル信号プロセッサ(DSP:Digital Signal Processer)またはマイクロコントローラで実施することができる。このデジタルフィルタ240は、内部バス245によってCPU(図2には図示せず)に接続される。
【0037】
エアバッグシステムのような車両サブシステムからの付加情報6は、一組の導線16を介して第2の信号変換部230に入力される。また、第2の信号変換部230は、CPUで処理するために信号を変換する。さらに、一組の導線16からの信号がアナログ信号の場合には、信号レベルはアナログ/デジタル変換器(ADC:Analog to Digital Converter)260のマージンの利用を最適にするように変換される。この第2の信号変換部230により変換された信号は、導線235を介してフィルタ250に入力され、エイリアス除去フィルタリングされる。このフィルタリングされた信号が導電線255を介して、ADC260に送られ、ここで信号はデジタル化される。デジタル化された信号は、ADC260から導線265を介してCPUに転送される。一方、導線275を介してCPUに直接伝達される信号は、ABS動作状態判別信号のように静的状態であり、フィルタリングもAD変換もする必要がない。
【0038】
図3は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の構成を示す図である。
【0039】
変換およびフィルタリングされた車載データネットワークからの信号が、内部バス245を介して車両モデル/推定装置30に供給される。この車両モデル/推定装置30は、1つまたは複数のマイクロプロセッサ上で実行されるアルゴリズムで動作する。結果を他のサブシステムで処理できるようなより良い概観を提供するために、図3に示す種々のサブシステム(310、315、320、325、330、335、340、345、350、360)が矢印によって接続される。これらの矢印は、車両モデル/推定装置30の構成要素(サブシステム)間の信号の流れを示す。実際には、これらのサブシステムは、レジスタおよび/または外部メモリから必要な入力を得て、かつ、それらの計算結果をレジスタおよび/または1つまたは複数のマイクロプロセッサ(図3には図示せず)によってアドレス指定された外部メモリに保存する。
【0040】
車速推定器310は、導線312を介して車両の重心(CoG:center of gravity)における対地速度を出力するために、内部バス245からの車輪速および前後方向(longitudinal)の加速度信号を用いる。これは1つのブロックで示され、車両における冗長情報を用いてより正確なデータを処理する。この車速推定器310は、図5でより詳細に説明する。
【0041】
垂直方向の力サブシステム315は、車両の重心(CoG:center of gravity)に近接して取り付けられた加速度センサからの前後方向(longitudinal)および横方向の加速度により、垂直方向の車輪の力を求めるアルゴリズムを含む。車輪に垂直方向に作用する力は、転がり摩擦サブシステム330、車輪モデルサブシステム340、及び摩擦推定器350の入力として必要であり、一組の導線316を介してこれらのサブシステムへ供給される。垂直方向の車輪の力に関する方程式の詳細な導出は、非特許文献3において見ることができる。
【0042】
駆動モデルサブシステム320は、内部バス245を介して提供されるエンジン回転数およびスロットル開度を処理し、車重推定器345において必要とされる被駆動車輪にかかる駆動トルクを、導線321を介して出力する。この駆動モデルサブシステム320は、図7により詳細に示される。
【0043】
ステアリングモデルサブシステム325は、運転者のステアリング入力(ハンドル角)と、導線326に供給される出力、すなわち車輪の操舵角との関係を示すサブシステムである。ステアリングモデルサブシステム325は、入力角と出力角との間の伝達率(transmission factor)を含む。このサブシステムをより正確に記述するためには、ハンドルの固さおよび車輪に横方向に作用する力も同様に考慮し得る。このステアリングモデルサブシステム325の詳細な記述は、基礎となる方程式の構成とともに、非特許文献5において見ることができる。
【0044】
転がり摩擦サブシステム330において、垂直方向の力サブシステム315において計算される車輪に垂直方向に作用する力、ならびに車速推定器310から推定されるCoG速度は、転がり摩擦抵抗を得るために処理され、導線331を介して車重推定器345に供給される。この転がり摩擦サブシステム330において転がり抵抗を求める基礎となる方程式は、非特許文献4に記載される近似方程式である。これは、車速、車輪に垂直方向に作用する力、および経験的実験において求めたいくつかの固有の転がり抵抗係数を含む四次多項式である。
【0045】
道路勾配および風力サブシステム335は、道路勾配および風力を求めるユニットである。これら両方の影響は、力の平衡を設定する際に考慮されなければならない。道路勾配は、車重推定器345へ、かつ導線336を介して非線形状態空間モデルサブシステム360へ転送される。導線337は、道路勾配および風力サブシステム335において計算された風力を車重推定器345と非線形状態空間モデルサブシステム360に供給する。その結果、道路勾配は、内部バス245を介して得られる車輪速および前後方向の加速度によって処理される。前後方向の加速度と、車輪速の導関数(すなわち車輪加速度)との間に導出される結果は、道路勾配によって発生する加速度センサの誤った向きによるものである。風速は、導線312を介して供給される力と車両のCoG速度と力の二次関係式によって計算される。
【0046】
車両モデル/推定装置30の重要なサブシステムの1つは、車輪モデルサブシステム340である。2組の導線341および342を介して出力される前後方向および横方向の車輪の力に加えて、タイヤの横すべり角(TSSA:tire side slip angle)が計算される。このタイヤの横すべり角は、事故再現において重要な測定不可能な基本変数である。このタイヤの横すべり角を解析することで、オーバーステアまたはアンダーステアの運転状況がわかる。タイヤの横すべり角が車輪モデルサブシステム340において求められ、一組の導線343を介して出力されるのはそのためである。この車輪モデルサブシステム340は、内部バス245を介して得ることができる前処理済ヨーレート信号、導線312を介して得られるCoG速度、一組の導線316からの車輪に垂直に作用する力、導線326を介して得られる車輪の操舵角、及び導線371を介して得られる車体横すべり角(VBSSA:vehicle body side slip angle)を用いる。この車体横すべり角は、非線形状態空間モデルサブシステム360の非線形観測器370において計算される。この車輪モデルサブシステム340は、図6において詳細に示される。
【0047】
ある車両が過荷重であるか否かを推論するためには、実際の車重を知ることが望ましい。車重推定器345において再帰的最小二乗(RLS:recursive least squares)推定が実行され、導線346に車重を出力するのはそのためである。このRLS推定法は、リアルタイムで変数を推定する一般的な手段であり、そのため、同定および推定理論に精通した者にはよく知られている。計算による可能性に応じて、基本的には操作変数法や最尤法などの他のパラメータ推定法も適用することができる。車重推定器345において実施される質量推定方程式は、CoGにおける力の平衡から導出される。駆動モデルサブシステム320によって導線321を介して供給される駆動力に加えて、抵抗力も処理されるのはそのためである。後者は、導線331を介して(転がり摩擦サブシステム330から)、ならびに導線336および337を介して(道路勾配および風力サブシステム335から)得ることができる。さらに、導線312を介するCoG速度が必要である。基礎となる力の平衡ならびにRLS法は、非特許文献3に記載される。
【0048】
摩擦推定器350は、垂直方向の力サブシステム315において計算され、導線316を介して供給される車輪に垂直方向に作用する力、ならびに内部バス245を介して得ることができる前後方向の加速度および車輪速を用いる。これらの情報は、2つの車軸上の個々の摩擦係数を得るために処理される。摩擦係数は、車輪におけるトルク平衡およびCoGでの力の平衡によって求められる。結果として得られる一組の一次方程式を解くことにより、摩擦係数が提供される(詳細な導出については非特許文献3を参照)。一対の導線351を介して供給される摩擦係数は、車両及び運転者の挙動解析システム10の基本変数を示すが、直接測定することはできない。摩擦係数を知ることは、車両及び運転者の挙動解析システム10について定義される必要性の1つである、車両の周辺環境についての情報を有することを意味する。
【0049】
車体横すべり角(VBSSA)を求めるために、状態空間モデル365は、非線形状態空間モデルサブシステム360において設定される。この車体横すべり角は、事故再現理論における基本変数でもある。これを標準的な車両で測定するためには、多大な努力が必要である。それゆえ、車体横すべり角を、導線371へ出力するとともに、車輪モデルサブシステム340へ出力する観測器370によって求めなければならない。この非線形状態モデルサブシステム360は、状態空間モデル365と非線形観測器370を含み、図4に詳細に記載される。
【0050】
図3に示す車両モデル/推定装置30は、測定不可能な変数を得るための1つの例である。基本的に、種々のサブシステムの他の実施態様も可能である。サブシステムの境界条件に応じて多様な可能性が存在する。モデルは、車体横すべり角(導線371)、タイヤ横すべり角(導線343)、摩擦係数(導線351)、重心速度(導線312)、および車両の重量(導線346)を求めることが可能である。さらに、横方向および垂直方向に作用する力ならびに駆動力が計算される。本モデルベースシステムが既存のイベントデータレコーダ(EDR:event data recorder)よりも多くの情報を提供可能であり、それによって事故をより正確に再現できるのは、このためである。
【0051】
図4は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の非線形状態空間モデルの構成を示す図である。
【0052】
非線形状態空間モデルサブシステム360の状態空間モデル365のプロセスは、次の式(1)及び(2)のように記述される。
【0053】
【数1】
Figure 0003763476
【0054】
状態ベクトル(−)xは、測定不可能な変数、例えば車体横すべり角(VBSSA)を含む。なお、(−)xは、式(1)に示すように、xにアンダーラインがあることを表し、以下同様である。
【0055】
初めに、関数(−)f((−)x,(−)u)は、重心における力およびトルクの平衡により設定される。それゆえ、測定可能な変数は測定不可能な変数に変換され、たとえばこの特定の場合では、前後方向および横方向の加速度が車体横すべり角(VBSSA)との関係に置換される。最後に、式(1)および(2)は、測定不可能な状態(−)x、入力ベクトル(−)u、および測定ベクトル(−)yの関係を含む。入力ベクトルは、たとえば前後方向の車輪の力および車輪操舵角からなる。
【0056】
非線形状態空間モデルサブシステム360の非線形観測器370は、次の式(3)及び(4)のように記述されるプロセスモデルによって構成される。
【0057】
【数2】
Figure 0003763476
【0058】
プロセスの開始状態(−)xとモデルの開始状態(−)(∧)xは異なるため、測定されたベクトル(−)yと推定測定ベクトル(−)(∧)yには差があるであろう。なお、(−)(∧)xは、式(3)に示すように、xにアンダーラインがあり、かつxの上にハットがあることを表し、以下同様である。
【0059】
差(−)y−(−)(∧)yは、行列(−)Lに加えられ、モデルにフィードバックされ、結果として式(3)および(4)が得られる。測定不可能な状態空間ベクトル(−)xを正確に推定するために、推定状態空間ベクトル(−)(∧)xは(−)xに等しくならねばならない。換言すると、(−)xと(−)(∧)xとの差は0にならねばならない。これは、推定誤差が長いタイムディスタンスで0になるように、行列(−)Lの要素を選択することによって達成され得る。(−)Lの要素を求める方法、制御技術分野の専門家には既知である。
【0060】
状態空間モデル365は、線形または非線形に設定され得る。前者の場合、単純なルーエンバーガー(LUENBERGER)観測器が、(−)Lの要素を提供する。後者の場合、たとえばZEITZによる非線形観測器を実施する必要がある。後者の方法は、非特許文献3に詳細に記載され、非線形4輪状態空間モデルが導入されている。測定不可能な状態ベクトル(−)xの推定ベクトル(−)(∧)xは、非線形状態空間モデルサブシステム360の出力である。非線形観測器370により、事故再現理論において重要な変数として車体横すべり角(VBSSA)を求めることができる。
【0061】
図5は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の車速推定器の構成を示す図である。
【0062】
車速推定器310は、マイクロプロセッサ上で実行されるアルゴリズムで動作する。したがって、サブシステム3100、3101、3104、3108、3110の出力は、1つまたは複数のマイクロプロセッサ内のレジスタ、および/または1つまたは複数のマイクロプロセッサにアドレス指定された外部メモリに書き込まれる。この車速推定器310に関するより良い構造および概観を提供するために、機能ブロック図により、システム内の信号の流れについて説明する。実際は、レジスタバンクおよび/または外部メモリの内容を意味する矢印は中間結果であり、最終の結果もまた次の処理のために書込まれる。
【0063】
前後方向の加速度および車輪速が内部バス245から得られる。車速推定器310は、より確実な車両の重心速度を、実際の運転状況によって決定するために、加速度および車輪速の冗長情報を用いる。ハードブレーキングや加速状況においては、スリップのため車輪速情報は確実ではない。信号不良が生じる場合もある。加速度信号は、通常、雑音を含み、オフセットおよび/またはオフセットドリフトを受けやすい。運転状況に従って、各センサは多少の重み付けがされる。信号をファジィシステムに加える前に、それら信号を前処理して、系統誤差を排除しなければならない。
【0064】
車輪速変換サブシステム3100は、個々の車輪位置の速度を車両の重心位置の速度に変換し、それゆえ外輪が内輪よりも速くなるカーブでの運転操作による誤差を排除する。加速度信号は、加速度修正サブシステム3101において修正される。この場合、道路勾配による系統誤差を除去しなければならない。一組の導線3102上に出力される変換済みの車輪速はサブシステム3104に入力され、ここでは、変換済みの車輪速間の最大偏差が計算され、同時に最新のサンプリングステップにおける車輪速と推定CoG速度との間の速度差ΔvRijが計算される。
【0065】
この結果、最大偏差Δmaxは、ある特定の時点における4つの車輪速値のうちの最高値と最低値との差となる。したがって、最大偏差Δmaxは概して車輪速の信頼度の測度であり、一方、速度差ΔvRijは、個々の車輪についてのスリップの測定量となる。速度差は、一組の導線3106を介して伝達され、最大偏差は導線3107を介してファジィシステム3108に転送される。
【0066】
修正された加速度は、ファジィシステム3108への別の入力であり、導線3103を介してファジィシステム3108に供給される。特定の運転状況についての発見的知識はルールベースで実施され、ファジィシステム3108内のIF...THENブロックとして示される。ルールベースは、ある特定の信号がその現在の信頼度に従ってどのように重み付けされなければならないかに関する知識を含む。たとえば、車輪速は、大きなブレーキスリップによるABSブレーキング状況では低く重み付けされるであろう。メンバーシップ関数の選択、ルールベース、および非ファジィ化法を伴う完全なファジィシステム3108は、非特許文献3において見ることができる。
【0067】
一組の導線3109を介するファジィシステム3108の出力は、4つの車輪速と加速度(k)に対する重み係数k、k、k、kを含む。重み平均サブシステム3110において、推定重心速度が求められる。
【0068】
【数3】
Figure 0003763476
【0069】
CoG速度は導線312を介して出力されるが、次のサンプリングステップについての速度差ΔvRijを求めるために、導線3105を介してサブシステム3104に戻される。
【0070】
車速推定器310で説明したファジィ速度推定方法は、CoG速度推定の一可能性を示すだけである。他のアプローチは、カルマンフィルタを用いるか、または全地球測位システム(GPS:Global Positioning System)のような衛星測位システムからの評価情報を用いる。非特許文献2を参照されたい。
【0071】
車輪の速度で再現を行うのは信頼度が低すぎる。そのため、CoGの速度が、事故再現における基本変数の1つでもある。
【0072】
図6は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の車輪モデルの構成を示す図である。
【0073】
車輪モデルサブシステム340は、マイクロプロセッサ上のアルゴリズムにおいて動作する。したがって、実際には、中間結果および最終結果が、マイクロプロセッサのレジスタおよび/またはマイクロプロセッサにアドレス指定された外部メモリに書き込まれる。車輪モデルサブシステム340における信号の流れについてよりわかりやすく説明するために、図6に信号の流れを矢印で表すブロック図で示す。
【0074】
車輪モデルサブシステム340は、他のサブモデルにおいて処理される情報、ならびに内部バス245からの情報を用いて、二組の導線341および342を介して前後方向および横方向に車輪に作用する力を、かつ一組の導線343を介してタイヤの横すべり角(TSSA)を出力する。
【0075】
TSSA計算サブシステム3401において、内部バス245を介して得ることができるヨーレート信号、導線312を介して車両の重心速度、導線326を介して車輪操舵角、導線371を介して車体横すべり角(VBSSA)が、タイヤ横すべり角(TSSA)を計算するために必要とされる。車輪の速度とシャーシの速度との関係の設定は、タイヤ横すべり角計算の方程式を導く。このタイヤ横すべり角計算の一組の導線343を介してスリップ計算サブシステム3404に入力されるが、車輪モデルサブシステム340の出力の1つでもある。
【0076】
車輪接地点の速度(VWGCP:velocity of wheel ground contact point)は、VWGCP計算サブシステム3402において求められる。この車輪接地点の速度(VWGCP)は、車輪における実速度ベクトルを示す。車輪接地点の速度は、タイヤ横すべり角(TSSA)により、車輪面(wheel plain)からそれる。
【0077】
このVWGCP計算サブシステム3402は、TSSA計算サブシステム3401と同じ入力を処理し、一組の導線3403を介して各車輪における車輪接地点の速度(VWGCP)を出力する。このVWGCPは、各車輪における異なるカーブ半径を考慮すること、およびそのヨー軸を中心とした車両の動きに対する関係を築くことによって求められる。
【0078】
導線343を介してタイヤ横すべり角(TSSA)と、導線3403における車輪接地点の速度(VWGCP)は、スリップ計算サブシステム3404に入力される。このスリップ計算サブシステム3404で、非特許文献1の著者であるBurckhardtの定義に基づいてスリップが計算され、導線3405を介して摩擦係数確定サブシステム3406に伝達される。
【0079】
この摩擦係数確定サブシステム3406において、車輪のスリップと摩擦係数との特性関係が、後者(すなわち摩擦係数)を得るために用いられる。この関係は経験的実験において見出され、車両工学に精通する者にとって既知の近似方程式で実施された。一組の導線3407における摩擦係数と、一組の導線316を介する車輪に垂直に作用する力は、車輪の力サブシステム3408に入力される。
【0080】
この車輪の力サブシステム3408は、一組の導線341を介して車輪に前後方向に作用する力を、かつ一組の導線342を介して車輪に横方向に作用する力を出力する。車輪の力サブシステム3408は、車輪に垂直に作用する力と摩擦係数の関係を表す方程式を含む。さらに、それらの力は車輪座標系からシャーシ座標系へと変換される。
【0081】
図7は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の駆動モデルの構成を示す図である。
【0082】
駆動モデルサブシステム320において、内部バス245を介して得ることができる前処理済エンジン回転数とスロットル開度は、エンジントルクマップサブシステム3201に入力される。エンジントルクマップは、マイクロプロセッサメモリにおけるルックアップテーブルとして実施される。このエンジントルクマップサブシステム3201は、導線3202を介して車両のエンジントルクを出力する。
【0083】
車重推定器345(図7には図示せず)において実行される車輪上のトルク平衡のためには、車輪の駆動トルクが直接必要である。したがって、動力伝達機構伝達関数サブシステム3203において車両の動力伝達機構(powertrain)がモデル化される。この動力伝達機構伝達関数サブシステム3203において実施される動力伝達機の伝達関数は、動力伝達機の剛性および減衰を考慮するために、たとえば二次システムであり得る。伝達関数の利得は、選択されるギアに応じて適合されなければならない。駆動モデルサブシステム320は、導線321を介して車輪駆動トルクを出力する。
【0084】
図8は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムのモデルベース再現装置の構成を示す図である。
【0085】
インタフェース45は、車両モデル/推定装置30において求められた付加情報と、車両で取得された情報2、4、6の全てを含む。モデルベース再現装置50の目的は、事故前後の物理変数を計算することと、当該装置内で処理されたこれら物理変数および他の全情報を、共通の時刻管理情報(ブロック550で示す)に統合することである。
【0086】
さらに、モデルベース再現装置50において、相対座標系で計算される車両の動きは、絶対座標系、いわゆる慣性座標系に、車両静止後の車両の最終位置を起点として固定される(anchored)。基本的に、再現は単独で入手可能なセンサ情報を用いて実行され得る。しかしながら、測定誤差により、再現された軌跡と「実」軌跡との間に大きな偏差が生じる結果となるであろう。
【0087】
本発明によるモデルベース再現装置50は、たとえばパーソナルコンピュータ上のマイクロプロセッサで実行されるアルゴリズムにおいて動作する。より良い概観を提供するために、モデルベース再現装置50内の信号の流れを矢印で表す。実際には、モデルベース再現装置50の入力信号、サブシステム510、520、530、540における計算の中間結果と最終結果、ならびに出力信号は、レジスタおよび/または1つまたは複数のマイクロプロセッサによりアドレス指定および管理される外部メモリに格納される。
【0088】
車載データネットワーク14(図8には図示せず)を介して送信されるか、またはローカルマイクロプロセッサにおいて生成され、かつ導線501を介して得ることができる時刻管理情報は、全サブシステムに流れる。これは、時刻管理情報がインタフェース45を介して得ることができる複数の情報と比較する基準であり、それゆえ「いつ何が起こったか」について推論することを可能にするからである。
【0089】
一組の導線502は、車両の複数の静的信号を含み、導線501を介する時刻管理情報に加えて、信号状態再現サブシステム510への第2の入力であり、サブシステム510において、スイッチや電源の入り切り信号を経時的な運転者の操作として再現され、導線512を介して入手可能にされる。
【0090】
軌跡再現サブシステム520において、経時的な車両の位置が計算される。したがって、導線501を介して取得された時刻管理情報に加えて、一組の導線503を介して取得されたデータ、複数の導線504を介して車両モデル/推定装置30(図8には図示せず)において得られた測定不可能な変数についての付加情報、および導線18を介して人による情報8が処理される。後者の人による情報8は、たとえば車両の動きの起点を固定するために必要である。この起点は、車両の停止場所の位置を示す。これは、車両及び運転者の挙動解析システム10におけるデータの記録が事故後に停止したときに車両が位置付けられている位置である。
【0091】
例えば、警察によって起点が固定されることは、慣性座標系の固定を意味する。この地点から開始して、車両の軌跡(−)x(n)=[x(n)y(n)]及び車両の向きは、二組の導線501、503、および504を介して得ることができる情報を用いて、たとえば以下の式(6)、(7)に従って再現される。
【0092】
【数4】
Figure 0003763476
【0093】
これは、導線501、503、および504を介する情報が、Δs(n)およびΔΘ(n)を正確に求めるために適切に処理されることを意味する。よってΔs(n)は、2つの連続した位置ベクトル(−)x(n)と(−)x(n+1)の間の距離を表し、一方、ΔΘ(n)は、2つの連続した時点nとn+1の間の車両のCoG速度ベクトルの回転方向の動きを表す。この軌跡再現サブシステム520の出力は、一組の導線522を介して内部バス65に送信される。この出力は、車両の停止位置から事故前の指定された時間における車両の位置までの、慣性座標系における車両の軌跡および向きを含む。
【0094】
車両状態変数再現サブシステム530は、軌跡再現サブシステム520と同じ入力を有し、共通の時刻管理情報に関する事故再現のために重要なパラメータおよび状態変数を一組の導線532を介して出力する。事故再現のために重要な変数は、たとえば、タイヤ横すべり角と、車体横すべり角と、前後方向、横方向、および車輪に垂直方向に作用する力と、対地速度である。
【0095】
軌跡再現サブシステム520及び車両状態変数再現サブシステム530と同じ入力を有する車両周辺環境再現サブシステム540において、経時的な摩擦係数が決定され、道路およびライト状況、またはカーブ半径のような付加情報が一組の導線542を介して出力される。このモデルベース再現装置50において生成される全情報は、内部バス65に入力され、図9に示す査定装置70に送信される。
【0096】
図9は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの査定装置の構成を示す図である。
【0097】
物理変数が計算されるモデルベース再現装置50とは対照的に、査定装置70では、物理変数の判定が実行される。再現および査定の分割が目的ごとに行なわれる。これは、モデルベース再現装置50(図9には図示せず)における物理変数の計算は、いかなる場合にも行なう必要があるが、査定はある特定の国における個別の法律に従って行なわれるからである。
【0098】
車両の軌跡、「複数の静的状態」、目撃者の証言等の事故のシナリオについての情報は、内部バス65を介して査定装置70に供給され、さらに運転者査定サブシステム710、車両周辺環境査定サブシステム720、車両挙動査定サブシステム730についての個々の査定サブシステムに供給される。
【0099】
運転者査定サブシステム710において、たとえば運転者の運転が速すぎたか否かの判定が下される。しかしながら、この判定は、警察による人の情報(たとえば、事故が速度制限区域で発生したか否か等)や、国の法律の特性(たとえば、郊外の道路全般に当てはまる速度制限があるか否か等)に左右される。
【0100】
さらに、判定は、車両挙動査定サブシステム730において査定される車両の挙動によって決まる。たとえば、ブレーキが故障していた場合、これは運転者の過失とはみなされない。したがって、車両と運転者の挙動との相互作用が考慮されなければならず、これは図9において双方向矢印713で示される。
【0101】
双方向矢印712で示される車両周辺環境査定サブシステム720(たとえば、「凍結道路」)と運転者の挙動との間、および車両周辺環境査定サブシステム720と車両の挙動(双方向矢印723)との間の付加的な相互作用も、考慮に入れなければならない。また、査定装置70は、運転者査定サブシステム710による運転者の操作及び挙動と、車両周辺環境査定サブシステム720による車両周辺環境と、車両挙動査定サブシステム730による車両の挙動との関連を評価する関連評価サブシステム740を含む。
【0102】
運転者または車両の誤った挙動を説明するための測度を得るために、車輪速に加えて、たとえば、車両がアンダーステアの状態であったかオーバーステアの状態であったかについて推論するために前輪と後輪の軌道上でのタイヤの横すべり角(TSSA)を得ることができる。運転者によるステアリングの入力についての情報とともに、カーブにおける運転者の挙動を査定することができる。車体横すべり角(VBSSA)および車速の比率、ならびに道路半径および車速の比率が、考慮され、かつ運転者の入力と比較されなければならない。たとえば、上記の比率についてのある特定の閾値を下回るか、または上回る場合、運転者の不適切な挙動が検出され得る。
【0103】
車両の状態を査定するための別の例は、車輪速の解析である。たとえば、車輪速の数値が、短期間での著しい車速の低下を示すものでない場合、ABSシステムの故障という結論を下すことができ、運転者の責任は免除され得る。種々の国の個別の条例により、査定装置70は異なる境界条件に適合されなければならない。査定装置70は、事故状況の説明において優れた役割を果たす。モデルベース再現装置50(図9には図示せず)におけるモデルベース再現と同様に、車両モデル/推定装置30(図9には図示せず)は、危機的な運転状況における運転者および車両の挙動を査定装置70において正確に査定するために必要な部分である。
【0104】
図10は、この発明の実施の形態1に係る車両及び運転者の挙動解析システムの3つの異なる構成(システムアーキテクチャ)を示す図である。
【0105】
図10(a)は、いわゆる「オンラインオンボードシステム(Online-Onboard-System)」のブロック図を示す。入手可能な車両データ2、4、および6は、複数の導線12、16と車載データネットワーク14を介して、車両及び運転者の挙動解析システム10のデータ前処理装置20に入力される。
【0106】
前処理済データは、一組の導線と内部バス25を介して車両モデル/推定装置30に伝達される。車両モデル計算および推定アルゴリズムは、車両モデル/推定装置30においてリアルタイムに(「オンライン」で)実行される。双方向内部バス35は、中間結果を保存およびロードするためにメモリ40とデータを交換するためのものである。
【0107】
事故の発生後、推定および計算されたデータ、ならびに前処理済測定信号が、インタフェース45を介してモデルベース再現装置50にアップロードされる。車両モデル/推定装置30内で取得された全情報は、モデルベース再現装置50において、共通の時刻管理情報に統合され、事故のシナリオが再現される。
【0108】
再現された事故のシナリオは、モデルベース再現装置50から、内部バス65を介して査定装置70に入力される。この査定装置70において、車両および運転者の挙動ならびに環境条件が査定される。
【0109】
図10(a)に示されるシステムアーキテクチャは、車両において実行されている一部分(「オンボード」)101により特徴付けられる。この部分は、サブシステム20におけるデータの前処理と、サブシステム30におけるリアルタイムでのモデル計算および推定を含む。また、サブシステム50において実行されるモデルベース再現と、サブシステム70における事故のシナリオの査定は、オフボードで試験所(ブロック102)において行なわれ、車内では行なわれない。
【0110】
図10(b)は、いわゆる「オフラインオンボードシステム(Offline-Onboard-System)」のブロック図である。車両データ2、4、および6は、複数の導線12、16と、車載データネットワーク14を介して、車両及び運転者の挙動解析システム10のデータ前処理装置20に入力される。
【0111】
データは、内部バス55を介してメモリ40に保存される。事故の発生後、CPU(図10(b)には図示せず)が自身をトリガして車両モデルの計算を開始し、サブシステム30においてモデルベース推定を実行する。「オフラインオンボードシステム」はリアルタイムシステムではない。双方向内部バス35により、中間結果および最終結果をメモリ40に書き込み、必要であればそれらの結果をリロードすることができる。
【0112】
数学的計算の終了後、インタフェース45によってメモリ40からモデルベース再現装置50に車両モデル計算および推定結果がアップロードされる。モデルベース再現装置50において、取得されたデータと、車両モデル/推定装置30において判定された付加的なデータが共通の時刻管理情報に統合され、事故のシナリオが再現される。
【0113】
再現された事故のシナリオは、内部バス65を介して査定装置70に入力される。この査定装置70において、環境条件と同じく、車両と運転者の挙動が査定される。
【0114】
「オフラインオンボードシステム」は、サブシステム20におけるデータの前処理と、サブシステム30における車輪モデルの計算および推定をオンボードで、すなわち車内で実行する(ブロック101)。図10(a)で説明した「オンラインオンボードシステム」とは対照的に、サブシステム30における計算はリアルタイムで計算されない(「オフライン」)。また、サブシステム50におけるモデルベース再現と、サブシステム70における事故のシナリオの査定は、試験所(ブロック102で示す「オフボード」)で行なわれ、車内では行なわれない。
【0115】
図10(c)のブロック図に示す「オフラインオフボードシステム(Offline-Offboard-System)」は、以下のように機能する。車両データ2、4、および6が複数の導線12、16と、車載データネットワーク14を介して車両及び運転者の挙動解析システム10のデータ前処理装置20に供給される。
【0116】
データは、内部バス55を介してメモリ40に保存される。事故の発生後、データはメモリ40からインタフェース35を介して試験所内のパーソナルコンピュータ(PC)にアップロードされる。このパーソナルコンピュータにおいて、車両モデル/推定装置30で各計算が実行される(「オフライン」)。
【0117】
モデル計算および推定が実行された後、その結果はインタフェース45を介してモデルベース再現装置50に入力される。このモデルベース再現装置50において、取得されたデータおよび車両モデル/推定装置30において求められた情報が時刻管理情報に統合され、事故のシナリオが再現される。
【0118】
再現された事故のシナリオは、内部バス65を介して査定装置70に入力される。この査定装置70において、車両と運転者の挙動、ならびに環境条件が査定される。
【0119】
「オフラインオフボードシステム」もまたオンボードブロック101に分割され、これはシステムの車両にインストールされた部分を示す。この「オフラインオフボードシステム」では、データ前処理装置20と、メモリ40への測定および前処理済データの保存行為のみがブロック101においてオンボードで実行される。また、車両モデルおよび推定計算と、モデルベース再現と、事故のシナリオの査定は試験所で行なわれ、これはブロック102における「オフボード」を意味するものであり、リアルタイムではない(「オフライン」)。
【0120】
これにより、図10(c)の「オフラインオフボードシステム」は、図10(a)に示す「オンラインオンボードシステム」および図10(b)に示す「オフラインオンボードシステム」と区別され、後者ではモデルベース再現および査定のみが車外(オフボード)で行なわれる。
【0121】
【発明の効果】
この発明に係る車両及び運転者の挙動解析システムは、以上説明したとおり、車両モデル/推定装置によって、車両の挙動を表す微分方程式の解が得られるだけでなく、推定手法を用いて車両の挙動状態、および車重など車両に搭載されるセンサでは測定不可能な変数を求めることができ、事故等の危険な運転状況を定量的に解析するために必要な情報を記録することができるという効果と、さらに、査定装置によって運転者の挙動、車両の挙動、環境条件の相互関係を明らかにでき、交通規則違反の証明、または必ずしも事故に至らなかった不安定な運転状況の解析の手段としても使用することができるという効果を奏する。
【図面の簡単な説明】
【図1】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの構成を示す図である。
【図2】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムのデータ前処理装置の構成を示す図である。
【図3】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の構成を示す図である。
【図4】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の非線形状態空間モデルの構成を示す図である。
【図5】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の車速推定器の構成を示す図である。
【図6】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の車輪モデルの構成を示す図である。
【図7】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの車両モデル/推定装置の駆動モデルの構成を示す図である。
【図8】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムのモデルベース再現装置の構成を示す図である。
【図9】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの査定装置の構成を示す図である。
【図10】 この発明の実施の形態1に係る車両及び運転者の挙動解析システムの3つの異なる構成を示す図である。
【符号の説明】
2 電気機器、4 車載データネットワーク、6 付加情報、8 人による情報、10 システム、20 データ前処理装置、30 車両モデル/推定装置、40 メモリ、50 モデルベース再現装置、70 査定装置、210 第1の信号変換部、220 信号復号器、230 第2の信号変換部、240 デジタルフィルタ、250 フィルタ、260 アナログ/デジタル変換器、310 車速推定器、315 垂直方向の力サブシステム、320 駆動モデルサブシステム、325 ステアリングモデルサブシステム、330 摩擦サブシステム、335 道路勾配および風力サブシステム、340 車輪モデルサブシステム、345 車重推定器、350 摩擦推定器、360 非線形状態空間モデルサブシステム。

Claims (7)

  1. 車載電気機器の動作情報、車載データネットワークを介して得られるデータ、及び車両を構成する電気装置から得られる付加情報を内部で演算可能な形式に変換を行うデータ前処理手段と、
    前記データ前処理手段からの前処理済データに基づいて車両搭載センサでは測定できない変数を推定する車両モデル/推定手段と、
    前記データ前処理手段からの前処理済データ、及び前記車両モデル/推定手段からの推定された変数を共通の時刻管理情報に統合して事故のシナリオを再現するモデルベース再現手段と、
    前記モデルベース再現手段により再現された事故のシナリオに基づいて車両、及び運転者の挙動を査定する査定手段とを備え
    前記車両モデル/推定手段は、
    前記データ前処理手段からの前処理済データに基づいて、車両に搭載されるセンサでは測定できない変数として、車両の重心速度を求める車速推定部と、
    前記データ前処理手段からの前処理済データに基づいて、車両の動力伝達機構を数式として構築した駆動モデルを用いて駆動トルク求める駆動モデル部と、
    前記前処理済データ、前記車両の重心速度、及び前記駆動トルクに基づいて、車両に搭載されるセンサでは測定できない変数として、車両の重量を推定する車重推定部と、
    前記前処理済データ、及び前記車両の重心速度に基づいて、車輪を数式として構築した車輪モデルを用い、車両に搭載されるセンサでは測定できない変数として、タイヤの横すべり角と車輪に作用する力を計算する車輪モデル部と、
    前記前処理済データ、前記車両の重心速度、前記車両の重量、及び前記車輪に作用する力に基づいて、走行時における車両の状態を数式として構築した状態空間モデルを用い、車両に搭載されるセンサでは測定できない変数として、車体の横すべり角を導出する状態空間モデル部とを含む
    ことを特徴とする車両及び運転者の挙動解析システム。
  2. 車載電気機器の動作情報、車載データネットワークを介して得られるデータ、及び車両を構成する電気装置から得られる付加情報を内部で演算可能な形式に変換を行うデータ前処理手段と、
    前記データ前処理手段からの前処理済データに基づいて車両搭載センサでは測定できない変数を推定する車両モデル/推定手段と、
    前記データ前処理手段からの前処理済データ、及び前記車両モデル/推定手段からの推定された変数を共通の時刻管理情報に統合して事故のシナリオを再現するモデルベース再現手段と、
    前記モデルベース再現手段により再現された事故のシナリオに基づいて車両、及び運転者の挙動を査定する査定手段とを備え、
    前記モデルベース再現手段は、
    前記データ前処理手段からの前処理済データのうち、経時的な運転者の操作による信号を再現する信号状態再現部と、
    前記データ前処理手段からの前処理済データ、前記車両モデル/推定手段からの推定された変数、及び外部から入力されたデータに基づいて経時的な車両の軌跡及び向きを再現する軌跡再現部と、
    前記データ前処理手段からの前処理済データ、前記車両モデル/推定手段からの推定された変数、及び外部から入力されたデータに基づいて車両状態を表す経時的な変数を求める車両状態変数再現部と、
    前記データ前処理手段からの前処理済データ、前記車両モデル/推定手段からの推定された変数、及び外部から入力されたデータに基づいて経時的な車両の周辺環境を再現する車両周辺環境再現部とを含む
    ことを特徴とする車両及び運転者の挙動解析システム。
  3. 前記データ前処理手段は、
    前記車載電気機器の動作情報、及び前記車両を構成する電気装置から得られる付加情報の信号レベルをマイクロプロセッサで使用する信号レベルに変換する信号変換部と、
    前記車載データネットワークから得られるデータを復号する信号復号部と、
    前記信号変換部及び前記信号復号部の少なくともいずれかからの信号のフィルタ処理を行うフィルタ処理部と、
    前記信号変換部からのアナログ信号をデジタル信号に変換するアナログ/デジタル変換部
    を含むことを特徴とする請求項1又は2記載の車両及び運転者の挙動解析システム。
  4. 前記査定手段は、
    前記モデルベース再現手段により再現された事故のシナリオに基づいて運転者の操作及び挙動を査定する運転者査定手段と、
    前記モデルベース再現手段により再現された事故のシナリオに基づいて車両の挙動を査定する車両挙動査定手段と、
    前記モデルベース再現手段により再現された事故のシナリオに基づいて車両周辺環境を査定する車両周辺環境査定手段と、
    前記運転者査定手段による運転者の操作及び挙動と、前記車両挙動査定手段による車両の挙動と、前記車両周辺環境査定手段による車両周辺環境との関連を評価する関連評価手段
    を含むことを特徴とする請求項1又は2記載の車両及び運転者の挙動解析システム。
  5. 前記データ前処理手段、及び前記車両モデル/推定手段は、車内において実行され、かつ
    前記車両モデル/推定手段は、リアルタイムで実行され、
    前記モデルベース再現手段、及び前記査定手段は、車外において実行される
    ことを特徴とする請求項1又は2記載の車両及び運転者の挙動解析システム。
  6. 前記データ前処理手段、及び前記車両モデル/推定手段は、車内において実行され
    記モデルベース再現手段、及び前記査定手段は、車外において実行される
    ことを特徴とする請求項1又は2記載の車両及び運転者の挙動解析システム。
  7. 前記データ前処理手段、車内において実行され、
    前記車両モデル/推定手段、前記モデルベース再現手段、及び前記査定手段は、車外において実行される
    ことを特徴とする請求項1又は2記載の車両及び運転者の挙動解析システム。
JP2003153184A 2003-05-29 2003-05-29 車両及び運転者の挙動解析システム Expired - Fee Related JP3763476B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2003153184A JP3763476B2 (ja) 2003-05-29 2003-05-29 車両及び運転者の挙動解析システム
US10/680,219 US7239945B2 (en) 2003-05-29 2003-10-08 System for analyzing vehicle and driver behavior
DE10353650A DE10353650A1 (de) 2003-05-29 2003-11-17 System zum Analysieren eines Fahrzeug- und Fahrerverhaltens

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003153184A JP3763476B2 (ja) 2003-05-29 2003-05-29 車両及び運転者の挙動解析システム

Publications (2)

Publication Number Publication Date
JP2004352116A JP2004352116A (ja) 2004-12-16
JP3763476B2 true JP3763476B2 (ja) 2006-04-05

Family

ID=33447824

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003153184A Expired - Fee Related JP3763476B2 (ja) 2003-05-29 2003-05-29 車両及び運転者の挙動解析システム

Country Status (3)

Country Link
US (1) US7239945B2 (ja)
JP (1) JP3763476B2 (ja)
DE (1) DE10353650A1 (ja)

Families Citing this family (57)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9563869B2 (en) 2010-09-14 2017-02-07 Zonar Systems, Inc. Automatic incorporation of vehicle data into documents captured at a vehicle using a mobile computing device
US10185455B2 (en) 2012-10-04 2019-01-22 Zonar Systems, Inc. Mobile computing device for fleet telematics
US7359821B1 (en) 2002-06-11 2008-04-15 Injury Sciences Llc Methods and apparatus for using black box data to analyze vehicular accidents
US7805228B2 (en) 2004-08-19 2010-09-28 Spx Corporation Vehicle diagnostic device
US7885739B2 (en) * 2004-08-19 2011-02-08 Spx Corporation Open-ended vehicle diagnostic device interface
US7826948B2 (en) * 2004-10-15 2010-11-02 Ford Global Technologies Vehicle loading based vehicle dynamic and safety related characteristic adjusting system
US7430465B2 (en) * 2004-11-17 2008-09-30 Spx Corporation Open-ended PC host interface for vehicle data recorder
EP1672389B1 (en) * 2004-12-20 2013-08-07 Volvo Car Corporation Method for estimating a traffic situation
DE102005004742A1 (de) * 2005-02-02 2006-08-10 Robert Bosch Gmbh Sicherheitssystem
US20070005404A1 (en) * 2005-06-09 2007-01-04 Drive Diagnostics Ltd. System and method for providing driving insurance
US7953536B2 (en) * 2005-07-29 2011-05-31 GM Global Technology Operations LLC Inertial sensor software architecture security method
WO2007094191A1 (ja) * 2006-02-16 2007-08-23 Pioneer Corporation 移動体保険サービス支援方法、移動体保険サービス支援プログラム、記録媒体、および移動体保険サービス支援装置
US7502772B2 (en) * 2006-02-27 2009-03-10 Injury Sciences Llc Method and apparatus for obtaining and using impact severity triage data
US7266468B1 (en) * 2006-03-03 2007-09-04 Perceptron, Inc. Structural data analysis system
US20070271079A1 (en) * 2006-05-17 2007-11-22 Kentaro Oguchi Simulator for Vehicle Radio Propagation Including Shadowing Effects
US8239220B2 (en) * 2006-06-08 2012-08-07 Injury Sciences Llc Method and apparatus for obtaining photogrammetric data to estimate impact severity
US7698086B2 (en) * 2006-06-08 2010-04-13 Injury Sciences Llc Method and apparatus for obtaining and using event data recorder triage data
US10056008B1 (en) 2006-06-20 2018-08-21 Zonar Systems, Inc. Using telematics data including position data and vehicle analytics to train drivers to improve efficiency of vehicle use
US20130164715A1 (en) 2011-12-24 2013-06-27 Zonar Systems, Inc. Using social networking to improve driver performance based on industry sharing of driver performance data
US20130164713A1 (en) 2011-12-23 2013-06-27 Zonar Systems, Inc. Method and apparatus for gps based slope determination, real-time vehicle mass determination, and vehicle efficiency analysis
DE102007002791A1 (de) * 2007-01-18 2008-07-24 Wabco Gmbh Verfahren und Vorrichtung zur Bestimmung der Geschwindigkeit eines Fahrzeugs
US8195364B2 (en) * 2007-02-12 2012-06-05 Deere & Company Perception model for trajectory following autonomous and human augmented steering control
US8498796B2 (en) * 2007-02-12 2013-07-30 Deere & Company Perception model for trajectory following autonomous and human augmented speed control
US8050856B2 (en) * 2007-04-18 2011-11-01 Chrysler Group Llc Methods and systems for powertrain optimization and improved fuel economy
DE102007034142B4 (de) * 2007-07-21 2016-03-10 Bayerische Motoren Werke Aktiengesellschaft Verfahren zur Plausibilisierung einer die Gierrate oder die Querbeschleunigung oder den Lenkwinkel eines Fahrzeugs repräsentierenden Zustandsgröße
DE102007034143B4 (de) * 2007-07-21 2016-03-10 Bayerische Motoren Werke Aktiengesellschaft Verfahren zum Bestimmen einer Modell-Gierrate und/oder einer Modell-Querbeschleunigung sowie eines Modell-Lenkwinkels
ES2341633B1 (es) * 2008-01-25 2011-04-08 Fundacion Para Progreso Soft Computing Detector personalizado de conduccion anomala de un automovil.
US8340855B2 (en) 2008-04-22 2012-12-25 Spx Corporation USB isolation for vehicle communication interface
JP5585229B2 (ja) * 2010-06-16 2014-09-10 富士通株式会社 車両の挙動推定装置、プログラム、及び、方法
WO2012018733A2 (en) 2010-08-03 2012-02-09 Spx Corporation Vehicle diagnostic, communication and signal delivery system
JP5146523B2 (ja) 2010-11-26 2013-02-20 株式会社デンソー 車両用挙動データ記録装置
US10431020B2 (en) 2010-12-02 2019-10-01 Zonar Systems, Inc. Method and apparatus for implementing a vehicle inspection waiver program
US9527515B2 (en) 2011-12-23 2016-12-27 Zonar Systems, Inc. Vehicle performance based on analysis of drive data
US8764126B2 (en) * 2011-05-03 2014-07-01 Robert Bosch Gmbh Fuzzy logic based brake control
DE102011084204A1 (de) * 2011-10-10 2013-04-11 Robert Bosch Gmbh Verfahren zum Ansteuern von Sicherheitsaktuatorik eines Kraftfahrzeugs
CN102592330B (zh) * 2012-01-19 2014-08-06 上海交通大学 基于电动车电池数据提取道路车辆行驶状况的方法
US10061745B2 (en) 2012-04-01 2018-08-28 Zonar Sytems, Inc. Method and apparatus for matching vehicle ECU programming to current vehicle operating conditions
US9424696B2 (en) 2012-10-04 2016-08-23 Zonar Systems, Inc. Virtual trainer for in vehicle driver coaching and to collect metrics to improve driver performance
CN104837403A (zh) * 2012-11-27 2015-08-12 佛吉亚汽车座椅有限责任公司 具有集成传感器的车辆座椅
KR101438606B1 (ko) 2012-12-20 2014-09-15 현대자동차 주식회사 이종 센서의 융복합을 위한 가상 센서 네트워크 시스템
WO2014149042A1 (en) * 2013-03-20 2014-09-25 International Truck Intellectual Property Company, Llc Smart cruise control system
CN106068097B (zh) 2014-02-20 2020-09-29 佛吉亚汽车座椅有限责任公司 整合有传感器的车辆座椅
US9600541B2 (en) * 2014-05-02 2017-03-21 Kookmin University Industry Academy Corporation Foundation Method of processing and analysing vehicle driving big data and system thereof
US10046671B2 (en) 2015-08-14 2018-08-14 Faurecia Automotive Seating, Llc Occupant-recognition system for vehicle seat
DE102015218810A1 (de) 2015-09-29 2017-03-30 Continental Teves Ag & Co. Ohg Verfahren zum Auswählen von Lokalisierungsalgorithmen in einem Fahrzeug
CN106323269B (zh) * 2015-12-10 2019-06-07 上海思岚科技有限公司 自主定位导航设备、定位导航方法及自主定位导航***
US10918337B2 (en) 2016-06-03 2021-02-16 Faurecia Automotive Seating, Llc Vehicle seat with integrated sensors
US11197637B2 (en) 2016-06-20 2021-12-14 Faurecia Automotive Seating, Llc Control system for a vehicle seat
US10730524B2 (en) 2016-09-07 2020-08-04 Faurecia Automotive Seating, Llc Vehicle seat
WO2018150888A1 (ja) * 2017-02-14 2018-08-23 パイオニア株式会社 情報記録装置、情報記録方法及び情報記録用プログラム
CN106953856B (zh) * 2017-03-16 2020-02-21 大连理工大学 一种结合车辆行驶状态分析车载网can总线协议中数据有效性的方法
US11083379B2 (en) 2017-08-02 2021-08-10 Faurecia Automotive Seating, Llc Health-monitoring seat cover
US10343695B2 (en) * 2017-11-14 2019-07-09 General Motors Llc Controlling vehicle functionality based on geopolitical region
JP7001009B2 (ja) * 2018-07-11 2022-01-19 トヨタ自動車株式会社 風データ推定装置
CN110070642A (zh) * 2019-03-22 2019-07-30 天津大学 一种基于深度学习的交通事故责任评估方法及装置
US20210342500A1 (en) * 2020-05-01 2021-11-04 Steering Solutions Ip Holding Corporation Systems and methods for vehicle modeling
CN112406767B (zh) * 2020-11-11 2022-08-16 广州小鹏汽车科技有限公司 车辆陷坑事故的监控方法、监控装置、服务器和存储介质

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4132981A1 (de) 1991-01-23 1992-08-06 Schimmelpfennig Karl Heinz Verfahren zur rekonstruktion der bewegungstrajektorie eines strassenfahrzeugs
DE4216301C2 (de) 1992-05-16 1997-05-22 Daimler Benz Ag Verfahren zur Bestimmung das Fahrverhalten charakterisierender Größen
DE19509711C2 (de) 1994-03-09 1997-10-16 Mannesmann Ag Unfalldatenschreiber und Verfahren zur Analyse eines Unfallgeschehens
US6223125B1 (en) * 1999-02-05 2001-04-24 Brett O. Hall Collision avoidance system
US6847892B2 (en) * 2001-10-29 2005-01-25 Digital Angel Corporation System for localizing and sensing objects and providing alerts
JP4569798B2 (ja) 2001-07-30 2010-10-27 マツダ株式会社 車両事故解析装置

Also Published As

Publication number Publication date
US20040243368A1 (en) 2004-12-02
JP2004352116A (ja) 2004-12-16
US7239945B2 (en) 2007-07-03
DE10353650A1 (de) 2004-12-23

Similar Documents

Publication Publication Date Title
JP3763476B2 (ja) 車両及び運転者の挙動解析システム
CN113748030B (zh) 用于车辆轮胎性能建模和反馈的***和方法
AU762134B2 (en) Monitoring system for determining and communicating a cost of insurance
Stephant et al. Virtual sensor: Application to vehicle sideslip angle and transversal forces
Doumiati et al. A method to estimate the lateral tire force and the sideslip angle of a vehicle: Experimental validation
JP6796527B2 (ja) 車両状態監視装置、車両状態監視システム及び車両状態監視方法
Allotta et al. Evaluation of odometry algorithm performances using a railway vehicle dynamic model
US10940865B2 (en) Method and device for determining a state of a roadway of a vehicle
CN106816020A (zh) 基于数据分析的交通事故信息处理方法
Kidambi et al. Accuracy and robustness of parallel vehicle mass and road grade estimation
Hiemer Model based detection and reconstruction of road traffic accidents
WO2020193860A1 (en) Modelling dynamics of a vehicle
Kissai et al. Importance of vertical dynamics for accurate modelling, friction estimation and vehicle motion control
US10414409B2 (en) Method for determining a dangerous driving indicator of a vehicle
Antonyan et al. The algorithm of diagnosing the development of a skid when driving a two-axle vehicle
Hong et al. An interacting multiple-model-based algorithm for driver behavior characterization using handling risk
Ghandour et al. Risk indicators prediction based on the estimation of tire/road forces and the maximum friction coefficient: Experimental validation
Bastiaan Physical Validation Testing of a Smart Tire Prototype for Estimation of Tire Forces
Albinsson Online and offline identification of Tyre model parameters
US20230391342A1 (en) Road surface damage detection device, road surface damage detection method, and storage device
Ghandour et al. Risk indicators anticipation based on the vehicle dynamics anticipation to avoid accidents
US20240118175A1 (en) System and method for identifying a tire contact length from radial acceleration signals
US20220219680A1 (en) Information processing device, information processing method, and non-transitory computer-readable medium storing information processing program
Ko et al. Robust Mass Estimation for Enhanced Braking Distance Calculation in Electric Vehicles with Autonomous Emergency Braking Systems
Jiang Real-time estimation and diagnosis of vehicle's dynamics states with low-cost sensors in different driving condition

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050905

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050920

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051108

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060110

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060112

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100127

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100127

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110127

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120127

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130127

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees