223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

state = subState0 + doState * dt 

this is the usual state update used in all other integrators. also
in-line with logic in crystallite_stress
This commit is contained in:
Martin Diehl 2019-01-22 23:23:48 +01:00
parent 5ed3707013
commit ac9d49f6be
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/crystallite.f90
View File

@ -2102,13 +2102,13 @@ eIter = FEsolving_execElem(1:2)
c = phasememberAt(g,i,e)
mySizePlasticDotState = plasticState(p)%sizeDotState
plasticState(p)%state( 1:mySizePlasticDotState,c) = &
plasticState(p)%state( 1:mySizePlasticDotState,c) &
plasticState(p)%subState0( 1:mySizePlasticDotState,c) &
+ plasticState(p)%dotState(1:mySizePlasticDotState,c) &
* crystallite_subdt(g,i,e)
do mySource = 1_pInt, phase_Nsources(p)
mySizeSourceDotState = sourceState(p)%p(mySource)%sizeDotState
sourceState(p)%p(mySource)%state( 1:mySizeSourceDotState,c) = &
sourceState(p)%p(mySource)%state( 1:mySizeSourceDotState,c) &
sourceState(p)%p(mySource)%subState0( 1:mySizeSourceDotState,c) &
+ sourceState(p)%p(mySource)%dotState(1:mySizeSourceDotState,c) &
* crystallite_subdt(g,i,e)
enddo