吳益輝,萬樂天

(中國船舶重工集團(tuán)有限公司第七一〇研究所,湖北宜昌443003)

0 引言

球柱組合殼結(jié)構(gòu)是單殼體潛艇耐壓船體結(jié)構(gòu)的主要形式，其極限承載能力是衡量結(jié)構(gòu)安全性的一個重要指標(biāo)。單殼體潛艇具有超長艙和大分艙特點(diǎn)[1]，在靜水壓力作用下球柱組合殼結(jié)構(gòu)的變形不再是單一的構(gòu)件受載荷作用下的變形，而是會受到相鄰結(jié)構(gòu)構(gòu)件變形的誘導(dǎo)，在連接處會產(chǎn)生應(yīng)力集中現(xiàn)象，使結(jié)合處成為整個結(jié)構(gòu)的薄弱環(huán)節(jié)，影響了整個結(jié)構(gòu)的承載能力。國內(nèi)外很多學(xué)者對柱殼或球殼的后屈曲變形和極限承載能力進(jìn)行了相關(guān)研究，而對球柱組合殼的研究主要集中在力學(xué)性能的分析方法上，但少涉及到球柱組合殼及其試驗(yàn)相關(guān)方法。白雪飛等[2]設(shè)計(jì)制作了3個凹型加肋錐—柱/錐—環(huán)—柱結(jié)合殼系列精車模型，靜水外壓試驗(yàn)和計(jì)算結(jié)果表明，凹型加肋錐—環(huán)—柱結(jié)合殼能有效降低錐、柱殼結(jié)合部的應(yīng)力峰值。郭日修等[3]通過分析加肋錐—環(huán)—柱結(jié)合殼中環(huán)殼塊的應(yīng)力分布特點(diǎn)，指出凹環(huán)殼塊需局部加強(qiáng)的必要性，并提出了凹環(huán)殼塊中部增設(shè)肋骨和凹環(huán)殼塊適當(dāng)加厚兩種加強(qiáng)方式。ABS[4]及DNV[5]規(guī)范中給出了相對準(zhǔn)確和可靠的水下航行器的相關(guān)設(shè)計(jì)方法和失效載荷；Radha和Rajagopalan[6]采用非線性有限元及蒙特卡洛法分析潛艇耐壓殼體的承載可靠性，并利用Johnson-Ostenfeld彈塑性修正法及非線性有限元逼近法研究了加筋圓柱殼的彈塑性失穩(wěn)，表明有限元法的精度高于理論計(jì)算。而關(guān)于初撓度的研究，往往通過導(dǎo)入其特征屈曲失穩(wěn)模態(tài)波形實(shí)現(xiàn)，雖取得了較好的結(jié)果，但其失穩(wěn)波形存在一定的假設(shè)，可能與實(shí)際情況不符，并且殼體結(jié)構(gòu)在建造過程中其初撓度是很難用解析函數(shù)表示的[7]。

本文根據(jù)實(shí)測三維初撓度數(shù)據(jù)建立含缺陷的球柱組合殼有限元模型，計(jì)算了外部均勻壓力作用下球柱組合殼的極限承載能力和變形，并與試驗(yàn)結(jié)果進(jìn)行了對比分析，得到了組合殼體的變形與受力特性，為球柱組合殼體的優(yōu)化設(shè)計(jì)提供了有效的依據(jù)。

1 球柱組合殼體結(jié)構(gòu)力學(xué)性能

球柱組合殼結(jié)構(gòu)在受均勻靜水壓力作用下的力學(xué)模型如圖1所示。球冠結(jié)構(gòu)邊緣半徑與軸線之間夾角為α，本文球柱組合殼體端面利用法蘭固定，加壓時力學(xué)模型可簡化如此。

圖1 球柱組合殼力學(xué)模型Fig.1 Mechanicalmodel of sphere-cylinder combined shell

2 球柱組合殼靜水壓力試驗(yàn)

球柱組合殼靜水壓力試驗(yàn)在模擬深海的壓力罐內(nèi)進(jìn)行，3個球柱組合殼結(jié)構(gòu)的焊接模型，分別命名為1#模型、2#模型和3#模型。這3個模型的圓柱殼結(jié)構(gòu)尺度相同，球形殼殼板厚度與圓柱殼相同。1#、2#模型為半球殼結(jié)構(gòu)，3#模型為球冠結(jié)構(gòu)，球冠殼的圓弧角為120°，采用逐步加載法進(jìn)行了強(qiáng)度試驗(yàn)和承載能力試驗(yàn)[8]。

分別利用回轉(zhuǎn)平臺、分度儀和高度尺等對3個球柱組合殼的壁厚及肋骨尺寸和肋骨間距進(jìn)行了測量，其中殼體厚度及肋骨間距測量結(jié)果如表1-2所示。表1殼體厚度測量中涉及周向6個角度下的測量均值，取平均為最終殼體厚度；表2中l(wèi)i(i=1，2，3，4，5)表示肋骨離球柱組合殼連接處及相鄰肋骨間的距離，其中2#模型有4根環(huán)肋。

表1 模型殼體厚度測量均值Table 1 Averagemeasurement values of models’hull thickness

表2 模型肋距測量均值Table 2 Averagemeasurement values ofmodels’frame space

3 非線性有限元分析

3.1 模型建立

根據(jù)模型結(jié)構(gòu)實(shí)測參數(shù)建立了相應(yīng)含缺陷幾何模型，殼體厚度及環(huán)肋加強(qiáng)筋由測量均值給出，組合殼半徑R=155mm，肋骨截面為矩形且腹板長為l=10 mm，法蘭半徑r=195 mm，法蘭板厚t=5.2 mm。球柱組合殼結(jié)構(gòu)焊接模型材料彈性模量E=1.97×105MPa，泊松比μ=0.3，屈服極限σs=418 MPa。

3.2 數(shù)值仿真結(jié)果

試驗(yàn)及有限元結(jié)果見圖2，由變形形狀及應(yīng)力集中區(qū)分布情況來看，兩者符合的很好。1#模型最大變形發(fā)生在球殼與柱殼連接處靠近第一個環(huán)肋處；2#模型最大變形發(fā)生在柱殼第二與第三環(huán)肋中間處；3#模型最大變形則發(fā)生在第一與第二環(huán)肋中間處；變形區(qū)域平均 Mises應(yīng)力分別達(dá)到539.5 MPa、476.9 MPa、442.7 MPa，可見，在外部均布水壓作用下，組合殼結(jié)構(gòu)均失穩(wěn)破壞。

圖2 變形及Mises應(yīng)力分布Fig.2 Deformations and Mises stress distributions

試驗(yàn)中3個模型均發(fā)生總體失穩(wěn)破壞，且失穩(wěn)的部分為環(huán)肋圓柱殼部分。最終測得1#、2#、3#試驗(yàn)?zāi)Ｐ推茐谋罎毫Ψ謩e為2.2 MPa、2.4 MPa、2.3 MPa，而有限元結(jié)果為 2.5 MPa、2.8 MPa、2.4 MPa?？梢?，組合殼2#模型有著較高的承載能力，與1#模型相比，靠近球柱連接處的環(huán)肋提高了結(jié)構(gòu)的整體穩(wěn)定性，其應(yīng)力集中區(qū)較??；而3#模型與1#模型的抗壓能力相當(dāng)，其球冠結(jié)構(gòu)使得變形區(qū)域下移，在連接處較大的結(jié)構(gòu)不連續(xù)性使其產(chǎn)生了比1#模型更大的塑性應(yīng)變；3個模型球殼部分應(yīng)力分布較為均勻。

4 對比分析研究

圖3 有限元與試驗(yàn)變形對比Fig.3 Results comparison of finite element and experiment deformation

1#與2#球柱組合殼體極限承載能力分別為554.9 kn、551.2 kn，觀察結(jié)構(gòu)達(dá)到極限強(qiáng)度之后的后極限強(qiáng)度行為變化過程可以發(fā)現(xiàn)，2#模型的承載能力變化較1#模型緩慢，且具有較高的剩余強(qiáng)度，可見在連接處的加強(qiáng)筋對于球柱組合殼的后極限強(qiáng)度行為影響較大。3#模型承載能力變化與1#模型相當(dāng)，但是最大抵抗反力相比1#模型較小，僅為539.8 kn，可見球冠結(jié)構(gòu)的不連續(xù)性減弱了結(jié)構(gòu)的極限強(qiáng)度。

輸出和測量了縱向沿半殼體經(jīng)線分別于球殼頂部、中部、球柱殼連接處及柱殼肋骨及肋骨中間11個位置點(diǎn)的軸向變形，橫向?qū)ΨQ沿半殼體緯線方向的9個位置點(diǎn)的軸向變形，均以最大變形點(diǎn)為基準(zhǔn)點(diǎn)，且變形方向均為沿殼體半徑方向。非線性有限元與實(shí)驗(yàn)結(jié)果如圖3所示。由圖3可知，1#與3#模型縱橫向變形較2#模型平緩，2#模型由小變形到極限變形發(fā)展速度很快。1#與3#球柱組合殼連接處的3號位置點(diǎn)變形量較大，而2#球柱組合殼3號位置點(diǎn)變形量較小，說明靠近連接處的柱殼環(huán)肋加強(qiáng)在很大程度上影響了連接處的變形與應(yīng)力分布。

表3為球柱組合殼最大變形數(shù)值模擬與試驗(yàn)結(jié)果對比，其誤差分別為0.35%、19.61%、9.75%。2#模型較其他2個模型的誤差偏大，可能是由于測量誤差和有限元處理大變形時迭代路徑的影響，但總體上可以看出非線性有限元結(jié)果與實(shí)驗(yàn)誤差在工程結(jié)構(gòu)物計(jì)算的可接受范圍內(nèi)，較為準(zhǔn)確地反映了結(jié)構(gòu)的實(shí)際變形。

表3 最大變形結(jié)果對比Table 3 Comparison ofmaximum deformation results

5 結(jié)束語

本文通過數(shù)值模擬與試驗(yàn)對比分析研究了3個不同形式球柱組合殼的極限承載能力和變形，得出了以下結(jié)論：

1)3個模型的破壞都屬于柱殼部分的總體失穩(wěn)破壞；

2)通過實(shí)測初始缺陷建立的球柱組合殼有限元模型可以有效提高數(shù)值模擬的有效性與準(zhǔn)確性，相比于通過屈曲失穩(wěn)波形引入缺陷的方法更加符合實(shí)際；

3)環(huán)肋提高了結(jié)構(gòu)的整體穩(wěn)定性，但其對球柱組合殼結(jié)構(gòu)的極限強(qiáng)度影響不大，而對后極限強(qiáng)度行為影響較大，可以通過在離球柱組合殼連接處不同位置布置相應(yīng)肋骨來減弱過渡區(qū)域?qū)Y(jié)構(gòu)的影響。