uelastic.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 muelastic
  implicit none
contains
 
  subroutine uelasticmatrix( matl, strain, D )
    use hecmw
    real(kind=kreal), intent(in)  :: matl(:)   
    real(kind=kreal), intent(in)  :: strain(6) 
    real(kind=kreal), intent(out) :: d(6,6)    
 
    ! following examples of linear elasticicty
    real(kind=kreal) :: ee, pp
 
    d(:,:)=0.d0
 
    ee = matl(1)
    pp = matl(2)
    d(1,1)=ee*(1.d0-pp)/(1.d0-2.d0*pp)/(1.d0+pp)
    d(1,2)=ee*pp/(1.d0-2.d0*pp)/(1.d0+pp)
    d(1,3)=d(1,2)
    d(2,1)=d(1,2)
    d(2,2)=d(1,1)
    d(2,3)=d(1,2)
    d(3,1)=d(1,3)
    d(3,2)=d(2,3)
    d(3,3)=d(1,1)
    d(4,4)=ee/(1.d0+pp)*0.5d0
    d(5,5)=ee/(1.d0+pp)*0.5d0
    d(6,6)=ee/(1.d0+pp)*0.5d0
  end subroutine
 
  subroutine uelasticupdate( matl, strain, stress )
    use hecmw
    real(kind=kreal), intent(in)  :: matl(:)   
    real(kind=kreal), intent(in)  :: strain(6) 
    real(kind=kreal), intent(out) :: stress(6) 
 
    ! following examples of linear elasticicty
    real(kind=kreal) :: d(6,6)
    call uelasticmatrix( matl(:), strain, d )
    stress = matmul( d, strain )
  end subroutine
 
end module