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.