mrotens.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/.
!||====================================================================
!||    mrotens         ../engine/source/materials/mat_share/mrotens.F
!||--- called by ------------------------------------------------------
!||    m25law          ../engine/source/materials/mat/mat025/m25law.F
!||    mmain           ../engine/source/materials/mat_share/mmain.F90
!||    mulaw           ../engine/source/materials/mat_share/mulaw.f90
!||    mulaw8          ../engine/source/materials/mat_share/mulaw8.F90
!||    usermat_solid   ../engine/source/materials/mat_share/usermat_solid.F
!||====================================================================
      SUBROUTINE mrotens (LFT,LLT,E1,E2,E3,
     .                   E4,DYZ,E6,R11,R12,
     .                   R13,R21,R22,R23,R31,
     .                   R32,R33)
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LFT,LLT
      my_real E1(*),E2(*),E3(*),E4(*),DYZ(*),E6(*),
     .   R11(*), R12(*), R13(*), 
     .   r21(*), r22(*), r23(*),
     .   r31(*), r32(*), r33(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   A11(MVSIZ), A12(MVSIZ), A13(MVSIZ), 
     .   a21(mvsiz), a22(mvsiz), a23(mvsiz),
     .   a31(mvsiz), a32(mvsiz), a33(mvsiz)
C-----------------------------------------------
      DO 200 i=lft,llt
      a11(i)=e1(i)*r11(i)+e4(i)*r21(i)+e6(i)*r31(i)
      a12(i)=e1(i)*r12(i)+e4(i)*r22(i)+e6(i)*r32(i)
      a13(i)=e1(i)*r13(i)+e4(i)*r23(i)+e6(i)*r33(i)
      a21(i)=e4(i)*r11(i)+e2(i)*r21(i)+dyz(i)*r31(i)
      a22(i)=e4(i)*r12(i)+e2(i)*r22(i)+dyz(i)*r32(i)
      a23(i)=e4(i)*r13(i)+e2(i)*r23(i)+dyz(i)*r33(i)
      a31(i)=e6(i)*r11(i)+dyz(i)*r21(i)+e3(i)*r31(i)
      a32(i)=e6(i)*r12(i)+dyz(i)*r22(i)+e3(i)*r32(i)
      a33(i)=e6(i)*r13(i)+dyz(i)*r23(i)+e3(i)*r33(i)
 200  CONTINUE
      DO 201 i=lft,llt
      e1(i)=r11(i)*a11(i)+r21(i)*a21(i)+r31(i)*a31(i)
      e2(i)=r12(i)*a12(i)+r22(i)*a22(i)+r32(i)*a32(i)
      e3(i)=r13(i)*a13(i)+r23(i)*a23(i)+r33(i)*a33(i)
      e4(i)=r11(i)*a12(i)+r21(i)*a22(i)+r31(i)*a32(i)
      dyz(i)=r12(i)*a13(i)+r22(i)*a23(i)+r32(i)*a33(i)
      e6(i)=r11(i)*a13(i)+r21(i)*a23(i)+r31(i)*a33(i)
 201  CONTINUE
C
      RETURN
      END
!||====================================================================
!||    mrotensns       ../engine/source/materials/mat_share/mrotens.f
!||--- called by ------------------------------------------------------
!||    mulaw           ../engine/source/materials/mat_share/mulaw.F90
!||    usermat_solid   ../engine/source/materials/mat_share/usermat_solid.F
!||====================================================================
      SUBROUTINE mrotensns (LFT,LLT,
     .               DXX,DXY,DXZ,
     .               DYX,DYY,DYZ,
     .               DZX,DZY,DZZ,
     .               R11,R12,R13,
     .               R21,R22,R23,
     .               R31,R32,R33)
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 LFT,LLT
      my_real
     .   DXY(*) ,DYX(*) ,DYZ(*) ,DZY(*) ,DZX(*) ,DXZ(*)
      my_real DXX(*),DYY(*),DZZ(*),
     .   R11(*), R12(*), R13(*), 
     .   R21(*), R22(*), R23(*),
     .   R31(*), R32(*), R33(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   a11, a12, a13, 
     .   a21, a22, a23,
     .   a31, a32, a33
C-----------------------------------------------
      DO i=lft,llt
        a11 = dxx(i)*r11(i) + dxy(i)*r21(i) + dxz(i)*r31(i)
        a12 = dxx(i)*r12(i) + dxy(i)*r22(i) + dxz(i)*r32(i)
        a13 = dxx(i)*r13(i) + dxy(i)*r23(i) + dxz(i)*r33(i)
 
        a21 = dyx(i)*r11(i) + dyy(i)*r21(i) + dyz(i)*r31(i)
        a22 = dyx(i)*r12(i) + dyy(i)*r22(i) + dyz(i)*r32(i)
        a23 = dyx(i)*r13(i) + dyy(i)*r23(i) + dyz(i)*r33(i)
 
        a31 = dzx(i)*r11(i) + dzy(i)*r21(i) + dzz(i)*r31(i)
        a32 = dzx(i)*r12(i) + dzy(i)*r22(i) + dzz(i)*r32(i)
        a33 = dzx(i)*r13(i) + dzy(i)*r23(i) + dzz(i)*r33(i)
 
        dxx(i) = r11(i)*a11 + r21(i)*a21 + r31(i)*a31
        dxy(i) = r11(i)*a12 + r21(i)*a22 + r31(i)*a32
        dxz(i) = r11(i)*a13 + r21(i)*a23 + r31(i)*a33
 
        dyx(i) = r12(i)*a11 + r22(i)*a21 + r32(i)*a31
        dyy(i) = r12(i)*a12 + r22(i)*a22 + r32(i)*a32
        dyz(i) = r12(i)*a13 + r22(i)*a23 + r32(i)*a33
 
        dzx(i) = r13(i)*a11 + r23(i)*a21 + r33(i)*a31
        dzy(i) = r13(i)*a12 + r23(i)*a22 + r33(i)*a32
        dzz(i) = r13(i)*a13 + r23(i)*a23 + r33(i)*a33
      ENDDO
C
      RETURN
      END