JP4005803B2 - 流体解析方法、及び、その流体解析方法を用いた流体解析装置 - Google Patents
流体解析方法、及び、その流体解析方法を用いた流体解析装置 Download PDFInfo
- Publication number
- JP4005803B2 JP4005803B2 JP2001377582A JP2001377582A JP4005803B2 JP 4005803 B2 JP4005803 B2 JP 4005803B2 JP 2001377582 A JP2001377582 A JP 2001377582A JP 2001377582 A JP2001377582 A JP 2001377582A JP 4005803 B2 JP4005803 B2 JP 4005803B2
- Authority
- JP
- Japan
- Prior art keywords
- subdivision
- calculation
- fluid analysis
- area
- analysis method
- 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 - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Description
【０００１】
【発明の属する技術分野】
本発明は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、上記設計領域を複数の演算領域に細分割して流体の解析を行う流体解析方法、及び、その流体解析方法を用いた流体解析装置に関する。
【０００２】
【従来の技術】
従来より、空調ダクト等の設計においては、限られたスペース内で配風要件を満足させ、しかも圧力損失の少ないダクト形状を検討することが必要である。このようなことから、本発明者は、社団法人自動車技術会 学術講演会前刷集９６４ Ｎｏ．９６３６８９９（１９９６年１０月）において、圧力損失の少ないダクトを最適ダクトと定義し、限られたダクト容積で出入り口間を連結する概念形状を求めるダクトレイアウト最適化の手法を提案している。
【０００３】
このダクトレイアウト最適化の手法では、設計領域として、図１２に示すように、内部に小さな固体が無数に散らばっている多孔質体を考え、この多孔質体の空隙率に相当する流動抵抗が流体に作用するものとする。そして、例えば図１３に示すような検討空間（設計領域）を、図１４（ａ）、図１４（ｂ）、或いは、図１４（ｃ）に示すような、複数の演算領域に細分割（複数の計算セルに細分割）して考え、設計領域内に占める物体の容積を一定に保ちながら（図１５におけるダクト容積：これが制約条件となる）、多孔質体の空隙率をコントロールすることで圧力損失の少ない概念形状を求める手法である。すなわち、圧力勾配の大きい、流れやすい計算セルの空隙率を大きく、逆に圧力勾配の小さい、あまり流れのない計算セルの空隙率を小さくコントロールする。こうして得られる結果から、例えば、空隙率が０．５の等値面を境界部位としてポスト処理することで物体形状（概念形状）を求めるようになっている。
【０００４】
【発明が解決しようとする課題】
しかしながら、上述のダクトレイアウト最適化の手法では、計算当初に設定しておいた計算セルによっては、図１４（ａ）に示すように計算セルのサイズが粗い場合には、図１６（ａ）に示すように物体形状が滑らかに表現することができないという問題があり、逆に図１４（ｂ）に示すように各計算セルのサイズが細かすぎると演算負荷が大きくなるという問題がある。また、図１４（ｃ）に示すように各計算セルの密度が不均一な場合には、形状表現にムラを生じてしまい、図１６（ｃ）に示すように、精度の良い滑らかな物体形状が得られないという問題がある。精度の良い滑らかな物体形状を得るために、予め必要な部位の計算セルのサイズを細かくしておくことも考えられるが、予め必要な部位を予想して他の部位と調整しながら細かく分割しておくことは困難で、このような作業を行うことは解析方法の汎用性を低下させる一因となってしまう。
【０００５】
本発明は上記事情に鑑みなされたもので、汎用性に優れ、また、演算負荷が必要以上に増加することなく、且つ、精度良く滑らかな物体形状を解析し得ることが可能な流体解析方法、及び、その流体解析方法を用いた流体解析装置を提供することを目的とする。
【０００６】
【課題を解決するための手段】
上記目的を達成するため、本発明の請求項１記載の流体解析方法は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、上記設計領域を複数の演算領域に細分割して流体の解析を行う流体解析方法であって、上記解析対象とする設計領域には、更なる細分割を実行する領域と更なる細分割を禁止する領域の少なくともどちらかを予め設定自在で、上記複数の演算領域の細分割を変更自在とすることを特徴とする。
【０００７】
すなわち、請求項１記載の流体解析方法は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、設計領域を複数の演算領域に細分割し、更に、必要に応じて複数の演算領域の細分割を変更して、流体の解析を行う。このため、様々な設計領域で応用が可能で汎用性に優れ、また、演算領域の細分割の変更は、精度良く滑らかな物体形状を解析する部位にのみ行えば良いため、演算負荷が必要以上に増加することがない。また、解析対象とする設計領域には、更なる細分割を実行する領域と更なる細分割を禁止する領域の少なくともどちらかを予め設定自在であり、ユーザがより深く解析したい部位や、或いは、解析を必要としない部位における演算の節約を図れるように満足して、ユーザにとっての利便性をより向上させる。
【０００８】
また、本発明の請求項２記載の流体解析方法は、請求項１記載の流体解析方法において、上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順を有し、上記所定のパラメータが予め設定しておいた範囲の値の演算領域を更に細分割することを特徴とする。
【０００９】
すなわち、請求項２記載の流体解析方法は、請求項１記載の流体解析方法において、複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算し、所定のパラメータが予め設定しておいた範囲の値の演算領域を更に細分割する。
【００１０】
更に、本発明の請求項３記載の流体解析方法は、請求項１又は請求項２に記載の流体解析方法において、上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順と、上記所定のパラメータが予め設定しておいた値の部位を境界部位として連続させる手順とを有し、上記境界部位の形状の変化量が予め設定する値を超える際に、該境界部位が存在する演算領域を更に細分割することを特徴とする。
【００１１】
すなわち、請求項３記載の流体解析方法は、請求項１又は請求項２に記載の流体解析方法において、複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算し、所定のパラメータが予め設定しておいた値の部位を境界部位として連続させる。そして、境界部位の形状の変化量が予め設定する値を超える際に、境界部位が存在する演算領域を更に細分割する。
【００１４】
更に、本発明の請求項４記載の流体解析方法は、請求項１乃至請求項３の何れか一つに記載の流体解析方法において、上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、元の領域の細分割を実行して形成した領域を元の領域に対する１世代後の子の領域とし、該子の領域の細分割を実行して形成した領域を上記元の領域に対する２世代後の孫の領域として複数世代で形成する場合、２世代以上離れて形成した領域間の隣接を禁止することを特徴とする。
【００１５】
すなわち、請求項４記載の流体解析方法は、請求項１乃至請求項３の何れか一つに記載の流体解析方法において、細分割を繰り返すことによる各演算領域間の接続が不自然となることを防止する。
【００１６】
また、本発明の請求項５記載の流体解析方法は、請求項１乃至請求項４の何れか一つに記載の流体解析方法において、上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、上記細分割した後に上記細分割の対象外となった演算領域は、上記細分割の前の大きさの演算領域に復元自在であることを特徴とする。
【００１７】
すなわち、請求項５記載の流体解析方法は、請求項１乃至請求項４の何れか一つに記載の流体解析方法において、細分割を繰り返すことにより、不要なデータが多く蓄積されることを未然に防止できるようにする。
【００１８】
更に、本発明の請求項６記載の流体解析装置は、上記請求項１乃至請求項５の何れか一つに記載の流体解析方法を備え、流体の解析を行うことを特徴とする。
【００１９】
【発明の実施の形態】
以下、図面に基づいて本発明の実施の形態を説明する。図１〜図１１は本発明の実施の一形態に係わり、図１は流体解析装置の機能ブロック図、図２は検討空間の初期計算セルの説明図、図３は細分割処理に係る２つの領域についての説明図、図４は図３の２つの領域の組み合わせについての説明図、図５は細分割する計算セルの形状の説明図、図６は計算セルの細分割処理の世代の説明図、図７は細分割処理を実行する条件の他の例の説明図、図８は細分割処理により新たに細分割処理が必要となる場合の説明図、図９は流体解析プログラムのフローチャート、図１０は２次元四角形セルによる流体解析結果の説明図、図１１は２次元三角形セルによる流体解析結果の説明図である。
【００２０】
図１において、符号１はマイクロコンピュータとその周辺回路を備えて構成された流体解析装置を示し、この流体解析装置１は、初期空隙率設定部２、流れ解析実行部３、圧力空間微分値演算部４、空隙率更新部５、計算セル細分割判断実行部６、更新判断部７、境界部位形成部８、データ出力部９、データ保存部１０から主要に構成されている。尚、本発明の実施の形態においては、検討空間（設計領域）として、例えば図１３に示す空間、すなわち、左側の壁に垂直吹き出しの入口を有し、右側の壁に自由流出の出口を備えた空間を例に説明する。
【００２１】
初期空隙率設定部２は、予め設定された後述の制約ボリュームが入力され、この制約ボリュームが初期空隙率として設定されて、流れ解析実行部３に出力される。制約ボリュームは、図１５に示すように、検討空間内において予め指定された、占有することのできるダクト容積を検討空間の全容積で除した（＝（ダクト容積）／（検討空間容積））ものである。
【００２２】
流れ解析実行部３は、検討空間の境界条件や、初期計算セルの形状等が入力され、以下の多孔質体内における流れの数値計算を実行する。そして、この数値計算により得られた各計算セルの圧力値ｐは、圧力空間微分値演算部４へと入力される。
・連続の式
∂（ε・ｕi）／∂ｘi＝０ …（１）
・流れの運動方程式
ｕj・（∂（ε・ｕi）／∂ｘj）＝−∂（ε・ｐ）／∂ｘi＋（１／Ｒｅ）・（∂２（ε・ｕi）／∂ｘj２）−Ｒi …（２）
ここで、計算セルの空隙率εを次のように定義する。
【００２３】
ε＝（流体領域の体積）／（計算セル体積） …（３）
ｕ：速度成分、ｘ：座標軸、ｐ：圧力、Ｒｅ：レイノルズ数、Ｒ：流動抵抗（＝ｆ（ε，ｕｉ））、添え字ｉ，ｊは、総和規則に従う（＝１，２，３）
【００２４】
また、初期計算セルの形状は、図２に示すように、２次元における解析では、複数の四角形の計算セル（図２（ａ））、或いは、複数の三角形の計算セル（図２（ｂ））となっている。
【００２５】
更に、図２の実線の閉空間で示すように、検討空間には、例えば長方形（図２（ａ））、楕円（図２（ｂ））、円及びこれらの組み合わせ形状による幾何学的領域が、ユーザにより任意に設定自在になっており、計算セルの中心が、この領域内にある計算セルを、強制的に更に細分割することや、或いは、細分割することを禁止することが可能になっている。これは、例えば、細分割する領域が図３（ａ）の斜線のように得られ、強制的に更に細分割する領域を図３（ｂ）の斜線のように設定しておいたものとすると、結果として細分割される領域は図４（ａ）に示す斜線領域となる。一方、細分割する領域が図３（ａ）の斜線のように得られ、細分割することを禁止する領域を図３（ｂ）の斜線のように設定しておいたものとすると、結果として細分割される領域は図４（ｂ）に示す斜線領域となる。このような、細分割を選択可能な領域を設けることにより、ユーザが、より高度で緻密な解析を行いたい部位での利便性を図ることができ、また、解析にあまり必要としないデータ部分での演算の削減を行うことが可能となっている。
【００２６】
尚、図示はしないが、３次元における解析では、初期計算セルの形状は、複数の６面体の計算セル、或いは、複数の四面体の計算セル等である。この場合においても、検討空間には、例えば円柱や、球及びこれらの組み合わせ形状による幾何学的領域が、ユーザにより任意に設定自在になっており、計算セルの中心が、この領域内にある計算セルを、強制的に更に細分割することや、或いは、細分割することを禁止することが可能になっている。
【００２７】
圧力空間微分値演算部４は、流れ解析実行部３から各計算セルの圧力値ｐが入力され、これら圧力値ｐの空間微分を演算し、空隙率更新部５に出力する。
すなわち、
∂ｐi／∂ｎ＝（Σ（∂ｐi／∂ｘi）２）（１／２） …（４）
【００２８】
空隙率更新部５は、圧力空間微分値演算部４から各計算セルの空間微分値（∂ｐi／∂ｎ）が入力され、空間微分値（∂ｐi／∂ｎ）の大きい方から計算セル容積の和を算出し、計算セル容積の和が制約ボリュームと等しくなった時点の空間微分値（∂ｐi／∂ｎ）（＝（∂ｐ／∂ｎ）ε0.5とする）を、空隙率＝０．５に更新する。
【００２９】
そして、空間微分値（∂ｐi／∂ｎ）が、（∂ｐ／∂ｎ）ε0.5より大きい計算セルの空隙率εを０．５以上に、小さい計算セルの空隙率εを０．５以下に更新し、これを繰り返し計算することで空隙率εの分布を更新していく。最終的に得られた空隙率分布から、ε＝０．５を物体表面形状（ε＞０．５を流体領域、ε＜０．５を物体領域）と定義するのである。尚、本実施の形態の例では、空隙率εが０．５で設定しているが、制約ボリューム等の条件によっては、他の値、例えば、０．３や、０．２として演算しても良い。
【００３０】
このように、空隙率更新部５で更新される各計算セルの空隙率の分布状態は、計算セル細分割判断実行部６へと出力される。
【００３１】
計算セル細分割判断実行部６は、空隙率更新部５から各計算セルの空隙率の分布状態が入力され、また、細分割する条件（初期条件：例えば、細分割を行う計算セルの空隙率の範囲εmin〜εmax）が入力されて、各計算セルの細分割の実行可否を判断し、細分割する計算セルを抽出して細分割処理し、細分割処理を終えたデータを更新判断部７に出力する。
【００３２】
ここで、計算セルの細分割は、例えば図５に示すように、計算セルの形状によって異なっている。例えば、計算セルが２次元の四角形であれば、図５（ａ）に示すように、細分割する計算セルは、四角形を４等分するものであり、これを更に細分割する計算セルも、細分割した計算セルを４等分するものとなっている。
【００３３】
また、計算セルが２次元の三角形であれば、図５（ｂ）に示すように、細分割する計算セルは、三角形を４等分するものであり、これを更に細分割する計算セルも、細分割した計算セルを４等分するものとなっている。
【００３４】
更に、計算セルが３次元の６面体であれば、図５（ｃ）に示すように、細分割する計算セルは、６面体を８等分するものであり、これを更に細分割する計算セルも、図示はしないが、細分割した計算セルを８等分するものとなっている。
【００３５】
また、計算セルが３次元の４面体であれば、図５（ｄ）に示すように、細分割する計算セルは、４面体を８等分するものであり、これを更に細分割する計算セルも、図示はしないが、細分割した計算セルを８等分するものとなっている。
【００３６】
細分割する計算セルは、図６に示すように、本実施の形態においては、世代で定義する。すなわち、初期の計算セルは、複数の親セルから構成され、この親セルから細分割されたセルは、親セルに対して１代下の代として子セルとして定義する。同様に、子セルから細分割されたセルは、親セルに対して２代下の代として孫セルとして定義する。更に、孫セルから細分割されたセルは、親セルに対して３代下の代としてひ孫セルとして定義する。
【００３７】
計算セル細分割判断実行部６では、空隙率が、予め細分割する条件、すなわち、予め設定しておいた初期条件の範囲にある計算セルを抽出し、細分割処理を行うが、以下のような条件としても良い。
【００３８】
空隙率が、予め設定した空隙率（例えば０．５）でポスト処理される境界部位（以降、等値面という）の形状変化率が、予め設定する値を超える際に、この等値面が存在する計算セルのみを細分割する。
【００３９】
すなわち、図７（ａ）に示すように、まず、ポスト処理にて空隙率が０．５の等値面を形成する。次いで、これら等値面の各計算セル毎の法線ベクトルｎ１〜ｎ１１を算出し、これら法線ベクトルｎ１〜ｎ１１の隣接するベクトル同士のなす角が、予め設定しておいた角度を超えている場合、その法線ベクトルの等値面が所属する計算セルを、図７（ｂ）に示すように細分割処理する。そして、これを繰り返していくことで、形状変化の大きな部分が滑らかに表現できるようになっている。尚、３次元においても、図７（ｃ）に示すように、上述の２次元の細分割処理と同様に、ポスト処理にて形成した空隙率０．５の各等値面に法線ベクトルｍ１〜ｍ４を算出し、これら各法線ベクトルの隣接する法線ベクトル同士のなす角を求め、設定値と比較することにより、等値面が所属する計算セルを細分割処理することで行うことができる。
【００４０】
また、計算セル細分割判断実行部６では、細分割する計算セルは、滑らかな連続性を維持するために、２世代以上離れた計算セルと隣接しないようにする。特定の計算セルの細分割により、２世代以上離れるような場合は、図８に示すように、例えば、子セルの細分割により新たな孫セルが親セルと隣接するような場合は、この孫セルと隣接する親セルも子セルに細分割する。
【００４１】
また、計算セル細分割判断実行部６では、計算セルの細分割を繰り返し、過去に細分割した計算セルにおいて、現在の物体形状表現に使われていない（計算に関与していない）計算セルは、適宜、１世代前の計算セルに戻すことが可能になっている。このようにすることで、全体の計算セル数の削減を図り、データ数を削減することが可能になっている。
【００４２】
更新判断部７は、当初入力された演算の更新回数（例えば、１００回等）に達しているか否か判定し、演算回数が更新回数に達している場合には、演算により得られたデータを境界部位形成部８に出力する。
【００４３】
境界部位形成部８では、入力された各計算セルの空隙率のデータに対し、予め設定しておいた物体表面形状とみなす空隙率、すなわち、０．５の等値面を表示させる。
【００４４】
データ出力部９では、境界部位形成部８からのデータを、図示しないモニタに出力させる。また、必要に応じて、等値面を構成する面群データ或いは点群データからＣＡＤ（Computer Aided Design）のSurfaceデータ又はVolumeデータを作成しておき、ハードディスク等のデータ保存部１０に、ＩＧＥＳ形式、ＳＴＥＰ形式、或いは、ＳＴＬ形式等のＣＡＤで多く用いられるデータ形式で保存する。
【００４５】
そして、ダクト等の設計の際に、ハードディスク等のデータ保存部１０から、直接、保存しておいた最適な物体表面形状のデータを読み込むことにより、流体解析の設計へのフィードバックが迅速に行われるようになっている。
【００４６】
次に、上記構成の流体解析装置１で実行される流体解析について、図９のフローチャートで説明する。
まず、ステップ（以下、「Ｓ」と略称）１０１では、解析に必要なデータの入力が行われる。これは、検討空間の形状や入口・出口の境界条件、初期計算セルの形状、ダクトの制約ボリューム、演算の更新回数、細分割する条件（細分割する計算セルの空隙率の範囲）等である。
【００４７】
次いで、Ｓ１０２では、初期空隙率設定部２において、初期空隙率の設定を行う。本実施の形態において、この初期空隙率は、制約ボリュームが初期空隙率として設定される。
【００４８】
次に、Ｓ１０３に進み、流れ解析実行部３において、上述の（１）、（２）、（３）式に基づき、各計算セルの流れ解析が実行され、各計算セルの圧力値ｐが演算される。
【００４９】
次いで、Ｓ１０４に進み、圧力空間微分値演算部４において、上述の（４）式に基づき、各計算セルの圧力値ｐの空間微分を演算する。
【００５０】
次に、Ｓ１０５に進み、空隙率更新部５において、各計算セルの空間微分値（∂ｐi／∂ｎ）の大きい方から計算セル容積の和を算出し、計算セル容積の和が制約ボリュームと等しくなった時点の空間微分値（∂ｐi／∂ｎ）（＝（∂ｐ／∂ｎ）ε0.5とする）を、空隙率＝０．５に更新する。
【００５１】
そして、空間微分値（∂ｐi／∂ｎ）が、（∂ｐ／∂ｎ）ε0.5より大きい計算セルの空隙率εを０．５以上に、小さい計算セルの空隙率εを０．５以下に更新し、これを繰り返し計算することで空隙率εの分布を更新していく。最終的に得られた空隙率分布から、ε＝０．５を物体表面形状（ε＞０．５を流体領域、ε＜０．５を物体領域）と定義する。
【００５２】
次いで、Ｓ１０６に進み、計算セル細分割判断実行部６において、計算セルを細分割するか否か判定する。この判定は、計算セルの空隙率が、予め細分割する条件、すなわち、予め設定しておいた初期条件の範囲（εmin〜εmax）にあるか否かで行い、計算セルが、初期条件の範囲にある場合には細分割を行うと判定し、初期条件の範囲にない場合には細分割を行わないと判定する。
【００５３】
尚、前述の如く、計算セル細分割判断実行部６では、空隙率が、予め設定した空隙率（例えば０．５）でポスト処理される等値面の形状変化率が、予め設定する値を超える際に、この等値面が存在する計算セルのみを細分割するように判定するものであっても良い。
【００５４】
そして、上述のＳ１０６の判定の結果、計算セルを細分割すると判定した場合はＳ１０７に進み、計算セルの細分割を実行し、Ｓ１０８へと進む。一方、Ｓ１０６の判定の結果、計算セルを細分割しないと判定した場合は直接Ｓ１０８へと進む。
【００５５】
Ｓ１０８では、初期入力した演算の更新回数に達したか否か判定し、更新回数に達していない場合は再びＳ１０３へと戻り演算を繰り返す。
【００５６】
一方、Ｓ１０８で演算の更新回数に達したと判定した場合は、Ｓ１０９に進み、境界部位形成部８において、ポスト処理にて空隙率０．５の等値面で最終的な境界部位、すなわち、表面形状を形成する。
【００５７】
その後、Ｓ１１０に進み、データ出力部９により、境界部位形成部８からのデータを、モニタ等に出力させる。また、必要に応じて、等値面を構成する面群データ或いは点群データからＣＡＤ（ComputerAided Design）のSurfaceデータ又はVolumeデータを作成しておき、ハードディスク等のデータ保存部１０に、ＩＧＥＳ形式、ＳＴＥＰ形式、或いは、ＳＴＬ形式で保存する。
【００５８】
このように、本発明の実施の形態により、実際に形成した表面形状の結果を図１０、及び、図１１に示す。
【００５９】
図１０は、２次元四角形セルによる流体解析結果の説明図であり、表面形状が存在する部分の計算セルが細分化され、ダクトの形状が滑らかに再現することができている。
【００６０】
また、図１１は２次元三角形セルによる流体解析結果の説明図であり、この図１１も図１０と同様に、表面形状が存在する部分の計算セルが細分化され、ダクトの形状が滑らかに再現することができている。
【００６１】
そして、図１０、図１１の図からも分かるように、本発明は、計算セルの大きさや、形状に影響されることがなく、汎用性に優れ、最適なダクト形状を得ることができる。また、計算を行う領域も、ダクト形状が存在する領域に限られるため、演算量が特に大きく増加することも無い。
【００６２】
【発明の効果】
以上、説明したように本発明によれば、汎用性に優れ、また、演算負荷が必要以上に増加することなく、且つ、精度良く滑らかな物体形状を解析し得ることが可能となる。
【図面の簡単な説明】
【図１】本発明の実施の一形態に係る、流体解析装置の機能ブロック図
【図２】同上、検討空間の初期計算セルの説明図
【図３】同上、細分割処理に係る２つの領域についての説明図
【図４】同上、図３の２つの領域の組み合わせについての説明図
【図５】同上、細分割する計算セルの形状の説明図
【図６】同上、計算セルの細分割処理の世代の説明図
【図７】同上、細分割処理を実行する条件の他の例の説明図
【図８】同上、細分割処理により新たに細分割処理が必要となる場合の説明図
【図９】同上、流体解析プログラムのフローチャート
【図１０】同上、２次元四角形セルによる流体解析結果の説明図
【図１１】同上、２次元三角形セルによる流体解析結果の説明図
【図１２】従来のダクトレイアウト最適化の手法を説明する、多孔質体の説明図
【図１３】同上、検討空間の説明図
【図１４】同上、当初より検討空間に設定される計算セルの説明図
【図１５】同上、検討空間内のダクト容積の説明図
【図１６】同上、図１４の各計算セルにて得られた物体形状の結果の説明図
【符号の説明】
１ 流体解析装置
２ 初期空隙率設定部
３ 流れ解析実行部
４ 圧力空間微分値演算部
５ 空隙率更新部
６ 計算セル細分割判断実行部
７ 更新判断部
８ 境界部位形成部
９ データ出力部
１０ データ保存部
【発明の属する技術分野】
本発明は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、上記設計領域を複数の演算領域に細分割して流体の解析を行う流体解析方法、及び、その流体解析方法を用いた流体解析装置に関する。
【０００２】
【従来の技術】
従来より、空調ダクト等の設計においては、限られたスペース内で配風要件を満足させ、しかも圧力損失の少ないダクト形状を検討することが必要である。このようなことから、本発明者は、社団法人自動車技術会 学術講演会前刷集９６４ Ｎｏ．９６３６８９９（１９９６年１０月）において、圧力損失の少ないダクトを最適ダクトと定義し、限られたダクト容積で出入り口間を連結する概念形状を求めるダクトレイアウト最適化の手法を提案している。
【０００３】
このダクトレイアウト最適化の手法では、設計領域として、図１２に示すように、内部に小さな固体が無数に散らばっている多孔質体を考え、この多孔質体の空隙率に相当する流動抵抗が流体に作用するものとする。そして、例えば図１３に示すような検討空間（設計領域）を、図１４（ａ）、図１４（ｂ）、或いは、図１４（ｃ）に示すような、複数の演算領域に細分割（複数の計算セルに細分割）して考え、設計領域内に占める物体の容積を一定に保ちながら（図１５におけるダクト容積：これが制約条件となる）、多孔質体の空隙率をコントロールすることで圧力損失の少ない概念形状を求める手法である。すなわち、圧力勾配の大きい、流れやすい計算セルの空隙率を大きく、逆に圧力勾配の小さい、あまり流れのない計算セルの空隙率を小さくコントロールする。こうして得られる結果から、例えば、空隙率が０．５の等値面を境界部位としてポスト処理することで物体形状（概念形状）を求めるようになっている。
【０００４】
【発明が解決しようとする課題】
しかしながら、上述のダクトレイアウト最適化の手法では、計算当初に設定しておいた計算セルによっては、図１４（ａ）に示すように計算セルのサイズが粗い場合には、図１６（ａ）に示すように物体形状が滑らかに表現することができないという問題があり、逆に図１４（ｂ）に示すように各計算セルのサイズが細かすぎると演算負荷が大きくなるという問題がある。また、図１４（ｃ）に示すように各計算セルの密度が不均一な場合には、形状表現にムラを生じてしまい、図１６（ｃ）に示すように、精度の良い滑らかな物体形状が得られないという問題がある。精度の良い滑らかな物体形状を得るために、予め必要な部位の計算セルのサイズを細かくしておくことも考えられるが、予め必要な部位を予想して他の部位と調整しながら細かく分割しておくことは困難で、このような作業を行うことは解析方法の汎用性を低下させる一因となってしまう。
【０００５】
本発明は上記事情に鑑みなされたもので、汎用性に優れ、また、演算負荷が必要以上に増加することなく、且つ、精度良く滑らかな物体形状を解析し得ることが可能な流体解析方法、及び、その流体解析方法を用いた流体解析装置を提供することを目的とする。
【０００６】
【課題を解決するための手段】
上記目的を達成するため、本発明の請求項１記載の流体解析方法は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、上記設計領域を複数の演算領域に細分割して流体の解析を行う流体解析方法であって、上記解析対象とする設計領域には、更なる細分割を実行する領域と更なる細分割を禁止する領域の少なくともどちらかを予め設定自在で、上記複数の演算領域の細分割を変更自在とすることを特徴とする。
【０００７】
すなわち、請求項１記載の流体解析方法は、解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、設計領域を複数の演算領域に細分割し、更に、必要に応じて複数の演算領域の細分割を変更して、流体の解析を行う。このため、様々な設計領域で応用が可能で汎用性に優れ、また、演算領域の細分割の変更は、精度良く滑らかな物体形状を解析する部位にのみ行えば良いため、演算負荷が必要以上に増加することがない。また、解析対象とする設計領域には、更なる細分割を実行する領域と更なる細分割を禁止する領域の少なくともどちらかを予め設定自在であり、ユーザがより深く解析したい部位や、或いは、解析を必要としない部位における演算の節約を図れるように満足して、ユーザにとっての利便性をより向上させる。
【０００８】
また、本発明の請求項２記載の流体解析方法は、請求項１記載の流体解析方法において、上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順を有し、上記所定のパラメータが予め設定しておいた範囲の値の演算領域を更に細分割することを特徴とする。
【０００９】
すなわち、請求項２記載の流体解析方法は、請求項１記載の流体解析方法において、複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算し、所定のパラメータが予め設定しておいた範囲の値の演算領域を更に細分割する。
【００１０】
更に、本発明の請求項３記載の流体解析方法は、請求項１又は請求項２に記載の流体解析方法において、上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順と、上記所定のパラメータが予め設定しておいた値の部位を境界部位として連続させる手順とを有し、上記境界部位の形状の変化量が予め設定する値を超える際に、該境界部位が存在する演算領域を更に細分割することを特徴とする。
【００１１】
すなわち、請求項３記載の流体解析方法は、請求項１又は請求項２に記載の流体解析方法において、複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算し、所定のパラメータが予め設定しておいた値の部位を境界部位として連続させる。そして、境界部位の形状の変化量が予め設定する値を超える際に、境界部位が存在する演算領域を更に細分割する。
【００１４】
更に、本発明の請求項４記載の流体解析方法は、請求項１乃至請求項３の何れか一つに記載の流体解析方法において、上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、元の領域の細分割を実行して形成した領域を元の領域に対する１世代後の子の領域とし、該子の領域の細分割を実行して形成した領域を上記元の領域に対する２世代後の孫の領域として複数世代で形成する場合、２世代以上離れて形成した領域間の隣接を禁止することを特徴とする。
【００１５】
すなわち、請求項４記載の流体解析方法は、請求項１乃至請求項３の何れか一つに記載の流体解析方法において、細分割を繰り返すことによる各演算領域間の接続が不自然となることを防止する。
【００１６】
また、本発明の請求項５記載の流体解析方法は、請求項１乃至請求項４の何れか一つに記載の流体解析方法において、上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、上記細分割した後に上記細分割の対象外となった演算領域は、上記細分割の前の大きさの演算領域に復元自在であることを特徴とする。
【００１７】
すなわち、請求項５記載の流体解析方法は、請求項１乃至請求項４の何れか一つに記載の流体解析方法において、細分割を繰り返すことにより、不要なデータが多く蓄積されることを未然に防止できるようにする。
【００１８】
更に、本発明の請求項６記載の流体解析装置は、上記請求項１乃至請求項５の何れか一つに記載の流体解析方法を備え、流体の解析を行うことを特徴とする。
【００１９】
【発明の実施の形態】
以下、図面に基づいて本発明の実施の形態を説明する。図１〜図１１は本発明の実施の一形態に係わり、図１は流体解析装置の機能ブロック図、図２は検討空間の初期計算セルの説明図、図３は細分割処理に係る２つの領域についての説明図、図４は図３の２つの領域の組み合わせについての説明図、図５は細分割する計算セルの形状の説明図、図６は計算セルの細分割処理の世代の説明図、図７は細分割処理を実行する条件の他の例の説明図、図８は細分割処理により新たに細分割処理が必要となる場合の説明図、図９は流体解析プログラムのフローチャート、図１０は２次元四角形セルによる流体解析結果の説明図、図１１は２次元三角形セルによる流体解析結果の説明図である。
【００２０】
図１において、符号１はマイクロコンピュータとその周辺回路を備えて構成された流体解析装置を示し、この流体解析装置１は、初期空隙率設定部２、流れ解析実行部３、圧力空間微分値演算部４、空隙率更新部５、計算セル細分割判断実行部６、更新判断部７、境界部位形成部８、データ出力部９、データ保存部１０から主要に構成されている。尚、本発明の実施の形態においては、検討空間（設計領域）として、例えば図１３に示す空間、すなわち、左側の壁に垂直吹き出しの入口を有し、右側の壁に自由流出の出口を備えた空間を例に説明する。
【００２１】
初期空隙率設定部２は、予め設定された後述の制約ボリュームが入力され、この制約ボリュームが初期空隙率として設定されて、流れ解析実行部３に出力される。制約ボリュームは、図１５に示すように、検討空間内において予め指定された、占有することのできるダクト容積を検討空間の全容積で除した（＝（ダクト容積）／（検討空間容積））ものである。
【００２２】
流れ解析実行部３は、検討空間の境界条件や、初期計算セルの形状等が入力され、以下の多孔質体内における流れの数値計算を実行する。そして、この数値計算により得られた各計算セルの圧力値ｐは、圧力空間微分値演算部４へと入力される。
・連続の式
∂（ε・ｕi）／∂ｘi＝０ …（１）
・流れの運動方程式
ｕj・（∂（ε・ｕi）／∂ｘj）＝−∂（ε・ｐ）／∂ｘi＋（１／Ｒｅ）・（∂２（ε・ｕi）／∂ｘj２）−Ｒi …（２）
ここで、計算セルの空隙率εを次のように定義する。
【００２３】
ε＝（流体領域の体積）／（計算セル体積） …（３）
ｕ：速度成分、ｘ：座標軸、ｐ：圧力、Ｒｅ：レイノルズ数、Ｒ：流動抵抗（＝ｆ（ε，ｕｉ））、添え字ｉ，ｊは、総和規則に従う（＝１，２，３）
【００２４】
また、初期計算セルの形状は、図２に示すように、２次元における解析では、複数の四角形の計算セル（図２（ａ））、或いは、複数の三角形の計算セル（図２（ｂ））となっている。
【００２５】
更に、図２の実線の閉空間で示すように、検討空間には、例えば長方形（図２（ａ））、楕円（図２（ｂ））、円及びこれらの組み合わせ形状による幾何学的領域が、ユーザにより任意に設定自在になっており、計算セルの中心が、この領域内にある計算セルを、強制的に更に細分割することや、或いは、細分割することを禁止することが可能になっている。これは、例えば、細分割する領域が図３（ａ）の斜線のように得られ、強制的に更に細分割する領域を図３（ｂ）の斜線のように設定しておいたものとすると、結果として細分割される領域は図４（ａ）に示す斜線領域となる。一方、細分割する領域が図３（ａ）の斜線のように得られ、細分割することを禁止する領域を図３（ｂ）の斜線のように設定しておいたものとすると、結果として細分割される領域は図４（ｂ）に示す斜線領域となる。このような、細分割を選択可能な領域を設けることにより、ユーザが、より高度で緻密な解析を行いたい部位での利便性を図ることができ、また、解析にあまり必要としないデータ部分での演算の削減を行うことが可能となっている。
【００２６】
尚、図示はしないが、３次元における解析では、初期計算セルの形状は、複数の６面体の計算セル、或いは、複数の四面体の計算セル等である。この場合においても、検討空間には、例えば円柱や、球及びこれらの組み合わせ形状による幾何学的領域が、ユーザにより任意に設定自在になっており、計算セルの中心が、この領域内にある計算セルを、強制的に更に細分割することや、或いは、細分割することを禁止することが可能になっている。
【００２７】
圧力空間微分値演算部４は、流れ解析実行部３から各計算セルの圧力値ｐが入力され、これら圧力値ｐの空間微分を演算し、空隙率更新部５に出力する。
すなわち、
∂ｐi／∂ｎ＝（Σ（∂ｐi／∂ｘi）２）（１／２） …（４）
【００２８】
空隙率更新部５は、圧力空間微分値演算部４から各計算セルの空間微分値（∂ｐi／∂ｎ）が入力され、空間微分値（∂ｐi／∂ｎ）の大きい方から計算セル容積の和を算出し、計算セル容積の和が制約ボリュームと等しくなった時点の空間微分値（∂ｐi／∂ｎ）（＝（∂ｐ／∂ｎ）ε0.5とする）を、空隙率＝０．５に更新する。
【００２９】
そして、空間微分値（∂ｐi／∂ｎ）が、（∂ｐ／∂ｎ）ε0.5より大きい計算セルの空隙率εを０．５以上に、小さい計算セルの空隙率εを０．５以下に更新し、これを繰り返し計算することで空隙率εの分布を更新していく。最終的に得られた空隙率分布から、ε＝０．５を物体表面形状（ε＞０．５を流体領域、ε＜０．５を物体領域）と定義するのである。尚、本実施の形態の例では、空隙率εが０．５で設定しているが、制約ボリューム等の条件によっては、他の値、例えば、０．３や、０．２として演算しても良い。
【００３０】
このように、空隙率更新部５で更新される各計算セルの空隙率の分布状態は、計算セル細分割判断実行部６へと出力される。
【００３１】
計算セル細分割判断実行部６は、空隙率更新部５から各計算セルの空隙率の分布状態が入力され、また、細分割する条件（初期条件：例えば、細分割を行う計算セルの空隙率の範囲εmin〜εmax）が入力されて、各計算セルの細分割の実行可否を判断し、細分割する計算セルを抽出して細分割処理し、細分割処理を終えたデータを更新判断部７に出力する。
【００３２】
ここで、計算セルの細分割は、例えば図５に示すように、計算セルの形状によって異なっている。例えば、計算セルが２次元の四角形であれば、図５（ａ）に示すように、細分割する計算セルは、四角形を４等分するものであり、これを更に細分割する計算セルも、細分割した計算セルを４等分するものとなっている。
【００３３】
また、計算セルが２次元の三角形であれば、図５（ｂ）に示すように、細分割する計算セルは、三角形を４等分するものであり、これを更に細分割する計算セルも、細分割した計算セルを４等分するものとなっている。
【００３４】
更に、計算セルが３次元の６面体であれば、図５（ｃ）に示すように、細分割する計算セルは、６面体を８等分するものであり、これを更に細分割する計算セルも、図示はしないが、細分割した計算セルを８等分するものとなっている。
【００３５】
また、計算セルが３次元の４面体であれば、図５（ｄ）に示すように、細分割する計算セルは、４面体を８等分するものであり、これを更に細分割する計算セルも、図示はしないが、細分割した計算セルを８等分するものとなっている。
【００３６】
細分割する計算セルは、図６に示すように、本実施の形態においては、世代で定義する。すなわち、初期の計算セルは、複数の親セルから構成され、この親セルから細分割されたセルは、親セルに対して１代下の代として子セルとして定義する。同様に、子セルから細分割されたセルは、親セルに対して２代下の代として孫セルとして定義する。更に、孫セルから細分割されたセルは、親セルに対して３代下の代としてひ孫セルとして定義する。
【００３７】
計算セル細分割判断実行部６では、空隙率が、予め細分割する条件、すなわち、予め設定しておいた初期条件の範囲にある計算セルを抽出し、細分割処理を行うが、以下のような条件としても良い。
【００３８】
空隙率が、予め設定した空隙率（例えば０．５）でポスト処理される境界部位（以降、等値面という）の形状変化率が、予め設定する値を超える際に、この等値面が存在する計算セルのみを細分割する。
【００３９】
すなわち、図７（ａ）に示すように、まず、ポスト処理にて空隙率が０．５の等値面を形成する。次いで、これら等値面の各計算セル毎の法線ベクトルｎ１〜ｎ１１を算出し、これら法線ベクトルｎ１〜ｎ１１の隣接するベクトル同士のなす角が、予め設定しておいた角度を超えている場合、その法線ベクトルの等値面が所属する計算セルを、図７（ｂ）に示すように細分割処理する。そして、これを繰り返していくことで、形状変化の大きな部分が滑らかに表現できるようになっている。尚、３次元においても、図７（ｃ）に示すように、上述の２次元の細分割処理と同様に、ポスト処理にて形成した空隙率０．５の各等値面に法線ベクトルｍ１〜ｍ４を算出し、これら各法線ベクトルの隣接する法線ベクトル同士のなす角を求め、設定値と比較することにより、等値面が所属する計算セルを細分割処理することで行うことができる。
【００４０】
また、計算セル細分割判断実行部６では、細分割する計算セルは、滑らかな連続性を維持するために、２世代以上離れた計算セルと隣接しないようにする。特定の計算セルの細分割により、２世代以上離れるような場合は、図８に示すように、例えば、子セルの細分割により新たな孫セルが親セルと隣接するような場合は、この孫セルと隣接する親セルも子セルに細分割する。
【００４１】
また、計算セル細分割判断実行部６では、計算セルの細分割を繰り返し、過去に細分割した計算セルにおいて、現在の物体形状表現に使われていない（計算に関与していない）計算セルは、適宜、１世代前の計算セルに戻すことが可能になっている。このようにすることで、全体の計算セル数の削減を図り、データ数を削減することが可能になっている。
【００４２】
更新判断部７は、当初入力された演算の更新回数（例えば、１００回等）に達しているか否か判定し、演算回数が更新回数に達している場合には、演算により得られたデータを境界部位形成部８に出力する。
【００４３】
境界部位形成部８では、入力された各計算セルの空隙率のデータに対し、予め設定しておいた物体表面形状とみなす空隙率、すなわち、０．５の等値面を表示させる。
【００４４】
データ出力部９では、境界部位形成部８からのデータを、図示しないモニタに出力させる。また、必要に応じて、等値面を構成する面群データ或いは点群データからＣＡＤ（Computer Aided Design）のSurfaceデータ又はVolumeデータを作成しておき、ハードディスク等のデータ保存部１０に、ＩＧＥＳ形式、ＳＴＥＰ形式、或いは、ＳＴＬ形式等のＣＡＤで多く用いられるデータ形式で保存する。
【００４５】
そして、ダクト等の設計の際に、ハードディスク等のデータ保存部１０から、直接、保存しておいた最適な物体表面形状のデータを読み込むことにより、流体解析の設計へのフィードバックが迅速に行われるようになっている。
【００４６】
次に、上記構成の流体解析装置１で実行される流体解析について、図９のフローチャートで説明する。
まず、ステップ（以下、「Ｓ」と略称）１０１では、解析に必要なデータの入力が行われる。これは、検討空間の形状や入口・出口の境界条件、初期計算セルの形状、ダクトの制約ボリューム、演算の更新回数、細分割する条件（細分割する計算セルの空隙率の範囲）等である。
【００４７】
次いで、Ｓ１０２では、初期空隙率設定部２において、初期空隙率の設定を行う。本実施の形態において、この初期空隙率は、制約ボリュームが初期空隙率として設定される。
【００４８】
次に、Ｓ１０３に進み、流れ解析実行部３において、上述の（１）、（２）、（３）式に基づき、各計算セルの流れ解析が実行され、各計算セルの圧力値ｐが演算される。
【００４９】
次いで、Ｓ１０４に進み、圧力空間微分値演算部４において、上述の（４）式に基づき、各計算セルの圧力値ｐの空間微分を演算する。
【００５０】
次に、Ｓ１０５に進み、空隙率更新部５において、各計算セルの空間微分値（∂ｐi／∂ｎ）の大きい方から計算セル容積の和を算出し、計算セル容積の和が制約ボリュームと等しくなった時点の空間微分値（∂ｐi／∂ｎ）（＝（∂ｐ／∂ｎ）ε0.5とする）を、空隙率＝０．５に更新する。
【００５１】
そして、空間微分値（∂ｐi／∂ｎ）が、（∂ｐ／∂ｎ）ε0.5より大きい計算セルの空隙率εを０．５以上に、小さい計算セルの空隙率εを０．５以下に更新し、これを繰り返し計算することで空隙率εの分布を更新していく。最終的に得られた空隙率分布から、ε＝０．５を物体表面形状（ε＞０．５を流体領域、ε＜０．５を物体領域）と定義する。
【００５２】
次いで、Ｓ１０６に進み、計算セル細分割判断実行部６において、計算セルを細分割するか否か判定する。この判定は、計算セルの空隙率が、予め細分割する条件、すなわち、予め設定しておいた初期条件の範囲（εmin〜εmax）にあるか否かで行い、計算セルが、初期条件の範囲にある場合には細分割を行うと判定し、初期条件の範囲にない場合には細分割を行わないと判定する。
【００５３】
尚、前述の如く、計算セル細分割判断実行部６では、空隙率が、予め設定した空隙率（例えば０．５）でポスト処理される等値面の形状変化率が、予め設定する値を超える際に、この等値面が存在する計算セルのみを細分割するように判定するものであっても良い。
【００５４】
そして、上述のＳ１０６の判定の結果、計算セルを細分割すると判定した場合はＳ１０７に進み、計算セルの細分割を実行し、Ｓ１０８へと進む。一方、Ｓ１０６の判定の結果、計算セルを細分割しないと判定した場合は直接Ｓ１０８へと進む。
【００５５】
Ｓ１０８では、初期入力した演算の更新回数に達したか否か判定し、更新回数に達していない場合は再びＳ１０３へと戻り演算を繰り返す。
【００５６】
一方、Ｓ１０８で演算の更新回数に達したと判定した場合は、Ｓ１０９に進み、境界部位形成部８において、ポスト処理にて空隙率０．５の等値面で最終的な境界部位、すなわち、表面形状を形成する。
【００５７】
その後、Ｓ１１０に進み、データ出力部９により、境界部位形成部８からのデータを、モニタ等に出力させる。また、必要に応じて、等値面を構成する面群データ或いは点群データからＣＡＤ（ComputerAided Design）のSurfaceデータ又はVolumeデータを作成しておき、ハードディスク等のデータ保存部１０に、ＩＧＥＳ形式、ＳＴＥＰ形式、或いは、ＳＴＬ形式で保存する。
【００５８】
このように、本発明の実施の形態により、実際に形成した表面形状の結果を図１０、及び、図１１に示す。
【００５９】
図１０は、２次元四角形セルによる流体解析結果の説明図であり、表面形状が存在する部分の計算セルが細分化され、ダクトの形状が滑らかに再現することができている。
【００６０】
また、図１１は２次元三角形セルによる流体解析結果の説明図であり、この図１１も図１０と同様に、表面形状が存在する部分の計算セルが細分化され、ダクトの形状が滑らかに再現することができている。
【００６１】
そして、図１０、図１１の図からも分かるように、本発明は、計算セルの大きさや、形状に影響されることがなく、汎用性に優れ、最適なダクト形状を得ることができる。また、計算を行う領域も、ダクト形状が存在する領域に限られるため、演算量が特に大きく増加することも無い。
【００６２】
【発明の効果】
以上、説明したように本発明によれば、汎用性に優れ、また、演算負荷が必要以上に増加することなく、且つ、精度良く滑らかな物体形状を解析し得ることが可能となる。
【図面の簡単な説明】
【図１】本発明の実施の一形態に係る、流体解析装置の機能ブロック図
【図２】同上、検討空間の初期計算セルの説明図
【図３】同上、細分割処理に係る２つの領域についての説明図
【図４】同上、図３の２つの領域の組み合わせについての説明図
【図５】同上、細分割する計算セルの形状の説明図
【図６】同上、計算セルの細分割処理の世代の説明図
【図７】同上、細分割処理を実行する条件の他の例の説明図
【図８】同上、細分割処理により新たに細分割処理が必要となる場合の説明図
【図９】同上、流体解析プログラムのフローチャート
【図１０】同上、２次元四角形セルによる流体解析結果の説明図
【図１１】同上、２次元三角形セルによる流体解析結果の説明図
【図１２】従来のダクトレイアウト最適化の手法を説明する、多孔質体の説明図
【図１３】同上、検討空間の説明図
【図１４】同上、当初より検討空間に設定される計算セルの説明図
【図１５】同上、検討空間内のダクト容積の説明図
【図１６】同上、図１４の各計算セルにて得られた物体形状の結果の説明図
【符号の説明】
１ 流体解析装置
２ 初期空隙率設定部
３ 流れ解析実行部
４ 圧力空間微分値演算部
５ 空隙率更新部
６ 計算セル細分割判断実行部
７ 更新判断部
８ 境界部位形成部
９ データ出力部
１０ データ保存部
Claims (6)
- 解析対象とする設計領域内を小さな固体が無数に散らばっている多孔質体と仮定し、上記設計領域を複数の演算領域に細分割して流体の解析を行う流体解析方法であって、
上記解析対象とする設計領域には、更なる細分割を実行する領域と更なる細分割を禁止する領域の少なくともどちらかを予め設定自在で、
上記複数の演算領域の細分割を変更自在とすることを特徴とする流体解析方法。 - 上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順を有し、
上記所定のパラメータが予め設定しておいた範囲の値の演算領域を更に細分割することを特徴とする請求項１記載の流体解析方法。 - 上記複数の演算領域毎に、連続の式と流れの運動方程式に基づき所定のパラメータを演算する手順と、
上記所定のパラメータが予め設定しておいた値の部位を境界部位として連続させる手順とを有し、
上記境界部位の形状の変化量が予め設定する値を超える際に、該境界部位が存在する演算領域を更に細分割することを特徴とする請求項１又は請求項２に記載の流体解析方法。 - 上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、元の領域の細分割を実行して形成した領域を元の領域に対する１世代後の子の領域とし、該子の領域の細分割を実行して形成した領域を上記元の領域に対する２世代後の孫の領域として複数世代で形成する場合、２世代以上離れて形成した領域間の隣接を禁止することを特徴とする請求項１乃至請求項３の何れか一つに記載の流体解析方法。
- 上記複数の演算領域の細分割の変更で、更なる細分割を繰り返して実行する際は、上記細分割した後に上記細分割の対象外となった演算領域は、上記細分割の前の大きさの演算領域に復元自在であることを特徴とする請求項１乃至請求項４の何れか一つに記載の流体解析方法。
- 上記請求項１乃至請求項５の何れか一つに記載の流体解析方法を備え、流体の解析を行うことを特徴とする流体解析装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001377582A JP4005803B2 (ja) | 2001-12-11 | 2001-12-11 | 流体解析方法、及び、その流体解析方法を用いた流体解析装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001377582A JP4005803B2 (ja) | 2001-12-11 | 2001-12-11 | 流体解析方法、及び、その流体解析方法を用いた流体解析装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003178100A JP2003178100A (ja) | 2003-06-27 |
| JP4005803B2 true JP4005803B2 (ja) | 2007-11-14 |
Family
ID=19185505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001377582A Expired - Fee Related JP4005803B2 (ja) | 2001-12-11 | 2001-12-11 | 流体解析方法、及び、その流体解析方法を用いた流体解析装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4005803B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4804189B2 (ja) * | 2006-03-29 | 2011-11-02 | 富士通株式会社 | 構造解析装置、構造解析方法、及び構造解析プログラム |
| US7921002B2 (en) * | 2007-01-04 | 2011-04-05 | Honda Motor Co., Ltd. | Method and system for simulating flow of fluid around a body |
| JP5324820B2 (ja) * | 2008-05-19 | 2013-10-23 | 住友ゴム工業株式会社 | 解析モデルの作成方法 |
| JP5782069B2 (ja) * | 2013-06-25 | 2015-09-24 | 大阪瓦斯株式会社 | 貯湯水温度分布計算方法及び貯湯式熱源装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0363880A (ja) * | 1989-08-02 | 1991-03-19 | Hitachi Ltd | オブジェクトモデリングシステム |
| JPH04679A (ja) * | 1990-04-18 | 1992-01-06 | Hitachi Ltd | 座標格子作成支援方法及びその装置 |
| JP2642070B2 (ja) * | 1994-11-07 | 1997-08-20 | インターナショナル・ビジネス・マシーンズ・コーポレイション | 四角形メッシュの生成方法及びシステム |
| JP2881389B2 (ja) * | 1995-05-16 | 1999-04-12 | 日本アイ・ビー・エム株式会社 | 自動メッシュ生成方法及びシステム |
| JPH08320947A (ja) * | 1995-05-25 | 1996-12-03 | Matsushita Electric Ind Co Ltd | 数値解析用メッシュ作成方法及び装置 |
-
2001
- 2001-12-11 JP JP2001377582A patent/JP4005803B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003178100A (ja) | 2003-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Williams et al. | Numerical simulations of laminar flow over a 3D backward‐facing step | |
| CN111125942B (zh) | 用于三维单元结构建模和拓扑优化的b样条高清晰度单元水平集方法和计算机存储介质 | |
| EP1840841B1 (en) | Evolutionary direct manipulation of free form deformation representations for design optimisation | |
| Gokce et al. | Branch and bound search to optimize injection gate locations in liquid composite molding processes | |
| US20070229544A1 (en) | Nurbs surface deformation apparatus and the method using 3d target curve | |
| Yaman et al. | Slice coherence in a query-based architecture for 3D heterogeneous printing | |
| US10121279B1 (en) | Systems and methods for generating a mesh | |
| EP1675068B1 (en) | Evolutionary optimisation and free form deformation | |
| Zhang et al. | A prototype mesh generation tool for CFD simulations in architecture domain | |
| CN115994410B (zh) | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 | |
| CN110852003A (zh) | 计算机辅助设计定义的几何形状中的对象间的间隙的检测 | |
| JP4005803B2 (ja) | 流体解析方法、及び、その流体解析方法を用いた流体解析装置 | |
| JPH08129636A (ja) | 地形モデル作成装置および方法 | |
| Zeng et al. | Modeling method for variable and isotropic permeability design of porous material based on TPMS lattices | |
| Letov et al. | A geometric modelling framework to support the design of heterogeneous lattice structures with non-linearly varying geometry | |
| Wang et al. | Topology optimization of self-supporting lattice structure | |
| EP1840842A1 (en) | Evolutionary design optimisation by means of extended direct manipulation of free form deformations | |
| CN115203997A (zh) | 一种基于多变量设计的点阵-实体复合结构拓扑优化方法 | |
| Gifford | Compensating for die swell in the design of profile dies | |
| Azman et al. | Creation of lattice structures with skeleton model for additive manufacturing | |
| CN113887095A (zh) | 一种基于等几何分析的渐进式结构拓扑优化方法 | |
| JP4639292B2 (ja) | ３次元メッシュ生成方法 | |
| JP2000340476A (ja) | 形状シミュレーション方法、装置および記録媒体 | |
| Li et al. | Topology optimization of multi-phase shell-infill composite structure for additive manufacturing | |
| JP3770991B2 (ja) | 解析モデルの生成方法および装置ならびに射出成形過程の解析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041124 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070605 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070709 |
|
| 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: 20070821 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070824 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100831 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110831 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110831 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120831 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120831 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130831 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |