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

bugfix: same tau was used for all systems

This commit is contained in:
Martin Diehl 2018-10-11 17:28:00 +02:00
parent 3a39c2d68c
commit 51a3308a3d
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/plastic_dislotwin.f90
View File

@ -1269,11 +1269,9 @@ subroutine plastic_dislotwin_dotState(Mp,Temperature,instance,of)
- sum(stt%stressTransFraction(1_pInt:prm%totalNtrans,of)) &
- sum(stt%strainTransFraction(1_pInt:prm%totalNtrans,of))
slipState1: do i = 1_pInt, prm%totalNslip
tau = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i))
enddo slipState1
call kinetics_slip(prm,stt,mse,of,Mp,temperature,gdot_slip)
slipState2: do i = 1_pInt, prm%totalNslip
slipState: do i = 1_pInt, prm%totalNslip
tau = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i))
DotRhoMultiplication = abs(gdot_slip(i))/(prm%burgers_slip(i)*mse%mfp_slip(i,of))
EdgeDipMinDistance = prm%CEdgeDipMinDistance*prm%burgers_slip(i)
@ -1315,7 +1313,7 @@ subroutine plastic_dislotwin_dotState(Mp,Temperature,instance,of)
dot%rhoEdge(i,of) = DotRhoMultiplication-DotRhoDipFormation-DotRhoEdgeEdgeAnnihilation
dot%rhoEdgeDip(i,of) = DotRhoDipFormation-DotRhoEdgeDipAnnihilation-DotRhoEdgeDipClimb
dot%accshear_slip(i,of) = abs(gdot_slip(i))
enddo slipState2
enddo slipState
twinState: do i = 1_pInt, prm%totalNtwin