ranim33.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/.
!||====================================================================
!||    ranim33   ../engine/source/elements/joint/ranim33.F
!||--- called by ------------------------------------------------------
!||    rgjoint   ../engine/source/elements/joint/rgjoint.F
!||--- uses       -----------------------------------------------------
!||    h3d_mod   ../engine/share/modules/h3d_mod.F
!||====================================================================
      SUBROUTINE ranim33(JFT,JLT,FORX,FORY,FORZ,XMOM,YMOM,ZMOM,XL,TANI,H3D_DATA)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE h3d_mod
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 JFT,JLT
      my_real
     .   forx(*), fory(*), forz(*), xmom(*),
     .   ymom(*), zmom(*),tani(15,*)
      DOUBLE PRECISION XL(MVSIZ,3)
      TYPE(h3d_database) :: H3D_DATA
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr14_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real m1,m2,m3
C=======================================================================
C
      IF(anim_ft(1)/=0 .OR. h3d_data%UND_FORC /= 0)THEN
       DO i=jft,jlt
        tani(1,i)= forx(i)
        tani(2,i)= fory(i)
        tani(3,i)= forz(i)
        m1 = half*(xl(i,2)*forz(i)-xl(i,3)*fory(i))
        m2 = half*(xl(i,3)*forx(i)-xl(i,1)*forz(i))
        m3 = half*(xl(i,1)*fory(i)-xl(i,2)*forx(i))
        tani(4,i)= xmom(i)  + m1
        tani(5,i)= ymom(i)  + m2
        tani(6,i)= zmom(i)  + m3
        tani(7,i)= -xmom(i) + m1
        tani(8,i)= -ymom(i) + m2
        tani(9,i)= -zmom(i) + m3
       ENDDO
      ENDIF
C
      RETURN
      END