标签：red div 学习笔记 当前时间 用两个 nstat ops user 迭代
参考：
1.http://ivt-abaqusdoc.ivt.ntnu.no:2080/v6.14/books/sub/default.htm
2.ABAQUS 用户材料子程序开发及应用-杨曼娟
3.费康老师abaqus书
4.论坛：http://www.mechbbs.com/forum.php?mod=viewthread&tid=104&highlight=abaqus
1.首先看费康老师书第12章邓肯张模型开发。
刚度矩阵就是下面的矩阵。其实刚度矩阵类似于一维里面胡克定律的k（F=k*l 这个胡克定律）
看懂了上面的内容看书上代码就容易的多了。
还需要明白的是，step是分析步，比如加边界条件是第0分析步，开挖第一部分是第1分析步（step1），开挖第二部分是第二分析步（step2）。而umat中有两个时间，一个是Time(1)，代表分析步内当前时间，是以当前分析步开始时间为零点参照。下一个分析步时又会从零开始。而Time(2)是当前总时间，是从step0开始计算的。abaqus在施加力时是一点一点的加上去的，比如某单元上表面施加100kPa压力，且这个步叫做step1，是加载荷步，比如这个分析步中有1000迭代步，那就是每个迭代步增加0.1kPa在单元上表面。
2.再看上面的硕士论文：
umat的一半形式：
1 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 2 1 RPL,DDSDDT,DRPLDE,DRPLDT, 3 2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 4 3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 5 4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC) 6 C 7 INCLUDE ‘ABA_PARAM.INC‘ 8 C 9 CHARACTER*80 CMNAME 10 DIMENSION STRESS(NTENS),STATEV(NSTATV), 11 1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), 12 2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 13 3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3), 14 4 JSTEP(4) 15 16 17 user coding to define DDSDDE, STRESS, STATEV, SSE, SPD, SCD 18 and, if necessary, RPL, DDSDDT, DRPLDE, DRPLDT, PNEWDT 19 20 21 RETURN 22 END
怎样在一个文档里用两个模型呢：
1 IF (CMNAME(1:4) .EQ. ‘MAT1‘) THEN 2 CALL UMAT_MAT1(argument_list) 3 ELSE IF(CMNAME(1:4) .EQ. ‘MAT2‘) THEN 4 CALL UMAT_MAT2(argument_list) 5 END IF
添加材料库：http://www.mechbbs.com/forum.php?mod=viewthread&tid=13
abaqus在线帮助文档：http://www.mechbbs.com/forum.php?mod=viewthread&tid=2
标签：red div 学习笔记 当前时间 用两个 nstat ops user 迭代
原文地址：http://www.cnblogs.com/zhubinglong/p/7765471.html
