解検証(Solution Verification)
解検証(Solution Verification)
理论与物理
概述
老师!今天要讲的是解验证(Solution Verification)对吧?具体是什么内容呢?
针对特定问题的离散化误差估计。通过Richardson外推和GCI进行不确定性量化。
控制方程
$$ \hat{f}_h = f_h + \frac{f_1 - f_2}{r^p - 1} $$
$$ \text{GCI}_{fine} = \frac{F_s|\varepsilon_{12}|}{r_{12}^p - 1} $$
离散化方法
这个方程,在计算机上实际是怎么求解的呢?
使用有限元法(FEM)进行空间离散化。组装单元刚度矩阵,构建整体刚度方程。
矩阵求解算法
矩阵求解算法,具体是指什么呢?
通过直接法(LU分解、Cholesky分解)或迭代法(CG法、GMRES法)求解联立方程。对于大规模问题,带预处理的迭代法非常有效。
| 解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小~中规模 |
| Cholesky分解 | 直接法(对称正定) | O(n²) | 小~中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 前处理 | O(n) | 超大规模 |
也就是说在有限元法这部分偷懒的话,后面会吃苦头对吧。我铭记在心!
商用工具中的实现
那么,要进行解验证(Solution Verification)可以用哪些软件呢?
| 工具名 | 开发商/现状 | 主要文件格式 |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software) | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| OpenFOAM | 开源(OpenCFD/ESI、OpenFOAM Foundation) | 字典文件(blockMeshDict等), .foam |
供应商谱系与产品整合历程
各个软件的诞生过程,是不是还挺有戏剧性的?
MSC Nastran / NX Nastran
接下来是MSC Nastran的内容对吧。具体是什么?
作为NASA结构分析(NASTRAN)于1960年代开发。MSC Software将其商业化,之后UGS(现Siemens)派生出NX Nastran。MSC于2017年被Hexagon AB收购。
当前所属:MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software)
Abaqus FEA (SIMULIA)
Abaqus FEA,具体是指什么呢?
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被达索系统收购,并整合至SIMULIA品牌下。
当前所属:Dassault Systèmes SIMULIA
等等等等,结构分析的话,也就是说这种案例也能用吗?
Ansys Mechanical (旧ANSYS Structural)
请讲讲「Ansys Mechanical」!
哦~,结构分析的话题,超级有趣!请再多讲一些。