UMATNeoHookean.docVIP

UMAT 分析Neo-Hookean材料 !由于本人学ABA不久, 到现在为止, 对umat还是处于半懂状态, 从ABA论坛下载了大大们提供的资料和读了一些大大的经验之谈, 受益非浅! 以下这段是我关于neo-hookean材料子程序的部分阐释,不当之处还望各为大大指正.(由于本人水平有限,但本着大家互相帮助的出发点,所以抛点砖头出来,^_^) Neo-Hookean材料模型是各向同性线性定律(Hookean定律)至大变形的扩展. 在关键词*ELASTIC中定义弹性模量不能有效的解决有限变形问题,必须定义一个合适的势能函数 SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL, 1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, TIME, DTIME, TEMP, DTEMP, 2 PREDEF, DPRED, CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, NPROPS, 3 COORDS, DROT, PNEWDT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, LAYER, 4 KSPT, KSTEP, KINC) C INCLUDE ’ABA_PARAM.INC’ C CHARACTER*8 CMNAME C DIMENSION STRESS(NTENS), STATEV(NSTATV), DDSDDE(NTENS, NTENS), 1 DDSDDT(NTENS), DRPLDE(NTENS), STRAN(NTENS), DSTRAN(NTENS), 2 PREDEF(1), DPRED(1), PROPS(NPROPS), COORDS(3), DROT(3, 3), 3 DFGRD0(3, 3), DFGRD1(3, 3) C LOCAL ARRAYS C ---------------------------------------------------------------- C EELAS - LOGARITHMIC ELASTIC STRAINS C EELASP - PRINCIPAL ELASTIC STRAINS C BBAR - DEVIATORIC RIGHT CAUCHY-GREEN TENSOR RIGHT 应为LEFT?? C BBARP - PRINCIPAL VALUES OF BBAR C BBARN - PRINCIPAL DIRECTION OF BBAR (AND EELAS) C DISTGR - DEVIATORIC DEFORMATION GRADIENT (DISTORTION TENSOR) C ---------------------------------------------------------------- C DIMENSION EELAS(6), EELASP(3), BBAR(6), BBARP(3), BBARN(3, 3), 1 DISTGR(3,3) C PARAMETER(ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, FOUR=4.D0, 1 SIX=6.D0) C C ---------------------------------------------------------------- C UMAT FOR COMPRESSIBLE NEO-HOOKEAN HYPERELASTICITY C CANNOT BE USED FOR PLANE STRESS C ---------------------------------------------------------------- C PROPS(1) - E C PROPS(2) – NU C ---------------------------------------------------------------- C C ELASTIC PROPERTIES C EMOD=PROPS(1) ENU=PROPS(2) C10=EMOD/(FOUR*(ONE+ENU)) D1=SIX*(ONE-TWO*ENU)/EMOD C C JACOBIAN AND DISTORTION TENSOR C 1.程序里面的DET=,当平面应变问题时,=0 2. DISTGR=Fdev=J-1/3.F 为偏量变形梯度, 体积变形梯度 FV0L=J1/3.I, 且FVOL.Fdev=F DET=DFGRD1(1, 1)*DFGRD1(2, 2)*DFGRD1(3, 3) 1