简介：螺旋桨流体仿真.zip 面向机械工程、航空航天与船舶推进方向的学习者和设计人员，围绕螺旋桨在水或空气中的推力、扭矩与阻力分析，提供一套可复现的 SolidWorks Flow Simulation 仿真资料。压缩包共105个文件、约92.06MB，包含sldprt与sldasm三维模型及装配文件、fld和par_val等仿真设置与结果数据、cv_mm等网格相关文件，以及txt说明、html报告和log、stdout运行记录，覆盖从几何模型到计算结果的完整链路。已有772人学习，说明其在实际项目与课程练习中具有一定参考价值。读者可据此对照螺旋桨叶片形状、角度和排列对性能的影响，查看流场与压力、速度结果，理解旋转边界及入口出口等边界条件设置，并借助日志和报告排查仿真过程。适合具备SolidWorks基础、希望系统练习CFD或优化螺旋桨设计的中级用户。

1. 一副桨叶卡住的三个星期

某同学把桨叶模型画得很漂亮，导入流体模块的第一天就卡住了：网格生成直接报错，改完网格算出来的推力是负的，再调边界条件又变成残差来回震荡。问题不在理论推导，而在流程——螺旋桨流体仿真本质上是几何、流体域、网格、旋转区域、后处理五件事相互关联，任何一环处理不当，结果都不可信。这套基于 SolidWorks 的螺旋桨流体仿真工程，把上述环节按可复现的顺序串了起来：参数化建模驱动桨叶形状，旋转域与静止域分开切，稳态用旋转参考系、瞬态换滑移网格，最后把桨叶表面的力和力矩换算成推力系数、扭矩系数与效率。它适合要交毕设课设的学生，也适合第一次把 CAD 与 CFD 接起来的从业者。

2. 需求与技术选型：一套螺旋桨仿真到底要交付什么

2.1 需求拆解：仿真真正要回答的是三个数加两张图

把"做个螺旋桨仿真"这句话拆开，交付物其实非常具体。三个数是推力 T、扭矩 Q、效率 η，两张图是桨叶表面压力分布和叶尖附近的涡量或速度云图。T 和 Q 直接由桨叶表面的压力积分与黏性力积分得到，η 则由前进比 J、推力系数 Kt、扭矩系数 Kq 换算：η = J·Kt/Kq。压力云图用来说明载荷分布是否合理，叶尖涡图用来说明网格是否捕捉到了真实的流动结构。

需求拆到这个粒度之后，技术选型的判断标准就清楚了：工具能不能稳定地给出这三个数，能不能在改动几何参数后快速复算，能不能导出足够细的表面数据做后处理。剩下的事情——求解器内核是有限体积还是格子玻尔兹曼、网格是贴体还是笛卡尔——都属于实现细节，不该成为选型的第一顺位。

2.2 为什么是 SolidWorks 加内嵌流体模块，而不是另起一套 CFD

做课设和毕设的典型约束是：时间只有一个学期，机器是一台普通工作站，人只有一个，而且这个人的强项通常在建模而不在数值方法。在这种约束下，SolidWorks 内嵌的流体仿真模块有几个实打实的优势。几何和仿真在同一个文件体系里，改一个螺距角，网格和边界条件不用重新导一遍；笛卡尔网格加部分单元的处理方式对复杂曲面的容错率高，桨叶这种扭曲面不容易生成失败；湍流模型里 k-ω SST 和 k-ε 都有，外流场和旋转机械够用。

代价也要说清楚。笛卡尔网格在薄壁和尖角处的分辨率不如贴体网格，叶梢这种曲率剧烈变化的位置必须手动加密，否则叶尖涡会被抹平。求解器对旋转机械的支持主要靠旋转区域，稳态算的是旋转参考系下的近似解，真正的时间效应要靠瞬态滑移网格，计算量会涨一个数量级。

提示：选型时先确认自己的交付物是"趋势正确"还是"数值精确"。课设通常前者就够，工业交付才需要后者。

2.3 能力边界对照：什么时候该换工具

把三条常见路线放在一起比，边界就清楚了。下表里的判断依据是网格类型、旋转机械处理方式、学习成本和典型精度，不涉及任何具体版本。

我一般会这样建议：如果目标是把 J-Kt-Kq 曲线画出来、和公开图谱对得上、误差控制在 10% 以内，SolidWorks 这条路线完全够；如果要研究空化、非定常尾涡演化或者做优化循环，尽早换通用求解器，别在工具上耗时间。

3. 总体设计：几何到推力曲线的五段式流程

3.1 五段式流程与模块职责划分

整套方案拆成五段，每段有明确的输入输出，段与段之间用文件而非口头约定衔接。这样做的直接好处是：出问题时能快速定位是哪一段的锅，而不是对着一个黑匣子反复试参数。

第一段几何建模，输出是参数化的桨叶实体和旋转轴，关键约束是直径、螺距、弦长分布、桨毂直径这几个驱动尺寸必须写成方程式，改动时全模型联动。第二段流体域构建，输出旋转域和静止域两个封闭实体，关键是两个域的交接面要贴合、不能有缝隙或重叠。第三段网格划分，输出是带局部加密的网格，关键是叶梢、前缘、桨毂根部三处的加密足够。第四段求解设置，输出是边界条件、旋转区域参数和目标函数，关键是入口速度、出口压力、转速三者与前进比 J 自洽。第五段后处理，输出是 T、Q、η 和云图，关键是从表面目标值到无量纲系数的换算不能错。

这五段里，最容易返工的是第二段和第三段。流体域尺寸给错，后面网格再怎么加密都是错的；网格没加密，后面收敛再好结果也是假的。

3.2 关键参数表：把所有输入钉死在一张表上

仿真结果不可复现，十有八九是因为参数散落在各个对话框里没记录。我的做法是先建一张参数表，所有输入以这张表为准，建模、流体域、求解设置都从表里取值。下面是一组船用螺旋桨的示例参数，换项目时整表替换即可。

雷诺数按 Re = n·D²/ν 估，这组参数下约 1.25×10⁶，属于湍流充分发展区，用 SST 是合适的。如果换成小直径低转速的空气桨，Re 会掉到 10⁴ 量级，这时候层流与转捩的处理比湍流模型选择更关键。

3.3 文件组织与版本管理：别让结果覆盖结果

一个螺旋桨项目跑下来，动辄十几个算例：不同 J、不同网格密度、稳态和瞬态各一套。文件管理混乱是隐形成本，改天想找回"那次收敛得特别好的算例"却找不到，只能重算。

我习惯按"模型 / 算例 / 结果"三层组织。模型层只放参数化实体和方程式，算例层按 J 值命名，结果层存目标值导出的 CSV 和云图截图。命名带上日期和网格级别，例如 J080_mesh3_steady，一看就知道是什么。每次改参数前先复制一份算例文件夹，不要在原算例上直接改，这样任何一次"改坏了"都有后悔药。

注意：结果文件不要和 CAD 文件混在同一目录，流体模块的结果文件夹体积增长很快，混放会让备份变得很痛苦。

4. 关键实现：从参数化建模到推力扭矩提取

4.1 参数化建模：用螺距角把几何和物理连起来

螺旋桨几何的核心不是画得像，而是让每个半径处的安装角符合螺距定义。等螺距条件下，半径 r 处的螺距角 β 满足 tanβ = P/(2πr)，半径越大角度越小，这就是桨叶看起来"越往外越平"的原因。把这个关系写成脚本，先算出各截面的角度和弦长，再回到 CAD 里用方程式或放样驱动。

import math D = 0.25 # 桨直径，m P = 0.20 # 螺距，m R = D / 2.0 r_hub = 0.02 # 桨毂半径，m n_sec = 6 # 截面数 for i in range(1, n_sec + 1): # 从桨毂到叶梢等距取截面 r = r_hub + (R - r_hub) * i / n_sec # 等螺距条件下的螺距角 beta = math.degrees(math.atan(P / (2 * math.pi * r))) # 弦长按椭圆分布衰减，0.35D 为根部弦长系数 chord = 0.35 * D * math.sqrt(max(0.0, 1 - (r / R) ** 2)) print(f"r={r:.4f} m beta={beta:.2f} deg chord={chord:.4f} m")

这段脚本输出的是建模用的截面表。beta 是小节 3.2 里螺距和直径的直接函数，改 P 或 D 后重跑一遍就能得到新表，把结果填进 CAD 的方程式或放样草图，几何就跟着变了。chord 那一行是经验分布，根部粗、梢部细，实际项目里如果有型值表就以型值表为准，脚本只做校核。值得一提的是，桨毂附近的 beta 会接近 80 度，这个位置的截面几乎垂直于旋转轴，建模时容易被切成自相交面，通常的做法是把第一个截面从 0.25R 开始取，根部用圆角过渡。

4.2 流体域构建：旋转域和静止域怎么切

流体域分两块。静止域是以桨盘为中心的圆柱或长方体，入口距桨盘面 3D 到 5D，出口 8D 到 10D，侧向边界 5D。这个比例不是玄学，入口太近会让来流被桨的诱导速度污染，出口太短会让尾流在出口面上还没充分发展，压力边界会把涡反射回来。旋转域是一个包住桨叶的小圆柱，直径约为 1.1D 到 1.2D，轴向长度略大于桨叶在该半径处的轴向投影。

D = 0.25 V = 4.0 # 来流速度 m/s n = 20.0 # 转速 rps J = V / (n * D) upstream = 4 * D # 入口到桨盘面 downstream = 10 * D # 桨盘面到出口 radial = 5 * D # 侧向边界 rot_dia = 1.15 * D # 旋转域直径 rot_len = 0.45 * D # 旋转域轴向长度 print(f"J={J:.3f}") print(f"静止域总长={upstream + downstream:.3f} m, 半径={radial:.3f} m") print(f"旋转域直径={rot_dia:.4f} m, 轴向长度={rot_len:.4f} m")

参数的含义很直白，upstream 和 downstream 是硬约束，别为了省网格压缩它们。rot_dia 取 1.15D 是为了让叶梢和旋转域壁面之间留出至少 0.05D 的间隙，间隙太小会让叶梢附近的流动被壁面人为加速，推力偏高；间隙太大则旋转域内的网格量浪费在无意义的区域。旋转域和静止域的交接面在建模阶段就要做成同一个圆柱面，导入流体模块后两域共面，不用再做布尔运算，这一点能省掉大量网格失败。

提示：旋转域轴向长度不要卡得太紧，桨叶有侧斜和纵倾时，轴向投影会比正投影长，留 0.45D 是给这类几何留余量。

4.3 网格划分：全局级别和局部加密的取舍

流体模块的网格是笛卡尔网格加部分单元，控制手段主要是初始网格级别、最小缝隙尺寸、最小壁面厚度和局部细化。初始网格级别一般取 3 到 5，级别每加一，网格量大约翻倍。最小缝隙尺寸决定它能分辨多窄的通道，最小壁面厚度决定壁面附近能铺几层，这两个值给小了网格量会爆，给大了叶梢直接被抹圆。

工程上的做法是：全局级别先取 4 跑一遍，看叶梢附近的网格是不是至少有三到四个单元跨越叶梢厚度，不够就在叶梢和前缘加局部细化区域，细化区域做成包住叶梢的小圆柱，半径给 0.1D 左右。桨毂根部也建议加一个细化区，那里曲率变化剧烈，不加密容易出现负体积或部分单元被切成薄片。

网格是否够用不靠感觉，靠网格无关性验证。同一个 J 下用粗、中、细三套网格各算一次，比较推力。三套网格的推力偏差如果落在 2% 以内，说明网格已经收敛，后面所有算例沿用中等那套就行，不用每算一次都加密。

表里的数值是示例量级，真实项目要以自己的算例为准。判断方法和趋势是一致的：粗网格明显偏低，是因为叶梢被抹平、吸力面的负压峰值捕捉不到；中网格和细网格差 1.4%，已经进入可接受区间。如果粗中细三套网格的结果一路往上飘没有收敛趋势，多半是流体域尺寸不够或者边界条件有问题，而不是网格还不够密。

4.4 旋转区域与边界条件：稳态与瞬态的两种走法

稳态算例用旋转区域的方式，在流体模块里指定旋转轴、旋转方向和角速度，角速度由转速换算：ω = 2πn/60，1200 rpm 对应 125.66 rad/s。稳态求解的是旋转参考系下的定常近似，桨叶相对流体静止，离心力和科氏力以源项形式加进去。这种方法算得快，对 J 不太小、流动分离不严重的工况足够准。

瞬态算例要打开时间相关选项，把旋转区域改成滑移网格模式，让旋转域真的转起来，时间步长按每步转 1 到 2 度给，1200 rpm 下每步约 1.4×10⁻⁴ s，转三到五圈看推力是否进入周期稳定。瞬态的计算量大概是稳态的二十倍以上，只在需要看非定常力波动或者叶尖涡脱落时才做。

边界条件方面，入口面给速度，大小等于来流速度 V，方向沿轴向；出口面给静压，表压 0；侧向边界给滑移壁面或对称条件，视流体域形状而定。桨叶表面和桨毂表面默认是无滑移壁面，粗糙度按加工情况给，课设里给光滑就行。

注意：入口速度、转速、直径三者要和 J 对上，J 算错了后面所有无量纲系数全错，而且错得很隐蔽——曲线形状看着正常，数值整体偏移。

4.5 后处理：从表面力到 Kt、Kq、η 的换算

求解前就要在桨叶表面设好目标：沿轴向的力（推力 T）和绕旋转轴的力矩（扭矩 Q）。求解过程中这两个值会被实时记录，收敛判断也靠它们的曲线。算完之后导出目标值历史，用下面这段脚本换算成无量纲系数，顺便做一次收敛性检查。

import csv rho = 998.2 # 介质密度 D = 0.25 # 桨直径 n = 20.0 # 转速 rps V = 4.0 # 来流速度 rows = [] with open("goals_history.csv", newline="") as f: reader = csv.DictReader(f) for row in reader: rows.append((float(row["iteration"]), float(row["force_x"]), # 轴向力，N float(row["torque_z"]))) # 绕轴力矩，N·m # 取最后 20% 迭代做平均，避开启动段的剧烈波动 tail = rows[int(len(rows) * 0.8):] T = sum(r[1] for r in tail) / len(tail) Q = sum(r[2] for r in tail) / len(tail) J = V / (n * D) Kt = T / (rho * n**2 * D**4) Kq = Q / (rho * n**2 * D**5) eta = J * Kt / Kq # 收敛性检查：尾部波动幅度 spread = (max(r[1] for r in tail) - min(r[1] for r in tail)) / abs(T) print(f"J={J:.3f} T={T:.2f} N Q={Q:.3f} N·m") print(f"Kt={Kt:.4f} Kq={Kq:.5f} eta={eta:.3f} 推力波动={spread:.1%}")

代码里几个地方值得多说一句。取尾部 20% 做平均而不是取最后一个值，是因为定常解的尾部仍有小幅振荡，单点取值随机性太大。Kt 和 Kq 的分母分别是 ρn²D⁴ 和 ρn²D⁵，指数不能错，D 的四次方和五次方在 D 只有 0.25 时数值很小，手算容易漏掉量级。eta 用 J·Kt/Kq 而不是 T·V/(2πnQ)，两个式子等价，但前者的无量纲形式便于和其他算例横向比较。最后那行推力波动是关键，超过 3% 就说明还没收敛，别看残差曲线已经降到 10⁻⁴ 就以为万事大吉。

5. 常见问题排查：五类最容易翻车的仿真现象

5.1 推力算出负值或者量级离谱

现象很直观：T 是负的，或者比预期大一个数量级。最常见的原因是桨叶的旋转方向与几何的螺旋方向反了，桨叶在"倒转"，产生的自然是反向力。检查方法是看压力云图，正常工况下吸力面应该是负压、压力面正压，如果两张面反了，就是旋向问题。第二个原因是力的方向取错，目标里取的是 X 方向力，但模型坐标系的轴向是 Z，取错了方向得到的是侧向力，数值小且方向漂移。第三个原因相对少见但更隐蔽：单位没统一，几何用毫米建模、流体域按米输入，D 在公式里差了 1000 倍。

解决方式很机械：先固定一个已知工况，把旋向、轴向、单位三项逐一核对，用理论估算值做参照，Kt 的合理区间在 0.1 到 0.3 之间，超出这个范围基本可以判定设置有误。

5.2 残差一路下降，目标值却一直在漂

残差曲线很漂亮，推力曲线还在缓慢爬升或者持续下探，这是典型的"假收敛"。原因通常是网格太粗，叶梢附近的流动结构一直在被网格重排，求解器在追一个不存在的定常解；或者流体域出口太近，尾流被压力边界反复反射。还有一种情况是时间步长或松弛因子给得太激进，稳态求解在两个解之间来回跳。

处理顺序是：先看流体域尺寸是否满足 4D 入口、10D 出口的要求；再把全局网格级别提一级，看推力曲线是否变平；最后才动松弛因子。如果三套网格的推力都在漂，那就是边界条件的问题，回头检查入口速度和出口压力是否自洽。

5.3 网格生成直接失败或者单元数爆炸

网格生成失败通常卡在两个位置：旋转域和静止域的交接面有缝隙或者重叠，以及桨叶根部和桨毂的过渡圆角被切成了薄片。前者是建模问题，两个域必须共用同一个圆柱面，不能各画一个再布尔；后者是圆角半径小于最小缝隙尺寸，部分单元被切成接近零体积的碎片。

单元数爆炸则是另一个极端，往往是局部细化区域给得太大或者最小壁面厚度给得太小。细化区域做成包住叶梢的小圆柱就够了，不要沿着整个桨叶铺一层细化盒。最小壁面厚度按实际加工粗糙度和边界层需求给，课设里给 0.1 到 0.2 mm 量级足够，再小就是在给壁面铺十几层网格，收益很低。

5.4 加了旋转区域，结果和没加一样

这种情况最让人怀疑人生。原因之一是旋转区域没有真正包住桨叶，桨叶有一截露在旋转域外面，那一截还在静止参考系里，等于部分桨叶没转。检查方法是把旋转域半透明显示，看桨叶是否完整落在里面，尤其是叶梢和根部。原因之二是旋转轴的方向和位置给错，轴心偏离桨毂中心，旋转域在"绕着别的地方转"，效果接近于没转。原因之三是角速度单位搞混，把 rpm 当成 rad/s 填进去，转速只有实际的十分之一，推力的变化小到看不出来。

核对顺序是：旋转域包住桨叶 → 轴线穿过桨毂中心且方向沿轴向 → 角速度用 rad/s 且数值等于 2πn/60。三项逐一确认，基本能一次排掉。

5.5 同一个模型两次算出来结果不一样

可复现性问题是隐性的，但答辩时最容易被问。常见来源有三个：网格是自动生成的，每次生成的单元数和分布有微小差异，尤其是部分单元的处理顺序；求解的初始条件没有固定，第二次算的时候用了上一次的结果做初值；以及并行求解的线程数不同导致浮点累加顺序变化。前两个影响大，第三个通常只在小数点后第三位。

做法是把网格文件保存下来复用，不要每次重新生成；每次算例都从统一的初始条件开始，别偷懒用上次的结果续算；并行线程数固定。做到这三条，同一个算例两次的推力偏差可以压到 0.5% 以内，答辩时被问到也能拿得出证据。

6. 验收与答辩：怎么证明这份仿真结果可信

仿真做完不算完，能证明它对才算完。验收的核心逻辑是"多重交叉验证"，任何单一证据都不够。我一般按下面这张表逐项过，每一项都要能拿出文件或者截图。

表格里"趋势正确"这一项最容易被忽略，但它恰恰是答辩时最有说服力的证据。单点数值再准，也可能是一次巧合；把 J 从 0.4 扫到 1.0，Kt 单调下降、η 在 J≈0.7 附近出现峰值，这个形状对了，说明整套流程的物理逻辑是通的，评委也更愿意相信你的数值。

答辩时讲什么，我的建议是把失败案例也放进去。只讲"我算出来 Kt=0.21"，评委会追问"你怎么知道不是 0.15"；如果你能补一句"第一版网格没加密，算出来是 0.18，加密后稳定在 0.21，粗中细三套网格的偏差 1.4%"，可信度立刻不一样。误差来源也要主动交代：网格分辨率、湍流模型对分离流动的近似、边界条件与真实水池的差异，各占多少说不清没关系，能指出方向就说明你想清楚了。

最后说个习惯。早些年我做第一个螺旋桨算例，收敛曲线漂亮得不行，兴冲冲拿去汇报，被问了一句"你入口速度是几"，才发现自己填的是 0.4 m/s 而不是 4 m/s，J 差了十倍，整条曲线全部作废。从那以后我每次开算之前都强制走一遍"参数表核对—单位核对—旋向核对"的三步检查，把这三行写在算例文件夹的 readme 里，改一次参数就重读一次。这套流程看着笨，但省下来的返工时间远比它多。希望帮到你。

本文还有配套的精品资源，点击获取