(中國(guó)核動(dòng)力研究設(shè)計(jì)院，四川 成都 610213)

燃料組件壓緊系統(tǒng)驗(yàn)證中，需分析確定在正常運(yùn)行期間壓緊組件能防止燃料組件跳起，以保證燃料組件的完整性。

目前，在核電工程上通常采用的壓緊力驗(yàn)證分析方法為確定論分析方法，并采用最保守的假設(shè)，所有的設(shè)計(jì)參數(shù)均取對(duì)最小壓緊力最不利的極限值。造成計(jì)算得到的壓緊力與堆內(nèi)的實(shí)際值差別很大。這種分析方法大大增加了燃料組件防跳起的安全性，但同時(shí)可能出現(xiàn)壓緊力裕量不足的情況。如果就此修改壓緊系統(tǒng)設(shè)計(jì)來(lái)增加壓緊力，使得驗(yàn)證結(jié)果滿(mǎn)足要求；此時(shí)真實(shí)的壓緊力可能會(huì)過(guò)大，反而會(huì)造成燃料組件發(fā)生彎曲，從而可能延長(zhǎng)控制棒的落棒時(shí)間，甚至出現(xiàn)不完全落棒，進(jìn)而影響反應(yīng)堆的安全運(yùn)行[1]。因此，為優(yōu)化壓緊系統(tǒng)的設(shè)計(jì)，在保證足夠壓緊力的同時(shí)減小對(duì)燃料組件彎曲的影響，需研究更為合理的壓緊力計(jì)算方法。

在保證壓緊的前提下減小施加到燃料組件上的軸向壓緊力成為壓緊系統(tǒng)設(shè)計(jì)新的關(guān)注重點(diǎn)，也對(duì)壓緊系統(tǒng)的壓緊力大小提出了更加苛刻的要求。壓緊系統(tǒng)設(shè)計(jì)時(shí)不能為了保證燃料組件的壓緊而僅從結(jié)構(gòu)設(shè)計(jì)上增加壓緊力，也需要通過(guò)更為精確的壓緊力計(jì)算方法來(lái)增加壓緊力值。為了在保證壓緊系統(tǒng)壓緊功能的同時(shí)，燃料組件不致承受過(guò)大的軸向壓緊力，本文研究采用蒙特卡洛方法計(jì)算了燃料組件板彈簧壓緊力對(duì)燃料組件壓緊影響，并對(duì)蒙特卡洛方法的精準(zhǔn)度進(jìn)行了研究。

1 壓緊系統(tǒng)分析驗(yàn)證方法

1.1 方法概要

燃料組件的壓緊狀態(tài)與反應(yīng)堆及燃料組件自身的眾多參數(shù)相關(guān)，無(wú)論采用哪一種方法進(jìn)行壓緊力的計(jì)算和驗(yàn)證，都需要考慮這些因素影響，將相關(guān)參數(shù)作為計(jì)算輸入。根據(jù)不同參數(shù)的特性，可將輸入?yún)?shù)劃分為可變參數(shù)和固定參數(shù)兩大類(lèi)。

對(duì)于給定的燃料組件、壓緊系統(tǒng)和堆芯，固定參數(shù)的取值不會(huì)隨計(jì)算方法發(fā)生變化。其取值的確定也會(huì)考慮相應(yīng)的保守性，如壓緊板彈簧剛度特性采用根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合得到的下限值(最小剛度)。

對(duì)于可變參數(shù)，數(shù)據(jù)組成包括名義值及其不確定度，不同的計(jì)算方法差異在于對(duì)于參數(shù)不確定的處理，對(duì)于整個(gè)計(jì)算所考慮的物理模型相同?？勺儏?shù)包括堆腔高度、組件軸向結(jié)構(gòu)參數(shù)、壓緊系統(tǒng)結(jié)構(gòu)參數(shù)及輻照特性等。

1.2 確定論方法

目前工程設(shè)計(jì)中一般采用保守假設(shè)的確定論方法(Deterministic Methodology)進(jìn)行壓緊系統(tǒng)驗(yàn)證。通過(guò)對(duì)可變參數(shù)的名義值加上或者減去不確定度以使平衡力達(dá)到最小，從而使計(jì)算結(jié)果達(dá)到最保守的程度，以包絡(luò)實(shí)際反應(yīng)堆中所有燃料組件的最危險(xiǎn)壓緊狀態(tài)[2]。

采用確定論方法計(jì)算壓緊力時(shí)，對(duì)各種影響平衡力的參數(shù)均取了極限值。但實(shí)際運(yùn)行中，各種影響平衡力的參數(shù)均同時(shí)處于極限值的概率是非常小的。因此，采用確定論方法驗(yàn)證計(jì)算得到的最小平衡力過(guò)于保守[3]。

1.3 蒙特卡洛方法

蒙特卡洛方法(Monte Carlo Algorithm)也稱(chēng)統(tǒng)計(jì)模擬方法，是一種以概率統(tǒng)計(jì)理論為指導(dǎo)的一類(lèi)非常重要的數(shù)值計(jì)算方法，可以使用隨機(jī)數(shù)(或常見(jiàn)的偽隨機(jī)數(shù))來(lái)解決很多計(jì)算問(wèn)題。

采用蒙特卡洛方法計(jì)算時(shí)，使用數(shù)學(xué)模擬的方法，進(jìn)行大量模擬實(shí)際情況的計(jì)算。每次計(jì)算時(shí)，可變參數(shù)按照一定的概率分布規(guī)律隨機(jī)生成數(shù)據(jù)，并經(jīng)模型計(jì)算得到相應(yīng)的輸出結(jié)果。最終對(duì)大量的計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì)處理，得到在一定置信度下壓緊力和平衡力的統(tǒng)計(jì)輸出結(jié)果[4]。為了保證蒙卡方法計(jì)算壓緊力的可靠性及適用性，需選擇適宜的輸入?yún)?shù)分布類(lèi)型及數(shù)據(jù)統(tǒng)計(jì)時(shí)的置信度。

2 壓緊驗(yàn)證結(jié)果對(duì)比

采用蒙特卡洛方法驗(yàn)證燃料組件的壓緊，需要開(kāi)發(fā)配套的軟件，目前中國(guó)核動(dòng)力研究設(shè)計(jì)院自主研發(fā)的HOFA軟件可以實(shí)現(xiàn)該功能[5]。經(jīng)分別采用確定論方法及蒙特卡洛方法對(duì)燃料組件進(jìn)行壓緊力驗(yàn)證，壓緊力及平衡力計(jì)算結(jié)果對(duì)比見(jiàn)表1，其中蒙特卡洛方法中的置信度為95%。

通過(guò)對(duì)采用確定論方法及蒙特卡洛方法計(jì)算的最小壓緊力、平衡力數(shù)據(jù)進(jìn)行對(duì)比，可以看出在未考慮最小設(shè)計(jì)裕量時(shí)，采用確定論方法計(jì)算的平衡力已經(jīng)不足，不滿(mǎn)足設(shè)計(jì)準(zhǔn)則。采用蒙特卡洛方法計(jì)算的壓緊力及平衡力都能獲得較大的提升。相較于確定論方法，采用蒙特卡洛方法最少獲取755 N的提升，在壓緊情況最?lèi)毫拥臒釕B(tài)邊緣區(qū)仍有707 N的裕量，滿(mǎn)足了設(shè)計(jì)準(zhǔn)則要求的最小250 N的裕量，從而避免由于計(jì)算方法的過(guò)度保守性導(dǎo)致壓緊力不足的問(wèn)題。

表1 最小力計(jì)算結(jié)果對(duì)比

3 壓緊力離散度對(duì)比分析

早期的燃料組件設(shè)計(jì)中，板狀壓緊彈簧設(shè)計(jì)中重點(diǎn)關(guān)注了壓緊系統(tǒng)的壓緊功能以保證燃料組件運(yùn)行過(guò)程中不發(fā)生跳起，未考慮過(guò)大壓緊力對(duì)燃料組件彎曲的影響。

采用確定論方法計(jì)算的壓緊力為最極限的值，離散度大，偏離燃料組件堆內(nèi)所受的實(shí)際壓緊力較大。以前在壓緊系統(tǒng)設(shè)計(jì)時(shí)，如果最小壓緊力不能滿(mǎn)足壓緊的要求，設(shè)計(jì)者往往會(huì)修改結(jié)構(gòu)設(shè)計(jì)來(lái)增加壓緊力。如果分析方法得到的結(jié)果不夠精確，燃料組件在堆內(nèi)承受的實(shí)際壓緊力會(huì)比計(jì)算得到的最小壓緊力大得多，可能會(huì)對(duì)燃料組件產(chǎn)生不利的影響。壓緊系統(tǒng)設(shè)計(jì)不能繼續(xù)按照保守壓緊力的方式進(jìn)行設(shè)計(jì)，需要更為準(zhǔn)確的壓緊力分析方法來(lái)計(jì)算壓緊力值，以避免過(guò)大壓緊力對(duì)組件的不利影響。通過(guò)研究壓緊力的離散度，來(lái)分析燃料組件壓緊力評(píng)價(jià)方法的準(zhǔn)確性。

3.1 離散度評(píng)價(jià)指標(biāo)

離散度是指變量各個(gè)取值之間的差異程度，可采用極差進(jìn)行衡量。極差,是觀測(cè)變量的最大取值與最小取值之間的離差，也就是觀測(cè)變量的最大觀測(cè)值Max(xi)與最小觀測(cè)值Min(xi)之間的區(qū)間跨度，采用下式進(jìn)行計(jì)算：

R=Max(xi)-Min(xi)

為評(píng)價(jià)壓緊力分析結(jié)果的離散度，采用極差進(jìn)行各方法計(jì)算結(jié)果的離散度分析。

3.2 離散度對(duì)比

從方法原理上講，在計(jì)算壓緊力時(shí)，確定論方法中對(duì)可變參數(shù)取極限值，因此壓緊力的最小值和最大值均比蒙特卡洛方法計(jì)算結(jié)果更保守。確定論方法計(jì)算的壓緊力變化范圍很大，離散度高。采用蒙特卡洛方法計(jì)算的壓緊力變化范圍較小，數(shù)據(jù)更精確，離散度低。

采用確定論方法及蒙特卡洛方法對(duì)燃料組件進(jìn)行壓緊力計(jì)算。對(duì)不同計(jì)算方法下各工況壓緊力的離散度進(jìn)行了研究分析，研究對(duì)象為壽期內(nèi)各時(shí)間點(diǎn)對(duì)應(yīng)工況下壓緊力的極差。

蒙特卡洛方法計(jì)算的壓緊力被包絡(luò)在確定論方法的最小值和最大值之間，與根據(jù)方法原理判斷的預(yù)期是相符的。蒙特卡洛方法計(jì)算的壓緊力的冷態(tài)及熱態(tài)下極差變化相對(duì)確定論方法減小45%～65%，見(jiàn)圖1、圖2。以上對(duì)比數(shù)據(jù)表明采用蒙特卡洛方法可有效降低壓緊力計(jì)算結(jié)果的離散度，提高計(jì)算結(jié)果的準(zhǔn)確性。

圖1 冷態(tài)壓緊力分布及極差變化率Fig.1 Distribution of hold-down force andchange of range in cold state

圖2 熱態(tài)壓緊力分布及極差變化率Fig.2 Distribution of hold-down forceand change of range in hot state

4 結(jié) 論

采用確定論方法計(jì)算得到的燃料組件壓緊力不夠精確，過(guò)于保守。在堆內(nèi)燃耗水平不斷提高的情況下，采用確定論方法對(duì)燃料組件進(jìn)行驗(yàn)證時(shí)，可能會(huì)出現(xiàn)壓緊力不滿(mǎn)足準(zhǔn)則要求的情況。本文基于蒙特卡洛方法對(duì)燃料組件壓緊力進(jìn)行了分析，結(jié)論如下：

1)采用蒙特卡洛方法計(jì)算能夠獲取更大的壓緊力裕量，可避免在實(shí)際工程中由于計(jì)算方法的過(guò)度保守性導(dǎo)致壓緊力不足的問(wèn)題，從而保證燃料組件設(shè)計(jì)滿(mǎn)足工程驗(yàn)證的要求。

2)采用蒙特卡洛方法替代確定論方法進(jìn)行燃料組件壓緊力分析驗(yàn)證，降低了壓緊力分析結(jié)果的離散度，提高了壓緊力計(jì)算的準(zhǔn)確性，能夠?yàn)樵O(shè)計(jì)者設(shè)計(jì)壓緊系統(tǒng)時(shí)避免燃料組件承受過(guò)大的壓緊力提供指導(dǎo)。