(57)【要約】【課題】解析モデル化最適制御に必要な時間的コストを削減する。【解決手段】代表形状モデル１０９の代表形状特徴パラメータを取得する手段（１０２）と、代表形状モデル１０９の代表形状最適メッシュ制御パラメータを決定する手段（１０３）と、代表形状特徴パラメータと代表形状最適メッシュ制御パラメータとの１以上の組み合わせ（１１２）から応答曲面を構築する手段（１０４）と、前記応答曲面の応答曲面構築パラメータ１１３を格納する応答曲面データベース１０５と、メッシュ生成対象形状モデル１１４のメッシュ生成対象形状特徴パラメータ１１５と応答曲面構築パラメータ１１３に基づいて対象形状メッシュ制御パラメータ１１６を決定する手段（１０７）を備える。これにより、解析モデル化最適制御を解析モデルごとに毎回行う必要が無くなるので、形状に応じてメッシュ生成を高速に行うことを可能とした。【選択図】図１
2011243016.tif 000002
【特許請求の範囲】
【請求項１】
代表形状のモデルから取得され、前記代表形状の形状に関する特徴を示す第１の形状特徴パラメータと、
前記代表形状に用いるメッシュに関する特徴を示す第１のメッシュ制御パラメータを最適化した最適メッシュ制御パラメータと、
前記第１の形状特徴パラメータおよび前記最適メッシュ制御パラメータに基づいて生成された応答曲面を構築する応答曲面構築パラメータと、を記憶する記憶部と、
メッシュ生成対象形状のモデルから、前記メッシュ生成対象形状の形状に関する特徴を示す第２の形状特徴パラメータを取得する制御と、
前記応答曲面構築パラメータから構築した前記応答曲面に、前記第２の形状特徴パラメータを入力することで、前記メッシュ生成対象形状に用いるメッシュに関する特徴を示す第２のメッシュ制御パラメータを決定する制御と、
前記第２のメッシュ制御パラメータに基づいて、前記メッシュ生成対象形状に用いるメッシュを生成する制御と、を実行する制御部と、を備える
ことを特徴とする解析用メッシュ生成装置。
【請求項２】
前記最適メッシュ制御パラメータは、メッシュの粗密を定めるパラメータである
ことを特徴とする請求項１に記載の解析用メッシュ生成装置。
【請求項３】
前記制御部は、
前記第１のメッシュ制御パラメータを変化させる制御と、
前記代表形状について、前記変化させた第１のメッシュ制御パラメータに基づいて生成したメッシュを用いて解析したときの解析結果を取得する制御と、
前記記憶部に記憶され、前記代表形状について所定の精度を示す参照結果と、前記取得した解析結果との差異を求める制御と、
前記差異を最小化する第１のメッシュ制御パラメータを前記最適メッシュ制御パラメータとする制御と、を実行する
ことを特徴とする請求項１に記載の解析用メッシュ生成装置。
【請求項４】
前記制御部は、
前記第１のメッシュ制御パラメータが複数ある場合、前記第１のメッシュ制御パラメータの組み合わせを入力すると、前記差異を出力する応答曲面を構築する制御と、
前記構築した応答曲面を用いて前記最適メッシュ制御パラメータを求める制御と、を実行する
ことを特徴とする請求項３に記載の解析用メッシュ生成装置。
【請求項５】
前記制御部は、
前記第１のメッシュ制御パラメータを変化させる制御と、
入力部からの入力により、前記代表形状について目的関数を設定する制御と、
前記設定された目標関数の値を最小化する第１のメッシュ制御パラメータを前記最適メッシュ制御パラメータとする制御と、を実行する
ことを特徴とする請求項１に記載の解析用メッシュ生成装置。
【請求項６】
前記制御部は、
前記第１のメッシュ制御パラメータが複数ある場合、前記第１のメッシュ制御パラメータの組み合わせを入力すると、前記目標関数の値を出力する応答曲面を構築する制御と、
前記構築した応答曲面を用いて前記最適メッシュ制御パラメータを求める制御と、を実行する
ことを特徴とする請求項５に記載の解析用メッシュ生成装置。
【請求項７】
代表メッシュから取得され、前記代表メッシュの形状に関する特徴を示す第１の形状特徴パラメータと、
前記代表メッシュに関する特徴を示す第１のメッシュ制御パラメータを最適化した最適メッシュ制御パラメータと、
前記第１の形状特徴パラメータおよび前記最適メッシュ制御パラメータに基づいて生成された応答曲面を構築する応答曲面構築パラメータと、を記憶する記憶部と、
変形対象メッシュに対し、前記変形対象メッシュの形状の変形に関する特徴を示す第２の形状特徴パラメータを取得する制御と、
前記応答曲面構築パラメータから構築した前記応答曲面に、前記第２の形状特徴パラメータを入力することで、前記変形対象メッシュに関する特徴を示す第２のメッシュ制御パラメータを決定する制御と、
前記第２のメッシュ制御パラメータに基づいて、前記変形対象メッシュを変形する制御と、を実行する制御部と、を備える
ことを特徴とする解析用メッシュ生成装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、有限要素法等の数値解析の技術に関する。
【背景技術】
【０００２】
特許文献１には、実験結果や厳密な解析結果とモデル化された解析結果との間で大きな相違を生じさせず、かつ解析メッシュ数を可能な限り低減し解析時間を短縮するように、解析対象を適切にモデル化する解析モデル化最適制御装置を提供する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００１−１３４７８３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記特許文献１に記載の従来技術は、実験結果や詳細な解析モデルを用いた解析結果と、何らかのモデル化が施された解析モデルを用いた解析結果との間での大きな差異が生じないよう、三次元形状モデルを構成する基本形状となるブロック、穴等に対し、最も解析結果の差異が小さくなるようなモデル化手法の組み合わせを、遺伝的アルゴリズムを用いて決定し、解析モデル作成支援を行う。しかしながら、上記従来技術での解析モデル化最適制御は、解析モデルごとに毎回行わなければならず、時間的コストを要する。
【０００５】
そこで本発明の目的は、解析モデル化最適制御に必要な時間的コストを削減することである。
【課題を解決するための手段】
【０００６】
前記目的を達成するため、本発明は、代表形状から形状特徴パラメータを取得する手段と、前記代表形状についてメッシュ制御パラメータを変化させてメッシュ生成を行い、これを用いて行った解析結果と前記代表形状の実験結果または詳細解析結果との差異、もしくは、装置利用者が定義する全体解析メッシュ数や解析時間といった目的関数を最適化した最適メッシュ制御パラメータを取得する手段と、上記代表形状の形状特徴パラメータと上記最適メッシュ制御パラメータの複数の組み合わせから、形状特徴パラメータとメッシュ制御パラメータに関する応答曲面を構築する手段と、上記構築された応答曲面の構築パラメータを格納する応答曲面データベースと、入力されたメッシュ生成対象となる形状から形状特徴パラメータを取得する手段と、上記メッシュ生成対象となる形状の形状特徴パラメータを入力し、上記応答曲面データベースに格納されている応答曲面の構築パラメータから構築される応答曲面上でのメッシュ制御パラメータを決定する手段と、上記応答曲面上で決定されたメッシュ制御パラメータにより、上記入力されたメッシュ生成対象となる形状に対しメッシュを生成する手段、を備えることにより、新たにメッシュ生成対象となる形状に対して、解析精度の良いメッシュ制御パラメータを高速に決定することができる。
詳細は、後記する。
【発明の効果】
【０００７】
本発明によれば、解析モデル化最適制御に必要な時間的コストを削減することができる。
【図面の簡単な説明】
【０００８】
【図１】本実施形態の解析用メッシュ生成装置の構成例である。
【図２】代表形状特徴パラメータ取得部における画面構成例である。
【図３】代表形状メッシュ制御パラメータ最適化部における誤差最小化対象設定の画面構成例である。
【図４】代表形状メッシュ制御パラメータ最適化の手順を示すフローチャートである。
【図５】代表形状メッシュ制御パラメータ最適化部における最適化問題定義の画面構成例である。
【図６】複数代表形状に対する形状特徴パラメータ、最適メッシュ制御パラメータの取得例である。
【図７】応答曲面構築パラメータの具体例である。
【図８】メッシュ生成対象形状の形状モデルと形状特徴パラメータの具体例である。
【図９】メッシュ制御パラメータの具体例である。
【図１０】対象形状メッシュの具体例である。
【図１１】本装置をメッシュモーフィングによるメッシュ生成に適用した場合の具体例である。
【発明を実施するための形態】
【０００９】
次に、本発明を実施するための形態（以下、「実施形態」という。）について、適宜図面を参照しながら説明する。
【００１０】
図１は、本実施形態の解析用メッシュ生成装置の構成例である。本装置は、入出力装置１０１をハードウェア構成として備えている。また、本装置は、代表形状特徴パラメータ取得部１０２と、代表形状メッシュ制御パラメータ最適化部１０３と、応答曲面構築部１０４と、応答曲面データベース１０５と、メッシュ生成対象形状特徴パラメータ取得部１０６と、メッシュ制御パラメータ決定部１０７と、メッシュ生成部１０８とをソフトウェア構成として備えている。以下、これらの構成の概要を説明する。
【００１１】
入出力装置１０１は、装置利用者が、データを入力、表示するためのキーボード、マウス、ディスプレイ等の入力部、表示部を備えるとともに、制御部、記憶部を備える。
前記制御部は、例えば、ＣＰＵ（Central Processing Unit）であり、メッシュ生成処理等の各種情報処理を実行する。
前記記憶部は、例えば、ＨＤＤ（Hard Disk Drive）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）等であり、前記ソフトウェア構成による機能を実現し、前記ＣＰＵが情報処理を実行するために読み出すプログラムを記憶する。また、記憶部は、代表形状モデル１０９と、解析条件１１０と、参照結果１１１と、代表形状特徴パラメータと代表形状最適メッシュ制御パラメータの組み合わせ１１２と、応答曲面構築パラメータ１１３と、メッシュ生成対象形状モデル１１４と、メッシュ生成対象形状特徴パラメータ１１５と、対象形状メッシュ制御パラメータ１１６と、対象形状メッシュ１１７を記憶する。
【００１２】
代表形状特徴パラメータ取得部１０２は、代表形状をＣＡＤ（Computer Aided Design）ソフトウェア等によりモデル化した代表形状モデル１０９を入力し、その代表形状の形状特徴パラメータを代表形状特徴パラメータ（第１の形状特徴パラメータ）として取得する。
【００１３】
代表形状メッシュ制御パラメータ最適化部１０３は、代表形状モデル１０９および解析条件１１０を入力するとともに、事前に取得した代表形状の実験結果または詳細解析結果を参照結果１１１として入力し、代表形状に用いるメッシュに関する特徴を示すメッシュ制御パラメータ（第１のメッシュ制御パラメータ）を変化させることで、上記解析条件１１０を満たす代表形状の解析結果と、参照結果１１１との差異が許容値以下となる最適なメッシュ制御パラメータを探索する。探索の結果、見つけたメッシュ制御パラメータを代表形状最適メッシュ制御パラメータ（最適メッシュ制御パラメータ）と称する。また、装置利用者が設定する目的関数を最適化する。
【００１４】
応答曲面構築部１０４は、上記代表形状特徴パラメータと上記代表形状最適メッシュ制御パラメータの１以上の組み合わせから、形状特徴パラメータとメッシュ制御パラメータに関する応答曲面を構築する。「応答曲面」とは、サンプリング点のデータを元に作成した近似関数をいい、近似関数上で最適化を行う方法は応答曲面法と呼ばれる。本実施形態では、前記サンプリング点のデータとは、代表形状特徴パラメータおよび代表形状最適メッシュ制御パラメータである。
【００１５】
応答曲面データベース１０５は、上記応答曲面構築部１０４で得られた応答曲面を規定する応答曲面構築パラメータを格納する。
【００１６】
メッシュ生成対象形状特徴パラメータ取得部１０６は、装置利用者が解析したいと希望するメッシュ生成対象形状を、ＣＡＤ（Computer Aided Design）ソフトウェア等によりモデル化したメッシュ生成対象形状モデル１１４を入力し、そのメッシュ生成対象形状の形状特徴パラメータをメッシュ生成対象特徴パラメータ１１５（第２の形状特徴パラメータ）として取得する。
【００１７】
メッシュ制御パラメータ決定部１０７は、上記メッシュ生成対象特徴パラメータ１１５と上記応答曲面データベース１０５に格納されている応答曲面構築パラメータを入力し、上記メッシュ生成対象形状に適したメッシュ制御パラメータを対象形状メッシュ制御パラメータ１１６（第２のメッシュ制御パラメータ）として決定する。
【００１８】
メッシュ生成部１０８は、上記メッシュ生成対象形状モデル１１４と上記決定された対象形状メッシュ制御パラメータ１１６を用いて、対象形状メッシュ１１７を生成する。
以下、図１および他の図も参照して、本装置の各部について詳細に説明する。 
【００１９】
（１）代表形状特徴パラメータ取得部１０２
代表形状特徴パラメータ取得部１０２では、本装置利用者が入出力装置１０１を使用して入力した代表形状モデル１０９から、形状特徴パラメータを代表形状特徴パラメータとして取得する。入力する代表形状モデル１０９のデータは、形状を定義する点列形式、ベジェ曲線・曲面データ、B-Spline曲線・曲面データ、NURBS曲線・曲面データ等の関数形式、二次元・三次元CADデータ等の形式を持つ。
【００２０】
また、「形状特徴パラメータ」とは、着目する製品または部品等の形状を定めるパラメータであり、本装置に入力する形状モデルから取得可能なものとする。代表形状モデル１０９は、代表形状の形状モデルである。形状特徴パラメータの具体例としては、例えば、形状モデル中部材の寸法や、円弧半径、曲線関数の定義パラメータ、曲率等がある。
【００２１】
図２は、代表形状特徴パラメータ取得部における画面構成例である。この図２には、代表形状特徴パラメータ取得部１０２に入力する代表形状モデル１０９を指定したときの入出力装置１０１の画面表示例が示されている。また、この際入力する代表形状モデルを２０１とし、それに対し取得した形状特徴パラメータが２０２であったとする。本例では、形状特徴パラメータ２０２として、「板厚」、「ベンド径」が採り上げられている。
【００２２】
（２）代表形状メッシュ制御パラメータ最適化部１０３
代表形状メッシュ制御パラメータ最適化部１０３では、本装置利用者が入出力装置１０１を使用して、上記代表形状特徴パラメータ取得部１０２で入力した代表形状モデル１０９と、解析条件１１０と、実験結果または詳細解析結果である参照結果１１１を入力し、上記代表形状の解析結果の誤差が許容値以下となるような代表形状最適メッシュ制御パラメータを探索する。
【００２３】
解析条件１１０とは、代表形状について、事前にメッシュを生成して、所定の計算機により有限要素法等の数値解析を行ったときに定めた条件であり、例えば、境界条件等である。
【００２４】
参照結果１１１の一つである「実験結果」は、代表形状を実際に測定したときの測定位置と、測定パラメータの組み合わせとして入力し、それらは複数あっても構わない。また、参照結果１１１の一つである「詳細解析結果」も同様であり、解析結果の位置と解析によって得られる何らかの値（物理量）とし、それらは複数あっても構わない。ただし、測定位置、解析結果の位置は各代表形状において等しいとする。なお、詳細解析結果は、一般的には、メッシュを用いた解析とは別の解析を行うことで求めるようにしてもよい。
【００２５】
また、参照結果１１１は、代表形状最適メッシュ制御パラメータを見つけるときの基準となるため、一般的には、高精度な結果であることが好ましい。しかし、実際にメッシュ生成対象形状についてメッシュを生成して数値解析を行う利用者の便宜を考慮して、必ずしも高精度でなくともよく、所定の精度を備えた結果であればよい。
【００２６】
なお、実験結果や詳細解析結果の使用不使用に関わらず、利用者が設定する目的関数を最適化することも可能とする。その際、扱う目的関数の例としては、全体の解析メッシュ数や解析時間、メッシュ品質を示すアスペクト比やスキューネス等とする。
【００２７】
図３は、代表形状メッシュ制御パラメータ最適化部における誤差最小化対象設定の画面構成例である。この図３には、図２の代表形状モデル２０１に対し、入力する参照結果１１１と、参照結果１１１に含まれる測定位置もしくは詳細解析結果の位置を指定する際の画面構成例が示されている。符号３０１、３０２、３０３、３０４は、代表形状の解析結果と参照結果１１１との差異を許容値以下としたい点を示している。この点は複数指定可能であり、参照結果１１１から自動で指定することも、画面上で代表形状モデル２０１の所望の箇所をピックアップして指定することも可能である。
【００２８】
誤差最小化対象設定３０５は、前記差異を許容値以下としたい点の情報を示す。その情報として、点３０１、３０２、３０３、３０４について、その座標値（例：「ｘ座標」、「ｙ座標」、「ｚ座標」）、誤差の最小化対象となる物理量の名称（例：「変位」）が指定されている。ここで座標値、物理量の名称を直接指定することも可能である。
【００２９】
ここで、図４を参照して、代表形状メッシュ制御パラメータ最適化部１０３による代表形状メッシュ制御パラメータ最適化の手順について説明する。
図４は、代表形状メッシュ制御パラメータ最適化の手順を示すフローチャートである。
【００３０】
まず、後記する応答曲面データベース１０５の構築に用いる代表形状モデル１０９を、入出力装置１０１に入力し、ステップＳ４０１において、例えば、実験計画法を利用してメッシュ制御パラメータを多種多様に変更して（変化させて）メッシュ生成を行う。生成したメッシュは、「代表形状メッシュ」と称する。
【００３１】
ここで、「メッシュ制御パラメータ」とは、メッシュ分割数、メッシュ粗密、メッシュサイズ、メッシュトポロジー等の解析用メッシュを生成する際に必要なパラメータであり、実験計画法を用いる際にはそれらの値を変化させる範囲を設定する。また、実験計画法としてはラテン超方格法等を利用し、上記メッシュ制御パラメータの変化の範囲内で偏りの無いメッシュ制御パラメータの組み合わせを決定し、それらに基づいて、１以上の代表形状メッシュを生成する。
【００３２】
次に、ステップＳ４０２において、上記実験計画法によりメッシュ制御パラメータを様々に変化させて生成した代表形状メッシュと、それらの解析に用いる解析条件１１０を用いて解析を実行し、それぞれのメッシュについての解析結果を取得する。取得した解析結果は、「代表形状解析結果」と称する。
【００３３】
次に、ステップＳ４０３において、代表形状解析結果と参照結果１１１との差異の最小化、もしくは、その他装置利用者が設定する目的関数値の最小化を行う。参照結果１１１となる実験結果または詳細解析結果と、上記得られた各代表形状解析結果を入力し、単数もしくは複数の実験での測定位置や詳細解析結果の比較位置での結果と、各代表形状解析結果との差の絶対値をとり、目的関数値とする。また、装置利用者が入力部から入力して設定するその他の目的関数、例えば全体解析メッシュ数や解析時間等についても、所望する全体解析メッシュ数もしくは解析時間との差の絶対値等をとる。もし、所望する値が無く、どんな値でも構わないから単純に最小化したいという場合にも、そのような値を目的関数値として取得する。
【００３４】
それら単数もしくは複数の目的関数値に対し、メッシュ分割数、メッシュ粗密、メッシュサイズ、メッシュトポロジー等のメッシュ制御パラメータを設計変数とし、目的関数値を最小化する最適化問題を解く。最適化問題を解くアルゴリズムには、例えば最急降下法、共役勾配法、焼きなまし法、遺伝的アルゴリズム（ＧＡ）、Particle Swarm Optimization法等を用いる。
【００３５】
目的関数が単数の場合は最適なメッシュ制御パラメータの組み合わせが一意に求められるが、目的関数が複数の場合は、多目的最適化問題として解いてメッシュ制御パラメータのパレート解を求め、装置利用者の意志により最適なメッシュ制御パラメータを決定することも可能である。その一方、複数の目的関数を線形結合して単目的最適化問題として解き、目的関数が単数の場合と同様に最適なメッシュ制御パラメータの組み合わせを一意に求めることも可能である。
【００３６】
また、設計変数としてのメッシュ制御パラメータが多数存在する場合、メッシュ制御パラメータの組み合わせを入力とし、上記単数もしくは複数の目的関数値を出力する応答曲面を構築することにより、代表形状の解析ケース数を減らすことも可能である。応答曲面の構築法としては、Kriging法、Radial Basis Function法、多項式近似、ニューラルネットワーク等を用いる。なお、ここで述べた応答曲面の構築は、応答曲面構築部１０４で行われる応答曲面の構築とは別である。
【００３７】
以上の手順により得られたメッシュ制御パラメータの組み合わせを、「代表形状最適メッシュ制御パラメータ」として出力する。また、代表形状特徴パラメータと最適メッシュ制御パラメータとの組み合わせを１１２とし、応答曲面構築部１０４の入力とする。
【００３８】
ここで、図５を参照して、図３で指定した代表形状モデル２０１、参照結果に対し、最適化問題を定義する方法について説明する。
図５は、代表形状メッシュ制御パラメータ最適化部における最適化問題定義の画面構成例である。この画面には、主に、目的関数定義５０１、許容誤差５０２、ユーザ定義目的関数５０３、メッシュ制御パラメータ定義５０４、最適化アルゴリズム５０５、最適化個体数５０６が表示されている。
【００３９】
まず、目的関数定義５０１では、誤差最小化対象設定３０５で設定した誤差最小化対象について、それぞれの最適化を行うか否かを指定する。その際、最小化対象となる物理量の誤差の許容範囲を許容誤差５０２で指定する。
また、ユーザ定義目的関数５０３では、解析結果誤差以外に、装置利用者が設定するその他の目的関数を指定する。
【００４０】
また、メッシュ制御パラメータ定義５０４では、最適な組み合わせを探索したいメッシュ制御パラメータの種類と、その最小値および最大値を指定する。
また、最適化アルゴリズム５０５では最適メッシュ制御パラメータ探索に用いる最適化アルゴリズムの種類を指定し、最適化個体数５０６では最適化に用いる個体数を指定する。
【００４１】
メッシュ制御パラメータ定義５０４、最適化アルゴリズム５０５、最適化個体数５０６の情報に従って、メッシュ制御パラメータ定義５０４の組み合わせを様々に変化させた最適化個体を、最適化個体数５０６で指定した個数分作成し、最適化アルゴリズム５０５で指定した最適化アルゴリズムによって、入力された代表形状モデル２０１に対し、代表形状最適メッシュ制御パラメータを探索する。
【００４２】
（３）応答曲面構築部１０４
応答曲面構築部１０４では、１以上の代表形状に対し、代表形状特徴パラメータ取得部１０２で得られた代表形状特徴パラメータと、代表形状メッシュ制御パラメータ最適化部１０３で得られた代表形状最適メッシュ制御パラメータの情報の組み合わせ１１２から、形状特徴パラメータを入力することで、それに応じたメッシュ制御パラメータを出力する機能を備える応答曲面の構築を行う。この際、応答曲面の構築法として、Kriging法、Radial Basis Function法、多項式近似、ニューラルネットワーク等を用いる。応答曲面の構築の結果、応答曲面構築パラメータ１１３を得る。得られた応答曲面構築パラメータ１１３は、応答曲面データベース１０５に格納される。
【００４３】
図６は、複数代表形状に対する形状特徴パラメータ、最適メッシュ制御パラメータの取得例である。また、図７は、応答曲面構築パラメータの具体例である。
【００４４】
図６には、（ａ）複数の代表形状６０１、６０２、６０３、６０４、６０５、６０６について、符号６０８、６０９、６１０、６１１で示した実験結果の測定位置と、（ｂ）実験結果６０７として示した実験結果の変位量と、（ｃ）代表形状特徴パラメータ６１２と、（ｄ）代表形状最適メッシュ生成パラメータ６１３とが、結果例として図示されている。
【００４５】
また、図７には、代表形状特徴パラメータ６１２および代表形状最適メッシュ生成パラメータ６１３について、Kriging法を用いた際に得られた応答曲面構築パラメータ７０１が一例として示されている。
【００４６】
（４）応答曲面データベース１０５
応答曲面データベース１０５には、応答曲面構築部１０４で得られた応答曲面構築パラメータ１１３を格納する。ここで例として、図７の符号７０１で示した応答曲面構築パラメータ７０１を格納したとする。
【００４７】
（５）メッシュ生成対象形状特徴パラメータ取得部１０６
メッシュ生成対象形状特徴パラメータ取得部１０６では、新たな解析対象となる形状モデル、つまり、メッシュ生成対象形状モデル１１４を入力し、代表形状特徴パラメータ取得部１０２と同じ仕組みにより、形状特徴パラメータとしてメッシュ生成対象形状特徴パラメータ１１５を取得する。取得したメッシュ生成対象形状特徴パラメータ１１５は、代表形状特徴パラメータ取得部１０２で取得した代表形状特徴パラメータと同じものとする。
【００４８】
図８は、メッシュ生成対象形状モデルとメッシュ生成対象形状特徴パラメータの具体例である。図８（ａ）には、メッシュ生成対象形状モデル８０１と、図８（ｂ）には、それに対し取得した形状特徴パラメータ８０２を示す。本例では、形状特徴パラメータ８０２として、「板厚」、「ベンド径」が採り上げられている。
【００４９】
（６）メッシュ制御パラメータ決定部１０７
メッシュ制御パラメータ決定部１０７では、メッシュ生成対象形状特徴パラメータ取得部１０６で取得した、メッシュ生成対象形状特徴パラメータ１１５を、応答曲面データベース１０５に格納されている応答曲面構築パラメータ１１３から構築した応答曲面の入力とし、メッシュ生成対象形状に対する応答曲面上の対象形状メッシュ制御パラメータ１１６を出力する。
【００５０】
図９は、対象形状メッシュ制御パラメータの具体例として、図８のメッシュ生成対象形状モデル８０１について応答曲面を利用して決定された対象形状メッシュ制御パラメータ９０１を示す。対象形状メッシュ制御パラメータ９０１の種類は、メッシュ制御パラメータ定義５０４に示した種類と同じになる。
【００５１】
（７）メッシュ生成部１０８
メッシュ生成部１０８では、新たにメッシュ生成対象となる形状のモデル、つまりメッシュ生成対象形状モデル１１４について、メッシュ制御パラメータ決定部１０７で決定された対象形状メッシュ制御パラメータ１１６を用いてメッシュ生成を行う。生成されたメッシュは、対象形状メッシュ１１７として、その後に行われる有限要素法等の数値解析に用いられる。
【００５２】
図１０は、対象形状メッシュの具体例として、メッシュ生成対象形状モデル８０１について、対象形状メッシュ制御パラメータ９０１を用いてメッシュ生成を行った結果物である対象形状メッシュ１００１を示す。
【００５３】
以上のメッシュ生成装置により、メッシュ生成対象形状の高精度、またはその他装置利用者が設定した全体解析メッシュ数や、解析時間等といった目的関数値が良好なメッシュを生成するためのメッシュ制御パラメータの組み合わせを、高速に決定しメッシュ生成を自動的に行うことが可能となる。代表形状について、高精度に解析されたときの応答曲面を事前に蓄えてあるため、メッシュの生成対象となる新たな形状について、その応答曲面を用いることで、高精度に解析できるメッシュを高速に生成することができる。よって、解析モデル化最適制御に必要な時間的コストを削減することができる。
【００５４】
また、前記した特許文献１には、解析モデルの作成支援を行うとはいえ、個々のモデル化に対し、より精度の高い解析に必要なメッシュ粗密の制御を実行する技術は含まれていない。しかし、本実施形態のメッシュ生成装置によれば、代表形状に用いるメッシュ制御パラメータを変更すれば、相応の応答曲面ができあがるので、メッシュ生成対象形状に用いるメッシュの粗密を実現することができる。
【００５５】
≪使用例≫
以下、その他の本装置の使用例について、メッシュモーフィングについて説明する。「メッシュモーフィング」とは、既存のメッシュを変形して新形状のメッシュを生成する方法をいい、本装置の一連の利用手順の１例である。
【００５６】
図１１は、本装置をメッシュモーフィングによるメッシュ生成に適用した場合の具体例である。図１１（ａ）は、代表形状１１０１、１１０２、１１０３である。図１１（ｂ）は、代表形状特徴パラメータ１１０９である。図１１（ｃ）は、代表形状最適メッシュ制御パラメータ１１０４である。図１１（ｄ）は、メッシュ生成対象形状１１０５である。図１１（ｅ）は、メッシュ生成対象形状特徴パラメータ１１０６である。図１１（ｆ）は、対象形状メッシュ制御パラメータ１１０７である。図１１（ｇ）は、対象形状メッシュ１１０８である。
【００５７】
本使用例は、平板に開けられた穴の位置を上下左右それぞれに0.0〜10mmの範囲で変更する場合のメッシュ制御パラメータを決定する例として説明する。この説明では、図１１（ａ）に示す代表形状１１０１〜１１０３およびメッシュ生成対象形状１１０５は、メッシュとして与えられたものとする（代表メッシュ、変形対象メッシュ）。
【００５８】
本装置では、まず代表形状特徴パラメータ取得部１０２において代表形状特徴パラメータ（第１の形状特徴パラメータ）を取得する。本例では形状特徴パラメータとして、「穴径」を取得できるものとし、１１０１〜１１０３に示すメッシュの代表形状に対して代表形状特徴パラメータ取得部１０２により取得した結果を、図１１（ｂ）の代表形状特徴パラメータ１１０９に示す。また、穴位置の「移動方向」および「移動量」も形状特徴パラメータとする。
【００５９】
ここで穴位置を上下左右それぞれに0.0〜10mmの範囲で変更する場合の代表形状最適メッシュ制御パラメータ（最適メッシュ制御パラメータ）を代表形状メッシュ制御パラメータ最適化部１０３で探索する。メッシュ制御パラメータ（第１のメッシュ制御パラメータ）としては、移動端側テンション、固定端側テンションを用いることとする。
【００６０】
「移動端側テンション」は、メッシュ変形時に移動端側に節点を寄せる割合であり、-1.0〜1.0の範囲で指定する。移動端側テンションが示す値が大きいほど移動端側に節点が寄る。「移動端側」とは、メッシュ変形時に移動する節点が向かう側をいう。
「固定端側テンション」は、メッシュ変形時に固定端側に節点を寄せる割合であり、-1.0〜1.0の範囲で指定する。固定端側テンションが大きい値ほど固定端側に節点がよる。「固定端側」とは、メッシュ変形時に移動する節点が、移動しない節点に向かう側をいう。
【００６１】
また、メッシュ変形においては、節点可動範囲が定められている。「節点可動範囲」とは、穴中心点を中心とする正方形の１辺の長さをいう。
また、本例では、目的関数としては、メッシュ品質として要素アスペクト比を用いることとする。
【００６２】
前記探索の結果、代表形状１１０１〜１１０３の代表形状最適メッシュ制御パラメータ１１０４として、図１１（ｃ）に示す値が得られる。図１１（ｃ）には、移動方向、上下左右それぞれに5mm、10mmとした際の代表形状最適メッシュ制御パラメータが表として示されている。
【００６３】
次に本装置は、応答曲面構築部１０４において、代表形状特徴パラメータと最適メッシュ制御パラメータに関する応答曲面を構築し、その応答曲面の応答曲面構築パラメータを応答曲面データベース１０５に登録する。
【００６４】
次に本装置は、メッシュ生成対象形状特徴パラメータ取得部１０６において、メッシュ生成対象形状特徴パラメータ（第２の形状特徴パラメータ）を取得する。本例では図１１（ｄ）に示すメッシュ生成対象形状１１０５がメッシュとして入力されたものとする。また、移動方向は右、移動量は7mmといった具合に形状特徴パラメータを与えた場合、メッシュ生成対象形状特徴パラメータ取得部１０６により取得した結果を、図１１（ｅ）のメッシュ生成対象形状特徴パラメータ１１０６に表として示す。
【００６５】
次に本装置は、メッシュ制御パラメータ決定部１０７において、メッシュ生成対象形状特徴パラメータ１１０６と応答曲面データベース１０５に格納されている応答曲面構築パラメータ１１３からメッシュ制御パラメータ（第２のメッシュ制御パラメータ）を決定する。その決定した結果を、図１１（ｆ）のメッシュ制御パラメータ１１０７に表として示す。
【００６６】
最後に本装置は、メッシュ生成部１０８において、メッシュ制御パラメータ１１０７を用いて対象形状メッシュを生成（モーフィング）する。その結果、図１１（ｆ）に示す対象形状メッシュ１１０８が生成される。
【００６７】
このように、本実施形態では、既存のメッシュを変形することで、新規のメッシュを高速に生成することができる。
【００６８】
≪その他≫
前記実施形態は、本発明を実施するために好適のものであるが、その実施形式はこれらに限定されるものでなく、本発明の要旨を変更しない範囲内において種々変形することが可能である。
【００６９】
例えば、本実施形態では、代表形状は、有形物（またはメッシュ）を取り扱ったが、これに限らず、例えば、流体解析や電磁場解析を行うために無形物を取り扱うこともできる。
【００７０】
また、本実施形態で説明した種々の技術を適宜組み合わせた技術を実現することもできる。
【００７１】
その他、ハードウェア、ソフトウェア、各フローチャート等の具体的な構成について、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。
【符号の説明】
【００７２】
１０１  入出力装置（入力部、制御部、記憶部を備える）
１０２  代表形状特徴パラメータ取得部
１０３  代表形状メッシュ制御パラメータ最適化部
１０４  応答曲面構築部
１０５  応答曲面データベース
１０６  メッシュ生成対象形状特徴パラメータ取得部
１０７  メッシュ制御パラメータ決定部
１０８  メッシュ生成部
１０９  代表形状モデル
１１０  解析条件
１１１  参照結果
１１２  代表形状特徴パラメータ（第１の形状特徴パラメータ）と代表形状最適メッシュ制御パラメータ（最適メッシュ制御パラメータ）の組み合わせ
１１３  応答曲面構築パラメータ
１１４  メッシュ生成対象形状モデル
１１５  メッシュ生成対象形状特徴パラメータ（第２の形状特徴パラメータ）
１１６  対象形状メッシュ制御パラメータ（第２のメッシュ制御パラメータ）
１１７  対象形状メッシュ
【図１】
2011243016.tif 000003
【図２】
2011243016.tif 000004
【図３】
2011243016.tif 000005
【図４】
2011243016.tif 000006
【図５】
2011243016.tif 000007
【図６】
2011243016.tif 000008
【図７】
2011243016.tif 000009
【図８】
2011243016.tif 000010
【図９】
2011243016.tif 000011
【図１０】
2011243016.tif 000012
【図１１】
2011243016.tif 000013