usage example
This commit is contained in:
parent
4bce087d3d
commit
0dbc6fb435
|
@ -14,6 +14,7 @@
|
||||||
#include "Lambert.f90"
|
#include "Lambert.f90"
|
||||||
#include "rotations.f90"
|
#include "rotations.f90"
|
||||||
#include "FEsolving.f90"
|
#include "FEsolving.f90"
|
||||||
|
#include "geometry_plastic_nonlocal.f90"
|
||||||
#include "element.f90"
|
#include "element.f90"
|
||||||
#include "mesh_base.f90"
|
#include "mesh_base.f90"
|
||||||
#ifdef Abaqus
|
#ifdef Abaqus
|
||||||
|
@ -28,7 +29,6 @@
|
||||||
#endif
|
#endif
|
||||||
#include "material.f90"
|
#include "material.f90"
|
||||||
#include "lattice.f90"
|
#include "lattice.f90"
|
||||||
#include "geometry_plastic_nonlocal.f90"
|
|
||||||
#include "source_thermal_dissipation.f90"
|
#include "source_thermal_dissipation.f90"
|
||||||
#include "source_thermal_externalheat.f90"
|
#include "source_thermal_externalheat.f90"
|
||||||
#include "source_damage_isoBrittle.f90"
|
#include "source_damage_isoBrittle.f90"
|
||||||
|
|
|
@ -25,4 +25,28 @@ module geometry_plastic_nonlocal
|
||||||
real(pReal),dimension(:,:,:,:), allocatable, public, protected :: &
|
real(pReal),dimension(:,:,:,:), allocatable, public, protected :: &
|
||||||
geometry_plastic_nonlocal_IPareaNormal !< area normal of interface to neighboring IP (initially!)
|
geometry_plastic_nonlocal_IPareaNormal !< area normal of interface to neighboring IP (initially!)
|
||||||
|
|
||||||
|
public :: &
|
||||||
|
geometry_plastic_nonlocal_set_IPneighborhood, &
|
||||||
|
geometry_plastic_nonlocal_set_IPvolume
|
||||||
|
|
||||||
|
contains
|
||||||
|
|
||||||
|
subroutine geometry_plastic_nonlocal_set_IPneighborhood(IPneighborhood)
|
||||||
|
|
||||||
|
integer, dimension(:,:,:,:), intent(in) :: IPneighborhood
|
||||||
|
|
||||||
|
geometry_plastic_nonlocal_IPneighborhood = IPneighborhood
|
||||||
|
|
||||||
|
end subroutine geometry_plastic_nonlocal_set_IPneighborhood
|
||||||
|
|
||||||
|
|
||||||
|
subroutine geometry_plastic_nonlocal_set_IPvolume(IPvolume)
|
||||||
|
|
||||||
|
real(pReal), dimension(:,:), intent(in) :: IPvolume
|
||||||
|
|
||||||
|
geometry_plastic_nonlocal_IPvolume = IPvolume
|
||||||
|
|
||||||
|
end subroutine geometry_plastic_nonlocal_set_IPvolume
|
||||||
|
|
||||||
|
|
||||||
end module geometry_plastic_nonlocal
|
end module geometry_plastic_nonlocal
|
||||||
|
|
|
@ -7,7 +7,8 @@
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
module mesh
|
module mesh
|
||||||
use, intrinsic :: iso_c_binding
|
use, intrinsic :: iso_c_binding
|
||||||
use prec, only: pReal, pInt
|
use prec
|
||||||
|
use geometry_plastic_nonlocal
|
||||||
use mesh_base
|
use mesh_base
|
||||||
|
|
||||||
implicit none
|
implicit none
|
||||||
|
@ -191,6 +192,7 @@ subroutine mesh_init(ip,el)
|
||||||
if (myDebug) write(6,'(a)') ' Built IP areas'; flush(6)
|
if (myDebug) write(6,'(a)') ' Built IP areas'; flush(6)
|
||||||
|
|
||||||
call mesh_spectral_build_ipNeighborhood
|
call mesh_spectral_build_ipNeighborhood
|
||||||
|
call geometry_plastic_nonlocal_set_IPneighborhood(mesh_ipNeighborhood)
|
||||||
|
|
||||||
if (myDebug) write(6,'(a)') ' Built IP neighborhood'; flush(6)
|
if (myDebug) write(6,'(a)') ' Built IP neighborhood'; flush(6)
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue