sigeps01gpinch.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sigeps01gpinch (jft, jlt, g, thk, off, gs, a1, a2, nu, thk0, nel, ssp, rho, depsxx, depsyy, depszz, depsxy, depsyz, depszx, depbxx, depbyy, depbxy, deppxz, deppyz, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, momoxx, momoyy, momoxy, momopxz, momopyz, signxx, signyy, signzz, signxy, signyz, signzx, momnxx, momnyy, momnxy, momnpxz, momnpyz)

Function/Subroutine Documentation

sigeps01gpinch()

subroutine sigeps01gpinch	(	integer	jft,
		integer	jlt,
			g,
			thk,
			off,
			gs,
			a1,
			a2,
			nu,
			thk0,
		integer	nel,
			ssp,
			rho,
			depsxx,
			depsyy,
			depszz,
			depsxy,
			depsyz,
			depszx,
			depbxx,
			depbyy,
			depbxy,
			deppxz,
			deppyz,
			sigoxx,
			sigoyy,
			sigozz,
			sigoxy,
			sigoyz,
			sigozx,
			momoxx,
			momoyy,
			momoxy,
			momopxz,
			momopyz,
			signxx,
			signyy,
			signzz,
			signxy,
			signyz,
			signzx,
			momnxx,
			momnyy,
			momnxy,
			momnpxz,
			momnpyz )

Definition at line 28 of file sigeps01gpinch.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   I N P U T   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT,NEL,ISMSTR
C     REAL
      my_real
     .  off(*),gs(*),thk(*)
      my_real
     .   nu(mvsiz),g(mvsiz),e(mvsiz),a1(mvsiz),a2(mvsiz),
     .   depsxx(nel),depsyy(nel),depszz(nel),
     .   depsxy(nel),depsyz(nel),depszx(nel),
     .   depbxx(nel),depbyy(nel),depbxy(nel),
     .   deppxz(nel),deppyz(nel),
     .   sigoxx(nel),sigoyy(nel),sigozz(nel),
     .   sigoxy(nel),sigoyz(nel),sigozx(nel),
     .   momoxx(nel),momoyy(nel),momoxy(nel),
     .   momopxz(nel),momopyz(nel)
C-----------------------------------------------
C   O U T P U T   A R G U M E N T S
C-----------------------------------------------
      my_real
     .    signxx(nel),signyy(nel),signzz(nel),
     .    signxy(nel),signyz(nel),signzx(nel),
     .    momnxx(nel),momnyy(nel),momnxy(nel),
     .    momnpxz(nel),momnpyz(nel),
     .    ssp(nel),rho(nel)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX,J
C     REAL
      my_real
     .   pa1(mvsiz),pa2(mvsiz),pa3(mvsiz),pa4(mvsiz),pa5(mvsiz),
     .   b1(mvsiz),b2(mvsiz),b3(mvsiz),b4(mvsiz),b5(mvsiz),
     .   thk0(mvsiz),thk08(mvsiz)
C-----------------------------------------------
C
      DO i=jft,jlt
C         [PM] reverse engineer E and ANU from A11 and A12
C         [PM] can't use A11/A12 since they correspond to shell (plane stress)
          e(i)   = a1(i)*(one-nu(i)**2)
C
C         PA1 = C1111, C2222, C3333
C         PA2 = C1122, C1133, C2233
C         PA3 = C1212, C1313, C2323
C         PA4 = C1111*, C2222* 
C         PA5 = C1122*
          pa1(i) = e(i)*(one-nu(i))/(one+nu(i))/(one-two*nu(i))          
          pa2(i) = e(i)*nu(i)/(one+nu(i))/(one-two*nu(i))         
          pa3(i) = e(i)/two/(one+nu(i))         
          pa4(i) = (one+nu(i))*(one-two*nu(i))/(one-nu(i)**2)/(one-nu(i))*pa1(i)         
          pa5(i) = (one+nu(i))*(one-two*nu(i))/(one-nu(i)**2)*pa2(i)
C
          thk08(i) = thk0(i)*one_over_12
          b1(i) = pa1(i)*thk08(i)
          b2(i) = pa2(i)*thk08(i)
          b3(i) = pa3(i)*thk08(i)
          b4(i) = pa4(i)*thk08(i)
          b5(i) = pa5(i)*thk08(i)
      ENDDO
C
      DO i=jft,jlt
C
        signxx(i)=sigoxx(i)+pa1(i)*depsxx(i)+pa2(i)*depsyy(i)+pa2(i)*depszz(i)
        signyy(i)=sigoyy(i)+pa2(i)*depsxx(i)+pa1(i)*depsyy(i)+pa2(i)*depszz(i)
        signzz(i)=sigozz(i)+pa2(i)*depsxx(i)+pa2(i)*depsyy(i)+pa1(i)*depszz(i)
C
        signxy(i)=sigoxy(i)+g(i)*depsxy(i)
        signyz(i)=sigoyz(i)+gs(i)*(depsyz(i)+zero*deppyz(i))
        signzx(i)=sigozx(i)+gs(i)*(depszx(i)+zero*deppxz(i))
C
        momnxx(i)=momoxx(i)+b4(i)*depbxx(i)+b5(i)*depbyy(i)
        momnyy(i)=momoyy(i)+b5(i)*depbxx(i)+b4(i)*depbyy(i)
        momnxy(i)=momoxy(i)+b3(i)*depbxy(i)
C
        momnpxz(i)=momopxz(i)+b3(i)*deppxz(i)
        momnpyz(i)=momopyz(i)+b3(i)*deppyz(i)
C
C       Sound speed
C       SSP(I) = SQRT(E(I)/RHO(I))                      ! historic
C       SSP(I) = SQRT(E(I)/RHO(I)/(ONE-ANU(I)*ANU(I)))   ! for true shells
        ssp(i) = sqrt(pa1(i)/rho(i))                    ! like solids
C
      ENDDO
      RETURN