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

fixed small bug in newstate

This commit is contained in:
Martin Diehl 2014-07-10 08:47:00 +00:00
parent f62ce37ba9
commit 0252fea3d7
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
code/crystallite.f90
View File

@ -978,14 +978,13 @@ subroutine crystallite_stressAndItsTangent(updateJaco,rate_sensitivity)
.and. ((e == debug_e .and. i == debug_i .and. g == debug_g) & .and. ((e == debug_e .and. i == debug_i .and. g == debug_g) &
.or. .not. iand(debug_level(debug_crystallite), debug_levelSelective) /= 0_pInt)) then .or. .not. iand(debug_level(debug_crystallite), debug_levelSelective) /= 0_pInt)) then
if (crystallite_todo(g,i,e)) then if (crystallite_todo(g,i,e)) then
write(6,'(a,f12.8,a,i8,1x,i2,1x,i3)') '<< CRYST >> cutback step in crystallite_stressAndItsTangent & write(6,'(a,f12.8,a,i8,1x,i2,1x,i3,/)') '<< CRYST >> cutback step in crystallite_stressAndItsTangent &
&with new crystallite_subStep: ',& &with new crystallite_subStep: ',&
crystallite_subStep(g,i,e),' at el ip g ',e,i,g crystallite_subStep(g,i,e),' at el ip g ',e,i,g
else else
write(6,'(a,i8,1x,i2,1x,i3)') '<< CRYST >> reached minimum step size & write(6,'(a,i8,1x,i2,1x,i3,/)') '<< CRYST >> reached minimum step size &
&in crystallite_stressAndItsTangent at el ip g ',e,i,g &in crystallite_stressAndItsTangent at el ip g ',e,i,g
endif endif
write(6,*)
endif endif
#endif #endif
endif endif
@ -1054,7 +1053,6 @@ subroutine crystallite_stressAndItsTangent(updateJaco,rate_sensitivity)
call crystallite_integrateStateRKCK45() call crystallite_integrateStateRKCK45()
end select end select
endif endif
where(.not. crystallite_converged .and. crystallite_subStep > subStepMinCryst) & ! do not try non-converged & fully cutbacked any further where(.not. crystallite_converged .and. crystallite_subStep > subStepMinCryst) & ! do not try non-converged & fully cutbacked any further
crystallite_todo = .true. crystallite_todo = .true.
@ -1094,7 +1092,6 @@ subroutine crystallite_stressAndItsTangent(updateJaco,rate_sensitivity)
enddo enddo
enddo elementLooping5 enddo elementLooping5
! --+>> STIFFNESS CALCULATION <<+-- ! --+>> STIFFNESS CALCULATION <<+--
computeJacobian: if(updateJaco) then computeJacobian: if(updateJaco) then
@ -3007,9 +3004,10 @@ subroutine crystallite_integrateStateFPI()
else else
e = eIter(1) e = eIter(1)
i = iIter(1,e) i = iIter(1,e)
do g = iIter(1,e), iIter(2,e) do g = gIter(1,e), gIter(2,e)
p = mappingConstitutive(2,g,i,e) p = mappingConstitutive(2,g,i,e)
c = mappingConstitutive(1,g,i,e) c = mappingConstitutive(1,g,i,e)
plasticState(p)%previousDotState (:,c) = 0.0_pReal plasticState(p)%previousDotState (:,c) = 0.0_pReal
plasticState(p)%previousDotState2(:,c) = 0.0_pReal plasticState(p)%previousDotState2(:,c) = 0.0_pReal
damageState(p)%previousDotState (:,c) = 0.0_pReal damageState(p)%previousDotState (:,c) = 0.0_pReal