🌿 A题代码全解 · 药材烘干四问

2026-09-22 08:59 生成 · 左侧目录选择文件 · 快捷键 H 隐藏/恢复注释

这是什么

2025 国赛 A 题「圆柱形药材烘干」的全部支撑代码浏览器：左侧是从工作区逐字节读取的原版代码，右侧是面向初学者的逐行详细讲解（含库函数说明与函数间跳转）。两侧同屏同步滚动：讲解块所在的行在左侧以空白占位，点击顶部「隐藏注释」后空白自动消失、右侧收起，只剩纯原版代码。

怎么用

· 讲解里的蓝色虚线链接（如 solver.py 的 solve）可点击，直接跳到对应函数定义行并高亮闪一下；
· 顶部「☰ 函数大纲」列出当前文件全部函数/类，点击任一个跳到定义处；
· 快捷键 H 隐藏/恢复注释；目录树顶部输入框可按文件名或用途筛选；地址栏会记住当前文件，刷新不丢。

建议阅读顺序

顺序	文件	对应论文
① 背景参数	`A/PROBLEM_FACTS.json`	题目附录数据
② 数据入口	`A/code/data_check.py`	第三章数据与假设
③ 核心求解器 ★	`A/code/solver.py`	第四章建模（全文心脏）
④ 四问入口	`A/code/problem1~4.py`、`main.py`	第四~七章
⑤ 结果输出	`A/code/results_io.py`、`portable_excel.py`	表1~6
⑥ 验证与敏感性	`verify_cylinder.py`、`convergence.py`、`endface_check.py`、`sensitivity_analysis.py` 等	第八章
⑦ 相平衡扩展	`code/phase_model.py` → `long_phase.py` → 其余	第九章

一句话模型

把药材看作半径 2 cm、长 25 cm 的圆柱，温度场 T(r,t) 与含水率场 C(r,t) 沿径向一维分布，用柱坐标有限体积法离散、BDF 刚性积分器推进时间；表面是第三类边界（对流换热+对流传质），烘干终点（最大含水率降到 0.15 kg/kg）由事件函数自动定位；第四问通过材料坐标 ξ=r/R(t) 处理实测半径收缩。核心结果：问题三烘干 57.47 h，问题四 51.09 h。

必须会讲的数字

水分收支守恒误差 < 1e-7；Bessel 解析解验证误差 ~3e-7；二维端面效应 < 9 s（一维简化成立）；附录3物性+收缩 = 25.2 h，附录4物性+固定半径 = 129.8 h（“收缩加快、附录4物性减慢”）。