林舟曦 何世平

（1.海軍工程大學動力工程學院 武漢 430033）（2.湖北工業(yè)大學機械工程學院 武漢 430068）

1 引言

關于浮體水中上浮的相關動力學模型國內(nèi)外已經(jīng)做了大量研究，也產(chǎn)生了很多設計應用，例如黃祥兵、黃興玲、董云飛[1]對水下航行器集體逃生艙上浮過程進行了分析，建立逃生艙上浮模型并用數(shù)學方法和仿真軟件驗證了模型的可靠性。孫石、宋兆麗[2]分析了水球在水中的受力情況和能量變化，引入浮力勢能概念，給出了水球上浮的理論模型。李志偉、崔維成[3]分析了深海載人潛水器的下潛上浮過程，推算出深海潛水器運動阻力估算公式。潘彬彬等[4]引入收縮率概念改進潛水器上浮與下潛的理論模型，并開發(fā)出了無動力潛浮過程的分析系統(tǒng)。徐萌萌[5]等對一種浮體自由上浮運動過程進行測量與研究，建立了浮體自由上浮的速度，位移與時間的理論模型與仿真結(jié)果并通過實驗驗證。Nur Afande AliHussain[6]等對水下航行器上浮下潛過程進行分析，建立了水下航行器平衡浮力數(shù)學模型。Michael Woolwright[7]設計了一種水下釋放裝置，用于救生筏的水下釋放。Wainman[8]采用類似的滾筒、繩索以及浮體的機械結(jié)構(gòu)，設計了一種應急浮標的水下釋放裝置等。

而對于攜帶繩子和滾筒的浮體上浮模型，在許多領域也可以見到。比如水雷的海底布置過程[9]，水下航行器應急通信浮標釋放，以及水下航行器上救生筏釋放的上浮過程等物理過程。這些物理過程基本上在目前已經(jīng)實現(xiàn)或是仍在研究中，但是適用于這些物理過程的理論模型以及相關的研究較少，雖然文獻[5]對浮體上浮同樣給出速度位移與時間的數(shù)學模型，但模型不考慮繩子和滾筒。因此，本文進行了進一步的研究，即對存在繩子和固定滾筒的浮體上浮過程進行研究，建立其理論模型，通過仿真，繪制出攜帶繩子和固定滾筒的浮體上浮過程速度，位移與時間的函數(shù)關系曲線圖。最后根據(jù)實際救生筏水下釋放實驗，驗證所建立的理論模型以及仿真結(jié)果具有一定的準確性，為這具有一定普遍性、多領域應用性的物理過程提供一個理論參考。

2 帶繩子滾筒的浮體上浮過程動力學分析

因為實際實驗的浮體近似球形，為了接近實驗的實際情況，浮體簡化為球體，建立模型后如圖1所示。

圖1 浮體、繩子、滾筒模型

對浮體模型作受力分析，如圖2所示。

有些情況繩子的重量可以忽略不計，例如救生筏繩索系統(tǒng)，繩子輕。而有些情況繩子的重量不可忽略不計，例如水雷的水下釋放，因為錨鏈的重量比較大，錨雷在上浮過程中，錨鏈的重量就不能忽略，但是這里本文統(tǒng)一考慮繩子重量，并用m（t）表示，讀者可根據(jù)實際繩索情況確定繩子模型的線密度，從而確定繩子質(zhì)量。

圖2 受力分析

所以浮體上升的總質(zhì)量為浮體質(zhì)量與繩子質(zhì)量之和：

環(huán)境是海水，因此浮體在海水中受到海水的流體阻力，因此：

而增加的質(zhì)量與繩子的線密度以及從滾筒輸送出來的長度有關，得：

式中：F浮為浮體受到的浮力，單位為N；M（t）為t時刻，浮體和繩子的總質(zhì)量，單位為kg；M浮體為球形的浮體、救生筏以及相關設備的質(zhì)量，單位為kg；F阻為海水阻力，單位為N；F拉為繩子拉力，單位為N；ps為海水密度，單位為kg/m3；pl為繩子線密度單位為kg/m；L（t）為t時刻繩子的長度，單位為m；c為水中阻力系數(shù)；A為垂直于來流方向的物體橫截面積，單位為m2；vt為t時刻浮體上升的速度，單位為m/s。

對球形浮體上浮過程中的滾筒作受力分析，因為本次實驗存在滾筒的制動裝置，因此存在制動力矩，另外還需考慮滾筒中間的軸零件與其他零件相對運動產(chǎn)生的摩擦，因此還存在軸的摩擦力矩，但是制動裝置制動也是采用摩擦力制動原理，并且兩者都是常數(shù)，因此這里進行簡化（制動裝置忽略不畫），統(tǒng)一用摩擦力矩表示，符號M摩，如圖3所示。

圖3 受力分析

因為繩子物理模型的拉伸量可忽略不計，所以浮體位移距離和滾筒外邊緣釋放的繩子長度，即滾筒外邊緣的總線位移相同，即：

另外：

所以可得到：

式中：s總為浮體上浮位移，單位為m；α總為浮體上浮過程中，滾筒轉(zhuǎn)過的總弧度，單位為rad；r為滾筒半徑，單位為m；a為浮體上浮加速度，單位為m/s2；α為浮體上浮過程中，滾筒轉(zhuǎn)動的角加速度，單位為rad/s2；M摩為滾筒滾動時受到的摩擦力矩，單位為N·m。J為轉(zhuǎn)動慣量，單位為kg·m2。

聯(lián)立式（1）～式（9），得到：

3 上浮過程仿真模擬及分析

3.1 仿真模擬圖像

為了繪制浮體上浮的速度、加速度、位移以及繩子拉力和時間的關系曲線，使用實際測量的數(shù)據(jù)：M浮體=200kg，F(xiàn)浮=320N，g=9.8N/kg，J=0.6kg·m2，r=0.045m，ps=1025kg/m3，pl=0.1kg/m（鋼絲繩），c=0.5，A=0.6，M摩=3N·m。并帶入式（10），通過仿真計算，繪制曲線圖如下。

圖4 速度與時間的關系

圖5 上升位移與時間的關系

3.2 仿真模擬結(jié)果分析

從位移圖來看，浮體上浮過程總體比較接近直線，說明浮體上浮狀態(tài)接近勻速運動。但從速度圖像來看，速度先從零上升到最大速度，再漸漸下降，而文獻[5]仿真的速度不會有明顯下降，初步說明繩子的存在對浮體上浮有一定影響。

4 浮體上浮過程實驗與分析

4.1 上浮實驗數(shù)據(jù)采集

以救生筏水下釋放為實驗背景，氣脹式救生筏一般用于水面船只上，并從水面釋放，但是為了能夠讓其應用在水下航行器上，釋放到水面后為遭遇事故后的水下航行器落水人員提供庇護所[10～12]，本課題組制作出工程樣機進行水下釋放實驗。本次實驗中，救生筏安裝在球形的浮體中，實驗時，在滾筒上每250mm的鋼絲繩處，綁上紅色的防水膠，鋼絲繩與浮體用尼龍繩綁定，白色尼龍繩上用藍色油性筆每隔150mm做上記號。實驗過程中用水下攝像機對準滾筒處進行錄像（圖6）。

圖6 水下攝像機畫面

本次總共進行了水下10m和水下25m的實驗，釋放實驗結(jié)束后，對滾筒錄像視頻進行數(shù)據(jù)采集，計算出浮體上浮速度、位移與時間的關系，繪制曲線圖（圖7和圖8）。實驗時，浮體上浮到水面，因此，速度最終降至零點。

圖7 速度與時間的關系

圖8 位移與時間的關系

4.2 實驗結(jié)果與理論的比較分析

將實驗曲線與理論曲線進行比較（圖9與圖10）。

圖9 速度比較

圖10 位移比較

對比結(jié)果后總結(jié)如下：

采樣的精度產(chǎn)生誤差，本次實驗中因為攝像機角度的原因，雖然同樣是25m的長度，但是離攝像機越遠，同樣的長度在畫面上所顯現(xiàn)的距離就會越短，會給測量估計造成困難，因此在數(shù)據(jù)采樣時會產(chǎn)生測量誤差。

浮體離開釋放裝置后所達到的最大速度比較接近，都在2.5m/s～2.8m/s范圍內(nèi)。因此最大速度方面仿真模擬結(jié)果接近實際。

從速度趨勢上來說，通過25m的實驗曲線和仿真曲線的比較，可以發(fā)現(xiàn)，浮體的速度在到達最大速度后，都開始逐漸下降（文獻5中實驗結(jié)果速度沒有明顯下降）。說明浮體上浮過程中繩子的存在確實對浮體上浮產(chǎn)生了影響。因此速度趨勢方面，所建立的模型與實際情況接近。

速度降低形式上，從25m釋放的實驗曲線結(jié)果可以看出，與文獻[5]類似，速度產(chǎn)生了波動，可能原因是所建立的模型與本次實際實驗有所差別。如圖11和圖12，在所建立的模型中，浮體從開始的釋放到上浮整個過程中，繩子與球形浮體的連接點始終被認為處于浮體的正下方，然而實際在進行水下實驗時，尼龍繩的一端是綁在浮體側(cè)面上的，也就是說，實際實驗時，浮體開始釋放瞬間，會產(chǎn)生繞繩子與浮體連接點的翻轉(zhuǎn)并很可能在之后的過程中形成周期性的擺動現(xiàn)象（圖13），而且這種翻轉(zhuǎn)或者擺動幅度現(xiàn)象可能在浮體上浮過程中，因為受到旋轉(zhuǎn)周圍的海水阻力而逐漸減?。愃茩C械阻尼結(jié)構(gòu)的效果）。因此呈現(xiàn)出繩子與浮體的連接點產(chǎn)生速度起伏并且幅度降低的現(xiàn)象。此現(xiàn)象需要在今后試驗中改進所建立的模型并繼續(xù)研究。

圖11 想象的安裝情況

圖12 實際安裝情況

圖13 浮體翻轉(zhuǎn)

5 結(jié)語

1）實驗結(jié)果浮體上浮過程所達到的最大速度與理論最大速度接近。

2）繩子的存在對浮體的上浮速度產(chǎn)生影響。

3）實驗結(jié)果中，浮體上浮過程速度總體上逐漸降低的趨勢特點與仿真結(jié)果相符合。

4）實驗模型具有一定的準確性，但存在誤差，并且所建立的模型需要改進。

5）仿真以及實驗結(jié)果可以為本領域或是其他領域的同類物理過程的工程問題提供參考。