特許 7085748 処理装置、ＣＡＤモデルの特徴部分検出方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-06-09

(45)【発行日】2022-06-17

(54)【発明の名称】処理装置、ＣＡＤモデルの特徴部分検出方法及びプログラム

(51)【国際特許分類】

   G06F 30/10 20200101AFI20220610BHJP

   G06F 30/15 20200101ALI20220610BHJP

   G06F 30/23 20200101ALI20220610BHJP

   G06F 113/28 20200101ALN20220610BHJP

【ＦＩ】

G06F30/10 100

G06F30/15

G06F30/23

G06F113:28

【請求項の数】 4

(21)【出願番号】P 2018071823

(22)【出願日】2018-04-03

(65)【公開番号】P2019185186

(43)【公開日】2019-10-24

【審査請求日】2021-02-01

(73)【特許権者】

【識別番号】503361400

【氏名又は名称】国立研究開発法人宇宙航空研究開発機構

(74)【代理人】

【識別番号】110003339

【氏名又は名称】特許業務法人南青山国際特許事務所

(74)【代理人】

【識別番号】100104215

【弁理士】

【氏名又は名称】大森 純一

(74)【代理人】

【識別番号】100196575

【弁理士】

【氏名又は名称】高橋 満

(74)【代理人】

【識別番号】100168181

【弁理士】

【氏名又は名称】中村 哲平

(74)【代理人】

【識別番号】100117330

【弁理士】

【氏名又は名称】折居 章

(74)【代理人】

【識別番号】100160989

【弁理士】

【氏名又は名称】関根 正好

(74)【代理人】

【識別番号】100168745

【弁理士】

【氏名又は名称】金子 彩子

(74)【代理人】

【識別番号】100176131

【弁理士】

【氏名又は名称】金山 慎太郎

(74)【代理人】

【識別番号】100197398

【弁理士】

【氏名又は名称】千葉 絢子

(74)【代理人】

【識別番号】100197619

【弁理士】

【氏名又は名称】白鹿 智久

(72)【発明者】

【氏名】金森 正史

(72)【発明者】

【氏名】石田 崇

(72)【発明者】

【氏名】橋本 敦

(72)【発明者】

【氏名】青山 剛史

【審査官】堀井 啓明

(56)【参考文献】

【文献】特開２００４－２６５０５０（ＪＰ，Ａ）

【文献】特開平１０－０１１６１２（ＪＰ，Ａ）

【文献】特開平１１－００３４３９（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ ３０／００－３０／３９８

(57)【特許請求の範囲】

【請求項1】

ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータから、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成し、
前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求め、
前記求められた電位分布に基づき前記物体形状の特徴部分を検出する
特徴部分検出部
を具備する処理装置。

【請求項2】

請求項１に記載の処理装置であって、
前記ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータは、ＳＴＬ(Standard Triangulated Language/Standard Tessellation Language)のデータである
処理装置。

【請求項3】

コンピュータが、ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を再現し、
コンピュータが、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成し、
コンピュータが、前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求め、
コンピュータが、前記求められた電位分布に基づき前記物体形状の特徴部分を検出する
ＣＡＤモデルの特徴部分検出方法。

【請求項4】

ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータから、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成するステップと、
前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求めるステップと、
前記求められた電位分布に基づき前記物体形状の特徴部分を検出するステップと
をコンピュータに実行させるプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、例えば航空機の機体の数値解析に必要な計算格子を得るための処理装置、ＣＡＤモデルの特徴部分検出方法及びプログラムに関する。

【背景技術】

【0002】

計算力学分野においては、計算領域を有限の離散点に分割する作業、いわゆる計算格子生成が重要な作業である。離散点が適切な位置に配置されていない場合、計算結果に非物理的な誤差が生じることがある。計算格子とは、すなわち物体周りにくまなく張り巡らされた網目であり、物体の形をきれいに覆う必要があり、きれいに覆うことができないと、不自然な「バリ」が出て計算精度が極端に悪化する。

【0003】

特に、物体形状の稜線や物体同士の交差部のような部分(これを「特徴部分」と呼ぶ)では、その特徴部分に離散点を配置しないと、形状の特徴を表した数値計算は実施できない。

【0004】

従来、特徴部分は人間の目により発見し、その部分を正確に再現できるように格子網を生成する作業が行われてきた。しかし、モデルが複雑になるにつれて人間の作業には時間的な限界を生じ、特徴部分の自動検出が望まれてきた。

【0005】

従来技術として、モデル表面の角度や曲率を図る手法がある。ＣＡＤモデルは小さな三角形や四角形などの面を無数に並べることで物体の形状を再現しており、隣り合う面同士の角度や曲率を計算し、それらがある閾値よりも大きい場合、その部分を特徴部分として検出する方法である。

【0006】

より具体的には、例えば隣り合う三角形のなす角を、各三角形の単位法線ベクトルから計算する。つまり、

として、θがある閾値θ_threshよりも大きい場合、特徴部分と判断する（非特許文献１参照）。

【先行技術文献】

【非特許文献】

【0007】

【文献】'Hexahedra Grid Generation Method for Flow Computation' a Paulus R. Lahur of Japan Aerospace Exploration Agency, 22nd Applied Aerodynamics Conference and Exhibit, 16-19 August 2004, Providence, Rhode Island

【発明の概要】

【発明が解決しようとする課題】

【0008】

上記の手法は簡便だが、単純な例でも検出されないようなパターンが存在する。例えば、θ_thresh＝60degと設定したときに、２平面が30deg程度と緩やかに交差するような特徴部分を検出することはできない。つまり、隣り合った三角形の関係だけでは、特徴部分を適切に検出することはできない。

【0009】

以上のような事情に鑑み、本発明の目的は、物体形状の特徴部分を適切に顕在化することができる処理装置、ＣＡＤモデルの特徴部分検出方法及びプログラムを提供することにある。

【課題を解決するための手段】

【0010】

上記目的を達成するため、本発明の一形態に係る処理装置は、ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータから、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成し、前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求め、前記求められた電位分布に基づき前記物体形状の特徴部分を検出する特徴部分検出部を具備する。

【0011】

本発明では、離散点によって囲まれる各面を各面の法線方向に移動させることで、膨張モデル又は収縮モデルを生成し、それをラプラス方程式で解くことで、物体形状の特徴部分を強調することが可能となるので、特徴部分を適切に顕在化することができる。

【0012】

ここで、前記ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータは、ＳＴＬ(Standard Triangulated Language/Standard Tessellation Language)のデータであってもよい。

【0013】

本発明の一形態に係るＣＡＤモデルの特徴部分検出方法は、ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を再現し、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成し、前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求め、前記求められた電位分布に基づき前記物体形状の特徴部分を検出する。

【0014】

本発明の一形態に係るプログラムは、ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる面を並べることで前記ＣＡＤモデルの物体形状を表現したデータから、各面を各面の法線方向に移動させることで、前記ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成するステップと、前記膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで前記膨張モデル又は収縮モデルの表面の電位分布を求めるステップと、前記求められた電位分布に基づき前記物体形状の特徴部分を検出するステップをコンピュータに実行させる。

【発明の効果】

【0015】

本発明により、物体形状の特徴部分を適切に顕在化することができる。これにより、従来手作業で実施することが必要であった特徴部分の取り扱いが自動化され、その結果、ＣＡＤモデルを用意するだけで格子網を生成し、数値解析を実施するところまでを自動的に進めることできる。このことで、数値解析全体にかかるターンアラウンドは減少し、最適化等の手間のかかる作業を大幅に削減することが可能になる。

【図面の簡単な説明】

【0016】

【図1】本発明の一実施形態に係る処理装置の構成を示すブロック図である。

【図2】ＳＴＬのデータによって表現されたＣＡＤモデルとしての航空機の機体を示す斜視図である。

【図3】図２の符号Ａの部分の拡大図である。

【図4】図２の符号Ｂの部分の拡大図である。

【図5】ＳＴＬのデータによって表現される三角形を示す図である。

【図6】ＳＴＬのデータの一例を示す図である。

【図7】図１に示した処理装置における処理の流れを示すフローチャートである。

【図8】図２～図４に示した三角形のうち一つの三角形に着目し、その三角形を法線ベクトル方向に移動する状態を示す図である。

【図9】機体の凸部に着目した膨張モデルを示す図である。

【図10】機体の凹部に着目した膨張モデルを示す図である。

【図11】特徴部分を検出したＣＡＤモデルの外観を示す斜視図である。

【図12】膨張させずに電位を求めた結果を示す図である。

【図13】膨張させて電位を求めた結果を示す図である。

【図14】図２～図４に示した三角形のうち一つの三角形に着目し、収縮に向かうようにその三角形を法線ベクトル方向に移動する状態を示す図である。

【発明を実施するための形態】

【0017】

以下、図面を参照しながら、本発明の実施形態を説明する。
図１は、本発明の一実施形態に係る処理装置を示す図である。
図１に示すように、処理装置１０は、データ変換部１１と、特徴部分検出部１２と、数値解析部１３とを有する。この処理装置１０は、典型的には、コンピュータシステムに各部を構成するプログラムをインストールすることによって構成される。

【0018】

データ変換部１１は、ＣＡＤデータからＳＴＬのデータに変換する。
ここで、ＳＴＬ(Standard Triangulated Language/Standard Tessellation Language)とは、３次元の任意の表面形状を、無数の三角形で表現するＣＡＤデータの形式の一つである。本実施形態に係るＳＴＬのデータは、ＣＡＤモデルの表面を格子網に離散化し、離散点によって囲まれる三角形を並べることでＣＡＤモデルの物体形状を表現したデータである。図２に、ＳＴＬのデータによって表現されたＣＡＤモデルとしての航空機の機体の形状を示し、図３及び図４にその一部拡大部分（図２のＡとＢ）を示す。
これらのデータは、図５に示すように、三角形ごとに、その三角形の単位法線ベクトル及びその三角形の３頂点の座標を記載したデータとなる。データフォーマットの一例を図６に示す。

【0019】

特徴部分検出部１２は、各三角形を各三角形の法線方向に移動させることで、ＣＡＤモデルの物体形状の膨張モデル又は収縮モデルを生成し、膨張モデル又は収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで膨張モデル又は収縮モデルの表面の電位分布を求め、求められた電位分布に基づき物体形状の特徴部分を検出する。

【0020】

数値解析部１３は、例えば機体表面の抵抗を算出する。より具体的には、上記の特徴部分検出部１２により検出された特徴部分を使って、物体の特徴を保ったまま計算格子を生成し、その計算格子を用いて数値流体力学解析を実施することで、機体表面の圧力分布が得られ、それを機体表面全体で積分することで抵抗を算出する。

【0021】

次に、処理装置１０における処理の具体例を図７に示すフローチャートに沿って説明する。
データ変換部１１において、ＣＡＤデータからＳＴＬのデータに変換し、ＳＴＬのデータを用意する（ステップ７１）。なお、本発明に係る処理装置では、ＳＴＬのデータが予め用意される場合には、データ変換部１１は不要である。

【0022】

次に、特徴部分検出部１２において、ＣＡＤモデルを拡大した膨張モデルを作成する（ステップ７２）。
具体的には、物体表面を構成する各三角形の頂点ｘ（ベクトル）_ｉに対して、新しい頂点の位置をｘ（ベクトル）_ｉ＋ｓｎ（ベクトル）とする。
ここで、ｎ（ベクトル）は頂点に対応する三角形の法線ベクトルであり、ｓは膨張量を表す。膨張量は、各ＣＡＤモデルによって切り替える必要があるが、概ね全三角形の辺の長さの最小値に対して５倍程度が望ましい。なお、後述の縮小量もこれと同程度と考えられる。

【0023】

図８は図２～図４に示した三角形のうち一つの三角形に着目し、その三角形を法線ベクトル方向に移動する状態を示す図である。
ここでは、ｓ＞０で、ｎ（ベクトル）が機体の表面から外方向に向かっているとすると、ｘ（ベクトル）_ｉ＋ｓｎ（ベクトル）によって、三角形は機体の表面から外方向に移動することになる。

【0024】

ここで、図９に機体の凸部に着目した膨張モデルを示し、図１０に機体の凹部に着目した膨張モデルを示す。機体の凸部とは例えば図３に示した翼の端部などであり、機体の凹部とは例えば図４に示した機体の胴体と翼の境目などである。
三角形が機体の表面から外方向に移動すると（図中矢印の方向）、機体の凸部においては、図９に示すように、機体Ａ表面で隙間なく並んでいた三角形Ｔはこの移動によって三角形Ｔ同士が離れるため隙間Ｇをもった膨張モデルＥとなる。機体の凹部においては、図１０に示すように、機体Ａ表面で隙間なく並んでいた三角形Ｔはこの移動によって三角形Ｔ同士が近づきあうため交差Ｉをもった膨張モデルＥとなる。

【0025】

以上の膨張モデルＥを細かく見ると、凹凸が大きいほど交差Ｉ及び隙間Ｇが大きくなることがわかる。

【0026】

次に、膨張モデルＥの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで膨張モデルＥの表面の電位分布を求める（ステップ７３）。

【0027】

ここで、以下の説明において、記号は以下の意味である。
φ：静電場

ｘ（ベクトル）：表面を指す位置ベクトル
ｎ：表面の法線方向
Ｄ：解析領域

ｒ：表面上の点ｉとｊの距離

【0028】

静電場と電位の関係は、ラプラス方程式によって支配されることが知られている。

【0029】

表面上に一様に静電場φ（ベクトル）が分布しているものとする。

【0030】

(2)式を境界条件として用い、(1)式はラプラス方程式のグリーン関数を用いて次のように変形できる。

【0031】

(3)式は境界要素法の基本式であり、これを表面三角形を要素として離散化し、ｑを求める。

【0032】

全ての三角形においてφは一定値であり、また各三角形においてｑが一定値をとると仮定すると（これは境界要素を１次要素とすることに相当)、（４）式は次のように変形できる。

【0033】

(5)式を行列形式に置き換えると次式が得られる。

【0034】

ここで、

【0035】

(6)式は未知数ベクトルｑに関する連立1次方程式であり、これを解けば各三角形上のｑが得られる。

【0036】

なお、以上に示した境界要素法の一例にすぎず、本発明はこれに限定されず、これ以外の境界要素法を用いても勿論かまわない。

【0037】

次に、境界要素法により求めた表面電位ｑ（電位分布）に対して、その絶対値｜ｑ｜が大きい部分を抽出する（ステップ７４）。つまり、これにより電位分布に基づき物体形状の特徴部分を検出する。

【0038】

絶対値｜ｑ｜を考えるのは、特徴部分が凹の部分で絶対値の大きな負の値を、凸の部分で絶対値の大きな正の値をとるためであり、凹凸の両者を同時に顕在化したいため、ｑ単体ではなく｜ｑ｜を用いる。

【0039】

ここで、図１１は特徴部分を検出したＣＡＤモデルの外観を示す斜視図であり、図１２に膨張させずに電位を求めた結果を示し、図１３に膨張させて電位を求めた結果を示す。
図１２の結果と図１３の結果を比較すると、膨張させずに電位を求めた結果は特徴部分が殆ど顕在化されていない（図１２）のに対して、膨張させて電位を求めた結果は凹凸部を顕在化できている（図１３）ことがわかる。

【0040】

本発明は上記の実施形態に限定されず、その発明の技術的思想の範囲内で様々な変形や応用が可能である。その変形や応用による実施も本発明の技術的範囲に属する。

【0041】

例えば、上記の実施形態では、ＣＡＤモデルを拡大した膨張モデルを作成していたが、ＣＡＤモデルを縮小した収縮モデルを作成し、収縮モデルの表面に一様な電荷が帯電しているものとみなし、境界要素法による表面上の静電場に対するラプラス方程式を解くことで収縮モデルの表面の電位分布を求め、特徴部分の検出を行ってもよい。

【0042】

図１４に図２～図４に示した三角形のうち一つの三角形に着目し、収縮に向かうようにその三角形を法線ベクトル方向に移動する状態を示す。例えば、ｓ＜０で、ｎ（ベクトル）が機体の表面から外方向に向かっているとすると、ｘ（ベクトル）_ｉ＋ｓｎ（ベクトル）によって、三角形は機体の表面から内方向に移動することになる。
そして、収縮させて電位を求めることで、膨張の場合と同様に、凹凸部である特徴部分を顕在化できる。

【0043】

さらに、上記の実施形態では、三角形の単位法線ベクトル及びその三角形の３頂点の座標の情報は、すべてＳＴＬのデータをそのまま利用していたが、三角形の単位法線ベクトルについては例えば次の手順で計算することも可能である。

【0044】

ここでは「×」ベクトル同士の外積を表し、ｒ（ベクトル）とｓ（ベクトル）の外積がｒ（ベクトル）とｓ（ベクトル）の両方に垂直な方向になることを利用した計算法である。
従って、本発明は、ＳＴＬのデータに三角形の単位法線ベクトルが含まれていない場合であっても実施が可能である。

【0045】

さらにまた、上記の本実施形態では、三角形を並べることでＣＡＤモデルの物体形状を表現したＳＴＬのデータを用いていたが、本発明は、四角形やそれ以上の多角形を並べることでＣＡＤモデルの物体形状を表現したデータであっても実施可能である。

【符号の説明】

【0046】

１０ 処理装置
１１ データ変換部
１２ 特徴部分検出部
１３ 数値解析部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】