(57)【要約】  【課題】被塗装ワークを浸漬槽に投入する際のエアーポケット発生の有無を簡単な演算で求めることができるようにする。  【解決手段】被塗装ワーク１０２を浸漬槽に投入する際に、被塗装ワーク１０２に発生するエアーポケットを予測するものにおいて、計算領域Ａ内に被塗装ワーク１０２に対応する形状データを入力し、計算領域Ａ内に発生させる複数の気泡粒子Ｂを移動させると共に、この気泡粒子Ｂと形状データとの干渉からエアーポケットの発生を予測する。【選択図】図３
2006116385.tif 000002
【特許請求の範囲】
【請求項１】
被塗装ワークを浸漬槽に投入する際に、該被塗装ワークに発生するエアーポケットを予測する浸漬塗装解析方法において、
計算領域内に上記被塗装ワークに対応する形状データを入力し、
上記計算領域内に発生させる少なくとも１つの気泡粒子を移動させると共に、該気泡粒子と上記形状データとの干渉から上記エアーポケットの発生を予測する
ことを特徴とする浸漬塗装解析方法。
【請求項２】
上記気泡粒子の発生位置は任意に設定することができる
ことを特徴とする請求項１記載の浸漬塗装解析方法。
【請求項３】
上記気泡粒子の発生位置は上記形状データの線上に設定する
ことを特徴とする請求項１或いは２記載の浸漬塗装解析方法。
【請求項４】
上記気泡粒子の属性は、該気泡粒子の発生位置及び上記形状データに応じて変更する
ことを特徴とする請求項１〜３の何れかに記載の浸漬塗装解析方法。
【請求項５】
上記形状データに周囲が閉ざされた内部空間領域が存在する場合、該内部空間領域を固体領域と設定して、該固体領域では上記気泡粒子を発生させない
ことを特徴とする請求項１〜４の何れかに記載の浸漬塗装解析方法。
【請求項６】
上記形状データを表す上記線に板厚データが設定されていない場合、上記気泡粒子を該線の外表面と内表面との双方から発生させる
ことを特徴とする請求項３記載の浸漬塗装解析方法。
【請求項７】
上記線で囲まれた閉空間領域が存在する場合、該閉空間領域を固体領域と設定して、該固体領域側には上記気泡粒子を発生させない
ことを特徴とする請求項６記載の浸漬塗装解析方法。
【請求項８】
上記気泡粒子の移動を計算するに際し、該気泡粒子の前回の位置座標と今回の位置座標とを比較し、両位置座標が同一の場合は、当該気泡粒子を停滞粒子として設定し、該停滞粒子に対するその後の移動計算を打ち切る
ことを特徴とする請求項１〜７の何れかに記載の浸漬塗装解析方法。
【請求項９】
移動計算が打ち切られた上記気泡粒子の属性を変更する
ことを特徴とする請求項８記載の浸漬塗装解析方法。
【請求項１０】
移動計算が打ち切られた上記気泡粒子に接触する上記形状データの線部分と該気泡粒子との少なくとも一方にマーキングを付す
ことを特徴とする請求項８或いは９記載の浸漬塗装解析方法。
【請求項１１】
移動する上記気泡粒子の軌跡を表示装置上に表示させる
ことを特徴とする請求項１〜１０の何れかに記載の浸漬塗装解析方法。
【請求項１２】
上記計算領域に、上記形状データの少なくとも一部を取り囲む外部領域を設定し、
上記気泡粒子が上記外部領域から外れている場合、該気泡粒子の移動計算は行わない
ことを特徴とする請求項１〜１１の何れかに記載の浸漬塗装解析方法。
【請求項１３】
全ての上記気泡粒子の移動計算が打ち切れられた場合、上記各気泡粒子の移動計算を全て終了する
ことを特徴とする請求項８〜１０の何れかに記載の浸漬塗装解析方法。
【請求項１４】
１つの座標領域に集まる上記停滞粒子の数をカウントし、該カウント値に基づき上記エアーポケットの発生し易い部位を特定する
ことを特徴とする請求項１〜１３の何れかに記載の浸漬塗装解析方法。
【請求項１５】
上記カウント値に応じて上記形状データの上記停滞粒子の干渉している線部分と該停滞粒子との少なくとも一方に付す上記マーキングの色を変化させる
ことを特徴とする請求項１４記載の浸漬塗装解析方法。
【請求項１６】
上記気泡粒子の移動方向は計算途中において任意に変更可能である
ことを特徴とする請求項１〜１５の何れかに記載の浸漬塗装解析方法。
【請求項１７】
上記気泡粒子の直径を設定し、該直径と上記形状データに形成される通気孔との大きさを比較して、エアーポケットの発生位置を検出する
ことを特徴とする請求項１〜１６の何れかに記載の浸漬塗装解析方法。
【請求項１８】
上記気泡粒子の移動は重力と反対方向に発生させる
ことを特徴とする請求項１〜１７の何れかに記載の浸漬塗装解析方法。
【請求項１９】
上記気泡粒子の移動計算は上記計算領域に設定した塗料の流速に基づいて設定する
ことを特徴とする請求項１〜１７の何れかに記載の浸漬塗装解析方法。
【請求項２０】
上記重力方向を上記計算領域の下方に設定する
ことを特徴とする請求項１８記載の浸漬塗装解析方法。
【請求項２１】
上記気泡粒子を該気泡粒子に設定する密度と上記計算領域に設定する密度との比に基づいて移動させることを特徴とする請求項１〜１７記載の浸漬塗装解析方法。
【請求項２２】
上記気泡粒子に設定する密度と上記計算領域に設定する密度とを逆に設定する
ことを特徴とする請求項２１記載の浸漬塗装解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、車体部品等のワークを浸漬槽に浸漬して塗装処理を施す浸漬塗装解析方法に関する。
【背景技術】
【０００２】
従来から、浸漬塗装として、ワークを電着槽に浸漬し、電着塗料とワークとの間に電流を流し、電気泳動作用により車体部品の表面に電着塗膜を析出させた後、乾燥させる電着塗装が知られている。電着塗装は、複数枚の板金部材を重ねた状態でスポット溶接することで形成されたワークに対しても施すことができるため、自動車の車体部品（ホワイトボディを含む）等に、広く採用されている。
【０００３】
又、最近の自動車の車体構造は、軽量化、及び耐衝撃性の向上から、パネル本体の裏面に溶接される補強部材の立体形状が複雑化する傾向にある。このような形状の複雑な補強部材が溶接されているワークに対して浸漬塗装を施した場合、塗料がパネル本体の裏面と補強部材との間等に形成される袋状部に移動込み難くなり、その内面には塗装が全く或いは充分に付着されず、塗装不良が発生し易くなる。この場合の塗装不良の代表としてエアーポケットがある。エアーポケットは、ワークを電着槽に浸漬する際に、補強部材とパネル本体との間等の内部に閉じこめられた空気が外部へ抜けきらずに残留し、空気溜り（エアーポケット）が発生するものである。
【０００４】
一般に、複雑な立体形状のワークに対してエアーポケットを発生させることなく、良好な塗装を行うためには、先ず、そのワークの試作品を作成し、様々な投入姿勢及び投入方向から電着槽に浸漬させて、その都度、エアーポケットの発生具合を調べ、エアーポケットの発生しないワークの投入姿勢及び投入方向を決定するようにしている。
【０００５】
しかし、エアーポケットの発生具合を、人為的な手法により一々調べることは、多大な手間と時間、及び費用を必要とするため非効率的である。又、このような手法では、浸漬塗装すべきワークを設計し、試作品を完成させてから浸漬塗装に適したワークの投入姿勢及び投入方向を決定することとなるため、製品の設計から実際の量産へ移行するまでの期間が長期化し、製品開発に要する時間を短縮化する上での障害となっている。
【０００６】
そのため、最近では、例えば特許文献１（特開２０００−５１７５０号公報）に開示されているように、設計段階において試作品を作成することなく、コンピュータによる解析プログラムを用いて、浸漬処理すべき製品についての、エアーポケット発生状況をディスプレイ上で予測するシステムが採用されている。
【０００７】
すなわち、同文献に開示されている技術は、ワークの各表面をその頂点座標データに基づいて任意数の解析用メッシュで形成し、各メッシュについて、ワークを浸漬槽に投入する姿勢を設定し、この姿勢でワークを浸漬槽に模擬浸漬した際にエアーポケットとなりうる表面部位か否かを判定し、その判定の結果をディスプレイ上に表示するものである。
【特許文献１】特開２０００−５１７５０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかし、上述した特許文献１に開示されている技術では、複数に分割されたメッシュ毎に、空気と浸漬槽内の塗料との状態変化を時々刻々と計算しなければならないため、演算が複雑で、その分コンピュータにかかる負荷が大きくなり、その結果、演算に膨大な時間を必要とし、パーソナルコンピュータ等の汎用コンピュータを用いて演算を行うには限界があり、実用性に欠ける問題がある。
【０００９】
本発明は、上記事情に鑑み、ワークを浸漬槽に投入する際のエアーポケット発生の有無を簡単な演算で求めることができ、コンピュータにかかる負荷を小さくして計算時間の短縮化を実現することができると共に、パーソナルコンピュータ等の汎用コンピュータを用いた演算を可能として、実用性に優れた浸漬塗装解析方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記目的を達成するため本発明は、被塗装ワークを浸漬槽に投入する際に、該被塗装ワークに発生するエアーポケットを予測する浸漬塗装解析方法において、計算領域内に上記被塗装ワークに対応する形状データを入力し、上記計算領域内に発生させる少なくとも１つの気泡粒子を移動させると共に、該気泡粒子と上記形状データとの干渉から上記エアーポケットの発生を予測することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、気泡粒子の移動からエアーポケットの発生を予測するようにしたので、計算が容易になり、その分、コンピュータにかかる負荷を小さくして、相対的に計算時間の短縮化を実現することができる。更に、コンピュータにかかる負荷が小さいため、汎用コンピュータを用いて演算を行うことが可能となり実用的である。
【発明を実施するための最良の形態】
【００１２】
以下、図面に基づいて本発明の一形態を説明する。図１〜図８に本発明の第１形態を示す。図１は電着塗装解析装置の機能ブロック図である。尚、以下においては、浸漬塗装として電着塗装を用いた場合を例に掲げて説明する。
【００１３】
電着塗装解析装置１は、ＲＯＭ，ＲＡＭ等のメモリ及びＣＰＵを有するマイクロコンピュータ等のコンピュータで構成されており、後述する各処理部１１，１２，１４の機能を実現するためのプログラムをメモリにロードして実行することにより浸漬塗装の解析処理が実行される。
【００１４】
図１に示すように、電着塗装解析装置１は、入力側にキーボード、マウス、外部接続端子等を含む入力装置２が接続され、出力側に液晶モニタ等の表示装置３が接続されている。又、電着塗装解析装置１は、解析モデル設定処理部１１、モデル解析処理部１２、解析モデルデータベース１３、及び表示処理部１４を備えている。
【００１５】
解析モデル設定処理部１１には解析に必要なデータが格納されている。解析に必要なデータとしては、仮想電着槽として規定した三次元の計算領域Ａ（図４参照）内の物性値（流体密度、粘性係数等）、計算領域Ａ内に発生させる空気等の気泡粒子Ｂの物性値（密度、粘性係数等）、重力（Ｇ）の方向、及び計算領域Ａ内に定義される被塗装ワーク１０２の形状データが設定されている。
【００１６】
尚、本形態では、仮想電着槽内に電着塗料が装填されていることを想定しており、従って、計算領域Ａ内の物性値は、電着塗料の有する物性値となる。又、気泡粒子Ｂとして本形態は空気を想定しており、その物性値は空気の有する物性値となる。但し、気泡粒子Ｂの移動計算中に、各気泡粒子Ｂが分裂したり、合体したりすることは考慮しない。
【００１７】
これら各データは、入力装置２から入力される。又、被塗装ワーク１０２としては、例えば自動車等の車両においては、ドアパネル、エンジンフードパネル、ホワイトボディー等、電着塗装を必要とする車体の各部品がある。この被塗装ワーク１０２の形状データは、設計者自身が直接作成し、或いは予め作成されているＣＡＤデータ、ＳＴＬ(Stereo Lithography)データ等に基づいて自動的に作成されるものであり、入力装置２を介して直接入力しても良いが、外部接続端子に、図示しないＬＡＮ等の通信手段を介して接続された、他のコンピュータに備えられている記憶装置（図示せず）から読込むようにしても良い。
【００１８】
モデル解析処理部１２は、解析モデル設定処理部１１で設定される計算領域Ａの物性値、気泡粒子Ｂの物性値、被塗装ワーク１０２との対応付けが行われた形状モデル、及び入力される計算条件等の解析条件に従って、被塗装ワーク１０２の形状モデル内部の気泡粒子Ｂの移動計算を行う。尚、モデル解析処理部１２で実行される気泡粒子Ｂの移動計算処理については後述する。
【００１９】
解析モデルデータベース１３は、モデル解析処理を行った結果得られた解析モデルを記録するデータベースであり、ハードディスク装置や光磁気ディスク装置、フラッシュメモリ等の書き込み可能な不揮発性記録媒体で構成されている。解析モデル設定処理部１１は、解析モデルデータベース１３から未設定の解析モデルを読込み、この解析モデルの構成要素に対して、計算領域Ａ、気泡粒子Ｂ、被塗装ワーク１０２の構成要素との対応付けを行い、解析モデルデータベース１３へ再度記録する。モデル解析処理部１２は、解析モデルデータベース１３に記録された、解析モデル設定処理部１１で設定済みの解析モデルを取得して気泡粒子Ｂの移動計算を行う。
【００２０】
表示処理部１４は、気泡粒子Ｂの移動計算処理を実行している際の制御画像や作業画像等の表示画像を表示装置３へ出力する。尚、制御画像には、解析モデル設定処理部１１で実行可能な設定機能のメニュー表示部と、モデル解析処理部１２で実行可能な解析機能のメニュー表示部とが含まれる。作業画像には、各メニュー表示部において選択された内容の条件設定画像を作成して表示する。表示処理部１４は、モデル解析処理部１２で実行した移動計算により取得した計算領域Ａ、及び被塗装ワーク１０２の形状モデル、気泡粒子Ｂの状態変化を表示装置３に表示させる。
【００２１】
尚、解析モデル設定処理部１１、モデル解析処理部１２、表示処理部１４は、コンピュータによって処理される必要はなく、それぞれ専用のハードウェアにより実現するようにしても良い。
【００２２】
上述したモデル解析処理部１２で実行される気泡粒子Ｂの移動計算は、具体的には、図２、図３に示す気泡粒子移動計算ルーチンに従って処理される。
【００２３】
このルーチンでは、先ず、ステップＳ１で、解析モデル設定処理部１１で対応付けが行われた被塗装ワーク１０２の形状データを入力する。すると、図４（ａ）に示すように、表示装置３に表示されている計算領域Ａに、被塗装ワーク１０２の形状が表示される。尚、実際の被塗装ワーク１０２は、表示装置３上に三次元表示されるが、本形態では便宜的に二次元で表示としている。
【００２４】
同図（ａ）に例示されている被塗装ワーク１０２は、その外周一側に開口部１０２ｂが形成され、内部に周囲が閉ざされた内部空間領域１０２ａが形成された角形部材が示されている。尚、内部空間領域１０２ａには中実状の内部領域も含まれる。
【００２５】
又、表示装置３に表示された被塗装ワーク１０２は、平面上で傾斜方向（三次元表示の場合は、これに奥行きを加えた方向への傾斜）を自在に設定することができる。本形態では、開口部１０２ｂを斜め下方へ指向させた状態が示されている。
【００２６】
そして、ステップＳ２で、計算領域Ａに充満されている流体の物性値（流体密度、粘性係数等）、及び気泡粒子Ｂの物性値（密度、粘性係数等）をそれぞれ入力する。
【００２７】
次いで、ステップＳ３で重力Ｇの方向を設定する。尚、図５に示すように、本形態では、重力Ｇの方向を縦軸に沿った下方に設定している。
【００２８】
次に、ステップＳ４へ進み、計算条件を入力する。計算条件は気泡粒子Ｂの発生個数Ｎ、規定回数Ｍ、気泡粒子Ｂの演算毎の移動距離ΔＺである。尚、発生個数Ｎ、規定回数Ｍの最小単位は０であり、Ｎ＝０、或いはＭ＝０の場合は、１個の気泡粒子Ｂを１回のみ計算し、移動距離ΔＺだけ上昇させることになる。
【００２９】
その後、ステップＳ５へ進み、Ｎ個の気泡粒子Ｂの発生位置を、例えば図４（ａ），（ｂ）の黒丸で示すように設定する。この気泡粒子Ｂの発生位置は、表示装置３を参照しながら、マウス等を用いて任意の位置を指定することで設定する。この場合、表示装置３にグリッドを表示し、このグリッドをマウスで指定することで、気泡粒子Ｂの発生位置を設定するようにしても良い。或いは開始行を指定して一括設定するようにしても良い。又、この発生位置は、被塗装ワーク１０２の内部に設定することも可能である。被塗装ワーク１０２の内部に設定することで、最も注目した部位を集中的に解析することができる。
【００３０】
次いで、ステップＳ６へ進み、解析が開始されるまで待機する。そして、[Enter]キーをＯＮする等して、解析開始の指令が入力されると、ステップＳ７へ進み、ステップＳ７以下で、気泡粒子Ｂの移動計算が開始される。先ず、ステップＳ７では、発生位置に設定されている気泡粒子Ｂを、ステップＳ４で設定した移動距離ΔＺ分、重力（Ｇ）と反対の方向、すなわち図４では上方へ移動させる。
【００３１】
その結果、気泡粒子Ｂは演算周期毎に移動距離ΔＺずつ上昇する。尚、図４（ａ），（ｂ）には、便宜的にＮ個の気泡粒子Ｂを横軸方向へ１列に配設した状態が示されているが、気泡粒子Ｂの発生位置は１列に設定する必要はない。又、図４（ａ），（ｂ）においては、気泡粒子Ｂの演算周期毎に移動する軌跡が一点鎖線で示され、表示装置３上に実際に表示される気泡粒子Ｂが黒丸で示されている。従って、本形態に示すフローチャートに従って、計算を繰り返すと、各列毎に１つの気泡粒子Ｂ（黒丸で示す）のみが移動（浮上）するように表示される。
【００３２】
次いで、ステップＳ８へ進み、計算領域Ａ内に配設した被塗装ワーク１０２の形状を表す線と気泡粒子Ｂとが干渉するか否かが調べられる。これは、同一座標上に、被塗装ワーク１０２の形状を表す線と気泡粒子Ｂとが存在するか否で判断し、同一座標上に無い場合は、ステップＳ１１へジャンプする。
【００３３】
一方、被塗装ワーク１０２の形状を表す線と気泡粒子Ｂとが同一座標上にある場合は、干渉しているため、ステップＳ９へ進む。ステップＳ９では、干渉している部位の表面１０２ｄ（図５参照）の被塗装ワーク１０２の気泡粒子Ｂの移動方向を基準とする傾斜角θを調べ、−９０°＜θ＜＋９０°のときは、ステップＳ１０へ進む。又、それ以外のとき（−９０°≧θ、或いはθ≧＋９０°）はステップＳ１１へジャンプする。
【００３４】
図５に示すように、本形態では気泡粒子Ｂの移動方向を基準として反時計回り方向の角度をプラス（＋）、時計回り方向の角度をマイナス（−）として設定されており、＋９０°と−９０°とは水平方向を示す。従って、被塗装ワーク１０２の表面１０２ｄの傾斜角θが、−９０°≧θ、或いはθ≧＋９０°の場合は、気泡粒子Ｂが被塗装ワーク１０２の表面１０２ｄに阻害されてそれ以上、上昇することはできない。一方、被塗装ワーク１０２の表面１０２ｄの傾斜角θが、−９０°＜θ＜＋９０°の範囲にあるときは、その表面１０２ｄに沿って気泡粒子Ｂを上昇させることが可能である。
【００３５】
従って、被塗装ワーク１０２の表面１０２ｄの傾斜角θが、−９０°＜θ＜＋９０°の範囲にあるときは、ステップＳ１０へ進み、気泡粒子Ｂを被塗装ワーク１０２を表す線に沿って移動（上昇）させる。
【００３６】
そして、ステップＳ８，Ｓ９或いはＳ１０からステップＳ１１へ進むと、気泡粒子Ｂの数（以下「気泡粒子数」と称する）ｎが、初期設定した発生個数Ｎに達したか否かを調べ、達していない場合（ｎ＜Ｎ）、ステップＳ１２へ分岐し、気泡粒子数ｎをインクリメントして（ｎ←ｎ＋１）、ステップＳ７へ戻る。
【００３７】
そして、ステップＳ１１において、気泡粒子数ｎが発生個数Ｎに達したと判定された場合（ｎ＝Ｎ）、気泡粒子Ｂの全てが移動距離ΔＺだけ上昇されたと判断し、ステップＳ１３へ進み、計算回数ｍが規定回数Ｍに達したか否かを調べる。そして、ｍ＜Ｍのときは、規定回数に達していないため、ステップＳ１４へ進み、計算回数ｍをインクリメントして（ｍ←ｍ＋１）、ステップＳ１５へ進み、気泡粒子数ｎをクリアして（ｎ←０）、ステップＳ７へ進み、気泡粒子Ｂを移動距離ΔＺだけ移動（上昇）させる。
【００３８】
そして、ステップＳ１３において、計算回数ｍが規定回数Ｍに達したとき判断したとき、ステップＳ１６へ進み、各列の気泡粒子Ｂの軌跡をポスト処理により、アニメーション化して可視化表示した後、ルーチンを終了する。
【００３９】
その結果、表示装置３には、図４（ａ）に示すように、各列に設定した気泡粒子Ｂが、気泡粒子移動計算ルーチンで、解析した軌跡に沿って発生位置から規定回数Ｍ達する位置まで上昇する様子が繰り返し表示される。
【００４０】
図４（ａ）に示す被塗装ワーク１０２では、その下方から発生した気泡粒子Ｂが上昇する過程で、被塗装ワーク１０２の表面１０２ｄに干渉した気泡粒子Ｂは、この表面１０２ｄに沿って上昇する。又、この被塗装ワーク１０２に開口されている開口部１０２ｂから内部に浸入した気泡粒子Ｂの内、車体用ワークの内表面１０２ｄに干渉した気泡粒子Ｂはその内表面１０２ｄに沿って上昇される。しかし、図４（ａ）に示す被塗装ワーク１０２内の上部は閉塞されているため、この被塗装ワーク１０２内に浸入した気泡粒子Ｂは、内部に停滞する。その結果、図４（ａ）に表示された被塗装ワーク１０２では、気泡粒子Ｂが停滞する部位にエアーポケットが発生することが視覚的に理解できる。
【００４１】
そして、今度は、図４（ｂ）に示すように、気泡粒子Ｂが停滞する部位、すなわちエアーポケットが発生する部位に通気孔１０２ｃを形成し、この通気孔１０２ｃが形成された被塗装ワーク１０２を用いて、図２、図３に示す気泡粒子移動計算ルーチンを再度実行する。すると、被塗装ワーク１０２内に停滞していた気泡粒子Ｂが通気孔１０２ｃから外部に抜ける様子が視覚的に認識されるので、エアーポケット対策を容易に検討することができる。
【００４２】
このように、本形態では、解析用メッシュを用いず、気泡粒子Ｂの移動計算から、浸漬槽に投入する被塗装ワーク１０２の内部にエアーポケットが発生するか否かの解析を行うようにしたので、計算時間を大幅に短縮することができるばかりでなく、対策すべき通気孔の形成位置を瞬時に検討することができるため作業性が良い。
【００４３】
又、コンピュータにかかる負荷が小さくなり、その結果、パーソナルコンピュータ等の汎用コンピュータを用いた計算が可能となり実用的である。
【００４４】
ところで、図２に示す気泡粒子移動計算ルーチンのステップＳ４において、計算条件として入力する項目に、気泡粒子Ｂの直径Ｄを追加することも可能である。計算条件に気泡粒子Ｂの直径Ｄを加えることで、図６に示すように、被塗装ワーク１０２に形成する通気孔１０２ｃの大きさＬをどの程度にすべきかを、気泡粒子Ｂとの関係で検討することができる。すなわち、気泡粒子Ｂの直径Ｄよりも通気孔１０２ｃの大きさＬが小さい場合、気泡粒子Ｂは、通気孔１０２ｃから抜けることができないため、内部にエアーポケットが発生する。従って、気泡粒子Ｂの直径Ｄを計算条件として入力することで、通気孔１０２ｃの幅Ｌ（或いは径）をどの程度にすべきかを容易に検討することができる。
【００４５】
又、図２に示すフローチャートのステップＳ３において、重力（Ｇ）方向を上方に設定すれば、図７、図８に示すように、気泡粒子Ｂは、計算領域Ａ内で下方へ移動させることができる。気泡粒子Ｂの下方への移動を検討することで、例えば被塗装ワーク１０２を電着槽から引き上げたときの液溜まり対策を講じることができる。
【００４６】
すなわち、図７（ａ）に示すように、図４（ａ）と同じ姿勢の状態で被塗装ワーク１０２を電着槽から引き上げた場合は、下部にエアーポケットが形成されるが、同図（ｂ）に示すように、通気孔１０２ｃが下方へ指向するように回転させることで、下部から電着塗料を排出させることができる。
【００４７】
尚、この場合、重力（Ｇ）方向は下方に設定した状態とし、ステップＳ４に示す計算条件として、計算領域Ａ内の密度ρｏと気泡粒子Ｂの密度ρとの比を、ρｏ＜ρとなるような値に設定しても同様の結果を得ることができる。
【００４８】
又、図９、図１０に本発明の第２形態を示す。本形態は第１形態の変形例である。第１形態では、気泡粒子移動計算ルーチン終了後に、表示装置３に対して気泡粒子Ｂの移動をアニメーション化して動的に表示するようにしたが、本形態では、計算途中においても、気泡粒子Ｂの移動移動を画像表示するようにしたものである。
【００４９】
図９に示すルーチンは、図２、図３に示す気泡粒子移動計算ルーチンのステップＳ１０とステップＳ１１との間に追加されるでもので、その他のステップは第１形態と同一の処理が行われるため、同一のステップ部分については説明を省略する。
【００５０】
すなわち、図３に示す気泡粒子移動計算ルーチンのステップＳ８，Ｓ９或いはＳ１０からステップＳ２１へ進むと、今回の気泡粒子Ｂが移動した位置座標を解析モデルデータベース１３（図１参照）に記憶させ、その後、ステップＳ２２へ進み、この気泡粒子Ｂを表示装置３に画像表示させて、ステップＳ１１へ進む。
【００５１】
その結果、図１０（ａ），（ｂ）に示すように、表示装置３には計算途中の時々刻々と変化する気泡粒子Ｂの挙動が表示される。操作者は計算途中から気泡粒子Ｂの挙動を視覚的に認識することができるため、計算結果が表示される迄のいらだちを解消することができる。
【００５２】
又、図１１、図１２に本発明の第３形態を示す。本形態は第１形態の変形例である。第１形態では、気泡粒子Ｂの開始位置を任意の位置から設定できるようにしたが、本形態では、気泡粒子Ｂの開始位置を、被塗装ワーク１０２を表す線に近接する位置に設定した場合、最も近い線上に気泡粒子Ｂの開始位置を自動的に設定するようにしたものである。
【００５３】
図１１に示すルーチンは、図２、図３に示す気泡粒子移動計算ルーチンのステップＳ５とステップＳ６との間に追加されるもので、その他のステップは第１形態と同一の処理が行われるため、同一のステップ部分については説明を省略する。
【００５４】
ステップＳ５で、Ｎ個の気泡粒子Ｂの発生位置を設定すると、ステップＳ３１へ進み、設定した各気泡粒子Ｂ中に、被塗装ワーク１０２を表示する線に近いものがあるか否かを調べ、この線に近い位置に設定した気泡粒子Ｂがある場合は、当該気泡粒子Ｂを被塗装ワーク１０２を表す線の中で、最も近い位置にある線上に貼り付ける。従って、操作者がステップＳ５において、被塗装ワーク１０２を表す線に近接する位置に意識的に気泡粒子Ｂを設定することで、図１２（ａ）に示すように、当該気泡粒子Ｂは最も位階位置の線上に自動的に貼り付けられる。尚、同図に示す気泡粒子Ｂ’のように、線から比較的離れた位置に設定した場合は、その位置に気泡粒子Ｂ’の開始位置が設定される。
【００５５】
又、通常、被塗装ワーク１０２を表す線には板厚データが設定されていないため、当該線上に貼り付けられた気泡粒子Ｂが、被塗装ワーク１０２の外表面に張り付いたのか、内表面に張り付いたのか判別できない。そのため、図３のステップＳ７において、当該線に貼り付けられている気泡粒子Ｂを移動距離ΔＺだけ移動させる際しては、図１２（ｂ）に示すように、当該気泡粒子Ｂを外表面と内表面との双方から発生させる。
【００５６】
本形態では、気泡粒子Ｂの開始位置を、被塗装ワーク１０２の線上から発生させるようにしたので、計算時間を短縮することができる。すなわち、気泡粒子Ｂが停滞する部位は被塗装ワーク１０２を表す線と干渉する位置に形成されるため、気泡粒子Ｂの発生位置を線上に設定することで、計算をより効率的に行うことができる。又、その際、当該線に板厚が設定されていない場合は、気泡粒子Ｂを外表面と内表面との双方から発生させるようにしたので、より精度の高い解析を行うことができる。
【００５７】
又、図１３に本発明の第４形態を示す。本形態は第１形態の変形例である。第１形態では、気泡粒子Ｂの開始位置を任意の位置から設定できるようしたが、本形態では、被塗装ワーク１０２の閉空間領域に気泡粒子Ｂの開始位置を設定した場合、その設定を自動的に取り消し、気泡粒子Ｂを閉空間領域に発生させないようにしたものである。
【００５８】
図１３に示すルーチンは、図２、図３に示す気泡粒子移動計算ルーチンのステップＳ５とステップＳ６との間に追加されるもので、その他のステップは第１形態と同一の処理が行われるため、同一のステップ部分については説明を省略する。
【００５９】
すなわち、本形態では、気泡粒子移動計算ルーチンを実行するに際し、ステップＳ１（図２参照）で被塗装ワーク１０２の形状データを入力したとき、この被塗装ワーク１０２の形状データに周囲が閉ざされた閉空間領域あるか否かを調べ、閉空間領域が存在する場合は、これを固体領域として設定する。尚、本形態による被塗装ワーク１０２では内部空間領域１０２ａがこれに相当する。
【００６０】
そして、図３に示す気泡粒子移動計算ルーチンのステップＳ５でＮ個の気泡粒子Ｂの発生位置を設定するした後、ステップＳ４１へ進むと、発生位置に設定した各気泡粒子Ｂの中に、固体領域に設定したものがあるか否かを調べ、存在していない場合は、そのままステップＳ６へ進み、存在している場合は、ステップＳ４２へ分岐し、当該気泡粒子Ｂの設定を取り消した後、ステップＳ６へ進む。
【００６１】
本形態では、気泡粒子Ｂの挙動を解析する必要のない固体領域に気泡粒子Ｂの発生位置が設定された場合、これを取り消すことで、計算効率が向上し、計算時間の短縮化を実現することができる。
【００６２】
又、図１４〜図１６に本発明の第５形態を示す。本形態は第１形態の変形例である。第１形態では、発生位置から発生させた各列の気泡粒子Ｂの一部が、例えば図４（ａ）に示すような被塗装ワーク１０２の内部に浸入した場合、気泡粒子Ｂが上部の閉塞されている空間に達したにも拘わらず、気泡粒子Ｂの計算回数ｍが規定回数Ｍに達するまでは繰り返し計算が行われる。本形態では、気泡粒子Ｂが前回の位置座標と同じ位置にある場合、すなわち、今回のルーチン実行時に移動距離ΔＺだけ移動させたにも拘わらず、移動していない場合は、当該気泡粒子Ｂの計算を打ち切るようにしたものである。
【００６３】
尚、図１４に示すルーチンは、図２、図３に示す気泡粒子移動計算ルーチンのステップＳ１１，１２に代えて適用するもので、その他のステップは第１形態と同一の処理が行われるため、同一のステップ部分については説明を省略する。
【００６４】
すなわち、図３に示す気泡粒子移動計算ルーチンのステップＳ８，Ｓ９或いはＳ１０からステップＳ５１へ進むと、今回の気泡粒子Ｂが移動した位置座標を解析モデルデータベース１３（図１参照）に記憶させ、ステップＳ５２へ進む。
【００６５】
ステップＳ５２では、今回の気泡粒子Ｂの座標位置と、当該気泡粒子Ｂの前回の座標位置とを比較し、同一か否かを調べる。そして、前回と今回との気泡粒子Ｂの座標位置が相違しているときは、そのままステップＳ１１へ進む。一方、前回と今回との気泡粒子Ｂの座標位置が同一のときは、ステップＳ５３へ分岐し、停滞粒子設定処理を実行して、ステップＳ１１へ進む。
【００６６】
ステップＳ５３で実行される停滞粒子設定処理は、図１５に示す停滞粒子設定処理サブルーチンに従って行われる。このルーチンでは、先ず、ステップＳ５６で停滞している気泡粒子Ｂが干渉している被塗装ワーク１０２の形状データの線（図１４の太線で示す部分）にマーキングを付す。マーキングは、例えば図１４に示すように、停滞している気泡粒子Ｂと干渉している部分の線の太さや色を変える等により行う。停滞している気泡粒子Ｂが干渉している線にマーキングを付すことで、エアーポケットの発生する部位を一目で認識することができる。
【００６７】
次いで、ステップＳ５７へ進み、停滞している気泡粒子Ｂの属性を変更し、この気泡粒子Ｂを停滞粒子Ｂ’として設定する。属性の変更は、停滞粒子Ｂ’の形状や色を、通常の気泡粒子Ｂと異なる形状にする等により行う。例えば、図１６では、通常の気泡粒子Ｂの形状を円形とした場合、停滞粒子Ｂの形状を三角形としている。
【００６８】
停滞粒子Ｂ’の形状を、通常の気泡粒子Ｂと異なる形状とすることで、計算中の気泡粒子Ｂと、計算が打ち切られた停滞粒子Ｂ’とを容易に区別することができる。
【００６９】
その後、ステップＳ５８へ進み、当該列ｎを停滞列ｎｓとして記憶し、ルーチンを抜け、図１４に示す気泡粒子移動計算ルーチンのステップＳ１１へ進む。
【００７０】
ステップＳ１１では、気泡粒子数ｎが発生個数Ｎに達したか否かを調べ、達していない場合（ｎ＜Ｎ）、ステップＳ５４へ分岐し、気泡粒子数ｎをインクリメントして（ｎ←ｎ＋１）、ステップＳ５５へ進む。又、達した場合（ｎ＝Ｎ）、ステップＳ１３へ進む。
【００７１】
ステップＳ５５へ進むと、現在の気泡粒子数ｎ、すなわち、今回の列ｎが、停滞粒子Ｂ’を表す列（停滞列）ｎｓか否かを調べ、ｎ≠ｎｓのときは、ステップＳ７へ戻り、気泡粒子Ｂを上昇させる。一方、ｎ＝ｎｓのときはステップＳ１１戻る。その結果、停滞列ｎｓにある停滞粒子Ｂ’は、結果的にその後の計算が打ち切られる。
【００７２】
このように、本形態では、１つの列ｎの気泡粒子Ｂの前回の位置座標と今回の位置座標とが同一のときは、気泡粒子Ｂが被塗装ワーク１０２の形状データを表す線に干渉して停滞していると判断し、当該気泡粒子Ｂのその後の計算を打ち切るようにしたので、計算効率が向上し、計算時間を短縮することができる。
【００７３】
更に、当該線の干渉部分にマーキングを付すと共に、停滞している気泡粒子Ｂの属性を変更して停滞粒子Ｂ’とするようにしたので、操作者は表示装置３に表示されたマーキング、及び停滞粒子Ｂ’を認識することで、エアーポケットの発生する部位を簡単に認識することができ、当該部位のエアーポケット対策が容易になる。
【００７４】
又、図１７〜図１９に本発明の第６形態を示す。本発明は上述した第５形態の変形例である。第５形態では、気泡粒子Ｂの座標位置が前回と今回とで同じ場合、当該気泡粒子Ｂを停滞粒子Ｂ’として設定することで、以後の計算を打ち切るようにしたが、本形態では、気泡粒子Ｂを球体として設定し、特定箇所に集まる停滞粒子Ｂ’の数をカウントすることで、エアーポケットの発生し易い部位を特定しようとするものである。
【００７５】
図１７に示すルーチンは、図１４に示す気泡粒子移動計算ルーチンのステップＳ１１とステップＳ１３との間に、ステップＳ６１，Ｓ６２を追加したもので、その他の処理は、図１４に示すフローチャートと同一である。又、図１７のステップＳ５３で実行される停滞粒子設定処理サブルーチンは、図１８に示すフローチャートに従って処理される。
【００７６】
又、本形態では、気泡粒子移動計算ルーチンを実行するに際し、ステップＳ４（図２参照）では、計算条件として気泡粒子Ｂの直径Ｄを追加入力する。尚、この直径Ｄは操作者が任意に設定することができる。
【００７７】
そして、図３に示す気泡粒子移動計算ルーチンのステップＳ８，Ｓ９或いはＳ１０からステップＳ５１へ進むと、今回の気泡粒子Ｂが移動した位置座標を解析モデルデータベース１３（図１参照）に記憶させ、ステップＳ５２へ進む。
【００７８】
ところで、ステップＳ１０では、気泡粒子Ｂが干渉している被塗装ワーク１０２の形状データの線に沿って気泡粒子Ｂを移動（上昇）させるが、この場合、図１９（ｂ），（ｃ）に示すように、本形態による気泡粒子Ｂは、直径Ｄにより領域が設定されているため、この気泡粒子Ｂの外表面が、形状データの線に沿って移動する。
【００７９】
そして、ステップＳ５１からステップＳ５２へ進むと、今回の気泡粒子Ｂの座標位置と、当該気泡粒子Ｂの前回の座標位置とを比較し、同一か否かを調べる。そして、前回と今回との気泡粒子Ｂの座標位置が相違しているときは、そのままステップＳ１１へ進む。一方、前回と今回との気泡粒子Ｂの座標位置が同一のとき、すなわち、気泡粒子Ｂが被塗装ワーク１０２の表面に干渉し、それ以上移動できない状態のときは、ステップＳ５３へ分岐し、停滞粒子設定処理を実行して、ステップＳ１１へ進む。
【００８０】
尚、図１９に示すように、本形態で採用する被塗装ワーク１０２は、内部空間の上部が閉塞されており、又、被塗装ワーク１０２の両側に、山状に屈曲されたフランジ部１０２ｆが形成されている。この場合、気泡粒子Ｂは、被塗装ワーク１０２は内部空間の上部と、フランジ部１０２ｆの下方へ指向する表面に集まり易い。
【００８１】
ステップＳ５３で実行される停滞粒子設定処理は、図１８に示す停滞粒子設定処理サブルーチンに従って行われる。このルーチンでは、先ず、ステップＳ６６で停滞している気泡粒子Ｂが干渉している被塗装ワーク１０２の形状データの線にマーキングを付す。マーキングは、例えば図１９（ｂ），（ｃ）に太線で示す干渉部分を他の線と異なる色で表示することで行う。
【００８２】
次いで、ステップＳ６７へ進み、当該列ｎを停滞列ｎｓとして記憶し、ルーチンを抜け、図１７に示すステップＳ１１へ進む。
【００８３】
ステップＳ１１では、気泡粒子数ｎが発生個数Ｎに達したか否かを調べ、達していない場合（ｎ＜Ｎ）、ステップＳ５４へ分岐し、気泡粒子数ｎをインクリメントして（ｎ←ｎ＋１）、ステップＳ５５へ進む。又、達した場合（ｎ＝Ｎ）、ステップＳ６１へ進む。
【００８４】
ステップＳ５５へ進むと、現在の気泡粒子数ｎ、すなわち、今回の列ｎが、停滞粒子Ｂ’を表す列（停滞列）ｎｓか否かを調べ、ｎ≠ｎｓのときは、ステップＳ７へ戻り、気泡粒子Ｂを、移動距離ΔＺだけ上昇させる。一方、ｎ＝ｎｓのときはステップＳ１１戻る。その結果、停滞列ｎｓにある停滞粒子Ｂ’は、結果的にその後の計算が打ち切られる。
【００８５】
又、ステップＳ１１からステップＳ６１へ進むと、停滞粒子Ｂ’の数を、各停滞している領域毎にカウントする。次いで、ステップＳ６２へ進み、各領域毎のカウント値に応じて、各領域に設定されているマーキングの色（或いは濃淡）を可変設定し、ステップＳ１３へ進む。
【００８６】
停滞粒子Ｂ’の数が多い領域は、エアーポケットが発生し易く領域であるため、当該部位のエアーポケット対策を他の部位に優先して検討する必要がある。本形態では、停滞粒子Ｂ’の数からエアーポケットが発生し易く領域に順番を付け、その順番に従ってマーキングの色（或いは濃淡）を可変設定するようにしたので、エアーポケット対策を行う順番を一目で把握することができ、検討に要する時間を短縮することができる。
【００８７】
尚、本形態においては、最初の発生位置に設定する気泡粒子Ｂの間隔を等間隔とすることで、停滞粒子Ｂ’が集まりやすい領域の順位をより高精度に検出することができる。
【００８８】
又、図２０〜図２３に本発明の第７形態を示す。本形態は第５形態の変形例である。本形態は、第５形態の構成に更に、被塗装ワーク１０２を取り囲む外部領域（任意に設定できる空間領域）を設定し、この外部領域に気泡粒子Ｂが到達し、或いは通過した場合、気泡粒子Ｂの移動計算を打ち切るようにした構成を加えたものである。尚、図２３に示す被塗装ワーク１０２は、第６形態と同一のものである。
【００８９】
図２０〜図２２に示す気泡粒子移動計算ルーチンのステップＳ７１〜Ｓ７５は、第１形態の図２に示すフローチャートのステップＳ１〜Ｓ５に対応している。すなわち、ステップＳ７１で、被塗装ワーク１０２の形状データを入力し、ステップＳ７２で計算領域Ａに充満されている流体の物性値（流体密度、粘性係数等）、及び気泡粒子Ｂの物性値（密度、粘性係数等）をそれぞれ入力する。次いで、ステップＳ７３で重力Ｇの方向を設定する。次いで、ステップＳ７４で計算条件を入力する。計算条件は気泡粒子Ｂの発生個数Ｎ、規定回数Ｍ、気泡粒子Ｂの演算毎の移動距離ΔＺ、及び気泡粒子Ｂの直径Ｄである。その後、ステップＳ５へ進み、Ｎ個の気泡粒子Ｂの発生位置を設定する。
【００９０】
そして、ステップＳ７６へ進み、被塗装ワーク１０２を取り囲む外部領域Ｓを入力装置２（図１参照）に設けられていマウス等を用いて設定する（図２３参照）。
【００９１】
次いで、ステップＳ７７へ進む。ステップＳ７７〜Ｓ８１は、図２、図３に示すフローチャートのステップＳ６〜１０に対応しており、ステップＳ７７で解析が開始されて、ステップＳ７８へ進むと、今回の気泡粒子Ｂを移動距離ΔＺ分、重力（Ｇ）と反対の方向へ移動させる。
【００９２】
そして、ステップＳ７９で、移動した気泡粒子Ｂが外部領域Ｓ内にあるか否かを調べる。気泡粒子Ｂが外部領域Ｓ内にあるか否かは、外部領域Ｓの外郭を設定する位置座標と、気泡粒子Ｂの中心を示す位置座標とを比較して行う。そして、気泡粒子Ｂが外部領域Ｓの外側にある場合は、ステップＳ８８へジャンプする。又、気泡粒子Ｂが外部領域Ｓ内にある場合は、ステップＳ８０へ進む。
【００９３】
ステップＳ８０では、被塗装ワーク１０２の形状を表す線と気泡粒子Ｂとが干渉するか否かが調べられる。気泡粒子Ｂには直径Ｄの領域が設定されているため、この気泡粒子Ｂの外周が、被塗装ワーク１０２の形状を表す線と接触した場合、具体的には、気泡粒子Ｂの外周の位置座標が被塗装ワーク１０２の形状を表す線の位置座標上にあるとき干渉と判断される。干渉と判断された場合、ステップＳ８１へ進み、干渉した気泡粒子Ｂが上昇するか否かを、線の傾斜角θから判定する。そして、気泡粒子Ｂが上昇すると判定した場合、ステップＳ８２へ進み、当該気泡粒子Ｂを被塗装ワーク１０２を表す線に沿って移動（上昇）させる。又、ステップＳ８０で被塗装ワーク１０２に対して非干渉状態にあり、又、ステップＳ８１で気泡粒子Ｂが上昇できないと判定したときは、ステップＳ８３へジャンプする。
【００９４】
そして、ステップＳ８０，Ｓ８１、或いはＳ８２からステップＳ８３へ進むと、今回の気泡粒子Ｂが移動した位置座標を解析モデルデータベース１３（図１参照）に記憶させ、ステップＳ８４へ進む。ステップＳ８４では、今回の気泡粒子Ｂの座標位置と、当該気泡粒子Ｂの前回の座標位置とを比較し、同一か否かを調べる。そして、前回と今回との気泡粒子Ｂの座標位置が相違しているときは、そのままステップＳ８８へ進む。一方、前回と今回との気泡粒子Ｂの座標位置が同一のとき、すなわち、気泡粒子Ｂが被塗装ワーク１０２の表面に干渉し、それ以上移動できない状態のときは、ステップＳ８５へ分岐し、停滞粒子設定処理を実行して、ステップＳ８８へ進む。
【００９５】
この停滞粒子設定処理は、第５形態の図１５に示す停滞粒子設定処理サブルーチンに従って行われる。尚、この停滞粒子設定処理サブルーチンは既に説明したので、ここでの説明は省略する。
【００９６】
ステップＳ８８では、気泡粒子数ｎが発生個数Ｎに達したか否かを調べ、達していない場合（ｎ＜Ｎ）、ステップＳ８９へ分岐し、気泡粒子数ｎをインクリメントして（ｎ←ｎ＋１）、ステップＳ９０へ進む。又、達した場合（ｎ＝Ｎ）、ステップＳ９１へ進む。
【００９７】
ステップＳ９０へ進むと、現在の気泡粒子数ｎ、すなわち、今回の列ｎが、停滞粒子Ｂ’を表す列（停滞列）ｎｓかが調べられ、ｎ≠ｎｓのときは、ステップＳ７８へ戻り、気泡粒子Ｂを移動距離ΔＺだけ上昇させる。一方、ｎ＝ｎｓのときはステップＳ８８戻る。
【００９８】
その結果、図２３に示すように、気泡粒子Ｂは外部領域Ｓ内でのみ移動計算が行われ、又、この外部領域Ｓ内でき気泡粒子Ｂが停滞粒子Ｂ’に設定された場合は、結果的に計算が打ち切られる。従って、図２３に示すように、気泡粒子Ｂの開始位置は外部領域Ｓの下端の外郭となり、終了位置は外部領域Ｓの上端の外郭となる。従って、本形態では、例えば気泡粒子Ｂが外部領域Ｓの上端側外郭に達した場合、計算回数ｍが規定回数Ｍに達しなくても移動計算は結果として打ち切られることになる。
【００９９】
このように、本形態では、第５形態による効果に加え、外部領域Ｓで計算領域を限定するようにしたので、エアーポケットの発生に関係しない気泡粒子Ｂの移動計算を行う必要が無くなり、その分、計算効率が向上し、計算時間を短縮することができる。
【０１００】
又、図２４、図２５に本発明の第８形態を示す。本形態は第１形態の変形例である。第１形態では、重力（Ｇ）方向を最初に設定するようにしているが、本形態では、最初に設定した重力（Ｇ）方向を、演算途中であっても変更できるようにしたものである。
【０１０１】
図２４に示すルーチンは、図３に示す気泡粒子移動計算ルーチンのステップＳ１５からステップＳ７へ戻るルーチン間に、ステップＳ１０１，Ｓ１０２を追加したもので、その他の処理は、図２、図３と同一である。
【０１０２】
すなわち、図３に示す気泡粒子移動計算ルーチンのステップＳ１５で気泡粒子数ｎをクリアした後（ｎ←０）、ステップＳ１０１へ進み、重力（Ｇ）方向が変更されているか否かを調べる。重力（Ｇ）方向は、キーボードやマウスを用い、表示装置３に表示されている画像を参照しながら変更する。
【０１０３】
そして、重力（Ｇ）方向が変更された場合、ステップＳ１０２へ分岐し、指定された方向に重力（Ｇ）方向を変更し、ステップＳ７へ戻る。一方、重力（Ｇ）方向が変更されていない場合は、そのままステップＳ７へ戻る。
【０１０４】
図２５（ａ）には、重力（Ｇ）方向を変更する前の状態が示されている。同図に示されている被塗装ワーク１０２は開口部１０２ｂのほぼ対角線上に他の開口部１０２ｇが形成されている。このような形状の被塗装ワーク１０２は、図２５（ａ）に示すように、開口部１０２ｂを斜め下方に指向させた姿勢で計算領域Ａに配設した場合、内部に浸入した気泡粒子Ｂは、上部に停滞し易くなる。そこで、同図（ｂ）に示すように、気泡粒子Ｂが上昇している過程で、例えばＰ行目の気泡粒子Ｂから、この気泡粒子Ｂの移動方向を他の開口部１０２ｇ側に設定する。換言すれば、Ｐ行目移行の気泡粒子Ｂに対する重力（Ｇ）方向を他の開口部１０２ｇと反対方向へ指向するように変更する。
【０１０５】
その結果、Ｐ行目の気泡粒子Ｂから、移動方向が他の開口部１０２ｇ側へ設定されるため、気泡粒子Ｂが被塗装ワーク１０２の内部に停滞せず、他の開口部１０２ｇから抜けさせることができる。以上のことから、図２５（ａ）に示すような形状の被塗装ワーク１０２のエアーポケット対策としては、被塗装ワーク１０２を電着槽に浸漬した後、この被塗装ワーク１０２を、他の開口部１０２ｇが上方へ指向するような姿勢に回転させことで、エアーポケットの発生を回避することができる。
【０１０６】
又、例えば被塗装ワーク１０２を電着槽に投入する際の前後左右方向への振れを、重力（Ｇ）方向を変更することで考慮することもできる。
【０１０７】
尚、本発明は、上述した各形態に限るものではなく、例えば、計算領域Ａに解析用メッシュを生成し、重力（Ｇ）方向に変えて、この重力（Ｇ）方向とは反対方向へ電着塗料の流速を設定するようにし、この電着塗料の流速から気泡粒子Ｂの移動計算を行うようにしても良い。
【図面の簡単な説明】
【０１０８】
【図１】第１形態による電着塗装解析装置の機能ブロック図
【図２】同、気泡粒子移動計算ルーチンを示すフローチャート（その１）
【図３】同、気泡粒子移動計算ルーチンを示すフローチャート（その２）
【図４】同、気泡粒子移動計算の結果を示す説明図
【図５】同、気泡粒子の移動状態を示す説明図
【図６】同、気泡移動計算結果を説明する要部拡大図
【図７】同、気泡粒子移動計算を用いて被塗装ワークを電着槽から引き上げる際の気泡粒子の移動計算の結果を示す説明図
【図８】同、気泡粒子の下方への移動状態を示す説明図
【図９】第２形態による気泡粒子移動計算ルーチンを示すフローチャート
【図１０】同、気泡粒子移動計算の結果を示す説明図
【図１１】第３形態による気泡粒子移動計算ルーチンを示すフローチャート
【図１２】同、（ａ）は気泡粒子の発生位置を設定する状態の説明図、（ｂ）は要部拡大図
【図１３】第４形態による気泡粒子移動計算ルーチンを示すフローチャート
【図１４】第５形態による気泡粒子移動計算ルーチンを示すフローチャート
【図１５】同、停滞粒子設定処理サブルーチンを示すフローチャート
【図１６】同、エアーポケットの発生位置を示す説明図
【図１７】第６形態による気泡粒子移動計算ルーチンを示すフローチャート
【図１８】同、停滞粒子設定処理サブルーチンを示すフローチャート
【図１９】同、（ａ）は気泡粒子移動計算の結果を示す説明図、（ｂ）はフランジ部の拡大図、（ｃ）は上部の拡大図
【図２０】第７形態による気泡粒子移動計算ルーチンを示すフローチャート（その１）
【図２１】同、気泡粒子移動計算ルーチンを示すフローチャート（その２）
【図２２】同、気泡粒子移動計算ルーチンを示すフローチャート（その３）
【図２３】同、気泡粒子移動計算の結果を示す説明図
【図２４】第８形態による気泡粒子移動計算ルーチンを示すフローチャート
【図２５】同、（ａ）は重力方向を下方に設定した状態の気泡粒子移動計算の結果を示す説明図、（ｂ）は途中で重力方向を変更した状態の気泡粒子移動計算の結果を示す説明図
【符号の説明】
【０１０９】
１…電着塗装解析装置、２…入力装置、３…表示装置、１１…解析モデル設定処理部、１２…モデル解析処理部、１３…解析モデルデータベース、１４…表示処理部、ΔＺ…移動距離、θ…傾斜角、ρ，ρｏ…密度、１０２…被塗装ワーク、１０２ａ…内部空間領域、１０２ｂ…開口部、１０２ｃ…通気孔、１０２ｄ…表面、１０２ｆ…フランジ部、１０２ｇ…開口部、Ａ…計算領域、Ｂ…気泡粒子、Ｂ’…停滞粒子、Ｄ…直径、Ｇ…重力、Ｌ…幅、Ｍ…規定回数、Ｎ…発生個数、Ｓ…外部領域、ｍ…計算回数、ｎ…気泡粒子数、ｎｓ…停滞列
代理人    弁理士    伊  藤    進
【図１】
2006116385.tif 000003
【図２】
2006116385.tif 000004
【図３】
2006116385.tif 000005
【図４】
2006116385.tif 000006
【図５】
2006116385.tif 000007
【図６】
2006116385.tif 000008
【図７】
2006116385.tif 000009
【図８】
2006116385.tif 000010
【図９】
2006116385.tif 000011
【図１０】
2006116385.tif 000012
【図１１】
2006116385.tif 000013
【図１２】
2006116385.tif 000014
【図１３】
2006116385.tif 000015
【図１４】
2006116385.tif 000016
【図１５】
2006116385.tif 000017
【図１６】
2006116385.tif 000018
【図１７】
2006116385.tif 000019
【図１８】
2006116385.tif 000020
【図１９】
2006116385.tif 000021
【図２０】
2006116385.tif 000022
【図２１】
2006116385.tif 000023
【図２２】
2006116385.tif 000024
【図２３】
2006116385.tif 000025
【図２４】
2006116385.tif 000026
【図２５】
2006116385.tif 000027