这个系列我会持续的说一些Abaqus求解方面的小问题,都是我日常遇到的,帮助大家避坑的同时,也让大家感受学习有限元算法的乐趣,大家一起来“祛魅”,Abaqus虽然是一款成熟的产品,但是也不免有些小问题,我等“小辈”可以发现也或可超越。
下图为一个实体单元组成的案例,案例中左端固定,右端施加点力,左侧案例采用C3D20,而右侧单元采用C3D20R (缩减积分),故意设置的单元长细比较大,但未达到报错的地步,计算出来的位移结果如下图所示:
可以看到右侧已经出现明显的沙漏现象,一般情况下二阶单元其实是对沙漏效应有一定抵抗能力的,然后当大家想调一调沙漏控制方法的时候,你会发现Abaqus对其二阶单元抵抗沙漏特别有自信,并没有给你提供沙漏控制的选项(一阶是有的),如下图所示: 所以弱小和无知不是生存的障碍,傲慢才是!(《三体》)。
然后我们进一步看这个案例,当我们查看其支座反力,会发现支座区域会有明显的渲染错误,其实这种情况也好解释,证明Abaqus后处理采用的还是最简单的三角面片渲染,并没有针对六面体特别是二阶单元进行一些修正,这就会导致在中心节点处经常会出现渲染错误,而由于支座反力只在支座处存在,其他地方为0,这种效应就特别明显。比如Femap和Hyperview就没有这种问题。
最后重申Abaqus是一款成熟的软件,我非常爱用,大多数问题上非常厉害(算法稳健、操作流程合理),有些小问题,大家看着当个乐,当然也不少,后续持续更新......