中国腐蚀与防护学报, 2025, 45(2): 307-318 DOI: 10.11902/1005.4537.2024.157

临氢关键材料服役行为研究专刊

重水堆压力管延迟氢化物开裂行为研究进展

潘春婷1,2, 明洪亮,1,2, 石秀强3, 鲍一晨3, 王俭秋1,2,4, 韩恩厚2,4

1.中国科学技术大学材料科学与工程学院 沈阳 110016

2.中国科学院金属研究所 沈阳 110016

3.上海核工程研究设计院股份有限公司 上海 200233

4.广东腐蚀科学与技术创新研究院 广州 510530

Research Progress on Delayed Hydrides Cracking Behavior of Heavy Water-Reactor Pressure Tube

PAN Chunting1,2, MING Hongliang,1,2, SHI Xiuqiang3, BAO Yichen3, WANG Jianqiu1,2,4, HAN En-Hou2,4

1.School of Materials of Science and Engineering, University of Science and Technology of China, Shenyang 110016, China

2.Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China

3.Shanghai Nuclear engineering Research and Design Institute Co., Ltd., Shanghai 200233, China

4.Institute of Corrosion Science and Technology, Guangzhou 510530, China

通讯作者: 明洪亮,E-mail:hlming12s@imr.ac.cn,研究方向为材料力学化学交互作用

收稿日期: 2024-05-18   修回日期: 2024-06-27  

基金资助: 中国科学院青年创新促进会.  2022187

Corresponding authors: MING Hongliang, E-mail:hlming12s@imr.ac.cn

Received: 2024-05-18   Revised: 2024-06-27  

Fund supported: Youth Innovation Promotion Association CAS.  2022187

作者简介 About authors

潘春婷,女,1997年生,博士生

摘要

压力管是重水堆中重要的结构部件,其服役性能是确保重水堆安全稳定运行的关键。延迟氢化物开裂(DHC)是压力管服役过程中最重要的潜在风险之一,因此研究压力管的DHC行为具有重要的意义。本文就DHC行为测试的实验方法、机理和模型以及DHC行为的影响因素等方面的研究进展进行综述,并指出目前研究存在的不足和未来的发展趋势。

关键词: 重水堆 ; 压力管 ; Zr-2.5Nb合金 ; 延迟氢化物开裂

Abstract

Zr-2.5Nb alloy pressure tubes are important structural components in heavy-water reactor. During operation of a heavy-water reactor, a large amount of hydrogen isotopes is produced by the corrosion reaction between the pressure tubes and the heavy-water coolant. Some of the hydrogen isotopes can be absorbed into the pressure tube. When the concentration of hydrogen atoms exceeds the solid solubility of hydrogen in Zr-2.5Nb alloy, hydrides are precipitated. The precipitation of hydrides will lead to deterioration of mechanical properties of Zr-2.5Nb alloy, and then leads to the expansion of microcracks inside the pressure tubes. The phenomenon is called Delayed hydrides cracking (DHC) which is one of the most important potential risks during the service of pressure tubes. Therefore, it is of great significance to study the DHC behavior of pressure tube. In this paper, the research progress on the testing methods for DHC behavior, the relevant mechanisms and models as well as the influencing factors of DHC behavior are reviewed, and the shortcomings of current researches and the future development trends are pointed out.

Keywords: heavy-water reactor ; pressure tubes ; Zr-2.5Nb alloy ; delayed hydrides cracking

PDF (2070KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

潘春婷, 明洪亮, 石秀强, 鲍一晨, 王俭秋, 韩恩厚. 重水堆压力管延迟氢化物开裂行为研究进展. 中国腐蚀与防护学报[J], 2025, 45(2): 307-318 DOI:10.11902/1005.4537.2024.157

PAN Chunting, MING Hongliang, SHI Xiuqiang, BAO Yichen, WANG Jianqiu, HAN En-Hou. Research Progress on Delayed Hydrides Cracking Behavior of Heavy Water-Reactor Pressure Tube. Journal of Chinese Society for Corrosion and Protection[J], 2025, 45(2): 307-318 DOI:10.11902/1005.4537.2024.157

核能发电以其燃料能量密度高、基本零排放和土地占用少等优势已经成为一些发达国家的首选和我国能源发展的重要趋势之一。根据国际原子能机构(IAEA)统计,截至2023年10月,中国的在运核电机组规模居世界第三,居于世界领先地位;《中国核能发展报告2023》指出,未来我国核能发电占比将会逐步提高,保证核电产出的高效和安全是当前国家发展战略之一。核电站的安全运行始终是扩大核电规模的重要前提,反应性控制、余热排出和包容放射性物质是核电厂三项基本安全功能,其中“包容放射性物质”这一项对核电关键性结构材料在苛刻工况中服役的稳定性和可靠性提出了极高的要求。锆合金具有力学性能优良、耐腐蚀强、中子吸收截面小和热导率高等优点,广泛应用于核反应堆的燃料包壳和堆芯结构部件的制造[1]。其中,重水堆中的重要结构部件压力管就由Zr-2.5Nb合金制成。压力管作为重水堆中燃料棒束和冷却剂的容器,是重水堆一回路燃料通道的重要组成部分[2]

19世纪70年代,研究人员首次在加拿大重水铀反应堆(CANDU)中发现了压力管内部氢化物引起压力管内容物泄漏现象。研究表明,压力管的破裂是由一种现在被称为延迟氢化物开裂(DHC)的裂纹扩展造成的[3,4]。在反应堆运行过程中,压力管与重水发生腐蚀会反应产生大量的氢原子,由于锆金属与氢结合力较强,氢原子极易进入压力管内部。在应力场的作用下,氢原子易向高应力区(如微裂纹尖端)扩散,当高应力区的氢含量超出其极限固溶度时,难以继续溶解的氢原子就会以氢化物脆性相的形式析出并稳定存在于锆合金基体中。氢化物脆性相的存在会降低压力管材料的断裂韧性,当氢化物在应力作用发生断裂时就会导致原有微裂纹的扩展或产生新的微裂纹,上述氢扩散、氢化物形核和生长、氢化物断裂过程的反复发生,就会导致宏观上DHC裂纹的扩展,最终致使重水冷却剂的泄露,对重水堆的安全运行造成影响。目前,DHC被视为压力管服役过程中最重要的潜在风险之一。来自加拿大、韩国、阿根廷和印度等国家的研究者均对重水堆压力管用Zr-2.5Nb合金的DHC行为开展了大量研究,其中既包括DHC机理、模型构建及可靠性验证,也较为系统的研究了温度、氢含量、热循环方式以及锆合金本身微观结构和力学性能等因素对DHC行为的影响。DHC扩展速率(DHCR)作为评估压力管是否符合“破前漏”(LBB)准则的重要指标之一[5,6],是评价压力管DHC行为的重点。此外,DHC作为一种受环境影响的亚临界裂纹扩展[7],其萌生条件之一为应力强度因子达到DHC的门槛值KIH[8],为了预防和抑制DHC的产生,KIH测定是评价压力管DHC失效的另一个重点。

当前我国仅有秦山三期建造的两座重水堆核电站。秦山三期重水堆核电站自2003年7月建成投产以来,已服役超过20年(CANDU6型重水堆核电机组设计寿命为40年),根据加拿大标准协会相关标准CSA N285.4的要求,应定期对压力管取样并进行服役性能测试(其中包括DHCR测试和KIH测试),从而为压力管的老化评估和寿命预测提供必要的数据支撑。此外,根据秦山三期重水堆设计文件手册及最终安全分析报告,机组运行中期需更换压力管。为提高核电站中国产零部件的比例,降低压力管更换的经济成本,减少进口依赖,近几年我国陆续开展了国产压力管的研制工作。在国产压力管正式投入使用前,需对其抗DHC性能进行测试和验证。目前国内对压力管DHC行为的相关实验研究较少,因此建立DHC行为相关测试方法,厘清各因素对DHC行为的影响,对压力管的寿命评估与延寿以及国产压力管制造和性能测试有重要意义。本文就常用DHC测试方法、DHC相关机理以及各种因素对DHC行为的影响等方面进行综述,指出当前的研究不足,并对未来的研究趋势进行了展望。

1 DHC行为研究方法

为模拟压力管服在富氢环境下大量吸氢的现象,需在DHC实验前对压力管材料进行预充氢实验,通常用的充氢方法有气体充氢、高压釜渗氢、电解充氢等[9~11]。其中,电解充氢所需设备简单易得、操作便捷且可重复性高,应用最为广泛。由于压力管材料的各向异性,研究者结合压力管的形状和实际服役工况,通常采用两种不同的试样分别用于研究压力管轴向和径向DHC行为。其中,测试径向DHC行为采用的是悬臂梁(CB)试样;测试轴向DHC行为采用的是弯曲紧凑拉伸(CCT)试样,其取样方式如图1所示。为促进DHC裂纹的萌生,缩短加载与裂纹萌生之间的孕育期,通常需要在CCT试样上预制疲劳裂纹;受限于CB试样的形状,在CB试样上引入疲劳裂纹较为困难,因此在加工试样时尽量降低缺口底部几何形状的直径。图2为测试过程中两种试样的受载示意图,除图2中所示的基本受载条件外,DHC实验还需要高温炉等能精确控制温度变化的设备加以辅助,常采用声发射法(AE)或直流电位降法(DCPD)监测裂纹的萌生与扩展[12~14]。由于DHCR和KIH是评估压力管服役性能的重要指标,目前中国科学院金属研究所已经建成了压力管轴向DHCR和KIH的测试方法,并联合上海核工程研究设计院股份有限公司、中核核电运行管理有限公司、广东腐蚀科学与技术创新研究院和核工业标准化研究所编制了针对压力管轴向KIH和DHCR测试的中国核学会团体标准3项,预计将于近期发布实施。

图1

图1   试样取样及尺寸示意图

Fig.1   Schematic diagrams of sampling (a) and dimensions of CCT specimen (b) and CB specimen (c). B is thickness of pressure tube in Fig.1b and 1c, the unit is mm


图2

图2   试样加载示意图

Fig.2   Schematic diagrams of loading for CB specimen (a) and CCT specimen (b). F is the external load


1.1 KIH 测试

常见的KIH测试方法包括降K法和升K法。其中降K法为先对试样施加较大的载荷,若此时试样发生裂纹扩展,则继续降低载荷值,使K减小。当K值在10 MPa·m-2及以上时,K值每次降幅为2~3 MPa·m-2,当K值降至10 MPa·m-2以下时,K值每次降幅为0.2~1 MPa·m-2,直至试样裂纹保持24 h以上不再扩展(即DCPD电位信号超过24 h未发生明显增长),则认为此时对应的K值即为该温度下的KIH[15,16]。升K法则是先对试样施加较小的载荷,若裂纹未发生扩展,则继续增加载荷,增至裂纹发生扩展,此时的K值即为该温度下的KIH

由于升K法测KIH时,先施加的较小载荷虽不会引起裂纹扩展,但可能会导致裂纹尖端发生钝化[16],从而降低了裂纹尖端的应力集中程度,因此一定程度上会造成所测的KIH值偏大;而降K法实验过程中,施加初始载荷初期就使裂纹尖端产生DHC扩展,可避免裂纹尖端的钝化,保证了KIH测试值的准确性。因此,通常采用降K法测试KIH。当采用降K法测试KIH时,判定所测的KIH是否有效,除需满足DHC裂纹保持24 h以上不发生扩展外,还需满足两个条件:(1) 试样断口表面,靠近试样两侧自由表面的韧性断裂区域的面积小于试样表面总的断裂区域面积的10%;(2) DHC扩展过程中,在试样的厚度方向,裂纹扩展的前沿近似均匀、平直,厚度方向大多数位置的裂纹扩展长度与利用九点法测得的平均裂纹长度之间的偏差小于10%。

1.2 DHCR测试

如前文所述,DHCR是评价材料抗DHC性能和是否符合“破前漏”要求的重要指标之一。如图3所示,DHCR随KI的变化可分为3个阶段,当位于Stage II时,DHCR几乎不随KI的变化而改变[17,18]。因此,在DHCR测试实验中,对试样施加适当恒载荷,使KI处于Stage II阶段。KIH一般小于10 MPa·m-2KIC大于30 MPa·m-2,因此初始KI通常设置在15~20 MPa·m-2的范围内。在此条件下,虽然KI随着DHC扩展而增加,但DHCR基本不受KI改变的影响。实验结束后,通过九点法或面积法测量断口的DHC裂纹长度,裂纹扩展时间则通过DCPD或AE信号的变化获得。由于在Stage II阶段,DHCR对KI无明显依赖性,因此通常将实验过程中的裂纹扩展视为匀速扩展,将裂纹长度与裂纹扩展时间的比值作为DHCR的测试值。

图3

图3   DHCR随KI变化规律[19]

Fig.3   Influence of KI on DHCR[19]. KIH is the threshold of stress intensity factor for DHC, KIC is the fracture toughness


2 锆合金中氢化物析出行为和氢固溶度

2.1 氢化物在锆合金中的析出行为

氢在锆合金中的扩散和聚集是发生DHC现象的首要原因,因此对氢化物在锆合金中析出行为的研究是明确DHC机理的基础。如上文所提到的,压力管在高温高压的服役过程中与重水发生腐蚀反应而产生大量的氢原子:

Zr+2H2O=ZrO2+4H

当进入材料内部的氢原子数量超出氢固溶度时,氢原子会以氢化物的形式析出,造成锆合金部件局部或整体力学性能的下降。Zr为密排六方结构,目前在Zr及其合金中观察到的氢化锆有δ-ZrH1.66γ-ZrH、ε-ZrH2ζ-Zr2H[20,21],最为常见的是前两种。其中δ-ZrH1.66结构最为稳定[22],为面心立方结构,在锆合金中沿特定的惯习面析出[23~25];近期,西安交通大学在Zr-Zn-Nb-Fe合金中观察到应力诱发δ-ZrH1.5在Zr晶体的柱面析出,取向关系为<101¯0> α-Zr//<110> δ[26]γ-ZrH被认为是亚稳相,通常在冷却速度较快的情况下析出[24]。研究认为ε-ZrH2δ-ZrH1.66经过马氏体相变转变形成的[27],在氢含量较高的条件下ε-ZrH2更加稳定,但对于该结论的实验支撑较少。ζ-Zr2H是最晚发现的氢化锆,有学者认为ζ-Zr2H是α-Zr向δ-ZrH1.5γ-ZrH转变时的中间产物[28]

在工程使用前,压力管经历了高温挤压、冷加工和去应力退火等过程,压力管用Zr-2.5Nb合金由α-Zr相和β-Zr相两相组成。α-Zr相被沿压力管轴向方向拉长,β-Zr相呈薄片状分布在两个相邻的α-Zr相晶粒之间[14]δ-ZrH1.66多在α-Zr和β-Zr相界处形核,并由相界处向α-Zr相内部生长形成穿晶氢化物片[29]。Zr-2.5Nb合金中观察到的δ-ZrH1.66析出惯习面多近似平行于(0001) α-Zr[30]。Perovic等[29]研究表明无外加应力条件下,氢化物片沿压力管切向堆叠,在宏观上呈现出氢化物长度方向与压力管切向平行的现象;有外加应力条件下,氢化物片沿压力管径向堆叠,在宏观上呈现出氢化物长度方向与压力管径向平行的现象。切向氢化物通常对压力管影响较小,而外加应力作用下产生的径向氢化物会显著降低锆合金管材的力学性能,其破裂直接导致DHC的萌生和扩展[31,32]

2.2 锆合金中的氢固溶度

锆合金中氢的极限固溶度(TSS)直接影响氢化物能否析出和溶解氢的含量,后者还是氢原子扩散速率的重要影响参数。获得准确的TSS有助于预测压力管的DHC行为,因此许多研究者用差示扫描量热法和红外热导法等测量了不同锆合金在不同温度下的TSS[33,34]。锆合金的TSS通常有两种表示方法:一种是在升温过程中锆合金中所能溶解的最高氢含量,称为溶解极限固溶度(TSSD);另一种是在降温过程中氢化物恰好不析出的最大溶解度,称为析出极限固溶度(TSSP)。Motta等[35]总结了以上两种条件下测得的极限固溶度数据。如图4所示,TSSD固溶曲线和TSSP固溶曲线并不重合,相同温度下,TSSP固溶度高于TSSD固溶度。造成TSSD和TSSP差异的原因主要是析出的氢化物与锆基体的晶格失配度较大。降温过程中,氢原子从基体中脱溶导致晶格应变能的降低并不足以抵消氢化物析出造成晶格应变能的增加,因此氢化物析出受到抑制,氢在锆中过饱和溶解,导致TSSP高于平衡条件下的TSS;在升温过程中,氢化物发生溶解,由于已经析出的氢化物在其周围造成了塑性变形,消耗了由于氢化物析出造成的多余晶格应变能,降低了氢化锆溶解进锆基体的倾向,使TSSD低于TSSP[36]

图4

图4   锆合金TSSP及TSSD汇总图[35]

Fig.4   Summary of TSSP and TSSD for Zirconium and its alloys[35]


对于Zr-2.5Nb等双相合金,其TSS高于纯锆,但其对β-Zr相的分解[37,38]和合金经历的热循环过程[39]较为敏感。研究结果表明,β-Zr相分解程度的增加会导致双相锆合金TSSD和TSSP的下降[36]。反复升温、降温使氢化物析出具有“记忆效应”,随着热循环次数的增加,氢化物更易析出,因而也会导致TSSP的下降[33]

3 DHC行为机理及模型

图5为DHC裂纹扩展的基本过程示意图。如5a所示,氢原子在应力诱导作用下向高应力区(裂纹尖端或缺口前端)扩散并富集;当裂纹尖端氢浓度高于氢的固溶度时,氢化物形核析出并长大(图5b);位于裂尖的氢化物在应力集中的作用下断裂,且裂纹与原有的裂纹尖端相连,引起原有裂纹扩展一小段距离后止裂(图5c);随后氢原子向新的裂纹尖端扩散,并形成新的氢化物(图5d)。以上过程不断重复进行就导致了DHC裂纹的扩展[40]。目前,不同研究者对上述的DHC裂纹扩展过程基本认同,但对氢原子向高应力区扩散的根本驱动力还存有争议。国内外研究者针对DHC裂纹扩展机理建立了多种模型,其中比较具有代表性的是Dutton-Puls模型[41,42]和Kim's模型[36],下面就两个模型进行讨论。

图5

图5   DHC裂纹扩展过程示意图[40]

Fig.5   Diagrams of the propagation of cracks during DHC[40]


3.1 Dutton-Plus模型

Dutton-Puls模型(D-P模型)认为发生DHC的第一步是氢原子向裂纹尖端等应力集中区域的扩散[43]。拉应力使锆的晶格常数增加,更易容纳氢原子,即拉应力的存在使氢原子在锆基体中的化学势降低[44]。在应力集中的作用下,裂纹尖端受拉应力最大,因而氢原子在裂纹尖端的化学势最低,相当于在裂纹尖端形成“氢阱”[43,45]。位于低应力区的氢原子在化学势梯度的作用下向裂纹尖端扩散,造成氢原子不断在裂纹尖端富集,直至达到平衡状态。由上述分析可知,氢原子扩散的直接原因是化学势梯度,而化学势梯度是由应力梯度造成的,因而D-P模型中认为氢原子扩散的根本驱动力是应力梯度。如果达到平衡时裂纹尖端的氢含量高于TSSP,则在氢扩散至平衡状态前,裂纹尖端处会因氢含量达到TSSP而析出氢化物,从而导致裂纹尖端的溶解氢浓度降低;由于未达到平衡状态,氢原子持续向裂纹尖端扩散,从而导致裂尖区氢化物的生长,生长至临界尺寸后发生断裂。反之,若达到平衡所需的氢含量低于TSSP,则裂纹尖端无氢化物析出,DHC现象也就不会发生[46]

为模拟H原子向高应力区扩散的过程,D-P理论模型中建立圆柱几何模型[7],并且提出DHCR基本上等同于裂纹尖端片状氢化物的生长速率的假设,这种假设使DHCR与进入高应力区的氢扩散通量建立了直接联系[43]

D-P模型中首先提出两个假设:(1) 氢原子稳态扩散,即▽J≡0;(2) 扩散到高应力区的氢原子全部都生成氢化物。基于以上假设,推导并得到了稳态条件下的氢扩散通量表达 式(2)[41],以及DHCR表达 式(4)[42,44]

JH=-DHrΦ(ζ, L)ΩZrCHSEL-Eζ
Φ(ζ, L)=ζLr-1expϕ(r)RTdr
DHCR=dadr=2πDHΩZrΦ(ζ, l)thydNHx(EL-Eζ)
EL=CHSexp-p(L)V¯HhRT
Eζ=CHSexp-p(ζ)V¯HhRT

其中,DH = α-Zr中氢扩散常数(m2/s);CHS =氢在α-Zr中的平衡固溶度(mg/kg);ΩZr = α-Zr中的Zr的原子体积(m3);Φ(ζ, L) =关于应力势能的积分常数;ϕ(r) =距裂纹尖端距离为r处的应力势能;a =微裂纹长度(m);t =时间(s);thyd =氢化物厚度(m);NH =氢化物中的原子密度(atoms/m3);x =每个氢化物分子中的氢原子数;V¯Hh =氢原子在氢化物中的摩尔体积(m3/mol);p(r) =静水应力(MPa)。

从表达 式(4)可知,当E(L) > E(ζ)时,DHC现象才会发生。其中,E(L)和E(ζ)分别为r = Lr = ζ处的表观氢浓度,其表达式如式(5)、(6)所示。从(5)、(6)可知,r = Lr = ζ处的静水应力差造成E(L)和E(ζ)差异,与D-P模型中提出的观点相一致,即应力梯度是氢原子向裂纹尖端扩散的驱动力。

图3可知,DHC裂纹扩展分为3个阶段。其中在Stage I和Stage III阶段,DHCR与KI呈明显正相关关系,而在Stage II阶段,KI对DHCR无明显影响,这与D-P模型认为的应力梯度是DHC的驱动力明显不符。有研究中将不同的K值带入D-P模型中进行拟合[47]。拟合结果表明,K值对最终的DHCR拟合值的影响不大,且带入不同K值后得到的拟合值均与实验值相近,但这只能说明以上表达式具有一定的准确性,但仍不能作为应力梯度是DHC发生的驱动力的直接证据。

3.2 Kim's模型

Kim's模型认为DHC现象发生的第一步是氢化物的析出[43],DHC裂纹扩展的驱动力是裂纹尖端处氢化物析出后导致裂纹尖端与其他位置的氢浓度差ΔC[36]。Kim等用TSSD和TSSP曲线对该模型进行了解释,如图6所示,当温度升高至300 ℃时,试样中的溶解氢含量处于A点,此时裂纹尖端和基体没有浓度差。当再降温至250 ℃时,由于TSS的“滞后效应”,基体的溶解氢含量会降至B点,但此时基体中的溶解氢含量未达到析出值。裂纹尖端处应力较大,诱发了氢化物的析出(Kim's模型认为[36],应力集中使裂纹尖端产生塑性变形,能降低氢化物析出产生的晶格应变能,使氢化物更易析出),溶解氢含量降至C点。至此,裂纹尖端和基体产生了氢的浓度差ΔC。为证明DHC裂纹扩展的驱动力是ΔC而非D-P模型中的应力梯度,研究中通过将试样加热到不同温度后再降温的方法获得具有不同ΔC的试样,并对这些试样施以相同的应力强度因子。结果表明,在应力强度因子相同的情况下,ΔC越大,试样的DHCR越大,ΔC可达到的最大值为测试温度下的TSSP减去TSSD[36]

图6

图6   基于TSS曲线的热循环过程中氢含量变化示意图[36]

Fig.6   Variations of hydrogen concentration with temperature during thermal cycling, obtained based on TSS curves[36]


基于以上分析,Kim还提出了DHCR计算公式:

v=kJH=kDHCX=kDHC(KIH/σYS)2

其中,k =与材料有关的系数;JH =氢扩散通量(kg/(m2/s));ΔC =基体和裂纹尖端的氢浓度差(mg/kg);DH =氢扩散系数(m2/s);ΔX =裂纹尖端塑性区尺寸(m);KIH =裂纹扩展临界应力强度因子(MPa·m-2);σYS =屈服强度(MPa)。

在Kim's模型的所有版本中,均强调基体处的氢含量高于裂纹尖端的氢含量,原因是应力作用下裂纹尖端的TSSP下降,使原本处于TSSP以下的实际氢含量(即图6中的B点)超过了应力作用后的TSSP,导致氢化物析出[36,43],进而导致裂纹尖端氢含量下降。但目前还没有实验结果能直接证明应力作用会使裂纹尖端的TSSP下降[43]。有研究测定了氢在锆合金中的偏摩尔体积,表明应力对溶解度的影响很小[48,49],无应力状态下测得的溶解曲线可以直接应用在裂纹尖端,这是Kim's模型存在的争议点。

3.3 其它模型

不同于以上两种模型对DHC过程进行大量简化和假设,Varias 模型将含有氢化物的锆合金视为弹性复合材料,引入连续的扩散-热力学-机械耦合场,模拟DHC过程中真实的应力和温度[50,51]。Jernkvist等在Varias模型中引入氢化物相转变动力学,以更精确的描述氢化物的析出和溶解行为[52,53]。Shmakov模型采用了与D-P模型相似的物理原理,不同的是该模型认为氢原子扩散是非稳态的,并用分步数值计算的方法获得氢原子扩散通量[54]。Shmakov模型能较好的预测DHCR随KI变化的3个阶段和接近临界温度时DHCR降低的现象。上述模型各具特色,但仍难以统一共识。随着模拟计算能力和微观分析检测能力的不断提升,相信未来的模型会越来越准确。

4 不同因素对DHC行为的影响

4.1 材料微观结构和力学性能的影响

压力管在服役过程中受到高温和辐照的双重影响,其微观结构和力学性能会发生变化,进而对DHC行为产生影响。

前文已经提到,Zr-2.5Nb合金中含有α-Zr相和β-Zr相两相,α-Zr相晶粒沿压力管轴向被拉长,β-Zr相呈连续或半连续状态沿轴向分布在α-Zr相晶粒之间[14]。由于微观结构的各向异性,DHC行为也表现出各向异性。Kim等[14]研究表明氢在β-Zr相中的扩散速率较快,当裂纹扩展方向与β相分布方向一致时,会形成氢扩散的快速通道,裂纹起裂要求较低,因此轴向DHCR是径向DHCR的1.8~2倍,轴向KIH也高于径向KIH。研究结果表明,400 ℃热处理会使原本在轴向上连续分布的β-Zr相发生分解,热处理时间越长,β-Zr相分解程度越高[55]。随着热处理时间延长,轴向和径向的DHCR均降低。Jovanovic等[12]认为β-Zr相的分解降低了氢原子在压力管轴向的扩散速率,这是热处理导致轴向DHCR降低的主要原因;此外,材料的屈服强度随着热处理时长的增加而降低,有效缓解了裂纹尖端的应力集中程度,降低了裂纹尖端氢化物的断裂倾向,一定程度上抑制了DHC裂纹的萌生和扩展,这是导致轴向DHCR随热处理时长的延长而降低的另一个原因。对于径向DHCR,轴向氢扩散速率的下降对其基本没有影响,导致径向DHCR降低的主要原因是屈服强度的降低,因此径向DHCR下降幅度也小于轴向DHCR,且随着热处理时间的延长,轴向DHCR与径向DHCR的比值减小[55],压力管DHC行为的各向异性程度降低。除热处理外,辐照也能显著提高材料的屈服强度,从而导致较高的DHCR,使代表氢化物断裂尺寸的断口条纹宽度降低[56]

4.2 KI 的影响

KI对DHC行为的影响可根据KIKIHKIC相对关系大致分为3个阶段(图3):在Stage I,KI略大于KIH,DHCR随KI的增加而急剧上升;Stage I的KI范围较小,随后进入DHCR对KI不敏感的Stage II;当KI继续增加至KIC附近时,DHCR再次随着KI的增加而急剧增长,即进入Stage III。在Stage III发生裂纹失稳扩展的现象较易理解,但对Stage I和Stage II所呈现的规律的机理以及这两个阶段的具体的分界点仍不清晰。

Yan等[57,58]重点关注了KIKIH附近时氢化物的生长及断裂情况。结果表明,当KI < KIH时,裂纹尖端的氢化物尺寸随KI的增加而增加且不发生断裂;当KI > KIH时,氢化物断裂,其断裂尺寸随KI的增加呈先降低后趋于稳定的趋势。氢化物的生长是一个自限性过程,随着氢化物尺寸的增加,氢化物生长的驱动力逐渐减弱,因此在达到临界断裂尺寸前,氢化物就已经停止生长,即DHC不发生。随着KI的增加,裂纹尖端氢化物所受的应力增加,这会导致氢化物更易断裂,即临界断裂尺寸减小;裂纹尖端塑性区尺寸随KI增加,这为氢化物提供了更大的生长空间和驱动力。因此,在KI增加的过程中,必然存在一个临界点,使氢化物的实际尺寸等于临界断裂尺寸,该临界点对应的KIKIH[41]。在KI略大于KIH的小范围内,KI的增加使氢化物临界断裂尺寸继续减小,氢化物更易断裂;并且KI的增加促进了氢原子向裂纹尖端扩散和氢化物的生长速率,因此在该范围内DHCR随KI的增加而急剧升高(Stage I)。

Shi等[46]认为应力改变了裂纹尖端氢的平衡固溶度,并在圆柱几何模型的基础上提出了裂纹尖端及其附近的应力分布(图7a)及氢浓度分布(图7b)示意图,并针对DHCR在Stage II的范围内对KI无明显依赖性的现象做出了解释。

图7

图7   裂纹尖端静水应力分布和稳态扩散条件下裂纹尖端氢浓度分布示意图[46]

Fig.7   Stress distribution at the crack tip (a) and hydrogen concentration distribution at the crack tip under steady state diffusion condition (b)[46]


图中,C0为无应力区原始氢浓度,rp为静水应力开始降低的转折点(在裂纹尖端附近),L ≫ rpC(L)为r=L处的氢浓度,Cmax为应力作用下裂纹尖端所能达到的最高氢浓度(并非裂纹尖端实际的氢浓度,若Cmax > TSSP,则在实际氢浓度到达Cmax前就有氢化物析出,导致氢含量降低)。从图中可知,当r < rp时,即在裂纹尖端及其附近的范围内,静水应力保持恒定,氢浓度也保持恒定。C0C(L)以及Cmax的关系,如 式(8)和(9)所示,

C(L)C01+13Φ(L)
Cmax=βC0
β=exp(Φ(rp)-Φ(L))1+13Φ(L)
Φ(r)=V¯HRTP(r)

其中,Φ(r)为应力势能;V¯H为氢溶解在锆合金中的摩尔体积(m3/mol);P(r)为距裂尖r处的静水应力(MPa);

裂纹足够尖锐的条件下,塑性区内P(r)的表达式如下:

P(r)=2(1+v)KI32πr

对于P(rp),可以根据HRR原理计算[59,60],但在锆合金中,P(rp)约等于P(rpz)[42],(rpz为塑性区尺寸)

r(pz)=(1-2v)22πKIσY2

其中,σy为屈服强度(MPa);KI为应力强度因子。

式(12)带入 式(11)则有

P(rp)P(rpz)=2(1+v)σY3(1-2v)

式(13)带入 式(11)可得

Φ(rp)=V¯HRT2(1+v)σY3(1-2v)

式(15)可知,裂纹尖端的应力势能Φ(rp)与施加的应力强度因子无关。且LrpP(L) ≪ P(rp),r =L处的应力势能Φ(L)可以忽略,因此CmaxC0的比值β也与应力强度因子无关,即裂纹尖端所能达到的最高氢固溶度不受KI的影响。综上所述,KI的变化无法改变r = rpr = L处的应力势能差(Φ(rp)-Φ(L))和氢浓度差(Cmax-C(L)),所以KI的变化对氢扩散通量的影响很小。根据D-P模型的假设可知,扩散进入高应力区的氢原子全部生成氢化物,且氢化物在X轴方向上的生长速率等同于DHCR。氢扩散通量直接决定了氢原子的扩散速度和DHCR,所以在Stage II范围内KI的变化对DHCR影响不大。

4.3 温度的影响

DHC裂纹扩展是一个热激活过程[41],其中,氢原子扩散、氢化物析出和氢化物断裂等均受到温度的影响。虽然压力管的实际服役温度较为稳定,但停堆检修期间涉及到温度的变化,因此明确温度对DHC行为的影响有助于寻找降低压力管泄漏风险的措施。

对于固定的氢浓度,温度过高会使氢浓度低于或略高于TSSP,从而导致氢化物不析出或析出困难,因此DHC现象的发生存在临界温度。Sagat等发现当温度升高至350 ℃时,压力管不产生DHC现象[61]

在DHC现象能发生的温度范围内,DHCR通常是随着温度的升高而增大。许多研究者[61~63]在100~300 ℃的范围内测量了Zr-2.5Nb压力管轴向的DHCR,结果表明,该温度范围内DHCR和温度呈现出正相关,并符合Arrhenius规律[63]。随着温度的升高,断口的条纹特征也变得更加明显[13],条纹宽度也增大[56]。根据Kim's模型,DHCR主要由氢扩散速率和氢浓度梯度决定,其中氢扩散速率占主导地位;当温度增大到300 ℃以上时,氢扩散速率不再是DHCR的决定性因素,DHCR会随着温度的升高而降低[61],DHCR和温度之间不符合Arrhenius规律。Sunil等[64]认为,DHCR的降低和存在DHC裂纹扩展临界温度的原因是温度升高使材料发生蠕变,导致裂纹尖端的应力松弛,使氢化物的断裂变得困难;也有研究认为[65],在加载的过程中只有锆合金基体的屈服强度大于氢化物的断裂应力,氢化物才会发生断裂,而随着温度的升高,基体的屈服强度降幅超过了氢化物断裂应力的降幅,使氢化物的断裂变得困难。

Holston和Stjärnsäter[66]采用降载荷法测量了不同温度下Zr-4包壳管DHC的KIH,研究表明随着温度的升高,KIH值缓慢升高,但当温度超过300 ℃时KIH会急剧增加,Kim等[67]也观察到了类似的现象。这与上述关于300 ℃以上DHCR降低的现象相符合,即KIH越高,DHC裂纹扩展越困难。但同样采用降载荷法测量KIH值,Sagat等[68]研究表明,KIH值与温度的关联性并不大。Kim等[67]认为,裂纹尖端塑性变形区中产生的变形孪晶与氢化物的相互作用使氢化物断裂。同时,Kim还分别采用降载荷法和升载荷法获得了300 ℃以下不同温度的KIH值,发现使用降载荷法时,裂纹前端的塑性变形区较大,变形孪晶数量多[16]。虽然在300 ℃以下时,变形孪晶数量随温度的升高而增加,但在较低温度时变形孪晶数量并不是氢化物断裂的限制条件,所以KIH值与温度关联性不大;使用升载荷法时,裂纹前端塑性变形区小,变形孪晶数量较少,因此成为氢化物断裂的限制条件,温度越高,变形孪晶数量越多,氢化物越容易断裂,KIH值随温度的升高而降低[16]。当温度在300 ℃以上时,变形孪晶数量随温度的升高而减少,导致氢化物不容易发生断裂,使KIH急剧上升,也导致了DHCR的下降。

4.4 热循环过程的影响

在重水堆检修期间,压力管道不可避免的会经历升温或降温过程。因此,明确热循环过程对DHC行为的影响有助于改进检修方案,最大限度的降低DHC裂纹扩展对压力管的损伤。

升温、降温过程的重复会使氢化物反复溶解和析出。当氢化物第一次析出后会使析出周围的基体发生塑性变形,因而造成氢化物析出的“记忆效应”,再次降温时氢化物会更易析出,导致DHC发生的临界温度上升[69,70]。DHCR受裂纹尖端(高应力区)和低应力区的氢浓度梯度影响,但由于溶解氢含量受温度变化过程影响较大,材料经历的热循环过程会导致其中的氢含量产生差异,进而造成DHCR的不同。实验结果表明,升温至实验温度时所测的DHCR明显低于降温至实验温度时所测的DHCR[47]。对于具有相同氢含量的试样,如果升温到实验温度,则DHCR发生转折的温度比降温至实验温度时所测的转折温度低[18,71]。Shi等[46]利用D-P模型预测了不同热循环条件下所得的临界温度,结果表明,降温至实验温度条件下计算得到临界温度取决于氢含量,而升温至实验温度条件下计算得到的临界温度不受氢含量影响[46],该预测结果得到了实验数据的验证。

4.5 H含量的影响

压力管在服役过程中与重水持续发生腐蚀反应并产生氢(氘),因此压力管内部的氢(氘)含量随着服役时间的延长而增加;在同一时间,压力管的不同位置的氢(氘)含量也会有所差异[72,73]。因此,探究氢含量对压力管DHC行为的影响对延长压力管使用寿命具有重要的意义。

研究认为,在一定温度下,裂纹尖端DHC行为发生的临界氢含量与TSSD接近,因此TSSD是裂纹尖端发生DHC现象的氢含量阈值[74]。Kim等[75]用升K法和降K法分别测量了280 ℃时不同氢含量试样的KIH值,结果表明,当实际氢含量高于该温度下的TSSD,且差值ΔCC =实际氢含量-TSSD)在25 mg/kg以下时,两种方式测得的KIH值均随氢含量的增加而降低;当ΔC超过25 mg/kg时,KIH值保持稳定[75]。Kim等[67]在随后的研究中还观察到氢化物尺寸越大,越容易与变形孪晶相互作用导致氢化物断裂,即ΔC越大,形成的氢化物尺寸越大,DHC越容易发生,这解释了KIH对氢含量的依赖性。同样,DHCR也是先随着氢含量的增加而增大,但是当氢含量增加至35 mg/kg后,DHCR趋于稳定[67]。Kim's模型[36]中DHC裂纹扩展的驱动力是裂纹尖端处与其他位置的溶解氢的浓度差,而裂尖处溶解氢含量为TSSD,因此在一定温度下,不同位置的溶解氢的差值最大为TSSP减去TSSD,所以当氢浓度差超过TSSP-TSSD的值时,氢含量的变化不会对DHCR和KIH产生影响。然而,也有研究发现当温度为250 ℃,氢含量在40~60 mg/kg的范围内,DHCR随氢含量的升高而增加[76]

5 总结与展望

针对锆合金在高温环境中的延迟氢化物开裂行为已开展了大量的研究,为压力管等核电关键结构的锆合金部件的安全运行和老化管理提供了大量理论和数据支撑。但是,在研究中仍存在以下不足:(1) 由于锆合金的延迟氢化物开裂行为受温度、氢含量、受力和材料微观结构等因素的影响,其发生机理存在争议,氢化物的形核、生长和断裂过程以及影响氢化物断裂尺寸的因素仍不明确,对于不同因素间的耦合作用的影响有待进一步探索。(2) 氢含量对DHCR的影响趋势尚存在争议。(3) 普遍的锆合金DHC行为测试方法较为单一,且测试过程耗时长,效率较低。(4) 缺乏能实时监测DHC裂纹扩展长度的方法。(5) 缺乏模拟核用锆合金实际服役工况的实验设备,尤其是中子辐照的影响,使实验测试结果可能与材料实际服役结果有较大的出入。

针对上述的不足之处,未来要继续深入研究不同因素对DHC行为的影响。结合压力管的实际服役工况,不断改进实验方法,开发更接近实际服役状态的实验设备,尤其是针对中子辐照样品的自动化操作设备及热室。随着研究的深入和计算机技术的发展,模拟计算 被越来越广泛的应用到DHC行为的研究中。将模拟计算的结果与实验结果相结合,会极大的提高实验效率,有助于进一步探索DHC机理;还可以模拟压力管实际服役环境,为降低压力管的失效风险和延长使用寿命提供数据支撑。

参考文献

Liao J P, Mao Y L, Jin D S, et al.

Laboratory simulation of crud deposition on Zr-alloy fuel cladding in simulated pressurized water reactor primary coolant

[J]. J. Chin. Soc. Corros. Prot., 2023, 43: 197

[本文引用: 1]

廖家鹏, 毛玉龙, 金德升 .

锆合金包壳在模拟压水堆一回路冷却剂中的表面污垢沉积行为研究

[J]. 中国腐蚀与防护学报, 2023, 43: 197

DOI      [本文引用: 1]

阐述了影响燃料包壳表面污垢沉积行为的主要因素,通过设计并开展针对国产锆合金燃料包壳的污垢沉积试验,成功制备出含多孔基体和烟囱结构主要元素组成的燃料包壳氧化物污垢,污垢的主要成分为NiFe<sub>2</sub>O<sub>4</sub>、Fe<sub>2</sub>O<sub>3</sub>和NiO,并在污垢内部检测到LiBO<sub>2</sub>的析出,实现了对PWR堆芯燃料包壳污垢沉积的实验室模拟。

Ru X, Staehle R W.

Historical experience providing bases for predicting corrosion and stress corrosion in emerging supercritical water nuclear technology: part 1—review

[J]. Corrosion, 2013, 69: 211

[本文引用: 1]

Grade A M.

Effects of irradiation and hydriding on the mechanical properties of zircaloy-4 at high fluence

[A]. Zirconium in the Nuclear Industry: Eighth International Symposium [M]. Philadelphia: ASTM, 1989

[本文引用: 1]

Daunys M, Dundulis R, Grybenas A, et al.

Hydrogen influence on mechanical and fracture mechanics characteristics of zirconium Zr-2.5Nb alloy at ambient and elevated temperatures

[J]. Nucl. Eng. Des., 2008, 238: 2536

[本文引用: 1]

Moan G D, Coleman C E, Price E G, et al.

Leak-before-break in the pressure tubes of CANDU reactors

[J]. Int. J. Pressure Vessels Piping, 1990, 43: 1

[本文引用: 1]

Wang M J, Qiu S Z, Su G H, et al.

Research on the leak-rate characteristics of leak-before-break (LBB) in pressurized water reactor (PWR)

[J]. Appl. Therm. Eng., 2014, 62: 133

[本文引用: 1]

Puls M P.

Delayed hydride cracking: theory and experiment

[A]. The Effect of Hydrogen and Hydrides on the Integrity of Zirconium Alloy Components. Engineering Materials [M]. London: Springer, 2012: 333

[本文引用: 2]

Shi S Q, Puls M P.

Criteria for fracture initiation at hydrides in zirconium alloys I. Sharp crack tip

[J]. J. Nucl. Mater., 1994, 208: 232

[本文引用: 1]

Liu Y Z, Zhao W J, Peng Q, et al.

Study on electrolytic hydrogen infiltration and determination of hydrogen content of Zr-Sn-Nb alloy

[A]. China Materials Seminar [C]. Beijing, 2002: 4

[本文引用: 1]

刘彦章, 赵文金, 彭 倩 .

Zr-Sn-Nb合金电解渗氢及氢含量确定研究

[A]. 中国材料研讨会 [C]. 北京, 2002: 4

[本文引用: 1]

Yao C, Chen J, Ming H L, et al.

Research progress on hydrogen permeability behavior of pipeline steel

[J]. J. Chin. Soc. Corros. Prot., 2023, 43: 209

姚 婵, 陈 健, 明洪亮 .

管线钢氢渗透行为的研究进展

[J]. 中国腐蚀与防护学报, 2023, 43: 209

Zhou X, Wu D K, Cheng X, et al.

Research progress of detection techniques for permeated hydrogen

[J]. J. Chin. Soc. Corros. Prot., 2023, 43: 1203

[本文引用: 1]

周 欣, 吴大康, 成 旭 .

渗透氢检测方法研究进展

[J]. 中国腐蚀与防护学报, 2023, 43: 1203

[本文引用: 1]

Mieza J I, Vigna G L, Domizzi G.

Evaluation of variables affecting crack propagation by delayed hydride cracking in Zr-2.5Nb with different heat treatments

[J]. J. Nucl. Mater., 2011, 411: 150

[本文引用: 2]

Pan C T, Zhao G N, Bao Y C, et al.

Effect of temperature on the delayed hydride cracking rate of Zr-2.5Nb alloy pressure tubes

[J]. J. Nucl. Mater., 2024, 588: 154778

[本文引用: 1]

Kim Y S, Cheong Y M.

Anisotropic delayed hydride cracking velocity of CANDU Zr-2.5Nb pressure tubes

[J]. J. Nucl. Mater., 2008, 373: 179

[本文引用: 4]

Sun C, Tan J, Ying S H, et al.

Threshold stress intensity factor for delayed hydride cracking of a recrystallized N18 alloy plate along the rolling direction

[J]. J. Nucl. Mater., 2010, 406: 212

[本文引用: 1]

Kim Y S, Park S S, Kwun S I.

Threshold stress intensity factor, KIH for delayed hydride cracking of a Zr-2.5Nb tube with loading mode

[J]. J. Alloy. Compd., 2008, 462: 367

[本文引用: 4]

Shmakov A A, Singh R N, Yan D, et al.

A combined SIF and temperature model of delayed hydride cracking in zirconium materials

[J]. Comput. Mater. Sci., 2007, 39: 237

[本文引用: 1]

Simpson L A, Puls M P.

The effects of stress, temperature and hydrogen content on hydride-induced crack growth in Zr-2.5 pct Nb

[J]. Metall. Trans., 1979, 10A: 1093

[本文引用: 2]

Coleman C E.

Simulating the behavior of zirconium-alloy components in nuclear reactors

[A]. MoanG D, RudlingP. Zirconium in the Nuclear Industry: Thirteenth International Symposium [M]. West Conshohocken: ASTM, 2002: 3

[本文引用: 2]

Lumley S C, Grimes R W, Murphy S T, et al.

The thermodynamics of hydride precipitation: the importance of entropy, enthalpy and disorder

[J]. Acta Mater., 2014, 79: 351

[本文引用: 1]

Ma Y, Lan Y N, Chen J W.

A novel cross-sectional metallography method for determining hydrogen absorption concentration and hydrogen absorption amount of Zr-Sn-Nb alloy cladding caused by high temperature water corrosion

[J]. J. Chin. Soc. Corros. Prot., 2024, 44: 261

[本文引用: 1]

马 雁, 蓝宇宁, 陈嘉威.

Zr-Sn-Nb包壳管腐蚀吸氢中氢浓度测算的截面金相法

[J]. 中国腐蚀与防护学报, 2024, 44: 261

[本文引用: 1]

Kim J S, Kim S D, Yoon J.

Hydride formation on deformation twin in zirconium alloy

[J]. J. Nucl. Mater., 2016, 482: 88

[本文引用: 1]

Une K, Nogita K, Ishimoto S, et al.

Crystallography of zirconium hydrides in recrystallized Zircaloy-2 fuel cladding by electron backscatter diffraction

[J]. J. Nucl. Sci. Technol., 2004, 41(7): 731

[本文引用: 1]

Nath B, Lorimer G W, Ridley N.

Effect of hydrogen concentration and cooling rate on hydride precipitation in α-zirconium

[J]. J. Nucl. Mater., 1975, 58(2): 153

[本文引用: 1]

Roy C, Jacques J G.

{1017} Hydride habit planes in single crystal zirconium

[J]. J. Nucl. Mater., 1969, 31: 233

[本文引用: 1]

Lin X H, Beyerlein I J, Han W Z.

Annealing cracking in Zr and a Zr-alloy with low hydrogen concentration

[J]. J. Mater. Sci. Technol., 2024, 182: 165

DOI      [本文引用: 1]

Hydrogen embrittlement, also known as hydrogen-induced cracking, is a longstanding concern for zirconium (Zr) and its alloys. The affinity for hydrogen in zirconium results in the formation of brittle hydrides. Hydrides are conventionally thought to facilitate crack propagation, and the prevailing recommendation to hinder hydride embrittlement involves reducing the hydrogen concentration in the material and applying slow cooling to reduce thermal stresses. Here, in Zr with a maximum hydrogen concentration of 43 wppm, we show that cracks initiate unexpectedly after annealing, a step that involves slow cooling and no external stress. With a suite of high-resolution microscopy, we show that annealing cracks initiate at prismatic &#x003b4;-hydrides, whose presence is unexpected. The microscopic annealing crack paths are consequently crystallographic, following the habit plane of these prismatic &#x003b4;-hydrides. The commonly reported &#x003b3;-hydride and &#x003b4;-hydrides that develop on the basal plane are found not to initiate cracks. We rationalize that prismatic &#x003b4; hydrides result from the intergranular thermal residual stresses that develop during slow cooling. This finding of uncommon annealing cracking further advances the mechanistic understanding of hydride embrittlement in zirconium and can lead to new strategies for design and fabrication of advanced zirconium alloys for harsh environments.

Maimaitiyili T, Steuwer A, Blomquist J, et al.

In-situ hydrogen charging of zirconium powder to study isothermal percipitation of hydrides and determination of Zr-hydride crystal structure

[A]. Proceedings of the 16th International Conference on Environmental Degradation of Materials in Nuclear Power Systems-Water Reactors [C]. Asheville, 2013: 8

[本文引用: 1]

Zhao Z, Blat-Yrieix M, Morniroli J P, et al.

Characterization of zirconium hydrides and phase field approach to a mesoscopic-scale modeling of their precipitation

[J]. J. ASTM Int., 2008, 5: 1

[本文引用: 1]

Perovic V, Weatherly G C, Simpson C J.

Hydride precipitation in α/β zirconium alloys

[A]. AshbyM F, HirthJ P. Perspectives in Hydrogen in Metals [M]. Oxford, New York: Pergamon, 1986: 469

[本文引用: 2]

Northwood D O, Gilbert R W.

Hydrides in zirconium-2.5wt.% niobium alloy pressure tubing

[J]. J. Nucl. Mater., 1978, 78: 112

[本文引用: 1]

Lee K W, Hong S I.

Zirconium hydrides and their effect on the circumferential mechanical properties of Zr-Sn-Fe-Nb tubes

[J]. J. Alloy. Compd., 2002, 346: 302

[本文引用: 1]

Xu C R, Zhao W J, Deng Z G, et al.

Review of research on stress reorientation of hydrides in zirconium alloy cladding tube

[J]. Hot Working Technol., 2016, 45(12): 19

[本文引用: 1]

徐春容, 赵文金, 邓治国 .

锆合金包壳管氢化物应力再取向研究概述

[J]. 热加工工艺, 2016, 45(12): 19

[本文引用: 1]

Parodi S A, Ponzoni L M E, De Las Heras M E, et al.

Study of variables that affect hydrogen solubility in α + β Zr-alloys

[J]. J. Nucl. Mater., 2016, 477: 305

[本文引用: 2]

Fang Q.

Characterization of hydrides and delayed hydride cracking in zirconium alloys

[D]. Ontario: Queen's University, 2016

[本文引用: 1]

Motta A T, Capolungo L, Chen L Q, et al.

Hydrogen in zirconium alloys: a review

[J]. J. Nucl. Mater., 2019, 518: 440

DOI      [本文引用: 3]

Hydrogen absorbed into zirconium alloy nuclear fuel cladding as a result of the waterside corrosion reaction can affect the properties of nuclear fuel, principally through the precipitation of brittle hydride particles. Multiple phenomena are involved in this overall process: after hydrogen pickup degradation of mechanical properties is controlled by hydrogen transport, hydride precipitation and dissolution kinetics and the formation of specific mesoscale hydride microstructures. The precipitation of hydrides especially affects cladding ductility and fracture toughness, but can also affect other phenomena, including via stress-induced hydride reorientation. These processes can affect cladding performance both during normal operation and during extended dry storage, as hydride morphology can be modified during the preparatory vacuum drying processes. We review the processes of hydrogen transport, hydride precipitation and dissolution and formation of mesoscale hydride microstructures, and highlight where more research is needed, both from an experimental and from a modeling point of view. Crown Copyright (C) 2019 Published by Elsevier B.V.

Kim Y S.

Driving force for delayed hydride cracking of zirconium alloys

[J]. Met. Mater. Int., 2005, 11: 29

[本文引用: 10]

Khatamian D.

Solubility and partitioning of hydrogen in metastable Zr-based alloys used in the nuclear industry

[J]. J. Alloy. Compd., 1999, 293-295: 893

[本文引用: 1]

Khatamian D.

Effect of β-Zr decomposition on the solubility limits for H in Zr-2.5Nb

[J]. J. Alloy. Compd., 2003, 356-357: 22

[本文引用: 1]

Pan Z L, Ritchie I G, Puls M P.

The terminal solid solubility of hydrogen and deuterium in Zr-2.5Nb alloys

[J]. J. Nucl. Mater., 1996, 228: 227

[本文引用: 1]

Chu W Y, Qiao L J, Li J X. Hydrogen Embrittlement and Stress Corrosion Cracking: the Base Component [M]. Beijing: Science Press, 2013: 174

[本文引用: 3]

褚武扬, 乔利杰, 李金许. 氢脆和应力腐蚀: 基础部分 [M]. 北京: 科学出版社, 2013: 174

[本文引用: 3]

Dutton R, Nuttall K, Puls M P, et al.

Mechanisms of hydrogen induced delayed cracking in hydride forming materials

[J]. Metall. Trans., 1977, 8A: 1553

[本文引用: 4]

Puls M P.

Effects of crack tip stress states and hydride-matrix interaction stresses on delayed hydride cracking

[J]. Metall. Trans., 1990, 21A: 2905

[本文引用: 3]

McRae G A, Coleman C E, Leitch B W.

The first step for delayed hydride cracking in zirconium alloys

[J]. J. Nucl. Mater., 2010, 396: 130

[本文引用: 6]

Puls M P.

Review of the thermodynamic basis for models of delayed hydride cracking rate in zirconium alloys

[J]. J. Nucl. Mater., 2009, 393: 350

[本文引用: 2]

Thompson A W, Bernstein I M.

Effect of hydrogen on behavior of materials

[A]. Proceedings of Conference on effects of hydrogen on behavior of materials [C]. Moran, 1976: 717

[本文引用: 1]

Shi S Q, Shek G K, Puls M P.

Hydrogen concentration limit and critical temperatures for delayed hydride cracking in zirconium alloys

[J]. J. Nucl. Mater., 1995, 218: 189

[本文引用: 6]

De Las Heras M E, Parodi S A, Ponzoni L M E, et al.

Effect of thermal cycles on delayed hydride cracking in Zr-2.5Nb alloy

[J]. J. Nucl. Mater., 2018, 509: 600

[本文引用: 2]

MacEwen S R, Coleman C E, Ells C E, et al.

Dilation of h.c.p. zirconium by interstitial deuterium

[J]. Acta Metall., 1985, 33: 753

[本文引用: 1]

Eadie R L, Coleman C E.

Effect of stress on hydride precipitation in zirconium-2.5% niobium and on delayed hydride cracking

[J]. Scr. Metall., 1989, 23: 1865

[本文引用: 1]

Varias A G, Massih A R.

Simulation of hydrogen embrittlement in zirconium alloys under stress and temperature gradients

[J]. J. Nucl. Mater., 2000, 279: 273

[本文引用: 1]

Feng J L, Varias A G, Sui Y K.

Finite element analysis for steady-state hydride-induced fracture in metals by composite model

[J]. Int. J. Solids Struct., 2006, 43: 2174

[本文引用: 1]

Jernkvist L O, Massih A R.

Multi-field modelling of hydride forming metals. Part I: model formulation and validation

[J]. Comput. Mater. Sci., 2014, 85: 363

[本文引用: 1]

Jernkvist L O.

Multi-field modelling of hydride forming metals part II: application to fracture

[J]. Comput. Mater. Sci., 2014, 85: 383

[本文引用: 1]

Shmakov A A, Kalin B A, Ioltukhovskii A G.

A theoretical study of the kinetics of hydride cracking in zirconium alloys

[J]. Met. Sci. Heat Treat., 2003, 45: 315

[本文引用: 1]

Jovanović M T, Eadie R L, Ma Y, et al.

The effect of annealing on hardness, microstructure and delayed hydride cracking in Zr-2.5Nb pressure tube material

[J]. Mater. Charact., 2001, 47: 259

[本文引用: 2]

Shah P K, Dubey J S, Kumar A, et al.

Delayed hydride crack growth study on irradiated Zr-2.5Nb pressure tube

[J]. J. Nucl. Mater., 2015, 460: 1

[本文引用: 2]

Yan D, Eadie R L.

An approach to explain the stage I/II behaviour of the delayed hydride cracking velocity vs. KI curve for Zr-2.5Nb

[J]. Scr. Mater., 2000, 43: 89

[本文引用: 1]

Yan D, Eadie R L.

The critical length of the hydride cluster in delayed hydride cracking of Zr-2.5wt%Nb

[J]. J. Mater. Sci., 2000, 35: 5667

[本文引用: 1]

Rice J R, Rosengren G F.

Plane strain deformation near a crack tip in a power-law hardening material

[J]. J. Mech. Phys. Solids, 1968, 16: 1

[本文引用: 1]

Hutchinson J W.

Singular behaviour at the end of a tensile crack in a hardening material

[J]. J. Mech. Phys. Solids, 1968, 16: 13

[本文引用: 1]

Sagat S, Puls M P.

Temperature limit for delayed hydride cracking in Zr-2.5Nb alloys

[A]. Proceedings of the 17th International Conference on Structural Mechanics in Reactor Technology [C]. Prague, 2003: 17

[本文引用: 3]

Markelov V A, Kotov P V, Zheltkovskaya T N.

Temperature dependence of the velocity of delayed hydride cracking in Zr-2.5%Nb alloy

[J]. Inorg. Mater. Appl. Res., 2010, 1: 217

Kim Y S.

Temperature dependency of delayed hydride cracking velocity in Zr-2.5Nb tubes

[J]. Mater. Sci. Eng., 2007, 468-470: 281

[本文引用: 2]

Sunil S, Bind A K, Khandelwal H K, et al.

Delayed hydride cracking behavior of Zr-2.5Nb alloy pressure tubes for phwr700

[J]. J. Nucl. Mater., 2015, 466: 208

[本文引用: 1]

Shi S Q, Puls M P.

Fracture strength of hydride precipitates in Zr-2.5Nb alloys

[J]. J. Nucl. Mater., 1999, 275: 312

[本文引用: 1]

Holston A M A, Stjärnsäter J.

On the effect of temperature on the threshold stress intensity factor of delayed hydride cracking in light water reactor fuel cladding

[J]. Nucl. Eng. Technol., 2017, 49: 663

[本文引用: 1]

Kim Y S, Ahn S B, Kim K S, et al.

Temperature dependence of threshold stress intensity factor, KIH in Zr-2.5Nb alloy and its effect on temperature limit for delayed hydride cracking

[J]. Key Eng. Mater., 2006, 326-328: 919

[本文引用: 4]

Sagat S, Coleman C E, Griffiths M, et al.

The effect of fluence and irradiation temperature on delayed hydride cracking in Zr-2

.5Nb [A]. Zirconium in the Nuclear Industry: Tenth International Symposium [M]. Philadelphia: ASTM, 1994: 35

[本文引用: 1]

Sun C, Tan J, Ying S H, et al.

Study of the critical temperatures for delayed hydride cracking in N18 zirconium alloy

[J]. Acta Metall. Sin., 2009, 45: 541

[本文引用: 1]

<p>Zirconium alloys are used extensively in nuclear reactor cores. During their service a part of hydrogen produced through the corrosion reaction of Zr with hot coolant is absorbed by materials. Hydride induced embrittlement significantly influences the in&ndash;service performance of the Zr&ndash;alloy components. Delayed hydride cracking (DHC) is a localized form of hydride embrittlement, consequently, hydrogen atoms in the solid solution will diffuse into this region ahead of the crack tip subjected to a triaxial state of stress, which may lower the chemical potential of the region. Once the hydrogen concentration in this region reaches the terminal solid solubility (TSS), hydrides will start to form and grow. When the hydrides at the crack tip reach a critical size, the main crack will propagate through this hydrided region. The crack front finally is arrested at the end of the hydrided region by the ductile zirconium matrix, and the whole process repeats itself.<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Most of the investigations on DHC in zirconium alloys are focused on Zr&ndash;Nb alloys. Few literatures were found on the subject of DHC in Zr&ndash;Sn alloys. The purpose of the present study was to investigate critical temperature for initiating and arresting delayed hydride cracking in Zr&ndash;Sn&ndash;Nb alloy.<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A critical temperature for DHC study was carried out to determine the critical temperature for initiating and arresting in N18 zirconium alloy (Zr&ndash;Sn&ndash;Nb alloy). For a given hydrogen concentration of a specimen, the two critical temperatures were observed&mdash;a DHC initiation temperature, <em>T</em><sub>c</sub>, at which DHC would initiate when approaching the test temperature from above the terminal solid solubility (Cd) temperature in hydride dissolution and a DHC arrest temperature, Th, obtained by heating the same specimen from <em>T</em><sub>c</sub> after DHC had started. <em>T</em><sub>c</sub> slightly below Th. Both <em>T</em><sub>c</sub> and <em>T</em><sub>h</sub> fall below the dissolution solvus temperature and above the precipitation solvus temperature. A theoretical analysis was carried out to quantitatively determine the hydrogen concentration limit and these critical temperatures using the method of Dutton and Plus, a key assumption in the method is that, while the local crack tip stress concentration causes a local enhancement of the hydrogen concentration in solution, the hydride precipitation solvus is unaffected by stress. Good agreements are obtained between measured and predicted values of critical temperatures, which support the Dutton--Plus theory.</p>

孙 超, 谭 军, 应诗浩 .

N18锆合金氢致裂纹延迟开裂临界温度研究

[J]. 金属学报, 2009, 45: 541

[本文引用: 1]

Sun C, Tan J, Ying S H, et al.

Prediction of critical temperature for delayed hydride cracking in irradiated N18 zirconium alloy

[J]. Acta Metall. Sin., 2010, 46: 805

DOI      [本文引用: 1]

Delayed hydride cracking (DHC) has been recognized as a potential threat to the structural integrity of zirconium alloy components in water-cooled nuclear reactors since the early 1970 s. Although DHC has been studied for a long time, most of the investigations were focused on unirradiated Zr-Nb and Zr-Sn alloys, no information was found on Zr-Sn-Nb alloy. Currently, several new Zr-based alloys have been developed from the point of view of enhancement of corrosion resistance and hydrogen pickup properties. In China, N18 (Zr-Sn-Nb) alloy has been developed based on Zr-Sn and Zr-Nb systems to meet the requirements of higher fuel burn-up. Therefore, it is necessary to investigate the DHC behavior of N18 alloy systematically. The objective of this study is to predict the critical temperatures for initiating and arresting DHC in irradiated N18 alloy. A model used to forecast DHC critical temperature of irradiated zirconium alloys was established based on the terminal solid solubility of hydrogen in irradiated zirconium alloys and stress-induced hydrogen diffusion. By using the model, the critical temperatures for irradiated N18 alloy were predicted. The result showed that the irradiated N18 alloy has high sensitivity for DHC initiation, because neutron irradiation increased the solubility of hydrogen and yield strength of N18 alloy. Compared with unirradiated N18 alloy, the critical temperature of irradiated N18 alloy rise nearly 20℃. The variation of critical temperatures of irradiated N18 alloy is similar to the unirradiated one. Tc is slightly less than Th.

孙 超, 谭 军, 应诗浩 .

辐照后N18锆合金氢致延迟开裂临界温度预测

[J]. 金属学报, 2010, 46: 805

[本文引用: 1]

Cheadle B A, Coleman C E, Ambler J F R.

Prevention of delayed hydride cracking in zirconium alloys

[A]. Zirconium in the Nuclear Industry [M]. Strasbourg: ASTM, 1987: 846

[本文引用: 1]

Choo K N, Kim Y S.

Hydrogen uptake and corrosion behavior of Zr-2.5Nb pressure tubes in wolsong unit 1

[J]. J. Nucl. Mater., 2001, 297: 52

[本文引用: 1]

Bao Y C, Shi X Q, Zhao C L.

Hydrogen corrosion-uptake analysis and modeling for heavy water reactor Zr-2.5Nb pressure tubes

[J]. Corros. Prot., 2020, 41(11): 22

[本文引用: 1]

鲍一晨, 石秀强, 赵传礼.

重水堆Zr-2.5Nb压力管腐蚀吸氢分析与建模

[J]. 腐蚀与防护, 2020, 41(11): 22

[本文引用: 1]

Coleman C E, Ambler J F R.

Solubility of hydrogen isotopes in stressed hydride-forming metals

[J]. Scripta Metallurgica, 1983, 17: 77

[本文引用: 1]

Kim Y S, Kim S J, Im K S.

Delayed hydride cracking in Zr-2.5Nb tube with the cooling rate and the notch tip shape

[J]. J. Nucl. Mater., 2004, 335: 387

[本文引用: 2]

Kim Y S, Grybenas A.

Effect of load ratio and hydrogen concentration on the crack growth rate in Zr-2.5Nb tubes

[J]. Mater. Sci. Eng., 2009, 520: 147

[本文引用: 1]

/