球的抗力（Stokes/Oseen解）

低雷诺数的球周围流动。斯托克斯阻力规律和Oseen修正。粒子追踪模型的验证。

采用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，建立整体刚度方程。

将微分方程转换成弱形式（变分形式），使用试函数和形状函数进行Galerkin方法定义。单元类型的选择（低阶单元 vs. 高阶单元、完全积分 vs. 减缩积分）直接影响解的精度和计算成本的权衡。

直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解线性方程组。大规模问题采用前处理的迭代法效率更高。

NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。

在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表示差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）往往不能在求解器之间直接转换。

网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析非常重要。

是的，不错的状态！亲自动手实践是最好的学习。有任何疑问随时来问我。

使用形状函数 $N_i$ 来近似未知量：

用公式表示就是这样。

用公式表示就是这样。

连续体的控制方程经过离散化后，得到下面这样的代数方程组：

这里 $[K]$ 是整体刚度矩阵（或等价的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。

收敛速度：二阶单元以 $O(h^2)$ 的量级减小误差（光滑解的情况）

说明球的抗力（Stokes/Oseen解）的实务解析流程和注意事项。

1. 前处理 (Pre-processing)

2. 求解 (Solving)

3. 后处理 (Post-processing)

NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

是的，不错的状态！亲自动手实践是最好的学习。有任何疑问随时来问我。

对应球的抗力（Stokes/Oseen解）的主要商用CAE工具的功能比较，以及各产品的历史背景进行详述。

NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，整合到SIMULIA品牌。

Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。