常元智棒狀燃料元件的穩(wěn)態(tài)熱工水力分教材

棒狀燃料元件的穩(wěn)態(tài)熱工水力

分析作者姓名：常兀智專業(yè)名稱：核工程與核技術(shù)指導(dǎo)教師：欒峰講師摘要本文為研究棒狀燃料元件所組成的堆芯的熱工水力，計(jì)算出堆芯穩(wěn)態(tài)運(yùn)行時(shí)的相關(guān)參數(shù)，以此估算反應(yīng)堆運(yùn)行時(shí)的安全性。在反應(yīng)堆中，核燃料裂變產(chǎn)生的熱量通過導(dǎo)熱傳遞給燃料元件的外表面，再通過間隙傳遞給包殼內(nèi)表面，再繼續(xù)傳遞給包殼外表面，最后通過對(duì)流換熱將熱量傳遞給冷卻劑。通過計(jì)算得出燃料元件的中心溫度和最高溫度、包殼外表面的最高溫度，還要計(jì)算燃料元件溫度場(chǎng)、冷卻劑溫度場(chǎng)和焓場(chǎng)，并對(duì)結(jié)果進(jìn)行了仔細(xì)的分析。分析結(jié)果表明本文估算所得的相關(guān)參數(shù)與實(shí)際運(yùn)行參數(shù)基本吻合。關(guān)鍵詞：壓水堆核電站熱工水力計(jì)算棒狀燃料元件AbstractInthispaper,thethermalpowerofthecoreoftherodlikefuelelementisstudied,andtherelatedparametersofthecoresteadystateoperationarecalculated,whichcanestimatethesafetyofthereactor.Inthereactor,thenuclearfissiontoproducefuelheatthroughtheheattransfertotheoutersurfaceofthefuelelement,againthroughthegapispassedtothebagontheinnersurfaceoftheshell,tocontinuetotransfertotheoutersurfaceoftheshellpackage,finallybyconvectionheattransfertotransferheattothecoolant.Thetemperaturefield,thetemperaturefieldandtheenthalpyfieldofthefuelelementarecalculatedbycalculatingthemaximumtemperatureandthemaximumtemperatureofthefuelelement,andthetemperaturefieldofthefuelelement.Theanalysisresultsshowthattherelatedparametersarebasicallyconsistentwiththeactualoperatingparameters.Keywords：Pressurizedwaterreactornuclearpowerstation,Thermalhydrauliccalculation,Rodlikefuelelement、八、■

刖言核能在人類生產(chǎn)和生活中的應(yīng)用的主要形式是核電。核電廠具有污染小，發(fā)電成本低等優(yōu)點(diǎn)。結(jié)合我國(guó)國(guó)情，我國(guó)主要依賴煤炭，大量的燃煤帶來了嚴(yán)重的污染問題。盡管采取了脫硫等環(huán)保措施，二氧化硫和氮氧化物的排放總量還是巨大的。加之國(guó)內(nèi)可開發(fā)的水電資源有限，可再生能源等新能源成本高、難以成規(guī)模，環(huán)境狀況非常嚴(yán)峻。在此形式下，發(fā)展核電對(duì)調(diào)整能源結(jié)構(gòu)，減少污染，實(shí)現(xiàn)經(jīng)濟(jì)和生態(tài)環(huán)境協(xié)調(diào)發(fā)展具有十分重要的戰(zhàn)略意義。當(dāng)今世界上運(yùn)行的機(jī)組中大多數(shù)都是壓水堆，在壓水堆堆芯中基本都是采用的棒狀燃料元件，因此對(duì)棒狀燃料元件所組成的堆芯進(jìn)行熱工水力分析是核電廠安全運(yùn)行所必需的。1棒狀燃料元件概述按燃料元件的幾何形狀分類，可將燃料元件分為棒狀、板狀、管狀和球狀等燃料元件形式，本文將主要介紹棒狀燃料元件。對(duì)于不同形狀的燃料元件有不同的分析方法。核燃料裂變產(chǎn)生的熱量通過導(dǎo)熱傳遞給燃料元件的外表面，再通過間隙傳遞給包殼內(nèi)表面，再繼續(xù)傳遞給包殼外表面，最后通過對(duì)流換熱將熱量傳遞給冷卻劑。1.1燃料元件的設(shè)計(jì)要求1.2棒狀燃料元件包殼的作用是防止核燃料與反應(yīng)堆冷卻劑接觸，以免裂變產(chǎn)物逸出對(duì)冷卻劑造成放射性污染。Zr-4合金的中子吸收截面小，在高溫下有較高的機(jī)械強(qiáng)度和抗腐蝕性能。值得提出的是，Zr-4包殼與水相容溫度不超過350C，包殼的熔點(diǎn)為1852C，包殼溫度達(dá)到820E后鋯與水反應(yīng)產(chǎn)生氫氣并放出大量的熱。在運(yùn)行中應(yīng)使燃料元件保持在適度的溫度，以免包殼與冷卻劑之間發(fā)生相互作用。1.3包殼材料在堆芯壽期內(nèi)，核燃料裂變所產(chǎn)生的能量將傳遞給冷卻劑。假如燃料芯塊直接與冷卻劑相接觸，那么由于核燃料裂變產(chǎn)生廢物將直接進(jìn)入冷卻劑，從而造成冷卻劑具有強(qiáng)放射性，對(duì)核電廠的維護(hù)和人員安全造成威脅，所以在燃料芯塊和冷卻劑之間加入一層非裂變材料即包殼是非常必要的。通過以上敘述可總結(jié)包殼的作用：使燃料芯塊不受冷卻劑的機(jī)械侵蝕和化學(xué)腐蝕。包容由于核燃料裂變產(chǎn)生的裂變氣體和其他裂變產(chǎn)物。封裝燃料芯塊，有支撐燃料芯塊的幾何形狀。對(duì)于包殼材料的選擇，因考慮一下因素：具有較好的核性能，即中子吸收截面小，感生放射性小。與燃料芯塊有較好的相容性，且耐高溫。具有良好的傳熱的性能。具有較好的力學(xué)性能。具有良好的抗腐能力，對(duì)冷卻劑是惰性的。具有穩(wěn)定的抗輻照性能。容易加工，成本低廉，便于后處理。綜合考慮上述要求，可用作包殼材料有：鎂、鋁、鋯、鈹、石墨、不銹鋼、鎳基合金等。由于本文所討論包殼選用Zr-4合金，下我們重點(diǎn)討論Zr-4的性質(zhì)。Zr-4中合金所含的微量元素如表1.1所示、、微量元素鋯合金\SnFeCrNiZr-41.5007表1.1Zr-4合金中微量元素（%）1.熱導(dǎo)率。Zr-4合金的熱導(dǎo)率在表1.2體現(xiàn)t/C100300500700K/WcmCZr-40.1360.1520.1800.225表1.2Zr-4的熱導(dǎo)率（試樣在750c下經(jīng)20h真空退火）Zr-4合金的熱導(dǎo)率關(guān)于溫度的函數(shù)表示為式(2-5):■?C=7.7310—3.1510%-2.8710Jt21.55210」°t3(2-5)式中：'c是Zr-4合金的熱導(dǎo)率，W(cm-C)；t是溫度，C。2.比定壓熱容。Zr-4合金的比定壓熱容關(guān)于溫度變化的關(guān)系表示為式(2-6):c=286.50.1t,0::t::750Cp(2-6cp二360,t750C式中：Cp是Zr-4的比定壓熱容(j/(kgC))；t是溫度(C)2棒狀燃料元件的傳熱分析燃料元件傳熱分析在反應(yīng)堆熱工分析中有著重要的地位。之所以要對(duì)燃料元件進(jìn)行傳熱分析，是出于對(duì)反應(yīng)堆的安全性和經(jīng)濟(jì)性兩個(gè)角度考慮的。首先，從安全性來講：第一，要確保在任何工況下不發(fā)生燃料元件熔化；第二，在燃料芯塊和結(jié)構(gòu)材料設(shè)計(jì)是要考慮由于溫度梯度會(huì)造成熱應(yīng)力；第三，包殼和冷卻劑的化學(xué)反應(yīng)也與溫度相關(guān)；第四，燃料和慢化劑由于溫度變化而引入反應(yīng)性的變化會(huì)影響到反應(yīng)堆的控制。再者，在保證反應(yīng)堆安全的前提下，要盡可能通過提高堆芯單位體積的熱功率，冷卻劑的溫度等方式來提高反應(yīng)堆的經(jīng)濟(jì)性。2.1熱傳導(dǎo)微分方程燃料元件的導(dǎo)熱是通過熱傳導(dǎo)的方式將燃料元件內(nèi)產(chǎn)生的裂變熱，從溫度高的芯塊傳遞到向溫度低的包殼外表面的過程。此過程的導(dǎo)熱方程便可描述成式(2-1):'Cpr,'r,tItr,.qvr,.(2-1)式中：，為密度氽9m3)；cp(r,)為.時(shí)刻r處的比定壓熱容(J(kg-C))；'(r,t)為tC時(shí)；處的熱導(dǎo)率(W(cm-C))；t(r*,)為.時(shí)刻r處的溫度(C)qv(r,)為.時(shí)刻r處的體積是熱率(kWm3)。本文對(duì)穩(wěn)定工況的燃料元件分析，得到式(2-1)穩(wěn)態(tài)形式(2-2):

(2-2)八r*,ttrqvr(2-2)我們已經(jīng)知道，熱導(dǎo)率■是關(guān)于溫度t的函數(shù)。由于在此處我們認(rèn)為導(dǎo)熱過程是主要矛盾，而熱導(dǎo)率■的變化只是次要矛盾，所以假設(shè)熱導(dǎo)率’是一個(gè)常數(shù)，得到式(2-3):-4才vt2rqr=0(2-3)2.2有內(nèi)熱源柱形芯塊的溫度場(chǎng)2.2有內(nèi)熱源柱形芯塊的溫度場(chǎng)先做兩點(diǎn)合理假設(shè)：忽略軸線方向的導(dǎo)熱，認(rèn)為導(dǎo)熱只在徑向進(jìn)行。認(rèn)為徑向的體積釋熱率是常數(shù)?；谝陨蟽牲c(diǎn)假設(shè)，將棒狀芯塊的導(dǎo)熱問題轉(zhuǎn)化為一維問題。棒狀燃料元件采用一維圓柱形坐標(biāo)系就有(2-4):(2-4)2孝」，4毛drrdru(2-4)如圖2.1，設(shè)燃料芯塊半徑為尺，包殼的外半徑為Rc,包殼的厚度為：，先不考慮氣隙的厚度，三者之間有關(guān)系式(2-5Rc二Ru、(2-5)圖2.1棒狀燃料元件橫截面現(xiàn)求解微分方程式（2-4）:2=0「dr2'u(2-6)由假設(shè)（1）可知芯塊的溫度分布具有對(duì)稱性:tit=0drr=0得到G=0與式（2-8）:史［蟲r=0dr2\(2-8)再對(duì)式（2-8）進(jìn)行一次積分，得到式（2-9）:1qvr2(2-9)將邊界條件，r=R,t=tu代入式（2-9）得到C2和式（2-10）：2.3無內(nèi)熱源筒形包殼的溫度場(chǎng)首先引入傅里葉導(dǎo)熱定律的一維形式:Q=cQ=cA—-J1dr（2-11）式中：Q為通過圓筒形包殼的總功率（W）；瓷c為包殼的熱導(dǎo)率（Wcm；C）；A為垂直于導(dǎo)熱方向的面積（m2）。雖然包殼會(huì)吸收來自燃料芯塊的射線、1射線等而產(chǎn)生熱量，但由于包殼的厚度太小，在其中產(chǎn)生的熱量相對(duì)于燃料芯塊傳遞給包殼的熱量可以忽略。故可以把包殼導(dǎo)的熱處理成無內(nèi)熱源的導(dǎo)熱過程，便將式（2-3）寫成式（2-12）:2彳八(2-13)八2tr=0（2-12）Jc回2將式（2-12）在圓柱形坐標(biāo)系中改寫成式（2-13）:(2-13)1dtr邛津卜。dr..dr.j—（2-14）求解微分方程式（2-13）:現(xiàn)將傅里葉導(dǎo)熱定律的一維形式根據(jù)圓柱形包殼改寫成式drG(2-15):(2-15)式中：A=2二|。故(2-16)得到(2-16)12二I將式（2-16）中的Q定義為線功率，用qj表示。故式（2-14）改寫成式（2-17）：1dtqii（2-17）dr2二cr再對(duì)式（2-17）進(jìn)行一次積分，得到式（2-18）:"J'C2川％eg）將邊界條件，r譏，5代入式（2-18）得到C2和式（2-19）:C2trtuin（2-19）U2叭Ru將包殼的半徑Rc代入式（2-19）得到包殼外表面的溫度表達(dá)式（2-20）:tc鞏n昱C(2-20)2.4氣隙熱傳導(dǎo)從圖2.2可以看到，燃料棒內(nèi)除了二氧化鈾，氣隙能夠引起明顯溫度變化。由此可見，氣隙傳熱計(jì)算的準(zhǔn)確程度，將會(huì)很大程度影響燃料芯塊溫度計(jì)算的準(zhǔn)確性。氣隙是由惰性氣體氦占據(jù)的環(huán)形空間。氣隙內(nèi)的氣體成分會(huì)隨燃耗加深以及裂變的氣體增加而改變。隨著燃耗的加深，燃料芯塊因?yàn)檩椪债a(chǎn)生腫脹和變形導(dǎo)致與包殼接觸，這就

型就比較合適本文穩(wěn)定工況條件下的氣隙傳熱求解。此處我們認(rèn)為燃料芯塊與包殼恰好接觸，且接觸壓力為零。經(jīng)過這樣近似處理后可以應(yīng)用接觸導(dǎo)熱的等效傳熱系，此=5768W,'(m2LG。對(duì)于接觸傳熱傳熱模型，我們應(yīng)用牛頓冷卻公式(Newton'slawofcooling)求解內(nèi)包殼的溫度：圖2.2典型輕水堆燃料棒的溫度分布-492W/om1760159314261260IG9392760■593q-492W/om=t?2土(2-21)Rh圖2.2典型輕水堆燃料棒的溫度分布2.5包殼外表面與冷卻劑的對(duì)流換熱換熱過程是一個(gè)熱量由包殼外表面?zhèn)鬟f給冷卻劑的過程，而本文所討論的壓水堆一回路是不允許出現(xiàn)相變的，致使冷卻劑流

過包殼的換熱過程是單向?qū)α鲹Q熱，因此包殼向冷卻劑傳遞的熱量可以應(yīng)用牛頓冷卻公式（Newton'slawofcooling求得：（2-22）Q=hAJo）式中：Q為熱功率（W）；h為對(duì)流傳熱系數(shù)（Wf（m2|_C））；A為傳熱面積（A為傳熱面積（m2）;tcs為包殼外表面溫度：co為冷卻劑主流溫度（C）。ql2二Rh(2-24)將式（2-23）改寫成線功率的形式：ql2二Rh式中：ql是線功率（Wm）；Rcs為包殼外徑（m）。若要求解包殼外表面流度，需要先確定式（2-25）中的對(duì)流換熱系數(shù)h，對(duì)于采用棒束燃料組件的壓水堆，采用Weisman的經(jīng)驗(yàn)關(guān)系式：TOC\o"1-5"\h\zNu=CRe0.8Pr13（2-25）式中C是一個(gè)常數(shù)，它取決于柵格的排列形式。本文討論的壓水堆使用的是正方形柵格：PP當(dāng)1.11.3時(shí)，C=0.0420.024。dd

2.6棒狀燃料元件的溫度分布如圖2.3,當(dāng)冷卻劑沿z方向流動(dòng)時(shí)，累積吸收來自燃料芯塊裂變產(chǎn)生的熱能，使其溫度不斷升高，故溫度t不僅要考慮徑向的分布，還要準(zhǔn)確描述軸向z的分布。圖2.3圖2.3燃料元件換熱示意圖首先給出輸熱量關(guān)于z的函數(shù):(2-26)(2-26)Qz二WqI。z式中：Q（z）為冷卻劑在堆芯z處所傳出的熱量，（W）；W為冷卻劑質(zhì)量流量（kgs）；Cp為冷卻劑的比定壓熱容（J；（kg上））；tcoz為冷卻劑從入口到堆芯z處的溫度變化（C）。將式（2-26）改寫成（2-27）:Q（z）LoZ^co,inWCP式中：tcoin為冷卻劑入口溫度（C）。如圖2.4所示，以棒狀燃料元件的半高處為原點(diǎn)，建立線功率q與z的直角坐標(biāo)系圖2.4圖2.4棒狀燃料元件上建立的q-z坐標(biāo)系(2-28)先給出輸熱量Qz與線功率qz的關(guān)系：

z

Qz一嚴(yán)(2-28)_2式中：H為棒狀燃料的高度（m）。若該棒狀元件的軸向發(fā)熱按照余弦分布：(2-29)_,/(2-29)q.=Q.（0）cosj<HReJ式中：q.0為在z=0處的線功率（Wm）。聯(lián)立式（2-27）、式（2-28）、式（2-29）得到在z的冷卻劑tz二t%樵*sinzsin兀Hr（2-30）co河n.WcpHRe2HR?！箿囟萾coz:將代入式（2-31），得到出口溫度tco,ex將式(2-32)將式(2-24)將式(2-29)將式(2-36)t■-七卜co,ex2q(°叫sinCP將式（2-31）寫成進(jìn)出口溫差.-：代入式（2-30）:cP改寫成關(guān)于z的函數(shù):代入式（2-31）:tcsZ-ttcoR2H二的形式:H—sin兀WCpR2Hr:」(2-31)(2-32)\o2cozcsz-tcoz二C02二2兀Rsh（z）在z=0處冷卻劑與外包殼的溫差為:cs0=tcs0(2-33)EeRcshzcos■:zHReJ(2-34)(2-35)tco0=ARcshz(2-36)、式（2-35）、式（2-33）聯(lián)立得：七今也譜竺sintcszino)coscs((2-3CPHReJ在推導(dǎo)式（2-20）時(shí)，未考慮氣隙的存在7）現(xiàn)將式其改寫成考tciz=tcsn八Rci(2-38)慮氣隙且關(guān)于z的函數(shù):2，c將式（2-29）代入式（2-38）:tciz=tReJ(2-39)(2-40)(2-41)將式tciz%z(2-40)(2-41)將式tciz%z怙0cos兀zARe(2-21)改寫成關(guān)于z的函數(shù):q(z)tuz口寸z(2-42)(2-40)代入式(2-39):將式將式(2-29)代入式(2-42):七吃亦q(0)h二(2-43)tzzcOsRe」Ui0二tci0-tcs0=q0n字070七。^1寸寸(2-44)g在z=0處燃料溫差為：將式(2-44)代入式(2-43):tu(z)=tci(z)+%(0)cos—/、jj<tu(z)=tci(z)+%(0)cos—/、jj<HRe)j1H(2-45)將r=0代入式(2-10):to=tR24R(2-46)將式(2-46)中體積釋熱率qv中線功率q替換并且寫成關(guān)于z的函數(shù):(2-47)tozUz+(2-47)q(0)cost。z7z"4匚將式（2-29）代入式（2-47）:在z=0處燃料芯塊中心和燃料芯塊表面溫差為:「（-0“00

u將式（2-49）代入式（2-48）:0—tu(2-48)(2-49)(2-50)toz譏z%zcos聯(lián)立式（2-37）（2-41）、（2-45）、（2-50）得到燃料芯t0zJJnT，sin+巨也(0)1cos(2-51)HRe)".)塊中心溫度：(2-51)其中：—l0「cs0*ci0Tg0"u03燃料元件的傳熱計(jì)算在反應(yīng)堆中，核燃料裂變產(chǎn)生的熱量通過導(dǎo)熱傳遞給燃料元件的外表面，再通過氣隙傳遞給包殼內(nèi)表面，再繼續(xù)傳遞給包殼外表面，最后通過對(duì)流換熱將熱量傳遞給冷卻劑。而熱量通過包殼外表面?zhèn)鬟f給冷卻劑時(shí)，可能發(fā)生多種傳熱方式，比如單相對(duì)流換熱，流動(dòng)沸騰傳熱。具體是哪種，要根據(jù)包殼外表面的各種參數(shù)確定。在本文中我們假定只發(fā)生單相對(duì)流換熱，而不考慮流動(dòng)沸騰傳熱。3.1各參數(shù)的獲得3.1.1最大線功率的獲得現(xiàn)已知某核電廠堆芯額定熱功率為1930MW,而在堆芯堆芯裂變產(chǎn)生熱的過程中，大約有2.6%的熱功率來自于反應(yīng)堆冷卻劑中射線與流體的相互作用，故額定熱功率的97.4%的熱功率來自于核燃料的裂變。燃料中產(chǎn)生的功率：193097.4%=1879.82MW燃料總長(zhǎng)度(燃料組件數(shù)每個(gè)燃料組件的元件數(shù)每根元件的長(zhǎng)度)：平均線功率121264368.1=1.176107cmq:1879.8210A-1.176107“10=1.60104Wm也’心z)z?匚占q(0)cosIdz(3-1)2(3-1)■RTRe那么由平均線功率可以解得最大線功率：變形得：代入數(shù)據(jù)，得到訂(3-2)說明：雖然堆芯」中的燃料所在區(qū)域不q0同富集度不相同，但這種分區(qū)裝料的方式是為了將功率展平，故在此處使用平均線功率是有意義的。3.1.2對(duì)流換熱系數(shù)的獲得該核電廠的正常運(yùn)行壓力為15.5MPa,反應(yīng)堆平均溫度310C,根據(jù)表3.1獲得該壓力和溫度下的比定壓熱容。t/Cp/MPa、\28032015.05.0646.12816.05.0446.056表3.1水在不同壓力和溫度下的比定壓熱容Cp.:(KJKgl_C)對(duì)表3.1使用內(nèi)插法計(jì)算，解得15.5MPa,310C下的比定壓熱容:5.833103J..KgLC。同理得到15.5MPa,310C下的比體積v=1.41910"m3「Kg;運(yùn)動(dòng)粘度j=8.8410*Pag;熱導(dǎo)率=0.552W.mUC。根據(jù)以上參數(shù)可以求得普朗特?cái)?shù)Pu0.934運(yùn)動(dòng)粘度=1.254310mis。已知該核電廠所用燃料包殼直徑為d=9.6mm，燃料棒間距P=12.5mm,冷卻劑流速V=3.5ms，現(xiàn)求得等效直徑：-4P2:

.TOC\o"1-5"\h\zD二(3-3)e一將數(shù)據(jù)代入式(3-3)，得到De=1.1M10‘m?，F(xiàn)已知V、De、、??，便可求得雷諾數(shù)Re:6VDeRee(3-4)將數(shù)據(jù)代入式(3-4)，得到Re=3.097勺05根據(jù)柵距和棒徑比有Pd=1.3，在式(2-25)的適用范圍內(nèi)，P故可求出C=0.42匚-0.024=0.0306。d現(xiàn)已知C、Re、Pr，代入式(2-25)得努塞爾數(shù)Nu=738.9(3-5)根據(jù)努塞爾數(shù)Nu的定義式：解得對(duì)流傳熱系數(shù)h=3.674104W「m2l_C。說明：相關(guān)資料表明，對(duì)流傳熱系數(shù)沿z方向大不，故將其做常數(shù)處理，即hz二h=3.67104W..；m2l_C。3.1.2包殼熱導(dǎo)率與燃料芯塊熱導(dǎo)率的獲得包殼熱導(dǎo)率■c關(guān)于溫度t的變化如圖3.1所示：圖3.1包殼熱導(dǎo)率與溫度的關(guān)系觀察圖3.1可知道熱導(dǎo)率在此溫度范圍內(nèi)變化不大，故作平均處理：t2t'cdt(3-6)K

c

t1聯(lián)立式（2-5）、式（3-6），并代入數(shù)據(jù)解得心=15.38W/mic燃料芯塊熱導(dǎo)率*關(guān)于溫度t的變化如圖3.2所示4a300400500SOQ7004a300400500SOQ700BOO9401DOO觀察圖3.2可知道熱導(dǎo)率在此溫度范圍內(nèi)變化不大，故同樣作平均處理：（3-7）聯(lián)立式（2-5）、式（3-6），并代入數(shù)據(jù)解得畸=3.71W/m°C。3.2燃料芯塊中心溫度及其最高溫度該核電廠堆芯入口溫度293C，出口溫度為327C。燃料元件的幾何參數(shù)如表3.2所示。燃料包殼外徑燃料包殼內(nèi)徑燃料包殼厚度燃料芯塊直徑燃料元件長(zhǎng)度58.353681表3.2燃料元件幾何參數(shù)表mm有等效直徑De、比體積v、流速V估算流過燃料元件的冷卻劑質(zhì)量流量W=0.318Kgs。將數(shù)據(jù)代入各式得到億0=536.、2.-:g0=165.2、cs0=22.6、0=31.5，故'0=755.5。ci代入數(shù)據(jù)，式（2-51）可以寫成:toz=31015.86sin0.85z755.5cos0.85z(3-8)現(xiàn)來求解燃料芯塊最高溫度:dtoz13.481cos0.85z-642.175sin0.85ztoz=31015.86sin0.85z755.5cos0.85z(3-8)現(xiàn)來求解燃料芯塊最高溫度:dtoz13.481cos0.85z-642.175sin0.85z=0dz