q4vis2.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	q4vis2 (pm, off, rho, y1, y2, y3, y4, z1, z2, z3, z4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, py1, py2, pz1, pz2, t11, t12, t13, t14, t21, t22, t23, t24, area, cxx, nel)

Function/Subroutine Documentation

q4vis2()

subroutine q4vis2	(		pm,
			off,
			rho,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			py1,
			py2,
			pz1,
			pz2,
			t11,
			t12,
			t13,
			t14,
			t21,
			t22,
			t23,
			t24,
			area,
			cxx,
		integer, intent(in)	nel )

Definition at line 28 of file q4vis2.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   pm(npropm,*),off(*), rho(*),
     .   y1(*),y2(*),y3(*),y4(*),z1(*),z2(*),z3(*),z4(*),
     .   vy1(*), vy2(*), vy3(*), vy4(*), vz1(*), vz2(*), vz3(*),
     .   py1(*), py2(*), pz1(*), pz2(*),
     .   t11(*), t12(*), t13(*), t14(*), t21(*), t22(*), t23(*), t24(*),
     .   vz4(*), area(*), cxx(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
      my_real
     .   caq(mvsiz),fcl(mvsiz) , 
     .   g11(mvsiz) , g21(mvsiz) , g31(mvsiz) , g41(mvsiz),
     .   hgy(mvsiz), hgz(mvsiz), 
     .   hy,hz,fac,px1h1,px2h1, are
C-----------------------------------------------
C 
c      will use CAQ if necessary CAQ=0.001 (Quad BELYTSCHKO CAQ=0.1)
C
      DO i=1,nel
       fac = em03*off(i)
       are = max(area(i),em20)
       fcl(i)=half*fac*rho(i)*sqrt(are)*cxx(i)
      END DO 
C
C
C HOURGLASS BELYTSCHKO
C
        DO i=1,nel
         hy=y1(i)-y2(i)+y3(i)-y4(i)
         hz=z1(i)-z2(i)+z3(i)-z4(i)
         fac=one/max(em20,area(i))
         px1h1=fac*(py1(i)*hy+pz1(i)*hz)
         px2h1=fac*(py2(i)*hy+pz2(i)*hz)
         g11(i)= one -px1h1
         g21(i)=-one -px2h1
         g31(i)= one +px1h1
         g41(i)=-one +px2h1 
        END DO 
        DO i=1,nel
         hgy(i)=fcl(i)*(
     &    g11(i)*vy1(i)+g21(i)*vy2(i)+g31(i)*vy3(i)+g41(i)*vy4(i))
         hgz(i)=fcl(i)*(
     &    g11(i)*vz1(i)+g21(i)*vz2(i)+g31(i)*vz3(i)+g41(i)*vz4(i))
        END DO 
        DO i=1,nel
         t11(i) =t11(i) - g11(i)*hgy(i)
         t12(i) =t12(i) - g21(i)*hgy(i)
         t13(i) =t13(i) - g31(i)*hgy(i)
         t14(i) =t14(i) - g41(i)*hgy(i)
         t21(i) =t21(i) - g11(i)*hgz(i)
         t22(i) =t22(i) - g21(i)*hgz(i)
         t23(i) =t23(i) - g31(i)*hgz(i)
         t24(i) =t24(i) - g41(i)*hgz(i)
        END DO 
C
      RETURN