s4fint3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s4fint3 (sig, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, vol, qvis, nel, svis)

Function/Subroutine Documentation

s4fint3()

subroutine s4fint3	(		sig,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			f11,
			f21,
			f31,
			f12,
			f22,
			f32,
			f13,
			f23,
			f33,
			f14,
			f24,
			f34,
			vol,
			qvis,
		integer	nel,
		intent(inout)	svis )

Definition at line 28 of file s4fint3.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
C     REAL
      my_real
     .   sig(nel,6),
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   f11(*),f21(*),f31(*),f12(*),f22(*),f32(*),
     .   f13(*),f23(*),f33(*),f14(*),f24(*),f34(*),
     .   vol(*),qvis(*)
      my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   s1, s2, s3,
     .   s4, s5, s6
C-----------------------------------------------
      DO i=1,nel
        s1=(sig(i,1)+svis(i,1)-qvis(i))*vol(i)
        s2=(sig(i,2)+svis(i,2)-qvis(i))*vol(i)
        s3=(sig(i,3)+svis(i,3)-qvis(i))*vol(i)
        s4=(sig(i,4)+svis(i,4))*vol(i)
        s5=(sig(i,5)+svis(i,5))*vol(i)
        s6=(sig(i,6)+svis(i,6))*vol(i)
C
        f11(i)=f11(i)-(s1*px1(i)+s4*py1(i)+s6*pz1(i))
        f21(i)=f21(i)-(s2*py1(i)+s4*px1(i)+s5*pz1(i))
        f31(i)=f31(i)-(s3*pz1(i)+s6*px1(i)+s5*py1(i))
        f12(i)=f12(i)-(s1*px2(i)+s4*py2(i)+s6*pz2(i))
        f22(i)=f22(i)-(s2*py2(i)+s4*px2(i)+s5*pz2(i))
        f32(i)=f32(i)-(s3*pz2(i)+s6*px2(i)+s5*py2(i))
        f13(i)=f13(i)-(s1*px3(i)+s4*py3(i)+s6*pz3(i))
        f23(i)=f23(i)-(s2*py3(i)+s4*px3(i)+s5*pz3(i))
        f33(i)=f33(i)-(s3*pz3(i)+s6*px3(i)+s5*py3(i))
        f14(i)=f14(i)-(s1*px4(i)+s4*py4(i)+s6*pz4(i))
        f24(i)=f24(i)-(s2*py4(i)+s4*px4(i)+s5*pz4(i))
        f34(i)=f34(i)-(s3*pz4(i)+s6*px4(i)+s5*py4(i))
      ENDDO
C
      RETURN