soltosph_hour.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sig_heph1 (jr0, js0, jt0, gsig, fhour, sig_heph, pm, ixs, ii, nel, lft, llt)

Function/Subroutine Documentation

sig_heph1()

subroutine sig_heph1	(		jr0,
			js0,
			jt0,
			gsig,
			fhour,
			sig_heph,
			pm,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(6)	ii,
		integer	nel,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt )

Definition at line 30 of file soltosph_hour.F.

      use element_mod , only : nixs
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER IXS(NIXS,*),II(6),NEL
      my_real
     .   gsig(*),fhour(nel,3,4),jr0(*),js0(*),jt0(*) ,
     .   sig_heph(mvsiz,6,7),pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   nu1,nu2,nu3
C-----------------------------------------------
 
C-----------------------------------------------
C
C       SIG_HEPH(,1) -> mean stress
C       SIG_HEPH(,2) -> hourglass stress component proportional to ZETA
C       SIG_HEPH(,3) -> hourglass stress component proportional to ETA 
C       SIG_HEPH(,4) -> hourglass stress component proportional to KSI 
C       SIG_HEPH(,5) -> hourglass stress component proportional to ZETA*ETA
C       SIG_HEPH(,6) -> hourglass stress component proportional to ZETA*KSI 
C       SIG_HEPH(,7) -> hourglass stress component proportional to ETA*KSI    
C
        DO i=lft,llt
           nu1=one/(one-pm(21,ixs(1,i)))
           nu2=pm(21,ixs(1,i))*nu1
           nu3=one+pm(21,ixs(1,i))
C
           sig_heph(i,1,1)=gsig(ii(1)+i)
           sig_heph(i,2,1)=gsig(ii(2)+i)
           sig_heph(i,3,1)=gsig(ii(3)+i)
           sig_heph(i,4,1)=gsig(ii(4)+i)
           sig_heph(i,5,1)=gsig(ii(5)+i)
           sig_heph(i,6,1)=gsig(ii(6)+i)
C
           sig_heph(i,1,2)=2*twenty4*(nu1*fhour(i,1,2)+nu2*fhour(i,2,1))
           sig_heph(i,2,2)=2*twenty4*(nu1*fhour(i,2,1)+nu2*fhour(i,1,2))
           sig_heph(i,4,2)=twenty4*(jr0(i)*fhour(i,1,1)/js0(i)+js0(i)*fhour(i,2,2)/jr0(i))
           sig_heph(i,1,3)=2*twenty4*(nu1*fhour(i,1,3)+nu2*fhour(i,3,1))
           sig_heph(i,3,3)=2*twenty4*(nu1*fhour(i,3,1)+nu2*fhour(i,1,3))
           sig_heph(i,6,3)=twenty4*(jr0(i)*fhour(i,1,1)/jt0(i)+jt0(i)*fhour(i,3,3)/jr0(i))
           sig_heph(i,2,4)=2*twenty4*(nu1*fhour(i,2,3)+nu2*fhour(i,3,2))
           sig_heph(i,3,4)=2*twenty4*(nu1*fhour(i,3,2)+nu2*fhour(i,2,3))
           sig_heph(i,5,4)=twenty4*(js0(i)*fhour(i,2,2)/jt0(i)+jt0(i)*fhour(i,3,3)/js0(i))
           sig_heph(i,1,5)=2*twenty4*nu3*fhour(i,1,4)        
           sig_heph(i,2,6)=2*twenty4*nu3*fhour(i,2,4)
           sig_heph(i,3,7)=2*twenty4*nu3*fhour(i,3,4)
        END DO
C
      RETURN