From 872228dd47e42d75e405249ec57797248c976b61 Mon Sep 17 00:00:00 2001 From: achalhp Date: Wed, 6 Mar 2024 11:37:41 +0530 Subject: [PATCH] frozen in kinetics_tw --- src/phase_mechanical_plastic_phenopowerlaw.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/phase_mechanical_plastic_phenopowerlaw.f90 b/src/phase_mechanical_plastic_phenopowerlaw.f90 index 2c93774fe..058e06933 100644 --- a/src/phase_mechanical_plastic_phenopowerlaw.f90 +++ b/src/phase_mechanical_plastic_phenopowerlaw.f90 @@ -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