s16sigs3.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s16sigs3 (dxx, dyy, dzz, d4, d5, d6, rx, ry, rz, sx, sy, sz, tx, ty, tz, eint, vol0, sigs, mxt, pm, dt1, vol0dp, nel)

Function/Subroutine Documentation

s16sigs3()

subroutine s16sigs3	(		dxx,
			dyy,
			dzz,
			d4,
			d5,
			d6,
			rx,
			ry,
			rz,
			sx,
			sy,
			sz,
			tx,
			ty,
			tz,
			eint,
			vol0,
			sigs,
		integer, dimension(*)	mxt,
			pm,
			dt1,
		double precision, dimension(*)	vol0dp,
		integer, intent(in)	nel )

Definition at line 28 of file s16sigs3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER MXT(*) 
      my_real
     1    dxx(*) ,dyy(*) ,dzz(*) ,d4(*) ,d5(*) ,d6(*) ,
     2    rx(*) ,ry(*) ,rz(*) ,sx(*) ,sy(*) ,sz(*) ,
     3    tx(*) ,ty(*) ,tz(*) ,eint(*),vol0(*),sigs(*), 
     4    pm(npropm,*), dt1   
      DOUBLE PRECISION 
     .        VOL0DP(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX
      my_real
     .   aa, dxy, dyz, dzx, drr, dss, dtt, drs, dst, dtr,
     .   shf,a11,a12,a13,a21,a22,a23,a31,a32,a33,e,nu,d11
C-----------------------------------------------
      mx = mxt(1)
      e  = pm(20,mx)
      nu = pm(21,mx)
      DO i=1,nel
        shf = five_over_6
        aa  = em01 * sigs(i)* ((one +nu) * (one -nu-nu)) / (e * (one -nu))
 
        dxy = half * d4(i)
        dyz = half * d5(i)
        dzx = half * d6(i)
 
        a11 = dxx(i)*rx(i)+dxy*ry(i)+dzx*rz(i)
        a12 = dxx(i)*sx(i)+dxy*sy(i)+dzx*sz(i)
        a13 = dxx(i)*tx(i)+dxy*ty(i)+dzx*tz(i)
        a21 = dxy*rx(i)+dyy(i)*ry(i)+dyz*rz(i)
        a22 = dxy*sx(i)+dyy(i)*sy(i)+dyz*sz(i)
        a23 = dxy*tx(i)+dyy(i)*ty(i)+dyz*tz(i)
        a31 = dzx*rx(i)+dyz*ry(i)+dzz(i)*rz(i)
        a32 = dzx*sx(i)+dyz*sy(i)+dzz(i)*sz(i)
        a33 = dzx*tx(i)+dyz*ty(i)+dzz(i)*tz(i)
 
        drr = rx(i)*a11+ry(i)*a21+rz(i)*a31
        dss = sx(i)*a12+sy(i)*a22+sz(i)*a32
        dtt = tx(i)*a13+ty(i)*a23+tz(i)*a33
        dtr = tx(i)*a11+ty(i)*a21+tz(i)*a31
        drs = rx(i)*a12+ry(i)*a22+rz(i)*a32
        dst = sx(i)*a13+sy(i)*a23+sz(i)*a33
 
        vol0(i) = vol0(i) * (one + aa)
        vol0dp(i) = vol0dp(i) * (one + aa)
        eint(i) = eint(i) / (1. + aa)
        dss     = dss - aa/dt1
 
        drs     = drs * shf
        dst     = dst * shf
 
        a11 = drr*rx(i)+drs*sx(i)+dtr*tx(i)
        a12 = drr*ry(i)+drs*sy(i)+dtr*ty(i)
        a13 = drr*rz(i)+drs*sz(i)+dtr*tz(i)
        a21 = drs*rx(i)+dss*sx(i)+dst*tx(i)
        a22 = drs*ry(i)+dss*sy(i)+dst*ty(i)
        a23 = drs*rz(i)+dss*sz(i)+dst*tz(i)
        a31 = dtr*rx(i)+dst*sx(i)+dtt*tx(i)
        a32 = dtr*ry(i)+dst*sy(i)+dtt*ty(i)
        a33 = dtr*rz(i)+dst*sz(i)+dtt*tz(i)
 
        dxx(i)= rx(i)*a11+sx(i)*a21+tx(i)*a31
        dyy(i)= ry(i)*a12+sy(i)*a22+ty(i)*a32
        dzz(i)= rz(i)*a13+sz(i)*a23+tz(i)*a33
        dzx   = rz(i)*a11+sz(i)*a21+tz(i)*a31
        dxy   = rx(i)*a12+sx(i)*a22+tx(i)*a32
        dyz   = ry(i)*a13+sy(i)*a23+ty(i)*a33
 
        d4(i) = two * dxy 
        d5(i) = two * dyz 
        d6(i) = two * dzx
      ENDDO
C-----------
      RETURN