(57)【要約】【課題】スプリングバックの要因となる部位を高精度かつ効率的に特定することができるプレス成形解析方法を提供する。【解決手段】本発明に係るプレス成形解析方法は、コンピュータが、複数の要素で構成されるプレス成形品モデルにおける一つ又は複数の応力変更領域を選定し、プレス成形品モデルの離型前の応力分布に対して応力変更領域の応力の値を変更してプレス成形品モデルのスプリングバック量を算出し、算出したスプリングバック量に基づいてスプリングバックの要因となる部位を特定するものであって、応力変更領域の選定の基準となる基準要素を選択し、選択した基準要素との間で所定の条件を満たす要素を領域要素として選出し、選出した領域要素に近接する要素のうち、領域要素及び基準要素との間に所定の条件を満たす要素を領域要素として追加して選出し、選択した基準要素と選出した領域要素とからなる領域を応力変更領域として選定するものである。【選択図】  図１
2020011277.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータが、複数の要素で構成されるプレス成形品モデルにおける一つ又は複数の領域を応力変更領域として選定し、前記プレス成形品モデルの離型前の応力分布に対して前記応力変更領域の応力の値を変更し、該変更した応力分布に基づいて前記プレス成形品モデルのスプリングバック量を算出し、該算出したスプリングバック量に基づいてスプリングバックの要因となる部位を特定するプレス成形解析方法であって、
前記応力変更領域は、以下のステップＳ１〜ステップＳ１５の手順に従って選定することを特徴とするプレス成形解析方法。
（Ｓ１）前記プレス成形品モデルの全ての要素を母集団とする。
（Ｓ３）前記応力変更領域の選定に用いる応力の方向を指定する。
（Ｓ５）前記プレス成形品モデルの離型前の応力分布における前記指定した方向の応力の値に基づいて、前記母集団の中から前記領域の選定の基準となる基準要素を選択する。
（Ｓ７）該選択した基準要素に近接する要素のうち、該基準要素との間で所定の条件を満たす要素を領域要素として選出する。
（Ｓ９）該選出した領域要素に近接する要素のうち、該領域要素との間に所定の条件を満たし、かつ前記基準要素との間にも所定の条件を満たす要素を領域要素として追加して選出する。
（Ｓ１１）ステップＳ９において追加して選出される領域要素がなくなるまで、ステップＳ９を繰り返し実行する。
（Ｓ１３）前記選択された基準要素と前記選出された領域要素とからなる領域を応力変更領域として選定する。
（Ｓ１５）前記基準要素及び前記領域要素を前記母集団から除外したものを新たな母集団とし、該新たな母集団に含まれる要素がなくなるまで、又は、選定される応力変更領域が所定の個数に達するまで、ステップＳ５からステップＳ１３を繰り返し実行する。
【請求項２】
ステップＳ７及びステップＳ９における前記基準要素との間の所定の条件は、前記基準要素との角度、重心間距離、法線方向距離及び応力比の全てがそれぞれについて予め定めた範囲内にあることとし、
ステップＳ９における前記領域要素との間の所定の条件は、該領域要素との角度が予め定めた範囲内にあることとすることを特徴とする請求項１記載のプレス成形解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プレス成形品の離型後に生じるスプリングバックの要因となる部位を特定するプレス成形解析方法に関する。
【背景技術】
【０００２】
金属板のプレス成形において、成形対象物であるプレス成形品に要求される形状精度を満足するためには、成形後の金型からプレス成形品を取り出す時に発生するスプリングバック量を低減させることが重要である。スプリングバックとは、成形下死点におけるプレス成形品の残留応力が金型からの離型時に解放されることで発生する弾性回復挙動である。そして、スプリングバックによる弾性回復量（スプリングバック量）を効果的に低減させるためには、プレス成形品の下死点における残留応力の内、どの部位の残留応力がどの程度、スプリングバック挙動に影響しているのかを知ることが重要である。
【０００３】
このような観点から、スプリングバックの発生要因となる部位を特定する方法がこれまでに提案されている。例えば、特許文献１には、離型前の成形対象物におけるある領域の残留応力を消去し、当該領域の残留応力の消去後のスプリングバックに関する量を算出し、当該領域の残留応力を消去する前後で前記スプリングバックに関する量がどのように変化するかを算出することにより、スプリングバックの発生要因となる部位を特定する方法が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許４８９４２９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に開示されている方法においては、成形対象物を碁盤の目状に分割し、該分割した部位を前記成形対象物の応力値に基づいて細分化することにより、残留応力を消去する領域を選定するようにしている。しかしながら、このような領域の選定においては成形対象物の全体形状に基づいて人手で碁盤の目状に分割することを前提とする。そのため、成形対象物を細分化した領域毎に応力値を消去してスプリングバック解析を行ってスプリングバックの要因となる部位を特定するためには、人（例えば、解析者）が手動で領域を選定し、該選定した領域における応力の値を変更するので、スプリングバックの要因となる部位を詳細かつ効率的に特定するのに多大な労力と時間を要する場合があった。
【０００６】
本発明は、上記のような課題を解決するためになされたものであり、プレス成形品の離型後に生じるスプリングバックの要因となる部位を特定するに際し、成形下死点における応力の値を変更する領域を効率的に選定することができるプレス成形解析方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
（１）本発明に係るプレス成形解析方法は、コンピュータが、複数の要素で構成されるプレス成形品モデルにおける一つ又は複数の領域を応力変更領域として選定し、前記プレス成形品モデルの離型前の応力分布に対して前記応力変更領域の応力の値を変更し、該変更した応力分布に基づいて前記プレス成形品モデルのスプリングバック量を算出し、該算出したスプリングバック量に基づいてスプリングバックの要因となる部位を特定するものであって、前記応力変更領域は、以下のステップＳ１〜ステップＳ１５の手順に従って選定することを特徴とするものである。
（Ｓ１）前記プレス成形品モデルの全ての要素を母集団とする。
（Ｓ３）前記応力変更領域の選定に用いる応力の方向を指定する。
（Ｓ５）前記プレス成形品モデルの離型前の応力分布における前記指定した方向の応力の値に基づいて、前記母集団の中から前記領域の選定の基準となる基準要素を選択する。
（Ｓ７）該選択した基準要素に近接する要素のうち、該基準要素との間で所定の条件を満たす要素を領域要素として選出する。
（Ｓ９）該選出した領域要素に近接する要素のうち、該領域要素との間に所定の条件を満たし、かつ前記基準要素との間にも所定の条件を満たす要素を領域要素として追加して選出する。
（Ｓ１１）ステップＳ９において追加して選出される領域要素がなくなるまで、ステップＳ９を繰り返し実行する。
（Ｓ１３）前記選択された基準要素と前記選出された領域要素とからなる領域を応力変更領域として選定する。
（Ｓ１５）前記基準要素及び前記領域要素を前記母集団から除外したものを新たな母集団とし、該新たな母集団に含まれる要素がなくなるまで、又は、選定される応力変更領域が所定の個数に達するまで、ステップＳ５からステップＳ１３を繰り返し実行する。
【０００８】
（２）上記（１）に記載のものにおいて、ステップＳ７及びステップＳ９における前記基準要素との間の所定の条件は、前記基準要素との角度、重心間距離、法線方向距離及び応力比の全てがそれぞれについて予め定めた範囲内にあることとし、ステップＳ９における前記領域要素との間の所定の条件は、該領域要素との角度が予め定めた範囲内にあることとすることを特徴とするものである。
【発明の効果】
【０００９】
本発明によれば、成形下死点における応力の値を変更する領域を効率的に選定することができ、スプリングバックの要因となる部位の特定を高精度かつ効率的に行うことができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の形態に係るプレス成形解析方法において、スプリングバックの要因となる部位を特定するために応力を変更する応力変更領域を選定する処理の流れを示す図である。
【図２】本実施の形態に係るプレス成形解析方法において、解析対象とするプレス成形品モデルを示す図である（x軸：長手方向、y軸：幅方向、z軸：高さ方向）。
【図３】本実施の形態に係るプレス成形解析方法において、解析対象とするプレス成形品モデルの離型前の応力分布を示す図である（x軸：長手方向、y軸：幅方向、z軸：高さ方向）。
【図４】本実施の形態に係るプレス成形解析方法において、応力変更領域の選定の基準とする基準要素に基づいて領域要素を選出する処理の流れを示す図である。
【図５】本発明の実施の形態に係るプレス成形解析方法において、基準要素に基づいて選出した領域要素についてさらに追加して領域要素を選出する処理の流れを示す図である。
【図６】本実施の形態に係るプレス成形解析方法において、スプリングバック（ＳＢ）の要因となる部位を特定するための処理の流れを示す図である。
【図７】本実施の形態に係るプレス成形解析方法において、応力変更領域の選定に用いたプレス成形品モデルの離型前における応力分布を示す図である（x軸：長手方向、y軸：幅方向）。
【図８】本実施の形態に係るプレス成形解析方法において、ステップＳ１からステップＳ１５の処理により選定された応力変更領域と、該応力変更領域の応力の値をコンター表示した図である（x軸：長手方向、y軸：幅方向）。
【図９】本実施の形態に係るプレス成形解析方法において、スプリングバック量の算出法を説明する図である（x軸：長手方向、y軸：幅方向、z軸：高さ方向）。
【図１０】本実施の形態に係るプレス成形解析方法において、各応力変更領域における応力の値を変更して算出したスプリングバック影響度の結果である（x軸：長手方向、y軸：幅方向）。
【発明を実施するための形態】
【００１１】
本発明の実施の形態に係るプレス成形解析方法は、コンピュータが解析するものであって、図２に一例として示すような、複数の要素で構成されるプレス成形品モデル１における一つ又は複数の領域を応力変更領域として選定し、図３に示すようなプレス成形品モデル１の離型前（成形下死点）の応力分布に対して前記応力変更領域の応力の値を変更し、該変更した応力分布に基づいてプレス成形品モデル１のスプリングバック量を算出し、該算出したスプリングバック量に基づいてスプリングバックの要因となる部位を特定するものであり、前記応力変更領域は、図１に示すようにステップＳ１〜ステップＳ１５の手順に従って選定するものである。
【００１２】
以下、本実施の形態において対象としたプレス成形品モデル１について説明した後、ステップＳ１〜ステップＳ１５の各ステップについて説明する。
【００１３】
＜プレス成形品モデル＞
図２に示すプレス成形品モデル１は、長手方向に沿って湾曲し、天板部３と縦壁部５とフランジ部７とを有する断面ハット形状のプレス成形品を平面要素で要素分割してモデル化したものである。ここで、図２におけるx軸方向、y軸方向及びz軸方向は、それぞれプレス成形品モデル１の長手方向、幅方向及び高さ方向に相当する。
そして、プレス成形品モデル１には、成形下死点となる離型前の形状と応力分布（残留応力）が与えられている。図３は、プレス成形品モデル１の長手方向の応力の値をコンター表示したものを示す。
【００１４】
プレス成形品モデル１の離型前の形状と応力分布は、プレス成形品モデル１の形状にプレス成形するプレス成形解析を事前に実施したり、あるいは、既に実施したプレス成形解析の解析結果を記憶した解析ファイルを読み込むことにより取得することができる。
なお、プレス成形解析には、市販の有限要素法解析ソフトウェアを用いることができる。
【００１５】
＜ステップＳ１＞
ステップＳ１は、プレス成形品モデル１の全ての要素を母集団とするステップである。母集団に含まれる各要素は、要素形状、要素を構成する節点に関する情報、節点座標や、応力やひずみ等に関する情報を有している。なお、各要素が有する応力に関する情報とは、通常、プレス成形解析により要素毎に求められた応力のことをいい、応力の値は、各要素の板厚方向に設定されている積分点に対して算出されている。
【００１６】
＜ステップＳ３＞
ステップＳ３は、応力変更領域の選定に用いる応力の方向を指定するステップである。本実施の形態では、前記応力変更領域の選定に用いる応力の方向としてプレス成形品モデル１の長手方向を指定する。もっとも、ステップＳ１は、プレス成形品モデル１の全体座標系における各軸方向（例えば、図２中のx軸方向、y軸方向及びz軸方向）や、プレス成形品モデル１を構成する要素の面内方向を選定に用いる応力の方向として指定してもよく、プレス成形品モデル１に生じるスプリングバック量の要因となりうる応力成分に応じて、前記選定に用いる応力の方向を指定すればよい。
【００１７】
＜ステップＳ５＞
ステップＳ５は、プレス成形品モデル１の離型前の応力分布における前記指定した方向の応力の値に基づいて、前記母集団の中から前記領域の選定の基準となる基準要素を選択するステップである。
【００１８】
本実施の形態では、図１に示すように、母集団の中から、ステップＳ３で指定した方向の応力の絶対値が最大となる要素を選択し（Ｓ５１）、該選択した要素の応力の絶対値が予め設定した基準値以上であるかどうかを判定し（Ｓ５３）、該基準値以上の要素を基準要素として選択する（Ｓ５５）。
なお、基準要素の選択に用いる応力の値は、各要素の積分点の平均値や最大値を用いることができるが、板厚方向の任意の位置における積分点の値を用いてもよい。
【００１９】
＜ステップＳ７＞
ステップＳ７は、ステップＳ５において選択した基準要素に近接する要素のうち、該基準要素との間で所定の条件を満たす要素を領域要素として選出するステップである。
【００２０】
ステップＳ７における具体的な処理の一例を図４に示す。
まず、領域要素の候補として、ステップＳ５で選択した基準要素に近接する要素を抽出する（Ｓ７１）。ここで、基準要素に近接する要素とは、要素が平面要素の場合、基準要素の辺を共有する要素のことをいい、以下、近接要素と表記する。
そして、抽出した近接要素が基準要素との間で所定の条件を満たすかどうかを判定する（Ｓ７５〜Ｓ８１）。
【００２１】
本実施の形態では、基準要素との間の所定の条件は、近接要素と基準要素との角度（Ｓ７５）、近接要素と基準要素の重心間距離（Ｓ７７）、近接要素と基準要素の法線方向距離（Ｓ７９）、及び基準要素における前記指定した方向の応力の値を基準とした近接要素における前記指定した方向の応力の値の比（応力比）（Ｓ８１）のそれぞれが予め定めた範囲内にあることとする。そして、これらの全てが予め定めた範囲内にある近接要素を領域要素として選出し（Ｓ８３）、前記所定の条件のいずれか一つでも満たさないものは領域要素として選出しないものとする。
【００２２】
このように、基準要素について抽出した各近接要素についてＳ７５〜Ｓ８１の判定を行い、全ての近接要素について判定を終えた場合（Ｓ７３）、当該基準要素についての領域要素の選出を完了する（Ｓ８５）。
【００２３】
本実施の形態において、基準要素に近接する近接要素のうち領域要素として選出するかどうかの判定に用いた前記所定の条件は、それぞれ以下の理由による。
【００２４】
まず、近接要素と基準要素との角度（Ｓ７５）は、プレス成形品モデル１において角度変化が大きく形状が変化する部位においては、スプリングバックに対する寄与が異なるとするものである。角度変化が大きく形状が変化する部位とは、プレス成形品モデル１における要素サイズにもよるが、例えば、図２に示す天板部３とパンチ肩Ｒ部９や、パンチ肩Ｒ部９と縦壁部５、等が挙げられる。
【００２５】
近接要素と基準要素の重心間距離（Ｓ７７）は、スプリングバックの要因となる部位を詳細に特定するため、基準要素から距離が離れすぎている要素は応力を変更する領域には含めないようにするものである。
【００２６】
近接要素と基準要素の法線方向距離（Ｓ７９）は、近接要素と基準要素との角度が予め定めた範囲内であっても、法線方向距離が所定の値を越えて形状が変化する部位は段差形状であり、スプリングバックに対する寄与が異なるとするものである。
【００２７】
近接要素の応力の値の比（Ｓ８１）は、基準要素と同程度の応力の値を有する部位は、スプリングバックに対して同じように寄与するとみなされることによる。なお、近接要素の応力の値については、前述の基準要素の応力の値と同様、各要素における積分点の平均値や最大値、あるいは板厚方向の任意の位置における積分点の値を用いることができる。
【００２８】
なお、上記のＳ７５〜Ｓ８１での各判定における予め定めた所定の条件の具体的な数値は、プレス成形品モデル１における要素のサイズに応じて適宜設定することができる。
【００２９】
＜ステップＳ９＞
ステップＳ９は、ステップＳ７において選出した領域要素に近接する要素のうち、該領域要素との間の所定の条件を満たし、かつステップＳ５で選択した基準要素との間の所定の条件を満たす要素を領域要素として追加して選出するステップである。領域要素に近接する要素とは、基準要素に近接する近接要素と同様、要素が平面要素の場合においては当該領域要素の辺を共通とする辺を有する要素のことをいう。
【００３０】
ステップＳ９における具体的な処理を図５に示す。
まず、基準要素に近接する近接要素と同様、ステップＳ７で選出した領域要素に近接する近接要素を抽出する（Ｓ９１）。このとき、ステップＳ７において基準要素の近接要素として抽出されて領域要素となった要素は、ステップＳ９１では近接要素として抽出しないものとする。
【００３１】
そして、抽出した各近接要素について、領域要素との間の所定の条件を満たし、かつ前記基準要素との間の所定の条件を満たす要素を領域要素として追加して選出する。
【００３２】
本実施の形態では、領域要素との間の所定の条件としては、近接要素と領域要素との角度が予め定めた範囲内にあることとする（Ｓ９５）。さらに、基準要素との間の所定の条件としては、近接する要素と基準要素との角度（Ｓ９７）、近接する要素と基準要素の重心間距離（Ｓ９９）、近接する要素と基準要素の法線方向距離（Ｓ１０１）、及び基準要素における前記指定した方向の応力の値を基準とした近接する要素における前記指定した方向の応力の値の比（応力比）（Ｓ１０３）のそれぞれが予め定めた範囲内にあることとする。そして、これらの全てが予め定めた範囲内にある近接要素を領域要素として追加して選出し（Ｓ１０５）、前記所定の条件のいずれか一つでも満たさないものは、領域要素として追加して選出しないものとする。
【００３３】
このように、領域要素について抽出した各近接要素についてＳ９５〜Ｓ１０３の判定を行い、全ての近接要素について判定を終えた場合（Ｓ９３）、当該領域要素について追加して領域要素を選出するのを完了する（Ｓ１０７）。
【００３４】
追加して選出する領域要素に関する上記所定の条件のうち、領域要素との角度については、前述の近接要素と基準要素との角度（Ｓ７５）と同様、プレス成形品モデル１において角度が変化して形状が変わる部位においては、スプリングバックに対する寄与が異なるとするものである。
【００３５】
これに対して、近接要素と基準要素との角度（Ｓ９７）については、近接要素と領域要素との角度が小さくても、基準要素から一つ又は複数の領域要素を挟んで連続する基準要素に対する角度が大きくなるように形状が変化する部位において、スプリングバックに対する寄与が異なるとするものである。このような部位として、例えば、図２に示すパンチ肩Ｒ部９を挟んで存在する天板部３と縦壁部５とが挙げられる。
【００３６】
ステップＳ９におけるその他の所定の条件（Ｓ９９、Ｓ１０１及びＳ１０３）については、前述の基準要素と近接要素との間の所定の条件（Ｓ７７、Ｓ７９及びＳ８１）と同様の理由である。また、領域要素に近接する近接要素の応力の値については、前述の基準要素に近接する近接要素の場合と同様、各要素の積分点の平均値や最大値を用いることができるが、板厚方向の任意の位置における積分点の値を用いてもよい。
【００３７】
＜ステップＳ１１＞
ステップＳ１１は、ステップＳ９において追加して選出された領域要素があるかどうかを判定し、追加して選出される領域要素がなくなるまで、ステップＳ９を繰り返し実行するステップである。
【００３８】
ステップＳ９において追加して選出された領域要素がある場合、当該追加して選出された領域要素に近接する近接要素を抽出する（Ｓ９１）。このとき、基準要素又は領域要素の近接要素として既に抽出されて領域要素となった要素は、追加して選出された領域要素に近接する近接要素として抽出しないものとする。そして、抽出された近接要素が領域要素かどうかを判定し（Ｓ９５〜Ｓ１０３）、当該判定を満たすものを領域要素としてさらに追加して選出する（Ｓ１０５）。
【００３９】
一方、ステップＳ９において追加して選出された領域要素がない場合、領域要素を追加して選出する処理を終了する。なお、ステップＳ１１は、追加して選出された領域要素の要素数が所定の個数に達したら、領域要素を追加して選出する処理を終了してもよい。
【００４０】
＜ステップＳ１３＞
ステップＳ１３は、ステップＳ７で選択された基準要素とステップＳ９で選出された領域要素とからなる領域を応力変更領域として選定するステップである。
ここで、応力変更領域として選定する際に、基準要素と領域要素の要素数に制限を設けてもよい。例えば、基準要素と領域要素の要素数が所定の個数以上のものは応力変更領域として選定する。
【００４１】
＜ステップＳ１５＞
ステップＳ１５は、前記基準要素及び前記領域要素を前記母集団から除外し、所定の個数の応力変更領域を選定するまで、又は、母集団に含まれる要素若しくはステップＳ５において選択される基準要素がなくなるまで、ステップＳ５からステップＳ１３の手順を繰り返し実行するステップである。
【００４２】
ステップＳ１５の具体的な処理としては、図１に示すように、選定された応力変更領域が所定の個数に達したか否かを判定し（Ｓ１５１）、所定の個数に達している場合、応力変更領域の選定を終了する。これに対し、応力変更領域が所定の個数に達していない場合、選択された基準要素と選出された領域要素を母集団から除外したものを新たな母集団とする（Ｓ１５３）。そして、新たな母集団に含まれる要素数が0であるかどうかを判定し（Ｓ１５５）、母集団に含まれる要素数が0でなければ新たな母集団についてステップＳ５からステップＳ１３を繰り返して実行し、そうでなければ応力変更領域の選定を終了する。
【００４３】
次に、本実施の形態に係るプレス成形解析方法により、スプリングバックの要因となる部位を特定する処理の流れを、図２に示すプレス成形品モデル１を解析対象とし、図１に示す手順に従って応力変更領域を選定する場合を例として、図６に基づいて説明する。
【００４４】
まず、解析者またはコンピュータにより、スプリングバックの要因部位の特定に用いる解析結果が記憶された解析ファイル（基本ＳＢ解析ファイル）を準備する（Ｓ２０１）。
基本ＳＢ解析ファイルには、要素分割されたプレス成形品モデル１の離型前（成形下死点）における形状や応力分布といった情報が記録されており、これらの情報は、例えば、プレス成形品モデル１のプレス成形解析により取得されているものとする。
【００４５】
次に、解析者またはコンピュータが、応力を変更する応力変更領域の選定及び当該応力変更領域において変更する応力の方向を指定する（Ｓ２０３）。本実施の形態では、図２に示すプレス成形品モデル１の長手方向の応力の値を指定する。図７に、プレス成形品モデル１の離型前における長手方向の応力分布を示す。
【００４６】
次に、コンピュータが、指定した方向の応力の値に基づいて、応力変更領域を選定する（Ｓ２０５）。応力変更領域の選定には、前述の図１に示す各処理（Ｓ１〜Ｓ１５）を実行する。
応力変更領域の選定においては、基準要素を選択するステップＳ５と、領域要素として選出する判定を行うステップＳ７及びステップＳ９における所定の条件は、以下に示す例のように設定した。
【００４７】
まず、ステップＳ５における基準要素の選択に対する所定の条件は、応力の絶対値が200MPa以上とし（Ｓ５３）、応力の絶対値がこの範囲外の要素は基準要素として選択しないものとする。
【００４８】
次に、基準要素に近接する要素について、該基準要素との所定の条件は、基準要素との角度が15°以内（Ｓ７５）、基準要素との重心間距離が20mm以内（Ｓ７７）、基準要素との法線方向距離が1.5mm以内（Ｓ７９）、応力の絶対値が100MPa以上であって基準要素の応力絶対値に対する比が0.4以上（Ｓ８１）とする。
【００４９】
さらに、選出された領域要素に近接する要素について、該領域要素との所定の条件は、領域要素との角度が5°以内（Ｓ９５）とし、基準要素との所定の条件は、上記と同様、基準要素との角度が15°以内（Ｓ９７）、基準要素との重心間距離が20mm以内（Ｓ９９）、基準要素との法線方向距離が1.5mm以内（Ｓ１０１）、応力の値が100MPa以上であって基準要素の応力絶対値に対する比が0.4以上（Ｓ１０３）とする。
【００５０】
さらに、領域要素の選出においては、選出された領域要素の最小要素数を10とし、選出された領域要素の数がこれに満たない場合、当該基準要素について選定される応力変更領域はないものとする。
さらに、一つの基準要素に対して選出される領域要素の最大要素数は100とし、選出される領域要素が最大要素数に達した時点で、その基準となる基準要素に対して領域要素の選出を終了する。
【００５１】
また、本実施の形態において、応力変更領域の選定にかかるステップＳ３で指定する方向は、スプリングバック要因部位の特定にかかるステップＳ２０３で指定する方向と同一とする。もっとも、本発明は、ステップＳ３において指定する方向とステップＳ２０３において指定する方向とが同一であることを必須の要件とするものではない。
【００５２】
なお、選定された応力変更領域は、コンピュータが、プレス成形品モデル１にマッピングして表示させてもよい（Ｓ２０７）。図８に、一例として、選定された応力変更領域をプレス成形品モデル１にマッピングして表示した結果を示す。
【００５３】
ステップＳ２０５で選定した各応力変更領域について変更した応力の値を入力する応力変更入力ファイルを作成する（Ｓ２０９）。ここでは、応力変更領域の各要素において、ステップＳ２０３で指定した方向の応力を消去（応力の値をゼロ）とする。
【００５４】
もっとも、ステップＳ２０５は、応力変更領域における要素の応力を消去することに限るものではなく指定した方向の応力の値を定数倍する、定数を加算する、定数乗する、板厚方向の平均値又は中央値に置き替える、などのいずれかであってもよい。また、応力は要素の板厚方向における複数の積分点について与えられているので、各要素のうち所定の積分点について応力を変更してもよい。
【００５５】
次に、コンピュータが、応力変更領域の応力を変更したプレス成形品モデル１についてスプリングバック解析を行い（Ｓ２１１）、スプリングバック解析結果を取得する（Ｓ２１３）。当該スプリングバック解析と解析結果の取得は、応力変更領域の応力を変更したものそれぞれについて実行する。
【００５６】
そして、コンピュータが、スプリングバック解析の解析結果を用いて、各応力変更領域の応力を変更したときのスプリングバック影響度を求める（Ｓ２１５）。
スプリングバック影響度として、例えば、図９に示すように、プレス成形品モデル１の長手方向における一端側を固定したときのスプリングバックによる他端側のねじれ角を求める。この場合、応力変更領域の応力を変更せずにスプリングバック解析を実行したときのねじれ角を基準とし、応力変更領域の応力を変更してスプリングバック解析したときのねじれ角の変化量がスプリングバック影響度を表す指標となる。
【００５７】
その他、スプリングバック影響度を求めるにあたり、プレス成形品モデル１において指定した節点の変位や、指定した範囲における変位の二乗和や絶対値平均を用いることができる。さらに、例えば、４つの節点座標から算出される角度変化により、パンチ肩Ｒ部９のスプリングバックによる縦壁部の壁開き、３つ節点座標から曲率変化により、縦壁部５の壁そり等に対するスプリングバック影響度を求めることができる。このように、節点座標（位置情報）からスプリングバック量を計算できるものであれば、スプリングバック影響度は上記のものに限られない。
【００５８】
続いて、コンピュータが、スプリングバック影響度を応力変更領域の要素データとして書き出し（Ｓ２１７）、その書き出した結果をプレス成形品モデル１にコンター表示する（Ｓ２１９）。図１０に、応力変更領域と、各応力変更領域のスプリングバック影響度をコンター表示した結果を示す。
【００５９】
以上、本発明に係るプレス成形解析方法においては、プレス成形品モデル１におけるスプリングバックの要因となる部位を特定するに際し、プレス成形品モデル１の形状と応力分布に基づいて応力を変更する領域の選定することができ、スプリングバックの要因となる部位の特定を高精度かつ効率的に行うことができる。
きる。
【００６０】
さらに、本発明に係るプレス成形解析方法によれば、解析者は、人手またはコンピュータを用いて前記基本ＳＢ解析ファイルを準備し、応力の方向を指定するだけで、応力変更領域の選定のみならず、該選定した応力変更領域毎にスプリングバック影響度を求め、該求めたスプリングバック影響度をプレス成形品モデルにマッピング表示したりグラフ作成に至るまでの一連の過程をコンピュータが自動でスクリプト処理することができる。これにより、スプリングバックの要因部位の特定のさらなる効率化と省力化を図ることができる。
【００６１】
なお、各応力変更領域において応力を変更したときのスプリングバック影響度（例えば、ねじれ角度の変化量）は、解析結果のポスト処理ソフトウェアを用い、マクロ処理により自動的に算出することができる。これらのスクリプト処理やマクロ処理には、市販の解析ソフトウェアを使用することができる。
【００６２】
また、応力変更領域の選定のために指定する方向は、前述のようにプレス成形品モデルの長手方向といった特定の方向に限定するものではなく、全体座標系における各軸方向や要素における面内方向を任意に指定し、該指定した方向の応力の値に基づいて応力変更領域を選定するものであってもよい。
【００６３】
なお、上記の説明において、プレス成形品モデルは平面要素でモデル化されたものを対象としていたが、本発明は平面要素に限るものではなく、立体要素でモデル化されたプレス成形品モデルを対象としてもよい。
【符号の説明】
【００６４】
１  プレス成形品モデル
３  天板部
５  縦壁部
７  フランジ部
９  パンチ肩Ｒ部
【図１】
2020011277.tif 000003
【図２】
2020011277.tif 000004
【図３】
2020011277.tif 000005
【図４】
2020011277.tif 000006
【図５】
2020011277.tif 000007
【図６】
2020011277.tif 000008
【図７】
2020011277.tif 000009
【図８】
2020011277.tif 000010
【図９】
2020011277.tif 000011
【図１０】
2020011277.tif 000012