From 07ccaf5fe76e05ddcaea988a27cf410bf06974ac Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 21 Dec 2020 12:14:09 +0100 Subject: [PATCH] some ideas ... --- src/constitutive.f90 | 13 +++++++++++++ src/constitutive_mech.f90 | 34 ++++++++++++++++++++++++++++++++++ 2 files changed, 47 insertions(+) diff --git a/src/constitutive.f90 b/src/constitutive.f90 index 58442283a..b5158be4e 100644 --- a/src/constitutive.f90 +++ b/src/constitutive.f90 @@ -140,6 +140,7 @@ module constitutive interface +! == cleaned:begin ================================================================================= module subroutine mech_init end subroutine mech_init @@ -149,6 +150,18 @@ module constitutive module subroutine thermal_init end subroutine thermal_init + + module subroutine mech_results(group,ph) + character(len=*), intent(in) :: group + integer, intent(in) :: ph + end subroutine mech_results + + module subroutine mech_restart_read(fileHandle) + integer(HID_T), intent(in) :: fileHandle + end subroutine mech_restart_read + +! == cleaned:end =================================================================================== + module function constitutive_collectDotState(S, FArray, Fi, FpArray, subdt, ipc, ip, el,phase,of) result(broken) integer, intent(in) :: & diff --git a/src/constitutive_mech.f90 b/src/constitutive_mech.f90 index dea5ed647..acfe3a23b 100644 --- a/src/constitutive_mech.f90 +++ b/src/constitutive_mech.f90 @@ -651,6 +651,40 @@ module function constitutive_deltaState(S, Fi, ipc, ip, el, phase, of) result(br end function constitutive_deltaState +module subroutine mech_results(group,ph) + + character(len=*), intent(in) :: group + integer, intent(in) :: ph + + select case(phase_plasticity(ph)) + + case(PLASTICITY_ISOTROPIC_ID) + call plastic_isotropic_results(phase_plasticityInstance(ph),group//'plastic') + + case(PLASTICITY_PHENOPOWERLAW_ID) + call plastic_phenopowerlaw_results(phase_plasticityInstance(ph),group//'plastic') + + case(PLASTICITY_KINEHARDENING_ID) + call plastic_kinehardening_results(phase_plasticityInstance(ph),group//'plastic') + + case(PLASTICITY_DISLOTWIN_ID) + call plastic_dislotwin_results(phase_plasticityInstance(ph),group//'plastic') + + case(PLASTICITY_DISLOTUNGSTEN_ID) + call plastic_dislotungsten_results(phase_plasticityInstance(ph),group//'plastic') + + case(PLASTICITY_NONLOCAL_ID) + call plastic_nonlocal_results(phase_plasticityInstance(ph),group//'plastic') + end select + +end subroutine mech_results + + module subroutine mech_restart_read(fileHandle) + integer(HID_T), intent(in) :: fileHandle + end subroutine mech_restart_read + + + !-------------------------------------------------------------------------------------------- !> @brief writes plasticity constitutive results to HDF5 output file !--------------------------------------------------------------------------------------------