本发明涉及耦合残余应力和拘束效应高温结构在塑性瞬态蠕变条件下的蠕变孕育期工程临界评价，就是确定当结构中存在表面裂纹且处在塑性瞬态蠕变条件下时，评价此高温结构的蠕变裂纹萌生寿命。
背景技术：
：：燃煤为主的能源结构是我国雾霾天气的主要成因之一，而燃煤发电是我国目前最主要的发电方式，该趋势将长期存在。因此，除了改变能源结构，发展高效洁净的超超临界(usc)机组是节能减排的重要途径之一。然而，蒸汽温度和压力等参数的提高导致机组关键高温管道的服役环境非常恶劣，特别是管道中存在裂纹、未焊透、焊接气孔和夹渣等各种缺陷，严重威胁机组的安全运行，需要对其进行科学精确的寿命评估。几十年来，对于高温下含裂纹构件，国外发展了多种高温蠕变寿命的评定规范和方法。蠕变孕育期是蠕变过程中经历时间最长的阶段，孕育期的准确预测对于高温结构的蠕变寿命预测具有重要意义；davies等人基于韧性耗散模型提出的孕育期预测模型，考虑了蠕变过程应力变化的完整性，但是结构的残余应力以及拘束效应对孕育期的影响并未得到研究；残余应力、拘束效应广泛的存在于加工制造的高温部件中，并且对部件的服役寿命造成了重大影响。大量针对于高温蠕变情况下残余应力和拘束效应的研究也广泛展开。因此建立耦合残余应力和拘束效应的蠕变孕育期预测模型，可以更加准确完整地评估复合加载结构的蠕变孕育期。技术实现要素：本发明的目的是针对现有技术中存在的技术缺陷，而提供一种发明名称。为实现本发明的目的所采用的技术方案是：本发明在davies工作的基础上，提出了耦合残余应力和拘束效应的蠕变孕育期预测模型。利用参考应立法，引入了弹性追随因子z计算考虑残余应力的蠕变孕育期。使用紧凑拉伸试样(ct)通过预压缩产生残余应力，并施加主载荷进行蠕变实验。为实现本发明的目的所采用的技术方案是：本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,包括以下步骤：s1：建立模型，包括ct试样本体，ct试样本体的中部前端设有槽，槽的后部设有缺口，ct试样本体上还设有上主载荷销孔、下主载荷销孔，上主载荷销孔、下主载荷销孔上下对称设置，分别设置在槽的上下两端。s2：首先利用上圆销、下圆销对ct试样本体进行预定大小的压缩加载，然后释放上圆销、下圆销，会在ct试样本体的缺口附近产生残余应力分布；s3：在含有残余应力的缺口处插入预制裂纹，以进行蠕变试验；s4：利用销子在上主载荷销孔、下主载荷销孔施加主载荷，进行高温蠕变试验；s5：通过蠕变有限元模拟可以获得计算含残余应力ct试样孕育期所需要的必要参数，在塑性瞬态蠕变条件下，如图4所示，研究点的初始应力是塑性应力状态，到达转换时间thrr-rr后进入到瞬态蠕变应力状态，计算孕育期主要包括以下步骤：(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是弹性模量，ν是泊松比，e和ν二者均参见文献：(zhaol,jingh,xul,hany,xiuj.evaluationofconstrainteffectsoncreepcrackgrowthbyexperimentalinvestigationandnumericalsimulation.engngfractmech2012；96:251–66.)，和js都利用有限元模拟结果提取；(ⅱ)中:lr是无量纲参量，描述主载荷幅度：其中：σy是屈服强度，单位为mpa，参见文献：(zhaol,jingh,xul,hany,xiuj.evaluationofconstrainteffectsoncreepcrackgrowthbyexperimentalinvestigationandnumericalsimulation.engngfractmech2012；96:251–66.)；是主载荷参考应力，单位为mpa，用下式计算：其中：nl为无量纲裂纹深宽比参数，通过下式计算:常数(ⅱ)中:其中：是弹性主载荷应力强度因子，单位为mpa·(m1/2)，是塑性主载荷应力强度因子，单位为mpa·(m1/2)；利用有限元模拟结果计算：(ⅱ)中:β描述残余应力的幅度，是无量纲参量；是二次载荷参考应力，单位为mpa，利用有限元模拟计算；(ⅱ)中:z为无量纲的弹性追随因子，从有限元模拟结果中提取出应力应变关系，取等效蠕变应变增量与等效弹性应变增量的比值：(2)计算稳态蠕变复合应力场下c*积分数值，其计算公式为：其中:a是蠕变硬化系数，单位为mpa-n·h-1，a参见文献：(zhaol,jingh,xul,hany,xiuj.evaluationofconstrainteffectsoncreepcrackgrowthbyexperimentalinvestigationandnumericalsimulation.engngfractmech2012；96:251–66.)，ki是复合应力强度因子，单位为mpa·(m1/2)，是初始参考应力，单位为mpa；(3)然后计算裂尖参量c(t)，它是反映瞬态蠕变过程的载荷环路积分值，单位为mpa·mm·(h-1)，利用参考应力法计算：(ⅴ)中：σref是总参考应力，单位为mpa，利用以下积分式计算：其中：是总参考应变速率，单位为h-1，是主载荷参考应变速率，单位为h-1，(ⅴ)中：εref是总参考应变，利用下式计算：εref＝ε0ref+a∫σnrefdt其中：ε0ref是初始参考应变，ε0ref通过有限元模拟提取。(4)计算瞬态蠕变条件下的拘束参量q*rr，其计算公式为：是利用有限元计算得出的裂纹前沿处的张开应力值，单位是mpa，σ0是材料的屈服强度，单位是mpa，参见文献：(zhaol,xul,hany,jingh.two-parametercharacterizationofconstrainteffectinducedbyspecimensizeoncreepcrackgrowth.engngfractmech2012；96:251–66.)，l是标量距离，取1mm；(ⅵ)中：σ22是利用hrr应力场计算得出的裂纹前沿的张开应力值，单位是mpa，其中：r是裂纹后部尖端到裂纹前沿研究点的间距，单位是mm，θ是裂纹尖端角度，是蠕变应变变化率，单位为h-1，与材料高温蠕变属性有关，n为无量纲的蠕变应力硬化指数，n和参见文献：(zhaol,jingh,xul,hany,xiuj.evaluationofconstrainteffectsoncreepcrackgrowthbyexperimentalinvestigationandnumericalsimulation.engngfractmech2012；96:251–66.)，in是与n有关的无量纲函数，是与θ和n有关的无量纲函数，in和具体值可以查阅文献得到：shih,c.f..1983.tablesofhutchinson-rice-rosengrensingularfieldquantities.brownuniversitytechnicalreport,mrle-147.(5)计算瞬态蠕变等效应力其计算公式为：其中：是与θ和n有关的无量纲函数，具体值可以查阅文献得到：shih,c.f..1983.tablesofhutchinson-rice-rosengrensingularfieldquantities.brownuniversitytechnicalreport,mrle-147.计算塑性等效应力其计算公式为：σp0是标准化应力，单位为mpa，εp0是标准化应变，单位为1，α为应变硬化系数，n为应变硬化指数，σp0，εp0，α和n参见文献：(zhaol,xul,hany,jingh.two-parametercharacterizationofconstrainteffectinducedbyspecimensizeoncreepcrackgrowth.engngfractmech2012；96:251–66.)，in是与n有关的无量纲函数,是与裂纹尖端角度θ和应变硬化指数n有关的无量纲函数，in和可查表获得shih,c.f..1983.tablesofhutchinson-rice-rosengrensingularfieldquantities.brownuniversitytechnicalreport,mrle-147.(6)利用matalab软件计算转换时间thrr-rr：在此时刻满足：塑性阶段损伤累计值：msfhrr为塑性条件下多轴应力因子，根据cocksandashby关系式计算：sinh是双曲正弦函数，hhrr为塑性应力三轴度，在塑性应力状态下：其中：和是与θ和n有关的无量纲函数，具体值可以查阅文献得到：shih,c.f..1983.tablesofhutchinson-rice-rosengrensingularfieldquantities.brownuniversitytechnicalreport,mrle-147.。(7)然后计算瞬态蠕变应力场下孕育期时间tihrr-rr，其计算公式为：(ⅷ)中：d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离，即蠕变萌生发生的临界距离。(ⅷ)中：msfrr为塑性条件下多轴应力因子，根据cocksandashby关系式计算：sinh是双曲正弦函数，hrr为瞬态蠕变应力三轴度，在塑性应力状态下：其中平均应力单位是mpa，其计算公式为：其中：σ11和σ33是利用rrss应力场计算得出的裂纹前沿的应力值，单位是mpa，其中：是与θ和n有关的无量纲函数，具体值可以查阅文献得到：shih,c.f..1983.tablesofhutchinson-rice-rosengrensingularfieldquantities.brownuniversitytechnicalreport,mrle-147.。优选的，d取所研究材料的晶粒尺寸。优选的，bn＝b。优选的，所述有限元模拟采用abaqus6.14进行计算模拟，js、ε0ref的提取过程包括以下步骤：(5)首先建立预压缩加载的ct试样的有限元模型，在材料属性模块里设置弹性塑性参数，在载荷模块里设置压缩载荷，以及拘束条件：包括对称条件和固定条件。在接触模块里设置压缩圆销与试样上下表面的的刚性接触，在分析步模块里设置好输出参量：应力值，在网格模块划分网格；(6)在作业模块提交任务计算，获得残余应力的计算结果，结果文件中，从场变量中可以直接提取二次载荷参考应力σrefs；(7)建立相同尺寸的试样模型，进行主载荷拉伸试验，在材料属性模块里设置高温下的弹性塑性蠕变参数，在网格模块划分网格，在接触模块里设置拉伸销子与销孔的刚性接触，并在模型中插入预制裂纹，在分析步模块里设置好输出参量：应力应变值、应力强度因子k值、断裂参量j积分值，在载荷模块里设置拉伸载荷，以及拘束条件：包括对称条件和固定条件，在预加载应力场中导入上步算好的残余应力；(8)在作业模块提交任务计算，获得含有残余应力的蠕变拉伸实验计算结果，结果文件中，在插入裂纹后还未施加拉伸载荷时刻，从场变量中可以获取初始参考应变从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量js，在施加拉伸载荷的初始时刻，可以获取塑性主载荷强度因子从历史变量中可以获取等效应力随总应变增量的变化曲线，从曲线中得到等效蠕变应变增量等效弹性应变增量进而得到弹性追随因子z计算方法。与现有技术相比，本发明的有益效果是：本发明提出了耦合残余应力和拘束效应的弹性瞬态蠕变条件下蠕变孕育期预测模型，与现有模型相比，本设计方法能够将原有的预测模型扩展到含残余应力的模型中，从而提出一种简化的瞬态蠕变条件下蠕变孕育期预测方法，因此能够在结构中简介有效的预测出弹性瞬态蠕变条件下蠕变孕育期。附图说明图1紧凑拉伸试样(ct)预压缩示意图；其中：1-上圆销，2-ct试样本体，3-上主载荷销孔，4-槽，5-缺口，6-预制裂纹，7-下主载荷销孔，8-下圆销。图2是蠕变裂纹萌生临界条件示意图。图3是应力应变关系曲线。图4初始应力向瞬态蠕变转换时间。具体实施方式以下结合附图和具体实施例对本发明作进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测模型，包括ct试样本体2，ct试样本体的中部前端设有槽4，槽的后部设有缺口5，槽4、缺口5在同一平面上，ct试样本体上还设有上主载荷销孔3、下主载荷销孔7，上主载荷销孔、下主载荷销孔上下对称设置，分别设置在槽的上下两端。选取p92高温耐热钢，以b＝20mm，w＝40mm，a/w＝0.5的ct试样作为研究对象，以预加载为12000n和主载荷p＝12000n作为研究载荷。其主要材料属性见下表：其中e-16为10的-16次方。本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,包括以下步骤：s1：首先利用上圆销1、下圆销8对ct试样本体进行预定大小的压缩加载，然后释放上圆销、下圆销，会在ct试样本体的缺口5附近产生残余应力分布；s2：在含有残余应力的缺口5处插入预制裂纹6，以进行蠕变试验；s3：利用销子在上主载荷销孔3、下主载荷销孔7施加主载荷，进行高温蠕变试验；s4：通过蠕变有限元模拟可以获得计算含残余应力ct试样孕育期所需要的必要参数，在塑性瞬态蠕变条件下，如图4所示，研究点的初始应力是塑性应力状态，到达转换时间thrr-rr后进入到瞬态蠕变应力状态。所述有限元模拟采用abaqus6.14进行计算模拟，js、ε0ref的提取过程包括以下步骤：(ⅰ)首先建立如图1所示的预压缩加载的ct试样的有限元模型，在材料属性模块里设置弹性塑性参数，在载荷模块里设置压缩载荷，以及拘束条件：包括对称条件和固定条件。在接触模块里设置压缩圆销与试样上下表面的的刚性接触，在分析步模块里设置好输出参量：应力值，在网格模块划分网格；(ⅱ)在作业模块提交任务计算，获得残余应力的计算结果，结果文件中，从场变量中可以直接提取二次载荷参考应力(ⅲ)建立如图1所示的试样模型，进行主载荷拉伸试验，在材料属性模块里设置高温下的弹性塑性蠕变参数，在网格模块划分网格，在接触模块里设置拉伸销子与销孔的刚性接触，并在模型中插入预制裂纹，在分析步模块里设置好输出参量：应力应变值、应力强度因子k值、断裂参量j积分值，在载荷模块里设置拉伸载荷，以及拘束条件：包括对称条件和固定条件，在预加载应力场中导入上步算好的残余应力；(ⅳ)在作业模块提交任务计算，获得含有残余应力的蠕变拉伸实验计算结果，结果文件中，在插入裂纹后还未施加拉伸载荷时刻，在插入裂纹后还未施加拉伸载荷时刻，从场变量中可以获取初始参考应变从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量js＝0.013mpa·m，可以计算得塑性残余应力强度因子：在施加拉伸载荷的初始时刻，可以获取塑性主载荷强度因子从历史变量提取出c(t)随时间的变化曲线，如图4所示，从历史变量中可以获取等效应力随总应变增量的变化曲线，如图3所示，从曲线中得到等效蠕变应变增量等效弹性应变增量进而得到弹性追随因子z计算方法。计算孕育期主要包括以下步骤：(1)首先计算各个参量：(a)弹性主载荷强度因子:(b)主载荷参考应力：(c)主载荷幅度:(d)残余应力参考应力：残余应力的幅度:(e)弹性追随因子：(f)塑性相关项：(2)所以，复合加载下的应力强度因子初始参考应力：稳态蠕变复合应力场下c*积分数值为：复合加载下的j积分值：(3)场变量可以获取应力值查表得：in＝4.99，p92钢的材料参数εcrit＝0.2；n＝5.23,在计算蠕变应力以及拘束时候，我们取裂尖前的距离r＝d＝0.05mm。查表得：塑性等效应力查表得：瞬态蠕变等效应力转换时间thrr-rr：利用以及matalab计算得：thrr-rr＝2h查表得：塑性阶段损伤累计值：(4)然后计算瞬态蠕变应力场下发生的萌生：查表得：平均应力：应力三轴度：多轴应力因子：是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离，即蠕变萌生发生的临界距离，一般取所研究材料的晶粒尺寸，如图2所示。瞬态蠕变条件下孕育期：利用matalab解积分得：tihrr-rr＝1822h。以上所述仅是本发明的优选实施方式，应当指出的是，对于本
技术领域：
：的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。当前第1页12当前第1页12