CN104785601A - Sheet stamping method - Google Patents
Sheet stamping method Download PDFInfo
- Publication number
- CN104785601A CN104785601A CN201510193879.0A CN201510193879A CN104785601A CN 104785601 A CN104785601 A CN 104785601A CN 201510193879 A CN201510193879 A CN 201510193879A CN 104785601 A CN104785601 A CN 104785601A
- Authority
- CN
- China
- Prior art keywords
- time point
- sensitization
- sheet
- wrinkling
- group
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21D—WORKING OR PROCESSING OF SHEET METAL OR METAL TUBES, RODS OR PROFILES WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21D22/00—Shaping without cutting, by stamping, spinning, or deep-drawing
- B21D22/02—Stamping using rigid devices or tools
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
The invention relates to the field of design of the sheet stamping technology, in particular to a sheet stamping method. The sheet stamping method can conduct variable blank holder force optimization on areas where sensitive changes happen and corresponding moments in the sheet stamping and forming process. According to the sheet stamping method, according to a stamping die and a three-dimensional model of a sheet, an finite element analyzing method is used for conducting meshing, the safe changing range of blank holder force of the sheet is determined through numerical simulation, the time needed for stamping and forming is equally divided into multiple parts, Latin square experiment design is applied, k sets of data including equally-dividing time points and blank holder force values corresponding to the dividing time points, k sets of blank holder force changing curves are obtained, k sets of experiment data are obtained through numerical simulation, k sets of sensitization time points of the sheet are obtained so that the final sensitization time point can be determined, a time point is added at the final sensitization time point to serve as a new sample point, the new sample point is optimized through a firefly algorithm and a Gaussian process model, the optimal blank holder force changing curve is obtained, and finally the sheet stamping process is completed.
Description
Technical field
The present invention relates to sheet metal stamping process design field, refer to a kind of sheet-metal press working method especially.
Background technology
In sheet metal forming process, pressure-pad-force is a very important technological parameter, is mainly used to produce frictional force, to increase the tension of plate, controls the flowing of plate.The size of pressure-pad-force affect plate produce the degree of plastic deformation and break, the problem such as wrinkling, the success or failure of drawing and forming and the quality of stamping parts are directly connected to the control of pressure-pad-force.
Pressure variable edge force technology is the technology being provided the pressure-pad-force with stroke or time variations by adjustable pressure or hydraulic means, effectively can control the Plastic Flow of plate.But determine that optimum pressure-pad-force change curve is still more difficult accurately at present.The method of existing acquisition pressure variable edge force has test method(s) and numerical simulation method.Test method(s) carries out tension test to multiple forecast model, is drawn by mass data described point with the figure obtaining this kind of construction instability region, and determines optimum pressure-pad-force change curve in conjunction with stamping parts actual size.This method needs to carry out the shock pressure experiments under a series of different flanging, and efficiency and precision are all lower.Numerical simulation method is combined optimized algorithm and finite element software, use the variation tendency of optimal theoretical determination pressure-pad-force, the relation of forming quality and pressure-pad-force is evaluated by Finite Element Numerical Simulation, the method of this acquisition pressure variable edge force is carried out time of whole stamping process or stroke, do not have specific aim.These method spended times obtaining pressure variable edge force are oversize at present, and efficiency is too low, and optimal way is single, and the pressure-pad-force change curve obtained also needs to be improved further in raising sheet forming precision, and effect of optimization is undesirable.
Summary of the invention
In view of this, the moment that the object of the invention is to propose a kind of region and correspondence for there is sensitivity change in sheet metal forming process carries out the sheet-metal press working method of variable Blank Holder Force Optimization.
Based on above-mentioned purpose sheet-metal press working method provided by the invention, comprise the following steps:
(1) according to the threedimensional model of diel and plate, finite element method is adopted to carry out stress and strain model, through the safe excursion of pressure-pad-force of numerical simulation determination plate;
(2) what step (1) numerical simulation obtained stamping is requiredly divided into many parts total time, obtain decile time point, use Latin square experimental design, obtain k group data, often organize the flanging force value that data comprise each decile time point and each decile time point of correspondence, according to described k group data acquisition k group pressure-pad-force change curve;
(3) k group pressure-pad-force change curve step (2) obtained obtains k group experimental data through numerical simulation, obtains the k group sensitization time point of plate in forming process, utilizes described k group sensitization time point to determine final sensitization time point;
Described sensitization time point comprises the time point that described plate wrinkling time point occurs and fracture tendency occurs;
The wrinkling time point of described generation is in forming limit diagram, and strain the time point being become corrugated regions from wrinkling trend region, the time point of described generation fracture tendency is in forming limit diagram, strains the time point become from safety zone close to burst region;
Described final sensitization time point comprises wrinkling sensitization time point and fracture tendency sensitization time point;
(4) time point is increased, as the new sample point corresponding to described decile time point with the residual quantity value determined at the one-sided or bilateral of described final sensitization time point;
(5) utilize the new sample point that step (4) obtains, adopt glowworm swarm algorithm and Gaussian process model to be optimized, obtain optimum pressure-pad-force change curve;
(6) Sheet Metal Stamping Process is completed according to described optimum pressure-pad-force change curve.
Preferably, described decile time point is 1-10, and the span of described k is 15-50.
Preferably, there is the mean value of wrinkling time point as the wrinkling sensitization time point in described final sensitization time point in the plate got in described k group sensitization time point, or gets plate in described k group sensitization time point and the maximum value of wrinkling time point occurrence number occurs as the wrinkling sensitization time point in final sensitization time point;
Get the mean value of the time point of the plate generation fracture tendency in described k group sensitization time point as the fracture tendency sensitization time point in described final sensitization time point, or the maximum value of the time point occurrence number of getting the plate generation fracture tendency in described k group sensitization time point is as the fracture tendency sensitization time point in final sensitization time point.
Preferably, increasing time point number at the one-sided or bilateral of each described sensitization time point is 1-3.
As can be seen from above, sheet-metal press working method provided by the invention, be compared with the method for optimized variable with traditional time based on whole punching course or stroke, to select in Sheet Metal Stamping Process the easily wrinkling and sensitizing range of breaking and time, and carry out variable Blank Holder Force Optimization at sensitization time, have more specific aim, more directly and effectively can optimize the wrinkling and fracture phenomena of plate in forming process.
The optimization precision of the method can also adjust as required, and the number of the sample point that can be increased in sensitization time point both sides by adjustment changes optimization precision, and wherein, the sample point forming pressure-pad-force change curve is more, optimizes precision higher.
Accompanying drawing explanation
Fig. 1 is the safe excursion schematic diagram of pressure-pad-force of the embodiment of the present invention 1 plate;
Fig. 2 is pressure-pad-force change curve schematic diagrames corresponding to the embodiment of the present invention 1 step (2) first groups of data;
Fig. 3 is the embodiment of the present invention 1 optimization method flow chart;
Fig. 4 is that the embodiment of the present invention 1 optimizes pressure-pad-force change curve schematic diagram corresponding to rear optimal solution;
Fig. 5 is that the embodiment of the present invention 2 optimizes pressure-pad-force change curve schematic diagram corresponding to rear optimal solution.
Detailed description of the invention
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.
In order to the problem that the deficiency and optimal way that overcome existing variable Blank Holder Force Optimization method effect of optimization are single, the present invention proposes sheet-metal press working method to obtain the time dependent curve of pressure-pad-force.Diel and plate import in DYNAFORM software by three-dimensional modeling by the method, adopt finite element method to carry out stress and strain model, through the rational change range of numerical simulation determination pressure-pad-force; Needed for stamping, be divided into some parts total time, the experimental design of utilization Latin square, obtains the flanging force value of many group time points and corresponding each time point, organizes pressure-pad-force change curve to obtain more; Obtained by numerical simulation and organize experimental data more, analyze simulation result to obtain, in forming process, the responsive region of change and the shaping moment of correspondence occur in plate; Increase sample point in sensitization time region, re-start experimental design, with Optimal Experimental result; Set up optimization system, in the pressure-pad-force scope of above-mentioned acquisition, hunt out optimum pressure-pad-force change curve.
Embodiment 1
The method is used for optimizing in the punch forming optimization of box part by the present embodiment, and wherein the material of plate is DP600, high strength steel.The size of plate is as shown in table 1:
Table 1 plate size
| Material | Size (mm) | Radius of corner (mm) | Sheet thickness (mm) | Size of mesh opening (mm) |
| DP600 | 300*250 | 100 | 1.0 | 4 |
Implementation step is as follows:
(1) model of diel and plate is set up by 3 d modeling software, and import in Special punching pressing formation simulation software DYNAFORM, multiple constant blank holder force value is selected to carry out stamping forming simulation to plate, analyze punching press effect and the shaping situation of plate, determine the safe excursion of the pressure-pad-force of plate, namely both can not there is the wrinkling pressure-pad-force scope also can not broken in plate, as shown in Figure 1; The safe excursion of pressure-pad-force that this plate is determined is 400kN-595kN.
(2) stamping process required time is divided into 3 parts, gets 4 time points, as shown in table 2; Then in the safe excursion of pressure-pad-force, use latin square experiment method for designing to extract the flanging force value of 30 groups of corresponding 4 time points, as shown in table 3; Form 30 pressure-pad-force change curves, wherein the pressure-pad-force change curve that obtains of first group of data as shown in Figure 2;
In the present embodiment, decile time point is 4, and can select the number waiting point time point in other cases as the case may be, range of choice is 1-10; The data utilizing latin square experiment design to obtain are 30 groups, and can select the data group number obtained in other cases as the case may be, range of choice is 15-50;
The result of latin square experiment design can be genetic to next group, and therefore, this step adopts latin square experiment method for designing when not affecting randomness, can avoid there are two groups of identical phenomenons in the curve occurring to obtain.
Tables 2 etc. divide time point
| t 1(s) | t 2(s) | t 3(s) | t 4(s) |
| 0 | 0.0352 | 0.0704 | 0.1056 |
Table 3 four sample point pressure-pad-force change curve
| Numbering | F 1(10 5N) | F 2(10 5N) | F 3(10 5N) | F 4(10 5N) |
| 1 | 5.2169 | 4.2488 | 4.5474 | 5.7778 |
| 2 | 5.2434 | 4.0178 | 5.9294 | 4.5007 |
| 3 | 5.8578 | 4.1623 | 4.6891 | 5.1334 |
| 4 | 5.7740 | 4.1851 | 5.0671 | 4.9656 |
| 5 | 5.9091 | 4.0123 | 4.6575 | 5.3498 |
| 6 | 4.7783 | 4.4280 | 5.8191 | 5.3368 |
| 7 | 5.3008 | 5.8758 | 4.3178 | 4.9157 |
| 8 | 4.5440 | 4.4415 | 5.0091 | 5.7555 |
| 9 | 4.9519 | 5.2956 | 4.1287 | 5.5626 |
| 10 | 4.2219 | 5.1280 | 4.5394 | 5.9233 |
| 11 | 5.0553 | 4.8178 | 5.6808 | 4.0102 |
| 12 | 5.1947 | 5.7889 | 4.1855 | 4.7992 |
| 13 | 4.6062 | 5.2607 | 5.7100 | 5.1022 |
| 14 | 5.4758 | 4.3277 | 4.5666 | 5.1022 |
| 15 | 4.0223 | 5.4469 | 5.7760 | 4.6520 |
| 16 | 4.8378 | 5.8377 | 5.1158 | 4.1830 |
| 17 | 4.1996 | 5.1405 | 5.9268 | 4.8050 |
| 18 | 4.2675 | 5.8326 | 5.1139 | 4.5576 |
| 19 | 5.8128 | 5.1061 | 4.4203 | 4.5671 |
| 20 | 4.4199 | 4.6882 | 5.7715 | 5.0692 |
| 21 | 4.2419 | 5.2238 | 5.5224 | 4.8028 |
| 22 | 5.2434 | 5.4803 | 4.9544 | 4.0132 |
| 23 | 4.8828 | 5.6248 | 4.2016 | 5.1334 |
| 24 | 4.3115 | 5.6476 | 5.0671 | 4.9656 |
| 25 | 4.9341 | 4.0123 | 5.1450 | 5.8373 |
| 26 | 5.2658 | 4.9155 | 4.3566 | 5.8243 |
| 27 | 5.7883 | 5.3883 | 4.3178 | 4.9157 |
| 28 | 5.0315 | 5.9040 | 4.0341 | 4.7805 |
| 29 | 5.9269 | 4.8081 | 5.1037 | 4.1001 |
| 30 | 4.2219 | 5.1280 | 4.5394 | 5.9233 |
(3) 30 pressure-pad-force change curves of step (2) gained are input in stamping forming simulation software DYNAFORM, sheet-metal press working effect is obtained through analogue simulation, obtain the forming limit diagram of plate, 30 groups of sensitization time points are found by stamping forming simulation result, often organize sensitization time point and comprise the time point that wrinkling time point and generation fracture tendency occur plate, simulation result is as shown in table 4;
In the forming limit diagram of plate, the change of strain comprises red (CRACK) region of breaking, yellow close (RISK OF CRACK) region of breaking, green safety (SAFE) region, blue wrinkling trend (WRINKLE TENDENCY) region and wrinkling (WRINKLE) region of pink colour, wrinkling time point occurs plate is in forming limit diagram, strain is become the time point in wrinkling (WRINKLE) region of pink colour from wrinkling trend (WRINKLETENDENCY) region of blueness, the time point of plate generation fracture tendency is in forming limit diagram, strain is become the yellow time point close to (RISK OF CRACK) region of breaking from safety (SAFE) region of green,
The list of table 4 simulation result
| Sample number | Wrinkling time of origin (s) | Fracture tendency time of origin (s) |
| 1 | 0.067887 | 0.052800 |
| 2 | 0.100000 | 0.052800 |
| 3 | 0.067887 | 0.052800 |
| 4 | 0.082972 | 0.052800 |
| 5 | 0.067887 | 0.052800 |
| 6 | 0.100000 | 0.052800 |
| 7 | 0.075429 | 0.044513 |
| 8 | 0.082972 | 0.052800 |
| 9 | 0.067887 | 0.052800 |
| 10 | 0.082972 | 0.052800 |
| 11 | 0.098059 | 0.052800 |
| 12 | 0.075429 | 0.045257 |
| 13 | 0.098059 | 0.052800 |
| 14 | 0.067887 | 0.052800 |
| 15 | 0.100000 | 0.052800 |
| 16 | 0.098059 | 0.045257 |
| 17 | 0.100000 | 0.052800 |
| 18 | 0.098059 | 0.045257 |
| 19 | 0.075429 | 0.045257 |
| 20 | 0.100000 | 0.052800 |
| 21 | 0.098059 | 0.052800 |
| 22 | 0.082972 | 0.045257 |
| 23 | 0.075429 | 0.045257 |
| 24 | 0.098059 | 0.045257 |
| 25 | 0.098059 | 0.052800 |
| 26 | 0.067887 | 0.052800 |
| 27 | 0.075429 | 0.045257 |
| 28 | 0.075429 | 0.045257 |
| 29 | 0.098059 | 0.045257 |
| 30 | 0.082972 | 0.052800 |
By time point comparative analysis wrinkling for the generation in these 30 groups of data, wrinkling time point distribution occurs and comparatively disperses, therefore can get its mean value as the wrinkling sensitization time point in final sensitization time point, then it is T
a=0.085s; By the time point comparative analysis of the generation fracture tendency in these 30 groups of data, there is the change that the time point of fracture tendency is too not large, concentrate on some set times, therefore can get this value as the fracture tendency sensitization time point in final sensitization time point, then it is T
b=0.0528s;
T
a, T
bbe final sensitization time point.
(4) after obtaining the final sensitization time point in step (3), the new sample point corresponding to decile time point is redefined, because final sensitive spot is the object of main research, so need to increase sample point at final sensitization time point place.Select difference amount Δ
t, then by T
a+ Δ
tand T
a-Δ
tdetermine two sample points, make to use the same method at T
bplace also increases by two sample points.In the operation of reality, can according to the number needing adjustment sample point of precision;
In the present embodiment, the time of serious wrinkling generation is at T
aplace's fluctuation is comparatively large, so choose Δ
tfor 0.006s; And fracture tendency time of origin Relatively centralized, so less Δ can be selected
t, get Δ
t=0.0006s; Then can obtain new time sampling point, as shown in table 5.
New sample point near table 5 sensitization time
| Sampled point | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
| Time (s) | 0 | 0.0522 | 0.0528 | 0.0534 | 0.0790 | 0.0850 | 0.0910 | 0.1056 |
(5) adopt firefly (FA) algorithm and Gaussian process (GP) model to be optimized, obtain optimum pressure-pad-force change curve.The flanging force value of the sample point that this step is determined with corresponding step (4) is for variable X
0, account for the percentage of total element number with wrinkling trend unit and fracture tendency unit for object function Y
0, use optimized algorithm to obtain optimum pressure-pad-force change curve;
Concrete steps are as follows:
[1] initial sample X is generated
0, Y
0; Wherein X
0for the flanging force value corresponding with time point, it forms pressure-pad-force curve; Y
0for object function, be the percentage that wrinkling and breaking unit accounts for the total number of unit of plate.Initial sample list is as shown in table 6;
The initial sample list of table 6
| Numbering | F 1(10 5N) | F 2(10 5N) | F 3(10 5N) | F 4(10 5N) | F 5(10 5N) | F 6(10 5N) | F 7(10 5N) | F 8(10 5N) |
| 1 | 5.8592 | 5.7876 | 5.6681 | 5.0178 | 5.8469 | 5.8205 | 4.0883 | 5.0719 |
| 2 | 4.4392 | 4.8436 | 4.7665 | 4.2480 | 5.6833 | 5.6270 | 5.4728 | 4.1495 |
| 3 | 4.8556 | 4.5611 | 4.0866 | 4.0658 | 4.8291 | 5.4368 | 5.8913 | 4.7517 |
| 4 | 5.6235 | 4.3250 | 4.6395 | 5.6224 | 5.3477 | 4.8261 | 4.8102 | 5.9497 |
| 5 | 5.2305 | 5.6628 | 5.4160 | 4.5889 | 5.0491 | 4.2977 | 5.3920 | 5.3498 |
| 6 | 4.0194 | 5.1875 | 5.1144 | 4.9487 | 4.3169 | 5.1155 | 4.3237 | 4.3841 |
| 7 | 5.2059 | 5.2190 | 5.8325 | 5.2416 | 4.5756 | 4.7091 | 4.5618 | 4.6189 |
| 8 | 4.5514 | 4.2020 | 4.4581 | 5.7356 | 4.2356 | 4.1788 | 5.2022 | 5.5185 |
| 9 | 4.1144 | 5.0066 | 4.9741 | 5.1337 | 4.2872 | 4.8968 | 4.4860 | 4.5820 |
| 10 | 4.1944 | 4.8308 | 4.8725 | 5.2357 | 4.2754 | 4.7221 | 4.6196 | 4.7670 |
| 11 | 4.2795 | 4.7073 | 4.8029 | 5.3149 | 4.2946 | 4.6209 | 4.7420 | 4.9187 |
| 12 | 4.2634 | 4.7242 | 4.7679 | 5.3470 | 4.2432 | 4.6108 | 4.7378 | 4.9436 |
| 13 | 4.3019 | 4.6464 | 4.6932 | 5.4057 | 4.2534 | 4.5164 | 4.8549 | 5.0849 |
| 14 | 4.3351 | 4.6120 | 4.7065 | 5.4243 | 4.2524 | 4.5550 | 4.8432 | 5.0281 |
| 15 | 4.2933 | 4.6324 | 4.6793 | 5.4423 | 4.2447 | 4.5414 | 4.8322 | 5.0730 |
| 16 | 4.3379 | 4.6330 | 4.7246 | 5.4030 | 4.2672 | 4.5439 | 4.8141 | 5.0766 |
| 17 | 4.3372 | 4.6043 | 4.7221 | 5.4303 | 4.2565 | 4.5262 | 4.8380 | 5.0769 |
| 18 | 4.3423 | 4.6294 | 4.7356 | 5.4540 | 4.2429 | 4.5544 | 4.8069 | 5.0473 |
Following setting is done to the element thickness obtained after sheet forming emulation: if element thickness is greater than l
1=1.1mm, then using this unit as wrinkling trend generating unit, and element thickness is less than l
2=0.88mm, then using this unit as fracture tendency generating unit; Wherein l
1, l
2value determine according to the performance of material;
[2] adopt FA algorithm and GP model to carry out iteration optimization, finally obtain the optimal value of above-mentioned sample; Method as shown in Figure 3, be exactly the point that wherein optimal solution target function value is minimum, total namely wrinkling and breaking unit accounts for the minimum point of the percentage of plate number of unit;
[3] the pressure-pad-force change curve corresponding to optimal solution is extracted, both optimum pressure-pad-force change curves, as shown in Figure 4.
Through said procedure optimization, obtain Y
0minimum of a value be 0.15157, the pressure-pad-force change curve data of its correspondence are as shown in table 7;
Flanging force value corresponding to each time point of table 7
| F 1(10 5N) | F 2(10 5N) | F 3(10 5N) | F 4(10 5N) | F 5(10 5N) | F 6(10 5N) | F 7(10 5N) | F 8(10 5N) | Y 0 |
| 4.1526 | 4.0333 | 4.6585 | 5.3916 | 4.0000 | 4.3631 | 4.9345 | 5.2659 | 0.15157 |
(6) Sheet Metal Stamping Process is completed according to the optimum pressure-pad-force change curve of gained.
Embodiment 2
Utilize as equipartition time point as sample point, adopt FA algorithm and GP model to obtain optimum pressure-pad-force change curve, as a comparison case.
This embodiment with flanging force value corresponding to the equipartition time point of step (2) gained in embodiment 1 for variable X
0, account for the percentage of total element number with wrinkling trend unit and fracture tendency unit for object function Y
0, use optimized algorithm to obtain optimum pressure-pad-force change curve;
Concrete steps are as follows:
[1] initial sample X is generated
0, Y
0; Be wherein the flanging force value corresponding with time point, it forms pressure-pad-force curve.Y
0for object function, be the percentage that wrinkling and breaking unit accounts for the total number of unit of plate.Initial sample list is as shown in table 8;
Table 8 contrasts the initial sample list of example
| Numbering | F 1(10 5N) | F 2(10 5N) | F 3(10 5N) | F 4(10 5N) |
| 1 | 4.7783 | 4.4280 | 5.8191 | 5.3368 |
| 2 | 5.3008 | 5.8758 | 4.3178 | 4.9157 |
| 3 | 4.5440 | 4.4415 | 5.0091 | 5.7555 |
| 4 | 4.9519 | 5.2956 | 4.1287 | 5.5626 |
| 5 | 4.2219 | 5.1280 | 4.5394 | 5.9233 |
| 6 | 4.2419 | 5.2238 | 5.5224 | 4.8028 |
| 7 | 5.2434 | 5.4803 | 4.9544 | 4.0132 |
| 8 | 4.8828 | 5.6248 | 4.2016 | 5.1334 |
| 9 | 4.3115 | 5.6476 | 5.0671 | 4.9656 |
| 10 | 4.9341 | 4.0123 | 4.6575 | 5.8373 |
| 11 | 5.7740 | 4.1851 | 5.0671 | 4.9656 |
| 12 | 5.1947 | 5.7889 | 4.1855 | 4.7992 |
| 13 | 4.8828 | 5.6290 | 4.2116 | 5.1334 |
| 14 | 5.3028 | 4.1587 | 4.2016 | 4.8028 |
| 15 | 4.1144 | 5.0066 | 4.9741 | 5.1337 |
| 16 | 5.3477 | 4.8261 | 4.8102 | 5.9497 |
| 17 | 4.2946 | 4.6209 | 4.7420 | 4.9187 |
| 18 | 4.2565 | 4.5262 | 4.8380 | 5.0769 |
[2] FA algorithm and GP model is adopted to carry out iteration optimization, after finally obtaining the optimal value of above-mentioned sample; Wherein optimal solution is exactly the minimum point of target function value;
[3] the pressure-pad-force change curve corresponding to optimal solution is extracted, both optimum pressure-pad-force change curves, as shown in Figure 5; Through said procedure optimization, obtain Y
0minimum of a value be that the pressure-pad-force change curve of its correspondence of 0.167796. is as shown in table 9;
Table 9 contrasts the corresponding flanging force value of the optimum pressure-pad-force curve of example
| F 1(10 5N) | F 1(10 5N) | F 1(10 5N) | F 1(10 5N) | Y 0 |
| 5.9091 | 4.0123 | 4.6585 | 5.3498 | 0.167796 |
Embodiment 1 analog result compared with this example 2 analog result, object function Y
0have dropped 0.016226, Y compared with embodiment 2 result
0reduce 9.67%, that is, after sheet metal forming, the wrinkling trend generating unit of embodiment 1 and the number sum of fracture tendency generating unit reduce, and the forming quality of plate is optimized.
Those of ordinary skill in the field are to be understood that: the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (4)
1. a sheet-metal press working method, is characterized in that, comprises the following steps:
(1) according to the threedimensional model of diel and plate, finite element method is adopted to carry out stress and strain model, through the safe excursion of pressure-pad-force of numerical simulation determination plate;
(2) what step (1) numerical simulation obtained stamping is requiredly divided into many parts total time, obtain decile time point, use Latin square experimental design, obtain k group data, often organize the flanging force value that data comprise each decile time point and each decile time point of correspondence, according to described k group data acquisition k group pressure-pad-force change curve;
(3) k group pressure-pad-force change curve step (2) obtained obtains k group experimental data through numerical simulation, obtains the k group sensitization time point of plate in forming process, utilizes described k group sensitization time point to determine final sensitization time point;
Described sensitization time point comprises the time point that described plate wrinkling time point occurs and fracture tendency occurs;
The wrinkling time point of described generation is in forming limit diagram, and strain the time point being become corrugated regions from wrinkling trend region, the time point of described generation fracture tendency is in forming limit diagram, strains the time point become from safety zone close to burst region;
Described final sensitization time point comprises wrinkling sensitization time point and fracture tendency sensitization time point;
(4) time point is increased, as the new sample point corresponding to described decile time point with the residual quantity value determined at the one-sided or bilateral of described final sensitization time point;
(5) utilize the new sample point that step (4) obtains, adopt glowworm swarm algorithm and Gaussian process model to be optimized, obtain optimum pressure-pad-force change curve;
(6) Sheet Metal Stamping Process is completed according to described optimum pressure-pad-force change curve.
2. sheet-metal press working method according to claim 1, is characterized in that, described decile time point is 1-10, and the span of described k is 15-50.
3. sheet-metal press working method according to claim 1, it is characterized in that, there is the mean value of wrinkling time point as the wrinkling sensitization time point in described final sensitization time point in the plate got in described k group sensitization time point, or gets plate in described k group sensitization time point and the maximum value of wrinkling time point occurrence number occurs as the wrinkling sensitization time point in final sensitization time point;
Get the mean value of the time point of the plate generation fracture tendency in described k group sensitization time point as the fracture tendency sensitization time point in described final sensitization time point, or the maximum value of the time point occurrence number of getting the plate generation fracture tendency in described k group sensitization time point is as the fracture tendency sensitization time point in final sensitization time point.
4. sheet-metal press working method according to claim 1, is characterized in that, increasing time point number at the one-sided or bilateral of each described sensitization time point is 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510193879.0A CN104785601B (en) | 2015-04-22 | 2015-04-22 | A kind of sheet-metal press working method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510193879.0A CN104785601B (en) | 2015-04-22 | 2015-04-22 | A kind of sheet-metal press working method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104785601A true CN104785601A (en) | 2015-07-22 |
| CN104785601B CN104785601B (en) | 2016-08-31 |
Family
ID=53551042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510193879.0A Expired - Fee Related CN104785601B (en) | 2015-04-22 | 2015-04-22 | A kind of sheet-metal press working method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104785601B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003048027A (en) * | 2001-08-03 | 2003-02-18 | Showa Denko Kk | Method for forming laminated aluminum foil container |
| JP2009279591A (en) * | 2008-05-20 | 2009-12-03 | Furukawa-Sky Aluminum Corp | Method of press-forming aluminum alloy sheet |
| JP2010064128A (en) * | 2008-09-12 | 2010-03-25 | Furukawa-Sky Aluminum Corp | Press drawing method for aluminum alloy sheet |
| CN101908090A (en) * | 2010-08-18 | 2010-12-08 | 湖南大学 | Optimization method of stamping based on space mapping of response function |
| CN102799735A (en) * | 2012-07-24 | 2012-11-28 | 湖南大学 | Springback compensation method based on technological parameter control |
| CN102950185A (en) * | 2012-11-16 | 2013-03-06 | 华中科技大学 | Variable speed and variable blank holder force linked sheet drawing process |
-
2015
- 2015-04-22 CN CN201510193879.0A patent/CN104785601B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003048027A (en) * | 2001-08-03 | 2003-02-18 | Showa Denko Kk | Method for forming laminated aluminum foil container |
| JP2009279591A (en) * | 2008-05-20 | 2009-12-03 | Furukawa-Sky Aluminum Corp | Method of press-forming aluminum alloy sheet |
| JP2010064128A (en) * | 2008-09-12 | 2010-03-25 | Furukawa-Sky Aluminum Corp | Press drawing method for aluminum alloy sheet |
| CN101908090A (en) * | 2010-08-18 | 2010-12-08 | 湖南大学 | Optimization method of stamping based on space mapping of response function |
| CN102799735A (en) * | 2012-07-24 | 2012-11-28 | 湖南大学 | Springback compensation method based on technological parameter control |
| CN102950185A (en) * | 2012-11-16 | 2013-03-06 | 华中科技大学 | Variable speed and variable blank holder force linked sheet drawing process |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104785601B (en) | 2016-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Qian et al. | 3D coupled macro–microscopic finite element modelling and simulation for combined blank-forging and rolling process of alloy steel large ring | |
| CN104573237B (en) | A kind of mould optimization method based on fretting wear CAE analysis | |
| CN106547999B (en) | Drawing die structure analysis method based on explicit dynamics finite element method | |
| CN102672059B (en) | The method of mould concave-convex mold amendment profile is determined according to emulation stamped workpieces thickness | |
| CN104729962A (en) | CH4169 alloy forging piece grain size analysis and predication method | |
| CN104077442A (en) | Method for controlling machining accuracy of large integrated thin-walled parts based on finite element analysis | |
| CN101537468B (en) | Design method of middle blank and prefabricating blank during conical shell forgeable piece forming process | |
| CN103357734A (en) | Taper blank-holding and deep-drawing stamping method | |
| CN104765912A (en) | Robustness optimizing method of aluminum plate punching process | |
| CN103871098B (en) | The structure of a kind of thickening plate (VRB) Three Dimensional Thermal forming limit diagram and using method | |
| CN104077439A (en) | Numerical simulation method of novel high-strength steel spoke drawing punching combined process | |
| CN109284539A (en) | U-shaped bellows is hydraulic or gas pressure compacting die size and process parameter optimizing algorithm | |
| CN103424318A (en) | Method for determination of forming limit curve of sheet metal | |
| CN109127945A (en) | The regulation method of the stamping precision of lightweight body outer skin | |
| CN105426629B (en) | A kind of Automobile Cover Drawing Die Strength Calculation method | |
| CN106503389A (en) | A kind of numerical method of quick determination drop stamping production line process parameter | |
| CN106202686B (en) | A kind of objective design method of turbine disk isothermal die forging process preform blank | |
| CN104636565A (en) | Magnesium alloy one-off die forging forming process optimizing method based on machinability analysis | |
| CN104785601B (en) | A kind of sheet-metal press working method | |
| Shamsi-Sarband et al. | Optimizitation of the pressure path in sheet metal hydroforming | |
| CN102722619A (en) | Method for determining material utilization rate of parts for stamping automobile covering parts | |
| CN101789031B (en) | Progressive die belt designing system based on pressure center optimization | |
| CN108971355A (en) | The corrugated method of large curved surface covering is eliminated based on gradual change gap mold | |
| Ugla et al. | An Effect of Wall Corner Radius of Punch on The Features of The Hexagonal Cup Fabricated by The Deep Drawing Process | |
| CN106294889A (en) | A kind of high-strength steel spoke centre hole flanging punch die angle of taper optimization method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160831 |