ccurv3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	ccurv3 (jft, jlt, vrl1, vrl2, vrl3, vrl4, px1, px2, py1, py2, area, rx1, rx2, rx3, rx4, ry1, ry2, ry3, ry4, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, kxx, kyy, kxy, exz, eyz, wxy, ismstr)

Function/Subroutine Documentation

ccurv3()

subroutine ccurv3	(	integer	jft,
		integer	jlt,
			vrl1,
			vrl2,
			vrl3,
			vrl4,
			px1,
			px2,
			py1,
			py2,
			area,
			rx1,
			rx2,
			rx3,
			rx4,
			ry1,
			ry2,
			ry3,
			ry4,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
			kxx,
			kyy,
			kxy,
			exz,
			eyz,
			wxy,
		integer	ismstr )

Definition at line 29 of file ccurv3.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, ISMSTR
      my_real
     .     rx1(mvsiz),  rx2(mvsiz),  rx3(mvsiz),  rx4(mvsiz),
     .     ry1(mvsiz),  ry2(mvsiz),  ry3(mvsiz),  ry4(mvsiz),
     .     px1(mvsiz),  px2(mvsiz),  py1(mvsiz),  py2(mvsiz),
     .     e1x(mvsiz),  e1y(mvsiz),  e1z(mvsiz),  e2x(mvsiz),
     .     e2y(mvsiz),  e2z(mvsiz),  e3x(mvsiz),  e3y(mvsiz),
     .     e3z(mvsiz),
     .     vrl1(mvsiz,3), vrl2(mvsiz,3), vrl3(mvsiz,3), vrl4(mvsiz,3),
     .     kxx(mvsiz),  kyy(mvsiz),
     .     kxy(mvsiz),  area(mvsiz), exz(mvsiz),  eyz(mvsiz),wxy(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real
     .     rxavt(mvsiz), rx13t(mvsiz), rx24t(mvsiz), ry13t(mvsiz),
     .     ryavt(mvsiz), ry24t(mvsiz),
     .     rxavta,rx13ta,rx24ta,ry13ta,ryavta,ry24ta
C=======================================================================
      DO i=jft,jlt
      rx1(i)=e1x(i)*vrl1(i,1)+e1y(i)*vrl1(i,2)+e1z(i)*vrl1(i,3)
      rx2(i)=e1x(i)*vrl2(i,1)+e1y(i)*vrl2(i,2)+e1z(i)*vrl2(i,3)
      rx3(i)=e1x(i)*vrl3(i,1)+e1y(i)*vrl3(i,2)+e1z(i)*vrl3(i,3)
      rx4(i)=e1x(i)*vrl4(i,1)+e1y(i)*vrl4(i,2)+e1z(i)*vrl4(i,3)
      ry1(i)=e2x(i)*vrl1(i,1)+e2y(i)*vrl1(i,2)+e2z(i)*vrl1(i,3)
      ry2(i)=e2x(i)*vrl2(i,1)+e2y(i)*vrl2(i,2)+e2z(i)*vrl2(i,3)
      ry3(i)=e2x(i)*vrl3(i,1)+e2y(i)*vrl3(i,2)+e2z(i)*vrl3(i,3)
      ry4(i)=e2x(i)*vrl4(i,1)+e2y(i)*vrl4(i,2)+e2z(i)*vrl4(i,3)
      ENDDO
C
      DO i=jft,jlt
        rx13ta =rx1(i)-rx3(i)
        rxavta =rx1(i)+rx2(i)+rx3(i)+rx4(i)
        rx24ta =rx2(i)-rx4(i)
C
        kyy(i)=-py1(i)*rx13ta-py2(i)*rx24ta
        kxy(i)= px1(i)*rx13ta+px2(i)*rx24ta
C
        ry13ta = ry1(i)-ry3(i)
        ryavta = ry1(i)+ry2(i)+ry3(i)+ry4(i)
        ry24ta = ry2(i)-ry4(i)
C
        kxx(i)= px1(i)*ry13ta+px2(i)*ry24ta
        kxy(i)= py1(i)*ry13ta+py2(i)*ry24ta
     +       -kxy(i)
C
        exz(i)=exz(i)+ryavta*(.25*area(i))
        eyz(i)=eyz(i)-rxavta*(.25*area(i))
      ENDDO
C
C      
       IF (ismstr==10) THEN
        DO i=jft,jlt
        rx13ta =rx1(i)-rx3(i)
        rx24ta =rx2(i)-rx4(i)
C
        ry13ta = ry1(i)-ry3(i)
        ry24ta = ry2(i)-ry4(i)
C
         wxy(i) = py1(i)*ry13ta+py2(i)*ry24ta + 
     +            px1(i)*rx13ta+px2(i)*rx24ta
        ENDDO 
       ENDIF
c-----------
      RETURN