mechgauss.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 mmechgauss
  use hecmw_util
  use mmaterial
  implicit none
 
  ! ----------------------------------------------------------------------------
  type tgaussstatus
    type(tmaterial), pointer  :: pmaterial => null()    
    real(kind=kreal)          :: strain(6)              
    real(kind=kreal)          :: stress(6)              
    integer, pointer          :: istatus(:) =>null()    
    real(kind=kreal), pointer :: fstatus(:) => null()   
    real(kind=kreal)          :: plstrain               
    real(kind=kreal)          :: strain_bak(6)          
    real(kind=kreal)          :: stress_bak(6)          
    real(kind=kreal)          :: nqm(12)                
    real(kind=kreal)          :: strain_out(6)          
    real(kind=kreal)          :: stress_out(6)          
    real(kind=kreal)          :: strain_energy          
    real(kind=kreal)          :: strain_energy_bak      
    real(kind=kreal)          :: plpotential            
  end type
 
  ! ----------------------------------------------------------------------------
  type telement
    integer                     :: etype
    integer                     :: iset
    real(kind=kreal), pointer   :: equiforces(:) => null()  
    type(tgaussstatus), pointer :: gausses(:) => null()  
    real(kind=kreal), pointer   :: aux(:,:) => null()    
    integer                     :: elemact_flag
    real(kind=kreal)            :: elemact_coeff           
  end type
 
contains
 
  subroutine fstr_init_gauss( gauss )
    use muyield, only: uelastoplasticnumstatus
    type( tgaussstatus ), intent(inout) :: gauss
    integer :: n
    gauss%strain=0.d0; gauss%stress=0.d0
    gauss%strain_bak=0.d0; gauss%stress_bak=0.d0
    gauss%strain_out=0.d0; gauss%stress_out=0.d0
    gauss%plstrain =0.d0
    gauss%nqm =0.d0
    gauss%strain_energy =0.d0
    gauss%strain_energy_bak =0.d0
    if( gauss%pMaterial%mtype==usermaterial ) then
      if( gauss%pMaterial%nfstatus> 0 ) then
        allocate( gauss%fstatus(gauss%pMaterial%nfstatus) )
        gauss%fstatus(:) = 0.d0
      endif
    else if( iselastoplastic(gauss%pMaterial%mtype) ) then
      allocate( gauss%istatus(1) )    ! 0:elastic 1:plastic
      if( getyieldfunction( gauss%pMaterial%mtype )==3 ) then  ! user defined
        n = uelastoplasticnumstatus( gauss%pMaterial%variables )
        if( n>0 ) allocate( gauss%fstatus(n) )
      elseif( iskinematicharden( gauss%pMaterial%mtype ) ) then
        allocate( gauss%fstatus(7+6) )  ! plastic strain, back stress
      else
        allocate( gauss%fstatus(2) )    ! plastic strain
      endif
      gauss%istatus = 0
      gauss%fstatus = 0.d0
    else if( isviscoelastic(gauss%pMaterial%mtype) ) then
      n = fetch_tablerow( mc_viscoelastic, gauss%pMaterial%dict )
      if( n>0 ) then
        allocate( gauss%fstatus(12*n+6) )    ! visco stress components
        gauss%fstatus = 0.d0
      else
        stop "Viscoelastic properties not defined"
      endif
    else if( gauss%pMaterial%mtype==norton ) then
      allocate( gauss%fstatus(2) )        ! effective stress, effective viscoplastic strain
      gauss%fstatus = 0.d0
      gauss%plstrain = 0.d0
    endif
  end subroutine fstr_init_gauss
 
  subroutine fstr_finalize_gauss( gauss )
    type( tgaussstatus ), intent(inout) :: gauss
    if( associated( gauss%istatus ) ) deallocate( gauss%istatus )
    if( associated( gauss%fstatus ) ) deallocate( gauss%fstatus )
  end subroutine
 
  subroutine fstr_copy_gauss( gauss1, gauss2 )
    type( tgaussstatus ), intent(in)    :: gauss1
    type( tgaussstatus ), intent(inout) :: gauss2
 
    gauss2%strain     = gauss1%strain
    gauss2%stress     = gauss1%stress
    gauss2%strain_bak = gauss1%strain_bak
    gauss2%stress_bak = gauss1%stress_bak
    gauss2%plstrain   = gauss1%plstrain
    gauss2%strain_energy = gauss1%strain_energy
 
    if( associated(gauss1%istatus) .and. associated(gauss2%istatus) ) then
      gauss2%istatus   = gauss1%istatus
    end if
    if( associated(gauss1%fstatus) .and. associated(gauss2%fstatus) ) then
      gauss2%fstatus   = gauss1%fstatus
    end if
  end subroutine fstr_copy_gauss
 
 
end module