Abaqus使用经验记录.doc

Abaqus使用经验记录1 Verification一闪而过，在你的 PATH 里面加上 %SystemRoot%system32;%SystemRoot%;%SystemRoot%System32Wbem 后，就不会 一闪而过了！就是不能进入verify的方法是查找verify.log文件上面会有信息在verify中可以找到需要的软件，fortran和visual等版本C:ABAQUSCommands2如何分析圆柱型模型的径向应力和环向应力各位：2 l z) Q1 b1 h- |我建立了一个圆柱模型，要分析一下它的径向应力和环向应力，但ABAQUS里面只有MISS应力，是不是要建立一个局部坐标系，如何建？各位帮忙指导一下！谢谢！建立柱面坐标系：! n$ T( D/ k3 b, S C$ z. V. 1 f8 y6 y3 m2 a9 / LToolsDatumCreate DAtumCylindrical( / 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 Babaqus job=Job-1 memory=5124求助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，的从面去掉上面一圈，用node10保证通过的地应力平衡方法（用excel加逗号）本人是个初学者，在论坛上看到地应力平衡法的贴子，发现有很多地方没讲清楚，现将地应力平衡法进行完善，调试通过，与大家分享。SC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent. # t6 | U6 H: gj地应力平衡方法：: 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( 第五步：按以下步骤,ModuleVisualization-Report-Report Field Output-下拉菜单里面选择0 I4 T0 M 6 D& w ucentriod，然后依次把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% be7 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 Ra2 ! s) R . , . , . , . , . , . ,9 . / q, A0 t3 O I这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2.这些这些数字这些数字。在CSV文件中，数字本身是用逗号隔开的，不需要另加逗号。6 M J9 7 J$ XT! E/ v7 W% P; q1 ! Z8 k/ J& z, X+ R* i( r2 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 sd4 h: q c! Y/ c% O% D* I p4 y0 ?! h4 W& n第八步：重新提交JOB，OK1 w; i% f3 A) A( Y: F a; u11材料属性变更Re:急：请问能不能对同一个part定义两种材料属性啊？以前在板上看到过大侠发过，具体的找不到了，方法应该是这样的，在cae中定义材料参数的时候，* Y- I% _( x. D6 R x8 U% S! G# Qproperty模块：定义材料性质时将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 Hstep模块的定义多个分析步 & t/ V+ ) g2 w- * W在step中插入3 B6 Z$ V5 M9 V: R. q. j) ?# b*field, variable=1Set-1,1则step1中Set-1的材料属性为Field 1为1时的数据在第二个*step之中添加*field, variable=1 % / D q0 B 9 ?7 b1 L1 N5 Y, pSet-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=localSOM, SOF注意中间要空开两行，其中name名字自己定义，surface的名字是自己定义的那个面14、输出面的弯矩方法在history添加*Section Print, name=a1, surface=a, axes=localsom,sof15. 水压力施加请问：输入静水压力时，zero pressure height和reference pressure height 是什么意思啊？该怎么输入？谢谢了比如说静液压力，p=agh（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 ZMagnitude = 2zero pressure height = 1E+10) |4 j9 v& P3 Q4 e4 kreference pressure height = 0那么这个load会在各处施加近似2MPa（因为近水压力的0点在很远很远的地方）的垂直表面的压力。16 求助！为什么孔隙水压力没有变化？你用visualization module 中的toolscreate field outputsfrom fields，将后续分析步中的孔压减去初始分析步中的孔压即可得超孔压分布。17移动legendviewportviewport annotation option选Gnerallegend 可以隐藏选LegentUpper left corner 可以移动18模型中接触对在空间中有位移,如何输出接触对面中从面和主面的相对滑动位移?在Step模块中设置Field output中择contact中的FSLIP，这个就是要求输出的接触面的相对滑移的% n0 j1 z2 d6 H0 A在后处理中选择相应的接触对，看结果就可以了/viewthread.php?tid=818579&highlight=%CF%E0%B6%D419. 相对输出响应减去输入即可。, N j! m _$ F |! W* ecae后处理可以通过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* ss,) 2 # 0 L: A这个命令加在以下位置，非常方便了版上很多帖子,受益菲浅!0 V2 # g+ 6 f/ t1 w, c, r5 T关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!1 f1 X9 C! q( u g8 UV5 D1.先施加重力荷载的作用,可以在cae中实现;5 * v6 J3 . . pt2.在inp文件中的output request中写上* * 3 t9 , Xa; o& 1 y$ L+ D*el print- M* i: o* J f, F# H8 C; 4 O9 Ks,- M4 r+ S5 o8 m6 q$ U这样就会将施加重力荷载后的应力输出到*.dat文件中了;S$ * T6 v& F6 . 5 v3 w: T% a5 M! M: z6 l3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如: n% 9 L! + , M4 D7 K3 X0 5 TELEMENTPT 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.% Ol+ 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 pE+ eN8 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- z5.接下来就在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/viewthread.php?tid=724821&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2/viewthread.php?tid=818747&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2边坡的地应力平衡/viewthread.php?tid=782714&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2地应力平衡贴/search.php?searchid=14&orderby=lastpost&ascdesc=desc&searchsubmit=yes&page=419输出桩轴力和弯矩/viewthread.php?tid=770244&highlight=%D6%E1%C1%A6/viewthread.php?tid=856959&highlight=%D6%E1%C1%A6/thread-763434-1-3.html具体的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 Step20关于接触收敛的经验总结/viewthread.php?tid=859723&highlight=%BD%D3%B4%A55 ? L) R6 g# ?4 iu! T8 I; K; J. R9 ?% z6 W, c% _# D# rby hitcondor4 T; ) O1 0 s$ c: yh- ?& xI% Y$ y& t S计算了一个有多处接触的变形问题，经过若干次的尝试终于成功，总结一点经验请高手斧正：. S% S% x( i5 s8 Y: , j5 n; m( 1 P: j5 g$ * t1 首先要判断你的问题是否是大变形问题，如果是，而你的Nlgeom选项又没有打开，那么你将为此付出沉重的代价。& m/ N- |& H# m2 载荷的施加要分多步进行，初始载荷要很小，可以尝试万分之一的量级。( t) Zm7 c7 T2 + W* m- * a) h3 接触建立起来之前，你可能需要增加很多额外的约束，总体的原则是要让接触的区域在初始载荷作用下平稳，不能产生刚体位移（注意：刚体位移远远不仅仅是指整个零件的大幅度移动，范围要宽的多，甚至是节点的稍微大一点点的位移都是）。( 4 R ! c/ e& j: o4 额外约束的去除可能会给你的模型带来较大的约束变动，这使得收敛变得困难，你可以尝试分多步去掉额外约束，绝不要吝惜子步的数量。子步的初始增量步要小，调试阶段你也许可以尝试0.0001，调试后发现收敛容易再放大不迟。1 H t/ w5 y8 O5 有时会发生这样的事情：当你去掉某个额外约束后发现，初始载荷的作用没有到达远离初始载荷位置的接触的地方，这可能使得那个约束并没有建立起来，所以你也许可以尝试在这个状态下再增加初始载荷的大小，让他的作用传播到更远一些的接触的地方。$ U/ R0 a% I/ J7 x$ 6 如果你的初始载荷是力，那么它对于接触的建立可能是一些不利因素，当你尝试过很多次之后若还不收敛，你也许可以把它换成位