umat到vumat转化程序_5.本构模型-应力更新专题-UMAT和VUMAT资料.pptx

计算固体力学应力更新专题

柳占立 庄茁liuzhanli@tsinghua.edu.cn ;导出有限元方程(完全，更新、共旋);应力更新;;;;;2.几种客观率;2.几种客观率;例4.1; 对于不同客观率采用了相同的材料常数，其差别是非常大的。事实上，这是误用了材料模型。材料模型必须根据不同的率转换。这是变形体，若是刚体转动，Jaumann率与Truesdell率是一致的。;Green-Naghdi率;切线模量之间的关系 ;2.几种客观率的关系;2.几种客观率的关系;2.几种客观率的关系;Comparison of different objective stress rate;2.几种客观率的关系;Solver;后拉、前推;后拉、前推;Lie导数 ;后拉、前推和Lie导数 ;4.常用客观率应力更新流程;速度梯度;4.常用客观率应力更新流程;4.常用客观率应力更新流程;刚体转动求解;4.常用客观率应力更新流程;4.常用客观率应力更新流程;4.常用客观率应力更新流程;4.常用客观率应力更新流程;1 ABAQUS的UMAT;ABAQUS用户子程序; 由于主程序与UMAT之间存在数据传递，甚至共用一些变量，因此必须遵守有关UMAT的书写格式，UMAT中常用的变量在文件开头予以定义，通常格式为：?SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,DRPLDE, 1 DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 2 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 3 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)C INCLUDE 'ABA_PARAM.INC'C CHARACTER*80 CMNAME DIMENSION 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, SCD and, if necessary, RPL, DDSDDT, DRPLDE, DRPLDT, PNEWDT? RETURN END;; UMAT中的应力矩阵、应变矩阵以及矩阵 DDSDDE，DDSDDT，DRPLDE 等，都是直接分量存储在前，剪切分量存储在后。直接分量有 NDI个，剪切分量有 NSHR 个。各分量之间的顺序根据单元自由度的不同有一些差异，所以编写UMAT时要考虑到所使用单元的类别。;用于存储状态变量的矩阵，在增量步开始时将数值传递到UMAT中。也可在子程序USDFLD或UEXPAN中先更新数据，然后在增量步开始时将??新后的数据传递到UMAT中。在增量步结束时必须更新状态变量矩阵中的数据。 和应力张量矩阵不同的是：对于有限应变问题，除了材料本构行为引起的数据更新以外，状态变量矩阵NSTATEV中的任何矢量或者张量都必须通过旋转来考虑材料的刚体运动。;2 ABAQUS的VUMAT;2 ABAQUS的VUMAT;UMat:For geometrically nonlinear analysis the strain increment and incremental rotation passed into the routine are based on the Hughes-Winget formulae.– Linearized strain and rotation increments are calculated in the mid-increment con?guration.– Approximations are made, particularly if rotation increments are large: more accurate measures can be obtained