umat自学整理.docx

什么是UMAT?v. l4 W* F9 T% rUMAT:[user defined material mechanical behavior]用户材料子程序,是ABAQUS提供给用户自定义材料属性的FORTRAN程序接口;注意是接口程序,真正的定义材料的力学行为即属性,是用户自己编译的FORTRAN程序来实现的UMAT通过与ABAQUS主求解程序的接口实现与ABAQUS的数据交流+ ^7 P Z1 ?. UMAT功能简介![摘自庄茁老师的书]UMAT子程序具有强大的功能，使用UMAT子程序：(1)可以定义材料的本构关系，使用ABAQUS材料库中没有包含的材料进行计算，扩充程序功能。ABAQUS软件2003年度用户年会论文集, I* A: (2)几乎可以用于力学行为分析的任何分析过程，几乎可以把用户材料属性赋予ABAQUS中的任何单元；1 K ? e# g3 l/ I! f/ X2 \(3)必须在UMAT中提供材料本构模型的雅可比（Jacobian）矩阵，即应力增量对应变增量的变化率。) z: w# Q. H/ W0 d+ M* J; f(4)可以和用户子程序“USDFLD”联合使用，通过“USDFLD”重新定义单元每一物质点上传递到UMAT中场变量的数值。UMAT开始的变量声明6 p4 q/ a C4 S- T由于主程序与UMAT之间存在数据传递，甚至共用一些变量，因此必须遵守有关书写格式，UMAT中常用的变量在文件开头予以定义，通常格式为： Q9 ?+ o6 Y- w3 Q4 i: h) CSUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,1 RPL,DDSDDT,DRPLDE,DRPLDT, U# \. h: W6 ?( w/ X+ {$ C; M2STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME3 u. _1 T v5 y7 J0 y Y) z8 K3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) E1 o B4 H! Y- y6 R1 F1 wCINCLUDEABA_PARAM.INC-----此处是将ABAQUS本身自带的参量精度定义的文件包含进来[后面详说]CCHARACTER*80 CMNAMEDIMENSION STRESS(NTENS),STATEV(NSTATV),1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)-------------------------------以上是变量声明吧!user coding to define DDSDDE,STRESS,STATEV,SSE,SPD,SCDand,if necessary,RPL,DDSDDT,DRPLDE,DRPLDT,PNEWDT( j9 l W( I n u-------------------------------此处,看来是将用户定义材料属性的fortran程序编入 J. O* d0 B! ORETURN------------------这是返回值END------------------------结束UMAT中各个变量的详细解释[凡是-注明david的,都是我自己猜的,仅供参考]DDSDDE (NTENS ,NTENS)是一个NTENS[Number of the Tensions----david]维的方阵，称作雅可比矩阵，应力增量/应变增量的偏导数，DDSDDE (I ,J)表示增量步结束时第J个应变分量的改变引起的第I个应力增量的变化!雅可比是一个对称矩阵，除非在“*USER MATERIAL”语句中加UNSYMM参数+ z) D0 \/ c$ N, E+ f2 x2 N) O Q. ` ~# S! B) |3 C1 E3 @6 d( Z1 I ]/ w6 ISTRESS (NTENS)1 t5 A, q- R, @* N4 K应力张量矩阵，对应NDI[Number of the Direct Components--david]个直接分量和NSHR[Number of the shear Components-david]个剪切分量.在增量步的开始,应力张量矩阵中的数值通