Abaqus学习资料摘取.docx
- 文档编号:4487126
- 上传时间:2022-12-01
- 格式:DOCX
- 页数:16
- 大小:174.23KB
Abaqus学习资料摘取.docx
《Abaqus学习资料摘取.docx》由会员分享,可在线阅读,更多相关《Abaqus学习资料摘取.docx(16页珍藏版)》请在冰豆网上搜索。
Abaqus学习资料摘取
材料的塑性性能可以用它的屈服点和屈服后硬化特性来描述。
从弹性到塑性
的转变发生在材料应力-应变曲线上的某个确定点,即通常所说的弹性极限或屈
服点(见图8-2)。
屈服点的应力叫做屈服应力。
大部分金属的屈服应力为材料弹
性模量的0.05%-0.1%。
金属在到达屈服点之前的变形只产生弹性应变,在卸载后可以完全恢复。
然
而,一旦应力超过屈服应力就会开始产生永久(塑性)变形。
与这种永久变形相关
的应变称为塑性应变。
在屈服后的区间里,弹性和塑性应变共同组成了金属的变
形。
金属的刚度在材料屈服后会显著下降(见图8-2)。
已屈服的延性金属在卸载
后将恢复它的初始刚度(见图8-2)。
通常,塑性变形会提高材料重新加载时的屈
服应力:
这一特性称为硬化。
金属塑性的另一个重要特点就是非弹性变形与材料几乎不可压缩的特性相
关,这一效应给用于模拟弹-塑性的单元类型的选用带来很大的限制。
承受拉力的金属在塑性变形时,可能会在材料失效时经历局部的高度伸长与
变细,称为颈缩(见图8-2)。
金属的工程应力(力除以变形前的截面面积)称为
名义应力,与之对应的为名义应变(伸长量除以原长)。
金属在发生颈缩时的名义
应力远低于材料的极限强度。
这种材料特性是由试件几何形状、实验本身特点以
及应力应变测量方法引起的。
例如,由相同材料的压缩实验所得到的应力应变曲
线就不会出现颈缩区域,因为试件在受压变形时不会变细。
因此,描述金属塑性
的数学模型应该能够考虑拉伸和压缩的不同特性,并与结构几何形状和加载特性
无关。
为了实现这一目的,应当把已十分熟悉的名义应力的定义,FA0,和名
义应变,Δll0(这里用下标0代表材料未变形状态下的值),替换为能考虑在有
限变形中面积的改变的应力与应变的新度量方法。
有限变形中的应力应变度量
只有考虑在极限Δl→dl→0的情况下,拉伸和压缩应变是相同的,例如:
其中是l当前长度,0l是原始长度,ε为真实应变或对数应变。
与真实应变对应的应力是真实应力,定义为:
其中F是材料承受的力,A是当前面积。
承受有限变形的延性材料的真实应力—
真实应变曲线对于拉伸和压缩是一致的。
在ABAQUS中定义塑性
在ABAQUS中必须用真实应力和真实应变定义塑性。
ABAQUS需要这些值并相
应地在输入文件中解释这些数据。
大多数实验数据常常是用名义应力和名义应变的值给出的。
这时,必须应用
下列公式把塑性材料的数据从名义应力/应变的值转换为真实应力/应变的值。
先把名义应力表达为如下公式:
表达式两边同时加上1,并取自然对数可以得到真实应力和名义应力之间的关系
为:
考虑塑性变形的不可压缩性,并假定弹性变形也是不可压缩的,由体积不变
得:
当前面积和原面积的关系为:
将A的定义代入到真实应力的定义式中,得到:
其中
也可以写为1+εnom,代入上式便得到真实应力和名义应力、名义应变之
间的关系:
ABAQUS中经典金属塑性模型定义了大部分金属的屈服后特性。
ABAQUS
用连接给定数据点的一系列直线的方法来逼近材料的光滑应力-应变曲线。
因为
可以用任意多个点来逼近实际的材料性质;所以,有可能非常逼真地模拟材料的
真实性质。
在塑性数据中将材料的真实屈服应力定义为真实塑性应变的函数。
第
一组数据定义材料的初始屈服应力,因此,该组数据中的塑性应变值应该为零。
在塑性性能的材料试验数据中所提供的应变不可能仅仅是材料的塑性应变,
而应当是材料的总应变。
所以必须将总应变分解为弹性和塑性应变分量。
弹性应
变等于真实应力与杨氏模量的比值,从总应变中减去弹性应变,就得到了塑性应
变(见图8-3)。
其关系式为:
其中
:
真实塑性应变,
:
总体真实应变,
el:
真实弹性应变,
σ:
真实应力,
E:
杨氏模量。
试验数据转换为ABAQUS输入数据的例题
图8-4中的名义应力-应变曲线可以作为一个例子用来示范如何将材料塑
性特性的试验数据转换为ABAQUS适用的输入格式。
名义应力-应变曲线上的
六个点用来确定塑性数据。
第一步是用前述公式将名义应力和名义应变转化为真实应力和应变。
一旦得
到这些值,就可以用上述塑性应变与总应变和弹性应变之间的关系式来确定与每
个屈服应力值相关的塑性应变。
表8.1中给出了转换后的数据。
当应变很小时,真实值和名义值之间差别很小,而当应变很大时,二者之间
就会有明显的差别;因此,如果模拟的应变比较大,就一定要向ABAQUS提供合
适的应力-应变数据,这是极为重要的。
生成等值线图的算法需要的是节点数据,而ABAQUS/Standard
计算的却是积分点的单元变量。
ABAQUS/CAE将积分点的数据外推从而得到节
点的数据。
外推算法的阶次由单元类型决定:
对二次减缩积分单元,
ABAQUS/CAE采用线性外推算法计算单元变量的节点数值。
为了显示Von
Mises应力的等值线图,ABAQUS/CAE在每个单元内将应力分量从积分点外推
到节点位置并计算VonMises应力。
如果由节点周围各个单元计算出来的Mises
应力之间的差值落在了所定义的平均阀值(specifiedaveragingthreshold)内,
则由各个单元的应力不变量计算节点平均Mises应力。
外推过程会导致不变量
超出弹性极限。
若绘制应力张量的每一个分量的等值线图(变量S11,S22,S33,S12,
S23和S13),可以发现在固定端截面上的单元应力有明显的变化。
这将导致外
推的节点应力高于积分点应力,从而由此计算出的vonMises应力也会偏高。
积分点的Meses应力不会超出单元材料的当前屈服应力,但应力等值线中
的外推节点应力可能会超出。
另外,独立的应力分量可以超出当前的屈服应力
值;因为只要求Mises应力小于等于当前的屈服应力值。
绘制等效塑性应变图
材料的等效塑性应变(PEEQ)是用来表示材料非弹性变形的标量。
如果该变
量大于零,则材料屈服。
可以在PEEQ等值线图中确定连接环已屈服的部分。
从主菜单中选择ResultFieldOutput,在弹出的对话框中的输出变量列表中
选择PEEQ可以绘出等效应变的等值线图。
激活ContourPlotOptions对话框,
设置等效塑性应变等值线的最小值为一个很小的数值(例如:
1.E-4),那么,
ABAQUS/CAE绘制的模型中的所有深兰色区域将仍保持弹性特性(见图8-11)。
从中可以清楚的发现,连接环与母体相连的部分有明显的屈服。
等值线图
例中给出的最大塑性应变为0.1082。
当然,这里面可能包括外推过程的误差,
但这种对应变的外推误差要比对应力的小。
利用可视化查询工具检查塑性应
变最大的单元积分点处的PEEQ值。
可以发现模型中积分点处最大等效塑性应
变约为0.067,这并不一定表明外推误差很大,因为在峰值塑性变形附近的应变
梯度很大。
应变势能
对于超弹性材料,ABAQUS不用杨氏模量和泊松比,而用应变势能(U)来表
达应力-应变关系。
有多种应变势能可利用:
多项式模型、Ogden模型、
Arruda-Boyce模型、vanderWaals模型。
也可以使用较简单的多项式模型,
包括Mooney-Rivlin模型、neo-Hookean模型、简约的多项式模型和Yeoh模
型。
多项式形式的应变势能是最常用的一种形式。
可以表达为:
其中U是应变势能;Jel是弹性体积比;I1和I2是材料的变形度量;N,Cij
和Di是与温度有关的材料参数。
参数Cij描述材料的剪切特性,参数Di引入可压缩性。
如果Di设置为0,ABAQUS将认为材料为完全不可压缩,并忽略上述
公式的第二项。
如果项数N为1,初始剪切模量μ0,体积模量K0的公式如下:
如果材料是不可压缩的,则应变能密度为:
该表达式即通常所说的Moongy-Rovlin材料模型。
如果C01为0,则称neo-Hookean材料模型。
图8–17用于定义材料超弹性的各种试验的变形模式
接触分析逻辑流程图
ABAQUS在每个增量步开始之前检察所有接触相互作用状态,以判断从属
节点是脱开还是闭合。
在图11.7中p表示从属节点上的接触压力,h表示从属节
点对主面的侵入距离。
如果一个节点是闭合的,ABAQUS确定它是在滑动还是
粘结。
ABAQUS对每个闭合节点加以约束,而对那些接触状态从闭合到脱离变
化的节点撤除约束。
然后ABAQUS再次进行迭代并用计算修正值来改变模型。
在检验力或力矩的平衡前,ABAQUS先检查从属节点上接触状态的变化。
若节点在迭代后间隙变为负的或零,则它的状态由脱离变为闭合。
若节点在迭代
后接触压力变为负的,则它的状态则由闭合变为脱开。
如果检测到当前迭代步的
接触状态有变化,ABAQUS将它标识为严重不连续迭代(severediscontinuity
iteration),且不进行平衡检验。
在第一次迭代结束后,ABAQUS通过改变接触约束来反映接触状态的改变,
然后进行第二次迭代。
ABAQUS重复这个过程,直到接触状态不再变化才结束
迭代。
接着的迭代为第一次平衡迭代,并且ABAQUS进行正常的平衡收敛检查。
如果收敛检查失败,ABAQUS将进行另一次迭代。
每当一个严重不连续迭代发
生时,ABAQUS将内部平衡迭代计数器重新置零。
这个平衡迭代的计数用于确
定是否因收敛慢而放弃这个增量步。
ABAQUS重复整个过程直至获得收敛的结
果,如图11.7所示。
在信息和状态文件中,每完成一个增量步就会总结显示有多少次严重不连续
迭代,和多少次平衡迭代。
增量步的总迭代数是这两者之和。
通过区分这两类迭代,可以看到ABAQUS非常适合处理接触计算和很恰当
地完成平衡迭代。
如果严重不连续迭代数很多,而只有很少的平衡迭代,那么
ABAQUS对确定合适的接触状态就会出现困难。
在默认情况下,ABAQUS会地
放弃那些超过12个严重不连续迭代的增量步,而改用更小的增量步。
如果没有
严重不连续迭代,接触状态从一个增量步到另一个增量步之间没有改变。
可能的刚体运动取决于模型的维数。
施加约束的规则如下:
•边界条件不应施加在从面上,它会与接触约束发生冲突。
•边界条件不应施加在有线性约束方程的自由度上。
如果在有约束方程的自由度上施加边界条件,ABAQUS将在数据文件中给出错误消息。
•为了容易检查约束处的反力,边界条件不应施加在载荷作用处。
*OUTPUT:
定义输出请求
该选项用来写入接触、单元、能量、节点或诊断输出到输出数据库中。
是一个统领的关键词,下面跟随的子关键字有*CONTACTOUTPUT、*ELEMENTOUTPUT、*NODEOUTPUT等。
Standard中的输出
必需的相互排斥的参数:
DIAGNOSTICS:
设置DIAGNOSTICS=YES(默认)表明详细的诊断信息将被写入输出数据库。
NO则抑制。
FIELD:
包含该参数表明与*OUTPUT联合使用的输出请求将以场的形式写入输出数据库。
HISTORY:
包含该参数表明与*OUTPUT联合使用的输出请求将以历程的形式写入输出数据库。
可选参数:
FREQUENCY:
设置该参数等于分期的输出频率。
输出将总是在每个分析步的最后被写入输出数据库。
设置FREQUENCY=0则抑制输出。
如果该参数和NUMBERINTERVAL、TIMEINTERVAL和TIMEPOINT被省略,输出将在每个增量步写入,除了*DYNAMIC和*MODALDYNAMIC,输出会每隔10个增量步写入。
FREQUENCY、NUMBERINTERVAL、TIMEINTERVAL和TIMEPOINT是相互排斥的。
MODELIST:
包含该参数表明期望输出的一系列模态会在数据行列表显示。
该参数仅在*FREQUENCY、*COMPLEXFREQUENCY或*BUCKLE有用,并且包含FIELD。
NUMBERINTERVAL:
设置该参数等于输出数据写入时间间隔的数量。
如果该参数和FREQUENCY、TIMEINTERVAL和TIMEPOINT被省略,输出将在每个增量步写入,除了*DYNAMIC和*MODALDYNAMIC,输出会每隔10个增量步写入。
FREQUENCY、NUMBERINTERVAL、TIMEINTERVAL和TIMEPOINT是相互排斥的。
TIMEMARKS:
设置TIMEMARKS=YES(默认)将在由NUMBERINTERVAL、TIMEINTERVAL或TIMEPOINT参数控制的准确时间写入。
=NO则将在由NUMBERINTERVAL、TIMEINTERVAL或TIMEPOINT参数控制的准确时间后的增量步结束时写入。
TIMEPOINTS:
设置该参数等于*TIMEPOINTS的名字来定义输出写入的时间点。
下面的参数是可选的,而且仅当包含FIELD和HISTORY参数时有效:
OP:
设置OP=NEW(默认)表明前一个分析步定义的所有输出数据请求都将被移除。
并定义新的输出请求。
设置OP=ADD表明所定义的输出请求被添加到签名分析步的输出请求中。
设置OP=REPLACE则当前输出请求会代替相同类型(比如field)的输出请求和频率。
如果没有匹配的请求,输出请求将被看作是ADD。
TIMEINTERVAL:
设置该参数等于要写入输出状态的时间间隔。
VARIABLE:
设置VARIBALE=ALL表明所有应用到该过程的变量和材料类型将被写入输出数据库。
设置VARIABLE=PRESELECT表明当前过程的默认输出变量被写入输出数据库。
额外的输出请求可以与*OPTION其他选项联合使用得到。
如果省略该参数,则只有那些各自输出选项的输出变量请求写入输出数据库。
包含MODELIST参数的数据行,列表期望的特征模态
第一行
1、执行一系列期望的特征模态
Explicit分析中的输出
必需的、相互排斥的参数:
FIELD:
包含该参数表明输出场类型
HISTORY:
表明输出历程类型
下列参数参数是可选的并只用于FIELD:
NUMBERINTERVAL:
设置该参数等于输出数据写入时间间隔的数量。
Explicit总是写分析步开始时的结果,例如,如果NUMBERINTERVAL=10,explicit会写入11个输出数据包含分析步开始时的值。
该参数默认值必须是正整数或0。
0表示抑制所有输出。
如果该参数被省略,值是20。
TIMEMARKS:
设置TIMEMARKS=NO(默认)在增量步结束立刻写出结果到输出数据库。
=YES则在准确时间写入结果。
TIMEPOINTS:
未完设置该参数等于*TIMEPOINTS的名字来定义输出写入的时间点。
基本形式:
*OUTPUT,FIELD,VARIABLE=PRESELECT
含义:
将ABAQUS默认的场变量写入ODB文件。
*OUTPUT,HISTORY,VARIABLE=PRESELECT
含义:
将ABAQUS默认的历史变量写入ODB文件。
*PREPRINT:
设置DAT文件中记录的内容。
可选参数:
CONTACT:
(该参数只用于standard分析)设置contact=YES,则打印由接触对定义数据产生的接触约束的详细信息;NO(默认),抑制该输出。
ECHO:
设置ECHO=YES打印输入数据的响应;NO(默认)抑制该输出。
History:
设置History=YES打印历史数据;NO(默认)抑制该输出。
Model:
设置Model=YES打印模型定义数据;NO(默认)抑制该输出。
Parsubstitution:
设置parsubstitution=YES打印输入模型参数化自由原始数据的更改版本;NO(默认)抑制该输出。
Parvalues:
设置parvalues=YES打印原始输入文件的更改版本并先死用来模型参数化的参数和相应数值;NO(默认)则抑制该输出。
该选项没有数据行。
Abaqus的默认设置:
*Preprint,echo=NO,model=NO,history=NO,contact=NO。
表示:
在DAT文件中不记录对INP文件的处理过程,以及详细的模型和历史数据。
*PRE-TENSIONSECTION:
关联一个预紧节点到一个预紧截面
必需参数:
NODE:
设置该参数等于预紧节点号码或包含预紧节点的节点集(只包含一个节点)名称。
必需的、相互排斥的参数:
ELEMENT:
设置该参数等于用来定义预紧截面的杆或梁单元号码,或单元集。
仅一个单元。
SURFACE:
设置该参数等于*SURFACE定义的面的名称,用来定义预紧截面。
数据行,定义截面法向,可选的
第一行(仅一行)
1、法向的第一个分量
2、法向的第二个分量
3、法向的第三个分量
如果该数据行被省略,则standard会为连续单元计算一个平均法向;而杆或梁单元的法向是两个端点方向。
*RESTART:
保存和重用数据和分析结果
警告:
该选项可能产生非常大的数据
使用在standard中
以下参数至少一个是必需的
READ:
包含该参数表示该分析是以前分析的重启动。
基本模型定义数据(单元、材料、节点)不能在这样的重启动被改变。
但是可以添加单元集、节点集和幅值表,而且与该零件有关的、已经分析的历史数据也能被更改。
WRITE:
包含该参数表示将写入重启动数据。
使用READ参数时可选的其它参数:
ENDSTEP:
该参数表明用户希望结束设置重启动的分析中的当前分析步。
当用户希望重定义载荷历史、输出选项等,可使用该参数;如果包含该参数,数据必须包含更多的分析步以定义分析如果继续。
如果省略该参数,软件会继续分析直到完成当前定义重启动的分析步。
INC:
设置该参数等于由STEP参数指定的分析步内的增量数,其后面的分析将会恢复。
如果省略该参数,重启动将在指定的分析步的结尾开始。
ITERATION:
如果新的分析是由前一个直接循环分析重启动,设置该参数等于指定分析步内的积分数。
由于新分析只能从前一个直接循环分析的载荷循环结尾重启动,所以如果INERATION参数已经设置,则INC参数可忽略。
STEP:
设置该参数等于重启动的分析步数。
如果省略,分析将重启动最后一个可用的分析步。
使用WRITE参数时可选的其它参数:
FREQUENCY:
该参数指定重启动将写入的增量。
例如,FREQUQNCY=2,则将在增量步2、4、6写入重启动信息。
没有数据行
基本形式:
*RESTART,WRITE,FREQUENCY=0
(含义:
不输出用于重启动分析的数据。
)
*RIGIDBODY:
定义刚体及其属性
必需参数:
REFNODE:
设置该参数等于刚体参考节点的节点号或节点集。
可选参数:
ANALYTICALSURFACE:
要分配到刚性体的解析面名称
ELSET:
要分配给刚性体的单元集合名称,一个单元不能分配给多个刚体
ISOTHERMAL:
只用于热力耦合分析。
设置=YES指定等温线刚体,默认是NO。
PINNSET:
设置该参数等于包含销钉类型的节点名来分配给刚性体。
该参数可用来添加节点到刚体或重定义节点类型。
销钉类型的节点只有移动自由度与刚体关联。
POSITION:
如果参考节点由用户定义,设置POSITION=INPUT(默认)。
如果参考点位于刚体的质心,则设置POSITION=CENTER
TIENSET:
设置该参数等于包含绑定类型的节点名来分配给刚体。
该参数可用来添加节点到刚体或重定义节点类型。
绑定类型的节点的移动和转动自由度都与刚体关联。
仅当指定单元集包含刚性单元时的可选参数
DENSITY:
仅用于explicit分析。
单元集内所有刚性单元的密度
NODALTHICKNESS:
仅用于explicit分析。
表明刚性单元的厚度不从数据行读入,而是由*NODALTHICKNESS选项指定。
OFFSET:
仅用于explicit分析。
定义从单元中面到参考面的距离。
由于对刚性单元来说,没有单元级别的计算,所以指定的偏移值只影响由刚性单元构成的刚性面的接触对的操作。
该参数接受正值或负值或标记SPOS和SNEG。
正值代表单元法向方向。
当OFFSET=0.5(或SPOS),刚性单元的顶面是参考面。
当当OFFSET=-0.5(或SNEG),刚性单元的底面是参考面。
默认是OFFSET=0,表明刚性单元的中面是参考面。
Standard分析中无数据行
对于explicit中的R2D2、RB2D2、RB3D2的数据行
第一行,仅一行
1、单元的横截面面积,默认是0
对于explicit中的RAX2、RB2D2、R3D4的数据行
第一行,仅一行
1、单元的厚度,默认是0
基本格式:
*RIGIDBODY,refnode=参考点集名称,analyticalsurface=解析面集名称
*RIGIDSURFACE:
定义解析刚性面
*ROTARYINERTIA:
定义刚体转动惯量
该选项用来定义与ROTARYI单元相关的刚体转动惯量的数值。
也用在standard分析中,定义质量比例阻尼(对于直接积分动力学分析)和复合阻尼(对于模态动力学分析)。
必需参数:
ELSET:
包含ROTARYI单元的单元集
可选参数:
ALPHA:
只用于standard分析。
当用于直接积分动力学,设置该参数等于因子来为ROTARYI单元创建质量比例阻尼。
在模态分析时忽略该参数。
默认是0
COMPOSITE:
只用于standard分析。
当用于模态动力学计算复合阻尼因子时,设置该参数等于临界阻尼。
在直接积分动力学时忽略该参数。
默认是0
定义转动惯量的数据行
第一行,仅一行
1、关于局部坐标轴1的转动惯量
2、关于局部坐标轴2的转动惯量
3、关于局部坐标轴3的转动惯量
4、
5、
6、
*SECTIONCONTROLS:
指定截面控制
警告:
对于沙漏控制,使用大于默认值会产生额外的刚度响应,甚至当值太大时有时导致不稳定。
默认沙漏控制参数下出现沙漏问题表明网格太粗糙,因此,更好的解决办法是细化网格而不是施加更大的沙漏控制。
该选项用来为减缩积分单元选择非默认的沙漏控制方法,和standard中的修正的四面体或三角形单元或缩放沙漏控制的默认系数;在explicit中,也为8节点块体单元选择非默认的运动方程:
为实体和壳选择二阶方程、为实体单元激活扭曲控制、缩放线性和二次体积粘度、设置当单元破损时是否删除他们、或为上述完全破损的单元指定一标量退化参数。
等
必需参数:
NAME:
名字
可选参数:
DISTORTIONCONTROL:
只用于explicit分析。
=YES激活约束防止负体积单元出现或其他可压缩材料的过度变形,这对超弹材料是默认的。
DISTORTIONCONTROL参数对线性动力学不可用并且不能防止单元由于时间不稳定、沙漏不稳定或不切实际的物理变形造成的扭曲。
=NO不激活约束,对于除了超弹材料等单元都是默认值。
HOURGLASS:
设置HOURGLASS=COMBINED定义沙漏控制的单元粘性-刚度形式;HOURGLASS=ENHANCED基于假定的增强应变方法来控制沙漏;HOURGLASS=RELAXSTIFFNES
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Abaqus 学习 资料 摘取