張 鍇

（上海核工程研究設(shè)計(jì)院 上海 200233）

蒸汽發(fā)生器傳熱管間隙對(duì)傳熱管動(dòng)態(tài)特性的影響分析

張 鍇

（上海核工程研究設(shè)計(jì)院 上海 200233）

傳熱管流致振動(dòng)是核電廠蒸汽發(fā)生器傳熱管失效的主要原因之一，在核電廠設(shè)計(jì)蒸汽發(fā)生器時(shí)，需對(duì)蒸汽發(fā)生器傳熱管流致振動(dòng)問題進(jìn)行分析。傳熱管與支撐板及抗振條之間存在小尺度間隙，這類間隙具有非線性效應(yīng)，在進(jìn)行流致振動(dòng)線性分析時(shí)應(yīng)考慮對(duì)間隙進(jìn)行線性化等效處理。本文從理論研究和模擬分析兩方面出發(fā)，對(duì)傳熱管與支撐板及抗振條之間間隙對(duì)傳熱管動(dòng)態(tài)特性的影響進(jìn)行分析。理論和模擬分析可知，傳熱管間隙對(duì)傳熱管整體振動(dòng)的作用接近于簡(jiǎn)支。在進(jìn)行流致振動(dòng)分析時(shí)，可采用簡(jiǎn)支代替間隙進(jìn)行線性分析。

傳熱管，流致振動(dòng)，間隙，簡(jiǎn)支

傳熱管流致振動(dòng)(Flow Induced Vibration，F(xiàn)IV)是核電廠蒸汽發(fā)生器傳熱管失效的主要原因之一。Enrico Fermi核電站試運(yùn)行13天后，蒸汽發(fā)生器有6根傳熱管失效[1]。之后該電廠又持續(xù)運(yùn)行42天，再次發(fā)現(xiàn)另39根傳熱管發(fā)生泄漏。瑞典Ringhals 3號(hào)反應(yīng)堆曾發(fā)現(xiàn)數(shù)十根蒸汽發(fā)生器傳熱管存在磨損，導(dǎo)致傳熱管厚度減至原有厚度的10%。這些事故均由蒸汽發(fā)生器二次側(cè)流動(dòng)誘發(fā)傳熱管振動(dòng)所致。由此可見，蒸汽發(fā)生器傳熱管流致振動(dòng)及由此引起的傳熱管磨損不僅會(huì)增加電廠維修工作量，造成經(jīng)濟(jì)損失，而且可能影響反應(yīng)堆安全和核電廠的可靠運(yùn)行。因此在進(jìn)行核電廠蒸汽發(fā)生器設(shè)計(jì)時(shí)，有必要對(duì)蒸汽發(fā)生器傳熱管流致振動(dòng)問題進(jìn)行研究和分析，以避免或減少由此帶來的危害。

在進(jìn)行傳熱管流致振動(dòng)分析時(shí)，注意到傳熱管與支撐板及抗振條之間存在間隙。間隙的尺度相對(duì)于整體結(jié)構(gòu)的尺寸非常小。當(dāng)傳熱管發(fā)生流致振動(dòng)時(shí)，傳熱管與支撐板及抗振條之間會(huì)反復(fù)碰撞。在這一從非接觸到接觸再到非接觸的過程中，間隙本身的剛度會(huì)在極大與極小之間反復(fù)突變。間隙的這種性質(zhì)使傳熱管的振動(dòng)出現(xiàn)顯著的非線性現(xiàn)象，由此導(dǎo)致傳統(tǒng)的線性分析方法不再適用于傳熱管的流致振動(dòng)。而采用完全非線性方法分析傳熱管的流致振動(dòng)，對(duì)計(jì)算資源和計(jì)算時(shí)間的耗費(fèi)又非常大。因此有必要在分析時(shí)對(duì)間隙進(jìn)行線性化等效處理，在簡(jiǎn)化分析過程的同時(shí)保留傳熱管的主要?jiǎng)討B(tài)特性。本文對(duì)傳熱管與支撐板及抗振條之間間隙對(duì)傳熱管動(dòng)態(tài)特性的影響進(jìn)行分析，由此確定在傳熱管流致振動(dòng)線性分析中傳熱管間隙的等效處理方法。經(jīng)分析可知，支撐板與抗振條對(duì)傳熱管可視作簡(jiǎn)支處理，以此為基礎(chǔ)進(jìn)行傳熱管線性振動(dòng)分析，可大幅簡(jiǎn)化分析過程。

1 傳熱管間隙的結(jié)構(gòu)

分析前首先給出核電廠蒸汽發(fā)生器傳熱管間隙的實(shí)際結(jié)構(gòu)。傳熱管間隙共有兩種，分別為傳熱管與支撐板(TSP)間的梅花孔間隙以及傳熱管與抗振條(AVB)之間的條狀間隙。圖1為傳熱管與支撐板間梅花孔間隙(左)及傳熱管與抗振條間條狀間隙(右)的示意圖，圖2為蒸汽發(fā)生器傳熱管的整體結(jié)構(gòu)。

2 理論分析

2.1間隙水膜的作用

傳熱管與支撐板間的間隙中存在水膜。當(dāng)傳熱管發(fā)生流致振動(dòng)時(shí)，水膜對(duì)傳熱管產(chǎn)生一定作用力。傳熱管與抗振條間的間隙中不存在水膜，可不考慮水膜的作用。根據(jù)傳熱管間隙的實(shí)際結(jié)構(gòu)及蒸汽發(fā)生器二次側(cè)的流場(chǎng)條件，傳熱管流致振動(dòng)中的間隙水膜問題可視作狹小間隙中的流體粘性流動(dòng)，即當(dāng)傳熱管沒有與支撐板碰撞時(shí)，可用經(jīng)典的雷諾方程(式1)對(duì)傳熱管在振動(dòng)過程中的受力情況進(jìn)行分析。

圖1 傳熱管與支撐板及抗振條間的間隙Fig.1 Gap between tube and TSP, AVB.

圖2 蒸汽發(fā)生器傳熱管的整體結(jié)構(gòu)Fig.2 Configuration of steam generator tubes.

式中，p為流體相對(duì)壓力(即相對(duì)靜壓的壓力變化值)，h為間隙厚度，η為間隙中水的動(dòng)力粘性系數(shù)。

簡(jiǎn)化式1，并考慮雷諾方程的邊界條件，得到定解問題如下：

式中，a和b分別為間隙中心到其長(zhǎng)度和寬度方向邊界的距離。根據(jù)文獻(xiàn)《Roark’s Formulas for Stress and Strain》[2]中的基本接觸公式，間隙實(shí)際寬度為：

其中，p為單位長(zhǎng)度接觸力，D為傳熱管外徑，EC為接觸柔度。間隙實(shí)際長(zhǎng)度等于支撐板厚度Kt，即：

求解定解問題(2)式，得到方程的解為：

其中，p0為通解，p*為特解。

進(jìn)一步對(duì)間隙水膜壓力在接觸面積上積分，得到未碰撞時(shí)間隙水膜對(duì)傳熱管的作用力為：

由式6可知，水膜對(duì)傳熱管的作用力與間隙厚度的變化速度，即傳熱管本身的振動(dòng)速度有關(guān)。這說明水膜作用力是一種阻尼力，水膜的存在導(dǎo)致傳熱管振動(dòng)的阻尼增大。

2.2傳熱管間隙振動(dòng)分析

考慮三個(gè)間隙之間兩段傳熱管的振動(dòng)，將兩端的間隙視作簡(jiǎn)支。對(duì)于中間的間隙，僅保留其對(duì)傳熱管的作用力，從而分析間隙作用力對(duì)傳熱管的影響。當(dāng)傳熱管流致振動(dòng)發(fā)生時(shí)，傳熱管分別受到流致振動(dòng)作用力Fq、間隙水膜作用力Fh及間隙的碰撞力Fk。將三個(gè)間隙之間兩段傳熱管的振動(dòng)問題等效為單自由度阻尼彈簧系統(tǒng)(如圖3)，其中，m為傳熱管第一階模態(tài)質(zhì)量，K為系統(tǒng)剛度，c為系統(tǒng)阻尼。

圖3 三個(gè)間隙間傳熱管振動(dòng)的單自由度阻尼彈簧系統(tǒng)等效圖Fig.3 Single DOF system with damp and spring equivalent to tube vibration among three support plates.

根據(jù)等效單自由度阻尼彈簧系統(tǒng)，列出平衡方程(式7)和各作用力的表達(dá)式(式8–10)：

其中，q為流致振動(dòng)作用力幅值，ω為流致振動(dòng)作用力頻率，l為兩塊支撐板或抗振條間距，σ為支撐板或抗振條厚度，D和d為傳熱管外徑及內(nèi)徑，z0為初始間隙厚度，K(z)為系統(tǒng)剛度(包括簡(jiǎn)支梁折算剛度及傳熱管的碰撞剛度，忽略間隙水膜的剛度)：

其中，E為材料彈性模量，Ks為碰撞剛度。

將式8–10代入式7，進(jìn)一步整理得到：

方程(12)是一個(gè)強(qiáng)非線性方程。將核電廠蒸汽發(fā)生器傳熱管的各物理量代入方程，并采用有限差分法求解該方程。引入中心差分格式改寫方程(12)，并將方程改寫為時(shí)間推進(jìn)形式：

根據(jù)式(13)編寫計(jì)算程序求解方程，得到振動(dòng)方程(12)的解。圖4(a)為流致振動(dòng)激勵(lì)下，傳熱管與支撐板或抗振條接觸位置的振動(dòng)情況，其中實(shí)線表示考慮水膜阻尼效應(yīng)的情況，虛線表示不考慮水膜阻尼效應(yīng)的情況。由圖可知，由于間隙的存在，傳熱管在間隙內(nèi)來回碰撞。若不考慮阻尼效應(yīng)，系統(tǒng)振動(dòng)十分復(fù)雜，表現(xiàn)出明顯的非線性振動(dòng)特征。圖4(b)為傳熱管間隙對(duì)傳熱管流致振動(dòng)的影響，其中實(shí)線表示存在間隙的情況，虛線表示沒有間隙的情況(即傳熱管與支撐板或抗振條緊密貼合)。由圖可知，間隙的存在使傳熱管在間隙中來回碰撞。由于間隙很小，間隙中傳熱管的振動(dòng)幅值受限，但間隙中傳熱管的轉(zhuǎn)動(dòng)卻未受影響。此時(shí)間隙的作用與簡(jiǎn)支類似，因此線性分析時(shí)，傳熱管間隙對(duì)傳熱管整體可視作簡(jiǎn)支處理。圖4(c)為不同的流致振動(dòng)激勵(lì)強(qiáng)弱對(duì)傳熱管振動(dòng)的影響，其中實(shí)線表示大強(qiáng)度激勵(lì)，虛線表示小強(qiáng)度激勵(lì)。由圖可見，激勵(lì)較小時(shí)，傳熱管經(jīng)若干次碰撞后很快衰減到不再碰撞的狀態(tài)。此時(shí)，傳熱管間隙的作用逐步由簡(jiǎn)支過渡到自由。激勵(lì)較大時(shí)，這種過渡則比較緩慢。

總之，流致振動(dòng)條件下，傳熱管間隙的存在使傳熱管與支撐板或抗振條發(fā)生反復(fù)碰撞，由此導(dǎo)致的微動(dòng)磨損是傳熱管流致振動(dòng)失效的主要原因。傳熱管間隙對(duì)傳熱管振動(dòng)的影響主要有以下幾點(diǎn)：

(1) 傳熱管間隙內(nèi)水膜是傳熱管振動(dòng)的阻尼，可有效減弱傳熱管的碰撞次數(shù)和頻率；(2) 傳熱管間隙對(duì)傳熱管整體振動(dòng)的作用接近于簡(jiǎn)支；(3) 傳熱管振動(dòng)經(jīng)一段時(shí)間衰減后，不再發(fā)生碰撞。此時(shí)，間隙對(duì)傳熱管振動(dòng)的影響由簡(jiǎn)支過渡到自由。

圖4 水膜阻尼(a)、傳熱管間隙(b)及流致振動(dòng)激勵(lì)(c)對(duì)傳熱管振動(dòng)的影響Fig.4 Influence of water film damp(a), tube gap(b), FIV excitation(c) on tube vibration.

3 模擬計(jì)算

3.1分析模型

采用有限元瞬態(tài)分析方法研究傳熱管間隙對(duì)傳熱管動(dòng)態(tài)特性的影響。建立三種單根傳熱管分析模型，以間隙單元模擬傳熱管與支撐板或抗振條間的間隙。三種模型中，分別對(duì)兩個(gè)間隙之間一段傳熱管的中點(diǎn)施加激勵(lì)，使其自由振動(dòng)。通過考察傳熱管激勵(lì)點(diǎn)的自由振動(dòng)狀態(tài)，分析傳熱管間隙對(duì)傳熱管動(dòng)態(tài)特性的影響。在分析間隙對(duì)單根傳熱管影響的同時(shí)，在三個(gè)模型中用簡(jiǎn)支替換間隙進(jìn)行再次分析，從而比對(duì)間隙與簡(jiǎn)支對(duì)傳熱管的影響，三種分析模型見圖5。模型1 底部固支、頂部間隙、管中點(diǎn)加載(模擬管板上方傳熱管)；模型2 兩端簡(jiǎn)支、中點(diǎn)間隙、1/4管長(zhǎng)點(diǎn)加載(模擬三塊支撐板間傳熱管)；模型3 兩端簡(jiǎn)支、三分點(diǎn)間隙、彎管中點(diǎn)加載(模擬與抗振條接觸的最外側(cè)傳熱管)。

圖5 模型1(a)，模型2(b)，模型3(c)Fig.5 Model 1(a), Model 2(b), Model 3(c).

3.2結(jié)果分析

考察間隙中的傳熱管振動(dòng)。當(dāng)傳熱管振動(dòng)尚未衰減時(shí)，傳熱管將在間隙中來回碰撞，兩側(cè)支撐板或抗振條將交替提供作用力使傳熱管回彈，如圖6(a)所示，這與理論分析的結(jié)果十分接近。

圖6 間隙中的傳熱管振動(dòng)(a)及模型1間隙中傳熱管振動(dòng)狀態(tài)的演變(b)Fig.6 Tube vibration in the gap(a) and tube vibration in the gap of Model 1(b).

分別考察三個(gè)分析模型的計(jì)算結(jié)果。對(duì)于模型1，當(dāng)傳熱管振動(dòng)尚未衰減時(shí)，傳熱管在間隙中來回碰撞。由圖6(b)看到，傳熱管在來回碰撞的過程中會(huì)多次碰撞同一側(cè)，表明該振動(dòng)的非線性特征。由于間隙很小，間隙中傳熱管的振動(dòng)幅值受限，但間隙中傳熱管的轉(zhuǎn)動(dòng)卻沒有受到影響。此時(shí)整體上傳熱管呈現(xiàn)底部固支頂部簡(jiǎn)支的動(dòng)態(tài)特性。經(jīng)過一段時(shí)間的衰減，間隙中的傳熱管振幅逐步減小到不再發(fā)生碰撞。此時(shí)，間隙的影響消失，傳熱管振動(dòng)頻率下降，逐步演化到頂部自由振動(dòng)的狀態(tài)。

模型2計(jì)算結(jié)果與模型1相似。當(dāng)傳熱管振動(dòng)尚未衰減時(shí)，間隙中的傳熱管在間隙中來回碰撞。此時(shí)整體上傳熱管呈現(xiàn)兩端簡(jiǎn)支中點(diǎn)也簡(jiǎn)支的動(dòng)態(tài)特性。經(jīng)過一段時(shí)間的衰減，間隙中的傳熱管振幅逐步減小到不再發(fā)生碰撞。傳熱管振動(dòng)頻率下降，逐步演化到中點(diǎn)沒有間隙的自由狀態(tài)下振動(dòng)。模型2中點(diǎn)存在間隙和簡(jiǎn)支兩種情況，傳熱管加載點(diǎn)振動(dòng)狀態(tài)比較及頻率見圖7。由圖可知，兩種情況下的振動(dòng)頻率基本相同，幅值比較接近，間隙可視作簡(jiǎn)支處理。

圖7 模型2傳熱管加載點(diǎn)振動(dòng)狀態(tài)比較(a)及其傅立葉變換(b)Fig.7 Tube load point vibration of Model 2(a) and FFT results(b).

模型3振動(dòng)狀態(tài)更為復(fù)雜。傳熱管中點(diǎn)的振動(dòng)雖然與簡(jiǎn)支條件下的振動(dòng)不完全一致，但兩者的振動(dòng)頻率仍然相同，幅值也比較接近(如圖8)，因此也可將間隙作為簡(jiǎn)支處理。

總之，通過分別模擬管板上方到第一塊支撐板、三塊支撐板之間以及U型彎管區(qū)域多條抗振條之間傳熱管的振動(dòng)狀態(tài)，分析了傳熱管間隙對(duì)傳熱管動(dòng)態(tài)特性的影響。經(jīng)過分析可知，較小的傳熱管間隙通?？勺鳛楹?jiǎn)支考慮。當(dāng)傳熱管振動(dòng)尚未衰減時(shí)，間隙中的傳熱管在間隙中來回碰撞，表現(xiàn)出一定的非線性特征。經(jīng)一段時(shí)間衰減后，間隙中的傳熱管振幅逐步減小到不再發(fā)生碰撞。此時(shí)，間隙的作用消失，傳熱管振動(dòng)頻率下降，間隙中的傳熱管逐步演化到自由振動(dòng)的狀態(tài)。

圖8 模型3傳熱管中點(diǎn)的振動(dòng)狀態(tài)比較(a)及其傅立葉變換(b)Fig.8 Tube apex vibration of Model 3(a) and FFT results(b).

4 結(jié)語

傳熱管流致振動(dòng)是核電廠蒸汽發(fā)生器傳熱管失效的主要原因之一，一般由蒸汽發(fā)生器二次側(cè)流動(dòng)誘發(fā)傳熱管振動(dòng)及由此引起的傳熱管磨損導(dǎo)致。為避免傳熱管流致振動(dòng)導(dǎo)致傳熱管失效，在進(jìn)行核電廠蒸汽發(fā)生器設(shè)計(jì)時(shí)，有必要對(duì)蒸汽發(fā)生器傳熱管流致振動(dòng)問題進(jìn)行研究和分析。注意到傳熱管與支撐板及抗振條之間存在小尺度間隙。當(dāng)流致振動(dòng)發(fā)生時(shí)，間隙中的傳熱管與支撐板及抗振條之間會(huì)反復(fù)碰撞，這意味著傳統(tǒng)的線性分析方法不適用于傳熱管流致振動(dòng)分析。而采用完全非線性方法又費(fèi)時(shí)費(fèi)力。因此，可考慮對(duì)間隙進(jìn)行線性化等效處理，以簡(jiǎn)化分析過程并保留傳熱管的主要?jiǎng)討B(tài)特性。

本文從理論研究和模擬分析兩方面出發(fā)，對(duì)傳熱管與支撐板及抗振條之間間隙對(duì)傳熱管動(dòng)態(tài)特性的影響進(jìn)行分析，以找出傳熱管流致振動(dòng)線性分析中傳熱管間隙的等效處理方法。經(jīng)理論分析可知，傳熱管間隙對(duì)傳熱管整體振動(dòng)的作用接近于簡(jiǎn)支，而傳熱管間隙內(nèi)水膜對(duì)傳熱管振動(dòng)起到阻尼效應(yīng)。當(dāng)傳熱管振動(dòng)經(jīng)一段時(shí)間的衰減后，傳熱管間隙的作用將逐步消失。模擬計(jì)算可得，較小的傳熱管間隙可作為簡(jiǎn)支處理。當(dāng)傳熱管振動(dòng)尚未衰減時(shí)，傳熱管在間隙中反復(fù)發(fā)生碰撞，并表現(xiàn)出一定的非線性特征。當(dāng)振動(dòng)經(jīng)一段時(shí)間衰減后，間隙的作用消失，傳熱管振動(dòng)頻率下降，間隙中的傳熱管逐步演化到自由振動(dòng)的狀態(tài)?？傊?，傳熱管間隙在局部上導(dǎo)致傳熱管振動(dòng)出現(xiàn)明顯的非線性效應(yīng)，但在整體上仍可視作簡(jiǎn)支處理。這意味著可采用線性方法分析傳熱管流致振動(dòng)，從而大幅簡(jiǎn)化流致振動(dòng)分析。

致謝本文在研究和撰寫過程中得到曲大莊和姚偉達(dá)老師的悉心指導(dǎo)以及張可豐和金樂的有力幫助，在此表示衷心感謝！

1 Chen S S, 馮振寧, 張希農(nóng). 圓柱結(jié)構(gòu)的流動(dòng)誘發(fā)振動(dòng)[M]. 北京: 石油工業(yè)出版社, 1993: 6–8 Chen S S, FENG Zhenning, ZHANG Xinong. Flow-induced vibration of circular cylindrical structures[M]. Beijing: Petroleum Industry Press, 1993: 6–8

2 Young W C, Budynas R G. Roark’s formulas for stress and strain[M]. 7thed. New York: McGraw-Hill Press, 2001: 702–705

Influence of steam generator tube gaps on the dynamic characteristics of tubes

ZHANG Kai
(Shanghai Nuclear Engineering Research and Design Institute, Shanghai 200233, China)

Background: Flow induced vibration (FIV) is one of the main cause of steam generator tubes failure. Flow induced vibration of steam generator tubes shall be analyzed with the design of steam generator. Small-scale gaps exist between tubes and tube support plates (TSP) or anti-vibration bars (AVB). These gaps make the system non-linear. Purpose: Linear equivalence to the gaps shall be performed in the linear FIV analysis. Methods: This paper analyzes the influence of tube gaps on the dynamic characteristics of tubes with theoretical and simulation methods. Results: According to the analyses, the influence of gaps on tubes is similar to that of simple supports. Conclusions: Simple supports can be used to substitute gaps in the linear analysis for flow induced vibration.

Tubes, Flow induced vibration, Gap, Simple support

TL35

10.11889/j.0253-3219.2013.hjs.36.040617

張鍇，男，1984年出生，2009年于清華大學(xué)獲碩士學(xué)位，工程師，專業(yè)：反應(yīng)堆結(jié)構(gòu)力學(xué)

2012-10-31，

2013-02-27

CLC TL35