(57)【要約】【課題】本発明は、単純化したダミー人形と車両のモデルを用いて、入力パラメータの数は少なく演算時間も短いにもかかわらず、車両後面からの衝突時にダミー人形の頚部が受ける傷害値指標値を、精度良く算出することが可能な、解析方法を提供する。【解決手段】本発明の解析方法は、乗員を模擬するダミー人形を用いた車両の後面衝突試験時の、シートに着座したダミー人形の挙動を、車両の前後方向の運動のみに単純化して、胸部と頭部の２自由度系バネマス非減衰モデルで簡易的に表現して運動方程式を立て、前記運動方程式は、車両座席のヘッドレストと、ダミー人形頭部と、が接触を開始するする時刻、即ちコンタクトタイムの前と、後と、についてそれぞれ用意し、前記２つの運動方程式を解くことにより、胸部と頭部の変位と、速度と、加速度と、を算出し、これらの解を公式に代入することによって、頚部障害値指標であるＮＩＣ値（Neck Injury Criterion）を算出している。【選択図】図１
2015087945.tif 000002
【特許請求の範囲】
【請求項１】
車両のシートに載せられるダミー人形の力学モデルを運動方程式により数式化して、車両衝突時のダミー人形の挙動を算出する車両衝突解析方法において、
前記力学モデルは２質点のバネマスダンパモデルとして構成された前記ダミー人形の
胸部と頭部の力学モデルを含むとともに、
乗員を模擬するダミー人形を用いた車両の後面衝突で試験に用いられる、定められた三角波の加速度の過渡入力が加わった時の、シートに着座した前記ダミー人形の挙動を、車両の前後方向の運動のみに単純化して、胸部と頭部の２自由度系バネマス非減衰モデルで簡易的に表現して運動方程式を立て、
前記運動方程式は、車両座席のヘッドレストと、ダミー人形頭部と、が接触を開始するする時刻、即ちコンタクトタイムの前と、後と、についてそれぞれ用意され、
胸部と頭部について立てた前記２つの運動方程式を、コンタクトタイムの前と、後と、各々について解くことにより、胸部と頭部の変位と、速度と、加速度と、を算出し、これらの解を公式に代入することによって、頚部障害値指標であるNIC値（Neck Injury Criterion）を算出する、
ことを特徴とする解析手法。
【請求項２】
請求項１の解析手法において、
ヘッドレストとダミー人形頭部との前後方向クリアランスは、後面衝突スレッド試験前の前記クリアランスに、車両座席のシートバックの動的変形によりヘッドレストが車両後方へ移動する量を加えて補正した値を算出し、
ダミー人形頭部がヘッドレストに衝突する時刻、即ちコンタクトタイムは、後面衝突スレッド試験での頭部加速度発生結果を参考の上、誤差を補正した値を算出し、
前記の補正された値を運動方程式に入力する、
ことを特徴とする解析手法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、乗員を模擬するダミー人形がシートを介して搭載された車両の、後面衝突時におけるダミー人形の挙動の解析方法に関する。
【背景技術】
【０００２】
これまで車室内乗員の衝突時の挙動や傷害値計算手法として、車両の全面衝突時の膝荷重や胸部が受ける加速度を、バネマスモデルを用いて予測する手法等は多く考えられてきた。
【０００３】
車両のシートに載せられた乗員を模擬するダミー人形を用いた車両衝突シミュレーションとして、例えば特許文献１に記載された技術が知られている。この車両衝突シミュレーション方法では、車両のシートに載せられたダミー人形と、当該ダミー人形を前記シートに拘束するシートベルトを含む拘束体と、から構成される力学モデルを、バネ特性を用いて比較的単純な力学モデルとしてモデル化する。そして、このモデルの運動方程式に対してダミー・スレッド試験（ダミー人形を用いた車両衝突試験）を通じて同定された上記の力学モデルのバネ特性を代入すると共に、車両に衝突減速特性を与えることにより、この運動方程式を解いて、車両衝突時のダミー人形の挙動を求めている。
【０００４】
詳述すれば、ダミー人形に関する、胸部と腰部と、の２質点で構成された力学モデルと、ダミー人形をシートに拘束する少なくともシートベルトを含む拘束体に関する、所定のバネ特性を付与される腰ベルトと、フォースリミッタ付き肩ベルトと、ニーボルスタ（衝突時に膝を支持し、ダミー人形の前方移動を抑制するもの）と、で力学モデルを構成している。尚、肩ベルトに対しては、フォースリミッタ作動前と作動後で、それぞれ異なった荷重特性を設定し、ニーボルスタに対しては、ソフトクッション領域と、非ソフトクッション領域と、に分けて、バネ特性を設定している。これにより、衝突傷害評価において重要な胸部と腰部の挙動に関して、これらの挙動に大きな影響を与える腰ベルトと肩ベルトとニーボルスタとが正確にモデル化され、胸部と腰部との力学的連係も十分に考慮されるので、制度良くシミュレーションが可能なものとなっている。
【０００５】
さらに、上記の力学モデルを運動方程式化する際、胸部のみの単独運動とした胸部運動方程式と、腰部のみの単独運動方程式とした腰部運動方程式とを作成し、胸部は腰部を中心とした回転運動と見なして、最終的な胸部の減速度を補正することも提案されている。
【０００６】
しかしながら、車両のシートに載せられた乗員を模擬するダミー人形を用いた車両衝突シミュレーションのうち、車両後面からの衝突時の胸部と腰部との運動差で生じる鞭打ち現象については、衝突中にシートバックがダミー人形からの慣性荷重を受けて略車両後方へ変形し、それに伴ってヘッドレストが後方へ移動し、頭部とヘッドレストとの距離が遷移する等の現象が発生するため、頭部とヘッドレストの衝突の表現には多自由度の複雑なモデルが必要になり、解析事例は、各質点間を弾性特性や各種ジョイント構造で接合して人体の解析を行う衝突解析ソフトウエア（例えばオランダ応用科学研究機構製 MADYMO）等の大規模解析を除いてほとんど報告されていない。
【０００７】
一方で、前記衝突解析ソフトウエアMADYMOを用いて、交通事故傷害の解析を行なったり、さらにはそれを発展させて、人体全身の有限要素モデルと車両構成部材の有限要素モデルを構成し、車両衝突時のシミュレーション結果から人体傷害を再現する等の技術が知られている。しかしながら、有限要素法を用いたシミュレーションの場合、車両の詳細な寸法レイアウトが要求されるのに対して、実際には試作車両などではそのようなデータを正確に与えることが困難である。また、仮にそのようなデータを与えたとしても、計算の前準備および計算そのものに時間およびコストがかかり、車体設計段階で頻繁に利用することは事実上不可能に近い。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００５−１１５７７０号公報（特許第４１４１３６７号広報）                  （〔０００９〕〜〔００１２〕段落、図２）
【特許文献２】特開２００２−１４９７１９号公報（〔００１８〕〜〔００１９〕段落、図１）
【発明の概要】
【発明が解決しようとする課題】
【０００９】
以上の問題に鑑みて本発明は創作されたものであり、本発明は、単純化したダミー人形と車両のモデルを用いて、入力パラメータの数は少なく演算時間も短いにもかかわらず、車両後面からの衝突時にダミー人形の頚部が受ける傷害値指標値を、精度良く算出することが可能な、解析方法の提供を目的とした。
【課題を解決するための手段】
【００１０】
上述した課題を解決すべく提供される本発明の解析方法は、車両のシートに載せられるダミー人形の力学モデルを運動方程式により数式化して、車両衝突時のダミー人形の挙動を算出する車両衝突解析方法において、前記力学モデルは２質点のバネマスダンパモデルとして構成された前記ダミー人形の胸部と頭部の力学モデルを含むとともに、乗員を模擬するダミー人形を用いた車両の後面衝突で試験に用いられる、定められた三角波の加速度の過渡入力が加わった時の、シートに着座した前記ダミー人形の挙動を、車両の前後方向の運動のみに単純化して、胸部と頭部の２自由度系バネマス非減衰モデルで簡易的に表現して運動方程式を立て、前記運動方程式は、車両座席のヘッドレストと、ダミー人形頭部と、が接触を開始するする時刻、即ちコンタクトタイムの前と、後と、について各々用意され、胸部と頭部について立てた前記２つの運動方程式を、コンタクトタイムの前と、後と、各々について解くことにより、胸部と頭部の、変位と、速度と、加速度と、を算出し、これらの解を公式に代入することによって、頚部障害値指標であるNIC値（Neck Injury Criterion）を算出することを特徴とするものである。
【００１１】
本発明においては、車両のシートに載せられるダミー人形について、後面衝突時における胸部と頭部の、変位と、速度と、加速度と、を算出可能な力学モデルを構成する。すなわち、この力学モデルを運動方程式により数式化した場合、その変数として胸部と頭部の、変位と、速度と、加速度と、を含むように力学モデルを構成する。そして、この運動方程式を解くことにより求めた、前記胸部と頭部の、変位と、速度と、加速度と、をNIC値を算出するための公式に代入することにより、後面衝突時のダミー人形のNIC値が算出される。
【００１２】
上述した構成とすることにより、入力パラメータの数が少ないにもかかわらず、後面衝突時におけるNIC値を精度良く算出可能な解析方法を提供することが可能となる。
【００１３】
また、前記解析方法の力学モデルにおいて使用される、ヘッドレストとダミー人形頭部との前後方向クリアランスは、後面衝突スレッド試験前の前記クリアランスに、車両座席のシートバックの動的変形によりヘッドレストが車両後方へ移動する量を加えて補正した値を算出し、また、ダミー人形頭部がヘッドレストに衝突する時刻、即ちコンタクトタイムは、後面衝突スレッド試験での頭部加速度発生結果を参考の上、誤差を補正した値を算出した値を算出し、前記の補正された値を前記運動方程式に入力することが望ましい。結果として、現象の解析精度が向上し、良好な解析結果を得ることができる。したがって、車両用シートの開発時における後面衝突鞭打ち性能の、開発初期段階での構造決定に際して有効なツールとなり、開発期間の短縮や、開発費の節減にも効力を発生する。
【００１４】
このように本発明は、運動方程式を、ダミー人形頭部がヘッドレストに衝突する時刻、即ちコンタクトタイムの前用と、後用と、の２つを用意して切り替えて使用し、さらに運動方程式の入力値として、ヘッドレストとダミー人形頭部との前後方向クリアランスと前記コンタクトタイムとの２項目について補正値を用いることにより、ダミー人形の後面衝突時のNIC値を精度良く算出できるという本発明者の知見に基づくものである。また、単純な力学モデルを用いることにより計算量を減少させることが可能となり、小規模なコンピュータシステムであっても迅速に計算でき、車両用シートの開発時における設計変更にも柔軟に対応し、解析することができ、開発効率を向上できる。
【発明の効果】
【００１５】
本発明によれば、単純化したダミー人形と車両シートのモデルを用いて、入力パラメータの数は少なく演算時間も短いにもかかわらず、車両後面からの衝突時にダミー人形の頚部が受ける傷害値指標値を、精度良く算出することが可能な車両の後面衝突時の解析方法を提供できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の車両衝突解析方法における力学モデルの説明図である。
【図２】後面衝突スレッド試験用装置の略図である。
【図３】後面衝突スレッド試験において車両に加える、定められた三角波の加速度特性を示す図である。
【図４】後面衝突スレッド試験におけるダミー胸部T1(第1胸椎) の前後方向荷重撓み曲線の測定結果、およびバネ定数：k1を示す直線を示すグラフである。
【図５】後面衝突スレッド試験におけるダミー頭部の前後方向荷重撓み曲線の測定結果、およびバネ定数：k3を示す直線を示すグラフである。
【図６】インパクタによるダミーのキャリブレーションテストの原理図である。
【図７】ダミー頚部バネ定数：k2を求める畳み込み積分手法を説明する図である。
【図８】計算に用いた胸部T1加速度入力値特性を示すグラフである。
【図９】ダミー頚部バネ定数：k2の算出結果と実験結果を示すグラフである。
【図１０】NIC値の解析に用いられる定められた三角波の加速度特性と、ダミーを用いた後面衝突スレッド試験で再現した実際の加速度と、を比較したグラフである。
【図１１】入力パラメータの一覧図である。
【図１２】試験結果と計算結果とによる胸部加速度特性グラフである。
【図１３】試験結果と計算結果とによる頭部加速度特性グラフである。
【図１４】試験結果と計算結果とによるNIC値の特性グラフである。
【発明を実施するための形態】
【００１７】
以下、本発明の実施形態を図面に基づいて説明する。
〔力学モデルの構成〕
図１〜図２に、後面衝突を想定した車両１のシート２に載せられたダミー人形３（以下、単に“ダミー”と称する）のバネマスダンパモデル（力学モデルの一例）の説明図を示す。このモデルでは、ダミー胸部の質点（胸部：質量m1）とダミー頭部の質点（頭部：質量m2）との２質点でモデル化している。このようにダミーの胸部と頭部をモデル化することで、胸部動的変位量（地上基準）χ1と、頭部動的変位量（地上基準）χ2と、を各々定義することができる。また、図３は本実施形態における車両１に加える、定められた三角波の加速度（車両１が後面から衝突を受けた際に車両１に加わる加速度）の特性を示す図である。
【００１８】
〔力学モデルの運動方程式〕
以下に、上記モデルを用いた際の運動方程式を説明する。尚、本実施形態において使用する記号は以下の通りである。
Y：車体動的変位量（地上基準）
χ0：ヘッドレストとダミー頭部の前後方向クリアランス
t１：ヘッドレストとダミー頭部のコンタクトタイム
χ1：胸部動的変位量（地上基準）
χ2：頭部動的変位量（地上基準）
m1：ダミー胸部質量
m2：ダミー頭部質量
k1：シートバック部バネ定数
k2：ダミー頚部バネ定数
k3：ヘッドレスト部バネ定数（シートバック部の影響も含む）
a0：車体入力加速度のピーク値
T：車体入力加速度の入力終了時間
kT：車体入力加速度のピーク値発生時間
X1：胸部の自由振動解の振幅
X2：頭部の自由振動解の振幅
T1：ダミー胸部の第１胸椎
【００１９】
まず、m2に関してはt１を境にしてバネマスモデルの形態が変わり、下記の運動方程式が得られる。
【００２０】
【数１】
2015087945.tif 000003
ここで車体入力加速度を積分し、境界条件を用いて強制速度、および強制変位を求めると、下式のようになる。
【００２１】
【数２】
2015087945.tif 000004
χ1,χ2の強制振動解は、式(1)〜(3)に式(10)〜(12)を代入し、微分方程式を解くことにより得られる。
χ1,χ2の強制振動解を、おのおのχfo1,χfo2とすると、
【００２２】
【数３】
2015087945.tif 000005
一方、自由振動解は、式(1)，(2)，(3)の中の強制変位成分を0とすることにより得られる。
χ1,χ2の自由振動解を、おのおのχfr1,χfr2とすると、
【００２３】
【数４】
2015087945.tif 000006
χfr1,χfr2を
【００２４】
【数５】
2015087945.tif 000007
とおき、式(21)，(22)，(23)に代入して整理すると、
【００２５】
【数６】
2015087945.tif 000008
X1，X2が共に0にならない解を得るためには、
【００２６】
【数７】
2015087945.tif 000009
これを展開してω2について解き、2根の小さい方をω12，大きい方をω22とすると、ωは±ω1と±ω2の4個の値をとるが、＋ω1と＋ω2の2個が意味を持ち、
【００２７】
【数８】
2015087945.tif 000010
2自由度系振動としては、m1とm2がω1で振動し、それぞれの振幅がX1，X2の場合と、m1とm2がω2で振動し、それぞれが別の振幅X1，X2で振動する2形態が存在する。式  (21)，(24)，(25)より
【００２８】
【数９】
2015087945.tif 000011
ω1に対するX1，X2，χfr1,χfr2をX1(1)，X2(1)，χfr1(1),χfr2(1)、ω2に対するX1，X2，χfr1,χfr2をX1(2)，X2(2)，χfr1(2),χfr2(2)とし、
【００２９】
【数１０】
2015087945.tif 000012
とおく。式(35)、(36)を式(24)、(25)に代入すると、
【００３０】
【数１１】
2015087945.tif 000013
【００３１】
式(37)、(38)は同位相、式(39)、(40)は逆位相で振動する。式(21)、(22)、(23)の微分方程式の解は、α、βを定数とすると、
【００３２】
【数１２】
2015087945.tif 000014
式(41)を置き換えると、
【００３３】
【数１３】
2015087945.tif 000015
これを展開すると、
【００３４】
【数１４】
2015087945.tif 000016
ここで、
【００３５】
【数１５】
2015087945.tif 000017
と置くと、式(44)、(45)は、
【００３６】
【数１６】
2015087945.tif 000018
E1〜E4は各時間帯にて異なった係数となるので、係数を変えた上で各時間帯別の自由振動解と強制振動解を加算し、全体の振動解を求めると、
【００３７】
【数１７】
2015087945.tif 000019
ここで式(52)〜(59)におけるχfo1,χfo2の値は、式(13)〜(20)の各時間帯に対応した値となる。また、式(52)〜(59)の各係数は未知数となり、合計16個存在するが、境界条件として、式(52)、(53)において、
【００３８】
【数１８】
2015087945.tif 000020
と、分割した車体入力加速度の隣り合う時間における
【数１９】
2015087945.tif 000021
を等しくし、連立方程式を解けば各係数が求まるが、ここでは紙面の関係により計算結果の記載を省略する。
【００３９】
以上によりχ1,χ2が求まったため
【数２０】
2015087945.tif 000022
も求まり、NIC値に関しても式(61)にて算出できる。
【００４０】
【数２１】
2015087945.tif 000023
また、ヘッドレストとダミー頭部のコンタクトタイムであるt１は、次式を満足した時間に一致する。
【００４１】
【数２２】
2015087945.tif 000024
【００４２】
〔計算に用いる入力値の選定〕
この衝突解析を行う際の入力パラメータの中でバネ定数：k1〜k3、ヘッドレストとダミー頭部のコンタクトタイム：t１に関しては、あらかじめ実験を行い、値を求めておく必要がある。
ここでk1、k3はシートバックやヘッドレストへの静的、又は動的荷重負荷試験により求めるべきものであるが、これらは実際の後面衝突スレッド試験により算出する。図４と図５に、後面衝突スレッド試験におけるダミー胸部T1(第1胸椎)と、頭部と、の前後方向荷重撓み曲線の測定結果、およびバネ定数：k1〜k3を示す直線が示されている。
【００４３】
また、ダミー頚部の前後方向バネ定数であるk2は、図６に示す、インパクタによるダミーのキャリブレーションテストにより算出する。前記テストでは、インパクタでダミーの背面に衝撃を加えたときのダミー胸部T1前後方向加速度とダミー頭部前後方向加速度が得られるため、図７に示す畳み込み積分手法によりk2を求める。図８に計算に用いた胸部T1加速度入力値、図９にk2の算出結果が示されている。
【００４４】
一方、ヘッドレストとダミー頭部のコンタクトタイム：t１は、本来ヘッドレストとダミー頭部のコンタクトタイムであるが，前記スレッド試験結果では、頭部がヘッドレストにコンタクト時発生する反発現象の影響を受けて実質的な頭部拘束開始時間がコンタクトタイムより遅れて発生することに着目し、前記スレッド試験での頭部加速度発生結果を参考の上、決定する。
【００４５】
また、χ0は本来、後面衝突スレッド試験前のダミー頭部とヘッドレストまでの前後方向のクリアランスであるが、頭部がヘッドレストとコンタクトする前にヘッドレストがシートバックの変形により後方へ移動することを考慮し、動的変形による移動量も加えた値とし、式(62)により算出する。
【００４６】
〔実験値と解析計算値との比較〕
本発明による簡易的な力学モデルによる解析による結果が、ダミーを用いた後面衝突スレッド試験による結果にどの程度合致しているかが調べられた。NIC値の解析に用いられる定められた三角波の加速度特性と、前記スレッド試験で再現した実際の加速度と、の比較グラフが図１０に示されている。この図から、スレッド試験における加速度特性は、解析に用いられる加速度特性とほぼ同等の波形が再現されていることが理解できる。
【００４７】
この比較テストの際、解析に用いられた入力パラメータが図１１に一覧化されている。このようなパラメータを、前述した式に代入して得られた、胸部と、頭部と、の加速度波形およびNIC値の計算結果と、実際にスレッド試験によって得られた試験結果と、計算結果と、が図１２〜図１４に示されている。図１２は時間に対する胸部加速度の変化を、図１３は時間に対する胸部加速度の変化を、図１４は時間に対するNIC値の変化を、それぞれ示している。これらの計算結果と試験結果とを比較した図から、本発明による後面衝突解析が、非常に簡易なモデル化から出発しているにもかかわらず、スレッド試験による試験結果にかなり類似していることが分かる。特にNIC値については、0〜150msまでの間のピーク値が採択されるが、計算値と実験値との結果比較では、ピーク値のみならずピーク値発生時間やNIC値の時間毎のグラフ推移もほぼ整合しており、この後面衝突解技術がシート設計時等での机上の頚部障害値指標の推測に十分利用できることを証明している。
【産業上の利用可能性】
【００４８】
本発明による車両後面衝突時の解析技術は、シートに載った乗員の、後面衝突事故における挙動を解析して乗員傷害の低減を図ろうとする自動車のみならず、鉄道やその他乗り物全般におけるシート設計の分野に応用することができる。
【００４９】
以上、本発明の解析方法についての実施形態について説明してきたが、本発明はこれに限定されるものではなく、特許請求の範囲および明細書等に記載の精神や教示を逸脱しない範囲で他の変形例、改良例が得られることが当業者は理解できるであろう。
【符号の説明】
【００５０】
１  車両
２  シート
３  ダミー人形
３a  胸部要素
３b  頭部要素
【図１】
2015087945.tif 000025
【図２】
2015087945.tif 000026
【図３】
2015087945.tif 000027
【図４】
2015087945.tif 000028
【図５】
2015087945.tif 000029
【図６】
2015087945.tif 000030
【図７】
2015087945.tif 000031
【図８】
2015087945.tif 000032
【図９】
2015087945.tif 000033
【図１０】
2015087945.tif 000034
【図１１】
2015087945.tif 000035
【図１２】
2015087945.tif 000036
【図１３】
2015087945.tif 000037
【図１４】
2015087945.tif 000038
【手続補正書】
【提出日】平成25年11月21日(2013.11.21)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００２５
【補正方法】変更
【補正の内容】
【００２５】
2015087945.tif 000039
X1，X2が共に0にならない解を得るためには、