CalculiX入門
理论与物理
（理论与物理部分）
数值解法与实现
数值方法详情
具体是用什么算法来求解CalculiX入门问题的呢？
我来讲解CalculiX入门的数值解法与实现要点。
原来如此。那是不是只要掌握了入门的数值解法与实现，就基本没问题了？
编译与构建
我听说过“编译与构建”，但可能理解得不是很透彻…
也就是说，在从源代码构建这一步偷懒的话，后面会吃苦头对吧。我记住了！
输入文件构成
在不同软件之间传递数据时有什么注意事项吗？
理解案例文件的结构和主要参数设置是实现的第一步。字典文件（dict）或命令文件的格式是各软件特有的，从官方教程的模板开始编辑是高效的方法。
脚本自动化
我听说过“脚本自动化”，但可能理解得不是很透彻…
使用Python或Bash脚本实现参数研究的自动化是提高生产力的关键。也应考虑利用PyFoam、cfMesh等封装工具。
调试与开发环境
使用GDB、Valgrind、AddressSanitizer进行内存泄漏检测和调试是有效的。利用IDE（VSCode, CLion）的远程调试功能，构建高效的开发环境。引入单元测试框架（Google Test, pytest），自动化回归测试。
求解器设置与算法
我想更详细地了解计算背后发生了什么！
OpenFOAM 求解器选择指南
CalculiX 输入文件结构
输入文件结构具体是指什么？
```
*NODE
1, 0.0, 0.0, 0.0
...
*ELEMENT, TYPE=C3D8
1, 1, 2, 3, 4, 5, 6, 7, 8
...
*MATERIAL, NAME=STEEL
*ELASTIC
210000., 0.3
*DENSITY
7.85e-9
*BOUNDARY
1, 1, 3
*CLOAD
100, 2, 1000.
啊，原来是这样！求解器选择指南原来是这样的机制啊。
Code_Aster 命令文件结构
接下来是命令文件结构的讲解吧。内容是什么？
```
DEBUT()
MAIL = LIRE_MAILLAGE()
MODELE = AFFE_MODELE(MAILLAGE=MAIL, ...)
RESULT = MECA_STATIQUE(MODELE=MODELE, ...)
FIN()
```
离散化格式选择
请讲解一下“离散化格式选择”！
OpenFOAM的离散化格式在 fvSchemes 文件中设置。对流项的离散化对精度和稳定性影响很大：
听到这里，我终于明白为什么求解器选择指南如此重要了！
- upwind: 一阶精度，稳定但数值扩散大
- linearUpwind: 二阶精度，受限
- limitedLinear: 二阶精度，带TVD限制器
- LUST: 混合格式，推荐用于LES
误差评估与精度验证
我听说过“误差评估与精度验证”，但可能理解得不是很透彻…
离散化误差评估
离散化误差评估具体是指什么？
基于Richardson外推法的离散化误差估计：
其中 $f_h$ 是网格尺寸 $h$ 下的解，$r$ 是网格比，$p$ 是离散化阶数。
GCI（网格收敛指数）
请讲解一下“GCI”！
基于 ASME V&V 20-2009的网格收敛性定量评估：
听到这里，我终于明白为什么离散化误差评估如此重要了！
用公式表示就是这样。
嗯…只看公式不太明白… 这表示什么？
安全系数 $F_s = 1.25$（用于三个或更多网格等级比较时）。GCI < 5% 可作为收敛的参考标准。
我明白前辈为什么说“离散化误差评估一定要认真做”了。
验证基准问题
请讲解一下“验证基准问题”！