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

output was screwed up

This commit is contained in:
Martin Diehl 2018-08-26 06:20:48 +02:00
parent 83740b5d7b
commit 0041d21777
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/plastic_phenopowerlaw.f90
View File

@ -666,14 +666,14 @@ subroutine plastic_phenopowerlaw_dotState(Mstar6,ipc,ip,el)
tau_slip_pos = tau_slip_pos + math_mul33xx33(Mstar,prm%nonSchmid_pos(1:3,1:3,k,j)) tau_slip_pos = tau_slip_pos + math_mul33xx33(Mstar,prm%nonSchmid_pos(1:3,1:3,k,j))
tau_slip_neg = tau_slip_neg + math_mul33xx33(Mstar,prm%nonSchmid_neg(1:3,1:3,k,j)) tau_slip_neg = tau_slip_neg + math_mul33xx33(Mstar,prm%nonSchmid_neg(1:3,1:3,k,j))
enddo nonSchmidSystems enddo nonSchmidSystems
gdot_slip(j) = prm%gdot0_slip*0.5_pReal* & gdot_slip(j) = prm%gdot0_slip*0.5_pReal* & !ToDo: save to dotState
( (abs(tau_slip_pos)/(stt%s_slip(j,of)))**prm%n_slip*sign(1.0_pReal,tau_slip_pos) & ( (abs(tau_slip_pos)/(stt%s_slip(j,of)))**prm%n_slip*sign(1.0_pReal,tau_slip_pos) &
+(abs(tau_slip_neg)/(stt%s_slip(j,of)))**prm%n_slip*sign(1.0_pReal,tau_slip_neg)) +(abs(tau_slip_neg)/(stt%s_slip(j,of)))**prm%n_slip*sign(1.0_pReal,tau_slip_neg))
enddo enddo
do j = 1_pInt, prm%totalNtwin do j = 1_pInt, prm%totalNtwin
tau_twin = math_mul33xx33(Mstar,prm%Schmid_twin(1:3,1:3,j)) tau_twin = math_mul33xx33(Mstar,prm%Schmid_twin(1:3,1:3,j))
gdot_twin(j) = (1.0_pReal-stt%sumF(of))*prm%gdot0_twin* (abs(tau_twin)/stt%s_twin(j,of))**prm%n_twin & gdot_twin(j) = (1.0_pReal-stt%sumF(of))*prm%gdot0_twin* (abs(tau_twin)/stt%s_twin(j,of))**prm%n_twin & !ToDo: save to dotState
* max(0.0_pReal,sign(1.0_pReal,tau_twin)) * max(0.0_pReal,sign(1.0_pReal,tau_twin))
enddo enddo
@ -733,7 +733,7 @@ function plastic_phenopowerlaw_postResults(Mstar6,ipc,ip,el)
plastic_phenopowerlaw_postResults plastic_phenopowerlaw_postResults
integer(pInt) :: & integer(pInt) :: &
ph, of, & of, &
o,c,j,k o,c,j,k
real(pReal) :: & real(pReal) :: &
tau_slip_pos,tau_slip_neg,tau_twin tau_slip_pos,tau_slip_neg,tau_twin
@ -741,6 +741,7 @@ function plastic_phenopowerlaw_postResults(Mstar6,ipc,ip,el)
type(tParameters) :: prm type(tParameters) :: prm
type(tPhenopowerlawState) :: stt type(tPhenopowerlawState) :: stt
of = phasememberAt(ipc,ip,el)
associate( prm => param(phase_plasticityInstance(material_phase(ipc,ip,el))), & associate( prm => param(phase_plasticityInstance(material_phase(ipc,ip,el))), &
stt => state(phase_plasticityInstance(material_phase(ipc,ip,el))) ) stt => state(phase_plasticityInstance(material_phase(ipc,ip,el))) )
Mstar = math_Mandel6to33(Mstar6) Mstar = math_Mandel6to33(Mstar6)