CN104928605B - 一种预测镍基合金高温流变应力和动态再结晶行为的方法 - Google Patents
一种预测镍基合金高温流变应力和动态再结晶行为的方法 Download PDFInfo
- Publication number
- CN104928605B CN104928605B CN201510424926.8A CN201510424926A CN104928605B CN 104928605 B CN104928605 B CN 104928605B CN 201510424926 A CN201510424926 A CN 201510424926A CN 104928605 B CN104928605 B CN 104928605B
- Authority
- CN
- China
- Prior art keywords
- base alloy
- dynamic recrystallization
- nickel
- strain
- stress
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
本发明公开了一种预测镍基合金高温流变应力和动态再结晶行为的方法,其方法的步骤包括:(1)通过高温压缩实验,获得镍基合金的真应力‑真应变数据;(2)建立预测镍基合金高温流变应力和动态再结晶行为的统一本构模型;(3)利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力‑真应变数据,确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数;(4)预测恒温恒应变速率和变温变应变速率条件下的镍基合金高温流变应力和动态再结晶行为。本发明方法可快速地准确预测恒温恒应变速率和变温变应变速率条件下镍基合金高温流变应力和动态再结晶行为,对合理制定镍基合金热加工工艺有重要的技术指导意义。
Description
技术领域:
本发明属于镍基合金加工工程技术领域,特别涉及一种预测镍基合金高温流变应力和动态再结晶行为的方法。
背景技术:
在镍基合金热加工过程中,镍基合金的热变形过程通常可以分为弹性变形和塑性变形两个阶段。镍基合金的弹性变形阶段通常可以通过胡克定律准确描述。当外加载荷超过了镍基合金的屈服应力,在宏观尺度上表现为镍基合金的塑性变形开始发生,在微观尺度上表现为镍基合金内部的位错运动过程。由于位错的产生和增殖导致的加工硬化行为,进一步促进了镍基合金真应力的增加;随着变形程度的增加,空位逐步扩散,位错滑移和攀移引起的位错相消和位错重排的动态回复过程开始发生,镍基合金真应力增长速度变慢。镍基合金作为典型低层错能金属,位错攀移和交滑移等动态回复机制发生缓慢,难以与加工硬化达到平衡,随着应变的增加,位错密度逐渐上升,当达到发生动态再结晶的临界位错密度时,动态再结晶开始发生,位错湮灭速率显著上升,流变软化现象明显,最终位错湮灭速率与增殖速率达到平衡,镍基合金的真应力达到稳态。研究表明镍基合金的热变形行为极为复杂,显著受到变形温度,应变速率和应变等宏观热变形条件的综合影响,又受到加工硬化,动态回复和动态再结晶等微观变形机制的影响。众多学者开展了大量实验和理论研究工作,发明了多种预测镍基合金高温流变应力和动态再结晶行为的方法。其中,Arrhenius模型、Cingara模型及相关修正模型能够准确预测恒温恒应变速率等理想热变形状态下的镍基合金高温流变应力,Avrami模型,Laasraoui-Jonas模型及相关修正模型能够准确预测恒温恒应变速率等理想热变形状态下的镍基合金动态再结晶行为,但都难以推广应用到具有时变变形特征(如,变温变应变速率)的工业实际热加工过程。然而国内外尚无见到过基于镍基合金热变形物理机理,提出能够预测时变变形条件下镍基合金高温流变应力和动态再结晶行为的方法。
因此,本发明从镍基合金热变形物理机理出发,发明了一种可快速地准确预测恒温恒应变速率和变温变应变速率条件下镍基合金高温流变应力和动态再结晶行为的方法,以解决现有预测方法应用范围狭窄,难以工程推广的弊端。该方法的发明和推广应用对合理制定镍基合金热加工工艺有重要的技术指导意义。
发明内容:
本发明的目的在于提供一种预测镍基合金高温流变应力和动态再结晶行为的方法,解决了现有预测方法应用范围狭窄,难以实现工程应用的弊端,对合理制定镍基合金热加工工艺有重要的技术指导意义。
为达到上述目的,本发明采用的技术方案是:一种预测镍基合金高温流变应力和动态再结晶行为的方法。该方法的具体步骤为:
步骤1:在变形温度为900℃~1100℃和应变速率为0.0005s-1~10s-1的热变形条件下,对初始晶粒尺寸为20μm~90μm的镍基合金进行高温压缩实验,获得镍基合金的真应力-真应变数据;
步骤2:建立预测镍基合金高温流变应力和动态再结晶行为的统一本构模型:
σ=σy+σi+σg (1)
其中σ为高温流变应力,σy为屈服应力,σi为位错密度引起的应力,σg为晶粒尺寸演化引起的应力;
建立镍基合金的屈服应力σy与变温温度T、应变速率初始晶粒尺寸d0之间的函数关系:
其中Ay、my、ny和Qy均为材料参数,R为普适气体常数,为应变速率,T为变形温度,d0为初始晶粒尺寸;根据高温压缩实验的变形条件和镍基合金的真应力-真应变数据,绘制镍基合金的屈服应力σy和变温温度T、应变速率初始晶粒尺寸d0之间的关系图,即lnσy-1/T和lnσy-ln d0关系图,并通过线性拟合的方法确定材料参数Ay、my、ny和Qy的具体数值;
建立镍基合金的位错密度与其引起的应力σi的函数关系:
其中ρi为位错密度,为位错密度演化速率,M为泰勒系数,α为位错交互作用常数,μ为材料剪切模量,b为柏氏矢量;位错增殖速度fw为加工硬化系数,动态回复导致的位错湮灭速度fv为动态回复系数,动态再结晶导致的位错湮灭速度σi-i(cr)为当前应变减去临界应变状态下的位错密度,fx为动态再结晶系数,Xi为动态再结晶份数,当变形大于动态再结晶临界应变εc时,ε0.5为当动态再结晶份数Xi达到0.5时的应变,fd为动态再结晶动力学指数,Aw、Av、Ax、Ad、mw、mv、mx、md、nw、nv、nx、nd、Qw、Qv、Qx和Qd均为材料参数,ρcr为动态再结晶临界位错密度,s为晶界能,s=μbθm/4π(1-v),θm为晶界取向角,v为泊松比,L为位错自由程,L=kdb(μ/σy)m,kd和m为材料常数,λ为位错线能量,λ=cμb2,c为材料常数,Mbm为晶界可动性,Mbm=vδDobexp(-Qdrf/RT)/kbT,δ为特征晶界厚度,Dob为晶界自扩散系数,kb为波尔茨曼常量,Qdrf晶界扩散激活能;
建立镍基合金的晶粒尺寸与其演化引起的应力σg之间的函数关系:
其中fg为晶粒尺寸演化系数,di为平均晶粒尺寸,ddrx为动态再结晶晶粒尺寸,Ag、Ar、mg、mr、ng、nr、Qg和Qr均为材料参数;根据高温压缩实验的变形条件和高温压缩实验后镍基合金的金相实验统计数据,绘制镍基合金的动态再结晶晶粒尺寸ddrx和变温温度T、应变速率初始晶粒尺寸d0之间的关系图,即ln ddrx-lnln ddrx-1/T和ln ddrx-ln d0关系图,并通过线性拟合的方法确定材料参数Ar、mr、nr和Qr的具体数值;
步骤3:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数;
利用数值差分原理,将与任意小应变增量Δε引起的位错密度增量Δρi表示为 将与任意小应变增量Δε引起的位错密度演化引起的应力增量Δσi表示为Δσi=(αMμbρi -1/2)Δρi/2,将与任意小应变增量Δε引起的晶粒尺寸演化引起的应力增量Δσg表示为将与任意小应变增量Δε引起的动态再结晶份数增量ΔXi表示为 将与任意小应变增量Δε引起的晶粒尺寸增量Δdi表示为Δdi=-(d0-ddrx)ΔXi,则任意小应变增量Δε引起的应力增量Δσ确定为Δσ=αMμbρi -1/2/2·Δρi-fg(1+Xi(d0-ddrx)/2di)di -1/2ΔXi;
编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,对预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数Aw,Av,Ax,Ag,Ad,mw,mv,mx,mg,md,nw,nv,nx,ng,nd,Qw,Qv,Qx,Qg和Qd进行优化求解,整个优化过程可以分为三步:第一步,根据位错密度演化规律,即确定与动态再结晶临界位错密度ρcr对应的临界应变εc;第二步,确定当动态再结晶份数Xi达到0.5时的应变ε0.5,动态再结晶份数同样可采用Xi=(σrec-σ)/(σsat-σss)进行计算,其中σsat为饱和应力,且 σss稳态应力,且 σrec为假定整个热变形过程中未有动态再结晶发生,只受加工硬化和动态回复机制控制情况下的应力,通过求解(σrec-σ)/(σsat-σss)=0.5,确定当动态再结晶份数Xi达到0.5时的应变ε0.5;第三步,优化确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数Aw,Av,Ax,Ag,Ad,mw,mv,mx,mg,md,nw,nv,nx,ng,nd,Qw,Qv Qx,Qg和Qd的具体数值;
步骤4:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,实现热变形参数及受热变形参数影响的材料参数在任意迭代步更新,进而预测恒温恒应变速率和变温变应变速率条件下的镍基合金高温流变应力和动态再结晶行为,其中热变形参数包括变形温度和应变速率,受热变形参数影响的材料参数包括屈服应力σy,动态再结晶晶粒尺寸ddrx,加工硬化系数fw,动态回复系数fv,动态再结晶系数fx、晶粒尺寸演化系数fg和动态再结晶动力学指数fd。
本发明通过镍基合金高温压缩实验,在位错密度理论和动态再结晶动力学基础上,建立了一种预测镍基合金高温流变应力和动态再结晶行为的模型,充分考虑了实时变形条件对镍基合金高温流变应力和动态再结晶行为的影响,可快速地准确预测恒温恒应变速率和变温变应变速率条件下镍基合金高温流变应力和动态再结晶行为。
本发明的有益效果为:本发明充分考虑了实时变形条件对镍基合金高温流变应力和动态再结晶行为的影响,可快速地准确预测恒温恒应变速率和变温变应变速率条件下镍基合金高温流变应力和动态再结晶行为,能够推广应用到变温变应变速率条件下的工业实际热加工过程中,解决了现有预测方法应用范围狭窄,难以工程推广的弊端。该方法的发明和推广应用对合理制定镍基合金热加工工艺有重要的技术指导意义。
附图说明:
图1优化过程示意图
图2恒温恒应变速率条件下GH4169合金预测结果:(a)高温流变应力;(b)动态再结晶行为;(c)平均晶粒尺寸演变;
图3变温变应变速率条件下GH4169合金预测结果:(a)高温流变应力;(b)动态再结晶行为;(c)平均晶粒尺寸演变;
具体实施方式:
下面结合附图和具体实施方式对本发明进行详细说明。
本发明是一种预测镍基合金高温流变应力和动态再结晶行为的方法,下面以GH4169合金(典型镍基合金)的高温流变应力和动态再结晶行为的预测为例,详细介绍本发明涉及的预测方法的具体实施细节,其方法包括:
步骤1:对GH4169合金进行高温压缩实验,初始晶粒尺寸分别为75μm、48μm和33μm,变形温度分别为920℃、950℃、980℃、1010℃和1040℃,应变速率分别为0.001s-1、0.01s-1、0.1s-1和1s-1,应变量为1.2。
步骤2:建立预测镍基合金高温流变应力和动态再结晶行为的统一本构模型:σ=σy+σi+σg,其中σ为高温流变应力,σy为屈服应力,σi为位错密度引起的应力,σg为晶粒尺寸演化引起的应力;
建立镍基合金的屈服应力σy与变温温度T、应变速率初始晶粒尺寸d0之间的函数关系:
其中Ay、my、ny和Qy均为材料参数,R为普适气体常数(8.314Jmol-1K-1),为应变速率,T为变形温度,d0为初始晶粒尺寸;利用高温压缩实验获得的GH4169合金的真应力-真应变数据,通过0.2%应变补偿法可以测得GH4169合金的屈服应力。根据GH4169合金的真应力-真应变数据,绘制屈服应力σy和变温温度T、应变速率初始晶粒尺寸d0之间的关系图,即lnσy-1/T和lnσy-ln d0关系图。通过数据线性拟合的方法,对lnσy-1/T和lnσy-ln d0关系图中的数据进行回归,确定材料参数Ay、my、ny和Qy的具体数值分别为0.708、-0.123、0.09和663.870kJ/mol。因此,屈服应力和变温温度T、应变速率初始晶粒尺寸d0之间的关系可表示为
建立镍基合金的位错密度与其引起的应力σi的函数关系:
其中ρi为位错密度,初始状态位错密度假定为1×1011m-2,为位错密度演化速率,M为泰勒系数,等于3.06;α为位错交互作用常数,等于0.3;μ为材料剪切模量,与温度显著相关,与温度T的关系可以表示为μ=86.94-0.027T;b为柏氏矢量(2.54×10-10m);位错增殖速度fw为加工硬化系数,动态回复导致的位错湮灭速度fv为动态回复系数,动态再结晶导致的位错湮灭速度ρi-i(cr)为当前应变减去临界应变状态下的位错密度,fx为动态再结晶系数,Xi为动态再结晶份数,当变形大于动态再结晶临界应变εc时,ε0.5为当动态再结晶份数Xi达到0.5时的应变,fd为动态再结晶动力学指数,Aw、Av、Ax、Ad、mw、mv、mx、md、nw、nv、nx、nd、Qw、Qv、Qx和Qd均为材料参数,ρcr为动态再结晶临界位错密度,s为晶界能,s=μbθm/4π(1-v),θm为晶界取向角,θm=15°,v为泊松比,v=0.33,L为位错自由程,L=kdb(μ/σy)m,kd和m为材料常数,(kd=10,m=1.33),λ为位错线能量,λ=cμb2,c为材料常数,c=0.5,Mbm为晶界可动性,Mbm=bδDcbexp(-Qdrf/RT)/kbT,δ为特征晶界厚度,Dob为晶界自扩散系数,(δDob=3.5×10-15m3/s),kb为波尔茨曼常量(1.381×10-23J/K),Qdrf晶界扩散激活能(115kJ/mol);
建立镍基合金的晶粒尺寸与其演化引起的应力σg之间的函数关系:
其中fg为晶粒尺寸演化系数,di为平均晶粒尺寸,ddrx为动态再结晶晶粒尺寸,Ag、Ar、mg、mr、ng、nr、Qg和Qr均为材料参数;
对高温压缩实验后镍基合金进行金相实验,根据截线法对镍基合金金相图片中的动态再结晶晶粒尺寸进行统计。绘制镍基合金的动态再结晶晶粒尺寸ddrx和变温温度T、应变速率初始晶粒尺寸d0之间的关系图,即ln ddrx-1/T和ln ddrx-ln d0关系图。通过数据线性拟合的方法,对ln ddrx-1/T和ln ddrx-ln d0关系图中的数据进行回归,确定材料参数Ar、mr、nr和Qr的具体数值分别为737.151,0.329,-0.139和474kJ/mol。因此,屈服应力和变温温度T、应变速率初始晶粒尺寸d0之间的关系可表示为
步骤3:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数;
利用数值差分原理,将与任意小应变增量Δε引起的位错密度增量Δρi表示为 将与任意小应变增量Δε引起的位错密度演化引起的应力增量Δσi表示为Δσi=(αMμbρi -1/2)Δρi/2,将与任意小应变增量Δε引起的晶粒尺寸演化引起的应力增量Δσg表示为将与任意小应变增量Δε引起的动态再结晶份数增量ΔXi表示为 将与任意小应变增量Δε引起的晶粒尺寸增量Δdi表示为Δdi=-(d0-ddrx)ΔXi,则任意小应变增量Δε引起的应力增量Δσ可表示为Δσ=αMμbρi -1/2/2·Δρi-fg(1+Xi(d0-ddrx)/2di)di -1/2ΔXi;
编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,对预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数Aw,Av,Ax,Ag,Ad,mw,mv,mx,mg,md,nw,nv,nx,ng,nd,Qw,Qv,Qx,Qg和Qd进行优化求解,整个优化过程可以分为三步:第一步,根据位错密度演化规律,即确定与动态再结晶临界位错密度ρcr对应的临界应变εc;第二步,确定当动态再结晶份数Xi达到0.5时的应变ε0.5,动态再结晶份数同样可采用Xi=(σrec-σ)/(σsat-σss)进行计算,其中,σsat为饱和应力,且 σss稳态应力,且 σrec为假定整个热变形过程中未有动态再结晶发生,只受加工硬化和动态回复机制控制情况下的应力,通过求解(σrec-σ)/(σsat-σss)=0.5,确定当动态再结晶份数Xi达到0.5时的应变ε0.5;第三步,优化确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数Aw,Av,Ax,Ag,Ad,mw,mv,mx,mg,md,nw,nv,nx,ng,nd,Qw,QvQx,Qg和Qd的具体数值,如表1所示。优化过程示意图如图1所示。
表1优化得到的GH4169合金参数值.
步骤4:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,实现热变形参数及受热变形参数影响的材料参数在任意迭代步更新,进而预测恒温恒应变速率和变温变应变速率条件下的镍基合金高温流变应力和动态再结晶行为,其中热变形参数包括变形温度和应变速率,受热变形参数影响的材料参数包括屈服应力σy,动态再结晶晶粒尺寸ddrx,加工硬化系数fw,动态回复系数fv,动态再结晶系数fx、晶粒尺寸演化系数fg和动态再结晶动力学指数fd。
图2所示分别为恒温恒应变速率条件下GH4169合金高温流变应力、动态再结晶行为和平均晶粒尺寸演变的预测结果。图3所示分别为变温变应变速率条件下GH4169合金高温流变应力、动态再结晶行为和平均晶粒尺寸演变的预测结果。从图中可以发现,GH4169合金的真应力-真应变的预测值和实验值吻合良好,表明本发明的方法能够准确地预测GH4169合金的高温流变应力和动态再结晶行为。
上面结合附图对本发明的实例进行了描述,但本发明不局限于上述具体的实施方式,上述的具体实施方式仅是示例性的,不是局限性的,任何不超过本发明权利要求的发明创造,均在本发明的保护之内。
Claims (1)
1.一种预测镍基合金高温流变应力和动态再结晶行为的方法,其特征在于:充分考虑了实时变形条件对镍基合金高温流变应力和动态再结晶行为的影响,基于位错密度理论和动态再结晶动力学,提出了一种预测镍基合金高温流变应力和动态再结晶行为的统一本构模型,可快速地准确预测恒温恒应变速率和变温变应变速率条件下镍基合金高温流变应力和动态再结晶行为,该方法包括以下步骤:
步骤1:在变形温度为900℃~1100℃和应变速率为0.0005s-1~10s-1的热变形条件下,对初始晶粒尺寸为20μm~90μm的镍基合金进行高温压缩实验,获得镍基合金的真应力-真应变数据;
步骤2:建立预测镍基合金高温流变应力和动态再结晶行为的统一本构模型:
其中 为高温流变应力, 为屈服应力, 为位错密度引起的应力, 为晶粒尺寸演化引起的应力;
建立镍基合金的屈服应力 与变温温度T 、应变速率初始晶粒尺寸d 0之间的函数关系:
其中A y、m y、n y和Q y均为材料参数,R为普适气体常数,为应变速率,T 为变形温度,d 0为初始晶粒尺寸;根据高温压缩实验的变形条件和镍基合金的真应力-真应变数据,绘制镍基合金的屈服应力和变温温度T 、应变速率初始晶粒尺寸d 0之间的关系图,即 和 关系图,并通过线性拟合的方法确定材料参数A y、m y、n y和Q y的具体数值;
建立镍基合金的位错密度与其引起的应力的函数关系:
其中ρ i 为位错密度,为位错密度演化速率,M 为泰勒系数, 为位错交互作用常数,μ为材料剪切模量,b 为柏氏矢量;位错增殖速度 f w为加工硬化系数,动态回复导致的位错湮灭速度 f v为动态回复系数,动态再结晶导致的位错湮灭速度 ρ i-i(cr)为当前应变减去临界应变状态下的位错密度,f x为动态再结晶系数,,X i 为动态再结晶份数,当变形大于动态再结晶临界应变 时,,为当动态再结晶份数X i 达到0.5时的应变,f d为动态再结晶动力学指数, A w、A v、A x、A d、m w、m v、m x、m d、n w、n v、n x、n d、Q w、Q v、Q x和Q d均为材料参数,ρ cr为动态再结晶临界位错密度, s为晶界能,s=μbθ m/4π(1-v),θ m为晶界取向角,v 为泊松比,L为位错自由程,L=k d b(μ /) m ,k d和m为材料常数,λ 为位错线能量,λ =cμb 2,c为材料常数,M bm为晶界可动性,M bm=bδD ob exp(-Q dif/RT)/k b T,δ为特征晶界厚度,D ob为晶界自扩散系数,k b为波尔茨曼常量,Q dif晶界扩散激活能;
建立镍基合金的晶粒尺寸与其演化引起的应力之间的函数关系:
其中f g为晶粒尺寸演化系数,d i为平均晶粒尺寸,d drx为动态再结晶晶粒尺寸,A g、A r、m g、m r、n g、n r、Q g和Q r均为材料参数;根据高温压缩实验的变形条件和高温压缩实验后镍基合金的金相实验统计数据,绘制镍基合金的动态再结晶晶粒尺寸d drx和变温温度T、应变速率初始晶粒尺寸d 0之间的关系图,即ln d drx-1/T 和ln d drx-ln d 0关系图,并通过线性拟合的方法确定材料参数A r、m r、n r和Q r的具体数值;
步骤3:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数;
利用数值差分原理,将与任意小应变增量Δε引起的位错密度增量Δρ i 表示为将与任意小应变增量Δε引起的位错密度演化引起的应力增量Δ表示为,将与任意小应变增量Δε引起的晶粒尺寸演化引起的应力增量Δ表示为 将与任意小应变增量Δε引起的动态再结晶份数增量ΔX i 表示为 将与任意小应变增量Δε引起的晶粒尺寸增量Δd i 表示为Δd i =-(d 0-d drx)ΔX i ,则任意小应变增量Δε引起的应力增量Δ确定为 ;
编写迭代累加算法程序,嵌入数值模拟软件,结合镍基合金的真应力-真应变数据,对预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数A w,A v,A x,A g,A d,m w,m v,m x,m g,m d,n w,n v,n x,n g,n d,Q w,Q v,Q x,Q g和Q d进行优化求解,整个优化过程可以分为三步:第一步,根据位错密度演化规律,即确定与动态再结晶临界位错密度ρ cr对应的临界应变 ;第二步,确定当动态再结晶份数X i 达到0.5时的应变,动态再结晶份数同样可采用进行计算,其中 为饱和应力,且 为稳态应力,且 为假定整个热变形过程中未有动态再结晶发生,只受加工硬化和动态回复机制控制情况下的应力,通过求解,确定当动态再结晶份数X i 达到0.5时的应变;第三步,优化确定预测镍基合金高温流变应力和动态再结晶行为的统一本构模型的材料参数A w,A v,A x,A g,A d,m w,m v,m x,m g,m d,n w,n v,n x,n g,n d,Q w,Q v,Q x,Q g和Q d的具体数值;
步骤4:利用数值差分原理,编写迭代累加算法程序,嵌入数值模拟软件,实现热变形参数及受热变形参数影响的材料参数在任意迭代步更新,进而预测恒温恒应变速率和变温变应变速率条件下的镍基合金高温流变应力和动态再结晶行为,其中热变形参数包括变形温度和应变速率,受热变形参数影响的材料参数包括屈服应力,动态再结晶晶粒尺寸d drx,加工硬化系数f w,动态回复系数f v,动态再结晶系数f x、晶粒尺寸演化系数f g和动态再结晶动力学指数f d。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510424926.8A CN104928605B (zh) | 2015-07-20 | 2015-07-20 | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510424926.8A CN104928605B (zh) | 2015-07-20 | 2015-07-20 | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104928605A CN104928605A (zh) | 2015-09-23 |
CN104928605B true CN104928605B (zh) | 2017-01-18 |
Family
ID=54116025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510424926.8A Expired - Fee Related CN104928605B (zh) | 2015-07-20 | 2015-07-20 | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104928605B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105426595A (zh) * | 2015-11-05 | 2016-03-23 | 内蒙古工业大学 | 一种铝合金热弹塑性变形模拟用本构模型的建立方法 |
CN105653822B (zh) * | 2016-01-29 | 2018-11-09 | 中南大学 | 一种模拟gh4169合金静态再结晶行为的元胞自动机方法 |
CN106053754B (zh) * | 2016-07-06 | 2017-11-14 | 中南大学 | 一种预测时变工况下高合金化材料动态再结晶分数的方法 |
CN107784149B (zh) * | 2016-09-05 | 2023-08-08 | 清华大学 | 一种bcj本构模型的参数识别方法及装置 |
CN106503397B (zh) * | 2016-11-16 | 2019-06-04 | 中国航空工业集团公司北京航空材料研究院 | 一种获得金属材料再结晶组织演化晶界可动性参数的方法 |
CN108107071B (zh) * | 2016-11-25 | 2020-06-16 | 中国科学院金属研究所 | 一种单晶高温合金再结晶倾向的评价方法 |
CN107220485B (zh) * | 2017-05-12 | 2018-02-23 | 华中科技大学 | 一种适用于多道次压缩的本构模型的建立方法 |
CN107245683B (zh) * | 2017-06-08 | 2018-09-04 | 太原钢铁(集团)有限公司 | 镍基合金径锻组织控制方法 |
CN108169019B (zh) * | 2017-12-26 | 2020-06-16 | 中国兵器工业第五九研究所 | 一种准静态塑性压缩应力应变参量的识别方法 |
CN108595827A (zh) * | 2018-04-20 | 2018-09-28 | 安徽工业大学 | 一种C-Mn-Al高强度钢热变形组织演变机制及热加工性能的确定方法 |
CN108563880B (zh) * | 2018-04-20 | 2022-01-28 | 中国兵器工业第五九研究所 | 一种计及应变硬化和软化的金属材料本构模型建立方法 |
CN109030496A (zh) * | 2018-06-27 | 2018-12-18 | 燕山大学 | 一种立方晶系中位错柏氏矢量沿任意晶向投影矢量的计算方法 |
CN109147873B (zh) * | 2018-07-13 | 2021-09-21 | 江西理工大学 | 一种预测微合金钢焊接粗晶区晶粒尺寸的方法 |
CN109255175B (zh) * | 2018-08-31 | 2021-04-02 | 西北工业大学 | 一种基于临界再结晶残余应力的再结晶控制方法 |
CN110008620B (zh) * | 2019-04-15 | 2023-06-16 | 中国科学院宁波材料技术与工程研究所 | 一种分析动态载荷条件下α-Fe应变率敏感系数的方法 |
CN110684938A (zh) * | 2019-08-28 | 2020-01-14 | 中南大学 | 一种预测变应变速率工况下金属或合金材料动态再结晶晶粒尺寸的方法 |
FR3100144B1 (fr) * | 2019-09-04 | 2021-10-01 | Safran Aircraft Engines | Procede de fabrication d’une piece metallique limitant l’apparition de grains recristallises dans ladite piece |
CN111380899A (zh) * | 2019-11-29 | 2020-07-07 | 中国科学院金属研究所 | 一种通过轧制模拟过程温升修正锆合金流变应力的方法 |
CN111337356A (zh) * | 2020-04-30 | 2020-06-26 | 西北工业大学 | 一种基于动态再结晶速率的金属热加工性能表征方法 |
CN112877628B (zh) * | 2021-01-13 | 2021-09-21 | 重庆大学 | 低能晶界密度与晶粒尺寸的协调优化方法与系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1718802A (zh) * | 2005-08-05 | 2006-01-11 | 武汉大学 | 镍基高温合金成分的优化设计方法 |
CN1851720A (zh) * | 2006-03-20 | 2006-10-25 | 上海发电设备成套设计研究所 | 一种汽轮机高可用性的设计方法及其评价方法 |
CN103246774A (zh) * | 2013-05-13 | 2013-08-14 | 天津大学 | 数值模拟p92钢管焊接热影响区宽度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140044951A1 (en) * | 2012-08-09 | 2014-02-13 | United Technologies Corporation | High strength-to-density nanocellular foam |
-
2015
- 2015-07-20 CN CN201510424926.8A patent/CN104928605B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1718802A (zh) * | 2005-08-05 | 2006-01-11 | 武汉大学 | 镍基高温合金成分的优化设计方法 |
CN1851720A (zh) * | 2006-03-20 | 2006-10-25 | 上海发电设备成套设计研究所 | 一种汽轮机高可用性的设计方法及其评价方法 |
CN103246774A (zh) * | 2013-05-13 | 2013-08-14 | 天津大学 | 数值模拟p92钢管焊接热影响区宽度的方法 |
Non-Patent Citations (2)
Title |
---|
42CrMo钢亚动态再结晶行为研究;蔺永诚等;《材料热处理学报》;20090430;第30卷(第2期);第70-74页 * |
形变温度对42CrMo钢塑性成形与动态再结晶的影响;蔺永诚等;《材料热处理学报》;20090228;第30卷(第1期);第71-75页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104928605A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104928605B (zh) | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 | |
Fournier et al. | Micromechanical model of the high temperature cyclic behavior of 9–12% Cr martensitic steels | |
CN106503397B (zh) | 一种获得金属材料再结晶组织演化晶界可动性参数的方法 | |
CN105930553A (zh) | 镍基高温合金缺口拉伸强度预测方法 | |
Wedberg et al. | Modelling high strain rate phenomena in metal cutting simulation | |
Qian et al. | Correction of flow stress and determination of constitutive constants for hot working of API X100 pipeline steel | |
Yan et al. | Effects of UMo irradiation creep on the thermo-mechanical behavior in monolithic UMo/Al fuel plates | |
CN105653822A (zh) | 一种模拟gh4169合金静态再结晶行为的元胞自动机方法 | |
Gan et al. | Atomic insights into the effects of Al element on the nanoindentation behavior of single-crystal FeNiCoCr-based multicomponent alloys | |
CN104951633B (zh) | 一种预测镍基合金加工硬化和动态回复行为的方法 | |
Rudd et al. | Multiscale strength (MS) models: their foundation, their successes, and their challenges | |
Jin et al. | Anisotropic cyclic plasticity modeling for additively manufactured nickel‐based superalloys | |
Stewart et al. | Characterization of the Creep Deformation and Rupture Behavior of DS GTD-111 Using the Kachanov–Rabotnov Constitutive Model | |
Xie et al. | Study of the effect of crystal orientation on the mechanical response and fatigue life of solder joints under thermal cycling by crystal plasticity | |
US11320353B2 (en) | Creep strength analysis and assessment method and computer device | |
Scholz et al. | Modeling of mechanical properties of alloy CMSX-4 | |
Zaletelj et al. | Numerical methods for TMF cycle modeling | |
CN111429981B (zh) | 针对单晶材料热-机疲劳本构模型的损伤计算方法 | |
Giordana et al. | Micromechanical modeling of the cyclic softening of EUROFER 97 steel | |
Mahajan et al. | Finite element analysis of printed circuit heat exchanger core for creep and creep-fatigue responses | |
Liang et al. | Time-dependent crack behavior in an integrated structure | |
Lurdos et al. | Empirical and physically based flow rules relevant to high speed processing of 304L steel | |
Yang et al. | Process emulation for predicting die shift and wafer warpage in wafer reconstitution | |
Ebner et al. | Methodology for advanced tool load analysis and lifetime prediction of tools | |
Bérard et al. | Viscoplastic behavior of a FeCrAl alloy for high temperature steam electrolysis (HTSE) sealing applications between 700° C and 900° C |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170118 Termination date: 20170720 |