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

added check to the "Plasticity_DetectChanges" test and removed commented-out parts for kinematic hardening

This commit is contained in:
Zhuowen Zhao 2018-05-03 22:39:49 -04:00
parent 362b958a8f
commit 76e6fa24fb
1 changed files with 3 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

48
src/plastic_kinematichardening.f90
View File

@ -51,12 +51,6 @@ module plastic_kinehardening
outputID !< ID of each post result output outputID !< ID of each post result output
real(pReal) :: & real(pReal) :: &
! F0, &
! mu, &
! mu0, &
! tau_hat0, &
! p1, &
! q1, &
gdot0, & !< reference shear strain rate for slip (input parameter) gdot0, & !< reference shear strain rate for slip (input parameter)
n_slip, & !< stress exponent for slip (input parameter) n_slip, & !< stress exponent for slip (input parameter)
aTolResistance, & aTolResistance, &
@ -325,25 +319,6 @@ subroutine plastic_kinehardening_init(fileUnit)
param(instance)%nonSchmidCoeff(j) = IO_floatValue(line,chunkPos,1_pInt+j) param(instance)%nonSchmidCoeff(j) = IO_floatValue(line,chunkPos,1_pInt+j)
enddo enddo
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! parameters independent of number of slip families
! case ('F0')
! param(instance)%F0 = IO_floatValue(line,chunkPos,2_pInt)
!
! case ('mu')
! param(instance)%mu = IO_floatValue(line,chunkPos,2_pInt)
!
! case ('mu0')
! param(instance)%mu0 = IO_floatValue(line,chunkPos,2_pInt)
!
! case ('tau_hat0')
! param(instance)%tau_hat0 = IO_floatValue(line,chunkPos,2_pInt)
!
! case ('p1')
! param(instance)%p1 = IO_floatValue(line,chunkPos,2_pInt)
!
! case ('q1')
! param(instance)%q1 = IO_floatValue(line,chunkPos,2_pInt)
case ('gdot0') case ('gdot0')
param(instance)%gdot0 = IO_floatValue(line,chunkPos,2_pInt) param(instance)%gdot0 = IO_floatValue(line,chunkPos,2_pInt)
@ -619,29 +594,6 @@ subroutine plastic_kinehardening_shearRates(gdot_pos,gdot_neg,tau_pos,tau_neg, &
state(instance)%crss(:,of))**param(instance)%n_slip & state(instance)%crss(:,of))**param(instance)%n_slip &
*sign(1.0_pReal,tau_neg-state(instance)%crss_back(:,of)) *sign(1.0_pReal,tau_neg-state(instance)%crss_back(:,of))
! gdot_pos = 0.5_pReal * param(instance)%gdot0 * &
! exp(-param(instance)%F0/(1.38e-23*298.15)* &
! (1-((abs(tau_pos-state(instance)%crss_back(:,of)) &
! -state(instance)%crss(:,of)*param(instance)%mu/param(instance)%mu) / &
! !----------------------------------------------------------------------------
! param(instance)%tau_hat0*param(instance)%mu/param(instance)%mu &
! )**param(instance)%p1 &
! )**param(instance)%q1 &
! )*sign(1.0_pReal,(tau_pos-state(instance)%crss_back(:,of)))
!
!
!
! gdot_neg = 0.5_pReal * param(instance)%gdot0 * &
! exp(-param(instance)%F0/(1.38e-23*298.15)* &
! (1-((abs(tau_neg-state(instance)%crss_back(:,of)) &
! -state(instance)%crss(:,of)*param(instance)%mu/param(instance)%mu) / &
! !----------------------------------------------------------------------------
! param(instance)%tau_hat0*param(instance)%mu/param(instance)%mu &
! )**param(instance)%p1 &
! )**param(instance)%q1 &
! )*sign(1.0_pReal,(tau_neg-state(instance)%crss_back(:,of)))
end subroutine plastic_kinehardening_shearRates end subroutine plastic_kinehardening_shearRates