详情

原文传递塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法

专利名称：	塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法
摘要：	本发明公开了一种塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，在Davies工作的基础上，提出了耦合残余应力和拘束效应的蠕变孕育期预测模型。利用韧性耗散损伤模型，引入了与载荷无关拘束参量Q*计算考虑拘束效应的蠕变孕育期。使用紧凑拉伸试样(CT)施加主载荷进行蠕变模拟实验，本发明的有益效果：能够在结构中简洁有效的预测出塑性条件下蠕变孕育期。
专利类型：	发明专利
国家地区组织代码：	天津;12
申请人：	天津大学
发明人：	徐连勇;邬栋权;荆洪阳;韩永典;赵雷;吕小青
专利状态：	有效
申请日期：	2017-12-15T00:00:00+0800
发布日期：	2019-06-25T00:00:00+0800
申请号：	CN201711352974.6
公开号：	CN109932241A
代理机构：	天津创智天诚知识产权代理事务所(普通合伙)
代理人：	王秀奎
分类号：	G01N3/08(2006.01);G;G01;G01N;G01N3
申请人地址：	300072 天津市南开区卫津路92号
主权项：	1.塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，其特征在于，包括以下步骤： S1：建立预测模型，所述预测模型包括CT试样本体，CT试样本体的中部前端设有槽，槽的后部设有缺口，槽、缺口在同一平面上，CT试样本体上还设有上主载荷销孔、下主载荷销孔，上主载荷销孔、下主载荷销孔上下对称设置，分别设置在槽的上下两端， S2：首先利用上圆销、下圆销对CT试样本体的上下两端进行预定大小的压缩加载，然后释放上圆销、下圆销，会在CT试样本体的缺口附近产生残余应力分布； S3：在含有残余应力的缺口处插入预制裂纹，以进行蠕变试验； S4：利用销子在上主载荷销孔、下主载荷销孔施加主载荷，进行高温蠕变试验； S5：通过蠕变有限元模拟可以获得计算含残余应力、拘束参量CT试样孕育期所需要的必要参数，在塑性条件下，计算孕育期主要包括以下步骤： (1)首先计算复合加载下的应力强度因子，其计算公式为： (Ⅰ)中：其中：是模拟计算的只含有残余应力下的应力强度因子，单位为MPa·(m1/2)；是主载荷应力强度因子，单位为MPa·(m1/2)；P是主载荷，单位为N；B是试样厚度，单位为mm，Bn是试样净厚度，单位为mm；a/W是预制裂纹长度比率，a是预制裂纹长度，采用上主载荷销孔圆心到预制裂纹后端的水平直线距离，单位为mm；W是名义试样宽度，采用上主载荷销孔圆心到CT试样本体后端的水平直线距离，单位为mm；f(a/W)是CT试样几何系数，只与a/W有关；V是无量纲的塑性相关项，计算如下： (Ⅱ)中:V0是无量纲参量，其中：是塑性残余应力强度因子，单位为MPa·(m1/2)；是弹性残余应力强度因子，单位为MPa·(m1/2)，利用JS计算，JS是残余应力场下断裂参量，单位为MPa·m；其中：E′是有效弹性模量：E'＝E/(1-ν2)，E是弹性模量，ν是泊松比，和JS都利用有限元模拟结果提取； (Ⅱ)中:Lr是无量纲参量，描述主载荷幅度：其中：σy是屈服强度，与材料属性有关，单位为MPa；σrefP是主载荷参考应力，单位为MPa，用下式计算：其中：nL为无量纲裂纹深宽比参数，通过下式计算: 常数 (Ⅱ)中: 是弹性主载荷应力强度因子，是塑性主载荷应力强度因子，单位为MPa·(m1/2)；利用有限元模拟结果计算： (Ⅱ)中:β描述残余应力的幅度，是无量纲参量；是二次载荷参考应力，利用有限元模拟计算， (Ⅱ)中:Z为无量纲的弹性追随因子，从有限元模拟结果中提取出应力应变关系，取等效蠕变应变增量与等效弹性应变增量的比值： (2)计算塑性复合应力场下J积分数值，其计算公式为：其中：KI是复合应力强度因子，E′是有效弹性模量：E'＝E/(1-ν2)，ν是泊松比，E是弹性模量，ν是泊松比； (3)然后计算塑性条件下的拘束参量Q*HRR，其计算公式为：是利用有限元计算得出的裂纹前沿处的张开应力值，单位是Mpa，σP0是标准化应力，单位为MPa，εP0是标准化应变，单位为1，α为应变硬化系数，N为应变硬化指数，IN是与N有关的无量纲函数，L是标量距离，取1mm； (Ⅲ)中：σ22是利用HRR应力场计算得出的裂纹前沿的张开应力值，单位是MPa，其中：r是裂纹后部尖端到裂纹前沿研究点的间距，单位是mm，θ是裂纹尖端角度，是与θ和N有关的无量纲函数， (4)计算等效应力单位是MPa，其计算公式为：其中：是与θ和N有关的无量纲函数； (5)然后计算塑性应力场下孕育期时间tiHRR，其计算公式为： (Ⅴ)中：n是无量纲的蠕变应力硬化指数，εcrit是单轴蠕变韧性，与材料属性有关，单位为1，是蠕变应变变化率，单位为h-1，与材料高温蠕变属性有关， (Ⅴ)中：MSFHRR为塑性条件下多轴应力因子，根据Cocks and Ashby关系式计算： sinh是双曲正弦函数，hHRR为塑性应力三轴度，在塑性应力状态下：其中平均应力单位是MPa，其计算公式为：其中：σ11,σ22和σ33是利用HRR应力场计算得出的裂纹前沿的应力值，单位是MPa，其中：是与θ和N有关的无量纲函数。 2.如权利要求1所述的塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，其特征在于，裂纹后部尖端到裂纹前沿研究点的间距r取d，d是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离，即蠕变萌生发生的临界距离。 3.如权利要求2所述的塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，其特征在于，d取所研究材料的晶粒尺寸。 4.如权利要求1所述的塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，其特征在于，Bn＝B。 5.如权利要求1所述的塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法，其特征在于，所述有限元模拟采用ABAQUS6.14进行计算模拟，JS、的提取过程包括以下步骤： (1)首先建立预压缩加载的CT试样的有限元模型，在材料属性模块里设置弹性塑性参数，在载荷模块里设置压缩载荷，以及拘束条件：包括对称条件和固定条件，在接触模块里设置压缩圆销与试样上下表面的的刚性接触，在分析步模块里设置好输出参量：应力值，在网格模块划分网格； (2)在作业模块提交任务计算，获得残余应力的计算结果，结果文件中，从场变量中可以直接提取二次载荷参考应力 (3)建立相同尺寸的试样模型，进行主载荷拉伸试验，在材料属性模块里设置高温下的弹性塑性蠕变参数，在网格模块划分网格，在接触模块里设置拉伸销子与销孔的刚性接触，并在模型中插入预制裂纹，在分析步模块里设置好输出参量：应力应变值、应力强度因子K值，在载荷模块里设置拉伸载荷，以及拘束条件：包括对称条件和固定条件，在预加载应力场中导入上步算好的残余应力； (4)在作业模块提交任务计算，获得含有残余应力的蠕变拉伸实验计算结果，结果文件中，在插入裂纹后还未施加拉伸载荷时刻，在场变量可以获取应力值从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量JS，在施加拉伸载荷的初始时刻，可以获取塑性主载荷强度因子从历史变量中可以获取等效应力随总应变增量的变化曲线，从曲线中得到等效蠕变应变增量，等效弹性应变增量进而得到弹性追随因子Z计算方法。
所属类别：	发明专利