JP4308239B2 - Springback influence display method, apparatus thereof, and program thereof - Google Patents

Springback influence display method, apparatus thereof, and program thereof Download PDF

Info

Publication number
JP4308239B2
JP4308239B2 JP2006268855A JP2006268855A JP4308239B2 JP 4308239 B2 JP4308239 B2 JP 4308239B2 JP 2006268855 A JP2006268855 A JP 2006268855A JP 2006268855 A JP2006268855 A JP 2006268855A JP 4308239 B2 JP4308239 B2 JP 4308239B2
Authority
JP
Japan
Prior art keywords
springback
molding
press
molded product
computer
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
JP2006268855A
Other languages
Japanese (ja)
Other versions
JP2008090481A (en
JP2008090481A5 (en
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
Priority to JP2006268855A priority Critical patent/JP4308239B2/en
Application filed by Nippon Steel Corp filed Critical Nippon Steel Corp
Priority to EP07806800.4A priority patent/EP2058062A4/en
Priority to TW096132255A priority patent/TW200824813A/en
Priority to CN200780032577.9A priority patent/CN101511502B/en
Priority to CA2662109A priority patent/CA2662109C/en
Priority to CN201210056854.2A priority patent/CN102708222B/en
Priority to CN201410204754.9A priority patent/CN103995927B/en
Priority to MX2009002150A priority patent/MX2009002150A/en
Priority to US12/439,126 priority patent/US9767234B2/en
Priority to BRPI0716142A priority patent/BRPI0716142B1/en
Priority to KR1020097004343A priority patent/KR101088115B1/en
Priority to PCT/JP2007/067358 priority patent/WO2008026777A1/en
Publication of JP2008090481A publication Critical patent/JP2008090481A/en
Publication of JP2008090481A5 publication Critical patent/JP2008090481A5/ja
Application granted granted Critical
Publication of JP4308239B2 publication Critical patent/JP4308239B2/en
Priority to US15/689,209 priority patent/US20180075177A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/80Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82Elements for improving aerodynamics

Description

本発明は、自動車用部材などを鋼板や棒鋼などの鋼材を始め、アルミ、マグネシウム、チタン等の金属からプレス成形した際にプレス成形製品に発生するスプリングバックの原因部位を特定するための方法、その装置、及びプログラムに関する。   The present invention is a method for identifying a cause part of a springback generated in a press-formed product when press-molding a member for an automobile such as a steel plate or a steel bar from a metal such as aluminum, magnesium, titanium, The present invention relates to the apparatus and the program.

ドアやバンパーなど多くの自動車用部材、家電部材、建材等は鋼板のプレス成形により製造されている。近年、それらの部材に対する軽量化の要求が高まっており、その要求を実現するために高強度を有する鋼材を使用することによって鋼材を薄手化するなどの対応が図られている。
しかしながら、鋼板の高強度化に伴い、プレス成形による製品形状の確保には厳しい管理が必要となっており、その管理において重要な項目の一つは、プレス成形時に鋼板に生じた残留応力を駆動力として、鋼板の弾性変形分が弾性回復する変形であるスプリングバックである。
Many automotive parts such as doors and bumpers, household electrical appliances, and building materials are manufactured by press forming of steel plates. In recent years, demands for weight reduction of these members have increased, and measures such as thinning the steel material by using a steel material having high strength in order to realize the demand have been attempted.
However, with the strengthening of steel sheets, strict management is required to secure the product shape by press forming, and one of the important items in the management is driving residual stress generated in the steel sheet during press forming. As a force, a spring back is a deformation in which the elastic deformation of the steel plate is elastically recovered.

特に、昨今は、自動車等の開発工数及びコスト削減のため、デザイン段階と同時に、その成形部材の成形方法を検討する設計段階が開始される傾向にあり、デザイン段階でのデザイン変更が、設計段階での成形部材の変更を生じさせるため、成形部材の成形方法を検討する設計段階での工数やコストは、自動車等の開発工程や開発費において、より大きな問題となっている。   In recent years, in order to reduce the development man-hours and costs of automobiles, etc., there is a tendency to start a design stage at the same time as the design stage to examine the molding method of the molded member. Therefore, the man-hours and costs at the design stage for examining the molding method of the molded member have become a bigger problem in the development process and development cost of automobiles and the like.

図1は、従来のスプリングバックへの対応を示す成形部材の断面模式図である。(a)は、成形製品の形状断面を示し、(b)は、(a)で示した成形製品と同型の金型で鋼板を冷間プレス成形後、成形製品に生じるスプリングバックを示し、(c)は、スプリングバックを想定して補正した金型の形状断面を示す。このように、(a)に示される成形製品を得るために、(c)に示すようなスプリングバックを予想した「見込み」の金型によって、所望の成形製品を得る対応がある。   FIG. 1 is a schematic cross-sectional view of a molded member showing correspondence to a conventional spring back. (A) shows the shape cross section of the molded product, (b) shows the spring back generated in the molded product after cold press-molding the steel plate with the same mold as the molded product shown in (a), ( c) shows the cross section of the mold shape corrected for spring back. As described above, in order to obtain the molded product shown in (a), there is a correspondence to obtain a desired molded product by using a “prospect” mold that anticipates spring back as shown in (c).

このようなスプリングバックを予想した金型を成形する方法として、有限要素法により金型にプレスされた下死点における鋼板の残留応力を解析し、その残留応力と反対向きの残留応力により生じる変形(スプリングフォワード)形状の金型を数値解析することにより、簡易にスプリングバックを考慮した金型を成形する方法が提案されている(特許文献1、非特許文献1)。   As a method of forming a mold that anticipates such springback, the residual stress of the steel plate at the bottom dead center pressed to the mold is analyzed by the finite element method, and deformation caused by the residual stress in the opposite direction to the residual stress There has been proposed a method of simply forming a mold in consideration of springback by numerically analyzing a (spring forward) mold (Patent Document 1, Non-Patent Document 1).

しかし、スプリングバックを完全に考慮した金型を数値解析により設計することは、非線形問題であり非常に困難であるため、提案されている方法は、あくまでも有限要素法によって簡易にスプリングバックを考慮した金型を成形するというものである。そのため、その金型によって、スプリングバックの許容量を満たさない場合にどのような対策が必要となるかについては、数値的に解析することが困難であるため、依然として、何の解決方法も提案されていない。   However, since it is very difficult to design a die that fully considers springback by numerical analysis, it is very difficult, so the proposed method simply considers springback by the finite element method. It is to mold a mold. For this reason, it is difficult to numerically analyze what measures are required when the mold does not satisfy the allowable amount of springback. Not.

そのため、簡易にスプリングバックを考慮した金型で、スプリングバックの許容量が満たされなかった場合、所望の成形製品を得るためにどのような対策をとるかは、技術者の経験に依存することになる。結局、その成形法による金型と実際の鋼板とによるトライアンドエラーテストが必要になる。   For this reason, a mold that easily considers springback, and if the allowable amount of springback is not satisfied, what measures are taken to obtain the desired molded product depends on the experience of engineers. become. Eventually, a trial-and-error test using a mold by the forming method and an actual steel plate is required.

また、金型形状ではなく、鋼材や成形製品の形状に残留応力を除去する修正を加えることでスプリングバックを少なくする方法が提案されている。
図2は、スプリングバックによる変形の発生原因となる部位を探る従来の方法を例示する斜視図である。(a)は、成形製品の形状を示し、(b)は、製品の一部1を切断除去した場合を示し、(c)は、製品に孔2を開けた場合を示し、(d)は、製品の一部にスリット3を付与する場合を示す。このような対策を講じることによってスプリングバックの挙動を観察すると共に、スプリングバックを少なくする対策が試みられている。
In addition, a method has been proposed in which the springback is reduced by modifying the shape of the steel material or molded product, not the mold shape, to remove residual stress.
FIG. 2 is a perspective view illustrating a conventional method for searching for a site that causes deformation due to springback. (A) shows the shape of the molded product, (b) shows the case where part 1 of the product is cut and removed, (c) shows the case where a hole 2 is made in the product, (d) shows The case where the slit 3 is given to a part of the product is shown. By taking such measures, an attempt has been made to observe the behavior of the springback and to reduce the springback.

また、スプリングバック発生部位への対策により、スプリングバックの原因となる残留応力は低減されても、切断除去や孔開け等のために部材自体の剛性が低下するので、わずかな残留応力でスプリングバックが発生してしまうという問題が生じ、根本的な原因究明には至らない。さらに、このような対策は、実際に試験金型と鋼板とによるテストを必要とするため、設計段階の工数とコストの増大という問題が生じる。   Even if the residual stress that causes the springback is reduced by measures against the springback occurrence site, the rigidity of the member itself decreases due to cutting and drilling, etc. This causes the problem that the problem will occur, and the root cause cannot be investigated. Furthermore, since such a measure actually requires a test using a test mold and a steel plate, there arises a problem of an increase in man-hours and costs at the design stage.

特開2003−33828号公報JP 2003-33828 A 三菱自動車テクニカルレビュー(2006 No.18 126〜131頁)Mitsubishi Motors Technical Review (2006 No.18, pages 126-131)

上述のような問題点に鑑み、本発明は、数値解析によりプレス成形製品のスプリングバック発生の原因となる部位を効率的にかつ経済的に特定し、それにより、成形部材の成形方法の検討時間を効率的にかつ経済的に短縮するために、スプリングバックの影響度の表示方法を提供することを課題とする。   In view of the above-described problems, the present invention efficiently and economically identifies a part that causes the occurrence of springback of a press-molded product by numerical analysis, thereby enabling a time for studying a molding member molding method. It is an object to provide a method for displaying the degree of influence of springback in order to efficiently and economically shorten the above.

また、本発明は、製品強度を維持しつつ、スプリングバック量を許容値以下とする成形製品を効率的にかつ経済的に特定するために、スプリングバックの影響度の表示方法を提供することを課題とする。   In addition, the present invention provides a method for displaying the degree of influence of springback in order to efficiently and economically identify a molded product having a springback amount that is less than or equal to an allowable value while maintaining product strength. Let it be an issue.

以上の課題を解決するために、本発明は、スプリングバック影響度を表示するコンピュータにおいて、プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析ステップ(手順:以下、カッコ内はプログラムの発明に対応する)と、プレス成形品の成形データのうち、プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理ステップ(手順)と、演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップ(手順)と、各領域毎に、算出したスプリングバック量を表示する表示ステップ(手順)とを有するスプリングバック影響度表示方法(プログラム)を提供する。   In order to solve the above-described problems, the present invention provides a press molding analysis step (procedure: the following) that obtains molding data of a press-molded product by numerically analyzing press molding molding conditions in a computer that displays a springback influence degree. (In parentheses correspond to the invention of the program) and among the molding data of the press-molded product, computation processing is performed on at least one of the physical property value and physical quantity data of each area into which the press-molded product is divided A step (procedure), a springback amount calculating step (procedure) for calculating the springback amount based on the result of the arithmetic processing, and a display step (procedure) for displaying the calculated springback amount for each region. Provided is a springback influence display method (program).

演算処理ステップ(手順)は、プレス成形品の成形データのうち、プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行うことが出来る。また、物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値としても良い。   The arithmetic processing step (procedure) performs arithmetic processing on at least one of physical property values and physical quantity data at some integration points in each region obtained by dividing the press-molded product among the molding data of the press-molded product. I can do it. The physical property value and physical quantity may be a plate thickness, an elastic coefficient, a plastic coefficient, a stress component value, or a strain component value.

表示ステップ(手順)は、各領域毎に、出したスプリングバック量を各領域の面積で割った値をコンタ表示するまた、表示ステップ(手順)は、各領域毎に、算出したスプリングバック量をプレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値を表示することができる。さらに、表示ステップ(手順)は、各領域毎に、算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を表示することもできる。 Display step (Step), for each region, a value obtained by dividing the amount of spring back which issued calculated by the area of each region contour display. In addition, in the display step (procedure), for each region, the calculated springback amount is divided by any one of the representative length, representative width, representative height, representative plate thickness, and tensile strength of the press-formed product. The value can be displayed. Furthermore, the display step (procedure) can also display a value obtained by dividing the calculated springback amount by the punch speed or the wrinkle pressing force for each region.

また、本発明は、プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析部と、成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、プレス成形品の成形データのうち、プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行ない、スプリングバック解析部に、演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、各領域毎に算出したスプリングバック量を表示する表示部とを有するスプリングバック影響度表示装置を提供する。   The present invention also includes a press molding analysis unit that numerically analyzes molding conditions for press molding to obtain molding data of a press molded product, and a spring back analysis unit that performs numerical analysis of the molding data and calculates a springback amount. The calculation processing is performed on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product, and the spring based on the result of the calculation processing is performed in the springback analysis unit. Provided is a springback influence display device having an arithmetic processing unit that calculates a back amount and a display unit that displays a springback amount calculated for each region.

演算処理部は、プレス成形品の成形データのうち、プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行なうことができ、物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値としても良い。   The arithmetic processing unit can perform arithmetic processing on at least one of the physical property value and physical quantity data of a part of the integration points of each region obtained by dividing the press molded product among the molding data of the press molded product, The physical property value and physical quantity may be a plate thickness, an elastic coefficient, a plastic coefficient, a stress component value, or a strain component value.

表示部は、各領域毎に、出したスプリングバック量を各領域の面積で割った値を、表示部にコンタ表示る。表示部は、各領域毎に、算出したスプリングバック量をプレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値を、表示部に表示することもできる。さらに、表示部は、各領域毎に、算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を、表示することもできる。 Display unit, for each region, a value obtained by dividing the amount of spring back which issued calculated by the area of each region, you Conta displayed on the display unit. The display unit displays the value obtained by dividing the calculated springback amount for each region by one of the representative length, representative width, representative height, representative plate thickness, and tensile strength of the press-formed product. Can also be displayed. Further, the display unit can display a value obtained by dividing the calculated springback amount by the punching speed or the wrinkle pressing force for each region.

本発明のプログラムは、上記に説明したプログラムの他に、プレス成形の成形条件を数値解析してプレス成形品の成形データを算出するプレス成形解析プログラムと、該成形データを数値解析してスプリングバック量を算出するスプリングバック解析プログラムと、該スプリングバック量をコンタ表示するポスト処理プログラムは市販のプログラムを使用し、これらのプログラムと、データ入出力が可能であるスプリングバック影響度表示プログラムであって、プレス成形解析プログラムから、成形データを取得する手順と、成形データのうち、プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、スプリングバック解析プログラムに、前記演算処理の結果を出力する手順と、スプリングバック解析プログラムが各領域毎に前記算出したスプリングバック量を、前記ポスト処理プログラムに、コンタ表示させるコンタ表示手順をコンピュータに実行させるものでも良い。市販のスプリングバック解析プログラムは、本発明のプログラムから出力された演算処理データに基づいてスプリングバック量を算出する。   In addition to the program described above, the program of the present invention includes a press molding analysis program for numerically analyzing press molding molding conditions to calculate molding data of a press molded product, and a numerical analysis of the molding data for springback. The springback analysis program for calculating the amount and the post-processing program for contour-displaying the amount of springback use commercially available programs, and these programs and the springback influence display program capable of data input / output. , A procedure for obtaining molding data from the press molding analysis program, and an arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region into which the press molded product is divided, among the molding data, A procedure for outputting the result of the arithmetic processing to a springback analysis program; The spring back amount of springback analysis program is the calculated for each of the regions, in the post-processing program may be one which perform the contour display procedure for contour displayed on the computer. A commercially available springback analysis program calculates the amount of springback based on the arithmetic processing data output from the program of the present invention.

本発明によれば、プレス成形製品のスプリングバック発生の原因として解析対象となる特定部位を、プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理し、その演算処理の結果に基づいて、スプリングバック量を各領域毎に、コンタ表示することにより、コンピュータ上での解析的なトライアンドエラーを可能とし、スプリングバック発生の原因部位の特定を経済的にかつ効率的に短縮し、かつ、成形製品の成形方法の検討時間を経済的にかつ効率的に短縮する方法を提供する。   According to the present invention, the specific part to be analyzed as a cause of the occurrence of springback of the press-molded product is calculated with respect to at least one of the physical property value and physical quantity data of each region obtained by dividing the press-molded product, Based on the result of the calculation process, the amount of springback is displayed for each region in a contoured manner, enabling analytical trial and error on a computer, and identifying the cause of the occurrence of springback economically. The present invention also provides a method for efficiently and efficiently reducing the examination time of a molding method of a molded product economically and efficiently.

さらに、本発明は、特定部位及びその特定部位の物性値及び物理量の少なくとも何れかを変化させることで、実際の成形製品に対して切断除去や孔開けを行わずに、スプリングバック変化量を許容値以下とする特定部位とその性状データを求めることが可能である。そのため、このように解析された成形製品は、スプリングバック変化量を許容値以下とすると同時に、製品の剛性を維持することが確認できるため、成形製品に求められる所定のスプリングバック変化量を維持するために、剛性を犠牲にするような成形製品の切断除去や孔開け等の対策を不要とすることができる。   Furthermore, the present invention allows the amount of springback change without cutting and punching the actual molded product by changing the specific part and / or the physical property value and physical quantity of the specific part. It is possible to obtain a specific part and its property data that are less than or equal to the value. Therefore, it can be confirmed that the molded product analyzed in this way maintains the rigidity of the product at the same time as the amount of change in the springback is less than the allowable value, and therefore the predetermined amount of change in springback required for the molded product is maintained. For this reason, it is possible to eliminate the need for measures such as cutting and removing the molded product and punching at the expense of rigidity.

図3は、本発明の一実施形態によるスプリングバック発生原因部位特定装置の機能構成図を示す。
スプリングバック発生原因部位特定装置10Aは、プレス成形解析部11、スプリングバック解析部12、物性値・物理量演算処理部14、ファイル格納部16、成形条件入力部18、及びスプリングバック量出力画面19を有する。
FIG. 3 is a functional configuration diagram of the springback occurrence cause identifying apparatus according to the embodiment of the present invention.
The springback generation cause identifying apparatus 10A includes a press molding analysis unit 11, a springback analysis unit 12, a physical property value / physical quantity calculation processing unit 14, a file storage unit 16, a molding condition input unit 18, and a springback amount output screen 19. Have.

成形条件入力部18は、後述するプレス成形解析部11及びスプリングバック解析部12で解析対象となる鋼板の形状データ(板厚、長さ、幅、曲率、歪等)、性状(強度、伸びなどの材質、板厚などの形状)、金型形状(ダイ(金型)及びパンチ形状、曲率、径、クリアランス、潤滑条件)、プレス条件(しわ押さえ荷重、パッド荷重、ビード張力、プレス圧力、温度)などの成形条件を入力する入力部である。または、成形解析におけるデータ領域、物性値・物理量演算処理部14におけるデータ領域、スプリングバック量出力画面19で画面表示する際の分割領域等を別個に設定することも可能である。   The forming condition input unit 18 includes shape data (plate thickness, length, width, curvature, strain, etc.) and properties (strength, elongation, etc.) of a steel plate to be analyzed by the press forming analysis unit 11 and the springback analysis unit 12 described later. Material, shape such as plate thickness), mold shape (die (die) and punch shape, curvature, diameter, clearance, lubrication condition), press condition (wrinkle holding load, pad load, bead tension, press pressure, temperature) ) And other molding conditions. Alternatively, a data area in the molding analysis, a data area in the physical property value / physical quantity calculation processing unit 14, a divided area when the screen is displayed on the springback amount output screen 19, and the like can be set separately.

プレス成形解析部11は、成形条件入力部18で入力された成形条件を入力情報として、弾塑性有限要素法、剛塑性有限要素法、ワンステップ有限要素法、境界要素法、初等解析等を用いて、プレス成形された被加工物である鋼板等の応力や歪や板厚などを求める数値解析を行う。プレス成形解析部11は、被加工物の板厚、応力の成分値、歪の成分値等の変数や、その変数の分布という形式で数値解析結果を出力する。その出力データは、例えば、ファイル「P org.k」として、スプリングバック解析部12、後述する物性値・物理量演算処理部14及びファイル格納部16に出力される。   The press molding analysis unit 11 uses an elastic-plastic finite element method, a rigid-plastic finite element method, a one-step finite element method, a boundary element method, an elementary analysis, etc., with the molding conditions input by the molding condition input unit 18 as input information. Then, numerical analysis is performed to obtain stress, strain, plate thickness and the like of a press-formed workpiece such as a steel plate. The press forming analysis unit 11 outputs a numerical analysis result in the form of variables such as a plate thickness of the workpiece, a stress component value, a strain component value, and a distribution of the variables. The output data is output, for example, as a file “P org.k” to the springback analysis unit 12, a physical property value / physical quantity calculation processing unit 14 and a file storage unit 16 described later.

このプレス成形解析部11による数値解析は、有限要素法(例えば市販のソフトウェアPAM-STAMP、LS-DYNA、Autoform、OPTRIS、ITAS-3D、ASU/P-FORM、ABAQUS、ANSYS、MARC、HYSTAMP、Hyperform、SIMEX、Fastform3D、Quikstamp)を用いて、プレス成形する製品の形状データ(板厚、長さ、幅、曲率、歪みなど)及び、使用する金属板の性状(強度、伸びなどの材質、板厚などの形状)に基づいて、必要であれば金型形状(ダイ及びパンチ形状、曲率、径、クリアランス、潤滑条件)、プレス条件(温度、圧力)などの成形条件を設定し、成形解析を行い、かつ、成形後の応力、歪値の分布を数値的に得ることができる。
また、プレス成形解析部11は、有限要素法による応力分布や曲率等の解析結果を、ポスト処理ソフトを用いて、スプリングバック量出力画面19にコンタ表示することが可能である。
Numerical analysis by the press forming analysis unit 11 is performed by a finite element method (for example, commercially available software PAM-STAMP, LS-DYNA, Autoform, OPTRIS, ITAS-3D, ASU / P-FORM, ABAQUS, ANSYS, MARC, HYSTAMP, Hyperform. , SIMEX, Fastform3D, Quikstamp) shape data (plate thickness, length, width, curvature, strain, etc.) of the product to be press-molded and properties of the metal plate used (strength, elongation, etc., material thickness) If necessary, set molding conditions such as die shape (die and punch shape, curvature, diameter, clearance, lubrication conditions), press conditions (temperature, pressure), etc., and perform molding analysis. In addition, the distribution of stress and strain values after molding can be obtained numerically.
Further, the press forming analysis unit 11 can display the analysis results such as stress distribution and curvature by the finite element method on the springback amount output screen 19 using post processing software.

スプリングバック解析部12は、プレス成形解析部11の出力データファイル「P org.k」及び後述する物性値・物理量演算処理部14の出力データファイル「P rem.casen.k」を入力データとして用いて、スプリングバック解析を行う。スプリングバック解析とは、プレス成形解析部11の出力結果である被加工物の板厚、応力の成分値、歪の成分値等の変数、及び変数分布に基づいて、弾性有限要素法,弾塑性有限要素法、ワンステップ有限要素法、初等解析等により、除荷過程の計算を行い、被加工物に生じる変形量であるスプリングバック量を数値解析することである。そのスプリングバック量は、被加工物を有限要素法等により分割し、3次元データ座標の要素毎に計算される。なお、被加工物に生じる変形量であるスプリングバック量としては、被加工物の任意点の変形量、被加工物の指定領域内の最大変位点もしくは最小変位点の変形量、被加工物内の複数の任意の面又は線同士がなす角度、または、被加工物の任意の面又は線のなす曲率等がある。   The springback analysis unit 12 uses an output data file “P org.k” of the press forming analysis unit 11 and an output data file “P rem.casen.k” of a physical property value / physical quantity calculation processing unit 14 described later as input data. To perform springback analysis. The springback analysis is an elastic finite element method, elastoplasticity, based on variables such as the plate thickness of the workpiece, the stress component value, the strain component value, and the variable distribution, which are output results of the press forming analysis unit 11. The unloading process is calculated by a finite element method, a one-step finite element method, an elementary analysis, and the like, and a springback amount that is a deformation amount generated in a workpiece is numerically analyzed. The springback amount is calculated for each element of the three-dimensional data coordinates by dividing the workpiece by the finite element method or the like. The amount of springback, which is the amount of deformation that occurs in the workpiece, includes the amount of deformation at an arbitrary point of the workpiece, the amount of deformation at the maximum or minimum displacement point within the specified area of the workpiece, There are an angle formed by a plurality of arbitrary surfaces or lines, or a curvature formed by an arbitrary surface or line of the workpiece.

スプリングバック解析部12のスプリングバック解析結果である出力データは、スプリングバック量、スプリングバック時の歪等の形状、性状、応力等を含み、スプリングバック量出力画面19に出力され、かつ、入力データファイル「P org.k」による数値解析結果出力データファイル「SB org.k」、又は、「P rem.casen.k」による数値解析結果出力データファイル「SB rem.casen.k」として、スプリングバック解析部12、後述する物性値・物理量演算処理部14及びファイル格納部16に出力される。   The output data that is the result of the springback analysis by the springback analysis unit 12 includes the springback amount, shape such as distortion at the time of springback, properties, stress, etc., and is output to the springback amount output screen 19 and input data. Springback as numerical analysis result output data file “SB org.k” by file “P org.k” or numerical analysis result output data file “SB rem.casen.k” by “Prem.casen.k” The data is output to the analysis unit 12, a physical property value / physical quantity calculation processing unit 14 and a file storage unit 16 which will be described later.

物性値・物理量演算処理部14は、プレス成形解析部11の出力データファイル「P org.k」、スプリングバック解析部12の数値解析結果出力データファイル「SB org.k」、又は、「SB rem.casen.k」を入力し、演算処理を行い、その結果として「P rem.casen.k」をスプリングバック解析部12に出力する。なお、ここで言う演算処理は、上記変数の少なくとも1個以上を係数倍すること、ゼロを含む一定値にすること、四則演算をすること、関数に基づいて計算すること、一定ではない任意の値に置換することをいう。物性値・物理量演算処理部14は、このような、演算処理により、スプリングバック量を減少させる変数及び特定部位を判別することを可能にする。   The physical property value / physical quantity calculation processing unit 14 outputs the output data file “P org.k” of the press forming analysis unit 11, the numerical analysis result output data file “SB org.k” of the springback analysis unit 12, or “SB rem”. .Casen.k ”is input, the calculation process is performed, and“ P rem.casen.k ”is output to the springback analysis unit 12 as a result. Note that the arithmetic processing referred to here includes multiplying at least one of the above variables by a coefficient, making a constant value including zero, performing four arithmetic operations, calculating based on a function, and any non-constant arbitrary It means replacing with a value. The physical property value / physical quantity calculation processing unit 14 makes it possible to determine a variable and a specific part that reduce the springback amount by such calculation processing.

物性値・物理量演算処理部14は、入力データから被加工体の1つであるプレス成形製品のデータを取得すると、そのプレス成形製品の一部から特定部位を発見するため、そのプレス成形製品のデータを複数の領域に分割する処理を行う必要がある。その領域分割方法の1つは、均等な寸法により被加工体の領域を分割することである。   When the physical property value / physical quantity calculation processing unit 14 acquires data of a press-molded product that is one of the workpieces from the input data, the physical property value / physical quantity calculation processing unit 14 finds a specific part from the part of the press-molded product. It is necessary to divide the data into a plurality of areas. One of the region dividing methods is to divide the region of the workpiece by uniform dimensions.

均等な寸法は、特定部位に対する対策を考慮することが望ましい。つまり、特定部位を発見しても、特定部位への有効な対策がとれない極めて微小な領域であっては、解析結果が有効に活かせない場合がある。例えば、特定部位への対策が20mm×100mmのビードを加えることである場合、分割領域は、20mm角となることが挙げられる。   For uniform dimensions, it is desirable to consider measures for specific parts. That is, even if a specific part is found, the analysis result may not be used effectively in a very small area where an effective measure cannot be taken for the specific part. For example, when the countermeasure for a specific part is to add a bead of 20 mm × 100 mm, the divided area may be 20 mm square.

成形品の分割領域を決定するための1つの方法として、曲率や要素に基づいて、分割領域を決定する方法もある。
プレス成形解析部11は、数値解析により、成形解析後の被加工物の幾何学的な情報として各節点の座標値を得、各要素の曲率や要素間の角度を客観的に算出することが可能である。プレス成形品の場合、変形後の被加工物はパンチ肩Rやダイ肩R等の曲げの稜線部では、曲げ稜線に垂直な方向で他の部位に比べ非常に大きな曲率となる。
ここでは各要素の最大曲率を計算し、ある閾値以上の曲率となる要素を非表示にすれば、パンチ肩Rやダイ肩R以外の部位(ウェブ面、縦壁部、フランジ部)を連続していない別々の領域として判別分離することが可能である。
As one method for determining a divided area of a molded product, there is a method of determining a divided area based on curvature and elements.
The press forming analysis unit 11 can obtain the coordinate value of each node as geometric information of the workpiece after the forming analysis by numerical analysis, and objectively calculate the curvature of each element and the angle between the elements. Is possible. In the case of a press-formed product, the deformed workpiece has a very large curvature in the bending ridge line portion such as the punch shoulder R and the die shoulder R in the direction perpendicular to the bending ridge line compared to other parts.
Here, if the maximum curvature of each element is calculated and the element having a curvature greater than a certain threshold value is hidden, the parts other than the punch shoulder R and die shoulder R (web surface, vertical wall portion, flange portion) are continuous. It is possible to discriminate and separate as separate areas that are not.

図4〜図8にプレス成形解析から演算処理までの工程の例を示す。
図4がプレス成形解析結果であり、応力分布がコンタ図で表示されている。プレス成形解析結果である各要素の曲率を表示したものが図5であり、曲率の大きい要素を非表示としたのが図6である。段差部で非常に大きな曲率となることから、縦壁部内である閾値以上の曲率となる要素を非表示にすれば、段差の稜線部を境界とした複数の領域に分離可能であり、個別の領域を抽出または選択することで特徴形状に基づいた領域を選択可能である。図6により、段差のある側壁は3つの面と、ウェブ、段差のない縦壁部(図から隠れている)、両側の各フランジ部に分割できることがわかる。図6により分割したそれぞれの領域について、図4に示した成形データ(応力分布)に基づいて、応力の高い領域について演算処理を行うと判定したのが、図7に示す領域A15〜A17である。判定工程では、曲率及び/又は角度に基づいて分割した領域について、後述の演算を行うかを成形データに基づいて判定する。例えば、分割した各領域のプレス成形解析後の応力レベルにより判定する。
4 to 8 show examples of steps from press forming analysis to calculation processing.
FIG. 4 shows the press forming analysis result, and the stress distribution is displayed in a contour diagram. FIG. 5 shows the curvature of each element as a result of the press forming analysis, and FIG. 6 shows the element with a large curvature hidden. Since the stepped portion has a very large curvature, it can be separated into a plurality of regions with the ridgeline portion of the step as a boundary by hiding elements that have a curvature greater than the threshold value in the vertical wall portion, A region based on the feature shape can be selected by extracting or selecting the region. From FIG. 6, it can be seen that the stepped side wall can be divided into three surfaces, a web, a vertical wall portion with no step (hidden from the figure), and each flange portion on both sides. The areas A15 to A17 shown in FIG. 7 are determined to be subjected to the arithmetic processing for the high stress areas based on the molding data (stress distribution) shown in FIG. . In the determination step, it is determined on the basis of the molding data whether the calculation described later is performed on the region divided based on the curvature and / or the angle. For example, the determination is made based on the stress level after press molding analysis of each divided area.

領域A15〜A17について、それぞれの領域の全ての要素に対して、例えば応力をゼロとする演算処理を行った状態が図8である。
ここでは要素の最大曲率を用いたが要素間の角度を用いても同様の領域分割を行うことができる。
また、段差のないハット断面形状のメンバー部材を考えると,分離された各領域として、ウェブ面、両側の縦壁部、両側のフランジ部を選択することができる。また、選択した領域が平面でない場合、同様の手法を繰り返すことで、さらに特徴的な領域を選択していくことが可能である。
このようにして、物性値・物理量演算処理部14は、プレス成形解析部11の解析結果より、曲率や要素に基づいて、分割領域を決定することが可能である。
For the areas A15 to A17, FIG. 8 shows a state in which, for example, a calculation process for setting the stress to zero is performed on all the elements in the respective areas.
Although the maximum curvature of the element is used here, the same region division can be performed using the angle between the elements.
Further, when considering a member member having a hat cross-sectional shape without a step, a web surface, vertical wall portions on both sides, and flange portions on both sides can be selected as separated regions. Further, when the selected region is not a plane, it is possible to select a more characteristic region by repeating the same method.
In this way, the physical property value / physical quantity calculation processing unit 14 can determine the divided region based on the curvature and the elements based on the analysis result of the press forming analysis unit 11.

スプリングバック解析部12は、有限要素法(例えば市販のソフトウェアPAM-STAMP、LS-DYNA、Autoform、OPTRIS、ITAS-3D、ASU/P-FORM、ABAQUS、ANSYS、MARC、HYSTAMP)を用いて、プレス成形解析部11で得られた「P org.k」内に記載の応力分布を入力し、スプリングバック解析を実施する。ソフトウェア内でのスプリングバックの計算は例えば、有限弾塑性変形の基礎式や離散化手法に沿った内容で計算される。ただし、スプリングバックの計算は弾性解析であっても、弾塑性解析であっても良い。   The springback analysis unit 12 uses a finite element method (for example, commercially available software PAM-STAMP, LS-DYNA, Autoform, OPTRIS, ITAS-3D, ASU / P-FORM, ABAQUS, ANSYS, MARC, HYSTAMP). The stress distribution described in “P org.k” obtained by the forming analysis unit 11 is input, and the springback analysis is performed. The calculation of the springback in the software is, for example, calculated based on the basic equation of finite elasto-plastic deformation and the content according to the discretization method. However, the springback calculation may be an elastic analysis or an elasto-plastic analysis.

また、スプリングバック解析部12は、有限要素法によるスプリングバック解析結果を、ポスト処理ソフトを用いて、スプリングバック量出力画面19にコンタ表示することが可能である。また、物性値・物理量演算処理部14において、スプリングバック量を、分割領域面積、プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さ等で除算することにより規格化し、除算した変数とスプリングバックの影響度をより判別し易くしてデータ出力し、その出力データを用いて、ポスト処理ソフトにより、コンタ表示することが可能である。このような規格化により、被加工体の物性値・物理量とスプリングバック量との関係を視覚的に捉え易くすることができる。   Further, the springback analysis unit 12 can display the result of the springback analysis by the finite element method on the springback amount output screen 19 using post processing software. Further, the physical property value / physical quantity calculation processing unit 14 standardizes by dividing the spring back amount by the divided area, the representative length of the press-formed product, the representative width, the representative height, the representative plate thickness, the tensile strength, and the like. It is possible to output the data by making it easier to discriminate between the divided variable and the influence degree of the springback, and use the output data to display the contour by post processing software. Such normalization makes it easy to visually grasp the relationship between the physical property value / physical quantity of the workpiece and the springback amount.

さらに、物性値・物理量演算処理部14において、スプリングバック量を、パンチ速度又はしわ押さえ力等のプレス成形条件で除算することで規格化した値を求め、ポスト処理ソフトにより、コンタ表示することで、プレス成形条件と、スプリングバック量との関係を視覚的に捉え易くすることもできる。このようなスプリングバック量の規格化とその数値のコンタ表示は、スプリングバック発生の原因部位の特定を経済的にかつ効率的に短縮し、かつ、成形製品の成形方法の検討時間を経済的にかつ効率的に短縮することが可能とする。   Furthermore, the physical property value / physical quantity calculation processing unit 14 obtains a normalized value by dividing the springback amount by a press molding condition such as a punch speed or a wrinkle holding force, and displays the contour by post processing software. The relationship between the press molding conditions and the amount of springback can be easily grasped visually. This standardization of the amount of springback and the contour display of the numerical value shorten the identification of the cause of the occurrence of springback economically and efficiently, and economically consider the molding method of molded products. And it is possible to shorten it efficiently.

ファイル格納部16は、プレス成形解析部11、スプリングバック解析部12、後述する物性値・物理量演算処理部14の出力結果であるデータファイル「P org.k」、「SB org.k」、「P rem.casen.k」、「SB rem.casen.k」、「P trim.casen.k」等を保存するための格納部である。しかし、それらのデータファイル等が、プレス成形解析部11と、スプリングバック解析部12と、物性値・物理量演算部14との間で直接入出力される場合は、このファイル格納部16は、必ずしも必要では無い。   The file storage unit 16 includes data files “P org.k”, “SB org.k”, “SB”, which are output results of the press forming analysis unit 11, the springback analysis unit 12, and a physical property value / physical quantity calculation processing unit 14 described later. A storage unit for storing “P rem.casen.k”, “SB rem.casen.k”, “P trim.casen.k”, and the like. However, when these data files are directly input / output between the press forming analysis unit 11, the springback analysis unit 12, and the physical property value / physical quantity calculation unit 14, the file storage unit 16 is not necessarily provided. Not necessary.

物性値・物理量演算処理部14は、プレス成形解析部11の出力データファイル「P org.k」内の一部の領域を選択し、その一部の領域の板厚、弾性係数、塑性係数、応力の成分値、歪の成分値の少なくとも一つ以上の変数に対して演算処理を実施する。そして、演算処理後の出力データファイル「P rem.case1.k」を生成して、スプリングバック解析部12に出力する。このデータファイルの送受信は、繰り返し行うことが可能であり、ファイル格納部16に逐次保存可能であり、また、ファイル形式ではなく実行プロセス又はスレッドのデータ入出力という形でデータのやりとりを行っても良い。   The physical property value / physical quantity calculation processing unit 14 selects a partial region in the output data file “P org.k” of the press forming analysis unit 11, and the plate thickness, elastic modulus, plasticity coefficient, An arithmetic process is performed on at least one of the stress component value and the strain component value. Then, an output data file “P rem.case1.k” after the arithmetic processing is generated and output to the springback analysis unit 12. This data file transmission / reception can be repeated, can be stored in the file storage unit 16 sequentially, and can be exchanged in the form of data input / output of the execution process or thread instead of the file format. good.

有限要素法で使用されるスプリングバック特定部位又は原因判断のための領域データは、有限要素法における1個以上の要素、複数の要素からなる微小領域、又は、要素が連続した集合体などの計算単位区分、有限要素法による数値積分手法の計算単位で応力や歪を計算する点のうちの一部である積分点(板厚方向における一部の点、面内の一分の点の双方を含む)等とすることができる。「1個以上の要素」とは有限要素法の領域分割単位である要素をいい、「計算単位区分」とは初等解析における計算単位となる微小領域の一つまたは連続した集合体をいう。また、「積分点」とは、一般的に有限要素法で行う近似的な積分を行う点である。「一部の積分点」とは、有限要素法による数値積分手法の計算単位で応力や歪を計算する点のうちの一部をいい、板厚方向における一部の点、及び面内の一部の点の双方を含むものとする。   The region data for determining the springback specific part or cause used in the finite element method is the calculation of one or more elements in the finite element method, a micro area consisting of a plurality of elements, or a collection of continuous elements. Integral points that are part of points for calculating stress and strain in the unit of calculation, unit of numerical integration by the finite element method (some points in the thickness direction, both points in the plane) Including). “One or more elements” refers to an element that is a region division unit of the finite element method, and “calculation unit classification” refers to one or a continuous set of minute regions that are calculation units in primary analysis. The “integration point” is a point at which approximate integration is generally performed by the finite element method. “Partial integration point” refers to a part of the points where stress and strain are calculated in the calculation unit of the numerical integration method based on the finite element method. Both points are included.

また、「板厚、弾性係数、塑性係数、応力の成分値、歪の成分値の少なくとも一つ以上の変数」とは、被加工物の一部の板厚、弾性係数(ヤング率、ポアソン比)、塑性係数(降伏応力、塑性係数、加工硬化指数)、解析する3次元座標系(x、y、z)としたときに、x、y、z方向の応力若しくは歪(それぞれ3成分)、xy平面、yz平面、zx平面の剪断応力若しくは剪断歪(それぞれ3成分)、の計18の因子のうち少なくとも一つ以上をいう。   In addition, “at least one variable of plate thickness, elastic modulus, plasticity coefficient, stress component value, and strain component value” means the plate thickness, elastic modulus (Young's modulus, Poisson's ratio) of a part of the workpiece. ), Plastic coefficient (yield stress, plastic coefficient, work hardening index), stress or strain in x, y, z direction (each of three components), when a three-dimensional coordinate system (x, y, z) is analyzed. It means at least one of a total of 18 factors including the xy plane, the yz plane, and the shear stress or shear strain (each of three components) in the zx plane.

図9は、有限要素法で使用される座標系を示す図である。(a)は、3次元全体の座標系を、(b)は、局所座標系を示す。
このとき、図9(a)で示すように、(x、y、z)座標軸での成分値の他に、例えば同図(b)に示すように、要素31の平面をxy平面とみなした局所座標系(x1、x2、x3)における応力、歪、剪断応力、剪断歪の計12成分も含むものとする。また、相当応力、相当塑性歪、弾性エネルギー(弾性仕事)、塑性エネルギー(塑性仕事)等の各応力成分や歪成分を使って得られる変数や、応力増分やひずみ増分等の各成分値の増分形で計算された成分値も含むものとする。
FIG. 9 is a diagram showing a coordinate system used in the finite element method. (A) shows the coordinate system of the whole three-dimensional, (b) shows a local coordinate system.
At this time, as shown in FIG. 9A, in addition to the component values on the (x, y, z) coordinate axes, for example, as shown in FIG. 9B, the plane of the element 31 is regarded as the xy plane. It is assumed to include a total of 12 components of stress, strain, shear stress, and shear strain in the local coordinate system (x1, x2, x3). Also, variables obtained using each stress component and strain component such as equivalent stress, equivalent plastic strain, elastic energy (elastic work), plastic energy (plastic work), and increment of each component value such as stress increment and strain increment It also includes component values calculated in the form.

一般に、鋼板を解析対象とする有限要素法解析では、その鋼板の物性値・物理量は、2mm〜4mm角の有限要素に分割される。しかし、要素の分割単位は、物性値・物理量が一定と近似される長さで分割されるため、必ずしも2mm〜4mm角に限定されるものでは無い。つまり、残留応力が大きい部位は、有限要素をより小さく限定することが必要な場合もある。また、各要素は、3次元座標面で規定され、有限要素の面の角度や曲率は、他の面との比較によって規定される。   In general, in a finite element method analysis using a steel plate as an analysis target, the physical property value / physical quantity of the steel plate is divided into 2 mm to 4 mm square finite elements. However, the division unit of the element is not necessarily limited to 2 mm to 4 mm square because it is divided by a length that approximates that the physical property value / physical quantity is constant. In other words, it may be necessary to limit the finite element to be smaller in a region where the residual stress is large. Each element is defined by a three-dimensional coordinate plane, and the angle and curvature of the plane of the finite element are defined by comparison with other planes.

このように、スプリングバック発生原因部位の特定のために、角度や曲率に基づいて部材を領域分割し、残留応力レベルの高い部分を演算対象とすることで、計算量を大幅に減らすことにより、より迅速なスプリングバック発生原因部位の特定が可能である。   Thus, in order to identify the springback occurrence cause part, the member is divided into regions based on the angle and the curvature, and by calculating the portion having a high residual stress level, the calculation amount is greatly reduced, It is possible to identify the cause of the occurrence of springback more quickly.

演算処理の例として、σx成分値を係数倍する方法を説明する。
選択された領域の積分点における演算前の応力成分を(σx0,σy0,σz0,τxy0,τyz0,τzx0)、歪成分を(εx0,εy0,εz0,γxy0,γyz0,γzx0)とすると、演算処理後の応力成分(σx,σy,σz,τxy,τyz,τzx)、歪成分(εx,εy,εz,γxy,γyz,γzx)は以下のとおりとなる。
σx=K×σx0,σy=σy0,σz=σz0,τxy=τxy0,τyz=τyz0,τzx=τzx0
εx=εx0,εy=εy0,εz=εz0,γxy=γxy0,γyz=γyz0,γzx=γzx0
ここでKは−1000〜1000までの範囲で変化させた演算を行い、板厚は演算処理前の値をそのまま使うこととし、弾性係数、塑性係数は、成形条件入力部18で入力した値を使う場合が一般的である。例示のため、Kは、σxに関してのみ示すが、他の成分についても、同様に変化させても良い。
As an example of arithmetic processing, a method of multiplying the σ x component value by a coefficient will be described.
The stress component before calculation at the integration point of the selected region is (σ x0 , σ y0 , σ z0 , τ xy0 , τ yz0 , τ zx0 ), and the strain component is (ε x0 , ε y0 , ε z0 , γ xy0 , γ yz0 , γ zx0 ), the stress components (σ x , σ y , σ z , τ xy , τ yz , τ zx ), strain components (ε x , ε y , ε z , γ xy , γ yz , γ zx ) is as follows.
σ x = K × σ x0, σ y = σ y0, σ z = σ z0, τ xy = τ xy0, τ yz = τ yz0, τ zx = τ zx0
ε x = ε x0 , ε y = ε y0 , ε z = ε z0 , γ xy = γ xy0 , γ yz = γ yz0 , γ zx = γ zx0
Here, K is calculated in a range of −1000 to 1000, the value before the processing is used as it is for the plate thickness, and the elastic coefficient and the plastic coefficient are values input by the molding condition input unit 18. It is common to use it. For illustration, K is shown only for σ x , but other components may be varied as well.

また、材料特性(板厚、弾性係数、塑性係数)を演算対象とすることも可能である。例えば、物性値・物理量演算処理部14において、プレス成形解析部11から得られた被加工物の一部の領域を選択し、選択された領域の例えば、ヤング率を演算処理において係数倍する。その場合は選択された領域の演算前の板厚t0、ヤング率E0、ポアソン比ν0、塑性係数F0、a0、n0(σ=F0(ε+a0n0の場合)とすると、演算後の板厚t、ヤング率E、ポアソン比ν、塑性係数F、a、n(σ=F(ε+a)nの場合)は以下のとおりとなる。
t=t0,E=K×E0,ν=ν0,F=F0,a=a0,n=n0
ここでKは−1000〜1000の範囲で変化させることができる。例示のため、Kは、Eに関してのみ示すが、他の材料特性についても、同様に変化させても良い。
In addition, material characteristics (plate thickness, elastic modulus, plasticity coefficient) can be calculated. For example, the physical property value / physical quantity calculation processing unit 14 selects a partial region of the workpiece obtained from the press forming analysis unit 11 and multiplies the selected region, for example, Young's modulus by a coefficient in the calculation processing. In that case, the thickness t 0 before calculation of the selected region, Young's modulus E 0 , Poisson's ratio ν 0 , plastic coefficients F 0 , a 0 , n 0 (when σ = F 0 (ε + a 0 ) n0 ) Then, the thickness t, Young's modulus E, Poisson's ratio ν, plastic coefficient F, a, n (when σ = F (ε + a) n ) after calculation are as follows.
t = t 0 , E = K × E 0 , ν = ν 0 , F = F 0 , a = a 0 , n = n 0
Here, K can be changed in the range of −1000 to 1000. For illustration, K is shown only for E, but other material properties may be varied as well.

また、物性値・物理量演算処理部14では、プレス成形解析部11から得られた状態変数を出力したファイルを直接編集により変更することで領域を選択して演算処理することも可能である。例えばファイルの内容をワードパット等のテキスト編集機能を有するソフトで表示し、成形条件入力部18の操作により演算処理を行いたい領域の成分値を直接書き換えたり、カットアンドペーストにより成分値を入れ替えても良い。
これにより、選択した領域の演算した変数(成分)がスプリングバックにどのくらい影響を与えているかを定量的に評価できる。
Further, the physical property value / physical quantity calculation processing unit 14 can also select and perform calculation processing by selecting a region by directly editing the file that outputs the state variable obtained from the press forming analysis unit 11. For example, the contents of the file are displayed by software having a text editing function such as word pad, and the component value of the area to be processed is directly rewritten by operating the molding condition input unit 18, or the component value is replaced by cut and paste. Also good.
Thereby, it is possible to quantitatively evaluate how much the variable (component) calculated in the selected region affects the springback.

例えばスプリングバック解析部12で有限要素節点の変位分布をコンタ図(等高線図)で表示した場合や変形形状の切断面を表示した場合は、選択した領域の演算した変数(成分)ごとの結果の図を、スプリングバック量出力画面19上に並べて比較することや図示しないプリンタに出力して比較することが可能である。また、任意の指定位置の節点の座表値(X、Y、Z)をソフト上またはファイル出力により確認した場合は、選択した領域の演算した変数(成分)ごとの数値を比較したり、表計算ソフト等を用いてグラフ化して比較することが可能である。   For example, when the displacement distribution of a finite element node is displayed in a contour diagram (contour map) or a deformed cut surface is displayed in the springback analysis unit 12, the result of each calculated variable (component) in the selected region is displayed. It is possible to compare the figures side by side on the springback amount output screen 19 or output them to a printer (not shown) for comparison. In addition, when the coordinate value (X, Y, Z) of the node at any specified position is confirmed on the software or by file output, the numerical value for each calculated variable (component) in the selected area is compared, It is possible to make a graph and compare using calculation software or the like.

スプリングバック量出力画面19上に表示されるコンタ図は、物性値・物理量演算処理部14による物性値・物理量演算処理、及び、スプリングバック解析部12によるスプリングバック量計算を繰り返し行いながら、スプリングバックの原因となる部位や変数を厳密に特定していくことが出来る。   The contour diagram displayed on the springback amount output screen 19 shows the springback while repeating the physical property value / physical quantity calculation processing by the physical property value / physical quantity calculation processing unit 14 and the springback amount calculation by the springback analysis unit 12. It is possible to strictly identify the parts and variables that cause the problem.

さらに、物性値・物理量演算処理部14の特定領域を選択により変更したり、選択された領域の成分値をσxだけではなくσyやτxyなどの成分値についても演算処理として係数倍してから、スプリングバック解析部12のスプリングバック解析を実施し、スプリングバック量出力画面18にスプリングバックの原因となる部位や変数を表示させ、スプリングバックによる変位の変化量を比較することで、スプリングバックの原因となる部位や変数を厳密に特定していくことが可能である。   Further, the specific region of the physical property value / physical quantity calculation processing unit 14 is changed by selection, or the component value of the selected region is multiplied not only by σx but also by component coefficients such as σy and τxy as calculation processing, The springback analysis unit 12 performs the springback analysis, displays the portion and variables that cause the springback on the springback amount output screen 18, and compares the amount of change in the displacement due to the springback, thereby causing the cause of the springback. It is possible to specify exactly the parts and variables that become.

また、物性値・物理量演算処理部14において、被加工物の一部の中で複数の領域を選択して、それぞれの選択した領域に対し演算処理を同時に実施することで原因特定のために必要な解析工数を低減することが可能である。   In addition, the physical property value / physical quantity calculation processing unit 14 selects a plurality of areas in a part of the workpiece, and performs the calculation process on each selected area at the same time, which is necessary for specifying the cause. It is possible to reduce the number of analysis steps.

図10は、本発明の一実施形態によるハット型断面形状部品の領域選択の例を示す斜視図である。
例えば、プレス成形解析部11により得られた被加工物が図5に示すようなハット型断面形状部品の場合に、物性値・物理量演算処理部14において図10に示すような成形後の応力が所定値より大きい複数の領域A1、A2を同時に選択し、選択した領域A1、A2に対してσxを係数倍にする演算処理を実施し、スプリングバック解析部12においてその演算処理結果を用いたスプリングバック解析を行う。
FIG. 10 is a perspective view showing an example of region selection of a hat-shaped cross-sectional component according to an embodiment of the present invention.
For example, when the workpiece obtained by the press molding analysis unit 11 is a hat-shaped cross-sectional part as shown in FIG. 5, the physical property value / physical quantity calculation processing unit 14 has a stress after molding as shown in FIG. A plurality of areas A1 and A2 that are larger than a predetermined value are simultaneously selected, a calculation process for multiplying σx by a factor is performed on the selected areas A1 and A2, and a spring using the calculation process result in the springback analysis unit 12 Perform back analysis.

そのスプリングバック解析結果から、スプリングバックによる変位の変化量を調べて、領域A1、A2に加わっている応力成分値σxがどのくらいスプリングバックに影響しているかを、片側ずつ演算処理をしなくても1回の演算処理のみで評価できる。   From the result of the springback analysis, the amount of change in displacement caused by the springback is examined, and it is possible to determine how much the stress component value σx applied to the regions A1 and A2 has an influence on the springback without performing calculation processing for each side. It can be evaluated with only one calculation process.

さらにスプリングバック解析部12のスプリングバック解析を行った後に最もスプリングバック量が変化するときの領域を再分割して、再分割した領域の大きさが所定値以下になるまで、判定して、再分割したそれぞれの領域について、物性値・物理量演算処理部14において演算処理を行い、再度スプリングバック解析部12でスプリングバック解析を行うことで、最もスプリングバック量が変化する領域をスプリングバック発生の原因を特定することができ、効率的な原因特定が可能である。なお、物性値・物理量演算処理部14は、スプリングバック変化量が閾値以下となるように、物性値及び物理量を操作変数として変化させるような収束計算を行うことができ、同時に、特定領域の分割領域も操作変数として変化させるような収束計算を行っても良い。   Further, after the springback analysis of the springback analysis unit 12 is performed, the region where the amount of springback changes the most is subdivided, and determination is performed until the size of the subdivided region becomes a predetermined value or less. For each divided area, the physical property value / physical quantity calculation processing unit 14 performs calculation processing, and the springback analysis unit 12 performs the springback analysis again. Can be identified, and efficient cause identification is possible. The physical property value / physical quantity calculation processing unit 14 can perform convergence calculation so that the physical property value and the physical quantity are changed as manipulated variables so that the springback change amount is equal to or less than the threshold, and at the same time, the division of the specific region is performed. Convergence calculation may be performed so as to change the region as an operation variable.

また、プレス成形解析部11において、例えば有限要素法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部14において領域を選択する際に1個以上の要素を選択して演算対象の領域とすることができる。場の問題を解くための解析手法である有限要素法においては、連続体を有限の要素に分割して扱うのが特徴であり、要素はその辺上に位置する有限個の節点において連結され、個々の要素に定義された形状関数と節点の変位から要素内の変位場が一意的に決定されるものである。例えば、プレス成形解析部11で得られた応力成分σxの分布を表示し、σx値が最も大きい要素を選択領域としたり、σxが最も大きい要素とその要素に接する要素の集団を選択領域とすることができる。   In addition, when the press forming analysis unit 11 performs forming analysis using, for example, a finite element method and numerically obtains the stress and strain value distribution after forming, the physical property value / physical quantity calculation processing unit 14 selects a region. In this case, one or more elements can be selected as a calculation target area. The finite element method, which is an analysis method for solving field problems, is characterized by dividing the continuum into finite elements, and the elements are connected at a finite number of nodes located on their sides. The displacement field in the element is uniquely determined from the shape function defined for each element and the displacement of the node. For example, the distribution of the stress component σx obtained by the press forming analysis unit 11 is displayed, and an element having the largest σx value is selected as a selected area, or a group of elements having the largest σx and elements in contact with the element is selected. be able to.

また、プレス成形解析部11において、例えば初等解析手法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部14において領域を選択する際に1個以上の計算単位区分を選択して演算対象の領域とすることができる。   In addition, when the press forming analysis unit 11 performs forming analysis using, for example, an elementary analysis method and numerically obtains the stress and strain value distribution after forming, the physical property value / physical quantity calculation processing unit 14 selects a region. In this case, one or more calculation unit sections can be selected as a calculation target area.

また、プレス成形解析部11において、例えば有限要素法を用いて成形解析を行い、成形後の応力、歪値の分布を数値的に得た場合、物性値・物理量演算処理部14において領域を選択する際に1個以上の積分点を選択して演算対象の領域とすることができる。例えば、プレス成形解析部11の有限要素解析において被加工物は板厚方向に5個の積分点(板裏面から表面に向けて1、2、3、4、5)を有するシェル要素で解析しておき、物性値・物理量演算処理部14において選択した要素の板裏面に最も近い1と板表面に最も近い5の積分点のみを演算対象の領域とすることができる。これにより、金型の曲率を有する部分に接触して変形した際に受ける曲げ変形の影響を分離した評価を行える。   In addition, when the press forming analysis unit 11 performs forming analysis using, for example, a finite element method and numerically obtains the stress and strain value distribution after forming, the physical property value / physical quantity calculation processing unit 14 selects a region. In this case, one or more integration points can be selected as a calculation target region. For example, in the finite element analysis of the press forming analysis unit 11, the workpiece is analyzed with a shell element having five integration points (1, 2, 3, 4, 5 from the back surface to the front surface) in the thickness direction. In addition, it is possible to set only 1 integration point closest to the plate back surface of the element selected in the physical property value / physical quantity calculation processing unit 14 and 5 integration points closest to the plate surface as a calculation target region. Thereby, the evaluation which isolate | separated the influence of the bending deformation received when contacting and deform | transforming the part which has a curvature of a metal mold | die can be performed.

図11は、本発明の一実施例によるスプリングバック原因部位特定処理のスプリングバック発生原因特定処理のフローチャートを示す。   FIG. 11 is a flowchart of the springback generation cause specifying process of the springback cause site specifying process according to the embodiment of the present invention.

ステップS101では、成形条件入力部18に成形条件が入力される。次に、ステップS102に進む。   In step S <b> 101, molding conditions are input to the molding condition input unit 18. Next, the process proceeds to step S102.

ステップS102では、プレス成形解析部11が、成形条件で規定された被加工物に対して数値解析処理を行い、プレス成形する製品の成形後の応力、歪値の分布を計算する。次に、ステップS103に進む。   In step S102, the press molding analysis unit 11 performs numerical analysis processing on the workpiece defined by the molding conditions, and calculates the distribution of stress and strain values after molding of the product to be press molded. Next, the process proceeds to step S103.

ステップS103では、物性値・物理量演算処理部14で、変換対象領域となる特定部位が決定される。この特定部位は、プレス成形品について曲率及び/又は角度を算定し、前記曲率及び/又は角度に基づいて成形品の領域を分割し、分割したそれぞれの領域について前記成形データに基づいて判定して、後述の演算処理を行う領域を決定する。次に、ステップS104に進む。   In step S <b> 103, the physical property value / physical quantity calculation processing unit 14 determines a specific part to be a conversion target region. The specific part is calculated by calculating the curvature and / or angle of the press-molded product, dividing the region of the molded product based on the curvature and / or angle, and determining each divided region based on the molding data. Then, an area for performing arithmetic processing described later is determined. Next, the process proceeds to step S104.

ステップS104では、S103で演算処理を行うと判定した一部の領域における物性値及び物理量の少なくとも一つを変換する上述の演算処理が行われる。次に、ステップS105に進む。   In step S104, the above-described arithmetic processing for converting at least one of the physical property value and the physical quantity in the partial area determined to be arithmetic processing in S103 is performed. Next, the process proceeds to step S105.

ステップS105では、プレス成形時の被加工体の変換処理された物性値及び物理量に基づいて、スプリングバック量が計算される。次に、ステップS106に進み、同時にステップS107でスプリングバック量の結果がコンタ表示などにより画面表示される。   In step S105, the springback amount is calculated based on the physical property value and the physical quantity subjected to the conversion processing of the workpiece during press molding. Next, the process proceeds to step S106, and at the same time, the result of the springback amount is displayed on the screen by a contour display or the like in step S107.

ステップS106では、スプリングバックの許容値が許容値以内か判定を行う。許容値以内であれば、本処理を終了する。許容値以上であれば、ステップS108に進む。なお、ステップS106及びステップS108は、人間が、特定部位の指定を、スプリングバック変化量を見ながら手動で行っても良い。   In step S106, it is determined whether the allowable value of springback is within the allowable value. If it is within the allowable value, this process is terminated. If it is equal to or greater than the allowable value, the process proceeds to step S108. In step S106 and step S108, a human may manually designate a specific part while watching the amount of change in springback.

ステップS108では、特定部位及び被加工体の物性値・物理量を変更する。ここでは、特定部位のみの変更を行ってもよく、物性値・物理量の変更だけを行っても良く、それら双方を同時に行っても良い。そして、ステップS105に戻る。ステップS105、ステップS106、ステップS108の処理は、スプリングバック変化量が閾値以下になるような収束計算により行うことが可能で、閾値以下になるまで繰り返される。   In step S108, the physical property value / physical quantity of the specific part and the workpiece are changed. Here, only a specific part may be changed, only the physical property value / physical quantity may be changed, or both of them may be simultaneously performed. Then, the process returns to step S105. The processing of step S105, step S106, and step S108 can be performed by a convergence calculation such that the amount of springback change is less than or equal to the threshold value, and is repeated until it becomes less than or equal to the threshold value.

なお、本処理の繰り返し回数を制限し、制限した繰り返し回数を超えた場合は、処理を終了しても良い。その場合、ファイル格納部16に格納されたスプリングバック解析部12の出力データファイル「SB rem.casen.k」を人間が検討し、人間がスプリングバック量出力画面19を見ながら、最適な特定部位を探索して求めることが可能となる。   Note that the number of repetitions of this process is limited, and the process may be terminated when the limited number of repetitions is exceeded. In that case, the human examines the output data file “SB rem.casen.k” of the springback analysis unit 12 stored in the file storage unit 16, and the person looks at the springback amount output screen 19 to determine the optimum specific part. Can be obtained by searching.

図12は、本発明の一実施形態によるスプリングバック発生原因部位特定装置のハードウェア構成図を示す。上述したプレス成形解析部11、スプリングバック解析部12、物性値・物理量演算処理部14における各処理は、プログラム100に規定され、そのプログラム100をコンピュータ90に実行させても良い。コンピュータ90は、必要な処理を実行するCPU91、処理結果を格納するメモリ92(例えば、RAM(Random Access Memory))、ディスプレイ93、例えば、キーボードやマウスのような入力装置94、ハードディスク95、CD/DVDドライブのような外部記憶装置96、NIC(ネットワーク・インタフェース・カード)97、プリンタ98を備える。コンピュータ90は、NIC97に接続されたイーサネット(登録商標)ケーブルにより構成されるネットワーク99を介して、他のコンピュータ90Aと接続可能である。   FIG. 12 shows a hardware configuration diagram of a springback generation cause identifying apparatus according to an embodiment of the present invention. Each processing in the press forming analysis unit 11, the spring back analysis unit 12, and the physical property value / physical quantity calculation processing unit 14 described above may be defined in the program 100, and the program 90 may be executed by the computer 90. The computer 90 includes a CPU 91 for executing necessary processing, a memory 92 for storing processing results (for example, RAM (Random Access Memory)), a display 93, for example, an input device 94 such as a keyboard and a mouse, a hard disk 95, a CD / CD An external storage device 96 such as a DVD drive, a NIC (network interface card) 97, and a printer 98 are provided. The computer 90 can be connected to another computer 90 </ b> A via a network 99 configured by an Ethernet (registered trademark) cable connected to the NIC 97.

プログラム100は、記録媒体に保存され、外部記憶装置96からローディングされ、もしくは、他のコンピュータ90Aからネットワーク99を介してダウンロードされ、CPU91の制御によって、コンピュータ90のハードディスク95に保存される。次に、保存されたプログラム100は、CPU91によって、実行され、実行プロセス又はスレッドとしてメモリ92に格納される。例えば、プレス成形解析部11、スプリングバック解析部12、物性値・物理量演算処理部14における各処理が、それぞれ実行プロセス又はスレッドととなり、各実行プロセス又はスレッド間で上述したデータファイル、又は、データが入出力される。また、それらの各実行プロセス又はスレッドは、他のコンピュータ90Aに分散して存在し、各処理がコンピュータ90と他のコンピュータ90Aによって分散処理されても良い。   The program 100 is stored in a recording medium, loaded from the external storage device 96, or downloaded from another computer 90A via the network 99, and stored in the hard disk 95 of the computer 90 under the control of the CPU 91. Next, the saved program 100 is executed by the CPU 91 and stored in the memory 92 as an execution process or thread. For example, each process in the press molding analysis unit 11, the springback analysis unit 12, and the physical property value / physical quantity calculation processing unit 14 is an execution process or thread, and the data file or data described above between the execution processes or threads. Are input and output. Each of these execution processes or threads may be distributed in the other computer 90A, and each process may be distributed by the computer 90 and the other computer 90A.

また、図3に示される成形条件入力部18、スプリングバック量出力画面19は、それぞれ、入力装置94、ディスプレイ93であっても良い。上述したファイル格納部16は、ハードディスク95であっても良い。プログラム100は、ハードディスク95に保存しても良い。そして、ディスプレイ93に出力された上述のコンタ図は、プリンタ98に出力することが可能である。   Further, the molding condition input unit 18 and the springback amount output screen 19 shown in FIG. 3 may be an input device 94 and a display 93, respectively. The file storage unit 16 described above may be a hard disk 95. The program 100 may be stored in the hard disk 95. The above contour map output to the display 93 can be output to the printer 98.

以下、実施例を用いて本発明を具体的に説明する。   Hereinafter, the present invention will be specifically described with reference to examples.

図13は、本発明の一実施例によるプレス成形品である被加工物の例を示す斜視図と断面図である。図13に示すように、解析対象となる被加工物は、ハット型断面形状部品である。
まずプレス成形解析処理として、有限要素法に基づく市販の板成形シミュレーション解析ソフトウェアLS−DYNAを使用して解析処理を行った。金属板の性状として、板厚1.4mm、引張強さ590MPa級の高強度鋼板のデータを用いた。また、金型(ダイ、パンチ、ホルダー)の形状をシェル要素によりモデリングし、剛体と仮定して解析した。金型クリアランスは0mmとした。摩擦係数は0.1を入力し、成形荷重として1400kNを設定した。
FIG. 13 is a perspective view and a cross-sectional view showing an example of a workpiece that is a press-formed product according to an embodiment of the present invention. As shown in FIG. 13, the workpiece to be analyzed is a hat-shaped cross-sectional component.
First, as press forming analysis processing, analysis processing was performed using commercially available plate forming simulation analysis software LS-DYNA based on the finite element method. As the properties of the metal plate, data on a high-strength steel plate having a thickness of 1.4 mm and a tensile strength of 590 MPa was used. The shape of the die (die, punch, holder) was modeled with shell elements and analyzed assuming a rigid body. The mold clearance was 0 mm. As the friction coefficient, 0.1 was input and 1400 kN was set as the forming load.

物性値・物理量演算処理を実行するプログラムは、プレス成形解析から得られた応力や歪が出力されたファイルを、入力情報として取り込み、演算処理を行うために被加工物の領域を分割する。   A program that executes physical property value / physical quantity calculation processing takes in a file in which stress or strain obtained from press forming analysis is output as input information, and divides the region of the workpiece to perform calculation processing.

図14は、本発明の一実施例による図13に示した被加工物をX方向に10等分、Y、Z方向に一定の曲率で分割した場合の分割領域を示す斜視図である。ここでは、各領域に属する要素の全積分点に演算手段により以下の演算を施した。
σx=0,σy=0,σz=0,τxy=0,τyz=0,τzx=0
εx=εx0,εy=εy0,εz=εz0,γxy=γxy0,γyz=γyz0,γzx=γzx0
FIG. 14 is a perspective view showing a divided area when the workpiece shown in FIG. 13 according to an embodiment of the present invention is divided into 10 equal parts in the X direction and constant curvature in the Y and Z directions. Here, the following calculation was performed by calculation means on all integration points of elements belonging to each region.
σ x = 0, σ y = 0, σ z = 0, τ xy = 0, τ yz = 0, τ zx = 0
ε x = ε x0 , ε y = ε y0 , ε z = ε z0 , γ xy = γ xy0 , γ yz = γ yz0 , γ zx = γ zx0

ここで、選択された領域の積分点における演算前の応力成分を(σx0,σy0,σz0,τxy0,τyz0,τzx0)、歪成分を(εx0,εy0,εz0,γxy0,γyz0,γzx0)と表し、演算処理後における応力成分を(σx,σy,σz,τxy,τyz,τzx)、歪成分を(εx,εy,εz,γxy,γyz,γzx)とした。板厚は、演算処理前の値をそのまま使うこととし、弾性係数、塑性係数はプレス解析処理におけるソフトウェアLS−DYNAに入力した値をそのまま使用した。演算された応力や歪値はファイル出力手段により演算結果ファイルとして出力した。 Here, the stress component before calculation at the integration point of the selected region is (σ x0 , σ y0 , σ z0 , τ xy0 , τ yz0 , τ zx0 ), and the strain component is (ε x0 , ε y0 , ε z0 , γ xy0 , γ yz0 , γ zx0 ), the stress components after the calculation processing are (σ x , σ y , σ z , τ xy , τ yz , τ zx ), and the strain components are (ε x , ε y , ε z , γ xy , γ yz , γ zx ). For the plate thickness, the value before the calculation processing is used as it is, and the value input to the software LS-DYNA in the press analysis processing is used as it is for the elastic coefficient and plasticity coefficient. The calculated stress and strain values were output as a calculation result file by the file output means.

スプリングバック解析処理として、上述のソフトウェアLS−DYNAを使用した。上記の物性値・物理量演算処理実行プログラムの出力結果を、ソフトウェアLS−DYNAに入力してスプリングバック解析を実施した。スプリングバック解析は静的陰解法による弾性解析を用いた。各領域での演算処理・スプリングバック解析を領域分割数繰り返した。   The above software LS-DYNA was used as the springback analysis process. The output result of the physical property value / physical quantity calculation processing execution program was input to the software LS-DYNA to perform a springback analysis. The springback analysis used the elastic analysis by the static implicit method. The calculation process and springback analysis in each area were repeated for the number of area divisions.

図15は、本発明の一実施例によるスプリングバック前後の断面A、及び、注目点BのZ方向変位であるスプリングバック量を示す図である。図16は、本発明の一実施例による演算処理を行った各領域でのスプリングバック量を示すコンタ図である。このコンタ図は、ソフトウェアLS−DYNAの出力データにより、ポスト処理ソフトLS−PREPOSTにより得たものである。図16より、スプリングバック量の大きい領域Cがスプリングバック発生原因部位であると特定することができた。   FIG. 15 is a diagram showing a cross section A before and after the spring back according to an embodiment of the present invention and a spring back amount that is a Z direction displacement of the point of interest B. FIG. 16 is a contour diagram showing the amount of springback in each region where the arithmetic processing according to one embodiment of the present invention has been performed. This contour diagram is obtained by the post-processing software LS-PREPOST from the output data of the software LS-DYNA. From FIG. 16, it was possible to specify that the region C where the amount of springback is large is the site causing the springback.

対象部品は図13に示す形状及び寸法であるハット型断面形状部品とする。
まずプレス成形解析には、上述のソフトウェアLS−DYNAを使用した。金属板の性状として、板厚1.4mm、引張強さ590MPa級の高強度鋼板のデータを用いた。また、金型(ダイ、パンチ、ホルダー)の形状をシェル要素によりモデリングし、剛体と仮定して解析した。金型クリアランスは0mmとした。摩擦係数は0.1を入力し、成形荷重として1400kNを設定した。
The target part is a hat-shaped cross-sectional part having the shape and dimensions shown in FIG.
First, the above-mentioned software LS-DYNA was used for press molding analysis. As the properties of the metal plate, data on a high-strength steel plate having a thickness of 1.4 mm and a tensile strength of 590 MPa was used. The shape of the die (die, punch, holder) was modeled with shell elements and analyzed assuming a rigid body. The mold clearance was 0 mm. As the friction coefficient, 0.1 was input and 1400 kN was set as the forming load.

物性値・物理量演算処理を実行するプログラムは、プレス成形解析から得られた応力や歪が出力されたファイルを、入力情報として取り込み、演算処理を行うために被加工物の領域を分割する。図14に示す各領域として分割し、各領域に属する要素の全積分点に演算手段により以下の演算を施した。
σx=0,σy=0,σz=0,τxy=0,τyz=0,τzx=0
εx=εx0,εy=εy0,εz=εz0,γxy=γxy0,γyz=γyz0,γzx=γzx0
A program that executes physical property value / physical quantity calculation processing takes in a file in which stress or strain obtained from press forming analysis is output as input information, and divides the region of the workpiece to perform calculation processing. Each region shown in FIG. 14 was divided, and the following calculation was performed by calculation means on all integration points of elements belonging to each region.
σ x = 0, σ y = 0, σ z = 0, τ xy = 0, τ yz = 0, τ zx = 0
ε x = ε x0 , ε y = ε y0 , ε z = ε z0 , γ xy = γ xy0 , γ yz = γ yz0 , γ zx = γ zx0

ここで、選択された領域の積分点における演算前の応力成分を(σx0,σy0,σz0,τxy0,τyz0,τzx0),歪成分を(εx0,εy0,εz0,γxy0,γyz0,γzx0)と表し,演算処理後における応力成分を(σx,σy,σz,τxy,τyz,τzx),歪成分を(εx,εy,εz,γxy,γyz,γzx)とした。板厚は、演算処理前の値をそのまま使うこととし、弾性係数、塑性係数は、プレス成形解析としてソフトウェアLS−DYNAで入力した値をそのまま使用した。演算された応力や歪値はファイル出力手段により演算結果ファイルとして出力した。 Here, the stress component before calculation at the integration point in the selected region is (σ x0 , σ y0 , σ z0 , τ xy0 , τ yz0 , τ zx0 ), and the strain component is (ε x0 , ε y0 , ε z0 , γ xy0 , γ yz0 , γ zx0 ), the stress components after the processing are (σ x , σ y , σ z , τ xy , τ yz , τ zx ), and the strain components are (ε x , ε y , ε z , γ xy , γ yz , γ zx ). For the plate thickness, the values before the arithmetic processing were used as they were, and the values input by the software LS-DYNA as the press molding analysis were used as they were for the elastic coefficient and plastic coefficient. The calculated stress and strain values were output as a calculation result file by the file output means.

スプリングバック解析処理では、物性値・物理量演算処理実行プログラムによる演算結果ファイルを、再度ソフトウェアLS−DYNAに入力してスプリングバック解析を実施した。スプリングバック解析は静的陰解法による弾性解析を用いた。
ここでは、各領域毎に、物性値・物理量演算処理プログラムによる演算処理、及び、ソフトウェアLS−DYNAによるスプリングバック解析を繰り返し行った。
In the springback analysis process, the calculation result file obtained by the physical property value / physical quantity calculation process execution program is input again to the software LS-DYNA to perform the springback analysis. The springback analysis used the elastic analysis by the static implicit method.
Here, for each region, calculation processing by a physical property value / physical quantity calculation processing program and springback analysis by software LS-DYNA were repeatedly performed.

図17は、本発明の一実施例による各領域にスプリングバック影響度をコンタ表示した斜視図である。
図15に示すスプリングバック前後の断面Aに属する注目点BのZ方向変位をスプリングバック量とした。ここで、正規化された値であるスプリングバック影響度として、スプリングバック影響度=スプリングバック量÷各分割領域面積×1000+10と定義し、演算処理を行った各領域でのスプリングバック影響度をポスト処理ソフトにて図17に示すコンタ図として表示した。図17より、スプリングバック量の大きい領域Dがスプリングバック発生原因部位であると特定することができた。
FIG. 17 is a perspective view showing a springback influence level in each area according to an embodiment of the present invention.
The Z-direction displacement of the point of interest B belonging to the cross-section A before and after the springback shown in FIG. Here, the springback influence degree, which is a normalized value, is defined as springback influence degree = springback amount / each divided area area × 1000 + 10, and the springback influence degree in each area where the calculation processing is performed is posted. It was displayed as a contour diagram shown in FIG. 17 by the processing software. From FIG. 17, it was possible to specify that the region D having a large amount of springback was the cause of the occurrence of springback.

なお、上記スプリングバック影響度は、分割領域による正規化した例であるが、プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値とすることもできる。代表寸法を用いて、スプリングバック影響度を求めることで、被加工体の寸法がスプリングバックに与える影響を定量的に判断でき、特定部位の判別と、その対策がより容易になる。例えば、長手方向の代表寸法で規格化したスプリングバック影響度が、横手方向のそれよりも大きな値を示した場合、ビードの形状を長手方向に伸ばした方がよりスプリングバック量を減らすことができる等の情報を容易に得ることが可能となる。   In addition, although the said springback influence degree is the example normalized by the division | segmentation area | region, it divided by any one of the representative length, representative width, representative height, representative board thickness, and tensile strength of a press-formed product. It can also be a value. By determining the influence of the springback using the representative dimensions, the influence of the dimensions of the workpiece on the springback can be quantitatively determined, and the identification of specific parts and the countermeasures can be facilitated. For example, when the influence of springback normalized by the representative dimension in the longitudinal direction shows a larger value than that in the lateral direction, the amount of springback can be reduced more by extending the bead shape in the longitudinal direction. Such information can be easily obtained.

また、スプリングバック影響度は、パンチ速度又はしわ押さえ力で割った値をコンタ表示することもできる。このような変数によりコンタ表示をさせることにより、スプリングバック量を減らすために、被加工体の形状に対する対策のみならず、プレス成形条件に対しての定量的な判断材料が提供され、特定部位へのプレス成形条件による対策をとることが可能である。   In addition, the springback influence degree can be displayed in a contour by dividing the value by the punching speed or the wrinkle holding force. By displaying contours with such variables, in order to reduce the amount of springback, not only measures against the shape of the workpiece but also quantitative judgment materials for press molding conditions are provided, and to specific parts It is possible to take measures according to the press molding conditions.

上述においては、スプリングバック量、及びスプリングバック影響度の表示方法としてコンタ表示を説明したが、スプリングバック量等をベクトル表示、棒等の3次元物体で表示したり、良く知られた表示方法によって表示することが可能である。
コンタ表示するポスト処理プログラムとして、LS−PREPOST,HyperView,Animator,Meta−Post、等を使用することができる。
In the above description, contour display has been described as a method of displaying the amount of springback and the degree of influence of springback. However, the amount of springback or the like is displayed as a vector display or a three-dimensional object such as a bar, It is possible to display.
LS-PREPOST, HyperView, Animator, Meta-Post, etc. can be used as a post-processing program for contour display.

以上に例を挙げて示したように、被加工物のスプリングバック影響度を、定量的に分析し、さらにその情報を表示して視覚的に捉えることにより、スプリングバック発生原因部位を容易に特定することが可能となる。
このように、本発明によるスプリングバック発生原因部位特定装置は、従来実際の装置を用いてスプリングバックの検討を、数値解析によりスプリングバック発生原因部位を特定できるので、成形部材の設計段階でのテスト工数及び費用を減少させるものである。また、このようなスプリングバック発生原因部位特定装置は、被加工物全般に適用されることが期待されることから、産業界において多大な利益をもたらす。
As shown in the examples above, the cause of the springback can be easily identified by quantitatively analyzing the degree of influence of the workpiece on the springback and displaying the information visually. It becomes possible to do.
As described above, the device for identifying the cause of the occurrence of springback according to the present invention can identify the cause of the occurrence of springback by numerical analysis by using the conventional device to examine the springback, so that the test at the design stage of the molded member can be performed. It reduces man-hours and costs. Further, such a device for identifying the cause of occurrence of springback is expected to be applied to all workpieces, and thus brings great benefits to the industry.

以上説明した実施形態は典型例として挙げたに過ぎず、その各実施形態の構成要素を組合せること、その変形及びバリエーションは当業者にとって明らかであり、当業者であれば本発明の原理及び請求の範囲に記載した発明の範囲を逸脱することなく上述の実施形態の種趣の変形を行えることは明らかである。   The embodiments described above are merely given as typical examples, and it is obvious to those skilled in the art to combine the components of each embodiment, and variations and variations thereof. Those skilled in the art will understand the principles and claims of the present invention. It is obvious that various modifications of the above-described embodiment can be made without departing from the scope of the invention described in the above.

従来のスプリングバックへの対応を示す成形部材の断面模式図である。It is a cross-sectional schematic diagram of a molded member showing correspondence to a conventional spring back. スプリングバックによる変形の発生原因となる部位を探る従来の方法を例示する斜視図である。It is a perspective view which illustrates the conventional method which searches the site | part which becomes the generation | occurrence | production of a deformation | transformation by a springback. 本発明の一実施形態によるスプリングバック発生原因部位特定装置の機能構成図を示す。The functional block diagram of the spring back generation | occurrence | production cause site | part identification apparatus by one Embodiment of this invention is shown. 本発明例によるプレス成形解析結果の応力分布を示すコンタ図である。It is a contour figure which shows the stress distribution of the press molding analysis result by the example of this invention. 本発明例による各要素の曲率を示すコンタ図である。It is a contour figure which shows the curvature of each element by the example of this invention. 本発明例による曲率の大きい要素を非表示としたコンタ図である。FIG. 6 is a contour diagram in which elements having a large curvature according to an example of the present invention are not displayed. 本発明例による曲率に基づいて領域分割したことを示す図である。It is a figure which shows having divided into areas based on the curvature by the example of this invention. 本発明例による一部の領域の応力をゼロとしたコンタ図である。It is the contour figure which made the stress of the one part area | region by the example of this invention zero. 本発明の一実施形態による有限要素法で使用される座標系を示す図である。It is a figure which shows the coordinate system used by the finite element method by one Embodiment of this invention. 本発明の一実施形態による被加工物の領域選択の例を示す斜視図である。It is a perspective view which shows the example of the area | region selection of the workpiece by one Embodiment of this invention. 本発明の一実施形態によるスプリングバック発生原因部位特定処理のフローチャートである。It is a flowchart of the spring back generation | occurrence | production cause site | part identification process by one Embodiment of this invention. 本発明の一実施形態によるスプリングバック原因部位特定処理のハードウェア構成図である。It is a hardware block diagram of the springback cause site | part identification process by one Embodiment of this invention. 本発明の一実施例による被加工物の例を示す斜視図と断面図である。It is the perspective view and sectional drawing which show the example of the workpiece by one Example of this invention. 本発明の一実施例による被加工物の領域分割の例を示す斜視図である。It is a perspective view which shows the example of the area | region division of the workpiece by one Example of this invention. 本発明の一実施例による被加工物の断面位置と注目点を示す斜視図及び断面図である。It is the perspective view and sectional drawing which show the cross-sectional position and attention point of the workpiece by one Example of this invention. 本発明の一実施例による被加工物の各領域にスプリングバック量をコンタ表示した例を示す斜視図である。It is a perspective view which shows the example which displayed the amount of springbacks in each area | region of the to-be-processed object by one Example of this invention. 本発明の一実施例による被加工物の各領域にスプリングバック影響度をコンタ表示した例を示す斜視図である。It is a perspective view which shows the example which carried out the contour display of the springback influence degree in each area | region of the workpiece by one Example of this invention.

符号の説明Explanation of symbols

10A,10B スプリングバック発生原因部位特定装置
11 プレス成形解析部
12 スプリングバック解析部
14 物性値・物理量演算処理部
16 ファイル格納部
18 成形条件入力部
19 スプリングバック量出力画面
10A, 10B Springback occurrence cause identifying device 11 Press molding analysis unit 12 Springback analysis unit 14 Physical property value / physical quantity calculation processing unit 16 File storage unit 18 Molding condition input unit 19 Springback amount output screen

Claims (22)

コンピュータが、プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析ステップと、
コンピュータが、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理ステップと、
コンピュータが、前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップと、
コンピュータが、前記各領域毎に、前記算出したスプリングバック量を前記各領域の面積で割った値を表示する表示ステップと、
を有することを特徴とするスプリングバック影響度表示方法。
A press molding analysis step in which a computer numerically analyzes press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing step in which the computer performs arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A computer calculates a springback amount based on a result of the arithmetic processing, and calculates a springback amount;
Computer, a display step in which the in each area, to display the value obtained by dividing the amount of spring back which is the calculated area of the respective regions,
A method for displaying the influence of springback, characterized by comprising:
コンピュータが、プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析ステップと、
コンピュータが、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理ステップと、
コンピュータが、前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップと、
コンピュータが、前記各領域毎に、前記算出したスプリングバック量を前記プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値を表示するステップと、
を有することを特徴とするスプリングバック影響度表示方法。
A press molding analysis step in which a computer numerically analyzes press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing step in which the computer performs arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A computer calculates a springback amount based on a result of the arithmetic processing, and calculates a springback amount;
The computer displays a value obtained by dividing the calculated springback amount by one of the representative length, representative width, representative height, representative plate thickness, and tensile strength of the press-formed product for each area. And steps to
A method for displaying the influence of springback, characterized by comprising:
コンピュータが、プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析ステップと、
コンピュータが、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理ステップと、
コンピュータが、前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出ステップと、
コンピュータが、前記各領域毎に、前記算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を表示するステップと、
を有することを特徴とするスプリングバック影響度表示方法。
A press molding analysis step in which a computer numerically analyzes press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing step in which the computer performs arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A computer calculates a springback amount based on a result of the arithmetic processing, and calculates a springback amount;
And step computer, wherein for each of the regions, to display the value obtained by dividing the amount of spring back which is the calculated punch speed or blank holding force,
A method for displaying the influence of springback, characterized by comprising:
前記演算処理ステップは、コンピュータが、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行うステップである請求項1〜3の何れか1項に記載のスプリングバック影響度表示方法。 In the calculation processing step, the computer performs calculation processing on at least one of the physical property value and physical quantity data of a part of integration points of each region obtained by dividing the press molded product among the molding data of the press molded product. The springback influence display method according to any one of claims 1 to 3, wherein the step of performing 前記表示ステップは、コンピュータが、前記各領域毎に、前記算出したスプリングバック量をコンタ表示するステップである請求項1〜4の何れか1項に記載のスプリングバック影響度表示方法。 Wherein the display step, the computer, the each region, the springback effect of the display method according to any one of claims 1 to 4, a step for contour display the amount of spring back which is the calculated. 前記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値である請求項1〜5のいずれか一項に記載のスプリングバック影響度表示方法。   6. The springback influence degree display method according to claim 1, wherein the physical property value and the physical quantity are a plate thickness, an elastic coefficient, a plastic coefficient, a stress component value, and a strain component value. プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析部と、
前記成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行ない、前記スプリングバック解析部に、前記演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、
前記各領域毎に前記算出したスプリングバック量を前記各領域の面積で割った値を表示する表示部と、
を有することを特徴とするスプリングバック影響度表示装置。
A press molding analysis unit that numerically analyzes the molding conditions of press molding and obtains molding data of the press molded product,
A numerical analysis of the molding data, and a springback analysis unit for calculating a springback amount;
Among the molding data of the press-molded product, perform arithmetic processing on at least one of the physical property value and physical quantity data of each region into which the press-molded product is divided, and the spring back analysis unit has the result of the arithmetic processing. An arithmetic processing unit for calculating a springback amount based on
A display unit for displaying a value obtained by dividing the calculated springback amount by the area of each region for each region;
A spring back influence display device characterized by comprising:
プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析部と、
前記成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行ない、前記スプリングバック解析部に、前記演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、
前記各領域毎に、前記算出したスプリングバック量を前記プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値を、前記表示部に表示する表示部と、
を有することを特徴とするスプリングバック影響度表示装置。
A press molding analysis unit that numerically analyzes the molding conditions of press molding and obtains molding data of the press molded product,
A numerical analysis of the molding data, and a springback analysis unit for calculating a springback amount;
Of the molding data of the press-molded product, perform arithmetic processing on at least one of the physical property value and physical quantity data of each region into which the press-molded product is divided, and the springback analysis unit has the result of the arithmetic processing. An arithmetic processing unit for calculating a springback amount based on
A value obtained by dividing the calculated springback amount by any one of the representative length, the representative width, the representative height, the representative plate thickness, and the tensile strength of the press-molded product for each region. A display section to be displayed on
A spring back influence display device characterized by comprising:
プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析部と、
前記成形データを数値解析して、スプリングバック量を算出するスプリングバック解析部と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行ない、前記スプリングバック解析部に、前記演算処理の結果に基づくスプリングバック量を算出させる演算処理部と、
前記各領域毎に、前記算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を、前記表示部にコンタ表示する表示部と、
を有することを特徴とするスプリングバック影響度表示装置。
A press molding analysis unit that numerically analyzes the molding conditions of press molding and obtains molding data of the press molded product,
A numerical analysis of the molding data, and a springback analysis unit for calculating a springback amount;
Of the molding data of the press-molded product, perform arithmetic processing on at least one of the physical property value and physical quantity data of each region into which the press-molded product is divided, and the springback analysis unit has the result of the arithmetic processing. An arithmetic processing unit for calculating a springback amount based on
A display unit that contour-displays a value obtained by dividing the calculated springback amount by a punching speed or a wrinkle holding force for each region, on the display unit;
A spring back influence display device characterized by comprising:
前記演算処理部は、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行なう請求項7〜9の何れか1項に記載のスプリングバック影響度表示装置。   The arithmetic processing unit performs arithmetic processing on at least one of physical property values and physical quantity data of a part of integration points in each region obtained by dividing the press molded product among the molding data of the press molded product. Item 10. The springback influence display device according to any one of Items 7 to 9. 前記表示部は、前記各領域毎に前記算出したスプリングバック量をコンタ表示する請求項7〜10の何れか1項に記載のスプリングバック影響度表示装置。   The springback influence display device according to any one of claims 7 to 10, wherein the display unit displays the calculated amount of springback for each of the regions. 前記物性値及び物理量は、板厚、弾性係数、塑性係数、応力の成分値、歪の成分値である請求項7〜11のいずれか一項に記載のスプリングバック影響度表示装置。   The spring back influence display device according to any one of claims 7 to 11, wherein the physical property value and the physical quantity are a plate thickness, an elastic coefficient, a plastic coefficient, a stress component value, and a strain component value. プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析手順と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出手順と、
前記各領域毎に、前記算出したスプリングバック量を前記各領域の面積で割った値を、表示部に表示する表示手順と、
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis procedure for numerically analyzing press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A springback amount calculation procedure for calculating a springback amount based on a result of the arithmetic processing;
A display procedure for displaying a value obtained by dividing the calculated springback amount by the area of each region on the display unit for each region,
A program for displaying the influence of springback, which causes a computer to execute the program.
プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析手順と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出手順と、
前記各領域毎に、前記算出したスプリングバック量を前記プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値を表示する表示手順と、
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis procedure for numerically analyzing press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A springback amount calculation procedure for calculating a springback amount based on a result of the arithmetic processing;
Display procedure for displaying the value obtained by dividing the calculated springback amount by one of the representative length, representative width, representative height, representative plate thickness, and tensile strength of the press-formed product for each region. When,
A program for displaying the influence of springback, which causes a computer to execute the program.
プレス成形の成形条件を数値解析して、プレス成形品の成形データを得るプレス成形解析手順と、
前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記演算処理の結果に基づいて、スプリングバック量を算出するスプリングバック量算出手順と、
前記各領域毎に、前記算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を表示する表示手順と、
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis procedure for numerically analyzing press molding molding conditions to obtain molding data of a press molded product,
An arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product among the molding data of the press molded product,
A springback amount calculation procedure for calculating a springback amount based on a result of the arithmetic processing;
A display procedure for displaying the value obtained by dividing the calculated springback amount by the punching speed or the wrinkle holding force for each region,
A program for displaying the influence of springback, which causes a computer to execute the program.
前記演算処理手順は、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う請求項13〜15の何れか1項に記載のスプリングバック影響度表示プログラム。   The calculation processing procedure performs calculation processing on at least one of physical property values and physical quantity data of a part of integration points in each region obtained by dividing the press-formed product out of the molding data of the press-formed product. Item 16. The springback influence degree display program according to any one of Items 13 to 15. 前記表示手順は、前記各領域毎に、前記算出したスプリングバック量をコンタ表示する請求項13〜16の何れか1項に記載のスプリングバック影響度表示プログラム。   The springback influence display program according to any one of claims 13 to 16, wherein the display procedure contour-displays the calculated springback amount for each of the regions. プレス成形の成形条件を数値解析してプレス成形品の成形データを算出する手順をコンピュータに実行させるプレス成形解析プログラム、該成形データを数値解析してスプリングバック量を算出する手順をコンピュータに実行させるスプリングバック解析プログラム、及び該スプリングバック量をコンタ表示する手順をコンピュータに実行させるポスト処理プログラム、と協働するスプリングバック影響度表示プログラムであって、
前記プレス成形解析プログラムを実行するコンピュータによって出力される前記成形データを取得する手順と、
前記成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記スプリングバック解析プログラムを実行するコンピュータに、前記演算処理の結果を出力する手順と、
前記スプリングバック解析プログラムが各領域毎に前記算出したスプリングバック量を前記各領域の面積で割った値を、前記ポスト処理プログラムを実行するコンピュータに、コンタ表示させるコンタ表示手順
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis program for causing a computer to execute a procedure for calculating the molding data of a press molded product by numerically analyzing the molding conditions of the press molding, and causing the computer to execute a procedure for calculating the springback amount by numerically analyzing the molding data . springback analysis program, and post-processing program for executing a procedure for contour displaying the amount of spring back to the computer, a spring-back effect of the display program for working capital together,
A procedure for acquiring the molding data output by a computer that executes the press molding analysis program;
Among the molding data, an arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product,
A procedure for outputting the result of the arithmetic processing to a computer that executes the springback analysis program;
The springback analysis program by dividing the amount of spring back which is the calculated for each region by the area of each region value, the computer running the post-processing program, and Contours display procedure for contour display,
A program for displaying the influence of springback, which causes a computer to execute the program.
プレス成形の成形条件を数値解析してプレス成形品の成形データを算出する手順をコンピュータに実行させるプレス成形解析プログラム、該成形データを数値解析してスプリングバック量を算出する手順をコンピュータに実行させるスプリングバック解析プログラム、及び該スプリングバック量をコンタ表示する手順をコンピュータに実行させるポスト処理プログラム、と協働するスプリングバック影響度表示プログラムであって、
前記プレス成形解析プログラムを実行するコンピュータによって出力される前記成形データを取得する手順と、
前記成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記スプリングバック解析プログラムを実行するコンピュータに、前記演算処理の結果を出力する手順と、
前記各領域毎に、前記算出したスプリングバック量を前記プレス成形品の代表長さ、代表幅、代表高さ、代表板厚、引張強さのいずれか1つで割った値、前記ポスト処理プログラムを実行するコンピュータに、コンタ表示させるコンタ表示手順と、
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis program for causing a computer to execute a procedure for calculating the molding data of a press molded product by numerically analyzing the molding conditions of the press molding, and causing the computer to execute a procedure for calculating the springback amount by numerically analyzing the molding data . springback analysis program, and post-processing program for executing a procedure for contour displaying the amount of spring back to the computer, a spring-back effect of the display program for working capital together,
A procedure for acquiring the molding data output by a computer that executes the press molding analysis program;
Among the molding data, an arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product,
A procedure for outputting the result of the arithmetic processing to a computer that executes the springback analysis program;
Wherein for each of the regions, the calculated amount of springback the press molded product of the representative lengths, representative width, representative height, representative thickness, the value obtained by dividing the one of tensile strength, the post-processing A contour display procedure for displaying a contour on a computer executing the program;
A program for displaying the influence of springback, which causes a computer to execute the program.
プレス成形の成形条件を数値解析してプレス成形品の成形データを算出する手順をコンピュータに実行させるプレス成形解析プログラム、該成形データを数値解析してスプリングバック量を算出する手順をコンピュータに実行させるスプリングバック解析プログラム、及び該スプリングバック量をコンタ表示する手順をコンピュータに実行させるポスト処理プログラム、と協働するスプリングバック影響度表示プログラムであって、
前記プレス成形解析プログラムを実行するコンピュータによって出力される前記成形データを取得する手順と、
前記成形データのうち、前記プレス成形品を分割した各領域の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う演算処理手順と、
前記スプリングバック解析プログラムを実行するコンピュータに、前記演算処理の結果を出力する手順と、
前記各領域毎に、前記算出したスプリングバック量をパンチ速度又はしわ押さえ力で割った値を、前記ポスト処理プログラムを実行するコンピュータに、コンタ表示させるコンタ表示手順と、
をコンピュータに実行させることを特徴とするスプリングバック影響度表示プログラム。
A press molding analysis program for causing a computer to execute a procedure for calculating the molding data of a press molded product by numerically analyzing the molding conditions of the press molding, and causing the computer to execute a procedure for calculating the springback amount by numerically analyzing the molding data . springback analysis program, and post-processing program for executing a procedure for contour displaying the amount of spring back to the computer, a spring-back effect of the display program for working capital together,
A procedure for acquiring the molding data output by a computer that executes the press molding analysis program;
Among the molding data, an arithmetic processing procedure for performing arithmetic processing on at least one of the physical property value and physical quantity data of each region obtained by dividing the press molded product,
A procedure for outputting the result of the arithmetic processing to a computer that executes the springback analysis program;
A contour display procedure for contour-displaying a value obtained by dividing the calculated springback amount by a punch speed or a wrinkle holding force for each area on a computer that executes the post-processing program,
A program for displaying the influence of springback, which causes a computer to execute the program.
前記演算処理手順は、前記プレス成形品の成形データのうち、前記プレス成形品を分割した各領域の一部の積分点の物性値及び物理量のデータの少なくとも一つに対して演算処理を行う手順である請求項18〜20の何れか1項に記載のスプリングバック影響度表示プログラム。 The calculation processing procedure is a procedure for performing calculation processing on at least one of physical property values and physical quantity data of a part of integration points in each region obtained by dividing the press molded product among the molding data of the press molded product. springback influence display program according to any one of claims 18 to 20 is. 請求項13〜21のいずれか一項に記載のコンピュータプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。   A computer-readable recording medium on which the computer program according to any one of claims 13 to 21 is recorded.
JP2006268855A 2006-08-31 2006-09-29 Springback influence display method, apparatus thereof, and program thereof Active JP4308239B2 (en)

Priority Applications (13)

Application Number Priority Date Filing Date Title
JP2006268855A JP4308239B2 (en) 2006-09-29 2006-09-29 Springback influence display method, apparatus thereof, and program thereof
KR1020097004343A KR101088115B1 (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs
CN200780032577.9A CN101511502B (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method thereof
CA2662109A CA2662109C (en) 2006-08-31 2007-08-30 Method of identification of cause of occurrence of springback, method of display of degree of effect of springback, method of identification of location of cause of occurence of springback, method of identification of position of measure against springback, apparatuses of these, and programs of these
CN201210056854.2A CN102708222B (en) 2006-08-31 2007-08-30 Method of identification of cause of occurrence of springback, method of identification of location of cause of occurence of springback, method of identification of position of measure against springback, and apparatuses of these
CN201410204754.9A CN103995927B (en) 2006-08-31 2007-08-30 Resilience countermeasure position specifying method and resilience countermeasure position specific device
MX2009002150A MX2009002150A (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs.
US12/439,126 US9767234B2 (en) 2006-08-31 2007-08-30 Method of identification of cause and/or location of cause of occurrence of springback
EP07806800.4A EP2058062A4 (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs
TW096132255A TW200824813A (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs
PCT/JP2007/067358 WO2008026777A1 (en) 2006-08-31 2007-08-30 Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs
BRPI0716142A BRPI0716142B1 (en) 2006-08-31 2007-08-30 method of identifying cause of occurrence of elastic return and methods of identifying location of cause of occurrence of elastic return
US15/689,209 US20180075177A1 (en) 2006-08-31 2017-08-29 Method of identification of cause and/or location of cause of occurrence of springback

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006268855A JP4308239B2 (en) 2006-09-29 2006-09-29 Springback influence display method, apparatus thereof, and program thereof

Publications (3)

Publication Number Publication Date
JP2008090481A JP2008090481A (en) 2008-04-17
JP2008090481A5 JP2008090481A5 (en) 2009-01-08
JP4308239B2 true JP4308239B2 (en) 2009-08-05

Family

ID=39374581

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006268855A Active JP4308239B2 (en) 2006-08-31 2006-09-29 Springback influence display method, apparatus thereof, and program thereof

Country Status (1)

Country Link
JP (1) JP4308239B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5821403B2 (en) * 2011-08-22 2015-11-24 Jfeスチール株式会社 Method and apparatus for confirming springback countermeasure effect of press-formed product
JP5797594B2 (en) * 2012-03-22 2015-10-21 株式会社豊田中央研究所 Press forming analysis system and program thereof

Also Published As

Publication number Publication date
JP2008090481A (en) 2008-04-17

Similar Documents

Publication Publication Date Title
JP4724626B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
KR101088115B1 (en) Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their devices, and their programs
JP4410833B2 (en) Springback generation cause analysis method, apparatus, program and recording medium
JP4418452B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
KR101368108B1 (en) Springback occurrence cause analyzing method, springback occurrence cause analyzing device, computer readable recording medium for recording springback occurrence cause analyzing program
JP4894294B2 (en) Press forming analysis method
JP4964988B2 (en) Molding simulation method, molding simulation apparatus, program, recording medium, and molding method based on simulation results
JP6060591B2 (en) Method and apparatus for identifying springback factor of press-formed product
JP4344375B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP5068783B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP4418451B2 (en) Springback countermeasure position specifying method, apparatus, and program thereof
JP4308239B2 (en) Springback influence display method, apparatus thereof, and program thereof
JP4418449B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
RU2424074C2 (en) Method of identifying elastic afterworking cause, method of displaying influence of elastic afterworking, method of identifying location of elastic afterworking cause, method of identifying position of elastic afterworking preventer, devices for their implementation and programmes
JP6303815B2 (en) Forging crack prediction method, forging crack prediction program, and recording medium
Haase et al. Enabling stamping processes through meticulous FE Modelling-segmented drawbeads and remesh criteria
Nakkhlai et al. NEURAL NETWORK APPLICATION FOR FINE-BLANKED EDGE QUALITY IN ROLLED STEEL SHEETS
JP2008149356A (en) Method and apparatus for designing highly rigid component, computer program, and computer readable recording medium

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080805

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20081114

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20081114

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20081208

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081216

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090216

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: 20090421

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090430

R151 Written notification of patent or utility model registration

Ref document number: 4308239

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120515

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130515

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130515

Year of fee payment: 4

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130515

Year of fee payment: 4

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130515

Year of fee payment: 4

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130515

Year of fee payment: 4

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140515

Year of fee payment: 5

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350