CN110348055A - Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method - Google Patents
Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method Download PDFInfo
- Publication number
- CN110348055A CN110348055A CN201910500501.9A CN201910500501A CN110348055A CN 110348055 A CN110348055 A CN 110348055A CN 201910500501 A CN201910500501 A CN 201910500501A CN 110348055 A CN110348055 A CN 110348055A
- Authority
- CN
- China
- Prior art keywords
- stress
- indicate
- parameter
- chaboche
- equation
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Biomedical Technology (AREA)
- Genetics & Genomics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Physiology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention provides a kind of acquisition of Chaboche Viscoplastic Constitutive Model material parameter and optimization methods, comprising the following steps: step 1) establishes the Chaboche Viscoplastic Constitutive Model of material；The stress-strain data that step 2 is obtained according to experiment, with Ramberg-Osgood equation, the material parameter initial value of quick obtaining Chaboche constitutive model；Step 3) utilizes the material parameter of genetic algorithm optimization Chaboche model.The present invention is imported genetic algorithm, is obtained Optimal Parameters, obtain preferable effect of optimization by verification experimental verification by Ramberg-Osgood equation based on the material parameter of self-compiling program quick obtaining Chaboche constitutive model.
Description
Technical field
The invention belongs to aerospace high-temperature material intensity technique fields, and in particular to a kind of this structure of Chaboche visco-plasticity
Cast material parameter obtains and optimization method.
Background technique
Constitutive model for the analysis of high temperature laser intensity can not only consider the influence of temperature, and can reasonably describe to follow
Various inelastic deformations under ring load, such as cyclic hardening, cyclic softening, nonlinear motion hardening, creep, relaxation and recovery effect
It should wait.Chaboche Viscoplastic Constitutive Model is it is verified that it can be considered that the influence of temperature and time, reasonably description recycles
A variety of materials inelastic deformation under load effect.Chaboche model is applied in actual analysis, it is necessary to have corresponding
Finite element method and program, and reasonably to determine model parameter.And contain a large amount of undetermined parameters in the model, by
Some more complicated material behavior responses can be described in the mathematic(al) representation of these nonlinearities, to parameter
Integral algorithm has certain sensibility.Therefore determination method for parameter is worth research in Chaboche Viscoplastic Constitutive Model.
Summary of the invention
Goal of the invention: Chaboche Viscoplastic Constitutive Model can be used for predicting material stress-strain-responsive, precondition
It is known to model parameter.The present invention proposes that a kind of Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method, obtains
To more accurate model parameter, the estimation for subsequent fatigue life provides foundation.
Technical solution: to achieve the goals above, the invention adopts the following technical scheme:
A kind of Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method, comprising the following steps:
Step 1) establishes the Chaboche Viscoplastic Constitutive Model of material；
Step 2) is based on Ramberg-Osgood equation, obtains Chaboche model initial parameter；
Step 3) utilizes genetic algorithm optimization Chaboche model parameter.
Specific step is as follows in the step 1):
Step 1.1), Chaboche Viscoplastic Constitutive Model is by flow rule, the equation of motion and evolution equation of internal variables three
Part forms；
Flow rule are as follows:
Flow rule is related to Von Mises yield function, the expression formula of yield function F are as follows:
F=F (σij,Vk)=J2(σij-Xij)-k0-R
In formula: VkIndicate internal variable tensor, XijIndicate that kinematic hardening internal variable tensor, R indicate isotropism variable, VkPacket
Include kinematic hardening internal variable tensor XijWith isotropism variable R, k0Indicate initial yield stress, σijIndicate stress tensor, J2Table
It is shown with the second invariant of efficacy deviator, indicates that the equivalents of tensor in bracket, size indicate in physical significance are as follows:
In formula, σ 'ijWith X 'ijRespectively indicate stress tensor σijWith kinematic hardening internal variable tensor XijDeviator form；
Inelastic strain rate equation:
Wherein:Indicate that inelastic strain rate tensor, Λ are expressed as the function of yield function；
Yield function is brought into, inelastic strain rate equation can be obtained:
The equation of motion are as follows:
Λ is expressed as to a power function of yield function:
In formula: K and n is the parameter for describing the power function of the equation of motion, only related with temperature level, is carried it into non-resilient
Strain-rate equation:
Evolution equation of internal variables are as follows:
Kinematic hardening internal variable is indicated using Armstrong-Frederick equation are as follows:
In formula:Indicate kinematic hardening internal variable tensor evolution rate, c, a indicate that material constant, c indicate control hardening
The inclined degree of variable, a indicate the maximum value of control hardening variable,Indicate accumulated plastic strain rate；
Isotropic hardening equation:
In formula,Indicate isotropic hardening internal variable evolution rate, b indicates the speed of stress decrease, and Q is indicated under stress
The amplitude of drop；
Triaxiality strain problems are reduced to one-dimensional stress strain on the basis of multiaxis constitutive model by step 1.2)
Problem is depicted as Chaboche Visco-plastic Model a series of following sides using the Temperature measurement of two component statement materials
Journey:
In formula,Indicate inelastic strain rate, σ indicates that stress, X indicate kinematic hardening parameter, and X=X1+X2, X1It indicates
First kinematic hardening parameter, X2Indicate second kinematic hardening parameter；Indicate the evolution rate of kinematic hardening parameter,WithRespectively indicate the evolution rate of first with second kinematic hardening parameter；Wherein, a1With c1Respectively indicate X1Corresponding a with
C, a2With c2Respectively indicate X2Corresponding a and c；
For the mechanical characteristic of different materials, the hardening component number needed in model is different, under normal circumstances, only needs two
A component can describe the Temperature measurement of material；Therefore the parameter being fitted is needed to have K, n, k0、a1、c1、a2、c2,b,Q；Wherein, K,
N is the parameter for describing the power function of the equation of motion, and for isotropic material, there are also elastic modulus Es and Poisson's ratio μ；It is single
In the case of axis, Poisson's ratio μ without the concern for；When obtaining initial parameter values, need establish it is assumed hereinafter that:
(1) isotropism scalar R describes cyclic hardening, and monotonic tension curve is put aside；
(2) kinematic hardening parameter X1And X2The stress-strain diagram within the scope of short distance and long-range, i.e. X are expressed respectively1Influence modeling
Property strain less than 0.1% when stress intensity, X2Always stress intensity is influenced；
Specifically: short distance is consistent with following size plastic strain with long-range range, is with some plastic strain
Value is divided into two parts as boundary, by nonelastic range, but this cut off value value is different to different materials, in this patent with
For GH4169, plastic strain takes 0.1%；Short distance and long-range represent two regions of different change rates in stress-strain relation,
It is divided with different plastic strain values, can be taken as 0.1% in this patent.
Obtain the relationship of stress σ and other parameters:
σ=X1+X2+k0+R+σv
In formula: σvIndicate viscous stress, andStress remaining value i.e. in above formula.
Specific step is as follows in the step 2):
Step 2.1) reads material stress-strain data, carries out data fitting, bullet using Ramberg-Osgood equation
Property modulus E and initial yield stress k0It is obtained by the ess-strain linear relationship of 1/4 circulation before cylic stress-strain curve；
Step 2.2) obtains the relationship of plastic strain and stress, is successively fitted c according to evolution equation of internal variables2And a2、c1
And a1Value, is finally fitted K, n；
When plastic strain is larger, and in the case where not considering cyclic hardening (softening), due to viscous stress item and plasticity
Strain unrelated, X1It is considered a constant under big plastic strain state, utilizes slope (this of stress when big plastic strain state
Text takes when plastic strain reaches 0.1%), at this time stress variation only and X2Correlation, by the relational expression of stress and each parameter to modeling
Property strain derivation after, and take logarithm, obtain:
A is obtained by curve matching2、c2；In formula, σ indicates stress, εpIndicate inelastic strain.
In small plastic strain, i.e., when plastic strain is less than 0.1%, if σX=σ-X2, represent stress and remove X2Part
Obtain the stress variation situation under small plastic strain state:
A is obtained by curve matching1、c1；
In uncertain time and the relationship of plastic strain, the value that following methods determine K and n can use: straining greatly
When situation, elastic strain increment is almost nil, and strain is completely converted into plastic strain, and stress no longer increases at this time, and K value indicates surplus
Residue stress maximum value, n indicate σvSituation of change before reaching K, by being fitted obtained a1、c1、a2、c2It can true defining K value
Size；Since rate relevant parameter is relatively difficult to measure, and it is related to load duration, therefore assume first that a range, finally
Optimal solution is obtained using parameter global optimization.Plastic strain rate is always less than total strain rate in positive loading procedure, can be in the hope of
The range of n value out；
It is the material parameter for describing the function of yield function, the national defence in 2013 write by Yang Xiaoguang, Shi Duoqi about K, n
Have recorded in the books of " visco-plasticity constitutive theoryr and its application " that industrial publishing house publishes, K, n do not have practical significance, be for
Fit inelastic strain rate and stress, internal variable relationship and the power function form that constructs, the books page 22 has specific note
Carry content.
Step 2.3), the expression formula of isotropism variable R are as follows:
R=Q (1-e-bp)
By to the maximum stress difference σ recycled every timeRThe size of Q and b are estimated with accumulated plastic strain p, in which:
Specific step is as follows in the step 3):
Step 3.1) determines the thresholding of parameter to be optimized on the basis of step 2) obtained initial parameter；
The stress-strain data that experiment obtains are imported into genetic algorithm program, Q by step 3.2)cFor Chaboche mould
The theoretical stress value σ that type calculatesthWith actual stress value σexplAbsolute value of the difference, the objective function ∑ Q of selectioncFor Chaboche
The theoretical stress value σ that model calculatesthWith actual stress value σexplThe sum of absolute value of the difference, expression is as follows:
∑Qc=∑ | σexpl-σth|
Step 3.3) starts to optimize, the Chaboche model parameter after obtaining genetic algorithm optimization.
Step 3.1) determines the thresholding of parameter to be optimized on the basis of step 2) obtained initial parameter；
The stress-strain data that experiment obtains are imported into genetic algorithm program, Q by step 3.2)cFor Chaboche mould
The theoretical stress value σ that type calculatesthWith actual stress value σexplAbsolute value of the difference, the objective function ∑ Q of selectioncFor Chaboche
The theoretical stress value σ that model calculatesthWith actual stress value σexplThe sum of absolute value of the difference, expression is as follows:
∑Qc=∑ | σexpl-σth|
Step 3.3) starts to optimize, the Chaboche model parameter after obtaining genetic algorithm optimization.
It is obtained and optimization method the utility model has the advantages that the present invention provides a kind of Chaboche Viscoplastic Constitutive Model material parameter,
Compared with the prior art by using the above technical solution, it has following technical effect that
(1) it object of the present invention is to obtain Chaboche Viscoplastic Constitutive Model material parameter and optimize, obtains more
Accurate model parameter.
(2) present invention uses Ramberg-Osgood equation, and quick obtaining Chaboche model initial parameter utilizes heredity
Algorithm optimization Chaboche model parameter, improves the accuracy of model.
Detailed description of the invention
Fig. 1 is genetic algorithm flow chart.
Fig. 2 is flow chart of the method for the present invention.
Fig. 3 is the source of parameter in the present invention.
Specific embodiment
Below with reference to embodiment, the present invention will be described in detail.
The Chaboche Viscoplastic Constitutive Model material parameter that the present invention obtains is by its reliability of verification experimental verification and accurately
Property.
A kind of Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method, comprising the following steps:
Step 1) establishes the Chaboche Viscoplastic Constitutive Model of material；Specific step is as follows in step 1):
Step 1.1), Chaboche Viscoplastic Constitutive Model is by flow rule, the equation of motion and evolution equation of internal variables three
Part forms；
Flow rule are as follows:
Flow rule is related to Von Mises yield function, the expression formula of yield function F are as follows:
F=F (σij,Vk)=J2(σij-Xij)-k0-R
In formula: VkIndicate internal variable tensor, XijIndicate that kinematic hardening internal variable tensor, R indicate isotropism variable, VkPacket
Include kinematic hardening internal variable tensor XijWith isotropism variable R, k0Indicate initial yield stress, σijIndicate stress tensor, J2Table
It is shown with the second invariant of efficacy deviator, indicates that the equivalents of tensor in bracket, size indicate in physical significance are as follows:
In formula, σ 'ijWith X 'ijRespectively indicate stress tensor σijWith kinematic hardening internal variable tensor XijDeviator form；
Inelastic strain rate equation:
Wherein:Indicate that inelastic strain rate tensor, Λ are expressed as the function of yield function；
Yield function is brought into, inelastic strain rate equation can be obtained:
The equation of motion are as follows:
Λ is expressed as to a power function of yield function:
In formula: K and n is the parameter for describing the power function of the equation of motion, only related with temperature level, is carried it into non-resilient
Strain-rate equation:
Evolution equation of internal variables are as follows:
Kinematic hardening internal variable is indicated using Armstrong-Frederick equation are as follows:
In formula:Indicate kinematic hardening internal variable tensor evolution rate, c, a indicate that material constant, c indicate control hardening
The inclined degree of variable, a indicate the maximum value of control hardening variable,Indicate accumulated plastic strain rate；
Isotropic hardening equation:
In formula,Indicate isotropic hardening internal variable evolution rate, b indicates the speed of stress decrease, and Q is indicated under stress
The amplitude of drop；
Triaxiality strain problems are reduced to one-dimensional stress strain on the basis of multiaxis constitutive model by step 1.2)
Problem is depicted as Chaboche Visco-plastic Model a series of following sides using the Temperature measurement of two component statement materials
Journey:
In formula,Indicate inelastic strain rate, σ indicates that stress, X indicate kinematic hardening parameter, and X=X1+X2, X1It indicates
First kinematic hardening parameter, X2Indicate second kinematic hardening parameter；Indicate the evolution rate of kinematic hardening parameter,WithRespectively indicate the evolution rate of first with second kinematic hardening parameter；Wherein, a1With c1Respectively indicate X1Corresponding a with
C, a2With c2Respectively indicate X2Corresponding a and c；
For the mechanical characteristic of different materials, the hardening component number needed in model is different, under normal circumstances, only needs two
A component can describe the Temperature measurement of material；Therefore the parameter being fitted is needed to have K, n, k0、a1、c1、a2、c2,b,Q；Wherein, K,
N is the parameter for describing the power function of the equation of motion, and for isotropic material, there are also elastic modulus Es and Poisson's ratio μ；It is single
In the case of axis, Poisson's ratio μ without the concern for；When obtaining initial parameter values, need establish it is assumed hereinafter that:
(1) isotropism scalar R describes cyclic hardening, and monotonic tension curve is put aside；
(2) kinematic hardening parameter X1And X2The stress-strain diagram within the scope of short distance and long-range, i.e. X are expressed respectively1Influence modeling
Property strain less than 0.1% when stress intensity, X2Always stress intensity is influenced；
Specifically: short distance is consistent with following size plastic strain with long-range range, is with some plastic strain
Value is divided into two parts as boundary, by nonelastic range, but this cut off value value is different to different materials, in this patent with
For GH4169, plastic strain takes 0.1%；Short distance and long-range represent two regions of different change rates in stress-strain relation,
It is divided with different plastic strain values, can be taken as 0.1% in this patent.
Obtain the relationship of stress σ and other parameters:
σ=X1+X2+k0+R+σv
In formula: σvIndicate viscous stress, andStress remaining value i.e. in above formula.
Step 2) is based on Ramberg-Osgood equation, obtains Chaboche model initial parameter；It is specific in step 2)
Steps are as follows:
Step 2.1) reads material stress-strain data, carries out data fitting, bullet using Ramberg-Osgood equation
Property modulus E and initial yield stress k0It is obtained by the ess-strain linear relationship of 1/4 circulation before cylic stress-strain curve；
Step 2.2) obtains the relationship of plastic strain and stress according to evolution equation of internal variables.
When plastic strain is larger, and in the case where not considering cyclic hardening (softening), due to viscous stress item and plasticity
Strain unrelated, X1It is considered a constant under big plastic strain state, utilizes slope (this of stress when big plastic strain state
Text takes when plastic strain reaches 0.1%), at this time stress variation only and X2Correlation, by the relational expression of stress and each parameter to modeling
Property strain derivation after, and take logarithm, obtain:
A is obtained by curve matching2、c2；In formula, σ indicates stress, εpIndicate inelastic strain.
In small plastic strain, i.e., when plastic strain is less than 0.1%, if σX=σ-X2, represent stress and remove X2Part
Obtain the stress variation situation under small plastic strain state:
A is obtained by curve matching1、c1；
In uncertain time and the relationship of plastic strain, the value that following methods determine K and n can use: straining greatly
When situation, elastic strain increment is almost nil, and strain is completely converted into plastic strain, and stress no longer increases at this time, and K value indicates surplus
Residue stress maximum value, n indicate σvSituation of change before reaching K, by being fitted obtained a1、c1、a2、c2It can true defining K value
Size；Since rate relevant parameter is relatively difficult to measure, and it is related to load duration, therefore assume first that a range, finally
Optimal solution is obtained using parameter global optimization.Plastic strain rate is always less than total strain rate in positive loading procedure, can be in the hope of
The range of n value out；
It is the material parameter for describing the function of yield function, the national defence in 2013 write by Yang Xiaoguang, Shi Duoqi about K, n
Have recorded in the books of " visco-plasticity constitutive theoryr and its application " that industrial publishing house publishes, K, n do not have practical significance, be for
Fit inelastic strain rate and stress, internal variable relationship and the power function form that constructs, the books page 22 has specific note
Content is carried, as shown in Figure 3.
Step 2.3), the expression formula of isotropism variable R are as follows:
R=Q (1-e-bp)
By to the maximum stress difference σ recycled every timeRThe size of Q and b are estimated with accumulated plastic strain p, in which:
Step 3) starts to optimize, the Chaboche model parameter after obtaining genetic algorithm optimization, i.e., excellent using genetic algorithm
Change Chaboche model parameter, specific step is as follows in step 3):
Step 3.1) determines the thresholding of parameter to be optimized on the basis of step 2) obtained initial parameter；
The stress-strain data that experiment obtains are imported into genetic algorithm program, genetic algorithm flow chart by step 3.2)
As shown in Figure 1.QcThe theoretical stress value σ calculated for Chaboche modelthWith actual stress value σexplAbsolute value of the difference, selection
Objective function ∑ QcThe theoretical stress value σ calculated for Chaboche modelthWith actual stress value σexplAbsolute value of the difference it
With expression is as follows:
∑Qc=∑ | σexpl-σth|
Step 3.3) starts to optimize, the Chaboche model parameter after obtaining genetic algorithm optimization.
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (4)
1. a kind of Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method, it is characterised in that: including following step
It is rapid:
Step 1) establishes the Chaboche Viscoplastic Constitutive Model of material；
Step 2) is based on Ramberg-Osgood equation, obtains Chaboche model initial parameter；
Step 3) utilizes genetic algorithm optimization Chaboche model parameter.
2. Chaboche Visco-plastic Model material parameter according to claim 1 obtains and optimization method, it is characterised in that:
Specific step is as follows in the step 1):
Step 1.1), Chaboche Viscoplastic Constitutive Model is by flow rule, the equation of motion and evolution equation of internal variables three parts
Composition；
Flow rule are as follows:
Flow rule is related to Von Mises yield function, the expression formula of yield function F are as follows:
F=F (σij,Vk)=J2(σij-Xij)-k0-R
In formula: VkIndicate internal variable tensor, XijIndicate that kinematic hardening internal variable tensor, R indicate isotropism variable, VkIncluding fortune
Dynamic hardening internal variable tensor XijWith isotropism variable R, k0Indicate initial yield stress, σijIndicate stress tensor, J2It indicates
The second invariant of efficacy deviator indicates that the equivalents of tensor in bracket, size indicate in physical significance are as follows:
In formula, σ 'ijWith X 'ijRespectively indicate stress tensor σijWith kinematic hardening internal variable tensor XijDeviator form；
Inelastic strain rate equation:
Wherein:Indicate that inelastic strain rate tensor, Λ are expressed as the function of yield function；
Yield function is brought into, inelastic strain rate equation can be obtained:
The equation of motion are as follows:
Λ is expressed as to a power function of yield function:
In formula: K and n is the parameter for describing the power function of the equation of motion, only related with temperature level, carries it into inelastic strain
Rate equation:
Evolution equation of internal variables are as follows:
Kinematic hardening internal variable is indicated using Armstrong-Frederick equation are as follows:
In formula:Indicate kinematic hardening internal variable tensor evolution rate, c, a indicate that material constant, c indicate control hardening variable
Inclined degree, a indicate control hardening variable maximum value,Indicate accumulated plastic strain rate；
Isotropic hardening equation:
In formula,Indicate isotropic hardening internal variable evolution rate, b indicates the speed of stress decrease, and Q indicates stress decrease
Amplitude；
Triaxiality strain problems are reduced to one-dimensional stress strain problems on the basis of multiaxis constitutive model by step 1.2),
Using the Temperature measurement of two component statement materials, Chaboche Visco-plastic Model is depicted as a series of following equations:
In formula,Indicate inelastic strain rate, σ indicates that stress, X indicate kinematic hardening parameter, and X=X1+X2, X1Indicate first
Kinematic hardening parameter, X2Indicate second kinematic hardening parameter；Indicate the evolution rate of kinematic hardening parameter,WithRespectively
Indicate the evolution rate of first with second kinematic hardening parameter；Wherein, a1With c1Respectively indicate X1Corresponding a and c, a2With c2
Respectively indicate X2Corresponding a and c；
For the mechanical characteristic of different materials, the hardening component number needed in model is different, under normal circumstances, only needs two points
The Temperature measurement of amount description material；Therefore the parameter being fitted is needed to have K, n, k0、a1、c1、a2、c2, b, Q, wherein K, n be description
The parameter of the power function of the equation of motion；For isotropic material, there are also elastic modulus Es and Poisson's ratio μ；Obtaining parameter
When initial value, need establish it is assumed hereinafter that:
(1) isotropism scalar R describes cyclic hardening, and monotonic tension curve is put aside；
(2) kinematic hardening parameter X1And X2The stress-strain diagram within the scope of short distance and long-range, i.e. X are expressed respectively1Plasticity is influenced to answer
The stress intensity to become smaller when 0.1%, X2Always stress intensity is influenced；
Obtain the relationship of stress σ and other parameters:
σ=X1+X2+k0+R+σv
In formula: σvIndicate viscous stress, andStress remaining value i.e. in above formula.
3. Chaboche Visco-plastic Model material parameter according to claim 1 obtains and optimization method, it is characterised in that:
Specific step is as follows in the step 2):
Step 2.1) reads material stress-strain data, carries out data fitting, springform using Ramberg-Osgood equation
Measure E and initial yield stress k0It is obtained by the ess-strain linear relationship of 1/4 circulation before cylic stress-strain curve；
Step 2.2) obtains the relationship of plastic strain and stress, is successively fitted c according to evolution equation of internal variables2And a2、c1And a1
Value, is finally fitted K, n；
Step 2.3), the expression formula of isotropism variable R are as follows:
R=Q (1-e-bp)
By to the maximum stress difference σ recycled every timeRThe size of Q and b are estimated with accumulated plastic strain p, in which:
4. Chaboche Visco-plastic Model material parameter according to claim 1 obtains and optimization method, it is characterised in that:
Specific step is as follows in the step 3):
Step 3.1) determines the thresholding of parameter to be optimized on the basis of step 2) obtained initial parameter；
The stress-strain data that experiment obtains are imported into genetic algorithm program, Q by step 3.2)cFor Chaboche model meter
The theoretical stress value σ of calculationthWith actual stress value σexplAbsolute value of the difference, the objective function ∑ Q of selectioncFor Chaboche model
The theoretical stress value σ of calculatingthWith actual stress value σexplThe sum of absolute value of the difference, expression is as follows:
∑Qc=∑ | σexpl-σth|
Step 3.3) starts to optimize, the Chaboche model parameter after obtaining genetic algorithm optimization.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910500501.9A CN110348055B (en) | 2019-06-11 | 2019-06-11 | Method for obtaining and optimizing material parameters of Chaboche viscoplasticity constitutive model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910500501.9A CN110348055B (en) | 2019-06-11 | 2019-06-11 | Method for obtaining and optimizing material parameters of Chaboche viscoplasticity constitutive model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110348055A true CN110348055A (en) | 2019-10-18 |
| CN110348055B CN110348055B (en) | 2021-01-08 |
Family
ID=68181786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910500501.9A Active CN110348055B (en) | 2019-06-11 | 2019-06-11 | Method for obtaining and optimizing material parameters of Chaboche viscoplasticity constitutive model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110348055B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112651089A (en) * | 2020-08-27 | 2021-04-13 | 同济大学 | Bending shear metal component hybrid reinforcement restoring force model and parameter determination method |
| CN112668167A (en) * | 2020-12-21 | 2021-04-16 | 合图智造科技(西安)有限公司 | Material parameter construction method based on small amount of experimental data |
| CN113742914A (en) * | 2021-09-02 | 2021-12-03 | 南京工业大学 | Method suitable for predicting cyclic load deformation behaviors of multiple control modes |
| CN117558381A (en) * | 2024-01-12 | 2024-02-13 | 四川大学 | Calculation method of plastic hardening model related to temperature and strain rate of metal material |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107704642A (en) * | 2017-07-07 | 2018-02-16 | 水利部交通运输部国家能源局南京水利科学研究院 | A kind of soil constitutive model parameter determination method based on intelligent algorithm |
| CN108931448A (en) * | 2018-05-07 | 2018-12-04 | 华南理工大学 | A kind of prediction technique of high chrome Material Thermodynamics response and spleen tissue extracts damage |
| CN109614715A (en) * | 2018-12-13 | 2019-04-12 | 电子科技大学 | A kind of lower Field strength method and its application for considering notch effect of multiaxial loading effect |
-
2019
- 2019-06-11 CN CN201910500501.9A patent/CN110348055B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107704642A (en) * | 2017-07-07 | 2018-02-16 | 水利部交通运输部国家能源局南京水利科学研究院 | A kind of soil constitutive model parameter determination method based on intelligent algorithm |
| CN108931448A (en) * | 2018-05-07 | 2018-12-04 | 华南理工大学 | A kind of prediction technique of high chrome Material Thermodynamics response and spleen tissue extracts damage |
| CN109614715A (en) * | 2018-12-13 | 2019-04-12 | 电子科技大学 | A kind of lower Field strength method and its application for considering notch effect of multiaxial loading effect |
Non-Patent Citations (3)
| Title |
|---|
| A.H. MAHMOUDI, ET AL.: "Parameter determination of Chaboche kinematic hardening model using a multi objective Genetic Algorithm", 《COMPUTATIONAL MATERIALS SCIENCE》 * |
| 吴志荣 等: "基于Chaboche模型的金属材料稳态循环应力-应变曲线的本构建模方法", 《机械工程材料》 * |
| 徐可君 等: "基于Chaboche理论的GH901合金本构模型改进", 《航空发动机》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112651089A (en) * | 2020-08-27 | 2021-04-13 | 同济大学 | Bending shear metal component hybrid reinforcement restoring force model and parameter determination method |
| CN112668167A (en) * | 2020-12-21 | 2021-04-16 | 合图智造科技(西安)有限公司 | Material parameter construction method based on small amount of experimental data |
| CN113742914A (en) * | 2021-09-02 | 2021-12-03 | 南京工业大学 | Method suitable for predicting cyclic load deformation behaviors of multiple control modes |
| CN113742914B (en) * | 2021-09-02 | 2023-08-08 | 南京工业大学 | Method suitable for predicting cyclic load deformation behaviors of multiple control modes |
| CN117558381A (en) * | 2024-01-12 | 2024-02-13 | 四川大学 | Calculation method of plastic hardening model related to temperature and strain rate of metal material |
| CN117558381B (en) * | 2024-01-12 | 2024-03-22 | 四川大学 | Calculation method of plastic hardening model related to temperature and strain rate of metal material |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110348055B (en) | 2021-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110334405A (en) | High temperature Multiaxial Low Cycle Fatigue Life Prediction method based on this structure of Chaboche and Lemaitre damage model | |
| CN110348055A (en) | Chaboche Viscoplastic Constitutive Model material parameter obtains and optimization method | |
| CN108333208A (en) | A kind of complete machine grade product storage-life accelerated test method | |
| CN101436057A (en) | Numerical control machining tool heat error Bayes network compensation method | |
| CN101706355A (en) | NExT/ARMA-based structural response analysis method | |
| Cortez-Vega et al. | A hybrid dynamic model of shape memory alloy spring actuators | |
| Makris et al. | The engineering merit of the “effective period” of bilinear isolation systems | |
| CN103853927B (en) | Based on the method that cluster global optimization approach predicts material behavior | |
| CN106442727B (en) | A kind of method and system recognizing hard coat material mechanics parameters | |
| Chen et al. | Training material models using gradient descent algorithms | |
| Yassin et al. | Extended analysis of bpso structure selection of nonlinear auto-regressive model with exogenous inputs (NARX) of direct current motor. | |
| Babaei et al. | Adaptive Neuro-Fuzzy Inference System (ANFIS) integrated with genetic algorithm to optimize piezoelectric cantilever-oscillator-spring energy harvester: Verification with closed-form solution | |
| Haddadi et al. | Determination of Material Parameters of In740H Under Different Experimental Situations Using Chaboche Model | |
| Dey et al. | Univariate extreme value analysis | |
| Hasanpour et al. | Finite element simulation of polymer behaviour using a three-dimensional, finite deformation constitutive model | |
| CN103488826A (en) | Experience acceleration model based degradation amount distribution parameter modeling and extrapolating method | |
| CN111721647A (en) | Low-cycle fatigue test data processing and internal stress evaluation method | |
| CN104516858A (en) | Phase diagram matrix method for nonlinear dynamic behavior analysis | |
| Zhan et al. | High‐efficiency algorithms for simulating metal failure effects under multiaxial repeated loadings | |
| Loukanov et al. | Multi-criteria identification of VibroBot dynamic characteristics | |
| CN114251075A (en) | Reservoir exploitation scheme determination method, device and equipment based on multi-target parameters | |
| CN117289337B (en) | Rayleigh surface wave inversion method and device, electronic equipment and storage medium | |
| Wang et al. | Operator based robust nonlinear control design to an ionic polymer metal composite with uncertainties and input constraints | |
| Politis et al. | A method of determining unified viscoplastic constitutive equations for hot forging simulations | |
| Tsalouchidis et al. | IDA-based seismic collapse patterns and their predictability by generalized linear models |
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 |