cdk6fint3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cdk6fint3 (jft, jlt, vol, thk0, for, mom, px2, py2, px3, py3, pb1, pb2, pb3, f11, f12, f13, f21, f22, f23, f31, f32, f33, f14, f15, f16, f24, f25, f26, f34, f35, f36, n4x, n4y, n4z, n5x, n5y, n5z, n6x, n6y, n6z, nel)

Function/Subroutine Documentation

cdk6fint3()

subroutine cdk6fint3	(	integer	jft,
		integer	jlt,
			vol,
			thk0,
			for,
			mom,
			px2,
			py2,
			px3,
			py3,
			pb1,
			pb2,
			pb3,
			f11,
			f12,
			f13,
			f21,
			f22,
			f23,
			f31,
			f32,
			f33,
			f14,
			f15,
			f16,
			f24,
			f25,
			f26,
			f34,
			f35,
			f36,
			n4x,
			n4y,
			n4z,
			n5x,
			n5y,
			n5z,
			n6x,
			n6y,
			n6z,
		integer	nel )

Definition at line 28 of file cdk6fint3.F.

C--------------------------------------------------------------------------------------------------
C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT,NEL
      my_real 
     .   for(nel,5),mom(nel,3), 
     .   px2(*), py2(*), px3(*), py3(*), 
     .   pb1(mvsiz,3,3),pb2(mvsiz,3,3),pb3(mvsiz,3,6),
     .   f11(*), f12(*), f13(*),
     .   f21(*), f22(*), f23(*), f31(*), f32(*), f33(*),
     .   n4x(*), n4y(*), n4z(*),n5x(*), n5y(*), n5z(*),
     .   n6x(*), n6y(*), n6z(*),
     .   f14(*), f15(*), f16(*),f24(*), f25(*), f26(*),
     .   f34(*), f35(*), f36(*),thk0(*),vol(*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I
      my_real 
     .    c2,fxx,fyy,fxy,fz4,fz5,fz6,mxx,myy,mxy,fxn,fyn,fxyn
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      DO i=jft,jlt
        f12(i) =vol(i)*(px2(i)*for(i,1)+py2(i)*for(i,3))
        f22(i) =vol(i)*(py2(i)*for(i,2)+px2(i)*for(i,3))
        f13(i) =vol(i)*(px3(i)*for(i,1)+py3(i)*for(i,3))
        f23(i) =vol(i)*(py3(i)*for(i,2)+px3(i)*for(i,3))
        f11(i)=-f12(i)-f13(i)
        f21(i)=-f22(i)-f23(i)
      ENDDO     
      DO i=jft,jlt
          c2=vol(i)*thk0(i)
        mxx=c2*mom(i,1)
        myy=c2*mom(i,2)
        mxy=c2*mom(i,3)
        f31(i) = pb1(i,1,1)*mxx+pb1(i,2,1)*myy+pb1(i,3,1)*mxy
        f32(i) = pb1(i,1,2)*mxx+pb1(i,2,2)*myy+pb1(i,3,2)*mxy
        f33(i) = pb1(i,1,3)*mxx+pb1(i,2,3)*myy+pb1(i,3,3)*mxy
C
        fz4    = pb2(i,1,1)*mxx+pb2(i,2,1)*myy+pb2(i,3,1)*mxy
        f14(i) = fz4*n4x(i)
        f24(i) = fz4*n4y(i)
        f34(i) = fz4*n4z(i)
        fz5    = pb2(i,1,2)*mxx+pb2(i,2,2)*myy+pb2(i,3,2)*mxy
        f15(i) = fz5*n5x(i)
        f25(i) = fz5*n5y(i)
        f35(i) = fz5*n5z(i)
        fz6    = pb2(i,1,3)*mxx+pb2(i,2,3)*myy+pb2(i,3,3)*mxy
        f16(i) = fz6*n6x(i)
        f26(i) = fz6*n6y(i)
        f36(i) = fz6*n6z(i)
C
        fxx    = pb3(i,1,1)*mxx + pb3(i,2,1)*myy + pb3(i,3,1)*mxy
        fyy    = pb3(i,1,2)*mxx + pb3(i,2,2)*myy + pb3(i,3,2)*mxy
        fxy    = pb3(i,1,3)*mxx + pb3(i,2,3)*myy + pb3(i,3,3)*mxy
        fxn  = n5x(i)*fxx
        fyn  = n6x(i)*fyy
        fxyn = n4x(i)*fxy
        f11(i) = f11(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
        f12(i) = f12(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
        f13(i) = f13(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
        fxn  = n5y(i)*fxx
        fyn  = n6y(i)*fyy
        fxyn = n4y(i)*fxy
        f21(i) = f21(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
        f22(i) = f22(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
        f23(i) = f23(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
        fxn  = n5z(i)*fxx
        fyn  = n6z(i)*fyy
        fxyn = n4z(i)*fxy
        f31(i) = f31(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
        f32(i) = f32(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
        f33(i) = f33(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
      ENDDO     
C    
      RETURN