田輝 袁志華 陳永 吳祎迪 杜偉 張駐軍 時(shí)永磊

摘要：小麥倒伏每年造成小麥產(chǎn)量損失約10%～30%，研究小麥倒伏對提高小麥莖稈品質(zhì)和抗倒伏特性具有重要意義。在小麥倒伏的影響因素中，莖稈的力學(xué)性能至關(guān)重要。在力學(xué)理論和有限元模擬的基礎(chǔ)上，研究小麥莖稈力學(xué)性能隨維管束分布的變化規(guī)律。結(jié)果表明，維管束分布對倒伏的影響與維管束和基礎(chǔ)組織的彈性模量有關(guān)。在維管束彈性模量大于基本組織彈性模量的情況下，維管束越多，越難倒伏;當(dāng)維管束數(shù)量相同時(shí)，維管束越靠近表皮，越難倒伏。增大維管束的彈性模量和維管束所在處的圓周直徑，可以增強(qiáng)小麥莖稈的抗倒伏能力。

關(guān)鍵詞：小麥倒伏;維管束;撓度;抗彎剛度;有限元模擬

中圖分類號(hào)： S512.101 ?文獻(xiàn)標(biāo)志碼： A ?文章編號(hào)：1002-1302（2021）13-0181-05

小麥?zhǔn)鞘澜缟现匾霓r(nóng)作物，每年的產(chǎn)量為5.945億t，可滿足世界35%～40%的人口需求。倒伏導(dǎo)致小麥產(chǎn)量的損失每年約占產(chǎn)量總損失的10%～30%[1]，以平均值20%作為評估基礎(chǔ)，每年倒伏造成的總損失約為1.189億t，是我國小麥總產(chǎn)量的1.09倍。因此，研究抗倒伏性，尋找提高小麥莖稈品質(zhì)和抗倒伏特性的新途徑具有重要意義。目前解決小麥倒伏問題的方法是使用矮稈品種。以矮稈品種為代表的綠色革命，解決了億萬人的糧食問題。但是，世界人口的增長對糧食生產(chǎn)提出了新的要求。傳統(tǒng)的矮稈品種由于生物量積累量低，進(jìn)一步提高產(chǎn)量受到限制。增加小麥株高可以提高生物產(chǎn)量，但是增加株高會(huì)帶來新的倒伏問題。為了提高生物產(chǎn)量，必須提高植株的抗倒伏能力。

影響小麥莖稈倒伏的因素有很多[2-3]，但歸根結(jié)底與小麥莖稈的力學(xué)性質(zhì)有關(guān)，小麥莖稈力學(xué)性質(zhì)又由其形態(tài)、結(jié)構(gòu)、組成成分等多種因素決定。對小麥莖稈各節(jié)間做拉伸、彎曲、剪切、扭轉(zhuǎn)等一系列力學(xué)試驗(yàn)，得到了小麥莖稈的拉伸強(qiáng)度、剪切強(qiáng)度、楊氏模量、彎曲剛度、剪切模量、彎曲強(qiáng)度等參數(shù)。這些參數(shù)隨品種、生長期、化學(xué)成分、節(jié)點(diǎn)間位置、含水率、形態(tài)和微觀組織結(jié)構(gòu)而變化[4-5]。Baker等建立了作物莖稈的力學(xué)模型，闡明了倒伏過程，并用簡單的圖形公式計(jì)算了倒伏風(fēng)險(xiǎn)[6]。Leblicq等利用數(shù)字高程模型（DEM）建立了作物莖稈模型，考慮了塑性變形和損傷的影響，模擬分析了單個(gè)稈件的彎曲性能[7]。Huang等研究表明，纖維束的平均拉伸模量和拉伸強(qiáng)度明顯低于表皮層，抗拉剛度沿稈的徑向和軸向呈梯度分布，這種梯度分布增加了莖稈基礎(chǔ)的剛度，這是影響玉米植株抗倒伏的關(guān)鍵因素[8]。關(guān)于小麥倒伏的影響因素，莖的力學(xué)性能至關(guān)重要，彎曲剛度可以用作抗倒伏性的評價(jià)指標(biāo)[9]。汪燦等在蕎麥莖稈解剖結(jié)構(gòu)及木質(zhì)素代謝與抗倒性關(guān)系的研究中指出，莖稈中木質(zhì)素含量、機(jī)械組織厚度、莖壁厚度、大維管束數(shù)目和維管束面積可以作為禾本科莖稈抗倒伏育種的重要指標(biāo)[10]。提高木質(zhì)素含量可以提高玉米莖稈的抗倒伏能力[11]。在溫室環(huán)境中應(yīng)用機(jī)械刺激可用于改變植物生長，實(shí)現(xiàn)更緊密的密實(shí)度、更大的周長[12]。有研究表明，小麥莖稈的彎曲剛度與成熟度、形態(tài)、結(jié)構(gòu)、木質(zhì)素和纖維素含量有關(guān)[13-14]。有文獻(xiàn)研究了不同品種小麥第2節(jié)間的抗壓強(qiáng)度，結(jié)果表明，不同小麥品種的抗壓性隨機(jī)械組織的增厚而增強(qiáng)。機(jī)械組織的纖維層越厚，維管束的數(shù)量越多，小麥莖的強(qiáng)度、剛度和穩(wěn)定性越好，實(shí)心小麥具有較強(qiáng)的抗倒伏性[15]。還有研究表明，倒伏指數(shù)與株高、重力高度和莖長呈極顯著正相關(guān)關(guān)系，而與莖徑、壁厚、莖密度、彎曲強(qiáng)度、維管束數(shù)量和面積呈負(fù)相關(guān)關(guān)系。相關(guān)研究表明，莖的機(jī)械組織與莖的直徑、壁厚、薄壁組織、彎曲強(qiáng)度以及維管束的數(shù)量和面積呈正相關(guān)關(guān)系[16]。相關(guān)研究表明，莖稈維管束結(jié)構(gòu)特征與穗粒數(shù)、結(jié)實(shí)率和產(chǎn)量成正相關(guān)，而與千粒質(zhì)量呈負(fù)相關(guān)，特別是上部莖節(jié)間小維管束和基部莖節(jié)間大維管束結(jié)構(gòu)特征顯著影響超級(jí)雜交稻產(chǎn)量及產(chǎn)量構(gòu)成要素[17]。

本研究基于力學(xué)理論，分析抗彎剛度與維管束數(shù)量和面積的關(guān)系。在力學(xué)理論和有限元分析的基礎(chǔ)上，進(jìn)行了抗彎剛度與維管束分布之間的關(guān)系研究。研究結(jié)果以期為小麥高產(chǎn)栽培和抗倒伏性研究提供參考。

1 材料與方法

1.1 理論分析

對小麥莖稈節(jié)間三點(diǎn)彎曲進(jìn)行分析。節(jié)點(diǎn)間長度為l，P是荷載，EI是抗彎剛度，δ為稈的最大撓度。根據(jù)材料力學(xué)理論，稈的最大撓度與抗彎剛度的關(guān)系為公式（1）：

1.2 有限元模擬

有限元模擬可以在任何時(shí)間、任何地點(diǎn)進(jìn)行，可以用來研究小麥莖稈的力學(xué)特性。ANSYS軟件是融結(jié)構(gòu)、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件，是目前國際上最流行的有限元分析軟件之一[18]。它操作簡單且功能強(qiáng)大，可預(yù)先模擬以確保實(shí)際工程應(yīng)用安全[19]。

根據(jù)小麥莖稈截面結(jié)構(gòu)的實(shí)際情況，利用有限元軟件ANSYS17.0建立小麥莖稈的ansys模型。橫截面可分為2層，外層代表表皮，內(nèi)層含有大的維管束。適當(dāng)擴(kuò)大試驗(yàn)數(shù)據(jù)，模型幾何尺寸如表1所示。模型1、模型2、模型3的維管束均勻分布在同一圓周上，模型4、模型5、模型6的維管束交替均勻分布在不同圓周上。模型7沒有維管束。在這7個(gè)模型中，只有維管束的數(shù)目和分布不同，其余參數(shù)都是相同的。表1中表皮材料的泊松比為v1，基本組織材料的泊松比為v2，維管束的泊松比為v3。E1為 22 930 MPa，E2為60 MPa，E3為29 400 MPa。單元類型是Solid186。模型1橫截面網(wǎng)格劃分如圖5所示，模型1截面變形如圖6所示，模型4橫截面網(wǎng)格劃分如圖7所示，模型4截面變形如圖8所示。三點(diǎn)彎曲試驗(yàn)的模擬結(jié)果如表2所示。

2 結(jié)果與分析

2.1 維管束數(shù)量對小麥倒伏的影響

由公式（5）可知，n越大，I3越大。

在E3>E2的情況下，由公式（3）可知，I3越大，EI值越大。由公式（1）可知，EI越大，撓度越小。撓度越小，莖稈就越難倒伏。因此，對于E3>E2，n越大，莖稈就越難倒伏。

在E3 ? ? 維管束數(shù)量對倒伏的影響與維管束彈性模量有關(guān)。如果維管束的彈性模量大于基本組織的彈性模量，則維管束越多，越難倒伏。如果維管束的彈性模量小于基本組織的彈性模量，維管束越少，就越難倒伏。有研究表明，在新鮮狀態(tài)下，玉米秸稈內(nèi)部的彈性模量為0.06 GPa[20];大維管束的平均彈性模量為4.44 GPa，小維管束的平均彈性模量為10.80 GPa[8]。根據(jù)這些研究結(jié)果，在有限元模擬中，假設(shè)維管束的彈性模量大于基本組織的彈性模量。在ANSYS模型中，小麥莖稈的維管束彈性模量大于基本組織的彈性模量。仿真結(jié)果與理論分析結(jié)果吻合。

2.2 維管束位置對小麥倒伏的影響

n是大維管束的數(shù)目。這些維管束均勻分布在同一圓周上，Dw為該圓的直徑。I3是這些所有維管束的慣性矩。

維管束數(shù)目仍為n，其中（n-1）個(gè)維管束仍在直徑Dw的圓周上，第n個(gè)維管束的中心位置在直徑為Dw1的圓周上，Dw1 ? ? 維管束數(shù)目仍為n，有一半均勻分布在直徑為Dw的周周上。另一半交替均勻分布在直徑為Dw1圓周上，Dw1 ? ? 由公式（7）、公式（6）、公式（4）可得I*3E2，I3越大，EI值越大。由公式（1）可知，EI越大，撓度越小。因此，可以得出，當(dāng)維管束均勻分布在同一圓周上時(shí)，莖稈較難倒伏。在維管束數(shù)目相同的情況下，Dw越大，莖稈越難倒伏。

Dw是維管束所在圓周的直徑。Dw越大，即維管束離表皮越近。因此，可以得到維管束越靠近表皮，在維管束數(shù)目相同的情況下，莖稈越難倒伏。

利用有限元軟件ANSYS17.0建立了小麥莖稈的7個(gè)模型，在這7個(gè)模型中，只有維管束數(shù)量和分布不同，其余參數(shù)均相同。結(jié)果表明，當(dāng)維管束數(shù)目相同時(shí)，在相同圓周上的撓度較小，在不同圓周上的撓度較大。仿真結(jié)果與理論結(jié)果基本一致。

3 討論

在影響小麥倒伏的因素中，莖稈的力學(xué)性能非常重要。莖稈的力學(xué)性能與小麥品種、生育期、組成、形態(tài)和解剖結(jié)構(gòu)有關(guān)。彎曲剛度可用作莖稈抗倒伏性的評估指標(biāo)。抗彎剛度越大，抗倒伏性越強(qiáng)?？箯潉偠仍酱螅o稈越難倒伏。許多研究表明，增加維管束的面積和數(shù)量可以顯著降低倒伏指數(shù)并提高莖的抗倒伏能力[15-16，21]，本研究結(jié)果與該結(jié)果相同。本研究結(jié)果還表明，維管束數(shù)量對倒伏的影響與維管束的彈性模量有關(guān)。如果維管束彈性模量大于基本組織的彈性模量，則維管束越多，莖稈越難倒伏;如果維管束彈性模量小于基本組織的彈性模量，則維管束越少，莖稈越難倒伏。維管束的彈性模量越大，彎曲剛度越大。增加維管束的彈性模量可以提高莖稈的抗倒伏能力。維管束位置對倒伏的影響與維管束的彈性模量有關(guān)。在維管束數(shù)量相同的情況下，如果維管束的彈性模量大于基本組織的彈性模量，則維管束均勻地分布在同一圓周上時(shí)，莖稈越難倒伏。維管束離表皮越近，莖稈越難倒伏。如果維管束的彈性模量小于基本組織的彈性模量，則維管束離表皮越遠(yuǎn)，莖稈越難倒伏。

根據(jù)復(fù)合材料力學(xué)理論，推導(dǎo)出小麥莖稈抗彎剛度的計(jì)算公式。該公式可用于分析維管束分布對小麥倒伏的影響。利用有限元軟件ANSYS建立了小麥莖稈的7個(gè)力學(xué)模型，對三點(diǎn)彎曲進(jìn)行了模擬。仿真結(jié)果與理論結(jié)果基本一致。說明該建模方法、單元類型選取是可行的。

4 結(jié)論

理論分析和模擬結(jié)果表明，小麥莖稈的力學(xué)性能與維管束分布等諸多因素有關(guān)。維管束分布對倒伏的影響與維管束和基礎(chǔ)組織的彈性模量有關(guān)。在維管束彈性模量大于基本組織彈性模量的情況下，維管束越多，莖稈越難倒伏。當(dāng)維管束的數(shù)量相同時(shí)，維管束越靠近表皮，莖稈越難倒伏。通過增加維管束的彈性模量和維管束所在處的圓周直徑，可以增強(qiáng)小麥莖稈的抗倒伏能力。

參考文獻(xiàn)：

[1]李金才，尹 鈞，魏鳳珍. 播種密度對冬小麥莖稈形態(tài)特征和抗倒指數(shù)的影響[J]. 作物學(xué)報(bào)，2005，31（5）：662-666.

[2]Martinez-Vazquez P. Crop lodging induced by wind and rain[J]. Article in Agricultural and Forest Meteorology，2016，228/229：265-275.

[3]Piera-Chavez F J，Berry P M，F(xiàn)oulkes M J，et al. Avoiding lodging in irrigated spring wheat. Ⅱ. Genetic variation of stem and root structural properties[J]. Field Crops Research，2016，196：64-74.

[4]Leblicq T，Vanmaercke S，Ramon H，et al. Mechanical analysis of the bending behaviour of plant stems[J]. Biosystems Engineering，2015，129：87-99.

[5]Peake A S，Huth N I，Carberry P S，et al. Quantifying potential yield and lodging-related yield gaps for irrigated spring wheat in sub-tropical Australia[J]. Field Crops Research，2014，158：1-14.

[6]Baker C J，Sterling M，Berry P. A generalised model of crop lodging[J]. Journal of Theoretical Biology，2014，363（7）：1-12.

[7]Leblicq T，Smeets B，Vanmaercke S，et al. A discrete element approach for modelling bendable crop stems[J]. Computers and Electronics in Agriculture，2016，124：141-149.

[8]Huang J L，Liu W Y，Zhou F，et al. Mechanical properties of maize fibre bundles and their contribution to lodging resistance[J]. Biosystems Engineering，2016，151：298-307.

[9]梁 莉，郭玉明. 不同生長期小麥莖稈力學(xué)性質(zhì)與形態(tài)特性的相關(guān)性[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24（8）：131-134.

[10]汪 燦，阮仁武，袁曉輝，等. 蕎麥莖稈解剖結(jié)構(gòu)和木質(zhì)素代謝及其與抗倒性的關(guān)系[J]. 作物學(xué)報(bào)，2014，40（10）：1846-1856.

[11]Kamran M，Ahmad I，Wang H Q，et al. Mepiquat chloride application increases lodging resistance of maize by enhancing stem physical strength and lignin biosynthesis[J]. Field Crops? Research，2018，224：148-159.

[12]Brnke F，Rocksch T. Thigmomorphogenesis-control of plant growth by mechanical stimulation[J]. Scientia Horticulturae，2018，234：344-353.

[13]潘 婷，胡文靜，李東升，等. 小麥莖稈實(shí)心度對莖稈強(qiáng)度的影響及相關(guān)性狀QTL分析[J]. 作物學(xué)報(bào)，2017，43（1）：9-18.

[14]郭維俊，王芬娥，黃高寶，等. 小麥莖稈力學(xué)性能與化學(xué)組分試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009，40（2）：110-114.

[15]Kong E Y，Liu D C，Guo X L，et al. Anatomical and chemical characteristics associated with lodging resistance in wheat[J]. The Crop Journal，2013，1（1）：43-49.

[16]文廷剛，王偉中，楊文飛，等. 水稻莖稈形態(tài)特征與抗倒伏能力對外源植物生長調(diào)節(jié)劑的響應(yīng)差異[J]. 南方農(nóng)業(yè)學(xué)報(bào)，2020，51（1）：48-55.

[17]彭 廷，可文靜，熊加豹，等. 超級(jí)雜交稻持續(xù)增產(chǎn)的莖稈維管束結(jié)構(gòu)基礎(chǔ)研究[J]. 華北農(nóng)學(xué)報(bào)，2019，34（1）：165-171.

[18]王亞利. ANSYS軟件在機(jī)械結(jié)構(gòu)分析中的應(yīng)用[J]. 價(jià)值工程，2014，33（3）：203-203，204.

[19]楊少凡. 淺談ANSYS在平面應(yīng)力問題分析中的應(yīng)用[J]. 四川建材，2020，46（8）：41-42.

[20]張李嫻. 玉米秸稈力學(xué)特性的離散元建模方法研究[D]. 楊凌：西北農(nóng)林科技大學(xué)，2017.

[21]田文濤，邵 平，王 燚，等. 超級(jí)雜交稻莖稈形態(tài)結(jié)構(gòu)及其與抗倒性的關(guān)系研究[J]. 雜交水稻，2017，32（2）：67-71，81.