基于ANSYS-ADAMS的立式油菜割曬機鋪放角形成機理

李海同，吳崇友，沐森林，關卓懷，江 濤

基于ANSYS-ADAMS的立式油菜割曬機鋪放角形成機理

李海同，吳崇友※，沐森林，關卓懷，江 濤

（農業(yè)農村部南京農業(yè)機械化研究所，南京 210004）

為提高油菜割曬作業(yè)性能，研究油菜植株參數(shù)和割曬機參數(shù)對鋪放角的影響規(guī)律，該研究利用ANSYS軟件和ProE軟件分別建立油菜植株的柔性體和割曬機模型，在ADAMS軟件中建立割曬機-油菜植株剛柔耦合模型，開展油菜植株鋪放過程單因素與多因素仿真試驗。單因素試驗結果表明：油菜植株高度在1.2～1.6 m范圍內，鋪放角先減小后增加且在植株高度為1.4 m時取得最小值；油菜植株質量密度在438～588 kg/m3范圍內，鋪放角先增加后減小且在質量密度為538 kg/m3處取得最大值；植株與植株的摩擦系數(shù)在0.1～0.5范圍內，摩系數(shù)越大鋪放角越小。多因素響應面試驗結果表明：輸送鏈與前進速度的速比、前進速度、割茬高度、輸送鏈速度與前進速度的速比及前進速度與割茬高度的交互作用對鋪放角的影響極顯著（<0.01），顯著性順序從大到小為輸送鏈速度與前進速度的速比、前進速度、割茬高度；根據多因素試驗結果建立鋪放角的三元二次回歸方程模型，并建立響應面分析多因素交互作用對鋪放角的影響規(guī)律；建立目標優(yōu)化方程組，確定割曬機最優(yōu)作業(yè)參數(shù)組合為前進速度0.6 m/s，速比1.37，割茬高度0.42 m。田間試驗結果表明，優(yōu)化作業(yè)參數(shù)條件下，鋪放角預測值與實測值的相對誤差小于6%；與對照相比，鋪放角由134.9°下降為115.8°，角度差由11.2°下降為9.6°，根差由0.22 m下降為0.14 m，割曬損失率由2.2%下降為1.5%，割曬機作業(yè)性能提高。研究結果可為油菜割曬作業(yè)的適宜條件確定和割曬機作業(yè)參數(shù)優(yōu)化提供依據。

農業(yè)機械；收獲；油菜；鋪放角；剛柔耦合模型；響應面分析；優(yōu)化

0 引 言

油菜機械化收獲主要有分段收獲和聯(lián)合收獲2種方式[1-3]，分段收獲采用先割曬鋪放，再撿拾脫粒的作業(yè)流程，可充分利用油菜的后熟特性，既保證了脫粒前籽粒成熟度一致，又實現(xiàn)了油菜籽充分完熟，保障了菜油的品質[4-7]，具有收獲期長、對油菜的成熟度及其一致性和株型等不敏感、收獲損失率小等特點[8-10]。鋪放角是評價油菜割曬鋪放質量的重要指標之一，可反映割曬機的作業(yè)性能，研究割曬參數(shù)對鋪放角的影響規(guī)律可為割曬機設計與優(yōu)化提供依據。

李平等[11-12]設計了集成橫向和縱向輸送的組合式莖稈輸送裝置的4SY-1.8改進型油菜割曬機，田間試驗表明該機收獲直播油菜時，上、下層油菜莖稈鋪放角分別為23.6°和17.1°，角度差小于10°，得到影響割曬機田間作業(yè)質量的主要因素，并確定了油菜割曬機的最佳參數(shù)范圍。金誠謙等[13-14]研制了與聯(lián)合收割機底盤配套的4SY–2型油菜割曬機，確定了單帶式輸送器、橫向輸送機構與縱向輸送機構的參數(shù)，田間試驗表明割曬作業(yè)損失率較低，鋪放角度差為6.77o。石增祥等[15]對油菜莖稈在割曬機鋪放輸送裝置上的運動情況和梳理推送拋送條件進行理論分析，得出鋪放角由油菜生物學特性和割曬機結構和運行參數(shù)決定。王修善等[16]通過分析油菜植株有序鋪放條件，從理論上得出立式割臺油菜割曬機前進速度和輸送速度對鋪放角的影響規(guī)律，田間試驗結果表明鋪放角為98.2°，角度差10.5°，割茬高度350 mm。現(xiàn)階段對油菜割曬機鋪放角的研究，主要是利用理論分析與田間試驗相結合方法，探究割曬機參數(shù)和油菜物料參數(shù)對鋪放角的影響，但由于田間試驗受收獲季節(jié)限制，國內外對鋪放角的研究少見報道，缺少割曬作業(yè)條件和參數(shù)優(yōu)化的理論依據。

隨著計算機技術和軟件水平的發(fā)展，計算機輔助設計及模擬仿真已作為一種重要的研究方法應用于相關領域[17-18]。利用計算機開展割曬機仿真試驗，可突破季節(jié)條件限制，提高試驗的速度和精度，有利于縮短研發(fā)周期，節(jié)約試驗成本。為此，本文通過分析立式油菜割曬機鋪放過程建立鋪放角的數(shù)學模型，分析鋪放角的影響因素；利用有限元軟件ANSYS建立油菜植株的柔性體模型，并導入ADAMS中建立割曬臺-油菜剛柔耦合模型；開展單因素與多因素響應面仿真試驗，分別研究油菜植株物料參數(shù)和割曬機作業(yè)參數(shù)對鋪放角的影響規(guī)律，并進行田間驗證試驗，以期為油菜割曬適宜條件選擇和割曬機作業(yè)參數(shù)優(yōu)化提供依據。

1 割曬機結構與工作原理

1.1 割曬機整機結構

手扶立式割臺油菜割曬機主要由機架、行走底盤、傳動系統(tǒng)、水平切割器、側邊切割器、橫向輸送鏈組和分行器等部件組成，整機結構如圖1所示。行走底盤為履帶式手扶拖拉機，割臺掛接在底盤前方并可調節(jié)割茬高度，側邊切割器安裝于機架左側將待割區(qū)與未割區(qū)交錯的油菜切割分離開來，橫向輸送鏈組平行安裝于機架上、中、下部，分行器總成安裝于機架前端并向前伸出。

1.機架 2.割茬高度調節(jié)機構 3.行走底盤 4.發(fā)動機 5.鏈條 6.變速箱 7.割臺機架 8.側邊切割器 9.橫向輸送鏈 10.水平切割器 11.分行器 12.地面

1.2 工作原理

立式割臺割曬機作業(yè)時，側邊切割器將割曬機正前方待割區(qū)與未割區(qū)交錯纏繞的分枝剪斷分開并推向兩側；機器正前方的油菜在分行器作用下運動至割臺，被水平切割器切斷后在分行器和上、中、下3條橫向輸送鏈的夾持輸送作用下向割臺右側輸送，脫離橫向輸送鏈后油菜植株在自重和橫向輸送鏈與機組縱向前進速度所形成的慣性力作用下向割臺右側連續(xù)鋪放，形成側邊條鋪，便于晾曬。

1.3 主要技術參數(shù)

立式割臺油菜割曬機為手扶懸掛式，參考現(xiàn)有割曬機相關標準及后續(xù)油菜撿拾機作業(yè)要求，其主要技術參數(shù)為：發(fā)動機功率7.5 kW，割副1.2 m，前進速度0.6～1.0 m/s，橫向輸送鏈速度0.8～1.4 m/s，割茬高度0.2～0.5 m。

2 鋪放過程與鋪放角分析

2.1 鋪放過程植株受力及運動過程分析

脫離割臺至鋪放成行的過程中，油菜植株的運動過程中與割臺的相對位置關系如圖2所示。

1.油菜植株 2.油菜割曬機

1.Rape plant 2.Rape windrower

注：為坐標原點，、為坐標軸；1為油菜植株落地點；為油菜植株的初始傾角，(°)；為割茬高度，mm；為油菜植株的轉動角速度，rad·s-1；為鋪放角，(°)；0為油菜植株的初始位置；1為油菜植株與地面接觸時的位置；2為油菜植株的最終鋪放位置。

Note:is coordinate origin;、are axis;1is the landing point of rape plant;is the initial dip angle of rape plant, (°);is stubble height, mm;is the angular velocity of rape plant, rad·s-1;is laying angle, (°);0is the initial position of rape plant;1isthe position of rape plant in contact with ground;2isthe final laying position of rape plant.

圖2 油菜植株與割曬機位置關系示意圖

Fig.2 Schematic diagram of position relationship between rape plant and windrower

根據油菜植株的運動過程和受力情況，鋪放過程分為2個連續(xù)的過程：傾角為的油菜植株脫離割曬機作用至底端接觸地面的三維平動，即由位置0至位置1的下落過程；油菜植株以落地點為旋轉點繞底端轉動至鋪放于地面的過程，即由位置1至位置2的轉動過程。

2.1.1 油菜植株下落過程分析

油菜植株在下落過程中的運動為三維空間內的平拋運動，以植株離開割臺瞬間最下端在地面的投影點為坐標原點建立笛卡爾坐標系，軸正方向為割曬機前進方向，軸垂直向上，軸正方向為橫向輸送速度方向。油菜植株平拋運動在3個方向的運動分別為：方向的與機器前進速度相等的勻速運動，方向的初速度為0的自由落體運動，方向的與橫向輸送速度相等的勻速運動。依據文獻[16]可知，油菜植株下落過程結束時植株底端在3個方向上的坐標和速度分別為

式中x、y、z為油菜植株底端在3個方向的坐標；v為割曬機前進速度，m/s；v為橫向輸送鏈速度，m/s；為油菜植株底端與地面的距離，即割茬高度，m；v為油菜植株底端在方向上的速度，m/s；v為油菜植株底端在方向上的速度，m/s；v為油菜植株底端在方向上的速度，m/s；為重力加速度，9.8 m/s2。

由此可得油菜下落過程的運動時間1、植株的合速度1與軸的夾角1分別為

式中1為油菜植株下落時間，s；

油菜植株底端與地面接觸時發(fā)生完全非彈性碰撞，碰撞后油菜植株由空間剛體平動變成繞落地點1的轉動。碰撞過程中油菜植株的動量守恒方程為

1-0sin=0 （6）

由此求得碰撞后油菜植株的角速度

式中為植株的質量，kg；0為質心到底端的距離，m，為植株對1點的轉動慣量，kg·m2。

由下落過程分析可以得到下落過程結束時油菜植株的位置、速度和角速度等狀態(tài)參數(shù)，確定植株轉動過程分析的初始條件。

2.1.2 油菜植株轉動過程分析

油菜植株與地面碰撞后繞落地點1轉動，割曬機連續(xù)作業(yè)時由于相互牽連的油菜植株存在相對運動而產生摩擦力，速度較小的油菜植株對速度較大的油菜植株產生的摩擦牽連力，已鋪放油菜植株對處于轉動過程中的油菜植株的摩擦力為阻力。因此轉動過程中作用于油菜植株的力包括重力、牽連力1和阻力2。

重力作用于油菜植株質心，豎直向下；牽連力1作用于油菜角果層交錯重疊區(qū)域，方向與植株運動速度相反，大小由重疊系數(shù)、摩擦系數(shù)等參數(shù)決定；阻力2作用于已鋪放油菜的分枝與植株分枝交錯重疊的區(qū)域，方向與機組前進方向相反，大小由重疊系數(shù)、壓力、摩擦系數(shù)等參數(shù)決定。轉動過程中油菜植株的受力分析如圖3所示。

油菜植株轉動過程中牽連力1和阻力2的大小由多個因素決定，因此利用式（8）的多元方程表達其大小。

式中為包含摩擦系數(shù)1、重疊系數(shù)2、時間等參數(shù)的多維矩陣；1為到1的映射法則；為包含摩擦系數(shù)1、壓力2、時間等參數(shù)的多維矩陣；2為到2的映射法則。

由油菜植株繞1點定軸轉動動能定理得

式中1為角果層與1的距離，m。

由角動量定理可得

由式（11）可得油菜植株在轉動過程中的運動時間(2)為

注：v為油菜植株運動速度，m·s-1；γ為油菜植株運動速度與植株主干線的夾角，(°)；α為油菜植株速度與y軸的夾角，(°)；G為油菜植株重力，N；F1為牽連力，N；F2為阻力，N。

由轉動過程分析得到油菜植株的受力、轉動角速度和運動時間等參數(shù)，為計算油菜植株鋪放角提供依據。

2.2 鋪放角及其影響因素分析

由油菜植株轉動過程動力學分析可知，油菜在重力和牽連力1作用下轉動，機組前進方向上在牽連力1和阻力2作用下做減速運動，設其當量加速度大小為a。定義鋪放角為割倒的油菜植株主莖稈與機器前進方向在水平面內的后夾角，則油菜植株的理論鋪放角t為

式中為割下部分的油菜植株高度，m。

將式（7）~（8）、式（10）、式（12）代入式（13），得到鋪放角的理論計算公式為

式中為鋪放角理論值t的正弦值，其表達式為

由式（15）可知，鋪放角的影響因素包括割后的油菜植株高度、質量、轉動慣量、牽連力1、阻力2、質心到底端的距離0等物料參數(shù)和機器前進速度v、橫向輸送鏈速度v、割茬高度等割曬機作業(yè)參數(shù)。

為研究油菜植株參數(shù)和割曬機參數(shù)對鋪放角的影響規(guī)律，運用ANSYS和ADAMS軟件進行油菜植株鋪放過程的模擬仿真試驗。根據研究結果選擇適宜的割曬條件和割曬機的最佳作業(yè)參數(shù)，提高油菜割曬作業(yè)性能。

3 剛柔耦合模擬仿真試驗

3.1 剛柔耦合多體動力學模型構建

為準確描述油菜植株的外形和力學特征，利用ANSYS軟件構建油菜的柔性體力學模型并生成中間文件，在ProE軟件中建立割曬機三維模型，將2個模型導入ADAMS中構建剛-柔耦合多體動力學仿真模型。

3.1.1 油菜植株柔性力學模型構建

根據文獻[19-23]中油菜植株表型參數(shù)和機械物理特性，建立油菜植株的三維模型，由主莖稈和3個圓周方向均勻分布的分枝組成。以青油14號油菜為例，依據文獻[24-25]，油菜植株的力學參數(shù)為質量密度438 kg/m3、彈性模量1.34×1011Pa，泊松比0.3，選用Solid（Brick8node45）為材料單元，出于計算精度和時間的考慮，定義最小網格邊長為2 mm對模型進行網格劃分，生成73 275個網格，分別在植株底端和頂端創(chuàng)建2個連接點。利用ANSYS與ADAMS的接口Solve and create export file to ADAMS輸出包含莖稈形態(tài)、力學特征等參數(shù)的柔性體模型并保存為.mnf文件。

3.1.2 剛柔耦合模型構建

在ProE中建立油菜割曬機的三維模型，經簡化處理后導入ADAMS中，通過ADAMS的Flexible Bodies選項在指定目錄下打開已生成的mnf文件，構建割曬機-油菜植株剛柔耦合模型。

3.1.3 定義約束、作用力和驅動

根據油菜割曬機工作的實際情況和鋪放過程理論分析，在割曬機-油菜植株組成的剛柔耦合動力學仿真模型中依次定義約束、作用力和驅動。

1）在油菜植株與輸送鏈之間添加Forces，接觸類型為Flex Boby to Solid（柔性體與實體），根據油菜植株參數(shù)和文獻[26-27]，接觸力參數(shù)設置為：Normal Force的類型選為impact，Stiffness定義為1.0E+5，F(xiàn)orce Exponent定義為2.2，Damping定義為100，Penetration Depth定義為0.1。摩擦力Friction Force 選取為None。

2）牽連力1和阻力2影響因素較為復雜，仿真中定義阻力作為兩者的合力，大小為植株質量和油菜植株與植株的摩擦系數(shù)v的乘積，即=v。

3）定義割臺與地面沿機器前進方向為滑移副，并添加直線驅動；定義橫向輸送撥齒在輸送鏈方向為滑移副，并添加直線驅動。

3.1.4 定義關鍵傳感器

依據鋪放過程分析，定義2個傳感器模擬植株的運動過程。

1）油菜底端與地面位置傳感器1：用于檢測植株底端與地面的距離，檢測值為0時通過腳本控制使該傳感器失效，植株與地面發(fā)生碰撞，下落運動過程結束。

2）油菜頂端與地面位置傳感器2：用于檢測植株頂端與地面的距離，檢測值為0時通過腳本控制使該傳感器失效，植株鋪放于地面上，轉動過程結束。

3.1.5 仿真過程與后處理

在ADAMS中創(chuàng)建通過2個連接點的線段，根據鋪放角的定義可知，線段與軸負方向的夾角即為鋪放角，如圖4所示。

注：A為油菜植株頂端連接點，B為油菜植株底端連接點。

由油菜植株運動過程分析可知，鋪放角仿真值β的計算公式為

式中x為點在軸上的坐標；x為點在軸上的坐標；z為點在軸上的坐標；z為點在軸上的坐標。

創(chuàng)建Measure測量仿真過程中線段與軸負方向的夾角。設置仿真時間和步長，點擊開始進行鋪放過程仿真計算。仿真完成后進入ADAMS的Post process后處理模塊提取仿真結果，繪制鋪放角隨時間的變化曲線，如圖5所示。

圖5 鋪放角隨時間變化曲線

3.2 試驗設計

3.2.1 單因素試驗

為探究適合割曬的油菜植株參數(shù)，以油菜植株的物料參數(shù)為因素開展單因素試驗，研究其對鋪放角的影響規(guī)律。根據油菜割曬機的技術參數(shù)和行業(yè)標準JB/T 7733—2007《割曬機技術條件》，結合試驗油菜植株特性，單因素試驗作業(yè)參數(shù)設置為：前進速度0.8 m/s、橫向輸送鏈與前進速度的速比1.2、割茬高度0.4 m。

根據文獻[28]和油菜植株參數(shù)實際測試結果，單因素試驗的因素水平為：油菜植株高度1.2、1.3、1.4、1.5、1.6 m，油菜植株質量密度388、438、488、538、588 kg/m3，油菜植株與植株的摩擦系數(shù)0.1、0.2、0.3、0.4、0.5。

3.2.2 多因素參數(shù)優(yōu)化試驗

為研究割曬機作業(yè)參數(shù)對鋪放角的影響規(guī)律并確定最優(yōu)參數(shù)組合，根據理論分析和文獻[16]，以鋪放角為響應值，以割曬機前進速度1、橫向輸送鏈與前進速度的速比（以下稱“速比”）2和割茬高度3作為影響因子開展響應面試驗。

根據Box-Behnken中心組合設計理論，以鋪放角為響應值，以前進速度1、速比2和割茬高度3為因素開展響應面試驗[29-31]。試驗因素水平編碼如表1所示。

表1 因素水平編碼表

4 試驗結果與分析

4.1 單因素試驗結果與分析

4.1.1 油菜植株高度對鋪放角的影響

圖6為油菜植株高度對鋪放角的影響曲線。由圖6可知，油菜植株高度為1.2～1.4 m時，鋪放角隨植株高度的增加而減少，植株高度在1.4～1.6 m范圍內鋪放角隨植株高度的增加而增加，植株高度為1.4 m時，鋪放角取得最小值120.4°。

圖6 油菜植株高度對鋪放角的影響

4.1.2 油菜植株質量密度對鋪放角的影響

圖7為油菜植株質量密度對鋪放角的影響曲線。由圖7可知，油菜植株質量密度小于538 kg/m3時鋪放角隨質量密度增加而增加，質量密度繼續(xù)增加時鋪放角減小，植株質量密度為538 kg/m3時，鋪放角取得最大值133.6°。

圖7 油菜植株質量密度對鋪放角的影響

4.1.3 油菜植株與植株的摩擦系數(shù)對鋪放角的影響

圖8為油菜植株與植株的摩擦系數(shù)對鋪放角的影響曲線。由圖8可知，油菜植株與植株的摩擦系數(shù)為0.1～0.5時，鋪放角隨油菜植株與植株的摩擦系數(shù)增加而減少，變化范圍為117.6°～134.3°。

圖8 油菜植株與植株的摩擦系數(shù)對鋪放角的影響

單因素仿真試驗結果表明，油菜植株高度為1.4 m時鋪放角取得最小值120.4°，質量密度為538 kg/m3時鋪放角取得最大值133.6°，油菜植株與植株的摩擦系數(shù)為0.5時鋪放角取得最小值117.6°。

4.2 多因素響應面試驗結果與分析

4.2.1 回歸模型與顯著性分析

多因素試驗設計方案與結果如表2所示，利用 Design-Expert 軟件[32-33]對試驗結果進行顯著性分析，結果如表3所示。

表2 多因素試驗設計方案與結果

由表3方差分析結果可知，方程評價指標鋪放角仿真值的值為0.000 2＜0.01，表明回歸模型極顯著，其決定系數(shù)2值為0.97，表明模型與試驗值擬合程度較好。其中，1、2、3、13、12的值小于0.01，說明對模型影響極為顯著；12、13、32的值大于0. 05，說明對模型影響不顯著，模型中可不予考慮。根據值可知各參數(shù)對分布變異系數(shù)由大到小順序為：2、1、3，即速比、前進速度、割茬高度。去除不顯著因素，可得鋪放角的回歸函數(shù)值r與前進速度1、速比2、割茬高度3的三元二次回歸方程為

根據回歸方程（16），分別對3個自變量求導可得

4.2.2 交互作用對鋪放角影響分析

通過 Design-Expert 軟件生成響應面曲線如圖9，根據響應面分析前進速度、速比和割茬高度的交互因素對鋪放角的影響。

表3 回歸方程方差分析

注：< 0.01（極顯著）；0.01≤< 0.05（顯著）；≥0.05（不顯著）。

Note:< 0.01（high significance）；0.01≤< 0.05（significant）；≥0.05（not significant）。

圖9 交互作用對鋪放角的影響

圖9a為割茬高度位于中心水平（0.35 m）時，前進速度與速比對鋪放角交互作用的響應面曲線。由圖9a可知，前進速度增大時，鋪放角先增大后減小，增加幅度先增大后減小；鋪放角隨速比的增加而減小，減小幅度較為平緩。在前進速度和速比對鋪放角的交互作用中，前進速度對交互作用的影響較大。

圖9b為速比位于中心水平（1.2）時，前進速度與割茬高度對鋪放角交互作用的響應面曲線。由圖9b可知，前進速度增大時，鋪放角先增大后減小，增加幅度先增大后減小；割茬高度增加時，鋪放角呈增加趨勢，但增加幅度較小。在前進速度和割茬高度對鋪放角的交互作用中，前進速度對交互作用的影響較大。

圖9c為前進速度位于中心水平（0.8 m/s）時，速比與割茬高度對鋪放角交互作用的響應面曲線。由圖9c可知，速比增加時鋪放角減小，減小幅度逐步變小；鋪放角隨割茬高度的增加而增加，增加幅度較為平緩。在速比和割茬高度對鋪放角的交互作用中，速比對交互作用的影響較大。

4.2.3 參數(shù)優(yōu)化

割曬機前進速度、速比和割茬高度及其交互作用對鋪放角的影響規(guī)律各不相同，為了尋求油菜割曬機最佳工作參數(shù)組合，提高割曬作業(yè)性能，進行參數(shù)優(yōu)化設計。

參考行業(yè)標準JB/T 7733—2007《割曬機技術條件》的技術要求，鋪放角越接近90°越好，建立目標函數(shù)與各參數(shù)變量的約束條件如式（21）所示：

運用Design-Expert 中的Optimization模塊對目標方程進行優(yōu)化求解，得出割曬機的最優(yōu)參數(shù)組合為：前進速度0.6 m/s，速比1.37，割茬高度0.42 m，此時的鋪放角為122.6°。

4.3 田間試驗

為驗證模型預測的準確性和油菜割曬機最優(yōu)參數(shù)組合條件下的作業(yè)性能，于2019年9月21日—23日在甘肅省張掖市民樂縣進行田間性能試驗。試驗油菜品種為青油14號，種植方式為直播，油菜種植密度33 株/m2，植株平均高度1.6 m，底莢高度0.56 m，無倒伏，油菜籽粒及莖稈含水率分別為19.8%和58.9%，籽粒千粒質量3.74 g。

根據多因素響應面試驗優(yōu)化結果設置對照組T1與優(yōu)化組T0，參考行業(yè)標準JB/T 7733—2007《割曬機技術條件》依據作業(yè)經驗，對照組試驗條件設定為：前進速度1.0 m/s、速比1.2、割茬高度0.5 m；根據參數(shù)優(yōu)化結果優(yōu)化組試驗條件設置為：前進速度0.6 m/s、速比1.4、割茬高度0.4 m，選擇地勢平坦、油菜長勢均勻、無雜草的田塊，在風力不大于3 m/s的情況下開展油菜割曬作業(yè)，每次作業(yè)長度30 m，重復3次。

割曬作業(yè)完成后，參考標準JB/T 7733—2007《割曬機技術條件》和GB/T 8097—2008《收獲機械聯(lián)合收割機試驗方法》測量割曬作業(yè)各項指標。具體測量方法為：將皮卷尺平行固定于鋪放條上，每隔2 m選取1 m的測試區(qū)，利用角度尺測量測試區(qū)內油菜主莖稈與皮卷尺的夾角，利用鋼卷尺測量油菜植株底端到皮卷尺的距離，記錄數(shù)據并求平均值。則油菜植株鋪放角為

式中β為每個測區(qū)內鋪放角的平均值，(°)。

角度差為

式中為Δ角度差，(°)；maxβ為每個測區(qū)內鋪放角的最大值，(°)；minβ為每個測區(qū)內鋪放角的最小值，(°)。

根差為

式中為Δ根差，m；maxρ為測區(qū)內油菜植株底端到皮卷尺距離的最大值，m；minρ為測區(qū)內油菜植株底端到皮卷尺距離的最小值，m。

每隔2 m撿拾1.2 m×1 m面積內的落粒、掉穗，脫粒后稱其籽粒質量，減去每平方米自然落粒得出每個測區(qū)內油菜割曬損失質量，則割曬損失率為

式中g為割曬損失率，%；w為每平方米內損失的籽粒質量，kg；為油菜產量，kg/hm2。

根據試驗測試與計算，鋪放角、角度差、根差和割曬損失率的試驗結果見表4。

表4 田間試驗結果

由表4可知，鋪放角的預測值與實測值的相對誤差小于6%，在參數(shù)優(yōu)化組合條件下，鋪放角由134.9°下降為115.8°，角度差由11.2°下降為9.6°，根差由0.22 m下降為0.14 m，損失率由2.2%下降為1.5%，說明優(yōu)化參數(shù)組合可使鋪放角更接近90°，同時降低了角度差、根差和損失率，有利于提高割曬作業(yè)質量。

試驗油菜品種的分枝較多，結穗位置較高，植株重心偏上，割曬作業(yè)中油菜頂端前傾較大，導致鋪放角偏大，但優(yōu)化參數(shù)條件下的鋪放角接近90°。實際割曬作業(yè)中，應根據植株高度和種植密度，調節(jié)割茬高度使割后植株高度接近1.4 m，種植密度較大時選用較小的前進速度和較大的速比，使鋪放角接近90°，以提高油菜割曬作業(yè)質量。

5 結 論

1）闡述了立式油菜割曬機結構和工作原理，將鋪放過程分為下落過程和轉動過程，分析油菜植株在2個過程中的運動狀態(tài)和受力情況，得出鋪放角的函數(shù)模型并分析影響因素。

2）利用ANSYS軟件和ProE軟件分別建立油菜植株的柔性體和割曬機模型，導入ADAMS中建立割曬機-油菜植株剛柔耦合模型，分別開展單因素試驗和多因素優(yōu)化試驗研究油菜植株參數(shù)和割曬機參數(shù)對鋪放角的影響規(guī)律。

3）單因素試驗結果表明，鋪放角隨油菜植株高度的增加先減小后增加，油菜植株質量密度增加時鋪放角先增加后減小，鋪放角隨油菜植株與植株的摩擦系數(shù)增加而減少。多因素試驗結果表明試驗因素對鋪放角影響顯著性順序大小為速比、前進速度、割茬高度；建立鋪放角的三元二次回歸模型，利用響應面法分析因素交互作用對鋪放角的影響規(guī)律。

4）建立鋪放角目標優(yōu)化方程，得到割曬機優(yōu)化參數(shù)組合：前進速度0.6 m/s，速比1.37，割茬高度0.42 m。田間試驗結果表明，鋪放角的預測值與實測值的相對誤差小于6%，在優(yōu)化參數(shù)組合條件下，鋪放角由134.9°下降為115.8°，角度差由11.2°下降為9.6°，根差由0.22 m下降為0.14 m，損失率由2.2%下降為1.5%，說明優(yōu)化參數(shù)組合可使鋪放角更接近90°，同時降低角度差、根差和損失率，提高了割曬作業(yè)質量。

由于仿真精度和時間限制，本文將油菜植株的相互牽連力添加為定值作用力，導致鋪放角預測值與實測值的相對誤差接近6%。實際作業(yè)中牽連作用力受油菜品種、重疊系數(shù)、含水率、摩擦系數(shù)等多個因素影響是變化的，后續(xù)需研究油菜植株與植株的相互作用力在鋪放過程中的變化及其對割曬作業(yè)質量的影響，進一步提高仿真結果的可靠性。

[1] 吳崇友，金誠謙，肖體瓊，等. 我國油菜全程機械化現(xiàn)狀與技術影響因素分析[J]. 農機化研究，2007，29(12)：207-210. Wu Chongyou, Jin Chengqian, Xiao Tiqiong, et al. Analysis on rape mechanization present situation and technical affection factors in whole productive course in China[J]. Journal of Agricultural Mechanization Research, 2007, 29(12): 207-210. (in Chinese with English abstract)

[2] 吳崇友，易中懿. 我國油菜全程機械化技術路線的選擇[J]. 中國農機化學報，2009，222(2)：3-6. Wu Chongyou, Yi Zhongyi. Selection of technology route of oilseed rape mechanization in entire production proceeding in China[J]. Journal of Chinese Agricultural Mechanization, 2009, 222(2): 3-6. (in Chinese with English abstract)

[3] 金誠謙，吳崇友，石磊. 油菜生產全程機械化技術體系關鍵技術研究[J]. 農機化研究，2010，32(5)：221-223. Jin Chengqian, Wu Chongyou, Shi Lei. Key technologies for the systems of rape mechanization in whole productive course[J]. Journal of Agricultural Mechanization Research, 2010, 32(5): 221-223. (in Chinese with English abstract)

[4] 羅海峰，湯楚宙，官春云，等. 油菜機械化收獲研究進展[J]. 農機化研究，2015，37(1)：1-8. Luo Haifeng, Tang Chuzhou, Guan Chunyun, et al. Research progress in rapeseed mechanized harvest[J]. Journal of Agricultural Mechanization Research, 2015, 37(1): 1-8. (in Chinese with English abstract)

[5] 吳崇友，王積軍，廖慶喜，等. 油菜生產現(xiàn)狀與問題分析[J]. 中國農機化學報，2017，38(1)：124-131. Wu Chongyou, Wang Jijun, Liao Qingxi, et al. Current status and problems of rapeseed production[J]. Journal of Chinese Agricultural Mechanization, 2017, 38(1): 124-131. (in Chinese with English abstract)

[6] 陳紅琳，陳尚洪，沈學善，等. 不同收獲方式對油菜子粒損失、含油量及種植效益的影響[J]. 作物雜志，2015，30(5)：74-79. Chen Honglin, Chen Shanghong, Shen Xueshan, et al. Effects of harvest methods on yield loss rate，oil content and planting efficiency for rape[J]. Crops, 2015, 30(5): 74-79. (in Chinese with English abstract)

[7] 石增祥，吳明亮，楊文敏，等. 我國油菜分段收獲機械研究現(xiàn)狀與發(fā)展對策[J]．農業(yè)工程，2015，5(5)：1-4． Shi Zengxiang, Wu Mingliang, Yang Wenmin, et al. Research status and development measures of rape segment harvester in China[J]．Agricultural Engineering, 2015, 5(5): 1-4．(in Chinese with English abstract)

[8] 吳崇友，丁為民，石磊，等．油菜分段收獲撿拾脫粒機撿拾損失響應面分析[J]．農業(yè)機械學報，2011，42(8)：89-93. Wu Chongyou, Ding Weimin, Shi Lei, et al. Response surface analysis of pickup losses in two-stage harvesting for rapeseed[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(8): 89-93. (in Chinese with English abstract)

[9] 吳崇友，肖圣元，金梅. 油菜聯(lián)合收獲與分段收獲效果比較[J]. 農業(yè)工程學報，2014，30(17)：10-16. Wu Chongyou, Xiao Shengyuan, Jin Mei. Comparation on rape combine harvesting and two-stage harvesting[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(17): 10-16. (in Chinese with English abstract)

[10] 金誠謙，尹文慶，吳崇友. 油菜割曬機撥指輸送鏈式輸送裝置研制與試驗[J]. 農業(yè)工程學報，2013，29(21)：11-18. Jin Chengqian, Yin Wenqing, Wu Chongyou. Development and experiment of rape windrower transportation device with poke finger conveyor chain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(21): 11-18. (in Chinese with English abstract)

[11] 李平，廖慶喜，李磊，等. 4SY-1.8改進型油菜割曬機主要裝置設計與試驗[J]. 農業(yè)機械學報，2014，45(1)：53-58. Li Ping, Liao Qingxi, Li Lei, et al. Design and experiment of the main device of 4SY-1.8 modified rape windrower[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(1): 53-58. (in Chinese with English abstract)

[12] 李平，廖慶喜，舒彩霞，等. 油菜割曬機莖稈鋪放質量的故障成因分析與參數(shù)匹配研究[J]. 應用基礎與工程科學學報，2016，24(1)：197-209. Li Ping, Liao Qingxi, Shu Caixia, et al. Fault analysis of stem windrowing quality and parameters match on the rape windrower[J]. Journal of Basic Science and Engineering, 2016, 24(1): 197-209. (in Chinese with English abstract)

[13] 金誠謙，吳崇友，金梅，等． 4SY-2 型油菜割曬機設計與試驗[J]. 農業(yè)機械學報，2010，41(10)：76-79． Jin Chengqian, Wu Chongyou, Jin Mei, et al. Design and experiment of 4SY-2 rape windrower[J]. Transactions of the Chinese Society for Agricultural Machinery, 2010, 41(10): 76-79．(in Chinese with English abstract)

[14] 金誠謙，尹文慶，吳崇友．4SY–2型油菜割曬機鋪放質量數(shù)學模型與影響因素分析[J]．農業(yè)工程學報，2012，28(2)：45-48． Jin Chengqian, Yin Wenqing, Wu Chongyou. Mathematical model and influencing factors analysis for windrow quality of 4SY-2 rape windrower[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(2): 45-48. (in Chinese with English abstract)

[15] 石增祥，謝偉，任述光，等. 4SY-2.2型油菜割曬機設計與試驗[J]. 中國農學通報，2017，33(10)：140-145. Shi Zengxiang, Xie Wei, Ren Shuguang, et al. Design and test of 4SY-2.2 rape windrower[J]. Chinese Agricultural Science Bulletin, 2017, 33(10): 140-145. (in Chinese with English abstract)

[16] 王修善，劉大為，李旭，等. 4SY–2. 0 型自走式油菜割曬機的設計與試驗[J]. 湖南農業(yè)大學學報：自然科學版，2016，42(4)：445-453. Wang Xiushan, Liu Dawei, Li Xu, et al. Design and experiment of 4SY-2.0 self-propelled rape windrower[J]. Journal of Hunan Agricultural University: Natural Sciences, 2016, 42(4): 445-453. (in Chinese with English abstract)

[17] 賀俊林，胡偉，郭玉富，等. 扶禾桿在不對行導入玉米莖稈中的運動仿真[J]. 農業(yè)工程學報，2007，23(6)：125-129. He Junlin, Hu Wei, Guo Yufu, et al. Kinematic simulation of no-row feed-in mechanism with guide-rod for corn harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(6): 125-129. (in Chinese with English abstract)

[18] 杜岳峰，毛恩榮，宋正河，等. 基于ADAMS的玉米植株收獲過程仿真[J]. 農業(yè)機械學報，2012，43(S)：106-111. Du Yuefeng, Mao Enrong, Song Zhenghe, et al. Simulation on corn plants in harvesting process based on ADAMS[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(S): 106-111. (in Chinese with English abstract).

[19] 羅海峰，湯楚宙，官春云，等. 適應機械化收獲的田間油菜植株特性研究[J]. 農業(yè)工程學報，2010，26(S1)：61-66. Luo Haifeng, Tang Chuzhou, Guan Chunyun, et al. Plant characteristic research on field rape based on mechanized harvesting adaptability[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(S1): 61-66. (in Chinese with English abstract)

[20] 董曉芳，田保明，姚永芳，等. 密度對油菜品種機械化收獲特性的影響[J]. 中國農學通報，2012，28(3)：71-74. Dong Xiaofang, Tian Baoming, Yao Yongfang, et al. Effects of the density on the characteristics of the mechanization harvest in brassica napus L.[J]. Chinese Agricultural Science Bulletin, 2012, 28(3): 71-74. (in Chinese with English abstract)

[21] 羅海峰，吳明亮，翁偉，等. 成熟期油菜莖稈物理模型試驗研究[J]. 湖南農業(yè)大學學報，2013，39(S1)：88-89，93. Luo Haifeng, Wu Mingliang, Weng Wei, et al. Experimental study on physical model of rape stem at mature stage[J]. Journal of Hunan Agricultural University, 2013, 39(S1): 88-89, 93. (in Chinese with English abstract)

[22] 金鑫. 4SY-1.8型油菜割曬機輸送系統(tǒng)的研制[D]. 武漢：華中農業(yè)大學，2012. Jin Xin. Design on Conveying System of 4SY-1.8 Rape Windrower[D]. Wuhan: Huazhong Agricultural University, 2012. (in Chinese with English abstract)

[23] 李平. 4SY-1.8型油菜割曬機關鍵部件工作機理與試驗研究[D]. 武漢：華中農業(yè)大學，2014. Li Ping. Experimental Study and Working Mechanism of Key Parts for 4SY-1.8 Rape Windrower[D]. Wuhan: Huazhong Agricultural University, 2014. (in Chinese with English abstract)

[24] 吳明亮，官春云，湯楚宙，等. 油菜莖稈切割力影響因素試驗[J]. 農業(yè)工程學報，2009，25(6)：141-144. Wu Mingliang, Guan Chunyun, Tang Chuzhou, et al. Experiments on influencing factors of cutting force of rape stem[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(6): 141-144. (in Chinese with English abstract).

[25] 伍文杰，吳崇友，江濤. 油菜莖稈切割試驗研究驗[J]. 農機化研究，2018，40(6)：141-144. Wu Wenjie, Wu Chongyou, Jiang Tao. Experimental study on cutting of rape salks[J]. Journal of Agricultural Mechanization Research, 2018, 40(6): 141-144. (in Chinese with English abstract).

[26] 崔濤，劉佳，張東興，等. 基于ANSYS和ADAMS的玉米莖稈柔性體仿真[J]. 農業(yè)機械學報，2012，43(S)：112-115. Cui Tao, Liu Jia, Zhang Dongxing, et al. Flexible body simulation for corn stem based on ANSYS and ADAMS[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(S): 112-115. (in Chinese with English abstract).

[27] 宗望遠，廖慶喜，黃鵬，等. 組合式油菜脫粒裝置設計與物料運動軌跡分析[J]. 農業(yè)機械學報，2013，44(S2)：41-46. Zong Wangyuan, Liao Qingxi, Huang Peng, et al. Design of combined rape threshing device and analysis of rape cane movement trail[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(S2): 41-46. (in Chinese with English abstract).

[28] 冉軍輝，沐森林，李海同，等. 油菜聯(lián)合收獲機往復式雙動割刀行星輪驅動器設計與試驗[J]. 農業(yè)工程學報，2020，36(9)：17-25. Ran Junhui, Mu Senlin, Li Haitong, et al. Design and test of planet gear driver of reciprocating double-acting cutter for rapeseed combine harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(9): 17-25. (in Chinese with English abstract)

[29] 施麗莉，胡志超，顧峰瑋，等. 耙齒式壟作花生殘膜回收機設計及參數(shù)優(yōu)化[J]. 農業(yè)工程學報，2017，33(2)：8-15. Shi Lili, Hu Zhichao, Gu Fengwei, et al. Design and parameter optimization on teeth residue plastic film collector of ridged peanut[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(2): 8-15. (in Chinese with English abstract)

[30] 謝建華，唐煒，張學軍，等. 擺桿驅動式殘膜回收機的設計與參數(shù)優(yōu)化[J]. 農業(yè)工程學報，2019，35(4)：56-63. Xie Jianhua, Tang Wei, Zhang Xuejun, et al. Design and parameter optimization of swing rod driven residual plastic film collector[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(4): 56-63. (in Chinese with English abstract)

[31] 楊風波，薛新宇，蔡晨，等. 航空專用離心噴頭霧化性能試驗與影響因子研究[J]. 農業(yè)機械學報，2019，50(9): 96-104． Yang Fengbo, Xue Xinyu, Cai Chen, et al. Atomization performance test and influence factors of aviation special centrifugal nozzle[J]. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(9): 96-104．(in Chinese with English abstract)

[32] 張佳喜，汪珽玨，陳明江，等. 齒盤式棉稈收獲機的設計[J]. 農業(yè)工程學報，2019，35(15)：1-8. Zhang Jiaxi, Wang Tingjue, Chen Mingjiang, et al. Design of toothed disc cotton stalk harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(15): 1-8. (in Chinese with English abstract)

[33] 沈公威，王公仆，胡良龍，等. 甘薯莖尖收獲機研制[J]. 農業(yè)工程學報，2019，35(19)：46-55. Shen Gongwei, Wang Gongpu, Hu Lianglong, et al. Development of harvesting mechanism for stem tips of sweet potatoes[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(19): 46-55. (in Chinese with English abstract)

Formation mechanism of laying angle of vertical rape windrower based on ANSYS-ADAMS

Li Haitong, Wu Chongyou※, Mu Senlin, Guan Zhuohuai, Jiang Tao

(,,210014,)

In order to improve the operation quality of vertical rape windrower, the effects of rape plant characteristics and windrower parameters on laying angle were studied in this study. The structure and working principle of the vertical rape mower were described. The laying process of rape plant was divided into falling process and rotation process according to the working principle of the rape windrower. The movement state of rape plant during the laying process was analyzed based on the dynamic theory, and the mathematical model of laying angle was established and the influencing factors were analyzed. Analysis results showed that the laying angle influenced by material parameters of rape plant and working parameters of windrower. ANSYS and ProE software were used to establish the flexible body of rape plant and windrower model, respectively. Then the rigid flexible coupling model of windrower and rape plant was established in ADAMS software, and the single factor and multifactor simulation experiments of laying process were carried out. The results of single factor experiments showed that with the increasing of the plant height in the range of 1.2-1.6 m, laying angle decreased and reached the minimum value at the plant height of 1.4m, and then increased. The laying angle increased and then decreased when the plant mass density was in the range of 438-588 kg/m3, and reached the maximum value at the plant mass density of 538 kg/m3. When the friction coefficient between plants increased in the range of 0.1-0.5, the laying angle decreased. The results of multifactor experiments showed that forward speed, speed ratio and the interaction between forward speed and stubble height had a significant influence on laying angle (< 0.01).The significant order of factors was speed ratio, forward speed and stubble height. The quadratic regression equation model was established according to the results of multifactor experiments, and the response surface methodology was utilized to analyze the influence of interaction on the laying angle. Optimization equations and constraints were established, and optimal parameters of windrower were determined that forward speed was 0.6 m/s, the speed ratio was 1.37, and stubble height was 0.42 m. The field experiments were carried out and the results indicated that the relative errors between predicted and measured values of laying angle were less than 6%。Compared with the control tests, the laying angle of rape plant decreased from 134.9° to 115.8°, the angle difference decreased from 11.2°to 9.6°, the root difference decreased from 0.22 m to 0.14 m, and the loss rate decreased from 2.2% to 1.5%, the working performance of rape windrower improved. The research results can provide the basis for determining the suitable conditions of rape cutting and optimizing the operation parameters of the windrower.

agricultural machinery; harvesting; rape; laying angle; rigid-flexible coupling model; response surface; optimization

李海同，吳崇友，沐森林，等. 基于ANSYS-ADAMS的立式油菜割曬機鋪放角形成機理 [J]. 農業(yè)工程學報，2020，36(14)：96-105.doi：10.11975/j.issn.1002-6819.2020.14.012 http://www.tcsae.org

Li Haitong, Wu Chongyou, Mu Senlin, et al. Formation mechanism of laying angle of vertical rape windrower based on ANSYS-ADAMS[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 96-105. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.14.012 http://www.tcsae.org

2020-04-20

2020-06-24

國家重點研發(fā)計劃項目（2016YFD0702101）；國家油菜產業(yè)體系專項資助項目（CARS-13-10B）

李海同，博士，助理研究員，主要從事收獲機械研究。Email：lhtdolphin@163.com

吳崇友，博士，研究員，博士生導師，主要從事收獲機械研究研發(fā)。Email：542681935@qq.com

10.11975/j.issn.1002-6819.2020.14.012

S225.2

A

1002-6819(2020)-14-0096-10