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

shortened

This commit is contained in:
Martin Diehl 2019-01-27 12:15:11 +01:00
parent 8c18b29420
commit aecb5f20bf
1 changed files with 10 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/plastic_dislotwin.f90
View File

@ -686,7 +686,7 @@ end subroutine plastic_dislotwin_init
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief returns the homogenized elasticity matrix !> @brief returns the homogenized elasticity matrix
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
function plastic_dislotwin_homogenizedC(ipc,ip,el) function plastic_dislotwin_homogenizedC(ipc,ip,el) result(homogenizedC)
use material, only: & use material, only: &
material_phase, & material_phase, &
phase_plasticityInstance, & phase_plasticityInstance, &
@ -694,7 +694,7 @@ function plastic_dislotwin_homogenizedC(ipc,ip,el)
implicit none implicit none
real(pReal), dimension(6,6) :: & real(pReal), dimension(6,6) :: &
plastic_dislotwin_homogenizedC homogenizedC
integer(pInt), intent(in) :: & integer(pInt), intent(in) :: &
ipc, & !< component-ID of integration point ipc, & !< component-ID of integration point
ip, & !< integration point ip, & !< integration point
@ -713,17 +713,18 @@ function plastic_dislotwin_homogenizedC(ipc,ip,el)
- sum(stt%stressTransFraction(1_pInt:prm%totalNtrans,of)) & - sum(stt%stressTransFraction(1_pInt:prm%totalNtrans,of)) &
- sum(stt%strainTransFraction(1_pInt:prm%totalNtrans,of)) - sum(stt%strainTransFraction(1_pInt:prm%totalNtrans,of))
plastic_dislotwin_homogenizedC = f_unrotated * prm%C66 homogenizedC = f_unrotated * prm%C66
do i=1_pInt,prm%totalNtwin do i=1_pInt,prm%totalNtwin
plastic_dislotwin_homogenizedC = plastic_dislotwin_homogenizedC & homogenizedC = homogenizedC &
+ stt%twinFraction(i,of)*prm%C66_twin(1:6,1:6,i) + stt%twinFraction(i,of)*prm%C66_twin(1:6,1:6,i)
enddo enddo
do i=1_pInt,prm%totalNtrans do i=1_pInt,prm%totalNtrans
plastic_dislotwin_homogenizedC = plastic_dislotwin_homogenizedC & homogenizedC = homogenizedC &
+(stt%stressTransFraction(i,of)+stt%strainTransFraction(i,of))*& +(stt%stressTransFraction(i,of)+stt%strainTransFraction(i,of))*prm%C66_trans(1:6,1:6,i)
prm%C66_trans(1:6,1:6,i)
enddo enddo
end associate end associate
end function plastic_dislotwin_homogenizedC end function plastic_dislotwin_homogenizedC