楊 望，梁 磊，楊 堅(jiān)

(1.廣西大學(xué) 機(jī)械工程學(xué)院，南寧 530004；2.廣西蔗糖產(chǎn)業(yè)協(xié)同創(chuàng)新中心，南寧 530004)

0 引言

我國(guó)甘蔗主產(chǎn)區(qū)受臺(tái)風(fēng)影響，易出現(xiàn)大面積倒伏和折斷，對(duì)甘蔗產(chǎn)量和機(jī)械收獲造成嚴(yán)重影響，嚴(yán)重阻礙甘蔗收獲機(jī)械化的推廣和甘蔗種植業(yè)的發(fā)展[1]。而目前風(fēng)作用甘蔗的動(dòng)力學(xué)仿真研究處于起步階段[2]，相關(guān)研究匱乏，因此進(jìn)行風(fēng)-甘蔗流固耦合動(dòng)力學(xué)仿真研究具有重要意義。本文研究的風(fēng)馬赫數(shù)較低，屬于不可壓縮流體，故采用LS-DYNA不可壓縮流ICFD求解器[3]和強(qiáng)耦合計(jì)算方法，建立風(fēng)-甘蔗流固耦合動(dòng)力學(xué)仿真模型，并驗(yàn)證其構(gòu)建方法的可行性。本研究對(duì)臺(tái)風(fēng)-甘蔗流固耦合系統(tǒng)動(dòng)力學(xué)仿真研究具有重要意義。

1 甘蔗仿真模型

1.1 甘蔗幾何模型

甘蔗主要由莖稈、蔗葉和根系構(gòu)成，本文研究甘蔗與風(fēng)之間的作用，不涉及甘蔗根系，即甘蔗底部采用質(zhì)量塊固定。甘蔗的結(jié)構(gòu)復(fù)雜，且各部分材料參數(shù)差別較大[4]，因此建模時(shí)需要對(duì)甘蔗先進(jìn)行分段，具體如下：

1)將莖稈分為上稈、中稈和下稈。由于上稈直徑變化較大，中下稈直徑變化較小，去除未緊密包裹莖稈的葉鞘后，莖稈頂部往下到出現(xiàn)第1個(gè)蔗節(jié)之間為上稈，剩下的中、下稈等分。

2)根據(jù)蔗葉所處位置，對(duì)蔗葉進(jìn)行分類。莖稈最頂端蔗葉為頂葉，長(zhǎng)在頂葉旁邊的兩片蔗葉和頂葉為上葉，其余位于上稈的蔗葉為中葉，位于莖稈中下稈的蔗葉為下葉。

甘蔗的各段分布如圖1所示。

圖1 甘蔗分段圖Fig.1 The fractal of sugarcane

在保證計(jì)算精度的情況下，為了便于建模、網(wǎng)格劃分及減少計(jì)算時(shí)間，對(duì)甘蔗各部分進(jìn)行如下適當(dāng)簡(jiǎn)化：

1)葉片和葉中脈建模時(shí)簡(jiǎn)化為面，定義單元屬性時(shí)設(shè)定厚度，如圖2所示。葉鞘需與葉片和葉中脈連接起來(lái)，其剖面與圖2(b)類似。

2)頂葉由多張蔗葉包裹一起，在劃分流體網(wǎng)格時(shí)，易出現(xiàn)細(xì)小網(wǎng)格，故在建模時(shí)將其簡(jiǎn)化成一片橫截面為半圓的葉。

3)由圖1可知：上稈的葉鞘緊密包裹莖稈，很難與莖稈分開，而中、下稈的葉鞘相對(duì)與莖稈沒有那么緊密。因此，上稈的蔗葉考慮為通過葉舌與莖稈直接相連，中、下稈的蔗葉考慮為通過葉鞘與莖稈相連，同時(shí)上稈的節(jié)間長(zhǎng)度取葉舌之間的距離，中、下稈節(jié)間取甘蔗節(jié)之間距離。甘蔗葉舌和鞘基直接與莖稈相連，劃分流體網(wǎng)格時(shí)易產(chǎn)生小網(wǎng)格，導(dǎo)致計(jì)算時(shí)間過長(zhǎng)；而葉舌和鞘基主要起連接作用，所以建模時(shí)把葉舌和鞘基向莖稈外水平延伸一些，如圖3和圖4所示。

圖2 蔗葉剖面對(duì)比圖Fig.2 Comparative picture of actual and model leaf section

圖3 甘蔗葉舌對(duì)比圖Fig.3 Comparison between the actual and the model middle leaf

圖4 甘蔗鞘基對(duì)比圖Fig.5 Comparison between the actual and the model bottom leaf

4)由圖1可知：上稈為新生稈，蔗節(jié)和節(jié)間被葉鞘緊密包裹，蔗節(jié)和節(jié)間材料參數(shù)相近，所以建模時(shí)把上稈的蔗節(jié)和節(jié)間作為整體建模；而中、下稈的蔗節(jié)與節(jié)間材料參數(shù)差別大，所以需要建立蔗節(jié)，如圖4所示。

甘蔗平躺后，參考文獻(xiàn)[5]的方法測(cè)量其幾何尺寸，蔗葉從上到下編號(hào)。蔗葉各部分尺寸示意圖及尺寸標(biāo)注如圖5所示。

圖5中建模時(shí)葉鞘底部包裹莖稈，b3根據(jù)莖稈尺寸確定。由于葉中脈從底端到頂端逐漸變細(xì)，所以只測(cè)量葉中脈底端寬f，其厚度取平均2mm。由于葉片厚度變化不大，取平均0.32mm，葉鞘厚度取平均為1.4mm，其余尺寸如表1所示。

圖5 蔗葉示意圖Fig.5 Diagrammatic sketch of sugarcane leaf表1 蔗葉尺寸Table 1 Size of sugarcane leaf

編號(hào)b1/mmb2/cmL1/cmL2/cmL3/cme/(°)f/mm1285.560130-30142405.460130-30163425.860130-30154326.060125-150145445.867.5135-170156405.868130-150167285.471.513034150168354.4511353614014

上稈頂端直徑為20mm，中、下稈直徑取平均值為30mm，建模時(shí)蔗節(jié)與莖稈直徑一致，蔗節(jié)高取平均為20mm，甘蔗莖稈總高為1.64m，上稈長(zhǎng)度為0.74m。莖稈各部分尺寸標(biāo)注示意圖，如圖6所示。

圖6 甘蔗莖稈示意圖Fig.6 Diagrammatic sketch of sugarcane stalk

莖稈的節(jié)間從上到下編號(hào)，節(jié)間長(zhǎng)度h如表2所示。蔗葉與莖稈之間夾角i、葉鞘與蔗葉之間夾角j和葉鞘與莖稈之間夾角k如表3所示。

表2 甘蔗莖稈節(jié)間長(zhǎng)度Table 2 Length of each segment

表3 蔗葉與莖稈的角度Table 3 Angle of leaf to sugarcane stalk

通過CREO建立甘蔗三維模型如圖7所示。其中，甘蔗底部質(zhì)量塊長(zhǎng)×寬×高為0.2m×0.2m×0.2m。

圖7 甘蔗模型Fig.7 Sugarcane model

1.2 甘蔗材料模型

甘蔗各部材料模型均采用彈塑性材料(MAT_PLASTIC_KINEMATIC)[6-7]，參考文獻(xiàn)[4]材料測(cè)試方法，得甘蔗各部材料參數(shù)，如表4所示。

表4 甘蔗各部材料參數(shù)Table 4 Material parameters of each part of sugarcane

甘蔗底部質(zhì)量塊材料模型選用剛體材料(MAT_RIGID)，其密度為7 800kg/m3，彈性模量為2.0×1011Pa，泊松比為0.27。

1.3 甘蔗網(wǎng)格劃分及約束條件

采用HyperMesh軟件對(duì)甘蔗模型進(jìn)行網(wǎng)格劃分。甘蔗莖稈和質(zhì)量塊為實(shí)體單元，采用掃掠方法劃分為六面體網(wǎng)格，網(wǎng)格大小取0.02m，甘蔗莖稈和質(zhì)量塊的單元數(shù)分別為972和204。葉片、葉中脈、葉鞘、葉舌、鞘基為殼單元，厚度在單元屬性中定義，采用自由劃分，為四邊形和三角形單元混合劃分，網(wǎng)格大小取0.02m。葉片、葉中脈、葉鞘、葉舌和鞘基單元數(shù)目分別為1 210、481、114、46、10。劃分完網(wǎng)格后，甘蔗各部共節(jié)點(diǎn)，質(zhì)量塊全約束且甘蔗和質(zhì)量塊采用共節(jié)點(diǎn)相連，通過關(guān)鍵字LOAD_BODY_Z對(duì)甘蔗施加重力。

2 流場(chǎng)仿真模型

2.1 流場(chǎng)幾何模型

流場(chǎng)模型為長(zhǎng)方體流場(chǎng)，如圖8所示。其中，流場(chǎng)風(fēng)源豎直平面與甘蔗最大迎風(fēng)面平行，距離La為0.5m。流場(chǎng)尺寸Lb×Lc×Ld為2.1m×3.4m×3.4m。

圖8 流場(chǎng)模型Fig.8 Fluid model

2.2 流場(chǎng)材料模型

空氣材料模型采用不可壓縮流體材料ICFD_MAT，一個(gè)標(biāo)準(zhǔn)大氣壓下，25℃的干空氣密度為1.184kg/m3，動(dòng)力粘度為1.849×10-5Pa·s。

2.3 流場(chǎng)網(wǎng)格劃分

流場(chǎng)域由風(fēng)源面Sin、出口面Sout、自由邊界面Sfree-slip和障礙物(甘蔗和質(zhì)量塊)表面Snon-slip構(gòu)成，如圖8所示。流場(chǎng)采用流體網(wǎng)格，風(fēng)源面、出口面、自由邊界面劃分為等邊三角形單元，網(wǎng)格大小為0.1m；同時(shí)，為實(shí)現(xiàn)耦合，流場(chǎng)中甘蔗表面和甘蔗結(jié)構(gòu)網(wǎng)格大小相同，封閉流場(chǎng)殼單元數(shù)目為16 883個(gè)。流場(chǎng)體網(wǎng)格通過MESH_VOLUME和MESH_EMBEDSHELL結(jié)合流場(chǎng)封閉的邊界殼單元，自動(dòng)劃分為四面體單元。

2.4 流場(chǎng)邊界條件

利用ICFD_BOUNDARY_PRESCRIBED_VEL在流場(chǎng)的入口定義沿X方向風(fēng)速，即風(fēng)源；流場(chǎng)的出口通過ICFD_BOUNDARY_PRESCRIBED_PRE定義，流場(chǎng)出口壓力為0；自由邊界通過ICFD_BOUNDARY_FREESLIP定義自由滑移面；甘蔗表面通過ICFD_BOUNDARY_NONSLIP定義流固耦合邊界。

3 控制參數(shù)和風(fēng)-甘蔗仿真模型

甘蔗與風(fēng)流固耦合采用強(qiáng)耦合的計(jì)算方法，采用CONTROL_IMPLICIT_GENERAL激活隱式求解實(shí)現(xiàn)強(qiáng)耦合?？紤]流場(chǎng)中湍流的存在，采用ICFD_CONTROL_TURBULENCE選取大渦模擬(LES)湍流模型[8]，通過采用關(guān)鍵字ICFD_CONTROL_FSI控制實(shí)現(xiàn)雙向流固耦合[9]。建立的風(fēng)-甘蔗仿真模型，如圖9所示。

圖9 風(fēng)-甘蔗仿真模型Fig.9 Wind-sugarcane simulation model

4 仿真模型驗(yàn)證

4.1 試驗(yàn)方法及設(shè)備

通過在仿真模型中加載風(fēng)源風(fēng)速后測(cè)定莖稈一位置的加速度曲線及記錄甘蔗受風(fēng)作用過程與在物理試驗(yàn)中加載同一風(fēng)源風(fēng)速獲得的加速度曲線和受風(fēng)作用過程進(jìn)行對(duì)比，驗(yàn)證模型的精度。

在無(wú)自然風(fēng)條件下通過工業(yè)風(fēng)扇(桌子墊高45cm)對(duì)甘蔗施加風(fēng)載荷，扇葉豎直平面與甘蔗最大迎風(fēng)面平行，試驗(yàn)過程中采用動(dòng)態(tài)測(cè)試儀記錄位于莖稈高度83cm處的加速度曲線，同時(shí)采用高速攝像機(jī)記錄甘蔗的受風(fēng)作用過程。圖10為驗(yàn)證試驗(yàn)圖。

圖10 驗(yàn)證試驗(yàn)示意圖Fig.10 Schematic diagram of verification testing

在無(wú)自然風(fēng)條件下，在甘蔗迎風(fēng)面正前方0.5m處放置自制測(cè)風(fēng)架，測(cè)定風(fēng)源的剖面風(fēng)速，從無(wú)風(fēng)的位置向最大風(fēng)速中心每間隔20cm水平移動(dòng)一次，并記錄每次7個(gè)數(shù)字測(cè)風(fēng)儀數(shù)據(jù)；接著繼續(xù)水平移動(dòng)，一直移動(dòng)到另一無(wú)風(fēng)位置結(jié)束。試驗(yàn)示意圖如圖11所示。

圖11 測(cè)風(fēng)試驗(yàn)示意圖Fig.11 Schematic diagram of measuring wind speed

試驗(yàn)設(shè)備：特強(qiáng)型工業(yè)風(fēng)扇(佛山市德勝電器廠FS850，扇葉直徑為850mm，功率為480W)、動(dòng)態(tài)測(cè)試儀(東華測(cè)試有限公司DH5902，采樣頻率50Hz)、加速度傳感器(東華測(cè)試有限公司DH131E)、高速攝像機(jī)(日本Photron 公司FASTCAM Mini UX100)、筆記本電腦2臺(tái)、數(shù)字測(cè)風(fēng)儀7個(gè)(深圳市聚茂源科技有限公司GM8902)、自制鋼制測(cè)風(fēng)架(由底座和豎直桿組成，在豎直桿下方往上每間隔30cm放置1個(gè)測(cè)風(fēng)儀，共放置7個(gè)測(cè)風(fēng)儀，底座高37cm，豎直桿高3m)及白色幕布等。試驗(yàn)地點(diǎn)為廣西大學(xué)農(nóng)機(jī)實(shí)驗(yàn)室，甘蔗品種為園林17號(hào)，試驗(yàn)時(shí)間為2017年7月。

4.2 試驗(yàn)結(jié)果

風(fēng)源面風(fēng)速測(cè)量時(shí)，測(cè)風(fēng)儀編號(hào)按從高到低排列，如表5所示。當(dāng)測(cè)風(fēng)架水平移動(dòng)距離為65cm時(shí)，測(cè)風(fēng)儀位于風(fēng)扇中心正前方，編號(hào)7的測(cè)風(fēng)儀距地面高度為67cm。經(jīng)過對(duì)表5中的風(fēng)速進(jìn)行平均合并處理，獲得風(fēng)源面風(fēng)速。

表5 風(fēng)速測(cè)量結(jié)果Table 5 Wind speed

通過加載風(fēng)源面風(fēng)速，仿真計(jì)算獲得的加速度曲線與物理試驗(yàn)實(shí)測(cè)曲線進(jìn)行對(duì)比，如圖12所示。由圖12可知：仿真曲線與實(shí)測(cè)曲線變化趨勢(shì)一致。圖13為實(shí)際甘蔗和仿真甘蔗受風(fēng)載作用截圖。由圖13可知：仿真和實(shí)際風(fēng)作用甘蔗過程較一致。結(jié)果表明：所建立風(fēng)-甘蔗流固耦合仿真模型構(gòu)建方法可行，可用于模擬風(fēng)作用甘蔗的動(dòng)力學(xué)過程。

圖12 加速度曲線對(duì)比圖Fig.12 Contrastive diagram of acceleration curve

圖13 甘蔗受風(fēng)作用對(duì)比圖Fig.13 Contrastive diagram of sugarcane affected by wind

5 結(jié)論

首先對(duì)實(shí)際甘蔗進(jìn)行了簡(jiǎn)化，然后測(cè)定甘蔗幾何參數(shù)建立了甘蔗幾何模型，測(cè)定了甘蔗材料參數(shù)，建立了甘蔗材料模型。同時(shí)，建立了流場(chǎng)幾何模型，定義材料參數(shù)，對(duì)甘蔗和流場(chǎng)分別網(wǎng)格劃分，并定義邊界條件和控制參數(shù)進(jìn)行計(jì)算。最后，對(duì)甘蔗風(fēng)載仿真模型進(jìn)行了物理驗(yàn)證，結(jié)果表明：風(fēng)-甘蔗流固耦合仿真模型構(gòu)建方法可行，可以用于模擬風(fēng)作用甘蔗的動(dòng)力學(xué)過程，為臺(tái)風(fēng)-甘蔗和自然風(fēng)-作物流固耦合系統(tǒng)仿真提供了依據(jù)。