Abaqus无限单元的建立方法.docx

Abaqus无限单元的建立方法
什么是无限单元无限单元是Abaqus单元库中的一种单元，它的单元形式如下图所示
无限单元的作用
使用无限单元作为反射边界，将无反射，防止在边界上产生的应力波反射，重新进入模型，从而导致结果不正确。
怎么建立无限单元（1）先建立CAE模型，通过partition将需要设置为无限单元的部分分割，在mesh模块中设置该部分单元类型，这样在inp文件中需要设置为无限单元的部分就会集中在一起，方便修改其节点的编号顺序。
（2）在inp文件里对无限单元进行单元属性及编号的改变，需要注意无限单元的方向，在二维无限单元中，前两个节点所组成直线中点A与后两个节点所组成直线中点B，无限单元的方向就是A指向B的方向，如图1，二维无限元的方向是朝下；在三维无限单元中，前四个节点所组成平面的中心点C与后四个节点所组成平面的中心点D，无限单元的方向就是C指向D的方向，如图1，三维无限单元的方向朝右。
（3）将修改后的inp导入，建立job提交就可以了。
实例讲解-钢丸撞击金属板本人在做机械喷丸的模拟，其中设置金属板边界部分为无限单元：（1）建立CAE，将金属板的边界partition切割，全部设置为C3D8R单元类型，修改边界部分为C3D8I，建立job-writeinput生成inp文件。
（2）修改inp文件的单元属性及编号将C3D8修改为CIN3D8，下一步修改无限单元节点编号的排列顺序在inp的datalines中第一列是单元编号，后面的是节点编号，在本文中，无限单元的方向是从金属板的中心向外。View-AssemblyDisplayOption-Mesh-shownodelabels（显示节点编号），showelementlabels（显示单元编号），在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:ErrorinjobJob-6:THEFACES6ISNOTVALIDFORTHEUNDERLYINGELEMENT56INSTANCEPART-5-1，即某个面包含了无限单元，导致不能运算。
修改好的无限单元为：运算后
修改前*Heading
**Jobname:Job-1Modelname:Model-1
**Generatedby:Abaqus/CAE6.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.9