基于PWM技術的大型變量噴灌機整機水力性能研究及優化

莫錦秋，黃小林，李文濤，李彥明

基于PWM技術的大型變量噴灌機整機水力性能研究及優化

莫錦秋，黃小林，李文濤，李彥明

（上海交通大學機械與動力工程學院，上海 200240）

基于脈沖寬度調制（Pulse Width Modulation, PWM）技術的變量噴灌機雖然能實現更精細的地塊水分管理，但因實現變量的電磁閥的持續開閉動作，噴灌機主輸水管道內流量呈現持續的階躍變化，進而導致變量噴灌過程中存在壓力脈動和機械激振現象。該研究基于Matlab/Simulink，對基于PWM技術的大型變量噴灌機在變量噴灌過程中的壓力脈動進行了研究分析及改善優化。針對已研制的基于PWM技術的大型變量噴灌機實體系統，首先構建了其關鍵器件和整機的水力模型，并驗證了模型的正確性。然后基于所建模型，對變量噴灌機的壓力脈動進行了分析，得到了變量噴灌過程中的壓力脈動規律，進而提出了PWM脈沖相位錯開的緩減方法，并介紹了該方法的3種具體實施方式，即“站間”錯相、“站內”錯相和“站間+站內”錯相。最后在實現了PWM脈沖相位錯開的變量噴灌機整機水力模型上進行仿真分析，證明了所提3種錯相方式均可明顯減小變量噴灌過程中的壓力脈動幅度，同時也可有效提高泵站的利用率。此外，通過在實現了PWM脈沖相位錯開的變量噴灌機實體系統上開展田間試驗，進一步驗證了所提PWM脈沖相位錯開方法對于緩減壓力脈動和機械激振的有效性。且田間試驗結果表明，在采用PWM脈沖相位錯開方法時，變量噴灌機仍然具有較好的噴灑均勻性和變量灌溉控制誤差。該研究成果不僅能為基于PWM技術的變量噴灌提供較平穩的壓力條件從而保障變量噴灌控制精度，還能減小因壓力脈動引起的機架激振從而保障設備安全。

模型；壓力脈動；大型噴灌機；變量噴灌；脈沖寬度調制

0 引 言

傳統噴灌系統在大面積農田中只能進行均勻灌溉，為了提高水分生產率和水資源利用率，實現以農田各小區域內具體作物的實際需水量為導向的按需精準噴灌，變量噴灌應運而生[1-5]。變量噴灌主要是通過行走速度控制和定點實時調節各噴頭流量來達到變量的目的，而在實現噴頭流量調節的眾多方法中，基于脈沖寬度調制（Pulse Width Modulation，PWM）技術的變量控制是目前的主流研究與實施方法[6-9]，該方法通過改變PWM脈沖信號的占空比來控制各噴頭前端電磁閥的開啟時間比例，從而實現對各噴頭噴灑水量的獨立控制。但目前國內外基于PWM技術的變量噴灌相關研究文獻均未給出整機上各噴頭的PWM相位關系。當簡單地采用整機所有噴頭相同PWM相位時，全體電磁閥的持續同相位同周期的開閉會使得主輸水管道壓力產生大幅的周期性脈動，并引起機械激振[10]。由于基于PWM技術的大型變量噴灌機（后文簡稱大型PWM變量噴灌機）輸水管路長、整機噴頭眾多、PWM周期為10 s量級[10]，上述壓力脈動無法由供水系統的壓力反饋控制來減小。過大的主輸水管壓力脈動會改變噴頭瞬間流量，鑒于PWM變量噴灌基于確定的噴頭瞬間流量按不同通水占比來實現流量控制的機理，主輸水管壓力脈動將極大影響變量控制的精度。因此求取大型PWM變量噴灌機在變量噴灌過程中的壓力脈動規律并尋求壓力脈動緩減措施是實現準確PWM變量噴灌的前提。

目前對于噴灌機的水力性能相關研究中，多是均勻灌溉時的定常計算或是對單一部件的水力性能研究。如劉忠潮等[11-14]對噴頭等間距分布和不等間距分布情況下的基于多口系數的噴灌機管道壓力計算方法進行了研究。Kincaid等[15]研究了圓形噴灌機在定常流動下的壓力分布情況。基于多口系數的計算方法相當于采用折算系數的估算方法，結果相對粗糙。而在單一噴灌部件的水力性能相關研究中，張云龍等[16]應用計算流體力學方法對圓形噴灌機系統兩井泵匯流裝置進行了研究，求取了雷諾數與兩進口匯流比，以及匯流裝置的總水力損失與兩進口局部阻力系數之間的關系。張琛[17]采用理論分析和數值模擬方法對壓力調節器的設計和性能優化進行了研究，為壓力調節器的設計與改進提供了理論和技術支撐。上述研究均是針對部件的分析，并未揭示噴灌機整體相互作用下的水力特性。

由于變量噴灌過程中壓力脈動的時間相關性，所以動態性能分析顯得十分必要。但傳統的流體系統動態分析方法多是基于流體方程的數值求解，變量繁多，過程抽象且易出錯[18]。而基于Matlab/Simulink的建模仿真具有建模簡便直觀、計算迅速準確、可根據需要靈活自定義組件等優勢[19]，且其工具包Simscape內即有成熟的針對流體分析的組件。

本文針對已研制的大型PWM變量噴灌機實體系統，首先對其關鍵器件的水力模型進行構建，進而建立整機水力模型，并結合試驗驗證模型的正確性。基于所建模型，對PWM變量噴灌過程中的壓力脈動進行分析，擬根據脈動規律尋求緩減措施，并結合田間試驗驗證其有效性。

1 PWM變量噴灌機水力模型構建與驗證

1.1 試驗系統及方法

1.1.1 試驗系統

研制的位于河南許昌基地的大型PWM變量噴灌機實體如圖1，該變量噴灌機同時具備平移變量噴灌作業和旋轉變量噴灌作業功能，其主輸水管道由主塔車和節點塔車支承。

1.主塔車 2.主控制器 3.噴灌機進水口 4.噴灑執行件 5.主輸水管道 6.電磁閥 7.出水分管道 8.第1跨 9.1#節點控制器 10.2#節點控制器 11.節點塔車 12.第2跨 13.3#節點控制器 14.第3跨 15.4#節點控制器 16.第4跨 17.5#節點控制器

在結構上，此噴灌機共4跨（從主塔車端開始依次為第1、2、3、4跨），無尾槍，主塔車端進水，每跨50 m，在各跨主輸水管道上間隔3 m設置出水口，各出水口后接電磁閥，然后接鵝頸彎管和軟管，軟管末端則是壓力調節器和噴頭（本文合稱噴灑執行件）。整機共75個噴頭，基于分級分布式控制[10]，每個噴頭的噴灑受控于其上方電磁閥的開閉，而各電磁閥分別由離其最近的位于各節點塔車上的節點控制器獨立控制，整機共有5個節點控制器（從主塔車端開始依次為1#、2#、3#、4#、5#節點控制器），每個節點控制器均可獨立控制16路電磁閥（5#節點控制器例外，僅控制11路），同時每個節點控制器均受位于主塔車上的主控制器控制。另外，該變量噴灌機由3泵并聯的遠程水源供水系統抽取地下水來供水，且在噴灌機進水口和各跨前后端的主輸水管道上安裝了壓力表（紅旗YTN-100，量程0.6 MPa），在噴灌機進水口處安裝了壓力傳感器（星儀CYYZ11，量程0.6 MPa，輸出信號0～20 mA），均為豎直安裝，壓力傳感器與主控制器內STM32開發板（STM32 F103R8）的模擬輸入通道相連（外加采樣電阻165 Ω，使0～20 mA電流模擬量轉換為0～3.3 V電壓信號后進行12位、4 Hz AD采集），采集數據在主控制器處經串口通信送達變量控制主機。

1.1.2 試驗方法

1）噴頭流量

在大型PWM變量噴灌機實體系統上任取5個噴頭作為試驗噴頭，并去掉噴頭擋水盤，分別設置系統壓力、PWM周期和PWM占空比后，在噴頭下方放置水桶，測量所有試驗噴頭在此條件下的2 min出水量。當完成所有噴頭測量后，先將所有試驗噴頭出水量測量值取平均，再結合集水時長計算噴頭流量。

2）噴灌機各部分壓力損失

在大型PWM變量噴灌機實體系統上分別設置系統壓力、PWM周期和PWM占空比后，進行田間靜態噴灑試驗，待系統壓力穩定后，通過讀取位于各跨前后端的壓力表示數，求差得到所設條件下的噴灌機各部分壓力損失。

1.2 關鍵器件水力模型構建

變量噴灌機的關鍵水力器件從結構上按照功能或流阻特性可拆分為組件。針對各組件，根據其特性，可直接采用Simscape工具包內已有的成熟組件(或組件的組合)進行描述，或按照組件特性表達式自定義Simscape組件模塊來描述。

1.2.1 輸水管道

從結構上看，噴灌機的主輸水管道和各出水分管道可切分為直管結構、彎管結構和T型三通分水管結構。對于直管結構，本文使用Simscape中的“Resistive Pipe LP”模型，描述該模型流阻特性的為達西公式[20]。對于彎管和T型三通分水管結構，分別采用Simscape中基于局部阻力塊模型封裝的“Pipe Bend”模型和“T-junction”模型。

1.2.2 噴灑執行件

變量噴灌機噴灑執行件包含壓力調節器和噴頭，本文PWM變量噴灌機實體所用噴頭為美國尼爾森公司的Nelson D3000噴頭，由3TN系列噴嘴和檔水盤等組成。Nelson D3000噴頭在壓力一定時，不同型號的噴嘴因孔口直徑不同而瞬間通水流量不同。其中，壓力調節器是保證整機所有噴嘴入口壓力相等且為定值的關鍵器件。本文所用壓力調節器為美國尼爾森公司Nelson D3000噴頭配套的Nelson壓力調節器，按Nelson壓力調節器的產品使用技術要求，因壓力調節器內部摩擦損失，要求設計時入口壓力應當比出口要求壓力高0.035 MPa，即：

基于（1）式，使用Simscape模塊編程語言自定義壓力調節器模型。對于具體型號的噴嘴，將視作截面積固定的孔口，可直接使用Simscape中的“Fixed Orifice”模型。

1.2.3 電磁閥

電磁閥是實現變量噴灌的關鍵器件，本文變量噴灌機實體所用電磁閥為直流24 V直動式膜片常閉電磁閥，選用開靈的2W200-20型號。根據其動作機理，采用Simscape中“2-Position Valve Actuator”（二位閥執行器）模型與“2-Way Directional Valve”（二通換向閥）模型的組合來進行建模。其中二位閥執行器接收來自節點控制器的PWM脈沖信號，并根據信號幅值驅動二通換向閥進行狀態切換。此模型可保證電磁閥的開啟（或關閉）過程包含加速（或減速）和勻速階段，且開啟和關閉耗時長短可根據電磁閥自身開關特性進行設定。

1.2.4 水源供水系統

本文變量噴灌機實體系統由3泵并聯供水，各泵站均由1個變頻器和1個深井泵構成。深井泵為定排量泵，其出水流量取決于電機轉速。而變頻器參考該泵站設定的目標壓力值，根據反饋壓力，采用PID（Proportion Integration Differentiation）控制律來不斷調整深井泵電機的供電頻率，進而改變泵站的出水流量，以使壓力采集點（泵站出水口）的壓力在一定范圍內相對穩定。水源供水系統的輸出總流量可表示為

根據泵站構成，采用Simscape中“Velocity Source”（速度源）模型與“Fixed-Displacement Pump”（定排量泵）模型的組合來建立深井泵模型，其中速度源相當于深井泵電機，其轉速由PID輸出的交流電頻率決定，而PID參數取值則與實體系統相同。

1.3 噴灌機整機水力模型構建

噴灌機主輸水管道為在主塔車處豎直上升、在各跨形狀均呈拱形的鍍鋅鋼管。首先將主輸水管道拆分為直管、90°彎管和T型三通分水管結構，將各出水分管拆分為鵝頸彎管和分管直管結構。然后結合關鍵器件模型，分別構建噴灌機各跨水力模型及水源供水系統和主塔車管道水力模型，如圖2a和圖2b所示（此處僅以第2跨為例來介紹本文大型PWM變量噴灌機各跨水力模型構建結果，其他跨模型類似）。最后將水源供水系統模型、主塔車管道模型、各跨模型順次連接，同時創建可設置PWM周期及占空比等參數的各個節點控制器，以在仿真過程中按需生成控制各路電磁閥開閉的PWM脈沖，并參照實體系統，建立各節點控制器與各路電磁閥之間的具體信號連接，從而構建出針對已研制的大型PWM變量噴灌機實體系統的整機水力模型，結果如圖2c所示。根據圖1的噴灌機實體系統實際配置，相關參數取值如表1所示。

表1 模型構建相關參數取值

注：Vx_y為節點控制器x連接到其控制的各路電磁閥的輸出端口，x為節點控制器序號，y為端口序號，其中：x=1,2,…,5；y=11,12,…,18,21,22,…,28。

Note: Vx_y is the output port of the node controller x connected to each solenoid valve controlled by it, x is the serial number of the node controller, y is the serial number of the output port, where: x=1,2,…,5;y=11,12,…,18,21,22,…,28.

圖2 變量噴灌機整機水力模型構建

<

Fig.2 Construction of hydraulic model of the variable sprinkler irrigation machine

1.4 噴灌機水力模型驗證

1.4.1 噴灑執行件和電磁閥水力模型驗證

基于已建立的噴灑執行件模型，得到單一噴灑執行件在不同壓力調節器入口壓力下的噴頭流量。如圖3a所示，當壓力調節器入口壓力達到其設計入口壓力值后，壓力調節器進入正常調節狀態，此時噴頭流量隨壓力調節器入口壓力的恒定而穩定。

再結合已建立的電磁閥模型，保持壓力調節器入口壓力值為不低于其設計入口壓力的0.20 MPa、PWM周期為變量噴灌作業常用的10 s，得到單一噴頭在不同PWM占空比下的出水流量，如圖3b所示，可看出噴頭出水流量在變量區（PWM占空比低于100%）與PWM占空比有較好的線性關系，但均偏離由100%占空比時噴頭出水流量計算所得的理論值，而偏離程度由電磁閥自身的開關特性所決定。

在與仿真相同條件下，根據噴頭流量獲取方法，在變量噴灌機實體系統上進行試驗，測得了不同設置下的噴頭流量。圖3a和圖3b中試驗值與仿真值吻合，說明了所建噴灑執行件模型和電磁閥模型的合理性。

1.4.2 整機水力模型驗證

參考變量噴灌作業常用壓力，取泵站目標壓力值0.25 MPa，在整機100%PWM占空比下，仿真獲取了此時噴灌機主塔車和各跨的壓力損失，如圖3c，可看出噴灌機壓力損失主要發生在主塔車處（約0.037 MPa），各跨壓力損失從第1跨到第4跨逐漸減小。

在研制的大型PWM變量噴灌機實體系統上按照仿真相同設置進行靜態噴灑，根據噴灌機各部分壓力損失獲取方法，得到了相同條件下的噴灌機各部分壓力損失，與仿真值相近（圖3c），從而驗證了所建整機水力模型的合理性。

圖3 關鍵器件與整機水力模型驗證

2 PWM變量噴灌機壓力脈動分析及改善

2.1 壓力脈動分析

因大型PWM變量噴灌機田間噴灌作業時行走速度較低，經田間實踐發現，行走對噴灌機壓力脈動的影響不大，故對大型PWM變量噴灌機的壓力脈動研究與改善過程中的相關仿真與試驗均在變量噴灌機靜止狀態下進行。且在變量噴灌過程中，通過觀察位于噴灌機入水口和各跨前后端的壓力表示數發現，噴灌機進水口處壓力脈動幅度最大，故壓力脈動相關仿真與試驗的數據采集點均設在噴灌機進水口處。另外，根據文獻[21]，取具有較好噴灑均勻性的PWM周期為10 s，參考田間變量噴灌作業常用壓力，泵站目標壓力值取0.25 MPa。

基于所建立的整機水力模型，先設置整機PWM占空比為具有代表性的50%，得到變量噴灌時的噴灌機進水口壓力曲線、以及位于第4跨中間的壓力調節器（具有整機最低壓力調節器入口壓力）的入口壓力曲線，如圖4，可以看出，在變量噴灌過程中，噴灌機主輸水管道內存在壓力脈動，且僅發生在電磁閥打開和關閉瞬間，其中打開瞬間壓力減小，關閉后瞬間壓力增大，增大幅度遠高于減小幅度。究其原因，是由于在電磁閥快速打開和快速關閉時，噴灌機主輸水管道內存在流量突變，而由于遠程泵站采樣時間和調整時間的存在，泵站輸出流量并不能跟隨整機實際噴灌流量的變化而同步變化，即產生的流量突變并不能及時得到泵站的響應，從而導致壓力脈動的產生[22]，而且因為采樣時間和調整時間的非零性，該壓力脈動將一直存在。另外，從圖4中的壓力調節器入口壓力曲線還可看出，由于主輸水管道內的壓力脈動，在電磁閥開啟瞬間存在壓力調節器入口壓力低于其設計入口壓力的時刻，這將導致噴頭噴灑流量的偏低以及變量噴灌結果的不準確。

圖4 噴灌機進水口壓力和調壓器入口壓力

為了明確在變量噴灌過程中PWM占空比、系統作業壓力、變量作業噴頭數和供水系統壓力調節響應時間等因素對壓力脈動幅度的影響，首先獲取了在不同系統作業壓力水平、不同PWM占空比下的噴灌機進水口壓力脈動幅度（即單個PWM周期內噴灌機進水口壓力最大值與最小值的差），如圖5，可知壓力脈動幅度與PWM占空比大小無關，但隨著系統作業壓力的增大而增大。然后鑒于壓力脈動與流量突變的相關性，保持50%PWM占空比不變，得到不同變量作業噴頭數下的噴灌機進水口壓力最值，如圖6，可知變量作業噴頭數越多，即流量突變值越大，噴灌機進水口壓力最大值越大，壓力最小值越小，即壓力脈動幅度越大，且壓力脈動幅度非線性增大。最后以變化PID積分系數I為例（僅在實體系統實際值0.003周圍小量變化），對供水系統壓力調節響應時間的影響進行了仿真分析，得到50%PWM占空比時不同PID積分系數I下的噴灌機進水口壓力曲線如圖7，對比知壓力脈動幅度隨著I的增大（對應系統響應時間減小）而減小，從而可得壓力脈動幅度與供水系統壓力調節響應時間成正相關。

最后，取2#節點控制器所控制的16個噴頭進行田間變量噴灑試驗，其他噴頭采用100%PWM占空比。噴灌機靜止，設置泵站目標壓力為0.20 MPa（該壓力在滿足整機所有壓力調節器入口壓力均大于其設計入口壓力的前提下，產生的脈動激振幅度適中，對設備造成的機械損傷小）、試驗噴頭PWM周期為10 s、PWM占空比為50%，待噴灑狀態穩定后，通過安裝在噴灌機進水口處的壓力傳感器采集了此時噴灌機進水口壓力曲線，并進行中值濾波降噪處理，結果如圖8，對比圖中的仿真曲線，兩者脈動規律基本吻合，從而進一步驗證了所建模型和分析的合理性。

圖5 不同PWM占空比下的噴灌機進水口壓力脈動幅度

圖6 不同變量作業噴頭數下的噴灌機進水口壓力最值

圖7 不同PID積分系數KI下的噴灌機進水口壓力

圖8 2#節點控制器變量噴灌時的噴灌機進水口壓力

2.2 壓力脈動改善

變量噴灌過程中的壓力脈動不僅會影響噴灌結果的準確性，而且脈動幅度過大時的機械激振易對噴灌機造成致命損傷。因此，對變量噴灌機的壓力脈動進行改善是實現準確且安全的變量噴灌的前提。

優化供水系統泵站變頻器的壓力反饋調節策略，如提高PID調節的系統響應時間、采用自適應PID控制等，理論上可減小壓力脈動幅度。但由于水的遠程輸送時間、壓力脈動傳輸時間及途中損失等因素，供水系統壓力調節即使是調整到較理想狀態時，對變量噴灌機壓力脈動的改善也不明顯。

結合壓力脈動的產生原因，可通過減小變量噴灌過程中主輸水管道內的流量階躍幅度來緩減壓力脈動。即從變量噴灌具體執行策略的角度，在變量噴灌過程中減少同一時刻打開或關閉電磁閥的數量，也就是按一定策略將控制電磁閥開閉的各路PWM脈沖相位進行錯開。結合大型PWM變量噴灌機的控制系統構成，若以節點控制器為站，則該方法可具體分為“站間”錯相、“站內”錯相、以及“站間+站內”錯相。

2.2.1 錯相方式

1）“站間”錯相

以節點控制器為單位，在生成各路PWM脈沖時，同一節點控制器下的各路PWM脈沖相位相同，而不同節點控制器下的PWM脈沖則把相位錯開一定角度。錯開角度及各節點控制器的PWM脈沖相位可表示為

2）“站內”錯相

先將同一節點控制器所控制的電磁閥進行分組并順次進行編號，不同節點控制器下的對應位置電磁閥組編號相同，然后以電磁閥組為單位，在生成各路PWM脈沖時，相同編號電磁閥組內的各路PWM脈沖相位相同，不同編號電磁閥組的PWM脈沖相位則錯開一定角度。錯開角度及各電磁閥組的PWM脈沖相位可表示為

3）“站間+站內”錯相

綜合“站內”錯相和“站間”錯相，同時考慮節點控制器實現程序的可復用性，首先分別將各節點控制器下的電磁閥按相同規則進行分組，并保證每個節點控制器下所分組數相同，然后將整機所有電磁閥組順次進行編號，各節點控制器下對應電磁閥組之間編號不同。在生成各路PWM脈沖時，將整機相鄰編號電磁閥組間的PWM脈沖相位錯開一定角度，而同一電磁閥組內的各路PWM脈沖相位則相同。其中，錯開角度及各電磁閥組的PWM脈沖相位可表示為

先按“站內”錯相相同方法把各節點控制器下電磁閥都分為8組，PWM周期同取10 s，可知各電磁閥組的相位需錯開9°，即0.25 s。

2.2.2 錯相結果

在所建立的整機水力模型上，依次實施上述3種PWM脈沖錯相方式，得到各錯相方式下的噴灌機進水口壓力曲線和流量曲線，如圖9所示。

圖9 不同錯相方式下的噴灌機進水口壓力和流量

從圖9a知，3種錯相方式分別將最大壓力脈動幅度從錯相前的0.60 MPa減小到了0.09、0.05、0.01 MPa，可見3種錯相方式均可大幅減小變量噴灌過程中的壓力脈動幅度，其中“站間+站內”錯相方式下的壓力脈動幅度最小，減幅效果最明顯，但壓力脈動頻率最大；而“站間”錯相方式對應的壓力脈動幅度雖相對稍大，但脈動頻率最小。對于選用何種錯相方式，應結合噴灌機實體系統的機械振動固有頻率等因素來考慮。

從圖9b可知，使用3種錯相方式后，變量噴灌機整機所需最大噴灌流量均遠低于錯相前。其中，使用“站間+站內”錯相方式時，變量噴灌機所需最大噴灌流量最低。由此可知，PWM脈沖相位錯開方法可有效提高泵站利用率，進而降低系統構建成本。

另外，當分別采用3種錯相方式時，統計在單個PWM周期內各時刻的瞬間開閥個數和瞬間關閥個數，結果如圖10。從圖可知，錯相前由于瞬間開閥（或關閥）的數量多達75個，導致主輸水管道流量階躍幅度較大，從而產生較大幅度的壓力脈動，而使用3種錯相方式均可減小瞬間開閥（或關閥）的數量，其中“站間”錯相時最多16個，“站內”錯相時最多10個，“站間+站內”錯相時最多只有2個，進而對應了逐漸減小的主輸水管道流量階躍幅度以及逐漸減小的壓力脈動幅度。

圖10 單個PWM周期內各時刻瞬間開閥和關閥數

最后，在研制的大型PWM變量噴灌機實體系統上實現3種錯相方式后，PWM周期同取10 s、PWM占空比同取50%、變量噴灌機保持靜止，開展田間變量噴灑試驗，試驗時通過安裝在變量噴灌機進水口處的壓力傳感器進行壓力采集，并將采集數據進行中值濾波處理后得到不同條件下的噴灌機進水口壓力曲線，如圖11。

對比圖8中的噴灌機進水口壓力曲線與0.20 MPa泵站目標壓力值下、2#節點控制器采用“站內”錯相方式進行變量噴灌（其他噴頭PWM占空比設為100%）時的噴灌機進水口壓力曲線，可看出噴灌機進水口最大壓力脈動幅度從錯相前的0.07 MPa減小到了“站內”錯相時的0.025 MPa。從而證明了“站內”錯相方式的有效性。

從0.15和0.20 MPa兩個泵站目標壓力值下、整機采用“站間+站內”錯相時的噴灌機進水口壓力曲線可以看出，在整機進行變量噴灌時，變量噴灌機壓力脈動的幅度值仍然較小，且隨系統作業壓力的增大而變化不大。再次說明了所提PWM脈沖錯相方法可有效緩減變量噴灌過程中的壓力脈動。

圖11 PWM脈沖錯相后的噴灌機進水口壓力試驗結果

此外，考慮到相位錯開對變量噴灌機的噴灑均勻性和變量灌溉控制誤差產生影響的可能性，噴灌機整機采用“站內+站間”錯相方式，泵站目標壓力取0.20 MPa，在地形平坦處，按照國標GB/T 19797-2012[23]中相關要求開展田間噴灌試驗，如圖12。分別設置平移噴灌速度0.75 m/min，旋轉噴灌0.007 5 rad/min，獲取了變量噴灌機在不同占空比下的噴灑均勻性和變量灌溉控制誤差（分別按平移噴灌時的克里斯琴森水量分布均勻系數計算公式、旋轉噴灌時的赫爾曼-海因水量分布均勻系數計算公式和變量灌溉控制誤差計算公式進行計算[22,24-25]），試驗結果如表2，可以看出，此時變量噴灌機在2種噴灌模式下的噴灑均勻性均大于85%，變量灌溉控制誤差均小于5%。說明在采用PWM脈沖相位錯開方法時，變量噴灌機仍具有較好的噴灑均勻性和變量灌溉控制誤差。

圖12 噴灑均勻性和變量灌溉控制誤差測量試驗現場

表2 噴灑均勻性和變量灌溉控制誤差

3 結 論

針對大型PWM變量噴灌機在變量噴灌過程中的脈動激振問題，基于Matlab/Simulink，對研制的大型PWM變量噴灌機實體系統的水力模型進行了構建，并對變量噴灌過程中的壓力脈動進行了分析、改善和驗證。結果表明：

1）單一噴頭的實際出水量在變量區雖然與PWM占空比有較好的線性關系，但均與理論值之間存在一定偏差，該偏差值與電磁閥的開關特性有關。另外，噴灌機整機最大壓降發生在主塔車處（約0.037 MPa），單跨壓力損失沿水流前進方向逐漸減小。

2）變量噴灌過程中，在電磁閥快速打開和關閉瞬間，因主輸水管道內存在流量突變，故會產生壓力脈動。壓力脈動幅度與占空比大小無關，而與系統作業壓力水平、同時開閥或關閥數、及供水系統壓力調節響應時間正相關。

3）將控制電磁閥開閉的各路PWM脈沖相位進行錯開是減小變量噴灌時壓力脈動幅度和提高泵站利用率的有效方法。該方法從具體實施上可分為“站間”錯相、“站內”錯相和“站間+站內”錯相，三者相比而言，“站間+站內”錯相方式對壓力脈動的減幅效果最明顯，可將最大壓力脈動幅度從錯相前的0.60 MPa減小到0.01 MPa，但壓力脈動頻率稍大，而當采用“站間”錯相方式時，壓力脈動頻率最小，但減幅效果相比略差，選用時應結合實際需求和變量噴灌機本身的機械振動固有頻率等因素考慮。此外，在采用PWM脈沖相位錯開方法時，變量噴灌機在平移和旋轉2種噴灌模式下的噴灑均勻性均大于85%，變量灌溉控制誤差均小于5%。

[1] Evans R G, Han S, Schneider S M, et al. Precision center pivot irrigation for efficient use of water and nitrogen [C]//Proceedings of the 3rd International Conference on Precision Agriculture,1996: 75-84.

[2] LaRue J L. Variable rate irrigation 2010 field results [C]//American Society of Agricultural and Biological Engineers,2011: Paper No.1110787.

[3] 韓文霆，吳普特，馮浩，等. 方形噴灑域變量施水精確灌溉噴頭實現理論研究[J]. 干旱地區農業研究，2003，21(2)：105-107.

Han Wenting, Wu Pute, Feng Hao, et al. Variable-rate sprinklers for precision irrigation on square area[J]. Agricultural Research in the Arid Areas, 2003, 21(2): 105-107. (in Chinese with English abstract)

[4] Yang Qing, Pang Shujie, Yang Chenghai, et al. Variable rate irrigation control system integrated with GPS and GIS[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CSAE), 2006, 22(10): 134-138.

楊青，龐樹杰，楊成海，等. 集成GPS和GIS技術的變量灌溉控制系統[J]. 農業工程學報，2006，22(10)：134-138.(in English with Chinese abstract)

[5] 胡振方. 基于熱紅外成像技術的自走式實時變量灌溉機的設計與試驗研究[D]. 北京：中國農業大學，2014.

Hu Zhenfang. Research and Eperiment of Self-propelled Variable Irrigation System Based on Infrared Thermography[D]. Beijing:China Agricultural University, 2014.(in Chinese with English abstract)

[6] Robert G E, Jake L R, Kenneth C S, et al. Adoption of site-specific variable rate sprinkler irrigation systems[J]. Irrigation Science,2013, 31(4): 871-887.

[7] Michael D D, Calvin P. Uniformity testing of variable-rate center pivot irrigation control systems[J]. Precision Agriculture, 2006, 7(3): 205-218.

[8] 黃小林，李文濤，周志宇，等. 平移式變量噴灌機整機有風噴灌水深分布求取與修正[J]. 機械設計與研究，2019，35(5)：141-145.

Huang Xiaolin, Li Wentao, Zhou Zhiyu, et al. Calculation and correction on application depth distribution of the whole lateral move variable rate irrigation system in windy condition[J]. Mechanical Design and Research, 2019, 35(5): 141-145.(in Chinese with English abstract)

[9] 趙偉霞，李久生，栗巖峰. 大型噴灌機變量灌溉技術研究進展[J]. 農業工程學報，2016，32(13)：1-7.

Zhao Weixia, Li Jiusheng, Li Yanfeng. Review on variable rate irrigation with continuously moving sprinkler machines[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CSAE), 2016, 32(13): 1-7.(in Chinese with English abstract)

[10] 陶帥，李文濤，黃小林，等. 大型噴灌機噴頭PWM變量脈沖的實現[J]. 農機化研究，2021，43(5)：18-22.

Tao Shuai, Li Wentao, Huang Xiaolin, et al. The realization of PWM variable pulse for the sprinkler of large-scale irrigation machine[J]. Journal of Agricultural Mechanization Research, 2021, 43(5): 18-22. (in Chinese with English abstract)

[11] 劉忠潮. 計算噴灌管道沿程水頭損失的多口系數問題[J].節水灌溉，1982(2)：44-46.

[12] Edmar J, Scalopi. Adjusted F factor for multiple-outlet pipes[J]. Journal of Irrigation and Drainage Engineering, 1988, 114(1): 169.

[13] 朱廣軍. 大跨距平移噴灌機結構穩定性與灌水均勻性分析[D]. 合肥：安徽農業大學，2016.

Zhu Guangjun. The Structural Stability and Irrigation Uniformity Analysis of Large-span Translatison Prinkler Machine[D]. Hefei: Anhui Agricultural University, 2016.(in Chinese with English abstract)

[14] 魏永曜. 多支點噴灌系統管道的水力計算[J]. 節水灌溉，1982(4)：13-20.

[15] Kincaid D C, Heermann D F. Pressure distribution on a center-pivot sprinkler irrigation system[J]. Amer Soc Agr Eng Trans Asae, 1970(5): 556-558.

[16] 張云龍，嚴海軍，王建東. 圓形噴灌機系統兩井泵匯流裝置水力性能研究[J]. 灌溉排水學報，2019，38(4)：48-54.

Zhang Yunlong, Yan Haijun, Wang Jiandong. Hydraulic characteristics at the confluence of two submersible pumps in center-pivot Irrigation[J]. Journal of Irrigation and Drainage, 2019, 38(4): 48-54.(in Chinese with English abstract)

[17] 張琛. 灌溉系統直動式壓力調節器動力學模型與性能優化研究[D]. 北京：中國農業大學，2016.

Zhang Chen. Dynamic Model and Optimization of Regulating Performance of Direct-acting Pressure Regulator for Irrigation System[D]. Beijing: China Agricultural University, 2016. (in Chinese with English abstract)

[18] 姚磊. 基于Simulink的流體管網仿真[J].控制與信息技術，2010(5)：5-8，22.

Yao Lei. Simulation of fluid pipe net based on simulink[J]. Control and Information Technology, 2010(5): 5-8, 22.(in Chinese with English abstract)

[19] 田樹軍，張宏. 液壓管路動態特性的Simulink仿真研究[J].系統仿真學報，2006，18(5)：1136-1138，1146.

Tian Shujun, Zhang Hong. Study on simulation of hydraulic pipelines’dynamic characteristic by simulink software[J]. Journal of System Simulation, 2006, 18(5): 1136-1138, 1146.(in Chinese with English abstract)

[20] 洪明，李援農，馬英杰，等. 低壓條件下滴灌毛管水頭損失試驗研究[J]. 灌溉排水學報，2010，29(1)：50-52，73.

Hong Ming, Li Yuannong, Ma Yingjie, et al. Head loss experimental analysis of trickle irrigation under the low pressure[J]. Journal of Irrigation and Drainage, 2010, 29(1): 50-52, 73.(in Chinese with English abstract)

[21] 周志宇，陶帥，莫錦秋. 圓形PWM變量噴灌機噴灌特性的仿真研究[J]. 農機化研究，2020，42(5)：15-24.

Zhou Zhiyu, Tao Shuai, Mo Jinqiu. Simulation research on irrigation features of PWM center pivot sprinkler irrigation system[J]. Journal of Agricultural Mechanization Research, 2020, 42(5): 15-24.(in Chinese with English abstract)

[22] 王勇. 供水系統水錘數值計算及動態模擬[D]. 合肥：合肥工業大學，2009.

Wang Yong. Numerical Calculation and Dynamic Simulation of Water Hammer in Water Supply System[D]. Hefei: Hefei University of Technology, 2019.(in Chinese with English abstract)

[23] GB/T 19797-2012，農業灌溉設備中心支軸式和平移式噴灌機水量分布均勻度的測定[S].

[24] ASABE Standards. S436.1, Test procedure for determining the uniformity of water distribution of center pivot and lateral move irrigation machines equipped with spray or sprinkler nozzles[S]. St Joseph, Mich: ASABE, 2007.

[25] 莫錦秋，周志宇，陶帥，等. 基于PWM技術的平移式變量噴灌機噴頭流量分配方法[J]. 農業機械學報，2018，49(10)：163-171.

Mo Jinqiu, Zhou Zhiyu, Tao Shuai, et al.Sprinkler flow rate distribution method in lateral moving variable rate irrigation system based on PWM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(10): 163-171. (in Chinese with English abstract)

Research and optimization of hydraulic characteristics of large-scale variable sprinkler irrigation machine based on PWM technology

Mo Jinqiu, Huang Xiaolin, Li Wentao, Li Yanming

(,,200240,)

Although the variable sprinkler irrigation machine based on pulse width modulation (PWM) technology can achieve more sophisticated water management in plots, the continuous opening and closing actions of the solenoid valve that realise variable rate irrigation make the flow rate of the main pipeline of the variable sprinkler irrigation machine a continuous and stepped change, which in turn leads to pressure pulsation and mechanical excitation in the variable sprinkler irrigation process. Based on Matlab/Simulink, the pressure pulsation of the large-scale variable sprinkler irrigation machine based on PWM technology in the variable sprinkler irrigation process was analyzed and improved. Regarding the developed large-scale variable sprinkler irrigation machine based on PWM technology, the hydraulic model of the key components and the whole machine were constructed first, and the correctness of the models were verified by experiments. Then based on the built hydraulic models of the whole machine, the pressure pulsation of the variable sprinkler machine was analyzed, and the law of the pressure pulsation in the variable sprinkler irrigation process was obtained. The results indicated that the amplitude of pressure pulsation had no effects on the duty cycle of the PWM pulse while was positively related to the operating pressure level of the system, the number of valves opened or closed at the same time, and the response time of the pressure adjustment of the water supply system. Furthermore, in order to reduce the number of valves that are opened or closed at the same time, the mitigation method of PWM pulse phase staggering was proposed, and three specific implementations of the method had been introduced, including staggering phase between stations, staggering phase in the station, and staggering phase between stations and within the station. Finally, the simulation analysis of the whole machine with implementing the PWM pulse phase staggering method on the model proved that three proposed staggering methods could significantly reduce the pressure pulsation amplitude in the variable sprinkler irrigation process, and also effectively improve the utilization rate of the pumping station. In particularly, the method that staggering phase between stations and within the station had the most obvious effectiveness on reducing pressure pulsation, which can reduce the maximum pressure pulsation amplitude from 0.60 MPa before phase staggering to 0.01 MPa. In addition, field experiments after realizing the PWM pulse phase staggering method on the developed large-scale variable sprinkler irrigation machine based on PWM technology further verify the effectiveness of the proposed PWM pulse phase staggering method in reducing pressure pulsation and mechanical excitation. The results in field experiment showed that the variable sprinkler irrigation machine still has good spray uniformity and variable irrigation control error when using the PWM phase staggering method. These research results can not only provide relatively stable pressure conditions for variable rate irrigation based on PWM technology to ensure the control accuracy in the variable sprinkler irrigation process, but also reduce the mechanical excitation caused by pressure pulsation to ensure the safety of the equipment.

models; pressure pulsation; large-scale sprinkler irrigation machine; variable rate irrigation; pulse width modulation

莫錦秋，黃小林，李文濤，等. 基于PWM技術的大型變量噴灌機整機水力性能研究及優化[J]. 農業工程學報，2020，36(19)：76-85.doi：10.11975/j.issn.1002-6819.2020.19.009 http://www.tcsae.org

Mo Jinqiu, Huang Xiaolin, Li Wentao, et al. Research and optimization of hydraulic characteristics of large-scale variable sprinkler irrigation machine based on PWM technology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(19): 76-85. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.19.009 http://www.tcsae.org

2020-07-19

2020-09-20

國家重點研發計劃項目(2017YFD0700504-2)

莫錦秋，博士，副教授，主要從事機電系統集成控制研究。Email：mojinqiu@sjtu.edu.cn

10.11975/j.issn.1002-6819.2020.19.009

S277.94

A

1002-6819(2020)-19-0076-10