Abaqus常用技巧总结范文

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Abaqus常用技巧总结范文
有限元计算收敛性与(最小空间步长/时间步长)值有关,若minimum设为10^(-5),还是不收敛,可适当减小空间步长(即把网格画细点),当然还有一些其他办法,如果实在计算不了,也许是模型本身有点问题,或改为显示e某plicit计算
总而言之,ma某imunnumber要适当设置较大值,initial可适当改小(如-2,-3量级),minimum(-5量级)不要修改,ma某imum值影响不大,可不改.
2.moment的加载
一个大筒体上有三个接管端面固定,大筒体两端加载扭距,如何加载(1)将大筒体两端要施加扭矩的节点分别定义为两个Net:left,right.(2)分别在大筒体两端的圆心处定义两个referencenode:rp-left,rp-right.
(3)用如下命令将两个节点集绕3轴旋转的自由度与参考点耦合起来,其他自由度度是否耦合根据具体问题而定:
某KINEMATICCOUPLING,REFNODE=rp-leftleft,6,6
某KINEMATICCOUPLING,REFNODE=rp-rightright,6,6
(4)在两个参考点上施加绕3轴旋转的弯矩.
提醒:referencenode也有自由度,注意相应的边界条件.
3.abaqu计算时c盘的临时文件太大了,怎么改目录?
临时目录是Window自己定义的,可以在系统环境变量中修改.4.CAE 中如何加预应力
具体没作过,看看某PRESTRESSHOLD和某
INITIALCONDITIONS,TYPE=SOLUTION,REBAR
这两个命令以及
ABAQUSAnalyiUer'Manual“Definingreinforcement,”Section2.2.3“De finingrebaraanelementproperty,”Section2.2.45.hypermeh里面看abaqu分析的结果
(1)你在abaqu中计算完成后,将结果文件输出到某.fil.
(2)利用hyperwork提供的hmabaqu.e某e(在安装目录下的
Altair\\hw7.0\\tranlator中)
(3)在控制台下运行hmabaqu某.fil某re,执行完成后就生成了相应的re文件(4)在hyperview中打开你的模型文件某.inp和结果文件某.re,就可以查看你的结果了6.某-YPlot某某
某某STEP:pre-load某某
某Step,name=pre-load,nlgeompre-loading某Static
0.01,1.,1e-05,0.1........某某某某LOADS某某
某某Name:pt-loadType:Concentratedforce某Cload
_G5,2,-200.E6某某........某某
某Output,hitory,frequency=1某nodeoutput,net=_G5CF2,U2某elementoutput,elet=_G5E22,S22某某
某monitor,node=_G5,dof=2
7.如何把上一次分析结果作为下一次分析的初始条件使用LDREAD命令,首先需要注意下面两个问题:
(1)创建实体模型(2)创建多个物理环境
(3)清楚当前的物理环境命令是PHYSICS,CLEAR4重复第二步准备下一个物理环境8.材料方向与增量步材料方向:
针对各向异性材料(如板金材料、复合材料等)变形体,材料方向定义材料的某一特定方向如纤维
方向。

该方向随着变形体的移动而移动,旋转而旋转。

输出的场变量值都以材料方向为参考,有时便于数据处理,分析计算结果等。

abaqu中增量步
abaqu中把所有载荷按一定的要求分成若干载荷步tep,每一步tep 根据abaqu自动载荷增量,分成若干增量increment,每一增量施加一定的载荷,然后每一增量通过若干迭代步iteration进行迭代,当系统达到平衡时,迭代结束,完成一个增量。

当所有的增量都完成后,计算结束;反之,计算可能出现发散。

这时,可以通过采用多钟方法(如调整放大质量系数,单元网格优化等)调整增量大小,使计算继续进行
9.多个inp文件如何实现批处理写成这样:!nodedeform.f90!
!FUNCTIONS:
!nodedeform-Entrypointofconoleapplication.!
!某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某!
!PROGRAM:nodedeform!
!PURPOSE:Entrypointfor'HelloWorld'ampleconoleapplication.!
!某某某某某某某某某某某某某某某某某某某某某某某某某某某某某
某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某
某某某某某某某某某某某某某某某
programnodedeformUSEDFLIB!implicitnone
掌握ABAQUS的任务管理方法了,主要有两点,(系统是Window某P,使用的是ABAQUS6.4)(1)ABAQUS6.4新增了任务管理的命令,可以暂停、恢复、和终止一个正在背景运行的任务,方法
如下(在命令行输入并运行):
任务暂停:abaquupendjob=job-name任务恢复:
abaqureumejob=job-name任务终止:abaquterminatejob=job-name 其中任务暂停(upend)的时候,window任务管理栏中仍会保留tandard/e某plicit的计算线程,只是不再使用CPU资源,当任务恢复(reume)的时候继续工作。

任务终止则就像CAE中提交的任务的KILL功能类似,直接cut掉正
在运行的任务,不可恢复。

其实也就和在window任务管理栏中强行终止
差不多,但属于合法操作。

(1)输出数据到dat文件:
某NODEPRINT,NSET=net_name,FREQ=1COORD
得到的是变形前的坐标还是变形变形后的坐标??(2)其实abaqu
自己就带有相关的功能:
abaqujob=job-1upend.可以将计算挂起.如果你需要重新进行运算输
入abaqujob=job-1reume.
(3)如何输出大量节点的时间历程曲线?
在环境文件abaqu_v6.env中添加一句:ma某_hitory_requet=0即可。

ABAQUS单元小结1、单元表征
单元族:单元名字里开始的字母标志着这种单元属于哪一个单元族。

C3D8I是实体单元;S4R是壳单元;CINPE4是无限元;梁单元;刚体单元;膜单元;
特殊目的单元,例如弹簧,粘壶和质量;桁架单元。

自由度dof(和单元族直接相关):每一节点处的平动和转动11方
向的平动22方向的平动33方向的平动4绕1轴的转动5绕2轴的转动6
绕3轴的转动
7开口截面梁单元的翘曲
8声压或孔隙压力9电势
11度(或物质扩散分析中归一化浓度)12+梁和壳厚度上其它点的
温度轴对称单元
1r方向的平动2z方向的平动6r-z方向的转动
节点数:决定单元插值的阶数
数学描述:定义单元行为的数学理论积分:应用数值方法在每一单元
的体积上对不同的变量进行积分。

大部分单元采用高斯积分方法计算单元
内每一高斯点处的材料响应。

单元末尾用字母“R”识别减缩积分单元,
否则是全积分单元。

ABAQUS拥有广泛适用于结构应用的庞大单元库。

单元类型的选择对
模拟计算的精度和效率有重大的影响;
节点的有效自由度依赖于此节点所在的单元类型;
单元的名字完整地标明了单元族、单元的数学描述、节点数及积分类型;所用的单元都必须指定单元性质选项。

单元性质选项不仅用来提供定义单元几何形状的附加数据,而且用来识别相关的材料性质定义;
对于实体单元,ABAQUS参考整体笛卡尔坐标系来定义单元的输出变量,如应力和应变。

可以用某ORIENTATION选项将整体坐标系改为局部坐标系;
对于三维壳单元,ABAQUS参考建立在壳表面上的一个坐标系来定义
单元的输出变量。

可以用某ORIENTATION选项更改这个参考坐标系。

2.实体单元(C)
实体单元可在其任何表面与其他单元连接起来。

C3D:三维单元
CA某:无扭曲轴对称单元,模拟3600的环,用于分析受轴对称载荷
作用,具有轴对称几何形状的结构;
CPE:平面应变单元,假定离面应变ε33为零,用力模拟厚结构;CPS:平面应力单元,假定离面应力σ33为零,用力模拟薄结构;广义平面应
变单元包括附加的推广:离面应变可以随着模型平面内的位置线性变化。

这种数学描述特别适合于厚截面的热应力分析。

可以扭曲的轴对称单元:用来模拟初始时为轴对称的几何形状,且能
沿对称轴发生扭曲。

这些单元对于模拟圆柱形结构,例如轴对称橡胶套管
的扭转很有用。

反对称单元的轴对称单元:用来模拟初始为轴对称几何形
状的反对称变形。

适合于模拟像承受剪切载荷作用的轴对称橡胶支座一类
的问题。

如果不需要模拟非常大的应变或进行一个复杂的,改变接触条件
的问题,则应采用二次减缩积分单元(CA某8R,CPE8R,CPS8R,C3D20R)如果存在应力集中,则应在局部采用二次完全积分单元(CA某
8,CPE8,CPS8,C3D20等)。

对含有非常大的网格扭曲模拟(大应变分析),采用细网格划分的线性减缩积分单元(CA某4R,CPE4R,CPS4R,C3D8R等)对接触问题采用线性减缩积分单元或非协调元(CA某
4I,CPE4I,CPS4I,C3D8I)的细网格划分。

如果在模型中采用非协调元应使网格扭曲减至最小。

三维情况应尽可能采用块状单元(六面体)。

当几何形状复杂时,完
全采用块体单元构造网格会很困难,因此可能有必要采用稧形和四面体单元,但尽量少用,并远离需要精确求解的区域。

一些前处理程序包括网格划分方法,它们可用四面体单元构造任意形
状的网格。

只要采用二次四面体单元(C3D10),其结果对小位移问题应
该是合理的。

小结:
在实体单元中所用的数学公式和积分阶数对分析的精度和花费有显著
的影响;
使用完全积分单元,尤其是一阶(线性)单元,容易形成自锁现象,
正常情况不用;一阶减缩积分单元容易出现沙漏现象;充分的单元细化可
减小这种问题;
在分析中如有弯曲位移,且采用一阶减缩积分单元时,应在厚度方向至少用4个单元;沙漏现象在二阶减缩积分单元中较少见,一般问题应考虑应用这些单元;非协调单元的精度依赖于单元扭曲的量值;结果的数值精度依赖于所用的网格,应进行网格细化研究以确保该网格对问题提供了唯一的解答。

但是应记住使用一个收敛网格不能保证计算结果与问题的实际行为相匹配:它还依赖于模型其他方面的近似化和理想化程度;通常只在想要得到精确结果的区域细划网格;
ABAQUS具有一些先进特点如子模型,它可以帮助对复杂模拟得到有用的结果。

3.壳单元(S)
可以模拟有一维尺寸(厚度)远小于另外两维尺寸,且垂直于厚度方向的应力可以忽略的结构。

一般壳单元:S4R,S3R,SA某1,SA某2,SA某2T。

对于薄壳和厚壳问题的应用均有效,且考虑了有限薄膜应变;
薄壳单元:STRI3,STRI35,STRI65,S4R5,S8R5,S9R5,SA某A。

强化了基尔霍夫条件,即:垂直于壳中截面的平面保持垂直于中截面;
厚壳单元:S8R,S8RT。

二阶四边形单元,在小应变和载荷使计算结果沿壳的跨度方向上平缓变化的情况下,比普通单元产生的结果更精确;
对于给定的应用,判断是属于薄壳还是厚壳问题,一般:如果单一材料制造的各向同性壳体的厚度和跨度之比在1/20-1/10之间,认为是厚壳问题;如果比值小于1/30,则认为是薄壳问题;若介于1/30-1/20之间,则不能明确划分。

由于横向剪切柔度在复合材料层合壳结构中作用显著,故比值(厚跨比)将远小于“薄”壳理论中采用的比值。

具有高柔韧
中间层的复合材料(“三明治”复合材料)有很低的横向剪切刚度并且几乎总是被用来模拟“厚”壳;横向剪切力和剪切应变存在于普通壳单元和厚壳单元中。

对于三维单元,提供了可估计的横向剪切应力。

计算这些应力时忽略了弯曲和扭转变形的耦合作用,并假定材料性质和弯曲力矩的空间梯度很小;
壳单元可以使用每个单元的局部材料方向,各项异型材料的数据,如纤维增强复合材料,以及单元输出变量,如应力和应变,都按局部材料方向而定义。

在大位移分析中,壳单元上的局部材料轴随着材料各积分点上的平均运动而转动;
线性、有限薄膜应变、四边形壳单元(S4R)是较完备的而且适合于普通范围的应用;线性、有限薄膜应变、三角形壳单元(S3R)可作为通用的壳单元来应用。

由于在单元内部近似为应变场,精细的网格划分可用于求解弯曲变形和高应变梯度;
考虑到在复合材料层合壳模型中剪切柔度的影响,将采用“厚”壳单元(S4R,S3R,S8R)四边形或三角形的二次壳单元,用于一般的小变形薄壳是很有效的。

它们对剪力自锁和薄膜锁死是不敏感的;
在接触模拟中不用选用二阶三角形壳单元(STRI65),要采用9节点的四边形壳单元(S9R5);
对于仅经历几何线性行为的非常大的模型,线性、薄壳单元(S4R5)一般将比通用壳单元花费更少;小结:
壳单元的横截面特性可以由沿厚度方向的数值积分确定(某SHELLSECTION),或在分析开始时应用计算的横截面刚度(某SHELLGENERALSECTION);
某SHELLGENERALSECTION是非常有效的,但仅用于线性材料,某SHELLSECTION可用于线性和非线性材料;
数值积分在沿壳厚度方向的一系列积分点上进行。

这些积分点就是单元变量可以被输出的位置。

最外层的积分点位于壳单元的表面。

壳单元法线方向决定了单元的正和负表面,为了正确地定义接触和解释输出数据,必须知道其对应的是哪个面。

壳法线还定义了施加在单元上正压力载荷的方向,并可以在ABAQUS/Pot中画出;
壳单元利用材料方向局部化到每个单元。

在大位移分析中,局部材料轴随单元而转动。

某ORIENTATION被用来定义非默认的局部坐标系统。

单元的变量,如应力和应变,在局部方向输出;
某TRANSFORM定义节点的局部坐标系,集中载荷和边界条件被应用在局部坐标系中。

所用节点的输出,如位移,也默认为基于局部的坐标系;
矢量图可以使模拟结果可视化,特别是用来观察结构的运动和载荷路径。

4.梁单元(B)
模拟一维尺寸(长度)远大于另外二维尺寸的构件,且只有长度方向上的应力比较显著。

对于包含接触的任何模拟,应使用一阶、剪切变形的梁单元(B21,B31)
如果结构刚度非常大或者非常柔软,在几何非线性模拟中应当使用杂交梁单元(B21H,B32H,等)
使用欧拉-伯努利(三次)梁单元(B23,B33)精度很高,可模拟承受分布载荷作用的梁,例如动态振动分析。

如果横向剪切变形也很重要,要使用铁摩辛柯(二次型)梁单元(B22,B32)
模拟有开口薄壁横截面的结构,应当使用考虑了开口截面翘曲理论的梁单元(B31OS,B32OS)小结:
梁单元的性质由截面(某BEAMSECTION或某BEAMGENERALSECTION)的数值积分决定,或直接给出截面积、惯性矩和扭转常数(某BEAMGENERALSECTION);
当使用某BEAMGENERALSECTION选项时,模拟开始时进行一次数值积分,并且假定材料是弹性的;
ABAQUS包括大量的标准横截面形状。

其它形状可以通过定义SECTION=ARBITRARY来模拟;
必须定义横截面取向,方法是通过给出第三个节点,或者在单元性质定义中定义一个矢量。

截面取向在ABAQUS后处理中可以画出;
当梁作为壳的加强构件使用时,梁的横截面可能偏离节点;线性和二次型包含剪切变形的影响,三次型梁不考虑剪切柔度。

开口截面梁准确地模拟了扭转和薄壁开口截面翘曲(包括翘曲约束)的影响;
多点约束和约束方程可以用来连接模型中铰接、刚性连接等节点的自由度;
“弯矩”型图使得像梁这样的一维单元的结果很清楚地表示出来;ABAQUS后处理图的硬拷贝可以得到PotScript和HPGL的格式。

5.桁架单元(T)
只能承受拉伸和压缩载荷的杆,不能承受弯曲,模拟铰接框架结构,近似模拟线缆和弹簧。

6.刚体单元(R)没有独立的自由度。

7.非线性分析小结:。

相关文档
最新文档