From 14c0503a7e1ffc0141553cc8728edcce97808b40 Mon Sep 17 00:00:00 2001 From: Zhuowen Zhao Date: Thu, 14 Dec 2017 16:51:58 -0500 Subject: [PATCH] change tag to all lowercase; make pheno hardeing only depend on its own system --- src/plastic_kinematichardening.f90 | 8 ++++---- src/plastic_phenopowerlaw.f90 | 16 +++++++++------- 2 files changed, 13 insertions(+), 11 deletions(-) diff --git a/src/plastic_kinematichardening.f90 b/src/plastic_kinematichardening.f90 index a82d9066a..3bd35aaea 100644 --- a/src/plastic_kinematichardening.f90 +++ b/src/plastic_kinematichardening.f90 @@ -1,7 +1,7 @@ !-------------------------------------------------------------------------------------------------- !> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH !> @author Zhuowen Zhao, Michigan State University -!> @brief Introducing Voce-type kinematic hardening rule into crystal phenopowerlaw plasticity +!> @brief Introducing Voce-type kinematic hardening rule into crystal plasticity !! formulation using a power law fitting !-------------------------------------------------------------------------------------------------- module plastic_kinehardening @@ -318,7 +318,7 @@ subroutine plastic_kinehardening_init(fileUnit) do j = 1_pInt, Nchunks_SlipSlip param(instance)%interaction_slipslip(j) = IO_floatValue(line,chunkPos,1_pInt+j) enddo - case ('nonSchmidCoeff') + case ('nonschmidcoeff') if (chunkPos(1) < 1_pInt + Nchunks_nonSchmid) & call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_KINEHARDENING_label//')') do j = 1_pInt,Nchunks_nonSchmid @@ -350,10 +350,10 @@ subroutine plastic_kinehardening_init(fileUnit) case ('n_slip') param(instance)%n_slip = IO_floatValue(line,chunkPos,2_pInt) - case ('aTolResistance') + case ('atol_resistance') param(instance)%aTolResistance = IO_floatValue(line,chunkPos,2_pInt) - case ('aTolShear') + case ('atol_shear') param(instance)%aTolShear = IO_floatValue(line,chunkPos,2_pInt) case default diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90 index 319f1b585..facc76ff1 100644 --- a/src/plastic_phenopowerlaw.f90 +++ b/src/plastic_phenopowerlaw.f90 @@ -987,13 +987,15 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el) index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family slipSystems1: do i = 1_pInt,plastic_phenopowerlaw_Nslip(f,instance) j = j+1_pInt - left_SlipSlip(j) = 1.0_pReal + plastic_phenopowerlaw_H_int(f,instance) ! modified no system-dependent left part - left_SlipTwin(j) = 1.0_pReal ! no system-dependent left part - right_SlipSlip(j) = abs(1.0_pReal-plasticState(ph)%state(j,of) / & - (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) & - **plastic_phenopowerlaw_a_slip(instance)& - *sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / & - (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) + left_SlipSlip(j) = (1.0_pReal + plastic_phenopowerlaw_H_int(f,instance)) & + *abs(1.0_pReal-plasticState(ph)%state(j,of) / & ! no system-dependent left part + (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) & + **plastic_phenopowerlaw_a_slip(instance)& + *sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / & + (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) + left_SlipTwin(j) = 1.0_pReal + right_SlipSlip(j) = 1.0_pReal ! system-dependent part (beta summation) + right_TwinSlip(j) = 1.0_pReal ! no system-dependent part !--------------------------------------------------------------------------------------------------