CN114496124A - Method for measuring parameters of GISSMO material failure model under high-speed working condition - Google Patents
Method for measuring parameters of GISSMO material failure model under high-speed working condition Download PDFInfo
- Publication number
- CN114496124A CN114496124A CN202210112825.7A CN202210112825A CN114496124A CN 114496124 A CN114496124 A CN 114496124A CN 202210112825 A CN202210112825 A CN 202210112825A CN 114496124 A CN114496124 A CN 114496124A
- Authority
- CN
- China
- Prior art keywords
- speed
- gissmo
- strain
- curve
- test
- 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.)
- Pending
Links
- 239000000463 material Substances 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012360 testing method Methods 0.000 claims abstract description 46
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 230000003068 static effect Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004080 punching Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 229910000831 Steel Inorganic materials 0.000 abstract description 6
- 239000010959 steel Substances 0.000 abstract description 6
- 239000010960 cold rolled steel Substances 0.000 abstract description 4
- 238000009864 tensile test Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to a method for measuring parameters of a GISSMO material failure model under a high-speed working condition, which takes full consideration of the influence of high-speed strain generated under the action of impact on ultrahigh-strength cold-rolled steel and hot-formed steel on the mechanical property of the material in the collision process, measures a tension failure front force-displacement curve of the material at a certain strain rate by using a high-speed stretching machine, and obtains failure equivalent plastic strain of the material under different stress triaxial degrees and lode angles by combining simulation results and test data. The invention mainly solves the problem that the existing material failure model has lower precision under high-speed load generated in the collision process, does not need a complex clamp, only needs to punch a hole on a sample, and is matched with a simple bolt clamp, so that the test difficulty and the cost are reduced.
Description
Technical Field
The invention relates to the technical field of materials, in particular to a method for measuring parameters of a GISSMO material failure model under a high-speed working condition.
Background
In recent years, the requirement of consumers on the collision performance of passenger cars is higher and higher, so that the application proportion of high-strength low-ductility steel such as ultrahigh-strength cold-rolled steel, hot-formed steel and the like in car bodies is greatly improved, and the materials can be broken in the collision process due to poor plasticity. In order to accurately predict the failure behavior of the automobile material in the collision process by using a finite element simulation technology, a generalized incremental stress state related damage model (GISSMO) is developed, and the fracture behavior of the isotropic ductile material under a plurality of load paths can be accurately reproduced by applying the criterion.
Currently, the GISSMO failure model parameters are determined by measuring the force-displacement curves of the material under different stress states in a static test mode. In the actual collision process, the concerned part is usually in a high-speed deformation working condition, and the test under the static load does not consider the change of the mechanical property of the material under the high-speed load, so that the precision of a material failure model is reduced.
Disclosure of Invention
Aiming at the problems, the invention aims to provide a method for measuring parameters of a GISSMO material failure model under a high-speed working condition, and particularly aims to solve the problem that the existing material failure model is low in precision under a high-speed load generated in a collision process. This patent need not complicated anchor clamps, only needs punch the processing on the sample, cooperates simple and easy bolt anchor clamps, reduces experimental degree of difficulty and cost.
The technical scheme adopted by the invention is as follows:
the invention provides a method for measuring parameters of a GISSMO material failure model under a high-speed working condition, which specifically comprises the following steps:
s1, measuring the static stretching and high-speed stretching performance of the material, obtaining a force-displacement curve, converting the obtained force-displacement curve into a true stress-strain curve, and extrapolating the true stress-strain curve to 1 at different strain rates; inputting the true stress-strain curves of the material at different strain rates into a simulation model, and correcting the material curve by a simulation benchmarking method;
s2, punching a hole on the required tensile sample, wherein the hole is required to pass through a clamp bolt so as to be loaded on a test device; measuring force-displacement curves of various fracture tests of the material under corresponding high-speed load working conditions by using a high-speed stretcher and a cup-shaped protrusion testing machine;
s3, inputting the corrected material data into a fracture mechanics test simulation model, and combining the simulation result and the test data to obtain the corresponding equivalent plastic strain of the material under different stress triaxial degrees and lode angles; fitting data points by using an interpolation method to obtain a GISSMO failure model parameter curve of the material under a high-speed working condition, and manufacturing a material failure card by combining a tensile curve of the material;
s4, loading the broken material card in the simulation model of the steps S1 and S3 for calculation, and correcting material data by combining the simulation result and the test data until the precision of the material failure model meets the requirement.
Further, in step S2, the diameter range of the hole on the tensile sample is
Further, in the step S2, the strain rate of the high-speed stretching material is 0.001-1000/S.
Compared with the prior art, the invention has the following beneficial effects:
the GISSMO material failure model is mostly applied to collision simulation calculation to predict the fracture failure of high-strength low-ductility steel such as ultrahigh-strength cold-rolled steel or hot-formed steel in collision. Compared with the traditional GISSMO material failure model parameter measurement method, the method provided by the invention fully considers the high-speed strain generated by impact during the collision process of the ultra-high-strength cold-rolled steel and the hot-formed steel. And measuring a force-displacement curve before the material fails under tension at a certain strain rate by using a high-speed stretcher, and acquiring failure equivalent plastic strain of the material under different stress triaxial degrees and lode angles by combining simulation results and test data.
The data of the material shear test, the tensile shear test, the R5 tensile test, the R20 tensile test, the center hole tensile test, the cupping test and the like obtained by the method are more fit with the actual working condition, the accuracy of the drawn material failure curve is higher on the basis, and the collision simulation analysis accuracy is improved.
Drawings
FIG. 1 is a flow chart of a method for measuring parameters of a GISSMO material failure model under a high-speed working condition, which is provided by the invention;
FIG. 2 is a schematic illustration of a center hole tensile specimen;
FIG. 3 is a schematic representation of an R5 tensile specimen;
FIG. 4 is a schematic representation of an R20 tensile specimen;
FIG. 5 is a schematic view of a shear tensile specimen;
FIG. 6 is a schematic drawing of a tension-shear tensile specimen;
FIG. 7 is a schematic diagram of the cupping test.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
The method for measuring the parameters of the GISSMO material failure model under the high-speed working condition, disclosed by the invention, comprises the following steps of:
step S1, measuring the static stretching and high-speed stretching performance of the material, converting the obtained force-displacement curve into a true stress-strain curve in a data processing mode, and fitting and extrapolating to 1; inputting the true stress-strain curves of the material at different strain rates into a simulation model, calculating, calibrating the simulation calculation result and the test data, and if the calibration result is not good, correcting the material curve and repeating the step;
step S2, punching a hole on the required tensile sample, wherein the hole can pass through a clamp bolt so as to be loaded on the test device, and the diameter range of the hole is
Measuring force-displacement curves of various fracture tests of the material under corresponding high-speed load working conditions by using a high-speed stretcher and a cup-shaped protrusion testing machine; strain rate of high speed stretch materialThe specific value of the rate is between 0.001 and 1000/s;
step S3, inputting the corrected material data into various fracture mechanics test finite element simulation models, and combining the simulation result and the test data to obtain the corresponding equivalent plastic strain of the material under different stress triaxial degrees and lode angles; fitting data points by using an interpolation method to obtain a GISSMO failure model parameter curve of the material under a high-speed working condition, and establishing a corresponding material card by combining a tensile curve and a fracture failure curve of the material;
and step S4, loading the material card into the simulation model used in the steps S1 and S3 for calculation, and correcting the material data by combining the simulation result and the test data until the precision meets the requirement.
The tensile simulation sample and the test sample need to be identical, and the simulation boundary condition and the working condition are subjected to a reference test.
The simulation model required in step S1 at least includes a static tensile test model, and if there is a high-speed tensile data requirement, a high-speed tensile test model needs to be added.
The samples required for step S2 include at least standard uniaxial tensile, high speed tensile, shear, tensile shear, R5 tensile, R20 tensile, center hole tensile and cupping.
The finite element simulation model used in step S3 includes a shear specimen, a tension-shear specimen, an R5 tensile specimen, an R20 tensile specimen, a center hole tensile specimen, and a cupping specimen.
The strain rate of the high-speed stretching material is selected according to actual requirements and is not higher than 1000/s.
Determining failure force-displacement curves of different failure modes of the material under a certain strain rate and positions of failure on the sample through high-speed tensile tests of different samples, extracting stress triaxial degrees, lode angles and equivalent plastic strains of corresponding failure time and positions in a simulation model, and drawing the failure curve of the material by using an interpolation method.
Preparing a uniaxial tensile sample and a high-speed tensile sample of the material to be tested, measuring true stress-strain curves of the material at different strain rates through a static tensile testing machine and a high-speed tensile testing machine, and extrapolating to finally obtain an extrapolated mechanical property curve of the material for collision simulation analysis. And establishing a simulation model according to the static and high-speed tensile sample and the test conditions of the material. And inputting the extrapolated true stress-strain curves of the material at different strain rates into a simulation model, comparing the simulation result with the test data, and continuously correcting the curve extrapolation part until the simulation result is basically consistent with the test data.
As shown in FIGS. 2 to 6, a shear specimen, a tension-shear specimen, an R5 tensile specimen, an R20 tensile specimen and a center hole tensile specimen were prepared, respectively, and the specimens were drilled with a drill
The holes of (a) are matched with the bolts of the fixture, so that key data points in a fracture curve required for establishing a GISSMO material failure model can be covered. Performing a tensile test on the sample, selecting a proper strain rate for testing, wherein the strain rate range is 0.001-1000/s, and obtaining a force-displacement curve before failure and a failure part of the material at a certain strain rate; simultaneously preparing a cupping test sample, and finishing the cupping test as shown in fig. 7; establishing a simulation model according to the test sample and the test conditions, and requiring the simulation conditions to be aligned with the standard test conditions; extracting the time and the position of the displacement equal to the test failure displacement in the simulation result, and checking the stress triaxial degree, the lode angle and the corresponding equivalent plastic strain of the corresponding unit; and processing the obtained data points by using an interpolation method to prepare a material fracture characteristic curve required by the GISSMO material failure model.
Adding MAT _ ADD _ EROSION, namely a GISSMO failure model, into a data model of a material static uniaxial tensile test, a high-speed tensile test, a shearing test, a pulling-shearing test, an R5 tensile test, an R20 tensile test, a center hole tensile test and a cup drawing test, comparing errors of a force-displacement curve output by simulation with test data, and calculating the mean square error of the simulation and test force-displacement curves. And (3) correcting the fracture characteristic curve of the material, optimizing the GISSMO model, and repeating the process. And continuously iterating parameters to reduce the mean square error of the simulation and test force-displacement curve until the precision requirement is met.
The above-mentioned embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solution of the present invention by those skilled in the art should fall within the protection scope defined by the claims of the present invention without departing from the spirit of the present invention.
Claims (3)
1. A method for measuring parameters of a GISSMO material failure model under a high-speed working condition is characterized by comprising the following steps:
s1, measuring the static stretching and high-speed stretching performance of the material, obtaining a force-displacement curve, converting the obtained force-displacement curve into a true stress-strain curve, and extrapolating the true stress-strain curve to 1 at different strain rates; inputting the true stress-strain curves of the material at different strain rates into a simulation model, and correcting the material curve by a simulation benchmarking method;
s2, punching a hole on the required tensile sample, wherein the hole is required to pass through a clamp bolt so as to be loaded on a test device; measuring force-displacement curves of various fracture tests of the material under corresponding high-speed load working conditions by using a high-speed stretcher and a cup-shaped protrusion testing machine;
s3, inputting the corrected material data into a fracture mechanics test simulation model, and combining the simulation result and the test data to obtain the corresponding equivalent plastic strain of the material under different stress triaxial degrees and lode angles; fitting data points by using an interpolation method to obtain a GISSMO failure model parameter curve of the material under a high-speed working condition, and manufacturing a material failure card by combining a tensile curve of the material;
s4, loading the broken material card in the simulation model of the steps S1 and S3 for calculation, and correcting material data by combining the simulation result and the test data until the precision of the material failure model meets the requirement.
2. The method for measuring the parameters of the GISSMO material failure model under the high-speed working condition according to claim 1, wherein: in step S2, the diameter range of the hole on the tensile sample is
3. The method for measuring the parameters of the GISSMO material failure model under the high-speed working condition according to claim 1, wherein: in the step S2, the strain rate of the high-speed stretching material is 0.001-1000/S.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210112825.7A CN114496124A (en) | 2022-01-29 | 2022-01-29 | Method for measuring parameters of GISSMO material failure model under high-speed working condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210112825.7A CN114496124A (en) | 2022-01-29 | 2022-01-29 | Method for measuring parameters of GISSMO material failure model under high-speed working condition |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114496124A true CN114496124A (en) | 2022-05-13 |
Family
ID=81479419
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210112825.7A Pending CN114496124A (en) | 2022-01-29 | 2022-01-29 | Method for measuring parameters of GISSMO material failure model under high-speed working condition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114496124A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115828459A (en) * | 2022-11-30 | 2023-03-21 | 常州市新创智能科技有限公司 | R angle failure mode control method for interlayer eccentric structure joint |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112800645A (en) * | 2021-01-05 | 2021-05-14 | 唐山钢铁集团有限责任公司 | GISSMO material failure model parameter measurement method |
-
2022
- 2022-01-29 CN CN202210112825.7A patent/CN114496124A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112800645A (en) * | 2021-01-05 | 2021-05-14 | 唐山钢铁集团有限责任公司 | GISSMO material failure model parameter measurement method |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115828459A (en) * | 2022-11-30 | 2023-03-21 | 常州市新创智能科技有限公司 | R angle failure mode control method for interlayer eccentric structure joint |
| CN115828459B (en) * | 2022-11-30 | 2024-04-30 | 常州市新创智能科技有限公司 | R-angle failure mode control method for joint of sandwich eccentric structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111060396B (en) | Material mechanical property calibration method based on Ls-Dyna No. 187 material card | |
| US20020077795A1 (en) | System, method and storage medium for predicting impact performance of thermoplastic | |
| CN109933925B (en) | Method for predicting stamping forming performance of metal plate | |
| CN110987621A (en) | Method for establishing three-dimensional fracture model of metal material in complex stress state | |
| JP2020046184A (en) | Stress-strain relationship estimation method | |
| CN114496124A (en) | Method for measuring parameters of GISSMO material failure model under high-speed working condition | |
| CN101487780A (en) | Method and apparatus for detecting non-work directional performance of medium plate | |
| CN110749510B (en) | Method for detecting bending property of metal material based on finite element simulation | |
| CN110793853B (en) | Tension-torsion steady-state cyclic stress-strain modeling method based on basic mechanical parameters | |
| CN112800645A (en) | GISSMO material failure model parameter measurement method | |
| CN109716100B (en) | Material sample and method for determining the geometry of a sample | |
| CN109490334B (en) | Nondestructive testing method for T-shaped forge piece by using residual stress prediction model | |
| Evans et al. | Numerical modelling of small disc creep test | |
| CN113033040B (en) | Accurate modeling method for vehicle flexible connection | |
| CN112763345B (en) | Airplane airfoil load calibration area loading test method | |
| CN113218604B (en) | Device and method for measuring local rigidity of door limiter hole | |
| Uemori et al. | Elasto-plasticity behavior of high strength steel sheet in biaxial stress path change | |
| CN114018733B (en) | High-temperature fatigue crack propagation test method and device based on flexibility principle | |
| KR20160055630A (en) | Apparatus for safety assessement of glass materials and evaluation method thereof | |
| Panwitt et al. | Automated crack length measurement for mixed mode fatigue cracks using digital image correlation | |
| KR102521067B1 (en) | Evaluation system of directionality of material yield strength using spherical indentation and surface displacement analysis | |
| CN116659894A (en) | Method and system for testing loading force of automobile chassis part bench endurance test | |
| CN215492316U (en) | Door limiter hole local rigidity measuring device | |
| CN115512797A (en) | Method for analyzing strength of part based on material failure model | |
| Aydin et al. | On the study of constitutive parameter identification of advanced yield criteria |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |