蔗段离散元仿真建模方法与参数标定

【目的】蔗段作为大长径比秆状物料，其离散元模型的构建方法与仿真参数的设定尚不清楚，仿真模型精度对颗粒间的动力学响应特性有较大影响，需通过参数标定提高仿真参数的准确性。【方法】以蔗段物理堆积角为响应值，采用仿真试验方法优化标定离散元参数。首先，采用物理试验测定蔗段的基本物性参数，并基于多球聚合模型和XML的方法构建蔗段仿真模型；然后，应用Plackett-Burman试验对蔗段离散元仿真中的8个初始参数进行显著性筛选，并对显著性参数进行最陡爬坡试验，确定最优参数区间；最后，基于Box-Behnken试验建立显著性参数与堆积角的二阶回归方程，以物理试验堆积角42.70°为目标值，对回归方程进行优化求解。【结果】显著性筛选试验得出蔗段泊松比、蔗段-蔗段静摩擦系数、蔗段-蔗段滚动摩擦系数对仿真堆积角影响显著；最优参数组合为：蔗段泊松比0.35、蔗段-蔗段静摩擦系数0.53、蔗段-蔗段滚动摩擦系数0.04。最优参数组合的仿真试验结果表明，仿真堆积角与物理试验堆积角无显著性差异，两者相对误差为0.99%，进一步验证了蔗段离散元标定参数的可靠性。【结论】蔗段离散元模型与最优仿真参数可用于蔗段离散元仿...     

基于响应面法的玉米籽粒离散元参数标定

【目的】离散元法普遍用于模拟农业物料在农业生产过程中的运动状态。离散元模拟的准确度取决于所使用的离散元参数,对颗粒离散元参数进行标定,为其在仿真试验中的应用提供可靠基础。【方法】以玉米籽粒为研究对象,以堆积角为响应值,将采用不同参数组合的仿真试验结果与实体试验结果对比,提出一种基于响应面法的离散元参数标定方法。采用Plackett-Burman试验对玉米籽粒离散元参数进行显著性检验,筛选对响应值影响显著的因素;根据最陡爬坡试验确定显著性因素的最佳水平范围;采用中心组合设计进行3因素5水平响应面优化试验。【结果】玉米泊松比与玉米–玉米静摩擦系数交互作用显著,得到玉米籽粒离散元参数最佳组合为玉米泊松比0.438、玉米–玉米静摩擦系数0.182、玉米–玉米滚动摩擦系数0.051。采用标定好的参数进行仿真试验得到堆积角平均值为26.89°,接近实体试验的堆积角值。【结论】基于响应面法的玉米籽粒离散元参数标定方法是可行的,可以提高玉米籽粒离散元仿真试验的准确性。     

前胡种子物性参数测定及其离散元仿真模型参数标定

【目的】为确定前胡排种器离散元仿真所需的物性参数,对前胡种子物性参数进行测量及标定.【方法】采用空心圆筒堆积试验法与EDEM仿真结合的方法,结合图像处理技术测定前胡种子堆积角,基于响应面优化法对前胡种子物性参数进行标定,采用Plackett-Burman法从7组因素中筛选出了2组前胡种子堆积角影响显著的因素,再根据星点设计-效应面法试验原理,对前胡种子与前胡种子静摩擦系数、前胡种子与前胡种子滚动摩擦系数进行响应面分析试验.用效应面法预测影响显著因素的最佳参数组合.【结果】前胡种子密度为0.129 g/cm~3,前胡种子与尼龙板静摩擦系数为0.55,前胡种子与前胡种子静摩擦系数0.53,前胡种子与前胡种子滚动摩擦系数为0.05,将最佳参数组合进行仿真试验得到的堆积角与真实试验所得堆积角进行对比,相对误差在2.05%以内,在允许误差范围内.【结论】通过EDEM快速填充的方法建立了前胡种子离散元颗粒模型,为后续前胡的研究奠定了基础.通过空心圆筒法结合图像处理技术,获得了前胡种子堆积角均值为40.93°.     

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

【目的】确定鲜莲籽机械化加工过程中离散元仿真模型参数,为鲜莲籽机械化加工仿真试验提供数据参考。【方法】本研究利用EDEM仿真软件开展鲜莲籽离散元仿真参数标定。以产自湖北洪湖的‘太空莲36号’为试验对象,通过落种试验测定鲜莲籽实际落种的堆积角和休止角。基于Hertz-Mindlin (no slip)接触模型进行鲜莲籽落种仿真试验,以鲜莲籽堆积角和休止角的实测值与仿真值之间的误差为试验指标,通过PlackettBurman试验确定对堆积角和休止角影响显著的接触参数,通过最陡爬坡试验确定鲜莲籽离散元模型最优接触参数组合。采用料斗进行实际落种验证试验,以莲籽落种速率为试验指标,对比实际与仿真落种验证试验莲籽落种速率,验证最优参数组合可靠性。【结果】莲籽间静摩擦系数、莲籽间滚动摩擦系数对堆积角影响极显著(P<0.01);莲籽间滚动摩擦系数对休止角影响极显著(P<0.01),莲籽间静摩擦系数、莲籽-有机玻璃静摩擦系数对休止角影响显著(P<0.05)。最优接触参数组合为莲籽间静摩擦系数0.4、莲籽间滚动摩擦系数0.02、莲籽-有机玻璃静摩擦系数0.4。落料验证试验结果表明,实际试...     

新疆农田粉土离散元仿真参数标定

【目的】标定新疆地区粉土离散元仿真模型,量化分析农业机械触土部件与土壤的相互作用。【方法】采用EDEM离散元仿真软件,运用斜坡试验标定土壤与触土材料（65 Mn）之间的接触参数,堆积角试验标定土壤与土壤之间的接触参数,以斜坡滚动距离为优化目标进行回归分析。【结果】土壤与65 Mn之间离散元模型参数的一种优化组合为恢复系数0.51,静摩擦系数0.56,动摩擦系数0.08,JKR表面能4.12;以堆积角为优化目标,土壤之间离散元模型参数的一种优化组合为恢复系数0.57,静摩擦系数0.65,动摩擦系数0.23,JKR表面能4.49。滚动距离和堆积角的误差分别为6.05%、1.28%。【结论】在建立的仿真模型条件下,在8、9、10 km/h作业速度下,犁体仿真试验与实际试验的工作阻力相对误差为5.69%、5.95%、6.49%,建立的模型真实可靠。     

基于离散元法的水稻种子含水率表征方法

【目的】针对离散元仿真软件EDEM中的输入参数缺少含水率,无法对水稻种子的物理特性进行准确描述的问题,探究表征种子含水率的参数并提出标定方法.【方法】利用圆筒试验法和EDEM仿真试验,找到了1种用颗粒间静摩擦系数来表征水稻含水率的方法,并进行了验证.首先进行了5种不同含水率的水稻种子物理特性试验,验证含水率对水稻种子物理特性的影响,为仿真所需的输入参数提供依据.随后利用圆筒法测量水稻种子堆积角,得到堆积角与含水率的关系,并在EDEM软件中重复该过程,通过单因素重复试验找到对堆积角影响最显著的参数为颗粒间静摩擦系数,得到堆积角与颗粒间静摩擦系数的关系.【结果】颗粒间静摩擦系数为水稻种子含水率的表征参数,其误差小于5%.【结论】在EDEM软件中用颗粒间静摩擦系数表征水稻种子含水率的方法可行,并为其他种子在EDEM软件中含水率的表征提供了思路.     

水稻机械直播的土壤表层泥浆离散元参数标定

【目的】针对水稻机械直播稻田播前土壤与触土部件的相互作用机理不明确等问题,可利用离散元分析法研究影响水稻精量穴直播质量的因素。由于播前土壤的结构复杂,其表层泥浆与耕层存在一定差异,因此,需分层构建复合土壤模型。为获取精确的水稻机械直播的表层泥浆离散元参数,开展试验研究和仿真分析,并标定表层泥浆的离散元仿真参数。【方法】按照水稻精量穴直播的整地要求处理,并在播种前去田间获取表层泥浆,采用漏斗法进行泥浆堆积角试验,利用数显倾角仪多次测量泥浆堆积角并取平均值。选择EDEM离散元软件中Hertz-Mindlin with JKR接触模型,开展泥浆堆积角仿真试验。以堆积角为响应值,通过Plackett-Burman试验筛选出对堆积角影响显著的3个参数,进一步开展最陡爬坡试验缩小显著性参数的取值范围。采用Box-Behnken试验建立表层泥浆堆积角与筛选的显著性参数的回归模型,并以物理试验测得的泥浆堆积角为目标值,对显著性参数寻优得到最佳参数组合。将最优参数代入仿真软件验证表层泥浆颗粒离散元参数的准确性。【结果】泥浆堆积角物理试验获取表层泥浆堆积角为40.20°。Plackett-Burman试验...     

考虑颗粒间黏结力的黏性土壤离散元模型参数标定

【目的】实现黏性土壤离散元模型的接触参数与接触模型参数标定。【方法】基于土壤堆积角物理试验结果,采用考虑颗粒间黏结力的"Hertz-Mindlin with JKR"接触模型进行土壤堆积角仿真试验,借助GEMM(Generic EDEM material model database)数据库获得离散元模型关键参数(包括JKR表面能、恢复系数、静摩擦系数与动摩擦系数),进一步运用Box-Behnken试验方法进行堆积角仿真试验。【结果】通过对试验结果进行多元回归拟合分析获得了堆积角回归模型,回归模型的方差分析表明该模型极显著,试验因素对堆积角的影响为二次多项式,且存在复杂的一次与二次交互作用。以堆积角40.45°为目标对回归模型进行寻优,得到了优化解:JKR表面能7.91J·m-2;恢复系数0.66;静摩擦系数0.83;动摩擦系数0.25。以此优化解进行仿真试验获得的堆积角为39.73°。堆积角仿真试验与物理试验在堆积角度和形状上具有较高的相似性。【结论】可利用该优化参数对样品土壤进行进一步的黏性土壤与触土部件间的离散元仿真,从而揭示黏性土壤在触土部件作用下的运动规律。     

基于离散元的竹粉颗粒接触参数标定

【目的】建立竹粉颗粒材料的离散元模型并标定其相关接触参数。【方法】基于毛竹Phyllostachys edulis粉末物理堆积角实验结果,采用离散元法(DEM)模拟仿真和实验设计(DOE)相结合的方法,通过Plackett-Burman (P-BD)实验、爬坡实验和响应面实验获得竹粉堆积角和仿真参数之间的二次多项式回归模型,以物理堆积角为目标值预测得到接触参数的最优组合。【结果】实测竹粉物理堆积角为49.29°;P-BD实验结果表明了竹粉-竹粉滚动摩擦系数、竹粉-竹粉恢复系数和竹粉-不锈钢板静摩擦系数对堆积角影响显著(P<0.05);响应面实验结果进一步证实了P-BD实验的结果,同时表明竹粉-竹粉滚动摩擦系数自身交互作用、竹粉-竹粉滚动摩擦系数和竹粉-竹粉恢复系数的交互作用对堆积角有显著影响(P<0.05);最优因素组合为竹粉-竹粉恢复系数0.44,竹粉-竹粉滚动摩擦系数0.22,竹粉-不锈钢板静摩擦系数0.45。【结论】本研究在41.20°～55.27°堆积角范围内建立的竹粉颗粒堆积角定量模拟标定与预测平台具有较高的可靠性,DEM标定方法可为后续建立评估生物质颗粒运动行为...     

基于堆积试验的玉米包衣种子离散元参数标定

为获得基于黏结颗粒模型（BPM）玉米包衣种子离散元仿真所需的精确接触参数，基于堆积试验对玉米包衣种子的仿真参数进行标定。对郑单958玉米包衣种子分类筛选后，通过激光扫描仪对轮廓较好的种子进行扫描，得到点云数据，并通过CATIA软件对点云数据进行处理，最终得到玉米种子仿真模型。设计Plackett-Burman试验，通过Isight软件的近似模型与DOE联合模块对试验结果进行分析，筛选出对堆积角影响显著的参数：玉米种子-玉米种子静摩擦系数、法向刚度与切向刚度。基于Isight软件RSM优化模块，根据Box-Behnken试验结果建立堆积角与显著性参数的二阶回归模型，得到参数的最佳组合：玉米种子-玉米种子静摩擦系数0.269、法向刚度2.54×10⁸ N·m^-3、切向刚度5.93×10⁷ N·m^-3。将标定参数仿真所得堆积角与真实试验值进行对比，二者相对误差为0.98%。上述结果表明，响应面分析可用于标定玉米包衣种子的离散元仿真参数，为玉米气力式排种器的结构设计与参数优化提供参考。     

白萝卜种子颗粒模型离散元接触参数标定与试验

在利用EDEM与FLUENT耦合模拟白萝卜排种器排种过程时,白萝卜种子颗粒参数设置以及模型的选择直接影响仿真结果的可靠性。结合白萝卜种子真实颗粒和仿真堆积试验,标定了不同填充球颗粒半径白萝卜种子模型的2个主要接触参数：白萝卜种子-有机玻璃静摩擦系数和白萝卜种子间静摩擦系数。采用自动填充方式创建了不同填充球颗粒半径白萝卜种子的离散元模型。利用Plackett-Burman试验对物料特性影响参数进行分析,发现白萝卜种子-有机玻璃静摩擦系数、白萝卜种间静摩擦系数对堆积角的影响极显著。结合台架与仿真堆积试验,建立了2个主要接触参数与堆积角的二元回归模型,以白萝卜种子实际堆积角为目标对参数进行寻优,得到不同填充球颗粒半径白萝卜种子模型的白萝卜种子-有机玻璃静摩擦系数和白萝卜种子间静摩擦系数。结合堆积角相对误差率与仿真时间分析最佳球颗粒填充半径,当白萝卜种子离散元填充颗粒半径为0.25 mm时,其仿真精度和仿真时间最优。得到的最佳填充球颗粒模型以及标定的接触参数,为白萝卜种子排种器研究提供一定的参考。     

蓝亚麻蒴果离散元模型仿真参数标定

【目的】对蓝亚麻蒴果离散元模型仿真参数进行标定.【方法】通过查阅文献法、物理试验法测定了蓝亚麻蒴果的基本物理参数(千粒质量、含水率、三轴尺寸、密度、泊松比、剪切模量)和接触力学参数(碰撞恢复系数、静摩擦系数和滚动摩擦系数)的取值或取值范围.以物理试验法测定的各参数为仿真试验参数选择依据,利用Minitable18软件进行Plackett-Burman试验设计,对物理试验法粗测得的参数进行显著性筛选;进一步以实际休止角与仿真试验休止角的相对误差作为评价指标,对显著性参数进行最陡爬坡试验,通过仿真逼近预测法对显著性参数进行寻优.最后,以较优参数组合进行仿真试验验证,测得对应仿真试验休止角,并与实际休止角进行对比.【结果】Plackett-Burman试验结果表明:蓝亚麻蒴果间的滚动摩擦系数对仿真试验蒴果休止角影响显著;由最陡爬坡试验结果可得蓝亚麻蒴果间的滚动摩擦系数为0.18;仿真试验验证结果表明:仿真试验休止角与实际休止角相对误差为0.95%.【结论】应用上述各方法及优化试验来标定蓝亚麻蒴果离散元模型仿真参数是可行的.     

基于堆积试验的沙蓬种子离散元仿真参数标定

为了获得沙蓬种子仿真播种过程所需的参数,通过物理堆积试验与仿真堆积试验相结合的方法对沙蓬种子的离散元仿真参数进行标定。首先,对沙蓬种子相关物性参数进行了测定,建立其离散元模型,并利用物理堆积试验获得沙蓬种子的堆积角为45.995°。其次,采用Plackett-Burman试验进行仿真堆积试验,筛选出对堆积角影响显著的参数为：沙蓬种子之间恢复系数、沙蓬种子之间滚动摩擦系数、沙蓬种子与ABS板间静摩擦系数。然后,根据最陡爬坡试验确定显著参数的取值范围;设计二次正交旋转组合试验,建立堆积角相对误差与显著参数的二阶回归方程,得到离散元仿真最佳参数组合：沙蓬种子之间恢复系数为0.529、沙蓬种子之间滚动摩擦系数为0.057、沙蓬种子与ABS板间静摩擦系数为0.629。最后,将最佳参数组合的值进行仿真堆积试验验证,得到仿真试验堆积角为46.740°,与物理堆积试验相对误差为1.62%,表明最佳参数组合选取合理,可为沙蓬种子播种机械离散元仿真参数设置及优化设计提供参考依据。     

胶园植保硫磺粉颗粒离散元模型参数标定

【目的】白粉病是影响橡胶树天然橡胶产量的重要病害,常使用硫磺粉进行防治。为确定胶园植保过程中喷施硫磺粉颗粒的离散元仿真参数,对硫磺粉粉体的离散元参数进行标定。【方法】通过BT-1000型粉体综合特性测试仪进行物理试验获取硫磺粉休止角,测量多次取平均值。以休止角为响应值,应用Plackett-Burman试验筛选出对休止角影响显著的3个参数,进一步基于最陡爬坡试验确定显著性参数的最优取值区间,并通过BoxBehnken试验建立并优化休止角与显著性参数的二阶回归模型,以物理试验休止角为目标值,对显著性参数寻优得到最佳组合。【结果】经过物理试验测得的硫磺粉休止角大小为48.56°。通过Plackett-Burman试验筛选出的3个影响显著参数为：硫磺粉-硫磺粉滚动摩擦系数、硫磺粉-不锈钢滚动摩擦系数和JKR表面能,其余影响不显著的参数取值：硫磺粉-硫磺粉恢复系数为0.2、硫磺粉-硫磺粉静摩擦系数为0.6、硫磺粉-不锈钢恢复系数为0.2、硫磺粉-不锈钢静摩擦系数为0.5。最陡爬坡试验确定的最优取值区间：硫磺粉-硫磺粉滚动摩擦系数为0.1～0.2、硫磺粉-不锈钢滚动摩擦系数为0.15～0.35、...     

西南区坡耕地紫色土离散元模型参数标定

土壤离散元模型是离散元法研究土壤与触土部件相互作用的基础。为了获得西南区坡耕地紫色土离散元模型,本文基于堆积试验,采用EDEM2020软件中的Hertz-Mindlin with JKR接触模型,以休止角为响应值,结合土壤本征参数实测试验,对模型参数进行标定。实测参数包括土壤密度、含水率、休止角、粒径分布、颗粒-颗粒静摩擦系数、颗粒-45钢静摩擦系数,利用Origin2021b软件图像数字化工具处理图像,得到实际休止角为34.11°,变异系数为0.452%。利用Design-Expert11软件设计了Plackett-Burman试验,对影响拟合休止角的9个参数进行筛选,得到4个影响显著参数,即颗粒表面能、颗粒-颗粒滚动摩擦系数、颗粒-45钢静摩擦系数、颗粒-45钢滚动摩擦系数。设计最陡爬坡试验进一步缩小模型参数取值范围,设计Box-Behnken试验,对模型进行分析优化,得到4个参数最优组合,即颗粒表面能0.134 J/m²,颗粒-颗粒滚动摩擦系数0.231,颗粒-45钢静摩擦系数1.045,颗粒-45钢滚动摩擦系数0.228。结果表明,仿真休止角与实际休止角相...     

揉碎玉米秸秆螺旋输送仿真离散元模型参数标定

为进一步揭示揉碎玉米秸秆饲料在螺旋输送器中的输送机理，提高揉碎玉米秸秆螺旋输送过程中离散元仿真研究所用参数的准确度，将揉碎玉米秸秆分为秸秆穰、秸秆叶、秸秆皮，采用物理试验和仿真优化设计相结合的方法对离散元仿真参数进行标定。首先，通过物理试验得到穰、叶、皮各项本征参数和接触参数的平均值和范围，将范围相对较大的接触参数通过Plackett-Burman试验利用其显著性进行筛选。结果表明，对仿真堆积角影响显著的因素有穰-叶滚动摩擦系数、皮-叶滚动摩擦系数、叶-叶静摩擦系数。对这3个显著性参数进行最陡爬坡试验，将其取值范围进行优化缩小，通过Box-Behnken试验建立堆积角与显著性参数的二阶回归模型，堆积角的目标值设置为物理试验测得的30.4°，寻优得到最佳的显著性参数的取值：叶-穰滚动摩擦系数为0.325，皮-叶滚动摩擦系数为0.377，叶-叶静摩擦系数为0.411，在此组合下仿真堆积角为30.77°。为进一步检验参数取值的准确性，通过t检验得到P>0.05，表明仿真与物理试验堆积角的值无明显差异，验证了最佳参数取值的可靠性。由此表明，应用上述各物理试验和优化试验来测定及标定离散元仿...     

典型蓖麻蒴果离散元参数标定与分析

蓖麻蒴果离散元参数难以通过试验或文献直接获得,为获得蓖麻蒴果脱壳过程离散元参数需进行参数标定。以堆积试验为依据,对典型蓖麻蒴果离散元参数进行标定。选用通蓖11号和品种1987进行试验,以蓖麻休止角为响应值,应用PlackettBurman试验对物料离散元参数筛选出显著因素。利用最陡爬坡试验逼近最优参数组合区间,设计中心组合试验,选取最优参数组合。结果表明:通蓖11号的显著因素为蓖麻-蓖麻静摩擦系数、蓖麻-蓖麻动摩擦系数与蓖麻-Q235钢静摩擦系数,3个系数最优组合值分别为0.59,0.04,0.33,虚拟试验休止角平均值与真实值比值为94.48%;1987品种3个系数最优组合值分别为0.44,0.02,0.14,虚拟试验休止角平均值与真实值比值为96.44%。对通蓖11号3个参数进行优化,优化后各系数为0.6,0.05,0.4,此时虚拟试验休止角平均值与真实值比值为98%,仿真程度提高;通蓖11号各系数均大于1987品种,其中蓖麻-蓖麻静摩擦系数占比为0.27、蓖麻-蓖麻动摩擦系数占比为0.60、蓖麻-Q235钢静摩擦系数占比为0.65。3个参数越大,脱壳时的最大速度越小,仿真时间越长,最大压力及脱壳率越大。标定后的参数可为蓖麻蒴果离散元建模及分析提供参考。     

稻谷颗粒模型离散元接触参数标定与试验

EDEM与FLUENT耦合模拟稻谷清选时,稻谷颗粒参数设置的准确性直接影响仿真结果的可信度。结合稻谷颗粒堆积角的实测试验与仿真试验,标定了稻谷颗粒模型间静摩擦系数、滚动摩擦系数2个主要接触参数。设计了堆积角形成装置,该装置采用坍塌法和注入法同时形成2种堆积角,可以减小非同时形成所引起的测量误差;创建了稻谷颗粒的离散元模型,结合堆积角的实测试验与仿真试验,建立了2个主要接触参数与2种堆积角之间的二元回归方程;以稻谷颗粒2种堆积角的实测结果作为目标值,对回归方程进行数值求解,得到颗粒模型间的静摩擦系数、滚动摩擦系数分别为0.433 6、0.161 5;设计了抛洒试验对标定后的参数进行验证,抛洒长度、宽度和近抛洒端堆积角等特征参数的实测结果与仿真结果的相对误差小于5.7%,这表明颗粒的模拟运动轨迹与实际运动轨迹基本一致,所创建的稻谷颗粒离散元模型及标定得到的接触参数可以为EDEM-FLUENT耦合仿真稻谷清选过程提供参考。     

离散元法中所用土壤参数测量及标定方法研究

[目的]针对离散元仿真软件涉及的土壤参数进行测量及标定。[方法]基于Hertz-Mindlin with JKR黏结模型，通过直接测量法测量土壤的固体密度、弹性模量和泊松比，并用堆积角和滑动摩擦角来标定土壤接触参数。通过中心组合试验，采用Design-Expert 8.0.6软件，以土壤休止角、土壤与65Mn钢滑动摩擦角的仿真值与实测值的相对误差为优化目标进行回归分析。[结果]通过分析获得最优的离散元接触参数组合为土壤间恢复系数0.28、静摩擦系数0.49、滚动摩擦系数0.24、土壤表面能0.04 J/m²,土壤与65Mn钢间恢复系数0.59、静摩擦系数0.67、滚动摩擦系数0.13。在所测土壤参数及最优标定参数下，采用离散元仿真模拟探针入土行为，获得探针在8 mm/s的贯入速度下，贯入20、40、60、80和100mm处仿真试验和土槽试验探针阻力相对误差分别为8.59%、9.88%、9.72%、0.15%、6.98%,误差在可接受范围内。[结论]参数测量和标定方法准确可靠性，可为松软土壤的离散元仿真提供参考。     

基于离散元法的柔性片烟建模及仿真参数标定

为确定片烟在离散元仿真模拟中的合理模型参数取值，更好地计算出片烟在松散润叶设备中的复杂运动，基于离散元软件（Rocky DEM），创建了与以往球形黏结刚性颗粒模型不同的柔性片烟壳体模型，并分析了柔性片烟壳体模型形变特点。以物理和仿真试验的堆积角为对比对象，选取“Hysteretic Linear Spring”法向接触模型、“Linear Spring Coulomb Limit”切向接触模型以及“Constant Adhesive Force”黏附力模型，对仿真接触参数进行标定。再通过对片烟本征参数测定、查阅文献及Rocky DEM说明手册，分析得到了片烟的离散元参数取值范围，设计Plackett-Burman、最陡爬坡以及Box-Behnken试验，确定了影响柔性片烟堆积角的显著性参数以及堆积角的拟合方程。以堆积角物理试验35.83°为目标值，得到显著性参数的最优解：片烟与片烟间静摩擦系数为0.67，片烟与片烟间滑动摩擦系数为0.31，片烟与片烟黏间黏附力系数为0.2，其余非显著性参数取平均值。仿真试验与物理试验相对误差为0.78%，从而验证了模型的准确性与真实性，可为工业生产进一...