• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

鲜莲籽离散元仿真参数标定

张国忠, 张清洪, 焦俊, 陈亚欣, 梁升, 刘浩蓬

张国忠, 张清洪, 焦俊, 等. 鲜莲籽离散元仿真参数标定[J]. 华南农业大学学报, 2023, 44(1): 170-178. DOI: 10.7671/j.issn.1001-411X.202204007
引用本文: 张国忠, 张清洪, 焦俊, 等. 鲜莲籽离散元仿真参数标定[J]. 华南农业大学学报, 2023, 44(1): 170-178. DOI: 10.7671/j.issn.1001-411X.202204007
ZHANG Guozhong, ZHANG Qinghong, JIAO Jun, et al. Parameters calibration of discrete element simulation for fresh lotus seeds[J]. Journal of South China Agricultural University, 2023, 44(1): 170-178. DOI: 10.7671/j.issn.1001-411X.202204007
Citation: ZHANG Guozhong, ZHANG Qinghong, JIAO Jun, et al. Parameters calibration of discrete element simulation for fresh lotus seeds[J]. Journal of South China Agricultural University, 2023, 44(1): 170-178. DOI: 10.7671/j.issn.1001-411X.202204007

鲜莲籽离散元仿真参数标定

基金项目: 国家特色蔬菜产业技术体系专项(CARS-24-D-02)
详细信息
    作者简介:

    张国忠,教授,博士,主要从事农业机械化与自动化研究,E-mail: zhanggz@mail.hzau.edu.cn

    通讯作者:

    刘浩蓬,工程师,主要从事农业机械化与自动化研究,E-mail: liuhaopeng@mail.hzau.edu.cn

  • 中图分类号: S226

Parameters calibration of discrete element simulation for fresh lotus seeds

  • 摘要:
    目的 

    确定鲜莲籽机械化加工过程中离散元仿真模型参数,为鲜莲籽机械化加工仿真试验提供数据参考。

    方法 

    本研究利用EDEM仿真软件开展鲜莲籽离散元仿真参数标定。以产自湖北洪湖的‘太空莲36号’为试验对象,通过落种试验测定鲜莲籽实际落种的堆积角和休止角。基于Hertz-Mindlin (no slip)接触模型进行鲜莲籽落种仿真试验,以鲜莲籽堆积角和休止角的实测值与仿真值之间的误差为试验指标,通过Plackett-Burman试验确定对堆积角和休止角影响显著的接触参数,通过最陡爬坡试验确定鲜莲籽离散元模型最优接触参数组合。采用料斗进行实际落种验证试验,以莲籽落种速率为试验指标,对比实际与仿真落种验证试验莲籽落种速率,验证最优参数组合可靠性。

    结果 

    莲籽间静摩擦系数、莲籽间滚动摩擦系数对堆积角影响极显著(P<0.01);莲籽间滚动摩擦系数对休止角影响极显著(P<0.01),莲籽间静摩擦系数、莲籽−有机玻璃静摩擦系数对休止角影响显著(P<0.05)。最优接触参数组合为莲籽间静摩擦系数0.4、莲籽间滚动摩擦系数0.02、莲籽−有机玻璃静摩擦系数0.4。落料验证试验结果表明,实际试验与仿真试验的鲜莲籽落种速率最大相对误差不超过3.65%。

    结论 

    标定的鲜莲籽离散元仿真模型接触参数准确可靠,研究结果可为莲籽加工机械的结构优化设计提供数据参考。

    Abstract:
    Objective 

    To determine the parameters of discrete element simulation model in the mechanized processing of fresh lotus seeds, and provide data references for the mechanized processing simulation test of fresh lotus seeds.

    Method 

    The calibration of discrete element simulation parameters of fresh lotus seeds was carried out by EDEM simulation software. The accumulation angle and repose angle of the actual fresh lotus seeds were measured by seed drop test with ‘Space lotus 36’ from Honghu, Hubei. Based on the Hertz-Mindlin (no slip) contact model, a simulation test of fresh lotus seed drop was conducted, and the error between the measured and simulated values of fresh lotus seed accumulation angle and repose angle was used as the test index to determine the contact parameters with significant effects on accumulation angle and repose angle through the Plackett-Burman test. The steepest climb test was conducted to determine the optimal contact parameter combinations in discrete element model for fresh lotus seeds. The actual seed drop verification test was carried out using the hopper with seed drop rate as the test index. The seed drop rates in actual and simulated seed drop verification tests were compared to verify the reliability of the optimal parameter combination.

    Result 

    The effects of static friction coefficient between lotus seeds and rolling friction coefficient between lotus seeds on the accumulation angle were highly significant (P<0.01); The effects of rolling friction coefficient between lotus seeds on the repose angle were highly significant (P<0.01), and the effects of static friction coefficient between lotus seeds and static friction coefficient between lotus seeds and plexiglass on the repose angle were significant (P<0.05). The optimal combination of contact parameters was 0.4 for static friction coefficient between lotus seeds, 0.02 for rolling friction coefficient between lotus seeds, and 0.4 for static friction coefficient between lotus seeds and plexiglass. The results of the drop verification test showed that the maximum relative error of fresh lotus seed drop rate between the actual test and the simulation test was not exceeding 3.65%.

    Conclusion 

    The contact parameters of the calibrated discrete element simulation model for fresh lotus seeds are accurate and reliable, and the findings can provide data references for the structural design optimization of lotus seed processing machinery.

  • 滴灌作为有效的节水灌溉和一体化管理方式在很多地区得到应用。滴灌设计时,为了保证灌水质量和提高水分利用率,要求滴灌毛管网中滴头出水相对均匀,而水流动过程中受毛管边壁及滴头阻力影响,毛管压力沿水流方向逐渐衰减。如何在低压条件下最大化地减少管网水流水头损失及提高灌水均匀度成为滴灌设计的关键[1-3]。滴灌毛管水头损失包括滴头段的局部水头损失和非滴头段的沿程水头损失。针对滴灌毛管沿程水头损失的研究已经较为成熟,主要基于Darcy-Weisbach公式确定[4]。对于滴头局部水头损失,依据滴灌设计规范,在缺乏资料情况下,滴头局部水头损失按毛管沿程水头损失的10%~20%计算[5]。不少学者针对滴头局部水头损失做了相关研究,虽然取得了一些成果,但是仍未给出被行业一致认可的计算模型[6]

    毛管局部水头损失因其滴头种类、尺寸差异较大,目前仍没有形成统一认识。Bagarello等[7]和Provenzano等[8]对毛管局部水头损失做了相关研究,通过回归试验结果得出了滴头局部水头损失系数与滴头相对形状有关的结论,并给出计算表达式。仇振杰等[9]通过试验研究了滴头插入式毛管水头损失规律,认为毛管局部水头损失既与滴头形状有关又与雷诺数有关,但由于只是对特定尺寸毛管进行了研究,其公式适用范围受限。上述研究都是通过试验回归得出的经验公式。此外,Demir等[10]、Ferro[11]和Yildirim[12]利用量纲分析法对毛管综合摩擦阻力做了相关研究,认为试验与量纲分析理论结合可得到更合理、准确的结果。Provenzano等[13]通过计算流体动力学(Computational fluid dynamics, CFD)模拟对内镶式毛管流速分布进行相关研究,认为CFD具有很高的准确性,且节省时间和成本。石喜等[14]通过CFD数值模拟与试验研究了三通管局部水头损失,认为模拟在研究边壁形状不规则、存在流体分离等强旋流时有较高的精度。

    以上关于毛管局部水头损失的研究,大多是根据试验得到经验公式,在实际应用时计算精度仍有待继续提高。本文选用2种内镶式及1种外包式滴灌毛管,并以CFD数值模拟为基础,采用试验与模拟相结合的方法,对圆柱式滴灌毛管水力特性进行了研究。

    试验在三峡大学水工厅进行,试验装置主要包括水箱、离心泵、变频柜、高精度压力表、游标卡尺、内镶式滴灌毛管、外包式滴灌毛管、电子秤等,如图1所示。试验选用外径分别为12和16 mm的2种内镶式滴灌毛管和一种外径为16 mm的外包式滴灌毛管,利用游标卡尺测量滴头尺寸、滴头间距等,精度为0.01 mm,利用称重法测定毛管流量,精度为0.01 g。

    图  1  试验装置示意图
    1:泵;2:水箱;3:进水管;4:回水管;5:阀门;6:变频稳压器;7:电脑;8:压力表;9:回水箱;10:毛管;11:滴头;12:外包式滴头;13:内镶式滴头
    Figure  1.  Schematic diagram of experimental device
    1: Pump; 2: Water box; 3: Water inlet pipe; 4: Return pipe; 5: Valve; 6: Variable frequency stabilizer; 7: Computer; 8: Pressure gauge; 9: Return tank; 10: Trickle lateral; 11: Emitter; 12: Out-line emitter; 13: In-line emitter

    试验首先将滴头封堵,在毛管首末两端安装精度为0.1 kPa的压力计,通过首末端阀门控制毛管流量,测量周期为100 s,首部压力范围为0~20 m(压力单位用水头高度表示)。滴头间距为0.33 m,毛管长度为8.42 m,毛管详细尺寸参数见表1。温度计测得水体温度为18~23 ℃。根据所测毛管流量计算毛管流速及雷诺数,毛管首末端压力差即为总水头损失。内镶式滴灌毛管采用秦川节水公司生产的12和16 mm毛管,外包式滴灌毛管采用三峡大学自研的16 mm毛管。外包圆柱式滴头主要由滴头入口、外包迷宫流道、滴头外壁、滴头出口组成。滴头流道采用矩齿形迷宫流道,迷宫流道螺旋环绕于毛管外壁。

    表  1  滴灌毛管几何参数
    Table  1.  Geometric parameters of drip irrigation pipes
    编号
    No.
    毛管类型
    Type of trickle lateral
    毛管外径/mm
    External diameter of
    trickle lateral
    毛管内径/mm
    Internal diameter of
    trickle lateral
    滴头内径/mm
    Emitter internal
    diameter
    滴头长度/mm
    Emitter length
    滴头间距/mm
    Emitter spacing
    E1 内镶式 In-line emitter 12 10.68 8.50 40 260
    E2 内镶式 In-line emitter 16 14.34 11.60 40 260
    E3 外包式 Out-line emitter 16 14.34 14.34 40 260
    下载: 导出CSV 
    | 显示表格

    通过绘图软件CAD建立毛管滴头几何模型,采用GAMBIT构建网格模型,将网格文件导入FLUENT进行模拟计算,进口边界采用速度进口边界,出口边界采用压力边界条件,以实测温度均值20 ℃为模拟水体温度,对应运动黏滞系数为1.01×10−6 m/s2,毛管采用混合网格形式,非滴头段网格单元以六面体网格元素进行划分,滴头段网格单元采用四面体网格形式,增加局部形变处网格个数和离散程度,网格数范围为1.2×105~2.8×105。模拟计算采用Realizablek-ε模型[15],方程如下:

    $$\frac{{\partial (\rho k{u_i})}}{{\partial {x_i}}} = \frac{\partial }{{\partial {x_j}}}\left( {{\alpha _k}{\mu _{{\rm{eff}}}}\frac{{\partial k}}{{\partial {x_j}}}} \right) + {G_k} + \rho \varepsilon ,$$ (1)
    $$ \begin{split}&\;\\ \frac{{\partial (\rho \varepsilon {u_i})}}{{\partial {x_i}}} = &\frac{\partial }{{\partial {x_j}}}\left( {{\alpha _\varepsilon }{\mu _{{\rm{eff}}}}\frac{{\partial \varepsilon }}{{\partial {x_j}}}} \right) +\\ & \frac{\varepsilon }{k}(C_{1\varepsilon }^ * {G_k} - {C_{2\varepsilon }}\rho \varepsilon ), \end{split} $$ (2)

    式中:k为湍动能;ε为耗散速率;ui为流速在i方向的分量,m/s;ρ为液体的密度,kg/m3μeff为修正后流体的黏滞系数,Pa∙s。C1ε*、C2ε为经验常数;Gk为平均速度梯度引起的湍动能产生项;修正系数αk=αε=1.39;xixj表示ij方向符号。

    滴灌毛管出口与大气相接,出口条件选择自由出流,采用标准壁面函数法处理近壁区流动,使用SIMPLE算法进行求解,外包式滴头实体模型和网格模型如图2所示。

    图  2  外包式滴头实体和网格模型
    Figure  2.  Entity and network models of trickle laterals equipped with integrated out-line emitters

    根据表1中的几何参数,分别构建3组滴灌毛管几何模型及网格模型,通过Fluent对模型数值模拟。图3列出了3组滴灌毛管总水头损失实测值与模拟值对比关系。从图3中可以看出,随流速增大,毛管总水头损失逐渐增大;毛管总水头损失实测值与模拟值略有偏差,但趋势基本吻合。E1内镶式滴灌毛管实测值与模拟值最小相对误差为6.9%,最大相对误差为12.8%;E2内镶式滴灌毛管实测值与模拟值最小相对误差为5.9%,最大相对误差为17.5%。E3外包式滴灌毛管实测值与模拟值最小相对误差为3.1%,最大相对误差为7.2%。通过实测值与模拟值对比结果可以看出,2种内镶式滴灌毛管模拟值略低于实测值,主要由于试验中毛管入口、毛管出口与支管通过管道直接相连,水流经过此处存在突然扩大和突然收缩的局部水头损失,而Gambit几何模型没有考虑此段管径变化。外包式滴灌毛管模拟值与实测值相对误差更小,分析原因认为实测内镶式滴灌毛管圆柱滴头存在制造偏差,导致滴头圆柱空腔与标准圆形存在偏差,而此偏差对毛管局部水头损失影响很大,外包式滴灌毛管圆柱滴头镶嵌于毛管外部,滴头偏差不会对毛管水力特性造成影响,因此外包式滴灌毛管总水头损失模拟值与实测值更吻合。与试验相比,数值模拟具有操作方便、省工省时,具有可避免试验中制造偏差、读数误差等不可控因素的影响等优点。

    图  3  滴灌毛管水头损失实测值与模拟值
    Figure  3.  Measured and simulated head loss of drip irrigation pipes

    为比较外包式滴灌毛管与内镶式滴灌毛管总水头损失差异,按照标准灌溉小区毛管铺设长度60 m计算,设置5组滴头流量:1、3、5、7和9 L/h,并假定滴头流量均匀一致,对内镶式滴灌毛管E2和外包式滴灌毛管E3进行模拟分析。图4列出了2组毛管总水头损失对比结果。从图4中可以看出,内镶式滴灌毛管总水头损失明显大于外包式滴灌毛管,3、5、7和9 L/h流量下,外包式滴灌毛管只有内镶式滴灌毛管总水头损失的0.454、0.423、0.402和0.387;1 L/h流量下外包式滴灌毛管与内镶式滴灌毛管总水头损失的比值为0.508。图4还显示了外包式滴灌毛管相对内镶式滴灌毛管减阻比例变化关系,由图4可以看出,相对于内镶式滴灌毛管,外包式滴灌毛管减阻比例随滴头流量的增大不断增大,在灌溉流量较大时,外包式滴灌毛管相对内镶式滴灌毛管水头损失更小,减阻优势更明显,外包式滴灌毛管相较内镶式滴灌毛管减阻比例最小值为49.2%,最大值为61.3%。

    图  4  内镶式与外包式滴灌毛管总水头损失比较
    Figure  4.  Comparison of total head loss of trickle laterals equipped with integrated in-line and out-line emitters

    为比较外包式与内镶式滴灌毛管最大铺设长度,以内镶式滴灌毛管E2和外包式滴灌毛管E3为研究对象,设置首部初始工作压力水头20 m,假定滴头在5~20 m压力水头范围内流量恒定,毛管末端的压力水头刚好为5 m,设置4组流量,分别为1、3、5和7 L/h,沿程水头损失采用Darcy-Weisbach公式和Blasius公式计算[16],公式如下:

    $$ {h_f} = \lambda \frac{l}{d} \cdot \frac{{{v^2}}}{{2g}}, $$ (3)
    $$ \lambda = \frac{{64}}{{{Re}}}, $$ (4)
    $$ \lambda = \frac{{0.3164}}{{Re^{0.25}}}, $$ (5)

    式中:hf为毛管沿程水头损失,m;λ为沿程阻力系数;d为毛管内径,m;l为滴头长度,m;v为毛管水流速度,m/s;g为重力加速度,m/s2Re为毛管水流雷诺数。在Re<2 000的层流范围内,沿程阻力系数λ主要采用式(4)计算;在3000<Re<100000的光滑紊流范围内,λ采用式(5)计算;在2 000<Re<3 000的过渡区,将沿程阻力系数λ视为常数0.04[17]。局部水头损失 ${h_j} $ 的计算采用内镶式滴灌毛管局部水头损失系数法,计算表达式采用式(6)和式(7):

    $$ {h_j} = \alpha \frac{{{v^2}}}{{2g}}, $$ (6)
    $$ \alpha = 0.056\left[ {{{\left( {\frac{{{D_i}}}{{{D_g}}}} \right)}^{17.83}} - 1} \right], $$ (7)

    式中:α为局部水头损失系数;Dg为圆柱式滴头内径,m;Di为毛管内径,m。

    内镶式与外包式滴灌毛管最大铺设长度如图5所示。由图5可知,在初始压力水头为20 m,滴头设计流量为1 L/h时,内镶式滴灌毛管最大铺设长度为178 m,相同首部压力水头、毛管内径和设计流量条件下,外包式滴灌毛管极限铺设长度为240 m,极限铺设长度增加了34.8%。当设计流量为3 L/h时,内镶式滴灌毛管最大铺设长度为79 m,同等条件下,外包式滴灌毛管极限铺设长度为124 m,极限铺设长度增加了56.9%。当设计流量为5 L/h时,内镶式滴灌毛管最大铺设长度为63 m,同等条件下,外包式滴灌毛管极限铺设长度为92 m,极限铺设长度增加了46.0%。当设计流量为7 L/h时,内镶式滴灌毛管最大铺设长度为51 m,同等条件下,外包式滴灌毛管极限铺设长度为62 m,极限铺设长度增加了21.6%。结果表明:相同条件下,外包式滴灌毛管相较内镶式滴灌毛管可增大铺设长度、增加灌溉面积。

    图  5  各流量下的内镶式与外包式滴灌毛管最大铺设长度
    Figure  5.  Maximum laying length of trickle laterals equipped with integrated in-line emitters and out-line emitters at each flow rate

    图5可以看出,随着流量的增加,内镶式和外包式滴灌毛管极限铺设长度逐渐减少,同时,外包式滴灌毛管极限铺设长度明显优于内镶式滴灌毛管,外包式滴灌毛管在满足灌水均匀性的同时,能够有效地减少滴头局部水头损失,更适合较长距离的灌溉小区。相较于内镶式滴灌毛管,外包式滴灌毛管能有效地提高铺设长度,在相同铺设长度条件下,外包式滴灌毛管灌水均匀度提高。

    根据试验,外包式滴灌毛管压力−流量关系为q = 1.86 H0.54,其中,q代表滴头流量,H代表压力。图6给出了外包式滴头流道流速及压力分布规律。压力最大值出现在流道入口处,沿水流方向流道压力逐渐减小。低压区主要集中在直角弯道凸面,由于直角转弯导致边壁处水流边界层分离出现低压回旋区。高压区主要分布在弯道凹面拐角处,由于流道凹面拐角水流发生滞留回旋,流速相对较低,由能量守恒原理可知,此处压力较高,由图6a可得到验证。流速较大处主要发生在流道主流区,这与部分学者研究结果一致[18-19]。在流道直角弯折处流速较低,边壁直角弯折引起水流流向突变,水流以涡流回旋的形式存在于折角处,导致折角近壁流速降低。

    图  6  滴头流道流速和压力分布
    Figure  6.  Velocity and pressure distribution in emitter channel

    图7a可知,外包式滴灌毛管过水断面沿程均匀光滑,毛管至滴头过渡处不存在断面突变,水流沿程平顺,因此外包式滴灌毛管局部水头损失很小,压力衰减主要由管道沿程阻力引起。最大流速发生在管轴线位置,主流区流道较大,靠近边壁流速逐级减小。内镶式滴灌毛管由于滴头镶嵌于毛管内部,滴头厚度引起毛管断面形状突变,导致滴灌毛管出现局部阻力,产生局部水头损失,由图7a流速分布可知,内镶圆柱滴头处断面平均流速明显大于非滴头的毛管断面,水流从毛管到滴头段以及从滴头段到毛管段,流速大小和方向发生激变,流速在滴头段重新分布,此过程消耗大量水流能量,产生局部阻力,引起局部水头损失。由图7b可知,外包式滴灌毛管压力梯度基本均匀变化且压力衰减较慢,滴头段压力未发生急剧变化;内镶式滴灌毛管在毛管至滴头过渡段,压力梯度急剧变化,压力值迅速减小,滴头至毛管段压力逐渐增大,在滴头下游恢复至均匀衰减。

    图  7  滴灌毛管流速和压力分布
    Figure  7.  Velocity and pressure distribution in drip irrigation pipes

    数值模拟因操作方便、省工省时、可精确表达管网流场特性等优点在滴灌领域得到广泛应用。数值模拟在构建滴头外包式滴灌毛管几何模型时,通常采用标准几何体,而实际滴灌毛管存在制造偏差和运行偏差,导致实际管道几何形状与数值模拟建立的几何模型存在一定偏差,引起模拟结果与实际产生一定差异[4]。滴灌毛管总水头损失包括毛管非滴头段沿程水头损失和滴头段局部水头损失,滴头内镶改变毛管原断面形状,产生额外局部水头损失,毛管断面变化规律对滴灌毛管总水头损失影响显著,这与部分学者研究结果一致[20-22]。在流量一定条件下,毛管沿程水头损失主要与毛管管径和相对形状有关,对于首部运行压力水头大于10 m的滴灌毛管,沿程管径可视为标准圆形,而对于低压灌溉毛管沿程管径并非标准圆形,应根据实际形状进行换算[23-24]。在管径和管长一定情况下,毛管沿程水头损失通常无法减小,但通过改变滴头形状或改变镶嵌方式可降低滴头局部水头损失,本文借助CFD数值模拟滴灌毛管压力梯度衰减和精确显示流场特性的优势,研究了滴头结构形式对滴灌毛管阻力变化的影响。

    滴头流态指数是衡量滴头水力性能的一个重要指标,流态指数越小滴头流量受压力变化影响越小,而流态指数主要与流道水流紊动特性有关,流道制紊能力越强水流在流道中消耗能量越大,滴头流量受压力影响越小。因此流道结构形式对滴头水力性能影响很大[25-26]。在流道宽度和深度固定的前提下,流道锯齿形状、间距和高度偏差等是影响流道水力性能的重要因素,主要可通过增大主流区长度、增加主航道弯折频率和弯折幅度的形式增加流道水流沿程和局部水头损失,实现压力补偿,提高滴头水力性能。矩齿形迷宫流道最大流速主要发生在流道主流区,这与部分学者研究结果一致[27]。在流道直角弯折处流速较低,边壁直角弯折引起水流流向突变,水流以涡流回旋的形式存在于折角处,导致折角近壁流速降低。魏正英等[28]提出最大流速发生在流道主航道内,低速区主要发生在主航道以外,这与本文模拟结果相符,通过优化流道形状可降低阻滞区范围,增大流道过流能力和抗堵能力。

    对内镶式和外包式滴灌毛管进行灌水试验,同时采用CFD方法对不同滴头流量下内镶式及外包式滴灌毛管水力特性模拟计算,结果如下:

    1)外包式滴灌毛管及内镶式滴灌毛管总水头损失模拟值略低于实测值,但变化趋势基本吻合。外包式滴灌毛管模拟值与实测值相对误差更小。内镶式滴灌毛管实测值与模拟值最小相对误差为5.9%,最大相对误差为17.5%,外包式滴灌毛管实测值与模拟值最小相对误差为3.1%,最大相对误差为7.2%。

    2)外包式滴灌毛管相较于内镶式滴灌毛管总水头损失更小,总水头损失减小比例随滴头流量增加而增大,减阻比例大于49.2%。相同条件下,外包式滴灌毛管最大铺设长度比内镶式滴灌毛管更长,试验所选工况下,最大铺设长度增加范围为21.6%~56.9%。内镶式滴灌毛管滴头处,流速重分布引起压力梯度急变,消耗水流能量,这是内镶式滴灌毛管总水头损失增大的主要原因。

    3)外包式滴灌毛管相较于内镶式滴灌毛管滴头段水流平顺,滴灌毛管局部水头损失更小,在允许流量偏差范围内,可提高滴灌毛管铺设长度,更适合较长距离的灌溉小区,在相同铺设长度条件下,外包式滴灌毛管灌水均匀度更高。

  • 图  1   鲜莲籽短轴压缩试验装备

    Figure  1.   Equipment for short axis compression of fresh lotus seeds

    图  2   实际落种试验

    1:台架,2:种箱,3:鲜莲籽,4:挡板,5:落种盘,a1:台架长度,a2:种箱长度,a3:落种盘直径,a4:挡板长度,b1:台架宽度,b2:种箱宽度、挡板宽度,h1:台架高度,h2:种箱高度,h3:落种盘高度,h4:挡板到落种盘底面的距离,θ:休止角,φ:堆积角

    Figure  2.   Actual drop test

    1: Bench, 2: Seed box, 3: Fresh lotus seeds, 4: Baffle, 5: Seed drop tray, a1: Length of benchtop, a2: Length of seed box, a3: Diameter of seed drop tray, a4: Length of bench, b1: Width of bench, b2: Width of seed box and baffle, h1: Height of bench, h2: Height of seed box, h3: Height of seed drop tray, h4: Distance from baffle to bottom of seed drop tray, θ: Repose angle, φ: Accumulation angle

    图  3   种堆图像处理

    Figure  3.   Processing of seed pile image

    图  4   鲜莲籽离散元模型

    Figure  4.   Discrete element model of fresh lotus seed

    图  5   落种仿真试验

    Figure  5.   Simulated seed dropping test

    图  6   验证试验

    Figure  6.   Validation test

    图  7   实测与仿真结果对比

    Figure  7.   Comparison of actual and simulated result

    表  1   鲜莲籽实际落种试验堆积角和休止角测量结果

    Table  1   Measurement results of accumulation angle and repose angle of fresh lotus seeds in actual seed drop test (°)

    试验编号 Test No. 堆积角 Accumulation angle 休止角 Repose angle
    1 2 3 4 1 2
    1 31.47 28.97 29.96 26.50 38.73 35.42
    2 29.70 28.22 26.27 29.62 31.78 37.98
    3 27.91 31.49 23.97 33.84 32.98 36.84
    4 30.26 35.43 29.42 29.36 36.66 33.09
    5 28.52 24.97 32.76 33.83 33.25 36.99
    6 34.32 27.34 27.56 28.50 38.39 37.86
    7 32.89 29.84 29.43 32.89 36.23 34.54
    8 32.92 33.37 32.34 29.62 33.29 37.03
    9 31.90 32.54 26.32 30.14 41.89 35.59
    10 31.85 31.44 27.03 33.19 34.37 36.00
    平均值 Average value 30.20 35.95
    下载: 导出CSV

    表  2   Plackett-Burman试验设计与结果1)

    Table  2   Design and results of Plackett-Burman test scheme

    试验编号 Test No. x1 x2 x3 x4 x5 x6 堆积角/(°) Accumulation angle 休止角/(°) Repose angle
    1 0.5 0.3 0.10 0.5 0.3 0.01 32.24 35.10
    2 0.1 0.9 0.01 0.1 0.9 0.10 35.78 36.03
    3 0.5 0.9 0.10 0.5 0.9 0.10 44.56 50.07
    4 0.5 0.9 0.10 0.1 0.9 0.01 47.02 49.41
    5 0.5 0.3 0.01 0.1 0.9 0.01 25.04 27.62
    6 0.5 0.3 0.10 0.1 0.3 0.10 34.48 35.01
    7 0.5 0.9 0.01 0.5 0.3 0.01 31.39 22.80
    8 0.5 0.9 0.01 0.1 0.3 0.10 32.66 29.52
    9 0.1 0.9 0.01 0.5 0.9 0.01 28.58 32.92
    10 0.5 0.3 0.01 0.5 0.9 0.10 24.02 32.17
    11 0.1 0.3 0.01 0.5 0.3 0.10 25.53 27.68
    12 0.1 0.3 0.10 0.5 0.9 0.01 39.48 40.63
    13 0.1 0.3 0.10 0.1 0.9 0.10 37.71 38.43
    14 0.1 0.9 0.10 0.5 0.3 0.10 44.42 40.10
    15 0.1 0.9 0.10 0.1 0.3 0.01 46.27 44.42
    16 0.1 0.3 0.01 0.1 0.3 0.01 30.09 28.64
     1)x1:莲籽间碰撞恢复系数,x2:莲籽间静摩擦系数,x3:莲籽间滚动摩擦系数,x4:莲籽−有机玻璃碰撞恢复系数,x5:莲籽−有机玻璃静摩擦系数,x6:莲籽−有机玻璃动摩擦系数  1)x1: Coefficient of collision recovery between lotus seeds, x2: Coefficient of static friction between lotus seeds, x3: Coefficient of rolling friction between lotus seeds, x4: Coefficient of collision recovery between lotus seeds and plexiglass, x5: Coefficient of static friction between lotus seeds and plexiglass, x6: Coefficient of dynamic friction between lotus seeds and plexiglass
    下载: 导出CSV

    表  3   Plackett-Burman 试验参数显著性分析1)

    Table  3   Significance analysis of Plackett-Burman test parameters

    方差来源 Source of variance 自由度 Degree of freedom 均方和 Sum of square F P
    堆积角 Accumulation angle 模型 Model 6 137.22 21.49 <0.000 1**
    x1 1 16.91 2.65 0.138 0
    x2 1 240.95 37.74 0.000 2**
    x3 1 541.61 84.84 <0.000 1**
    x4 1 22.16 3.47 0.095 3
    x5 1 1.63 0.26 0.625 3
    x6 1 0.06 0.01 0.927 2
    残差 Residual 9 6.38
    总和 Total 15
    休止角 Repose angle 模型 Model 6 134.13 9.90 0.001 6**
    x1 1 3.20 0.24 0.638 9
    x2 1 99.95 7.37 0.023 8*
    x3 1 573.48 42.31 0.000 1**
    x4 1 3.62 0.27 0.617 8
    x5 1 121.06 8.93 0.015 2*
    x6 1 3.49 0.26 0.624 2
    残差 Residual 9 13.55
    总和 Total 15
     1)x1:莲籽间碰撞恢复系数,x2:莲籽间静摩擦系数,x3:莲籽间滚动摩擦系数,x4:莲籽−有机玻璃碰撞恢复系数,x5:莲籽−有机玻璃静摩擦系数,x6:莲籽−有机玻璃动摩擦系数;“*”和“**”分别表示在P<0.05和P<0.01水平差异显著(F检验)  1) x1: Coefficient of collision recovery between lotus seeds, x2: Coefficient of static friction between lotus seeds, x3: Coefficient of rolling friction between lotus seeds, x4: Coefficient of collision recovery between lotus seeds and plexiglass, x5: Coefficient of static friction between lotus seeds and plexiglass, x6: Coefficient of dynamic friction between lotus seeds and plexiglass; “*” and “**” indicate significant differences at P<0.05 andP<0.01 levels, respectively (F test)
    下载: 导出CSV

    表  4   最陡爬坡试验I设计方案及相对误差1)

    Table  4   Steepest climb test I design scheme and relative error

    试验编号 Test No. x2 x3 x5 堆积角/% Accumulation angle 休止角/% Repose angle
    1 0.3 0.010 0.3 11.51 6.68
    2 0.4 0.025 0.4 9.79 6.27
    3 0.5 0.040 0.5 13.39 18.62
    4 0.6 0.055 0.6 24.29 29.31
    5 0.7 0.070 0.7 31.54 38.57
    6 0.8 0.085 0.8 38.73 41.71
    7 0.9 0.100 0.9 43.89 46.34
     1)x2:莲籽间静摩擦系数,x3:莲籽间滚动摩擦系数,x5:莲籽−有机玻璃静摩擦系数  1) x2: Coefficient of static friction between lotus seeds, x3: Coefficient of rolling friction between lotus seeds, x5: Coefficient of static friction between lotus seeds and plexiglass
    下载: 导出CSV

    表  5   最陡爬坡试验II设计方案及相对误差1)

    Table  5   Design scheme and relative error of steepest climb test II

    试验编号 Test No. x2 x3 x5 堆积角/% Accumulation angle 休止角/% Repose angle
    1 0.30 0.01 0.30 12.43 13.49
    2 0.35 0.02 0.35 7.07 8.27
    3 0.40 0.02 0.40 3.66 1.83
    4 0.45 0.03 0.45 4.48 8.49
    5 0.50 0.03 0.50 10.38 12.53
     1) x2:莲籽间静摩擦系数,x3:莲籽间滚动摩擦系数,x5:莲籽−有机玻璃静摩擦系数  1) x2: Coefficient of static friction between lotus seeds, x3: Coefficient of rolling friction between lotus seeds, x5: Coefficient of static friction between lotus seeds and plexiglass
    下载: 导出CSV

    表  6   料斗排种试验实测值与仿真值误差表

    Table  6   Error table of measurement value and simulation value in hopper seeding test

    试验 Test 试验编号 Test No. 斜率/排种速率(v)/(kg·s−1) Slope/seed discharge rate 截距/莲籽总质量(m0)/kg Intercept/total mass of lotus seeds R2 相对误差/% Relative Error
    v m0
    仿真 Simulated 1.672 10.167 0.999 9
    实测 Actual 1 1.634 10.098 0.999 4 2.30 0.68
    2 1.659 10.174 0.999 1 0.75 0.07
    3 1.670 10.063 0.999 3 0.12 1.02
    4 1.636 10.027 0.999 4 2.18 0.39
    5 1.652 10.169 0.999 1 1.17 0.02
    6 1.643 10.209 0.999 2 1.73 0.41
    7 1.663 10.068 0.998 6 0.52 0.97
    8 1.632 10.199 0.999 3 2.40 0.31
    9 1.611 10.019 0.999 3 3.65 1.46
    10 1.642 9.982 0.999 4 1.81 1.82
    下载: 导出CSV
  • [1] 郑宝东, 郑金贵, 曾绍校. 我国主要莲子品种中三种功效成分的研究[J]. 营养学报, 2004, 26(2): 158-160. doi: 10.3321/j.issn:0512-7955.2004.02.021
    [2] 郑宝东. 莲子科学与工程[M]. 北京: 科学出版社, 2010: 21-47.
    [3]

    ZHANG Y, ZENG H L, WANG Y, et al. Structural characteristics and crystalline properties of lotus seed resistant starch and its prebiotic effects[J]. Food chemistry, 2014, 155: 311-318. doi: 10.1016/j.foodchem.2014.01.036

    [4]

    BHAT R, SRIDHAR K R. Nutritional quality evaluation of electron beam-irradiated lotus (Nelumbo nucifera) seeds[J]. Food Chemistry, 2008, 107(1): 174-184. doi: 10.1016/j.foodchem.2007.08.002

    [5] 侯群喜, 曾荣, 张国忠, 等. 不同成熟度莲籽力学特性试验[J]. 华中农业大学学报, 2020, 39(5): 159-166. doi: 10.13300/j.cnki.hnlkxb.2020.05.022
    [6] 马秋成, 郭耿君, 马婕, 等. 莲仁力学特性参数测定及挤压破碎特性试验[J]. 农业工程学报, 2018, 34(6): 263-271. doi: 10.11975/j.issn.1002-6819.2018.06.033
    [7] 谢丽娟. 莲子破裂最佳位置的有限元模拟及试验研究[D]. 武汉: 华中农业大学, 2005.
    [8] 谢丽娟, 宗力. 莲子受力有限元分析[J]. 农业机械学报, 2006, 37(6): 94-97. doi: 10.3969/j.issn.1000-1298.2006.06.024
    [9] 徐谐庆, 饶洪辉, 李涛, 等. 全自动莲子剥壳去皮机的设计与试验[J]. 农业工程学报, 2014, 30(13): 28-34. doi: 10.3969/j.issn.1002-6819.2014.13.004
    [10] 梁诗华, 林毅鑫, 范泽泽, 等. 鲜莲剥壳技术研究进展与探讨[J]. 湖北农业科学, 2017, 56(4): 603-607. doi: 10.14088/j.cnki.issn0439-8114.2017.04.002
    [11] 廖宜涛, 廖庆喜, 周宇, 等. 饲料油菜薹期收获茎秆破碎离散元仿真参数标定[J]. 农业机械学报, 2020, 51(6): 73-82. doi: 10.6041/j.issn.1000-1298.2020.06.008
    [12] 马文鹏, 尤泳, 王德成, 等. 基于RSM和NSGA-Ⅱ的苜蓿种子离散元模型参数标定[J]. 农业机械学报, 2020, 51(8): 136-144. doi: 10.6041/j.issn.1000-1298.2020.08.015
    [13]

    UCGUL M, FIELKE J M, SAUNDERS C. Three-dimensional discrete element modelling of tillage: Determination of a suitable contact model and parameters for a cohesionless soil[J]. Biosystems Engineering, 2014, 121: 105-117. doi: 10.1016/j.biosystemseng.2014.02.005

    [14] 彭才望, 许道军, 贺喜, 等. 黑水虻处理的猪粪有机肥离散元仿真模型参数标定[J]. 农业工程学报, 2020, 36(17): 212-218. doi: 10.11975/j.issn.1002-6819.2020.17.025
    [15] 郭访全, 钟艳. 太空莲36号高产栽培技术[J]. 农家致富顾问, 2016(8): 4. doi: 10.3969/j.issn.1003-9902.2016.08.004
    [16] 张刘炜, 王海滨, 郭艳玲, 等. 基于响应面法的蓝莓离散元参数标定[J]. 沈阳农业大学学报, 2020, 51(5): 540-548.
    [17] 马秋成, 卢安舸, 高连兴, 等. 莲子物料空气动力学特性与壳仁分离装置试验[J]. 农业工程学报, 2015, 31(6): 297-303.
    [18] 马秋成, 雷林韬, 卢安舸, 等. 基于CFD-DEM耦合法的莲子壳仁分离装置优化与试验分析[J]. 食品与机械, 2016, 32(12): 87-91. doi: 10.13652/j.issn.1003-5788.2016.12.019
    [19] 刘文政, 何进, 李洪文, 等. 基于离散元的微型马铃薯仿真参数标定[J]. 农业机械学报, 2018, 49(5): 125-135. doi: 10.6041/j.issn.1000-1298.2018.05.014
    [20] 侯群喜. 小型鲜食莲籽剥壳通芯一体机设计及剥壳性能试验研究[D]. 武汉: 华中农业大学, 2021.
    [21] 侯占峰, 戴念祖, 陈智, 等. 冰草种子物性参数测定与离散元仿真参数标定[J]. 农业工程学报, 2020, 36(24): 46-54. doi: 10.11975/j.issn.1002-6819.2020.24.006
    [22] 何为, 薛卫东, 唐斌. 优化试验设计方法及数据分析[M]. 北京: 化学工业出版社, 2012: 218-231.
    [23] 吴孟宸, 丛锦玲, 闫琴, 等. 花生种子颗粒离散元仿真参数标定与试验[J]. 农业工程学报, 2020, 36(23): 30-38. doi: 10.11975/j.issn.1002-6819.2020.23.004
    [24] 于庆旭, 刘燕, 陈小兵, 等. 基于离散元的三七种子仿真参数标定与试验[J]. 农业机械学报, 2020, 51(2): 123-132. doi: 10.6041/j.issn.1000-1298.2020.02.014
  • 期刊类型引用(3)

    1. 丁永军,杨文涛,赵一龙. 基于一致性半监督学习的苹果叶片病斑分割模型研究. 农业机械学报. 2024(12): 314-321 . 百度学术
    2. 梁毅. 山东东营市小麦条锈病的发生与防治. 农业工程技术. 2023(22): 34-35 . 百度学术
    3. 强沥文,贾广慧,王伟,周莉. 图像分析技术在农作物污染损害鉴定中的应用. 农业环境科学学报. 2023(12): 2851-2859 . 百度学术

    其他类型引用(7)

图(7)  /  表(6)
计量
  • 文章访问数: 
  • HTML全文浏览量:  0
  • PDF下载量: 
  • 被引次数: 10
出版历程
  • 收稿日期:  2022-04-07
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2023-01-09

目录

/

返回文章
返回