[转帖]Abaqus 使用问答

TNC

　　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 ?

　　PS1

　　PS1>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 输出。 。

TNC

　　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 Inf 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中写的[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中调试才能事半功倍。

yejet

　　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建立模型划分单元后再用小程序变成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兆，佩服佩服 个人认为hypermesh是最好的前处理软件，做复杂的大模型再好不过 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格式的图片比较好，而且文件也比较小，可以插入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. Sincethe 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
　　[ 本帖最后由 vibshare 于 2006-9-2 14:56 编辑 ]