szstrainhg.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	szstrainhg (jr_1, js_1, jt_1, strhg, nel, hgx1, hgx2, hgx3, hgx4, hgy1, hgy2, hgy3, hgy4, hgz1, hgz2, hgz3, hgz4, nu, nu1)

Function/Subroutine Documentation

szstrainhg()

subroutine szstrainhg	(		jr_1,
			js_1,
			jt_1,
			strhg,
		integer	nel,
			hgx1,
			hgx2,
			hgx3,
			hgx4,
			hgy1,
			hgy2,
			hgy3,
			hgy4,
			hgz1,
			hgz2,
			hgz3,
			hgz4,
			nu,
			nu1 )

Definition at line 31 of file szstrainhg.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-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   hgx1(mvsiz), hgx2(mvsiz), hgx3(mvsiz), hgx4(mvsiz),
     .   hgy1(mvsiz), hgy2(mvsiz), hgy3(mvsiz), hgy4(mvsiz),
     .   hgz1(mvsiz), hgz2(mvsiz), hgz3(mvsiz), hgz4(mvsiz),
     .   jr_1(*),js_1(*),jt_1(*) ,strhg(nel,3,6),nu(*),nu1(*)
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
     .   hgx(mvsiz,4), hgy(mvsiz,4), hgz(mvsiz,4),
     .   ksi,eta,zeta,coef,njr_1(mvsiz),njs_1(mvsiz),njt_1(mvsiz),
     .   jsr_1(mvsiz),jrs_1(mvsiz),jtr_1(mvsiz),jrt_1(mvsiz),jts_1(mvsiz),
     .   jst_1(mvsiz),f1,f2,f3,f1x,f2y,f3z
C-----------------------------------------------
      coef = eight
      DO i=1,nel
       njr_1(i) = coef*jr_1(i)*dt1
       njs_1(i) = coef*js_1(i)*dt1
       njt_1(i) = coef*jt_1(i)*dt1
      ENDDO
      DO i=1,nel
       hgx(i,1) = njr_1(i)*hgx1(i)
       hgx(i,2) = njr_1(i)*hgx2(i)
       hgx(i,3) = njr_1(i)*hgx3(i)
       hgx(i,4) = njr_1(i)*hgx4(i)
       hgy(i,1) = njs_1(i)*hgy1(i)
       hgy(i,2) = njs_1(i)*hgy2(i)
       hgy(i,3) = njs_1(i)*hgy3(i)
       hgy(i,4) = njs_1(i)*hgy4(i)
       hgz(i,1) = njt_1(i)*hgz1(i)
       hgz(i,2) = njt_1(i)*hgz2(i)
       hgz(i,3) = njt_1(i)*hgz3(i)
       hgz(i,4) = njt_1(i)*hgz4(i)
      ENDDO
C------ STRHG(NEL,3,j) : J=1,3: term ksi,eta,zeta of Eii; J=4: shear;   
C------                  J=5: mode 4 HGi4; J=6: nu*HGi4;   
         DO i=1,nel
          f1x = hgx(i,4)
          f2y = hgy(i,4)
          f3z = hgz(i,4)
          strhg(i,1,1)= strhg(i,1,1)-nu1(i)*(hgy(i,3)+hgz(i,2))
          strhg(i,1,2)= strhg(i,1,2)+hgx(i,3)
          strhg(i,1,3)= strhg(i,1,3)+hgx(i,2)
          strhg(i,1,4)= strhg(i,1,4)+hgx1(i)*njs_1(i)+hgy2(i)*njr_1(i)
          strhg(i,1,5)= strhg(i,1,5)+f1x
          strhg(i,1,6)= strhg(i,1,6)-nu(i)*f1x
          strhg(i,2,1)= strhg(i,2,1)+hgy(i,3)
          strhg(i,2,2)= strhg(i,2,2)-nu1(i)*(hgx(i,3)+hgz(i,1))
          strhg(i,2,3)= strhg(i,2,3)+hgy(i,1)
          strhg(i,2,4)= strhg(i,2,4)+hgx1(i)*njt_1(i)+hgz3(i)*njr_1(i)
          strhg(i,2,5)= strhg(i,2,5)+f2y
          strhg(i,2,6)= strhg(i,2,6)-nu(i)*f2y
          strhg(i,3,1)= strhg(i,3,1)+hgz(i,2)
          strhg(i,3,2)= strhg(i,3,2)+hgz(i,1)
          strhg(i,3,3)= strhg(i,3,3)-nu1(i)*(hgx(i,2)+hgy(i,1))
          strhg(i,3,4)= strhg(i,3,4)+hgy2(i)*njt_1(i)+hgz3(i)*njs_1(i)
          strhg(i,3,5)= strhg(i,3,5)+f3z
          strhg(i,3,6)= strhg(i,3,6)-nu(i)*f3z
         END DO
C      
      RETURN