Author Topic: Regarding storage of History variables  (Read 6529 times)

ad10

  • Jr. Member
  • **
  • Posts: 39
Regarding storage of History variables
« on: July 06, 2021, 01:16:49 AM »
Dear all

I am formulating a User element in which I required to store previous step strains

Though I am assigning the history variable hr(nh1) for the previous and hr(nh2) for the current step. But, during simulation values are not getting stored!!  any suggestions about what probably going wrong?

do i = 1,3
                  en(i) = hr(nh1+nn)  ! getting previously stored strains
 end do
And updating the history with,

 do i = 1,3
                  hr(nh2+nn)     = ee(i)
end do

For each step, I am getting en as a zero vector.

Thanks in advance!

JStorm

  • Sr. Member
  • ****
  • Posts: 250
Re: Regarding storage of History variables
« Reply #1 on: July 06, 2021, 01:51:10 AM »
Do you want to have the values of the last time step (equilibrium step) or from the last Newton iteration?

You are right with your procedure if you need the values of the last time step.
Maybe you have missed to inform FEAP about the required number of history components?
It is explained in pmanual85.pdf in chapter 5.7.1.

If you want to store your variable in each iteration, then look for nh3 in the manual.

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: Regarding storage of History variables
« Reply #2 on: July 06, 2021, 04:39:24 AM »
Set nh1 = 3 in isw=1 part of element

Should be hr(nh1+nn-1) = ee(nn). Etc

ad10

  • Jr. Member
  • **
  • Posts: 39
Re: Regarding storage of History variables
« Reply #3 on: July 21, 2021, 02:32:58 PM »
Thanks for replying back.

I have formulated an element for Linear Viscoelastic material,

The strange issue I am facing here,
My element works very perfectly for single-arm but the solution unconverges (UN) as I introduce another arm(Increases the viscoelastic terms nv)

I am thinking the issue is with storing the qi(1:ntm,2) values because I am getting results for qi(1:ntm,1).

I am storing the elastic strains and shear strains in ISW 3 or 6 only. I have added my element here. Can you please comment?

Thanks in advance

« Last Edit: July 26, 2021, 10:37:31 AM by ard10 »

Prof. S. Govindjee

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 1164
Re: Regarding storage of History variables
« Reply #4 on: July 21, 2021, 09:52:10 PM »
nh1 should be set to the number of (scalar) history variables times the number of element gauss points.

So for example if you are storing a single symmetric tensor per arm of the VE model, and you had 4 Gauss points in the element, then nh1 = 6 * 4 * number_arms_in_VE_model.

ad10

  • Jr. Member
  • **
  • Posts: 39
Re: Regarding storage of History variables
« Reply #5 on: July 25, 2021, 08:29:16 AM »
Thanks a lot Professor!!

I able to formulate one. I am making mistake in sizing of nh1.

ad10

  • Jr. Member
  • **
  • Posts: 39
Re: Regarding storage of History variables
« Reply #6 on: November 03, 2021, 08:04:14 AM »
Dear all,

I was able to use nh1: previous step variables at step (tn) and nh2: current step variables (tn+1). I want to know, can I get data of all previous values (time data) of variables available for calculating a value of the variable for the current step. Is that possible to do?

Prof. S. Govindjee

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 1164
Re: Regarding storage of History variables
« Reply #7 on: November 03, 2021, 10:40:44 AM »
We only store t_n and t_n+1.  If you would like to have t_{n-1}, t_{n-2}, etc., then you will need to store it yourself.  Two obvious options are to write a macro to write the data to disk or to write a macro to store the data in memory.

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: Regarding storage of History variables
« Reply #8 on: November 04, 2021, 11:34:07 AM »
Storing ALL previous values is not a good idea as it would involve too much data storage for a problem with many elements.   We always try to do solutions from t_n to t_n+1.