q4fint2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4fint2 (sig, ay, fay, faz, py1, py2, py3, py4, pz1, pz2, pz3, pz4, byz1, byz2, byz3, byz4, bzy1, bzy2, bzy3, bzy4, fy1, fz1, fy2, fz2, fy3, fz3, fy4, fz4, r22, r23, r32, r33, air, vol, qvis, icp, nel, jhbe, jcvt, svis)

Function/Subroutine Documentation

q4fint2()

subroutine q4fint2	(		sig,
			ay,
			fay,
			faz,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			byz1,
			byz2,
			byz3,
			byz4,
			bzy1,
			bzy2,
			bzy3,
			bzy4,
			fy1,
			fz1,
			fy2,
			fz2,
			fy3,
			fz3,
			fy4,
			fz4,
			r22,
			r23,
			r32,
			r33,
			air,
			vol,
			qvis,
		integer	icp,
		integer	nel,
		integer, intent(in)	jhbe,
		integer, intent(in)	jcvt,
		intent(inout)	svis )

Definition at line 28 of file q4fint2.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      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JHBE
      INTEGER, INTENT(IN) :: JCVT
      INTEGER ICP,NEL
C     REAL
      my_real
     .    sig(nel,6),ay(*),fay(*),faz(*),
     .    py1(*),py2(*),py3(*),py4(*),
     .    pz1(*),pz2(*),pz3(*),pz4(*),
     .    byz1(*),byz2(*),byz3(*),byz4(*),
     .    bzy1(*),bzy2(*),bzy3(*),bzy4(*),
     .    fy1(*),fz1(*),fy2(*),fz2(*),
     .    fy3(*),fz3(*),fy4(*),fz4(*),
     .    r22(*),r23(*),r32(*),r33(*),
     .    air(*),vol(*),qvis(*)
      my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
C-----------------------------------------------
C FUNCTION: 
C ARGUMENTS:  (I: input, O: output, IO: input & output, W: workspace)
C TYPE NAME                FUNCTION
C  I   SIG(NEL,6)           - STRESS COMPONENTS
C  I   AY(*)              - Ni/r AT CENTER FOR AXISYMMETRIC CASE
C  I   FAY(*),FAZ(*)      - COMPONENT OF INTERNAL FORCE FOR AXISYMMETRIC CASE
C  I   PY1(*)~PZ4(*)      - SHAPE DERIVATIVES (dNi/dY, dNi/dZ)
C  O   FY1(*)~FZ4(*)      - INTERNAL FORCE IN LOCAL STORAGE
C  I   R22(*)~R33(*)      - TRANSFORMATION MATRIX
C  I   AIR(*)             - W*|J|
C  I   VOL(*)             - W*|J| FOR PLAIN CASE; r'*W*|J| FOR AXISYMMETRIC CASE
C  I   QVIS(*)            - VISCOUS PRESSURE
C  I   ICP                - FLAG FOR CONSTANT PRESURE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .    s1(mvsiz),s2(mvsiz),s3(mvsiz),s4(mvsiz),
     .    fac(mvsiz),p(mvsiz),t1,t2,t3,t4
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      IF(n2d==1 .AND. jhbe==17) THEN
        DO i=1,nel
          fac(i) = air(i)
        ENDDO
      ELSE
        DO i=1,nel
          fac(i) = vol(i)
        ENDDO
      ENDIF
C
C     ADD VISCOUS STRESSES "SVIS" & "QVIS" INTO THE TOTAL ONES
C     REMOVE PRESSURE FROM STRESSES FOR THE CONSTANT PRESSURE CASE
      IF(icp==1) THEN
        DO i=1,nel
          p(i)=third*(sig(i,1)+sig(i,2)+sig(i,3)+
     .                svis(i,1)+svis(i,2)+svis(i,3))
          s1(i)=(sig(i,1)+svis(i,1)-p(i))*fac(i)
          s2(i)=(sig(i,2)+svis(i,2)-p(i))*fac(i)
        ENDDO
      ELSE
        DO i=1,nel
          s1(i)=(sig(i,1)+svis(i,1)-qvis(i))*fac(i)
          s2(i)=(sig(i,2)+svis(i,2)-qvis(i))*fac(i)
        ENDDO
      ENDIF
C
      DO i=1,nel
        fy1(i) = fy1(i) - s1(i)*py1(i)- s2(i)*bzy1(i)
        fy2(i) = fy2(i) - s1(i)*py2(i)- s2(i)*bzy2(i)
        fy3(i) = fy3(i) - s1(i)*py3(i)- s2(i)*bzy3(i)
        fy4(i) = fy4(i) - s1(i)*py4(i)- s2(i)*bzy4(i)
        fz1(i) = fz1(i) - s2(i)*pz1(i)- s1(i)*byz1(i)
        fz2(i) = fz2(i) - s2(i)*pz2(i)- s1(i)*byz2(i)
        fz3(i) = fz3(i) - s2(i)*pz3(i)- s1(i)*byz3(i)
        fz4(i) = fz4(i) - s2(i)*pz4(i)- s1(i)*byz4(i)
      ENDDO
C
      IF(n2d==1.AND.jhbe==17) THEN
         IF(icp==1) THEN
          DO i=1,nel
            s3(i)=(sig(i,3)+svis(i,3)-p(i))*fac(i)
            s4(i)=(sig(i,4)+svis(i,4))*fac(i)
          ENDDO
           ELSE
          DO i=1,nel
            s3(i)=(sig(i,3)+svis(i,3)-qvis(i))*fac(i)
            s4(i)=(sig(i,4)+svis(i,4))*fac(i)
          ENDDO
           END IF
         IF(jcvt/=0) THEN
            DO i=1,nel
              t1=s1(i)*r22(i)+s4(i)*r23(i)
              t2=s4(i)*r32(i)+s2(i)*r33(i)
              t3=s1(i)*r32(i)+s4(i)*r33(i)
              t4=s4(i)*r22(i)+s2(i)*r23(i)
              s1(i)=r22(i)*t1+r23(i)*t4
              s2(i)=r32(i)*t3+r33(i)*t2
              s4(i)=r22(i)*t3+r23(i)*t2
            ENDDO
          ENDIF
          DO i=1,nel
            fay(i) = fay(i) + (s1(i)-s3(i))*ay(i)
            faz(i) = faz(i) + s4(i)*ay(i)
          ENDDO
      ENDIF
C
      RETURN