JP6852426B2 - 成形性評価方法、プログラム及び記録媒体 - Google Patents

成形性評価方法、プログラム及び記録媒体 Download PDF

Info

Publication number
JP6852426B2
JP6852426B2 JP2017018964A JP2017018964A JP6852426B2 JP 6852426 B2 JP6852426 B2 JP 6852426B2 JP 2017018964 A JP2017018964 A JP 2017018964A JP 2017018964 A JP2017018964 A JP 2017018964A JP 6852426 B2 JP6852426 B2 JP 6852426B2
Authority
JP
Japan
Prior art keywords
strain
fracture
hole expansion
hole
different
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017018964A
Other languages
English (en)
Other versions
JP2017140653A (ja
Inventor
繁 米村
繁 米村
淳 新田
淳 新田
吉田 亨
亨 吉田
聡 白神
聡 白神
隆 安富
隆 安富
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Steel Corp
Original Assignee
Nippon Steel Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Steel Corp filed Critical Nippon Steel Corp
Publication of JP2017140653A publication Critical patent/JP2017140653A/ja
Application granted granted Critical
Publication of JP6852426B2 publication Critical patent/JP6852426B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Shaping Metal By Deep-Drawing, Or The Like (AREA)

Description

本発明は、成形性評価方法、プログラム及び記録媒体に関するものである。
近年、衝突安全性と軽量化の要請から、自動車車体への高強度鋼板の適用が急速に進展しつつある。これら高強度鋼板は板厚を増加させることなく、衝突時の吸収エネルギーや強度を高めることができる。しかしながら、鋼板の高強度化に伴う延性の低下は、プレス成形時の破断の危険性を高めるため、有限要素法による材料の破断予測とその高精度化のニーズは高まってきている。
成形時の破断に対する余裕度は、一般に、板厚減少率や成形限界線図(FLD)を用いて判断される。FLDは破断限界を与える最大主歪みを最小主歪みごとに示した図であり、成形解析や衝突解析での破断評価に用いられている。実験によるFLDの測定方法は、一般に、予め金属板の表面にエッチング等によりサークル状或いは格子状の模様を描いておき、液圧成形や剛体工具での張出し成形で破断させた後に、サークルの変形量から破断限界歪みを測定する。破断限界線は、様々な面内歪み比について金属板を比例負荷し、それぞれの歪み比での破断限界歪みを主歪み軸上にプロットして線で結ぶことで得られる(図1)。
一方、FLDの理論予測としては、Hillの局部くびれモデルとSwiftの拡散くびれモデルの併用、Marciniak-Kuczynski法、Storen-Riceモデルなど、種々がある。材料の延性破壊は局部くびれにより変形が局所化した位置で発生する。この局部くびれが発生すると極めて短時間で破断に至るため、実用上、破断限界は局部くびれ発生限界と考えることが多く、破断限界予測は塑性不安定の枠組みで取り扱うことが多い。このようにして得た破断限界線と有限要素法による数値シミュレーションの結果から得られる各部位の歪み状態との位置関係を比較することで破断の危険性を評価し、変形過程の歪みがこの限界歪みに達したときに破断、若しくはその危険性が高いと判断する。
特願2014−137185 特開2011−140046号公報
実験や理論予測から得られるFLDは、一様な応力状態のもとで材料が分離するとき、若しくは局部くびれが生じるときを対象にしたものである(図2)。しかしながら、鋼板の端部から亀裂が生じる伸びフランジ成形では、フランジ端部から内側に向けて歪みが小さくなるため、材料端部は内側の拘束を受け、くびれの発生が抑制される(図3)。即ち、伸びフランジ端部が一様分布における破断条件を満足しても内側ではまだその条件に達していないので、内側の支持効果により全体としては塑性不安定状態とはなり得ず破断に至らない。この点が単軸引張りや張出しや深絞りのような一様応力場での局所くびれ発生と異なる点であり、伸びフランジ破断のようにフランジ端部から内側に向けて歪み勾配が存在する場合の不安定くびれの発生条件はまだ解明されていない。また、せん断加工時の鋼板端部に導入される微視的損傷の影響により破断機構は複雑であり、これと前述した歪み勾配の影響により従来のFLDによる破断予測では予測精度を確保することができない。
本発明は、上記の課題に鑑みてなされたものであり、高強度鋼板適用時の成形課題の一つである伸びフランジ破断を未然に回避し、高強度で軽量な部品のプレス成形を実現することができる成形性評価方法、プログラム及び記録媒体を提供することを目的とする。
上記の課題を解決するため、本発明者らは、鋭意検討の結果、以下に示す発明の諸様態に想到した。本発明の要旨は、次の通りである。
1.材料のプレス成形の可否を評価する方法であって、
穴広げ試験から得られる2つ以上の異なる穴広げ率を入力する第1のステップと、
前記2つ以上の異なる穴広げ率のデータから、前記2つ以上の異なる破断限界歪み及び前記2つ以上の異なる半径方向の歪み勾配を計算する第2のステップと、
前記2つ以上の破断限界歪みと前記2つ以上の歪み勾配との関係から破断クライテリアを計算する第3のステップと、
有限要素法による数値解析から得られる最大主歪み及び隣接する要素間における歪み勾配が前記破断クライテリアに達したときに破断したと評価する第4のステップと
を含み、
前記第1のステップ〜前記第4のステップを一連の工程として実行して得られた解析結果に基づいて、前記材料の成形性を評価することを特徴とする成形性評価方法。
2.前記第1のステップにおいて、前記2つ以上の異なる頂角の円錐穴広げ試験から得られる夫々の前記穴広げ率を入力することを特徴とする1.に記載の成形性評価方法。
3.前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円錐工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする1.に記載の成形性評価方法。
4.前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円筒工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする1.に記載の成形性評価方法。
5.前記第2のステップにおいて、前記穴広げ率をλ(%)、前記板端部の破断限界歪みを(εθr=0、前記歪み勾配をdεθ/dr、円錐工具又は円筒工具の半頂角をφ(deg.)、素板の穴直径をd0(mm)、変形後の位置を表す半径座標をr(mm)とすると、
前記材料の板端部で破断が生じたときの穴の直径dは、(1+λ/100)d0として得られ、
(εθr=0=ln(d/d0
Figure 0006852426
として計算することを特徴とする1.〜4.のいずれか1項に記載の成形性評価方法。
6.前記歪み勾配であるdεθ/drを計算する際に、有限要素法の解析に用いる要素サイズを基準長さdrとすることを特徴とする1.〜5.のいずれか1項に記載の成形性評価方法。
7.前記2つ以上の破断限界歪み及び歪み勾配から前記破断クライテリアεcr=f(dεθ/dr)を計算し、有限要素法による数値解析から、隣接する要素間での前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11を求め、前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11が前記破断クライテリアに達しているか否かの指標としてε11/εcrを計算し、その結果をコンター表示することを特徴とする1.〜6.のいずれか1項に記載の成形性評価方法。
8.材料のプレス成形の可否を評価するためのプログラムであって、
穴広げ試験から得られる2つ以上の異なる穴広げ率を入力する第1のステップと、
前記2つ以上の異なる穴広げ率のデータから、前記2つ以上の異なる破断限界歪み及び前記2つ以上の異なる半径方向の歪み勾配を計算する第2のステップと、
前記2つ以上の破断限界歪みと前記2つ以上の歪み勾配との関係から破断クライテリアを計算する第3のステップと、
有限要素法による数値解析から得られる最大主歪み及び隣接する要素間における歪み勾配が前記破断クライテリアに達したときに破断したと評価する第4のステップと
をコンピュータに実行させ、
前記第1のステップ〜前記第4のステップを一連の工程として実行して得られた解析結果に基づいて、前記材料の成形性を評価することを特徴とする成形性評価プログラム。
9.前記第1のステップにおいて、前記2つ以上の異なる頂角の円錐穴広げ試験から得られる夫々の前記穴広げ率を入力することを特徴とする8.に記載の成形性評価プログラム。
10.前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円錐工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする8.に記載の成形性評価プログラム。
11.前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円筒工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする8.に記載の成形性評価プログラム。
12.前記第2のステップにおいて、前記穴広げ率をλ(%)、前記板端部の破断限界歪みを(εθr=0、前記歪み勾配をdεθ/dr、円錐工具又は円筒工具の半頂角をφ(deg.)、素板の穴直径をd0(mm)、変形後の位置を表す半径座標をr(mm)とすると、
前記材料の板端部で破断が生じたときの穴の直径dは、(1+λ/100)d0として得られ、
(εθr=0=ln(d/d0
Figure 0006852426
として計算することを特徴とする8.〜11.のいずれか1項に記載の成形性評価プログラム。
13.前記歪み勾配であるdεθ/drを計算する際に、有限要素法の解析に用いる要素サイズを基準長さdrとすることを特徴とする8.〜12.のいずれか1項に記載の成形性評価プログラム。
14.前記2つ以上の破断限界歪み及び歪み勾配から前記破断クライテリアεcr=f(dεθ/dr)を計算し、有限要素法による数値解析から、隣接する要素間での前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11を求め、前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11が前記破断クライテリアに達しているか否かの指標としてε11/εcrを計算し、その結果をコンター表示することを特徴とする8.〜13.のいずれか1項に記載の成形性評価プログラム。
15.8.〜14.のいずれか1項に記載の成形性評価プログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
本発明によれば、板端部からの歪み勾配の影響と打抜き状態を考慮した伸びフランジ破断のクライテリアを推定することができ、これを成形シミュレーションに適用することで伸びフランジ破断の危険性を定量評価することができる。これにより、高強度鋼板適用時の成形課題の一つである伸びフランジ破断を未然に回避し、高強度で軽量な部品のプレス成形を実現することが可能となる。
従来技術の説明に用いる成形限界線図(FLD)を示す特性図である。 一様応力状態における局部くびれを説明するための特性図である。 伸びフランジ部の板端部から内側に向けた歪み勾配を説明するための特性図である。 穴広げ試験を説明するための模式図である。 穴広げ試験における穴縁の歪み及び歪み勾配を求めるための模式図である。 穴広げ試験における穴縁の歪み及び歪み勾配の計算結果と実験との比較を示す特性図である。 破断クライテリアである歪み勾配と破断限界歪みとの関係を説明するための特性図である。 本実施形態の一つであり、異なる穴径の素材を円錐工具により試験したときの穴広げ率から穴縁の歪み及び歪み勾配を計算した結果とこれを用いた破断クライテリアを説明する特性図である。 穴広げ試験のときに用いる円筒工具を説明する模式図である。 本実施形態の一つであり、異なる穴径の素材を円筒工具により試験したときの穴広げ率から穴縁の歪み及び歪み勾配を計算した結果とこれを用いた破断クライテリアを説明する特性図である。 本実施形態の破断予測方法を示すフローチャートである。 本実施形態により、破断危険性を評価した結果を示す模式図である。 パーソナルユーザ端末装置の内部構成を示す模式図である。
以下、図面を用いて、本実施形態を詳細に説明する。
(破断予測方法の基本骨子)
鋼板の端部から亀裂が生じる伸びフランジ成形では、フランジ端部から内側に向けて歪みが小さくなるため、材料端部は内側の拘束を受け、くびれの発生が抑制される。また、せん断加工時の鋼板端部に導入される微視的損傷の影響により、破断機構は複雑であり、従来のFLDによる破断予測では予測精度を確保することができない。そこで、板端部からの歪み勾配と端部の打抜き状態を考慮した穴広げ試験の結果を活用した伸びフランジ破断予測技術を想到した。
エッジ部の成形限界はフランジ端の周方向の伸び歪みが限界値に達したときに発生するため、フランジ端の周方向伸び歪みが成形難易度を示す指標として考えられ、通常、円錐パンチ穴広げ試験により得られる穴広げ率λで評価される。穴広げ試験は、一般に、素材に例えば直径10mmの円形の穴をあけ、例えば頂角60゜の円錐パンチにより穴径を広げ、穴縁端部で亀裂が生じて板厚を貫通したときを破断限界とする。その時の直径d及び初期径d0から、穴広げ率{(d−d0)/d0}×100を計算することで得られる(図4)。このときの穴は、一般に、パンチとダイスとの間のクリアランスが例えば板厚の12%になるように設定してせん断される。このため、フランジ端の破断面にはディンプル、亀裂、ボイド等の微視的損傷が観察される。従って、穴広げ試験は、せん断加工時に導入される微視的損傷を考慮した試験方法でもある。
この穴広げ率は、穴径、パンチ径等の試験寸法やパンチ底形状によって大きく変化することが知られている。これは、伸びフランジ成形限界が、フランジ端部から内側へ向けて小さくなる歪み勾配の影響を強く受けることによる。従って、実部品の伸びフランジ成形で導入される歪み状態(板端部から内側への歪み勾配)と穴広げ試験での歪み勾配とを合わせた評価試験で破断限界を評価する必要がある。しかしながら、実部品で問題となる伸びフランジ破断では、歪み分布は部品形状やブランクの形状によって大きく変化することが多い。実験から得られる穴広げ率を用いて破断を評価する場合、歪み勾配に応じて穴径やパンチ径等を無数に変化させ破断限界を準備せざるを得ない。
一方、実部品の伸びフランジ成形で導入される歪み状態(板端部から内側への歪み勾配、端部に沿った歪み勾配)と破断限界の関係を検討する方法として、サイドベンド試験法が提案されている(例えば、特許文献1,2を参照)。その特徴は、歪みの計測に画像処理システムを導入すると共に、端面の亀裂発生の挙動を観察するためのカメラを備えているところにある。この方法から得られた破断限界歪みと歪み分布の関係から破断限界曲面を定義し、伸びフランジ破断の予測に活用する方法が提案されている。
しかしながら、この方法により端部から半径方向の歪み勾配を得るためには、予め素材の表面にエッチング等によりサークル状或いは格子状の模様を描いておき、変形後にサークルの変形量から歪み勾配を測定する必要がある。他には、数値シミュレーションにより種々の歪み勾配の破断限界を得る方法もあるが、実部品で問題となる歪み分布に対応した無数の試験片形状に対して破断限界を得るには煩雑であり、実用上、難しい。
本実施形態では、穴広げ試験、具体的には円錐穴広げ試験から得られる穴広げ率の試験値から板端部の破断限界歪みと半径方向の歪み勾配を容易に得られる方法を検討する。円錐ポンチによる穴広げ試験の概要を図5に示す。ここで、円錐ポンチの円錐面の法線とz軸とのなす角をφ(円錐ポンチの半頂角は90°−φとなる)とする。すると、軸対称に変形中の材料の要素のつりあい方程式は以下で与えられる。
Figure 0006852426
ここで、rは変形後の要素の位置を示す半径座標、tは変形後の板厚、σθ及びσφはそれぞれ円周方向及び半径方向の応力、μは摩擦係数を表す。
続いて、円周方向の歪みεφは、半径方向の歪みεθは、それぞれ以下で与えられる。
Figure 0006852426
但し、sは板の要素の変形前の半径座標を表す。これらから次式の歪みの適合条件式が得られる。
Figure 0006852426
更に、変形は全歪み理論によるものとし、材料の加工硬化特性はn乗硬化則で近似すると仮定する。
Figure 0006852426
すると、円周方向、半径方向の歪み分布はそれぞれ次式で与えられる。
Figure 0006852426
ここで、1次の微係数(歪み勾配)、2次の微係数はそれぞれ次式で与えられる。
Figure 0006852426
ここで、例えば、初期穴径d0の素板が穴広げ試験により直径dとなったとき端部の歪みは次式により得られ、
Figure 0006852426
これと式(6),(9)〜(11)から半径方向歪み分布σθが計算できる。
以下、1.4mm厚で980MPa級の高強度鋼板を例に採り、破断限界のクライテリアを計算する方法を説明する。
先ず、オフラインで異なる頂角の円錐穴広げ試験を行う。ここでは、その例として30゜円錐及び60゜円錐のパンチを用いた。このときのダイスとパンチとの間のクリアランスを板厚の12%に設定し、素板中央に直径10mmの穴を打抜いた。この素板を穴広げ試験に供し、穴縁端部で亀裂が板厚を貫通したときを破断限界とし、その時の直径d及び初期径d0から穴広げ率{(d−d0)/d0}×100を計算した。その結果、30゜円錐及び60゜円錐の穴広げ率は、それぞれ42%及び35%であった。
この結果及び式(6),(9)〜(11)から、それぞれの円錐パンチを用いたときの半径方向歪み分布σθを計算することができる。その結果を図6に示す。実験値は、評点間距離0.5mmの同心円のスクライブドサークルを素材にエッチングし、試験後のサークルの変形量から求めた。その結果、計算結果は実験を良好な精度で再現することを確認した。即ち、30゜円錐、60゜円錐いずれの結果も、穴縁で周方向歪みεθは最大値を観測し、穴縁から内側に向けて単調に減少した。また、穴縁では60゜円錐よりも30゜円錐の方が周方向歪みεθは大きいが、穴縁から内側に向けてその傾向は逆転し、穴縁から十分離れた位置では60゜円錐の方が大きな値を示した。
続いて、上記の計算結果から、歪み勾配に対応した周方向破断歪みを計算する。周方向歪み分布は穴縁から離れると共に緩やかに減少しているため、基準となる大きさdrにより歪み勾配は変化する。そのため、歪み勾配dεθ/drを計算する際に、有限要素法の解析に用いる要素サイズを基準長さdrとすることで破断予測の精度が高まる。なお、それほど高い精度での評価が必要なければ、式(9)の1次微係数を用いても良い。
このようにして得られた2つ以上の歪み勾配|dε/dr|と穴縁での破断限界歪みεcr=(εθr=0との関係から、直線近似又はそれらを結ぶ多直線データを破断クライテリアεcr=f(dεθ/dr)とする(図7)。
続いて、成形性を評価する部品の成形解析を有限要素法による数値解析により行い、その結果の歪み分布から隣接する要素間での歪み勾配dε11/dr及び最大主歪みε11を求め、これらが破断クライテリアに達していれば破断の危険性が高いと判断する。更に、破断の危険性を定量評価するための指標としてε11/εcrを計算し(図7中のOR/OAの大きさに相当する。)、その結果をコンター表示する。
(破断予測方法のその他の例)
破断クライテリアを求めるときの穴広げ試験としては、2つ以上の異なる頂角の円錐穴広げの他に、異なる穴径の素材を単一形状の円錐パンチ又は円筒パンチにより拡大する穴広げ試験を行うようにしても良い。異なる穴径の場合には、式(13)の初期穴径d0をそれぞれの穴径とすれば良い。
[1]円錐パンチを用いた穴広げ試験
円錐パンチを用いた穴広げ試験の例として、60゜の円錐パンチを用い、素板中央に直径10mmの穴と直径50mmの穴を打抜いた素材を穴広げ試験に供した。その結果、直径10mmの穴及び直径50mmの素材を用いた穴広げ率は、それぞれ42%及び18%であった。この結果及び式(6),(9)〜(11),(13)から、それぞれの初期穴径を用いたときの半径方向ひずみ分布を計算することができる。その結果を図8(a)に示す。実験値は、標点間距離0.5mmの同心円のスクライブドサークルを素材にエッチングし、試験のサークルの変形量から求めた。その結果、計算値は実験を良好な精度で再現できることを確認した。ここで得られた2つ以上のひずみ勾配と穴縁での破断限界ひずみとの関係から決定した破断クライテリアを図8(b)に示す。
[2]円筒パンチを用いた穴広げ試験
円筒パンチを用いる場合には、式(1)〜(12)に対して2φ=180゜とすれば同様に、破断クライテリアを計算することができる。ここでは、図9に示すように、異なる穴径の素材を単一形状の円筒パンチ(パンチ肩半径10mm、ダイス肩半径10mm)を用い、素板中央に直径10mmの穴(d=10mm)と直径50mm(d=50mm)の穴を打抜いた素材を穴広げ試験に供した。その結果として、式(6),(9)〜(11),(13)に対してφ=180゜として計算した穴縁の破断限界ひずみと穴縁からの距離との関係を図10(a)に、破断限界ひずみとひずみ勾配との関係を図10(b)に示す。
(本実施形態の具体的構成)
図11を用いて、自動車部品の成形性を評価する例として、本実施形態の具体的構成を説明する。
自動車部品の成形性を評価するにあたり、先ず当該自動車の構造を設定した後(ステップS1)、CADを用いて自動車部品の形状を設定し(ステップS2)、3次元の部品形状をコンピュータ上に記録する(ステップS3)。ここで、金型を用いてプレス加工ができるかどうか評価するために、金型CADにて金型を設計し(ステップS4)、目的に応じてソフトを選択してコンピュータ上に記録する。次に、成形解析により成形性の評価を行う(ステップS8)。そのために先ず、成形性評価に供する部品の材料パラメータ、板厚、成形条件を設定する(ステップS5)。
続いて、オフラインで試験した2つ以上の異なる頂角の円錐穴広げ試験の結果を入力する(ステップS6)。ここで、2つ以上の異なる頂角の円錐穴広げ試験の代わりに、2つ以上の異なる初期穴径の素材を用いた単一形状の頂角の円錐工具による穴広げ試験を行う場合には、その試験から得られた穴広げ率をステップS6において入力する。また、2つ以上の異なる頂角の円錐穴広げ試験の代わりに、2つ以上の異なる初期穴径の素材を用いた単一形状の頂角の円筒工具による穴広げ試験を行う場合には、その試験から得られた穴広げ率をステップS6において入力する。
続いて、この2つ以上の穴広げ率の値を用いて前述した方法で破断クライテリアεcr=f(dεθ/dr)を計算する(ステップS7)。更に、破断の評価は、このようにして求めた破断クライテリアと変形過程の有限要素法によるシミュレーションの結果から得られる各要素の最大主歪み、隣接する歪み勾配との位置関係を比較することで評価する。この歪みが破断クライテリアに達したときに破断、若しくは、その危険性が高いと判断する(ステップS8〜S10)。具体的には、図7に記載の歪み勾配−破断限界歪みの関係において、有限要素法により得られた要素の歪み状態をR、歪み0の状態OとRを結ぶ直線とクライテリアとの交点をAとしたとき、破断危険率はOR/OAとして定量化することができる。
最後に、破断危険率をコンター表示し(ステップS11)、成形可否を判断する。具体的な実施として、フランジ高さH=30mm、コーナーR=30mm、開き角120°の鞍型形状部品を、1.4mm厚の980MPa級高強度鋼板により成形したときの破断危険性を評価した結果を図12に示す。
以上説明したように、本実施形態によれば、板端部からの歪み勾配の影響と打抜き状態を考慮した伸びフランジ破断のクライテリアを推定することができ、これを成形シミュレーションに適用することで伸びフランジ破断の危険性を定量評価することができる。これにより、高強度鋼板適用時の成形課題の一つである伸びフランジ破断を未然に回避し、高強度で軽量な部品のプレス成形を実現することが可能となる。
(その他の実施形態)
上述した本実施形態による成形性予測評価方法の各ステップ(図11のステップS1〜S11等)は、コンピュータのRAMやROM等に記録されたプログラムが動作することによって実現できる。このプログラム及び当該プログラムを記録したコンピュータ読み取り可能な記録媒体は本実施形態に含まれる。
具体的に、上記のプログラムは、例えばCD−ROMのような記録媒体に記録し、或いは各種伝送媒体を介し、コンピュータに提供される。上記のプログラムを記録する記録媒体としては、CD−ROM以外に、フレキシブルディスク、ハードディスク、磁気テープ、光磁気ディスク、不揮発性メモリカード等を用いることができる。他方、上記のプログラムの伝送媒体としては、プログラム情報を搬送波として伝搬させて供給するためのコンピュータネットワークシステムにおける通信媒体を用いることができる。ここで、コンピュータネットワークとは、LAN、インターネットの等のWAN、無線通信ネットワーク等であり、通信媒体とは、光ファイバ等の有線回線や無線回線等である。
また、本実施形態に含まれるプログラムとしては、供給されたプログラムをコンピュータが実行することにより本実施形態の機能が実現されるようなもののみではない。例えば、そのプログラムがコンピュータにおいて稼働しているOS(オペレーティングシステム)或いは他のアプリケーションソフト等と共同して本実施形態の機能が実現される場合にも、かかるプログラムは本実施形態に含まれる。また、供給されたプログラムの処理の全て或いは一部がコンピュータの機能拡張ボードや機能拡張ユニットにより行われて本実施形態の機能が実現される場合にも、かかるプログラムは本実施形態に含まれる。
例えば、図13は、パーソナルユーザ端末装置の内部構成を示す模式図である。この図13において、1200はCPU1201を備えたパーソナルコンピュータ(PC)である。PC1200は、ROM1202またはハードディスク(HD)1211に記憶された、又はフレキシブルディスクドライブ(FD)1212より供給されるデバイス制御ソフトウェアを実行する。このPC1200は、システムバス1204に接続される各デバイスを総括的に制御する。
PC1200のCPU1201、ROM1202またはハードディスク(HD)1211に記憶されたプログラムにより、本実施形態の図11におけるステップS1〜S11の手順等が実現される。
1203はRAMであり、CPU1201の主メモリ、ワークエリア等として機能する。1205はキーボードコントローラ(KBC)であり、キーボード(KB)1209や不図示のデバイス等からの指示入力を制御する。
1206はCRTコントローラ(CRTC)であり、CRTディスプレイ(CRT)1210の表示を制御する。1207はディスクコントローラ(DKC)である。DKC1207は、ブートプログラム、複数のアプリケーション、編集ファイル、ユーザファイル、ネットワーク管理プログラム等を記憶するハードディスク(HD)1211、及びフレキシブルディスク(FD)1212とのアクセスを制御する。ここで、ブートプログラムとは、パソコンのハードやソフトの実行(動作)を開始する起動プログラムである。
1208はネットワーク・インターフェースカード(NIC)であり、LAN1220を介して、ネットワークプリンタ、他のネットワーク機器、或いは他のPCと双方向のデータのやり取りを行う。
なお、パーソナルユーザ端末装置を用いる代わりに、成形性予測評価方法に特化された所定の計算機等を用いても良い。

Claims (13)

  1. 材料のプレス成形の可否を評価する方法であって、
    穴広げ試験から得られる2つ以上の異なる穴広げ率を入力する第1のステップと、
    前記2つ以上の異なる穴広げ率のデータから、前記2つ以上の異なる破断限界歪み及び前記2つ以上の異なる半径方向の歪み勾配を計算する第2のステップと、
    前記2つ以上の破断限界歪みと前記2つ以上の歪み勾配との関係から破断クライテリアを計算する第3のステップと、
    有限要素法による数値解析から得られる最大主歪み及び隣接する要素間における歪み勾配が前記破断クライテリアに達したときに破断したと評価する第4のステップと
    を含み、
    前記第2のステップにおいて、前記穴広げ率をλ(%)、板端部の前記破断限界歪みを(ε θ r=0 、前記歪み勾配をdε θ /dr、円錐工具又は円筒工具の半頂角をφ(deg.)、素板の穴直径をd 0 (mm)、変形後の位置を表す半径座標をr(mm)とすると、
    前記材料の板端部で破断が生じたときの穴の直径dは、(1+λ/100)d 0 として得られ、
    (ε θ r=0 =ln(d/d 0
    Figure 0006852426
    として計算し、
    前記第1のステップ〜前記第4のステップを一連の工程として実行して得られた解析結果に基づいて、前記材料の成形性を評価することを特徴とする成形性評価方法。
  2. 前記第1のステップにおいて、前記2つ以上の異なる頂角の円錐穴広げ試験から得られる夫々の前記穴広げ率を入力することを特徴とする請求項1に記載の成形性評価方法。
  3. 前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円錐工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする請求項1に記載の成形性評価方法。
  4. 前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円筒工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする請求項1に記載の成形性評価方法。
  5. 前記歪み勾配であるdεθ/drを計算する際に、有限要素法の解析に用いる要素サイズを基準長さdrとすることを特徴とする請求項1〜のいずれか1項に記載の成形性評価方法。
  6. 前記2つ以上の破断限界歪み及び歪み勾配から前記破断クライテリアεcr=f(dεθ/dr)を計算し、有限要素法による数値解析から、隣接する要素間での前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11を求め、前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11が前記破断クライテリアに達しているか否かの指標としてε11/εcrを計算し、その結果をコンター表示することを特徴とする請求項1〜のいずれか1項に記載の成形性評価方法。
  7. 材料のプレス成形の可否を評価するためのプログラムであって、
    穴広げ試験から得られる2つ以上の異なる穴広げ率を入力する第1のステップと、
    前記2つ以上の異なる穴広げ率のデータから、前記2つ以上の異なる破断限界歪み及び前記2つ以上の異なる半径方向の歪み勾配を計算する第2のステップと、
    前記2つ以上の破断限界歪みと前記2つ以上の歪み勾配との関係から破断クライテリアを計算する第3のステップと、
    有限要素法による数値解析から得られる最大主歪み及び隣接する要素間における歪み勾配が前記破断クライテリアに達したときに破断したと評価する第4のステップと
    をコンピュータに実行させ、
    前記第2のステップにおいて、前記穴広げ率をλ(%)、板端部の前記破断限界歪みを(ε θ r=0 、前記歪み勾配をdε θ /dr、円錐工具又は円筒工具の半頂角をφ(deg.)、素板の穴直径をd 0 (mm)、変形後の位置を表す半径座標をr(mm)とすると、
    前記材料の板端部で破断が生じたときの穴の直径dは、(1+λ/100)d 0 として得られ、
    (ε θ r=0 =ln(d/d 0
    Figure 0006852426
    として計算し
    前記第1のステップ〜前記第4のステップを一連の工程として実行して得られた解析結果に基づいて、前記材料の成形性を評価することを特徴とする成形性評価プログラム。
  8. 前記第1のステップにおいて、前記2つ以上の異なる頂角の円錐穴広げ試験から得られる夫々の前記穴広げ率を入力することを特徴とする請求項に記載の成形性評価プログラム。
  9. 前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円錐工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする請求項に記載の成形性評価プログラム。
  10. 前記第1のステップにおいて、前記2つ以上の異なる初期穴径の前記材料を用いた単一形状の頂角の円筒工具による穴広げ試験から得られる前記穴広げ率を入力することを特徴とする請求項に記載の成形性評価プログラム。
  11. 前記歪み勾配であるdεθ/drを計算する際に、有限要素法の解析に用いる要素サイズを基準長さdrとすることを特徴とする請求項10のいずれか1項に記載の成形性評価プログラム。
  12. 前記2つ以上の破断限界歪み及び歪み勾配から前記破断クライテリアεcr=f(dεθ/dr)を計算し、有限要素法による数値解析から、隣接する要素間での前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11を求め、前記歪み勾配であるdε11/dr及び前記最大主歪みであるε11が前記破断クライテリアに達しているか否かの指標としてε11/εcrを計算し、その結果をコンター表示することを特徴とする請求項11のいずれか1項に記載の成形性評価プログラム。
  13. 請求項12のいずれか1項に記載の成形性評価プログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
JP2017018964A 2016-02-05 2017-02-03 成形性評価方法、プログラム及び記録媒体 Active JP6852426B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016021311 2016-02-05
JP2016021311 2016-02-05

Publications (2)

Publication Number Publication Date
JP2017140653A JP2017140653A (ja) 2017-08-17
JP6852426B2 true JP6852426B2 (ja) 2021-03-31

Family

ID=59628165

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017018964A Active JP6852426B2 (ja) 2016-02-05 2017-02-03 成形性評価方法、プログラム及び記録媒体

Country Status (1)

Country Link
JP (1) JP6852426B2 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7110976B2 (ja) * 2017-12-27 2022-08-02 日本製鉄株式会社 成形性評価方法、プログラム及び記録媒体
CN108614945B (zh) * 2018-05-10 2022-01-28 哈尔滨电气股份有限公司 一种背部承压平盖设计方法
JP6819832B1 (ja) 2019-03-14 2021-01-27 Jfeスチール株式会社 伸びフランジ割れ評価方法、金属板の選定方法、プレス金型の設計方法、部品形状の設計方法、及びプレス部品の製造方法
CN110376059B (zh) * 2019-07-11 2021-11-19 上海理工大学 基于静强度的结构全场轻量化定量评价方法
JP7382913B2 (ja) 2020-08-20 2023-11-17 株式会社神戸製鋼所 数値シミュレーションの欠陥発生リスク評価方法および装置
KR102459164B1 (ko) * 2021-06-23 2022-10-25 울산과학기술원 강재의 신장 플랜지 성형성 평가방법
JP7217322B1 (ja) * 2021-09-01 2023-02-02 Jfeスチール株式会社 金属薄板の成形限界判定方法および成形限界判定システムならびにコンピュータプログラム

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4935713B2 (ja) * 2008-02-27 2012-05-23 Jfeスチール株式会社 プレス品のせん断縁における成形可否判別方法
JP5472518B1 (ja) * 2012-11-19 2014-04-16 Jfeスチール株式会社 伸びフランジの限界ひずみ特定方法およびプレス成形可否判定方法
JP6314626B2 (ja) * 2014-04-21 2018-04-25 新日鐵住金株式会社 プレス成形性の評価方法、装置、プログラム及びコンピュータ読み取り可能な記憶媒体
EP3165298A4 (en) * 2014-07-02 2018-03-07 Nippon Steel & Sumitomo Metal Corporation Stretch-flange crack prediction method, stretch-flange crack prediction device, computer program, and recording medium

Also Published As

Publication number Publication date
JP2017140653A (ja) 2017-08-17

Similar Documents

Publication Publication Date Title
JP6852426B2 (ja) 成形性評価方法、プログラム及び記録媒体
JP6828476B2 (ja) エッジ部破断予測方法、プログラム及び記録媒体
JP6897413B2 (ja) 成形性評価方法、プログラム及び記録媒体
JP5445381B2 (ja) 材料の曲げ破断予測方法および装置、ならびにプログラムおよび記録媒体
JP5146395B2 (ja) 歪勾配を考慮した伸びフランジ割れの推定方法およびプレス成形シミュレーションの伸びフランジ割れ判定システム
EP3016009B1 (en) Method for determining bending fracture in metal plate, program, and storage medium
JP4865036B2 (ja) 破断予測方法、演算処理装置、プログラム及び記録媒体
JP4621217B2 (ja) 破断予測方法及び装置、並びにプログラム及び記録媒体
JP7110976B2 (ja) 成形性評価方法、プログラム及び記録媒体
US8990028B2 (en) Fracture prediction method, device, a program arrangement and computer-accessible medium therefor
Huang et al. The formability limitation of the hole-flanging process
JP4980499B2 (ja) 破断判定方法、破断判定装置、プログラムおよびコンピュータ読み取り可能な記録媒体
JP4621216B2 (ja) 破断限界取得方法及び装置、並びにプログラム及び記録媒体
Singh et al. Prediction of earing defect and deep drawing behavior of commercially pure titanium sheets using CPB06 anisotropy yield theory
JP4880043B2 (ja) 破断予測方法、破断予測システム、プログラム及び記録媒体
JP6314626B2 (ja) プレス成形性の評価方法、装置、プログラム及びコンピュータ読み取り可能な記憶媒体
WO2016002880A1 (ja) 伸びフランジ割れ予測方法、伸びフランジ割れ予測装置、コンピュータープログラム、及び記録媒体
JP6819832B1 (ja) 伸びフランジ割れ評価方法、金属板の選定方法、プレス金型の設計方法、部品形状の設計方法、及びプレス部品の製造方法
JP6176410B2 (ja) 破断予測方法、プログラム、記録媒体及び演算処理装置
Roy et al. Experimental and numerical investigation of deformation characteristics during tube spinning
JP7206902B2 (ja) 成形性評価方法、プログラム及び記録媒体
Chen An analysis of forming limit in the elliptic hole-flanging process of sheet metal
JP5900751B2 (ja) 曲げ内側割れの評価方法および予測方法
JP2020159834A (ja) 破断クライテリア解析方法、破断クライテリア解析プログラム、及び破断クライテリア解析システム
JP6303815B2 (ja) 鍛造割れ予測方法、鍛造割れ予測プログラム、および記録媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191007

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200820

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201006

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201202

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210209

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210222

R151 Written notification of patent or utility model registration

Ref document number: 6852426

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151