刘正奇,高付海
(中国原子能科学研究院 核工程设计研究所,北京 102413)
核压力容器高温设计规范的提出是为了核设施部件在高温下使用时有统一的设计分析评定准则,它的完善和发展与科学技术的进步、人们对失效模式的认识联系紧密,理想情况下核压力容器高温设计规范应包括:高温的定义,考虑的设计载荷与失效标准,设计准则及评价。
目前国内核承压设备设计使用的主流规范是ASME,其在1915年首次颁布,该规范在力学评价弹性分析中首次提出了应力分类法的思想,为后续压力容器力学评价的弹性分析方法奠定了基础,影响深远。2015年正式发布的ASME规范第Ⅲ卷第5册(简称ASME-Ⅲ-5)[1]是研究核容器中关于高温环境运行下的规范,综合了ASME自1976年以来高温设备设计规范的主要成果[2-3],涉及蠕变、热老化等失效机制,提供了核1级部件设计规则。此外随着时间的发展,国际上许多其他国家核工业因为不同的需求和侧重,也诞生了各自的核承压设备高温设计规范,且各有所长,著名的如:英国R5规程[4],注重参考应力法和极限分析的应用;法国的RCC-MRx规范[5],与ASME一脉相承,但区分蠕变显著/不显著、辐照显著/不显著,更加细化;日本的MONJU规范[6-7],注重地震引起的失效和长/短期载荷作用对不同失效机制的影响。
随着当下计算方法、计算机高性能计算技术、商用有限元软件和高温非弹性分析方法技术的快速发展进步,相较而言,ASME规范传统采用的作为主流的基于弹性分析的应力分类评价方法本身存在结果过于保守、设计经济性不够等局限性。本文通过对比ASME-Ⅲ-5与国际上其他主流设计规范或规程,基于损伤模式下载荷控制的应力限制和变形控制限制分析方法及准则,探讨ASME-Ⅲ-5规范在分析方法准则中的优势与局限性,吸收其他规范技术强项,提出若干项改进ASME-Ⅲ-5高温1级部件设计规范的方法,指出ASME-Ⅲ-5高温规范改进的方向,以进一步提高规范发展使用的经济性和安全性。
理想情况下,核反应堆运行环境中的高温指的是温度的一个工作范围,在此范围内,即使在低于屈服强度的理论稳定载荷作用下,与时间有关的热激活变形和损伤过程也可以成为影响承载部件稳定性的重要因素。
由于时间相关性是温度的平滑函数,没有可以清晰区分时间相关和时间无关行为的界限,因此对于特定的材料和使用条件,可以定义一个温度限制即高温条件,低于该限制时,温度影响可以忽略不计。此温度限制可以确定为时间相关性首次达到显著标准的阈值温度。高温条件的定义不唯一,需依据规范考虑的失效模式和设计寿命综合决定,达到或高于高温条件需要考虑时间相关的热激活变形和失效,如热老化、蠕变等,低于时可以忽略时间相关性影响。
ASME规范有高温使用温度和不显著蠕变温度的定义,也给出了与时间、温度相关的材料许用应力限制Smt。对铁素体钢温度超过700 ℉(371 ℃)或奥氏体温度超过800 ℉(427 ℃)时,要考虑蠕变[1],否则可以忽略。ASME-Ⅲ-5关于不显著蠕变的定义是模糊的,交由设计师自行判断,但可以参考附录T中Test-A3试验给定的明确标准,即当蠕变断裂的时间使用分数为0.1或更小,且累积蠕变应变小于限制的0.2%时,蠕变可认为不显著。
ASME-Ⅲ-5有关于不同载荷的分类和处理,包含了设计载荷、正常工作载荷、经常出现的异常载荷、罕见异常载荷、极限故障载荷(分别对应规范的设计,A、B、C、D级使用限制)以及试验载荷,载荷类型可分为机械载荷(如自重、内压、外压、地震、流体动力载荷、接管载荷等)、热载荷(如温度引起的热应力)和热-机循环载荷,不同载荷或组合引起不同的失效机制,如机械载荷一般引起载荷控制的一次应力失效(即载荷控制的应力限制),热载荷一般引起变形控制的二次应力,热-机循环载荷引起变形控制的安定棘轮和疲劳失效(即变形控制限制)。
在ASME-Ⅲ-5设计准则中涉及的失效机制分为单一载荷、循环载荷、蠕变相关的失效,包括:单一载荷下的极限垮塌,单一载荷下的结构不稳定或屈曲;循环载荷下的疲劳失效,棘轮作用导致的渐进坍塌;高温运行下的特定失效包括蠕变断裂、增强蠕变、蠕变屈曲、蠕变疲劳、弹性跟随效应。通过免责声明进行说明的有:单一载荷时的过度变形;稳定载荷下的蠕变产生过度变形;由于腐蚀、传质现象等导致压力边界坍塌或破裂;辐照效应。
ASME-Ⅲ-5以详细表格的形式总结提供了如何使用设计规范来处理各种失效机制,其中的重点是Ⅲ-5册不允许对载荷控制的应力限制(即一次应力限制)进行高温非弹性分析,D级载荷除外。其载荷控制的应力限制分析中采用一种弹性方法来解释高温蠕变引起的一次弯曲应力的重新分布,通过一次膜应力加上考虑或不考虑蠕变影响的主弯曲应力限制,有效形成了一次应力的破坏标准。ASME-Ⅲ-5附录T提供了变形控制的限制,即关于沿厚度方向平均应变、表面平均加弯曲应变和峰值应变的最大累计非弹性应变不能超过1%、2%和5%,焊接件中的应变集中则通过将许用应变减少为母材的1/2来解决[1]。上述应变限制是通用的,不考虑给定材料的蠕变延性。ASME-Ⅲ-5没有提供解决腐蚀、辐照和传质现象的方法。
比较从载荷控制的应力限制和变形控制限制两方面进行,每个方面分为时间相关和时间无关领域,进而可以直观看出ASME设计规范的弹性分析方法、非弹性分析方法的侧重点和优劣势。ASME-Ⅲ-5时间无关领域直接引用ASME第Ⅲ卷第1册NB分卷。
1) 载荷控制的应力限制
(1) 时间无关
在弹性分析中忽略蠕变的情况下,ASME规范关于一次、二次应力以及峰值应力的应力分类的定义在世界压力容器设计领域具有突破性意义,对国际上其他不同设计标准的后续发展和要求都有重大影响。
关于载荷控制的应力限制的计算,ASME允许使用非弹性方法的极限载荷分析和弹塑性分析。
(2) 时间相关
弹性分析时,ASME承认随蠕变发生的应力重分布现象,同时其使用应力折减因子Kt来反映高温蠕变引起的沿壁厚的外层纤维弯曲应力减小。
此处ASME规范不允许在载荷控制的应力限制分析中使用非弹性分析和极限载荷分析[8]。
目前ASME规范使用存在的一项主要困难是面对复杂结构或工况条件下如何进行合理的应力分类,在对应力分类把握不准时,设计师一般采用保守处理方式,这往往可能导致不恰当的设计结果。
2) 变形控制限制
变形控制限制的介绍和对比主要从弹性跟随现象的处理和弹性分析、非弹性分析方法3个方面进行。
弹性跟随指的是不同刚度零部件连接在一起且柔性部分处于高应力状态时,柔性部分可能产生很高的应变集中。由于弹性分析方法只是进行了弹性计算,并不能真实反映出结构实际存在的弹性跟随这一重要效应,这是弹性分析方法的一个固有不足,规范必须加以考虑解决。在变形控制限制中,进行安定棘轮分析时ASME规范简化非弹性方法中直接将弹性跟随引起的膜应力视为一次应力(即载荷控制的应力)进行保守评估。完全非弹性分析可以对弹性跟随进行更进一步的评估,ASME规范是允许的,但依赖高温本构模型[9]。
ASME规范在变形控制限制中主要采用的弹性分析和简化非弹性分析方法是指使用弹性应力分析和应力分类,其结果应用于弹性、安定和棘轮分析的Bree图(适用于理想弹塑性结构和小变形假设,以影响热应力棘轮分析的3个无量纲参数为坐标轴的图,用以分析安定棘轮状态),即附录T中的A &B试验。
ASME规范附录T部分允许对变形控制限制的计算进行非弹性分析,基于材料的弹黏塑性本构模型和有限元分析软件,在A、B、C使用载荷、使用寿期内得到的最大累计非弹性应变应满足:沿厚度方向平均应变、表面平均加弯曲应变和峰值应变的应变不超过1%、2%和5%。ASME规范Ⅲ-5(2021版)提供了若干种高温材料的弹塑性、蠕变本构模型,可以据此开展高温结构的基于分离型本构模型的弹性跟随、安定棘轮、蠕变-疲劳分析,但提供的模型参数正确性仍存在一些问题[10]。高温非弹性本构模型的应用方法在ASME规范中只是简要提及,并未深入涉及设计师操作层面的细则规定[9],由于现今计算机高性能计算技术的巨大进步,这会是ASME规范未来进一步细化改进发展的一个方向。
1) 载荷控制的应力限制
(1) 时间无关
弹性分析方面,R5直接借鉴了ASME对一次应力、二次应力以及峰值应力的定义。
非弹性分析方面,与ASME不同,R5允许并鼓励使用非弹性分析,在第2/3卷附录A12节中提供了非弹性分析和各种形式本构方程的使用指南。同时R5允许并在很大程度上鼓励使用极限载荷方法,此类方法的使用以及与设计标准的结合非常出色,充分体现了R5前沿性、学术性更强的技术特点。
(2) 时间相关
弹性分析方面,R5并不提供一种用弹性方法解释蠕变变形引起的应力重分布[4],与ASME相反,R5更加依赖极限载荷分析和参考应力方法。极限载荷分析实际上是极限分析的一种类型,极限分析可被视为寻找极端条件的分析,如最大允许载荷或最小不安全负荷;这还包括循环状态,如安定和棘轮极限。极限载荷分析通常通过使用有限元分析来实现,其假定结构由理想塑性材料制成,采用极限分析法来计算结构最大承载能力。R5规范将极限载荷与参考应力联系起来,以分析一次载荷限制,其中参考应力指给定外载荷条件下复杂构件的稳态蠕变行为可以跟单轴蠕变试验结果建立关联的应力[11]。1968年,Sim[12]在其博士论文中建立了参考应力和极限载荷、所受外载荷的关联式:
σref=P0/Puσy
(1)
式中:σref为参考应力;Pu为结构的极限载荷;σy为结构材料的屈服极限;P0为结构的外加载荷。
此方法于1991年被R5规范采纳后用于蠕变强度分析。同时,参考应力法基于结构整体失稳而不是局部破坏的极限载荷,其本质上要求材料必须具有足够的高温蠕变延性才能使用[11],不过R5规范同时提供了蠕变脆性材料的参考应力法。
非弹性分析方面,R5同样允许并鼓励使用非弹性分析。
与ASME规范不同,R5规范尽管也采用了应力分类的技术路线,但极限载荷分析及其在参考应力方法中的应用是R5的时间相关载荷控制的应力限制分析方法的核心。通过限制参考应力小于断裂参考应力判据,保障满足蠕变强度和设计寿命要求[11],根据蠕变延性/脆性材料有不同的断裂参考应力计算形式。此方法避免了ASME规范经常面对的复杂结构如何进行应力分类的问题,提高了方法应用的客观性和设计结果准确度。
2) 变形控制的限制
R5认识到弹性跟随对结构性能分析的重要性。关于循环载荷,有3种方案来解决变形限制问题,分析程度和保守程度各不相同,分别是弹性分析、简化非弹性分析和完全非弹性分析。这里介绍基于简化非弹的弹性跟随因子法。R5为了量化结构中不连续部位的弹性跟随效应定义了弹性跟随因子,提供了其求取流程[11]:
(2)
ASME没有对弹性跟随进行量化处理。
R5允许对变形限制进行弹性分析,并将结果应用于蠕变不显著的快速循环,同样也允许考虑增量时间步长使用非弹性分析来分析结构的弹性、塑性和蠕变行为。与载荷控制的限制类似,直接使用非弹性分析方法的计算效率不高。此外R5若要在棘轮状态下运行,则必须进行完全非弹性分析。
相比ASME规范使用时间分数法计算蠕变损伤,而R5规范更推荐使用延性耗竭法,延性耗竭法基于延性耗竭理论,认为当局部区域累积蠕变应变达到断裂延性值时,损伤达到临界值,即失效[13]。此方法计算的是蠕变损伤。蠕变延性数据为断裂时蠕变应变,可通过单轴蠕变试样截面收缩率换算或延伸率获取。
R5还提供了分析安定棘轮边界的直接法,即线性匹配法(Linear Matching Method, LMM),其通过一系列迭代修正弹性模量的线弹性分析来拟合结构的塑性行为[14],是一种求解材料非线性行为响应的快速直接方法,其不同于R5规程中的其他方法以校核为主,该方法计算载荷极限,是真正的设计方法。
RCC-MRx对失效机制分为P型损伤和S型损伤两大类,其描述了在两种主要失效机制下1级部件的一般设计规则。P型损伤是由稳定单调增加的载荷或恒定载荷造成的,这包括瞬时过度变形、塑性失稳、时间相关过度变形、塑性、断裂和弹性或弹塑性失稳。S型损伤是指由于循环往复加载而产生的所有失效机制,如渐进变形(棘轮)和疲劳开裂。在区分P型和S型损伤后,还区分了蠕变显著/不显著、辐照显著/不显著[15],再进行后续强度分析。RCC-MRx关于辐照能否忽略有明确标准,通过结构使用期间接收的最大剂量小于限值和每个时间段中累计剂量小于此时间段内最高温下的允许累计剂量来区分,在辐照不可忽略时许用应力强度加入了辐照项,两者都不可忽略时蠕变使用分数限制降低到考虑蠕变、不考虑辐照的10%[5]。而ASME规范则完全未考虑辐照的影响。下面的对比为不考虑辐照的情况。
1) 载荷控制的应力限制
(1) 时间无关
弹性分析方面,RCC-MRx以与ASME相同的方式对一次应力、二次应力和峰值应力进行分类。
非弹性分析方面,如极限载荷法,与ASME相同,RCC-MRx允许使用非弹性分析,它提供了描述可以使用的非弹性分析方法和本构材料模型的附录。
(2) 时间相关
RCC-MRx可以使用弹性分析,与ASME相同,RCC-MRx使用应力折减因子解释蠕变变形引起的应力重新分布。同时RCC-MRx使用修正因子修正了弹性分析中弹性跟随对膜应力的影响[5]。
RCC-MRx允许对时间相关的应力限制进行非弹性和极限分析(不可与其他方法混用),ASME规范则不允许。
2) 变形控制的限制
RCC-MRx同样认识到需要在弹性分析方法中解决弹性跟随问题,但处理不充分,只是提供一个免责声明。
弹性分析方面,RCC-MRx有完全不同于ASME的方法,两者差异较大,该方法仅限于布里管型问题,未对一般结构和复杂载荷的适用性进行说明,其需要定义效率指数和二次比,很难评估是否比ASME规范提供的方法更有效。
非弹性分析方面,RCC-MRx允许使用弹塑性和黏弹塑性类型的非弹性分析。RCC-MRx提供了该规范所涉及的高温使用结构材料的塑性、蠕变本构模型方程形式和参数。
MONJU规范很大程度上基于ASME规范,主要区别在于MONJU区分了长/短期载荷,尤其是对于日本经常发生的地震载荷,提出了一种与地震载荷相关的新型失效机制,称为棘轮疲劳[6]。其是一种由渐进循环变形加速疲劳的失效机制,并主要发生在管道中,该机制是短期的,与蠕变不直接相关,因此它在高温或低温下都可能发生。
1) 载荷控制的应力限制
(1) 时间无关
弹性分析方面,应力分类的定义与ASME相同,区别是区分了短期和长期应力,限制两者之和小于限制,这反映了日本对地震载荷的特殊关注。
非弹性分析方面与ASME相同。
(2) 时间相关
弹性分析方面与ASME一致。非弹性分析方面,有新提出类似R5的基于参考应力的无应力分类设计,不过无进一步的详细信息,说明其正处于发展之中。
2) 变形控制的限制
MONJU更多地考虑了弹性跟随现象,与ASME规范类似,将热应力视情况分为二次应力或纳入一次应力分析,在低应力、长期运行的情况下简单几何结构的弹性跟随应变计算参照提供的附录。
弹性分析方面,MONJU有一种处理循环载荷和相关变形极限的简化方法,类似于ASME附录T的B试验。
非弹性分析方面没有详细资料。
总的来说,R5允许使用各种分析类型,如使用简化非弹性分析以计算载荷控制的应力限制和变形控制极限,但当使用弹性分析时,ASME采用整体厚度的平均温度,而R5则采用所选相关截面的最高温度,R5相对更保守。参考应力法是其简化非弹性分析的核心。使用参考应力法进行蠕变强度分析可以有效避免ASME如何合理进行应力分类的问题,提高了设计准确度,同时弹性跟随因子应用也很好地解决了弹性跟随问题。
RCC-MRx与ASME不同,区分了载荷类型及蠕变是否显著/不显著,考虑了辐照的影响。其允许对P型载荷的蠕变强度进行极限分析,提出了修正因子修正弹性跟随对膜应力的影响,也允许对载荷控制的应力限制进行非弹性分析。
MONJU规范大体上基于ASME,但它重点区分了长/短期载荷和应力限制,提供了如何判断弹性跟随程度的详细指南,并基于地震强震载荷引起的失效机理新提出了一种失效模式——棘轮疲劳。
通过上述对比,可知ASME-Ⅲ-5对高温环境下载荷控制的应力限制仅允许进行弹性分析,不允许使用非弹性分析和与极限载荷有关的参考应力法,在变形控制限制的弹性分析方面对弹性跟随效应进行简单保守处理,非弹性分析中没有高温非弹性蠕变本构模型使用的细则规定,没有考虑辐照相关失效等。所以本文尝试总结面向损伤模式下各个设计方法的优势和劣势并进行改进讨论。
基于R5多年实践经验,可以看出在特定情况下蠕变相关的载荷控制的应力限制分析中使用参考应力法替代传统的应力分类法是可行的,传统应力分类方法在面对部分复杂结构的蠕变强度分析过于保守。MONJU规范中针对地震的短期、长期载荷区分,RCC-MRx关于辐照限制的影响也值得参考。
变形控制限制中R5提出了弹性跟随因子及其求取流程,在简化非弹性分析中处理了弹性跟随效应,而ASME规范处理得不充分。
已被R5规范采纳的线性匹配法是一种先进的数值算法,适用于复杂工况下的棘轮极限载荷分析,具有较高的计算效率。
1) ASME-Ⅲ-5规范主体部分可以不变,对载荷控制的应力限制的分析方法中除了强制性的弹性分析外可以引入非强制的可替代的非弹性分析方法,如参考应力法,并提供配套的更详尽的材料数据,如蠕变延性。参考应力法基于极限分析,避免了应力分类,给蠕变强度分析带来了极大地便利,也更加合理。在一次载荷分析领域中,参考应力法已有了足够的数据支持,在其他国际ETD规范和指南中被采纳或拟议采纳,可作为一种可选方法纳入ASME准范中。
2) 在载荷控制的应力限制中,可以加入考虑蠕变显著/不显著、辐照显著/不显著的影响以及地震条件下对于长/短期载荷明确采用不同应力限制。
3) 除ASME现有的弹性分析、简化非弹性分析和非弹性分析3种方法外,可以研究考虑将目前国际上先进的直接法设计方法——线性匹配法纳入,单独作为一种与上述3种方法并列的设计师可选的设计方法。
4) 高温非弹性本构模型是高温结构非弹性应力应变响应分析的基础,现有ASME规范提供了若干材料的分离型本构模型和参数,但对于循环塑性本构、适用于9Cr1MoV材料的统一型非弹性本构模型以及非弹性分析方法的使用细则等内容ASME都未明确,有待细化规定。
5) ASME规范中变形控制限制的弹性分析对弹性跟随处理不充分,可以在可选替代方案中考虑引入弹性跟随因子及其求取流程,量化考虑弹性跟随效应。
6) ASME规范中变形控制限制的蠕变疲劳损伤交互作用判据对于蠕变损伤计算采用考虑时间分数法,结果往往过于保守,在一定条件下可以考虑研究引入保守性适度的延性耗竭法作为可选替代方案。
ASME规范历史悠久,是实践了很久且行之有效的高温设计准则。为了更好地认识和使用规范,本文将ASME-Ⅲ-5与其他国际主流规范或规程对比,探讨改进方向,得出如下主要结论。
1) ASME-Ⅲ-5基于应力分类的弹性分析方法思想先进,表现优异,但存在复杂结构应力难以合理分类问题,且在高温使用时不允许对载荷控制的应力限制进行非弹性分析。鉴于以上原因以及现今计算技术的巨大进步,ASME可以尝试引入非弹性分析法如R5规程的参考应力法作为一种替代方法,在保证安全性前提下提高设计方法的经济性。也可以考虑引入线性匹配法作为一种独立的设计方法供设计师选择作为可选方案。
2) ASME-Ⅲ-5对于载荷控制的应力限制可以考虑针对服役条件和地震条件进行进一步的明确细化,加入考虑蠕变显著/不显著、辐照显著/不显著的影响以及地震条件下对于长/短期载荷明确采用不同应力限制。
3) 高温非弹性本构模型是高温结构非弹性分析的基础。对于循环塑性本构、适用于9Cr1MoV材料的统一型非弹性本构模型以及非弹性分析方法的使用细则等内容ASME可以进一步细化规定。
4) ASME规范中变形控制限制的蠕变疲劳损伤交互作用判据对于蠕变损伤计算方法过于保守,在一定条件下可以研究考虑使用保守性适度的延性耗竭法作为可选的替代方案。
本文对比研究指出了若干项ASME-Ⅲ-5高温1级部件规范设计方法准则技术改进的方向,有望通过进一步的技术基础研究正式纳入ASME,进而提高规范发展使用的先进性、经济性和安全性,同时以上结论对开发构建我国自主化的核承压高温使用1级部件设计规范亦具有重要指导价值。