DAMASK_EICMD/examples/ConfigFiles/Phase_Phenopowerlaw_Magnesi...

#-------------------#
<phase>
#-------------------#
/echo/
[Mg]
plasticity            phenopowerlaw
elasticity            hooke

(output)                resistance_slip
(output)                shearrate_slip
(output)                resolvedstress_slip
(output)                resistance_twin
(output)                shearrate_twin
(output)                resolvedstress_twin

lattice_structure       hex
c/a                     1.62350                          # from Tromans 2011, Elastic Anisotropy of HCP Metal Crystals and Polycrystals
c11                     59.3e9                           #  - " -
c33                     61.5e9                           #  - " -
c44                     16.4e9                           #  - " -
c12                     25.7e9                           #  - " -
c13                     21.4e9                           #  - " -

#                      basal  prism  prism pyr(a) pyr(c+a) pyr(c+a)
Nslip                  3      3      0      6        0      6      # from Agnew et al 2006, Validating a polycrystal model for the elastoplastic response of mg alloy AZ32 using in situ neutron diffraction
#                      T1 C1 T2 C2         
Ntwin                  6  0  0  6                                  #  - " - 
#                      basal  prism  prism pyr(a) pyr(c+a) pyr(c+a)
tau0_slip               10.0e6  55.0e6  0   60.0e6  0.0   60.0e6     #  - " - table 1, pyr(a) set to pyr(c+a)
tausat_slip             40.0e6 135.0e6  0  150.0e6  0.0  150.0e6     #  - " - table 1, pyr(a) set to pyr(c+a)
#                      T1    C1   T2   C2 
tau0_twin              40e6  0.0  0.0  60.0e6                        #  - " - table 1, compressive twin guessed by Steffi, tensile twin modified to match experimental results

h0_twintwin           50.0e6                                         #  - " - table 1, same range as theta_0
h0_slipslip          500.0e6                                         #  - " - table 1, same range as theta_0
h0_twinslip          150.0e6                                         # guessing

interaction_slipslip    1 1 1 1 1 1 1 1 1 1   1 1 1 1 1 1 1 1 1 1    1 1 1 1 1 1 1 1 1 1    1 1 1 1 1 1 1 1 1 1   1 1  # just guessing
interaction_twintwin    1 1 1 1 1 1 1 1 1 1   1 1 1 1 1 1 1 1 1 1                                                      #  - " - 
interaction_sliptwin    1 1 1 1 1 1 1 1 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 1 1 1 1    1 1 1 1                                           #  - " - 


####################################################
# open for discussion
####################################################
atol_resistance         1
atol_shear              1e-6
atol_twinfrac           1e-6
n_twin                  20
n_slip                  20

gdot0_twin              0.001
gdot0_slip              0.001

twin_b                  0
twin_c                  0
twin_d                  20
twin_e                  20

a_slip                  2.25
s_pr                    10.0                  # push-up factor for slip saturation due to twinning
added parameter set for Magensium 2014-01-15 17:59:44 +05:30			`#-------------------#`
			`<phase>`
			`#-------------------#`
			`/echo/`
			`[Mg]`
			`plasticity phenopowerlaw`
			`elasticity hooke`

			`(output) resistance_slip`
			`(output) shearrate_slip`
			`(output) resolvedstress_slip`
			`(output) resistance_twin`
			`(output) shearrate_twin`
			`(output) resolvedstress_twin`

			`lattice_structure hex`
c/a ratio matches the lattice_init name 2019-03-29 17:00:00 +05:30			`c/a 1.62350 # from Tromans 2011, Elastic Anisotropy of HCP Metal Crystals and Polycrystals`
added parameter set for Magensium 2014-01-15 17:59:44 +05:30			`c11 59.3e9 # - " -`
			`c33 61.5e9 # - " -`
			`c44 16.4e9 # - " -`
			`c12 25.7e9 # - " -`
			`c13 21.4e9 # - " -`

			`# basal prism prism pyr(a) pyr(c+a) pyr(c+a)`
			`Nslip 3 3 0 6 0 6 # from Agnew et al 2006, Validating a polycrystal model for the elastoplastic response of mg alloy AZ32 using in situ neutron diffraction`
			`# T1 C1 T2 C2`
			`Ntwin 6 0 0 6 # - " -`
			`# basal prism prism pyr(a) pyr(c+a) pyr(c+a)`
			`tau0_slip 10.0e6 55.0e6 0 60.0e6 0.0 60.0e6 # - " - table 1, pyr(a) set to pyr(c+a)`
			`tausat_slip 40.0e6 135.0e6 0 150.0e6 0.0 150.0e6 # - " - table 1, pyr(a) set to pyr(c+a)`
			`# T1 C1 T2 C2`
			`tau0_twin 40e6 0.0 0.0 60.0e6 # - " - table 1, compressive twin guessed by Steffi, tensile twin modified to match experimental results`

			`h0_twintwin 50.0e6 # - " - table 1, same range as theta_0`
			`h0_slipslip 500.0e6 # - " - table 1, same range as theta_0`
			`h0_twinslip 150.0e6 # guessing`

			`interaction_slipslip 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 # just guessing`
			`interaction_twintwin 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 # - " -`
			`interaction_sliptwin 1 1 1 1 1 1 1 1 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 1 1 1 1 1 1 1 1 # - " -`



			`####################################################`
			`# open for discussion`
			`####################################################`
			`atol_resistance 1`
atol_shear and atol_twinfrac are not allowed to be zero 2019-03-29 17:06:30 +05:30			`atol_shear 1e-6`
			`atol_twinfrac 1e-6`
added parameter set for Magensium 2014-01-15 17:59:44 +05:30			`n_twin 20`
			`n_slip 20`

			`gdot0_twin 0.001`
			`gdot0_slip 0.001`

			`twin_b 0`
			`twin_c 0`
			`twin_d 20`
			`twin_e 20`

			`a_slip 2.25`
			`s_pr 10.0 # push-up factor for slip saturation due to twinning`