籽棉壓縮與應力松弛力學特性及模型構建

孔凡婷，吳 騰，陳長林，孫勇飛，謝 慶，石 磊

籽棉壓縮與應力松弛力學特性及模型構建

孔凡婷，吳 騰，陳長林，孫勇飛，謝 慶，石 磊※

（農業農村部南京農業機械化研究所，南京 210014）

為深入研究籽棉壓縮及應力松弛過程的力學特性并構建其本構模型，該研究以籽棉為研究對象開展試驗，分別利用改進西原模型和五元件廣義Maxwell模型對壓縮和應力松弛的應力應變曲線進行描述。通過定義的本構模型對不同含水率及喂入量的籽棉力學特性試驗數據結果進行參數辨識，得到相關模型參數，并探究不同因素對壓縮及松弛過程應力的影響規律。研究結果表明：曲線擬合法求解籽棉壓縮及應力松弛過程本構模型參數的決定系數均大于0.9，改進西原模型、五元件廣義Maxwell模型可以較好地描述籽棉壓縮力學特性和應力松弛特性。對籽棉壓縮及應力松弛力學特性解析顯示，模型參數表現出明顯的應力規律：通過Duncan均值比較結果可知：籽棉壓縮本構模型參數中不同含水率組間彈性模量差異均顯著（<0.05），不同喂入量組間串聯虎克體彈性模量差異顯著（<0.05），壓縮應力值與含水率、喂入量呈正相關關系；籽棉應力松弛模型參數中不同含水率及喂入量組間彈性模量及粘性模量差異均顯著（<0.05），彈性模量、粘性模量與含水率及喂入量均呈正相關關系。籽棉壓縮與應力松弛力學特性的研究可為籽棉壓縮過程機理研究、模擬仿真提供理論依據，縮短機具研發過程。

壓縮；應力松弛；本構模型；試驗研究；籽棉

0 引 言

棉花是世界上重要的經濟作物和紡織工業原料，在經濟發展中占據重要地位[1]。中國棉花種植面積超300萬hm2，無論是種植面積還是產量，都位于世界前列[2]。在當前城市化進程不斷加快、大量農村勞動力向城市轉移的背景下，自2007年以來，機采棉的發展速度逐年加快。隨著機采棉面積的不斷增加，傳統、落后的籽棉儲運技術暴露出來的問題愈發明顯。

目前，籽棉儲運主要有散裝、打垛和打包3種方式。散裝儲運耗費勞力多、運輸效率低、占用場地大、存放期間易造成棉花霉爛。打垛可提高籽棉運輸效率，但需專用的打垛設備和拉垛車輛，工序繁瑣，不便于移動。籽棉打包隨著美國迪爾CP690型摘棉打包機的引進率先在新疆地區應用，其技術能顯著提高儲運效率，降低籽棉運輸成本；棉花品質不因天氣變化影響，也可杜絕火災發生，還有效防止異性纖維混入；圓形包捆移動和堆放方便，占用場地空間小并無特殊儲運要求[3]。籽棉打包儲運技術優勢突出，先進性和實用性已得到生產實踐認可，是籽棉儲運方式發展的一個重要方向。在中美貿易戰的大背景下，研發具有自主知識產權的打包裝備，可解決中國棉花儲運難的問題。

籽棉壓縮過程的力學特性是籽棉收獲、打包等裝備優化設計的基礎，尤其是打包機研發過程中，籽棉壓縮力學特性對打包機預壓過程中預壓輥的結構及工藝參數、打包過程中喂入量與打包帶轉速優化匹配、打包輥結構、布置與工作參數及包裝過程中包膜機構的結構與工藝參數、打包膜物理特性要求等影響較大。但目前關于農業物料壓縮力學特性的研究多集中于玉米、小麥或牧草秸稈等物料。Sun等[4]、邢獻軍等[5]、馬彥華等[6]、Huang等[7]分別對玉米秸稈進行了不同的力學特性試驗研究；馬瑞峻等[8]、Nona等[9]、 Zhang等[10]、吳敏等[11]、房佳佳等[12]學者分別對小麥、水稻秸稈及苜蓿類牧草等物料進行了壓縮力學試驗及分析；朱凱等[13]進行了粉體棉花秸稈應力松弛特性試驗并將所得數據擬合分析；Yang等[14]、劉慶庭等[15]學者對竹子、甘蔗等硬質秸稈物料的壓縮力學特性展開研究與分析；部分學者對去籽皮棉壓縮特性[16-18]進行研究，Hardin等[19-20]對籽棉壓縮和蠕變開展研究，其主要研究打垛過程中籽棉力學特性，壓縮密度范圍較小（64～128 kg/m3），不能為打包裝備研發提供理論參考。目前，各種物料壓縮或成型的建模以及對成型塊質量的影響因素分析已取得一定進展，霍麗麗等[21]、陳嘯等[22]、馬方等[23]、Kaliyan等[24]進行了生物質物料的本構模型構建研究及分析，并為生產實踐提供了良好的理論支撐。

基于此，本文以籽棉為研究對象開展壓縮力學特性試驗，以含水率及喂入量為試驗因素利用萬能試驗機進行壓縮及應力松弛試驗，將得到的應力應變數據結果進行曲線擬合，并對定義的本構模型進行參數辨識，進而進行籽棉壓縮力學特性研究。擬為籽棉的壓縮過程機理研究、模擬仿真及籽棉打包機構的優化設計提供理論參考。

1 材料與方法

1.1 試驗材料與設備

1.1.1 試驗材料

試驗籽棉種植于農業農村部南京農業機械化研究所白馬基地，品種為中6913。試驗樣品為于2019年10月28日農業農村部南京農業機械化研究所研制的4MZ-3型刷輥式采棉機收獲籽棉，收獲時籽棉含水率為10.8%。

1.1.2 試驗設備

試驗使用的主要設備為WDW-10 微機控制電子萬能試驗機及自制壓縮模具，如圖1所示。萬能試驗機試驗力示值精度在±1%內，最大試驗力為10 kN；壓縮模具主要包括壓縮腔及壓縮桿：壓縮腔由PVC管制作與底板螺紋連接，便于拆卸；壓縮桿與壓縮端銷軸連接。其他輔助儀器包括南京試驗儀器廠生產的DGH30-1A型電熱鼓風干燥箱；中國上海立晨儀器技術有限公司生產的YP300001D電子秤（量程0～30kg，精度0.1 g）；游標卡尺、直尺和秒表等。

1.2 試驗方法

本文對籽棉進行單向施力擠壓，以籽棉在受壓狀態下的應力應變變化作為表征籽棉壓縮性能的重點指標。預試驗結果顯示，籽棉的棉纖維屬于應變敏感材料，靜態壓縮條件下加載速度對籽棉壓縮應力無明顯影響；而喂入量及含水率對壓縮及應力松弛力學特性值具有較大影響，因此選擇喂入量及含水率作為籽棉壓縮及應力松弛力學特性試驗的因素。收獲時因采棉機采摘原理不同及季節氣候影響，收獲籽棉含水率范圍在10%～30%，因此將試樣進行預調濕處理，將每組試樣放置在溫度為50 ℃、濕度為10%的環境下預調濕2h，然后在溫度20 ℃、相對濕度6%的恒溫恒濕條件下調濕24h后進行測試[25]。

籽棉壓縮及應力松弛試驗均在萬能試驗機上進行，為減少誤差，基于試驗條件將0.1、0.15、0.2 kg不同喂入量籽棉自由堆放至壓縮腔，為保證籽棉壓縮效果并兼顧效率，以30 mm/min壓縮速度分別對10%、20%、30%含水率的籽棉進行試驗。每組試驗重復3次。對不同含水率的籽棉試驗數據，分別以喂入量-序號進行重復性試驗編號。

1.2.1 壓縮試驗

約翰迪爾打包采棉一體機打包密度為220 kg/m3，預試驗將不同含水率的籽棉壓縮至約400 kg/m3時棉籽產生破壞。因此，設定壓縮桿從籽棉自由堆積密度87.73 kg/m3壓縮至指定密度250 kg/m3后停止加載，由萬能試驗機記錄加載過程的負荷值，由此表征籽棉在整個壓縮過程中的受力情況。壓縮過程中將籽棉視為一個集合體，在萬能試驗機輸入壓縮截面積及壓縮位移量，試驗軟件將測得負荷值及形變量自動將坐標變換可得到應力及應變值。

1.2.2 應力松弛試驗

每組壓縮試驗停止加載后，壓縮桿不動保持籽棉變形，由萬能試驗機記錄其應力變化值，90s后結束試驗，得到籽棉壓縮后應力松弛曲線。

1.3 試驗結果

1.3.1 壓縮試驗數值

以30mm/min壓縮速度分別對10%、20%、30%含水率下0.1、0.15、0.2kg喂入量的籽棉進行壓縮，每組試驗重復3次。由圖2可知，各組重復試驗得到的應力應變曲線特征基本一致，且離散性不大。

注：圖中0.1kg-、0.15 kg-、0.2 kg-(=1,2,3)分別代表0.1 kg、0.15 kg、0.2 kg喂入量的3組重復試驗編號，下同。

Note: The 0.1 kg-, 0.15 kg-, 0.2 kg-(=1,2,3) stand for test No. were repeated in three groups of 0.1 kg, 0.15 kg, 0.2 kg feeding quantity, The same below.

圖2 不同含水率籽棉壓縮應力應變曲線

Fig.2 Compression stress-strain curve of seed cotton with different moisture content

棉籽上均布棉纖維的籽棉與其他農業物料相比，壓縮過程中的力學變化比較復雜，主要表現為粘彈特性。剛開始壓縮時，籽棉纖維間的間隙逐漸減小，主要為棉纖維內部纖維束發生彎折彈性變形及纖維間內摩擦的粘性，此階段應力近似線性緩慢增加以排除纖維間空隙為主，定義為排空階段；該階段曲線表現為近似線性，此過程物料松散變形雖大，但應力變化很小。隨著壓縮空間減小，單向壓力從纖維傳遞至棉籽，棉籽在力的作用下產生滑移及變形，因棉籽排列方式的隨機性，此階段應力變化加速且可能有階躍式變化，內部纖維體受壓穿插、滑移，進行內部調整，定義為滑移變形階段，應力應變曲線呈逐漸上升趨勢。當籽棉壓縮密度達到約150 kg/m3后繼續壓縮，此階段內部應力隨應變增加而急劇上升，集合體內部纖維排列緊密，主要以棉纖維自身受壓發生部分塑性變形、物料間內摩擦的粘性及棉籽受壓的彈性變形為主，該階段定義為非線性粘彈塑性變形階段；此階段曲線表現為快速上升，即壓縮應變較小時應力快速上升。

1.3.2 應力松弛試驗數值

以30mm/min壓縮速度分別對10%、20%、30%含水率；0.1、0.15、0.2 kg不同喂入量的籽棉壓縮停止加載后，壓縮桿保持90 s變形得到的應力松弛曲線，每組試驗重復3次。由圖3可知，各組重復試驗得到的曲線特征基本一致，即應力隨時間呈先快后慢的非線性減少，最后曲線趨于恒定值，這與秸稈、草料等農業物料應力松弛取得結果一致[8,12]。

2 模型構建

研究借鑒相關生物質壓縮成型及秸稈物料壓縮建模方法進行籽棉壓縮本構模型構建。

2.1 壓縮本構模型

西原體模型[26]是由虎克體、開爾文體和粘塑性體串聯組成，適合反映軟材料的力學特性，其力學模型如圖4a所示。西原模型的本構方程為

但西原體模型在描述壓縮第1階段僅為彈性，無法真實地反映出籽棉壓縮過程中的粘彈特性。為更好的與實際材料曲線相吻合，建立統一反映籽棉壓縮3個階段的模型，本文添加一個牛頓體與原西原體模型進行串聯，來替代傳統西原模型中理想粘-彈-塑性體，得到改進西原模型來反映籽棉應力的加速特征，其力學模型如圖4b所示。

注：為總應力，MPa；0為串聯牛頓體的粘性模量，MPa·s；1為串聯虎克體的彈性模量，MPa；2為開爾文體中虎克體的彈性模量，MPa；2為開爾文體中牛頓體的粘性模量，MPa·s；σ為庫倫體的初始應力，MPa；3為粘塑體中牛頓體的粘性模量，MPa·s。

Note:is total stress, MPa;0is viscosity modulus of tandem Newtonian body, MPa·s;1is elasticity modulus of tandem Hookean body, MPa;2is elasticity modulus ofHookean body in Kelvin body, MPa;2is viscosity modulus of Newtonian body in Kelvin body, MPa·s;σis initialstress of Cullen body, MPa;3is viscosity modulus of Newtonian body in Viscoplastic body, MPa·s.

圖4 壓縮力學模型

Fig.4 Compression mechanical model

在推導籽棉壓縮本構方程時，可將牛頓體、虎克體、開爾文體和粘塑體看成單個元件，然后應用串聯運算規則[26]，即可求得改進西原模型體的本構方程。

當<σ時，

0=1=2（2）

=0+1+2（3）

將式（3）對時間求一階導、二階導分別為（以下公式求導均為對時間求導）

式中為總應力，MPa；為總應變；0、1、2分別為串聯元件牛頓體、虎克體和開爾文體的應力，MPa；0、1、2分別為串聯元件牛頓體、虎克體和開爾文體的應變。

對于牛頓體有

求一階導為

對于虎克體有

求一階導、二階導分別為

對于開爾文體有

求一階導為

將公式（6）、（9）代入式（4）解得：

將公式（7）、（10）、（12）、（13）代入式（5）解得：

當>σ時，

式中3為串聯粘塑體的應力，MPa；3為串聯元件粘塑體的應變。

將式（15）求一階導、二階導分別為

對于粘塑性體有：

求一階導為

式中，3為開爾文元件中牛頓體的粘性模量，MPa·s。

將公式（6）、（9）、（18）代入式（16）解得：

將公式（7）、（10）、（12）、（19）、（20）代入式（17）解得：

以含水率20%、喂入量0.15kg第2次重復試驗數據進行擬合，擬合效果如圖5所示。由擬合曲線觀察到，西原模型無法擬合出籽棉壓縮曲線加速階段特征，其決定系數2小于0.85；改進西原模型決定系數2大于0.98，與西原模型相比，在應力水平超過某一應變值后應力由緩慢增加變為非線性快速增加，表現出明顯的粘彈性。因此，選擇改進西原模型可較全面地表征籽棉的壓縮特性，可選用該模型來構建籽棉壓縮階段本構模型。

2.2 應力松弛本構模型

籽棉壓縮的應力松弛階段應變不變，應力逐漸減小，由于籽棉內部存在彈性恢復力，使彈性變形逐漸變成非彈性變形。應力松弛階段一般采用廣義Maxwell模型描述[27-29]，從圖3知不同壓縮條件下的籽棉應力松弛曲線的變化規律基本一致，且曲線的形狀相似，所以應力松弛曲線可以用同一種模型進行描述。

三元件廣義Maxwell力學模型主要由1個彈簧和1個Maxwell元件并聯組成圖6a，三元件廣義Maxwell模型松弛階段方程為[30]

式中()為時間在時刻的應力，MPa；ε為初始應變；為時間s，T為Maxwell元件的應力松弛時間T=η/E，s；η為牛頓體的粘性模量，MPa·s。

五元件廣義Maxwell力學模型主要由1個彈簧和2個Maxwell元件并聯組成，見圖6b，五元件廣義Maxwell模型應力松弛方程為[30]

式中T、T分別為Maxwell元件的應力松弛時間，T=η/E，T=η/E，s；η、η分別為牛頓體的粘性模量，MPa·s。

注：E為虎克體的彈性模量，MPa；為Maxwell體中牛頓體的粘性模量，MPa·s；E為Maxwell體中虎克體的彈性模量，MPa；為并聯Maxwell體中牛頓體的粘性模量，MPa·s；E為并聯Maxwell體中虎克體的彈性模量，MPa。

Note:Eis elasticity modulus of Hookean body, MPa;is viscosity modulus of Newtonian body in Maxwell body, MPa·s;Eis elasticity modulus ofHookean body inMaxwell body, MPa;is viscosity modulus of Newtonian body in parallel Maxwell body, MPa·s;Eis elasticity modulus ofHookean body in parallel Maxwell body, MPa.

圖6 應力松弛力學模型

Fig.6 Stress relaxation mechanical model

公式（23）和公式（24）可以看出五元件廣義Maxwell模型中既包含第1應力松弛時間，又包含第2應力松弛時間；而三元件廣義Maxwell模型中僅包含一個應力松弛時間T。因此，選擇五元件廣義Maxwell模型可較全面地表征籽棉的應力松弛特性。以含水率20%、喂入量為0.15kg第2次重復試驗數據進行擬合，擬合效果如圖7所示。五元件廣義Maxwell模型擬合決定系數2大于0.99，因此選用五元件來構建松弛階段本構模型。

3 參數確定與分析

3.1 參數確定

基于籽棉在壓縮過程中經歷的彈性變形和塑性變形及籽棉纖維間固有的粘阻性和受壓后表現出的粘結性，將籽棉壓縮應力應變曲線分為近似線性、逐漸上升和快速上升3個階段。為較好的擬合籽棉壓縮力學特性曲線采用改進西原模型進行建模，將應力應變試驗數據導入Origin Pro 8軟件并輸入籽棉本構模型，以曲線擬合法求解壓縮本構模型的參數值；將求解的模型參數值輸入SPSS 22軟件進行數據處理得到不同因素水平下籽棉壓縮的模型參數，利用Duncan均值多重比較法對結果進行對比分析，結果如表1所示。

表1 含水率和喂入量對籽棉壓縮本構模型參數的影響

注：采用Duncan均值多重比較法進行比較，同列肩標上字母表示不同含水率（小寫字母）及不同喂入量（大寫字母）組間差異顯著（<0.05）。下同。

Note: The results are compared by Duncan’s multiple comparisons, different letters on shoulder indicate significant difference among difference moisture content groups (lowercase letters) and difference feeding quantity groups (uppercase letter) (<0.05). The same below.

從表1可知采用改進西原模型擬合籽棉壓縮曲線其決定系數2均大于0.9，結果表明改進西原模型可以很好的描述籽棉壓縮過程力學特性。

以曲線擬合法分別求解10%、20%、30%含水率及0.10、0.15、0.20 kg喂入量的試驗條件下籽棉應力松弛模型參數值，再利用Duncan均值多重比較法進行對比分析，分析結果見表2。

從表2可知采用五元件廣義Maxwell模型擬合籽棉應力松弛曲線其決定系數2均大于0.99，結果表明五元件廣義Maxwell模型可以很好的描述籽棉壓縮松弛力學特性。

改進西原模型、應力松弛模型及其對應參數值，利用通過ANSYS/LS-DYNA或ABAQUS等仿真軟件二次開發，可進行籽棉壓縮與相關機構間的仿真模擬，以縮短機具研發過程。

表2 含水率和喂入量對籽棉應力松弛模型參數的影響

3.2 各因素對壓縮力學特性影響

3.2.1 含水率

3.2.2 喂入量

由表1擬合參數結果可知，彈性模量隨喂入量的增加而增加，Duncan均值比較結果顯示不同喂入量組間串聯虎克體彈性模量1差異顯著（<0.05），可能由于棉纖維的卷曲使纖維在壓縮過程中縱向收縮并具有彈性伸長，從而增加籽棉整體縱向可變形性，提高其彈性模量。喂入量對三組粘性模量及粘塑性初始應力的正相關影響規律表明，在相同的壓縮腔內自然堆積的籽棉，喂入量越大籽棉堆積高度越高，壓力傳輸路徑越長，籽棉壓縮過程中纖維體彎曲變形，纖維體內部越密實，纖維體之間相互接觸、擠壓產生的局部應力越大，即粘性模量增加。

3.3 各因素對應力松弛力學特性影響

3.3.1 含水率

表2中結果顯示，Duncan均值比較結果顯示不同含水率組間彈性模量及粘性模量差異均顯著（<0.05），含水率與籽棉彈性模量、粘性模量均呈正相關關系。主要原因為棉的天然轉曲使纖維相互糾纏、掛鉤和粘著，使得籽棉應力松弛過程中不易分散，而隨含水率加大，水分子滲透到棉纖維間，棉纖維發生膨潤現象，彈性模量增大[32]。且水分使纖維間的抱合力大大增加，增大了運動阻力，粘性模量值隨之快速增長。

3.3.2 喂入量

由表2擬合參數結果可知，模型系數表現出明顯的應力規律，Duncan均值比較結果顯示不同喂入量組間彈性模量及粘性模量差異均顯著（<0.05）。隨著喂入量的增加，彈性模量數值增大，其原因主要為籽棉應力松弛過程中壓力的衰減速率并非均勻，喂入量越大籽棉內部卷曲接觸點較多，籽棉間自下而上傳遞松弛壓力，彈性模量越大；隨著喂入量的增加，粘性模量明顯增大，說明隨棉纖維的體積質量和纖維間相互接觸點數的增加，纖維間抱合力變大，運動阻力增加從而粘性模量增大。

4 結 論

1）采用改進西原模型描述籽棉壓縮過程的粘彈塑特性，以曲線擬合法求解籽棉壓縮過程中非線性粘彈塑性本構模型的參數值，其決定系數2均大于0.9，結果表明模型可以很好的描述籽棉壓縮過程；采用五元件廣義Maxwell模型擬合籽棉應力松弛曲線，其決定系數2均大于0.99，結果表明五元件廣義Maxwell模型可以全面地表征籽棉的應力松弛特性。

2）籽棉壓縮試驗擬合參數結果顯示，壓縮應力及應力松弛與含水率、喂入量均呈正相關關系。通過Duncan均值比較結果可知：不同含水率組間彈性模量差異均顯著（<0.05），而粘塑體中牛頓體粘性模量3影響不顯著，不同喂入量組間串聯虎克體彈性模量1差異顯著。

3）籽棉應力松弛擬合參數結果顯示，含水率對籽棉應力松弛彈性模量、粘性模量均呈正相關影響。Duncan均值比較結果表明：不同含水率及喂入量組間的彈性模量及粘性模量差異均顯著（<0.05）。

通過研究籽棉壓縮各階段的力學模型和本構模型，將籽棉壓縮本構模型通過軟件二次開發進行仿真模擬，可縮短機具研發過程；也可為籽棉壓縮相關裝備的設計及技術參數的確定提供理論依據。

[1] 武建設，陳學庚. 新疆兵團棉花生產機械化發展現狀問題及對策[J]. 農業工程學報，2015，31(18)：5-10.

Wu Jianshe, Chen Xuegeng. Present situation, problems and countermeasures of cotton production mechanization development in Xinjiang production and construction corps[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(18): 5-10. (in Chinese with English abstract)

[2] 喻樹迅. 中國棉花產業百年發展歷程[J]. 農學學報，2018，8(1)：85-91.

Yu Shuxun. The development of cotton production in the recent hundred years of China[J]. Journal of Agriculture, 2018, 8(1): 85-91. (in Chinese with English abstract)

[3] Bennett J, Woodhouse N, Keller T, et al. Advances in cotton harvesting technology: A review and implications for the John Deere round baler cotton picker[J]. The Journal of Cotton Science, 2015, 19: 225-249.

[4] Sun Z Z, Jiang H X, Cai H P, et al. The viscoelasticity model of corn straw under the different moisture contents[J]. Mathematical Problems in Engineering, 2013, 2013: 1-5

[5] 邢獻軍，李濤，馬培勇，等. 生物質固體成型燃料熱壓成型實驗研究[J]. 太陽能學報，2016，37(10)：2660-2667.

Xing Xianjun, Li Tao, Ma Peiyong, et al. The typical biomass cold-pressing constitutive model and the influence factors analysis of viscous-elastic-plasticity[J]. Acta Energiae Solaris Sinica, 2016, 37(10): 2660-2667.

[6] 馬彥華，宣傳忠，武佩，等. 玉米秸稈振動壓縮過程的應力松弛試驗[J]. 農業工程學報，2016，32(19)：88-94.

Ma Yanhua, Xuan Chuanzhong, Wu Pei, et al. Experiment on stress relaxation of corn stover during compression with assisted vibration[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(19): 88-94. (in Chinese with English abstract)

[7] 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.

[8] 馬瑞峻，蕭金慶，鄭普峰，等. 穴盤水稻秧苗莖稈蠕變與應力松弛特性的試驗研究[J]. 農業工程學報，2018，34(13)：43-53.

Ma Ruijun, Xiao Jinqing, Zheng Pufeng, et al. Experimental study on characteristics of creep and stress relaxation for rice seedling stem raised in cell tray[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 43-53. (in Chinese with English abstract)

[9] Nona K, Lenaerts B, Kayacan E, et al. Bulk compression characteristics of straw and hay[J]. Biosystems Engineering. 2014, 118: 194-202.

[10] Zhang X, Peng W, Han L, et al. Effects of different pretreatments on compression molding of wheat straw and mechanism analysis[J]. Bioresource Technology, 2018, 251: 210-217.

[11] 吳敏，潘虹宇，趙東林，等. 燕麥莖稈理化組組分分析與動態力學特性研究[J]. 農業機械學報，2018，49（增刊）：447-455.

Wu Min, Pan Hongyu, Zhao Donglin, et al. Analysis of physicochemical composition and dynamic mechanical properties of avenasativa stalk[J]. Transaction of the Chinese Society for Agricultural Machinery, 2018, 49(Supp.): 447-455. (in Chinese with English abstract)

[12] 房佳佳，張永，楊明韶，等. 紫花苜蓿圓捆機卷壓過程的應力松弛特性與模型建立[J]. 農業工程學報，2018，34(16)：50-56.

Fang Jiajia, Zhang Yong, Yang Mingshao, et al. Stress relaxation behavior and modeling of alfalfa during rotary compression[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(16): 50-56. (in Chinese with English abstract)

[13] 朱凱，牛智有. 粉體棉花秸稈的應力松弛特性與模型建立[J]. 華中農業大學學報，2014，33(5)：130-134.

Zhu Kai, Niu Zhiyou. Compression stress relaxation characteristic of cotton stalk particles and its model building of cot[J]. Journal of Huazhong Agricultural University, 2014, 33(5): 130-134. (in Chinese with English abstract)

[14] Yang D, Li H T, Xiong Z H, et al. Mechanical properties of laminated bamboo under off-axis compression[J]. Composites Part A, 2020, 138: 106042.

[15] 劉慶庭，區穎剛，王萬章，等.甘蔗莖稈的材料模型與本構方程[C]//中國農業工程學會.中國農業工程學會2007年學術年會論文匯編，2007.

[16] 李勇，李健，吳蓓，等. 棉纖維集合體壓縮力傳遞規律分析與研究[J]. 上海紡織科技，2019，47(10)：8-11，14.

[17] 李勇，張宏，張有強，等. 棉纖維集合體壓縮力傳遞與密度關系[J]. 紡織學報，2016，37(11)：19-25.

[18] Privas E, Gawrysiak G, Lapeyre M, et al. Influence of cotton variety on compression and destructuration abilities under elevated pressure[J]. Cellulose, 2013, 20(3): 1013-1022.

[19] Hardin R, Searcy S.Viscoelastic properties of seed cotton[J]. Transactions of the ASABE, 2008, 51(3): 803-810.

[20] Hardin R, Searcy S. Modeling of seed cotton viscoelastic properties[J]. Transactions of the ASABE, 2009, 52(3): 707-714.

[21] 霍麗麗，趙立欣，田宜水，等. 生物質顆粒燃料成型的黏彈性本構模型[J]. 農業工程學報，2013，29(9)：200-206.

Huo Lili, Zhao Lixin, Tian Yishui, et al. Viscoelastic constitutive model of biomass pellet[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(9): 200-206. (in Chinese with English abstract)

[22] 陳嘯，孔丹丹，王紅英，等. 基于本構模型的顆粒飼料成型特性研究[J]. 農業工程學報，2017，33(23)：267-275.

Chen Xiao, Kong Dandan, Wang Hongying, et al. Analysis of forming properties based on pellet feed forming constitutive model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(23): 267-275. (in Chinese with English abstract)

[23] 馬方，白雪衛，劉德軍，等. 典型生物質冷壓本構模型及黏彈塑性影響因素[J]. 太陽能學報，2017，38(1)：98-105.

Ma Fang, Bai Xuewei, Liu Dejun, et al. The typical biomass cold pressing constitutive model and the influence factors analysis of viscous-elastic-plasticity[J]. Acta Energiae Solaris Sinica, 2017, 38(1): 98-105. (in Chinese with English abstract)

[24] Kaliyan N, Morey R V. Constitutive model for densification of corn stover and switchgrass[J]. Biosystems Engineering, 2009, 104(1): 47-63.

[25] 中國纖維檢驗局，江蘇省纖維檢驗所，農業部棉花纖維品質監督檢驗測試中心. HVI棉纖維物理性能試驗方法：GB/T 20392-2006[S]. 北京：中國標準出版社，2006.

[26] 朱兆祥. 材料本構關系理論講義[M].北京：科學出版社，2015：56-72.

[27] Ozturk O K, Takhar P S. Stress relaxation behavior of oat flakes[J]. Cereal Science, 2017, 77: 84-89.

[28] Hernandez Z, Figueroa J, Rayas-Duarte P, et al. Influence of high and low molecular weight glutenins on stress relaxation of wheat kernels and the relation to sedimentation and rheological properties[J]. Journal of Cereal Science, 2012, 55(3): 344-350.

[29] 閆翠珍，武國峰，王效華，等. 稻秸高密度壓縮過程中應力松弛的試驗研究[J]. 農機化研究，2015，(11)：181-187.

Yan Cuizhen, Wu Guofeng, Wang Xiaohua, et al. The study of rice straws’ stress-relaxation in high density compression[J]. Journal of Agricultural Mechanization Research, 2015, (11): 181-187.

[30] 周祖鍔.農業物料學[M]. 北京：中國農業出版社，1994：25-48.

[31] Long R, Bange M, Gordon S, et al. Fiber quality and textile performance of some Australian cotton genotypes[J]. Crop Science, 2010, 50(4): 1509-1518.

[32] Hearle J. Cotton-science and technology[M]. England: Woodhead Publishing Limited, 2007: 35-67.

Mechanical properties and construction of constitutive model for compression and stress relaxation of seed cotton

Kong Fanting, Wu Teng, Chen Changlin, Sun Yongfei, Xie Qing, Shi Lei※

(210014,)

Cotton is one of the most important cash crops and textile materials in the world. The cotton production of China is also ranked at the top of the world, in terms of planting area and total yield. Specifically, the growing area of cotton is annually more than 3 million hectares. Packing storage and transportation of seed cotton have become two main determinants of production capacity in the cotton industry. However, the traditional storage and transportation of seed cotton cannot meet the high demand for the machine-picked cotton area in recent years. The packing facilities inevitably need to fit into the specific producing conditions in China, such as local topography, climate, and vegetation. Mechanical properties of seed cotton are thus fundamental, including compression and stress relaxation, to the optimal design of equipment for harvesting, stacking, and packing. Most previous researches were focused mainly on the mechanical properties of corn, wheat, and herbage straw in agricultural materials. Few studies were reported on the mechanical properties of seed cotton, particularly on compression and stress relaxation. Taking the seed cotton as the research object in this study, a systematic test of mechanical properties was performed on the compression and stress relaxation of seed cotton. An improved Nishihara model and a generalized Maxwell 5-element mechanical model were selected to represent the stress-strain curves of compression and stress relaxation. Various level factors were utilized to verify the constitutive models for the compression and stress relaxation of seed cotton. Relevant model parameters were obtained to determine the influence rules with different factors. The results showed that the determination coefficients of parameters in the constitutive model were beyond 0.9 using the curve-fitting in the compression and stress relaxation of seed cotton. An obvious regularity of coefficients indicated that two models were better suitable for the compression and stress relaxation of seed cotton. There was also a significant influence of water content and feeding quantity on the mechanical properties of seed cotton. Furthermore, the compression stress was positively correlated to the water content and feeding quantity. In addition, the viscosity coefficient and elastic modulus increased significantly at a higher level of water content, due mainly to the porosity of cotton fiber. The cell wall of fiber became stronger and tougher, because the hydrophilic groups on the cellulose macromolecule absorbed water from the external environment. The curl of cotton fiber also made the fiber shrink longitudinally and elastic elongation during the compression under the larger amount of feeding quantity, where there was a significant increase in the overall longitudinal deformability of seed cotton and the elastic modulus. The viscosity coefficient and packing height both rose up significantly, while the pressure transmission path was longer for the naturally stacked seed cotton in the same compression chamber, as the feeding quantity was larger. The main reason was that the bending deformation, contact and extrusion were induced to generate the greater local stress between the fiber bodies during the compressing process of seed cotton. Moreover, the water content was positively correlated to the elastic modulus and viscosity coefficient of seed cotton. This positive correlation was possible because the natural twist of cotton made the fibers entangled, linked, and adhered, difficult to disperse during the stress relaxation of seed cotton. Correspondingly, the elastic modulus and viscosity coefficient increased significantly, as the feeding quantity increased. This improvement was due mainly to a non-uniform attenuation rate of pressure during the stress relaxation of seed cotton. There were much more contact points of fiber curl subjected to the greater feeding amount of seed cotton. The elastic modulus increased, while the relaxation pressure was transferred from the bottom to the top of seed cotton. The adhesion force between fibers and the viscosity coefficient rapidly rose with the motion resistance increased. This finding can provide an insightful theoretical basis to simulate the compression of seed cotton in most machine-picked cotton areas.

compression; stress relaxation; constitutive model; experimental study; seed cotton

2021-01-04

2021-03-19

江蘇省農業科技自主創新資金項目（CX(18)3046）；中央級公益性科學研究院所基本科研業務費專項（S202014）；中央級公益性科學研究院所重點任務-糧棉油收獲（15創新-2002）

孔凡婷，助理研究員，主要研究方向為大田作物機械化收獲。Email：kongfanting1989@163.com

石磊，研究員，主要研究方向為大田作物機械化收獲。 Email：shileijsnj@126.com

10.11975/j.issn.1002-6819.2021.07.007

S222.1

A

1002-6819(2021)-07-0053-07

孔凡婷，吳騰，陳長林，等. 籽棉壓縮與應力松弛力學特性及模型構建[J]. 農業工程學報，2021，37(7)：53-60. doi：10.11975/j.issn.1002-6819.2021.07.007 http://www.tcsae.org

Kong Fanting, Wu Teng, Chen Changlin, et al. Mechanical properties and construction of constitutive model for compression and stress relaxation of seed cotton[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(7): 53-60. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.07.007 http://www.tcsae.org