王 墨，馬 駿

（大連理工大學(xué)船舶工程學(xué)院，遼寧大連 116024）

0 引言

影響船舶上層建筑首階固有頻率的因素有很多，如上層建筑自身總體剛度、上層建筑根部支撐的剛度、上層建筑重量大小及重心位置等。人們試圖分析和總結(jié)各種因素對上層建筑固有頻率影響的大小，在此過程中逐漸引入了靈敏度分析的方法。靈敏度分析是研究模型的輸入和輸出間的關(guān)系的一種分析方法，用于定性或定量地研究不同的模型輸入對模型的輸出的影響以及模型對輸入的依賴程度。過去在對上層建筑首階固有頻率進行估算時根據(jù)上層建筑的主要尺度使用基于統(tǒng)計數(shù)據(jù)和回歸分析的經(jīng)驗公式估算[1-3]。但是這些經(jīng)驗公式提出時間較早，適用船型均為排水量較小的船型，沒有完全體現(xiàn)出上層建筑整體剛性和質(zhì)量分布影響，對現(xiàn)代超大型船舶的計算精確度仍需探討[4]。隨著有限元計算方法的完善和計算機計算能力的突飛猛進，人們逐漸利用有限元計算軟件代替經(jīng)驗公式對大型結(jié)構(gòu)進行分析計算。本文利用Ansys有限元計算軟件對某 300,000 DWT級原油船上層建筑進行建模計算，討論船舶上層建筑的重量和中心對上層建筑首階橫向和縱向固有頻率影響的靈敏度。

1 分析模型

本文對一個300,000噸級的原油船建模并進行分析計算。建立的模型包括整個上層建筑結(jié)構(gòu)、上層建筑后方煙囪結(jié)構(gòu)以及兩者下方機艙雙層底以上的全部結(jié)構(gòu)。

模型采用直角坐標系，其中：

原點取為基線的0站處；

X軸沿船長方向且指向艏為正；

Y軸沿船寬方向且指向左舷為正；

Z軸沿垂向且向上為正。

圖1 有限元計算模型

模型中尺寸較大的橫、縱向構(gòu)件及所有板材以SHELL63單元建立，尺寸較小的梁及骨材等加強結(jié)構(gòu)以BEAM188單元建立。上述兩種單元具有拉壓、扭轉(zhuǎn)和彎曲能力。每種單元的每個節(jié)點均有六個自由度，包括沿x、y、z方向的線位移和繞x、y、z軸的轉(zhuǎn)角。機艙和上層建筑內(nèi)的設(shè)備以MASS21單元建立。

模型的材料選為屈服模量為楊氏模量為211GPa、泊松比為0.3、默認密度為7850kg/m3的普通鋼。

該船的肋距為0.85m，故為方便劃分網(wǎng)格并兼顧計算精度，定義網(wǎng)格大小也為0.85m。在計算上層建筑首階固有頻率時，對模型主船體的前、后截斷面及雙層底位置進行簡支約束。模型的主要數(shù)據(jù)見表1和表2。

表1 各層甲板主要數(shù)據(jù)

表2 模型的其他主要數(shù)據(jù)

2 上層建筑重量對固有頻率的影響

在對船舶上層建筑建模時，對線纜、舾裝、空冷管材等物件的重量的處理一般為均布在所在處的甲板上。通過修改各層甲板的密度來調(diào)節(jié)重量和重心，可以避免改變上層建筑的結(jié)構(gòu)剛度。

2.1 方案設(shè)計

討論上層建筑重量對首階固有頻率的影響，一般通過調(diào)節(jié)上層建筑整體結(jié)構(gòu)的密度的方法達到調(diào)節(jié)上層建筑的重量的目的。用這種方法可以僅改變重量而不影響重心和結(jié)構(gòu)的剛度。以整體相當(dāng)密度為10338kg/m3的方案作為標準對比方案(即原船數(shù)據(jù))，另設(shè)計6組相當(dāng)密度與標準方案相當(dāng)密度不同的方案（3組小于相當(dāng)密度、3組大于相當(dāng)密度）。各方案詳情見表3。

2.2 計算結(jié)果

根據(jù)上述方案對模型進行修改并計算，得到的計算結(jié)果匯總見表4。

從表4中的結(jié)果可以看出，隨著上層建筑整體重量的線性增大，上層建筑首階橫向和縱向固有頻率逐漸減小。

表3 各方案詳情

以方案4作為標準方案時，其首階橫向頻率7.87Hz與首階縱向頻率7.40Hz作為對比的標準頻率。其余方案計算得到的結(jié)果均同方案4中的結(jié)果對比。為了能更直觀簡潔地觀察，將表4中計算結(jié)果匯總為圖2和圖3。

為了能更清楚地對比首階固有頻率隨上層建筑總重的變化情況，把得到的計算結(jié)果同重量變化情況繪制在一張圖中，見圖2。其中：橫向、縱向固有頻率按主坐標軸（垂直）繪制，單位Hz；上層建筑總重的變化情況按次坐標軸（垂直）繪制，單位kg。圖2中橫軸為方案序號。

從圖2中可以看出，當(dāng)上層建筑總重呈線性增加時，其首階橫向、縱向固有頻率均按接近線性的方式逐漸減小，固有頻率變化的斜率約為重量變化斜率的一半。把上層建筑固有頻率變化情況以百分比表示，繪制成圖3。

表4 各方案結(jié)果

圖2 上層建筑固有頻率隨重量變化情況

圖3 上層建筑固有頻率隨重量變化情況（%）

圖3中縱軸為變化量，以百分數(shù)表示，橫軸為方案號。從這張圖中可以較直觀地看出上層建筑首階橫向、縱向固有頻率隨重量變化的情況。重量是影響上層建筑的主要因素之一，當(dāng)重量改變較大時，其固有頻率也有較大的變化。

3 上層建筑重心位置對固有頻率的影響

除了重量會影響上層建筑的固有頻率，重心位置的不同也會上層建筑的固有頻率有所影響。對比重心位置不同而對固有頻率產(chǎn)生的影響，要保證上層建筑的總重和結(jié)構(gòu)剛度不變。本節(jié)中使用的方法為調(diào)節(jié)上層建筑各層甲板的密度，使總重保持不變，但重心位置的高度有所變化。

3.1 方案設(shè)計

通過調(diào)節(jié)上層建筑中各層甲板的密度來調(diào)節(jié)上層建筑整體重心位置。為此共設(shè)計了6種方案，并與方案1作對比。各方案的詳情見表5。

表5 各方案詳情

3.2 計算結(jié)果

根據(jù)以上各方案對模型進行修改并計算，得到的結(jié)果匯總到表6。

表6 各方案結(jié)果

將表6中的結(jié)果繪制成曲線圖，見圖4。

圖4 上層建筑固有頻率隨重心變化情況

圖4中橫向、縱向頻率曲線對應(yīng)主坐標軸，單位為Hz；重心高度曲線對應(yīng)次坐標軸，單位為m；橫軸為方案號。

為了更清楚地對比首階頻率隨重心高度位置變化的情況，把頻率改變量和重心高度改變量寫為百分比的形式，見圖5。

圖5 上層建筑固有頻率隨重心變化情況

圖5中圖中重心高度改變曲線對應(yīng)主坐標軸，單位為m；橫向、縱向頻率改變曲線對應(yīng)次坐標軸，單位為Hz；橫軸為方案號。

從以上兩圖中可以看出，隨著上層建筑整體重心的提高，上層建筑的橫向和縱向首階固有頻率均有下降的趨勢。橫向首階固有頻率相對縱向首階固有頻率改變更大。兩固有頻率的改變量同上層建筑重心的改變量大致在同一個數(shù)量級。

當(dāng)船舶的上層建筑較高時，可以簡化地視之為一個懸臂梁。如果懸臂梁的剛度和重量不變，其重心位置越遠離支點端，其固有頻率越低。

4 結(jié)論

改變船舶上層建筑的重量和重心均可有效改變其首階固有頻率的大小。增大上層建筑的重量、提高上層建筑的重心均可降低其首階固有頻率，反之減小上層建筑的重量、降低上層建筑的重心可以增大其首階固有頻率。相對而言，上層建筑的固有頻率對重心高度的變化更為敏感。

修改上層建筑的重心較簡單，可以修改總布置挪動一些不重要的機器設(shè)備。改變上層建筑重量就相對麻煩一些，在不改變結(jié)構(gòu)的前提下可以修改舾裝布置，但是對總重影響不大。當(dāng)通過其他方法對上層建筑重量改變較大時對整船影響也較大。所以改變上層建筑的重心是一個簡單有效又可行的方法。

[1]曹迪, 趙德有, 馬廣宗.船舶上層建筑縱向振動固有頻率估算[J].大連理工大學(xué)學(xué)報, 1989年第4期: 463-469.

[2]広渡智雪, 松本亙平, 加道博章.上部構(gòu)造の前後振動に関する研究（続）[J].日本造船學(xué)會論文集, 1979, (125號): 147-155.

[3]大沼覺, 山本鷹司, 中野元博等.船尾船橋樓の船樓の振動につぃ[J].西部造船會會報, 1969, (38號): 93-120.

[4]趙德有, 林哲.船舶上層建筑整體縱向固有頻率算法研究[J].中國造船.2001(03): 22-26.