grelem_sav.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/.
!||====================================================================
!||    grelem_sav   ../engine/source/output/th/grelem_sav.F
!||--- called by ------------------------------------------------------
!||    c3bilan      ../engine/source/elements/sh3n/coque3n/c3bilan.F
!||    cbilan       ../engine/source/elements/shell/coque/cbilan.F
!||    ige3dbilan   ../engine/source/elements/ige3d/ige3dbilan.F
!||    pbilan       ../engine/source/elements/beam/pbilan.F
!||    qbilan       ../engine/source/elements/solid_2d/quad/qbilan.F
!||    r23bilan     ../engine/source/elements/spring/r23bilan.F
!||    r3bilan      ../engine/source/elements/spring/r3bilan.F
!||    r5bilan      ../engine/source/elements/spring/r5bilan.F
!||    rbilan       ../engine/source/elements/spring/rbilan.F
!||    rbilan33     ../engine/source/elements/joint/rbilan33.F
!||    s10bilan     ../engine/source/elements/solid/solide10/s10bilan.F
!||    s16bilan     ../engine/source/elements/thickshell/solide16/s16bilan.F
!||    s20bilan     ../engine/source/elements/solid/solide20/s20bilan.F
!||    s4bilan      ../engine/source/elements/solid/solide4/s4bilan.F
!||    s6cbilan     ../engine/source/elements/thickshell/solide6c/s6cbilan.F
!||    s8bilan      ../engine/source/elements/solid/solide8/s8bilan.F
!||    sbilan       ../engine/source/elements/solid/solide/sbilan.f
!||    spbilan      ../engine/source/elements/sph/spbilan.F
!||    sr8bilan     ../engine/source/elements/solid/solide8/sr8bilan.F
!||    srbilan      ../engine/source/elements/solid/solide/srbilan.F
!||    tbilan       ../engine/source/elements/truss/tbilan.F
!||====================================================================
      SUBROUTINE grelem_sav(JFT   ,JLT   ,GRESAV,IGRTH ,GRTH  ,
     2                      OFF   ,EI    ,EK    ,XM    ,YM    ,
     3                      ZM    ,XMAS  ,XCG   ,YCG   ,ZCG   ,
     4                      XXM   ,YYM   ,ZZM   ,IXX   ,IYY   ,
     5                      IZZ   ,IXY   ,IYZ   ,IZX   ,REI   ,
     6                      REK   ,FLAG)
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-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT,JLT,FLAG
      INTEGER IGRTH(*),GRTH(*)
      my_real
     .   GRESAV(NPSAV,*),OFF(*),
     .   XMAS(*), EI(*),REI(*),REK(*),
     .   ek(*), xm(*), ym(*), zm(*),
     .   xxm(*), yym(*), zzm(*), 
     .   xcg(*), ycg(*), zcg(*),
     .   ixx(*), iyy(*), izz(*),
     .   ixy(*), iyz(*), izx(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
C-----------------------------------------------
      DO I=jft,jlt
        IF (igrth(i) /= igrth(i+1)) THEN
          IF(off(i) /= zero)THEN
            DO j = igrth(i),igrth(i+1)-1
              gresav(1,grth(j)) = gresav(1,grth(j)) + ei(i)
              gresav(2,grth(j)) = gresav(2,grth(j)) + ek(i)
              gresav(3,grth(j)) = gresav(3,grth(j)) + xm(i)
              gresav(4,grth(j)) = gresav(4,grth(j)) + ym(i)
              gresav(5,grth(j)) = gresav(5,grth(j)) + zm(i)
              gresav(6,grth(j)) = gresav(6,grth(j)) + xmas(i)
            ENDDO
          ENDIF
        ENDIF
      ENDDO
      IF (flag /= 0)THEN
        DO i=jft,jlt
          IF (igrth(i) /= igrth(i+1)) THEN
            DO j = igrth(i),igrth(i+1)-1
              gresav(9,grth(j))  = gresav(9,grth(j))  + xcg(i)
              gresav(10,grth(j)) = gresav(10,grth(j)) + ycg(i)
              gresav(11,grth(j)) = gresav(11,grth(j)) + zcg(i)
              gresav(12,grth(j)) = gresav(12,grth(j)) + xxm(i)
              gresav(13,grth(j)) = gresav(13,grth(j)) + yym(i)
              gresav(14,grth(j)) = gresav(14,grth(j)) + zzm(i)
              gresav(15,grth(j)) = gresav(15,grth(j)) + ixx(i)
              gresav(16,grth(j)) = gresav(16,grth(j)) + iyy(i)
              gresav(17,grth(j)) = gresav(17,grth(j)) + izz(i)
              gresav(18,grth(j)) = gresav(18,grth(j)) + ixy(i)
              gresav(19,grth(j)) = gresav(19,grth(j)) + iyz(i)
              gresav(20,grth(j)) = gresav(20,grth(j)) + izx(i)
              gresav(21,grth(j)) = gresav(21,grth(j)) + rei(i)
              gresav(22,grth(j)) = gresav(22,grth(j)) + rek(i)
            ENDDO
          ENDIF
        ENDDO
      ENDIF
      RETURN
      END