圆管内层流（哈根-泊肃叶流）

完全发达的圆管内层流的严密解。速度分布为抛物线，流量与R⁴成正比。CFD的基本验证问题。

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

进行弱形式（变分形式）的变换，用试验函数和形状函数进行Galerkin法的定式化。单元类型的选择（低阶单元 vs. 高阶单元、完全积分 vs. 减缩积分）直接关系到解的精度和计算成本的权衡。

用直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解线性方程组。大规模问题中预处理迭代法很有效。

1960年代作为NASA构造解析（NASTRAN）开发。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）为基础。

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

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

圆管（R=0.01m, L=0.1m）内层流。μ=0.001Pa·s, dp/dx=-100Pa/m。Re=10。

u_max = R²/(4μ)·|dp/dx| = 0.025 m/s, Q = πR⁴/(8μ)·|dp/dx| = 3.927×10⁻⁶ m³/s

二阶方案比一阶明显高精度。边界层网格（壁面附近的棱柱层）精度直接相关。

上表展示了系统改变单元类型和网格密度的结果。二阶单元相比线性单元显示出明显更快的收敛，且粗网格也能获得实用的精度。应该用GCI（Grid Convergence Index）计算离散化误差的95%置信区间定量评估。

嗯，不错啊！ 其实最好的学习就是动手去做。有问题的话随时来问。

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

用公式表示如下。

用公式表示如下。

连续体支配方程进行离散化后，得到下列代数方程组：

其中 $[K]$ 是全体刚度矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量矢量，$\{F\}$ 是外力矢量。

收敛速度：二阶单元达到 $O(h^2)$ 阶误差减少（光滑解的情况）

解说圆管内层流（哈根-泊肃叶流）的实务解析流程和注意事项。

1. 预处理 (Pre-processing)

2. 求解 (Solving)

3. 后处理 (Post-processing)

1960年代作为NASA构造解析（NASTRAN）开发。MSC Software进行商用化，之后UGS（现Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。

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

嗯，不错啊！ 其实最好的学习就是动手去做。有问题的话随时来问。

对应圆管内层流（哈根-泊肃叶流）的主要商用CAE工具的功能比较，以及各产品的历史背景详述。