劉 基，金誠謙,2，梁蘇寧，倪有亮，王 昕，李澤峰

(1.農(nóng)業(yè)部南京農(nóng)業(yè)機械化研究所，南京 210014；2.山東理工大學，山東 淄博 255000)

0 引言

機械化收獲在大豆生產(chǎn)過程中起著重要作用，目前黃淮海地區(qū)的大豆收獲大多使用稻麥聯(lián)合收割機[1]，大豆專用收割機較少，且收獲效果并不理想，損失率和破碎率都比較高。改進大豆收獲機械的工作性能可以有效提高收獲質(zhì)量，減少收獲環(huán)節(jié)的損失，推進大豆產(chǎn)業(yè)更好的發(fā)展[2]。大豆收割機的作用對象是大豆莖稈、大豆籽粒及豆莢，從作用對象的物理特性入手，掌握物料的各項力學參數(shù)[3-4]，可以更加準確地對收割機的機械部件進行改進。在收割機的脫粒滾筒中，莖稈與大豆籽粒、豆莢以及釘齒和凹版篩之間有相互作用力[5]，掌握莖稈的物理特性[6]可以更加準確地對收割機的脫粒系統(tǒng)進行研究[7-8]。

國內(nèi)對作物莖稈的力學特性進行了一些研究。閆以勛[9]等人從免耕播種機設計的角度對成熟期的大豆莖稈進行了試驗研究，發(fā)現(xiàn)底部節(jié)間能承受的彎曲力較大。吳杰、王艷云[10-11]等人對棉稈的壓縮特性進行了試驗研究。何曉莉[12]等人從大豆莖稈倒伏的角度出發(fā)，通過軸向壓縮試驗對干濕大豆莖稈做了對比研究。李紅波[13]等人從抗倒伏的角度出發(fā)，通過折彎試驗對小麥莖稈進行了研究，發(fā)現(xiàn)生長期會引起抗倒伏特性的變化。袁志華[14]等人建立了適用于小麥、水稻等作物莖稈的與倒伏相關的力學模型。劉慶庭[15]等人認為，作物莖稈力學特性的研究應該考慮生物材料的各向異性及復雜結構。王芬娥[16]等人對小麥莖稈的微觀結構進行了研究。周揚[17]等人對大麻莖稈的力學特性進行研究，獲得了大麻莖稈的剪切模量等參數(shù)。目前，還沒有從優(yōu)化脫粒裝置的角度對大豆莖稈進行的多品種對比研究。

考慮到莖稈是作為一個整體與滾筒部件及大豆籽粒產(chǎn)生相互作用力，本文對大豆莖稈的整體結構進行分析。先把莖稈按照圓柱體[18]模型進行力學試驗，再結合莖稈的內(nèi)部結構對試驗結果進行分析。由于黃淮海地區(qū)的大豆品種多樣，各個品種的大豆莖稈物理特性不盡相同，為了使研制的大豆收割機通用性較強，有必要對多個品種的莖稈做對比分析[19]，找到共性及各自的特點，為之后大豆收割機上脫粒裝置的研究提供理論依據(jù)。

1 試驗材料與方法

1.1 試驗材料

試樣選用成熟期安徽柳豐種業(yè)良繁基地的5個大豆品種(見圖1)，每個品種隨機取30株，試樣均貼地表面截取，平均含水率65%。阜豆15的平均株高為764mm，荷豆19的平均株高為613mm，柳豆109的平均株高為782mm，皖豆28的平均株高為687mm，皖豆33的平均株高為684mm。

1.2 試驗設備

UTM6503電子萬能試驗機(見圖2)、測試夾具、DM.3數(shù)字電子稱、滬制01130048游標卡尺、烘干箱、直尺、卷尺及鐮刀等。

圖1 試驗材料

圖2 試驗設備

1.3 試驗方法

5個大豆品種分別取10株進行三點彎曲試驗、10株進行剪切試驗，量取每株莖稈的高度。每株莖稈去掉頂端最細一節(jié)(直徑過小不具代表性)，其余部分均勻分成6段，每段量取最小及最大直徑，求得其平均值作為該段的直徑。做彎曲試驗時跨距為40mm，加載速度為10mm/min。

2 三點彎曲試驗

2.1 彎曲試驗測定原理

利用三點彎曲夾具在萬能材料試驗機上進行試驗，跨距40mm。試驗前輸入直徑數(shù)據(jù)，試驗結束后慣性矩I、彈性模量E都可以直接從電腦中讀出，彎曲載荷隨位移的變化曲線在萬能材料試驗機配套軟件中自動生成，抗彎剛度由彈性模量E和慣性矩I的乘積計算得出。

2.2 彎曲試驗過程

試驗過程中發(fā)現(xiàn)，莖稈斷裂過程為(見圖3)：在半徑為5mm的刀刃的壓力下，被壓點逐漸被壓縮，載荷逐漸增大；之后接觸點周圍發(fā)生縱向裂開，達到最大載荷；隨著莖稈被壓彎，其他部位也發(fā)生縱向裂開現(xiàn)象，載荷力逐漸減小，直到下表面橫向斷裂；此時整個莖稈瞬間斷開，載荷急劇下降。

圖3 三點彎曲試驗莖稈斷裂過程

2.3 彎曲試驗數(shù)據(jù)分析

彎曲試驗數(shù)據(jù)如表1～表3所示。

表1 莖稈不同高度處的最大彎曲載荷

Table 1 The maximum bending load of different heights of the stem N

試樣編號彎曲力F阜豆15荷豆19柳豆109皖豆28皖豆331228.60198.69125.52189.60225.752105.82128.0695.1674.11125.61389.7795.0979.1457.8468.88475.2459.1066.7955.5269.51547.3738.2135.7639.0339.67613.7126.7513.7428.4514.68

表2 莖稈不同高度處的慣性矩

Table 2 The moment of inertia of different heights of the stem mm4

試樣編號慣性矩I阜豆15荷豆19柳豆109皖豆28皖豆331124.74140.93221.93136.37165.502104.9588.71155.3169.77112.56372.5360.2987.6269.7785.49467.9751.0967.9742.3653.29533.2126.4928.3039.3226.49613.8612.828.2414.149.40

表3 莖稈不同高度處的彈性模量

Table 3 The elastic modulus of different heights of the stem N/mm2

試樣編號彈性模量E阜豆15荷豆19柳豆109皖豆28皖豆331171.00402.12288.84306.52334.632307.00476.79421.75476.99354.853477.00616.38545.91525.05384.724394.00660.87775.06553.10607.325472.00837.911052.36598.32892.016458.001146.661366.57878.06896.36

2.3.1 彎曲載荷

圖4中(a)、(b)、(c)、(d)、(e)各自顯示的6條曲線分別表示每個品種的莖稈從底部到頂部的6段試樣受彎曲載荷的情況。觀察這5個圖，可發(fā)現(xiàn)如下共同點：莖稈從底部到頂部所承受最大彎曲載荷逐漸降低；莖稈底部的承載能力明顯大于其他部位；莖稈頂部的承載能力明顯小于其他部位。不同點是：荷豆19、柳豆109、皖豆33的莖稈承受最大彎曲載荷從底部到頂部均勻減小，阜豆15和皖豆28的莖稈承受最大彎曲載荷從底部到第2段迅速減小，從第2段到頂部緩慢減小。

通過分析，最大彎曲載荷從莖稈底部到頂部逐漸減小的原因如下：一是莖稈底部直徑最大，且隨著高度的增加直徑逐漸減??；二是從頂部到底部，莖稈的木質(zhì)化程度越來越高。5個品種變化規(guī)律的不同可能是由直徑大小的和木質(zhì)部生長的不規(guī)則變化引起的。

由表1可知：除柳豆109的最大載荷是125.52N外，阜豆15、荷豆19、皖豆28、皖豆33最大載荷都在189～228N之間，最小載荷在13～29N之間。

2.3.2 慣性矩

試驗所得各種類莖稈不同高度的慣性矩如表2所示，變化趨勢如圖5所示。

圖5 慣性矩隨莖稈高度的變化曲線

由圖5可知：試驗中5個品種的大豆莖稈的慣性矩都是從底部到頂部逐漸減小。

2.3.3 彈性模量

彈性模量是衡量物體抵抗彈性變形能力大小的尺度。圖6為5個品種的大豆莖稈的彈性模量隨莖稈高度變化的趨勢。

圖6 彈性模量隨莖稈高度的變化曲線

由圖6可知：從莖稈底部到頂部彈性模量逐漸增大。這一結論與閆以勛[9]等人提出的從莖稈底部到頂部彈性模量先減小、后增大的結論不一致，原因可能是測量方法不同。本試驗是從整體考慮，樣式使用莖稈原有的形態(tài)，未對莖稈進行切片[18]。即本文中的彈性模量是莖稈整體的彈性模量，是受莖稈結構的變化影響的，而不是傳統(tǒng)意義上材料的彈性模量。分析莖稈的原始結構可知：莖稈底部的木質(zhì)層較厚，并且隨莖稈高度的增加木質(zhì)化程度呈降低的趨勢，因此從底部到頂部管狀結構越來越明顯，抵抗變形的能力越來越強[20]。導致這一現(xiàn)象發(fā)生的其他原因有待進一步研究。

2.3.4 抗彎剛度

抗彎剛度是指物體抵抗其彎曲變形的能力，用彈性模量E和慣性矩I的乘積EI表示。抗彎剛度隨莖稈高度的變化曲線如圖7所示。

圖7 抗彎剛度隨莖稈高度的變化曲線

由圖7可知：除阜豆15莖稈的抗彎剛度是從底部到頂部先增大、后減小之外，荷豆19、柳豆109、皖豆28、皖豆33等4個品種的抗彎剛度都是從莖稈的底部到頂部呈逐漸減小的趨勢，且柳豆109的抗彎剛度最大，阜豆15的抗彎剛度最小。

3 剪切試驗

3.1 測定原理

試驗利用剪切夾具在萬能材料試驗機上進行。直徑在試驗前測得并輸入到與試驗機配套的軟件中，剪切力隨位移變化的曲線隨著試驗的進行在軟件中自動生成，最大剪切力在試驗結束時自動生成，切應力由剪切力和截面積計算求得。

3.2 試驗過程

剪切試驗(見圖8)過程中發(fā)現(xiàn)：剪切第1段試樣(莖稈最底端)時，刀刃逐漸切入莖稈中，當切到半徑處時，剪切力達到最大值；隨后割刀繼續(xù)向下深入，并在很短的時間內(nèi)斷裂，整個過程中莖稈沒有明顯的被壓縮現(xiàn)象。剪切第2段試樣時，刀刃接觸到莖稈后切入莖稈中，在刀刃向下深入的同時，伴隨著莖稈的徑向壓縮和縱向開裂，剪切力產(chǎn)生很大波動；當?shù)度猩钊氲桨霃轿恢脮r達到最大剪切力，隨后剪切力以比之前較快的速度減小，直到莖稈被剪斷，剪切力瞬間降低。剪切第3段及更高部位的莖稈時，刀刃接觸莖稈后，首先對莖稈徑向壓縮，壓縮過程伴隨著莖稈的縱向開裂，直到把莖稈壓縮到接近片狀時，刀刃開始深入到莖稈中對莖稈進行剪切，剪切力會有明顯增大；此時的莖稈相當于雙層的板材結構，當上層被切斷時剪切力達到最大值，隨后切斷下層，剪切力瞬間減小。

通過對莖稈的斷面觀察可知：莖稈由表皮、木質(zhì)部、中心髓部組成，且髓部呈泡沫狀，質(zhì)地柔軟，分析受力時可忽略。從莖稈第2段開始一直到頂部，莖稈的木質(zhì)化程度較底部有明顯降低，木質(zhì)部厚度逐漸減小，莖稈管狀化結構越來越明顯。當?shù)度袑ηo稈徑向施力時，莖稈被徑向壓縮，壓縮到一定程度時會出現(xiàn)莖稈的縱向開裂；當莖稈被壓縮到接近板材形狀時，刀刃才開始切入到莖稈中。

圖8 剪切試驗過程

3.3 剪切試驗數(shù)據(jù)分析

剪切試驗數(shù)據(jù)如表4所示。

表4 莖稈不同高度處的直徑、最大剪切力

3.3.1 剪切載荷

觀察5個品種的大豆莖稈特性曲線可知，其所呈現(xiàn)的規(guī)律大致相同，因此做統(tǒng)一分析。

圖9中每個圖中的6條曲線分別為莖稈底部到頂部的6段試樣所受剪切力隨位移的變化情況。由圖9可知:莖稈底部所需的剪切力遠大于其他部位，且隨著莖稈的高度增加，所需剪切力逐漸減小。其原因是：一方面大豆莖稈底部較為粗壯，直徑比較大；另一方面莖稈底部木質(zhì)化程度較大，莖稈的強度較大，在剪切過程中突然斷裂。從第2條曲線開始，剪切力隨刀刃位移的波動較大，開始階段剪切力均勻增大，之后開始有較大波動，且增大的速率比之前小；當達到一個峰值后開始緩慢減小，減小到一定程度后瞬間降為零。這與在試驗過程中觀察到的刀刃與莖稈之間的作用力情況一致。起初階段刀刃對莖稈進行壓縮，之后出現(xiàn)莖稈的縱向開裂；當壓縮到一定程度時刀刃切入莖稈中，把莖稈上層切斷后達到最大值，刀刃繼續(xù)向下，把下層切斷時剪切力瞬間降為零。

圖9 剪切力-位移曲線

3.3.2 直徑

圖10中的5條曲線分別為5個品種的大豆莖稈直徑隨莖稈高度的變化趨勢。由圖10可知：從莖稈底部到頂部直徑整體呈下降的趨勢，但并不是嚴格的下降，中部會有微小浮動。這說明，大豆莖稈不是規(guī)則的錐形圓柱體，也符合作物隨機生長的特性，但這種浮動對莖稈整體力學性能趨勢的研究影響不大。

圖10 直徑隨莖稈高度的變化曲線

3.3.3 切應力

物體由于外因(受力、濕度變化等)而變形時，在物體內(nèi)各部分之間產(chǎn)生相互作用的內(nèi)力，以抵抗這種外因的作用，并力圖使物體從變形后的位置恢復到變形前的位置。在所考察的截面某一點單位面積上的內(nèi)力稱為應力，同截面相切的稱為剪應力或切應力。切應力為

(1)

圖11中分別為5個品種的大豆莖稈切應力隨莖稈高度的變化趨勢。由圖11可知：隨莖稈高度的增加，切應力呈逐漸減小的趨勢，從底部到頂部變化較大。柳豆109從底部到第5段變化較為緩慢，最后一段迅速減?。黄渌?個品種從底部到第3段迅速減小，第3段到最后一段減小較為緩慢。

圖11 切應力隨莖稈高度的變化曲線

4 結論

1)通過對成熟期5個不同大豆品種的三點彎曲和剪切試驗，得到了最大彎曲力、最大剪切力、慣性矩、最大切應力、彈性模量及抗彎剛度等力學參數(shù)，掌握了大豆莖稈不同高度承受彎曲力及剪切力的一般規(guī)律，為今后大豆收獲機械的研究提供了參考依據(jù)，對研制不同品種通用的大豆收獲機械有重要的意義。

2)大豆受彎曲力及剪切力表現(xiàn)出來的特性和莖稈的直徑變化、不同部位的木質(zhì)化程度相關。

3)阜豆15、荷豆19、柳豆109、皖豆28、皖豆33莖稈整體的彈性模量都是從底部到頂部逐漸增大，造成這一現(xiàn)象的一個原因是莖稈的管狀化程度從底部到頂部越來越明顯，其他方面的原因有待進一步研究。

4)隨大豆莖稈高度的增加，荷豆19、柳豆109、皖豆28、皖豆33的抗彎剛度都是逐漸降低，阜豆15莖稈的抗彎剛度是從底部到頂部先增大、后減小，在中段靠下部位達到最大。

5)隨莖稈高度的增加，5個品種的切應力都呈逐漸減小的趨勢，且從底部到頂部變化較大。

參考文獻：

[1] 梁蘇寧,沐森林,金誠謙,等.黃淮海地區(qū)大豆生產(chǎn)機械化現(xiàn)狀與發(fā)展趨勢[J].農(nóng)機化研究,2015,37(1):261-264，268.

[2] 劉基,金誠謙,梁蘇寧,等.大豆機械收獲損失的研究現(xiàn)狀[J]. 農(nóng)機化研究,2017,39(7):1-9，15.

[3] O 'Dogherty M J. A review of research on forage chopping[J]. J. Agric. Engng. Res., 1982,27(4): 267-289.

[4] Chatopadhyay P S, Pandey K P. Impact cutting behavior of sorghum stalk using a flail-cutter-a mechanical model and its experimental verification[J]. J. Agric. Engng. Res., 2001, 78(4): 369-376.

[5] Hirai Y, Inove E, Mori K. Application of a quasi-static stalk bending analysis to the dynamic response of rice and wheat stalks gathered by a combine harvester reel[J]. Biosystem Engineering, 2004, 88(3): 281-294.

[6] 梁莉,郭玉明. 作物莖稈生物力學性質(zhì)與形態(tài)特性相關性研究[J].農(nóng)業(yè)工程學報,2008(7):1-6.

[7] Hirai Y, Inove E, Mori K, et al. Investigation of mechanical interaction between a combine harvester reel and crop stalk[J].Biosystems Engineering, 2002, 83(3):307-317.

[8] O 'Dogherty M J, Huber J A, Dyson J, et al. A study of the physical and mechanical properties of wheat straw[J].J. Agric. Engng Res., 1995,62(2):133-142.

[9] 閆以勛,趙淑紅,楊悅乾,等.成熟期大豆莖稈力學特性研究[J].東北農(nóng)業(yè)大學學報,2012(5):46-49.

[10] 吳杰，黃勇.棉稈軸向壓縮特性的試驗研究[J].農(nóng)機化研究，2004(4):148-149.

[11] 王艷云，蘇萍.棉稈壓縮力學特性的研究[J].農(nóng)機化研究，2006(6):171-172.

[12] 何曉莉,吳曉強,張立峰,等.大豆莖稈壓縮力學特性的研究[J].農(nóng)機化研究,2010,32(11):164-169.

[13] 李紅波,郭玉明,陳維毅.冬小麥莖稈的力學性質(zhì)研究[J].太原理工大學學報,2006(1):31-34.

[14] 袁志華,馮寶萍,趙安慶,等.作物莖稈抗倒伏的力學分析及綜合評價探討[J].農(nóng)業(yè)工程學報,2002(6):30-31.

[15] 劉慶庭,區(qū)穎剛,卿上樂,等.農(nóng)作物莖稈的力學特性研究進展[J].農(nóng)業(yè)機械學報,2007(7):172-176.

[16] 王芬娥,黃高寶,郭維俊,等.小麥莖稈力學性能與微觀結構研究[J].農(nóng)業(yè)機械學報,2009(5):92-95.

[17] 周楊,李顯旺,沈成,等.工業(yè)大麻莖稈力學模型的試驗分析(英文)[J].農(nóng)業(yè)工程學報,2016(9):22-29.

[18] 陳海濤,頓國強. 基于虛擬樣機動力學仿真的大豆扶禾器參數(shù)優(yōu)化[J].農(nóng)業(yè)工程學報,2012(18):23-29.

[19] 郭玉明,袁紅梅,陰妍,等.莖稈作物抗倒伏生物力學評價研究及關聯(lián)分析[J].農(nóng)業(yè)工程學報,2007(7):14-18.

[20] 劉慶潭.材料力學[M].北京:機械工業(yè)出版社，2003:8-27,101.