推力軸承基座帶預(yù)應(yīng)力諧響應(yīng)分析

丁德勇 張 鈞 賀雙元

1 大連船舶重工有限責(zé)任公司軍事代表室，遼寧大連 116021 2 中國(guó)艦船研究設(shè)計(jì)中心，湖北 武漢 430064

1 引 言

靜力分析也許能確定一個(gè)結(jié)構(gòu)可承受穩(wěn)定載荷的條件，但這還遠(yuǎn)遠(yuǎn)不夠，尤其在載荷隨時(shí)間變化時(shí)更是如此。工程中存在較多的受迫振動(dòng)，即振動(dòng)系統(tǒng)在外界干擾力或干擾位移作用下產(chǎn)生的振動(dòng)。由于干擾力形式不同，可將受迫振動(dòng)分為簡(jiǎn)諧激振、周期激振、脈沖激振、階躍激振和任意激振。

諧響應(yīng)分析［1-3］是用于分析持續(xù)的周期性載荷在結(jié)構(gòu)系統(tǒng)中產(chǎn)生的持續(xù)的周期響應(yīng)，以及確定線性結(jié)構(gòu)在承受隨時(shí)間按正弦（簡(jiǎn)諧）規(guī)律變化的載荷時(shí)，穩(wěn)態(tài)響應(yīng)的一種技術(shù)，一般用于確保一個(gè)給定的結(jié)構(gòu)能經(jīng)受住不同頻率的各種正弦載荷，探測(cè)其共振響應(yīng)，并在必要時(shí)避免其發(fā)生。

推力軸承［4］是連接螺旋槳軸與船舶的軸承，軸的力作用在推力軸承上，通過推力軸承推動(dòng)船舶運(yùn)動(dòng)。在航行的過程中，推力軸承基座實(shí)際上受力很大且很復(fù)雜，包括衡定的軸承的壓力、軸力以及軸的回轉(zhuǎn)對(duì)基座施加穩(wěn)態(tài)的、交變的作用力［5］。由于這些力都施加在基座的頂部面板上，使得基座頂端發(fā)生較大變形，但推力軸承對(duì)變形有較高要求，這就要求基座的設(shè)計(jì)必須保證其剛度。本文擬采用ANSYS軟件，對(duì)某實(shí)船的推力軸承基座及經(jīng)過加強(qiáng)處理的基座進(jìn)行帶預(yù)應(yīng)力的諧響應(yīng)分析［6-8］，驗(yàn)證其強(qiáng)度、剛度及振動(dòng)特性，分析試驗(yàn)時(shí)推力軸承基座在倒車時(shí)存在比較明顯的晃動(dòng)原因。

2 計(jì)算原理

系統(tǒng)受迫振動(dòng)的微分方程即運(yùn)動(dòng)方程［9］一般可以表述為：

當(dāng)系統(tǒng)受簡(jiǎn)諧激振力作用時(shí)，［K］矩陣和｛u｝矩陣是簡(jiǎn)諧的，其頻率為ω，因此：

式中，［M］是結(jié)構(gòu)質(zhì)量矩陣；［C］是結(jié)構(gòu)阻尼矩陣；［K］是結(jié)構(gòu)剛度矩陣；｛F｝是隨時(shí)間變化的載荷函數(shù)；｛u｝是節(jié)點(diǎn)位移矢量。

ANSYS軟件提供三種方法用以解決結(jié)構(gòu)諧響應(yīng)分析，分別是完整法、縮減法和模態(tài)疊加法。

完整法是三種方法中最容易使用的方法，采用完整的系統(tǒng)矩陣計(jì)算諧響應(yīng)。系統(tǒng)矩陣不存在矩陣縮減，它可以是對(duì)稱的，也可以是不對(duì)稱的。

縮減法通過使用縮減矩陣，比完整法計(jì)算更快。它需要選擇主自由度，據(jù)主自由度的位移值擴(kuò)展得到完整的DOF集。同時(shí)還采用縮減矩陣來(lái)壓縮問題的規(guī)模。

模態(tài)疊加法通過對(duì)模態(tài)分析得到的振形乘上參與因子并求和來(lái)計(jì)算出結(jié)構(gòu)的響應(yīng)，是所有求解方法中最快的。

本文采用模態(tài)疊加法對(duì)受預(yù)應(yīng)力的推力軸承基座進(jìn)行諧響應(yīng)分析，確定其穩(wěn)態(tài)動(dòng)力響應(yīng)曲線。

3 計(jì)算模型

3.1 模型范圍

為了避免邊界條件對(duì)計(jì)算結(jié)果的影響，在實(shí)際建模過程中，將推力軸承基座連同船體外板及相關(guān)構(gòu)件作為整體建模，模型縱向選取為基座兩端外推兩檔肋骨間距。具體模型見圖1和圖2。

3.2 單元選取及網(wǎng)格劃分

在建模過程中，船體外板、基座的各組面板和腹板以及船底縱桁及強(qiáng)肋骨等T型材均采用4節(jié)點(diǎn)殼單元SHELL63模擬；肋骨為12號(hào)球扁鋼，采用梁?jiǎn)卧狟EAM188模擬。模型中采用國(guó)際單位制，單元尺寸設(shè)定為以0.05 m為基準(zhǔn)。原模型共有17 000個(gè)殼單元、670個(gè)梁?jiǎn)卧?7 969個(gè)節(jié)點(diǎn)；加強(qiáng)后模型共有22 365個(gè)殼單元、690個(gè)梁?jiǎn)卧?3 163個(gè)節(jié)點(diǎn)。

3.3 材料參數(shù)

材料為B級(jí)鋼，其彈性模量為Ε＝2.06×1011N／m2，泊松比 υ＝0.3，密度為 7 850 kg／m3。

3.4 載荷分布及邊界條件

1）靜力分析

工作狀態(tài)下，基座頂部面板同時(shí)承受軸承帶來(lái)的1.5 t的垂直向下的重力，以及7 t的縱向力。這些力分布在基座面板的6個(gè)節(jié)點(diǎn)上，見圖3、圖4。此時(shí)邊界條件為：模型兩端所有節(jié)點(diǎn)剛固，模型舷頂所有節(jié)點(diǎn)簡(jiǎn)支。

2）模態(tài)計(jì)算

進(jìn)行模態(tài)計(jì)算時(shí)邊界條件與靜力分析一致。

3）諧響應(yīng)分析

諧響應(yīng)分析時(shí)，給基座頂部施加幅值為1 t的縱向簡(jiǎn)諧激振力，邊界條件與靜力分析時(shí)一致。

4 計(jì)算過程及結(jié)果

4.1 靜力分析

要對(duì)推力軸承基座進(jìn)行帶預(yù)應(yīng)力的諧響應(yīng)分析，首先需要進(jìn)行靜力分析。原推力軸承基座靜力分析應(yīng)力云圖見圖5，其位移云圖見圖6。加強(qiáng)后的推力軸承基座靜力分析應(yīng)力云圖見圖7，位移云圖見圖8。

從靜力計(jì)算結(jié)果可知，在7 t的縱向和1.5 t的垂向力的作用下，原基座整體模型應(yīng)力很大，且底部肋骨處存在應(yīng)力集中，基座局部應(yīng)力最大達(dá)到192 MPa，應(yīng)力也偏大；加強(qiáng)后基座的最大應(yīng)力值為59.5 MPa。可見，加強(qiáng)后基座力能有效傳遞出去，整體應(yīng)力值小，強(qiáng)度儲(chǔ)備充足且不存在應(yīng)力集中。原基座整體模型縱向位移達(dá)到10.4 mm，其變形值偏大；加強(qiáng)后基座最大位移值為0.594 mm，可見加強(qiáng)后基座剛度也足夠。

4.2 帶預(yù)應(yīng)力模態(tài)分析

對(duì)于考慮預(yù)應(yīng)力影響的模態(tài)分析，需要指定模態(tài)求解方法、提取的模態(tài)數(shù)以及要擴(kuò)展的模態(tài)數(shù)。對(duì)于本文，指定模態(tài)求解的方法為分塊蘭索斯法（Block Lanczos），并指定提取模態(tài)的階數(shù)為40。由于模型相對(duì)復(fù)雜，提取的模態(tài)中有很多階屬于局部的。對(duì)模型模態(tài)進(jìn)行分析和提取，得到模型前5階的固有頻率見表1。

表1 基座固有頻率（Hz）

4.3 諧響應(yīng)分析

諧響應(yīng)分析中需要指定分析的方法。本文將采用模態(tài)疊加法來(lái)進(jìn)行分析。另外，還需要對(duì)結(jié)果輸出形式進(jìn)行設(shè)置，以便程序?qū)⑿枰那蠼饨Y(jié)果寫到結(jié)果文件中。

在進(jìn)行諧響應(yīng)分析時(shí)還需要定義載荷步選項(xiàng)，包括諧響應(yīng)分析的頻率范圍、載荷子步數(shù)、載荷施加方式等選項(xiàng)。

對(duì)推力軸承基座進(jìn)行諧響應(yīng)分析得到的曲線如圖9、圖10所示。可見，對(duì)于原基座，在激振力頻率為 33.846 Hz和頻率為 47.957 Hz時(shí)，會(huì)發(fā)生諧響應(yīng)共振。對(duì)于加強(qiáng)后的推力軸承基座在激振力頻率為 79.974 Hz、85.334 Hz 和 142.024 Hz 時(shí)發(fā)生諧響應(yīng)共振。

由于本例螺旋槳的最大轉(zhuǎn)動(dòng)頻率為25 Hz，無(wú)論是原推力軸承基座還是加強(qiáng)后的基座發(fā)生諧響應(yīng)共振的頻率值都大于25 Hz。但是從圖11和圖12可知，對(duì)于原基座，在頻率為25 Hz的簡(jiǎn)諧激振力作用時(shí)，基座的諧響應(yīng)在3.3 mm左右；對(duì)于加強(qiáng)后基座，則減小至0.07 mm以下，這個(gè)值相當(dāng)小，可見加強(qiáng)后基座結(jié)構(gòu)在運(yùn)行時(shí)安全，不會(huì)再出現(xiàn)基座晃動(dòng)的問題。

5 結(jié) 論

1）從文中的靜力分析的計(jì)算結(jié)果來(lái)看，對(duì)比原推力軸承基座，加強(qiáng)后的基座在承受相同的靜載荷時(shí)，應(yīng)力和位移值都明顯小很多，而且不會(huì)在基座底部發(fā)生應(yīng)力集中的問題，基座結(jié)構(gòu)的強(qiáng)度和剛度儲(chǔ)備都較為充分。

2）從諧響應(yīng)分析計(jì)算結(jié)果來(lái)看，原推力軸承基座和加強(qiáng)后的基座在軸承轉(zhuǎn)速頻率范圍內(nèi)的簡(jiǎn)諧激振力作用下不會(huì)發(fā)生共振。但比較而言，加強(qiáng)后基座在激勵(lì)頻率作用下的響應(yīng)幅值大大降低。

3）類似于推力軸承基座這種旋轉(zhuǎn)設(shè)備的支座，一般來(lái)說(shuō)，對(duì)剛度的要求比對(duì)強(qiáng)度的要求高。首先，軸承的安裝精度對(duì)基座的剛度提出了要求。其次，基座剛度越大，其固有頻率越高，有利于避免在簡(jiǎn)諧激振力作用下的共振。

［1］朱從兵.基于有限元法的輕型港口起重機(jī)結(jié)構(gòu)動(dòng)態(tài)分析［J］.中國(guó)工程機(jī)械學(xué)報(bào)，2008，6（3）：333-339.

［2］姚學(xué)詩(shī)，周傳榮.基于預(yù)應(yīng)力法的轉(zhuǎn)子不平衡響應(yīng)研究［J］.振動(dòng)與沖擊，2005，24（2）：84-86.

［3］嚴(yán)大考.DF50／160 架橋機(jī)主梁的動(dòng)力諧響應(yīng)分析［J］.山西建筑，2009，35（10）：52-53.

［4］周渭鎬，包傳福，黃忠信.艦船汽輪機(jī)結(jié)構(gòu)和強(qiáng)度［M］.北京：國(guó)防工業(yè)出版社，1985.

［5］王維信.某型輔助船推力軸承故障研究與維修技術(shù)淺析［J］.中國(guó)修船，2009，22（2）：27-29.

［6］WANG Y Y.Finite element analysis and structural optimization of a vibration feeder based on ANSYS ［J］.Intern J.Plant Engineering&Management，2009，14（1）：43-47.

［7］馬宇，楊彩紅，馬永杰.基于ANSYS的連續(xù)梁橋的動(dòng)態(tài)特性分析［J］.武漢工程職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2008，20（3）：17-19.

［8］康方，范晉偉.基于ANSYS的數(shù)控機(jī)床動(dòng)態(tài)特性分析［J］.機(jī)械設(shè)計(jì)與制造，2008（7）：181-182.

［9］鄒經(jīng)湘.結(jié)構(gòu)動(dòng)力學(xué)［M］.哈爾濱：哈爾濱工業(yè)大學(xué)出版社，1996.