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

added kmc example material.config 

added elastic test for mkc, titanmod, and J2
This commit is contained in:
Martin Diehl 2014-08-08 13:51:11 +00:00
parent 0ecdd2692b
commit f55dc4b639
2 changed files with 45 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
code/config/Phase_Dislokmc_Tungsten.config Normal file
View File

@ -0,0 +1,43 @@
### $Id$ ###
[Tungsten]
elasticity hooke
plasticity dislokmc
### Material parameters ###
lattice_structure bcc
C11 523.0e9 # From Marinica et al. Journal of Physics: Condensed Matter(2013)
C12 202.0e9
C44 161.0e9
grainsize 2.0e-5 # Average grain size [m] 2.0e-5
SolidSolutionStrength 0.0 # Strength due to elements in solid solution
### Dislocation glide parameters ###
#per family
Nslip 12 0
slipburgers 2.72e-10 # Burgers vector of slip system [m]
rhoedge0 1.0e12 # Initial edge dislocation density [m/m**3] 1.0e12
rhoedgedip0 1.0 # Initial edged dipole dislocation density [m/m**3]
Qedge 2.725e-19 # Activation energy for dislocation glide [J]
v0 3560.3 # Initial glide velocity [m/s] 1.0e-4 (kmC)
p_slip 0.16 # p-exponent in glide velocity
q_slip 1.00 # q-exponent in glide velocity
u_slip 2.47 # u-exponent of stress pre-factor (kmC)
s_slip 0.97 # self hardening in glide velocity (kmC)
tau_peierls 2.03e9 # peierls stress [Pa]
#v_slip 1.58 # v-exponent in glide velocity (kmC)
#hardening
dipoleformationfactor 0 # to have hardening due to dipole formation off
CLambdaSlip 10.0 # Adj. parameter controlling dislocation mean free path
D0 4.0e-5 # Vacancy diffusion prefactor [m**2/s]
Qsd 4.5e-19 # Activation energy for climb [J]
Catomicvolume 1.0 # Adj. parameter controlling the atomic volume [in b]
Cedgedipmindistance 1.0 # Adj. parameter controlling the minimum dipole distance [in b]
#interaction_slipslip 0 0 0 0 0 0
#interaction_slipslip 0.1 0.1 0.14 0.14 0.14 0.14
#interaction_slipslip 1 0.1 0.9 0.5 0.1 0.7
interaction_slipslip 0.2 0.11 0.19 0.15 0.11 0.17
#interaction_slipslip 1 1 1.4 1.4 1.4 1.4

4
code/constitutive_titanmod.f90
View File

@ -669,7 +669,7 @@ subroutine constitutive_titanmod_init(fileUnit)
sanityChecks: do phase = 1_pInt, size(phase_plasticity)
myPhase: if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then
instance = phase_plasticityInstance(phase)
if (sum(constitutive_titanmod_Nslip(:,instance)) <= 0_pInt) &
if (sum(constitutive_titanmod_Nslip(:,instance)) < 0_pInt) &
call IO_error(211_pInt,el=instance,ext_msg='nslip ('//PLASTICITY_TITANMOD_label//')')
if (sum(constitutive_titanmod_Ntwin(:,instance)) < 0_pInt) &
call IO_error(211_pInt,el=instance,ext_msg='ntwin ('//PLASTICITY_TITANMOD_label//')')
@ -718,7 +718,7 @@ subroutine constitutive_titanmod_init(fileUnit)
if (constitutive_titanmod_twinhpconstant(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinhpconstant ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_aTolRho(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='aTolRho ('//PLASTICITY_TITANMOD_label//')')
call IO_error(211_pInt,el=instance,ext_msg='aTol_rho ('//PLASTICITY_TITANMOD_label//')')
!--------------------------------------------------------------------------------------------------
! determine total number of active slip or twin systems