c3stra3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr14_c.inc"
#include "param_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	c3stra3 (jft, jlt, pm, mat, area, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, geo, pid, nu, shf, gstr, ssp, rho, epsdot, nft, istrain, ismstr, ux1, ux2, ux3, uy1, uy2, uy3, px1, py1, py2, mtn, f_def, wxy, gstrw, nel)

Function/Subroutine Documentation

c3stra3()

subroutine c3stra3	(	integer	jft,
		integer	jlt,
			pm,
		integer, dimension(mvsiz)	mat,
			area,
			exx,
			eyy,
			exy,
			exz,
			eyz,
			kxx,
			kyy,
			kxy,
			geo,
		integer, dimension(mvsiz)	pid,
			nu,
			shf,
			gstr,
			ssp,
			rho,
			epsdot,
		integer	nft,
		integer	istrain,
		integer	ismstr,
			ux1,
			ux2,
			ux3,
			uy1,
			uy2,
			uy3,
			px1,
			py1,
			py2,
		integer	mtn,
			f_def,
			wxy,
			gstrw,
		integer	nel )

Definition at line 29 of file c3stra3.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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr14_c.inc"
#include      "param_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT,NFT,ISTRAIN,ISMSTR,NEL
      INTEGER MAT(MVSIZ),PID(MVSIZ),MTN
C     REAL
      my_real
     .   gstr(nel,8), pm(npropm,*),geo(npropg,*),
     .   area(mvsiz),shf(mvsiz),
     .   exx(mvsiz), eyy(mvsiz), exy(mvsiz), exz(mvsiz), eyz(mvsiz),
     .   kxx(mvsiz), kyy(mvsiz), kxy(mvsiz),nu(mvsiz),ssp(*),rho(*),
     .   epsdot(6,*),
     .   ux1(*),ux2(*),ux3(*),uy1(*),uy2(*),uy3(*),
     .   px1(*),py1(*),py2(*),wxy(*),f_def(mvsiz,*),gstrw(*)
      my_real
     .   ux12,ux13,ux23,uy13,uy23,exxt,eyyt,exyt
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ISH(MVSIZ),
     .        I, MX, J
C     REAL
      my_real
     .   fsh(mvsiz),
     .   thk2 ,fac1, fac2, fc1p
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C
      IF(iepsdot/=0)THEN
        DO i=jft,jlt
          fac1  =one/area(i)
          j = i + nft
          epsdot(1,j)=exx(i)*fac1
          epsdot(2,j)=eyy(i)*fac1
          epsdot(3,j)=exy(i)*fac1
          epsdot(4,j)=kxx(i)*fac1
          epsdot(5,j)=kyy(i)*fac1
          epsdot(6,j)=kxy(i)*fac1
        ENDDO
      ENDIF
C
      DO i=jft,jlt
        fac1  =dt1/area(i)
        exx(i)=exx(i)*fac1
        eyy(i)=eyy(i)*fac1
        exy(i)=exy(i)*fac1
        eyz(i)=eyz(i)*fac1
        exz(i)=exz(i)*fac1
        kxx(i)=kxx(i)*fac1
        kyy(i)=kyy(i)*fac1
        kxy(i)=kxy(i)*fac1
      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)*dt1/area(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  =one/area(i)
          ux12=ux1(i)-ux2(i)
          ux13=ux1(i)-ux3(i)
          ux23=ux2(i)-ux3(i)
          uy13=uy1(i)-uy3(i)
          uy23=uy2(i)-uy3(i)
          exxt= px1(i)*ux12*fac1
          eyyt=(py1(i)*uy13+py2(i)*uy23)*fac1
          exyt=(py1(i)*ux13+py2(i)*ux23)*fac1
          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