"ワークを処理液に沈めたときの空気溜まりを予測するシミュレーション方法であって、前記ワークの解析モデルを構成する複数の要素から、端に位置する要素を処理液に接する境界要素として設定する境界設定ステップと、前記ワークの解析モデルを構成する複数の要素から、前記境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定する気液判定ステップとを有し、前記気液判定ステップは、前記境界要素と前記隣接要素との間のワーク面に、所定の付着力で付着する気泡モデルを設定する気泡設定ステップと、前記気泡モデルに作用する浮力に基づいて、前記ワーク面に平行となる方向の浮力分力を算出する浮力分力算出ステップと、前記気泡モデルに作用する処理液の流動推力に基づいて、前記ワーク面に平行となる方向の推力分力を算出する推力分力算出ステップと、前記浮力分力と前記推力分力とを合算し、前記気泡モデルに対して前記ワーク面に平行となる方向に作用する合力を算出する合力算出ステップと、前記合力が前記付着力を上回る場合には前記隣接要素が処理液に接すると判定する一方、前記合力が前記付着力を下回る場合には前記隣接要素が空気に接すると判定する判定ステップとを有することを特徴とするシミュレーション方法。請求項1記載のシミュレーション方法において、前記隣接要素が処理液に接すると判定された場合に、この隣接要素を新たな境界要素として設定する境界更新ステップを有し、前記気液判定ステップは、新たに設定された境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定することを特徴とするシミュレーション方法。請求項1または2記載のシミュレーション方法において、前記気泡モデルに作用する処理液の流動推力は、前記解析モデルの座標毎に処理液の流れ方向および流速が格納されるデータベースを参照して求められることを特徴とするシミュレーション方法。請求項1~3のいずれか1項に記載のシミュレーション方法において、前記解析モデルを構成する要素は、前記ワークの表面形状を分割する複数の表面要素、または前記表面要素の頂点に設けられる複数の節点要素であることを特徴とするシミュレーション方法。ワークを処理液に沈めたときの空気溜まりを予測するシミュレーションプログラムであって、前記ワークの解析モデルを構成する複数の要素から、端に位置する要素を処理液に接する境界要素として設定する境界設定ステップと、前記ワークの解析モデルを構成する複数の要素から、前記境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定する気液判定ステップとを有し、前記気液判定ステップは、前記境界要素と前記隣接要素との間のワーク面に、所定の付着力で付着する気泡モデルを設定する気泡設定ステップと、前記気泡モデルに作用する浮力に基づいて、前記ワーク面に平行となる方向の浮力分力を算出する浮力分力算出ステップと、前記気泡モデルに作用する処理液の流動推力に基づいて、前記ワーク面に平行となる方向の推力分力を算出する推力分力算出ステップと、前記浮力分力と前記推力分力とを合算し、前記気泡モデルに対して前記ワーク面に平行となる方向に作用する合力を算出する合力算出ステップと、前記合力が前記付着力を上回る場合には前記隣接要素が処理液に接すると判定する一方、前記合力が前記付着力を下回る場合には前記隣接要素が空気に接すると判定する判定ステップとを有することを特徴とするシミュレーションプログラム。請求項5記載のシミュレーションプログラムにおいて、前記隣接要素が処理液に接すると判定された場合に、この隣接要素を新たな境界要素として設定する境界更新ステップを有し、前記気液判定ステップは、新たに設定された境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定することを特徴とするシミュレーションプログラム。請求項5または6記載のシミュレーションプログラムにおいて、前記気泡モデルに作用する処理液の流動推力は、前記解析モデルの座標毎に処理液の流れ方向および流速が格納されるデータベースを参照して求められることを特徴とするシミュレーションプログラム。請求項5~7のいずれか1項に記載のシミュレーションプログラムにおいて、前記解析モデルを構成する要素は、前記ワークの表面形状を分割する複数の表面要素、または前記表面要素の頂点に設けられる複数の節点要素であることを特徴とするシミュレーションプログラム。特許文献1のシミュレーション方法は、処理槽内における処理液の流動状況と空気溜まりの発生状況とをデータベース化し、このデータベースを参照して空気溜まりの発生状況を予測している。このように、処理液の流動状況を考慮したものであるため、高精度に空気溜まりを予測することが可能であるが、個々のワーク形状に合わせて処理液の流動状況と空気溜まりの発生状況とをデータベース化することは、作業コストを増大させる要因となっていた。本発明の目的は、作業コストを抑制しつつ空気溜まりの予測精度を高めることにある。本発明のシミュレーション方法は、ワークを処理液に沈めたときの空気溜まりを予測するシミュレーション方法であって、前記ワークの解析モデルを構成する複数の要素から、端に位置する要素を処理液に接する境界要素として設定する境界設定ステップと、前記ワークの解析モデルを構成する複数の要素から、前記境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定する気液判定ステップとを有し、前記気液判定ステップは、前記境界要素と前記隣接要素との間のワーク面に、所定の付着力で付着する気泡モデルを設定する気泡設定ステップと、前記気泡モデルに作用する浮力に基づいて、前記ワーク面に平行となる方向の浮力分力を算出する浮力分力算出ステップと、前記気泡モデルに作用する処理液の流動推力に基づいて、前記ワーク面に平行となる方向の推力分力を算出する推力分力算出ステップと、前記浮力分力と前記推力分力とを合算し、前記気泡モデルに対して前記ワーク面に平行となる方向に作用する合力を算出する合力算出ステップと、前記合力が前記付着力を上回る場合には前記隣接要素が処理液に接すると判定する一方、前記合力が前記付着力を下回る場合には前記隣接要素が空気に接すると判定する判定ステップとを有することを特徴とする。本発明のシミュレーション方法は、前記隣接要素が処理液に接すると判定された場合に、この隣接要素を新たな境界要素として設定する境界更新ステップを有し、前記気液判定ステップは、新たに設定された境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定することを特徴とする。本発明のシミュレーション方法は、前記気泡モデルに作用する処理液の流動推力は、前記解析モデルの座標毎に処理液の流れ方向および流速が格納されるデータベースを参照して求められることを特徴とする。本発明のシミュレーション方法は、前記解析モデルを構成する要素は、前記ワークの表面形状を分割する複数の表面要素、または前記表面要素の頂点に設けられる複数の節点要素であることを特徴とする。本発明のシミュレーションプログラムは、ワークを処理液に沈めたときの空気溜まりを予測するシミュレーションプログラムであって、前記ワークの解析モデルを構成する複数の要素から、端に位置する要素を処理液に接する境界要素として設定する境界設定ステップと、前記ワークの解析モデルを構成する複数の要素から、前記境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定する気液判定ステップとを有し、前記気液判定ステップは、前記境界要素と前記隣接要素との間のワーク面に、所定の付着力で付着する気泡モデルを設定する気泡設定ステップと、前記気泡モデルに作用する浮力に基づいて、前記ワーク面に平行となる方向の浮力分力を算出する浮力分力算出ステップと、前記気泡モデルに作用する処理液の流動推力に基づいて、前記ワーク面に平行となる方向の推力分力を算出する推力分力算出ステップと、前記浮力分力と前記推力分力とを合算し、前記気泡モデルに対して前記ワーク面に平行となる方向に作用する合力を算出する合力算出ステップと、前記合力が前記付着力を上回る場合には前記隣接要素が処理液に接すると判定する一方、前記合力が前記付着力を下回る場合には前記隣接要素が空気に接すると判定する判定ステップとを有することを特徴とする。本発明のシミュレーションプログラムは、前記隣接要素が処理液に接すると判定された場合に、この隣接要素を新たな境界要素として設定する境界更新ステップを有し、前記気液判定ステップは、新たに設定された境界要素に隣接する隣接要素を抽出し、この隣接要素が処理液と空気とのいずれに接するかを判定することを特徴とする。本発明のシミュレーションプログラムは、前記気泡モデルに作用する処理液の流動推力は、前記解析モデルの座標毎に処理液の流れ方向および流速が格納されるデータベースを参照して求められることを特徴とする。本発明のシミュレーションプログラムは、前記解析モデルを構成する要素は、前記ワークの表面形状を分割する複数の表面要素、または前記表面要素の頂点に設けられる複数の節点要素であることを特徴とする。本発明によれば、ワーク面に所定の付着力で付着する気泡モデルを設定し、気泡モデルに作用する浮力からワーク面に平行となる浮力分力を算出し、気泡モデルに作用する処理液の流動推力からワーク面に平行となる推力分力を算出する。そして、浮力分力と推力分力とに基づき気泡モデルに作用する合力を算出し、この合力と付着力とを比較することで気液判定を行っている。これにより、処理液の流動推力を考慮することができ、空気溜まりの予測精度を向上させることが可能となる。しかも、処理液の流動推力による空気溜まりへの影響について、気泡モデルを用いて演算するようにしたので、大規模なデータベースを準備する必要がなく、作業コストを大幅に引き下げることが可能となる。以下、本発明の実施の形態を図面に基づいて詳細に説明する。図1は電着塗装工程を示す概略図である。図1に示すように、ワークである車体10に対して電着塗装を施すため、電着塗装工程には電着液を溜めた処理槽11が設置されている。車体10に対して脱脂や表面調整等の前処理を施した後に、車体10を処理槽11に沈めて通電を施すことにより、車体10のアウタ面やインナ面には電着塗膜が形成されることになる。この電着塗装工程において良好な塗装品質を得るためには、車体10を処理槽11に沈めたときに空気溜まりを発生させることのない車体形状が求められている。このような車体形状を効率良く設計するため、設計段階において空気溜まりの発生状況を予測することが重要となっている。また、処理槽11には図示しないノズルが複数設置されており、これらのノズルから電着液が吐出されている。このため、空気溜まりの発生状況を高精度に予測するためには、処理槽11における電着液の流れを考慮することが必要となっている。以下、本発明の一実施の形態であるシミュレーション方法およびシミュレーションプログラムについて説明する。ここで、図2はシミュレーション装置12を示すブロック図であり、このシミュレーション装置12によってシミュレーション方法やシミュレーションプログラムが実行される。図2に示すように、シミュレーション装置12は、CPUやメモリ等によって構成される演算装置13、キーボード等の入力装置14、液晶ディスプレイ等の表示装置15、磁気ディスク等の記憶装置16を備えている。このシミュレーション装置12は、マイクロコンピュータやパーソナルコンピュータ等の単一のコンピュータを用いて構成しても良く、ネットワークを介して相互に接続される複数のコンピュータを用いて構成しても良い。シミュレーション装置12の記憶装置16には、車体全体を表す車体解析モデルが格納されている。ここで、図3は車体解析モデルを示す概略図であり、図4は記憶装置16に格納される車体解析モデルのデータを示す説明図である。図3の拡大部分に示すように、車体解析モデルは、車体10の表面形状を分割する複数の表面要素と、表面要素の頂点に設けられる節点要素とによって構成されている。この車体解析モデルについては、各節点要素が位置する座標データ等の形で記憶装置16に格納されている。図4に示すように、記憶装置16には、車体解析モデルを構成するパネル部材の番号データ、各パネル部材が備える節点の番号データ、各節点の座標データ等が格納されている。なお、XYZ直交座標系におけるZ軸の方向は、重力方向Gに対して逆向きに設定されている。また、車体解析モデルとしては、衝突変形シミュレーション等に用いられる車体解析モデルを流用することが可能である。また、演算装置13には気液判定部17が設けられており、この気液判定部17は車体解析モデルの各節点について気液判定処理を実行する。この気液判定処理とは、車体解析モデルを構成する各節点について、電着液に接する節点であるか、或いは空気に接する節点であるかを判定する処理である。すなわち、気液判定処理とは、各節点の属性データが電着液または空気の何れであるかを判定する処理である。ここで、図5およびは車体解析モデルIを気液判定する際の手順を示す説明図である。なお、説明を容易にするため、図5には簡略化された車体解析モデルIが示されている。また、図6は図5のA-A線に沿って車体解析モデルIを示す説明図であり、図6は図5のA-A線に沿って車体解析モデルIを示す説明図である。図5および図6に示すように、まず、気液判定処理においては、車体解析モデルIを構成する全節点i1~i15の属性データが空気に設定される。続いて、端に位置する節点i1~i6,i10~i15が境界節点として設定され、図5および図6に示すように、境界節点i1~i6,i10~i15の属性データが電着液に変更される。なお、端に位置する節点とは、車体10の縁の部分に位置する節点、すなわちパネル部材の縁部やパネル部材に形成される穴の縁部に位置する節点を意味し、車体10を処理槽11に沈めたときに無条件で電着液に接する節点である。また、換言すれば、端に位置する節点とは、表面要素によって囲まれていない節点、つまり表面要素が四角形で構成される場合に4つの表面要素によって共有されていない節点である。続いて、図6に示すように、複数の節点から解析対象となる境界節点が設定されるとともに、この境界節点に隣接する隣接節点が設定される。すなわち、節点i6が解析対象となる境界節点として設定されたときには、節点i6に隣接する節点i7が隣接節点として設定される。また、節点i10が解析対象となる境界節点として設定されたときには、節点i10に隣接する節点i9が隣接節点として設定される。このように隣接節点が設定されると、この隣接節点について気液判定が実行される。以下、隣接節点i7の気液判定について説明した後に、隣接節点i9の気液判定について説明する。図7~は隣接節点i7について気液判定を行う際の手順を示す説明図である。図7に示すように、境界節点i6と隣接節点i7との間のワーク面に付着する気泡モデルが設定され、この気泡モデルに作用する浮力Fuが算出される。また、境界節点i6と隣接節点i7との座標値に基づいて、水平方向に対するワーク面の傾斜角αが算出される。そして、浮力Fuと傾斜角αとに基づいて、ワーク面に平行となる方向の浮力分力Fupが算出される。この浮力分力Fupとは、気泡モデルを隣接節点i7に向けて移動させようとする力である。また、図7に示すFuvとは、ワーク面に垂直となる方向の浮力分力である。なお、浮力Fuを算出する際には、例えば、以下の式を用いることが可能であり、浮力分力Fupを算出する際には、例えば、以下の式を用いることが可能である。ρwは電着液の密度であり、ρaは空気の密度であり、Vは気泡モデルの体積であり、gは重力加速度である。Fu=×V×g...Fup=Fu×sinα...続いて、図7に示すように、隣接節点i7の座標値に基づき所定のデータベースを参照することにより、隣接節点i7に作用する電着液の流れ方向および流速を読み込み、気泡モデルに作用する電着液の流動推力Flが算出される。ここで、図8は流動推力Flを算出する際に参照されるデータベースの一例を示す概略図である。図8に示すように、車体解析モデルおよびその周辺が複数の領域に区画されるとともに、これらの領域毎に電着液の流れ方向および流速を示すデータD11~DnNが格納されている。このようなデータベースを参照し、電着液が気泡モデルを押し込む流動推力Flが算出されると、流動推力Flと傾斜角αとに基づいてワーク面に平行となる方向の推力分力Flpが算出される。この推力分力Flpとは、気泡モデルを境界節点i6に向けて移動させようとする力である。また、図7に示すFlvとは、ワーク面に垂直となる方向の推力分力である。なお、流動推力Flを算出する際には、例えば、以下の式を用いることが可能であり、推力分力Flpを算出する際には、例えば、以下の式を用いることが可能である。Sは気泡モデルの受圧面積断面積であり、vwは電着液の流速である。Fl=S×ρw×vw2/2...Flp=Fl×cosα...次いで、図7に示すように、浮力分力Fupと推力分力Flpとに基づいて、気泡モデルに作用する合力が算出される。なお、図示する場合には、推力分力Flpと浮力分力Fupとの作用方向が逆向きであり、浮力分力Fupが推力分力Flpよりも大きいことから、気泡モデルには隣接節点i7に向かう方向に合力が作用している。気泡モデルに作用するスカラー量としての合力は|Fup-Flp|となる。続いて、隣接節点i7の属性データが電着液であるか否かを判定するため、気泡モデルに作用する合力が、ワーク面に付着する気泡モデルの付着力Fsを上回るか否かが判定される。なお、気泡モデルの付着力Fsとは、気泡モデルを移動させる際に移動方向とは反対方向に作用する抗力であり、電着液の表面張力や分子間力等に基づいて予め設定される力である。そして、気泡モデルに作用する合力が付着力を上回る場合には、気泡モデルが移動して抜けることから、隣接節点i9の属性データは電着液であると判定される。一方、気泡モデルに作用する合力が付着力を上回らない場合には、気泡モデルがそのまま停止することから、隣接節点i9の属性データは空気であると判定される。なお、図示する場合には、合力が付着力Fsを下回ることから、隣接節点i7の属性データは空気のまま維持されることになる。続いて、隣接節点i9の気液判定について説明する。図9~は隣接節点i9について気液判定を行う際の手順を示す説明図である。図9に示すように、境界節点i10と隣接節点i9との間のワーク面に付着する気泡モデルが設定され、この気泡モデルに作用する浮力Fuが算出される。また、境界節点i10と隣接節点i9との座標値に基づいて、水平方向に対するワーク面の傾斜角αが算出される。そして、浮力Fuと傾斜角αとに基づいて、ワーク面に平行となる方向の浮力分力Fupが算出される。この浮力分力Fupとは、気泡モデルを境界節点i10に向けて移動させようとする力である。なお、図9に示すFuvとは、ワーク面に垂直となる方向の浮力分力である。続いて、図9に示すように、隣接節点i9の座標値に基づいて図8のデータベースを参照することにより、隣接節点i9に作用する電着液の流れ方向および流速を読み込み、気泡モデルに作用する電着液の流動推力Flが算出される。このように流動推力Flが算出されると、流動推力Flと傾斜角αとに基づいてワーク面に平行となる方向の推力分力Flpが算出される。この推力分力Flpとは、気泡モデルを隣接節点i9に向けて移動させようとする力である。なお、図9に示すFlvとは、ワーク面に垂直となる方向の推力分力である。次いで、図9に示すように、浮力分力Fupと推力分力Flpとに基づいて、気泡モデルに作用する合力が算出される。なお、図示する場合には、推力分力Flpと浮力分力Fupとの作用方向が逆向きであり、推力分力Flpが浮力分力Fupよりも大きいことから、気泡モデルには隣接節点i9に向かう方向に合力が作用している。気泡モデルに作用するスカラー量としての合力は|Fup-Flp|となる。続いて、隣接節点i9の属性データが電着液であるか否かを判定するため、気泡モデルに作用する合力が、ワーク面に付着する気泡モデルの付着力Fsを上回るか否かが判定される。なお、図示する場合には、合力が付着力Fsを上回り、気泡モデルが移動して抜けることから、隣接節点i9の属性データは空気から電着液に変更されることになる。このように、隣接節点i9の属性データが空気から電着液に変更されると、この隣接節点i9は新たな境界節点として設定される。そして、境界節点i9に隣接する節点i8が隣接節点として設定されるとともに、この隣接節点i8について同様の気液判定が実行される。このような気液判定は、未解析の隣接節点が存在しなくなるまで繰り返され、車体解析モデルを構成する全節点の属性データが電着液または空気に分けられることになる。なお、全節点に関する気液判定の結果は、演算装置13のポスト処理部18を経て表示装置15に出力される。ポスト処理部18においては、例えば、電着液に接する節点や空気に接する節点を、色相や濃淡等によって区分して表現する処理が実行される。ここで、図10は流動推力Flを考慮していない従来のシミュレーション方法やシミュレーションプログラムによる車体解析モデルIの気液判定結果を示す説明図である。また、図10は流動推力Flを考慮した本発明のシミュレーション方法やシミュレーションプログラムによる車体解析モデルIの気液判定結果を示す説明図である。また、図11は図10のA-A線に沿って車体解析モデルIを示す説明図であり、図11は図10のA-A線に沿って車体解析モデルIを示す説明図である。まず、図10および図11に示すように、電着液の流れを考慮していない場合、つまり節点間における上下の位置関係だけを考慮した場合には、車体解析モデルIが片側に傾斜する構造であるため、全節点i1~i15の属性データが電着液と判定される。すなわち、車体解析モデルIが下方に開口する凹形状を有していないことから、車体解析モデルIには空気溜まりが発生しないと判定されることになる。一方、図11および図11に示すように、電着液の流れを考慮した場合には、浮力によって上方に抜けようとする空気が流動する電着液によって押し戻されることから、節点i7において空気溜まりが発生するとの正しい結果を導くことが可能となる。以下、前述した空気溜まりの予測手順をフローチャートに沿って簡単に説明する。図12は空気溜まりの予測手順を示すフローチャートである。図12に示すように、ステップS1では、車体解析モデルを構成する全節点の属性データが空気に設定される。続いて、ステップS2では、端に位置する節点が境界節点として設定され、ステップS3では、境界節点の属性データが電着液に変更される。続くステップS4では、後述する隣接節点の有無について未解析の境界節点が存在するか否かが判定される。ステップS4において、未解析の境界節点が存在する場合には、ステップS5に進み、解析対象となる境界節点が設定される。続くステップS6では、境界節点に隣接する隣接節点の有無が判定される。ステップS6において、隣接節点が存在する場合には、ステップS7に進み、隣接節点について気液判定が実施される。ここで、図13はステップS7における気液判定の手順を詳細に示すフローチャートである。図13に示すように、ステップS11では境界節点の座標値が読み込まれ、ステップS12では隣接節点の座標値が読み込まれる。ステップS13では、隣接節点の座標値に基づきデータベースを参照し、電着液の流れ方向および流速が読み込まれる。続くステップS14では、境界節点と隣接節点との間のワーク面に付着する気泡モデルの付着力Fsが読み込まれる。また、ステップS15では、境界節点と隣接節点との座標値に基づきワーク面の傾斜角αが算出される。ステップS16では、気泡モデルに作用する浮力Fuと傾斜角αとに基づき、ワーク面に平行となる方向の浮力分力Fupが算出される。ステップS17では、気泡モデルに作用する流動推力Flと傾斜角αとに基づき、ワーク面に平行となる方向の推力分力Flpが算出される。続いて、ステップS18では、浮力分力Fupと推力分力Flpとに基づき、気泡モデルに作用する合力が算出される。続くステップS19において、気泡モデルに作用する合力が付着力Fsを上回るか否かが判定される。そして、ステップS19において、合力が付着力Fsを上回ると判定された場合には、ステップS20に進み、隣接節点の属性データが電着液に変更される。一方、ステップS19において、合力が付着力Fsを下回ると判定された場合には、ステップS21に進み、隣接節点の属性データが空気のまま維持される。このように、隣接節点についての気液判定が完了すると、図12のステップS8に進み、隣接節点の属性データが電着液であるか否かが判定される。ステップS8において、隣接節点の属性データが電着液であると判定された場合には、ステップS9に進み、電着液に変更された隣接節点が新たな境界節点として設定される。そして、再びステップS4に戻り、未解析の境界節点について気液判定が実行されることになる。一方、ステップS8において、隣接節点の属性データが空気であると判定された場合には、そのままステップS4に戻り、未解析の境界節点について気液判定が実行されることになる。このように、未解析の境界節点が存在しなくなるまで、ステップS4からの気液判定が繰り返され、車体解析モデルを構成する全節点の属性データが電着液と空気とに分けられることになる。これまで説明したように、本発明においては、電着液の流れを考慮して空気溜まりの有無を判断するようにしたので、空気溜まりの予測精度を高めることが可能となる。しかも、電着液の流れによる空気溜まりへの影響について、気泡モデルを用いて演算するようにしたので、大規模なデータベースを準備する必要がなく、作業コストを大幅に引き下げることが可能となる。本発明は前記実施の形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能であることはいうまでもない。例えば、前述の説明では、車体解析モデルの節点要素について気液判定を実行しているが、これに限られることはなく、車体解析モデルの表面要素について気液判定を実行しても良い。この場合には、表面要素の代表点を設定した上で、代表点の座標データを用いて傾斜角αを算出するとともに、座標データを用いて電着液の流れ方向や流速が求められることになる。また、前述の説明では、気泡モデルの直径寸法を2mmに設定しているが、この値に限られることはなく、車体仕様や処理槽仕様等に応じて適宜変更しても良い。また、前述の説明では、ワークとして車体10を挙げているが、これに限られることはなく、ワークとしてケース等の他の部品を用いても良い。さらに、図示する場合には、四角形の表面要素によってパネル部材の表面形状を表しているが、これに限られることはなく、三角形や五角形等の表面要素を用いて解析モデルを構成しても良い。なお、電着塗装工程を例に挙げて説明しているが、本発明のシミュレーション方法やシミュレーションプログラムを、単なるディップ方式の塗装処理に対して適用しても良く、ワークに金属皮膜を形成するメッキ処理に対して適用しても良い。"