鋯合金輻照蠕變和生長的宏介觀模型研究進展

馬慧彤，劉桂森，崔嚴光，趙冠楠，沈 耀,*

(1.上海交通大學 材料科學與工程學院，上海 200240；2.上海核工程研究設計院，上海 200233)

反應堆燃料組件是反應堆堆芯的核心部件，其尺寸穩定性直接關系到核反應堆運行的安全性和經濟性。在反應堆運行過程中，燃料組件所處的工作條件十分苛刻，它不僅受到強烈的中子輻照，還會受到高溫高速冷卻劑流的侵蝕以及熱和機械應力的作用[1]。鋯合金因具有低熱中子吸收截面、較好的力學性能及高溫抗氧化、抗腐蝕能力，已成為反應堆堆芯燃料組件的主要結構材料，如壓水堆用燃料包殼管(Zr-2、Zr-4、M5、Zirlo等)及重水堆用壓力管(Zr-2.5Nb)等[2]。由于堆內服役環境嚴苛，鋯合金構件在服役過程中不僅會被氧化、腐蝕，更重要的是會輻照變形，這將嚴重影響其使用可靠性。預測鋯合金的堆內變形行為對評估其安全使用壽命、提供事故下的安全裕量具有重要意義，對開發核燃料預測系統具有重要價值。

鋯合金的輻照變形主要由材料在高溫、輻照、應力作用下缺陷(包括輻照缺陷和晶體固有缺陷)之間的相互作用及演化引起[3-5]，包括熱蠕變、輻照蠕變及輻照生長。輻照生長指各向異性材料在無外加應力作用下，輻照時發生的體積守恒的變形，是由輻照產生的點缺陷及其團簇與金屬中各向異性阱(位錯、位錯環及晶界)的擇優相互作用引起的[6-8]。蠕變則是在外加應力低于屈服應力時發生的緩慢塑性變形行為。熱蠕變是指高溫下，沒有輻照作用時材料的蠕變行為。輻照蠕變則是在輻照作用下，材料內部產生的大量點缺陷促進位錯的運動，導致的熱蠕變基礎上的額外蠕變增量[9]。這3種變形的疊加導致鋯合金構件產生了宏觀尺寸和形狀的改變。

近年來，國外有很多綜述總結了鋯合金的輻照變形行為和影響因素，以及相關變形機制的研究進展，但很少有文獻直接立足于其輻照變形的預測模型，國內更是缺乏相關的研究報道[9-12]。本文擬分別從宏介觀尺度詳細綜述兩類鋯合金構件的輻照變形預測模型。針對包殼管的模型有宏觀經驗模型及介觀力學模型。經驗模型基于堆內蠕變和生長數據，描述蠕變和生長應變與環境參量(應力、溫度、中子注量)之間的經驗關系，是當前應用在壓水堆包殼管設計及安全許可分析中的主要模型，但其特定的經驗參數值只適用于具有特定材料特性并處于特定服役環境下的鋯合金。鑒于經驗模型的局限性，最新發展起來的介觀力學模型變得尤為重要。介觀力學模型基于鋯合金輻照變形的物理機制，在微觀尺度考慮輻照過程中微觀結構的演化，在介觀尺度通過粘塑性自洽(VPSC)程序考慮織構、晶粒尺寸及形狀對輻照變形的影響，模型具有廣泛適用性。針對壓力管，主要是宏觀經驗模型，包括C6方程及秦山方程。C6方程考慮了材料的微觀結構及織構對輻照變形各向異性的影響，具備一定的物理基礎，是當前應用在CANDU重水堆壓力管壽命預估中的主要模型，但模型關鍵各向異性參數及計算方法并未詳細公開；秦山方程基于秦山CANDU重水堆壓力管輻照變形的統計規律，是一種新的壓力管輻照變形預測方程。因此，通過對壓水堆用鋯合金包殼管及重水堆用Zr-2.5Nb壓力管輻照變形預測模型的詳細綜述，不僅對評估鋯合金的輻照蠕變和生長行為具有參考價值，對Zr-2.5Nb壓力管國產化研究也具有重要意義。

1 鋯合金包殼管的輻照變形預測模型

宏觀經驗模型及介觀力學模型都可用于模擬包殼管的蠕變和生長行為。宏觀經驗模型基于包殼管熱蠕變、輻照蠕變、輻照生長的經驗變形規律，通常通過有限元方法求解包殼管變形過程中的力學平衡問題。介觀力學模型基于鋯合金輻照變形的物理機制，主要通過VPSC程序耦合鋯合金熱蠕變、輻照蠕變、輻照生長的單晶模型，考慮晶體織構以及晶粒尺寸和形狀效應，預測鋯合金多晶的蠕變和生長行為。

1.1 包殼管的輻照蠕變和生長行為的基本規律

不同類型的鋯合金包殼管的蠕變應變隨中子注量變化的曲線示于圖1[2]。由于鋯合金包殼管織構(Kearns織構因子通常為fr=0.636,ft=0.317,fa=0.047[2]，r、t、a分別代表包殼管的3個宏觀特征方向(徑向、環向和軸向))，包殼管的環向應變通常只由蠕變貢獻，而沒有輻照生長的貢獻[13]。鋯合金的蠕變變形由熱蠕變和輻照蠕變組成，但當前測試技術很難將兩者分離開來。對不同類型的鋯合金，蠕變曲線都可分為初始的瞬態階段和穩態階段。初始階段由于輻照劑量較小、輻照硬化不明顯，通常有很高的蠕變率。隨著輻照損傷的累積，材料被充分硬化，熱蠕變被抑制，蠕變速率逐漸下降至穩定，進入穩態變形階段。

圖1 不同類型鋯合金包殼管的輻照蠕變曲線[2]Fig.1 Irradiation creep curves of different types of zirconium alloy cladding tubes[2]

圖2 不同類型鋯合金包殼管輻照生長曲線[2]Fig.2 Irradiation growth curves of different types of zirconium alloy cladding tubes[2]

不同類型的鋯合金包殼管軸向生長應變隨輻照劑量(中子注量) 的變化示于圖2，大致分為3個階段[2]：輻照初期的瞬態階段，由于〈a〉型位錯環大量成核，生長應變迅速增加；隨著輻照劑量的增加，〈a〉型位錯環逐漸飽和，生長率下降至穩定，此時進入輻照生長第2階段；當輻照劑量超過一定值時，〈c〉型位錯環開始成核，生長率明顯增加，進入“加速生長”階段。RXA(再結晶退火)鋯合金相比CW(冷加工)及SRA(消應力退火)鋯合金的“加速生長”趨勢更加明顯。

宏觀經驗模型及介觀力學模型都可描述包殼管蠕變和生長不同階段的變形規律。對不同類型的包殼管材料，其蠕變和生長應變隨中子注量變化的趨勢是相似的(如圖1～2)，因此可通過調整經驗模型及力學模型的模型參數預測不同類型鋯合金包殼管材料的蠕變和生長行為。

1.2 宏觀經驗預測模型

宏觀經驗模型假設輻照變形由獨立的、可疊加的3部分變形構成[3]，包括熱蠕變、輻照蠕變、輻照生長，總變形率為：

(1)

Matsuo[14]根據Zr-4包殼管在內壓作用下的熱蠕變特征，提出了適用于鋯合金包殼管的熱蠕變模型，該模型通過雙曲正弦函數將穩態熱蠕變率與外加應力聯系起來。根據應力水平，指數n的取值范圍在2～4之間[15-16]。

(2)

文獻[3,14,17]在式(2)的基礎上，同時考慮材料的初始熱蠕變及穩態熱蠕變，描述總的熱蠕變應變隨時間的變化：

(3)

Hoppe[18]假設輻照蠕變主要是由非熱機制引起的，提出了適用于鋯合金包殼管的輻照蠕變經驗模型，輻照蠕變速率與外加應力和中子注量率呈冪律關系：

(4)

Luscher等[16]在Hoppe等的基礎上，考慮了輻照蠕變對溫度的依賴，同時修正應力項為等效應力，并根據Franklin等[19]的試驗數據擬合修正了模型參數：

(5)

其中，f(T)為溫度依賴系數。對于SRA及RXA鋯合金包殼管，應力指數C2=1.0、快中子注量率指數C1=0.85、溫度依賴系數f(T)=0.72～1.47(服役溫度在300～350 ℃之間)。當等效應力形式基于屈服準則以及應變方向由屈服面法線決定時，這種修正形式類似于塑性流動的常規表達式。

Geelhood等[16,20]提出了鋯合金包殼管的輻照生長模型，在輻照生長的不同階段，輻照生長應變以與時間有關的增量形式表示：

εt+Δt=εt+Δε

(6)

Δε=AGd(φΔt)m

(7)

其中：εt為t時刻的生長應變；Δε為Δt時間內的生長應變增量；φ為快中子注量率；A、Gd及m為材料常數，A和m由合金成分和加工方式決定，Gd為生長各向異性因子，在d特征方向上，Gd=1-3fd，fd為Kearns各向異性織構因子。

1.3 介觀力學預測模型

Montgomery等[3]提出了以VPSC框架為基礎的介觀模型思想。即通過VPSC程序耦合鋯合金的單晶蠕變和生長模型，在多晶尺度上綜合考慮多晶織構、晶粒間相互作用、晶粒尺寸和晶粒形狀效應，模擬鋯合金多晶的輻照變形行為。

1) 輻照變形物理機制與單晶模型

熱蠕變機制主要有擴散控制及位錯控制兩種類型。擴散控制的蠕變機制假設熱蠕變率與應力呈線性依賴關系[21-22]，包括Nabarro-Herring蠕變(體擴散)[23]及Coble蠕變(晶界擴散)[24]。位錯攀移控制的蠕變機制假設熱蠕變率與應力滿足冪律關系，冪指數隨應力的變化在4～9之間分布[25-27]。Montgomery等[3]采用黏塑性冪律規則描述不同滑移系對熱蠕變的貢獻，在鋯合金構件正常服役的情況下，冪指數n取4，單晶熱蠕變率為：

j=1,2,…,5

(8)

(9)

(10)

輻照蠕變機制主要有應力誘導的優先吸收(SIPA)機制[29-30]、攀移增強滑移機制[31]、應力誘導的優先成核(SIPN)機制[32]及擴散傳質機制[33]，其中前兩者較為突出[11]。SIPA機制假設自間隙原子優先被柏氏矢量平行于應力方向的刃位錯捕獲，空位優先被柏氏矢量垂直于應力方向的刃位錯捕獲，這種捕獲偏差導致材料產生沿應力方向的伸長現象[34-35]。攀移增強的滑移機制假設刃位錯捕獲輻照產生的點缺陷發生攀移，克服障礙后在另一個滑移面上繼續滑移[36-38]，從而進一步產生變形，最具代表性的是Gittus的I-creep模型[39-40]。目前學術界對鋯合金的輻照蠕變機制仍無明確定論，包殼管的輻照蠕變計算仍以經驗模型為主。Patra等[41]在經驗模型基礎上，假設蠕變率與應力及中子劑量率呈正比，在單晶尺度上考慮每個滑移系對變形的貢獻，開發了輻照蠕變晶體學模型，單晶輻照蠕變率為：

j=1,2,…,5

(11)

輻照生長的物理模型主要有輻照損傷理論(RDT)模型、點缺陷各向異性擴散差異(DAD)模型及反應-擴散模型。RDT模型假設輻照損傷過程產生等量的空位和自間隙原子，它們的運動都是各相同性的，微觀結構演化的驅動力來自刃位錯對自間隙原子的優先捕獲[42-43]，這并不能解釋試驗觀察到的輻照生長現象[44-45]。在此基礎上，DAD模型考慮到自間隙原子優先沿基面a方向遷移，易被柱面位錯捕獲；空位遷移近似為各向同性，易被基面位錯捕獲，這可解釋鋯合金晶胞沿c軸的收縮行為，以及〈c〉位錯環在鋯合金“加速生長”的行為中扮演至關重要的角色[46-47]。但DAD模型沒有考慮到輻照使材料內部產生大量可動缺陷團簇的情況，這與試驗和分子動力學模擬結果不符[34,48-49]。Golubov等[44-45]將模擬bcc和fcc晶體輻照腫脹的生產偏差模型(PBM)思想應用于鋯合金輻照生長模型中，稱為反應-擴散模型。該模型考慮最初的輻照損傷由點缺陷及自間隙原子團簇組成，忽略了點缺陷遷移的各向異性，但考慮自間隙原子團簇沿基面一維遷移。Patra等[41]在Golubov等的基礎上，進一步考慮了晶界捕獲輻照點缺陷發生收縮或膨脹對輻照生長的貢獻，具體如圖3所示。

圖3 鋯合金HCP晶胞的晶體學特征[2] Fig.3 Crystallographic characteristics of zirconium alloy HCP elementary cell[2]

Patra等考慮輻照生長率由位錯捕獲點缺陷發生攀移及晶界吸收點缺陷兩部分貢獻。位錯捕獲點缺陷對輻照生長的貢獻(忽略位錯線與位錯環捕獲強度的差異)為：

(12)

j=a1,a2,a3

(13)

(14)

晶界吸收點缺陷發生收縮或膨脹由晶界捕獲點缺陷的凈通量決定，對輻照生長的貢獻為：

j=a1,a2,a3,c；m=x′,y′,z′

(15)

由以上公式可得總輻照生長率：

(16)

Christiaen等[50]采用蒙特卡羅方法計算分析得出，鋯合金中空位運動也是各向異性的，沿平行于基面的方向擴散較快。Patra等[41]的反應-擴散模型假設空位和自間隙原子運動都是各相同性的，雖然模擬輻照生長應變隨輻照劑量的變化趨勢與試驗結果較為接近，但其在理論上仍存在進一步發展空間。

2) VPSC框架

VPSC程序由Tomé等于1993年提出，基于單晶變形時滑移和孿生的物理剪切機制，考慮材料織構及晶粒間相互作用，模擬多晶在外部約束(應力、應變、高溫、輻照等)下的各向異性的塑性變形[51]。VPSC將每個晶粒看作一個嵌入在有效黏塑性介質中的橢球形黏塑性夾雜，夾雜和介質都是各向異性的(圖4)。當對有效介質施加應力、應變率或輻照時，介質與晶粒發生相互作用，在晶粒及其附近產生局部應力。當晶粒為橢球時，晶粒中的應力是均勻的，但可能與介質中應力存在偏差。VPSC的自洽性即通過循環迭代強制實現所有晶粒上的平均應力和應變率等于施加在介質上的宏觀應力或應變率[51-52]。

圖4 VPSC機制示意圖Fig.4 Schematic diagram of VPSC mechanism

式(17)基于非線性速率敏感性方程，通過對所有滑移系s上的變形速率累加，描述了VPSC中的單晶變形速率(VPSC在計算過程中應力及應變率均用5維矢量形式表示[51])，并假設變形速率與應力呈線性關系：

i=1,2，…,5

(17)

(18)

(19)

其中：I為單位矩陣；S為Eshelby張量，與晶粒形狀及宏觀粘塑性蠕變柔度張量有關。

(20)

(21)

3) 輻照蠕變和生長在VPSC中的實現

Patra等以VPSC框架為基礎，在VPSC中實現了鋯合金輻照生長的反應-擴散模型與輻照蠕變的晶體學模型的耦合[41]，并根據冷加工Zr-2合金的微觀結構演化及輻照變形數據確定了模型參數，初步預測了單軸拉伸應力作用下鋯合金的軸向變形行為。其預測結果較好地反映了試驗規律[41,53-54]，表明該介觀模型整體上是可靠的。但Patra等忽略了輻照初期熱蠕變的貢獻，而堆內試驗結果表明，在輻照初期的瞬態階段，熱蠕變的貢獻不可忽略(圖1)。本文在Patra 等模型的基礎上，增加了熱蠕變修正項，假設單晶變形速率為：

i=1，2，…，5

(22)

(23)

當前的介觀力學模型展現了較好的預測能力，但目前輻照蠕變機制尚無明確定論，單晶蠕變模型采用了以經驗模型為基礎的晶體學模型，而輻照生長的反應-擴散模型中也忽略了點缺陷各向異性運動這一重要物理現象。因此，介觀力學模型仍存在進一步的發展空間，其發展核心應是鋯合金蠕變和生長物理機制的持續探索以及單晶模型的開發。

2 Zr-2.5Nb壓力管的輻照變形預測模型

目前，Zr-2.5Nb壓力管的輻照變形預測模型主要為宏觀經驗模型，雖然1.3節介紹的介觀力學模型框架理論上也適用于壓力管，但具體應用實例并未見文獻報道。宏觀經驗模型包括C6方程及秦山方程。C6方程是Christodoulou等[55]于1996年提出的適用于CANDU重水堆壓力管的模型，可綜合考慮壓力管的服役條件、織構、各向異性微觀結構特征(位錯、晶粒形狀)，預測壓力管各向異性的輻照變形，但相關各向異性參數并沒有公開。秦山方程是唐迥然等[56]通過對秦山壓力管的在役檢查數據進行統計分析，在統計規律的基礎上建立的一種新的壓力管輻照變形計算方程，該方程符合秦山重水堆壓力管變形規律，但針對不同的重水堆，方程參數需要重新擬合。

圖5 SRA Zr-4包殼管輻照變形預測結果與試驗對比Fig.5 Comparison between predicted result and experimental result of SRA Zr-4 alloy cladding tube

圖6 RXA Zr-4包殼管輻照變形預測結果與試驗對比Fig.6 Comparison between predicted result and experimental result of RXA Zr-4 cladding tube

在重水堆國產化進程中，壓力管材料也需國產化，應選擇一種合適的方程預測及評估國產壓力管的壽命。秦山方程中的參數不能直接應用于國產化的壓力管材料，當前也缺少堆內變形數據支持擬合，而C6方程中有相關文獻公開的部分參數，僅涉及到微觀結構及織構的各向異性參數沒有明確給出(各向異性參數可基于國產化壓力管的織構及微觀結構計算)。鑒于C6方程有在CANDU重水堆(包括秦山重水堆)中的實際應用基礎[57]，C6方程預測國產化壓力管的輻照變形行為應是較為合適且可行的選擇，也需相應地國產化。

2.1 壓力管的輻照蠕變和生長行為的基本規律

由于壓力管具有與包殼管明顯不同的織構(Kearns織構因子通常為fr=0.33、ft=0.61、fa=0.06)，其蠕變和生長應變也體現出和包殼管明顯不同的各向異性。OSIRIS測試堆(中子注量率為2×1018m-2·s-1)內微壓力管(MPT，織構和微觀結構與標準壓力管類似，尺寸較小)的蠕變應變隨中子注量變化的曲線示于圖7a[58]，雙軸應力(環向應力為160 MPa)作用下，軸向和環向的蠕變應變均為正值(包殼管的軸向蠕變應變為負值)，且都沒有體現出明顯的瞬態階段，可能是因為中子注量率較高，輻照初期熱蠕變即被明顯抑制。微壓力管的輻照生長應變隨中子注量的變化曲線示于圖7b[59]，與包殼管不同的是，壓力管的環向生長應變不為0，軸向生長應變約為環向生長應變的2倍，同時也沒有明顯的“加速生長”階段。

圖7 OSIRIS堆內Zr-2.5Nb微壓力管的輻照蠕變和生長曲線Fig.7 Irradiation creep and growth curve of Zr-2.5 Nb micropressure tube in OSIRIS reactor

微壓力管的蠕變和生長應變隨中子注量的變化規律能直接反映標準壓力管蠕變和生長行為的宏觀規律，但與微壓力管不同，標準壓力管長6 m、直徑104 mm、厚度4.2 mm，沿軸向不同位置，其微觀結構與織構略有不同，服役時溫度、應力及中子注量分布也有一定差異[12]。這些因素導致壓力管產生沿軸向的不均勻變形(圖8)，進一步增加了預測模型開發的難度。

圖8 典型的CANDU 壓力管平均內徑沿軸向的分布[12]Fig.8 Average inner diameter distribution along axis of typical CANDU pressure tube[12]

2.2 宏觀經驗預測模型

1) CANDU 壓力管C6方程

加拿大原子能公司(AECL)對Zr-2.5Nb壓力管進行了多年材料輻照試驗研究(包括輻照生長、輻照蠕變和熱蠕變試驗)，建立了Zr-2.5Nb壓力管在CANDU反應堆中的輻照變形計算方程(C6方程)，并不斷進行改進，于1996年最終定型。總變形率為3種變形模式的疊加，其中熱蠕變率、輻照蠕變率、輻照生長率[55]分別為：

(24)

(25)

(26)

σi=[Fi(σa-σt)2+Gi(σt-σr)2+

Hi(σr-σa)2]1/2i=1,2

(27)

其中，F、G和H為Hill各向異性屈服常數。對于熱蠕變，i=1或2 (應力小于120 MPa，i=1起主導作用；應力大于200 MPa，i=2起主導作用)。由于熱蠕變對總變形的貢獻較小，忽略Fi、Gi、Hi對位置x的依賴，而輻照蠕變考慮了壓力管不同位置處織構變化對變形各向異性的影響，各向異性屈服常數隨位置x的變化為：

F(x)=Fb+(Ff-Fb)x/6

(28)

G(x)=Gb+(Gf-Gb)x/6

(29)

H(x)=1.5-F(x)-G(x)

(30)

其中：Ff、Gf、Hf分別為壓力管前端的輻照蠕變各向異性屈服常數；Fb、Gb、Hb分別為壓力管后端的輻照蠕變各向異性屈服常數。

熱蠕變及輻照蠕變的各向異性參數可表示為：

(31)

(32)

(33)

輻照生長各向異性參數為：

(34)

(35)

(36)

Christodoulou等[55]根據壓力管堆內輻照變形數據及壓力管的織構和微觀結構特征擬合并計算了上述方程的參數，預測結果與堆內變形數據的對比示于圖9，表明該模型能較好地預測壓力管軸向及環向的變形行為。

2) 秦山壓力管輻照變形計算方程

唐迥然等[56]對大量秦山CANDU重水堆壓力管的在役檢查數據進行統計分析，得出輻照蠕變速率及輻照生長速率分別與快中子注量率呈冪律關系，且熱蠕變是不可忽略、非穩態的行為，這與C6方程明顯不同。在此基礎上，他們進一步建立了秦山重水堆壓力管熱蠕變率、輻照蠕變率及輻照生長率的計算方程，總變形率為3種變形模式的疊加：

(37)

(38)

(39)

秦山方程中沿壓力管不同方向的輻照蠕變率及輻照生長率的參數需分別擬合，大量的參數擬合需要大量的堆內運行數據支持。在缺少堆內運行數據支持的情況下，C6方程預測國產化壓力管的輻照變形行為應是較為合適的選擇。而秦山方程描述的壓力管的變形規律與C6方程明顯不同，這也應作為下一步對C6方程進行改進的參考。

3) C6方程的國產化

CANDU C6方程考慮了微觀結構及織構對輻照變形各向異性的影響，具備一定的物理基礎，這主要體現在C6方程的各向異性參數中，而在1996年的公開文獻中，這些參數并沒有被公開，其計算方法也沒有被詳細闡述。本文將壓力管熱蠕變、輻照蠕變、輻照生長的單晶模型和VPSC程序耦合，輸入國產化壓力管的初步樣品織構，計算了壓力管熱蠕變及輻照蠕變的各向異性屈服常數F、G、H，以及輻照生長各向異性因子，并將計算結果與Causey等計算的OSIRIS微壓力管[59]進行了對比。

圖9 Zr-2.5Nb壓力管在CANDU反應堆中實測和預測結果對比[55]Fig.9 Comparison of measured and predicted results of Zr-2.5Nb pressure tube in CANDU reactor[55]

單晶熱蠕變模型基于黏塑性冪律準則，冪指數n取2：

j=1,2,…,5

(40)

單晶輻照蠕變模型基于I-creep機制[61]：

k1,k2=k5,k3=k4,k6=0

(41)

其中：ki為第i類變形模式的蠕變柔度，k1為錐面滑移貢獻，k2為柱面和錐面滑移貢獻，k3為基面和錐面滑移貢獻；σi為相應變形模式下的晶粒應力。k1、k2、k3由材料微觀結構決定。

將上述單晶模型與VPSC程序耦合，采用Causey等[59]的單晶模型參數，輸入國產壓力管初步樣品織構(織構因子fr=0.34、ft=0.57、fa=0.09)，分別計算熱蠕變各向異性屈服常數Fth、Gth、Hth，輻照蠕變各向異性屈服常數F、G、H，以及輻照生長各向異性因子Ga/Gt，計算結果與Causey等[59]計算的OSIRIS微壓力管(織構因子fr=0.33、ft=0.61、fa=0.06)結果的對比列于表1。

表1 國產壓力管樣品與OSIRIS微壓力管的計算結果對比Table 1 Comparison of calculated results between localized pressure tube and OSIRIS micro-pressure tube

在織構相近的情況下，計算所得Hill各向異性屈服常數與文獻結果較接近，且軸向與環向的生長率比值也較接近，可進一步基于式(28)～(36)計算各向異性參數。表明可通過該方法定量評估織構及微觀結構對國產壓力管輻照變形各向異性的影響。如果國產化壓力管織構和微觀結構與加拿大原產品接近，則這些關鍵的各向異性參數也應接近，且其差異可定量評估。

C6方程中壓力管的輻照蠕變率及輻照生長率均與快中子注量率呈線性關系，這與秦山方程的乘冪形式明顯不同，下一步可考慮參考秦山方程的形式，結合已發表的CANDU壓力管蠕變和生長的堆內試驗數據對C6方程進行進一步改進。但無論是C6方程還是秦山方程，都是經驗性的，其經驗參數值受材料類型及環境參量影響，具有一定的局限性。而理論上1.3節介紹的介觀力學模型框架也適用于壓力管，鑒于介觀力學模型具有更廣的適用性，應是更進一步的發展方向。

3 總結與展望

鋯合金構件在堆內服役過程中發生的蠕變及生長行為嚴重影響了其使用可靠性。預測鋯合金的輻照蠕變和生長對保障反應堆安全、有效運行具有重要意義。

對壓水堆用鋯合金包殼管，宏觀經驗模型及介觀力學模型都可預測其蠕變和生長行為。介觀力學模型基于鋯合金輻照變形物理機制，理論上可預測具有任意微觀結構特征的鋯合金多晶在任意環境參量作用下的蠕變及生長行為，相比宏觀經驗模型，其理論性更強，適用范圍更廣，將是未來的發展趨勢。當前介觀模型已初步展現了較好的預測能力，但其蠕變和生長物理機制及相應的單晶模型仍存在進一步的發展空間。

對重水堆用Zr-2.5Nb壓力管，主要的預測模型為宏觀經驗模型，包括CANDU C6方程及秦山方程。C6方程相比秦山方程具備一定的物理基礎，是預測國產化壓力管材料輻照變形行為較為合適的選擇。當前C6方程的國產化研究已取得了階段性進展，但用C6方程描述的壓力管的蠕變和生長速率對中子注量率的依賴規律與秦山方程有明顯區別，下一步應考慮參考秦山方程的形式，結合壓力管蠕變和生長的堆內試驗數據對C6方程進行進一步修正。同時，考慮到C6方程和秦山方程均為經驗模型，其適用范圍受限，因此，結合介觀力學方法發展機理性的輻照蠕變和生長模型應是未來的一個重要發展方向。