q4defo2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4defo2 (py1, py2, py3, py4, pz1, pz2, pz3, pz4, byz1, byz2, byz3, byz4, bzy1, bzy2, bzy3, bzy4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, eyz, eyy, ezz, exx, wxx, r22, r23, ay, off, offg, volo, eint, dsv, icp, fac, nel, jcvt)

Function/Subroutine Documentation

q4defo2()

subroutine q4defo2	(		py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			byz1,
			byz2,
			byz3,
			byz4,
			bzy1,
			bzy2,
			bzy3,
			bzy4,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			eyz,
			eyy,
			ezz,
			exx,
			wxx,
			r22,
			r23,
			ay,
			off,
			offg,
			volo,
			eint,
			dsv,
		integer	icp,
			fac,
		integer, intent(in)	nel,
		integer, intent(in)	jcvt )

Definition at line 28 of file q4defo2.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      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JCVT
      INTEGER ICP
C     REAL
      my_real
     .    py1(*),py2(*),py3(*),py4(*),
     .    pz1(*),pz2(*),pz3(*),pz4(*),
     .    byz1(*),byz2(*),byz3(*),byz4(*),
     .    bzy1(*),bzy2(*),bzy3(*),bzy4(*),
     .    vy1(*),vy2(*),vy3(*),vy4(*),
     .    vz1(*),vz2(*),vz3(*),vz4(*),
     .    eyz(*),eyy(*),ezz(*),exx(*),wxx(*),
     .    r22(*),r23(*),
     .    ay(*),off(*),offg(*),volo(*),
     .    eint(*),dsv(*),fac(*)
C-----------------------------------------------
c FUNCTION: 
c ARGUMENTS:  (I: input, O: output, IO: input & output, W: workspace)
c TYPE NAME                FUNCTION
c  I   PY1(*)~PZ4(*)      - SHAPE DERIVATIVES
c  I   VY1(*)~VZ4(*)      - NODAL VELOCITIES
c  O   EYY(*),EZZ(*),EXX(*) - RATE OF NORMAL STRAIN
c  O   WXX(*)             - COMPONENT OF ROTATION TENSOR
c  I   AY(*)              - Ni/r AT CENTER FOR AXISYMMETRIC CASE
c  IO  OFF(*)             - ELEMENT D/A FLAG, OF POINT
c  I   OFFG(*)            - ELEMENT D/A FLAG, OVERALL
c  IO  VOLO(*)            - EV(NB6)
c  IO  EINT(*)            - EV(NB3)
c  I   DSV(*)             - VOLUME STRAIN RATE AT CENTER
c  I   ICP                - FLAG FOR CONSTANT PRESURE
c  I   FAC(*)             - FRACTION ABOUT PLASTICITY STATE
c  I   R22(*)~R33(*)      - TRANSFORMATION MATRIX
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .    dyz(mvsiz),dzy(mvsiz),
     .    dvc(mvsiz),dv1,tol,
     .    dt1d2,vyg
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C     STRAIN RATE BEFORE CORRECTION
      DO i=1,nel
        eyy(i)=py1(i)*vy1(i)+py2(i)*vy2(i)+py3(i)*vy3(i)+py4(i)*vy4(i)
        ezz(i)=pz1(i)*vz1(i)+pz2(i)*vz2(i)+pz3(i)*vz3(i)+pz4(i)*vz4(i)
        dyz(i)=pz1(i)*vy1(i)+pz2(i)*vy2(i)+pz3(i)*vy3(i)+pz4(i)*vy4(i)
        dzy(i)=py1(i)*vz1(i)+py2(i)*vz2(i)+py3(i)*vz3(i)+py4(i)*vz4(i)
      ENDDO
C
Ccw>>>
C
C     STRAIN RATE AFTER CORRECTION AND INTRINSIC ROTATION
      dt1d2=half*dt1
      IF(jcvt==0) THEN
        DO i=1,nel
          wxx(i)=dt1d2*(dzy(i)-dyz(i))
        ENDDO
      ELSE
        DO i=1,nel
          wxx(i)=zero
          eyy(i)=eyy(i)-dt1d2*(eyy(i)*eyy(i)+dzy(i)*dzy(i))
          ezz(i)=ezz(i)-dt1d2*(ezz(i)*ezz(i)+dyz(i)*dyz(i))
        ENDDO
      ENDIF
C
Ccw<<<
      DO i=1,nel
        dvc(i)=eyy(i)+ezz(i)+exx(i)
      ENDDO
Ccw>>> ----- only assumed strain
        DO i=1,nel
         eyy(i)=eyy(i)+
     .     byz1(i)*vz1(i)+byz2(i)*vz2(i)+byz3(i)*vz3(i)+byz4(i)*vz4(i)
         ezz(i)=ezz(i)+
     .     bzy1(i)*vy1(i)+bzy2(i)*vy2(i)+bzy3(i)*vy3(i)+bzy4(i)*vy4(i)
        ENDDO
C
C
Ccw<<<
C     MODIFY NORMAL STRAIN RATE FOR CONSTANT PRESSURE CASE
c      IF(ICP>0) THEN
c       IF(N2D==1) THEN
c        IF(ICP==2) THEN
c          DO I=1,NEL
c            DVC(I)=FAC(I)*THIRD*(DSV(I)-DVC(I))
c            EYY(I)=EYY(I)+DVC(I)
c            EZZ(I)=EZZ(I)+DVC(I)
c            EXX(I)=EXX(I)+DVC(I)
c          ENDDO
c        ELSEIF(ICP==1) THEN
c          DO I=1,NEL
c            DVC(I)=THIRD*(DSV(I)-DVC(I))
c            EYY(I)=EYY(I)+DVC(I)
c            EZZ(I)=EZZ(I)+DVC(I)
c            EXX(I)=EXXC(I)+DVC(I)
c          ENDDO
c        ENDIF
c       ELSE
c        IF(ICP==2) THEN
c          DO I=1,NEL
c            DVC(I)=FAC(I)*HALF*(DSV(I)-DVC(I))
c            EYY(I)=EYY(I)+DVC(I)
c            EZZ(I)=EZZ(I)+DVC(I)
c          ENDDO
c        ELSEIF(ICP==1) THEN
c          DO I=1,NEL
c            DVC(I)=HALF*(DSV(I)-DVC(I))
c            EYY(I)=EYY(I)+DVC(I)
c            EZZ(I)=EZZ(I)+DVC(I)
c          ENDDO
c        ENDIF
c       END IF!(N2D==1) THEN
C----For hypoelastic laws-----       
c       DO I=1,NEL
c        IF(OFFG(I)==TWO) CYCLE
c   DV1 = ONE- DVC(I)*DT1
c         VOLO(I) = VOLO(I)*DV1
c         EINT(I) = EINT(I)/DV1
c       ENDDO
c      END IF!(ICP>0) THEN
Ccw>>>
C
      RETURN