桑园土壤非等径颗粒离散元仿真模型参数标定与试验

为获取土壤离散元仿真模型的土壤颗粒物理参数和接触参数,本文采用试验与仿真相结合的方法,以桑园土壤为例,对土壤颗粒的接触参数进行了仿真标定。首先利用粉体仪、斜面仪、等应变直剪仪等,分析了试验地不同深度土壤的粒径分布,测量了试验地不同深度土壤休止角、滑动摩擦角、剪应力、内聚力、内摩擦角;然后,根据实测土壤粒径分布,利用EDEM软件建立了非等直径土壤球形颗粒模型。在此基础上,以土壤颗粒间及土壤与65Mn钢间的静摩擦因数、滚动摩擦因数、恢复系数为试验因素,土壤休止角、土壤-65Mn钢滑动摩擦角为目标值,建立了基于中心组合试验设计（CCD）方案,并利用Design-Expert软件对仿真试验结果进行了分析,得到了仿真标定的土壤-土壤间静摩擦因数、滚动摩擦因数和恢复系数的最优值分别为0.89、0.45和0.43;标定的土壤-65Mn钢间静摩擦因数、滚动摩擦因数和恢复系数的最优值分别为1.15、0.05和0.4。利用以上标定的最优参数对桑园土壤进行了休止角与滑动摩擦角仿真试验,试验结果表明,休止角仿真值与试验值相对误差为1.69%,土壤-65Mn钢的滑动摩擦角仿真值与试验值相对误差为2.88%。在此...     

桑园土壤非等径颗粒离散元仿真模型参数标定与试验

为获取土壤离散元仿真模型的土壤颗粒物理参数和接触参数,本文采用试验与仿真相结合的方法,以桑园土壤为例,对土壤颗粒的接触参数进行了仿真标定。首先利用粉体仪、斜面仪、等应变直剪仪等,分析了试验地不同深度土壤的粒径分布,测量了试验地不同深度土壤休止角、滑动摩擦角、剪应力、内聚力、内摩擦角;然后,根据实测土壤粒径分布,利用EDEM软件建立了非等直径土壤球形颗粒模型。在此基础上,以土壤颗粒间及土壤与65Mn钢间的静摩擦因数、滚动摩擦因数、恢复系数为试验因素,土壤休止角、土壤-65Mn钢滑动摩擦角为目标值,建立了基于中心组合试验设计（CCD）方案,并利用Design-Expert软件对仿真试验结果进行了分析,得到了仿真标定的土壤-土壤间静摩擦因数、滚动摩擦因数和恢复系数的最优值分别为0.89、0.45和0.43;标定的土壤-65Mn钢间静摩擦因数、滚动摩擦因数和恢复系数的最优值分别为1.15、0.05和0.4。利用以上标定的最优参数对桑园土壤进行了休止角与滑动摩擦角仿真试验,试验结果表明,休止角仿真值与试验值相对误差为1.69%,土壤-65Mn钢的滑动摩擦角仿真值与试验值相对误差为2.88%。在此基础上,依据实测的土壤剪应力,采用试错法,以实测土壤内摩擦角为目标值,优化标定了土壤-土壤颗粒Hertz-Mindlin with Bonding接触模型中的粘结参数,标定法向粘结刚度、切向粘结刚度分别为1×108、5×107N/m3,临界法向应力和临界切向应力均为10kPa,接触半径为1.1倍颗粒半径,直剪仿真得到内摩擦角为30.24°,仿真值与直剪试验内摩擦角平均值相对误差为5.53%。本文提出的土壤颗粒建模方法、标定方法及其所标定的参数值,可用于砂质壤土桑园耕作机械触土部件与土壤相互作用的离散元仿真分析及其结构优化。     

基于离散元的赤芍移栽接触参数标定

为确定赤芍在移栽过程中的仿真参数以进行离散元仿真试验,以物理试验测量土壤堆积角为基础,选取土壤间的静摩擦因数、滚动摩擦因数、碰撞恢复系数以及土壤的表面能作为试验因素,以土壤堆积角为试验指标,进行四因素三水平离散元仿真试验设计,建立了土壤接触参数与堆积角的回归模型。以实际堆积角值为目标进行了回归模型参数优化。通过实际物理试验完成了赤芍与土壤间的静摩擦系数、碰撞恢复系数、滚动摩擦系数以及碰撞恢复系数的标定。试验结果表明,土壤接触参数近似组合为碰撞恢复系数为0.223、静摩擦因数为0.630、滚动摩擦因数为0.373、土壤表面能为0.340 J/m³,该组合离散元仿真所得土壤堆积角为31.57°,与实际测量结果误差仅为1.2%。赤芍与土壤接触参数为：碰撞恢复系数为0.25、静摩擦因数为0.84、滚动摩擦因数为0.37。研究结果可以为后期移栽离散元仿真分析奠定基础。     

蚕豆物性参数测定与离散元仿真参数标定

采用离散元法对蚕豆与收获机具间的相互作用进行数值仿真研究,实现对蚕豆收获机具关键部件相关参数的设计与优化,以提高研究效率、减少研究成本。为此,通过物理试验对蚕豆基本物性参数、接触力学参数和休止角进行测定,参考物理试验测定结果选取仿真试验参数取值范围。以仿真试验休止角为响应值,由Plackette-Burman(PB)试验得到显著影响休止角的试验参数(蚕豆-蚕豆静摩擦因数、蚕豆-钢板静摩擦因数、蚕豆-钢板滚动摩擦因数),并利用最陡爬坡试验在PB试验基础上得到显著参数的取值范围。采用Box-Behnken(BB)试验在休止角与显著参数之间建立二阶回归模型,以休止角物理试验测定值为目标值,对显著参数进行寻优,获得最优参数组合,即蚕豆-蚕豆静摩擦因数0.25、蚕豆-钢板静摩擦因数0.47、蚕豆-钢板滚动摩擦因数0.35。最后,通过T检验得到P>0.05,表明采用最优参数组合得到的休止角仿真值与实测值相差无几,验证了最优参数组合的可靠性。     

果荚初期饲料油菜茎秆离散元接触模型参数标定

针对饲料油菜与不同材料的接触参数实测难度大、机械化收获离散元仿真模拟缺乏接触模型参数的问题，以果荚初期饲料油菜为对象，基于EDEM开展了饲料油菜茎秆颗粒离散元接触模型参数标定。测定了果荚初期饲用油菜茎秆本征参数，茎秆平均直径为20.4mm，密度为809kg/m3，茎秆弹性模量、剪切模量和泊松比平均值分别为115.73MPa、47.04MPa和0.23；以休止角为评价指标，应用Hertz-Mindlin基本模型和圆筒提升堆积法开展了饲料油菜茎秆颗粒堆积的虚拟二水平因子试验，结果表明饲料油菜茎秆与钢之间的碰撞恢复系数和滚动摩擦因数以及茎秆之间的碰撞恢复系数对休止角的影响较小，其值分别为0.60、0.10和0.60；通过最陡爬坡试验和响应面分析，确定了饲料油菜茎秆颗粒间静摩擦因数、滚动摩擦因数和饲料油菜茎秆-钢静摩擦因数的取值范围，建立了颗粒休止角的回归模型，以实测休止角与仿真试验休止角之间相对误差最小进行响应面分析和优化求解，确定其参数值分别为0.36、0.03和0.23。在接触参数最优组合条件下，根据回归模型计算得出的休止角理论值与实测值误差为2.15%，仿真试验得出休止角模拟值与实测值误差为1.83%，表明标定方法正确，标定参数准确。研究可为饲料油菜机械化收获过程的离散元仿真分析提供基本参数。     

基于EDEMpy的枸杞离散元仿真参数标定与试验

为探究枸杞离散元仿真边界参数的最优组合,提高离散单元仿真方法研究枸杞与采收机械作用机理的准确性,提升后续枸杞机械采收设备的研发效率及机械作业性能,选取EDEM中“Hertz-Mindlin with Johnson-Kendall-Roberts”凝聚力接触模型,通过物理试验测量枸杞堆积角、枸杞-钢板静摩擦因数、枸杞-钢板碰撞恢复系数。在物理试验的基础上,采用基于EDEMpy的仿真试验方法对枸杞的接触参数进行试验研究,以枸杞颗粒的静摩擦因数、滚动摩擦因数、碰撞恢复系数、表面能JKR为影响因素,以堆积角为目标进行多因素试验,并建立回归模型。试验结果表明：对枸杞仿真堆积角影响显著的参数为枸杞-枸杞静摩擦因数、枸杞-枸杞滚动摩擦因数、表面能JKR。以实际堆积角为目标值,对回归模型进行优化分析,得到显著性参数最优组合值分别为：枸杞-枸杞静摩擦因数0.506,枸杞-枸杞滚动摩擦因数0.064,表面能JKR0.048,该组合下的仿真堆积角平均值为28.82°,与实际堆积角的偏差小于4%,边界参数标定结果可靠。研究结果可为枸杞离散元模型边界参数的选取提供参考。     

黑土区玉米秸秆-土壤混料离散元模型参数标定

由于黑土区保护性耕作中关键农机部件设计优化过程中缺乏准确的离散元仿真模型参数,在一定程度上制约了机具的优化改进。以黑土区玉米秸秆-土壤混料为研究对象,构建玉米秸秆-土壤混料离散元仿真模型,采用物理试验与EDEM仿真试验相结合的方法,选用Hertz-Mindlin with JKR接触模型进行离散元仿真接触参数标定。首先,采用圆筒提升的方法确定玉米秸秆-土壤混料的实际堆积角,利用Design-Expert软件中Plackett-Burman试验筛选出对堆积角有显著影响的参数为：土壤-土壤滚动摩擦因数、土壤-钢静摩擦因数、秸秆-土壤滚动摩擦因数、土壤JKR表面能;进一步通过最陡爬坡试验确定4个参数的最优取值范围,根据Box-Behnken试验原理以堆积角为响应值,建立堆积角与显著参数的二次回归模型;以实际堆积角为目标,利用软件寻优功能对显著参数进行优化并得到最优参数组合：秸秆-土壤滚动摩擦因数0.16、土壤-土壤滚动摩擦因数0.24、土壤-钢静摩擦因数0.75、土壤JKR表面能0.67J/m2。通过仿真试验对最优参数组合进行对比验证,仿真堆积角与物理试验堆积角相对误差为1.64%。研究结果表明标定的参数真实可靠,可为黑土区玉米秸秆-土壤混料的离散元仿真提供理论参考。     

基于注入截面法的颗粒饲料离散元模型参数标定

采用离散元法进行颗粒饲料后喷涂、冷却、输送、仓储、饲喂等关键环节的工作过程仿真分析时,需要建立颗粒饲料离散元模型并通过测定休止角来标定离散元参数,用以测定休止角的传统装置与方法,存在样品用量多、测定繁琐的情况。本文提出了一种基于注入截面法的休止角测定装置与方法,通过颗粒堆积体截面的轮廓线直接获取休止角,从而进行休止角的模拟与测定。建立与颗粒饲料形态相近的离散元模型,借用GEMM(Generic EDEM material model database)数据库获得离散元模型参数范围。以滑动摩擦因数X1、碰撞恢复系数X2和滚动摩擦因数X3为试验因素,以颗粒饲料堆积休止角Y1为评价指标,按照3因素5水平正交旋转组合设计试验方法,利用DesignExpert 8.0.6软件回归分析法和响应面分析法,建立了3个因素以堆积休止角为评价指标的数学模型。以颗粒饲料休止角真实测得值为目标,对回归模型进行寻优,得到优化后的标定参数组合:颗粒饲料间滑动摩擦因数为0.41,碰撞恢复系数为0.53,滚动摩擦因数为0.08。以此优化解进行仿真试验,结果显示预测休止角为29.43°±0.70°,误差为3.1%,休止角仿真和试验在堆积角度和形态上相似度较高。结果表明了基于注入截面法的颗粒饲料离散元建模与休止角测定试验的有效性和可行性。     

蒜种颗粒离散元模型参数标定

为确定蒜种离散元模型仿真参数,采用多球聚合的方法建立了蒜种离散元模型,并对可实现测量的模型参数(蒜种-钢板恢复系数、蒜种-蒜种恢复系数、蒜种-钢板静摩擦因数、蒜种-蒜种静摩擦因数)进行了物理试验测定。对不易测量的模型参数(蒜种-钢板滚动摩擦系数、蒜种-蒜种滚动摩擦系数)进行了仿真试验标定,应用、建立了两种不易测量的接触参数与蒜种堆休止角、密相区域圆半径的二次回归模型,最终求解寻优得到参数最优值。通过上述试验方法,得到蒜种离散元模型参数:蒜种-钢板恢复系数为0.511,蒜种-蒜种恢复系数为0.487,蒜种-钢板静摩擦因数为0.473,蒜种-蒜种静摩擦因数为0.503,蒜种-钢板滚动摩擦因数为0.111,蒜种-蒜种滚动摩擦因数为0.108。最后,在得到的接触参数下进行仿真验证试验,结果表明:蒜种堆休止角、密相区域圆半径的仿真值与实际值的相对误差分别为1.36%、1.50%,无明显差异,排种功率、排种速度的仿真与试验值在整体变化趋势上具有一致性,验证了蒜种离散元模型与仿真试验的有效性,可为大蒜播种机械仿真设计与优化提供参考。     

酒用高粱种子离散元仿真参数标定

为获得酒用高粱种子离散元仿真参数,以仁怀市大坝镇地区的茅梁糯2号高粱种子为研究对象,完成了物性参数的测量种子-尼龙、种子-ABS及种子颗粒间接触参数的标定。首先,采用物理试验测出高粱种子的三轴尺寸、密度、泊松比、弹性模量、剪切模量及休止角等参数,且获取高粱种子-尼龙与高粱种子-ABS的静摩擦因数、滚动摩擦因数和碰撞恢复系数;其次,采用Box-Behnken法设计三因素三水平试验,得到高粱种子颗粒间接触参数与休止角的二阶回归模型;最后,以实测休止角29.14°为优化目标值,求得高粱种子颗粒间的碰撞恢复系数、静摩擦因数、滚动摩擦因数分别为0.422、0.504、0.074,代入EDEM中得到仿真试验结果为28.39°,与实测值的相对误差为2.57%,进一步验证了仿真试验的可靠性和有效性。酒用高粱种子物性参数测定和接触参数标定可为高粱播种机仿真分析及优化设计提供一定参考和依据。     

黄土高原坡地土壤与旋耕部件互作离散元仿真参数标定

针对黄土高原坡地土壤-旋耕部件互作机理研究以及坡地专用旋耕机具设计缺乏准确可靠离散元仿真参数的问题,以典型坡地粘壤土(含水率13.4％±1％)为研究对象,选取EDEM中Hertz-Mindlin with JKR Cohesion接触模型,对相关仿真参数进行标定.首先,对土壤颗粒间接触参数进行了标定,以土壤颗粒的仿真堆...     

茎瘤芥种子离散元模型参数标定

为解决茎瘤芥种子离散元模拟中缺乏准确仿真参数的问题,通过直接测量和虚拟标定方法确定“涪杂2号”茎瘤芥种子离散元参数。利用数显游标卡尺、电子秤等物理机械特性测定工具,得到茎瘤芥种子的三轴尺寸、等效半径、千粒重和密度;通过斜面碰撞法、三点滑动法和斜面法分别得到茎瘤芥种子与有机玻璃之间的恢复系数、静摩擦系数和滚动摩擦系数分别为0.561、0.420、0.155;采用旋转鼓试验和圆筒提升试验结合MATLAB图像处理技术获得茎瘤芥种子的动态堆积角和静态堆积角分别为35.235 2°和14.483 0°。基于Plackett Burman试验筛选出对动态堆积角和静态堆积角影响显著的参数因子均为种子与种子间的滚动摩擦系数、种子与种子间的静摩擦系数;通过中心复合设计分别获得显著性因子与动态堆积角、静态堆积角的二次回归方程,并以动态堆积角和静态堆积角的实测值为目标,基于回归模型对因子参数组合进行优化,得到两组种子与种子间的静摩擦系数、种子与种子间的滚动摩擦系数最优参数组合为0.666 1、0.023 7以及0.341 6、0.036 0。对两组最优参数组合进行仿真验证,结果显示两组参数组合下动态堆积角、静态堆积角的仿真值与实测值相对误差分别为0.477%、3.590%和1.820%、4.950%。结果表明,仿真中种子与种子间的静摩擦系数、种子与种子间的滚动摩擦系数分别为0.666 1、0.023 7时,更符合茎瘤芥种子实际的物理特性。本研究能为茎瘤芥播种过程及相关装置设计优化提供参考。     

颗粒肥料离散元仿真边界参数系统化研究

离散元边界参数对仿真精确性具有重要的影响。为提高离散元仿真精度,以大颗粒尿素颗粒为研究对象,利用Plackett-Burman休止角仿真试验进行了重要边界参数的筛选,确定影响显著的参数依次为尿素颗粒间滚动摩擦因数、颗粒间静摩擦因数和颗粒与ABS板间静摩擦因数,且3种边界参数影响尿素颗粒堆积特性,休止角随着3种边界参数的增大而增大。利用自制静摩擦因数测量仪和虚拟仿真标定方法分别对颗粒与ABS板间的静摩擦因数、颗粒间静摩擦因数和颗粒间滚动摩擦因数进行研究,并对确定值进行了堆积过程的仿真和试验验证,仿真休止角与实际试验休止角相对误差仅为0.36%,不同含水率下的实际试验休止角与标定参数下的仿真休止角相对误差均不大于3.25%,表明仿真确定的边界参数和仿真模型的有效性。     

基于EDEM的青稞接触参数仿真标定

针对青稞品种多样,且部分接触参数难以直接测量的问题,本研究以青稞的休止角为响应值,对青稞的接触参数进行标定。采用响应面法:通过Plackett-Burman试验对青稞的8个离散元仿真参数进行筛选,得到对青稞休止角有显著性影响的三个因素:青稞—青稞的静摩擦系数、青稞—青稞的滚动摩擦系数、青稞—pla平板的静摩擦因数;通过最陡爬坡试验确定三个显著因素的最优取值范围;在此基础上进行Box-Behnken试验,建立休止角与显著性参数的二阶回归模型,并对方程进行优化,求得显著性参数的最优参数值:青稞—青稞的静摩擦系数0.19、青稞—青稞的滚动摩擦系数0.01、青稞与pla平板的静摩擦因数0.43。最后进行仿真验证试验,结果表明仿真试验所得休止角与实际休止角值误差为0.64%,误差很小表明仿真标定的参数是可信的,可以为今后的青稞排种器的研究提供参考。     

离散元模拟中沙土参数标定方法研究

农业机械与土壤相互作用仿真时,选用颗粒相互作用参数的准确度将直接影响仿真结果。本文提出一种通过试验与模拟相结合系统地标定沙土颗粒相互作用参数的方法。通过堆积角测试装置、三轴剪切试验仪、图像颗粒分析系统等设备测量计算沙土的堆积角、剪切模量、粒径分布和外观形貌等参数,为颗粒或工厂建模提供参考。使用标准球和非标准球对沙土颗粒的碰撞恢复因数、静摩擦因数、滚动摩擦因数进行标定。研究不同质量和不同标定方法(抽板法和漏斗法)是否对堆积角产生影响。模拟结果表明,选用标准球标定时,碰撞恢复因数是0.15,静摩擦因数是0.8,滚动摩擦因数是0.2,抽板法得到的堆积角是33.99°,相对误差是4.74%;漏斗法得到的堆积角是33.59°,相对误差是3.51%。同时,选用非标准球进行标定时,碰撞恢复因数是0.15,静摩擦因数是0.2,滚动摩擦因数是0.3,抽板法得到的堆积角是32.06°,相对误差是1.20%。由此看出,颗粒外观形貌对颗粒间静摩擦因数影响相对较大。     

颗粒饲料破损离散元仿真参数标定

为了探究气力输送中颗粒饲料的破损机理,针对当前缺乏颗粒饲料准确破损仿真模型的问题,利用EDEM仿真软件进行颗粒饲料破损离散元仿真参数标定研究。以粒径为2.50 mm混养成鱼颗粒饲料为研究对象,通过基础试验测定了颗粒饲料本征参数;通过颗粒饲料休止角试验、碰撞恢复系数标定试验和落料时间,结合试验优化设计方法,确定了饲料间的碰撞恢复系数、静摩擦因数、滚动摩擦因数为0.58、0.23、0.12,饲料和软塑料（软PVC）间的碰撞恢复系数、静摩擦因数、滚动摩擦因数为0.69、0.22、0.18;通过颗粒饲料单轴压缩破碎试验和仿真试验,结合响应面优化确定了单位面积法向刚度、单位面积切向刚度、临界法向应力、临界切向应力,分别为2.25×10⁹ N/m³、8.05×10⁸ N/m³、455 MPa、305 MPa。以确定的参数进行休止角仿真试验、单轴压缩仿真试验,结果表明,休止角、破碎力、落料时间的仿真值与实测值相对误差分别为0.35%、1.43%、2.81%;通过自由落料、斜面滑动、斜面滚动试验对粘结模型接触参数进...     

颗粒肥料离散元仿真摩擦因数标定方法研究

对分体圆筒法、倾斜法、抽板法和斜面法4种颗粒特性测试方法进行Plackett-Burman多因素显著性筛选试验,试验方差分析结果表明,不同的测试方法影响测量结果的显著因素与因素显著程度。根据分体圆筒法、倾斜法和斜面法的方差分析结果,提出一种基于颗粒物料整体特性的摩擦因数标定方法,将仿真试验与真实试验相结合,依次标定出尿素颗粒与PVC材料间静摩擦因数为0. 41,颗粒间静摩擦因数为0. 36,颗粒间滚动摩擦因数为0. 15。将所标定的摩擦因数采用无底圆筒法进行验证试验,休止角仿真试验结果为30. 57°,真实试验结果为31. 74°,相对误差为3. 69%,不同含水率下的实际试验休止角与所标定摩擦因数下的仿真休止角相对误差均不大于4. 59%,仿真试验结果与真实试验结果无显著差异,验证了所标定摩擦因数的有效性。本方法可用于其他颗粒状物料间摩擦因数的标定试验。