潘彬彬，崔維成，葉 聰，劉正元

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

1 引 言

基于氣球原理設(shè)計(jì)的第一代和第二代載人潛水器，均是采用無動(dòng)力下潛和上浮的方式[1-3]。這樣做獲得的最大好處就是可以大大地節(jié)省潛水器上的能源，但需要解決無動(dòng)力潛浮過程中的水動(dòng)力學(xué)問

題[4]。

在“蛟龍”號(hào)設(shè)計(jì)階段，劉正元等就開始研究載人潛水器的無動(dòng)力下潛上浮運(yùn)動(dòng)[5-6]，給出了一套簡化分析方法并應(yīng)用于“蛟龍”號(hào)載人潛水器的無動(dòng)力潛浮設(shè)計(jì)。為了更精確和方便地進(jìn)行海試和今后應(yīng)用下潛的配載，我們?cè)凇膀札垺碧?hào)載人潛水器3 000米級(jí)海試數(shù)據(jù)的基礎(chǔ)上，曾經(jīng)建立了一個(gè)理論上更為精細(xì)的載人潛器無動(dòng)力上浮下潛運(yùn)動(dòng)模型，并根據(jù)該模型開發(fā)了一個(gè)運(yùn)動(dòng)分析系統(tǒng)[7]，用于預(yù)報(bào)“蛟龍”號(hào)5 000米級(jí)海試的配載和速度。但是由于3 000米海試所有潛次的海試數(shù)據(jù)中都存在開動(dòng)推進(jìn)器和調(diào)節(jié)壓載水艙的操作，沒有真正意義上的無動(dòng)力上浮下潛數(shù)據(jù)用于標(biāo)定模型參數(shù)，因此，所建立的無動(dòng)力潛浮運(yùn)動(dòng)分析系統(tǒng)并不十分可靠。在2011年7-8月間進(jìn)行的5 000米級(jí)海試過程中發(fā)現(xiàn)該模型預(yù)報(bào)值比實(shí)際值偏大了約80kg。在這次的5 000米級(jí)海試中，D44潛次首次實(shí)現(xiàn)了全程無動(dòng)力上浮下潛運(yùn)動(dòng)，為改進(jìn)載人潛器無動(dòng)力上浮下潛運(yùn)動(dòng)模型提供了寶貴的數(shù)據(jù)。本文將對(duì)D44潛次的數(shù)據(jù)進(jìn)行分析，在此基礎(chǔ)上改進(jìn)潛器無動(dòng)力上浮下潛運(yùn)動(dòng)模型，期望開發(fā)出真正可靠的“蛟龍”號(hào)載人潛水器的無動(dòng)力潛浮運(yùn)動(dòng)分析系統(tǒng)，為2012年的7 000米級(jí)海試以及今后的應(yīng)用提供可靠的手段，同時(shí)，本文所建立的精細(xì)分析方法也可以應(yīng)用于我國正在進(jìn)行的4 500米載人潛器無動(dòng)力潛浮運(yùn)動(dòng)的設(shè)計(jì)。

2 載人潛水器無動(dòng)力潛浮運(yùn)動(dòng)模型

在文獻(xiàn)[7]中詳細(xì)給出了載人潛水器的無動(dòng)力潛浮運(yùn)動(dòng)分析精細(xì)模型，在該模型中，認(rèn)為無動(dòng)力潛浮運(yùn)動(dòng)過程中潛器所受的外力只有重力W、浮力B及流體阻力Fd。此外考慮到穩(wěn)態(tài)運(yùn)動(dòng)時(shí)其角速度為零等因素，模型建立在固定坐標(biāo)系內(nèi)，如下圖1所示。

潛器在固定坐標(biāo)系中oz0方向的力平衡方程為：

其中M=潛器的質(zhì)量與z0方向附加質(zhì)量之和；0=潛器垂向（z0方向）加速度；W=重力；B=浮力；Fd=潛器在鉛垂方向運(yùn)動(dòng)時(shí)受到的 （z0方向）的阻力。潛器無動(dòng)力下潛上浮的大部分時(shí)間中潛器的運(yùn)動(dòng)可以近似為平穩(wěn)運(yùn)動(dòng)，即近似認(rèn)為平穩(wěn)階段潛器的加速度為零：

圖1 坐標(biāo)系Fig.1 Coordinate system

但（2）式中各部分力與文獻(xiàn)[7]中的模型有所變化，所以本文將對(duì)每個(gè)部分力重新介紹。

2.1 阻力項(xiàng)Fd

如圖1所示，潛器運(yùn)動(dòng)時(shí)存在一定的縱傾角θ，文獻(xiàn)[6]推導(dǎo)了潛器速度方向（攻角）α和縱傾角θ的關(guān)系為：

其中：Cx（α）為潛器縱向（x軸方向）運(yùn)動(dòng)阻力系數(shù)，Cz（α）為潛器垂向（z軸方向）運(yùn)動(dòng)阻力系數(shù)，兩個(gè)方向的阻力試驗(yàn)結(jié)果見文獻(xiàn)[7]。

3 000 米和5 000米級(jí)海試的數(shù)據(jù)顯示縱傾角θ接近零度，而由上式（3）可知，當(dāng)縱傾角θ接近零度時(shí)攻角α可表示為θ的函數(shù)，采用數(shù)值算法可近似求出該函數(shù)，與文獻(xiàn)[7]類似，阻力項(xiàng)寫為：

0

2.2 浮力項(xiàng)B

隨著潛器所處深度不同，海水密度和壓力都會(huì)不同。一方面密度的變化會(huì)使?jié)撈鞲×Πl(fā)生相應(yīng)變化；另一方面潛器的各耐壓部件、密封部位和壓力補(bǔ)償部件的排水體積也會(huì)隨著海水壓力變化而變化。所以為了較準(zhǔn)確地計(jì)算潛水器在各個(gè)深度的浮力，需要采用對(duì)應(yīng)深度處的密度代入計(jì)算，而且需要較準(zhǔn)確地計(jì)算潛器排水體積的變化量，潛器的浮力項(xiàng)可寫為：

B=ρg （V-ΔV）

(5)其中：ρ=ρ （depth）為深度depth處的海水密度；ΔV為深度depth處潛器排水體積減少量，ΔV的變化分為兩個(gè)大的階段，第一階段是潛器壓載水箱注水或者排水階段，此階段ΔV變化速度較快；第二階段是潛器壓載水箱注滿水后的階段，此階段潛器的排水體積變化主要是由于海水壓力導(dǎo)致的，本文主要分析第二階段的ΔV；所以V為潛器的壓載水箱注滿水后的排水體積。

根據(jù)水池均衡實(shí)驗(yàn)結(jié)果[7]，潛器在水池平衡時(shí)總重20 807kg，而此時(shí)除了潛器本身的排水體積外，還有27kg鉛塊（密度11 700kg/m3）。此時(shí)潛水器在水池中平衡，有：

其中：Vm=潛器本體的排水體積。由（6）式可求得：

除了潛器本體，潛器攜帶的固定壓載、可拋壓載和作業(yè)工具等也產(chǎn)生一定的排水體積，這部分排水體積稱為可變排水體積Vc。則潛器剛潛入水面以下時(shí)的排水體積（未考慮體積變化量）為：

對(duì)于排水體積變化量ΔV，如前所述，在海水壓力下潛器排水體積發(fā)生變化的主要是耐壓結(jié)構(gòu)、密封部位和壓力補(bǔ)償部件，下面分兩部分進(jìn)行分析。

2.2.1 耐壓結(jié)構(gòu)排水體積變化

“蛟龍”號(hào)載人潛器上的主要耐壓結(jié)構(gòu)按外形可以分為兩類：

（1）球殼：載人艙、可調(diào)壓載水艙、高壓氣罐；

（2）圓柱殼：計(jì)算機(jī)罐、水聲通信機(jī)罐、測深側(cè)掃聲納罐、配電罐。

根據(jù)設(shè)計(jì)規(guī)范，可知潛器的耐壓結(jié)構(gòu)主體的材料在工作壓力下還屬于彈性范圍內(nèi)，所以可以采用理論公式計(jì)算排水體積變化量。

2.2.1.1 球殼耐壓結(jié)構(gòu)

表1 “蛟龍”號(hào)主要球殼耐壓結(jié)構(gòu)Tab.1 Large spherical pressure hull of “JiaoLong”

球殼在下潛開始時(shí)內(nèi)壓已經(jīng)存在，此時(shí)球殼的排水體積為：

其中Δr0為球殼徑向變形量：

（10）式中：P1=球罐的內(nèi)壓，Ri=球殼內(nèi)半徑，t=殼厚，E=材料楊氏模量，μ=材料泊松比，各參數(shù)值可見表1。

當(dāng)潛器下潛到某深度時(shí)，設(shè)海水壓力為P2，此時(shí)潛器上的各個(gè)球殼同時(shí)受外壓P2和內(nèi)壓P1，球殼的排水體積為：

其中Δr1為球殼徑向變形量：

則各個(gè)球殼的排水體積減小量為：

2.2.1.2 圓柱殼耐壓結(jié)構(gòu)

表2 “蛟龍”號(hào)主要圓柱殼耐壓結(jié)構(gòu)Tab.2 Large cylindrical pressure hull of “JiaoLong”

圓柱殼在下潛開始時(shí)內(nèi)壓已經(jīng)存在，此時(shí)球殼的排水體積為：

其中Δr0和Δl0分別是圓柱殼的徑向和軸向變形量：

（15）式和（16）式中：r2=圓柱殼外半徑；r1=球殼內(nèi)半徑；l為圓柱殼軸線長度；其他參數(shù)的定義同（10）式。各參數(shù)的取值見表2。

當(dāng)潛器下潛到某深度時(shí)，設(shè)海水壓力為P2，此時(shí)潛器上的各個(gè)圓柱殼同時(shí)受外壓P2和內(nèi)壓P1，圓柱殼的排水體積為：

其中Δr1和Δl1分別是圓柱殼的徑向和軸向變形量：

則各個(gè)圓柱殼的排水體積減小量為：

則整個(gè)潛水器的耐壓結(jié)構(gòu)排水體積減小量為：

需要注意的是，圓柱筒形耐壓結(jié)構(gòu)的封蓋可能是平封頭，也可能是半球封頭，且封頭上一般有穿艙件和密封結(jié)構(gòu)，所以封蓋的變形量不能用理論公式方便地求出，這部分的變形量將并入下一類。

2.2.2 其他排水體積變化

除了主要的耐壓結(jié)構(gòu)，潛器上存在大量的密封部件和壓力補(bǔ)償部件，例如載人艙觀察窗有機(jī)玻璃和觀察窗基座之間存在摩擦，需要通過有限元接觸分析才能計(jì)算觀察窗有機(jī)玻璃的變形量。又例如潛器的耐壓結(jié)構(gòu)上的封頭和耐壓結(jié)構(gòu)主體的裝備處往往采用O型圈等密封措施。在這些密封部件處，裝配體之間往往存在間隙，當(dāng)海水壓力小于某個(gè)值時(shí)，裝配部位通過密封圈接觸而裝備體互相之間不接觸，由于密封圈材料比較軟（通常密封圈材料為橡膠），此時(shí)變形量隨海水壓力增加而有較大的增加；當(dāng)海水壓力大于某個(gè)值時(shí)，密封圈變形量已經(jīng)足夠大，使得裝配體之間開始互相直接接觸，此后變形量隨海水壓力增加的速度明顯減小。這樣的現(xiàn)象在壓力筒試驗(yàn)中已經(jīng)觀察到了數(shù)次。例如模型球壓力筒試驗(yàn)中密封端蓋的變形量在約25MPa處存在明顯的拐點(diǎn)。此外潛器上存在的許多壓力補(bǔ)償?shù)牟考驮O(shè)備，例如充油電池箱、充油管路和電纜線等，這些設(shè)備和部件的內(nèi)部承壓部件往往外形復(fù)雜且包含多種材料。例如充油電池箱中承壓的是一塊塊的電池，而這些電池的外形比較復(fù)雜，電池的等效楊氏模量也未知，所以采用結(jié)構(gòu)分析來計(jì)算變形量變得不切實(shí)際。本文將所有不能通過球殼和圓柱殼理論公式計(jì)算變形量的部件的體積變化歸為“其他排水體積變化”，且假設(shè)這一變化量為潛器所受的海水壓力的函數(shù)，即：

則潛器總的體積變化量為：

2.3 重力項(xiàng)W

根據(jù)水池均衡試驗(yàn)和海試的經(jīng)驗(yàn)，潛器的重量可分為潛器本體重量和可變重量兩部分組成。

根據(jù)水池均衡試驗(yàn)，潛器在水池中均衡時(shí)重量為20 807kg[7]，此時(shí)除了潛器本體重量外，還包含有潛航員重量240kg、艙內(nèi)耗材（包括二氧化碳吸收劑）13kg和可調(diào)壓載水艙中140kg的淡水，則潛器本體重量為：

需要注意的是正常情況下潛器上安裝的鈦制采樣籃是不會(huì)變動(dòng)的，采樣籃將算作潛器本體的一部分，故M0中包括了20kg的鈦合金制采樣籃，而D44潛次中改為重量約25kg的鋼制采樣籃，所以將D44潛次中鋼制采樣籃比原裝鈦制采樣籃多出的5kg歸入作業(yè)工具。因此，在計(jì)算作業(yè)工具的排水體積時(shí)，需要扣除更換采樣籃帶來的排水體積變化。潛器的重量可以劃分為如表3所示的幾個(gè)部分。將除了M0外的潛器可變重量之和記為MC，則潛器下潛時(shí)的重量M為：

表3 重量組成Tab.3 Weight components

2.4 潛器無動(dòng)力潛浮運(yùn)動(dòng)平衡方程

潛器無動(dòng)力潛浮運(yùn)動(dòng)遵守牛頓第二定律，運(yùn)動(dòng)方程如（1）式所示，在文獻(xiàn)[7]中，潛器的下潛和上浮運(yùn)動(dòng)劃分為以下幾個(gè)階段：

（1）壓載水箱注水階段；

（2）平穩(wěn)下潛階段；

（3）拋掉P1后減速階段；

（4）潛器靠底階段；

（5）拋掉P2后加速階段；

（6）平穩(wěn)上浮階段；

（7）壓載水箱排水階段。

在文獻(xiàn)[9]中，對(duì)本文模型的7個(gè)階段均進(jìn)行了詳細(xì)的分析，出于篇幅考慮，此處只給出平穩(wěn)下潛階段和平穩(wěn)上浮階段的分析結(jié)果。

2.4.1 平穩(wěn)下潛階段

壓載水箱注水完成后潛器排水體積不再有大的變化 （除了由于海水壓力導(dǎo)致的排水體積減小），此階段潛器的運(yùn)動(dòng)也是動(dòng)態(tài)平衡的，即隨著重力與浮力阻力之差程減速—加速—減速—加速……的過程，但是通過分析3 000米海試和5 000米海試的數(shù)據(jù)可見每個(gè)減速—加速循環(huán)的周期不長，且潛器的速度隨著下潛深度的增加會(huì)逐漸降低，D44潛次從潛航員關(guān)閉壓載水箱上閥門到拋載P1的下潛過程的速度—深度曲線如圖2所示。

圖2 D44平穩(wěn)下潛階段速度—深度曲線Fig.2 The velocity curve of the steady descent stage of dive 44

為了簡化分析，本文將這一階段潛器的運(yùn)動(dòng)近似為平穩(wěn)運(yùn)動(dòng)，即認(rèn)為加速度項(xiàng)為0，加速度為零則不考慮附加質(zhì)量，此階段潛器的各項(xiàng)力為：

2.4.2 潛器平穩(wěn)上浮階段

潛器的重力、浮力和阻力達(dá)到動(dòng)態(tài)平衡后，潛器的上浮運(yùn)動(dòng)可看成平穩(wěn)運(yùn)動(dòng)，此階段前期的運(yùn)動(dòng)稱為平穩(wěn)上浮階段，此階段潛器的各項(xiàng)受力為：

則潛器這一階段的運(yùn)動(dòng)可以近似描述為平衡方程：

則潛器這一階段的運(yùn)動(dòng)可以近似描述為平衡方程：

3 D44潛次數(shù)據(jù)分析

D44潛次實(shí)現(xiàn)了完全的無動(dòng)力下潛上浮運(yùn)動(dòng)，其速度曲線相對(duì)沒有實(shí)現(xiàn)無動(dòng)力上浮下潛運(yùn)動(dòng)的潛次來說規(guī)律性比較明顯，為了確定已經(jīng)建立的無動(dòng)力運(yùn)動(dòng)模型中的未知參數(shù)，本文對(duì)D44潛次的配載情況、航行記錄和潛浮運(yùn)動(dòng)數(shù)據(jù)進(jìn)行了詳細(xì)的分析。

3.1 D44配載情況

D44潛次實(shí)現(xiàn)了無動(dòng)力上浮下潛運(yùn)動(dòng)，其配載較合理，根據(jù)水池均衡試驗(yàn)結(jié)果和D44潛次的記錄，D44潛次潛器的配載情況如下表4所示。

表4 D44配載情況Tab.4 Weight of dive 44

3.2 D44無動(dòng)力潛浮運(yùn)動(dòng)數(shù)據(jù)

與無動(dòng)力潛浮運(yùn)動(dòng)相關(guān)的D44潛次的海水密度—深度曲線、深度—時(shí)間曲線、速度—深度曲線如下圖3-5所示。

圖3 海水密度—深度曲線Fig.3 Density curve

圖4 D44下潛和上浮的深度—時(shí)間曲線Fig.4 The depth curve of submerge and ascent stage of dive 44

需要說明以下幾點(diǎn)：（a）從圖3中可看到深度大于1 000米后，海水密度隨深度的變化呈線性，所以本文經(jīng)將D44潛次的海水密度曲線外推至7 500米，且將該曲線和3 000米海試的海水密度曲線比較后發(fā)現(xiàn)，1000米以深的海水密度曲線差別較小。（b）圖4和圖5中的下潛曲線包括了從壓載水艙開始注水到拋掉P1的整個(gè)下潛過程；而上浮曲線包括了拋掉P2到壓載水艙開始排水的整個(gè)上浮過程。（c）與上浮至3 600米左右控制潛器上浮的操作記錄相對(duì)應(yīng)，可看到圖5中平穩(wěn)上浮階段的速度曲線在3 600米附近有一個(gè)“尖刺”。

而本文主要分析潛器下潛和上浮的平穩(wěn)運(yùn)動(dòng)階段，故將加速度值較大的曲線段去除，僅對(duì)平穩(wěn)階段的曲線進(jìn)行分析，由于潛器是動(dòng)態(tài)平衡的，所以速度曲線存在較多小幅的加速—減速振蕩，本文對(duì)D44平穩(wěn)下潛階段和平穩(wěn)上浮階段的速度曲線進(jìn)行了光順處理，以抓住平穩(wěn)運(yùn)動(dòng)的規(guī)律，D44平穩(wěn)下潛階段的速度曲線和光順后的曲線見圖6所示，同理D44平穩(wěn)上浮階段的速度曲線和光順后的曲線見圖7所示。

圖5 D44下潛和上浮的速度—深度曲線Fig.5 Velocity curve of submerge and ascent stage of dive dive 44

圖6 D44平穩(wěn)下潛階段速度—深度曲線及光順后曲線Fig.6 The smoothed velocity curve of the steady descent stage of dive 44

圖7 D44平穩(wěn)上浮階段速度—深度曲線及光順后曲線Fig.7 The smoothed velocity curve of the steady asscent stage of dive 44

4 模型參數(shù)辨識(shí)

4.1 平穩(wěn)下潛階段

由平穩(wěn)下潛階段的平衡方程式（26）可知，無動(dòng)力運(yùn)動(dòng)模型中的“其他排水體積變化量—ΔVother”是未知的，需要由5 000米海試D44潛次的平穩(wěn)下潛階段的曲線求得。D44下潛的ΔVother對(duì)海水壓力P的曲線見圖8所示，由光順后曲線可看到在約20-26MPa之間潛器的排水體積變化量存在明顯的拐點(diǎn)，可認(rèn)為潛器上大部分設(shè)備的密封部位和補(bǔ)償部位在20-26MPa時(shí)達(dá)到緊實(shí)，此后ΔVother隨P增長的速度降低（表現(xiàn)為圖上的曲線斜率降低）。

圖8 D44平穩(wěn)下潛階段ΔVother隨海水壓力P變化曲線Fig.8 The ΔVothercurve of the steady descent stage of dive 44

4.2 平穩(wěn)上浮階段

平穩(wěn)上浮階段的平衡方程見（27）式，ΔVother同樣是未知的，需要由5 000米海試D44潛次的平穩(wěn)上浮階段的曲線求得，需要注意的是由于壓縮量恢復(fù)有一定的滯后，上浮階段的ΔVother和下潛階段的ΔVother可能不同。ΔVother隨P變化的曲線見圖9，可見當(dāng)海水壓力大于6MPa時(shí)潛器ΔVother隨P的變化可近似為直線，而P在4.5MPa至6MPa之間曲線存在明顯的拐點(diǎn)，當(dāng)P小于4.5MPa時(shí)ΔVother減小的速度明顯加快，可認(rèn)為在4.5～6MPa之間潛器密封部件的密封面脫離閉合狀態(tài)，當(dāng)然，僅由D44潛次一個(gè)潛次的數(shù)據(jù)還不能下充分的結(jié)論。

圖9 D44平穩(wěn)上浮階段ΔVother隨海水壓力P變化曲線Fig.9 The ΔVothercurve of the steady ascent stage of dive 44

4.3 參數(shù)取值

將圖8和圖9中的D44下潛和上浮求出的ΔVother隨海水壓力P變化的光順后的曲線取出放在同一張圖中比較，如圖10所示?？梢姡弘S著下潛深度增大，下潛和上浮的ΔVother值逐漸接近，這個(gè)現(xiàn)象和實(shí)際情況較符合，因?yàn)镈44潛次最大潛深為5180米，在較深的深度處，例如5 000米處，潛器下潛和上浮至該深度處的壓縮量應(yīng)該差別不大，因?yàn)闈撈飨聺摮^5 000米深度只有180米，在最深的5 180米處作業(yè)接近2小時(shí)，這2個(gè)小時(shí)內(nèi)潛器相當(dāng)于處在保壓狀態(tài)，ΔVother將隨著作業(yè)時(shí)間增加而略有增加，但是從5 180米上浮至5 000米時(shí)的壓縮量和下潛至5 000米時(shí)的壓縮量不應(yīng)有大的差別。

圖10 D44航次下潛和上浮ΔVother隨P變化曲線（光順后）比較Fig.10 Comparison of smoothed ΔVothercurve of the steady descent stage and ascent stage of dive 44

為了將5 000米海試的結(jié)果外推至7 000米，需要對(duì)圖10中的曲線進(jìn)行分析和外插。由圖10，下潛的ΔVother-P曲線可近似為三段：

（1）P＜13MPa時(shí)，此階段潛器的ΔVother隨P的變化可近似為直線；

（2）13MPa＜P＜33MPa時(shí)，此階段潛器的ΔVother隨P的變化出現(xiàn)拐點(diǎn)，直接取由D44下潛數(shù)據(jù)反求出的光順后的曲線段；

（3）P>33MPa時(shí)，ΔVother隨P的變化也可近似為直線，將該直線段外插至76.7MPa（約7 500米），以便于估算7 000米處的ΔVother。

圖11 下潛ΔVother-P曲線的處理和外插Fig.11 Extrapolation of ΔVothercurve of the steady descent stage of dive 44

近似結(jié)果如上圖11所示。

類似的，上浮曲線也可分為三段：

（1）P＜3MPa時(shí)，此階段潛器的ΔVother隨P的變化可近似為直線；

（2）3MPa＜P＜7.6MPa時(shí)，此階段潛器的ΔVother隨P的變化出現(xiàn)拐點(diǎn)，直接取由D44上浮數(shù)據(jù)反求出的光順后的曲線段；

（3）P>7.6MPa時(shí)，ΔVother隨P的變化也可近似為直線，將該直線段外插至76.7MPa（約7 500米）。

上浮近似結(jié)果如下圖12所示。

將處理并外插后的下潛和上浮ΔVother-P曲線置于同一張圖中，如圖13所示，可看到下潛和上浮曲線的交點(diǎn)在約5 736米處，與最大下潛深度5 180米差556米，一方面可能是由于潛器在高壓下作業(yè)的過程中ΔVother會(huì)存在少量增加，另一方面也可能是由于上浮和下潛的數(shù)據(jù)存在誤差等導(dǎo)致。

圖12 上浮ΔVother-P曲線的處理和外插Fig.12 Extrapolation of ΔVothercurve of the steady ascent stage of dive 44

圖13 上浮和下潛的ΔVother-P曲線Fig.13 Used ΔVothercurve

可見，在前述模型的基礎(chǔ)上，為了更準(zhǔn)確地描述潛器上浮階段的運(yùn)動(dòng)，需要增加新的變量：作業(yè)階段（保壓階段）體積壓縮量ΔVwork—反應(yīng)潛器在底部作業(yè)過程中，即保壓狀態(tài)下排水體積減小量。對(duì)指定的潛器，ΔVwork是作業(yè)壓力和作業(yè)時(shí)間的函數(shù)，由于“蛟龍”號(hào)的作業(yè)時(shí)間基本固定，且還沒有足夠的數(shù)據(jù)支持，本文在此暫時(shí)假設(shè)ΔVwork只是作業(yè)壓力的線性函數(shù)：

則平穩(wěn)上浮階段平衡方程改為：

4.4 程序修改

根據(jù)以上計(jì)算模型和參數(shù)取值，對(duì)原有程序進(jìn)行了修改，界面也進(jìn)行了改進(jìn)，如圖14所示。計(jì)算輸入分為兩個(gè)部分，第一部分是預(yù)計(jì)作業(yè)深度Z0、上浮拋載、下潛拋載及密度曲線選擇；第二部分是潛器重量和排水情況。計(jì)算結(jié)果將輸出Z0處的均衡情況、下潛速度曲線和上浮速度曲線，同時(shí)也給出下潛/上浮的速度均值以及時(shí)間估算。

程序的大部分操作和原程序類似，但是存在以下不同：

（1）程序要求較準(zhǔn)確地輸入潛器各部分重量以及排水體積情況。注意“艙內(nèi)耗材”、“艙內(nèi)人員”和“壓載水箱注水”三項(xiàng)的排水體積為0；浮拋載和下潛拋載需要根據(jù)計(jì)算均衡結(jié)果人工進(jìn)行調(diào)整；

圖14 程序界面圖14 GUI of the program

（2）程序的“參考密度曲線”為5 000米海試得到的密度—深度曲線拓展至10 000米的密度曲線；“實(shí)測密度曲線”為下潛海區(qū)的實(shí)測密度—深度曲線，數(shù)據(jù)格式見程序目錄下的 “深度密度曲線.csv”文件，且要求“深度密度曲線.csv”中的第一個(gè)深度是0，最后一個(gè)深度大于程序輸入的Z0（大深度下密度-深度曲線可近似為直線，所以操作者可將實(shí)測密度-深度曲線的后半段 （例如對(duì)于Z0=7 000米，而實(shí)測曲線只有6 000米的數(shù)據(jù)，則操作者可以對(duì)3 000米以深的密度—深度曲線進(jìn)行線性擬合，操作者可以根據(jù)具體曲線形狀在excel中方便地進(jìn)行線性擬合和外插操作，本程序不再提供該功能）線性外插拓展至需要的深度 （一般可外插至10 000米），并將拓展后的曲線保存至 “深度密度曲線.csv”中）。

（3）程序的輸出調(diào)整為下潛深度Z0處的均衡情況。如果要求Z0處潛器達(dá)到均衡，則應(yīng)調(diào)整下潛拋載P2（因?yàn)闈撈飨聺撝罿0處時(shí)P1已經(jīng)被拋掉，所以P1不影響均衡位置）使結(jié)果接近于0；如果要求Z0處潛器具有足夠的坐底力，則應(yīng)調(diào)整下潛拋載P2使結(jié)果接近于設(shè)計(jì)坐底力（“蛟龍 ” 號(hào) 設(shè) 計(jì) 坐 底 力 ≈40）。

4.5 D44潛次計(jì)算結(jié)果

采用修改后的程序，將D44潛次的載荷數(shù)據(jù)填入后采用程序進(jìn)行計(jì)算，計(jì)算結(jié)果如圖15所示。

可見5 180處潛器重力浮力差（即坐底力）為35kg，與40kg比較接近，且速度曲線和圖6及圖7類似。

5 結(jié) 語

本文根據(jù)“蛟龍”號(hào)5 000米海試的D44潛次無動(dòng)力潛浮海試數(shù)據(jù)對(duì)“蛟龍”號(hào)無動(dòng)力潛浮運(yùn)動(dòng)模型進(jìn)行了再分析，改進(jìn)了3 000米海試后提出的分析模型，主要是采用了可變的體積壓縮率的概念，且對(duì)程序進(jìn)行了修改，開發(fā)出了更為可靠的蛟龍?zhí)栞d人潛器無動(dòng)力潛浮運(yùn)動(dòng)分析系統(tǒng)，使得計(jì)算結(jié)果更接近實(shí)際情況，程序界面更加符合實(shí)際應(yīng)用，為2012年的7 000米級(jí)海試做好了準(zhǔn)備，也為4 500米載人潛器無動(dòng)力潛浮運(yùn)動(dòng)的設(shè)計(jì)打下了基礎(chǔ)。

[1]Busby,Frank.Manned Submersibles[M].Arlington,VA:Office of the Oceanographer of the Navy,1976.

[2]Allmendinger,Eugene E.Submersible Vehicle Systems Design[M].ill.Jersey City,NJ:Society of Naval Architects and Marine Engineers,1990.

[3]朱繼懋主編.潛水器設(shè)計(jì)[M].上海:上海交通大學(xué)出版社,1992.

[4]崔維成,馬 嶺.潛水器設(shè)計(jì)中所要解決的水動(dòng)力學(xué)問題[C].第九屆全國水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨第二十二屆全國水動(dòng)力學(xué)研討會(huì)文集,吳有生,劉 樺,許唯臨,周連第,楊顯成主編,海洋出版社,北京,2009:9-29.

[5]劉正元.7000米載人潛器無動(dòng)力潛浮運(yùn)動(dòng)研究[R].無錫:中國船舶科學(xué)研究中心科技報(bào)告,2004.

[6]劉正元.潛水器大攻角范圍內(nèi)運(yùn)動(dòng)的仿真[J].船舶力學(xué),2005,9(2):54-59.Liu Zhengyuan.Simulation of submersible motion in large attack angle[J].Journal of Ship Mechanics,2005,9(2):54-59.

[7]劉正元,潘彬彬,崔維成.7000m載人潛器海試試驗(yàn)數(shù)據(jù)整理和分析[R].無錫:中國船舶科學(xué)研究中心科技報(bào)告,2010.

[8]范欣珊.壓力容器的應(yīng)力分析與強(qiáng)度設(shè)計(jì)[M].北京:原子能出版社,1979.

[9]潘彬彬,崔維成,葉 聰,劉正元.7000m載人潛器無動(dòng)力潛浮運(yùn)動(dòng)分析[R].無錫:中國船舶科學(xué)研究中心科技報(bào)告,2011.