ABAQUS计算常见问题系列二

2017-03-29 by:CAE仿真在线来源:互联网

模型不能算或不收敛,都需要去monitor,msg文件查看原因,如何分析这些信息呢?这个需要具体问题具体分析,但是也存在一些共性。这里只是尝试做一个一般性的大概的总结。如果你看见此贴就认为你的warning以为迎刃而解了,那恐怕令你失望了。不收敛的问题千奇万状,往往需要头疼医脚。接触、单元类型、边界条件、网格质量以及它们的组合能产生许多千奇百怪的警告信息。企图凭一个警告信息就知道问题所在,那就只有神仙有这个本事了。一个warning出现十次能有一回参考这个汇总而得到解决了,我们就颇为欣慰了。

类似于:
Fixed time is too large
Too many attamps have been made
THE SOLUTION APPEARS TO BE DIVERGING.
CONVERGENCE ISJUDGED UNLIKELY.

Time increment required is less than the minimum specified

这样的信息几乎是无用信息(除了告诉你的模型分析失败以外,没有告诉你任何有用的东西)。宜再查找别的信息来考察。根据经验,改小增量步也不一定能收敛,虽然也有人报告过改好的先例,我是从来没有遇到过,也从来没有那个奢望。所以我一般从模型的设置入手。
必须说明的是:Error和warning的性质是完全不同的。Error意味着运算失败,but出现warning可能还能算,而且有些运算必定会出现warning(比如接触分析必定出“负特征值”,下有详述)。很多警告只是通知性质的,或者只是说明一下而已,不一定都是模型有问题。比如以下warning完全可以忽略:

xxxxx will (not) printed,这种只是通知你一声,某些玩意儿不输出了。

还有:

The parameter frequency cannot be used with the parameter field. It will be ignored(都说某某被ignored了).

A系列

如果模型能算,且结果合理,那么大部分警告信息可以不管。但是以下除外:

1numericalsigularity(数值奇异):刚体位移(欠约束)

solver problem. numericalsigularity when processing node105 instance

pile D.O.F. 1ratio=1.735e13

2Zero pivot(零主元):过约束或者欠约束。

这2个问题一般都意味着模型约束存在问题。1)、2)都会伴随着产生大量负特征值。解决方案当然第一步是检查约束了。

B系列

有一些直接导致计算aborted,那就得仔细分析了,比如:

1xxxxx is not a valid in ABAQUS/Standard(告诉你这种计算standard不支持了,换别的)
2missing property
在perperty步检查材料属性是不是都加上了。如果有梁单元,看看梁法向定义对了没有。
3Detected lock file Job-1.lck. Please confirm that no other applications are attempting to write to the output database associated with
this job before removing the lock file and resubmitting.
删除.lck文件就可以了,它是一个自动生成的文件。你也可以另存为(另取名),再运算。
4The rigid part xx is missing a refernce point
刚体(or刚体约束)都必须通过stools--reference point给它定义一个参考点 (RP),载荷都加在这个RP上。

5The area of 54 elements is zero, small, or negative. Check coordinates or node numbering, or modify the mesh seed.The elements 8have been identified in element set ErrElemAreaSmallNegZero

这个一般是节点编号不对的问题。必须是逆时针方向。

6The value of 256 MB that has been specified for standard_memory is too small to run the analysis and must be increased. Theminimum possible value for standard_memory is 470 MB

7HM to ABA的问题:集合和面的几何的名称最好不要用特殊符号和数值 (特别是从hm转过来),全部用英文字母是最安全的。

8令很多人抓狂的error code 5
1)使用了子程序, 子程序有问题, (例如数组定义跟实际赋值不一致,这个我也遇到过)
2 )模型有问题, 通常模型很大,很复杂(这个我没遇到过)
Please make sure that the mesh density of the slave surface in the tie
pair( mbly__pickedset37_cns_,assembly__pickedsurf36) is finer than the master surface.The analysis may run slower, may yield inaccurate results, and may require more memory if this is not the case
3 )硬盘没空间了(这个其实不会引起error code 5,但是出错是肯定的了), 或者是内存太小.或者产生的文件太大.
4 )关闭杀毒软件试试(特别是卡巴)

5 ) 有人认为边界条件不正确,也会引起这个错误.

9system error code 29539
关闭杀毒软件and try。

10**ERROR: Issue cannot be deleted Not all data Released
在windows中,单击“控制面板”--“系统”--“高级”-- “性能设置”- -“数据执行保护”命令 ,把pre.exe和standard.exe添加进去。重起动后尝试

11Surfaces associated with analytical rigid part MANDREL may have their orientation flipped
刚体相连的接触面方向定义反了,在接触定义的地方edit--flip

12CONTACT PAIR (ASSEMBLY_BLANKBOT,ASSEMBLY_TIE- 1_DIEDURF)NODE BLANK-1.5 ISOVERCLOSED
BY 0.0512228 WHICH IS TOO SEVERE

这往往是因为接触面的法线方向定义反了。定义刚体和shell的surface 时, 要注意选择外侧

13123456 elements are distorted。Excessive distortion of element number 5 of instance PART-1-1
如果有子程序,一般不是材料设置有问题,就是边界条件的问题

14XML parsing failure for job 1.Shutting down socket and terminating all further messages.Please check the .log, .dat, .sta, or .msg files for information about the status of the job.

http://forum.simwe.com/viewthrea ... 26amp;typeid=68

15The number of history output requests in this ABAQUS analysis (>5000)

may cause SIGNIFICANT performance problems during analysis and postprocessing
输出项太多,恐硬件资源不够。要是你确保硬件够,这条也不怕了。一般的,应该减少History中的输出项,尽量输出你最感兴趣的内容。

16Value for parameter nset will be truncated to 80 characters
nset名字取太长了,80字符限制

17compilation - ifort.exe 问题
Problem during compilation - ifort.exe not found in PATH.
安装的时候没有装好或是二次开发版本冲突。检查环境变量的设置;然后

verify一下,看看是子程序功能否能通过?

C系列

如上所说,有很多warning并一定意味着你的模型存在问题。常被问起的有:

1负特征值问题
THE SYSTEM MATRIX HAS 8 NEGATIVE EIGENVALUES.
负特征值是非线性分析的必然产物。所以不必大惊小怪,甚至久而久之,对于你熟悉的问题,你都会视而不见了。若出了问题,可先检查下有没有伴随的numericalsigularity(数值奇异)和Zero pivot(零主元)产生。如果没有,可以参考这几个方面:
1).刚体位移
2).单元异常,过度变形、过度扭曲等
3).应力应变关系有负斜率
4) 如果有流体的话,在容器发生形变的话,也可能出现negative eigenvalue 的情况,不过不会出现警告,这是被允许的
5) 失稳发生2The ratio of deformation speed to wave speed exceeds 1.0000 这个警告是指单元形变速度V(单元最大形变率/特征尺寸)和膨胀波速C(通过材料本构关系求得)的比例超过1。
解决这个问题的方案有以下几种:
(1)检查单位是否封闭(参数设置有数量级的错误),此错误新手常犯;
(2)检查网格质量 ;
(3)检查加载速度,如果条件允许的话就降低速度,该方法也很有效,但在很多情况下无法降低速度;
(4)调整STEP中的TIME SCALING FACTOR;调整STEP中的
MASS SCALING FACTOR;
(5)加*SECTION CONTROLS,NAME=SC,DISTORTION CONTROL,LENGTH RATION=0.1
或者YES也可以,加在MATERIAL 前面;或加* DIAGNOSTICS,
DEFORMATION SPEED CHECK=OFF;

或者加*DIAGNOSTICS, CUTOFF RATIO=RATIO(具体数值),在其他方法修改后还有问题的的情况下使用增加.关键字的方法见http://forum.simwe.com/thread-862510-1-1.html

3zero force/ZERO MOMENT问题

THERE IS ZERO MOMENT EVERYWHERE IN THE MODEL BASED ON THE DEFAUL CRITERION. PLEASE CHECK THE VALUE OF THE AVERAGE MOMENT DURING THECURRENT ITERATION TO VERIFY THAT THE MOMENT IS SMALL ENOUGH TO BE
TREATED AS ZERO. IF NOT, PLEASE USE THE SOLUTION CONTROLS TO RESET
THE CRITERION FOR ZERO MOMENT.
这个警告是告诉你模型中没有弯矩,没问题的,可以继续计算。
如果提示中出现特征值奇异的时候才是计算有可能出现不收敛的问题。4Degree of freedom 4 is not active in this model and can not be restrained
有限元软件计算对于实体步考虑转动自由度,所以你在边界条件中限制了 456的自由度后,软件会忽略的啊.

5The option *boundary,type=displacement has been used; check status file between steps for warnings on any jumps
prescribed across the steps in displacement values of translational dof. For rotational dof make sure that there are
no such jumps.All jumps in displacements across steps are ignored.
你采用了位移边界条件,但在平动自由度上,可能在不同的分析步骤里面有突变(你可以从sta文件里面查看),
并且应保证转动自由度无突变。通知性质的warning,一般是因为你采用位移加载方式,都出这个。

6The strain increment has exceeded fifty times the strain to cause first yield at 377 points

检查下约束够不够,约束够了就不用管了,这只是通知你,你的模型塑性应变很大,一般没多大问题。

7123 nodes are used more than once as a slave node in *TIE keyword.One of the *TIE constraints at each of these

nodes have been identified in node set WarNodeOverconTieSlave
定义接触的时候,公共节点重复定义了好几次,这样可能会出现过约束问题(只是可能影响)..

8There are 2 unconnected regions in the model.
可能是接触面由空隙,最好在接触属性中定义一个容差范围。一般各个 parts之间定义接触,aba都会这样通知用户的,只要接触设置对了, 一般没事。

9Boundary conditions are specified on inactive dof of 124 nodes.
The nodes have been identified in node set WarnNodeBCIactiveDof
边界条件定义的有问题:在124个节点的非自由度上有边界加载

10The plasticity/creep/connector friction algorithm did not converg

一般是塑性应变太大,单元扭曲导致的。可以先改为弹性模型看看是否收敛;

11The ratio of the maximum incremental adjustment to the average characteristiclength is 1.82846e-02 at node 10868 instance jiti1 on the surface pair (assembly_jq22,assembly_q22).

可以通过调大预设值消除该提示and检查网格质量。

12ELEMENT 42 INSTANCE SOIL3-1 IS DISTORTING SO MUCH THAT IT TURNS
应改进单元质量

13650 nodes are either missing intersection with their respective master surface or outside the adjust zone.
改改tie里的tolarance试试

14Dependent part instances cannot be edited or assigned mesh attributes
模型树--assembly-打击part 右键--make independent。也可以到模型树part步展开点mesh。

15The aspet ratio for nnn elements exceeds 100 to 1.
单元划分网格长宽比不合适。如果这些单元在不重要的区域(对结果肯定有些影响,

影响大小取决于这三个单元的位置,在模型中的作用等),而且能计算, 那就没问题了

16123elements are distorted

存在单元扭曲,如果这些单元在不重要的区域(对结果肯定有些影响,
影响大小取决于这三个单元的位置,在模型中的作用等),而且能计算, 那就没问题了(同15)

17***WARNING: DEGREE OF FREEDOM 1 IS NOT ACTIVE ON NODE 6 - THIS BOUNDARY
CONDITION IS IGNORED
约束了单元没有得自由度对求解没有影响,可以查看下

18热分析时出现了这样的警告“THERE IS ZERO HEAT FLUX EVERYWHERE
There is zero HEAT FLUX everywhere in the model based on the default criterion. please check the value of the
averageHEAT FLUX during the current iteration to verify that the HEAT FLUX is small enough to be treated as zero.
if not, please use the solution controls to reset the criterion for zero HEAT FLUX.
试试:
(1)是不是热源定义的问题,错误信息是说热源量几乎为零。

(2)定义热源的子程序调用命令流应该为*HEAT GENERATION,在材料模块中定义,子程序为HETVAL。

以上内容转自博客陌上良人

http://blog.sina.com.cn/u/1699672322

开放分享：优质有限元技术文章,助你自学成才