szsigpara.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	szsigpara (jr0, js0, jt0, fhour, sig0, sig, ksi, eta, zeta, nu, nel, i)

Function/Subroutine Documentation

szsigpara()

subroutine szsigpara	(		jr0,
			js0,
			jt0,
			fhour,
			sig0,
			sig,
			ksi,
			eta,
			zeta,
			nu,
		integer	nel,
		integer	i )

Definition at line 30 of file szsigpara.F.

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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL,I
      my_real
     .   fhour(nel,3,4),jr0(*),js0(*),jt0(*) ,
     .   sig(mvsiz,6) ,sig0(nel,6),nu(*),ksi  ,eta  ,zeta
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   coef,nu1,nu2,nu3,sxx,syy,szz,sxy,sxz,syz
C---JR0=-X1+X2+X3-X4-X5+X6+X7-X8 (local coordinates)     
C---JS0=-Y1-Y2+Y3+Y4-Y5-Y6+Y7+Y8     
C---JT0=-Z1-Z2-Z3-Z4+Z5+Z6+Z7+Z8
C---NU : Poisson's coefficient     
C-----------------------------------------------
        coef=twenty4
c        DO I=LFT,LLT
         nu1=one/(one-nu(i))
         nu2=nu(i)*nu1
         nu3=one+nu(i)
         sxx=nu1*(fhour(i,1,2)*zeta+fhour(i,1,3)*eta)+
     .       nu2*(fhour(i,2,1)*zeta+fhour(i,3,1)*eta)+         
     .       nu3*fhour(i,1,4)*zeta*eta         
         syy=nu1*(fhour(i,2,1)*zeta+fhour(i,2,3)*ksi)+
     .       nu2*(fhour(i,1,2)*zeta+fhour(i,3,2)*ksi)+         
     .       nu3*fhour(i,2,4)*zeta*ksi         
         szz=nu1*(fhour(i,3,1)*eta+fhour(i,3,2)*ksi)+
     .       nu2*(fhour(i,1,3)*eta+fhour(i,2,3)*ksi)+         
     .       nu3*fhour(i,3,4)*eta*ksi         
         sig(i,1)=sig0(i,1)+two*sxx*twenty4
         sig(i,2)=sig0(i,2)+two*syy*twenty4
         sig(i,3)=sig0(i,3)+two*szz*twenty4
c        ENDDO
C        
c        DO I=LFT,LLT
         sxy=zeta*(jr0(i)*fhour(i,1,1)/js0(i)+
     .             js0(i)*fhour(i,2,2)/jr0(i))         
         sxz=eta*(jr0(i)*fhour(i,1,1)/jt0(i)+
     .            jt0(i)*fhour(i,3,3)/jr0(i))
         syz=ksi*(js0(i)*fhour(i,2,2)/jt0(i)+
     .            jt0(i)*fhour(i,3,3)/js0(i))
         sig(i,4)=sig0(i,4)+sxy*twenty4         
         sig(i,5)=sig0(i,5)+syz*twenty4
         sig(i,6)=sig0(i,6)+sxz*twenty4
c        ENDDO
C
      RETURN