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

non-Schmid contribution were ignored in test 

reverted, just to ensure that the results have not changed: Never
polishing and change physics at the same time
This commit is contained in:
Martin Diehl 2018-09-14 05:21:05 +02:00
parent 8ecb019566
commit b165e48d9f
1 changed files with 4 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/plastic_phenopowerlaw.f90
View File

@ -699,7 +699,6 @@ subroutine kinetics_slip(prm,stt,of,S,gdot_slip_pos,gdot_slip_neg, &
enddo enddo
enddo enddo
gdot_slip_pos = 0.5_pReal*prm%gdot0_slip & gdot_slip_pos = 0.5_pReal*prm%gdot0_slip &
* sign(abs(tau_slip_pos/stt%s_slip(:,of))**prm%n_slip, tau_slip_pos) * sign(abs(tau_slip_pos/stt%s_slip(:,of))**prm%n_slip, tau_slip_pos)
gdot_slip_neg = 0.5_pReal*prm%gdot0_slip & gdot_slip_neg = 0.5_pReal*prm%gdot0_slip &
@ -828,10 +827,10 @@ function plastic_phenopowerlaw_postResults(Mstar6,ipc,ip,el) result(postResults)
do i = 1_pInt, prm%totalNslip do i = 1_pInt, prm%totalNslip
tau_slip_pos = math_mul33xx33(S,prm%Schmid_slip(1:3,1:3,i)) tau_slip_pos = math_mul33xx33(S,prm%Schmid_slip(1:3,1:3,i))
tau_slip_neg = tau_slip_pos tau_slip_neg = tau_slip_pos
do j = 1,size(prm%nonSchmidCoeff) !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_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)) ! tau_slip_neg = tau_slip_neg + math_mul33xx33(S,prm%nonSchmid_neg(1:3,1:3,j,i))
enddo !enddo
postResults(c+i) = 0.5_pReal*(tau_slip_pos+tau_slip_neg) postResults(c+i) = 0.5_pReal*(tau_slip_pos+tau_slip_neg)
enddo enddo
c = c + prm%totalNslip c = c + prm%totalNslip