← 返回总览王龙宝 · 详细佐证

详细佐证(详细评审意见)

论文:《基于单粒分离-稳定投种的小麦宽幅匀播技术与装置研究》
作者:王龙宝 培养单位:中国农业大学 学位:工学博士 专业:农业机械化工程


一、总体评价

本文面向现行小麦宽幅播种由外槽轮排种器配合开沟器撒播作业后种沟内种子分布均匀性差、断条漏播重播突出的问题,提出基于单粒分离与稳定投种的宽幅匀播方法,研制了双气室小麦精量排种器与气流附壁式投种装置。作者建立小麦种子多面体离散元模型并以休止角为目标标定接触参数,对排种与投种各过程作理论受力分析,搭建 CFD-DEM 耦合仿真,并完成排种台架、投种台架与整机田间四级试验。结果给出各环节最优参数组合(排种最优单粒率80.7%、投种最优粒距合格率76.7%、分布变异系数33.1%、田间一类苗占比86.7%),核心装置获授权发明专利与实用新型专利支撑。论文章节结构完整,理论-仿真-台架-田间链条清晰,图表与中英文摘要较规范,工作量饱满,工艺路线可行,具一定工程参考价值。论文也存在若干不宜忽略的问题:理论播种粒距5cm这一设计基准与式(2-1)回算结果不符并传至下游结构参数;标定响应面最优落约束边界、预测误差低于全部实测点;三章寻优最优解均落设计空间边界、验证试验为重测已测点;投种达标声明混用了规范无阈值的自定指标;田间核心对照四因素同时变化致归因不可分离、且全程无统计检验。上述问题多属数值核对、口径界定与试验设计完善,可通过修改与补充说明解决,答辩前须较多修改完善。


二、主要问题与修改建议

意见1:理论播种粒距5cm与式(2-1)回算的8.63cm不符,并传播为全文设计基准

定位:式(2-1)(印p22);式(3-4)型孔数Z=36(印p42);第4章台架“保持理论播种粒距为5cm”(印p100)

问题:第2章用式(2-1)推算理论播种粒距,正文写“最小理论播种粒距约为5cm”,但把论文自给的同一组输入数据代回同一印刷公式,得到的是8.63cm,二者相差约73%;这个5cm 又被式(3-4)与第4章台架当作硬基准沿用。

证据

同一个式(2-1)、同一组输入,正文结论与回算结果对不上。

  • 式(2-1):Ns=1000²×Q/q,Ss=666.7×4×10000/(Sr×Ns)
  • 输入(正文自给):Q=9kg/亩、q=43.7g(鲁原502千粒重)、Sr=15cm
  • 回算:Ns=205950粒/亩,Ss=666.7×4×10000/(15×205950)=8.63cm
  • 正文却写:约5cm
  • 差距:8.63cm 对 5cm,约73%

为什么是问题:这个5cm 是后续结构参数链的起点,起点的数与公式对不上,整条链就跟着不确定。式(3-4)用 Ss=5cm 定出每列型孔数Z=36(回算35.9,与设计值一致),第4章投种台架又把“保持理论播种粒距5cm”作为试验前提;由于式(2-1)结果、式(3-4)的Z=36 与第4章5cm 基准三者共用一个5cm,而该5cm 与式(2-1)直接回算的8.63cm 不符,型孔数与投种粒距判据链的起点随之存疑——若 Ss 取8.63cm,则Z≈21,排种轮型孔数与列间距结构需随之改变。这是同一公式、同一数据自相差73%的内部不一致,不需要任何外部假设即可看出。

修改建议


意见2:休止角标定的响应面最优落约束下界、预测误差低于全部实测点

定位:式(2-20)优化回归(印p37);式(2-21)约束与表2-7 BBD结果(印p37-39);优化结果D=0.15/E=0.11694/H=0.0797、预测δ=1.96%(印p39)

问题:第2章以休止角相对误差δ最小为目标,在式(2-20)回归模型上求最优接触参数,得到一组最优组合并称其预测误差为1.96%、验证仿真2.3%。但该最优组合中摩擦参数D 落在约束区间的下边界、不是模型内部的驻点,且预测的1.96% 低于全部15个BBD实测点的最小实测误差2.4%。

证据

最优点位置与预测误差水平,均落在数据约束之外。

  • 式(2-20)对D 只有一次项、无D²项,对D 是线性的
  • 最优 D=0.15,恰为式(2-21)约束区间[0.15, 0.275]的下边界(非内部驻点)
  • 预测δ=1.96%
  • 全部15个BBD实测点的最小实测误差=2.4%(run7)
  • 即:预测值1.96% < 全部实测点最小值2.4%

为什么是问题:一个优化“最优”如果落在因素水平的边界上,往往提示所设的因素区间没有覆盖到真正的最优区,而不是在区间内找到了驻点;预测误差低于所有实测点,则说明这个“最优”是外推到了数据之下。由于式(2-20)对D 无二次项、响应对D 单调,极小被推到约束下界D=0.15,因此最优落边界是模型结构决定的几何事实,而预测1.96% 低于全部实测点的2.4% 同样来自论文自身的两个数。至于式(2-20)回代得1.19%、与论文1.96% 不一致这一点,可能源于 Design-Expert 求解时使用了未截断的全精度系数、而印刷公式只保留了2至3位有效数字;若作者给出全精度系数或复算过程,这一项或可解释。上述差异使第2章“响应面优化得到可靠最优参数”的论证强度有待加强,但不直接证伪标定参数本身。

修改建议


意见3:三章响应面最优解均落设计空间边界,验证试验为重测已测点

定位:式(3-48)排种寻优与表3-4第11组(印p58);式(4-36)投种寻优与表4-7第14组(印p109);第2章式(2-21)边界最优

问题:第2、3、4章都用BBD响应面求多目标最优并各做一次验证试验。三章的最优解分别落在设计空间的边界顶点或棱边上,而非内部驻点;其中第3、4章的最优组合恰好等于某一已经测过的BBD试验行,所谓“验证试验”实际是对已测边界点的重测,而不是对优化新解的独立校验。

证据

三章最优解的编码位置与对应的已测试验行。

  • 第2章:最优D 编码=−1(约束下界)
  • 第3章:最优编码(X1=+1, X2=0, X3=+1),回代回归式得Y1=81.29%、Y2=4.52%、Y3=14.16%(与论文吻合);该编码=表3-4第11组(实测80.3/5.2/14.5),X1、X3 均在+1边界
  • 第4章:最优编码(x1=−1, x2=−1, x3=0),回代得y1=76.03%、y2=8.41%、y3=15.50%、y4=32.64%(与论文吻合);该编码=表4-7第14组(实测75.8/8.3/15.8/32.8),x1、x2 均在−1边界
  • 第3章“验证试验”80.7/5.8/13.5、第4章“验证试验”76.7/9.2/14.2/33.1 均为对应已测点的同点重测

为什么是问题:响应面优化的目的,是在所测试验点之外的设计空间内部找到一个更优的工作点,再用一个新的、不等于任何已测行的工况去验证它,验证才具有独立校验的意义。三章的最优都落在因素水平的边界上、且第3、4章最优恰为已测行,使“响应面优化获得最优工作参数”这一贯穿三章的方法学结论缺少独立验证支撑;若真实最优确在边界,则提示三章的因素水平区间可能都没有覆盖到最优区。这是全文优化环节系统重现的同一类特征,不是某一章的偶然。

修改建议


意见4:接触参数标定以“10次取极值”定因素水平、休止角仅单均值,输出却给5位有效数字

定位:本征与接触参数测定“重复10次试验取极值”与表2-2(印p27-29);休止角均值32.9°(印p29);优化输出动摩擦0.11694、小麦-树脂动摩擦0.0797(印p39)

问题:第2章测定杨氏模量、泊松比、堆积密度、各类摩擦与碰撞恢复系数时,均以“重复10次试验取极值”的上下限直接作为PB试验的高低水平;休止角作为唯一寻优目标只给了均值32.9°、未给离散度;而标定输出的动摩擦系数却给到5位有效数字。输入的精度交代方式与输出的有效数字位数不相称。

证据

为什么是问题:用极值上下限作因素水平,区间宽度会被极端样本支配,得到的不是均值附近的常规波动范围;而休止角只有一个均值、没有标准差,δ=2.3% 这个误差相对实测本身的分散度是大是小就无从判断。由于标定结果是第3、4章全部CFD-DEM仿真的接触参数输入,输入端用极值定区间、输出端却报到5位有效数字,使“参数标定可靠”这一结论的论证完整性有所欠缺。这是论证交代方式的问题,并不直接否定标定数值本身。

修改建议


意见5:投种“各项指标满足NY/T1143-2006规范”混用了规范无阈值的自定指标

定位:表4-5(引NY/T1143-2006,印p101);式(4-31)分布变异系数定义(印p101);4.5本章小结(4)(印p111)

问题:第4章投种台架对标NY/T1143-2006,表4-5 列出该规范理论粒距≤10cm 档的四项阈值,其中第四项是“合格粒距变异系数≤40%”;论文 p101 自己注明该项需去除重播与漏播粒距。但论文用于评价的 y4 是“分布变异系数”(含全部粒距、二维),最优值33.1%,与规范第四项不是同一个量;小结却称投种“各项指标满足规范要求”。

证据

为什么是问题:“满足规范”要成立,前提是所报指标与规范规定的指标是同一个量、按同一口径计算。合格率、重播率、漏播率三项确可对标≤10cm 档且达标,但33.1% 是含全部粒距的分布变异系数,规范规定的是剔除重漏后的合格粒距变异系数,两者口径不同,规范对前者并无阈值。把33.1% 计入“满足规范”,使达标声明的覆盖范围被高估——严格说只覆盖前三项;而“满足规范”又上行到摘要与结论的“稳定投种可行”。这一口径界定问题由论文自身文本即可确认,不依赖外部引文能否取得。

修改建议


意见6:田间“匀播与撒播”对照四因素同时变化,核心主张归因不可分离

定位:样机配置(印p112-113);播种试验小结(2)(3)(印p128);撒播侧镇压轮驱动(印p124)

问题:第5章田间以“宽幅匀播与宽幅撒播”对比作业质量,结论称匀播种沟内分布均匀性相对较高。但两组之间,开沟器类型、排种器类型、有无投种装置、驱动方式四者同时不同(匀播为芯铧开沟器加双气室精量排种器加气流附壁投种装置加电机驱动;撒播为圆盘或翼铲开沟器加外槽轮排种器加镇压轮驱动加软管撒播),没有任何一个只切换排种或投种系统、其余固定的解耦对照单元。

证据

为什么是问题:要把“种沟内分布均匀性相对较高”归因于本文核心的排种与投种系统,需要把开沟器、驱动方式等其它变量固定下来、只切换排种或投种系统。第5章四个变量同时改变,使均匀性的改善无法单独归因于排种或投种系统,它也可能来自芯铧开沟器或恒速电机驱动。第5章是全文核心主张(单粒分离-稳定投种提高种沟内均匀性)的整机验证环节,归因不可分离会削弱该主张在整机层面的实证支撑。

修改建议


意见7:田间播深合格率在0.9m/s时73.3%未达标,小结只呈现有利指标

定位:播深合格率随速度83.3→81.7→73.3%(印p123);播深合格率“核心指标”自述(印p121);播种试验小结(2)(印p128)

问题:第5章宽幅匀播播深合格率在0.6、0.75、0.9m/s 三个速度下为83.3%、81.7%、73.3%,论文 p121 自述播深合格率是评价核心指标、p123 自认0.9m/s 时未满足行业标准(NY/T996-2006 阈值≥80%)。但同时圆盘式撒播在播深上整体优于匀播;小结(2)只提匀播播幅与播种均匀性变异系数相对较低,未把播深合格率随速度递减、0.9m/s 不达标、且反劣于撒播的事实写入对比结论。

证据

为什么是问题:播深合格率是论文自述的核心评价指标,0.9m/s 已低于行业阈值、且在该指标上撒播反而更优,但这两点都未在小结体现,“满足标准”也未限定速度范围,读者会以为整体可行。由于核心质量指标在高速档不达标却未在结论呈现,第5章“提高播种质量”的总主张与播深这一关键指标的实测水平之间出现张力。

修改建议


意见8:田间试验全程无方差分析、显著性检验与重复离散度

定位:播种试验“每组重复3次”(印p121);植株破坏性取样n=60(印p125-126);田间各结果(印p123-128)

问题:第5章田间播种质量各指标每组重复3次、植株生长每模式取60株,给出一类苗占比86.7%/48.3%/53.3%、地上部干物质0.43/0.36/0.32g 以及株高、分蘖、叶龄等组间差异并直接作结论。但全章没有方差分析、显著性检验或置信区间,也未给3次重复的离散度,而台架试验(第3、4章)是配有ANOVA 的。

证据

为什么是问题:组间差异是否超出试验误差,需要方差分析、显著性检验或置信区间来支撑;田间是核心验证环节,却在统计推断上弱于台架,读者无法判断一类苗占比、干物质等差异是真实差异还是落在试验波动之内。这使“苗齐苗壮”“分布均匀性相对较高”等田间结论的证据强度低于台架。

修改建议


意见9:排种与投种台架最优性能偏低且无部件级对照基线

定位:排种最优单粒率80.7%、漏吸5.8%、重吸13.5%(印p73);投种最优分布变异系数33.1%、粒距合格率76.7%(印p109)

问题:排种台架最优单粒率80.7%、投种台架最优分布变异系数33.1%与粒距合格率76.7%,均为各自最优参数组合下的结果。两台架都没有与外槽轮撒播或其它精量排种器、投种方式的部件级对照;80.7% 只与自身回归预测81.29% 比较,33.1% 只与自身回归预测32.64% 比较。

证据

为什么是问题:“精量”“稳定”“提高均匀性”都是相对性的判断,要成立需要一个可比的基线(如外槽轮撒播或同类精量排种器、投种方式);两台架的最优指标只与自身回归预测比较,属自比而非他比。部件级缺对照,使排种与投种环节的相对性结论无法单独建立,而整机环节又有意见6 所述四因素混淆,核心机制的有效性主要靠自比支撑。

修改建议


意见10:投种仿真所示均匀性(CV约3%)与台架实测离散度(33.1%)量级悬殊,外推未说明

定位:CFD-DEM仿真投种速度CV 3.1~3.9%(表4-2/4-3,印p92);台架最优分布变异系数33.1%(印p109);仿真“8/9/10号型孔压力突变致漏吸”(图3-23/3-25,印p61-62)

问题:第4章CFD-DEM仿真(单管50粒)给出投种速度变异系数仅3.1至3.9%,同一装置的台架最优分布变异系数却为33.1%,二者量级相差约10倍。第3章还有“隔板平板扰流致8/9/10号携种区型孔总压突变、漏吸增大”这一来自仿真的机理结论,但台架以整体单粒率、漏吸率、重吸率为指标,未对特定型孔列做分列统计印证。

证据

仿真与台架的两个离散度数值,度量对象并不相同。

  • 仿真投种速度变异系数:3.1~3.9%(表4-2/4-3,单管50粒投种速度的离散度)
  • 台架分布变异系数:33.1%(种沟内粒距分布离散度,含重播漏播、二维)
  • 量级差:约10倍
  • 注:表4-1至4-4 各数值回算均自洽,无算术错误

为什么是问题:仿真所示“投种速度很均匀”与台架所示“种沟分布较离散”度量的不是同一个量,因此10倍的差距本身不构成数值矛盾,但仿真结论要向台架、田间外推,需要以同口径指标或对照来说明其预测力。由于CFD-DEM仿真是第3、4章机制论证的主要手段,而仿真所示的均匀性与特定型孔漏吸都未被台架同口径指标隔离印证,若不补充说明,仿真结论向台架的外推效度尚不充分;这一项指向外推说明的补充,并不否定仿真本身的正确性。

修改建议


三、次要问题

以下为校对、勘误、取值交代与完善性问题,逐条不单列,集中说明,建议作者全文通校核对。校对与勘误方面:第3章 p64 段79 变量名标注串行,“分别将种层高度设置为2、4、6、8cm,并将种层高度设置为中间水平27.5r/min”,其中被赋转速值的应为“排种转速设置为中间水平27.5r/min”;式(3-48)寻优约束“−3kPa≤X1≤−5kPa”数值上下限方向异常(−3大于−5),应为“−5≤X1≤−3”或注明负压绝对值口径;表3-2/3.3.3 引用“固相参数如3.3.21章节所述”疑为“3.3.2.1”编号脱点,内容回勾对应DEM模型节;同一量两名与命名不统一两处,“休止角”与“堆积角”(第2章正文与2.3.4.1标题并存)、“气流附壁式投种装置”与“气力式小麦投种装置”(题名摘要与结论6.2(1)并用),建议全文各取一名统一。取值依据与试验交代方面:型孔直径2mm 落经验式 d=(0.64~0.66)W=2.41~2.48mm 的下界(取2 而非中值2.45 的依据未给,与自述“孔径过小致漏吸”风险方向一致),投种通孔Φ1=11mm 非式(4-1) 1.8d′=8.80mm 的取整值(取整应为9mm,11mm 高于下限约25%,当量直径取谐均、系数1.8 无引文),建议补两处取值依据;BBD编码中心与单因素中间水平两处不一致(排种BBD中心−4kPa/25/6cm 对单因素中间−3.5/27.5/5,投种BBD中心0.75/9/2 对单因素中间0.8/8/2.5),影响“中心点重复即纯误差”的代表性,建议统一或说明;投种台架工作气压固定取最低值2kPa、依据“投种速度具随机性”交代偏笼统,而仿真示气压2 至8kPa 投种速度2.74 至4.12m/s(差约50%),建议补充固定取最低气压的依据并说明气压与转速对投种速度的耦合;投种最优作业速度0.6m/s 即最低水平、合格率随作业速度单调降(回归一次项主导),与宽幅播种高效作业取向之间的权衡未在文中讨论,建议补一句权衡说明;投种台架共测120个粒距(排种台架为600个型孔),未说明分几个测区、是否含3次重复的合计,建议补样本构成说明;越冬期匀播株高18.9cm 小于撒播21.3、20.7cm(匀播最小),小结结合叶龄、分蘖、干物质较高解释为“苗齐苗壮”,但株高较小是否为正面指标未给农学依据,建议补依据或调整表述;蒋杉(2025)“一器四行”既有气吸方案、丁启朔(2024)与王淞(2024)“二维分布变异系数”两处近新引文,因检索限制未独立核到原文(其存在性合理、不判作者有误),建议作者补全可核出处,并就创新点2 相对既有“一器四行”方案的差异点作一句辨析。


四、评审结论

较大修改后答辩。

依据:论文工作量饱满、理论-仿真-台架-田间链条完整、核心装置有授权专利支撑,工艺路线可行,不存在颠覆整体的硬伤;但理论播种粒距基准与公式回算不符并传播至下游、响应面最优落约束边界且预测低于全部实测点、三章寻优落边界且验证为重测、投种达标声明口径错配、田间核心对照四因素混淆且无统计检验等问题,分布于设计基准、标定可靠性、寻优方法、达标声明与核心验证多个环节,需作者较多修改与补充说明后方可答辩。

← 返回总览本篇提交版 →