fstr_contact_param.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2021 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
module mcontactparam
 
  use hecmw
 
  implicit none
 
  private
  public :: tcontactparam
  public :: init_contactparam
  public :: tcontactinterference
  public :: init_contact_if
 
  type tcontactparam
    character(HECMW_NAME_LEN) :: name
    real(kind=kreal) :: clearance        
    real(kind=kreal) :: clr_same_elem    
    real(kind=kreal) :: clr_difflpos     
    real(kind=kreal) :: clr_cal_norm     
    real(kind=kreal) :: distclr_init     
    real(kind=kreal) :: distclr_free     
    real(kind=kreal) :: distclr_nocheck  
    real(kind=kreal) :: tensile_force    
    real(kind=kreal) :: box_exp_rate     
  end type tcontactparam
 
  type tcontactinterference
    integer(kind=kint) :: if_type
    real(kind=kreal)   :: etime
    real(kind=kreal)   :: initial_pos
    real(kind=kreal)   :: end_pos
    character(HECMW_NAME_LEN) :: cp_name
  end type tcontactinterference
 
contains
 
  subroutine init_contactparam( cparam )
    type(tcontactparam), intent(out) :: cparam
    cparam%name            = ''
    cparam%CLEARANCE       = 1.d-4
    cparam%CLR_SAME_ELEM   = 5.d-3
    cparam%CLR_DIFFLPOS    = 1.d-2
    cparam%CLR_CAL_NORM    = 1.d-4
    cparam%DISTCLR_INIT    = 1.d-6
    cparam%DISTCLR_FREE    =-1.d-6
    cparam%DISTCLR_NOCHECK = 1.d0
    cparam%TENSILE_FORCE   =-1.d-8
    cparam%BOX_EXP_RATE    = 1.05d0
  end subroutine init_contactparam
 
  subroutine init_contact_if( contact_if )
    type(tcontactinterference), intent(out) :: contact_if
    contact_if%if_type     = 1
    contact_if%etime       = 1.0d0
    contact_if%initial_pos = 0.d0
    contact_if%end_pos     = 0.d0
    contact_if%cp_name = ''
  end subroutine init_contact_if
 
end module mcontactparam