abaqus常用技巧总结 - 图文(12)

: 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) c

include 'aba_param.inc' c

character*80 ciname,slname,msname

dimension tau(nfdir),ddtddg(nfdir,nfdir),ddtddp(nfdir), 1 dslip(nfdir),ddtddt(nfdir,2),statev(*),

2 dgam(nfdir),taulm(nfdir),slip(nfdir),time(2), 3 coords(mcrd),rcoord(mcrd),drot(2,2),temp(2), 4 predef(2,*),props(nprops) c

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

if (lm .eq. 2) return c

a=props(1) b=props(2) c

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) C

INCLUDE 'ABA_PARAM.INC' C

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(*) c c

FIELD(1)=time(2) c print*,field(1) c c

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. …… 此处隐藏：1509字，全部文档内容请下载后查看。喜欢就下载吧 ……