核電廠疲勞監(jiān)測(cè)系統(tǒng)管壁導(dǎo)熱反演數(shù)值計(jì)算方法

劉　浪，孟阿軍，凌　君，何大宇

（深圳中廣核工程設(shè)計(jì)有限公司，廣東深圳　518172）

核電廠疲勞監(jiān)測(cè)系統(tǒng)管壁導(dǎo)熱反演數(shù)值計(jì)算方法

劉浪，孟阿軍，凌君，何大宇

（深圳中廣核工程設(shè)計(jì)有限公司，廣東深圳518172）

核電廠疲勞監(jiān)測(cè)系統(tǒng)中的導(dǎo)熱反演計(jì)算是關(guān)鍵步驟，首先推演得到一維非穩(wěn)態(tài)導(dǎo)熱反演問(wèn)題的數(shù)值計(jì)算公式，然后基于該公式，利用試驗(yàn)獲取的管道外壁實(shí)測(cè)溫度時(shí)程數(shù)據(jù)計(jì)算得到管道內(nèi)壁溫度，最后通過(guò)與試驗(yàn)實(shí)測(cè)管道內(nèi)壁溫度進(jìn)行比較，驗(yàn)證了一維非穩(wěn)態(tài)導(dǎo)熱反演問(wèn)題的數(shù)值計(jì)算方法的正確性。

疲勞監(jiān)測(cè)系統(tǒng)；一維非穩(wěn)態(tài)導(dǎo)熱反演；數(shù)值計(jì)算

符號(hào)說(shuō)明：

Φλ（i-1）——控制容積i-1導(dǎo)入i的熱量，W

d U——控制容積i熱力學(xué)能的增加，W

A——控制容積橫截面面積，m2

λ——熱導(dǎo)率，W/（m·℃）

Δx——控制容積長(zhǎng)度，即空間離散步長(zhǎng)，m

Δτ——時(shí)間離散步長(zhǎng)，s

ρ——密度，kg/m3

c——比熱容，J/（kg·℃）

ti

s——控制容積i在s時(shí)刻的溫度，℃

t∞

s——環(huán)境流場(chǎng)在s時(shí)刻的溫度，℃

a——熱擴(kuò)散率，m2/s，a=λ/（ρc）

To（i，j）——管道外壁的溫度，℃

Ti（i，j）——管道內(nèi)壁的溫度，℃

0　引言

在核電廠老化與壽命管理的理論體系中，關(guān)鍵設(shè)備老化降質(zhì)的監(jiān)督監(jiān)測(cè)結(jié)果是后續(xù)開(kāi)展設(shè)備老化狀態(tài)和壽命評(píng)估的基礎(chǔ)。核電廠關(guān)鍵設(shè)備存在的主要老化機(jī)理包括疲勞、腐蝕、磨損、輻照脆化以及熱老化等，其中疲勞是造成核電廠關(guān)鍵設(shè)備失效最主要的原因之一。目前，在世界范圍內(nèi)核電廠已經(jīng)發(fā)生了一系列反應(yīng)堆冷卻劑管道與疲勞相關(guān)的事件，已影響到電廠安全、經(jīng)濟(jì)運(yùn)行，嚴(yán)重時(shí)甚至引起非計(jì)劃停堆。

通常，在核電廠的工程設(shè)計(jì)階段，會(huì)按照規(guī)范要求對(duì)核一級(jí)管道進(jìn)行疲勞分析，但是該分析結(jié)果不能真實(shí)反映管道的實(shí)際損耗，因?yàn)樵摲治鲞^(guò)程中使用的設(shè)計(jì)瞬態(tài)進(jìn)行計(jì)算，然而電廠實(shí)際運(yùn)行瞬態(tài)與設(shè)計(jì)瞬態(tài)在數(shù)值上相比有一定的差距。為了獲取實(shí)際運(yùn)行瞬態(tài)數(shù)據(jù)，掌握管道的真實(shí)疲勞狀態(tài)，基于設(shè)備的實(shí)際損傷狀態(tài)，優(yōu)化運(yùn)行規(guī)程及在役檢查大綱，合理挖掘關(guān)鍵管道疲勞設(shè)計(jì)的安全裕度，為電廠定期安全審查（PSR）或電廠延壽提供真實(shí)的數(shù)據(jù)支撐，部分核電國(guó)家開(kāi)展了核電廠疲勞監(jiān)測(cè)系統(tǒng)的研發(fā)工作［1］。

一般的管道計(jì)算軟件如SYSPIPE，PIPESTRESS可以用設(shè)計(jì)瞬態(tài)數(shù)據(jù)開(kāi)展疲勞分析得到使用系數(shù)，但是上述軟件無(wú)法應(yīng)用到疲勞監(jiān)測(cè)系統(tǒng)，原因有兩個(gè):一方面，實(shí)際運(yùn)行瞬態(tài)與設(shè)計(jì)瞬態(tài)相比在數(shù)據(jù)量上差別非常大，SYSPIPE，PIPESTRESS只能計(jì)算有限的瞬態(tài)個(gè)數(shù)；另一方面，設(shè)計(jì)瞬態(tài)描述的是工質(zhì)溫度變化，而疲勞監(jiān)測(cè)系統(tǒng)實(shí)測(cè)到的是管道外壁溫度變化，SYSPIPE，PIPESTRESS軟件不具備導(dǎo)熱反演計(jì)算的能力。綜上所述，需要專門(mén)開(kāi)發(fā)導(dǎo)熱反演的計(jì)算方法，滿足疲勞監(jiān)測(cè)系統(tǒng)的需要［2-6］。

本文基于試驗(yàn)數(shù)據(jù)，在研究導(dǎo)熱數(shù)值反演計(jì)算過(guò)程中，重點(diǎn)分析了傳熱數(shù)值模型、邊界條件、內(nèi)部節(jié)點(diǎn)和邊界節(jié)點(diǎn)的溫度差分方程、空間離散度對(duì)計(jì)算精度的影響。

1　分析流程

通過(guò)疲勞監(jiān)測(cè)系統(tǒng)計(jì)算得到疲勞使用系數(shù)需要以下4個(gè)步驟，本文重點(diǎn)研究第2步工作。

第1步:通過(guò)疲勞監(jiān)測(cè)系統(tǒng)在管道外壁安裝的熱電偶（如圖1，2所示），實(shí)測(cè)獲取管道外壁溫度；

第2步:根據(jù)上一步實(shí)測(cè)管道外壁溫度，導(dǎo)熱數(shù)值反演計(jì)算得到管壁徑向各點(diǎn)溫度；

第3步:計(jì)算管壁徑向溫度梯度，包括平均溫度、線性溫度梯度和非線性溫度梯度；

第4步:將溫度梯度代入RCCM公式計(jì)算應(yīng)力，最終求得疲勞使用系數(shù)。

圖1　監(jiān)測(cè)點(diǎn)示意

圖2　監(jiān)測(cè)點(diǎn)橫截面圖

本文重點(diǎn)研究第2步工作，基于試驗(yàn)（詳見(jiàn)第3節(jié)）數(shù)據(jù)（包括實(shí)際測(cè)量的管道外壁和內(nèi)壁溫度時(shí)程數(shù)據(jù)），首先驗(yàn)證試驗(yàn)數(shù)據(jù)的準(zhǔn)確性（詳見(jiàn)第3.3節(jié)），即用ANSYS軟件依據(jù)實(shí)測(cè)內(nèi)壁溫度計(jì)算得到外壁溫度，并與實(shí)測(cè)外壁溫度比較，判斷試驗(yàn)數(shù)據(jù)是否準(zhǔn)確；其次，驗(yàn)證一維非穩(wěn)態(tài)導(dǎo)熱反演問(wèn)題的數(shù)值計(jì)算方法的正確性（詳見(jiàn)第3.4節(jié)），即依據(jù)實(shí)測(cè)外壁溫度數(shù)值反演計(jì)算得到內(nèi)壁溫度，并與實(shí)測(cè)內(nèi)壁溫度比較，判斷計(jì)算精確度；然后分析壁厚的最小空間離散度（詳見(jiàn)第3.5節(jié)），即用ANSYS軟件計(jì)算不同空間離散度下壁厚的溫度梯度，得到最小空間離散度，也就是說(shuō)壁厚徑向最少需要幾個(gè)節(jié)點(diǎn)；隨后用一維非穩(wěn)態(tài)導(dǎo)熱反演問(wèn)題的數(shù)值計(jì)算方法類似第2步求的各節(jié)點(diǎn)溫度（詳見(jiàn)第3.6節(jié)），最后求得壁厚溫度梯度（詳見(jiàn)第3.7節(jié)），并與ANSYS軟件計(jì)算的溫度梯度比較，論證數(shù)值計(jì)算結(jié)果的正確性。

2　分析方法

2.1導(dǎo)熱數(shù)值反演

2.1.1基本假設(shè)

核電廠被監(jiān)測(cè)的熱載荷主要有兩種類型:熱沖擊和熱分層，如圖3，4所示。

圖3　熱沖擊示意

圖4　熱分層示意

在熱沖擊的情況下管道周向溫度沒(méi)有差異，在熱分層的情況下管道周向溫度有差異，但是對(duì)于薄壁管道而言，監(jiān)測(cè)點(diǎn)處的導(dǎo)熱計(jì)算可以基本假設(shè)如下:

（1）管道軸向傳熱忽略；

（2）管道為薄壁軸對(duì)稱模型，將分析對(duì)象簡(jiǎn)化為大平壁一維熱傳導(dǎo)問(wèn)題。

2.1.2溫度差分方程

（1）內(nèi)部節(jié)點(diǎn)溫度差分方程［7］。

如圖5所示，對(duì)于無(wú)內(nèi)熱源的大平壁的一維非穩(wěn)態(tài)導(dǎo)熱，內(nèi)部節(jié)點(diǎn)i所代表的控制容積（灰色）的熱平衡可表述為:在s時(shí)刻，單位時(shí)間內(nèi)從相鄰控制容積i-1與i+1分別導(dǎo)入的熱量Φλ（i-1）與Φλ（i+1）之和等于該控制容積熱力學(xué)能的增加d U，即:

圖5　一維非穩(wěn)態(tài)導(dǎo)熱內(nèi)部節(jié)點(diǎn)溫度差分方程

節(jié)點(diǎn)i的溫度對(duì)時(shí)間的變化率采用向前差分，則熱平衡方程為:

令FOΔ=AΔτ/Δx2，則上式可以整理為:

針對(duì)上式說(shuō)明兩點(diǎn):

1）內(nèi)部節(jié)點(diǎn)i在s時(shí)刻的溫度，可以由該節(jié)點(diǎn)及相鄰節(jié)點(diǎn)在s-1時(shí)刻的節(jié)點(diǎn)溫度求出；

2）上式1-2FOΔ不能為負(fù)，稱為一維非穩(wěn)態(tài)導(dǎo)熱內(nèi)部節(jié)點(diǎn)溫度方程的穩(wěn)定性條件。

（2）邊界節(jié)點(diǎn)溫度差分方程［7］。

對(duì)于邊界節(jié)點(diǎn)（如圖6所示）的導(dǎo)熱計(jì)算，同上述過(guò)程，可以得到邊界節(jié)點(diǎn)i-1的熱平衡方程為:

令FOΔ=aΔτ/Δx2，BiΔ=hΔx/λ，則上式可以整理為:

針對(duì)上式說(shuō)明兩點(diǎn):

1）邊界節(jié)點(diǎn)i在s+1時(shí)刻的溫度，可以由該節(jié)點(diǎn)及相鄰節(jié)點(diǎn)在s時(shí)刻的節(jié)點(diǎn)溫度求出；

2）上式1-2BiΔFOΔ-2FOΔ不能為負(fù)，稱為一維非穩(wěn)態(tài)導(dǎo)熱邊界節(jié)點(diǎn)溫度方程的穩(wěn)定性條件。

圖6　一維非穩(wěn)態(tài)導(dǎo)熱邊界節(jié)點(diǎn)溫度差分方程

（3）邊界節(jié)點(diǎn)溫度差分方程反演。

根據(jù)式（5）可以得到如下方程:

針對(duì)上式說(shuō)明兩點(diǎn):

1）內(nèi)部節(jié)點(diǎn)i在s時(shí)刻的溫度，可以由相鄰邊界節(jié)點(diǎn)在s，s+1時(shí)刻的節(jié)點(diǎn)溫度求出；

2）上式1-2BiΔFOΔ-2FOΔ不能為負(fù)，稱為一維非穩(wěn)態(tài)導(dǎo)熱邊界節(jié)點(diǎn)溫度反演方程的穩(wěn)定性條件。2.2溫度梯度計(jì)算［8］

對(duì)管道進(jìn)行一維非穩(wěn)態(tài)導(dǎo)熱分析后得到管壁徑向的溫度分布（見(jiàn)圖7），然后根據(jù)RCC-M規(guī)范對(duì)管壁徑向溫度分布結(jié)果進(jìn)行處理以得到溫度梯度分布，即求得沿管壁徑向平均溫度分布T（i，j）、線性溫度分布ΔT1（i，j）以及非線性溫度分布ΔT2（i，j）（其中i，j分別表示載荷i與j）。本文計(jì)算溫度梯度的目的是分析在一維非穩(wěn)態(tài)導(dǎo)熱數(shù)值計(jì)算過(guò)程中空間離散尺度對(duì)疲勞計(jì)算結(jié)果的影響。

圖7　管壁徑向溫度分布

（1）平均溫度T（i，j）。

（2）線性分布的溫度ΔT1（i，j）。

（3）非線性分布的溫度ΔT2（i，j）。

3　分析實(shí)例

3.1試驗(yàn)?zāi)Ｐ徒榻B

試驗(yàn)?zāi)Ｐ腿鐖D8所示。試驗(yàn)過(guò)程:第1步，001VP，002VP關(guān)閉，001BA內(nèi)流質(zhì)加熱至75℃，同時(shí)003VP，004VP打開(kāi)，002PO運(yùn)轉(zhuǎn)，用002BA中的常溫流質(zhì)冷卻監(jiān)測(cè)點(diǎn)管段至常溫；第2步，002PO停運(yùn)，003VP，004VP關(guān)閉，001VP，002VP打開(kāi)，001PO運(yùn)轉(zhuǎn)，001BA中的加熱流質(zhì)進(jìn)入監(jiān)測(cè)點(diǎn)管段，形成熱沖擊，觀察探測(cè)到的管道外壁溫度數(shù)據(jù)，待溫度監(jiān)測(cè)數(shù)據(jù)進(jìn)入穩(wěn)態(tài)后，001PO停運(yùn)，計(jì)算并分析監(jiān)測(cè)結(jié)果。

圖8　試驗(yàn)?zāi)Ｐ?/p>

圖9　監(jiān)測(cè)點(diǎn)橫截面圖

如圖9所示，監(jiān)測(cè)點(diǎn)處沿管壁徑向，在管道內(nèi)壁、外壁設(shè)置（焊接）熱電偶，測(cè)量管壁溫度變化；管道帶保溫，管壁厚9mm，管道材料Z2CN18.10，材料性能數(shù)據(jù)見(jiàn)表1。

表1　材料性能

3.2試驗(yàn)監(jiān)測(cè)數(shù)據(jù)

如上所述，試驗(yàn)可以監(jiān)測(cè)到管道內(nèi)壁、外壁的溫度數(shù)據(jù)，如圖10所示。以下將根據(jù)一組試驗(yàn)監(jiān)測(cè)到的管道外壁溫度數(shù)據(jù)，進(jìn)行一維非穩(wěn)態(tài)導(dǎo)熱反演的數(shù)值計(jì)算，算出管道內(nèi)壁的溫度，然后與試驗(yàn)監(jiān)測(cè)的數(shù)據(jù)（如圖10所示）進(jìn)行對(duì)比，論證計(jì)算方法的正確性。但是在開(kāi)展反演計(jì)算前，首先需要分析監(jiān)測(cè)數(shù)據(jù)（如圖10所示）的準(zhǔn)確度（詳見(jiàn)第3.3節(jié)），因?yàn)楸O(jiān)測(cè)點(diǎn)處，焊接于管道外壁和內(nèi)壁的熱電偶會(huì)因?yàn)楹附淤|(zhì)量問(wèn)題導(dǎo)致監(jiān)測(cè)數(shù)據(jù)失真。

圖10　試驗(yàn)監(jiān)測(cè)溫度數(shù)據(jù)

3.3試驗(yàn)監(jiān)測(cè)數(shù)據(jù)準(zhǔn)確性分析

用ANSYS軟件進(jìn)行熱傳導(dǎo)計(jì)算，根據(jù)管道截面實(shí)際尺寸（外徑和壁厚）確定矩形截面，依據(jù)實(shí)測(cè)內(nèi)壁溫度計(jì)算得到外壁溫度，外壁有保溫層，所以邊界條件按絕熱處理，ANSYS軟件計(jì)算模型見(jiàn)圖11，計(jì)算結(jié)果見(jiàn)圖12，可以發(fā)現(xiàn)計(jì)算的外壁溫度與實(shí)測(cè)外壁溫度吻合，說(shuō)明試驗(yàn)獲取的管道內(nèi)壁和外壁溫度數(shù)據(jù)準(zhǔn)確［9］。

圖11　有限元模型

圖12　試驗(yàn)監(jiān)測(cè)溫度數(shù)據(jù)驗(yàn)證

3.4一維非穩(wěn)態(tài)導(dǎo)熱反演數(shù)值計(jì)算方法正確性分析

基于試驗(yàn)實(shí)測(cè)外壁溫度，根據(jù)第2.1.2（3）節(jié)中邊界點(diǎn)溫度差分方程反演計(jì)算管道內(nèi)壁溫度，與試驗(yàn)實(shí)測(cè)內(nèi)壁溫度比較（如圖13所示），發(fā)現(xiàn)最大相差小于10℃。由于溫度場(chǎng)數(shù)據(jù)在疲勞計(jì)算中是一個(gè)過(guò)程數(shù)據(jù)，該偏差對(duì)疲勞計(jì)算結(jié)果的影響暫時(shí)還無(wú)法判斷，需要求出壁厚溫度梯度后再作判斷。根據(jù)第3.7節(jié)壁厚溫度梯度結(jié)果可以認(rèn)為，一維非穩(wěn)態(tài)導(dǎo)熱反演數(shù)值計(jì)算方法滿足工程需要。

圖13　反演計(jì)算內(nèi)壁溫度

3.5數(shù)值計(jì)算最小空間離散度分析

壁厚空間離散度對(duì)壁厚溫度梯度（直接影響疲勞計(jì)算結(jié)果）有影響，本文用ANSYS軟件計(jì)算不同離散尺度下壁厚溫度數(shù)據(jù)，進(jìn)而數(shù)值計(jì)算出溫度梯度，發(fā)現(xiàn)平均溫度、線性溫度梯度受離散尺度影響較小，非線性溫度梯度受離散尺度影響較大，分析可知（見(jiàn)圖14，15），試驗(yàn)用的管壁至少9等分可以保證計(jì)算精度。

圖14　ANSYS計(jì)算管壁徑向線性溫度梯度

圖15　ANSYS計(jì)算管壁徑向非線性溫度梯度

3.6一維非穩(wěn)態(tài)導(dǎo)熱反演數(shù)值計(jì)算壁厚溫度場(chǎng)

根據(jù)第3.5節(jié)可知，壁厚方向至少要有10個(gè)節(jié)點(diǎn)的溫度數(shù)據(jù)，才能保證計(jì)算精度。但是根據(jù)第2.1.2（3）節(jié)中邊界點(diǎn)溫度差分方程反演的穩(wěn)定性條件，時(shí)間離散尺度1 s（溫度檢測(cè)裝置間隔1 s采集數(shù)據(jù)）的情況下，空間離散尺度最小3mm。于是管壁9等分后（如圖16所示），點(diǎn)1，2，3，4，5，6，7的溫度可以直接通過(guò)點(diǎn)10的溫度導(dǎo)熱數(shù)值反演計(jì)算得到，點(diǎn)8，9與管道外壁點(diǎn)10間距都小于3mm，不滿足反演方程的穩(wěn)定性條件，無(wú)法直接通過(guò)點(diǎn)10的溫度導(dǎo)熱反演計(jì)算，但是點(diǎn)8，9的溫度既可以通過(guò)其他點(diǎn)反演計(jì)算得到，也可以通過(guò)線性插值簡(jiǎn)化處理，因?yàn)辄c(diǎn)7，8， 9，10溫度十分接近。最終管壁上10個(gè)點(diǎn)溫度全部得到。

圖16　管壁徑向離散示意

3.7壁厚溫度梯度計(jì)算

利用第3.6節(jié)中一維非穩(wěn)態(tài)導(dǎo)熱反演計(jì)算得到的各節(jié)點(diǎn)溫度，根據(jù)式（7）數(shù)值計(jì)算得到管壁徑向平均溫度（見(jiàn)圖17），與ANSYS軟件計(jì)算結(jié)果相比一致（ANSYS軟件用試驗(yàn)監(jiān)測(cè)的管道內(nèi)壁溫度數(shù)據(jù)計(jì)算得到各層溫度，然后求得管壁徑向平均溫度）。

圖17　管壁徑向平均溫度

利用第3.6節(jié)中一維非穩(wěn)態(tài)導(dǎo)熱反演計(jì)算得到的各節(jié)點(diǎn)溫度，根據(jù)式（8）數(shù)值計(jì)算得到管壁徑向線性溫度梯度（見(jiàn)圖18），與ANSYS軟件計(jì)算結(jié)果相差小于5℃（ANSYS軟件用試驗(yàn)監(jiān)測(cè)的管道內(nèi)壁溫度數(shù)據(jù)計(jì)算得到各層溫度，然后求得管壁徑向線性溫度梯度）。

利用第3.6節(jié)中一維非穩(wěn)態(tài)導(dǎo)熱反演計(jì)算得到的各節(jié)點(diǎn)溫度，根據(jù)式（9）數(shù)值計(jì)算得到管壁徑向非線性溫度梯度（見(jiàn)圖19），與ANSYS軟件計(jì)算結(jié)果相差小于2℃（ANSYS軟件用試驗(yàn)監(jiān)測(cè)的管道內(nèi)壁溫度數(shù)據(jù)計(jì)算得到各層溫度，然后求得管壁徑向非線性溫度梯度）。

圖18　管壁徑向線性溫度梯度

圖19　管壁徑向非線性溫度梯度

4　結(jié)語(yǔ)

（1）核電廠疲勞監(jiān)測(cè)系統(tǒng)所涉及的管道，在熱沖擊和熱分層兩種載荷下，可以簡(jiǎn)化用一維數(shù)值模型；一維非穩(wěn)態(tài)數(shù)值反演計(jì)算方法可以得到管壁各節(jié)點(diǎn)溫度值，計(jì)算精度滿足工程要求。

（2）對(duì)于一個(gè)具體的監(jiān)測(cè)點(diǎn)，在進(jìn)行導(dǎo)熱數(shù)值反演計(jì)算前，首先需要分析最小空間離散度，目的是確保疲勞計(jì)算結(jié)果的精度，本文試驗(yàn)中用到的9mm厚管道，經(jīng)分析至少需要10個(gè)節(jié)點(diǎn)。

［1］李崗，梁兵兵，馬志才.穩(wěn)壓器波動(dòng)管溫度監(jiān)測(cè)及壽命評(píng)估［J］.壓力容器，2012，29（6）:63-67.

［2］唐永進(jìn).壓力管道應(yīng)力分析［M］.北京:中國(guó)石化出版社，2003.

［3］石媛媛，王明毓，湯鳳，等.某研發(fā)項(xiàng)目穩(wěn)壓器排放管線應(yīng)力分析［J］.壓力容器，2015，32（7）:47-51.

［4］［美］Suresh S.材料的疲勞［M］.北京:國(guó)防工業(yè)出版社，1993.

［5］署恒木，黃朝琴，李翠偉.瞬態(tài)熱傳導(dǎo)問(wèn)題的加權(quán)最小二乘無(wú)網(wǎng)格法［J］.計(jì)算力學(xué)學(xué)報(bào)，2008，25（6）:904-908.

［6］張征明，王敏稚，何樹(shù)延.HTR-10核安全一級(jí)管道的力學(xué)分析［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2000，40（12）:14-17.

［7］張學(xué)學(xué)，李桂馥.熱工基礎(chǔ)［M］.北京:高等教育出版社，2009.

［8］AFCEN，RCC-M，壓水堆核島機(jī)械設(shè)備設(shè)計(jì)和建造規(guī)則（2000版+2002補(bǔ)遺）［S］.

［9］署恒木.軸對(duì)稱熱傳導(dǎo)有限元格式［J］.計(jì)算力學(xué)學(xué)報(bào)，2000，17（1）:123-126.

A Numerical Calculation of Inverse Heat Conduction for Fatigue M onitoring System of Nuclear Power Plant

LIU Lang，MENG A-jun，LING Jun，HE Da-yu
（China Nuclear Power Design Co.，Ltd.（Shenzhen），Shenzhen 518172，China）

The unsteady inverse heat conduction calculation is very important for the fatigue monitoring system，and is the focus of the thesis.The study is based on a test，the temperature of inner wall and outer wall ismeasured directly.First，the equation of numerical calculation of one dimensional unsteady inverse heat conduction is studied.Second，with the help of the equation，the temperature of innerwall is calculated by themeasured temperature of outer wall.Last，the calculated and themeasured temperature of inner wall is compared.The numerical calculation of one dimensional unsteady inverse heat conduction could be applied to the fatiguemonitoring system.

fatiguemonitoring system；one dimensional unsteady inverse heat conduction；numerical calculation

TH123；TL353

A

1001-4837（2015）11-0020-07

10.3969/j.issn.1001-4837.2015.11.004

2015-09-14

2015-10-27

劉浪（1983-），男，工程師，主要從事工程力學(xué)方面的工作，通信地址:518172廣東省深圳市龍崗區(qū)黃閣北路天安數(shù)碼城5號(hào)樓深圳中廣核工程設(shè)計(jì)有限公司，E-mail:liulang@cgnpc.com.cn