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

frozen in kinetics_tw

This commit is contained in:
achalhp 2024-03-06 11:37:41 +05:30
parent e38ad80077
commit 872228dd47
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/phase_mechanical_plastic_phenopowerlaw.f90
View File

@ -553,6 +553,7 @@ associate(prm => param(ph), stt => state(ph), dot => dotState(ph), dlt => deltas
Success_Nucleation: if (random <= stt%f_twin(twin_var,en)) then ! Instead of sum take max
twinJump = .true.
deltaFp = prm%CorrespondanceMatrix(:,:,twin_var)
write(6,*)'deltaFp',deltaFp,'element',en
dlt%f_twin(:,en) = 0.0_pReal - stt%f_twin(:,en)
dlt%fmc_twin(:,en) = 0.0_pReal - stt%fmc_twin(:,en)
dlt%frozen(en) = 1.0_pReal - stt%frozen(en)
@ -739,7 +740,7 @@ associate(prm => param(ph), stt => state(ph))
tau_tw = [(math_tensordot(Mp,prm%P_tw(1:3,1:3,i)),i=1,prm%sum_N_tw)]
where(tau_tw > 0.0_pReal)
where(tau_tw > 0.0_pReal .and. stt%frozen(en) < 0.9_pReal)
dot_gamma_tw = (1.0_pReal-sum(stt%gamma_tw(:,en)/prm%gamma_char)) & ! only twin in untwinned volume fraction
* prm%dot_gamma_0_tw*(abs(tau_tw)/stt%xi_tw(:,en))**prm%n_tw
fdot_twin = (0.05_pReal*(abs(tau_tw)/stt%xi_tw(:,en))**prm%n_tw)/prm%gamma_char !Achal 0.05 is constant