JP2008087035A - スプリングバック発生原因部位特定方法、その装置、及びそのプログラム - Google Patents

Abstract

【課題】プレス成形製品のスプリングバックの原因となる部位を特定するという課題があった。
【解決手段】第１の成形条件を数値解析して、第１の成形データを算出する第１のプレス成形解析ステップと、第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出する第２のプレス成形解析ステップと、第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行うステップと、演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップとを有してスプリングバック発生原因部位特定方法を構成する。
【選択図】図３

Description

本発明は、自動車用部材などを鋼板や棒鋼などの鋼材を始め、アルミ、マグネシウム、チタン等の金属からプレス成形した際にプレス成形製品に発生するスプリングバックの原因部位を特定するための方法、その装置、及びプログラムに関する。

ドアやバンパーなど多くの自動車用部材、家電部材、建材等は鋼板のプレス成形により製造されている。近年、それらの部材に対する軽量化の要求が高まっており、その要求を実現するために高強度を有する鋼材を使用することによって鋼材を薄手化するなどの対応が図られている。
しかしながら、鋼板の高強度化に伴い、プレス成形による製品形状の確保には厳しい管理が必要となっており、その管理において重要な項目の一つは、プレス成形時に鋼板に生じた残留応力を駆動力として、鋼板の弾性変形分が弾性回復する変形であるスプリングバックである。

特に、昨今は、自動車等の開発工数及びコスト削減のため、デザイン段階と同時に、その成形部材の成形方法を検討する設計段階が開始される傾向にあり、デザイン段階でのデザイン変更が、設計段階での成形部材の変更を生じさせるため、成形部材の成形方法を検討する設計段階での工数やコストは、自動車等の開発工程や開発費において、より大きな問題となっている。

図１は、従来のスプリングバックへの対応を示す成形部材の断面模式図である。（ａ）は、成形製品の形状断面を示し、（ｂ）は、（ａ）で示した成形製品と同型の金型で鋼板を冷間プレス成形後、成形製品に生じるスプリングバックを示し、（ｃ）は、スプリングバックを想定して補正した金型の形状断面を示す。このように、（ａ）に示される成形製品を得るために、（ｃ）に示すようなスプリングバックを予想した「見込み」の金型によって、所望の成形製品を得る対応がある。

このようなスプリングバックを予想した金型を成形する方法として、有限要素法により金型にプレスされた下死点における鋼板の残留応力を解析し、その残留応力と反対向きの残留応力により生じる変形（スプリングフォワード）形状の金型を数値解析することにより、簡易にスプリングバックを考慮した金型を成形する方法が提案されている（特許文献１、非特許文献１）。

しかし、スプリングバックを完全に考慮した金型を数値解析により設計することは、非線形問題であり非常に困難であるため、提案されている方法は、あくまでも有限要素法によって簡易にスプリングバックを考慮した金型を成形するというものである。そのため、その金型によって、スプリングバックの許容量を満たさない場合にどのような対策が必要となるかについては、数値的に解析することが困難であるため、依然として、何の解決方法も提案されていない。

そのため、簡易にスプリングバックを考慮した金型で、スプリングバックの許容量が満たされなかった場合、所望の成形製品を得るためにどのような対策をとるかは、技術者の経験に依存することになる。結局、その成形法による金型と実際の鋼板とによるトライアンドエラーテストが必要になる。

また、金型形状ではなく、鋼材や成形製品の形状に残留応力を除去する修正を加えることでスプリングバックを少なくする方法が提案されている。
図２は、スプリングバックによる変形の発生原因となる部位を探る従来の方法を例示する斜視図である。（ａ）は、成形製品の形状を示し、（ｂ）は、製品の一部１を切断除去した場合を示し、（ｃ）は、製品に孔２を開けた場合を示し、（ｄ）は、製品の一部にスリット３を付与する場合を示す。このような対策を講じることによってスプリングバックの挙動を観察すると共に、スプリングバックを少なくする対策が試みられている。

また、スプリングバック発生部位への対策により、スプリングバックの原因となる残留応力は低減されても、切断除去や孔開け等のために部材自体の剛性が低下するので、わずかな残留応力でスプリングバックが発生してしまうという問題が生じ、根本的な原因究明には至らない。さらに、このような対策は、実際に試験金型と鋼板とによるテストを必要とするため、設計段階の工数とコストの増大という問題が生じる。

特開２００３−３３８２８号公報 三菱自動車テクニカルレビュー（２００６ Ｎｏ．１８ １２６〜１３１頁）

上述のような問題点に鑑み、本発明は、数値解析によりプレス成形製品のスプリングバック発生の原因となる部位を特定することで、成形部材の成形方法の検討時間を効率的にかつ経済的に短縮する方法を提供することを課題とする。

また、本発明は、製品強度を維持しつつ、スプリングバック量を許容値以下とする成形製品を数値解析により提供することを課題とする。

以上の課題を解決するために、本発明に係るスプリングバック発生原因部位特定方法(プログラム：以下、カッコ内はプログラムの発明に対応する)は、第１の成形条件を数値解析して、第１の成形データを算出する第１のプレス成形解析ステップ(手順)と、第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出する第２のプレス成形解析ステップ(手順)と、第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行うステップ(手順)と、演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップ(手順)と、を有する。

また、上記成形条件は、鋼板の形状及び性状、金型形状、プレス条件であっても良く、上記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値であっても良い。

さらに、本発明に係るスプリングバック発生原因部位特定装置は、第１の成形条件を数値解析して、第１の成形データを算出し、第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出するプレス成形解析部と、成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当領域の第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行い、スプリングバック解析部に、演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、を提供する。

また、上記成形条件は、鋼板の形状及び性状、金型形状、プレス条件であっても良く、上記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値であっても良い。
本発明のプログラムは、上記に説明したプログラムの他に、プレス成形品の成形データと複数領域毎の物性値及び物理量のデータを算出するプレス成形解析プログラムと、該成形データと複数領域毎の物性値及び物理量のデータを数値解析してスプリングバック量とスプリングバック後の複数領域毎の物性値及び物理量のデータを算出する数値解析プログラムは市販のプログラムを使用し、これらのプログラムと、データ入出力が可能であるスプリングバック発生原因部位特定プログラムであって、プレス成形解析プログラムから、第１の成形データを取得する手順と、プレス成形解析プログラムから、第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを取得する手順と、第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の前記第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行う手順と、スプリングバック解析プログラムに、前記演算処理の結果を出力する手順をコンピュータに実行させるものでも良い。市販のスプリングバック解析プログラムは、本発明のプログラムから出力された演算処理データに基づいてスプリングバック量を算出する。

本発明によれば、プレス成形製品のスプリングバック発生の原因として解析対象となる特定部位を変更し、かつ、その特定部位の物性値・物理量を数値演算しながら、スプリングバック量を最小化することにより、スプリングバック発生の原因部位を特定し、かつその原因部位の物性値・物理量を正確に導き出すことを可能とし、又は、解析的なトライアンドエラーによりそれらを確認することを可能とし、これにより、成形部材の成形方法の検討時間を経済的にかつ効率的に短縮する方法を提供する。

さらに、本発明は、特定部位及びその特定部位の物性値及び物理量の少なくとも何れかを変化させることで、実際の成形製品に対して切断除去や孔開けを行わずに、スプリングバック変化量を許容値以下とする特定部位とその物性値及び物理量を求めることが可能である。そのため、このように解析された成形製品は、スプリングバック変化量を許容値以下とすると同時に、剛性等の製品の品質を維持することが確認できるため、成形製品に求められる所定のスプリングバック変化量を維持するために、製品品質を犠牲にするような成形製品の特定部位への切断除去や孔開け等の対策を不要とすることができる。

図３は、本発明の一実施形態によるスプリングバック発生原因部位特定装置の機能構成図を示す。
スプリングバック発生原因部位特定装置１０Ａは、プレス成形解析部１１、スプリングバック解析部１２、物性値・物理量演算処理部１４、ファイル格納部１６、成形条件入力部１８、及びスプリングバック量出力画面１９を有する。

成形条件入力部１８は、後述するプレス成形解析部１１及びスプリングバック解析部１２で解析対象となる鋼板の形状データ（板厚、長さ、幅、曲率、歪等）、性状（強度、伸びなどの材質、板厚などの形状）、金型形状（ダイ（金型）及びパンチ形状、曲率、径、クリアランス、潤滑条件）、プレス条件（しわ押さえ荷重、パッド荷重、ビード張力、プレス圧力、温度）などの成形条件を入力する入力部である。または、成形解析におけるデータ領域、物性値・物理量演算処理部１４におけるデータ領域、スプリングバック量出力画面１９で画面表示する際の分割領域等を別個に設定することも可能である。

プレス成形解析部１１は、成形条件入力部１８で入力された成形条件を入力情報として、弾塑性有限要素法、剛塑性有限要素法、ワンステップ有限要素法、境界要素法、初等解析等を用いて、プレス成形された被加工物である鋼板等の応力や歪や板厚などを求める数値解析を行う。プレス成形解析部１１は、被加工物の板厚、応力の成分値、歪の成分値等の変数や、その変数の分布という形式で数値解析結果を出力する。その出力データは、例えば、ファイル「Ｐ ｏｒｇ．ｋ」として、スプリングバック解析部１２、後述する物性値・物理量演算処理部１４及びファイル格納部１６に出力される。

このプレス成形解析部１１による数値解析は、有限要素法（例えば市販のソフトウェアPAM-STAMP、LS-DYNA、Autoform、OPTRIS、ITAS-３D、ASU/P-FORM、ABAQUS、ANSYS、MARC、HYSTAMP、Hyperform、SIMEX、Fastform3D、Quikstamp）を用いて、プレス成形する製品の形状データ（板厚、長さ、幅、曲率、歪みなど）及び、使用する金属板の性状（強度、伸びなどの材質、板厚などの形状）に基づいて、必要であれば金型形状（ダイ及びパンチ形状、曲率、径、クリアランス、潤滑条件）、プレス条件（温度、圧力）などの成形条件を設定し、成形解析を行い、かつ、成形後の応力、歪値の分布を数値的に得ることができる。
また、プレス成型成形解析部１１は、有限要素法による応力分布や曲率等の解析結果を、ポスト処理ソフトを用いて、スプリングバック量出力画面１９にコンタ表示することが可能である。

スプリングバック解析部１２は、プレス成形解析部１１の出力データファイル「Ｐ ｏｒｇ．ｋ」及び後述する物性値・物理量演算処理部１４の出力データファイル「Ｐ ｒｅｍ．ｃａｓｅｎ．ｋ」を入力データとして用いて、スプリングバック解析を行う。スプリングバック解析とは、プレス成形解析部１１の出力結果である被加工物の板厚、応力の成分値、歪の成分値等の変数、及び変数分布に基づいて、弾性有限要素法，弾塑性有限要素法、ワンステップ有限要素法、初等解析等により、除荷過程の計算を行い、被加工物に生じる変形量であるスプリングバック量を数値解析することである。そのスプリングバック量は、被加工物を有限要素法等により分割し、３次元データ座標の要素毎に計算される。なお、被加工物に生じる変形量であるスプリングバック量としては、被加工物の任意点の変形量、被加工物の指定領域内の最大変位点もしくは最小変位点の変形量、被加工物内の複数の任意の面又は線同士がなす角度、または、被加工物の任意の面又は線のなす曲率等がある。

スプリングバック解析部１２のスプリングバック解析結果である出力データは、スプリングバック量、スプリングバック時の歪等の形状、性状、応力等を含み、スプリングバック量出力画面１９に出力され、かつ、入力データファイル「Ｐ ｏｒｇ．ｋ」による数値解析結果出力データファイル「ＳＢ ｏｒｇ．ｋ」、又は、「Ｐ ｒｅｍ．ｃａｓｅｎ．ｋ」による数値解析結果出力データファイル「ＳＢ ｒｅｍ．ｃａｓｅｎ．ｋ」として、スプリングバック解析部１２、後述する物性値・物理量演算処理部１４及びファイル格納部１６に出力される。

物性値・物理量演算処理部１４は、プレス成形解析部１１の出力データファイル「Ｐ ｏｒｇ．ｋ」、スプリングバック解析部１２の数値解析結果出力データファイル「ＳＢ ｏｒｇ．ｋ」、又は、「ＳＢ ｒｅｍ．ｃａｓｅｎ．ｋ」を入力し、演算処理を行い、その結果として「Ｐ ｒｅｍ．ｃａｓｅｎ．ｋ」をスプリングバック解析部１２に出力する。なお、ここで言う演算処理は、上記変数の少なくとも１個以上を係数倍すること、ゼロを含む一定値にすること、四則演算をすること、関数に基づいて計算すること、一定ではない任意の値に置換することをいう。物性値・物理量演算処理部１４は、このような、演算処理により、スプリングバック量を減少させる変数及び特定部位を判別することを可能にする。

物性値・物理量演算処理部１４は、入力データから被加工体の１つであるプレス成型成形製品のデータを取得すると、そのプレス成形製品の一部から特定部位を発見するため、そのプレス成形製品のデータを複数の領域に分割する処理を行う必要がある。その領域分割方法の１つは、均等な寸法により被加工体の領域を分割することである。

均等な寸法は、特定部位に対する対策を考慮することが望ましい。つまり、特定部位を発見しても、特定部位への有効な対策がとれない極めて微小な領域であっては、解析結果が有効に活かせない場合がある。例えば、特定部位への対策が２０ｍｍ×１００ｍｍのビードを加えることである場合、分割領域は、２０ｍｍ角となることが挙げられる。

成形品の分割領域を決定するための１つの方法として、曲率や要素に基づいて、分割領域を決定する方法もある。
プレス成形解析部１１は、数値解析により、成形解析後の被加工物の幾何学的な情報として各節点の座標値を得、各要素の曲率や要素間の角度を客観的に算出することが可能である。プレス成形品の場合、変形後の被加工物はパンチ肩Ｒやダイ肩Ｒ等の曲げの稜線部では、曲げ稜線に垂直な方向で他の部位に比べ非常に大きな曲率となる。
ここでは各要素の最大曲率を計算し、ある閾値以上の曲率となる要素を非表示にすれば、パンチ肩Ｒやダイ肩Ｒ以外の部位（ウェブ面、縦壁部、フランジ部）を連続していない別々の領域として判別分離することが可能である。

図４〜図８にプレス成形解析から演算処理までの工程の例を示す。
図４がプレス成形解析結果であり、応力分布がコンタ図で表示されている。プレス成形解析結果である各要素の曲率を表示したものが図５であり、曲率の大きい要素を非表示としたのが図６である。段差部で非常に大きな曲率となることから、縦壁部内である閾値以上の曲率となる要素を非表示にすれば、段差の稜線部を境界とした複数の領域に分離可能であり、個別の領域を抽出または選択することで特徴形状に基づいた領域を選択可能である。図６により、段差のある側壁は３つの面と、ウェブ、段差のない縦壁部（図から隠れている）、両側の各フランジ部に分割できることがわかる。図６により分割したそれぞれの領域について、図４に示した成形データ（応力分布）に基づいて、応力の高い領域について演算処理を行うと判定したのが、図７に示す領域Ａ１５〜Ａ１７である。判定工程では、曲率及び／又は角度に基づいて分割した領域について、後述の演算を行うかを成形データに基づいて判定する。例えば、分割した各領域のプレス成形解析後の応力レベルにより判定する。

領域Ａ１５〜Ａ１７について、それぞれの領域の全ての要素に対して、例えば応力をゼロとする演算処理を行った状態が図８である。
ここでは要素の最大曲率を用いたが要素間の角度を用いても同様の領域分割を行うことができる。
また、段差のないハット断面形状のメンバー部材を考えると，分離された各領域として、ウェブ面、両側の縦壁部、両側のフランジ部を選択することができる。また、選択した領域が平面でない場合、同様の手法を繰り返すことで、さらに特徴的な領域を選択していくことが可能である。
このようにして、物性値・物理量演算処理部１４は、プレス成形解析部１１の解析結果より、曲率や要素に基づいて、分割領域を決定することが可能である。

スプリングバック解析部１２は、有限要素法（例えば市販のソフトウェアPAM-STAMP、LS-DYNA、Autoform、OPTRIS、ITAS-3D、ASU/P-FORM、ABAQUS、ANSYS、MARC、HYSTAMP）を用いて、プレス成形解析部１１で得られた「Ｐ ｏｒｇ．ｋ」内に記載の応力分布を入力し、スプリングバック解析を実施する。ソフトウェア内でのスプリングバックの計算は例えば、有限弾塑性変形の基礎式や離散化手法に沿った内容で計算される。ただし、スプリングバックの計算は弾性解析であっても、弾塑性解析であっても良い。

また、スプリングバック解析部１２は、有限要素法によるスプリングバック解析結果を、ポスト処理ソフトを用いて、スプリングバック量出力画面１９にコンタ表示することが可能である。また、物性値・物理量演算処理部１４において、スプリングバック量を、分割領域面積、プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さ等で除算することにより規格化し、除算した変数とスプリングバックの影響度をより判別し易くしてデータ出力し、その出力データを用いて、ポスト処理ソフトにより、コンタ表示することが可能である。このような規格化により、被加工体の物性値・物理量とスプリングバック量との関係を視覚的に捉え易くすることができる。

さらに、物性値・物理量演算処理部１４において、スプリングバック量を、パンチ速度又はしわ押さえ力等のプレス成形条件で除算することで規格化した値を求め、ポスト処理ソフトにより、コンタ表示することで、プレス成形条件と、スプリングバック量との関係を視覚的に捉え易くすることもできる。このようなスプリングバック量の規格化とその数値のコンタ表示は、スプリングバック発生の原因部位の特定を経済的にかつ効率的に短縮し、かつ、成形製品の成形方法の検討時間を経済的にかつ効率的に短縮することが可能とする。

ファイル格納部１６は、プレス成形解析部１１、スプリングバック解析部１２、後述する物性値・物理量演算処理部１４の出力結果であるデータファイル「Ｐ ｏｒｇ．ｋ」、「ＳＢ ｏｒｇ．ｋ」、「Ｐ ｒｅｍ．ｃａｓｅｎ．ｋ」、「ＳＢ ｒｅｍ．ｃａｓｅｎ．ｋ」、「Ｐ ｔｒｉｍ．ｃａｓｅｎ．ｋ」等を保存するための格納部である。しかし、それらのデータファイル等が、プレス成形解析部１１と、スプリングバック解析部１２と、物性値・物理量演算部１４との間で直接入出力される場合は、このファイル格納部１６は、必ずしも必要では無い。

物性値・物理量演算処理部１４は、プレス成形解析部１１の出力データファイル「Ｐ ｏｒｇ．ｋ」内の一部の領域を選択し、その一部の領域の板厚、弾性係数、塑性係数、応力の成分値、歪の成分値の少なくとも一つ以上の変数に対して演算処理を実施する。そして、演算処理後の出力データファイル「Ｐ ｒｅｍ．ｃａｓｅ１．ｋ」を生成して、スプリングバック解析部１２に出力する。このデータファイルの送受信は、繰り返し行うことが可能であり、ファイル格納部１６に逐次保存可能であり、また、ファイル形式ではなく実行プロセス又はスレッドのデータ入出力という形でデータのやりとりを行っても良い。

有限要素法で使用されるスプリングバック特定部位又は原因判断のための領域データは、有限要素法における１個以上の要素、複数の要素からなる微小領域、又は、要素が連続した集合体などの計算単位区分、有限要素法による数値積分手法の計算単位で応力や歪を計算する点のうちの一部である積分点（板厚方向における一部の点、面内の一分の点の双方を含む）等とすることができる。「１個以上の要素」とは有限要素法の領域分割単位である要素をいい、「計算単位区分」とは初等解析における計算単位となる微小領域の一つまたは連続した集合体をいう。また、「積分点」とは、一般的に有限要素法で行う近似的な積分を行う点である。「一部の積分点」とは、有限要素法による数値積分手法の計算単位で応力や歪を計算する点のうちの一部をいい、板厚方向における一部の点、及び面内の一部の点の双方を含むものとする。

また、「板厚、弾性係数、塑性係数、応力の成分値、歪の成分値の少なくとも一つ以上の変数」とは、被加工物の一部の板厚、弾性係数（ヤング率、ポアソン比）、塑性係数（降伏応力、塑性係数、加工硬化指数）、解析する３次元座標系（ｘ、ｙ、ｚ）としたときに、ｘ、ｙ、ｚ方向の応力若しくは歪（それぞれ３成分）、ｘｙ平面、ｙｚ平面、ｚｘ平面の剪断応力若しくは剪断歪（それぞれ３成分）、の計１８の因子のうち少なくとも一つ以上をいう。

図９は、有限要素法で使用される座標系を示す図である。（a）は、３次元全体の座標系を、（ｂ）は、局所座標系を示す。
このとき、図９（ａ）で示すように、(ｘ、ｙ、ｚ)座標軸での成分値の他に、例えば同図（ｂ）に示すように、要素３１の平面をｘｙ平面とみなした局所座標系（ｘ１、ｘ２、ｘ３）における応力、歪、剪断応力、剪断歪の計１２成分も含むものとする。また、相当応力、相当塑性歪、弾性エネルギー（弾性仕事）、塑性エネルギー（塑性仕事）等の各応力成分や歪成分を使って得られる変数や、応力増分やひずみ増分等の各成分値の増分形で計算された成分値も含むものとする。

一般に、鋼板を解析対象とする有限要素法解析では、その鋼板の物性値・物理量は、２ｍｍ〜４ｍｍ角の有限要素に分割される。しかし、要素の分割単位は、物性値・物理量が一定と近似される長さで分割されるため、必ずしも２ｍｍ〜４ｍｍ角に限定されるものでは無い。つまり、残留応力が大きい部位は、有限要素をより小さく限定することが必要な場合もある。また、各要素は、３次元座標面で規定され、有限要素の面の角度や曲率は、他の面との比較によって規定される。

このように、スプリングバック発生原因部位の特定のために、角度や曲率に基づいて部材を領域分割し、残留応力レベルの高い部分を演算対象とすることで、計算量を大幅に減らすことにより、より迅速なスプリングバック発生原因部位の特定が可能である。

演算処理の例として、σ_x成分値を係数倍する方法を説明する。
選択された領域の積分点における演算前の応力成分を（σ_x0，σ_y0，σ_z0，τ_xy0，τ_yz0，τ_zx0）、歪成分を（ε_x0，ε_y0，ε_z0，γ_xy0，γ_yz0，γ_zx0）とすると、演算処理後の応力成分（σ_x，σ_y，σ_z，τ_xy，τ_yz，τ_zx）、歪成分（ε_x，ε_y，ε_z，γ_xy，γ_yz，γ_zx）は以下のとおりとなる。
σ_x＝Ｋ×σ_x0，σ_y＝σ_y0，σ_z＝σ_z0，τ_xy＝τ_xy0，τ_yz＝τ_yz0，τ_zx＝τ_zx0
ε_x＝ε_x0，ε_y＝ε_y0，ε_z＝ε_z0，γ_xy＝γ_xy0，γ_yz＝γ_yz0，γ_zx＝γ_zx0
ここでＫは−１０００〜１０００までの範囲で変化させた演算を行い、板厚は演算処理前の値をそのまま使うこととし、弾性係数、塑性係数は、成形条件入力部１８で入力した値を使う場合が一般的である。例示のため、Ｋは、σ_xに関してのみ示すが、他の成分についても、同様に変化させても良い。

また、材料特性（板厚、弾性係数、塑性係数）を演算対象とすることも可能である。例えば、物性値・物理量演算処理部１４において、プレス成形解析部１１から得られた被加工物の一部の領域を選択し、選択された領域の例えば、ヤング率を演算処理において係数倍する。その場合は選択された領域の演算前の板厚ｔ₀、ヤング率Ｅ₀、ポアソン比ν₀、塑性係数Ｆ₀、ａ₀、ｎ₀（σ＝Ｆ₀（ε＋ａ₀）ⁿ⁰の場合）とすると、演算後の板厚ｔ、ヤング率Ｅ、ポアソン比ν、塑性係数Ｆ、ａ、ｎ（σ＝Ｆ（ε＋ａ）ⁿの場合）は以下のとおりとなる。
ｔ＝ｔ₀，Ｅ＝Ｋ×Ｅ₀，ν＝ν₀，Ｆ＝Ｆ₀，ａ＝ａ₀，ｎ＝ｎ₀
ここでＫは−１０００〜１０００の範囲で変化させることができる。例示のため、Ｋは、Ｅに関してのみ示すが、他の材料特性についても、同様に変化させても良い。

また、物性値・物理量演算処理部１４では、プレス成形解析部１１から得られた状態変数を出力したファイルを直接編集により変更することで領域を選択して演算処理することも可能である。例えばファイルの内容をワードパット等のテキスト編集機能を有するソフトで表示し、成形条件入力部１８の操作により演算処理を行いたい領域の成分値を直接書き換えたり、カットアンドペーストにより成分値を入れ替えても良い。
これにより、選択した領域の演算した変数（成分）がスプリングバックにどのくらい影響を与えているかを定量的に評価できる。

例えばスプリングバック解析部１２で有限要素節点の変位分布をコンタ図（等高線図）で表示した場合や変形形状の切断面を表示した場合は、選択した領域の演算した変数（成分）ごとの結果の図を、スプリングバック量出力画面１９上に並べて比較することや図示しないプリンタに出力して比較することが可能である。また、任意の指定位置の節点の座表値（Ｘ、Ｙ、Ｚ）をソフト上またはファイル出力により確認した場合は、選択した領域の演算した変数（成分）ごとの数値を比較したり、表計算ソフト等を用いてグラフ化して比較することが可能である。

スプリングバック量出力画面１９上に表示されるコンタ図は、物性値・物理量演算処理部１４による物性値・物理量演算処理、及び、スプリングバック解析部１２によるスプリングバック量計算を繰り返し行いながら、スプリングバックの原因となる部位や変数を厳密に特定していくことが出来る。

さらに、物性値・物理量演算処理部１４の特定領域を選択により変更したり、選択された領域の成分値をσｘだけではなくσｙやτｘｙなどの成分値についても演算処理として係数倍してから、スプリングバック解析部１２のスプリングバック解析を実施し、スプリングバック量出力画面１８にスプリングバックの原因となる部位や変数を表示させ、スプリングバックによる変位の変化量を比較することで、スプリングバックの原因となる部位や変数を厳密に特定していくことが可能である。

また、物性値・物理量演算処理部１４において、被加工物の一部の中で複数の領域を選択して、それぞれの選択した領域に対し演算処理を同時に実施することで原因特定のために必要な解析工数を低減することが可能である。

図１０は、本発明の一実施形態によるハット型断面形状部品の領域選択の例を示す斜視図である。
例えば、プレス成形解析部１１により得られた被加工物が図５に示すようなハット型断面形状部品の場合に、物性値・物理量演算処理部１４において図１０に示すような成形後の応力が所定値より大きい複数の領域Ａ１、Ａ２を同時に選択し、選択した領域Ａ１、Ａ２に対してσｘを係数倍にする演算処理を実施し、スプリングバック解析部１２においてその演算処理結果を用いたスプリングバック解析を行う。

そのスプリングバック解析結果から、スプリングバックによる変位の変化量を調べて、領域Ａ１、Ａ２に加わっている応力成分値σｘがどのくらいスプリングバックに影響しているかを、片側ずつ演算処理をしなくても１回の演算処理のみで評価できる。

さらにスプリングバック解析部１２のスプリングバック解析を行った後に最もスプリングバック量が変化するときの領域を再分割して、再分割した領域の大きさが所定値以下になるまで、判定して、再分割したそれぞれの領域について、物性値・物理量演算処理部１４において演算処理を行い、再度スプリングバック解析部１２でスプリングバック解析を行うことで、最もスプリングバック量が変化する領域をスプリングバック発生の原因を特定することができ、効率的な原因特定が可能である。なお、物性値・物理量演算処理部１４は、スプリングバック変化量が許容値以下となるように、物性値及び物理量を操作変数として変化させるような収束計算を行うことができ、同時に、特定領域の分割領域も操作変数として変化させるような収束計算を行っても良い。

また、プレス成形解析部１１において、例えば有限要素法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部１４において領域を選択する際に１個以上の要素を選択して演算対象の領域とすることができる。場の問題を解くための解析手法である有限要素法においては、連続体を有限の要素に分割して扱うのが特徴であり、要素はその辺上に位置する有限個の節点において連結され、個々の要素に定義された形状関数と節点の変位から要素内の変位場が一意的に決定されるものである。例えば、プレス成形解析部１１で得られた応力成分σｘの分布を表示し、σｘ値が最も大きい要素を選択領域としたり、σｘが最も大きい要素とその要素に接する要素の集団を選択領域とすることができる。

また、プレス成形解析部１１において、例えば初等解析手法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部１４において領域を選択する際に１個以上の計算単位区分を選択して演算対象の領域とすることができる。

また、プレス成形解析部１１において、例えば有限要素法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部１４において領域を選択する際に１個以上の積分点を選択して演算対象の領域とすることができる。例えば、プレス成形解析部１１の有限要素解析において被加工物は板厚方向に５個の積分点（板裏面から表面に向けて１、２、３、４、５）を有するシェル要素で解析しておき、物性値・物理量演算処理部１４において選択した要素の板裏面に最も近い１と板表面に最も近い５の積分点のみを演算対象の領域とすることができる。これにより、金型の曲率を有する部分に接触して変形した際に受ける曲げ変形の影響を分離した評価を行える。

図１１は、本発明の一実施例によるスプリングバック原因部位特定処理のスプリングバック発生原因特定処理のフローチャートを示す。

ステップＳ１０１では、成形条件入力部１８に成形条件が入力される。次に、ステップＳ１０２に進む。

ステップＳ１０２では、プレス成形解析部１１が、成形条件で規定された被加工物に対して数値解析処理を行い、プレス成形する製品の成形後の応力、歪値の分布を計算する。次に、ステップＳ１０３に進む。

ステップＳ１０３では、物性値・物理量演算処理部１４で、変換対象領域となる特定部位が決定される。この特定部位は、プレス成形品について曲率及び／又は角度を算定し、前記曲率及び／又は角度に基づいて成形品の領域を分割し、分割したそれぞれの領域について前記成形データに基づいて判定して、後述の演算処理を行う領域を決定する。次に、ステップＳ１０４に進む。

ステップＳ１０４では、Ｓ１０３で演算処理を行うと判定した一部の領域における物性値及び物理量の少なくとも一つを変換する上述の演算処理が行われる。次に、ステップＳ１０５に進む。

ステップＳ１０５では、プレス成形時の被加工体の変換処理された物性値及び物理量に基づいて、スプリングバック量が計算される。次に、ステップＳ１０６に進み、同時にステップＳ１０７でスプリングバック量の結果がコンタ表示などにより画面表示される。

ステップＳ１０６では、スプリングバックの許容値が許容値以内か判定を行う。許容値以内であれば、本処理を終了する。許容値以上であれば、ステップＳ１０８に進む。なお、ステップＳ１０６及びステップＳ１０８は、人間が、特定部位の指定を、スプリングバック変化量を見ながら手動で行っても良い。

ステップＳ１０８では、特定部位及び被加工体の物性値・物理量を変更する。ここでは、特定部位のみの変更を行ってもよく、物性値・物理量の変更だけを行っても良く、それら双方を同時に行っても良い。そして、ステップＳ１０５に戻る。ステップＳ１０５、ステップＳ１０６、ステップＳ１０８の処理は、スプリングバック変化量が許容値以下になるような収束計算により行うことが可能で、閾値以下になるまで繰り返される。

なお、本処理の繰り返し回数を制限し、制限した繰り返し回数を超えた場合は、処理を終了しても良い。その場合、ファイル格納部１６に格納されたスプリングバック解析部１２の出力データファイル「ＳＢ ｒｅｍ．ｃａｓｅｎ．ｋ」を人間が検討し、人間がスプリングバック量出力画面１９を見ながら、最適な特定部位を探索して求めることが可能となる。

図１２は、本発明の一実施例によるスプリングバック原因部位特定処理のスプリングバック発生原因特定処理のフローチャートを示す。

ステップＳ２０１では、成形条件入力部１８に第１の成形条件が入力される。次に、ステップＳ２０２に進む。

ステップＳ２０２では、プレス成形解析部１１が、第１の成形条件で規定されたプレス成形品に対して数値解析処理を行い、プレス成形する製品の成形後の応力、歪値の分布を計算する。さらに、プレス成形解析部１１は、上記第１の成形条件のすくなくとも１つを変更した第２の成形条件で規定された被加工物に対して数値解析を行い、プレス成形する製品の成形後の応力、歪値の分布（物性値及び物理量）を計算する。次に、ステップＳ２０３に進む。

ステップＳ２０３では、物性値・物理量演算処理部１４で、変換対象領域となる特定部位が決定される。この特定部位は、第１の成形条件で算出されたプレス成形品の領域毎の物性値及び物理量と、第２の成形条件で算出されたプレス成形品の領域毎の物性値及び物理量との差分を計算し、その差分値が、所定値より大きいか否かを判断する。所定値より大きい差分値が有った場合、その差分値を有する領域を特定部位と決定し、ステップＳ２０４に進む。

ステップＳ２０４では、Ｓ２０３で決定した領域における物性値及び物理量の少なくとも一つを変換する演算処理が行われる。次に、ステップＳ２０５に進む。

ステップＳ２０５では、演算処理された物性値及び物理量に基づいて、スプリングバック量が計算される。次に、ステップＳ２０６に進み、同時にステップＳ２０７でスプリングバック量の結果がコンタ表示などにより画面表示される。

ステップＳ２０６では、スプリングバックの許容値が許容値以内か判定を行う。許容値以内であれば、本処理を終了する。許容値以上であれば、ステップＳ２０８に進む。

ステップＳ２０８では、ステップＳ２０２の第２の成形条件において、第１の成形条件と相違させる条件を変更する場合、ステップＳ２０２に戻り、ステップＳ２０３で特定部位を決定するための物性値・物理量や、判断のための所定値を変更する場合は、ステップＳ２０３に戻り、ステップＳ２０４の演算処理を変更する場合は、ステップＳ２０４に戻る。なお、ステップＳ２０８は、画面出力等の結果を見ながら、人間が、判断して行っても良い。ステップＳ２０２〜Ｓ２０８の処理は、スプリングバック変化量が許容値以下になるような収束計算により、スプリングバック許容値以下になるまで繰り返しても良い。

なお、本処理の繰り返し回数を制限し、制限した繰り返し回数を超えた場合は、処理を終了しても良い。その場合、ファイル格納部１６に格納されたスプリングバック解析部１２の出力データファイル「ＳＢ ｒｅｍ．ｃａｓｅｎ．ｋ」を人間が検討し、人間がスプリングバック量出力画面１９を見ながら、最適な特定部位を探索して求めても良い。

図１３は、本発明の一実施形態によるスプリングバック発生原因部位特定装置のハードウェア構成図を示す。上述したプレス成形解析部１１、スプリングバック解析部１２、物性値・物理量演算処理部１４、部分残留応力除去処理部２１、における各処理は、プログラム１００に規定され、そのプログラム１００をコンピュータ９０に実行させても良い。コンピュータ９０は、必要な処理を実行するＣＰＵ９１、処理結果を格納するメモリ９２（例えば、ＲＡＭ（Random Access Memory））、ディスプレイ９３、例えば、キーボードやマウスのような入力装置９４、ハードディスク９５、ＣＤ／ＤＶＤドライブのような外部記憶装置９６、ＮＩＣ（ネットワーク・インタフェース・カード）９７、プリンタ９８を備える。コンピュータ９０は、ＮＩＣ９７に接続されたイーサネット（登録商標）ケーブルにより構成されるネットワーク９９を介して、他のコンピュータ９０Ａと接続可能である。

プログラム１００は、記録媒体に保存され、外部記憶装置９６からローディングされ、もしくは、他のコンピュータ９０Ａからネットワーク９９を介してダウンロードされ、ＣＰＵ９１の制御によって、コンピュータ９０のハードディスク９５に保存される。次に、保存されたプログラム１００は、ＣＰＵ９１によって、実行され、実行プロセス又はスレッドとしてメモリ９２に格納される。例えば、プレス成形解析部１１、スプリングバック解析部１２、物性値・物理量演算処理部１４又は部分残留応力除去処理部２１における各処理が、それぞれ実行プロセス又はスレッドととなり、各実行プロセス又はスレッド間で上述したデータファイル、又は、データが入出力される。また、それらの各実行プロセス又はスレッドは、他のコンピュータ９０Ａに分散して存在し、各処理がコンピュータ９０と他のコンピュータ９０Ａによって分散処理されても良い。

また、図３に示される成形条件入力部１８、スプリングバック量出力画面１９は、それぞれ、入力装置９４、ディスプレイ９３であっても良い。上述したファイル格納部１６は、ハードディスク９５であっても良い。プログラム１００は、ハードディスク９５に保存しても良い。そして、ディスプレイ９３に出力された上述のコンタ図は、プリンタ９８に出力することが可能である。

以下、実施例を用いて本発明を具体的に説明する。

図１４は、本発明の実施例によるスプリングバック解析の対象となるハット型断面形状部品を示す図であり、（ａ）はその斜視図、（ｂ）はその断面図であり、ウェブ面Ｗ０、側壁Ｗ１、Ｗ２、フランジＦ１、Ｆ２から構成される。

プレス成形解析部１１及びスプリングバック解析部１２が行う数値解析には有限要素法に基づく市販の板成形シミュレーション解析ソフトウェアＰＡＭ−ＳＴＡＭＰを使用した。スプリングバック解析は静的陰解法による弾性解析を用いた。

成形条件は、金属板の性状として、板厚１．６ｍｍ、引張強さ７８０ＭＰａ級の高強度鋼板のデータを用いた。また、金型（ダイ、パンチ、ホルダー）の形状をシェル要素によりモデリングし、剛体と仮定して解析した。ダイとパンチのクリアランスは板厚と同じ１．６ｍｍとした。摩擦係数は０．１５を入力し、しわ押さえ荷重として６００ｋＮを設定した。

図１５は、本発明の実施例によるハット型断面形状部品の断面位置とねじれ角度を示す図であり、（ａ）はその斜視図、（ｂ）はその断面図である。
得られた結果をポスト処理ソフトにて図１１に示す位置Ａ、Ｂでの断面を構成する点の座標値を取得し、断面のウェブ面のなす角度θを求めた。スプリングバック量とは、ここでは、角度θとする。

図１６は、第１及び第２の被加工物のプレス成形後の状態量分布を示す図である。第１の被加工物は、第１の成形条件によるプレス成形品であり、第２の被加工物は、第１の成形条件の少なくとも１つの要素が異なる第２の成形条件によるプレス成形品である。図１６では、第２の被加工物の異なる成形条件要素として、形状に変化を加えるために、ビードを加えた。

図１７は、第１及び第２の被加工物のプレス成形後の状態量差分を示す図である。図１８は、状態量差分に基づいた領域選択を示す図である。図１７に示すように、物性値・物理量演算処理プログラムにより、全領域にわたって、第１及び第２の被加工物の応力差分を求め、図１８に示すように、その応力差分の高位１０％を選択領域とした。

物性値・物理量演算処理を実行する物性値・物理量演算処理実行プログラムは、上記選択領域を演算対象領域として分割し、第１の被加工物の選択した各領域に属する要素の全積分点について以下の演算処理を行った。
σ_x＝０，σ_y＝０，σ_z＝０，τ_xy＝０，τ_yz＝０，τ_zx＝０
ε_x＝ε_x0，ε_y＝ε_y0，ε_z＝ε_z0，γ_xy＝γ_xy0，γ_yz＝γ_yz0，γ_zx＝γ_zx0

ここで、選択された領域の積分点における演算前の応力成分を（σ_x0，σ_y0，σ_z0，τ_xy0，τ_yz0，τ_zx0）、歪成分を（ε_x0，ε_y0，ε_z0，γ_xy0，γ_yz0，γ_zx0）と表し、演算処理後における応力成分を（σ_x，σ_y，σ_z，τ_xy，τ_yz，τ_zx）、歪成分を（ε_x，ε_y，ε_z，γ_xy，γ_yz，γ_zx）とした。板厚は、演算処理前の値をそのまま使うこととし、弾性係数、塑性係数はプレス解析処理におけるソフトウェアＰＡＭ−ＳＴＡＭＰに入力した値をそのまま使用した。演算された応力や歪値はファイル出力手段により演算結果ファイルとして出力した。

スプリングバック解析処理として、上述のソフトウェアＰＡＭ−ＳＴＡＭＰを使用した。上記の物性値・物理量演算処理実行プログラムの出力結果を、ソフトウェアＰＡＭ−ＳＴＡＭＰに入力してスプリングバック解析を実施した。スプリングバック解析は、静的陰解法による弾性解析を用いた。

図１５に示す断面Ａに対する断面Ｂのねじれ角をスプリングバック量とした。図１９は、演算処理を行わない場合のスプリングバック量と、演算処理を行った場合のスプリングバック量の比較を示す図である。これにより、図１８で選択した領域への演算処理によりスプリングバック量が低減したことが明らかであり、図１８に示す選択領域がスプリングバック発生原因部位であると特定できた。

以上に例を挙げて示したように、被加工物のスプリングバック影響度を、定量的に分析し、さらにその情報を表示して視覚的に捉えることにより、スプリングバック発生原因部位を容易に特定することが可能となる。
このように、本発明によるスプリングバック発生原因部位特定装置は、従来実際の装置を用いてスプリングバックの検討を、数値解析によりスプリングバック発生原因部位を特定できるので、成形部材の設計段階でのテスト工数及び費用を減少させるものである。また、このようなスプリングバック発生原因部位特定装置は、被加工物全般に適用されることが期待されることから、産業界において多大な利益をもたらす。

以上説明した実施形態は典型例として挙げたに過ぎず、その各実施形態の構成要素を組合せること、その変形及びバリエーションは当業者にとって明らかであり、当業者であれば本発明の原理及び請求の範囲に記載した発明の範囲を逸脱することなく上述の実施形態の種趣の変形を行えることは明らかである。

従来のスプリングバックへの対応を示す成形部材の断面模式図である。 スプリングバックによる変形の発生原因となる部位を探る従来の方法を例示する斜視図である。 本発明の一実施形態によるスプリングバック発生原因部位特定装置の機能構成図を示す。 本発明例によるプレス成形解析結果の応力分布を示すコンタ図である。 本発明例による各要素の曲率を示すコンタ図である。 本発明例による曲率の大きい要素を非表示としたコンタ図である。 本発明例による曲率に基づいて領域分割したことを示す図である。 本発明例による一部の領域の応力をゼロとしたコンタ図である。 本発明の一実施形態による有限要素法で使用される座標系を示す図である。 本発明の一実施形態による被加工物の領域選択の例を示す斜視図である。 本発明の一実施形態によるスプリングバック発生原因部位特定処理のフローチャートである。 本発明の一実施形態によるスプリングバック発生原因部位特定装置の機能構成図を示す。 本発明の一実施形態によるスプリングバック原因部位特定処理のハードウェア構成図である。 本発明の一実施例による被加工物の例を示す斜視図と断面図である。 本発明の一実施例による被加工物の領域分割の例を示す斜視図である。 本発明の一実施例による第１及び第２の被加工物のプレス成形後の状態量分布を示す図である。 本発明の一実施例による第１及び第２の被加工物のプレス成形後の状態量差分を示す図である。 本発明の一実施例による状態量差分に基づいた領域選択を示す図である。 本発明の一実施例による演算処理を行わない場合のスプリングバック量と、演算処理を行った場合のスプリングバック量の比較を示す図である。

符号の説明

１０Ａ，１０Ｂ スプリングバック発生原因部位特定装置
１１ プレス成形解析部
１２ スプリングバック解析部
１４ 物性値・物理量演算処理部
１６ ファイル格納部
１８ 成形条件入力部
１９ スプリングバック量出力画面

Claims (11)

1. 第１の成形条件を数値解析して、第１の成形データを算出する第１のプレス成形解析ステップと、
   前記第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出する第２のプレス成形解析ステップと、
   前記第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の前記第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行うステップと、
   前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップと、
   を有することを特徴とするスプリングバック発生原因部位特定方法。
2. 前記成形条件は、鋼板の形状及び性状、金型形状、プレス条件である請求項１に記載のスプリングバック発生原因部位特定方法。
3. 前記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値である請求項１又は２に記載のスプリングバック発生原因部位特定方法。
4. 第１の成形条件を数値解析して、第１の成形データを算出し、該第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出するプレス成形解析部と、
   前記成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、
   前記第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の前記第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行い、前記スプリングバック解析部に、該演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、
   を有することを特徴とするスプリングバック発生原因部位特定装置。
5. 前記成形条件は、鋼板の形状及び性状、金型形状、プレス条件である請求項４に記載のスプリングバック発生原因部位特定装置。
6. 前記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値である請求項４又は５に記載のスプリングバック発生原因部位特定装置。
7. 第１の成形条件を数値解析して、第１の成形データを算出する第１のプレス成形解析手順と、
   前記第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを算出する第２のプレス成形解析手順と、
   前記第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の前記第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行う手順と、
   前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出手順と、
   をコンピュータに実行させることを特徴とするスプリングバック発生原因部位特定プログラム。
8. プレス成形の成形条件を数値解析して、プレス成形品の成形データと複数領域毎の物性値及び物理量のデータを算出するプレス成形解析プログラムと、該成形データと複数領域毎の物性値及び物理量のデータを数値解析してスプリングバック量とスプリングバック後の複数領域毎の物性値及び物理量のデータを算出する数値解析プログラムと、データ入出力が可能であるスプリングバック発生原因部位特定プログラムであって、
   前記プレス成形解析プログラムから、第１の成形データを取得する手順と、
   前記プレス成形解析プログラムから、前記第１の成形条件の少なくとも１つが異なる第２の成形条件を数値解析して、第２の成形データを取得する手順と、
   前記第１及び第２の成形データの複数領域毎の物性値及び物理量のデータのうち、物性値及び物理量のデータの少なくとも１つの差分が、所定値より大きい領域がある場合、当該領域の前記第１の成形データの物性値及び物理量のデータの少なくとも１つに対して演算処理を行う手順と、
   前記スプリングバック解析プログラムに、前記演算処理の結果を出力する手順
   をコンピュータに実行させることを特徴とするスプリングバック発生原因部位特定プログラム。
9. 前記成形条件は、鋼板の形状及び性状、金型形状、プレス条件である請求項９に記載のスプリングバック発生原因部位特定プログラム。
10. 前記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値である請求項７〜９の何れか１項に記載のスプリングバック発生原因部位特定プログラム。
11. 請求項７〜１０の何れか１項に記載のコンピュータプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。