鄭琳珠，胡京招，梁鑫，謝立全

(1.中交疏浚技術(shù)裝備國(guó)家工程研究中心有限公司，上海 201208;2.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092)

?

絞吸挖泥船絞刀切削過(guò)程數(shù)值模擬及刀臂曲線優(yōu)化

鄭琳珠1，胡京招1，梁鑫2，謝立全2

(1.中交疏浚技術(shù)裝備國(guó)家工程研究中心有限公司，上海 201208;2.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092)

為了研究絞吸挖泥船切削過(guò)程土體與絞刀相互作用的功率消耗和切削產(chǎn)量情況，基于離散元理論，結(jié)合實(shí)船施工數(shù)據(jù)，采用PFC3D程序，將土體離散為有限單元，模擬絞刀切削土體過(guò)程，監(jiān)測(cè)絞刀功率消耗和切削產(chǎn)量，與實(shí)船數(shù)據(jù)對(duì)比分析，驗(yàn)證數(shù)值模擬的可靠性，在此基礎(chǔ)上對(duì)現(xiàn)有絞刀刀臂曲線進(jìn)行改進(jìn)優(yōu)化。

絞吸式挖泥船；絞刀；離散元；PFC3D；數(shù)值模擬

絞刀是絞吸挖泥船的主要挖掘部件，在疏浚施工過(guò)程中，絞刀直接與土體相接觸，通過(guò)橫移與自旋轉(zhuǎn)達(dá)到切削挖掘土體的效果，因此，其挖掘性能是影響絞吸挖泥船產(chǎn)量的主要因素。影響絞刀挖掘性能的因素非常多且復(fù)雜，包括刀臂包角、刀齒切削角、刀齒剪切角和刀齒-土體摩擦角等，以及切削過(guò)程中的橫移速度、旋轉(zhuǎn)速度，這就決定了絞刀頭理論研究的復(fù)雜性[1-2]。目前，國(guó)內(nèi)外對(duì)絞刀的研究主要是基于二維切削理論與巖土力學(xué)的巖土破壞理論，通過(guò)建立刀齒與土體切削破壞模型，用方程來(lái)簡(jiǎn)化分析絞刀切削過(guò)程中的刀齒受力問(wèn)題。這在一定程度上是可行的，但是絞刀作為三維形態(tài)結(jié)構(gòu)，且與土體的切削過(guò)程也是一個(gè)三維過(guò)程，用二維模型很難全面描述整個(gè)切削過(guò)程[3]。同時(shí)，由于各個(gè)刀齒的空間位置不同，以及前后刀齒的相互影響，這就導(dǎo)致各刀齒在施工中切削過(guò)程各不相同，對(duì)產(chǎn)量的挖掘貢獻(xiàn)也各不相同，同一絞刀頭上的不同刀齒磨損以及斷裂程度各不相同。因此，為了分析絞吸挖泥船施工過(guò)程中絞刀的真實(shí)切削情況，對(duì)絞刀與土體的切削過(guò)程進(jìn)行三維描述研究就顯得很有必要了。

散體單元方法屬于一種動(dòng)態(tài)分析方法，將物體假設(shè)為不連續(xù)的離散的單元體，通過(guò)定義單元顆粒與單元顆粒之間的相互作用特征，來(lái)體現(xiàn)物體的物理屬性。考慮到塊體或者顆粒受力后的運(yùn)動(dòng)狀態(tài)以及由此導(dǎo)出受力狀態(tài)隨時(shí)間的變化，以不連續(xù)體(散粒體)力學(xué)的方法研究各個(gè)單元之間的相互接觸和作用，多種類型的單元適應(yīng)了不同問(wèn)題的需要。這種離散單元的理論方法在土體、水流等物體的模擬和研究中得到了較廣泛的使用和認(rèn)可[4-5]。因此，考慮將離散元理論引入到絞刀與土體的切削過(guò)程研究中，在計(jì)算機(jī)的輔助下，采用PFC3D程序?qū)g刀與土體相互作用建立三維數(shù)值模型，模擬并分析絞吸挖泥船絞刀挖掘土體過(guò)程中的功率消耗、產(chǎn)量分布情況。

1 基于離散元理論的三維絞刀建模及土體細(xì)觀參數(shù)確定

1.1 采用落雨法建立海床結(jié)構(gòu)模型

由于砂土等散粒體介質(zhì)在荷載作用下的本構(gòu)關(guān)系十分復(fù)雜，通過(guò)現(xiàn)場(chǎng)或室內(nèi)實(shí)驗(yàn)來(lái)得到反映其真實(shí)力學(xué)規(guī)律的本構(gòu)模型相當(dāng)困難，并且耗時(shí)耗錢。土體的宏觀表現(xiàn)必然與其內(nèi)在的細(xì)、微觀特性相關(guān)。隨著計(jì)算機(jī)功能的不斷增強(qiáng)，從細(xì)觀角度對(duì)土體進(jìn)行模擬研究成為可能。另外由于顆粒流理論是基于顆粒的細(xì)觀力學(xué)特征建立的，其有關(guān)屬性和參數(shù)，包括顆粒大小，可以在整個(gè)模擬過(guò)程中進(jìn)行調(diào)整，而且顆粒流理論克服了傳統(tǒng)連續(xù)介質(zhì)力學(xué)的宏觀連續(xù)性假設(shè)，因此對(duì)于模擬砂土等散粒體介質(zhì)的力學(xué)問(wèn)題具有其特殊的優(yōu)勢(shì)[6-7]。

土體顆粒的制備采用落雨法。所謂落雨法，就是在區(qū)域上方一定高度的位置利用GENERATE命令生成無(wú)規(guī)則排列的小球(balls)，通過(guò)重力作用使小球自由下落，通過(guò)下落高度的調(diào)整，來(lái)控制最終區(qū)域內(nèi)砂子的孔隙率等一系列參數(shù)。因?yàn)槁溆攴梢暂^為真實(shí)地?cái)M合現(xiàn)實(shí)中沙土顆粒的沉降過(guò)程，在數(shù)值模擬中得到廣泛的使用。

1.2 三維絞刀的數(shù)值建模

以通用型絞刀為試驗(yàn)對(duì)象，絞刀配備寬型常規(guī)刀齒，對(duì)其進(jìn)行三維數(shù)值建模。為了方便計(jì)算，將絞刀齒簡(jiǎn)化為一個(gè)矩形平面，取A、B及C點(diǎn)作為矩形的3個(gè)頂點(diǎn)，作刀齒矩形，見圖1。

圖1 絞刀刀齒及齒座安裝示意

根據(jù)設(shè)計(jì)圖紙尺寸，建立通用型絞刀PFC3D模型，見圖2。

圖2 通用型絞刀PFC3D模型

1.3 土體數(shù)值模型細(xì)觀參數(shù)的確定

土體在完成三維建模后需要通過(guò)賦予與實(shí)際土體相一致的微觀參數(shù)來(lái)確定土體的物理性質(zhì)，保證所形成的土體數(shù)值模型與實(shí)際土體具有相同的力學(xué)性質(zhì)。微觀參數(shù)在數(shù)值上并不等于宏觀參數(shù)，但是其力學(xué)性質(zhì)的表征卻可以一致。采用PFC3D三軸試驗(yàn)來(lái)調(diào)試微觀參數(shù)的選取[8-9]。

以細(xì)粉砂土質(zhì)為模擬對(duì)象，利用實(shí)驗(yàn)室土體三軸試驗(yàn)來(lái)調(diào)試數(shù)值模擬土體模型的微觀參數(shù)，達(dá)到模型在宏觀的力學(xué)性質(zhì)上符合室內(nèi)試驗(yàn)的土體模型。通過(guò)不斷調(diào)整土體數(shù)值模型微觀參數(shù)，最終形成與實(shí)驗(yàn)室物理土體三軸試驗(yàn)結(jié)果相一致的土體模型，圖3為模擬試驗(yàn)和室內(nèi)試驗(yàn)的應(yīng)力-應(yīng)變對(duì)比。

圖3 模擬試驗(yàn)和室內(nèi)三軸試驗(yàn)的應(yīng)力-應(yīng)變對(duì)比

由圖3可見，數(shù)值模型土體的模擬試驗(yàn)得出的模型土體應(yīng)力應(yīng)變關(guān)系曲線與實(shí)驗(yàn)室物理土體三軸試驗(yàn)得出的實(shí)際土體應(yīng)力應(yīng)變曲線很接近。此時(shí)所得到的一組細(xì)觀參數(shù)便可以認(rèn)為是與實(shí)際土體宏觀參數(shù)相一致的數(shù)值模擬土體參數(shù)。

最終選定的顆粒參數(shù)可見表1。表1顯示的僅僅是微觀參數(shù)，并不直接等于宏觀參數(shù)。

表1 顆粒的微觀參數(shù)

2 絞刀切削過(guò)程數(shù)值模擬及分析

將生成的“通用型絞刀”模型置入土體當(dāng)中，對(duì)土體進(jìn)行分層染色，結(jié)果見圖4。

由圖4可見，在船行方向(絞刀推進(jìn)方向)存在一個(gè)斜坡，這是為了模擬真實(shí)生產(chǎn)穩(wěn)定切削過(guò)程中的絞刀推進(jìn)狀態(tài)。僅對(duì)其中的一個(gè)絞刀擺移過(guò)程進(jìn)行模擬，絞刀在向左平移過(guò)程中的豎直坐標(biāo)(z坐標(biāo))不改變。為監(jiān)測(cè)數(shù)據(jù)方便，絞刀進(jìn)入土體之前為完全懸空，絞刀從右到左移動(dòng)，逐漸進(jìn)入土體。

圖4 “絞刀-土體”相互作用過(guò)程模擬

“絞刀-土體”相互作用過(guò)程參數(shù)設(shè)置為：絞刀模型前進(jìn)速度為0.26 m/s，轉(zhuǎn)速為25 r/min，與工程實(shí)際的前進(jìn)和轉(zhuǎn)動(dòng)速度保持一致。切削時(shí)長(zhǎng)為23 s，切削長(zhǎng)度為6 m。在計(jì)算機(jī)輔助下運(yùn)行程序，絞刀模擬施工過(guò)程切削土體。在模擬過(guò)程中，疊加絞刀頭上48片刀齒的功率之和作為絞刀頭的消耗功率[10]。

其中，單個(gè)絞刀刀齒的功率計(jì)算如下。

(1)

式中:Pi為第i個(gè)絞刀刀齒的功率，kW;F為絞刀刀齒所受的切向力，kN;n為絞刀的自轉(zhuǎn)速度，r/min;di為絞刀刀齒距離絞刀中軸線的距離，m。

因此，絞刀的消耗功率計(jì)算如下。

(2)

施工過(guò)程中，吸泥管在泥泵的作用下將絞刀切削下的土體混合海水吸走。為了模擬泵吸效應(yīng)，在絞刀切削過(guò)程中將進(jìn)入絞刀刀臂外輪廓范圍內(nèi)的土體顆粒移除，統(tǒng)計(jì)移除土體顆粒的體積，以此作為計(jì)算絞刀工作時(shí)挖方量的基礎(chǔ)，記錄其體積v，作為絞刀的切削量[11-12]?？偼诜搅繛?/p>

(3)

式中：V為總挖方量；vi為每個(gè)土體顆粒的體積；P為土體的孔隙率。

在模擬過(guò)程中，監(jiān)測(cè)“通用型絞刀”的消耗功率和產(chǎn)量，對(duì)比數(shù)值模擬結(jié)果與同樣工況下實(shí)際施工數(shù)據(jù)，見圖5。

圖5 “通用型”絞刀功率消耗數(shù)值模擬與實(shí)船數(shù)據(jù)對(duì)比

由圖5可見，“通用型”絞刀數(shù)值模擬中功率消耗隨時(shí)間變化與實(shí)船施工中功率消耗隨時(shí)間變化基本吻合。絞刀由開始的空轉(zhuǎn)到逐漸切削到土體；因此絞刀消耗功率從0逐漸變大。絞刀在轉(zhuǎn)動(dòng)切削土體過(guò)程中，刀臂與刀臂之間存在切削間隔，且土體模型顆粒為70 mm，為均勻分布土體，因此，絞刀功率消耗曲線出現(xiàn)有規(guī)律的震蕩，且振幅較小。針對(duì)實(shí)船施工過(guò)程，在開始階段(0～1 s范圍)，絞刀從水體進(jìn)入土體，參與切削的刀臂以及刀齒最多，此時(shí)絞刀功率最大，隨著絞刀的橫移，絞刀切削逐漸穩(wěn)定，因此，最大功率從開始的630 kW逐漸降低為350 kW。在實(shí)際工程中，刀臂與刀臂之間存在間隙，刀臂切削土體為非連續(xù)性切削，且海底土體不均勻，并含有大尺寸礫石；因此，絞刀功率消耗曲線出現(xiàn)周期性震蕩的同時(shí)會(huì)在某些時(shí)刻出現(xiàn)較大的功率突變。

從功率消耗對(duì)比來(lái)看，計(jì)算機(jī)數(shù)值模擬的結(jié)果與實(shí)船數(shù)據(jù)較為吻合，平均功率偏差8 kW，偏差值占實(shí)際功率的2.9%，屬于可以接受的偏差范圍，驗(yàn)證了計(jì)算機(jī)數(shù)值模擬在功率監(jiān)測(cè)和計(jì)算方面的有效性。圖6為數(shù)值模擬過(guò)程中，“通用型”絞刀切削土體的體積隨時(shí)間的變化。

圖6 “通用型”絞刀數(shù)值模擬切削量隨時(shí)間的變化

由圖6可見，從0～2 s過(guò)程中，切削量為0，從2 s以后，絞刀切削量逐漸增加，穩(wěn)定后大致成線性增長(zhǎng)。由圖6可以計(jì)算出絞刀單位時(shí)間切削量為4 320 m3/h。在實(shí)際施工中，同樣工況下“通用型”絞刀的單位時(shí)間挖掘量平均為3 500 m3/時(shí)。數(shù)值模擬計(jì)算結(jié)果為4 320 m3/時(shí)，比實(shí)際挖方量偏差820 m3，偏差率約為23%，這是由于海上施工時(shí)實(shí)際工況較為惡劣，海底土體顆粒不均勻，含有礫石顆粒；而在數(shù)值模擬過(guò)程中，土體顆粒直徑均為70 mm小顆粒，且均勻分布，絞刀切削過(guò)程比較平穩(wěn)。因此，絞刀在數(shù)值模擬過(guò)程中計(jì)算出的產(chǎn)量偏大，偏差率達(dá)到23%。

總體來(lái)看，數(shù)值模擬的過(guò)程與實(shí)際施工過(guò)程較為吻合，偏差在可接受的范圍內(nèi)。進(jìn)而驗(yàn)證了基于離散元方法的絞刀計(jì)算機(jī)數(shù)值模擬程序?qū)Ξa(chǎn)量的計(jì)算和監(jiān)測(cè)是有效的。

3 絞刀改進(jìn)優(yōu)化

絞刀在切削土壤過(guò)程中，刀臂間隙直接影響著切削連續(xù)性，進(jìn)而影響絞刀的切削產(chǎn)量，但是，刀臂間隙過(guò)小亦會(huì)引起刀臂堵塞。因此，針對(duì)不同土質(zhì)，最優(yōu)的刀臂間隙和刀臂輪廓曲線是不同的。基于此，在原有一款典型冠形絞刀基礎(chǔ)上優(yōu)化刀臂輪廓形狀和刀臂間隙，使其刀臂輪廓形狀更傾向于錐形，更適用于細(xì)沙土質(zhì)的挖掘，形成優(yōu)化型錐形絞刀。

優(yōu)化前后絞刀刀臂的三維立體圖形如圖7、圖8。

圖8 優(yōu)化后錐形絞刀刀臂模型

根據(jù)優(yōu)化前后不同絞刀的設(shè)計(jì)數(shù)據(jù)，利用PFC軟件分別建立絞刀PFC3D模型。

以實(shí)船施工工況數(shù)據(jù)為基礎(chǔ)，在數(shù)值模擬中，設(shè)定絞刀的自轉(zhuǎn)速度為26 r/min，前進(jìn)速度為0.26 m/s，傾斜角度為27°，與實(shí)際工況保持一致。

以相同參數(shù)設(shè)置的土體模型為基礎(chǔ)，將優(yōu)化前冠形絞刀和優(yōu)化后絞刀分別與土體作用，在計(jì)算機(jī)輔助下利用PFC3D程序進(jìn)行數(shù)值模擬，切削時(shí)間均設(shè)置為20 s。監(jiān)測(cè)在整個(gè)模擬切削過(guò)程中的功率消耗情況和切削量變化情況。

將監(jiān)測(cè)所得數(shù)據(jù)繪成圖9、10。

圖9 優(yōu)化前后絞刀切削過(guò)程消耗功率對(duì)比

從功率消耗來(lái)看，冠形絞刀在與土體切削作用過(guò)程中所消耗的功率大于錐形絞刀所消耗的功率值。

圖10 優(yōu)化前后絞刀切削過(guò)程土體切削量對(duì)比

對(duì)比發(fā)現(xiàn)，優(yōu)化前后絞刀在前期切削過(guò)程中產(chǎn)量相差不大，曲線非常接近，在達(dá)到穩(wěn)定狀態(tài)之后，優(yōu)化前后絞刀曲線開始出現(xiàn)分叉，產(chǎn)量相差逐漸增大。單位時(shí)間產(chǎn)量，錐形絞刀比冠形絞刀切削產(chǎn)量高大約10%。

4 結(jié)論

通過(guò)將實(shí)船施工數(shù)據(jù)與采用PFC3D的三維動(dòng)態(tài)數(shù)值模擬方法所得模擬結(jié)果進(jìn)行對(duì)比，驗(yàn)證了PFC3D三維動(dòng)態(tài)數(shù)值模擬方法的可靠性，說(shuō)明將基于離散元理論的數(shù)值模擬方法引入到絞吸挖泥船絞刀切削土體過(guò)程的研究中是可行的。這與之前采用ANSYS有限元方法對(duì)絞刀孤立的進(jìn)行靜態(tài)切削受力研究有著較大的不同，基于離散元理論的PFC3D動(dòng)態(tài)數(shù)值模擬將絞刀與土體相結(jié)合進(jìn)行動(dòng)態(tài)研究，實(shí)時(shí)監(jiān)測(cè)土體對(duì)刀齒產(chǎn)生的切削反作用力變化，更接近實(shí)際施工工況。將離散元理論的數(shù)值模擬方法引入到絞刀的開發(fā)研究中，有效解決了絞刀優(yōu)化研發(fā)過(guò)程中的消耗功率監(jiān)測(cè)、產(chǎn)量預(yù)報(bào)問(wèn)題。同時(shí)，利用計(jì)算機(jī)輔助模擬，避免了實(shí)驗(yàn)室物模實(shí)驗(yàn)成本高、時(shí)間長(zhǎng)的缺陷，具有較好的應(yīng)用價(jià)值。

[1] 陳曉華，楊啟.大型絞吸挖泥船挖巖絞刀切削力計(jì)算分析[J].武漢理工大學(xué)學(xué)報(bào)，2013,37(1):106-110.

[2] 歐陽(yáng)義平.巖石疏浚用刀齒的切削機(jī)理研究[D].上海：上海交通大學(xué)，2013.

[3] 朱文亮,倪福生,張德新.挖泥船絞刀的三維建模方法[J].船海工程，2007,36(1):45-48.

[4] 楊楨毅.絞吸式挖泥船絞刀結(jié)構(gòu)與性能優(yōu)化的研究[D].武漢：武漢理工大學(xué)，2010.

[5] 王泳嘉,刑紀(jì)波.離散元法及其在巖土力學(xué)中的應(yīng)用[M].沈陽(yáng)：東北大學(xué)出版社，1991.

[6] 焦玉勇.三維離散單元法及其應(yīng)用[D].武漢：中國(guó)科學(xué)院武漢巖土力學(xué)研究所，1998.

[7] 徐泳,孫其誠(chéng),張凌,等.顆粒離散元法研究進(jìn)展[J].力學(xué)進(jìn)展，2003,33(2)：251-260.

[8] 周健,池毓蔚,池永,等.砂土雙軸試驗(yàn)的顆粒流模擬[J].巖土工程學(xué)報(bào)，2000,22(6):701-704.

[9] 周健,池永,池毓蔚,等.顆粒流方法及PFC2D程序[J].巖土力學(xué),2000,21(3):271-274.

[10] KWAPINSKA M, SAAGE G, TSOTSAS E. Mixing of particles in rotary drums: A comparison of discrete element simulations with experimental results and penetration models for thermal processes[J]. Powder technology,2006(161):69-78.

[11] ALUKO O B, CHANDLER H W. Characterization and modeling of brittle fracture in two-dimensional soil cutting[J]. Biosystems engineering,2004，88(3):369-381.

[12] HE J. VLASBLOM W. Modeling of saturated sand cutting with large rake angle[C]. 15thworld dredging congress, June 1998, Las Vegas, Nevada,USA,1998:905-920.

Numerical Simulation of Cutting Process and Optimization of Cutter Arm Curve for the Cutter Suction Dredger

ZHENG Lin-zhu1, HU Jing-zhao1, LIANG Xin2, XIE Li-quan2

(1.CCCC National Engineering Research Center of Dredging Technology and Equipment Co.,Ltd., Shanghai 201208, China;2.College of Civil Engineering, Tongji University, Shanghai 200092 China)

In order to research the power dissipation and cutting production in cutting process, based on distinct element method, the cutting process of cutter head was numerically simulated in PFC3D program to analyze the cutting production and the power dissipation of cutter head. Comparing the numerical results with that of the working dredger showed that the numerical simulation is valid. Based on the numerical results, the cutter arm curve of the cutter suction dredger was optimized.

cutter suction dredger; cutter head; distinct element method; PFC3D program; numerical simulation

10.3963/j.issn.1671-7953.2017.04.016

2016-09-27

鄭琳珠(1965—)，女，學(xué)士，高級(jí)工程師

研究方向：疏浚船舶特種設(shè)備研發(fā)

U674.31

A

1671-7953(2017)04-0072-05

修回日期：2016-10-23