geometry_plastic_nonlocal.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/geometry_plastic_nonlocal.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/geometry_plastic_nonlocal.f90"
!--------------------------------------------------------------------------------------------------
! plasticity model
!--------------------------------------------------------------------------------------------------
module geometry_plastic_nonlocal
  use prec
  use results
 
  implicit none
  public
  
  integer, protected :: &
    geometry_plastic_nonlocal_nipneighbors
  
  integer,     dimension(:,:,:,:), allocatable, protected :: &
    geometry_plastic_nonlocal_ipneighborhood
 
  real(preal), dimension(:,:),     allocatable, protected :: &
    geometry_plastic_nonlocal_ipvolume0
 
  real(preal), dimension(:,:,:),   allocatable, protected :: &
    geometry_plastic_nonlocal_iparea0
  
  real(preal), dimension(:,:,:,:), allocatable, protected :: &
    geometry_plastic_nonlocal_ipareanormal0
  
 
contains
 
!---------------------------------------------------------------------------------------------------
!            face ID (second index) gives the element ID (1 @ first index), IP ID (2 @ first index)
!            and face ID (3 @ first index).
!            A triangle (2D) has 3 faces, a quadrilateral (2D) had 4 faces, a tetrahedron (3D) has
!            4 faces, and a hexahedron (3D) has 6 faces.
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_setipneighborhood(IPneighborhood)
 
  integer, dimension(:,:,:,:), intent(in) :: IPneighborhood
 
  geometry_plastic_nonlocal_ipneighborhood = ipneighborhood
  geometry_plastic_nonlocal_nipneighbors   = size(ipneighborhood,2)
  
 
end subroutine geometry_plastic_nonlocal_setipneighborhood
 
 
!---------------------------------------------------------------------------------------------------
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_setipvolume(IPvolume)
 
  real(pReal), dimension(:,:), intent(in) :: IPvolume
 
  geometry_plastic_nonlocal_ipvolume0 = ipvolume
 
end subroutine geometry_plastic_nonlocal_setipvolume
 
 
!---------------------------------------------------------------------------------------------------
!         encompassing an integration point
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_setiparea(IParea)
 
  real(pReal), dimension(:,:,:), intent(in) :: IParea
 
  geometry_plastic_nonlocal_iparea0 = iparea
 
end subroutine geometry_plastic_nonlocal_setiparea
 
 
!---------------------------------------------------------------------------------------------------
!         encompassing an integration point
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_setipareanormal(IPareaNormal)
 
  real(pReal), dimension(:,:,:,:), intent(in) :: IPareaNormal
 
  geometry_plastic_nonlocal_ipareanormal0 = ipareanormal
 
end subroutine geometry_plastic_nonlocal_setipareanormal
 
 
!---------------------------------------------------------------------------------------------------
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_disable
 
  if(allocated(geometry_plastic_nonlocal_ipneighborhood)) &
    deallocate(geometry_plastic_nonlocal_ipneighborhood)
    
  if(allocated(geometry_plastic_nonlocal_ipvolume0)) &
    deallocate(geometry_plastic_nonlocal_ipvolume0)
    
  if(allocated(geometry_plastic_nonlocal_iparea0)) &
    deallocate(geometry_plastic_nonlocal_iparea0)
    
  if(allocated(geometry_plastic_nonlocal_ipareanormal0)) &
    deallocate(geometry_plastic_nonlocal_ipareanormal0)
  
end subroutine geometry_plastic_nonlocal_disable
 
 
!---------------------------------------------------------------------------------------------------
!---------------------------------------------------------------------------------------------------
subroutine geometry_plastic_nonlocal_results
  
  integer,     dimension(:),   allocatable :: shp
 
  call results_openjobfile
 
  writevolume: block
    real(pReal), dimension(:), allocatable :: temp
    shp = shape(geometry_plastic_nonlocal_ipvolume0)
    temp = reshape(geometry_plastic_nonlocal_ipvolume0,[shp(1)*shp(2)])
    call results_writedataset('geometry',temp,'v_0',&
                              'initial cell volume','m³')
  end block writevolume
 
  writeareas: block
    real(pReal), dimension(:,:), allocatable :: temp
    shp = shape(geometry_plastic_nonlocal_iparea0)
    temp = reshape(geometry_plastic_nonlocal_iparea0,[shp(1),shp(2)*shp(3)])
    call results_writedataset('geometry',temp,'a_0',&
                              'initial cell face area','m²')
  end block writeareas
 
  writenormals: block
    real(pReal), dimension(:,:,:), allocatable :: temp
    shp = shape(geometry_plastic_nonlocal_ipareanormal0)
    temp = reshape(geometry_plastic_nonlocal_ipareanormal0,[shp(1),shp(2),shp(3)*shp(4)])
    call results_writedataset('geometry',temp,'n_0',&
                              'initial cell face normals','-',transposed=.false.)
  end block writenormals
 
 
  call results_closejobfile
  
end subroutine geometry_plastic_nonlocal_results
 
end module geometry_plastic_nonlocal