王永濤，劉堅，李家春，李繼學(xué)，蔡家斌

(1. 湖南大學(xué)，汽車車身先進(jìn)設(shè)計制造國家重點(diǎn)實驗室，長沙市，410082； 2. 貴州省水利科學(xué)研究院， 貴陽市，550002； 3. 貴州大學(xué)機(jī)械工程學(xué)院，貴陽市，550025)

0 引言

水肥一體化技術(shù)以水肥一體化施肥機(jī)為核心，采用現(xiàn)代控制方法通過灌溉施肥管網(wǎng)，實現(xiàn)精準(zhǔn)灌溉、精量施肥。該技術(shù)具有灌溉與施肥同步、節(jié)水、降施和減污等優(yōu)點(diǎn)，對于我國發(fā)展現(xiàn)代農(nóng)業(yè)具有重要意義。國外如荷蘭、以色列等國家對于水肥一體化施肥機(jī)的研究較早，且水平較為領(lǐng)先。尤以FERTIKIT 3G水肥一體化為代表，其常用的PL型工作壓力范圍250 kPa～650 kPa，文丘里射流器是水肥一體化施肥器的關(guān)鍵部件，其壓力差由吸肥泵提供。我國在施肥機(jī)領(lǐng)域原創(chuàng)性技術(shù)偏少，國際競爭力不強(qiáng)，多以整機(jī)進(jìn)口為主[1]。同時，因各區(qū)域水文地質(zhì)條件、灌溉習(xí)慣、種植結(jié)構(gòu)等差異，亟需對該類施肥機(jī)進(jìn)行結(jié)構(gòu)優(yōu)化設(shè)計，提高適應(yīng)性。

在壓差方面：鄧蘭生等[2]采用壓差施肥罐進(jìn)行施肥條件下壓差、流量以及肥料的品種、形態(tài)、用量等因素對施肥時間的影響試驗研究，得出了對于液體肥料流量是影響施肥時間的直接因素。范軍亮等[3]通過田間試驗綜合評價了施肥罐兩端壓差(50 kPa，100 kPa，150 kPa，200 kPa和250 kPa)和管道布置方式(縱向一端，縱向中間，橫向一端和橫向中間供水)對系統(tǒng)灌水與施肥均勻性的影響，研究表明施肥罐壓差是決定施肥均勻性的直接因素，應(yīng)減低壓差，且管道布置應(yīng)采用橫向布置。孟一斌等[4]認(rèn)為施肥罐壓差是影響肥液濃度變化的最直接因素，影響微灌系統(tǒng)施肥均勻性。施肥罐的流量與隨壓差的增大呈冪函數(shù)關(guān)系，隨時間增加施肥罐的出口肥液濃度降低。韓啟彪等[5]以9.42 L圓柱體模型為例，使用CFD技術(shù)對進(jìn)口流量1.425 m3/h和0.71 m3/h時的施肥罐質(zhì)量濃度衰減規(guī)律進(jìn)行了模擬研究，并與前人理論比較，初步探討了CFD技術(shù)在壓差施肥罐質(zhì)量濃度衰減研究中的可行性。結(jié)果顯示：質(zhì)量濃度衰減過程與封俊等理論公式的計算結(jié)果基本一致，模擬得出的施肥結(jié)束時間與阿莫斯·泰奇經(jīng)驗公式偏差僅3%左右，壓差施肥罐濃度衰減研究中使用CFD模擬是可行的。在施肥器結(jié)構(gòu)方面：嚴(yán)海軍等[6]基于雷諾時均Navier-Stokes方程和標(biāo)準(zhǔn)k-ε湍流模型，利用Fluent軟件對文丘里施肥器的內(nèi)部流動進(jìn)行了數(shù)值模擬，并對數(shù)值計算方法的正確性進(jìn)行了驗證。根據(jù)效率最高原則，得出喉管直徑比應(yīng)為1.2～1.3。邱振宇等在試驗基礎(chǔ)上對并聯(lián)式文丘里施肥器進(jìn)行模擬仿真研究，研究表明收縮角α最佳取值為23°～29°，擴(kuò)大角β最佳取值為7°～11°；劉永華等[7]通過三因素三水平的正交試驗方案，得出基于CFD數(shù)值計算的文丘里吸肥器最優(yōu)結(jié)構(gòu)參數(shù)組合,漸縮角α為20°，漸擴(kuò)角β為8°，喉部直徑d0為6 mm，且當(dāng)吸肥管與文丘里主管道呈40°傾角時吸肥性能最優(yōu)。模擬數(shù)據(jù)與水肥一體化灌溉施肥機(jī)運(yùn)行數(shù)據(jù)表明，CFD數(shù)值優(yōu)化后的文丘里吸肥器吸肥流量提高約38.6%，與模糊自動控制系統(tǒng)相配合的水肥一體化灌溉施肥機(jī)的吸肥流量總體提高約47.6%，節(jié)能效果顯著。在T型出口結(jié)構(gòu)方面：李家春、田莉等通過一通道、二通道、三通道文丘里射流器的仿真對比分析，提出了最優(yōu)結(jié)構(gòu)為三通道[8-9]。李繼學(xué)等通過改變主管道長度和射流泵間距對吸肥器結(jié)構(gòu)進(jìn)行優(yōu)化并分析每條射流泵吸肥量，確定吸肥器結(jié)構(gòu)為主管道所有三通向左側(cè)同時縮短15 mm 作為吸肥器新結(jié)構(gòu)，吸肥量趨于相同且吸肥量高于優(yōu)化前。在控制模型及方法方面：李加念等[10]通過調(diào)節(jié)文丘里施肥機(jī)吸入管段的電磁閥PWM信號，實現(xiàn)了射流器的變量施肥，綜合試驗結(jié)果提出了電磁閥最佳頻率為6 Hz，施肥裝置最佳入口壓力范圍為150 kPa～250 kPa；袁洪波等[11]設(shè)計了一種營養(yǎng)液制備的裝備及數(shù)學(xué)模型，采用基于增量式PID算法和改進(jìn)Smith預(yù)估器的營養(yǎng)液調(diào)控算法，加快了系統(tǒng)的響應(yīng)速度，提高了控制精度。以上研究均具有較強(qiáng)的借鑒意義，但對并聯(lián)四通道文丘里射流器吸肥結(jié)構(gòu)的優(yōu)化、工程實際中四通道水肥一體化施肥機(jī)適宜工作壓力和安裝位置等方面研究較為缺乏。

本文針對四通道文丘里射流器Z型布置條件下，各通道不同安裝位置導(dǎo)致吸肥能力和振動差異較大等問題，采用理論分析與仿真相結(jié)合的方法，得出了并聯(lián)四通道文丘里射流器的優(yōu)化結(jié)構(gòu)，實現(xiàn)了文丘里射流器吸肥特性的改善。為進(jìn)一步探討山地地勢起伏、水肥一體化管網(wǎng)壓差變化大，施肥機(jī)適宜的安裝位置問題，采用四通道文丘里射流器不同入口壓力下吸肥量的仿真，得出在本文優(yōu)化結(jié)構(gòu)下施肥機(jī)適宜的工作入口壓力，為并聯(lián)四通道文丘里施肥機(jī)的安裝位置設(shè)計、施工等工作提供參考。

1 工作原理

文丘里射流器吸肥原理是具有一定壓力的水流通過錐形的噴嘴時，水流速度發(fā)生快速增大的同時在腔內(nèi)產(chǎn)生負(fù)壓，肥料溶液在外界大氣壓的作用下通過吸肥管道被吸入到文丘里射流器內(nèi)。出水段通過吸肥泵將水肥混合液有壓輸至田地，實現(xiàn)水肥一體化，圖1為文丘里施肥器結(jié)構(gòu)示意圖。

圖1 文丘里施肥器結(jié)構(gòu)示意圖

為便于研究，假設(shè)文丘里射流器中的流體為不可壓縮流體，遵循質(zhì)量和能量守恒定律，則其工作原理滿足伯努利方程和連續(xù)性方程。

伯努利方程

(1)

式中：P1——入口壓力，Pa；

P0——喉部壓力，Pa；

V1——入口斷面液體的平均速度，m/s；

V0——喉部斷面液體的平均速度，m/s；

Z1——入口斷面上任一點(diǎn)相對于基準(zhǔn)面的高程，m；

Z0——喉部斷面上任一點(diǎn)相對于基準(zhǔn)面的高程，m；

g——重力加速度，m/s2；

ρ——液體的密度，kg/m3；

hw1-0——入口斷面至喉部斷面管路的各種損失，m。

連續(xù)性方程

(2)

式中：Q1——入口流體流量，m3/s；

Q0——喉部流體流量，m3/s；

d1——入口斷面直徑，m；

d0——喉部斷面直徑，m。

質(zhì)量守恒

Q0+Q=Q2

(3)

式中：Q——吸肥管肥料母液流量，m3/s。

若不計吸入管路的各種損失，且文丘里射流器為水平放置，通過上述方程可以推導(dǎo)出流速及流量計算公式，如式(4)所示。

(4)

式中：v——吸肥管液體的平均速度，m/s；

h——以肥液面為起點(diǎn)，文丘里射流器喉部中心距肥液面的垂直高度，m，文丘里射流器喉部中心高于肥液面為負(fù)，低于肥液面為正；

d——吸肥管直徑，m。

2 計算模型

基于FERTIKIT 3G旁路吸肥式水肥一體化施肥機(jī)，通過增加一條吸肥通道構(gòu)建并聯(lián)四通道旁路吸肥式水肥一體化施肥機(jī)，并對其并聯(lián)四通道文丘里射流器吸肥結(jié)構(gòu)優(yōu)化及仿真。該并聯(lián)四通道文丘里射流器結(jié)構(gòu)可配置為3路施肥通道，1路施酸通道，實現(xiàn)對不同類型單元素液體肥料的吸取，較好滿足農(nóng)業(yè)生產(chǎn)實際需要。采用Solidworks完成主管進(jìn)水管、文丘里射流器、營養(yǎng)液出水管的三維建模，按照空間相對位置關(guān)系進(jìn)行裝配完成四通道文丘里射流器旁路吸肥結(jié)構(gòu)三維建模，如圖2所示。根據(jù)實際應(yīng)用條件，模擬分析中水源管道入口靜壓設(shè)定為250 kPa，文丘里射流器出口與吸肥泵連接，體積流量設(shè)定為0.002 2 m3/s，吸肥泵出口靜壓為450 kPa，文丘里射流器吸肥口及營養(yǎng)液出口環(huán)境壓力設(shè)定為101.325 kPa。

圖2 并聯(lián)四通道文丘里射流器旁路吸肥結(jié)構(gòu)物理模型

為便于分析，將圖2物理模型進(jìn)一步簡化為四通道文丘里射流器的T型結(jié)構(gòu)布置。依據(jù)微尺度理論，在結(jié)構(gòu)彎曲多變的流道內(nèi)，除部分靠近壁面位置外，其他流體基本為湍流，故采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行數(shù)值模擬。為提高網(wǎng)格劃分效率，采用非結(jié)構(gòu)化網(wǎng)格對模型劃分網(wǎng)格，在變化大的復(fù)雜位置(吸肥口、拐角和斜面)，進(jìn)行網(wǎng)格加密處理[12-13]。四通道文丘里射流器網(wǎng)格劃分結(jié)果為：總網(wǎng)格數(shù)13 949個，流體網(wǎng)格13 949 個，接觸固體的流體網(wǎng)格9 605個，圖3為四通道文丘里射流器網(wǎng)格劃分圖。

圖3 四通道文丘里射流器網(wǎng)格劃分圖

3 仿真及結(jié)構(gòu)優(yōu)化

3.1 仿真分析

采用有限體積法對四通道文丘里射流器進(jìn)行仿真分析。當(dāng)滿足最大行程為4.0或分析間隔時間為0.5 s時結(jié)束仿真。將三維模型導(dǎo)入道FloEFD軟件中，完成網(wǎng)格劃分后，軟件運(yùn)行迭代105次后結(jié)束仿真，輸出仿真分析結(jié)果。其中，圖4為四通道文丘里射流器的靜壓云圖，圖5為四通道文丘里射流器的速度云圖。通過靜壓云圖和速度云圖可以看出流體在喉部處速度變大，壓力變小，該仿真分析結(jié)果進(jìn)一步驗證了伯努利方程。

圖4 四通道文丘里射流器的靜壓云圖

四通道文丘里射流器的吸肥量仿真數(shù)據(jù)見圖6，由圖6可知通道1、通道2、通道3和通道4的吸肥量不均勻。通道1和通道4吸肥量接近且吸肥量維持在700～750 L/h；通道2吸肥量維持在550 L/h左右，通道3吸肥量維持在462.5 L/h左右，各通道吸肥量差值最大為287.5 L/h。表明此時射流器出現(xiàn)不均勻的振動導(dǎo)致各通道吸肥量不平衡，主要原因為流體經(jīng)過通道1、通道2、通道3和通道4時產(chǎn)生的能量損失不同。

圖5 四通道文丘里射流器的速度云圖

圖6 四通道文丘里射流器吸肥量

進(jìn)一步對段管內(nèi)1-2斷面的能量損失hw1-2、2-3 斷面的能量損失hw2-3，3-4斷面的能量損失hw3-4，4-5斷面的能量損失hw4-5展開分析。四通道文丘里射流器采用Z型布置，其結(jié)構(gòu)可簡化為圖7。

圖7 四通道文丘里射流器Z型布置

射流器能量損失見式(5)，由式(5)可以看出射流器能量損失與吸肥器斷面管徑、距離和流速呈非線性關(guān)系。射流器的能量損失可分為1-2斷面的沿程水頭損失hf1-2、2-3斷面的沿程水頭損失hf2-3、3-4斷面的沿程水頭hf3-4、4-5斷面的hf4-5沿程水頭損失和1-2斷面的局部水頭損失hm1-2，2-3斷面的局部水頭損失hm2-3，3-4斷面的局部水頭損失hm3-4，4-5 斷面的局部水頭損失hm4-5局部水頭損失[14-16]。

(5)

式中：λ——沿程阻力系數(shù)；

ε——局部阻力系數(shù)；

l1-2——文丘里射流器1-2斷面距離；

l2-3——文丘里射流器2-3斷面距離；

l3-4——文丘里射流器3-4斷面距離l3-4；

l4-5——文丘里射流器4-5斷面距離；

v1——文丘里射流器1-2斷面液體的平均速度，m/s；

v2——文丘里射流器2-3斷面液體的平均速度，m/s；

v3——文丘里射流器3-4斷面液體的平均速度，m/s；

v4——文丘里射流器4-5斷面液體的平均速度，m/s。

進(jìn)一步，文丘里射流器流速與壓力滿足式(6)。

(6)

式中：d2——出口斷面直徑，m。

由式(5)、式(6)可以看出，對于具體某一文丘里射流器來說，管徑d1、d2保持不變，通過改變文丘里射流器斷面距離l1-2、l2-3、l3-4、l4-5和壓力P1，進(jìn)而改變射流器的能量損失實現(xiàn)動態(tài)調(diào)節(jié)其吸肥量趨于一致和最大的目的[17-18]。該分析為進(jìn)一步改善文丘里射流器吸肥特性，實現(xiàn)流量與能力趨于相同和振動平衡提供了理論依據(jù)。

3.2 結(jié)構(gòu)優(yōu)化

本文構(gòu)建四通道文丘里射流器如圖8所示，它是由4個相同的T型通道組成，其部分參數(shù)為：La=105 mm，內(nèi)徑為50 mm，4通T型出口兩側(cè)長度L1=L2=L3=L4=L5=L6=L7=L8=36 mm。由四通道文丘里射流器仿真分析可知，通過對四通道文丘里射流器L1，L2，L3，L4，L5，L6，L7，L8的長度，實現(xiàn)射流器的能量損失調(diào)節(jié)、振動均衡和吸肥量均勻目標(biāo)。

圖8 四通道文丘里射流器尺寸圖

3.3 結(jié)構(gòu)優(yōu)化

同上所述，仿真分析中設(shè)定水源管道入口靜壓為250 kPa，文丘里射流器出口與吸肥泵連接，體積流量設(shè)定為0.002 2 m3/s，吸肥泵出口靜壓為450 kPa，文丘里射流器吸肥口及營養(yǎng)液出口環(huán)境壓力設(shè)定為101.325 kPa。基于四通道文丘里射流器對稱性的要求，設(shè)定T型三通縮減梯度分別為ΔL1=5 mm，ΔL2=10 mm，ΔL3=15 mm，ΔL4=20 mm。縮進(jìn)方向分別為L1、L3、L5、L7同時向左縮進(jìn)，L2、L4、L6、L8同時向右縮進(jìn)，L2、L3、L4、L5、L6、L7同時向中間縮進(jìn)三種情況。利用FloEFD軟件分析不同結(jié)構(gòu)下通道1、通道2、通道3、通道4射流器的吸肥量。

3.3.1 同時向左縮進(jìn)L1、L3、L5、L7

四通道文丘里射流器在總體機(jī)架結(jié)構(gòu)尺寸固定的情況下，通過向左減小通道1、通道2、通道3、通道4射流器的安裝位置，實現(xiàn)射流器整體安裝位置的左移，縮減梯度分別為向左縮進(jìn)ΔL1=5 mm，ΔL2=10 mm，ΔL3=15 mm，ΔL4=20 mm。利用FloEFD軟件分析四通道文丘里射流器不同位置的吸肥量見圖9～圖12。

圖9 左縮進(jìn)5 mm吸肥量

圖10 左縮進(jìn)10 mm吸肥量

圖11 左縮進(jìn)15 mm吸肥量

圖12 左縮進(jìn)20 mm吸肥量

3.3.2 同時向右縮進(jìn)L2、L4、L6、L8

四通道文丘里射流器在總體機(jī)架結(jié)構(gòu)尺寸固定的情況下，通過向右減小通道1、通道2、通道3、通道4射流器的安裝位置，實現(xiàn)射流器整體安裝位置的右移，縮減梯度分別為向右縮進(jìn)ΔL1=5 mm，ΔL2=10 mm，ΔL3=15 mm，ΔL4=20 mm。利用FloEFD軟件分析四通道文丘里射流器不同位置的吸肥量見圖13～圖16。

圖13 右縮進(jìn)5 mm吸肥量

圖14 右縮進(jìn)10 mm吸肥量

圖15 右縮進(jìn)15 mm吸肥量

圖16 右縮進(jìn)20 mm吸肥量

3.3.3 同時向中間縮進(jìn)L2、L3、L4、L5、L6、L7

四通道文丘里射流器在總體機(jī)架結(jié)構(gòu)尺寸固定的情況下，通過向中間減小通道1、通道2、通道3、通道4射流器的安裝位置，實現(xiàn)射流器整體安裝位置的中間移動，縮減梯度分別為向中間縮進(jìn)ΔL1=5 mm，ΔL2=10 mm，ΔL3=15 mm，ΔL4=20 mm。利用FloEFD軟件分析四通道文丘里射流器不同位置的吸肥量見圖17～圖20。

圖17 中間縮進(jìn)5 mm吸肥量

圖18 中間縮進(jìn)10 mm吸肥量

圖19 中間縮進(jìn)15 mm吸肥量

圖20 中間縮進(jìn)20 mm吸肥量

3.3.4 結(jié)果分析

向左縮進(jìn)、右縮進(jìn)和中間縮進(jìn)三個縮進(jìn)方向，5 mm、10 mm、15 mm、20 mm四個縮進(jìn)梯度共12個水平，通道1、通道2、通道3、通道4的吸肥量仿真分析結(jié)果見表1。為提高數(shù)據(jù)的準(zhǔn)確性，重復(fù)仿真10次并取平均值，采用方差表征通道1、通道2、通道3、通道4吸肥量的均勻程度。由表1可知向右縮進(jìn)20 mm的試驗組RI4，吸肥量為2 317.20 L/h，基本持平或高于其他試驗組吸肥量水平。且RI4的方差為218.43，低于其他試驗組。仿真結(jié)果表明RI4試驗組吸肥量大且均勻性最好。故采用試驗組RI4向右縮進(jìn)20 mm為優(yōu)化后的新結(jié)構(gòu)OS，如圖21所示。

表1 不同縮進(jìn)條件下吸肥量Tab. 1 Fertilizer absorption under different indentation conditions

圖22～圖25為向右縮進(jìn)20 mm新結(jié)構(gòu)OS的流體跡線圖和云圖。

圖21 新結(jié)構(gòu)OS尺寸圖

圖22 新結(jié)構(gòu)OS的速度跡線圖

圖23 新結(jié)構(gòu)OS的靜壓跡線圖

圖24 新結(jié)構(gòu)OS的速度云圖

圖25 新結(jié)構(gòu)OS的靜壓云圖

4 新結(jié)構(gòu)OS適宜工作壓力分析

鑒于山地地勢起伏、水肥一體化管網(wǎng)壓差變化大，四通道文丘里射流器水源入口壓力多為150 kPa～650 kPa 的實際情況，本文在新結(jié)構(gòu)OS靜壓250 kPa的研究基礎(chǔ)上，進(jìn)一步開展水源入口壓力分別為150 kPa，350 kPa，450 kPa，550 kPa，650 kPa的吸肥量對比研究，探尋新結(jié)構(gòu)OS的適宜工作壓力，為水肥一體化工程設(shè)計、施工提供技術(shù)依據(jù)。

4.1 仿真結(jié)果

分別重復(fù)仿真4次，取各通道吸肥量平均值見表2。當(dāng)壓力為150 kPa時，吸肥量最大為4 300.55 L/h；當(dāng)壓力為250 kPa時，吸肥量降為2 317.20 L/h；當(dāng)壓力為350 kPa時，吸肥量僅為849.03 L/h；當(dāng)壓力為450 kPa時，倒吸量為17 534.71 L/h；當(dāng)壓力為550 kPa時，倒吸量為1 772.52 L/h；當(dāng)壓力為650 kPa時，倒吸量為2 771.27 L/h。

表2 最優(yōu)結(jié)構(gòu)OS仿真分析結(jié)果Tab. 2 Optimal structure OS simulation analysis results

4.2 結(jié)果分析

由表2可知，當(dāng)入口壓力為150 kPa時，新結(jié)構(gòu)OS的吸肥量最大為4 300.55 L/h。當(dāng)入口壓力從150 kPa增至350 kPa時，吸肥量降低至849.03 L/h，且通道1～通道4吸肥量方差增大，表明吸肥量與入口壓力呈負(fù)相關(guān)關(guān)系。當(dāng)入口壓力為350 kPa時，吸肥通道已開始出現(xiàn)不均勻的震蕩。當(dāng)壓力進(jìn)一步增大至450 kPa時，開始出現(xiàn)倒吸現(xiàn)象，吸肥通道出現(xiàn)劇烈震蕩。隨著壓力進(jìn)一步增大，倒吸現(xiàn)象逐漸變小，震蕩現(xiàn)象逐漸減弱。研究表明新結(jié)構(gòu)OS的適宜入口壓力范圍為150 kPa～250 kPa，吸肥量最大的入口壓力為150 kPa。該研究進(jìn)一步在四通道文丘里射流器中驗證了李加念等提出的施肥裝置最佳入口壓力范圍為150 kPa～250 kPa的結(jié)論。

5 結(jié)論

1) 在文丘里射流器出口與吸肥泵連接旁路吸肥模式下，水源管道入口靜壓為250 kPa，文丘里射流器出口與吸肥泵連接，體積流量設(shè)定為0.002 2 m3/s，吸肥泵出口靜壓為450 kPa，文丘里射流器吸肥口及營養(yǎng)液出口環(huán)境壓力設(shè)定為101.325 kPa。通過向左、右和中間三個方向縮進(jìn)5 mm、10 mm、15 mm、20 mm共12個水平，改變并聯(lián)四通道文丘里射流器T型出口安裝位置達(dá)到調(diào)節(jié)能量損失的仿真分析結(jié)果表明：RI4向右縮進(jìn)20mm的結(jié)構(gòu)，吸肥量為2 317.20 L/h，基本持平或高于其他試驗組吸肥量水平，且RI4的方差為218.43低于其他試驗組，均勻性最好。故采用RI4向右縮進(jìn)20 mm為優(yōu)化后的新結(jié)構(gòu)OS。

2) 在實際應(yīng)用水肥一體化管網(wǎng)入口壓力多為150 kPa～650 kPa的條件下，對新結(jié)構(gòu)OS進(jìn)行仿真分析表明：文丘里射流器吸肥量與入口壓力呈負(fù)相關(guān)關(guān)系。當(dāng)壓力增至350 kPa時，通道1～通道4已開始出現(xiàn)不均勻的震蕩；當(dāng)壓力進(jìn)一步增大至450 kPa時，吸肥通道出現(xiàn)劇烈震蕩并伴隨出現(xiàn)倒吸現(xiàn)象。研究表明四通道文丘里射流器適宜的入口壓力為150 kPa～250 kPa。吸肥量最大的入口壓力為150 kPa，此時吸肥量為4 300.55 L/h。該仿真結(jié)果可有效指導(dǎo)水肥一體化設(shè)計及施工。