From 414050303b918d269b6407c5f8a61ebab1bccedb Mon Sep 17 00:00:00 2001 From: Denny Tjahjanto Date: Wed, 24 Mar 2010 13:28:46 +0000 Subject: [PATCH] forgotten: the change in the input config --- code/material.config | 37 +++++++++++++++++++++++-------------- code/numerics.config | 19 ++++++++++++------- 2 files changed, 35 insertions(+), 21 deletions(-) diff --git a/code/material.config b/code/material.config index 5cb068517..a475ad641 100644 --- a/code/material.config +++ b/code/material.config @@ -13,13 +13,22 @@ Ngrains 1 [Grain_Cluster] type RGC ngrains 8 -clustersize 2 2 2 # product of these numbers must be equal to ngrains(!) -grainsizeparameter 1.0e+3 1.0e+3 1.0e+3 # in [nm] -overproportionality 1.0e+1 1.0e+1 1.0e+1 # typical range between 0.1 (very large grain) to 100 (very small grain) +clustersize 2 2 2 # product of these numbers must be equal to ngrains(!) +clusterorientation 0.0 0.0 0.0 # orientation of cluster in terms of 313 Euler-angles in degree (random if not present) +# clusterorientation 0.0 26.57 0.0 # [012] +# clusterorientation 0.0 45.00 0.0 # [011] +# clusterorientation 0.0 26.57 24.10 # [112] +# clusterorientation 0.0 45.00 19.47 # [122] +# clusterorientation 0.0 45.00 35.26 # [111] +grainsize 4.0e-6 4.0e-6 2.0e-6 # in [m] +overproportionality 2.0e+0 # typical range between 0.001 to 1000 +scalingparameter 1.0e+1 # typical range between 0.001 to 1000 (output) constitutivework -(output) penaltyenergy (output) magnitudemismatch +(output) penaltyenergy (output) volumediscrepancy +(output) averagerelaxrate +(output) maximumrelaxrate [Taylor2] type isostrain @@ -235,7 +244,7 @@ c44 118.0e9 gdot0_slip 0.001 n_slip 20 tau0_slip 88.0e6 0 0 0 # per family -tausat_slip 201.0e6 0 0 0 # per family +tausat_slip 205.0e6 0 0 0 # per family gdot0_twin 0.001 n_twin 20 tau0_twin 31.0e6 0 0 0 # per family @@ -244,7 +253,7 @@ twin_b 0 twin_c 0 twin_d 0 twin_e 0 -h0_slipslip 391.0e6 +h0_slipslip 495.0e6 h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 @@ -260,13 +269,13 @@ constitution phenopowerlaw lattice_structure bcc Nslip 12 0 0 0 # per family Ntwin 0 0 0 0 # per family -c11 511.1e9 -c12 296.9e9 -c44 258.5e9 +c11 417.4e9 +c12 242.4e9 +c44 211.1e9 gdot0_slip 0.001 n_slip 20 -tau0_slip 396.0e6 0 0 0 # per family -tausat_slip 1120.0e6 0 0 0 # per family +tau0_slip 575.00e6 0 0 0 # 460.0e6 # per family +tausat_slip 1280.0e6 0 0 0 # 1020.0e6 # per family gdot0_twin 0.001 n_twin 20 tau0_twin 31.0e6 0 0 0 # per family @@ -275,7 +284,7 @@ twin_b 0 twin_c 0 twin_d 0 twin_e 0 -h0_slipslip 215000.0e6 +h0_slipslip 53500.0e6 # 35000.0e6 h0_sliptwin 0 h0_twinslip 0 h0_twintwin 0 @@ -283,7 +292,7 @@ interaction_slipslip 1 1 1.4 1.4 1.4 1.4 interaction_sliptwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twinslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 interaction_twintwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -w0_slip 4.0 +w0_slip 2.0 relevantResistance 1 @@ -362,4 +371,4 @@ symmetry orthotropic # or monoclinic [fiber example] # fiber axis in spherical coordinates: alpha crystal system, beta sample system -(fiber) alpha1 123 alpha2 123 beta1 12 beta2 45 scatter 15 fraction 0.333 \ No newline at end of file +(fiber) alpha1 123 alpha2 123 beta1 12 beta2 45 scatter 15 fraction 0.333 diff --git a/code/numerics.config b/code/numerics.config index 95a16bcc8..a0998aa76 100644 --- a/code/numerics.config +++ b/code/numerics.config @@ -31,12 +31,17 @@ rTol_RGC 1.0e-3 # relative ... aMax_RGC 1.0e+10 # absolute upper-limit of RGC residuum (in Pa) rMax_RGC 1.0e+3 # relative ... perturbPenalty_RGC 1.0e-7 # perturbation for computing penalty tangent -relevantMismatch_RGC 1.0e-5 # minimum threshold of mismatch -viscosityRate_RGC 1.0e+0 # power (sensitivity rate) of numerical viscosity in RGC scheme -viscosityModulus_RGC 0.0e+0 # stress modulus of RGC numerical viscosity (zero = without numerical viscosity) maxRelaxation_RGC 1.0e+0 # threshold of maximum relaxation vector increment (if exceed this then cutback) -maxvoldiscrepancy_RGC 1.0e-5 # threshold of maximum volume discrepancy allowed -voldiscrepancymod_RGC 0.0e+8 # energy modulus of RGC volume discrepancy (zero = without volume discrepancy constraint) -discrepancypower_RGC 5.0 # powerlaw penalty for volume discrepancy -fixed_seed 1234 # put any number larger than zero, integer, if you want to have a pseudo random distribution +relevantMismatch_RGC 1.0e-5 # minimum threshold of mismatch + +viscosityPower_RGC 1.0e+0 # power (sensitivity rate) of numerical viscosity in RGC scheme +viscosityModulus_RGC 1.0e+7 # stress modulus of RGC numerical viscosity (zero = without numerical viscosity) + # suggestion: larger than the aTol_RGC but still far below the expected flow stress of material +refRelaxationRate_RGC 1.0e-3 # reference rate of relaxation (about the same magnitude as straining rate, possibly a bit higher) + +maxVolDiscrepancy_RGC 1.0e-6 # maximum allowable relative volume discrepancy +volDiscrepancyMod_RGC 1.0e+5 +discrepancyPower_RGC 2.0 + +fixed_seed 10468 # put any number larger than zero, integer, if you want to have a pseudo random distribution