SUBOFF模型阻力、变攻角、变侧滑角（斜航运动）水动力计算

学习目标：SUBOFF阻力计算、变攻角计算、变侧滑角计算（斜航运动）

学习内容：

1. 使用HEXPRESS网格划分，面细化、体细化、设置边界层。
2. 计算SUBOFF body阻力，对比不同湍流模型对结果的影响。
3. 网格无关性验证。
4. 变攻角alpha计算，使用Fluent批处理技术，使用matlab编写变攻角计算脚本、编写导出水动力系数脚本、水动力系数绘图。
5. 无量纲的轴向力X系数、垂向力Z系数、纵倾力矩M系数与试验值对比。
6. 变侧滑角beta计算（斜航运动），同样使用Fluent批处理技术，使用matlab编写脚本。
7. 无量纲的轴向力X系数、横向力系数、偏航力矩N系数与试验值对比。

学习时间：

• 2023.10-11

---

学习产出：

1.网格划分

suboff body网格：
计算域：半模计算，4L2LL
面细化：head，tail，tail cap：0.01、0.005、0.005
paralle：0.02、0.01、0.01
体细化box：0.4、0.2、0.2
边界层：0.001、1.54或1.2、10
网格数：202万左右

2.Fluent计算

湍流模型：standard k-espilong
壁面函数：standard wall function
边界条件：前面和周面设置为速度入口，指定velocity magnitude和x/y/z flow direction；
出口为压力出口，指定压力出口相对压力为0；
航行器表面为为无滑移壁面，no slip；
求解方法：simple方法，二阶迎风格式。
监视阻力。

3.阻力计算结果对比

改变面细化网格和体细化尺寸，得到不同的网格数量，分别计算阻力，然后和实验值进行对比误差，找到最合适的网格数。

4.SUBOFF full变攻角alpha计算

使用matlab编写Fluent 批处理脚本，包含：变攻角（-16度—16度）、导出水动力系数txt
然后使用matlab编写从水动力系数txt中收集系数到txt，绘图。
alpha=16度的速度云图和压力云图：

alpha=0度的速度云图和压力云图：

alpha=-16度的速度云图和压力云图：

计算结果：无量纲化的力、力矩系数随攻角的变化曲线，并进行拟合。

5.SUBOFF full变侧滑角计算

重新画网格，采用全模计算，网格数翻一倍。
侧滑角：0—16度
beta=0度的速度云图和压力云图（水平面）：

beta=16度的速度云图和压力云图（水平面）：

计算结果：无量纲化的力、力矩系数随侧滑角的变化曲线，并进行拟合。


与下面这张图进行计算结果对比
SUBOFF 阻力试验结果：