青贮玉米粉碎物料压缩特性响应面优化试验.pdf

引文格式：赵玉清,齐文涛,赵勇,等.青贮玉米粉碎物料压缩特性响应面优化试验J.云南农业大学学报(自然科学),2023,38(4):704713.DOI:10.12101/j.issn.1004-390X(n).202303039青贮玉米粉碎物料压缩特性响应面优化试验*赵玉清1,2,3，齐文涛1,3，赵勇1,3，马兴彪4，朱树早4，陈再兴1,3，熊坚2，赵伟1,3*(1.云南农业大学机电工程学院，云南昆明650201；2.昆明理工大学交通工程学院，云南昆明650093；3.云南省高原特色农业产业研究院，云南昆明650201；4.云南省寻甸县农业机械技术推广站，云南昆明655200)摘要:【目的】研究青贮玉米粉碎物料压缩特性，解决青贮玉米压缩机械作业能耗大、压缩密度低的问题。【方法】以喂入量、压缩密度和压缩速度为试验因素，以比能耗、压缩时间和最大压缩力为试验指标，开展二次回归正交试验，进行试验因素显著性分析，建立回归模型，对试验结果进行响应面优化。【结果】各试验因素对比能耗影响顺序为压缩密度喂入量压缩速度，对压缩时间影响顺序为压缩速度压缩密度喂入量，对最大压缩力影响顺序为压缩密度压缩速度喂入量。响应面优化后的最佳因素水平组合为喂入量0.90kg、压缩密度 550kg/m3、压缩速度 245.50mm/min。【结论】本研究得到的回归模型误差小，优化结果可信度高，可用于青贮玉米饲料制备工艺参数选取，也可为青贮玉米压缩装置设计提供理论参考。关键词:青贮玉米粉碎物料；压缩特性；响应面优化中图分类号:S513文献标志码:A文章编号:1004390X(2023)04070410Experiment Study on Compression Characteristics of Silage CornCrushed Material Based on Response Surface OptimizationZHAOYuqing1,2,3，QIWentao1,3，ZHAOYong1,3，MAXingbiao4，ZHUShuzao4，CHENZaixing1,3，XIONGJian2，ZHAOWei1,3(1.CollegeofMechanicalandElectricalEngineering,YunnanAgriculturalUniversity,Kunming650201,China;2.SchoolofTransportationEngineering,KunmingUniversityofScienceandTechnology,Kunming650093,China;3.YunnanPlateauCharacteristicAgriculturalIndustryResearchInstitute,Kunming650201,China;4.AgriculturalMachineryTechnologyExtensionStationofXundianCounty,Kunming655200,China)Abstract:PurposeTostudythecompressioncharacteristicsofcrushedsilagecornmaterialsandsolvetheproblemsofhighenergyconsumptionandlowcompressiondensityofsilagecorncompres-sionmachinery.MethodsUsingfeedingamount,compressiondensity,andcompressionspeedasexperimentalfactors,andspecificenergyconsumption,compressiontime,andmaximumcompres-sionforceasexperimentalindicators,aquadraticregressionorthogonalexperimentwasconductedtoanalyzethesignificanceofexperimentalfactors,establisharegressionmodel,andoptimizethere-sponsesurfaceoftheexperimentalresults.ResultsTheorderofinfluenceonspecificenergycon-云南农业大学学报（自然科学），2023，38(4)：704713http:/JournalofYunnanAgriculturalUniversity(NaturalScience)E-mail:收稿日期：2023-03-22修回日期：2023-07-01网络首发日期：2023-08-02*基金项目：云南省重大科技专项绿色食品国际合作研究中心项目(2019ZG009)；云南省教育厅科学研究基金项目(2020Y198)。作者简介：赵玉清(1975)，男，云南大理人，硕士，副教授，主要从事智能农业机械装备研究。E-mail：*通信作者 Correspondingauthor：赵伟(1983)，男，山西长治人，博士，讲师，主要从事农业收获技术与装备研究。E-mail：网络首发地址：https:/ speed,the order of influence on compression time was compression speedcompressiondensityfeedingamount,andtheorderofinfluenceonmaximumcompressionforcewascompressiondensitycompressionspeedfeedingamount.Theoptimalfactorlevelcombinationafterresponsesur-faceoptimizationwasthefeedingamountof0.90kg,compressiondensityof550kg/m3,andcom-pressionspeedof245.50mm/min.ConclusionTheregressionmodelobtainedinthisstudyhassmallerrorandhighreliabilityofoptimizationresults,whichcanbeusedforselectingprocesspara-metersforthepreparationofsilagecornfeedandalsoprovidetheoreticalreferenceforthedesignofsilagecorncompressionequipment.Keywords:crushedsilagecornmaterial;compressioncharacteristics;responsesurfaceoptimization2015 年中国“粮改饲”政策实施以来，玉米因含有丰富的糖类(73.2%)、粗蛋白(8.5%)和脂肪(4.3%)等营养物质而成为反刍动物的重要饲料来源之一，被广泛运用于农牧产业，逐步成为现代农业发展的重点农作物1-5。青贮玉米饲料是将处于乳熟末期至蜡熟前期的全株玉米经收获机械切割、揉搓、输送、压缩和打捆后密封储藏，再经厌氧发酵后制备而成的一种营养丰富的饲料6-8。粉碎物料的压缩是青贮玉米饲料制备过程中的关键环节，压缩效果直接对后续青贮玉米饲料发酵质量产生重要影响9-11。因此，青贮玉米粉碎物料压缩特性研究对青贮玉米饲料机械化工艺参数的优化和压缩装置的设计开发具有重要意义。目前，较多学者从优化玉米秸秆切割和粉碎等装备的工作参数角度出发，对玉米秸秆的力学特性开展了系列研究12-16，但对经切碎和揉搓后的青贮玉米粉碎物料开展压缩特性试验的研究相对较少。青贮玉米植株收获时的含水率达到 60%70%，其力学特性也与籽粒玉米收获后的秸秆差异较大。本研究以红单 10 号青贮玉米粉碎物料为研究对象，利用万能试验机和自制压缩装置设计3 因素 3 水平二次回归正交试验，测定青贮玉米粉碎物料压缩过程中的力学特性参数，可为后续青贮饲料制备工艺和青贮玉米压缩打捆装置的设计提供理论参考。1 材料与方法1.1试验设备及材料1.1.1试验设备主要设备为美斯特 CMT6104 力学万能试验机(图 1)和自制可观测压缩装置(图 2)。万能试验机主要工作参数为：工作电压 220V，最大压力 10kN，功率 0.4kW，精度误差范围在0.1%0.1%之间，能通过调节不同的压缩速度测定压缩状态下青贮玉米粉碎物料的力学特性数据。自制可观测压缩装置由压缩底座、压缩仓、压缩杆和透明压缩挡板组成，压缩仓尺寸为 100mm图 1 CMT6104 力学万能试验机Fig.1CMT6104mechanicaluniversaltestingmachine1234注：1.压缩底座；2.透明压缩挡板；3.压缩杆；4.压缩仓。Note:pressionbase;2.transparentcompressionbaffle;pres-sionrod;pressionbin.图 2 自制可观测压缩装置Fig.2Self-madeobservablecompressiondevice第4期赵玉清，等：青贮玉米粉碎物料压缩特性响应面优化试验705100mm350mm(长 宽 高)，压缩杆底面积为98mm98mm。试验时将压缩杆与万能试验机相连接，通过控制万能试验机升降速度，让压缩杆获得相应的升降速度从而实现压缩杆对压缩仓内青贮玉米粉碎物料的压缩，可实时对压缩仓内的压缩过程进行观察。其他试验设备包括 101-2A电热鼓风干燥箱和 LT501 电子天平称等。1.1.2试验材料试验用青贮玉米粉碎物料试样(图 3)于 2020年 9 月在云南省寻甸县青贮玉米种植基地收获，其初始密度为 312kg/m3。品种为红单 10 号，种植行距 6070cm，株距 1520cm，种植方式为单株密植，产量为 46kg/m2。按照农业机械试验条件测定方法的一般规定(GB/T52562008)中关于含水率的测定方法，测得红单 10 号青贮玉米粉碎物料平均含水率为 68.42%，满足青贮农艺技术要求。1.2二次回归正交试验设计1.2.1因素水平选取根据青贮农艺技术要求，农业纤维物料在压缩青贮时，压缩密度应大于 450kg/m3，且压缩密度对后期青贮饲料发酵品质有重要影响17-18。课题组前期对青贮玉米粉碎物料进行的单因素压缩试验结果表明：喂入量和压缩速度对压缩特性有显著影响。因此，选取喂入量(A)、压缩密度(B)和压缩速度(C)为试验因素，设计青贮玉米粉碎物料压缩特性的正交试验(表 1)。1.2.2试验指标选取比能耗是指单位质量青贮玉米粉碎物料被压缩成型时所消耗的能量，是压缩过程中评价能耗高低的重要技术指标。利用 MATLAB 软件对压缩过程中压缩力与压缩量曲线图进行梯形数值积分，梯形总面积即为比能耗15。压缩力和压缩量可由万能试验机测得，比能耗计算公式为=Wm=103mwF(x)dx。WmF(x)x式中：为压缩过程比能耗，J/kg；为压缩过程总能耗，J；为物料质量，kg；为压缩过程中随压缩量变化的压力，N；为压缩过程中的压缩量，mm。压缩时间是指压缩时青贮玉米粉碎物料状态由蓬松到压实所需的时间，能反映压缩过程的压缩效率。最大压缩力是指将青贮玉米粉碎物料压缩成型时所需的压力，能反映青贮玉米粉碎物料的可压缩程度，对压缩装置的动力选择有重要指导意义。因此，选取比能耗()、压缩时间(T)和最大压缩力(P)为试验指标。1.2.3试验方案设计响应面分析法(Box-Behnken)是一种综合试验设计和数学建模的优化方法，采用多元二次方程来拟合因素和响应值之间的函数关系，通过对回归方程的分析寻求最优工艺参数。响应面分析法能够以最少试验次数优化试验过程，通过方差分析、响应面分析和响应面优化等试验数据分析得出因素的最佳组合和响应的最优值。本研究利用 Design-Expert10.0.7 中的 Box-Behnken 模块对青贮玉米粉碎物料试样进行 3 因素 3 水平二次回归正交压缩试验(表 2)，建立回归模型以表达比能耗、压缩时间和最大压缩力3 个试验指标与喂入量、压缩密度和压缩速度3 个因素取值之间的变化规律，共 17 次试验。试验过程中，各因素水平组合重复 3 次，结果取平均值。因素间交互效应可用响应面和等高线直观表示。响应面坡度越陡，因素间交互作用越显著，反之则不显著19；等高线则用于反映因素间交互作用的强弱20。为获得压缩青贮玉米粉碎物图 3 青贮玉米粉碎物料Fig.3Crushedsilagecornmaterial表 1 二次回归正交试验因素水平表Tab.1Factorleveltableofquadraticregressionorthogonaltest水平level因素factor喂入量/kgfeedingamount压缩密度/(kgm3)compressiondensity压缩速度/(mmmin1)compressionspeed10.555010000.760020010.9650300706云南农业大学学报第38卷料过程中的最佳因素水平组合，采用 Box-Behn-ken 模块对二次回归正交试验结果进行优化处理；为验证优化后最佳因素水平组合的可信度，按照最佳因素水平组合进行 10 次重复试验。2 结果与分析2.1压缩特性试验结果由表 2 可知：青贮玉米粉碎物料压缩时，比能耗、压缩时间和最大压缩力范围分别为 27.8045.26J/kg、29.1099.00s、50.3387.91kPa。在试验过程中，当青贮玉米粉碎物料被压缩到一定程度时，会有汁液从压缩底板与压缩仓的间隙流出(图 4a)。压缩后的青贮玉米粉碎物料见图 4b。2.2各指标回归模型与方差分析利用 Design-Expert10.0.7 软件建立的回归方程分别为：=36.551.75A+7.40B+1.04C；a)b)注：a)粉碎物料汁液流出；b)压缩后的粉碎物料。Note:a)crushedmaterialjuiceflowsout;b)crushedmaterialaftercompression.图 4 压缩试验过程及压缩试验后粉碎物料情况Fig.4Compressiontestprocessandcrushedmaterialaftercompressiontest表 2 Box-Behnken 试验设计方案及结果Tab.2Box-Behnkenexperimentaldesignandresults试验序号testserialnumber试验因素testfactors比能耗/(Jkg1)specificenergyconsumption压缩时间/scompressiontime最大压缩力/kPamaximumcompressionforceABC111030.8946.1051.77211027.8042.0052.15311044.2050.7085.89411045.2649.6086.04510137.3996.9068.12610131.2092.3059.00710141.4132.3067.74810135.6530.8072.72901129.6587.4050.331001143.5399.0082.661101129.2329.1052.321201143.7933.2087.911300032.4046.9068.511400037.1946.8066.611500037.0945.9068.351600037.0246.7069.991700037.6846.8068.66第4期赵玉清，等：青贮玉米粉碎物料压缩特性响应面优化试验707T=46.531.41A+3.49B31.27C+0.75AB+0.78AC1.87BC+0.68A2+15.75C2；P=68.160.45A+16.99B+2.57C0.06AB+3.53AC+0.81BC。比能耗、压缩时间和最大压缩力的综合决定系数分别为 0.9090、0.9997 和 0.9910，校正决定系数分别为 0.8880、0.9993 和 0.9850，各回归模型 P 值均小于 0.0001，失拟项分别为 0.7247、0.1277 和 0.2738，表明试验误差较小，以上 3 个指标的回归模型与实际过程相符，均具有较高可靠性，可用于预测比能耗、压缩时间和最大压缩力为最小的最佳因素水平组合。由表 3 可知：对比能耗影响极显著的因素是压缩密度(B)，影响显著的因素是喂入量(A)，各因素对比能耗的影响顺序为压缩密度(B)喂入量(A)压缩速度(C)。由表 4 可知：对压缩时间影响极显著的因素是 A、B、C、BC 和 C2，对压缩时间影响显著的因素是 AB 和 AC，各因素对压缩时间的影响顺序为压缩速度(C)压缩密度(B)喂入量(A)。由表 5 可知：对最大压缩力影响极显著的因素是 B、C 和 AC，各因素对最大压缩力的影响顺序为压缩密度(B)压缩速度(C)喂入量(A)。2.3响应面分析压缩时间随喂入量的增加而缓慢减少，随压缩密度的增加而缓慢增加，随压缩速度的增加而表 3 比能耗回归模型方差分析Tab.3Varianceanalysisofspecificenergyconsumptionregressionmodel方差来源variancesource方差和sumofvariances自由度degreeoffreedom均方差meansquaredeviationF值F-valueP值P-value模型model471.293157.1043.260.0001*A24.43124.436.730.0223*B438.231438.23120.690.0001*C8.6318.632.380.1471残差residualerror47.20133.63失拟项misfitterm28.1693.130.660.7247纯误差pureerror19.0544.76总和sum518.4916注：“*”表示影响极显著(P0.01)，“*”表示影响显著(P0.05)；下同。Note:“*”indicatesanextremelysignificantimpact(P0.01),“*”indicatesasignificantimpact(P0.05);thesameasbelow.表 4 压缩时间回归模型方差分析Tab.4Varianceanalysisofcompressiontimeregressionmodel方差来源variancesource方差和sumofvariances自由度degreeoffreedom均方差meansquaredeviationF值F-valueP值P-value模型model9015.1181126.893018.940.0001*A15.96115.9642.760.0002*B97.30197.30260.670.0001*C7825.0117825.0120963.220.0001*AB2.2512.256.030.0396*AC2.4012.406.440.0349*BC14.06114.0637.670.0003*A21.9511.965.210.0518C21047.9011047.902807.330.0001*残差residualerror2.9980.37失拟项misfitterm2.3240.583.470.1277纯误差pureerror0.6740.17总和sum9018.1016708云南农业大学学报第38卷快速减小；喂入量和压缩密度的交互作用以及喂入量和压缩速度的交互作用均对压缩时间的影响显著(图 56)，压缩密度和压缩速度的交互作用对压缩时间的影响极显著(图 7)。当压缩速度在 100200mm/min 区间内，最大压缩力随喂入量的增加而减少；当压缩速度在200300mm/min 区间内，最大压缩力随喂入量的增加而增加、随压缩速度的增加而增加；喂入量和压缩速度的交互作用对最大压缩力的影响极显著(图 8)。2.4响应面参数优化与验证试验2.4.1响应面参数优化根据响应面分析结果，以 3 个试验指标最小为优化目标对试验结果进行响应面优化。各因素和各指标间的约束条件设置见表 6。通过分析，得到青贮玉米粉碎物料压缩时的最佳因素水平组合为喂入量 0.90kg、压缩密度 550kg/m3、压缩速度 245.50mm/min，在此因素水平组合下对应的各试验指标分别为比能耗 27.88J/kg、压缩时间 31.67s、最大压缩力 53.18kPa(图 9)。2.4.2验证试验由表 7 可知：在青贮玉米粉碎物料压缩的最佳因素水平组合条件下，比能耗、压缩时间和最大压缩力的平均值分别为 28.50J/kg、32.23s 和52.32kPa。对比优化结果，比能耗、压缩时间和最大压缩力相对误差分别为 2.09%、1.36%和5500.5204060801000.60.70.80.9 5505705906106306500.50.6喂入量/kgfeeding amount喂入量/kgfeeding amount压缩时间/scompression time压缩时间/scompression time压缩密度/(kgm3)compression density压缩密度/(kgm3)compression density0.70.80.9570590610630650图 5 喂入量与压缩密度对压缩时间的交互效应Fig.5Interactioneffectoffeedamountandcompressiondensityonthecompressiontime表 5 最大压缩力回归模型方差分析Tab.5Varianceanalysisofthemaximumcompressiveforceregressionmodel方差来源variancesource方差和sumofvariances自由度degreeoffreedom均方差meansquaredeviationF值F-valueP值P-value模型model2416.566402.76178.220.0001*A1.6311.630.720.4157B2309.6212309.621021.99缠网松紧度喂入次数；宫泽奇28研究发现：在不同含水率条件下，各因素对螺旋压缩装置总功耗的影响排序为压缩密度行进速度切断长度，压缩密度为 400kg/m3、切断长度为 30mm、含水率为65%、行进速度为 0.6m/s 时压缩装置总功耗最小；吴海辉30以黄贮玉米秸秆粉碎物料的喂入量为试验因素，对压缩室内压缩力与压缩密度之间的关系进行研究，得到了压缩力随压缩密度的增加而不断增加的结论。本研究中，各试验因素对比能耗的影响顺序为压缩密度喂入量压缩速度，对压缩时间的影响顺序为压缩速度压缩密度喂入量，对最大压缩力的影响顺序为压缩密度压缩速度喂入量。从青贮玉米饲料品质角度考虑，压缩密度越大，发酵效果越好，更有利于保证青贮饲料质量。但是，从青贮玉米压缩打捆机的设计角度而言，压缩密度越大，所需机械动力越大，结构越复杂，其使用的经济性和适用性会变差。在一定的作业条件下，通过试验研究得到相对较优的作业参数组合，对青贮玉米压缩打捆机的设计具有较大的理论意义和实际应用价值。4 结论青贮玉米粉碎物料压缩过程中，压缩密度对比能耗的影响极显著，喂入量对比能耗的影响显著；喂入量与压缩密度交互作用、喂入量与压缩速度交互作用对压缩时间的影响显著；喂入量与压缩速度交互作用对最大压缩力的影响极显著。青贮玉米粉碎物料压缩过程的最佳因素水平组合为喂入量 0.90kg、压缩密度 550kg/m3、压缩速度 245.50mm/min，对应的比能耗、压缩时间和最大压缩力分别为27.88J/kg、31.67s 和53.18kPa。表 7 优化结果验证试验Tab.7Validationtestofoptimizationresults试验号testnumber比能耗/(Jkg1)specificenergyconsumption压缩时间/scompressiontime最大压缩力/kPamaximumcompressionforce123.5932.9048.73224.9333.6350.17324.1632.1849.75433.9729.9458.86532.3230.5856.42631.9530.7557.38725.9332.2351.23829.1335.4046.82930.2729.8855.331028.3434.4948.46平均值average28.5032.2352.32712云南农业大学学报第38卷参考文献侯燕红,张一桦,周凯,等.不同青贮玉米品种农艺性状及营养价值评价J.云南农业大学学报(自然科学),2022,37(4):604.DOI:10.12101/j.issn.1004-390X(n).202201037.1叶瑞卿,薛世明,杨国荣,等.云南适种高产优质青贮玉米品种筛选试验研究J.云南农业大学学报(自然科学),2012,27(4):467.DOI:10.3969/j.issn.1004-390X(n).2012.04.002.2宗亚倩,鲁洪智,段新慧,等.青贮玉米中乳酸菌的分离和鉴定J.云南农业大学学报(自然科学),2021,36(6):1071.DOI:10.12101/j.issn.1004-390X(n).202009016.3耿端阳,孙延武,牟孝栋,等.基于差速辊的青贮玉米籽粒破碎仿真试验及优化J.吉林大学学报(工学版),2022,52(3):693.DOI:10.13229/ki.jdxbgxb20200812.4王国伟,朱庆辉,于海业,等.基于数字化农机装备的青贮饲料可追溯系统J.吉林大学学报(工学版),2022,52(1):242.DOI:10.13229/ki.jdxbgxb20210392.5段珍,张红梅,张建华,等.浅谈青贮玉米品种的选择J.饲料研究,2017,40(2):4.DOI:10.13557/ki.issn1002-2813.2017.09.002.6刘组钊,刘杰,何静.青贮玉米的概念分类及高产栽培技术J.农业与技术,2019,39(2):98.DOI:10.19754/j.nyyjs.20190130055.7张钰,田海清,赵志宇,等.基于改进离散粒子群算法的青贮玉米原料含水率高光谱检测J.农业工程学报,2019,35(1):285.DOI:10.11975/j.issn.1002-6819.2019.01.035.8韩战强,宋艳画,王志方.影响青贮玉米品质因素的研究进展J.饲料研究,2020,43(1):106.DOI:10.13557/ki.issn1002-2813.2020.01.029.9王帅,张彬,李显旺,等.玉米鲜秸秆力学特性试验研究J.农机化研究,2018,40(3):171.DOI:10.13427/ki.njyi.2018.03.034.10刘富强,李景彬,坎杂,等.粉碎玉米秸秆压缩特性试验研究J.农机化研究,2019,41(1):169.DOI:10.13427/ki.njyi.2019.01.030.11张锋伟,宋学锋,张雪坤,等.玉米秸秆揉丝破碎过程力学特性仿真与试验J.农业工程学报,2019,35(9):58.DOI:10.11975/j.issn.1002-6819.2019.09.007.12杜晓雪.青贮甜高粱秸秆压缩及应力松弛特性试验研究D.呼和浩特:内蒙古农业大学,2019.13陈美舟,徐广飞,宋志才,等.全株青贮玉米径向局部压缩特性试验J.吉林大学学报(工学版),2021,6(1):1.DOI:10.13229/ki.jdxbgxb20201008.14邢芳.玉米秸秆物理机械特性试验研究D.沈阳:沈阳农业大学,2017.15朱振伟,张开飞,李赫,等.玉米秸秆力学特性的研究J.江苏农业科学,2017,45(12):178.DOI:10.15889/j.issn.1002-1302.2017.12.047.16王旭哲,张凡凡,马春晖,等.压实度对玉米青贮开窖后营养品质及有氧稳定性的影响J.农业工程学报,2018,34(6):300.DOI:10.11975/j.issn.1002-6819.2018.06.038.17盛凯,冉毅,艾平,等.不同压缩比和青贮时间下玉米秸秆厌氧发酵性能研究J.农业工程学报,2020,51(7):338.DOI:10.6041/j.issn.1000-1298.2020.07.038.18王秀文,柴智,王颖莉.响应面优化复方黄芪饲料添加提取工艺与体外抗氧化作用研究J.饲料研究,2021,44(3):53.DOI:10.13557/ki.issn1002-2813.2021.03.013.19王亚玲,陶然,梁峰源,等.响应面法优化中药降脂复方多糖提取工艺及抗氧化研究J.当代化工,2020,49(3):555.DOI:10.13840/21-1457/tq.2020.03.012.20李艳,李彤,张丹娜,等.发酵饲料在畜牧养殖业中的应用J.中国畜禽种业,2021,17(12):86.DOI:10.3969/j.issn.1673-4556.2021.12.040.21李心平,熊师,耿令新,等.含水率对玉米果穗抗压特性的影响J.农业工程学报,2018,34(2):25.DOI:10.11975/j.issn.1002-6819.2018.02.004.22刘立意,邴媛媛,董欣,等.玉米秸秆压缩加载条件对其裂纹和机器功耗的影响J.农业工程学报,2014,30(3):163.DOI:10.3969/j.issn.1002-6819.2014.03.022.23刘禹辰,张锋伟,宋学锋,等.基于离散元法玉米秸秆双层粘结模型力学特性研究J.东北农业大学学报,2022,53(1):45.DOI:10.19720/ki.issn.1005-9369.2022.01.006.24赵玉清,赵勇,杨志,等.青贮玉米秸秆径向和轴向剪切特性试验研究J.饲料研究,2021,44(17):93.DOI:10.13557/ki.issn1002-2813.2021.17.022.25赵玉清,方佳梦,杨志,等.青贮玉米植株剪切力学特性试验研究J.饲料研究,2021,44(3):93.DOI:10.13557/ki.issn1002-2813.2021.03.022.26刘晓东.青贮玉米秸秆圆捆应力松弛特性试验研究D.呼和浩特:内蒙古农业大学,2020.27宫泽奇.青贮玉米秸秆螺旋致密成型工艺及压缩过程流变学试验研究D.北京:中国农业大学,2017.28马斌,韩少平,贺成柱,等.不同含水率苜蓿干草的高密压缩成型特性的试验研究J.农机化研究,2022,44(5):152.DOI:10.13427/ki.njyi.2022.05.028.29吴海辉.玉米秸秆高密度压缩过程中的压缩及应力松弛特性研究D.呼和浩特:内蒙古农业大学,2010.30责任编辑：何謦成第4期赵玉清，等：青贮玉米粉碎物料压缩特性响应面优化试验713