Abaqus高效使用手册:常用小技巧集锦

1.   故障诊断基本方法

问:

提交分析后,应该

1)  察看dat文件,出现的warning一般可以不去理会,但如果出现error, 则分析无法进行;

2)  如果分析能够进行,则察看msg文件(ABAQUS/Standard分析)或sel文件(ABAQUS/Explicit分析)。如果msg文件出现error, 则分析同样无法进行;msg文件中的大部分warning都可以不去理会,但以下warning意味着模型有问题:

● Numerical Singularity

● Zero Pivot

● Too Many Attempts Made for This Increment:分析无法收敛,而不收敛的原因则可能是多种多样的,常见原因是接触或非线性材料定义不当,或网格太粗,需要具体问题具体分析

答:

1)根据上述error或warning信息找出问题。

2)简化模型,去掉模型中的复杂参数(接触、塑性材料等),做一个简单模型,成功了再做复杂的。

3)约束刚体位移 / NUMERICAL SINGULARITY

2.   Abaqus显示

问:

ABAQUS/CAE初始化草图后,无法显示栅格

答:

主要是显卡问题

在 显示——属性——设置——高级——疑难问题中的硬件加速调到最小或一格

3.   磁盘空间不足

问:

在进行接触分析遇到了一个问题,错误提示为:ERROR: UNABLE TO COMPLETE FILE WRITE. CHECK THAT SUFFICIENT DISK SPACE IS

AVAILABLE. FILE IN USE AT FAILURE IS subcasejob-1_metsp.uft.

可我计算机C盘剩余8G,D盘45G,abaqus安装在D盘,内存为8G 单元数量为26万体单元,可为什么还说磁盘空间不够呢。

答:

分析过程中除了在当前工作目录下生成odb, msg, dat等文件外,还会在默认的临时文件目录下生成很大的临时文件(此目录好像是c:/temp),分析中止后,这些临时文件会被自动删掉,所以尽管你看到c:有很大空间,但在分析过程中可能会空间不够。你的d盘很大,可以把临时文件的目录修改到d:, 方法是:

在ABAQUS安装目录下的文件夹site中可以找到环境文件abaqus_v6.env,在其开头添加scratch='D:\scratch'保存此文件,在D盘创建目录scratch. 重新启动ABAQUS/CAE.另外,在abaqus_v6.env中还可以看到关于内存的设置:

#对INP文件进行预处理所允许使用的最大内存

pre_memory = "256 mb"

#

# 使用ABAQUS/Standard求解时所允许使用的最大内存

standard_memory = "256 mb"

可以把这两个值改大些(但不要超过你计算机内存的实际大小).

6.5 Linux系统

问:

abaqus在linux系统上要比winxp节省很多时间…………不知有没有哪位也试验过?

答:

linux要比winxp的分析速度快15%, 且linux能够利用的内存上限要比winxp大。

磁盘阵列没用过。如果分析时内存的使用量超出了物理内存的大小,买一条内存肯定是有用的。

4.   不要在塑性材料上施加点载荷

问:

在我的模拟一个桶土接触问题中,我想加一个和坐标轴成一定角度的集中力,但是直接加集中力会导致不收敛,因此我改用加位移荷载的方法,求出反力,再反推合力,以及合力与坐标轴的夹角,请问这个合力方向是不是与所加的合位移方向相同,本人认为是相同,不知道大家有什其他的看法,请教一个?

答:

有些情况下一样,有些情况下不一样,你试试就知道一样不一样了。施加位移荷载时的反力就等于力载荷。

“直接集中力会导致不收敛”:

一是加载之前要有一个只施加很小载荷的分析步,使各个接触关系都建立起来。

二是如果受集中力的材料是塑性,会因为局部变形过大而不收敛,应该将力施加在一个小区域上,让多个节点共同受力。

5.   屈服准则

问:

在塑性加工中,不同的材料用如下变形抗力模型表示。请问:这个公式指的是整个弹塑性过程还是忽略弹性的?如果包括弹性,那么弹性模量又与该公式有何关系?屈服极限应力在该式中有体现吗?该公式反应的是真实应力应变关系吗?它的适用范围是什么?


答:

塑性加工中的具体材料模型有成千上万个,大类可以分为理想刚塑性(忽略弹性、忽略硬化,主要应用于弹性变形相对于塑性变形较小,可以忽略的情况下,如:锻造等)、理想弹塑性(忽略硬化)、弹塑性(在材料成形时其弹性变形不可忽略,例如板料成形)、等等。

你的这个模型属于弹塑性模型,没有忽略弹性。

屈服极限是材料发生塑性变形时的临界应力,是材料的本身性质,在本公式中不能体现。

纵向应力与纵向应变的比例常数就是材料的弹性模量,也叫杨氏模量。也不能在此公式中反映出来。

任何材料模型都不可能完全反映真实的应力应变关系,当然次模型也不例外,只是最大近似。

一般只有根据材料提出材料模型,其适应范围也只能适用于特定材料以及特定场合。

6.   后处理 / 判断材料是否屈服

问:

DOCUMENT上说PEEQ是等效塑性应变(Equivalent plastic strain),PEMAG是塑性应变值(Plastic strain magintude),感觉前者主要是判断材料是否屈服,而后者呢??

答:

在比例加载时(即加载过程中主应力方向和比值不变),大多数材料的PEMAG和PEEQ相等。这两个量的区别在于,PEMAG描述的是变形过程中某一时刻的塑性应变,与加载历史无关,而PEEQ是整个变形过程中塑性应变的累积结果。例如,单向拉伸一个圆柱体,使其发生塑性应变,再通过单向压缩使其恢复初始长度,则最终的PEMAG为0,而PEEQ是拉伸和压缩过程中塑性应变之和。


8.1 接触分析不收敛的常见现象和解决方法

做接触分析时,应观察msg文件中的以下信息

(接触状态的变化:1个节点由不接触变为接触,10个节点由接触变为不接触)

CONTACT CHANGE SUMMARY: 1 CLOSURES 10 OPENINGS.

SEVERE DISCONTINUITY ITERATION 2 ENDS

CONTACT CHANGE SUMMARY: 0 CLOSURES 4 OPENINGS.

如果分析能够收敛,CLOSURES和OPENINGS的数目会逐渐减少,否则一般是接触定义不当,可以考虑以下解决方法:

● 检查所定义的接触面、接触参数和边界条件是否正确;

● 在静力分析中,必须定义足够的约束条件,以消除各个平移和转动自由度上的刚体位移;

● 避免过约束(overconstraint);

● 合理地定义接触面、接触参数和过盈配合;

● 使用足够细化的网格;

● 不要在接触面上使用C3D20、C3D20R和C3D10等单元;

● 在接触对上设置微小的过盈量,以保证在分析的一开始就已经建立起接触关系;

● 施加临时边界条件,以保证在接触关系建立之前,模型也不会出现刚体位移;

● 将分析过程分解为多个分析步来完成,让各个载荷分别在不同的分析步中逐步施加到模型上,避免使接触状态发生剧烈的改变。


8.2 板料成形的接触问题

问:

请问厚度变化怎么看

在板料冲压过程中,冲压过后,我想看看板料的厚度变化情况,怎么看?我在输出变量设置里面看到有thickness,可是设置后,好像那个不是指板料的厚度

谢谢大侠指教

答:

体单元没有法子直接输出厚度。其实体单元已经很直观了。

可以查看厚度方向上下两个节点间的DISTANCE当然可能对不正,要采用一定的策略配准。厚度随时间变化的曲线某一点处的还是可以输出的。很复杂。。。。

板料成形还是壳单元比较划算些,一来可省大量的计算成本:如果要得到比较精确的结果,实体单元要划得比较细才行。如果用减缩单元,厚度方向要至少划四层才行,而别的方向上单元的尺寸要相对差不多。就是差不多的一块板子,划实体单元要多用十几倍的单元,再加上最小单元长度变小了好多,又是十几倍的时间,如果用EXPLICIT,相同的问题,要多耗掉100多倍的时间,10多倍的磁盘空间和内存,是很不现实的。二来用实体单元好象模拟板材起皱有困难。

7.   改变接触状态

问:

杀死后加part这个命令我已经掌握,现在还有一个比较严重的问题不能解决:如何在杀死单元的位置加上另一个part,同时这个part还要与周围的物体有接触作用。怎么办??

答:

一开始就在同一个位置定义重合的part 1 和 part 2, 并定义part 2 的接触。

开始时杀死part 2 的单元和接触。杀死part 1 后,激活part 2 的单元和接触。

8.   时间/增量步的含义

问:

Abaqus的History Output中的时间是真实的时间吗?还是时间增量步的概念?

答:

如果是和时间无关的静态分析就不是真实事件,如果是和时间相关的蠕变分析、动力学分析等就是真实时间

9.   用户子程序和FORTRAN

问:

我是第一次调用subroutine.调试的*.inp和*.f也是从abaqus中的help找来的.

在调试过程中,已通过了fortran的编译,生成了*.obj文件,我已将以下的文件dp.exe,aba_para_sp,aba_para_dp,并将aba_para_sp改为aba_para,但调试仍出现以下问题(如图),请大家看看是什么问题?谢谢了

答:

在本版搜索“valueerror”,可以找到很多。

我运行没问题。ABAQUS要求的FORTRAN版本是Compaq Visual Fortran 6.0。

不用自己在Fortran里编译*.f。直接把帮助文件里的 uhetvalx.inp 和 uhetvalx.f 保存在默认的工作目录下,把 uhetvalx.f 改名为uhetvalx.for. 然后输入:

abaqus job=uhetvalx user=uhetvalx interactive

答:

我总算把这个问题解决了,是fortran软件忍得祸。因为我在另一个机器也装了abaqus,和compaq visual fortan 6.5,还是出现同样的问题。我就开始想,是不是软件的问题,找了另一个形式的(里的内容一样,只是serial不同),安装完后,verify 成功,真高兴呀!

问:

INP文件和.for文件一起执行加入umat...

偶用document里的文件来运行:abaqus job=rotcyl2d user=rotcyl2dvumat interactive结果报错 roblem during compilation-df.exe not found in PATH.

答:

首先要保证安装了正确的fortran版本,见ABAQUS Installation and Licensing Guide A.1 Software。 Windows下的ABAQUS 6.5 要求是Compaq Visual Fortran 6.0。

错误信息Problem during compilation- df.exe not found in PATH, 意思是找不到fortran的df.exe , 可以设定一下windows的path.

10. 单元选择

问:

我对同一有限元模型在相同的冲击载荷下,分别用这两种软件进行了计算,都选用显式算法,可结果相差比较大,请问各位热心人,是什么使他们的计算结果有如此的差异呢(网格划分一样)?

答:

估计你的ABAQUS模型用的是reduced integration单元,比如C3D8R,这种单元只能察看单元积分点上的应力,节点应力则很不准确。可以自己做一个悬臂梁线性弯曲实例试试,和解析解或者C3D20R,C3D8I的结果比较一下。详见 Getting Started with ABAQUS 4.1.2 Reduced integration。LS-DYNA 的单元不知是不是这样。

问:

请教一下各位,我在Z向加载局部的均布荷载,采用一般静态分析,而计算出来的Z向位移比理论计算值小得多,在8204节点计算为0.17971mm,而理论计算应略为0.4449mm,相差太大,从应力分布看,没什么问题。请问是什么原因?(文件比较大,分成了两部分)

答:

检查一下你的边界条件是否和理论计算完全相同(比如理论解是无限大区域,而你的模型还不够大)。在应力集中部位应该再细化网格。C3D8R单元只有一个积分点,精度较差(尤其是应力结果)。用C3D8I, 得到8204节点U3=0.18. 模型中没有塑性材料和hard接触属性,可以用C3D20R。位移很小,没必要用NLGEOM, 就可以加快速度。对称问题,可以只取1/4建模。大部分区域应力很小,网格可以划粗些。

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空