数值仿真为什么需要生成网格

CAE十级退堂鼓

12月前浏览1737

数值仿真需要网格的主要原因在于其核心思想：将连续的物理问题离散化为有限数量的计算单元，以便用计算机进行求解。

网格在这一过程中扮演了关键角色，具体原因如下：

—

离散化需求

连续到离散的转换：

现实中的物理场（如流体、温度、应力等）是连续分布的，但计算机只能处理离散数据。网格将求解域划分为有限数量的单元（如三角形、四边形、六面体等），每个单元上的物理量通过节点或积分点上的值来近似表示。

控制方程的离散：

微分方程（如Navier-Stokes方程、热传导方程）需要在离散的网格上通过有限元法（FEM）、有限体积法（FVM）或有限差分法（FDM）转化为代数方程。

—

数值积分的实现

积分计算：

在有限元法中，刚度矩阵和载荷向量的计算需要数值积分（如高斯积分），而积分需要在每个网格单元内进行。网格提供了积分域的几何定义。

物理量的插值：

网格节点或单元上的变量（如位移、速度）通过形函数（Shape Functions）插值得到连续分布的解，形函数的定义依赖于网格的几何形状。

—

边界条件的施加

边界标识：

网格的边界单元或节点用于精确施加边界条件（如固定约束、热通量、压力入口等）。

无网格时，边界条件的数学描述和数值实现会变得复杂。

—

复杂几何的适配

几何逼近：

网格可以拟合复杂形状（如飞机翼型、血管结构），通过局部细化（如自适应网格）捕捉几何细节或物理场梯度（如边界层、激波），得到更加精准的边界计算结果。

多尺度问题：

通过非均匀网格（如加密网格）在关键区域（如裂纹尖端）提高计算精度，同时减少其他区域的计算量。

—

计算效率并行化

局部性：

网格单元之间的拓扑关系明确，便于构建稀疏矩阵（降低内存占用）和设计高效的迭代算法（如多重网格法）。

并行计算：

网格可分块处理（如区域分解法），适合分布式计算。

—

例外情况

无网格方法（Meshfree Methods）

虽然传统数值仿真依赖网格，但无网格方法（如光滑粒子流体动力学SPH、移动最小二乘法MLS）通过粒子或散点近似解，避免了显式网格生成。

这种技术已经开始出现，目前还没有大规模使用。

然而，这类方法目前来说计算成本更高或稳定性较差，因此网格法仍是主流。

—

总结

网格的本质是“离散化的工具”，它平衡了计算精度、效率与工程实用性。

对于大多数实际问题（如结构力学、CFD、电磁仿真），网格提供了可控制、可验证的数值框架，是连接连续物理世界与离散计算模型的桥梁。

来源：CAE十级退堂鼓

仿真模型不收敛咋办啊

在ANSYS结构仿真中，判断结果是否收敛以及解决不收敛问题是实际工程分析的核心挑战。以下是针对ANSYSWorkbench的实操指南，涵盖收敛判断方法、监控工具和不收敛时的解决方案。01—收敛性判断方法在workbench中，有如下三种方法来查看计算是否收敛：1.监控求解器输出信息在计算开始后，在solutioninformation中查看命令行消息，默认的选项为solveroutput：默认容差：力残差通常需小于0.5%（也可自定义）。出现如下消息表示收敛成功，计算完成。出现Error表示不收敛，具体可查看报错信息，判断模型的改进方向。2.图形化监测工具在计算开始后，在solutioninformation下拉菜单选择ForceConvergence，可观察力的收敛残差曲线。静力学计算的不收敛的情况多半如下图：力的收敛曲线在残差标准线之上，红色虚线表示子步计算失败，继续将子步做二等分计算。静力学计算的收敛的情况如下图：计算时间稳步增加，绿色的虚线表示子步收敛。蓝色虚线表示加载步收敛。对于显式动力学，通常观察能量平衡来判断内能（StrainEnergy）、动能（KineticEnergy）、外力功（ExternalWork）应平衡。能量误差（|内能+动能-外力功|/外力功）应<5%，下图这个超过5%了，可判断为收敛性不佳。3.收敛标准的设置：通过AnalysisSettings>NonlinearControls调整收敛容差，一般采用程序默认，在自己有一定的工程经验后可根据自己的判断来适当收严或放宽。默认的收敛容差为0.5%02—解决不收敛的办法1.网格优化问题：报错信息提示网格有畸变。这种报错有一部分原因是网格的问题，也可能是边界有问题，总之第一步是先检查网格质量。操作：在提示的网格单元附近区域局部加密网格。提高单元阶数检查网格质量关于如何查找高畸变网格的位置，可以通过报错的单元ID来查找2.接触设置调整问题：接触力不收敛或穿透过大。操作：修改接触算法：在接触对的Details中，将Formulation从PurePenalty改为AugmentedLagrange。勾选NormalStiffness=Aggressive（提高接触刚度）。启用自动时间步：在AnalysisSettings中，设置AutoTimeStepping=On，InitialSubsteps=100。调整接触探测：在接触Details中，设置DetectionMethod=Nodal-NormaltoTarget（减少穿透）。3.载荷步与子步控制问题：材料非线性或大变形导致发散。操作：分步加载：在AnalysisSettings中，设置NumberofSteps=2，第一个步施加50%载荷，第二步加载至100%。增加子步数：在AnalysisSettings中，设置InitialSubsteps=50MinimumSubsteps=20MaximumSubsteps=1000启用大变形：在AnalysisSettings中勾选LargeDeflection=On（适用于几何非线性）。4.约束与稳定性问题：刚体位移或数值不稳定。操作：添加弱弹簧：在AnalysisSetting中找到Stabilization，选择WeakSprings为On。检查约束：确保所有自由度的刚体运动被限制（如对称模型需施加对称约束）。5.数值阻尼（适用于振荡）操作：调整阻尼系数：在Stabilization下设置DampingFactor=0.1（默认0.2，值越小越严格）。以上方法为，确定边界条件无误的前提下进行的。在不收敛时第一步应该先检查边界条件简化的是否合理！以上为上周调试大变形模型而引发的思考。如有不足，评论区指正！来源：CAE十级退堂鼓