adding damping mechanism

This commit is contained in:
zhangc43 2016-04-15 10:28:02 -04:00
parent 00abdc34c1
commit 1e406169ec
3 changed files with 54 additions and 52 deletions

View File

@ -480,7 +480,7 @@ subroutine constitutive_microstructure(orientations, Fe, Fp, ipc, ip, el, F0s,Fe
Fes, & Fes, &
Fps Fps
real(pReal), intent(in), dimension(6,:,:,:) :: & real(pReal), intent(in), dimension(:,:,:,:) :: &
Tstar_vs !< crystal orientations as quaternions Tstar_vs !< crystal orientations as quaternions
ho = material_homog(ip,el) ho = material_homog(ip,el)

View File

@ -437,10 +437,10 @@ subroutine crystallite_init
call constitutive_microstructure(crystallite_orientation, & ! pass orientation to constitutive module call constitutive_microstructure(crystallite_orientation, & ! pass orientation to constitutive module
crystallite_Fe(1:3,1:3,c,i,e), & crystallite_Fe(1:3,1:3,c,i,e), &
crystallite_Fp(1:3,1:3,c,i,e), & crystallite_Fp(1:3,1:3,c,i,e), &
c,i,e, c,i,e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo enddo
enddo enddo
@ -1718,10 +1718,10 @@ subroutine crystallite_integrateStateRK4()
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP ENDDO
@ -2048,10 +2048,10 @@ subroutine crystallite_integrateStateRKCK45()
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP ENDDO
@ -2272,10 +2272,10 @@ subroutine crystallite_integrateStateRKCK45()
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP ENDDO
@ -2511,10 +2511,10 @@ subroutine crystallite_integrateStateAdaptiveEuler()
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP ENDDO
@ -2859,10 +2859,10 @@ eIter = FEsolving_execElem(1:2)
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP ENDDO
@ -3108,10 +3108,10 @@ subroutine crystallite_integrateStateFPI()
call constitutive_microstructure(crystallite_orientation, & call constitutive_microstructure(crystallite_orientation, &
crystallite_Fe(1:3,1:3,g,i,e), & crystallite_Fe(1:3,1:3,g,i,e), &
crystallite_Fp(1:3,1:3,g,i,e), & crystallite_Fp(1:3,1:3,g,i,e), &
g, i, e, g, i, e, &
crystallite_F0, crystallite_F0, &
crystallite_Fe, crystallite_Fe, &
crystallite_Fp, crystallite_Fp, &
crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states crystallite_Tstar_v) ! update dependent state variables to be consistent with basic states
p = phaseAt(g,i,e) p = phaseAt(g,i,e)
c = phasememberAt(g,i,e) c = phasememberAt(g,i,e)

View File

@ -793,6 +793,7 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
ph, & !my phase ph, & !my phase
of, & !my spatial position in memory (offset) of, & !my spatial position in memory (offset)
textureID, & !my texture textureID, & !my texture
index_myFamily, &
Nneighbors, & !number of neighbors (<= 6) Nneighbors, & !number of neighbors (<= 6)
vld_Nneighbors, & !number of my valid neighbors vld_Nneighbors, & !number of my valid neighbors
n, & !neighbor index (for iterating through all neighbors) n, & !neighbor index (for iterating through all neighbors)
@ -814,21 +815,21 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
f,i,& !loop counter for me f,i,& !loop counter for me
f_ne, i_ne !loop counter for neighbor f_ne, i_ne !loop counter for neighbor
real(pReal) kappa_max, & ! real(pReal) mprime_cut, & !m' cutoff to consider neighboring effect
tmp_myshear_slip, & !temp storage for accumulative shear for me
mprime_cut, & !m' cutoff to consider neighboring effect
dtaylor_cut, & !threshold for determine high contrast interface using Taylor factor dtaylor_cut, & !threshold for determine high contrast interface using Taylor factor
avg_acshear_ne, & !the average accumulative shear from my neighbor tau_slip, & !the average accumulative shear from my neighbor
taylor_me, & !Taylor factor for me taylor_me, & !Taylor factor for me
taylor_ne, & !Taylor factor for my current neighbor taylor_ne, & !Taylor factor for my current neighbor
d_vonstrain !von Mises delta strain (temp container) d_vonstrain, & !von Mises delta strain (temp container)
sum_gdot !total shear rate for given material point
real(pReal), dimension(3,3) :: & real(pReal), dimension(3,3) :: &
F0_me, & !my deformation gradient from last converged increment F0_me, & !my deformation gradient from last converged increment
Fe_me, & !my elastic deformation gradient Fe_me, & !my elastic deformation gradient
Fp_me, & !my plastic deformation gradient Fp_me, & !my plastic deformation gradient
dF_me, & !my deformation gradient change (delta) dF_me, & !my deformation gradient change (delta)
dE_me, & !my Green Lagrangian strain tensor (delta) dE_me, & !my Green Lagrangian strain tensor (delta)
F0_ne, & !
Fe_ne, & !elastic deformation gradient of my current neighbor Fe_ne, & !elastic deformation gradient of my current neighbor
Fp_ne, & !plastic deformation gradient of my current neighbor Fp_ne, & !plastic deformation gradient of my current neighbor
dF_ne, & !deformation gradient of my current neighbor dF_ne, & !deformation gradient of my current neighbor
@ -874,7 +875,7 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
!***check if all my neighbors have the same phase as me !***check if all my neighbors have the same phase as me
vld_Nneighbors = 0 vld_Nneighbors = 0
PHASECHECK DO n_phasecheck = 1_pInt, Nneighbors PHASECHECK: DO n_phasecheck = 1_pInt, Nneighbors
!******for each of my neighbor !******for each of my neighbor
neighbor_el = mesh_ipNeighborhood( 1, n_phasecheck, ip, el ) neighbor_el = mesh_ipNeighborhood( 1, n_phasecheck, ip, el )
neighbor_ip = mesh_ipNeighborhood( 2, n_phasecheck, ip, el ) neighbor_ip = mesh_ipNeighborhood( 2, n_phasecheck, ip, el )
@ -905,7 +906,7 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
tau_slip = dot_product(Tstar_v(1:6, ipc, ip, el),lattice_Sslip_v(1:6,1,index_myFamily+i,ph)) tau_slip = dot_product(Tstar_v(1:6, ipc, ip, el),lattice_Sslip_v(1:6,1,index_myFamily+i,ph))
sum_gdot = sum_gdot + & sum_gdot = sum_gdot + &
plastic_phenoplus_gdot0_slip(instance)* & plastic_phenoplus_gdot0_slip(instance)* &
((abs(tau_slip)/(state(instance)%s_slip(j,of))) & ((abs(tau_slip)/(plasticState(ph)%state(j,of))) &
**plastic_phenoplus_n_slip(instance))*sign(1.0_pReal,tau_slip) **plastic_phenoplus_n_slip(instance))*sign(1.0_pReal,tau_slip)
ENDDO slipSystems ENDDO slipSystems
ENDDO slipFamilies ENDDO slipFamilies
@ -934,14 +935,14 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
j = 0_pInt j = 0_pInt
slipFamiliesNeighbor: DO f_ne = 1_pInt,lattice_maxNslipFamily slipFamiliesNeighbor: DO f_ne = 1_pInt,lattice_maxNslipFamily
index_myFamily = sum(lattice_NslipSystem(1:f_ne-1_pInt,neighbor_ph)) ! at which index starts my family index_myFamily = sum(lattice_NslipSystem(1:f_ne-1_pInt,neighbor_ph)) ! at which index starts my family
slipSystemsNeighbor: DO i_ne = 1_pInt,plastic_phenopowerlaw_Nslip(f_ne,neighbor_instance) slipSystemsNeighbor: DO i_ne = 1_pInt,plastic_phenoplus_Nslip(f_ne,neighbor_instance)
j = j+1_pInt j = j+1_pInt
tau_slip = dot_product(Tstar_v(1:6, neighbor_ipc, neighbor_ip, neighbor_el), tau_slip = dot_product(Tstar_v(1:6, neighbor_ipc, neighbor_ip, neighbor_el), &
lattice_Sslip_v(1:6,1,index_myFamily+i_ne,neighbor_ph)) lattice_Sslip_v(1:6,1,index_myFamily+i_ne,neighbor_ph))
sum_gdot = sum_gdot & sum_gdot = sum_gdot &
+plastic_phenopowerlaw_gdot0_slip(neighbor_instance) & +plastic_phenoplus_gdot0_slip(neighbor_instance) &
*((abs(tau_slip)/(state(neighbor_instance)%s_slip(j,neighbor_of))) & *((abs(tau_slip)/(plasticState(neighbor_ph)%state(j,neighbor_of))) &
**plastic_phenopowerlaw_n_slip(neighbor_instance))*sign(1.0_pReal,tau_slip) **plastic_phenoplus_n_slip(neighbor_instance))*sign(1.0_pReal,tau_slip)
ENDDO slipSystemsNeighbor ENDDO slipSystemsNeighbor
ENDDO slipFamiliesNeighbor ENDDO slipFamiliesNeighbor
taylor_ne = d_vonstrain / sum_gdot taylor_ne = d_vonstrain / sum_gdot
@ -950,11 +951,11 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
ENDDO LOOPCALCTAYLOR ENDDO LOOPCALCTAYLOR
!***Only perform necessary calculation if high contrast interface is detected !***Only perform necessary calculation if high contrast interface is detected
IF (max(d_taylors) > dtaylor_cut) THEN IF (maxval(d_taylors) > dtaylor_cut) THEN
!*****calculate kappa per slip system base !*****calculate kappa per slip system base
LOOPMYSLIP DO me_slip = 1_pInt, ns LOOPMYSLIP: DO me_slip = 1_pInt, ns
ne_mprimes = 0.0_pReal !initialize max m' to 0 for all neighbors ne_mprimes = 0.0_pReal !initialize max m' to 0 for all neighbors
LOOPMYNEIGHBORS DO n=1_pInt, Nneighbors LOOPMYNEIGHBORS: DO n=1_pInt, Nneighbors
!*******only consider neighbor at the high contrast interface !*******only consider neighbor at the high contrast interface
IF (d_taylors(n) > dtaylor_cut) THEN IF (d_taylors(n) > dtaylor_cut) THEN
neighbor_el = mesh_ipNeighborhood( 1, n_calcTaylor, ip, el ) neighbor_el = mesh_ipNeighborhood( 1, n_calcTaylor, ip, el )
@ -975,14 +976,15 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
*abs(math_mul3x3(slipDirect(1:3,me_slip), & *abs(math_mul3x3(slipDirect(1:3,me_slip), &
math_qRot(absMisorientation, slipDirect(1:3,ne_slip)))) math_qRot(absMisorientation, slipDirect(1:3,ne_slip))))
ENDDO LOOPNEIGHBORSLIP ENDDO LOOPNEIGHBORSLIP
ne_mprimes(n) = max(m_primes) ne_mprimes(n) = maxval(m_primes)
ENDIF
!*******check if one of the neighbor already can provide a kick for this slip system
IF ( max(ne_mprimes) > mprime_cut ) THEN
plasticState(ph)%state(index_kappa+me_slip, of) = 1.5_pReal
EXIT
ENDIF ENDIF
ENDDO LOOPMYNEIGHBORS ENDDO LOOPMYNEIGHBORS
!*******check if one of the neighbor already can provide a kick for this slip system
IF ( maxval(ne_mprimes) > mprime_cut ) THEN
plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal + 0.2_pReal*maxval(ne_mprimes)
ELSE
plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal - 0.1_pReal*maxval(ne_mprimes)
ENDIF
ENDDO LOOPMYSLIP ENDDO LOOPMYSLIP
ENDIF ENDIF