2015-05-28 22:32:23 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Pratheek Shanthraj, Max-Planck-Institut für Eisenforschung GmbH
!> @brief material subroutine for thermal source due to plastic dissipation
!> @details to be done
!--------------------------------------------------------------------------------------------------
module source_thermal_dissipation
use prec , only : &
pReal , &
pInt
implicit none
private
integer ( pInt ) , dimension ( : ) , allocatable , public , protected :: &
source_thermal_dissipation_sizePostResults , & !< cumulative size of post results
source_thermal_dissipation_offset , & !< which source is my current thermal dissipation mechanism?
source_thermal_dissipation_instance !< instance of thermal dissipation source mechanism
integer ( pInt ) , dimension ( : , : ) , allocatable , target , public :: &
source_thermal_dissipation_sizePostResult !< size of each post result output
character ( len = 64 ) , dimension ( : , : ) , allocatable , target , public :: &
source_thermal_dissipation_output !< name of each post result output
integer ( pInt ) , dimension ( : ) , allocatable , target , public :: &
source_thermal_dissipation_Noutput !< number of outputs per instance of this source
real ( pReal ) , dimension ( : ) , allocatable , private :: &
source_thermal_dissipation_coldworkCoeff
2018-12-31 01:28:38 +05:30
type , private :: tParameters !< container type for internal constitutive parameters
real ( pReal ) :: &
coldworkCoeff
end type tParameters
type ( tParameters ) , dimension ( : ) , allocatable , private :: param !< containers of constitutive parameters (len Ninstance)
2015-05-28 22:32:23 +05:30
public :: &
source_thermal_dissipation_init , &
source_thermal_dissipation_getRateAndItsTangent
contains
!--------------------------------------------------------------------------------------------------
!> @brief module initialization
!> @details reads in material parameters, allocates arrays, and does sanity checks
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_dissipation_init ( fileUnit )
2018-02-02 17:06:09 +05:30
#if defined(__GFORTRAN__) || __INTEL_COMPILER >= 1800
2017-10-05 20:05:34 +05:30
use , intrinsic :: iso_fortran_env , only : &
compiler_version , &
compiler_options
#endif
2015-05-28 22:32:23 +05:30
use debug , only : &
debug_level , &
debug_constitutive , &
debug_levelBasic
use IO , only : &
IO_read , &
IO_lc , &
IO_getTag , &
IO_isBlank , &
IO_stringPos , &
IO_stringValue , &
IO_floatValue , &
IO_intValue , &
IO_warning , &
IO_error , &
IO_timeStamp , &
IO_EOF
use material , only : &
phase_source , &
phase_Nsources , &
phase_Noutput , &
SOURCE_thermal_dissipation_label , &
SOURCE_thermal_dissipation_ID , &
material_phase , &
2018-06-10 21:31:52 +05:30
sourceState
2018-06-14 10:09:49 +05:30
use config , only : &
2018-12-31 01:28:38 +05:30
config_phase , &
2018-06-10 21:31:52 +05:30
material_Nphase , &
2015-05-28 22:32:23 +05:30
MATERIAL_partPhase
use numerics , only : &
numerics_integrator
implicit none
integer ( pInt ) , intent ( in ) :: fileUnit
2015-08-28 13:08:48 +05:30
integer ( pInt ) , allocatable , dimension ( : ) :: chunkPos
2018-12-31 01:28:38 +05:30
integer ( pInt ) :: Ninstance , phase , instance , source , sourceOffset
2015-06-01 21:32:27 +05:30
integer ( pInt ) :: sizeState , sizeDotState , sizeDeltaState
2018-12-31 01:28:38 +05:30
integer ( pInt ) :: NofMyPhase , p
2015-05-28 22:32:23 +05:30
character ( len = 65536 ) :: &
tag = '' , &
line = ''
2016-07-25 23:42:00 +05:30
write ( 6 , '(/,a)' ) ' <<<+- source_' / / SOURCE_thermal_dissipation_label / / ' init -+>>>'
write ( 6 , '(a15,a)' ) ' Current time: ' , IO_timeStamp ( )
2015-05-28 22:32:23 +05:30
#include "compilation_info.f90"
2018-12-31 01:28:38 +05:30
Ninstance = int ( count ( phase_source == SOURCE_thermal_dissipation_ID ) , pInt )
if ( Ninstance == 0_pInt ) return
2015-05-28 22:32:23 +05:30
if ( iand ( debug_level ( debug_constitutive ) , debug_levelBasic ) / = 0_pInt ) &
2018-12-31 01:28:38 +05:30
write ( 6 , '(a16,1x,i5,/)' ) '# instances:' , Ninstance
2015-05-28 22:32:23 +05:30
allocate ( source_thermal_dissipation_offset ( material_Nphase ) , source = 0_pInt )
allocate ( source_thermal_dissipation_instance ( material_Nphase ) , source = 0_pInt )
do phase = 1 , material_Nphase
source_thermal_dissipation_instance ( phase ) = count ( phase_source ( : , 1 : phase ) == SOURCE_thermal_dissipation_ID )
do source = 1 , phase_Nsources ( phase )
if ( phase_source ( source , phase ) == SOURCE_thermal_dissipation_ID ) &
source_thermal_dissipation_offset ( phase ) = source
enddo
enddo
2018-12-31 01:28:38 +05:30
allocate ( source_thermal_dissipation_sizePostResults ( Ninstance ) , source = 0_pInt )
allocate ( source_thermal_dissipation_sizePostResult ( maxval ( phase_Noutput ) , Ninstance ) , source = 0_pInt )
allocate ( source_thermal_dissipation_output ( maxval ( phase_Noutput ) , Ninstance ) )
2015-05-28 22:32:23 +05:30
source_thermal_dissipation_output = ''
2018-12-31 01:28:38 +05:30
allocate ( source_thermal_dissipation_Noutput ( Ninstance ) , source = 0_pInt )
allocate ( source_thermal_dissipation_coldworkCoeff ( Ninstance ) , source = 0.0_pReal )
do p = 1 , size ( config_phase )
if ( all ( phase_source ( : , p ) / = SOURCE_thermal_dissipation_ID ) ) cycle
enddo
2015-05-28 22:32:23 +05:30
rewind ( fileUnit )
phase = 0_pInt
do while ( trim ( line ) / = IO_EOF . and . IO_lc ( IO_getTag ( line , '<' , '>' ) ) / = MATERIAL_partPhase ) ! wind forward to <phase>
line = IO_read ( fileUnit )
enddo
parsingFile : do while ( trim ( line ) / = IO_EOF ) ! read through sections of phase part
line = IO_read ( fileUnit )
if ( IO_isBlank ( line ) ) cycle ! skip empty lines
if ( IO_getTag ( line , '<' , '>' ) / = '' ) then ! stop at next part
line = IO_read ( fileUnit , . true . ) ! reset IO_read
exit
endif
if ( IO_getTag ( line , '[' , ']' ) / = '' ) then ! next phase section
phase = phase + 1_pInt ! advance phase section counter
cycle ! skip to next line
endif
if ( phase > 0_pInt ) then ; if ( any ( phase_source ( : , phase ) == SOURCE_thermal_dissipation_ID ) ) then ! do not short-circuit here (.and. with next if statemen). It's not safe in Fortran
instance = source_thermal_dissipation_instance ( phase ) ! which instance of my source is present phase
2015-08-28 13:08:48 +05:30
chunkPos = IO_stringPos ( line )
tag = IO_lc ( IO_stringValue ( line , chunkPos , 1_pInt ) ) ! extract key
2015-05-28 22:32:23 +05:30
select case ( tag )
case ( 'dissipation_coldworkcoeff' )
2015-08-28 13:08:48 +05:30
source_thermal_dissipation_coldworkCoeff ( instance ) = IO_floatValue ( line , chunkPos , 2_pInt )
2015-05-28 22:32:23 +05:30
end select
endif ; endif
enddo parsingFile
initializeInstances : do phase = 1_pInt , material_Nphase
if ( any ( phase_source ( : , phase ) == SOURCE_thermal_dissipation_ID ) ) then
NofMyPhase = count ( material_phase == phase )
instance = source_thermal_dissipation_instance ( phase )
sourceOffset = source_thermal_dissipation_offset ( phase )
sizeDotState = 0_pInt
2015-06-01 21:32:27 +05:30
sizeDeltaState = 0_pInt
2015-05-28 22:32:23 +05:30
sizeState = 0_pInt
sourceState ( phase ) % p ( sourceOffset ) % sizeState = sizeState
sourceState ( phase ) % p ( sourceOffset ) % sizeDotState = sizeDotState
2015-06-01 21:32:27 +05:30
sourceState ( phase ) % p ( sourceOffset ) % sizeDeltaState = sizeDeltaState
2015-05-28 22:32:23 +05:30
sourceState ( phase ) % p ( sourceOffset ) % sizePostResults = source_thermal_dissipation_sizePostResults ( instance )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % aTolState ( sizeState ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % state0 ( sizeState , NofMyPhase ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % partionedState0 ( sizeState , NofMyPhase ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % subState0 ( sizeState , NofMyPhase ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % state ( sizeState , NofMyPhase ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % dotState ( sizeDotState , NofMyPhase ) , source = 0.0_pReal )
2015-06-01 21:32:27 +05:30
allocate ( sourceState ( phase ) % p ( sourceOffset ) % deltaState ( sizeDeltaState , NofMyPhase ) , source = 0.0_pReal )
2015-05-28 22:32:23 +05:30
if ( any ( numerics_integrator == 1_pInt ) ) then
allocate ( sourceState ( phase ) % p ( sourceOffset ) % previousDotState ( sizeDotState , NofMyPhase ) , source = 0.0_pReal )
allocate ( sourceState ( phase ) % p ( sourceOffset ) % previousDotState2 ( sizeDotState , NofMyPhase ) , source = 0.0_pReal )
endif
if ( any ( numerics_integrator == 4_pInt ) ) &
allocate ( sourceState ( phase ) % p ( sourceOffset ) % RK4dotState ( sizeDotState , NofMyPhase ) , source = 0.0_pReal )
if ( any ( numerics_integrator == 5_pInt ) ) &
allocate ( sourceState ( phase ) % p ( sourceOffset ) % RKCK45dotState ( 6 , sizeDotState , NofMyPhase ) , source = 0.0_pReal )
endif
enddo initializeInstances
end subroutine source_thermal_dissipation_init
!--------------------------------------------------------------------------------------------------
!> @brief returns local vacancy generation rate
!--------------------------------------------------------------------------------------------------
2018-12-31 01:28:38 +05:30
subroutine source_thermal_dissipation_getRateAndItsTangent ( TDot , dTDOT_dT , Tstar_v , Lp , phase , constituent )
2015-05-28 22:32:23 +05:30
use math , only : &
math_Mandel6to33
implicit none
integer ( pInt ) , intent ( in ) :: &
2018-12-31 01:28:38 +05:30
phase , &
constituent
2015-05-28 22:32:23 +05:30
real ( pReal ) , intent ( in ) , dimension ( 6 ) :: &
Tstar_v !< 2nd Piola Kirchhoff stress tensor (Mandel)
real ( pReal ) , intent ( in ) , dimension ( 3 , 3 ) :: &
Lp
real ( pReal ) , intent ( out ) :: &
TDot , &
dTDOT_dT
integer ( pInt ) :: &
2018-12-31 01:28:38 +05:30
instance
2015-05-28 22:32:23 +05:30
instance = source_thermal_dissipation_instance ( phase )
TDot = source_thermal_dissipation_coldworkCoeff ( instance ) * &
sum ( abs ( math_Mandel6to33 ( Tstar_v ) * Lp ) )
dTDOT_dT = 0.0_pReal
end subroutine source_thermal_dissipation_getRateAndItsTangent
end module source_thermal_dissipation