s8edericto3.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/.
!||====================================================================
!||    s8edericto3   ../engine/source/elements/solid/solide8e/s8edericto3.F
!||--- called by ------------------------------------------------------
!||    s8eforc3      ../engine/source/elements/solid/solide8e/s8eforc3.F
!||====================================================================
      SUBROUTINE s8edericto3(
     1   OFF,     DET,     PX1,     PX2,
     2   PX3,     PX4,     PY1,     PY2,
     3   PY3,     PY4,     PZ1,     PZ2,
     4   PZ3,     PZ4,     JAC_I,   NEL)
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
C   C o m m o n   B l o c k 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   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
C     REAL   
      my_real
     .   OFF(*),DET(*),  
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   jac_i(10,mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J ,ICOR
C     REAL
C                                                                     12
      my_real
     .   aji1(mvsiz), aji2(mvsiz), aji3(mvsiz),
     .   aji4(mvsiz), aji5(mvsiz), aji6(mvsiz),
     .   aji7(mvsiz), aji8(mvsiz), aji9(mvsiz),
     .   aj12(mvsiz), aj45(mvsiz), aj78(mvsiz)
C-----------------------------------------------
C  
      DO i=1,nel
      aji1(i)=jac_i(1,i)
      aji4(i)=jac_i(4,i)
      aji7(i)=jac_i(7,i)
      aji2(i)=jac_i(2,i)
      aji5(i)=jac_i(5,i)
      aji8(i)=jac_i(8,i)
      aji3(i)=jac_i(3,i)
      aji6(i)=jac_i(6,i)
      aji9(i)=jac_i(9,i)
      det(i) =jac_i(10,i)
      END DO
C
      DO 160 i=1,nel
      aj12(i)=aji1(i)-aji2(i)
      aj45(i)=aji4(i)-aji5(i)
  160 aj78(i)=aji7(i)-aji8(i)
      DO 170 i=1,nel
      px2(i)= aj12(i)-aji3(i)
      py2(i)= aj45(i)-aji6(i)
      pz2(i)= aj78(i)-aji9(i)
      px4(i)=-aj12(i)-aji3(i)
      py4(i)=-aj45(i)-aji6(i)
  170 pz4(i)=-aj78(i)-aji9(i)
      DO 180 i=1,nel
      aj12(i)=aji1(i)+aji2(i)
      aj45(i)=aji4(i)+aji5(i)
  180 aj78(i)=aji7(i)+aji8(i)
      DO 190 i=1,nel
      px1(i)=-aj12(i)-aji3(i)
      py1(i)=-aj45(i)-aji6(i)
      pz1(i)=-aj78(i)-aji9(i)
      px3(i)=aj12(i)-aji3(i)
      py3(i)=aj45(i)-aji6(i)
  190 pz3(i)=aj78(i)-aji9(i)
C       
      RETURN
      END