(57)【要約】【課題】実際の浸漬工程に近い状態での空気溜りの評価を可能とする。【解決手段】対象物の表面形状のデータを複数の要素に分割して数値計算のための解析モデルを構築し（Ｓ１）、重力方向、初期条件を設定すると（Ｓ２，Ｓ３）、エアポケット判定計算処理を実施し、気液境界線を検出する（Ｓ４）。次に、気液境界線の高さＨに修正高さΔＨを加算して境界線を修正し（Ｓ５）、ポスト処理を実行する（Ｓ６）。このポスト処理では、シミュレーション結果を容易に把握可能とするための可視化処理を行い、エアポケットの高さやエアポケットの高さと境界面積との比に応じて序列付けや数値化を行い、これらの序列あるいは数値に応じてエアポケットの色相を変化させる等して識別表示する。これにより、実際の浸漬工程で生じるであろうエアポケットに対する対策の難易度を知ることができる。【選択図】図１０
2008049269.tif 000002
【特許請求の範囲】
【請求項１】
浸漬処理を施すワークの形状を複数の要素のデータでモデル化し、各要素の属性が気体であるか液体であるかを決定して空気溜りをシミュレーションする空気溜りのシミュレーション方法において、
属性が気体であると決定された要素と属性が液体であると決定された要素との境界線を、上記ワークの実際の浸漬条件に基づいて修正する境界線修正ステップと、
上記修正後の境界線に基づいて空気溜りを評価した評価データを算出し、該評価データに基づいて上記空気溜りを可視化処理するポスト処理ステップと
を備えることを特徴とする空気溜まりのシミュレーション方法。
【請求項２】
上記境界線修正ステップにおいては、
上記境界線を形成する上記複数の要素の高さを平滑化して直線化した上で、上記浸漬条件としてワークの揺れに基づく修正高さを加えるにより、上記境界線を修正することを特徴とする請求項１記載の空気溜りのシミュレーション方法。
【請求項３】
上記境界線修正ステップにおいては、
上記修正高さを、上記ワークの釣下げ部における揺れ角度と上記ワークの釣下げ部から端部までの水平方向距離との相関関係について予め求めた実験式に基づいて算出することを特徴とする請求項２記載の空気溜りのシミュレーション方法。
【請求項４】
上記ポスト処理ステップにおいては、
上記評価データを、上記空気溜りの深さと面積との少なくとも一方に基づいて算出することを特徴とする請求項１〜３の何れか一に記載の空気溜りのシミュレーション方法。
【請求項５】
上記ポスト処理ステップにおいては、
上記空気溜りの可視化処理を、上記評価データに基づく序列或いは数値の大きさに応じて上記空気溜りを識別表示する処理とすることを特徴とする請求項１〜４の何れか一に記載の空気溜りのシミュレーション方法。
【請求項６】
浸漬処理を施す対象物の形状を複数の要素のデータでモデル化し、各要素の属性が気体であるか液体であるかを決定して空気溜りをシミュレーションするコンピュータが実行可能な空気溜りのシミュレーションプログラムにおいて、
属性が気体であると決定された要素と属性が液体であると決定された要素との境界線を、上記ワークの実際の浸漬条件に基づいて修正する境界線修正ステップと、
上記修正後の境界線に基づいて空気溜りを評価した評価データを算出し、該評価データに基づいて上記空気溜りを可視化処理するポスト処理ステップと
を備えることを特徴とする空気溜まりのシミュレーションプログラム。
【請求項７】
上記境界線修正ステップにおいては、
上記境界線を形成する上記複数の要素の高さを平滑化して直線化した上で、上記浸漬条件としてワークの揺れに基づく修正高さを加えるにより、上記境界線を修正することを特徴とする請求項６記載の空気溜りのシミュレーションプログラム。
【請求項８】
上記境界線修正ステップにおいては、
上記修正高さを、上記ワークの釣下げ部における揺れ角度と上記ワークの釣下げ部から端部までの水平方向距離との相関関係について予め求めた実験式に基づいて算出することを特徴とする請求項７記載の空気溜りのシミュレーションプログラム。
【請求項９】
上記ポスト処理ステップにおいては、
上記評価データを、上記空気溜りの深さと面積との少なくとも一方に基づいて算出することを特徴とする請求項６〜８の何れか一に記載の空気溜りのシミュレーションプログラム。
【請求項１０】
上記ポスト処理ステップにおいては、
上記空気溜りの可視化処理を、上記評価データに基づく序列或いは数値の大きさに応じて上記空気溜りを識別表示する処理とすることを特徴とする請求項６〜９の何れか一に記載の空気溜りのシミュレーションプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、浸漬処理を施す際にワークに生じる空気溜りをシミュレーションする空気溜りのシミュレーション方法及びシミュレーションプログラムに関する。
【背景技術】
【０００２】
一般に、電着塗装やメッキ等におけるワークの浸漬工程においては、均一な塗膜やメッキ膜を形成させるために、ワークに空気溜り（エアポケット）を生じさせないことが重要である。このエアポケットは、複雑な形状のワークを浸漬した際に、その隙間にできる凹部、例えば、自動車等の車体であれば、フード内面、ルーフ内面、フロア下面等の凹部等の空間に空気が残留する現象であり、このエアポケットが発生すると、塗装不良やメッキ不良が発生してしまう。
【０００３】
そこで、ワークを浸漬処理する状態をコンピュータを用いてシミュレーションし、エアポケットの発生状況を予測する手法が提案されている。エアポケットのコンピュータシミュレーションとしては、ワークの表面形状を、メッシュやノード等の複数の要素のデータでモデル化し、各要素の属性が気体であるか液体であるかを、要素の特徴点（例えば、重心点、頂点、内心点、外心点、垂心点）同士の重力方向の位置を比較して判定する幾何学的な手法が知られている。
【０００４】
例えば、特許文献１には、ワークの各表面をその頂点座標データに基づいて任意数の三角ポリゴン（三角形の要素）に分割し、各々の三角ポリゴンについて空気と浸漬槽内の液体との状態変化を時々刻々と計算することにより、ワークを入槽角の姿勢にて浸漬槽に模擬浸漬した際にエアポケットとなりうる表面部位であるか否かの判定を行う手法が開示されている。
【特許文献１】特開２０００−５１７５０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
従来の幾何学的な手法によるエアポケットのシミュレーションは、塗料等の浸漬液や被塗装物等のワークが静止した状態での評価を行っている。しかしながら、実際の浸漬工程においては、塗料が流れていたり、被塗装物が揺れていたりする等、幾何学的なシミュレーションでは評価が不十分な点が多い。
【０００６】
例えば、図１３に示す電着塗装されたワーク（カットモデル）１００で検証すると、このワーク１００においては、塗料で塗装された領域Ｒ１に対して、図中の上部の領域Ｒ２でエアポケットの発生により素材の地肌が露出している。すなわち、このワーク１００の例では、図中の左側の穴１０１から上の部分では内部の空気が抜けず、穴１０１の上端より上にエアポケットが発生することが予想されるが、実際には、塗料の流れやワークの揺れ等により、穴より所定高さだけ高い位置にエアポケットが発生し、エアポケットと塗装部との境に気液境界線が現れている。
【０００７】
一方、このワーク１００に対して、従来の幾何学的な手法でシミュレーションを実施した結果は、図１４に示される。図１４のシミュレーション結果では、塗料に浸漬されたメッシュ群Ｍｌとエアポケットに相当するメッシュ群Ｍａ（図中の濃い部分）とが識別表示される。しかしながら、このシミュレーションでは、実際の塗装ラインにおける塗料の流れや被塗装物の揺れ等の影響を考慮していないため、エアポケットの領域Ｒは穴１０１の上端位置まで広がっており、図１３の領域Ｒ２で示される実際のエアポケット形状よりも大きくなってしまう。
【０００８】
このように、従来の幾何学的な手法では、エアポケットが発生していることはわかるが、どこにどの程度のエアポケットが発生しているかを的確に評価することは困難である。すなわち、従来の幾何学的な手法によるエアポケットのシミュレーションでは、実際に発生するエアポケットの対応策を考える上での情報が十分とは言えず、ワークの浸漬条件を変更すればエアポケットを無くせるのか、或いはワークの設計変更による形状変更で対応せざるを得ないのかといった判断まで踏み込むことはできない。
【０００９】
本発明は上記事情に鑑みてなされたもので、実際の浸漬工程に近い状態での空気溜りの評価を可能とする空気溜りのシミュレーション方法及びシミュレーションプログラムを提供することを目的としている。
【課題を解決するための手段】
【００１０】
上記目的を達成するため、本発明による空気溜りのシミュレーション方法は、浸漬処理を施すワークの形状を複数の要素のデータでモデル化し、各要素の属性が気体であるか液体であるかを決定して空気溜りをシミュレーションする空気溜りのシミュレーション方法において、属性が気体であると決定された要素と属性が液体であると決定された要素との境界線を、上記ワークの実際の浸漬条件に基づいて修正する境界線修正ステップと、上記修正後の境界線に基づいて空気溜りを評価した評価データを算出し、該評価データに基づいて上記空気溜りを可視化処理するポスト処理ステップとを備えることを特徴とする。
【００１１】
本発明による空気溜りのシミュレーションプログラムは、浸漬処理を施す対象物の形状を複数の要素のデータでモデル化し、各要素の属性が気体であるか液体であるかを決定して空気溜りをシミュレーションするコンピュータが実行可能な空気溜りのシミュレーションプログラムにおいて、属性が気体であると決定された要素と属性が液体であると決定された要素との境界線を、上記ワークの実際の浸漬条件に基づいて修正する境界線修正ステップと、上記修正後の境界線に基づいて空気溜りを評価した評価データを算出し、該評価データに基づいて上記空気溜りを可視化処理するポスト処理ステップとを備えることを特徴とする。
【発明の効果】
【００１２】
本発明によれば、実際の浸漬工程に近い状態での空気溜りの評価が可能となり、ワークの浸漬条件を変更すればエアポケットを無くせるのか、或いはワークの形状変更で対応せざるを得ないのかといった判断まで行うことが可能となり、単なる目安に止まらない実用的なシミュレーションへと発展させることができる。
【発明を実施するための最良の形態】
【００１３】
以下、図面を参照して本発明の実施の形態を説明する。図１〜図１２は本発明の実施の一形態に係り、図１はシミュレーション装置の基本構成図、図２は車体の塗装ラインの概略説明図、図３はフロアパネルの概略斜視図、図４はフロアパネルを複数の要素に分割した解析モデルを示す説明図、図５はメッシュによる境界線の高さを示す説明図、図６はノードによる境界線の高さを示す説明図、図７はワークの揺れを示す説明図、図８は境界線の変更を示す説明図、図９は境界線の表示とメッシュ分割を示す説明図、図１０はエアポケットシミュレーションプログラムのフローチャート、図１１はエアポケット断面を示す説明図、図１２はエアポケット高さ／境界面積の表示例を示す説明図である。
【００１４】
図１に示すように、本実施の形態におけるシミュレーション装置１は、自動車のボディシェル等のワークに電着塗装やメッキ処理等の浸漬処理を施す際に、ワークに生じる空気溜り（エアポケット）をシミュレーションするものである。このエアポケットのシミュレーション装置１は、マイクロコンピュータやパーソナルコンピュータ等の単一のコンピュータ、或いはネットワークを介して相互に接続される複数のコンピュータを用いて構成される。
【００１５】
以下では、便宜上、シミュレーション装置１を単一のコンピュータで構成する例について説明する。このコンピュータとしてのシミュレーション装置１は、演算装置１０、キーボードやマウス等の入力装置１１、ＣＲＴや液晶ディスプレイ等の表示装置１２、磁気ディスクや光ディスク等の外部記憶装置１３等を備えている。
【００１６】
演算装置１０は、ＣＰＵ、ＲＯＭ及びＲＡＭ等の内部メモリ、入出力インターフェース等を備えており、内部のＲＯＭ、外部記憶装置１３、外部の記憶媒体に記憶させたシミュレーションプログラム、或いは、図示しないネットワークや通信装置を介して外部からロードしたシミュレーションプログラムをＣＰＵで実行し、入力装置１１を介して指示された解析対象のワーク（対象物）を擬似的に浸漬槽内に浸漬させ、その浸漬によって対象物に発生するエアポケットをシミュレーションする。
【００１７】
例えば、本発明によるシミュレーションの適用例として、車体ボディの電着塗装時における電着液及び空気の分布を数値解析し、被塗装物におけるエアポケットの発生を予測するシミュレーションが挙げられる。ここで、車体ボディの塗装ラインについて、図２を用いて簡単に説明する。
【００１８】
図２に示すように、溶接等により複数の車体パネルを互いに接合して構成される自動車の車体ボディ２０は、搬送装置２１のハンガに搭載された状態で塗装ラインにて略水平方向へ搬送される。塗装ラインでは、電着塗装の前処理として、車体パネルに脱脂、水洗、表面調整、皮膜化成、水洗等の処理が施される。
【００１９】
これらの処理の後、車体ボディ２０は電着槽２２に向かって降下し、電着液２３に浸漬された状態で略水平に移動する。この状態で、車体ボディ２０と電着槽２２内の電極（図示せず）に電圧を加えることにより、車体パネルに塗料が析出するようになっている。この後、搬送装置２１により車体ボディ２０は電着槽２２から引き上げられ、水洗により車体パネルに電着せずに付着している電着液２３が除去される。
【００２０】
演算装置１０において実行されるシミュレーションプログラムの機能は、モデル構築部１０ａ、エアポケット判定部１０ｂ、境界線修正部１０ｃ、ポスト処理部１０ｄによって表現することができる。演算装置１０は、対象物の形状を複数の要素のデータでモデル化した解析モデルによる幾何学的手法を用いてエアポケットの発生をシミュレーションし、シミュレーション結果を表示装置１２に出力して表示させる。尚、表示装置１２には、シミュレーション結果のみならずシミュレーション過程を表示させるようにしても良い。
【００２１】
本発明によるエアポケットのシミュレーションは、従来の幾何学的な手法によって求めた気体と液体との気液境界線（以下、単に「境界線」と記載する場合もある）に対して、実際の塗料等の浸漬液の状態や浸漬槽内のワークの状態を考慮した修正を行い、より現実に近いエアポケット形状を求めることができる。以下、本実施の形態におけるエアポケットシミュレーションについて説明する。
【００２２】
モデル構築部１０ａは、対象物の形状データを複数の要素に分割して解析モデルを構築する。解析モデルを構成する手法としては、対象物の表面形状を複数の所定の形状に分割して表現するメッシュを用いる手法や、対象物の表面に複数の節点（ノード）を配置する手法等を用いることができる。
【００２３】
例えば、車体ボディ２０を構成する図３に示すようなフロアパネル２６を例に挙げて説明すると、図３に示すフロアパネル２６を、図４に示すように、フロアパネル２６の表面を複数のメッシュ２９に分割して数値計算のための解析モデルを構築する。以下では、対象物の表面形状を三角形のメッシュで表現したサーフェスモデルを用い、このサーフェスモデルのメッシュデータを三次元空間に配置して解析モデルを構築する例を主として説明するが、解析モデルの要素としてノードを用いる場合においても処理は同様である。
【００２４】
尚、解析モデルの要素としてのメッシュの形状は、三角形の形状に限られず、矩形や多角形の形状であっても良い。個々のメッシュには、例えば、自己を識別するためのメッシュ番号、メッシュの節点（ノード）の所定の基準点に対する座標値（三次元空間のＸＹＺ座標軸における座標値）、メッシュの属性等のデータが付与される。以下では、便宜上、ＸＹＺ直交座標系におけるＺ軸の方向は、重力と反対方向であるものとする。
【００２５】
また、外部記憶装置１３には、エアポケットのシミュレーション処理に必要な各種のデータが格納されている。例えば、外部記憶装置１３は、対象物毎に個別の識別番号（レコード番号）が付された属性レコード群で構成されるデータベースを有し、個々の属性レコードに、各メッシュのメッシュ番号、ノード点座標値、属性データ等が対応付けて記述される。
【００２６】
エアポケット判定部１０ｂは、解析モデルを形成する個々のメッシュに対して、エアポケットの箇所に相当する気体の属性と電着液等の浸漬個所に相当する液体の属性との何れかを設定する。このメッシュの属性は、当初、全メッシュが気体の属性に設定され（浸漬前）、その中で所定の初期条件を満足するメッシュの属性が液体に設定される。そして、気体の属性を有するメッシュに対して、メッシュの属性が液体となるか気体のままであるか（つまり、エアポケットとなるか否か）を隣接メッシュとの比較によって判定し、液体と判定されたとき、そのメッシュの属性を液体に変更し、属性データを書き換える。
【００２７】
尚、隣接メッシュとは、判定対象のメッシュに対して節点或いは辺を共有するメッシュである。
【００２８】
このエアポケットの判定は、基本的に、浸漬液の比重が空気の比重よりも大きいことを前提としている。すなわち、或るメッシュの属性を判定する際に、判定対象のメッシュと比較対象のメッシュとを、Ｚ方向（重力と反対方向）の高さで比較する。そして、比較対象メッシュが判定対象メッシュよりも高く液体の属性を有していれば、相対的に低い判定対象メッシュにも浸漬液が空気を押しのけて浸入し、浸漬液により満たされるものと判断して、判定対象メッシュの属性データを液体とする。これに対して、相対的に高い比較対象メッシュの属性が気体である場合には、判定対象メッシュの属性データは変更せず、比較対象を変更して判定を繰り返す。そして、この判定を繰り返して全メッシュの属性が液体と気体とに分離・決定され、属性が気体と決定されたメッシュの存在によって示されるエアポケットの気液境界線が導き出される。この気液境界線は、複数本導き出される場合もある。
【００２９】
境界線修正部１０ｃは、エアポケットの気液境界線を修正し、実際の浸漬工程に近い状態でエアポケットを評価するための基準とする。すなわち、ワークの形状をメッシュ等の複数の要素に分割して各要素の属性を気体と液体とに分離して得られる気液境界線Ｌは、図５に示すように凹凸がある。従って、各境界線のＺ方向の高さＨをそれぞれ求め、この高さＨを、実際の浸漬工程に対応した修正高さΔＨで修正することにより、気液境界線Ｌを変更してエアポケットを評価するための基準とする。
【００３０】
本形態においては、境界線Ｌの高さＨは、境界線を構成するメッシュの辺（エッジ）の長さにより加重平均して求める。ここで、属性が気体であるメッシュを「液外メッシュ」、属性が液体であるメッシュを「液内メッシュ」とし、境界線を構成する液内メッシュのエッジの数がｎであるものとすると、図５に示すように、液内メッシュのエッジの長さをｌｉ、また、このエッジの中点の高さをｈｉとしたとき、境界線の高さＨは、以下の（１）式で表現することができる。尚、（１）式におけるΣは、ｉ〜ｎについての総和を示すものとする。 
Ｈ＝Σ（ｌｉ×ｈｉ）／Σｌｉ  …（１）
この場合、図６に示すように、解析モデルの要素としてノード（図中の黒丸で示す）を用いる場合には、液体の属性のノード（液内ノード）のうち最も高い位置にあるノードと、隣接する気体の属性のノード（液外ノード）との中間を結んだラインが境界線Ｌとして導き出される。この場合には、エッジの長さｌｉは、境界線と、この境界線に隣接するノードとの距離、エッジの中点の高さｈiは、境界線と、この境界線に隣接するノードとの間の中点の高さとして、（１）式を適用する。
【００３１】
一方、実際の浸漬工程に対応した修正高さΔＨは、予め実際の浸漬条件を考慮して定めた一定値や、気液境界線を形成する液内メッシュの中の最も高い位置等でも良いが、一般的には、塗料等の浸漬液の流れ、ワークの揺れ、ワークのフック部（釣下げ部）から端部までの距離、ワークの内部か外部か等の条件によって変化する。従って、本形態においては、これらの条件と修正高さΔＨとの相関関係について予め基礎実験を行い、実験式を取得する。
【００３２】
例えば、図７に示すように、或るエアポケットＡＰにおいて、フックＨＫからの水平方向距離が最も長い境界線位置までの距離（端部までの距離）Ｌと修正高さΔＨとの関係は、ワークの揺れ角度をθとして、幾何学的に以下の(２)式のようになる。但し、(２)式においては、揺れ角度θが微小な角度（略１°程度）であるため、フックＨＫからの垂直方向距離を無視している。  
ΔＨ＝Ｌ×ｔａｎ２θ  …（２）
【００３３】
次に、以下の（３）式に示すように、境界線Ｌの高さＨに修正高さΔＨを加算し、新たな境界線Ｌ’の高さＨｓとして修正する。そして、図８に示すように、高さＨｓ以下の液外メッシュを液内メッシュに変更し、境界線の高さを変更する。この境界線の高さ変更は、全ての境界線について実施する。 
Ｈｓ＝Ｈ＋ΔＨ  …（３）
ポスト処理部１０ｄは、シミュレーション結果を容易に把握可能とするための可視化処理やシミュレーション結果のデータ出力を行う。シミュレーション結果の可視化処理においては、図９に示すように、修正前の境界線Ｌと修正後の境界線の高さＨｓ（＝Ｈ＋ΔＨ）とを同時にライン表示させ、更には、修正前後の境界線を境に、液外から液内に属性を変更されたメッシュと、修正前の境界線を形成する液内メッシュとを色相や濃淡等によって識別表示する。
【００３４】
また、エアポケットの高さや境界線の内側のエアポケットの境界面積を求めることにより、実際に発生するエアポケットの評価を可能とする指標を提供する。ここで、エアポケット高さは、修正後の境界線からエアポケットの最高点までの高さであり、エアポケット境界面積は、修正した境界線で形成される平面に接するエアポケットの面積である。
【００３５】
エアポケットを評価する指標としては、エアポケット高さ、エアポケットの高さと境界面積との比、エアポケット形状の座標値（エアポケットの位置及び方向；浸漬液が浸入し易いか否か）等を指標として導入し、これらに応じて序列付けや数値化を行う。そして、シミュレーション結果を表示する際にも、エアポケットに相当する気体の属性のメッシュを、これらの序列あるいは数値に応じた色相や濃淡で区別して表示することにより、実際の浸漬工程で生じるであろうエアポケットに対する対策の難易度を知ることができる。
【００３６】
例えば、エアポケットの高さを用いてシミュレーション結果を評価する場合には、エポケット高さの値が大きい程、つまりエアポケットの深さが深い程、エアポケットを強調表示することで、エアポケットの対策が容易か否かを知ることができる。更に、詳細にエアポケットのシミュレーション結果を評価する場合には、エアポケットの高さと境界面積との比を用い、この値が大きくなる程、つまりエアポケットが細く尖った形状になればなる程、エアポケットを強調表示することで、対策の難易度を知ることができる。
【００３７】
次に、演算装置１０のＣＰＵで実行されるエアポケットのシミュレーションプログラムについて、図１０のフローチャートを用いて説明する。
【００３８】
図１０のシミュレーションプログラムでは、先ず、ステップＳ１において、対象物の表面形状のデータを複数の要素に分割して数値計算のための解析モデルを構築する。本形態においては、三角形のメッシュを用いた解析モデルを構築する。
【００３９】
次に、ステップＳ２へ進んで解析モデルにおける重力方向を設定した後、ステップＳ３へ進み、対象物が電着液等に浸漬される前の周辺が空気で満たされた状態のシミュレーションとして、一旦、全てのメッシュを気体の属性（液外メッシュ）に設定する。さらに、対象物を浸漬した初期状態のシミュレーションとして、メッシュを液体の属性と気体の属性とに分離するための初期条件を設定する。
【００４０】
本実施の形態においては、初期状態でのメッシュの属性を液体（液内メッシュ）に決定する条件として、ワークのフリーエッジに接するメッシュを初期メッシュに設定する。そして、この初期メッシュの属性を気体から液体に変更する。ここで、フリーエッジとは、ワークの縁の部分、即ちワークの端部及びワークに形成された穴の縁の部分を意味する。初期メッシュとは、ワークのフリーエッジに接する要素であって、フリーエッジに接するメッシュとは、そのメッシュを構成する辺がフリーエッジに接していることを言う。
【００４１】
尚、初期に属性を液体として決定する条件として、ワークに対してその周囲よりかつ水平より上側の視線方向から見て認識可能なメッシュを初期メッシュとしても良い。この場合においても、初期メッシュ決定後の処理は同じである。
【００４２】
次に、ステップＳ３からステップＳ４へ進み、エアポケット判定計算処理を実施して、エアポケットの発生をシミュレーションする。この判定計算処理は、各メッシュの重力と反対方向の高さの比較により、それぞれの属性を決定してエアポケットの発生を判定する処理であり、従来の手法を用いることができる。各メッシュの高さの比較は、例えば、それぞれのメッシュの重心点、内心点、外心点、垂心点等のＺ方向（重力と反対方向）における位置や、それぞれのメッシュが有する各頂点のうちの最も高い頂点（Ｚ方向の位置が最も高い頂点）等の特徴点を比較することで行う。
【００４３】
例えば、重心点の位置を比較する場合には、判定対象とするメッシュと隣接メッシュとのそれぞれの重心点のＺ方向における高さを比較し、隣接メッシュの重心点が判定対象メッシュの重心点よりも高いまたは同じ場合、隣接メッシュの属性が液体であれば、判定対象メッシュの属性を液体に決定する。一方、隣接メッシュの重心点が判定対象メッシュの重心点よりも低い場合には、判定対象メッシュの属性は変更せず、比較対象を変更して判定を繰り返す。この判定を初期メッシュを含む全メッシュについて実行し、液内メッシュと液外メッシュとの境界を気液境界線として検出する。
【００４４】
以上のステップＳ４のエアポケット判定計算処理により気液境界線が検出されると、ステップＳ５へ進む。ステップＳ５は、本発明の境界線修正ステップに相当し、ステップＳ４のエアポケット判定計算処理で得られる境界線に対して、前述の（１）式によって境界線の高さＨを求めた後、境界線の高さＨに前述の（２）式による修正高さΔＨを加算し、新たな境界線の高さＨｓを求めることにより、境界線を修正する。
【００４５】
尚、境界線の修正処理は、ステップＳ４で一度計算した結果に対して、修正高さΔＨだけ境界線を高くすることが重要である。判定計算の途中にて評価指標をΔＨだけ増分して判定することはＮＧである。
【００４６】
その後、ステップＳ６へ進み、本発明のポスト処理に相当する処理を実行する。すなわち、ステップＳ６のポスト処理では、シミュレーション結果を可視化処理して表示装置１２へ出力したり、シミュレーション結果のデータのファイルへの出力を行う。このポスト処理は、エアポケットの有無を表示（図９参照）するに止まらず、エアポケットの高さＨａｐ、或いはエアポケットの高さと境界面積との比［Ｈａｐ／Ｒａｐ］に応じて序列付けや数値化を行い、これらの序列あるいは数値に応じてエアポケットの色相を変化させる等して識別表示する。
【００４７】
例えば、図１１に示すような断面形状を有する部材を例に取ると、この部材のＡ部，Ｂ部，Ｃ部に予想されるエアポケットをシミュレーションした結果を、各部のエアポケットの高さＡＨｉ（ｉ＝１，２，３）と境界面積ＡＳｉ（ｉ＝１，２，３）との比［ＡＨｉ／ＡＳｉ］にて評価し、図１２に示すように表示する。図１２においては、［ＡＨ1／ＡＳ1］の数値が最も小さいＡ部のエアポケットを例えば青色で表示し、次に、［ＡＨ3／ＡＳ3］の数値が大きいＣ部のエアポケットを例えば黄色で表示し、［ＡＨ2／ＡＳ2］の数値が最も大きいＢ部のエアポケットを例えば赤色で表示する。
【００４８】
これにより、エアポケットの配色からエアポケットの対策の難易度を即座に判断することができる。つまり、青色のエアポケットは、部材の釣り角度等を調整することで対策できる可能性があるが、赤色のエアポケットは対策が困難であり、空気穴を設ける等の構造変更が必要であると言ったように、エアポケットの危険性をインタラクティブに判断することが可能となる。
【００４９】
以上のように、本実施の形態によれば、実際の塗装ライン等の浸漬工程の状態を考慮してエアポケットの気液境界線を修正するため、実際の浸漬工程に近い状態でのエアポケットの評価を行うことができる。また、エアポケットの高さや高さ／境界面積等によって序列或いは数値化することで、エアポケットの危険度を知ることができ、単なる目安に止まることなく、より実用的なシミュレーションへと発展させることができる。
【図面の簡単な説明】
【００５０】
【図１】シミュレーション装置の基本構成図
【図２】車体の塗装ラインの概略説明図
【図３】フロアパネルの概略斜視図
【図４】フロアパネルを複数の要素に分割した解析モデルを示す説明図
【図５】メッシュによる境界線の高さを示す説明図
【図６】ノードによる境界線の高さを示す説明図
【図７】ワークの揺れを示す説明図
【図８】境界線の変更を示す説明図
【図９】境界線の表示とメッシュ分割を示す説明図
【図１０】エアポケットシミュレーションプログラムのフローチャート
【図１１】エアポケット断面を示す説明図
【図１２】エアポケット高さ／境界面積の表示例を示す説明図
【図１３】エアポケットが発生した被塗装物のカットモデルを示す説明図
【図１４】従来の幾何学的手法によるシミュレーション結果を示す説明図
【符号の説明】
【００５１】
１      シミュレーション装置
１０    演算装置
１０ａ  モデル構築部
１０ｂ  エアポケット判定部
１０ｃ  境界線修正部
１０ｄ  ポスト処理部
【図１】
2008049269.tif 000003
【図２】
2008049269.tif 000004
【図３】
2008049269.tif 000005
【図４】
2008049269.tif 000006
【図５】
2008049269.tif 000007
【図６】
2008049269.tif 000008
【図７】
2008049269.tif 000009
【図８】
2008049269.tif 000010
【図９】
2008049269.tif 000011
【図１０】
2008049269.tif 000012
【図１１】
2008049269.tif 000013
【図１３】
2008049269.tif 000014
【図１２】
2008049269.tif 000015
【図１４】
2008049269.tif 000016