孫 倩 廖建彬 蔡振雄

集美大學(xué)輪機(jī)工程學(xué)院 廈門361021

潛艇耐壓液艙區(qū)域有限元應(yīng)力計算

孫 倩 廖建彬 蔡振雄

集美大學(xué)輪機(jī)工程學(xué)院 廈門361021

針對目前潛艇耐壓液艙結(jié)構(gòu)應(yīng)力解析公式精度不高的情況，根據(jù)其受力特點建立了同心圓和準(zhǔn)同心圓式耐壓液艙結(jié)構(gòu)有限元分析的力學(xué)模型，編程通過實例計算表明自編有限元程序可靠性很高，可單獨作為一個模塊取代傳統(tǒng)的近似解析法進(jìn)行應(yīng)力分析及后續(xù)的優(yōu)化設(shè)計。

潛艇 船體結(jié)構(gòu) 有限元 程序

潛艇耐壓液艙結(jié)構(gòu)受力復(fù)雜，試驗表明近似解析法不能反映結(jié)構(gòu)內(nèi)的真實應(yīng)力狀態(tài)。盡管近年來國內(nèi)外對此展開了廣泛的分析研究，先后提出了一些改進(jìn)的計算方法1-5，但都不夠完善。因此在對受力特性最好的同心圓式耐壓液艙結(jié)構(gòu)進(jìn)行了大量的有限元分析的基礎(chǔ)上自編了有限元應(yīng)力計算的源程序，作為一個模塊進(jìn)行最終的優(yōu)化設(shè)計和應(yīng)力校核。

1 結(jié)構(gòu)型式選擇

傳統(tǒng)耐壓液艙結(jié)構(gòu)（見圖1 a））耐壓船體應(yīng)力集中嚴(yán)重、軸向應(yīng)力很大，同心圓式耐壓液艙結(jié)構(gòu)（見圖1 b））可以克服這些缺點，但是這種結(jié)構(gòu)型式在布置和通行上會造成不便，以致工程上不太實用。研究表明，如果將同心圓式耐壓液艙結(jié)構(gòu)上部的圓弧形液艙殼板改成平板，構(gòu)成準(zhǔn)同心圓式耐壓液艙結(jié)構(gòu)（見圖1 c）），通過計算，除平臺區(qū)局部外，基本不改變同心圓式耐壓液艙結(jié)構(gòu)的應(yīng)力分布狀態(tài)，卻可克服上述缺點。相關(guān)試驗也證明［6］：在彈性階段，準(zhǔn)同心圓式耐壓液艙與傳統(tǒng)耐壓液艙相比，前者船體殼板的應(yīng)力沿圓周方向趨于均勻化，最大應(yīng)力明顯低于后者相應(yīng)位置處的同類應(yīng)力，在外載增長過程中，前者達(dá)到局部屈服的時間比后者長，從而驗證了準(zhǔn)同心圓式耐壓液艙結(jié)構(gòu)是一種優(yōu)良的結(jié)構(gòu)型式。

圖1 耐壓液艙結(jié)構(gòu)簡圖

2 力學(xué)模型的建立

考慮到以文獻(xiàn)［1-5］為代表的近似解析法一般均采用軸對稱假設(shè)，即將傳統(tǒng)的耐壓液艙結(jié)構(gòu)簡化為一個中間由實肋板相連的雙層圓柱殼來分析潛艇耐壓液艙的應(yīng)力分布情況，同心圓和準(zhǔn)同心圓式耐壓液艙結(jié)構(gòu)更符合軸對稱實際，因此也可采用此假設(shè)。

2．1 力學(xué)模型的選取

以計算精度論，應(yīng)采用全剖面力學(xué)模型，但是它的計算量太大，不利于編程，所以采用中心對稱力學(xué)模型見圖2，即縱向取一檔實肋板間距，周向取一個縱骨間距夾角內(nèi)的結(jié)構(gòu)。此外，通過ANSYS軟件計算，驗證了平臺板的應(yīng)力如采用解析公式，即應(yīng)力按四周剛性固定板計算［7］，與有限元結(jié)果吻合良好。因此對于同心圓柱殼結(jié)構(gòu)自身的對稱性，完全可采用這一力學(xué)模型進(jìn)行計算，而準(zhǔn)同心圓式結(jié)構(gòu)的平臺板可用解析公式計算，其他部位的應(yīng)力可采用本文的有限元法的計算結(jié)果。

圖2 全實肋板加縱骨式耐壓液艙結(jié)構(gòu)中心對稱力學(xué)模型

2．2 有限元網(wǎng)格劃分

將耐壓液艙和耐壓船體殼板周向劃分成4個單元，縱向劃分成8個單元，其他結(jié)構(gòu)均劃分成相當(dāng)大小的單元?？紤]到結(jié)構(gòu)為圓柱形薄殼結(jié)構(gòu)，故采用平面四邊形單元，除實肋板采用四邊形線性等參元外，其余構(gòu)件均采用矩形單元。

2．3 邊界條件的處理

位移邊界條件的選?。簝啥藱M剖面上沿圓周的切線和法線方向轉(zhuǎn)角為零，不受外力一端約束縱向位移；在夾角扇面上沿殼板平面切線和法線轉(zhuǎn)角為零；沿夾角剖面的法線位移為零；在耐壓液艙縱骨和船體縱骨上三個方向的轉(zhuǎn)角均為零；約束各板單元面內(nèi)的第六自由度。

在外力方面，假設(shè)液艙殼板受法向壓力，縱向力則以集中力的形式按力的等效作用原理分配到端面內(nèi)的各結(jié)點上。計算公式參考文獻(xiàn)［1］。

3 應(yīng)力有限元計算源程序編制

根據(jù)上述有限元力學(xué)模型，用C語言編制有限元應(yīng)力計算程序。程序流程見圖3。

圖3 采用平板單元分析耐壓液艙結(jié)構(gòu)流程圖

3．1 程序編制目的

程序編制的主要目的是為以后的優(yōu)化設(shè)計打基礎(chǔ)，因為盡管大型有限元軟件如ANSYS、ADINA和NASTRAN等精度功能都比較突出，但明顯不適用于在優(yōu)化過程中的多次頻繁調(diào)用，而近似解析法的精度又不能令人滿意，所以有必要編制一個獨立性強(qiáng)并具有較好實用性的有限元程序嵌入到約束函數(shù)子程序中取代解析法來計算關(guān)鍵點應(yīng)力。另外還可以在初步設(shè)計階段可用于校核耐壓液艙結(jié)構(gòu)應(yīng)力。

3．2 程序介紹

該程序由1個主程序、32個子程序組成，可以進(jìn)行同心圓和準(zhǔn)同心圓式耐壓液艙結(jié)構(gòu)的有限元應(yīng)力計算，使用者只需輸入計算載荷以及外形尺寸就可以得出各關(guān)鍵點的應(yīng)力及位移，另外各子程序功能獨立，便于移植，使用方便，只要在約束條件子函數(shù)中調(diào)用該程序就可以進(jìn)行耐壓液艙結(jié)構(gòu)的優(yōu)化設(shè)計。程序流程見圖3。

3．3 程序主要特點

程序?qū)φw剛度矩陣采用一維變帶寬存儲，可放松對結(jié)點編號的要求，且適宜用LDLT解法，對受約束的結(jié)點方向不建立方程。因此，整體剛度矩陣的階數(shù)為結(jié)構(gòu)的自由度數(shù)，自由度序號按結(jié)構(gòu)的結(jié)點編號順序排列。

由于采用了四結(jié)點二維等參元，所得到的應(yīng)力值是單元內(nèi)高斯積分點處的值，而實際需要的是結(jié)點值，因此程序采用了E．Hinton和D．R．J．Owen所提出的局部應(yīng)力平滑技術(shù)［8］，就是一單元的2×2個高斯點的雙線性外插應(yīng)力值，最后按光滑過的應(yīng)力再求繞結(jié)點的平均值，從而得到一個唯一的而且是較滿意的結(jié)點應(yīng)力。

對第六自由度采用如下處理：對于局部坐標(biāo)系的單元剛陣，不管周圍單元是否共面，可在其對應(yīng)于θzi的剛陣對角元素上都加上一個非常小的值ε［9］。這樣，當(dāng)i結(jié)點周圍單元都共面時，對應(yīng)于θzi的對角元素也不為零，消除了奇異性；而當(dāng)i結(jié)點周圍單元不共面時，由于ε值非常小，也不會影響正常剛度元素迭加結(jié)果。這雖然是一種簡化的近似處理方法，但實際應(yīng)用是很方便有效的。

4 實例驗證及計算

為驗證程序的可靠性，進(jìn)行若干實例計算。

算例1：全實肋板雙層加縱骨同心圓柱殼式耐壓液艙結(jié)構(gòu)有限元分析。程序計算結(jié)果和采用ANSYS軟件的計算結(jié)果見表1，兩種方法均采用中心對稱力學(xué)模型。

算例2：全實肋板不加縱骨同心圓柱殼式耐壓液艙結(jié)構(gòu)有限元分析。

原始數(shù)據(jù)同算例1，只是船體及液艙殼板均不加縱骨，ANSYS采用三擋肋位全剖面有限元力學(xué)模型，近似解析法采用文獻(xiàn)［1］中的方法，計算結(jié)果見表2。

算例3：全實肋板液艙殼板加縱骨同心圓柱殼式耐壓液艙結(jié)構(gòu)有限元分析。

原始數(shù)據(jù)同算例1，只有液艙殼板加縱骨，ANSYS采用三擋肋位全剖面有限元模型，計算結(jié)果見表3。

算例4：全實肋板液艙殼板加縱骨準(zhǔn)同心圓柱殼式耐壓液艙結(jié)構(gòu)有限元分析。

此算例計算結(jié)果與ANSYS計算結(jié)果見表4，ANSYS采用三擋肋位全剖面有限元模型。

表1 全實肋板雙層加縱骨典型點結(jié)構(gòu)應(yīng)力MPa

表2 液艙殼板均不加縱骨典型點結(jié)構(gòu)應(yīng)力MPa

表3 液艙殼板加縱骨典型點結(jié)構(gòu)應(yīng)力MPa

表4 全實肋板液艙殼板加縱骨典型點結(jié)構(gòu)應(yīng)力MPa

5 結(jié)論

對照表1～4可以看出：

1）僅在耐壓液艙殼板加縱骨對耐壓船體的應(yīng)力影響較小，但可以降低耐壓液艙的應(yīng)力；對于雙層加縱骨的情況，耐壓船體和耐壓液艙殼板的軸向應(yīng)力都有非常明顯的下降，耐壓液艙的周向應(yīng)力也有所下降，但耐壓船體的周向應(yīng)力則有所升高。由此說明，殼板加縱骨可顯著降低殼板的軸向應(yīng)力，但有可能造成周向應(yīng)力惡化。因此，對于同心圓柱殼耐壓液艙結(jié)構(gòu)，在耐壓液艙上加縱骨是比較好的措施，但在耐壓船體上加縱骨就得慎重：如果耐壓船體殼板上軸向應(yīng)力不高，則不必在耐壓船體殼板上加縱骨，以免造成周向應(yīng)力惡化。

2）近似解析法計算出來的耐壓船體殼板軸向應(yīng)力明顯偏高而周向應(yīng)力與有限元結(jié)果還比較吻合，文中有限元程序的計算結(jié)果與ANSYS的結(jié)果比較吻合，兩者均采用中心對稱力學(xué)模型時結(jié)果非常接近。誤差主要是由于文中基于計算量考慮選用了中心對稱力學(xué)模型，而ANSYS軟件計算考慮到實際情況和世界影響采用的多為較密網(wǎng)格多跨距全剖面模型，另外線性方程的求解、等參元的選擇以及第六自由度處理方法等方面都會造成一定的偏差。

綜上所述，采用文中編制的有限元程序進(jìn)行耐壓液艙結(jié)構(gòu)的強(qiáng)度分析，誤差小于4%，在初步設(shè)計階段可用于校核耐壓液艙結(jié)構(gòu)應(yīng)力，還可以作為單獨一個模塊取代解析公式進(jìn)行優(yōu)化設(shè)計。

［1］GJB／Z21-2001．潛艇結(jié)構(gòu)設(shè)計計算方法［S］．國防科學(xué)技術(shù)工業(yè)委員會批準(zhǔn)發(fā)布，2002：8-27．

［2］沈 豐．潛艇縱骨式全實肋板耐壓液艙結(jié)構(gòu)強(qiáng)度計算方法（Ⅰ）［J］．艦船性能研究，1994（1）：33-38．

［3］沈 豐．潛艇縱骨式全實肋板耐壓液艙結(jié)構(gòu)強(qiáng)度計算方法（Ⅱ）［J］．艦船性能研究，1994（3）：42-46．

［4］羅 斌，陳 強(qiáng)，馬 騁．縱骨式全實肋板耐壓液艙的縱骨計算方法探討［J］．艦船科學(xué)技術(shù)，2004（2）：19-20．

［5］羅 斌，陳 強(qiáng)，徐云椿．潛艇縱骨式全實肋板耐壓液艙殼板強(qiáng)度計算方法研究［J］．中國造船，1999（2）：74-80．

［6］謝祚水，孫 倩．耐壓液艙結(jié)構(gòu)的模型試驗研究［J］．華東船舶工業(yè)學(xué)院學(xué)報，2000（5）：1-4．

［7］陳鐵云，陳伯真．船舶結(jié)構(gòu)力學(xué)［M］．上海：上海交通大學(xué)出版社，1991（7）：197-210．

［8］楊菊生，攬生瑞．有限元法程序設(shè)計［M］．西安：西安交通大學(xué)出版社，1990（6）：21-103．

［9］趙經(jīng)文，王宏鈺．結(jié)構(gòu)有限元分析［M］．哈爾濱：哈爾濱工業(yè)大學(xué)出版社，1998（4）：23-39．

FE stress analysis of the pressure tank structure of submarine

SUN Qian LIAO Jian-bin CAI Zhen-xiong
Marine Engineering Institute Jimei University Xiamen 361021

Because of lower accuracy of analytic method，the FE calculation models of the homocentric and quasi-homocentric pressure tank structures were made alternatively according to their properties of loading．The numerical results showed that the programs have high reliability and are suitable for stress analyses optimum calculations in stead of traditional methods．

submarines hull structure finite element method program

U663

A

1671-7953（2007）02-0017-04

2006-08-02

修回日期2006-11-06

孫 倩（1974—），女，碩士，講師。