CN101511502A - Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their - Google Patents

Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method, their Download PDF

Info

Publication number
CN101511502A
CN101511502A CN200780032577.9A CN200780032577A CN101511502A CN 101511502 A CN101511502 A CN 101511502A CN 200780032577 A CN200780032577 A CN 200780032577A CN 101511502 A CN101511502 A CN 101511502A
Authority
CN
China
Prior art keywords
mentioned
springback
calculation process
zone
physical quantity
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.)
Granted
Application number
CN200780032577.9A
Other languages
Chinese (zh)
Other versions
CN101511502B (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 claimed from JP2006244351A external-priority patent/JP4418449B2/en
Priority claimed from JP2006269253A external-priority patent/JP4344375B2/en
Priority claimed from JP2006268855A external-priority patent/JP4308239B2/en
Priority claimed from JP2006271097A external-priority patent/JP4418452B2/en
Priority claimed from JP2006270693A external-priority patent/JP4418451B2/en
Priority to CN201410204754.9A priority Critical patent/CN103995927B/en
Application filed by Nippon Steel Corp filed Critical Nippon Steel Corp
Priority claimed from PCT/JP2007/067358 external-priority patent/WO2008026777A1/en
Publication of CN101511502A publication Critical patent/CN101511502A/en
Publication of CN101511502B publication Critical patent/CN101511502B/en
Application granted granted Critical
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82—Elements for improving aerodynamics

Landscapes

  • Shaping Metal By Deep-Drawing, Or The Like (AREA)

Abstract

A method of identification of a cause of occurrence of springback comprising a press forming analysis step of numerically analyzing forming conditions of press forming to obtain forming data of a press formed part, a processing step of processing at least one of a physical property value and physical property quantity data of part of the regions of said press formed part among forming data of said press formed part, and a springback value calculation step of calculating a springback value based on the results of said processing.

Description

Springback occurrence cause ad hoc approach, springback influence degree display method, springback occurrence cause portion ad hoc approach, resilience countermeasure position specifying method, their device and their program
Technical field
The present invention relates in occurrence cause ad hoc approach, influence degree display method, occurrence cause portion ad hoc approach, countermeasure position specifying method, their device and their program by the resilience that formed article by press working takes place during with member etc. with metal stamping shaping automobiles such as the aluminium headed by the steel such as steel plate, bar steel, magnesium, titaniums.
Background technology
A lot of automobiles such as door and bumper are to make by the drawing of steel plate with member, tame electric components, building materials etc.In recent years, improve,, just seeking to have the countermeasure that high-intensity steel make steel thin-walled property etc. by use in order to realize this requirement for the light-weighted requirement of these members.
Yet, be accompanied by the high strength of steel plate, article shape in order to ensure drawing, need strict control, one of more important project is in this management, the residual stress that produces in the steel plate during with drawing is a driving force, and the distortion that the strain of steel plate is partly carried out elastic recovery is resilience.
Particularly recently in order to cut down the exploitation man-hour and the cost of automobile etc., has the trend of design phase that beginning and design phase side by side discuss the manufacturing process of this formed parts, the change of the formed parts in the design phase takes place in the design alteration meeting in the design phase, therefore man-hour and the cost in the design phase of the manufacturing process of investigation formed parts becomes bigger problem in the development technology of automobile etc. and fee to develop.
Fig. 1 is the cross section ideograph to the formed parts of the countermeasure of resilience that shows in the past.(a) show shaping shape of products cross section, (b) show adopt with (a) shown in the mould of molded article homotype steel plate is carried out cold stamping and is shaped after, the resilience that molded article produces (c) shows the shaped cross of the mould that imagination resilience and having carried out is revised.Like this, obtain the molded article shown in (a), utilize the mould of " predictability " of the anticipation resilience shown in (c), obtain the countermeasure of desirable molded article with regard to promising.
Method as the mould of such anticipation resilience that is shaped, once proposed: utilize FInite Element to resolve residual stress by the steel plate of the lower dead center of mould punching press, (resilience is advanced in the distortion that produces by numeric value analysis because with the reverse residual stress of this residual stress; Spring forward) mould of shape is considered the method (" TOHKEMY 2003-33828 communique ", " テ of Mitsubishi Motors Network ニ カ Le レ PVC ユ-(126～131 pages of 2006 No.18) ") of the shaping of carrying out mould of resilience simply.
, design the linear problem of mould right and wrong of considering resilience fully by numeric value analysis, difficulty, so the method that has proposed in the extreme is to adopt FInite Element to come the method that the mould of considering resilience is simply formed eventually.Therefore, do not satisfy what kind of countermeasure of occasion needs of the allowance of resilience for utilizing this mould,, therefore still do not propose what solution because difficulty is resolved on numerical value ground.
Therefore,, take what kind of countermeasure, depend on technical staff's experience in order to obtain desirable molded article utilizing the occasion of considering the mould of resilience simply and not satisfying the allowance of resilience.As a result, the tentative error checking of the mould that need obtain by this forming process and actual steel plate.
In addition, once having proposed is not to mold shape, but the shape of steel, molded article is applied the correction of removing residual stress, reduces the method for resilience thus.
Fig. 2 is the stereogram of the previous methods at the illustration position of seeking the occurrence cause become the distortion that is caused by resilience.(a) show the shaping shape of products, (b) show the situation of cutting off a part 1 of removing goods (c) to be presented at the situation of having opened hole 2 on the goods, (d) show the situation of giving otch 3 part of goods.By when taking such countermeasure to observe the behavior of resilience, attempted reducing the countermeasure of resilience.
; even the residual stress that makes the reason that becomes resilience by the resilience happening part is taken some countermeasures reduces since cut off remove, perforate etc., the member self rigidity reduces; therefore also produce the problem that resilience takes place owing to the residual stress of a little, do not find out the underlying cause as yet.And in fact such countermeasure needs the test of adopting test mould and steel plate to carry out, therefore produces the man-hour of design phase and the problem that cost increases.
Summary of the invention
In view of the above-mentioned problems, problem of the present invention is, come the position that becomes springback occurrence cause of specific formed article by press working by numeric value analysis, and the proterties of this privileged site is carried out numeric value analysis, thus high efficiency and shorten investigation time of the manufacturing process of formed parts economically.
In addition, problem of the present invention is to provide by numeric value analysis that to keep product strength and make springback capacity be molded article below the permissible value.
In addition, problem of the present invention is to provide the countermeasure that makes the springback capacity minimizing at the position that becomes springback occurrence cause of formed article by press working by numeric value analysis.
In order to solve above problem, the invention provides a kind of springback occurrence cause ad hoc approach, this method has: the molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product; The physics value in parts among the one-tenth graphic data of drawing product, the drawing product zone and at least one of physical quantity data are carried out the calculation process step of calculation process; Calculate the springback capacity of springback capacity with result and calculate step based on calculation process.
Above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain, said method, a part of zone repeats above-mentioned calculation process step and springback capacity is calculated step by changing, and also can become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving zone, physics value and the physical quantity when becoming maximum carried out specific to springback capacity.
In addition, calculate in the step in above-mentioned calculation process step and springback capacity, can have following step: repeat by at least one and/or the calculation process that changes physics value and physical quantity, come springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving zone, physics value and the physical quantity when becoming maximum carried out specific.
And then, above-mentioned part zone is a plurality of zones, for each zone, also can carry out calculation process simultaneously, above-mentioned calculation process step also can be: divide the zone that springback capacity becomes minimum, about the size in the zone that is divided to division is each zone below the setting, and at least one of physics value and physical quantity carried out the step of calculation process.In addition, above-mentioned part zone can be more than one unit or unit of account section (calculation unit sections), and in addition, above-mentioned a part of zone also can be more than one point.
In addition, the invention provides a kind of springback occurrence cause specific device, it has: the molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product; Numeric value analysis becomes graphic data, calculates the resilience analysis unit of springback capacity; With at least a calculation process of carrying out, make the resilience analysis unit calculate arithmetic processing section based on calculation process result's springback capacity for the data of the physics value that becomes parts among the graphic data, the drawing product zone of drawing product and physical quantity.
Above-mentioned physics value and physical quantity also can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain, change a part of zone by above-mentioned arithmetic processing section and repeat the result of calculation process and calculating of springback capacity, can become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving zone, physics value and the physical quantity when becoming maximum carried out specific springback capacity.
In addition, above-mentioned springback occurrence cause specific device, repeat the result of calculation process and calculating of springback capacity by utilizing above-mentioned arithmetic processing section to change physics value and physical quantity and/or calculation process, can become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving zone, physics value and the physical quantity when becoming maximum carried out specific springback capacity.
And then, above-mentioned part zone is a plurality of zones, can carry out calculation process simultaneously for each zone, in addition, the springback occurrence cause specific device can utilize above-mentioned arithmetic processing section, divides the zone of springback capacity minimum, about the size in the zone that is divided to division is each zone below the setting, and at least one of physics value and physical quantity carried out calculation process.
In addition, in order to solve above-mentioned problem, the invention provides a kind of springback occurrence cause ad hoc approach, it has: the molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product; Calculate among the one-tenth graphic data decide the drawing product, about the curvature and/or the angle of drawing product, zone based on curvature and/or angular divisions formed products, about each zone of dividing, for carrying out the calculation process step of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity based on above-mentioned shaping data judging; With result, calculate the springback capacity of springback capacity and calculate step based on calculation process.
Above-mentioned calculation process step also can be carried out calculation process at least one at least one of physics value physical quantity variable of the point of the big relatively zoning of curvature and/or angle.In addition, above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
Above-mentioned springback occurrence cause ad hoc approach is calculated in the step and can be had in calculation process step and springback capacity: repeats by changing a part of zone, come zone, physics value and physical quantity that springback capacity becomes hour are carried out specific step.
And then, calculate in the step in above-mentioned calculation process step and springback capacity, can have: repeat by at least one and/or the calculation process that changes physics value and physical quantity, come zone, physics value and physical quantity that springback capacity becomes hour are carried out specific step.
In addition, the invention provides a kind of springback occurrence cause specific device, it has: the molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product; Calculate among the one-tenth graphic data decide the drawing product, about the curvature and/or the angle of drawing product, zone based on curvature and/or angular divisions formed products, for each zone of dividing, to carrying out the arithmetic processing section of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity variable based on above-mentioned shaping data judging; With result, calculate the resilience analysis unit of springback capacity based on above-mentioned calculation process.
Physics value and physical quantity in the big relatively zoning of above-mentioned curvature and/or angle also can be at least one the physics value and the physical quantitys of point of the big relatively zoning of curvature and/or angle.In addition, above-mentioned physics value and physical quantity also can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
And then, above-mentioned springback occurrence cause specific device, by utilizing above-mentioned arithmetic processing section, change a part of zone and repeat the result of calculation process and calculating of springback capacity, can carry out specific to zone, physics value and physical quantity that springback capacity becomes hour, and, repeat the result of calculation process and calculating of springback capacity by changing physics value and physical quantity and/or calculation process, can carry out specific zone, physics value and physical quantity that springback capacity becomes hour.
In addition, in order to solve above-mentioned problem, the invention provides a kind of springback influence degree display method (program), it has: in the computer that shows springback influence degree, the molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product; At least one of each regional physics value among the one-tenth graphic data of drawing product, that divided the drawing product and physical quantity data carried out the calculation process step (process) of calculation process; Based on the result of calculation process, calculate the springback capacity of springback capacity and calculate step (process); Step display (process) with the springback capacity of calculating for each demonstration in each zone.
Calculation process step (process) can be carried out calculation process to the physics value of a part of point in each zone among the one-tenth graphic data of drawing product, that divided the drawing product and at least one of physical quantity data.In addition, physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
Step display (process) also can show springback capacity for each profile in each zone, and profile shows the value that the springback capacity of calculating obtains divided by each regional area.In addition, step display (process), the springback capacity that can calculate for each demonstration in each zone divided by the representative length of drawing product, represent width, represent height, represent any of thickness of slab, tensile strength and the value that obtains.And the springback capacity that step display (process) also can be calculated for each demonstration in each zone is divided by drift speed or blank holder pressure and the value that obtains.
In addition, the invention provides a kind of springback influence degree display unit, it has: the molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product; Numeric value analysis becomes graphic data, calculates the resilience analysis unit of springback capacity; Each regional physics value among the one-tenth graphic data of drawing product, that divided the drawing product and at least one of physical quantity data are carried out calculation process, make the resilience analysis unit calculate arithmetic processing section based on calculation process result's springback capacity; Display part with the springback capacity of calculating for each demonstration in each zone.
Arithmetic processing section can be carried out calculation process to the physics value of a part of point in each zone among the one-tenth graphic data of drawing product, that divided the drawing product and at least one of physical quantity data, and physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
Display part can show springback capacity for each profile in each zone, and the value profile that the springback capacity of calculating can be obtained divided by each regional area is shown in display part.Display part also can for each zone each with the springback capacity of calculating divided by the representative length of drawing product, represent width, represent height, represent any of thickness of slab, tensile strength and the value that obtains is shown in display part.And the springback capacity that display part also can be calculated for each demonstration in each zone is divided by drift speed or blank holder pressure and the value that obtains.
Program of the present invention, except the program of above-mentioned explanation, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product drawing analysis program, numeric value analysis this become graphic data, calculating resilience analysis program, the post processor that profile shows this springback capacity of springback capacity and can use existing program to carry out the springback influence degree display routine of these programs and the defeated input and output of data, can be to make computer carry out the program of following process: the process that is obtained into graphic data by the drawing analysis program; Each regional physics value among the one-tenth graphic data, that divided the drawing product and at least one of physical quantity data are carried out the calculation process process of calculation process; Export the result of above-mentioned calculation process in the resilience analysis program process; Make above-mentioned post processor profile show that the resilience analysis program carries out aforementioned calculation for each of each zone and the profile procedure for displaying of the springback capacity that obtains.Existing resilience analysis program is based on calculating springback capacity by the calculation process data of program output of the present invention.
In addition, in order to solve above-mentioned problem, springback occurrence cause portion ad hoc approach of the present invention (program: following bracket in corresponding with the invention of program) has: the molding condition of numeric value analysis drawing obtains the drawing analyzing step (process) of the one-tenth graphic data of drawing product; Based on the one-tenth graphic data of drawing product, calculate after the resilience a plurality of zones separately physics value and the 1st springback capacity of the data of physical quantity calculate step (process); Occasion at least one zone bigger than setting with physics value and physical quantity is carried out the calculation process step (process) of calculation process at least one of this regional physics value and physical quantity data; With result based on calculation process, further calculate after the resilience a plurality of zones separately physics value and the 2nd springback capacity of the data of physical quantity calculate step (process).
In addition, above-mentioned calculation process step (process), occasion in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the resilience and physical quantity data, also can carry out calculation process at least one of this regional physics value and physical quantity data, above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
In addition, springback occurrence cause portion specific device of the present invention provides: the molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product; Numeric value analysis becomes graphic data, calculates the resilience analysis unit of springback capacity; Occasion at least one zone bigger than setting with physics value and physical quantity, at least one of carrying out this regional physics value and physical quantity data carried out the calculation process of calculation process, makes the resilience analysis unit further calculate the arithmetic processing section of a plurality of zones springback capacity separately after the resilience based on the result of calculation process.
And then, occasion in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the resilience and physical quantity data, arithmetic processing section also can be carried out calculation process at least one of this regional physics value and physical quantity data, and above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
Program of the present invention, except the program of above-mentioned explanation, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product and a plurality of zone physics value separately and the drawing analysis program of the data of physical quantity, this becomes graphic data and a plurality of zone physics value separately and the data of physical quantity numeric value analysis, calculate a plurality of zones physics value separately and the resilience analysis program of the data of physical quantity after springback capacity and the resilience, can use existing program to carry out the springback occurrence cause portion specific program of the defeated input and output of these programs and data, can be to make computer carry out the program of following process: obtain a plurality of zones physics value separately and the process of the data of physical quantity after the above-mentioned resilience by the resilience analysis program; Occasion at least one zone bigger than setting with above-mentioned physics value and physical quantity is carried out the calculation process process of calculation process at least one of this regional physics value and physical quantity data; The result of above-mentioned calculation process is outputed to process in the resilience analysis program.Existing resilience analysis program based on the calculation process data by program output of the present invention, calculates springback capacity again.
In addition, in order to solve above-mentioned problem, springback occurrence cause portion ad hoc approach of the present invention (program: following bracket in corresponding with the invention of program) has: numeric value analysis the 1st molding condition, calculate the 1st drawing analyzing step (process) of the 1 one-tenth graphic data; At least one the 2nd different molding condition of numeric value analysis the 1st molding condition is calculated the 2nd drawing analyzing step (process) of the 2 one-tenth graphic data; The a plurality of zones with the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, the occasion in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out the step (process) of calculation process; With result, calculate the springback capacity of springback capacity and calculate step (process) based on calculation process.
In addition, above-mentioned molding condition can be the shape of steel plate and proterties, mold shape, punching press condition, and above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
In addition, springback occurrence cause portion specific device of the present invention provides: numeric value analysis the 1st molding condition, calculate the 1 one-tenth graphic data, at least one the 2nd different molding condition of numeric value analysis the 1st molding condition is calculated the drawing analysis unit of the 2 one-tenth graphic data; Numeric value analysis becomes graphic data, calculates the resilience analysis unit of springback capacity; The a plurality of zones with the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, the occasion in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out calculation process, make the resilience analysis unit calculate arithmetic processing section based on calculation process result's springback capacity.
In addition, above-mentioned molding condition can be the shape of steel plate and proterties, mold shape, punching press condition, and above-mentioned physics value and physical quantity can be the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
Program of the present invention, except the program of above-mentioned explanation, be to calculate the one-tenth graphic data of drawing product and a plurality of zone physics value separately and the drawing analysis program of the data of physical quantity, this becomes graphic data and a plurality of zone physics value separately and the data of physical quantity with numeric value analysis, calculate a plurality of zones physics value separately and the numeric value analysis program of the data of physical quantity after springback capacity and the resilience, and can use existing program to carry out the springback occurrence cause portion specific program of these programs and the defeated input and output of data, can be to make computer carry out the program of following process: the process that obtains the 1 one-tenth graphic data by the drawing analysis program; By at least one the 2nd different molding condition of drawing analysis program numeric value analysis the 1st molding condition, obtain the process of the 2 one-tenth graphic data; The a plurality of zones with the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, the occasion in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of above-mentioned the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out the process of calculation process; Output to process in the resilience analysis program with result with above-mentioned calculation process.Existing resilience analysis program is based on calculating springback capacity by the calculation process data of program output of the present invention.
In addition, in order to solve above-mentioned problem, resilience countermeasure position specifying method of the present invention has: the molding condition of numeric value analysis drawing product, calculate the drawing analyzing step (process) of the one-tenth graphic data of drawing product; At least a portion zones among the one-tenth graphic data of selected drawing product, drawing product, generate be replaced as with the variform shape of this selection area and/or, with the stress distribution different with the stress distribution of this selection area replace be replaced as graphic data be replaced as graphic data generation step (process); Be replaced as graphic data with numeric value analysis, calculate the springback capacity of springback capacity and calculate step (process).
Also can change the position and the quantity of above-mentioned selection area, repeat to be replaced as graphic data generation step and springback capacity and calculate step (process), can further have and judge whether springback capacity is the following step (process) of setting.
In addition, also can at least one predetermined data displacement with this displacement of the stress distribution of the shape in one or more zones of the following springback capacity of setting and this selection area will be given, and then change the position or the quantity of selection area, repeat above-mentionedly to be replaced as that graphic data generates step (process) and above-mentioned springback capacity is calculated step (process).
And then, can further have the step (process) that the one-tenth graphic data of drawing product is divided with the zone of regulation, be replaced as graphic data and generate step (process), can be for the zone of dividing with the zone of regulation whole, calculate and be replaced as graphic data.In addition, the data of afore mentioned rules are selected the zone for each, and at least one of shape and stress distribution can be different.
In addition, resilience countermeasure of the present invention position specific device has: the molding condition of numeric value analysis drawing product, calculate the drawing analysis unit of the one-tenth graphic data of drawing product; At least a portion zones among the one-tenth graphic data of selected drawing product, drawing product, generation with the stress distribution of the shape of this selection area and this selection area at least one with predetermined data replace form be replaced as graphic data be replaced as the graphic data generating unit; Be replaced as graphic data with numeric value analysis, calculate the resilience analysis unit of springback capacity.
Change the position and the quantity of above-mentioned selection area, be replaced as the graphic data generating unit and calculate and be replaced as graphic data, and the resilience analysis unit repeats springback capacity and calculates also passable.In addition, be replaced as the graphic data generating unit and can judge also whether springback capacity is below the setting.
In addition, at least one predetermined data displacement with this displacement of the stress distribution of the shape in one or more zones of the springback capacity below the setting and this selection area will be given, and then the position or the quantity of change selection area, be replaced as the graphic data generating unit and calculate the above-mentioned graphic data that is replaced as, and above-mentioned resilience analysis unit repetition springback capacity is calculated also passable.
Be replaced as the graphic data generating unit, can divide the one-tenth graphic data of drawing product with the zone of regulation, and whole for the zone of the zone division of regulation, calculate and be replaced as graphic data.In addition, the data of afore mentioned rules are selected the zone corresponding to each, and at least one of shape and stress distribution can be different.
Program of the present invention, except the program of above-mentioned explanation, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product drawing analysis program, numeric value analysis this become graphic data to calculate the resilience analysis program of springback capacity and can use existing program to carry out the resilience countermeasure position specific program of these drawing analysis programs and resilience analysis program and the defeated input and output of data, be to make computer carry out the program of following process: the process that obtains into graphic data by the drawing analysis program; At least a portion zones among the one-tenth graphic data of selected drawing product, the drawing product, calculate with predetermined data replace that at least one of stress distribution of the shape of this selection area and this selection area form be replaced as graphic data be replaced as the graphic data generative process; Export the above-mentioned graphic data that is replaced as in the above-mentioned resilience analysis program process.Existing resilience analysis program is based on calculating springback capacity by the graphic data that is replaced as of program output of the present invention.
According to the present invention, a kind of method is provided, change becomes the privileged site of analysis object as the reason of the resilience generation of formed article by press working, and, physics value physical quantity by one side this privileged site of numerical operation, on one side springback capacity is minimized, come to carry out specific to the cause portion that resilience takes place, and making correctly derives the physics value of this cause portion physical quantity and becomes possibility, or can confirm them by the tentative mistake of resolving, economical thus and shorten investigation time of the manufacturing process of formed parts expeditiously.
In addition, according to the present invention, a kind of method is provided, the reason that will take place as the resilience of formed article by press working becomes the privileged site of analysis object, each regional physics value of having divided the drawing product and at least one of physical quantity data are carried out calculation process, result based on this calculation process, springback capacity is carried out profile for each of each zone to be shown, make the tentative mistake of parsing on computers become possibility thus, economical and shorten the specific of cause portion that resilience takes place expeditiously, and economical and shorten investigation time of the manufacturing process of molded article expeditiously.
In addition, according to the present invention, a kind of method is provided, become the selection zone of analysis object for the reason that takes place as the resilience of formed article by press working, by being replaced as and this different shape and stress distribution in selection zone, and springback capacity minimized, and make the cause portion that resilience is taken place carry out specific and provide countermeasure to become possibility at this cause portion, economical thus and shorten investigation time of the manufacturing process of formed parts expeditiously.
In addition, the present invention by making privileged site and this privileged site physics value and physical quantity any changes at least, do not cut off for the molded article of reality remove, perforate just can obtain that to make the resilience variable quantity be privileged site and physics value and physical quantity below the permissible value.Therefore, the molded article of Xie Xiing in this wise, can confirm make the resilience variable quantity be below the permissible value in, kept the product property of rigidity etc., the countermeasure with perforate etc. is removed in the cut-out to the molded article privileged site of therefore can not need sacrifice product property for the resilience variable quantity of keeping the desired regulation of molded article.
In addition, the present invention does not cut off by the molded article for reality and removes, just at least one at least one or both sides that select different shape of regional replacement and stress distribution, can make springback capacity is below the permissible value in perforate.Therefore, the molded article of Xie Xiing in this wise, can confirm make springback capacity be below the permissible value in, kept the product property of rigidity etc., the countermeasure with perforate etc. is removed in the cut-out to the molded article privileged site of therefore can not need sacrifice product property for the springback capacity of keeping the desired regulation of molded article.
Description of drawings
Following with reference to description of drawings the present invention.
Fig. 1 is the cross section ideograph to the formed parts of the countermeasure of resilience that shows in the past.
Fig. 2 is the stereogram of the previous methods at the illustration position of seeking the occurrence cause become the distortion that is caused by resilience.
Fig. 3 is the function constitution map of the springback occurrence cause portion specific device that relates to of an embodiment of the invention.
Fig. 4 is the profile diagram of the stress distribution of the drawing analysis result that relates to of expression example of the present invention.
Fig. 5 is the profile diagram of the curvature of each unit of relating to of expression example of the present invention.
Fig. 6 is that the big unit of curvature that makes that example of the present invention relates to is the profile diagram of non-demonstration.
Fig. 7 is the figure that has carried out regional division based on curvature that expression example of the present invention relates to.
Fig. 8 is that the stress that makes part zone that example of the present invention relates to is zero profile diagram.
Fig. 9 is the figure of the coordinate system that uses in FInite Element that shows that an embodiment of the invention relate to.
Figure 10 is the stereogram of the regioselective example of the machined object that shows that an embodiment of the invention relate to.
Figure 11 is the flow chart of the springback occurrence cause portion particular procedure that relates to of an embodiment of the invention.
Figure 12 is the flow chart of the springback occurrence cause portion particular procedure that relates to of an embodiment of the invention.
Figure 13 is the flow chart of the springback occurrence cause portion particular procedure of the resilience cause portion specific device that relates to of an embodiment of the invention.
Figure 14 is the function constitution map of the springback occurrence cause portion specific device that relates to of an embodiment of the invention.
Figure 15 is the flow chart of the springback occurrence cause portion particular procedure that relates to of one embodiment of the invention.
Figure 16 is the function constitution map of the resilience countermeasure position specific device that relates to of an embodiment of the invention.
Figure 17 is the flow chart of the springback occurrence cause portion particular procedure that relates to of an embodiment of the invention.
Figure 18 is the hardware structure diagram of the springback occurrence cause portion particular procedure that relates to of an embodiment of the invention.
Figure 19 (a) and (b) are figure of the hat section shaped member that becomes the object that resilience resolves that shows that embodiments of the invention relate to.
Figure 20 (a) and (b) are the sectional position of the hat section shaped member that shows that embodiments of the invention relate to and the figure of windup-degree.
Figure 21 is the stereogram of the specific region A3 that shows that embodiments of the invention relate to.
Figure 22 shows the specific region A4 that embodiments of the invention relate to, the stereogram of A5.
Figure 23 is the specific region A6～A9 stereogram that shows that embodiments of the invention relate to.
Figure 24 (a) is the stereogram of specific region A10～A14 of showing that embodiments of the invention relate to, (b) is the sectional view of specific region A10～A14 of showing that embodiments of the invention relate to.
Figure 25 is the stereogram of specific region A15～A17 of showing that embodiments of the invention relate to.
Figure 26 is the stereogram of specific region A18～A21 of showing that embodiments of the invention relate to.
Figure 27 is the figure of the windup-degree of the specific region A3 that shows that embodiments of the invention relate to.
Figure 28 is the figure of the windup-degree of the specific region A4, the A5 that show that embodiments of the invention relate to.
Figure 29 is the figure of the windup-degree of specific region A6～A9 of showing that embodiments of the invention relate to.
Figure 30 is the figure of the windup-degree of specific region A10～A14 of showing that embodiments of the invention relate to.
Figure 31 is the figure of the windup-degree of specific region A15～A17 of showing that embodiments of the invention relate to.
Figure 32 is the figure of the windup-degree of specific region A18～A21 of showing that embodiments of the invention relate to.
Figure 33 is the stereogram of the specific region A22 that shows that embodiments of the invention relate to.
The figure of the windup-degree when Figure 34 is the change operation method of the specific region A22 that shows that embodiments of the invention relate to.
Figure 35 is the figure of the windup-degree of specific region A23～A24 of showing that embodiments of the invention relate to.
Figure 36 (a) and (b) are the sectional position of the hat bending part that shows that embodiments of the invention relate to and the figure of opening.
Figure 37 is the result's of the investigation example 9 that shows that embodiments of the invention relate to figure.
Figure 38 is the result's of the investigation example 10 that shows that embodiments of the invention relate to figure.
Figure 39 is the stereogram and the sectional view of the hat section shaped member that shows that embodiments of the invention relate to.
Figure 40 is the sectional position of the hat section shaped member that shows that embodiments of the invention relate to and the figure of windup-degree.
Figure 41 is the stereogram of specific region A35～A37 of showing that embodiments of the invention relate to.
Figure 42 is the figure of the windup-degree (all points) of specific region A35～A37 of showing that embodiments of the invention relate to.
Figure 43 is the figure of the windup-degree (point at the most surperficial and back side of thickness of slab direction) of specific region A35～A37 of showing that embodiments of the invention relate to.
Figure 44 is the stereogram and the sectional view of the example of the machined object that shows that one embodiment of the invention relate to.
Figure 45 is the stereogram of the example divided of the zone of the machined object that shows that one embodiment of the invention relate to.
Figure 46 is the sectional position of the machined object that shows that one embodiment of the invention relate to and the stereogram and the sectional view of focus.
Figure 47 is the stereogram that each region contour of showing the machined object that one embodiment of the invention are related to shows the example of springback capacity.
Figure 48 is the stereogram that each region contour of showing the machined object that one embodiment of the invention are related to shows the example of springback influence degree.
Figure 49 is the stereogram and the sectional view of the example of the machined object that shows that one embodiment of the invention relate to.
Figure 50 is the sectional position of the hat section shaped member that shows that embodiments of the invention relate to and the figure of windup-degree.
Figure 51 is the figure of the stress distribution after the resilience of the machined object that shows that one embodiment of the invention relate to.
Figure 52 is the figure in the zone that shows that the stress distribution after the resilience of the machined object that relates to based on one embodiment of the invention is selected.
Figure 53 be comparison one embodiment of the invention relate to do not carry out calculation process the time springback capacity and the figure of the springback capacity when having carried out calculation process.
Figure 54 is the sectional position of the hat section shaped member that shows that embodiments of the invention relate to and the figure of flange spring amount.
Figure 55 is before the resilience of the machined object that shows that one embodiment of the invention relate to and the figure of the stress distribution after the resilience.
Figure 56 is the distribution of the stress difference before and after the resilience of the machined object that shows that one embodiment of the invention relate to and based on the figure in the zone of stress difference component selections.
Figure 57 be comparison one embodiment of the invention relate to do not carry out calculation process the time springback capacity and the figure of the springback capacity when having carried out calculation process.
Figure 58 is the stereogram and the sectional view of the example of the machined object that shows that one embodiment of the invention relate to.
Figure 59 is the sectional position of the machined object that shows that one embodiment of the invention relate to and the stereogram of windup-degree.
Figure 60 is the figure that shows that the quantity of state after the drawing of the 1st machined object that one embodiment of the invention relate to and the 2nd machined object distributes.
Figure 61 is the figure that shows the quantity of state difference after the drawing of the 1st machined object that one embodiment of the invention relate to and the 2nd machined object.
Figure 62 is the regioselective figure that shows the quantity of state difference that relates to based on one embodiment of the invention.
Figure 63 be show that one embodiment of the invention relate to do not carry out calculation process the time springback capacity and the figure of the comparison of the springback capacity when having carried out calculation process.
Figure 64 is the figure of the hat section shaped member that becomes the object that resilience resolves that shows that embodiments of the invention relate to.
Figure 65 is the figure of the data drawing list of the rib (bead) of the displacement that shows that embodiments of the invention relate to and stress distribution.
Figure 66 is the figure that the zone of the processed body that shows that embodiments of the invention relate to is divided.
Figure 67 is the figure that has disposed shape and stress on the appointed area of the processed body that embodiments of the invention relate to region-wide.
Figure 68 has carried out the figure that defines to sectional position and torsional capacity that embodiments of the invention relate to.
Figure 69 is that the appointed area of the processed body that relates at embodiments of the invention has disposed the figure in rib and zone.
Figure 70 is that the appointed area of the processed body that relates at embodiments of the invention has disposed the figure in rib and zone.
Figure 71 is the figure that has disposed 1 rib etc. of the springback capacity minimum that embodiments of the invention relate to.
Figure 72 is the figure that has disposed the 2nd rib etc. on the basis of 1 configuration such as rib of the springback capacity minimum that embodiments of the invention relate to.
Figure 73 is the figure that has disposed the 2nd rib etc. on the basis of 1 configuration such as rib of the springback capacity minimum that embodiments of the invention relate to.
Figure 74 is the figure that the appointed area of the processed body that relates at embodiments of the invention has disposed n rib etc.
Figure 75 is the figure that the appointed area of the processed body that relates at embodiments of the invention has disposed actual rib.
The specific embodiment
Fig. 3 shows the function constitution map of the springback occurrence cause portion specific device that an embodiment of the invention relate to.
Springback occurrence cause portion specific device 10A has: drawing analysis unit 11, resilience analysis unit 12, physics value physical quantity arithmetic processing section 14, file storage part 16, molding condition input part 18 and springback capacity output screen 19.
Molding condition input part 18 is input parts of the input shape data (thickness of slab, length, width, curvature, strain etc.), the proterties (shapes such as materials such as intensity, percentage elongation, thickness of slab) that become the steel plate of analysis object in drawing analysis unit 11 described later and resilience analysis unit 12, mold shape (punch die (mould) and punch shape, curvature, diameter, gap, lubricating condition), punching press condition molding conditions such as (blank holder load, liner (pad) load, rib tension force, stamping pressure, temperature).Perhaps, also can set data area in be shaped resolving, the data area in the physics value physical quantity arithmetic processing section 14, the zoning when carrying out screen display etc. respectively by springback capacity output screen 19.
Drawing analysis unit 11, will be by the molding condition of molding condition input part 18 input as input information, use elasto-plastic finite element method, rigid-plastic finite element method, a step FInite Element, boundary element method, elementary parsing etc., obtain numeric value analysis as stress, strain and the thickness of slab etc. of the steel plate of machined object etc. through drawing.Drawing analysis unit 11 is with this form output numeric value analysis result of distribution of variable such as the thickness of slab of machined object, the component value of stress, the component value of strain and this variable.This exports data, for example exports resilience analysis unit 12, physics value physical quantity arithmetic processing section 14 described later and file storage part 16 to as file " P org.k ".
The numeric value analysis of this drawing analysis unit 11, use FInite Element (for example existing software PAM-STAMP, LS-DYNA, Autoform, OPTRIS, ITAS-3D, ASU/P-FORM, ABAQUS, ANSYS, MARC, HYSTAMP, Hyperform, SIMEX, Fastform3D, Quikstamp), based on the shape of products data (thickness of slab that carries out drawing, length, width, curvature, strain etc.) and the proterties (intensity of the metallic plate that uses, materials such as percentage elongation, shapes such as thickness of slab), then set mold shape (punch die and punch shape if desired, curvature, diameter, the gap, lubricating condition), punching press condition (temperature, pressure) molding condition such as, thereby form parsing, and, the stress after can numerically obtaining being shaped, the distribution of strain value.
In addition, drawing analysis unit 11 can use the analysis result of the stress distribution that obtained by FInite Element, curvature etc. the poster processing soft profile to be shown in springback capacity output screen 19.
Resilience analysis unit 12 is used the output data file " P org.k " of drawing analysis unit 11 and output data file " P rem.casen.k " the conduct input data of physics value physical quantity arithmetic processing section 14 described later, carries out resilience and resolves.So-called resilience is resolved, be based on variable and variable distributions such as the thickness of slab output result, machined object as drawing analysis unit 11, the component value of stress, the component value of strain, calculating by elastic finite, elasto-plastic finite element method, a step FInite Element, elementary parsing etc. remove the lotus process comes the springback capacity of numeric value analysis as the deflection that machined object produced.This springback capacity is to adopt FInite Element etc. to divide machined object, to calculating each unit of three-dimensional data coordinate.Moreover, as the springback capacity of the deflection that machined object produced, have: the curvature that the angle that a plurality of faces arbitrarily in the maximum displacement point in the deflection of the arbitrfary point of machined object, the appointed area of machined object or the deflection of least displacement point, the machined object or line constitute mutually or the face arbitrarily of machined object or line constitute etc.
Resilience analysis unit 12, use FInite Element (for example existing software PAM-STAMP, LS-DYNA, Autoform, OPTRIS, ITAS-3D, ASU/P-FORM, ABAQUS, ANSYS, MARC, HYSTAMP), input is implemented resilience and is resolved by " P org.k " interior stress distribution of being put down in writing that drawing analysis unit 11 obtains.Resilience in software is calculated, and for example can adopt the basic formula of having followed the 71st～127 page of limited elastic-plastic deformation of being put down in writing of " non-linear FInite Element " (コ ロ Na society, distribution on December 20th, 1994) and the content of discretization method to calculate.But the calculating of resilience can be that elasticity is resolved, and also can be that elastoplasticity is resolved.
Output data as the resilience analysis result of resilience analysis unit 12, shapes such as the strain when comprising springback capacity, resilience, proterties, stress etc., export springback capacity output screen 19 to, and, export resilience analysis unit 12, physics value physical quantity arithmetic processing section 14 described later and file storage part 16 to as numeric value analysis result output data file " SB org.k " that obtains by input data file " P org.k " or the numeric value analysis result output data file " SB rem.casen.k " that obtains by " P rem.casen.k ".
Physics value physical quantity arithmetic processing section 14, the output data file " P org.k " of input drawing analysis unit 11, numeric value analysis result output data file " SB org.k " or " the SB rem.casen.k " of resilience analysis unit 12, carry out calculation process, as its result, " Prem.casen.k " exported to resilience analysis unit 12.Moreover said here calculation process is meant with at least more than one coefficient of above-mentioned variable doubly, makes to become to comprise zero steady state value, carry out arithmetic, calculate, be replaced as inconstant value arbitrarily based on function.Physics value physical quantity arithmetic processing section 14 by such calculation process, can be differentiated variable and privileged site that springback capacity is reduced.
Physics value physical quantity arithmetic processing section 14, when the data that obtain by the input data as the formed article by press working of one of processed body, find privileged site for a part, must carry out the data of this formed article by press working are divided into the processing in a plurality of zones from this formed article by press working.One of this region partitioning method is to utilize impartial size to divide the zone of processed body.
Impartial size, the preferred countermeasure of considering at privileged site.That is,, then can not make full use of analysis result effectively sometimes if, also do not take atomic little zone to the effective countermeasure of privileged site even find privileged site.For example can enumerate: in the countermeasure at privileged site is to apply under the situation of rib of 20mm * 100mm, and the zoning is that 20mm is square.
A method as the zoning that is used for definite formed products also has the method for determining the zoning based on curvature, unit.
Drawing analysis unit 11 by numeric value analysis, obtains the coordinate figure of each node as the geometric information of the machined object after be shaped resolving, and can calculate the curvature of each unit, the angle between the unit objectively.Under the situation of drawing product, the machined object after the distortion in the crest line portion of the bending of drift shoulder R, punch die shoulder R etc., is compared very large curvature becoming with the vertical direction of crooked crest line with other positions.
At this, calculate the maximum curvature of each unit, if making the unit that becomes the above curvature of certain threshold value is non-demonstration, then can differentiate separation as each zone of drift not being takeed on R, punch die shoulder R position (web face, longitudinal wall part, flange part) connection in addition.
Fig. 4～Fig. 8 shows the example that is resolved to the process of calculation process from drawing.
Fig. 4 is the drawing analysis result, and stress distribution is shown by profile diagram.The stress level of ladder side sidewall (absolute value) is higher as can be known.Fig. 5 demonstrates each unit, the curvature as the drawing analysis result, and white～grey display part is represented the position that curvature is big.Fig. 6 is that to make the big unit of curvature be the figure of non-demonstration.Because at end difference is very large curvature, if therefore making the unit for the above curvature of threshold value in the longitudinal wall part is non-demonstration, then can be separated into the crest line portion of ladder is a plurality of zones on border, by extracting out or selecting other zone, can select the zone based on character shape.As shown in Figure 6, can be divided into stepped sidewall 3 faces, webs, the longitudinal wall part (hiding among the figure) that does not have ladder, each flange part of both sides.Regional A105～A107 shown in Figure 7 is each zone of dividing about according to Fig. 6, is judged to be the zone of carrying out calculation process based on the high zone of one-tenth graphic data (stress distribution) counter stress shown in Figure 4.In determination step,, whether carry out computing described later based on the shaping data judging about zone based on curvature and/or angular divisions.For example, judge according to the stress level after each regional drawing parsing of dividing.
Fig. 8 is about regional A105～A107, and it is the state of zero calculation process that whole unit in each zone has been carried out for example making stress.
At this, used the maximum curvature of unit, even but use the angle between the unit also can carry out same zone division.
In addition, consider the structural elements that does not have the hat section of ladder shape,, can select the longitudinal wall part of web face, both sides, the flange part of both sides as each separated zone.In addition, when being not the plane,, can further select distinctive zone by repeating same method in the zone of selecting.
Like this, physics value physical quantity arithmetic processing section 14 can be determined the zoning based on curvature, unit according to the analysis result of drawing analysis unit 11.
Resilience analysis unit 12, adopt FInite Element (for example existing software PAM-STAMP, LS-DYNA, Autoform, OPTRIS, ITAS-3D, ASU/P-FORM, ABAQUS, ANSYS, MARC, HYSTAMP), input is implemented resilience and is resolved by " P org.k " interior stress distribution of being put down in writing that drawing analysis unit 11 obtains.Resilience in software is calculated, and for example can adopt the basic formula of having followed limited elastic-plastic deformation, the content of discretization method to calculate.But the calculating of resilience can be that elasticity is resolved, and also can be that elastoplasticity is resolved.
In addition, resilience analysis unit 12 can use the resilience analysis result that is obtained by FInite Element the poster processing soft profile to be shown in springback capacity output screen 19.In addition, in physics value physical quantity arithmetic processing section 14, with springback capacity by with the representative length of zoning area, drawing product, represent width, represent highly, represent removal such as thickness of slab, tensile strength and standardization, make variable and the easier differentiation of influence springback degree of having carried out division arithmetic, carry out data output, use this output data, utilize the poster processing soft can carry out profile and show.By such standardization, can visually catch the physics value physical quantity of processed body and the relation of springback capacity easily.
And, in physics value physical quantity arithmetic processing section 14, springback capacity is removed by using drawing conditions such as drift speed or blank holder pressure, obtain normalized value, utilize the poster processing soft to carry out profile and show, also can visually catch the relation of drawing condition and springback capacity thus easily.Such springback capacity standardization and the profile of numerical value thereof show, can be economical and shorten the specific of cause portion that resilience takes place expeditiously, and can be economical and shorten investigation time of the manufacturing process of molded article expeditiously.
File storage part 16 is the storage parts that are used to preserve data file " P org.k " as the output result of drawing analysis unit 11, resilience analysis unit 12, physics value physical quantity arithmetic processing section 14 described later, " SB org.k ", " P rem.casen.k ", " SB rem.casen.k ", " P trim.casen.k " etc., the occasion of direct input and output between drawing analysis unit 11, resilience analysis unit 12 and physics value physical quantity operational part 14 such as these data files, this document storage part 16 may not need.
Physics value physical quantity arithmetic processing section 14, select the part zone in the output data file " P org.k " of drawing analysis unit 11, at least more than one variable of the component value of the component value of the thickness of slab in this part zone, coefficient of elasticity, plastic coefficient, stress, strain is implemented calculation process.Then, the output data file " P rem.case1.k " after the generation calculation process exports resilience analysis unit 12 to.Sending and receiving of this data file can repeat, and can be stored in file storage part 16 successively, in addition, can not be with document form also, but carries out the transfer of data with the form of the data input and output of implementation or thread (thread).
The resilience privileged site that uses in FInite Element or be used for the area data that reason is judged can be: the unit of account section of the aggregate that the more than one unit in the FInite Element, the tiny area that comprises a plurality of unit or unit are continuous etc., as with a part of point among the point of the unit of account calculated stress of the numerical integration method of FInite Element, strain (comprise a part point and the interior part of face on the thickness of slab direction put the two) etc.So-called " more than one unit " is meant the unit as the regional dividing unit of FInite Element, and so-called " unit of account section " is meant one of tiny area that becomes unit of account in the elementary parsing or continuous aggregate.In addition, so-called " point " is the point that generally adopts the approximate integration that FInite Element carries out.So-called " a part of point " is meant the part among the point of unit of account calculated stress with the numerical integration method of FInite Element, strain, comprise the interior part point of part point on the thickness of slab direction and face both.
In addition, so-called " thickness of slab; coefficient of elasticity; plastic coefficient; the component value of stress; the variable of at least more than one of the component value of strain ", be meant the thickness of slab of the part of machined object, coefficient of elasticity (Young's modulus, Poisson's ratio), plastic coefficient (yield stress, plastic coefficient, work hardening index), at three-dimensional system of coordinate (x for resolving, y, z) time, be meant x, y, the stress of z direction or strain (3 components separately), the xy plane, the yz plane, among 18 factors of total of the shear stress on zx plane or shear strain (3 components separately) at least more than one.
Fig. 9 is the figure that is presented at the coordinate system that uses in the FInite Element.(a) represent three-dimensional global coordinate, (b) the expression local coordinate system.
At this moment, shown in Fig. 9 (a), except the component value on (x, y, z) reference axis, for example, shown in Fig. 9 (b), also comprise 12 components of total of the plane of unit 31 being regarded as stress in the local coordinate system (x1, x2, x3) on xy plane, strain, shear stress, shear strain.In addition, also comprise each components of stress of using equivalent stress, quite plastic strain, elastic energy (elastic work), plasticity can (plastic work done) to wait, variable that the components of strain obtain, increase the component value of increase portion-form calculating that part, strain increase each component value of part etc. with stress.
Usually, in the FInite Element that with the steel plate is analysis object was resolved, the physics value physical quantity of this steel plate can be divided into the square finite element of 2mm～4mm., the dividing unit of unit is to be approximately constant length with the physics value physical quantity to divide, and it is square therefore may not to be limited to 2mm～4mm.That is, the position that residual stress is big also needs finite element is limited to smaller sometimes.In addition, each unit is by three-dimensional coordinate face regulation, the angle of the face of finite element, curvature by with relatively the stipulating of other faces.
Like this, specific for springback occurrence cause portion is carried out, based on angle, curvature member to be carried out the zone and divide, the part that residual stress level is high is as operand, thereby amount of calculation is reduced significantly, can more promptly carry out specific thus springback occurrence cause portion.
Example as calculation process illustrates σ xComponent value coefficient method doubly.
When the components of stress before the computing on the point in selected zone are made as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are made as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0) time, the components of stress (σ after the calculation process x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx) as follows.
σ x＝K×σ x0，σ y＝σ y0，σ z＝σ z0，τ xy＝τ xy0，τ yz＝τ yz0，τ zx＝T zx0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, generally be following situation: carry out the computing that K changes in-1000～1000 scope, thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient use the value by 18 inputs of molding condition input part.Though for illustration, only and σ xK is shown relatively, but, K is similarly changed about other components.
In addition, can be operand also with material behavior (thickness of slab, coefficient of elasticity, plastic coefficient).For example, in physics value physical quantity arithmetic processing section 14, select the part zone of the machined object that obtained by drawing analysis unit 11, coefficient is doubly in calculation process with for example Young's modulus in selected zone.Under this occasion, the thickness of slab t before the computing that is made as selected zone 0, Young's modulus E 0, Poisson's ratio v 0, plastic coefficient F 0, a 0, n 0(σ=F 0(ε+a 0) N0The time) time, thickness of slab t, the Young's modulus E after the computing, Poisson's ratio v, plastic coefficient F, a, n (σ=F (ε+a) nThe time) as follows.
t＝t 0，E＝K×E 0，v＝v 0，F＝F 0，a＝a 0，n＝n 0
At this, K can change in-1000～1000 scope.Though for illustration, only K is shown relatively,, K is similarly changed about other material behavior with E.
In addition,, also the file of having exported the state variable that is obtained by drawing analysis unit 11 can be changed by direct editor, select the zone to carry out calculation process thus in physics value physical quantity arithmetic processing section 14.For example file content can be adopted the software of Wordpad etc. to show, will carry out the component value in the zone of calculation process directly to rewrite by the operation of molding condition input part 18, or replace component value by cut and paste with text editing function.
Thus, can estimate quantitatively selection the zone computing variable (component) to resilience give what the influence.
For example occasion that the Displacements Distribution of finite element node is shown with profile diagram (contour map) in resilience analysis unit 12 and the occasion that shows the section of deformed shape, can be with the computing in the zone selected result's the figure of each variable (component) on springback capacity output screen 19, compare side by side or export to and do not have illustrated printer to compare.In addition, the coordinate figure of the node of assigned address (X, Y, Z) is arbitrarily exported the occasion of confirming on software or by file, the computing in the zone that can relatively select the numerical value of each variable (component), or carry out pictorialization with table software for calculation etc. and come comparison.
The profile diagram that on springback capacity output screen 19, shows, can repeat the physics value physical quantity calculation process of physics value physical quantity arithmetic processing section 14 and the springback capacity of resilience analysis unit 12 on one side and calculate, to the position, the variable that become resilience reason closely carry out specific on one side.
And, by the specific region of selecting to change physics value physical quantity arithmetic processing section 14, the component value in the zone that will select maybe, σ x not only, about branch values such as σ y, τ xy also as calculation process carry out coefficient doubly after, implement the resilience of resilience analysis unit 12 and resolve, make springback capacity output screen 18 show position, the variable that becomes the resilience reason, compare the variable quantity of the displacement that causes by resilience, can carry out specific closely to position, the variable that becomes the resilience reason thus.
In addition,, among the part of machined object, select a plurality of zones, calculation process is side by side implemented in the zone of each selection, can be reduced to reason specific required parsing man-hour thus in physics value physical quantity arithmetic processing section 14.
Figure 10 is the stereogram of the regioselective example of the hat section shaped member that shows that an embodiment of the invention relate to.
For example, at the machined object that is obtained by drawing analysis unit 11 is the occasion of hat section shaped member shown in Figure 5, in physics value physical quantity arithmetic processing section 14, side by side select stress after the shaping shown in Figure 10 greater than a plurality of regional A1, the A2 of setting, regional A1, the A2 that selects implemented σ x is carried out coefficient calculation process doubly, in resilience analysis unit 12, use this calculation process result's resilience to resolve.
The variable quantity of the displacement that is caused by resilience according to the investigation of this resilience analysis result, even each one-sided calculation process of not carrying out, how many influences only adopt 1 time calculation process also can estimate the components of stress value σ x that puts on regional A1, A2 has to resilience.
And, after resilience analysis unit 12 is carried out the resilience parsing, zone when dividing springback capacity again and changing most, judge that the size up to the zone of dividing again is below the setting, about each zone of dividing again, carry out calculation process in physics value physical quantity arithmetic processing section 14, carry out resilience by resilience analysis unit 12 once again and resolve, can carry out the specific of springback occurrence cause to the zone that springback capacity changes most thus, it is specific to carry out high efficiency reason.Moreover, physics value physical quantity arithmetic processing section 14, can carry out physics value and physical quantity changed as performance variable and make that the resilience variable quantity is that convergence below the permissible value is calculated, simultaneously, also can carry out the zoning of specific region is also calculated as the convergence that performance variable changes.
In addition, in drawing analysis unit 11, for example use FInite Element to form parsing, the occasion of the stress after numerically obtaining being shaped, the distribution of strain value, in physics value physical quantity arithmetic processing section 14, when selecting the zone, can select the zone of more than one unit as operand.The FInite Element that is used for the analytic method of a problem that solves in conduct, it is its feature that non-individual body is divided into that limited unit handles, the unit connects at limited the node that is positioned on its limit, by the displacement field in the displacement uniqueness ground determining unit of shape function that defines in the unit one by one and node.For example, show the distribution of the components of stress σ x that is obtained by drawing analysis unit 11, unit that can σ x value is maximum is as selecting the zone, or the group of the unit of σ x maximum and the unit that is connected with this unit is regional as selection.
In addition, in drawing analysis unit 11, for example adopt elementary analytic method to form parsing, the occasion of the stress after numerically obtaining being shaped, the distribution of strain value, in physics value physical quantity arithmetic processing section 14, when selecting the zone, can select the zone of more than one unit of account section as operand.
In addition, in drawing analysis unit 11, for example adopting, FInite Element forms parsing, the occasion of the stress after numerically obtaining being shaped, the distribution of strain value, in physics value physical quantity arithmetic processing section 14, when selecting the zone, can select the zone of more than one point as operand.For example, in the finite element of drawing analysis unit 11 is resolved, the machined object employing is resolved in the shell unit (shell elements) that the thickness of slab direction has 5 points (the slave plate back side is followed successively by 1,2,3,4,5 to the surface), in physics value physical quantity arithmetic processing section 14, can with the unit selected with the backboard face nearest 1 and with the zone of nearest 5 the point in plate surface as operand.Thus, can carry out the evaluation that will separate in the diastrophic influence that is subjected to during distortion with the part contact with curvature of mould.
Figure 11 represents the flow chart of the springback occurrence cause particular procedure that one embodiment of the invention relate to.
In step S101, molding condition input forming condition entry portion 18.Then, enter step S102.
In step S102,11 pairs of drawing analysis units are carried out numeric value analysis by the machined objects of molding condition regulation and are handled, and calculate the stress after the shaping of the goods that carry out drawing, the distribution of strain value.Then enter step S103.
In step S103, determine to become the privileged site in transforming object zone by physics value physical quantity arithmetic processing section 14.At this privileged site, curvature and/or angle are decided in calculation for the drawing product, based on the zone of above-mentioned curvature and/or angular divisions formed products, judge based on above-mentioned one-tenth graphic data about each zone of dividing, determine to carry out the zone of calculation process described later.Then enter step S104.
In step S104, carry out above-mentioned calculation process, carry out conversion in S103, being judged to be physics value in a part of zone that will carry out calculation process and at least one of physical quantity.Then enter step S105.
In step S105, the conversion process of the processed body during based on drawing physics value and physical quantity, calculate springback capacity.Then enter step S106, in step S107, wait screen display springback capacity result simultaneously by the profile demonstration.
In step S106, whether the permissible value of judging springback capacity is in permissible value.If at permissible value with interior then finish this processing.If beyond permissible value, then enter step S108.Moreover, at step S106 and step S108, while the people can observe the appointment that the resilience variable quantity manually carries out privileged site.
In step S108, the physics value physical quantity of change privileged site and processed body.At this, can only change privileged site, also can only carry out the change of physics value physical quantity, the change of the two of can also side by side carrying out them.Then, turn back to step S105.The processing of step S105, step S106, step S108 can be that convergence below the permissible value is calculated and carried out by the resilience variable quantity, repeats to become below the threshold value.
Moreover, also can limit the number of repetition of this processing, end process when surpassing the number of repetition of restriction.Under this occasion, can be stored in the output data file " SB rem.casen.k " of the resilience analysis unit 12 of file storage part 16 by people's investigation, on one side the people observes springback capacity output screen 19, Yi Bian explore, obtain best privileged site.
Figure 12 represents the flow chart of the resilience cause portion particular procedure that one embodiment of the invention relate to.
In step S201, the 1st molding condition input forming condition entry portion 18.Then enter step S202.
In step S202,11 pairs of drawing analysis units are carried out numeric value analysis by the drawing product of the 1st molding condition regulation and are handled, and the stress after the shaping of goods of drawing, the distribution of strain value are carried out in calculating.And, drawing analysis unit 11, to carrying out numeric value analysis, calculate the stress after the shaping of the goods that carry out drawing, the distribution (physics value and physical quantity) of strain value by the machined object that at least one of above-mentioned the 1st molding condition is changed the 2nd molding condition regulation that forms.Then enter step S203.
In step S203, determine to become the privileged site in transforming object zone by physics value physical quantity arithmetic processing section 14.This privileged site, each regional physics value of the drawing product that calculating is calculated with the 1st molding condition and physical quantity, with each the regional physics value of the drawing product of calculating and the difference of physical quantity with the 2nd molding condition, judge that whether this difference value is greater than setting.Have the occasion greater than the difference value of setting, the zone that will have this difference value is defined as privileged site, enters step S204.
In step S204, the calculation process of the physics value in the zone of carrying out in S203, determining and at least one conversion of physical quantity.Then enter step S205.
In step S205, based on calculation process physics value and physical quantity, calculate springback capacity.Then enter step S206, in step S207, wait the screen display springback capacity simultaneously by the profile demonstration.
In step S206, whether the permissible value of judging resilience is in permissible value.If in permissible value, then finish this processing.If beyond permissible value, then enter step S208.
In step S208, the occasion of the condition different in 2nd molding condition of change with the 1st molding condition at step S202, be back to step S202, be used for determining the physics value physical quantity of privileged site, the occasion of the setting that is used to judge in change at step S203, be back to step S203, occasion in the calculation process that changes step S204 is back to step S204.Moreover,, the result of people's view screen output on one side etc. carries out step S208 on one side judging.The processing of step S202～S208 can be that convergence below the permissible value is calculated by the resilience variable quantity, repeats to become below the resilience permissible value.
Moreover, limit the number of repetition of this processing, surpassed the occasion of the number of repetition of restriction, also can end process.Under this occasion, the people can discuss the output data file " SB rem.casen.k " of the resilience analysis unit 12 that is stored in file storage part 16, Yi Bian observe springback capacity output screen 19, Yi Bian explore, obtain best privileged site.
Figure 13 represents the flow chart of the resilience cause portion particular procedure that an embodiment of the invention relate to.
In step S301, molding condition input forming condition entry portion 18.Then enter step S302.
In step S302,11 pairs of drawing analysis units are carried out numeric value analysis by the machined objects of molding condition regulation and are handled, and calculate physics value and physical quantity after the shaping of the goods that carry out drawing.Then enter step S303.
In step S303, resilience analysis unit 12 dissection process are calculated the 1st springback capacity through the machined object of drawing.Then enter step S304.
In step S304, judge whether to have physics value after the resilience of drawing product and at least one zone bigger of physical quantity than setting, having this regional occasion, this zone is defined as the specific region.Perhaps, judge whether to have the zone of the difference of physics value before and after the resilience and physical quantity,, this zone is defined as the specific region having this regional occasion greater than setting.Then enter step S305.
In step S305, the physics value of definite privileged site and at least one of physical quantity data are carried out calculation process.Then enter step S306.
In step S306, based on calculation process the amount and physical quantity data computation the 2nd springback capacity of privileged site.Then enter step S307, in step S308, wait simultaneously the result of screen display the 1st springback capacity or the 2nd springback capacity by the profile demonstration.
In step S307, judge that the 2nd springback capacity is whether in permissible value.If in permissible value, then finish this processing.If beyond permissible value, then enter step S309.
In step S309, the physics value that is used for determining privileged site and physical quantity variable, exclusive disjunction are handled and changed.Determine the occasion of method at the change privileged site, change becomes the physics value of determining object and the variable of physical quantity, and change afore mentioned rules value turns back to step S304.The occasion of the variable that carries out in the change calculation process etc. turns back to step S305.Moreover, on one side judging, the result of people's view screen output on one side etc. carries out step S309.The processing of step S304～S309 can repeat to become below the resilience permissible value by convergence calculating.
Moreover, fastening, can not get the occasion of convergence solution with the pass of resilience permissible value, can limit this processing number of repetition, with the number of repetition end process of restriction.Even under this occasion, the people also can discuss the output data file " SBUVC.casen.k " of the resilience analysis unit 12 that is stored in file storage part 16, Yi Bian observe springback capacity output screen 19, Yi Bian explore, obtain best privileged site.
Figure 14 represents the function constitution map of the springback occurrence cause portion specific device that another embodiment of the invention relates to.
Springback occurrence cause portion specific device 10B shown in Figure 14 compares with springback occurrence cause portion specific device 10A shown in Figure 3, replaces physics value physical quantity arithmetic processing section 21, has the part residual stress and removes handling part 21.
The part residual stress is removed handling part 21, is for removing the handling part of the processing of residual stress as the privileged site that is considered to the position of springback occurrence cause.
This privileged site, with above-mentioned same, curvature and/or angle are decided in calculation for the drawing product, based on the zone of above-mentioned curvature and/or angular divisions formed products, judge based on above-mentioned one-tenth graphic data for each zone of dividing, thereby can determine to carry out the zone of calculation process described later.In addition, select the zone, can observe the curvature and/or the angle that are shown in springback capacity output screen 19 on one side,, also can carry out with specified scope such as coordinate figure on one side based on becoming graphic data (stress distribution etc.) to be undertaken by molding condition input part 18 (for example mouse and keyboard etc.).And this privileged site also can be set a plurality of.
In addition, privileged site, also can perhaps, obtain based on the residual stress distribution of data file " P org.k " by calculating based on springback capacity and/or the residual stress distribution of the data file that obtains by resilience analysis unit 12 " SB org.k " and " SB UVC.casen.k ".
That is, can be with the residual stress concentrations position at the lower dead center place of drawing as privileged site, the residual stress concentrations position in the time of also can be with resilience is as privileged site, can also be with the position of springback capacity maximum as privileged site.
And, the part residual stress is removed handling part 21, can not make among the one-tenth graphic data of formed products, physics value changes and just differentiates the resilience variable quantity by removing residual stress, therefore can obtain the state of the molded article of the intensity of keeping molded article and inhibition resilience.
The judgement of these privileged sites in the occasion that does not satisfy the resilience permissible value, also can be carried out various combinations, generates example automatically and calculates springback capacity.Usually, the springback capacity the best part is that the situation of molded article peripheral part is more, therefore preferably with the residual stress concentrations position preferentially as privileged site.For example, can stipulate to make the position that residual stress is the most concentrated to select candidate as the override of privileged site, secondly, with residual stress be not that maximum but residual stress are big to a certain degree, the springback capacity the best part is as the 2nd preferential candidate of selecting of privileged site, with residual stress with 3rd preferential the select candidate etc. of the 2nd concentrated position as privileged site.
The part residual stress is removed handling part 21, will determine privileged site as described above from output file " SBorg.k " the conduct input data of resilience analysis unit 12, and the residual stress that makes this privileged site is zero or lower value.Then, the data file that has changed this privileged site residual stress is defined as " P trim case1.k ", with the corresponding privileged site in the data replacement " P org.k " of " P trim case1.k ", the file after the displacement is made as " P UVCcase1.k ".
The part residual stress is removed handling part 21 and " P UVCcase1.k " exported to resilience analysis unit 12 is made resilience analysis unit 12 calculate springback capacities.At this moment, resilience analysis unit 12 will be stored in file storage part 16 as the data file " SB UVC.Case1.k " of result of calculation, and the springback capacity shown in this data file is shown in resilience output screen 19.
In addition, the part residual stress is removed handling part 21, residual stress in " SBorg.k " that be stored in file storage part 16 is removed the springback capacity that springback capacity before handling and residual stress of the 1st time in " SB UVC.Case1.k " remove after the processing compare, judge whether that springback capacity reduces to become below the threshold value.
At springback capacity in this occasion below permissible value, end process.At springback capacity in certain occasion more than the permissible value, for the processing of the residual stress that further reduces the residual stress position, the part residual stress remove handling part 21 also can other privileged site and/or, carry out the part residual stress with other residual-stress value and handle.Other privileged site, as mentioned above, residual stress concentrations position on the drawing lower dead center, residual stress residual stress concentrations position more than setting and when the position of springback capacity maximum and resilience etc. becomes candidate, can automatically calculate springback capacity for the situation of single position in these positions or various combinations.At this moment, the position coordinate data of the privileged site of generation and the file of stress data are each generations " P trimcasen.k " for n example.
The part residual stress is removed handling part 21, can serve as input data and " SB UVC.Casen.k " that resiliometer is let it pass is stored in file storage part 16 with " the P UVC.casen.k " that generates accordingly with " P trim casen.k " that each was generated with a plurality of examples.
Figure 15 represents the flow chart of the resilience cause portion particular procedure that an embodiment of the invention relate to.
In step S401, molding condition input forming condition entry portion 18.Then enter step S402.
In step S402,11 pairs of drawing analysis units are carried out numeric value analysis by the machined objects of molding condition regulation and are handled, and calculate the stress after the shaping of the goods that carry out drawing, the distribution of strain value.Then enter step S403.
In step S403, carry out dissection process through the springback capacity of the machined object of drawing by resilience analysis unit 12.Then enter step S404.
In step S404, determine the privileged site that residual stress is removed.This privileged site, curvature and/or angle are decided in calculation for the drawing product, based on the zone of above-mentioned curvature and/or angular divisions formed products, judge based on above-mentioned one-tenth graphic data for each zone of dividing, thereby determine to carry out the zone that residual stress described later is removed (S405).When this determines to be based on drawing or the residual stress during resilience or the springback capacity of residual stress when setting is above wait to determine.Then enter step S405.
In step S405, the residual stress of privileged site is reduced or be zero processing.Then enter step S406.
In step S406, reduce or shape or proterties data when being zero drawing based on the residual stress that makes privileged site, calculate springback capacity.Then enter step S407, simultaneously the result of screen display springback capacity in step S408.
In step S407, whether the permissible value of judging resilience is in permissible value.If in permissible value, then finish this processing.If beyond permissible value, then enter step S409.Moreover step S407 and step S409 are while the people can observe the appointment that the resilience variable quantity manually carries out privileged site.
In step S409, the change privileged site.The change of privileged site, when being based on drawing or the judgement material of residual stress concentrations portion during resilience or the springback capacity maximum portion of residual stress when setting is above etc., selected more than one privileged site.In this step, privileged site can be selected through manually coming by means of molding condition input part 18.Then, turn back to step S405, above-mentioned processing repeats to reach in the resilience permissible value of step S407.
Moreover, fastening the occasion that can not get convergence solution, can limit the number of repetition of this processing, with the number of repetition end process that has limited with the pass of resilience permissible value.Even this occasion, the people can discuss the output data file " SBUVC.casen.k " of the resilience analysis unit 12 that is stored in file storage part 16, on one side the people observes springback capacity output screen 19, Yi Bian explore, obtain best privileged site.
Figure 16 represents the function constitution map of the resilience countermeasure position specific device that an embodiment of the invention relate to.
Resilience countermeasure position specific device 10C has: drawing analysis unit 11, resilience analysis unit 12, be replaced as graphic data generating unit 15, file storage part 16, molding condition input part 18 and springback capacity output screen 19.
The output data of drawing analysis unit 11 for example export resilience analysis unit 12, graphic data generating unit 15 and the file storage part 16 of being replaced as described later to as file " P org.k ".
Resilience analysis unit 12, the output data file " P rem.casen.k " that uses the output data file " P org.k " of drawing analysis unit 11 and be replaced as graphic data generating unit 15 carries out resilience and resolves as the input data.
The output data of resilience analysis unit 12, as the numeric value analysis result output data file " SB org.k " that obtains according to input data file " P org.k " or according to the numeric value analysis result output data file " SB rem.casen.k " that " P rem.casen.k " described later obtains, export resilience analysis unit 12 to, be replaced as graphic data generating unit 15 and file storage part 16.
Be replaced as graphic data generating unit 15, obtain the output data file " Porg.k " of drawing analysis unit 11, the zoning of the scope of from output data file, determine among the position data of contained object member to take some countermeasures (rib etc.) and the scope that can take some countermeasures.Be replaced as graphic data generating unit 15, further obtain the shape of rib and the data of stress distribution " P trim.casen.k " from the rib of the wide variety of shapes of definition in file storage part 16 and the shape stress data chart T1 (not having diagram) of definition stress distribution, above-mentioned zoning and rib shape thereof and stress distribution are replaced, generate and be replaced as graphic data " P rem.casen.k ".
Be replaced as graphic data generating unit 15, repeat to make springback capacity become in the setting with the processing of different shapes, stress distribution data replacement zoning.Such reprocessing, by for whole zonings, shape and stress distribution data that each data of the total data in the shape stress data chart T1 are different are replaced, and are replaced as graphic data and generate.Therefore, about being replaced as graphic data, be to generate a plurality of graphic data " P rem.casen.k " that are replaced as for each example number n.
In addition, be replaced as the graphic data generating unit, also can be with after providing at least one specified data displacement (reducing the countermeasure of the 1st springback capacity) of stress distribution of the shape in one or more zones of the following springback capacity of setting and this selection area with this displacement, and then change the position or the quantity of above-mentioned selection area, calculate and be replaced as graphic data.By being replaced as graphic data based on this, utilize the resilience analysis unit to repeat springback capacity and calculate, on the countermeasure basis that reduces the 1st springback capacity, also can take to reduce the countermeasure of the 2nd springback capacity, therefore can further reduce springback capacity.
In addition, above-mentioned zoning also changes the size of zoning also according to selected shape, stress distribution data, makes springback capacity in setting, therefore can generate the graphic data that is replaced as of all examples.
And, be replaced as graphic data generating unit 15, springback capacity passed through divided by drawing conditions such as drift speed or blank holder pressure, and obtain standardized value, utilize the poster processing soft to carry out profile and show, also can visually catch the relation of drawing condition and springback capacity thus easily.The standardization of such springback capacity and the profile of numerical value thereof show, can be economical and shorten the specific of cause portion that resilience takes place expeditiously, and economical and shorten investigation time of the manufacturing process of molded article expeditiously.
File storage part 16 is the storage parts that are used to preserve as drawing analysis unit 11, resilience analysis unit 12, the data file " P org.k " that is replaced as the output result of graphic data generating unit 15 described later, " SBorg.k ", " P rem.casen.k ", " SB rem.casen.k ", " P trim.casen.k " etc., these data files etc. are in drawing analysis unit 11, resilience analysis unit 12 and be replaced as between the graphic data generating unit 15 the directly occasion of input and output, and this document storage part 16 may not need.
The profile diagram that shows on the springback capacity output screen 19, can be on one side based on by being replaced as the graphic data that is replaced as that graphic data generating unit 15 obtains, repeat springback capacity by resilience analysis unit 12 and calculate, Yi Bian visually judge the result of this resilience countermeasure of displacement of shape, stress distribution.
Figure 17 represents the flow chart of the springback occurrence cause portion particular procedure that one embodiment of the invention relate to.
In step S501, molding condition input forming condition entry portion 18.Then enter step S502.
In step S502,11 pairs of drawing analysis units are carried out numeric value analysis by the machined objects of molding condition regulation and are handled, and calculate the stress after the shaping of the goods that carry out drawing, the distribution of strain value.Then enter step S503.
In step S503, be specified to the zone that becomes the displacement subject area in the machined object of graphic data by being replaced as graphic data generating unit 15.This replaces subject area, is the zone with sufficient size of replaceable shape, stress distribution.And, be replaced as graphic data generating unit 15, select shape, the stress distribution data replace from shape stress distribution data drawing list T1, become the displacement subject area of graphic data, will be replaced as graphic data 15 and pass to resilience analysis unit 12 with this shape, stress distribution data replacement.Then enter step S504.
In step S504,, calculate springback capacity for being replaced as graphic data 15.Then enter step S505, in step S507, wait simultaneously the result of screen display springback capacity by the profile demonstration.
In step S505, judge that springback capacity is whether in permissible value.If in permissible value, then finish this processing.If beyond permissible value, then enter step S507.
In step S507, displacement subject area, the shape data of displacement and at least one of stress distribution are changed.At this, only change displacement subject area also can only change shape data or stress distribution, they can also be changed simultaneously the two.Then, turn back to step S503.The processing of step S503～S507 can become convergence below the permissible value by springback capacity and calculate and carry out, and repeats to reach below the permissible value.
Moreover this processing can not limit the permissible value of springback capacity in advance, by the whole example combinations for displacement subject area, the shape data of replacing, stress distribution, carries out the resilience of step S504 and resolves, and calculates springback capacity and is minimum example.
In addition, in the above description, though be to be target variable with the springback capacity, the problem that displacement subject area, shape data, stress distribution are calculated as the convergence of performance variable, but if the quantity of these performance variables is very big, then solution space is very big, therefore is difficult to find optimum solution sometimes or is following the separating of permissible value.Such occasion, also the algorithm that can use genetic algorithm, neutral net etc. to know is obtained optimum solution etc.
Figure 18 represents the hardware structure diagram of the springback occurrence cause portion specific device that an embodiment of the invention relate to.Above-mentioned drawing analysis unit 11, resilience analysis unit 12, physics value physical quantity arithmetic processing section 14, be replaced as graphic data generating unit 15, part residual stress remove handling part 21, in each handle, can in program 100, stipulate, make computer 90 carry out this program 100.Computer 90 possesses: external memory 96, NIC (NIC) 97, the printer 98 of the input unit 94 of the CPU91 of execution necessary processing, stores processor result's memory 92 (for example RAM (Random Access Memory)), display 93, for example keyboard, mouse and so on, hard disk 95, CD/DVD driver and so on.Computer 90 can be connected with other computer 90A by means of network 99, and network 99 is made of the Ethernet that is connected with NIC97 (registration mark) cable.
Program 100 is stored on the recording medium, packs into from external memory 96, or downloads by network 99 from other computer 90A, by the control of CPU91, is stored on the hard disk 95 of computer 90.Then, the program 100 of preservation is carried out by CPU91, as implementation or threads store in memory 92.For example, in drawing analysis unit 11, resilience analysis unit 12, physics value physical quantity arithmetic processing section 14, be replaced as that graphic data generating unit 15 or part residual stress remove in the handling part 21 each handle, become implementation or thread respectively, above-mentioned data file or data are transfused to output in each implementation or cross-thread.In addition, described each implementation or thread can be present among the other computer 90A dispersedly, everywhere reason computer 90 and other computer 90A dispersion treatment.
In addition, molding condition input part 18, springback capacity output screen 19 shown in Figure 3 also can be respectively input unit 94, display 93.Above-mentioned file storage part 16 also can be a hard disk 95.Program 100 also can be stored in hard disk 95.And, export the above-mentioned profile diagram of display 93 to, may be output to printer 98.
[embodiment 1]
Specify the present invention with embodiment.
Figure 19 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) be that its stereogram, (b) are its sectional view, is made of web face W0, sidewall W1, W2, flange F 1, F2.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Molding condition has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level as the proterties of metallic plate.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.Punch die is identical with the gap and the thickness of slab of drift, is 1.6mm.Coefficient of friction input 0.15 is set 60 tons as blank holder load.
Figure 20 is the sectional position of the hat section shaped member that relates to of expression embodiments of the invention and the figure of windup-degree, (a) be that its stereogram, (b) are its sectional view.
The result that obtains is obtained the coordinate figure of the point in the cross section that constitutes position A shown in Figure 20, B place with the poster processing soft, obtain the angle θ of the web face formation in cross section.So-called springback capacity is angle θ at this.
Figure 21 is the stereogram of the specific region A3 that relates to of expression embodiments of the invention.
Figure 22 is the specific region A4 that relates to of expression embodiments of the invention, the stereogram of A5.
Figure 23 is the stereogram of specific region A6～A9 of relating to of expression embodiments of the invention.
Figure 24 (a) is the stereogram of specific region A10～A14 of relating to of expression embodiments of the invention, and Figure 24 (b) is the sectional view of specific region A10～A14 of relating to of expression embodiments of the invention.
Figure 25 is the stereogram of specific region A15～A17 of relating to of expression embodiments of the invention.
Figure 26 is the stereogram of specific region A18～A21 of relating to of expression embodiments of the invention.
Figure 33 is the stereogram of the specific region A22 that relates to of expression embodiments of the invention.
The following describes for the specific region A3～A21 shown in Figure 21～25, components of stress change is carried out the result of resilience calculation process.
The components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, T Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).
Thickness of slab, former state ground uses the value before the calculation process, coefficient of elasticity, plastic coefficient, former state ground uses the input value of drawing analysis unit 11.Stress, strain value through computing are exported as the operation result file.
Table 1 expression is used to resolve the investigation example through the relation of the variable of calculation process and windup-degree of specific region.Table 1 expression specific region and the variable that in calculation process, has been changed.
Table 1
Figure A200780032577D00561
To belonging to whole unit in each zone, use the intact variable of calculation process that obtains by physics value physical quantity arithmetic processing section 14 to carry out resilience and resolve.
The result that obtains is obtained the coordinate figure of the point in the cross section that constitutes position A shown in Figure 20, B place with the poster processing soft, obtain the angle θ of the web face formation in cross section.
Figure 27 is the figure of the windup-degree of the specific region A3 that relates to of expression embodiments of the invention.
The σ x of regional A3 is carried out the result of calculation process, compare with the resilience analysis result of the occasion of not carrying out calculation process, θ is big.The σ x that distinguishes regional A3 thus is big to the influence that θ gives, distinguishes by this result, and be below the threshold value in order to make the resilience variable quantity, as long as the σ x value of regional A3 is changed.
Figure 28 is the figure of the windup-degree of expression embodiments of the invention the specific region A4, the A5 that relate to.
Regional A4 is carried out the result of calculation process, compare with the result who regional A5 is carried out calculation process, θ is big, even compare with the resilience analysis result of the occasion of not carrying out calculation process, θ is also big.Distinguish that thus the σ x of regional A4 compares with the σ x of regional A5, the influence that θ is given is bigger, distinguish by this result, and be below the threshold value in order to make the resilience variable quantity, as long as the σ x value from regional A5 to regional A4 is changed.
Figure 29 is the figure of the windup-degree of specific region A6～A9 of relating to of expression embodiments of the invention.
Regional A6, A7, A8, A9 are carried out the result of calculation process, compare with the resilience analysis result of the occasion of not carrying out calculation process, θ is little and approach 0.The stress σ x that distinguishes regional A6, A7, A8, A9 thus is bigger to the influence that θ gives.
Figure 30 is the figure of the windup-degree of specific region A10～A14 of relating to of expression embodiments of the invention.Distinguish thus, regional A10～A14 is being carried out in the occasion of regional A11 being carried out calculation process, compare the θ minimum among the result of calculation process with the resilience analysis result of the occasion of not carrying out calculation process.The stress of distinguishing regional A11 thus is bigger to the influence that θ gives.
Figure 31 is the figure of the windup-degree of specific region A15～A17 of relating to of expression embodiments of the invention.In view of this, regional A15～A17 is being carried out in the occasion of regional A16 being carried out calculation process, compare the θ minimum among the result of calculation process with the resilience analysis result of the occasion of not carrying out calculation process.The stress of distinguishing regional A16 thus is bigger to the influence that θ gives.
Figure 32 is the figure of the windup-degree of specific region A18～A21 of relating to of expression embodiments of the invention.Regional A16 is divided into regional A18, A19, A20, A21, the resilience parsing has been carried out in each zone of dividing.In view of this, regional A18～A21 is being carried out in the occasion of regional A20 being carried out calculation process, compare the θ minimum among the result of calculation process with the resilience analysis result of the occasion of not carrying out calculation process.The stress of distinguishing regional A20 thus is bigger to the influence that θ gives.Like this, at first the big zone of degree of influence in big zone is carried out specificly, under this scope, further carry out zone and divide and calculation process, investigate degree of influence, can carry out specific to the generation essential factor of resilience expeditiously thus.
Figure 34 be the specific region A22 that relates to of expression embodiments of the invention change the figure of windup-degree of occasion of operation method.In view of this, in the occasion of operation method 2, compare the θ minimum with the resilience analysis result of the occasion of not carrying out calculation process and the occasion of operation method 1.Distinguish thus,,, can seek θ and reduce by becoming and the approaching stress distribution of σ x of replacing with operation method 2 for the regional A22 that selects.
Figure 35 is the figure of the specific region A23 that relates to of expression embodiments of the invention, 24 windup-degree.At this, regional A23, A24 for the zone of each point appointment of pressing the thickness of slab direction, at this, discuss by adopting count 5 parsing of integration.The point numbering to outer surface (punch die side), has been attached symbol 1～5 from inner surface (drift side) successively.Among regional A11, with the whole points ( point 1,2,4,5) beyond the point 3 of the median plane of thickness of slab direction as regional A23.In addition, among regional A11, with the point 3 of the median plane of thickness of slab direction as regional A24.According to the result of Figure 35, the resilience analysis result among regional A23 and the regional A24 is roughly equal, and the influence of the point 3 of the point (1,2,4,5) beyond the median plane of thickness of slab direction and the median plane of thickness of slab direction is roughly equal.Like this, by pressing each point appointed area of thickness of slab direction, can distinguish thickness of slab direction (the surperficial back side) thus bending stress influence and the in-plane stress influence of thickness of slab direction median plane investigate degree of influence.
[embodiment 2]
Figure 36 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) be that its stereogram, (b) are its sectional view, and the hat section is the left-right symmetry shape.The width of the drift that uses is 80mm, and forming height is 80mm, and the width of processed material is 100mm.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.In addition, consider the symmetric shape of hat section, shown in Figure 27 (b), resolve with symmetric condition for 1/2nd zone.Utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.Punch die is identical with the gap and the thickness of slab of drift, is 1.6mm.Coefficient of friction input 0.15 is set 10 tons as blank holder load.
The result that obtains is obtained coordinate figure before and after the resilience of the point of the position shown in Figure 36 (b) is resolved with the poster processing soft, obtain the opening △ d in cross section.So-called springback capacity is opening △ d at this.
Table 2 illustrates the investigation example through the relation of the variable of calculation process and opening that is used to resolve the specific region.The variable that table 2 illustrates the specific region and changed in calculation process.
Table 2
The investigation example The specific region The variable that calculation process is intact
9 Whole unit σ x＝5×σ x0，σ y＝5×σ y0，σ z＝5×σ z0， τ xy＝5×τ xy0，τ yz＝5×τ yz0，τ zx＝5×τ zx0
10 The point 5 of whole unit σ x＝5×σ x0，σ y＝5×σ y0，σ z＝5×σ z0， τ xy＝5×τ xy0，τ yz＝5×τ yz0，τ zx＝5×τ zx0
In example 9, at first with a unit of element number minimum as the specific region, change the components of stress by computing, carry out resilience and resolve, repeat same parsing with the little order of element number, resolve for whole unit.
In example 10, at first with the point 5 of a unit of element number minimum as the specific region, change the components of stress by computing, carry out resilience and resolve, repeat same parsing with the little order of element number, resolve for the point 5 of whole unit.
The components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, T Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).
Thickness of slab, former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use the input value of drawing analysis unit 11.Computing stress, strain value export as the operation result file.
To belonging to whole unit in each zone, use the intact variable of calculation process that obtains by physics value physical quantity arithmetic processing section 14, carry out resilience and resolve.
Figure 37 is among the opening △ d of whole unit of obtaining in the investigation example 9 that embodiments of the invention relate to, shows the figure of 5% unit of whole unit number with the big order of △ d.Hence one can see that, and the big unit of △ d influence is concentrated near the drift shoulder R zone.
Figure 38 is among the opening △ d of whole unit of obtaining in the investigation example 10 that embodiments of the invention relate to, shows the figure of 5% unit of whole unit number with the big order of △ d.Hence one can see that, and the bigger unit with point 5 of △ d influence is concentrated near the longitudinal wall part central authorities zone.
[embodiment 3]
Figure 39 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) is its stereogram, (b) is its sectional view, is made of web face W0, sidewall W1, W2, flange F 1, F2.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of using static implicit expression solution and is resolved.
Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.Punch die is identical with the gap and the thickness of slab of drift, is 1.6mm.Coefficient of friction input 0.15 is set 600kN as blank holder load.
Figure 40 is the sectional position of the hat section shaped member that relates to of expression embodiments of the invention and the figure of windup-degree, (a) be that its stereogram, (b) are its sectional view.
The result that obtains is obtained the coordinate figure of the point in the cross section that is formed in the position A shown in Figure 40 (a), B place with the poster processing soft, obtain the angle θ that the web face in cross section constitutes.So-called springback capacity is angle θ at this.
In the step of determining privileged site, select to divide the zone of implementing calculation process based on curvature by the finite element data of machined object.Use curvature to calculate and decide program, read the unit information node coordinate value of machined object from the drawing analysis result, the curvature that machined object is decided in calculation, making curvature is that unit more than 0.01 is non-demonstration, thus each zone of flange web face longitudinal wall is divided.In addition, stepped longitudinal wall is divided into 3 regional A35～A37 as shown in figure 41.Each regional stress level after the observation drawing, therefore the stress level height of stepped longitudinal wall is judged regional A35～A37 shown in Figure 41, select as the calculation process object.
In the calculation process step, the operation processing program that employing has a file input medium is obtained output to be had and resolves the stress that obtains, the file of strain by being shaped, about regional A35～A37 shown in Figure 41, utilize the computing means to implement following computing to the point of unit that belongs to the zone for each of each zone.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use at i) value imported in the step.Stress, strain value through computing utilize the file output means to export as the operation result file.
In the resilience analyzing step, the operation result file that will obtain in the calculation process step inputs to existing plate finishing die quasi-solution and analyses software PAM-STAMP, implements resilience and resolves.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.Repeat to resolve with regional number of partitions in each regional calculation process resilience.Each result of obtaining is obtained the coordinate figure of the point in the cross section that constitutes position A shown in Figure 40, B place with the poster processing soft, obtain the windup-degree θ of web face formation in the cross section at position A, B place as springback capacity.
Figure 42 represents the resilience analysis result is carried out that the cross section is estimated and the result of the windup-degree θ that obtains.In view of this,, and regional A35～A37 is being carried out among the result of calculation process, regional A35 has been carried out the occasion of calculation process, windup-degree θ minimum in the occasion of not carrying out calculation process.The stress of distinguishing regional A35 thus is bigger to the influence that windup-degree θ gives.
[embodiment 4]
Except calculation process, carry out calculating similarly to Example 3.In the calculation process step, obtain output with operation processing program the stress that is obtained by the parsing that is shaped, the file of strain are arranged with file input medium, about regional A35～A37 shown in Figure 41, utilize the following computing of computing means enforcement for the point at each of each zone the most surperficial and back side to the thickness of slab direction of unit that belongs to the zone.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，T zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use at i) value imported in the step.Stress, strain value through computing utilize the file output means to export as the operation result file.
Figure 43 represents the resilience analysis result is carried out that the cross section is estimated and the result of the windup-degree that obtains (point at the most surperficial and back side of thickness of slab direction) θ.In view of this,, and regional A35～A37 is being carried out among the result of calculation process, regional A35 has been carried out the occasion of calculation process, windup-degree θ minimum in the occasion of not carrying out calculation process.The stress of distinguishing regional A35 thus is bigger to the influence that windup-degree θ gives.
[embodiment 5]
Figure 44 is stereogram and the sectional view as the example of the machined object of drawing product that expression one embodiment of the invention relate to.If shown in Figure 44, the machined object that becomes analysis object is the hat section shaped member.
At first, as the drawing dissection process, use is analysed software LS-DYNA based on the existing plate finishing die quasi-solution of FInite Element and is carried out dissection process.As the proterties of metallic plate, the data of the high-strength steel sheet of thickness of slab 1.4mm, tensile strength 590MPa level have been used.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.The die clearance is 0mm.Coefficient of friction input 0.1 is set 1400kN as shaping load.
Carry out the program of physics value physical quantity calculation process, for obtain output have the file of resolving the stress that obtains, strain by drawing as input information the row operation of going forward side by side handle, divide the zone of machined object.
The stereogram of the zoning of the occasion that to be expression machined object shown in Figure 44 that one embodiment of the invention are related to carry out 10 five equilibriums at directions X to Figure 45, divide with certain curvature at Y, Z direction.At this, utilize the computing means to implement following computing to the whole points that belong to each regional unit.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use the value that inputs to software LS-DYNA in the punching press dissection process.Stress, strain value through computing utilize the file output means to export as the operation result file.
As the resilience dissection process, use above-mentioned software LS-DYNA.The output result of above-mentioned physics value physical quantity calculation process performing a programme is inputed to software LS-DYNA enforcement resilience parsing.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.Repeat to resolve with regional number of partitions in each regional calculation process resilience.
Figure 46 is the figure as the springback capacity of Z direction displacement of section A before and after the resilience that relates to of expression one embodiment of the invention and focus B.Figure 47 is the profile diagram of each regional springback capacity that carries out calculation process of relating to of expression one embodiment of the invention.This profile diagram is to utilize the output data of software LS-DYNA to be obtained by the poster processing soft LS-PREPOST.By Figure 47 can specificly be that the big zone C of springback capacity is a springback occurrence cause portion.
[embodiment 6]
The object parts are the hat section shaped member of shape and size shown in Figure 44.
At first, in resolving, drawing adopted above-mentioned software LS-DYNA.As the proterties of metallic plate, the data of the high-strength steel sheet of thickness of slab 1.4mm, tensile strength 590MPa level have been used.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.The die clearance is 0mm.Coefficient of friction input 0.1 is set 1400kN as shaping load.
Carry out the program of physics value physical quantity calculation process, the file of resolving the stress that obtains, strain by drawing is arranged as input information in order to obtain output, the row operation of going forward side by side is handled, and divides the zone of machined object.Be divided into each zone shown in Figure 45, utilize the computing means to implement following computing the whole points that belong to each regional unit.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Value before the thickness of slab former state ground use calculation process, coefficient of elasticity, plastic coefficient former state ground use as drawing and resolve the value of being imported by software LS-DYNA.Stress, strain value through computing utilize the file output means to export as the operation result file.
In the resilience dissection process, will input to software LS-DYNA once more by the operation result file that physics value physical quantity calculation process performing a programme obtains, implement resilience and resolve.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
At this, resolve for each resilience that repeats to utilize the calculation process of physics value physical quantity operation processing program and utilize software LS-DYNA in each zone.
Figure 48 is the stereogram that each region contour that one embodiment of the invention relate to is shown springback influence degree.
With the Z direction displacement of the focus B that belongs to section A before and after the resilience shown in Figure 46 as springback capacity.At this, as the springback influence degree of standardized value, be defined as each zoning area * 1000+10 of springback influence degree=springback capacity ÷, each the regional springback influence degree that carries out calculation process shows with the poster processing soft becomes profile diagram shown in Figure 48.By Figure 48 can specificly be that the big region D of springback capacity is a springback occurrence cause portion.
Moreover above-mentioned springback influence degree is to have carried out standardized example according to the zoning, but also can be representative length divided by the drawing product, represents width, represents height, represents any of thickness of slab, tensile strength and the value that obtains.Represent size to obtain springback influence degree by adopting, the size that can judge processed body quantitatively makes the differentiation of privileged site and countermeasure thereof easier to the influence that resilience gives.For example, represent size to carry out standardized springback influence degree longitudinally in employing, demonstrate than adopting horizontal representative size to carry out the occasion of the big value of standardized springback influence degree, the shape of rib can be reduced springback capacity more at extending longitudinally, can easily obtain above-mentioned or the like information.
In addition, for springback influence degree, also can profile show the value that obtains divided by drift speed or blank holder pressure.Show by utilizing such variable to carry out profile, not only can be provided as and reduce springback capacity and at the countermeasure of processed shape, and quantitative judgement material for the drawing condition can be provided, can obtain countermeasure based on the drawing condition to privileged site.
In above-mentioned, as the display packing of springback capacity and springback influence degree, the profile demonstration has been described, show but springback capacity etc. can be carried out vector, or show, or the display packing that employing is known shows with the three-dimensional body of rod etc.
As carrying out the post processor that profile shows, can use LS-PREPOST, HyperView, Animator, Meta-Post etc.
[embodiment 7]
Figure 49 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) is its stereogram, (b) is its sectional view, is made of web face W0, sidewall W1, W2, flange F 1, F2.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.Punch die is identical with the gap and the thickness of slab of drift, is 1.6mm.Coefficient of friction input 0.15 is set 600kN as blank holder load.
Figure 50 is the sectional position of the hat section shaped member that relates to of expression embodiments of the invention and the figure of windup-degree, (a) be that its stereogram, (b) are its sectional view.
The result that obtains is obtained the coordinate figure of the point in the cross section that constitutes the position A shown in Figure 50 (a), B place with the poster processing soft, obtain the angle θ of the web face formation in cross section.So-called springback capacity is windup-degree θ at this.
Figure 51 is the figure of stress distribution after the resilience of expression machined object.Physics value physical quantity arithmetic processing section 14 is selected the zone of the high position 5% of stress level based on the stress distribution that is obtained by resilience analysis unit 12 shown in Figure 51.
Figure 52 is the figure of expression based on the zone of the selection of the stress distribution after the resilience of machined object.The region R of selecting with the high position 5% of stress level shown in Figure 52 10.
Physics value physical quantity arithmetic processing section 14 is to belonging to the following computing of whole points execution of the unit in the zone of selection in this wise.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use the value that inputs to software PAM-STAMP in the punching press dissection process.The stress that calculates, strain value utilize the file output means to export as the operation result file.
As the resilience dissection process, used above-mentioned software PAM-STAMP.The output result of above-mentioned physics value physical quantity calculation process performing a programme is inputed to software PAM-STAMP, implement resilience and resolve.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Figure 53 is the springback capacity when not carrying out calculation process and the figure of the springback capacity when having carried out calculation process.It is clear and definite thus that springback capacity reduces by the calculation process to the region R 10 selected in Figure 52, can be specific for the selection region R 10 shown in Figure 52 are springback occurrence cause portions.
[embodiment 8]
Figure 54 is the sectional position of the hat section shaped member that relates to of expression embodiments of the invention and the figure of flange spring amount, (a) be that its stereogram, (b) are its sectional view.Springback capacity is the displacement of the end points of the section A before and after the resilience, i.e. the spring amount of flange.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.Punch die is identical with the gap and the thickness of slab of drift, is 1.6mm.Coefficient of friction input 0.15 is set 600kN as blank holder load.
Stress distribution before the resilience of (a) expression machined object of Figure 55, (b) stress distribution after the resilience of expression machined object.(a) of Figure 56 is before the resilience of expression machined object and stress difference, (b) after the resilience are the figure of expression based on the zone of stress difference component selections.Figure 56 illustrates the region R of selecting with the high position 5% of stress level 11.
Physics value physical quantity arithmetic processing section 14 based on the stress distribution that is obtained by resilience analysis unit 12 shown in Figure 55, is selected the region R 11 of the high position 5% of the stress level shown in Figure 56.
Physics value physical quantity arithmetic processing section 14 is to belonging to the following computing of whole points execution of the unit in the zone of selection in this wise.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use the value that inputs to software PAM-STAMP in the punching press dissection process.The stress that calculates, strain value utilize the file output means to export as the operation result file.
As the resilience dissection process, used above-mentioned software PAM-STAMP.The output result of above-mentioned physics value physical quantity calculation process performing a programme is inputed to software PAM-STAMP, implement resilience and resolve.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.Springback capacity be the end points of the section A before and after the resilience displacement, be the spring amount of flange.
Figure 57 be comparison one embodiment of the invention relate to do not carry out calculation process the time springback capacity and the figure of the springback capacity when having carried out calculation process.It is clear and definite thus that springback capacity reduces by the calculation process to the region R 11 selected in Figure 56, can be specific for the selection region R 11 shown in Figure 56 are springback occurrence cause portions.
[embodiment 9]
Figure 58 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) is its stereogram, (b) is its sectional view, is made of web face W0, sidewall W1, W2, flange F 1, F2.
In the numeric value analysis that drawing analysis unit 11 and resilience analysis unit 12 are carried out, used based on the existing plate finishing die quasi-solution of FInite Element and analysed software PAM-STAMP.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.In addition, utilize shell unit to come modelling the shape of mould (punch die, drift, anchor clamps), be assumed to rigid body and resolve.The gap of punch die and drift, identical with thickness of slab, be 1.6mm.Coefficient of friction input 0.15 is set 600kN as blank holder load.
Figure 59 is the sectional position of the hat section shaped member that relates to of expression embodiments of the invention and the figure of windup-degree, (a) be that its stereogram, (b) are its sectional view.
The result that obtains is obtained the coordinate figure of the point in the cross section that constitutes the position A shown in Figure 59, B place with the poster processing soft, obtain the angle θ of the web face formation in cross section.So-called springback capacity is angle θ at this.
Figure 60 is the figure that the quantity of state after the drawing of expression the 1st machined object and the 2nd machined object distributes.(a) the 1st machined object shown in is based on the drawing product that the 1st molding condition obtains, and (b) the 2nd machined object shown in is based on the drawing product that the 2nd different molding condition of at least one unit of the 1st molding condition obtains.In Figure 60, the different molding condition unit as the 2nd machined object for shape is changed, has applied rib.
Figure 61 is the figure of the quantity of state difference after the drawing of expression the 1st machined object and the 2nd machined object.Figure 62 is the regioselective figure of expression based on the quantity of state difference.Shown in Figure 61, utilize physics value physical quantity operation processing program, obtain the stress difference of the 1st machined object and the 2nd machined object with spreading all over Zone Full, shown in Figure 62, that the high position 10% of its stress difference is regional as selecting.
Carry out the physics value physical quantity calculation process performing a programme of physics value physical quantity calculation process, above-mentioned selection zone is divided as the operand zone, for the selection that belongs to the 1st machined object whole points of each regional unit carry out following calculation process.
σ x＝0，σ y＝0，σ z＝0，τ xy＝0，τ yz＝0，τ zx＝0
ε x＝ε x0，ε y＝ε y0，ε z＝ε z0，γ xy＝γ xy0，γ yz＝γ yz0，γ zx＝γ zx0
At this, the components of stress before the computing on the point in selected zone are expressed as (σ X0, σ Y0, σ Z0, τ Xy0, τ Yz0, τ Zx0), the components of strain are expressed as (ε X0, ε Y0, ε Z0, γ Xy0, γ Yz0, γ Zx0), the components of stress after the calculation process are made as (σ x, σ y, σ z, τ Xy, τ Yz, τ Zx), the components of strain are made as (ε x, ε y, ε z, γ Xy, γ Yz, γ Zx).Thickness of slab former state ground uses the value before the calculation process, and coefficient of elasticity, plastic coefficient former state ground use the value that inputs to software PAM-STAMP in the punching press dissection process.The stress that calculates, strain value utilize the file output means to export as the operation result file.
As the resilience dissection process, used above-mentioned software PAM-STAMP.The output result of above-mentioned physics value physical quantity calculation process performing a programme is inputed to software PAM-STAMP, implement resilience and resolve.Resilience is resolved the elasticity of having used static implicit expression solution and is resolved.
With the cross section B shown in Figure 59 with respect to the torsion angle of section A as springback capacity.Figure 63 is the figure of the comparison of springback capacity when representing not carry out calculation process and the springback capacity when having carried out calculation process.It is clear and definite thus that springback capacity reduces by the calculation process to the zone selected among Figure 62, can be specific for the selection zone shown in Figure 62 be springback occurrence cause portion.
[embodiment 10]
Figure 64 is the figure of the hat section shaped member that becomes the resilience analysis object that relates to of expression embodiments of the invention, (a) is its stereogram, (b) is its sectional view, is made of web face, sidewall 1,2, flange 1,2.Molding condition as the proterties of metallic plate, has used the data of the high-strength steel sheet of thickness of slab 1.6mm, tensile strength 780MPa level.The numeric value analysis of drawing analysis unit 11 uses FInite Element LS-DYNA, and the resilience of resilience analysis unit 12 is resolved, and uses FInite Element LS-DYNA.
Figure 65 is the figure of the data drawing list of the expression rib of replacing that relates to of embodiments of the invention and stress distribution.In the present embodiment, adopted rib and stress distribution with 3 shapes.
Figure 66 is the figure that the zone of the processed body that relates to of expression embodiments of the invention is divided.As shown in the figure, meeting the mode of the shape shown in Figure 65, on the zone of the shape that can dispose rib and stress distribution of machined object, the zone is divided into the quantity of configurable shape.In addition, on the zone of being divided, be furnished with Position Number make as be numbered n=1,2 ..., n.
Figure 67 is the figure that the appointed area whole district of the processed body that relates at embodiments of the invention has disposed shape and stress.Figure 68 is the definition sectional position that relates to of embodiments of the invention and the figure of torsional capacity.At this, with this torsional capacity as springback capacity.Moreover in the A-stage of not carrying out the replacement Treatment of rib and stress etc., springback capacity is 5.60mm.
Figure 69 is the appointed area configuration rib of the processed body that relates at embodiments of the invention and the figure in zone.As investigation example 1, for processed body by the position configuration shown in the n=1@1 (initial position numbering) situation of rib, the situation that has disposed stress, the situation that has disposed rib and stress carried out numeric value analysis.The shape of using is the shape No.A12w3d shown in Figure 65 (degree of depth 3mm), and springback capacity is replaced at rib displacement, stress displacement, rib and stress and to be respectively 5.58mm, 5.64mm, 5.62mm under each situation.
Figure 70 is that the appointed area of the processed body that relates at embodiments of the invention has disposed the figure in rib and zone.As investigation example 2, for processed body by the position configuration shown in the n=2@1 (the 2nd Position Number) situation of rib, the situation that has disposed stress, the situation that has disposed rib and stress carried out numeric value analysis.Springback capacity is replaced at rib displacement, stress displacement, rib and stress and to be respectively 5.57mm, 5.64mm, 5.62mm under each situation.Like this, do not satisfy the situation of permissible value etc. at springback capacity, can increase Position Number in turn yet and discuss.
Figure 71 is the figure (reducing the countermeasure of the 1st springback capacity) that has disposed 1 rib etc. of the springback capacity minimum that embodiments of the invention relate to.As investigation example 3, for processed body by the position configuration shown in the n=best@1 situation of rib, the situation that has disposed stress, the situation that has disposed rib and stress carried out numeric value analysis.Springback capacity is respectively 5.22mm, 4.26mm, 4.03mm in each situation of rib displacement, stress displacement, rib and stress displacement.Like this, under the situation of the springback capacity minimum when having carried out replacement Treatment a position, 28% resilience is enhanced.
Figure 72 is on the basis of configuration such as 1 rib of the springback capacity minimum that embodiments of the invention relate to, and disposed the figure of the 2nd rib etc.As investigation example 4, for processed body by the position configuration shown in n=1@2, the n=best@1 situation of rib, the situation that has disposed stress, the situation that has disposed rib and stress carried out numeric value analysis.Springback capacity is respectively 5.21mm, 4.34mm, 4.11mm in each situation of rib displacement, stress displacement, rib and stress displacement.Like this, on the preferred example basis of configuration 1 rib etc., for disposing the occasion that the 2nd rib etc. carries out replacement Treatment, 27% resilience is enhanced (countermeasure that reduces the 2nd springback capacity).
Figure 73 is the figure that has disposed the 2nd rib etc. on the basis of 1 configuration such as rib of the springback capacity minimum that embodiments of the invention relate to.As investigation example 5, for processed body by the position configuration shown in n=2@2, the n=best@1 situation of rib, the situation that has disposed stress, the situation that has disposed rib and stress carried out numeric value analysis.Springback capacity is respectively 5.23mm, 4.32mm, 4.09mm in each situation of rib displacement, stress displacement, rib and stress displacement.Like this, on the preferred example basis that has disposed 1 rib etc., carry out the 2nd rib etc. is disposed at the occasion of the 2nd replacement Treatment on the Position Number, 27% resilience is enhanced.Like this, do not satisfy the occasion of permissible value etc. at springback capacity, can increase Position Number in turn yet and discuss.
Figure 74 is the figure that the appointed area of the processed body that relates at embodiments of the invention has disposed n rib etc.As investigation example 6, be after the position n=best@1 that springback capacity becomes minimum when having disposed 1 rib etc., this replacement Treatment that the 2nd rib etc. is disposed at the position n=best@2 that springback capacity becomes minimum is carried out n time example.Springback capacity is respectively 3.22mm, 1.60mm, 1.51mm in each situation of rib displacement, stress displacement, rib and stress displacement.Like this, the occasion of the replacement Treatment of repetition is carried out in the position of exploring the n suboptimum, and 73% resilience is enhanced.
And, as investigation example 7, carried out n position in Figure 74, finding, be replaced as the investigation of the situation of the shape No.A12w3d different (degree of depth 4mm) with Figure 74.Under this situation, springback capacity is respectively 3.20mm, 1.54mm, 1.46mm in each situation of rib displacement, stress displacement, rib and stress displacement.Like this, explore the position of n suboptimum, carry out under the situation of replacement Treatment of repetition, 74% resilience is enhanced.
Figure 75 is the figure that the appointed area of the processed body that relates at embodiments of the invention has disposed actual rib.As investigation example 8, the rib in the position configuration reality shown in Figure 74 carries out numeric value analysis.Even at the big rib of the position configuration reality that has disposed small rib strip shown in Figure 74 like that, springback capacity also is 1.60mm, and 72% resilience is enhanced.Like this, even, disposed under the situation of actual rib, also can seek equal resilience and cut down based on the investigation of the allocation position of above-mentioned rib etc.
Table 3 illustrates the resilience result about above-mentioned example 1～8.
Table 3
Figure A200780032577D00721
Like this, the present invention, except the resilience occurrence positions is carried out specific, can also calculate its countermeasure.
As above for example shown in, the resilience happening part specific device that the present invention relates to can come the investigation of the resilience of in the past using actual device to carry out to carry out specific to the resilience happening part by numeric value analysis.In addition, carry out quantitative analysis, and show its information, seizure visually, can easily carry out specific springback occurrence cause portion by springback influence degree to machined object.And the present invention can provide the countermeasure at the machined object resilience.
Like this, the springback occurrence cause portion specific device that the present invention relates to, the investigation of the resilience of in the past using actual device to carry out can be come springback occurrence cause portion is carried out specific by numeric value analysis, test man-hour and the expense of the design phase of formed parts are reduced.In addition, though used actual device to carry out the investigation of resilience in the past, therefore but the resilience countermeasure position specific device that the present invention relates to can carry out specificly to resilience countermeasure position by numeric value analysis, test man-hour and the expense of the design phase of formed parts are very reduced.
In addition, such device involved in the present invention is expected to be applied to all machined objects, therefore brings huge interests to industrial circle.
More than Shuo Ming embodiment is only listed as typical case; make up inscape, its distortion and the variation of this each embodiment; concerning those skilled in the art, all can clearly know, clearly carry out the various distortion of above-mentioned embodiment under the situation of the invention scope that can in the scheme that does not break away from principle of the present invention and ask for protection, be put down in writing of those skilled in the art.

Claims (105)

1. springback occurrence cause ad hoc approach is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
Carry out the calculation process step of calculation process for the physics value in parts among the one-tenth graphic data of above-mentioned drawing product, above-mentioned drawing product zone and at least one of physical quantity data; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
2. springback occurrence cause ad hoc approach according to claim 1, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
3. springback occurrence cause ad hoc approach according to claim 1 and 2, have following step: above-mentioned a part of zone repeats above-mentioned calculation process step and springback capacity is calculated step by changing, and comes above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific.
4. according to each described springback occurrence cause ad hoc approach of claim 1～3, has following step:, come above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific by changing at least one of above-mentioned physics value and physical quantity and/or above-mentioned calculation process repeats above-mentioned calculation process step and springback capacity is calculated step.
5. according to each described springback occurrence cause portion ad hoc approach of claim 1～4, above-mentioned part zone is a plurality of zones, side by side carries out above-mentioned calculation process for each zone.
6. according to each described springback occurrence cause portion ad hoc approach of claim 1～5, above-mentioned calculation process step, be to divide the zone that above-mentioned springback capacity becomes minimum, for each zone of this division, at least one of above-mentioned physics value and physical quantity carried out the size of calculation process up to the zone of dividing become step below the setting.
7. the ad hoc approach of the described springback occurrence cause of each of claim 1～6, above-mentioned a part of zone is more than one unit or unit of account section.
8. according to the ad hoc approach of each described springback occurrence cause of claim 1～7, above-mentioned part zone is more than one point.
9. springback occurrence cause specific device is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product;
The above-mentioned one-tenth graphic data of numeric value analysis is calculated the resilience analysis unit of springback capacity; With
The physics value in a part of zones among the one-tenth graphic data of above-mentioned drawing product, above-mentioned drawing product and at least one of physical quantity data are carried out calculation process, make above-mentioned resilience analysis unit calculate arithmetic processing section based on the result's of above-mentioned calculation process springback capacity.
10. springback occurrence cause specific device according to claim 9, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
11. according to claim 9 or 10 described springback occurrence cause specific devices, change above-mentioned a part of zone by above-mentioned arithmetic processing section and repeat the result of above-mentioned calculation process and calculating of springback capacity, come above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific.
12. each described springback occurrence cause specific device according to claim 9～11, change above-mentioned physics value and physical quantity and/or above-mentioned calculation process by above-mentioned arithmetic processing section and repeat the result of above-mentioned calculation process and calculating of springback capacity, come above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific.
13. according to each described springback occurrence cause specific device of claim 9～12, above-mentioned part zone is a plurality of zones, side by side carries out above-mentioned calculation process for each zone.
14. each described springback occurrence cause specific device according to claim 9～13, above-mentioned arithmetic processing section, divide the zone that above-mentioned springback capacity becomes minimum, for each zone of this division, at least one of above-mentioned physics value and physical quantity carried out the size of calculation process up to the zone of dividing become below the setting.
15. a springback occurrence cause specific program is characterized in that, makes computer carry out following steps:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
The physics value in parts among the one-tenth graphic data of above-mentioned drawing product, above-mentioned drawing product zone and at least one of physical quantity data are carried out the calculation process step of calculation process; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
16. springback occurrence cause specific program according to claim 15, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
17. according to claim 15 or 16 described springback occurrence cause specific programs, have following step: above-mentioned a part of zone repeats above-mentioned calculation process step and above-mentioned springback capacity is calculated step by changing, and comes above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific.
18. each described springback occurrence cause specific program according to claim 15～17, has following step:, come above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific by changing at least one of above-mentioned physics value and physical quantity and above-mentioned calculation process repeats above-mentioned calculation process step and above-mentioned springback capacity is calculated step.
19. according to each described springback occurrence cause portion specific program of claim 15～18, above-mentioned part zone is a plurality of zones, side by side carries out above-mentioned calculation process for each zone.
20. each described springback occurrence cause portion specific program according to claim 15～19, above-mentioned calculation process step, it is the zone of above-mentioned division springback capacity minimum, for each zone of this division, at least one of above-mentioned physics value and physical quantity carried out the size of calculation process up to the zone of dividing become step below the setting.
21. springback occurrence cause specific program, it is the molding condition of numeric value analysis drawing, calculate one-tenth graphic data and the numeric value analysis program of springback capacity and the springback occurrence cause portion specific program that can carry out the input and output of this numeric value analysis program and data of drawing, it is characterized in that, make computer carry out following steps:
Obtain the step of above-mentioned one-tenth graphic data by above-mentioned numeric value analysis program;
The physics value in parts among the one-tenth graphic data of above-mentioned drawing product, above-mentioned drawing product zone and at least one of physical quantity are carried out the step of calculation process; With
Export the result of above-mentioned calculation process to above-mentioned numeric value analysis program, and make it calculate the step of above-mentioned springback capacity.
22. springback occurrence cause specific program according to claim 21, have following step: above-mentioned a part of zone repeats above-mentioned calculation process step and above-mentioned springback capacity is calculated step by changing, and comes above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific.
23. according to claim 21 or 22 described springback occurrence cause specific programs, has following step:, come above-mentioned springback capacity is become hour or the difference of springback capacity when carrying out resilience with not carrying out above-mentioned calculation process and resolving above-mentioned zone, physics value and the physical quantity when becoming maximum carried out specific by changing at least one of above-mentioned physics value and physical quantity and above-mentioned calculation process repeats above-mentioned calculation process step and above-mentioned springback capacity is calculated step.
24. a computer-readable recording medium is characterized in that, records each described springback occurrence cause portion specific program of claim 15～23.
25. a springback occurrence cause ad hoc approach is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
Calculate among the one-tenth graphic data of fixed above-mentioned drawing product, for the curvature and/or the angle of above-mentioned drawing product, zone based on above-mentioned curvature and/or angular divisions formed products, for each zone of dividing, to carrying out the calculation process step of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity based on above-mentioned shaping data judging; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
26. springback occurrence cause ad hoc approach according to claim 25, above-mentioned calculation process step, to the point of the big relatively above-mentioned zoning of above-mentioned curvature and/or angle at least one above-mentioned physics value and at least one of physical quantity variable carry out calculation process.
27. according to claim 25 or 26 described springback occurrence cause ad hoc approach, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
28. each described springback occurrence cause ad hoc approach according to claim 25～27, have following step: above-mentioned a part of zone repeats above-mentioned calculation process step and springback capacity is calculated step by changing, and comes to carry out specific to above-mentioned zone, physics value and physical quantity that above-mentioned springback capacity becomes hour.
29. each described springback occurrence cause ad hoc approach according to claim 25～28, has following step:, come to carry out specific to above-mentioned zone, physics value and physical quantity that above-mentioned springback capacity becomes hour by changing at least one of above-mentioned physics value and physical quantity and/or above-mentioned calculation process repeats above-mentioned calculation process step and springback capacity is calculated step.
30. a springback occurrence cause specific device is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product;
Calculate among the one-tenth graphic data of fixed above-mentioned drawing product, for the curvature and/or the angle of above-mentioned drawing product, zone based on above-mentioned curvature and/or angular divisions formed products, for each zone of dividing, to carrying out the arithmetic processing section of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity variable based on above-mentioned shaping data judging; With
Based on the result of above-mentioned calculation process, calculate the resilience analysis unit of springback capacity.
31. springback occurrence cause specific device according to claim 30, physics value and physical quantity in the big relatively above-mentioned zoning of above-mentioned curvature and/or angle are the above-mentioned physics value and the physical quantity of at least one point of the big relatively above-mentioned zoning of above-mentioned curvature and/or angle.
32. according to claim 30 or 31 described springback occurrence cause specific devices, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
33. each described springback occurrence cause specific device according to claim 30～32, change above-mentioned a part of zone by above-mentioned arithmetic processing section and above-mentioned resilience analysis unit and repeat the result of above-mentioned calculation process and calculating of springback capacity, above-mentioned arithmetic processing section is carried out specific to above-mentioned zone, physics value and physical quantity that above-mentioned springback capacity becomes hour.
34. each described springback occurrence cause specific device according to claim 30～33, change above-mentioned physics value and physical quantity and/or above-mentioned calculation process by above-mentioned arithmetic processing section and above-mentioned resilience analysis unit and repeat the result of above-mentioned calculation process and calculating of springback capacity, above-mentioned arithmetic processing section is carried out specific to above-mentioned springback capacity above-mentioned zone, physics value and physical quantity hour.
35. a springback occurrence cause specific program is characterized in that, makes computer carry out following steps:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
Calculate among the one-tenth graphic data of fixed above-mentioned drawing product, for the curvature and/or the angle of above-mentioned drawing product, zone based on above-mentioned curvature and/or angular divisions formed products, for each zone of dividing, to carrying out the calculation process step of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity based on above-mentioned shaping data judging; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
36. springback occurrence cause specific program according to claim 35, above-mentioned calculation process step is carried out calculation process to the above-mentioned physics value of the point of the big relatively above-mentioned zoning of above-mentioned curvature and/or angle and at least one of physical quantity variable.
37. springback occurrence cause specific program, it is the molding condition of numeric value analysis drawing, calculate one-tenth graphic data and the numeric value analysis program of springback capacity and the springback occurrence cause portion specific program that can carry out the input and output of this numeric value analysis program and data of drawing, it is characterized in that, make computer carry out following steps:
Obtain the step of above-mentioned one-tenth graphic data by above-mentioned numeric value analysis program;
Calculate among the one-tenth graphic data of fixed above-mentioned drawing product, for the curvature and/or the angle of above-mentioned drawing product, zone based on above-mentioned curvature and/or angular divisions formed products, for each zone of dividing, to carrying out the calculation process step of calculation process for the physics value in a part of zone that will carry out calculation process and at least one of physical quantity based on above-mentioned shaping data judging; With
Export the result of above-mentioned calculation process to above-mentioned numeric value analysis program, and make it calculate the step of above-mentioned springback capacity.
38. according to the described springback occurrence cause specific program of claim 37, above-mentioned calculation process step is carried out calculation process to the above-mentioned physics value of the point of the big relatively above-mentioned zoning of above-mentioned curvature and/or angle and at least more than one variable of physical quantity variable.
39. a computer-readable recording medium is characterized in that, records the described computer program of claim 35～37.
40. a springback influence degree display method is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
At least one of each regional physics value among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and physical quantity data carried out the calculation process step of calculation process;
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step; With
For above-mentioned each zone each, show the step display of above-mentioned springback capacity of calculating.
41. according to the described springback influence degree display method of claim 40, above-mentioned calculation process step is carried out calculation process to the physics value of a part of point in each zone among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and at least one of physical quantity data.
42. according to claim 40 or 41 described springback influence degree display methods, above-mentioned step display is for the above-mentioned springback capacity of calculating of each profile demonstration in above-mentioned each zone.
43. according to each described springback influence degree display method of claim 40～42, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
44. according to each described springback influence degree display method of claim 40～43, above-mentioned step display for above-mentioned each zone each, shows the value that above-mentioned springback capacity of calculating obtains divided by above-mentioned each regional area.
45. each described springback influence degree display method according to claim 40～43, above-mentioned step display, for above-mentioned each zone each, show above-mentioned springback capacity of calculating divided by the representative length of above-mentioned drawing product, represent width, represent height, represent any in thickness of slab, the tensile strength and the value that obtains.
46. according to each described springback influence degree display method of claim 40～43, above-mentioned step display for above-mentioned each zone each, shows that above-mentioned springback capacity of calculating is divided by drift speed or blank holder pressure and the value that obtains.
47. a springback influence degree display unit is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product;
The above-mentioned one-tenth graphic data of numeric value analysis is calculated the resilience analysis unit of springback capacity;
Each regional physics value among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and at least one of physical quantity data are carried out calculation process, make above-mentioned resilience analysis unit calculate arithmetic processing section based on the result's of above-mentioned calculation process springback capacity; With
For above-mentioned each zone each, show the display part of above-mentioned springback capacity of calculating.
48. according to the described springback influence degree display unit of claim 47, above-mentioned arithmetic processing section is carried out calculation process to the physics value of the point of each regional part among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and at least one of physical quantity data.
49. according to claim 47 or 48 described springback influence degree display unit, above-mentioned display part shows above-mentioned springback capacity of calculating for each profile in above-mentioned each zone.
50. according to each described springback influence degree display unit of claim 47～49, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
51. according to each described springback influence degree display unit of claim 47～50, above-mentioned display part, for above-mentioned each zone each, the value that above-mentioned springback capacity of calculating is obtained divided by above-mentioned each regional area is shown in above-mentioned display part.
52. each described springback influence degree display unit according to claim 47～50, above-mentioned display part, for each of above-mentioned each zone, with above-mentioned springback capacity of calculating divided by the representative length of above-mentioned drawing product, represent width, represent height, represent any in thickness of slab, the tensile strength and the value that obtains is shown in above-mentioned display part.
53. each described springback influence degree display unit according to claim 47～50, above-mentioned display part, for each of above-mentioned each zone, with above-mentioned springback capacity of calculating divided by drift speed or blank holder pressure and the value profile that obtains is shown in above-mentioned display part.
54. a springback influence degree display routine is characterized in that, makes computer carry out following steps:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
At least one of each regional physics value among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and physical quantity data carried out the calculation process step of calculation process;
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step; With
For above-mentioned each zone each, above-mentioned springback capacity of calculating is shown in the step display of display part.
55. according to the described springback influence degree display routine of claim 54, above-mentioned calculation process step is carried out calculation process to the physics value of a part of point in each zone among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and at least one of physical quantity data.
56. according to claim 54 or 55 described springback influence degree display routines, above-mentioned step display, for each of above-mentioned each zone, profile shows above-mentioned springback capacity of calculating.
57. according to each described springback influence degree display routine of claim 54～56, above-mentioned step display is to show the step of the value that above-mentioned springback capacity of calculating obtains divided by above-mentioned each regional area for each of above-mentioned each zone.
58. according to each described springback influence degree display routine of claim 54～56, above-mentioned step display be for each of above-mentioned each zone show above-mentioned springback capacity of calculating divided by the representative length of above-mentioned drawing product, represent width, represent highly, represent any in thickness of slab, the tensile strength and the step of the value that obtains.
59. according to each described springback influence degree display routine of claim 54～56, above-mentioned step display is the step that shows the value that above-mentioned springback capacity of calculating obtains divided by drift speed or blank holder pressure for each of above-mentioned each zone.
60. springback occurrence cause display routine, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product drawing analysis program, numeric value analysis this become graphic data, calculate the resilience analysis program of springback capacity, the post processor that profile shows this springback capacity and the springback influence degree display routine that can carry out the data input and output, it is characterized in that, make computer carry out following steps:
Obtain the step of above-mentioned one-tenth graphic data by above-mentioned drawing analysis program;
At least one of each regional physics value among the above-mentioned one-tenth graphic data, that divided above-mentioned drawing product and physical quantity data carried out the calculation process step of calculation process;
The result of above-mentioned calculation process is outputed to the step of above-mentioned resilience analysis program; With
Make above-mentioned post processor profile show that above-mentioned resilience analysis program carries out aforementioned calculation for each of each zone and the profile step display of the springback capacity that obtains.
61. according to the described springback occurrence cause display routine of claim 60, above-mentioned calculation process step is carried out calculation process to the physics value of a part of point in each zone among the one-tenth graphic data of above-mentioned drawing product, that divided above-mentioned drawing product and at least one of physical quantity data.
62. a computer-readable recording medium is characterized in that, records each described computer program of claim 54～61.
63. a springback occurrence cause portion ad hoc approach is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
Based on the one-tenth graphic data of above-mentioned drawing product, calculate after the resilience a plurality of zones separately physics value and the 1st springback capacity of the data of physical quantity calculate step;
Under the situation at least one big zone, at least one of this regional physics value and physical quantity data carried out the calculation process step of calculation process than setting with above-mentioned physics value and physical quantity; With
Based on the result of above-mentioned calculation process, further calculate after the resilience a plurality of zones separately physics value and the 2nd springback capacity of the data of physical quantity calculate step.
64. according to the described springback occurrence cause portion ad hoc approach of claim 63, above-mentioned calculation process step, under the situation in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the resilience and physical quantity data, at least one of this regional physics value and physical quantity data carried out calculation process.
65. according to claim 63 or 64 described springback occurrence cause portion ad hoc approach, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
66. a springback occurrence cause portion specific device is characterized in that having:
The molding condition of numeric value analysis drawing obtains the drawing analysis unit of the one-tenth graphic data of drawing product;
The above-mentioned one-tenth graphic data of numeric value analysis is calculated the resilience analysis unit of springback capacity;
Under the situation at least one big zone than setting with above-mentioned physics value and physical quantity, at least one of carrying out this regional physics value and physical quantity data carried out the calculation process of calculation process, makes above-mentioned resilience analysis unit further calculate the arithmetic processing section of a plurality of zones springback capacity separately after the resilience based on the result of above-mentioned calculation process.
67. according to the described springback occurrence cause portion specific device of claim 66, above-mentioned arithmetic processing section, under the situation in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the resilience and physical quantity data, at least one of this regional physics value and physical quantity data carried out calculation process.
68. according to claim 66 or 67 described springback occurrence cause portion specific devices, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
69. a springback occurrence cause portion specific program is characterized in that, makes computer carry out following steps:
The molding condition of numeric value analysis drawing obtains the drawing analyzing step of the one-tenth graphic data of drawing product;
Based on the one-tenth graphic data of above-mentioned drawing product, calculate after the resilience a plurality of zones separately physics value and the 1st springback capacity of the data of physical quantity calculate step;
Under the situation at least one big zone, at least one of this regional physics value and physical quantity data carried out the calculation process step of calculation process than setting with above-mentioned physics value and physical quantity; With
Based on the result of above-mentioned calculation process, further calculate after the resilience a plurality of zones separately physics value and the 2nd springback capacity of the data of physical quantity calculate step.
70. according to the described springback occurrence cause portion specific program of claim 69, above-mentioned calculation process step, under the situation in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the resilience and physical quantity data, at least one of this regional physics value and physical quantity data carried out calculation process.
71. springback occurrence cause portion specific program, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product and a plurality of zone physics value separately and the drawing analysis program of the data of physical quantity, this becomes graphic data and a plurality of zone physics value separately and the data of physical quantity numeric value analysis, calculate a plurality of zones physics value separately and the resilience analysis program of the data of physical quantity after springback capacity and the resilience, with the springback occurrence cause portion specific program that can carry out the data input and output, it is characterized in that, make computer carry out following steps:
Obtain a plurality of zones physics value separately and the step of the data of physical quantity after the above-mentioned resilience by above-mentioned resilience analysis program;
Under the situation at least one big zone, at least one of this regional physics value and physical quantity data carried out the calculation process step of calculation process than setting with above-mentioned physics value and physical quantity; With
The result of above-mentioned calculation process is outputed to the step of above-mentioned resilience analysis program.
72. according to the described springback occurrence cause portion specific program of claim 71, also have by above-mentioned drawing analysis program and obtain a plurality of zones physics value separately and the step of the data of physical quantity before the above-mentioned resilience,
Above-mentioned calculation process step, under the situation in the big zone of at least one differential ratio setting of physics value with the drawing product before and after the above-mentioned resilience and physical quantity data, at least one of this regional physics value and physical quantity data carried out calculation process.
73. according to each described springback occurrence cause portion specific program of claim 69～72, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
74. a computer-readable recording medium is characterized in that, records each described computer program of claim 69～73.
75. a springback occurrence cause portion ad hoc approach is characterized in that having:
Numeric value analysis the 1st molding condition is calculated the 1st drawing analyzing step of the 1 one-tenth graphic data;
At least one the 2nd different molding condition of above-mentioned the 1st molding condition of numeric value analysis is calculated the 2nd drawing analyzing step of the 2 one-tenth graphic data;
The a plurality of zones with above-mentioned the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, under the situation in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of above-mentioned the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out the step of calculation process; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
76. according to the described springback occurrence cause portion ad hoc approach of claim 75, above-mentioned molding condition is the shape of steel plate and proterties, mold shape, punching press condition.
77. according to claim 75 or 76 described springback occurrence cause portion ad hoc approach, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
78. a springback occurrence cause portion specific device is characterized in that having:
Numeric value analysis the 1st molding condition is calculated the 1 one-tenth graphic data, and at least one the 2nd different molding condition of numeric value analysis the 1st molding condition is calculated the drawing analysis unit of the 2 one-tenth graphic data;
The above-mentioned one-tenth graphic data of numeric value analysis is calculated the resilience analysis unit of springback capacity; With
The a plurality of zones with above-mentioned the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, under the situation in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of above-mentioned the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out calculation process, make above-mentioned resilience analysis unit calculate arithmetic processing section based on the result's of this calculation process springback capacity.
79. according to the described springback occurrence cause portion specific device of claim 78, above-mentioned molding condition is the shape of steel plate and proterties, mold shape, punching press condition.
80. according to claim 78 or 79 described springback occurrence cause portion specific devices, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
81. a springback occurrence cause portion specific program is characterized in that, makes computer carry out following steps:
Numeric value analysis the 1st molding condition is calculated the 1st drawing analyzing step of the 1 one-tenth graphic data;
At least one the 2nd different molding condition of above-mentioned the 1st molding condition of numeric value analysis is calculated the 2nd drawing analyzing step of the 2 one-tenth graphic data;
The a plurality of zones with above-mentioned the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, under the situation in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of above-mentioned the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out the step of calculation process; With
Based on the result of above-mentioned calculation process, calculate the springback capacity of springback capacity and calculate step.
82. springback occurrence cause portion specific program, it is the molding condition of numeric value analysis drawing, calculate the one-tenth graphic data of drawing product and a plurality of zone physics value separately and the drawing analysis program of the data of physical quantity, this becomes graphic data and a plurality of zone physics value separately and the data of physical quantity numeric value analysis, calculate a plurality of zones physics value separately and the numeric value analysis program of the data of physical quantity after springback capacity and the resilience, with the springback occurrence cause portion specific program that can carry out the data input and output, it is characterized in that, make computer carry out following steps:
Obtain the step of the 1 one-tenth graphic data by above-mentioned drawing analysis program;
By above-mentioned drawing analysis program, at least one the 2nd different molding condition of above-mentioned the 1st molding condition of numeric value analysis is obtained the step of the 2 one-tenth graphic data;
The a plurality of zones with above-mentioned the 1 one-tenth graphic data and the 2 one-tenth graphic data separately physics value and the data of physical quantity among, under the situation in zone that at least one differential ratio setting of physics value and physical quantity data is big, the physics value of above-mentioned the 1 one-tenth graphic data in this zone and at least one of physical quantity data are carried out the step of calculation process; With
The result of above-mentioned calculation process is outputed to the step of above-mentioned resilience analysis program.
83. 2 described springback occurrence cause portion specific programs according to Claim 8, above-mentioned molding condition are the shape of steel plate and proterties, mold shape, punching press condition.
84. the described springback occurrence cause portion specific program of each of 1～83 according to Claim 8, above-mentioned physics value and physical quantity are the component value of thickness of slab, coefficient of elasticity, plastic coefficient, stress, the component value of strain.
85. a computer-readable recording medium is characterized in that, records each described computer program of claim 81～84.
86. a resilience countermeasure position specifying method is characterized in that having:
The molding condition of numeric value analysis drawing product is calculated the drawing analyzing step of the one-tenth graphic data of drawing product;
At least a portion zones among the one-tenth graphic data of selected drawing product, drawing product, calculate with specified data replace that at least one of stress distribution of the shape of this selection area and this selection area form be replaced as graphic data be replaced as graphic data generation step; With
The above-mentioned graphic data that is replaced as of numeric value analysis is calculated the springback capacity of springback capacity and is calculated step.
87. 6 described resilience countermeasure position specifying methods according to Claim 8 change the position or the quantity of above-mentioned selection area, repeat above-mentionedly to be replaced as that graphic data generates step and above-mentioned springback capacity is calculated step.
88. 6 or 87 described resilience countermeasure position specifying methods also have and judge whether above-mentioned springback capacity becomes the following step of setting according to Claim 8.
89. 8 described resilience countermeasure position specifying methods according to Claim 8, the above-mentioned graphic data that is replaced as generates step, at least one specified data displacement with this displacement of the stress distribution of the shape in one or more zones of the springback capacity below the setting and this selection area will be provided, and then change the position or the quantity of above-mentioned selection area, repeat above-mentionedly to be replaced as that graphic data generates step and above-mentioned springback capacity is calculated step.
90. the described resilience countermeasure of each of 6～89 position specifying method also has the step that the one-tenth graphic data of above-mentioned drawing product is divided with the zone of regulation according to Claim 8,
The above-mentioned graphic data that is replaced as generates step, for the Zone Full in the zone of dividing with the zone of afore mentioned rules, calculates the above-mentioned graphic data that is replaced as.
91. the described resilience countermeasure of each of 6～90 position specifying method according to Claim 8, the afore mentioned rules data, for above-mentioned selection zone each, at least one of shape and stress distribution is different.
92. a resilience countermeasure position specific device is characterized in that having:
The molding condition of numeric value analysis drawing product is calculated the drawing analysis unit of the one-tenth graphic data of drawing product;
At least a portion zones among the one-tenth graphic data of selected drawing product, drawing product, generation with the stress distribution of the shape of this selection area and this selection area at least one with specified data replace form be replaced as graphic data be replaced as the graphic data generating unit; With
The above-mentioned graphic data that is replaced as of numeric value analysis is calculated the resilience analysis unit of springback capacity.
93. according to the described resilience countermeasure of claim 92 position specific device, change the position or the quantity of above-mentioned selection area, the above-mentioned graphic data generating unit that is replaced as is calculated the above-mentioned graphic data that is replaced as, and above-mentioned resilience analysis unit repeats springback capacity and calculates.
94. according to claim 92 or 93 described resilience countermeasure position specific devices, the above-mentioned graphic data generating unit that is replaced as judges whether above-mentioned springback capacity becomes below the setting.
95. according to the described resilience countermeasure of claim 94 position specific device, the above-mentioned graphic data generating unit that is replaced as, at least one specified data displacement with this displacement of the stress distribution of the shape in one or more zones of the springback capacity below the setting and this selection area will be provided, and then change the position or the quantity of above-mentioned selection area, the above-mentioned graphic data generating unit that is replaced as is calculated the above-mentioned graphic data that is replaced as, and above-mentioned resilience analysis unit repeats springback capacity and calculates.
96. each described resilience countermeasure position specific device according to claim 92～95, the above-mentioned graphic data generating unit that is replaced as, the one-tenth graphic data of above-mentioned drawing product is divided with the zone of regulation, and the Zone Full for the zone of dividing with the zone of afore mentioned rules generates the above-mentioned graphic data that is replaced as.
97. according to each described resilience countermeasure position specific device of claim 92～96, the afore mentioned rules data, for above-mentioned selection zone each, at least one of shape and stress distribution is different.
98. a resilience countermeasure position specific program makes computer carry out following steps:
The molding condition of numeric value analysis drawing product is calculated the drawing analyzing step of the one-tenth graphic data of drawing product;
At least a portion zones among the one-tenth graphic data of selected drawing product, drawing product, calculate with the stress distribution of the shape of this selection area and this selection area at least one with specified data replace form be replaced as graphic data be replaced as graphic data generation step; With
The above-mentioned graphic data that is replaced as of numeric value analysis is calculated the springback capacity of springback capacity and is calculated step.
99. according to the described resilience countermeasure of claim 98 position specific program, change the position and the quantity of above-mentioned selection area, repeat above-mentionedly to be replaced as that graphic data generates step and above-mentioned springback capacity is calculated step.
100., also have and judge whether above-mentioned springback capacity becomes the following step of setting according to claim 98 or 99 described resilience countermeasure position specific programs.
101. according to the described resilience countermeasure of claim 100 position specific program, the above-mentioned graphic data that is replaced as generates step, at least one specified data displacement with this displacement of the stress distribution of the shape in one or more zones of the springback capacity below the setting and this selection area will be provided, and then change the position or the quantity of above-mentioned selection area, repeat above-mentionedly to be replaced as that graphic data generates step and above-mentioned springback capacity is calculated step.
102. according to each described resilience countermeasure position specific program of claim 98～101, also have the step that the one-tenth graphic data of above-mentioned drawing product is divided with the zone of regulation,
The above-mentioned graphic data that is replaced as generates step, and the Zone Full for the zone of dividing with the zone of afore mentioned rules generates the above-mentioned graphic data that is replaced as.
103. according to each described resilience countermeasure position specific program of claim 98～102, the afore mentioned rules data, for above-mentioned selection zone each, at least one of shape and stress distribution is different.
104. resilience countermeasure position specific program, it is the molding condition of numeric value analysis drawing, calculate the drawing product the one-tenth graphic data drawing analysis program and numeric value analysis this become graphic data, the resilience countermeasure position specific program of calculating the resilience analysis program of springback capacity and can carry out the data input and output, it is characterized in that, make computer carry out following steps:
Obtain into the step of graphic data by above-mentioned drawing analysis program;
At least a portion zones among the one-tenth graphic data of selected drawing product, selected drawing product, calculate with the stress distribution of the shape of this selection area and this selection area at least one with specified data replace form be replaced as graphic data be replaced as graphic data generation step; With
Be replaced as the step that graphic data exports above-mentioned resilience analysis program to above-mentioned.
105. a computer-readable recording medium is characterized in that, records each described program of claim 98～104.
CN200780032577.9A 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 Active CN101511502B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410204754.9A CN103995927B (en) 2006-08-31 2007-08-30 Resilience countermeasure position specifying method and resilience countermeasure position specific device

Applications Claiming Priority (13)

Application Number Priority Date Filing Date Title
JP2006236314A JP4724626B2 (en) 2006-08-31 2006-08-31 Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP236314/2006 2006-08-31
JP244351/2006 2006-09-08
JP2006244351A JP4418449B2 (en) 2006-09-08 2006-09-08 Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP269253/2006 2006-09-29
JP268855/2006 2006-09-29
JP2006268855A JP4308239B2 (en) 2006-09-29 2006-09-29 Springback influence display method, apparatus thereof, and program thereof
JP2006269253A JP4344375B2 (en) 2006-09-29 2006-09-29 Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP2006270693A JP4418451B2 (en) 2006-10-02 2006-10-02 Springback countermeasure position specifying method, apparatus, and program thereof
JP271097/2006 2006-10-02
JP270693/2006 2006-10-02
JP2006271097A JP4418452B2 (en) 2006-10-02 2006-10-02 Springback occurrence cause identifying method, apparatus thereof, and program thereof
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

Related Child Applications (2)

Application Number Title Priority Date Filing Date
CN201410204754.9A Division CN103995927B (en) 2006-08-31 2007-08-30 Resilience countermeasure position specifying method and resilience countermeasure position specific device
CN201210056854.2A Division 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

Publications (2)

Publication Number Publication Date
CN101511502A true CN101511502A (en) 2009-08-19
CN101511502B CN101511502B (en) 2013-03-27

Family

ID=39238834

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200780032577.9A Active 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

Country Status (2)

Country Link
JP (1) JP4724626B2 (en)
CN (1) CN101511502B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102968524A (en) * 2012-11-06 2013-03-13 西北工业大学 Modeling method for two-dimensional variable-curvature process model of section bar part
CN103260781A (en) * 2010-12-14 2013-08-21 杰富意钢铁株式会社 Method of designing material for cylinder formation process, and cylinder formation processed product
CN103459058A (en) * 2011-03-30 2013-12-18 杰富意钢铁株式会社 Press forming analysis method
CN103747890A (en) * 2011-08-22 2014-04-23 杰富意钢铁株式会社 Method and device for confirming press molded product springback countermeasure effect
CN104245174A (en) * 2012-03-26 2014-12-24 杰富意钢铁株式会社 Press-forming analysis method
TWI492797B (en) * 2013-03-13 2015-07-21 China Steel Corp Cap Type Stamping Forming Die Mold Device and Cap Type Stamping Forming Mold
CN111432952A (en) * 2017-12-07 2020-07-17 杰富意钢铁株式会社 Rebound quantity variation reason position determination method
CN111801175A (en) * 2018-03-28 2020-10-20 杰富意钢铁株式会社 Method for designing press-molded article, press-molding die, press-molded article, and method for manufacturing press-molded article
CN113365751A (en) * 2019-02-01 2021-09-07 杰富意钢铁株式会社 Method for determining springback value variation factor part
CN113727790A (en) * 2019-04-25 2021-11-30 杰富意钢铁株式会社 Method and device for determining main cause part of springback value deviation
CN114505374A (en) * 2022-02-25 2022-05-17 北京计算机技术及应用研究所 Real-time following bending and springback eliminating method applied to automatic bending

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4410833B2 (en) 2007-12-25 2010-02-03 新日本製鐵株式会社 Springback generation cause analysis method, apparatus, program and recording medium
TWI438041B (en) 2008-09-30 2014-05-21 Nippon Steel & Sumitomo Metal Corp Forming simulation method, forming simulation apparatus, forming simuration program and media for storing the program, and forming method using the simulation result
CN102264486B (en) 2008-12-25 2013-11-06 新日铁住金株式会社 Springback occurrence cause analyzing method, springback occurrence cause analyzing device, springback occurrence cause analyzing program, and recording medium
JP5737059B2 (en) * 2011-08-22 2015-06-17 Ｊｆｅスチール株式会社 Press forming simulation analysis method and apparatus
JP6060591B2 (en) 2012-09-26 2017-01-18 Ｊｆｅスチール株式会社 Method and apparatus for identifying springback factor of press-formed product
US10161892B2 (en) 2013-02-08 2018-12-25 Jfe Steel Corporation Method of analyzing press forming
JP6888703B1 (en) * 2020-02-14 2021-06-16 Ｊｆｅスチール株式会社 Method for predicting shape change of press-molded products
JP7243657B2 (en) * 2020-02-14 2023-03-22 Ｊｆｅスチール株式会社 Shape change prediction method for press-formed products
JP7342921B2 (en) * 2021-09-29 2023-09-12 Ｊｆｅスチール株式会社 Method for predicting shape change of press-formed products

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100390017B1 (en) * 1994-07-08 2004-02-11 가부시키가이샤 아마다 Bending process of press brake and press brake used in the method
JP2003033828A (en) * 2001-07-23 2003-02-04 Toyota Motor Corp Method and program for forming die model
JP4000941B2 (en) * 2001-09-10 2007-10-31 スズキ株式会社 Die shape optimization support system and optimization support program in plastic working
JP4352658B2 (en) * 2002-05-29 2009-10-28 Ｊｆｅスチール株式会社 Springback analysis method for press-formed products
CN1181332C (en) * 2003-03-07 2004-12-22 江苏大学 Edge pressing force and punching speed changeable single moving sheet hydraulic test machine
JP4894294B2 (en) * 2006-02-27 2012-03-14 Ｊｆｅスチール株式会社 Press forming analysis method
JP4739147B2 (en) * 2006-08-28 2011-08-03 トヨタ自動車株式会社 Shape defect factor identification method, apparatus and program

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103260781B (en) * 2010-12-14 2015-07-15 杰富意钢铁株式会社 Method of designing material for cylinder formation process, and cylinder formation processed product
CN103260781A (en) * 2010-12-14 2013-08-21 杰富意钢铁株式会社 Method of designing material for cylinder formation process, and cylinder formation processed product
CN103459058A (en) * 2011-03-30 2013-12-18 杰富意钢铁株式会社 Press forming analysis method
CN103459058B (en) * 2011-03-30 2015-11-25 杰富意钢铁株式会社 Analysis of stamping forming method
CN103747890B (en) * 2011-08-22 2016-08-17 杰富意钢铁株式会社 The recoil countermeasure effect confirmation method of stamping product and device
CN103747890A (en) * 2011-08-22 2014-04-23 杰富意钢铁株式会社 Method and device for confirming press molded product springback countermeasure effect
CN104245174B (en) * 2012-03-26 2016-04-20 杰富意钢铁株式会社 Compressing analytic method
CN104245174A (en) * 2012-03-26 2014-12-24 杰富意钢铁株式会社 Press-forming analysis method
CN102968524B (en) * 2012-11-06 2016-02-24 西北工业大学 A kind of modeling method of two-dimentional variable curvature section bar part process model
CN102968524A (en) * 2012-11-06 2013-03-13 西北工业大学 Modeling method for two-dimensional variable-curvature process model of section bar part
TWI492797B (en) * 2013-03-13 2015-07-21 China Steel Corp Cap Type Stamping Forming Die Mold Device and Cap Type Stamping Forming Mold
US11221272B2 (en) 2017-12-07 2022-01-11 Jfe Steel Corporation Springback variation cause analysis method
CN111432952A (en) * 2017-12-07 2020-07-17 杰富意钢铁株式会社 Rebound quantity variation reason position determination method
CN111432952B (en) * 2017-12-07 2022-01-18 杰富意钢铁株式会社 Rebound quantity variation reason position determination method
CN111801175A (en) * 2018-03-28 2020-10-20 杰富意钢铁株式会社 Method for designing press-molded article, press-molding die, press-molded article, and method for manufacturing press-molded article
CN113365751A (en) * 2019-02-01 2021-09-07 杰富意钢铁株式会社 Method for determining springback value variation factor part
US11911816B2 (en) 2019-02-01 2024-02-27 Jfe Steel Corporation Method for identifying variation factor portion of springback amount
CN113727790A (en) * 2019-04-25 2021-11-30 杰富意钢铁株式会社 Method and device for determining main cause part of springback value deviation
CN113727790B (en) * 2019-04-25 2023-03-28 杰富意钢铁株式会社 Method and device for determining main factor part of springback value deviation
CN114505374A (en) * 2022-02-25 2022-05-17 北京计算机技术及应用研究所 Real-time following bending and springback eliminating method applied to automatic bending
CN114505374B (en) * 2022-02-25 2024-08-02 北京计算机技术及应用研究所 Real-time follow-up bending and rebound elimination method applied to automatic bending

Also Published As

Publication number Publication date
JP2008055476A (en) 2008-03-13
CN101511502B (en) 2013-03-27
JP4724626B2 (en) 2011-07-13

Similar Documents

Publication Publication Date Title
CN101511502B (en) Springback occurrence cause identifying method, springback influence degree display method, springback occurrence cause portion identifying method, springback measure position specifying method thereof
CN102708222B (en) 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
CN101657712B (en) breakage prediction method, calculation processing device
JP4894294B2 (en) Press forming analysis method
KR101368108B1 (en) Springback occurrence cause analyzing method, springback occurrence cause analyzing device, computer readable recording medium for recording springback occurrence cause analyzing program
US20110246150A1 (en) Method, device, program, and recording medium of analyzing cause of springback
KR20090046963A (en) Shape defect factor identification method, device, and program
JP4418452B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
JP4344375B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
Zhang et al. Selective element fission approach for fast FEM simulation of incremental sheet forming based on dual-mesh system
Ingarao et al. Optimization methods for complex sheet metal stamping computer aided engineering
JP2009148838A (en) Method of specifying springback occurrence cause portion, its device and its program
JP2011183417A (en) Method of evaluating stability of spring back
Corrado et al. Methods of influence coefficients to evaluate stress and deviation distributions of parts under operating conditions-a review
JP4308239B2 (en) Springback influence display method, apparatus thereof, and program thereof
Naranje et al. Knowledge-based system for design of deep drawing die for axisymmetric parts
CN105487488B (en) A kind of progressive die Complex Panels Sheet Metal Forming Technology sequential program(me) method
JP4418449B2 (en) Springback occurrence cause identifying method, apparatus thereof, and program thereof
Chang et al. A novel optimal working process in the design of parametric progressive dies
JP2002530197A (en) A Modeling Method for Anisotropic Metal Sheet Forming
Jonsson et al. A case-based reasoning method including tooling function for case retrieval and reuse in stamping tooling design
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
Gakwaya et al. ALE formulation and simulation techniques in integrated computer aided design and engineering system with industrial metal forming applications
Neto et al. Applying Nagata patches in the description of smooth tool surfaces used in sheet metal forming simulations
Venkatachalam et al. Influence of hole size, hole shape and hole pattern on spring-back effect in perforated sheet metals using FEM

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee

Owner name: NIPPON STEEL + SUMITOMO METAL CORPORATION

Free format text: FORMER NAME: SHIN NIPPON STEEL LTD.

CP01 Change in the name or title of a patent holder

Address after: Tokyo, Japan, Japan

Patentee after: Nippon Steel Corporation

Address before: Tokyo, Japan, Japan

Patentee before: Nippon Steel Corporation

CP01 Change in the name or title of a patent holder

Address after: Tokyo, Japan, Japan

Patentee after: Nippon Iron & Steel Corporation

Address before: Tokyo, Japan, Japan

Patentee before: Nippon Steel Corporation

CP01 Change in the name or title of a patent holder