Added ideas for Voigt<->Mandel transformations (PE)

This commit is contained in:
Luc Hantcherli 2007-03-28 17:41:54 +00:00
parent 4b69c1d738
commit 3ba46ac7d7
1 changed files with 74 additions and 54 deletions

View File

@ -31,6 +31,20 @@
(/1.0_pReal,1.0_pReal,1.0_pReal, 1.414213562373095_pReal, 1.414213562373095_pReal, 1.414213562373095_pReal/)
real(pReal), dimension(6), parameter :: invnrmMandel = &
(/1.0_pReal,1.0_pReal,1.0_pReal,0.7071067811865476_pReal,0.7071067811865476_pReal,0.7071067811865476_pReal/)
! *** Voigt notation ***
integer(pInt), dimension (2,6), parameter :: mapMandel = &
reshape((/&
1,1, &
2,2, &
3,3, &
2,3, &
1,3, &
1,2 &
/),(/2,6/))
real(pReal), dimension(6), parameter :: nrmVoigt = &
(/1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal/)
real(pReal), dimension(6), parameter :: invnrmVoigt = &
(/1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal,1.0_pReal/)
CONTAINS
@ -326,23 +340,6 @@
END SUBROUTINE Gauss
!********************************************************************
! calculate v6 stress from M6x6 stiffness and v6 strain
!********************************************************************
FUNCTION math_Hooke(C,e)
use prec, only: pReal,pInt
implicit none
real(pReal) C(6,6), e(6), factoredE(6)
real(pReal), dimension(6), parameter :: preFactor = (/1.0_pReal,1.0_pReal,1.0_pReal,2.0_pReal,2.0_pReal,2.0_pReal/)
real(pReal) math_Hooke(6)
factoredE = dot_product(preFactor,e)
math_Hooke = matmul(C,factoredE)
return
END FUNCTION
!********************************************************************
@ -444,6 +441,29 @@
END FUNCTION
!********************************************************************
! convert Voigt matrix 6x6 back to symmetric 3x3x3x3 tensor
!********************************************************************
FUNCTION math_Voigt66to3333(m66)
use prec, only: pReal,pInt
implicit none
real(pReal), dimension(6,6) :: m66
real(pReal), dimension(3,3,3,3) :: math_Voigt66to3333
integer(pInt) i,j
forall (i=1:6,j=1:6)
math_Voigt66to3333(mapVoigt(1,i),mapVoigt(2,i),mapVoigt(1,j),mapVoigt(2,j)) = invnrmVoigt(i)*invnrmVoigt(j)*m66(i,j)
math_Voigt66to3333(mapVoigt(2,i),mapVoigt(1,i),mapVoigt(1,j),mapVoigt(2,j)) = invnrmVoigt(i)*invnrmVoigt(j)*m66(i,j)
math_Voigt66to3333(mapVoigt(1,i),mapVoigt(2,i),mapVoigt(2,j),mapVoigt(1,j)) = invnrmVoigt(i)*invnrmVoigt(j)*m66(i,j)
math_Voigt66to3333(mapVoigt(2,i),mapVoigt(1,i),mapVoigt(2,j),mapVoigt(1,j)) = invnrmVoigt(i)*invnrmVoigt(j)*m66(i,j)
end forall
return
END FUNCTION
!********************************************************************