打印面上的层纹和你要研究的凸起一样大：仿生叶片的公差敏感性与风洞阻塞修正后的真实 Cp

推荐优先级：高 · 族 B（自建低成本装置）／参赛子类：Energy — Sustainable Materials and Design · 资源需求：总器材 ¥1,500–2,500，最长采购周期为自制风洞整流段约 2 周 · 技能取向：机械制造 / 流体测量 / 误差传播

1 · 研究问题

当 3D 打印叶片的表面层纹幅值与被研究的仿生凸起处在同一量级时，风洞里测到的性能差异有多少来自设计、多少来自打印工艺？以及：一个小型风洞里测出的功率系数 Cp = 0.48，做完阻塞修正之后还剩多少？

2 · 研究背景与空白

背景。 垂直轴风力机（Vertical Axis Wind Turbine, VAWT）的性能用功率系数（power coefficient）Cp = P / (½ρAV³) 衡量，其中 A 是叶轮扫掠面积、V 是来流风速，理想上限是贝兹极限（Betz limit）0.593。仿生凸起（tubercle）源自鲸类鳍前缘的波状结构，被认为能推迟失速、改善大攻角下的升阻比（lift-to-drag ratio, L/D）。但有两件与"叶片形状好不好"完全无关的事会直接进入这个数字。其一，熔融沉积（Fused Deposition Modeling, FDM）打印件的表面由逐层堆叠形成，层纹（layer lines）的峰谷幅值与层高同量级（常用 0.1–0.3 mm），而凸起、凹坑这类仿生特征在一片小尺寸模型叶片上的特征高度往往也只有零点几毫米——制造噪声与被研究的信号落在同一尺度上。其二，小型风洞里叶轮占据了试验段可观比例的截面，被叶轮阻挡的气流绕流加速，使叶轮实际感受到的来流快于试验段测得的风速；不做阻塞（blockage）修正，而 Cp 与 V 的三次方成反比，误差会被放大三倍地灌进结果。

已有工作到哪一步。 两个直接相关的获奖项目构成一条连续的脉络。EGSD032（2025, Grand Award, "Hammerhead-Inspired Wind Energy Generation"） 把双髻鲨头翼的凸起沿弦长布置到 VAWT 叶片上，报告 "The CFD simulations produced a 67% improved lift-to-drag ratio for the tubercle inspired blade, while physical experiments in real world operations showcased a 12% RPM gain and a 20% improvement in generated voltage"。EGSD028T（2026, Grand Award + Special Award, "WinDTurbinePower: A Novel Dimple-Tubercle Blade"） 在 NACA0012 上叠加凸起与高尔夫球凹坑，报告 "a power coefficient (Cp) of 0.48 for the hybrid, 60% higher than 0.30 for the smooth blade"、"the hybrid achieved a Cp of about 80% of the Betz limit (0.593)"。两篇都用 OnShape 建模、SimScale 仿真、3D 打印实测，两篇都没有报告打印层高、打印方向或表面粗糙度，也都没有提到风洞阻塞修正。 另可查 2024 EGSD016（改良叶片的螺旋风机，Second Award）。更大的背景是：本赛道两届 35 篇获奖摘要中，calibrat 出现 0 次、uncertaint 出现 0 次、reproducib 出现 0 次；作为横向参照，同两届四个工程赛道全部 173 个获奖项目里 error bar 与 sensitivity analys 合计出现 0 次。检索还显示 wind turbine blade 五届累计 6 项获奖，"仿生凸起"这一支连续两年获奖（EGSD032 → EGSD028T），再叠加一种仿生特征几乎必然撞车；但测量方法侧（阻塞修正、打印公差）完全空白。

空白在于：没有人量化过"制造公差散布"与"仿生特征效应"这两个量的相对大小，也没有人对这一支项目普遍使用的小型风洞报告过阻塞比与修正结果。这个空白适合一年期学生课题，因为它不需要任何新理论——阻塞修正是标准做法，层高与打印方向是切片软件里的两个下拉框——技术门槛全在"肯不肯把这两件事当自变量做"；而且结论正负都成立：若工艺散布远小于设计效应，那是对前两年结论的一次独立支持，同样值得报告。

3 · 可检验假设

H1 几何完全相同、仅打印工艺不同（3 种层高 × 2 种打印方向）的一组"光滑基准叶片"，其 Cp 的组间极差 ≥ 0.05，即达到 EGSD028T 报告的 0.30 → 0.48 提升幅度（0.18）的四分之一以上。

数值依据的推算过程：对弦长约 50 mm 的模型叶片，层高 0.1 mm 与 0.3 mm 对应的相对粗糙度 k/c 约为 0.002 与 0.006，相差三倍；在模型尺度的低雷诺数下，这个尺度与层流边界层厚度同量级，足以改变转捩（transition）位置。转捩位置移动对 L/D 的影响量级需核实（须查低雷诺数翼型粗糙度效应的公开风洞数据，并核实其雷诺数区间是否覆盖本装置）。保守取 L/D 变化 10%，并近似认为该区间内 Cp 与 L/D 单调同步，则 Cp 变化约 0.048 ≈ 0.05。

H2（机制假设，H1 被否证时仍有内容）若工艺散布显著小于 0.05，则 Cp = 0.48 这一偏高读数的主因是风洞阻塞未修正：本装置实测阻塞比若为 ε，一阶修正会把等效来流速度抬高约 τ 倍，由 Cp ∝ V⁻³ 得 Cp 下降约 3τ。修正后 Cp 相对下降 ≥ 20%。 依据：若阻塞修正给出 τ ≈ 7%，Cp 下降约 21%，0.48 → 0.38。具体采用 Maskell 还是 Glauert 修正、以及它们对 VAWT 这类介于致动盘与钝体之间的对象是否适用，需核实——这本身是本课题必须交代清楚的一个方法学选择，不能默认。

4 · 量化验收标准

1. 方法学校验（硬门槛，不过关则后续全部结论无效）。 两件事都要过。（a）风洞标定：在同一试验段内测一个已知压力系数分布的标准体（例如标准圆柱或标准翼型剖面）的表面压力系数分布，与公开参考值的偏差 ≤ 10%。该参考值需核实——须核实所选标准体的公开 Cp 分布数据出处、其雷诺数区间是否覆盖本装置的工作点、以及原始文献是否报告了不确定度。注意术语撞车：此处的 Cp 是压力系数，与主结论里的功率系数 Cp 是两个不同的量，全文必须用不同符号并在图注中说明。（b）层纹幅值实测：用轮廓仪（校内是否有需核实）或标定过的 USB 显微照片实测每种工艺组合的层纹峰谷幅值与 Ra，并证明在"精修基准组"（最小层高 + 最优方向 + 统一后处理）上，层纹幅值 ≤ 凸起特征高度的 1/10。这一组是整个实验的零点：没有一组低噪声试件，"设计效应"就没有可比的参照系。 这两条任一不过关，后续所有叶片间比较都无意义。
2. 统计口径预先写死，且必须区分两个 n。 每种工艺组合打 3 件（不同试件，n_试件 = 3），每件在同一工况下重测 5 次（同一试件重测，n_重测 = 5）。这两个 n 不是一回事：同一试件重测只反映装夹与气流的短期变异，试件间差异才是本课题要测的制造公差。 全部方差分析必须按嵌套结构做，报告"试件间方差 / 试件内方差"的分解比例；任何把 3×5 = 15 当作 n = 15 的做法都构成伪重复（pseudoreplication），直接判本条不达标。不确定度按 A 类（重复测量标准差）与 B 类（风速仪、功率测量、空气密度、扫掠面积的仪器与读数分量）分别列出后合成，并写明 Cp 的相对不确定度按 3×u_r(V) 与 u_r(P) 的平方和合成——风速的分量先乘 3 再进平方和，这一步漏掉会把主要误差源低估一个数量级。
3. 修正前后必须同口径逐点对比。 报告实测阻塞比 ε（实测叶轮扫掠面积 ÷ 实测试验段截面积），并对每一个尖速比（tip speed ratio, TSR）工况给出"未修正 Cp"与"修正后 Cp"两个值并列成表，附修正式与其适用条件的说明。只报修正后数字、不给修正前，等于让读者无法判断修正的量级。
4. 参数点数量。 光滑基准组 3 层高 × 2 打印方向 × 3 件 = 18 件；仿生特征组至少覆盖 2 种层高 × 3 件 = 6 件（用于回答"设计效应是否穿透了工艺噪声"）。每件扫 ≥ 6 个 TSR 点。叶轮平面风速剖面在 ≥ 9 个位置测量，报告不均匀度；若不均匀度 > 5%，全部 Cp 必须使用面平均风速而非单点风速，并在方法中写明。
5. 阴性结果预案（开题时即写好）。 若工艺散布显著小于仿生效应，主结论改为"在本装置分辨能力下，制造公差不足以解释已报告的提升，且阻塞修正把 Cp 从 0.48 拉低到 X"——这仍是完整结论，但必须给出误差棒证明有能力分辨 0.05 这个量级的差异，否则等于没做。
6. 可复现性。 CAD/STL 源文件、切片参数（层高、方向、填充率、打印速度、挤出温度）、风洞几何图纸、采集脚本、分析脚本、全部原始时序数据开源。

5 · 数据与工具

用途	来源 / 工具
叶片与整流段制造	校内 FDM 打印机。必须记录并作为自变量报告的参数：层高（0.1 / 0.2 / 0.3 mm）、打印方向（2 种）、填充率、喷嘴温度。机型与耗材批次需核实并记录
层纹幅值与表面粗糙度	轮廓仪（校内是否具备需核实；若无，用 USB 显微镜 + 标定尺按显微照片测峰谷幅值，源文件已列此替代方案）。测量方法的重复性必须先自测并报告
风速	风速仪，¥300–800，现货。须核实其响应时间与在叶轮尾流湍流场中的适用性；叶轮平面多点平均，不得单点代表
输出功率 / 转矩	电磁制动 + 电流法，或直接测发电功率，¥200，现货
风洞	校内若无，用变频轴流风机 + 蜂窝整流段自制，¥600–1,200，约 2 周。试验段截面尺寸必须实测并用于阻塞比计算
阻塞修正方法	Maskell 或 Glauert 一类标准修正。具体选用哪一支、其对 VAWT 的适用性、以及原始出处，全部需核实——本文件不指定标准号
对照基准（仅用于校验与对比，不计入本项目的数据贡献）	EGSD032 的 67% L/D 提升 / 12% RPM 增益 / 20% 电压提升；EGSD028T 的 Cp = 0.48 vs 0.30、Betz limit 0.593

须核实而非引用的量： ① 低雷诺数下粗糙度对 L/D 影响的量级（H1 的推算依赖它）；② EGSD028T 与 EGSD032 所用凸起与凹坑的实际特征高度——两篇摘要都没有给出，若拿不到，本课题的"同量级"论断只能基于自己复刻几何的实测尺寸，须在局限性中明确声明；③ 所选风洞标定标准体的公开 Cp 分布参考值与其雷诺数适用区间；④ 阻塞修正式的出处与适用条件。以上四项在方案里一律标注为待核实，不得引用任何未经查证的数字。

6 · 方法路径

1. 搭风洞 + 用标准体标定 + 完成 §4 第 1 条两项校验。 建整流段、实测试验段截面、跑标准体压力系数分布、实测层纹幅值并确认精修基准组达到 1/10 判据。这一步不过不许继续。
2. 核实器材能力边界。 逐条验证以下静默失败点（不报错但给错结果）：（a）打印层纹尺度与被研究的表面特征同量级——这是本课题的立论点，必须先量化再谈其他；（b）风洞阻塞比未修正——修正前后 Cp 差多少，必须先用一个已知工况算出来；（c）风速仪放在叶轮尾流或壁面边界层内，读数偏低而仪器一切正常，Cp 因此被三次方地抬高；（d）光电测速的倍频与混叠——叶轮上贴了多个反射标记时转速会成倍读高，采样率低于转频两倍时会读出一个稳定但完全错误的低频；（e）空气密度未按当日温压修正，秋冬室温波动 10 °C 对应 ρ 变化约 3%，直接线性进入 Cp。
3. 打满 18 件光滑基准 + 6 件仿生特征试件，逐件测层纹幅值与 Ra，建立"工艺参数 → 表面形貌"的对照表。此表独立于气动结果，即使后续风洞出问题也是一个可交付产物。
4. 跑满光滑基准组的风洞测试（18 件 × 6 个 TSR × 5 次重测），做嵌套方差分解，得到工艺造成的 Cp 组间极差——这是 H1 的直接判据。
5. 跑仿生特征组，在同一工艺条件下与光滑基准比较，判断设计效应是否穿透工艺噪声；并对每个工况同时报告未修正与阻塞修正后的 Cp。
6. 独立交叉校验。 用另一种原理测同一个量：用自由旋转减速法（spin-down，从断电后的转速衰减曲线反推轴系机械功率与摩擦损失）独立估计机械输出，与电功率法所得对照；风速端另用皮托管与风速仪两种原理交叉验证同一测点。两条链若系统性不一致，说明测量链而非叶片有问题。

7 · 新颖性边界

本课题不声称： - 不声称首次把仿生凸起用于风力机叶片。EGSD032（2025, Grand Award）与 EGSD028T（2026, Grand Award + Special Award）已连续两年获奖，且这一支五届累计 6 项。 - 不声称提出新的仿生几何。再叠加一种仿生特征必然撞车，本课题明确不走这条路。 - 不声称 EGSD028T 的 Cp = 0.48 是错的。本课题只声称它未报告阻塞修正与打印公差——这是摘要原文可核对的事实，不是对结论的否定。 - 不声称提出新的阻塞修正方法。Maskell / Glauert 一类修正是既有标准做法，本课题只是把它补上。 - 不声称本装置的 Cp 绝对值可与工业风机比较。模型尺度雷诺数远低于全尺寸，绝对值不外推。

已有工作做到哪： 设计侧，凸起 + 凹坑的复合几何已用 OnShape 建模、SimScale 仿真、3D 打印实测，报出 67% 的 L/D 提升与 0.48 的 Cp；测量侧，两篇均未报告层高、打印方向、表面粗糙度与阻塞比，本赛道两届 35 篇获奖摘要中 calibrat 与 uncertaint 各出现 0 次。

本项目的贡献（属「方法可复现性贡献」型）： 已有工作评估的是换一个几何能提升多少；本项目评估的是换一次打印工艺会波动多少，以及未修正的阻塞把这个数字抬高了多少——即给出一条"制造公差 → Cp 散布"的传递关系和一份阻塞修正前后的对照表。这是主结论，不是附加内容：交付物就是这两张表本身，而不是"我的叶片更好"。

为什么有价值： 如果工艺散布可比于仿生特征的提升，那么这一支连续两年获奖的结论就需要重新定界——不是被推翻，而是被加上一个此前缺失的置信区间；如果散布远小于提升，那是对前作的一次独立支持，同样有价值。而阻塞修正是一个任何后来者都能直接复用的固定动作。

若结论不显著： 若工艺组间极差的置信区间跨过 0.05，即"在本装置分辨能力下无法判定制造公差是否可比于设计效应"，这同样是有效结论，但必须给出误差棒证明有能力分辨 0.05 这个量级的差异，否则该结论没有信息量。

8 · 决策门槛（go / no-go）

🟢 2026-08-25（第 3 周末）：风洞建成，标准体标定通过，层纹幅值实测完成。 若标准体压力系数分布偏差 > 10%，说明整流段品质或压力测量链有问题。降级路径：放弃 Cp 绝对值的声称，把主结论框架从"修正后的真实 Cp 是多少"收缩为"同一装置上工艺组间的相对散布有多大"——相对比较不要求绝对标定，主结论（制造公差 vs 设计效应）完整保留。这是本课题最重要的一条保险。

🟡 2026-09-14（第 6 周末）：18 件光滑基准试件打印与形貌测量全部完成。 若精修基准组的层纹幅值做不到凸起高度的 1/10，说明校内打印机的分辨能力不足以构造零点。降级路径：改用统一后处理（打磨 + 涂层）制造低噪声基准组，并把"后处理"作为第三个工艺自变量纳入——主结论框架不变，反而多了一个可操作的工程建议。

🟢 2026-10-12（第 10 周末）：光滑基准组风洞测试跑满，H1 判据出结果。 此后才允许把时间投入仿生特征组。

2026-12-15：数据采集截止，留 3 周分析与展板。

⚠️ 与「附属赛后不得修改」的冲突点： ① 阻塞修正式的选择（Maskell / Glauert / 不修正三选一）；② 层高档位（0.1 / 0.2 / 0.3 mm）与打印方向的定义；③ Cp 计算中风速取单点还是面平均。这三项一旦在附属赛后想改就不被允许。必须在 2026 年 11 月底前全部定死，并在研究计划中写明选择理由与一次敏感性检验结果。 事后改修正式等价于事后改结论量级，是评委最敏感的问题。

须提前核实而非边做边发现的事项： - 轮廓仪的校内可及性。若没有且显微照片法的重复性不达标，§4 第 1 条的（b）过不了，整个课题的前提塌掉。这件事在第 1 周就要问清楚，不要拖到打完 18 件试件才发现没法测。 - EGSD032 与 EGSD028T 的凸起/凹坑实际尺寸。两篇摘要都没有给，建议 2026 年 8 月底前设法取得全文或向作者询问；取不到就用自己复刻几何的实测尺寸，并在局限性中声明。 - 阻塞修正式对 VAWT 的适用性。这是一个真实的方法学争议点，不是一句"按标准做"能带过的，必须自己查证并写明依据。 - 本文件的证据基础是本地五届语料、未做外部文献核查，"没人做过"不等于不存在前作。开题前必须自己补一轮独立文献检索（重点查低雷诺数翼型粗糙度效应与小型风洞阻塞修正两条线）。

已知困难及其定位： Cp 正比于风速的三次方，风速测量 3% 的误差就是 Cp 的 9%。这个困难本身就是研究内容——本课题问的正是"测到的差异里有多少不是叶片给的"，所以误差传播不是要绕开的障碍，而是要量化的对象。

安全与合规： 高速旋转叶轮必须加防护网罩，工作电压 < 12 V，需常规 Form 3（风险评估），列明护目镜、防护罩与成人监督。不涉及化学品、激光、辐射，无 SRC 事前审批需求。

选择前提： 适合愿意做重复制造与重复测量的学生。最耗时的一环是打满 24 件试件并逐件测形貌，它既不出成果也不好看，但没有它就没有试件间方差，整个主结论不成立。 不适合希望"做一个更好的叶片"的学生——本课题的定位是标定，不是超越。

预算裁剪顺序： 先砍自制风洞（¥600–1,200），改用校内现有风机 + 自制整流段，代价是阻塞比更大、修正更重要；再砍仿生特征组（6 件），只做光滑基准的工艺散布；最后砍打印方向这一维，只做 3 种层高 × 3 件。砍到最后主结论"制造公差造成的 Cp 散布有多大"仍然成立。

🔴 展台能否展示实物：可以。 叶片与打印试件是固体塑料件，非液体、非粉末、非化学品、非电池组，可以带上展位；风洞整体过大不能带，用照片与视频替代。这是本文件中少数能带实物的条目之一，但注意展板的说服力仍主要由方差分解图与修正前后对照表承担，不要把版面让给叶片照片。

评分： O=8, W=9, F=7, S=8, Fit=9 → base=81.8，h=9 → [73, 91]，置信度：中

---