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

109 lines
4.4 KiB
Fortran
Raw Normal View History

added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
!--------------------------------------------------------------------------------------------------
set svn Id property for recently added files
2014-10-15 19:18:19 +05:30
! $Id$
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief material subroutine for purely elastic material
!--------------------------------------------------------------------------------------------------
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
module thermal_isothermal
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
use prec, only: &
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
pInt, &
pReal
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
implicit none
private
integer(pInt), dimension(:), allocatable, public, protected :: &
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
thermal_isothermal_sizePostResults
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
integer(pInt), dimension(:,:), allocatable, target, public :: &
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
thermal_isothermal_sizePostResult !< size of each post result output
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
real(pReal), dimension(:), allocatable, public :: &
thermal_isothermal_temperature
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
public :: &
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
thermal_isothermal_init
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
contains
!--------------------------------------------------------------------------------------------------
!> @brief module initialization
!> @details reads in material parameters, allocates arrays, and does sanity checks
!--------------------------------------------------------------------------------------------------
no need to pass fileunit during init of none types
2014-11-07 17:45:28 +05:30
subroutine thermal_isothermal_init(temperature_init)
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
use debug, only: &
debug_level, &
debug_constitutive, &
debug_levelBasic
use IO, only: &
IO_timeStamp
use numerics, only: &
only output from the root processor for parallel runs
2014-10-10 01:53:06 +05:30
worldrank, &
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
numerics_integrator
use material, only: &
phase_thermal, &
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
LOCAL_THERMAL_ISOTHERMAL_label, &
LOCAL_THERMAL_ISOTHERMAL_ID, &
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
material_phase, &
thermalState, &
MATERIAL_partPhase
implicit none
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
real(pReal), intent(in) :: temperature_init !< initial temperature
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
integer(pInt) :: &
maxNinstance, &
phase, &
NofMyPhase, &
sizeState, &
sizeDotState
polished field state related changes and brittle damage
2014-09-22 23:45:19 +05:30
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
mainProcess: if (worldrank == 0) then
write(6,'(/,a)') ' <<<+- thermal_'//LOCAL_THERMAL_ISOTHERMAL_label//' init -+>>>'
set svn Id property for recently added files
2014-10-15 19:18:19 +05:30
write(6,'(a)') ' $Id$'
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
#include "compilation_info.f90"
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
endif mainProcess
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
maxNinstance = int(count(phase_thermal == LOCAL_THERMAL_ISOTHERMAL_ID),pInt)
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
if (maxNinstance == 0_pInt) return
if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 0_pInt) &
write(6,'(a16,1x,i5,/)') '# instances:',maxNinstance
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
initializeInstances: do phase = 1_pInt, size(phase_thermal)
NofMyPhase=count(material_phase==phase)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
if (phase_thermal(phase) == LOCAL_THERMAL_ISOTHERMAL_ID) then
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
sizeState = 0_pInt
thermalState(phase)%sizeState = sizeState
sizeDotState = sizeState
thermalState(phase)%sizeDotState = sizeDotState
for consistency, pushed assignment of sizePostResults from aggregators to individual thermal/damage label. Also fixed wrong index while use of mapping
2014-07-04 14:52:57 +05:30
thermalState(phase)%sizePostResults = 0_pInt
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
allocate(thermalState(phase)%state0 (sizeState,NofMyPhase))
allocate(thermalState(phase)%partionedState0(sizeState,NofMyPhase))
allocate(thermalState(phase)%subState0 (sizeState,NofMyPhase))
allocate(thermalState(phase)%state (sizeState,NofMyPhase))
allocate(thermalState(phase)%state_backup (sizeState,NofMyPhase))
allocate(thermalState(phase)%aTolState (NofMyPhase))
allocate(thermalState(phase)%dotState (sizeDotState,NofMyPhase))
allocate(thermalState(phase)%dotState_backup(sizeDotState,NofMyPhase))
if (any(numerics_integrator == 1_pInt)) then
allocate(thermalState(phase)%previousDotState (sizeDotState,NofMyPhase))
allocate(thermalState(phase)%previousDotState2 (sizeDotState,NofMyPhase))
endif
if (any(numerics_integrator == 4_pInt)) &
allocate(thermalState(phase)%RK4dotState (sizeDotState,NofMyPhase))
if (any(numerics_integrator == 5_pInt)) &
allocate(thermalState(phase)%RKCK45dotState (6,sizeDotState,NofMyPhase))
endif
enddo initializeInstances
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
allocate(thermal_isothermal_sizePostResults(maxNinstance), source=0_pInt)
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
allocate(thermal_isothermal_temperature(maxNinstance), source=temperature_init)
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
end subroutine thermal_isothermal_init
added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing. damage_none: does nothing damage_gradient: interacts with solver to solve gradient damage problems thermal_none: does nothing thermai_adiabatic: local heating only thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-25 04:29:16 +05:30
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
end module thermal_isothermal