Abaqus使用经验记录.doc

Abaqus使用经验记录
1 Verification一闪而过，在你的 PATH 里面加上 %SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem 后，就不会 一闪而过了！！就是不能进入verify的方法是查找verify.log文件上面会有信息
在verify中可以找到需要的软件，fortran和visual等版本
C:\ABAQUS\Commands
2如何分析圆柱型模型的径向应力和环向应力
各位：2 l" z) Q1 b1 h- |
我建立了一个圆柱模型，要分析一下它的径向应力和环向应力，但ABAQUS里面只有MISS应力，是不是要建立一个局部坐标系，如何建？各位帮忙指导一下！谢谢！
建立柱面坐标系：! n$ T( D/ k3 b, S" C
$ z. V. ~1 f8 y6 y3 m2 a9 ^/ L
Tools－－－Datum－－－－Create DAtum－－Cylindrical( {" {/ T: N/ Z( i% F;
3怎么在command里来修改内存呢（abaqus）
***ERROR: THE VALUE OF 256 MB THAT HAS BEEN SPECIFIED FOR STANDARD_MEMORY IS $ `1 w# t' N) S1 v9 r: i
TOO SMALL TO RUN THE ANALYSIS AND MUST BE INCREASED. THE MINIMUM
POSSIBLE VALUE FOR STANDARD_MEMORY IS 375 MB.
: l5 Q W8 D& ]+ U' `
怎么在command里来修改内存呢？
for instance: }% ]# E1 G: Q( H, H6 B
abaqus job=Job-1 memory=512
4求助abaqus内存使用极限
abaqus计算时需要2g以上内存，配有4g的内存，尽管不出现内存不够的error，但还是不算，log里说需要设定3gb开关，设了也不行，请各位指教

64位不会出现这种问题，现在我用的就是64位系统，ABAQUS也有为64位系统的安装盘。

5画网时容易死机，内存使用过高，怎么办
我的模型不大，但画网的时候老容易死机，内存使用近100%，怎么回事啊
你多少网格啊，我在hypermsh中划几百万是没有问题的，abaqus没怎么试过，也有可能是版本的问题。
abaqus 新增的bottom-up的网格划分功能很强大，在向HM靠拢了
6、 3G开关的启动
间的分配。无论系统的物理内存量是多少，Windows 都使用 4 GB 的虚拟地址空间，其中 2 GB 分配给用户模式进程（例如应用程序），另外 2 GB 分配给内核模式进程（例如操作系统和内核模式驱动程序）。在物理内存等于或大于 1 GB 的系统中，可以使用上述两个启动开关，为应用程序分配较多内存 (3 GB)，为操作系统分配较少内存 (1 GB)。

当你的内存超过2G时，xp系统认为已达上限，但实际上系统可以承载最多3G的内存，如果你需要内存超过2G这时就需打开3G开关。如果要求超过3G，那估计就得换64位系统了
这里的内存包括RAM(物理内存)和VM（虚拟内存），有的朋友可能内存只有2G但是弹出要求扩展通道的错误，应该是你的虚拟内存太大引起的，可以把它调到3G一下就没事了
所以在32位下，不管是物理内存还是虚拟内存最好不要超过3G,多了应该也没用
打开3G开关的方法：
我的电脑-属性-高级-启动和故障恢复系统-启动-编辑，进入boot
在multi(0)disk(0)rdisk(0)partition(1)\WINDOWS="Microsoft Windows XP Professional" /noexecute=optin /fastdetect /usepmtimer

添加multi(0)disk(0)rdisk(0)partition(1)\WINDOWS="Microsoft Windows XP Professional 3GB"/3GB /noexecute=optin /fastdetect /usepmtimer
保存即可
每个电脑的上面那段东西可能不一样，没关系只要在上句的基础上添加红色部分再黏贴就行了，要紧挨着上一句
重启电脑后，选择带3G的那个系统进入即可
具体可以搜索3G开关，上面有具体的说明
我这里提供一个英文版的说明，也是从论坛上下的，谢谢上传的朋友
7、 查找单元数目及发现计算中的问题需要近dat文件

8\后处理网格变坏，怎么调回来
做的是单桩静载试验模拟，后处理时网格变坏了，如图，怎么调回来呢
不影响计算结果但太难看了也不好，希望指点一下
你是用View Cut来看截面应力，网格没变坏，是正常的。可以有不同解决方法，常用的方法是去掉网格线，Option-> Common Option, 在弹出来的菜单，选择'Free Edge'或‘No Edge'。
9\接触的次面选择情况
一律用node-surface，大面选择surface，小面用node去掉两边
下部tie，的从面去掉上面一圈，用node
10\保证通过的地应力平衡方法（用excel加逗号）
本人是个初学者，在论坛上看到地应力平衡法的贴子，发现有很多地方没讲清楚，现将地应力平衡法进行完善，调试通过，与大家分享。SC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
. ^# t6 |" U6 H: g  j
地应力平衡方法：: M1 W! @8 W& Z/ }
第一步：建立模型，材料，分析步（GEOSTATIC）5 C( l4 [$ s* Y. X! ~% ]  Z2 B5 `
第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8SimWe5 ~3 x6 ~) I* I9 g
第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)  (请严格按照这个格式，注意大小写的字母)仿9 W7 [, z8 J! O; e4 Q! d
第四步：提交J0B，完成后
1 x( t4 _: h7 v" _$ ?% H( @第五步：按以下步骤,[Module]Visualization-Report---Report Field Output---下拉菜单里面选择0 I4 T0 M' `6 D& w' u
centriod，然后依次把s11、s22、s33、s12、s13、s23点选上，setup页面选择报告文件的名字---***.txt---Write中选择Field Output-------------ok!!!$ J3 L. ^9 w/ i7 ?# `+ ~
第六步：# }: c# L. l8 [' l" k$ y
先打开excel，在excel中打开***.txt，选择“分隔符号”—选择“导入的起始行”---点下一步-----在出来的界面上点“完成”------保存为CSV逗号分隔文件***.csv（文件内容格式如下）. m2 c$ D8 ~; z% b  e7 g- c; R
单元号     S11    S22   S33   S12   S13   S23 （请注意，在保存内容中没有这一行的）0 o$ }8 }  W+ T, \: C9 Y
   1        ,     .    ,   .    ,   .   ,    .   ,   .    ,
& g' j1 w5 ?( s/ |; d4 e% X, H   2        ,     .    ,   .    ,   .   ,    .   ,   .    ,  
# g/ {  W! \' P* T! l9 S2 z3 X   .         ,     .    ,   .    ,   .   ,    .   ,   .    ,
2 k! Q/ b! b7 R  a2 [! s) R   .         ,     .    ,   .    ,   .   ,    .   ,   .    ,  
9 ]. {/ q, A0 t3 O' I这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些这些数字这些数字。在CSV文件中，数字本身是用逗号隔开的，不需要另加逗号。
6 M" J9 \7 J$ X  T! E/ v7 W% P; q1 @
! Z8 k/ J& z, X+ R* i( r
2 s3 C: t4 h+ P( S! s/ C. ~第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上： *initial conditions,type=stress,input=***.txt
完成仿% ~8 H: S! N' ~sd

4 h: q' c! Y/ c
% O% D* I' p4 y0 ?! h4 W& n第八步：重新提交JOB，OK1 w; i% f3 A) A( Y: F" a; u
11材料属性变更
Re:急：请问能不能对同一个part定义两种材料属性啊？
以前在板上看到过大侠发过，具体的找不到了，方法应该是这样的，
在cae中定义材料参数的时候，* Y- I% _( x. D6 R' x8 U% S! G# Q
property模块：定义材料性质时将Number of field variables设为1。
Field 1可以设为1，2，每个数字对应不同材料数据。 # M7 ]8 U6 g3 a# k* E/ ^
设置几何体为一个set，如Set-1。 6 `" d( Q* }3 M9 Q6 H
step模块的定义多个分析步 & t/ V+ \) g2 w- \* W
在step中插入3 B6 Z$ V5 M9 V: R. q. j) ?# b
*field, variable=1
Set-1,1
则step1中Set-1的材料属性为Field 1为1时的数据
在第二个*step之中添加
*field, variable=1 % ]/ D q0 B" `9 ?7 b1 L1 N5 Y, p
Set-1,2
则step2中Set-1的材料属性为Field 1为2时的数据
这样就可以了，当然在开始的时候要把物体附上材料常数的哦。
12\一定要从弹性模型看起
13、如何输出实体截面的弯矩啊，急！！
在INP文件中进行编辑;在*Output, history, variable=PRESELECT后面添加以下语句0 j X! T% {' f0 h& J5 G
*Section Print, name=a1, surface=asd, axes=local

SOM, SOF
注意中间要空开两行，其中name名字自己定义，surface的名字是自己定义的那个面
14、输出面的弯矩方法
在history添加
*Section Print, name=a1, surface=a, axes=local
som,sof
15. 水压力施加
请问：输入静水压力时，zero pressure height和reference pressure height 是什么意思啊？该怎么输入？谢谢了
比如说静液压力，p=a×g×h（a是密度，g重力加速度，h是深度），那么zero pressure height就是规定了静液压力为0的z坐标值，reference pressure height 设定了具有你比所给定的压力值位置的z坐标值。' [; r" u& Y6 ]* ]6 n7 O# W0 b
. R5 d1 a+ N5 I. Y; {. O
还有必须要注意的是，ABAQUS里面静液压力都是随z坐标升高或降低的。
比如你要在模型上施加相等的静水压力2MPa，那么你可以( v# J7 ~0 @/ j' v- b3 b( T3 I
把整个模型放在靠近z=0的位置，设置
Distribution = Hydrostatic$ @ f9 K0 |! p5 Z
Magnitude = 2
zero pressure height = 1E+10) { |4 j9 v& P3 Q4 e4 k
reference pressure height = 0

那么这个load会在各处施加近似2MPa（因为近水压力的0点在很远很远的地方）的垂直表面的压力。
16 求助！为什么孔隙水压力没有变化？

你用visualization module 中的tools→create field outputs→from fields，将后续分析步中的孔压减去初始分析步中的孔压即可得超孔压分布。

17移动legend
viewport－－－－－viewport annotation option

选Gneral－－legend 可以隐藏
选Legent－－Upper left corner 可以移动

18模型中接触对在空间中有位移,如何输出接触对面中从面和主面的相对滑动位移?
在Step模块中设置
Field output中择contact中的FSLIP，这个就是要求输出的接触面的相对滑移的% n0 j1 z2 d6 H0 A
在后处理中选择相应的接触对，看结果就可以了

19. 相对输出
响应减去输入即可。, N' j! m _$ F" |! W* e
cae后处理可以通过tools----->XY Date提供的相减算法实现，很简单。
后处理中能输出相对位移云图吗，就是只输出在某一级荷载下产生的位移云图，可以吗？
你可以试试 tools/create filed output，可以对结果进行运算，应该可以满足你的要求
20. 平衡地应力的问题
# Z( x2 Z) g
我两种方法我都试了，我觉得还是下面的命令简单，当然每个人不一样，习惯了都行：* R* X/ I5 H1 D; N) q4 V3 [
*el print) N: O& V+ u& \% E3 n* s
s,
) }  `2 \# @0 L: A这个命令加在以下位置，非常方便

了版上很多帖子,受益菲浅!0 V2 `# g+ `6 f/ t1 w, c, r5 T
关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!1 f1 X9 C! q( u' g8 U  V5 D
1.先施加重力荷载的作用,可以在cae中实现;5 }* v6 J3 @. }. p  t
2.在inp文件中的output request中写上
* ~* @3 t9 }, X  a; o& `1 y$ L+ D*el print- M* i: o* J' f, F# H8 C; ]4 O9 K
s,- M4 r+ S5 o8 m6 q$ U
这样就会将施加重力荷载后的应力输出到*.dat文件中了;  S$ ^* T6 v& F6 {
. @5 v3 w: T% a5 M! M: z6 l
3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如: n% ~9 L! \+ ^, M4 D7 K3 X0 {  {5 T
ELEMENT  PT FOOT-       S11         S22         S33         S12     & t1 d2 O% l. E/ @/ f1 I' I
                NOTE
6 ^# [( U% L/ B& O% v# B- R2 I( e8 Y, A0 l
           1   1        -1.2598E+05  -1434.     -3.1852E+04   892.7) H$ }; N3 |- X7 O8 D
           2   1        -1.2249E+05  -6287.     -3.2194E+04   1223.! D2 c0 e2 ]/ _9 o- z4 y
           3   1        -1.1795E+05  -497.7     -2.9611E+04   1664.2 g! z5 b$ w) y5 t
           4   1        -1.1210E+05  -7240.     -2.9834E+04   1992.
+ e! Y/ v# n- U. Y0 G+ V* g           5   1        -1.0485E+05   579.0     -2.6068E+04   2600.
( H% r* V" p2 W2 V7 f# R           6   1        -9.5803E+04  -8272.     -2.6019E+04   3031.
% O  l+ M/ z( m$ }+ r           7   1        -8.4709E+04   1915.     -2.0698E+04   4083.
# O4 s" O/ h, x; n0 T4 u& @$ ?           8   1        -7.0634E+04  -9746.     -2.0095E+04   4339.) O0 n1 W" |+ m( |' O) O- j! E
           9   1        -5.1088E+04   5401.     -1.1422E+04   8519.
' I7 u$ j6 L# d7 a% o$ J3 L          10   1        -2.4353E+04 -1.1150E+04  -8876.      1.2126E+04
( ^' Y6 H! E' E9 U4 B          11   1        -1.2847E+05   268.1     -3.2050E+04   738.1) z: f( W$ `- }- c; g. g
          12   1        -1.2786E+05  -9868.     -3.4433E+04   629.1
$ D5 p  E+ e  N8 g' y; Q8 k0 r7 ~, e          13   1        -1.2938E+05  -4224.     -3.3402E+04   502.5! R, R6 n/ H' Y6 M2 I& M, e; o
          14   1        -1.3039E+05  -3458.     -3.3461E+04   165.92 k3 i! T- @1 L
单独存为一个*.dat文件,
. f$ |/ X7 j7 u8 _- [  g+ s4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为:  instance名称.序号  ;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!" `- _. J9 ?/ J- s- z
5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),
  B/ r+ S7 s' O& c6.最后在inp文件的step之前写上' h/ m! |5 n) x7 |( a1 a9 g- K
*initial conditions,type=stress,input=文件名.csv即可!
* @3 U: p: U: c* C- @这种方法不需要用python,比较简单,希望能对大家有用!9 u& _' m4 e* j4 t' J$ q

边坡的地应力平衡

地应力平衡贴

19输出桩轴力和弯矩

具体的sof是什么结果我也没搞明白，还望哪位弄弄看
1先定义体内的界面
用1, 2
*Elset, elset=SS2, instance=PART-1-1, generate
49, 51, 1
*surface,type=cutting surface,name=s1
0,1,0,0,1,0
ss1
*surface,type=cutting surface,name=s2
6.4,6,1,0,1,0
ss2
*End Assembly
**
** MATERIALS
**
*Material, name=MATERIAL-1
*Elastic
3e+11, 0.167
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1; W& E3 d0 H. o8 T: n1 q
在每部最后然后写输出语句
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*section print,name=f1,surface=s1,axes=local,frequency=1,update=yes
sof,som（单元集）
*section print,name=f2,surface=s2,axes=local,frequency=1,update=yes
sof,som
*End Step

20关于接触收敛的经验总结

5 ?" L) R6 g# ^  ?4 i  u! T8 I; K; J. R9 ?% z6 W, c% _# D# r
by hitcondor
4 T; ]) O1 ^0 s$ c: y  h
- ?& x  I% Y$ y& t" S计算了一个有多处接触的变形问题，经过若干次的尝试终于成功，总结一点经验请高手斧正：
. S% S% x( i5 s8 Y: ~, j5 n; m( \1 P: j5 g$ @* t
1           首先要判断你的问题是否是大变形问题，如果是，而你的Nlgeom选项又没有打开，那么你将为此付出沉重的代价。& m/ N- |& H# m
2            载荷的施加要分多步进行，初始载荷要很小，可以尝试万分之一的量级。
( t) Z  m7 c7 T2 }+ W* m- {* a) h3           接触建立起来之前，你可能需要增加很多额外的约束，总体的原则是要让接触的区域在初始载荷作用下平稳，不能产生刚体位移（注意：刚体位移远远不仅仅是指整个零件的大幅度移动，范围要宽的多，甚至是节点的稍微大一点点的位移都是）。( ~4 R" ]! c/ e& j: o
4          额外约束的去除可能会给你的模型带来较大的约束变动，这使得收敛变得困难，你可以尝试分多步去掉额外约束，绝不要吝惜子步的数量。子步的初始增量步要小，调试阶段你也许可以尝试0.0001，调试后发现收敛容易再放大不迟。
1 H" t/ w5 y8 O5          有时会发生这样的事情：当你去掉某个额外约束后发现，初始载荷的作用没有到达远离初始载荷位置的接触的地方，这可能使得那个约束并没有建立起来，所以你也许可以尝试在这个状态下再增加初始载荷的大小，让他的作用传播到更远一些的接触的地方。
$ U/ R0 a% I/ J7 x$ ]6          如果你的初始载荷是力，那么它对于接触的建立可能是一些不利因素，当你尝试过很多次之后若还不收敛，你也许可以把它换成位移，小位移，比如过压0.001mm，当然你需要防止刚体位移，比如增加一个反向的力。/ r% [' ~3 }+ [
7            往往我们非常关注整个固定零件系统的约束，却忽略了运动部件的约束，你是否忽略了这个问题呢？
5 ^; |; c1 X& v' h7 Z4 l8          你若判断能用tie连接的地方，就不要选择接触；能用small sliding 的地方就不要用finite sliding；在没有什么办法的时候，你还可以定义一下torlerance of adjustment zone，0.02mm是你的一个选择；或者尝试一下surface to surface 与nodes to surface的区别，这个选项我还没有形成什么经验。
- i/ N3 a1 c4 U! N# g. b, W9            Interaction --- contact control 下面可以定义一个接触控制，如果选择了stabilization将会有助于收敛，不过建议这个选择在确实没有办法时候使用，他的理论太复杂了，超出了我的理解能力，请高手解释一下吧。  f" b  T& X) _- O! b
10     据说版本的不同，会导致收敛的不同，随着abaqus的发展，可能会有这种变化吧。2 o- H/ H- D  s2 r, F; Q4 g* F
11     祝所有的问题都能收敛。
21、过盈

22、刚体

a离散刚体怎么赋予质量
是在property --- special ----- inertia  --point mass/inertia  里面设置  
B
刚体约束是将某个区域上的单元定义为刚体，因此需要mesh；如果直接在part里面建刚体部件，如果是属性为解析刚体则不需要mesh，如果属性为离散刚体则需要mesh。不过解析刚体和离散刚体part定义是都要设置RP。
! F# X4 p8 U8 s3 v- r+ \( a其中解析刚体可以通过设置属性变为变形体或离散刚体，而变形体则只可以变成离散刚体。
7 X8 {/ D- H7 m& W1 t) o3 U  A不过刚体不会发生变形，所以无需定义材料和截面属性。