CN104809273B - 一种描述蠕变变形的方法 - Google Patents
一种描述蠕变变形的方法 Download PDFInfo
- Publication number
- CN104809273B CN104809273B CN201510158190.4A CN201510158190A CN104809273B CN 104809273 B CN104809273 B CN 104809273B CN 201510158190 A CN201510158190 A CN 201510158190A CN 104809273 B CN104809273 B CN 104809273B
- Authority
- CN
- China
- Prior art keywords
- msub
- creep
- mrow
- eta
- curve
- 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)
Abstract
一种描述蠕变变形的方法,它是一种能完整描述3个阶段蠕变变形的方法,步骤如下:一、将试验得到的不同温度和应力下的蠕变曲线归一化;二、对曲线的第2阶段进行拟合,得到η1,η2,η3;三、拟合得到η4,η5;四、将ηi该i=1,2,3,4,5,表示为温度和应力的函数,将ηi该i=1,2,3,4,5表示为无量纲应力σ/σ0.2与无量纲温度T/Tm的函数;五、将该方法与有限元软件结合,编写usercreep子程序;六、选择适当的适用于变载情况下的蠕变模型;七、建立实际结构的有限元模型,进行蠕变变形计算和应力松弛行为分析;所发明的蠕变变形方法,解决了当前常用的蠕变模型描述蠕变曲线的不足,同时针对变载情况,在子程序中实现3种变载模型。
Description
一.技术领域
本发明提供一种能完整描述3个阶段蠕变变形的模型及将其应用于实际结构的蠕变计算的方法,属于高温结构蠕变模拟技术领域。
二.背景技术
蠕变是指在温度、载荷不变的条件下,材料的变形也会随时间增长而逐渐增大的现象,而且这种变形即使在应力小于屈服极限时仍具有不可逆的变形性质。典型的蠕变变形可以分为3个阶段:第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低,称为初始蠕变阶段;第2阶段为直线,这是由于加工硬化与回复软化过程达到动态平衡而导致蠕变速率保持不变,即为稳态蠕变阶段;第3阶段蠕变速率随时间增大直至断裂,称为加速蠕变阶段。
当前工程中常用的蠕变模型有时间硬化模型和应变硬化模型等,但这类模型只能模拟蠕变第1阶段或者前两个阶段(有限元软件ANSYS和ABAQUS提供的蠕变模型均不能模拟蠕变第3阶段);1985年Evans和Wilshire在其著作中提出一种能模拟完整的蠕变变形的模型—θ投射法,虽然θ投射法可以模拟蠕变变形的3个阶段,但它并不注重描述稳态蠕变阶段,即它所描述的整个蠕变变形中,并不存在蠕变速率为常数的阶段;一些粘塑性本构模型虽能较好地描述循环载荷下的塑形变形,而在描述蠕变第3阶段方面必须引入损伤参量对模型进行修正。这些模型用于试样在试验载荷下变形的描述较精准,但由于计算较为复杂,目前还不易用于实际结构(如涡轮叶片和轮盘)的蠕变分析。
因此,仍有必要发展能够完整描述蠕变曲线的3个阶段的方法,同时能与有限元结合用于实际结构的蠕变变形计算及应力松弛问题的分析。
三.发明内容
1.发明目的
本发明提供一种能完整描述3个阶段蠕变变形的模型及将其应用于实际结构的蠕变计算的方法,达到能计算实际结构3个阶段的蠕变变形和模拟应力松弛行为的目的,解决现有技术的不足。
2.技术方案
本发明提供一种描述蠕变变形的方法,它能完整描述3个阶段(第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低,称为初始蠕变阶段;第2阶段蠕变速率不随时间变化,即稳态蠕变阶段;第3阶段蠕变速率随时间增大直至断裂,称为加速蠕变阶段)蠕变变形的方法,同时通过将该方法与有限元软件ANSYS结合,编写usercreep子程序,用于实际结构的蠕变变形计算。同时提出3种适用于变载情况的模型,用于计算应力松弛行为。本发明用归一化参数描述蠕变变形,将蠕变应变表示为:εc=η1ζη4+η2ζ+η3ζη5或分别称之为表达式1和表达式2,其中3项分别表示蠕变的3个阶段。表达式中各参数含义如下:ζ=t/tc为无量纲时间,tc为给定温度和应力下的持久寿命,则ζ∈[0,1];ηi(i=1,2,3,4,5)为材料参数,η1,η2,η3,分别为蠕变3个阶段的蠕变量,η4,η5分别控制蠕变第1阶段和第3阶段的变化快慢,且η5>1。
本发明一种描述蠕变变形的方法,它是一种能完整描述3个阶段蠕变变形的方法,其具体步骤如下:
步骤一:将试验得到的不同温度和应力下的蠕变曲线归一化;
用试验得到的不同温度和应力下的蠕变曲线的时间坐标除以该温度和应力下的持久断裂时间,则所有曲线的横坐标均为归一化时间坐标ζ=t/tc,ζ∈[0,1];
步骤二:对曲线的第2阶段进行拟合,得到η1,η2,η3;
对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程;得到的一组k值则为不同应力和温度下的一组η2值,即归一化坐标下的稳态蠕变应变率或最小蠕变率;拟合直线中的b值即η1,继而可得到η3=εr-η1-η2,εr为断裂时的蠕变应变;
步骤三:拟合得到η4,η5;
然后由已得到的参数η1,η2,η3,根据表达式或对每条曲线进行拟合(可采用matlab等软件进行最小二乘拟合),得到η4,η5;
步骤四:将ηi(i=1,2,3,4,5)表示为温度和应力的函数;
将ηi(i=1,2,3,4,5)表示为无量纲应力σ/σ0.2与无量纲温度T/Tm的函数;
步骤五:将该方法与有限元软件结合,编写usercreep子程序;
利用通用有限元软件ANSYS中usercreep子程序,将步骤一~步骤四得到的表达式写入子程序中,以达到计算实际结构蠕变变形的目的;对所编写的子程序进行编译连接后,在ANSYS主程序中调用子程序即可用所发明的方法对实际结构进行蠕变变形计算;
步骤六:选择适当的适用于变载情况下的蠕变模型;
在usercreep子程序中实现3种不同的适用于变载情况下蠕变模型:时间硬化模型,相对时间硬化模型,应变硬化模型;时间硬化模型为当载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下t时刻后的蠕变应变曲线上下平移得到;相对时间硬化模型为当载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下t2=t×tc,1/tc,2时刻后的蠕变曲线平移得到,tc,1和tc,2分别表示σ1,T1和σ2,T2状态下的持久寿命;应变硬化模型为载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下产生与前一状态相同的蠕变应变所对应的时间之后曲线左右平移得到;
步骤七:建立实际结构的有限元模型,进行蠕变变形计算和应力松弛行为分析;
对实际结构在CAE前处理软件中建立有限元模型;由于应力松弛行为属于变载条件下的蠕变行为,能通过所编写的usercreep子程序进行计算分析;在ANSYS主程序中调用所编写的usercreep子程序,同时输入所求得模型参数值,对实际结构进行有限元数值模拟,计算蠕变变形和应力松弛行为。
其中,在步骤二中所述的“对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程”,其方法如下:利用第二阶段的试验数据点,可采用excel或matlab等软件对该数据点进行一次多项式拟合;若蠕变曲线的第2阶段并不明显,可采用多次多项式对蠕变曲线进行拟合,然后对多次多项式求导,导数值最小的点即为最小蠕变率点,导数值即η2。
其中,在步骤四中所述的“将ηi(i=1,2,3,4,5)表示为温度和应力的函数”,具体做法如下:选取特定函数将无量纲应力σ/σ0.2与无量纲温度T/Tm表示为ηi=f(σ/σ0.2,T/Tm),例如:
或
其中,ai,bi,ci,di(i=1,2,3,4,5)为材料相关系数,可由进行最小二乘拟合(可采用matlab等软件)得到;通过一组ai,bi,ci,di的系数值,则可求得任意温度和应力下的ηi,继而得到该温度和应力下的蠕变曲线。
其中,在步骤五中所述的“usercreep子程序”,其编写所需的输出变量为:蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1),蠕变应变增量对蠕变应变的导数dcrda(2);表达式2的各输出变量为:
蠕变应变增量delcr:
蠕变应变增量对等效应力的导数dcrda(1):
式中,
蠕变应变增量对蠕变应变的导数dcrda(2):
其中,在步骤六中所述的“变载情况下的蠕变模型”,在子程序中实现方法为:
在子程序中,根据t2=t×tc,1/tc,2,可计算得到相对时间硬化模型中的t2;通过得到主程序返回当前时刻的应变值ε,根据表达式2可迭代求得在σ2,T2状态下的蠕变曲线产生应变ε的对应的时间t2′。分别用t2和t2′替换usercreep子程序中输出变量中的t2,即可完成符合相对时间硬化模型和应变硬化模型的usercreep子程序。
3.优点功效
所发明的基于归一化参数的蠕变模型能够模拟完整的3个阶段的蠕变变形,解决当前常用的蠕变模型描述蠕变曲线的不足。通过编写usercreep子程序,能够将该模型用于实际结构的蠕变变形分析,同时针对变载情况,在子程序中实现3种变载模型,能用于计算应力松弛行为。
四.附图说明:
图1:本发明所述方法流程图;
图2:典型的蠕变曲线;
图3:归一化参数蠕变模型各部分的意义;
图4:蠕变变形表达式中各系数的意义;
图5:直接时效GH4169G的蠕变曲线(680℃下不同应力);
图6:直接时效GH4169G的蠕变曲线(650MPa下不同温度);
图7:用表达式1对试验曲线的拟合结果(680℃下不同应力);
图8:用表达式1对试验曲线的拟合结果(650MPa下不同温度);
图9:用表达式2对试验曲线的拟合结果(680℃下不同应力);
图10:用表达式2对试验曲线的拟合结果(650MPa下不同温度);
图11:适用于变载情况的时间硬化模型;
图12:适用于变载情况的相对时间硬化模型;
图13:适用于变载情况的应变硬化模型;
图中符号说明如下:
η1,η2,η3为模型中参数,分别表示蠕变3个阶段的蠕变量;
ζ为无量纲时间,ζ=t/tc,tc为给定温度和应力下的持久寿命,ζ∈[0,1]。
五.具体实施方案
本发明一种描述蠕变变形的方法,它是一种能完整描述3个阶段蠕变变形的方法,所述的“3个阶段”是指:第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低,称为初始蠕变阶段;第2阶段蠕变速率不随时间变化,即稳态蠕变阶段;第3阶段蠕变速率随时间增大直至断裂,称为加速蠕变阶段,其具体步骤如下:
步骤一:将试验得到的不同温度和应力下的蠕变曲线归一化;
用试验得到的不同温度和应力下的蠕变曲线的时间坐标除以该温度和应力下的持久断裂时间,则所有曲线的横坐标均为归一化时间坐标ζ=t/tc,ζ∈[0,1]。
步骤二:对曲线的第2阶段进行拟合,得到η1,η2,η3;
对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程;得到的一组k值则为不同应力和温度下的一组η2值,即归一化坐标下的稳态蠕变应变率或最小蠕变率;拟合直线中的b值即η1,继而可得到η3=εr-η1-η2,εr为断裂时的蠕变应变。
步骤三:拟合得到η4,η5;
然后由已得到的参数η1,η2,η3,根据表达式或对每条曲线进行拟合(可采用matlab等软件进行最小二乘拟合),得到η4,η5。
步骤四:将ηi(i=1,2,3,4,5)表示为温度和应力的函数;
将ηi(i=1,2,3,4,5)表示为无量纲应力σ/σ0.2与无量纲温度T/Tm的函数。
步骤五:将该方法与有限元软件结合,编写usercreep子程序;
利用通用有限元软件ANSYS中usercreep子程序,将步骤一~步骤四得到的表达式写入子程序中,以达到计算实际结构蠕变变形的目的;对所编写的子程序进行编译连接后,在ANSYS主程序中调用子程序即可用所发明的方法对实际结构进行蠕变变形计算。
步骤六:选择适当的适用于变载情况下的蠕变模型;
在usercreep子程序中实现3种不同的适用于变载情况下蠕变模型:时间硬化模型,相对时间硬化模型,应变硬化模型。时间硬化模型为当载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下t时刻后的蠕变应变曲线上下平移得到;相对时间硬化模型为当载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下t2=t×tc,1/tc,2时刻后的蠕变曲线平移得到,tc,1和tc,2分别表示σ1,T1和σ2,T2状态下的持久寿命;应变硬化模型为载荷在t时刻变化(由σ1,T1变到σ2,T2),t时刻后的蠕变曲线由σ2,T2状态下产生与前一状态相同的蠕变应变所对应的时间之后曲线左右平移得到。
步骤七:建立实际结构的有限元模型,进行蠕变变形计算和应力松弛行为分析;
对实际结构在CAE前处理软件中建立有限元模型;由于应力松弛行为属于变载条件下的蠕变行为,能通过所编写的usercreep子程序进行计算分析;在ANSYS主程序中调用所编写的usercreep子程序,同时输入所求得模型参数值,对实际结构进行有限元数值模拟,计算蠕变变形和应力松弛行为。
其中,在步骤二中所述的“对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程”,其方法如下:利用第二阶段的试验数据点,可采用excel或matlab等软件对该数据点进行一次多项式拟合;若蠕变曲线的第2阶段并不明显,可采用多次多项式对蠕变曲线进行拟合,然后对多次多项式求导,导数值最小的点即为最小蠕变率点,导数值即η2。
其中,在步骤四中所述的“将ηi(i=1,2,3,4,5)表示为温度和应力的函数”,具体做法如下:选取特定函数将无量纲应力σ/σ0.2与无量纲温度T/Tm表示为ηi=f(σ/σ0.2,T/Tm),例如:
或
其中,ai,bi,ci,di(i=1,2,3,4,5)为材料相关系数,可由进行最小二乘拟合(可采用matlab等软件)得到;通过一组ai,bi,ci,di的系数值,则可求得任意温度和应力下的ηi,继而得到该温度和应力下的蠕变曲线。
其中,在步骤五中所述的“usercreep子程序”,其编写所需的输出变量为:蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1),蠕变应变增量对蠕变应变的导数dcrda(2);表达式2(即)的各输出变量为:
蠕变应变增量delcr:
蠕变应变增量对等效应力的导数dcrda(1):
式中,
蠕变应变增量对蠕变应变的导数dcrda(2):
其中,在步骤六中所述的“变载情况下的蠕变模型”,在子程序中实现方法为:
在子程序中,根据t2=t×tc,1/tc,2,可计算得到相对时间硬化模型中的t2;通过得到主程序返回当前时刻的应变值ε,根据表达式2(即),可迭代求得在σ2,T2状态下的蠕变曲线产生应变ε的对应的时间t2′;分别用t2和t2′替换usercreep子程序中输出变量中的t2,即可完成符合相对时间硬化模型和应变硬化模型的usercreep子程序。
下面将结合附图和实例对本发明作进一步的详细说明。图1为本发明所述方法的流程图。图2为典型的蠕变蠕变,蠕变3个阶段较为明显。第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低,称为初始蠕变阶段;第2阶段为直线,这是由于加工硬化与回复软化过程达到动态平衡而导致蠕变速率保持不变,即为稳态蠕变阶段;第3阶段蠕变速率随时间增大直至断裂,称为加速蠕变阶段。图3为所采用的蠕变应变表达式的各组成部分的曲线,说明其具有描述蠕变3个阶段的能力。
图4说明蠕变应变表达式各参数的物理意义,η1,η2,η3,分别为蠕变3个阶段的蠕变量,η4,η5分别控制蠕变第1阶段和第3阶段的变化快慢,且η5>1。步骤二拟合的直线即图中的斜线,根据各参数物理意义可得η1,η2,η3的值,从而拟合得到η4,η5。图5和图6为直接时效GH4169G的蠕变试验曲线,图7、图8和图9、图10分别为根据技术方案步骤一到步骤四,用表达式1和表达式2描述蠕变曲线得到的结果,可看出模型曲线和试验曲线的各部分均较为接近,说明该方法能够较好地完整描述3个阶段的蠕变曲线。
图11、图12和图13分别为变载下的时间硬化模型、相对时间硬化模型和应变硬化模型,分别由σ2,T2状态下的曲线通过平移得到。
Claims (5)
1.一种描述蠕变变形的方法,其特征在于:它是一种能完整描述3个阶段蠕变变形的方法,所述的“3个阶段”是指:第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低,称为初始蠕变阶段;第2阶段蠕变速率不随时间变化,即稳态蠕变阶段;第3阶段蠕变速率随时间增大直至断裂,称为加速蠕变阶段,其具体步骤如下:
步骤一:将试验得到的不同温度和应力下的蠕变曲线归一化;
用试验得到的不同温度和应力下的蠕变曲线的时间坐标除以该温度和应力下的持久断裂时间,则所有曲线的横坐标均为归一化时间坐标ζ=t/tc,ζ∈[0,1];ζ为无量纲时间,tc为给定温度和应力下的持久寿命;
步骤二:对曲线的第2阶段进行拟合,得到η1,η2,η3;η1、η2、η3分别为蠕变3个阶段的蠕变量;
对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程;得到的一组k值则为不同应力和温度下的一组η2值,即归一化坐标下的稳态蠕变应变率及最小蠕变率;拟合直线中的b值即η1,继而得到η3=εr-η1-η2,εr为断裂时的蠕变应变;
步骤三:拟合得到η4,η5;η4、η5分别控制蠕变第1阶段和第3阶段的变化快慢,且η5>1;
然后由已得到的参数η1,η2,η3,根据表达式蠕变应变量及对每条曲线进行拟合得到η4,η5;
步骤四:将ηi表示为温度和应力的函数;i=1,2,3,4,5;
将ηi表示为无量纲应力σ/σ0.2与无量纲温度T/Tm的函数;i=1,2,3,4,5;
步骤五:将该方法与有限元软件结合,编写usercreep子程序;
利用通用有限元软件ANSYS中usercreep子程序,将步骤一~步骤四得到的表达式写入子程序中,以达到计算实际结构蠕变变形的目的;对所编写的子程序进行编译连接后,在ANSYS主程序中调用子程序即能用所发明的方法对实际结构进行蠕变变形计算;
步骤六:选择预定的适用于变载情况下的蠕变模型;
在usercreep子程序中实现3种不同的适用于变载情况下蠕变模型:时间硬化模型,相对时间硬化模型,应变硬化模型;时间硬化模型为当载荷在t时刻变化,由σ1,T1变到σ2,T2,t时刻后的蠕变曲线由σ2,T2状态下t时刻后的蠕变应变曲线上下平移得到;相对时间硬化模型为当载荷在t时刻变化,由σ1,T1变到σ2,T2,t时刻后的蠕变曲线由σ2,T2状态下t2=t×tc,1/tc,2时刻后的蠕变曲线平移得到,tc,1和tc,2分别表示σ1,T1和σ2,T2状态下的持久寿命;应变硬化模型为载荷在t时刻变化,由σ1,T1变到σ2,T2,t时刻后的蠕变曲线由σ2,T2状态下产生与前一状态相同的蠕变应变所对应的时间之后曲线左右平移得到;
步骤七:建立实际结构的有限元模型,进行蠕变变形计算和应力松弛行为分析;
对实际结构在CAE前处理软件中建立有限元模型;由于应力松弛行为属于变载条件下的蠕变行为,能通过所编写的usercreep子程序进行计算分析;在ANSYS主程序中调用所编写的usercreep子程序,同时输入所求得模型参数值,对实际结构进行有限元数值模拟,计算蠕变变形和应力松弛行为。
2.根据权利要求1所述的一种描述蠕变变形的方法,其特征在于:在步骤二中所述的“对每条曲线第2阶段进行拟合,若曲线第2阶段较为明显,则得到直线方程y=kx+b,否则,根据最小蠕变率点得到直线方程”,其方法如下:利用第二阶段的试验数据点,采用excel及matlab软件对该数据点进行一次多项式拟合;若蠕变曲线的第2阶段并不明显,则采用多次多项式对蠕变曲线进行拟合,然后对多次多项式求导,导数值最小的点即为最小蠕变率点,导数值即η2。
3.根据权利要求1所述的一种描述蠕变变形的方法,其特征在于:在步骤四中所述的“将ηi表示为温度和应力的函数”,具体做法如下:选取特定函数将无量纲应力σ/σ0.2与无量纲温度T/Tm表示为ηi=f(σ/σ0.2,T/Tm),
及
其中,ai,bi,ci,di为材料相关系数,i=1,2,3,4,5;由进行最小二乘拟合得到;通过一组ai,bi,ci,di的系数值,则求得任意温度和应力下的ηi,继而得到该温度和应力下的蠕变曲线。
4.根据权利要求1所述的一种描述蠕变变形的方法,其特征在于:在步骤五中所述的“usercreep子程序”,其编写所需的输出变量为:蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1),蠕变应变增量对蠕变应变的导数dcrda(2);表达式的各输出变量为:
蠕变应变增量delcr:
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>r</mi>
<mo>=</mo>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>c</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<mi>&zeta;</mi>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<msup>
<mi>&zeta;</mi>
<mrow>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
蠕变应变增量对等效应力的导数dcrda(1):
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>r</mi>
<mi>d</mi>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>c</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&sigma;</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
</mfrac>
<mfrac>
<mn>1</mn>
<msub>
<mi>&sigma;</mi>
<mn>0.2</mn>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<mi>&zeta;</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>4</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<mi>&zeta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<msup>
<mi>&zeta;</mi>
<mrow>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>5</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>5</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mi>ln</mi>
<mi>&zeta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,ci,di为材料相关系数;
蠕变应变增量对蠕变应变的导数dcrda(2):
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>r</mi>
<mi>d</mi>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>c</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&epsiv;</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mfrac>
<mrow>
<msub>
<mi>d&epsiv;</mi>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>&eta;</mi>
<mn>4</mn>
<mn>2</mn>
</msubsup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<mi>&zeta;</mi>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>&zeta;</mi>
<mrow>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>4</mn>
</msub>
<mi>&zeta;</mi>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<msup>
<mi>&zeta;</mi>
<mrow>
<msub>
<mi>&eta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>.</mo>
</mrow>
5.根据权利要求1所述的一种描述蠕变变形的方法,其特征在于:在步骤六中所述的“变载情况下的蠕变模型”,在子程序中实现方法为:
在子程序中,根据t2=t×tc,1/tc,2,计算得到相对时间硬化模型中的t2;通过得到主程序返回当前时刻的应变值ε,根据表达式迭代求得在σ2,T2状态下的蠕变曲线产生应变ε的对应的时间t2′;分别用t2和t2′替换usercreep子程序中输出变量中的t2,即完成符合相对时间硬化模型和应变硬化模型的usercreep子程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510158190.4A CN104809273B (zh) | 2015-04-03 | 2015-04-03 | 一种描述蠕变变形的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510158190.4A CN104809273B (zh) | 2015-04-03 | 2015-04-03 | 一种描述蠕变变形的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104809273A CN104809273A (zh) | 2015-07-29 |
CN104809273B true CN104809273B (zh) | 2017-12-15 |
Family
ID=53694090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510158190.4A Expired - Fee Related CN104809273B (zh) | 2015-04-03 | 2015-04-03 | 一种描述蠕变变形的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104809273B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529017A (zh) * | 2016-11-08 | 2017-03-22 | 华南理工大学 | 一种高铬钢构件高温蠕变变形预测与蠕变损伤分析方法 |
CN106777479B (zh) * | 2016-11-18 | 2020-04-21 | 中国航空动力机械研究所 | 基于梁理论的涡轮叶片非线性蠕变分析方法 |
CN109359330B (zh) * | 2018-09-05 | 2021-02-26 | 南京航空航天大学 | 一种描述高温材料蠕变变形三阶段的方法及模型 |
CN109142083B (zh) * | 2018-09-05 | 2020-04-24 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法 |
CN110967245B (zh) * | 2018-09-28 | 2023-01-20 | 中国航发商用航空发动机有限责任公司 | 材料蠕变时间及寿命实验方法以及实验系统 |
CN109117604B (zh) * | 2018-10-25 | 2020-03-31 | 北京航空航天大学 | 一种镍基单晶合金基于滑移系的Walker模型材料常数的获取方法 |
CN110688788B (zh) * | 2019-08-28 | 2021-06-22 | 南京航空航天大学 | 一种高温材料蠕变变形和寿命预测方法 |
CN112730092B (zh) | 2020-12-24 | 2022-07-12 | 南京航空航天大学 | 一种基于失效模式一致的蠕变载荷等效加速方法 |
CN114493012B (zh) * | 2022-01-28 | 2023-08-08 | 武汉科技大学 | 考虑硬化损伤协同效应的胶结充填体蠕变失稳预测方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003167929A (ja) * | 2001-09-18 | 2003-06-13 | Fujitsu Nagano Systems Engineering Ltd | 構造解析プログラム、構造解析方法、構造解析装置および半導体集積回路の製造方法 |
CN102607951A (zh) * | 2012-03-08 | 2012-07-25 | 北京航空航天大学 | 一种无附加弯矩双头榫连拉伸试验结构 |
CN103324805A (zh) * | 2013-07-04 | 2013-09-25 | 沈阳华创风能有限公司 | 风机底架极限强度校核的计算方法 |
CN104361167A (zh) * | 2014-11-04 | 2015-02-18 | 南京航空航天大学 | 一种基于相场法分析的含电极铁电单晶有限元预测方法 |
-
2015
- 2015-04-03 CN CN201510158190.4A patent/CN104809273B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003167929A (ja) * | 2001-09-18 | 2003-06-13 | Fujitsu Nagano Systems Engineering Ltd | 構造解析プログラム、構造解析方法、構造解析装置および半導体集積回路の製造方法 |
CN102607951A (zh) * | 2012-03-08 | 2012-07-25 | 北京航空航天大学 | 一种无附加弯矩双头榫连拉伸试验结构 |
CN103324805A (zh) * | 2013-07-04 | 2013-09-25 | 沈阳华创风能有限公司 | 风机底架极限强度校核的计算方法 |
CN104361167A (zh) * | 2014-11-04 | 2015-02-18 | 南京航空航天大学 | 一种基于相场法分析的含电极铁电单晶有限元预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104809273A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104809273B (zh) | 一种描述蠕变变形的方法 | |
CN103745114B (zh) | 一种钛合金应力松弛数值及回弹计算方法 | |
CN105628511B (zh) | 一种高温合金蠕变寿命预测方法 | |
US9283695B1 (en) | Computer-implemented simulation method and non-transitory computer medium capable of predicting fiber orientation for use in a molding process | |
CN103886125A (zh) | 一种钛合金热复合成形数值模拟方法 | |
CN109470553B (zh) | 应用Laplace变换实现沥青混合料材料参数快速转换的方法 | |
CN110532727B (zh) | 可用于常见非牛顿流体的数值模拟方法 | |
CN106525612B (zh) | 基于拉伸及压缩实验的聚碳酸酯统一本构模型的构建方法 | |
CN104316388A (zh) | 一种对各向异性材料结构件进行疲劳寿命测定的方法 | |
CN104182585B (zh) | 硅胶弹性体手感分析及疲劳寿命预测的有限元模拟方法 | |
CN104615840B (zh) | 一种数字仿真模型的修正方法及系统 | |
CN105737881B (zh) | 一种电子设备加速可靠性增长试验方法 | |
CN111062162B (zh) | 一种岩土材料精确本构模型的数值建模与应用方法 | |
WO2013042600A1 (ja) | Chabocheモデルを用いた応力-ひずみ関係シミュレーション方法、応力-ひずみ関係シミュレーションシステム、応力-ひずみ関係シミュレーションプログラム | |
CN104614995A (zh) | 一种二阶系统有限时间滑模控制器的通用设计方法 | |
CN104239625A (zh) | 一种基于修正流体运动方程线性迭代的稳态求解方法 | |
CN109271655B (zh) | 一种基于非对称有限元算法的材料尺度效应分析方法 | |
CN105844025B (zh) | 一种针对高超声速舵面的非概率热气动弹性可靠性设计方法 | |
CN105046030A (zh) | 基于有限元法的三维传热条件下的铝合金构件淬火过程换热系数的获得方法 | |
CN109117604A (zh) | 一种镍基单晶合金基于滑移系的Walker模型材料常数的获取方法 | |
CN106682284A (zh) | 复合材料构件热校形工艺的仿真模拟方法 | |
CN115034121B (zh) | 一种基于组织性能智能预报模型的带钢工艺调控方法 | |
CN115310311B (zh) | 一种金属材料在谱载荷下疲劳寿命计算方法 | |
CN110955999B (zh) | 轮胎三维温度场仿真方法 | |
CN109558669A (zh) | 一种基于有限元模型的汽轮机转子疲劳损伤在线计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171215 Termination date: 20190403 |
|
CF01 | Termination of patent right due to non-payment of annual fee |