特開 2022-134169 軌道状態推定方法、軌道状態推定装置及び車両

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

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

(11)【公開番号】P2022134169

(43)【公開日】2022-09-15

(54)【発明の名称】軌道状態推定方法、軌道状態推定装置及び車両

(51)【国際特許分類】

   B61K 9/08 20060101AFI20220908BHJP

【ＦＩ】

B61K9/08

【審査請求】未請求

【請求項の数】6

【出願形態】ＯＬ

(21)【出願番号】P 2021033126

(22)【出願日】2021-03-03

(71)【出願人】

【識別番号】000006655

【氏名又は名称】日本製鉄株式会社

(71)【出願人】

【識別番号】000203977

【氏名又は名称】日鉄テックスエンジ株式会社

(74)【代理人】

【識別番号】110002044

【氏名又は名称】弁理士法人ブライタス

(72)【発明者】

【氏名】加藤 宗

(72)【発明者】

【氏名】久積 和正

(72)【発明者】

【氏名】君塚 清

(72)【発明者】

【氏名】坂田 哲郎

(72)【発明者】

【氏名】寺地 平

(57)【要約】

【課題】稼働車両を利用して簡易に、かつ、精度よく軌道検測することの可能な、軌道状態推定方法を提供する。
【解決手段】車両の軸箱に設置された慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定ステップと、輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出ステップと、を含む、軌道状態推定方法が提供される。
【選択図】図８

【特許請求の範囲】

【請求項1】

車両の軸箱に設置された慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定ステップと、
前記輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出ステップと、
を含む、軌道状態推定方法。

【請求項2】

前記傾斜角推定ステップでは、
前記慣性計測装置により計測された３軸加速度から、車両走行時の並進加速度成分をフィルタリング処理により抽出し、
抽出された前記並進加速度成分に基づいて、カルマンフィルタの観測モデルの誤差共分散行列を設定し、
前記誤差共分散行列が設定された観測モデルにより構成される状態空間モデルから定式化されたカルマンフィルタを用いて、前記３軸加速度及び前記３軸角速度から前記輪軸の水平傾斜角及び前後傾斜角を推定し、
前記軌道状態算出ステップでは、前記輪軸の水平傾斜角に基づいて、前記軌道状態情報として平面性変位が算出される、請求項１に記載の軌道状態推定方法。

【請求項3】

前記並進加速度成分は、前記３軸加速度のうち左右加速度及び前後加速度に対して、軌道変位の卓越する周波数成分を除去するバンドパスフィルタ処理により抽出される、請求項２に記載の軌道状態推定方法。

【請求項4】

前記誤差共分散行列は、対角成分に、抽出された前記並進加速度成分に対して所定の係数を乗じた値を加算して設定される、請求項２または３に記載の軌道状態推定方法。

【請求項5】

車両の軸箱に設置された慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定部と、
前記輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出部と、
を備える、軌道状態推定装置。

【請求項6】

車両の軸箱に設置された慣性計測装置と、
前記慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定部、及び、
前記輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出部、を有する軌道状態推定装置と、
を備える、車両。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、軌道状態推定方法、軌道状態推定装置及び車両に関する。

【背景技術】

【0002】

レール上を走行する車両の脱線要因の１つとして、元位置に対するレールの不整（以下、「軌道変位」ともいう。）がある。軌道変位は脱線に直接的に影響することから、軌道変位を適切に点検し整備することが、脱線のリスクを低減するために重要である。

【0003】

軌道変位の点検（以下、「軌道検測」ともいう。）は、例えば、手押し式の軌道検測装置を用いて行われている。しかし、軌道検測対象となるレールが長くなると、高い頻度で検測することは難しい。また、手押し式の軌道検測装置による軌道検測は、脱線時の状況とは異なる無負荷時状態での静的検測である。このため、車両走行時の動的検測を高い頻度で実施できることが望まれている。

【0004】

車両走行時の動的検測を高い頻度で実施するための一手法として、稼働車両を利用して軌道検測を行うことが考えられる。車両を利用した軌道検測装置として、例えば、特許文献１には、車両の台車枠に各種検出器から構成される検出器ユニットを取り付け、車両走行時に慣性正矢法に基づき軌道変位５項目を検測する技術が開示されている。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特許第３４１１８６１号公報

【発明の概要】

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

【0006】

しかし、特許文献１に記載の技術では、検出器ユニットにより軌道変位５項目を検測することはできるが、検出器ユニットによる検測機構が複雑であり高価となる。また、大型の検出器ユニットを装架する必要があるため、車両の構造によっては計測位置が車輪から離れる場合があり、車両載荷時のレールの変形を十分にとらえられない可能性がある。

【0007】

そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、稼働車両を利用して簡易に、かつ、精度よく軌道検測することの可能な、軌道状態推定方法、軌道状態推定装置及び車両を提供することにある。

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

【0008】

上記課題を解決するために、本発明のある観点によれば、車両の軸箱に設置された慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定ステップと、輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出ステップと、を含む、軌道状態推定方法が提供される。

【0009】

傾斜角推定ステップでは、慣性計測装置により計測された３軸加速度から、車両走行時の並進加速度成分をフィルタリング処理により抽出し、抽出された並進加速度成分に基づいて、カルマンフィルタの観測モデルの誤差共分散行列を設定し、誤差共分散行列が設定された観測モデルにより構成される状態空間モデルから定式化されたカルマンフィルタを用いて、３軸加速度及び３軸角速度から輪軸の水平傾斜角及び前後傾斜角を推定し、軌道状態算出ステップでは、輪軸の水平傾斜角に基づいて、軌道状態情報として平面性変位が算出されるようにしてもよい。

【0010】

また、並進加速度成分は、３軸加速度のうち左右加速度及び前後加速度に対して、軌道変位の卓越する周波数成分を除去するバンドパスフィルタ処理により抽出してもよい。

【0011】

誤差共分散行列は、対角成分に、抽出された並進加速度成分に対して所定の係数を乗じた値を加算して設定してもよい。

【0012】

また、上記課題を解決するために、本発明の別の観点によれば、車両の軸箱に設置された慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定部と、輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出部と、を備える、軌道状態推定装置が提供される。

【0013】

さらに、上記課題を解決するために、本発明の別の観点によれば、車両の軸箱に設置された慣性計測装置と、慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する傾斜角推定部、及び、輪軸の水平傾斜角及び前後傾斜角に基づいて、軌道状態を表す軌道状態情報を算出する軌道状態算出部、を有する軌道状態推定装置と、を備える、車両が提供される。

【発明の効果】

【0014】

以上説明したように本発明によれば、稼働車両を利用して簡易に、かつ、精度よく軌道検測することができる。

【図面の簡単な説明】

【0015】

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

【図2】平面性変位を説明するための模式図である。

【図3】慣性計測装置による計測値を説明するための説明図である。

【図4】慣性計測装置による計測値を説明するための輪軸断面模式図である。

【図5】並進加速度を表す指標Ｒ_ｘの一算出例を示すグラフである。

【図6】並進加速度を表す指標Ｒ_ｙの一算出例を示すグラフである。

【図7】同実施形態に係る軌道状態推定装置の構成を示す機能ブロック図である。

【図8】同実施形態に係る軌道状態推定方法の一例を示すフローチャートである。

【図9】実施例における平面性変位の推定結果を示すグラフである。

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

【0016】

以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。

【0017】

［１．車両構成］
まず、図１に基づいて、本発明の一実施形態に係る車両１の概略構成について説明する。図１は、本実施形態に係る車両１の概略構成を示す模式図である。

【0018】

本実施形態に係る車両１は、延設された一対のレール５上を走行可能に構成されている。車両１は、図１に示すように、車両の主構造部分である構体からなる車体部１０と、車体部１０を支持し、走行機構を有する台車部２０とにより構成される。なお、以下では、車両進行方向（前後方向）をＸ方向、幅方向（左右方向）をＹ方向、高さ方向（上下方向）をＺ方向とする。

【0019】

図１に示す車両１の台車部２０には、車両進行方向に４つの車輪対が配置されている。各車輪対は、幅方向に対となる２つの車輪２３が輪軸２５によって連結されることにより構成されている。各輪軸２５の両端には、輪軸２５の軸受部分である軸箱２７が設けられている。軸箱２７は、コイルばね２９により台車枠２１に接続されている。なお、図１では、輪軸２５を示すために、１つの車輪２３にのみ軸箱２７及びコイルばね２９を記載しているが、すべての車輪２３が軸箱２７及びコイルばね２９を備えている。

【0020】

本実施形態に係る車両１は、少なくとも１つの車輪２３の軸箱２７に慣性計測装置（ＩＭＵ：Inertial Measurement Unit）３０を備えている。慣性計測装置３０は、３次元の慣性運動（直行３軸方向の並進運動及び回転運動）を検出する装置である。慣性計測装置３０は、加速度センサにより並進運動を検出し、角速度センサにより回転運動を検出する。慣性計測装置３０による計測値は、車両１が走行するレール５の軌道状態を推定する軌道状態推定装置（図３の軌道状態推定装置１００）へ出力される。

【0021】

また、本実施形態に係る車両１は、レール５上の車両１の絶対位置を検出する位置検出装置（図示せず。）を備えている。位置検出装置は、例えばＲＦＩＤ（Radio Frequency Identification）リーダであってもよい。この場合、一対のレール５の幅方向中間位置に予め設定されているＲＦＩＤタグをＲＦＩＤリーダによって読み取ることにより、車両１の絶対位置を検出することができる。なお、車両１の絶対位置を検出する方法はかかる例に限定されず、周知の技術を用いてもよい。例えば、車両１の絶対位置として、人工衛星からの電波を用いて得られる位置情報を利用してもよい。

【0022】

［２．軌道状態の推定］
レールの軌道状態としては、平面性変位、通り変位、高低変位、水準変位、軌間変位がある。以下では、軌道状態として、脱線への寄与の大きい平面性変位を推定する場合について説明する。平面性変位とは、図２に示すように、軌道の平面に対するねじれの状態をいい、車両進行方向に所定の間隔ｄだけ離れた２地点（例えば、点Ｐ１及び点Ｐ３の地点と点Ｐ２と点Ｐ４の地点）における左右のレール５の水準差ＣＬ（ｘ）により表される。具体的には、一般車両においては固定軸距が５ｍ程度であることから、平面性変位Ｔｗ（ｘ）は、下記式（１）により表される。

【0023】

Ｔｗ（ｘ）＝ＣＬ（ｘ）－ＣＬ（ｘ－５） ・・・（１）

【0024】

ここで、車輪２３の浮きが存在しないと仮定したとき、輪軸２５の鉛直変位が軌道の鉛直変位に等しいとみなすことができる。このため、軸箱２７に設置された慣性計測装置３０の計測値に基づき輪軸２５のロール角θｘを推定することで、下記式（２）より水準差ＣＬ（ｘ）を推定することができる。なお、式（２）において、Ｌは基準軌間長を表す。

【0025】

ＣＬ（ｘ）＝Ｌｔａｎθｘ ・・・（２）

【0026】

本実施形態では、軸箱２７に設置された慣性計測装置３０の計測値のみを用いて、車輪直下の軌道変位を検測する。慣性計測装置３０を軸箱２７に設置することで、車輪直下の軌道変位を検測することができるが、車両走行に伴う衝撃的な振動等のノイズも計測されやすくなる。このため、慣性計測装置３０の計測値から推定された軌道変位には、慣性計測装置３０が計測軸方向に並進運動することで計測される、車両走行に伴う衝撃的な振動応答をはじめとした加速度成分（以下、「並進加速度成分」ともいう。）による誤差が発生しやすくなる。そこで、本実施形態では、慣性計測装置３０の計測値に含まれる車両走行時の並進加速度成分による誤差影響を低減して、平面性変位を推定する手法を提示する。

【0027】

以下、軌道状態として平面性変位を推定する軌道状態推定手法の詳細と、軌道状態推定手法に基づく軌道状態推定処理を実行する軌道状態推定装置、及び、軌道状態推定方法について説明する。

【0028】

［２－１．軌道状態推定手法］
まず、図３～図６に基づいて、本実施形態に係る軌道状態を推定する軌道状態推定手法について説明する。図３は、慣性計測装置３０による計測値を説明するための説明図である。図４は、慣性計測装置３０による計測値を説明するための輪軸断面模式図である。図５は、並進加速度を表す指標Ｒ_ｘの一算出例を示すグラフである。図６は、並進加速度を表す指標Ｒ_ｙの一算出例を示すグラフである。

【0029】

以下の説明においては、慣性計測装置３０により計測された時刻ｋにおける３軸加速度の計測値をａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ、３軸角速度の計測値をω_ｘｋ、ω_ｙｋ、ω_ｚｋとする。また、図３に示すように、時刻ｋにおけるロール角をθ_ｋ、ピッチ角をφ_ｋ、ヨー角をψ_ｋとする。

【0030】

本実施形態では、カルマンフィルタ（Kalman Filter）を用いて、慣性計測装置３０の計測値から輪軸２５のロール角θｘを推定する。カルマンフィルタは、データ同化手法の１つであり、ベイズ推定を利用し観測量との誤差を最小化する状態量を逐次的に推定する手法である。３軸加速度の計測値から得られるロール角推定値及びピッチ角推定値を観測量とし、絶対座標におけるロール角速度及びピッチ角速度の積分による時間更新式から、カルマンフィルタを定式化することができる。

【0031】

状態ベクトルｘ_ｋは、下記式（３）で定義される。

【0032】

【数1】

【0033】

また、観測ベクトルｙ_ｋは、図４に示す、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの低周波成分を抽出した重力加速度ｇの３軸方向の射影の関係から、下記式（４）のように表すことができる。観測ベクトルｙ_ｋにおいて、θ_０ｋは時刻ｋにおけるロール角の推定値、φ_０ｋは時刻ｋにおけるピッチ角の推定値である。

【0034】

【数2】

【0035】

３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの重力加速度射影成分は、加速度低周波成分であり、例えば０～１Ｈｚの成分とする。傾斜角推定部１１０は、例えばローパスフィルタにより、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから重力加速度射影成分を抽出し得る。

【0036】

時間更新式は、絶対座標系における角度の積分から導かれ、状態ベクトルについて非線形の関係式となる。そこで、本実施形態では、状態ベクトルｘ_ｋ－１の周辺で線形化を行う拡張カルマンフィルタ（Extended Kalman Filter）を用いる。このとき、システム方程式は、下記式（５）で表される。ｖ_ｋはプロセスノイズである。

【0037】

【数3】

【0038】

また、観測ノイズｗ_ｋ、観測行列Ｈ_ｋをとしたとき、観測方程式は下記式（７）で表される。

【0039】

【数4】

【0040】

上記式（５）～式（７）で表される状態空間モデルに基づくカルマンフィルタにより、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ、及び、３軸角速度の計測値ω_ｘｋ、ω_ｙｋ、ω_ｚｋからロール角θ_ｋ及びピッチ角φ_ｋを推定することができる。

【0041】

ここで、本実施形態では、慣性計測装置３０の計測値に含まれる車両走行時の並進加速度成分による誤差影響を低減するため、傾斜角推定部１１０は、誤差因子となる車両走行時の３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから並進加速度成分をフィルタリング処理により抽出し、抽出された並進加速度成分に基づいてカルマンフィルタの観測誤差共分散行列を設定する。

【0042】

本願発明者は、角度推定値の誤差要因を調べた結果、急曲線部分において車輪フランジがレールに接触する際に発生する衝撃応答が、角度推定値の誤差要因となっていることを知見した。角速度の積分から角度を算出する上記式（５）、及び、加速度低周波成分（すなわち重力加速度の射影成分）から角度を算出する上記式（７）において、プロセスノイズｖ_ｋ及び観測ノイズｗ_ｋは平均値ゼロのガウスノイズと仮定される。このノイズ設定は各関係式の誤差成分を表す。プロセスノイズｖ_ｋは、角速度データの計測誤差、及び、積分誤差を意味する。観測ノイズｗ_ｋは、加速度データの計測誤差、及び、加速度から抽出された重力加速度射影成分の推定誤差を意味する。

【0043】

急曲線部分において、車輪フランジがレールに接触する際には、角速度計測値及び加速度計測値のいずれにも衝撃的な応答が見られる場合があるが、衝撃的な応答は全周波数帯において卓越するという特徴がある。ここで、上記式（５）は衝撃応答の有無にかかわらず成立するが、上記式（７）の観測ベクトルｙ_ｋは、重力加速度射影成分から算出された角度であり、慣性計測装置３０に並進加速度が作用する場合には誤差要因となる。すなわち、このような走行区間においては、観測ノイズｗ_ｋを大きな値に設定することにより、上記式（７）の確からしさを低減させる（換言すると、上記式（５）の確からしさを増大させる）ことで、実応答を適切に表現できると考えられる。

【0044】

そこで、本実施形態では、慣性計測装置３０に作用する並進加速度の大きさに応じて観測ノイズｗ_ｋを設定し、角度推定値の精度向上を図る。具体的には、並進加速度を表す指標Ｒ_ｘ、Ｒ_ｙを導入し、並進加速度成分の大きさに応じた観測ノイズｗ_ｋの共分散行列Ｒを逐次的に設定する。

【0045】

観測ベクトルｙ_ｋのロール角θ_０ｋ及びピッチ角φ_０ｋの算出精度は、それぞれ左右加速度ａ_ｙｋ、前後加速度ａ_ｘｋに支配され、いずれも並進加速度により誤差が生じうる。左右加速度ａ_ｙｋについては急曲線通過時の衝撃応答が並進加速度として現れ、前後加速度ａ_ｘｋについては車両速度を一定に保つための加減速が並進加速度として現れる。そこで、それぞれの並進加速度を表す指標Ｒ_ｘ、Ｒ_ｙを定義する。

【0046】

指標Ｒ_ｘは、前後加速度の測定値ａ_ｘからバンドパスフィルタにより高周波成分のみを抽出し、抽出した値を二乗した後、所定区間（移動平均サンプル数）の移動平均を取った値とする。指標Ｒ_ｙは、左右加速度の測定値ａ_ｙからバンドパスフィルタにより高周波成分のみを抽出し、抽出した値を二乗した後、所定区間（移動平均サンプル数）の移動平均を取った値とする。指標Ｒ_ｘ、Ｒ_ｙは、低周波の応答（すなわち、重力加速度射影成分から推定される傾斜角）を補正する指標であることから、二乗して移動平均をとることにより高周波を除いた値とすることで、推定値を安定化させる。例えば、抽出する高周波成分は１～７Ｈｚの周波数成分とし、移動平均サンプル数は２５０サンプル時間としてもよい。

【0047】

衝撃応答あるいは加減速による応答は、外部から作用する強制的な応答であり、全周波数成分において卓越していると仮定できる。一方、例えば１Ｈｚ未満の低周波数帯では角度変化による成分が卓越しており、例えば７Ｈｚ超の高周波数帯では計測ノイズが卓越していると考えられる。これより、本実施形態では、バンドパスフィルタを用いてこれらの帯域の周波数成分を除去している。

【0048】

詳細には、車両速度がほぼ一定状態であると仮定すると、フィルタリングを施す計測データ中の車両の平均速度をＶ_ａｖｅ［ｍ／ｓ］、軌道変位の最大空間周波数Ｆ_{ｓｐａｃｅ}［ｃｙｃｌｅ／ｍ］としたとき、時刻歴波形のカットオフ周波数Ｆ_ｔｉｍｅ［Ｈｚ］は下記式（８）により表すことができる。

【0049】

Ｆ_ｔｉｍｅ＝Ｆ_{ｓｐａｃｅ}×Ｖ_ａｖｅ ・・・（８）

【0050】

平面性変位の特性は、おおよそ、軌道変位の最大空間周波数Ｆ_{ｓｐａｃｅ}＝１、車両速度Ｖ_ａｖｅ＝１である。これより、平面性変位の推定においては、低周波数帯と高周波数帯との境界とするカットオフ周波数Ｆ_ｔｉｍｅを１Ｈｚとしている。

【0051】

高周波のカットオフ周波数Ｆ_{ｔｉｍｅ＿ｈ}については、経験式として、下記式（９）のような一般性を持った形として表すことができる。

【0052】

Ｆ_{ｔｉｍｅ＿ｈ}＝７×Ｆ_{ｓｐａｃｅ}×Ｖ_ａｖｅ ・・・（９）

【0053】

指標Ｒ_ｘ、Ｒ_ｙの移動平均サンプル数ｓについては、経験式として、計測値のサンプリング周波数ｆｓ［Ｈｚ］から下記式（１０）により表すことができる。

【0054】

ｓ＝ｆｓ／（２×Ｆ_ｔｉｍｅ） ・・・（１０）

【0055】

例えば、サンプリング周波数ｆｓ＝５００Ｈｚ、カットオフ周波数Ｆ_ｔｉｍｅ＝１Ｈｚであるとき、移動平均サンプル数ｓは２５０サンプル時間となる。

【0056】

図５は並進加速度を表す指標Ｒ_ｘの一例であり、図６は並進加速度を表す指標Ｒ_ｙの一例である。図５において、破線枠で示した指標Ｒ_ｘの値が大きい位置では、車両速度を一定に保つための加減速の影響を受けていることを表している。また、図６において、破線枠で示した指標Ｒ_ｙの値が大きい位置では、急曲線通過時の衝撃の影響を受けていることを表している。

【0057】

次に、算出された指標Ｒ_ｘ、Ｒ_ｙを用いて、並進加速度成分の大きさに応じた観測ノイズｗ_ｋを逐次的に設定する。記誤差共分散行列Ｒ（ｋ）は、対角成分に、抽出された並進加速度成分に対して所定の係数を乗じた値を加算して設定してもよい。具体的には、時刻ｋにおける観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）は、並進加速度成分を表す指標Ｒ_ｙ（ｋ）、Ｒ_ｘ（ｋ）、及び、指標Ｒ_ｙ（ｋ）、Ｒ_ｘ（ｋ）に対する係数Ｋ_ｙ、Ｋ_ｘを用いて、下記式（１１）のように設定される。

【0058】

【数5】

【0059】

また、プロセスノイズＱは、時刻によらず一定とし、例えば下記式（１２）のように設定してもよい。

【0060】

【数6】

【0061】

なお、係数Ｋ_ｙ、Ｋ_ｘは、指標Ｒ_ｙ（ｋ）、Ｒ_ｘ（ｋ）から上記式（７）で表される観測方程式のノイズ分散を設定するための値であり、理論的に導出することはできない。このため、係数Ｋ_ｙ、Ｋ_ｘには、平面性変位の推定精度が最も高くなる値が最適値として設定される。プロセスノイズＱも、平面性変位の推定精度が最も高くなる値に基づき、適宜設定される。

【0062】

このように、本実施形態では、時刻ｋにおける観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）を上記式（１１）のように設定し、上記式（７）に示した観測方程式の観測ノイズｗ_ｋを並進加速度成分の大きさに応じて逐次的に設定する。これにより、慣性計測装置３０を軸箱２７に設けることによって計測値に含まれる、車両走行時の並進加速度成分による誤差影響を低減することができる。その結果、時刻ｋにおける観測ベクトルｙ_ｋ、すなわち、ロール角の推定値θ_０ｋ及びピッチ角の推定値φ_０ｋを精度よく求めることができる。式（５）～式（７）で表され、式（１１）により観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）が設定された状態空間モデルに基づくカルマンフィルタにより、ロール角θ_ｋ及びピッチ角φ_ｋが推定されれば、上記式（１）及び式（２）を用いて、平面性変位Ｔｗ（ｘ）を算出することができる。

【0063】

［２－２．軌道状態推定装置］
図７に基づいて、上述の軌道状態推定手法に基づき軌道状態を推定する軌道状態推定装置１００の構成について説明する。図７は、本実施形態に係る軌道状態推定装置１００の構成を示す機能ブロック図である。軌道状態推定装置１００は、図７に示すように、傾斜角推定部１１０と、軌道状態算出部１２０とを有する。

【0064】

傾斜角推定部１１０は、車両１の軸箱２７に設置された慣性計測装置３０により計測された３軸加速度ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ及び３軸角速度ω_ｘｋ、ω_ｙｋ、ω_ｚｋに基づいて、輪軸２５の水平傾斜角（ロール角θ_ｋ）及び前後傾斜角（ピッチ角φ_ｋ）を推定する。

【0065】

まず、傾斜角推定部１１０は、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから重力加速度射影成分と並進加速度成分とを抽出する。また、傾斜角推定部１１０は、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから抽出された並進加速度成分から、並進加速度を表す指標Ｒ_ｘ、Ｒ_ｙを算出し、上記式（１１）を用いて時刻ｋにおける観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）を設定する。そして、傾斜角推定部１１０は、上記式（７）において、設定した分散共分散行列Ｒ（ｋ）を用いて、観測ベクトルｙ_ｋ（ロール角の推定値θ_０ｋ及びピッチ角の推定値φ_０ｋ）を算出する。さらに、傾斜角推定部１１０は、式（５）～式（７）で表され、式（１１）により観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）が設定された状態空間モデルに基づくカルマンフィルタにより、ロール角θ_ｋ及びピッチ角φ_ｋを算出する。傾斜角推定部１１０は、算出した輪軸２５のロール角θ_ｋ及びピッチ角φ_ｋを、軌道状態算出部１２０へ出力する。

【0066】

軌道状態算出部１２０は、輪軸２５の水平傾斜角（ロール角θ_ｋ）及び前後傾斜角（ピッチ角φ_ｋ）に基づいて、軌道状態を表す軌道状態情報を算出する。本実施形態では、軌道状態情報として、平面性変位が算出される。平面性変位を軌道状態情報として求める場合、軌道状態算出部１２０は、傾斜角推定部１１０により算出された輪軸２５のロール角θ_ｋに基づき、上記式（１）及び式（２）を用いて、平面性変位Ｔｗ（ｘ）を算出する。軌道状態算出部１２０は、算出した平面性変位Ｔｗ（ｘ）を、例えば出力装置２００へ出力する。出力装置２００は、例えばディスプレイ等の表示装置であってもよい。

【0067】

本実施形態に係る軌道状態推定装置１００は、例えばＣＰＵ（Central Processing Unit）やＤＳＰ（Digital Signal Processor）等の各種のプロセッサによって構成され、当該プロセッサが所定のプログラムにしたがって動作されることにより上記機能を実現し得る。

【0068】

本実施形態に係る軌道状態推定装置１００は、車両１に設置されてもよく、別途の場所に設置されていてもよい。慣性計測装置３０による３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ及び３軸角速度の計測値ω_ｘｋ、ω_ｙｋ、ω_ｚｋは、位置検出装置によって検出されたレール５上の車両１の絶対位置に関連付けて、記憶装置（図示せず。）に記録されている。軌道状態推定装置１００は、記憶装置に記録された３軸加速度の計測値ａ_ｙｋ、ａ_ｚｋ、３軸角速度の計測値ω_ｘｋ、ω_ｙｋ、ω_ｚｋ、及び、車両１の絶対位置を取得することにより、軌道状態を推定し得る。

【0069】

［２－３．軌道状態推定方法］
図８に基づいて、本実施形態に係る軌道状態推定装置１００により実施される軌道状態推定方法を説明する。図８は、本実施形態に係る軌道状態推定方法の一例を示すフローチャートである。

【0070】

（Ｓ１０、Ｓ２０：３軸加速度及び３軸角速度の取得）
まず、車両１の軸箱２７に設置された慣性計測装置３０により、３軸加速度ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ及び３軸角速度ω_ｘｋ、ω_ｙｋ、ω_ｚｋが計測される（Ｓ１０、Ｓ２０）。慣性計測装置３０により計測された３軸加速度及び３軸角速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋ、ω_ｘｋ、ω_ｙｋ、ω_ｚｋは、軌道状態推定装置１００の傾斜角推定部１１０へ出力される。

【0071】

（Ｓ３０：重力加速度射影成分の抽出）
次いで、傾斜角推定部１１０は、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから、重力加速度射影成分を抽出する（Ｓ３０）。３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの重力加速度射影成分は、加速度低周波成分であり、例えば０～１Ｈｚの成分とする。傾斜角推定部１１０は、例えばローパスフィルタにより、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから重力加速度射影成分を抽出し得る。

【0072】

（Ｓ４０：並進加速度成分の抽出）
また、傾斜角推定部１１０は、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから、並進加速度成分を抽出する（Ｓ４０）。３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの並進加速度成分は、加速度高周波成分であり、例えば１～７Ｈｚの成分とする。傾斜角推定部１１０は、例えばバンドパスフィルタにより、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋから並進加速度成分を抽出し得る。上述したように、平面性変位について、車両１の曲線走行時に卓越する衝撃に起因した加速度成分である並進加速度成分が、ロール角の推定値θ_ｋ及びピッチ角の推定値φ_ｋの推定精度に大きく影響しているとの知見が得られている。そこで、傾斜角推定部１１０は、並進加速度成分を抽出し、並進加速度を表す指標Ｒ_ｘ、Ｒ_ｙを算出する。

【0073】

（Ｓ５０：傾斜角の推定）
次いで、傾斜角推定部１１０は、輪軸２５の傾斜角θ（ここでは、ロール角θｘ）を推定する（Ｓ５０）。輪軸２５のロール角θｘは、上記式（５）～式（７）で表される状態空間モデルに基づくカルマンフィルタを用いて算出し得る。ここで、傾斜角推定部１１０は、ステップＳ４０にて算出した指標Ｒ_ｘ、Ｒ_ｙを用いて、上記式（１１）から時刻ｋにおける観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）を設定する。これにより、上記式（７）に示した観測方程式の観測ノイズｗ_ｋが、並進加速度成分の大きさに応じて逐次的に設定される。その結果、慣性計測装置３０を軸箱２７に設けることによって計測値に含まれる、車両走行時の並進加速度成分による誤差影響が低減され、観測ベクトルｙ_ｋ（ロール角の推定値θ_０ｋ及びピッチ角の推定値φ_０ｋ）を精度よく求めることができる。

【0074】

傾斜角推定部１１０は、式（５）～式（７）で表され、式（１１）により観測ノイズｗ_ｋの分散共分散行列Ｒ（ｋ）が設定された状態空間モデルに基づくカルマンフィルタにより、ロール角θ_ｋ及びピッチ角φ_ｋを算出する。そして、傾斜角推定部１１０は、算出した輪軸２５のロール角θ_ｋ及びピッチ角φ_ｋを、軌道状態算出部１２０へ出力する。

【0075】

（Ｓ６０：軌道状態情報の推定）
その後、軌道状態算出部１２０は、輪軸２５の水平傾斜角（ロール角θ_ｋ）及び前後傾斜角（ピッチ角φ_ｋ）に基づいて、軌道状態を表す軌道状態情報を算出する（Ｓ６０）。本実施形態では、軌道状態情報として、平面性変位が算出される。平面性変位を軌道状態情報として求める場合、軌道状態算出部１２０は、ステップＳ５０にて算出された輪軸２５のロール角θ_ｋに基づき、上記式（１）及び式（２）を用いて、平面性変位Ｔｗ（ｘ）を算出する。軌道状態算出部１２０は、算出した平面性変位Ｔｗ（ｘ）を、出力装置２００へ出力する。

【0076】

［３．まとめ］
以上、本発明の一実施形態に係る車両と、軌道状態として平面性変位を推定する軌道状態推定装置及び軌道状態推定方法とについて説明した。本実施形態によれば、車両の軸箱に慣性計測装置を設置し、当該慣性計測装置により計測された３軸加速度及び３軸角速度に基づいて、輪軸の水平傾斜角及び前後傾斜角を推定する。輪軸の水平傾斜角及び前後傾斜角はカルマンフィルタにより推定し得る。この際、慣性計測装置の計測値に含まれる車両走行時の並進加速度成分による誤差影響を低減するため、車両走行時の３軸加速度の計測値から誤差因子となる並進加速度成分をフィルタリング処理により抽出し、抽出された並進加速度成分に基づいてカルマンフィルタの観測誤差共分散行列を逐次的に設定する。これにより、輪軸の水平傾斜角及び前後傾斜角を精度よく推定することができる。高精度に輪軸の水平傾斜角及び前後傾斜角を推定することで、平面性変位の推定精度を高めることができる。また、車両を走行させれば３軸加速度及び３軸角速度を測定することができるため、簡易に検測することができる。

【実施例0077】

上記実施形態に係る軌道状態推定手法の有効性を検証するため、当該軌道状態推定手法による平面性変位の推定値を評価した。検証に使用した車両は、図１に示した構成の車両であり、車両進行方向（Ｘ方向）最前の一対の車輪の軸箱にそれぞれ慣性計測装置を設けた。軌道条件として、走行区間は、一部に最大４０ｍｍ程度の平面性変位が設けられた全長２００ｍの曲線区間（Ｒ＝１００ｍ）を設定した。

【0078】

本検証では、手押し式の軌道検測装置により測定された平面性測定値を真値（Ｔｒｕｅ）とした。実施例として、上記実施形態に係る軌道状態推定手法を用いて平面性測定値を推定した。このとき、上記式（１１）における係数Ｋ_ｙ、Ｋ_ｘについて、平面性変位の推定精度が最も高くなる値を探索した。その結果、係数Ｋ_ｙ、Ｋ_ｘはともに１．０×１０^－６となり、かかる係数Ｋ_ｙ、Ｋ_ｘに基づき分散共分散行列Ｒを設定した。プロセスノイズＱは上記式（１２）とした。

【0079】

比較例として、上記式（７）に示した観測方程式の観測ノイズｗ_ｋを固定値として観測ベクトルｙ_ｋを算出し、平面性変位を算出した。すなわち、比較例では、車両の曲線走行時に卓越する衝撃に起因した加速度成分（並進加速度成分）がノイズとして残っている状態で、平面性変位が推定される。比較例において、分散共分散行列Ｒ及びロセスノイズＱは下記の通りとした。

【0080】

【数7】

【0081】

図９に、実施例及び比較例において推定された平面性変位を示す。図９に示すように、比較例では、推定された平面性変位は真値からずれが生じていたが、実施例では、推定された平面性変位は真値とほぼ同値となった。これより、上記実施形態に係る軌道状態推定手法を適用することで、車両走行時の衝撃応答による影響をはじめとする並進加速度による誤差影響が低減され、平面性変位の推定精度を高めることができることが示された。

【0082】

以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。

【0083】

例えば、上記実施形態では、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの重力加速度射影成分である加速度低周波成分を０～１Ｈｚの成分とし、３軸加速度の計測値ａ_ｘｋ、ａ_ｙｋ、ａ_ｚｋの並進加速度成分である加速度高周波成分を１～７Ｈｚの成分としたが、本発明はかかる例に限定されない。上記加速度低周波成分及び加速度高周波成分の周波数帯域は、例えば、製鉄所において重量の大きい資材等を搬送するために使用される、低速走行車両を想定して設定されたものである。加速度低周波成分及び加速度高周波成分の周波数帯域は、例えば上記式（８）及び式（９）に基づき、車両の走行速度に応じて適宜設定すればよい。指標Ｒ_ｘ、Ｒ_ｙの移動平均サンプル数ｓについても、上記式（１０）に基づき設定すればよい。

【符号の説明】

【0084】

１ 車両
５ レール
１０ 車体部
２０ 台車部
２１ 台車枠
２３ 車輪
２５ 輪軸
２７ 軸箱
２９ コイルばね
３０ 慣性計測装置
１００ 軌道状態推定装置
１１０ 傾斜角推定部
１２０ 軌道状態算出部
２００ 出力装置

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】