板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行うことを特徴とする円筒深絞りの成形シミュレーション方法。Hillの異方性降伏条件式における異方性パラメータを求めることを特徴とする請求項1に記載の円筒深絞りの成形シミュレーション方法。板厚が4mm以上の鋼板の複数の面それぞれから複数の傾きで立方体を切出し、圧縮して、その歪みに基づいて、前記Hillの異方性降伏条件式における異方性パラメータを求めることを特徴とする請求項2に記載の円筒深絞りの成形シミュレーション方法。板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行う手段を備えたことを特徴とする円筒深絞りの成形シミュレーション装置。板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行う処理をコンピュータに実行させるためのプログラム。しかしながら、厚手の鋼板の円筒深絞りのFEM解析においては、従来の薄板のFEM解析に用いられている面内異方性を考慮した解析では形状の予測精度が不十分であった。本発明は上記のような事情に鑑みてなされたもので、円筒深絞りの成形シミュレーションを行う際に、形状の予測精度を高めることを目的とする。本発明の円筒深絞りの成形シミュレーション方法は、板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行うことを特徴とする。また、本発明の円筒深絞りの成形シミュレーション方法の他の特徴とするところは、Hillの異方性降伏条件式における異方性パラメータを求める点にある。その場合、板厚が4mm以上の鋼板の複数の面それぞれから複数の傾きで立方体を切出し、圧縮して、その歪みに基づいて、前記Hillの異方性降伏条件式における異方性パラメータを求める。本発明の円筒深絞りの成形シミュレーション装置は、板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行う手段を備えたことを特徴とする。本発明のプログラムは、板厚が4mm以上の鋼板の円筒深絞りの成形シミュレーションを行う際に、面内異方性に加えて、板厚断面内の異方性も考慮して成形シミュレーションを行う処理をコンピュータに実行させる。本発明によれば、円筒深絞りの成形シミュレーションを行う際に、形状の予測精度を高めることが可能になる。以下、添付図面を参照して、本発明の好適な実施形態について説明する。板厚が10mmのSPH590の円板の円筒深絞り行い、それを成形シミュレーションであるFEM有限要素法:FiniteElementMethodにて解析した。まず、面内異方性のみを考慮したFEM解析を行った。降伏条件式は式1に示すHillの異方性降伏条件式を用いた。面内異方性を考慮するため、式1における異方性パラメータF、G、H及びNは、例えば非特許文献2にあるように以下のように求めることができる。図2に示すように、鋼板11から、圧延方向に対する傾角α=0°、45°、90°を有する試験片12を作り、引張試験を行い、式2で示すLankfordのr値を求める。上述のようにして求めたr値を式3、式4、式5に代入することにより、それぞれ異方性パラメータG、F、Nを得ることができる。その値を表1に示す。面内異方性を考慮したFEM解析では、図3に示すように、圧延方向から45°方向のカップ高さ図4に示すHが大きくなる傾向は再現できたが、カップ高さが実験結果より大きくなった。また、図4に示すように、成形後の絞り品21の4ヶ所a~dで板厚を調査した結果、図5に示すように、パンチ肩部での板厚が実験結果よりかなり小さくなる結果となった。本願発明者は鋭意検討を重ねた結果、面内異方性のみならず、板厚断面内の異方性が形状変化に影響することを突き止めた。面内異方性に加えて、板厚断面内の異方性も考慮するため、式1における異方性パラメータF、G、H、L、M及びNを、例えば非特許文献3にあるように以下のように求める。図6に示すように、辺長が6mmの立方体を加工し、圧延方向に対する傾角β=0°、45°、90°の方向に圧縮し、式6で示すZ値を求める。図6は、Z値を求めるために用いる立方体試験片の切出し方向、圧縮方向及び歪の測定方向を示す図である。31は鋼板である。32は圧縮方向が圧延方向に対する傾き角β=0°になるように加工した辺長が6mmの立方体試験片である。33は圧縮方向が圧延方向に対する傾き角β=45°になるように加工した辺長が6mmの立方体試験片である。34は圧縮方向が圧延方向に対する傾き角β=90°になるように加工した辺長が6mmの立方体試験片である。また、図7に示すように、辺長が6mmの立方体を加工し、圧延方向に対する傾角γ=0°、45°、90°の方向に圧縮し、式7で示すY値を求める。図7は、Y値を求めるために用いる立方体試験片の切出し方向、圧縮方向及び歪の測定方向を示す図である。41は鋼板である。42は圧縮方向が圧延方向に対する傾き角γ=0°になるように加工した辺長が6mmの立方体試験片である。43は圧縮方向が圧延方向に対する傾き角γ=45°になるように加工した辺長が6mmの立方体試験片である。44は圧縮方向が圧延方向に対する傾き角γ=90°になるように加工した辺長が6mmの立方体試験片である。また、図8に示すように、辺長が6mmの立方体を加工し、圧延方向に対する傾角θ=0°、45°、90°の方向に圧縮し、式8で示すX値を求める。図8は、X値を求めるために用いる立方体試験片の切出し方向、圧縮方向及び歪の測定方向を示す図である。51は鋼板である。52は圧縮方向が幅方向に対する傾き角θ=0°になるように加工した辺長が6mmの立方体試験片である。53は圧縮方向が幅方向に対する傾き角θ=45°になるように加工した辺長が6mmの立方体試験片である。54は圧縮方向が幅方向に対する傾き角θ=90°になるように加工した辺長が6mmの立方体試験片である。上述のようにして求めたZ値、Y値、X値を式9、式10、式11、式12、式13に代入して、それぞれ係数G、F、N、M、Lを得ることができる。その値を表2に示す。図9及び図10に、面内異方性に加えて、板厚断面内の異方性も考慮したFEM解析結果を示す。図9は、被加工材の初期板厚が10mmの場合の実験及び計算によるカップ高さを示す特性図である。また、図10は、被加工材の初期板厚が10mmの場合の実験及び計算による板厚を示す特性図である。被加工材の初期板厚が10mmの場合、面内異方性に加えて、板厚断面内の異方性を考慮した解析結果は、面内異方性のみを考慮した解析結果と比較して実験結果に近く、解析精度が向上している。一方、図11は、被加工材の初期板厚が3.5mmの場合の実験及び計算によるカップ高さを示す特性図である。また、図12は、被加工材の初期板厚が3.5mmの場合の実験及び計算による板厚を示す特性図である。図11及び図12に示すように、板厚が4mm未満になると、面内異方性に加えて、板厚断面内の異方性を考慮した解析結果と、従来の面内異方性のみを考慮した解析結果とは略同等となる。従って、板厚が4mm以上の鋼板の円筒深絞りのFEM解析を行う場合に、面内異方性に加えて、板厚断面内の異方性を考慮した解析を行うのが好適である。板厚が10mm、8mm、4mmのSPH590の鋼板から円板に切削加工したものを被加工材とした。図13に示す金型パンチ61、ダイ62を用い、表3に示す条件にて被加工材63の絞り成形を行った。板厚が厚いため、しわ押さえが無くても皺が発生しないため、しわ押さえを使用しなかった。FEM解析条件を表3に示す。符号No.1及びNo.2は、板厚10mmのFEM解析であり、その結果を図14及び図15に示す。図14は、被加工材の初期板厚を10mmとしたNo.1、No.2及び実験によるカップ高さを示す特性図である。また、図15は、被加工材の初期板厚を10mmとしたNo.1、No.2及び実験による板厚を示す特性図である。本発明を適用したNo.1では、面内異方性のみを考慮した比較例No.2と比較して解析精度が良い。No.3及びNo.4は、板厚8mmのFEM解析であり、その結果を図16及び図17に示す。図16は、被加工材の初期板厚を8mmとしたNo.3、No.4及び実験によるカップ高さを示す特性図である。また、図17は、被加工材の初期板厚を8mmとしたNo.3、No.4及び実験による板厚を示す特性図である。本発明を適用したNo.3では、面内異方性のみを考慮した比較例No.4と比較して解析精度が良い。No.5及びNo.6は、板厚4mmのFEM解析であり、その結果を図18及び図19に示す。図18は、被加工材の初期板厚を4mmとしたNo.5、No.6及び実験によるカップ高さを示す特性図である。また、図19は、被加工材の初期板厚を4mmとしたNo.5、No.6及び実験による板厚を示す特性図である。本発明を適用したNo.5では、面内異方性のみを考慮した比較例No.6と比較して解析精度が良い。上述したような円筒深絞りのFEM解析は、CPU、ROM、RAM等を備えたコンピュータ装置からなる成形シミュレーション装置により実行可能である。また、上述した円筒深絞りのFEM解析を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ又はCPUやMPU等がプログラムを読み出して実行することも可能である。