4、*static,<adiabatic>,<direct>,<fully plastic>,<riks>,<stabilize>
数据行
(a) 此关键词的作用为声明此荷载步为静态分析,在这里同时要声明所做的分析类型是属于哪一种类型(稳定、屈曲、热等不同类型)。与ansys中的命令为antype,static相同,但此命令同时封装了时间步及时间子步的定义,即ansys中的命令;
(b) <adiabatic> 此参数设置分析为非线性热力学类型
(c) <direct> 当采用此参数时,说明使用者对荷载步的大小进行了设置,当该参数没有被定义时,程序会自动选择求解时间步长。
(d) <fully plastic> 采用此参数表明将采用塑性变形理论分析模型,即对单元进行全塑性行为监控。
(e) <riks> 此参数定义程序采用修正的riks法分析比例加载情况;
(f)<stabilize> 此参数定义程序采用自动稳定算法。
(g) 数据行的内容分两种不同情况:
**通用静力分析选项:
数据行(初始时间增量,荷载步的时间大小,允许的最小时间增量,最大时间增量)
**当关键词的选项为riks时,数据行的内容为
数据行(沿静力平衡路径的最小弧长增量,总的弧长比例系数,最小弧长增量,最大弧长增量,最大荷载增量系数,监控结点的位移限值,被监控的结点自由度,结束增量的结点整体位移值)
(h) 从上面可以看出,当采用通用静力分析选项时,数据行定义了荷载步以及荷载子步的最大与最小值;采用*static,riks关键词后表明分析模型时采用弧长法,
5、*control,analysis/parameters/reset/type,<field>
数据行
(a) 此关键词为定义求解控制选项;在ansys中有各个不同领域的求解模块,如力学模块,结构模块,流体模块,这些是在进入ansys程序中就进行选择,但是在abaqus中没有这些特定的模块,因此在此项关键词中进行区分,以确定不同领域的求解要求与符合的域平衡方程;
(b) 当关键词的必选项不同时,数据行的内容也有所不同;
(c) analysis=discontinuous 选项将使得高度非线性问题(如接触面间的滑动,混凝土的开裂等)的求解更加有效;个人的理解是当设置此参数时,程序在求解时会自动在非线性部分增加迭代次数或者增加荷载子步,使得高度非线性的求解趋于精确。
(d) parameters=field
parameters=constraints 设置约束方程的误差
parameters=line search 设置线性搜索选项
parameters=time incremention 设置时间增量选项
(e) field可选参数仅仅与parameters=field配套使用;
field=concentration 设置集中质量场平衡方程参数;
field=displacment 设置位移场与扭转角平衡方程参数;
field=electrical 设置电势场平衡方程参数;
field=global 设置一整套激活场变量平衡方程参数;
field=hydrostatic fluid pressure 设置静水流体压力平衡方程参数;
field=pore fluid pressure 设置孔隙水压力平衡方程参数;
field=rotation 设置转角平衡方程参数;
field=temperature 设置温度场平衡方程参数;
(f) reset 选项将各种设置复原,恢复到系统缺省的设置;
(h) type=direct cyclic 当设置此参数后,说明将对结构进行直接循环分析,并在分析过程中控制稳定状态;
(i) 当必选参数不同时,其可选参数以及数据行的内容也将会不同
*control,parameters=line search
数据行(搜寻到零点的最大迭代次数(缺省为0,建议值为4),最大校正系数(缺省为4),最小校正系数(缺省为0.25),线性搜索结束时残余缩减系数(缺省为0.25),新搜索步与上一搜索步校正系数之比(缺省为0.15))
*control,parameters=
**对模型施加荷载的类型有分布荷载,集中荷载等。ansys中可以对几何体进行加载,也可以对单元结点进行加载,当对几何体加载后,需要通过转换将荷载从几何体转至结点上;在abaqus中荷载只能加在结点上,而不能加载于几何实体上(不知道理解对否,有待确认),每种不同形式的荷载由不同的关键词来定义,此处仅介绍常用的两种荷载类型:均布荷载与集中荷载
6、*dsload,<amplitude>,<constant resultant>,<op>,
数据行(表面名称,荷载类型,荷载大小)
(a) 此关键词表明对模型施加表面分布荷载;ansys中采用sf,sfa等命令定义均布荷载;
(b) <amplitude>此参数定义荷载变化类型,是渐变荷载,还是跃阶荷载;如果没有定义则参照*step中的设置;
(c) <constant resultant>
(c) 当op=mod(缺省)时告诉程序保持从前定义的荷载,当op=new时以前定义的荷载将被删除,重新定义新的荷载;
7、*cload,<amplitude>,<op>,,,,,
数据行(结点号,自由度,荷载大小)
(a) 此关键词定义对模型施加集中力或者集中弯矩;ansys中采用f定义结点上的集中力;
(b) <amplitude> 指定加载曲线,当此选项省略,程序根据*step中对amplitude选项的定义对结构施加荷载。
(c) <op> 选项定义荷载施加情况:
op=mod 表明前面所施加的荷载继续存在,在本荷载步中仅对荷载进行修改或增加
op=new 表明前面对模型施加的荷载全部删除,新的荷载在本荷载步重新施加
(d) *cload的其余可选关键词如<follower>,<load case>,<region type>由于较少用到,所以在这里不进行介绍
8、*restart,read/write,<>,<>
无数据行
(a) 保存与重新使用分析结果的关键词,以下为*restart关键词在abaqus/standard分析中的内容,其在abaqus/explicit分析中的参数参见abaqus帮助文件;
(b) read/write 是必选项,两者选其一;
当为read时指明是重启动分析,此时模型的基本信息不能改动(单元,材料,结点),但是结点集与单元集可以增加。
当为write时表明重新启动数据将被记录。
(c) 可选参数根据read与write的不同而不同
*restart,read,<end step>,<inc>,<step>
*restart,write,<frequency>,<overlay>
(d) <end step> 指明使用者希望在哪一荷载停止有限元分析,当有限元分析终止后,可以重新设置荷载、输出选项、迭代精度等内容;
(e) <inc> 从启动分析中荷载子步的设置,当此项没有被定义时,采用*step中的设置;
(f) <frequency> 指明写入数据文件中的频率,frequency=2表示在荷载子步为2.4.6时将结果写入文件,缺省设置为frequency=1,如果frequency=0,则表示结果不写入文件;
(g) <overlay> 指明每一步只能保存一个荷载步数据,这样可以减小存储空间。
*****文件输出部分
文件输出部分包括对结点计算结果的输出,单元计算结果的输出,历史分析结果的输出等内容。
9、*node output,<参数>
数据行(定义需要输出的结点计算结果)
(a) 定义结点结果输出选项,且必须与*output选项匹配使用;
(b) 当与*output,field同时使用时,此项关键词的格式为 *node output,<nset>/<tracer set>
当与*output,history同时使用时,此项关键词的格式为 *node output,<nset>,<tracer set>,<variable>
(c) <nset>为结点集的名称,
(d) <tracer set>仅仅在explicit分析中定义,且只适用于输出结点位移结果;
(e) 当variable=all,程序将结点的所有变量以及材料类型均输出至结果文件中,
当variable=preselect,程序将以默认的结点结果输入至结果文件中,其余的结果将根据数据行定义的结果变量来输出;当这项参数没有定义时,结点输出变量必须在数据行进行定义。
10、*element output,<参数>
数据行(定义需要输出的单元计算结果)
(a) 定义单元结果输出选项,且必须与*output选项匹配使用;
当与*output,field同时使用时,此项关键词的格式为 *element output,<directions>,<elset>,<position>,<varialbe>
(b) directions=yes指定输出单元的材料类型,directions=no时不输出单元的材料类型。
(c) <elset>为单元集的名称;
(d) position=centroidal 是指明在单元重心点输出单元的计算结果;
position=integration points(缺省),指明在单元积分点输出计算结果;
position=node,指明通过外推法输出单元结点上的计算结果。
(e) 当variable=all,程序将结点的所有变量以及材料类型均输出至结果文件中,
当variable=preselect,程序将以默认的结点结果输入至结果文件中,其余的结果将根据数据行定义的结果变量来输出;当这项参数没有定义时,结点输出变量必须在数据行进行定义。
11、*output,diagnostics/filed/history,<>,<>,<>
(a) 定义写入输出文件中的计算结果内容,如结点计算结果,单元计算结果等,类似于ansys中的output,命令;
(b) 此关键词是一个统领关键词,可以把它想象为一个菜单,下面有子菜单,点击了菜单之后就必须点击其中的一个子菜单。子菜单的内容是确定什么计算结果需要写入结果文件中;
(c) 子菜单中常用的关键词有:
*node output:结点计算结果的提取; *element output:单元计算结果的提取; *contact output:接触结算结果的提取:
(d) 根据必选参数的不同,*output的可选参数<>也随之不同:
*output,diagnosic=yes
*output,fild,<op>,<variable>
*output,history,<frequency>,<mode list>,<op>,,<variable>
(e) diagnostics=yes(缺省)指明是否将诊断说明写入结果文件中;
filed指明计算结果是某一确定荷载步的计算结果,在ansys中为在/post26中定义变量并输出计算结果,输出的结果是某一time的计算结果;
history指明计算结果是历史记录的计算结果,在ansys中为在/中定义变量并输出计算结果,输出的内容是某一个结点在整个加载过程中的计算结果。
(f) <frequency>指定输出频率,即每几步荷载步的计算结果被输出;如果省略此选项则在每一荷载步之后输出计算结果;
<mode list> 指定输出结构的特征曲线;
<op> op=new说明在前面荷载步中定义的数据将被删除,结果文件中输入数据行定义的结果;op=add表明数据行定义的结果将追加在结果文件中;op=replace表明将代替原来的计算结果;
<variable>同*node output中的说明;
(g) 其余参数见abaqus帮助。
另外有个小疑问,*node output, *node file, *node print 三者相同的地方均为定义输出结点计算结果,那么它们的区别是什么?有点不太理解。
三者均属于History data中setp步输出文件的命令。
*node output 输出数据至.odb文件,可在view模块中观察;
*node file 输出数据至.fil文件,并不默认生成且无法在CAE中实现;
*node print 输出数据至.dat文件。
********************************
接触分析步骤与相关关键词
********************************
1、采用*surface 关键词定义接触面, 这是进行接触分析的第一步。ansys与abaqus都相同,ansys中通过cm,,node命令定义接触面的点集合;
2、采用*contact pair 关键词定义接触对,这一步是确认哪两个面是接触对,以及接触对的特性。ansys通过实常数识别接触对,abaqus通过接触名称确定接触对,但同时又对接触的一些特征进行了说明;
3、采用*surface interaction 关键此定义接触面特性
4、采用*friction 定义接触面间的摩擦系数。 ansys中通过材料将摩擦赋予接触对;
5、采用*surface behavior 定义接触面间的特征(大滑移,小滑移,有分离,无分离);ansys中接触的这些特性在接触单元的关键词中进行了设置,如对于174号单元,其keyopt(12)中定义了接触摩擦的类型;
6、在ansys中没有定义接触压应力与接触间隙的关系曲线,而在abaqus中有此方面的内容,用户可以根据需要选择线性、指数型、列表型接触压应力—间隙曲线,因此如果想分析关系复杂的接触问题,abaqus将会是一个比较有效的软件。
从上面的基本步骤可以看出,abaqus与ansys中接触的定义步骤大同小异,虽然对接触特性的定义不一样,但是其定义的内容却是一样的,从这个角度去看问题,可能会使你在学习abaqus的过程中理解与接受起来更加快。还需要说明的是在abaqus中没有专门的接触单元,接触过程是通过接触特性程序根据实际情况自动生成接触面上的各种算法来进行分析。
a、*surface,name,<combine>,<crop>,<internal>,<region type>,<trim>,<type>
数据行
(a) 此关键词定义模型中可能接触到的面,以下为在abaqus/standard分析中的内容;
(b) <combine> 定义接触面是由几个接触面共同组成的,且此关键词只能在辐射分析中使用
combine=union 指明一个或多个面组成一个整体接触界面;
combine=intersection 指出两个类型相同的面相交组成接触界面;
combine=difference 指出接触界面是由两个类型相同的面相减而得;
(c) <crop> 也是在辐射分析中使用,在此不做过多介绍;
(d) <internal> 指明界面在内部产生,默认的设置是不对这个参数进行声明;
(e) <trim> 对界面进行修正,
trim=yes(缺省设置) 指对开放的边界进行修正;
trim=no 指出不对界面进行修正。
(f) <type> 设置接触面的类型
type=element(缺省设置) 指接触界面是由单元来定义的;
type=node 定义接触界面是通过结点来定义的;
b、*contact pair,interaction,<adjust>,<extension zone>,<no thickness>,<small sliding>,<smooth>,<type>——————在abaqus/standard下的形式
从面,主面,从面切向滑动的的可能方向,主面切向滑动的的可能方向
(1) adjust 对接触面间进行调整,使得接触面闭合,这种调整不产生任何应变;此选项适用于tied接触,不适用于self-contact接触类型;
(2) extension 指出在接触分析中将主面扩大一定范围,防止从面滑出主面以外;
(3) no thickness 表明在接触计算中不考虑厚度的影响,此选项只适用于small-sliding以及surface to surface情况;
(4) small sliding 指出接触面间为小滑移情况,这个参数不能用于自接触问题
(5) smooth 表明接触面间是光滑的,缺省值为0.2,取值范围在0.0~0.5之间;
(6) type=node to surface
type=surface to surface 指定接触分析为面-面分析,并产生相应的约束系数,这种接触类型只适用于tied或者small sliding参数存在的选项;
c、*contact pair,<cpset>,<interaction>,<mechanical constraint>,<op>,<small sliding>,<weight>——————在abaqus/explicit 下的形式
第一个面,第二个面(如果两个面的名称相同,程序认为是自接触)
interaction 同surface interaction 所定义的特性相同
(1)、mechanical constraint=kinematic(缺省)选择随动接触法分析接触
mechanical constraint=penalty(缺省)选择随动接触法分析接触
(2)、op=add(缺省)在原有的接触对中加入新的接触对,
op=delet指定删除接触对
(3)、small sliding 指出接触面间为小滑动接触,而不是有限滑动,这种情况只能采用随动接触分析法;
(4)、weigth 权重系数
d、*surface interaction,name,<pad thickness>,<user>
数据行(接触面厚度或者以点为基础的接触面的截面面积,默认为1.0)
(1)、这个关键词是用来定义接触面间相互作用特性;
(2)、<pad thickness> 仅用于abaqus/explicit分析中接触对算法中;
(3)、<user> 为用户自定义接触面间的特性;
e、*surface behavior,augmented lagrange/pressure_overcloser,<no separation>
数据行
(1)、此关键词用来修改力学接触分析中默认的接触压应力—接触间隙曲线关系,且必须同关键词*surface interaction同时使用;
(2)、augmented lagrange 选项表明程序将选择增强的拉各朗日法求解接触问题;
(3)、pressure_overcloser 选项指明程序将根据命令选择相应的压应力—接触间隙曲线。应当指出的是此项关键词与<no separation>选项不能同时使用;
pressure_overcloser=hard (缺省)指明接触类型为硬接触;
pressure_overcloser=exponential 将定义压应力—接触间隙曲线为指数型曲线;
pressure_overcloser=linear 压应力—接触间隙曲线为线型关系;
pressure_overcloser=tabular 列表方式定义压应力—闭合曲线;
(4)、数据行根据参数的不同而不同
*surface behavior,augmented largrange
数据行(罚刚度(必须为正值,如果没有定义程序会根据计算来确定),接触压应力为零时接触面间的空隙(默认为零))
*surface behavior,pressure_overcloser=hard
无数据行
*surface behavior,pressure_overcloser=exponential
数据行(接触压应力为零时的接触间隙,接触间隙为零时的接触压应力)
*surface behavior,pressure_overcloser=linear
数据行(接触压应力—接触间隙曲线斜率)
*surface behavior,pressure_overcloser=tabular
数据行(接触压应力,接触间隙)
f、*friction,elastic slip/lagrange/rough/sliptolerance/user,<anisotropic>,<dependencies>,<depvar>,<exponential decay>,<properties>,<shear traction slope>,<taumax>,<test data>
(1)、此关键词定义摩擦类型,必须与*surface interaction一起使用;
(2)、elastic slip 此选项仅用于abaqus/standard分析中,定义了最大允许弹性滑移绝对值,如果这项没有定义,则弹性滑移选取slip tolerance所定义的值;
(4)、lagrange 接触间的摩擦满足拉各朗日乘法公式;
(5)、rough 指定接触为完全粗糙的,接触面间没有滑移;
(6)、slip tolerance 此选项仅用于abaqus/standard分析中,默认为0.005。此项定义了允许的最大弹性滑速度同接触面大小的比有关简单连接单元以及预应力单元的关键词
在有限元分析中常常会碰到使用连接单元的情况:接触单元分析不收敛而采用连接单元来代替;模拟两个模型间通过点对点的连接情况等。在ansys中有各种各样的弹簧单元,如最常用的39号单元,带有阻尼的14号单元等。这些单元的存在为我们分析有限元模型增添了许多灵活易变的方法。同样在abaqus中也有弹簧单元,由于这个软件所包含的内容丰富,因此对于弹簧单元将划分的更为详细。在abaqus中,模拟两者之间联系情况的单元主要有以下几种单元: spring单元,dashpot单元,joint单元以及connector等连接单元。以下将介绍这几个单元的用法与使用范围。
* spring单元
(1) 特性描述
a、可以描述力与相对位移的情况
b、可以描述弯矩与相对转角的情况
c、具有线性与非线性特性
(2) 单元类型(有三种可供选择的弹簧单元类型)
a、spring1
b、spring2
c、springA
(3)单元使用说明
a、spring1 仅用于abaqus/standard分析中,用于模拟结点与地面之间的关系。也就是说,如果选用了spring1单元,那么在你的模型中,必须有一个是固定的,被约束的;
b、spring2 仅用于abaqus/standard分析中,用于模拟某一个方向的具有弹簧特性的连接,即单元只有一个自由度;
c、springA 在abaqus/standard 与abaqus/explicit分析中都可以运用。
(4)inp文件定义形式
*element,type=spring1/spring2/springA
单元号,单元结点1,单元结点2
*spring,elset,<dependencies>,<nonlinear>,<orientation>,<rtol>
数据行
a、*element与*spring必须一起使用。在ansys中*spring相当于弹簧单元的实常数,定义弹簧单元的力与位移的关系;
b、<dependencies> 定义弹簧单元是否与场变量有关;
c、<nonlinear> 声明弹簧单元为非线性单元,缺省默认为线性单元;
d、<orientation> 声明弹簧单元约束哪一个方向的自由度;
e、<rtol> 在abaqus/explicit分析中使用
f、数据行的内容根据各个单元以及参数的变化而变化
*element,type=springA (线性)
*spring,elset=,
空行
弹簧刚度,频率,温度,第一场变量,第二场变量......
*element,type=springA (非线性)
*spring,elset=,nonlinear
空行
力,相对位移,温度,第一场变量,第二场变量......
*element,type=spring1/spring2 (线性)
*spring,elset=,
第一点的结点自由度,第二点的结点自由度
弹簧刚度,频率,温度,第一场变量,第二场变量......
*element,type=spring1/spring2 (非线性)
*spring,elset=,nonlinear
第一点的结点自由度,第二点的结点自由度
力,相对位移,温度,第一场变量,第二场变量......
*dashpot单元
(1) 特性描述
a、可以描述力与相对速度的情况
b、可以描述弯矩与相对转角变化的情况
c、具有线性与非线性特性
(2) 单元类型(有三种可供选择的弹簧单元类型)
a、dashpot1
b、dashpot2
c、dashpotA
(3)单元使用说明
a、dashpot1 仅用于abaqus/standard分析中,用于模拟结点与地面之间的关系。
b、dashpot2 仅用于abaqus/standard分析中,用于模拟某一个方向的具有缓冲特性的连接,即单元只有一个自由度;
c、dashpotA 在abaqus/standard 与abaqus/explicit分析中都可以运用。
(4)inp文件定义形式
*element,type=dashpot1/dashpot2/dashpotA
*dashpot,elset,<dependencies>,<nonlinear>,<orientation>,<rtol>
数据行
(5) *dashpot中各个参数的含义与*spring单元中参数的含义相同,两者不同的只是数据行定义的内容不同;
*element,type=dashpotA (线性)
*dashpot,elset=,
空行
单元缓冲系数,频率,温度,第一场变量,第二场变量......
*element,type=dashpotA (非线性)
*dashpot,elset=,nonlinear
空行
力,相对速度,温度,第一场变量,第二场变量......
*element,type=dashpot1/dashpot2 (线性)
*dashpot,elset=,
第一点的结点自由度,第二点的结点自由度
单元缓冲系数,频率,温度,第一场变量,第二场变量......
*element,type=dashpot1/dashpot2 (非线性)
*dashpot,elset=,nonlinear
第一点的结点自由度,第二点的结点自由度
力,相对速度,温度,第一场变量,第二场变量......
*jointc单元
(1) 特性描述
a、用来模拟连接相互作用的单元;
b、由平动弹簧以及转动弹簧(*spring单元)与阻尼器(*dashpot单元)联合组成。
(2) 单元形式
a、 jointc
(3) inp文件形式
*element,type=jointc
*joint,elset,<orientation>
*dashpot
数据行
*spring
数据行