damage_local.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/damage_local.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/damage_local.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module damage_local
  use prec
  use material
  use config
  use numerics
  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_local_init, &
    damage_local_updatestate, &
    damage_local_results
 
contains
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_local_init
 
  integer :: ninstance,nofmyhomog,h
 
  write(6,'(/,a)') ' <<<+-  damage_'//damage_local_label//' init  -+>>>'; flush(6)
 
  ninstance = count(damage_type == damage_local_id)
  allocate(param(ninstance))
 
  do h = 1, size(config_homogenization)
    if (damage_type(h) /= damage_local_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_local_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function damage_local_updatestate(subdt, ip, el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal),   intent(in) :: &
    subdt
  logical,    dimension(2)  :: &
    damage_local_updatestate
  integer :: &
    homog, &
    offset
  real(preal) :: &
    phi, phidot, dphidot_dphi
 
  homog  = material_homogenizationat(el)
  offset = material_homogenizationmemberat(ip,el)
  phi = damagestate(homog)%subState0(1,offset)
  call damage_local_getsourceanditstangent(phidot, dphidot_dphi, phi, ip, el)
  phi = max(residualstiffness,min(1.0_preal,phi + subdt*phidot))
 
  damage_local_updatestate = [     abs(phi - damagestate(homog)%state(1,offset)) &
                                <= err_damage_tolabs &
                              .or. abs(phi - damagestate(homog)%state(1,offset)) &
                                <= err_damage_tolrel*abs(damagestate(homog)%state(1,offset)), &
                              .true.]
 
  damagestate(homog)%state(1,offset) = phi
 
end function damage_local_updatestate
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_local_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_local_getsourceanditstangent
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine damage_local_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_local_results
 
 
end module damage_local