czstra3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	czstra3 (jft, jlt, nft, vdef, gstr, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, dt1c, epsdot, iepsdot, istrain, ux1, ux2, ux3, ux4, uy1, uy2, uy3, uy4, px1, px2, py1, py2, area, ismstr, mtn, wxy, f_def, gstrw, nel)

Function/Subroutine Documentation

czstra3()

subroutine czstra3	(	integer	jft,
		integer	jlt,
		integer	nft,
			vdef,
			gstr,
			exx,
			eyy,
			exy,
			exz,
			eyz,
			kxx,
			kyy,
			kxy,
			dt1c,
			epsdot,
		integer	iepsdot,
		integer	istrain,
			ux1,
			ux2,
			ux3,
			ux4,
			uy1,
			uy2,
			uy3,
			uy4,
			px1,
			px2,
			py1,
			py2,
			area,
		integer	ismstr,
		integer	mtn,
			wxy,
			f_def,
			gstrw,
		integer	nel )

Definition at line 29 of file czstra3.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, NFT, IEPSDOT, ISTRAIN, ISMSTR, MTN,NEL
      my_real
     .   vdef(mvsiz,8), gstr(nel,8),
     .   exx(mvsiz),eyy(mvsiz),exy(mvsiz),exz(mvsiz),eyz(mvsiz),
     .   kxx(mvsiz),kyy(mvsiz),kxy(mvsiz),dt1c(mvsiz),epsdot(6,*),
     .   ux1(*),ux2(*),ux3(*),ux4(*),uy1(*),uy2(*),uy3(*),uy4(*),
     .   px1(*),px2(*),py1(*),py2(*),area(*),f_def(mvsiz,8),wxy(*),gstrw(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real
     .        fac1(mvsiz),ux13,ux24,uy13,uy24,
     .        exxt,eyyt,exyt
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C
      IF(iepsdot/=0)THEN
       DO i=jft,jlt
        j = i + nft
        epsdot(1,j) = vdef(i,1)
        epsdot(2,j) = vdef(i,2)
        epsdot(3,j) = vdef(i,3)
        epsdot(4,j) = vdef(i,6)
        epsdot(5,j) = vdef(i,7)
        epsdot(6,j) = vdef(i,8)
       ENDDO
      ENDIF
C
      DO i=jft,jlt
        exx(i) = vdef(i,1) * dt1c(i)
        eyy(i) = vdef(i,2) * dt1c(i)
        exy(i) = vdef(i,3) * dt1c(i)
        eyz(i) = vdef(i,5) * dt1c(i)
        exz(i) = vdef(i,4) * dt1c(i)
        kxx(i) = vdef(i,6) * dt1c(i)
        kyy(i) = vdef(i,7) * dt1c(i)
        kxy(i) = vdef(i,8) * dt1c(i)
      ENDDO
C
      IF (istrain /= 0.OR.ismstr == 10) THEN
       IF(ismstr == 10)THEN
        DO i=jft,jlt
         gstr(i,1)=gstr(i,1)+exx(i)
         gstr(i,2)=gstr(i,2)+eyy(i)
         gstr(i,3)=gstr(i,3)+exy(i)
         gstr(i,4)=gstr(i,4)+eyz(i)
         gstr(i,5)=gstr(i,5)+exz(i)
         gstr(i,6)=gstr(i,6)+kxx(i)
         gstr(i,7)=gstr(i,7)+kyy(i)
         gstr(i,8)=gstr(i,8)+kxy(i)
C-------- WXY         
         gstrw(i)=gstrw(i)+wxy(i)*dt1c(i)
        ENDDO
C-----   
        DO i=jft,jlt 
          f_def(i,6) = gstr(i,6)
          f_def(i,7) = gstr(i,7)
          f_def(i,8) = (gstr(i,8)+gstrw(i))*half
          f_def(i,5) = (gstr(i,8)-gstrw(i))*half
        ENDDO
       ELSEIF(ismstr /= 11)THEN
        DO i=jft,jlt
         gstr(i,1)=gstr(i,1)+exx(i)
         gstr(i,2)=gstr(i,2)+eyy(i)
         gstr(i,3)=gstr(i,3)+exy(i)
         gstr(i,4)=gstr(i,4)+eyz(i)
         gstr(i,5)=gstr(i,5)+exz(i)
         gstr(i,6)=gstr(i,6)+kxx(i)
         gstr(i,7)=gstr(i,7)+kyy(i)
         gstr(i,8)=gstr(i,8)+kxy(i)
        ENDDO
       ELSE
        DO i=jft,jlt
         fac1(i)  =one/area(i)
         ux13=ux1(i)-ux3(i)
         ux24=ux2(i)-ux4(i)
         uy13=uy1(i)-uy3(i)
         uy24=uy2(i)-uy4(i)
         exxt=(px1(i)*ux13+px2(i)*ux24)*fac1(i)
         eyyt=(py1(i)*uy13+py2(i)*uy24)*fac1(i)
         exyt=(py1(i)*ux13+py2(i)*ux24
     .        +px1(i)*uy13+px2(i)*uy24)*fac1(i)
         gstr(i,1)=exxt
         gstr(i,2)=eyyt
         gstr(i,3)=exyt
         gstr(i,4)=gstr(i,4)+eyz(i)
         gstr(i,5)=gstr(i,5)+exz(i)
         gstr(i,6)=gstr(i,6)+kxx(i)
         gstr(i,7)=gstr(i,7)+kyy(i)
         gstr(i,8)=gstr(i,8)+kxy(i)
        ENDDO
       ENDIF
      ENDIF
C
      RETURN