User login

Navigation

You are here

initialization of initial condition in state variable in vumat

i am still struggling with  defining an initial condition in vumat. i have defined mine like this:

STATEV(1:nstatev)=stateOld(KM,1:nstatev) 

          if (steptime .eq. zero) then

          STATEV(1) = 1.0d0

            STATEV(2) = 1.0d0

            STATEV(3) = 1.0d0

            STATEV(4) = 0.0d0

            STATEV(5) = 0.0d0

            STATEV(6) = 0.0d0

            STATEV(7) = 0.0d0

            end if

      BBAR_DE_N(1,1)= STATEV(1)

      BBAR_DE_N(2,2)= STATEV(2)

      BBAR_DE_N(3,3)= STATEV(3)

      BBAR_DE_N(1,2)= STATEV(4)

      BBAR_DE_N(2,3)= STATEV(5)

      BBAR_DE_N(1,3)= STATEV(6)

      BBAR_DE_N(2,1)=BBAR_DE_N(1,2)

      BBAR_DE_N(3,2)=BBAR_DE_N(2,3)

      BBAR_DE_N(1,3)=BBAR_DE_N(3,1)

BBAR_TRIAL_NN= MATMUL(MATMUL(Fbarrel,BBAR_DE_N),TRANSPOSE(Fbarrel))

 STATEV(1)=BBAR_TRIAL_NN(1,1)

      STATEV(2)=BBAR_TRIAL_NN(2,2)

      STATEV(3)=BBAR_TRIAL_NN(3,3)

      STATEV(4)=BBAR_TRIAL_NN(1,2)

      STATEV(5)=BBAR_TRIAL_NN(2,3)

      STATEV(6)=BBAR_TRIAL_NN(1,3)

     STATEV(7)=Beta_zero

      i dont seem to be able to update bbar_trial_nn. when i write bbar_trial_nn all the terms are zero. 

by the way i also used *initial condition, type=solution

elem,1,1,1,0,0,0,0 with elem being the name of my element.

any help would be highly appreciated. 

Subscribe to Comments for "initialization of initial condition in state variable in vumat"

Recent comments

More comments

Syndicate

Subscribe to Syndicate