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

added missing allocation of slopeLnRate array. 

changed parsing to lowercase name.
This commit is contained in:
Philip Eisenlohr 2012-04-20 12:18:38 +00:00
parent d2c4874db3
commit 83e22c00a2
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
code/constitutive_j2.f90
View File

@ -175,6 +175,8 @@ subroutine constitutive_j2_init(myFile)
constitutive_j2_n = 0.0_pReal constitutive_j2_n = 0.0_pReal
allocate(constitutive_j2_h0(maxNinstance)) allocate(constitutive_j2_h0(maxNinstance))
constitutive_j2_h0 = 0.0_pReal constitutive_j2_h0 = 0.0_pReal
allocate(constitutive_j2_h0_slopeLnRate(maxNinstance))
constitutive_j2_h0_slopeLnRate = 0.0_pReal
allocate(constitutive_j2_tausat(maxNinstance)) allocate(constitutive_j2_tausat(maxNinstance))
constitutive_j2_tausat = 0.0_pReal constitutive_j2_tausat = 0.0_pReal
allocate(constitutive_j2_a(maxNinstance)) allocate(constitutive_j2_a(maxNinstance))
@ -226,7 +228,7 @@ subroutine constitutive_j2_init(myFile)
constitutive_j2_n(i) = IO_floatValue(line,positions,2_pInt) constitutive_j2_n(i) = IO_floatValue(line,positions,2_pInt)
case ('h0') case ('h0')
constitutive_j2_h0(i) = IO_floatValue(line,positions,2_pInt) constitutive_j2_h0(i) = IO_floatValue(line,positions,2_pInt)
case ('h0_slope','slopeLnRate') case ('h0_slope','slopelnrate')
constitutive_j2_h0_slopeLnRate(i) = IO_floatValue(line,positions,2) constitutive_j2_h0_slopeLnRate(i) = IO_floatValue(line,positions,2)
case ('tausat') case ('tausat')
constitutive_j2_tausat(i) = IO_floatValue(line,positions,2_pInt) constitutive_j2_tausat(i) = IO_floatValue(line,positions,2_pInt)
@ -503,7 +505,7 @@ pure function constitutive_j2_dotState(Tstar_v, Temperature, state, g, ip, el)
! hardening coefficient ! hardening coefficient
if (abs(gamma_dot) > 1e-12_pReal) then if (abs(gamma_dot) > 1e-12_pReal) then
if (constitutive_j2_tausat_SinhFitA(matID) == 0) then if (constitutive_j2_tausat_SinhFitA(matID) == 0.0_pReal) then
saturation = constitutive_j2_tausat(matID) saturation = constitutive_j2_tausat(matID)
else else
saturation = ( constitutive_j2_tausat(matID) & saturation = ( constitutive_j2_tausat(matID) &