含水率对黑水虻生物转化猪粪有机肥黏结流动的影响.pdf

中国农业科技导报，2023，25（10）：126-136Journal of Agricultural Science and Technology含水率对黑水虻生物转化猪粪有机肥黏结流动的影响周婷，孙松林，朱海英，彭才望*（湖南农业大学机电工程学院，长沙 410128）摘要：黑水虻生物转化猪粪后获得的有机肥含水率介于40%60%之间，其黏结流动参数难以通过常规手段获取。为探究含水率对有机肥黏结流动的影响，通过堆积角测定与离散元仿真相结合的方法，从数值上量化分析不同含水率（41.21%60.52%）对有机肥颗粒流动与黏结能力的影响。结果表明，有机肥-不锈钢碰撞恢复系数、有机肥-不锈钢静摩擦系数、有机肥JKR（Johnson-Kendall-Roberts）表面能显著影响有机肥堆积角；有机肥含水率在41.21%60.52%范围内，随着含水率的增加，有机肥-不锈钢静摩擦系数先增大后减小，在含水率56.56%时达到最大值，为0.39，而JKR表面能由0.22 J m-2逐渐升高到0.47 J m-2；有机肥直剪试验表明，当含水率大于40%时，有机肥颗粒间的内摩擦角逐渐减小，但颗粒间的内聚力随JKR表面能增大而增大，有机肥颗粒间易发生团聚现象，一定程度上阻碍颗粒流动，堆积角增大。研究结果可为转移、输送有机肥的机械设备研究提供基础数据支撑。关键词：黑水虻；有机肥；含水率；黏结流动；离散元doi：10.13304/j.nykjdb.2022.0389 中图分类号：S216；X713文献标志码：A文章编号：10080864（2023）10012611Effect of Moisture Content on Bond Flows of Black Soldier Fly Larvae Biotransformation Pig Manure Organic FertilizerZHOU Ting，SUN Songlin，ZHU Haiying，PENG Caiwang*（College of Mechanical and Electrical Engineering，Hunan Agricultural University，Changsha 410128，China）Abstract：The moisture content of black soldier fly biotransformation pig manure organic fertilizer can reach 40%60%under natural condition，but the flow parameters are difficult to obtain by conventional means.In order to explore the effect of moisture content on the bonding flow of organic fertilizer，the effects of different water content（41.21%60.52%）on the flow and bonding ability of organic fertilizer particles were quantitatively analyzed by the method of accumulation angle measurement and discrete element simulation.The results showed that the organic fertilizer-steel collided restitution coefficient，organic fertilizer-steel static friction coefficient and organic fertilizer Johnson-Kendall-Roberts（JKR）surface energy significantly affected the accumulation angle of organic fertilizer.With the rising of moisture content from 41.21%to 60.52%，the organic fertilizer-steel static friction coefficient firstly decreased and then increased，especially，which reached the maximum value 0.39 when the moisture content was 56.56%，the surface energy changed from 0.22 to 0.47 J m-2.The direct shear test of organic fertilizer showed that when the moisture content was greater than 40%，the internal friction angle between the organic fertilizer particles gradually decreased，the cohesion between particles increases with the increase of JKR surface energy，organic fertilizer particles were easy to agglomerate，which would restrict the flow movement of organic fertilizer to some extent，and the repose angle increased.The study was expected to provide numerical parameter reference for the black soldier fly biotransformation pig manure organic fertilizer application and development of mechanization.Key words：black soldier fly；organic fertilizer；moisture content；bond flow；discrete element收稿日期：20220510；接受日期：20220906基金项目：湖南省大学生创新训练计划项目（S202110537054）；湖南省农机装备与技术创新项目（湘财农指 2020 0107）。联系方式：周婷 E-mail：；*通信作者 彭才望 E-mail：10 期周婷等：含水率对黑水虻生物转化猪粪有机肥黏结流动的影响畜禽养殖中，生猪养殖及其养殖废弃物的规模都比较大，生猪养殖废弃物处理面临的形势非常严峻1。相对于传统的猪粪通过堆肥、发酵或直接还田等处理方法，黑水虻生物转化猪粪可有效降低猪粪累积程度、解决恶臭问题，2周内猪粪堆积减少56%，猪粪残渣中氮、磷、钾含量分别降低55.1%、44.1%、52.5%，同时可有效降低大肠杆菌、获得可作鱼粉的高蛋白黑水虻幼虫，成为猪粪资源化与肥料化利用的重要手段2。因黑水虻生物转化猪粪工艺的差异，黑水虻生物转化猪粪后有机肥含水率在40%60%，经后期进一步的有机肥工艺处理，其具备抗植物病害的活性物质且肥效高，广泛应用于牧草、果蔬种植业，可促进农作物生长发育与代谢过程3。近年来，国内外学者针对黑水虻幼虫生长发育4、生物转化效率5、重金属迁移6、幼虫营养价值及利用7、有机肥肥效与应用3等开展了大量研究，但是黑水虻生物转化猪粪后的虫沙有机肥流动特性及其黏结能力与含水率间的变化规律鲜有研究，导致机械设备设计缺乏基础数据支撑。本团队前期针对黑水虻生物转化猪粪后形成含水率为43.6%的虫沙有机肥（简称“有机肥”）进行了接触参数标定8，并在此基础上设计了斗式取料机9、双向螺旋集料装置10，一定程度上实现了有机肥收集、提升与输送。但是，试验发现机械收集、输送不同含水率（40%60%）的有机肥过程中机械作业稳定性与可靠性降低。有机肥流动性参数（内聚力、内摩擦角、接触面的摩擦系数、表面能等）是衡量有机肥流动过程的重要指标11，但不同含水率有机肥颗粒黏结流动相关接触参数难以通过常规测试方法量化获得。目前，针对黏湿颗粒物料接触参数难以获取的问题，通过堆积角测定与离散元仿真试验相结合的方法对黏重黑土12、蚯蚓粪1314、猪粪15、一般有机肥16、土壤1718等物料进行“虚拟标定”获取物料参数，具有较强的可行性，为探究含水率对黑水虻生物转化猪粪后的有机肥黏结流动影响提供了借鉴。本研究以黑水虻生物转化猪粪后的有机肥作为研究对象，对不同含水率的有机肥进行堆积角测定、离散元仿真及直剪试验；通过 Plackett-Burman Design、爬坡、Box-Burman Design 等试验方法，建立堆积角数学模型，量化获得不同含水率的有机肥流动特性参数，综合分析含水率与堆积角之间的变化关系及其对有机肥黏结流动特性的影响，以期为黑水虻生物转化猪粪后有机肥在收集、输送、分离等不同阶段的机械化作业提供有效的理论参数。1材料与方法1.1试验材料与堆积角测定1.1.1试验材料以黑水虻幼虫生物转化猪粪后形成的有机肥为试验材料，取自湖南农业大学与湖南艾布鲁环保科技股份有限公司黑水虻科研基地。将新鲜猪粪中投放一定数量的4日龄黑水虻幼虫（料虫质量比为8 1），经过812 d的生物转化获得有机肥原料，后期通过进一步的有机肥处理工艺（水分蒸发、虫卵灭杀等），形成植物用有机肥。试 验 采 用 卤 素 水 分 测 定 仪（HC-SFY001，0.5%，上海花潮实业有限公司）测定不同处理工艺下的有机肥含水率，为40%60%。结合实际生产中黑水虻生物转化猪粪工艺条件，获得5组不同 含 水 率 有 机 肥 样 品，分 别 为 A1（41.21%）、A2（44.46%）、A3（51.37%）、A4（56.56%）、A5（60.52%）。含水率低于50%时，有机肥中大部分呈疏松状的近似球体颗粒；含水率高于50%时，有机肥中出现部分黏结团状颗粒块。经筛分测量，粒径分布范围为1.42.8 mm，平均粒径为2.0 mm。1.1.2有机肥堆积角测定本试验结合圆筒提升法测量有机肥堆积角19，试验所用钢质圆筒高120 mm、内径40 mm，钢质圆筒内填满有机肥，利用 CMT5105 型万能试验机夹住钢质圆筒并以0.03 m s-1的速度向上提升，如图 1所示，参照文献8和18，用摄像机垂直拍照，照片导入计算机CAD软件中标注堆积角，分别测量5组（A1A5）不同含水率有机肥的堆积角，每组重复3次，取平均值作为最后结果。1.2仿真模型构建含水率较高的有机肥颗粒之间或有机肥对接触部件材料之间存在较大的粘附力910。本试验选取 EDEM 软件中 Hertz-Mindlin with JKR（Johnson-Kendall-Roberts）模型作为有机肥颗粒间接触模型8，JKR 表面能可表征有机肥颗粒之间或颗粒与接触部件材料之间的黏附力。在127中 国 农 业 科 技 导 报25 卷EDEM2.7软件中建立圆筒堆积角的仿真试验模型，圆筒内径 40 mm、高120 mm，接料底板为200 mm200 mm的方形钢板。建立平均粒径为2 mm的球体模型，采用标准差为0.1 mm的标准正态分布方式生成有机肥颗粒，颗粒模型粒径变化范围设置为满足平均值，圆筒以0.03 m s-1的速度垂直底板方向匀速提升，仿真模型如图2所示。有机肥的各个待标定参数如表1所示，其他参数包括不锈钢的泊松比为0.320、剪切模量为7.91010 Pa、密度为7 865 kg m-3、重力加速度为9.81 m s-2、仿真步长时间为瑞利时间的22%、数据保存间隔时间为0.01 s、仿真时间为5.5 s。1.3试验方法1.3.1Plackett-Burman（PB）试验设计Plackett-Burman筛选试验设计中10个待标定参数（T1T10）均以1和-1形式分别代表各参数高低2个水平，注：1万能试验机；2钢质圆筒；3肥堆；4摄像机；5计算机。Note：1Universal testing machine；2Steel cylinder；3Fertilizer heap；4Video camera；5Computer.图1有机肥堆积试验Fig.1Repose angle of organic fertilizer 图2离散元仿真模型Fig.2Discrete element simulation model表1离散元仿真标定参数Table 1Parameters required for discrete element simulation参数符号Parameter notationT1T2T3T4T5T6T7T8T9T10参数Parameter泊松比Poisson ratio剪切模量Shear modulus/MPa密度Particle density/（kg m-3）有机肥-有机肥碰撞恢复系数Organic fertilizer-organic fertilizer collided restitution coefficient有机肥-有机肥静摩擦系数Organic fertilizer-organic fertilizer static friction coefficient有机肥-有机肥滚动摩擦系数Organic fertilizer-organic fertilizer rolling friction coefficient有机肥-不锈钢碰撞恢复系数Organic fertilizer-steel collided restitution coefficient有机肥-不锈钢静摩擦系数Organic fertilizer-steel static friction coefficient有机肥-不锈钢滚动摩擦系数Organic fertilizer-steelr rolling friction coefficientJKR表面能JKR surface energy/（J m-2）参数水平 Parameter level-10.101.001 600 0.400.100.050.150.100.100.0500.305.502 000 0.600.550.250.400.500.400.4010.5010.002 400 0.801.000.450.650.900.700.7512810 期周婷等：含水率对黑水虻生物转化猪粪有机肥黏结流动的影响并选择1个中心点，试验次数共25次。1.3.2最陡爬坡试验将Plackett-Burman试验得到的3个显著性参数有机肥-不锈钢碰撞恢复系数、有机肥-不锈钢静摩擦系数和有机肥 JKR表面能进行最陡爬坡试验，3个显著性参数均呈现正向效应，参数值按选定步长逐渐增加，其余对堆积角的影响不显著的参数采用中间水平值。以3个显著性参数最低值为起点进行爬坡，并与含水率为41.21%和60.52%的有机肥堆积角试验结果进行极端值验证对比。1.3.3Box-Behnken（BB）试验根据最陡爬坡试验结果，Box-Behnken Design响应面分析试验中的非显著性参数选取表1中所示中间水平值，有机肥-不锈钢碰撞恢复系数为0.250.45，有机肥-不锈钢静摩擦系数为0.250.55，JKR表面能为0.200.50 J m-2。选择1个中心点进行误差估计，共进行15次试验，包括对中心点的3次重复试验。利用Minitab软件对试验结果进行拟合，获得回归方程模型并进行分析。1.4有机肥直剪试验方法内摩擦角、抗剪强度等力学参数反映散体物料的摩擦特性13，是用于衡量物料内部颗粒克服自身阻力发生滚动、滑动、流动的重要参数。本试验通过有机肥直剪试验测量物料抗剪强度，获得物料的内摩擦角，分析有机肥含水率与流动性之间的关系。试验所用设备为河北路兴安达仪器有限公司生产的BJZ型应变控制式直剪仪。1.4.1含水率测定试验用有机肥先进行过筛处理，去除杂物后，将其在遮阴条件下自然风干，使有机肥含水率缓慢下降，不定时取样，用卤素水分测定仪（HC-SFY001，0.5%，上海花潮实业有限公司）测定，获得与物理堆积角试验时相同含水率（含水率误差为2%）的有机肥。试验采用样品盒、电子天平（UTP313-4、精度为0.01 g）等仪器对不同粒度有机肥的含水率进行测试。测试过程如下：对部分含水率较高呈现黏结团状的有机肥颗粒，先通过搅拌杆打散团状有机肥至粒径为23 mm的小颗粒。然后将有机肥放入干燥洁净的样品盒中，通过电子天平称重，质量记为 m1;将盛放有机肥的样品盒放到提前预热好的卤素水分测定仪中，采用自动干燥模式，设置干燥温度为105，干燥率达到99%时自动停止干燥过程，有机肥烘干至没有水分状态，再将其放到电子天平上称重，质量记为m2；通过2次平行试验计算获得有机肥含水率(1-m2/m1)的均值，作为其最后的含水率测定值。1.4.2直剪试验针对不同含水率有机肥，分别用环刀切取4个有机肥试样。试验前将配好的有机肥样件置于直剪仪上下金属盒之间，调整好有机肥试样与垂直加压框架的接触距离，保证其接触到位。试验开始，压力和位移为0，设置垂直加压分别为50、100、200、300 kPa，手轮转速为4 r min-1，使剪切速度控制在 0.8 mm min-1，直到有机肥试样构件破裂为止，记录应力的变化值。试验采用4个试样为1组，分别在不同的垂直法向压力（）下，施加水平剪应力进行剪切，直至有机肥试样受到剪切破坏，按照试验要求加载4种不同垂直压力（），每组试验重复3次，分别得到剪切破坏时抗剪强度（），取平均值。根据库伦定律进一步确定有机肥的抗剪强度参数，内摩擦角（）和内聚力（c）。不同垂直压力下的抗剪强度参照下面公式计算。=CR（1）式中，为有机肥的抗剪强度，kPa；C为测力计率定系数，166.7 kPa mm-1；R为测力计量表读数，精度为0.01 mm。以垂直压力（）为横坐标，以抗剪强度为纵坐标，根据试验测得点（，），绘制-关系曲线，其拟合直线的倾角为有机肥的内摩擦角，直线在纵轴上的截距为有机肥的内聚力c。泊松比计算公式如下。v=K01+K0K0=1-sin（2）式中，v为有机肥颗粒的近似泊松比，K0为侧压力系数，无量纲，根据式（2），可由内摩擦角求得有机肥的近似泊松比。2结果与分析2.1Plackett-Burman（PB）试验结果与分析不同参数在不同水平组合下测定的堆积角值如表2所示。方差分析结果（表3）表明，有机肥-不锈钢碰撞恢复系数、有机肥-不锈钢静摩擦系129中 国 农 业 科 技 导 报25 卷数、有机肥的JKR表面能显著影响有机肥堆积角。对有机肥堆积角的影响程度依次为：有机肥JKR表面能有机肥-不锈钢静摩擦系数有机肥-不锈钢碰撞恢复系数，且3个显著性参数与有机肥堆积角均呈正相关，有机肥堆积角随3个显著性参数的增大而增大。其中，JKR表面能是用于衡量有机肥颗粒间接触表面黏结能力的重要数值指标15,2122，有机肥的运动状态主要取决于颗粒之间表面的黏附能力大小。2.2最陡爬坡试验结果与分析从表4可以看出，2组堆积角仿真试验与物理试验相对误差先变小后变大，且在3号试验组的相对误差均最小。因此，最陡爬坡试验可快速逼近最优参数区域，试验最优值区间位于3号组附近，选取 3 号组试验的 3 个因素值作为中心点，2号、4号试验3个因素值分别为低、高水平进行后续响应面设计。2.3Box-Behnken（BB）试验结果与分析根据PB试验结果及最陡爬坡试验结果，重新选取各参数的水平值，进行三因素三水平中心组合试验，试验结果如表5所示。对表5试验结果进行拟合分析，获得堆积角与3个显著性参数间的表2PB筛选试验结果Table 2Pile angle result of Plackett-Burman design序号No.12345678910111213141516171819202122232425参数水平 Parameter levelT111111-11-111-1-111-1-11-11-1-1-1-1-10T2-111111-11-111-1-111-1-11-11-1-1-1-10T3-1-111111-11-111-1-111-1-11-11-1-1-10T4-1-1-111111-11-111-1-111-1-11-11-1-10T5-1-1-1-111111-11-111-1-111-1-11-11-10T61-1-1-1-111111-11-111-1-111-1-11-1-10T7-11-1-1-1-111111-11-111-1-111-1-11-10T81-11-1-1-1-111111-11-111-1-111-1-1-10T9-11-11-1-1-1-111111-11-111-1-111-1-10T10-1-11-11-1-1-1-111111-11-111-1-111-10堆积角Repose angle/（）25.7714.0754.2911.3021.6922.2914.0652.3135.8267.4850.5150.9474.6250.0623.8957.0727.6459.0156.5931.5231.9832.0351.4611.9441.8913010 期周婷等：含水率对黑水虻生物转化猪粪有机肥黏结流动的影响二次多项式回归方程如下。Y=43.68+0.29T7+4.12T8+3.14T10+0.17T7T8+0.35T7T10-1.63T8T10+0.42T72-3.32T82-1.46T210（3）该回归模型方差分析结果如表6所示，在给定的试验因素水平范围内，T8（有机肥-不锈钢静摩擦系数）、T10（有机肥JKR表面能）对堆积角影响极显著；T7（有机肥-不锈钢碰撞恢复系数）对堆积角影响不显著；T82对堆积角影响极显著，T8T10对堆积角影响显著，其余对堆积角影响不显著。模型决定系数 R2=0.976 9，与实际数据的拟合程度高；模型的P=0.001，说明该堆积角回归模型因变量与自变量之间的关系表现为极显著；失拟项表4最陡爬坡试验结果Table 4Results of steep climbing test序号No.12345有机肥-不锈钢碰撞恢复系数Organic fertilizer-steel collided restitution coefficient0.150.250.350.450.55有机肥-不锈钢静摩擦系数Organic fertilizer-steel static friction coefficient0.100.250.400.550.700.050.200.350.500.65JKR表面能JKR surface energy/（J m-2）仿真结果Simulation result/（）12.8126.2244.3667.2371.44含水率41.21%有机肥堆积角Repose angle of 41.21%moisture organic fertilizer试验值Test value/（）37.25相对误差Relative error/%65.6129.6116.0344.5947.86含水率60.52%有机肥堆积角Repose angle of 60.52%moisture organic fertilizer试验值Test value/（）45.53相对误差Relative error/%71.8642.412.5732.2836.27表3PB筛选试验方差分析Table 3Analysis of variance of Plackett-Burman design来源Resources模型Model线性LinerT1T2T3T4T5T6T7T8T9T10弯曲Bend误差Error合计Sum自由度Df1110111111111111324效应Effect-1.80-0.96-5.62-0.204.554.3510.8711.872.5226.93离均差平方和Adi SS6 406.136 396.2419.375.51189.730.25124.03113.36709.16845.3838.104 351.359.891 699.328 105.45均方Adi MS582.38639.6219.375.51189.730.25124.03113.36709.16845.3838.104 351.359.89130.72F值F value4.464.890.150.041.450.000.950.875.436.470.2933.290.08P值P value0.006*0.005*0.7070.8410.2500.9660.3480.3690.037*0.025*0.5980.000*0.788注：*和*分别表示在P0.05和P0.01水平显著。Note：*and*indicate significance at P0.05 and P0.01 levels，respectively.131中 国 农 业 科 技 导 报25 卷表5Box-Behnken Design方案及结果Table 5Scheme and results of Box-Behnken design序号NO.123456789101112131415T7：有机肥-不锈钢碰撞恢复系数Organic fertilizer-steel collided restitution coefficient0.250.450.250.450.250.450.250.450.350.350.350.350.350.350.35T8：有机肥-不锈钢静摩擦系数Organic fertilizer-steel static friction coefficient0.250.250.550.550.400.400.400.400.250.550.250.550.400.400.40T10：JKR表面能 JKR surface energy/（J m-2）0.350.350.350.350.200.200.500.500.200.200.500.500.350.350.35Y：堆积角Repose angle/（）32.8937.9442.2845.9943.5142.5938.9741.4729.6840.7133.3543.6439.9238.6142.51表6Box-Behnken Design二次回归模型方差分析Table 6ANOVA of Box-Behnken Design quadratic model来源Resources模型ModelT7T8T10T27T28T210T7T8T7T10T8T10误差Error失拟项Lack of fit纯误差Pure error合计Sum自由度Df911111111153214均方SSAdj275.1100.684135.79579.0020.64440.7597.8970.1090.50410.6936.5034.2112.291281.613离均差平方和MSAdj30.5680.684135.79579.0020.64440.7597.8970.1090.50410.6931.3011.4041.146F值F value23.500.53104.4160.750.4931.346.070.080.398.221.23P值P value0.001*0.5010.000*0.001*0.5130.003*0.0570.7840.5610.035*0.479注：*和*分别表示在P0.05和P0.01水平显著。Note：*and*indicate significance at P0.05 and P0.01 levels，respectively.13210 期周婷等：含水率对黑水虻生物转化猪粪有机肥黏结流动的影响P=0.479，表明方程拟合良好。因此，可通过堆积角结果，反推导有机肥物料的有机肥-不锈钢碰撞恢复系数、有机肥-不锈钢静摩擦系数、有机肥JKR表面能，从而进一步推导出不同含水率有机肥的接触参数。为验证模型的有效性，对不同含水率有机肥堆积角进行试验验证，结果如表7所示。由表7可知，不同含水率(41.21%60.52%)有机肥物理与仿真堆积角试验的相对误差均在5%以内，吻合结果较好。有机肥的JKR表面能随含水率的增加而不断增大，JKR表面能对堆积角影响极其显著，其主要原因是有机肥颗粒为表面非圆滑颗粒，随着含水率的增加，导致有机肥颗粒表面黏附性能提高，发生颗粒团聚现象，表现出其表面能随含水率的增加而提高；有机肥-不锈钢静摩擦系数先增大后减小，在含水率为56.56%时有机肥颗粒与不锈钢静摩擦系数达到最大值，为0.39，此后随含水率的继续增大，其静摩擦系数迅速减小，其主要原因是有机肥含水率由小增大的过程中，有机肥中越来越多的弱结合水逐渐在强结合水外围吸附，有机肥颗粒外围结合水分子数量增加，因此，加大了有机肥与不锈钢之间的静摩擦系数。但是，当含水率超过56.56%，有机肥颗粒间的弱结合水逐渐形成较厚的结合水膜，颗粒间距增大，尤其是在结合水膜厚度增大到一定程度时，可起到润滑作用，使有机肥与不锈钢之间的静摩擦系数减小。2.4含水率与有机肥堆积角关系分析有机肥的堆积角随含水率的增加而增加，与一般散体物料的研究结果相似14,23。不同含水率（41.21%60.52%）有机肥的堆积角变化如图3所示，根据实际有机肥含水率与堆积角的变化关系，含水率降低时，堆积角减小的速率呈放缓趋势。本研究分别采用多项式方程与线性方程对散点图进行拟合，得到有机肥物料的堆积角（y，）和有机肥物料的含水率（x，%）关系的数学模型如下所示。y=0.001 4x3-0.219 6x2+11.382x-160.44，R2=0.993 2（4）y=0.406 1x+20.541，R2=0.984 3（5）拟合方程（4）和（5）的拟合度都相对较高，通过对比图3中测量点误差范围与趋势，可知多项式方程（4）更符合含水率与堆积角两者之间的变化趋势。因此，可通过测量有机肥的含水率，并进一步利用式（5）计算预测得到有机肥的堆积角。2.5有机肥黏结流动分析通过直剪试验得到的不同含水率有机肥法向应力-抗剪强度（-）关系曲线如图 4 所示。在50300 kPa的法向应力下发生剪切破坏，不同含水率的有机肥-关系拟合直线的回归方程决定系数R2均在0.96以上，线性关系较好。从图4分析和式（1）计算可知，有机肥内聚力和内摩擦角随含水率的增大而减小，主要原因是受有机肥颗粒间的摩擦和咬合作用影响，含水率的增加使有机肥颗粒间的结合水膜变厚，颗粒间的间距不断增大，颗粒间的吸引力和摩擦作用减弱，其相互胶结作用受到破坏。同时，根据式（2），通过侧压力系数Ko，计算获得A1A5组有机肥的平均泊松比分别为0.33、0.35、0.38、0.40、0.44。表7模型验证试验结果Table 7Model verification test results组别GroupA1A2A3A4A5含水率Moisture content/%41.2144.4651.3756.5660.52堆积角试验值Repose angle test value/（）37.2538.6741.6242.8345.53堆积角仿真值Repose angle simulation value/（）37.8439.5342.7341.5246.81相对误差Relative error/%1.562.182.603.062.73有机肥-不锈钢碰撞恢复系数Organic fertilizer-steel collided restitution coefficient0.430.290.350.260.43有机肥-不锈钢静摩擦系数Organic fertilizer-steel static friction coefficient0.340.350.360.390.38JKR表面能JKR surface energy/（J m-2）0.220.240.320.320.47133中 国 农 业 科 技 导 报25 卷通过直剪试验得到有机肥内聚力和内摩擦角与不同含水率的关系如图5所示，有机肥内聚力随含水率的增加先减小后增大，内摩擦角随含水率的增加而减小，主要原因是有机肥内聚力主要受颗粒间相互吸引力作用影响，内摩擦角主要受颗粒间的摩擦作用影响；有机肥中水分主要起润滑与黏合作用，当含水率远小于40%时，因有机肥颗粒间的水分缺乏，颗粒间润滑作用减弱，此时内聚力较低，摩擦效应较强而导致内摩擦角较大；当含水率逼近 40%时，有机肥颗粒间可能存在“咬合关系”最强的某空间位置，此时，有机肥抗剪强度指标内聚力和内摩擦角值均较高或处于减小趋势中；随着含水率增大并大于40%时，有机肥颗粒间因水膜的联结作用增强而导致颗粒间孔隙增多，有机肥颗粒间的关联性降低，有机肥颗粒间的内摩擦角不断减小，而内聚力因颗粒间表面能的增大而增大，易发生黏结团聚，一定程度上阻碍有机肥的翻滚运动。3讨论黑水虻生物转化猪粪是畜禽废弃物资源化利用研究的热点方向之一，针对黑水虻生物转化猪粪后的有机肥含水率处于40%60%，而不同含水率有机肥的物理接触参数及其流动特性并不清图4不同含水率有机肥法向应力-抗剪强度关系Fig.4Relationship between normal stress and shear strength of organic fertilizers with different moisture content图3有机肥含水率与堆积角的关系曲线Fig.3Relationship between water content and physical accumulation angle of organic fertilizer13410 期周婷等：含水率对黑水虻生物转化猪粪有机肥黏结流动的影响晰，导致机械设备开发缺乏对应的基础数据支撑。因此，方便、准确获得不同含水率的黑水虻生物转化猪粪有机肥仿真参数及流动特性具有重要意义。与其他已有研究仅对单一含水率的物料进行参数标定相比，本文应用堆积角测定与离散元仿真试验相结合的方法，量化不同含水率有机肥流动性参数指标并进行验证分析，比较全面的分析了显著影响有机肥颗粒堆积角的因素，对比验证试验相对误差在5%以内，试验结果较为理想，为后续试验提供了基础数据支撑；建立了含水率与堆积角之间的数学模型，并采用多项式拟合模型获得了相对更高拟合精度，可及时准确地对有机肥堆积角进行预测，有利于进一步推导物料其他接触参数；本研究进一步通过有机肥直剪试验反映了有机肥含水率与内聚力、内摩擦角之间的关联，研究结果为机械收集、转移、输送不同含水率有机肥提供了有效的基础数据。总体来说，本研究以黑水虻生物转化猪粪后的有机肥为研究对象，获得了不同含水率与堆积角之间的变化关系及其对有机肥黏结流动特性的影响规律，为后续设备开发与试验提供了基础数据支撑。参 考 文 献1 陈秋红,张宽.新中国70年畜禽养殖废弃物资源化利用演进J.中国人口资源与环境,2020,30(6):166-176.CHEN Q H,ZHANG K.The evolution of resource utilization of livestock and poultry breeding waste in the past 70 years since the founding of P.R.China J.China Population,Resour.Environ.,2020,30(6):166-176.2 BESKINA K V,HOLCOMBA C D,CAMMACK