絞吸式挖泥船絞刀的流場(chǎng)數(shù)值模擬*

周棟彬 熊 庭 張新卓 江 攀

(武漢理工大學(xué)能源與動(dòng)力工程學(xué)院 武漢 430063)

0 引 言

絞吸式挖泥船在港灣建設(shè)、航道湖泊清淤、填海造地、水域環(huán)保等領(lǐng)域有重要應(yīng)用，安裝于其橋架前端的絞刀是絞吸式挖泥船的主要工作部件，絞刀旋轉(zhuǎn)時(shí)由刀臂帶動(dòng)刀齒進(jìn)行破碎切削工作.切削得到的疏浚泥漿的主要成分以淤泥和黏土為主，具有固體顆粒粒度小，分布范圍寬，含水量高的特點(diǎn)[1].絞刀的流場(chǎng)特性對(duì)挖掘性能有重要影響，而絞刀的挖掘性能是影響挖泥船挖泥效率的關(guān)鍵因素.在絞刀工作時(shí)，刀齒及刀臂和介質(zhì)直接接觸，刀齒安裝于齒座上，當(dāng)?shù)洱X出現(xiàn)磨損時(shí)可即時(shí)更換，刀臂則無(wú)法更換，而刀臂為絞刀的主要連接及支撐結(jié)構(gòu)，以往的相關(guān)研究對(duì)刀臂上的壓力變化情況的關(guān)注較少[2-3]，對(duì)絞刀的流場(chǎng)分析一般以轉(zhuǎn)速為主要研究對(duì)象且多為簡(jiǎn)化的清水流場(chǎng)[4-5]. 鑒于此，文中采用流場(chǎng)仿真軟件fluent建立了絞刀的三維有限元模型，將絞刀工作時(shí)的下傾角及入口流速作為研究對(duì)象，采用數(shù)值模擬方法研究絞刀在最佳轉(zhuǎn)速工況下，不同下傾角和所處水域流速對(duì)絞刀流場(chǎng)的壓力分布和出口泥沙體積分?jǐn)?shù)的影響，得到減小絞刀所受壓力并提高出口泥沙體積分?jǐn)?shù)的工況參數(shù).

1 絞刀流場(chǎng)的數(shù)值仿真模型

1.1 絞刀三維造型

絞刀主體結(jié)構(gòu)由大環(huán)、輪轂、刀臂、刀齒組成，輪轂與絞刀軸連接，傳遞力和力矩，刀齒附著在刀臂上，大環(huán)和輪轂把刀臂連接起來(lái)，實(shí)現(xiàn)連續(xù)切削過(guò)程.絞刀的設(shè)計(jì)尺寸對(duì)切削性能有重要影響，而刀臂的幾何形狀是絞刀設(shè)計(jì)中的關(guān)鍵參數(shù)[6].數(shù)值仿真中使用的鉸刀三維模型包括六個(gè)刀臂，48個(gè)刀齒，刀臂由內(nèi)外輪廓曲線及截面生成空間曲面.絞刀尺寸見(jiàn)表1,三維模型見(jiàn)圖1.刀臂的內(nèi)外輪廓曲線為

(1)

(2)

式中：d1為輪轂外徑；Ω1為刀臂外輪廓線包角；Ω2為刀臂內(nèi)輪廓線包角；D為絞刀直徑；k為絞刀形狀系數(shù)，與刀臂外展角有關(guān)；N為土質(zhì)的標(biāo)準(zhǔn)貫入擊數(shù)；H為刀臂外輪廓線頂點(diǎn)高度；h為刀臂內(nèi)輪廓線頂點(diǎn)高度；φ為刀臂安裝角.

表1 絞刀尺寸參數(shù)

1-大環(huán)；2-刀臂；3-刀齒；4-輪轂圖1 絞刀三維造型

1.2 絞刀流場(chǎng)分析的數(shù)學(xué)模型

絞吸式挖泥船在進(jìn)行切削工作時(shí)，橋架將絞刀下放至待挖土層，由電機(jī)驅(qū)動(dòng)絞刀對(duì)土質(zhì)進(jìn)行擠壓破碎，并隨著橋架繞定位樁旋轉(zhuǎn)，由于絞刀直徑相比繞定位樁旋轉(zhuǎn)的直徑較小，因此，可以近似于直線橫移運(yùn)動(dòng).土質(zhì)剝離破碎后一部分?jǐn)U散至周?chē)蚧蛘叱练e在水底，另一部分在泵吸作用下經(jīng)由吸口進(jìn)入吸泥管，這部分的土質(zhì)就是挖泥船的最終產(chǎn)量，絞刀的切削原理見(jiàn)圖2.

圖2 絞刀切削原理圖

計(jì)算使用的Mixture模型是由離散型和連續(xù)相構(gòu)成的多相流模型，允許相之間互相滲透，它求解的是各相的連續(xù)方程和混合相的動(dòng)量方程，與拉格朗日模型相比，mixture模型考慮的變量更少，計(jì)算更有優(yōu)勢(shì).其控制方程描述如下.

連續(xù)方程：

(3)

式中：ρm為第m相密度，kg/m3；vm為第m相矢量速度；φm為m相的體積分?jǐn)?shù).

Mixture模型的動(dòng)量方程為

(4)

Mixture模型的能量方程為

(5)

計(jì)算的湍流模型使用k-ε方程：

(6)

(7)

式中：ρ為水的密度；ui為速度矢量；μt為湍流黏度；Gk為平均速度梯度引起的湍流k的生產(chǎn)項(xiàng)；sk=1.0，sε=1.3，s1ε=1.44，s2ε=1.92.

2 絞刀流場(chǎng)數(shù)值模擬

2.1 流場(chǎng)計(jì)算域的設(shè)置和網(wǎng)格劃分

計(jì)算中使用的鉸刀直徑為1 750 mm、高度為1 125 mm，在geometry中將鉸刀所在流場(chǎng)域劃分為旋轉(zhuǎn)域和靜止域.旋轉(zhuǎn)域設(shè)置為包圍鉸刀的圓柱，直徑1.85 m、高1.325 m，實(shí)際施工中下傾角的一般大于20°[7]，因此，將下傾角分別設(shè)為20°，25°，30°，35°，40°；靜止域設(shè)置為長(zhǎng)方體，長(zhǎng)5.85 m、寬5.325 m、高4.15 m，設(shè)置完成后使用布爾運(yùn)算去除固體設(shè)備.絞刀計(jì)算域設(shè)置見(jiàn)圖3.

圖3 絞刀流場(chǎng)計(jì)算域的劃分

鉸刀在旋轉(zhuǎn)過(guò)程中，越靠近鉸刀附近的流場(chǎng)所受影響越大，網(wǎng)格需進(jìn)行局部加密，對(duì)旋轉(zhuǎn)域及旋轉(zhuǎn)域和靜止域交界面上的網(wǎng)格進(jìn)行加密，對(duì)靜止域的網(wǎng)格進(jìn)行適當(dāng)?shù)南∈杌?在mesh中不同下傾角的劃分網(wǎng)格數(shù)分別為：20°網(wǎng)格數(shù)2 571 225，25°網(wǎng)格數(shù)2 573 376,30°網(wǎng)格數(shù)2 573 340，35°網(wǎng)格數(shù)2 574 386，40°網(wǎng)格數(shù)2 579 742.經(jīng)過(guò)網(wǎng)格修改和優(yōu)化，計(jì)算中99%以上的網(wǎng)格質(zhì)量在0.5以上，符合計(jì)算要求.

2.2 邊界條件的設(shè)定

將5組下傾角不同的mesh文件分別導(dǎo)入fluent中，對(duì)模型進(jìn)行邊界條件的設(shè)置.

1) 材料的設(shè)置 計(jì)算中液固兩相分別為水和泥沙顆粒，具體參數(shù)設(shè)置見(jiàn)表2.

2) 邊界的設(shè)置 將入口設(shè)置為速度入口，泥沙和水的流入速度分別取2，5，8 m/s，方向垂直于入口面，出口設(shè)置為自由出口.

3) 壁面的設(shè)置 將絞刀軸與靜止域的接觸面設(shè)置為無(wú)滑移壁面，其他壁面設(shè)置為開(kāi)放式壁面.

4) 旋轉(zhuǎn)域的設(shè)置 有研究顯示，泥沙吸入百分比及顆粒分布均勻度隨絞刀轉(zhuǎn)速呈先升后降的現(xiàn)象，合理轉(zhuǎn)速范圍在30～40 r/min之間[8-9]，計(jì)算中設(shè)置旋轉(zhuǎn)域的轉(zhuǎn)速為35 r/min,方向?yàn)轫槙r(shí)針?lè)较蛐D(zhuǎn).

5) 模型采用mixture模型，湍流方程采用k-ε方程[10-11].

表2 泥沙參數(shù)

3 結(jié)果分析

下傾角α分別為20°，30°，40°時(shí)絞刀刀臂的壓力分布云圖見(jiàn)圖4，由圖4可知，絞刀上最大壓力的位置在刀臂上接近輪轂的地方，此處是刀臂延輪廓線發(fā)展時(shí)外徑最大的位置，也是刀臂與泥沙接觸最充分的地方，比較三種下傾角的壓力云圖發(fā)現(xiàn)，隨著下傾角的增大，最大壓力的范圍向大環(huán)處發(fā)展，因此角度越小最大壓力的擴(kuò)散范圍也就越小.

圖4 不同下傾角刀臂壓力云圖

為了分析不同工況組合下絞刀上的壓力變化情況，在流場(chǎng)分析后處理中從絞刀上取三個(gè)位置進(jìn)行分析，分別是絞刀最大壓力處、刀臂與輪轂連接處、刀臂與大環(huán)連接處.這三處是絞刀的主要承力處及結(jié)構(gòu)支撐點(diǎn)，對(duì)絞刀的結(jié)構(gòu)強(qiáng)度有重要影響，將這三個(gè)點(diǎn)標(biāo)記為①，②，③，三點(diǎn)的具體位置見(jiàn)圖5.

圖5 壓力分析所取三點(diǎn)位置

在絞刀流場(chǎng)分析的后處理中不同下傾角、不同進(jìn)口速度工況下的絞刀上三處位置壓強(qiáng)的變化情況見(jiàn)圖6.

圖6 不同入口速度下絞刀壓力隨下傾角變化趨勢(shì)

絞刀上的壓力變化對(duì)絞刀的使用壽命和結(jié)構(gòu)強(qiáng)度有較大影響，進(jìn)而影響絞刀的切削性能.由圖6可知，當(dāng)下傾角較小時(shí)，絞刀上的壓力較大，隨著下傾角增加，壓力持續(xù)降低，在角度為35°時(shí)壓力達(dá)到最小值，繼續(xù)增大下傾角，壓力反而出現(xiàn)增大的情況.當(dāng)下傾角一定時(shí)，隨著流體的入口速度增加，絞刀上的壓力會(huì)增大，且幅度很明顯.

因此，流體的入口速度對(duì)絞刀壓力的影響比下傾角大，入口速度越小，鉸刀所處流域就越穩(wěn)定，泥沙對(duì)絞刀的沖擊越小，絞刀所承受的載荷也就越小.下傾角過(guò)小時(shí)，絞刀所承受的載荷很大，且此時(shí)的挖泥深度很小，對(duì)實(shí)際工況及地形的要求比較嚴(yán)格；下傾角過(guò)大時(shí)，由于吸泥管布置于橋架內(nèi)，吸管的下傾角也會(huì)很大，此時(shí)泥沙比小角度更易沉降，使泥泵的功耗增加，導(dǎo)致工作效率降低.綜合分析可以得出，在絞吸式挖泥船實(shí)際施工過(guò)程中，需要在較平緩的海況中設(shè)置合理的下傾角來(lái)降低絞刀所承受的載荷，以優(yōu)化其挖掘受力性能.

不同下傾角出口泥沙的平均體積分?jǐn)?shù)變化情況見(jiàn)圖7.由圖7可知，隨著下傾角的增大，出口泥沙的平均體積分?jǐn)?shù)逐漸增大，當(dāng)下傾角增加到35°時(shí)平均體積分?jǐn)?shù)達(dá)到較大值，角度繼續(xù)增大，平均體積分?jǐn)?shù)會(huì)繼續(xù)增加，但是幅度趨于平穩(wěn)，且隨著入口流速的增大，出口泥沙平均體積分?jǐn)?shù)逐漸降低，在入口流速為2 m/s時(shí)出口泥沙平均體積分?jǐn)?shù)最大且隨著下傾角的變化波動(dòng)較大，因此，挖泥船施工時(shí)所處流域越穩(wěn)定，產(chǎn)量就會(huì)越大，挖掘效率也就越高.

圖7 出口泥沙平均體積分?jǐn)?shù)

4 結(jié) 論

1) 在絞刀流場(chǎng)的數(shù)值模擬過(guò)程中，旋轉(zhuǎn)域尺寸小，是流場(chǎng)分析的主要研究對(duì)象，網(wǎng)格需加密.對(duì)于10%以下體積分?jǐn)?shù)的泥沙顆粒，適合使用雙歐拉模型計(jì)算，而對(duì)于較高體積分?jǐn)?shù)，應(yīng)采用mixture模型.

2) 絞刀刀臂為絞刀的連接及支撐結(jié)構(gòu)，其載荷變化對(duì)絞刀的受力性能影響較大，當(dāng)絞刀下傾角在20°～40°范圍內(nèi)時(shí)，刀臂的主要連接和支撐處的壓力及最大壓力隨下傾角的增加有先增后降的趨勢(shì)，且當(dāng)入口流速在2～8 m/s范圍內(nèi)時(shí)，隨著入口流速的增大，壓力升高幅度很大，最大壓力的范圍也隨下傾角的增大向大環(huán)擴(kuò)散；出口泥沙體積分?jǐn)?shù)會(huì)隨著下傾角的增大而增大，且隨著入口流速的增大而降低，因此，合適的下傾角及入口流速較低的工況會(huì)有利改善絞刀的受力性能，并提高工作效率.

3) 影響絞刀壓力及速度分布的影響因素除了下傾角和入口流速，還有絞刀的橫移速度、前進(jìn)步距、泥沙粒徑范圍等因素，本文的數(shù)值模擬簡(jiǎn)化了工況，忽略了這些因素的影響，進(jìn)一步的研究需考慮這些因素的影響.