damage_nonlocal.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/damage_nonlocal.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/damage_nonlocal.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module damage_nonlocal
  use prec
  use material
  use config
  use numerics
  use crystallite
  use lattice
  use source_damage_isobrittle
  use source_damage_isoductile
  use source_damage_anisobrittle
  use source_damage_anisoductile
  use results
 
  implicit none
  private
 
  type :: tparameters
    character(len=pStringLen), allocatable, dimension(:) :: &
      output
  end type tparameters
 
  type(tparameters),             dimension(:), allocatable :: &
    param
 
  public :: &
    damage_nonlocal_init, &
    damage_nonlocal_getsourceanditstangent, &
    damage_nonlocal_getdiffusion, &
    damage_nonlocal_getmobility, &
    damage_nonlocal_putnonlocaldamage, &
    damage_nonlocal_results
 
contains
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_nonlocal_init
 
  integer :: ninstance,nofmyhomog,h
 
  write(6,'(/,a)') ' <<<+-  damage_'//damage_nonlocal_label//' init  -+>>>'; flush(6)
 
  ninstance = count(damage_type == damage_nonlocal_id)
  allocate(param(ninstance))
 
  do h = 1, size(config_homogenization)
    if (damage_type(h) /= damage_nonlocal_id) cycle
    associate(prm => param(damage_typeinstance(h)),config => config_homogenization(h))
 
    prm%output = config%getStrings('(output)',defaultval=emptystringarray)
 
    nofmyhomog = count(material_homogenizationat == h)
    damagestate(h)%sizeState = 1
    allocate(damagestate(h)%state0   (1,nofmyhomog), source=damage_initialphi(h))
    allocate(damagestate(h)%subState0(1,nofmyhomog), source=damage_initialphi(h))
    allocate(damagestate(h)%state    (1,nofmyhomog), source=damage_initialphi(h))
 
    nullify(damagemapping(h)%p)
    damagemapping(h)%p => material_homogenizationmemberat
    deallocate(damage(h)%p)
    damage(h)%p => damagestate(h)%state(1,:)
 
    end associate
  enddo
 
end subroutine damage_nonlocal_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_nonlocal_getsourceanditstangent(phiDot, dPhiDot_dPhi, phi, ip, el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal),   intent(in) :: &
    phi
  integer :: &
    phase, &
    grain, &
    source, &
    constituent
  real(preal) :: &
    phidot, dphidot_dphi, localphidot, dlocalphidot_dphi
 
  phidot = 0.0_preal
  dphidot_dphi = 0.0_preal
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    phase = material_phaseat(grain,el)
    constituent = material_phasememberat(grain,ip,el)
    do source = 1, phase_nsources(phase)
      select case(phase_source(source,phase))
        case (source_damage_isobrittle_id)
         call source_damage_isobrittle_getrateanditstangent  (localphidot, dlocalphidot_dphi, phi, phase, constituent)
 
        case (source_damage_isoductile_id)
         call source_damage_isoductile_getrateanditstangent  (localphidot, dlocalphidot_dphi, phi, phase, constituent)
 
        case (source_damage_anisobrittle_id)
         call source_damage_anisobrittle_getrateanditstangent(localphidot, dlocalphidot_dphi, phi, phase, constituent)
 
        case (source_damage_anisoductile_id)
         call source_damage_anisoductile_getrateanditstangent(localphidot, dlocalphidot_dphi, phi, phase, constituent)
 
        case default
         localphidot = 0.0_preal
         dlocalphidot_dphi = 0.0_preal
 
      end select
      phidot = phidot + localphidot
      dphidot_dphi = dphidot_dphi + dlocalphidot_dphi
    enddo
  enddo
 
  phidot = phidot/real(homogenization_ngrains(material_homogenizationat(el)),preal)
  dphidot_dphi = dphidot_dphi/real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end subroutine damage_nonlocal_getsourceanditstangent
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function damage_nonlocal_getdiffusion(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal), dimension(3,3) :: &
    damage_nonlocal_getdiffusion
  integer :: &
    homog, &
    grain
 
  homog  = material_homogenizationat(el)
  damage_nonlocal_getdiffusion = 0.0_preal
  do grain = 1, homogenization_ngrains(homog)
    damage_nonlocal_getdiffusion = damage_nonlocal_getdiffusion + &
      crystallite_push33toref(grain,ip,el,lattice_damagediffusion(1:3,1:3,material_phaseat(grain,el)))
  enddo
 
  damage_nonlocal_getdiffusion = &
    charlength**2*damage_nonlocal_getdiffusion/real(homogenization_ngrains(homog),preal)
 
end function damage_nonlocal_getdiffusion
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
real(pReal) function damage_nonlocal_getmobility(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  integer :: &
    ipc
 
  damage_nonlocal_getmobility = 0.0_preal
 
  do ipc = 1, homogenization_ngrains(material_homogenizationat(el))
    damage_nonlocal_getmobility = damage_nonlocal_getmobility + lattice_damagemobility(material_phaseat(ipc,el))
  enddo
 
  damage_nonlocal_getmobility = damage_nonlocal_getmobility/&
                                real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end function damage_nonlocal_getmobility
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_nonlocal_putnonlocaldamage(phi,ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal),   intent(in) :: &
    phi
  integer :: &
    homog, &
    offset
 
  homog  = material_homogenizationat(el)
  offset = damagemapping(homog)%p(ip,el)
  damage(homog)%p(offset) = phi
 
end subroutine damage_nonlocal_putnonlocaldamage
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_nonlocal_results(homog,group)
 
  integer,          intent(in) :: homog
  character(len=*), intent(in) :: group
 
  integer :: o
 
  associate(prm => param(damage_typeinstance(homog)))
  outputsloop: do o = 1,size(prm%output)
    select case(prm%output(o))
      case ('damage')
        call results_writedataset(group,damage(homog)%p,'phi',&
                                  'damage indicator','-')
    end select
  enddo outputsloop
  end associate
 
end subroutine damage_nonlocal_results
 
end module damage_nonlocal