SpaceTimeDeformClass.f90 Source File

module SpaceTimeDeformClass
   use, intrinsic :: iso_fortran_env
   use MathClass
   use LinearSolverClass
   use FEMDomainClass
   use PostProcessingClass
   use ConstitutiveModelClass

   type :: SpaceTimeFEM_
      type(STFEMDomain_), pointer ::STFEMDomain
      real(real64), allocatable ::DeformStress(:, :, :)
      real(real64), allocatable ::DeformStrain(:, :, :)
      real(real64), allocatable ::DeformStressInit(:, :, :)
      real(real64), allocatable ::DeformStressMat(:, :, :)
      real(real64), allocatable ::DeformStressRHS(:, :)
      real(real64), allocatable ::DeformVecEBETot(:, :)
      real(real64), allocatable ::DeformVecEBEInc(:, :)

      real(real64), allocatable ::DeformVecGloTot(:)
      real(real64), allocatable ::DeformVecGloInc(:)

      real(real64), allocatable ::TractionVecGlo(:)
      real(real64), allocatable ::ResidualVecGlo(:)
      real(real64), allocatable ::InternalVecGlo(:)

      real(real64), allocatable ::VolInitCurrEBE(:, :)
      real(real64)             ::dt, error

      logical :: MeshMove
   contains
      procedure :: SetMovingMesh => SetMovingMesh
      procedure :: Init => Init
      procedure :: SetMatrix => SetMatrix
      procedure :: SetRHS => SetRHS
   end type

contains

! ############################################################
   subroutine SetMovingMesh(obj, on_off)
      class(SpaceTimeFEM_), intent(inout)  :: obj
      logical, intent(in)::on_off

      obj%MeshMove = on_off

   end subroutine
! ############################################################

! ############################################################
   subroutine Init(obj)
      class(SpaceTimeFEM_), intent(inout)  :: obj

      obj%MeshMove = .false.

   end subroutine
! ############################################################

! ############################################################
   subroutine SetMatrix(obj)
      class(SpaceTimeFEM_), intent(inout)  :: obj

   end subroutine
! ############################################################

! ############################################################
   subroutine GetSTShapeFunc(obj, ElemID, GpID)
      class(SpaceTimeFEM_), intent(inout)  :: obj
      integer(int32), optional, intent(in)::ElemID
      integer(int32), optional, intent(in)::GpID

      if (obj%MeshMove .eqv. .false.) then
         call obj%STFEMDomain%ShapeFunction%GetAll(elem_id=ElemID, nod_coord=obj%STFEMDomain%Mesh%NodCoord, &
                                                   elem_nod=obj%STFEMDomain%Mesh%ElemNod, OptionalGpID=GpID)
         call obj%STFEMDomain%TimeShapeFunction%GetAll(elem_id=ElemID, nod_coord=obj%STFEMDomain%Mesh%NodCoord, &
                                                       elem_nod=obj%STFEMDomain%Mesh%ElemNod, OptionalGpID=GpID)
      elseif (obj%MeshMove .eqv. .true.) then
         call obj%STFEMDomain%ShapeFunction%GetAll(elem_id=ElemID, nod_coord=obj%STFEMDomain%Mesh%NodCoord, &
                                                   elem_nod=obj%STFEMDomain%Mesh%ElemNod, OptionalGpID=GpID)
         call obj%STFEMDomain%TimeShapeFunction%GetAll(elem_id=ElemID, nod_coord=obj%STFEMDomain%Mesh%NodCoord, &
                                                       elem_nod=obj%STFEMDomain%Mesh%ElemNod, OptionalGpID=GpID)
      else
         stop "ERROR :: SetElementMatrix >> please run constructor obj%Init to create instance"
      end if

   end subroutine
! ############################################################

! ############################################################
   subroutine SetElementMatrix(obj)
      class(SpaceTimeFEM_), intent(inout)  :: obj

   end subroutine
! ############################################################

! ############################################################
   subroutine SetSpaceTimeShapeFunction(Nobj, Tobj)
      class(ShapeFunction_), intent(inout)  :: Nobj, Tobj

   end subroutine

! ############################################################

! ############################################################
   subroutine SetRHS(obj)
      class(SpaceTimeFEM_), intent(inout)  :: obj

   end subroutine
! ############################################################

end module