animig3d.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/.
!||====================================================================
!||    animig3d       ../engine/source/output/anim/generate/animig3d.F
!||--- called by ------------------------------------------------------
!||    genani         ../engine/source/output/anim/generate/genani.F
!||--- calls      -----------------------------------------------------
!||    projecig3d     ../engine/source/elements/ige3d/projecig3d.f
!||--- uses       -----------------------------------------------------
!||    elbufdef_mod   ../common_source/modules/mat_elem/elbufdef_mod.F90
!||    message_mod    ../engine/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE animig3d(ELBUF_TAB,IPARG    ,X     ,D        , V       ,A     ,
     .                    WIGE     ,KXIG3D   ,IXIG3D,IG3DSOLID,NANIM3D_L,X_TEMP,
     .                    D_TEMP   ,V_TEMP   ,A_TEMP,TABSTRESL,IGEO     ,KNOT  ,
     .                    ITAB     ,IPARTIG3D,IPART , CONT, CONT_TEMP,
     .                    FINT, FINT_TEMP, FEXT, FEXT_TEMP, FREAC, FREAC_TEMP,KNOTLOCPC,KNOTLOCEL)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod  
      USE message_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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
#include      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER KXIG3D(NIXIG3D,*),IXIG3D(*),IGEO(NPROPGI,*),
     .        NANIM3D_L, IPARG(NPARG,*), ITAB(*)  ! DISCUTER AVEC SEB DE COMMENT RENVOYER NANIM3D_L
      INTEGER IG3DSOLID(8,27,*),IPARTIG3D(*),IPART(LIPART1,*)
      TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
      my_real
     .   x(3,*),d(3,*),v(3,*),a(3,*),wige(*),knot(*),
     .   x_temp(3,*),d_temp(3,*),v_temp(3,*),
     .   a_temp(3,*),tabstresl(6,*),cont(3,*),cont_temp(3,*),
     .   fint(3,*),fint_temp(3,*),fext(3,*),fext_temp(3,*),
     .   freac(6,*),freac_temp(3,*),knotlocpc(*),knotlocel(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NISOELCUT, NG, NBG, NCTRL, ITEST, NEL_P, I, PX,PY,PZ
C----------------------------------------------------------
      NISOELCUT = 0
      ng=1
      DO WHILE (ng<=ngroup)
       nel_p = 0
       itest = 0
       IF (iparg(5,ng)==101) THEN 
        nft=iparg(3,ng)
        IF (ng==ngroup) THEN
         nel_p = iparg(2,ng)
        ELSE
         i=nft+1
         nel_p=nel_p+1
         DO WHILE (i<numelig3d)
          IF (ipartig3d(i+1)==ipartig3d(i)) THEN
            nel_p=nel_p+1
            i=i+1
          ELSE
            EXIT
          END IF
         END DO
        ENDIF
        nbg  = (nel_p-1)/nvsiz + 1
        nctrl = iparg(75,ng)
        px = igeo(41,iparg(62,ng))
        py = igeo(42,iparg(62,ng))
        pz = igeo(43,iparg(62,ng))
        CALL projecig3d(elbuf_tab, iparg    , x     , d, v,a,
     .                  wige     , kxig3d   , ixig3d,
     .                  ig3dsolid, nanim3d_l, x_temp, d_temp, v_temp,a_temp,
     .                  tabstresl, igeo     , knot  ,
     .                  ng       , nbg      , nisoelcut,
     .                  nctrl    , nel_p    , itab  ,cont , cont_temp,
     .                  fint, fint_temp, fext, fext_temp, freac, freac_temp,
     .                  px, py, pz,knotlocpc,knotlocel)
        ng=ng+nbg
       ELSE
        ng=ng+1
       ENDIF
      ENDDO
C-----------
      RETURN
      END