从abaqus的inp文件编写节点坐标系的例子
abaqus中的局部坐标系2013
abaqus中的局部坐标系在abaqus中，在不同的模块中使用不同的局部坐标系类型，主要有下面几种局部坐标系：1、定义节点的局部坐标系2、定义节点自由度的局部坐标系3、单元局部坐标系1.定义节点的局部坐标系（*SYSTEM）有时在一个局部坐标系中定义节点坐标是方便的，然后再转换这些坐标值到总体坐标系中。
用户可以定义一个局部节点坐标系，Abaqus将转换和旋转局部坐标值（）到总体坐标系。
在input完成后立刻完成该转换，并且该转换将施加到在定义了节点坐标系后的所有输入的节点坐标上。
【1】关键字*SYSTEM关键字*SYSTEM使用规则：先定义坐标系，再使用。
先用*SYSTEM 定义局部坐标系，然后在node definition（如*NODE关键字）里面选择定义的坐标系。
返回到整体的坐标系使用*SYSTEM+ 空白的数据行。
*SYSTEM第一行1、, global X-coordinate of the origin of the local coordinate system (point a in Figure 18.58–1).2、, global Y-coordinate of the origin of the local coordinate system.3、, global Z-coordinate of the origin of the local coordinate system.The following entries are not needed for a pure translation:4、, global X-coordinate of a point on the -axis of the local coordinate system (point b in Figure 18.58–1).5、, global Y-coordinate of a point on the -axis of the local coordinate system.6、, global Z-coordinate of a point on the -axis of the local coordinate system.第二行(optional; if not provided, the Z-axis direction remains unchanged, and the -axis is projected onto the plane):1、, global X-coordinate of a point in the plane of the localcoordinate system, on the side of the positive -axis (for example, point c in Figure 18.58–1).2、, global Y-coordinate of a point in the plane of the local coordinate system, on the side of the positive -axis.3、, global Z-coordinate of a point in the plane of the local coordinate system, on the side of the positive -axis.如果*SYSTEM中只给出一个点（原点），Abaqus假设用户只需要平移。
abaqus 提取节点坐标参考资料
1、umat中如何得到单元节点的信息大家好：最近在学习umat（standard），因为程序中要用到单元节点的信息，特别是当前增量步的单元各个节点上的节点位移值和节点的整体坐标值。
据我所知，UMAT中由abqus主程序传入的STESS，COORDS等值分别是当前积分点上的信息。
请大家指点，如何得到当前增量步中积分点所在单元各个节点上信息，特别是节点位移与节点整体坐标。
有没有类似STESS，COORDS的变量可以直接传入UMAT中。
谢谢！望大家出手指点确实UMAT是在每个material point也就是高斯积分点调用，对于每一增量步，第一次迭代调用两次，以后每步迭代调用一次。
所有的操作都是针对material point来说的。
COORDS确实只是是当前积分点上的信息。
然而，我们可以应用utility subroutine GETPARTINFO来得到此时传入UMAT的节点和单元号，有了节点和单元号当然就知道他们的信息了。
2、如何输出后处理中部分节点的坐标？谢谢如果只需要一个节点的坐标的话可以通过Query来得到，但是如果要得到大量的节点坐标就不好弄了。
如何才能得到后处理中大量节点的坐标，如下图，我要得到节点150-285的坐标该如何弄啊？多谢各位！！！Re:如何输出后处理中部分节点的坐标？谢谢建立一个set然后后处理的时候用python编个程序把节点号码和坐标都读出来Re:如何输出后处理中部分节点的坐标？谢谢设置节点SET，使用如下命令：*NODE PRINT, NSET=？？, FREQUENCY=？？COORD到*.dat文件中查看。
Re:如何输出后处理中部分节点的坐标？谢谢*NODE PRINT, NSET=？？, FREQUENCY=？？COORD谢谢，上面那些该放在inp文件的什么位置阿？Re:如何输出后处理中部分节点的坐标？谢谢nset设置放到相应的part中，*node print放到完成你要的step后部。
abaqus中的局部坐标系2013
abaqus中的局部坐标系在abaqus中，在不同的模块中使用不同的局部坐标系类型，主要有下面几种局部坐标系：1、定义节点的局部坐标系2、定义节点自由度的局部坐标系3、单元局部坐标系1.定义节点的局部坐标系（*SYSTEM）有时在一个局部坐标系中定义节点坐标是方便的，然后再转换这些坐标值到总体坐标系中。
用户可以定义一个局部节点坐标系，Abaqus将转换和旋转局部坐标值（）到总体坐标系。
在input完成后立刻完成该转换，并且该转换将施加到在定义了节点坐标系后的所有输入的节点坐标上。
【1】关键字*SYSTEM关键字*SYSTEM使用规则：先定义坐标系，再使用。
先用*SYSTEM 定义局部坐标系，然后在node definition（如*NODE关键字）里面选择定义的坐标系。
返回到整体的坐标系使用*SYSTEM+ 空白的数据行。
*SYSTEM第一行1、, global X-coordinate of the origin of the local coordinate system (point a in Figure 18.58–1).2、, global Y-coordinate of the origin of the local coordinate system.3、, global Z-coordinate of the origin of the local coordinate system.The following entries are not needed for a pure translation:4、, global X-coordinate of a point on the -axis of the local coordinate system (point b in Figure 18.58–1).5、, global Y-coordinate of a point on the -axis of the local coordinate system.6、, global Z-coordinate of a point on the -axis of the local coordinate system.第二行(optional; if not provided, the Z-axis direction remains unchanged, and the -axis is projected onto the plane):1、, global X-coordinate of a point in the plane of the localcoordinate system, on the side of the positive -axis (for example, point c in Figure 18.58–1).2、, global Y-coordinate of a point in the plane of the local coordinate system, on the side of the positive -axis.3、, global Z-coordinate of a point in the plane of the local coordinate system, on the side of the positive -axis.如果*SYSTEM中只给出一个点（原点），Abaqus假设用户只需要平移。
abaqus_inp文件精讲
如何写input文件一、输入文件的组成和结构：1．一个输入文件由模型数据和历史数据两部分组成.模型数据的作用:定义一个有限元模型.包括单元,节点,单元性质,定义材料等等有关说明模型自身的数据.模型数据可被组织到零件中(零件可以被组装成一个有意义的模型).历史数据的定义是模型发生了什么----事情的进展,模型响应的荷载,历史被分成一系列的时步层序.每一步就是一个响应(静态加载,动态响应等)，时步的定义包括过程类型(比如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出要求.At a minimum the model consists of the following information: geometry, element section properties, material data, loads and boundary conditions, analysis type, and output requests.2．ABAQUS输入文件的结构形式。
1) 必须有一个*HEADING开头。
2)接下来就是模型数据部分，定义节点，单元，材料，初始条件等。
模型数据的层次为：部件，组装，模型。
必须的模型数据：（1）几何数据：模型的几何形状是用单元和节点来定义的，结构性单元的截面是必须定义的。
比如梁单元。
特殊的特征也可以用特殊的单元来定义，比如弹簧单元，阻尼器，点式群体等。
（2）材料的定义：材料必须定义比如使用的是钢啊，岩石，土啊等材料。
可选的模型数据：（1）零件和组合：一个模型可以用几个零件来定义有可以把几个零件组合成一个集来定义。
（2）初始条件：比如初始应力，温度，或者速度等（3）边界条件：（4）运动约束（5）相互作用（6）振幅定义（7）输出控制（8）环境特性（9）用户子程序（10）分析附属部分3)接下来就是历史数据：定义分析的类型，荷载，输出要求等。
分析的目的就是预测模型对某些外部荷载或者某些初始条件的反映。
abaqus中定义管道局部坐标系、切线载荷
1.CAE操作对应inp文件：*Orientation, name="Datum csys-2", system=CYLINDRICAL-20., -2.5, 20., -20., -2.5, 21. ** Name: Load-1 Type: Surface traction*Dsload, orientation="Datum csys-2", constant resultant=YESSurf-1, TRSHR, 10., -0.707107, 0.707107, 0.2.操作原理：using a local coordinate system to define shear directionsIt is sometimes convenient to give shear and general traction directions with respect to a local coordinate system. The following two examples illustrate the specification of the direction of a shear traction on a cylinder using global coordinates in one case and a local cylindrical coordinate system in the other case. The axis of symmetry of the cylinder coincides with the global z-axis. A surface named SURFA has been defined on the outside of the cylinder.In the first example the direction of the shear traction,, is given in global coordinates. The sense of theresulting shear tractions using global coordinates is shown in Figure 34.4.3–3(a).Figure 34.4.3–3 Shear tractions specified using global coordinates (a) and a local cylindrical coordinate system (b).*STEPStep 1 - Specify shear directions in global coordinates...*DSLOADSURFA, TRSHR, 1., 0., 1., 0. ！suefa是自定义名称；trshr为载荷类型标签，见附表一；1.为载荷量，见关键字*DSLOAD； 0., 1., 0.为投影前向量。
abaqus inp文件的详解
对ABAQUS例子的理解很多人学习ABAQUS很长时间但是却不能编写一个INP文件，在论坛中有位朋友编写了一个INP文件，但是依然有朋友问是怎么编写的，下面是我对那个例子的解释，也许会对有的朋友有些帮助，当然我的理解也可能不对，那就请斑竹和大虾指点。
我不明白的我已经在里面注明。
参见原文件可以看出，一个好的INP文件的顺序应该是这样的（本人的理解）首先定义节点，然后定义单元，再定义材料，然后定义边界条件，这是模型数据。
接下来就是历史数据，关键就是步骤的定义，当然我们需要的那些数据的输出是我们下一步进行工作的资源是一定要定义好的。
其实一个好的INP文件中在模型数据的工作中的目的就是为了得到好的网格，历史数据的目的就是得到我们想要得到的数据。
当然了有了CAE我们不需要编写INP来工作，但是对刚刚接触和学习ABAUQS的朋友来说，编写一个好的INP文件既能有一种学习的成就感也能很好的对问题有个好的理解，对学习和使用CAE来分析大型的模型是有帮助的。
附件是我对原文件的理解，请对指教。
不建议手写数据文件，可以用CAE生成，用HM生成这不是我们研究的核心，但强烈建议用手写Hitory Data，有助于加深对问题的理解*HEADINGTHE PLANAR(TWO DIMENSIONAL PROBLEM)UNITS:LENGTH-MM FORCE-N STRESS-N/MM2****THE DEFINITION OF NODE （节点的定义）****LEFT EDGE （左边的定义或者说是产生一条左边）*NODE （节点的定义，*NODE关键句定义的其实是一些独立的节点；下面的解释：1,0,0 （节点1，坐标是（0，0））30,30,0 （节点30，坐标是（30，0））*NGEN,NSET=BOTTOMEDGE （*NGEN，关键句产生一个节点集，在这个节点集中所1,30,1 使用的节点中1是初始节点，30是终点，第三个数字1是它们之间的增量。
abaqus系列之inp文件讲解
如何写input文件一、输入文件的组成和结构：1．一个输入文件由模型数据和历史数据两部分组成.模型数据的作用:定义一个有限元模型.包括单元,节点,单元性质,定义材料等等有关说明模型自身的数据.模型数据可被组织到零件中(零件可以被组装成一个有意义的模型).历史数据的定义是模型发生了什么----事情的进展,模型响应的荷载,历史被分成一系列的时步层序.每一步就是一个响应(静态加载,动态响应等)，时步的定义包括过程类型(比如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出要求.At a minimum the model consists of the following information: geometry, element section properties, material data, loads and boundary conditions, analysis type, and output requests.2．ABAQUS输入文件的结构形式。
1) 必须有一个*HEADING开头。
2)接下来就是模型数据部分，定义节点，单元，材料，初始条件等。
模型数据的层次为：部件，组装，模型。
必须的模型数据：（1）几何数据：模型的几何形状是用单元和节点来定义的，结构性单元的截面是必须定义的。
比如梁单元。
特殊的特征也可以用特殊的单元来定义，比如弹簧单元，阻尼器，点式群体等。
（2）材料的定义：材料必须定义比如使用的是钢啊，岩石，土啊等材料。
可选的模型数据：（1）零件和组合：一个模型可以用几个零件来定义有可以把几个零件组合成一个集来定义。
（2）初始条件：比如初始应力，温度，或者速度等（3）边界条件：（4）运动约束（5）相互作用（6）振幅定义（7）输出控制（8）环境特性（9）用户子程序（10）分析附属部分3)接下来就是历史数据：定义分析的类型，荷载，输出要求等。
分析的目的就是预测模型对某些外部荷载或者某些初始条件的反映。
abaqus-inp编程
1, 1, 2 (单元类型的参数)
fig.3 对称边界条件的施加及载荷的分布
input 文件祥解：
*HEADING
STRESS ANALYSIS FOR A PLATE WITH A HOLE **文件 名
的前面的是模型数据，后面的就是历史数据。 必须的历史数据：
响应类型：必须立刻出现在*STEP 选项后面。ABAQUS 中 有两种响应步，一种是总体分析响应步，可以是线形和非线 形的，另一种是线形扰动步。
可选历史数据： （1）荷载：通常定义某种加的荷载类型和大小。荷载可以 被描述成时间的函数。 （2）边界条件输出控制 （3）辅助控制 （4）再生单元和曲面 二、书写 input 文件的语法和规则： 1．关键词行： 1）必须以*开始，后面接的是选项的名字,然后随之定义选项 的内容.如: * MATERIAL NAME=STEEL 注释行是以**开始的. 2）如果有参数，则参数和关键词之间必须用“，”格开。 3）在参数之间必须用“，”格开。 4）关键词行中的空格可以忽略。
也可以用特殊的单元来定义，比如弹簧单元，阻尼器，点式 群体等。
（2）材料的定义：材料必须定义比如使用的是钢啊，岩石， 土啊等材料。
可选的模型数据：
（1）零件和组合 ：一个模型可以用几个零件来定义有可以 把几个零件组合成一个集来定义。
（2）初始条件 ：比如初始应力，温度，或者速度等
（3）边界条件：
【ABAQUS笔记】ABAQUS中如何提取变形后节点的空间坐标？后附在matlab中处理.。。。
【ABAQUS笔记】ABAQUS中如何提取变形后节点的空间坐标？后附在matlab中处理.。
1.新建⼀个ABAQUS⼯作界⾯，点选顶部菜单栏File-Import-Part。
2.弹出对话框下⾯的⽂件类型选“.odb”，选中需要提取节点坐标的结果⽂件。
3.弹出的对话框中勾选“Import deformed configuation”，选择需要输出的step和Frame，得到需要输出坐标的⽬标。
4.在Assembly ⾥将该Part实例化。
5.新建⼀个Job，不需要提交，点击右上⾓“Write Input”，即可成功在⼯作⽬录下输出⼀个.inp⽂件。
6.⽤记事本打开.inp⽂件，如下图所⽰，数据栏第⼀列为节点编号，第⼆三四列为空间坐标。
将该⽂件另存为.txt格式，在matlab⼯作⽬录下保存⼀份.txt，运⾏后⾯的程序即可在matlab中绘制出提取的节点。
（附matlab程序段，程序段⾮本⼈⼿写，来源于同课题组师弟）clc;clear all;for num=1fid1 = fopen('Job151030IN.txt'); %该⽂件来源于INP⽂件另存str = fread(fid1,'*char')';fclose(fid1);idstart = strfind(str,'*Node');idend = strfind(str,'*Element');fid2 = fopen('node_get_process1.txt','wt'); %新建⼀个保存节点编号和空间坐标的⽂件 fprintf(fid2,'%s',str(idstart+7:idend-3));fclose(fid2);data1 = importdata('node_get_process1.txt');nodelist_done_1 = data1;nodenumber1 = size(data1,1);nodelist11(1:nodenumber1,1)=nodelist_done_1(1:nodenumber1,2);nodelist11(1:nodenumber1,2)=nodelist_done_1(1:nodenumber1,3);nodelist11(1:nodenumber1,3)=nodelist_done_1(1:nodenumber1,4);figurescatter3(nodelist11(:,1),nodelist11(:,2),nodelist11(:,3),'o')axis equalxlabel('X轴');ylabel('Y轴');zlabel('Z轴');hold on;end;。
abaqus的inp文件经验总结
例子一：悬臂梁输入文件的开始就是文件头，以HEADING开始，随后是模型的名字，如下所示：*HEADINGBEAM然后是网格定义：现在就是模型数据的开始了。
一般选择从网格的定义开始，网格包括（单元和节点）假如我们的悬臂梁有五个单元，六个节点，下面我们首先详细说明节点：*NODE , NSET=ENDS1，0.6, 100.NGEN1, 6节点组集，NSET其值（名字）为ENDS.下面的就是这样理解的，第一个节点是从0开始的，第六个节点是在100结束的. 同样我们来定义单元：*ELEMENT, TYPE=B21(单元类型)1， 1， 2 （单元类型的参数）*ELGEN, ELSET=BEAM(产生单元集，及其名称)1， 5 （一个单元集，包括5个单元）现在定义单元的性质：*BEAM SECTION, SECTION=RECTANGULAR, ELSET=BEAM, MATERIAL=STEEL1.,2.梁截面，截面的形状是矩形，单元集的名称是梁单元，材料是钢。
截面尺寸是1*2.下面定义材料的性质：*MATERIAL, NAME=STEEL*ELASTIC,材料是钢，弹性模量是30E6。
下面定义边界：*BOUNDARY6， ENCASTRE边界也可以用下面的形式来定义：*BOUNDARY6, 1, 6ABAQUS对结构单元的节点的自由度使用常规的编号。
1，2，3代表的是位移分量；4，5，6代表的是旋转分量。
以上是模型数据的定义，下面开始历史数据的定义：（加载的次序，事件的发生，还有我们想看到的变量的响应）时间（步骤）的定义。
*STEP, PERTURBATION(步骤的开始，扰动是其名称)*STATIC (静态分析)*CLOAD (集中载荷)1， 2， -2000. （在节点1，y方向施加载荷，载荷大小是-20000.也就是向下施加载荷。
[2代表的是方向，假如是1则代表的是x方向，3代表的是z方向]）*END STEP（步骤的结束）下面来解释输出要求：*EL PRINT, POSITION=AVERAGED AT NODES, (节点的平均值) SUMMARY=YES （在表的下部求和）S11, (积分点的应力分量在X方向) E11，（在积分点的应变分量在X方向）SF, (在积分点的截面力)*NODE FILE,NSET=ENDSU, CF, RF------------------文件的执行1、数据的检查abaqus job=tutorial datacheck interactiveor abaqus datacheck job=frame interactive2、运行abaqus job=tutorial interactiveor abaqus job=tutorial continue interactive3、后处理运行abaqus viewer启动后处理模块。
调整abaqus中节点坐标的方法
Abaqus Analysis User's Manual2.1.6 Adjusting nodal coordinatesProducts: Abaqus/Standard Abaqus/Explicit Abaqus/CAEReferences∙ *ADJUST∙ “Defining adjust points constraints,” Section 15.15.5 of the Abaqus/CAE User's Manual OverviewNodal adjustment is used for:∙ adjusting user-specified nodal coordinates so that the nodes lie on a given surface; and∙ specifying the direction along which the nodes are moved. Adjusting nodal coordinatesIn general, user-specified nodal coordinates are not modified during input file processing. However, there are some situations where mesh coordinates are known only in a generic way and it is inconvenient to determine their coordinates for their actual usage. For example, when using fasteners the specified reference node should be positioned at its projection point on the associated surface. Since that location may be known only approximately, you can use nodal adjustment to move the reference node to that location automatically. For typical usage of the nodal adjustment feature, refer to “About assembled fasteners,” Section 29.1.3 of the Abaqus/CAE User's Manual .When using this feature, the nodes are adjusted to lie on the specified surface without regard for shell thickness or shell offsets. Therefore,it is not advisable to use this feature as a way of correcting initial overclosures for contact or for tie constraints. In addition, care should be taken when choosing the nodes to be adjusted because the feature does not respect any constraints relating the relative position of the adjusted node with other nodes (e.g., rigid body definitions).Input File Usage: U se the following option to identify thenodes to be moved and the surface onto whichthe nodes are to be moved:*ADJUST, NODE SET=name, SURFACE=nameAbaqus/CAE Usage: U se the following option to move the control pointof a coupling constraint onto the coupling surface:Interaction module: Constraint Create: Coupling;Adjust control point to lie on surfaceUse the following option to move any point or pointsonto any surface:Interaction module: Constraint Create: AdjustpointsSpecifying the nodal adjustment directionA node can be moved to the surface using a normal adjustment or a directed adjustment. By default, the node is adjusted to the closest point on the specified surface along the normal to the surface. You can specify an orientation to move the node to the surface along a given direction rather than along the normal to the surface. The vector along the localZ-direction from the orientation definition is used to move the node to the surface (see “Orientations,” Section 2.2.5). If no projection can be found, the nodal coordinates are left unmodified.Input File Usage: *ADJUST, ORIENTATION=nameAbaqus/CAE Usage: T he orientation projection option is not supportedin Abaqus/CAE.。
ABAQUS中的坐标系
1、定义节点自由度的局部坐标系*TRANSFORM，NSET=<节点集>,TYPE=<坐标系类型>x a，y a，z a，x b，y b，z b，数据行指定了a、b两点的坐标，可以指定坐标系为直角（TYPE=R），柱坐标（TYPE=C），球坐标（TYPE=S）；用NSET指定应用在这个局部坐标系的节点集。
若一个节点基于*TRANSFORM关键词定义的局部坐标系，则所有属于该节点的数据，如边界条件、集中荷载、节点输出变量（如位移、速度、反力等）也被定义在该局部坐标系中。
在大位移分析中，此局部坐标系的方向不会随着材料的旋转而旋转。
2、定义单元局部坐标系*ORIENTATION, NAME=<局部坐标系名称>，SYSTEM=<局部坐标系类型> This option is used to define a local coordinate system for definition of material properties; for material calculations at integration points; for element property definitions (e.g., connector elements); for output of components of stress, strain, and element section forces; and for kinematic and distributing coupling constraints.用于定义材料特性，定义单元应力应变的输出方向、单元截面力输出方向；实体单元默认的材料方向为整体直角坐标系；壳单元、膜单元默认的材料方向为整体直角坐标系在壳或膜单元表面的投影；在大位移分析中，此局部坐标系的方向会随着材料的旋转而旋转。
*TRANSFORM与*ORIENTATION定义的两个局部坐标系是相互独立的。
abaqus中定义管道局部坐标系、切线载荷
1.CAE操作对应inp文件：*Orientation, name="Datum csys-2", system=CYLINDRICAL-20., -2.5, 20., -20., -2.5, 21. ** Name: Load-1 Type: Surface traction*Dsload, orientation="Datum csys-2", constant resultant=YESSurf-1, TRSHR, 10., -0.707107, 0.707107, 0.2.操作原理：using a local coordinate system to define shear directionsIt is sometimes convenient to give shear and general traction directions with respect to a local coordinate system. The following two examples illustrate the specification of the direction of a shear traction on a cylinder using global coordinates in one case and a local cylindrical coordinate system in the other case. The axis of symmetry of the cylinder coincides with the global z-axis. A surface named SURFA has been defined on the outside of the cylinder.In the first example the direction of the shear traction,, is given in global coordinates. The sense of theresulting shear tractions using global coordinates is shown in Figure 34.4.3–3(a).Figure 34.4.3–3 Shear tractions specified using global coordinates (a) and a local cylindrical coordinate system (b).*STEPStep 1 - Specify shear directions in global coordinates...*DSLOADSURFA, TRSHR, 1., 0., 1., 0. ！suefa是自定义名称；trshr为载荷类型标签，见附表一；1.为载荷量，见关键字*DSLOAD； 0., 1., 0.为投影前向量。
abaqusinp文件的详解
abaqusinp文件的详解对ABAQUS例子的理解很多人学习ABAQUS很长时间但是却不能编写一个INP文件，在论坛中有位朋友编写了一个INP文件，但是依然有朋友问是怎么编写的，下面是我对那个例子的解释，也许会对有的朋友有些帮助，当然我的理解也可能不对，那就请斑竹和大虾指点。
我不明白的我已经在里面注明。
参见原文件可以看出，一个好的INP文件的顺序应该是这样的（本人的理解）首先定义节点，然后定义单元，再定义材料，然后定义边界条件，这是模型数据。
接下来就是历史数据，关键就是步骤的定义，当然我们需要的那些数据的输出是我们下一步进行工作的资源是一定要定义好的。
其实一个好的INP文件中在模型数据的工作中的目的就是为了得到好的网格，历史数据的目的就是得到我们想要得到的数据。
当然了有了CAE我们不需要编写INP 来工作，但是对刚刚接触和学习ABAUQS的朋友来说，编写一个好的INP文件既能有一种学习的成就感也能很好的对问题有个好的理解，对学习和使用CAE 来分析大型的模型是有帮助的。
附件是我对原文件的理解，请对指教。
不建议手写数据文件，可以用CAE生成，用HM生成这不是我们研究的核心，但强烈建议用手写Hitory Data，有助于加深对问题的理解*HEADINGTHE PLANAR(TWO DIMENSIONAL PROBLEM)UNITS:LENGTH-MM FORCE-N STRESS-N/MM2****THE DEFINITION OF NODE （节点的定义）****LEFT EDGE （左边的定义或者说是产生一条左边）*NODE （节点的定义，*NODE关键句定义的其实是一些独立的节点；下面的解释：1,0,0 （节点1，坐标是（0，0））30,30,0 （节点30，坐标是（30，0））*NGEN,NSET=BOTTOMEDGE （*NGEN，关键句产生一个节点集，在这个节点集中所1,30,1 使用的节点中1是初始节点，30是终点，第三个数字1是它们之间的增量。
ABAQUS读取节点集坐标
在abaqus中使用python实现的功能(2010-02-06 15:04:10)转载▼标签：杂谈功能一：实行提交多个job的功能。
对象：Job object使用：在源文件开始写上import job，源程序用mdb.jobs[name] 使用名字为name的job对象。
建立一个job对象的方法：λ利用已有的inp文件中建立job：mdb.JobFromInputFile（）λ利用已有的cae中建立job： Job(...)建议用第一种方法。
设定参数的方法：⎫利用第一种方法建立job的时候，可以设定很多的参数，比如type，queue，userSubroutine等。
格式：mdb.JobFromInputFile（name＝，inputFile＝，type＝，queue＝，userSubroutine＝，…….）。
⎫也可以先建立一个job，然后利用job对象的setValues来设定参数，格式：job.setValues(type＝，queue＝，userSubroutine＝，…….)。
一个简单的例子：文件：job.pyfrom abaqusConstants import *import jobmdb.JobFromInputFile(name='job-1-1',inputFileName='Job-1.inp')＃基于inp文件Job-1.inp建立名称为job-1-1的jobmdb.jobs['job-1-1'].setValues(waitMinutes=1)＃设定参数mdb.jobs['job-1-1'].submit()＃提交任务mdb.jobs['job-1-1'].waitForCompletion()运行：在cmd下面运行：Abaqus cae nogui=job.py如果是多个job，同样道理了，不多说了。
abaqus-inp参数化编程
对入门来说还是相当不错的。
abaqus inp文件精讲如何写input文件一、输入文件的组成和结构：1．一个输入文件由模型数据和历史数据两部分组成.模型数据的作用:定义一个有限元模型.包括单元,节点,单元性质,定义材料等等有关说明模型自身的数据.模型数据可被组织到零件中(零件可以被组装成一个有意义的模型).历史数据的定义是模型发生了什么----事情的进展,模型响应的荷载,历史被分成一系列的时步层序.每一步就是一个响应(静态加载,动态响应等)，时步的定义包括过程类型(比如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出要求.At a minimum the model consists of the following information: geometry, element section properties, material data, loads and boundary conditions, analysis type, and output requests.2． ABAQUS输入文件的结构形式。
1) 必须有一个*HEADING开头。
2) 接下来就是模型数据部分，定义节点，单元，材料，初始条件等。
模型数据的层次为：部件，组装，模型。
必须的模型数据：（1）几何数据：模型的几何形状是用单元和节点来定义的，结构性单元的截面是必须定义的。
比如梁单元。
特殊的特征也可以用特殊的单元来定义，比如弹簧单元，阻尼器，点式群体等。
（2）材料的定义：材料必须定义比如使用的是钢啊，岩石，土啊等材料。
可选的模型数据：（1）零件和组合：一个模型可以用几个零件来定义有可以把几个零件组合成一个集来定义。
（2）初始条件：比如初始应力，温度，或者速度等（3）边界条件：（4）运动约束（5）相互作用（6）振幅定义（7）输出控制（8）环境特性（9）用户子程序（10）分析附属部分3) 接下来就是历史数据：定义分析的类型，荷载，输出要求等。
abaqus二次开发输出结点坐标
abaqus⼆次开发输出结点坐标from abaqus import *from abaqusConstants import *import __main__import sectionimport regionToolsetimport displayGroupMdbToolset as dgmimport partimport materialimport assemblyimport stepimport interactionimport loadimport meshimport optimizationimport jobimport sketchimport visualizationimport xyPlotimport displayGroupOdbToolset as dgoimport connectorBehaviorfrom driverUtils import executeOnCaeStartupfrom decimal import *from operator import *from string import *import datetimeimport mathimport meshEditimport osimport randomimport shutilimport timeimport numpy################################################################################### CREATE PART##################################################################################Mdb()s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=1.0)g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraintss.setPrimaryObject(option=STANDALONE)s.rectangle(point1=(0.5, 0.5), point2=(-0.5, -0.5))s.CircleByCenterPerimeter(center=(0.0, 0.0), point1=(0.0, 0.3))p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=TWO_D_PLANAR,type=DEFORMABLE_BODY)p = mdb.models['Model-1'].parts['Part-1']p.BaseShell(sketch=s)s.unsetPrimaryObject()del mdb.models['Model-1'].sketches['__profile__']################################################################################## # DEFINE MATERIAL AND SECTION################################################################################## mdb.models['Model-1'].Material(name='Material-1')mdb.models['Model-1'].materials['Material-1'].Elastic(table=((200.0, 0.3), ))mdb.models['Model-1'].HomogeneousSolidSection(name='Section-1',material='Material-1', thickness=None)f = p.facesregion = p.Set(faces=f, name='Set-1')p.SectionAssignment(region=region, sectionName='Section-1', offset=0.0,offsetType=MIDDLE_SURFACE, offsetField='',thicknessAssignment=FROM_SECTION)################################################################################## # ASSEMBLE################################################################################## a = mdb.models['Model-1'].rootAssemblya.DatumCsysByDefault(CARTESIAN)a.Instance(name='Part-1-1', part=p, dependent=ON)################################################################################## # SEED PART################################################################################## # create partitionp = mdb.models['Model-1'].parts['Part-1']f1, e, d1 = p.faces, p.edges, p.datumst = p.MakeSketchTransform(sketchPlane=f1[0], sketchPlaneSide=SIDE1, origin=(0.0, 0.0, 0.0))s1 = mdb.models['Model-1'].ConstrainedSketch(name='__profile__',sheetSize=2.82, gridSpacing=0.07, transform=t)g, v, d, c = s1.geometry, s1.vertices, s1.dimensions, s1.constraintss1.setPrimaryObject(option=SUPERIMPOSE)p.projectReferencesOntoSketch(sketch=s1, filter=COPLANAR_EDGES)s1.Line(point1=(-0.5, 0.0), point2=(0.5, 0.0))s1.Line(point1=(0.0, 0.5), point2=(0.0, -0.5))e1, d2 = p.edges, p.datumsp.PartitionFaceBySketch(faces=f, sketch=s1)s1.unsetPrimaryObject()del mdb.models['Model-1'].sketches['__profile__']# mesh controlp.seedPart(size=0.05, deviationFactor=0.1, minSizeFactor=0.1)p.setMeshControls(regions=p.faces, elemShape=TRI, technique=STRUCTURED) p.generateMesh()# it is importanta = mdb.models['Model-1'].rootAssemblya.regenerate()Nodeset = mdb.models[NameModel].rootAssembly.instances[Nameinst].nodes j = 0x = []y = []for i in Nodeset:x.insert(j, i.coordinates[0])y.insert(j, i.coordinates[1])j = j + 1。
Abaqus INP输出接触力-接触力矩-摩擦力-接触面积-接触中心坐标
Abaqus输出接触对的接触力(CFN)、接触力矩(CMN)、摩擦力（CFS）等
关键字： *Output, history *Contact Output CFN1, CFN2, CFN3, CFNM, CMN1, CMN2, CMN3, CMNM, CFS1 说明： 1. 接触力、摩擦力只能在历史变量中输出 2. 如果按照上述方法书写INP，则CFN输出的是包含所有接触对的接触力 3. 如果输出指定接触对的接触力，则后边需要指定接触对的主从面： *Output, history *Contact Output, Mater=A-Cnt-B, Slave=B-Cnt-A CFN1, CFN2
4. 接触力、接的接触力力矩： CFM， 全局坐标系下输出的,且是对全局坐标系原点的力矩 即：包括了该接触对处本身的力矩，以及该点（XN）的力对原点的力矩 所以要求该处本身的力矩，需要用输出的CFM 减去力对原点的力矩
Abaqus关于接触对的常用输出变量