Author Topic: Initialization of history variables in materials routines (UMATL)  (Read 4613 times)

Aamir Dean

  • Jr. Member
  • **
  • Posts: 15
Dear All,
I hope this is not a repeated topic. I actually went through all the topics that discuss history variables initialization, but I just couldn't find an answer to my question since all these topics are related to initialization of history variables in Element routines.
Anyway, I'm just wondering about how to initialize history variables in Material routines?
Is it simply just I have to add these lines at the beginning of the umatl routine if I would like to initialize the first history variable to 1.0 for example?
c... INITIALIZE HISTORY VARIABLES
        IF (ISW.EQ. 14) THEN
           HN(1) = 1.0D0
        END IF

Best regards,

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: Initialization of history variables in materials routines (UMATL)
« Reply #1 on: June 11, 2015, 01:39:08 PM »
Your form should work,

You could make it
     if(isw.eq.14) then
       hn(1) = 0.0d0
       h1(1) = 0.0s0   ! For safety
     else
       ...  compute your model here
     endif

Aamir Dean

  • Jr. Member
  • **
  • Posts: 15
Re: Initialization of history variables in materials routines (UMATL)
« Reply #2 on: June 11, 2015, 03:07:15 PM »
Many thanks dear Prof. Tylor for the quick reply. But do you mean
       hn(1) = 0.0d0
       h1(1) = 0.0d0 
or (I suppose)
       hn(1) = 1.0d0
       h1(1) = 1.0d0 ?

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: Initialization of history variables in materials routines (UMATL)
« Reply #3 on: June 12, 2015, 03:31:13 AM »
Sorry, yes 1.0d0