"製品の多目的工学設計最適化において最も影響力のある設計変数を識別するためのコンピュータシステムにおいて実行される方法であって、コンピュータシステムにおけるアプリケーションモジュールによって、最適化される製品の記述を受け取るステップであって、該記述が第1セットの設計変数と工学設計最適化に対する目的および制約としての1セットの応答関数とを含んでいるステップと、前記アプリケーションモジュールによって、第1セットの設計変数によって定義された第1設計空間から選択される1セットの候補製品の前記1セットの応答関数の結果を取得するステップと、前記アプリケーションモジュールによって、前記取得された結果に基づいて前記第1セットの設計変数のそれぞれ対前記1セットの応答関数のグローバル感度解析から前記第1セットの設計変数のうちのどれが最も影響力のある設計変数であるかを決定し、よって第2セットの設計変数を選択するステップであって、前記グローバル感度解析は、それぞれの重み因子が与えられた、前記1セットの応答関数に対する前記第1セットの設計変数のそれぞれの影響を表すトータルグローバル感度指標を示すよう構成されたグラフィカルプレゼンテーションツールを含んでおり、前記第2セットの設計変数は、製品の工学設計最適化がその後行われる第2設計空間を定義するために用いられるステップと、を備える方法。請求項1に記載の方法であって、前記第1セットの設計変数のそれぞれは、上下限境界値によって定義される範囲を含んでいる方法。請求項1に記載の方法であって、前記第1設計空間は前記第2設計空間より大きい方法。請求項1の方法であって、さらに、前記アプリケーションモジュールによって、前記取得された応答を適合させることにより1セットのメタモデルを生成するステップを備える方法。請求項1に記載の方法であって、前記第2セットの設計変数のそれぞれは、所定のしきい値より高いトータル感度指標を有する方法。請求項1に記載の方法であって、第2セットの設計変数のそれぞれは第2セットの設計変数を選択する前記ステップが、前記第1セットの設計変数をグローバル感度指標の降順で並べ換えるステップと、第1セットの設計変数のうち最も高い感度指標を有する設計変数を、前記選択された設計変数の累積感度指標が予め定義されたしきい値より大きくなるまで、選択するステップと、を備えている方法。請求項1に記載の方法であって、前記重み因子は、グラフィカルプレゼンテーションツールにおける対話型調整機構を用いてユーザによって調整される方法。請求項1に記載の方法であって、前記グラフィカルプレゼンテーションツールは、累積トータル感度指標を示す積み重ね棒グラフである方法。請求項8に記載の方法であって、前記積み重ね棒グラフは、前記第1セットの設計変数の前記それぞれに対して1つのバーを有しており、該バーは応答関数のうちの1つにそれぞれに対応する複数の成分を有している方法。請求項9に記載の方法であって、前記複数の成分は前記それぞれの重み因子に応じて変化する方法。請求項1の方法であって、さらに、前記アプリケーションモジュールによって、前記応答関数の重み因子を対話式に調整する手段を図表式に表現するステップを備える方法。請求項11に記載の方法であって、重み因子を対話式に調整する前記手段は、ユーザが対話式に操作できるシミュレート用スライドメーターを備えている方法。請求項1に記載の方法であって、前記図表式表現は円グラフである方法。請求項1に記載のコンピュータシステムにおいて実行される方法によって、製品の多目的工学設計最適化において最も影響力のある設計変数を識別する命令を有するコンピュータ可読媒体。製品の多目的工学設計最適化において最も影響力のある設計変数を識別するシステムであって、少なくとも1つのアプリケーションモジュールに関するコンピュータ読取り可能なコードを記憶しているメインメモリと、前記メインメモリに連結される少なくとも1つのプロセッサであって、該少なくとも1つのプロセッサがメインメモリ内のコンピュータが読み取り可能なコードを実行して、前記少なくとも1つのアプリケーションモジュールに、方法に基づいてオペレーションを実行させる少なくとも1つのプロセッサと、を備えるシステムであって、該方法が、最適化される製品の記述を受け取るステップであって、該記述が第1セットの設計変数と工学設計最適化に対する目的および制約としての1セットの応答関数とを含んでいるステップと、第1セットの設計変数によって定義された第1設計空間から選択される1セットの候補製品の前記1セットの応答関数の結果を取得するステップと、前記取得された結果に基づいて前記第1セットの設計変数のそれぞれ対前記1セットの応答関数のグローバル感度解析から前記第1セットの設計変数のうちのどれが最も影響力のある設計変数であるかを決定することによって第2セットの設計変数を選択するステップであって、前記グローバル感度解析は、それぞれの重み因子によって変倍された、前記第1セットの設計変数のそれぞれの前記1セットの応答関数に対する影響を表すトータルグローバル感度指標を示すよう構成されたグラフィカルプレゼンテーションツールを含んでおり、前記第2セットの設計変数は、製品の工学設計最適化がその後行われる第2設計空間を定義するために用いられるステップと、を備えているシステム。非常に大きな数の連結部分を有する複雑な製品例えば自動車、航空機などの最適化設計を行なう場合、多数のパラメータあるいは設計変数が必要となり、ひいては設計空間のサイズが増加する。その結果、工業的な問題あるいは現実的な問題例えば車耐衝撃性シミュレーションの応答関数の評価を必要とする最適化のトータルコスト工学的観点および演算的観点からは、全設計空間に対する調査である設計が多数あるために、非常に高くなる。しかしながら、一般に、いくつかの主要なパラメータのみが応答関数に最も影響力のある。その多数のパラメータを最も影響力のあるいくつかに絞り込むことによって、設計空間を小さくすることができ、ひいては最適化を行なうためのコストを低減することができることになる。したがって、効率的かつ有効かつ簡単かつ直感的な方法で工学設計最適化において最も影響力のある設計変数を識別することが望まれる。製品の多目的工学設計最適化において最も影響力のある設計変数を識別する方法を開示する。本発明の一の面では、製品が、目的および制約として1セットの設計変数および1セットの応答関数を用いて最適化される。代表的な製品設計選択肢サンプルが、設計空間設計変数によって定義されるから選択され、応答に関して評価される。そして、メタモデルが、サンプルの応答を設定し、すべての設計変数対応答関数のグローバル感度解析を容易にするよう、用いられる。グラフィカルプレゼンテーションツールは、それぞれの重み因子を応答関数に対話式に与えるスケールするscaleことによってユーザが起こりうる事態what-ifのシナリオを実行し、最も影響力のある設計変数を容易に識別できるよう、構成される。最も影響力のある設計変数が識別されると、その場合、工学設計最適化を低減された設計空間において行うことができる。低減された設計空間は、最も影響力のある設計変数の全範囲によって、それほど重要でない残りの設計変数が特定の値に固定された状態で、定義される。その特定の値は、最適解から得られた値あるいは設計変数の基準値とすることができる。一の態様では、本発明は、製品の多目的工学設計最適化において最も影響力のある設計変数を識別するコンピュータシステムにおいて実行される方法に対してなされたものであり、最適化される製品の記述を受け取るステップと、1セットの候補製品の1セットの応答関数の結果を取得するステップと、第2セットの設計変数を選択するステップと、を有する。最適化される製品の記述を受け取るステップでは、該記述が第1セットの設計変数と工学設計最適化に対する目的および制約としての1セットの応答関数とを含んでいる。1セットの候補製品の1セットの応答関数の結果を取得するステップでは、1セットの候補製品の1セットの応答関数は、第1セットの設計変数によって定義された第1設計空間から選択される。第2セットの設計変数を選択するステップでは、取得された結果に基づいて第1セットの設計変数のそれぞれ対1セットの応答関数のグローバル感度解析から第1セットの設計変数のうちのどれが最も影響力のある設計変数であるかを決定することによって第2セットの設計変数を選択する。また、グローバル感度解析は、それぞれの第1セットの設計変数のそれぞれの1セットの応答関数に対する影響を表す重み因子によって変倍されたトータルグローバル感度指標を示すよう構成されたグラフィカルプレゼンテーションツールを含んでいる。また、第2セットの設計変数は、製品の工学設計最適化がその後行われる第2設計空間を定義するために用いられる。本発明の他の目的、特徴および利点は、添付した図面を参照し、以下の本発明の実施の形態の詳細な説明を考察することによって明らかとなろう。本発明のこれらおよび他の特徴、面および利点は、以下の説明、添付の特許請求の範囲および添付した図面を考慮してより理解されよう。図面は次の通りである。まず図1を参照して、管状の構造部材102つまり例示的な工学製品が、ある設計荷重条件の下で、設計目的が与えられた材料例えば通常の強さ鋼の重さの最小化、したがって、コストの最小化である工学最適化において、最適化される。厚さ104が薄いほど、構造物の重量が小さくなることは明らかである。しかしながら、ある点では、その構造が設計荷重に耐えるには弱くなりすぎることになる例えば、材料降伏および/または座屈による構造的な破壊。したがって、この管状の構造体の工学最適化は、強さを最大限にし、その結果より安全な構造体とする他の設計目的を必要とする。この例示的なケースにおいて、厚さ104は、設計空間としての上下限境界間の範囲例えば8分の1インチから1/2インチを有することができる設計変数である。あらゆる設計選択肢がこの空間から選択される。多目的進化型アルゴリズムにおいて、それぞれの世代における母集団すなわち設計選択肢は設計空間から選択される。1つの設計変数だけである場合、設計空間は一次元例えば線である。2つの変数に対しては、設計空間は二次元領域になる、などである。3つを超える設計変数に対しては、設計空間は図示することができない超空間ハイパースペースhyperspaceである。図2に示すフローチャートは、1セットの最も影響力のある設計変数を識別することによって、製品例えば車、航空機などの工学設計最適化を行なう例示的なプロセス200を示す。プロセス200は、好ましくはソフトウェアで実行される。プロセス200は、ステップ202において、最適化される製品の記述を受け取ることによってスタートする。工学設計最適化は、設計目的および制約として1セットの設計変数つまり可能な限り多くの設計変数および1セットの応答関数を用いて実行される。耐衝撃性に関する自動車設計の最適化において、設計目的および制約には、例えば、車の質量を最小限にすること、予め決められた衝撃を受ける防護壁の陥入を最小限にすることが含まれる。それぞれの設計変数は、初期値と、上下限境界値の範囲と、を有する。工学設計最適化の目標は、すべての設計変数の範囲によって境界付けられた設計空間つまりn-次元の四次元以上の空間内にある最適解あるいは1セットのパレート最適解を得ることである。非常に多くの設計変数がある場合、工学設計最適化を実行するのに必要なコストと時間は相当なものになる。例えば、自動車耐衝撃性など現実世界の応用において、何百もの設計変数が想定される。次に、ステップ204において、最も影響力のある設計変数が、対話型起こりうる事態判断ツールつまりユーザがあるパラメータを対話により調整することを可能にする図表式表現アプリケーションと共に、グローバル感度解析を用いて、1セット全体の設計変数をスクリーニングすることによって選択される。他の態様において、ユーザは、また、最も影響力のある設計変数を判断するために、重み因子を予め定義することができる。グローバル感度解析の理論的定式化は次のように要約される。領域[0,1]における独立した一様に分布するランダム入力変数xのベクトルの関数としての自乗積分応答squareintegrableresponseの関数fxが想定される。関数は、次元数を増加させる関数の和として分解することができる。ここでである。もし以下の条件がに対して課せられる場合、方程式1に記載した分解は一義的である。グローバル感度解析のコンテキストにおいて、Vfによって表される関数fの全分散totalvarianceは、次の方程式と同等なものとしてできる。ここで、であり、方程式3におけるそれぞれの項は、全分散に対するそれぞれの独立変数Viあるいは1セットの変数例えばVijの部分的寄与partialcontributionあるいは偏分散partialvarianceを表し、それらの相対的な重要性を表示している。偏分散は以下の表現を用いて計算することができる。等がそれであり、VおよびEは分散および期待値をそれぞれ示している。なお、である。この定式化により、独立変数および1セットの変数に対応する感度指標の計算が容易になる。例えば、第1順序感度指標および第2順序感度指標は次のように計算することができる。独立モデル入力仮定の下では、すべての感度指標の和は1に等しい。所定の変数に対する第1順序感度指標Siは、変数の主な効果を表すが、変数の相互作用の効果を考慮していない。全分散に対する変数のトータルな寄与totalcontributionは、すべての相互作用の和として定められ、変数の主な効果、対応して、トータル感度指標Sitotalが次のように定義される。上記の参考表現は、あらゆる応答関数例えば放射基底関数radialbasisfunctionのメタモデルmetamodelを用いて容易に評価できる。分散に基づいたノンパラメトリックアプローチvariancebasednon-parametricapproachにおいて、グローバル感度は、モンテカルロ法を用いて、設計変数のあらゆる組み合わせに対して推定される。あらゆる設計変数xiのトータル感度を計算するために、設計変数セットは、xiおよびの2つ相補的な部分セットsubsetに分割される。これらの部分セットを用いる目的は、xiの影響を、部分セットZに含まれる残りの設計変数の影響から遮断することにある。そして、xiに対するトータル感度指標は次のように定義される。ここで、Viはxiに関する応答の偏分散であり、ViZは、xiとZとの間の相互作用に依存する応答分散の程度themeasureoftheresponsevarianceである。同様に、Zに対する偏分散はVzとして定義することができる。したがって、トータル応答可変性totalresponsevariabilityは次のように記述できる。これらの表現は、分析的にあるいは基底関数の態様に応じたモンテカルロ法を介して演算できる。感度指標が無次元のエンティティnon-dimensionalentitiesであるので、あらゆる変数のシステムレベル重要性は、以下のように対応する指標を単純にすべて加えることによって容易に推定することができる。ここでStotalijおよびSijは、j番目の応答に対するi番目の変数のトータル指標および主な感度指標をそれぞれ表し、wjは、j番目の応答に関連した重みを表す。グラフィカルベースの、あるいはグラフィカルプレゼンテーションツール例えば対話型グラフィカルユーザインタフェースは、ユーザがそれぞれの重み因子wjをそれぞれのStotalijおよび/またはSijに対話式に適用することを可能にするように構成される。複雑な数値的結果つまりグローバル感度指標のそのような図表式表現によって、1セットの最も影響力のある設計変数を直観的に容易に判断することができる。グローバル感度指標の例示的なグラフィカル表示を図4に示す。それぞれの設計変数DV-1、DV-2、DV-3、502の目的および制約に対する累積的トータルグローバル感度指標511~515がプロットされている。図4の例において、異なる斜線掛けによって5つの応答関数511~515を示している。応答関数のそれぞれの寄与を、重み因子を適用することによって調整することができる。一の態様では、図5は、複数のスライディングスケール520を示しており、複数のスライディングスケールの1つ1つは応答関数511~515のそれぞれと対応している。スライディングスケール520は、多数の起こりうる事態のシナリオを生成することができるよう、ユーザによって対話式で調整されるように構成される。起こりうる事態のシナリオの結果、1セットの最も影響力のある設計変数が判断される。一般に、最も影響力のある設計変数の数は、1セットの全体の数より著しく少ない。自動車耐衝撃性の例において、総計何百もの設計変数のうちで最も影響力のある設計変数の数は10個ほどである。どの設計変数が重要なあるいは最も影響力のある設計変数として保持すべきかを判断する基準には、以下の方法を挙げることができる。方法1:予め定義されたしきい値より大きな感度指標を有する設計変数を選択する。方法2:次の2段階プロシージャによって設計変数を選択する。a設計変数を感度指標を用いて降順に並べ換える。b累積的な感度指標が予め定義されたしきい値より大きくなるまで、最も高い感度指標を有する設計変数を抽出する。最も影響力のある設計変数が識別されると、ステップ206において、低減された設計空間において1セットの最適設計つまりパレート最適解セットを取得するよう工学設計最適化が実行される。低減された設計空間は、最も影響力のある設計変数の範囲によって定義される。残りの設計変数は、特定の値に固定される。特定の値は、基準値あるいは論理的であると考えられる他の値例えばスクリーニング段階の後に最適解から得られた値とできる。次に図3を参照して、本発明の一の態様にかかる最も影響力のある設計変数を識別する例示的なプロセスを示すフローチャートを示す。ステップ304において、1セットのサンプルあるいは製品設計選択肢が、1セットの全体の設計変数によって定義される設計空間から選択される。選択されたサンプルのそれぞれが、ステップ306において、設計目的および制約に基づいた応答に対して評価される。この評価は、コンピュータ支援工学解析例えば有限要素解析法を用いて達成することができる。言いかえれば、それぞれのサンプルの応答はコンピューターシミュレーションを介して取得される。次に、ステップ308において、近似モデルを用いてすべての設計変数対すべての応答関数のグローバル感度解析が行なわれる例えば、放射基底関数を用いて取得されたサンプル応答を適合フィッティングfittingさせる。累積的トータルグローバル感度指標は、ユーザに図表式で提示され例えば図4、ステップ310において応答の重み因子を調整することによって、ユーザが対話式に起こりうる事態のシナリオを行なうことを可能にする。最後に、ステップ312において、最も影響力のある設計変数が、予め定義された基準例えばパーセンテージしきい値に応じて選択される。自動車耐衝撃性の例を用いれば、10個の最も影響力のある設計変数が、可変性の70以上に寄与している。13個の変数が、最適化に1パーセント以上の影響を与える。なお、グローバル感度解析を、既に低減されている設計空間をさらに低減するために複数回適用することができることを付言しておく。一の側面において、本発明は、ここに説明した機能を実行可能な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を、メインメモリ608に、ロードすることもできる。動作においては、少なくとも1つのプロセッサ604がアプリケーションモジュール606のうちの1つが実行されると、結果が演算されて二次メモリ610つまりハードディスクドライブ612に記憶される。CAE解析あるいは工学設計最適化の状況例えば組み合わせ前後のパレート最適解は、テキストあるいはグラフィック表現で、I/Oインタフェース630を介してユーザに報告される。本発明を具体的な実施形態を参照しながら説明したが、これらの実施形態は単なる例示であって、本発明を限定するものではない。開示した例示的な実施形態に対する種々の変更あるいは変形を、当業者は思いつくであろう。例えば、応答関数の数および設計変数を比較的少数の例示的な数を用いて説明したが、非常に多数の設計変数および目的を用いることができる。また、積み重ね棒グラフをグローバル感度解析の結果を図表式で表すのに示し説明したが、例えば、他のタイプ例えば円グラフを代わりに用いることができる。つまり、発明の範囲は、ここで開示した具体的で例示的な実施形態に限定されず、当業者が容易に想到するあらゆる変更が、本願の精神および認識範囲そして添付の特許請求の範囲の権利範囲に含まれる。"