srota6_s8s.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	srota6_s8s (kcvt, tens, gama, khbe, ityp, frame, iint, isorth)

Function/Subroutine Documentation

srota6_s8s()

subroutine srota6_s8s	(	integer	kcvt,
			tens,
			gama,
		integer	khbe,
		integer	ityp,
			frame,
		integer, intent(in)	iint,
		integer, intent(in)	isorth )

Definition at line 30 of file srota6_s8s.F.

      use element_mod , only : nixs
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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: IINT
      INTEGER, INTENT(IN) :: ISORTH
C     REAL
      my_real
     .   tens(6),gama(6),frame(3,3)
      INTEGER  KCVT, KHBE, ITYP
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C     REAL
      my_real
     .   l11,l12,l13,l22,l23,l33,
     .   r11,r12,r13,r21,r22,r23,r31,r32,r33,
     .   g11,g22,g33,g12,g21,g23,g32,g13,g31,
     .   t11,t22,t33,t12,t21,t23,t32,t13,t31,
     .   s11,s12,s21,s13,s31,s22,s23,s32,s33
C-----------------------------------------------
C---- should not rotate w/ global iso system
!       IF (KCVT==0.AND.ISORTH==0) RETURN
C----------------------------
      IF (khbe /= 17 .OR. iint /= 3) RETURN
         r11 = frame(1,1)
             r21 = frame(2,1)
             r31 = frame(3,1)
             r12 = frame(1,2)
             r22 = frame(2,2)
             r32 = frame(3,2)
             r13 = frame(1,3)
             r23 = frame(2,3)
             r33 = frame(3,3)
!       ENDIF
C-----------
C     ORTHOTROPIC ELEMENTS : Rotation from ORHTOROPIC FRAME TO LOCAL FRAME
C-----------
      IF (kcvt==2) THEN
        IF (isorth > 0) THEN
        g11=gama(1)
        g21=gama(2)
        g31=gama(3)
        g12=gama(4)
        g22=gama(5)
        g32=gama(6)
        g13=g21*g32-g31*g22
        g23=g31*g12-g11*g32
        g33=g11*g22-g21*g12
C       MATRICE DE PASSAGE GLOBAL -> ORTHOTROPE.
        t11=r11*g11+r12*g21+r13*g31
        t12=r11*g12+r12*g22+r13*g32
        t13=r11*g13+r12*g23+r13*g33
        t21=r21*g11+r22*g21+r23*g31
        t22=r21*g12+r22*g22+r23*g32
        t23=r21*g13+r22*g23+r23*g33
        t31=r31*g11+r32*g21+r33*g31
        t32=r31*g12+r32*g22+r33*g32
        t33=r31*g13+r32*g23+r33*g33
        r11=t11
        r12=t12
        r13=t13
        r21=t21
        r22=t22
        r23=t23
        r31=t31
        r32=t32
        r33=t33
       ENDIF
      ENDIF
C-----------
C      Rotation from LOCAL FRAME TO GLOBAL FRAME
C-----------
C     TENSOR ROTATION.
C        R passage Global -> Local.
C        R S Transpose(R)
C    1     SX
C    2     SY
C    3     SZ
C    4     SXY
C    5     SYZ
C    6     SXZ
C-----------
        l11    =tens(1)
        l22    =tens(2)
        l33    =tens(3)
        l12    =tens(4)
        l23    =tens(5)
        l13    =tens(6)
        s11    =l11*r11+l12*r12+l13*r13 
        s12    =l11*r21+l12*r22+l13*r23 
        s13    =l11*r31+l12*r32+l13*r33
        s21    =l12*r11+l22*r12+l23*r13
        s22    =l12*r21+l22*r22+l23*r23
        s23    =l12*r31+l22*r32+l23*r33 
        s31    =l13*r11+l23*r12+l33*r13
        s32    =l13*r21+l23*r22+l33*r23
        s33    =l13*r31+l23*r32+l33*r33
        tens(1)=r11*s11+r12*s21+r13*s31
        tens(2)=r21*s12+r22*s22+r23*s32
        tens(3)=r31*s13+r32*s23+r33*s33
        tens(4)=r11*s12+r12*s22+r13*s32
        tens(5)=r21*s13+r22*s23+r23*s33
        tens(6)=r11*s13+r12*s23+r13*s33
C-----------
      RETURN