首页
发现
课程
培训
文章
案例
问答
需求
服务
行家
赛事
热门搜索
发布
消息
注册
|
登录
首页
/
文章
/
详情
ABAQUS显式动力学入门案例-丢硬币/掷骰 子
团长
1年前
浏览5523
关注
大家好,我是团长。
从【CAE知识地图】公 众号发现一个很不错的学习案例,
用两个刚体显式动力学案例入门Abaqus/Explicit
,转载一下,有兴趣的可以看看,也可以下载模型自己尝试哦(见文末)。
以下原文(略有改动):
最近科研项目遇到一点小小的挫折,涉及接触的算例怎么修改边界条件都难以收敛。我猜想如果换用显式动力学可能可以让问题更快实现收敛。之前虽然玩过一两次,但一直没有很认真的正式学习显式动力学。
今天这篇文章,就用USim邓怡超工程师几年前分享的两个inp文件,
学习一下Abaqus/Explicit显式动力学分析
。
既然用了别人的模型,那就在开头显著位置写一下致谢吧。感谢USim,在我还在犹豫没想好是在仿真领域一条道走到黑还是毕业转行当码农的时刻,是USim邓工分享的一波Abaqus算例动画让我燃起了对有限元仿真的极大热情(emmm这么说来也不知是福是祸,但我确实发自内心的喜欢力学和有限元仿真)。
01、显式动力学
先翻译一些不那么有趣的东西。
如果想直接看丢硬币和投骰 子的小伙伴可以往下滑(P.S.还是按惯例翻译自帮助文档。写到这里发现我的彩云小译半年会员居然已经到期了,果断再续一年~)
什么是显式动力学分析?
它和隐式分析的区别在何处?
使用直接积分法的隐式动力学在Abaqus中是由Abaqus/Standard求解器提供的,而Abaqus/Explicit关于时间部分使用了中心差分算子。在隐式动力学分析中,必须对积分算子矩阵求逆,对每个时间增量步都需要求解一组非线性的平衡方程。而在显式动力学中,节点的位移和速度是根据时间增量步开始时刻已知的量计算出来的,不需要组装成整体质量阵和刚度阵并求逆。这样显式分析中每一个增量步的计算成本相对隐式来说更低。然而,在显式动力学分析中,稳定时间增量步的大小是有限的,因为它使用的中心差分算子是条件稳定的,而Abaqus/Standard中使用的隐式算子则是无条件稳定的。因此,在Abaqus/Standard中的大多数分析里,对时间增量步都没有限制(能限制时间增量的只有求解精度)。
中心差分法的稳定极限与应力波在模型中最小单元中传播的时间有关。因此,如果模型的最小网格尺寸较小,或者材料中应力波速度较高,那么显式分析的稳定时间增量步就非常短。因此,显式分析适合用于建模总动态响应时长较短的分析,例如波的传播等问题。同时,显式分析的许多优点也适用于准静态过程,在此情况下,需要使用质量缩放功能来降低材料中应力波的传播速度。
Abaqus/Explicit有如下特点:
a、分析成本随问题规模线性增加。而隐式算法计算成本随问题规模增加的速度比线性更快。因此Abaqus/Explicit对于规模非常大的问题是很有吸引力的。
b、显式积分算法在求解极不连续的短期事件或过程时往往比隐式算法更有效。
c、Abaqus/Explicit在求解与应力波传播有关的问题时比Abaqus/Standard更有效。
这个显式时间积分,公式也很简洁明晰。下一步的速度通过前一步的速度和加速度求得,而下一时间步的位移通过速度乘以时间增量得到。
而这个u\ddot_i^N,这个加速度又是用牛顿第二定律计算得到的。由于使用了集中质量矩阵,所以可以把整个矩阵解耦合,不需要组装总体刚度阵即可将时间步向下推进。集中质量阵把所有的质量都放在对角线元素上,所以非常容易求逆,而且计算复杂度只会随自由度线性增加。显式过程不需要在同一步内进行迭代,也不需要计算切线刚度阵。
因为显式分析无需迭代也无需切线刚度阵,这也同样简化了接触问题的处理。因此,对于接触问题占主导的分析,显式动力学分析也具有优势,例如金属板成型分析。
02、丢个硬币
这个案例如前所述,来自USim公 众号分享的inp文件。模型本身倒是非常简单。硬币是体单元,下面的地面分别为体单元和面单元。
所有进行显式动力学分析的材料都需要指定密度。因为对于动力学分析来说,不管是刚体还是柔性体,都有惯性效应。
分析步使用动力,显式分析步。
这里如果使用自动的话,会自动按照材料中应力波传播的速率计算稳定时间增量。本案例作为显式动力学的入门,全部零件设置为刚体,因此由用户指定时间增量大小。
在Abaqus里,指定材料为刚体有多种方式。可以在部件模块处就将材料定义为刚体,也可以先定义成柔性体,然后在相互作用这里使用刚体约束将材料定义为刚体。
而显式动力学中定义接触的方式也比隐式分析要简单。在显式分析中,只需要定义一个全局的通用接触,general_contact即可。软件会自动在所有单元的表面都覆盖上一层接触单元。
载荷模块,底部固定,全局施加重力载荷,然后在硬币的刚体约束参考点处,施加一个角速度的预定义场。
可以看到,虽然用了30000个增量步,但总计算时长只有不到1分钟。显式动力学每个时间增量步的计算开销确实非常低。
结果——
夭寿啦w(゚Д゚)w!硬币居然摔下去了
掉落无尽深渊……
我们来把转速改小一点再试试。顺便修改一下场输出和历程输出的时间间隔。
嗯,作为一个动画来说,这回帧数够用了,而且也成功地没有掉下去。学会Abaqus显式,真的可以把它当大型通用动画软件来玩✧(≖ ◡ ≖✿)
03、掷个
骰 子
刚才那个硬币模型建立的不够细致,没分出正反面。下面这个
骰 子
模型就更进一步,把点数一起画了出来。
和丢硬币模型一样,全部零件都约束为刚体,在整个模型上添加通用接触。载荷只有一个重力,让
骰 子
自由下落。只是它们的初始下落角度比较随意。
这个点数是——小蓝2分,小紫5分。那如果我们其他参数什么都不改,只把时间步长从1e-4改为2e-5会发生什么事呢?
小蓝变成了4分,小紫则变成了2分!
……再把求解从6核修改成2核试试?
结果又不一样了。
这完全是个玄学啊。。。。
可能因为整个模型中都施加了刚体约束,而且时间增量步不同,求解使用的核心数不同,确实会对数值计算的误差产生影响。
——看来,即使是今天的计算机,想要模拟一只拉普拉斯妖预测
骰 子
掷出来的点数,也还是很难的事情呀。
04、后记
学会了显式动力学,挺多生活中有趣的物理现象就可以用有限元来模拟了。可以说要想把Abaqus玩成动画软件,显式动力学分析是必经之路。其实本文中的几个分析还可以更进一步,用Altair的工具给odb文件做个转换,用HyperView读取以后,还可以做更真实的物理渲染。(就更像动画了。。。)
本文完。
下载方式:
见附件
来源:CAEer
登录后免费查看全文
立即登录
附件
免费
ABAQUS显式动力学入门案例-丢硬币/掷骰 子.txt
下载
HyperView
Abaqus
显式动力学
碰撞
非线性
通用
汽车
ANSA
材料
渲染
Altair
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-06-29
最近编辑:1年前
团长
硕士
|
白车身结构设...
专注ANSA使用技巧-微信公众号『C...
关注
获赞 115
粉丝 970
文章 141
课程 6
点赞
收藏
0/200
清空
提交
还没有评论
课程
培训
服务
行家
准静态拉伸试验与质量缩放的应用
Abaqus材料库讲解【附插件和材料参考手册】
航空橡胶轮胎材料及制品疲劳失效的热力耦合仿真分析
分离式霍普金森压杆(SHPB)数值模拟
相关推荐
材料非线性问题matlab有限元编程求解
CoSim 系列软件
ABAQUS UMAT 材料子程序开发59讲:经典材料本构UMAT编程实例详解
从零开始学Abaqus到接触非线性分析高阶应用(赠送知识圈答疑和订阅用户交流群)
ABAQUS-复合材料工程应用案例二-玻璃纤维复合材料泡沫夹层板弯曲变形损伤失效模拟
最新文章
笃行致远 智启新程 | 热烈祝贺世冠科技深圳子公司乔迁新址
车企裁员风波席卷全球!本科生保研升学高薪就业如何“先人一步”
经典课程指导:《数字处理系统课程设计与实践》指导(1)无线信道及其估计
探秘LS-DYNA ALE法流固耦合-从原理到实战案例分享会
仿真人才库:苏州芯慧联半导体科技力学仿真工程师高新内推
热门文章
几种常见的热仿真软件
ABAQUS中Cohesive粘聚力模型的2种定义方式(附案例操作步骤)
Abaqus分析常见问题及解决方法(2):零主元和过约束
Abaqus分析常见问题及解决方法(3):负特征值(Negative Eigenvalue)
仿真笔记——ANSYS APDL命令汇总(收藏备用)
其他人都在看
STAR CCM 案例|电池包散热
工业互联网中数字孪生系统的机理+数据融合建模方法
一种基于多通道信号Transformer(MSiT)的数控机床主轴电机故障诊断方法
2024R1来了,来看看新功能
勃发奔涌!HELLO 2024 邀你发现“仿真好内容”
VIP会员
学习计划
福利任务
下载APP
联系我们
微信客服
联系客服
人工服务时间为周一至周五的9:30-19:30
非工作时间请在微信客服留言
客服热线:
4000-969-010
邮箱:
service@fangzhenxiu.com
地址:
北京市朝阳区莱锦创意园CN08座
帮助与反馈
返回顶部