有限要素法に基づいたシミュレーション解析であって、板厚方向応力を考慮したシェル要素によるプレス成形解析に引き続き、スプリングバック解析を実施する際に、シェル要素の節点は固定した上で板厚方向応力を変化させてその値を0にするとともに、その変化に伴う板厚方向以外の応力の変化量を求め、この求めた応力の変化量で、下死点での応力状態を補正し、この応力状態をもとにスプリングバック解析を行うことを特徴とするプレス成形シミュレーション解析方法。成形対象物に対してプレス成形解析を行い、離型前の形状の応力状態やひずみなどを計算するプレス成形解析演算処理部と、前記プレス成形解析演算処理部の結果をもとに、離型前の形状の応力状態のうち、板厚方向応力を0にするとともに面内応力の変化量を計算し、この計算結果をもとに離型前の応力状態を補正する板厚方向応力解放演算処理部と、前記板厚方向応力解放演算処理部の結果をもとに、スプリングバック解析を行い、離型後の形状の応力状態やひずみなどを計算するスプリングバック解析演算処理部と、前記スプリングバック解析演算処理部の結果をもとに、離型後の形状の応力状態などを表示装置に表示する表示処理部とを備えたことを特徴とするプレス成形シミュレーション解析装置。しかし、プレス成形解析の結果を前提としてスプリングバック解析を行う際の、プレス成形解析とスプリングバック解析間における板厚方向応力の取り扱いについて具体的に示したものはない。一般に、スプリングバック解析は離型後の応力状態についての解析であるから、プレス成形解析からスプリングバック解析への橋渡しの際に、下死点における応力状態における板厚方向応力を0にすることが考えられる。例えば、特許文献3において、なお、ここで言う演算処理は、上記変数の少なくとも1個以上を係数倍すること、ゼロを含む一定値にすること、四則演算をすること、関数に基づいて計算すること、一定ではない任意の値に置換することをいう。特許文献3[0033]参照と記載されていることからも、板厚方向応力を0とすることが伺える。また、実際に、ある市販ソフトでは、板厚方向応力を考慮したシェル要素によるプレス成形解析に引き続きスプリングバック解析を実施する際に、板厚方向応力を単に0とする処理がなされているものがある。しかし、このような処理方法では、スプリングバック解析の計算精度が低下するという問題があった。本発明はかかる課題を解決するためになされたものであり、板厚方向応力を考慮した計算精度の高いプレス成形シミュレーション解析方法及び装置を提供することを目的としている。発明者は、このスプリングバック解析の精度が低くなる原因を追究し、その結果、市販ソフトでは橋渡しの際に、プレス成形解析の下死点での応力状態について、板厚方向応力のみを0としただけで、その変化に伴う板厚方向以外の応力の変化を考慮していないことが原因であると知見した。本発明はかかる知見をもとになされたものであり、具体的には以下のような構成を備えてなるものである。1本発明に係るプレス成形シミュレーション解析方法は、有限要素法に基づいたシミュレーション解析であって、板厚方向応力を考慮したシェル要素によるプレス成形解析に引き続き、スプリングバック解析を実施する際に、シェル要素の節点は固定した上で板厚方向応力を変化させてその値を0にするとともに、その変化に伴う板厚方向以外の応力の変化量を求め、この求めた応力の変化量で、下死点での応力状態を補正し、この応力状態をもとにスプリングバック解析を行うことを特徴とする。2本発明に係るプレス成形シミュレーション解析装置は、成形対象物に対してプレス成形解析を行い、離型前の形状の応力状態やひずみなどを計算するプレス成形解析演算処理部と、前記プレス成形解析演算処理部の結果をもとに、離型前の形状の応力状態のうち、板厚方向応力を0にするとともに面内応力の変化量を計算し、この計算結果をもとに離型前の応力状態を補正する板厚方向応力解放演算処理部と、前記板厚方向応力解放演算処理部の結果をもとに、スプリングバック解析を行い、離型後の形状の応力状態やひずみなどを計算するスプリングバック解析演算処理部と、前記スプリングバック解析演算処理部の結果をもとに、離型後の形状の応力状態などを表示装置に表示する表示処理部とを備えたことを特徴とする。本発明に係るプレス成形シミュレーション解析方法によれば、計算精度の高いプレス成形シミュレーション解析を行うことができる。本発明の実施の形態を図面に基づいて説明する。プレス成形シミュレーション解析装置まず、図1に示すブロック図に基づき、プレス成形におけるプレス成形シミュレーション解析装置1以下、単にプレス成形シミュレーション解析装置1というの構成について説明する。本実施の形態に係るプレス成形シミュレーション解析装置1は、PCパーソナルコンピュータによって構成され、表示装置3と入力装置5と主記憶装置7と補助記憶装置9および演算処理部11とを有している。また、演算処理部11には、表示装置3と入力装置5と主記憶装置7および補助記憶装置9が接続され、演算処理部11の指令によって各機能を行う。<表示装置>表示装置3は計算結果の表示等に用いられ、液晶モニター等で構成される。<入力装置>入力装置5は構造体モデルファイルの表示指示、オペレータの条件入力等に用いられ、キーボードやマウス等で構成される。<主記憶装置>主記憶装置7は、演算処理部11で使用するデータの一時保存や演算等に用いられ、RAM等で構成される。<補助記憶装置>補助記憶装置9は、ファイルの記憶等に用いられ、ハードディスク等で構成される。補助記憶装置9内には、少なくとも、成形対象物情報ファイル13等の各種の情報が格納される。<演算処理装置>演算処理部11はPCのCPUによって構成され、以下に説明する各処理部はCPUが所定のプログラムを実行することによって実現される。演算処理部11内には、プレス成形解析演算処理部15と、板厚方向応力解放演算処理部17と、スプリングバック解析演算処理部19と、表示処理部21とを有する。以下、演算処理部11内の各処理部についてさらに詳細に説明する。プレス成形解析演算処理部15は、オペレータから成形対象物情報ファイル13の指定を受け付ける成形対象物情報ファイル13選択画面と、プレス成形解析条件を受け付けるプレス成形解析条件受付画面を有する。また、プレス成形解析演算処理部15は、オペレータから指定された成形対象物情報ファイル13から読み込んだ成形対象物情報と、オペレータから入力されたプレス成形解析条件をもとに、プレス成形解析を行い、演算結果をプレス成形解析データとして主記憶装置7に記憶させる。板厚方向応力解放演算処理部17は、プレス成形解析データを主記憶装置7から読み出し、板厚方向応力解放演算処理を行い、演算結果を板厚方向応力解放補正データとして、主記憶装置7に記憶させる。以下、板厚方向応力解放演算処理についてさらに詳細に説明する。本発明では、スプリングバック解析を行う前に、板厚方向応力解放演算処理において、板厚方向応力を0にするとともに、板厚方向のひずみを変化させ、その変化に伴う板厚方向以外の応力の変化量を求め、その求めた応力の変化量と、板厚方向のひずみの変化量で下死点における応力状態を補正する板厚方向応力解放演算処理を本演算処理部において行う。具体的には、シェル要素の節点を並進回転ともに固定し、下死点状態のZ方向応力を変化させてその値を0にするとともに、板厚方向のひずみを変化させ、その変化に伴うZ方向以外の応力の変化量を求め、この求めた応力の変化量で下死点での応力状態を補正するものである。ここで、X方向、Y方向は面内方向であり、Z方向は板厚方向である。以下、Z方向以外の応力の変化量を求める方法を説明する。まず、Z方向以外の応力の変化量を求めるために、弾性体の応力-ひずみ関係について計算する。ここで、3次元弾性体の一般的な応力-ひずみ関係は、下記の式1に示す通りである。式1におけるεはひずみ、γはせん断ひずみ、σは応力、Eは縦弾性係数、νはポアソン比であり、添え字の1はX方向、2はY方向、3はZ方向を表す。本発明においては、Z方向の応力の解放前後における、ひずみと応力の変化量について求めることから、式1のひずみεと応力σに差分を表すΔをつけて表現すると下記に示す式2のようになる。Z方向以外の応力の変化量は、弾性変形を仮定して、Z方向の応力σ33を絶対値が等しく向きが反対の応力を代入して0にリセットし、かつZ方向のひずみε33のみを変化させて求める。つまり、式2に対して、Z方向のひずみε33のみを変化させ、ε33以外は一定とすることから、Δε11=Δε22=Δγ23=Δγ31=Δγ12=0とし、またZ方向の応力の変化量は、下死点で作用している応力σ33を解放して0にするために、Δσ33=-σ33を代入すれば下記の式3となる。式3をΔσ11、Δσ22、Δσ23、Δσ31、Δσ12、Δε33について解くと、式4~6のようになる。式4および式5は、Z方向の応力を0にリセットすべくZ方向のひずみのみを変化させた場合のZ方向以外の応力の変化量である。式6は、Z方向のひずみの変化量である。上記のように、プレス成形解析データのある点において、σ33を0にリセットするのに伴い、下死点の応力状態σ11およびσ22を式4のΔσ11およびΔσ22で補正し、スプリングバック解析を行うことにより、スプリングバック解析を正確に行うことができる。ちなみに、Δσ23、Δσ31、Δσ12の値は0なのでσ23、σ31、σ12に対して補正は不要である。スプリングバック解析演算処理部19は、上記のようにして求められた板厚方向応力解放補正データを主記憶装置7から読み出し、読み出した板厚方向応力解放補正データをもとにスプリングバック解析を行い、その結果をスプリングバック解析演算結果データとして主記憶装置7に記憶させる。表示処理部21は、主記憶装置7からスプリングバック解析結果データを読み出し、読み出したスプリングバック解析結果データをもとに、表示装置3を用いて、プレス成形シミュレーション解析結果を表示する。プレス成形シミュレーション解析方法次に、図2に示すフローチャートに基づいてプレス成形シミュレーション解析方法の処理の流れを、図3に示す矩形状の板をゆるく山なりに折り曲げた形状の成形対象物31についてプレス成形シミュレーション解析をする場合を例に挙げて説明する。なお、必要な図は適宜参照する。まず、オペレータは入力装置5を用いて、成形対象物情報ファイル選択画面を介して、プレス成形解析演算処理部15に対し、成形対象物情報ファイル13を指示する。プレス成形解析演算処理部15は、オペレータから指示された成形対象物情報ファイル13から成形対象物情報を読み込む。また、オペレータは入力装置5を用いて、プレス成形解析条件受付画面を介して、プレス成形解析演算処理部15に対し、プレス成形解析条件を入力する。プレス成形解析条件としては、例えば、後述する実施例に示すように、図3の成形対象物31の決め押し領域33に対し、300tonfの決め押し荷重をかけるというようなものである。プレス成形解析演算処理部15は、成形対象物情報とプレス成形解析条件をもとに、プレス成形解析を行い、プレス成形解析データデータD1を算出し、主記憶装置7に記憶させるステップS1。次に、板厚方向応力解放演算処理部17は、プレス成形解析データデータD1を主記憶装置7から読み出し、前述したような板厚方向応力解放演算処理を行い、演算結果を板厚方向応力解放補正データデータD2として、主記憶装置7に記憶させるステップS2。次に、スプリングバック解析演算処理部19は、プレス成形解析データデータD1と板厚方向応力解放補正データデータD2を主記憶装置7から読み出し、プレス成形解析データデータD1を板厚方向応力解放補正データデータD2で補正し、その結果に対してスプリングバック解析を行い、演算結果をスプリングバック解析演算結果データデータD3として、主記憶装置7に記憶させるステップS3。次に、表示処理部21は、スプリングバック解析演算結果データデータD3を主記憶装置7から読み出し、表示装置3を用いてスプリングバック解析演算結果を表示するステップS4。本発明の効果を確認するシミュレーションを行ったので、これについて説明する。本実施例におけるプレス成形シミュレーション解析は、成形対象物31の決め押し領域33に対し、300tonfの決め押し荷重をかけるような成形条件によってプレス成形解析を行い、離型前の形状を求める図4b参照。さらに、プレス成形金型離型後の形状をスプリングバック解析を行って求める図4c参照。このような単純な形状である成形対象物31の場合にはスプリングバック後の形状が想定できるので、このような成形対象物31について、市販ソフトを用いてプレス成形シミュレーション解析行った場合と、本発明を適用したプレス成形シミュレーション解析を行った場合と、いずれが想定される形状を再現しているかを比較する。市販ソフトの場合、決め押し荷重をかけるプレス成形解析の後、スプリングバック解析を行うにあたり、板厚方向応力を単に0にしたのみである。プレス成形シミュレーション解析を行った結果について図5を用いて説明する。図5は、各プレス成形シミュレーション解析におけるスプリングバック解析結果について、成形対象物の中心を基準として、スプリングバック量をコンター表示で示した図である。図5aは市販ソフトを用いてプレス成形シミュレーション解析を行った結果であり、図5bは本発明を適用してプレス成形シミュレーション解析を行った結果である。比較例まず、比較例として市販ソフトを用いて前記プレス成形シミュレーション解析を行った結果について考察する。比較例で行ったスプリングバック解析の結果は、図5aに示すように、スプリングバック量は成形対象物31の中心点から複数の楕円を描くように分布している。これは、成形対象物31がスプリングバックによって、成形対象物31の中心点に対して四隅が最も下がり、中心に向かって楕円状に盛り上がっていることを示している。しかしながら、実際に成形対象物31の決め押し領域33に決め押し荷重をかけるような、単純なプレス成形を行えば、離型後、成形対象物31は幅方向に均一にスプリングバック変形することから、比較例のプレス成形シミュレーション解析の結果は、明らかに間違った結果となっている。実施例これに対して、本発明を適用して前記プレス成形シミュレーション解析を行った結果であるスプリングバック量をコンター表示した図5bをみると、全体的に幅方向に均一にスプリングバックが発生しており、実際のプレス成形に近い結果が出ており、前記比較例より明らかに計算精度が向上したことがわかる。また、本実施の形態で行ったプレス成形シミュレーション解析では、板厚を考慮せずにプレス成形シミュレーション解析を行った場合に比べて計算時間は1割増し程度であり、3次元ソリッド要素を使用して計算した場合のように計算時間が膨大となるようなことはなかった。なお、上記実施の形態および実施例では、板厚方向応力σ33を0にする過程では弾性的な挙動が支配的であるため、弾性体の応力-ひずみ関係の式を用いたが、塑性変形を考慮して弾塑性体の応力-ひずみ関係の式を用いてもよい。