如果满足稳定性约束,FEM 和 FVM 代码收敛到相同物理问题(相同的 BC、几何、属性等)的略微不同的解决方案是否合理?甚至,不同的 FEM 方法是否有可能收敛到不同的解决方案?
这应该是一个更普遍的问题,但它受到我几周前测试的一个具体示例的启发。
几周前,我尝试在 Abaqus、ANSYS 和 FVM 求解器中建模一个简单的热膨胀问题,所有这些都在线性弹性下运行。我很确定所有 3 个求解器中的条件都相同,但似乎 3 个求解器收敛到略有不同的结果,差异约为 10-20%(Abaqus 与 ANSYS 也存在差异,2 个 FEM 求解器)。
我自己编写了 FVM 求解器,并且验证了线性弹性方程得到了正确求解,并且我假设 ANSYS 和 Abaqus 等已建立的求解器也正确求解了它们的方程。这让我很困惑 3 个求解器如何各自收敛到不同的解决方案,我无法直观地想到这怎么可能,在数字上?