吊籃振動特性分析研究

黃旋， 吳萬軍， 齊歡歡， 馮志鵬， 葉獻輝

（中國核動力研究設計院 核反應堆系統設計技術重點實驗室，成都610041）

0 引言

1 結構簡介

吊籃組件由吊籃法蘭、吊籃筒體、出水口管嘴、堆芯下板、堆芯支撐柱、堆芯支撐板和二次支撐構件組成。吊籃組件通過吊籃法蘭固定，吊籃法蘭在預緊螺栓和壓緊環的作用下，壓緊固定在壓力容器支撐臺階上，吊籃下端自由。吊籃法蘭、堆芯支撐板與吊籃筒體通過焊接連接；堆芯下板與吊籃筒體通過螺栓連接，堆芯支撐柱通過螺栓與堆芯下板和堆芯支撐板連接。

2 分析模型

2.1 模型簡化

為能準確模擬吊籃組件固有振動特性，又能使模型盡量簡化，忽略堆芯支撐板上的流水孔，簡化為實心圓板；忽略吊籃出水口的局部倒角；模型中詳細模擬了吊籃法蘭、吊籃筒體、圍板-成型板組件等部件，上述部件是吊籃組件主體，代表了其主要的剛性和質量特征。其余部件只考慮質量的影響，以附件質量的形式附加到結構上。吊籃組件簡化模型如圖1所示。

圖1 吊籃簡化模型

2.2 網格劃分

本文對吊籃及其主要附屬結構在空氣和安裝于反應堆壓力容器內水下環境的情況分別進行模擬。

2.2.1 空氣中

為降低求解規模并提高求解精度，對影響結構剛度的主要部件采用映射網格劃分，包括吊籃法蘭、吊籃筒體和圍板；對主要體現為質量的次要部件采用自由網格劃分，包括出水口管嘴、成型板和堆芯支撐板。為了準確模擬吊籃法蘭在壓緊彈簧的作用下，壓緊固定在壓力容器支撐臺階上的邊界條件，吊籃法蘭與壓力容器支撐臺接觸的部分用實體單元SOLID185模擬，其余部分均用殼單元SHELL63模擬，空氣中吊籃網格模型見圖2。

2.2.2 靜水中

2.1四組患者心臟結構相關指標比較:OSAHS患者組中,中度組與重度組患者的LVDD、LVDS、LA、IVS、LVPWD指標與健康對照組比較明顯升高,差意有統計學意義(P<0.05),見表1。[3]

吊籃及其主要附屬結構安裝于反應堆壓力容器內水下環境中，其結構模型與空氣中相同，但是需要考慮吊籃內外水對其振動特性的影響?？紤]吊籃處于反應堆壓力容器內部水下環境時，使用聲流體Fluid30單元模擬水環境。由于使用殼單元模擬吊籃筒體，只能與吊籃單側的水體進行流固耦合。因此，吊籃外側水采用Fluid30單元模擬，內側水通過附加質量的方式折算到結構上。靜水中的吊籃整體網格模型如圖2所示。

圖2 吊籃在空氣和靜水中的網格模型

3 彈簧等效剛度計算

壓緊彈簧位于壓緊部件法蘭與吊籃法蘭之間，將吊籃組件壓緊并定位于反應堆壓力容器內。在對吊籃組件邊界條件進行模擬時，將壓緊彈簧簡化為垂向線性彈簧，彈簧剛度通過靜力計算確定。

3.1 壓緊彈簧結構

壓緊彈簧是一個鍛環，位于上支承法蘭和吊籃法蘭之間，當壓力容器頂蓋蓋緊時，頂蓋與上支承法蘭上表面貼合，從而壓縮彈簧。通過壓緊彈簧的作用，可阻止上部堆內和下部堆內構件因縱向載荷作用引起的豎直方向的移動。

3.2 壓緊彈簧模型

計算采用二維軸對稱模型。壓緊彈簧在兩個接觸面處均呈圓弧形，它與上部組件法蘭及吊籃法蘭之間的接觸圓在徑向的位置會隨壓緊力的變化而有微小的移動，應作為接觸問題來處理。在計算時，采用8節點平面軸對稱單元進行離散，在接觸面上采用面對面接觸單元，并使初始狀態為壓緊彈簧與上部組件法蘭及吊籃法蘭剛好接觸。上部組件法蘭與吊籃法蘭均采用剛性接觸面來處理。模型的有限元網格劃分見圖3。

圖3 壓緊彈簧網格模型

3.3 計算方法

采用有限元分析法進行計算，在計算剛度時，給定接觸點的位移，計算出接觸點的反力，從而求出垂向剛度。

3.4 載荷及邊界條件

模型的邊界條件見圖3。計算剛度時，壓緊彈簧的下部，吊籃法蘭接觸面的約束條件為Uy=0 mm，上部組件法蘭接觸面施加逐步增大的位移載荷：Uy=-1 mm,-2 mm,-3 mm,-4 mm。

3.5 壓緊彈簧剛度計算結果

表1 壓緊彈簧在不同載荷下的垂向剛度值

4 邊界條件

吊籃放置在反應堆壓力容器內，其下端自由，上端通過吊籃法蘭下緣掛在壓力容器支撐臺上，上緣通過壓緊彈簧壓緊。建模時將壓緊彈簧簡化為線彈簧單元；約束法蘭下緣外側節點的Z向位移；約束吊籃法蘭X=0處節點X向位移，約束Y=0處節點Y向位移。在彈性變形假設下，為模擬吊籃法蘭與壓力容器支撐臺階之間的摩擦力[15]，在法蘭下緣XY平面內0°、45°、90°、135°、180°、225°、270°、315°方向上建立8個徑向拉壓彈簧和8個切向彈簧，單個彈簧剛度為壓緊環等效總剛度的1/8。

5 分析方法

報告采用模態分析方法進行振動特性分析，計算吊籃組件的梁式模態和殼式模態?？紤]反應堆為線性n個自由度的結構系統，其無阻尼系統的自由振動方程為

其中：M為n×n維質量矩陣；K為 n×n維剛度矩陣；u¨為n維節點加速度矢量；u為n維節點位移矢量；0為 n維零元素矢量。

令u=Φsin（ωt+α），Φ為n維振幅矢量，代入式（1）得：

由線性代數理論可知，若使式（2）有非零解，必須有：

求解式（3）特征方程，得出n個特征值ωi（i=1,2,...,n），則系統的固有頻率為

將ωi代入式（2），可求解出相應的振幅矢量Φi，稱為第i階固有振型矢量。

固有振型矢量具有正交性，即當i≠j時：

而當i=j時：

為進行后續響應譜分析，將Φi進行正則化處理，使

得：

式中，Φ′i稱為正則振型。

6 計算結果

6.1 空氣中

吊籃在空氣中振動特性的分析結果見表2，梁式模態和殼式模態如圖4~圖6所示。6.2 靜水中

圖4 吊籃在空氣中的一階梁式振型

圖5 吊籃在空氣中的殼式振型（1,2）

圖6 吊籃在空氣中的殼式振型（1,3）

吊籃在靜水中振動特性的分析結果見表3，梁式模態和殼式模態如圖7~圖9所示。

7 結論

圖7 吊籃在靜水中的一階梁式振型

圖8 吊籃在靜水中的殼式振型（1,2）

圖9 吊籃在靜水中的殼式振型（1,3）

表3 吊籃在靜水中的固有頻率 Hz

本文建立吊籃及其主要附屬結構在空氣和安裝于反應堆壓力容器內水下環境中的三維有限元模型，并分別進行了振動特性分析。將動態特性分析結果與試驗結果進行比較可知，通過有限元模型計算得到的吊籃組件在空氣和靜水中的前三階頻率與相應試驗測量結果的偏差均在5之內，有限元模型可以準確模擬吊籃組件固有振動特性。因此，靜水中的吊籃組件有限元模型可用于后續的吊籃組件流致振動響應分析。