czstrah3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	czstrah3 (jft, jlt, stra_h, vhg, a_i, mx23, mx34, my23, my34, dt1c, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, nel)

Function/Subroutine Documentation

czstrah3()

subroutine czstrah3	(	integer, intent(in)	jft,
		integer, intent(in)	jlt,
		intent(inout)	stra_h,
		intent(in)	vhg,
		intent(in)	a_i,
		intent(in)	mx23,
		intent(in)	mx34,
		intent(in)	my23,
		intent(in)	my34,
		intent(in)	dt1c,
		intent(in)	exx,
		intent(in)	eyy,
		intent(in)	exy,
		intent(in)	exz,
		intent(in)	eyz,
		intent(in)	kxx,
		intent(in)	kyy,
		intent(in)	kxy,
		integer, intent(in)	nel )

Definition at line 28 of file czstrah3.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, INTENT(IN) :: JFT, JLT,NEL
      my_real, DIMENSION(MVSIZ),INTENT(IN) :: a_i,
     .                    dt1c,mx23,mx34,my23,my34,
     .                    exx,eyy,exy,exz,eyz,  
     .                    kxx ,kyy ,kxy 
      my_real, DIMENSION(MVSIZ,6),INTENT(IN) :: vhg
      my_real, DIMENSION(NEL,8,4),INTENT(INOUT) :: stra_h
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER PT,PT00,PT0,I,J,EP,IUN,NG,K,NPG
      my_real 
     .   hx(mvsiz),hy(mvsiz),dhg(mvsiz,6),pg1,pg,vpg(2,4)
      my_real 
     .   a4_1(mvsiz),a_4(mvsiz),strain(8),hxz(mvsiz),hyz(mvsiz)
       parameter(pg=.577350269189626)
       parameter(pg1=-.577350269189626)
      DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/  ! VPG(1:2,4):ksi,eta
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      npg =4
      a4_1(jft:jlt) = four*a_i(jft:jlt)
      a_4(jft:jlt) = fourth*a_i(jft:jlt)
      DO j=1,6
        dhg(jft:jlt,j) = vhg(jft:jlt,j) * dt1c(jft:jlt)
      ENDDO
      DO ng = 1,npg
        DO i=jft,jlt          
          hx(i) = a4_1(i)*(my34(i)*vpg(2,ng)- my23(i)*vpg(1,ng))
          hy(i) = a4_1(i)*(-mx34(i)*vpg(2,ng)+mx23(i)*vpg(1,ng))
          hxz(i)= a_4(i)*(my34(i)*vpg(2,ng)*dhg(i,5)- my23(i)*vpg(1,ng)*dhg(i,6))
          hyz(i)= a_4(i)*(mx23(i)*vpg(1,ng)*dhg(i,5)- mx34(i)*vpg(2,ng)*dhg(i,6))
          strain(1) = exx(i) + hx(i)*dhg(i,1)
          strain(2) = eyy(i) + hy(i)*dhg(i,2)
          strain(3) = exy(i)  
          strain(6) = kxx(i) + hx(i)*dhg(i,3)
          strain(7) = kyy(i) + hy(i)*dhg(i,4)
          strain(8) = kxy(i)  
          strain(5) = exz(i) + hxz(i)
          strain(4) = eyz(i) + hyz(i)
          stra_h(i,1,ng) = stra_h(i,1,ng) + strain(1) 
          stra_h(i,2,ng) = stra_h(i,2,ng) + strain(2) 
          stra_h(i,3,ng) = stra_h(i,3,ng) + strain(3) 
          stra_h(i,4,ng) = stra_h(i,4,ng) + strain(4) 
          stra_h(i,5,ng) = stra_h(i,5,ng) + strain(5) 
          stra_h(i,6,ng) = stra_h(i,6,ng) + strain(6) 
          stra_h(i,7,ng) = stra_h(i,7,ng) + strain(7) 
          stra_h(i,8,ng) = stra_h(i,8,ng) + strain(8) 
        ENDDO
      ENDDO
C
      RETURN