223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

remove outdated code in comment

This commit is contained in:
zhangc43 2016-04-21 14:07:20 -04:00
parent 45872e840f
commit 4029346aed
1 changed files with 0 additions and 71 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

71
src/plastic_phenoplus.f90
View File

@ -990,77 +990,6 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
ENDIF ENDIF
! !***calculate kappa between me and all my neighbors
! LOOPMYSLIP: DO me_slip=1_pInt,ns
! vld_Nneighbors = Nneighbors
! tmp_myshear_slip = me_acshear(me_slip)
! tmp_mprime = 0.0_pReal !highest m' from all neighbors
! tmp_acshear = 0.0_pReal !accumulative shear from highest m'
! !***go through my neighbors to find highest m'
! LOOPNEIGHBORS: DO n=1_pInt,Nneighbors
! neighbor_el = mesh_ipNeighborhood(1,n,ip,el)
! neighbor_ip = mesh_ipNeighborhood(2,n,ip,el)
! neighbor_n = 1 !It is ipc
! neighbor_of = phasememberAt( neighbor_n, neighbor_ip, neighbor_el)
! neighbor_ph = phaseAt( neighbor_n, neighbor_ip, neighbor_el)
! neighbor_tex = material_texture(1,neighbor_ip,neighbor_el)
! neighbor_orientation = orientation(1:4, neighbor_n, neighbor_ip, neighbor_el) !ipc is always 1.
! absMisorientation = lattice_qDisorientation(my_orientation, &
! neighbor_orientation, &
! 0_pInt) !no need for explicit calculation of symmetry
! !***find the accumulative shear for this neighbor
! LOOPFINDNEISHEAR: DO ne_slip_ac=1_pInt, ns
! ne_acshear(ne_slip_ac) = plasticState(ph)%state(offset_acshear_slip+ne_slip_ac, &
! neighbor_of)
! ENDDO LOOPFINDNEISHEAR
! !***calculate the average accumulative shear and use it as cutoff
! avg_acshear_ne = SUM(ne_acshear)/ns
! !***
! IF (ph==neighbor_ph) THEN
! !***walk through all the
! LOOPNEIGHBORSLIP: DO ne_slip=1_pInt,ns
! !***only consider slip system that is active (above average accumulative shear)
! IF (ne_acshear(ne_slip) > avg_acshear_ne) THEN
! m_primes(ne_slip) = abs(math_mul3x3(slipNormal(1:3,me_slip), &
! math_qRot(absMisorientation, slipNormal(1:3,ne_slip)))) &
! *abs(math_mul3x3(slipDirect(1:3,me_slip), &
! math_qRot(absMisorientation, slipDirect(1:3,ne_slip))))
! !***find the highest m' and corresponding accumulative shear
! IF (m_primes(ne_slip) > tmp_mprime) THEN
! tmp_mprime = m_primes(ne_slip)
! tmp_acshear = ne_acshear(ne_slip)
! ENDIF
! ENDIF
! ENDDO LOOPNEIGHBORSLIP
! ELSE
! ne_mprimes(n) = 0.0_pReal
! vld_Nneighbors = vld_Nneighbors - 1_pInt
! ENDIF
! ENDDO LOOPNEIGHBORS
! !***check if this element close to rim
! IF (vld_Nneighbors < Nneighbors) THEN
! !***rim voxel, no modification allowed
! plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal
! ELSE
! !***patch voxel, started to calculate push up factor for gamma_dot
! IF ((tmp_mprime > mprime_cut) .AND. (tmp_acshear > tmp_myshear_slip)) THEN
! plasticState(ph)%state(index_kappa+me_slip, of) = 1.5_pReal
! ELSE
! !***minimum damping factor is 0.5
! plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal
! ENDIF
! ENDIF
! ENDDO LOOPMYSLIP
end subroutine plastic_phenoplus_microstructure end subroutine plastic_phenoplus_microstructure