[发明专利]针对蠕变损伤及体积型缺陷的承压结构的安全评定方法

摘要

本发明针对蠕变损伤及体积型缺陷的承压结构的安全评定方法，包含以下步骤⑴对长期工作的承压结构建立弹塑性本构模型，用以预测所述承压结构的弹塑性响应，为估算其承载极限做准备；⑵运用有限元分析方法研究带各种凹坑及局部减薄缺陷的承压结构在内压或弯矩载荷作用下塑性区域的扩展过程和失效模式，研究各种凹坑及局部减薄所致的不同蠕变损伤对承压结构承载极限的影响，定量地算出承压结构的承载极限值；⑶确定安全评定步骤，对含有蠕变损伤和/或体积型缺陷的承压结构进行安全评估。本发明首次将蠕变损伤与体积型缺陷尺寸的影响紧密结合，制定出适合高温高压下长期工作的承压结构的安全评定标准，具有明确的物理意义，适合在工程中使用。

主权项

针对蠕变损伤及体积型缺陷的承压结构的安全评定方法，其特征在于，包含以下步骤：(1)对长期工作的承压结构建立弹塑性本构模型，用以预测所述承压结构的弹塑性响应；在本步骤中，基于Ramberg‑Osgood模型构建一个耦合材料各向同性的蠕变损伤的弹塑性本构模型，将所述承压结构所受到的蠕变损伤耦合到弹塑性本构模型中，预测长期服役后的承压结构的弹塑性响应，进而为估算其承载极限做准备；所述弹塑性本构模型的构建方法为：①不考虑氧化、腐蚀等环境因素的影响，从微观组织的角度分析，所述蠕变损伤受三种因素的综合影响：粒子粗化、固溶元素贫化以及微空洞形核长大；设：由粒子粗化导致的蠕变损伤为DP；由固溶元素贫化导致的蠕变损伤为DS；由微空洞形核长大导致的蠕变损伤为DN；DP的损伤速率为：D·P=KP(1-DP)43,DP(t0)=0,DP(tr)=1---(1-1)]]>式中，t0,tr为蠕变开始时间和断裂时间；KP为粒子粗化的速率常数；DS的损伤速率为：D·S=KSDS1/3(1-DS),Ds(t0)=0,Ds(tr)=1---(1-2)]]>式中，KS为粒子析出的速率常数；DN的损伤速率为：式中，A，υ，是与最小蠕变应变速率以及断裂相关的材料常数；②将上述蠕变损伤的三种因素进行综合考虑，设总的蠕变损伤为D；D的损伤速率为： D＝1‑(1‑Dp)(1‑Ds)(1‑DN)        (1‑4)式中，Dp,Ds,DN通过求解公式(1‑1)、(1‑2)、(1‑3)的微分方程并根据相应的蠕变时间来确定；③基于Ramberg‑Osgood模型构建弹塑性本构模型步骤(1)所述的Ramberg‑Osgood模型与弹塑性本构基本的关系为：ϵ=σE+α(σσp)n---(1-5)]]>式中，E为弹性模量，ε为应变，σ为应力，σp为参考应力，n为非线性项的硬化参数，α为屈服偏移量；将求解公式(1‑4)中的D代入Ramberg‑Osgood模型中，得到长期工作的承压结构蠕变损伤的弹塑性本构模型：ϵ=σE(1-D)+α(σ(1-D)σp)n---(1-6);]]>(2)在步骤(1)的基础上，得到含有蠕变损伤材料的力学特性，运用有限元分析方法，研究带各种形状、尺寸的凹坑及局部减薄缺陷的承压结构在内压或弯矩载荷作用下塑性区域的扩展过程和失效模式，研究各种凹坑及局部减薄缺陷所致的不同蠕变损伤对所述承压结构承载极限的影响，定量地算出所述承压结构的承载极限值；所述的有限元分析方法是指基于增量有限元理论的分析方法，所述方法注重载荷加载过程，能关注到加载过程中结构应力‑应变场的变化；其确定极限载荷的原理是：在有限元分析中将载荷不再增大或者有微量增大而位移或应变增量很大时的载荷确定为极限载荷，其步骤包括：①载荷加载采用自动步长，迭代至有限元计算发散为止；②通过已加载到发散的步长，得到对应的总的加载步长，然后确定发散时加载的载荷值；③分析重点部位的加载步和应变的关系，通过载荷‑位移曲线，确定结构的极限载荷；(3)根据步骤(2)有限元的分析和计算结果及拟合公式，确定安全评定步骤，对含有蠕变损伤和/或体积型缺陷的承压结构进行安全评估，从而评估该承压结构的安全状态；所述的安全评估包括：①对含凹坑的承压结构的极限载荷的评估Ⅰ、参数的无量纲化设：缺陷尺寸的无量纲参数为g0，其表达式为g0=CT·AB·BRm·(RmT)12=CT·ARmT---(1-7)]]>式中，A,B,C分别代表缺陷轴向、环向和深度方向的尺寸；Rm，T表示圆筒的中径和壁厚；因含凹坑的承压结构的极限载荷对上述无量纲参数g0的平方根更加敏感，故G0=g0---(1-8)]]>Ⅱ、蠕变损伤的无量纲化根据步骤(1)，蠕变损伤与蠕变时间的有一一对应的关系，故用相对的蠕变时间tC表示蠕变损伤，tr为蠕变断裂时间；tC＝t/tr                        (1‑9)Ⅲ、极限载荷的无量纲化极限载荷PLRC用PLRC＝PL/PL0来进行归一化，式中，PL和PL0分别为有缺陷和无缺陷结构的极限载荷；Ⅳ、极限载荷与凹坑缺陷尺寸G0及蠕变损伤tC的关系根据所述步骤(2)的计算结果，对于不同tC下的数据进行拟合，出于保守估计，考虑蠕变损伤的含凹坑圆筒的极限内压载荷拟合公式为PLRC=e-tC-0.3G0tC≤0.5PLRC=e-tC-0.2G0tC>0.5---(1-10)]]>②对含局部减薄缺陷的承压结构极限载荷的评估Ⅰ、参数的无量纲化Ⅰ‑1、缺陷尺寸的无量纲化设无量纲参数AeAe=c·abc3---(1-11)]]>其中，a,b,c分别代表缺陷轴向、环向和深度方向的相对尺寸。Ⅰ‑2、蠕变损伤的无量纲化步骤同①Ⅱ；Ⅰ‑3、极限载荷的无量纲化步骤同①Ⅲ；③极限载荷与减薄缺陷尺寸Ae及蠕变损伤tC的关系Ⅰ、纯内压下，根据所述步骤(2)的计算结果，对不同tC下的数据进行拟合，出于保守估计，考虑蠕变损伤的含凹坑圆筒的极限内压载荷拟合公式为：PLRC=e-tC(0.95-0.85Ae)a/b≤7PLRC=e-tC(0.95-1.04Ae)7<a/b≤25PLRC=e-tC(0.95-1.47Ae)a/b>7---(1-12)]]>Ⅱ、纯弯矩载荷作用下蠕变损伤与极限载荷的关系，拟合得到公式，为区别于内压，弯矩作用下的极限载荷用mLs来表示：mLs=e-tC(cos(cπb2)-csin(πb)2)c<1-bbmLs=e-tC(1-c)sin[π(1-bc)2(1-c)]+csin(πb)2c>1-bb---(1-13).]]>