srota6_s8s.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    srota6_s8s    ../engine/source/output/anim/generate/srota6_s8s.F
!||--- called by ------------------------------------------------------
!||    tensors       ../engine/source/output/anim/generate/tensor6.F
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||====================================================================
      SUBROUTINE srota6_s8s(
     1   KCVT,    TENS,    GAMA,    KHBE,
     2   ITYP,    FRAME,   IINT,    ISORTH)
      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
      END