fe_close.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/.
!||====================================================================
!||    fe_close   ../engine/source/elements/solid/solide/fe_close.F
!||--- called by ------------------------------------------------------
!||    sforc3     ../engine/source/elements/solid/solide/sforc3.F
!||    szforc3    ../engine/source/elements/solid/solidez/szforc3.F
!||====================================================================
      SUBROUTINE fe_close(
     1   HH,      RHO,     VOLN,    VX1,
     2   VX2,     VX3,     VX4,     VX5,
     3   VX6,     VX7,     VX8,     VY1,
     4   VY2,     VY3,     VY4,     VY5,
     5   VY6,     VY7,     VY8,     VZ1,
     6   VZ2,     VZ3,     VZ4,     VZ5,
     7   VZ6,     VZ7,     VZ8,     F11,
     8   F21,     F31,     F12,     F22,
     9   F32,     F13,     F23,     F33,
     A   F14,     F24,     F34,     F15,
     B   F25,     F35,     F16,     F26,
     C   F36,     F17,     F27,     F37,
     D   F18,     F28,     F38,     NEL)
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      "mvsiz_p.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   F11(MVSIZ),F21(MVSIZ),F31(MVSIZ),
     .   F12(MVSIZ),F22(MVSIZ),F32(MVSIZ),
     .   F13(MVSIZ),F23(MVSIZ),F33(MVSIZ),
     .   F14(MVSIZ),F24(MVSIZ),F34(MVSIZ),
     .   F15(MVSIZ),F25(MVSIZ),F35(MVSIZ),
     .   F16(MVSIZ),F26(MVSIZ),F36(MVSIZ),
     .   F17(MVSIZ),F27(MVSIZ),F37(MVSIZ),
     .   F18(MVSIZ),F28(MVSIZ),F38(MVSIZ),
     .   HH(*),RHO(*),VOLN(*),
     .   VX1(*),VX2(*),VX3(*),VX4(*),VX5(*),VX6(*),VX7(*),VX8(*),
     .   VY1(*),VY2(*),VY3(*),VY4(*),VY5(*),VY6(*),VY7(*),VY8(*),
     .   vz1(*),vz2(*),vz3(*),vz4(*),vz5(*),vz6(*),vz7(*),vz8(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     . FAC  
C      
      IF(DT1==zero)RETURN
      DO i=1,nel
        fac=-one_over_8*hh(i)*rho(i)*voln(i)/dt1
          f11(i) = f11(i) + fac*vx1(i)
          f12(i) = f12(i) + fac*vx2(i)
          f13(i) = f13(i) + fac*vx3(i)
          f14(i) = f14(i) + fac*vx4(i)
          f15(i) = f15(i) + fac*vx5(i)
          f16(i) = f16(i) + fac*vx6(i)
          f17(i) = f17(i) + fac*vx7(i)
          f18(i) = f18(i) + fac*vx8(i)
 
          f21(i) = f21(i) + fac*vy1(i)
          f22(i) = f22(i) + fac*vy2(i)
          f23(i) = f23(i) + fac*vy3(i)
          f24(i) = f24(i) + fac*vy4(i)
          f25(i) = f25(i) + fac*vy5(i)
          f26(i) = f26(i) + fac*vy6(i)
          f27(i) = f27(i) + fac*vy7(i)
          f28(i) = f28(i) + fac*vy8(i)
 
          f31(i) = f31(i) + fac*vz1(i)
          f32(i) = f32(i) + fac*vz2(i)
          f33(i) = f33(i) + fac*vz3(i)
          f34(i) = f34(i) + fac*vz4(i)
          f35(i) = f35(i) + fac*vz5(i)
          f36(i) = f36(i) + fac*vz6(i)
          f37(i) = f37(i) + fac*vz7(i)
          f38(i) = f38(i) + fac*vz8(i)
      ENDDO
      RETURN
      END