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

Changed to SI conformity. Mattex file changed bg->burgers

This commit is contained in:
Luc Hantcherli 2007-11-15 17:04:30 +00:00
parent ad3ccf22c2
commit 3c0d51542e
2 changed files with 23 additions and 25 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
trunk/constitutive_dislo.f90
View File

@ -22,8 +22,8 @@ character(len=300), parameter :: mattexFile = 'mattex.mpie'
!*************************************
!* Physical parameter, attack_frequency != Debye frequency
real(pReal), parameter :: attack_frequency = 1.0e10_pReal
!* Physical parameter, Boltzman constant in mJ/Kelvin
real(pReal), parameter :: Kb = 1.38e-20_pReal
!* Physical parameter, Boltzman constant in J/Kelvin
real(pReal), parameter :: Kb = 1.38e-23_pReal
!*************************************
!* Definition of material properties *
@ -301,16 +301,16 @@ do while(.true.)
material_C33(section)=IO_floatValue(line,positions,2)
case ('c44')
material_C44(section)=IO_floatValue(line,positions,2)
case ('rho0') !* conversion in 1/mm²
material_rho0(section)=IO_floatValue(line,positions,2)/1.0e6_pReal
case ('rho0')
material_rho0(section)=IO_floatValue(line,positions,2)
case ('interaction_coefficients')
do i=1,6
material_SlipIntCoeff(i,section)=IO_floatValue(line,positions,i+1)
enddo
case ('bg') !* conversion in mm
material_bg(section)=IO_floatValue(line,positions,2)*1.0e3_pReal
case ('Qedge') !* conversion in mJ/Kelvin
material_Qedge(section)=IO_floatValue(line,positions,2)*1.0e3_pReal
case ('burgers')
material_bg(section)=IO_floatValue(line,positions,2)
case ('Qedge')
material_Qedge(section)=IO_floatValue(line,positions,2)
case ('tau0')
material_tau0(section)=IO_floatValue(line,positions,2)
case ('c1')
@ -886,14 +886,10 @@ matID = constitutive_matID(ipc,ip,el)
do i=1,constitutive_Nstatevars(ipc,ip,el)
tau_slip = dot_product(Tstar_v,crystal_Sslip_v(:,i,material_CrystalStructure(matID)))
gdot_slip = constitutive_g0_slip(i)*sinh((abs(tau_slip)*constitutive_activation_volume(i))/(Kb*Tp))*&
sign(1.0_pReal,tau_slip)
if (abs(tau_slip)>1.0e-20_pReal) then
lock(i)=(material_c4(matID)*sqrt(constitutive_rho_f(i))*abs(gdot_slip))/material_bg(matID)
recovery(i)=material_c5(matID)*state(i)*abs(gdot_slip)
constitutive_dotState(i)=lock(i)-recovery(i)
else
constitutive_dotState(i)=0.0_pReal
endif
sign(1.0_pReal,tau_slip)
lock(i)=(material_c4(matID)*sqrt(constitutive_rho_f(i))*abs(gdot_slip))/material_bg(matID)
recovery(i)=material_c5(matID)*state(i)*abs(gdot_slip)
constitutive_dotState(i)=lock(i)-recovery(i)
enddo
return

20
trunk/mattex.mpie
View File

@ -2,17 +2,19 @@
[Aluminium]
crystal_structure 1
Nslip 12
C11 106.75e3
C12 60.41e3
C44 28.34e3
## Elastic constants
# Unit in [Pa]
C11 106.75e9
C12 60.41e9
C44 28.34e9
## Parameters for phenomenological modeling
# Unit in [MPa]
s0_slip 31.0
# Unit in [Pa]
s0_slip 31.0e6
gdot0_slip 0.001
n_slip 20
h0 75
s_sat 63
h0 75.0e6
s_sat 63.0e6
w0 2.25
# Self and latent hardening coefficients
hardening_coefficients 1.0 1.4
@ -21,10 +23,10 @@ hardening_coefficients 1.0 1.4
# Initial dislocation density [m]²
rho0 1.5e11
# Burgers vector [m]
bg 2.86e-10
burgers 2.86e-10
# Activation energy for dislocation glide [J/K]
Qedge 3.0e-19
# Reference for passing stress [MPa]
# Reference for passing stress [Pa]
tau0 0.0
# Passing stress adjustment
c1 0.1