abaqus程序编写.doc

如何写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分析是建立在STEP的概念上的，（在历史数据中描述）在分析中可以定义多种STEPS每个STEP用*STEP开始，用*END STEP结束STEP是历史数据和模型数据的分界点，第一次浮现*STEP的前面的是模型数据，背面的就是历史数据 必须的历史数据：响应类型：必须立即出目前*STEP选项背面。

ABAQUS中有两种响应步，一种是总体分析响应步，可以是线形和非线形的，另一种是线形扰动步可选历史数据：（1）荷载：一般定义某种加的荷载类型和大小荷载可以被描述成时间的函数2）边界条件输出控制（3）辅助控制（4）再生单元和曲面二、书写input文献的语法和规则：1．核心词行：1）必须以*开始，背面接的是选项的名字,然后随之定义选项的内容.如:* MATERIAL NAME=STEEL注释行是以**开始的.2）如果有参数，则参数和核心词之间必须用“，”格开3）在参数之间必须用“，”格开4）核心词行中的空格可以忽视5）每行的长度不能超过256个字符6）核心词和参数对大小写是不辨别的7）参数值一般对大小写也是不辨别的，但是唯一的例外是文献名辨别大小写8）核心词和参数，以及大多数状况下的参数值是不需要全拼写出来的，只要她们之间可以互相辨别就可以了9) 如果参数有响应的值，则付值号是“=”10) 核心词行可以延续，例如参数的名字很长，要在下一行继续这个核心词行的话就可以这样做，用“，”来连接例如：*ELASTIC, TYPE=ISOTROPIC,DEPENDENCIES=111) 有些选项容许INPUT和FILE的参数作为一种输入文献名，这样的文献名必须涉及一种完整的途径名或者是一种相对途径名。

2．数据行：（数据行如果和核心词相联系必须紧跟核心词行）1） 一种数据行涉及空格在内不能超过256个字符2） 所有的数据条目之间必须用“，”格开3） 一行中必须涉及指定阐明的数据条目的数字4） 每行结尾的空数据域可以省略5） 浮点数最多可以占用20个字符6） 整数可以是10个7） 字符串可以是80个8） 延续行可以被用到特定的状况3．标签：什么是标签呢，例如集名，曲面名，钢筋名，她们是辨别大小写的，长度可以有80个字符长标签中的空格是可以省略的，除非用“”来标示那就不能省略了没有用“”来标示的标签必须用字母来开头如果一种标签用“”来定义那么“”也是标签的一部分标签的开始和结束不能用双重“＿”下面是一种使用了引号和没有使用引号的例子：*ELEMENT, TYPE=SPRINGA, ELSET="One element"1,1,2*SPRING, ELSET="One element"1.0E-5,*NSET, ELSET="One element", NSET=NODESET*BOUNDARYnodeset,1,24．数据行反复:数据行可以反复,也就是说每行数据可以有一行响应的变量也可以有几行。

同样也可以有多行数据行,相应各自的变量行.如：*ELASTIC, TYPE=ISOTROPIC 200.E3, 0.3, 20. 150.E3, 0.35, 400. 80.E3, 0.42, 700.定义了一种零件的材料性质,均质,线弹性,在不同应力下的扬式模量和泊松比.三、从外存储器中引入模型或者历史数据The *INCLUDE option can be used to specify an external file that contains a portion of the ABAQUS input file. This file can include model and history definition data, comment lines, and other *INCLUDE references. When the option is encountered, ABAQUS will immediately process the input data within the file specified by the INPUT parameter. When the end-of-file is reached, ABAQUS will return to processing the original file.A maximum of five levels of nested *INCLUDE references can be used. The option supports UNIX environment variables.Usage: *INCLUDE, INPUT=file_name四、举例：下面简介两个简朴的例子:例一：悬臂梁.输入文献的开始就是文献头,以HEADING 开始,随后是模型的名字,如下所示:*HEADING CANTILEVER BEAM然后是网格定义:目前就是模型数据的开始了.一般选择从网格的定义开始,网格涉及(单元和节点)如果我们的悬臂梁有五个单元,六个节点,下面我们一方面具体阐明节点:*NODE, NSET=ENDS 1, 0. 6, 100.*NGEN 1, 6节点组集,NSET其值(名字)为ENDS.下面的就是这样理解的,第一种节点是从0开始的,第六个节点是在100结束的.同样我们来定义单元:*ELEMENT, TYPE=B21(单元类型) 1, 1, 2 (单元类型的参数)*ELGEN, ELSET=BEAM (产生单元集,及其名称) 1, 5 (一种单元集,涉及5个单元)目前定义单元的性质：*BEAM SECTION, SECTION=RECTANGULAR, ELSET=BEAM, MATERIAL=STEEL 1., 2.梁截面，截面的形状是矩形，单元集的名称是梁单元，材料是钢。

截面的尺寸是1*2下面定义材料的性质：*MATERIAL, NAME=STEEL*ELASTIC 30.E6,材料是钢，弹性，弹性模量是30E6下面定义边界：*BOUNDARY 6, ENCASTRE边界是在6节点，通过ENCASTRE来描述边界也可以用下面的形式来定义：*BOUNDARY 6, 1, 6ABAQUS对构造单元的中的节点的自由度使用常规的编号方式1，2，3 代表的是位移分量；4，5，6代表的是旋转分量以上是模型数据的定义，下面开始历史数据的定义：（加载的顺序，事件的发生，尚有我们想看到的变量的响应）时步（环节）的定义STEP, PERTURBATION（环节的开始，扰动是其名称）*STATIC（静态分析）*CLOAD（集中荷载） 1, 2, -0.（在节点1，y[2代表的是y方向，如果是1则代表的是x方向，3代表的是z方向]方向施加荷载，荷载的大小是－0也就是向下施加荷载END STEP（环节的结束）下面来解释输出规定：*EL PRINT, POSITION=AVERAGED AT NODES, （节点的平均值）SUMMARY=YES在表的下部求和） S11, （积分点的应力分量在X方向）E11（在积分点的应变分量在X方向） SF,（在积分点的截面力）*NODE FILE, NSET=ENDS U,（节点的空间位移） CF,（节点的点荷载）RF节点的反作用力）*OUTPUT, FIELD, VARIABLE=PRESELECT*ELEMENT OUTPUT SF,fig.1 应力应变曲线*OUTPUT, HISTORY*NODE OUTPUT, NSET=ENDS U, CF, RF------------------------------例二、孔平板分析A thin plate, dimensions 10cm×10cm, contains a hole of radius 1cm at its center. The plate is made from steel, which is idealized as an elastic—strain hardening plastic solid, with Young’s modulus E=210GPa and Poisson’s ratioν=0.3. The uniaxial stress—strain curve for steel is idealized as a series of straight line segments, as shown fig.1.The plate is loaded in the horizontal direction by applying tractions to its boundary.The magnitude of the loading increases linearly with time, as shown fig.2.fig.2几何及加载状况由于板内圆孔的应力集中因子大概为3，因此在时间 t＝1时 A 点将达到屈服 (板的初始屈服应力为200MPa)；在时间 t＝3时，平板都将达到屈服。

下面我们运用abaqus分析并输出时间t=1、 t=2 和t=3时板的塑。