s4derito3.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/.
!||====================================================================
!||    s4derito3   ../engine/source/elements/solid/solide4/s4derito3.F
!||--- called by ------------------------------------------------------
!||    s4forc3     ../engine/source/elements/solid/solide4/s4forc3.F
!||====================================================================
      SUBROUTINE s4derito3(
     1   OFF,     DET,     X1,      X2,
     2   X3,      X4,      Y1,      Y2,
     3   Y3,      Y4,      Z1,      Z2,
     4   Z3,      Z4,      PX1,     PX2,
     5   PX3,     PX4,     PY1,     PY2,
     6   PY3,     PY4,     PZ1,     PZ2,
     7   PZ3,     PZ4,     RX,      RY,
     8   RZ,      SX,      SY,      SZ,
     9   TX,      TY,      TZ,      JAC_I,
     A   NEL)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   OFF(*),DET(*),
     .   PX1(*), PX2(*), PX3(*), PX4(*),  
     .   PY1(*), PY2(*), PY3(*), PY4(*),  
     .   PZ1(*), PZ2(*), PZ3(*), PZ4(*), 
     .   RX(*), RY(*), RZ(*), SX(*), SY(*), SZ(*),TX(*), TY(*), TZ(*),
     .   JAC_I(10,*)
      DOUBLE PRECISION
     .   X1(*), X2(*), X3(*), X4(*),
     .   y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*)
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     REAL
C                                                                     12
      my_real
     .   X41, Y41, Z41, X42, Y42, Z42, X43, Y43, Z43
      my_real
     .   JI1, JI2, JI3,  
     .   JI4, JI5, JI6,  
     .   JI7, JI8, JI9,  
     .   A1, A2, A3, A4, D, AREAMX2
C-----------------------------------------------
      DO i=1,nel
       x43 = x4(i) - x3(i)
       y43 = y4(i) - y3(i)
       z43 = z4(i) - z3(i)
       x41 = x4(i) - x1(i)
       y41 = y4(i) - y1(i)
       z41 = z4(i) - z1(i)
       x42 = x4(i) - x2(i)
       y42 = y4(i) - y2(i)
       z42 = z4(i) - z2(i)
C
       rx(i) =  -x41
       ry(i) =  -y41
       rz(i) =  -z41
       sx(i) =  -x42
       sy(i) =  -y42
       sz(i) =  -z42
       tx(i) =  -x43
       ty(i) =  -y43
       tz(i) =  -z43
C
      ENDDO
C
      DO i=1,nel
       ji1 = jac_i(1,i)
       ji2 = jac_i(2,i)
       ji3 = jac_i(3,i)
       ji4 = jac_i(4,i)
       ji5 = jac_i(5,i)
       ji6 = jac_i(6,i)
       ji7 = jac_i(7,i)
       ji8 = jac_i(8,i)
       ji9 = jac_i(9,i)
C
        px1(i)= ji1
        py1(i)= ji4
        pz1(i)= ji7
        px2(i)= ji2
        py2(i)= ji5
        pz2(i)= ji8
        px3(i)= ji3
        py3(i)= ji6
        pz3(i)= ji9
        px4(i)=-(ji1+ji2+ji3)
        py4(i)=-(ji4+ji5+ji6)
        pz4(i)=-(ji7+ji8+ji9)
      ENDDO
C
      RETURN
C
      END