Abaqus常用技巧总结.docx

Abaqus常用技巧总结.docx

《Abaqus常用技巧总结.docx》由会员分享，可在线阅读，更多相关《Abaqus常用技巧总结.docx（10页珍藏版）》请在冰豆网上搜索。

Abaqus常用技巧总结

Abaqus常用技巧总结

1.对timeincrement的根本理解

      abaqus的step里有maximumnumberofincrement、initialincrement、minimumincrement、maximumincrement四个量许多网友不知怎样设置合理，合理设置是建立在深刻理解基础上的。

     要理解这个问题，首先需要了解abaqus的计算过程和有限元计算收敛性问题，abaqus首先用initial值输入进行叠代计算，如果计算结果收敛，则继续以这个值代入计算下一步，如果不收敛，则自动减小时间步长（timeincrement）重新计算直到收敛然后计算下一步。

      但是如果时间步长减小到最小值minimum时计算结果还是不收敛，则abaqus将停止计算，由此可知maximum值和minimum值分别是abaqus在收敛计算时时间步长的上下限，同时totaltime=求和（timeincrement*number），当时间步长很小时，需要计算的步数number相应增大（电脑计算花的时间也随之增大），因此number一般要设置较大值。

minimum并不是越小越好，因为1）number即计算时间增大2）abaqus计算精度约在10^（-5）,当时间步长小于这个值，计算结果已经没什么意义了。

     有限元计算收敛性与（最小空间步长/时间步长）值有关，若minimum设为10^（-5）,还是不收敛，可适当减小空间步长（即把网格画细点），当然还有一些其他办法，如果实在计算不了，也许是模型本身有点问题，或改为显示explicit计算

     总而言之，maximunnumber要适当设置较大值，initial可适当改小（如-2，-3量级），minimum（-5量级）不要修改，maximum值影响不大，可不改.

2.moment的加载

      一个大筒体上有三个接管端面固定,大筒体两端加载扭距,如何加载?

（1）将大筒体两端要施加扭矩的节点分别定义为两个Nset:

left,right. 

（2）分别在大筒体两端的圆心处定义两个referencenode:

rp-left,rp-right. 

    （3）用如下命令将两个节点集绕3轴旋转的自由度与参考点耦合起来,其他自由度度是否耦合根据具体问题而定:

     *KINEMATICCOUPLING,REFNODE=rp-left

    left,6,6

    *KINEMATICCOUPLING,REFNODE=rp-right 

     right,6,6 

    （4）在两个参考点上施加绕3轴旋转的弯矩.

     提醒:

referencenode也有自由度,注意相应的边界条件.

3.abaqus计算时c盘的临时文件太大了，怎么改目录？

     临时目录是Windows自己定义的,可以在系统环境变量中修改.

4.CAE中如何加预应力

      具体没作过，看看*PRESTRESSHOLD和*INITIALCONDITIONS,TYPE=SOLUTION,REBAR

这两个命令以及ABAQUSAnalysisUser'sManual“Definingreinforcement,”Section2.2.3“Definingrebarasanelementproperty,”Section2.2.4

5.hypermesh里面看abaqus分析的结果

（1）你在abaqus中计算完成后，将结果文件输出到*.fil.

（2）利用hyperworks提供的hmabaqus.exe（在安装目录下的Altair\hw7.0\translators中）

     （3）在控制台下运行hmabaqus*.fil*res，执行完成后就生成了相应的res文件

     （4）在hyperview中打开你的模型文件*.inp和结果文件*.res，就可以查看你的结果了

6.X-YPlots

       **

      **STEP:

pre-load

      **

      *Step,name=pre-load,nlgeom

      pre-loading

      *Static

      0.01,1.,1e-05,0.1

      ........

      **

      **LOADS

      **

      **Name:

pt-loadType:

Concentratedforce

      *Cload

       _G5,2,-200.E6

      **

      ........

      **

      *Output,history,frequency=1

      *nodeoutput,nset=_G5

      CF2,U2

      *elementoutput,elset=_G5

      E22,S22

      **

      *monitor,node=_G5,dof=2

7.如何把上一次分析结果作为下一次分析的初始条件

     使用LDREAD命令，首先需要注意下面两个问题：

（1）每一个ANSYS的实体模型的面或体都要定义对应的单元类型编号材料属性.编号实常数编号单元坐标系编号这些参数在整个分析过程中保持不变而这些编号对应的属性在各个步骤中是不同的. 

（2）网格划分要满足所有步骤的要求单元类型必须兼容 

     步骤：

（1）创建实体模型

（2） 创建多个物理环境  

     设定一个物理环境中的单元类型材料属性实常数坐标系等,将这些参数的编号赋给实体模型的面或体 施加基本物理载荷和边界条件.设定求解选项:

选择一个标题使用PHYSICS,WRITE命令将物理环境存入文件中 

     （3）清楚当前的物理环境命令是PHYSICS,CLEAR 4重复第二步准备下一个物理环境 

8.材料方向与增量步

     材料方向：

     针对各向异性材料（如板金材料、复合材料等）变形体，材料方向定义材料的某一特定方向如纤维方向。

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

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

     abaqus中增量步

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

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

这时，可以通过采用多钟方法（如调整放大质量系数，单元网格优化等）调整增量大小，使计算继续进行

9.多个inp文件如何实现批处理

      写成这样：

!

nodedeform.f90

!

!

FUNCTIONS:

!

 nodedeform-Entrypointofconsoleapplication.

!

!

 Exampleofdisplaying'HelloWorld'atexecutiontime.

!

!

****************************************************************************

!

!

PROGRAM:

nodedeform

!

!

PURPOSE:

Entrypointfor'HelloWorld'sampleconsoleapplication.

!

!

****************************************************************************

 programnodedeform

USEDFLIB

!

 implicitnone

result=SYSTEMQQ（'abaqusjob=3d05ctJinteractive'）

 result=SYSTEMQQ（'abaqusjob=3d05ctJmodel-2interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-1interactive'）

!

result=SYSTEMQQ（'abaqusjob=test-2interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-3interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-4interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-5interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-6interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-7interactive'）

!

 result=SYSTEMQQ（'abaqusjob=test-8interactive'）

 endprogramnodedeform

10.关于ABAQUS的任务管理

     掌握ABAQUS的任务管理方法了，主要有两点，（系统是WindowsXP，使用的是ABAQUS6.4）

（1）ABAQUS6.4新增了任务管理的命令，可以暂停、恢复、和终止一个正在背景运行的任务，方法如下（在命令行输入并运行）：

     任务暂停：

 abaqussuspendjob=job-name

     任务恢复：

 abaqusresumejob=job-name

     任务终止：

abaqusterminatejob=job-name

     其中任务暂停（suspend）的时候，windows任务管理栏中仍会保留standard/explicit的计算线程，只是不再使用CPU资源，当任务恢复（resume）的时候继续工作。

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

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

（2）就是利用windos的任务管理栏了，使用它的processes管理功能时，不仅仅可以强行关掉一个正在运行的任务，还可以通过改变它的优先级别来改变其对CPU的占用程度。

共有Realtime，High，AboveNormal，Normal，BelowNormal，Low五个可选等级，windows默认的等级是Normal，此时所有的任务都在随机抢占CPU资源，一般ABAQUS在运行时想要运行别的程序就比较困难了，特别是一些大程序。

在觉得机器使用时有明显的延迟时，就可以把ABAQUS任务的优先级别设的低一些，就可以腾出CPU资源给别的级别高的任务了，不用机器的时候再把ABAQUS任务的级别调上去，这样就可以娱乐、工作两不误了。

 11.关于数据的输入输出

（1）输出数据到dat文件：

      *NODEPRINT,NSET=nset_name，FREQ=1

      COORD

得到的是变形前的坐标还是变形变形后的坐标？

？

（2）其实abaqus自己就带有相关的功能:

      abaqusjob=job-1suspend.可以将计算挂起.如果你需要重新进行运算

      输入abaqusjob=job-1resume.

      （3）如何输出大量节点的时间历程曲线？

      在环境文件abaqus_v6.env中添加一句：

max_history_requests=0即可。

-

ABAQUS单元小结

1、    单元表征

单元族：

单元名字里开始的字母标志着这种单元属于哪一个单元族。

C3D8I是实体单元；

S4R是壳单元；

CINPE4是无限元；

梁单元；

刚体单元；

膜单元；

特殊目的单元，例如弹簧，粘壶和质量；

桁架单元。

自由度dof（和单元族直接相关）：

每一节点处的平动和转动

1    1方向的平动

2    2方向的平动

3    3方向的平动

4    绕1轴的转动

5    绕2轴的转动

6    绕3轴的转动

7    开口截面梁单元的翘曲

8    声压或孔隙压力

9    电势

11    度（或物质扩散分析中归一化浓度）

12＋梁和壳厚度上其它点的温度

轴对称单元

1    r方向的平动

2    z方向的平动

6  r-z方向的转动

节点数：

决定单元插值的阶数

数学描述：

定义单元行为的数学理论

积分：

应用数值方法在每一单元的体积上对不同的变量进行积分。

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

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

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

单元类型的选择对模拟计算的精度和效率有重大的影响；

节点的有效自由度依赖于此节点所在的单元类型；

单元的名字完整地标明了单元族、单元的数学描述、节点数及积分类型；

所用的单元都必须指定单元性质选项。

单元性质选项不仅用来提供定义单元几何形状的附加数据，而且用来识别相关的材料性质定义；

对于实体单元，ABAQUS参考整体笛卡尔坐标系来定义单元的输出变量，如应力和应变。

可以用*ORIENTATION选项将整体坐标系改为局部坐标系；

对于三维壳单元，ABAQUS参考建立在壳表面上的一个坐标系来定义单元的输出变量。

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

2．实体单元（C）

实体单元可在其任何表面与其他单元连接起来。

C3D:

三维单元

CAX:

无扭曲轴对称单元，模拟3600的环，用于分析受轴对称载荷作用，具有轴对称几何形状的结构；

CPE:

平面应变单元，假定离面应变ε33为零，用力模拟厚结构；

CPS:

平面应力单元，假定离面应力σ33为零，用力模拟薄结构；

广义平面应变单元包括附加的推广：

离面应变可以随着模型平面内的位置线性变化。

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

可以扭曲的轴对称单元：

用来模拟初始时为轴对称的几何形状，且能沿对称轴发生扭曲。

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

反对称单元的轴对称单元：

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

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

如果不需要模拟非常大的应变或进行一个复杂的，改变接触条件的问题，则应采用二次减缩积分单元（CAX8R,CPE8R,CPS8R,C3D20R）

如果存在应力集中，则应在局部采用二次完全积分单元（CAX8,CPE8,CPS8,C3D20等）。

对含有非常大的网格扭曲模拟（大应变分析），采用细网格划分的线性减缩积分单元（CAX4R,CPE4R,CPS4R,C3D8R等）

对接触问题采用线性减缩积分单元或非协调元（CAX4I,CPE4I,CPS4I,

C3D8I）的细网格划分。

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

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

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

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

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

小结：

在实体单元中所用的数学公式和积分阶数对分析的精度和花费有显著的影响；

使用完全积分单元，尤其是一阶（线性）单元，容易形成自锁现象，正常情况不用；

一阶减缩积分单元容易出现沙漏现象；充分的单元细化可减小这种问题；

在分析中如有弯曲位移，且采用一阶减缩积分单元时，应在厚度方向至少用4个单元；

沙漏现象在二阶减缩积分单元中较少见，一般问题应考虑应用这些单元；

非协调单元的精度依赖于单元扭曲的量值；

结果的数值精度依赖于所用的网格，应进行网格细化研究以确保该网格对问题提供了唯一的解答。

但是应记住使用一个收敛网格不能保证计算结果与问题的实际行为相匹配：

它还依赖于模型其他方面的近似化和理想化程度；

通常只在想要得到精确结果的区域细划网格；

ABAQUS具有一些先进特点如子模型，它可以帮助对复杂模拟得到有用的结果。

3．壳单元（S）

可以模拟有一维尺寸（厚度）远小于另外两维尺寸，且垂直于厚度方向的应力可以忽略的结构。

一般壳单元：

S4R,S3R,SAX1,SAX2,SAX2T。

对于薄壳和厚壳问题的应用均有效，且考虑了有限薄膜应变；

薄壳单元：

STRI3,STRI35,STRI65,S4R5,S8R5,S9R5,SAXA。

强化了基尔霍夫条件，即：

垂直于壳中截面的平面保持垂直于中截面；

厚壳单元：

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/Post中画出；

壳单元利用材料方向局部化到每个单元。

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

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

单元的变量，如应力和应变，在局部方向输出；

*TRANSFORM定义节点的局部坐标系，集中载荷和边界条件被应用在局部坐标系中。

所用节点的输出，如位移，也默认为基于局部的坐标系；

矢量图可以使模拟结果可视化，特别是用来观察结构的运动和载荷路径。

4.梁单元（B）

模拟一维尺寸（长度）远大于另外二维尺寸的构件，且只有长度方向上的应力比较显著。

对于包含接触的任何模拟，应使用一阶、剪切变形的梁单元（B21,B31）

如果结构刚度非常大或者非常柔软，在几何非线性模拟中应当使用杂交梁单元（B21H,B32H,等）

使用欧拉－伯努利（三次）梁单元（B23,B33）精度很高，可模拟承受分布载荷作用的梁，例如动态振动分析。

如果横向剪切变形也很重要，要使用铁摩辛柯（二次型）梁单元（B22,B32）

模拟有开口薄壁横截面的结构，应当使用考虑了开口截面翘曲理论的梁单元（B31OS,B32OS）

小结：

梁单元的性质由截面（*BEAMSECTION或*BEAMGENERALSECTION）的数值积分决定，或直接给出截面积、惯性矩和扭转常数（*BEAMGENERALSECTION）；

当使用*BEAMGENERALSECTION选项时，模拟开始时进行一次数值积分，并且假定材料是弹性的；

ABAQUS包括大量的标准横截面形状。

其它形状可以通过定义SECTION=ARBITRARY来模拟；

必须定义横截面取向，方法是通过给出第三个节点，或者在单元性质定义中定义一个矢量。

截面取向在ABAQUS后处理中可以画出；

当梁作为壳的加强构件使用时，梁的横截面可能偏离节点；

线性和二次型包含剪切变形的影响，三次型梁不考虑剪切柔度。

开口截面梁准确地模拟了扭转和薄壁开口截面翘曲（包括翘曲约束）的影响；

多点约束和约束方程可以用来连接模型中铰接、刚性连接等节点的自由度；

“弯矩”型图使得像梁这样的一维单元的结果很清楚地表示出来；

ABAQUS后处理图的硬拷贝可以得到PostScript和HPGL的格式。

5．桁架单元（T）

只能承受拉伸和压缩载荷的杆，不能承受弯曲，模拟铰接框架结构，近似模拟线缆和弹簧。

  6．刚体单元（R）

没有独立的自由度。

  7．非线性分析

小结：

结构问题中存在着三种非线性来源：

材料、几何和边界（接触）。

这些因素的任意组合都可以出现在ABAQUS的分析中；

几何非线性发生在位移量值影响结构响应的情况下。

这包括大位移和转动效应、突然翻转和载荷硬化；

非线性问题是利用牛顿－拉弗森方法来进行迭代求解的。

非线性问题比线性问题所需要的计算机资源要高许多倍；

非线性分析步被分为许多增量步。

ABAQUS通过迭代，在新的载荷增量结束时近似地达到静力学平衡。

ABAQUS在整个模拟计算中完全控制载荷的增量和收敛性；

状态文件允许在分析运行时监控分析过程的进展。

信息文件包含了载荷增量和迭代过程的详细信息；

在每个增量步结束时可以保存计算结果，这样结构响应的演化就可以用ABAQUS/Post显示出来。

计算结果也可以用x-y图的形式绘出。

  8．材料

小结：

ABAQUS包含一个广泛的材料库，可模拟各种工程材料的性质。

其中包括金属塑性和橡胶弹性模型；

金属塑性模型的应力－应变数据必须用真实应变定义；

金属塑性模型假定材