CN109408969A - A method of constitutive model is established using finite element software identification rubber viscoelastic parameter - Google Patents

A method of constitutive model is established using finite element software identification rubber viscoelastic parameter Download PDF

Info

Publication number
CN109408969A
CN109408969A CN201811263229.9A CN201811263229A CN109408969A CN 109408969 A CN109408969 A CN 109408969A CN 201811263229 A CN201811263229 A CN 201811263229A CN 109408969 A CN109408969 A CN 109408969A
Authority
CN
China
Prior art keywords
rubber
viscoelastic
parameter
ansys
model
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
CN201811263229.9A
Other languages
Chinese (zh)
Other versions
CN109408969B (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201811263229.9A priority Critical patent/CN109408969B/en
Publication of CN109408969A publication Critical patent/CN109408969A/en
Application granted granted Critical
Publication of CN109408969B publication Critical patent/CN109408969B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a kind of modes that the method present invention for establishing constitutive model using finite element software identification rubber viscoelastic parameter is combined using experiment with software, carry out rubber viscoelastic mode parameter identification, the physical model of rubbery sample is established using ANSYS, specification is designed according to size of sample used is tested.It establishes rubber superlastic and adds viscoelastic constitutive model.Experimental data combination ANSYS carries out viscoelastic parameter identification.Parameters of constitutive model is inputted, viscoelastic mode is filled according to step 3 result, carries out grid dividing, and conditions setting is arranged solving condition, carries out the calculating of the kinetic characteristics of rubber.Method before relatively, the method for the present invention is simple, and as a result accurately and reliably, the result error obtained using this method simulation analysis is less than 3%.Positive meaning is played to rubber viscoelastic parameter identification and rubber dynamical property analysis.

Description

A method of constitutive model is established using finite element software identification rubber viscoelastic parameter
Technical field
The present invention relates to a kind of methods for establishing constitutive model using finite element software identification rubber viscoelastic parameter, especially relate to And a kind of method for establishing constitutive model progress rubber kinetic characteristics calculating using finite element software identification rubber viscoelastic parameter, Belong to dynamics of solids calculating field.
Background technique
Rubber is widely applied to the row such as automobile, aviation, precision instrument because having good absorbing, damping, sealing performance In the vibration insulation structure of the machine of industry, such as: engine vibration isolation piece, fuel tank sealing ring.Rubber has complicated mechanical property, a side Super-elasticity is embodied in face, with very big adaptability to changes (ratio that stress, strain have very little), the stress-strain relation of material Indicated by strain energy density function, material proximate it is incompressible;On the one hand viscoplasticity, under ringing power, material are embodied Material response shows rate dependence, there are hysteresis effects for ess-strain.Different composition proportions can produce the different rubber of property Glue product, the reflected physical property of institute differ widely.Therefore, before rubber product is really applied in equipment, rationally Selection rubber, to rubber carry out finite element dynamic analysis, to rubber mechanical property carry out Accurate Prediction, determine product essence Degree even life of product.
Accurately establishing rubber constitutive model is the basis for carrying out rubber dynamical property analysis, and the constitutive model established is determined Can determine the emulation of limited rubber member characterize rubber kinetic characteristics comprehensively.Superlastic model is used not merely in previous research The mechanical property of the accurate comprehensive representation rubber of energy, it is meaningful indefinite using viscoelastic parameter when superlastic plus viscoelastic constitutive model, Discrimination method is complicated and applies theory deduction more, and experiment condition is confined to unidirectionally drawing, pressure more, and plane shear etc. is static and standard is quiet The disadvantages of state is tested.Which results in not comprehensive enough to rubber dynamic analysis, the rubber kinetic characteristics of prediction do not meet reality Border leads to unnecessary loss.
Summary of the invention
For the problem present on, the present invention provides a kind of application finite element software identification rubber viscoelastic parameters to establish The method that constitutive model carries out the calculating of rubber kinetic characteristics, the form that method is combined using experiment with emulation, by rubber reality It tests data and imported into progress viscoelastic parameter fitting in finite simulation element analysis software, the discrimination method is simple, and practical, it is objective to meet It is practical.A kind of new method is provided for rubber viscoelastic mode parameter identification, so that rubber dynamical property analysis is more comprehensive.
The technical solution adopted by the present invention is that a kind of application finite element software identification rubber viscoelastic parameter establishes constitutive model The method for carrying out the calculating of rubber kinetic characteristics, comprising the following steps:
1) rubbery sample physical model is established in finite element software.
2) it establishes rubber superlastic and adds viscoelastic constitutive model, wherein superlastic model uses bis- parameter model of Mooney-Rivlin It is indicated, viscoelastic mode is characterized using general Maxwell model, and general Maxwell model is expressed as three-dimensional flow in Ansys Varying model, as shown in Equation 1:
In formula: σ is cauchy stress, Pa;G (t) is shear relaxation kernel function;K (t) is volume relaxation kernel function;E is to cut to answer Become deviator;Δ is body strain;T is current time, s;τ is process time, s.Wherein:
In formula: G0And K0The respectively shearing at 0 moment and bulk relaxation module, Pa;nG、nKFor PRONY term of a series number;For relative modulus;For slack time, s.WhereinAlso referred to as shear loading coefficient; Also by for volume response coefficient.
3) experimental data that rubbery sample relaxation modulus changes over time is extracted
According to institute's DMA (Dynamic thermomechanical analysis) dynamic thermomechanical analysis system to be applied The standby rubbery sample with the same ingredient of studied rubber carries out the setting of experiment condition in conjunction with actual condition, comprising: experiment temperature Degree, test frequency scanning range, plastic strain amplitude etc..Select DMA temperature/frequency sweep experiment, obtain rubber loss modulus and Storage modulus versus time curve, by DMA tester software, the relaxation modulus data for obtaining rubber change with time Curve, the curve have fully demonstrated rubber viscous-elastic behaviour under the conditions of frequency and temperature change.
4) viscoelastic parameter identification is carried out in finite element software
The data that the rubber relaxation modulus of extraction is changed over time are input to the prony curve fitting of ANSYS In, adjustment coeff value initial parameter is 0,1, is fitted, obtains relative modulus and slack time i.e. viscoelastic in formula 2,3 Parameter.
5) it is based on finite element ANSYS, emulation project is subjected to grid dividing, boundary condition is established, sets solver Parameter carries out dynamical property analysis.
The present invention is using the mode that combines with software of experiment, progress rubber viscoelastic mode parameter identification, relatively before side For method, the method for the present invention is simple, and as a result accurately and reliably, the result error obtained using this method simulation analysis is less than 3%.It is right Rubber viscoelastic parameter identification and rubber dynamical property analysis play positive meaning.
Detailed description of the invention
Fig. 1 is method flow diagram according to the present invention.
Fig. 2 is that experimental data combination ANSYS recognizes rubber viscoelastic parameter flow chart.
Fig. 3 is to obtain rubber relaxation modulus to change over time data flowchart.
Fig. 4 is that ANSYS fitting experimental data obtains viscoelastic parameter flow chart.
Fig. 5 is application example of the present invention, and obtained rubber viscoelastic parameter is fitted by ANSYS.
Fig. 6 is rubbery sample used in present example and emulation setting example.
Fig. 7 is C01/C10With the relation curve of shore hardness;
Specific embodiment
Specific measures for implementation are provided with reference to the accompanying drawing.
The method of the invention is completed by DMA tester and ANSYS software.
The flow chart of the method for the invention as shown in Figure 1, specifically includes the following steps:
Step 1: establishing the physical model of rubbery sample using ANSYS, and specification is designed according to size of sample used is tested.
Step 2: it establishes rubber superlastic and adds viscoelastic constitutive model.
Step 3: experimental data combination ANSYS carries out viscoelastic parameter identification, and method is as shown in Figure 2.
Step 4: input Parameters of constitutive model, wherein superlastic model parameter is obtained by Fig. 7 relationship, and viscoelastic mode is according to step Rapid three result is filled in, and grid dividing is carried out, and conditions setting (as shown in Figure 6) is arranged solving condition, carries out the power of rubber Learn the calculating of characteristic.
The method that experimental data combination ANSYS carries out viscoelastic parameter identification, as shown in Figure 2:
Step 1: it obtains rubber relaxation modulus and changes with time data, method is as shown in Figure 3;
Step 2: ANSYS fitting experimental data obtains viscoelastic parameter, and method is as shown in Figure 4;
It obtains rubber relaxation modulus to change with time data, method is as shown in Figure 3:
Step 1: rubbery sample is formulated according to selected DMA experiment instrument model
Step 2: the actual working environment according to applied by rubber, the condition and range of choice experiment, including temperature, frequency, Plastic strain amplitude.
Step 3: selection DMA temperature frequency scanning experiment obtains rubber loss modulus, storage modulus, fissipation factor.
Step 4: choosing reference temperature appropriate, analyzes software by DMA and obtains the change of the relaxation modulus of rubber at any time Change data.
ANSYS fitting experimental data obtains viscoelastic parameter, and method is as shown in Figure 4:
Step 1: being deposited into notepad according to rule for the rubber relaxation modulus data of extraction, name should using English into Row name;
Step 2: entering ANSYS PRONY CURVE FITTING module, and the data that name is completed import, and adjustment is just Beginning parameter carries out data fitting;
Step 3: check fitting as a result, initial parameter is adjusted again, until the error of fitting data and experimental data exists The range of permission;
Step 4: saving and applies, and obtains rubber viscoelastic parameter.
The example that an application present invention calculates certain rubber structure kinetic characteristics is given below.
Fig. 5 is given by ANSYS fitting, experiment relaxation data and fitting relaxation data correlation curve.
Fig. 6 gives the dynamics simulation example that rubbery sample is carried out according to experiment condition.
1 viscoelastic parameter list of table
2 experimental data of table and emulation data comparison
Note: HS represents the hardness of rubber
Table 1 gives the viscoelastic parameter that rubber relaxation modulus data are obtained through ansys Fitting Analysis, and combines Fig. 5 can Show that the precision of fitting is very high.Dynamical property analysis is carried out using this group of parameter as rubber viscoelastic parameter, analysis obtains The fissipation factor and experiment gained fissipation factor of rubber are within 3%, as shown in table 2.Therefore the method for the present invention is substantially better than it Preceding method.

Claims (5)

1. a kind of method for establishing constitutive model using finite element software identification rubber viscoelastic parameter, it is characterised in that: specific packet Include following steps:
Step 1: establishing the physical model of rubbery sample using ANSYS, and specification is designed according to size of sample used is tested;
Step 2: it establishes rubber superlastic and adds viscoelastic constitutive model;
Step 3: experimental data combination ANSYS carries out viscoelastic parameter identification;
Step 4: input Parameters of constitutive model, viscoelastic mode are filled according to step 3 result, carry out grid dividing, set boundary Condition is arranged solving condition, carries out the calculating of the kinetic characteristics of rubber.
2. a kind of side for establishing constitutive model using finite element software identification rubber viscoelastic parameter according to claim 1 Method, it is characterised in that: this approach includes the following steps,
1) rubbery sample physical model is established in finite element software;
2) it establishes rubber superlastic and adds viscoelastic constitutive model, wherein superlastic model is carried out using bis- parameter model of Mooney-Rivlin It indicates, viscoelastic mode is characterized using general Maxwell model, and general Maxwell model is expressed as three-dimensional rheology mould in Ansys Type, as shown in Equation 1:
In formula: σ is cauchy stress, Pa;G (t) is shear relaxation kernel function;K (t) is volume relaxation kernel function;E is that shear strain is inclined Amount;Δ is body strain;T is current time, s;τ is process time, s;Wherein:
In formula: G0And K0The respectively shearing at 0 moment and bulk relaxation module, Pa;nG、nKFor PRONY term of a series number;For relative modulus;For slack time, s;WhereinAlso referred to as shear loading coefficient; Also by for volume response coefficient;
3) experimental data that rubbery sample relaxation modulus changes over time is extracted
According to DMA dynamic thermomechanical to be applied analysis preparation and the same ingredient of studied rubber rubbery sample, in conjunction with reality The setting of border operating condition progress experiment condition, comprising: experimental temperature, test frequency scanning range, plastic strain amplitude;Selection DMA temperature/ Frequency sweep experiment obtains the loss modulus and storage modulus versus time curve of rubber, by DMA tester software, The relaxation modulus data versus time curve of rubber is obtained, which has fully demonstrated rubber in frequency and temperature change item Viscous-elastic behaviour under part;
4) viscoelastic parameter identification is carried out in finite element software
The data that the rubber relaxation modulus of extraction is changed over time are input in the prony curve fitting of ANSYS, Adjusting coeffvalue initial parameter is 0,1, is fitted, obtains relative modulus and slack time i.e. viscoelastic parameter in formula 2,3;
5) it is based on finite element ANSYS, emulation project is subjected to grid dividing, boundary condition is established, sets solver parameter, Carry out dynamical property analysis.
3. a kind of side for establishing constitutive model using finite element software identification rubber viscoelastic parameter according to claim 1 Method, it is characterised in that: step 1: the method that experimental data combination ANSYS carries out viscoelastic parameter identification obtains rubber relaxation modulus Change with time data;
Step 2: ANSYS fitting experimental data obtains viscoelastic parameter.
4. a kind of side for establishing constitutive model using finite element software identification rubber viscoelastic parameter according to claim 1 Method, it is characterised in that: obtain rubber relaxation modulus and change with time data:
Step 1: rubbery sample is formulated according to selected DMA experiment instrument model
Step 2: the actual working environment according to applied by rubber, the condition and range of choice experiment, including temperature, frequency, strain Amplitude;
Step 3: selection DMA temperature frequency scanning experiment obtains rubber loss modulus, storage modulus, fissipation factor;
Step 4: choosing reference temperature appropriate, analyzes software by DMA and show that the relaxation modulus of rubber changes with time number According to.
5. a kind of side for establishing constitutive model using finite element software identification rubber viscoelastic parameter according to claim 1 Method, it is characterised in that: ANSYS fitting experimental data obtains viscoelastic parameter,
Step 1: the rubber relaxation modulus data of extraction are deposited into notepad according to rule, name should be ordered using English Name;
Step 2: entering ANSYS PRONY CURVE FITTING module, and the data that name is completed import, and adjusts initial ginseng Number carries out data fitting;
Step 3: check fitting as a result, initial parameter is adjusted again, until the error of fitting data and experimental data is allowing Range;
Step 4: saving and applies, and obtains rubber viscoelastic parameter.
CN201811263229.9A 2018-10-28 2018-10-28 Method for identifying viscoelastic parameters of rubber by using finite element software to establish constitutive model Active CN109408969B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811263229.9A CN109408969B (en) 2018-10-28 2018-10-28 Method for identifying viscoelastic parameters of rubber by using finite element software to establish constitutive model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811263229.9A CN109408969B (en) 2018-10-28 2018-10-28 Method for identifying viscoelastic parameters of rubber by using finite element software to establish constitutive model

Publications (2)

Publication Number Publication Date
CN109408969A true CN109408969A (en) 2019-03-01
CN109408969B CN109408969B (en) 2022-12-23

Family

ID=65469281

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811263229.9A Active CN109408969B (en) 2018-10-28 2018-10-28 Method for identifying viscoelastic parameters of rubber by using finite element software to establish constitutive model

Country Status (1)

Country Link
CN (1) CN109408969B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110008552A (en) * 2019-03-26 2019-07-12 北京工业大学 Consider the reed-type space deployable structure rapid modeling analysis and optimization method of material viscoelasticity
CN110006644A (en) * 2019-04-20 2019-07-12 北京工业大学 A method of judging steel-rubber roll structure rubber viscoplasticity influence degree under dynamic rotary operating condition
CN110472307A (en) * 2019-07-29 2019-11-19 惠州市德赛西威汽车电子股份有限公司 A kind of display system optical cement removing emulation mode
CN111951970A (en) * 2020-08-07 2020-11-17 山东大学 Dentin material biological tissue mechanical property calibration method, system and device
CN112182928A (en) * 2020-09-17 2021-01-05 中国海洋大学 Software platform for mechanical property analysis and optimization design of rubber products
CN112199823A (en) * 2020-09-17 2021-01-08 中国海洋大学 Numerical prediction method for stress relaxation and damage effects of rubber material
CN112199873A (en) * 2020-10-09 2021-01-08 青岛科技大学 Rubber dynamic heat generation calculation method and device
CN113239476A (en) * 2021-04-01 2021-08-10 中国第一汽车股份有限公司 Rubber bushing constitutive parameter reverse calibration method based on binomial Mooney-Rivlin model
CN113553706A (en) * 2021-07-21 2021-10-26 青岛博锐智远减振科技有限公司 Method for judging inflection point of conical rubber spring bonding experiment
CN112199873B (en) * 2020-10-09 2024-06-11 青岛科技大学 Dynamic heat generation calculation method and device for rubber

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184332A (en) * 2011-05-17 2011-09-14 中国船舶重工集团公司第七二五研究所 Method for expanding dynamic mechanical test data of viscous-elastic material out of test frequency
CN105069241A (en) * 2015-08-19 2015-11-18 山东大学 Step-by-step analysis and prediction method for dynamic performances of rubber material structure
CN106202683A (en) * 2016-07-04 2016-12-07 南京理工大学 A kind of elastomeric material considering correlation effect glues the modeling method of superlastic constitutive model
JP2017224098A (en) * 2016-06-14 2017-12-21 株式会社ブリヂストン Simulation method and simulation program

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184332A (en) * 2011-05-17 2011-09-14 中国船舶重工集团公司第七二五研究所 Method for expanding dynamic mechanical test data of viscous-elastic material out of test frequency
CN105069241A (en) * 2015-08-19 2015-11-18 山东大学 Step-by-step analysis and prediction method for dynamic performances of rubber material structure
JP2017224098A (en) * 2016-06-14 2017-12-21 株式会社ブリヂストン Simulation method and simulation program
CN106202683A (en) * 2016-07-04 2016-12-07 南京理工大学 A kind of elastomeric material considering correlation effect glues the modeling method of superlastic constitutive model

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于海富 等: "有限变形下橡胶材料非线性高弹-粘弹性本构模型", 《橡胶工业》 *
张银喜等: "Abaqus中橡胶粘弹性能参数识别", 《特种橡胶制品》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110008552B (en) * 2019-03-26 2023-06-02 北京工业大学 Reed type space expandable structure rapid modeling analysis and optimization method considering material viscoelasticity
CN110008552A (en) * 2019-03-26 2019-07-12 北京工业大学 Consider the reed-type space deployable structure rapid modeling analysis and optimization method of material viscoelasticity
CN110006644B (en) * 2019-04-20 2021-05-25 北京工业大学 Method for judging rubber viscoelasticity influence degree of steel-rubber roller structure under dynamic rotation working condition
CN110006644A (en) * 2019-04-20 2019-07-12 北京工业大学 A method of judging steel-rubber roll structure rubber viscoplasticity influence degree under dynamic rotary operating condition
CN110472307A (en) * 2019-07-29 2019-11-19 惠州市德赛西威汽车电子股份有限公司 A kind of display system optical cement removing emulation mode
CN111951970A (en) * 2020-08-07 2020-11-17 山东大学 Dentin material biological tissue mechanical property calibration method, system and device
CN112182928B (en) * 2020-09-17 2023-04-07 中国海洋大学 Software system for mechanical property analysis and optimization design of rubber products
CN112199823B (en) * 2020-09-17 2023-04-07 中国海洋大学 Numerical prediction method for stress relaxation and damage effects of rubber material
CN112199823A (en) * 2020-09-17 2021-01-08 中国海洋大学 Numerical prediction method for stress relaxation and damage effects of rubber material
CN112182928A (en) * 2020-09-17 2021-01-05 中国海洋大学 Software platform for mechanical property analysis and optimization design of rubber products
CN112199873A (en) * 2020-10-09 2021-01-08 青岛科技大学 Rubber dynamic heat generation calculation method and device
CN112199873B (en) * 2020-10-09 2024-06-11 青岛科技大学 Dynamic heat generation calculation method and device for rubber
CN113239476A (en) * 2021-04-01 2021-08-10 中国第一汽车股份有限公司 Rubber bushing constitutive parameter reverse calibration method based on binomial Mooney-Rivlin model
CN113239476B (en) * 2021-04-01 2022-07-19 中国第一汽车股份有限公司 Rubber bushing constitutive parameter reverse calibration method based on binomial Mooney-Rivlin model
CN113553706A (en) * 2021-07-21 2021-10-26 青岛博锐智远减振科技有限公司 Method for judging inflection point of conical rubber spring bonding experiment
CN113553706B (en) * 2021-07-21 2024-01-19 青岛博锐智远减振科技有限公司 Method for judging inflection point of conical rubber spring adhesion experiment

Also Published As

Publication number Publication date
CN109408969B (en) 2022-12-23

Similar Documents

Publication Publication Date Title
CN109408969A (en) A method of constitutive model is established using finite element software identification rubber viscoelastic parameter
Carroll A strain energy function for vulcanized rubbers
CN105372136B (en) A kind of fatigue limit method for quick predicting based on strain increment
Manes et al. Experimental and numerical investigations of low velocity impact on sandwich panels
CN103886125B (en) A kind of titanium alloy hot combined shaping method for numerical simulation
CN109470553B (en) Method for realizing rapid conversion of material parameters of asphalt mixture by applying Laplace transformation
CN105653754A (en) Solid rocket engine virtual test base system and testing method
CN107885930A (en) Relay elastic metallic material accelerated storage degradation failure mechanism changes method of discrimination
CN105677973A (en) Test spectrum clipping method achieving monaxial and triaxial random vibration stress equivalence
Needleman et al. Indentation of elastically soft and plastically compressible solids
Suman et al. Development of a multiaxial fatigue damage parameter and life prediction methodology for non-proportional loading
Bentil et al. Viscoelastic properties of inert solid rocket propellants exposed to a shock wave
Dalrymple et al. Elastomer rate-dependence: A testing and material modeling methodology
CN113553739B (en) Mixed explosive explosion output characteristic calculation method
CN107782616A (en) Indentation creep stress exponent method of testing and system
Bayat et al. XFEM analysis of cracked orthotropic media under non-classic thermal shock
CN110750926A (en) Particle swarm algorithm-based high-speed tensile curve processing and predicting method
Zhou Experimental and numerical investigation of soft impact loading on aircraft materials
Su et al. Dynamic mechanical response and a constitutive model of Fe-based high temperature alloy at high temperatures and strain rates
Dou et al. Experimental calibration of ISV damage model constants for pure copper for high-speed impact simulation
Koprowski-Theiß et al. Compressible rubber materials: experiments and simulations
Kim et al. Determination of the viscoelastic properties of agar/agar-gelatin gels based on finite element method optimization
Xu et al. A Rate‐Dependent Constitutive Model of HTPB Propellant Based on Parallel Rheological Framework
Kindrachuck et al. Simulation of non-isothermal mechanical tests on a single crystal nickel-basis superalloy
CN109829222A (en) A kind of composite propellant Compressive Mechanical this structure method of the wide temperature range of width strain rate

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
GR01 Patent grant
GR01 Patent grant