JP2011252905A - パラメータ同定のため計算された曲線を目標曲線に一致させる方法 - Google Patents
パラメータ同定のため計算された曲線を目標曲線に一致させる方法 Download PDFInfo
- Publication number
- JP2011252905A JP2011252905A JP2011118253A JP2011118253A JP2011252905A JP 2011252905 A JP2011252905 A JP 2011252905A JP 2011118253 A JP2011118253 A JP 2011118253A JP 2011118253 A JP2011118253 A JP 2011118253A JP 2011252905 A JP2011252905 A JP 2011252905A
- Authority
- JP
- Japan
- Prior art keywords
- curve
- points
- target
- computer
- calculated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 238000013507 mapping Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims description 11
- 230000006399 behavior Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 description 15
- 238000004590 computer program Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 11
- 238000013461 design Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Complex Calculations (AREA)
Abstract
【解決手段】パラメータ同定の最適化は、2つの曲線302,304の間の相違が最小限となるようにシミュレーションモデルのパラメータ入力を調整することによって達成される。一致させる2つの曲線302,304上の点312,316をペアにするので、ヒステリシス曲線を含む2つのあらゆる開いた曲線の一致に対応できる。元の座標において完全に離れている曲線を、計算上の不安定性が問題になることなく、パラメータ同定に対する共通の座標系に合成することができる。部分的に一致させる手法(410g)が、2つの曲線302,304のうちより短い曲線を定義する点312をより長い曲線上の1セットのマップ点316にマッピングするよう用いられる。長い曲線304の第1点から1以上のオフセット320が、複数回の試行401gに用いられて、最もよい一致を見つける。
【選択図】図4A
Description
座標)によって定義される。「曲線G」304上には、第3セットの点すなわちマップ点316(円で表示し、(ξ”、η”)はそれぞれの座標)が配置される。マップ点316
のそれぞれは、第1点312の各点からマップされる(割り当てられる)。両方の曲線は、短い曲線の最大制限範囲(maximum bounding ranges)に規格化される。第1および第3セット間の各対応するペアの点は、距離di310で表示される。
1)目標曲線のm個の点座標iを制限ボックス(bounding box)すなわち範囲が最も小さくなるように規格化し、図3Aに示す「曲線F」302を生成する。
また、個々のセグメント長さ322を計算する。
4)曲線全長さSに対する各セグメント長さを比で表す(スケールする)。
6)S>Tの場合、2つ曲線を入れ換えて、短い曲線を「曲線F」302とする。
7)「曲線G」404上の全長Sの曲線部分のスタート点として、オフセット320を定義する。オフセットを、短い曲線を長い曲線に沿って「スライドさせる」ために、多くの異なる位置(P)上で変化させる。オフセットの範囲は、ゼロと、2つの曲線の長さ差(つまり、|T−S|)との間とする。部分的に一致させる試行の回数がPであると仮定すると、各オフセットインクリメントの大きさは|T−S|/Pである。
8)オフセットを設定し、長い曲線の新たな一部分を生成するとともに、短い曲線を定義する各点を長い曲線にマップすることによって点座標ペアを生成する。オフセットは、最初、第1回目の部分的に一致させる試行ではゼロからスタートし、以降の各試行においてはオフセットインクリメントずつインクリメントされる。その結果、最後の試行におけるオフセットは|T−S|すなわち2つの曲線の長さの差と等しくなる。短い曲線上のセグメントiに対応する長い曲線上の典型的な曲線セグメントkの長さは、δUk=δSiである(図3A参照)。これにより、第3セットの点すなわち1セットのマップ点が生成される。長い曲線のマップされた部分の長さは短い曲線の長さと等しい。
9)2つの曲線間の相違(不整合エラー)を計算する。これは、個々のセグメントエラーを表わすボリューム成分vi(面積)を合計することによって行われる。まず、点ペア間の距離を計算する。
11)相違
を見つける。
これは、2つの曲線間の最もよい一致である(つまりそれがすべての部分的に一致させる試行の中の最小の相違である)。
502 バス
504 プロセッサ
506 アプリケーションモジュール
508 メインメモリ
510 二次メモリ
512 ハードディスクドライブ
514 リムーバブルストレージドライブ
518 リムーバブルストレージユニット
520 インタフェース
522 リムーバブルストレージユニット
524 通信インタフェース
530 I/Oインタフェース
Claims (14)
- コンピュータ支援工学シミュレーションにおいて計算曲線(304)を材料挙動(100)を表す目標曲線(302)に一致させるコンピュータシステム(500)において実行される方法(400)であって、
コンピュータシステム(500)にインストールされたアプリケーションモジュール(506)によって目標曲線(302)の定義を受け取るステップ(402)と、
複数の制御パラメータ(404)を介して目標曲線(302)を生成するよう構成されたコンピュータ支援工学解析モデル(404)を用いて計算曲線(304)を取得するステップ(408)と、
前記目標曲線および計算曲線(302、304)のうちより短い曲線を短曲線(302)としてそしてより長い曲線を長曲線(304)として指定するステップ(410e)であって、前記短曲線(302)が第1セットの点(312)によって定義されるとともに前記長曲線(304)は第2セットの点(314)によって定義されるステップと、
1回以上の部分的に一致させる試行(410g)を含むマッピング法を用いて前記第1セットの点(312)を前記長い曲線(304)上にマップすることにより1セットのマップ点(316)を生成するステップ(410)であって、すべての前記部分的に一致させる試行の中で最もよい一致を示した前記部分的に一致させる試行のうちの1回が前記1セットのマップ点(316、410h)として指定され、前記最もよい一致が前記マップ点(316)のそれぞれと前記部分的に一致させる試行のそれぞれに対する前記第1セットの点(312)のうちの対応する1つとの間で計算される相違測定を用いて判定されるステップと、
前記計算曲線(302)と前記目標曲線(304)との間の前記相違測定が許容範囲内となるまで前記制御パラメータを調整することによって前記計算曲線(304)を反復的に更新するステップ(408、410、412、414)と、
前記制御パラメータを前記コンピュータ支援工学解析モデル(404)とともにユーザの指示に応じてコンピュータシステム(500)に連結された記憶装置(512)のファイルに保存するステップと、
を備える方法。 - 請求項1に記載の方法であって、前記部分的に一致させる試行(410g)のそれぞれは、ゼロから短曲線および長曲線(302、304)の長さの差と等しい最大値の範囲にわたるオフセット(320)に対応しており、該オフセット(320)は前記第2セットの点(314)の第1点から測定される方法。
- 請求項2に記載の方法であって、部分的に一致させる試行の前記それぞれは、隣接した2つのマップ点(316)によって形成される各セグメント(326)の長さが前記短い曲線の対応するセグメントのセグメント長さと等しいことを保証することによって暫定的セットのマップ点(316)を生成することをさらに備えており、前記暫定的セットのマップ点(316)の第1点は対応するオフセット(320)からスタートしている方法。
- 請求項3に記載の方法であって、前記短曲線(322)の前記セグメントは、前記短曲線を定義する2つの隣接点(312)間の直線によって定義される方法。
- 請求項3に記載の方法であって、2つの隣接するマップ点(316)によって形成される各セグメント(326)の長さは、その2つの隣接するマップ点(316)の間に配置されたセグメント全体あるいはセグメントの一部分の長さの和である方法。
- 請求項2に記載の方法であって、前記1回以上の部分的に一致させる試行の数は、前記オフセット(320)の最大値と前記長い曲線の長さとの比に基づいている方法。
- 請求項1に記載の方法であって、前記目標および計算曲線(302、304)のうちより短い曲線を短い曲線としてそしてより長い曲線を長い曲線として指定する前記ステップ(410e)は、さらに、
前記目標曲線(302)の最小制限範囲を用いて前記第1および第2セットの点(312、314)の座標を規格化するステップ(410b)と、
前記目標曲線(302)の全長と前記計算曲線(304)の全長を計算するステップ(410c)と、
を備えている方法。 - 請求項1に記載の方法であって、前記最もよい一致の相違測定が前記すべての部分的に一致させる試行(410g)の中で最も低い方法。
- 請求項8に記載の方法であって、前記短曲線と前記長曲線(302、304)との間の相違測定は、前記短曲線(302)上の各セグメントと前記長曲線(304)上のその対応するセグメントとの間の面積の和である方法。
- 請求項1に記載の方法であって、前記目標曲線(302)は、材料試料試験において得られる材料(100)の歪み−応力関係曲線である方法。
- コンピュータ支援工学シミュレーションにおいて計算曲線(304)を材料挙動(100)を表す目標曲線(302)に一致させるシステムであって、
コンピュータ可読コードを記憶するメモリ(508)と、
前記メモリ(508)に連結される少なくとも1つのプロセッサ(504)と、
を備えており、
前記少なくとも1つのプロセッサ(504)が前記メモリ(508)内の前記コンピュータ可読コードを実行して、これにより、1つ以上のアプリケーションモジュール(506)に、
システムにインストールされた前記1つ以上のアプリケーションモジュール(506)によって目標曲線(302)の定義を受け取るオペレーション(400、402)と、
複数の制御パラメータ(404)を介して計算曲線(302)を生成するよう構成されたコンピュータ支援工学解析モデル(404)を用いて計算曲線(304)を取得するオペレーション(408)と、
目標および計算曲線(302、304)のうちより短い曲線を短曲線(302)としてそしてより長い曲線を長曲線(304)として指定するオペレーション(410e)であって、短曲線(302)が第1セットの点(312)によって定義されるとともに長曲線(304)は第2セットの点(314)によって定義されるオペレーションと、
1回以上の部分的に一致させる試行(410g)を含むマッピング法を用いて前記第1セットの点(312)を前記長曲線(304)上にマップすることにより1セットのマップ点(316)を生成するオペレーション(410)であって、すべての前記部分的に一致させる試行の中で最もよい一致を示した前記部分的に一致させる試行のうちの1回が前記1セットのマップ点(316、410h)として指定され、前記最もよい一致が前記マップ点(316)のそれぞれと前記部分的に一致させる試行のそれぞれに対する前記第1セットの点(312)のうちの対応する1つとの間で計算される相違測定を用いて判定されるオペレーションと、
前記計算曲線(304)と前記目標曲線(302)との間の前記相違が許容範囲内となるまで前記制御パラメータを調整することによって前記計算曲線(304)を反復的に更新するオペレーション(408、410、412、414)と、
前記制御パラメータを前記コンピュータ支援工学解析モデル(404)とともにユーザの指示に応じてシステムに連結された記憶装置(512)のファイルに保存するオペレーションと、
を実行させるシステム。 - 請求項11に記載のシステムであって、前記部分的に一致させる試行のそれぞれは、ゼロから短曲線および長曲線(302、304)の長さの差と等しい最大値の範囲にわたるオフセット(320)に対応しており、該オフセット(320)は前記第2セットの点(314)の第1点から測定される方法。
- 請求項12に記載のシステムであって、部分的に一致させる試行(410g)の前記それぞれは、隣接した2つのマップ2点(316)によって形成される各セグメントの長さが前記短曲線(302)の対応するセグメントのセグメント長さと等しいことを保証することによって暫定的セットのマップ点(316)を生成することをさらに備えており、前記暫定的セットのマップ点(316)の第1点は対応するオフセット(320)からスタートしている方法。
- 請求項1、2あるいは3に記載の方法(400)を実行することによって、コンピュータ支援工学シミュレーションにおいて計算曲線(304)を材料挙動(100)を表す目標曲線(302)に一致させるコンピュータシステム(500)を制御するコンピュータが実行可能な命令を有するコンピュータ可読媒体。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/792,700 | 2010-06-02 | ||
US12/792,700 US8521484B2 (en) | 2010-06-02 | 2010-06-02 | Curve matching for parameter identification |
US13/048,322 | 2011-03-15 | ||
US13/048,322 US8515714B2 (en) | 2010-06-02 | 2011-03-15 | Curve matching for parameter identification |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2011252905A true JP2011252905A (ja) | 2011-12-15 |
JP2011252905A5 JP2011252905A5 (ja) | 2014-05-29 |
JP5823170B2 JP5823170B2 (ja) | 2015-11-25 |
Family
ID=44654626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011118253A Active JP5823170B2 (ja) | 2010-06-02 | 2011-05-26 | パラメータ同定のため計算された曲線を目標曲線に一致させる方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8515714B2 (ja) |
EP (1) | EP2393027A3 (ja) |
JP (1) | JP5823170B2 (ja) |
CN (1) | CN102314533B (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019507889A (ja) * | 2015-12-10 | 2019-03-22 | メカニカル テスティング サービシズ, エルエルシーMechanical Testing Services, Llc | インテリジェント自動負荷制御システムおよび方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8521484B2 (en) * | 2010-06-02 | 2013-08-27 | Livermore Software Technology Corp. | Curve matching for parameter identification |
US9576088B2 (en) | 2013-01-23 | 2017-02-21 | Toyota Motor Engineering & Manufacturing North America, Inc. | Methods for orienting material physical properties using constraint transformation and isoparametric shape functions |
US9576084B2 (en) * | 2013-08-27 | 2017-02-21 | Halliburton Energy Services, Inc. | Generating a smooth grid for simulating fluid flow in a well system environment |
CN106826813B (zh) * | 2016-12-12 | 2019-08-23 | 杭州宇芯机器人科技有限公司 | 一种六足机器人稳定运动控制方法 |
CN108460113B (zh) * | 2018-02-09 | 2021-09-24 | 水利部交通运输部国家能源局南京水利科学研究院 | 可用于比较水文领域中理论计算曲线过拟合程度的方法 |
CN113342881B (zh) * | 2021-05-21 | 2023-07-07 | 中广核工程有限公司 | 核电厂仪控系统测试曲线构建方法、装置、设备、介质 |
TWI802363B (zh) * | 2022-04-08 | 2023-05-11 | 閎康科技股份有限公司 | 曲線對齊方法及曲線對齊裝置 |
CN116227351A (zh) * | 2023-03-06 | 2023-06-06 | 陕西海鹰汽车部件有限公司 | 基于主+副双运输管道的新型高位进气管设计方法、装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09152392A (ja) * | 1995-11-30 | 1997-06-10 | Sony Corp | 非線形構造解析方法 |
JP2002296164A (ja) * | 2001-03-30 | 2002-10-09 | Sumitomo Chem Co Ltd | 樹脂成形品の衝撃解析方法及び設計方法 |
US20030025712A1 (en) * | 2001-08-06 | 2003-02-06 | Corr William E. | Interpolation error minimization for data reduction |
US20030093235A1 (en) * | 2001-11-15 | 2003-05-15 | Ibm Corporation | Method and apparatus for detecting and correcting inaccuracies in curve-fitted models |
JP2007200290A (ja) * | 2005-12-27 | 2007-08-09 | Semiconductor Energy Lab Co Ltd | パラメータ抽出方法及び回路動作検証方法並びに当該パラメータ抽出方法を実行させるためのプログラムを具備する記憶媒体 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5694535A (en) * | 1995-03-24 | 1997-12-02 | Novell, Inc. | Direct interactive, constant-time curve apparatus and method |
CN101436254B (zh) * | 2007-11-14 | 2013-07-24 | 佳能株式会社 | 图像处理方法和图像处理设备 |
CN101266645B (zh) * | 2008-01-24 | 2011-01-19 | 电子科技大学中山学院 | 一种基于多分辨率分析的虹膜定位方法 |
-
2011
- 2011-03-15 US US13/048,322 patent/US8515714B2/en active Active
- 2011-05-26 EP EP11167629.2A patent/EP2393027A3/en not_active Withdrawn
- 2011-05-26 JP JP2011118253A patent/JP5823170B2/ja active Active
- 2011-05-27 CN CN201110140142.4A patent/CN102314533B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09152392A (ja) * | 1995-11-30 | 1997-06-10 | Sony Corp | 非線形構造解析方法 |
JP2002296164A (ja) * | 2001-03-30 | 2002-10-09 | Sumitomo Chem Co Ltd | 樹脂成形品の衝撃解析方法及び設計方法 |
US20030025712A1 (en) * | 2001-08-06 | 2003-02-06 | Corr William E. | Interpolation error minimization for data reduction |
US20030093235A1 (en) * | 2001-11-15 | 2003-05-15 | Ibm Corporation | Method and apparatus for detecting and correcting inaccuracies in curve-fitted models |
JP2007200290A (ja) * | 2005-12-27 | 2007-08-09 | Semiconductor Energy Lab Co Ltd | パラメータ抽出方法及び回路動作検証方法並びに当該パラメータ抽出方法を実行させるためのプログラムを具備する記憶媒体 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019507889A (ja) * | 2015-12-10 | 2019-03-22 | メカニカル テスティング サービシズ, エルエルシーMechanical Testing Services, Llc | インテリジェント自動負荷制御システムおよび方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5823170B2 (ja) | 2015-11-25 |
US8515714B2 (en) | 2013-08-20 |
EP2393027A2 (en) | 2011-12-07 |
CN102314533B (zh) | 2015-04-15 |
US20110301920A1 (en) | 2011-12-08 |
EP2393027A3 (en) | 2017-06-28 |
CN102314533A (zh) | 2012-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5823170B2 (ja) | パラメータ同定のため計算された曲線を目標曲線に一致させる方法 | |
JP5411057B2 (ja) | 有限要素解析法におけるビーム−面接触をシミュレートする改善方法およびシステム | |
CN107016155B (zh) | 非线性pde和线性求解器的收敛估计 | |
US11281824B2 (en) | Authoring loading and boundary conditions for simulation scenarios | |
Witowski et al. | An effective curve matching metric for parameter identification using partial mapping | |
JP6737630B2 (ja) | 金属くびれ破損の発生が予想される構造の時間進行数値シミュレーションを行う方法およびシステム | |
US20170255724A1 (en) | Enhanced Global Design Variables Used In Structural Topology Optimization Of A Product In An Impact Event | |
US9977993B2 (en) | System and method for constructing a statistical shape model | |
WO2020031802A1 (ja) | 学習方法、学習装置、モデル生成方法及びプログラム | |
JP6560577B2 (ja) | 適応分類を用いた多目的設計最適化 | |
JP2011159288A (ja) | 有限要素法において8ノード六面体エレメントの剪断ロッキングを低減する方法 | |
JP2019220156A (ja) | 物理システムのハイブリッドモデル生成 | |
US8521484B2 (en) | Curve matching for parameter identification | |
WO2021117162A1 (ja) | 検査装置、検査方法及び検査プログラム | |
Jensen | Statistical inference for partially observed diffusion processes | |
Stankiewicz et al. | Genetic algorithm-based calibration of reduced order Galerkin models | |
Liatsikouras et al. | Aerodynamic shape optimization under flow uncertainties using non-intrusive polynomial chaos and evolutionary algorithms | |
US11347911B2 (en) | Systems and methods of determining a numerical material model that optimally emulates physical material test results | |
JP2010181082A (ja) | モデル関数処理装置および方法 | |
Pagani et al. | A sampling strategy based on B-wavelets decomposition | |
US20150025859A1 (en) | Computer-implemented method for performing simulation | |
Gralla et al. | Inverting prediction models in micro production for process design | |
Bellucci et al. | A Taylor series approach to the robust parameter design of computer simulations using kriging and radial basis function neural networks | |
Mosegaard | Parameter optimisation for the behaviour of elastic models over time | |
US20230142773A1 (en) | Method and system for real-time simulations using convergence stopping criterion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140414 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140414 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150210 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150403 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150929 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151007 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5823170 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |