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

change tag to all lowercase; make pheno hardeing only depend on its own system

This commit is contained in:
Zhuowen Zhao 2017-12-14 16:51:58 -05:00
parent dcf9e139d0
commit 14c0503a7e
2 changed files with 13 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/plastic_kinematichardening.f90
View File

@ -1,7 +1,7 @@
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH !> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @author Zhuowen Zhao, Michigan State University !> @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 !! formulation using a power law fitting
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
module plastic_kinehardening module plastic_kinehardening
@ -318,7 +318,7 @@ subroutine plastic_kinehardening_init(fileUnit)
do j = 1_pInt, Nchunks_SlipSlip do j = 1_pInt, Nchunks_SlipSlip
param(instance)%interaction_slipslip(j) = IO_floatValue(line,chunkPos,1_pInt+j) param(instance)%interaction_slipslip(j) = IO_floatValue(line,chunkPos,1_pInt+j)
enddo enddo
case ('nonSchmidCoeff') case ('nonschmidcoeff')
if (chunkPos(1) < 1_pInt + Nchunks_nonSchmid) & if (chunkPos(1) < 1_pInt + Nchunks_nonSchmid) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_KINEHARDENING_label//')') call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_KINEHARDENING_label//')')
do j = 1_pInt,Nchunks_nonSchmid do j = 1_pInt,Nchunks_nonSchmid
@ -350,10 +350,10 @@ subroutine plastic_kinehardening_init(fileUnit)
case ('n_slip') case ('n_slip')
param(instance)%n_slip = IO_floatValue(line,chunkPos,2_pInt) param(instance)%n_slip = IO_floatValue(line,chunkPos,2_pInt)
case ('aTolResistance') case ('atol_resistance')
param(instance)%aTolResistance = IO_floatValue(line,chunkPos,2_pInt) param(instance)%aTolResistance = IO_floatValue(line,chunkPos,2_pInt)
case ('aTolShear') case ('atol_shear')
param(instance)%aTolShear = IO_floatValue(line,chunkPos,2_pInt) param(instance)%aTolShear = IO_floatValue(line,chunkPos,2_pInt)
case default case default

16
src/plastic_phenopowerlaw.f90
View File

@ -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 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) slipSystems1: do i = 1_pInt,plastic_phenopowerlaw_Nslip(f,instance)
j = j+1_pInt j = j+1_pInt
left_SlipSlip(j) = 1.0_pReal + plastic_phenopowerlaw_H_int(f,instance) ! modified no system-dependent left part left_SlipSlip(j) = (1.0_pReal + plastic_phenopowerlaw_H_int(f,instance)) &
left_SlipTwin(j) = 1.0_pReal ! no system-dependent left part *abs(1.0_pReal-plasticState(ph)%state(j,of) / & ! 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_tausat_slip(f,instance)+ssat_offset)) & **plastic_phenopowerlaw_a_slip(instance)&
**plastic_phenopowerlaw_a_slip(instance)& *sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / &
*sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / & (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset))
(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 right_TwinSlip(j) = 1.0_pReal ! no system-dependent part
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------