JP4639292B2 - 3次元メッシュ生成方法 - Google Patents
3次元メッシュ生成方法 Download PDFInfo
- Publication number
- JP4639292B2 JP4639292B2 JP2002168538A JP2002168538A JP4639292B2 JP 4639292 B2 JP4639292 B2 JP 4639292B2 JP 2002168538 A JP2002168538 A JP 2002168538A JP 2002168538 A JP2002168538 A JP 2002168538A JP 4639292 B2 JP4639292 B2 JP 4639292B2
- Authority
- JP
- Japan
- Prior art keywords
- mesh
- dimensional
- model
- shape
- analysis
- 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.)
- Expired - Lifetime
Links
Images
Description
【発明の属する技術分野】
本発明は、3次元メッシュ生成方法に関し、より詳しくは、3次元形状のオリジナル・データを、所定のアスペクト比の範囲内で、かつ、所定の形状近似誤差範囲内となるように、四辺形メッシュ及び三角形メッシュで自動的に再分割する3次元メッシュ生成方法に関する。
【0002】
【従来の技術】
数値データに基づく数値解析計算により様々な物理現象や物性を解明し、これにより製品の性能や信頼性を効率良く高かめる支援技術としてCAEシステムが知られているが、3次元ソリッドモデルCADソフトによるオリジナル・形状モデルデータや、非接触3次元計測オリジナル・点群データを、そのまま、直接、CAEの数値解析ソフトに入力データとして入力できないので、FEM(Finite Element Method)等の数値解析ソフトのために解析用メッシュ・データを新規に生成(Healing)する必要がある。かかる解析用メッシュとしては、数値計算の精度を向上させるために、3次元モデルには六面体メッシュの採用が好ましいが、解析用メッシュの総要素個数を増加させたくない場合には、四面体メッシュを併用するのが好ましい。
従来、上記六面体(四角形)メッシュのみを生成する場合、曲線座標変換法により直交格子状の写像モデルの格子点を解析対象形状に写像し解析メッシュを生成する写像法が実用になっており、解析対象の形状モデルから解析メッシュ生成のための写像モデルを自動生成する方法として、具体的には以下の技術が知られている。
a1) 形状モデルの面や稜線の接続情報(位相情報)を用いて、形状モデルに適合した写像モデルを自動生成する方法では、特開平1−311373号公報(日立)に、「メッシング方法及び装置」が記載されている。又、
a2) 自動生成した写像モデルを端末画面上に表示し、これを対話的に修正しメッシュの要素形状や粗密を制御する方法では、特開平7−127754号公報(日立)に、「数値解析用メッシュ生成方法及び生成装置」が記載されている。更に、
a3) 3次元シェル形状モデルを構成する各面のメッシュ自動分割を行う方法では、特開平6−180738号公報(日立)に、「有限要素生成装置、有限要素生成方法」が記載されている。又、
a4) 歪みの少ない要素で構成される六面体(四角形)メッシュを自動生成できる解析メッシュ生成方法では、特開平11−144093号公報(日立)に、「解析メッシュ生成方法及び装置」が記載されている。
【0003】
【発明が解決しようとする課題】
上記a1)〜a4)の従来方法では、曲線座標変換法を利用するため、いずれの手法も、上記入力した3次元空間等の形状モデルを、直接、曲線座標変換法の基礎となる単位正方形又は単位立方体で分割する必要があり、かかる形状モデルを単位立方体(単位正方形を含む)で直接分割する手法では、単位立方体の生成/分割に非常に演算時間及び作業メモリを浪費し、高速なメッシュ生成はとても望めないという根本的な問題点があった。又、曲線座標変換法では、非常に長時間の演算時間を必要とし、実用に耐える高速なメッシュ生成方法又は装置は、とても望めなかった。かかる問題点は、解析対象物が大きくなると、指数関数的に、演算時間や修正作業が増加するので、複雑な構造物は、実際には、そもそも有限要素解析の対象とできないという問題点ともなっていた。又、上記曲線座標変換法では、写像モデルは、必ず、六面体要素(2次元の場合には、四角形要素)からなる解析用メッシュで構成されなければならず、更に、対向する四角形の対辺の格子点数は、必ず同一の格子点数となるように生成されなければならないので、円錐形状等では、形状誤差が、上下の長さの異なる対辺の組の分割メッシュに自動的に含まれてしまうといった問題点もあった。
更に、3次元測定機等から入力した3次元計測オリジナル・データ点群からなる複雑な形状モデルに対しては、上記曲線座標変換法は、直接適用できなかった。
よって、本発明は上述の様な事情に鑑みて成されたものであり、本発明の目的は、従来の四面体要素のみの自動メッシュ生成方法、六面体要素のみの自動メッシュ生成方法では実現不可能であった3次元測定機等から入力した3次元計測オリジナル・データ点群からなる不規則/複雑な形状でも、計算精度を十分保った状態で、解析メッシュ要素の作成が実用に耐える高速な速度で可能であり、3次元計測オリジナル・データ点群で規定される面の形状・大きさ等面の品質に左右されることなく、品質の良いメッシュを生成し、精度の高い数値計算用のメッシュモデルを作成する方法/装置を提供することにあり、歪みの少ない要素で構成される六面体(四角形)メッシュ/四面体(三角形)メッシュを自動生成できる実用的な高速解析メッシュ生成方法及び装置を提供することにある。
【0004】
【課題を解決するための手段】
本発明は、解析対象の3次元形状モデルと解析メッシュ生成用基礎データとを入力する形状入力部と、写像モデルの格子点を3次元形状モデルに写像して解析メッシュを生成し、生成された解析メッシュをデータベースに格納するメッシュ生成部と、を含む解析メッシュ生成装置におけるメッシュ生成方法に関し、本発明の上記目的は、
前記解析メッシュ生成装置が、
前記3次元形状モデルを、所定の投影方向に対応して予め設定されている投影基準面に対して当該投影方向から投影し、当該基準面に前記3次元形状モデルの最外郭線の内部からなる基準面モデルを生成する処理を実行する第1工程と、前記基準面モデルの内部に、所定のドット生成規則により生成された格子点からなる2次元写像モデルを生成する処理を実行する第2工程と、を含み、
前記解析メッシュ生成装置が、
前記2次元写像モデルの格子点を、前記投影基準面から、仮想平行光線により前記3次元形状モデルに投影して解析メッシュを生成する処理を実行する第3工程と、
前記解析メッシュを後記分割パターンにより仮想的に分割したときの分割メッシュを構成する各仮想分割辺にて隣り合う2辺の比が最大の比をアスペクト比とし、前記分割メッシュを構成する平面と前記3次元形状モデルの曲面との間隔の最大値を3次元形状の近似誤差として、前記分割メッシュに対して、前記アスペクト比、前記仮想分割辺の長さのバランス、ワープ歪、及び、要素形状の最大角度と最小角度がそれぞれ所定の範囲内で、かつ、前記3次元形状の近似誤差が所定の範囲内であることを条件として、予め用意した四辺形メッシュ及び/又は三角形メッシュの分割パターンに基づいて、前記条件を満たしていない要素メッシュを分割する処理を実行する第4工程と、
前記第3工程での前記2次元写像モデルの格子点の処理対象を前記第4工程による分割後の要素メッシュの格子点として前記第3工程に戻り、前記条件を満たすまで前記第3工程と前記第4工程とを繰り返し実行し、前記条件を満たしたと判定したときに、前記第3工程により生成された解析メッシュを前記投影基準面に対する第1の最終3次元メッシュとして前記データベースに格納すると共に、前記所定の投影方向として複数の投影方向が設定されている場合には、前記第1工程から前記第4工程の処理によって前記第1の最終3次元メッシュを前記複数の投影方向毎に生成し、生成された第1の最終3次元メッシュ同士を相互に連結して1つの解析メッシュに統合する処理を実行し、第2の最終3次元メッシュとして前記データベースに格納する第5工程と、を更に含むことによって達成される。
本発明では、対象となる3次元形状モデルの非接触計測オリジナル・データやソリッド・モデルを、ユーザーが指定した基準面又は自動的に計算した貼り付け方向を法線とする無限の投影基準平面に投影し、当該基準面に投影された最外郭線の境界及びその内部に、投影用の所定の間隔で構成された格子点パターン(初期メッシュ)を生成する。次に、基準面上の格子点パターンを3次元形状モデルの表面に投影し、モデル表面の形状に合わせて貼り付ける。もし、モデル表面の要素形状に所定の閾値以上の歪が検出された場合、要素を再分割して形状を整えることで、数値解析可能なメッシュ生成が可能となる。複雑な3次元対象物の場合には、基準面の投影方向を複数個設定し、それぞれの基準面毎に、格子点パターンを生成し、その後、元の3次元形状モデルに、それぞれ基準面毎の格子点パターンを投影し、それぞれ生成されたメッシュ同士を相互に連結、統合することで最終的な3次元解析用メッシュが生成できる。
【0005】
【発明の実施の形態】
以下、図面に基づいて、本発明の好適な実施例について、詳細に説明する。
図1は、本実施例の解析メッシュ生成装置の全体システム構成の説明図である。
図2は、メッシュ生成処理の流れの説明図である。図3は、形状モデル、写像モデル及び解析メッシュの一例の説明図である。
本発明の解析メッシュ生成装置は、図1に示すように、3次元測定機1と、形状入力部2と、写像用2次元初期メッシュ生成部(写像モデル生成部)4と、3次元メッシュ生成部6と、要素形状評価部8と、写像モデル修正部10と、メッシュ出力部12と、メッシュ統合部14と、データベース16と、入出力装置18と、から構成される。尚、3次元測定機1は、3次元ソリッド・モデルCADソフトからオリジナル・3次元形状データを入力する場合には不要である。又、形状入力部2〜メッシュ統合部14は、コンピュータを含み、装置内部の制御を行うCPU(Central Processing Unit:中央演算装置)/MPUや、主記憶装置、固体メモリからなる補助記憶装置を含み、ソフトウエア・プログラムとして、一般には、実現するのが好ましい。
形状入力部2は、解析対象の3次元形状モデルのオリジナル・データ、解析メッシュの目標要素寸法及び要素の歪みを判定するパラメータ等を入力する。入力した形状モデルデータ等はデータベース16に格納される。
解析メッシュの目標要素寸法及び要素の歪みを判定するパラメータは、解析メッシュ生成用基礎データである。解析メッシュ生成用基礎データとして解析メッシュの格子点の個数等を使用することも可能である。
初期メッシュ生成部4は、入力した形状モデルデータと解析メッシュの目標要素寸法とからメッシュ生成のための四角形メッシュ及び三角形メッシュの混在した2次元平面状の写像モデルを投影基準面に生成する。3次元メッシュ生成部6は、仮想平行光線により投影基準面に生成した写像モデルの2次元平面状格子点を3次元オリジナル・形状モデルに投影して平面状の3次元解析用要素メッシュを生成する。従って、基準面の三角形メッシュからは、三角形要素形状を含む四面体/五面体メッシュが生成され、基準面の四角形メッシュからは、四角形形状要素を含む五面体/六面体メッシュが生成される。要素形状評価部8は、モデル表面に生成された形状要素の歪みを判定するパラメータを使用して生成された三次元メッシュの各要素の歪みを算出し、形状歪みが所定の範囲内か否かチェックし、形状歪の大きい要素を検出する。写像モデル修正部10は、要素形状評価部8で歪みの大きい形状要素を検出すると、解析メッシュの要素歪みの大きい要素に相当する写像モデルを再分割することにより修正する。メッシュ出力部12は、生成された解析メッシュを表示する。メッシュ統合部14は、複数の投影方向から生成された解析メッシュ・データを1つの解析メッシュ・データに統合する。データベース16は、形状モデル、写像モデル及び解析メッシュのデータを格納する。入出力装置18は、表示装置(CRT,液晶表示手段等)18a、キーボード、マウス等のポインティング手段18b、LAN、インターネット等のネットワークを介した通信手段18c、CD,MO、DVD等の大容量記憶媒体18d等の読込/書込装置を含む。
【0006】
次に、3次元オリジナル・形状モデルの入力から3次元解析用要素メッシュ出力までの処理の流れを図2を用いて説明する。本発明では、
解析対象の3次元形状モデルと解析メッシュ生成用基礎データとを入力する工程と、写像モデルの格子点を3次元形状モデルに写像して解析メッシュを生成する工程と、を含むメッシュ生成方法において、
前記写像モデルの生成工程が、
前記形状モデルを、所定の投影基準面に投影し、当該基準面に前記形状モデルの最外郭線の内部からなる基準面モデルを生成する工程と、前記基準面モデルの内部に、所定のドット生成規則により生成された格子点からなる2次元写像モデルを生成する工程と、を含み、
前記2次元写像モデルの格子点を、前記投影基準面から、仮想平行光線により前記3次元形状モデルに投影して解析メッシュを生成する工程と、
前記生成された解析メッシュに対して、所定のアスペクト比の範囲内で、かつ、所定の近似誤差範囲内となるまで、四辺形メッシュ及び/又は三角形メッシュで再分割する工程と、を更に含んでいる。具体的には、
c1) 入出力装置18を用いて、解析対象の3次元オリジナル形状モデルデータ30a(図3(A))と、生成する解析メッシュの目標要素寸法及び要素歪みを判定するパラメータ等の要素形状基準値のデータを含む解析メッシュ生成用基礎データとを入力する(S2)。解析メッシュ生成用基礎データについては後述する。
c2) 入力された形状モデル30aから写像モデルの2次元初期メッシュ32a(図3(C))を、所定の投影方向の投影基準面毎に生成する(S4)。ここでは、先ず、
d1) 形状モデルを、所定の投影基準面方向毎に投影し、最外郭輪郭線を生成する。尚、投影基準面は、±X方向・±Y方向・±Z方向の6つの方向を含む複数の方向に分けることが好ましいが、−Z方向等の単一方向とすることも可能である。
d2) 次に、各投影基準面方向毎に、最外郭輪郭線の内部に、いわゆる、シルエット領域を生成し、このシルエット領域に、目標要素寸法から指定された間隔で、正方格子状に、四角形メッシュ及び/又は三角形メッシュからなる写像用2次元初期メッシュを生成する(S4)。
e1) この2次元初期メッシュにより生成された投影基準面上の格子点を、3次元オリジナル形状モデルに投影し、六面体、五面体及び/又は四面体からなる解析メッシュ34を生成する(S6)。上記投影は、仮想平行光線により、投影基準面上の格子点を3次元形状モデル表面に投影し、3次元の形状要素を生成する。
e2) 次に、要素形状の歪みを算出する(S8)。要素形状の歪みは、予め四角形メッシュ又は三角形メッシュに対して、複数の歪算出モデルを用意しておき、要素メッシュ形状を構成する平面と一般には曲面からなるオリジナル形状モデル表面との適合度(バランス具合)を基本とする。
e3) 要素歪みが要素形状基準値以下であるか否か判定し(S10)、要素形状基準値以上の形状歪が大きい要素が検出されると、かかる形状要素を、投影基準面の写像モデルに戻って、形状歪が最小となる分割モデルで再分割し(S12)、上記ステップe1)に戻る。
e4) 全ての要素歪みが要素形状基準値以下であると判断されたとき、生成されたメッシュを、所定の投影基準面に対する第1の最終メッシュとして記憶する(S14)。
e5) 上記e1)〜e4)の処理を、全ての投影基準面に対して行う(S16)。
c3) 投影基準面が複数ある場合には、上記e1)〜e5)の処理により生成された各投影基準面に対する第1の最終メッシュを、合成・統合し、第2の最終メッシュとして記憶する(S18)。
【0007】
以下、各処理について詳細に説明する。先ず、入出力装置18を用いて、図3(A)に示すような解析対象の3次元オリジナル形状モデル30a、生成する解析メッシュの目標要素寸法及び要素歪みを判定するパラメータ等の要素形状基準値のデータを含む解析メッシュ生成用基礎データとを入力する(S2)。
尚、本発明の3次元オリジナル形状モデルは、3次元レンジファインダ、3次元レーザスキャナや断層画像撮像装置等の3次元非接触計測器を用いて対象物体を計測した非接触計測点データからなる点群データや、3次元接触式測定機を用いて物体を計測した計測点データからなる点群データを含み、更に、3次元ソリッド・モデルCADソフトで作成された中実のCADデータを含み、CADデータの形式としては、IGES(Initial Graphics Exchange Specification)が主に用いられているがこれに限られない。
解析メッシュ生成用基礎データとしては、次のものがある。
f2) 投影基準面の数g1、及び、方向g2
f4) 目標要素寸法g4
f6) 最小要素寸法g6
f8) 隣接し、隙間がある立体間を統合する場合の許容間隔g8
f10) 採用する初期メッシュ作成モードの種類g10
f12) 四角形要素形状の辺のバランス具合最大g12/最小g13
f14) 三角形要素形状の辺のバランス具合最大g14/最小g15
f16) ワープ歪の角度g16
f18) 四角形要素形状の最大角度g18/最小角度g19
f20) 三角形要素形状の最大角度g20/最小角度g21
f22) 要素メッシュ形状を構成する平面と形状モデル表面との最大許容距離g22
【0008】
次に、入力された3次元オリジナル形状モデル30aから写像モデル32a用平面状2次元初期メッシュを、所定の投影方向の投影基準面毎に生成する(S4及び図3)。ここでは、先ず、
d1) オリジナル形状モデルを、所定の投影基準面方向毎に投影し、3次元形状モデルのシルエット領域を生成する。尚、投影基準面は、±X方向・±Y方向・±Z方向の6つの方向を含む複数の方向に分割することが好ましいが、−Z方向等の単一方向に設定することも可能である。則ち、本発明では、投影基準面の数g1は、1乃至6の基準面が少なくとも採用可能であり、その方向g2は、−Z方向や、±X方向・±Y方向・±Z方向等が利用可能である。例えば、図3(A)に示すような3次元形状モデル30aを+Z方向の投影基準面に投影すると、同図(B)に示すように、投影基準面31上に、最外郭線30abが生成され、これを上方から観察すると,同図(C)のシルエット領域となる。
【0009】
次に、3次元形状モデルが投影された投影基準面のシルエット領域に、目標要素寸法から指定された格子間隔で、3次元形状モデルに適合した四角形メッシュ及び/又は三角形メッシュからなる写像用2次元初期メッシュを生成する処理について説明する(S4)。かかる初期メッシュの生成方法には、複数の生成方法が利用可能であるが、具体的には、本発明では、初期メッシュ作成モードの種類g10は、少なくとも以下の3種類の初期メッシュ生成モードがある。
(h1) 正方格子式初期メッシュ生成モード
図4(A)に示すような投影基準面31に投影した3次元計測点群データからなる形状モデルに対し、目標要素寸法g4の間隔で、投影基準面上で、その2次元シルエットを格子点で分割し、投影基準面内に形状モデルに適合した四角形メッシュ及び/又は三角形メッシュからなる写像モデルの正方格子式初期メッシュを生成する処理について説明する。図4(A)では、最外郭境界線30cb1の内部に窪んだ領域30cb2が形成されている。
h11) 先ず、x-y方向の所定の間隔g4で、正方格子状に、無限平面に広がるドット・パターンで仮想メッシュを生成する。
h12) 又、形状モデルを所定の投影基準面上に投影し、最外郭境界線30cb1を抽出し、この最外郭境界線の内部に、中空部と中実部とを区別しないシルエット領域(図4(B))を生成する。
h13) 次に、上記シルエット領域に対して、正方格子状のドット生成規則により生成された格子点からなる仮想メッシュとシルエット領域との論理積演算を行い、シルエット領域(シルエット領域の内部領域、及び、最外郭境界線の近傍領域を含む)に存在する正方格子点(図4(B)に示すような内部が空白の格子点)を生成する。
h14) 又、上記シルエット領域の中で、特に、最外郭境界線上には、所定の目標要素寸法g4の間隔で、制御節点用格子点(図4(B)の点vcb21、vcb22,vcb23,...等)を生成してもよい。更に、シルエット領域の最外郭境界線の近傍領域に存在する格子点を検出し、かかる格子点を総て境界線上に移動させるようにしてもよい。
尚、格子点相互の間隔は、最小要素寸法g6以上、相互に離れているのが好ましく、最小要素寸法g6は、目標要素寸法g4*0.6以上の値が好ましく、最小要素寸法g6は、目標要素寸法g4*0.7以上の値であると、より好ましい。図4(A)の例で格子点を生成すると、同図(B)のような内部を黒く塗り潰した格子点(vcb21,vcb22,vcb23,...)で2次元境界線が分割される。
h15) 上記境界線上に生成した制御節点用格子点と、上記論理積演算によりシルエットの内部領域に対して生成した正方格子点とを論理和演算により合成し(図4(C))、その後、比較・統合して、全ての格子点が所定の距離(最小要素寸法g6)以上に離れているか否かチェックする。
上記外形線及び特徴線の境界線上に生成した制御節点用格子点と、上記論理積演算により生成した内部正方格子点(vgb21,vgb22,vgb23,...)とが、所定の距離以内に接近して存在していた場合には、境界線上に生成した制御節点用格子点を好ましくは優先させて残し、又は、境界線上で移動させ、上記論理積演算によりシルエットの内部領域に対して生成し、かつ、境界線上の制御格子点に接近した正方格子点は、削除する。
h16) 上記境界線上に生成した制御節点用格子点と、近傍チェックをして内部領域に残った正方格子点とを基に、四角形メッシュ、及び/又は、三角形メッシュからなる初期メッシュを生成する(図4(C))。
かくして、正方格子式初期メッシュ生成処理は終了する。則ち、図4(A)に示すような2次元形状モデルに対して、正方格子式初期メッシュ生成モード処理した2次元平面状の初期メッシュの1例は、図4(C)に示すような形状となる。
(h2) その他の初期メッシュ生成モード
h21) 従来から知られているデローニ三角分割法による初期メッシュ生成モード
従来から知られているデローニ三角分割法による初期メッシュを、投影基準面に対して、境界線から2次元三角形メッシュの集合として生成することも当業者には、容易である。
h22) 公知のアドバンシングフロント法を用いた初期メッシュ生成モード
又、公知のアドバンシングフロント法を用いて、投影基準面に、2次元境界線から、四角形メッシュ/三角形メッシュの混在した初期メッシュを生成することも、可能である。
【0010】
続いて、c2)で生成した初期メッシュからなる2次元平面状写像モデルにより生成された投影基準面上の格子点を、3次元形状モデルの表面に投影し、形状モデルの外表面に、六面体、五面体及び/又は四面体からなる3次元解析用要素メッシュ34を生成する(S6)。上記投影処理は、図5(A)に示すように、仮想平行光線により、各格子点を3次元形状モデルの外表面に投影して行う。
図4(A)では、最外郭境界線30cb1の内部に窪んだ領域30cb2が形成されているので、図4(C)のような2次元平面状初期写像モデルにより生成された投影基準面上の格子点を、図4(A)の3次元形状モデルの表面に投影すると、形状モデルの外表面で最外郭境界線30cb1の内部に図4(D)に示すような3次元解析用要素メッシュが生成され、更にその内側の窪んだ領域30cb2の内部に3次元形状モデルの底部要素メッシュms30cが生成される。
本発明では、3次元形状モデルを、一度、投影基準面に投影して四角形メッシュ、及び/又は、三角形メッシュからなる2次元初期メッシュを生成し、この初期メッシュ生成過程で発生され位置の特定された2次元写像モデルの各メッシュを構成する格子点を、再び3次元形状モデルに投影しフィードバックして3次元解析用要素メッシュを生成する。従って、3次元形状モデルの表面では、要素メッシュの表面形状は、最初に、2次元初期メッシュの四角形メッシュ又は三角形メッシュに基づいて生成された四角柱モデル又は三角柱モデルにより切り取られた3次元形状モデルの非接触計測オリジナル・データ点群又は接触計測オリジナル・データ点群やソリッド・モデルの点群を平面近似すると共に、可能であれば所定の範囲内のアスペクト比を満たす六面体モデル、五面体モデル、又は、四面体モデルを生成する。
次に、平面近似した六面体モデル、五面体モデル、又は、四面体モデルの外表面を更にチェックして、3次元形状要素の歪みを算出する。かかる要素形状の歪みは、予め四角形メッシュ又は三角形メッシュに対して、複数の歪算出モデルを用意しておき、要素メッシュ形状を構成する平面と一般には曲面からなるオリジナル形状モデル表面との適合度(バランス具合)を基本とし、要素歪みが要素形状基準値以下であるか否か判定し、要素形状基準値以上の形状歪が大きい要素が検出されると、かかる形状要素を、投影基準面の写像モデルに戻って、形状歪が最小となる分割モデルで再分割し、形状歪が所定の範囲内となるまで再分割を繰り返す。
従って、本発明では、従来の六面体のみからなる写像モデル等は一切不要であり、又、本発明の写像モデルの格子点は、単に平面上に分散して形成されているだけなので、写像モデルの記憶容量は必要最小限ですむと共に、従来の曲線座標変換法のような複雑で、演算時間が膨大に要求される処理が不要となり、形状モデルの表面に、いわゆるシェル・メッシュを生成することが、高速かつ、ワーク・メモリが最小限の状態で、可能である。
尚、上記2次元平面状写像モデルとして生成された投影基準面上の初期メッシュの格子点を、仮想平行光線により、3次元形状モデルに投影した場合、隣接したオリジナルの形状モデルに隙間や重なりがあり、特に、かかる隙間部分に基準面モデルの格子点を投影しても、該当するモデル外表面が存在しないので、最初の投影段階では、隙間部分の格子点位置を確定できず、3次元形状モデル要素を直ちに作成できない場合がある。かかる隙間空間で、かつ、3次元モデル外表面のない領域で、3次元格子点の座標を演算し決定する方法を、次に、図6を参照して説明する。先ず、図6に示すように、形状モデルの面3f1と3f2との間に、間隔L3を隔てて、隙間が形成されている場合、
j1) 2次元格子点を、3次元モデル外表面に投影しても、投影先に外表面がなく、当該2次元格子点(図6のbjk1〜bjk5の元の投影基準面上の格子点)が隙間空間に投影されたことを記録し、元の2次元格子点は動かさず、対応する3次元格子点の座標が、未定であることも、記録する。
j2) j1)で記録された2次元格子点に関して、その格子点の周囲で、3次元モデル外表面上に投影され3次元モデル座標の確定した各格子点(abj1〜abj5及びabk1〜abk5)の座標値から、平均化処理、又は、重心位置演算処理等により、上記未定の3次元座標を演算して求め(例えば、周囲の3次元座標の確定した1〜8個の格子点座標から平均化処理して、隙間に投影された格子点(bjk2)の座標値とする処理では、周囲の3次元座標の確定した単数又は複数個の格子点座標(abj1〜abj3及びabk1〜abk3)から、重心演算や平均化演算処理により、隙間に投影された格子点(bjk2)の座標値とし、3次元格子点座標を移動させる。かかる処理を、j1)で記録された全ての2次元格子点に関して行う。
j3) j1)で記録された2次元格子点が隣接している場合、移動量が所定の範囲内で小さく移動し収束するまで、j1)で記録された全ての2次元格子点に関して、j2)の操作を繰り返す。
かかる空隙部に対する投影例のように、投影に対応する3次元形状モデル表面が無い場合の仮想3次元座標演算処理は、面の存在しない領域(空隙部)の幅L3がメッシュサイズ(g4)より、小さく設定するのが好ましい。
【0011】
次に、3次元形状モデル外表面に沿って生成された3次元平面状の要素メッシュ形状の歪みを、予め設定した複数種類の形状誤差評価パターンに基づいて、それぞれ評価し、算出する(S8)。かかる3次元要素形状の歪みは、予め投影する2次元四角形メッシュ又は2次元三角形メッシュに対して、複数の3次元歪算出モデルを用意しておき、3次元要素メッシュ形状を構成する平面と一般には曲面からなる3次元形状モデル表面との間で、複数の観点から歪の程度を演算し評価する。尚、ここでは理解し易くするために、2次元形状モデルで説明しているが3次元形状モデルへの拡張は容易である。
k1) 要素形状の表面が四角形メッシュの場合
2次元四角形メッシュの投影により形成された3次元形状要素に対しては、予め登録した以下の複数種類の形状誤差評価パターンを用意し、それぞれ形状誤差を評価し、ユーザが予め設定した許容誤差の範囲内で、総合的に、最も形状誤差が少なくなるように、必要であれば、当該四角形メッシュを再分割する。
k11) ワープ歪
3次元形状モデル表面上の四角形要素メッシュmsh1の曲がり具合を、図7に示すように、三角形メッシュk11aーk11b、k11c-k11dのような2組の対角線状に四角形を分割し形成した三角形形状要素により仮想的に分割し、これら2つの仮想三角形形状要素の平面が形成する角度で、それぞれ3次元形状モデルとの歪を評価し、ユーザが予め設定した許容誤差の範囲外であれば、大きい歪評価角度を形成する分割線b12-b14又はb11-b13に沿って、当該四角形メッシュを、三角形メッシュに再分割する。
上記仮想三角形形状要素の平面が形成するワープ歪の角度g16は、0.0度が最適であり、三角形形状要素の間で形成される歪角度が大きいほど、解析精度が低下するが、例えば、60度以内に設定するのが好ましく、55度以内に設定するとより好ましく、50度以内に設定すると更に好ましい。
かかる再分割処理により、近傍節点が四角形メッシュの対角線上にある場合、四角形メッシュが、2つの三角形メッシュに分割され、境界部分等に三角形メッシュが含まれるので、境界形状に、歪んだり潰れた四角形要素メッシュが生成されず、非常に高精度のメッシュモデルが生成可能となる。
k12) 四角形メッシュの仮想分割辺の長さのバランス及びアスペクト比
3次元形状モデルへの投影により生成された四角形要素メッシュの仮想分割辺の長さのバランスは、図8〜図11に示すように、基本的には、3種類の仮想分割パターンk121〜k123により評価する。尚、上記仮想分割パターンは、仮想分割する方向により、大きく変動するため、それぞれ、仮想分割パターンk121は、図8に示す2方向(k1211,k1212等)の対辺(好ましくは、対辺中央の中点)に分割点を設定して、それぞれ四角形形状要素を2つの小四角形に仮想分割し、表面が四角形メッシュからなる3次元形状要素の仮想分割辺の3次元長さのバランスをチェックするようになっており、具体的には、仮想分割パターンk1211では、図8(A)に示す上下の対辺に設定した分割点により、対辺を分割し、3次元長さで分割比ed1:ed2,ed3:ed4,edi:ed5(i=1〜4)等を求め、左右の3次元四角形形状要素の長さのバランスを評価するようになっている。又、仮想分割パターンk1212は、仮想分割パターンk1211を90度回転せしめて、上下方向に四角形メッシュを小さく分割し、図8(B)に示す左右の対辺に設定した分割点により、対辺を分割し、分割比ed6:ed7,ed8:ed9,edi:ed10(i=6〜9)等で、上下の3次元四角形形状要素の3次元長さのバランスを評価するようになっている。
更に、仮想分割パターンk122は、図9に示す4方向(k1221〜k1224等)の各分割パターンに対して、4つの各辺(好ましくは、各辺中央の中点)に、それぞれ四角形の分割点を設定して、各四角形形状要素を、上記中点及び対向する辺の頂点を連結して形成される3つの三角形形状要素に仮想分割し、四角形メッシュの仮想分割辺の長さのバランスをチェックするようになっており、具体的には、仮想分割パターンk1221では、図9(A)に示す上辺に設定した分割点により、1辺を分割し、分割比ed51:ed52,ed51:ed53,ed52:ed54等で、3つの三角形形状要素の長さのバランスを評価するようになっている。又、仮想分割パターンk1222は、仮想分割パターンk1221を90度時計回りに回転せしめて、右辺中央で四角形メッシュを3つの三角形形状要素に分割し、図9(B)に示す右辺に設定した分割点により、分割比ed55:ed56,ed55:ed57,ed56:ed58等で、3つの三角形形状要素の長さのバランスを評価するようになっている。又、仮想分割パターンk1223,k1224は、図9(A)に示す上辺に設定した分割点を更に時計回りに90度づつ回転させて形成したもので、1辺を中点で分割し、この中点と対辺の両端点とをそれぞれ連結して3つの三角形形状要素を仮想的に作成し、上記と同様な分割比演算により、3つの三角形形状要素の3次元長さのバランスを評価するようになっている。
更に又、仮想分割パターンk123は、図10に示す4方向(k1231〜k1234等)の分割パターンに対して、四角形形状要素の中央中心点及び1つの頂点を共有する2つの各辺(好ましくは、各辺中央の中点)に、それぞれ四角形の分割点を設定して、元の四角形形状要素を、それぞれ3つの小四角形形状要素に仮想分割し、小四角形メッシュの仮想分割辺の長さのバランスをチェックするようになっており、具体的には、仮想分割パターンk1231では、図10(A)に示す上辺及び右辺に設定した分割点と中央中点を連結して、3つの四角形形状要素に仮想分割し、かかる再分割された小四角形に生成された各辺を利用して、四角形メッシュの仮想分割辺の長さのバランスをチェックするようになっており、分割比ed11:ed12,ed13:ed14,ed15:ed17,ed16:ed17等で、3つの再分割四角形の長さのバランスを評価するようになっている。又、仮想分割パターンk1232は、仮想分割パターンk1231を90度時計回りに回転せしめて、四角形メッシュを3つの四角形に再分割したもので、図10(B)に示す右辺及び下辺に設定した分割点により、分割比ed21:ed22,ed23:ed24,ed25:ed27,ed26:ed27等で、3つの四角形の長さのバランスを評価するようになっている。又、仮想分割パターンk1233,k1234は、図10(A)の分割パターンを更に時計回りに90度づつ回転させて、それぞれ形成したものである。
尚、3次元形状モデル上の四角形要素メッシュの仮想分割辺の長さのバランスbg12が、ユーザが予め設定した許容誤差(最大g12/最小g13)の範囲外であれば、最もバランスの劣化した評価値を形成する仮想分割線に沿って、当該四角形メッシュを、再分割する。
上記仮想分割辺比のバランス値は、1.0が最適であり、この値から遠ざかるほど、解析精度が低下するが、例えば、バランスbg12が、最大g12=10以下、最小g13=0.1以上であるのが好ましく、バランスbg12が、最大g12=8以下、最小g13=0.15以上であるとより好ましい。
又、3次元形状モデル上の四角形要素メッシュの3次元アスペクト比のバランスは、図11に示すように、基本的には、各辺の中点と中央中心点とを連結して形成した1種類の仮想分割パターンk124により、4つの小さな仮想小四角形メッシュ(3次元)に再分割して、それぞれ、4つの再分割四角形メッシュに関して、アスペクト比ed71:ed73,ed72:ed79,ed71:ed76,ed76:ed75,ed75:ed73,ed76:3ed77,ed77:ed79等により、それぞれアスペクト比をチェックするようになっている。尚、3次元形状モデル上の四角形要素メッシュの仮想分割辺のアスペクト比が、ユーザが予め設定した許容誤差の範囲外であれば、図11に示す仮想分割線に沿って、当該四角形メッシュを、4つの四角形メッシュに再分割する。
上記仮想アスペクト比は、1.0が最適であり、この値から遠ざかるほど、解析精度が、低下し、例えば、アスペクト比が、最大10以下、最小0.1以上であるのが好ましく、アスペクト比が、最大8以下、最小0.15以上であるとより好ましい。
【0012】
k13) 四角形要素メッシュと、形状モデル表面との距離
次に、上記k11)ワープ歪及びk12)の長さのバランスに基づいた形状誤差評価基準では、十分に3次元形状モデル表面の曲面と要素メッシュの平面との形状誤差を近似できず、更に詳細な形状近似処理を実行したい場合には、図12(B)に示すように、実際の3次元形状モデル表面の曲面30gと、四角形要素メッシュms2との空間的3次元距離L6を、例えば、各要素メッシュの辺に沿って演算し、上記距離L6が所定の間隔以上に離れていた場合には、当該四角形メッシュを、より小さな面積の四角形メッシュで再分割するようにする。かかる再分割パターンとしては、基本的には、図13に示すような2種類(方向の違いを考慮すると、5種類)の仮想分割パターンk131、k132が利用可能であり、仮想分割パターンは、上記k12)における仮想分割パターンk122及び、k124と、それぞれ、同等であり、上記仮想分割パターンk131は、仮想分割する方向により、大きく変動するため、仮想分割パターンk131は、図13(A)〜(D)に示すような4方向(k1311〜k1314等)で、それぞれ仮想分割して、実際の3次元形状モデル表面の曲面30gと、四角形要素メッシュms2との空間的距離L6を評価し、又、仮想分割パターンk132は、図13(E)に示すような4つの正方格子状再分割四角形メッシュ・パターンで、それぞれ仮想分割して、実際の3次元形状モデル表面の曲面30gと、四角形要素メッシュms2との空間的距離L6を、合計5種類の3次元距離で評価し、形状誤差が少なくなるように、最大の形状誤差を示す再分割パターンで、該当する四角形メッシュを再分割するのが好ましい。
具体的には、図13(A)に示す仮想分割パターンk1311の場合では、辺ed84〜ed91に沿って、実際の3次元形状モデル表面の曲面30gと、3次元四角形要素メッシュms2との空間的3次元距離L6を評価し、3次元形状誤差を演算する。又、図13(E)に示す仮想分割パターンk132の場合では、辺ed92〜ed103に沿って、実際の3次元形状モデル表面の曲面30gと、四角形要素メッシュms2との空間的3次元距離L6を評価し、3次元近似形状誤差を演算する。
尚、上記空間的距離は、0.0が最適であり、この値から遠ざかるほど、解析精度が、低下し、例えば、要素メッシュ形状を構成する平面と形状モデル表面との再分割判定用距離g22は、モデルの最大長の1/50以内であるのが好ましく、モデルの最大長の1/100以内であるとより好ましく、モデルの最大長の1/200以内であると更に好ましい。
k14) 四角形要素形状の最大角度g18/最小角度g19
更に、四角形要素形状の形状誤差評価基準としては、個々の四角形要素の4隅の角度の大きさを利用することも可能であり、例えば、形状要素の4つの角度の大きさが、最大角度g18=130度〜最小角度g19=40度の範囲内であるのが好ましく、最大角度g18=120度〜最小角度g19=45度の範囲内であるとより好ましく、最大角度g18=115度〜最小角度g19=50度の範囲内であると更に好ましい。
又、上記4種類の形状誤差評価基準は、それぞれ相互に、包含関係を有するので、予め優先順位を設定して演算するのが好ましく、具体的には、k12)の長さのバランス > k11)ワープ歪 > k13)形状モデル表面との距離 > k14)4つの角度 の順番で、それぞれ形状誤差を評価し、ユーザが予め設定した許容誤差の範囲内で、順次、形状誤差が少なくなるように、必要であれば該当する四角形メッシュを、四角形メッシュ、及び/又は、三角形メッシュで再分割するのが好ましい。
【0013】
k2) 要素形状が三角形メッシュの場合
2次元三角形メッシュの投影により形成された3次元形状要素に対しては、予め登録した以下の複数種類の形状誤差評価パターンを用意し、それぞれ3次元形状誤差を評価し、ユーザが予め設定した許容誤差の範囲内で、総合的に、最も形状誤差が少なくなるように、必要であれば、当該2次元三角形メッシュを再分割して、3次元形状モデルに再投影する。
k21) 三角形メッシュの仮想分割辺の長さのバランス及びアスペクト比
3次元形状モデルへの投影により生成された3次元三角形要素メッシュの仮想分割辺の3次元長さのバランスは、図14(B)〜(H)に示すように、基本的には、3種類の仮想分割パターンk211〜k213により評価する。尚、上記仮想分割パターンは、仮想分割する方向により、大きく変動するため、それぞれ、仮想分割パターンk211は、図14(B)〜(D)に示す3方向(k2111〜k1213等)の1つの頂角を挟む辺(好ましくは、2つの挟辺のそれぞれ中点)に仮想分割点を設定して、それぞれ3次元三角形状要素を1つの三角形メッシュと1つの四角形メッシュとに仮想分割し、三角形及び四角形メッシュの仮想分割辺の3次元長さのバランスをチェックするようになっており、具体的には、仮想分割パターンk211では、図14(B)に示す左右の斜辺に設定した分割点により、当該斜辺を分割し、分割比ed110:ed111,ed110:ed112,ed114:ed115,ed112:ed114,ed112:ed113等で、上下の3次元三角形形状要素及び四角形形状要素の3次元長さのバランスを評価するようになっている。又、仮想分割パターンk212は、仮想分割パターンk211を時計周りに90度回転せしめて、右下方向に三角形メッシュを分割し、図14(C)に示す挟辺に設定した分割点により、挟辺を分割し、複数の分割比で、3次元三角形/四角形形状要素の3次元長さのバランスを評価するようになっている。
尚、3次元形状モデル上の三角形要素メッシュの仮想分割辺の3次元長さのバランスが、ユーザが予め設定した許容誤差の範囲外であれば、最もバランスの劣化した評価値を形成する仮想分割線に沿って、当該四角形メッシュを、再分割する。上記仮想分割辺比のバランス値は、1.0が最適であり、この値から遠ざかるほど、解析精度が、低下するが、例えば、バランスbg14が、最大g14=8以下、最小g15=0.15以上であるのが好ましく、バランスbg14が、最大g12=7以下、最小g13=0.2以上であるとより好ましい。
k22) 又、3次元形状モデル上の三角形要素メッシュのアスペクト比のバランスは、図14(E)〜(G)に示すように、基本的には、各辺の中点と当該辺に対応する対角頂点とを連結して形成した3種類の仮想分割パターンk221〜k223により、3種類の小さな仮想3次元三角形メッシュに再分割して、それぞれ、3つの再分割三角形メッシュに関して、分割パターンk221の場合には、アスペクト比ed120:ed121,ed120:ed122,ed121:ed122等により、それぞれアスペクト比をチェックし、分割パターンk222及びk223は、分割パターンk221を、反時計周りに、順次、90度づつ回転させて、アスペクト比を演算する方向を変化せしめ、3次元形状モデル表面と3次元形状要素との形状歪を評価するようになっている。
尚、3次元形状モデル上の三角形要素メッシュの仮想分割辺のアスペクト比が、ユーザが予め設定した許容誤差の範囲外であれば、図14に示す仮想分割線に沿って、当該三角形メッシュを、2つの三角形メッシュに再分割する。上記仮想アスペクト比は、1.0が最適であり、この値から遠ざかるほど、解析精度が低下するが、好ましいアスペクト比の範囲は、バランスbg14と同様である。
【0014】
k23) 三角形要素メッシュと、形状モデル表面との距離
次に、上記k21)分割辺のバランス比及びk22)のアスペクト比に基づいた3次元形状誤差評価基準では、十分に3次元モデル表面の曲面と要素メッシュの平面との3次元形状誤差を近似できず、更に詳細な3次元形状近似処理を実行したい場合には、図14(H)に示す分割パターンk23のように、三角形の重心と各辺の中点を連結して形成した3つの四角形メッシュの各分割辺に沿って、曲面と平面との3次元距離を演算評価し、実際の3次元形状モデル表面の曲面30gと、三角形要素メッシュとの空間的3次元距離L6を、例えば、各三角形要素メッシュの辺ed131〜ed139に沿って演算し、上記3次元距離L6が所定の間隔以上に離れて得いた場合には、当該三角形メッシュを、図14(H)に示すような、より小さな面積の3つの四角形メッシュで再分割するようにする。尚、上記空間的距離は、0.0が最適であり、この値から遠ざかるほど、解析精度が、低下し、例えば、要素メッシュ形状を構成する平面と形状モデル表面との再分割判定用距離は、モデルの最大長の1/50以内であるのが好ましく、モデルの最大長の1/100以内であるとより好ましく、モデルの最大長の1/200以内であると更に好ましい。
k24) 三角形要素形状の最大角度g20/最小角度g21
更に、三角形要素形状の形状誤差評価基準としては、個々の三角形要素の3隅の角度の大きさを利用することも可能であり、例えば、形状要素の3つの角度の大きさが、最大角度g20=170度〜最小角度g21=20度の範囲内であるのが好ましく、最大角度g20=160度〜最小角度g21=25度の範囲内であるとより好ましく、最大角度g20=155度〜最小角度g21=30度の範囲内であると更に好ましい。
又、上記4種類の3次元形状誤差評価基準は、それぞれ相互に、包含関係を有するので、予め優先順位を設定して演算するのが好ましく、具体的には、k21)長さのバランス > k22)アスペクト比 > k23)形状モデル表面との距離 > k24)3つの角度の大きさ の順番で、それぞれ形状誤差を評価し、ユーザが予め設定した許容誤差の範囲内で、順次、形状誤差が少なくなるように、必要であれば該当する三角形メッシュを、四角形メッシュ及び/又は三角形メッシュで再分割するのが好ましい。
【0015】
e3) 要素歪みが要素形状基準値以下であるか否か判定し(S10)、要素形状基準値以上の3次元形状歪が大きい要素が検出されると、かかる形状要素を、投影基準面の2次元写像モデルに戻って、形状歪が最大の分割モデルで再分割し(S12)、上記ステップe1)に戻る。則ち、再分割メッシュに所定の範囲外の形状歪みの大きい要素が無くなる迄、上記の処理を1回又は複数回繰返す。
例えば、図5(A)に示すような仮想平行光線によって生成された図5(C)の最外郭輪郭線の周辺の要素形状ms32aの場合、表面が三角形メッシュ又は四角形メッシュいずれの場合にも、最初の各辺のバランス比の範囲をチェックすると、いずれも、所定の範囲外なので、直ちに、再分割処理を行うと、図5(D)に示す要素形状ms32a1,ms32a2のような高精度の3次元メッシュが、高速に生成可能である。
尚、本発明では、再分割処理をする場合、初期メッシュとして生成した2次元写像モデルの各格子点は、再分割処理の時点では、新しく最適な格子点位置を求めるために移動/変動させる操作は最小限ですみ、更に、初期メッシュの各格子点に再分割用の2次元格子点を、所定の分割パターンに従って追加するだけで再分割可能なので、再分割の演算時間が非常に高速化可能である。
e4) 全ての3次元要素歪みが、要素形状基準値以下であると判断されたとき、生成された3次元メッシュを、所定の投影基準面に対する第1の最終3次元メッシュとして記憶する(S14)。
e5) 複数の投影基準面が設定されている場合には、上記e1)〜e4)の処理を、全ての投影基準面に対して行う(S16)。
c3) 投影基準面が複数ある場合には、上記e1)〜e5)の処理により生成された各投影基準面に対する第1の最終3次元メッシュを、上記d1),d2)で作成した境界線に基づいて合成し、第2の最終3次元メッシュとして記憶する(S18)。
【0016】
【発明の効果】
本発明によれば、3次元測定機等から入力した3次元計測オリジナル・データ点群からなる不規則/複雑な形状でも、計算精度を十分保った状態で、解析メッシュ要素の作成が実用に耐える高速な速度で可能であり、3次元計測オリジナル・データ点群で規定される面の形状・大きさ等面の品質に左右されることなく、品質の良いメッシュを生成し、精度の高い数値計算用のメッシュモデルを作成する方法を提供することができ、ユーザが設定した要素形状基準値を超える形状歪みの大きい要素は、効率的に検出され、該要素に対応する2次元写像モデルを修正し、形状歪みの大きい要素を再分割し修正する処理を高速かつ効率的に行うことができるので、従来の曲線座標変換法のような要素の形状歪みを低減させるための3次元六面体モデルの変更作業やメッシュ作成後の3次元メッシュ修正作業が不要となり、1〜20分以内で、大抵の3次元形状モデルに対して、3次元メッシュを生成可能である。従って、3次元形状モデルの修正、変更が容易となり、実際の設計変更に柔軟に対応し、尚かつ、有限要素解析結果を、直ちに、3次元形状モデルにフィードバックでき、精度が高く、信頼性の高い解析モデルのメッシュ作成が実現できる。
又、本発明のメッシュ生成法では、形状データの面の形状・大きさ・隣り合った面の間隔等、3次元形状モデルの面の品質に影響されずに品質の良いメッシュを自動生成可能であり、六面体メッシュを中心に精度の高い数値計算用の3次元メッシュモデルを作成することができる。更に、投影先を複数にすることで、より複雑な3次元形状モデルに対してもメッシュを生成することができる。
【図面の簡単な説明】
【図1】本発明のシステム全体の構成の1例を示すブロック図である。
【図2】本発明のメッシュ生成工程の1例を示すフローチャートである。
【図3】本発明の2次元初期メッシュ生成工程の概念を示す図である。
【図4】本発明の正方格子式初期メッシュ生成工程の1例を示す図である。
【図5】本発明のシェル・メッシュ生成工程の概念を説明する図である。
【図6】本発明のメッシュ生成工程において、格子点の投影先に形状モデルが無い場合に、新しく格子点の座標を演算する過程を説明する図である。
【図7】本発明の四角形形状要素のワープ歪の大きさを演算する例を説明する図である。
【図8】本発明の四角形形状要素の分割辺の長さのバランスを、2つの分割四角形メッシュから演算する例を示す図である。
【図9】本発明の四角形形状要素の分割辺の長さのバランスを、3つの分割三角形メッシュから演算する例を示す図である。
【図10】本発明の四角形形状要素の分割辺の長さのバランスを、3つの分割四角形メッシュから演算する例を示す図である。
【図11】本発明の四角形形状要素のアスペクト比を、4つの分割四角形メッシュから演算する例を示す図である。
【図12】本発明の形状モデルと3次元形状要素との近似誤差を、3次元距離により演算する例を示す図である。
【図13】本発明の形状モデルと3次元四角形形状要素との近似誤差を、3つの分割四角形メッシュの辺に沿って演算した3次元距離により評価する例を示す図である。
【図14】本発明の形状モデルと3次元三角形形状要素との近似誤差を、分割辺の長さのバランス及びアスペクト比により評価する例を示す図である。
【符号の説明】
2 形状入力部
4 初期メッシュ生成部(写像モデル生成部)
6 3次元メッシュ生成部
8 要素形状評価部
10 写像モデル修正部
12 メッシュ出力部
14 メッシュ統合部
Claims (3)
- 解析対象の3次元形状モデルと解析メッシュ生成用基礎データとを入力する形状入力部と、写像モデルの格子点を3次元形状モデルに写像して解析メッシュを生成し、生成された解析メッシュをデータベースに格納するメッシュ生成部と、を含む解析メッシュ生成装置におけるメッシュ生成方法において、
前記解析メッシュ生成装置が、
前記形状入力部から入力した前記3次元形状モデルを、所定の投影方向に対応して予め設定されている投影基準面に対して当該投影方向から投影し、当該基準面に前記3次元形状モデルの最外郭線の内部からなる基準面モデルを生成する処理を実行する第1工程と、前記基準面モデルの内部に、所定のドット生成規則により生成された格子点からなる2次元写像モデルを生成する処理を実行する第2工程と、を含み、
前記解析メッシュ生成装置が、
前記2次元写像モデルの格子点を、前記投影基準面から、仮想平行光線により前記3次元形状モデルに投影して解析メッシュを生成する処理を実行する第3工程と、
前記解析メッシュを後記分割パターンにより仮想的に分割したときの分割メッシュを構成する各仮想分割辺にて隣り合う2辺の比が最大の比をアスペクト比とし、前記分割メッシュを構成する平面と前記3次元形状モデルの曲面との間隔の最大値を3次元形状の近似誤差として、前記分割メッシュに対して、前記アスペクト比、前記仮想分割辺の長さのバランス、ワープ歪、及び、要素形状の最大角度と最小角度がそれぞれ所定の範囲内で、かつ、前記3次元形状の近似誤差が所定の範囲内であることを条件として、予め用意した四辺形メッシュ及び/又は三角形メッシュの分割パターンに基づいて、前記条件を満たしていない要素メッシュを分割する処理を実行する第4工程と、
前記第3工程での前記2次元写像モデルの格子点の処理対象を前記第4工程による分割後の要素メッシュの格子点として前記第3工程に戻り、前記条件を満たすまで前記第3工程と前記第4工程とを繰り返し実行し、前記条件を満たしたと判定したときに、前記第3工程により生成された解析メッシュを前記投影基準面に対する第1の最終3次元メッシュとして前記データベースに格納すると共に、前記所定の投影方向として複数の投影方向が設定されている場合には、前記第1工程から前記第4工程の処理によって前記第1の最終3次元メッシュを前記複数の投影方向毎に生成し、生成された第1の最終3次元メッシュ同士を相互に連結して1つの解析メッシュに統合する処理を実行し、第2の最終3次元メッシュとして前記データベースに格納する第5工程と、
を更に含むことを特徴とする3次元メッシュ生成方法。 - 前記3次元形状のオリジナル・データが、3次元計測オリジナル・データ、又は、3次元ソリッドモデルCAD形状オリジナル・データである請求項1に記載の3次元メッシュ生成方法。
- 解析対象の3次元形状モデルデータと解析メッシュの目標要素寸法とを入力する手段と、該形状モデルデータと該目標要素寸法とから解析メッシュ生成のための2次元写像モデルを生成する手段と、該2次元写像モデルの格子点を3次元形状モデルに写像する手段と、から構成される解析メッシュ生成装置において、
前記3次元形状モデルを、所定の投影方向に対応して予め設定されている投影基準面に対して当該投影方向から投影し、当該基準面に前記3次元形状モデルの最外郭線の内部からなる基準面モデルを生成し、この基準面モデルの内部に、所定のドット生成規則により生成された格子点からなる2次元写像モデルを生成する手段と、前記2次元写像モデルの格子点を、前記投影基準面から、仮想平行光線により前記3次元形状モデルに投影して解析メッシュを生成するメッシュ生成手段と、
前記解析メッシュを後記分割パターンにより仮想的に分割したときの分割メッシュを構成する各仮想分割辺にて隣り合う2辺の比が最大の比をアスペクト比とし、前記分割メッシュを構成する平面と前記3次元形状モデルの曲面との間隔の最大値を3次元形状の近似誤差として、前記分割メッシュに対して、前記アスペクト比、前記仮想分割辺の長さのバランス、ワープ歪、及び、要素形状の最大角度と最小角度がそれぞれ所定の範囲内で、かつ、前記3次元形状の近似誤差が所定の範囲内であることを条件として、予め用意した四辺形メッシュ及び/又は三角形メッシュの分割パターンに基づいて、前記条件を満たしていない要素メッシュを分割する分割手段と、
前記メッシュ生成手段での前記2次元写像モデルの格子点の処理対象を前記分割手段による分割後の要素メッシュの格子点として、前記条件を満たすまで、前記メッシュ生成手段による前記解析メッシュの生成と前記分割手段による前記要素メッシュの分割とを繰り返し実行し、前記条件を満たしたと判定したときに、前記第3工程により生成された解析メッシュを前記投影基準面に対する第1の最終3次元メッシュとしてデータベースに格納すると共に、前記所定の投影方向として複数の投影方向が設定されている場合には、前記第1工程から前記第4工程の処理によって前記第1の最終3次元メッシュを前記複数の投影方向毎に生成し、生成された第1の最終3次元メッシュ同士を相互に連結して1つの解析メッシュに統合する処理を実行し、第2の最終3次元メッシュとして前記データベースに格納する手段と、
を有することを特徴とする解析メッシュ生成装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002168538A JP4639292B2 (ja) | 2002-06-10 | 2002-06-10 | 3次元メッシュ生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002168538A JP4639292B2 (ja) | 2002-06-10 | 2002-06-10 | 3次元メッシュ生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004013672A JP2004013672A (ja) | 2004-01-15 |
JP4639292B2 true JP4639292B2 (ja) | 2011-02-23 |
Family
ID=30435421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002168538A Expired - Lifetime JP4639292B2 (ja) | 2002-06-10 | 2002-06-10 | 3次元メッシュ生成方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4639292B2 (ja) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006318275A (ja) * | 2005-05-13 | 2006-11-24 | Toyota Motor Corp | 数値解析用メッシュ作成装置及び方法 |
JP4592101B2 (ja) * | 2006-11-20 | 2010-12-01 | 関東自動車工業株式会社 | 有限要素解析のメッシュデータ自動変換装置 |
JP5561975B2 (ja) * | 2009-03-09 | 2014-07-30 | カルソニックカンセイ株式会社 | 表面加工データの作成方法および装置 |
JP5227436B2 (ja) * | 2011-03-18 | 2013-07-03 | 住友ゴム工業株式会社 | フィラー配合ゴムの有限要素モデルの作成方法 |
JP5834317B2 (ja) | 2011-05-11 | 2015-12-16 | 株式会社ミツトヨ | 誤差画像生成方法 |
CN104525668B (zh) * | 2014-11-12 | 2016-06-08 | 常州工学院 | 一种多点成形加工曲面的冲头调形方法 |
CN108228962B (zh) * | 2017-11-27 | 2024-04-26 | 中建华东投资有限公司 | 一种基于bim的基坑周边墙体测斜热力图生成方法 |
CN110117640A (zh) * | 2019-05-15 | 2019-08-13 | 重庆市公安局刑事警察总队 | 一种物证的网格标记方法及装置 |
CN110135067B (zh) * | 2019-05-16 | 2023-06-16 | 杭州电子科技大学 | 一种双时间步方法下的直升机流场重叠混合网格并行方法 |
JP7324726B2 (ja) * | 2020-03-03 | 2023-08-10 | 大成建設株式会社 | メッシュモデル生成装置及びメッシュモデル生成方法 |
CN112528541A (zh) * | 2020-12-09 | 2021-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种基于节点重构的双曲率座舱盖有限元建模方法 |
CN116597109B (zh) * | 2023-01-13 | 2023-12-22 | 东莘电磁科技(成都)有限公司 | 一种复杂三维曲面共型网格生成方法 |
-
2002
- 2002-06-10 JP JP2002168538A patent/JP4639292B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2004013672A (ja) | 2004-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7212205B2 (en) | Curved surface image processing apparatus and curved surface image processing method | |
JP4639292B2 (ja) | 3次元メッシュ生成方法 | |
JP4185698B2 (ja) | メッシュ生成方法 | |
Marchandise et al. | CAD and mesh repair with radial basis functions | |
WO2007083602A1 (ja) | 補間処理方法、補間処理装置、形状評価方法、および形状評価装置 | |
Bronson et al. | Lattice cleaving: A multimaterial tetrahedral meshing algorithm with guarantees | |
US20240153123A1 (en) | Isogeometric Analysis Method Based on a Geometric Reconstruction Model | |
CN116229015B (zh) | 一种基于2N-Tree带附面层的贴体笛卡尔网格生成方法 | |
US7348976B2 (en) | Three-dimensional shape display program, three-dimensional shape display method, and three-dimensional shape display apparatus | |
Wang et al. | Feature-sensitive tetrahedral mesh generation with guaranteed quality | |
JP3265879B2 (ja) | 3次元直交格子データの生成装置 | |
US6914601B2 (en) | Method, apparatus, and computer program for generating three-dimensional shape data or volume data | |
CN115758938A (zh) | 面向粘性边界流场数值模拟的附面层网格生成方法 | |
CN106649992A (zh) | 舰船与尾迹的网格模型的融合与优化方法 | |
CN110796735B (zh) | Nurbs曲面有限元板壳网格划分方法及计算机实现系统 | |
JPH0636013A (ja) | 地形データの作成方法および装置 | |
Xu et al. | Hexahedral meshing with varying element sizes | |
CN115272379A (zh) | 一种基于投影的三维网格模型外轮廓提取方法及系统 | |
CN111797551B (zh) | 一种面向试验应变数据三维云图显示的插值方法 | |
JP2007156604A (ja) | 曲面作成方法及び曲面作成プログラム並びに3次元形状処理装置 | |
JP3463843B2 (ja) | 自由曲面生成装置および自由曲面生成方法 | |
JPH11195139A (ja) | 自由曲面生成装置および方法、並びに提供媒体 | |
JP2655056B2 (ja) | テクスチャデータ生成装置 | |
Tse et al. | Tin meets CAD-extending the TIN concept in GIS | |
JP2003044528A (ja) | 物体の表面格子生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050603 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20050603 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20080303 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20080303 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080715 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080916 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20081111 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090113 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20090604 |
|
A912 | Re-examination (zenchi) completed and case transferred to appeal board |
Free format text: JAPANESE INTERMEDIATE CODE: A912 Effective date: 20090703 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20101020 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4639292 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131210 Year of fee payment: 3 |
|
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 |
|
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 |
|
EXPY | Cancellation because of completion of term |