付建明 王 濤 于廣錄 王 冰 閆 俊

（1.遼寧紅沿河核電有限公司，遼寧 大連 116001；2.大連鴻澤泵業(yè)有限公司，遼寧 大連 116036；3. 大連理工大學(xué)，海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024）

0 引言

由于大量海生物的出現(xiàn)，核電站經(jīng)常發(fā)生海生物和漂浮物進(jìn)入冷源入水口的事件，從而造成核電廠冷源系統(tǒng)堵塞[1-2]。目前，核電站采用傳統(tǒng)“網(wǎng)兜攔截+人工海上漁船駁運(yùn)清理”的方式來防護(hù)冷源系統(tǒng)，這種方式效率低且夜間作業(yè)風(fēng)險(xiǎn)高，會(huì)直接對(duì)核電站安全運(yùn)行造成影響[3]。因此，遼寧某公司和大連某公司聯(lián)合設(shè)計(jì)、制造了一種冷源取水口攔污網(wǎng)網(wǎng)兜海生物抽吸裝置，將其安裝在工作船上，將被攔截入網(wǎng)兜的海生物和漂浮物輸送到海面，并運(yùn)送至指定地點(diǎn)，代替?zhèn)鹘y(tǒng)人工海上漁船駁運(yùn)清理冷源入水口攔截網(wǎng)兜中海生物的作業(yè)方式，該裝置可以在設(shè)計(jì)工況下全天候、無間斷地開展清理作業(yè)，徹底解決人工海上漁船駁運(yùn)清理存在的安全隱患，提高經(jīng)濟(jì)效益，保障冷源安全。冷源取水口攔污網(wǎng)網(wǎng)兜抽吸裝置工作船在進(jìn)行工作時(shí)，啟動(dòng)錨機(jī)系統(tǒng)，拋錨索固定，將海生物切割裝置調(diào)整至工作狀態(tài)，連接網(wǎng)兜，即可在操作臺(tái)上按要求模式進(jìn)行抽吸作業(yè)。因此，需要根據(jù)具體工作的海洋環(huán)境條件對(duì)工作船的運(yùn)動(dòng)響應(yīng)進(jìn)行分析，獲取關(guān)鍵位置的運(yùn)動(dòng)響應(yīng)結(jié)果，評(píng)估其抽吸工作的安全性能。

1 工作船設(shè)計(jì)需求

1.1 工作海況條件

根據(jù)船體的對(duì)稱性，考慮風(fēng)、浪和流環(huán)境荷載入射方向與船長方向的夾角為0°～180°（風(fēng)、浪和流共向），如圖1所示。選取2 級(jí)海況波浪，3 級(jí)風(fēng)風(fēng)速，具體的環(huán)境參數(shù)值見表1。

圖1 工作船環(huán)境荷載入射方向

表1 作業(yè)工況環(huán)境參數(shù)

1.2 運(yùn)動(dòng)響應(yīng)限制

工作船結(jié)構(gòu)示意圖如圖2 所示。工作船能夠承載提升裝置、輸送泵系統(tǒng)、滿負(fù)荷為1 m3的過濾箱、錨機(jī)以及就地控制柜等設(shè)備和4 名操作人員，工作臺(tái)設(shè)計(jì)有4 個(gè)錨絞機(jī)和導(dǎo)纜孔，以保證工作船水下連接網(wǎng)兜位置處的耦合裝置的漂移范圍不超過1.5 m。

圖2 工作船結(jié)構(gòu)示意圖

2 運(yùn)動(dòng)響應(yīng)分析模型

根據(jù)工作船的運(yùn)行情況，建立其對(duì)應(yīng)的簡化運(yùn)動(dòng)響應(yīng)分析模型，如圖3 所示。工作船會(huì)受到風(fēng)、浪以及流環(huán)境的聯(lián)合作用。

圖3 工作船運(yùn)動(dòng)響應(yīng)分析模型

2.1 環(huán)境荷載計(jì)算

2.1.1 波浪荷載

采用勢流理論對(duì)規(guī)則波作用下浮式結(jié)構(gòu)物的輻射和繞射進(jìn)行分析。假設(shè)浮式結(jié)構(gòu)物所受的一階波浪力由波浪激勵(lì)力和輻射力組成，二者可以線性疊加。其中，波浪激勵(lì)力可看作是由入射力和繞射力組成的，輻射力是由浮式結(jié)構(gòu)物運(yùn)動(dòng)引起的流體反作用力，由附加質(zhì)量力和輻射阻尼力組成。

對(duì)受到簡諧波浪作用的浮式結(jié)構(gòu)物來說，其周圍流場速度勢可以分解為入射勢、輻射勢和繞射勢3 個(gè)部分[4]，如公式（1）所示。

2.1.2 風(fēng)荷載

對(duì)非超大型原油運(yùn)輸船（Very Large Crude Carrier，VLCC）船型浮式結(jié)構(gòu)來說，可以通過模塊法來計(jì)算風(fēng)力和風(fēng)力矩。根據(jù)API RP 2SK 規(guī)范[5]，風(fēng)力計(jì)算如公式（2）所示。

式中：FW為風(fēng)力；CW為風(fēng)力系數(shù)；CS為構(gòu)件形狀系數(shù)；Ch為構(gòu)件高度系數(shù)；A為構(gòu)件迎風(fēng)面的垂向投影面積；VW為設(shè)計(jì)風(fēng)速。

2.1.3 流荷載

對(duì)非VLCC 船型浮式結(jié)構(gòu)來說，可以根據(jù)模塊法來計(jì)算流力和力矩。根據(jù)API RP 2SK 規(guī)范[5]，流力計(jì)算如公式（3）所示。

式中：FC為流力；CC為流力系數(shù)；Cd為拖曳力系數(shù)，對(duì)圓形單元取0.5，對(duì)平面單元，可參考圖 4 取值；Ac為水線面下圓柱形構(gòu)件的投影面積總和；Af為水線面下平板構(gòu)件的投影面積；Vc為來流速度。

圖4 平面單元的流力系數(shù)

2.2 時(shí)域耦合運(yùn)動(dòng)方程

浮式結(jié)構(gòu)物在海上作業(yè)時(shí)會(huì)受到各種外界力的作用，通常包括風(fēng)、浪以及流等所產(chǎn)生的海洋環(huán)境載荷和定位系統(tǒng)所產(chǎn)生的系泊力等，會(huì)導(dǎo)致浮式結(jié)構(gòu)物產(chǎn)生6 個(gè)自由度的運(yùn)動(dòng)。浮式結(jié)構(gòu)及其系泊系統(tǒng)總體運(yùn)動(dòng)如公式（4）所示[6]。

式中：x˙˙、x˙和x分別為浮式結(jié)構(gòu)物的6 個(gè)自由度運(yùn)動(dòng)的加速度、速度和位移；M為浮式結(jié)構(gòu)質(zhì)量矩陣；Ma為附加質(zhì)量矩陣；C為線性阻尼矩陣；K為靜水剛度矩陣；Fws為作用在浮式結(jié)構(gòu)物上的平均波浪漂移力；Fwf為波頻力；Fsd為波浪低頻慢漂力；Fwind為風(fēng)力；Fcurent為流力；Fmooring為系泊力。

3 計(jì)算結(jié)果與分析

3.1 水動(dòng)力模型

首先，利用船體建造組立圖，通過三維建模軟件CATIA 建立工作船的三維質(zhì)量模型，包括船體外殼、船體內(nèi)部艙壁、船體內(nèi)部型材以及船體駕駛室等，如圖5 所示。計(jì)算船體重心、質(zhì)量以及轉(zhuǎn)動(dòng)慣量等參數(shù)，結(jié)果見表2，其中，Gx、Gy和Gz分別為工作船重心的x坐標(biāo)、y坐標(biāo)和z坐標(biāo)，坐標(biāo)原點(diǎn)設(shè)置在船艉底部船中處，Ixx、Iyy和Izz分別為工作船對(duì)x軸、y軸和z軸的轉(zhuǎn)動(dòng)慣量。

表2 清潔船質(zhì)量、重心和轉(zhuǎn)動(dòng)慣量

圖5 工作船三維質(zhì)量模型

其次，采用ANSYS APDL 建立船體的幾何外形并剖分網(wǎng)格，得到船體的三維水動(dòng)力模型，如圖6 所示。利用宏命令A(yù)NSTOAQWA將ANSYS 中的面元模型導(dǎo)為AQWA-LINE 的計(jì)算文件，采用AQWA-LINE 計(jì)算該船體的水動(dòng)力參數(shù)。選取周期分布為3 s～20 s 的單位振幅的規(guī)則波作為荷載，考慮船體的對(duì)稱性，計(jì)算不同波浪入射角度（與X軸的夾角為0°、45°、90°、135°和180°）時(shí)該船體的水動(dòng)力系數(shù)，部分水動(dòng)力結(jié)果如圖7 所示。其中，圖7（a）～圖7（c）中繞射力包括不同頻率單位振幅規(guī)則波入射勢和繞射勢引起的作用力，圖7（d）～圖7（f）中RAO 表示不同頻率單位振幅規(guī)則波引起的浮體運(yùn)動(dòng)響應(yīng)幅值。

圖6 工作船三維水動(dòng)力模型

圖7 工作船水動(dòng)力系數(shù)

3.2 時(shí)域響應(yīng)

考慮表1 中的環(huán)境荷載以及圖1 中的不同入射方向，分析工作船重心位置縱蕩、橫蕩、垂蕩、縱搖、橫搖和艏搖6 個(gè)自由度位移響應(yīng)。同時(shí)，分析工作船水下連接網(wǎng)兜位置處的位移響應(yīng)，運(yùn)動(dòng)響應(yīng)統(tǒng)計(jì)值見表3（0°方向環(huán)境荷載作用的結(jié)果如圖8 所示）。

從圖8 和表3 可知，工作船在系泊作業(yè)工況且入射角度不同的特定環(huán)境荷載作用下，其縱蕩自由度最大響應(yīng)發(fā)生在環(huán)境荷載0°入射時(shí)，最大值為0.413 m；橫蕩自由度最大響應(yīng)發(fā)生在環(huán)境荷載90°入射時(shí)，最大值為-0.433 m；垂蕩自由度最大響應(yīng)發(fā)生在環(huán)境荷載90°入射時(shí)，最大值為0.319 m；橫搖自由度最大響應(yīng)發(fā)生在環(huán)境荷載90°入射時(shí)，最大值為4.19°；縱搖自由度最大響應(yīng)發(fā)生在環(huán)境荷載0°入射時(shí)，最大值為2.33°；艏搖自由度最大響應(yīng)發(fā)生在環(huán)境荷載135°入射時(shí)，最大值為0.953°；環(huán)境荷載0°入射時(shí)，網(wǎng)口X向位移最大，為0.572 m；環(huán)境荷載90°入射時(shí)，網(wǎng)口Y向位移最大，為-0.653 m；網(wǎng)口Z向最大位移為0.488 m，發(fā)生在環(huán)境荷載90°入射時(shí)。

圖8 工作船位移響應(yīng)（0 度方向環(huán)境荷載作用）

表3 作業(yè)工況的響應(yīng)統(tǒng)計(jì)值

4 結(jié)論

對(duì)核電站冷源取水口攔污網(wǎng)網(wǎng)兜抽吸裝置工作船在系泊作業(yè)工況條件下的運(yùn)動(dòng)響應(yīng)進(jìn)行數(shù)值分析可以得出以下2個(gè)結(jié)論：1）建立的工作船運(yùn)動(dòng)響應(yīng)計(jì)算模型可以準(zhǔn)確計(jì)算整體裝置的響應(yīng)并評(píng)估其設(shè)計(jì)安全。2）在給定的環(huán)境工況條件下，設(shè)計(jì)的工作船系泊狀態(tài)下能夠滿足水下連接網(wǎng)兜位置處的運(yùn)動(dòng)響應(yīng)約束限制，保證整體裝置良好地運(yùn)行。