JPH04309183A - 有限要素メッシュ生成方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、有限要素解析を使用可
能にする方法に関し、より詳細には、すべての有限要素
が全体的に解析される体の内部、又は外部に存在するの
を確保するために、生成されたメッシュを自動的にトリ
ミングする方法に関する。

【０００２】

【従来の技術】自動モデル化及び設計手順には、体の幾
何学的モデルの初期展開、及び当該モデルと様々な境界
条件、材料の特性、及び後続の構造解析における許容誤
差の指示との関連性が含まれる。この解析は、モデルを
後で個別に解析される複数の有限要素に分割するメッシ
ュを生成する適応性のある方法を利用することにより、
自動的に実行される。

【０００３】メッシュの生成は、完全に自動的に実行さ
れることが理想的であり、即ち、入力として、関連する
メッシュ制御情報と一緒にモデルの幾何学的な表示を取
り、次に自動的に有効なメッシュを生成するような方法
である。メッシュを生成させる一つのアプローチとして
は、ドローネー（Ｄｅｌａｕｎａｙ）三角形化（後で引
用されるフレイ（Ｆｒｅｙ）、及びガベンディシュ（Ｃ
ａｖｅｎｄｉｓｈ）の論文を参照のこと）に依存するも
のである。ドローネー方法は、三次元の場合に一般的に
一組の四面体要素を意味する有限要素メッシュを生成す
るために、モデル空間内に点集合を三角形化することに
基づいている。そこで、ドローネー方法は、初期におい
て、点集合を生成し、次に、四面体を含む外接球がその
体積内にメッシュ点を全く含まないという特性を満たす
ような点の三角形化を実行する。

【０００４】このようなメッシュを生成する上での問題
点の一つに、有限要素は完全にモデルの境界の内部、ま
たは外部に存在するのを確認することがあり、従って、
境界外にある有限要素は後続の解析では分離、及び無視
される可能性が出てくる。

【０００５】図１は、実質的にＣ形を有するモデルの一
例を示す斜視図である。このモデルは、複数の頂点によ
り定義され、モデルの各面を定義している頂点の集合（
例えば、頂点１０、１２、１４、１６、１８、２０、２
２、及び２４はモデルの面２６を定義している）を有す
る。周知のように、各頂点は、三次元空間における位置
を定義する一組の座標として入力される。

【０００６】図１のモデル上で有限要素解析を実行する
には、有名なドローネー・メッシュ生成プログラムを用
いて、ドローネー四面体化がモデルの各頂点と各面に実
行されることが必要である。図２では、三次元において
ドローネー方法の基本出力を形成する四面体が示されて
いる。図３では、ドローネーのモデル四面体化に続いて
、面２６の平面図が示されている。ドローネー手順を実
行中は、形成された四面体がモデルの境界を交叉するエ
ッジを表す事象が頻繁に発生する。このような状況の下
では、「全体が内部／外部にある」の規格を満たすため
に、そのような四面体を見つけ、細分化する必要がる。 図３に示されたような場合には、頂点１２、１４、２０
によって、及び１２、１６、２０によって定義された四
面体は頂点１４と１６の間の境界を不明瞭にする。 この例外的な事象は、有限要素解析を実行する前に訂正
されなくてはならない。

【０００７】図４の（ａ）と（ｂ）では、ドローネー四
面体化の初期の段階で発生するかもしれない共通する「
違反」を幾つか示している。図４（ａ）では、オブジェ
クトのエッジは、三角形の頂点以外の点においてドロー
ネー三角形の内部を通過するように示されている。図４
（ｂ）では、ドローネー三角形のエッジは、エッジを定
義している頂点上の面Ｆの頂点以外の点において、モデ
ルの面Ｆを通過するように示されている。これらの二つ
の状況は，「全体が内部／外部にある」の規格に違反す
る四面体を生じさせる。

【０００８】従来技術では、様々な方法で上記定義の違
反を取り扱おうとしてきた。全体が内部／外部の規格を
達成するために用いられた先の方法のうちの一つは、モ
デル上でドローネー手順を実行させ、ユーザに対して、
モデルを調べさせ、次に、違反している四面体を取り除
くために、メッシュの中に追加頂点を挿入させることで
あった。幾つかの従来技術の引用例では、初期の点メッ
シュを「もっと機敏に」生成することによって、この問
題を回避しようとしている。例えば、次のような引用例
を参照しよう。１９８４年の数値方法と工学インターナ
ショナルジャーナル（Ｉｎｔｅｒｎａｔｉｏｎａｌ　　
Ｊｏｕｒｎａｌ　　Ｆｏｒ　　Ｎｕｍｅｒｉｃａｌ　　
Ｍｅｔｈｏｄｓ　　ａｎｄ　　Ｅｎｇｉｎｅｅｒｉｎｇ
）のコラム（Ｃｏｌｕｍｎ），第２０巻（Ｖｏｌ．２０
）」の１９６５ページから１９９０ページに発表された
エリー（Ｙｅｒｒｙ）らによる「修正されたオクトリー
（Ｏｃｔｒｅｅ）技術を用いた自動三次元メッシュ生成
」、１９８５年の同誌第２１巻（Ｖｏｌ．２１）の３２
９ページから３４７ページに発表されたガベンディシュ
（Ｃａｖｅｎｄｉｓｈ）らによる「自動三次元有限要素
メッシュ生成へのアプローチ」、１９８３年１１月発行
の磁気学ＩＥＥＥ会報（ＩＥＥＥ　　Ｔｒａｎｓａｃｔ
ｉｏｎｓ　　ｏｎ　　Ｍａｇｎｅｔｉｃｓ），第Ｍａｇ
−１９巻（Ｖｏｌ．Ｍａｇ−１９），Ｎｏ．６の２５５
２ページから２５５４ページに発表されたセンデス（Ｃ
ｅｎｄｅｓ）らによる「ドローネー三角形化及び相補有
限要素方法を用いた磁界計算」、更に、１９８７年の数
値方法と工学インターナショナルジャーナル（Ｉｎｔｅ
ｒｎａｔｉｏｎａｌ　　Ｊｏｕｒｎａｌ　　Ｆｏｒ　　
Ｎｕｍｅｒｉｃａｌ　　Ｍｅｔｈｏｄｓ　　ａｎｄ　　
Ｅｎｇｉｎｅｅｒｉｎｇ）のコラム（Ｃｏｌｕｍｎ），
第２４巻（Ｖｏｌ．２４）」の２１８３ページから２２
００ページに発表されたフレイ（Ｆｒｅｙ）による「選
択的精製：段三角形メッシュにおける自動ノード配置の
ための新しい戦略」等の論文がある。

【０００９】１９８８年の数値方法と工学インターナシ
ョナルジャーナル（Ｉｎｔｅｒｎａｔｉｏｎａｌ　　Ｊ
ｏｕｒｎａｌ　　ｏｆ　　Ｎｕｍｅｒｉｃａｌ　　Ｍｅ
ｔｈｏｄｓａｎｄ　　Ｅｎｇｉｎｅｅｒｉｎｇ），第２
６巻（Ｖｏｌ．２６）」の２５０３ページから２５１５
ページにフィールド（Ｆｉｅｌｄ）によって発表された
「三次元におけるワトソン・アルゴリズムの実行」と題
する論文では、立体の境界付近に四面体を有し、更にそ
の境界上に四面体の一部分が延びるという問題点が認識
される。但し、著者は、そのような四面体の細分化を可
能するのに追加点が必要であることを指示しているだけ
であって、その追加される点を挿入する方法については
述べていない。

【００１０】シュローダー（Ｓｃｈｒｏｅｄｅｒ）らに
よって工学における数値方法インターナショナルジャー
ナル（Ｉｎｔｅｒｎａｔｉｏｎａｌ　　Ｊｏｕｒｎａｌ
　　ｆｏｒ　　Ｎｕｍｅｒｉｃａｌ　　Ｍｅｔｈｏｄｓ
　　ｉｎ　　Ｅｎｇｉｎｅｅｒｉｎｇ），第２６巻（Ｖ
ｏｌ．２６）の２５０３ページから２５１５ページに発
表された「幾何学に基づく完全自動化メッシュ生成及び
ドローネー三角形化」と題する論文において、ドローネ
ー生成四面体のトリミングを可能にする解決法が提示さ
れている。シュローダー（Ｓｃｈｒｏｅｄｅｒ）らは、
追加される点がメッシュ要素の境界とモデルの表面との
正確な交叉点に配置されるべきあると述べている。この
解決法は、「全体が内部／外部にある」の規格を達成す
る一方、数多くの点がメッシュに追加される方向に導く
のを可能にする。これは、必ずしもそのような詳細な細
分化を必要としないモデル点に多くの異なる四面体を生
成させる。 その結果、分析される有限要素の数が多くなるために、
後続の有限要素分析を実行している最中の計算の速度が
実質的に遅くなる。

【００１１】メッシュ生成とその改善に多少関係がある
他の従来技術が、米国特許Ｎｏ．４，９１２，６６４、
米国特許Ｎｏ．４，８８８，７１３，米国特許Ｎｏ．４
，６９７，１７８において、更に、１９８９年６月のＩ
ＢＭ技術公開公報第３２巻（Ｖｏｌ．３２），Ｎｏ．１
の３４０ページから３４２ページに記載されたコヤマダ
（Ｋｏｙａｍａｄａ）の記事、及び、１９７５年９月の
同報第１８巻（Ｖｏｌ．１８），Ｎｏ．１の１１６３ペ
ージから１１７５ページに記載されたシュレイバー（Ｓ
ｃｈｒｅｉｂｅｒ）による記事に見つけることができる
。

【００１２】

【発明が解決しようとする課題】本発明の目的は、生成
された四面体に対して、「全体が内部／外部にある」の
規格が自動的に満たされるように改善された有限要素生
成システムを提供することである。

【００１３】更にもう一つの目的としては、前記規格を
満たすように生成された追加メッシュ点の数が最小化さ
れる三次元オブジェクトのための自動メッシュ生成シス
テムを提供することがある。

【００１４】

【課題を解決するための手段及び作用】解析されるべき
エッジや面を表わす体の内部にある有限要素のメッシュ
を生成する方法が記載されている。この方法は、体に対
して、各々が複数の頂点とエッジを有する要素から成る
初期メッシュを生成するステップと、有限要素エッジと
一致しない体エッジ・セグメントを見つけるステップと
、該体エッジ・セグメント上に、必要とされる新しい頂
点の数を最小限にするように、頂点を一つ追加するステ
ップとを含む。次に、新しい頂点を隣接する頂点に接続
することにより、追加有限要素が作成される。体のエッ
ジ・セグメントと一致するエッジと体の一つの面を通過
する別のエッジ・セグメントとの両方を有する有限要素
が識別される。交叉している要素エッジの付近にある体
表面上に、追加されるべき頂点の数を最小限にするよう
に、新しい頂点が配置され、次に、本発明では、追加さ
れた頂点を隣接する頂点に接続することにより、追加の
有限要素を生成する。

【００１５】

【実施例】上述したように、立体のオブジェクトのため
のドローネー・メッシュ生成手順は、四面体メッシュを
生成する。このメッシュ生成の間、ドローネー手順では
、四面体毎に外接球を計算する。生成されたメッシュは
、「全体が内部／外部にある」の規格を満たすかもしれ
ないし、又、満たさないかもしれない。後で明白になる
ように、既に計算された外接球は、頂点の数を最小限に
し、更に四面体化構造を改善するために、本発明によっ
て利用されている。

【００１６】さて、図５を考察すると、トリミング手順
の全体の高度なフロー図が示されている。本質的に、ト
リミング手順では、すべての四面体は、「全体が内部／
外部にある」の規格を満たすことを確保するものである
。記載される手順には、汎用コンピュータの操作を制御
する複数のアルゴリズムが含まれる。このようなアルゴ
リズムは、図形表示機能を備えているワークステーショ
ン・レベルのシステム上で実行される。

【００１７】一般的なトリミング手順に対する入力は、
分析されるオブジェクト、及び、ドローネー四面体化ア
ルゴリズム（ボックス３０）によって生成された、既に
完成している有限要素メッシュのための境界定義である
。オブジェクトの面Ｆ毎に、該面のエッジは、エッジ・
トリミング手順（ボックス３２）を呼び出すことによっ
て復元される。このサブルーチンは、いずれかのオブジ
ェクト・エッジ、即ち、四面体エッジと完全には一致し
ないいずれかのセグメントに対して、該オブジェクト・
エッジ・セグメントが後に四面体エッジと一致するのを
可能にする追加された頂点が提供されることを確保する
。

【００１８】境界上のすべてのエッジ・セグメントは、
一旦復元されると、更に処理を進めるためにリストＬに
配置される（ボックス３４）。この機能がすべてのオブ
ジェクト面に対して実行された後、リストＬの各オブジ
ェクト・エッジ・セグメントは、それが該オブジェクト
の面と交叉しているか否かを判定するために更に先の処
理を実行され、仮に交叉する場合には、そのような交叉
しているセグメントを一部分とする四面体を細分化する
ために、頂点が１個追加される。

【００１９】先ず初めに、セグメントＡＢがリストＬか
ら除去され（ボックス３６）、面トリミング手順が呼び
出される。これにより、頂点Ｘが見つけだされ、三角形
ＡＢＸが面の一部分として復元される（ボックス３８）
。頂点Ｘが新しい場合は、エッジ・トリミング手順が反
復される（ボックス４０）。一方、頂点Ｘが新しくなく
て、更に、セグメントＡＸとＢＸがリストＬに記載され
ていない場合には、これらのセグメントはリストＬに追
加される（ボックス４２）。

【００２０】次に図６を考察すると、図７から図１１ま
での図と関連させて、エッジ・トリミング手順の詳細が
記述されている。図６の流れ図を考察する前に、エッジ
・セグメントＡＢが図示され、それがモデル境界の一部
分を構成する図７を参照することにしよう。モデルは既
にドローネー四面体化されているので、複数の四面体が
セグメントＡＢの周囲に配置されている。但し、セグメ
ントＡＢは四面体のいずれの境界とも一致しないことは
判定できる。従って、「全体が内部／外部にある」の規
格に対して違反が存在する。

【００２１】図８は、８−８線に沿って切断され、四面
体境界ＡＥを中心に置いて描かれた、図７の四面体構造
を示す図である。この図は、モデル線分セグメントＡＢ
を囲んでいる四面体がどの様にして見つけられるのかを
図示するために、図９においても利用されている。探索
手順は、図６の流れ図において図示されている。

【００２２】ボックス５０で示されるように、エッジ・
セグメントＡＢの点Ａに頂点を１個持っている四面体が
識別される。（この探索は、四面体メッシュの初期の生
成時において、各頂点毎に入射四面体が格納されていた
場合には、単純化される。）一旦入射四面体が識別され
ると、エッジ・セグメントＡＢが四面体Ｔの基底を通過
する四面体Ｔを見つけるために、頂点Ａを共有する四面
体の間で探索処理が実行される（ボックス５２）。

【００２３】この探索は、入射四面体（即ち、四面体Ａ
ＥＩＪ）の各面を平面として考えることによって、図９
に示されるように完了する。四面体の内部に最も接近し
ている平面の側面はプラス側として分類され、反対側の
側面はマイナス側として分類される。図９に見ることが
できるように、平面８０は面ＡＩＪと一致し、平面８２
は面ＡＩＥと一致し、更に、平面８４は面ＡＪＥと一致
する。頂点Ｅから開始して、左回りに処理を進めていく
と、平面毎の等式が見つけられ、頂点Ｂが個々の平面の
プラスまたはマイナス側にあるか否かが判定される。頂
点Ｂが平面のいずれかのマイナス側にある場合は、手順
は、該平面を共有する四面体に進む。図９では、上記規
格を満たしている唯一の平面は、平面８４であり、四面
体ＡＥＪＣが次に考察される四面体であることが判る。 該四面体の解析は、該四面体が（頂点Ｂに対応してマイ
ナス値を有する）四面体平面ＡＥＣと一致する平面であ
るので、四面体ＡＥＣＤに導いていく。

【００２４】一つ以上の平面が頂点Ｂに対応して「負の
側に」配置されていることが判った場合には、点Ｂに対
して最小の絶対変位値を有する平面が選択される。この
探索は、頂点Ｂがすべての平面構造の正の側にある時に
終了する。これは、四面体ＡＥＣＤが見つけられた時に
発生することに注意されたい。頂点Ｂ（中心に置いて）
の優位な点から図９を眺めてみると、上述した平面関係
がはっきりと眼に見えてくる。

【００２５】図６に戻って考察すると、四面体Ｔが一旦
見つけられると、エッジ・セグメントＡＢが四面体Ｔの
基底の頂点Ｘを通過しているか否かが判定される。通過
している場合には、該エッジ・セグメント（即ち、ＡＸ
）が格納される（ボックス５６）。（ＡＸは四面体の境
界であり、上記規格が少なくともこの線分セグメントに
対して満たされているので、これ以上のアクションは必
要ではない。）

【００２６】該線分セグメントが頂点Ｘを通過しないこ
とが判明すると（ボックス５４）、元のドローネー・メ
ッシュ生成時に計算されて格納された外接球Ｓが見つけ
られる（ボックス５８）。図１０に示されるように、外
接球Ｓは線分セグメントＡＢに付随する２個の点、すな
わち、頂点Ａと点Ｐを有する。外接球Ｓの等式だけでは
なく、頂点ＡとＢの座標も知られているので、線分セグ
メントＡＢがＳと交叉する点が見つけられる（図６のボ
ックス６０）。次に、線分セグメントＡＢの中間点Ｍが
見つけられ（ボックス６２）、点ＰとＭの位置は頂点Ａ
に関連させて比較される（ボックス６４）。点Ｐまたは
Ｍのいずれかが頂点Ａから離れていると判っても、その
点は新しい頂点Ｘに等しく設定され（ボックス６６、６
８）、更に、その値は、新しいセグメントＡＸを含む修
正された四面体を得るために、四面体に挿入される（ボ
ックス７０）。Ａの値は、Ｘに等しく割り当てられ（ボ
ックス７２）、プログラムは、モデルのすべてのエッジ
が考察されるまで、図６に示されるように継続される。

【００２７】上記の構造は、線分セグメントＡＢ上の点
ＰとＭの相対位置を示す図１０だけではなく、図１１に
おいても見られる。図１１においては、新しい点Ｘを示
している最終的な四面体が描かれている。上記の２個の
点が頂点Ａから離れている距離を選択することにより、
新しい四面体の集合が生成され、更に、モデル・エッジ
に沿った点の数が最小化されることが理解できる。

【００２８】この時点では、現行面のすべてのエッジが
四面体のエッジと一致することが知られている。但し、
四面体のいずれかのエッジがモデル面と交叉しているか
否かは判らない。これについては、判定される必要があ
り、そのような交叉が見つけられた場合には常に、「違
反している」四面体はトリミングされ、新しい三角形が
、エッジ／面の交叉を除去するモデルの面に付加されな
くてはならない。そのようなアクションを達成する手順
が、図１２と図１３のフロー図において示され、更に、
図１４から図２５に示される幾何学的構造によって描か
れている。

【００２９】図１２を考察すると、面トリミング手順は
、モデル・エッジ・セグメントＡＢを含むいずれかの四
面体を見つけることによって、開始され（ボックス１０
０）、次にそのような四面体がモデルの面Ｆと交叉する
エッジを有するか否かを調べる処理に進む。そのような
四面体は、モデルが面Ｆを有することを示している図１
４を調べることにより、明確になり、その四面体のエッ
ジＣＤは面Ｆと交叉している。この方法は、エッジＣＤ
が面Ｆを通過し、頂点Ｃが面Ｆ上にあり、更に、頂点Ｄ
が面Ｆの下にある四面体ＡＢＣＤを見つけるために、モ
デル・エッジ・セグメントＡＢを共有する四面体の間を
探索する。図１５では、モデル・エッジ・セグメントＡ
Ｂを共有するすべての四面体が斜視図で示され、図１６
では、図１５の１６−１６線に沿う断面図が表面の面Ｆ
と四面体エッジＣＤとの関係を示している。

【００３０】図１２に戻って考察するとき、頂点Ｄが面
Ｆに一致している場合、新しい頂点ＸはＤに等しく設定
され（ボックス１０６）、更に、新しい三角形ＡＢＸが
面Ｆに割り当てられる。図１５と図１６に示されるよう
に、頂点Ｄが面Ｆ上にないので、プログラムはエッジＣ
Ｄを共有する四面体Ｔ１−Ｔｎを探索する処理に進む（
ボックス１１０）。図１７では、四面体エッジＣＤが共
通四面体ＡＢＣＤ、ＢＥＣＤ、ＥＦＣＤ及びＦＡＣＤと
共に示されている。四面体Ｔ１−Ｔｎの各四面体毎にあ
るすべての外接球Ｓ１−Ｓｎが識別され（図１２のボッ
クス１１２）、面Ｆを有するＳ１−Ｓｎの交点の円Ｃ１
−Ｃｎも見つけられる（ボックス１１４）。次に、交点
の円Ｃ１−Ｃｎのすべてがオーバラップする領域Ｒが見
つけられる（ボックス１１６）。

【００３１】図１８から図２２では、ボックス１１１、
１１４、１１６で参照されている幾何学的構造が示され
ている。図１８から図２２で示される構造は、通常、四
面体メッシュに見られるような事象上に実質的に単純化
されることが理解されるだろう。但し、これらの構造は
、本発明の方法を描くのには充分である。

【００３２】図１８（ａ）では、四面体ＡＢＣＤを囲む
外接球Ｓ１が図示されている。更に、面Ｆ上の４個の四
面体Ｔ１−Ｔ４すべてのフットプリントＡＢＥＦが点線
で示されている。面Ｆの各頂点と外接球Ｓ１の等式が判
っているので、交点の円Ｃ１を見つけることができる。 交点の円Ｃ１と、及び、フットプリントＡＢＥＦとＣ１
と、の関係が図１８（ｂ）に示されている。

【００３３】図１９（ａ）では、四面体ＢＥＣＤのため
の外接球Ｓ２が、面Ｆを有する交点の円Ｃ２と共に示さ
れている。図１９（ｂ）では、交点の円Ｃ２が、フット
プリントＡＢＥＦと関連させて示されている。図２０（
ａ）と図２１（ａ）では、外接球Ｓ３とＳ４が、各々、
四面体ＥＦＣＤとＡＦＣＤに対して示されている。 更に、面Ｆを有する個々の交点の円Ｃ３とＣ４が、各々
、図２０（ａ）と図２１（ｂ）に示されている。図２２
は、交点の円Ｃ１−Ｃ４のすべての複合図であり、これ
らの交点の円がオーバラップする領域Ｒを示している。

【００３４】図１３を考察するとき、領域Ｒが一旦決定
されると、最初にこの領域がモデルの面Ｆのいずれかの
境界エッジと交叉しているか否かが判定される。交叉し
ている場合は、交点はＭとＮと指定される。次に、新し
い頂点Ｘの値は、これらの点ＭとＮの間の中間点座標に
等しく設定される。図２３を参照すると、そのような領
域Ｒ’が面Ｆのエッジ２００と交叉しているように示さ
れている。そこでは、エッジ２００に沿って存在する点
ＭとＮは交叉しているＲ’によって定義されることを見
ることができる。新しい頂点Ｘ’の値は、ＭとＮのあい
だの中間点に等しく設定される。一方、（図２３の）領
域Ｒが、セグメントＭ’Ｎ’を生成するために、面Ｆの
対称軸と交叉しているならば、ＸはＭ’Ｎ’の中間点に
等しく設定される。図２３では、部分的に対称軸２０２
が示されているが、図２４はより良く示されている。こ
の軸は、面Ｆの局所対称軸でなくてはならない。

【００３５】この場合、領域Ｒが対称軸２０２と交叉し
ていないので、プログラムはボックス１２０から図１３
のボックス１２２に処理を進め、新しい頂点Ｘの値を、
対称軸２０２に最も近い領域Ｒの境界上の点に等しく設
定する。点Ｘが一旦決定されると、この点は四面体に挿
入され（ボックス１２４）、三角形ＡＢＸが面Ｆの一部
分として復元される。

【００３６】この段階において、プログラムは、頂点Ｘ
が新しい点であるか否かを判定するために、図５のボッ
クス４０に復帰する。この後、プログラムは、再び、エ
ッジ・トリミング手順を実行するために再循環し、すべ
てのモデル・エッジ・セグメントが考察されるまで、処
理を継続する。上記の手順に従って、（図２５に示され
てあるように）個々の新しい三角形は、見つけられ、テ
ストされた後に、リストから除去されるので、残りの四
面体は後続の考察のために残しておく。最終的に、手順
は、すべての三角形が復元された後に終了する。

【００３７】新しい頂点の配置を決定するためにオーバ
ラップ領域Ｒを利用すると、メッシュのドローネー特性
を維持しつつ、頂点をメッシュ三角形の他の２個の頂点
からできる限り離して配置することが可能になる。従っ
て、本質的には、新しい頂点の数は最小化されるのであ
る。

【００３８】

【発明の効果】上記のように構成されるので、生成され
た四面体に対して、「全体が内部／外部にある」の規格
が自動的に満たされるように改善された有限要素システ
ムが提供されるという優れた効果を有する。

【図面の簡単な説明】

【図１】この方法を説明する上で用いられるモデルの斜
視図である。

【図２】基本有限要素、即ち、四面体を示す図である。

【図３】図１のモデルの面における初期のドローネー三
角形化の結果を示す図である。

【図４】（ａ）乃至（ｂ）は、四面体のための全体が内
部／外部の規格に対する「違反」を示す図である。

【図５】この発明を具体化するトリミング手順の高度な
フロー図である。

【図６】この発明のエッジ・トリミング手順を図示する
高度なフロー図である。

【図７】隣接するいずれかのメッシュ四面体のエッジと
一致しないモデル・エッジ・セグメントＡＢの斜視図で
ある。

【図８】８−８線に沿って切断された図７を示す図であ
る。

【図９】メッシュ四面体の各平面と一致する平面を利用
することによって、線分セグメントＡＢが交叉する四面
体面がどのようにして見つけられるのかを示す図である
。

【図１０】外接球が作成された後の図９において見つけ
られた四面体と、外接球と見つけられた線分セグメント
ＡＢとの間の交点を示す図である。

【図１１】追加頂点が一つ線分セグメントＡＢ上に配置
された後の図１０の四面体を示す図である。

【図１２】この発明の面トリミング手順を図示する高度
な流れ図である。

【図１３】この発明の面トリミング手順を図示する高度
な流れ図である。

【図１４】面Ｆと交叉する線分セグメントＣＤを有する
四面体ＡＢＣＤを示している、モデルの面Ｆの一部分の
斜視図である。

【図１５】四面体ＡＢＣＤの他に、該セグメントを共用
する他の四面体を示している線分セグメントＡＢの斜視
図である。

【図１６】１６−１６線に沿って切断された図１５を示
す図である。

【図１７】全部が内部／外部の規格に違反する４個の四
面体を有するモデル面の斜視図である。

【図１８】（ａ）は、　　四面体ＡＢＣＤの周囲に描か
れた外接球、及び、該モデル面と交叉するその円周を示
す図である。 （ｂ）は、外接球ＡＢＣＤがモデル面と交叉する円を示
す平面図である。

【図１９】（ａ）は、四面体ＢＥＣＤの周囲に描かれた
外接球を示す図である。 （ｂ）は、外接球ＢＥＣＤがモデル面と交叉する円を示
す図である。

【図２０】（ａ）は、四面体ＥＦＣＤの周囲に描かれた
外接球を示す図である。 （ｂ）は、外接球ＥＦＣＤがモデル面と交叉する円を示
す図である。

【図２１】（ａ）は、四面体ＡＦＣＤの周囲に描かれた
外接球を示す図である。 （ｂ）は、外接球ＡＦＣＤがモデル面と交叉する円を示
す図である。

【図２２】オーバラップする領域で、上記の交叉する円
のすべてを示す複合図である。

【図２３】追加された頂点の配置点を示し、その上に重
ねられたオーバラップ領域を有するモデル面を示す図で
ある。

【図２４】対象の中心線を示すモデルの斜視図である。

【図２５】後続の処理を実行する前に、追加される頂点
がモデル面上に配置され、作成された三角形が除去され
た後の該モデル面の平面図である。

【符号の説明】

Ａ　　　　　　　　点 ＡＢ　　　　　　エッジ・セグメント Ｆ　　　　　　　　面 ８４　　　　　　平面 Ｓ　　　　　　　　外接球 Ｘ　　　　　　　　頂点 Ｓ１　　　　　　外接球 Ｓ２　　　　　　外接球 Ｓ３　　　　　　外接球 Ｃ１　　　　　　交点の円 Ｃ２　　　　　　交点の円 Ｃ３　　　　　　交点の円 ２００　　　　エッジ ２０２　　　　対称軸

Claims (9)

【特許請求の範囲】
1. 【請求項１】　　エッジ及び面を表し、解析の対象であ
   る体の内部において、有限要素メッシュを完全に生成し
   、ソフトウェア・プログラムと組み合わされたコンピュ
   ータ・システムによって実行される有限要素メッシュ生
   成方法であって、ａ．複数の頂点とエッジを有する個々
   の有限要素の初期メッシュを前記体に対して生成するス
   テップと、ｂ．要素エッジと一致しないいずれかの体エ
   ッジ・セグメントを見つけ、更に、前記体エッジ・セグ
   メント上に数を最小限にするように位置づけて頂点を追
   加するステップと、ｃ．前記追加された頂点を隣接する
   頂点に接続することにより、追加の有限要素を生成する
   ステップと、ｄ．前記体のエッジ・セグメントと一致す
   るエッジを有する有限要素を見つけ、更に、体の一つの
   面が別のエッジを通過するように前記別のエッジを定義
   する頂点の間に配置されるステップと、ｅ．前記別のエ
   ッジと前記体面との交点付近で前記体面上に数を最小限
   にするように位置づけて新しい頂点を配置し、この新し
   く追加された頂点を近接の頂点に接続することにより、
   更に、有限要素を生成するステップと、から成る有限要
   素メッシュ生成方法。
2. 【請求項２】　　ｆ．個々の有限要素が完全に前記体の
   面内、又は完全に前記体の面外に存在するようになるま
   で、ステップｅにおいて新しい頂点が配置及び接続され
   る度に、ステップａからステップｅまでを繰り返し実行
   するステップを、更に含む請求項１に記載の有限要素メ
   ッシュ生成方法。
3. 【請求項３】　　個々の前記有限要素が四面体である請
   求項２に記載の有限要素メッシュ生成方法。
4. 【請求項４】　　前記ステップｂが、ｂ１．前記体エッ
   ジ・セグメントの第一の端と一致する頂点を有する四面
   体を見つけるステップであって、前記体エッジ・セグメ
   ントは、前記四面体の基底と前記基底を定義する頂点以
   外の位置で交叉するステップ、ｂ２．前記四面体によっ
   て定義される外接球を見つけるステップと、ｂ３．前記
   体エッジ・セグメントの前記外接球との交点Ｐを見つけ
   るステップと、ｂ４．前記体エッジ・セグメントの中間
   点Ｍを見つけるステップと、ｂ５．前記体エッジ・セグ
   メントの両端の間で前記体エッジ・セグメント上に新し
   い頂点を挿入するステップであって、前記新しい頂点は
   点Ｍに位置され、前記点Ｐが前記体エッジの第一の端か
   らかけ離れていないならば、前記新しい頂点が点Ｐに位
   置されるステップ、から成る請求項３に記載の有限要素
   メッシュ生成方法。
5. 【請求項５】　　ステップｄとステップｅが、ｄ１．体
   エッジと一致するエッジを有する第一の四面体を見つけ
   るステップと、ｄ２．前記体面によって交叉されている
   エッジを有する四面体を見つけるために、前記エッジを
   共有する四面体の間を探索するステップと、ｅ１．前記
   体面が前記四面体の頂点において前記エッジと交叉して
   いるか否かを判定するステップであって、交叉している
   場合には、前記頂点を、前記共有されたエッジの各頂点
   と接続し、生成された三角形を前記面の一部分として指
   定するステップと、を更に含む請求項３に記載の有限要
   素メッシュ生成方法。
6. 【請求項６】　　前記ステップｅが、ｅ２．前記共有さ
   れたエッジを有する四面体に対する外接球の前記体面と
   交叉する各円を見つけるステップと、ｅ３．前記交叉す
   る各円のオーバラップする領域を判定するステップと、
   ｅ４．前記オーバラップ領域上の前記新しい頂点を前記
   体面に配置するステップと、を更に含む請求項４に記載
   の有限要素メッシュ生成方法。
7. 【請求項７】　　前記オーバラップ領域が前記体面のエ
   ッジと交叉すると共に、前記エッジに沿って交叉するエ
   ッジ・セグメントＭＮを生成し、前記新しい頂点は前記
   エッジ・セグメントＭＮの中間点に配置される請求項６
   に記載の有限要素メッシュ生成方法。
8. 【請求項８】　　前記オーバラップ領域が前記体面のエ
   ッジと交叉していないが、前記体面の対称的なアクセス
   と交叉して、前記対称的アクセス上の前記オーバラップ
   領域内にセグメントＭＮを生成するならば、前記新しい
   頂点をセグメントＭＮの中間点に配置する請求項６に記
   載の有限要素メッシュ生成方法。
9. 【請求項９】　　前記オーバラップ領域が前記体面のエ
   ッジとも、更に、前記体面の対称軸おも交叉していない
   ならば、前記頂点が、前記対称軸に最も近接している前
   記領域上の点に配置される請求項６に記載の有限要素メ
   ッシュ生成方法。