呂俊男，楊 爍，趙 毅，李 群，龍沖生,*

(1.中國(guó)核動(dòng)力研究設(shè)計(jì)院 反應(yīng)堆燃料及材料重點(diǎn)實(shí)驗(yàn)室，四川 成都 610213；2.西安交通大學(xué) 機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710049)

彌散燃料是將易裂變?nèi)剂舷?UO2、U3O8、U3Si2等)以小顆粒的形式彌散在相對(duì)惰性的非裂變材料基體(Al、不銹鋼、Mg等)中的一種混合物核燃料。彌散燃料芯體中的每個(gè)燃料顆粒都可單獨(dú)看作一個(gè)微型燃料元件。由于燃料顆粒被基體材料包圍和束縛，有利于包容裂變產(chǎn)物，抑制燃料腫脹，且導(dǎo)熱率高，從而允許達(dá)到較大塊燃料更高的燃耗，因此廣泛應(yīng)用于高安全可靠性、高比功率小堆芯的研究堆[1-2]。

輻照條件下，裂變氣體的快速生成會(huì)使燃料顆粒內(nèi)產(chǎn)生超壓和大量納米級(jí)的氣孔，致使燃料顆粒開(kāi)始腫脹，且隨著燃耗不斷加深，納米級(jí)氣孔聚合粗化形成微米級(jí)大氣孔，并逐漸演變成微米級(jí)氣孔占據(jù)主導(dǎo)的核燃料高燃耗結(jié)構(gòu)，且高燃耗結(jié)構(gòu)中氣孔的圓整度相對(duì)低燃耗情況更高[3-4]。在高溫和高燃耗條件下，燃料顆粒內(nèi)部氣孔中裂變氣體的壓力會(huì)達(dá)到上百兆帕[3]，事故工況下，燃料芯體內(nèi)部燃料相破損嚴(yán)重，大量裂變氣體釋放導(dǎo)致各顆粒之間的金屬基體發(fā)生開(kāi)裂連通，形成長(zhǎng)達(dá)數(shù)毫米的宏觀裂紋，最終導(dǎo)致燃料元件發(fā)生起泡現(xiàn)象?？梢哉f(shuō)燃料顆粒的開(kāi)裂失效是彌散燃料失效的關(guān)鍵因素[5-6]。

針對(duì)彌散燃料顆粒的開(kāi)裂行為研究，國(guó)內(nèi)外學(xué)者在這方面已開(kāi)展了一些有意義的工作[1-10]。但是，一方面前期研究是基于燃料顆粒內(nèi)氣孔尺寸、位置分布均勻性假設(shè)條件下開(kāi)展的[8-11]，與實(shí)際不符，忽視了氣孔尺寸、位置分布非均勻性對(duì)燃料內(nèi)部應(yīng)力場(chǎng)的隨機(jī)干涉效應(yīng)對(duì)其開(kāi)裂行為的影響；另一方面，前期工作未系統(tǒng)研究基體材料的力學(xué)性能、燃料顆粒間距、運(yùn)行環(huán)境靜水壓力及裂變碎片損傷層對(duì)燃料顆粒開(kāi)裂行為的影響。

本文將針對(duì)以上兩方面的難題提供解決方案，研究?jī)?nèi)容將有益于高效獲取高燃耗燃料顆粒的內(nèi)部應(yīng)力狀態(tài)，達(dá)到快速分析燃料顆粒強(qiáng)度力學(xué)性能和預(yù)測(cè)燃料開(kāi)裂溫度的目的，并為彌散型核燃料的輻照穩(wěn)定性評(píng)估和優(yōu)化設(shè)計(jì)提供技術(shù)支撐。

1 彌散燃料芯體的抽象與簡(jiǎn)化

1.1 彌散燃料芯體代表性體元

彌散燃料顆粒中含有大量的非均勻分布的高內(nèi)壓氣孔，氣孔平均尺寸和氣孔率隨燃耗加深而快速增長(zhǎng)，是導(dǎo)致燃料顆粒和基體開(kāi)裂行為的直接或間接驅(qū)動(dòng)力[10-12]。為系統(tǒng)研究基體材料對(duì)高燃耗燃料顆粒的約束效應(yīng)及強(qiáng)度影響規(guī)律，本文首先構(gòu)建了彌散燃料芯體的代表性體元，如圖1所示。

圖1 彌散燃料芯體代表性體元Fig.1 Representative volume element of dispersion fuel core

圖1中，彌散燃料芯體抽象單元由未損傷的基體層、裂變碎片損傷層和高燃耗燃料顆粒組成，相應(yīng)的特征尺寸半徑分別為Rm、Rl和Rf；彌散燃料芯體抽象單元外表面受到均布的環(huán)境壓應(yīng)力pm作用，需要注意的是，在數(shù)值計(jì)算分析燃料顆粒開(kāi)裂影響因素及規(guī)律時(shí)，pm的取值范圍為燃料元件運(yùn)行時(shí)所受的環(huán)境靜水壓力變化范圍；高燃耗燃料顆粒由致密的燃料相和彌散在其中的氣孔兩相組成，氣孔內(nèi)壓為pg。

1.2 高燃耗燃料顆粒的參數(shù)化建模

高燃耗燃料顆粒中，氣孔的尺寸和分布均會(huì)改變?nèi)剂项w粒內(nèi)部應(yīng)力場(chǎng)的分布，進(jìn)而影響開(kāi)裂行為。開(kāi)展考慮燃料顆粒內(nèi)部氣孔非均勻性的數(shù)值模擬研究，能更好地反映顆粒的強(qiáng)度特性。考慮這些參數(shù)必須參數(shù)化建模，通過(guò)參考輻照后UO2陶瓷燃料高燃耗結(jié)構(gòu)微觀形貌的電鏡照片(圖2)可知，高燃耗燃料內(nèi)彌散小氣孔的形狀近似為球形[3-6,13]，將氣孔形狀等效為球形也是相關(guān)研究的一貫處理方式[8-12]。

圖2 UO2陶瓷燃料高燃耗結(jié)構(gòu)的電鏡照片[4]Fig.2 SEM image of high burnup UO2 fuel structure[4]

為分析燃料顆粒內(nèi)部氣孔尺寸和位置分布的非均勻性對(duì)其開(kāi)裂行為的影響，本研究開(kāi)發(fā)了高燃耗燃料顆粒的參數(shù)化建模程序?；陔S機(jī)順序增加(RSA)的建模思想，結(jié)合Python編程工具，向燃料顆粒內(nèi)投放小氣孔，并記錄氣孔的坐標(biāo)和尺寸信息(xi，yi；ri)，通過(guò)循環(huán)判定算法確保任意氣孔相互不重疊，且氣孔半徑ri和相鄰氣孔間距di都服從Weibull分布。高燃耗燃料顆粒建模效果如圖3所示。

圖3 氣孔尺寸和位置非均勻分布的高燃耗燃料顆粒模擬效果圖Fig.3 Simulation of high burnup fuel particle with heterogeneous distribution of pore size and position

2 數(shù)值計(jì)算及結(jié)果分析

基于前文開(kāi)發(fā)的高燃耗燃料顆粒建模程序，進(jìn)一步使用Python腳本語(yǔ)言對(duì)Abaqus有限元二次開(kāi)發(fā)建立彌散燃料芯體的代表性體元幾何模型，并實(shí)現(xiàn)程序化加載、網(wǎng)格劃分和靜力學(xué)計(jì)算。根據(jù)實(shí)際燃耗、溫度及運(yùn)行環(huán)境條件，給每個(gè)裂變氣孔施加內(nèi)壓并在彌散燃料芯體的代表性體元外表面施加環(huán)境壓應(yīng)力，即可系統(tǒng)研究基體材料對(duì)高燃耗燃料顆粒的約束效應(yīng)及強(qiáng)度影響規(guī)律。

2.1 有限元計(jì)算基本輸入?yún)?shù)

如前文所述，高燃耗燃料顆粒中氣孔尺寸r和間距d均采用適應(yīng)性更強(qiáng)的Weibull概率密度分布函數(shù)模擬，其具體表達(dá)式為：

(1)

式中：x為隨機(jī)變量，本文中代表氣孔半徑r或氣孔間距d；λ為尺度參數(shù)；k為形狀參數(shù)。其中k反映變量的非均勻性，k越大則離散性越??；λ反映變量的平均特性。當(dāng)k取不同值時(shí)，可得到正偏、負(fù)偏和對(duì)稱的概率密度函數(shù)。值得注意的是，Spino等[4]統(tǒng)計(jì)的不同燃耗深度下UO2燃料內(nèi)部氣孔尺寸的分布實(shí)驗(yàn)數(shù)據(jù)已證實(shí)了Weibull概率密度分布模型的合理性。

為方便定性分析，本文有限元計(jì)算中采用的燃料顆粒材料為UO2陶瓷，燃料顆粒的半徑為160 μm，顆粒內(nèi)氣孔平均半徑為2 μm，且半徑r的Weibull概率密度分布函數(shù)雙參數(shù)設(shè)定為λ=2、k=5；氣孔平均間距為4 μm，且間距d的Weibull概率密度分布函數(shù)雙參數(shù)設(shè)定為λ=4、k=4。

此外，計(jì)算輸入的高燃耗燃料顆粒服役狀態(tài)參數(shù)列于表1，其中燃耗深度BU=0.3、環(huán)境溫度T=500 ℃。參考龍沖生等[10]的關(guān)于燃料顆粒開(kāi)裂數(shù)理模型研究，即可計(jì)算此服役狀態(tài)下燃料顆粒內(nèi)部氣孔率ρg和氣孔內(nèi)壓pg。

表1 燃料顆粒數(shù)值計(jì)算輸入?yún)?shù)Table 1 Calculation parameter of fuel particle model

2.2 參數(shù)分析

由于UO2燃料顆粒是脆性材料，可采用最大主應(yīng)力強(qiáng)度理論作為其開(kāi)裂判據(jù)，因此燃料顆粒內(nèi)部最大張應(yīng)力σmax即可反映其強(qiáng)度特性。在本文的有限元模擬研究中，將系統(tǒng)分析基體材料力學(xué)性能、燃料顆粒間距2(Rm-Rf)及環(huán)境壓應(yīng)力pm等因素對(duì)燃料顆粒最大張應(yīng)力σmax的影響規(guī)律。

1) 燃料顆粒間距及環(huán)境壓應(yīng)力對(duì)最大張應(yīng)力的影響規(guī)律

輸入表1所列基本參數(shù)，并輸入環(huán)境壓應(yīng)力pm=0.1 MPa、基體材料為304不銹鋼(彈性模量為146 GPa、泊松比為0.325)。由于UO2陶瓷熔點(diǎn)高達(dá)2 840 ℃，在高燃耗和高溫(500 ℃)服役條件下，燃料顆粒的塑性變形可忽略。因此本文將研究燃料顆粒間距2(Rm-Rf)和環(huán)境壓應(yīng)力pm變化對(duì)燃料內(nèi)部最大張應(yīng)力σmax的影響規(guī)律。

核裂變所產(chǎn)生的能量80%以上會(huì)轉(zhuǎn)變成裂變碎片的動(dòng)能，導(dǎo)致燃料顆粒近鄰的基體材料會(huì)受到高能裂變碎片的輻射，造成基體材料力學(xué)性能下降。其中，235U的裂變碎片在Al基體中的理論平均射程可達(dá)13.8 μm，在Fe基體中的理論平均射程為6.6 μm[6]。

本研究同時(shí)考慮了有無(wú)裂變碎片損傷層兩種情況下，燃料顆粒間距2(Rm-Rf)、環(huán)境壓應(yīng)力pm變化對(duì)燃料內(nèi)部最大張應(yīng)力σmax的影響規(guī)律，結(jié)果示于圖4，相應(yīng)有限元計(jì)算的應(yīng)力云圖如圖5a、b所示。

由圖4a可發(fā)現(xiàn)，高燃耗燃料顆粒內(nèi)部的最大張應(yīng)力σmax隨燃料顆粒間距的增大而減小，且減小的速率越來(lái)越慢。不考慮裂變碎片損傷層且環(huán)境壓應(yīng)力pm=0.1 MPa情況下，顆粒間距為20 μm時(shí)，σmax為192.4 MPa，顆粒間距為320 μm時(shí)，σmax為157.3 MPa，下降幅度達(dá)18.24%?？紤]裂變碎片損傷層影響后，即使損傷層厚度取20 μm的理論上限值(圖5b)，對(duì)σmax的影響也較小。具體地，顆粒間距為20 μm時(shí)，σmax小幅增大了3.69%；顆粒間距為320 μm時(shí)，σmax小幅增大了2.42%。

圖4 燃料顆粒間距及運(yùn)行環(huán)境壓應(yīng)力對(duì)最大張應(yīng)力的影響Fig.4 Influence of separation distance between fuel particles and ambient hydro-static pressure on maximum tensile stress

a——不含裂變碎片損傷層；b——含裂變碎片損傷層；c——輻照后UO2燃料顆粒斷裂行為的顯微分析[5]

上述研究結(jié)果說(shuō)明，燃料顆粒周圍起約束效果的基體材料層越厚，燃料顆粒越不易開(kāi)裂；裂變碎片損傷層的存在會(huì)小幅削弱基體材料對(duì)燃料顆粒的約束效應(yīng)，這種削弱效應(yīng)主要體現(xiàn)在損傷層彈性模量的下降，但這種效應(yīng)不明顯，不是影響燃料顆粒開(kāi)裂行為的主要因素?？紤]到燃料顆粒間距是彌散燃料顆粒微觀分布均勻性的重要特征參數(shù)，而均勻性關(guān)系到彌散燃料元件堆內(nèi)運(yùn)行時(shí)的溫度分布、輻照腫脹及起泡現(xiàn)象[14]。因此，本研究還間接說(shuō)明，彌散燃料芯體中局部燃料相體積分?jǐn)?shù)越大，燃料顆粒就越易開(kāi)裂，局部芯體出現(xiàn)宏觀裂紋的概率也越大。

此外，由圖4b可發(fā)現(xiàn)，高燃耗燃料顆粒內(nèi)σmax隨pm的增大而線性減小。具體地，當(dāng)燃料顆粒間距為160 μm時(shí)，pm從0.1 MPa增大到15 MPa，相應(yīng)的σmax從167.1 MPa減小至126 MPa，下降幅度為24.6%。說(shuō)明pm會(huì)抑制燃料顆粒的開(kāi)裂行為，且這種抑制作用較顯著。

當(dāng)燃料顆粒間距為20 μm時(shí)，pm從0.1 MPa增大到15 MPa，σmax下降幅度為23.2%(圖4a)；而燃料顆粒間距為320 μm時(shí)，相應(yīng)的σmax下降幅度為25.2%。本研究結(jié)果表明，燃料顆粒間距較小會(huì)削弱pm的影響程度，針對(duì)彌散燃料芯體中存在的更加明顯的燃料顆粒局部團(tuán)聚現(xiàn)象，pm對(duì)燃料顆粒開(kāi)裂行為的影響會(huì)更小。

從圖5a、b可直觀看出，由于高燃耗燃料顆粒內(nèi)氣孔尺寸和位置非均勻分布，顆粒內(nèi)應(yīng)力場(chǎng)存在復(fù)雜的疊加/干涉效應(yīng)，致使顆粒的開(kāi)裂危險(xiǎn)區(qū)從多處出現(xiàn)，且大多出現(xiàn)在靠近燃料顆粒的表層，本文數(shù)值計(jì)算結(jié)果與伍曉勇等[5]的實(shí)驗(yàn)觀察現(xiàn)象(圖5c)相吻合。

關(guān)于燃料顆粒表層開(kāi)裂現(xiàn)象，本研究認(rèn)為，燃料顆粒表層氣孔僅受到基體材料的約束，而靠近燃料顆粒中心區(qū)域的氣孔，受到基體材料和燃料相本身的共同約束，使得顆粒外層氣孔的平均張應(yīng)力水平較內(nèi)層氣孔更高，燃料顆粒在外層開(kāi)裂的概率更大。此外，對(duì)高燃耗顆粒內(nèi)開(kāi)裂危險(xiǎn)區(qū)觀察還發(fā)現(xiàn)，開(kāi)裂危險(xiǎn)區(qū)的氣孔尺寸普遍較大，且通常都是數(shù)個(gè)較大氣孔串聯(lián)在一起形成開(kāi)裂源。本研究佐證了Romano等[13]電鏡觀測(cè)到的高燃耗UO2燃料微觀結(jié)構(gòu)從氣孔較大處起裂的實(shí)驗(yàn)現(xiàn)象。

2) 燃料顆粒表面約束壓應(yīng)力pf數(shù)值計(jì)算

彌散燃料元件在堆內(nèi)運(yùn)行過(guò)程中，包殼外表面直接受到冷卻劑靜水壓力作用，靜水壓力ph可達(dá)到15 MPa。ph通過(guò)包殼和金屬基體的傳遞作用到燃料顆粒表面。但這種對(duì)燃料顆粒開(kāi)裂行為起到抑制作用的約束壓應(yīng)力pf的數(shù)值大小未知，亟待開(kāi)展數(shù)值計(jì)算研究。本研究構(gòu)建了如圖6所示的彌散燃料元件代表性體元，以系統(tǒng)分析外部靜水壓力ph對(duì)燃料顆粒表面約束壓應(yīng)力pf的影響規(guī)律。

圖6 燃料元件代表性體元有限元計(jì)算模型邊界條件Fig.6 Representative volume element of dispersion fuel and its boundary condition

對(duì)靠近燃料元件包殼的燃料顆粒受到的約束壓應(yīng)力pf進(jìn)行計(jì)算，相應(yīng)的燃料顆粒徑向應(yīng)力云圖如圖7所示。

取外部靜水壓力ph分別為0.1、5、10、15 MPa，提取基體與燃料顆粒表面一周節(jié)點(diǎn)上的徑向應(yīng)力。靜水壓力ph對(duì)pf數(shù)值分布的影響規(guī)律如圖8所示。由圖8a可見(jiàn)：燃料顆粒表面受到的金屬基體的約束壓應(yīng)力pf有波動(dòng)性，這是由于燃料顆粒內(nèi)部含有尺寸和位置非均勻分布的裂變氣孔所致，這一發(fā)現(xiàn)有益于對(duì)燃料顆粒開(kāi)裂溫度和開(kāi)裂位置的精準(zhǔn)預(yù)測(cè)；隨著ph的增大，燃料顆粒表面一周受到的pf波動(dòng)性降低。具體地，ph分別為0.1、5、10、15 MPa時(shí)，燃料顆粒表面一周pf數(shù)值波動(dòng)的標(biāo)準(zhǔn)差s分別為8.15、6.25、4.83、4.65。由圖8b可見(jiàn)，燃料顆粒表面一周受到的pf均值、最大值及最小值隨ph的增大而近乎線性增大，ph對(duì)pf最大值的影響程度明顯低于最小值。具體地，ph從0.1 MPa增大到15 MPa過(guò)程中，pf最大值增大了14.41 MPa，pf最小值增大了25.35 MPa。

圖8 靜水壓力ph對(duì)pf數(shù)值分布的影響Fig.8 Impact of hydro-static pressure ph on distribution of pf

3) 基體材料力學(xué)性能對(duì)最大張應(yīng)力的影響規(guī)律

輸入表1所列基本參數(shù)，并輸入環(huán)境壓應(yīng)力pm=0.1 MPa、燃料顆粒間距2(Rm-Rf)=160 μm，研究基體材料力學(xué)性能(彈性模量Em為變量)變化對(duì)燃料顆粒內(nèi)部最大張應(yīng)力σmax的影響規(guī)律，結(jié)果如圖9所示。由圖9可發(fā)現(xiàn)，高燃耗燃料顆粒內(nèi)部的最大張應(yīng)力σmax隨基體材料彈性模量的增大而線性降低。具體地，當(dāng)基體彈性模量Em從40 GPa增大到150 GPa時(shí)，相應(yīng)的σmax從188.2 MPa減小至165.1 MPa，下降幅度為12.3%。由圖9可直觀看出，若不考慮中子經(jīng)濟(jì)性，以304不銹鋼作為基體材料較銅、鋁/鎂合金對(duì)燃料顆粒的防護(hù)效果更好。本研究說(shuō)明，增大彌散燃料基體材料的彈性模量可在一定程度上抑制燃料顆粒的開(kāi)裂行為。

圖9 基體材料彈性模量變化對(duì)最大張應(yīng)力的影響Fig.9 Influence of elasticity modulus of matrix on maximum tensile stress

3 結(jié)論

本文通過(guò)有限元技術(shù)結(jié)合Python語(yǔ)言二次開(kāi)發(fā)功能，使得數(shù)值分析高燃耗燃料顆粒的強(qiáng)度影響因素及規(guī)律變?yōu)楝F(xiàn)實(shí)，得到如下主要結(jié)論。

1) 彌散燃料顆粒間距越大，燃料顆粒越不易開(kāi)裂；裂變碎片損傷層的存在會(huì)小幅削弱基體材料對(duì)燃料顆粒的約束效應(yīng)，使燃料顆粒開(kāi)裂風(fēng)險(xiǎn)小幅增大。

2) 高燃耗燃料顆粒內(nèi)氣孔尺寸、位置分布的非均勻性，會(huì)致使燃料顆粒的開(kāi)裂危險(xiǎn)區(qū)從多處出現(xiàn)；開(kāi)裂危險(xiǎn)區(qū)普遍具有氣孔尺寸較大且大氣孔串聯(lián)的特征；由于基體材料和燃料相的耦合約束效應(yīng)，燃料顆粒從外層開(kāi)裂的概率更大。

3) 燃料顆粒表面受到的金屬基體的約束壓應(yīng)力pf有波動(dòng)性；隨著靜水壓力ph的增大，燃料顆粒表面一周pf均值、最大值及最小值隨ph的增大而近乎線性增大，且ph對(duì)pf最大值的影響程度明顯低于最小值。

4) 增大彌散燃料基體材料的彈性模量，可在一定程度上抑制燃料顆粒的開(kāi)裂行為；彌散燃料元件所受環(huán)境靜水壓力越大，燃料顆粒越不易開(kāi)裂，而燃料顆粒間距縮短會(huì)小幅削弱靜水壓力的影響程度。