From 28ea09050d806b1a6cada0bcc5e0bb4aaafde3a2 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 25 Oct 2022 23:18:22 +0200 Subject: [PATCH] polishing --- .../config/phase/mechanical/plastic/dislotwin_TWIP-TRIP.yaml | 4 ++-- src/phase_mechanical_plastic_dislotwin.f90 | 4 ++++ 2 files changed, 6 insertions(+), 2 deletions(-) diff --git a/examples/config/phase/mechanical/plastic/dislotwin_TWIP-TRIP.yaml b/examples/config/phase/mechanical/plastic/dislotwin_TWIP-TRIP.yaml index 962f157cb..3e6087140 100644 --- a/examples/config/phase/mechanical/plastic/dislotwin_TWIP-TRIP.yaml +++ b/examples/config/phase/mechanical/plastic/dislotwin_TWIP-TRIP.yaml @@ -33,7 +33,7 @@ h_sl-sl: [0.122, 0.122, 0.625, 0.07, 0.137, 0.137, 0.122] # Twin N_tw: [12] -b_tw: [1.47e-10] # a/sqrt(6) +b_tw: [1.47e-10] # a_cF/sqrt(6) L_tw: 1.91e-7 # 1300 *b_tw i_tw: 10.0 t_tw: [5.0e-8] @@ -44,7 +44,7 @@ h_sl-tw: [1.0, 1.0, 1.0] # Transformation N_tr: [12] -b_tr: [1.47e-10] # a/sqrt(6) +b_tr: [1.47e-10] # a_cF/sqrt(6) L_tr: 2.21e-7 # 1500 *b_tr i_tr: 10.0 # adjusted t_tr: [1.0e-7] diff --git a/src/phase_mechanical_plastic_dislotwin.f90 b/src/phase_mechanical_plastic_dislotwin.f90 index 984992ba7..ffa6d50b9 100644 --- a/src/phase_mechanical_plastic_dislotwin.f90 +++ b/src/phase_mechanical_plastic_dislotwin.f90 @@ -371,6 +371,10 @@ module function plastic_dislotwin_init() result(myPlasticity) if (prm%fccTwinTransNucleation .and. size(prm%N_tr) /= 1) extmsg = trim(extmsg)//' N_tr: nucleation' end if slipAndTransActive + twinAndTransActive: if (prm%sum_N_tw * prm%sum_N_tr > 0) then + if (dNeq(prm%b_tw(1),prm%b_tr(1))) extmsg = trim(extmsg)//' b_tw != b_tr' + end if twinAndTransActive + !-------------------------------------------------------------------------------------------------- ! allocate state arrays Nmembers = count(material_phaseID == ph)