《abaqus输出文件fil-dat-刚度矩阵等关键词及相关例子.doc》由会员分享，可在线阅读，更多相关《abaqus输出文件fil-dat-刚度矩阵等关键词及相关例子.doc（12页珍藏版）》请在三一文库上搜索。
1、输出fil文件关键词*Output, history, variable=PRESELECT, time in terval=0.01*EI file, freque ncy=999*Node file, freque ncy=999*El fileS,e*Node fileCF, TF, Ufrequency=1表示每一个时间增量步输出一次frequency=999表示每999个时间增量步输出一次当frequency值足够大时，只在分析步的末端结束时输出数据输出dat文件关键词*Output, history, variable=PRESELECT, time in terval=0.01*
2、El print, freque ncy=999*Node print, freque ncy=999*El printS11,e11*Node printCF, TF, UEl print El file将单元上的分析结果（应力、应变和截面力等输出）Node file Node print将单元上的分析结果（位移、反力输出）输出单元刚度矩阵的方法！*Element Matrix Output,Elset=AII, frequency=nFile Name=user defined, file name= abc, Stiffness=Yes mass=yes， DLOAD=YESfreque
3、ncy=n每隔_n个增量步输出一次单元矩阵，_Elset=All all为单元集合名称集合可以在 assembly instanee中进行定义Stiffness=Yes 表示输出刚度矩阵的的算子矩阵mass=yes表示输出质量矩阵DLOAD=YES表示输出载荷向量no表示输出输出单元刚度矩阵的 最终正确方法！*File Format,Ascii*Element N 力atri x Output,Elset=E1,File Name=abc,Frequency=50, -Output File=User Defi ned,Stiffness=YesE:xch-cae-nsoft'jied
4、ianli-yanjiuABC.mtx生成文件.mtx 所在路径"FIELD OUTPUT: F-Output-2B-fr*Node OutputCFd RE TE UHELD C JTPUT: F-Output- 3Tlement Output, directions= YESLEd PE. PE EQ PEMAG, 5* HISTORY OUTPUT; H-Output-1Output, history, ¥0 ria bl 亡三 PRE SELECTFile FormatAsniElement Matrix Output,Else-t= El,Ale Mame=abc
5、rFrequency = 50,Output File-User &frfinedrStiffn<ss=Y&s"End S-tep十ABAQUS输出实体单元某一截面的轴力和弯矩一、方法在INP文件中进行编辑；在*Output, history, variable=PRESELECT 后面添加以下语句 *Sect ion Print, n ame=a1, surface=asd, axes=local（空两行）SOM, SOF1、中间要空开两行；2、 name名字自己定义，surface的名字是自己定义的那个面；3、算完后在DAT文件里面有SOM和SOF的输出数值
6、；4、cutting surface需要事先定义好才可以输出该位置的轴力和弯矩可以用以下方法：方法一：*定义切平面*surface,type=cutti ng surface ,n ame=fifth-ri ng-up*定义输出一一（SOM为截面总弯矩，SOF为截面合力）*secti on prin t, name=forcemome nt,surface=fifth-ri ng-up,update=yesSOM, SOF注意：关键词 surface和 section print 的用法可以参见 abaqus帮助 /Abaqus Keywords Referrenee Manual 以及 Ab
7、aqus/CAE User's Manual 等。方法二：首先需要导出节点应力， 节点应力可以：（1 ）从Visualization/Report/Field output里选择Unique Nodal导出；（2）由于节点应力不能直接在 inp中关键词输出，但是可以在 inp文件中将节 点附近单元应力导出，然后取平均值即可；语法为：*el print, elset=fifth-ring-outside-els将截面内外节点的应力值求出后，用材料力学方法求解。ABAQUS中如何通过cutting surface和section print输出桩的轴力经过两个星期的摸索与学习，今天终于学会
8、了桩轴力的输出。现总结如下：1.主要步骤是先定义截面 cutting surface，然后用section print输出轴力sof。2所有操作均是在inp文件中进行修改的，而不是 ABAQUS/CAE 中的编辑关键词（editkeywords）。原因：在CAE中编辑关键词是可以修改inp文件，但CAE并不能识别所有的inp文件关键词，下面将举例说明。3.最后提交的inp文件也不是在 CAE中导入模型文件（import model），然后提交job进行运 算的，而是在 ABAQUS命令窗口（小黑屏）进行的。原因同2中的一样，CAE并不能识别关键词*section print。好了，下面开始详细
9、的步骤讲解吧！第一步：定义截面（cutting surface）,具体的关键语句为：*surface,type=cutt ing surface ,n ame=cutsurf-10.6,25,0,0,1,0Set-pile解读：第一行，定义 surface、surface类型以及名称。第二行，定义截面上的一点（0.6,25,0）以及截面的法向量（0,1,0）。法向量不一定是单位向 量。第三行，截面所在的单元或集合。这个集合可以是事先在CAE里定义好的。需要注意的是，这个 cutting surface是垂直于桩径方向的一个桩截面，而不是桩的侧表面。 我一开始理解错了！此关键句在inp文件中的位
10、置是在SECTIONPRINTThe model "Job-20101215" has bee n imported from an in put file.Please scroll up to check for error and warning messages.这再次验证了， ABAQUS/CAE是不能识别关键词*section print的。强制提交job是可以运 算的，但是在最后的 dat结果文件里没用任何输出结果。第三步：保存并提交inp文件。在保存之前需要仔细检查上面关键语句的位置及拼写，以及是否含有中文字符。Assembly, name=Assembly
11、这一行之后，我试了下在放在 材料定义之后，运算不成功。估计是因为我的Set-pile是在assembly里定义的，而不是在part里定义的单元集。若是在part里定义了桩的集合，是不是可以将此关键句移动到材料定义之后，这个我倒没试。如果是在CAE中通过编辑关键词来添加上述语句，将会有下面的错误提示：第二步：定义输出（section print），具体的关键语句为：*secti on prin t, name=secpri nt-1,surface=cutsurf-1sof,som解读：第一行，定义输出的名称及截面。第二行，定义输出选项，sof为截面合力，som为截面弯矩。注意，此处第一行 su
12、rface=然后就可以在 ABAQUS命令窗口提交inp文件进行运算了。具体操作是：开始宀所有程序必须是上一步中 cutting surface的名称（name=*），因为要输 出的就是第一步定义的截面轴力及弯矩。如果此处乱写surface=？的话，估计到最后什么也输出不了了。有些人还在第一句后面加上了axes=local,freque ncy=1,update=yes，由于我的模型较简单，也就没管这几个东西。除了命名之外，inp文件中的大小写英文字符，是没有区别的。所以大家不要再犹豫写*SECTION PRINT 还是*section print此关键句在inp文件中的位置是最后一个* HI
13、STORY OUTPUT: H-Output-3 之后，*Output,history, variable=PRESELECT之前。也就是我的inp文件中的倒数第四和第三行。如果放错 了位置，估计也会引起不识别等问题，最终导致无输出结果。如果是在CAE中通过导入inp文件来进行运算，在导入时，CAE下方窗口将会有下述警告信息：WARNING:Thefollowi ng keywords/parameters are not yet supported by the in put file reader:t Abaqus 6.9.1 Abaqus Comma nd输入：abaqus job=*.
14、inp （就是你刚才修改的inp文件名）回车即可。然后黑屏窗口会弹出一个警告语句：Abaqus Warnin g:The .inp exte nsion has bee n removed from the job inden tifier大概意思是说，后缀.inp被省略。由此可判断，在上面输入的时候，完全可以不写后缀名字。我之所以写了，是因为我的默认工作路径 D:Temp下还有个跟inp文件完全同名的.cae文件， 而如果这两个文件信息不是完全一样的话，建议还是加上.inp后缀为好。到此，计算完毕。黑屏窗口并不会像 CAE那样提示job completed或ANSYS中的solution i
15、s done !（害得我还以为没反应 -_-汗！）默认工作路径下已经生产了结果文件.com .dat .log .msg .odb .prt .sta 等一系列后缀文件。第四步：提取数据用记事本打开刚才生产的.dat文件，在最下面会找到想要的输出结果THE FOLLOWING TABLES ARE PRINTEDSOFMSOF1SOF2SOF31.0556E+064.0515E+041.0548E+060.000三个方向的sof,以及合力SOFM。需要注意的是，此 sof值是合应力，单位为 N。无需乘以截面面积至此，大功告成！后话：（1）本人查看了仿真论坛里的一些帖子，发现输出轴力的方法有三：
16、?后处理模块中运用free body cut，此法较为方便简单，可直接得到截面合力，无奈我的模型是二维轴对称的，不能进行自由体切割。?输出截面上的应力值进行积分运算，主要方法是创建路径（path ）或者拾取节点应力值，桩身的网格越密集，最后的积分值越精确，此法我已学会。?cutting surface和section print联合使用，由于手动修改inp文件，此法最复杂，而且不知道最终的sof值，软件是如何算出来的。如果需要得到桩深方向N个截面的轴力，就需要在cutting surface里定义N个截面。Abaqus中in put文件介绍1. 本INP文件以及解释参考石亦平博士所著的ABAQ
17、US有限元分析实 例详解；2. 注释行以INP文件总是以*Heading开头，接下来可以用一行或多行来写下此模型的标 题和相关信息. 口头，以保持与INP文件的格式一致（为方便起见以绿色 显示）。*Head ing* Job name: Plate-CPS8 Model name: Model-1*Prepri nt, echo=NO, model=NO, history=NO, co ntact=NO*Preprint可设置在DAT文件（*.dat）中记录的内容。上述为 ABAQUS默认，内 容为：在DAT文件*中不记录对INP文件的处理过程，以及详细的模型和历史数据。* PARTS*Par
18、t, name=Plate1*定义Part的标准格式：*Part，name=E件名称*对于非独立实体，比如此例，要在下面数据块中定义详细的节点，单元，截面 属性等数据。以便 Mesh时，节点定义以*Node开始*节点的基本表示方法：* 一维模型：节点编号，节点坐标*二维模型：节点编号，节点坐标1，节点坐标2*三维模型：节点编号，节点坐标 1，节点坐标2，节点坐标3*注意：1.节点编号可以不从1开始也可以是不连续的。*注意：2.在INP文件中，不同的部件或实体可以有相同的节点或单元编号。比如，部件Part-A的节点编号*是1, 2, 3,部件Par-B的节点编号也可以I ，厶，xO ， *注意：
19、3.如果在定义荷载Load，边界条件BC或约束Constraint是需要引用这 些节点编号，需要加上相应*的实体名称作为前缀。例如，部件Part-A和*部件Part-B的相应 实体 分别名为Part A-1和Part-B-1,则实体Part-A-1的节 点记作Part-A-1.5,而实体*Part-B-1的节点5记作Part-B-1. 5Mesh在Part上面；*对于独立实体，则Part数据块只包含*Part和*End Part两行，而没有实质性的 数据。节点，单元，截面属*性等数据将在Instanee数据块中定义，而*Mesh时也只能在Instanee上面做。*如果INP文件是由Abaqus
20、/CAE生成的，其结构会包括部件*Part，装配件 *Assembly，实体*ln sta nee等数*据块，比如此例。*如果INP文件是由其他前处理器如 MSC, PATRAN，FEMAP等生成的，其结 构将不包含部件，装配件和实*体等数据块，而是直接定义节点和单元等数据。*Node1,0.,27.52,0.,5.3,略若干3.53553391,二节点3.53553391279,10.1391716,28.3766441280,4.99425983,29.4102631281,0.,30.3125*Eleme nt, type=CPS8*单元的定义方法：*Element, type=单元类型
21、*单元编号，节点1编号，节点2编号，节点三编号1,2,3,12,13,14,1,12,13,省略若干单元57,58,59,23, 102, 103, 104, 10557, 106, 107, 108, 10358, 109, 110, 111, 10778, 99, 100,79, 100, 101,80, 101, 56,22,23,1,21,273, 279, 126, 27822, 275, 280, 117, 27923, 277, 281, 105, 280*下面为集合的定义：集合分为节点集合*Nset和单元集合*Elset，又分为下面两 种：*1.定义在Part或Instanee
22、数据块中的集合：这类集合出现在 *Part和*丘门4 par 之间（比如此例），或出现*在*lnstanee和*丘门4 instanee之间，一般用来定义*界面属性。*2.定义在Assembly数据块的集合：这类集合出现在 *End Instanee之后，*End assembly之前，一般用来定*义荷载，边界条件，接触或约束等。*节点集合和单元集合的表示方法：号是连续的，可以表示为*Nset， Nset= 点集合名称，Ge nerate起始节点编号，结束节点编号，节点编号增量*Elset， Elset= 元集合名称，Gen erate起始单元编号，结束单元编号，单元编号增量80个字符，必须以
23、字母或下划线*1.如果集合中的节点或单元编*单元集合:*节点集合：号*注意：节点集合和单元集合的名称不得超过 开始。*2.如果集合中的节点或单元编号是不连续的，表示方法为：依次列出集合中的 所有节点或单元（每个数据行接点或单元编号不得超过16个），格式：*节点集合：*Nset， Nset=点集合名称节点编号1,节点编号2，节点编号16*单元集合：Elset, Elset=元集合名称单元编号1,单元编号2,单元编号16*Nset, n set=_PickedSet2, intern al, gen erate1, 281, 1*Elset, elset=_PickedSet2, intern a
24、l, gen erate1, 80, 1* Regi on: (Sectio n-1:Picked)*Elset, elset=_PickedSet2, intern al, gen erate1, 80, 1*截面属性的基本表示方法为：*Solid Section , Elset=元集合名称，Material材料名称*界面参数*界面参数可以是二维模型的厚度或一维模型的截面面积等。材料名称不超过 80个字符，必须以字母开头。* Sectio n: Sectio n-1*Solid Sectio n, elset=_PickedSet2, material=Steel1.,*End Part*定
25、义Assembly数据块的格式为：*Assembly, Name=S配件名称*End Assembly*省略号代表在Assembly数据块中的Instanee数据块，以及定义在Assembly数 据块中的几何数据块，以及面和约束有关的数据块。* ASSEMBLY*Assembly, n ame=Assembly*ln sta nee, n ame=Plate1-1, part=Plate1 *End In sta nee*定义在Assembly数据块中的集合表示方法与定义在Part或Instanee数据块中的集合基本相同，只是需要在其后面加上参数In sta nce体名称*Nset, n se
26、t=_PickedSet5, i nternal, i nsta nce=Plate1-11,2,11, 12,13,14, 54, 55,56,102,106,109, 112, 260, 268, 276281,*Elset, elset=_PickedSet5, i nternal, i nsta nce=Plate1-11, 2, 3, 4, 68, 72, 76, 80*Nset, n set=_PickedSet6, i nternal, i nsta nce=Plate1-15,6,7, 27, 28, 29, 33,34, 35, 152, 160, 168, 176, 178
27、, 181, 184187,*Elset, elset=_PickedSet6, i nternal, i nsta nce=Plate1-120, 24, 28, 32, 33, 34, 35, 36*Elset, elset=_PickedSurf4_S2, i nternal, i nsta nce=Plate1-1, gen erate36, 48,4*Elset, elset=_PickedSurf4_S1, intern al, i nsta nce=Plate1-1, gen erate49, 52,1*定义Surface数据块的格式为：*Surface， Type=H的类型，N
28、ame得名称*构成面的集合1，名称1*构成面的集合2，名称2*像定义节点和单元集合一样，不需要使用*End Surface等语法来结束。*Surface, type=ELEMENT, name=_PickedSurf4, i nternal _PickedSurf4_S2, S2_PickedSurf4_S1, S1*End Assembly* 定义Material数据块的格式为：*Material， ”&口6=材料名称*Elastic*弹性模量，泊松比*Plastic*屈服应力（列），塑性应变（列）*还可以定义*Density等*注意：所有使用Abaqus/Explicit的分析，都
29、要使用*Density来定义密度* MATERIALS*Material, n ame=Steel *Elastic210000., 0.3*定义Boundary Condition数据块的格式（以边界条件为例，速度等以此类推）：*Bou ndary*节点编号或节点集合，约定的边界条件类型*注意：如果一个边界条件定义在初始步（Initial Step）中则相应的Boundary数 据块出现在*Step之前；*如果一个边界条件定义在后续分析步中，则相应的Bou ndary数据块出现在此后续分析步的*Step* 和*End Step 之间。* BOUNDARY CONDITIONS* Name:
30、Fix-X Type: Symmetry/A ntisymmetry/E ncastre *Bo un dary_PickedSet5, XSYMM* Name: Fix-Y Type: Symmetry/A ntisymmetry/E ncastre *Bo un dary_PickedSet6, YSYMM*定义Step数据块的格式为（以静力分析为例）：*Step，”&口6=分析部名称*Static*初始增量步，分析时间，最小增量步，最大增量步 STEP: Apply Load*Step, n ame="Apply Load" *Static1., 1., 1e
31、-05, 1.*定义载荷Load数据块的格式：1.集中载荷：Output, history, variable=PRESELECT Cload*节点编号或节点集合，自由度编号，载荷值*2.定义在单元上的分布荷载：*Dload*单元编号或单元集合，载荷类型的代码*3.定义在面上的分布载荷：*DSload*面的名称，载荷类型的代码，载荷值*注意：描述载荷的关键词为：集中载荷 *Cload ( Concentrated Load)，定义在 单元上的分布载荷*Dload* (Distributing Load)和定义在面上* 的分布载荷 *DSload (Distributing Surface Lo
32、ad)*在Dload和DSload中的载荷类型的代码，从 ABAQUS Analysis User'sManual的Distributing loads查*找。下例的P代表均布面荷载。* LOADS* Name: Load-1 Type: Pressure *Dsload_PickedSurf4, P, -100.* OUTPUT REQUESTS*Restart, write, freque ncy=0*不输出用于重启动分析的数据。* FIELD OUTPUT: F-Output-1*Output, field, variable=PRESELECT*将Abaqus默认的场变量写入
33、ODB文件。* HISTORY OUTPUT: H-Output-1*将Abaqus默认的历史变量写入 ODB文件。*End Step从Abaqus/Sta ndard结果文件中读取单元矩阵为了从Abaqus/Standard结果文件中读取单元矩阵，必须在先前的分析中将刚度 和质量矩阵写入结果文件中作为单 元矩阵输出或子结构矩阵输出。必须指定单元号n或子结构标识符Zn,来知道与哪个矩阵对应。对于通过部件 装配得到的模型，输出到结果文件中的单元号是Abaqus/Sta ndard中产生的内部号。另外，在输出的单元矩阵中必须为其指定子步号和增量号。在子结构内部矩阵输 出时这些都不需要。输入格式：
34、*USER ELEMENT, FILE=name, OLD ELEMENT=n or Zn, STEP=n, INCREMENT=n输出单元刚度矩阵的方法！1用命令：* *ELEMENT MA TRIX OUTPUT只设定Required parameter： ELSET的话，由于结果文件（*.fil ）是二进制文件，用文本编辑器打开是一堆我们看不明白的乱码，所以有必要设置一下文件格式。2用命令：*FILE FORMAT设定文件格式将文件格式设为 ASCII，那么结果文件我们就能够看明白了。3命令：*ELEMENT MA TRIX OUTPUT 的其他参数设定：首先FREQUENCY要设定，否
35、则会生成每一步的结果文件，这个文件会很大，举例来说，1800个单元，共50步，全部输出的单元刚度矩阵有600M，这是一个很大的数据，所以要把不需要的删除。并且在打开这个600M大的文件时发现，UntraEdit比记事本和写字板的打开速度快得多！假如共50步，我只想要最后一步的单元刚度矩阵，那么可以设定FREQUENCY = 50 ,设定FREQUENCY = 0表示不输出，默认为每一步都输出。其次可以设定结果文件名， 设定时要使用 OUTPUT FILE=USER DEFINED 这个选项，将FILENAME的值设定为想要的文件名，注意不要扩展名。例如：*Elset, Elset=AII, G
36、en erate1, 1800, 1*FiIe Format,Ascii*EIeme ntMatrix Output,Elset=AII,File Name=abc,Freque ncy=50,Output FiIe=UserDefin ed,Stiffness=Yes计算完后可以在工作目录下找到文件名为ABC.mtx的单元刚度矩阵的结果文件，结果数据为十进制，大小仅为12M。*定义Instanee数据块的格式为：*Instanee, Name=体名称，部件名称*End In sta nee*在下面的格式中，不包含任何节点，单元，集合和截面属性等数据，因为文件 中的实体是非独立实体，必*须在Part模块中定义。*