CN103605830B  Establishing method of constructive model for rolling simulation of GH4169 materials  Google Patents
Establishing method of constructive model for rolling simulation of GH4169 materials Download PDFInfo
 Publication number
 CN103605830B CN103605830B CN201310499514.1A CN201310499514A CN103605830B CN 103605830 B CN103605830 B CN 103605830B CN 201310499514 A CN201310499514 A CN 201310499514A CN 103605830 B CN103605830 B CN 103605830B
 Authority
 CN
 China
 Prior art keywords
 stress
 hardening
 strain
 alloy
 parameter
 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.)
 Active
Links
 239000000463 material Substances 0.000 title claims abstract description 84
 238000004088 simulation Methods 0.000 title claims abstract description 17
 238000005096 rolling process Methods 0.000 title claims abstract description 15
 239000000956 alloy Substances 0.000 claims abstract description 32
 229910045601 alloy Inorganic materials 0.000 claims abstract description 32
 REDXJYDRNCIFBQUHFFFAOYSAN aluminium(3+) Chemical class data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5BPC90ZXh0Pgo8dGV4dCB4PScxNjUuNicgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+bDwvdGV4dD4KPHRleHQgeD0nMTc0Ljc4NicgeT0nMTU0JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+MzwvdGV4dD4KPHRleHQgeD0nMTg4LjY1JyB5PScxNTQnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID4rPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMzUuMDQ1NScgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPkE8L3RleHQ+Cjx0ZXh0IHg9JzUxLjA0MDknIHk9JzUzLjU5MDknIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjNweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5sPC90ZXh0Pgo8dGV4dCB4PSc1Ni4zNjQ3JyB5PSc0NC4zMTgyJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjE1cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+MzwvdGV4dD4KPHRleHQgeD0nNjQuMzk5NicgeT0nNDQuMzE4MicgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToxNXB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPis8L3RleHQ+Cjwvc3ZnPgo= [Al+3] REDXJYDRNCIFBQUHFFFAOYSAN 0.000 claims abstract description 32
 238000000034 method Methods 0.000 claims abstract description 7
 238000004458 analytical method Methods 0.000 claims abstract description 6
 238000005242 forging Methods 0.000 claims abstract description 4
 238000000465 moulding Methods 0.000 claims description 3
 230000003068 static Effects 0.000 claims description 3
 238000007906 compression Methods 0.000 claims description 2
 241000208340 Araliaceae Species 0.000 claims 1
 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
 238000005553 drilling Methods 0.000 claims 1
 235000005035 ginseng Nutrition 0.000 claims 1
 235000008434 ginseng Nutrition 0.000 claims 1
 239000002184 metal Substances 0.000 abstract description 2
 238000005457 optimization Methods 0.000 abstract 1
 238000009864 tensile test Methods 0.000 abstract 1
 206010039897 Sedation Diseases 0.000 description 1
 238000004364 calculation method Methods 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
Abstract
The invention discloses an establishing method of a constructive model for rolling simulation of Gh4169 materials. The method comprises the following steps of 1) GH4169 alloy monotonic tensile test; 2) GH4169 alloy lowcycle strain control symmetric loading test; 3) constructive model parameter identification including selecting a YoshidaUemori model and obtaining constructive model parameters through software computation; 4) model verification including after finishing secondary development of a material constructive model, verifying the material constructive model through a cuboid forging process. The establishing method of the constructive model for rolling simulation of the GH4169 materials has the advantages that the constructive model established through the method can accurately represent the stressstrain of deformed GH 4169 materials and really reflect the flow stress change of metal. The method can provide accurate material models for numerical simulation software, improve the accuracy of simulation results and have significance to material deformation characteristic analysis and forming process optimization.
Description
Technical field
The present invention relates to aviation processing technique field, when specifically belonging to roll forming, material constitutive relation model builds
Cube method.
Background technology
Constitutive relation, the relation of stress and strain when that is, material deforms.The constitutive model of material is in plastic history
Indispensable basic theory model, the computational accuracy of model directly affects the result of calculation of numerical simulation.Most limited
Meta software provides flow stress variation model in the range of uniform temperature, strain and strain rate for the material, then passes through interpolation
Calculate flow stress model in the range of other temperature, strain and strain rate for the material with the method for extrapolation.Actually most
Flow stress model in the range of different temperatures, strain and strain rate for the material changes very greatly, in different distortion condition model
Enclose the interior method using interpolation and extrapolation and calculate that flow stress can have a strong impact on the precision of finite element simulation.Meanwhile, at a certain group
Flow stress curve under deformation condition is made up of some experimental data points, and in finite element software, for determining every stream
Several picks of moving curve are counted out very few, and the flow stress curve that so obtains has very big gap with experiment gained, causes emulation
Result precision substantially reduces.
Because the GH4169 alloy flow stress curve in simulation softward material database is the side by interpolation and extrapolation
Method obtains, and very few for determining the data point of a flow stress curve, causes result and the experiment institute of numerical simulation
Obtain result and differ larger, with the increase of strain, the error of flow stress curve is also gradually increased.General Plastic Forming
Journey, deformation ratio is larger, will have a strong impact on the precision of numerical simulation using inaccurate material model.Therefore, set up accurate material
Material model is to ensure that the precondition of finite element simulation precision.Accurately material model can improve the accuracy of analog result,
For analysis of material deformation behaviour, Optimize the forming scheme technique is significant.
Content of the invention
In order to solve abovementioned technical problem, the invention aims to set up accurate GH4169 material rolling simulation using
Constitutive model it is ensured that rolling simulation accuracy.Concrete technical scheme is as follows：
A kind of method for building up of GH4169 material rolling simulation constitutive model, comprises the steps：
1) GH4169 alloy monotonic tension experiment：
Purpose is to obtain GH4169 alloy monotonic tension performance at ambient temperature, is GH4169 material rolling this structure mould
The determination of type material parameter provides data；
2) GH4169 alloy low circulation strain controlling asymmetrical load experiment：
Purpose is the tension and compression circulation experiment obtaining GH4169 alloy at room temperature, in the range of low circulation number of times, symmetric strained
Stress one strain retardant curve, determine the material of isotropic hardening and kinematic hardening in GH4169 material rolling constitutive model
Material parameter；
3) Parameters of constitutive model identification
From Yoshida Uemori model, constitutive relation model parameter is calculated by software；
4) model checking
After completing secondary development to material constitutive model, the forging process using cuboid is verified.
The method of described GH4169 alloy monotonic tension experiment is：
Feeding from GH4169 alloy cold bar, is processed intoWithStandard tensile specimen carry out
Static(al) stretching experiment.
The method of described GH4169 alloy low circulation strain controlling asymmetrical load experiment is：By GH4169 alloy cold bar
Blank is processed intoStandard lowcycle fatigue sample, carry out low circulation strain controlling asymmetrical load experiment.Range of strain
It is respectively 0.4%, 0.8%, 1.2%, 1.6%, 2.0%, the stressstrain curve and the circulation that obtain the first loaded cycle are steady
Stress one strain lag loop after fixed.
Described Parameters of constitutive model knows method for distinguishing：
From Yoshida Uemori model mainly comprise the following steps：Dynamics, yield function and related streams
Dynamic criterion, isotropic hardening and kinematic hardening；The particular content of each several part is respectively：
1. kinetics equation
In formula,For total strain rate,WithIt is respectively elastic strain speed and moulding strain rate；
2. yield function and related flowing criterion
Comprise isotropism one kinematic hardening and the yield function of related flowing criterion is
In formula, S and α is respectively cauchy stress deviator and back stress deviator, Y and R is respectively initial yield stress and each
To same sex hardening stress；Stress be defined as material unit and longpending on active force, cauchy stress is exactly true in analysis on Large Deformation
Stress, is different from nominal stress；And in the case of small deformation, nominal stress and trus stress numerically equal；
Cauchy stress
σ=F/A
Nominal stress σ '=F/Ao
In formula, A is the instantaneous cross section of deformable material；Ao is the sectional area before material deformation；
3. isotropic hardening equation
The isotropic hardening property of material to be represented by internal variable isotropic hardening stress, and isotropic hardening should
The EVOLUTION EQUATION of power is
In formula,For equivalent plastic strain speed, Rsat is the saturation value of R in the case of big plastic strain, m be control each to
The material parameter of same sex hardening rate；Above EVOLUTION EQUATION is integrated, obtains
R=R_{sat}(1e^{mp})
Kinematic hardening
In Yoshida Uemori model, back stress α is made up of two parts：α 1 is nonlinear kinematic hardening back stress, α 2
For linear kinematic hardening back stress；
α=α_{1}+α_{2}
In formula, c, α are and nonlinear servoactuated related material parameter, H '_{∞}For linear kinematic hardening coefficient, α 1, α 2 are respectively
For nonlinear kinematic hardening back stress and linear kinematic hardening back stress；
In Yoshida Uemori model, material parameter undetermined includes：
Y is yield surface radius, i.e. initial yield stress；
R_{sat}Saturation value for isotropic hardening stress；
M is the material parameter controlling isotropic hardening speed；
C is the material parameter controlling nonlinear kinematic hardening speed：
α is the material parameter related to nonlinear kinematic hardening；
H′_{∞}For linear kinematic hardening coefficient；
Using constitutive model as 6 material parameters function f (Y, R_{sat}, m, a, C, H '_{∞}), using Matlab software pair
The stressstrain curve that gained tested by GH4169 alloy carries out curve fitting, you can obtain Y, R_{sat}, m, a, C, H '_{∞}Value, profit
Try to achieve the material parameter of GH4169 alloy state before rolling with the method.
Beneficial effects of the present invention：This patent can accurately reproduce the essstrain after the deformation of GH4169 material, truly
The flow stress change of reflection metal.There is provided accurate material model for numerical simulation software, the essence of analog result can be improved
Exactness, for analysis of material deformation behaviour, Optimize the forming scheme technique is significant.
Brief description
Fig. 1 is GH4169 material stressstrain curve (experimental result and analog result).
Specific embodiment：
1) GH4169 alloy monotonic tension experiment
Draw Zhong Xingneng to obtain GMI69 alloy dullness at ambient temperature, take from GH4169 alloy cold bar
Material, is processed intoWithStandard tensile specimen carry out static(al) stretching experiment.Experiment condition is shown in Table I.
The experiment condition of table 1 stretching experiment
Experiment condition  Stretching experiment one  Stretching experiment two 
Specimen size  Φ8mm  Φ5mm 
Strain rate  4.8mm/min, 3.2mm/min, 1.6mm/min  3mm/min, 2mm/min, 1mm/min 
Temperature  Room temperature  Room temperature 
Range of strain  Stretching is up to sample fracture  Stretching is up to sample fracture 
Experimental facilities  WDS100 electronic type universal experimental machine  AGIS 50kN universal testing machine 
2) GH4169 alloy low circulation strain controlling asymmetrical load experiment
GH4169 alloy cold pieces of bar stock is processed intoStandard lowcycle fatigue sample, carry out low circulation
Strain controlling asymmetrical load is tested.Range of strain is respectively 0.4%, 0.8%, 1.2%, 1.6%, 2.0%, obtains the first loading
Stress one strain lag loop after the stressstrain curve of circulation and stable circulation.
3) Parameters of constitutive model identification
YoshidaUemori (Jitian is upper gloomy) kinematic hardening material model can accurate description strain paths change
When material property change, thus preferably reflection complex loading in the case of material anisotropy.
Yoshida Uemori model mainly comprises the following steps：Dynamics, yield function are accurate to related flowing
Then, isotropic hardening and kinematic hardening.The particular content of each several part is respectively：
1. kinetics equation
In formula,For total strain rate,WithIt is respectively elastic strain speed and moulding strain rate；
2. yield function and related flowing criterion
Comprise isotropism one kinematic hardening and the yield function of related flowing criterion is
In formula, S and a is respectively cauchy stress deviator and back stress deviator, Y and R is respectively initial yield stress and each
To same sex hardening stress；Stress be defined as material unit and longpending on active force, cauchy stress is exactly true in analysis on Large Deformation
Stress, is different from nominal stress；And in the case of small deformation, nominal stress and trus stress numerically equal；
Cauchy stress
σ=F/A
Nominal stress σ '=F/Ao
In formula, A is the instantaneous cross section of deformable material；Ao is the sectional area before material deformation.
3. isotropic hardening equation
The isotropic hardening property of material to be represented by internal variable isotropic hardening stress, isotropic hardening stress
EVOLUTION EQUATION be
In formula,For equivalent plastic strain speed, Rsat is the saturation value of R in the case of big plastic strain, m be control each to
The material parameter of same sex hardening rate.Above EVOLUTION EQUATION is integrated, obtains
R=R_{sat}(1e^{mp})
Kinematic hardening
In Yoshida Uemori model, back stress α is made up of two parts：α_{1}For nonlinear kinematic hardening back stress, α_{2}
For linear kinematic hardening back stress.
α=α_{1}+α_{2}
In formula, c, α are and nonlinear servoactuated related material parameter, H '_{∞}For linear kinematic hardening coefficient, α_{1}、α_{2}It is respectively
Nonlinear kinematic hardening back stress and linear kinematic hardening back stress.
In Yoshida Uemori model, material parameter undetermined includes：
Y is yield surface radius, i.e. initial yield stress；
R_{sat}Saturation value for isotropic hardening stress；
M is the material parameter controlling isotropic hardening speed；
C is the material parameter controlling nonlinear kinematic hardening speed：
α is the material parameter related to nonlinear kinematic hardening；
H′_{∞}For linear kinematic hardening coefficient.
Using constitutive model as 6 material parameters function f (Y, R_{sat}, m, a, C, H '_{∞}), public using U.S. MathWorks
The stressstrain curve that the Matlab software that department produces tests gained to GH4169 alloy carries out curve fitting, you can obtain Y,
R_{sat}, m, a, C, H '_{∞}Value, the GH4169 alloy tried to achieve using the method material parameter of state before rolling is as shown in table 2.
The material parameter of table 2GH4169 alloy
Materials behavior  a/MPa  C  H′_{∞}/MPa  m  R_{sat}/MPa  Y/MPa 
Original blank  300  7.8  200  7.75  200  560 
4) model checking
After completing secondary development to material constitutive model, the forging process using cuboid is verified.Using
GH4169 alloy original blank material parameter carries out numerical simulation, the result of numerical simulation and experimental result is contrasted, sees
Fig. 1, finds that analog result can be good at coincideing with experimental result, demonstrates the reliability of model.
Claims (1)
1. a kind of method for building up of GH4169 material rolling simulation constitutive model is it is characterised in that comprise the steps：
1) GH4169 alloy monotonic tension experiment：
Purpose is to obtain GH4169 alloy monotonic tension performance at ambient temperature, is GH4169 material rolling constitutive model material
The determination of material parameter provides data；
2) GH4169 alloy low circulation strain controlling asymmetrical load experiment：
Purpose be the tension and compression circulation experiment that obtains GH4169 alloy at room temperature, in the range of low circulation number of times, symmetric strained should
Power one strains retardant curve, determines the material ginseng of isotropic hardening and kinematic hardening in GH4169 material rolling constitutive model
Number；
3) Parameters of constitutive model identification
From YoshidaUemori model, constitutive relation model parameter is calculated by software；
4) model checking
After completing secondary development to material constitutive model, the forging process using cuboid is verified；
The method of described GH4169 alloy monotonic tension experiment is：
Feeding from GH4169 alloy cold bar, is processed intoWithStandard tensile specimen carry out static(al)
Stretching experiment；
The method of described GH4169 alloy low circulation strain controlling asymmetrical load experiment is：By GH4169 alloy cold pieces of bar stock
It is processed intoStandard lowcycle fatigue sample, carry out low circulation strain controlling asymmetrical load experiment；Range of strain is respectively
For 0.4%, 0.8%, 1.2%, 1.6%, 2.0%, after obtaining stressstrain curve and the stable circulation of the first loaded cycle
Stress one strain lag loop；
Described Parameters of constitutive model knows method for distinguishing：
From YoshidaUemori model mainly comprise the following steps：Dynamics, yield function are accurate to related flowing
Then, isotropic hardening and kinematic hardening；The particular content of each several part is respectively：
1. kinetics equation
In formula,For total strain rate,WithIt is respectively elastic strain speed and moulding strain rate；
2. yield function and related flowing criterion
Comprise isotropism one kinematic hardening and the yield function of related flowing criterion is
In formula, S and a is respectively cauchy stress deviator and back stress deviator, Y and R be respectively initial yield stress and each to
Property hardening stress；Stress be defined as material unit and longpending on active force, cauchy stress is exactly the trus stress in analysis on Large Deformation,
It is different from nominal stress；And in the case of small deformation, nominal stress and trus stress numerically equal；
Cauchy stress
σ=F/A
Nominal stress σ '=F/Ao
In formula, A is the instantaneous cross section of deformable material；Ao is the sectional area before material deformation；
3. isotropic hardening equation
The isotropic hardening property of material to be represented by internal variable isotropic hardening stress, the drilling of isotropic hardening stress
Changing equation is
In formula,For equivalent plastic strain speed, Rsat is the saturation value of R in the case of big plastic strain, and m is to control isotropism
The material parameter of hardening rate；Above EVOLUTION EQUATION is integrated, obtains
R=R_{sat}(1e^{mp})
Kinematic hardening
In YoshidaUemori model, back stress α is made up of two parts：α_{1}For nonlinear kinematic hardening back stress, α_{2}For linear
Kinematic hardening back stress；
α=α_{1}+α_{2}
In formula, c, α are and nonlinear servoactuated related material parameter, H '_{∞}For linear kinematic hardening coefficient, α_{1}、α_{2}It is respectively nonthread
Property kinematic hardening back stress and linear kinematic hardening back stress；
In YoshidaUemori model, material parameter undetermined includes：
Y is yield surface radius, i.e. initial yield stress；
R_{sat}Saturation value for isotropic hardening stress；
M is the material parameter controlling isotropic hardening speed；
C is the material parameter controlling nonlinear kinematic hardening speed：
α is the material parameter related to nonlinear kinematic hardening；
H′_{∞}For linear kinematic hardening coefficient；
Using constitutive model as 6 material parameters function f (Y, R_{sat}, m, a, C, H '_{∞}), using Matlab software to GH4169
The stressstrain curve that gained tested by alloy carries out curve fitting, you can obtain Y, R_{sat}, m, a, C, H '_{∞}Value, using this side
Method tries to achieve the material parameter of GH4169 alloy state before rolling.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201310499514.1A CN103605830B (en)  20131019  20131019  Establishing method of constructive model for rolling simulation of GH4169 materials 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201310499514.1A CN103605830B (en)  20131019  20131019  Establishing method of constructive model for rolling simulation of GH4169 materials 
Publications (2)
Publication Number  Publication Date 

CN103605830A CN103605830A (en)  20140226 
CN103605830B true CN103605830B (en)  20170208 
Family
ID=50124051
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201310499514.1A Active CN103605830B (en)  20131019  20131019  Establishing method of constructive model for rolling simulation of GH4169 materials 
Country Status (1)
Country  Link 

CN (1)  CN103605830B (en) 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN109372826A (en) *  20181213  20190222  中国工程物理研究院化工材料研究所  A kind of linear continuous static loading control method 
Families Citing this family (12)
Publication number  Priority date  Publication date  Assignee  Title 

CN104376201B (en) *  20141107  20170531  沈阳黎明航空发动机(集团)有限责任公司  A kind of determination method of the key process parameter of influence alloy material mechanical property 
CN104732012B (en) *  20150213  20180814  西北工业大学  The method for building up of GH4169 alloy forged pieces grain size and forging thermal parameter relationship 
CN105468829B (en) *  20151119  20190322  东华理工大学  A kind of calculating and detection method of Strength of Metallic Materials and index based on broad sense Hosford yield function 
CN107305174B (en) *  20160420  20201211  中国特种设备检测研究院  Numerical representation method and system for material stressstrain constitutive relation 
CN106227928A (en) *  20160720  20161214  福州大学  The method for numerical simulation of anisotropy of wood plastic yielding constitutive model 
CN106610357B (en) *  20161202  20190426  燕山大学  A kind of measuring method of metal material yield surface 
CN106834642B (en) *  20170118  20210223  抚顺特殊钢股份有限公司  Optimized forging process of GH6783 alloy bar 
CN106769555B (en) *  20170119  20190419  北京工业大学  A kind of high temperature Multiaxial stress strain stress relation modeling method under tensiontorsion load 
CN107633120B (en) *  20170907  20180717  东南大学  A kind of construction method of fibre reinforced composites dynamic shearing constitutive model 
CN107748817B (en) *  20171019  20210409  北京工业大学  Hightemperature multiaxis constitutive relation determination method considering nonproportional additional reinforcement 
CN108254250B (en) *  20171230  20200508  北京工业大学  Heat engine multiaxis stressstrain relation determination method considering dynamic strain aging 
CN110826285A (en) *  20191122  20200221  华东理工大学  Fatigue life prediction method for geometric discontinuous structure 

2013
 20131019 CN CN201310499514.1A patent/CN103605830B/en active Active
NonPatent Citations (4)
Title 

Multiaxial fatigue damage parameter and life prediction under low cycle loading for GH4169 alloy and other structural materials;GuoQin Sun,et al.;《International Journal of Fatigue》;20091221;全文 * 
X12CrMoWVNbN1011转子钢室温低周疲劳特性;吴海利 等;《北京科技大学学报》;20110731;第33卷(第7期);摘要，第1，3节 * 
以ZenerHollomon参数表示的GH169合金的本构关系;刘东，罗子健;《塑性工程学报》;19950331;全文 * 
高强钢回弹预测中材料模型的适用性研究及回弹补偿的自;张璐;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20120715;摘要，第2.1，2.4，3.13.2节 * 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN109372826A (en) *  20181213  20190222  中国工程物理研究院化工材料研究所  A kind of linear continuous static loading control method 
Also Published As
Publication number  Publication date 

CN103605830A (en)  20140226 
Similar Documents
Publication  Publication Date  Title 

CN103605830B (en)  Establishing method of constructive model for rolling simulation of GH4169 materials  
US10089422B2 (en)  Stressstrain relation simulation method, springbackamount prediction method, and springback analyzer  
CN102221503B (en)  Singleshaft tensile overall true stresstrue strain curve testing technique  
Taherizadeh et al.  Finite element simulation of springback for a channel draw process with drawbead using different hardening models  
CN102175511B (en)  Method and system for estimating material property  
CN102455263B (en)  Method for obtaining mechanical property of metal material based on loaddepth curve  
CN103886125B (en)  A kind of titanium alloy hot combined shaping method for numerical simulation  
CN103745114B (en)  Method for computing stress relaxation numerical values and resilience of titanium alloy  
DucToan et al.  A modified Johnson–Cook model to predict stressstrain curves of boron steel sheets at elevated and cooling temperatures  
Yang et al.  Flow behavior and processing maps of a lowcarbon steel during hot deformation  
Campos et al.  Hydraulic bulge test for stressstrain curve determination and damage calibration for ItoGoya model  
CN106644711A (en)  Test method for uniaxial constitutive relation of ductile material  
Hochholdinger et al.  Determination of flow curves by stack compression tests and inverse analysis for the simulation of hot forming  
CN107292029A (en)  A kind of determination method that sheet forming technological parameter is predicted based on forming defects  
CN107271273B (en)  A kind of method of cross tensile predeformation load measurement yield surface  
Chen  A general and unified theory of the kinematic relationships in bulge tests  
Barzegar et al.  Finite element modeling of damage evolution in cold pilgering process  
Huang et al.  Fracture Prediction and Correlation of AlSi Hot Stamped Steels with Different Models in LS DYNA  
CN107609223A (en)  A kind of coldrolled biphase steel dynamic deformation constitutive model method for building up of tensile strength 1200MPa  
Wang et al.  Probabilistic design of uncertainty for aluminum alloy sheet in rubber fluid forming process  
CN106525612A (en)  Construction method of polycarbonate uniform constitutive model based on tensile tests and compression tests  
Altenbach et al.  Short Term Transversally Isotropic Creep of Plates Under Static and Periodic Loading  
CN107818184A (en)  A kind of method for building material deformation resistance model  
Du et al.  Effect of Different Yield Criteria and Material Parameter Identification Methods on the Description Accuracy of the Anisotropic Behavior of 5182O Aluminum Alloy  
Chen et al.  Analysis of the Microtype Circular Rod Screw Thread Machining Process 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
C14  Grant of patent or utility model  
CP01  Change in the name or title of a patent holder 
Address after: 110043 Dong TA street, Dadong District, Shenyang, Liaoning Province, No. 6 Patentee after: Chinese Hangfa Shenyang Liming Aero engine limited liability company Address before: 110043 Dong TA street, Dadong District, Shenyang, Liaoning Province, No. 6 Patentee before: Liming Aeroplane Engine (Group) Co., Ltd., Shenyang City 

CP01  Change in the name or title of a patent holder 