223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

added math_Hooke to calc stress from stiffness and strain

This commit is contained in:
Philip Eisenlohr 2007-03-27 15:57:18 +00:00
parent 4bd3bbd2f5
commit 31d5ccf16d
1 changed files with 22 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
trunk/math.f90
View File

@ -326,6 +326,25 @@
END SUBROUTINE END SUBROUTINE
!********************************************************************
! 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
!******************************************************************** !********************************************************************
! calculate the determinant of a (3x3) ! calculate the determinant of a (3x3)
!******************************************************************** !********************************************************************
@ -345,7 +364,7 @@
!******************************************************************** !********************************************************************
! covert a symmetric 3,3 matrix into an array of 6 ! convert a symmetric 3,3 matrix into an array of 6
!******************************************************************** !********************************************************************
FUNCTION math_33to6(m33) FUNCTION math_33to6(m33)
@ -836,7 +855,7 @@
integer(pInt) i,dimen integer(pInt) i,dimen
real(pReal) math_identity2nd(dimen,dimen) real(pReal) math_identity2nd(dimen,dimen)
math_identity = 0.0_pReal math_identity2nd = 0.0_pReal
forall (i=1:dimen) math_identity2nd(i,i) = 1.0_pReal forall (i=1:dimen) math_identity2nd(i,i) = 1.0_pReal
return return
@ -850,7 +869,7 @@
use prec, only: pReal, pInt use prec, only: pReal, pInt
implicit none implicit none
integer(pInt) i,dimen integer(pInt) i,j,k,l,dimen
real(pReal) math_identity4th(dimen,dimen,dimen,dimen) real(pReal) math_identity4th(dimen,dimen,dimen,dimen)
forall (i=1:dimen,j=1:dimen,k=1:dimen,l=1:dimen) math_identity4th(i,j,k,l) = & forall (i=1:dimen,j=1:dimen,k=1:dimen,l=1:dimen) math_identity4th(i,j,k,l) = &