朱永梅，戴永健，朱俊臣，張 建,2

(1. 江蘇科技大學(xué)，江蘇 鎮(zhèn)江 212003；2. 中國(guó)船舶科學(xué)研究中心，江蘇 無(wú)錫 214082)

0 引 言

深海潛水器是進(jìn)行大洋勘查與深?？茖W(xué)研究的重要運(yùn)載工具，受到世界各國(guó)的普遍重視[1]。耐壓殼作為潛水器的關(guān)鍵組成部分，在潛水器的設(shè)計(jì)開(kāi)發(fā)進(jìn)程中舉足輕重，因此在對(duì)海洋的認(rèn)知、開(kāi)發(fā)過(guò)程中耐壓殼的研究占有重要位置[2 - 3]。

球形耐壓殼體具有強(qiáng)度高、穩(wěn)定性好等特點(diǎn)，已成為潛水器耐壓殼體的首選，但是其空間利用率較低，人們?cè)谂搩?nèi)工作時(shí)多數(shù)空間不易被利用，造成了較高的空間資源浪費(fèi)；其次由于其橫截面積較大從而水阻系數(shù)較高，不便于水下流動(dòng)作業(yè)[4]。蛋殼屬于薄殼結(jié)構(gòu)，它的母線是一種正高斯曲線，這種結(jié)構(gòu)不僅能給人一種美學(xué)的藝術(shù)感，還能節(jié)約材料，其次蛋形結(jié)構(gòu)還符合圓頂原理，能夠很好地分散來(lái)自外在的均布載荷[5]。

混凝土具有抗壓、抗腐蝕、抗凍以及防滲透等特性，為其運(yùn)用于海底建筑建設(shè)中提供了保障。通過(guò)改變混凝土的水灰比，可以得到不同性能的混凝土材料，這更多的體現(xiàn)了其運(yùn)用的靈活性[6 - 8]。混凝土結(jié)構(gòu)在海洋建筑中的運(yùn)用已經(jīng)非常普遍，如橋梁、堤壩、海上工作臺(tái)，但將混凝土運(yùn)用于潛水器耐壓殼的建設(shè)中，至今仍未見(jiàn)報(bào)道。因此，綜合考慮蛋殼的耐壓性能和混凝土的耐久性，研究并設(shè)計(jì)能在水下工作的耐壓殼結(jié)構(gòu)，可大大節(jié)約海洋開(kāi)發(fā)、研究的成本，這不僅對(duì)未來(lái)海洋資源的開(kāi)發(fā)，還對(duì)海軍軍事發(fā)展都顯得尤為重要。

本文結(jié)合蛋形結(jié)構(gòu)和混凝土的耐壓特性，完成2 000 m水深耐壓殼的結(jié)構(gòu)設(shè)計(jì)，分析其在承受均布載荷時(shí)，理想狀態(tài)、考慮初始幾何缺陷以及考慮蠕變時(shí)的失穩(wěn)形態(tài)以及極限載荷，得出了2年內(nèi)蛋形混凝土耐壓殼的蠕變曲線。

1 蛋形混凝土耐壓殼的結(jié)構(gòu)設(shè)計(jì)

1.1 混凝土配料的選擇

國(guó)家混凝土結(jié)構(gòu)建設(shè)規(guī)范根據(jù)混凝土的強(qiáng)度劃分若干等級(jí)，常用的有C20-C50，設(shè)計(jì)時(shí)根據(jù)不同的工作環(huán)境，選擇合適的混凝土強(qiáng)度等級(jí)，其中基礎(chǔ)在海水里混凝土標(biāo)號(hào)不能低于C30[9]。本文設(shè)計(jì)2 000 m水深蛋形耐壓殼，工程主體建設(shè)材料選用混凝土，計(jì)算載荷為：

通過(guò)計(jì)算可得需承受靜水壓力32.67 MPa，因此，混凝土強(qiáng)度等級(jí)初選C35，配料比為：水泥209 kg/m3；砂731 kg/m3；石 1 144 kg/m3；水 152 kg/m3；粉煤灰 95 kg/m3；減水劑 4.490 kg/m3。

常規(guī)建筑中混凝土澆筑之后，有必要進(jìn)行養(yǎng)護(hù)，其原因是剛剛澆筑好的混凝土仍處于熔漿質(zhì)，不具備較高的強(qiáng)度等力學(xué)特性，待其硬化后，才能充分發(fā)揮混凝土的力學(xué)特性。C35混凝土不同齡期抗壓強(qiáng)度見(jiàn)表1。C35混凝土養(yǎng)護(hù)28 d后，其抗壓強(qiáng)度可達(dá)到44.62 MPa，用于 2 000 m 水深的耐壓殼符合要求，其材料基本參數(shù)為：彈性模量，泊松比，材料密度kg/m3。

1.2 蛋形耐壓殼的結(jié)構(gòu)設(shè)計(jì)

國(guó)內(nèi)外球形耐壓殼直徑一般為2 m，本文采用等體積的設(shè)計(jì)方法，設(shè)計(jì)蛋形耐壓殼。通過(guò)前期的試驗(yàn)研究分析[10]，采用N-R函數(shù)[11]作為蛋形耐壓殼輪廓設(shè)計(jì)的數(shù)學(xué)模型，蛋形耐壓殼的體積公式[12]以及蛋形系數(shù)公式[12]為：

表1 C35 混凝土不同齡期抗壓強(qiáng)度Tab.1 Compressive strength of concrete at different age

式中：V為體積；L為蛋形耐壓殼長(zhǎng)軸尺寸；B為短軸尺寸；SI為蛋形系數(shù)。

根據(jù)對(duì)蛋形系數(shù)的統(tǒng)計(jì)分析，蛋形系數(shù)取0.69為最佳。通過(guò)計(jì)算可得到蛋形耐壓殼的長(zhǎng)軸尺寸為2.56 m，短軸尺寸為1.77 m。

根據(jù)薄殼理論[13-14]，可確定蛋形耐壓殼的經(jīng)向應(yīng)力和緯向應(yīng)力，分別為：

式中：Ps為計(jì)算載荷；為第一曲率半徑；為第二曲率半徑；t為蛋形耐壓殼的厚度。

蛋形耐壓殼的屈曲臨界載荷，采用Mushtari方程[15]表述：

式中：PQ為屈曲臨界載荷；μ為泊松比；E為彈性模量；為第一曲率半徑均值；為第二曲率半徑均值。

綜合考慮極限載荷和屈曲載荷，蛋形耐壓殼的厚度可由式（8）確定。通過(guò)等厚設(shè)計(jì)方法確定厚度為t=55.29 mm。

2 蛋形混凝土耐壓殼力學(xué)特性研究

對(duì)于耐壓殼的力學(xué)特性研究，主要集中在強(qiáng)度和屈曲兩方面的研究，通過(guò)研究2 000 m水深下的蛋形耐壓殼最終失穩(wěn)形式為屈曲失穩(wěn)[4]。可以認(rèn)為當(dāng)耐壓殼發(fā)生屈曲破壞時(shí)，其強(qiáng)度仍未達(dá)到極限值，因此本文只對(duì)混凝土耐壓殼的屈曲展開(kāi)研究。

2.1 理想狀態(tài)下蛋形混凝土耐壓殼屈曲分析

在混凝土薄殼結(jié)構(gòu)力學(xué)穩(wěn)定性驗(yàn)算方法上，國(guó)際殼體和空間結(jié)構(gòu)協(xié)會(huì)（IASS）推薦了半經(jīng)驗(yàn)的驗(yàn)算方法，其公式如下：

使用三維建模軟件根據(jù)N-R方程建立蛋形耐壓殼輪廓曲線，抽取模型內(nèi)部保留規(guī)定的厚度，得到蛋形殼體結(jié)構(gòu)，使用Ansa軟件進(jìn)行有限元計(jì)算前處理，采用網(wǎng)球劃分法進(jìn)行網(wǎng)格劃分。由于耐壓殼在2 000 m水深下工作，因此其承受均布載荷，理論上模型是不受任何約束力的，但是為了消除模型的剛性位移，在模型中任意選取3點(diǎn)限制模型6個(gè)自由度的剛性位移。通過(guò)有限元計(jì)算可得屈曲臨界載荷為51.91 MPa。

圖1為蛋形耐壓殼一階線彈性屈曲失穩(wěn)模式圖。

圖1 屈曲失穩(wěn)模式圖Fig.1 Buckling model

由圖1可知，蛋形耐壓殼在外載荷情況下，失穩(wěn)部位主要發(fā)生在中部區(qū)域，且呈現(xiàn)波浪狀失穩(wěn)模式，具有5個(gè)波峰值，是一種具有正高斯曲線殼的典型失穩(wěn)模式，沿緯向呈現(xiàn)有規(guī)律的凹進(jìn)與凸出，波峰與波谷數(shù)均為5，最大位移出現(xiàn)在如圖所示波谷處，值為1.006 mm。

2.2 缺陷狀態(tài)下蛋形混凝土耐壓殼屈曲分析

殼類結(jié)構(gòu)的穩(wěn)定性對(duì)初始缺陷十分敏感，在求解殼類結(jié)構(gòu)穩(wěn)定性時(shí)，初始缺陷必須予以考慮，它能引起殼體臨界載荷的下降，無(wú)論是在求解線性還是非線性時(shí)，初始幾何缺陷的特征將十分明顯。對(duì)于有限元模型，結(jié)構(gòu)初始缺陷的分布趨勢(shì)，一般可以假設(shè)為結(jié)構(gòu)的某一模態(tài)或者某幾個(gè)模態(tài)的組成。特征值屈曲模態(tài)，是目前非常普遍的計(jì)算殼體初始缺陷的一種方法。即假設(shè)初始缺陷按照最低階特征值屈曲模態(tài)形式分布。

從開(kāi)始的理想狀態(tài)逐漸增大模型的初始缺陷展開(kāi)研究，分別設(shè)置缺陷值為厚度的0.8，0.6，0.4以及0.2倍，研究隨著缺陷的增加，模型的屈曲臨界載荷的變化情況。

圖2 幾何缺陷對(duì)屈曲的影響Fig.2 The effect of geometric imperfections on buckling

表2 缺陷狀態(tài)下的臨界屈曲載荷Tab.2 Critical buckling load under the condition of imperfections

通過(guò)圖2可得隨著缺陷的增加，耐壓殼結(jié)構(gòu)承載能力逐漸降低，屈曲臨界載荷不斷下降，只是下降的速率不斷減小，說(shuō)明初始幾何缺陷對(duì)混凝土耐壓殼屈曲有一定的影響，相比而言缺陷發(fā)展初期對(duì)屈曲的影響較大，后期影響較小，這也就是混凝土在裂縫初期對(duì)構(gòu)件的穩(wěn)定性影響較為顯著的原因。

圖3 平衡路徑圖Fig.3 Equilibrium path

圖4 為圖3中平衡路徑末端點(diǎn)的后屈曲失穩(wěn)模式圖，由圖可見(jiàn)，混凝土蛋形耐壓殼后屈曲失穩(wěn)模式不同于其線彈性失穩(wěn)模式，其主要呈現(xiàn)局部凹坑失穩(wěn)，且失穩(wěn)部位位于混領(lǐng)土蛋形耐壓殼的中部區(qū)域，表明蛋形耐壓殼中部為危險(xiǎn)區(qū)域，極易發(fā)生屈曲失穩(wěn)。

圖4 屈曲失穩(wěn)模式圖Fig.4 Buckling model

2.3 混凝土蠕變對(duì)耐壓殼屈曲的影響

本文使用Abaqus模擬混凝土耐壓殼在2年時(shí)間內(nèi)恒定溫度下，受到恒定載荷的情況下的蠕變?nèi)^(guò)程，觀察其應(yīng)力分布、應(yīng)力應(yīng)變過(guò)程以及屈曲臨界載荷。在Abaqus中模型的前處理，如單元類型的選擇、網(wǎng)格的劃分等都與上文相同，但是材料設(shè)置添加了混凝土塑性參數(shù)的定義。在Abaqus中，混凝土本構(gòu)關(guān)系由三維模型簡(jiǎn)化為平面模型，因此在計(jì)算時(shí)只需要定義混凝土單軸受壓的應(yīng)力應(yīng)變曲線即可，混凝土的應(yīng)力應(yīng)變曲線如圖5所示。

在計(jì)算蠕變應(yīng)變時(shí)，容差的設(shè)置十分關(guān)鍵，它與增量步的大小有著密切的關(guān)系，容差越小，增量步也跟著降低，因此在計(jì)算時(shí)容差值太大會(huì)影響精度，太小則會(huì)使計(jì)算時(shí)間變大導(dǎo)致系統(tǒng)崩潰。一般取1E-4～1E-6左右，容差可根據(jù)具體計(jì)算作調(diào)整。本文沿著經(jīng)線，取蛋形耐壓殼結(jié)構(gòu)赤道與鈍端之間中點(diǎn)處取一點(diǎn)A，中間赤道部位取點(diǎn)B，赤道與尖端之間中點(diǎn)處取點(diǎn)C。研究隨時(shí)間的發(fā)展各點(diǎn)處的變形情況，結(jié)果如圖6所示。

圖5 混凝土本構(gòu)關(guān)系曲線Fig.5 The stress-strain curve of concrete

圖6 蠕變曲線Fig.6 Creep curve

從圖6可看出，隨著時(shí)間的發(fā)展，耐壓殼各點(diǎn)處均發(fā)生蠕變，但是在相同時(shí)間，同等載荷的情況下耐壓殼B點(diǎn)處的蠕變發(fā)展比A點(diǎn)的大，A點(diǎn)處蠕變比C點(diǎn)大，這也證明了耐壓殼內(nèi)部應(yīng)力分布規(guī)律為赤道處應(yīng)力比兩端大，鈍端比尖端應(yīng)力大；混凝土破壞分為3個(gè)階段，第1階段為蠕變發(fā)展階段，第2階段為加速階段，第3階段為蠕變破壞階段，本例在加載2年時(shí)間內(nèi)，混凝土蠕變未發(fā)展為第3階段，因此可以認(rèn)為未達(dá)到破壞極限。最終計(jì)算蠕變后屈曲極限載荷為 48.34 MPa。

3 結(jié) 語(yǔ)

結(jié)合蛋形結(jié)構(gòu)和混凝土的耐壓特性，完成了2 000 m水深耐壓殼的結(jié)構(gòu)設(shè)計(jì)，分析了其在承受均布載荷時(shí)，理想狀態(tài)、考慮初始幾何缺陷以及考慮蠕變時(shí)的失穩(wěn)形態(tài)以及極限載荷，可得到如下結(jié)論：

1）耐壓殼結(jié)構(gòu)在承受均布載荷時(shí)，緯向應(yīng)力變化幅度大于經(jīng)向應(yīng)力，因此設(shè)計(jì)耐壓殼時(shí)應(yīng)重點(diǎn)考慮緯向應(yīng)力；蛋形結(jié)構(gòu)兩端承受的應(yīng)力較小，中間段承受的應(yīng)力較大，相比而言尖端承受應(yīng)力比鈍端大，此結(jié)果說(shuō)明失穩(wěn)最先發(fā)生于耐壓殼中部，和實(shí)驗(yàn)結(jié)果相符合。

2）隨著缺陷的增加，耐壓殼結(jié)構(gòu)承載能力逐漸降低，屈曲臨界載荷下降，但是下降的速率逐漸減小，說(shuō)明在缺陷發(fā)展初期對(duì)屈曲的影響較大，后期影響較小，此結(jié)果說(shuō)明混凝土在裂縫初期對(duì)構(gòu)件穩(wěn)定性影響更為顯著。

3）隨著時(shí)間的增大，耐壓殼各點(diǎn)處均發(fā)生蠕變，但是在相同時(shí)間，同等載荷的情況下耐壓殼蠕變進(jìn)展依次為：赤道處，鈍端處，尖端處。發(fā)生蠕變后耐壓殼屈曲極限載荷與理想狀態(tài)相比下降了6.88%。