abaqus技巧汇总1

abaqus技巧汇总1
abaqus技巧汇总1

浙江大学bbs

采用abaqus的cae进行力学问题的分析,其对模型的处理存在很多的技巧,对abaqus的一些分析技巧进行一些概述,希望对大家有所帮助。

1.abaqus的多图层绘图

abaqus的cae默认一个视区仅仅绘出一个图形,譬如contor图,变形图,x-y 曲线图等,其实在abaqus里面存在一个类似于origin里面的图层的概念,对于每个当前视区里面的图形都可以建立一个图层,并且可以将多个图层合并在一个图形里面,称之为Overlay Plot,譬如你可以在同一副图中,左边绘出contor图,右边绘出x-y图等等,并且在abaqus里面的操作也是很简单的。

1.首先进入可视化模块,当然要先打开你的模型数据文件(。odb)

2.第一步要先创建好你的图形,譬如变形图等等

3.进入view里面的overlay plot,点击creat,创建一个图层,现在在viewport layer 里出现了你创建的图层了

4.注意你创建的图层,可以看到在visible 下面有个选择的标记,表示在视区里面你的图层是否可见,和autocad里面是一样,取消则不可见current表示是否是当前图层,有些操作只能对当前图层操作有效,同cad name是你建立图层的名称,其他的属性值和你的模型数据库及图形的类型有关,一般不能改动的。

5.重复2-4步就可以创建多个图层了

6.创建好之后就可以选择plot/apply,则在视区显示出所有的可见的图层

1.什么是子结构

子结构也叫超单元的(两者还是有点区别的,文后会谈到),子结构并不是abaqus 里面的新东东,而是有限元里面的一个概念,所谓子结构就是将一组单元组合为一个单元(称为超单元),注意是一个单元,这个单元和你用的其他任何一种类型的单元一样使用。

2.为什么要用子结构

使用子结构并不是为了好玩,凡是建过大型有限元模型的兄弟们都可能碰到过计算一个问题要花几个小时,一两天甚至由于单元太多无法求解的情况,子结构正是针对这类问题的一种解决方法,所以子结构肯定是对一个大型的有限元模型的,譬如在求解非线性问题的时候,因为对于一个非线性问题,系统往往经过多次迭代,每次这个系统的刚度矩阵都会被重新计算,而一般来说一个大型问题往往有很大一部分的变形是很小的,把这部分作为一个子结构,其刚度矩阵仅要计算一次,大大节约了计算时间。

3.哪些情况可以使用子结构

前面提到的非线性问题,包括了很小变形的或者线弹性部分可以使用子结构,特别是当模型中有很多相同的部分时,提到的最多的一个例子就是桌子的四条腿,四条腿作为子结构(因为基本时弹性变形)可以包括了很多的实体单元,可以大大提高效率再一个就是问题确实太大,只有采用子结构将问题分成很多块,计算出结果后再次采用子结构分块计算,一直到能对每块单独计算为止。

4.abaqus中子结构的特点及要注意的问题

子结构是一组单元的集合,但是在子结构中仅仅只有你指定的那些节点的自由度会保留下来而其他节点的自由度都被消除了,其他的节点均是通过线性插值的方式获得求解;

子结构是通过你指定的节点与其他的单元建立联系的;在abaqus的6.4版本中只有子结构这个概念没有超单元了,其区别就是子结构可以求得单元内部准确的解而超单元不行;当你定义子结构的时候不要包含太多的单元,因为单元的刚度矩阵集成的时候会花掉太多的时间,可以用更多的含有较少单元的子结构代替。5.在abaqus中子结构的用法

一般包括如下部分,首先你要利用substructure generate和retained nodal dofs定义一个子结构,然后你可以定义子结构的内部荷载,边界条件一旦你已经定义了子结构以后你就可以象利用一般的单元一样使用子结构了,譬如输出请求等等。因为内容太多只能另外的文章再述了,大家也可以直接看看手册,要注意的是,abaqus cae是不支持子结构的。

如何在不同的分析步改变材料的参数

我所了解的大概有三种不同的方法:

1.最强大的当然是采用umat的方式,不过需要有深厚的有限元基础,一般人不推荐使用

2.采用场变量,不过功能相对简单

3.采用abaqus的import命令将前面分析的结果传递到新的分析之中

这里介绍下第二种方法

*什么是场变量

所谓场变量,我的理解就是一个环境变量,它建立了一个与材料参数之间的中介,虽然不能直接指定材料参数在不同的分析步具有不同的值,但是通过场变量,间接的达到了目的。

*怎样使用场变量

其实场变量用的较多的实在热力学和流体力学的分析种,这里介绍的仅仅是在固体力学中的用法

1.定义场变量

*你可以在initial中指定场变量的值,格式如下

initial conditions,type=field,variable=n(场变量的编号)

Set-1(你定义的结点集),1.0(场变量的值)

场变量是通过编号来识别的,一次只能定义一个场变量

*你也可以直接在分析步中指定场变量的值,格式如下

*field, variable=1

Set-1,1

当然也可以同时使用initial和field,当你指定的场变量改变时,默认材料的参数是在增量步间线性变化的。

2.建立材料参数和场变量之间的联系

如果你用cae,在prop模块里面的材料参数一般都有Number of field variables,场变量都是从1开始的,你也可以选择多个场变量。填入场变量的值和材料参数间的关系,譬如

杨是模量泊松比field1

200.E9, 0.3, 1.

180.E9, 0.3, 2.

如果你用的是命令格式,则在inp文件里面键入:

*ELASTIC, DEPENDENCIES=1

200.E9, 0.3, , 1.

180.E9, 0.3, , 2.

第四个参数表示场变量的值

3.注意

场变量在不同分析步中的值有你在不同的分析步中指定,如果没有指定,材料参数默认微是场变量1的值,例如

*STEP,name=step1

*STATIC......

*FIELD, V ARIABLE=1

NALL, 1.

*END STEP

*STEP,name=step2

*STATIC......

*FIELD, V ARIABLE=1

NALL, 2.

*END STEP

cae步支持场变量,所以你必须自己更改inp文件,

4.技巧

如果材料的参数变化比较复杂,一般是利用幅值曲线来定义场变量值的变化,*FIELD, V ARIABLE=1,amplitude=???

对幅值曲线步清楚的自己可以看手册

*模型的重启动分析-restart

按理说restart不应该算是一个分析的技巧,而是一个常识,不过呢可能有很多朋友没有建过大型模型导致restart也用的较少,所以也介绍下

1.什么是restart

你的job可能包含多个step,可是如果你的模型很大,可能会有这样一种情况,当你花了几天几夜,终于分析好的时候,你发现the first step的边界条件设置的有问题,这对于你真是晴天霹雳,于是你只好重新来过,可是低二天你发现你的电脑restart,这时的你可能只能问上帝了,how can i do?

*restart,就是将一个复杂的模型分析过程分成很多的阶段,甚至是一个increatment step一个阶段,你可以对每个阶段的结果进行检验,然后进入下一个阶段进行分析。

2.重启动需要那些文件

对于standard来说,.res,.mdl,.stt,.prt,.odb,这些文件是用于重启动的,explict 是.abq,.stt,.prt,.odb.

3.如何在一个分析中设置重启动来生成以上文件。

这里只介绍下在standard的用法,其实很简单?

inp文件里面加入*RESTART, WRITE, FREQUENCY=N就可以了

cae默认加入了重启选项,不过可以在step->output->restart request里面设置输出的频率,也就是frequency。

*技巧:因为res文件包含了模型的几乎全部信息,所以非常大,你可以设置overlay 参数使后面的数据覆盖吊前面的数据,不过restart的话你也只能从最后一个增量步开始

4.如何重启

你要指定一个重启点,

inp文件里面加上*RESTART, READ, STEP=step, INC=increment就可以了cae中更简单,首先在model->edit attribute里面选择restart,指定前面分析的job名和你想重启动的开始分析步和增量步就可以了,然后在job里面指定重新创建的工作类型,restart,that's all.

5.注意

重启动不能改变你的原始分析中的任何参数,也就是说,你的启动点的模型必须和原始分析中的模型完全一致的,所以不要企图采用restart的方法来改变边界条件,材料参数或者网格的密度等等。这些需要另外的技巧来实现。

Q: 第一次安装abaqus6.2.1,装完exceed6.2&3D后,开始装 6.2.1的Product installation for network licensing,结果每到进度为75%时,弹出界面: An error ocurred during the move data process: -115 Component:Complete Inslattation File Group: fg_common

File:c:\\abaqus\\6.2-1\\cae\\External\\ebt\\adi3xcol.dll 就停止安装了,重新安装依然如此,现在感到头都大了,哪位好心人帮帮我,不深感激?

A: 没有关闭防火墙!

Q: WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENV ALUES.

一般在什么情况下会发生系统矩阵出现负特征值?

A: 很多情况。比如边界条件定义不对;接触定义不对。

Q: 就是手册的umatmst3.for

1)Mises各向同性的子程序中,调用硬化曲线函数中的Table(2,nvalue)这个二维数值保存是应力应变曲线的数据,可是ABAQUS并没有告诉UMAT中Table 数值这些应力应变的数值啊,UMAT里面Table怎么知道的呢?

2)RHS=SMISESS-EG3*DEQPL-SYIELD这个式子是什么意思?尤其是3G*DEQPL是什么意思?

A: 就是Constant=8中的从第三个开始的数值,数组名传递。。

Q: 我在计算时MSG文件出现如下错误,是不是由于节点数太多啦?该如何处理这个问题?

ERROR: SPECIFIED STANDARD_MEMORY V ALUE OF 8000000 IS TOO SMALL TO RUN THE

NALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE V ALUE IS 23477555. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHER

INFORMATION ?

A: 对abaqus_v6.env文件中的STANDARD_MEMORY的值进行修改。

修改Site文件夹下的abaqus_v6.env中的配置,如下:

#

# System-Wide ABAQUS Environment File

# -------------------------------------

pre_memory = 33554432

standard_memory = 33554432 ##*********修改这里******************

#

# NT specific settings 。

Q: WARNING: THE SYSTEM MATRIX HAS 148 NEGATIVE EIGENV ALUES

1 ABAQUS VERSION 6.3-1 DATE 27-NOV-200

2 TIME 22:08:00 PAGE..

For use by None user license from HKS Inc. ..

STEP 1 INCREMENT 1 STEP TIME 0.00

STEP 1 S T E A D Y S T A T E S ..

AUTOMATIC TIME CONTROL WITH -

A SUGGESTED INITIAL TIME INCREMENT OF 0.300?

A: 将你的INITIAL TIME INCREMENT 改小些试试看,不过,一般出现此类问题,多半是你的模型有问题,欠约束或者其它什么的。

Q: ZERO PIVOT是什么意思?

A: zero pivot可以理解为刚度矩阵出了问题,例如奇异。

可能有不同的原因,如:你所模拟的是软化性质材料,该点因破坏等原因而软化至不能吃劲;模型有问题,如约束不够,或者是单元拓扑出错等;或者是你的UMAT中写的[ddsdde]有错;。

Q: abaqus-uamt的老问题,缺少\'ABA_PARAM.INC\'文件?

A: 在cvf6.5调试时,显示缺少\'ABA_PARAM.INC\'文件!

这个没有任何关系的,这个错误将在ABAQUS调用UMAT的时候自动会找到,仅仅有这个错误将没有任何影响的。也就是说,ABAQUS中调用的时候,实际并不存在这个错误。FT,忘了说一句了,你把ABA_PARAM.INC.dp或ABA_PARAM.INC.sp拷到你的程序工作空间后,应该将把ABA_PARAM.INC.dp或ABA_PARAM.INC.sp的后缀.sp或.dp去掉,即将ABA_PARAM.IN C.dp或ABA_PARAM.INC.sp改名为ABA_PARAM.INC。

呵呵,他的意思是在Visual Fortran中调试其子程序,我觉得这是一个好办法,我当时也是这麽办的,毕竟在ABAQUS中调试是非常麻烦的,只有当你的UMAT 没有语法或者明显的逻辑错误,你在ABAQUS中调试才能事半功倍。

Q: abaqus中有没有布尔运算?能否在abaqus中将两个体粘成一个体?

A: 可以,用constraint和interaction可以将两个part进行定位的。

Q: abaqus里做静力分析时怎么加重力加速度?在inp文件里的card是什么?A: 如:*DLOAD

ELSET, GRA V, 9.81, 0., 0., 1. 有疑问查一些*DLOAD卡就清楚了。

Q: 如何在ABAQUS CAE如何设定在某一节点上受集中力呢

请问一下在ABAQUS CAE如何设定在某一节点上受集中力呢??

我用CAE中的Set好象不能设定一个节点丫….

请问在CAE要如何设定,或则祇能在inp下设定呢??

A: 用set可以定义一个节点的。

另外如果只是在一个节点上加集中力就没有必要用set了,在Load模块里creat a load,选择集中力,然后点中该节点(别忘了done)就可以输入力的数值和方向了设定一个节点当然可以了,在cae里一点的集中荷载还是挺容易加了。partition, datumn point, 都可以。用CAE建的模型好象不能直接加在接点上,不过可以加在vertex上。

Q: 如何将两个part接在一起?

几个人将一个工作分开来做,完成后要将造好型的part接在一起,再加边界条件和载荷。该如何把这几个part接起来?

A: 先导入File->import->part,然后形成assembly,这样做试试看如何

Q: abaqu cae中作梁的分析时可否定义不规则的profile形状?

只提供了box,pipe,circle,角钢,方钢,T型钢,l型钢,工字钢,我可不可以定义像三角形空心钢管,或者梯形空心钢管这样的profile?

A: 直接给出I 和G

Q: 在CAE中建立PART时为什么只能先建立在X-Y平面上的模型

在CAE中建立PART时为什么只能先建立在X-Y平面上的模型,可不可以直接建立空间的点,结点线,面,体呀

A: 在图纸上画投影图当然是2维的了。

这是它的建模思路,也是一种专利了。直接在空间里建立模型当然也是一种思路,哪个好呢?各有千秋。你要使用CAE,就要熟悉并掌握它的思路,把整体结果分解为各个小part,把每个part都简化到何以在平面上画出投影图并extrude或回转就得到part的空间实体。然后把这些part给instance,指定其section特性(就是其材料呀,厚度呀等在图形上显示不了的数据)装配起来构成复杂的模型。同一个part的草图可以有许多instance,你可以尽量利用装配功能来减少你建立part 的工作量。对机械图比较熟悉的,可以想象为制作零件然后装配成机器的过程,对于土木类工程师可能开始会不太习惯,但是多做几个例子就会熟悉并喜欢上这种操作了。如果你实在是喜欢空间直接建模法,你可以用ANSYS建立模型划分单元后再用小程序变成ABAQUS的输入文件。

Q: abaqus的所有功能在cae中都能完成吗?还是需要在命令行中输入A: 大部分工作可以在CAE中实现,有些功能还需要在INPUT文件中修改,将CAE转化成*.INP,就可。abaqus的前置确实不太好。cae的功能一直没有长进。如板料成形方面,回弹的分析就是不能在cae里完成,必须写inp文件来做。写inp 文件的效率可一直不高啊因为一般成形我们最好用动态显示算法,而回弹最好静态隐式算法。

不过在cae里就是不行。只要建立了动态算法的step,后面就不能建立静态算法的step。还有啊,abaqus的cae里面是先设置边界条件和输出,然后剖分网格,这个不适合大多数有限元前置处理的思路。很少前置软件是这样的。比起其他的前置软件,如femap,patran,还有工业用软件autoform,pam-stamp,差不少啊!

Q: 如何给一个面的一部分加载

A: Are you using /CAE? If so, please partition your surface and apply loading.

对的,将矩形和圆区分开来,将圆定义为一个块就行了!

将要加载的面上所有的节点设为一个nset就行了

在矩形面上画个circle 然后就在tool-》中partion

Q: 怎么样使用Abaqus最好呢?

1) 用CAE

2)直接填写Keyword文件

3)用Python script

4)用其它前处理及后处理(那么用什么前处理最好呢?Hypermesh还是?怎么样使用呢)

A: 根据具体情况而定,不复杂几何模型模型最好在CAE中建模这是最容易对模型控制的方式,无论是几何还是网格写keyword会累死人的,除非很多相同的部分,简单的装配,当年师兄用C++编了一个程序写斜拉桥input文件,input 文件有N兆,佩服佩服个人认为hype rmesh是最好的前处理软件,做复杂的大模型再好不过python edit国外用的很多,当然用它来继续开发CAE的功能还是很方便!有限元要充分利用各种工具的特点和长处。比如造型,比如剖发网格,比如前置处理。我的习惯思路:

ug 造型,导出模型;femap或者hypermesh 剖发网格;femap作前置,导出inp 文件;手工检查inp文件,部分性的修改;最后用abaqus计算。简单的模型可以在cae里作。最近看着别人用patran,我也临时学点patran作前置,挺好。

Q: 想用patran调用abaqus,如何在ini文件中设置那?

A: 这个问题我也碰到过。安装patran的时候,没有装abaqus的介质。只要重新安装patran选择用户自定义安装。注意abaqus的介质选项即可。

Q: 在Abaqus/CAE中如何施加三角形线荷载?

例如,一根悬臂梁上作用有三角形线荷载。

只找到了均布荷载Uniform,还有就是User Subroutine了,

这么简单的问题应该不用如此复杂吧?

Inp文件中似乎也只有Dload。

在Ansys中就很好处理此问题~~~~~~~~

A: you can use either AMPLITUDE or TABULAR to define load values. If it is a triangle, just give a few data from the triangle load line. Then when you define load, be sure to specify the name of the AMPLITUDE or TABULAR.

Q: Abaqus的背景只能黑色吗?怎么修改?

A: 是个问题呀!我是把图片硬拷贝下来后用其他图形软件改变背景的。还有啊,图片怎么

导出呢?导出是什么格式的呀?

利用Abaqus的print功能可以将图导出到图形文件中,具体路径如下:

File>Print,在弹出的菜单中,不选print viewport background,

rendition可以是white&black,也可以是color.

将Destination由printer改为file,自己定义file name,文件格式由PS改为Tiff图形格式。此文件可以用画笔打开,也可以插入到word中。可惜,图片文件质量太差,还不如抓图得到的质量呢。抓图有个毛病,就是黑色北京需要反色,而反色后,"热"色变成"凉"色,看上去不大好。改变Contour的色谱,用相反排列色谱,反色后,应力最高的反而成了黄色,看着也不爽。大家有好办法的话,我就不用画笔程序再一点点的改了,虽然工作量不大,就是觉得不舒服!我试了一下,觉得还是直接保存为.png格式的图片比较好,而且文件也比较小,可以插入word。的确不错!不过需要注意,保存前最好把那些legend,title block和state block 的字体颜色改为黑色,背景改为其他你喜欢的颜色,这样效果会更好,不至于看不清图上的文字。

Q: ABAQUS建模问题:如何建立复合材料的模型呢?

我看了CAE的文档,但是还是不知道如何建立下列模型。

我要建立一个圆柱体,高度H,半径R=X1,其中x2>R>0为材料1,X1>R>x2为材料2,在CAE中如何实现?不考虑接触,因为两种材料性能相差不大。请高手指点。

A: 我想你应该首先建立两种材料的定义。然后,使用assign命令对于不同的模型部分,分配不同的材料。那要如何操作?首先建立一个圆柱PART,然后建立2种材料,那如何将2种材料分配给一个PART INSTANCE的不同部分?或者如何将2个不同直径的PART组装为一个整体?请指教。

would like to answer your question, although I have never played around Abaqus/cae. Since the geometry of the model is simple (a cylinder), my solution is you do not have to use /cae, and you manually write your input file or rewrite the input file generated by abaqus/cae. I guess I am kind of cheating in answering your question. However, I think some hints could be helpful.

1) Again, take advantage of the symmetry if possible. If this is an axisymmetrical problem, you need to model only a plane. Maybe a quarter model?

2) When you define the elements, use "elset" to group two different elements: such as: elset, elset=material1, generate

1,100 200,300

elset, elset=material2, generate

101,199

which means the elements from #1 to #100, and from #200 to #300 are one material, while elements from #101 to #199 are the other. Then, you could define material properties by material, name=material1 ...... material, name=material2

Q: abaqus的图形如何copy?

A: file>print>file格式为png,可以用Acdsee打开。

Q: 用Abaqus能否计算[Dep]不对称的问题?

A: 可以,并且在step里面的edit step对话框other里面的matrix solver有个选项。Q: 弹塑性矩阵【D】与ddsdde有何联系?

A: stress=D*stran;d(stress)=ddsdde*d(stran)。

Q: 在abaqus中,如果采用umat,利用自己的本构,如何让abaqus明白这种材料的弹塑性应变,也就是说,如何让程序返回弹性应变与塑性应变,好在output 中输出,我曾想用最笨地方法,在uvarm中定义输出,利用getvrm获取材料点的值,但无法获取增量应力,材料常数等,研究了帮助中的例子,umatmst3.inp,umatmst3.for,他采用mises J2 流动理论,我在output history 显示他已进入塑性状态,但他的PE仍然为0!!?

A: 用uvar( )勉强成功。

Q: 本人在用umat作本构模型时,*static, 1,500,0.000001,0.1 此时要求的增量步很多,即每次增量要很小,*static 1,500 时,在弹性向塑性过度时,出现错误,增量过大,出现尖点?

A:这样尝试:

*STEP,EXTRAPOLATION=NO,INC=2000000

*STATIC 0.001,500.0,0.00001,0.1。

Q: 模型中存在两个物体的接触,计算过程中报错,怎么回事?

A: 接触问题不收敛有两个方面不妨试试:

一、在*CONTACT PAIR 里调试ADJUST参数;

二、调一些模型参数,比如FRICTION等。

Q: 在边界条件和加载时,总是有initial这个步,然后是我们自己定义的加载步,请问这个initial步,主要作用是什么?能不能去掉?

A: 不能去掉,所有的分析都有,是默认的步。

Q: A solid extrusion base feature 这句话是什么意思?extrusion、revolution等是什么意思?

A: 这两的是三维建模时候,在画完二维图形,如何来生成三维图形, extrusion意思是你给定一个厚度,然后二维图形第三个方向上面伸展这么多形成三维图形revolution意思是你给定一个旋转轴,二维图形绕其旋转后形成三维轴对称图形。Q: 偶在umat中调用求主应力函数CALL SPRINC(STRESS,PS,LSTR,NDI,NSHR) 后,存储主应力得数组PS中各个主应力排列顺序是什么?

PS1>PS2>PS3 ?

PS1PS1>PS3>PS2 ?

A: 第二个。

Q: 在*USER MATERIAL的定义中,Δσ对Δε的偏倒数,即

DDSDDE被称为:

"material stiffness matrix";

而在UMAT中,DDSDDE被称为:

"material\'s Jacobian matrix"。

请问DDSDDE和材料的切线刚度矩阵的关系是什么?二者是一个概念么?

A: 一般说可以这样理解:

σ=ddsdde*Δε

有点像我们常说的弹塑性矩阵:Dep。

Q: 请问field output和history output什么区别?

关键字*node output和*node print有什么区别??

A: field output和history output在viewer模块中很明显。field output是场量输出,history output为历史数据输出,会记录场量随时间的曲线。至于*node output与*node print的区别在于他们写入不同的后缀文件,比如*.odb,*.dat等,在help中有详细介绍*node output是给CAR或VIEWER做后处理用的,是二进制文件,而*node print则是写到DAT文件中,你可以自己看的文本文件。

field output用来输出模型中较大部分(a large portion of the model)的那些输出频率较小的变量,如模型的等值线,变形图等;而history oupput用来输出模型中

较小部分的那些输出频率较大的变量,如荷载作用点处的荷载--位移图等。所以,应力,应变,位移,反力等既可以作为field output也可以作为history oupput 输出。

Q: 我在学习ABAQUS时,遇到以下专有名词想请教各位。

1、orphan mesh instance,

2、self-contact,

3、elastic foundation,

4、convective interaction,

5、amplitude,

6、solution-dependent state variable,

7、datum,

8、thermal film condition ?

A:

1、独立网格实体;

2、自接触;

3、弹性基础,如弹性地基梁;

4、对流相互作用;

5、幅度,数值大小,如定义随时间变化的荷载;

6、依赖于解的状态变量,在UMAT中经常用到,它的大小取决于某一增量步收敛解;

7、辅助数据,在CAE中经常用于定义数据点、刚体参考点、辅助平面等;

8、

Q: (1),请问如何得到M,C,K矩阵?用什么命令?

(2)ABAQUS中能实现这样的东西吗?就是我需要平滑ABAQUS产生的位移场,还是这个平滑只能在其他环境中完成?

A: (1)试一试:*element matrix output

(2)试试这个输出选项:

*El print,Position=average at Nodes。

Q: solution-dependent variable 和time-dependent variable 这两种变量有什么差别?

A: solution-dependent variable从字面上理解,即为"与解答有关的变量",它的值与每一增量结束时的变量有关,常用于UMAT中,即SDV,可以由*DEPDV来定义其个数。

time-dependent variable从字面上理解,即为"与时间有关的变量",常用于定义随时间变化的量如地震荷载等,可以由*APPLITUDE来定义。。

Q: S4R单元可以输出sth(就是section thickness单元厚度变化),但C3D8R却不能输出sth。该怎么看单元层厚度方向的变化?

A: 个人觉得定义了局部坐标方向,特别是对于各向异性问题,有限元在材料方向上计算应力,应变,输出也就在材料方向,假如第三方向为厚度方向,那第三方向的塑性应变就是他的厚度变化。

Q: ABAQUS 在运行过程中,对于一些比较大的问题,经常会出现以下的信息:ABAQUS Info: License Timeout set to 70(可能是不同的值) minutes. ?

A: 个人感觉是没有问题的,你放心吧。

Q: 在计算固结过程中,msg文件总是有一下信息***ERROR: TOO MANY ATTEMPTS MADE FOR TTHIS INCREMENT: ANALYSIS TERMINATED

改变其最小时间变量也不行,太郁闷了,是不能收敛吗?请指教可能是什么原因? A: 很多原因,比如边界条件不对;约束不够;接触定义不对;单元划分畸变等等。

Q:

A: 个人觉得:

umat实现自己的本构没有固定的方法,对于不同的本构有可能必须采用不同的方法。这要靠自己不断地摸索。有可能一种方法对于简单加载问题还行,但有可能对于复杂问题并不收敛。最重要一点,就是umat中采用的算法必须consistent.再就是ddsdde必须正确,(如果采用back_Euler 方法等一些算法,ddsdde错误有时不影响结果(对于简单加载问题没有影响,能收敛,),但对于复杂问题不收敛。

uptonow,你这个算法对于Mises,hill,J2,J2d等一类的屈服函数是正确的,但具体的本构还要灵活运用,这我也正学习,正在摸索。

有时,umat需要很强的有限元基础,并且对采用的本构要很熟悉,不要在一颗树上吊死才好。首先要确认自己的umat没有错误,如果没有,但就是不收敛(在不断减小加载步长的情况下,当然最好对步长不敏感,特别是对于粘弹性,粘塑性,内变量一类的材料,有的本构取决于背应力的计算)。那就应该考虑换一种算法。

下面续贴本人摘自剑桥大学BBS上的FAQ:

1. General Questions

Q1.1 : How much disk space do I need?

Minimum of 10 MBytes to run simple examples if using PATRAN for Pre or Post processing. Probably about 20 MBytes for medium sized problems.

The amount of disk space required depends on the numbers of nodes/elements present in the mesh and the output frequency requested for outputting data to the *.fil, *.dat and *.res files.

The following ABAQUS files are created in the user's home directory.

*.dat, *.fil, *.res, *.msg, *.sta, *.log, *.job, *.023 (deleted at the completion of the ABAQUS job).

________________________________________

Q1.2 : What type of analysis can I do?

The procedures available in ABAQUS are listed below (in alphabetical order) : ?BUCKLE

?COUPLED TEMPERATURE-DISPLACEMENT (steady state and transient)

?COUPLED THERMAL-ELECTRICAL (steady state and transient) ?DYNAMIC

?FREQUENCY

?GEOSTATIC

?HEAT TRANSFER (steady state and transient)

?MASS DIFFUSION (steady state and transient)

?MODAL DYNAMIC

?RANDOM RESPONSE

?RESPONSE SPECTRUM

?SOILS, CONSOLIDATION

?SOILS, STEADY STATE

?STATIC

?STEADY STATE DYNAMICS

?VISCO

________________________________________

Q1.3 : What output options (hard copy) are there?

With ABAQUS 5.7 any size plots (including A3 and A4) can be obtained.

If using PATRAN of FEMGV hard copies of any size (includng A4 and A3) can be obtained.

________________________________________

Q1.4 : What pre and post-processing programs are available other than ABAQUS/Post? PATRAN and FEMGV are available in the teaching system.

Both these programs read the results file *.fil for post processing.

If you are post processing using FEMGV and have written results at element Gauss points then you need to write the gauss point co-ordinates as well. Include the following lines in the ABAQUS input file :

*EL FILE, POSITION=INTEGRATION POINTS

S

E

COORD

Otherwise it will not be possible to use the FEMGV program to post process the element gauss point results ie FEMGV would not know the location of the gauss points.

________________________________________

Q1.5 : Is ABAQUS/Pre available?

No. ABAQUS/Pre is not available in the teaching system. However PATRAN which provides the same functionality and which has the same look and feel of ABAQUS/Pre is available in the teaching system. Type patran to run the PATRAN program.

Q1.6 : Is ABAQUS/Explicit available?

ABAQUS/Explicit is only licenced on tw900. You need to rlogin to tw900 to run this version of ABAQUS. ABAQUS/Standard is also available on tw900.

________________________________________

Q1.7 : What on-line documentation are available and how do I access it?

Type abaqus57 doc in the CUED teaching system to access the full set of the ABAQUS Users' manual (V olumes I, II and III).

At present this is the only on-line documentation available except for ABAQUS Release Notes (100 pages) and ABAQUS Site Guide (80 pages).

________________________________________

Q1.8 : What units are used in ABAQUS?

There is no inherent set of units used in ABAQUS. It is up to the user to decide on a consistent set of units and use that units. Typical sets of units :

1 2

Length - metres mm

Force - Newtons Newtons

Time - second second

Mass - Kg tonne (**)

Density - Kg/m3 tonne/mm3

Stress - N/m2 N/mm2 (= MPa)

Young's Modulus - N/m2 N/mm2 (= MPa)

** 1 tonne = 1000 kilograms

Decide on the units before you start preparing your data. This is critical. If you start typing in the nodal co-ordinates that means you have already decided on what units to user for the Length parameter. This only leaves the choice for the units of Force.

It is not a good idea to choose mm for length, Newton for force and then specify the Young's modulus in KN/m2.

2. Jobs

Q2.1 : How do I run small jobs?

Use the following command

abaqus job=job-id interactive

Example : abaqus job=plate

________________________________________

Q2.2 : How do I go about running many small jobs?

Create a batch file (say aba.run) with one line per analysis as shown below :

abaqus job=analysis-a interactive

abaqus job=analysis-b interactive

abaqus job=analysis-c interactive

abaqus job=analysis-d interactive

Then make the file an executable using the following unix command :

chmod u+x aba.run

Then type aba.run to execute the ABAQUS jobs one at a time while you are logged ON. This is only suitable for small jobs which only take a few minutes to run. These jobs will run one at a time and in sequence.

This is preferable to submitting all the jobs at the same time (for example typing the above commands directly at the terminal without the interactive parameter). This will

put a strain on the server and its resources and inconvenience the other users as well. For medium to large jobs use the batch command available in the CUED teaching system.

________________________________________

Q2.3 : How do I run large jobs using batch?

In the CUED teaching system use the batch command. To run the job on tw500 or tw900 servers use :

batch -QX -mbao "abaqus job=job-id interactive"

To run a job specifically on the tw900 server use :

batch -QN -mbao "abaqus job=job-id interactive"

Example :

batch -QN -mbao "abaqus job=plate interactive"

The progress of the submitted job can be monitored using the batchq command. Use the batchrm command to delete any batch jobs you had submitted before these are run, if you change your mind. See the man pages on batchq, batchrm for more details. Example : Type man batchq.

________________________________________

Q2.4 : How do I run ABAQUS/Post?

Using the following command :

abaqus post job=job-id

Example : abaqus post job=cantilever

________________________________________

Q2.5 : How do I run ABAQUS/Plot?

Using the following command :

abaqus plot job=job-id device=cps or hgl or x11

Use device=x11 to view the plots on the screen. Use one of the other options (cps for colour postscript or hgl for hpgl) to create a hard copy of the plot.

Example : abaqus post device=cps job=cantilever

Q2.6 : How do I get copies of the ABAQUS examples input data files (*.inp)?

The ABAQUS datafiles used in the examples manual can be found in the /export/abaqus/samples/exastd directory. Similarly the datafiles in the verification manual can be found in /export/abaqus/samples/verstd directory. These files will have the extension name inp. All file names are in lowercase. use the cp command to copy the relevant file.

Example : cp /export/abaqus/samples/exastd/1010101.inp .

This will copy the 1010101.inp file to the current directory.

If these directories don't exist then can use the abaqus fetch command.

Example : abaqus fetch job=1010101

________________________________________

Q2.7 : How do I run a ABAQUS job which uses a user subroutine?

Using the following command :

abaqus job=job-id

As can be seen this is no different from running a standard abaqus job. The user subroutine itself can be embedded in the abaqus input file. Here it is illustrated with the umat subroutine.

<....part of the abaqus input file ....>

........

........

*END STEP

*USER SUBROUTINES

SUBROUTINE UMAT(........)

........

........

END

........

Alternatively the user subroutine can be in a separate file (say my_material.f) and the INPUT parameter is set to that file name.

<....part of the abaqus input file ....>

........

........

*USER SUBROUTINES, INPUT=my_material.f

........

........

________________________________________

Q2.8 : How do I run a user written post processing program which accesses the *.fil file?

Using the following command :

abaqus make job=job-id user=name-of-file

Example : abaqus make job=cantilever user=disp1

This will compile the user program in a file called disp1.f and then create an executable called cantilever.x. Type cantilever.x to run this program.

________________________________________

Q2.9 : How do I find out about the different execution procedures that are available with ABAQUS?

Type abaqus help and this will list all the abaqus execution procedures. These are listed below :

Execution Procedure for ABAQUS/Standard and ABAQUS/Explicit

abaqus job=job-name [ analysis | datacheck | continue | help | recover |

convert={restart|select|all} |

information={environment|local|memory|release|status} ]

[ input=input-file ] [ user=source-file ]

[ oldjob=oldjob-name ] [ fil={append|new} ]

[ globalmodel=results file-name ] [ double ]

[ memory=memory-size ] [ buffer=buffer-size ]

[ interactive | background | queue=queue-name ]

[ cpus=number-of-cpus ] [scratch=scratch-dir]

[ subcomplex=subcomplex-name]

Note: subcomplex is only valid on the Convex Exemplar

Execution Procedure for ABAQUS/Post

abaqus post [ job=job-name ]

[ restart=restart-name ] [ input=input-file ]

[ device=device-name ] [ display=display-name ]

[ geometry=widthXheight+xpos+ypos ]

[ memory=memory-size ] [ buffer=buffer-size ]

Execution Procedure for ABAQUS/Abares

abaqus abares job=job-name [ restart=restart-name ]

[ beginstep=step-number ] [ endstep=step-number ]

[ increment={all|endstep|final|none|integer-list} ] Execution Procedure for ABAQUS/Plot

abaqus plot job=job-name [ input=input-file ] [ device=device-name ]

[ frame={all|integer-list} ] [ options=options-file ]

Execution Procedure for ASCII translation of results (.fil) files

abaqus ascfil job=job-name [ input=input-file ]

Execution Procedure for on-line documentation

abaqus doc

Execution Procedure for ABAQUS/Append

abaqus append job=job-name oldjob=oldjob-name input=input-file

Execution Procedure for ABAQUS/Fetch

abaqus fetch job=job-name [ input=input-file ]

Execution Procedure for ABAQUS/Findkeyword

abaqus findkeyword [ job=job-name ] [ maximum=maximum-matches ]

Execution Procedure for ABAQUS/Make

abaqus make job=job-name [ user={source-file|object-file} ]

Q2.10 : How do I find what the current settings are for the environment variables? Type abaqus info=environment and this will list all the current setting of the ABAQUS environmental variables.

memory, local, release, status are other options on which you can get more information on.

________________________________________

Q2.11 : How do I change the current settings of the environment variables?

Create a file called abaqus.env in the directory from which ABAQUS is run which contains lines of environment variables you want to change set equal to new values. However make sure that the computer on which you are running ABAQUS can support the changes. For example you can increase the memory used by ABAQUS. But you cannot increase this beyond what is available in the computer.

Example : Include the following lines in the abaqus.env file to increase the size of post_buffer and the post_memory used by ABAQUS/Post.

post_buffer="1000000"

post_memory="3000000"

________________________________________

Q2.12 : Can I run long jobs on the twgs?

No. is is not a good idea to use the twgs to run long ABAQUS jobs. It is better to use the batch system on one of the faster servers (example : tw900) for long jobs.

Also becasue of the limited resources (memory, swap space, scratch disk space) large jobs also cannot be run on the twgs.

3. Elements

Q3.1 : How do I find the positive normal of a shell element?

For shells the positive normal is given by the right-hand rule going around the nodes of the element in the order they are given in the element-nodal connectivity data line (datalines which follow the keyword *ELEMENT line).

Q3.2 : What is the difference between a shell element and a 2D solid element?

The 2D Planar elements (For example : CPE4, CPE8 - plane strain analysis and CP4, CPS8 - plane stress analysis) are only used in situations where the loading is confined to the plane of the elements. The elements only have planar variables (d.o.f) ux, uy. Shell elements are needed for out-of-plane loading. Consider a square plate subjected to a loading normal to the plane of the plate. This requires shell elements and use of plane stress/plane strain type of elements would be inappropriate under these circumstances.

Q3.3 : How do I specify the local-1 direction for a beam in space?

See the section 6.1.2 of the Getting Started with ABAQUS/Standard manual for a description of this.

There are a number of ways of specifying this. In the following figure n1 represents the local-1 direction. Consider a beam element of rectangular cross-section. Then local-1 direction is the direction parallel to the width (or base of the element). See Figure 15.3.4-1 and the figure in Page 15.3.9-15 of the ABAQUS User's manual(Ver 5.7).

Q3.4 : Why do you need to specify the local-1 direction for a beam in space? Consider a non-circular cross-section (example : rectangular section). The bending stiffness is affected by which way the width is oriented. local-1 direction fixes the orientation without any ambiguity.

________________________________________

Q3.5 : I am having a problem interpreting the stress output from a shell element?

If you are using the shell elements (example : S4, S8R5) then you need to be aware that the stresses are defined in local directions which in turn are dependent on the orientation of each element w.r.t the global axes.

The Shell element sign convention explains the sign convention for the local directions.

4. ABAQUS - Mesh

Q4.1 : I have created a mesh in two parts separately and have ended up with two sets of node numbers along a common edge. Is there any alternative to editing the *.inp file to replace one set of node numbers by the other set in the element-nodal connectivity list ie element entries?

Yes. Use the TIE option available under MPC (Multi Point Constraint) to tie the respective nodes.

5. ABAQUS - Materials

Q5.1 : How do I find what material properties are needed for a particular analysis ? Read the relevant section in Chapter 6 : Analysis Types (User's manual V ol. I). This gives an overview about the analysis and has more information about the material properties.

Read also the following sections in Chapter 9 : Materials Introduction of the ABAQUS User's manual.

?Section 9.1.1 - Material Library : Overview

?Section 9.1.2 - Material Data Definition

?Section 9.1.3 - Combining Material Properties

Section 9.1.3 lists the material model combination tables. Several models are available to define the mechanical behaviour (elastic, plastic).

Some material options require the presence of other material options. Some exclude the use of the other material options. For example *DEFORMATION PLASTICITY completely defines the material's mechanical behaviour and should not be used with *ELASTIC.

Once you have all the relevant keywords to define the material properties consult the

keyword section - Chapter 23 (User's manual V olume III) for each of the keywords. This will explain what data is required for each of the keyword.

________________________________________

Q5.2 : What material properties need to be specified in a thermal-electrical analysis ? Referring to Section 9.1.3 of the ABAQUS User's manual you will require the heat transfer properties as well as the electrical properties. These are listed below : ?Heat Transfer properties

o *CONDUCTIVITY

o *LATENT HEAT

o *SPECIFIC HEAT

o *HEAT GENERATION

?Electrical properties

o *DIELECTRIC

o *ELECTRICAL CONDUCTIVITY

o *JOULE HEAT FRACTION

o *PIEZOELECTRIC

This forms the complete set of properties. If Piezoelectric elements are not used then *PIEZOELECTRIC and *DIELECTRIC properties will not be required.

If only the steady state heat transfer response is of interest then *SPECIFIC HEAT properties are not required. Similarly if there are no phase changes involved then *LATENT HEAT is not required.

*JOULE HEAT FRACTION is used to specify the fraction of electrical energy that will be released as heat.

Example problem 3.2.24 - thermal-electrical modeling of an automotive fuse illustrates the thermal-electrical analysis.

ABAQUS allows for redundant material properties to be specified. It will simply ignore the material properties not required for the current analysis.

Typical example of material properties :

*MATERIAL, NAME=ZINC

*CONDUCTIVITY

0.1121, 20.0

0.1103, 100.0

*ELECTRICAL CONDUCTIVITY

16.75E3, 20.0

12.92E3, 100.0

*JOULE HEAT FRACTION

相关主题
相关文档
最新文档