FEAP User Forum
FEAP => Programming => Topic started by: 36051121 on September 11, 2014, 02:28:43 AM
-
Dear ALL,
I am using user element, which has 3 degrees of freedom, corresponding to ul(1,nen), ul(2,nen), and ul(3,nen). In FEAP, the values of ul(ndf, nen) are output automatically and act as the initial values for the next time step. However, in my problem, I must normalize ul(ndf,nen) by
ul*(j,nen)=ul(j,nen)/sqrt(ul(1,nen)^2+ul(2,nen)^2+ul(3,nen)^2),
and use ul*(j,nen) as the initial values for the next time step.
So my question is that how I can normalize ul(j,nen) as ul*(j,nen), output ul*(j,nen), and automatically transform ul*(j,nen) as the initial values for the next time step ?
Many thanks !
-
Look in program/plocal.f but beaware that plocal is called for every iteration within a time step.