今天来聊一个仿真工程师们最不想看到的一个情况:不收敛。明明跟着做,怎么就在求解时出现了不收敛呢?这是什么原因导致的?

其实这背后还是有很多小技巧的,即使你在网上搜,也不一定能够针对你遇到的问题给出解答。下面我们元王仿真就来详细讲解,搞明白“收敛”,才知道为啥不收敛。
01、为什么非得"迭代"
有限元的本质,是把一个解不出来的问题拆成一堆能解的小问题,之后再拼成一整块。而在操作中,麻烦就在"拼"这一步。
每个单元都连着周围的单元,A的变形取决于B,B又反过来受A影响。这种互相牵制让方程变成非线性的,一次算不出来。
而如果不熟悉,就需要一轮又一轮的调整迭代,将差值缩小,“差”越来越小就叫收敛。
举个生活中很常见的例子,就是有时我们看乐手在调乐器,他都会先弹一下,听着差半个音,拧一点;再听,再拧,这个过程就是将音调正的过程。
02、三大准则:三个角度的"验收"
①力收敛:内外力平衡了没
软件算出的内力(结构内部各单元扛住的力)要和外部施加的载荷对得上。两者的差值叫不平衡力,也叫残差。不平衡力相对于载荷总量小到某个比例以下,就算收敛。
这是多数结构分析软件的默认准则,我们要解的本质就是力的平衡方程,力不平衡,结果就没有物理意义。

看曲线时注意:它应该往下走,而且持续往下走。如果在某个水平上反复震荡,说明模型里有东西在跟求解器"对着干",多迭代几轮也解决不了。
②位移收敛:结果是不是稳住了
比较这一轮和上一轮的位移,各节点变化量都很小,说明结果已经稳定。有个很实用的经验:位移比应力更容易收敛。
位移是积分量,是整个结构变形的宏观体现,局部波动会被平均掉,所以它很容易掌控。应力是位移求导来的量,等于把局部差异放大了,天生比位移难控制。这就是为什么有时候位移曲线早平了,应力曲线还在晃。
这时候其实不是软件出问题,是这两类物理量的数学性质决定的。
③能量收敛:整体还有没有改进空间

从"功"的角度看,每轮迭代都会改变系统能量状态,能量增量小到可忽略,说明这轮没有实质改进。涉及大变形时用得较多,不好直观理解,但对判断整体是否稳定相当敏感。
这三准则该怎么选?不用纠结。三者并不互斥,很多软件允许组合勾选。实际工作中默认的力收敛准则已经能满足绝大多数静力分析,仿真工程师要做的是看懂曲线。
03、收敛公差
公差越小,结果越精确,但迭代变多、时间变长,太小了还可能永远达不到,直接报不收敛。
新手最容易犯的错就是把公差使劲去调小,觉得"越严越专业",结果算一晚上没出来,还以为是模型有问题。够用就好,结构分析里默认的1%到0.1%量级,对绝大多数工程判断已经足够。
如果结果对公差极其敏感,从1%改到0.1%答案就变了,那问题不在公差,在模型:可能有应力奇异,可能接触没调好。
继续收紧公差,只是把问题往后推。
04、不收敛了,按这五步查
接下来这段是本篇文章最重要且需要记下来的一段,遇到不收敛该怎么办?别慌!一步步来自查,也能轻松把问题解决。多遇到几次,就熟能生巧。
如果在日常工作中,想提升仿真技能,也可以找我们元王仿真咨询,培训这块我们一直在做,无论是个人还是企业,都有适合的课程。下面接着继续分享步骤。
第一步:看曲线,判断卡在哪。第一轮就不动,大概率模型有硬伤(约束不足、刚体位移、单位制错);算到一半崩掉,多半是非线性因素触发了问题(材料屈服、接触跳变、大变形);曲线反复震荡,通常是接触或摩擦设置有问题。曲线形态本身就是最直接的诊断线索。
第二步:退回线性做对照。材料改纯弹性、去掉接触、关掉大变形,先跑一个最干净的线弹性分析。能过,说明问题在你加的非线性因素上;不能过,说明基础就错了,别浪费时间调参数。
第三步:一项项加回来。线弹性通过后,别一次全开:先加塑性,收敛了再加接触,再加大变形。哪一步开始不收敛,问题就在哪一步。是最笨的办法也是最可靠。
第四步:排查五大高频诱因。
网格质量太差:扭曲畸变、长宽比离谱的单元会让刚度矩阵条件数变得极差,求解器算不动。四面体网格在薄壁尖角处自动生成的畸变单元,数量不多但杀伤力很大。
约束不足:结构能整体平移或转动,平衡方程没有唯一解。
接触设置不当:接触状态在"接触"和"分离"之间跳变,接触刚度过大易震荡、过小会穿透。先确认接触对定义对了,再谈调参数。
材料参数不合理:最常见的是只给弹性模量就去做塑性分析,线弹性能过,一进塑性就崩。
单位制没统一:毫米当成米、力用成吨,刚度离谱到求解器无法处理。看似离谱,但对于粗心的人来说都有可能。在开头统一单位,就能省掉一半莫名其妙的问题。
第五步:最后才动求解器。时间步长、自动时间步、线性搜索、弧长法。顺序不能反。这些设置是用来"帮模型算下去"的,不是用来"修模型错误"的。硬算一个本来就错的模型,只会得到一份看起来很漂亮的错误报告。
05、收敛了,不代表算对了
其实你还要知道一点,就是收敛只是数学上的成功,不是物理上的正确。
软件说不平衡力已经足够小,它说的是"我按你给的方程算到位了",不是"你给的方程是对的"。边界条件设错的模型可以收敛得很好看,漏了载荷的模型照样跑出漂亮云图......而且这份垃圾还带着收敛报告,格外有说服力。
所以每次算完,问自己三个问题:变形趋势符合力学常识吗?应力集中的位置跟以往失效案例对得上吗?如果拿去做试验,我有几成把握?
答不上来,就回去再查模型。收敛曲线不骗人,但它只对你给的那个模型负责。
收敛就是软件在问你一句:"我这样算,你觉得对不?"
力收敛问"平衡了没",位移问"稳住了没",能量问"还有没有改进空间"。
作为仿真工程师,能做的就三件:模型建对(往往八成不收敛是这个原因),公差设得合理,出问题按顺序查而不是乱动求解器。
下次要是再遇到不收敛问题,别急着直接关软件。按照上面的方法来进行,一步步排查,学会了就不再觉得是难题了。
如果这篇内容有用就点个赞或者分享给更多的人,一起学习有限元分析知识点,一起提升有限元分析技能。元王仿真(有限元科技)您身边的CAE仿真顾问,有问题就找我们。
下一篇:没有了