kinematics_thermal_expansion.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/kinematics_thermal_expansion.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/kinematics_thermal_expansion.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module kinematics_thermal_expansion
  use prec
  use io
  use config
  use debug
  use math
  use lattice
  use material
 
  implicit none
  private
 
  integer, dimension(:), allocatable :: kinematics_thermal_expansion_instance
 
  type :: tparameters
    real(preal) :: &
      t_ref
    real(preal), dimension(3,3,3) :: &
      expansion = 0.0_preal
  end type tparameters
 
  type(tparameters), dimension(:), allocatable :: param
 
  public :: &
    kinematics_thermal_expansion_init, &
    kinematics_thermal_expansion_initialstrain, &
    kinematics_thermal_expansion_lianditstangent
 
contains
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine kinematics_thermal_expansion_init
 
  integer :: ninstance,p,i
  real(preal), dimension(:), allocatable :: temp
 
  write(6,'(/,a)') ' <<<+-  kinematics_'//kinematics_thermal_expansion_label//' init  -+>>>'; flush(6)
 
  ninstance = count(phase_kinematics == kinematics_thermal_expansion_id)
  if (iand(debug_level(debug_constitutive),debug_levelbasic) /= 0) &
    write(6,'(a16,1x,i5,/)') '# instances:',ninstance
 
  allocate(kinematics_thermal_expansion_instance(size(config_phase)), source=0)
  allocate(param(ninstance))
 
  do p = 1, size(config_phase)
    kinematics_thermal_expansion_instance(p) = count(phase_kinematics(:,1:p) == kinematics_thermal_expansion_id)
    if (all(phase_kinematics(:,p) /= kinematics_thermal_expansion_id)) cycle
 
    associate(prm => param(kinematics_thermal_expansion_instance(p)), &
              config => config_phase(p))
 
    prm%T_ref = config%getFloat('reference_temperature', defaultval=0.0_preal)
 
    ! read up to three parameters (constant, linear, quadratic with T)
    temp = config%getFloats('thermal_expansion11')
    prm%expansion(1,1,1:size(temp)) = temp
    temp = config%getFloats('thermal_expansion22',defaultval=[(0.0_preal, i=1,size(temp))],requiredsize=size(temp))
    prm%expansion(2,2,1:size(temp)) = temp
    temp = config%getFloats('thermal_expansion33',defaultval=[(0.0_preal, i=1,size(temp))],requiredsize=size(temp))
    prm%expansion(3,3,1:size(temp)) = temp
    do i=1, size(prm%expansion,3)
      prm%expansion(1:3,1:3,i) = lattice_applylatticesymmetry33(prm%expansion(1:3,1:3,i),config%getString('lattice_structure'))
    enddo
 
    end associate
  enddo
 
end subroutine kinematics_thermal_expansion_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
pure function kinematics_thermal_expansion_initialstrain(homog,phase,offset)
 
  integer, intent(in) :: &
    phase, &
    homog, &
    offset
 
  real(preal), dimension(3,3) :: &
    kinematics_thermal_expansion_initialstrain
 
  associate(prm => param(kinematics_thermal_expansion_instance(phase)))
  kinematics_thermal_expansion_initialstrain = &
    (temperature(homog)%p(offset) - prm%T_ref)**1 / 1. * prm%expansion(1:3,1:3,1) + &               ! constant  coefficient
    (temperature(homog)%p(offset) - prm%T_ref)**2 / 2. * prm%expansion(1:3,1:3,2) + &               ! linear    coefficient
    (temperature(homog)%p(offset) - prm%T_ref)**3 / 3. * prm%expansion(1:3,1:3,3)                   ! quadratic coefficient
  end associate
 
end function kinematics_thermal_expansion_initialstrain
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine kinematics_thermal_expansion_lianditstangent(Li, dLi_dTstar, ipc, ip, el)
 
  integer, intent(in) :: &
    ipc, &                                                                                          !< grain number
    ip, &                                                                                           !< integration point number
    el
  real(preal),   intent(out), dimension(3,3) :: &
    li
  real(preal),   intent(out), dimension(3,3,3,3) :: &
    dli_dtstar
 
  integer :: &
    phase, &
    homog
  real(preal) :: &
    t, tdot
 
  phase = material_phaseat(ipc,el)
  homog = material_homogenizationat(el)
  t = temperature(homog)%p(thermalmapping(homog)%p(ip,el))
  tdot = temperaturerate(homog)%p(thermalmapping(homog)%p(ip,el))
 
  associate(prm => param(kinematics_thermal_expansion_instance(phase)))
  li = tdot * ( &
                prm%expansion(1:3,1:3,1)*(t - prm%T_ref)**0 &                                       ! constant  coefficient
              + prm%expansion(1:3,1:3,2)*(t - prm%T_ref)**1 &                                       ! linear    coefficient
              + prm%expansion(1:3,1:3,3)*(t - prm%T_ref)**2 &                                       ! quadratic coefficient
              ) / &
       (1.0_preal &
             + prm%expansion(1:3,1:3,1)*(t - prm%T_ref)**1 / 1. &
             + prm%expansion(1:3,1:3,2)*(t - prm%T_ref)**2 / 2. &
             + prm%expansion(1:3,1:3,3)*(t - prm%T_ref)**3 / 3. &
       )
  end associate
  dli_dtstar = 0.0_preal
 
end subroutine kinematics_thermal_expansion_lianditstangent
 
end module kinematics_thermal_expansion