30-7-2020

*.avi

*.pyc

*.pdf

import connectorBehavior

a='/home'

b='apadhy'

c='tryout'

c='projects'

d='PhD'

e='ABAQUS'

f='UMAT_implementation'

!DIR$ NOFREEFORM

      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,JSTEP,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),

     4 JSTEP(4)

C     ELASTIC USER SUBROUTINE

      PARAMETER (ONE=1.0D0, TWO=2.0D0)

	  E=PROPS(1)

	  ANU=PROPS(2)

      ALAMBDA=E/(ONE+ANU)/(ONE-TWO*ANU)

	  BLAMBDA=(ONE-ANU)

         CLAMBDA=(ONE-TWO*ANU)

	     DO I=1,NTENS

	      DO J=1,NTENS

	      DDSDDE(I,J)=0.0D0

	      ENDDO

	     ENDDO

		DDSDDE(1,1)=(ALAMBDA*BLAMBDA)

		DDSDDE(2,2)=(ALAMBDA*BLAMBDA)

		DDSDDE(3,3)=(ALAMBDA*BLAMBDA)

		DDSDDE(4,4)=(ALAMBDA*CLAMBDA)

		DDSDDE(5,5)=(ALAMBDA*CLAMBDA)

		DDSDDE(6,6)=(ALAMBDA*CLAMBDA)

		DDSDDE(1,2)=(ALAMBDA*ANU)

		DDSDDE(1,3)=(ALAMBDA*ANU)

		DDSDDE(2,3)=(ALAMBDA*ANU)

		DDSDDE(2,1)=(ALAMBDA*ANU)

		DDSDDE(3,1)=(ALAMBDA*ANU)

		DDSDDE(3,2)=(ALAMBDA*ANU)

         DO I=1,NTENS

	      DO J=1,NTENS

	      STRESS(I)=STRESS(I)+DDSDDE(I,J)*DSTRAN(J)

	      ENDDO

	     ENDDO

      RETURN

      END

geometry:

    film_length : 1.0 #x length of film

    film_height : 0.5 # y height of film

    film_width : 0.5 #z width of film

    mesh: 0.25               #Mesh size in mm

UMAT: "3D_SMA_UM.f"

materials:

    Density: 6450

    Depvar: 100

    Conductivity: 22.0

    Specific_heat: 329.0

    User_material:

        # E: 70e9

        # NU: 0.33

        IPHASE: 1   # 1: Austenite, 2: Martensite

        MODEL: 2  # Constitutive model 1: Tanaka's exponential, 2: Bo and Lagoudas, 3: Liand and Roger's cosine model

        TOL: 1.0e-8  # Convergence Tolerance

        xi0: 0.0  #Initial volume of the martensitic volume fraction

        NELMTP: 16  #Number of integration points in all SMA finite elements

        EA: 70e9   # Young's modulus for austenite, Pa

        EM: 30e9 # YOung's modulus for martensite, Pa

        nu: 0.33 # Poissons ratio

        alphaA: 22e-6    #Thermal expansion coefficient of austenite K-1

        alphaM:  22e-6  #Thermal expansion coefficient of martensite K-1

        Mos: 291.0   #Martensite start temperature

        Mof: 271.0   #Martensite finish temperature

        Aos: 295.0  #Austenite start temperature

        Aof: 315.0   #Austenite finish temperature

        H: 0.05    #Maximum current transformation strains

        rSoMH: -0.35e6     #Martensite stress influence coefficient, Pa/K

        rSoAH:  -0.35e6    #Austenite stress influence coefficient, Pa/K

        epstr11: 0.0  #Initial value of transformation strain tensor 11 component

        epstr22: 0.0   #Initial value of transformation strain tensor 22 component

        epstr33: 0.0   #Initial value of transformation strain tensor 33 component

        2epstr23: 0.0    #2X Initial value of transformation strain tensor 23 component

        2epstr13: 0.0   #2X Initial value of transformation strain tensor 13 component

        2epstr12: 0.0   #2X Initial value of transformation strain tensor 12 component

        FRULE: 1  #Flag for the form of the transformation tensor 1- Direction dependent 2- Direction independent

#Loading condition in mm,N

loading:

    init_temp: 0.0

    steps: 2  #Number of steps

    type: 1 #1- Force, 2 - Displacement, 3 - Temperature

    value: 150000