功能梯度材料UMAT

chenyue

C    THE ABAQUS USER SUBROUTINE UMAT FOR GRADED ELEMENTS IS PRESENTED HERE. IT IS WRITTEN IN FORTRAN LANGUAGE.
C
C    PROBLEM DESCRIPTION
C    THIS UMAT USER SUBROUTINE IS FOR ISOTROPIC AND ISOTHERMAL MATERIALS.
C    THIS SUBROUTINE HANDLES FOR 3-D AND PLANE STRAIN PROBLEMS.
C    THIS SUBROUTINE NEEDS TO BE MODIFED FOR PLANE STRESS.
C    BOTH HOMOGENEOUS AND GRADED PROPERTIES ARE POSSIBLE BY SELECTING
C    APPROPRIATE FUNCTIONAL FORM OF "PROPS".
C    THIS UMAT IS BASED ON "DIRECT INTEGRATION METHOD".
C    IMPORTANT NOTICE
C    FOR NEW PROBLEMS, CHOOSE THE FUNCTION FOR
C    THE MATERIAL VARIATION AND THE VALUES OF "PROPS".
C    AFTER THAT, GO TO THE SECTION -CHANGE- AND CHANGE THE FUNCTION AND
C    VALUES ACCORDINGLY.

C ------------------------------------------------------------------------------------------------------------------------C
      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
      INCLUDE 'C:\SIMULIA\ABAQUS\6.9-1\SITE\ABA_PARAM_DP.INC'
C
      CHARACTER*8 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
1DDSDDT(NTENS),DRPLDE(NTENS),
     2STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
     3PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
C

C    DETERMINE MATERIAL PROPERTIES BASED ON GLOBAL COORDINATES OF GAUSS POINTS.
C    COORDS (1) IS X-COORDINATE OF GAUSS POINTS.
C    COORDS (2) IS Y-COORDINATE OF GAUSS POINTS.
C    COORDS (3) IS Z-COORDINATE OF GAUSS POINTS.
C    PROPS IS DESIGNED BY USERS.
C    THE FUNCTION CAN BE ALSO DESIGNED BY USERS.

C----CHANGE ------------------------------------------------------------------------------------------------------------------------------------------------------
C    PROPS(1)-----E(0)
C    PROPS(2)-----(E(h)-E(0))/h
C    V=0.3
      E=PROPS(1)+PROPS(2)*(COORDS(2)-1)
C ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

C    DETERMINE LAME'S CONSTANTS
C    MU IS MU
C    ALAMBDA IS LAMBDA
      AMU=E/2.0D0/(1.0D0+0.3)
      ALAMBDA=E*0.3/(1.0D0+0.3)/(1.0D0-2.0D0*0.3)
C    DETERMINE THE TANGENT(JACOBIAN) MATRIX
C    NDI IS NUMBER OF NORMAL STRESSES (E.G. SXX, SYY, SZZ)
C    NSHR IS NUMBER OF SHEAR STRESSES(E.G. SXY)
C    NTENS IS DIMENSION OF CONSTITUTIVE MATRIX (NTENS=NDI+NSHR)
      DO I=1, NDI
      DO J=1, NDI
      IF (I.EQ.J) THEN
C    定义主对角线元素
      DDSDDE(I,I)=ALAMBDA+2.0D0*AMU
      ELSE
      DDSDDE(I,J)=ALAMBDA
      ENDIF
      ENDDO
      ENDDO
      DO I=NDI+1, NTENS
      DDSDDE(I,I)=AMU
      ENDDO
C    DETERMINE THE STRESS AND UPDATE THE STRESS
      DO I=1, NTENS
      DO J=1, NTENS
      STRESS(I)=STRESS(I)+DDSDDE(I,J)*DSTRAN(J)
      ENDDO
      ENDDO
      RETURN
      END

这个子程序试用于动力计算分析吗？希望各位大侠指点！

chenyue

我现在要做关于功能梯度材料的动力分析，但是材料的子程序不正确，希望大侠指点！

zsq-w

lz新来需要先学习下如何发帖
1 已有资料不要在发
2 自己问题不要重发

看看以下2贴
http://forum.simwe.com/thread-867179-1-1.html
http://forum.simwe.com/thread-923999-1-1.html

北鹰南飞

动力计算中，你的材料沿厚度方向的密度怎么处理呢？ 2# chenyue

北鹰南飞

2# chenyue
这段子程序怎么这么眼熟？LZ你是从哪里搞来的？。（当然如果是楼主自创的话就另当别论了，它并不复杂）

北鹰南飞

2# chenyue
这段子程序没有密度的定义，并不符合功能梯度材料的密度分布，LZ觉得做动力计算合适吗？