volint.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/.
!||====================================================================
!||    volint     ../starter/source/interfaces/inter3d1/volint.F
!||--- called by ------------------------------------------------------
!||    i11sti3    ../starter/source/interfaces/inter3d1/i11sti3.F
!||    i20sti3    ../starter/source/interfaces/inter3d1/i20sti3.F
!||    i20sti3e   ../starter/source/interfaces/inter3d1/i20sti3.F
!||    i24gapm    ../starter/source/interfaces/inter3d1/i24sti3.F
!||    i24ll_kg   ../starter/source/interfaces/inter3d1/i24sti3.F
!||    i25gapm    ../starter/source/interfaces/inter3d1/i25sti3.F
!||    i2buc1     ../starter/source/interfaces/inter3d1/i2buc1.F
!||    i2cor3     ../starter/source/interfaces/inter3d1/i2cor3.F
!||    i3sti3     ../starter/source/interfaces/inter3d1/i3sti3.F
!||    i7sti3     ../starter/source/interfaces/inter3d1/i7sti3.F
!||====================================================================
      SUBROUTINE volint(VOL)
C-----------------------------------------------
C   COMPUTE HEXAHEDRON VOLUME
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      "scr08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real, INTENT(OUT) :: vol
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real x17, x28, x35, x46,
     .        y17, y28, y35, y46,
     .        z17, z28, z35, z46
      my_real a17, a28, b17, b28, c17, c28 
      my_real jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9 
      my_real jac_59_68, jac_67_49, jac_48_57 
C-------------------------------------------------------
      x17 = xc(7)-xc(1)
      x28 = xc(8)-xc(2)
      x35 = xc(5)-xc(3)
      x46 = xc(6)-xc(4)
C
      y17 = yc(7)-yc(1)
      y28 = yc(8)-yc(2)
      y35 = yc(5)-yc(3)
      y46 = yc(6)-yc(4)
C
      z17 = zc(7)-zc(1)
      z28 = zc(8)-zc(2)
      z35 = zc(5)-zc(3)
      z46 = zc(6)-zc(4)
C
      jac4 = x17 + x28 - x35 - x46 
      jac5 = y17 + y28 - y35 - y46 
      jac6 = z17 + z28 - z35 - z46 
C
      a17 = x17 + x46 
      a28 = x28 + x35 
      b17 = y17 + y46 
      b28 = y28 + y35 
      c17 = z17 + z46 
      c28 = z28 + z35 
C
      jac7 = a17 + a28 
      jac8 = b17 + b28 
      jac9 = c17 + c28 
      jac1 = a17 - a28 
      jac2 = b17 - b28 
      jac3 = c17 - c28 
C
      jac_59_68 = jac5 * jac9 - jac6 * jac8 
      jac_67_49 = jac6 * jac7 - jac4 * jac9 
      jac_48_57 = jac4 * jac8 - jac5 * jac7 
C
      vol =  zep015625 * (jac1 * jac_59_68 + jac2 * jac_67_49 + jac3 * jac_48_57)
 
      RETURN
      END