drag coefficient B introduced and read in
This commit is contained in:
parent
e7b8adadd0
commit
5d1e648d17
|
@ -109,7 +109,8 @@ module plastic_dislotwin
|
||||||
q, & !< q-exponent in glide velocity
|
q, & !< q-exponent in glide velocity
|
||||||
r, & !< r-exponent in twin nucleation rate
|
r, & !< r-exponent in twin nucleation rate
|
||||||
s, & !< s-exponent in trans nucleation rate
|
s, & !< s-exponent in trans nucleation rate
|
||||||
shear_twin !< characteristic shear for twins
|
shear_twin, & !< characteristic shear for twins
|
||||||
|
B !< drag coefficient
|
||||||
real(pReal), dimension(:,:), allocatable, private :: &
|
real(pReal), dimension(:,:), allocatable, private :: &
|
||||||
interaction_SlipSlip, & !< coefficients for slip-slip interaction for each interaction type and instance
|
interaction_SlipSlip, & !< coefficients for slip-slip interaction for each interaction type and instance
|
||||||
interaction_SlipTwin, & !< coefficients for slip-twin interaction for each interaction type and instance
|
interaction_SlipTwin, & !< coefficients for slip-twin interaction for each interaction type and instance
|
||||||
|
@ -326,6 +327,9 @@ subroutine plastic_dislotwin_init(fileUnit)
|
||||||
if (size(prm%burgers_slip) /= size(prm%Nslip)) call IO_error(150_pInt,ext_msg='slipburgers')
|
if (size(prm%burgers_slip) /= size(prm%Nslip)) call IO_error(150_pInt,ext_msg='slipburgers')
|
||||||
prm%burgers_slip = math_expand(prm%burgers_slip,prm%Nslip)
|
prm%burgers_slip = math_expand(prm%burgers_slip,prm%Nslip)
|
||||||
|
|
||||||
|
prm%B = config_phase(p)%getFloats('B',defaultVal=[(0.0_pReal, i=1,size(prm%Nslip))])
|
||||||
|
prm%B = math_expand(prm%B,prm%Nslip)
|
||||||
|
|
||||||
prm%Qedge = config_phase(p)%getFloats('qedge')
|
prm%Qedge = config_phase(p)%getFloats('qedge')
|
||||||
prm%Qedge = math_expand(prm%Qedge,prm%Nslip)
|
prm%Qedge = math_expand(prm%Qedge,prm%Nslip)
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue