(57)【特許請求の範囲】
【請求項１】
剛体壁と衝突する車両の時間進行シミュレーションにおいて詳細な剛体壁力一覧を提供する方法であって、
有限要素解析（ＦＥＡ）アプリケーションモジュールがインストールされたコンピュータシステムが、ＦＥＡモデルおよび剛体壁定義を受け取るステップであって、前記ＦＥＡモデルは、車両と前記剛体壁定義とを表しており、前記ＦＥＡモデルおよび前記剛体壁定義は、前記剛体壁と衝突する前記車両の時間進行シミュレーションにおいて用いられ、前記ＦＥＡモデルは、１以上のグループを編成する複数の有限要素によって接続される複数のノードを有し、前記剛体壁定義は、前記剛体壁上の対応する空間的位置を表す１以上のセグメントを有するステップと、
詳細な剛体壁力（ＲＷＦ）一覧を得ることが所望される前記グループのリストが定義されているか否かを判断するステップと、
前記リストが定義されていない場合に、前記車両の全ての構造部品を含めるよう前記リストを生成するステップと、
前記リスト内の１以上の前記グループに対して、各前記ノードについて１以上の寄与重み係数をそれぞれ計算するステップと、
時間に関する複数のソリューションサイクルを含む、前記時間進行シミュレーションを開始するステップと、
各前記ソリューションサイクルにおいて、ノードタイプに応じて各前記ノードのノード力寄与を計算し、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、前記セグメントの対応する１つについて、前記グループの対応する１つの詳細な前記ＲＷＦ一覧へと集計するステップと、
前記時間進行シミュレーションの完了後、詳細な前記ＲＷＦ一覧を提供するステップと、
を備える方法。
【請求項２】
請求項１に記載の方法であって、
前記寄与重み係数は、各前記ノードを所定のスキームで共有する前記グループに対するそれぞれの重み係数を含んでいる、
方法。
【請求項３】
請求項２に記載の方法であって、
前記所定のスキームは、各前記ノードを共有する前記グループの数に基づいている、
方法。
【請求項４】
請求項１に記載の方法であって、
前記の、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、詳細な前記ＲＷＦ一覧へと集計するステップは、さらに、
ノード寄与重み係数で修正された、計算された前記ノード力寄与を、各前記セグメントに対し、前記剛体壁と接触する１以上の前記ノードを含む各前記グループに関して合計するステップと、
特定セグメントの剛体壁力を、前記特定セグメントに接触する全ての前記ノードについて、ノード寄与重み係数で修正された、計算された前記ノード力寄与を合計することにより取得するステップと、
前記セグメントの全てについて、前記剛体壁力の全てを合計することにより、総剛体壁力を取得するステップと、
特定グループによって前記剛体壁に与えられる力の合計を、前記特定グループによって前記セグメントの全てに与えられる力を合計することにより取得するステップと、
を含んでいる、
方法。
【請求項５】
請求項１に記載の方法であって、
前記ノードタイプが、各前記ノードが適切なノード質量でモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝ＭNODE（Ａ＋−Ａ−）
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
ＭNODEはスカラーのノード質量であり、
Ａ＋は前記剛体壁との接触後のノード加速度ベクトルを表し、
Ａ−は前記剛体壁との接触前のノード加速度ベクトルを表している、
方法。
【請求項６】
請求項５に記載の方法であって、
摩擦オプションが前記時間進行シミュレーションにおいて有効である場合、前記ノード力寄与は摩擦力を含んでいる、
方法。
【請求項７】
請求項１に記載の方法であって、前記ノードタイプが、ペナルティ法によって、各前記ノードがノード質量なしでモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝Ｆｎ＋Ｆｔ
Ｆｎ＝（ＫΔｕ）ｎ
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
Ｆｎは前記剛体壁に垂直な前記ノード力寄与の法線部分を表し、
Ｋは前記ペナルティ法においてノード質量なしでモデル化されている各前記ノードに割り当てられた剛性を表し、
Δｕは、前記剛体壁を通る、ノード質量なしでモデル化されている各前記ノードのノード貫通を表し、
ｎは前記剛体壁に対する単位法線ベクトルを表し、
Ｆｔは、摩擦オプションが前記時間進行シミュレーションにおいて有効である場合の、前記法線部分に直交する摩擦力を表している、
方法。
【請求項８】
剛体壁と衝突する車両の時間進行シミュレーションにおける詳細な剛体壁力一覧を提供するシステムであって、
有限要素解析法（ＦＥＡ）アプリケーションモジュールに関するコンピュータ可読コードを記憶するメインメモリと、
前記メインメモリに接続される少なくとも１つのプロセッサと、
を備えるシステムであり、
前記少なくとも１のプロセッサは、前記メインメモリ内の前記コンピュータ可読コードを実行して、前記ＦＥＡアプリケーションモジュールに、方法に基づいてオペレーションを実行させ、
前記方法が、
有限要素解析（ＦＥＡ）モデルおよび剛体壁定義を受け取るステップであって、前記ＦＥＡモデルは車両と前記剛体壁定義とを表しており、前記ＦＥＡモデルおよび前記剛体壁定義は、前記剛体壁と衝突する前記車両の時間進行シミュレーションにおいて用いられ、前記ＦＥＡモデルは、１以上のグループを編成する複数の有限要素によって接続される複数のノードを有し、前記剛体壁定義は、前記剛体壁上の対応する空間的位置を表す１以上のセグメントを有するステップと、
詳細な剛体壁力（ＲＷＦ）一覧を得ることが所望されている前記グループのリストが定義されているか否かを判断するステップと、
前記リストが定義されていない場合に、前記車両の全ての構造部品を含めるよう前記リストを生成するステップと、
前記リスト内の１以上の前記グループに対して、各前記ノードについて１以上の寄与重み係数をそれぞれ計算するステップと、
時間に関する複数のソリューションサイクルを含む、時間進行シミュレーションを開始するステップと、
各前記ソリューションサイクルにおいて、ノードタイプに応じて各前記ノードのノード力寄与を計算し、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、前記セグメントの対応する１つについて、前記グループの対応する１つの詳細なＲＷＦ一覧へと集計するステップと、
前記時間進行シミュレーションの完了後、詳細な前記ＲＷＦ一覧を提供するステップと、
を有するシステム。
【請求項９】
請求項８に記載のシステムであって、
前記寄与重み係数は、各前記ノードを所定のスキームで共有する前記グループに対するそれぞれの重み係数を含んでいる、
システム。
【請求項１０】
請求項９に記載のシステムであって、
前記所定のスキームは、各前記ノードを共有する前記グループの数に基づいている、
システム。
【請求項１１】
請求項８に記載のシステムであって、
前記の、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、詳細な前記ＲＷＦ一覧へと集計するステップは、さらに、
ノード寄与重み係数で修正された、計算された前記ノード力寄与を、各前記セグメントに対し、前記剛体壁と接触する１以上の前記ノードを含む各前記グループに関して合計するステップと、
特定セグメントの剛体壁力を、前記特定セグメントに接触する全ての前記ノードについて、ノード寄与重み係数で修正された、計算された前記ノード力寄与を合計することにより取得するステップと、
前記セグメントの全てについて、前記剛体壁力の全てを合計することにより、総剛体壁力を取得するステップと、
特定グループによって前記剛体壁に与えられる力の合計を、前記特定グループによって前記セグメントの全てに与えられる力を合計することにより取得するステップと、
を含んでいる、
システム。
【請求項１２】
請求項８に記載のシステムであって、前記ノードタイプが、各前記ノードが適切なノード質量でモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝ＭNODE（Ａ＋−Ａ−）
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
ＭNODEはスカラーのノード質量であり、
Ａ＋は前記剛体壁との接触後のノード加速度ベクトルを表し、
Ａ−は前記剛体壁との接触前のノード加速度ベクトルを表している、
システム。
【請求項１３】
請求項８に記載のシステムであって、前記ノードタイプが、ペナルティ法によって、各前記ノードがノード質量なしでモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝Ｆｎ＋Ｆｔ
Ｆｎ＝（ＫΔｕ）ｎ
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
Ｆｎは前記剛体壁に垂直な前記ノード力寄与の法線部分を表し、
Ｋは前記ペナルティ法においてノード質量なしでモデル化されている各前記ノードに割り当てられた剛性を表し、
Δｕは、前記剛体壁を通る、ノード質量なしでモデル化されている各前記ノードのノード貫通を表し、
ｎは前記剛体壁に対する単位法線ベクトルを表し、
Ｆｔは、摩擦オプションが前記時間進行シミュレーションにおいて有効である場合の、前記法線部分に直交する摩擦力を表している、
システム。
【請求項１４】
請求項８に記載のシステムであって、
前記セグメントは前記剛体壁上に配置されるロードセルに対応する、
システム。
【請求項１５】
方法に基づいて、剛体壁と衝突する車両の時間進行シミュレーションにおいて詳細な剛体壁力一覧を提供する命令を有する非一時的なコンピュータ可読記憶媒体であって、
前記方法が、
有限要素解析（ＦＥＡ）アプリケーションモジュールがインストールされたコンピュータシステムが、ＦＥＡモデルおよび剛体壁定義を受け取るステップであって、前記ＦＥＡモデルは、車両と前記剛体壁定義とを表しており、前記ＦＥＡモデルおよび前記剛体壁定義は、前記剛体壁と衝突する前記車両の時間進行シミュレーションにおいて用いられ、前記ＦＥＡモデルは、１以上のグループを編成する複数の有限要素によって接続される複数のノードを有し、前記剛体壁定義は、前記剛体壁上の対応する空間的位置を表す１以上のセグメントを有するステップと、
詳細な剛体壁力（ＲＷＦ）一覧を得ることが所望される前記グループのリストが定義されているか否かを判断するステップと、
前記リストが定義されていない場合に、前記車両の全ての構造部品を含めるよう前記リストを生成するステップと、
前記リスト内の１以上の前記グループに対して、各前記ノードについて１以上の寄与重み係数を計算するステップと、
時間に関する複数のソリューションサイクルを含む、前記時間進行シミュレーションを開始するステップと、
各前記ソリューションサイクルにおいて、ノードタイプに応じて各前記ノードのノード力寄与を計算し、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、前記セグメントの対応する１つについて、前記グループの対応する１つの詳細な前記ＲＷＦ一覧へと集計するステップと、
前記時間進行シミュレーションの完了後、詳細な前記ＲＷＦ一覧を提供するステップと、
を有する、非一時的なコンピュータ可読記憶媒体。
【請求項１６】
請求項１５に記載の非一時的なコンピュータ可読記憶媒体であって、
前記寄与重み係数は、各前記ノードを所定のスキームで共有する前記グループに対するそれぞれの重み係数を含んでいる、
非一時的なコンピュータ可読記憶媒体。
【請求項１７】
請求項１６に記載の非一時的なコンピュータ可読記憶媒体であって、
前記所定のスキームは、各前記ノードを共有する前記グループの数に基づいている、
非一時的なコンピュータ可読記憶媒体。
【請求項１８】
請求項１５に記載の非一時的なコンピュータ可読記憶媒体であって、
前記の、対応する前記寄与重み係数で修正された、計算された前記ノード力寄与を、詳細な前記ＲＷＦ一覧へと収集するステップは、さらに、
ノード寄与重み係数で修正された、計算された前記ノード力寄与を、各前記セグメントに対し、前記剛体壁と接触する１以上の前記ノードを含む各前記グループに関して合計するステップと、
特定セグメントの剛体壁力を、前記特定セグメントに接触する全ての前記ノードについて、ノード寄与重み係数で修正された、計算された前記ノード力寄与を合計することにより取得するステップと、
前記セグメントの全てについて、前記剛体壁力の全てを合計することにより、総剛体壁力を取得するステップと、
特定グループによって前記剛体壁に与えられる力の合計を、前記特定グループによって前記セグメントの全てに与えられる力を合計することにより取得するステップと、
を含んでいる、
非一時的なコンピュータ可読記憶媒体。
【請求項１９】
請求項１５に記載の非一時的なコンピュータ可読記憶媒体であって、前記ノードタイプが、各前記ノードが適切なノード質量でモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝ＭNODE（Ａ＋−Ａ−）
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
ＭNODEはスカラーのノード質量であり、
Ａ＋は前記剛体壁との接触後のノード加速度ベクトルを表し、
Ａ−は前記剛体壁との接触前のノード加速度ベクトルを表している、
非一時的なコンピュータ可読記憶媒体。
【請求項２０】
請求項１５に記載の非一時的なコンピュータ可読記憶媒体であって、前記ノードタイプが、ペナルティ法によって、各前記ノードがノード質量なしでモデル化されていることを示す場合、前記ノード力寄与は、
ＦNODE＝Ｆｎ＋Ｆｔ
Ｆｎ＝（ＫΔｕ）ｎ
で計算され、
ＦNODEは前記ノード力寄与を表すベクトルであり、
Ｆｎは前記剛体壁に垂直な前記ノード力寄与の法線部分を表し、
Ｋは前記ペナルティ法においてノード質量なしでモデル化されている各前記ノードにおいて割り当てられた剛性を表し、
Δｕは、前記剛体壁を通る、ノード質量なしでモデル化されている各前記ノードのノード貫通を表し、
ｎは前記剛体壁に対する単位法線ベクトルを表し、
Ｆｔは、摩擦オプションが前記時間進行シミュレーションにおいて有効である場合の、前記法線部分に直交する摩擦力を表している、
非一時的なコンピュータ可読記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概して、コンピュータ支援工学解析に関する。本発明は、特に、剛体壁と衝突する車両（例えば自動車）の時間進行シミュレーションにおいて詳細な剛体壁力一覧を提供する方法およびシステムに関する。
【背景技術】
【０００２】
コンピュータ支援工学（ＣＡＥ）が、多くのタスクにおいてエンジニアを支援するために用いられている。例えば、構造体あるいは製品の設計の処理手続において、ＣＡＥ解析、特に有限要素解析法（ＦＥＡ）が、種々の荷重条件（例えば、静的あるいは動的な荷重条件）下での応答（例えば応力、変位など）を評価するために、しばしば用いられている。最も重要なコンピューターシミュレーションのうちの１つに、剛体壁と衝突する車両のシミュレーションがある。剛体壁は、静止している場合も、動いている場合もある。一般的に、車両は、空間に固定された剛体壁に向かって移動する。しかしながら、自動車の屋根が押し潰されるシミュレーションの際には、剛体壁が静止している自動車に向かって移動する。このようなシミュレーションは、ユーザが安全性の要件を満たす自動車を設計するのを支援するために利用できる。収集すべき重要な情報の一つは、自動車と剛体壁との衝突によって生じる剛体壁力（ＲＷＦ：Rigid Wall Force）である。剛体壁力一覧の合計は、従来技術アプローチに示されているような数値シミュレーションにおいて提供される。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、よりよい自動車を設計および／または構築するために、自動車技術者は総剛体壁力の空間的分布を知りたいと考えるであろう。物理的な模型衝突試験においては、一般的に、剛体壁のセグメントにそれぞれ取り付けられた多数のロードセルが、空間的力分布（つまり剛体壁上の種々の位置すなわちセグメントにおけるＲＷＦ）を決定するために取り付けられている。しかしながら、この力分布情報からは、各車両構造部品がどのように総剛体壁力に寄与するかに関して詳細な情報が得られない。車両構造は複雑であるため、そのような詳細な情報は、実験では取得できない。設計技術者は、この情報を、構造部品の設計および変更において、車両における荷重経路を決定するために用いる。したがって、剛体壁と衝突する車両の時間進行数値シミュレーションにおいて、詳細な剛体壁力一覧を提供する方法およびシステムが望まれる。
【課題を解決するための手段】
【０００４】
この章は、本発明のいくつかの面を要約することを目的とし、いくつかの好適な実施形態を簡潔に紹介するものである。この章の目的が曖昧になることを避けるため、この章においては、要約およびここでの発明の名称と同様に、簡略化あるいは省略が行われている場合がある。このような簡略化あるいは省略は、本発明の範囲を制限するものではない。
【０００５】
剛体壁と衝突する車両の時間進行シミュレーションにおいて詳細な剛体壁力一覧を提供する方法およびシステムを提供する。本発明の１の態様では、車両の演算モデル（例えば有限要素解析（ＦＥＡ）モデル）および剛体壁定義が、有限要素解析アプリケーションモジュールがインストールされたコンピュータシステムに受け取られる。ＦＥＡモデルは、グループを編成する有限要素によって接続される多数のノードを有する。剛体壁は、剛体壁に取り付けられたロードセルに対応する、１以上のセグメントを有する。詳細な剛体壁力（ＲＷＦ）一覧を得ることが所望されるグループのリスト（つまり、ノードのリスト、有限要素のリスト、車両の構造部品（例えばバンパーまたはバンパーの一部）など）が、ユーザ（つまり、自動車の設計および／または製造に携わる技術者および／または科学者）によって定義される。グループのそれぞれは、任意の数の剛体壁セグメントと相互に作用できる。定義されたグループのリストがない場合、デフォルト（初期設定）リストを生成できる（例えば、リストにＦＥＡモデルの全ての構造部品を含める）。ＦＥＡモデルにおけるそれぞれのノードに対して、寄与重み係数が計算される。ノードがリストの１つのグループのみに属する場合、寄与重み係数は１である（つまり、ノード力のＲＷＦへの寄与は、全てそのグループへのものとなる）。ノードが２以上のグループによって共有されている場合、各グループに対応する寄与重み係数が計算される。一実施形態において、寄与重み係数はノードを共有するグループの数によって決定される。他の実施形態においては、寄与重み係数は、グループ間の比（例えば面積比、質量比または体積比）によって決定される。
【０００６】
剛体壁と衝突する自動車の時間進行シミュレーションは、ＦＥＡモデルを用いて行われる。各ソリューションサイクルにおいて、ノードタイプに応じてそれぞれのノードについてノード力寄与（nodal force contribution）が計算される。その後、剛体壁のそれぞれのセグメントについて、ノード寄与重み係数で修正された、計算されたノード力寄与が、セグメントと接触する１以上のノードを有する各グループに対して合計される。総セグメント力は、セグメントに接触している全てのノードの、ノード寄与重み係数で修正されたノード力寄与を合計することにより決定される。剛体壁に作用する力の合計は、セグメント力の合計である。１のグループによって剛体壁に与えられる力の合計は、１のグループによってそれぞれのセグメントに与えられる力を合計することにより求められる。全体の詳細なＲＷＦ一覧（例えば、特定セグメントに対する１の特定グループのＲＷＦ時間履歴）が提供される。
【０００７】
本発明の一の目的は、自動車衝突の数値シミュレーションにおいて、車両のそれぞれの部品（つまり詳細なＲＷＦ一覧を得ることが望まれている有限要素のグループ）から剛体壁のそれぞれの空間的位置（つまりセグメント）に対して寄付する力の詳細な内訳を得ることにある。
【０００８】
本発明の他の、目的、特徴および利点は、添付した図面を参照し、以下の本発明の一実施形態の詳細な説明を考察することで、明らかとなろう。
【図面の簡単な説明】
【０００９】
本発明のこれらおよび他の、特徴、態様および利点は、以下の説明、添付の特許請求の範囲および添付した図面によって、より一層理解されるであろう。図面は次の通りである。
【図１Ａ】剛体壁と衝突する例示的な自動車を示す概略平面図である。
【図１Ｂ】剛体壁と衝突する例示的な自動車を示す概略側面図である。
【図１Ｃ】剛体柱と衝突する例示的な自動車を示す概略図である。
【図１Ｄ】剛体角部と衝突する例示的な自動車を示す概略図である。
【図２Ａ】本発明の一実施形態に係る、多数のロードセルが配置された例示的な剛体壁を示す図である。
【図２Ｂ】本発明の一実施形態に係る、図２Ａのロードセルの一つにそれぞれ対応する多数のセグメントを有する剛体壁定義を示す図である。
【図３】本発明の一実施形態に係る、多数のグループで構成される車両の例示的な有限要素解析モデルを示す概略図である。
【図４Ａ】本発明の一実施形態に係る、２以上のグループの有限要素によって共有される種々の例示的なノードを示す概略図である。
【図４Ｂ】本発明の一実施形態に係る、２以上のグループの有限要素によって共有される種々の例示的なノードを示す概略図である。
【図４Ｃ】本発明の一実施形態に係る、２以上のグループの有限要素によって共有される種々の例示的なノードを示す概略図である。
【図５】本発明の一実施形態に係る、例示的な詳細な剛体壁力一覧を示す図である。
【図６Ａ】図６Ａは、図６Ｂと共に、本発明の一実施形態に係る、剛体壁と衝突する車両の時間進行シミュレーションにおける詳細な剛体壁力一覧を提供する例示的なプロセスを示すフローチャートである。
【図６Ｂ】図６Ｂは、図６Ａと共に、本発明の一実施形態に係る、剛体壁と衝突する車両の時間進行シミュレーションにおける詳細な剛体壁力一覧を提供する例示的なプロセスを示すフローチャートである。
【図７】本発明の一実施形態を実現可能である、例示的なコンピュータシステムの主要な構成を示す機能図である。
【発明を実施するための形態】
【００１０】
以下の説明においては、本発明の完全な理解を提供するために、多数の具体的な詳細が記述される。しかしながら、本発明をこれらの具体的な詳細によらなくとも実現できることは、当業者にとって明らかである。ここでの説明および表現は、当該技術で経験のある者すなわち当業者が彼らの仕事を他の当業者に最も効果的に伝えるために用いる通常の手段である。他の事例で、よく知られている方法、手順、および構成は、本発明の態様が不必要に不明瞭となることを回避するために、詳細には説明しない。
【００１１】
ここで言及する「一実施形態」あるいは「ある実施形態」は、その実施形態に関して説明される特定の特徴、構造あるいは特性を、本発明の少なくとも１つの実施形態に含めることができるということを意味する。本明細書において記載される「一実施形態において」という語は、全てが同一の実施形態を言及しているものである必要はなく、また、他の実施形態とは相互に排他的な、別のあるいは代替的な実施形態である必要もない。さらに、本発明の１つ以上の実施形態を表わすプロセスフローチャートあるいはプロセス図におけるブロックの順序は、本質的に、いかなる特定の順序も示すものではなく、また、本発明におけるいかなる限定も意味するものでもない。
【００１２】
本発明の実施形態を、図１Ａから図７を参照して、ここに説明する。しかしながら、ここで行われるこれらの図面を参照した詳細な説明は、例示を目的とするものであって、当業者には、本発明がこれらの限られた実施形態よりも広いことを容易に認識するであろう。
【００１３】
剛体壁と衝突する車両の時間進行シミュレーションにおける詳細な剛体壁力（ＲＷＦ）一覧を提供する方法およびシステムを開示する。図１Ａ〜図１Ｂは、剛体壁１０４と衝突する例示的な自動車１０２の、概略的な、平面図および側面図である。このような衝突イベントの時間進行シミュレーションは、アプリケーションモジュール（例えば有限要素解析ソフトウェアパッケージ）がインストールされたコンピュータシステムにおいて行われる。剛体壁１０４を平面壁として示す。当業者であれば、他のタイプの剛体壁、例えば、図１Ｃに示す剛体柱１１４、図１Ｄに示す剛体角部１１６などを利用可能であることは理解されるであろう。図２Ａは、多数のロードセル２１２が配置された例示的な剛体壁２１０を示す。それぞれのセグメント２２２ａ〜２２２ｎを含む対応する剛体壁定義（rigid wall definition）２２０を図２Ｂに示す。また、剛体壁は静止している場合と動いている場合のいずれの場合もある。一般的に、正面衝突のシミュレーションにおいては、車両は、空間に固定された剛体壁へと運動する。しかしながら、屋根の押し潰しシミュレーションにおいては、運動する剛体壁が用いられる。
【００１４】
図６Ａを参照して、本発明の一実施形態にかかる、剛体壁と衝突する車両の時間進行シミュレーションにおける詳細な剛体壁力一覧を提供する例示的なプロセス６００を示すフローチャートを説明する。プロセス６００は、好ましくはソフトウェアで実行される。
【００１５】
プロセス６００は、ステップ６０２で、コンピュータシステム（例えば図７のコンピュータシステム７００）が、車両（例えば自動車）を表す有限要素解析（ＦＥＡ）モデルと剛体壁定義とを受けとることによりスタートする。ＦＥＡモデルは、時間進行シミュレーションを用いた衝突イベントにおける車両の構造挙動を取得するために用いられる。コンピュータシステムには、有限要素解析アプリケーションモジュールがインストールされている。ＦＥＡモデルは、グループを編成する複数の有限要素によって接続される複数のノードを含む。それぞれのグループは、ノードのリスト、有限要素のリスト、構造部品または部品の一部を含むことができる。一実施形態において、それぞれのグループは車両の部品を表す。例えば、図３に示すように、例示的な車両のＦＥＡモデル３００は、有限要素のグループ（明示的には描画せず）によって編成される。このような分割スキームにより、詳細なＲＷＦ一覧に、各グループからの寄与を含めることができる。さらに、図３に、バンパー３１１、フェンダ３１２、車体３１３、ホイール３１４、およびエンジンフード３１５を示されている。同じことを達成するために、他のスキームを利用することも可能である。例えば、バンパー３１１を、２つ以上のグループ（つまり構造部品の一部）へと分割することもできる。
【００１６】
次に判断ステップ６０４において、詳細な剛体壁力一覧を得ることが所望されるグループのリストがユーザによって定義されているか否かが判断される。定義されていない場合、ステップ６０６において、グループのリストが所定のスキームで生成される（例えば、ＦＥＡにおける各部品を１のグループに含める）。定義されている場合、プロセス６００は、「ｙｅｓ」に従ってステップ６０８へと分岐し、ＦＥＡモデルの全てのノードに対して、寄与重み係数が割り当てられる、又は、計算される。１つのグループのみに属するノードには、単一寄与重み係数（１．０）が割り当てられる。リスト内の２つ以上のグループに共有されるノードに対しては、所定のスキームを用いて、例えば、ノードを共有するグループの数に基づくスキーム、ノードを共有するグループ間の比に基づくスキーム等により、対応する重み係数が計算される。
【００１７】
複数のグループによって共有されるノードでは、各グループがそのノードに割り当てた重み係数の合計が１となる。図４Ａは、２つのノード４１２ａ〜４１２ｂを共有している２つのグループ４１０ａ〜４１０ｂを示している。図４Ｂは、いくつかのノード４２２ａ〜４２２ｄを共有している３つのグループ４２０ａ〜４２０ｃを示している。図４Ｃは、ノード４３２ａ〜４３２ｅを共有している４つのグループ４３０ａ〜４３０ｄを示している。基準としてノードを共有するグループの数を用いる場合、隣接するグループ４１０ａ〜４１０ｂへの寄与がそれぞれ５０％であるので、ノード４１２ａ〜４１２ｂの寄与重み係数は５０％である。図４Ｂに示す例では、それぞれの隣接グループへの寄与重み係数は、ノード４２２ａ〜４２２ｃでは５０％であり、ノード４２２ｄでは３３％である。図４Ｃに示す例では、それぞれの隣接グループへの寄与重み係数は、ノード４３２ａ〜４３２ｄでは５０％であり、ノード４３２ｅでは２５％である。
【００１８】
面積比または体積比が計算の基礎として用いられる場合、ノード４２２ｄの寄与重み係数は、グループ４２０ａおよび４２０ｂについては２５％とでき、グループ４２０ｃについて５０％とできる。同様に、ノード４２２ｂおよび４２２ｃの寄与重み係数は、グループ４２０ａおよび４２０ｂについてはそれぞれ３３％とでき、グループ４２０ｃについては６７％とできる。
【００１９】
次に、ステップ６１０において、剛体壁と衝突する車両の時間進行シミュレーションが、ＦＥＡモデルおよび剛体壁定義を用いて行われる。時間進行シミュレーションは、時間に関する複数のソリューションサイクルを含んでいる。各ソリューションサイクルでは、ステップ６１２において、ＦＥＡモデルにおける各ノードのノード力寄与が、ノードタイプに応じて計算される。その後、計算されたノード力寄与は、対応するグループおよびセグメントに関する、詳細なＲＷＦ一覧へと集計される前に、対応する重み係数で修正される。一実施形態では、図６Ｂに示すように、ステップ６１２ａにおいて、剛体壁のそれぞれのセグメントについて、ノード寄与重み係数で修正された、計算されたノード力寄与が、セグメントと接触する１以上のノードを含む各グループに対して合計される。ステップ６１２ｂにおいて、総セグメント力は、セグメントに接触している全てのノードのノード寄与重み係数で修正されたノード力寄与を合計することにより決定される。ステップ６１２ｃにおいて、剛体壁に作用する力の合計は、セグメント力の合計である。ステップ６１２ｄにおいて、あるグループによって剛体壁に与えられる力の合計は、グループによって各セグメントに与えられる力を合計することで求められる。
【００２０】
ノード力寄与ベクトルは、時間進行シミュレーションにおいて摩擦オプションが有効である（つまり、摩擦がシミュレーションに含まれている）場合には、剛体壁に垂直な力と、摩擦成分を表す２つの接線力と、を有する。以下の式が、ノード力寄与の計算に用いられる。
【００２１】
適切なノード質量でモデル化された通常のノードに対して、
ＦNODE＝ＭNODE（Ａ＋−Ａ−）
である。
ここでは：
ＦNODEはノード力寄与ベクトルである。
ＭNODEはスカラーのノード質量である。
Ａ＋は剛体壁との接触後のノード加速度ベクトルである。
Ａ−は剛体壁との接触前のノード加速度ベクトルである。
【００２２】
ペナルティ法(penalty formulation)において正確なノードの質量なしでモデル化されているノード(例えば、剛体におけるノード)に対して、
ＦNODE＝Ｆｎ＋Ｆｔ
Ｆｎ＝（ＫΔｕ）ｎ
である。
ここでは：
FNODEはノード力寄与ベクトルである。
Fnは剛体壁に垂直なノード力寄与である。
Ｋはペナルティ法において割り当てられた剛性（つまり、ペナルティ）である。
Δｕは、剛体壁を通るノード(つまり、正確なノードの質量なしでモデル化されたノード)の貫通(penetration)である。
ｎは剛体壁に対する単位法線ベクトルである。
Ｆｔは法線力に直交する摩擦力である（シミュレーションにおいて摩擦オプションがシミュレーションにおいて有効である場合に含められる）。
【００２３】
一実施形態において、剛体におけるノードに対して割り当てられる剛性は、剛体の材料特性によって決定される。
【００２４】
最後に、ステップ６１４において、時間進行シミュレーションの終了後、詳細なＲＷＦ一覧が提供される。詳細なＲＷＦ一覧を種々の形式で提供できる。１の例として、図５に時間に対する剛体壁力を示す。異なる３つのグループを、スタックカーブ（stacked curves）として示す。総剛体壁力は曲線５１０として示され、曲線５１２および５１４は剛体壁の２つの異なるセグメント（つまり２つの別の空間的位置）における力寄与を示す。他の実施形態において、特定のセグメントに対する有限要素の特定グループのＲＷＦ時間履歴を用いることもできる（図示せず）。さらに他の実施形態において、全てのセグメントと相互作用する単一のグループのＲＷＦ時間履歴を用いることもできる（図示せず）。そして、プロセス６００は、終了する。
【００２５】
一の態様において、本発明は、ここに説明した機能を実行可能な１以上のコンピュータシステムに対してなされたものである。コンピュータシステム７００の一例を、図７に示す。コンピュータシステム７００は、プロセッサ７０４のようなプロセッサを１つ以上有する。プロセッサ７０４は、コンピュータシステム内部通信バス７０２に接続されている。種々のソフトウェアの実施形態を、この例示的なコンピュータシステムに関して説明する。この説明を見れば、関連する技術分野に習熟している者には、いかにして他のコンピュータシステムおよび／またはコンピューターアーキテクチャーを用いて本発明を実行するかが明らかである。
【００２６】
コンピュータシステム７００は、また、メインメモリ７０８好ましくはランダムアクセスメモリ（ＲＡＭ））を有しており、そして補助メモリ７１０を有することもできる。補助メモリ７１０は、例えば、１つ以上のハードディスクドライブ７１２、および／または、フレキシブルディスクドライブ、磁気テープドライブ、光ディスクドライブなどを表わす１つ以上のリムーバブルストレージドライブ７１４を含むことができる。リムーバブルストレージドライブ７１４は、よく知られている方法で、リムーバブルストレージユニット７１８を読み出し、および／または、リムーバブルストレージユニット７１８に書き込む。リムーバブルストレージユニット７１８は、リムーバブルストレージドライブ７１４によって読み出し・書き込みされるフレキシブルディスク、磁気テープ、光ディスクなどを表わす。上記の内容からわかるように、リムーバブルストレージユニット７１８は、コンピューターソフトウェアおよび／またはデータを内部に記憶している、コンピュータ可読媒体を有している。
【００２７】
代替的な実施形態では、補助メモリ７１０は、コンピュータプログラムあるいは他の命令をコンピュータシステム７００にロードすることを可能にする他の同様な手段を有することもできる。そのような手段は、例えば、リムーバブルストレージユニット７２２とインタフェース７２０とを有してもよい。そのようなものの例には、プログラムカートリッジおよびカートリッジのインタフェース（ビデオゲーム機に見られるようなものなど）や、リムーバブルメモリチップ（消去可能なプログラマブルＲＯＭ（ＥＰＲＯＭ（Erasable Programmable Read-Only Memory））、ＵＳＢ（Universal Serial Bus）フラッシュメモリ、あるいはＰＲＯＭ（Programmable Read-Only Memory）など）および関連するソケットや、リムーバブルストレージユニット７２２からコンピュータシステム７００にソフトウェアおよびデータを転送することを可能にする他のリムーバブルストレージユニット７２２およびインタフェース７２０が含まれうる。一般に、コンピュータシステム７００は、プロセススケジューリング、メモリ管理、ネットワーキングおよびＩ／Ｏ（Input/Output）サービスなどのタスクを行なうオペレーティングシステム（ＯＳ）ソフトウェアによって、制御され連係される。
【００２８】
通信インタフェース７２４も、バス７０２に接続することができる。通信インタフェース７２４は、コンピュータシステム７００と外部装置との間で、ソフトウェアおよびデータを転送することを可能にする。通信インタフェース７２４の例には、モデム、ネットワークインタフェイス（イーサネット（登録商標）・カードなど）、コミュニケーションポート、ＰＣＭＣＩＡ(Personal Computer Memory Card International Association)スロットおよびカードなど、が含まれうる。
【００２９】
コンピュータシステム７００は、特定のセットの規則（つまりプロトコル）に基づいて、データネットワーク上の他の演算装置と通信する。一般的なプロトコルのうちの１つは、インターネットにおいて一般に用いられているＴＣＰ／ＩＰ（伝送コントロール・プロトコル／インターネット・プロトコル）である。一般に、通信インタフェース７２４は、データファイルをアセンブルしてデータネットワーク上で伝達される小さいパケットに変換することを管理し、あるいは受信したパケットを逆アセンブルして元のデータファイルに変換することを管理する。さらに、通信インタフェース７２４は、パケットが正しい宛先に届くよう、それぞれのパケットのアドレス部を処理し、あるいはコンピュータシステム７００が宛先となっているパケットをインターセプトする。
【００３０】
本明細書において、「コンピュータが記録可能な記憶媒体」、「コンピュータが記録可能な媒体」および「コンピュータ可読媒体」という語は、通常、リムーバブルストレージドライブ７１４、および／または、ハードディスクドライブ７１２に組み込まれたハードディスクなどの媒体の意味で用いられる。これらのコンピュータプログラム製品は、コンピュータシステム７００にソフトウェアを提供する手段である。本発明は、このようなコンピュータプログラム製品に対してなされたものである。
【００３１】
コンピュータシステム７００は、また、コンピュータシステム７００をモニタ、キーボード、マウス、プリンタ、スキャナ、プロッタなどにアクセスさせる入出力（Ｉ／Ｏ）インタフェース７３０を有してもよい。
【００３２】
コンピュータプログラム（コンピュータ制御ロジックともいう）は、メインメモリ７０８および／または補助メモリ７１０にアプリケーションモジュール７０６として記憶される。コンピュータプログラムは、通信インタフェース７２４を介して受け取られてもよい。このようなコンピュータプログラムが実行された時、コンピュータプログラムによって、コンピュータシステム７００がここに説明される本発明の特徴を実行することが可能になる。詳細には、コンピュータプログラムが実行された時、コンピュータプログラムによって、プロセッサ７０４が本発明の特徴を実行することが可能になる。したがって、このようなコンピュータプログラムは、コンピュータシステム７００のコントローラを表わしている。
【００３３】
ソフトウェアを用いて発明が実行される実施形態においては、ソフトウェアをコンピュータプログラム製品に記憶でき、リムーバブルストレージドライブ７１４、ハードディスクドライブ７１２あるいは通信インタフェース７２４を用いてコンピュータシステム７００へとロードすることができる。アプリケーションモジュール７０６は、プロセッサ７０４によって実行された時、プロセッサ７０４に、ここに説明される本発明の機能を実行させる。
【００３４】
所望のタスクを達成するために、Ｉ／Ｏインタフェース７３０を介したユーザ入力によって、あるいは、ユーザ入力によることなしに、１つ以上のプロセッサ７０４によって実行することができる１つ以上のアプリケーションモジュール７０６を、メインメモリ７０８に、ロードすることもできる。動作中、少なくとも１つのプロセッサ７０４がアプリケーションモジュール７０６のうちの１つを実行すると、結果が演算されて補助メモリ７１０（つまりハードディスクドライブ７１２）に記憶される。時間進行シミュレーションの状況（例えば衝突の結果など）が、テキスト表現あるいはグラフィック表現で、Ｉ／Ｏインタフェース７３０を介してユーザに報告される。
【００３５】
本発明を、具体的な実施形態を参照しながら説明したが、これらの実施形態は単なる例示であって、本発明を限定するものではない。具体的に開示された例示的な実施形態に対する種々の変更あるいは変形を、当業者は思いつくであろう。例えば、グループの全てをデフォルトリストに含める所定の方法を説明し例示したが、他の同等な方法、例えば関心のある特定のグループを含める方法を用いて同じことを達成することもできる。さらに、１の例示的なグループ化方法のみを図３において示したが、ＦＥＡモデルを無数の方法で分割できる。さらに、いくつかの例示的な寄与重み係数を説明し例示した。しかし、同じことを達成するために、他の方法を用いることもできる。最後に、詳細なＲＷＦ一覧の一例のみを図５において示した。しかし、詳細なＲＷＦ一覧を提示する他の手段を、テキスト形式またはグラフィック形式で用いることもできる。つまり、本発明の範囲は、ここで開示した具体的で例示的な実施形態に限定されず、当業者が容易に想到するあらゆる変更が、本願の精神および範囲、そして添付の特許請求の範囲の権利範囲に含まれる。
【符号の説明】
【００３６】
１０２  自動車
１０４  剛体壁
１１４  剛体柱
１１６  剛体角部
２１０  剛体壁
２１２  ロードセル
２２０  剛体壁定義
２２２ａ〜２２２ｎ  セグメント
３００  車両のＦＥＡモデル
３１１  バンパー
３１２  フェンダ
３１３  車体
３１４  ホイール
３１５  エンジンフード
４１０ａ〜４１０ｂ  グループ
４１２ａ〜４１２ｂ  ノード
４２０ａ〜４２０ｃ  グループ
４２２ａ〜４２２ｄ  ノード
４３０ａ〜４３０ｄ  グループ
４３２ａ〜４３２ｅ  ノード
７００  コンピュータシステム
７０２  バス
７０４  プロセッサ
７０６  モジュール
７０８  メインメモリ（ＲＡＭ）
７１０  補助メモリ
７１２  ハードディスクドライブ
７１４  リムーバブルストレージドライブ
７１８  リムーバブルストレージユニット
７２０  インタフェース
７２２  リムーバブルストレージユニット
７２４  通信インタフェース
７３０  Ｉ／Ｏインタフェース
【図１Ａ】
0006267573.tif 000002
【図１Ｂ】
0006267573.tif 000003
【図１Ｃ】
0006267573.tif 000004
【図１Ｄ】
0006267573.tif 000005
【図２Ａ】
0006267573.tif 000006
【図２Ｂ】
0006267573.tif 000007
【図３】
0006267573.tif 000008
【図４Ａ】
0006267573.tif 000009
【図４Ｂ】
0006267573.tif 000010
【図４Ｃ】
0006267573.tif 000011
【図５】
0006267573.tif 000012
【図６Ａ】
0006267573.tif 000013
【図６Ｂ】
0006267573.tif 000014
【図７】
0006267573.tif 000015