udam24.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	udam24 (sig, ang, nel)
subroutine	udam24n (sig, ang, nbdama, damai, nel)

Function/Subroutine Documentation

udam24()

subroutine udam24	(		sig,
			ang,
		integer	nel )

Definition at line 28 of file udam24.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 sig(nel,6), ang(nel,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real anga
      my_real s11(nel), s12(nel), s13(nel), s21(nel), s22(nel), s23(nel), s31(nel),
     .        s32(nel), s33(nel), a11(nel), a12(nel), a13(nel), a21(nel), a22(nel),
     .        a23(nel), a31(nel), a32(nel), a33(nel)
C=======================================================================
      anga = zero
      DO i = 1,nel
        anga = anga + abs(ang(i,1)-one) + abs(ang(i,5)-one)
      ENDDO
      IF (anga <= em3) RETURN
C
      DO i = 1,nel
        s11(i) = ang(i,1)
        s21(i) = ang(i,2)
        s31(i) = ang(i,3)
        s12(i) = ang(i,4)
        s22(i) = ang(i,5)
        s32(i) = ang(i,6)
        s13(i) = s21(i)*s32(i)-s31(i)*s22(i)
        s23(i) = s31(i)*s12(i)-s11(i)*s32(i)
        s33(i) = s11(i)*s22(i)-s21(i)*s12(i)
      ENDDO
C . . . . . . . . . . . . . . . . . .
C     ROTATION INVERSE DES CONTRAINTES
C . . . . . . . . . . . . . . . . . .
      DO i = 1,nel
        a11(i)=sig(i,1)*s11(i)+sig(i,4)*s12(i)+sig(i,6)*s13(i)
        a12(i)=sig(i,1)*s21(i)+sig(i,4)*s22(i)+sig(i,6)*s23(i)
        a13(i)=sig(i,1)*s31(i)+sig(i,4)*s32(i)+sig(i,6)*s33(i)
        a21(i)=sig(i,4)*s11(i)+sig(i,2)*s12(i)+sig(i,5)*s13(i)
        a22(i)=sig(i,4)*s21(i)+sig(i,2)*s22(i)+sig(i,5)*s23(i)
        a23(i)=sig(i,4)*s31(i)+sig(i,2)*s32(i)+sig(i,5)*s33(i)
        a31(i)=sig(i,6)*s11(i)+sig(i,5)*s12(i)+sig(i,3)*s13(i)
        a32(i)=sig(i,6)*s21(i)+sig(i,5)*s22(i)+sig(i,3)*s23(i)
        a33(i)=sig(i,6)*s31(i)+sig(i,5)*s32(i)+sig(i,3)*s33(i)
      ENDDO
c
      DO i = 1,nel
        sig(i,1)=s11(i)*a11(i)+s12(i)*a21(i)+s13(i)*a31(i)
        sig(i,2)=s21(i)*a12(i)+s22(i)*a22(i)+s23(i)*a32(i)
        sig(i,3)=s31(i)*a13(i)+s32(i)*a23(i)+s33(i)*a33(i)
        sig(i,4)=s11(i)*a12(i)+s12(i)*a22(i)+s13(i)*a32(i)
        sig(i,5)=s21(i)*a13(i)+s22(i)*a23(i)+s23(i)*a33(i)
        sig(i,6)=s11(i)*a13(i)+s12(i)*a23(i)+s13(i)*a33(i)
      ENDDO
c-----------
      RETURN

udam24n()

subroutine udam24n	(		sig,
			ang,
		integer	nbdama,
		integer, dimension(nel)	damai,
		integer	nel )

Definition at line 97 of file udam24.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 NBDAMA,NEL
      INTEGER DAMAI(NEL)
      my_real sig(nel,6), ang(nel,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,K
      my_real anga,
     .   s11(nel), s12(nel), s13(nel), s21(nel), s22(nel), s23(nel), s31(nel),
     .   s32(nel), s33(nel), a11(nel), a12(nel), a13(nel), a21(nel), a22(nel), 
     .   a23(nel), a31(nel), a32(nel), a33(nel)
C=======================================================================
      DO k=1,nbdama
        i = damai(k)
        s11(i)=ang(i,1)
        s21(i)=ang(i,2)
        s31(i)=ang(i,3)
        s12(i)=ang(i,4)
        s22(i)=ang(i,5)
        s32(i)=ang(i,6)
        s13(i)=s21(i)*s32(i)-s31(i)*s22(i)
        s23(i)=s31(i)*s12(i)-s11(i)*s32(i)
        s33(i)=s11(i)*s22(i)-s21(i)*s12(i)
      ENDDO
C . . . . . . . . . . . . . . . . . .
C     ROTATION INVERSE DES CONTRAINTES
C . . . . . . . . . . . . . . . . . .
      DO k=1,nbdama
       i = damai(k)
       a11(i)=sig(i,1)*s11(i)+sig(i,4)*s12(i)+sig(i,6)*s13(i)
       a12(i)=sig(i,1)*s21(i)+sig(i,4)*s22(i)+sig(i,6)*s23(i)
       a13(i)=sig(i,1)*s31(i)+sig(i,4)*s32(i)+sig(i,6)*s33(i)
       a21(i)=sig(i,4)*s11(i)+sig(i,2)*s12(i)+sig(i,5)*s13(i)
       a22(i)=sig(i,4)*s21(i)+sig(i,2)*s22(i)+sig(i,5)*s23(i)
       a23(i)=sig(i,4)*s31(i)+sig(i,2)*s32(i)+sig(i,5)*s33(i)
       a31(i)=sig(i,6)*s11(i)+sig(i,5)*s12(i)+sig(i,3)*s13(i)
       a32(i)=sig(i,6)*s21(i)+sig(i,5)*s22(i)+sig(i,3)*s23(i)
       a33(i)=sig(i,6)*s31(i)+sig(i,5)*s32(i)+sig(i,3)*s33(i)
      ENDDO
c
      DO k=1,nbdama
       i = damai(k)
       sig(i,1)=s11(i)*a11(i)+s12(i)*a21(i)+s13(i)*a31(i)
       sig(i,2)=s21(i)*a12(i)+s22(i)*a22(i)+s23(i)*a32(i)
       sig(i,3)=s31(i)*a13(i)+s32(i)*a23(i)+s33(i)*a33(i)
       sig(i,4)=s11(i)*a12(i)+s12(i)*a22(i)+s13(i)*a32(i)
       sig(i,5)=s21(i)*a13(i)+s22(i)*a23(i)+s23(i)*a33(i)
       sig(i,6)=s11(i)*a13(i)+s12(i)*a23(i)+s13(i)*a33(i)
      ENDDO
c-----------
      RETURN