朱永梅,戴永健,朱俊臣,張 建,2
(1. 江蘇科技大學(xué),江蘇 鎮(zhèn)江 212003;2. 中國(guó)船舶科學(xué)研究中心,江蘇 無(wú)錫 214082)
深海潛水器是進(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)蛋形混凝土耐壓殼的蠕變曲線。
國(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。
國(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。
對(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)研究。
在混凝土薄殼結(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。
殼類結(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
本文使用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。
結(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%。