From 6c92e8d2cc533c9eafb8bd385a6479f262f79139 Mon Sep 17 00:00:00 2001 From: Sharan Roongta Date: Thu, 29 Apr 2021 16:16:51 +0200 Subject: [PATCH] belongs to elastic submodule --- src/phase_mechanical.f90 | 20 -------------------- src/phase_mechanical_elastic.f90 | 19 +++++++++++++++++++ 2 files changed, 19 insertions(+), 20 deletions(-) diff --git a/src/phase_mechanical.f90 b/src/phase_mechanical.f90 index 00b8e3667..29e9e8ada 100644 --- a/src/phase_mechanical.f90 +++ b/src/phase_mechanical.f90 @@ -1015,26 +1015,6 @@ module subroutine mechanical_forward() end subroutine mechanical_forward - -!-------------------------------------------------------------------------------------------------- -!> @brief returns the homogenize elasticity matrix -!> ToDo: homogenizedC66 would be more consistent -!-------------------------------------------------------------------------------------------------- -module function phase_homogenizedC(ph,en) result(C) - - real(pReal), dimension(6,6) :: C - integer, intent(in) :: ph, en - - plasticType: select case (phase_plasticity(ph)) - case (PLASTICITY_DISLOTWIN_ID) plasticType - C = plastic_dislotwin_homogenizedC(ph,en) - case default plasticType - C = lattice_C66(1:6,1:6,ph) - end select plasticType - -end function phase_homogenizedC - - !-------------------------------------------------------------------------------------------------- !> @brief calculate stress (P) !-------------------------------------------------------------------------------------------------- diff --git a/src/phase_mechanical_elastic.f90 b/src/phase_mechanical_elastic.f90 index 0e11c3888..9c18bfa7f 100644 --- a/src/phase_mechanical_elastic.f90 +++ b/src/phase_mechanical_elastic.f90 @@ -113,4 +113,23 @@ module subroutine phase_hooke_SandItsTangents(S, dS_dFe, dS_dFi, & end subroutine phase_hooke_SandItsTangents +!-------------------------------------------------------------------------------------------------- +!> @brief returns the homogenized elasticity matrix +!> ToDo: homogenizedC66 would be more consistent +!-------------------------------------------------------------------------------------------------- +module function phase_homogenizedC(ph,en) result(C) + + real(pReal), dimension(6,6) :: C + integer, intent(in) :: ph, en + + plasticType: select case (phase_plasticity(ph)) + case (PLASTICITY_DISLOTWIN_ID) plasticType + C = plastic_dislotwin_homogenizedC(ph,en) + case default plasticType + C = lattice_C66(1:6,1:6,ph) + end select plasticType + +end function phase_homogenizedC + + end submodule elastic