ハイブリッドエレメント110、120、130、140、150、160を用いて有限要素法FEMにおけるソリッドエレメントと粒子平滑化法SPHにおける粒子との間のカップリング効果を与えるコンピュータシステム600において実行される方法500であって、構造を表すグリッドモデル220を定義するステップ502であって、該グリッドモデルは、複数のFEMソリッドエレメントと、少なくとも1つのSPH粒子と、前記FEMソリッドエレメントと前記少なくとも1つのSPH粒子との間の少なくとも1つのレイヤのハイブリッドエレメント110、120、130、140、150、160と、を有しており、該ハイブリッドエレメント110、120、130、140、150、160のそれぞれがソリッドレイヤと粒子レイヤとを有しているステップと、前記グリッドモデル220を用いて荷重条件における前記構造の時間進行シミュレーションを実行するステップ504、506、508、510、512、516、518であって、該シミュレーションが前記ソリッドレイヤにおいて計算される1セットのノード量と前記粒子レイヤにおいて計算される対応する内力とを交換して532、533、534、535、536カップリング効果510を得るステップと、各ソリューションサイクルにおいて時間進行シミュレーションの結果を取得するステップであって506、536、該取得された結果が記憶装置にファイルとして保存されるあるいはユーザの指示に応じてモニタに図表によって表示されるステップと、を備える方法。請求項1に記載の方法であって、前記ソリッドレイヤにおいて計算される前記1セットのノード量と前記粒子レイヤにおいて計算される前記対応する内力との交換532、533、534、535、536が、FEMに基づいて前記ソリッドレイヤにおける前記1セットのノード量と要素応力状態を計算するステップ532と、前記1セットのノード量と要素応力状態を前記ソリッドレイヤから前記粒子レイヤにマップするステップ533と、SPHに基づいて前記粒子レイヤにおける前記対応する内力を計算するステップ534と、前記対応する内力を前記ソリッドレイヤに伝達して、FEMに基づいて時間進行シミュレーションの次のソリューションサイクルでのノードの力を計算するステップ535と、によって達成される方法。請求項2に記載の方法であって、前記1セットのノード量533は、カップリング効果506、510の開始時におけるノード変位と、ノード速度と、ノード加速とを含んでいる方法。請求項3に記載の方法であって、前記ノード量533は、カップリング効果510の開始後のソリューションサイクルにおけるノード変位を含んでいる方法。請求項4に記載の方法であって、前記ノード変位533は、制約条件として前記対応するSPH粒子にマップされる方法。請求項4に記載の方法であって、前記対応する内力を計算する前記ステップは、前記交換532、533、534、535、536の際には中止される方法。請求項2に記載の方法であって、前記粒子レイヤには1つ以上のSPH粒子が含まれる方法。請求項2に記載の方法であって、前記要素応力状態には、前記それぞれのハイブリッドエレメントの材料状態400を追跡する応力状態履歴変数が含まれる方法。請求項8に記載の方法であって、前記それぞれのハイブリッドエレメントの材料状態400には、ポスト降伏424と、弾性領域406と、塑性領域408と、が含まれる方法。ハイブリッドエレメント110、120、130、140、150、160を用いて有限要素法FEMにおけるソリッドエレメントと粒子平滑化法SPHにおける粒子との間のカップリング効果を与えるコンピュータシステム600であって、1つ以上のアプリケーションモジュール606に関するコンピュータ可読コードを記憶しているメインメモリ608と、前記メモリ608に連結される少なくとも1つのプロセッサ604と、を備えており、前記少なくとも1つのプロセッサ604が前記メモリ608内の前記コンピュータ可読コードを実行して、これにより、前記1つ以上のアプリケーションモジュール606に、構造を表すグリッドモデル220を定義するオペレーション502であって、該グリッドモデルは、複数のFEMソリッドエレメントと、少なくとも1つのSPH粒子と、前記FEMソリッドエレメントと前記少なくとも1つのSPH粒子との間の少なくとも1つのレイヤのハイブリッドエレメント110、120、130、140、150、160と、を有しており、該ハイブリッドエレメント110、120、130、140、150、160のそれぞれがソリッドレイヤと粒子レイヤとを有しているオペレーションと、前記グリッドモデル220を用いて荷重条件における前記構造の時間進行シミュレーションを実行するオペレーション504、506、508、510、512、516、518であって、該シミュレーションがソリッドレイヤにおいて計算される1セットのノード量と粒子レイヤにおいて計算される対応する内力とを交換して532、533、534、535、536カップリング効果510を与えるオペレーションと、各ソリューションサイクルにおいて時間進行シミュレーションの結果を取得するオペレーションであって506、536、ユーザの指示に応じて該取得された結果が記憶装置にファイルとして保存されるあるいはモニタに図表によって表示されるオペレーションと、を実行させるシステム。請求項10に記載のシステムであって、前記ソリッドレイヤにおいて計算される前記1セットのノード量と前記粒子レイヤにおいて計算される前記対応する内力との交換532、533、534、535、536が、FEMに基づいて前記ソリッドレイヤにおけるノード量を計算532し、前記ノード量と前記要素応力状態を前記ソリッドレイヤから前記粒子レイヤにマップする計算し533、SPHに基づいて前記粒子レイヤにおける内力を計算し534、内力を前記ソリッドレイヤに伝達して、FEMに基づいて時間進行シミュレーションの次のソリューションサイクルでのノードの力を計算する535ことによって実行されるシステム。請求項1あるいは2に記載の方法500によって、ハイブリッドエレメント110、120、130、140、150、160を用いて有限要素法FEMにおけるソリッドエレメントと粒子平滑化法SPHにおける粒子との間のカップリング効果を有するコンピュータシステム600を制御する命令を有するコンピュータ可読媒体。最も興味深いFEAのタスクの1つは、非常に大きい変形例えば車の破砕、爆発シミュレーションを受ける構造を含んでいる衝突イベントをシミュレートすることである。最近のコンピュータ性能の向上につれ、技術者は、構造的な破壊あるいは破損を伴う衝突イベントにおける挙動を単にシミュレートすることを望むだけでなく、衝突イベントから全破壊前であって物理的降伏後の構造的挙動をシミュレートすることをも求めている。しかしながら、ソリッドエレメントを用いるFEAによってそのような現象をシミュレートすることは難しい。例えば、ソリッドエレメントが、発泡フォーム材料のバンパーが押しつぶされあるいは圧縮されて、過度に変形されあるいはぐしゃぐしゃにつぶれてしまうと、ゼロあるいはマイナスの体積となり、シミュレーションにおいて数値的な問題が起こることがある例えば、デジタルコンピューターにおいて無効な数のためにシミュレーションが中止される。ゼロあるいはマイナスの体積になる問題を解くために、破壊されたソリッドエレメントは粒子平滑化法smoothedparticlehydrodynamicsSPHにより粒子で置き換えられる。しかしながら、FEMとSPHの数学的定式化マセマティカルフォーミュレーションは異なる。粒子とソリッドエレメントとを同じモデルに共存させるためには、粒子とソリッドエレメントとを接続するある種類の接続を確立しなければならない。従来技術アプローチは、ある粒子をソリッドエレメントに剛体的にリジッドに接続するタイド連結インタフェースtiedinterfaceを用いている。しかしながら、このアプローチでは、多くは、タイドインタフェースつまり剛体的リンクの配置が不定であるために、さらに非常に非現実的なシミュレーション結果となる。例えば、粒子とソリッドエレメントを一体的に連結することは、初期には妥当なこともある。しかし、粒子およびソリッドエレメントの変形が予測不能となるにつれ、粒子およびソリッドエレメントの剛体的リンクの配置が不定であるため、非常に非現実的な接続となってしまう。したがって、SPH粒子およびFEMソリッドが従来技術アプローチの問題および欠点を回避するよう共存できるコンピュータ支援工学解析モデルにおけるより現実的なインタフェースが望まれよう。SPH粒子とFEMソリッドとの間のカップリング効果を有するハイブリッドエレメントを開示する。本発明の一の面では、ハイブリッドエレメントは、有限要素法FEMに基づいたソリッドエレメントと粒子平滑化法SPHに基づいた1つ以上の対応する粒子とのカップリング効果を有するように構成される。ハイブリッドエレメントは、SPH粒子とFEMソリッドとの間の緩衝バッファーあるいはインタフェースとしてコンピュータ支援工学CAEグリッドモデルにおいて定義される。例えば、グリッドモデルの一部はSPH粒子からなる。大きい変形が持続する可能性のためである。一方、残りのモデルはFEMソリッドエレメントからなる。ハイブリッドエレメントはソリッドと粒子との間に配置される。各ハイブリッドエレメントは2つのレイヤ層からなる。ソリッドレイヤと粒子レイヤとである。第一に、ハイブリッドエレメントのカップリング効果は、下記のように達成される:FEMに基づいてソリッドレイヤにおいてノード加速、速度および変位を要素応力とともに計算し、計算されたこれらノード量を要素応力状態とともに粒子レイヤへとマップし要素応力状態には、要素の応力値および現在の材料状態例えば、弾性、塑性、降伏、歪み硬化などが含まれる、SPHに基づいて粒子レイヤにおける内力を計算し、その内力をソリッドレイヤに転送して次のソリューションサイクルにおけるノードの力を計算する。以降のソリューションサイクルにおいては、ソリッドレイヤのノード変位が更新されて、対応するSPH粒子にマップされる。内力は、SPHに基づいて計算され、そして次のソリューションサイクルにおけるノードの力を計算するためにソリッドレイヤに転送される。本発明の他の面では、コンピュータ支援解析モデルは、固体有限要素を外周部に有することができる、一方、SPH粒子をその他の部分に有することができる。外周部あるいは周辺部上のソリッドエレメントは、境界条件を与えるよう構成される。さらなる他の面では、SPH粒子は、腐食して降伏限界を超えたソリッドエレメントと置き換わるよう用いられる。SPH粒子を、材料の歪み硬化効果を表す、より柔らかい材料モデルによってモデル化することができる。本発明の他の目的、特徴および利点は、添付した図面を参照し、以下の本発明の実施の形態の詳細な説明を考察することによって明らかとなろう。本発明のこれらおよび他の特徴、面および利点は、以下の説明、添付の特許請求の範囲および添付した図面を考慮してより理解されよう。図面は次の通りである。まず図1を参照して、本発明の一の実施形態にかかる種々の例示的なハイブリッドエレメントの図を示す。ハイブリッドエレメントは、2つの部分からなる。ソリッドレイヤと、対応する粒子レイヤと、である。ソリッドレイヤは、FEMに基づいたソリッドエレメントからなり、対応する粒子レイヤはSPHに基づいた1つ以上の粒子からなる。ソリッドエレメントは、限定するものではないが、六面体、くさび形や四面体からなる。ハイブリッドエレメント110は、一つの対応する粒子を有する六面体である。要素エレメント120は、1つの粒子を有するくさび形要素である。また、要素130は、1つの粒子を有する四面体である。2つ以上の粒子を有する例示的なハイブリッドエレメントは、8つの粒子を有する要素140であり、6つの粒子を有する要素150であり、4つの粒子を有する要素160である。本発明の他の実施形態では、他の数の粒子も、実現することができる。例えば、27個の粒子を有する六面体である図示せず。ハイブリッドエレメントのカップリング効果は、ソリッドレイヤを粒子レイヤと関連付けることによって達成される。内力計算プロシージャの詳細を、図4Cに示すフローチャートに示す。例えば、ソリッドレイヤは、対応する粒子レイヤにおいて粒子に対する制約条件として機能する。一の実施形態において、ソリッドレイヤの体積は、粒子を取り巻く領域として構成される。図2A~図2Bは、構造つまりグリッドモデルとして部分的に示すプレート220に比較的高速度矢印215で示すで当たる物体つまり剛性球体の態様の発射体210、したがって衝撃荷重、のシーケンスを示す。衝撃荷重を受ける構造の部分を、要素225点線で示すによってモデル化している。要素225は、初期にはFEMソリッドエレメントによってモデル化できる。球体210からの衝突がプレート220に衝撃を与えるとき、要素225は破壊されるあるいは降伏する場合がある材料破壊あるいは降伏の定義は図4および対応する説明を参照。壊れた要素はSPH粒子と置き換えられ、シミュレーションが継続される。SPH粒子およびFEMソリッドエレメントは異なる定式化を用いるので、ハイブリッドエレメントの態様のインタフェースがSPH粒子とFEMソリッドエレメントとの間に生成され、カップリング効果を有する。さらに上記の実施例を説明するために、図3A~図3Dは、構造プレート220の平面図のシーケンスを示している。開始では、図3Aにおいてプレート220はすべてソリッドエレメントとして示される。次に、図3Bにおいて、中心のソリッドエレメントが壊れ、SPH粒子中心にドット描き斜線を入れた円で示すと置き換えられる。これは、発射体/球体210によって起きる場合があり、プレート220との強い接触を意味する。少なくとも1つのレイヤのハイブリッドエレメント網掛けの要素で示すが、SPH粒子とソリッドエレメントとの間のカップリング効果のインタフェースとして生成される。その後、図3Cにおいて、中心の要素の周囲のより多くの要素が壊れ、SPH粒子と置き換えられる。見ての通り、ハイブリッドエレメントのインタフェースは、常にSPH粒子とソリッドエレメントとの間に位置するよう、動的に調整される。さらに、ハイブリッドエレメントをCAEモデルの境界に配置し、残りのモデルをSPH粒子とすることもできる。図3Dに示す構成は、本発明のこの面を例示している。図4は、本発明の一の実施形態にかかる、降伏後構造的挙動post-yieldingstructuralbehaviorsを決定するために用いることができる例示的な応力-歪曲線を示している。曲線400では、垂直軸が応力402を表し、水平軸が歪み404を表している。材料は、2つの領域を有する。弾性領域406と塑性領域408とである。塑性領域408は、3つのカテゴリーにさらに分類される。降伏yieldingカテゴリー424と、歪み硬化strainhardeningカテゴリー426と、くびれneckingカテゴリー428と、である。応力-歪み曲線400の弾性領域の頂端部は、降伏応力に対応する降伏点414である。破壊応力ultimatestressは、極限強さ点ultimatestrengthpoint416に対応し、また、破損すなわち破壊応力は破壊位置418に対応する。一の実施形態では、FEMソリッドエレメントは材料の弾性挙動をモデル化するよう用いられる。材料が降伏を越えるとすぐに、SPH粒子がソリッドエレメントと置き換わるよう生成される。置き換えられたSPH粒子を、歪み硬化効果をより現実的にシミュレートすることができるよう、より柔らかい材料モデルによってモデル化することができる。次に図5Aを参照して、衝撃荷重に応じて降伏後構造的挙動および大きい変形を数値的にシミュレートするようハイブリッドエレメントを用いる例示的なプロセス500を示す。プロセス500は、好ましくはソフトウェアで実行される。プロセス500は、ステップ502において、構造例えば車、飛行機のコンピュータ支援分析グリッドモデル例えばFEMグリッドモデルを定義することによってスタートする。グリッドモデルは、大きい変形をよく受けやすい構造例えば衝突シミュレーションにおける車のバンパーの部分を表す1つ以上のハイブリッドエレメントを有する。グリッドモデルは時間進行シミュレーションにおいて用いられる。次に、ステップ504においては、すべての要素およびパラメータは、時間進行シミュレーションにおいて開始つまり時間がゼロ、あるいは第1ソリューションサイクル時に初期化される。そして、プロセス500は、判断506においてハイブリッドエレメントのカップリング効果が有効にされたか否かをチェックする。noであれば、プロセス500は、FEMにおけるソリッドエレメントのようにハイブリッドエレメントを処理することによって、シミュレーションを行うステップ508に移行する。ステップ508の詳細を、図5Bおよび対応する説明において説明する。言いかえれば、カップリング効果が有効にされていない場合、時間進行シミュレーションがFEMを用いて行われる。そうでなく、yesであれば、プロセス500は、カップリング効果を含めるようハイブリッドエレメントによって時間進行シミュレーションを行うステップ510に移行する。図5Cおよび関連する説明はステップ510に関するものである。プロセス500は、次のソリューションサイクルのシミュレーション時間へとインクリメントするステップ516に移行する。そして、判断518において、時間進行シミュレーションが終了に達したか否かが判定される。例えば、シミュレーション時間を所定の総シミュレーション時間と照合する。終了に達してなければ、プロセス500は、判断518が真になるまで、判断506に戻って次のソリューションサイクルの残りのステップを繰り返す。そして、プロセス500は、終了する。図5Bは、ステップ508のさらなる詳細を示す。ステップ522において、プロセス500が、ハイブリッドエレメントを有する各要素のノード加速、速度および変位を取得する。一の実施形態において、ノード量は、FEMにおける明示的求解器explicitsolverにおいて取得される例えばf=m×a、ここでfはノードの力であり、mはノード質量であり、aはノード加速である。次に、ステップ524において、要素の内力がFEMにおいてソリッド定式化に応じて計算される。最後に、ステップ526において、次のソリューションサイクルのノードの力を、要素の内力からの寄与を含めるよう計算できる。どのノードも、すべての接続している要素からの寄与を受けることができる。図5Cは、ステップ510のさらなる詳細を示す。ステップ532においては、プロセスが、FEMにおけるステップ522とほぼ同様に、ハイブリッドエレメントのソリッドレイヤにおいてノード量つまりノード加速、速度および変位、そして要素応力状態を取得する。次に、ステップ533においては、取得されたノード量および要素状態が、接続の最初に、対応する粒子レイヤへとマップ対応関係指定される。そして、以降のサイクルにおいては、ノード変位のみが対応する粒子レイヤにマップされる。ソリッドレイヤの変位が更新されて、制約条件として、対応するSPH粒子にマップされる。言いかえれば、SPH粒子は、FEMに基づいて計算されるノード変位によって制限される。内力は、粒子レイヤにおいて計算され、FEMに基づいて次のソリューションサイクルにおけるノードの力を計算するためにソリッドレイヤに転送される。言いかえれば、内力はSPHを用いて取得され、ソリッドレイヤにおける内力計算はハイブリッドエレメントにおいて交換の際には中止される。要素応力状態には、その要素に対して計算された応力値が少なくとも含まれる。さらに、要素応力状態には、要素応力状態履歴変数に要素の現在状態が含まれる。この変数あるいは他の同等な手段が、時間進行シミュレーションにわたって要素の状態の追跡のために用いられる。言いかえれば、要素の降伏後状態を、その要素の履歴変数から判定できる。そして、ハイブリッドエレメントの内力が、ステップ534において、SPH定式化に基づいて対応する粒子レイヤにおける粒子において計算される。次に、ステップ535においては、内力がソリッドレイヤに転送される。言いかえれば、ソリッドレイヤの要素の内力は、対応する粒子レイヤにおける粒子から計算されたものと置き換えられる。最後に、ステップ536においては、次のソリューションサイクルにおけるノードの力が、ステップ526とほぼ同様に要素の内力の寄与を含めて計算される。一の面において、本発明は、ここに説明した機能を実行可能な1つ以上のコンピュータシステムに対してなされたものである。コンピュータシステム600の一例を、図6に示す。コンピュータシステム600は、プロセッサ604など1つ以上のプロセッサを有する。プロセッサ604は、コンピュータシステム内部通信バス602に接続されている。種々のソフトウェアの実施形態を、この例示的なコンピュータシステムの点から説明する。この説明を読むと、いかにして、他のコンピュータシステムおよび/またはコンピューターアーキテクチャーを用いて、本発明を実行するかが、関連する技術分野に習熟している者には明らかになるであろう。コンピュータシステム600は、また、メインメモリ608好ましくはランダムアクセスメモリRAMを有しており、そして二次メモリ610を有することもできる。二次メモリ610は、例えば、1つ以上のハードディスクドライブ612、および/またはフレキシブルディスクドライブ、磁気テープドライブ、光ディスクドライブなどを表わす1つ以上のリムーバブルストレージドライブ614を有することができる。リムーバブルストレージドライブ614は、よく知られている方法で、リムーバブルストレージユニット618を読み取りおよび/またはリムーバブルストレージユニット618に書き込む。リムーバブルストレージユニット618は、リムーバブルストレージドライブ614によって読み取り書き込みされるフレキシブルディスク、磁気テープ、光ディスクなどを表わす。以下にわかるように、リムーバブルストレージユニット618は、コンピューターソフトウェアおよび/またはデータを内部に記憶しているコンピュータで使用可能な記憶媒体を有している。代替的な実施形態において、二次メモリ610は、コンピュータプログラムあるいは他の命令をコンピュータシステム600にロードすることを可能にする他の同様な手段を有することもできる。そのような手段は、例えば、リムーバブルストレージユニット622とインタフェース620とを有することができる。そのようなものの例には、プログラムカートリッジおよびカートリッジのインタフェースビデオゲーム機に見られるようなものなどと、リムーバブルメモリチップ消去可能なプログラマブルROMEPROM、ユニバーサルシリアルバスUSBフラッシュメモリ、あるいはPROMなどおよび関連するソケットと、ソフトウェアおよびデータをリムーバブルストレージユニット622からコンピュータシステム600に転送することを可能にする他のリムーバブルストレージユニット622およびインタフェース620と、が含まれうる。一般に、コンピュータシステム600は、プロセススケジューリング、メモリ管理、ネットワーキングおよびI/Oサービスなどのタスクを行なうオペレーティングシステムOSソフトウェアによって、制御され連係される。通信用インタフェース624も、また、バス602に接続することができる。通信用インタフェース624は、ソフトウェアおよびデータをコンピュータシステム600と外部装置との間で転送することを可能にする。通信用インタフェース624の例には、モデム、ネットワークインターフェイスイーサネット登録商標カードなど、コミュニケーションポート、PCMCIAPersonalComputerMemoryCardInternationalAssociationスロットおよびカードなど、が含まれうる。ソフトウェアおよびデータは通信インタフェース624を介して転送される。コンピュータ600は、専用のセットの規則つまりプロトコルに基づいて、データネットワーク上の他の演算装置と通信する。一般的なプロトコルのうちの1つは、インターネットにおいて一般に用いられているTCP/IP伝送コントロールプロトコル/インターネットプロトコルである。一般に、通信インタフェース624は、データファイルをデータネットワーク上で伝達される小さいパケットへのアセンブリングを管理し、あるいは受信したパケット元のデータファイルへと再アセンブルする。さらに、通信インタフェース624は、正しい宛先に届くようそれぞれのパケットのアドレス部分に対処し、あるいはコンピュータ600が宛先となっているパケットを他に向かわせることなく受信する。この書類において、コンピュータプログラム媒体、コンピュータが読取り可能な媒体、コンピュータが記録可能な媒体およびコンピュータが使用可能な媒体という語は、リムーバブルストレージドライブ614例えば、フラッシュストレージドライブおよび/またはハードディスクドライブ612に組み込まれたハードディスクなどの媒体を概して意味して用いられている。これらのコンピュータプログラム製品は、コンピュータシステム600にソフトウェアを提供する手段である。本発明は、このようなコンピュータプログラム製品に対してなされたものである。コンピュータシステム600は、また、コンピュータシステム600をアクセスモニタ、キーボード、マウス、プリンタ、スキャナ、プロッタなどに提供する入出力I/Oインタフェース630を有することができる。コンピュータプログラムコンピュータ制御ロジックともいうは、メインメモリ608および/または二次メモリ610にアプリケーションモジュール606として記憶される。コンピュータプログラムを、通信用インタフェース624を介して受け取ることもできる。このようなコンピュータプログラムが実行された時、コンピュータプログラムによって、コンピュータシステム600がここに説明した本発明の特徴を実行することが可能になる。詳細には、コンピュータプログラムが実行された時、コンピュータプログラムによって、プロセッサ604が本発明の特徴を実行することが可能になる。したがって、このようなコンピュータプログラムは、コンピュータシステム600のコントローラを表わしている。ソフトウェアを用いて発明が実行される実施形態において、ソフトウェアをコンピュータプログラム製品に記憶でき、リムーバブルストレージドライブ614、ハードドライブ612あるいは通信用インタフェース624を用いてコンピュータシステム600へとロードすることができる。アプリケーションモジュール606は、プロセッサ604によって実行された時、アプリケーションモジュール606によって、プロセッサ604がここに説明した本発明の機能を実行する。所望のタスクを達成するために、I/Oインタフェース630を介したユーザ入力によってあるいはよることなしに、1つ以上のプロセッサ604によって実行することができる1つ以上のアプリケーションモジュール606例えばFEMおよび/またはSPHアプリケーションモジュールを、メインメモリ608に、ロードすることもできる。動作においては、少なくとも1つのプロセッサ604がアプリケーションモジュール606のうちの1つが実行されると、結果が演算されて二次メモリ610つまりハードディスクドライブ612に記憶される。解析の状況が、ユーザの指示に応じてテキストあるいはグラフィック表現でI/Oインタフェース630を介してユーザに報告される。本発明を具体的な実施形態を参照しながら説明したが、これらの実施形態は単なる例示であって、本発明を限定するものではない。開示した例示的な実施形態に対する種々の変更あるいは変形を、当業者は思いつくであろう。例えば、大きい変形を受ける例示的な構造を発射体が衝突するプレートとして例示し説明したが、衝撃荷重下にある他の構造例えば車の衝突時の自動車のバンパーを本発明の特許請求の範囲に記載した方法によって数値的にシミュレートすることができる。さらに、ソリッドエレメントを六面体、くさび形および四面体として例示し説明したが、代わりに他のタイプのソリッドエレメント例えば五面体を用いることもできる。つまり、発明の範囲は、ここで開示した具体的で例示的な実施形態に限定されず、当業者が容易に想到するあらゆる変更が、本願の精神および認識範囲そして添付の特許請求の範囲の権利範囲に含まれる。