王 丹，萬(wàn)正權(quán)

（中國(guó)船舶科學(xué)研究中心，江蘇 無(wú)錫 214082）

考慮初始形狀影響的耐壓球殼臨界載荷簡(jiǎn)化計(jì)算公式

王 丹，萬(wàn)正權(quán)

（中國(guó)船舶科學(xué)研究中心，江蘇 無(wú)錫 214082）

文章研究了在載人深潛器結(jié)構(gòu)中承受高壓、保障潛員安全的耐壓球殼。文中通過(guò)對(duì)經(jīng)典理論公式的分析和合理假設(shè)，引入幾何修正系數(shù)k細(xì)化了球殼臨界載荷計(jì)算表達(dá)式，該幾何修正系數(shù)體現(xiàn)了球殼本身的初始形狀，如初撓度、局部不圓度及板厚減薄等因素對(duì)其的影響。并通過(guò)回歸分析得到計(jì)算數(shù)據(jù)與試驗(yàn)值吻合較好。同時(shí)，運(yùn)用該公式對(duì)鈦合金球殼臨界載荷進(jìn)行預(yù)報(bào)，同樣有較高精度。

耐壓球殼；臨界載荷；載人潛水器；初撓度

1 引 言

球殼是深潛裝備、航天航空、壓力容器及核工業(yè)等領(lǐng)域常用的一種結(jié)構(gòu)形式。這與球殼自身具備的物理特性有著密切的關(guān)系。與相同質(zhì)量的其他復(fù)雜結(jié)構(gòu)相比較，球殼具有容重比高的特點(diǎn)，只要布置合理就能夠最大可能利用空間來(lái)承載；理想球殼的外表面受力均勻，在周向和徑向的應(yīng)力值相同，計(jì)算時(shí)只需關(guān)注膜應(yīng)力數(shù)值即可控制設(shè)計(jì)參數(shù)；因其結(jié)構(gòu)簡(jiǎn)單，在理論計(jì)算或數(shù)值模擬過(guò)程中選擇的參數(shù)較少，變量?jī)H剩厚度與半徑。無(wú)論是深潛還是航天航空領(lǐng)域，選擇球形結(jié)構(gòu)作為載人艙都是比較適當(dāng)。

球形載人艙肩負(fù)保障水下潛員安全、提供生命保障及作業(yè)條件方面的重任，球殼一旦發(fā)生失穩(wěn)、破壞，后果將不堪設(shè)想。一直以來(lái)，球殼的安全性受到科研人員的廣泛關(guān)注，各國(guó)科研人員對(duì)各種材料制成的球殼結(jié)構(gòu)進(jìn)行了大量的理論與試驗(yàn)研究，取得了豐碩的科研成果。而大多數(shù)科研是針對(duì)球殼的安全性進(jìn)行的，如將球殼的安全性參數(shù)化，可用極限承載能力來(lái)衡量，而承載能力一般可以用強(qiáng)度和穩(wěn)定性這兩個(gè)特征來(lái)表示。隨著科學(xué)技術(shù)的發(fā)展，材料學(xué)、晶像學(xué)在微觀和宏觀領(lǐng)域的研究，用于加工球殼的材料其本身強(qiáng)度均已達(dá)到深潛要求。然而，球殼成形的過(guò)程與許多不利因素有著緊密的關(guān)聯(lián)，如材料的非線性變形、殼體與理想形狀之間的初始偏差（初撓度）、板材滾軋減薄、焊接殘余應(yīng)力及其他工藝內(nèi)應(yīng)力等。這些因素在一定程度上都會(huì)影響到球殼的穩(wěn)定性，它們的存在會(huì)使球殼損失一定的承載能力。

到目前為止，我國(guó)科研人員對(duì)球殼問(wèn)題進(jìn)行了長(zhǎng)期不懈的研究，如羅培林[1]從實(shí)際工程使用中球殼的歐拉應(yīng)力一般高于材料的屈服極限出發(fā)，應(yīng)用兩個(gè)經(jīng)驗(yàn)常數(shù)來(lái)反映生產(chǎn)工藝過(guò)程對(duì)球殼承載能力的影響；在文獻(xiàn)[2]中，作者根據(jù)梁柱失效原理，用極限應(yīng)力函數(shù)來(lái)分析球殼的穩(wěn)定性，獲得較好的計(jì)算結(jié)果；劉濤[3]通過(guò)對(duì)材料非彈性模量的數(shù)據(jù)擬合分析，用直線擬合切線彈性模量，獲得高精度耐壓球殼失穩(wěn)臨界壓力計(jì)算公式；陸蓓、劉濤、崔維成[4]通過(guò)有限元分析方法，提出合理的有限元分析模型并研究了初撓度的影響，給出了完善的耐壓球殼極限強(qiáng)度與經(jīng)典值的比值曲線以及初撓度效應(yīng)曲線；俞銘華，王自力等[5]對(duì)潛水器耐壓球殼結(jié)構(gòu)研究進(jìn)展進(jìn)行了綜述，并利用有限元對(duì)不同深度潛水器進(jìn)行分析[6]，獲得初始缺陷對(duì)耐壓球殼承載能力的影響曲線；潘彬彬，崔維成[7]對(duì)深潛器載人球殼各國(guó)船級(jí)社規(guī)范設(shè)計(jì)公式進(jìn)行了比較，指出公式計(jì)算的球殼厚度大多偏厚，需要制定統(tǒng)一的球殼設(shè)計(jì)公式；王丹，邱昌賢等[8]通過(guò)球殼模型試驗(yàn)及有限元數(shù)值模擬與泰勒水池公式、俄羅斯公式計(jì)算結(jié)果進(jìn)行比較，分析初撓度對(duì)球殼的影響，提出這兩種公式計(jì)算結(jié)果較為保守，有必要對(duì)其修正系數(shù)進(jìn)行調(diào)整。

由前述研究可以看出穩(wěn)定性問(wèn)題求解大多通過(guò)理論計(jì)算和有限元分析來(lái)研究，而最有效的檢驗(yàn)方法是進(jìn)行模型試驗(yàn)，通過(guò)模型試驗(yàn)可以獲得球殼實(shí)際失穩(wěn)壓力、失穩(wěn)形態(tài)以及臨界應(yīng)力，有助于對(duì)修正系數(shù)進(jìn)行調(diào)整，得到更精準(zhǔn)的計(jì)算公式。本文從球殼的經(jīng)典理論公式與經(jīng)驗(yàn)公式著手分析，簡(jiǎn)化了球殼臨界載荷計(jì)算公式，并利用模型試驗(yàn)結(jié)果進(jìn)行回歸分析，驗(yàn)證幾何修正系數(shù)的臨界載荷計(jì)算公式。

2 球殼的臨界載荷計(jì)算

2.1 理想球殼理論臨界載荷

在進(jìn)行球殼極限承載能力的理論分析時(shí)一般會(huì)做基本假定，假設(shè)球殼材料無(wú)強(qiáng)化，具有基本的拉伸曲線圖；初撓度形狀為圓形，失穩(wěn)形狀與初撓度形狀相同，可以通過(guò)求解軸對(duì)稱問(wèn)題來(lái)解決初撓度對(duì)結(jié)構(gòu)承載能力的影響問(wèn)題。

在實(shí)際分析方法中，由理想材料制成的具有正確形狀的球殼理論臨界載荷一般采用著名的崔林（Zoelly）關(guān)系式表示，即

其中：v為泊松比；E為彈性模量，MPa；R為球殼中面半徑，mm；t為殼板平均板厚，mm；對(duì)于鋼制球殼，其泊松比v=0.3，公式（1）可以表示為：

2.2 非理想球殼理論臨界載荷

對(duì)非理想球殼進(jìn)行分析時(shí)，需要考慮球殼的加工工藝、殼板真實(shí)厚度與額定厚度之間的偏差、球殼局部半徑變化、真實(shí)球殼與理想球殼之間的偏差、材料性能、焊接殘余應(yīng)力等因素。厚度、半徑等參數(shù)均可通過(guò)測(cè)量工具獲得較為精確的數(shù)值，一般選擇初撓度進(jìn)行理論分析。球殼表面產(chǎn)生的初撓度對(duì)球殼本身產(chǎn)生附加彎曲作用，隨著壓力的不斷增加，附加彎曲的作用也逐漸增強(qiáng)，引起塑性變形的力與彎矩也會(huì)隨之增大。當(dāng)球殼殼板不能抵抗初撓度引起的附加力和力矩時(shí)，就會(huì)發(fā)生失穩(wěn)。因此，要獲得真實(shí)的非理想球殼理論臨界載荷就需要對(duì)理想球殼理論臨界載荷進(jìn)行修正。

Paliy教授在1991年就發(fā)表了關(guān)于殼體承載能力的文章[9]，其中帶有修正系數(shù)的公式如下：

將（2）式代入（3）式中得到相應(yīng)的理論臨界應(yīng)力關(guān)系式為

該公式已經(jīng)應(yīng)用在俄羅斯載人潛水器及深潛裝備的分類和建造規(guī)范之中[10]。其中，α取值為：柱殼取0.95，球殼取1.2。β系數(shù)與很多不利影響有關(guān)，主要考慮了材料的機(jī)械性能、初撓度、焊接內(nèi)應(yīng)力及其他工藝內(nèi)應(yīng)力等。

這里，β1=1/（1+（2.8+f ′）f2/3），f′為相對(duì)初撓度，f′=f/t，f為球殼偏離理想形狀的最大初撓度；σ′為相對(duì)應(yīng)力，σ′=σ/σs，σ 為臨界載荷情況下對(duì)應(yīng)的臨界應(yīng)力。

將α取值代入（5）式中并結(jié)合（2）式得到

2.3 帶幾何修正系數(shù)k的臨界載荷計(jì)算公式

由（4）式和（7）式可以看出，兩種情況下的理論臨界應(yīng)力比值中隱藏了修正系數(shù)β，則實(shí)際情況下的非理想球殼的承載能力關(guān)系式（6）可以變換為：

此時(shí)，利用（6）式，實(shí)際臨界載荷情況下的臨界應(yīng)力可以寫作：

如何能將球殼承載能力計(jì)算公式繼續(xù)細(xì)化，不僅能直觀地看出球殼半徑與殼板厚度之間的關(guān)系，還能更加準(zhǔn)確地計(jì)算出球殼的承載能力，這需要進(jìn)行合理假設(shè)。引入材料屈服極限σs，假設(shè)其與σ之間的關(guān)系為

將（3）式和（9）式代入（8）式可以得到：

其中：k為幾何修正系數(shù)，與球殼的厚度和半徑變化量以及初撓度值有著密切的關(guān)系。將厚度和半徑變化量的影響記作kt、kr，其值可以通過(guò)測(cè)量得到真實(shí)球殼的半徑和板厚，并用（11）式和（12）式；將初撓度在臨界弧長(zhǎng)內(nèi)的變化量的影響記作kf，kf=（1- fmax′ ），fmax′是最大相對(duì)初撓度，fmax′=fmax/t ，fmax為球殼實(shí)測(cè)最大初撓度值，有正負(fù)之分，負(fù)值表示凹陷，正值表示凸起，一般情況下認(rèn)為最危險(xiǎn)的狀態(tài)是凹陷，在（12）式中fmax只取最大凹陷值。

kt、kr按照下式進(jìn)行計(jì)算：

其中：tc為殼板平均厚度，通過(guò)初始形狀測(cè)量獲得，這樣可以獲得帶有幾何修正系數(shù)的球殼臨界載荷計(jì)算公式：

3 球殼模型臨界載荷試驗(yàn)值與計(jì)算值的比較

3.1 球殼的靜水外壓試驗(yàn)

中國(guó)船舶科學(xué)研究中心水下工程結(jié)構(gòu)試驗(yàn)室針對(duì)2個(gè)（Q1-Q2）參數(shù)（直徑、厚度、材料）相同，3個(gè)參數(shù)不同的球殼（Q3-Q5）進(jìn)行初始形狀測(cè)量，得到其初撓度值、局部凹凸度和殼板平均厚度等，之后對(duì)這些球殼分別進(jìn)行了靜水外壓試驗(yàn)。

Q1-Q2模型直徑較小，為半球焊接，最大相對(duì)初撓度值小于0.15，失穩(wěn)壓力較高；Q3-Q4模型直徑較大，為瓣片式焊接（赤道面無(wú)對(duì)接焊縫），焊縫較多，最大相對(duì)初撓度達(dá)到0.43，失穩(wěn)壓力偏低；Q5模型直徑較大，亦瓣片式焊接（赤道面有對(duì)接焊縫），焊縫與焊縫錯(cuò)開(kāi)對(duì)接，最大相對(duì)初撓度僅0.15，在20 MPa靜水外壓情況下未發(fā)生失穩(wěn)。球殼模型相關(guān)參數(shù)及失穩(wěn)壓力[11-13]列入表1，模型失穩(wěn)后形狀見(jiàn)圖1和圖2。

表1 球殼模型主要參數(shù)及失穩(wěn)壓力Tab.1 Main parameters and buckling loads of the spherical hull models

圖1 Q1、Q2失穩(wěn)形態(tài)Fig.1 The buckling shapes of Q1 and Q2

圖2 Q3、Q4失穩(wěn)形態(tài)Fig.2 The buckling shapes of Q3 and Q4

3.2 計(jì)算值與試驗(yàn)值的比較

應(yīng)用本文給出的計(jì)算公式，根據(jù)5個(gè)球殼模型的實(shí)際測(cè)量數(shù)據(jù)計(jì)算臨界載荷。將計(jì)算值與試驗(yàn)值進(jìn)行比較，表2中列出了臨界載荷計(jì)算值與試驗(yàn)值之間的偏差。

表2 利用本文給出方法對(duì)球殼模型臨界載荷的計(jì)算Tab.2 Calculation of the critical loads for spherical hull models by the present method

圖3是計(jì)算數(shù)據(jù)與試驗(yàn)值的趨勢(shì)線。分析了在k值一定的情況下，最大相對(duì)初撓度與臨界載荷之間的關(guān)系。利用回歸分析對(duì)試驗(yàn)結(jié)果與計(jì)算結(jié)果進(jìn)行比對(duì)，其相關(guān)系數(shù)達(dá)到99%以上。

表3 三種方法計(jì)算結(jié)果與試驗(yàn)結(jié)果的比較Tab.3 Comparison between calculation results of the three methods and test results

圖3 本文方法計(jì)算值與試驗(yàn)值趨勢(shì)線Fig.3 Calculation and test values trend lines

圖4 當(dāng)kt/kr=0.95時(shí)的臨界載荷曲線Fig.4 Critical load curves for kt/kr=0.95

本文提供的計(jì)算方法與試驗(yàn)值之間的偏差較小，最大偏差為6.7%；用俄羅斯方法的最大偏差為15.4%，而用泰勒水池方法的偏差為30%，通過(guò)算例比較，本文方法更接近實(shí)際情況。

圖4中以5個(gè)球殼模型為例，給出利用本文所給計(jì)算公式計(jì)算得到當(dāng)kt/kr=0.95時(shí)的臨界載荷曲線。可以看出在薄殼范圍內(nèi)，球殼的臨界載荷隨著半徑厚度比的增加而減小；在相同R/t的情況下，臨界載荷隨著初撓度的增加而降低。

4 應(yīng)用本方法對(duì)鈦合金耐壓球殼的臨界壓力進(jìn)行計(jì)算

根據(jù)鈦合金球耐壓球殼模型的主要參數(shù)，運(yùn)用本文方法進(jìn)行臨界載荷預(yù)報(bào)，與試驗(yàn)值對(duì)比結(jié)果見(jiàn)表4。

表4 鈦合金球殼模型主要參數(shù)及臨界載荷計(jì)算結(jié)果Tab.4 Main parameters and calculation results of critical loads for titanium alloy spherical hull models

可以看出運(yùn)用本文方法預(yù)報(bào)的鈦合金球殼模型臨界載荷與其靜水外壓試驗(yàn)值比較吻合，最大偏差不超過(guò)7.3%。該方法可以運(yùn)用到鈦合金球殼臨界載荷的預(yù)報(bào)中。

5 結(jié) 論

本文未考慮焊接殘余應(yīng)力及其他裝配內(nèi)應(yīng)力對(duì)球殼臨界載荷的影響，下面給出根據(jù)本文簡(jiǎn)化公式計(jì)算及分析得到的結(jié)果：

（1）簡(jiǎn)化公式需要得知精確的初始形狀測(cè)量數(shù)據(jù)，如板厚、半徑、初撓度或不圓度等才能得到更接近實(shí)際情況的計(jì)算結(jié)果；當(dāng)球殼為精加工模型時(shí)（球殼的半徑變化不大，且初始缺陷較小），影響承載能力的因素主要是殼板厚度；

（2）當(dāng)球殼的半徑、厚度變化量的關(guān)系kr/kt變化不大時(shí)，初撓度對(duì)球殼的影響更加顯著；當(dāng)初撓度值增加5%左右，臨界載荷會(huì)隨之降低，泰勒水池方法降低約11%[8]，俄羅斯方法降低約5%，本文方法隨初撓度的增加表現(xiàn)出非線性變化，初撓度每增加5%，臨界載荷降低量從5%～7%逐漸增加；

（3）在薄殼范圍內(nèi)（R/t＞20），臨界載荷隨著R/t的增大而減小，本文方法適用于鈦合金制造的球殼結(jié)構(gòu)臨界載荷計(jì)算；

（4）本文方法預(yù)報(bào)的臨界載荷與試驗(yàn)值的某些偏差超過(guò)5%，需要對(duì)簡(jiǎn)化計(jì)算公式進(jìn)行進(jìn)一步的修正，使得計(jì)算結(jié)果與試驗(yàn)值更加接近；

（5） 根據(jù)強(qiáng)度標(biāo)準(zhǔn)（σ≥0.85σs）的要求反推得出：k=kf·kt/kr≥0.85；此時(shí)，即從強(qiáng)度方面限制了球殼厚度、半徑及初撓度的變化量參與求解球殼穩(wěn)定性問(wèn)題。

[1]羅培林.厚球殼的極限載荷-均勻壓力-計(jì)算與模擬試驗(yàn)方法的改進(jìn),1976,手寫稿.1994-2012 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net.[J/OL].

[2]羅培林.用極限分析方法分析梁柱和球殼的穩(wěn)定性[J].哈爾濱工程大學(xué)學(xué)報(bào),1981,2:21-50.

[3]劉 濤.大深度潛水器耐壓殼體彈塑性穩(wěn)定性簡(jiǎn)易計(jì)算方法[J].中國(guó)造船,2001,42(3):8-14.

[4]陸 蓓,劉 濤,崔維成.深海載人潛水器耐壓球殼極限強(qiáng)度研究[J].船舶力學(xué),2004,8(1):51-58.Lu Bei,Liu Tao,Cui Weicheng.Ultimate strength of pressure spherical hull in deep-sea manned submersibles[J],Journal of Ship Mechanics,2004,8(1):51-58.

[5]俞銘華,王自力,李良碧,王仁華.大深度載人潛水器耐壓殼結(jié)構(gòu)研究進(jìn)展[J].華東船舶工業(yè)學(xué)院學(xué)報(bào)(自然科學(xué)版),2004,18(4):1-6.

[6]王自力,王仁華,俞銘華,李良碧.初始缺陷對(duì)不同深度載人潛水器耐壓球殼極限承載能力的影響[J].中國(guó)造船.2007,48(2):45-50.

[7]潘彬彬,崔維成.深潛器載人球殼規(guī)范設(shè)計(jì)公式的比較[J].船舶力學(xué),2011,15(3):276-285.Pan Binbin,Cui Weicheng.A comparison of different rules for the spherical pressure hull of deep manned submersibles[J].Journal of Ship Mechanics,2011,15(3):276-285.

[8]Wang Dan,Qiu Changxian,Wang Yongjun,Wan Zhengquan.Theoretical and experimental research of spherical pressure hull[C].船舶與海洋結(jié)構(gòu)力學(xué)學(xué)術(shù)會(huì)議,2011:68-76.

[9]Paliy O M.Weight characteristics,reliability and operational safety of deep-sea submersible hulls[C]//International Symposium on Marine Structures(ISM’91),September 13-14,1991.Shanghai,China,1991:197-199.

[10]Правила классификации и постройки обитаемых[Z].подводных аппаратов и глубоководных водолазных комплексов[Z].1993.

[11]王 丹,楊華偉,王永軍.超大潛深耐壓球殼開(kāi)孔模型試驗(yàn)[R].無(wú)錫:中國(guó)船舶科學(xué)研究中心科技報(bào)告10696,2010.

[12]王 丹,王永軍.超大潛深耐壓球殼無(wú)開(kāi)孔模型試驗(yàn)[R].無(wú)錫:中國(guó)船舶科學(xué)研究中心科技報(bào)告10171,2010.

[13]王 丹,楊華偉,王永軍.超大潛深耐壓球殼模型理論與試驗(yàn)分析報(bào)告[R].無(wú)錫:中國(guó)船舶科學(xué)研究中心科技報(bào)告10351,2010.

Simplified calculation formula for critical load of spherical pressure shell considering the influence of initial shape

WANG Dan,WAN Zheng-quan
(China Ship Scientific Research Center,Wuxi,214082,China)

The spherical pressure shell main structure of manned submersible which is under great pressure and provides the security of human is researched.A new formula for calculating critical load of the spherical shell which bases on classic theory and considers the correction coefficient k is introduced.The influences of initial shape are shown,for example,initial deflections and out-of-roundness and thickness thinning in the manufactory technical,etc.The calculated data agree well with the experimental data by the regression analysis.The critical load of titanium spherical shell predicted by this formula is very precisely.

spherical pressure shell;critical load;manned submersible;initial deflections

U661.4

A

10.3969/j.issn.1007-7294.2014.05.010

1007-7294（2014）05-0557-08

2013-10-28

王 丹（1981-），女，中國(guó)船舶科學(xué)研究中心工程師，E-mail：wangdan21@hotmail.com；

萬(wàn)正權(quán)（1962-），男，研究員，博士生導(dǎo)師。

book=564,ebook=412