From 3c0d51542ead8686a5d86761aa7e763a92c9047b Mon Sep 17 00:00:00 2001 From: Luc Hantcherli Date: Thu, 15 Nov 2007 17:04:30 +0000 Subject: [PATCH] Changed to SI conformity. Mattex file changed bg->burgers --- trunk/constitutive_dislo.f90 | 28 ++++++++++++---------------- trunk/mattex.mpie | 20 +++++++++++--------- 2 files changed, 23 insertions(+), 25 deletions(-) diff --git a/trunk/constitutive_dislo.f90 b/trunk/constitutive_dislo.f90 index fa4fc3ce2..6a62a3e59 100644 --- a/trunk/constitutive_dislo.f90 +++ b/trunk/constitutive_dislo.f90 @@ -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 diff --git a/trunk/mattex.mpie b/trunk/mattex.mpie index 5ec920d9a..5ffd56fe6 100644 --- a/trunk/mattex.mpie +++ b/trunk/mattex.mpie @@ -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