starter_2source_2output_2anim_2srota6_8F_source.html

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     ../starter/source/output/anim/srota6.F

!||--- called by ------------------------------------------------------

!||    tensors    ../starter/source/output/anim/tensors.F

!||--- calls      -----------------------------------------------------

!||    sortho31   ../starter/source/output/anim/sortho31.F

!||====================================================================


      SUBROUTINE srota6(X,IXS,KCVT,

     .                   TENS,GAMA)

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-----------------------------------------------

      my_real

     .   x(3,*),tens(6),gama(6)

      INTEGER IXS(NIXS), KCVT

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      my_real

     .   x1, x2, x3, x4,

     .   x5, x6, x7, x8,

     .   y1, y2, y3, y4,

     .   y5, y6, y7, y8,

     .   z1, z2, z3, z4,

     .   z5, z6, z7, z8,

     .   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

      INTEGER NC1, NC2, NC3, NC4,

     .        nc5, nc6, nc7, nc8

C-----------------------------------------------

        nc1=ixs(2)

        nc2=ixs(3)

        nc3=ixs(4)

        nc4=ixs(5)

        nc5=ixs(6)

        nc6=ixs(7)

        nc7=ixs(8)

        nc8=ixs(9)

C----------------------------

C     NODAL COORDINATES     |

C----------------------------

        x1=x(1,nc1)

        y1=x(2,nc1)

        z1=x(3,nc1)

        x2=x(1,nc2)

        y2=x(2,nc2)

        z2=x(3,nc2)

        x3=x(1,nc3)

        y3=x(2,nc3)

        z3=x(3,nc3)

        x4=x(1,nc4)

        y4=x(2,nc4)

        z4=x(3,nc4)

        x5=x(1,nc5)

        y5=x(2,nc5)

        z5=x(3,nc5)

        x6=x(1,nc6)

        y6=x(2,nc6)

        z6=x(3,nc6)

        x7=x(1,nc7)

        y7=x(2,nc7)

        z7=x(3,nc7)

        x8=x(1,nc8)

        y8=x(2,nc8)

        z8=x(3,nc8)

C-----------

C     REPERE CONVECTE.

C-----------

      CALL sortho31(

     .   x1, x2, x3, x4, x5, x6, x7, x8,

     .   y1, y2, y3, y4, y5, y6, y7, y8,

     .   z1, z2, z3, z4, z5, z6, z7, z8,

     .   r11, r12, r13, r21, r22, r23, r31, r32, r33)

C-----------

C     REPERE ORTHOTROPE.

C-----------

      IF (kcvt==2) 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

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

my_real
#define my_real
Definition cppsort.cpp:32

sortho31
subroutine sortho31(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition sortho31.F:34

srota6
subroutine srota6(x, ixs, kcvt, tens, gama)
Definition srota6.F:32