首页 > 心得体会 > 工作体会 / 正文

abaqus常用技巧总结

2020-07-02 18:40:36 ℃

　Abaqus常用技巧总结

　建议阅读方式：

　ABAQUS常用技巧总结 1

　目录 2

　1.对TIME INCREMENT的根本理解 4

　2.ABAQUS 请问 MOMENT的加载 4

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

　4.CAE中如何加预应力 5

　5.HYPERMESH里面看到ABAQUS分析的结果 5

　6.X-Y PLOTS 5

　7.把上一次的分析结果作为下一次分析的初始条件该怎么做 6

　8.材料方向与增量步 7

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

　10.关于ABAQUS的任务管理 8

　11.关于数据的输入输出 9

　12.后处理积分 10

　13.接触分析激活杀死 10

　14.利用QUEUE的功能由本地机器向远程UNIX机器提交ABAQUS作业的方法 [精华] 10

　15.利用命令进行计算时如何设置调用内存量 13

　16.清华大学BBS的ABAQUS精华 13

　17.请问怎么实现双曲线 36

　18.取消坐标系等的显示 37

　19.如何在计算中修改材料特性 37

　20.输出计算过程中的总质量和总刚度矩阵 40

　21.先张预应力：

　22.用户子程序的使用 41

　23.怎样设定用双CPU机器进行ABAQUS计算 41

　24.中途停止正在运算的JOB 42

　25.自适应网格技术 42

　26.ABAQUS计算与内存 42

　27.质量缩放 43

　28.ABAQUS多处理器进行并行计算的效果研究 54

　29.YAHOO讨论组摘录--CONTACT+OVERCLOSURE 55

　30.原创：无限元建立方法，希望得到加分 70

　31.[分享]ABAQUS 使用问答 76

　32.[转帖]ABAQUS6.4导入外来模型的几点小经验！ 87

　33.ABAQUS的多图层绘图 88

　34.子结构 89

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

　36.模型的重启动分析－RESTART 91

　37. ABAQUS的单位心得 92

　1.对time increment的根本理解

　abaqus的step里有maximum number of increment、initial increment、minimum increment 、maximum increment四个量许多网友不知怎样设置合理，合理设置是建立在深刻理解基础上的。

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

　但是如果时间步长减小到最小值minimum时计算结果还是不收敛，则abaqus将停止计算，由此可知maximum值和minimum值分别是abaqus在收敛计算时时间步长的上下限，同时total time=求和(time increment*number)，当时间步长很小时，需要计算的步数number相应增大(电脑计算花的时间也随之增大)，因此number一般要设置较大值。minimum并不是越小越好，因为1)number即计算时间增大2)abaqus计算精度约在10^(-5),当时间步长小于这个值，计算结果已经没什么意义了。

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

　MF.c*b.i6＼#s#fwww.simwe.com 总而言之，maximun number要适当设置较大值，initial可适当改小（如-2，-3量级），minimum（-5量级）不要修改，maximum值影响不大，可不改

　2.moment的加载

　请教各位大牛

　我最近在算一个题目,一个大筒体上有三个接管端面固定,大筒体两端加载扭距,但是不知道如何加载.请各位高手能赐教

　谢谢

　1. 将大筒体两端要施加扭矩的节点分别定义为两个Nset:left, right.

　 2. 分别在大筒体两端的圆心处定义两个reference node: rp-left,rp-right.

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

　*KINEMATIC COUPLING, REF NODE=rp-left

　left, 6,6

　*KINEMATIC COUPLING, REF NODE=rp-right

　right, 6,6

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

　友情提醒:reference node也有自由度,注意相应的边界条件.

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

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

　4.CAE中如何加预应力

　具体没作过，看看*PRESTRESS HOLD和*INITIAL CONDITIONS, TYPE=SOLUTION, REBAR

　这两个命令以及ABAQUS Analysis User＇s Manual

　 “Defining reinforcement,” Section 2.2.3

　 “Defining rebar as an element property,” Section 2.2.4

　这两节吧。

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

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

　2、利用hyperworks提供的hmabaqus.exe(在安装目录下的Altair＼hw7.0＼translators中)

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

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

　gook luck!

　6.X-Y Plots

　** STEP: pre-load

　*Step, name=pre-load, nlgeom

　pre-loading

　*Static

　0.01, 1., 1e-05, 0.1

　........

　** LOADS

　** Name: pt-load Type: Concentrated force

　*Cload

　_G5, 2, -200.E6

　........

　*Output, history, frequency=1

　*node output, nset=_G5

　CF2,U2

　*element output, elset=_G5

　E22, S22

　*monitor, node=_G5, dof=2

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

　Q如果想把上一次的分析结果作为下一次分析的初始条件该怎么做

　A使用LDREAD命令

　首先需要注意下面两个问题

　1每一个ANSYS的实体模型的面或体都要定义对应的单元类型编号材料属性

　编号实常数编号单元坐标系编号这些参数在整个分析过程中保持不变而这些

　编号对应的属性在各个步骤中是不同的

　2网格划分要满足所有步骤的要求单

　元类型必须兼容

　步骤

　1创建实体模型

　创建多个物理环境

　设定一个物理环境中的单元类型材料属性实常数坐标系等

　将这些参数的编号赋给实体模型的面或体

　施加基本物理载荷和边界条件

　设定求解选项

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

　3清楚当前的物理环境命令是PHYSICS, CLEAR

　4重复第二步准备下一个物理环境

　8.材料方向与增量步

　材料方向：

　针对各向异性材料（如板金材料、复合材料等）变形体，材料方向定义材料的某一特定方向如纤维方向。该方向随着变形体的移动而移动，旋转而旋转。输出的场变量值都以材料方向为参考，有时便于数据处理，分析计算结果等。

　 abaqus中增量步

　abaqus 中把所有载荷按一定的要求分成若干载荷步step，每一步step根据abaqus自动载荷增量，分成若干增量increments，每一增量施加一定的载荷，然后每一增量通过若干迭代步iteration 进行迭代，当系统达到平衡时，迭代结束，完成一个增量。当所有的增量都完成后，计算结束；反之，计算可能出现发散。这时，可以通过采用多钟方法（如调整放大质量系数，单元网格优化等）调整增量大小，使计算继续进行

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

　我喜欢用fortran

　写成这样：

　 ! nodedeform.f90

　! FUNCTIONS:

　nodedeform - Entry point of console application.

　Example of displaying ＇Hello World＇ at execution time.

　!****************************************************************************

　! PROGRAM: nodedeform

　! PURPOSE: Entry point for ＇Hello World＇ sample console application.

　!****************************************************************************

　program nodedeform

　USE DFLIB

　implicit none

　result = SYSTEMQQ (＇abaqus job=3d05ctJ interactive＇)

　result = SYSTEMQQ (＇abaqus job=3d05ctJmodel-2 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-1 interactive＇)

　! result = SYSTEMQQ (＇abaqus job=test-2 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-3 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-4 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-5 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-6 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-7 interactive＇)

　result = SYSTEMQQ (＇abaqus job=test-8 interactive＇)

　end program nodedeform

　10.关于ABAQUS的任务管理

　关于ABAQUS的任务管理 [精华]

　 chibiko

　发帖: 26

　积分: 1

　于2004-07-03 03:56

　 --------------------------------------------------------------------------------

　前些日子的计算量很大，自己的机器又不是很快，一个任务要算上几天。一开始，在算题的时候由于CPU占用很大，基本上做不了别的。捉摸了好几天，现在终于有点掌握ABAQUS的任务管理方法了，主要有两点，这里和大家共享。我的系统是WindowsXP，使用的是ABAQUS6.4。

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

　任务暂停：

　 abaqus suspend job=job-name

　任务恢复：

　 abaqus resume job=job-name

　任务终止：

　 abaqus terminate job=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任务的级别调上去，这样就可以娱乐、工作两不误了。

　You may try this method which was posted at simwe forum, but not my original idea.

　 Anyway, I will feel very happy even if it probably just helps you just a little bit.

　11.关于数据的输入输出

　1.输出数据到dat文件：

　 *NODE PRINT, NSET=nset_name，FREQ=1

　COORD

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

　偶在前面已发文问过，偶有时得到的是变形前的，有时得到是变形后的，一直没有弄明白这个问题。偶想得到的是变形后的坐标

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

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

　输入abaqus job=job-1 resume.

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

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

　max_history_requests=0 即可。

　12.后处理积分

　CAE自带此功能，比如对速度时程积分得到位移时程，可以这样操作：

　利用XY Data---->ODB History Output将速度时程存为V，然后XY Data---->Operate on XY Data---->Operations

　看到那个 integrate (X)

　这个可能还是只是实现对时间的积分，如果对dx，dy的积分

　如何实现呢？

　比方说，Q=v*A

　已经知道流速的矢量分布，怎么得到流量的分布了！

　这也是可以的，关键在于你的x-y data数据的两轴数据是什么？其积分就是y数据值对x值的积分！

　看看：CAE user manual：33.4.4Overview of X–Y data operations

　13.接触分析激活杀死

　Usage: Use both of the following options:

　*MODEL CHANGE, TYPE=CONTACT PAIR, ADD

　surface_1, surface_2

　*CONTACT INTERFERENCE

　surface_1, surface_2,

　 Usage: *MODEL CHANGE, TYPE=CONTACT PAIR, REMOVE

　14.利用queue的功能由本地机器向远程unix机器提交abaqus作业的方法 [精华]

　利用queue的功能由本地机器向远程unix机器提交abaqus作业的方法

　假设：

　 1. 远程unix机器的安装目录为/abaqus，

　2. 远程计算机IP地址为192.168.0.166，主机名为ibmlinux

　3. 本地计算机IP地址为192.168.0.18，主机名为training

　下面的讨论涉及的内容相应改变

　Step 1 分别设置本地计算机和远程计算机之间的主机名和IP地址对应。

　 1. 修改本地hosts文件，windows上一般位于c:＼windows＼system32＼drivers＼etc下，unix上位于/etc下，加入远程计算机主机名和IP地址对应行，如：

　 192.168.0.166 ibmlinux

　2. 修改远程hosts文件，位于/etc下，加入本地计算机主机名和IP地址对应行，如：

　 192.168.0.18 training

　Step 2 设置本地计算机对远程计算机运行rsh和rcp命令正确。

　 1. 两台计算机登陆名和密码一致。举例而言，若远程计算机登陆名和密码为root/root，则设置本地计算机也用同样的用户名和密码登陆。

　 2. 修改远程计算机.rhosts文件，该文件位于对应用户名登陆后的主目录下，如/home/root，加入本地计算机的IP地址使得本地计算机可以访问到远程计算机。

　 3. 在本地计算机的提示符下键入

　rsh ibmlinux ls

　测试rsh是否能够正常使用，如果可以列出远程计算机主目录下文件，代表rsh和rcp工作正常。

　 Step 3 设置远程计算机该用户的默认登陆程序为csh，修改passwd文件，位于/etc下，如：

　 root:Ty91eFGzybEE2:0:3::/:/usr/bin/csh

　Step 4 修改远程计算机的abaqus环境文件abaqus_v6.env

　在环境文件最后加入

　queue_name =[“nice_queue“]

　after_prefix = “-b“

　queue_prefix = “-q“

　nice_queue = “/bin/sh -c ＇nice ./%S 1>%L 2>&1 &＇“

　保存

　在提示符下面键入

　abaqus fetch job=beam.inp

　abaqus job=beam queue=nice_queue

　如果可以顺利获得beam.odb，则表示远程计算机queue设置正确。

　 Step 5 修改本地计算机的abaqus环境文件abaqus_v6.env

　在环境文件最后加入

　def onCaeStartup():

　def makeQueues(*args):

　import os, driverUtils

　driverName = ＇./abaqus＇

　scratchDir = ＇/abaqus/Commands＇

　session.Queue(name=＇nice＇,

　queueName=＇nice_queue＇,

　hostName=＇ibmlinux＇,

　driver=driverName,

　localPlatform = NT,

　fileCopy=ON,

　directory=scratchDir)

　addImportCallback(＇job＇, makeQueues)

　保存

　Step 6 在本地计算机上启动abaqus/cae，创建或打开一个模型，在Job模块创建一个新任务，选择Run Mode为Queue，选择nice，OK，提交运算，观察是否任务被自动提交到了远程计算机的/abaqus/Commands下运行并自动返回odb文件到本地计算机。如有错误，会在启动cae的提示符窗口中出现，检查错误并找到相应解决办法。

　15.利用命令进行计算时如何设置调用内存量

　abaqus job=job-name memory=memory-size

　16.清华大学bbs的abaqus精华

　【在 luxz (panda) 的大作中提到: 】

　: 哪位有什么好办法？

　: 感觉abaqus自己的前处理用得很不爽

　: patran可以么？

　我感觉abaqus/cae不错的呀,实体建模颇有pro/e的风格

　感觉patran也好不到哪去,太麻烦,

　也许patran功能更强一些(要不然也对不起它那硕大的身躯呀),

　但对于不是很庞大的模型,根本体现不出来什么优势

　 : air1大侠，本人给你作广告，为何不可？

　: 请问：弹塑性矩阵【D]与ddsdde有何联系，

　: 你用过板壳单元吗？

　 stress=D*stran?

　d(stress)=ddsdde*d(stran)

　那应该就是一样的，因为全量理论，

　Sij=DijklEkl（满足张量求和约定）

　即Stress=D*Strain;

　而在增量理论中，

　△S=D*△E(在有限变形中，△其实应该为应力的客观率）

　似乎不对吧

　大变形下此D非彼D

　你看过黄克智的固体本构关系这本书么

　如果你从全量理论和增量理论的角度上讲

　那似乎第一个Digkl就不对

　你有第一个式子么

　如果有，求导不久完了？

　: 那应该就是一样的，因为全量理论，

　: Sij=DijklEkl（满足张量求和约定）

　: 即Stress=D*Strain;

　: 而在增量理论中，

　: △S=D*△E(在有限变形中，△其实应该为应力的客观率）

　是啊，大变形下的[D]与普通意义下的[D]在构型上是不一样的，

　毕竟[D]大是变形历史的函数，而[D]小则不是，

　我推导一种新的本构关系，

　△Sij=Dijkl△Ekl (其中△为Jaumann率）

　假设材料一开始就屈服（即屈服面为0）

　想用壳单元，

　: 似乎不对吧

　: 大变形下此D非彼D

　: 你看过黄克智的固体本构关系这本书么

　: 如果你从全量理论和增量理论的角度上讲

　: 那似乎第一个Digkl就不对

　: 你有第一个式子么

　: 如果有，求导不久完了？

　: 是啊，大变形下的[D]与普通意义下的[D]在构型上是不一样的，

　: 毕竟[D]大是变形历史的函数，而[D]小则不是，

　: 我推导一种新的本构关系，

　: △Sij=Dijkl△Ekl (其中△为Jaumann率）

　: 假设材料一开始就屈服（即屈服面为0）

　: 想用壳单元，

　唉，别提了，问题就出在，在abaqus中，明明写着可以考虑剪切效应，

　可我打印出剪切力个数是，nshr=1,即只有S12,

　那我的S13,S23就不知怎么计算，

　（DDSDDE(5,5)无法计算，因为ntens=3,最多只能计算DDSDDE(3,3))

　你编umat编进去不久行了

　他让用NDI,NSHR,NTENS表示变量，你就用这些表示变量

　这样他就可以任意的计算了呀，而不在意实际计算的变量数

　: 唉，别提了，问题就出在，在abaqus中，明明写着可以考虑剪切效应，

　: 可我打印出剪切力个数是，nshr=1,即只有S12,

　: 那我的S13,S23就不知怎么计算，

　: （DDSDDE(5,5)无法计算，因为ntens=3,最多只能计算DDSDDE(3,3))

　因为我的UMAT从abaqus传来的变量（当我选用壳单元时）ntens=3,

　ndi=2,nshr=1,这样的话，STRESS为3个，STRESS(3),不可能计算

　STRESS(4),STRESS(5),同理，STRAIN,DDSDDE也存在同样的情况

　在中厚度板壳元中(MINDLIN)中，DDSDDE为5×5

　你的FOR文件收到，可惜研究的对象相差太大，看不懂！

　不过你在文件中定义的变量有的根本就没有用到，

　还有有的地方似乎是毫无用处的代码，比如：SMISES一段

　我觉得你不应该太急，第一步应该保证代码的正确性。

　另外看了你上面的讨论，我觉得你应该换一种单元试一试。

　比如用梁单元。很有可能NDI和NSHR会不一样。

　不知道你是不是一定要用到壳单元上？

　如果必须用在壳单元上，可能S23，S31就必须自己在UAMT

　中输出了。具体的做法我也不清楚。

　还有对DDSDDE，我的理解也是: dS=DDSDDE*dE

　我的意思是不管他是几乘几的矩阵

　一般的D矩阵就是6*6

　你编的umat需要适应于任何单元

　所以当然6*6的矩阵应该给完全

　当算体的时候

　如果是壳单元，他只有ntens=3

　那就等于三呗，他用哪个就是哪个

　反正你的D6*6的矩阵放在哪了

　如果是三维实体单元，ntens=6

　所以你的D矩阵就全用上啦

　记住你编的umat应该适应于所有的单元，所以不用考虑哪几个参数的值

　而且这个stress(3)也许是人家只用到了你D6*6的第一行，第二行，和的四行哪

　sigh

　不知道我说明白了没有

　: 因为我的UMAT从abaqus传来的变量（当我选用壳单元时）ntens=3,

　: ndi=2,nshr=1,这样的话，STRESS为3个，STRESS(3),不可能计算

　: STRESS(4),STRESS(5),同理，STRAIN,DDSDDE也存在同样的情况

　: 在中厚度板壳元中(MINDLIN)中，DDSDDE为5×5

　对于DDSDDE，其实就是给了一个应力应变的斜率，我感觉

　真正计算的是利用更新的应力，DDSDDE我觉得就是为了确保收敛

　真开心有这么多人讨论ABAQUS，特别是umat

　我国外的师兄说Umaterial在国外应用得也很广

　: 对于DDSDDE，其实就是给了一个应力应变的斜率，我感觉

　: 真正计算的是利用更新的应力，DDSDDE我觉得就是为了确保收敛

　: 真开心有这么多人讨论ABAQUS，特别是umat

　: 我国外的师兄说Umaterial在国外应用得也很广

　各位大侠，看到我提出的uamt,及statev有这么多人感兴趣，实在很

　高兴，希望大家能把摸索到的东西介绍给大家，现在我的程序能通了，

　（因没有屈服面的概念),虽然结果是错误的，但总可以静下心来，不像以前

　一运行自己的umat程序就standard.exe error.(并且程序.dat.msg.中无

　任何错误信息，log说error see msg,but msg中无）。

　现在我明白了，在用户程序中，WRITE(*,*)则写入log,write(6,*)则写入

　dat,等，

　在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，

　还是为上一步的值，因为我每次都为0，怀疑有误

　还望大侠提醒！

　: dat,等，

　: 在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，

　应该不是0，如果在每个增量步的最后你重新更新过他

　输入文件*.inp中，*depvar不要忘了写statev()的维数哦！

　: 还是为上一步的值，因为我每次都为0，怀疑有误

　: 还望大侠提醒！

　你只有一个单元么？

　【在 saturn (好好工作) 的大作中提到: 】

　: ~!~~~~~~~~~~~~~~~~~~~~~~~~~~~~这我以前还不知道

　: 应该不是0，如果在每个增量步的最后你重新更新过他

　: 输入文件*.inp中，*depvar不要忘了写statev()的维数哦！

　【在 air1 (无影如风※勤奋羊羊※课题在即) 的大作中提到: 】

　: 可是这个写的命令有用么

　: 对于输出而言

　: 你只有一个单元么？

　已经写了，定义了statev的维数，

　*depvar

　20,

　看来我的错误不在statev,

　因为我的输出显示在第三增量步中，statev不为0，

　但我还是有些不清楚uamt程序的执行过程，

　比如，

　一板一次壳元s4r，面内采用缩减积分，即1×1，厚度方向为5点simpson

　积分，按道理，

　没一增量步的迭代步，对每一积分点都会运行一次umat,

　即上问题，每一迭代该运行5次，

　那上面的statev,是每次迭代都会改变（不是每次运行）？

　请air1,beauty_life,saturn,......等各位相助

　看过书上例子，

　如果定义statev为（nblock，*）的数组，

　应该对每个积分点都进行迭代计算，从而更新变量值

　（写stateNew(nblock,*)

　请问Mikejwg，变量statev怎么输出？

　可以象应力应变那样进行绘图、动画等后处理吗？

　不好意思，我只在log文件中输出，还没有在odb中输出，

　想同element output等那样输出，只能在UVARM中定义，但

　同statev是不同的，老兄你可以试试。

　发信人: Mikejwg (mikejwg), 信区: FEA

　标题: 用abaqus中的UMAT的大侠看过来。。。

　发信站: BBS 水木清华站 (Tue Dec 25 15:36:27 2001)

　各位大侠，

　在加入自己弹塑性本构模型时，如何加入区分弹性部分与塑性部分（或

　非弹性部分）。也就是说，在输出时，EE,LE,IE([PE])，因为在umat中，

　我们只发现dstran,stress,除此并没有发现其它有关应力应变得变量。

　你所说的区分弹性部分和塑性部分是指应变的弹塑性部分：E=EE+PE吗？

　还是指加载过程的弹塑性？

　E=EE+PE

　我觉得想从umat的输入参数里面区分是不可能的

　可不可以你先算出弹性的，然后用stran减去弹性的，得到塑性的

　我猜的哦

　我自己是明白的，但不知怎么让程序明白

　比如，在一定条件下可解出位移增量，abaqus得到应变增量dstran,stran,

　而弹性部分，dEEij=dSij/(2G),dEPij=dstranij-dEEij

　但怎么控制程序了？

　默认得情况又是怎样计算的

　谢谢

　我不清楚你的本构，

　我大概觉得你如果是用屈服面的话

　你是不是可以用statev（）记录每一步的EE，EP

　然后通过应力偏量来计算DEE,DEP，然后自己更新EE，EP，就不与stran扯上关系

　我以前见过用Umat算J2流动理论的时候，是这么做的

　首先我的本构不能采用屈服面的概念（一种新型的弹塑性本构模型），

　我也看了那个采用j2流动法则的例子，我可以在umat中的statev重定义

　ee,pe,但statev不能像EE,PE,IE那样在odb中输出。

　再说在umat中的ee,pe与inp中的ee,pe是不同的。

　谢谢

　看来你只能试试uvarm了

　呵呵，

　发信人: Mikejwg (mikejwg), 信区: FEA

　标题: Re: ABAQUS中的弹塑性应变的计算问题？？

　发信站: BBS 水木清华站 (Thu Dec 27 14:19:27 2001)

　用uvar()勉强成功

　标题: 在abaqus中，DSTRAN,STRAN为工程剪应变吗？

　Sat Dec 29 09:44:03 2001)

　同题，这涉及到ddsdde,peeq的计算

　对，我也是最近才发现的

　peeq怎么计算？

　在umat中计算peeq,(因用umat后，系统缺省peeq=0)然后再uvarm输出

　uvarm可以输出到哪里？.odb可以么？

　另外那个strav？？就是自己定义用于umat的那个数组里的数可以输出到.odb里么？

　谢了，你这个成长中的大侠

　在umat中，statev是不能用在output中的，statev只是作为一个解的

　状态变量，说来惭愧，我是在umat中定义peeq,(peeq）的求解一般在

　弹塑性力学书上有（等效塑性应变），用write()写入一个临时文件，

　（（切记：这个文件unit号不要与abaqus中的重合，因为他有一些系统默认

　的文件号，））

　然后在uvarm中读取，以uvarm输出，因为uvarm可以以odb的形式输出，

　支持output,field,output,history

　不知谁还有更好的方法。

　听说清华有人abaqus用的很熟了。

　大家过奖了，我是没办法，要出论文，总得学习，上午版主还在qq中指点了我一番

　YOU CAN TRY AS FOLLOWS:

　*STEP,EXTRAPOLATION=NO,INC=2000000

　*STATIC

　0.001,500.0,0.00001,0.1

　可以

　并且在step里好像有选项选择

　你是手写的程序么

　新建一个step

　然后定义边条

　在上一step定义边条的节点上重新定义边条

　而且，好像，是不是有个选项是new

　这样就可以取消上一step的边条设置了

　在以后的载荷步里重新定义边界条件，加参数new

　比如在第二个载荷步里：

　 *STEP

　......

　*BOUNDARY CONDITIONS,OP=NEW

　......

　*END STEP

　Command方式：abaqus job=job_name user=source_file（或者obj_file);

　CAE方式：在最后submit的时候可以选择source_file（或者obj_file).

　不太明白你的意思。

　如果只是要得到所需要的如应力，应变等，那么在数据行请求

　输出就行了啊，而且个人觉得abaqus的后处理功能基本上也够用了，直接输出到*.odb文件

　中，用CAE或者Viewer后处理即可。

　而*el print输出的dat文件包含其他的信息，如版本号、处理内容等。

　若想要直接输出到其他后处理软件的话，可以输出到result file文件*.fil，不过，

　这可是二进制文件或者ASCII文件，还是需要处理的，或者直接用Hypermesh等后处理

　程序打开做图。

　好像是在tool->display->create－>再选定一个单元集或者是节点集

　然后操作是“取代”还是其他等等（就是里面的两个圆圈相交，

　还是想减）。你试试看看。

　这个是要到后处理才行吧

　非要计算以后形成*.odb吗

　前处理里面怎么才可以啊

　 Python在abaqus中可以运行scrpit文件，你在CAE manual中就

　可以将其附录中的例子存为*.py文件，然后在CAE的菜单中

　run script就可以运行了。

　 BTW:请问你的是abaqus的那个版本？6.2 or 6.3？6.3的有盗版了？

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

　试一试：*element matrix output

　: (2),想要用Newmark方法求解.求解问题,在ABAQUS中如何实现?或者

　: 那里能找到类似的例子?

　看看：standard manuals:analysis的dynamics analysis和explicit manuals

　: (3),ABAQUS中能实现这样的东西吗?就是我需要平滑ABAQUS产生的位移场

　: 还是这个平滑只能在其他环境中完成?

　试试这个输出选项：

　 *El print,Position=average at Nodes

　: abaqus/CAE好像不能生成节点集，只能定义某个实体，后处理用这个实体上的节点。

　 : 看到CAE的例题都是这样做的。

　 NODE SET是abaqus的基本功能，也是其方便使用之处，

　在CAE中是支持的。

　这个例子我做过，没有问题。

　 : : 我用tool》set》create居然不能选择由CAE生成的节点（不是顶点）

　这是因为你没有正确操作Partition Edge这一步。

　【在 sunhaolan 的大作中提到: 】

　: abaqus/CAE好像不能生成节点集，只能定义某个实体，后处理用这个实体上的节点。

　 ABAQUS divides all of its analysis procedures into

　two main groups: linear perturbation and general.

　ABAQUS makes a very clear distinction between

　perturbation and general analysis procedures because

　loading conditions and “time“ are defined differently

　for the two cases. Furthermore, the results from each

　type of procedure should be interpreted differently.

　The response of the model during a general analysis

　procedure, known as a general step, may be either

　nonlinear or linear. In a step that uses a perturbation

　procedure, which is called a perturbation step,

　the response can only be linear. ABAQUS treats

　such steps as a linear perturbation about the preloaded,

　predeformed state (known as the base state) created by

　any previous general steps; therefore, its capability

　for doing linear simulations is rather more general

　than that of a purely linear

　analysis program.

　各位大侠，我在定义step的时候，前三个用的是standard中的相关命令.在

　step 4的时候，我想用explicit来求解，并且定义了amplitude,来模拟一个

　位移随时间的简单变化，命令如下：

　 *step name=“～～～“，nlgeom,inc=1000(好像explicit里没有inc这个参数）

　*dynamic,explicit

　*amplitude name=ramp

　0,0,1,1,2,1

　*boundary amplitude=ramp

　refpunch(刚体参照点）2，2，-0.03

　*end step

　有问题吗？

　另外standard的命令和explicit的命令能出现在同一个.inp文件中吗？

　要转换的.

　ABAQUS/Standard and ABAQUS/Explicit are separate program modules with

　different data structures; hence, the

　explicit dynamics procedure cannot be used in the same analysis as any of

　the procedures in ABAQUS/Standard.

　However, ABAQUS provides a capability to import a deformed mesh and associated

　 material state from

　ABAQUS/Explicit into ABAQUS/Standard and vice versa. This procedure is

　described in ``Transferring results between

　ABAQUS/Explicit and ABAQUS/Standard,＇＇ Section 7.6.2 of the ABAQUS/Standard

　User＇s Manual.

　【在 linyuwho (苍龙有悔) 的大作中提到: 】

　: 各位大侠，我在定义step的时候，前三个用的是standard中的相关命令.在

　: step 4的时候，我想用explicit来求解，并且定义了amplitude,来模拟一个

　: 位移随时间的简单变化，命令如下：

　 : *step name=“～～～“，nlgeom,inc=1000(好像explicit里没有inc这个参数）

　: *dynamic,explicit

　: ,2

　: *amplitude name=ramp

　: 0,0,1,1,2,1

　: *boundary amplitude=ramp

　: refpunch(刚体参照点）2，2，-0.03

　: *end step

　: ...................

　feature菜单下有删除一项。基准点是一种feature

　: 工具栏不有个像铅笔一样的的东东么

　我想把计算出来的应力，应变，位移等数据提取出来，在其他软件里处理，

　请问该怎么办？我在。dat里没有找到这些数据阿？。res不是文本格式，

　其他的好像也没有啊。

　计算结果，没有文本格式，只能自己把数据倒出来。

　在显示结果的模块中，名字忘了，

　就是sketch前面那个，

　在report菜单中，可以到处自己需要的文件。

　文件格式是：rpt

　谢谢，但是这个为什么只有应力和应变呢？计算出的位移的数据怎么会没有呢？

　【在 santong (supernova) 的大作中提到: 】

　: 计算结果，没有文本格式，只能自己把数据倒出来。

　 : 在显示结果的模块中，名字忘了，

　: 就是sketch前面那个，

　: 在report菜单中，可以到处自己需要的文件。

　 : 文件格式是：rpt

　发信人: wjytj (赖子), 信区: FEA

　你要在STEP模块field output和history output改变默认

　请求输出位移U

　或者在*step中加上：

　 *NODE PRINT

　计算完成后在*.DAT文件中就可以提取节点位移。

　请问field output和history output什么区别？

　关键字*node output和*node print有什么区别？谢谢！

　我的。inp数据文件里有这么几行：

　 *Output, field, frequency=99999

　*Node Output

　U, RF

　*Element Output

　S, E, EE

　*Output, history, variable=PRESELECT

　*El Print, freq=999999

　*Node Print, freq=999999

　其中在history output中有*node print关键字，但是没有指定输入U，

　最后导出的数据没有U这一选项，请问是什么原因？

　【在 wjytj (赖子) 的大作中提到: 】

　: 你要在STEP模块field output和history output改变默认

　: 请求输出位移U

　: 或者在*step中加上：

　 : *NODE PRINT

　: U

　: 计算完成后在*.DAT文件中就可以提取节点位移。

　发信人: Mikejwg (mikejwg), 信区: FEA

　field output和history output这在viewer模块中很明显。

　 field output是场量输出，history output为历史数据输出，会记录场量随时间的曲线

　至于*node output与*node print的区别在于他们写入不同的后缀文件，

　比如*.odb,*.dat等，在help中有详细介绍。

　【在 wsun 的大作中提到: 】

　: 请问field output和history output什么区别？

　: 关键字*node output和*node print有什么区别？谢谢！

　: 我的。inp数据文件里有这么几行：

　 : *Output, field, frequency=99999

　: *Node Output

　: U, RF

　: *Element Output

　: S, E, EE

　: *Output, history, variable=PRESELECT

　: *El Print, freq=999999

　: ...................

　发信人: wjytj (赖子), 信区: FEA

　然也！！补充一句：

　 *node output是给CAR或VIEWER做后处理用的，是二进制文件；

　而*node print则是写到DAT文件中，你可以自己看的文本文件。

　【在 Mikejwg 的大作中提到: 】

　: field output和history output这在viewer模块中很明显。

　 : field output是场量输出，history output为历史数据输出，会记录场量随时间的曲..

　 : 至于*node output与*node print的区别在于他们写入不同的后缀文件，

　: 比如*.odb,*.dat等，在help中有详细介绍。

　 : 【在 wsun 的大作中提到: 】

　: : 请问field output和history output什么区别？

　: : 关键字*node output和*node print有什么区别？谢谢！

　: : 我的。inp数据文件里有这么几行：

　 : : *Output, field, frequency=99999

　: : *Node Output

　: ...................

　发信人: octopuszy (猴哥), 信区: FEA

　*node output是在odb文件中输出?br /> *node print是在dat文件中输出

　【在 Mikejwg (mikejwg) 的大作中提到: 】

　: field output和history output这在viewer模块中很明显。

　 : field output是场量输出，history output为历史数据输出，会记录场量随时间的曲..

　 : 至于*node output与*node print的区别在于他们写入不同的后缀文件，

　: 比如*.odb,*.dat等，在help中有详细介绍。

　 file>print>file

　文件类型选择：png或tiff

　【在 Dshmy 的大作中提到: 】

　: 谢谢

　一般说可以这样理解：

　 Δσ＝ddsdde*Δε

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

　【在 wjytj 的大作中提到: 】

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

　: DDSDDE被称为：

　 : “material stiffness matrix“;

　: 而在UMAT中，DDSDDE被称为：

　 : “material＇s Jacobian matrix“。

　 : 请问DDSDDE和材料的切线刚度矩阵的关系是什么？

　: 二者是一个概念么？

　第二个

　【在 octopuszy (猴哥) 的大作中提到: 】

　: 偶在umat中调用求主应力函数

　: CALL SPRINC(STRESS,PS,LSTR,NDI,NSHR)

　: 后，存储主应力得数组PS中

　: 各个主应力排列顺序是什么？

　: PS1>PS2>PS3 ?

　: PS1<PS2<PS3 ?

　: PS1>PS3>PS2 ?

　: 谢谢！

　 abaqus job=jobname parametercheck

　abaqus job=jobname datacheck

　abaqus job=jobname continue

　【在 gentlewtj (gentlewtj) 的大作中提到: 】

　: 那有没有检查错误的方法？

　【在 uptonow (到目前为止) 的大作中提到: 】

　: 1、在边界条件和加载时，总是有initial这个步，然后是我们自己定义的

　: 加载步，请问这个initial步，主要作用是什么？能不能去掉？

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

　: 2、A solid extrusion base feature 这句话是什么意思？

　: extrusion、revolution、等是什么意思？

　这两的是三维建模时候,在画完二维图形,如何来生成三维图形,

　extrusion意思是你给定一个厚度,然后二维图形第三个方向上面伸展这么多形成三维图形

　 revolution意思是你给定一个旋转轴,二维图形绕其旋转后形成三维轴对称图形

　: 谢谢！

　发信人: fangxj (fxj), 信区: FEA

　标题: Re: abaqus中接触问题请教

　发信站: BBS 水木清华站 (Fri Jun 14 15:51:20 2002)

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

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

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

　【在 octopuszy (猴哥) 的大作中提到: 】

　: 偶的模型中存在两个物体的接触，计算过程中报错：

　 : ***WARNING: OVERCLOSURE OF CONTACT SURFACES FOUNDATIONOUTSIDE and

　: BUCKETOUTSIDE IS TOO SEVERE -- CUTBACK WILL RESULT. YOU MAY WANT

　: TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION.

　: 。。。。。。。。。。。。。。。。。。。。。。。。。

　 : ***WARNING: CONVERGENCE JUDGED UNLIKELY. INCREMENT WILL BE ATTEMPTED AGAIN

　: WITH A TIME INCREMENT OF 9.76563E-05

　: ***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT: ANALYSIS TERMINATED

　: 不知该如何解决？

　standard使用隐式作动力分析，explicit用显式作动力分析，你该用explicit。

　【在 uptonow (到目前为止) 的大作中提到: 】

　: 请问这两本书有什么不同，

　: 我做的是圆形杯子的深冲的模拟过程，我这里只有standard一书

　: 那位有explicit书的大虾请留个言，我想复印，谢谢：）

　一般说来，你那样是正确的，如果不行，你可以在

　dat,log,msg中找到中断的增量步，假如为60，不妨可以

　在res_pulse中，

　*restart,read,step=1,inc=59

　【在 ypm (红塔山) 的大作中提到: 】

　: 我的一个abaqus程序(结构动力学,用的/standard模块,假设名

　: 为pulse.inp,一个step, inc=1000, restart命令的参数为

　: *restart, write, frequency=10),在算了一半之后人为停掉,然

　: 后我想从中止的地方继续运算.

　: 我定义一个res_pulse.inp,其中就一个命令

　: *restart,read

　: 然后在command窗口,输入 abaqus job=res_pulse oldjob=pulse

　: 结果程序从起始点开始运算,而不是从中断的地方.请问是什么

　: 原因? 刚学不久,希望大侠们指点指点!thanks

　发信人: Mikejwg (mikejwg), 信区: FEA

　uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难一个一个

　出现，永远有解不完的问题，就本人体会，写一点东西：

　在umat中，

　(我总结三点：

　 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵）

　2，据应力势函数(屈服函数）求解满足条件的应力，应变，

　在此我们要注意试算应力与真实应力的关系：sigma_true=Sigma_tri-De*dpstran: dpstr

　 an为塑性应变增量，据此满足屈服条件，F（Sigmaij）＝0，求得真实的应力

　更新应力 stress=stress-De*dpstran

　3，（与2联合，弹塑性区域的转换）判断是否屈服，

　4，组集ddsdde(与个人的本构有关）

　这4步都要在UMAT实现

　采用上一迭代步构造ddsdde,然后dsigma=ddsdde*dstran不可取，以我个人经验，收敛很慢

　并且敏感与步长。

　希望大家交流！

　发信人: uptonow (为什么受伤的总是我?), 信区: FEA

　衷心感谢Mikejwg大侠的该出手时，就出手！希望Mikejwg大侠以后经常出手！

　造福于UMAT的万民！　

　【在 Mikejwg (mikejwg) 的大作中提到: 】

　: uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难一?.

　 : 出现，永远有解不完的问题，就本人体会，写一点东西：

　 : 在umat中，

　: (我总结三点：

　 : 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵）

　: 2，据应力势函数(屈服函数）求解满足条件的应力，应变，

　: 在此我们要注意试算应力与真实应力的关系：sigma_true=Sigma_tri-De*dpstran: d..

　 : an为塑性应变增量，据此满足屈服条件，F（Sigmaij）＝0，求得真实的应力

　: 更新应力 stress=stress-De*dpstran

　: 3，（与2联合，弹塑性区域的转换）判断是否屈服，

　: 4，组集ddsdde(与个人的本构有关）

　: ...................

　发信人: wjytj (赖子), 信区: FEA

　精华区有一些帖子，很有帮助，大家感兴趣可以去看看

　【在 Mikejwg 的大作中提到: 】

　: uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难一?..

　 : 出现，永远有解不完的问题，就本人体会，写一点东西：

　 : 在umat中，

　: (我总结三点：

　 : 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵）

　: 2，据应力势函数(屈服函数）求解满足条件的应力，应变，

　: 在此我们要注意试算应力与真实应力的关系：sigma_true=Sigma_tri-De*dpstran: d..

　 : an为塑性应变增量，据此满足屈服条件，F（Sigmaij）＝0，求得真实的应力

　: 更新应力 stress=stress-De*dpstran

　: 3，（与2联合，弹塑性区域的转换）判断是否屈服，

　: ...................

　发信人: uptonow (为什么受伤的总是我?), 信区: FEA

　精华区的帖子太少了，还需要你补充呢！不要推辞了！：）

　推辞就是逃避责任呵呵，开玩笑的

　【在 wjytj (赖子) 的大作中提到: 】

　: 精华区有一些帖子，很有帮助，大家感兴趣可以去看看

　发信人: edodo (大漠孤沙), 信区: FEA

　老哥可以挺身而出呀！哈哈！

　期待中！

　【在 uptonow (为什么受伤的总是我?) 的大作中提到: 】

　: 精华区的帖子太少了，还需要你补充呢！不要推辞了！：）

　: 推辞就是逃避责任呵呵，开玩笑的

　发信人: air1 (无影如风※勤奋36小时※永不言败), 信区: FEA

　我曾经犯的错误

　1,发现456方向的应变是工程应变

　2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵.

　3,ddsdde何时更新对于收敛性很有影响.

　【在 Mikejwg (mikejwg) 的大作中提到: 】

　: uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难一?.

　 : 出现，永远有解不完的问题，就本人体会，写一点东西：

　 : 在umat中，

　: (我总结三点：

　 : 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵）

　: 2，据应力势函数(屈服函数）求解满足条件的应力，应变，

　: 在此我们要注意试算应力与真实应力的关系：sigma_true=Sigma_tri-De*dpstran: d..

　 : an为塑性应变增量，据此满足屈服条件，F（Sigmaij）＝0，求得真实的应力

　: 更新应力 stress=stress-De*dpstran

　: 3，（与2联合，弹塑性区域的转换）判断是否屈服，

　: 4，组集ddsdde(与个人的本构有关）

　: ...................

　发信人: wjytj (赖子), 信区: FEA

　: 1,发现456方向的应变是工程应变

　说白了，因为在UMAT中应变是以列向量的形式存储的，而初始刚度DDSDDE的形式为：

　 λ＋μ μ μ 0 0 0

　 μ λ＋μ μ 0 0 0

　 μ μ λ＋μ 0 0 0

　 0 0 0 μ 0 0

　 0 0 0 0 μ 0

　 0 0 0 0 0 μ

　所以，应变是工程应变，而如果DDSDDE的形式为（直接用张量表示）：

　 λ＋μ μ μ 0 0 0

　 μ λ＋μ μ 0 0 0

　 μ μ λ＋μ 0 0 0

　 0 0 0 2μ 0 0

　 0 0 0 0 2μ 0

　 0 0 0 0 0 2μ

　那么，应变就应该是应变张量中的应变了，也就是说，采用工程应变只是为了方便应用而

　已，比如，求弹性应变能，可以直接将应力列向量的转置乘以应变列向量，显然，这在编

　程中很容易实现。

　 : 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵.

　同感，应该可以减少不少工作量。

　 : 3,ddsdde何时更新对于收敛性很有影响.

　希望air1能够说得更详细些，比如在开始是可以采用初始刚度，在后来可以采用切线刚度

　，可是如何控制，编程如何实现呢？

　【在 air1 的大作中提到: 】

　: 我曾经犯的错误

　: 1,发现456方向的应变是工程应变

　: 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵.

　: 3,ddsdde何时更新对于收敛性很有影响.

　: 【在 Mikejwg (mikejwg) 的大作中提到: 】

　: : uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难?..

　 : : 出现，永远有解不完的问题，就本人体会，写一点东西：

　 : : 在umat中，

　: : (我总结三点：

　 : : 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵）

　: ...................

　发信人: wjytj (赖子), 信区: FEA

　: 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵.

　在我的本构中，DDSDDE的确定需要知道STRESS中的某些量，那么

　请问，按照你的意思和经验，在UMAT中是否可以先得到STRESS(比如，

　利用全量法），然后再得到DDSDDE？

　如果可以的话，就应该可以减少迭代求解的次数了。

　谢谢！！

　假定采用mises屈服准则。

　 1 在调用UMAT之前，ABAQUS传递给UMAT本次增量开始时的应力sigma(0)，总应变E,应变增

　量delta(E)。

　状态变量保存：弹性应变，塑性应变，等效塑性应变。

　（请问，状态变量保存的弹性应变＋塑性应变是否等于ABAQUS传递给UMAT的总应变？？？

　？？）

　2，然后在UMAT中利用上述的ABAQUS传递的量和状态变量得到DDSDDE矩阵，然后返回给

　ABAQUS，ABAQUS根据delta(sigma)=ddsdde*delta(E)，并且得到本次增量结束时的应力

　sigma(1)=sigma(0)+delta(sigma)=sigma(0)+ddsdde*delta(E)

　3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下次调

　用UMAT

　请问手册上UMAT必须更新应力，可是根据上述我的理解好像是ABAQUS根据UMAT提供的本

　次增量的DDSDDE在ABAQUS中更新，请问到底是怎么回事？谢谢！

　个人觉得，在Abaqus中更新应力是必须的，而且必须准确，这是你的本构的核心；DDSDDE

　则可以是不准确的，比如，你可以保持DDSDDE一直为初始刚度矩阵，当然，这样做的问题

　会是收敛速度较慢或者根本不收敛。

　在Manual中也说过，DDSDDE不影响结果，前提是，只要能得到结果。

　与大家共同讨论。

　【在 uptonow 的大作中提到: 】

　: 假定采用mises屈服准则。

　 : 1 在调用UMAT之前，ABAQUS传递给UMAT本次增量开始时的应力sigma(0)，总应变E,应..

　 : 量delta(E)。

　 : 状态变量保存：弹性应变，塑性应变，等效塑性应变。

　 : （请问，状态变量保存的弹性应变＋塑性应变是否等于ABAQUS传递给UMAT的总应变?.

　: ？？）

　: 2，然后在UMAT中利用上述的ABAQUS传递的量和状态变量得到DDSDDE矩阵，然后返回?.

　 : ABAQUS，ABAQUS根据delta(sigma)=ddsdde*delta(E)，并且得到本次增量结束时的应..

　 : sigma(1)=sigma(0)+delta(sigma)=sigma(0)+ddsdde*delta(E)

　: 3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下..

　 : ...................

　期盼了好久的abaqus6.3终于安装成功

　寄希望于此来成功import我的.igs模型（unigraphics v18生成的）

　但是，结果却不尽如人意，现将我的感受说一下：

　 1. 6.3在import进模型后，它的形状基本上与真实形状差不多,比6.2更像原始模型

　但是，依然是invalid geometry

　2. 经过delete invalid geometry后，它确实成了imprecise geometry(意味着你

　以划分网格了），但是此时的形状已经面目全非了，删除了太多的东西了，得到的

　imprecise geometry 就毫无意义了

　3. 如果我用automated repair，它还是invalid geometry

　我现在急需把文件import进后来进行下一步工作，

　而老板出国前定的任务是11月底出结果

　我现在连模型都是个问题了

　急急急！！！

　请大侠帮忙！！！！！

　救救我吧！给个好的建议吧

　用PATRAN导入，生成网格后用ABAQUS计算

　试试parasolid的文件格式，iges格式先天不足，尤其是表达曲面。

　对了俺不知道abaqus支持parasolid文件不？不过一般情况下比igs格式

　好点。

　 State variables

　Since the redefinition of field variables in USDFLD is local to the current

　increment (field variables are restored to the values

　interpolated from the nodal values at the start of each increment),

　any history dependence required to update material properties by

　using this subroutine must be introduced with user-defined state variables.

　The state variables can be updated in USDFLD and then passed into other user

　subroutines that can be called at this material point,

　such as those listed above. The number of such state variables is

　specified with the *DEPVAR option...

　要引入材料的非线性，需要usdfld;而要引入历史相关，就得要靠statev帮忙了！

　(我的理解不全面，可以自己参考manual中的这一段.版主g上吧，我以后忘了看文摘区

　就行了

　感谢！

　现在弄明白了，原来我的文件里多了几行空行，没想到abaqus的inp文件不允许有空行的。

　感谢wjytj的帮助！刚开始用abaqus，太弱了。

　【在 wjytj 的大作中提到: 】

　: :*NGEN

　: :1,111,11

　: :11,121,11

　: **此时，12，22，23……110节点并未生成，所以生成12，22以后的节点是没有意义?..

　 : 看看这样行否？

　: *NODE

　: 1,0,10

　: 11,10,10

　: 111,0,0

　: 121,10,0

　: ...................

　: solution-dependent variable

　从字面上理解，即为“与解答有关的变量”，它的值与每一增量结束时的变量有关，常用

　于UMAT中，即SDV，可以由*DEPDV来定义其个数。

　 : time-dependent variable

　从字面上理解，即为“与时间有关的变量”，常用于定义随时间变化的量如地震荷载等，

　，可以由*APPLITUDE来定义。

　【在 octopuszy 的大作中提到: 】

　: solution-dependent variable

　: time-dependent variable

　: 这两种变量有什么差别？

　: 谢谢！

　【在 octopuszy 的大作中提到: 】

　: UFIELD和USDFLD定义的场量有什么不同

　: 另外solution-dependent variable

　: 和time-dependent variable

　: 有什么不同？

　: 谢谢！

　标题: 给大家几点关于UMAT的建议！

　今天UMAT调试通过，非常高兴！！

　就我的经验来说，给大家几点建议：

　 1、UMAT一般应该适合各种应力状态，也就是无论NDI＝2或3

　都应该适用；

　2、UMAT中应力和应变列向量的存储方式（下标）是：

　 11，22，33，12，13，23，而不是11，22，33，12，23，31；

　3、ABAQUS第一次计算中传递给UMAT的应变STRAN和应变增量

　DSTRAN的值均为零（因为需要得到DDSDDE），此时UMAT中的

　应力也应该为零，所以，UMAT中分母有可能为零，这种情况

　应该注意！！（我就是这样陷在里面好几天）

　17.请问怎么实现双曲线

　subroutine fric(lm,tau,ddtddg,ddtddp,dslip,sed,spd,

　 1 ddtddt,pnewdt,statev,dgam,taulm,press,dpress,ddpddh,

　 2 slip,kstep,kinc,time,dtime,noel,ciname,slname,

　 3 msname,npt,node,npatch,coords,rcoord,drot,temp,

　 4 predef,nfdir,mcrd,npred,nstatv,chrlngth,props,nprops)

　 include ＇aba_param.inc＇

　 character*80 ciname,slname,msname

　 dimension tau(nfdir),ddtddg(nfdir,nfdir),ddtddp(nfdir),

　 dslip(nfdir),ddtddt(nfdir,2),statev(*),

　 dgam(nfdir),taulm(nfdir),slip(nfdir),time(2),

　 coords(mcrd),rcoord(mcrd),drot(2,2),temp(2),

　 predef(2,*),props(nprops)

　 parameter(zero=0.0D0,two=2.0D0)

　 if (lm .eq. 2) return

　 a=props(1)

　 b=props(2)

　 lm=0

　 tau(1)=dslip(1)/(a+b*dslip(1))

　 tau(2)=dslip(2)/(a+b*dslip(2))

　 ddtddg(1,1)=tau(1)/slip(1)

　 ddtddg(2,2)=tau(2)/slip(2)

　 ddtddg(1,2)=zero

　 ddtddg(2,1)=zero

　 dslip(1)=dgam(1)

　dslip(2)=dgam(2)

　 return

　 end

　18.取消坐标系等的显示

　CAE进行后处理画的等值线图挺不错，只是如果通过Print to File得到的图形质量很差，远远不如屏幕截图的效果。但通过屏幕截图也有个问题，就是坐标系——特别是有关ODB中当前step，increment等不需要的信息也一起给截来了。请问如何关闭这些信息？

　Module Visualization-Canvas-Viewport Annotation Options

　19.如何在计算中修改材料特性

　If you like to use USDFLD, in the subroutine, you can get step number and just use a if to

　 determine the field variable value due to the step number. Say, step 1 is 1 and step 2 field

　 variable become 2. In the material definition, just define material properties as a function of

　 your field variables. Here I use elastic as example. I＇m not sure about the grammar since I do

　 not have manual at hand. You need to check them. something like this:

　 *elastic, dependencies

　 2e5,0.3,,1

　 1e5,0.3,,2

　 Here I assume your material does not have any temperature dependency and you need to make

　 sure field variable is the 4th value. Then, in step 1, you have 2e5 as Young＇s modulus and in step

　 you have 1e5.

　 yogayoga，我还是有问题，我修改材料特性后，计算不收敛。请您帮我看看。

　在inp文件中，按你讲的使用了下面语句定义材料特性：

　*elastic, dependencies＝1

　 15e7,0.3,,1

　 2e10,0.2,,2

　然后采用usdfld子程序，现将子程序附在下面请您帮我检查检查，并提出宝贵意见，谢谢。

　SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,

　1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,

　2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)

　 INCLUDE ＇ABA_PARAM.INC＇

　 CHARACTER*80 CMNAME,ORNAME

　 CHARACTER*3 FLGRAY(15)

　 DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),

　1 T(3,3),TIME(2)

　 DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)

　 FIELD(1)=time(2)

　print*,field(1)

　 RETURN

　 END

　 First, if you want a sudden change, you can not use total time. Abaqus will linearly interpolate

　 all the value. If that＇s what you want, OK. If you want a sudden change at step two, use kstep

　 instead.

　 Second, consider the real physical situation: you got equilibrium before you changed material

　 properties. After you changed material properties, you naturally lose it. Stress need to be

　 redistributed and changed to get force balance. This is the reason for non convergence.

　 You have several way to deal with it: if you have static problem and rate effect is minimal, you

　 can keep the bc and let the model balance itself. The other way is to decrease the increment

　 length. You can also try damping (you can just apply stabilize directly) or adding some viscoplastic

　 or plastic effect. Dynamic analysis will also help to do the trick. Basicly, the method normally you

　 deal with unstable problem.

　改变材料的参数甚至模型，可以通过以下的办法：

　（1）在不同的分析之间传递结果

　细节请见：ABAQUS Analysis User＇s Manual /7.7.1 Transferring results between ABAQUS

　 analysis products: overview

　该方法的主要缺点就是建立模型时不能使用多个parts或assembly。

　但在某些情况下非常方便，如想在第一步变形的基础上做后续分析。

　（2）使用场变量，然后定义材料时指明参数依赖于该变量

　如果场变量（不一定要有物理意义）很明确，比如目前讨论的问题，可以定义一个变量，在

　step1的时候为1.，step2为2.(可以随便定义了)：

　 *STEP，name=step1

　*STATIC......

　*FIELD, VARIABLE=1

　NALL, 1.

　*END STEP

　*STEP，name=step2

　*STATIC......

　*FIELD, VARIABLE=1

　NALL, 2.

　*END STEP

　NALL是包含所有节点的集合，可以根据需要改变成欲改变参数的节点集。

　然后再定义材料时，指明某个参数依赖于该场变量：

　 *MATERIAL

　*ELASTIC, DEPENDENCIES=1

　200.E9, 0.3, , 1.

　 180.E9, 0.3, , 2.

　这样计算到step1时，所有节点的场变量1为1.0，abaqus由此查表确定模量为200.E9，到step2，

　因为场变量变为2.，所以模量为180.E9。

　如果场变量连续变化或依赖于其他计算结果，如温度，应变等，可以使用USDFLD子程序（25.2.38

　 USDFLD

　User subroutine to redefine field variables at a material point.

　）来连续改变该变量。在定义材料时同样只需要指定一个参数随该变量变化的表，ABAQUS在计算

　时自动插值计算材料参数。

　对你的情况，不需要使用该子程序。

　If you like to use USDFLD, in the subroutine, you can get step number and just use a if

　 to determine the field variable value due to the step number. Say, step 1 is 1 and step

　 2 field variable become 2. In the material definition, just define material properties

　 as a function of your field variables. Here I use elastic as example. I＇m not sure about

　 the grammar since I do not have manual at hand. You need to check them. something like

　 this:

　 *elastic, dependencies

　 2e5,0.3,,1

　 1e5,0.3,,2

　 Here I assume your material does not have any temperature dependency and you need to make

　 sure field variable is the 4th value. Then, in step 1, you have 2e5 as Young＇s modulus

　 and in step you have 1e5.

　 20.输出计算过程中的总质量和总刚度矩阵

　element matrix output,elset,file name=user defined,MASS=YES,OUTPUT FILE=USER DEFINED

　21.先张预应力：

　先张预应力：

　INITIAL CONDITIONS, TYPE=STRESS, REBAR

　22.用户子程序的使用

　假设你的输入文件为

　a.inp

　b.for

　那么在ABAQUS command 中的命令应该是这样的：

　 abaqus job=a user=b

　对于abaqus64pr11，command 中输入：

　 abq64pr11 job=a user=b

　就可以了。

　当然首先你要用 cd 命令进入输入文件所在的当前文件目录。

　强烈建议使用command来操作。

　子程序文件名后缀应为 .for，而不是 .f

　23.怎样设定用双CPU机器进行ABAQUS计算

　try ABAQUS j=your-job CPUS=2

　你是什么操作系统系统。我只试过2000和xp。

　在上述系统下abaqus仅支持domain并行计算。CMD为

　 abaqus job=? parallel=domains domains=2 cpus=2

　直接用CPUS=2不行的说。不过6.41我没试过。

　如果直接用abaqus job=?

　此时看任务管理器，应该2个cpu都是50%在跑：explicite。我没试过stantard下会如何。

　如果在xp下，单击进程名explicite,然后右键。可找到选项分配cpu。

　可以指定进程在那一个下跑或同时在两个上跑。还可以指定优先级。

　如果同时算两道题，可选一个进城为低于标准。这样还能干别的事。

　如果只是上网、下载就用默认就行了。

　我是2 X xp 1800+的机器。

　 24.中途停止正在运算的JOB

　If you run Abaqus with interactive mode, you can use Ctrl+C to kill it.

　If you run in backgroud, you can use:

　abaqus suspend job=jobname to suspend your job;

　abaqus resume job=jobname to resume your suspended job;

　abaqus terminate job=jobname to terminate your job.

　25.自适应网格技术

　一个模板：

　 *HEADING

　...

　*ELEMENT, TYPE=..., ELSET=ACOUSTIC

　Data lines to define acoustic elements

　*ELEMENT, TYPE=..., ELSET=SOLID

　Data lines to define structural elements

　*SURFACE, NAME=TIE_ACOUSTIC

　Data lines to define the acoustic surface interface with the structural mesh

　*SURFACE, NAME=TIE_SOLID

　Data lines to define the solid surface interface with the acoustic mesh

　*TIE, NAME=COUPLING

　TIE_ACOUSTIC, TIE_SOLID

　...

　*STEP

　*STATIC

　*ADAPTIVE MESH, ELSET=ACOUSTIC, MESH SWEEP=10

　...

　*END STEP

　*STEP

　*STEADY STATE DYNAMICS, DIRECT

　...

　*END STEP

　26.abaqus计算与内存

　求解器会估算内存使用量，ABAQUS计算时估算的是物理内存加上硬盘共享内存，如果所

　需超出二者之和，就停止计算。如果一次物理内存装不下，就先放在硬盘内存区上，然后逐步读入物理内存，这就是为啥物理内存扩大后，解算精度可以提高的原因了。同时隐式不能超过3G，显式不能超过2G ，就像增量步不能超过1000000步一样。（硬软件决定的，不可能无限的）而且，显示、隐式对CPU和内存的要求不同，这要看你算的问题的类型了。对于同一个算题，隐式计算，CPU的作用没有内存大，就是说CPU的数量和快慢没有内存的改变对计算时间的影响大。而显式则不同，CPU的个数和速度对计算时间有很大的影响。

　 8Q＇Bm-M6SL,x9P*F仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentabaqus内存调整

　U)I7U^}K.uQH仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM如果你的内存比较小，可以在这里调整下

　9|!^WgA ?|Simwe.com|仿真|设计|有限元|虚拟仪器1，edit job 对话框中点击memory 然后把里面的内存值设为比你现有的内存稍大些，至于memory policy 可以选moderate 或者maximum

　x0co8?“q(G1U|Simwe.com|仿真|设计|有限元|虚拟仪器2，可以将虚拟内存调大或者直接固定虚拟内存大小，还可以在abaqus_v6.env这个文件里面修改（ABAQUS＼6.5-1＼site）

　27.质量缩放

　ABAQUS Analysis User’s Manual 11.7.1

　概要

　准静态分析或某些动态分析中，少数尺寸较小的单元控制稳态时间增量，为提高计算效率，ABAQUS/Explicit常采用质量缩放的方法。质量缩放可用于：

　 § 缩放整个模型，单个单元或单元组的质量

　§ 多步分析中，缩放每个分析步中的质量

　§ 分析步起始或整个过程中进行质量缩放

　质量比例缩放可通过以下方式执行：

　 § 采用给定的常数因子对特定单元进行质量缩放

　§ 对所有指定的单元采用相同比例因子进行质量缩放，使单元组内任意单元的最小稳态时间增量等于用户给定的时间增量

　§ 仅对单元组内稳态时间增量小于用户给定时间增量的单元进行质量缩放，使这些单元的稳态时间增量等于用户给定的时间增量

　§ 缩放所有指定单元的质量，使每个单元的稳态时间增量等于用户给定的时间增量

　§ 对于金属成形分析，基于网格的几何形状和初始条件，自动进行质量缩放

　简介

　显示动态过程常用于解决以下两类问题：瞬时动态响应计算和含复杂非线性效应（最常见的是复杂的接触条件）的准静态模拟。由于求解动态方程时采用了显示中心差分法，平衡方程中离散的质量矩阵对计算效率和精度都起到了关键性的影响。如果恰当地运用质量缩放方法，可以在保证计算精度的情况下，大大提高计算效率。然而，最适合于准静态模拟的质量缩放技术与动态分析中必须采用的质量缩放方法存在很大差异。

　准静态分析

　对于应变率无关材料的准静态分析，自然时间并不重要。为节省计算时间，有效的办法是有两种：减少分析的时间步长或人为地增加模型的质量（质量缩放）。对于率无关材料，这两种方法产生的效果相同；但如果模型中含有率相关材料，首选质量缩放方法，因为该方法保留了自然时间。

　准静态分析的质量缩放方法通常用于整个模型上执行。然而，当模型各部分的刚度和质量不同时，常选中模型的某部分进行质量缩放或对每部分分别进行缩放。任何情况下，都没有必要减小模型质量的实际值，并且随意地增加质量通常都会影响到计算精度。对于大多数准静态问题，一定程度的质量缩放可以增加ABAQUS/Explicit时间增量，从而减小计算时间。然而，必须保证质量的改变和随之增加的惯性力对计算结果没有显著影响。

　诚然，修改材料密度也可以达到质量缩放的效果，但本节描述的方法灵活性更大，特别是在多步分析中。

　有关准静态过程中质量缩放的讨论，可参看”Rolling of thick plates,” Section 1.3.6 of the ABAQUS Example Problems Manual

　动态分析

　动态分析中，自然时间度量非常重要，为了获得瞬态响应，必须精确地表示模型的实际质量和惯性。然而，许多复杂的动态模型包含了一些尺寸极小的单元，使显示动态分析采用很小的时间增量。这些小尺寸单元通常是在生成复杂网格时形成的。通过在分析步起始时对这些控制单元的质量进行缩放，可以显著地增加稳态时间增量，而对整个模型的动态行为的影响可以忽略不计。

　对真正的动态问题，只能对少数单元进行质量缩放，并且不允许明显地增加整个模型的质量，否则会降低动态结果的精度。

　有关动态分析中采用质量缩放的讨论，可参看“Impact of a copper rod,” Section 1.3.10 of the

　 ABAQUS Benchmarks Manual,

　稳态时间增量

　本节中，“单元稳态时间增量”是指单个单元的稳态时间增量；“单元－单元稳态时间增量”是指某个单元组内单元稳态时间增量的最小值；“稳态时间增量”是指整个模型的稳态时间增量。

　模型中引入质量缩放

　显示动态分析中有两种质量缩放方法：定比例质量缩放和变比例质量缩放。两种方法可以分开使用，也可以结合起来使用。质量缩放可用于整个模型，也可以用在单元组上。

　定比例质量缩放

　定比例缩放方法是对组装成全局节点质量阵的单元质量进行缩放，该方法在分析步起始时执行。缩放后的质量阵接着用于该分析步的每个增量步，除非同时采用了变比例缩放。如果接下来的分析步中没有重新定义质量缩放，定比例缩放方法将延续下去。

　定比例缩放的两种基本方法：直接定义质量缩放因子，或者用户定义最小的稳态时间增量，由显示动态分析过程来确定质量缩放因子。

　定比例缩放方法简单，在分析步起始时修改准静态模型的质量，或修改动态模型少数单元的质量，使它们不控制稳态时间增量的大小。由于只在分析步起始时执行一次质量缩放，因此该方法的计算效率很高。

　 Input File Usage:

　*FIXED MASS SCALING

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: At beginning of step

　变比例质量缩放

　在分析步中变比例缩放方法同期性地缩放单元质量。当采用此类型的质量缩放方法时，需定义最小的稳态时间增量：质量缩放比例因子自动计算，并按要求施加到单元上。

　当分析步中控制稳态时间增量的刚度变化剧烈时，变比例缩放非常有用。准静态体积成形分析和单元压缩量很大的动态分析中常会出现这样的情形。

　 Input File Usage:

　*VARIABLE MASS SCALING

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below:

　 Create: Semi-automatic mass scaling, Scale: Throughout step

　直接定义质量缩放因子

　对于动能必须保持很小的准静态分析中，直接定义质量缩放因子很有用。用户可以对指定单元组内的所有单元定义一个固定的质量缩放因子。这些单元的质量在分析步开始时被缩放将在整个分析步中保持不变，除非通过变比例质量因子进一步修改质量。

　 Input File Usage:

　*FIXED MASS SCALING, FACTOR=scale_factor

　例如，将单元组内的单元质量放大10倍

　*FIXED MASS SCALING, FACTOR=10., ELSET=elset

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: At beginning of step, Scale by factor: scale_factor

　定义理想的单元－单元稳态时间增量

　对于定比例或变比例质量缩放方法，用户可以对单元组定义单元－单元稳定时间增量。ABAQUS/Explicit将确定必要的质量缩放因子。定义单元－单元稳态时间增量时，有三种互斥的质量缩放方法可供选择。下面对各种方法进行详细阐述。

　为了确定各增量步的稳态时间增量，ABAQUS/Explicit首先以单元－单元为基础确定最小的稳态时间增量。然后，根据模型的最高频率用全局估计的算法确定稳态时间增量。选用两种估计值中较大者作为稳态时间增量。一般来说，全局估计器确定的稳态时间增量大于单元－单元估计器确定的值。当采用定比例或变比例质量缩放方法，并且对单元组指定单元－单元稳态时间增量时，直接影响到单元－单元稳态时间增量的估计值。如果模型中所有单元采用单一的质量缩放定义，则单元－单元估计值将等于单元－单元稳态时间增量给定值，除非采用了罚方法强加接触约束。罚接触会导致单元－单元估计值比单元－单元稳态时间增量值的给定值略小。由于使用了全局估计器，实际使用的稳态时间增量值可能大于单元－单元稳态时间增量给定值。如果仅对模型的一部分执行质量缩放，没有经过质量缩放的单元，它们的单元稳态时间增量可能小于单元－单元稳态时间增量的给定值，这些单元将控制单元－单元稳态时间增量的估计值。结果，仅对部分模型进行质量缩放时，时间增量通常都不等于单元－单元稳态时间增量。

　如果显示动态分析步的定比例时间增量尺度是以初始的单元－单元稳态极限为基础或者直接指定的，将按6.6.3节描述的规则计算使用的时间增量。

　均匀缩放质量

　对于要求模型的动能保持很小的准静态分析，均匀缩放质量很有用。这种方法与直接指定比例因子相似。两种情况下，所有单元的质量都统一地根据单一比例因子进行缩放。然而，用均匀质量缩放方法，质量缩放因子由ABAQUS/Explicit确定，而不是由用户指定。对所有单元施加均匀的、相同的质量缩放因子，使得这些单元中的最小稳态时间增量等于单元－单元时间增量的给定值,dt。

　 Input File Usage:

　Use either of the following options:

　*FIXED MASS SCALING, TYPE=UNIFORM, DT=dt

　*VARIABLE MASS SCALING, TYPE=UNIFORM, DT=dt

　ABAQUS/CAE Usage:

　仅对单元稳态时间增量低于给定值的单元进行质量缩放

　仅缩放单元稳态时间增量低于给定值的单元，适用于准静态分析和动态分析。对于增加关键单元的稳态时间增量非常有用。

　分析步开始时网格包含了控制稳态时间增量的小单元时，采用定比例的方式缩放这些单元的质量，并且以理想的时间增量开始分析。只增加控制单元的质量意味着可显著提高稳态时间增量，但对整个模型的影响可以忽略不计。

　对于累积变形产生一定数量的小单元的分析过程，采用变比例的方式对这些单元进行质量缩放，从而限制稳态时间增量的下降。

　 Input File Usage:

　Use either of the following options:

　*FIXED MASS SCALING, TYPE=BELOW MIN, DT=dt

　*VARIABLE MASS SCALING, TYPE=BELOW MIN, DT=dt

　ABAQUS/CAE Usage:

　通过质量缩放使所有单元具有相等的单元稳态时间增量

　对所有单元进行质量缩放，致使它们具有相同的稳态时间增量，有效地影响到模型的特征频谱。由于会引起质量属性的剧烈变化，所以这种方法只适用于准静态分析。并且它意味着某些单元的比例缩放因子可能会小于

　1。

　 Input File Usage:

　Use either of the following options:

　*FIXED MASS SCALING, TYPE=SET EQUAL DT, DT=dt

　*VARIABLE MASS SCALING, TYPE=SET EQUAL DT, DT=dt

　ABAQUS/CAE Usage:

　全局质量缩放和局部质量缩放

　对单元组指定定比例或变比例的质量缩放，用于对模型的局部区域进行质量缩放。对于指定的单元组，重复定义质量缩放方法时，局部定义将覆盖全局定义。

　 Input File Usage:

　Use either of the following options:

　*FIXED MASS SCALING, ELSET=elset

　*VARIABLE MASS SCALING, ELSET=elset

　ABAQUS/CAE Usage:

　分析步开始时进行质量缩放

　定比例质量缩放仅用在分析步开始时指定质量缩放比例，并且通常缩放原始单元的质量。直接定义缩放因子时，用给定的缩放因子进行质量缩放。如果指定单元－单元稳态时间增量dt,基于该值进行质量缩放。如果缩放因子和单元－单元稳态时间增量都指定，首先用质量缩放因子进行缩放，再依据单元－单元稳态时间增量和所选的定比例缩放类型，决定是否再进行缩放。

　局部的质量缩放针对特定的单元组定义。如果没有指定单元组，将对模型中所有单元采用定比例质量缩放。每个单元组只允许使用一个定比例质量缩放。多个定比例缩放包含的单元组不能有重叠部分。对指定单元组定义的局部质量缩放会覆盖全局质量缩放。

　 Input File Usage:

　*FIXED MASS SCALING, FACTOR=factor, DT=dt, TYPE=type, ELSET=elset

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: At beginning of step,

　 Scale by factor: factor, Scale to target time increment of: dt

　例：

　对某个准静态分析过程，对模型的所有单元，定义质量缩放因子为50.0。而且，假定定义了该缩放因子后，少数极小的或形状极差的单元仍影响稳态时间增量，使之小于理想值。为了进一步加大稳态时间增量，用以下选项：

　 *FIXED MASS SCALING, FACTOR=50., TYPE=BELOW MIN, DT=.5E-6

　指定的缩放因子使模型中所有单元的质量放大50倍。如果质量放大50倍之后，仍有单元的稳定时间增量小于0.5*10E-6，这些单元的质量将进一步放在，使得单元的稳态时间增量等于0.5*10E-6.

　分析步中进行质量缩放

　指定单元－单元稳态时间增量的变比例质量缩放方法，可用于分析步起始和整个分析步内。同时，必须指定频率或间隔数，用来定义质量缩放执行的频繁程度。不进行质量缩放的增量步中，使用的时间增量一般都不同于单元－单元稳态时间增量的给定值。

　可以对特定的单元组定义局部质量缩放，并且每组单元只能定义一次。如果没有定义单元组，将针对所有单元进行变比例缩放。多个变比例质量缩放所涵盖的单元组不能重叠。局部质量缩放覆盖全局质量缩放。

　 Input File Usage:

　*VARIABLE MASS SCALING, DT=dt, TYPE=type, ELSET=elset

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: Throughout step, Scale to target time increment of: dt

　等增量步间隔进行质量缩放

　用户可以指定两次质量缩放之间的增量步数。例如，指定频率为5，表示执行质量缩放的时刻分别为分析步开始、第5、第10、第15个增量步。值得注意的是，选择过小的频率会会增加计算时间。

　 Input File Usage:

　*VARIABLE MASS SCALING, TYPE=type, DT=dt, FREQUENCY=n

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: Throughout step, Scale to target time increment of: dt, Scale: Every n increments

　等时间间隔进行质量缩放

　另外，可以指定执行质量缩放的时间间隔。例如，对历时1.0秒的分析步指定时间间隔为5，表示执行质量缩放的时刻分别为分析步开始、0.2秒、0.4秒……、1.0秒。

　 Input File Usage:

　*VARIABLE MASS SCALING, TYPE=type, DT=dt, NUMBER INTERVAL=n

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: Throughout step, Scale to target time increment of: dt, Scale: At n equal intervals

　分析步的开始和过程中采用不同的质量缩放

　对于有些情况，理想的做法是在分析步之初采用一种质量缩放方法，而分析过程中进行修改。

　 Input File Usage:

　采用以下两个选项：

　 *FIXED MASS SCALING, FACTOR=factor, TYPE=type, DT=dt_init

　*VARIABLE MASS SCALING, TYPE=type, DT=dt_min, FREQUENCY=n

　or NUMBER INTERVAL=n

　ABAQUS/CAE Usage:

　举例

　动态冲击分析中，网格中存在少数尺寸极小或形状极差的单元，这些单元控制了稳态时间增量。分析步之初，可对这些单元进行质量缩放。此外，冲击导致部分区域内的网格发生严重扭曲。稳态时间增量可能受冲击区内的单元控制。

　实质上，冲击区内的单元相对于刚性表面是稳态的，选择性地使用质量缩放方法可以保证整个动态响应不受影响。用指定的时间增量对这些单元进行质量缩放，可以有效地地减少计算时间。

　例如，采用定比例质量缩放，指定模型中所有单元的稳态时间增量下限值为1.0E-6。此外，采用变比例质量缩放，指定冲击区单元(elset1)的稳态时间下限值为0.5E-6。这样，分析步开始时，检查所有单元，如果单元的稳态时间增量小于1.0E-6，分别对这些单元进行质量缩放，使之达到1.0E-6。除单元组elset1之外，所有单元在随后的分析过程中保持该质量缩放。在分析步中，变比例缩放影响单元组

　elset1，使其稳态时间增量不小于0.5E-6。由于分析过程中只对单元组elset1进行了质量缩放，可能会出现整个模型的稳态时间增量小于0.5E-6的情况。

　多分析步中的质量缩放

　从一个分析步转到另一个分析步，质量缩放可以保留，也可以删除，已经缩放过的单元质量也可以重新初始化。跨分析步应用质量缩放方法，应遵循以下规则：

　 § 如果新分析步中没有重新定义变比例质量缩放方法，前一步定义的变比例质量缩放自动保留。

　 § 如果新分析步中没有重新定义定比例质量缩放方法，前一步分析结束时，无论单元的质量是否经过缩放，都将保留下来。

　多步分析中，除了第一步之外，分析步开始时由于使用了质量缩放，单元质量变化较大，可能会影响到质量的计算精度。当单元质量变化较大时，建议新分析步开始时，先用定比例缩放的方法把单元质量重新初始化（使之回到原始值），然后再定义必要的质量缩放方法。

　删除质量缩放

　在当前分析步中定义变比例缩放方法，将删除前一步中所有的变比例缩放。因此，为了将保留前一步中的变比例缩放，当前分析步中应重新对其进行定义。

　 Input File Usage:

　*VARIABLE MASS SCALING

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Semi-automatic mass scaling, Scale: Throughout step

　举例

　假定在准静态分析的第一步，单元经历变形导致稳态时间增量急剧下降。此外，假定第二步的变形对稳态时间没有明显影响。

　 *VARIABLE MASS SCALING, TYPE=BELOW MIN, DT=1.E-5, FREQUENCY=10

　*VARIABLE MASS SCALING

　第一步，定义每10个增量步对模型中所有单元进行一次质量缩放，单元－单元稳态时间增量为1 × 10–5。第二步，没有进一步采用质量缩放，沿用第一步经过缩放后的单元质量。

　重新初始化

　对于第一个分析步之外的其他分析步，默认采用定比例缩放重新初始化单元质量，使之回到原始值。这样，定比例质量缩放可用于防止缩放后的质量用于新的分析步。这种方法适用于从准静态分析步（需要进行质量缩放）转到动态分析步（无需进行质量缩放）。

　当缩放新分析步中的质量时，可以指定合适的比例因子，或者指定合适的单元－单元稳态时间增量和缩放类型。当仅对部分单元进行质量缩放时，需定义单元组。

　 Input File Usage:

　Use the following option to reinitialize the mass of the entire model to its original value for a new step:

　*FIXED MASS SCALING

　ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Reinitialize mass

　举例

　假定某个分析过程依次包括准静态分析和动态分析两步。在准静态分析步中执行质量缩放，在动态分析步中关闭该功能。

　 *FIXED MASS SCALING, TYPE=BELOW MIN, DT=1.E-5

　*FIXED MASS SCALING

　第一个分析步中，定义分析步开始时整个模型的单元－单元稳态时间增量为1.0 × 10–5.第二步中采用定比例缩放的方法将质量矩阵重新初始化。

　质量缩放适用场合

　质量缩放不会影响下列项目：

　 § 热－力耦合中的热响应结果

　§ 重力载荷，粘滞压力载荷

　§ 绝热分析

　§ 材料的状态方程

　§ 流体单元和流体连接器单元

　§ 弹簧和阻尼器单元

　与以上项目相关的密度保持不变。质量元、旋转惯量元、无限体和刚性单元也可以进行质量缩放。然而，由于没有这些单元没有稳态时间增量，所以它们只能通过两种方式进行质量缩放，一种是用户指定比例缩放因子，另一种是施加非均匀的单元－单元稳态时间增量。如果采用后者，则至少有一个单元是有稳态时间增量的。壳和梁的旋转惯性元就是基于质量缩放。

　无限元可进行质量缩放，但是，与变形体相邻的单元密度必须进行相同比例的缩放，否则边界不安定。

　金属辊轧成形分析中的自动质量缩放

　金属辊轧成形通常是准静态过程，但常采用ABAQUS/Explicit进行模拟。为了使求解过程省时，常人为增加产品质量。然而，必须优选质量缩放因子，使质量变化和对应的惯性力的变化不明显改变计算结果。比例因子选择过高会产生非准静态解。选择过低的比例因子，会导致运行时间过长。“辊轧变比例质量缩放”这一选项可以自动优先该过程的比例因子。

　自动策略的基础是半自动方法，缩放所有的单元质量使它们具有相等的稳态时间增量。该方法从辊轧过程的几个参数自动计算合适的稳态时间增量目标值。该目标值，，由以下几个参数确定：辊轧方向单元的平均长度，和送料速率V，以及产品截面上节点的个数n。送料速率是辊轧过程中稳态情况下的平均速度。分析过程中调整目标值，使它与实际的送料速率相匹配。因此，用户必须指定平均速度的估计值，辊轧方向上单元的平均长度和截面节点的个数。

　任何单元的质量都不得小于原始质量。这一点与等单元稳态时间增量的质量缩放方法有所不同。该强制条件意味着，把有惯性效果显著的辊轧问题当作准静态问题分析时，不自动调整单元质量。

　为了获得好的计算结果，建议：

　采用二维截面拉伸的方法进行网格划分

　辊轧方向上单元长度变化尽可能小

　初始的送料速率尽可能接近稳态送料速率

　横截面上的单元尺寸尽小于或等于辊轧方向的单元尺寸

　除了辊轧自动变比例质量缩放之外，没有其它的质量缩放方法

　Input File Usage:

　*VARIABLE MASS SCALING, ELSET=elset1, FREQUENCY=n,

　 TYPE=ROLLING, FEED RATE=V, EXTRUDED LENGTH=,

　CROSS SECTION NODES=n

　 ABAQUS/CAE Usage:

　Step module: Create Step: General, Dynamic, Explicit or Dynamic, Temp-disp, Explicit: Mass scaling: Use scaling definitions below: Create: Automatic mass scaling, Feed rate: V, Extruded element length: , Nodes in cross section: n

　28.ABAQUS多处理器进行并行计算的效果研究

　环境：

　 ABAQUS6.3+8 IntelXeon 700MHz CPU+4G Ram+Win2k AdvServer SP3

　在命令行模式下，abaqus命令的下面三个参数进行并行计算的控制：

　 [cpus=number-of-cpus]

　[parallel={loop | domain | supernode | tree}]

　[domains=number-of-domains]

　ABAQUS/Explicit:

　parallel参数可选domain,loop

　domain进行拓扑域并行，loop进行循环级并行（默认）

　但在NT系统下，不支持loop参数

　cpus数要可以整除domains数，也就是一个cpu可以进行多个domain的计算

　以Getting Started with ABAQUS/Explicit

　6.5

　Example: circuit board drop test

　circuit.inp分析为例(standard_memory = “256 Mb“)

　1cpu:

　abaqus job=circuit int

　cup利用率100%，运行时间506s

　2cpu:

　abaqus job=circuit parallel=domain domains=4 cpus=2 int

　每个cup利用率接近100%，运行时间402s

　4cpu:

　abaqus job=circuit parallel=domain domains=4 cpus=4 int

　每个cup利用率为80%左右，运行时间297s

　8cpu:

　abaqus job=circuit parallel=domain domains=8 cpus=8 int

　每个cup利用率为40%左右，运行时间364s

　ABAQUS/Standard:

　parallel参数可选supernode,tree

　supernode(默认)对单波前进行并行处理，tree对多波前同时进行并行处理

　domains参数无效

　对于线性方程并有稀疏刚度矩阵的模型并行计算有效

　以ABAQUS Release Notes

　2.7

　Parallel sparse solver

　t1-std.inp分析为例(standard_memory = “1000 Mb“)

　1cpu:

　abaqus j=t1-std int

　cup利用率为70%左右，运行时间390s

　4cpu (supernode):

　abaqus j=t1-std parallel=supernode cpus=4 int

　所有系统cpu均被使用，每个cpu利用率运行过程中不断大幅变化，运行时间454s

　4cpu (tree):

　abaqus j=t1-std parallel=tree cpus=4 int

　所有系统cpu均被使用，每个cpu利用率运行过程中不断大幅变化，运行时间352s

　8cpu (supernode):

　abaqus j=t1-std parallel=supernode cpus=8 int

　每个cpu利用率接近100%，运行40多分钟后还无法结束

　8cpu (tree):

　abaqus j=t1-std parallel=tree cpus=8 int

　每个cpu利用率接近90%，运行时间545s

　结论：

　多CPU并行处理对Explicit效果显著，对Standard在很多情况下效果不明显，甚至可能

　使运算更慢，只使用Standard的同学基本可以不用考虑买多cpu的机器

　29.yahoo讨论组摘录--contact+overclosure

　8776 Urgent help: overclosure of contact problems

　I have a contact problem which the mesh of the master Surface and the slave surface are of the same size. Acturally the Master nodes and the slave nodes are sharing the same coordinates. They have the same location but different node numbers. But when i run the abaqus. The warning message of “ the Overclosure of the master surface and slave surface is too much Severe.“.

　 In abaqus manual, one way to slove it is to use “ Adjust=0.0 “ in the “contact pair“ to establish contact at the Initial configuration?

　 Or to use “ contact interference“ ?

　How to manange contact when master surface and slave Surface are the same location, but will seperate or penetrate when Loading??

　Liang, I had solved a similar problem by increasing the mesh density of the Slave surface by about 2 times and increasing the distance between master and Slave surface so that there is a small gap between them. I am not sure if you Can do that in your problem but this is an option.

　These are the warnings you are seeing. What are the error messages when it aborts as seen from the message file? What exactly asre you simulating? Usually your slave surface (more deformable?) Must be more finely meshed than the master surface.

　5996 Exponential contact pressure-overclosure relationship & separation

　I am dealing with an interface where I have defined a softened interface following an exponential relationship. It appears that this configuration does not allow the separation of the two surfaces. Would you know how to disable this behavior? I have found in the doc how to not allow separation, but I have here the oposite case, I need the separation.

　You can just toggle off the box for “allow for separation“ :)

　Yes, just erase NO SEPARATION on the *CONTACT PAIR card !!

　Softened contact DO allow separation, the particularity is that you will have a contact pressure even if the surfaces are opened.

　5176frction and contact

　I have a co-axial cylindrical shell structure in which a polymeric cylinder has been placed between two steel cylinders and then pressure is applied on the external surface of the structure. I am looking for the collapse pressure of the structure when de-bonding between polymeric layer and steel layers is allowed. So far I have a model which works with low and moderate friction coefficient (0.1 to 0.6).

　 However, to make sure that my model is behaving well, I decided to run a fully bonded model (layers sharing node with each other) and then compare the result with a de-bonded model with high friction coefficient (e.g. 0.9-1.0). I expect that the buckling pressure should be close.

　 Now the problem is as friction coefficient is going up (>0.7), analysis fails without a single increment progress with messages regarding overclosure and contact opening, etc. Looking at the ABAQUS manual, it seems that there are two different methods to solve the problem. One is using NO SEPARATION contact behaviour and the other one is to use UNSYMM=NO. Both of them are working in my case when applied but I do not know which one is the correct one or the most appropriate one.

　 My questions are: 1- Is my understanding of the high friction and the way to avoid convergence problem is right? 2- How can I decide what should I choose (NO SEPARATION or UNSYMM=NO)? Any help would be greatly appreciated.

　4473 how to model the two plate contact together?

　I model the two steel plates connected, the sides of them are modeled by connector element, and the contact surface, I use contact command. But ABAUQUS, give me the warning below:

　 ***WARNING: THE SYSTEM MATRIX HAS 239 NEGATIVE EIGENVALUES.

　 ***WARNING: OVERCLOSURE OF CONTACT SURFACES ANGLEPLANE and GUSSETPLANE IS TOO SEVERE -- CUTBACK WILL RESULT. YOU MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION. SLAVE NODE 3513 OF CONTACT PAIR (ANGLEPLANE,GUSSETPLANE) HAS JUST CLOSED. A NEW CONTACT PATCH NEEDS TO BE DEFINED FOR THE SLAVE NODE.

　And I use this commond *Surface Behavior, no separation, pressure-overclosure=HARD and connenctor element is 3D Connector element : beam to located at the edge nodes of two plates I only want to model the two contact surface deform with each other without friction,They can＇t penetrate into each other＇s inside. Would you give me some suggestion! Thanks in advance

　I am not all that experienced with Abaqus, but I think what the error means is that there are 239 negative volumes / elemenst in ur model. It might be due to meshing very close surfaces together that your elements between the plates have turned inside out. Ur model has to be free of all the negative elements before you can do any analysis on it. If u are modeling two simple plates, u can try meshing them on another meshing program and then make sure that there is no negative elements in the mesh before u import it in Abaqus to do any analysis. Guys, feel free to correct me if I have mistaken.

　May be ur initial penetration is too much (modeling) or ur initial time step is too high. ABAQUS is trying to resolve overclosure there.

　Either adjust model so that it will be “just contact“ or “just before contact“ distance, or reduce ur initial time step so that u can avoid overclosure thingy.

　But it may be because of normal issue ?

　If it it plate contact, I think it will be working fine with few adjustment...

　 By the way, how are u dealing with edge line contact ???

　I get the same messages when running a contact problem. Is it realted to

　 ABAQUS CAE, did anyone tried to use the “Contact Pair“ option in ABAQUS Explicit?

　 If it worked some other way please advice, thanks

　**4400contact

　I am currently doing some analyses including contacts for a composite shell. To have a better understanding about contact, I started by modelling a composite beam consisting of three layers (Steel/Polymer/Steel). The steel plates were master surface and polymeric layer was slave. I was expecting that when I am pressing the steel plates, the polymeric layer should slide and squeeze out (the Poisson抯 effect).

　I was using solid elements (brick element). On the contact surfaces, the slave and mater nodes had the same coordinate but different node numbers. When I was using reduced integration point formulation for 1st or 2nd order elements, the solution did not converge and I had warning messages regarding numerical singularity and overclosure suggesting to change HCRIT, etc. Even when I refined the mesh for slave surface (twice as master surface), the problem did not disappear until mesh refinement in slave surface was in a way that the nodes in slave and master were not 揷oincide? Furthermore, I had to use full integration elements rather than reduced integration elements -even in 2nd order elements- to allow the analysis converges, which is strange to me. First of all, my question is why this is a case (mesh refinement in slave surface and avoiding reduced integration point formulation for the 2nd order elements)? Having succeeded in analysis of composite beam (although I have problem to understand why this is a case), I modelled a long shell (L/D=10) with three layers (Steel/Polymer/Steel), allowing the same sort of contact describing above. I ended up with memory shortage (I have modified the pre-memory and standard memory as high as I could). A simple solution to memory is to add RAM, but I am hopefully trying to solve the problem in other ways. Now my second question is, how one can 搊ptimise?The size of the model when you are dealing with contact. Any help will be greatly appreciated.

　It＇s very hard to define “optimal contact definition“ in general. It depends on many variables...but usually... 1. Define contact where only contact is occuring. Usually I run without contact (disp loading) to see where it penetrates, then grow contact from that region... 2. If it really is not serious contact (large sliding etc), try to

　 use simplified version (small sliding, etc).

　 So, I think the less contact I have it＇s better situation...

　I am not sure about full-integration, reduced integration shell...I thought ABAQUS is using kinematic relation type contact (penalty in explicit), surface is generated based on nodes/surface (master), then keep calculating relative distance between master/slave (nodes) etc...not sure where integration point thing comes in.

　 More node (refined) will generate smoother surface..depends on what contact u are using (finite sliding, small sliding, etc), the number of nodes to check distance will increase.....

　 Solution will converge if u use small step, displacement type loading until it reaches problem points, apply very small initial increment and check in what instance ABAQUS is having trouble...(assuming no serious overclosure)..

　4337 modeling an arch dam

　I＇m trying to model an arch dam by CAE scripting. As far as i＇ve understood, the solid modeling is based on extrusion operation for 3D solids. But this may not be applicable for generation of the arch dam blocks. Any idea on this?

　This may be create in others cad/cam/cae & then should be convert to *.inp for abaqus.

　Thanks for your reply. But I think ABAQUS should consider adding some more features for creating base solid models in next versions. Feature based modeling is quite a strong tool and I beleive adding further capabilities to it, CAE can be used as a very strong mesh generator.

　Do U have Abaqus 6.3.3?

　Yes, truly, all depends on terms. Up to an output of the new version it is possible to make, probably, only so with required quality.

　4141 Possible contact problem

　Hi group, I am running an impact simulation using ABAQUS/EXPLICIT between an analytical rigid body and a deformable element based solid. At the start of the step the rigid body is in contact (no overclosure) with the node based surface of the element based solid. A force is then applied to the reference node of the rigid body. I am interested in the displacement of the surface nodes on the deformable body, but they are consistently too low in comparison to the

　 experimental results I have.

　 I have read of lots of people having contact problems and was wondering if it is possible that I am having a contact related problem. I have checked the total contact force between the surfaces and it seems to be about right. Is it possible that the contact is not working properly if I am getting the right contact force between the surfaces?

　You haven＇t said anything about the nature of your deformable body. Is the real object undergoing plastic deformation or creep that you are not modeling?

　In rereading your original question , I see that you are doing impact. Creep is probably not an issue and I＇m not sure about plasticity or viscoelasticity.

　How great is the discrepancy between FEA and experiment?

　Hi Danny, thanks for your reply. Impact occurs over a short time (~1ms) so as you say creep, viscoelasticity shouldn＇t be a factor. I expect a displacement of 2.8mm based on experiment, and FEA yields a value of 1.5mm. There are other possible sources of error in the simulation that I am investigating, but I wanted to eliminate the possibility that the contact wasn＇t working properly.

　*4004 simulate micro indentation

　I am using ABAQUS 6.2/Standard to simulate micro indentation procedure. (ABAQUS/Explicit is not available) At first, I wanna to do a very simple simulation. A sphere indenter is put into aluminum block material. I set the indenter as rigid and constrain it fully. Then the Step is static and distributed load on the bottom surface of the block. Contace pair is the top surface of the block (slave) and the rigid sphere surface (master). Unfortunately, I always got the warning information: overclosure of contact surfaces is too severe. What has gone wrong? What shall I do to get rid of it and reach the right solution?

　(1) check, if the surfaces do not touch each other

　(2) use a small increment for the calculation start

　(3) check if your rigid surface definition is o.k., sometimes the surface normal points to the wrong direction, see manuals for right definition

　--- > the program computes the distance of the slave nodes from the master surface and takes into account this surface normal - so you will

　 get always a overclosure for separated surfaces if you choose the wrong direction

　I have noticed that you are modelling micro-indentation in ABAQUS. I wonder if you know of any references that can help me to get started with a micro-indentation model? I have no idea on how to do it, since all the literature that I have found says that classical plasticity does not apply but that I should use strain gradient plasticity theory, which is not available in ABAQUS. Any help will be greatly appreciated.

　Strain gradient plasticity is just a new kind of theory to consider the size effect of material. It means the material properties in small size scale (nano or micro) is different from large scale. The theory is called to be a bridge between continuum mechanics and dislocation theory. It is mainly proposed by Gao Huajian and Nix in Stanford univeristy and Huang in UIUC. They have done indentation simulation based on the theory. (check their papers). The constitutive equation should be changed in the simulation. I believe it is a quite complicated job. I am using classical plasticity theory because currently my scale is much larger. (Dozens of microns)

　**3666 CONTACT PROBLEMS

　I am facing problem defining the Contact.

　I must define whats the Geometry, Its 2D problem, Axis Symmetric. 2 Rigids body are holding Material Sheet, In the CONTACT Region of Rigid bodies and material, The material is not allowed to flow when the Deformation takes place, The rest of material part is facing Fluid Pressure, I hope you understand its Forming Process.

　 The problem is overclosure, The Nodes of material body with in Rigid bodies get the problem and Convergence is not reached so it aborts. Some body please guide me what options i should use to overcome this problem

　Are you applying load or displacement to your rigid bodies?

　*Does your problem start out with an overclosure? If it is your intention to have the rigid bodies contact the deformable solid with an initial line-to-line contact, you may want to use the ADJUST parameter on the *CONTACT PAIR command.

　I am not applying ne force on rigids bodies, they are just holding the material, the nodes of material that are under rigid bodies are fixed,

　*Good, ABAQUS has trouble resolving contact with load applied to rigid bodies; applying a displacement to the rigid bodies rather than load is much

　 easier for ABAQUS to resolve. As Danny suggested, ADJUST might help. Also, it will make your message file large, but you can try using *PRINT, CONTACT=YES to give you a better idea of what is going on at the contact interface (this will tell you which nodes are overclosing and by how much).

　This is puzzling. Do you mean that:

　a) The rigid bodies are controlled by displacements (zero or non-zero) and

　b) The nodes on the surface of the deformable body which are in contact with the rigid body are fixed????

　 Then nothing will move and no deformation will occur.

　This geometry is almost the same as an example in “ABAQUS Example Problems Manual - Volume 1“, the Example is CYLINDRICAL CUP DEEP DRAWING.

　 Now I want a little change in this.

　1-- Instead of Punch I want to use the FLUID PRESSURE.

　2-- Material held between the HOLDER and DIE should not be moved, when the deformation takes place meaning only part of material under FLUID PRESSURE is deformed. I hope this will clear what i want to do.

　I will be greatly looking for the suggestion.

　3535 composite structure and contact

　I would like to model 搒liding?Behaviour of middle layer of a three-layer composite structure. For simplicity, imaging there is a solid rectangular beam (C3D8RH), consisting of three layers. The middle layer (core) is much softer than the exposing (face) layers. Then structure is going under external pressure from one of the face layers. In reality, the nodes of core have the same coordinate than faces (i.e. Sharing the same nodes). But when you want to study the sliding behaviour of core, you need to have, I suppose, different series of node sets (for face and core respectively), but with same coordinates and then allow the contact between them. My problem is I will get the following warning message in the .msg file, which it is understandable:

　 ***WARNING: OVERCLOSURE OF CONTACT SURFACES 揂NN-I?Br> and 揊LINE?IS TOO SEVERE -- CUTBACK WILL RESULT. YOU

　MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION.

　 Now if I want to change the HCRIT, it will not be a case and in the case of 搒ame coordinates? I do not know how much it should be.

　Have anybody crossed such problems? Any help would be

　 greatly appreciated.

　***I have faced similar problems earlier and the following worked for me:

　 Use ADJUST parameter to specify the depth of adjustment zone. This might allow you to avoid changing the HCRIT value.In most contact problems, overclosure, due to smaller than required HCRIT, occurs right at the beginning of the simulation, which means that the slave nodes are penetrating the master surface in the initial configuration itself. ADJUST parameter will reposition the nodes and enable modeling of the sliding behvior with good accuracy.

　Thank you very much for your email. Actually, I have tried it and still could not perform the analysis successfully. I also had problem with memory shortage which means that I should make my model smaller.

　If it does not make any problem, could you please forward me the contact part of input file you have written for your own model to my email address

　3517 Contact Simulation...Severe Cut Back Problem

　I am having problem in my simulation my message file is giving this warning

　 ***WARNING: OVERCLOSURE OF CONTACT SURFACES SA000005 and SP000009 IS TOO SEVERE -- CUTBACK WILL RESULT. YOU MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION.

　I have reduced the value of HCRIT till 0.0001 m. But still having the same message. Is it due to the high load? And how can I get rid of this message as it sometimes stop after a no. Of iteration step due to this problem. Waiting...

　Check your surface normals. If they are in the wrong direction, ABAQUS thinks your surfaces have already penetrated each other.

　Thanks Rich for the reply...

　My *.dat file shows that the normals have been assumed by ABAQUS. The message is

　 ***WARNING: NODE 4 ON SURFACE SP000009 HAS FACETS WITH NORMAL VECTORS DIFFERING BY MORE THAN 30 DEGREES. CONVERGENCE DIFFICULTIES MAY OCCUR AT THIS NODE WITH FINITE-SLIDING CONTACT. THE NORMAL CONTACT DIRECTION AT THIS NODE WILL BE (0.0000,0.70622,0.70799).

　 Any remedial action to avoid this.

　I think this topic of Contact Penetration has been discussed in the previous mails. Try a word search with contact.

　*3504 Values for softened contact

　I am considering using the pressure-overclosure=linear option on the *SURFACE BEHAVIOUR card. I am interested in using this in conjuntion with the NO SEPARATION parameter and hopefully to have tension between the contacting bodies also (I don＇t think hard contact with no separation does this?).

　 My problem is what values to pick for the pressure-overclosure relationship. There is no particular data I need to use so is there numbers which people tend to use in general to get the advantages of softened contact (I am having problems with hourglassing) and to have tension available?

　I believe that NO SEPARATION parameter can＇t be used in conjuntion with the modified surface behavior such as pressure-overclosure=linear at least in ABAQUS/standard. My understanding to this is that once contact was established the slave nodes will be constrained fully to the master surface.

　***3200 Contact overclosure

　I have a model with 3D-solid elements (quad. Tetrahedral) and a static step in abaqus/standard. I defined a contact between two round surfaces, initially separated. During the second step the surfaces get in contact with each other. I allways have an overclosure of the surfaces and the increments are getting smaller and smaller untill the analysis stops. I changed master/slave and tried every value of hcrit, but it didn＇t work. It＇s allways the same. What could I do else?

　Perhaps you should use C3D10M instead of C3D10. C3D10M are designed for contact problems.

　Are you using ABAQUS/Standard? And is it a high enough speed impact? If so, like I had to, Explicit may be a better option. Just my thoughts,

　Thanks for your answer, but I use C3D10M elements.

　*Yes, I am using abaqus/standart. Can I do in explicit a static analysis?

　--Not really I think, it is only for dynamic analysis. Sorry, I thought from your previous mail you were running a dynamic analysis. Didn＇t read properly!

　A couple of suggestions. Try adding command *CONTACT CONTROLS, AUTOMATIC TOLERANCES in your steps. This helps to reduce ＇chattering＇ at the contact face (see the manual for detail.

　 Alternatively change your element type to 4 nodes (especially on the slave surface side) maybe C3D4 ?

　*The element choice plays a crucial part in contact element. Many of the problems like overclosure, Hourglassing atc can be reduced to a great extent by the right choice of elements. I would suggest you try to change the element type and run it once and the try observing at what point of the second time step is the overclosure starting. Keep checking for these in between the solution. Also if the geometry is not too complex try the brick elements.

　**Try to impose a small displacement instead of a load in order to initialized the contact.

　 First step: small displacement

　Second step: you replace the displacement by the load

　Becarfull: do not forget your real boundaries.

　 Second method: use the parameter APPROACH in *Contact Controls command

　Do you have SMALL SLIDING or finite sliding ? How do you move your second body? Just a force load ? If this is true, try a step with a small prescribed motion to close the contact, then run a second step, remove the boundary and add your force. This will run. Use *CONTACT CONTROLS, ANALYSIS=DISCONTINUOUS.

　1378 Explicit and Contact

　I＇m trying to do a simulation of a axisymmetric bolt thread being formed. I followed an example in the examples manual as a guide, but I＇m still having problems with the contact.

　 1) The die is made of a rigid surface.

　2) The bolt material starts out as a rectangular slab.

　3) The rigid surface is the master surface.

　4) Friction is used (mu=0.2).

　5) Material of bolt is steel

　6) The bolt mat＇l is fixed along the line of symmetry.

　7) The die moves towards the bolt at a velocity of 1600 mm/s.

　 The problem I encounter is that the forming result is totally bizarre.

　 First, I get an initial warning that my overclosure is 1.33 (mm?) And the ABAQUS will compensate for it in the first step. I think my problem is here.

　What I end up with is that the material “moves“ to contact the die without the die moving in the first increment. At least that is what is looks like to me. I checked it in ABAQUS/POST, and it does look like the elements first form around the die (which doesn＇t move in that increment), and then the when the die does move, it doesn＇t deform

　 anything.

　Any ideas or suggestions???

　What I end up with is that the material “moves“ to contact the die without the die moving

　 Abaqus post and viewer rescale the graphis to keep things centered on each increment. I＇m guessing this is causing the illusion of motion.

　Have you checked if the outward normal at the contact surface is in the right sense? Perhaps it is the reason of having an initial overclosure.

　Yep, that was the problem. The model works just fine now. Thanks for all help. Since i＇ve never used rigid analytical surfaces, I didn＇t realize this was an issue. It is.

　*1294 contact overclosure

　We have a problem involving too many contacts. Is there any other option other than ADJUST (*CONTACT PAIR - keyword) parameter to curb penetration / overclosure. We understand this ADJUST is to take care of geometry modelling / meshing errors. While we donot give this the solution terminates at the beginning itself. Giving ADJUST = 0.0 makes the solution progress. But we have been told not to use this parameter. Is there any other way in place. Thanks in advance.

　Between what type of bodies are you modeling contact? The problem depends to some extent on this. If you are modeling contact between deformable bodies, you might want to try the WEIGHT parameter. Also, try adjusting the mesh size of the elements belonging to master and slave surfaces.

　You might try *CONTACT INTERFERENCE,SHRINK

　**593 Couple of Questions

　I＇m wondering nobody couldn＇t give me guide(s) about my questions? Could you please reply me if you can?

　Soheil

　(below is my previous message)

　Dear Group

　I＇m a phd student at Heriot-Watt University, UK.

　I want to model the bending of pipe on rigid circular surface. I faced with problem in “contact“ of pipe which has been modelled by “S4“ element and the rigid surface which modelled by “R3D4“ elemnt.

　I would like to know:

　 1- What is the best way for modelling of a 3-D circular rigid surface (e.g. An exterior surface of a quarter of cylinder)?

　 2-What is the command name in ABAQUS for “merging“ of nodes? In other softwares like LUSAS or NISA, the “merge“ command is used. I

　 guess it is MPC in ABAQUS, isn＇t it?

　 3-Which parameters should be taken into account when you want to choose a magnitude for increment, i.e. INC=....?

　 4- I faced the following WARNING message in .msg file. What should I do to overcome it? What does it mean “VALUE OF HCRIT“?

　 ***WARNING: OVERCLOSURE OF CONTACT SURFACES SHELSUR and RGDSUR IS TOO SEVERE -- CUTBACK WILL RESULT. YOU MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION.

　(SHELSURF is positive side of pipe and RGDSUR is rigid surface).

　5-Finally, about print of .odb files, what is the procedure? I couldn＇t find any rough described path in ABAQUS/Viewer manual. I want to know what should I type in “Print Command“ box when you open “print“ window. I seek help of everybody and I do appreciate your help in advance.

　1.Best way of modeling rigid surface is building up the surface with analytical surface. Analytical one is always better than numerically discretized one in view of analysis costs, and accuracy. It＇s like mises is better than tresca in terms of easiness of numerical implementation.

　 2.It＇s MPC TIED. It does not merge the nodes as you want, it simply eliminates dof of specified nodes in stiffness matrix and makes this node inactive during current analysis.

　 3.Im not sure what is the keypoint of this question. If you want to use user defined time increment , use DIRECT keyword in *STATIC.

　 4. This message generally arise when there is large amount of penetration during current sdi. HCRIT means ... If there developed contact penetration larger than this value, there＇ll be cut-back rather than more iteration.

　305 Problems of Overclosure

　I am trying to model a machining process, where the workpiece is clamped onto the machining platform. Different sections of the workpiece are machined at different stages and this requires new set of clamps to be introduced into the model for each stage. The clamps are idealised by using contact analysis. This is because I do not want some sections of the workpiece to be rigidly constrained.

　 The problem is that I have severe overclosure in the model for some contact pairs when I activate new contact pairs in the restart analysis, to introduce the new set of clamps for machining the next section of the workpiece.

　 I am not very familiar with contact analysis, but I have read that we can use commands such as *contact interference, shrink or *contact pair,adjust to solve such problems (I was not aware of these command initially). However, as these comands are defined at the first step of the analysis, would the overclosure problems be rectified for the reactivated contact pairs that were previously inactive in the restart analyses? I would very much appreciate any other suggestions to solving the problem too.

　The clamps are idealised by using contact analysis. This is because I do not want some sections of the workpiece to be rigidly constrained.

　I＇d say dont use contact pairs at all. You are adding a complicated nonlinearity where all you want is a stiff-but-not-rigid boundary condition. Imagine your clamp is jsut a stiff spring and use *EQUATION.

　 If you must model the clamping, start with the clamp slightly away from the part and push it into contact as a step. Soft contact helps quite alot with convergance as well.

　*5994 assembly parts instance options in ABAQUS

　I hope everybody is doing fine. I am using the assembly, parts instance options in ABAQUS/Explicit to build my model. I am facing problems at the common interface of two parts of a model. Does anybody know--how abaqus joins the two parts or how they define the commmon nodes for two parts? When I am checking my data file I found that globally abaqus is definig two nodes at the same location. And if it is doing that then how the nodes are coneected between parts? Do I need to define the connection at the common interface or assembly option is supposed to do it automatically?

　 If any body has any idea or any experience with the assembly parts instance options in ABAQUS please, reply, ASAP.

　Abaqus will not join the two instances together. You will need to define some interaction between the two instances using interaction module. If u want to merge two parts then you can merge using assembly module itself. If u want to merge the nodes of the two instances then you can do the following, if your part is a geometrical representation then, mesh it and create input file for that. Open it using abaqus, which will open it with only nodes and elements (orphan mesh)...then you can opt for merging the two orphan mesh part instances in assembly module (but make sure to you dont merge the nodes)...then that creates third part with the

　 nodes at the interface unmerged. Then in part module you can click the node on the interface you want to merge and then merge it,

　Thanks a lot for your reply. But actually i am not using ABAQUS/CAE. We don＇t have the license for that. So I am building my model trough text editior. I am definig the parts first and meshing them. Each part contains tie constraints is some regions. The tie constarints are defined with in each part. Then I am definig the assembly by positioning the parts using *parts instance option. When the whole assembly is built I am defing the boundary conditions and contact regions. I have contact between steel and concrete at one face. And i am definig that (in the step data) afetr definig the assembly.

　 My model is running and there is no error message in any file. But the only thing is that each part is behaving individually.

　 As u have have mentioned about the merging option, how can I do that by text editor?

　The solution is to define a *CONTACT,TIED card between the two parts. But the stresses in this region won＇t be accurate...

　584 Couple of Questions about “CONTACT“

　I＇m a phd student at Heriot-Watt University, UK.

　 I would like to know:

　 1- What is the best way for modelling of a 3-D circular rigid surface (e.g. An exterior surface of a quarter of cylinder)?

　 2-What is the command name in ABAQUS for “merging“ of nodes? In other softwares like LUSAS or NISA, the “merge“ command is used. I guess it is MPC in ABAQUS, isn＇t it?

　 3-Which parameters should be taken into account when you want to choose a magnitude for increment, i.e. INC=....?

　 4- I faced the following WARNING message in .msg file. What should I do to overcome it? What does it mean “VALUE OF HCRIT“?

　 ***WARNING: OVERCLOSURE OF CONTACT SURFACES SHELSUR and RGDSUR IS TOO SEVERE -- CUTBACK WILL RESULT. YOU MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE *CONTACT PAIR OPTION.

　(SHELSURF is positive side of pipe and RGDSUR is rigid surface).

　 5-Finally, about print of .odb files, what is the procedure? I couldn＇t find any rough described path in ABAQUS/Viewer manual.

　 I seek help of everybody and I do appreciate your help in advance.

　30.原创：无限元建立方法，希望得到加分

　loveestboy1981

　这是我探索出来的关于无限元单元建立的过程，其实无限元建立最难的不是INP文件怎么通过*element去改变单元类型，而是如何把单元节点转换为与无限元单元节点一致的情况，也就是一个单元的最后几个节点要在同一侧（外侧）。www.simwe.com f

　以下是我建立的方法，如果能得到各位的同意的话，希望能得到加分肯定。

　 &S v/Z@8B7zA/yXr:C|Simwe.com|仿真|设计|有限元|虚拟仪器以一个简单例子作为探讨对象，建立模型，划分网格，通过mesh--creat mesh part建立网格单元PART：Part-1-mesh-1,然后定义装配。这些都是基础的步骤。

　 2006-8-19 22:14

　2.jpg (21.9 KB)

　 2006-8-19 22:14

　3.jpg (60.19 KB)

　然后可以通过TOOLS-----Query看到你所选单元的NODE排列号码（首先在VIEW--PART DISPLAY OPTIONS--MESH下选择show node labels）可以看到，所选单元node号为5407, 5408, 5434, 5433, 5121, 5122, 5148, 5147.node 5148, 5147,5434,5433这几个单元靠后的节点是在靠外的一侧。但是你需要在INP文件里把NODE靠后的单元号放在后面。

　 ~+V9x(HDD!o

　R[_/m＇o|Simwe.com|仿真|设计|有限元|虚拟仪器

　]N1ox3N.ewww.simwe.com[ 本帖最后由 loveestboy1981 于 2006-8-19 22:52 编辑 ]

　附件

　2006-8-19 22:23

　4.JPG (22.27 KB)

　 2006-8-19 22:23

　5.JPG (48.33 KB)

　当不是这种情况，你想把单元node号反过来，即把nodalconectivit改为 5121, 5122, 5148, 5147,5407, 5408, 5434, 5433.则可以通过修改inp文件。首先建立JOB点write input,然后用ULtraEdit打开建立的XXX.INP。如图。可以看到我想建立无限元单元集合为单元4251到4500.则利用ULtraEdit把单元4251到4500所有的节点编号复制下来新建一个文件，然后把NODE号改变，再粘贴覆盖原来这一部分，则已经把NODE号改变完成。可以看图7和图8的区别

　附件

　2006-8-19 22:36

　6.JPG (41.11 KB)

　 2006-8-19 22:36

　7.JPG (68.33 KB)

　 2006-8-19 22:36

　8.JPG (62.13 KB)

　最后别忘记了，对我想要建立的无限元部分进行单元类型修改。如下图。然后再通过以上查看节点的方法可以看到同样一个单元的NODE号把靠后大的号调了一个面。到此改变已经成功，单元类型也已经改成CIN3D8（三维无限单元）。谢谢，有不对的地方希望大家指正。

　 9DG0v*W

　k-{+vb|Simwe.com|仿真|设计|有限元|虚拟仪器

　aS }7＼%~/@p$C*Y“Rwww.simwe.com[ 本帖最后由 loveestboy1981 于 2006-8-19 22:57 编辑 ]

　附件

　2006-8-19 22:47

　9.JPG (56.61 KB)

　 2006-8-19 22:47

　10.JPG (90.55 KB)

　 2006-8-19 22:47

　11.JPG (76.76 KB)

　 31.[分享]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 ?

　 PS1<PS2

　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 输

　出。

　。

　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 ?

　 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

　 ......

　 32.[转帖]abaqus6.4导入外来模型的几点小经验！

　abaqus6.4的CAE建模功能仍不是太强大，许多情况下，还得借助于第二、第三方软件完成模型的建立，但在导入过程中，总会或多或少的存在问题，近来做了几个这方面的试验，有几点发现，解释如下：

　1.对于三维模型，pro/e等CAD软件可以建出很好的模型，存成iges、sat、step等格式导入abaqus可以直接用，个人感觉能用sat或step格式出现的模型最好不用iges文件导入！

　 2.复杂模型，当导入abaqus时，模型可能已经枝离破碎了，无法通过几何修补将模型改好，而且分网也相当困难，但通过第三方前处理软件hypermesh，就可以很好的实现模型在导入时出错的问题，hypermesh可导入的几何模型格式相当多，一般的CAD软件的格式都可以很好的导入，并且对小毛病可以通过几何清理完成！然后加上abaqus的模板，在hypermesh分出高质量的网后导出inp直接为abaqus所用，目前是相当好的一个强强联合的CAE思路！

　 3.有时也可以以这种方式导入模型：ansys现在在国内的受众群很大，我想许多人也会用ansys，起码最基本的建模都能做到，可以发现，ansys的建模功能相比abaqus要强的多，尤其现在有了workbench后，这点就更为突出了！所以不习惯用abaqus建模的朋友而又暂时不会用专业CAD软件的朋友可以直接在ansys里建好你所想要的模型，然后再导入abaqus里用它强大的计算功能吧！

　 a.曲线救国：ansys里建的模型通过hypermesh导入分网后导进abaqus！（说细情况可以到hypermesh版找到）

　 b.对于2d草图，abaqus也可以接收iges格式，所以在ansys做好2d模型后，通过file---export导出iges文件，然后abaqus导入sketch也就可以用了！（不过，在ansys建这个模型时，不要出现area及其以上的entity，因为sketch只认线和点，所以你只能在ansys通过line、keypoint完成你要建模型的轮廓，导入应该是没什么问题！有时导入效果比autocad导入效果都好！

　 4.说说2d模型的导入问题，这点版里有多人问到了，由于CAE无法做出诸如椭圆、正弦、余弦曲线等复杂平面图形来，所以通常通过autocad作出图形，存为dxf格式然后导入abaqus，有时会出现导入sketch时什么也没有或只有几条张而已，总之导入不成功！其实，abaqus导入sketch后的图形仅有线和点两种形式，而且最重要的一点就是abaqus的导入dxf时也认为cad的版本，目前autocad版本已到2004，这时如果你存dxf格式时选的是2004版本承认的dxf格式，当你导入abaqus时，有很大可能不会成功，但如果此时你选择较低版本的dxf格式存储时，问题就会解决！（存为低版本的dxf格式是可以在保存时选择的），我做过试验，存为目前来说最低的autocad R12版本和R13版本的dxf格式时再导入abaqus里，几乎可以很好地导入，不信大家试试！不过，要说明的一点，限于abaqus导入dxf格式功能的限制，它对于常见的几种几何图形如直接、圆、矩形等可以以线的形式出现（这在abaqus的documentation里有介绍），而对于复杂曲线，则是以多个点的形式出现（拟合成那个曲线），所以复杂一点的曲线导入后会是N个点！当然，有了点，曲线还愁出现不了嘛？

　两种方法，用直接一个点一点的连起来，但这样会产生很多点且曲线的拟合度不太好，用spline方法即多样线按点的顺序连成的曲线会与你所想要的曲线相差无几！（推荐用spline方法，因为这样仅产生一条线，对于以后的分网有好处）试试吧！（当然，如果点相当多，如果你没有耐心，我就没办法了！）

　还有就是版友iloveni所说的在autocad下作图时最好不用便捷图形直接画，能用直接直线画的图形就不要用诸如矩形、正六边形等快捷作图方式作图，这样导入后会出现直线，否则只能出现体现图形形状的定位顶点，你要做的工作就只能是自己在abaqus里用线将点连起来了！

　几点体会，都试过，反正在我这都通过，应该没什么问题！大家可以补充！！

　33.abaqus的多图层绘图

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

　 S$TcwZB!b“^:}L7HSimWe仿真论坛1.首先进入可视化模块，当然要先打开你的模型数据文件（。odb）

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

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

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

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

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

　34.子结构

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

　@6S u9J5O/T3d

　2.为什么要用子结构

　$oo:YS%M4N7K-N仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM 使用子结构并不是为了好玩，凡是建过大型有限元模型的兄弟们都可能碰到过计算一个问题要花几个小时，一两天甚至由于单元太多无法求解的情况，子结构正是针对这类问题的一种解决方法，所以子结构肯定是对一个大型的有限元模型的，譬如在求解非线性问题的时候，因为对于一个非线性问题，系统往往经过多次迭代，每次这个系统的刚度矩阵都会被重新计算，而一般来说一个大型问题往往有很大一部分的变形是很小的，把这部分作为一个子结构，其刚度矩阵仅要计算一次，大大节约了计算时间。

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

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

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

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

　M k[D“W;K|Simwe.com|仿真|设计|有限元|虚拟仪器子结构是通过你指定的节点与其他的单元建立联系的；在abaqus的6.4版本中只有子结构这个概念没有超单元了，其区别就是子结构可以求得单元内部准确的解而超单元不行；当你定义子结构的时候不要包含太多的单元，因为单元的刚度矩阵集成的时候会花掉太多的时间，可以用更多的含有较少单元的子结构代替。

　9Q“hE|yG仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent5.在abaqus中子结构的用法

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

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

　3dYe6wk9Uy Hr:CSimWe仿真论坛我所了解的大概有三种不同的方法：

　 J^lL1w/]e:y4A“A仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1.最强大的当然是采用umat的方式，不过需要有深厚的有限元基础，一般人不推荐使用

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

　%S|sA)h＼T*wRFwww.simwe.com 3.采用abaqus的import命令将前面分析的结果传递到新的分析之中

　＼]!q8t c2j*X*W

　 dSimWe仿真论坛这里介绍下第二种方法

　%R@j/K/?[I

　ki3m3USimWe仿真论坛*什么是场变量|Simwe.co

　＼

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

　8S4q!mq7[6Z)pwww.simwe.com*怎样使用场变量

　r0N6KhVj

　s|Simwe.com|仿真|设计|有限元|虚拟仪器其实场变量用的较多的实在热力学和流体力学的分析种，这里介绍的仅仅是在固体力学中的用法

　 1.定义场变量

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

　initial conditions，type＝field，variable＝n（场变量的编号）

　f&X:y^s

　Inp{仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent Set－1（你定义的结点集），1.0（场变量的值）

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

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

　*field, variable=1u-V^ e

　Set-1,1

　2G?$M8g:P FA4x|Simwe.com|仿真|设计|有限元|虚拟仪器当然也可以同时使用initial和field，当你指定的场变量改变时，默认材料的参数是在增量步间线性变化的。

　 z-N+i&]6B5gl仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent

　 8N＇Ex

　iCwww.simwe.com2.建立材料参数和场变量之间的联系9tHs&V3|(T

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

　T c ^/{8Ng7wtP

　杨是模量泊松比 field1

　-U3Cq(mI200.E9, 0.3, 1.|S

　180.E9, 0.3, 2.

　.RRTEaSimWe仿真论坛如果你用的是命令格式，则在inp文件里面键入：

　 *ELASTIC, DEPENDENCIES=1 www.simwe.com-TMZ4a3eH`

　200.E9, 0.3, , 1

　180.E9, 0.3, , 2

　7K2J+~n+U3v(a)Wq仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent第四个参数表示场变量的值＇}

　3.注意

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

　(G*?4A-xVo:{仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent*STEP，name=step1

　*STATIC

　B0Z&Ci6V/js Y#{SimWe仿真论坛*FIELD, VARIABLE=1

　NALL, 1.

　6＼ S}rB仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM *END STEP

　*STEP，name=step2

　1yiz3}X

　{SimWe仿真论坛*STATIC

　k+q*KD+kA仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent *FIELD, VARIABLE=1SimW

　NALL, 2.

　*END STEP

　`.^ W/zPQ$hSimWe仿真论坛cae步支持场变量，所以你必须自己更改inp文件，

　P:R6@HU*||Simwe.com|仿真|设计|有限元|虚拟仪器4.技巧

　＇Xm

　c.jb_+A9E!O仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM 如果材料的参数变化比较复杂，一般是利用副职曲线来定义场变量值的变化，

　*FIELD, VARIABLE=1，amplitude＝？？？

　+?x7^

　@K8I$w仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent 对幅值曲线步清楚的自己可以看手册仿真

　36.模型的重启动分析－restart

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

　bwww.simwe.com1.什么是restart

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

　o lmM)Ii:IaLlH5k ＼仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent*restart,就是将一个复杂的模型分析过程分成很多的阶段，甚至是一个increatment step一个阶段，你可以对每个阶段的结果进行检验，然后进入下一个阶段进行分析。

　 5e:d~&dz9H*HSimWe仿真论坛2.重启动需要那些文件

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

　l(W/J^|+B仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent仿

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

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

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

　1{g4Z＇w!Qf_仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentcae默认加入了重启选项，不过可以在step－>output->restart request里面设置输出的频率，也就是frequency。

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

　x#D%Kfh#@bYH仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent 4.如何重启

　TnQm,MSimWe仿真论坛你要指定一个重启点，

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

　 8udt o_ ?6t`B8T仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM

　 p#K#Q;|-{uSimWe仿真论坛 5.注意

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

　 37. ABAQUS的单位心得

　Quantity

　SI (mm)

　US Unit (ft)

　US Unit (inch)

　(Lxcad)cm

　Length

　Force

　lbf

　Mass

　tonne (103 kg)吨

　slug

　lbf s2/in

　(102 kg)

　Time

　Stress

　Pa (N/m2)

　MPa (N/mm2)

　lbf/ft2

　psi (lbf/in2)

　N/cm2

　Energy

　mJ (10–3 J)

　ft lbf

　in lbf

　Density

　kg/m3

　tonne/mm3(1012)

　slug/ft3

　lbf s2/in4

　102 kg/cm3(108)

　Ace

　m/s2

　mm/s2

　cm/s2

　单位:建议采用国际单位制

　采用m、kg、N、s国际单位制时，重力加速度9.8m/s2,质量为kg

　密度为7850

　kg/m3

　MP,EX,1,210e9

　MP,DENS,1,7850

　ACEL,9.8

　采用mm、Ton 、N、s时，重力加速度9800mm/s2, 密度为7850e-12

　Ton/mm3，应力MPa

　MPDATA,EX,1,,206e3

　 MPDATA,PRXY,1,,0.3

　 MPDATA,DENS,1,,7.85e-9

　ACEL,0,0,9800,

　力=质量*加速度

　关系为：N=kg.m/s2=0.001t.1000mm/s2=t.mm/s2

　采用cm、 Mcm、N、s时，重力加速度980cm/s2,力为N

　密度为7850e-8

　Mcm/cm3

　MPDATA,EX,1,,206e5

　 MPDATA,PRXY,1,,0.3

　 MPDATA,DENS,1,,7850e-8

　ACEL,0,0,980,

　力=质量*加速度

　关系为：N=kg.m/s2=0.01Mcm.100cm/s2=Mcm.mm/s2

　其中Mcm=100kg

Tags：

上一篇：2018年度节能宣传周及低碳日活动情况总结
下一篇：高考作文例文及写作指导六篇

首页 > 心得体会 > 工作体会 / 正文

abaqus常用技巧总结

猜你喜欢