首页/文章/ 详情

我看Abaqus的一些小问题-1

1天前浏览8

这个系列我会持续的说一些Abaqus求解方面的小问题,都是我日常遇到的,帮助大家避坑的同时,也让大家感受学习有限元算法的乐趣,大家一起来“祛魅”,Abaqus虽然是一款成熟的产品,但是也不免有些小问题,我等“小辈”可以发现也或可超越。

下图为一个实体单元组成的案例,案例中左端固定,右端施加点力,左侧案例采用C3D20,而右侧单元采用C3D20R (缩减积分),故意设置的单元长细比较大,但未达到报错的地步,计算出来的位移结果如下图所示:

可以看到右侧已经出现明显的沙漏现象,一般情况下二阶单元其实是对沙漏效应有一定抵抗能力的,然后当大家想调一调沙漏控制方法的时候,你会发现Abaqus对其二阶单元抵抗沙漏特别有自信,并没有给你提供沙漏控制的选项(一阶是有的),如下图所示: 所以弱小和无知不是生存的障碍,傲慢才是!(《三体》)。

然后我们进一步看这个案例,当我们查看其支座反力,会发现支座区域会有明显的渲染错误,其实这种情况也好解释,证明Abaqus后处理采用的还是最简单的三角面片渲染,并没有针对六面体特别是二阶单元进行一些修正,这就会导致在中心节点处经常会出现渲染错误,而由于支座反力只在支座处存在,其他地方为0,这种效应就特别明显。比如Femap和Hyperview就没有这种问题。

最后重申Abaqus是一款成熟的软件,我非常爱用,大多数问题上非常厉害(算法稳健、操作流程合理),有些小问题,大家看着当个乐,当然也不少,后续持续更新......


来源:大狗子说数值模拟
HyperViewAbaqusFemap控制渲染
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-02-20
最近编辑:1天前
大狗子说数值模拟
博士 传播国际一流的数值模拟算法
获赞 6粉丝 6文章 45课程 0
点赞
收藏
作者推荐

有限元基础知识:线性化

接着上一次的偏导与方向导数,微分与变分,有限元基础知识:偏导、方向导数、微分与变分这次讲一下线性化,线性化顾名思义其实其实就是一种”以直代曲“的方法,英文叫Linearization,对于非线性问题的求解极为重要。大家上大学学高等数学的时候都学过泰勒展开,当时可能就在想“这玩意有啥用?”,那其实线性化就是一个对于复杂非线性问题的泰勒级数展开,只是这里只取其一次项进行线性近似,忽略其他高阶项。首先我们来说一下牛顿法的基本流程,牛顿法的基本思想就是根据当前的状态求出切线矩阵,并求出残余向量,进行位移增量的求解,即:这本质上就是一种“以直代曲”的思想去逐步迭代接近真实值,如下图所示:那么何为切线矩阵呢?其实这里我们已经隐含了线性化的概念,现在我们看如下的一个非线性方程这里在非线性分析中是的函数,也是的函数(如几何非线性的情况下),这个式子非常难以求解,那么就需要我们对其进行线性化,具体做法为:那这就是一个典型对于复杂在附近的泰勒展开,但我们忽略高阶项,以上公式则变为这里为了公式简便,先暂时不考虑随着变化的情况,那么进一步展开上边的式子我们可以得到而我们可以认为,可以认为是在初始位移向量对应的内力,那么通过移项,我们就得到了最初的Newton法的基本公式,所以这也就是经常所说的牛顿法和泰勒级数展开与线性化都有很强的关系,也就是牛顿法所谓2阶收敛性的原因,因为他近似只取到泰勒级数的一阶,误差是这样对于材料非线性、接触非线性、几何非线性的问题我们就通过以下公式求解其切线矩阵现在可以回想下我所常用的材料:Mises塑性2中进行的那些应力增量的求解大家一定特别理解,因为本质上内力和应力有如下的关系那么对于求解切线刚度矩阵的时候对于内力求导,在不变的情况下则可以转化为最后大家需要注意的是,这里的线性化是非线性方程求解的线性化,与应力线性化是完全不同的概念。来源:大狗子说数值模拟

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈