abaqus程序编写.docx

该【abaqus程序编写 】是由【guoxiachuanyue013】上传分享，文档一共【21】页，该文档可以免费在线阅读，需要了解更多关于【abaqus程序编写 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。**Medel
+NODE
101,102f103,104f105,
0-.
玄」
,
,
de£inition
0-.
0-.
0-.
.
+EIEMENT,
f101,102
’102,103
#101,104
t102,104
t102,105
「103,105
’104,105
1234567
1111111
Comment
-
TYPE=T2D2f
*SOLIDSECTION
Modeldata
■-Keywordline
Datalines
Optionblock
如何写input文件
一、输入文件的组成和结构:
:,节点,单元性质,(零件可以被组装成一个有意义的模型).
历史数据的定义是模型发生了什么事情的进展,模型响应的荷载,历史被分成一系
(静态加载,动态响应等),时步的定义包括过程类型(比如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出要求.
Ataminimumthemodelconsistsofthefollowinginformation:geometry,elementsectionproperties,materialdata,loadsandboundaryconditions,analysistypeandoutputrequests.
-5,
占MATERIAL#NAME=STEEL*EIASTIC
,
**±a
占古

STATIC+BOUNTABY101HENCASTRE103,2
+CL0AD
102,2t-*NODEPRINT
Z
HF”*ELPRINT
S,
占ENDSTEP
。
1)必须有一个HEADING
Two—dimenssiona1overheadhoist;frame
SI:UnitE
1—axishoiri^~a_Kisveirtical
+PREPRINT,ECHOYES,MDDEL=YES,HISTaRY=YES
HEADING开头。
2)接下来就是模型数据部分,定义节点,单元,材料,初始条件等。模型数据的
层次为:部件,组装,模型。
必须的模型数据:
(1)几何数据:模型的几何形状是用单元和节点来定义的,结构性单元的截面是必须定义的。比如梁单元。特殊的特征也可以用特殊的单元来定义,比如弹簧单元,阻尼器,点式群体等。
(2)材料的定义:材料必须定义比如使用的是钢啊,岩石,土啊等材料可选的模型数据:
(1)零件和组合:一个模型可以用几个零件来定义有可以把几个零件组合成一个集来定义。
(2)初始条件:比如初始应力,温度,或者速度等
(3)边界条件:
(4)运动约束
(5)相互作用
(6)振幅定义
(7)输出控制
(8)环境特性
(9)用户子程序
(10)分析附属部分
3)接下来就是历史数据:定义分析的类型,荷载,输出要求等。分析的目的就是预测模型对某些外部荷载或者某些初始条件的反映。一个ABAQUS分析是建立在STEP的概念上的,(在历史数据中描述)在分析中可以定义多个STEPS。每个STEP用*STEP开始,用*ENDSTEP结束。*STEP是历史数据和模型数据的分界点,第一次出现*STEP的前面的是模型数据,后面的就是历史数据。必须的历史数据:
响应类型:必须立刻出现在*STEP选项后面。ABAQUS中有两种响应步,一种是总体分析响应步,可以是线形和非线形的,另一种是线形扰动步。
可选历史数据:
(1)荷载:通常定义某种加的荷载类型和大小。荷载可以被描述成时间的函数。
(2)边界条件输出控制
(3)辅助控制
(4)再生单元和曲面
、书写input文件的语法和规则:
:
1)必须以*开始,后面接的是选项的名字,:
*MATERIALNAME=STEEL
注释行是以**开始的.
2)如果有参数,则参数和关键词之间必须用“,”格开。
3)在参数之间必须用“,”格开。
4)关键词行中的空格可以忽略。
5)每行的长度不能超过256个字符
6)关键词和参数对大小写是不区分的。
7)参数值通常对大小写也是不区分的,但是唯一的例外是文件名区分大小写。
8)关键词和参数,以及大多数情况下的参数值是不需要全拼写出来的,只要他们之间
可以相互区分就可以了。
假如参数有响应的值,则付值号是“=”。
关键词行可以延续,比如参数的名字很长,要在下一行继续这个关键词行的话就可以这样做,用“,”来连接。比如:
*ELASTIC,TYPE=ISOTROPIC,
DEPENDENCIES=1
有些选项允许INPUT和FILE的参数作为一个输入文件名,这样的文件名必须包括一个完整的路径名或者是一个相对路径名。
:(数据行如果和关键词相联系必须紧跟关键词行)
一个数据行包括空格在内不能超过256个字符。
所有的数据条目之间必须用“,”格开。
一行中必须包括指定说明的数据条目的数字。
每行结尾的空数据域可以省略。
浮点数最多可以占用20个字符。
整数可以是10个
字符串可以是80个
延续行可以被用到特定的情况。
:什么是标签呢,比如集名,曲面名,钢筋名,他们是区分大小写的,长度可以有80个字符长。标签中的空格是可以省略的,除非用“”来标示。那就不能省略了。没有用“”来标示的标签必须用字母来开头。如果一个标签用“”来定义那么“”也是标签的一部分。标签的开始和结束不能用双重“_”下面是一个使用了引号和没有使用引号的例子:
*ELEMENT,TYPE=SPRINGA,ELSET="Oneelement"
1,1,2
*SPRING,ELSET="Oneelement"
-5,
*NSET,ELSET="Oneelement",NSET=NODESET
*BOUNDARY
nodeset,1,2
:数据行可以重复,也就是说每行数据可以有一行响应的变量也可以有几行。同样也可以有多行数据行,:
*ELASTIC,TYPE=ISOTROPIC
,,20.
,,400.
,,,均质,线弹性,在不同应力下的扬式模量和泊松比.
三、从外存储器中引入模型或者历史数据
The*,commentlines,andother*,-of-fileisreached,ABAQUSwillreturntoprocessingtheoriginalfile.
Amaximumoffivelevelsofnested*.
Usage:*INCLUDE,INPUT=file_name
四、举例:下面介绍两个简单的例子:例一:悬臂梁.
输入文件的开始就是文件头,以HEADING开始,随后是模型的名字,如下所示:
*HEADING
CANTILEVERBEAM然后是网格定义:,网格包括(单元和节点)假如我们的悬臂梁有五个单元,六个节点,下面我们首先详细说明节点:
*NODE,NSET=ENDS
1,0.
6,100.
*NGEN
1,6
节点组集,NSET其值(名字),第一个节点是从0开始的,第六个节点是在100结束的.
同样我们来定义单元:
*ELEMENT,TYPE=B21(单元类型)
1,1,2(单元类型的参数)
*ELGEN,ELSET=BEAM(产生单元集,及其名称)
1,5(一个单元集,包括5个单元)现在定义单元的性质:
*BEAMSECTION,SECTION=RECTANGULAR,ELSET=BEAM,MATERIAL=STEEL
1.,,截面的形状是矩形,单元集的名称是梁单元,材料是钢。截面的尺寸是1*2。下面定义材料的性质:
*MATERIAL,NAME=STEEL
*ELASTIC
,
材料是钢,弹性,弹性模量是30E6。
下面定义边界:
*BOUNDARY
6,ENCASTRE
边界是在6节点,通过ENCASTRE来描述。
边界也可以用下面的形式来定义:
*BOUNDARY
6,1,6
ABAQUS对结构单元的中的节点的自由度使用常规的编号方式。1,2,3代表的是位移分量;4,5,6代表的是旋转分量。
以上是模型数据的定义,下面开始历史数据的定义:(加载的次序,事件的发生,还有我们想看到的变量的响应)
时步(步骤)的定义。
*STEP,PERTURBATION(步骤的开始,扰动是其名称)
*STATIC(静态分析)
*CLOAD(集中荷载)
1,2,-20000.(在节点1,y[2代表的是y方向,假如是1则代表的是x方向,3代表的是z方向]方向施加荷载,荷载的大小是-20000。也就是向下施加荷载。)
*ENDSTEP(步骤的结束)下面来解释输出要求:
*ELPRINT,POSITION=AVERAGEDATNODES,(节点的平均值)SUMMARY=YES在表的下部求和)
S11,(积分点的应力分量在X方向)E11(在积分点的应变分量在X方向)
SF,(在积分点的截面力)
*NODEFILE,NSET=ENDS
U,(节点的空间位移)CF,(节点的点荷载)RF节点的反作用力)*OUTPUT,FIELD,VARIABLE=PRESELECT*ELEMENTOUTPUT

andPoisson'sratiov=
SF,*OUTPUT,HISTORY*NODEOUTPUT,NSET=ENDS
U,CF,RF例二、孔平板分析
Athinplate,dimensionslOcmxlOcm,,whichisidealizedasanelastic—strainhardeningplastic
solid,withYoung'smodulusE=210GPastress—straincurveforsteelisidealizedasaseriesofstraightlinesegments,.
,,所以在时间t=l时A点将达到屈服(板的初始屈服应力为200MPa);在时间t=3时,平板都将达到屈服。下面我们利用abaqus分析并输出时间t=l、t=2和t=3时板的塑性变化情况。
注意到板及载荷的对称性,只需要考虑板的1/4,并在底部和边缘施加对称边界条件。采用4节点二次平面应力单元。。
5cm
ABAQUS
荷的分布

件的施加及
input文件祥解:
*HEADING
STRESSANALYSISFORAPLATEWITHAHOLE**文件名*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES**指定什么输出到NGEN,LINE=C,NSET=HOLE
.dat文件中
*RESTART,WRITE,FREQ=1**,提供给后处理*FILEFORMAT,ZEROINCREMENT**,后处理用,这里要求输出分析前
的数据(即载荷施加前)
**下面产生网格(节点、单元)
**节点的产生
*NODE101,,
119,-02,
1919,,-02
131,-02,
1031,-02,-02
1931,,-02
**
**
**
119,1919,100,101
**
**
**
*NGEN,NSET=OUTER
131,1031,100
**
**
**
**将节点加到节点集outer中
*NGEN,NSET=OUTER
1031,1931,100
**
**
**
**沿着径向线填充
*NFILL,NSET=PLATE
HOLE,OUTER,12,1
**
**
**
**单元的产生
**
*ELEMENT,TYPE=CPS4
19,119,120,220,219
**
**
*ELGEN,ELSET=PLATE
19,12,1,1,18,100,100
**
**
****
331
231
131
1031
1919
1931
1930
1929
1928
119120>129130131
31
231
219H
119
-Row1
227128129130131
**给单元赋予属性(即单元由什么材料组成),将名为steel的材料赋给单元集plate
*SOLIDSECTION,MATERIAL=STEEL,ELSET=PLATE
**
**定义材料性质
**
*MATERIAL,NAME=STEEL
*ELASTIC
,
*PLASTIC
,
,
,
,
,
,
**
**
**指定与时间无关的边界条件(任何与时间无关的边界条件都必须在分析步前指定)
*NSET,NSET=BOTTOM,GENERATE
119,131,1
*NSET,NSET=LEFT,GENERATE
1919,1931,1
***BOUNDARY
BOTTOM,YSYMM
LEFT,XSYMM
****定义载荷施加的位置和随时间的变化
**definingloads
*ELSET,ELSET=EDGE,GENERATE
30,830,100
***AMPLITUDE,NAME=HIST,TIME=TOTALTIME
,,,,,,,
stress,MPa
200
166
83
2
Step1
El23456781
i1234^uicrementsS
Uicrementjfsec
►
Step2>
Step3
**这里的关键词AMPLITUDE为整个分析中的加载历史
**定义载荷步
**
**LoadStepDefinition
**FirstLoadStep
*STEP,AMPLITUDE=RAMP
**关键词step表示开始分析,这里的关键词AMPLITUDE表示在这步中如何施加载荷,
AMPLITUDE=RAMP为平滑加载,AMPLITUDE=STEP为立即加载*STATIC**,**前面的数表示估计在这步中计算变形的初始时间增量,由于可以预计在这步变形为弹性变形,ABAQUS可以在这步结束时直接得到结果,不必更小的时间增量步,所以设定时间增量等于步长。后面的数指定这个载荷步的时间间隔,这里分析开始时t=0,结束时t=1,故间隔为1。
*DL0AD,AMPLITUDE=HIST**关键词DLOAD表示压力,其值为正,当为拉力时为负
EDGE,P2,-**施加82MPA的拉力在单元集EDGE的面2上
**为后处理输出过程数据文件*ELFILE,P0SITI0N=AVERAGEDATNODES**其中POSITION=AVERAGEDATNODES目的是光滑曲线
S,E