(57)【特許請求の範囲】
【請求項１】
金属板材をプレス成形する成形型に配設されるドロービードのモデリングを省略し該ドロービードが設けられ得るドロービード位置に該ドロービードのモデリングの代替となる等価ドロービードを設定して該金属板材がプレス成形される様子をシミュレーションするプレス成形解析システムであって、
前記等価ドロービードを構成するドロービードパラメータは、前記ドロービードの入口側における該金属板材の状態を指標する入口側指標値を考慮して設定され、
該入口側指標値は、該金属板材に作用し得る断面力である入口側断面力、該金属板材の板厚である入口側板厚、該金属板材に生じ得る相当塑性ひずみである入口側相当塑性ひずみ、および該金属板材と前記成形型との間に生じ得る間隙である入口側板隙量を含むことを特徴とするプレス成形解析システム。
【請求項２】
前記金属板材は、前記ドロービードを通過する際に伸びフランジ変形または縮みフランジ変形を伴う請求項１に記載のプレス成形解析システム。
【請求項３】
前記ドロービードパラメータと前記入口側指標値の関係は、関数またはテーブルにより対応付けられる請求項１または２に記載のプレス成形解析システム。
【請求項４】
前記ドロービードパラメータは、前記金属板材に作用するバウシンガ効果、ひずみ速度依存性またはクーロン摩擦係数の面圧依存性の少なくとも一つ以上を考慮して算出される請求項１または３に記載のプレス成形解析システム。
【請求項５】
前記ドロービードパラメータは、前記ドロービードの出口側で前記金属板材に作用する引抜力、該ドロービード位置で該金属板材に作用する押付力、該ドロービードを通過した該金属板材の板厚減少量、または該ドロービードを通過した該金属板材の相当塑性ひずみ増加量の少なくとも一つ以上である請求項１〜４のいずれかに記載のプレス成形解析システム。
【請求項６】
請求項１〜５のいずれかに記載のプレス成形解析システムに用いられることを特徴とするプレス成形解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、等価ドロービードを設定して金属板材のプレス成形をシミュレーション可能なプレス成形解析システムおよびそのプログラムに関する。
【背景技術】
【０００２】
金属板をプレス成形した部材（プレス品）が多用されている。このうち、比較的大型で複雑な曲面形状をもつ自動車ボディ等のプレス品は、金型の成形キャビティの外周囲（ダイフェース上）に、突起状のドロービード（絞りビード）を設けてプレス成形される。このドロービードにより、成形キャビティへの材料の流入量等の成形条件が調整される。特に、割れ、しわ、表面ひずみ、スプリングバック（そり、ねじれ）等の成形不良が発生し易い複雑で大型なプレス品ほど、ドロービードの設定が重要となる。
【０００３】
このドロービードを含むプレス成形用金型（成形型）は非常に高価である。このため、通常、実際に金型を製作する前に、プレス成形のシミュレーションが十分になされ、ドロービードを含む最適な金型形状が検討される。このシミュレーションの際、実際の三次元形状に沿ってドロービードをモデリングし、金属板材がそのドロービードを通過する過程を十分に再現されると好ましい。しかし、そのような方法では、ドロービードに細かい有限要素メッシュを形成する必要があり、プレス品が大型化するほど、シミュレーションに長時間を要して現実的ではない。そこで、金型モデルからドロービードのモデリングを省略し、金属板材がドロービードを通過した場合と同様な境界条件を与える等価ドロービードを代替的に設定したシミュレーションがなされている。
【０００４】
この等価ドロービードの設定により、プレス成形のシミュレーション結果は大きく影響を受ける。このシミュレーション精度を高めるために、等価ドロービードの設定方法に関する提案がいくつかなされている。例えば、それに関連する記載が下記の特許文献にある。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−２９０９６１号公報
【特許文献２】特開２０００−３０１２６３号公報
【特許文献３】特開２００５−１９３２９２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１は、ドロービード通過後における材料のひずみとドロービード引抜き抵抗の関係を示す第１設計線図と、そのドロービード引抜き抵抗とドロービードの等価曲げ半径の関係を示す第２設計線図とを用いたドロービードの設計方法について記載している。
【０００７】
特許文献２は、ダイフェースとプレス方向のなす角度に応じて、ドロービードの形成箇所における材料流入抵抗（引抜力またはビード力）を補正するプレス成形解析方法について記載している。
【０００８】
特許文献３は、金型の押付力と板材の引抜力の対応関係を予め実験的に求めておき、その押付力にプレス成形時のブランクホールドフォースを対応させると共にその引抜力を等価ドロービードのドロービード張力に対応させるプレス成形シミュレーション方法について記載している。
【０００９】
しかし、これら特許文献に記載の方法はいずれも、材料の流動方向に作用する力のみを考慮してドロービード部分における材料流入抵抗（引抜力）を求めているに過ぎない。このような方法では、角部のようにプレス成形時に伸びフランジ変形や縮みフランジ変形を生じて材料の幅方向にも力が作用するような場合や、ドロービードを複数平行に配置する場合に、材料流入抵抗を的確に表すことは困難である。
【００１０】
本発明はこのような事情に鑑みて為されたものであり、材料の幅方向にも力が作用するような場合や、ドロービードを連続で通過する場合でも、高精度なシミュレーション結果を得ることができるプレス成形解析システムおよびそのプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明者はこの課題を解決すべく鋭意研究し、試行錯誤を重ねた結果、ドロービードの入口側における板材の状況を考慮して、ドロービードを通過した板材に作用する引抜力等を定めることを思いついた。この成果を発展させることにより、以降に述べる本発明を完成するに至った。
【００１２】
《プレス成形解析システム》
（１）本発明のプレス成形解析システムは、金属板材をプレス成形する成形型に配設されるドロービードのモデリングを省略し該ドロービードが設けられ得るドロービード位置に該ドロービードのモデリングの代替となる等価ドロービードを設定して該金属板材がプレス成形される様子をシミュレーションするプレス成形解析システムであって、前記等価ドロービードを構成するドロービードパラメータは、前記ドロービードの入口側における該金属板材の状態を指標する入口側指標値を考慮して設定され、該入口側指標値は、該金属板材に作用し得る断面力である入口側断面力、該金属板材の板厚である入口側板厚、該金属板材に生じ得る相当塑性ひずみである入口側相当塑性ひずみ、および該金属板材と前記成形型との間に生じ得る間隙である入口側板隙量を含むことを特徴とする。
【００１３】
（２）本発明のプレス成形解析システムでは、等価ドロービードのドロービードパラメータ（適宜、単に「パラメータ」ともいう。）が、ドロービードの入口側における金属板材の状態を指標する入口側指標値を考慮して決定される。このため本発明によれば、従来のシステムとは異なり、伸びフランジ変形や縮みフランジ変形を伴う状態でも、金属板材がドロービードを通過するときの状態が正確に再現され、プレス成形性（材料流入量、板厚減少量等押付力）の予測精度が向上し、ひいては、より高精度なプレス成形シミュレーション結果が得られる。
【００１４】
《その他》
（１）本発明は、上述したプレス成形解析システムとしてのみならず、プレス成形解析プログラム、プレス成形解析方法、さらにはドロービードパラメータの決定方法、ドロービードパラメータの解析方法やそのプログラム等としても把握できる。
【００１５】
（２）特に断らない限り本明細書でいう「ｘ〜ｙ」は下限値ｘおよび上限値ｙを含む。本明細書に記載した種々の数値または数値範囲に含まれる任意の数値を新たな下限値または上限値として「ａ〜ｂ」のような範囲を新設し得る。
【図面の簡単な説明】
【００１６】
【図１Ａ】ドロービード解析に用いたモデリングを示す説明図である。
【図１Ｂ】そのモデリングで、成形型に設けたドロービードでワークを挟み込んだ状態を示す説明図である。
【図２】ドロービード解析して得られた後方張力と引抜力の関係を示すグラフである。
【図３】その解析で得られた入口側相当塑性ひずみと引抜力の関係を示すグラフである。
【図４】その解析で得られた入口側板厚と引抜力の関係を示すグラフである。
【図５】その解析で得られた入口側板隙量と引抜力の関係を示すグラフである。
【図６】実プレス品の一例を示す立体図である。
【図７】プレス成形解析して得られたストロークと引抜力の関係を示すグラフである。
【図８】プレス成形解析して得られたストロークと押付力の関係を示すグラフである。
【発明を実施するための形態】
【００１７】
本明細書で説明する内容は、本発明のプレス成形解析システムのみならず、それに関するプログラム等にも該当し得る。上述した本発明の構成要素に、本明細書中から任意に選択した一つまたは二つ以上の構成要素を付加し得る。いずれの実施形態が最良であるか否かは、対象、要求性能等によって異なる。
【００１８】
《入口側指標値》
ドロービードの入口側における金属板材の状態を指標する入口側指標値は、例えば、金属板材に作用し得る断面力である入口側断面力、金属板材の板厚である入口側板厚、金属板材に生じ得る相当塑性ひずみである入口側相当塑性ひずみ、金属板材と成形型との間に生じ得る間隙である入口側板隙量などがある。本発明に係る入口側指標値は、これらの少なくとも一つ以上であると好ましい。
【００１９】
《ドロービードパラメータ》
プレス成形シミュレーションを行う場合に設定される等価ドロービードを構成するドロービードパラメータは、例えば、ドロービードの出口側で金属板材に作用する引抜力、ドロービード位置で金属板材に作用する押付力、ドロービードを通過した金属板材の板厚減少量、ドロービードを通過した金属板材の相当塑性ひずみ増加量などがある。本発明のドロービードパラメータは、これらの少なくとも一つ以上であると好ましい。
【００２０】
《ドロービード解析》
入口側指標値はドロービードパラメータに影響を与える。両者の関係は市販されている解析ソフト等を用いて、ドロービード部分に特化して解析することが可能である。この際、金属板材に作用するバウシンガ効果、ひずみ速度依存性またはクーロン摩擦係数の面圧依存性等の少なくとも一つ以上を考慮して解析されると好ましい。なお、このドロービード部分の解析を行う場合、有限要素メッシュのサイズは、金属板材の板厚（ｔ）の（０．５〜１．０）ｔ程度で行うと好ましい。
【００２１】
こうして得られた入口側指標値とドロービードパラメータの対応関係は、関数またはテーブル（表）により対応付けられると好ましい。これにより入口側指標値の変化に適切に対応したドロービードパラメータを用いることが可能となる。
【００２２】
両者間の関係を示す関数は、例えば、一次関数または多次関数で表される。また、その関数は、一変数関数で表されても、多変数関数で表されてもよい。もっとも、一変数の一次関数または二次関数で表すと、効率的な処理が可能となり、好ましい。
《成形型》
本発明に係る成形型は、成形キャビティの形態を問わない。もっとも、ドロービードを通過する金属板材に、伸びフランジ変形や縮みフランジ変形による幅方向の力が作用するような複雑なプレス成形をシミュレーションする場合や、ドロービードを複数連続で通過させるような場合に本発明は特に有効である。
【００２３】
また本発明に係るドロービードも、その形態を問わない。ドロービードは直線状でも曲線状でもよい。またドロービードの断面形状も、角ビード、丸ビード（シングルビード、ダブルビート）、ステップビード等でもよい。いずれの場合でも、本発明によれば高精度なプレス成形シミュレーション結果が得られる。
【００２４】
《その他》
本発明によりプレス成形解析の対象となるプレス品として、例えば、自動車用部材（ドアパネル、ボンネット、フェンダー、サイドパネル等）がある。
【００２５】
金属板材は、解析対象となるプレス品に応じて適宜選択されるが、例えば、鉄系板材（鋼板等）、アンモニア系板材、チタン系板材、マグネシウム系板材などがある。なお「〜系板材」は純金属板材と合金板材を含む意味である。
【実施例】
【００２６】
実施例を挙げて本発明をより具体的に説明する。
【００２７】
《ドロービード解析》
図１Ａに示すように、ドロービード部分のみに着目したモデリングを行い、ワークＷ０（金属板材）に関する入口側指標値とドロービードパラメータの関係を、ＦＥＭ解析した。用いた解析ソフトはＪＳＴＡＭＰ／ＮＶ（株式会社ＪＳＯＬ製）およびＪＶＩＳＩＯＮ（株式会社ＪＳＯＬ製）である。
【００２８】
具体的には、鋼板を想定したプレス成形前のワークＷ０を、上型Ｄ１と下型Ｄ２（両者を併せて「成形型Ｄ」という。）にそれぞれ設けた凸型ドロービードｂ１と凹型ドロービードｂ２（両者を併せて「ドロービードｂ」という。）で挟み込む場合を仮定した。
【００２９】
そして図１Ａに示す状態から、凸型ドロービードｂ１と凹型ドロービードｂ２の間に所定の押付力Ｆｐを作用させて、両型間にワークＷ０を挟み込んだ様子を図１Ｂに示した。この状態で、ドロービードｂの入口側のワークＷ１を、出口側へ所定のストロークで引き抜いた場合に、ドロービードｂの出口側のワークＷ２に作用するドロービードパラメータ（引抜力Ｆｄ等）を求めた。
【００３０】
この際、入口側のワークＷ１には、後方張力Ｆｂを設定し、入口側断面力および入口側相当塑性ひずみを考慮した。またワークＷ１の入口側板厚とワークＷ１と成形型Ｄのクリアランスである入口側板隙量も考慮した。さらに、ワークＷ１、Ｗ２に作用するバウシンガ効果、ひずみ速度依存性、クーロン摩擦係数の面圧依存性も考慮した。
【００３１】
このような前提条件の下で、ワークＷ１に作用する入口側指標値を種々変更し、それぞれに対応するドロービードパラメータを算出した。得られた結果に基づき、入口側指標値とドロービードパラメータの相関を多変数関数で近似した。このような解析結果の一例として、ドロービードパラメータの一つである引抜力Ｆｄと各入口側指標値との相関関係を図２〜図５にそれぞれ示した。
【００３２】
図２に示した後方張力（入口側断面力）と引抜力（ドロービードパラメータ）の相関は、入口側相当塑性ひずみ＝０、入口側板厚＝初期板厚（ワークＷ０の板厚）、入口側板隙量＝０．０５として求めた。こうして得た４つの算出結果に基づいて、図２に示した二次関数を得た。
【００３３】
図３に示した入口側相当塑性ひずみと引抜力の相関は、後方張力＝０、入口側板厚＝初期板厚、入口側板隙量＝０．０５として求めた。こうして得た７つの算出結果に基づいて図３に示した二次関数を得た。
【００３４】
図４に示した入口側板厚と引抜力の相関は、後方張力＝０、入口側相当塑性ひずみ＝０、入口側板隙量＝０．０５として求めた。こうして得た４つの算出結果に基づいて図４に示した一次関数を得た。
【００３５】
図５に示した入口側板隙量と引抜力の相関は、後方張力＝０、入口側相当塑性ひずみ＝０、入口側板厚＝初期板厚として求めた。こうして得た３つの算出結果に基づいて図４に示した一次関数を得た。
【００３６】
《プレス成形解析》
上述した解析により得られたドロービードパラメータを用いて、図６に示すような実プレス品をプレス成形する際のシミュレーションを行った。この際用いたソフト（プレス成形解析システム）は、ＪＳＴＡＭＰ／ＮＶ（株式会社ＪＳＯＬ製）である。
【００３７】
ここでは、従来のシステムでは高精度な解析が困難であったドロービードがカーブした部分（図６のカーブ部）を想定し、その模擬形状（曲率半径１５０）をモデリングして、プレス成形シミュレーションを行った。その結果の一例を図７および図８示した。
【００３８】
図７は、ドロービードを通過したワークの変位量（ストローク）と引抜力の関係を示す。図８は、そのストロークとドロービードにおける押付力の関係を示す。両図中に示した比較例は、ドロービードの入口側におけるワーク状態（入口側指標値）を考慮していない従来のドロービードパラメータ（ストレート部のドロービードパラメータ）を用いて算出した結果である。また参照例は、等価ドロービードを設定せずに、ドロービードの三次元形状をモデリングしてＦＥＭ解析した結果である。
【００３９】
《評価》
図７および図８から明らかなように、本発明に係る実施例のシミュレーション結果は、参照例とほぼ同様な結果となり、従来のシミュレーション結果（比較例）に対して、著しく高精度化されたことがわかる。
【００４０】
例えば、図８に示すように、従来はストロークに伴って押付力が変化していなかったが、実施例ではストロークに応じて押付力が変化しており、参照例と同様な結果が忠実に再現されていることがわかる。
【００４１】
従って、上述したような方法で求めたドロービードパラメータを用いれば、ドロービードの形態がストレート状でない場合でも、等価ドロービードを用いたシミュレーションにより、高精度なプレス成形の予測が可能となる。特に、図６に示したＡ部のように、ドロービード形状自体はストレートでも、材料端でカーブが現れており、流入材料がドロービードの位置によって不均一となるような部位を含む場合にも、本実施例によればプレス成形を高精度に予測できる。
【符号の説明】
【００４２】
ｂ１、ｂ２  ドロービード
Ｄ１        上型（成形型）
Ｄ２        下型（成形型）
Ｗ０        ワーク（プレス成形前）
Ｗ１        ワーク（ドロービードの入口側）
Ｗ２        ワーク（ドロービードの出口側）
【図１Ａ】
0005797594.tif 000002
【図１Ｂ】
0005797594.tif 000003
【図２】
0005797594.tif 000004
【図３】
0005797594.tif 000005
【図４】
0005797594.tif 000006
【図５】
0005797594.tif 000007
【図６】
0005797594.tif 000008
【図７】
0005797594.tif 000009
【図８】
0005797594.tif 000010