李 凱，朱鎮(zhèn)揚(yáng)，連雪海

(1.江蘇科技大學(xué) 機(jī)械工程學(xué)院，江蘇 鎮(zhèn)江 212000；2.中電科(寧波)海洋電子研究院有限公司，浙江 寧波 315040)

0 引言

長(zhǎng)江是我國(guó)重要的運(yùn)輸河流，但目前面臨繁重的生態(tài)環(huán)境修復(fù)任務(wù)，急需使用潛水器代替人工完成高危的水下工作。而耐壓殼作為潛水器的關(guān)鍵結(jié)構(gòu)，發(fā)揮著承受水壓、保護(hù)內(nèi)部元件的作用，占潛水器總質(zhì)量的1/4～1/2。

現(xiàn)有耐壓殼多為球形，雖然強(qiáng)度高，但空間利用率低，且在受載時(shí)對(duì)缺陷敏感。PAN等發(fā)現(xiàn)球形殼易在理論值壓力的1/4～1/3處發(fā)生失穩(wěn)，因此潛水器只能通過(guò)增大半徑拓展空間，但又會(huì)帶來(lái)水阻力增大的問(wèn)題?；诤｝斖鈿ど锾匦栽O(shè)計(jì)的耐壓殼與球形殼相比，具有以下優(yōu)勢(shì)：外形呈流線型，具有較小水阻力；扁平狀結(jié)構(gòu)使其具有更高的穩(wěn)定性，有助于在湍急的江水中保持平衡；曲率相對(duì)小，空間利用率高，可使?jié)撍鞔钶d更多功能模塊。

為此，本文進(jìn)行龜型耐壓殼的設(shè)計(jì)，研究等屈曲設(shè)計(jì)下的殼體厚度變化規(guī)律，最后基于數(shù)值法進(jìn)行屈曲分析。

1 耐壓殼設(shè)計(jì)

海龜外殼扁平，中間寬、兩端窄，截面上呈橢圓弧狀。在此基礎(chǔ)上，考慮用橢球形殼與球形殼的結(jié)合相擬合。在直角坐標(biāo)系下，本耐壓殼由方程(1)繞軸旋轉(zhuǎn)180°而成，見(jiàn)圖1。圖中：為橢圓半長(zhǎng)軸；為橢圓半短軸；為圓心與圓心′之間的距離；為圓的半徑；′為圓心；為橢圓與圓的切點(diǎn)；為橢圓上點(diǎn)，坐標(biāo)為(cos，sin)；為以為半徑的圓上的點(diǎn)，坐標(biāo)為(cos，sin)；為以為半徑的圓上的點(diǎn)，坐標(biāo)為(cos，sin)；與互為垂線，為垂足；為∠′。

圖1 耐壓殼幾何參數(shù)

(1)

式中：為橢圓與圓方程相交的橫坐標(biāo)；為所截取殼體部分距原點(diǎn)的縱坐標(biāo)。

選擇耐壓殼基本參數(shù)為：為200 mm，為120 mm，為50 mm，為110 mm，為80 mm，為40 mm。

耐壓殼等厚設(shè)計(jì)會(huì)導(dǎo)致安全裕度過(guò)剩，同時(shí)浪費(fèi)原材料，成本高。等屈曲變厚設(shè)計(jì)可以減少殼體材料的使用，在降低成本的同時(shí)，又保證良好的抗壓和抗沖擊能力。殼體由橢球殼部分與球殼部分共同組成，模型見(jiàn)圖2。

圖2 耐壓殼模型

對(duì)于橢球殼部分，橢球體第一曲率半徑、第二曲率半徑的表達(dá)式在通過(guò)橢圓參數(shù)方程：

(2)

簡(jiǎn)化后，可表示為

(3)

在線性范圍內(nèi)，根據(jù)Mushtari的屈曲理論公式，屈曲載荷可以用下式計(jì)算：

(4)

對(duì)于球殼部分，Zoelly的經(jīng)典屈曲理論公式，屈曲載荷可以用下式計(jì)算：

(5)

式中：為彈性模量；為泊松比；為殼厚。

現(xiàn)對(duì)橢球殼部分等屈曲變厚，將式(3)代入式(4)，反解出，得等屈曲橢球殼厚度表達(dá)式：

(6)

式中：=sin+cos；為殼體所受到的水壓，由于殼體自身的垂直高度(-)與水深相比可以忽略不計(jì)，因此取定值。

通過(guò)求導(dǎo)，發(fā)現(xiàn)以下厚度變化規(guī)律：

(7)

現(xiàn)對(duì)球殼部分等屈曲變厚，對(duì)式(5)反解出，得等屈曲橢球殼厚度表達(dá)式：

(8)

在水壓、材料和半徑確定時(shí)，()為常數(shù)。

2 屈曲數(shù)值分析

2.1 數(shù)值模型

在設(shè)計(jì)中，考慮到長(zhǎng)江流域水深，要求耐壓殼能夠承受50 m深的水壓，相應(yīng)水壓=0.49 MPa，并取安全系數(shù)=1.5。選擇樹(shù)脂材料，其性能參數(shù)為：=2 400 MPa，屈服強(qiáng)度[]=33.22 MPa，泊松比=0.41，=3.6 mm。采用前處理軟件ANSA對(duì)網(wǎng)格進(jìn)行劃分。應(yīng)力分析模型多采用四邊形殼單元，共獲得13 791個(gè)網(wǎng)格，見(jiàn)圖3。

圖3 網(wǎng)格劃分

考慮到安全系數(shù)，對(duì)殼體外表面施加均布載荷=0.735 MPa。約束定義：限制底部邊界所有節(jié)點(diǎn)的平移和旋轉(zhuǎn)的自由度，即===0且===0。

2.2 數(shù)值計(jì)算結(jié)果分析

利用ABAQUS軟件對(duì)耐壓殼模型進(jìn)行線性屈曲求解，得到前2階屈曲模態(tài)。前2階模態(tài)云圖見(jiàn)圖4。在線性屈曲分析下，耐壓殼模型臨界屈曲載荷為0.79 MPa，高于水壓0.735 MPa，材料符合要求。通過(guò)云圖發(fā)現(xiàn)，屈曲容易發(fā)生在殼體的上下及中間部分。由于材料制造過(guò)程存在初始缺陷，真實(shí)部件可能在低于理論值的載荷下發(fā)生屈曲，因此線性屈曲分析下的臨界載荷相對(duì)偏大。

圖4 前2階模態(tài)云圖

采用ABAQUS軟件的Riks算法進(jìn)行非線性屈曲分析，其應(yīng)變?cè)茍D見(jiàn)圖5。由圖5可知：在非線性屈曲分析下，其臨界屈曲載荷為0.76 MPa，高于水壓0.735 MPa，材料符合要求。其他結(jié)論與線性屈曲分析相同。

圖5 應(yīng)變?cè)茍D(弧長(zhǎng)=1.41)

3 結(jié)論

(1)由殼體厚度表達(dá)式可知，無(wú)論橢圓率在哪個(gè)范圍，橢球形殼體部分的厚度在=π/2處取最大值，球形殼體部分的厚度為常數(shù)。

(2)基于數(shù)值法的屈曲分析結(jié)果表明，龜型耐壓殼的上下及中間部分為殼體的薄弱環(huán)節(jié)，最容易發(fā)生破壞。