source_thermal_dissipation.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/source_thermal_dissipation.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/source_thermal_dissipation.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module source_thermal_dissipation
  use prec
  use debug
  use discretization
  use material
  use config
 
  implicit none
  private
 
  integer,           dimension(:),   allocatable :: &
    source_thermal_dissipation_offset, &                                                            !< which source is my current thermal dissipation mechanism?
    source_thermal_dissipation_instance
 
  type :: tparameters                                                                               
    real(preal) :: &
      kappa
  end type tparameters
 
  type(tparameters), dimension(:),   allocatable :: param
 
 
  public :: &
    source_thermal_dissipation_init, &
    source_thermal_dissipation_getrateanditstangent
 
contains
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_dissipation_init
 
  integer :: ninstance,sourceoffset,nipcmyphase,p
 
  write(6,'(/,a)') ' <<<+-  source_'//source_thermal_dissipation_label//' init  -+>>>'; flush(6)
 
  ninstance = count(phase_source == source_thermal_dissipation_id)
  if (iand(debug_level(debug_constitutive),debug_levelbasic) /= 0) &
    write(6,'(a16,1x,i5,/)') '# instances:',ninstance
 
  allocate(source_thermal_dissipation_offset(size(config_phase)), source=0)
  allocate(source_thermal_dissipation_instance(size(config_phase)), source=0)
  allocate(param(ninstance))
 
  do p = 1, size(config_phase)
    source_thermal_dissipation_instance(p) = count(phase_source(:,1:p) == source_thermal_dissipation_id)
    do sourceoffset = 1, phase_nsources(p)
      if (phase_source(sourceoffset,p) == source_thermal_dissipation_id) then
        source_thermal_dissipation_offset(p) = sourceoffset
        exit
      endif
    enddo
 
    if (all(phase_source(:,p) /= source_thermal_dissipation_id)) cycle
    associate(prm => param(source_thermal_dissipation_instance(p)), &
              config => config_phase(p))
 
    prm%kappa = config%getFloat('dissipation_coldworkcoeff')
 
    nipcmyphase = count(material_phaseat==p) * discretization_nip
    call material_allocatesourcestate(p,sourceoffset,nipcmyphase,0,0,0)
 
    end associate
  enddo
 
end subroutine source_thermal_dissipation_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_dissipation_getrateanditstangent(TDot, dTDot_dT, Tstar, Lp, phase)
 
  integer, intent(in) :: &
    phase
  real(preal),  intent(in), dimension(3,3) :: &
    tstar
  real(preal),  intent(in), dimension(3,3) :: &
    lp
 
  real(preal),  intent(out) :: &
    tdot, &
    dtdot_dt
 
  associate(prm => param(source_thermal_dissipation_instance(phase)))
  tdot = prm%kappa*sum(abs(tstar*lp))
  dtdot_dt = 0.0_preal
  end associate
 
end subroutine source_thermal_dissipation_getrateanditstangent
 
end module source_thermal_dissipation