FEAP User Forum

FEAP => Programming => Topic started by: Aamir Dean on June 11, 2015, 09:55:13 AM

Title: Initialization of history variables in materials routines (UMATL)
Post by: Aamir Dean on June 11, 2015, 09:55:13 AM

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,

Title: Re: Initialization of history variables in materials routines (UMATL)
Post by: Prof. R.L. Taylor 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

Title: Re: Initialization of history variables in materials routines (UMATL)
Post by: Aamir Dean 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 ?

Title: Re: Initialization of history variables in materials routines (UMATL)
Post by: Prof. R.L. Taylor on June 12, 2015, 03:31:13 AM

Sorry, yes 1.0d0