s10jacobt.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/.
!||====================================================================
!||    s10jacobt   ../engine/source/elements/solid/solide10/s10jacobt.F
!||====================================================================
      SUBROUTINE s10jacobt(
     1   ALPH,    BETA,    W,       PX1,
     2   PX2,     PX3,     PX4,     PX5,
     3   PX6,     PX7,     PX8,     PX9,
     4   PX10,    PY1,     PY2,     PY3,
     5   PY4,     PY5,     PY6,     PY7,
     6   PY8,     PY9,     PY10,    PZ1,
     7   PZ2,     PZ3,     PZ4,     PZ5,
     8   PZ6,     PZ7,     PZ8,     PZ9,
     9   PZ10,    NX1,     NX2,     NX3,
     A   NX4,     NX5,     NX6,     NX7,
     B   NX8,     NX9,     NX10,    VOL,
     C   JAC_I,   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, INTENT(IN) :: NEL
       my_real
     .   PX1(MVSIZ),PX2(MVSIZ),PX3(MVSIZ),PX4(MVSIZ),PX5(MVSIZ),
     .   PX6(MVSIZ),PX7(MVSIZ),PX8(MVSIZ),PX9(MVSIZ),PX10(MVSIZ),
     .   PY1(MVSIZ),PY2(MVSIZ),PY3(MVSIZ),PY4(MVSIZ),PY5(MVSIZ),
     .   PY6(MVSIZ),PY7(MVSIZ),PY8(MVSIZ),PY9(MVSIZ),PY10(MVSIZ),
     .   PZ1(MVSIZ),PZ2(MVSIZ),PZ3(MVSIZ),PZ4(MVSIZ),PZ5(MVSIZ),
     .   PZ6(MVSIZ),PZ7(MVSIZ),PZ8(MVSIZ),PZ9(MVSIZ),PZ10(MVSIZ),
     .   NX1(MVSIZ),NX2(MVSIZ),NX3(MVSIZ),NX4(MVSIZ),NX5(MVSIZ),
     .   NX6(MVSIZ),NX7(MVSIZ),NX8(MVSIZ),NX9(MVSIZ),NX10(MVSIZ),
     .   VOL(MVSIZ),ALPH,BETA,W,JAC_I(10,MVSIZ)
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
      my_real
     .   D,A4MB4,A4,B4,BB,AB,B2,
     .   DXDR,DXDS,DXDT,DYDR,DYDS,DYDT,DZDR,DZDS,DZDT
     
      DOUBLE PRECISION 
     .   AA,A4M1,B4M1
      
      my_real
     .   DRDX, DSDX, DTDX,
     .   DRDY, DSDY, DTDY,
     .   DRDZ, DSDZ, DTDZ
C-----------------------------------------------
c      AA = (TWO*ALPH - ONE)*ALPH
c      BB = (TWO*BETA - ONE)*BETA
c      B2 = FOUR*BETA*BETA
c      AB = FOUR*ALPH*BETA      
C-----------------------------------------------
C
      a4 = four * alph
      b4 = four * beta
      a4m1  = a4 - one
      b4m1  = b4 - one
      a4mb4  = a4 - b4
      DO i=1,nel
       drdx=jac_i(1,i)
       dsdx=jac_i(4,i)
       dtdx=jac_i(7,i)
C
       drdy=jac_i(2,i)
       dsdy=jac_i(5,i)
       dtdy=jac_i(8,i)
C
       drdz=jac_i(3,i)
       dsdz=jac_i(6,i)
       dtdz=jac_i(9,i)
C
       vol(i) = jac_i(10,i)
C
        px1(i) = b4m1 * drdx
        py1(i) = b4m1 * drdy
        pz1(i) = b4m1 * drdz
C
        px2(i) = b4m1 * dsdx
        py2(i) = b4m1 * dsdy
        pz2(i) = b4m1 * dsdz
C
        px3(i) = b4m1 * dtdx
        py3(i) = b4m1 * dtdy
        pz3(i) = b4m1 * dtdz
C
        px4(i) =-a4m1 * (drdx+dsdx+dtdx)
        py4(i) =-a4m1 * (drdy+dsdy+dtdy)
        pz4(i) =-a4m1 * (drdz+dsdz+dtdz)
C
        px5(i) = b4 * (drdx+dsdx)
        py5(i) = b4 * (drdy+dsdy)
        pz5(i) = b4 * (drdz+dsdz)
C
        px6(i) = b4 * (dsdx+dtdx)
        py6(i) = b4 * (dsdy+dtdy)
        pz6(i) = b4 * (dsdz+dtdz)
C
        px7(i) = b4 * (dtdx+drdx)
        py7(i) = b4 * (dtdy+drdy)
        pz7(i) = b4 * (dtdz+drdz)
C
        px8(i) = a4mb4  * drdx - px6(i)
        py8(i) = a4mb4  * drdy - py6(i)
        pz8(i) = a4mb4  * drdz - pz6(i)
C
        px9(i) = a4mb4  * dsdx - px7(i)
        py9(i) = a4mb4  * dsdy - py7(i)
        pz9(i) = a4mb4  * dsdz - pz7(i)
C
        px10(i)= a4mb4  * dtdx - px5(i)
        py10(i)= a4mb4  * dtdy - py5(i)
        pz10(i)= a4mb4  * dtdz - pz5(i)
C
      ENDDO
C
      RETURN
      END