s10volj.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/.
!||====================================================================
!||    s10volj      ../engine/source/elements/solid/solide4_sfem/s10volj.F
!||--- called by ------------------------------------------------------
!||    s10volnod3   ../engine/source/elements/solid/solide4_sfem/s10volnod3.F
!||====================================================================
      SUBROUTINE s10volj(W,
     .   X1B,X2B,X3B,X4A,X5B,X6B,X7B,X8B,X9B,X10B,X8A,X9A,X10A,
     .   Y1B,Y2B,Y3B,Y4A,Y5B,Y6B,Y7B,Y8B,Y9B,Y10B,Y8A,Y9A,Y10A,
     .   Z1B,Z2B,Z3B,Z4A,Z5B,Z6B,Z7B,Z8B,Z9B,Z10B,Z8A,Z9A,Z10A,
     .   VOLDP,NEL)
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
C     REAL
      DOUBLE PRECISION
     .   X1B(MVSIZ),X2B(MVSIZ),X3B(MVSIZ),X4A(MVSIZ),X5B(MVSIZ),
     .   x6b(mvsiz),x7b(mvsiz),x8b(mvsiz),x9b(mvsiz),x10b(mvsiz),
     .   x8a(mvsiz),x9a(mvsiz),x10a(mvsiz),
     .   y1b(mvsiz),y2b(mvsiz),y3b(mvsiz),y4a(mvsiz),y5b(mvsiz),
     .   y6b(mvsiz),y7b(mvsiz),y8b(mvsiz),y9b(mvsiz),y10b(mvsiz),
     .   y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
     .   z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
     .   z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
     .   z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(mvsiz)
      my_real
     .   w
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      DOUBLE PRECISION
     .   A4M1,B4M1,AA
      DOUBLE PRECISION
     .   d,a4mb4,a4,b4,bb,ab,b2,
     .   drdx, drdy, drdz,det,
     .   dxdr,dxds,dxdt,dydr,dyds,dydt,dzdr,dzds,dzdt
C-----------------------------------------------
C
      DO i=1,nel
       aa = x5b(i) + x6b(i) + x7b(i)
     .    - x4a(i) - x8b(i) - x9b(i) - x10b(i)
       dxdr = x1b(i) +x8a(i) - x6b(i) + aa
       dxds = x2b(i) +x9a(i) - x7b(i) + aa
       dxdt = x3b(i) +x10a(i)- x5b(i) + aa
C
       aa = y5b(i) + y6b(i) + y7b(i)
     .    - y4a(i) - y8b(i) - y9b(i) - y10b(i)
       dydr = y1b(i) +y8a(i) - y6b(i) + aa
       dyds = y2b(i) +y9a(i) - y7b(i) + aa
       dydt = y3b(i) +y10a(i)- y5b(i) + aa
C
       aa = z5b(i) + z6b(i) + z7b(i)
     .    - z4a(i) - z8b(i) - z9b(i) - z10b(i)
       dzdr = z1b(i) +z8a(i) - z6b(i) + aa
       dzds = z2b(i) +z9a(i) - z7b(i) + aa
       dzdt = z3b(i) +z10a(i)- z5b(i) + aa
C
       drdx =dyds*dzdt-dzds*dydt
       drdy =dzds*dxdt-dxds*dzdt
       drdz =dxds*dydt-dyds*dxdt
C
       det = dxdr * drdx + dydr * drdy + dzdr * drdz
       voldp(i) = w * det * one_over_6
C
      ENDDO
C
      RETURN
      END