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"
    endif

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