基于SimV&Ver Static的网格收敛性研究
问题背景
随着有限元技术的发展以及计算机技术的发展，有限元分析已经走向成熟，但是随着面对的问题越来越复杂，对仿真模型的计算精度要求越来越高，计算模型的可信度确认也会越来越困难，特别是在经验基础上的一些建模方式可能不再试用。数值误差包括建模误差、离散误差以及程序误差，其中建模误差是从物理现象转为数学问题而引入的，离散误差主要取决于网格的疏密程度，网格的疏密程度成为了仿真的精度的关键影响因素，也包括了静力学仿真分析的精度。
目的
国际上流行的方法是由Roache提出的，即网格收敛指数GCI（grid convergence index）。安怀信公司的静力学验模软件SimV&Ver static中集成这一方法，目的就是为了用户使用该方法快速的找到满足用户精度要求的网格尺寸，即保证精度的同时也最大程度的节省的计算的资源。
GCI方法
GCI，就是用于3套及以上系统细化网格，基于Richardson外推法估算网格收敛误差，是一种误差范围估计。
假定有3套网格及对应的输出响应量值如下：
粗网格 | 中网格 | 细网格 | |
结构域体积 | V | V | V |
单元数量 | N1 | N2 | N3 |
响应量值 | w1 | w2 | w3 |
I) 首先计算特征网格尺寸
数值仿真
采用的模型是一个存在在一定假设条件下的解析解的悬臂梁模型采用的模型，悬臂梁上部受压力10MPa，一边端面采用固定端，梁模型是正方形截面，边长为20mm，长为300mm，如图1所示。本次利用数值仿真的悬臂梁最大位移值来作为网格收敛性检验的情况，即随着网格细化是否数值计算的悬臂梁最大位移值会有一个收敛的趋势。根据梁自由端面的计算公式，悬臂梁最大挠度计算公式为：
网格收敛性分析
为了评估数值计算误差，对数值解进行网格收敛性分析，通过SimV&Ver static验模软件来对有限元模型网格进行细化，通过软件来计算p和GCI的值。
软件操作流程如下图所示：
具体的SimV&Ver static软件界面上的步骤如图4至图10所示。
将最终计算的结果如下表所示。
表1 悬臂梁网格收敛性分析
尺寸 mm | 网格 节点数 | 安全系数p | GCI （10-5） | 测点仿真 结果 mm | 挠度最大解析值mm | 误差 % | 外插值 mm |
10 | 782 | 2.4791 | 6.9644 | 72.0816 | 72.3214 | 0.331575 | 72.2067 |
5 | 5754 | 72.1843 | 0.189570 | ||||
2.5 | 60161 | 72.2027 | 0.164128 | ||||
注：测点选择应该尽量靠近悬臂梁自由端截面最大位移处的位置，为了与解析解方便比较。 | |||||||
结论
a) 悬臂梁施加均布力载荷随着网格尺寸的减小，测点结果越趋近于一个值，这个值越接近于解析解的数值；
b) 静力学验模软件SimV&Ver static可有效地运用GCI的方法为静力学数值仿真模型快速找到一个满足精度要求得网格尺寸。
欢迎给“安怀信正向设计研发港”***标星🌟，点赞。
--------------------------------------------------
如有内容来源于互联网,版权归原作者所有,供学习交流使用,严禁商用,如有侵权请联系我们删除。
--------------------------------------------------
