麦玉轮作区小麦种床制备工艺及床体成型-播深调节-镇压塑形装置研究

白圣贺 · 工学博士（机械设计及理论 · 中国农业大学，导师刘立晶） · 评审提交版

一、总体评价

本文面向华北平原一年两熟麦玉轮作壤土区小麦精量播种中种床不实、沟内落种位置不定、覆土镇压向内挤压致使播深稳定性与播种均匀性偏差的问题，研究高质量种床制备工艺及配套的床体成型、播深调节、镇压塑形成套装置。作者从镇压沉陷与多源切削阻力的解析建模入手，依次完成旋耕切削抛土装置、播深调节装置、播前镇压辊三套部件的设计与离散元仿真定参，再集成为整机并配播深闭环控制系统，最后在田间做对比与速度梯度试验。研究给出一组最优作业参数及覆土深度合格率、变异系数等量化指标，并报出较传统模式亩增产35kg（增幅5.6%）的田间结果；论文实验工作较完整，工艺路线可行，提供了可查的定量基础。选题属保护性耕作与精量播种装备的工程方向，对麦玉轮作区种床制备配套机具有一定参考价值；论文章节完整，工艺机理、装置设计、仿真定参、整机集成、田间验证一一对应，论点链不断裂，工作量饱满。论文结构完整，工艺机理、装置设计、仿真与田间各章衔接清楚，图表、公式与中英文摘要总体可循，但在关键公式数值复现、章节内部一致性、仿真向田间外推的口径、对照与统计严谨性以及符号与标准号规范等方面存在若干不宜忽略的问题，部分已进入摘要与结论，答辩前须较多修改完善。

二、主要问题与修改建议

问题1

式（3-75）给出的镇压沉陷量上限值与所列公式、参数对不上，且该结论已进入摘要与结论。第3章按光滑圆柱辊推出沉陷公式hq并取Dq=0.3m、q=800至2200N/m、k0=1.8×10⁵至4.2×10⁵N/m³，最小工况标30.1mm可独立复现，但最大工况按式（3-75）与所列参数代入算得约103.9mm，正文及摘要（3）、结论7.1（4）却标49.9mm（区间30至50mm），相差约一倍，且所列区间内任何参数组合都得不到49.9mm；另外实物为均布刀齿式镇压辊，与光滑圆柱接触模型力学本质不同，论文未给等效或修正说明。若上限工况确按q=2200、k0=1.8×10⁵取，则该被写入摘要与结论的设计落点偏离约一倍，下游辊体单位线载荷论证亦随之失据。建议作者：1. 复算式（3-75）的hq上限并订正数值，或交代上限工况实际所取的q与k0且与所列区间自洽；2. 同步修订摘要（3）与结论7.1（4）的30至50mm；3. 补充刀齿辊与光滑圆柱辊的沉陷等效或修正说明。

问题2

式（3-62）抛土匹配系数kQ的定义式两端数值对不上，恰好相差一个效率系数ηt。该式写kQ=Qb/Qa=ηt+Sb0·L/Qa约等于1.428，但按定义左端Qb/Qa=Sb/Sa=0.026/0.045=0.578，右端ηt+Sb/Sa=0.85+0.578=1.428，两端差恰为ηt=0.85，等号不成立；要使式（3-59）的Qb=Sb·L与式（3-62）隐含定义同时成立须Sb0=Sb−ηt·Sa=−0.012m²即负面积，且Sa大于Sb本应前轴抛土多、结论却为后/前=1.428。kQ究竟取0.578还是1.428取决于Qb是否计入接收的前轴料，两条定义并存使支撑上虚下实分层重构、并为第四章播深调节提供依据的关键比值失去自洽来源。建议作者：1. 厘清Qb定义并二选一订正——若Qb含接收前轴料则改写式（3-59）并明确Sb0含义与取值，若Qb=Sb·L为准则订正式（3-62）右端与1.428；2. 补充ηt=0.85与κ约3.3的来源，说明捕获比例大于1的物理含义。

问题3

第四章响应面优化的因素编码、目标函数与方差分析存在多处不自洽，削弱最优参数的论证过程。这些不自洽集中在第四章方法骨架，读者据此复核最优参数如何被推出会受影响；惟回归方程用正确中心编码能复现预测值，最优值602/323/28大概率不变，问题在表述与方差分析而非优化整体不成立。

1. 其一，因素编码与单因素试验对不上。表4-5将转速0级编为300r/min，与单因素三表的中心330r/min不一致，并使步长20与80不等距、字面上成为非法编码；而用330代入回归式才能复现预测值，反向印证0级应为330。
2. 其二，优化目标函数方向写反。式（4-28）将目标写为Qc取最小、Us取最大，与上方文字最小变异系数、最大合格率方向全反。
3. 其三，单因素与响应面趋势描述相反。单因素表4-3、表4-4显示Qc随转速、角度单调下降，与响应面文字先升后降相反，而最优解却取内部点。
4. 其四，方差分析平方和加不齐，但未必构成矛盾。表4-7各源项平方和之和与模型平方和差约4%与8.7%；若软件报告的是Type III偏平方和，则因二次项不正交本不必加和到模型平方和，故此处不径下矛盾之断。

建议作者：1. 订正表4-5转速0级为330，或说明非等距设计依据；2. 改正式（4-28）目标方向为Qc取最大、Us取最小；3. 核对表4-7平方和类型并标明，若为Type I应复核是否加齐、若为Type III宜在表注说明；4. 弥合单因素与响应面的趋势表述。

问题4

第五章镇压辊齿高、入土深度与防跳振速度多处自相矛盾，并违反本章自推约束。齿高与入土深度是镇压辊核心几何，临界速度更是高速作业稳定性结论，取值打架会让读者据此做出错误判断。

1. 其一，所取齿高超过本章自推的防挑土上限。式（5-8）防挑土须ht小于R·δ/(1−δ)，按δ取8至15%、R=133.5mm算得上限11.6至23.6mm，论文最终取ht=42mm，超上限约2至3.6倍，等价λ=1.12至1.21大于1，落入式（5-6/5-7）定义的齿尖反向切削、挑土工况。
2. 其二，入土深度按公式与正文标注对不上，并暴露齿高两个取值。式（5-13）按H=hc+ht应为62mm，正文却标52mm，反推H=52得ht=32mm，与p.127刚定的42mm冲突。
3. 其三，防跳振临界速度在正文与小结相差六倍多。5.3.4节正文6.8m/s可由式（5-45）复现并与图5-12量程一致，5.5小结却写42.9m/s（合154km/h，播种机不可能），二者差约6.3倍，而较大值进了小结作为对外结论。

建议作者：1. 复核ht=42与防挑土上限的关系，或说明实际滑移率达约24%以上的依据并与δ=8至15%弥合，或调整ht与δ使λ小于1；2. 统一ht取值并据此订正H及齿圈轴向厚度bt的流量守恒复算；3. 以正文6.8m/s为准订正小结42.9m/s，并标注式（5-45）长度统一取米的量纲约定。

问题5

表4-8的最优参数验证结果在仿真与田间之间口径不一，相对误差分母不统一，且标定与验证工况发生迁移。验证数据性质未定，直接关系到最优参数是否已被真实田间验证这一外推效度能否判定，分母两套口径削弱误差小于5%的可靠性论证。

1. 其一，同一张表的数据是仿真还是田间前后不一。表4-8标题为仿真验证结果、列名仿真次数，正文却叙述为田间验证3次重复，另一处又称与田间试验相比仍存在误差。
2. 其二，两个指标的相对误差用了不同分母。Qc的1.97%以实测均值作分母、Us的2.99%以预测值作分母，文字却称与预测值的相对误差，口径自相矛盾，惟两值仍均小于5%、头条结论不被推翻。
3. 其三，标定与验证工况发生迁移。最优解在菏泽标定（含水率14.5%、坚实度单位kPa），整机在临邑验证（含水率17.3%、紧实度单位N/cm²），地点、含水率、单位均不同，且紧实度单位全篇不统一。

建议作者：1. 明确表4-8数据性质并订正标题与列名使其与正文一致，若仿真与田间皆有则分表给出；2. 统一相对误差分母（建议一律用预测值）并复算Qc、Us两项；3. 说明菏泽标定最优解迁至临邑验证的工况可比性，或补同工况验证，并统一全篇紧实度单位。

问题6

第六章对比与速度梯度试验的对照公平性与统计严谨不足。该章田间试验是全文唯一的真实田间证据，承载增产、合格率、容重分层等对外结论；对照机型号与调参不明使优于传统模式的比较基础不透明，缺样本量、检验方法与速度间检验使增产5.6%等差异的统计可信度无从核验。

1. 其一，对照公平性不足。作对比的传统播种机未给型号，其播深与转速等调参状态亦未控。
2. 其二，统计交代不足。对比与速度梯度试验均仅重复3次，无样本量n标注、无检验方法说明，速度之间无显著性检验，趋势全凭文字定性。
3. 其三，显著性阈值标注前后冲突。正文称株高整齐度极显著（p<0.05）、表6-5标p<0.01，且极显著配p<0.05本身矛盾。
4. 其四，提升幅度一列混了三种口径。表6-5中出苗率、均匀性、株高是百分点差却标成%，茎秆强度+1.6%实为绝对差+1.6N·mm⁻²，仅产量+5.6%是真正的相对%。

建议作者：1. 补充传统对照机的型号与调参状态，说明两机各自调至适宜工况的公平性；2. 补样本量n、检验方法与速度间显著性检验，明确统计单元与重复定义；3. 统一显著性阈值与措辞（极显著对应p<0.01）；4. 统一表6-5提升幅度列口径并订正茎秆强度。

问题7

第二章沉陷与多源阻力解析模型的几个核心常数只给符号、无取值和来源，式（2-31）的速度修正项无出处。第二章建立的镇压沉陷与多源阻力模型是后三章装置参数的源头；若关键常数取值与来源缺失，该解析模型便无法被独立复算、向后三章的参数传导也难追溯，（1+λv）依据不明则该动态修正缺出处。

1. 其一，核心常数只给符号。kc、kφ、n、λ、Kp、Kpc在第二章仅给符号定义，无取值、无数据来源，k0要到第三章才以《Terramechanics》加课题组实测给出区间。
2. 其二，速度修正项无文献出处。式（2-31）的Bekker主体式与标准式一致，但其中（1+λv）线性修正在此处无任何文献出处，公开文献未见同形对应，更像作者自设。
3. 其三，2.5节未定量回算三大机理模型。该节工艺可行性验证仅给定性运移、平整度18mm、秸秆比1∶1.24∶0.40，未回算2.4节沉陷量hq、未定量验证2.3节播深映射。

建议作者：1. 在第二章给出kc、kφ、n、λ、Kp、Kpc的取值与来源，或指明其数值首次出现的章节与数据出处；2. 在式（2-31）处补（1+λv）的来源文献，或声明为自设并给出λ的标定方法、量纲依据与v的取值范围；3. 说明2.5节仿真是否、如何覆盖沉陷量hq与播深映射的定量验证。

问题8

多个符号跨章代表不同物理量、同一物理量跨节又换符号，影响公式与参数的可追溯性。全文存在多处符号漂移：δ三义（旋耕刀根部截面厚度10mm、滑移率8至15%、刮削角30°），γ三义（正切刃楔角35°、刀齿楔角55°、入土包角约90°），θ两义（相位角、悬臂下摆角31.6°），ω两体（旋耕角速度、镇压辊角速度），R符号族最乱、前后旋耕回转半径Ra/Rb与辊体半径R=133.5、齿尖半径Rt、镇压比R′并存，且镇压辊半径同时写作Ra、Rc、Rq；第六章另有Hr、Dr一符两义与Ht、Ha混用。同一符号跨章代表不同量、同一量跨节换符号，读者在跨章核对公式与参数时容易认错。建议作者：1. 通篇核查并统一符号，为复用符号加章节下标或换字母（如镇压辊半径统一为Rq、滑移率改用s）；2. 在论文前加主要符号表，标明各符号的章节归属与定义，重点厘清R符号族与δ、γ的三义。

问题9

旋耕刀引用的国家标准号用错，应为GB/T 5669-2017。论文在第三章旋耕刀选型及第四章相关段多处把规格依据引为GB/T 19901-2017《旋耕机械旋耕刀》，但经国家标准全文公开系统核对，GB/T 19901实为《温度计检测元件的金属套管 实用尺寸》、库中仅有2005版、无2017版，与旋耕刀无关；旋耕刀和刀座的现行国标是GB/T 5669-2017《旋耕机械 刀和刀座》。关键选型依据引错标准号，读者按GB/T 19901查到的是温度计套管标准，会影响规格依据的可追溯性与规范性。建议作者：1. 把全文GB/T 19901-2017《旋耕机械旋耕刀》更正为GB/T 5669-2017《旋耕机械 刀和刀座》，涉旋耕机整机另引GB/T 5668-2017《旋耕机》；2. 将弹簧钢GB/T 1222-2007更新为现行GB/T 1222-2016（65Mn的三项力学值正确、无须改动，仅更新年版）。

三、次要问题

以下为校对、记号、数值标签、术语、单位与口径类问题，建议随全文通校一次订正，不与上述方法问题等长。代表例如下：公式与记号方面，式（2-22）抛物线末项分母漏平方（对照同构式（2-20）有v0x²、掉了上标2），式（2-17）、（2-27）、（2-28）根号内y差的符号方向不一、v′py与v′Py大小写混用；数值与标签方面，式（3-53）邻近前刀轴旋转周期Ta=0.214s属值对名错（0.214s实为后刀轴周期），起始极径ρkmin文述180mm但式（3-33）实代160mm，式（6-17）y0=−390mm与文字机架下方340mm差50mm；术语与单位方面，变异系数被第四章正文一处称稳定性系数（语义相反），上虚下实、上松下实、上松下紧混用，沉陷量、下陷量、下沟量混用，Box-Behnken与中心组合回归、二次正交旋转组合混用，kg/亩与kg/mu混用，显著性p与P大小写混用，紧实度单位kPa与N/cm²全篇不统一。建议作者全文通校并统一符号、单位、标准号引用，并核校摘要、结论与正文之间的数值、量纲与命名一致性。

四、评审结论

较大修改后答辩。依据：论文工作量饱满、工艺路线可行、论点链完整，所列问题以数值矛盾、章节内部不一致、仿真向田间外推口径、对照与统计严谨性、符号与标准号规范为主，其中沉陷量上限、覆土深度合格率与变异系数、防跳振速度等已进入摘要与结论或属核心装置参数，但这些问题多可通过复算、勘误、补充说明与统一口径修改完善，属可修改、非根本性问题，故定为较大修改后答辩，须在答辩前较多修改完善。

附：凡科博士 rubric（如走凡科系统填此表，否则可忽略）