rcum33.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	rcum33 (jft, jlt, xl, nc1, nc2, f, forx, fory, forz, xm, xmom, ymom, zmom, sti, stir, stifn, stifr, fx1, fx2, fy1, fy2, fz1, fz2, mx1, mx2, my1, my2, mz1, mz2, nuvar, uvar)

Function/Subroutine Documentation

rcum33()

subroutine rcum33	(	integer	jft,
		integer	jlt,
		double precision, dimension(mvsiz,3)	xl,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
			f,
			forx,
			fory,
			forz,
			xm,
			xmom,
			ymom,
			zmom,
			sti,
			stir,
			stifn,
			stifr,
			fx1,
			fx2,
			fy1,
			fy2,
			fz1,
			fz2,
			mx1,
			mx2,
			my1,
			my2,
			mz1,
			mz2,
		integer	nuvar,
			uvar )

Definition at line 28 of file rcum33.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 JFT, JLT,NUVAR,NC1(*),NC2(*)
C     REAL
      my_real
     .   f(3,*), forx(*), fory(*), forz(*), xm(3,*), xmom(*),
     .   ymom(*), zmom(*), sti(3,*),stir(3,*),stifn(*),stifr(*),
     .   fx1(*) ,fx2(*)  ,fy1(*),  fy2(*) ,fz1(*)  ,
     .   fz2(*), mx1(*) ,mx2(*)  ,my1(*),  my2(*) ,mz1(*)  ,
     .   mz2(*),uvar(nuvar,*)
      DOUBLE PRECISION XL(MVSIZ,3)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real f11(mvsiz), f21(mvsiz), f31(mvsiz),
     .        m11(mvsiz), m21(mvsiz), m31(mvsiz),
     .        m12(mvsiz), m22(mvsiz), m32(mvsiz),
     .        xmom1, xmom2, ymom1, ymom2, zmom1, zmom2, m1,m2,m3,
     .        ex(9), xx
C=======================================================================
C-------------------------
C     FORCES
C-------------------------
      DO i=jft,jlt
        DO j=1,9
          ex(j) = uvar(21+j,i)
        ENDDO
        f11(i)=ex(1)*forx(i)+ex(4)*fory(i)+ex(7)*forz(i)
        f21(i)=ex(2)*forx(i)+ex(5)*fory(i)+ex(8)*forz(i)
        f31(i)=ex(3)*forx(i)+ex(6)*fory(i)+ex(9)*forz(i)
      ENDDO
   
C----  Assembly
      DO i=jft,jlt
        f(1,nc1(i))=f(1,nc1(i))+f11(i)
        f(2,nc1(i))=f(2,nc1(i))+f21(i)
        f(3,nc1(i))=f(3,nc1(i))+f31(i)
        fx1(i) = -f11(i)
        fy1(i) = -f21(i)
        fz1(i) = -f31(i)
        stifn(nc1(i))=stifn(nc1(i))+two*sti(1,i)
        f(1,nc2(i))=f(1,nc2(i))-f11(i)
        f(2,nc2(i))=f(2,nc2(i))-f21(i)
        f(3,nc2(i))=f(3,nc2(i))-f31(i)
        fx2(i) = f11(i)
        fy2(i) = f21(i)
        fz2(i) = f31(i)
        stifn(nc2(i))=stifn(nc2(i))+two*sti(2,i)
      ENDDO
C
C-------------------------
C     MOMENTS
C-------------------------
C
      DO i=jft,jlt
C      
        DO j=1,9
          ex(j) = uvar(21+j,i)
        END DO
        m1 = half*(xl(i,2)*forz(i)-xl(i,3)*fory(i))
        m2 = half*(xl(i,3)*forx(i)-xl(i,1)*forz(i))
        m3 = half*(xl(i,1)*fory(i)-xl(i,2)*forx(i))
        xmom1 = xmom(i) + m1
        ymom1 = ymom(i) + m2
        zmom1 = zmom(i) + m3
        m11(i)=ex(1)*xmom1+ex(4)*ymom1+ex(7)*zmom1
        m21(i)=ex(2)*xmom1+ex(5)*ymom1+ex(8)*zmom1
        m31(i)=ex(3)*xmom1+ex(6)*ymom1+ex(9)*zmom1
C
        xmom2 = xmom(i) - m1
        ymom2 = ymom(i) - m2
        zmom2 = zmom(i) - m3
        m12(i)=ex(1)*xmom2+ex(4)*ymom2+ex(7)*zmom2
        m22(i)=ex(2)*xmom2+ex(5)*ymom2+ex(8)*zmom2
        m32(i)=ex(3)*xmom2+ex(6)*ymom2+ex(9)*zmom2
C
        xx = xl(i,1)*xl(i,1) + xl(i,2)*xl(i,2) + xl(i,3)*xl(i,3)
        stir(1,i) = stir(1,i) + sti(2,i)*xx
        stir(2,i) = stir(2,i) + sti(1,i)*xx
      ENDDO
 
C----   Assembly
      DO i=jft,jlt
        xm(1,nc1(i))=xm(1,nc1(i))+m11(i)
        xm(2,nc1(i))=xm(2,nc1(i))+m21(i)
        xm(3,nc1(i))=xm(3,nc1(i))+m31(i)
        mx1(i) = -m11(i)
        my1(i) = -m21(i)
        mz1(i) = -m31(i)
        stifr(nc1(i))=stifr(nc1(i))+stir(1,i)
        xm(1,nc2(i))=xm(1,nc2(i))-m12(i)
        xm(2,nc2(i))=xm(2,nc2(i))-m22(i)
        xm(3,nc2(i))=xm(3,nc2(i))-m32(i)
        mx2(i) = m12(i)
        my2(i) = m22(i)
        mz2(i) = m32(i)
        stifr(nc2(i))=stifr(nc2(i))+stir(2,i)
      ENDDO
C-----------------------------------------------
      RETURN