fstr_FiniteRotationKinematics.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
 
module m_fstr_finiterotationkinematics
  use hecmw
  use elementinfo, only: fe_mitc4_shell
  use m_fstr, only: fstr_solid
  use mmaterial, only: tmaterial, totallag, iselastic
  implicit none
 
  private
 
  public :: fstr_uses_finite_rotation_kinematics
  public :: fstr_has_finite_rotation_kinematics
  public :: fstr_mark_finite_rotation_nodes
 
contains
 
  logical function fstr_uses_finite_rotation_kinematics( etype, nn, material )
    integer(kind=kint), intent(in) :: etype
    integer(kind=kint), intent(in) :: nn
    type(tmaterial), intent(in)    :: material
 
    fstr_uses_finite_rotation_kinematics = ( etype == fe_mitc4_shell .and. nn == 4 &
      .and. material%nlgeom_flag == totallag .and. iselastic( material%mtype ) )
  end function fstr_uses_finite_rotation_kinematics
 
  logical function fstr_has_finite_rotation_kinematics( hecMESH, fstrSOLID )
    type (hecmwst_local_mesh), intent(in) :: hecmesh
    type (fstr_solid), intent(in)         :: fstrsolid
 
    integer(kind=kint) :: itype, is, ie, ic_type, icel, iis, nn
 
    fstr_has_finite_rotation_kinematics = .false.
    if( hecmesh%n_dof < 6 ) return
    if( .not. associated( fstrsolid%elements ) ) return
 
    do itype = 1, hecmesh%n_elem_type
      is = hecmesh%elem_type_index(itype-1) + 1
      ie = hecmesh%elem_type_index(itype)
      ic_type = hecmesh%elem_type_item(itype)
      if( ic_type /= fe_mitc4_shell ) cycle
 
      do icel = is, ie
        iis = hecmesh%elem_node_index(icel-1)
        nn = hecmesh%elem_node_index(icel) - iis
        if( .not. associated( fstrsolid%elements(icel)%gausses ) ) cycle
        if( fstr_uses_finite_rotation_kinematics( ic_type, nn, &
            fstrsolid%elements(icel)%gausses(1)%pMaterial ) ) then
          fstr_has_finite_rotation_kinematics = .true.
          return
        endif
      end do
    end do
  end function fstr_has_finite_rotation_kinematics
 
  subroutine fstr_mark_finite_rotation_nodes( hecMESH, fstrSOLID, ndof, shell_node_mode )
    type (hecmwst_local_mesh), intent(in) :: hecmesh
    type (fstr_solid), intent(in)         :: fstrsolid
    integer(kind=kint), intent(in)        :: ndof
    integer(kind=kint), intent(out)       :: shell_node_mode(:)
 
    integer(kind=kint) :: itype, is, ie, ic_type, icel, iis, nn, j, node_id
 
    shell_node_mode(:) = 0
    if( ndof < 6 ) return
    if( .not. associated( fstrsolid%elements ) ) return
 
    do itype = 1, hecmesh%n_elem_type
      is = hecmesh%elem_type_index(itype-1) + 1
      ie = hecmesh%elem_type_index(itype)
      ic_type = hecmesh%elem_type_item(itype)
      if( ic_type /= fe_mitc4_shell ) cycle
 
      do icel = is, ie
        iis = hecmesh%elem_node_index(icel-1)
        nn = hecmesh%elem_node_index(icel) - iis
        if( .not. associated( fstrsolid%elements(icel)%gausses ) ) cycle
        if( .not. fstr_uses_finite_rotation_kinematics( ic_type, nn, &
            fstrsolid%elements(icel)%gausses(1)%pMaterial ) ) cycle
 
        do j = 1, nn
          node_id = hecmesh%elem_node_item(iis+j)
          if( node_id > 0 .and. node_id <= size(shell_node_mode) ) shell_node_mode(node_id) = 1
        end do
      end do
    end do
  end subroutine fstr_mark_finite_rotation_nodes
 
end module m_fstr_finiterotationkinematics