SUBROUTINE UMAT(stress,statev,ddsdde,sse,spd,scd, 1 rpl, ddsddt, drplde, drpldt, 2 stran,dstran,time,dtime,temp,dtemp,predef,dpred,cmname, 3 ndi,nshr,ntens,nstatv,props,nprops,coords,drot,pnewdt, 4 celent,dfgrd0,dfgrd1,noel,npt,layer,kspt,kstep,kinc) c WRITE (6,*) ' c NOTE: MODIFICATIONS TO *UMAT FOR ABAQUS VERSION 5.3 (14 APR '94) c c (1) The list of variables above defining the *UMAT subroutine, c and the first (standard) block of variables dimensioned below, c have variable names added compared to earlier ABAQUS versions. c c (2) The statement: include 'aba_param.inc' must be added as below. c c (3) As of version 5.3, ABAQUS files use double precision only. c The file aba_param.inc has a line "implicit real*8" and, since c it is included in the main subroutine, it will define the variables c there as double precision. But other subroutines still need the c definition "implicit real*8" since there may be variables that are c not passed to them through the list or common block. c c (4) This is current as of version 5.6 of ABAQUS. c ... ...

立方晶体，CPFE-UMAT子程序，需要MS+IVF以及Abaqus

构建晶体塑性模型，晶粒之间的动力学关系，可以预测晶体取向、晶粒形态、织构等

详细的晶体塑性有限元理论及算法，固体力学，金属学等研究生，教师可用。国外经典教材。

黄永刚单晶塑性有限元分析umat子程序，abaqus运行的用户材料子程序

通常用规则的六边形或四边形等来表示晶粒，不能够反映出晶界的不规则性，本文使用Voronoi方法在大型有限元软件ABAQUS中建立了多晶材料的几何模型，能够表达出晶粒的几何形状与晶界的不规则性；在此基础上提出了一种控制品粒大小分布及织构的方法，通过调整参数能够建立具有不同晶粒尺寸分布及织构的模型。

对黄永刚老师单晶塑型UMAT子程序进行修改，使之能够进行二维晶体塑性分析

1Dplasticity.pdf

黄永刚单晶塑性有限元umat子程序

晶体塑性有限元经典资料，详细介绍了晶体塑性有限元基础知识。