s8cumu3.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/.
!||====================================================================
!||    s8cumu3   ../engine/source/elements/solid/solide8/s8cumu3.F
!||--- called by ------------------------------------------------------
!||    s8forc3   ../engine/source/elements/solid/solide8/s8forc3.F
!||====================================================================
      SUBROUTINE s8cumu3(
     1   OFFG,    E,       NC,      STIFN,
     2   STI,     F11,     F21,     F31,
     3   F12,     F22,     F32,     F13,
     4   F23,     F33,     F14,     F24,
     5   F34,     F15,     F25,     F35,
     6   F16,     F26,     F36,     F17,
     7   F27,     F37,     F18,     F28,
     8   F38,     NEL)
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, INTENT(IN) :: NEL
      INTEGER NC(8,*)
      my_real
     .   OFFG(*),E(3,*),STIFN(*),STI(*),
     .   F11(*),F21(*),F31(*),F12(*),F22(*),F32(*),
     .   F13(*),F23(*),F33(*),F14(*),F24(*),F34(*),
     .   F15(*),F25(*),F35(*),F16(*),F26(*),F36(*),
     .   F17(*),F27(*),F37(*),F18(*),F28(*),F38(*)
CMasParINCLUDE 'scumu3.intmap.inc'
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
C-----------------------------------------------
      my_real
     .   OFF_L
C-----------------------------------------------
      OFF_L = 0.
      do i=1,nel
Cf small3b        IF(OFF(I)<1.)OFFG(I) = OFF(I)
        off_l = min(off_l,offg(i))
      ENDDO
      IF(off_l<0.)THEN
        DO i=1,nel
          IF(offg(i)<0.)THEN
            f11(i)=0.
            f21(i)=0.
            f31(i)=0.
            f12(i)=0.
            f22(i)=0.
            f32(i)=0.
            f13(i)=0.
            f23(i)=0.
            f33(i)=0.
            f14(i)=0.
            f24(i)=0.
            f34(i)=0.
            f15(i)=0.
            f25(i)=0.
            f35(i)=0.
            f16(i)=0.
            f26(i)=0.
            f36(i)=0.
            f17(i)=0.
            f27(i)=0.
            f37(i)=0.
            f18(i)=0.
            f28(i)=0.
            f38(i)=0.
            sti(i)=0.
          ENDIF
        ENDDO
      ENDIF
C
C     because 2*Mnodal = 2*Melement/8
      DO i=1,nel
        sti(i)=fourth*sti(i)
      END DO
C
      DO 100 i=1,nel
C
      e(1,nc(1,i))=e(1,nc(1,i))+f11(i)
      e(2,nc(1,i))=e(2,nc(1,i))+f21(i)
      e(3,nc(1,i))=e(3,nc(1,i))+f31(i)
      stifn(nc(1,i))=stifn(nc(1,i))+sti(i)
C
      e(1,nc(7,i))=e(1,nc(7,i))+f17(i)
      e(2,nc(7,i))=e(2,nc(7,i))+f27(i)
      e(3,nc(7,i))=e(3,nc(7,i))+f37(i)
      stifn(nc(7,i))=stifn(nc(7,i))+sti(i)
C
      e(1,nc(2,i))=e(1,nc(2,i))+f12(i)
      e(2,nc(2,i))=e(2,nc(2,i))+f22(i)
      e(3,nc(2,i))=e(3,nc(2,i))+f32(i)
      stifn(nc(2,i))=stifn(nc(2,i))+sti(i)
C
      e(1,nc(8,i))=e(1,nc(8,i))+f18(i)
      e(2,nc(8,i))=e(2,nc(8,i))+f28(i)
      e(3,nc(8,i))=e(3,nc(8,i))+f38(i)
      stifn(nc(8,i))=stifn(nc(8,i))+sti(i)
C
      e(1,nc(3,i))=e(1,nc(3,i))+f13(i)
      e(2,nc(3,i))=e(2,nc(3,i))+f23(i)
      e(3,nc(3,i))=e(3,nc(3,i))+f33(i)
      stifn(nc(3,i))=stifn(nc(3,i))+sti(i)
C
      e(1,nc(5,i))=e(1,nc(5,i))+f15(i)
      e(2,nc(5,i))=e(2,nc(5,i))+f25(i)
      e(3,nc(5,i))=e(3,nc(5,i))+f35(i)
      stifn(nc(5,i))=stifn(nc(5,i))+sti(i)
C
      e(1,nc(4,i))=e(1,nc(4,i))+f14(i)
      e(2,nc(4,i))=e(2,nc(4,i))+f24(i)
      e(3,nc(4,i))=e(3,nc(4,i))+f34(i)
      stifn(nc(4,i))=stifn(nc(4,i))+sti(i)
C
      e(1,nc(6,i))=e(1,nc(6,i))+f16(i)
      e(2,nc(6,i))=e(2,nc(6,i))+f26(i)
      e(3,nc(6,i))=e(3,nc(6,i))+f36(i)
      stifn(nc(6,i))=stifn(nc(6,i))+sti(i)
C
  100 CONTINUE
      RETURN
      END