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

correct names 

whether C66 is homogenized or not is a decision of the caller
This commit is contained in:
Martin Diehl 2021-11-18 21:31:08 +01:00
parent dfe6d0a195
commit 038dd1fc40
2 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/phase.f90
View File

@ -234,10 +234,10 @@ module phase
integer, intent(in) :: ph, en
real(pReal), dimension(6,6) :: C
end function phase_homogenizedC66
module function phase_damage_C66(C66_homogenized,ph,en) result(C66)
real(pReal), dimension(6,6), intent(in) :: C66_homogenized
module function phase_damage_C66(C66,ph,en)
real(pReal), dimension(6,6), intent(in) :: C66
integer, intent(in) :: ph,en
real(pReal), dimension(6,6) :: C66
real(pReal), dimension(6,6) :: phase_damage_C66
end function phase_damage_C66
module function phase_f_phi(phi,co,ce) result(f)

10
src/phase_damage.f90
View File

@ -151,18 +151,18 @@ end function phase_damage_constitutive
!--------------------------------------------------------------------------------------------------
!> @brief returns the degraded/modified elasticity matrix
!--------------------------------------------------------------------------------------------------
module function phase_damage_C66(C66_homogenized,ph,en) result(C66)
module function phase_damage_C66(C66,ph,en)
real(pReal), dimension(6,6), intent(in) :: C66_homogenized
real(pReal), dimension(6,6), intent(in) :: C66
integer, intent(in) :: ph,en
real(pReal), dimension(6,6) :: C66
real(pReal), dimension(6,6) :: phase_damage_C66
damageType: select case (phase_damage(ph))
case (DAMAGE_ISOBRITTLE_ID) damageType
C66 = C66_homogenized * damage_phi(ph,en)**2
phase_damage_C66 = C66 * damage_phi(ph,en)**2
case default damageType
C66 = C66_homogenized
phase_damage_C66 = C66
end select damageType
end function phase_damage_C66