Abaqus无限单元的建立方法

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. 实例讲解-钢丸撞击金属板

本人在做机械喷丸的模拟，其中设置金属板边界部分为无限单元：

（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，即某个面包含了无限单元，导致不能运算。

修改好的无限单元为：

运算后

修改前 *Heading

** Job name: Job-1 Model name: Model-1 ** Generated by: Abaqus/CAE 6.10-1

*Preprint, echo=NO, model=NO, history=NO, contact=NO **

** PARTS **

*Part, name=Part-2 *Node

1, 0., 0., 1.5 2, 0., 2., 1.5 3, -2., 0., 1.5 4, 0., 2., 0. 5, 0., 0., 0. 6, -2., 0., 0. 7, 0., 2.5, 1.5 8, 0., 2.5, 0. 9, 2., 0., 1.5 10, 2., 0., 0. 11, 2.5, 0., 0. 12, 2.5, 0., 1.5 13, -2.5, 0., 1.5 14, -2.5, 0., 0. 15, 0., 0.0786545798, 1.5 16, 0., 0.172710776, 1.5 17, 0., 0.285184383, 1.5 18, 0., 0.419681847, 1.5 19, 0., 0.580515683, 1.5 20, 0., 0.772842824, 1.5 21, 0., 1.00283027, 1.5 22, 0., 1.27785218, 1.5 23, 0., 1.60672712, 1.5 24, -0.312868923, 1.97537673, 1.5 25, -0.618034005, 1.90211308, 1.5 26, -0.907980978, 1.78201306, 1.5 27, -1.17557049, 1.61803401, 1.5 28, -1.41421354, 1.41421354, 1.5 29, -1.61803401, 1.17557049, 1.5 30, -1.78201306, 0.907980978, 1.5 31, -1.90211308, 0.618034005, 1.5 32, -1.97537673, 0.312868923, 1.5

33, -1.60672712, 0., 1.5 34, -1.27785218, 0., 1.5 35, -1.00283027, 0., 1.5 36, -0.772842824, 0., 1.5 37, -0.580515683, 0., 1.5 38, -0.419681847, 0., 1.5 39, -0.285184383, 0., 1.5 40, -0.172710776, 0., 1.5 41, -0.0786545798, 0., 1.5 42, 0., 2., 0.150000006 43, 0., 2., 0.300000012 44, 0., 2., 0.449999988 45, 0., 2., 0.600000024 46, 0., 2., 0.75 47, 0., 2., 0.899999976 48, 0., 2., 1.04999995 49, 0., 2., 1.20000005 50, 0., 2., 1.35000002 51, -2., 0., 1.35000002 52, -2., 0., 1.20000005 53, -2., 0., 1.04999995 54, -2., 0., 0.899999976 55, -2., 0., 0.75 56, -2., 0., 0.600000024 57, -2., 0., 0.449999988 58, -2., 0., 0.300000012 59, -2., 0., 0.150000006 60, 0., 0., 0.150000006 61, 0., 0., 0.300000012 62, 0., 0., 0.449999988 63, 0., 0., 0.600000024 64, 0., 0., 0.75 65, 0., 0., 0.899999976 66, 0., 0., 1.04999995 67, 0., 0., 1.20000005 68, 0., 0., 1.35000002 69, 0., 1.60672712, 0. 70, 0., 1.27785218, 0. 71, 0., 1.00283027, 0. 72, 0., 0.772842824, 0. 73, 0., 0.580515683, 0. 74, 0., 0.419681847, 0. 75, 0., 0.285184383, 0. 76, 0., 0.172710776, 0. 77, 0., 0.0786545798, 0. 78, -0.0786545798, 0., 0. 79, -0.172710776, 0., 0. 80, -0.285184383, 0., 0. 81, -0.419681847, 0., 0. 82, -0.580515683, 0., 0. 83, -0.772842824, 0., 0. 84, -1.00283027, 0., 0. 85, -1.27785218, 0., 0. 86, -1.60672712, 0., 0. 87, -1.97537673, 0.312868923, 0. 88, -1.90211308, 0.618034005, 0. 89, -1.78201306, 0.907980978, 0. 90, -1.61803401, 1.17557049, 0. 91, -1.41421354, 1.41421354, 0. 92, -1.17557049, 1.61803401, 0. 93, -0.907980978, 1.78201306, 0. 94, -0.618034005, 1.90211308, 0.