qfint2.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	qfint2 (sig, py1, py2, pz1, pz2, area, vol, qvis, f11, f12, f21, f22, ax1, ax2, r22, r23, r32, r33, nel, jcvt, svis)

Function/Subroutine Documentation

qfint2()

subroutine qfint2	(		sig,
			py1,
			py2,
			pz1,
			pz2,
			area,
			vol,
			qvis,
			f11,
			f12,
			f21,
			f22,
			ax1,
			ax2,
			r22,
			r23,
			r32,
			r33,
		integer	nel,
		integer, intent(in)	jcvt,
		intent(inout)	svis )

Definition at line 29 of file qfint2.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, INTENT(IN) :: JCVT
      INTEGER NEL
      my_real
     .   sig(nel,6),
     .   py1(*), py2(*), pz1(*), pz2(*), area(*), vol(*), qvis(*),
     .   f11(*), f12(*), f21(*), f22(*), ax1(*),ax2(*),
     .   r22(*), r23(*), r32(*), r33(*)
      my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   fac(mvsiz), s1(mvsiz), s2(mvsiz), s4(mvsiz),
     .   t1,t2,t3,t4
C-----------------------------------------------
      IF (n2d==1) THEN
      DO i=1,nel
        fac(i)=area(i)*area(i)/ max(em15,vol(i))/four
      ENDDO
      ELSE
      DO i=1,nel
        fac(i)=zero
      ENDDO
      ENDIF
 
      IF(jcvt==0)THEN
       DO i=1,nel
         ax1(i)=(sig(i,3)+svis(i,3)-sig(i,1)-svis(i,1))*fac(i)
         ax2(i)=(sig(i,4)+svis(i,4))*fac(i)
       END DO
      ELSE
       DO i=1,nel
         s1(i)=sig(i,1)+svis(i,1)
         s2(i)=sig(i,2)+svis(i,2)
         s4(i)=sig(i,4)+svis(i,4)
         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
       END DO
       DO i=1,nel
         ax1(i)= (sig(i,3)+svis(i,3)-s1(i))*fac(i)
         ax2(i)= s4(i)*fac(i)
       END DO
      END IF
 
      DO i=1,nel
        s1(i) = sig(i,1)+svis(i,1)-qvis(i)
        s2(i) = sig(i,2)+svis(i,2)-qvis(i)
        s4(i) = sig(i,4)+svis(i,4)
        f11(i)=s1(i)*py1(i)+s4(i)*pz1(i)
        f21(i)=s2(i)*pz1(i)+s4(i)*py1(i)
        f12(i)=s1(i)*py2(i)+s4(i)*pz2(i)
        f22(i)=s2(i)*pz2(i)+s4(i)*py2(i)
      ENDDO
 
      RETURN