私密分享 · 本报告为送审前自审参考,含未公开论文内容,仅供作者本人与导师查阅,请勿外传 · 不代表盲审或答辩委员会意见
学位论文自审报告
坡地播深精准调控方法与装置研究
E 研究内容与方法 / E2 数据可靠性与可溯源

E2 数据可靠性与可溯源

正文数据与图表是否一致、数据来源是否可查、有无查无对应的数字。

必须修改 3建议修改 5提醒 3
必须修改 E2-001

同一工况(θ=25°, s=30, v=5)振动加速度峰值在 §3.2.2 与 §3.2.3 仍给两套数(1561.97 vs 1264.90)

📍 位置物理 p65 / 印刷 p51(§3.2.2(2) 图3-5);物理 p72 / 印刷 p58(§3.2.3(4) 表3-9/图3-7 narrative)
问题
同一台坡地播种机、同一组参数(地表坡度 25°、弹簧初始增量 30 mm、作业速度 5 km·h⁻¹),振动加速度峰值在两节里仍报两个差很多的数,读者不知道哪个是真值。
为什么是问题
- 数据对比(独立成行): - θ=25°, s=30, v=5:§3.2.2 图3-5 = 1561.97 mm·s⁻² - θ=25°, s=30, v=5:§3.2.3 表3-5 序号 19(S19)经 §3.2.3(4) 复述 = 1264.90 mm·s⁻² - 差值 ≈ 297 mm·s⁻²(相对差 ≈ 23%) - 后果:两节参数完全相同却差 23%,盲审会质疑是同一仿真的两次抄录矛盾,还是两批独立仿真(单因素扫描 vs 中心组合设计)而正文未说明口径差异(随机种子、统计时间窗、稳定段选取不同)。MBD-DEM 为随机离散元模拟,批间波动可解释,但正文未交代,构成可溯源疑点。
修改建议
  • 复核两处 θ=25° 峰值的原始仿真记录,确认 1561.97 与 1264.90 各自来源。
  • 若属两批独立仿真,加一句说明二者为不同试验、取值口径不同;若属抄录/取数错误,统一为同一数值。
证据 / 原文锚点
  • 物理 p65「上坡时,随着地表坡度 θ 的增加,振动加速度峰值呈上升趋势,在地表坡度 θ = 25° 时达到振动加速度峰值局部最大值为 1561.97 mm·s⁻²」。
  • 物理 p72「当地表坡度 θ 逐渐增大时(S1、S3 和 S19)……在 θ = 25° 时,振动加速度峰值达到 1264.90 mm·s⁻²」。
必须修改 E2-011

控制模型只标定到坡度 15°,却用到 18.8°~20.1°;且两章把田间坡度上限写成不同的两个数(NEW)

📍 位置物理 p104 / 印刷 p90(式4-9~4-12 分段);物理 p100 / 印刷 p86(§4.3.1 田间坡度 0~20.1°);物理 p107 / 印刷 p93(§5.1 田间坡度 0~18.8°);物理 p109 / 印刷 p95(§5.2.2 出现 θ=19°)
问题
决定输出下压力的核心数学模型(式4-9~4-12)按地表坡度分段,分段区间只到 **0~6°、6~10°、10~15°**——也就是模型只在 |θ|≤15° 范围内拟合、定义。但田间试验和验证的坡度远超 15°:第四章 §4.3.1 说试验田坡度「0~20.1°」,第五章 §5.1 又说试验基地坡度「0~18.8°」,§5.2.2 正文还直接出现「地表坡度 θ 为 19° 时」的实测点。模型在 15° 以上没有对应的拟合段,却被用到 19°~20° 的工况,属把模型外推到标定范围之外;更直接的问题是同一片试验田的坡度上限,两章给了 20.1° 和 18.8° 两个对不上的数。
为什么是问题
- 范围对比(独立成行): - 控制模型式4-9~4-12 分段上限:|θ| ≤ 15°(最后一段 10~15°) - 第四章 §4.3.1 田间坡度范围:0~20.1° - 第五章 §5.1 田间坡度范围:0~18.8° - 第五章 §5.2.2 实测点:θ = 19° - 后果:① 19°~20° 工况下,控制器要么把模型外推(拿 10~15° 段的回归往外推 5°),要么没有对应段可用——正文未交代 15° 以上怎么算输出下压力,控制有效性在这段缺依据;② 同一试验田坡度上限 20.1° 与 18.8° 两章打架,读者无法确定真实坡度上界、也无法判断 19° 实测点是否在范围内;③ 这直接削弱「精准调控覆盖丘陵山区大坡度」的整体结论。
修改建议
  • 二选一:要么把控制模型分段补到实际用到的坡度上限(增加 15°~20° 段并给该段回归),要么在结论/方法里把模型与控制策略的适用坡度明确限定为 |θ|≤15°,并说明 15° 以上工况如何处理(外推/截断/另算)。
  • 统一第四章与第五章田间坡度范围上限(核对到底是 18.8° 还是 20.1°),全文一处取值;说明 §5.2.2 的 19° 实测点落在哪个范围内。
证据 / 原文锚点
  • 物理 p104「按照地表坡度类型分为上坡和下坡,地表坡度区间划分为 0~6°、6~10° 和 10~15°」(式4-9~4-12 分段只到 15°)。
  • 物理 p100「采用数显倾角仪……测定地表坡度 θ 范围为 0~20.1°」(§4.3.1)。
  • 物理 p107「该试验基地耕地地表坡度 θ 在 0~18.8° 范围内分布」(§5.1)。
  • 物理 p109「在输出下压力为 300 N 的条件下,地表坡度 θ 为 19° 时,实际下压力均值仅为 26.55 N」(§5.2.2)。
  • **另见**:数据范围/外推的诚信定性见 G_integrity.md 对应条目(互链,不重复合并)。
必须修改 E2-012

三章正文都改成「下坡下压力增大」,但第三章受力推导(式3-16/3-20)仍按 cosθ 推出「下坡减小」,公式与结论、与图3-14 数据三方打架(NEW)

📍 位置物理 p76–77 / 印刷 p62–63(§3.3.1 受力分析 式3-16/3-20/3-21);对照 物理 p80–81 / 印刷 p66–67(§3.3.4 图3-14 narrative)、物理 p102 / 印刷 p88(§4.3.3)、物理 p109 / 印刷 p95(§5.2.2)
问题
本轮把全文「下压力随坡度怎么变」的文字结论统一改成了「上坡减小、下坡增大」(并改用 |θ|),第三章图3-14 一节、第四章 §4.3.3、第五章 §5.2.2 三处叙述都一致,也和图3-14 实测形态(最大下压力在 θ=−20° 出现、上坡 θ=25° 最小)对得上——这部分改对了。但第三章 §3.3.1 的受力分析没跟着改:它给出的下压力公式式(3-16) 写的是 **F_D = G·cosθ + F₃·cosα**,cosθ 对上坡、下坡都是「角度越大越小」(cos 是偶函数),所以这条公式只能推出「上坡、下坡下压力都减小」。式(3-20) 处理下坡时也仍白纸黑字写着「下坡作业时……重力分力 G·cosθ 随 |θ| 增大而减小,导致对地下压力 F_D 减小」——结论是减小;可紧接着的下一句却又说「同时导致播深调节装置实际对地下压力 F_D 的增大」。同一页里,下压力公式说减小、下一句说增大,自己跟自己打架,而且和你已经改好的「下坡增大」结论也对不上。
为什么是问题
- 三方对比(独立成行): - 图3-14 实测数据:下坡 θ=−20° 下压力均值 921.6 N(最大)、θ=25° 640.1 N(最小)→ 下坡下压力大、上坡小 - 正文新结论(图3-14 一节 / §4.3.3 / §5.2.2):上坡减小、下坡增大(用 |θ|)→ 与图一致 - 受力推导式(3-16)/(3-20):F_D = G·cosθ + F₃·cosα,cosθ 随 |θ| 增大而减小 → 推出「下坡也减小」,且式(3-20) 明写「导致对地下压力 F_D 减小」→ 与图、与正文新结论相反 - 后果:① 受力分析是「为什么下坡下压力会增大」的机理依据,现在这条依据(式3-16/3-20)反而推出相反方向,等于结论没有受力模型支撑;② 同一页式(3-20) 先说 F_D 减小、后说 F_D 增大,自相矛盾;③ 「下坡增大」目前只靠式(3-21) 的前进牵引合力 F_Q / 前进阻力 F_Z 和「重心前移」一句话带过,而这套是水平方向受力,并不是 F_D 那条垂直方向公式——用水平力的变化去解释垂直下压力增大,逻辑链没接上;④ 盲审一旦顺着公式推,立刻能看出「公式给的方向和你结论、和你图都反着」。
修改建议
  • 重做下坡段受力推导,让 F_D 的表达式能解释「下坡 |θ| 增大时 F_D 增大」:要么补上 G·sinθ 经平行四杆几何传到垂直方向的分量、要么显式写出牵引角 β 随坡度变化如何反过来抬高 F_D,使式(3-16) 一类的最终 F_D 表达式与图3-14、与正文结论同方向。
  • 把式(3-20) 处「导致对地下压力 F_D 减小」这句按订正后的推导改为「增大」(或删去这句矛盾表述),消除同页自相矛盾。
  • 若「下坡增大」本质来自前进阻力/重心前移而非 G·cosθ 项,请把这条机理写进 F_D 的受力平衡式里、并在图3-11 受力图上标出对应的力,不要只用一句话定性带过。
证据 / 原文锚点
  • 物理 p76 式(3-16)「F_D = G·cosθ + F₃·cosα」(下压力按 cosθ)。
  • 物理 p77 式(3-20) 上文「下坡作业时,垂直地面方向上重力分力 G·cosθ 随 |θ| 增大而减小,导致对地下压力 F_D 减小」(结论:减小)。
  • 物理 p77 式(3-21) 下文「同时导致播深调节装置实际对地下压力 F_D 的增大,使实际播深变深」(同页结论:增大,自相矛盾)。
  • 物理 p80「下坡时,随着地表坡度 |θ| 的增大,下压力分布总体呈增大趋势,在 θ = -20° 时达到最大下压力均值 921.6 N……在 θ = 25° 时达到最小下压力均值 640.1 N」(图3-14 数据与正文新结论)。
建议修改 E2-003

§3.3.3 节研究「振动对播深」,但设计段落记录对象/窗口疑似沿用下压力段

📍 位置物理 p80–81 / 印刷 p66–67(§3.3.3/§3.3.4 设计段)
问题
本章振动/下压力相关节的设计段落与统计窗口仍存在与本节评价对象不一致的表述(记录对象、稳定窗口在相邻两节口径不同),读者无法判断该节究竟统计了播深还是下压力。
为什么是问题
- 现象:相邻两节稳定窗口 2~4 s(5 s 内)与 1~3 s(4 s 内)不一致(见 E1-009),且评价对象在「下压力」与「播深」之间表述不齐。 - 后果:读者无法判断该节统计量;更像从下压力段落复制后未改全,构成可溯源疑点。
修改建议
  • 核对各节「记录对象」与「统计时间窗」与本节评价指标一致;窗口不一致处给出说明。
证据 / 原文锚点
  • 物理 p80–81 §3.3.4「采集 RecurDyn 模拟试验过程中 4 s 内的 1~3 s 稳定范围内的下压力数据」与振动节窗口口径不一。
  • **另见**:段落复制残留视角见 C_language.md 对应条目。
建议修改 E2-004

式3-28/3-29 输出的 H、V 为小数(0.9060、0.0431),与式3-30 的「%」约束口径不一致

📍 位置物理 p86 / 印刷 p72(式3-28/3-29);物理 p87 / 印刷 p73(式3-30 约束)
问题
回归模型式3-28(合格率 H)常数项 0.9060、式3-29(变异系数 V)常数项 0.0431,量级表明 H、V 以「小数/分数」形式输出(0.9060≈90.60%、0.0431≈4.31%);但紧随的优化式3-30 约束写作「90%≤H≤100%、0%≤V≤5%」,H、V 用百分数;同一量三处口径不统一。
为什么是问题
- 现象:同一量 H 在式3-28 输出 0.9060(小数)、在式3-30 约束 90%(百分数)。 - 数据:若 H 单位为 %,式3-28 常数项应为 90.60 而非 0.9060;现 0.9060 只能解释为分数。 - 后果:直接代入式3-28/3-29 求得的 H、V 是小数,与式3-30 的 90%/5% 约束不在同一量纲,求解时若不换算会得到错误结论。
修改建议
  • 统一 H、V 量纲:要么回归模型也用百分数(常数项写 90.60、4.31),要么式3-30 约束改为小数(0.90≤H≤1.00、0≤V≤0.05),全章一致。
  • 在式中明确 H、V 是百分数还是分数,避免代入歧义。
证据 / 原文锚点
  • 物理 p86 式3-28「H = 0.9060 − 0.0372θ + …」、式3-29「V = 0.0431 + 0.0151θ − …」;物理 p87 式3-30 约束「90%≤H≤100%、0%≤V≤5%」。
建议修改 E2-005

式(4-8) 稳态误差用 m_b·g,但 m_b=99 kg 与表4-3 实际所用 990 N(≈101 kg)不一致

📍 位置物理 p98 / 印刷 p84(式4-8 + 表4-3);牵涉 物理 p94 / 印刷 p80(G=990 N)、物理 p96 / 印刷 p82(m_b=99 kg)
问题
稳态误差公式 e_ss = F(w)−(F_t+m_b·g) 里用的是播种单体质量 m_b。§4.2.1 把 m_b 定为 99 kg(重力约 970 N),但表4-3 那 11 行 e_ss 实际是用 990 N 算出来的(990 N 正是 §4.1.4 给的播种单体重力 G,对应质量约 101 kg)。论文给了两个不一样的单体重量(99 kg 与 990 N=101 kg),而公式符号 m_b 与表里实际用的数 990 N 不是同一个。
为什么是问题
- 复算对比(独立成行): - 用 990 N:θ=−25° 1290.2−(300+990)=0.2 ✓;θ=25° 2213.9−(1223+990)=0.9 ✓(11 行命中表值) - 用 970.2 N(m_b=99 kg):θ=−25° 得 20.0、θ=25° 得 20.7(11 行全部偏离表值约 +20 N) - 单体重量两处取值:G=990 N(物理 p94,≈101 kg) vs m_b=99 kg(物理 p96,≈970 N) - 后果:① 严格按式4-8 与 m_b=99 kg 算,e_ss 应全在 19~24 N,而非表中小值,公式与数据自相矛盾;② 同一单体重量在 §4.1.4 与 §4.2.1 差约 20 N;③ 读者按公式无法复现表4-3,影响可溯源。
修改建议
  • 统一播种单体质量/重力:确认是 99 kg 还是 101 kg(990 N),全章一处取值。
  • 修订式4-8:若表4-3 用的是 990 N(=G),公式应写 e_ss=F(w)−(F_t+G) 并把 G 定义为单体重力 990 N;不要公式用 m_b·g 而数据用 G。
  • 在式4-8 符号说明补 m_b(或 G)取值,使「公式符号—参数取值—表4-3 数据」三者闭合可复算。
证据 / 原文锚点
  • 物理 p98 式4-8「e_ss=F(w)−(F_t+m_b·g)」、表4-3 e_ss 列 11 行均可由 F(w)−(F_t+990) 精确复算;物理 p94「播种单体重量 G 为 990 N」;物理 p96「坡地播种单体质量 mb = 99 kg」。
建议修改 E2-006

正文「目标下压力 300~1223 N」与表4-3 列口径含糊(F_t 补充 vs F(w) 总)

📍 位置物理 p98 / 印刷 p84(§4.2.2(2) 表4-3「目标下压力 F_t/N」列)
问题
正文说「在地表坡度对应的目标下压力 300~1223 N 范围内,液压缸的调整时间在 0.16 s 至 0.65 s 内」。表4-3 里 F_t 列确是 300~1223 N,但同表「稳态值 F(w)」是 1290~2213 N。叙述「目标下压力」时只引 300~1223,容易让读者把系统实际工作的下压力误读为只有几百牛,而实际稳态下压力都在 1290 N 以上(含单体自重 990 N)。
为什么是问题
- 现象:F_t(目标下压力,即液压缸补充量)300~1223 N;F(w)(稳态实际下压力)1290.2~2213.9 N,二者相差约一个单体自重 990 N。 - 后果:正文未点明 F_t 是「液压缸输出/补充下压力」而非「单体总下压力」,读者易混淆 F_t(补充)与 F(w)(总)。 - 注:调整时间 0.16~0.65 s 与表4-3 t_g 列首末值一致,数值本身无误,仅口径表述需明确。
修改建议
  • 在正文/表注点明 F_t 为「液压缸目标输出(补充)下压力」、F(w) 为「含单体自重的实际总下压力」,二者差约 G=990 N。
  • 与 §4.3 节「输出下压力 300/800/1300 N」术语统一,避免口径漂移。
证据 / 原文锚点
  • 物理 p98「目标下压力 300~1223 N 范围内……调整时间在 0.16 s 至 0.65 s 内」;表4-3 F_t 列 300~1223、F(w) 列 1290.2~2213.9。
建议修改 E2-007

「50±0.88 mm」「50±3.22 mm」仍以对称 ± 概括非对称实测区间(PARTIAL)

📍 位置物理 p111 / 印刷 p97(§5.2.2(3) 图5-7);物理 p113–116 / 印刷 p99–102(§5.4、第六章 结论)
问题
本轮已把分项实测区间写清——上坡播深均值 49.36~50.34 mm、下坡 50.06~50.88 mm(物理 p111),这是进步;但仍把控制系统概括成「50±0.88 mm」、传统弹簧概括成「50±3.22 mm」。「50±0.88」=[49.12, 50.88],下界 49.12 与实测下界 49.36 仍对不上——这个 ±0.88 只取了上界一侧最大偏差(50.88−50=0.88)对称套用,下坡侧偏差被抹平。属「实测区间已给、对称 ± 概括口径未改」的半完成状态。
为什么是问题
- 区间对比(独立成行): - 控制系统分项实测:上坡 49.36~50.34 mm、下坡 50.06~50.88 mm,合并 [49.36, 50.88] - 概括写法:50±0.88 mm → [49.12, 50.88];下界 49.12 ≠ 实测 49.36 - 后果:对称 ±N 概括非对称区间,会让读者以为偏差对称分布,实际下偏更小;±3.22/±0.88 给了精确感,但其来源(单侧最大偏差)未说明。
修改建议
  • 改用真实区间表述:控制系统「播深均值 49.36~50.88 mm」、弹簧「46.78~53.22 mm」(或分上下坡列);若坚持 ±,注明「± 为偏离目标 50 mm 的最大单侧偏差」。
证据 / 原文锚点
  • 物理 p111「上坡……播深均值在 49.36~50.34 mm……下坡……播深均值在 50.06~50.88 mm……播深均值范围 50 ± 3.22 mm 降低至 50 ± 0.88 mm」。
  • **另见**:摘要/结论沿用视角见 D 类对应条目(互链,不重复合并)。
提醒 E2-008

限深轮橡胶剪切模量 3448 Pa 量级存疑,与 65Mn 的 7.9×10¹⁰ Pa 形成异常反差

📍 位置物理 p80 / 印刷 p66(§3.3.4(1) 坡地播种机 RecurDyn 模型材料参数)
问题
限深轮材料设为橡胶,剪切模量仍写「3448 Pa」(≈3.4 kPa);同处圆盘开沟器 65Mn 剪切模量为 7.9×10¹⁰ Pa。橡胶剪切模量通常在 10⁵~10⁶ Pa(0.3~2 MPa)量级,3448 Pa 偏低约两到三个数量级。
为什么是问题
- 现象:65Mn 用规范科学计数法 7.9×10¹⁰ Pa;橡胶用裸数 3448 Pa,写法与量级反差明显。 - 后果:若 3448 Pa 系漏单位/漏量级(应为 kPa 或 ×10⁶ Pa),会影响 RecurDyn 接触刚度与下压力仿真结果可信度;即便取自所引文献,也应核对原值与量纲。
修改建议
  • 核对橡胶剪切模量原始取值与单位,更正后统一计数法写法;与 65Mn 参数统一用科学计数法或统一用 Pa 裸数。
证据 / 原文锚点
  • 物理 p80「圆盘开沟器……剪切模量为 7.9×10¹⁰ Pa……限深轮的材料性能设定为橡胶……剪切模量为 3448 Pa,密度为 960 kg·m⁻³」。
提醒 E2-009

图3-12(b) 变异系数趋势描述与图不符(需据图复核)

📍 位置物理 p82 附近 / 印刷 p68 附近(§3.3.5 图3-12 变异系数面板)
问题
振动频率/坡度对播深变异系数的趋势描述(如「逐渐减小」)与图3-12(b) 形态可能不完全一致(图中存在先升后降段),属正文叙述与图查无对应。
为什么是问题
- 现象:图中变异系数存在非单调段,正文若概括为单调「逐渐减小」会漏掉上升段。 - 后果:趋势描述与所引图不一致,虽不改变结论方向,仍是可核查的图文偏差。
修改建议
  • 据图重述变异系数随频率/坡度的真实走向(先升后降处明写拐点)。
证据 / 原文锚点
  • §3.3.5 图3-12(b) 变异系数面板(物理 p82 附近,渲染核读)。
提醒 E2-010

图5-8 上/下坡播深均值「27.29~72.64 mm」「29.63~77.36 mm」把三个目标播深合并成一个区间,掩盖分档精度

📍 位置物理 p112–113 / 印刷 p98–99(§5.3.2 图5-8)
问题
控制策略试验目标播深为 30/50/70 mm 三档,但正文仍把播深均值笼统写成「27.29~72.64 mm(上坡)」「29.63~77.36 mm(下坡)」一个大区间。这区间跨了三个目标值(27.29 是 30 mm 档最浅、72.64 是 70 mm 档最深),混在一起后读者看不出每个目标播深各自的控制偏差。
为什么是问题
- 现象:27.29~72.64 mm 横跨 30/50/70 三个目标;合格率/变异系数同样合并报(88.36~96.44%、2.63~4.99%)。 - 后果:分目标精度被合并区间掩盖;30 mm 档(浅播)的控制偏差无法从合并区间判断;属报告口径笼统。
修改建议
  • 按目标播深 30/50/70 mm 分别给均值±偏差与合格率,或在图5-8 已分档基础上正文也分档点出关键档(尤其 30 mm 浅播)的控制偏差。
证据 / 原文锚点
  • 物理 p112–113「不同目标播深下的播深均值、播深合格率和播深变异系数分别在 27.29~72.64 mm、88.36%~96.44% 和 2.63%~4.99% 范围内……(上坡);……29.63~77.36 mm、87.67%~96.31% 和 3.43%~4.96%……(下坡)」。