船舶抵抗预测
船舶抵抗预测的理论基础
概要
老师！今天是船舶抵抗预测的话题吧？这是什么东西呢？
ITTC1957摩擦抵抗式。造波抵抗的CFD预测。模型船和实船的相关性。
控制方程
哇～，船舶抵抗预测的描述话题，超级有意思！请告诉我更多。
离散化方法
这个方程在计算机上到底怎么求解呢？
使用有限元法(FEM)进行空间离散化。组装单元刚性矩阵，构建整体刚性方程。
矩阵求解算法
矩阵求解算法，具体是怎么回事呢？
通过直接法(LU分解、Cholesky分解)或迭代法(CG法、GMRES法)求解联立方程组。对于大规模问题，预条件迭代法很有效。
| 求解法 | 分类 | 内存用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小~中规模 |
| Cholesky分解 | 直接法(对称正定) | O(n²) | 小~中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，在有限元法的地方偷工减料的话，之后就会很吃亏啊。我会牢记的！
商用工具中的实现
那么，做船舶抵抗预测可以用什么软件呢？
| 工具名称 | 开发商/现在 | 主要文件格式 |
|---|---|---|
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| OpenFOAM | 开源(OpenCFD/ESI、OpenFOAM Foundation) | 词典文件(blockMeshDict等), .foam |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
厂商血统与产品整合历史
各个软件的发展历程，好像都挺有戏剧性的呢？
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题吧。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是特点。
现在的所属: Siemens Digital Industries Software
Ansys Fluent
接下来是Ansys Fluent的话题吧。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。非结构网格型通用CFD求解器。
现在的所属: Ansys Inc.
听到这里，终于明白了为什么开发这么重要！
OpenFOAM
OpenFOAM，具体是怎么回事呢？
起源于Imperial College London的开源CFD。OpenCFD Ltd（ESI Group旗下）和The OpenFOAM Foundation并行开发。
现在的所属: 开源(OpenCFD/ESI、OpenFOAM Foundation)
哇～，开发话题，超级有意思！请告诉我更多。
文件格式与互操作性
不同软件之间交换数据时有什么要注意的吗？
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| CGNS | .cgns | CFD数据 | CFD通用符号系统。CFD结果的标准交换格式。 |
| VTK | .vtk/.vtu | 可视化 | Visualization Toolkit格式。供ParaView等使用。 |
在不同求解器间转换模型时，要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表示差异。特别是高阶单元或特殊单元(内聚单元、用户定义单元等)通常无法直接在求解器间转换。
我明白了…格式看似简单，但其实深度很大呢。
实务注意事项
教科书里没有的"现场智慧"这样的东西有吗？
网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析特别重要。
嗯，你进展不错！实际动手是最好的学习方式。有什么不明白的地方就随时问我。
弗劳德的"相似律"——19世纪的天才为什么至今仍支配着船舶设计
威廉·弗劳德(William Froude)在1870年代根据模型船水槽试验建立了预测实船抵抗的"弗劳德相似律"。他的想法很简单："如果将摩擦抵抗和造波抵抗分离，就可以分别用不同的相似律进行尺度变换"。这个概念至今仍活在CFD验证框架中，因为模型试验无法同时满足雷诺数和弗劳德数的相似性，所以由CFD分析来补充的方法论的背景就源于此。150年前的物理至今仍在规定着计算设计——这是令人感受到技术连续性的有趣话题。
船舶抵抗预测的数值计算方法
数值方法的详细内容
具体用什么算法来求解船舶抵抗预测呢？
离散化的公式化
使用形状函数 $N_i$ 近似未知量:
用公式表示就是这样。
基本方程的离散形式
用公式表示就是这样。
嗯…只看公式还是有点摸不着头脑…这代表什么呢？
将连续体的控制方程离散化后，得到以下代数方程组:
其中 $[K]$ 是整体刚性矩阵(或等价的系统矩阵)，$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
啊，这样啊！连续体的控制方程变成这样的结构…我终于明白了。
单元技术
我听说过"单元技术"，但可能没有真正理解…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低(剪切锁定) | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中~高 | 中 |
积分方案
积分方案，具体是怎么回事呢？
听到这里，终于明白了为什么单元类型这么重要！
收敛性与稳定性
如果不收敛，首先要检查什么？
收敛速度: 二阶单元的误差按 $O(h^2)$ 的数量级减小(对光滑解)
我明白了…网格细化看似简单，但其实深度很大呢。
求解器设置建议
具体用什么算法来求解船舶抵抗预测呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) or AMG | 根据问题规模 |
| 最大迭代次数 | 1000 | 不收敛时应重新审视设置 |
| 内存模式 | In-core | 尽可能使用 |
风上差分（Upwind）
一阶风上: 数值扩散较大但稳定。二阶风上: 精度提高但有振荡风险。高雷诺数流动中必需。
中心差分（Central Differencing）
二阶精度，但当Pe数 > 2时出现数值振荡。适用于低雷诺数的扩散控制流。
TVD方案（MUSCL、QUICK等）
用限制函数抑制数值振荡同时保持高精度。对冲击波和陡梯度的捕捉有效。
有限体积法 vs 有限元法
FVM: 自然满足守恒律。CFD的主流。FEM: 复杂形状、多物理耦合更有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式法: CFL ≤ 1为稳定条件。隐式法: CFL > 1也稳定，但影响精度和迭代次数。LES: 建议 CFL ≈ 1。物理意义: 一个时间步长内信息不传播超过一个单元。
残差监测
连续方程、动量、能量的各残差下降3~4位数视为收敛。质量守恒的残差特别重要。
松弛因子
压力: 0.2~0.3、速度: 0.5~0.7是一般初始值。发散时降低松弛因子。收敛后提高以加速。
非定常计算的内部迭代
每个时间步内迭代到定常解。内部迭代数: 5~20次左右为目安。时间步间残差波动时应重新检视时间步。
船舶抵抗预测的实务应用
实践指南
老师，"实践指南"请教一下！
讲解船舶抵抗预测的实务分析流程和注意事项。
我明白了…实务的说法看似简单，但其实深度很大呢。
分析流程
从最初的第一步开始请教我！应该从什么开始呢？
1. 预处理 (Pre-processing)
- CAD数据的导入和形状简化
- 材料特性的定义
- 网格生成(单元类型·尺寸的确定)
- 边界条件和荷载条件的设置
2. 求解 (Solving)
- 求解器设置(解法、收敛基准、输出控制)
- 作业提交和计算执行
- 收敛监测
3. 后处理 (Post-processing)
- 结果可视化(位移、应力、其他物理量)
- 结果验证和妥当性确认
- 报告作成
网格生成最佳实践
网格的好坏怎么判断呢？
单元品质指标
"单元品质指标"请教一下！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度低下 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度低下 |
| 歪斜度 | 0° | < 45° | 收敛性恶化 |
| 锥度比 | 0 | < 0.5 | 精度低下 |
网格密度的确定
网格密度的确定，具体是怎么回事呢？
边界条件设置指南
据说边界条件如果搞错，全部都会坏掉…
啊，这样啊！过度约束的注意…这样的结构啊…我终于明白了。
各商用工具的实现步骤
有很多软件吧？各个的特点请教一下！
| 工具名称 | 开发商/现在 | 主要文件格式 |
|---|---|---|
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| OpenFOAM | 开源(OpenCFD/ESI、OpenFOAM Foundation) | 词典文件(blockMeshDict等), .foam |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题吧。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是特点。
现在的所属: Siemens Digital Industries Software
Ansys Fluent
接下来是Ansys Fluent的话题吧。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。非结构网格型通用CFD求解器。
现在的所属: Ansys Inc.
老师的说明很清楚！工具名的疑惑消散了。
常见失败与对策
初学者容易犯什么样的失败？事先想要知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不适当的边界条件 | 网格改善、约束条件重新审视 |
| 应力异常大 | 应力奇点、网格依赖 | 奇点回避、局部网格细化 |
| 位移不现实 | 材料定数误差、单位系混在 | 输入数据确认 |
| 计算时间过长 | 不必要的细化、非效率的解法 | 网格最优化、并行计算 |
质量保证检查清单
教科书里没有的"现场智慧"这样的东西有吗？
嗯，你进展不错！实际动手是最好的学习方式。有什么不明白的地方就随时问我。
水槽试验 vs CFD——"选择相信哪一个"的现场冲突
在船舶抵抗实践分析中必然出现的讨论就是"与水槽试验的偏离"。在良好条件下CFD和水槽试验的全抵抗系数误差在2~5%以内，但在附加物(舵、舭鼍)安装的接近实船的模型上能出现超过10%的偏差。有趣的是"哪个是正确的"的问题并不重要，追踪"为什么不同"会深化对模型的理解。CFD有边界条件和乱流模型的影响，水槽试验有尺度效应和波反射的影响——了解了各自的弱点后如何取舍使用就是实践最佳实践。
船舶抵抗预测的软件比较
商用工具比较
有很多软件吧？各个的特点请教一下！
详述支持船舶抵抗预测的主要商用CAE工具的功能比较，以及各产品的历史背景。
也就是说，船舶抵抗预测的支持地方偷工减料的话，之后就会很吃亏啊。我会牢记的！
支持工具列表
那么，做船舶抵抗预测可以用什么软件呢？
| 工具名称 | 开发商/现在 | 主要文件格式 |
|---|---|---|
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| OpenFOAM | 开源(OpenCFD/ESI、OpenFOAM Foundation) | 词典文件(blockMeshDict等), .foam |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题吧。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是特点。
现在的所属: Siemens Digital Industries Software
Ansys Fluent
接下来是Ansys Fluent的话题吧。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。非结构网格型通用CFD求解器。
现在的所属: Ansys Inc.
听到这里，终于明白了为什么开发这么重要！
OpenFOAM
OpenFOAM，具体是怎么回事呢？
起源于Imperial College London的开源CFD。OpenCFD Ltd（ESI Group旗下）和The OpenFOAM Foundation并行开发。
现在的所属: 开源(OpenCFD/ESI、OpenFOAM Foundation)
Ansys CFX
"Ansys CFX"请教一下！
由AEA Technology (UK)开发的CFX。2003年被Ansys收购。耦合型求解器是特点。
现在的所属: Ansys Inc.
啊，这样啊！开发…这样的结构…我终于明白了。
功能比较矩阵
预算也有限时间也有限，性价比最强的是哪个呢？
| 功能 | Star-CCM+ | Fluent | OpenFOAM | CFX |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险，具体是怎么回事呢？
啊，这样啊！不同工具间的模型…这样的结构…我终于明白了。
许可形式
我听说过"许可形式"，但可能没有真正理解…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持是有偿的 |
| COMSOL | 节点锁定/浮动 | 按模块单位购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
最后到底选哪个，判断基准请教我？
船舶抵抗预测工具选择时考虑以下方面:
嗯，你进展不错！实际动手是最好的学习方式。有什么不明白的地方就随时问我。
船体抵抗CFD工具对比——SHIPFLOW vs StarCCM+ vs OpenFOAM的ROI
船体抵抗(造波抵抗+粘性摩擦抵抗)的CFD分析中存在海洋·造船专用工具与通用CFD的分工。SHIPFLOW(FLOWTECH)是专为船体抵抗分析优化的专用工具，具有自动化表面波解析和摩擦抵抗修正，多用于新造船设计验证。StarCCM+和ANSYS Fluent作为通用CFD，用自由表面(VOF法)+RANS分析来分析船体抵抗，强项在于详细流场分析(包括推进器)和附加物(舵、舭鼍)的影响评价。OpenFOAM的shipHydroFoam在学术研究中有实绩，但商业项目的认证实绩还很少。造船所通常采用SHIPFLOW用于向监管部门提交计算，StarCCM+用于研究和优化的2工具并用方式。
船舶抵抗预测的先进研究
前沿课题与研究动向
船舶抵抗预测领域接下来怎么发展呢？
看一下船舶抵抗预测领域的最新研究动向和先进手法。
我明白了…船舶抵抗预测中的说法看似简单，但其实深度很大呢。
最新的数值方法
接下来是最新数值方法的话题吧。内容是什么呢？
嗯…只看公式还是有点摸不着头脑…这代表什么呢？
对高性能计算(HPC)的支持
| 并行化方法 | 概述 | 适用求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 许多求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别在显式法中有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
船舶抵抗预测故障处理
故障处理
常见错误与对策
老师也在船舶抵抗预测中熬过夜间调试吗？(笑)
1. 收敛失败
收敛失败，具体是怎么回事呢？
现象: 求解器在指定迭代次数内不收敛，异常结束
可能的原因:
- 网格品质不足(过度扭曲的单元)
- 材料参数设置不当
- 初始条件不当
- 非线性性太强(荷载步不足)
对策:
- 进行网格品质检查(纵横比、雅可比比)
- 确认材料参数的单位系
- 将荷载分成多个步(增加子步数)
- 放宽收敛判定基准(但要注意精度)
也就是说，在收敛失败的地方偷工减料的话，之后就会很吃亏啊。我会牢记的！
2. 非物理的结果
接下来是非物理结果的话题吧。内容是什么呢？
现象: 应力/位移/温度等不符合物理规律的非现实值
可能的原因:
- 边界条件设置误差
- 单位系混在(SI单位和工程单位混用)
- 单元类型选择不当
- 应力奇点的存在
对策:
- 确认反力总和(力的平衡)
- 确认单位系的一