WO2012133867A1 - プレス成形解析方法 - Google Patents

プレス成形解析方法 Download PDF

Info

Publication number
WO2012133867A1
WO2012133867A1 PCT/JP2012/058774 JP2012058774W WO2012133867A1 WO 2012133867 A1 WO2012133867 A1 WO 2012133867A1 JP 2012058774 W JP2012058774 W JP 2012058774W WO 2012133867 A1 WO2012133867 A1 WO 2012133867A1
Authority
WO
WIPO (PCT)
Prior art keywords
press
stress distribution
stress
analysis
mold release
Prior art date
Application number
PCT/JP2012/058774
Other languages
English (en)
French (fr)
Inventor
亮伸 石渡
裕隆 狩野
Original Assignee
Ｊｆｅスチール株式会社
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 Ｊｆｅスチール株式会社 filed Critical Ｊｆｅスチール株式会社
Priority to EP12763272.7A priority Critical patent/EP2692454B1/en
Priority to MX2013011153A priority patent/MX2013011153A/es
Priority to KR1020137026029A priority patent/KR101544464B1/ko
Priority to US14/007,974 priority patent/US9410855B2/en
Priority to CN201280016493.7A priority patent/CN103459058B/zh
Publication of WO2012133867A1 publication Critical patent/WO2012133867A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/04—Measuring force or stress, in general by measuring elastic deformation of gauges, e.g. of springs
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21D—WORKING OR PROCESSING OF SHEET METAL OR METAL TUBES, RODS OR PROFILES WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21D22/00—Shaping without cutting, by stamping, spinning, or deep-drawing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82—Elements for improving aerodynamics

Definitions

  • the present invention relates to a press molding analysis method, and more particularly, to a method for accurately analyzing in a short time how and in which region of a press molded product the stress affects the springback.
  • a high-strength steel sheet has a large elastic recovery after press forming (hereinafter referred to as springback, and may be abbreviated as SB), and therefore it is difficult to ensure the dimensional accuracy of parts. For this reason, it is necessary to repeatedly adjust the shape of the press mold and adjust the mold several times until the normal dimensional accuracy is achieved.
  • step (1) Analysis of deformation, stress and strain of material under restraint by press mold
  • step (2) Analysis of springback in a state released from restraint by press mold
  • step (1) is finite
  • step (2) Analysis of springback in a state released from restraint by press mold
  • step (1) is finite
  • step (2) the analysis is performed by the springback theoretical formula or the finite element method.
  • Patent Documents 1 and 2 have studied the influence of the shape of a molding object and the molding conditions on the springback.
  • Patent Document 1 calculates how a defined amount of springback before and after changing the residual stress distribution changes by changing the residual stress distribution in a certain area of the press molding object, A method for predicting the impact on a region of springback is disclosed.
  • the present invention has been made to solve the above-described problems, and the stress change before release is performed in consideration of the stress (residual stress) distribution after the press-molded product is released from the press die.
  • the stress change before release is performed in consideration of the stress (residual stress) distribution after the press-molded product is released from the press die.
  • the gist of the present invention is as follows.
  • the shape of the press-formed product before release and the data of the residual stress distribution are calculated, and the stress in the local coordinate system before release by coordinate transformation based on the data before release.
  • a springback analysis is performed based on the data before release, the residual stress distribution of the press-formed product after release is calculated, and the residual stress distribution (b) of the press-formed product in the local coordinate system by coordinate transformation )
  • shape data (d) of the press-formed product after mold release (3)
  • the difference (ab) between the stress distributions (a) and (b) is calculated, and the difference is used as a spring bag effective stress (SB effective stress).
  • a press forming analysis method characterized by comprising:
  • the shape of the press-molded product before release and the data of residual stress distribution are calculated, and the stress in the local coordinate system before release by coordinate transformation based on the data before release.
  • a springback analysis is performed based on the data before release, the residual stress distribution of the press-formed product after release is calculated, and the residual stress distribution (b) of the press-formed product in the local coordinate system by coordinate transformation )
  • shape data (d) of the press-formed product after mold release (3)
  • the difference (ab) between the stress distributions (a) and (b) is calculated, and the difference is used as a spring bag effective stress (SB effective stress).
  • FIG. 4 is a flowchart showing the processing sequence of the press forming analysis method according to the first embodiment of the press forming analyzing method according to the present invention.
  • the springback analysis is performed after changing or removing the springback effective stress in the springback analysis target region, and the press-molded product after the mold release is performed. This is an analysis method for evaluating the influence of stress in the springback analysis target region by comparing the shapes.
  • the analysis calculation in the press forming analysis method shown in FIG. 4 is performed using a finite element method (FEM) analysis system.
  • FEM finite element method
  • step S1 First, in step S1, the stress distribution of the press-formed product before release is obtained.
  • step S2 the stress expressed in the global coordinate system obtained in step S1 is defined as a local coordinate system that does not change before and after mold release, and the stress is coordinate-transformed to obtain a stress distribution in the local coordinate system ( a) is determined. Since press-molded products are rotated by springback after mold release, if the stress expressed in the global coordinate system is directly differentiated, the stress contributing to the deformation cannot be obtained accurately. Therefore, for each finite element that gives a stress distribution, a local coordinate system that does not change before and after mold release, which is a method generally used in finite element analysis, is defined, and the stress difference is calculated before and after mold release after coordinate transformation of stress. Do it.
  • the coordinate transformation of the stress is performed by converting the three orthogonal unit vectors of the local coordinate system to (l 1 , m 1 , n 1 ), (l 2 , m 2 , n 2 ), (l 3 , m 3 , n 3 ) Then, it is given by equation (1).
  • the local coordinate system used in the press forming simulation generally uses four or three shell elements with no nodes in the thickness direction. Therefore, it is common to select a plate thickness direction, that is, a surface normal direction formed by four or three nodes as one third axis of the coordinate system.
  • a plate thickness direction that is, a surface normal direction formed by four or three nodes as one third axis of the coordinate system.
  • the four points are not necessarily on the same plane, so (1) the direction obtained from the outer product of the diagonal lines is the plate thickness direction, and (2) the surface normal direction of the approximate plane is the plate There are methods such as the thickness direction.
  • the first axis and the second axis include (1) a method in which the direction from the first node to the second node of the shell element is the first axis, and (2) a method in which the first axis and the second axis are made from a diagonal vector.
  • S6 The shape data (d) of the press-molded product after mold release is obtained from the springback analysis of S3.
  • SB effective stress The difference (ab) between the stress distribution (a) before mold release and the residual stress distribution (b) after mold release is called springback effective stress (hereinafter referred to as SB effective stress), and the distribution is obtained. .
  • step S11 A springback analysis is performed from the stress distribution obtained in step S10.
  • FIG. 5 is a flowchart showing the processing sequence of the press forming analysis method according to the second embodiment of the press forming analyzing method according to the present invention.
  • the stress distribution obtained by adding the residual stress of the local coordinate system to the spring bag effective stress distribution of the local coordinate system obtained by partially changing or removing the springback effective stress obtained in the first embodiment is an analysis method for evaluating the influence of the stress in the region subject to springback analysis by performing a springback analysis after converting to a global coordinate system and comparing it with the shape of the press-formed product after mold release.
  • FEM finite element method
  • step S21 First, in step S21, the stress distribution of the press-formed product before mold release is obtained.
  • step S22 the local coordinate system that does not change before and after the mold release is defined for the stress expressed in the global coordinate system obtained in step S21, and the stress is coordinate-transformed to obtain a stress distribution ( a) is determined. Since the method for obtaining the local coordinate system is the same as that in the first embodiment, it is omitted in this section.
  • SB effective stress The difference (ab) between the stress distribution (a) before mold release and the residual stress distribution (b) after mold release is called springback effective stress (hereinafter referred to as SB effective stress), and the distribution is obtained. .
  • S28 The coordinate transformation from the local coordinate system to the global coordinate system is performed to obtain the SB effective stress distribution in the global coordinate system.
  • press molding of a front pillar inner part will be described.
  • Figure 1 shows the general shape of this part. Using a 980 MPa class high-tensile cold-rolled steel sheet having a thickness of 1.6 mm, press molding and trim processing were performed by a mold, and the front pillar inner part was finished through bend processing. Bending is performed along the bend line shown in FIG.
  • the press forming analysis analyzed how much the stresses in the analysis target areas 1 to 6 shown in FIG. 2 are affected when the front pillar inner part shown in FIG. 1 is formed.
  • the analysis was performed along the steps shown in FIGS. 4 and 5 using a commercially available finite element method (FEM) analysis system. Specifically, after bending, the stress in each region (1 to 6) shown in FIG. 2 was removed, and the influence on the Z-direction displacement at points A and B was investigated. In addition, the displacement in the Z direction of the point A after the usual springback analysis (the stress is not removed in each region) is ⁇ 16.5 mm, and the point B is ⁇ 9.7 mm. In the conventional example (Patent Document 1), the springback calculation is performed by removing all the stresses regarding the analysis target region (any one of 1 to 6). Invention Example 1 takes out only the SB effective stress of the present invention (S8 in FIG.
  • FEM finite element method
  • Inventive example 2 is obtained by removing only the SB effective stress in the analysis target region (any one of 1 to 6) (S29 in FIG. 5) and performing a springback calculation (S35 in FIG. 5).
  • the horizontal axis represents the number of the analysis target area
  • the vertical axis represents the displacement (mm) in the Z direction of the point A.
  • the influence of the analysis target area 6 is the largest, but in the first and second invention examples, the influence of the analysis target area 2 is the largest.
  • the influence of the analysis target area 6 is increased because the analysis target area 6 has a relatively complicated shape and a large residual stress after the spring bag, and therefore, the influence is eliminated. . Further, in the bend process, the analysis target region 4 to the analysis target region 6 are far from the bend line to which plastic deformation is applied. Therefore, it is considered that the influence of the analysis target region 6 is excessively evaluated in the conventional example.
  • 1 to 6 are codes representing the analysis target area.

Abstract

（１）離型前のプレス成形品の形状，残留応力分布のデータを算出し，離型前の局所座標系での応力分布（ａ）を算出する工程，（２）前記離型前のデータに基づいてスプリングバック解析を行い，離型後のプレス成形品の残留応力分布を算出し，局所座標系でのプレス成形品の残留応力分布（ｂ）と，離型後のプレス成形品の形状データ（ｄ）を算出する工程，（３）前記応力分布（ａ）と（ｂ）との差分（ａ－ｂ）を算出し，該差分をスプリングバッグ有効応力（ＳＢ有効応力）として全体座標系におけるＳＢ有効応力分布を算出し，このＳＢ有効応力分布のうち解析対象領域のＳＢ有効応力を変更または除去したＳＢ有効応力分布に基づいてスプリングバック解析を行い，離型後のプレス成形品の形状データ（ｃ）を算出する工程，（４）前記形状データ（ｃ）と形状データ（ｄ）との差分（ｃ－ｄ）を算出し，解析対象領域の全体形状に対する影響度合を判定する工程，とからなるプレス成形解析方法。

Description

プレス成形解析方法
　本発明は、プレス成形解析方法に関し、特にプレス成形品のどの領域の応力がどのようにスプリングバックに影響しているかを短時間で且つ、正確に解析する方法に関するものである。
　近年、自動車の軽量化を実現させるため、自動車部品への高強度鋼板の適用を拡大することが進められている。高強度鋼板は、プレス成形後の弾性回復（以下スプリングバックと呼ぶ。またＳＢと略記することもある）が、軟鋼板と比較して大きいため、部品の寸法精度を確保することが難しい。そのため、正規の寸法精度を達成するまで幾度となく、プレス金型の形状修正を繰返し、金型を調整する作業が必要となる。
　そこで、このようなプレス金型の修正作業を軽減するために、金型設計段階でスプリングバック量を予測する技術が望まれており、その解析システムが開発されてきている。
　それらの解析方法は、次の２つの段階よりなる。
（１）プレス金型による拘束の下での、材料の変形、応力、歪の解析
（２）プレス金型による拘束から開放された状態でのスプリングバックの解析
　ここで、段階（１）は有限要素法により解析され、段階（２）はスプリングバック理論式あるいは有限要素法による解析が行われている。
　こうしたスプリングバックは、プレス金型から成形対象物を離型する前に成形対象物に残留する残留応力の不均一分布に起因することが知られており、従来より、有限要素法等の数値解析方法を用いた予測がなされてきた。しかし、成形対象物の残留応力の内、どの領域におけるものがスプリングバックに支配的なのかを予測するのが困難であった。
　従来、このスプリングバックに及ぼす成形対象物の形状及び成形条件の影響について検討を行ったものとして、例えば、特許文献１及び２がある。
　特許文献１には、プレス成形対象物のある領域の残留応力分布を変更することにより、残留応力分布を変更する前後におけるスプリングバックに関するある定義された量がどのように変化するかを算出し、ある領域のスプリングバックへの影響を予測する方法が開示されている。
特開２００７−２２９７２４号公報
　特許文献１の応力の変更方法は、プレス成形対象物の金型からの離型前の応力のみを考慮して、応力の除去および変更を行い、その状態でスプリングバック解析を行って、応力を変更しないものと比較することで応力変更による影響を検討するものである。
　しかし、プレス成形対象物（以下プレス成形品と呼ぶ）にはプレス金型からの離型後も、即ちスプリングバック後にも応力は残留している。このため、プレス金型からの離型前の応力のみに注目してその影響を検証しても正確な各領域の応力の影響を評価することが困難な場合がある。
　そこで、本発明は、上述の問題点を解決するためになされたものであり、プレス成形品がプレス金型から離型後の応力（残留応力）分布を考慮して、離型前の応力変更により、どの領域における応力がスプリングバックにどう影響しているかを明らかにすることを課題とする。
　本発明の要旨は以下の通りである。
　第一の発明は、（１）離型前のプレス成形品の形状、残留応力分布のデータを算出し、前記離型前のデータに基づいて座標変換により離型前の局所座標系での応力分布（ａ）を算出する工程、
　（２）前記離型前のデータに基づいてスプリングバック解析を行い、離型後のプレス成形品の残留応力分布を算出し、座標変換により局所座標系でのプレス成形品の残留応力分布（ｂ）を算出するとともに、離型後のプレス成形品の形状データ（ｄ）を算出する工程、
　（３）前記応力分布（ａ）と（ｂ）との差分（ａ−ｂ）を算出し、該差分をスプリングバッグ有効応力（ＳＢ有効応力）として座標変換により全体座標系におけるＳＢ有効応力分布を算出し、このＳＢ有効応力分布のうち解析対象領域のＳＢ有効応力を変更または除去して、プレス成形品のＳＢ有効応力分布を算出し、スプリングバック解析を行い、離型後のプレス成形品の形状データ（ｃ）を算出する工程、
　（４）前記形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、解析対象領域の全体形状に対する影響度合を判定する工程、
とからなることを特徴とするプレス成形解析方法である。
　第二の発明は、（１）離型前のプレス成形品の形状、残留応力分布のデータを算出し、前記離型前のデータに基づいて座標変換により離型前の局所座標系での応力分布（ａ）を算出する工程、
　（２）前記離型前のデータに基づいてスプリングバック解析を行い、離型後のプレス成形品の残留応力分布を算出し、座標変換により局所座標系でのプレス成形品の残留応力分布（ｂ）を算出するとともに、離型後のプレス成形品の形状データ（ｄ）を算出する工程、
　（３）前記応力分布（ａ）と（ｂ）との差分（ａ−ｂ）を算出し、該差分をスプリングバッグ有効応力（ＳＢ有効応力）として座標変換により全体座標系におけるＳＢ有効応力分布を算出し、このＳＢ有効応力分布のうち解析対象領域のＳＢ有効応力を除去して、プレス成形品のＳＢ有効応力分布（Ｘ）を算出する工程、
　（４）前記局所座標系での残留応力分布（ｂ）を離型前のプレス成形品の形状に付与して局所座標系での残留応力相当分を算出し、座標変換により全体座標系でのプレス成形品の残留応力相当分（Ｙ）を算出する工程、
　（５）前記プレス成形品のＳＢ有効応力分布（Ｘ）に前記プレス成形品の残留応力相当分（Ｙ）を加算した応力分布を算出してスプリングバック解析を行い、離型後のプレス成形品の形状データ（ｃ）を算出する工程、
　（６）前記形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、解析対象領域の全体形状に対する影響度合を判定する工程、
とからなることを特徴とするプレス成形解析方法である。
　実際にスプリングバッグに寄与している応力よりも過剰または過少に評価することなく適切に対象の領域の応力のスプリングバッグに及ぼす影響を評価することができる。また、局所座標系を設定し、特定の応力成分を除去することにより、その応力の影響のみを評価することができ、スプリングバッグ低減対策を検討する際に役に立つ。
ＳＢ解析を行ったフロントピラーインナーの概観を示す図である。 フロントピラーインナーのＳＢ解析対象領域を説明する図である。 解析結果を示す図である。 第１の実施の形態に係るフロー図である。 第２の実施の形態に係るフロー図である。
　（第１の実施の形態）
　本発明の第１の実施の形態を図面を参照して説明する。
　図４は本発明に係るプレス成形解析方法の第１の実施の形態に係るプレス成形解析方法の処理の順序を示すフロー図である。
第１の実施の形態は、スプリングバック有効応力にのみ着目して、スプリングバック解析対象領域のスプリングバック有効応力を変更または除去を行った後にスプリングバック解析を行い、離型後のプレス成形品の形状を比較することで、スプリングバック解析対象領域の応力の影響を評価する解析方法である。
　図４に示すプレス成形解析方法における解析計算は、有限要素法（ＦＥＭ）解析システムを用いて行われる。
　（１）
　Ｓ１：先ず、ステップＳ１において離型前のプレス成形品の応力分布を求める。
　Ｓ２：ステップＳ２において、ステップＳ１で求められた全体座標系で表わされている応力を、離型前後で変化しない局部座標系を定義し、応力を座標変換し局部座標系での応力分布（ａ）を求める。プレス成形品は離型後にスプリングバックにより回転が発生するため、全体座標系で表わされている応力を直接差分すると変形に寄与した応力を正確に求めることができない。そこで、応力分布を与える各有限要素について、一般に有限要素解析内で用いられている手法である離型前後で変化しない局所座標系を定義し応力を座標変換した後に離型前後で応力の差分をおこなう。
　応力の座標変換は局所座標系の直交する３軸の単位ベクトルをそれぞれ（ｌ１，ｍ１，ｎ１）、（ｌ２，ｍ２，ｎ２）、（ｌ３，ｍ３，ｎ３）とすると式（１）で与えられる。
Figure JPOXMLDOC01-appb-M000001
　プレス成形シミュレーションでの局所座標系の採り方は一般に節点が４つまたは３つのシェル要素が用いられ板厚方向には節点を持たない。そこで、一般に板厚方向、即ち４つまたは３つの節点で作られる面法線方向を座標系の１つの第３軸に選ぶことが普通である。４節点シェル要素の場合、４点が同一平面上に必ずあるわけではないので、（１）対角線の外積から求められる方向を板厚方向とする、（２）近似平面の面法線方向を板厚方向とするなどの方法がある。また第１軸、第２軸は、（１）シェル要素の第１節点から第２節点への方向を第１軸とする方法、（２）対角線のベクトルから作る方法などがある。
　（２）
　Ｓ３：ステップＳ１の応力分布に基づきスプリングバック解析を行う。
　Ｓ４：離型後のプレス成形品の残留応力分布を求める。
　Ｓ５：全体座標系から局所座標系への座標変換を行い、離型後の局所座標系でのプレス成形品の残留応力分布（ｂ）を求める。
　Ｓ６：Ｓ３のスプリングバック解析から離型後のプレス成形品の形状データ（ｄ）を求める。
　（３）
　Ｓ７：離型前の応力分布（ａ）と離型後の残留応力分布（ｂ）との差分（ａ−ｂ）をスプリングバック有効応力（以下ＳＢ有効応力と呼ぶ）と呼び、その分布を求める。
　Ｓ８：局所座標系から全体座標系への座標変換を行い全体座標系でのＳＢ有効応力分布を求める。
　Ｓ９：解析対象領域のＳＢ有効応力を変更または除去する。
　Ｓ１０：解析対象領域のＳＢ有効応力の変更または除去後のプレス成形品のＳＢ応力分布を求める。
　Ｓ１１：ステップＳ１０で求めた応力分布からスプリングバック解析を行う。
　Ｓ１２：離型後のプレス成形品の形状データ（ｃ）を求める。
　（４）
　Ｓ１３：形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、ＳＢ解析対象領域の全体形状への影響度を評価する。
　（第２の実施の形態）
　本発明の第２の実施の形態を図面を参照して説明する。
　図５は本発明に係るプレス成形解析方法の第２の実施の形態に係るプレス成形解析方法の処理の順序を示すフロー図である。
　第２の実施の形態は、第１の実施の形態で求めた部分的にスプリングバック有効応力を変更または除去した局所座標系のスプリングバッグ有効応力分布に局所座標系の残留応力を加えた応力分布を全体座標系に変換したのちにスプリングバック解析を行い、離型後のプレス成形品の形状と比較することで、スプリングバック解析対象領域の応力の影響を評価する解析方法である。
　図５に示すプレス成形解析方法における解析計算は、有限要素法（ＦＥＭ）解析システムを用いて行われる。
　（１）
　Ｓ２１：先ず、ステップＳ２１において離型前のプレス成形品の応力分布を求める。
　Ｓ２２：ステップＳ２２において、ステップＳ２１で求められた全体座標系で表わされている応力を、離型前後で変化しない局所座標系を定義し、応力を座標変換し局所座標系での応力分布（ａ）を求める。局所座標系の求め方は第１の実施の形態に同じであるので本項では省略する。
　（２）
　Ｓ２３：ステップＳ２１の応力分布に基づきスプリングバック解析を行う。
　Ｓ２４：離型後のプレス成形品の残留応力分布を求める。
　Ｓ２５：全体座標系から局所座標系への座標変換を行い、離型後の局所座標系でのプレス成形品の残留応力分布（ｂ）を求める。
　Ｓ２６：ステップＳ２３のスプリングバック解析から離型後のプレス成形品の形状データ（ｄ）を求める。
　（３）
　Ｓ２７：離型前の応力分布（ａ）と離型後の残留応力分布（ｂ）との差分（ａ−ｂ）をスプリングバック有効応力（以下ＳＢ有効応力と呼ぶ）と呼び、その分布を求める。
　Ｓ２８：局所座標系から全体座標系への座標変換を行い全体座標系でのＳＢ有効応力分布を求める。
　Ｓ２９：解析対象領域のＳＢ有効応力を変更または除去する。
　Ｓ３０：解析対象領域のＳＢ有効応力の変更または除去後のプレス成形品のＳＢ有効応力分布（Ｘ）を求める。
　（４）
　Ｓ３１：離型前のプレス成形品の形状に応力（ｂ）を付与する。
　Ｓ３２：局所座標系での残留応力相当分を求める。
　Ｓ３３：局所座標系を全体座標系に変換し全体座標系でのプレス成形品の残留応力相当分（Ｙ）を求める。
　（５）
　Ｓ３４：ＳＢ有効応力分布（Ｘ）に残留応力相当分（Ｙ）を加算した応力分布を求める。
　Ｓ３５：スプリングバック解析を行う。
　Ｓ３６：離型後のプレス成形品の形状データ（ｃ）を求める。
　（６）
　Ｓ３７：形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、ＳＢ解析対象領域の全体形状への影響度を評価する。
　本発明に係るプレス成形解析方法の実施例として、フロントピラーインナー部品のプレス成形加工について述べる。
　本部品の概観形状を図１に示す。板厚１．６ｍｍの９８０ＭＰａ級高張力冷延鋼板を用いて金型によりプレス成形及びトリム加工を行い、ベンド加工を経てフロントピラーインナー部品に仕上げた。ベンド加工は、図１に示すベンドラインに沿って折り曲げ加工がなされている。
　プレス成形解析は、図１に示すフロントピラーインナー部品を成形する場合において、図２に示す解析対象領域１から６までの応力が、それぞれ、どの程度影響しているかを解析した。
　解析は、市販の有限要素法（ＦＥＭ）解析システムを用いて、図４及び図５に示すステップに沿って行った。具体的には、ベンド加工後において図２に示す各領域（１から６）の応力を除去して、点Ａ、点ＢでのＺ方向変位に及ぼす影響を調査した。なお、通常（各領域での応力除去を行っていないもの）のスプリングバック解析後の点ＡのＺ方向変位は、−１６．５ｍｍ、点Ｂは−９．７ｍｍである。従来例（特許文献１）は、解析対象領域（１から６の何れか１箇所）に関して全ての応力を除去してスプリングバック計算を行ったものである。発明例１は、本発明のＳＢ有効応力のみを取り出し（図４のＳ８）、更に解析対象領域（１から６の何れか１箇所）のみＳＢ有効応力を除去（図４のＳ９）してスプリングバック計算（図４のＳ１１）を行ったものである。発明例２は解析対象領域（１から６の何れか１箇所）のＳＢ有効応力のみを除去（図５のＳ２９）してスプリングバック計算（図５のＳ３５））を行ったものである。
　上記の結果を図３に示す。横軸は解析対象領域の番号を、縦軸は点ＡのＺ方向変位（ｍｍ）を表したものである。従来例では解析対象領域６の影響がもっとも大きいが、発明例１及び２では、解析対象領域２の影響が最も大きい。
　次に、要因分析結果を具体的に検証するために、解析対象領域２と解析対象領域６のそれぞれのフランジに△形状の切り込みを入れた場合について比較を行った。切り込みを入れることにより、その周囲の変形が他の部分に影響しなくなり、その結果、応力を除去した結果と同様の効果を実現できる為である。
　その結果、解析対象領域６に切り込みを入れた場合はＡ点が１．２ｍｍ上昇し、解析対象領域２の場合では、４．３ｍｍ上昇が見られ、解析対象領域２の影響が大きいことがわかった。このことから、発明例の要因分析が従来例に比べて優れていることがわかる。
　従来例において解析対象領域６の影響が大きくなったのは、解析対象領域６が比較的形状が複雑でスプリングバッグ後の残留応力が大きいため、その分を除去した影響が出ていると考えられる。また、ベンド工程では、解析対象領域４~解析対象領域６は塑性変形を加えるベンドラインよりも遠いことからも、従来例では解析対象領域６の影響を過剰に評価しているものと考えられる。
　１~６　　解析対象領域を表す符号である。

Claims (2)

  1. 　（１）離型前のプレス成形品の形状、残留応力分布のデータを算出し、前記離型前のデータに基づいて座標変換により離型前の局所座標系での応力分布（ａ）を算出する工程、
    　（２）前記離型前のデータに基づいてスプリングバック解析を行い、離型後のプレス成形品の残留応力分布を算出し、座標変換により局所座標系でのプレス成形品の残留応力分布（ｂ）を算出するとともに、離型後のプレス成形品の形状データ（ｄ）を算出する工程、
    　（３）前記応力分布（ａ）と（ｂ）との差分（ａ−ｂ）を算出し、該差分をスプリングバッグ有効応力（ＳＢ有効応力）として座標変換により全体座標系におけるＳＢ有効応力分布を算出し、このＳＢ有効応力分布のうち解析対象領域のＳＢ有効応力を変更または除去して、プレス成形品のＳＢ有効応力分布を算出し、スプリングバック解析を行い、離型後のプレス成形品の形状データ（ｃ）を算出する工程、
    　（４）前記形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、解析対象領域の全体形状に対する影響度合を判定する工程、
    とからなることを特徴とするプレス成形解析方法。
  2. 　（１）離型前のプレス成形品の形状、残留応力分布のデータを算出し、前記離型前のデータに基づいて座標変換により離型前の局所座標系での応力分布（ａ）を算出する工程、
    　（２）前記離型前のデータに基づいてスプリングバック解析を行い、離型後のプレス成形品の残留応力分布を算出し、座標変換により局所座標系でのプレス成形品の残留応力分布（ｂ）を算出するとともに、離型後のプレス成形品の形状データ（ｄ）を算出する工程、
    　（３）前記応力分布（ａ）と（ｂ）との差分（ａ−ｂ）を算出し、該差分をスプリングバッグ有効応力（ＳＢ有効応力）として座標変換により全体座標系におけるＳＢ有効応力分布を算出し、このＳＢ有効応力分布のうち解析対象領域のＳＢ有効応力を除去して、プレス成形品のＳＢ有効応力分布（Ｘ）を算出する工程、
    　（４）前記局所座標系での残留応力分布（ｂ）を離型前のプレス成形品の形状に付与して局所座標系での残留応力相当分を算出し、座標変換により全体座標系でのプレス成形品の残留応力相当分（Ｙ）を算出する工程、
    　（５）前記プレス成形品のＳＢ有効応力分布（Ｘ）に前記プレス成形品の残留応力相当分（Ｙ）を加算した応力分布を算出してスプリングバック解析を行い、離型後のプレス成形品の形状データ（ｃ）を算出する工程、
    　（６）前記形状データ（ｃ）と形状データ（ｄ）との差分（ｃ−ｄ）を算出し、解析対象領域の全体形状に対する影響度合を判定する工程、
    とからなることを特徴とするプレス成形解析方法。
PCT/JP2012/058774 2011-03-30 2012-03-26 プレス成形解析方法 WO2012133867A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP12763272.7A EP2692454B1 (en) 2011-03-30 2012-03-26 Press forming analysis method
MX2013011153A MX2013011153A (es) 2011-03-30 2012-03-26 Metodo para el analisis de conformacion por presion.
KR1020137026029A KR101544464B1 (ko) 2011-03-30 2012-03-26 프레스 성형 해석 방법
US14/007,974 US9410855B2 (en) 2011-03-30 2012-03-26 Method for press forming analysis
CN201280016493.7A CN103459058B (zh) 2011-03-30 2012-03-26 冲压成形分析方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011-075778 2011-03-30
JP2011075778A JP5765014B2 (ja) 2011-03-30 2011-03-30 プレス成形解析方法

Publications (1)

Publication Number Publication Date
WO2012133867A1 true WO2012133867A1 (ja) 2012-10-04

Family

ID=46931555

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/058774 WO2012133867A1 (ja) 2011-03-30 2012-03-26 プレス成形解析方法

Country Status (7)

Country Link
US (1) US9410855B2 (ja)
EP (1) EP2692454B1 (ja)
JP (1) JP5765014B2 (ja)
KR (1) KR101544464B1 (ja)
CN (1) CN103459058B (ja)
MX (1) MX2013011153A (ja)
WO (1) WO2012133867A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258095A (zh) * 2013-05-14 2013-08-21 西北工业大学 平底立铣刀通用铣削力建模方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106295055A (zh) * 2016-08-22 2017-01-04 南京埃斯顿自动化股份有限公司 一种对上模进深进行折弯回弹补偿的折弯机折弯方法
JP6519639B1 (ja) 2017-12-07 2019-05-29 Ｊｆｅスチール株式会社 スプリングバック量変動要因部位特定方法
JP6683269B1 (ja) 2019-02-01 2020-04-15 Ｊｆｅスチール株式会社 スプリングバック量変動要因部位特定方法
JP6852750B2 (ja) 2019-04-25 2021-03-31 Ｊｆｅスチール株式会社 スプリングバック量乖離要因部位特定方法および装置
JP6841295B2 (ja) 2019-05-22 2021-03-10 Ｊｆｅスチール株式会社 スプリングバック量乖離要因部位特定方法および装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007229724A (ja) 2006-02-27 2007-09-13 Jfe Steel Kk プレス成形解析方法
JP2008049389A (ja) * 2006-08-28 2008-03-06 Toyota Motor Corp 形状不良要因特定方法、装置及びプログラム
WO2008026777A1 (fr) * 2006-08-31 2008-03-06 Nippon Steel Corporation Procédé d'identification de cause d'occurrence de retour élastique, procédé d'affichage du degré d'influence du retour élastique, procédé d'identification de la partie à l'origine de l'occurrence de retour élastique, procédé de spécif

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4631335B2 (ja) * 2004-07-21 2011-02-16 Ｊｆｅスチール株式会社 プレス成形用金型の形状決定方法
EP1908539B1 (en) * 2005-06-30 2021-09-22 Nippon Steel Corporation Method and device for designing member, computer program, and computer-readable recording medium
JP4724626B2 (ja) * 2006-08-31 2011-07-13 新日本製鐵株式会社 スプリングバック発生原因部位特定方法、その装置、及びそのプログラム
JP4410833B2 (ja) * 2007-12-25 2010-02-03 新日本製鐵株式会社 スプリングバック発生原因分析方法、その装置、そのプログラム及び記録媒体
CN101339574B (zh) * 2008-08-12 2010-06-16 江苏大学 基于回弹补偿的混凝土搅拌叶片模具型面设计系统及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007229724A (ja) 2006-02-27 2007-09-13 Jfe Steel Kk プレス成形解析方法
JP2008049389A (ja) * 2006-08-28 2008-03-06 Toyota Motor Corp 形状不良要因特定方法、装置及びプログラム
WO2008026777A1 (fr) * 2006-08-31 2008-03-06 Nippon Steel Corporation Procédé d'identification de cause d'occurrence de retour élastique, procédé d'affichage du degré d'influence du retour élastique, procédé d'identification de la partie à l'origine de l'occurrence de retour élastique, procédé de spécif

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258095A (zh) * 2013-05-14 2013-08-21 西北工业大学 平底立铣刀通用铣削力建模方法

Also Published As

Publication number Publication date
JP2012206158A (ja) 2012-10-25
US20140019071A1 (en) 2014-01-16
EP2692454B1 (en) 2017-10-18
KR101544464B1 (ko) 2015-08-13
US9410855B2 (en) 2016-08-09
JP5765014B2 (ja) 2015-08-19
EP2692454A4 (en) 2014-10-08
MX2013011153A (es) 2013-11-01
CN103459058A (zh) 2013-12-18
EP2692454A1 (en) 2014-02-05
KR20130129457A (ko) 2013-11-28
CN103459058B (zh) 2015-11-25

Similar Documents

Publication Publication Date Title
WO2012133867A1 (ja) プレス成形解析方法
JP4894294B2 (ja) プレス成形解析方法
Firat et al. Sheet metal forming analyses with an emphasis on the springback deformation
Firat Computer aided analysis and design of sheet metal forming processes:: Part III: Stamping die-face design
US20100241366A1 (en) Shape defect factor identification method, device, and program
JP2009172677A (ja) スプリングバック発生原因分析方法、その装置、そのプログラム及び記録媒体
JP2006263788A (ja) 見込み金型形状設計システム
JP4352658B2 (ja) プレス成形品のスプリングバック解析方法
JP2005266894A (ja) 金型設計支援システム及び方法並びに金型設計支援用プログラム
JP6350498B2 (ja) プレスしわ発生判定方法
JP2011183417A (ja) スプリングバック安定性評価方法
JP5419284B2 (ja) プレス成型金型のスプリングバック見込み形状生成方法及び装置
US11731187B2 (en) Press forming method, rigidity-improvement-position specifying method, press forming system, and press-formed product
JP2007229761A (ja) プレス成形における被成形材の皺の発生の有無の予測方法、予測装置、予測プログラム、ならびにプレス成形方法
JP5151652B2 (ja) インバース法を用いた解析による特性解析方法
JP5757224B2 (ja) 構造部材の設計方法
JP2021126680A (ja) プレス成形品の形状変化予測方法
Thipprakmas et al. Wiping Z-bending die design for precise part fabrication
JP7342921B2 (ja) プレス成形品の形状変化予測方法
JP2014241129A (ja) スプリングバック量評価方法
JP6236297B2 (ja) 表面形状定量化方法
Safaei et al. Finite element analysis of influence of material anisotropy on the springback of advanced high strength steel
JP7484862B2 (ja) 金型の設計方法、金型の製造方法、及びプログラム
JP2012011458A (ja) プレス成形シミュレーションにおける割れ判定方法およびそれを用いたプレス成形部品の製造方法
Slota et al. Experimental and numerical analysis of springback prediction in U-bendings of anisotropic sheet metals

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12763272

Country of ref document: EP

Kind code of ref document: A1

REEP Request for entry into the european phase

Ref document number: 2012763272

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 14007974

Country of ref document: US

Ref document number: MX/A/2013/011153

Country of ref document: MX

Ref document number: 2012763272

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 20137026029

Country of ref document: KR

Kind code of ref document: A