[分享]Abaqus 使用问答

发信人: dengnch (dengnch), 信区: CAD
标 题: Abaqus 使用问答
发信站: 哈工大紫丁香 (2003年06月01日10:47:48 星期天), 站内信件

Abaqus 使用问答

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: YOU CAN TRY AS FOLLOWS:
*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
从字面上理解，即为"与解答有关的变量"，它的值与每一增量结束时的变量有关，常用于U
MAT中，即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没有错误，如果没有，但就是不收敛（在不断减小加载步长的
情况下，当然最好对步长不敏感，特别是对于粘弹性，粘塑性，内变量一类的材料，有的
本构取决于背应力的计算）。 那就应该考虑换一种算法。

一点体会，请大家探讨。

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 EIGENVALUES.
一般在什么情况下会发生系统矩阵出现负特征值？
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 VALUE OF 8000000 IS TOO SMALL TO RUN THE
NALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUE 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 EIGENVALUES
1 ABAQUS VERSION 6.3-1 DATE 27-NOV-2002 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中写的[dds
dde]有错；。

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, GRAV, 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建立模型划分单元后再用小程序变成AB
AQUS的输入文件

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的介质。 只要重新安装patr
an选择用户自定义安装。注意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格式的图片比较好，而且文件也比较小，可以插入w
ord。的确不错！不过需要注意，保存前最好把那些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种材料分配给一个P
ART 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
......