下载文档原格式
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假设用户只需要平移。
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中的局部坐标系在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假设用户只需要平移。
如何写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)接下来就是历史数据：定义分析的类型，荷载，输出要求等。
分析的目的就是预测模型对某些外部荷载或者某些初始条件的反映。
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例子的理解很多人学习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是它们之间的增量。
如何写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中如何提取变形后节点的空间坐标？后附在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;。
例子一：悬臂梁输入文件的开始就是文件头，以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启动后处理模块。
