(57)【要約】【課題】ＣＡＥ解析と実際の成形品のスプリングバック量乖離要因部位を特定するスプリングバック量乖離要因部位特定方法および装置を提供する。【解決手段】本発明は、プレス成形品の駆動応力分布を取得する成形品駆動応力分布取得ステップと、スプリングバック解析の駆動応力分布を取得する解析駆動応力分布取得ステップと、解析駆動応力分布と成形品駆動応力分布の応力差分分布を下死点の成形品形状に設定する応力差分分布設定ステップと、該応力差分分布によるスプリングバック量を取得する応力差分スプリングバック量取得ステップと、応力差分分布の一部の領域の値を変更してスプリングバック量を取得する変更応力差分スプリングバック量取得ステップと、取得したスプリングバック量を比較して、スプリングバック量に乖離が生ずる要因となる部位を特定するスプリングバック量乖離要因部位特定ステップとを備えていることを特徴とする。【選択図】  図１
2020179409.tif 000002
【特許請求の範囲】
【請求項１】
実パネルをプレス成形したプレス成型品に生ずるスプリングバック量と、前記プレス成形品と同形状の解析モデルについてスプリングバック解析を行った際のスプリングバック量に乖離が生ずる場合において、該乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定方法であって、
前記プレス成形品の離型後における表面形状を測定して取得した三次元形状測定データからプレス成形品モデルを作成し、該プレス成形品モデルを金型モデルによって下死点まで挟み込んだ状態の力学的解析を行い、成形下死点における応力分布を前記プレス成形品のスプリングバックに寄与した駆動応力分布として取得する成形品駆動応力分布取得ステップと、
前記スプリングバック解析における下死点応力分布及び離型後の残留応力分布を取得し、該下死点応力分布と離型後の残留応力分布の差分をスプリングバック解析における駆動応力分布として取得する解析駆動応力分布取得ステップと、
該解析駆動応力分布取得ステップで取得した解析駆動応力分布と前記成形品駆動応力分布取得ステップで取得した成形品駆動応力分布の差分から応力差分分布を取得して、前記スプリングバック解析における下死点の成形品形状に前記応力差分分布を設定する応力差分分布設定ステップと、
該設定した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する応力差分スプリングバック量取得ステップと、
前記応力差分分布設定ステップで設定した前記応力差分分布のうち、一部の領域の応力差分の値を変更して、該変更した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する変更応力差分スプリングバック量取得ステップと、
該変更応力差分スプリングバック量取得ステップ及び前記応力差分スプリングバック量取得ステップで取得したスプリングバック量を比較して、前記プレス成形品と前記スプリングバック解析のスプリングバック量に乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定ステップと、を備えていることを特徴とするスプリングバック量乖離要因部位特定方法。
【請求項２】
前記変更応力差分スプリングバック量取得ステップにおける応力差分の値の変更は、前記応力差分分布に対して少なくとも１方向の成分を、消去する、定数倍する、定数を加算する、定数乗する、被加工材の板厚方向の平均値に置き換える、被加工材の板厚方向の中央値に置き換える、のいずれかであることを特徴とする請求項１記載のスプリングバック量乖離要因部位特定方法。
【請求項３】
実パネルをプレス成形したプレス成型品に生ずるスプリングバック量と、前記プレス成形品と同形状の解析モデルについてスプリングバック解析を行った際のスプリングバック量に乖離が生ずる場合において、該乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定装置であって、
前記プレス成形品の離型後における表面形状を測定して取得した三次元形状測定データからプレス成形品モデルを作成し、該プレス成形品モデルを金型モデルによって下死点まで挟み込んだ状態の力学的解析を行い、成形下死点における応力分布を前記プレス成形品のスプリングバックに寄与した駆動応力分布として取得する成形品駆動応力分布取得手段と、
前記スプリングバック解析における下死点応力分布及び離型後の残留応力分布を取得し、該下死点応力分布と離型後の残留応力分布の差分をスプリングバック解析における駆動応力分布として取得する解析駆動応力分布取得手段と、
該解析駆動応力分布取得ステップで取得した解析駆動応力分布と前記成形品駆動応力分布取得ステップで取得した成形品駆動応力分布の差分から応力差分分布を取得して、前記スプリングバック解析における下死点の成形品形状に前記応力差分分布を設定する応力差分分布設定手段と、
該設定した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する応力差分スプリングバック量取得手段と、
前記応力差分分布設定ステップで設定した前記応力差分分布のうち、一部の領域の応力差分の値を変更して、該変更した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する変更応力差分スプリングバック量取得手段と、
該変更応力差分スプリングバック量取得ステップ及び前記応力差分スプリングバック量取得ステップで取得したスプリングバック量を比較して、前記プレス成形品と前記スプリングバック解析のスプリングバック量に乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定手段と、を備えていることを特徴とするスプリングバック量乖離要因部位特定装置。
【請求項４】
前記変更応力差分スプリングバック量取得手段における応力差分の値の変更は、前記応力差分分布に対して少なくとも１方向の成分を、消去する、定数倍する、定数を加算する、定数乗する、被加工材の板厚方向の平均値に置き換える、被加工材の板厚方向の中央値に置き換える、のいずれかであることを特徴とする請求項３記載のスプリングバック量乖離要因部位特定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、実パネルのプレス成形品のスプリングバック量とＣＡＥ解析のスプリングバック量に乖離が生じる要因となる前記プレス成形品における部位を特定するスプリングバック量乖離要因部位特定方法および装置に関する。
【背景技術】
【０００２】
金属板のプレス成型における成形品に対しては、高い形状精度が求められている。要求される形状精度を満足させるために、プレス成形後に金型から取り出した成形品が弾性変形することによって生じるスプリングバックを低減させることが重要である。
成形品の下死点における内部応力がスプリングバックの挙動に影響を及ぼすため、成形品のどの部位における応力がスプリングバックに対してどのような影響を与えるのかを把握することはスプリングバック対策を講じるのに有効である。
その手法として、特許文献１に開示されるプレス成形解析方法では、有限要素法を用いた解析によって解析を行っている。
【０００３】
特許文献１のプレス成形解析方法によれば、「プレス成形品である成形対象物のある領域についての残留応力分布を変更する前後において、スプリングバックに関するある定義された量がどのように変化するかを算出するので、この結果に基づいて、離型前の成形対象物のある領域の残留応力のスプリングバックへの影響を予測することができる」（発明の効果参照）としている。
特許文献１に開示されたような手法を用いることで、実際の金型を作成する前にスプリングバック対策を検討することができ、形状精度を確保するための金型調整作業を大幅に低減することができる。
【０００４】
その他のＣＡＥ解析を用いたスプリングバックの要因分析方法として、特許文献２では離型前の残留応力と離型後の残留応力からスプリングバック（ＳＢ）有効応力を算出し、該ＳＢ有効応力を用いて要因分析を行うことで、より適切な評価を行う方法が開示されている。
【０００５】
上述した方法はＣＡＥ解析におけるデータ設定から算出した応力状態よりスプリングバック解析を行っている。これに対し、特許文献３では、実際にプレス成形した成形品の表面形状を測定して作成した測定三次元形状をＣＡＥ解析に取り込み、該測定三次元形状を金型モデルによって成形下死点状態まで挟み込んだ状態の力学的解析を行って応力分布状態を取得し、該応力分布状態を用いて要因分析を行うことで、より正確な評価を行う方法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００７−２２９７２４号公報
【特許文献２】特開２０１２−２０６１５８号公報
【特許文献３】特開２０１３−７１１２０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、ＣＡＥ解析に基づいたスプリングバック対策を講じた金型を作製してプレス成形を行っても、実際のプレス成形品に講じた対策で期待された形状とは異なる形状となることがある。
これは、プレス成形に使用した金型の形状や種々の成形条件によって、金属板に対してＣＡＥ解析で想定した成形荷重がかからなかった場合や、ＣＡＥ解析が離型前の応力状態を正確に再現することができなかった場合などに起きうるもので、そのような場合にはプレス成形に用いる金型や成形条件を調整したり、ＣＡＥ解析上の設定や金型形状を見直したりする必要がある。
【０００８】
そのような調整作業においては、ＣＡＥ解析によるスプリングバック解析結果と実際の成形品のスプリングバック量の乖離要因が、成形品のどの部位で発生しているのかを特定することが求められるが、特許文献１乃至３に開示されている方法はスプリングバックそのものの発生要因となる部位を特定するものであり、ＣＡＥ解析と実スプリングバック量の乖離要因となる部位を特定するものではなかった。
【０００９】
また、ＣＡＥ解析を用いたスプリングバック対策の有用性を担保するためには、スプリングバック解析の再現性を向上させて、実際の成形品のスプリングバックと乖離しないようにすることが求められる。そのためにも、成形品において乖離要因となる部位を特定することが必要である。
【００１０】
本発明は、上記のような課題を解決するためになされたものであり、実パネルのプレス成形品のスプリングバック量とＣＡＥ解析のスプリングバック量に乖離が生じる要因となる前記プレス成形品における部位を特定するスプリングバック量乖離要因部位特定方法および装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
（１）本発明に係るスプリングバック量乖離要因部位特定方法は、実パネルをプレス成形したプレス成型品に生ずるスプリングバック量と、前記プレス成形品と同形状の解析モデルについてスプリングバック解析を行った際のスプリングバック量に乖離が生ずる場合において、該乖離が生ずる要因となる成形品形状における部位を特定するものであって、前記プレス成形品の離型後における表面形状を測定して取得した三次元形状測定データからプレス成形品モデルを作成し、該プレス成形品モデルを金型モデルによって下死点まで挟み込んだ状態の力学的解析を行い、成形下死点における応力分布を前記プレス成形品のスプリングバックに寄与した駆動応力分布として取得する成形品駆動応力分布取得ステップと、前記スプリングバック解析における下死点応力分布及び離型後の残留応力分布を取得し、該下死点応力分布と離型後の残留応力分布の差分をスプリングバック解析における駆動応力分布として取得する解析駆動応力分布取得ステップと、該解析駆動応力分布取得ステップで取得した解析駆動応力分布と前記成形品駆動応力分布取得ステップで取得した成形品駆動応力分布の差分から応力差分分布を取得して、前記スプリングバック解析における下死点の成形品形状に前記応力差分分布を設定する応力差分分布設定ステップと、該設定した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する応力差分スプリングバック量取得ステップと、前記応力差分分布設定ステップで設定した前記応力差分分布のうち、一部の領域の応力差分の値を変更して、該変更した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する変更応力差分スプリングバック量取得ステップと、該変更応力差分スプリングバック量取得ステップ及び前記応力差分スプリングバック量取得ステップで取得したスプリングバック量を比較して、前記プレス成形品と前記スプリングバック解析のスプリングバック量に乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定ステップと、を備えていることを特徴とするものである。
【００１２】
（２）また、上記（１）に記載のものにおいて、前記変更応力差分スプリングバック量取得ステップにおける応力差分の値の変更は、前記応力差分分布に対して少なくとも１方向の成分を、消去する、定数倍する、定数を加算する、定数乗する、被加工材の板厚方向の平均値に置き換える、被加工材の板厚方向の中央値に置き換える、のいずれかであることを特徴とするものである。
【００１３】
（３）本発明に係るスプリングバック量乖離要因部位特定装置は、実パネルをプレス成形したプレス成型品に生ずるスプリングバック量と、前記プレス成形品と同形状の解析モデルについてスプリングバック解析を行った際のスプリングバック量に乖離が生ずる場合において、該乖離が生ずる要因となる成形品形状における部位を特定するものであって、前記プレス成形品の離型後における表面形状を測定して取得した三次元形状測定データからプレス成形品モデルを作成し、該プレス成形品モデルを金型モデルによって下死点まで挟み込んだ状態の力学的解析を行い、成形下死点における応力分布を前記プレス成形品のスプリングバックに寄与した駆動応力分布として取得する成形品駆動応力分布取得手段と、前記スプリングバック解析における下死点応力分布及び離型後の残留応力分布を取得し、該下死点応力分布と離型後の残留応力分布の差分をスプリングバック解析における駆動応力分布として取得する解析駆動応力分布取得手段と、該解析駆動応力分布取得ステップで取得した解析駆動応力分布と前記成形品駆動応力分布取得ステップで取得した成形品駆動応力分布の差分から応力差分分布を取得して、前記スプリングバック解析における下死点の成形品形状に前記応力差分分布を設定する応力差分分布設定手段と、該設定した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する応力差分スプリングバック量取得手段と、前記応力差分分布設定ステップで設定した前記応力差分分布のうち、一部の領域の応力差分の値を変更して、該変更した応力差分分布に基づいてスプリングバック解析を行ってスプリングバック量を取得する変更応力差分スプリングバック量取得手段と、該変更応力差分スプリングバック量取得ステップ及び前記応力差分スプリングバック量取得ステップで取得したスプリングバック量を比較して、前記プレス成形品と前記スプリングバック解析のスプリングバック量に乖離が生ずる要因となる成形品形状における部位を特定するスプリングバック量乖離要因部位特定手段と、を備えていることを特徴とするものである。
【００１４】
（４）また、上記（３）に記載のものにおいて、前記変更応力差分スプリングバック量取得手段における応力差分の値の変更は、前記応力差分分布に対して少なくとも１方向の成分を、消去する、定数倍する、定数を加算する、定数乗する、被加工材の板厚方向の平均値に置き換える、被加工材の板厚方向の中央値に置き換える、のいずれかであることを特徴とするものである。
【発明の効果】
【００１５】
本発明によれば、ＣＡＥ解析によるスプリングバック解析と実際の成形品に生じたスプリングバック量の乖離の要因となる部位を特定することができるので、実際のプレス成形金型や成形条件を調整する作業の負荷を低減し、ＣＡＥ解析を用いたスプリングバック対策の有用性を向上させることができる。
【図面の簡単な説明】
【００１６】
【図１】本実施の形態１に係るスプリングバック量乖離要因部位特定方法の処理の流れを示すフロー図である。
【図２】本実施の形態１におけるプレス成形品を説明する図である。
【図３】本実施の形態１において、ＣＡＥ解析における駆動応力分布を示す図である。
【図４】本実施の形態１において、実パネルにおける駆動応力分布を示す図である。
【図５】本実施の形態１において、ＣＡＥ解析における駆動応力分布（ａ）と、実パネルにおける駆動応力分布（ｂ）と、ＣＡＥ解析と実パネルの駆動応力差分分布（ｃ）を示す図である。
【図６】ＣＡＥ解析と実パネルの駆動応力差分分布（ａ）と、該駆動応力差分分布に基づくスプリングバック解析により算出した変位（ｂ）を示す図である。
【図７】本実施の形態１におけるスプリングバック解析における固定点と該スプリングバック解析により求めた変位の解析結果を示す図である。
【図８】本実施の形態１においてスプリングバック量である首振り量とはね量の算出方法を説明する図である。
【図９】本実施の形態１において、成形品の領域分割を説明する図である。
【図１０】本実施の形態１において、応力差分を変更する領域（ａ）と、該領域の応力差分を変更してスプリングバック解析により算出した変位（ｂ）を示す図である（その１）。
【図１１】本実施の形態１において、応力差分を変更する領域（ａ）と、該領域の応力差分を変更してスプリングバック解析により算出した変位（ｂ）を示す図である（その２）。
【図１２】本実施の形態１において、応力差分を変更する領域（ａ）と、該領域の応力差分を変更してスプリングバック解析により算出した変位（ｂ）を示す図である（その３）。
【図１３】本実施の形態１において、応力差分を変更する領域（ａ）と、該領域の応力差分を変更してスプリングバック解析により算出した変位（ｂ）を示す図である（その４）。
【図１４】本実施の形態１において、応力差分を消去した領域と、該応力差分を消去したときのスプリングバックにより生じる首振り量の関係を示すグラフである。
【図１５】本実施の形態１において、応力差分を消去した領域と、該応力差分を消去したときのスプリングバックにより生じるはね量の関係を示すグラフである。
【図１６】ＣＡＥ解析の成形下死点における応力分布を実パネルの応力分布に置き換える領域と、該応力分布を変更したときのＣＡＥ解析と実パネルの首振りの乖離量を示すグラフである。
【図１７】ＣＡＥ解析の成形下死点における応力分布を実パネルの応力分布に置き換える領域と、該応力分布を変更したときのＣＡＥ解析と実パネルのはねの乖離量を示すグラフである。
【図１８】本実施の形態２に係るスプリングバック量乖離要因部位特定装置の構成を説明するためのブロック図である。
【発明を実施するための形態】
【００１７】
本発明の実施の形態に係るスプリングバック量乖離要因部位特定方法を説明する前に、まず本発明に至った経緯について説明する。
【００１８】
＜本発明に至った経緯＞
図２に一例として示すような成形品１のプレス成形においては、プレス成形前に行うスプリングバック解析（ＣＡＥ解析ともいう）によって算出されたスプリングバック量と実際にプレス成形された成形品（実パネルともいう）のスプリングバック量に差（乖離）が生じる場合があった。そして、このようなＣＡＥ解析と実パネルにスプリングバック量の乖離が生じる要因となる部位は、スプリングバックそのものが発生する要因となる部位とは異なる場合がある。そのため、スプリングバックが発生する要因となる部位に何らかの対策を施してスプリングバックを低減したとしても、ＣＡＥ解析と実パネルの間に生じるスプリングバック量の乖離を低減するには至らないという問題があった。
【００１９】
そこで発明者は、このような問題を解決するために鋭意検討した。その結果、スプリングバック量の乖離の要因が、ＣＡＥ解析におけるスプリングバックに寄与した応力と、実パネルにおけるスプリングバックに寄与した応力の違いに起因するのではないかと着想するに至った。
【００２０】
そこで、この着想の妥当性を検証するために、図２に示すハット断面形状の成形品１について、ＣＡＥ解析及び実パネルのそれぞれについてスプリングバックに寄与した応力（以降、駆動応力分布ともいう）を求め、ＣＡＥ解析における駆動応力分布と実パネルにおける駆動応力分布の差分（以降、応力差分分布ともいう）を算出し、該応力差分分布をＣＡＥ解析の成形下死点における応力分布に置き換えてスプリングバック解析を行った（駆動応力算出方法については後述する実施の形態で詳しく説明する）。
【００２１】
スプリングバック量として、長手方向端部の板幅方向への変位量である首振り量と、成形ストローク方向への変位量であるはね量（図８参照）を算出し、ＣＡＥ解析と実パネルのスプリングバック量の差（乖離量）と比較した結果を表１に示す。
【００２２】
【表１】
2020179409.tif 000003
【００２３】
表１に示すように、ＣＡＥ解析と実パネルのスプリングバック量の差と応力差分分布を設定したスプリングバック解析により得られたスプリングバック量とはほぼ同じ値になることが確認された。つまり、ＣＡＥ解析と実パネルの駆動応力分布の違いである応力差分が、スプリングバック量の乖離の原因であると考えてよいという知見を得た。
【００２４】
本発明に係るスプリングバック量乖離要因部位特定方法は、上記の経緯により完成されたものであり、その具体的な方法を以下の実施の形態にて説明する。
【００２５】
［実施の形態１］
本発明の実施の形態１に係るスプリングバック量乖離要因部位特定方法は、実パネルをプレス成形したプレス成形品とスプリングバック解析のスプリングバック量に乖離が生じる場合において、該スプリングバック量に乖離が生じる要因となる成形品形状における部位を特定するものであって、図１に示すように、成形品駆動応力分布取得ステップＳ１と、解析駆動応力分布取得ステップＳ３と、応力差分分布設定ステップＳ５と、応力差分スプリングバック量取得ステップＳ７と、変更応力差分スプリングバック量取得ステップＳ９と、スプリングバック量乖離要因部位特定ステップＳ１１とを備えたものである。
【００２６】
図２に示すようなハット断面形状の成形品１をプレス成形する場合を例とし、上記の各ステップを説明する。
本実施の形態におけるＣＡＥ解析では、図２に示すように、ダイ５とパンチ７からなる金型モデル３により被加工材（鋼板）であるブランクモデル９を挟むプレス成形解析を実施する。なお、プレス成形解析においては、図２に示すように位置決めピンにより成形過程においてブランクモデル９を固定するものとし、ブランクモデル９の要素サイズを約1mm、解析条件として、ブランクモデル９と金型モデル３との間の摩擦係数を0.15、成形下死点位置を上下金型のモデルのスキが1.45mmとなる位置とした。また、被加工材は、板厚1.4mmの980MPa級GA鋼板とした。
なお、本実施の形態における実パネルはＣＡＥ解析に設定した成形条件と同じ条件下でプレス成形を行って成形したものとする。
【００２７】
＜成形品駆動応力分布取得ステップ＞
成形品駆動応力分布取得ステップＳ１は、実パネルの駆動応力分布を取得するステップである。
具体的には、上述したようなＣＡＥ解析に設定した成形条件と同じ条件でプレス成形を行って実パネルを成形し、該実パネルの離型後における表面形状を測定して取得した三次元形状測定データからプレス成形品モデルを作成し、該プレス成形品モデルを図２に示した金型モデル３によって成形下死点まで挟み込んだ状態の力学的解析を行って図４に示したような応力分布を取得するものである。
上記力学的解析として弾性有限要素解析を行い、該弾性有限要素解析により得られた応力分布は実パネルのスプリングバックに寄与した応力、すなわち、実パネルの駆動応力に相当する。
【００２８】
ここで、実パネルの三次元形状の測定、プレス成形品モデルの作成および弾性有限要素解析の具体的な方法としては、例えば、特許文献３に記載されている方法を用いることができる。
【００２９】
＜解析駆動応力分布取得ステップ＞
解析駆動応力分布取得ステップＳ３は、ＣＡＥ解析（スプリングバック解析）における下死点応力分布及び離型後の残留応力をそれぞれ取得して、その差分からＣＡＥ解析における駆動応力を取得するステップである。
プレス成形品の下死点における応力がスプリングバックの挙動に影響を及ぼすものではあるが、金型から離型した後のプレス成形品にも応力が残留している。すなわち、離型前の応力のすべてがスプリングバックに寄与するわけではなく、離型後のスプリングバックが生じたプレス成形品に残留する応力はスプリングバックに寄与しなかったものと考えられる。
成形品駆動応力分布取得ステップＳ１で取得した駆動応力（スプリングバックに寄与した応力）との正確な比較分析を行うためにはＣＡＥ解析における駆動応力を算出する必要がある。
【００３０】
そこで、解析駆動応力分布取得ステップＳ３では、図２に示した金型モデル３及びブランクモデル９を用いてＣＡＥ解析を行い、離型前（下死点）における成形品の応力及び離型後（スプリングバック後）における成形品の残留応力を取得し、下死点の応力から離型後の残留応力を差し引くことで図３に示したようなＣＡＥ解析における駆動応力分布（以降、単に応力分布ともいう）を算出するようにした。
【００３１】
なお、解析駆動応力分布取得ステップＳ３は、コンピュータがＣＡＥ解析を行うものであり、ＣＡＥ解析には、例えば、有限要素法解析ソフトウェアを用いることができる。本実施の形態では、市販の有限要素法解析ソフトウェアであるLS−DYNA  Ver.971をコンピュータ上で実行することにより解析を行い、ソルバーには動的陽解法を適用した。
【００３２】
また、本発明は、図２に示すような金型モデル３を用いるものや、ハット断面形状の成形品１を成形対象とするものに限らず、成形対象に応じて金型モデルや成形品などを適宜設定することができる。
さらに、離型後の残留応力がほぼ無視できるような小さな値であった場合には、解析下死点の応力分布の全てがスプリングバックに寄与するとみなし、下死点における応力分布を解析駆動応力分布としてもよい。
【００３３】
＜応力差分分布設定ステップ＞
応力差分分布設定ステップＳ５は、図５に示すように、解析駆動応力分布取得ステップＳ３で算出した応力分布（図３、図５（ａ））と成形品駆動応力分布取得ステップＳ１で算出した応力分布（図４、図５（ｂ））との差分を応力差分分布として算出し（図５（ｃ））、算出した応力差分分布を、ＣＡＥ解析の下死点形状における応力分布に設定するステップである。
なお、図５に示す各応力分布は応力の大きさの違いが色の濃淡で示されているが、図５（ｃ）は、図５（ａ）及び図５（ｂ）よりも色表示のレンジ幅を小さく設定したものである。このようにすることで、応力差が大きい部分（黒または白に近い部分）と応力差が小さい部分（図中右側のグレースケール中間色に近い部分）をより区別しやすくしている。
【００３４】
＜応力差分スプリングバック量取得ステップ＞
応力差分スプリングバック量取得ステップＳ７は、応力差分分布設定ステップＳ５で設定した応力差分分布を用いてスプリングバック解析を行い、そこで生じるスプリングバック量を算出するステップである。
【００３５】
本実施の形態では、応力差分分布が設定された下死点形状の成形品１（図６（ａ））に対してスプリングバック解析を行い、該スプリングバック解析によりスプリングバック後の変位（図６（ｂ））を算出する。スプリングバック解析においては、図７に示すように、成形品の一端側に設けた３箇所を固定点で固定し、スプリングバックによる変位を算出した。
【００３６】
次に、図８に示すように、スプリングバック解析により算出した変位に基づいて、首振り量（図８（ａ））およびはね量（図８（ｂ））を算出した。本実施の形態において、首振り量は、図８（ａ）に示すように、成形品１の他端側に設けた２箇所の評価点どちらかのスプリングバックによる図中矢印方向への移動量（矢印の向きを正）とし、はね量は、図８（ｂ）に示すように、２つの評価点の中点における成形ストローク方向の変位量（ダイ５から離れる向きを正）とした。
【００３７】
＜変更応力差分スプリングバック量取得ステップ＞
変更応力差分スプリングバック量取得ステップＳ９は、応力差分分布設定ステップＳ５で設定した応力差分分布のうち、成形品１における一部分のある部位の応力差分の値を変更し、変更した応力差分分布（変更応力差分分布）から、スプリングバック解析を行い、そこで生じるスプリングバック量を算出するステップである。
【００３８】
本実施の形態では、図９に示すように、成形品１を複数の領域に分割（長手方向にＡ〜Ｆの６分割、幅方向に１〜３の３分割）し、各領域における応力差分を消去、すなわち、応力差分の値をゼロにする。なお、応力差分を消去する領域は、必ずしも１領域である必要はなく、複数の領域の応力差分を消去するものでもよい。
ある領域の値を消去した後、スプリングバック解析を行って、スプリングバック量として図８に示した首振り量とはね量をそれぞれ算出する。
【００３９】
図１０に、応力差分分布設定ステップＳ５で設定した応力差分分布から領域１（図１０（ａ）参照）における応力差分を消去（応力差分の値をゼロ）し、スプリングバック解析を行って算出した変位（図１０（ｂ）参照）を示す。算出した変位に基づいてスプリングバック量（首振り量、はね量）を算出した。領域２、領域３についても同様にスプリングバック量を算出した（図示は省略）。
【００４０】
同様に、図１１に、応力差分分布設定ステップＳ５で設定した応力差分分布から領域Ａ（図１１（ａ）参照）における応力差分を消去（応力差分の値をゼロ）し、スプリングバック解析を行って算出した変位（図１１（ｂ）参照）を示す。算出した変位に基づいてスプリングバック量（首振り量、はね量）を算出した。領域Ｂ〜Ｆについても同様にスプリングバック量を算出した（図示は省略）。
【００４１】
また、図１２、図１３に、応力差分分布設定ステップＳ５で設定した応力差分分布から領域Ａ−１、領域Ａ−３（図１２（ａ）、図１３（ａ）参照）における応力差分を消去（応力差分の値をゼロ）し、スプリングバック解析を行って算出した変位（図１２（ｂ）、図１３（ｂ）参照）を示す。算出した変位に基づいてスプリングバック量（首振り量、はね量）を算出した。領域Ａ−１〜領域Ｆ−３についても同様にスプリングバック量を算出した（図示は省略）。
【００４２】
＜スプリングバック量乖離要因部位特定ステップ＞
スプリングバック量乖離要因部位特定ステップＳ１１は、応力差分スプリングバック量取得ステップＳ７及び変更応力差分スプリングバック量取得ステップＳ９で算出したスプリングバック量に基づいて、ＣＡＥ解析と実パネルでスプリングバック量に乖離が生じる要因となる部位を特定するステップである。
【００４３】
図１４に、応力差分スプリングバック量取得ステップＳ７及び変更応力差分スプリングバック量取得ステップＳ９で算出した首振り量の結果を示す。図１４において、baseは、応力消去を行っていない応力差分分布による首振り量であり、１〜３、Ａ〜Ｆ、Ａ−１〜Ｆ−３は、各領域の応力差分を消去した変更応力差分分布による首振り量である。
【００４４】
応力差分を変更しない場合の首振り量（base）は、-10.7mmである。これに対し、応力消去領域が領域１〜３の首振り量をみると、領域２のみ領域１、３に比べてbaseからの変化が大きいことがわかる。
同様に、応力消去領域が領域Ａ〜Ｆの首振り量をみると、領域Ｄ及び領域Ｅが他の領域に比べてbaseからの変化が大きいことがわかる。
【００４５】
次に、さらに細分化した領域である領域Ａ−１〜領域Ｆ−３を応力消去した場合の首振り量をみると、領域Ｄ−２及び領域Ｅ−２の首振り量が応力消去を行っていないbaseに比べて低減していることがわかる。このように、ある領域の応力差分を消去したことによって首振り量が低減するのは、当該領域が首振り量の乖離に大きく影響することを示している。
上述した領域２及び領域Ｄ、Ｅにおいても変化が大きかったことからもこれは妥当な結果であり、領域Ｄ−２及び領域Ｅ−２が、ＣＡＥ解析と実パネルの首振り量に乖離を生じさせる部位であると特定することができる。
【００４６】
なお、領域Ｄ−２及び領域Ｅ−２以外の領域における首振り量は、応力差分を変更しない場合の首振り量（base）と同程度であることから、これらの領域は、ＣＡＥ解析と実パネルの首振り量の乖離に対する影響が小さい部位であると判断される。
【００４７】
はね量に関しても同様に、算出した結果を図１５に示す。
応力差分を変更しない場合のはね量（base）は、4.2mmである。これに対し、応力消去領域が領域１〜３のはね量をみると、領域２（-8.3mm）のみ領域１（4.8mm）、領域３（7.5mm）に比べてbaseからの変化が大きいことがわかる。
同様に、応力消去領域が領域Ａ〜Ｆのはね量をみると、領域Ｄ（-2.9mm）及び領域Ｆ（-0.9mm）が他の領域に比べてbaseからの変化が大きいことがわかる。
【００４８】
次に、さらに細分化した領域である領域Ａ−１〜領域Ｆ−３を応力消去した場合のはね量をみると、領域Ｆ−２のはね量（0.1mm）が応力消去を行っていないbaseに比べて大きく低減しており、はねがほとんど生じていないことがわかる。首振り量と同様、ある領域の応力差分を消去したことによってはね量が低減するのは、当該領域がはね量の乖離に大きく影響することを示している。
領域２及び領域Ｆにおいてもはね量の変化が大きかったことからもこれは妥当な結果であり、領域Ｆ−２が、ＣＡＥ解析と実パネルのはね量に乖離を生じさせる部位であると特定することができる。
【００４９】
領域２及び領域Ｄもbaseからの変化が大きかったが、細分化した領域Ｄ−２の応力を消去した場合の結果をみるとはね量が-4.6mmとなっており、応力消去を行わなかったbaseのはね量4.2mmと比べると、はね方向が逆向きになり、かつ、はね量が増加している。
これは、領域Ｄ−２における応力差分を消去することで、当該領域がはね量をさらに乖離させやすい部位であることを示すものである。すなわち、領域Ｄ−２は応力差分があることによってはね量の乖離を抑制する部位であり、乖離を生じさせる部位ではないと判断される。
【００５０】
なお、本実施の形態では図１４、図１５に示すように、細分化した領域Ａ−１〜Ｆ−３のすべてにおいて応力消去を行った場合のスプリングバック量を算出したが、前述したように領域１〜３、領域Ａ〜Ｆの結果から、影響が大きいと想定される領域を絞ることができるため、その領域のみさらに細分化した場合の結果を求めるようにしてもよい。複雑な形状の成形品は分割領域を多くしてさらに細分化する必要があり、すべての細分化領域において結果を求めることは時間を要する。そのため、大きい領域分割から乖離に影響がありそうな部位を絞り込みながら分析を行うことで、乖離要因部位特定までの時間を短縮することができる。
【００５１】
次に、本実施の形態に係るスプリングバック量乖離要因部位特定方法により特定された部位が、ＣＡＥ解析と実パネルのスプリングバック量に乖離が生じる要因となる部位であると判断することの妥当性について説明する。
【００５２】
上述のとおり、図２に示す成形品１を成形するにあたって、成形品１における領域Ｄ−２及び領域Ｅ−２が、ＣＡＥ解析と実パネルの首振り量に乖離が生じる要因となる部位であると特定された。また、領域Ｆ−２が、はね量に乖離が生じる要因となる部位であると特定された。
そこで、ＣＡＥ解析によって得られた駆動応力分布（図３参照）における、前記特定された各部位の駆動応力を、これに対応する実パネルの各部位の駆動応力分布に置換して、スプリングバック解析を行い、スプリングバック量の乖離が低減するかどうかについて検証した。
【００５３】
まず、ＣＡＥ解析と実パネルの首振りの乖離量が低減するかどうかを検証した結果を図１６に示す。
図１６はＣＡＥ解析における首振り量から実パネルにおける首振り量を差し引いた結果（乖離量）をグラフにしたものであり、縦軸が乖離量を、横軸が応力置換領域をそれぞれ示している。横軸における「なし」は、応力置換領域が無い場合であり、この場合には、乖離量はＣＡＥ解析と実パネルの差である-11.6mmとなっている。また、「ＡＬＬ」は、ＣＡＥ解析の駆動応力分布の全ての領域を実パネルの駆動応力分布に置換した場合であり、この場合には、乖離量はほぼゼロとなる。
これらに対して、本実施の形態で特定された領域Ｄ−２または領域Ｅ−２のみに実パネルの応力分布を置き換えた場合は、図１６の「Ｄ−２」「Ｅ−２」に示されるように、いずれにおいても実パネルとの首振り乖離量が減少した。
さらに、「Ｄ−２／Ｅ−２」に示されるように、領域Ｄ−２、Ｅ−２の両方に実パネルの応力分布を置き換えることで、実パネルとの首振り乖離量をより低減できることが示された。
【００５４】
同様に、ＣＡＥ解析と実パネルのはねの乖離量が低減するかどうかを検証した結果を図１７に示す。はね量の乖離についても、「Ｆ−２」に示されるように、成形下死点における領域Ｆ−２の応力分布を実パネルの応力分布に置き換えることで、実パネルとのはね乖離量が減少した。
【００５５】
この結果は、本実施の形態に係るスプリングバック量乖離要因部位特定方法により特定された部位が、ＣＡＥ解析と実パネルのスプリングバック量乖離の要因となる部位として妥当であることを示すものである。
【００５６】
以上のように、本実施の形態に係る方法によれば、ＣＡＥ解析と実パネルのスプリングバック量乖離の要因となる部位を精度良く特定できることが示された。
このようにスプリングバック量乖離の要因となる部位を特定することで、ＣＡＥ解析と実パネルの離型後形状を近づけるために行っていた金型や成形条件などの調整作業を効率的に行うことができる。
また、特定した部位を中心にＣＡＥ解析上の設定や金型形状を見直すことで、スプリングバック解析の再現性を向上させることができ、ＣＡＥ解析を用いたスプリングバック対策の有用性を担保することができる。
【００５７】
なお、本実施の形態で行う応力差分分布の変更は、少なくともある一部の領域の応力差分を消去（すべての成分をゼロに）するものであったが、応力差分分布を変更する方法はこれに限るものではなく、応力差分分布に対して少なくとも１方向の成分を、消去する、定数倍する、定数を加算する、定数乗する、被加工材の板厚方向の平均値に置き換える、被加工材の板厚方向の中央値に置き換える、のいずれかにより行えばよい。
【００５８】
［実施の形態２］
実施の形態１で説明したスプリングバック量乖離要因部位特定方法は、予め設定されたプログラムをＰＣ（パーソナルコンピュータ）に実行させることで実現できる。そのような装置の一例を本実施の形態にて説明する。
本実施の形態に係るスプリングバック量乖離要因部位特定装置１１は、図１８に一例を示すような、表示装置１３と、入力装置１５と、主記憶装置１７と、補助記憶装置１９と、演算処理部２１を有している。演算処理部２１には、表示装置１３、入力装置１５、主記憶装置１７及び補助記憶装置１９が接続され、演算処理部２１の指令によって各機能を行う。
【００５９】
表示装置１３は実行結果の表示等に用いられ、液晶モニター等で構成される。入力装置１５はオペレータからの入力等に用いられ、キーボードやマウス等で構成される。主記憶装置１７は演算処理部２１で使用するデータの一時保存や演算等に用いられＲＡＭ等で構成される。補助記憶装置１９はデータの記憶等に用いられ、ハードディスク等で構成される。
補助記憶装置１９には少なくとも、三次元形状測定データ２３、金型モデル２５等のＣＡＥ解析に必要な各種データが記憶されている。
【００６０】
演算処理部２１はＰＣなどのＣＰＵ等によって構成されており、演算処理部２１で予め設定されたプログラムが実行されることで、成形品駆動応力分布取得手段２７と、解析駆動応力分布取得手段２９と、応力差分分布設定手段３１と、応力差分スプリングバック量取得手段３３と、変更応力差分スプリングバック量取得手段３５と、スプリングバック量乖離要因部位特定手段３７とが実現される。
【００６１】
成形品駆動応力分布取得手段２７は実施の形態１で説明した成形品駆動応力分布取得ステップＳ１と同様の処理を実現するものである。同様に、解析駆動応力分布取得手段２９は解析駆動応力分布取得ステップＳ３を、応力差分分布設定手段３１は応力差分分布設定ステップＳ５を、応力差分スプリングバック量取得手段３３は応力差分スプリングバック量取得ステップＳ７を、変更応力差分スプリングバック量取得手段３５は変更応力差分スプリングバック量取得ステップＳ９を、スプリングバック量乖離要因部位特定手段３７はスプリングバック量乖離要因部位特定ステップＳ１１をそれぞれ実現するものである。
【００６２】
上述したような本実施の形態によれば、実施の形態１と同様にＣＡＥ解析と実パネルのスプリングバック量乖離の要因となる部位を精度良く特定することができる。
よって、特定した部位に基づいてＣＡＥ解析と実パネルの離型後形状を近づけるための各種調整作業を効率的に行うことができ、さらにスプリングバック対策の有用性を担保することができる。
【００６３】
本発明は最終成形品（製品形状）に限らず、成形工程を分けて行うような成形品の場合には途中成形品にも適用可能である。また、上述した実施の形態では鋼板を被加工材とする例を説明したが、アルミにも適用可能である。
【符号の説明】
【００６４】
１  成形品
３  金型モデル
５  ダイ
７  パンチ
９  ブランクモデル
１１  スプリングバック量乖離要因部位特定装置
１３  表示装置
１５  入力装置
１７  主記憶装置
１９  補助記憶装置
２１  演算処理部
２３  三次元形状測定データ
２５  金型モデル
２７  成形品駆動応力分布取得手段
２９  解析駆動応力分布取得手段
３１  応力差分分布設定手段
３３  応力差分スプリングバック量取得手段
３５  変更応力差分スプリングバック量取得手段
３７  スプリングバック量乖離要因部位特定手段
【図１】
2020179409.tif 000004
【図２】
2020179409.tif 000005
【図３】
2020179409.tif 000006
【図４】
2020179409.tif 000007
【図５】
2020179409.tif 000008
【図６】
2020179409.tif 000009
【図７】
2020179409.tif 000010
【図８】
2020179409.tif 000011
【図９】
2020179409.tif 000012
【図１０】
2020179409.tif 000013
【図１１】
2020179409.tif 000014
【図１２】
2020179409.tif 000015
【図１３】
2020179409.tif 000016
【図１４】
2020179409.tif 000017
【図１５】
2020179409.tif 000018
【図１６】
2020179409.tif 000019
【図１７】
2020179409.tif 000020
【図１８】
2020179409.tif 000021