szstraingps.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	szstraingps (strain, str_pi, strhg, nel)

Function/Subroutine Documentation

szstraingps()

subroutine szstraingps	(		strain,
			str_pi,
			strhg,
		integer	nel )

Definition at line 30 of file szstraingps.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   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 NEL
      my_real
     .   strain(nel,6) ,str_pi(mvsiz,6,8),strhg(nel,3,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,  J,IR,IS,IT,IP
      my_real
     .   ksi,eta,zeta,coef,f1,f2,f3,f1x,f2y,f3z
      my_real
     .  a_gauss(9,9)
      DATA a_gauss / 
     1 0.               ,0.               ,0.               ,
     1 0.               ,0.               ,0.               ,
     1 0.               ,0.               ,0.               ,
     2 -.577350269189626,0.577350269189626,0.               ,
     2 0.               ,0.               ,0.               ,
     2 0.               ,0.               ,0.               ,
     3 -.774596669241483,0.               ,0.774596669241483,
     3 0.               ,0.               ,0.               ,
     3 0.               ,0.               ,0.               ,
     4 -.861136311594053,-.339981043584856,0.339981043584856,
     4 0.861136311594053,0.               ,0.               ,
     4 0.               ,0.               ,0.               ,
     5 -.906179845938664,-.538469310105683,0.               ,
     5 0.538469310105683,0.906179845938664,0.               ,
     5 0.               ,0.               ,0.               ,
     6 -.932469514203152,-.661209386466265,-.238619186083197,
     6 0.238619186083197,0.661209386466265,0.932469514203152,
     6 0.               ,0.               ,0.               ,
     7 -.949107912342759,-.741531185599394,-.405845151377397,
     7 0.               ,0.405845151377397,0.741531185599394,
     7 0.949107912342759,0.               ,0.               ,
     8 -.960289856497536,-.796666477413627,-.525532409916329,
     8 -.183434642495650,0.183434642495650,0.525532409916329,
     8 0.796666477413627,0.960289856497536,0.               ,
     9 -.968160239507626,-.836031107326636,-.613371432700590,
     9 -.324253423403809,0.               ,0.324253423403809,
     9 0.613371432700590,0.836031107326636,0.968160239507626/
      
C-----------------------------------------------
      DO ir=1,2
       DO is=1,2
        DO it=1,2
         ip = ir + ( (is-1) + (it-1)*2 )*2
         ksi = a_gauss(ir,2)
         eta = a_gauss(is,2)
         zeta = a_gauss(it,2)
         f1 = eta*zeta
         f2 = ksi*zeta
         f3 = eta*ksi
C
         DO i=1,nel
          str_pi(i,1,ip)=strain(i,1)+
     +                   strhg(i,1,1)*ksi+strhg(i,1,2)*eta+strhg(i,1,3)*zeta+
     +                   strhg(i,1,5)*f1+strhg(i,2,6)*f2+strhg(i,3,6)*f3
          str_pi(i,2,ip)=strain(i,2)+
     +                   strhg(i,2,1)*ksi+strhg(i,2,2)*eta+strhg(i,2,3)*zeta+
     +                   strhg(i,2,5)*f2+strhg(i,1,6)*f1+strhg(i,3,6)*f3
          str_pi(i,3,ip)=strain(i,3)+
     +                   strhg(i,3,1)*ksi+strhg(i,3,2)*eta+strhg(i,3,3)*zeta+
     +                   strhg(i,3,5)*f3+strhg(i,1,6)*f1+strhg(i,2,6)*f2
          str_pi(i,4,ip)=strain(i,4)+strhg(i,1,4)*zeta
          str_pi(i,6,ip)=strain(i,5)+strhg(i,2,4)*eta
          str_pi(i,5,ip)=strain(i,6)+strhg(i,3,4)*ksi
         END DO
        ENDDO 
       ENDDO 
      ENDDO
C      
      RETURN