(57)【要約】【課題】インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを構築する方法およびシステムを提供する。【解決手段】一実施形態は、パラメトリック状態モデルと、設計変数ベクトルとを含み、現実世界のシステムを表すモデルを定義することによって開始する。次いで、第１および第２の実験が実行され、パラメトリック状態ベクトルの経時的な応答を決定し、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成する。データセットは、次いで、パラメトリック状態ベクトルの１または複数の導関数で修正され、サロゲート微分方程式のセットが、データセット内のそれに関してパラメトリック状態ベクトルの高次導関数を近似するように構築され、サロゲート微分方程式のセットは、サロゲートモデルとして格納される。サロゲートモデルは、次いで、モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリから提供される。【選択図】図１
2016146169.tif 000002
【特許請求の範囲】
【請求項１】
インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを生成するコンピュータ実施方法であって、
現実世界のシステムを表すモデルを定義することであって、前記モデルは、振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含むことと、
前記パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験を実行し、かつ、前記設計変数ベクトルに対する第２の実験を実行することであって、前記第１および第２の実験を実行することは、前記パラメトリック状態ベクトルおよび前記設計変数ベクトルのデータセットを時間の関数として生成することと、
前記パラメトリック状態ベクトルの１または複数の導関数で前記データセットを修正することと、
前記データセット内の前記パラメトリック状態ベクトルに関して前記パラメトリック状態ベクトルの高次導関数を近似するサロゲート微分方程式のセットを構築することであって、前記構築することは、前記サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納することを含むことと、
前記モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリから前記サロゲートモデルを応答可能に提供することであって、前記サロゲートモデルを提供することは、前記サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む前記設計変数ベクトルの所与の値に対する前記サロゲート微分方程式を解くことと
を含む方法。
【請求項２】
前記サロゲート微分方程式のセットを構築することは、
前記サロゲート微分方程式の１または複数の候補セットを反復して構築することであって、前記サロゲート微分方程式のセットは、少なくとも前記データセットのサブセットに対する最小交差検証エラーを伴う前記セットとして前記１または複数の候補セットから選択されることを含む、請求項１に記載の方法。
【請求項３】
前記パラメトリック状態ベクトルの１または複数の導関数で前記データセットを修正することは、前記データセットから時間を除去することを含む、請求項１に記載の方法。
【請求項４】
前記第１および第２の実験は、
物理的実験、および、
数値高忠実度実験
のうちの少なくとも１つである、請求項１に記載の方法。
【請求項５】
前記定義されたモデルは、
離散事象、および、
境界条件のうちの少なくとも１つを含む、請求項１に記載の方法。
【請求項６】
前記定義されたモデルは、高忠実度の高次元数値モデルであり、
前記定義されたモデルの次元数を低減すること
をさらに含む、請求項１に記載の方法。
【請求項７】
前記定義されたモデルの次元数を低減することは、
主成分分析、
ｋ−最近傍法、および、
部分空間学習
のうちの少なくとも１つを採用する、請求項６に記載の方法。
【請求項８】
前記サロゲートモデルをファンクショナルモックアップインターフェイス（ＦＭＩ）にエクスポートすること
をさらに含む、請求項１に記載の方法。
【請求項９】
前記サロゲートモデルをハードウェアシステムに採用することであって、前記ハードウェアシステムは、飛行シミュレータまたは自動車電子制御装置（ＥＣＵ）を含むこと
をさらに含む、請求項１に記載の方法。
【請求項１０】
前記パラメトリック状態ベクトルの前記高次導関数は、
放射基底関数、
ニューラルネットワーク、
チェビシェフの多項式、
応答曲面法、
多項式応答曲面法、
任意項回帰、
サポートベクトルマシン、および、
空間マッピング
のうちの少なくとも１つを使用することによって近似される、請求項１に記載の方法。
【請求項１１】
前記サロゲートモデルを構築することは、
分散およびバイアスエラーを低減するために前記サロゲートモデルを最適化すること
をさらに含む、請求項１に記載の方法。
【請求項１２】
インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを生成するシステムであって、
現実世界のシステムを表すモデルを定義するように構成されたモデル定義モジュールであって、前記モデルは、振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含むモデル定義モジュールと、
前記モデル定義モジュールに動作可能に結合され、前記パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験、および、前記設計変数ベクトルに対する第２の実験を実行するように構成される実験モジュールであって、前記第１および第２の実験は、前記パラメトリック状態ベクトルおよび前記設計変数ベクトルのデータセットを時間の関数として生成する実験モジュールと、
前記実験モジュールから前記データセットを受信して、前記パラメトリック状態ベクトルの１または複数の導関数により前記データセットを修正するように構成されたデータセットモジュールと、
前記データセット内の前記パラメトリック状態ベクトルに関して前記パラメトリック状態ベクトルの高次導関数を近似するサロゲート微分方程式のセットを構築するように構成されたサロゲートモジュールであって、前記構築することは、前記サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納することを含む、サロゲートモジュールと、
前記モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリから前記サロゲートモデルを応答可能に提供するように構成されたインタラクションモジュールであって、前記サロゲートモデルを提供することは、前記サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む前記設計変数ベクトルの所与の値に対する前記サロゲート微分方程式を解く、インタラクションモジュールと
を備えたシステム。
【請求項１３】
前記サロゲートモジュールが、
前記サロゲート微分方程式の１または複数の候補セットを反復して構築することであって、前記サロゲート微分方程式のセットは、少なくとも前記データセットのサブセットに対する最小交差検証エラーを伴う前記セットとして前記１または複数の候補セットから選択されること
によって、前記サロゲート微分方程式のセットを構築するように構成された請求項１２に記載のシステム。
【請求項１４】
前記データセットモジュールが、前記パラメトリック状態ベクトルの１または複数の導関数で前記データセットを修正することにおいて、前記データセットから時間を除去するように構成された請求項１２に記載のシステム。
【請求項１５】
前記第１および第２の実験は、
物理的実験、および、
数値高忠実度実験
のうちの少なくとも１つである、請求項１２に記載のシステム。
【請求項１６】
前記定義されたモデルは、
離散事象、および、
境界条件のうちの少なくとも１つを含む、請求項１２に記載のシステム。
【請求項１７】
前記定義されたモデルは、高忠実度の高次元数値モデルであり、前記モデル定義モジュールが、前記定義されたモデルの次元数を低減するように構成された請求項１２に記載のシステム。
【請求項１８】
前記モデル定義モジュールが、
主成分分析、
ｋ−最近傍法、および、
部分空間学習
のうちの少なくとも１つを採用することによって、前記定義されたモデルの次元数を低減するように構成された請求項１７に記載のシステム。
【請求項１９】
前記サロゲートモデルをファンクショナルモックアップインターフェイス（ＦＭＩ）にエクスポートするように構成されたエクスポートモジュール
をさらに備えた請求項１２に記載のシステム。
【請求項２０】
飛行シミュレータまたは自動車電子制御装置（ＥＣＵ）のうちの少なくとも１つに組み入れられる、請求項１２に記載のシステム。
【請求項２１】
前記サロゲートモデルが、
放射基底関数、
ニューラルネットワーク、
チェビシェフの多項式、
応答曲面法、
多項式応答曲面法、
任意項回帰、
サポートベクトルマシン、および、
空間マッピング
のうちの少なくとも１つを使用することによって、前記パラメトリック状態ベクトルの前記高次導関数を近似するように構成された、請求項１２に記載のシステム。
【請求項２２】
１または複数のクライアントとのネットワークにわたる通信においてサーバによって実行されるコンピュータプログラム製品であって、
プログラム命令を備えたコンピュータ可読媒体を備え、前記プログラム命令はプロセッサに実行されると、
現実世界のシステムを表すモデルを定義することであって、前記モデルは、振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含むことと、
前記パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験を実行し、かつ、前記設計変数ベクトルに対して第２の実験を実行することであって、前記第１および第２の実験を実行することは、前記パラメトリック状態ベクトルおよび前記設計変数ベクトルのデータセットを時間の関数として生成することと、
前記パラメトリック状態ベクトルの１または複数の導関数により前記データセットを修正することと、
前記データセット内の前記パラメトリック状態ベクトルに関して前記パラメトリック状態ベクトルの高次導関数を近似するサロゲート微分方程式のセットを構築することであって、前記構築することは、前記サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納することを含むことと、
前記モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリから前記サロゲートモデルを応答可能に提供することであって、前記サロゲートモデルを提供することは、前記サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む前記設計変数ベクトルの所与の値に対する前記サロゲート微分方程式を解くことを含むことと
を実行させる、コンピュータプログラム製品。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は一般に、コンピュータプログラムおよびシステムの分野に関し、具体的には製品設計およびシミュレーションの分野に関する。本発明の実施形態はまた、ビデオゲーム、エンジニアリングシステム設計、協調意志決定、およびたとえば映画のようなエンターテイメントに採用されてもよい。
【背景技術】
【０００２】
いくつかの既存の製品およびシミュレーションシステムが、パーツまたはパーツのアセンブリの設計およびシミュレーション向けに市場に提供されている。そのようなシステムは通常、コンピュータ支援設計（ＣＡＤ）および／またはコンピュータ支援エンジニアリング（ＣＡＥ）プログラムを採用する。これらのシステムは、ユーザが、オブジェクトまたはオブジェクトのアセンブリの複合３次元モデルの構築、操作、およびシミュレーションを行うことができるようにする。したがって、これらのＣＡＤおよびＣＡＥシステムは、エッジまたは辺、特定の場合には面を使用してモデル化されたオブジェクトの表現を提供する。辺、エッジ、面、またはポリゴンは、たとえば非一様有理Ｂスプライン（ＮＵＲＢＳ）のような、さまざまな方法で表現されうる。
【０００３】
これらのＣＡＤシステムは、モデル化されたオブジェクトのパーツまたはパーツのアセンブリを管理するが、これらは主として形状の仕様である。特に、ＣＡＤファイルは仕様を含み、その仕様からジオメトリが生成される。ジオメトリから、表現が生成される。仕様、ジオメトリ、および表現は、単一のＣＡＤファイルまたは複数のＣＡＤファイルに格納されてもよい。ＣＡＤシステムは、モデル化されたオブジェクトを設計者に表現するためのグラフィックツールを含み、これらのツールは、複合オブジェクトの表示に専用である。たとえば、アセンブリは、数千のパーツを含むことができる。ＣＡＤシステムは、オブジェクトのモデルを管理するために使用されてもよく、これらは電子ファイルに格納される。
【０００４】
ＣＡＤおよびＣＡＥシステムの出現は、オブジェクトの広範な表現の可能性を実現させることができる。１つのそのような表現は、有限要素解析（ＦＥＡ）モデルである。ＦＥＡモデル、有限要素モデル（ＦＥＭ）、有限要素メッシュ、およびメッシュという用語は、本明細書において同義的に使用される。ＦＥＭは通常、ＣＡＤモデルを表現し、したがって、１または複数のパーツ、またはアセンブリ全体を表現することができる。ＦＥＭは、メッシュと称されるグリッドを構成するために相互に連結されるノードと呼ばれる点のシステムである。ＦＥＭは、ＦＥＭが基礎となるオブジェクトまたはそれが表すオブジェクトのプロパティを有するようにプログラムされてもよい。ＦＥＭまたは他のそのような、当技術分野において知られているオブジェクト表現が、そのような方法でプログラムされる場合、これはそれが表現するオブジェクトのシミュレーションを実行するために使用されてもよい。たとえば、ＦＥＭは、車両の内部空間、構造体を取り巻く音響流体、および任意の数の現実世界のオブジェクトを表現するために使用されてもよい。さらに、ＣＡＤおよびＣＡＥシステムはＦＥＭと共に、エンジニアリングシステムをシミュレートするために利用されてもよい。たとえば、ＣＡＥシステムは、車両のノイズおよび振動をシミュレートするために採用されてもよい。
【０００５】
しかし、そのような既存のシミュレーション方法にも、欠点がないわけではない。内容が参照により本明細書に援用される、ｈｔｔｐ：／／ｗｗｗ．３ｄｓ．ｃｏｍ／ｐｒｏｄｕｃｔｓ−ｓｅｒｖｉｃｅｓ／ｓｉｍｕｌｉａ／ｏｖｅｒｖｉｅｗ／において入手可能な非特許文献１、および非特許文献２において説明されているような、現実的な試験的および多数の実験は、多大な時間を要し、実行するために膨大な専門的知識を必要とする。そのようなシミュレーション／実験からの典型的な結果は、ビデオ、イメージ、またはグラフの形態で提供される。非特許文献３において提示されている調査のような、行動学的調査は、リアルタイムの実験を通じて得られることを示している。アクションが短時間間隔で直接の結果を生じる場合、「エクスペリエンス」が得られる。長期的な結果を有するアクションは、小規模なグループの人々に対する「エクスペリエンス」しかもたらさない。同様に、対話性が全くない（つまり、実験／シミュレーションを交換する機会がない）場合、ほとんどエクスペリエンスは得られない。したがって、たとえば、イメージを通じて科学的な結果を提示する専門家の現在のパラダイムは、オーディエンスのエクスペリエンスを増大させるための手段としては有効ではない。
【０００６】
既存のシミュレーション／実験に関するもう１つの問題は、非特許文献２におけるように、シミュレーション自体を実行し、次いで結果を再使用するために要する時間であり、このことは非特許文献４において説明されている。そのような時間を浪費するシミュレーションは、数百万もの検証実験を必要とするシステム検証または最適化には使用され得ない。そのようなシステム研究を行うために、現在の現実的なシミュレーションツール（たとえば、非線形有限要素解析）よりも、少なくとも３桁から４桁も高速な技法が必要となる。
【０００７】
シミュレーション／実験のユーザ「エクスペリエンス」を高めるための試みが行われてきたが、既存の技法では不十分である。特許文献５において説明されているような、広範に使用される一連の作業は、特許文献２において説明されているシミュレーションされた振る舞いを作り出すための多数のモデルを使用する実験、および／または特許文献４において説明されている実験的試験に関して使用可能である。これらの実験／シミュレーションは、多くの場合、応答曲面法、チェビシェフの多項式、クリギング、または放射基底関数のような技法を使用して補間される。これらの方法は、固定入力、１次元（１Ｄ）、パラメータのセットと、固定出力、１Ｄ、出力パラメータのセットとの間の正確な伝達関数を作り出す。１次元データは、極めて抽象的であり、抽象的情報は、オーディエンスのエクスペリエンスを増大させるための効率的な手段ではない。そのような方法は、シミュレーションの状態空間を縮小することに基づいている。そのような例において、完全３Ｄ計算（すなわち、メッシュ）モデルが先見的に知られているが、知られている振る舞いのモデルへの良好な適合を維持しつつ、モデルの複雑さと解決時間が短縮されている。これは通常、特許文献６において説明されている主成分分析のような技法を通じて、自由度を低減すること（たとえば、数百万のメッシュノードから数十のパラメータへ）により行われる。
【０００８】
近年、適正な直交分解によるシロビッチの「スナップショット法」を使用してフィールドデータの特殊補間を作成する取り組みが行われた。この方法は、エクスペリエンスの抽象性を減じる（基底）フィールドの小規模なセットの（線形）関数として、２次元（２Ｄ）または３Ｄ（フィールド）データを複製する。しかし、この方法は、所与のモデルの時系列シミュレーションの圧縮にしか適しておらず、したがって、シナリオの一部としてモデルの変動を考慮することはない。加えて、個々のセンサーデータ（フィールド内の場所）の予測値は、フィールド全体を近似する際のトレードオフにより、あまり正確ではない。
【０００９】
ｈｔｔｐ：／／ｗｗｗ．ｓｃａｉ．ｆｉ−ａｕｎｈｏｆｅｒ．ｄｅｌｇｅｓｃｈａｅｆｔｓｆｅｌｄｅｒｉｎｕｍｅｒｉｓｃｈｅ−ｓｏｆｔｗａｒｅ／ｐｒｏｄｕｋｔｅ／ｆｅｍｚｉｐ．ｈｔｍｌ？＆Ｌ＝１において入手可能な、特許文献７において説明されているような、シミュレーションフィールドデータの圧縮のための既存のツールがある。これらのツールは、（シロビッチにより説明される）基底フィールドが表示のために計算される速度を向上させる。そのようなツールはまた、より広い配布を可能にするためにデータストレージ要件を軽減するであろう。しかし、単独では、これらのツールは、幅広いオーディエンスに「エクスペリエンス」を生成するために、シナリオおよびモデルが「オンザフライ」で変更されうるようにはしない。
【００１０】
もう１つの既存の方法は、パラメータ推定に依拠するが、これはＫｌｅｉｎにおいて説明されている。そのような方法の実行においては、微分方程式が知られているが、方程式の定数の一部は知られていない。たとえば航空機力学の場合、未知の空気力学導関数は、モーションの操作の既知のセットに対する飛行データと航空機モデルデータとの間のエラーを最小化することにより見出されうる。
【００１１】
特許文献８
【００１２】
たとえば、クラインの方程式（３１）は、以下のとおりである。
次いで、調和運動について
【００１３】
【数１】
2016146169.tif 000003
【００１４】
ここでは、項ＣＬα、ＣＬｑ、
【００１５】
【数２】
2016146169.tif 000004
【００１６】
、およびｋは、ｔ、ＣＬ、ｑ、α、
【００１７】
【数３】
2016146169.tif 000005
【００１８】
の実験データに基づいて式に合わせて較正または適合されることになる。
【００１９】
多くの新規の製品設計シミュレーションに対して、モーションの方程式は、本明細書において示される航空機の方程式ほどよく知られてはいない。したがって、多くのシミュレーションは、上記で説明される方法に依拠することはできない。パラメータ推定に依拠する方法は、パラメトリック状態の数を、高階数学的サロゲートについて管理可能であるレベルまで減じるので、Ｊｏｌｌｉｆｆｅにおいて論じられる手法に補完的である。
【００２０】
本発明の実施形態は、上記において説明されるような既存の実験／シミュレーション技法の問題を克服する。
【先行技術文献】
【非特許文献】
【００２１】
【非特許文献１】Ｎｕｍｅｒｉｃａｌ  Ｍｏｄｅｌｓ  ｔｏ  Ｃｒｅａｔｅ  Ｓｉｍｕｌａｔｅｄ  Ｂｅｈａｖｉｏｒ
【非特許文献２】Ｈｅｌｍ、Ｅｘｐｅｒｉｍｅｎｔａｌ  Ｔｅｃｈｎｉｑｕｅｓ
【非特許文献３】Ｋａｈｎｅｍａｎ、Ｔｈｉｎｋｉｎｇ  Ｆａｓｔ  ａｎｄ  Ｓｌｏｗ
【非特許文献４】Ｎｕｍｅｒｉｃａｌ  Ｍｏｄｅｌｓ  ｔｏ  Ｃｒｅａｔｅ  Ｓｉｍｕｌａｔｅｄ  Ｂｅｈａｖｉｏｒ  ａｎｄ  Ｅｘｐｅｒｉｍｅｎｔａｌ  Ｔｅｃｈｎｉｑｕｅｓ
【非特許文献５】Ｂｏｘ  ｅｔ  ａｌ．、Ｓｔａｔｉｓｔｉｃ  ｆｏｒ  Ｅｘｐｅｒｉｍｅｎｔｅｒｓ：Ｄｅｓｉｇｎ，Ｉｎｎｏｖａｔｉｏｎ，ａｎｄ  Ｄｉｓｃｏｖｅｒｙ
【非特許文献６】Ｊｏｌｌｉｆｆｅ  Ｉ．Ｔ．、Ｐｒｉｎｃｉｐａｌ  Ｃｏｍｐｏｎｅｎｔ  Ａｎａｌｙｓｉｓ、Ｓｅｒｉｅｓ：Ｓｐｒｉｎｇｅｒ  Ｓｅｒｉｅｓ  ｉｎ  Ｓｔａｔｉｓｔｉｃｓ、２ｎｄ  ｅｄ．、Ｓｐｒｉｎｇｅｒ、ＮＹ
【非特許文献７】Ｓｉｒｏｖｉｃｈ、Ｔｕｒｂｕｌｅｎｃｅ  ａｎｄ  ｔｈｅ  Ｄｙｎａｍｉｃｓ  ｏｆ  Ｃｏｃｈｅｒｅｎｔ  Ｓｔｒｕｃｔｕｒｅｓ  Ｉ−Ｉｌｌ  ａｎｄ  Ｆｒａｕｎｈｏｆｅｒ
【非特許文献８】Ｋｌｅｉｎ  Ｖ．  ”Ａｅｒｏｄｙｎａｍｉｃ  Ｐａｒａｍｅｔｅｒｓ  ｏｆ  Ｈｉｇｈ  Ｐｅｒｆｏｒｍａｎｃｅ  Ａｉｒｃｒａｆｔ  Ｅｓｔｉｍａｔｅｄ  ｆｒｏｍ  Ｗｉｎｄ  Ｔｕｎｎｅｌ  ａｎｄ  Ｆｌｉｇｈｔ  Ｔｅｓｔ  Ｄａｔａ”、ＮＡＳＡ−９８−ＡＧＡＲＤ
【非特許文献９】Ｖａｎ  ｄｅｒ  Ｖｅｌｄｅｎ  ｅｔ  ａｌ．、”Ｐｒｏｂａｂｉｌｉｓｔｉｃ  Ｃｅｒｔｉｆｉｃａｔｅ  ｏｆ  Ｃｏｒｒｅｃｔｎｅｓｓ  ｆｏｒ  Ｃｙｂｅｒ  Ｐｈｙｓｉｃａｌ  Ｓｙｓｔｅｍｓ、”ＡＳＭＥ  ２０１２  ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｄｅｓｉｇｎ  Ｅｎｇｉｎｅｅｒｉｎｇ  Ｔｅｃｈｎｉｃａｌ  Ｃｏｎｆｅｒｅｎｃｅｓ  ＆  Ｃｏｍｐｕｔｅｒｓ  ａｎｄ  Ｉｎｆｏｒｍａｔｉｏｎ  ｉｎ  Ｅｎｇｉｎｅｅｒｉｎｇ  Ｃｏｎｆｅｒｅｎｃｅ、Ａｕｇｕｓｔ  １２−１５、Ｃｈｉｃａｇｏ、ＩＬ、ＵＳＡ、ＤＥＴＣ２０１２−７０１３５
【非特許文献１０】Ｈａｒｄｙ、”Ｍｕｌｔｉｑｕａｄｒａｔｉｃ  ｅｑｕａｔｉｏｎｓ  ｏｆ  ｔｏｐｏｇｒａｐｈｙ  ａｎｄ  ｏｔｈｅｒ  ｉｒｒｅｇｕｌａｒ  ｓｕｒｆａｃｅｓ”、Ｊ．Ｇｅｏｐｈｙｓｉｃｓ  Ｒｅｓ．７６、１９０５−１５（１９７１）
【発明の概要】
【課題を解決するための手段】
【００２２】
本発明の実施形態は、インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを構築するための方法およびシステムを提供する。少なくとも１つの実施形態によれば、インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを生成するための方法は、現実世界のシステムを表すモデルを定義するステップによって開始し、モデルは振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含む。そのような方法は、パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験（シミュレーション、分析など）を実行するステップと、設計変数ベクトルに対する第２の実験を実行するステップとをさらに含む。第１および第２の実験は、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成する。方法は、パラメトリック状態ベクトルの１または複数の導関数でデータセットを修正するステップによって続行し、次いで、データセット内のパラメトリック状態ベクトルに関してパラメトリック状態ベクトルの高次導関数を近似するサロゲート微分方程式のセットを構築する。そのような実施形態において、構築するステップは、サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納するステップを含む。そのような方法は、モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する（accelerating）方法で、メモリからサロゲートモデルを応答可能に提供することにより、インタラクティブエクスペリエンスにおいて使用される。これは、サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む設計変数ベクトルの所与の値に対するサロゲート微分方程式を解くことによって実施されてもよい。
【００２３】
方法の代替的実施形態において、サロゲート微分方程式のセットを構築するステップは、サロゲート微分方程式の１または複数の候補セットを反復して構築するステップと、サロゲートモデルとして格納されるべき１または複数の候補セットからサロゲート微分方程式のセットを選択するステップとを備え、選択されるセットは、データセットの少なくともサブセットに対する最小交差検証エラーを伴うセットである。そのような実施形態において、データセットのサブセットは、本明細書において検証データセットと称されてもよい。もう１つの実施形態によれば、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正するステップは、データセットから時間を除去するステップを含む。さらにもう１つの実施形態において、第１および第２の実験は、物理的実験および数値高忠実度実験のうちの少なくとも１つである。もう１つの実施形態は、定義済みモデルの次元数を低減するステップをさらに備えることができる。そのような実施形態において、定義済みモデルの次元数を低減するステップは、主成分分析、ｋ−最近傍法、および部分空間学習のうちの少なくとも１つを採用する。
【００２４】
さらに、さらなるもう１つの実施形態において、定義済みモデルは、離散事象および境界条件のうちの少なくとも１つを含む。方法の実施態様は、サロゲートモデルをファンクショナルモックアップインターフェイスにエクスポートするステップをさらに備えることができる。さらになお、実施形態は、サロゲートモデルをハードウェアシステムに採用することができ、ハードウェアシステムは飛行シミュレータまたは自動車電子制御装置である。
【００２５】
本明細書において説明される原理に従って実行される、本発明の１つの実施形態によれば、パラメトリック状態ベクトルの高次導関数は、放射基底関数、ニューラルネットワーク、チェビシェフの多項式、応答曲面法、多項式応答曲面法、任意項回帰（arbitrary term regression）、サポートベクトルマシン、および空間マッピングのうちの少なくとも１つを使用することによって近似されてもよい。本発明の方法によれば、サロゲートモデルを構築するステップは、分散およびバイアスエラーを低減するためにサロゲートモデルを最適化するステップをさらに含む。
【００２６】
本発明の代替的実施形態は、インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを生成するためのシステムを対象とする。そのようなシステムは、現実世界のシステムを表すモデルを定義するように構成されたモデル定義モジュールを備え、定義済みモデルは、振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含む。本発明の実施形態によれば、システムは、モデル定義モジュールに動作可能に結合され、パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験、および設計変数ベクトルに対する第２の実験を実行するように構成される実験モジュールをさらに含む。実験は、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成する。実施形態において、システムは、実験モジュールからデータセットを受信して、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正するように構成されたデータセットモジュールをさらに含む。サロゲートモジュールは、システムにさらに含まれ、サロゲートモジュールは、パラメトリック状態ベクトルの高次導関数をデータセット内のそれに関して近似するサロゲート微分方程式のセットを構築するように構成される。さらに、サロゲートモジュールは、サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納することができる。システムは、ユーザインタラクションに応答してモデルのシミュレートされる振る舞いを促進する方法で、メモリからサロゲートモデルを応答可能に提供するように構成されたインタラクションモジュールをさらに含む。前記サロゲートモデルを提供するステップは、サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む設計変数ベクトルの所与の値に対するサロゲート微分方程式を解く。
【００２７】
システムの実施形態によれば、サロゲートモジュールは、サロゲート微分方程式の１または複数の候補セットを反復して構築することによって、サロゲート微分方程式のセットを構築するように構成される。そのような実施形態において、サロゲートモデルとして格納されるサロゲート微分方程式のセットは、少なくともデータセットのサブセットに対する最小交差検証エラーを伴うセットとして１または複数の候補セットから選択される。実施形態において、データセットモジュールは、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正するように構成される。さらにもう１つの実施形態において、第１および第２の実験は、物理的実験および数値高忠実度実験のうちの少なくとも１つである。システムのさらなるもう１つの実施形態において、定義済みモデルは、離散事象および境界条件のうちの少なくとも１つを含む。
【００２８】
さらに、定義済みモデルが高忠実度の高次元数値モデルである実施形態において、モデル定義モジュールは、定義済みモデルの次元数を低減するようにさらに構成される。そのような実施形態において、モデル定義モジュールは、主成分分析、ｋ−最近傍法、および部分空間学習のうちの少なくとも１つを採用することによって、定義済みモデルの次元数を低減するように構成されてもよい。システムの代替的実施形態は、サロゲートモデルをファンクショナルモックアップインターフェイスにエクスポートするように構成されたエクスポートモジュールをさらに備える。さらにもう１つの実施形態において、システムは、飛行シミュレータまたは自動車電子制御装置に組み入れられる。
【００２９】
サロゲートモジュールは、放射基底関数、ニューラルネットワーク、チェビシェフの多項式、応答曲面法、多項式応答曲面法、任意項回帰、サポートベクトルマシン、および空間マッピングのうちの少なくとも１つを使用してパラメトリック状態ベクトルの高次導関数を近似するように構成されてもよい。
【００３０】
本発明のもう１つの実施形態は、サロゲートモデルを生成するためのクラウドコンピューティング実施を対象とする。そのような実施形態は、１または複数のクライアントとのネットワークにわたる通信においてサーバによって実行されるコンピュータプログラム製品を対象とし、コンピュータプログラム製品はコンピュータ可読媒体を備える。コンピュータ可読媒体は、プロセッサに実行されるとき、現実世界のシステムを表すモデルを定義させ、モデルは振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含み、パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験を実行させ、設計変数ベクトルに対して第２の実験を実行させ、第１および第２の実験を実行させることは、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成し、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正させ、パラメトリック状態ベクトルの高次導関数をデータセット内のそれに関して近似するサロゲート微分方程式のセットを構築させ、方程式のセットを構築させることは、サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納することを含み、モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリからサロゲートモデルを応答可能に提供させ、サロゲートモデルを提供させることは、サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む設計変数ベクトルの所与の値に対するサロゲート微分方程式を解くことを含む、プログラム命令を備える。
【図面の簡単な説明】
【００３１】
前述の内容は、さまざまな図表全体を通じて類似した参照符号が同じ部分を指す添付の図面に示されているように本発明の例示的な実施形態のさらに詳細な以下の説明から明らかとなろう。図面は、必ずしも一定の縮尺では示されていないが、本発明の実施形態を説明することに重点が置かれている。
【図１】本発明の実施形態によるインタラクティブ３Ｄエクスペリエンスにおいて使用するサロゲートモデルを構築するための方法を示す流れ図である。
【図２】本発明の原理を実施する翼フラッタのＮコードシミュレーションにおいて使用されるさまざまなコンポーネントを示す図である。
【図３】実施形態の原理を利用して決定される制御システム作動中のフラップのＣＦＤフローフィールドを示す図である。
【図４】本発明の実施形態において採用されうる、航空機のＣＡＤモデル、ならびに迎え角および制御サーボ角を説明する関連するグラフを示す図である。
【図５】時間の関数としてシミュレートされる翼ピッチを示すグラフである。
【図６】例示的実施形態により決定される２階常微分方程式を表す線形サロゲートモデルの検証を示す図である。
【図７】例示的実施形態により決定される非線形２階常微分方程式を表す線形サロゲートモデルの検証を示す図である。
【図８】実施形態の原理によるサロゲートモデルを生成する方法を示す流れ図である。
【図９】コントローラオフの４階ＲＳＭおよび４．８秒および３３ＨｚにおけるコントローラオンのＲＢＦで近似される、実施形態の原理により決定される翼ピッチおよび翼垂直変位を示す図である。
【図１０】コントローラオフの４階ＲＳＭおよび３．１秒および４２ＨｚにおけるコントローラオンのＲＢＦで近似される、翼フラッタのＮコードシミュレーションからの時系列を示すグラフである。
【図１１ａ】本発明の原理によるサロゲートモデルを生成するための方法を示す流れ図である。
【図１１ｂ】図１１ａに示されるサロゲートモデルを生成するための方法の続きを示す図である。
【図１２】シミュレーションにおいて本発明の原理を採用するための方法を示す流れ図である。
【図１３】実施形態の原理による３Ｄエクスペリエンスのサロゲートモデルを生成するためのシステムを示す簡略化ブロック図である。
【図１４】本発明の実施形態が実施されうるコンピュータネットワーク環境を示す簡略図である。
【発明を実施するための形態】
【００３２】
これ以降、本発明の例示的実施形態の説明が続く。
【００３３】
すべての特許の教示、公開された出願、および本明細書において引用される文献は、内容全体が参照により本明細書に組み入れられる。
【００３４】
本発明の実施形態は、シミュレーション／インタラクティブ３Ｄエクスペリエンスの前述の制限を克服する。本発明の１つのそのような実施形態は、数学的サロゲートを使用して現実的な振る舞いのシミュレーションを促進することによって、ほぼリアルタイムのエクスペリエンスを作り出すための方法を対象とする。これらのサロゲートモデルは、物理テストデータまたは高忠実度シミュレーションデータに基づく。高忠実度Ｎコード（ＦＥＡ、計算流体力学（ＣＦＤ）、または論理）協調シミュレーションは、振る舞い予測のリアルタイム秒ごとにＣＰＵ時間の１時間にも相当する時間を要することもある。本発明の実施形態は、リアルタイムまたはほぼリアルタイムのアプリケーションに使用されるべく十分な速さになるよう、このプロセスを何桁も高速化する。これは、それを、インタラクティブ設計、多くの専門分野にわたる最適化、およびサイバーフィジカルシステムの検証のための実践的技術にする。これらのサロゲートモデルは、ファンクショナルモックアップユニットとしてラップされてもよく、システムシミュレーションツールへの配備および／またはハードウェアへの組込みが行われてもよい。ＦＭＵは、ｈｔｔｐ：／／ｗｗｗ．ｆｕｎｃｔｉｏｎａｌ−ｍｏｃｋｕｐ−ｉｎｔｅｒｆａｃｅ．ｏｒｇ／ｉｎｄｅｘ．ｈｔｍｌにおいて説明されており、その内容は参照により本明細書に組み入れられる。
【００３５】
実施形態は、最初に、物理または数値モデルに物理または数値高忠実度実験を実行し、関心対象の振る舞いを表すパラメトリック状態ベクトルＰの経時的な応答（たとえば、時間、Ｐ）を観察することによって、現実的な振る舞いのモデル化を促進する。このモデルは、離散事象（たとえば、コントローラのオン／オフ）、ならびに道路励振のような強制境界条件を含むことができる。あるいは、そのような実験を実行する前に、実施形態は、最初に、主成分分析、または類似する技法を使用して、高忠実度の、高次元数値モデルの次元数を管理可能なパラメータの数（たとえば、約１００）まで低減することができる。続いて、実験を実行した後、そのような実施形態は次に、この実験を、関心対象のモデル設計変数ベクトルｖについて構造化された（たとえば、実験の設計）、または構造化されていない方法で反復して、初期条件ｐ（ｔｉｍｅ＝０）の所与のセットについて（ｔｉｍｅ、ｖ、ｐ）の関数としてデータセットを生成する。
【００３６】
方法は、データセットを結合し、後退微分方式を使用して状態変数をそれぞれ時間に関して微分し、次いでデータセットから時間を除去することによって続行し、
【００３７】
【数４】
2016146169.tif 000006
【００３８】
をもたらす。次いで、低次導関数のサロゲート関数としての高次導関数
【００３９】
【数５】
2016146169.tif 000007
【００４０】
は、放射基底関数、ニューラルネットワーク、チェビシェフの多項式、および応答曲面法のような技法を、
【００４１】
【数６】
2016146169.tif 000008
【００４２】
の完全セットの一部であるトレーニングセットにわたる項低減と共に使用して近似される。正しい予測サロゲートが選択され、その技法設定は、その分散およびバイアスエラーを低減するように最適化され、トレーニングに使用されないデータセットの部分で検証される。このようにして、サロゲートモデルを形成する。実施形態はまた、所与の値ｖについてｔｉｍｅ＝０からｐ（ｔ）を数値前進積分し、初期条件ｐ（ｔｉｍｅ＝０）および
【００４３】
【数７】
2016146169.tif 000009
【００４４】
を与えることができる。状態ベクトルおよび設計値のこれらの初期条件は、元の時系列の実験のものと同じであるように制限されることはない。このモデルは、次には、ＦＭＵにエクスポートされ、ほぼリアルタイムで現実的な振る舞いをシミュレートするように３Ｄモデリングソフトウェアに直接採用されてもよい。あるいは、これは、飛行シミュレータまたは自動車電子制御装置（ＥＣＵ）のようなハードウェアに直接使用されてもよい。
【００４５】
図１は、インタラクティブエクスペリエンスにおいて使用するサロゲートモデルを生成するための方法１００を示す流れ図である。方法１００は、現実世界のシステムを表すモデルを定義することによって開始する（１０１）。定義済みモデルは、振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含む。たとえば、パラメトリック状態ベクトルは、オブジェクトの位置および回転の角度（およびそれらの導関数）を含むことができ、設計変数ベクトルは、オブジェクトの質量および力（およびその方向）を備えることができる。ステップ１０１において定義されるモデルは、現実世界のシステムを表すことができ、その例については本明細書において後段で説明される。さらに、モデルは、たとえばＣＡＤシステムの使用を通じてなど、当技術分野において知られている原理に従って定義されてもよい。方法１００は、パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験を実行し、設計変数ベクトルに対する第２の実験を実行することによって続行する（１０２）。例示的な実験は、さまざまな質量、力、および力の方向についてパラメトリック状態の時系列を記録することを含む。実行１０２は、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成する。実験１０２を実行した後、方法１００は次に、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正する（１０３）。実施形態によれば、ステップ１０３は、有限差分法を使用して、時間ステップから直接計算されるパラメトリック状態ベクトルおよび時間ステップにおける状態ベクトル値の１または複数の導関数を決定するステップを含む。次に、パラメトリック状態ベクトルの高次導関数をデータセット内のそれに関して近似するサロゲート微分方程式のセットが構築される（１０４）。方法１００のコンピュータ実施の実施形態において、サロゲート微分方程式のセットを構築するステップ１０４は、サロゲート微分方程式のセットをサロゲートモデルとしてメモリに動作可能に格納するステップを備える。方法１００は、モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリからサロゲートモデルを応答可能に提供するステップによって終了する（１０５）。サロゲートモデルを提供するステップ１０５において、サロゲート微分方程式は、サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む設計変数ベクトルの所与の値に対して解かれる。サロゲートモデルを提供するこのステップ１０５は、したがって、改良されたユーザエクスペリエンスを提供する。
【００４６】
方法１００の実施形態によれば、サロゲート微分方程式のセットを構築するステップ１０４は、サロゲート微分方程式の１または複数の候補セットを反復して構築し、次いで、データセットの少なくともサブセットに対する最小交差検証エラーを伴う１または複数の候補セットからセットとしてサロゲートモデルとして格納されているサロゲート微分方程式のセットを選択することによって実行されてもよい。方程式構築を実行する際に使用されるデータセットのこのサブセットは、検証データセットと見なされてもよい。方法１００のさらにもう１つの実施形態によれば、データセットを修正するステップ１０３は、データセットから時間を除去するステップを含むことができる。さらにもう１つの実施形態において、方法１００のステップ１０２において実行される第１および第２の実験は、物理的実験および／または数値高忠実度実験である。
【００４７】
方法１００のさらにもう１つの実施形態において、定義済みモデルは、離散事象および境界条件をさらに含む。定義済みモデルが高忠実度の高次元数値モデルである実施形態において、方法１００は、定義済みモデルの次元数を低減するステップをさらに備えることができる。そのような実施形態によれば、次元数は、たとえば、主成分分析、ｋ−最近傍法、および／または部分空間学習を採用することによって、当技術分野において知られている原理に従って低減されてもよい。
【００４８】
方法１００はまた、サロゲートモデルをファンクショナルモックアップインターフェイスにエクスポートするステップを含むことができる。代替的実施態様において、方法１００は、飛行シミュレータまたは自動車電子制御装置のような、ハードウェアシステムに採用される。方法１００の実施形態において、パラメトリック状態ベクトルの高次導関数は、放射基底関数、ニューラルネットワーク、チェビシェフの多項式、応答曲面法、多項式応答曲面法、任意項回帰、サポートベクトルマシン、および空間マッピングのうちの少なくとも１つを使用してステップ１０４において決定される。さらにもう１つの実施形態において、ステップ１０４においてサロゲートモデルを構築するステップは、分散およびバイアスエラーを低減するためにサロゲートモデルを最適化するステップをさらに含む。
【００４９】
図２は、本発明の原理に従って実行される、インタラクティブエクスペリエンスが基づきうる個々の応答時系列データを作成するための方法２１０の概略を示す図である。後段において説明され、図２において示される例は、内容が参照により本明細書に組み入れられる非特許文献９において説明されるような４コードＣＦＤ−ＦＥＡ−論理組み込みソース協調シミュレーションのデータセットを採用する。これらのデータベースは、いずれもモデル抽象化２１１から生じる常微分方程式／１次元データセット２１３ａおよび偏微分方程式／３Ｄデータセット２１３ｂ、ならびに計算のモデル２１２から生じる連続時間データセット２１４ａおよび離散時間データセット２１４ｂを含む。データベース２１３ａ、２１３ｂ、２１４ａ、および２１４ｂは、非限定的な一例として、Ｄｙｍｏｌａ２１６ａ、Ａｂａｑｕｓ／ＣＦＤ２１６ｂ、Ａｂａｑｕｓ／Ｓｔａｎｄａｒｄ２１６ｃ、およびＣｏｎｔｒｏｌＢｕｉｌｄ２１６ｄを含むことができるドメインシミュレーションソフトウェア２１５によって処理される。これらのデータセット２１３ａ、２１３ｂ、２１４ａ、および２１４ｂは、本明細書において説明される原理に従ってインタラクティブに解決されうるモデルに変換されてもよい。図２に示される概要は、図３に示される有限要素モデル３２１についての、陰影付けによって示される、ＣＦＤフローフィールド３２０のデータをユーザに提供することができる。そのような例において、フローフィールドは、制御システムの作動中に部分的に分離されたフラップにナビエ−ストークスを採用することによって決定されてもよい。
【００５０】
図４は、翼４３０を備えるグライダー４３１を示す。翼４３０は、フラッタと呼ばれるねじれおよび垂直方向変位のリミットサイクル運動を呈する。所与の時間ｔにおいて、翼４３０の後縁フラップを動かすアクティブフラッタ制御システムが導入される。このシステムは、フラップ機構を作動させるために、センサー、組込み型ソフトウェア、および電気モーターで構成される。本明細書において説明される例示の方法は、翼４３０のフラッタをモデル化するインタラクティブエクスペリエンスを提供する。翼４３０上の所与の翼長方向の場所について、制御ソフトウェア作動４３３、および迎え角４３４を指示する、制御サーボ角４３２のような、高いアスペクト比データが、翼４３０の状態に対して提供される。その他のデータは、垂直方向変位を含むことができる。
【００５１】
図５は、時間の関数として翼ピッチαを示すグラフ５４０である。グラフ５４０は、フリーリミットサイクル運動の高忠実度協調シミュレーション（フラッタ抑止作動なし）により時間の関数としてシミュレートされるような翼ピッチを表すデータ５４１を含む。完全協調シミュレーションデータ５４１は、良好な精度をもたらすが生成には著しく時間を浪費して計算処理上高価となりうる、既存のバッチ処理方法を使用して決定されるデータを表す。したがって、本発明の実施形態は、完全協調シミュレーションデータ５４１を近似するサロゲートを計算する。
【００５２】
翼ピッチデータα、および垂直変位データｚから開始し、実施形態は、続いて後退数値ステップを使用してデータセットを微分し、状態パラメータおよびそれらの時間導関数のデータのセット
【００５３】
【数８】
2016146169.tif 000010
【００５４】
をもたらす。次に、時間はパラメータとして除去され、連立方程式（system of equations）が、協調シミュレーションデータの部分を使用して状態ベクトル（α，ｚ）の関数としてαおよびｚの高次導関数およびその低次導関数を表すために、サロゲートＦに基づいて作成される。残りのデータは、サロゲートを検証するために確保される。結果として得られる連立方程式は、以下に示される。
【００５５】
【数９】
2016146169.tif 000011
【００５６】
引き続き、Ｆ１およびＦ２は、よく知られているスプリングダンパーモデルを表す２階線形常微分方程式を得るための線形回帰サロゲートモデルとして表されてもよい。このスプリングダンパーモデルは、これが翼フラッタに適切なモデル形式であるという予備知識なくして導かれることに留意されたい。したがって、本発明の実施形態は、システムを表すための適切なモデルの予備知識を必要としない。この例のデータを使用して導かれた２階線形常微分方程式の例は、以下に示される。
【００５７】
【数１０】
2016146169.tif 000012
【００５８】
２階微分方程式を表す、結果として得られる１階サロゲートモデルは、図５の破線５４２によって示される。図５に示されるように、２階モデルは、完全協調シミュレーションデータ５４１と類似する周波数を有するが、振幅は全く異なる。グラフ５４０は、１階サロゲートで近似された翼フラッタＮコードシミュレーションを表すデータ５４２、および４階サロゲートを表すデータ５４３をさらに含む。同様に、フラッタ制御システムの所与の帯域幅ｗに対する時間の関数として、翼垂直変位ｚの図がある。そのようなデータは、図５において示されるデータと同様に表されうる。
【００５９】
本発明の実施形態は、前述の除外されているデータによりモデルを検証することによって、そのようなエラーのソースをさらに決定することができる。そのような評価は、実際の
【００６０】
【数１１】
2016146169.tif 000013
【００６１】
データと線形サロゲートから構築された２階ＯＤＥによって予測されるものとの間に大きいエラーがあることを指示する。このエラーは図６に示され、ここでグラフ６５０は
【００６２】
【数１２】
2016146169.tif 000014
【００６３】
データを示し、線６５２がモデルの完全な検証を表し、データポイント６５３が除外された実際のデータを示す。この検証から、線形サロゲートがリミットサイクル協調シミュレーションへの良好な近似ではないことが演繹的に知られる。これは、高角度における翼フロー分離の非線形に起因しうる。図６は、グラフ６５１において、
【００６４】
【数１３】
2016146169.tif 000015
【００６５】
近似が極めて良好であることをさらに示す。線６５４が予測される
【００６６】
【数１４】
2016146169.tif 000016
【００６７】
データである場合、データポイント６５５は、以前除外された残りのデータポイントである。
【００６８】
このエラーを訂正するために、実施形態は、サロゲートモデルを２階モデルから４階モデルに切替え、４階モデルを使用して検証を繰り返すことができる。この検証の結果は図７において示され、ここで
【００６９】
【数１５】
2016146169.tif 000017
【００７０】
データおよび
【００７１】
【数１６】
2016146169.tif 000018
【００７２】
データはそれぞれ、グラフ７６０および７６１によって示されている。グラフ７６０および７６１はいずれも、予測されたデータ７６２および７６４が、対応するデータポイント７６３および７６５と極めて良好に一致することを示している。この階モデルの結果は、元の協調シミュレーション振る舞い結果５４１と一致するドットポイント５４３によって図５において示され、良好に外挿する。したがって、そのような実施形態は、完全協調シミュレーションデータの使用による計算上のコストを生じることなくインタラクティブエクスペリエンスに使用されうるサロゲートモデルを展開してきた。
【００７３】
そのような方法の実施態様は図８に示され、ここで図７に示される検証モデルが、事象ｔｉｍｅ＝Ｔａｓｋ．ｃｏｎｔｒｏｌｌｅｒ＿ｏｎ＿ｔｉｍｅまでリミットサイクルモーションに対して方法８７０において使用される。方法８７０は、選択されたパラメータの値のシーケンスを、この場合は、タスクプロセス８７２のコンポーネントであるコントローラオン時間８７４まで、反復するループ８７１を備える。ループ８７１の各反復において、タスクプロセス８７２は、時間が、コントローラオフ時間８７３よりも小さく、コントローラがオフであることを示すか、またはコントローラオン時間８７４よりも大きく、コントローラがオンであることを指示するかどうかを識別する。時間がコントローラオン時間よりも小さい、つまりコントローラがオフである場合８７３、上記で説明される４階サロゲートが、プロセス８７５においてαおよびｚを計算するために使用される。上記で説明されるＦ１およびＦ２を備える連立方程式のインタラクティブモデル８７６は、連立方程式に異なる入力を使用して
【００７４】
【数１７】
2016146169.tif 000019
【００７５】
および
【００７６】
【数１８】
2016146169.tif 000020
【００７７】
を計算するために、システム８７０に採用されてもよい。
【００７８】
本発明の実施形態は、設計変数の効果をさらに導入することができる。上記で説明される例において、これは、ｗ＝２０、２５、３３、５０、および１００Ｈｚの帯域幅を備えるサンプル、ならびにｔｉｍｅ＝ｃｏｎｔｒｏｌｌｅｒ＿ｏｎ＿ｔｉｍｅにおいてフラップを作動するデジタルコントローラを追加することによって行われる。データセット
【００７９】
【数１９】
2016146169.tif 000021
【００８０】
をもたらす。次いで、上記で説明されるプロセスは、このデータセットで反復されて、ｔｉｍｅ＝ｃｏｎｔｒｏｌｌｅｒ＿ｏｎ＿ｔｉｍｅにおいて作動される新しいサロゲートモデルを作成する。そのような実施形態において、４階ＲＭＳ方法は機能しない。しかし、特許文献１０において説明されるように、ＲＭＳ方法を放射基底関数サロゲート方法に置換えることにより。
【００８１】
この修正されたサロゲートモデルの結果は図９において示され、ここでグラフ９８０はα対時間、およびグラフ９８１はｚ対時間を示す。グラフにおいて、より淡く陰影付けされた線９８２および９８４は、完全協調シミュレーション結果を表し、より濃い線９８３および９８５は、コントローラオフの４階ＲＳＭ、および４．８秒および３３ＨｚにおけるコントローラオンのＲＢＦ９８６で近似されたαおよびｚを示す。
【００８２】
上記で説明されるのは、トレーニングデータを使用する本発明の原理によるサロゲートモデルを構築するための方法である。しかし、本発明の実施形態は、それほど限定されてはおらず、構築されるサロゲートモデルは、ユーザによって提供されるデータのような、その他のデータを使用して、シミュレーション／インタラクティブエクスペリエンスに使用されてもよい。そのような実施形態において、翼システムの振る舞いは、異なるシステム帯域幅、異なるシステムコントローラオン時間、および異なる初期条件を伴うサロゲートモデルを使用することによって決定されてもよい。図１０は、異なるデータによる、コントローラオフの４階ＲＳＭ、および３．１秒および４２ＨｚにおけるコントローラオンのＲＢＦで近似される、システム原理の検証を示す。グラフ１０９０は、前述の条件の下の
【００８３】
【数２０】
2016146169.tif 000022
【００８４】
対時間の応答を示す。グラフの部分１０９１は、図１０においてさらに拡大されて、コントローラがオンにされ１０９２、ＲＢＦシミュレーション技法が採用される１０９４場合の応答のさらなる詳細を示す。グラフ１０９０は、サロゲートデータ１０９５が、完全協調シミュレーションデータ１０９３に密接に沿っていることを指示する。
【００８５】
図１０の拡大された領域１０９１で観察されうるように、数値積分エラー（オイラー法の時間ステップ）および時系列間のオフセットに至る可能性のある近似エラーにより、時間の経過に伴って一部の状態エラーが累加する。コントローラが、３．１ｓにおいてオンにされると、完全協調シミュレーションは、フラップ偏向により非常に高いピッチ加速度
【００８６】
【数２１】
2016146169.tif 000023
【００８７】
（−１２０ｒａｄ／ｓ２）への滑らかな移行を生み出す。しかし、サロゲートは、単一の時間ステップでコントローラオフ（０Ｈｚ）から４２Ｈｚへとジャンプし、そのため半径方向加速度にステップを作り出す。しかし、そのポイントからはサロゲートによって予測される加速度は滑らかであるが、デジタルコントローラによる完全協調シミュレーションは、ピッチ加速度に高頻度の離散的修正を有する。しかし、これは実践的なアプリケーションに対しては容認可能なエラーである。そのような実施形態は、リアルタイムのインタラクティブエクスペリエンスにおいて採用する正確なリアルタイム（サロゲート）モデルを提供する。たとえば、サロゲートモデルは、システムが不安定になった場合に、帯域幅または作動時間／条件があるかどうかを決定するために使用されうる。さらにもう１つの実施形態において、サロゲートモデルは、シミュレーションループの大気条件に起因する迎え角の変動のような、時間依存の境界条件を追加することによって、さらに機能強化される。
【００８８】
上記において紹介され、図１１ａおよび図１１ｂに関連してさらに説明されるように、本発明の実施形態は、システムをシミュレートする物理的微分方程式の知識を必要とはしない。必要なのは、データセットが関連性のある状態／設計ベクトルを含むことのみである。そのような関連性のあるベクトルが含まれない場合、方法は、高い交差検証エラーをもたらし、したがってモデルに関連ベクトルデータが不足していることを通知する。そのような方法は、いずれもモデル形式（ベクトルの方程式）が事前に定義されなければならないことを必要とするカルマンフィルタリングおよびパラメータ推定のような、既存の技法より優位にある。実施形態によれば、サロゲートモデルを構築するために、モデル形式は、検証データセットに対する最小交差検証エラーを伴うサロゲート微分方程式モデルとして反復方法を通じて構築され識別される。そのような方法は、図１１ａおよび図１１ｂに関連して本明細書において後段で説明される。
【００８９】
図１１ａおよび図１１ｂは、本発明の実施形態によるサロゲートモデルを構築するための方法１１００（図１１ａの１１００ａおよび図１１ｂの１１００ｂに分割）を示す。方法１１００は、オーサリングを開始し１１０１、３Ｄシステムセンサー位置におけるシステム設計変数ｖｉ、システム内部状態ｐｉ、関心対象のシステム状態ｑｉ、およびシステム境界条件ｂｉをメモリ領域に定義する（１１０２）。次に、時間系列実験が実行され、結果として得られるシステム内部状態結果ｐｉおよび関心対象のシステム状態ｑｉが、ｍｉｎ＜ｖｉ＜ｍａｘの範囲、および時間ｔ＝０からｔ＝ｔまでの初期条件状態ｑｔ＝０、ｐｔ＝０にわたるｖｉの値について決定される（１１０３）。方法１１００は、設計変数および時間ｔ対状態変数のデータ行を連結することによって続行し、ｒｏｗｉ＝ｖｉ，ｔ，ｐｉ，ｑｉ，ｂｉをもたらす（１１０４）。データ行を生成した後、タイムスタンプｔは、すべての関連する状態ｐｉおよび境界条件ｂｉ時間導関数を計算するために使用され、次いで時間はデータセット内でパラメータとして除去されて、結果としてｒｏｗｉ＝ｖｉ、ｐｉ、ｐ’ｉ、ｐ”ｉ、．．．ｑｉ、ｑ’ｉ、ｑ”ｉ、．．．ｂｉ、ｂ’ｉ、ｂ”ｉが得られ、ここでＸ’＝ｄｘ／ｄｔ、ｘ”＝ｄ２ｘ／ｄｔ２である（１１０５）。プロセス１１０５を実行した後、データ行ｒｏｗｉは、トレーニングデータセットおよび検証データセットにランダムに割り当てられる（１１０６）。この後、関心対象の状態ｐｉ、ｐ’ｉ、ｐ”ｉの時間導関数について、ｉは、たとえば、ｉ＝０から３までの間反復される（１１０７）。方法１１００において、ｉは、後段において説明されるテスト１１０９に従って反復される。各反復ごとに、関数Ｆは、多項式またはその他の回帰技法のようなサロゲート技術を使用して構築され、ｉ＝０である場合、ｑｉ＝Ｆ（ｖｉ，ｐｉ，ｐ’ｉ，ｐ”ｉ，．．．ｂｉ，ｂ’ｉ，ｂ”ｉ．．．）、ｉ＝１である場合、ｑ’ｉ＝Ｆ（ｖｉ，ｐｉ，ｐ’ｉ，ｐ”ｉ，．．．ｂｉ，ｂ’ｉ，ｂ”ｉ．．．）、ｉ＝２である場合、ｑ”ｉ＝Ｆ（ｖｉ，ｐｉ，ｐ’ｉ，ｐ”ｉ，．．．ｂｉ，ｂ’ｉ，ｂ”ｉ）の式がもたらされる（１１０８）。各関数Ｆが、サロゲート技術を使用して作成された後（１１０８）、これは、ステップ１１０３からの検証セットを減算した予測されるデータのｑｉ、ｑ’ｉ、ｑ”ｉがεよりも大きいかどうかを決定するためにテストされ１１０９、ここでεは、たとえばモデルの完全検証を表す線６５２と実際のデータポイント６５３との差のような、モデル検証エラーである。ｑｉ、ｑ’ｉ、ｑ”ｉ（予測−検証セット）がεよりも大きい場合、方法１１００は終了する１１１１。しかし、ｑｉ、ｑ’ｉ、ｑ”ｉ（予測−検証セット）がεよりも小さい場合、ｉは増分され１１１０、方法はプロセス１１０７に戻って続行する。
【００９０】
図１２は、方法１１００、または本明細書において説明される任意の方法を、インタラクティブエクスペリエンスにおいて採用するための方法１２２０を示す。方法１２２０は、ランタイムエクスペリエンスを介して開始し１２２１、次いで、図１１ａおよび図１１ｂに関連して説明されるサロゲート微分方程式Ｆのような方程式が、ｔｉｍｅ＝０における変数値ｋでの新しい初期条件ｍで、およびメモリ領域に格納される新しい時間依存境界条件または事象ｂｎ（時間）で、状態ｐｉおよびｑｉについてオフラインで解かれる（１２２２）。プロセス１２２２は、後処理のために、解の時系列をメモリ領域（ｔｉｍｅ、ｖｋ、ｔ、ｐｉ、ｑｉ、ｂｎ）に格納することによって続行する。次に、リアルタイム使用は、デジタルモックアップまたはハードウェアシミュレータの状態ｐｉおよびｑｉを、ｔｉｍｅ＝０における変数ｋの新しい初期条件ｍ、ならびに境界条件または事象値ｂｎにより駆動することによって開始する（１２２３）。方法１２２０は、プロセスの時点１２２４におけるユーザ、コントローラ、またはコンテキスト値のような、新しい事象値ｂｎの導入にさらに対応する。方法１２２０は、新しい事象値ｂｎが受信されるかどうかを監視し１２２５、受信される場合、リアルタイムで、方程式を解き、クロックと同期するｔｉｍｅ＝ｔｉｍｅ＋ｄｔに時間を進める（１２２６）。新しい事象値が１２２５において検出されない場合、ｐｉ、ｑｉ、およびｂｎは、ユーザのリアルタイムエクスペリエンスを作り出すために、時間ｔにおいてデジタルモックアップまたはハードウェアシミュレータの状態を更新するために使用される（１２２７）。
【００９１】
図１３は、本発明の実施形態によるサロゲートモデルを生成するために使用されうるコンピュータベースのシステム１３３０を示す簡略化されたブロック図である。システム１３３０は、バス１３３５を備える。バス１３３５は、システム１３３０のさまざまなコンポーネント間の相互接続としての役割を果たす。バス１３３５に接続されているのは、キーボード、マウス、ディスプレイ、スピーカーなどのようなさまざまな入力および出力デバイスをシステム１３３０に接続するための入出力デバイスインターフェイス１３３８である。中央演算処理装置（ＣＰＵ）１３３２は、バス１３３５に接続され、コンピュータ命令の実行を提供する。メモリ１３３７は、コンピュータ命令を実行するために使用されるデータの揮発性ストレージを提供する。ストレージ１３３６は、オペレーティングシステム（図示せず）のような、ソフトウェア命令の不揮発性ストレージを提供する。システム１３３０はまた、ワイドエリアネットワーク（ＷＡＮ）およびローカルエリアネットワーク（ＬＡＮ）を含む、当技術分野において知られている任意のさまざまなネットワークに接続するためのネットワークインターフェイス１３３１を備える。
【００９２】
モデル定義モジュール１３３３モジュールは、バス１３３５にさらに接続される。モデル定義モジュール１３３３は、現実世界のシステムを表すモデルを定義するように構成され、定義されるモデルは振る舞いを表すパラメトリック状態ベクトルと、設計変数ベクトルとを含む。
【００９３】
実験モジュール１３３４は、バス１３３５を介して、モデル定義モジュール１３３３に動作可能に結合され、パラメトリック状態ベクトルの経時的な応答を決定するための第１の実験、および設計変数ベクトルに対する第２の実験を実行するように構成される。システム１３３０のそのような実施形態において、第１および第２の実験は、パラメトリック状態ベクトルおよび設計変数ベクトルのデータセットを時間の関数として生成する。システム１３３０は、実験モジュール１３３４からデータセットを受信して、パラメトリック状態ベクトルの１または複数の導関数によりデータセットを修正するように構成されるデータセットモジュール１３４１をさらに含む。バス１３３５に接続されるのは、パラメトリック状態ベクトルの高次導関数をデータセットモジュール１３４１によって決定されたデータセット内のそれに関して近似するサロゲート微分方程式のセットを構築するように構成されるサロゲートモジュール１３３９である。サロゲートモジュール１３３９によるこの構築するステップは、サロゲート微分方程式のセットをサロゲートモデルとして、メモリ１３３７またはストレージデバイス１３３６に動作可能に格納するステップを含む。システム１３３０に動作可能に結合されるのは、モデルとのユーザインタラクションに応答してシミュレートされる振る舞いを促進する方法で、メモリ１３３７またはストレージ１３３６からサロゲートモデルを応答可能に提供するように構成されるインタラクションモジュール１３４０である。サロゲートモデルを提供するこのステップは、サロゲートモデルの状態ベクトルを時間の関数として解く時系列を使用することを含む設計変数ベクトルの所与の値に対するサロゲート微分方程式を解く。
【００９４】
実施形態において、促進されたシミュレートされる振る舞いは、通常のシミュレーション技法と比較して、促進されるものと見なされうる。さらに、促進されたシミュレートされる振る舞いは、たとえば、シミュレーションが翼フラッタの効果を、現実世界において生じるであろうそのようなフラッタよりも速い速度で示しうるなど、現実世界の事象と比較して、促進または減速（decelerated）されてもよい。
【００９５】
本明細書において説明される例示的な実施形態は、多くの異なる方法で実施されてもよいことが理解されるべきである。場合によっては、本明細書において説明されるさまざまな方法およびマシンは、各々、コンピュータシステム１３３０のような、物理、仮想、もしくはハイブリッド汎用コンピュータ、または図１４に関して本明細書において後段で説明される、コンピュータ環境１４００のようなコンピュータネットワーク環境によって実施されてもよい。コンピュータシステム１３３０は、たとえば、ソフトウェア命令をＣＰＵ１３３２により実行するためにメモリ１３３７または不揮発性ストレージ１３３６のいずれかにロードすることによって、本明細書において説明される方法を実行するマシンに変えられてもよい。さらに、モデル定義モジュール１３３３、実験モジュール１３３４、データセットモジュール１３４１、サロゲートモジュール１３３９、およびインタラクションモジュール１３４０は、例示的な実施形態において、別個のモジュールとして示されているが、これらのモジュールはさまざまな構成を使用して実施されてもよい。当業者は、システム１３３０およびそのさまざまなコンポーネントが、本明細書において説明される本発明の任意の実施形態を実行するように構成されうることを、さらに理解されたい。
【００９６】
図１４は、本発明の実施形態が実施されうるコンピュータネットワーク環境１４００を示す。コンピュータネットワーク環境１４００において、サーバ１４０１は、通信ネットワーク１４０２を通じてクライアント１４０３ａ〜１４０３ｎにリンクされる。環境１４００は、クライアント１４０３ａ〜１４０３ｎが、単独で、またはサーバ１４０１との組合せで、本明細書において説明される方法のいずれかを実行できるようにするために使用されてもよい。
【００９７】
その実施形態または態様は、ハードウェア、ファームウェア、またはソフトウェアの形態において実施されてもよい。ソフトウェアにおいて実施される場合、ソフトウェアは、プロセッサがソフトウェアまたはその命令のサブセットをロードできるように構成される任意の非一時的コンピュータ可読媒体上に格納されてもよい。次いで、プロセッサは、命令を実行し、本明細書において説明される方法で、動作するか、または装置に動作させるように構成される。
【００９８】
さらに、ファームウェア、ソフトウェア、ルーチン、または命令は、データプロセッサの特定のアクションおよび／または機能を実行するものとして本明細書において説明されてもよい。しかし、本明細書に含まれるそのような説明は、便宜上のものにすぎず、実際にそのようなアクションが、ファームウェア、ソフトウェア、ルーチン、命令などを実行するコンピューティングデバイス、プロセッサ、コントローラ、またはその他のデバイスによりもたらされることを理解されたい。
【００９９】
流れ図、ブロック図、およびネットワーク図は、異なって配置されるかまたは異なって表される、より多くまたはより少ない要素を含むことができることを理解されたい。しかし、特定の実施態様が、ブロックおよびネットワーク図を指示することができ、実施形態の実行を説明するブロックおよびネットワークの数が特定の方法で実施されてもよいことをさらに理解されたい。
【０１００】
したがって、さらなる実施形態はまた、さまざまなコンピュータアーキテクチャ、物理、仮想、クラウドコンピュータ、および／またはその一部の組合せにおいて実施されてもよく、したがって、本明細書において説明されるデータプロセッサは、例示のみを目的とし、実施形態を限定するものではないことが意図される。
【０１０１】
本発明は、その例示的な実施形態を参照して詳細に示され、説明されたが、添付の特許請求の範囲によって包含される本発明の範囲を逸脱することなく、形態および詳細のさまざまな変更がそこに行われてもよいことが、当業者には理解されるであろう。
【図１】
2016146169.tif 000024
【図５】
2016146169.tif 000025
【図９】
2016146169.tif 000026
【図１１ａ】
2016146169.tif 000027
【図１１ｂ】
2016146169.tif 000028
【図１２】
2016146169.tif 000029
【図１３】
2016146169.tif 000030
【図１４】
2016146169.tif 000031
【図２】
2016146169.tif 000032
【図３】
2016146169.tif 000033
【図４】
2016146169.tif 000034
【図６】
2016146169.tif 000035
【図７】
2016146169.tif 000036
【図８】
2016146169.tif 000037
【図１０】
2016146169.tif 000038
【外国語明細書】
2016146169000001.pdf