Abaqus无限单元的建立方法
1. 什么是无限单元
无限单元是Abaqus单元库中的一种单元,它的单元形式如下图所示
2. 无限单元的作用
使用无限单元作为反射边界,将无反射,防止在边界上产生的应力波反射,重新进入模型,从而导致结果不正确。
3. 怎么建立无限单元
(1)先建立CAE模型,通过partition将需要设置为无限单元的部分分割,在mesh模块中设置该部分单元类型,这样在inp文件中需要设置为无限单元的部分就会集中在一起,方便修改其节点的编号顺序。
(2)在inp文件里对无限单元进行单元属性及编号的改变,需要注意无限单元的方向,在二维无限单元中,前两个节点所组成直线中点A与后两个节点所组成直线中点B,无限单元的方向就是A指向B的方向,如图1,二维无限元的方向是朝下;在三维无限单元中,前四个节点所组成平面的中心点C与后四个节点所组成平面的中心点D,无限单元的方向就是C指向D的方向,如图1,三维无限单元的方向朝右。
(3)将修改后的inp导入,建立job提交就可以了。
4. 实例CAE-钢丸撞击金属板
本人在做机械喷丸的模拟,其中设置金属板边界部分为无限单元:
(1)建立CAE,将金属板的边界partition切割,全部设置为C3D8R单元类型,修改边界部分为C3D8I,建立job-write input 生成inp文件。
(2)修改inp文件的单元属性及编号
将C3D8修改为CIN3D8,
下一步修改无限单元节点编号的排列顺序在inp的data lines中第一列是单元编号,后面的是节点编号,在本文中,无限单元的方向是从金属板的中心向外。View-Assembly Display Option-Mesh-show node labels(显示节点编号),show element labels(显示单元编号),在mesh模块中可以看到各个单元的编号和节点的编号,任意找一个单元查看,750, 710, 114, 113, 621, 158, 12, 9, 141,
在网格显示中,750单元的节点编号顺序如上图所示,在此,无限单元的方向是朝右,inp修改的内容如下图。
而inp中的编号排列是750, 710, 114, 113, 621, 158, 12, 9, 141,所以需要进行修改,将12,158,710,114(逆时针排列)四个节点编号放在后面,而且12与9在一条直线上,改为750, 621, 113, 9, 141, 710, 114, 12, 158。对于大范围节点编号顺序的修改,建议采用UltraEdit软件,列-列模式进行批量修改。
(3)导入inp文件,建立job提交就可以了。
注意事项:在interaction和load模块中,定义面的时候不要将无限单元定义在里面,因为在修改单元编号时,也把面重新定义了,在提交inp文件时出现错误abaqus: Error in job Job-6: THE FACE S6 IS NOT VALID FOR THE UNDERLYING ELEMENT 56 INSTANCE PART-5-1,即某个面包含了无限单元,导致不能运算。
修改好的无限单元为:运算后