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

not used anywhere 

was only required for ductile damage
This commit is contained in:
Martin Diehl 2021-07-04 08:51:53 +02:00
parent 62146d0264
commit 6dc70fa4e7
1 changed files with 2 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/phase_mechanical_plastic_isotropic.f90
View File

@ -31,8 +31,7 @@ submodule(phase:plastic) isotropic
type :: tIsotropicState
real(pReal), pointer, dimension(:) :: &
xi, &
gamma
xi
end type tIsotropicState
!--------------------------------------------------------------------------------------------------
@ -122,7 +121,7 @@ module function plastic_isotropic_init() result(myPlasticity)
!--------------------------------------------------------------------------------------------------
! allocate state arrays
Nmembers = count(material_phaseID == ph)
sizeDotState = size(['xi ','gamma'])
sizeDotState = size(['xi'])
sizeState = sizeDotState
call phase_allocateState(plasticState(ph),Nmembers,sizeState,sizeDotState,0)
@ -135,11 +134,6 @@ module function plastic_isotropic_init() result(myPlasticity)
plasticState(ph)%atol(1) = pl%get_asFloat('atol_xi',defaultVal=1.0_pReal)
if (plasticState(ph)%atol(1) < 0.0_pReal) extmsg = trim(extmsg)//' atol_xi'
stt%gamma => plasticState(ph)%state (2,:)
dot%gamma => plasticState(ph)%dotState(2,:)
plasticState(ph)%atol(2) = pl%get_asFloat('atol_gamma',defaultVal=1.0e-6_pReal)
if (plasticState(ph)%atol(2) < 0.0_pReal) extmsg = trim(extmsg)//' atol_gamma'
end associate
!--------------------------------------------------------------------------------------------------
@ -285,8 +279,6 @@ module subroutine isotropic_dotState(Mp,ph,en)
dot%xi(en) = 0.0_pReal
endif
dot%gamma(en) = dot_gamma ! ToDo: not really used
end associate
end subroutine isotropic_dotState