基于分層分壓結(jié)構(gòu)的新型潛水器耐壓殼結(jié)構(gòu)設(shè)計(jì)1)

張 吟 劉小明 雷現(xiàn)奇 孫成奇 方 新 魏宇杰

*(中國科學(xué)院力學(xué)研究所非線性力學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京100190)

?(中國科學(xué)院大學(xué)工程科學(xué)學(xué)院，北京100049)

基于分層分壓結(jié)構(gòu)的新型潛水器耐壓殼結(jié)構(gòu)設(shè)計(jì)1)

張 吟*,?,2)劉小明*,?雷現(xiàn)奇*,?孫成奇*方 新*魏宇杰*,?

*(中國科學(xué)院力學(xué)研究所非線性力學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京100190)

?(中國科學(xué)院大學(xué)工程科學(xué)學(xué)院，北京100049)

耐壓殼是深海潛水器中最關(guān)鍵的結(jié)構(gòu)，直接關(guān)系到潛水器安全性和總體性能.本文對當(dāng)今的單層耐壓殼結(jié)構(gòu)設(shè)計(jì)進(jìn)行了評述，并設(shè)計(jì)了一種基于分層/分壓的新型耐壓殼結(jié)構(gòu)，該結(jié)構(gòu)借鑒了自然界的兩種深潛動物的結(jié)構(gòu)：抹香鯨分層結(jié)構(gòu)和鸚鵡螺隔片分割螺殼亞結(jié)構(gòu).綜合了這兩種結(jié)構(gòu)特性的雙層殼結(jié)構(gòu)能夠有效提高抗壓能力，從而提升深潛能力.與以往單層球殼的耐壓殼結(jié)構(gòu)相比，該結(jié)構(gòu)不僅提高了強(qiáng)度，也提高了抗屈曲能力.同時(shí)，該結(jié)構(gòu)還兼具大容積、高可靠性、以及避免超厚殼制備上的難點(diǎn)等特征，使得深海潛水器的綜合性能得到顯著提高.新結(jié)構(gòu)中的桁架將圓殼分割為若干個(gè)柱殼亞結(jié)構(gòu)，本文針對此亞結(jié)構(gòu)嚴(yán)格推導(dǎo)了桁架增強(qiáng)殼體抗屈曲的公式.從實(shí)驗(yàn)數(shù)據(jù)中總結(jié)出來的泰勒水池公式是目前廣泛使用的潛水器殼體設(shè)計(jì)依據(jù)，新推導(dǎo)的公式與之相比只有6%的差別，這使得新型結(jié)構(gòu)設(shè)計(jì)有了更堅(jiān)實(shí)的理論基礎(chǔ).

耐壓殼，深潛器，結(jié)構(gòu)，強(qiáng)度，屈曲

引言

迄今為止,人類研究的海洋區(qū)域不足其總面積的10%，以至于尚不能精確繪制出一幅全景式海底地圖[1].由于深海中蘊(yùn)藏巨大的生物、礦產(chǎn)、石油、天然氣和可燃冰等資源，美、俄和日等海洋強(qiáng)國一直以來為確保其基于海洋資源開發(fā)利用而形成的巨大經(jīng)濟(jì)利益，都在大力推進(jìn)深海探測的研究，其中各國對新型深海載人和無人潛水器的研制和應(yīng)用最為重視.近年來，我國也加大了深海探測的研發(fā)力度，但與航空航天科技的發(fā)展速度相比，我國海洋研究的進(jìn)展是相對緩慢的.制約海洋研究進(jìn)展速度的一個(gè)主要瓶頸因素就是深海潛水器的研制.受限于材料與工程科學(xué)的發(fā)展，世界各國現(xiàn)在仍然無法制造能長期承受數(shù)百乃至上千大氣壓的載人深海探測設(shè)備.而載人深海潛水器核心部件耐壓殼的設(shè)計(jì)成為成敗關(guān)鍵所在.

耐壓殼，顧名思義，就是要承受海水巨大壓力的外層結(jié)構(gòu).深海潛水器耐壓殼結(jié)構(gòu)的設(shè)計(jì)最早可追溯到1928—1929年間由美國工程師奧梯斯·巴頓(Otis Barton)設(shè)計(jì)的“深球”號 (Bathysphere)深海潛水器[2].“深球”號于1930—1934年間由巴頓和博物學(xué)家威廉·比伯(William Beebe)共同駕駛在百慕大群島附近海域進(jìn)行了一系列深潛，于1934年8月15日創(chuàng)造了深度1023m的潛水記錄，幫助海洋生物學(xué)家實(shí)現(xiàn)了在棲息地對深海生物進(jìn)行直接觀測[2].早期另一個(gè)重要的深海潛水器耐壓殼結(jié)構(gòu)是瑞士人奧古斯特·皮卡德(Auguste Piccard)設(shè)計(jì)的“深船”號(Bathyscaphe)深海潛水器[3].“深船”號于1948年11月3日下水，下潛深度達(dá)到了1370m.“深球”號和“深船”號的成功研發(fā)開啟了人類向深海探索的一個(gè)嶄新紀(jì)元.至今世界主要國家的深潛器耐壓殼結(jié)構(gòu)均采用“深球”號和“深船”號的單層球殼的結(jié)構(gòu)設(shè)計(jì).表1是各國設(shè)計(jì)的深潛器耐壓殼結(jié)構(gòu)關(guān)鍵參數(shù).

表1 各國深潛器耐壓殼結(jié)構(gòu)關(guān)鍵參數(shù)Table 1 The key parameters of the pressure hulls of deep-sea submersibles designed by various countries

目前下潛最深的深潛器為美國的“深海挑戰(zhàn)者”號(Deepsea Challenger).在深度10908m的馬里亞納海溝，球殼外壁承受的海水壓力約為110MPa，球殼結(jié)構(gòu)不同地方承受的應(yīng)力遠(yuǎn)遠(yuǎn)超過110MPa.以深海挑戰(zhàn)者球殼為例，內(nèi)徑1.1m，殼厚66mm在外壓110MPa，內(nèi)壓0.1MPa的環(huán)境下，球殼的環(huán)向應(yīng)力最大值在球殼內(nèi)壁，用壓力容器的公式可以算出該最大值在1000MPa左右[4]；深海挑戰(zhàn)者球殼選用的材料是高強(qiáng)鋼，其屈服應(yīng)力為1400MPa，所以該結(jié)構(gòu)的安全系數(shù)為1.4.從表1可以看出，當(dāng)前可選擇的耐壓殼材料只有高強(qiáng)度鋼和鈦合金兩種材料，限于材料科學(xué)的發(fā)展水平，現(xiàn)在并沒有可替代這兩種合金的材料.

基于現(xiàn)有材料，在材料性能達(dá)到極限的情況下，各國的深海潛水器耐壓殼結(jié)構(gòu)設(shè)計(jì)的安全系數(shù)依然很低.若耐壓殼設(shè)計(jì)僅從材料角度進(jìn)行發(fā)掘，則必然再次降低安全系數(shù)，以至于深潛器整體安全性難以保障.更為合理的是從結(jié)構(gòu)設(shè)計(jì)的角度，以增加殼體的厚度來降低應(yīng)力，這實(shí)際上也是深海挑戰(zhàn)者球殼設(shè)計(jì)的主要思路.但增加殼體的厚度會大大增加整個(gè)深潛器的重量.從表1可以看到，深海挑戰(zhàn)者球殼的密度已經(jīng)遠(yuǎn)遠(yuǎn)超過海水的密度了，這樣深潛器本身是沒法上浮的.這對耐壓浮力材料的研制提出了要求[5].作為球殼結(jié)構(gòu)，耐壓殼還可能在應(yīng)力遠(yuǎn)遠(yuǎn)沒有達(dá)到屈服應(yīng)力之前就發(fā)生屈曲失穩(wěn).雖然在后屈曲狀態(tài)，耐壓殼還能承受一定的壓力，但一旦屈曲失穩(wěn)，球殼結(jié)構(gòu)局部區(qū)域會出現(xiàn)較大應(yīng)變，導(dǎo)致塑性變形發(fā)生，使得球殼總體承受壓力的能力急劇降低.球殼結(jié)構(gòu)屈曲的一個(gè)重要特性就是對初始缺陷非常敏感，初始缺陷使得球殼結(jié)構(gòu)屈曲的臨界壓力急劇降低[6-11].除了屈曲失穩(wěn)，球殼還有突跳失穩(wěn)；某些情況下，發(fā)生突跳失穩(wěn)的壓力甚至還遠(yuǎn)遠(yuǎn)低于屈曲壓力[8].從制造的角度上講，要提高球殼的屈曲的應(yīng)力就要求非常高的加工精度.球殼加工精度的一個(gè)重要指標(biāo)就是真球度，即實(shí)測的曲率半徑和標(biāo)準(zhǔn)曲率半徑的比.以日本“深海6500”號深潛器和俄羅斯“和平2號”深潛器為例，真球度已經(jīng)非常高了，分別為[10]1.004和1.003.當(dāng)前單層球殼的耐壓倉結(jié)構(gòu)的設(shè)計(jì)已經(jīng)把材料性能和加工精度用到極致了，留給材料性能和加工精度的冗余空間是非常小的.

圖1 深潛器耐壓艙設(shè)計(jì)最小厚度與密度隨設(shè)計(jì)潛深變化關(guān)系Fig.1 The minimum thickness and density of the pressure hull as the functions of the submerging depth

圖1顯示了鈦合金(屈服強(qiáng)度900MPa，安全系數(shù)1.4)深潛器的單層殼體所需要的最小厚度以及結(jié)構(gòu)密度與下潛深度的關(guān)系.圖中方塊區(qū)域?yàn)閲鴥?nèi)外主要深潛器的殼壁厚度和潛水深度范圍.可以看出，隨著潛水深度的增加，所需要?dú)んw的最小厚度增加;但是同時(shí)，結(jié)構(gòu)密度也在增加，在萬米深度以下，結(jié)構(gòu)的密度將會超出1t/m3很多，使得海水浮力不足以支撐整體結(jié)構(gòu)的上浮.這又將會對深潛器浮力材料的設(shè)計(jì)提出新要求.

此外，傳統(tǒng)單層殼體結(jié)構(gòu)的設(shè)計(jì)，對材料的要求比較高.對于圖2所示的單層殼體結(jié)構(gòu)的受力模式，深潛器的主要受力狀態(tài)是深海中的高壓應(yīng)力，在深潛器的實(shí)際工作中，殼體內(nèi)部壓力Pin為0.1MPa(即1個(gè)大氣壓)，殼體外部壓力Pout取決于下潛深度.根據(jù)模型的受力狀態(tài)，殼體中徑向的應(yīng)力σr和環(huán)向的應(yīng)力σt為[4]

其中R為球殼內(nèi)徑，t為殼厚度.計(jì)算得到的Mises應(yīng)力按照殼體厚度分布的示意圖見圖2.最大應(yīng)力分布于殼體的內(nèi)表面，其他部分的材料遠(yuǎn)沒有達(dá)到屈服.這種單層殼體的設(shè)計(jì)，使得結(jié)構(gòu)中只有很小部分的材料承受了大應(yīng)力，大部分材料只承擔(dān)小應(yīng)力，這樣的設(shè)計(jì)勢必造成材料的浪費(fèi).因此，改進(jìn)結(jié)構(gòu)設(shè)計(jì)是目前唯一可大幅提高結(jié)構(gòu)耐壓性的途徑.

圖2 單層殼體的力學(xué)模型Fig.2 The free body diagram of one-layer shell

本文設(shè)計(jì)一種分層的球殼結(jié)構(gòu)，層間采用液壓油提供可調(diào)控的壓力，這樣球殼外壁的海水壓力就被分散到兩層或者多層殼上.分層結(jié)構(gòu)能調(diào)控層間液壓，優(yōu)化層間的壓力分布.壓力分布優(yōu)化的目標(biāo)是在滿足大容積的前提下，降低殼體內(nèi)的等效(最大)應(yīng)力,使得整體密度盡量小.層和層之間則通過桁架連接起來.這樣做的好處是：(1)內(nèi)層殼與外層殼通過桁架固定，不會相對旋轉(zhuǎn);(2)桁架連接兩層殼，使得兩層殼成為一個(gè)整體,能改變單層殼的高諧波屈曲模態(tài)[8]，大幅提高兩層殼的抗屈曲和抗突跳能力.本文的理論分析和數(shù)值模擬將顯示,在承受同樣海水壓力的情況下，與單層球殼相比，雙層球殼的結(jié)構(gòu)應(yīng)力和密度都有較大幅度的降低，抗屈曲能力顯著提高.這為將來的深潛器耐壓殼的設(shè)計(jì)提供了一個(gè)新思路.

1 雙層耐壓殼總體結(jié)構(gòu)和應(yīng)力分析

圖3為雙層耐壓殼總體結(jié)構(gòu)的示意圖.結(jié)構(gòu)分為內(nèi)外兩層，由密排的桁架相連；兩層殼體密封之后，層間充滿液壓油，外部的海水和液壓油通過一個(gè)可調(diào)壓的壓力平衡器連接(該壓力平衡器未出現(xiàn)在圖3中);上下兩個(gè)雙層球殼通過法蘭盤加固連接.該結(jié)構(gòu)的力學(xué)模型見圖4.

圖3 雙層耐壓殼總體結(jié)構(gòu)示意圖Fig.3 The schematic design of a two–layer pressure hull

圖4 雙層殼設(shè)計(jì)原理的力學(xué)模型Fig.4 Free body diagram of a two–layer pressure hull

使用式(1)和式(2)，可以得到內(nèi)外殼體沿殼厚度方向的應(yīng)力分布.計(jì)算中，取內(nèi)部壓力Pin為人體正常工作壓力 0.1MPa(近似為 1個(gè)大氣壓)，外部壓力Pout取100MPa，近似為10000m海深，R取為750mm.將應(yīng)力分布沿厚度做歸一化處理(歸一化厚度取為x/t,x為殼厚度方向的坐標(biāo),內(nèi)壁處x=0)，可以得到如圖5所示的內(nèi)外層殼體沿著厚度分布的應(yīng)力，包括單層和多層殼體在10000m海深環(huán)境下的應(yīng)力分布.采用單層殼，殼體厚度從75mm增加到100mm，可以使得最大Mises應(yīng)力從約600MPa降低到470MPa，降低了約22%；若采用雙層的殼體結(jié)構(gòu)，層間壓力取50MPa，層間厚度取10mm，最大應(yīng)力還能在這個(gè)基礎(chǔ)上再降低約20MPa.這里的雙層殼結(jié)構(gòu)的每層殼的厚度只有單層殼的一半.工藝上講，厚殼的加工難度高于薄殼.而且精度更難以控制.采用雙層殼體結(jié)構(gòu)的另一好處就是降低了對殼體加工的要求.因此，這種結(jié)構(gòu)不僅能夠有效緩解內(nèi)層壁面的應(yīng)力集中，而且還降低了殼體結(jié)構(gòu)對加工工藝的要求.另外，層間距離δ也是一個(gè)設(shè)計(jì)參數(shù)，可以通過調(diào)節(jié)δ的大小，來實(shí)現(xiàn)對結(jié)構(gòu)最大應(yīng)力的調(diào)節(jié)，見圖6.

圖5 單層和雙層(層間距離δ=10mm)殼體結(jié)構(gòu)應(yīng)力分布Fig.5 The comparison between the von Mises stress distributions of one-layer and two-layer pressure hulls(the two layer hull is with interlayer distance of δ=10mm)

圖6 殼體最大應(yīng)力與層間距離的關(guān)系Fig.6 The maximum von Mises stress versus the interlayer distance

2 生物體中的分層分壓結(jié)構(gòu)

這里以抹香鯨為例，簡述生物為適應(yīng)深海潛水而進(jìn)化出來的結(jié)構(gòu)，以對比分層分壓的結(jié)構(gòu)設(shè)計(jì)思想.抹香鯨的生物結(jié)構(gòu)遠(yuǎn)比雙層殼結(jié)構(gòu)復(fù)雜，但通過對比，能看到抹香鯨適應(yīng)深潛的結(jié)構(gòu)和要實(shí)現(xiàn)的核心功能是可以和雙層殼結(jié)構(gòu)一一對應(yīng)的.抹香鯨是潛水最深的動物，最大潛深2000～2200m，時(shí)間超過一小時(shí)[12-15]，下潛速度最快達(dá)到5m/s[12,14].抹香鯨中樞神經(jīng)所在的頭部在下潛過程中更需要保護(hù)，于是鯨頭進(jìn)化出了鯨蠟器官(spermaceti organ).鯨蠟 (spermaceti)是位于抹香鯨頭部的蠟質(zhì)，凝固點(diǎn)41?C～49?C.在抹香鯨下潛和捕獵的過程中，肌肉運(yùn)動發(fā)熱，體溫升高鯨蠟熔化變?yōu)橐后w，作用跟耐壓殼層間的液壓油一樣，可以更均勻地分布壓力.圖7橙色的部分為抹香鯨的脂肪組織(melon).在圖8中可以看到，脂肪組織中的脂肪(標(biāo)為f)里布滿了動脈(標(biāo)為a)和血竇(圖中標(biāo)為s)，造成了所謂的奇異網(wǎng)(rete mirabile)[13].奇異網(wǎng)的主要作用是收集氮?dú)?高壓下，血液中的氮?dú)馊芙舛壬?低壓下，血液中的氮?dú)馊芙舛冉档?，這樣上浮時(shí)，原先溶解的氮?dú)鈺纬蓺馀?，形成栓?這種氣栓障礙的形成就是潛水病根源[13].奇異網(wǎng)中的血管通過擴(kuò)散作用收集了溶解在血液里的氮?dú)猓獨(dú)庖徊糠秩芙庠谥纠?，剩下的就留在奇異網(wǎng)中的血管里，而奇異網(wǎng)血管里血液幾乎是不流動的[13].由于大量氮?dú)獗黄娈惥W(wǎng)血管收集，這就避免了在抹香鯨腦部血管形成氮?dú)鈿馀輀13].同時(shí)由于吸收氮?dú)獾木壒剩娈惥W(wǎng)中的血管和血竇的壓力是非常大的，巨大的壓力解釋了為什么奇異網(wǎng)中血竇尺寸很小的機(jī)理[13].抹香鯨調(diào)節(jié)血壓的另一個(gè)機(jī)制是二氧化碳.抹香鯨潛水時(shí)新陳代謝產(chǎn)生的二氧化碳并不排出體外，而用于調(diào)節(jié)血壓，并使得血管擴(kuò)張，這樣可以讓鯨蠟和脂肪組織在深海壓力下體積只有少量減少，保證抹香鯨上浮時(shí)有足夠浮力[15].動脈和血竇的作用相當(dāng)于耐壓殼的壓力平衡器，調(diào)節(jié)脂肪組織內(nèi)的壓力.鯨蠟和脂肪組織有兩重功能：分壓和讓壓力更均勻地分布.除了鯨蠟和脂肪組織之外，還有其他讓壓力均勻化的結(jié)構(gòu).如圖7所示，鯨蠟實(shí)際并不和鯨魚顱骨直接接觸，之間還有一個(gè)叫做前囊(frontal sac)的結(jié)構(gòu)，前囊貼合顱骨，囊內(nèi)充滿液體,可將鯨蠟的壓力再次均勻化傳遞到顱骨.抹香鯨頭部的結(jié)構(gòu)為皮膚/肌肉/鯨蠟/脂肪組織/前囊/顱骨的分層結(jié)構(gòu)，實(shí)現(xiàn)了分壓和壓力均勻化的功能，使得顱骨受到的力較小且相對均勻，這樣降低了顱骨總體所受應(yīng)力，避免了局部出現(xiàn)較大應(yīng)力，保障了顱內(nèi)中樞神經(jīng)的正常功能.總而言之，分層分壓結(jié)構(gòu)的核心功能就是將外載導(dǎo)致的壓力更加均勻地分布到各層，或者更準(zhǔn)確地說是降低應(yīng)力的非均勻分布，避免局部應(yīng)力過大而導(dǎo)致整體結(jié)構(gòu)失效.

圖7 抹香鯨頭部正面和側(cè)面的剖面圖[15]Fig.7 Anatomy of the sperm whale’s head with the views from the transverse and sagittal planes[15]

圖8 抹香鯨脂肪組織的微結(jié)構(gòu)電鏡圖.脂肪(f)中布滿了動脈(a)和血竇(s)形成的奇異網(wǎng)[13].奇異網(wǎng)里的動脈和血竇通過對氮?dú)夂投趸嫉氖占?，?shí)現(xiàn)了增壓和抹香鯨體積少量減少[13,15]Fig.8 Scanning electron microscopy representation of the thoracic rete,after freeze-cracking of the sample in liquid nitrogen,showing arteries(a)and sinusoids(s)embedded in fat(f)[13].The arteries and sinusoids are major components of rete mirabile with the function of collecting nitrogen and carbon dioxide.As a result,the pressure inside rete mirabile increases with a minor decrease of the sperm whale body volume[13,15]

3 結(jié)構(gòu)的屈曲分析

3.1 單層球殼屈曲

對于球殼結(jié)構(gòu)而言，可以根據(jù)厚度與半徑比將球殼分為薄殼和厚殼兩類.一般將厚度與半徑比小于0.05的殼體定義為薄殼，大于0.05的則為厚殼[16].相對于薄殼理論/模型而言，厚殼理論要復(fù)雜得多，目前對厚球殼應(yīng)力尚無統(tǒng)一的理論解[16-17].劉濤[17]考慮徑向、環(huán)向應(yīng)力形成的應(yīng)力強(qiáng)度和材料的非線性，推導(dǎo)了厚殼的非彈性失穩(wěn)壓力的簡化計(jì)算公式.從表1可以看出，除美國“深海挑戰(zhàn)者”號的耐壓球殼外，其余的所有耐壓球殼厚度與半徑比均小于0.05，薄殼理論適用，所以這里僅重點(diǎn)介紹薄殼理論.早期經(jīng)典的球殼屈曲分析在鐵木辛柯專著[18]中有很好的綜述,對于厚殼的情形，可參考文獻(xiàn)[19-20]的分析.對于薄殼模型，存在平截面假定：即垂直于中面的線元在殼變形后仍然和中面垂直并且保持直線.該平截面假定導(dǎo)致的直接結(jié)果就是薄殼變形的應(yīng)變能只由兩部分組成：彎曲和中面拉伸的應(yīng)變能[18].基于該平截面假定，鐵木辛柯推導(dǎo)出單層球殼屈曲的臨界壓力為[18]

式中，E,t,R和v分別為球殼的楊氏模量、厚度、半徑和泊松比.忽略上式的第2項(xiàng)，式(3)變?yōu)?/p>

式(4)即為Zoelly[21]在1915年和Schwerin[22]在1922年分別獨(dú)立得到的單層球殼屈曲臨界壓力.以上兩個(gè)臨界壓力公式的得出都假設(shè)了屈曲模態(tài)對稱[18].但對稱結(jié)構(gòu)在均勻?qū)ΨQ的外載作用下是可以發(fā)生非對稱屈曲的.更具一般性的、考慮對稱和非對稱屈曲模態(tài)的解由Van der Neut[23]在1932年給出.通常認(rèn)為均勻?qū)ΨQ、無幾何缺陷的球殼在均勻壓力作用下是不會發(fā)生非對稱屈曲的[18,24].所以，假設(shè)對稱屈曲模態(tài)的式(3)和式(4)被認(rèn)為是可以用于(無幾何缺陷)球殼屈曲壓力的計(jì)算[18];但需要強(qiáng)調(diào)的是一旦有幾何缺陷，球殼在均勻壓力作用下也可以發(fā)生非對稱屈曲[25-26].

式(3)和式(4)的推導(dǎo)過程包含如下的假設(shè)：擾動無窮小，材料均質(zhì)，球殼無幾何缺陷,應(yīng)力與應(yīng)變滿足線性關(guān)系等.實(shí)際上這些假設(shè)很難一一滿足,式(3)和式(4)的預(yù)測值遠(yuǎn)遠(yuǎn)高于實(shí)驗(yàn)測量值.早期殼體屈曲研究的一個(gè)重要內(nèi)容就是要解釋經(jīng)典理論和實(shí)驗(yàn)測量值之間存在巨大差距的問題.在這方面，馮·卡門、錢學(xué)森和他們加州理工的同事們做出了系統(tǒng)性和代表性的工作[27-31].他們發(fā)展了大撓度非線性和有限撓度擾動的殼屈曲理論來解釋理論和實(shí)驗(yàn)之間存在巨大差別的問題.當(dāng)時(shí)加州理工大學(xué)的水池實(shí)驗(yàn)顯示球殼屈曲的壓力只有式(4)預(yù)測的1/4左右[27,31].馮·卡門等[28]發(fā)現(xiàn)基于小撓度變形的線性理論總是會過高估計(jì)臨界屈曲壓力；同時(shí)他們也發(fā)現(xiàn),對于薄殼而言，即便材料處于彈性范圍，其外載與撓度的關(guān)系也是強(qiáng)非線性的[27-29].式(4)對應(yīng)的是線性和無窮小撓度擾動的理論[27].這里的無窮小撓度擾動有兩個(gè)涵義：(1)結(jié)構(gòu)是完美、無幾何缺陷的；(2)其所受擾動是無窮小的.這樣的假設(shè)是為了保證結(jié)構(gòu)剛發(fā)生屈曲時(shí)的撓度為零，這也是我們通常遇到的屈曲的情形.馮·卡門等[27]的大撓度非線性和有限撓度擾動的理論顯示球殼的屈曲壓力有一大一小兩個(gè)：大的屈曲壓力對應(yīng)的是無窮小撓度擾動的情形，也就是式(4)對應(yīng)的情形.他們的分析也指出這種通過無窮小撓度擾動達(dá)到的大屈曲壓力的情形需要對樣品加工和實(shí)驗(yàn)過程提出極其苛刻的要求才能實(shí)現(xiàn).小的屈曲壓力對應(yīng)的是有限撓度擾動的情形，即球殼屈曲時(shí)已經(jīng)通過突跳獲得了一個(gè)有限撓度，物理上該擾動可以通過外載變化或者結(jié)構(gòu)的動態(tài)效應(yīng)實(shí)現(xiàn).而且小的屈曲壓力對應(yīng)的就是實(shí)際應(yīng)用的情形，真實(shí)結(jié)構(gòu)“傾向于”在較小壓力下發(fā)生有限撓度屈曲[27].他們的計(jì)算結(jié)果顯示該屈曲壓力大約為式(4)預(yù)測值的30%，與實(shí)驗(yàn)數(shù)據(jù)的25%基本相符.這樣就解釋了經(jīng)典線性無窮小撓度擾動理論和實(shí)驗(yàn)數(shù)據(jù)之間巨大差距的問題.這里需要指出的是這種大撓度非線性和有限撓度擾動的理論是不依賴于初始缺陷的[27]，在不包含初始缺陷的情形下，擾動就使屈曲壓力降低了70%.前面提到的初始缺陷也是球殼結(jié)構(gòu)屈曲壓力急劇降低[6-11]的最重要的因素之一.馮·卡門和錢學(xué)森的工作[28-31]也考慮了這一重要因素.他們的另一重要貢獻(xiàn)是研究了殼體的另一種失穩(wěn)：突跳[28-30].拱或殼結(jié)構(gòu)在壓力作用下,其撓度并不一直隨壓力單調(diào)增加，在某個(gè)臨界壓力值(即壓力隨撓度開始下降時(shí)的值)，結(jié)構(gòu)會經(jīng)歷一個(gè)突然的變化，構(gòu)型翻轉(zhuǎn)，跳躍到新的平衡態(tài)[32].雖然都是失穩(wěn)，但屈曲和突跳對應(yīng)的失穩(wěn)類型是不一樣的：屈曲是叉式分岔失穩(wěn)，而突跳是鞍點(diǎn)分岔失穩(wěn).對于薄殼而言，突跳發(fā)生的壓力可能遠(yuǎn)遠(yuǎn)小于屈曲壓力[8].因?yàn)橥惶殡S著大的位移變化，突跳失穩(wěn)一旦發(fā)生極可能導(dǎo)致整個(gè)殼結(jié)構(gòu)失效.錢學(xué)森還發(fā)展了側(cè)面掛有非線性彈簧的桿受壓失穩(wěn)的模型并給出了近似解析解[31].用側(cè)面掛有非線性彈簧的受壓桿模型來研究殼的屈曲最早是由Cox提出的[33],這個(gè)模型在殼屈曲的研究中有重要意義.由于殼屈曲的控制方程非常復(fù)雜，各量的作用并不清晰;而這個(gè)模型，特別是側(cè)面懸掛的彈簧模型，非常直接地給出殼非對稱屈曲的物理機(jī)制：看似對稱的結(jié)構(gòu)實(shí)際上是不對稱的.殼非對稱屈曲的一個(gè)直接后果是在變形區(qū)域(有限撓度擾動情形下)，非對稱屈曲的壓力會小于對稱屈曲的壓力[31].非對稱屈曲的臨界壓力由桿本身的彎曲剛度和側(cè)掛非線性彈簧的剛度共同決定，而彈簧的剛度受結(jié)構(gòu)的初始幾何缺陷影響[31].

馮·卡門和錢學(xué)森的殼體失穩(wěn)研究針對的是航空航天器的結(jié)構(gòu),如火箭、飛機(jī)的殼體.而對深潛器耐壓結(jié)構(gòu)研究有非常大影響的是20世紀(jì)60年代美國海軍水面作戰(zhàn)中心的以Krenzke為首的先進(jìn)結(jié)構(gòu)小組[25-26].他們以深潛系統(tǒng)耐壓艇體為對象，對多個(gè)球殼模型進(jìn)行了試驗(yàn)研究.Krenzke等進(jìn)行實(shí)驗(yàn)的場所叫做大衛(wèi)·泰勒模型水池(David Taylor Model Basin).Krenzke等從殼體失穩(wěn)實(shí)驗(yàn)得出的經(jīng)驗(yàn)公式在國內(nèi)通常稱為泰勒水池公式[16-17].泰勒水池公式如下[17]

這里Es,Et分別為材料的正割楊氏模量和正切楊氏模量，也就是說對于塑性情形，等效楊氏模量換為泰勒水池公式還有如下形式[16]

式中，Cz為制造效應(yīng)影響系數(shù).我國潛水器規(guī)范整球殼屈曲壓力計(jì)算公式為[16,34]

這里C為半徑修正系數(shù)，可通過比值t/R查曲線確定；Cs為材料物理非線性修正系數(shù).

各國潛水器耐壓殼體設(shè)計(jì)主要是參照潛艇的設(shè)計(jì)標(biāo)準(zhǔn)，但由于結(jié)構(gòu)形式、使用條件和材料與潛艇不同,要求亦有所不同.各國規(guī)范有很大差別，很多國家甚至沒有規(guī)范，比如歐盟的結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)ENV 1993-1-6[35]對圓殼屈曲就沒有設(shè)計(jì)規(guī)則[8].直接后果就是現(xiàn)今很多船級社均有耐壓殼體設(shè)計(jì)標(biāo)準(zhǔn)，但互相之間差別很大[16].現(xiàn)今，還有研究者不斷地提出類似式(5)～式(7)的屈曲壓力公式來作殼體設(shè)計(jì)的一個(gè)指標(biāo).比如針對鈦合金材料球殼，國內(nèi)就有數(shù)種不同的屈曲壓力公式[36-38].文獻(xiàn)[16]對這些屈曲壓力公式有較好的綜述和歸納,這些屈曲壓力公式都是考慮了初始缺陷的擬合公式.與式(5)～式(7)相比，這些新近提出的屈曲壓力公式都比較復(fù)雜，這里就不一一列出.我國潛水器規(guī)范的式(7)與式(5)和式(6)相比有一個(gè)共同的系數(shù)0.84，可見泰勒水池公式實(shí)際是我國潛水器規(guī)范的基礎(chǔ).我們這里著重討論一下泰勒水池公式以及以此為標(biāo)準(zhǔn)可能存在的問題.

泰勒水池公式實(shí)際是以式 (4)為基礎(chǔ)得出的.Krenzke和Kiernan[25]的17個(gè)“近乎完美”圓殼試件的系列實(shí)驗(yàn)顯示這些“近乎完美”圓殼的屈曲壓力最多只能達(dá)到式(4)預(yù)測值的70%.Krenzke和Kiernan“近乎完美”的定義包含兩個(gè)方面：(1)加工圓殼產(chǎn)生的殘余應(yīng)力小到可以忽略；(2)圓殼的幾何缺陷小到幾乎無法測量[25].取泊松比v=0.3[25],壓力取70%的上限值，式(4)變?yōu)?/p>

也就是說式(5)～式(7)中的0.84系數(shù)實(shí)際上接近上限值，這就存在過高估計(jì)屈曲壓力的問題.另外，式(6)和式(7)中的系數(shù)Cz和C的物理意義是缺陷的影響，問題更大，因?yàn)閷?shí)際中遇到的缺陷模式各種各樣，很難精確標(biāo)度不同缺陷對殼屈曲壓力的影響.理論、計(jì)算與實(shí)驗(yàn)三者之間差距巨大，實(shí)驗(yàn)的不可重復(fù)性更是困擾殼屈曲研究的瓶頸問題.通過對“近乎完美”圓殼的研究，Krenzke和Kiernan實(shí)際上避開了這個(gè)至今仍困擾殼屈曲研究者的問題.

3.2 桁架連接的雙層耐壓殼屈曲分析

對單層殼的分析表明，球殼這么一個(gè)“簡單”的結(jié)構(gòu)其實(shí)很復(fù)雜：非線性導(dǎo)致對稱和非對稱屈曲模態(tài)并存,突跳、有限擾動和缺陷導(dǎo)致屈曲壓力下降，以及實(shí)驗(yàn)不可重復(fù).為了提高深潛器耐壓球殼的屈曲壓力，加工精度已經(jīng)被用到極致來消除幾何缺陷，如前所述耐壓球殼的真球度已經(jīng)達(dá)到1.003[10].但真球度實(shí)質(zhì)上還是一個(gè)總體性衡量幾何缺陷的指標(biāo)，還有其他的幾何缺陷，比如球殼厚度的均勻性.一個(gè)總的趨勢是球殼幾何缺陷對屈曲壓力的影響隨R/t(半徑與厚度比)增加而增大：pe/pcr隨R/t增大而單調(diào)遞減[8,39](pe為實(shí)驗(yàn)測量的屈曲壓力，pcr為式(4)給出的屈曲壓力).沒有準(zhǔn)確標(biāo)度缺陷影響的模型以及實(shí)驗(yàn)不可重復(fù)性使得現(xiàn)今的耐壓球殼的設(shè)計(jì)嚴(yán)重依賴于式(5)～式(7)之類的經(jīng)驗(yàn)公式.需要強(qiáng)調(diào)的是式(5)～式(7)是基于“近乎完美”的圓殼實(shí)驗(yàn)數(shù)據(jù)建立的，而且這些“近乎完美”的圓殼屈曲壓力值是式(4)預(yù)測值的42%～70%,波動幅度；缺陷的存在甚至讓屈曲壓力值降到預(yù)測值的20%以下，使實(shí)驗(yàn)結(jié)果分散性大[25].除了幾何缺陷，還有來自材料本身的缺陷，更有加工過程產(chǎn)生的殘余應(yīng)力、空隙、裂紋等等.雖然這些因素在式(6)和式(7)中有所考慮(比如系數(shù)Cz),但坦率地說,這些經(jīng)驗(yàn)公式是沒有堅(jiān)實(shí)的理論和實(shí)驗(yàn)數(shù)據(jù)作支撐的.這也是深潛器單層耐壓球殼設(shè)計(jì)和使用的最大問題.例如，Samuelson和Eggwertz[39]給出的作為圓殼安全設(shè)計(jì)的屈曲壓力極限曲線.該極限曲線給出的是有初始缺陷的圓殼屈曲壓力的下限.鋁制圓殼的實(shí)驗(yàn)數(shù)據(jù)均在該極限曲線之上，用該極限曲線設(shè)計(jì)鋁制圓殼是安全的；但鋼制圓殼的實(shí)驗(yàn)數(shù)據(jù)則是另一番圖景了：有相當(dāng)多的鋼制圓殼的屈曲壓力低于該極限曲線[8]，因此用該極限曲線來設(shè)計(jì)鋼制圓殼是不安全的.另一個(gè)要強(qiáng)調(diào)的事實(shí)是式(5)～式(7)只考慮了屈曲，沒有考慮突跳.Wunderlich和Albertin對完美無幾何缺陷的圓殼的計(jì)算顯示，在無塑性變形的情形下，完美圓殼的屈曲壓力大約為式(4)預(yù)測值的80%，如果發(fā)生塑性變形的話，完美圓殼在壓力僅為式(4)預(yù)測值的9%時(shí)就會發(fā)生突跳[8].桁架連接的雙層耐壓殼結(jié)構(gòu)設(shè)計(jì)的目的在于降低缺陷特別是幾何缺陷對球殼屈曲壓力以及突跳壓力的影響.

圖9給出桁架連接的雙層耐壓殼結(jié)構(gòu)的示意圖.內(nèi)外層球殼之間用密排的桁架連接.桁架密排設(shè)計(jì)的主要目的就是提高球殼屈曲壓力和突跳壓力.(無缺陷)單層球殼屈曲臨界壓力正如式(3)～式(8)給出的一樣：幾何上是由t/R的比值決定的.屈曲模態(tài)也同樣由t/R(或者R/t)決定[8].R/t越大，球殼屈曲的波數(shù)越大，比如R/t=50時(shí)，球殼是以8個(gè)傅里葉諧波模態(tài)屈曲的；R/t=1000時(shí)，球殼是以41個(gè)傅里葉諧波模態(tài)屈曲的[8].圖10給出了R/t=1000時(shí)的單層球殼的41個(gè)傅里葉諧波模態(tài)的截面圖.對于薄殼結(jié)構(gòu)，根據(jù)定義R/t>200，對于大的R/t值，屈曲的波數(shù)就很大.密排桁架的介入相當(dāng)于加入了約束,阻止此模態(tài)的屈曲發(fā)生.如果桁架都位于反節(jié)點(diǎn),即上述屈曲諧波變形的波峰和波谷的位置時(shí)，向內(nèi)變形的諧波因?yàn)殍旒芗s束，波幅就會被大大抑制住，而不是像圖10那樣內(nèi)外對稱的，相應(yīng)的屈曲壓力就提高了.對于二維(無桁架支撐的)結(jié)構(gòu)，一些高(共振)頻率也會引發(fā)如圖10所示的大波數(shù)變形[40].密排桁架的加入從結(jié)構(gòu)力學(xué)上講相當(dāng)于給殼結(jié)構(gòu)增加了剛度，有效抵抗大波數(shù)的變形.下面對桁架連接的雙層耐壓殼結(jié)構(gòu)做屈曲分析.

圖9 桁架連接的雙層耐壓殼剖面和帶桁架內(nèi)層殼示意圖Fig.9 The cross-section of two-layer pressure hull connected by trusses and a schematic of the inner shell with trusses

圖10 R/t=1000的球殼41個(gè)傅里葉諧波屈曲模態(tài)[8]Fig.10 Buckling mode of a shell with R/t=1000 and 41 Fourier harmonics[8]

如圖11所示,由于密排桁架連接的緣故，球殼實(shí)際被分割為多個(gè)小的亞結(jié)構(gòu).該亞結(jié)構(gòu)為柱殼，a為柱殼(弧)邊長；R為球殼半徑；t為球殼厚度.球殼外的水壓是pw,殼內(nèi)壓力為po；pw和po的方向?yàn)闅で娣ㄏ?

只考慮薄殼的情形，即不考慮殼徑向的應(yīng)力，只考慮環(huán)向的應(yīng)力.環(huán)向應(yīng)力p和面外壓力pw和po有如下關(guān)系[4]

為了和式 (4)比較，將po設(shè)為零.因?yàn)楸そY(jié)構(gòu)t?R，這樣式(9)可以近似為

引入如下兩個(gè)無量綱數(shù)[41]

這里，E*=E/(1?ν)為雙軸模量；其中E為楊氏模量，ν為泊松比.臨界的無量綱屈曲壓力為[41]

因?yàn)槭?12)的第二項(xiàng)是一個(gè)小量，可以忽略，這樣式(12)就近似為

結(jié)合式(10)～式(12)，得到柱殼亞結(jié)構(gòu)的屈曲壓力為

與式(4)的屈曲壓力相比，得到

很明顯，桁架連接的雙層殼中的單個(gè)殼要比單層殼抗屈曲性能高的指標(biāo)就是pw/pcr>1,從式(15)可以看出，這完全是由幾何因素決定的：當(dāng)ν=0.3時(shí)，pw/pcr>1對應(yīng)的幾何條件就是這也就是選擇密排桁架的原因：a是由桁架數(shù)量決定的，數(shù)量越多，a越小.增加桁架數(shù)量是提高球殼抗屈曲能力的辦法.Krenzke和Kiernan[25]用實(shí)驗(yàn)方法系統(tǒng)地研究過圓殼亞結(jié)構(gòu)對圓殼整體抗屈曲能力的影響，他們定義了如下的無量綱數(shù)

Krenzke和Kiernan[25]發(fā)現(xiàn)當(dāng)θ>2.5時(shí)，將圓殼分割成亞結(jié)構(gòu)不但沒有增加反而降低了圓殼屈曲壓力.也就是說，只有當(dāng)θ<2.5時(shí)才能增加殼屈曲壓力.當(dāng)ν=0.3時(shí)，那么θ<2.5對應(yīng)的就是這和前面得到的是基本一致的.在公式推導(dǎo)的過程中，式(12)是按柱殼結(jié)構(gòu)四邊簡支計(jì)算，四邊簡支結(jié)構(gòu)的屈曲壓力是最小的，實(shí)際的結(jié)構(gòu)介于簡支和固支之間.但同時(shí)，式(12)沒有考慮幾何缺陷，這兩個(gè)互相競爭/抵消的因素造成值和Krenzke和Kiernan實(shí)驗(yàn)值[25]有6%的差別.理論上講，如果桁架的密排程度達(dá)到時(shí)，那么pw/pcr=2,也就是說這時(shí)桁架連接的雙層殼結(jié)構(gòu)的單殼與同樣尺寸、同樣材料的單層球殼結(jié)構(gòu)相比，抗屈曲能力提高了一倍,雖然實(shí)際加工過程中，增加桁架數(shù)量會增加焊接點(diǎn)，會引入一些缺陷從而降低結(jié)構(gòu)的屈曲壓力，也會增加結(jié)構(gòu)安裝難度.在理論分析層面上講，增加桁架的密排程度是有效提高雙層殼結(jié)構(gòu)屈曲壓力的有效途徑，同時(shí)Krenzke和Kiernan實(shí)驗(yàn)也是支持該方法的.

桁架連接的雙層殼結(jié)構(gòu)還有一個(gè)好處就是降低了甚至還可能消除了殼體突跳失穩(wěn).在一些情況下，突跳在遠(yuǎn)小于結(jié)構(gòu)屈曲壓力時(shí)就發(fā)生[8].突跳和屈曲是兩種不同的失穩(wěn)，有4種情形：(1)屈曲發(fā)生而突跳未發(fā)生；(2)突跳發(fā)生而屈曲未發(fā)生；(3)突跳后引發(fā)屈曲或屈曲后引發(fā)突跳；(4)屈曲和突跳同時(shí)發(fā)生.如果這兩種失穩(wěn)的壓力不一樣，就是(1)和(2)的情形;但(3)和(4)的情形目前尚無有效方法確定.實(shí)驗(yàn)上，區(qū)分情形(1)和(2)極其困難.Krenzke和Kiernan的球殼實(shí)驗(yàn)[25]就發(fā)生過各種屈曲理論都解釋不了的現(xiàn)象，這應(yīng)該是突跳和屈曲都發(fā)生了的情形.Krenzke和Kiernan多次使用的一個(gè)詞就是“崩塌”(collapse)[25],對于圓殼結(jié)構(gòu)而言，“崩塌”可指上述四種情形的任何一種.突跳的發(fā)生不但和結(jié)構(gòu)有關(guān)還和載荷有關(guān)[32].在同一載荷條件下，“深”(deep)拱和殼比“淺”(shallow)拱和殼要更容易發(fā)生突跳，“淺”拱和殼在某些載荷條件下甚至不發(fā)生突跳,而“深”拱和殼發(fā)生突跳[32].“深”和“淺”是由拱/殼的高度和跨度的比值決定的，比值大的為“深”,比值小的為“淺”[32].如圖11所示，“深”的球殼因?yàn)楸昏旒芊指钭兂闪恕皽\”的柱殼了.密排桁架連接在提高抗屈曲能力的同時(shí)，也可以消除突跳失穩(wěn).

同樣，桁架連接的雙層耐壓殼體的設(shè)計(jì)也是受到了生物結(jié)構(gòu)啟發(fā)的：桁架連接的雙層耐壓殼體的設(shè)計(jì)即與鸚鵡螺的結(jié)構(gòu)有異曲同工之處.鸚鵡螺的螺殼的剖面圖見圖12.整個(gè)螺殼被隔片分割為大約30個(gè)腔室,腔室之間由連室細(xì)管連接.連室細(xì)管是頭足綱目動物特有的長條狀組織，通過調(diào)節(jié)海水中的鈉、氯離子濃度產(chǎn)生滲透壓梯度，讓腔室里的水進(jìn)出，從而調(diào)節(jié)鸚鵡螺的浮力[42].鸚鵡螺下沉?xí)r，腔室進(jìn)水但不是充滿的，腔室內(nèi)的壓力還稍微低于一個(gè)大氣壓[42]，所以鸚鵡螺的螺殼和我們的耐壓殼一樣面臨著深潛時(shí)海水巨大壓力的問題.從式(14)可以看出腔室殼的(特征)尺寸(a)越小，那么殼抗屈曲的能力就越大.鸚鵡螺最外的腔室最大，是鸚鵡螺軀體都在這個(gè)腔室，而且該腔室是直接與海水相通的，所以不存在抗壓力的問題.而其他腔室是與海水隔絕的，存在抗壓的問題.所以和最大腔室相比，其他腔室都小得多.在活體鸚鵡螺的壓裂實(shí)驗(yàn)中，總是最外層的隔片首先失效,因?yàn)樗某叽缱畲骩42].隔片在這里的作用和桁架一樣：分割結(jié)構(gòu)，減小結(jié)構(gòu)尺寸，從而提高結(jié)構(gòu)抗屈曲能力.鸚鵡螺還面臨一個(gè)材料上的挑戰(zhàn)：螺殼主要是由碳酸鈣組成的.與鈦合金和高強(qiáng)度鋼相比，碳酸鈣楊氏模量較小而且基本上沒有延性.鸚鵡螺螺殼材料的微結(jié)構(gòu)類似桁架.圖13(a)是鸚鵡螺殼材料微結(jié)構(gòu)的共聚焦顯微鏡圖像，層間有柱狀結(jié)構(gòu)支撐；圖13(b)是針對圖13(a)提出的正方形網(wǎng)格模型[43].同樣的這些柱狀結(jié)構(gòu)跟我們的桁架一樣在材料的微結(jié)構(gòu)層次對結(jié)構(gòu)分割.根據(jù)正方形網(wǎng)格模型和鸚鵡螺的活體實(shí)驗(yàn)，Sherrard給出了鸚鵡螺殼體材料能下潛的最大深度為[43]

這里D是下潛深度,單位為m，t是微結(jié)構(gòu)層的厚度，r是柱狀支撐物間隔的一半[43].式(17)和式(15)具有相似性：結(jié)構(gòu)和材料的抗壓能力和間隔(a或r)平方成反比.實(shí)驗(yàn)測量的數(shù)據(jù)顯示深海鸚鵡螺的r/t在3.81到5.96之間，平均值為[43]5.43.將上述r/t的值代入式(17)，可算出鸚鵡螺的最大潛深在617m至1508m之間，平均值為743m，這與海洋學(xué)家觀察到的最大潛深800m是大致相符的[42].

圖12 鸚鵡螺殼的剖面圖[42]Fig.12 Cross-section of a nautilus shell[42]

圖13 鸚鵡螺殼微結(jié)構(gòu)Fig.13 The microstructure of a nantilus shell

4 結(jié)論和展望

桁架連接的分層殼的耐壓性能通過兩重手段來提高：(1)降低每層殼所受外載,(2)提高每層殼的抗屈曲能力.分層結(jié)構(gòu)是為了實(shí)現(xiàn)分壓的功能，層間可調(diào)控的液壓油壓將海水壓力傳遞到內(nèi)層，分別由內(nèi)外球殼來承擔(dān).密排桁架將球殼分解為多個(gè)亞結(jié)構(gòu)，以提高結(jié)構(gòu)整體抗屈曲和突跳的能力.對桁架連接的分層結(jié)構(gòu)開展的理論分析和初步的數(shù)值模擬表明,該結(jié)構(gòu)與現(xiàn)今的單層耐壓殼結(jié)構(gòu)相比有著巨大的抗壓優(yōu)勢，是耐壓殼結(jié)構(gòu)設(shè)計(jì)的一個(gè)新方向,擺脫了對材料和加工精度的依賴.本文還從結(jié)構(gòu)力學(xué)的角度推導(dǎo)了從實(shí)驗(yàn)數(shù)據(jù)總結(jié)出來的泰勒水池經(jīng)驗(yàn)公式，使得深潛器殼體的設(shè)計(jì)有了更加堅(jiān)實(shí)的力學(xué)基礎(chǔ).

1 閆潔.院士解讀：海洋研究為何需要超級計(jì)算機(jī)，科學(xué)網(wǎng)2016年8月23日報(bào)道(http://news.sciencenet.cn/htmlnews/2016/8/354387.shtm)

2 維基百科：https://en.wikipedia.org/wiki/Bathysphere

3 維基百科：https://en.wikipedia.org/wiki/Bathyscaphe

4 Cook RD,Young WC.Advanced Mechanics of Materials,2nd edn.Upper Saddle River,NJ:Prentice Hall,1999

5 維基百科：https://en.Wikipedia.org/wiki/Deepsea_Challenger

6 Bushnell D.Buckling of shells—pitfall for designers.AIAA Journal,1981,19:1183-1226

7 Blachcut J,Galletly GD.Buckling strength of imperfect steel hemispheres.Thin-Walled Structures,1995,23:1-20

8 Wunderlich W,Albertin U.Buckling behaviour of imperfect spherical shells.International Journal of Non-Linear Mechanics,2002,37:589-604

9 Pan BB,Cui WC.An overview of buckling and ultimate strength of spherical pressure hull under external pressure.Marine Structures,2010,23:227-240

10 劉濤.深海載人潛水器耐壓球殼設(shè)計(jì)特性分析.船舶力學(xué)，2007,11(2):214-22(LiuTao.Researchonthedesignofsphericalpressure hull in manned deep-sea submersible.Journal of Ship Mechanics,2007,11(2):214-222(in Chinese))

11 Teng J.Buckling of thin shells:Recent advances and trends.Applied Mechanics Review,1996,49:263-274

12 Tyack PL,Johnson M,Soto NA,et al.Extreme diving of beaked whales.Journal of Experimental Biology,2006,209:4238-4253

13 Blix A,Wallφe L,Messelt EB.On how whales avoid decompression sickness and why they sometimes strand.Journal of Experimental Biology,2013,216:3385-3387

14 Goldbogen JA,Calambokidis J,Shadwick RE,et al.Kinematics of foraging dives and lunge-feeding in fi n whales.Journal of Experimental Biology,2006,209:1231-1244

15 維基百科：https://en.wikipedia.org/wiki/Sperm_whale

16 李文躍，王帥，劉濤等.大深度載人潛水器耐壓殼結(jié)構(gòu)研究現(xiàn)狀及最新進(jìn)展.中國造船,2016,57(1):210-221(Li Wenyue,Wang Shuai,Liu Tao,et al.Current status and progress on pressure hull structure of manned deep submersible.Shipbuilding of China,2016,57(1):210-221(in Chinese))

17 劉濤.大深度潛水器耐壓殼體彈塑性穩(wěn)定性簡易計(jì)算方法.中國造船,2000,42(3):8-14(Liu Tao.Simpli fi ed method of inelastic buckling analysis of deep submersible pressure hull.Shipbuilding of China,2000,42(3):8-14(in Chinese))

18 Timoshenko S.Theory of Elastic Stability,Engineering Societies Monographs.New York:McGraw Hill,1936

19 殷有泉,陳朝偉.軟化材料厚壁筒的解析解及其穩(wěn)定性分析.力學(xué)學(xué)報(bào),2010,42(1):56-64(Yin Youquan,Chen Chaowei.The analytical solutions of thick-walled cylinder of softening material and its stability.Chinese Journal of Theoretical and Applied Mechanics,2010,42(1):56-64(in Chinese))

20 趙均海,姜志琳,張常光等.不同拉壓特性的厚壁圓筒極限內(nèi)壓統(tǒng)一解.力學(xué)學(xué)報(bào),2017,49(4):836-847(Zhao Junhai,Jiang Zhilin,Zhang Changguang,et al.Uni fi ed solutions of limit internal pressureforthick-walledcylinderwithdi ff erentbehaviourintensionand compression.Chinese Journal of Theoretical and Applied Mechanics,2017,49(4):836-847(in Chinese))

21 Zoelly R.überein Knickungsproblem am der Kugelschale,Dissertation,Zrich,1915

22 Scherwin E.Zur Stabilitat der dunnwan digen Hohlkugel unter gleichmaβigem Auβendruck.,1922,2:81-91

23 Van der Neut A.De Elastische Stabiliteit van den Dunwandigen Bol,H.J.Paris,Amsterdam,1932

25 Krenzke MA,Kiernan TJ.Elastic stability of near-perfect shallow spherical shells.AIAA Journal,1963,1:2855-2857

26 Krenzke MA.The elastic buckling strength of near-perfect deep spherical shells with ideal boundaries.David Taylor Model Basin Report,1963:1713

27 Von Karman T,Tsien HS.The buckling of spherical shells by external pressure.Journal of the Aeronautical Sciences,1939,7:43-50

28 Von Karman T,Dunn LG,Tsien HS.The In fl uence of curvature on the buckling characteristics of structures.Journal of the Aeronautical Sciences,1940(7):276-289

29 Von Karman T,Tsien HS.The buckling of thin cylindrical shells under axial compression.Journal of the Aeronautical Sciences,1941,8:303-312

30 Tsien HS.A theory for the buckling of thin shells.Journal of the Aeronautical Sciences,1942,9:373-384

31 Tsien HS.Buckling of a column with non-linear lateral supports.Journal of the Aeronautical Sciences,1942,9:119-132

32 Zhang Y,Wang YS,Li ZH,et al.Snap-through and pull-in instabilities of an arch-shaped beam under electrostatic loading.Journal of Microelectromechanical Systems,2007,16:684-693

33 Cox HL.Stress analysis of thin metal construction.Journal of the Royal Aeronautical Society,1940,44:231-282

34 中國船級社.潛水系統(tǒng)和潛水器入級與建造規(guī)范.北京：人民交通出版社,2013(China Classi fi cation Society:The construction and classi fi cation rules for the submersible systems and vehicles.Beijing:China Communications Press，2013(in Chinese))

35 Eurocode 3,ENV 1993-1-6,Design of Steel Structures—Parts 1–6:General Rules—Supplementary Rules for the Shell Structures,CEN,Brussels,1999

36 石佳睿,唐文勇.載人深潛器鈦合金耐壓球殼極限強(qiáng)度可靠性分析.船海工程,2014,43(2):114-118(Shi Jiarui,Tang Wenyong.Ultimate strength reliability analysis of Titanium alloy spherical pressure shell in HOV.Ship&Ocean Engineering,2014,43(2):114-118(in Chinese))

37 王丹,萬正權(quán).考慮初始形狀影響的耐壓球殼臨界載荷簡化計(jì)算公式.船舶力學(xué),2014,18(5):557-564(Wang Dan,Wan Zhengquan.Simpli fi edcalculationformulaforcriticalloadofspherical pressure shell considering the in fl uence of initial shape.Journal of Ship Mechanics,2014,18(5):557-564(in Chinese))

38 Pan BB,Cui WC,Shen YS,et al.Further study on the ultimate strength analysis of spherical pressure hulls.Marine Structures,2010,23:444-461

39 Samuelson LA,Eggwertz S.(Eds.)Shell Stability Handbook.London:Elsevier,1992

40 方奕忠,王鋼,沈韓等.環(huán)形薄板二維駐波的研究.力學(xué)學(xué)報(bào),2015,47(4):664-671(Fang Yizhong,Wang Gang,Shen Han,et al.The research of 2-dimensional standing waves on annular plate.ChineseJournalofTheoreticalandAppliedMechanics,2015,47(4):664-671(in Chinese))

41 Chajes A.Principles of Structural Stability Theory.Upper Saddle River,NJ:Prentice Hall,1974

42 維基百科：https://en.wikipedia.org/wiki/Nautilus

43 Sherrard K.Cuttlebone morphology limits habitat depth in eleven speciesof sepia(Cephalopoda:Sepiidae).BiologicalBulletin,2000,198:404-414

NEW STRUCTURAL DESIGN OF SPHERICAL PRESSURE HULL FOR DEEP-SEA SUBMERSIBLES:A MULTILAYER AND PRESSURE REDISTRIBUTION APPROACH1)

Zhang Yin*,?,2)Liu Xiaoming*,?Lei Xianqi*,?Sun Chengqi*Fang Xin*Wei Yujie*,?
*(State Key Laboratory of Nonlinear Mechanics，Institute of Mechanics，Chinese Academy of Sciences，Beijing100190，China)
?(School of Engineering Science，University of Chinese Academy of Sciences，Beijing100049，China)

The pressure hull is the most important key structure to a deep-sea submersible,whose safety and overall performance depend on the integrity of the pressure hull.The current single shell design of the pressure hull is reviewed and a new design based on the multilayer and pressure redistribution mechanism is also proposed in this study.The design is inspired by two deep-sea animals of sperm whale and nautilus.The two mechanisms,the multilayer structure of sperm whale and the shell substructure divided by the septa of nautilus,are synthesized in the design to improve the pressurebearing capacity of the hull.Compared with the one layer structure of the pressure hull,the new design increases both the structural strength and buckling load.Furthermore,the new design can also signi fi cantly improve the overall performance of the deep-sea submersible by enlarging the hull volume,enhancing its reliability and lessening the technical challenges of fabricatingultra-thick shell structure.Because the truss structure is introduced in the newshell design,which e ff ectively divides a spherical shell into several substructures of cylindrical shell,a formula is strictly derived to evaluate the truss e ff ect on the buckling load bearing capability of a cylindrical shell.The di ff erence between the newly derived formula and the formula of the Taylor basin,which is summarized from the experimental data and widely used as a standard for the current design of submersible shell,is only 6%.Furthermore,this new derivation also lays a solid theoretical foundation for our new shell design.

pressure hull,deep-sea submersible,structure,strength,buckling

U674.941

A doi：10.6052/0459-1879-17-156

2017–05–04 收稿，2017–09–20 錄用,2017–09–28 網(wǎng)絡(luò)版發(fā)表.

1)中國科學(xué)院B類戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDB22020201)和國家自然科學(xué)基金(11372321,11772335)資助項(xiàng)目.

2)張吟，研究員，主要研究方向：振動，固體力學(xué)中的表界面效應(yīng).E-mail:zhangyin@lnm.imech.ac.cn

張吟,劉小明,雷現(xiàn)奇,孫成奇,方新,魏宇杰.基于分層分壓結(jié)構(gòu)的新型潛水器耐壓殼結(jié)構(gòu)設(shè)計(jì).力學(xué)學(xué)報(bào),2017,49(6):1231-1242

Zhang Yin,Liu Xiaoming,Lei Xianqi,Sun Chengqi,Fang Xin,Wei Yujie.New structural design of spherical pressure hull for deep-sea submersibles:A multilayer and pressure redistribution approach.Chinese Journal of Theoretical and Applied Mechanics,2017,49(6):1231-1242