gloa24.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	gloa24 (nel, strain, dt1, arm1, arm2, arm3, d1, d2, d3, d4, d5, d6, r11, r12, r13, r21, r22, r23, r31, r32, r33, deps1, deps2, deps3, deps4, deps5, deps6)

Function/Subroutine Documentation

gloa24()

subroutine gloa24	(	integer	nel,
		intent(in)	strain,
			dt1,
		intent(inout)	arm1,
		intent(inout)	arm2,
		intent(inout)	arm3,
		intent(in)	d1,
		intent(in)	d2,
		intent(in)	d3,
		intent(in)	d4,
		intent(in)	d5,
		intent(in)	d6,
		intent(in)	r11,
		intent(in)	r12,
		intent(in)	r13,
		intent(in)	r21,
		intent(in)	r22,
		intent(in)	r23,
		intent(in)	r31,
		intent(in)	r32,
		intent(in)	r33,
		intent(out)	deps1,
		intent(out)	deps2,
		intent(out)	deps3,
		intent(out)	deps4,
		intent(out)	deps5,
		intent(out)	deps6 )

Definition at line 28 of file gloa24.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real dt1
      my_real, DIMENSION(NEL,6), INTENT(IN)  :: strain
      my_real, DIMENSION(NEL)  , INTENT(IN)  :: d1,d2,d3,d4,d5,d6,
     .                                          r11,r12,r13,r21,r22,r23,r31,r32,r33
      my_real, DIMENSION(NEL)  , INTENT(INOUT) :: arm1,arm2,arm3
      my_real, DIMENSION(NEL)  , INTENT(OUT) :: deps1,deps2,deps3,deps4,deps5,deps6
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real, DIMENSION(NEL) :: h1, h2, h3,h4, h5, h6,
     .                           a11, a12, a13, a21, a22, a23, a31,a32, a33
C=======================================================================
      DO i=1,nel
        deps1(i)=d1(i)*dt1
        deps2(i)=d2(i)*dt1
        deps3(i)=d3(i)*dt1
        deps4(i)=half*d4(i)*dt1
        deps5(i)=half*d5(i)*dt1
        deps6(i)=half*d6(i)*dt1
        a11(i)=deps1(i)*r11(i)+deps4(i)*r21(i)+deps6(i)*r31(i)
        a12(i)=deps1(i)*r12(i)+deps4(i)*r22(i)+deps6(i)*r32(i)
        a13(i)=deps1(i)*r13(i)+deps4(i)*r23(i)+deps6(i)*r33(i)
        a21(i)=deps4(i)*r11(i)+deps2(i)*r21(i)+deps5(i)*r31(i)
        a22(i)=deps4(i)*r12(i)+deps2(i)*r22(i)+deps5(i)*r32(i)
        a23(i)=deps4(i)*r13(i)+deps2(i)*r23(i)+deps5(i)*r33(i)
        a31(i)=deps6(i)*r11(i)+deps5(i)*r21(i)+deps3(i)*r31(i)
        a32(i)=deps6(i)*r12(i)+deps5(i)*r22(i)+deps3(i)*r32(i)
        a33(i)=deps6(i)*r13(i)+deps5(i)*r23(i)+deps3(i)*r33(i)
      ENDDO
      DO i=1,nel
        deps1(i)=r11(i)*a11(i)+r21(i)*a21(i)+r31(i)*a31(i)
        deps2(i)=r12(i)*a12(i)+r22(i)*a22(i)+r32(i)*a32(i)
        deps3(i)=r13(i)*a13(i)+r23(i)*a23(i)+r33(i)*a33(i)
        deps4(i)=r11(i)*a12(i)+r21(i)*a22(i)+r31(i)*a32(i)
        deps5(i)=r12(i)*a13(i)+r22(i)*a23(i)+r32(i)*a33(i)
        deps6(i)=r11(i)*a13(i)+r21(i)*a23(i)+r31(i)*a33(i)
        deps4(i)=two*deps4(i)
        deps5(i)=two*deps5(i)
        deps6(i)=two*deps6(i)
      ENDDO
C . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
      DO i=1,nel
        h1(i)=strain(i,1)
        h2(i)=strain(i,2)
        h3(i)=strain(i,3)
        h4(i)=half*strain(i,4)
        h5(i)=half*strain(i,5)
        h6(i)=half*strain(i,6)
        a11(i)=h1(i)*r11(i)+h4(i)*r21(i)+h6(i)*r31(i)
        a12(i)=h1(i)*r12(i)+h4(i)*r22(i)+h6(i)*r32(i)
        a13(i)=h1(i)*r13(i)+h4(i)*r23(i)+h6(i)*r33(i)
        a21(i)=h4(i)*r11(i)+h2(i)*r21(i)+h5(i)*r31(i)
        a22(i)=h4(i)*r12(i)+h2(i)*r22(i)+h5(i)*r32(i)
        a23(i)=h4(i)*r13(i)+h2(i)*r23(i)+h5(i)*r33(i)
        a31(i)=h6(i)*r11(i)+h5(i)*r21(i)+h3(i)*r31(i)
        a32(i)=h6(i)*r12(i)+h5(i)*r22(i)+h3(i)*r32(i)
        a33(i)=h6(i)*r13(i)+h5(i)*r23(i)+h3(i)*r33(i)
      ENDDO
      DO i=1,nel
        h1(i)=r11(i)*a11(i)+r21(i)*a21(i)+r31(i)*a31(i)
        h2(i)=r12(i)*a12(i)+r22(i)*a22(i)+r32(i)*a32(i)
        h3(i)=r13(i)*a13(i)+r23(i)*a23(i)+r33(i)*a33(i)
        arm1(i)=arm1(i) * max(zero,1.-h2(i)-h3(i))
        arm2(i)=arm2(i) * max(zero,1.-h1(i)-h3(i))
        arm3(i)=arm3(i) * max(zero,1.-h1(i)-h2(i))
      ENDDO
C
      RETURN