dummy data layout

This commit is contained in:
Martin Diehl 2020-12-30 12:00:47 +01:00
parent dd23bec9aa
commit 0dac5f84ef
4 changed files with 56 additions and 0 deletions

View File

@ -202,11 +202,17 @@ module constitutive
real(pReal) :: T real(pReal) :: T
end function constitutive_thermal_T end function constitutive_thermal_T
module subroutine constitutive_mech_setF(F,co,ip,el) module subroutine constitutive_mech_setF(F,co,ip,el)
real(pReal), dimension(3,3), intent(in) :: F real(pReal), dimension(3,3), intent(in) :: F
integer, intent(in) :: co, ip, el integer, intent(in) :: co, ip, el
end subroutine constitutive_mech_setF end subroutine constitutive_mech_setF
module subroutine constitutive_thermal_setT(T,co,ip,el)
real(pReal), intent(in) :: T
integer, intent(in) :: co, ip, el
end subroutine constitutive_thermal_setT
! == cleaned:end =================================================================================== ! == cleaned:end ===================================================================================
module function crystallite_stress(dt,co,ip,el) result(converged_) module function crystallite_stress(dt,co,ip,el) result(converged_)
@ -414,6 +420,7 @@ module constitutive
constitutive_restartRead, & constitutive_restartRead, &
integrateSourceState, & integrateSourceState, &
constitutive_mech_setF, & constitutive_mech_setF, &
constitutive_thermal_setT, &
constitutive_mech_getP, & constitutive_mech_getP, &
constitutive_mech_getLp, & constitutive_mech_getLp, &
constitutive_mech_getF, & constitutive_mech_getF, &

View File

@ -139,4 +139,11 @@ module function constitutive_thermal_T(co,ip,el) result(T)
end function constitutive_thermal_T end function constitutive_thermal_T
! setter for homogenization
module subroutine constitutive_thermal_setT(T,co,ip,el)
real(pReal), intent(in) :: T
integer, intent(in) :: co, ip, el
end subroutine constitutive_thermal_setT
end submodule constitutive_thermal end submodule constitutive_thermal

View File

@ -27,6 +27,8 @@ module homogenization
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! General variables for the homogenization at a material point ! General variables for the homogenization at a material point
real(pReal), dimension(:), allocatable, public :: &
homogenization_T
real(pReal), dimension(:,:,:), allocatable, public :: & real(pReal), dimension(:,:,:), allocatable, public :: &
homogenization_F0, & !< def grad of IP at start of FE increment homogenization_F0, & !< def grad of IP at start of FE increment
homogenization_F !< def grad of IP to be reached at end of FE increment homogenization_F !< def grad of IP to be reached at end of FE increment
@ -56,6 +58,9 @@ module homogenization
num_homog !< pointer to mechanical homogenization numerics data num_homog !< pointer to mechanical homogenization numerics data
end subroutine mech_init end subroutine mech_init
module subroutine thermal_init
end subroutine thermal_init
module subroutine mech_partition(subF,ip,el) module subroutine mech_partition(subF,ip,el)
real(pReal), intent(in), dimension(3,3) :: & real(pReal), intent(in), dimension(3,3) :: &
subF subF

View File

@ -0,0 +1,37 @@
!--------------------------------------------------------------------------------------------------
!> @author Martin Diehl, KU Leuven
!--------------------------------------------------------------------------------------------------
submodule(homogenization) homogenization_thermal
contains
!--------------------------------------------------------------------------------------------------
!> @brief Allocate variables and set parameters.
!--------------------------------------------------------------------------------------------------
module subroutine thermal_init()
print'(/,a)', ' <<<+- homogenization_thermal init -+>>>'
allocate(homogenization_T(discretization_nIPs*discretization_Nelems), source=0.0_pReal)
end subroutine thermal_init
!--------------------------------------------------------------------------------------------------
!> @brief Partition T onto the individual constituents.
!--------------------------------------------------------------------------------------------------
module subroutine thermal_partition(T,ip,el)
real(pReal), intent(in) :: T
integer, intent(in) :: &
ip, & !< integration point
el !< element number
call constitutive_thermal_setT(T,1,ip,el)
end subroutine thermal_partition
end submodule homogenization_thermal