基于INP文件的ABAQUS三维无限元边界快速建模方法与流程

本发明涉及到一种无限单元边界建模方法，具体涉及一种基于inp文件的abaqus三维无限元边界快速建模方法。

背景技术：

计算机的出现和迅速发展，为工程分析提供了强有力的工具，从而为使用数值模拟求解许多大型工程实际问题成为可能。有限单元法作为最常用的数值模拟手段，广泛用于岩土，结构，地震和水利等工程领域。采用有限元方法的面临的问题是如何采用“有限”的模型模拟实际“无限”区域，常用的解决办法是根据圣维南原理忽略“无限”区域边界的影响，利用截断边界取“足够大”的模型尺寸进行几何尺寸的有限元网格划分，而在人为边界上施加相应的近似约束边界条件。实际计算时，这种“足够大”的模型尺寸鉴定困难，区域较小时计算成本低，但结果精度差；区域较大时满足一定精度要求，但是计算成本较高。特别是对于波动计算问题，常常由于网格边界上波的反射和散射造成有限元计算结果的失真。

通过引入边界单元建立合理的人工边界，可以尽量减少人工截断导致的有限元边界对地基中各种波动在边界上的反射，从而极大的减少单元的数量，提高计算效率。目前常见的人工边界如黏性边界、旁轴近似边界、透射边界、黏弹性边界、无限元边界等。

作为模拟无穷远区域一种方式，无限元的出现为克服有限元计算方法的缺陷，解决无限域边界模拟的问题提供了重要途径。无限元在概念上属于有限元的延伸，其主要思想是通过几何上无限大的“有限”单元，再对其在物理上进行界定模拟无限域物理场。又由于无限元必须反映近场的边界特征或与模拟近场的有限元结合，它实际上只在一个方向趋于无限，因而又被称为半无限元。从广义有限元概念上讲无限元仍然属于有限元的范畴。总之，无限元为克服有限元在解决无界域问题时而提出，常常与常规有限元同时用来解决更复杂的无界问题，是对有限元方法的一种补充，因而它与有限元方法的“协调”与生俱在，比边界元等其它求解无界域问题的数值方法更具有优势.

大量求解无限域问题的经验表明：有限元与无限元耦合模型在求解工程实际问题方面有着广泛的实用性。特别是在内源波动问题的研究，诸如列车振动和爆破等，以及外源振动的地震动力学问题方面等，无限元在模拟和近似模拟无限域问题方面表现出明显的优越性。此外有限元和无限元耦合模型也广泛用于电磁学，热力学以及声学等领域，并取得了良好的模拟效果。大型有限元计算软件abaqus和ansys中均提高了多种无限单元类型，模拟有限元计算的无限域问题。其中abaqus以其强大的动力非线性计算能力广泛用于岩土工程领域的动力计算。

abaqus提供了一阶和二阶无限单元，包括平面应变、平面应力、轴对称和三维无限元单元，这是基于zienkiewicz等静力计算分析，以及lysmer等动力响应分析而开发的，可用于静力无限域问题求解，以及作为无限元动力人工边界求解域内局部源振动问题，即对从有限域穿过人工边界进入无限域的外行波的模拟有效，而对外源入射问题通过abaqus无限元基础上进行二次开发也成功解决了地震动力学问题。abaqus无限元动力人工边界常见的三维无限单元包括cin3d8，cin3d12r^(s)和cin3d18r(s)等，这种单元可以与标准有限单元结合，用有限元模拟近场区域，而用无限元模拟远场区域。

尽管abaqus边界无限单元为静力无限元和动力人工边界模拟提供了便利。但是由于无限单元中节点编号要保证单元的第一个面为有限元和无限元的交接面，以保证无限元单元的方向性，单元延伸方向是从近场到远场区。abaqus中不能直接定义无限单元，只能通过建模时预设建出无限元边界部分，并采用其他单元类型以区分，导出到abaqus输入文件(.inp)，然后对比模型手动修改abaqus的输入文件(.inp)的无限元边界的单元节点顺序。特别是对于三维边界无限单元，修改过程繁琐，尤其是分析任务量大时，会耗费更多的人力和时间，并且经常会出现手动操作失误而导致网格操作错误的情况。

技术实现要素：
：

本发明所要解决的技术问题是提供一种基于inp文件的abaqus三维无限元边界快速建模方法，克服abaqus现有边界无限元模型建模过程繁琐且建模效率低，考虑到abaqus无限单元特点及abaqus的输入inp文件模型定义的特点，操作简单，自动化程度高，能够在较短时间内生成无限单元网格，极大地提高了abaqus的建模效率。

为解决上述技术问题，本发明采用的技术方案是：

根据abaqus三维无限单元(3dinfiniteelement)中单元节点定义规则，利用abaqus输入文件(.inp)的单元和节点定义方法，并借助于python实现输出abaqus无限单元定义的数据行，按照inp文件的语法规则完成模型定义的输入文件，导入abaqus/cae生成边界无限单元。

本发明解决其技术问题所采用的技术方案包括以下步骤：

建立有限元模型并创建指定part边界节点集作为无限元的近场节点集，导出该abaqus输入文件(job-1.inp)；

修改inp文件利用关键字(*ncopy)定义无限单元的远场节点，并另存为新的inp文件(job-2.inp)；

abaqus/cae导入修改后inp文件(job-2.inp)，查看生成的无限元远场节点并导出新inp文件(job-3.inp)；

inp文件(job-3.inp)指定位置插入关键字行(*element)定义无限单元，并另存为job-4.inp；

python编程构建无限元边界单元定义数据行并输入到abaqus的输入文件(job-4.inp)或者指定的inp文件；

通过abaqus模型导入的方式打开上述inp文件(job-4.inp)，abaqus指定part三维无限单元边界自动生成完毕。

进一步的，建立有限元模型中无限元近场节点集的并导出abaqus输入文件(job-1.inp)具体步骤包括：

建立abaqus/cae三维有限元模型，根据研究问题在part模块下利用set工具建立指定part的边界面上的节点集作为对应边界无限元的近场节点集，即无限单元第一个面上的四个节点；

在job模块下为上述问题的模型建立分析任务job-1，并在jobmanager上通过writeinput导出该abaqus输入文件(job-1.inp)。

修改上述inp文件利用关键字(*ncopy)生成无限单元的远场节点，并另存为新的inp文件(job-2.inp)主要内容包括：

修改上述inp文件job-1.inp，根据不同的边界区域在指定part定义行内插入多行节点定义的关键字(*ncopy)；

利用关键字(*ncopy)将近场节点复制偏移生成不同区域边界无限元的远场节点；

将上述修改后的inp文件另存为新的inp文件(job-2.inp)；

abaqus/cae导入修改后inp文件job-2.inp，查看生成的无限元远场节点并导出新inp文件job-3.inp包括：

通过abaqus/cae的图形用户界面gui以file→import→model的方式导入修改后的inp文件job-2.inp，无限元的远场节点生成，但暂时无法显示；

在part模块条件下，通过view→partdisplayoptions→mesh下点选shownodelabels显示单元节点号，查看生成的不同区域无限元远场节点；

将所述导入模型在job模块内，新建job-3任务，通过jobmanager上的writeinput导出新的inp文件job-3.inp；

inp文件job-3.inp指定位置插入关键字行*element定义无限单元，并另存为job-4.inp包括：

编辑inp文件job-3.inp，在指定行插入关键字行*element进行无限单元定义；

插入*element关键字行和数据行后的inp文件保存为新的inp文件job-4.inp。

进一步的，利用关键字(*ncopy)将近场节点复制偏移生成不同区域边界无限元的远场节点主要内容包括：

节点定义的关键字行和数据行必须按照abaqus的inp文件的语法规则编辑；

通过*ncopy关键字行和数据行定义新节点。

进一步的，关键字行和数据行必须按照abaqus的inp文件的语法规则编辑主要内容是指：

1)远场节点定义关键字行(keywordline)为：

*ncopy,changenumber＝n,oldset＝near_field_nodesetname,newset＝far_field_nodesetname,shift

2)远场节点定义数据行(dataline)为：

shifting-x,shifting-y,shifting-z

上述节点定义的关键字行和数据行插入位置必须位于无限元边界所在part的定义结束行*endpart关键字之前，近场节点集合定义数据行之后。

通过*ncopy关键字行和数据行定义新节点主要是指：

通过近场节点复制偏移一定距离，定义无限元边界远场节点，新节点的节点号在原有节点号基础上增加整数值n；

实现了同一无限单元延无限元方向上近场节点和远场节点一一对应，且单元号相差n；

数据行定义了新创建的远场节点在原有近场节点基础上的坐标偏移值。

进一步的，编辑inp文件(job-3.inp)，在指定位置插入关键字(*element)行进行无限单元定义详细内容是指：

通过*element关键字行和数据行定义新无限单元；

无限单元定义的关键字行和数据行必须按照abaqus的inp文件的语法规则编辑；

上述单元定义的关键字行和数据行插入位置必须位于无限元边界所在part的定义结束行*endpart关键字之前，远场节点集合定义数据行和有限元单元定义数据行之后。

进一步的，通过*element关键字行和数据行定义边界无限单元具体内容包括：

abaqus的单元和节点的定义主要是通过其前处理器abaqus/cae，创建几何形状并划分单元网格，形成单元节点，这里的abaqus无限单元在inp文件中定义是通过手动编辑inp文件，插入*element关键字，给无限单元分配单元号，并通过指定单元节点号定义单元，同时可以给边界无限元分组。

考虑到在装配件(assembly)上分配的单元号不唯一，需要结合部件(part)区分，而在部件上分配的单元号唯一，所以边界无限元的定义只在指定的part定义；

主要通过*element关键字定义无限单元，指定单元类型为cin3d8，并建立相应边界区域上无限边界元集合；

无限单元定义的数据行是指定单元号和单元组成的节点号；

abaqus无限单元定义的数据行也可以从其他包含单元定义数据的inp文件中读取；

上述数据行单元定义的节点号排列顺序满足abaqus帮助文档中三维无限单元(cin3d8)定义的节点次序。

无限单元定义的关键字行和数据行必须按照abaqus的inp文件的语法规则编辑具体是指：

1)无限元边界无限单元定义的关键字行(keywordline)有两种方式，分别为：

①*element,type＝cin3d8,elset＝infinite_element_setname

②*element,type＝cin3d8,elset＝infinite_element_setname,input＝filename.inp

2)相应的无限单元定义数据行(dataline)可以由两种方式，分别与上述关键字行对应：

①无限单元定义的数据行直接插入对应区域无限元定义关键字所在位置的下一行，主要包括无限单元号和组成节点号，即：

infinite_element_num,nodep1,nodep2,nodep3,nodep4,nodep5,nodep6,nodep7,nodep8

②无限单元定义的数据行由filename.inp给定，该文件内容只是包含预先完成输入的当前边界区域无限单元定义的数据行—包括无限单元号和组成节点号，即：

infinite_element_num,nodep1,nodep2,nodep3,nodep4,nodep5,nodep6,nodep7,nodep8。

进一步的，无限单元定义的数据行infinite_element_num,nodep1,nodep2,nodep3,nodep4,nodep5,nodep6,nodep7,nodep8相关约定包括：

根据权利要求5和6所述的远场四个节点的节点号比对应沿无限方向单元边上近场节点(有限元边界上节点)的节点号值大n，即：

根据无限单元定义规则，无限单元第一面内近场四个节点nodep1,nodep2,nodep3和nodep4，要满足从无限单元远场向近场观测时近场四个节点为逆时针方向的顺序排列；

无限单元定义指定的单元号是由当前part上有限单元个数，已定义边界区域无限单元个数以及区域内无限单元定义的次序决定。

进一步的，利用python编程构建无限元边界单元定义数据行并输入到abaqus的输入文件(job-4.inp)或者指定的inp文件具体包括以下步骤：

1)利用python编程读取abaqus的输入文件(job-4.inp)中的节点号及相应坐标，单元号及组成节点号，并存储在指定的列表中；

2)利用python编程读取abaqus的输入文件(job-4.inp)中无限单元定义的关键字行号并存储在指定的变量中，或者读取无限单元定义数据行写入的input参数—abaqus读取的inp文件filename.inp；

3)根据节点坐标特点遍历上述有限元模型的节点列表，根据各边界节点坐标值特点查找不同边界上节点，并添加到指定边界节点列表；

4)通过遍历上述有限元模型的单元列表，依据单元组成节点坐标与边界上节点坐标特点的关系，查找边界单元，并添加到指定的边界单元列表；

5)遍历指定边界节点列表，并取出其中任意四个节点a,b,c和d，判断这四个节点是否为边界无限单元近场第一个面上的四个节点；

6)若a,b,c和d四个节点是边界无限单元近场第一个面上的四个节点，则根据向量的数量积求向量夹角关系确定四个节点的相对位置关系；

7)利用向量间的向量积确定四个节点的逆时针方向的排列顺序为q1,q2,q3,q4；

8)无限单元序号由有限模型中指定part单元个数a，已定义边界无限单元个数为b以及该区域内无限单元定义次序c决定，则当前无限单元号为q＝a+b+c；

9)给出指定边界单元无限元定义的数据行输出：

q,q1,q2,q3,q4,q1+n,q2+n,q3+n,q4+n

其中：q为无限元单元号，q1,q2,q3,q4为近场四个节点，q1+n,q2+n,q3+n,q4+n为远场四个节点；

10)两种无限单元定义数据行输出方式，前者主要是针对无限单元定义采用关键字行①和数据行①，即单元定义数据行直接插入在job-4.inp中关键字定义位置的下一行并保存；

后者主要是针对无限单元定义采用关键字行②和数据行②，将单元定义的数据行写入无限单元定义行中的input参数filename.inp文件中，abaqus在导入模型时会自动读取filename.inp文件中的数据行，生成无限单元。

进一步的，遍历指定边界节点列表，并取出其中任意四个节点，判断该节点是否为无限元模型近场第一个面上的四个节点具体步骤为：

利用python的成员测试功能，若边界节点列表上取出的四个节点为边界单元组成节点的其中四个节点，则认为这四个节点为边界无限单元第一个面上的四个近场节点。

若a,b,c和d四个节点是边界无限单元近场第一个面上的四个节点，根据向量数量积求两个向量夹角确定四个节点的相对位置关系，具体步骤包括：

以点a作为四边形的第一个点，利用向量数量积分别计算∠bac,∠bad和∠cad的大小并求和:

令sum＝∠bac+∠bad+∠cad：

1)若sum＝2∠bac，则点a和点d，点b和点c互为对角节点，则节点逆时针方向顺序为a,b,d,c或者a,c,d,b；

2)若sum＝2∠bad，则点a和点c，点b和点d互为对角节点，则节点逆时针方向顺序为a,b,c,d或者a,d,c,b；

3)若sum＝2∠cad，则点a和点b，点c和点d互为对角节点，则节点逆时针方向顺序为a,d,b,c或者a,c,b,d。

利用向量间的向量积确定四个节点的逆时针方向的排列顺序具体步骤包括：

确定边界无限单元第一面内的四个近场节点p1,p2,p3,p4相对位置关系后，通过四个节点中的任意三个相邻节点首尾相接形成两个向量，假设p1与p3,p2与p4互为角点，以p1点为起始点，相邻节点必为p2和p4，则上述四个节点的逆时针顺序为p1,p2,p3,p4或者p1,p4,p3,p2，计算向量与的向量积

若与当前边界面上点p1的法向量夹角为锐角时，即：

则点p1,p2,p3为逆时针方向排列，四个节点的逆时针排列顺序为p1,p2,p3,p4；

若与当前边界面上点p1的法向量夹角为钝角时，即：

则点p1,p2,p3为顺时针方向排列，四个节点的逆时针排列顺序为p1,p4,p3,p2。

进一步的，通过abaqus模型导入的方式打开上述inp文件(job-4.inp)，abaqus指定part三维无限单元边界自动生成完毕具体步骤包括：

打开abaqus/cae，通过file菜单下的模型导入方式(import→model)，打开目标文件夹下的inp文件(job-4.inp),完成边界无限单元的建模过程，可在指定part下查看无限元建模效果。

与现有技术相比，本发明的有益效果是：本方法主要是根据abaqus三维无限单元中单元节点定义规则，利用abaqus输入文件(.inp)的单元和节点定义方法，并辅以一定的python编程完成无限单元定义数据行的自动输出，按照inp文件的语法规则完成模型定义的输入文件，导入abaqus/cae生成边界无限单元。本发明实施方案操作简单，避免了传统方法中的手动修改inp文件的繁琐过程，自动化程度高，能够在较短时间内生成无限单元网格，极大地提高了abaqus的建模效率。

附图说明

图1是本发明技术方案流程图。

图2是abaqus三维有限元模型。

图3是利用*ncopy生成的无限元远场节点。

图4是生成无限单元效果图。

图5是模型底面边界区域无限单元。

图6是模型左侧边界区域无限单元。

图7是模型右侧边界区域无限单元。

图8是模型后端边界区域无限单元。

图9是模型前端边界区域无限单元。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细的说明。

如图1所示，本发明实施例中基于inp文件的abaqus三维无限元边界快速建模方法包括：

基于三维单元网格划分完成的三维模型，并且所述三维模型的后面中心点位于坐标系的原点上，坐标系的x轴，y轴和z轴为其对称轴；

101.建立有限元模型并创建指定part边界节点集作为无限元的近场节点集，导出该abaqus输入文件(job-1.inp)；

建立abaqus/cae三维有限元模型如图2所示，根据研究问题在part模块下利用set工具建立指定part的边界面上的节点集作为对应边界无限元的近场节点(无限单元第一个面上的四个节点)集，不同边界区域的节点集情况如表1所示；在job模块下为上述问题的模型建立分析任务job-1，并在jobmanager上通过writeinput导出该abaqus输入文件(job-1.inp)。

102.修改inp文件利用关键字(*ncopy)定义无限单元的远场节点，并另存为新的inp文件(job-2.inp)；

修改上述inp文件job-1.inp，根据不同的边界区域在指定part定义行内插入多行节点定义的关键字(*ncopy)，命令如下所示

*ncopy,changenumber＝2000,oldset＝node-y1,newset＝node-y2,shift

0.0,-40.0,0.0

*ncopy,changenumber＝4000,oldset＝node-x1,newset＝node-x3,shift

-40.0,0.0,0.0

*ncopy,changenumber＝6000,oldset＝node-x2,newset＝node-x4,shift

40.0,0.0,0.0

*ncopy,changenumber＝8000,oldset＝node-z1,newset＝node-z3,shift

0.0,0.0,40.0

*ncopy,changenumber＝10000,oldset＝node-z2,newset＝node-z4,shift

0.0,0.0,-40.0

*ncopy,changenumber＝12000,oldset＝node-x1z1,newset＝node-x3z3,shift

-40.0,0.0,40.0

*ncopy,changenumber＝14000,oldset＝node-x1z2,newset＝node-x3z4,shift

-40.0,0.0,-40.0

*ncopy,changenumber＝16000,oldset＝node-x2z1,newset＝node-x4z3,shift

40.0,0.0,40.0

*ncopy,changenumber＝18000,oldset＝node-x2z2,newset＝node-x4z4,shift

40.0,0.0,-40.0

利用关键字(*ncopy)将近场节点复制偏移指定距离生成不同区域边界无限元的远场节点集如表2，将上述修改后的inp文件另存为新的inp文件(job-2.inp)。

通过近场节点复制偏移一定距离，定义无限元边界远场节点，新节点的节点号在原有节点号基础上增加整数值n，实现同一无限单元延无限元方向上近场节点和远场节点一一对应，且单元号相差n，n值的选取要结合模型中节点总数选取，确保原有限元节点号与新生成节点号不重复。数据行定义了新创建的远场节点在原有近场节点基础上的坐标偏移值。

表1不同边界区域无限元近场和远场节点集

103.abaqus/cae导入修改后inp文件(job-2.inp)，查看生成的无限元远场节点并导出新inp文件(job-3.inp)；

通过abaqus/cae的图形用户界面gui以file→import→model的方式导入修改后的inp文件(job-2.inp)，无限元的远场节点生成，在part模块条件下，通过view→partdisplayoptions→mesh下点选shownodelabels显示单元节点号，可以查看生成的不同区域无限元远场节点，如图3所示，将上述导入模型在job模块内，新建job-3任务，通过jobmanager上的writeinput导出新的inp文件(job-3.inp)。

104.inp文件(job-3.inp)指定位置插入关键字行(*element)定义无限单元，并另存为job-4.inp；

按照abaqus的inp文件的语法规则编辑inp文件(job-3.inp)，通过*element关键字行和数据行定义新单元，且上述单元定义的关键字行和数据行插入位置必须位于无限元边界所在part的定义结束行*endpart关键字之前，远场节点集合定义数据行和有限元单元定义数据行之后。

abaqus的inp文件主要通过*element关键字定义无限单元，指定单元类型为cin3d8，并建立相应边界区域上无限边界元集合；其数据行是指定单元号和单元组成的节点号，对于三维无限单元abaqus帮助文档给出了其定义规则，该数据行也可以从其他包含单元定义数据的inp文件中读取。

上述数据行单元定义的节点号排列顺序满足abaqus帮助文档中三维无限单元(cin3d8)定义的节点次序。

abaqus的inp文件中利用*element进行无限单元定义两种方式：

1)无限单元定义的数据行直接插入对应区域无限元定义关键字所在位置的下一行，命令如下所示；

*element,type＝cin3d8,elset＝bottomy

*element,type＝cin3d8,elset＝leftx

*element,type＝cin3d8,elset＝rightx

*element,type＝cin3d8,elset＝frontz

*element,type＝cin3d8,elset＝backz

*element,type＝cin3d8,elset＝leftfrontcorner

*element,type＝cin3d8,elset＝leftbackcorner

*element,type＝cin3d8,elset＝rightfrontcorner

*element,type＝cin3d8,elset＝rightbackcorner

*element,type＝cin3d8,elset＝frontleftcorner

*element,type＝cin3d8,elset＝frontrightcorner

*element,type＝cin3d8,elset＝backleftcorner

*element,type＝cin3d8,elset＝backrightcorner

2)无限单元定义的数据行由外部文件filename.inp给定，该文件内容只是包含预先完成输入的当前边界区域无限单元定义的数据行—包括无限单元号和组成节点号，命令如下所示。

*element,type＝cin3d8,elset＝bottomy,input＝bottomy_element1.inp

*element,type＝cin3d8,elset＝leftx,input＝leftx_element1.inp

*element,type＝cin3d8,elset＝rightx,input＝rightx_element1.inp

*element,type＝cin3d8,elset＝frontz,input＝frontz_element1.inp

*element,type＝cin3d8,elset＝backz,input＝backz_element1.inp

*element,type＝cin3d8,elset＝leftfrontcorner,input＝leftfrontcorner_element.inp

*element,type＝cin3d8,elset＝leftbackcorner,input＝leftbackcorner_element.inp

*element,type＝cin3d8,elset＝rightfrontcorner,input＝rightfrontcorner_element.inp

*element,type＝cin3d8,elset＝rightbackcorner,input＝rightbackcorner_element.inp

*element,type＝cin3d8,elset＝frontleftcorner,input＝frontleftcorner_element.inp

*element,type＝cin3d8,elset＝frontrightcorner,input＝frontrightcorner_element.inp

*element,type＝cin3d8,elset＝backleftcorner,input＝backleftcorner_element.inp

*element,type＝cin3d8,elset＝backrightcorner,input＝backrightcorner_element.inp

105.python编程构建无限元边界单元定义数据行并输入到abaqus的输入文件(job-4.inp)或者指定的inp文件，其具体过程为：

1)利用python编程读取abaqus的输入文件(job-4.inp)中的节点号及相应坐标，单元号及组成节点号，并存储在指定的列表中；

2)利用python编程读取abaqus的输入文件(job-4.inp)中无限单元定义的关键字行号并存储在指定的变量中，或者读取无限单元定义数据行写入的input参数filename.inp；

3)根据节点坐标特点遍历上述有限元模型的节点列表，根据各边界节点坐标值特点查找不同边界上节点，并添加到指定边界节点列表；

4)通过遍历上述有限元模型的单元列表，依据单元组成节点坐标与边界上节点坐标特点的关系，查找边界单元，并添加到指定的边界单元列表；

5)遍历指定边界节点列表，并取出其中任意四个节点a,b,c和d，判断这四个节点是否为边界无限单元近场第一个面上的四个节点；

6)若a,b,c和d四个节点是边界无限单元近场第一个面上的四个节点，则根据向量的数量积求向量夹角关系确定四个节点的相对位置关系；

7)利用向量间的向量积确定四个节点的逆时针方向的排列顺序为q1,q2,q3,q4；

8)无限单元序号由有限模型中指定part单元个数a，已定义边界无限单元个数为b以及该区域内无限单元定义次序c决定，则当前无限单元号为q＝a+b+c；

9)给出指定边界单元无限元定义的数据行输出：

q,q1,q2,q3,q4,q1+n,q2+n,q3+n,q4+n

其中：q为无限元单元号，q1,q2,q3,q4为近场四个节点，q1+n,q2+n,q3+n,q4+n为远场四个节点；

10)两种无限单元定义数据行输出方式，前者主要是针对无限单元定义采用关键字行①和数据行①，即单元定义数据行直接插入在job-4.inp中关键字定义位置的下一行并保存；

后者主要是针对无限单元定义采用关键字行②和数据行②，将单元定义的数据行写入无限单元定义行中的input参数filename.inp文件中，abaqus在导入模型时会自动读取filename.inp文件中的数据行生成无限单元。

106.通过abaqus模型导入的方式打开上述inp文件(job-4.inp)，abaqus指定part三维无限单元边界自动生成完毕。

打开abaqus/cae，通过file菜单下的模型导入方式(import→model)，打开目标文件夹下的inp文件(job-4.inp),完成边界无限单元的建模过程，可在指定part下查看无限元建模效果，如图4所示。图5-图9为不同边界区域的生成的无限单元图。