董仁義，吳崇健

(中國(guó)艦船研究設(shè)計(jì)中心，湖北 武漢430064)

0 引 言

艦船管系在設(shè)計(jì)布置過(guò)程中由于受到艦、艇體內(nèi)部空間的限制，往往具有較為復(fù)雜的空間特性。具有復(fù)雜空間特性的液壓、疏水、海水冷卻等管系中流體流速較高、流量較大，在進(jìn)行管路流體截止和開(kāi)啟過(guò)程中，由于其改變流體的方向或速度，在慣性作用下，會(huì)產(chǎn)生明顯的瞬態(tài)振動(dòng)，這種流體瞬變導(dǎo)致的管道發(fā)生的劇烈強(qiáng)迫振動(dòng)，會(huì)造成管路與支撐結(jié)構(gòu)表面產(chǎn)生微動(dòng)磨損，使支撐剛度逐漸下降，從而降低管道的固有頻率，造成管路頻率與管系泵組的壓力脈動(dòng)頻率接近或者重合，產(chǎn)生流固耦合振動(dòng)，影響整個(gè)液壓系統(tǒng)的性能和降低設(shè)備的使用壽命，甚至可能導(dǎo)致艦船液壓系統(tǒng)功能失效;同時(shí)對(duì)于隱蔽性要求很高的潛艇，流體瞬變引起的管路振動(dòng)沖擊能量將通過(guò)管路支撐傳遞到耐壓殼體上，加劇殼體振動(dòng)產(chǎn)生的聲輻射，使?jié)撏У碾[蔽性遭受極大威脅。

因此，很有必要對(duì)具有復(fù)雜空間布置特性的艦船管路進(jìn)行流體瞬變導(dǎo)致的管系振動(dòng)分析，并且針對(duì)管系的振動(dòng)特性，給出空間管系的隔振措施。本文建立一個(gè)典型的具有復(fù)雜空間走向特性的艦船管系，分析在關(guān)閉管系閥件后造成的流體瞬變對(duì)管系振動(dòng)的影響，對(duì)管系的振動(dòng)特性進(jìn)行仿真。

1 管系流體瞬變

當(dāng)閥門(mén)突然關(guān)閉，管道中流速驟然降至0，閥前微元流段內(nèi)的水體首先受阻停止流動(dòng)并導(dǎo)致水壓升高，管壁膨脹，水體壓縮。近臨段處于上游側(cè)的水體則繼續(xù)流動(dòng)，并有水流入由于膨脹和水體壓縮而空出的這部分體積內(nèi)。空間填滿(mǎn)后，緊鄰著第2 個(gè)微元流段內(nèi)的水體又停止流動(dòng)，水壓升高并繼續(xù)壓縮，如此一層一層地向上游的水池方向傳去，形成逆?zhèn)鞯乃N升壓波動(dòng)。最后水體全部停止流動(dòng)。這時(shí)從管系末端開(kāi)始，降壓波按相反的方向傳播過(guò)來(lái)，直到抵達(dá)閥門(mén)處。如果忽略管道阻尼，這種升壓、降壓的過(guò)程將反復(fù)振動(dòng)進(jìn)行下去，如圖1所示［2］。

圖1 流體瞬變的壓力脈動(dòng)Fig.1 Pressure pulsation of transient flow

當(dāng)管系閥門(mén)突然關(guān)閉時(shí)，系統(tǒng)管路內(nèi)流體動(dòng)能會(huì)急劇變化。對(duì)于管系閥門(mén)上游管路，靠近閥門(mén)的一層厚度為dl 的流體流速突變?yōu)?，從而受到后面流體的壓縮，壓力驟增dp，管路也因受壓而膨脹，截面積增大。上游的流體逐步向下游流體過(guò)渡而轉(zhuǎn)為流體停止、管路受壓膨脹的狀態(tài)，形成以速度a向上游傳播的壓縮波。

這種流體壓力瞬變過(guò)程在工程上稱(chēng)為管系水錘現(xiàn)象。由于系統(tǒng)管路內(nèi)存在一定的沿程損失和局部損失，因此，水錘波在壓縮和膨脹過(guò)程中會(huì)逐步消耗能量。水錘波相長(zhǎng)的計(jì)算公式為T(mén)=2L/a?？梢?jiàn)，相同條件下，管路較長(zhǎng)的系統(tǒng)，水錘波能量消耗較慢，水錘幅值較大。水錘最大壓力升高可采用下式計(jì)算:

Δp=ρ·a·Δν。

式中:Δp 為管路最大壓力變化值;ρ 為流體密度;a為流體內(nèi)聲速;Δν 為閥關(guān)閉后流體流速的變化值。另外，流體內(nèi)聲速由下式來(lái)計(jì)算:

式中:ρ 為流體密度;K 為體積模量;E 為管材的彈性模量;D 為管路直徑;e 為管路壁厚。

2 管系仿真模型

2.1 管系描述

某段典型艦船空間管系的走向如圖2所示，管路公稱(chēng)通徑300 mm，壁厚為9.5 mm，最小屈服強(qiáng)度為206 N/mm2。A00 點(diǎn)和A19 點(diǎn)固定，流體假定為海水，密度為1 028 kg/m3，由A00 點(diǎn)流入，沿管路至A19 點(diǎn)，同時(shí)在圖2 中標(biāo)注出各管點(diǎn)的標(biāo)號(hào)和管路的長(zhǎng)度。在管系終點(diǎn)A19 處有一管路截止閥，A19 點(diǎn)截止閥突然關(guān)閉將產(chǎn)生沿管線(xiàn)傳遞的壓力沖擊波。

2.2 管系初步固定方式

圖3 給出了管系初步固定方式。在管點(diǎn)A02和A03 處用彈性隔震器與艙壁固定，設(shè)定彈性系數(shù)為119 N/mm，在管點(diǎn)A06 處設(shè)置導(dǎo)向支架，與管路上下左右的間隙為0，在管點(diǎn)A11，A16，A18 處設(shè)置支架，與管路上下的間隙為0。

圖2 典型的艦船空間管系Fig.2 Typical space pipeline

圖3 艦船空間管系的初步固定方式Fig.3 Basal fixed style of the space pipeline

3 流體瞬變過(guò)程中管系振動(dòng)特性仿真

文獻(xiàn)［3］進(jìn)行了艦船首尾移水系統(tǒng)水錘特性仿真與試驗(yàn)研究，通過(guò)對(duì)管路節(jié)點(diǎn)處流體壓力的分析，對(duì)關(guān)閥時(shí)間長(zhǎng)短與產(chǎn)生水錘的關(guān)系、水錘峰值產(chǎn)生的位置進(jìn)行討論。本文將從關(guān)閉閥時(shí)間以及管系的不同固定方式2 個(gè)方面，對(duì)流體瞬變過(guò)程中對(duì)管路振動(dòng)特性的影響進(jìn)行仿真。

3.1 閥門(mén)關(guān)閉時(shí)間對(duì)管路振動(dòng)特性的影響

按照初始給定的管路固定方式，設(shè)定不同的閥門(mén)關(guān)閉時(shí)間，對(duì)管路的振動(dòng)特性進(jìn)行對(duì)比分析。閥門(mén)關(guān)閉時(shí)間分別按照0.1 s，0.4 s，0.7 s，1 s 進(jìn)行設(shè)定。

前8 階振動(dòng)頻率如表1所示。

表1 管路系統(tǒng)模態(tài)頻率計(jì)算結(jié)果Tab.1 Modal frequency result of pipeline

不同閥門(mén)關(guān)閉時(shí)間對(duì)應(yīng)的管路振動(dòng)位移有很大不同，流體瞬變引起的管路位移如圖4所示，表2給出了各個(gè)管點(diǎn)在X，Y，Z 方向的水平位移，繞坐標(biāo)軸旋轉(zhuǎn)的角度可以忽略。

圖4 不同閥門(mén)關(guān)閉時(shí)間對(duì)應(yīng)的管路振動(dòng)位移Fig.4 Vibration displacement of pipeline under different time for closing valve

表2 管路振動(dòng)位移Tab.2 Vibration displacement of pipeline

由表2 給出的數(shù)據(jù)可以得出，離關(guān)閉閥門(mén)越近的管路，其由閥門(mén)關(guān)閉引起的管路振動(dòng)越強(qiáng)烈，并且管路整體振動(dòng)位移隨著閥門(mén)關(guān)閉時(shí)間的延長(zhǎng)而變小。

3.2 管系固定方式對(duì)管路振動(dòng)特性的影響

管系閥門(mén)關(guān)閉時(shí)所引起的流體瞬變會(huì)引起管路較大的位移變化，在閥門(mén)關(guān)閉時(shí)間t=0.1 s 工況，管系初步固定的方式下，增加管系的固定方式，對(duì)不同固定方式對(duì)管路振動(dòng)特性的影響進(jìn)行分析。圖5 給出了在初步固定方式下增加管系固定后的管路固定方式。

圖5 管系的固定方式Fig.5 Fixed style of pipeline

第1 種固定方式是在管系中管點(diǎn)A20，A21，A22 處增加了導(dǎo)向支架，與管路上下左右的間隙為0;第2 種固定方式是在管系中管點(diǎn)A20，A21，A22 處增加了導(dǎo)向支架，與管路上下的間隙為0。

3.2.1 第1 種固定方式下的管系動(dòng)力學(xué)參數(shù)分析

圖6 給出了第1 種固定方式下管路系統(tǒng)前6 階模態(tài)變形圖。

圖6 第1 種固定方式下管路系統(tǒng)前6 階模態(tài)變形圖Fig.6 Modal deformation figure of pipeline under the first fixed sytle

前8 階振動(dòng)頻率如表3所示。

表3 管路系統(tǒng)模態(tài)頻率計(jì)算結(jié)果Tab.3 Modal frequency result of pipeline

3.2.2 第2 種固定方式下的管系動(dòng)力學(xué)參數(shù)分析

圖7 給出了第2 種固定方式下管路系統(tǒng)前6 階模態(tài)變形圖。

圖7 第2 種固定方式下管路系統(tǒng)前6 階模態(tài)變形圖Fig.7 Modal deformation figure of pipeline under the second fixed style

圖8 兩種固定方式下的振動(dòng)位移Fig.8 Vibration displacement under two fixed style

前8 階振動(dòng)頻率如表4所示。

表4 管路系統(tǒng)模態(tài)頻率計(jì)算結(jié)果Tab.4 Modal frequency result of pipeline

表5 管路振動(dòng)位移Tab.5 Vibration displacement of pipeline

3.2.3 不同固定方式下的振動(dòng)位移

對(duì)于不同的固定方式，選取閥門(mén)關(guān)閉時(shí)間為0.1 s 的工況，分別計(jì)算2 種固定方式下的振動(dòng)位移，圖8 給出了振動(dòng)位移圖，表5 給出了管系中的管點(diǎn)A00 ～A19 以及新增固定點(diǎn)A20，A21，A22 處的水平位移，表中未給出管點(diǎn)繞坐標(biāo)軸的旋轉(zhuǎn)位移，旋轉(zhuǎn)位移可以忽略。

表5 中的振動(dòng)位移具體數(shù)據(jù)顯示，不同固定的方式中振動(dòng)位移變化的趨勢(shì)相同，對(duì)于固定自由度越多的管路約束，相比固定自由度小的管路約束，其振動(dòng)位移越小，同價(jià)振動(dòng)頻率降低。這與實(shí)際情況相符，可以證明對(duì)于管路振動(dòng)特性的仿真可以作為管路設(shè)計(jì)參考。

4 結(jié) 語(yǔ)

通過(guò)建立艦船管系的空間模型，設(shè)定管路的固定方式以及流體流動(dòng)方向、流量，制定管路中閥門(mén)的關(guān)閉時(shí)間，討論不同的閥門(mén)關(guān)閉時(shí)間所引起的流體瞬變導(dǎo)致的管路振動(dòng)特性，給出了空間管路的動(dòng)力學(xué)特性參數(shù)，并且詳細(xì)對(duì)比了管路節(jié)點(diǎn)在不同閥門(mén)關(guān)閉時(shí)間條件下的位移，得出閥門(mén)關(guān)閉時(shí)間與管路管點(diǎn)的平面振動(dòng)位移成反比，閥門(mén)關(guān)閉時(shí)間越短時(shí)，管路管點(diǎn)的振動(dòng)位移越大，閥門(mén)關(guān)閉時(shí)間越長(zhǎng)時(shí)，管路管點(diǎn)振動(dòng)位移越小，并且越靠近閥門(mén)關(guān)閉點(diǎn)的管線(xiàn)振動(dòng)位移越大，振動(dòng)位移的峰值通常出現(xiàn)在離閥門(mén)較近的地方。

同時(shí)，在相同的閥門(mén)關(guān)閉時(shí)間工況下，管路在流體瞬變過(guò)程中，對(duì)不同固定方式，空間管路的振動(dòng)位移變化趨勢(shì)相同，管路振動(dòng)位移與固定方式限制的管路自由度有關(guān)，這對(duì)于設(shè)計(jì)艦船管路空間管路過(guò)程中，考慮布置所采取的固定方式和固定位置有一定的指導(dǎo)意義。

［1］董仁義，吳崇健，張京偉，等．基于Modelica/MWorks 的艦船液壓操舵系統(tǒng)建模與仿真研究［J］．艦船科學(xué)技術(shù)，2011，33(11):76－80．

DONG Ren-yi，WU Chong-jian，ZHANG Jing-wei，et al．Modeling and simulation for ship hydraulic rudder system based on Modelica/MWorks［J］．Ship Science and Technology，2011，33(11):76－80．

［2］譚平，徐蕾，凌曉聰．動(dòng)力管道水錘激振分析［J］．南京理工大學(xué)學(xué)報(bào)，2006(30):2，182－185．

TAN Ping，XU Lei，LING Xiao-cong，et al．Vibration of waterhammer of power pipe system［J］．Journal of Nanjing University of Science and Technology，2006，30(2):182－185．

［3］蔡標(biāo)華．艦船首尾移水系統(tǒng)水錘特性仿真與試驗(yàn)［J］．艦船科學(xué)技術(shù)，2011，33(9):52－55．

CAI Biao-hua．The water hammer characteristic simulation and test study in warship system pipeline［J］．Ship Science and Technology，2011，33(9):52－55．