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

Cannot write out resolved stress including non schmid 

Would need to write out tau_pos and tau_neg separately, so limit to
resolved stress in the sense of Schmid
This commit is contained in:
Martin Diehl 2018-10-05 07:26:13 +02:00
parent e74a852ef7
commit 47a748831f
1 changed files with 1 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/plastic_phenopowerlaw.f90
View File

@ -769,8 +769,6 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults)
integer(pInt) :: & integer(pInt) :: &
o,c,i,j o,c,i,j
real(pReal) :: &
tau_slip_pos, tau_slip_neg
real(pReal), dimension(param(instance)%totalNslip) :: & real(pReal), dimension(param(instance)%totalNslip) :: &
gdot_slip_pos,gdot_slip_neg gdot_slip_pos,gdot_slip_neg
@ -797,13 +795,7 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults)
c = c + prm%totalNslip c = c + prm%totalNslip
case (resolvedstress_slip_ID) case (resolvedstress_slip_ID)
do i = 1_pInt, prm%totalNslip do i = 1_pInt, prm%totalNslip
tau_slip_pos = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) postResults(c+i) = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i))
tau_slip_neg = tau_slip_pos
!do j = 1,size(prm%nonSchmidCoeff)
! tau_slip_pos = tau_slip_pos + math_mul33xx33(S,prm%nonSchmid_pos(1:3,1:3,j,i))
! tau_slip_neg = tau_slip_neg + math_mul33xx33(S,prm%nonSchmid_neg(1:3,1:3,j,i))
!enddo
postResults(c+i) = 0.5_pReal*(tau_slip_pos+tau_slip_neg)
enddo enddo
c = c + prm%totalNslip c = c + prm%totalNslip