一种描述蠕变变形的方法技术

一种描述蠕变变形的方法，它是一种能完整描述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种变载模型。

【技术实现步骤摘要】

本专利技术提供一种能完整描述3个阶段蠕变变形的模型及将其应用于实际结构的 蠕变计算的方法，属于高温结构蠕变模拟
二.
技术介绍
蠕变是指在温度、载荷不变的条件下，材料的变形也会随时间增长而逐渐增大的 现象，而且这种变形即使在应力小于屈服极限时仍具有不可逆的变形性质。典型的蠕变变 形可以分为3个阶段：第1阶段由于变形引起加工硬化导致蠕变速率随时间不断降低，称为 初始蠕变阶段；第2阶段为直线，这是由于加工硬化与回复软化过程达到动态平衡而导致 蠕变速率保持不变，即为稳态蠕变阶段；第3阶段蠕变速率随时间增大直至断裂，称为加速 懦变阶段。 当前工程中常用的蠕变模型有时间硬化模型和应变硬化模型等，但这类模型只能 模拟蠕变第1阶段或者前两个阶段（有限元软件ANSYS和ABAQUS提供的蠕变模型均不能 模拟蠕变第3阶段）；1985年Evans和Wilshire在其著作中提出一种能模拟完整的蠕变变 形的模型一 0投射法，虽然0投射法可以模拟蠕变变形的3个阶段，但它并不注重描述稳 态蠕变阶段，即它所描述的整个蠕变变形中，并不存在蠕变速率为常数的阶段；一些粘塑性 本构模型虽能较好地描述循环载荷下的塑形变形，而在描述蠕变第3阶段方面必须引入损 伤参量对模型进行修正。这些模型用于试样在试验载荷下变形的描述较精准，但由于计算 较为复杂，目前还不易用于实际结构（如涡轮叶片和轮盘）的蠕变分析。 因此，仍有必要发展能够完整描述蠕变曲线的3个阶段的方法，同时能与有限元 结合用于实际结构的蠕变变形计算及应力松弛问题的分析。 三.
技术实现思路
1?专利技术目的 本专利技术提供一种能完整描述3个阶段蠕变变形的模型及将其应用于实际结构的 蠕变计算的方法，达到能计算实际结构3个阶段的蠕变变形和模拟应力松弛行为的目的， 解决现有技术的不足。 2?技术方案 本专利技术提供，它能完整描述3个阶段（第1阶段由于 变形引起加工硬化导致蠕变速率随时间不断降低，称为初始蠕变阶段；第2阶段蠕变速 率不随时间变化，即稳态蠕变阶段；第3阶段蠕变速率随时间增大直至断裂，称为加速蠕 变阶段）懦变变形的方法，同时通过将该方法与有限元软件ANSYS结合，编写usercreep 子程序，用于实际结构的蠕变变形计算。同时提出3种适用于变载情况的模型，用于 计算应力松弛行为。本专利技术用归一化参数描述蠕变变形，将蠕变应变表示为：【主权项】1. ，其特征在于：它是一种能完整描述3个阶段蠕变变形的 方法，所述的"3个阶段"是指：第1阶段由于变形引起加工硬化导致蠕变速率随时间不断 降低，称为初始蠕变阶段；第2阶段蠕变速率不随时间变化，即稳态蠕变阶段；第3阶段蠕 变速率随时间增大直至断裂，称为加速蠕变阶段，其具体步骤如下： 步骤一：将试验得到的不同温度和应力下的蠕变曲线归一化； 用试验得到的不同温度和应力下的蠕变曲线的时间坐标除以该温度和应力下的持久 断裂时间，则所有曲线的横坐标均为归一化时间坐标G=t/t。，GG; 步骤二：对曲线的第2阶段进行拟合，得到1，n2,n3; 对每条曲线第2阶段进行拟合，若曲线第2阶段较为明显，则得到直线方程y=kx+b， 否则，根据最小蠕变率点得到直线方程；得到的一组k值则为不同应力和温度下的一组n2 值，即归一化坐标下的稳态蠕变应变率及最小蠕变率；拟合直线中的b值即ni，继而得到 n3=ern「n2, \为断裂时的蠕变应变； 步骤三：拟合得到n4,n5; 然后由已得到的参数ru，n2,n3，根据表达式及 & -e对每条曲线进行拟合得到n4，n5; 步骤四：将n3亥i= 1，2, 3, 4, 5,表示为温度和应力的函数； 将ni该i= 1，2, 3, 4, 5表示为无量纲应力0 / 0a2与无量纲温度T/Tj^函数； 步骤五：将该方法与有限元软件结合，编写usercreep子程序； 利用通用有限元软件ANSYS中usercreep子程序，将步骤一~步骤四得到的表达式写 入子程序中，以达到计算实际结构蠕变变形的目的；对所编写的子程序进行编译连接后，在 ANSYS主程序中调用子程序即能用所专利技术的方法对实际结构进行蠕变变形计算； 步骤六：选择预定的适用于变载情况下的蠕变模型； 在usercreep子程序中实现3种不同的适用于变载情况下懦变模型：时间硬化模型， 相对时间硬化模型，应变硬化模型；时间硬化模型为当载荷在t时刻变化（由〇u变到 〇 2,T2)，t时刻后的蠕变曲线由〇 2,T2状态下t时刻后的蠕变应变曲线上下平移得到；相对 时间硬化模型为当载荷在t时刻变化，由〇1，1\变到〇 2，T2,t时刻后的蠕变曲线由〇2，T2 状态下心二七乂七^/一时刻后的蠕变曲线平移得到^^和七"分别表示^上和％，!^ 状态下的持久寿命；应变硬化模型为载荷在t时刻变化，由〇uT/变到〇 2,T2,t时刻后的 蠕变曲线由〇2，T2状态下产生与前一状态相同的蠕变应变所对应的时间之后曲线左右平 移得到； 步骤七：建立实际结构的有限元模型，进行蠕变变形计算和应力松弛行为分析； 对实际结构在CAE前处理软件中建立有限元模型；由于应力松弛行为属于变载条件下 的蠕变行为，能通过所编写的usercre印子程序进行计算分析；在ANSYS主程序中调用所编 写的usercreep子程序，同时输入所求得模型参数值，对实际结构进行有限元数值模拟，计 算蠕变变形和应力松弛行为。2. 根据权利要求1所述的，其特征在于：在步骤二中所述 的"对每条曲线第2阶段进行拟合，若曲线第2阶段较为明显，则得到直线方程y=kx+b， 否则，根据最小蠕变率点得到直线方程"，其方法如下：利用第二阶段的试验数据点，采用 excel及matlab软件对该数据点进行一次多项式拟合；若蠕变曲线的第2阶段并不明显， 则采用多次多项式对蠕变曲线进行拟合，然后对多次多项式求导，导数值最小的点即为最 小蠕变率点，导数值即n2。3. 根据权利要求1所述的，其特征在于：在步骤四中所述的 "将n3亥i= 1，2, 3, 4, 5表示为温度和应力的函数"，具体做法如下：选取特定函数将无量 纲应力〇/〇Q2与无量纲温度T/Tm表示为ni=f(〇/〇Q.2，T/Tm)，例如：其中，aybyq，屯该i= 1，2, 3, 4, 5为材料相关系数，由进行最小二乘拟合得到；通过 一组叫七，(^，屯的系数值，则求得任意温度和应力下的1^，继而得到该温度和应力下的蠕 变曲线。4. 根据权利要求1所述的，其特征在于：在步骤五中所述 的"usercreep子程序"，其编写所需的输出变量为：蠕变应变增量delcr、蠕变应变增量 对等效应力的导数dcrda(l)，蠕变应变增量对蠕变应变的导数dcrda(2);表达式2即 A，(l-e + 的各输出变量为： 蠕变应变增量delcr:蠕变应变增量对等效应力的导数dcrda(l):式中，乐=f('.=.丨，2A4,5); 上m 蠕变应变增量对蠕变应变的导数dcrda(2):5. 根据权利要求1所述的，其特征在于：在步骤六中所述的 "变载情况下的蠕变模型"，在子程序中实现方法为： 在子程序中，根据t2 =tXt。，Vt。，2，计算得到相对时间硬化模型中的t2;通过得到主本文档来自技高网...

【技术保护点】
一种描述蠕变变形的方法，其特征在于：它是一种能完整描述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，根据表达式及对每条曲线进行拟合得到η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子程序，同时输入所求得模型参数值，对实际结构进行有限元数值模拟，计算蠕变变形和应力松弛行为。...

【技术特征摘要】

【专利技术属性】
技术研发人员：王延荣，程域钊，刘香，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11