diff --git a/code/config/Phase_Phenopowerlaw_Aluminum.config b/code/config/Phase_Phenopowerlaw_Aluminum.config index 769e0cdfa..78972cfb4 100644 --- a/code/config/Phase_Phenopowerlaw_Aluminum.config +++ b/code/config/Phase_Phenopowerlaw_Aluminum.config @@ -36,7 +36,6 @@ twin_c 0 twin_d 0 twin_e 0 h0_slipslip 75e6 -h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 interaction_slipslip 1 1 1.4 1.4 1.4 1.4 diff --git a/code/config/Phase_Phenopowerlaw_BCC-Ferrite.config b/code/config/Phase_Phenopowerlaw_BCC-Ferrite.config index cb8447834..76bab500a 100644 --- a/code/config/Phase_Phenopowerlaw_BCC-Ferrite.config +++ b/code/config/Phase_Phenopowerlaw_BCC-Ferrite.config @@ -21,7 +21,6 @@ twin_c 0 twin_d 0 twin_e 0 h0_slipslip 495.0e6 -h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 interaction_slipslip 1 1 1.4 1.4 1.4 1.4 @@ -29,4 +28,4 @@ interaction_sliptwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twinslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twintwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 a_slip 1.0 -atol_resistance 1 \ No newline at end of file +atol_resistance 1 diff --git a/code/config/Phase_Phenopowerlaw_BCC-Martensite.config b/code/config/Phase_Phenopowerlaw_BCC-Martensite.config index 830761034..4de48ae91 100644 --- a/code/config/Phase_Phenopowerlaw_BCC-Martensite.config +++ b/code/config/Phase_Phenopowerlaw_BCC-Martensite.config @@ -20,7 +20,6 @@ twin_c 0 twin_d 0 twin_e 0 h0_slipslip 53500.0e6 # 35000.0e6 -h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 interaction_slipslip 1 1 1.4 1.4 1.4 1.4 @@ -28,4 +27,4 @@ interaction_sliptwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twinslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twintwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 a_slip 2.0 -atol_resistance 1 \ No newline at end of file +atol_resistance 1 diff --git a/code/config/Phase_Phenopowerlaw_Gold.config b/code/config/Phase_Phenopowerlaw_Gold.config index 90281fac8..7d9270c9f 100644 --- a/code/config/Phase_Phenopowerlaw_Gold.config +++ b/code/config/Phase_Phenopowerlaw_Gold.config @@ -43,7 +43,6 @@ twin_c 0 twin_d 0 twin_e 0 h0_slipslip 75e6 -h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 interaction_slipslip 1 1 1.4 1.4 1.4 1.4 diff --git a/code/config/Phase_Phenopowerlaw_cpTi-alpha.config b/code/config/Phase_Phenopowerlaw_cpTi-alpha.config index 1d0128ee5..5cd215fc0 100644 --- a/code/config/Phase_Phenopowerlaw_cpTi-alpha.config +++ b/code/config/Phase_Phenopowerlaw_cpTi-alpha.config @@ -38,7 +38,6 @@ twin_c 0 twin_d 0 twin_e 0 h0_slipslip 15e6 -h0_sliptwin 15e6 h0_twinslip 15e6 h0_twintwin 15e6 atol_resistance 1 @@ -49,4 +48,4 @@ interaction_slipslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_sliptwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twinslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twintwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -#nonschmid_coefficients 0 \ No newline at end of file +#nonschmid_coefficients 0 diff --git a/code/plastic_phenopowerlaw.f90 b/code/plastic_phenopowerlaw.f90 index 4d9bad3e2..d4b48a01c 100644 --- a/code/plastic_phenopowerlaw.f90 +++ b/code/plastic_phenopowerlaw.f90 @@ -46,7 +46,6 @@ module plastic_phenopowerlaw plastic_phenopowerlaw_twinD, & plastic_phenopowerlaw_twinE, & plastic_phenopowerlaw_h0_SlipSlip, & !< reference hardening slip - slip (input parameter) - plastic_phenopowerlaw_h0_SlipTwin, & !< reference hardening slip - twin (input parameter, no effect at the moment) plastic_phenopowerlaw_h0_TwinSlip, & !< reference hardening twin - slip (input parameter) plastic_phenopowerlaw_h0_TwinTwin, & !< reference hardening twin - twin (input parameter) plastic_phenopowerlaw_a_slip, & @@ -207,7 +206,6 @@ subroutine plastic_phenopowerlaw_init(fileUnit) allocate(plastic_phenopowerlaw_twinD(maxNinstance), source=0.0_pReal) allocate(plastic_phenopowerlaw_twinE(maxNinstance), source=0.0_pReal) allocate(plastic_phenopowerlaw_h0_SlipSlip(maxNinstance), source=0.0_pReal) - allocate(plastic_phenopowerlaw_h0_SlipTwin(maxNinstance), source=0.0_pReal) allocate(plastic_phenopowerlaw_h0_TwinSlip(maxNinstance), source=0.0_pReal) allocate(plastic_phenopowerlaw_h0_TwinTwin(maxNinstance), source=0.0_pReal) allocate(plastic_phenopowerlaw_interaction_SlipSlip(lattice_maxNinteraction,maxNinstance), & @@ -419,9 +417,6 @@ subroutine plastic_phenopowerlaw_init(fileUnit) plastic_phenopowerlaw_twinE(instance) = IO_floatValue(line,positions,2_pInt) case ('h0_slipslip') plastic_phenopowerlaw_h0_SlipSlip(instance) = IO_floatValue(line,positions,2_pInt) - case ('h0_sliptwin') - plastic_phenopowerlaw_h0_SlipTwin(instance) = IO_floatValue(line,positions,2_pInt) - call IO_warning(42_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_PHENOPOWERLAW_label//')') case ('h0_twinslip') plastic_phenopowerlaw_h0_TwinSlip(instance) = IO_floatValue(line,positions,2_pInt) case ('h0_twintwin') @@ -900,7 +895,7 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el) offset_accshear_slip,offset_accshear_twin, & of real(pReal) :: & - c_SlipSlip,c_SlipTwin,c_TwinSlip,c_TwinTwin, & + c_SlipSlip,c_TwinSlip,c_TwinTwin, & ssat_offset, & tau_slip_pos,tau_slip_neg,tau_twin @@ -928,7 +923,6 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el) c_SlipSlip = plastic_phenopowerlaw_h0_SlipSlip(instance)*& (1.0_pReal + plastic_phenopowerlaw_twinC(instance)*plasticState(ph)%state(index_F,of)**& plastic_phenopowerlaw_twinB(instance)) - c_SlipTwin = 0.0_pReal c_TwinSlip = plastic_phenopowerlaw_h0_TwinSlip(instance)*& plasticState(ph)%state(index_Gamma,of)**plastic_phenopowerlaw_twinE(instance) c_TwinTwin = plastic_phenopowerlaw_h0_TwinTwin(instance)*& @@ -999,7 +993,6 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el) c_SlipSlip * left_SlipSlip(j) * & dot_product(plastic_phenopowerlaw_hardeningMatrix_SlipSlip(j,1:nSlip,instance), & right_SlipSlip*abs(gdot_slip)) + & ! dot gamma_slip modulated by right-side slip factor - c_SlipTwin * left_SlipTwin(j) * & dot_product(plastic_phenopowerlaw_hardeningMatrix_SlipTwin(j,1:nTwin,instance), & right_SlipTwin*gdot_twin) ! dot gamma_twin modulated by right-side twin factor plasticState(ph)%dotState(index_Gamma,of) = plasticState(ph)%dotState(index_Gamma,of) + &