s8e_sig_a.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

s8e_sig_a (c1, g, off, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, pij, sig_a, eint, vol, d1, d2, d3, sig, s1, s2, s3, nel, offg)

Function/Subroutine Documentation

s8e_sig_a()

subroutine s8e_sig_a	(		c1,
			g,
			off,
			bxy1,
			bxy2,
			bxy3,
			bxy4,
			bxy5,
			bxy6,
			bxy7,
			bxy8,
			byx1,
			byx2,
			byx3,
			byx4,
			byx5,
			byx6,
			byx7,
			byx8,
			bxz1,
			bxz2,
			bxz3,
			bxz4,
			bxz5,
			bxz6,
			bxz7,
			bxz8,
			bzx1,
			bzx2,
			bzx3,
			bzx4,
			bzx5,
			bzx6,
			bzx7,
			bzx8,
			byz1,
			byz2,
			byz3,
			byz4,
			byz5,
			byz6,
			byz7,
			byz8,
			bzy1,
			bzy2,
			bzy3,
			bzy4,
			bzy5,
			bzy6,
			bzy7,
			bzy8,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vx7,
			vx8,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vy7,
			vy8,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			vz7,
			vz8,
			pij,
			sig_a,
			eint,
			vol,
			d1,
			d2,
			d3,
			sig,
			s1,
			s2,
			s3,
		integer	nel,
			offg )

Definition at line 28 of file s8e_sig_a.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"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), vx7(*), vx8(*),
     .   vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), vy7(*), vy8(*),
     .   vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), vz7(*), vz8(*),
     .   bxy1(*),bxy2(*),bxy3(*),bxy4(*),
     .   bxy5(*),bxy6(*),bxy7(*),bxy8(*),
     .   byx1(*),byx2(*),byx3(*),byx4(*),
     .   byx5(*),byx6(*),byx7(*),byx8(*),
     .   bxz1(*),bxz2(*),bxz3(*),bxz4(*),
     .   bxz5(*),bxz6(*),bxz7(*),bxz8(*),
     .   bzx1(*),bzx2(*),bzx3(*),bzx4(*),
     .   bzx5(*),bzx6(*),bzx7(*),bzx8(*),
     .   byz1(*),byz2(*),byz3(*),byz4(*),
     .   byz5(*),byz6(*),byz7(*),byz8(*),
     .   bzy1(*),bzy2(*),bzy3(*),bzy4(*),
     .   bzy5(*),bzy6(*),bzy7(*),bzy8(*),
     .   pij(nel,*),c1,g,sig_a(mvsiz,3),off(*),eint(*),
     .   vol(*),d1(*),d2(*),d3(*),sig(nel,6),s1(*),s2(*),s3(*),offg(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,K
      my_real
     .   dc1,dg2,dxx(mvsiz),dyy(mvsiz),dzz(mvsiz),ekk,dav,
     .   dt05,e1,e2,e3,de,p2
C PIJ(NEL,25-27) are used call stress with v-terms
C-----------------------------------------------
       IF (tt==zero) pij(1:nel,25:27)=zero
C----1,compute Dii with v-terms, compute incremental stress
       DO i=1,nel
               IF (offg(i) >one) cycle
           dyy(i) =bxy1(i)*vx1(i)+bxy2(i)*vx2(i)
     +              +bxy3(i)*vx3(i)+bxy4(i)*vx4(i)
     +              +bxy5(i)*vx5(i)+bxy6(i)*vx6(i)
     +              +bxy7(i)*vx7(i)+bxy8(i)*vx8(i)
     +              +bzy1(i)*vz1(i)+bzy2(i)*vz2(i)
     +              +bzy3(i)*vz3(i)+bzy4(i)*vz4(i)
     +              +bzy5(i)*vz5(i)+bzy6(i)*vz6(i)
     +              +bzy7(i)*vz7(i)+bzy8(i)*vz8(i)
           dzz(i) =bxz1(i)*vx1(i)+bxz2(i)*vx2(i)
     +              +bxz3(i)*vx3(i)+bxz4(i)*vx4(i)
     +              +bxz5(i)*vx5(i)+bxz6(i)*vx6(i)
     +              +bxz7(i)*vx7(i)+bxz8(i)*vx8(i)
     +              +byz1(i)*vy1(i)+byz2(i)*vy2(i)
     +              +byz3(i)*vy3(i)+byz4(i)*vy4(i)
     +              +byz5(i)*vy5(i)+byz6(i)*vy6(i)
     +              +byz7(i)*vy7(i)+byz8(i)*vy8(i)
           dxx(i) =byx1(i)*vy1(i)+byx2(i)*vy2(i)
     +              +byx3(i)*vy3(i)+byx4(i)*vy4(i)
     +              +byx5(i)*vy5(i)+byx6(i)*vy6(i)
     +              +byx7(i)*vy7(i)+byx8(i)*vy8(i)
     +              +bzx1(i)*vz1(i)+bzx2(i)*vz2(i)
     +              +bzx3(i)*vz3(i)+bzx4(i)*vz4(i)
     +              +bzx5(i)*vz5(i)+bzx6(i)*vz6(i)
     +              +bzx7(i)*vz7(i)+bzx8(i)*vz8(i)
       ENDDO
       dc1=c1*dt1
       dg2=two*g*dt1
       sig_a(1:nel,1:3)=pij(1:nel,25:27)
       DO i=1,nel
               IF (offg(i) >one) cycle
         ekk=(dxx(i)+dyy(i)+dzz(i))*dc1
c       IF (OFFG(I) >ONE) EKK=ZERO
         dav=-third*(dxx(i)+dyy(i)+dzz(i))
         pij(i,25)=pij(i,25)*off(i)+ekk+dg2*(dxx(i)+dav)
         pij(i,26)=pij(i,26)*off(i)+ekk+dg2*(dyy(i)+dav)
         pij(i,27)=pij(i,27)*off(i)+ekk+dg2*(dzz(i)+dav)
       END DO
       dt05 =half*dt1
      DO i=1,nel
              IF (offg(i) >one) cycle
       e1=(dxx(i)+d1(i))*(sig_a(i,1)+pij(i,25))+dxx(i)*(s1(i)+sig(i,1))
       e2=(dyy(i)+d2(i))*(sig_a(i,2)+pij(i,26))+dyy(i)*(s2(i)+sig(i,2))
       e3=(dzz(i)+d3(i))*(sig_a(i,3)+pij(i,27))+dzz(i)*(s3(i)+sig(i,3))
       de= (e1+e2+e3)*dt05*off(i) 
       eint(i)=eint(i)+de 
       sig_a(i,1:3)=pij(i,25:27)
      ENDDO
C-----------
      RETURN