小麥收割機虛擬樣機運動仿真研究
———基于離散元虛擬現實技術

范媛媛，米西峰

(焦作師范高等專科學校 計算機與信息工程學院，河南 焦作 454000)

小麥收割機虛擬樣機運動仿真研究
———基于離散元虛擬現實技術

范媛媛，米西峰

(焦作師范高等專科學校 計算機與信息工程學院，河南 焦作 454000)

為了克服傳統小麥收割機設計過程中反復試驗-制造的長周期與低效率的缺點，且達到降低機械設計成本的目的，提出了一種基于離散元虛擬技術的小麥收割機虛擬樣機設計方法，實現了小麥收割機的虛擬優化設計。采用Windows平臺下VS2010集成了離散元仿真軟件，并使用OpenGL軟件仿真實現了脫粒過程的三維顯示，利用并行計算的方法對數據進行處理，最后通過保存三維輸出結果，對脫粒的性能進行了數據統計和計算。虛擬仿真測試表明：所提出的小麥植株建模方法和收割機虛擬樣機的性能仿真是可行的，通過顆粒捕捉仿真得到了麥粒的運動軌跡，并根據對脫粒的性能進行了統計，最終得到了小麥收割機脫粒的損失率、破碎率和脫凈率，對現代化農機的設計制造具有重要的現實意義。

小麥收割機；虛擬現實；并行計算；OpenGL軟件；離散元

0 引言

隨著計算機仿真技術的不斷發展，計算機的性能得到了大幅度提升，在各個領域得到了廣泛的應用。在小麥收割機的設計過程中，如果通過虛擬技術來仿真小麥的脫粒過程，可以克服傳統機械設計過程中反復試驗-制造的長周期與低效率缺點，從而達到降低機械設計成本的目的。利用計算機離散元虛擬仿真技術，發現小麥收割機細節設計中存在的問題，然后對裝置的結構進行優化設計，通過對不同結構方案的比較，找出存在的問題，以提高小麥收割機的設計質量。這種方法既省時又省力，相對于傳統的設計方法，具有無可比擬的優勢，可以為現代化農機設計制造提供一種嶄新的方法。

1 小麥收割機虛擬樣機設計

在小麥收割機的傳統設計過程中，利用傳統物理樣機設計時，只有將真實的小麥收割機樣機制造出來，才能夠發現設計中存在的問題和缺陷。但是，樣機的試制需要較長的時間，導致設計成本也較高。傳統設計流程如圖1所示。

圖1 小麥收割機傳統設計流程

在樣機的試驗階段，如果發現問題，還需要對樣機進行重新設計和加工制造，會使周期和成本進一步加大，浪費了大量的人力和物力。

圖2為采用虛擬樣機設計方法對小麥收割機進行優化設計。該方法利用計算機虛擬仿真技術，發現小麥收割機細節設計中存在的問題，然后對裝置的結構進行優化設計；通過對不同結構方案的比較，找出存在的問題，以提高小麥收割機的設計質量，既省時又省力，相對于傳統設計方法具有無可比擬的優勢。

圖3為小麥收割機割臺的虛擬樣機設計示意圖。由圖3可以看出：采用虛擬樣機可以設計出和真實的試制樣機非常相似的模型。本次設計主要是在Windows平臺下的VS2010中進行的，該軟件采用面向對象的設計方法，可以實現界面的可視化和預覽功能。采用編程的方法將小麥脫粒仿真軟件集成到離散元軟件中，其總體框架如圖4所示。

圖2 小麥收割機虛擬樣機設計流程

圖3 小麥收割機割臺虛擬樣機

圖4 小麥收割機虛擬樣機總體設計框架

小麥收割機虛擬樣機設計的總體框架主要分為3部分，包括離散建模、仿真計算和結果處理等。使用Open GL軟件仿真實現了脫粒過程的三維顯示，最后通過保存三維輸出結果，對脫粒的性能進行數據統計和計算。

2 小麥離散元建模和力學虛擬仿真

利用小麥自身的結構特征，構建小麥的三維離散元模型。由小麥的植株的實體特征發現：在每個秸稈的鄰節之間，有一個凸出的部分將多個小麥秸稈節連接在一起，每個節的尺寸不相同。為了虛擬仿真方便，將小麥節間進行編號，根據當地小麥的尺寸特征，在生成小麥節間時，結合節間的長度L、節間球的最大和最小尺寸Rmax和Rmin、相鄰節間球的重合度C，以及節間的密度ρ進行離散元建模。設節間球的數量為n，則可以得到n的近似量為

(1)

利用該方法建立的離散元模型，其兩個相鄰接的小球半徑R1、R2的差值公式為

R1-R2=a-C

(2)

其中，α為一個常數，該值表示相鄰的鄰接球的半徑的差值。節間或相鄰節間球也是符合上述規律的，在進行質量求解時，可以將節間看成是圓柱體來進行處理，于是節間的質量公式為

(3)

(4)

其中，P為莖稈在外力的作用下產生的沿著該方向上的變形；YD為虛擬仿真體到地面的垂直距離；b為P到地面的距離；E為小麥秸稈的彈性模量；d為秸稈的內徑；l為小麥的秸稈在地面以上的長度距離。將小麥的植株看作是空心的圓柱體，則

(5)

建立完成小麥植株的離散元模型后，需要通過力學分析來仿真模擬小麥收割機的運動性能，其流程如圖5所示。

圖5 離散元力學并行虛擬計算流程

考慮到在實際仿真過程中小麥植株的顆數較多，因此采用并行計算的方法，將小麥植株離散成多個部分，將力學分析賦予給每個離散單元體。通過多線程的方法計算線程的受力，并將多線程的受力賦給小麥植株，每次計算后將臨時保存的力數據進行清空，從而完成多植株的離散元虛擬仿真。

3 小麥收割機虛擬樣機運動仿真

小麥收割機虛擬樣機的運動仿真平臺可以使用Windows平臺下的VS2010軟件，采用面向對象的方法進行開發，從而實現小麥收割機虛擬樣機脫粒仿真和脫粒性能的統計。為了提高離散元虛擬仿真的計算速度，可以采用OpenMP把串行程序改成并行程序，由32位向64位進行拓展，最后通過編程的方式將小麥脫粒的仿真程序集成到離散元仿真軟件中。

圖6為本次小麥收割機樣機運動虛擬仿真建立的離散元模型。通過虛擬仿真可以得到模型的三維動畫演示，將計算輸出的結果文件進行保存，最后利用數據統計軟件對其性能進行分析。

根據三維仿真得到的數據結果可對單個的離散單元顆粒進行運動軌跡的捕捉。圖7為通過顆粒捕捉仿真得到的麥粒運動軌跡，由此軌跡可以實現麥粒最終脫落位置的標記，根據脫落位置可以統計脫粒得到的個數及是否破碎，最終得到小麥收割機脫粒的損失率，進而驗證小麥收割機的性能指標。

圖6 離散元虛擬仿真模型

圖7 麥粒運動軌跡捕捉

小麥收割機性能指標的分析主要包括小麥植株的顆粒損失率、顆粒破碎率和小麥植株的脫凈率，下面分別通過公式來說明性能分析的指標。

1)脫凈率。小麥收割機虛擬樣機的脫凈率表示在虛擬仿真過程中小麥顆粒脫落的數目和仿真小麥顆粒的總數目比值。脫凈率越高，表示小麥脫粒機的性能和效果越好。其計算公式為

(6)

2)破碎率。小麥收割機虛擬樣機的破碎率為虛擬仿真過程中小麥顆粒破碎的數目和仿真小麥顆粒的總數目比值。破碎率越高，表示小麥脫粒機的性能和效果越差。其計算公式為

(7)

3)籽粒損失率。籽粒損失率為小麥顆粒總數去掉脫落小麥的顆粒個數最后與總的仿真小麥顆粒數目的比值。其計算公式為

(8)

根據這3項指標的定量，對離散元仿真的三維模擬結果進行數據統計，然后對虛擬樣機進行優化設計，最后再次統計性能數據，得到了如表1所示的性能測試結果。

表1 模型優化前后虛擬樣機的性能數據 %

由表1可以看出：優化后的虛擬樣機性能數據要優于優化之前。這表明，采用離散單元方法對小麥收割機虛擬樣機進行優化設計的方法可行。

4 結論

針對小麥收割機設計過程中存在的設計周期長和成本高的缺點，提出了一種基于離散元虛擬仿真技術的小麥收割機設計優化方法，并建立了小麥植株的離散元虛擬仿真模型，最后對小麥收割機的性能進行了虛擬仿真計算。通過虛擬仿真計算得到了模型的三維動畫演示，將計算輸出的結果文件進行保存，利用數據統計軟件對其性能進行了分析，捕捉到了單個麥粒的運動軌跡，并由大量麥粒的運動軌跡的統計得到了小麥收割機的損失率、破碎率和脫凈率等性能指標。最后，對優化前后的小麥收割機的性能指標進行了對比，結果表明：采用離散單元方法對小麥收割機虛擬樣機進行優化設計的方法可行，能夠有效地提高收割機的脫凈率指標，從而降低損失率和破碎率，提高小麥收割機的設計效率。

[1] 雷曉俊,湯亮,張永會.小麥麥穗幾何模型構建與可視化[J].農業工程學報,2011,27(3):179-184.

[2] 于亞軍,周海玲,付宏,等.基于顆粒聚合體的玉米果穗建模方法[J].農業工程學報,2012,28(8):167-174.

[3] 劉凌霞,宋強.基于簡化的質點-彈簧模型織物變形仿真研究[J].計算機仿真,2011,28(7):406-409.

[4] 許紅,許浩然.彈簧-質點模型中彈簧力計算問題[J].湖北工業大學學報,2011,26(4):82-83.

[5] 青克樂其其格,郭新宇.基于脈序骨架的小麥葉精細建模研究[J].農機化研究,2012,34(5):53-57.

[6] 孔曉莉,蔣宇翔,劉樹東,等.基于Bezier曲線的卷曲水草模擬[J].蘇州大學學報,2011,31(6):1-5.

[7] 郝小忠,AHMEDA.A.D,陳文亮.Bezier曲面建模方法研究與實現[J].中國制造業信息化,2011,40(23): 42-45.

[8] 程飛.基于OpenGL的雙三次Bezier曲面的繪制[J].安徽電子信息職業技術學院學報,2010,46(9): 33-34.

[9] 王振華,王德成,劉貴林,等.方草捆壓捆機撿拾器參數設計[J].農業機械學報,2010,41(S1):107-109.

[10] 王建政.4GL-1.80型多功能小雜糧收割機的研制[J].中國農業大學學報,2009,14(2):103-106.

[11] 鄭浩,馮毅雄,譚建榮,等.一類制造資源的協同建模、優化與求解技術[J].計算機集成制造系統，2012(7):1387-1395.

[12] 尹超,黃必清,劉飛,等.中小企業云制造服務平臺共性關鍵技術體系[J]. 計算機集成制造系統,2011,17(3): 495-503.

[13] 戰德臣,趙曦濱,王順強,等.面向制造及管理的集團企業云制造服務平臺[J].計算機集成制造系統，2011,17(3): 487-494.

[14] 馬志欣,譚峰,侯召龍.基于 Android 的農業生產遠程監控系統[J].農機化研究,2014,36(10):226-229.

[15] 高偉,康倩.淺談嵌入式系統[J].計算機光盤軟件與應用,2014,14(8):69-70.

[16] 王崴,張宇紅,徐曉東，等.面向中小企業的零部件資源云端化封裝技術研究[J].機械設計與制造,2012(8):260-262.

[17] 賀瑤, 王文慶, 薛飛.基于云計算的海量數據挖掘研究[J].計算機技術與發展, 2013, 23(2): 69-72.

[18] 丁巖, 楊慶平, 錢煜明. 基于云計算的數據挖掘平臺架構及其關鍵技術研究[J]. 中興通訊技術, 2013,19(1): 53-60.

[19] 錢建平,楊信廷,吳曉明,等.自然場景下基于混合顏色空間的成熟期蘋果識別方法[J].農業工程學報,2012,28(17):137-142.

[20] 顧寶興,姬長英,王海青,等.智能移動水果采摘機器人設計與試驗[J].農業機械學報,2012,43(6):153- 160.

[21] 王海青,姬長英,顧寶興,等.基于機器視覺和支持向量機的溫室黃瓜識別[J].農業機械學報,2012,43(3): 163-167.

[22] 王輝,毛文華,劉剛,等.基于視覺組合的蘋果作業機器人識別與定位[J].農業機械學報,2012,43(12):165-170.

[23] 魏澤鼎,賈俊國,王占永.基于視覺傳感器的棉花果實定位方法[J].農機化研究,2012,34(6):66-68,112.

[24] 李立君,李昕,高自成,等.基于偏好免疫網絡的油茶果采摘機器人圖像識別算法[J].農業機械學報,2012,43(11):209-213.

[25] 王輝,毛文華,劉剛,等.基于視覺組合的蘋果作業機器人識別與定位[J].農業機械學報,2012,43(12): 165-170.

[26] 袁挺,紀超,陳英,等.基于光譜成像技術的溫室黃瓜識別方法[J].農業機械學報,2011,42(S1):172-176.

[27] 熊俊濤,鄒湘軍,陳麗娟,等.采摘機械手對擾動荔枝的視覺定位[J].農業工程學報,2012,28(14):36-41.

[28] 呂繼東,趙德安,姬偉,等.采摘機器人振蕩果實動態識別[J].農業機械學報,2012,43(5):173-196.

Motion Simulation of Wheat Harvester Virtual Prototype—Based on Reality Technology of Discrete Element Virtual

Fan Yuanyuan, Mi Xifeng

(School of Computer and Information Engineering, Jiaozuo Teachers College ,Jiaozuo 454000, China)

In order to overcome the disadvantages of long period and low efficiency shortcomings in traditional wheat harvester design process of repeated experiments - manufacture, save the cost of the mechanical design, it put forward a wheat harvester virtual prototype design of discrete element method based on virtual technology, the virtual design and optimization of wheat harvester. The use of Windows platform VS2010 integrated discrete element simulation software, and use the Open GL software simulation to achieve a three-dimensional display of threshing process by using parallel computing method, data processing, and finally through the preservation of the three-dimensional output, the threshing performance of data statistics and calculation. Through the simulation test shows that the performance simulation of the proposed modeling method and wheat harvester virtual prototype is feasible, the particle capture are simulated trajectories of wheat, and the track performance of threshing were collected, finally got the wheat threshing loss rate, broken rate and off the net rate of design of modern agricultural machinery manufacturing,which has important practical significance.

wheat harvester; virtual reality; parallel computing; Open GL software; discrete element

2016-12-11

河南省自然科學基金項目(2015GZC133)

范媛媛(1976-)，女，河南焦作人，副教授，碩士。

米西峰(1976-)，男，河南南陽人，副教授，碩士，(E-mail)mixifeng@163.com。

S225.3；TS112.2

A

1003-188X(2018)02-0024-05