s8eselecsh.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/.
!||====================================================================
!||    s8eselecsh   ../engine/source/elements/solid/solide8e/s8eselecsh.F
!||--- called by ------------------------------------------------------
!||    s8eforc3     ../engine/source/elements/solid/solide8e/s8eforc3.F
!||====================================================================
      SUBROUTINE s8eselecsh(
     1   NPTR,    NPTS,    NPTT,    NNPT,
     2   PX1,     PX2,     PX3,     PX4,
     3   PX5,     PX6,     PX7,     PX8,
     4   PY1,     PY2,     PY3,     PY4,
     5   PY5,     PY6,     PY7,     PY8,
     6   PZ1,     PZ2,     PZ3,     PZ4,
     7   PZ5,     PZ6,     PZ7,     PZ8,
     8   P0XY1,   P0XY2,   P0XY3,   P0XY4,
     9   P0XY5,   P0XY6,   P0XY7,   P0XY8,
     A   P0YX1,   P0YX2,   P0YX3,   P0YX4,
     B   P0YX5,   P0YX6,   P0YX7,   P0YX8,
     C   P0XZ1,   P0XZ2,   P0XZ3,   P0XZ4,
     D   P0XZ5,   P0XZ6,   P0XZ7,   P0XZ8,
     E   P0ZX1,   P0ZX2,   P0ZX3,   P0ZX4,
     F   P0ZX5,   P0ZX6,   P0ZX7,   P0ZX8,
     G   P0YZ1,   P0YZ2,   P0YZ3,   P0YZ4,
     H   P0YZ5,   P0YZ6,   P0YZ7,   P0YZ8,
     I   P0ZY1,   P0ZY2,   P0ZY3,   P0ZY4,
     J   P0ZY5,   P0ZY6,   P0ZY7,   P0ZY8,
     K   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   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
      INTEGER NPTR,NPTS,NPTT,NNPT 
C     REAL
      my_real
     .  PX1(MVSIZ,NNPT),PX2(MVSIZ,NNPT),PX3(MVSIZ,NNPT),PX4(MVSIZ,NNPT),
     .  PX5(MVSIZ,NNPT),PX6(MVSIZ,NNPT),PX7(MVSIZ,NNPT),PX8(MVSIZ,NNPT),
     .  PY1(MVSIZ,NNPT),PY2(MVSIZ,NNPT),PY3(MVSIZ,NNPT),PY4(MVSIZ,NNPT),
     .  PY5(MVSIZ,NNPT),PY6(MVSIZ,NNPT),PY7(MVSIZ,NNPT),PY8(MVSIZ,NNPT),
     .  PZ1(MVSIZ,NNPT),PZ2(MVSIZ,NNPT),PZ3(MVSIZ,NNPT),PZ4(MVSIZ,NNPT),
     .  PZ5(MVSIZ,NNPT),PZ6(MVSIZ,NNPT),PZ7(MVSIZ,NNPT),PZ8(MVSIZ,NNPT),
     .   P0XY1(MVSIZ,2),P0XY2(MVSIZ,2),P0XY3(MVSIZ,2),P0XY4(MVSIZ,2),
     .   P0XY5(MVSIZ,2),P0XY6(MVSIZ,2),P0XY7(MVSIZ,2),P0XY8(MVSIZ,2),
     .   P0YX1(MVSIZ,2),P0YX2(MVSIZ,2),P0YX3(MVSIZ,2),P0YX4(MVSIZ,2),
     .   P0YX5(MVSIZ,2),P0YX6(MVSIZ,2),P0YX7(MVSIZ,2),P0YX8(MVSIZ,2),
     .   P0XZ1(MVSIZ,2),P0XZ2(MVSIZ,2),P0XZ3(MVSIZ,2),P0XZ4(MVSIZ,2),
     .   P0XZ5(MVSIZ,2),P0XZ6(MVSIZ,2),P0XZ7(MVSIZ,2),P0XZ8(MVSIZ,2),
     .   P0ZX1(MVSIZ,2),P0ZX2(MVSIZ,2),P0ZX3(MVSIZ,2),P0ZX4(MVSIZ,2),
     .   P0ZX5(MVSIZ,2),P0ZX6(MVSIZ,2),P0ZX7(MVSIZ,2),P0ZX8(MVSIZ,2),
     .   P0YZ1(MVSIZ,2),P0YZ2(MVSIZ,2),P0YZ3(MVSIZ,2),P0YZ4(MVSIZ,2),
     .   P0YZ5(MVSIZ,2),P0YZ6(MVSIZ,2),P0YZ7(MVSIZ,2),P0YZ8(MVSIZ,2),
     .   p0zy1(mvsiz,2),p0zy2(mvsiz,2),p0zy3(mvsiz,2),p0zy4(mvsiz,2),
     .   p0zy5(mvsiz,2),p0zy6(mvsiz,2),p0zy7(mvsiz,2),p0zy8(mvsiz,2)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J ,IR,IS,IT,IP,I1,I2,I3,I4
      INTEGER IXY(4,2),IXZ(4,2),IYZ(4,2)
      DATA IXY/1,5,3,7,2,6,4,8/
      DATA IXZ/1,3,2,4,5,7,6,8/
      DATA IYZ/1,5,2,6,3,7,4,8/
C------ moyen in 3 directions
      DO IR=1,nptr
        i1 = ixy(1,ir)
        i2 = ixy(2,ir)
        i3 = ixy(3,ir)
        i4 = ixy(4,ir)
       DO i=1,nel
       p0yx1(i,ir)=fourth*(px1(i,i1)+px1(i,i2)+px1(i,i3)+px1(i,i4))
       p0yx2(i,ir)=fourth*(px2(i,i1)+px2(i,i2)+px2(i,i3)+px2(i,i4))
       p0yx3(i,ir)=fourth*(px3(i,i1)+px3(i,i2)+px3(i,i3)+px3(i,i4))
       p0yx4(i,ir)=fourth*(px4(i,i1)+px4(i,i2)+px4(i,i3)+px4(i,i4))
       p0yx5(i,ir)=fourth*(px5(i,i1)+px5(i,i2)+px5(i,i3)+px5(i,i4))
       p0yx6(i,ir)=fourth*(px6(i,i1)+px6(i,i2)+px6(i,i3)+px6(i,i4))
       p0yx7(i,ir)=fourth*(px7(i,i1)+px7(i,i2)+px7(i,i3)+px7(i,i4))
       p0yx8(i,ir)=fourth*(px8(i,i1)+px8(i,i2)+px8(i,i3)+px8(i,i4))
       ENDDO
       DO i=1,nel
       p0xy1(i,ir)=fourth*(py1(i,i1)+py1(i,i2)+py1(i,i3)+py1(i,i4))
       p0xy2(i,ir)=fourth*(py2(i,i1)+py2(i,i2)+py2(i,i3)+py2(i,i4))
       p0xy3(i,ir)=fourth*(py3(i,i1)+py3(i,i2)+py3(i,i3)+py3(i,i4))
       p0xy4(i,ir)=fourth*(py4(i,i1)+py4(i,i2)+py4(i,i3)+py4(i,i4))
       p0xy5(i,ir)=fourth*(py5(i,i1)+py5(i,i2)+py5(i,i3)+py5(i,i4))
       p0xy6(i,ir)=fourth*(py6(i,i1)+py6(i,i2)+py6(i,i3)+py6(i,i4))
       p0xy7(i,ir)=fourth*(py7(i,i1)+py7(i,i2)+py7(i,i3)+py7(i,i4))
       p0xy8(i,ir)=fourth*(py8(i,i1)+py8(i,i2)+py8(i,i3)+py8(i,i4))
       ENDDO
      ENDDO
C      
      DO it=1,nptt
        i1 = ixz(1,it)
        i2 = ixz(2,it)
        i3 = ixz(3,it)
        i4 = ixz(4,it)
       DO i=1,nel
       p0zx1(i,it)=fourth*(px1(i,i1)+px1(i,i2)+px1(i,i3)+px1(i,i4))
       p0zx2(i,it)=fourth*(px2(i,i1)+px2(i,i2)+px2(i,i3)+px2(i,i4))
       p0zx3(i,it)=fourth*(px3(i,i1)+px3(i,i2)+px3(i,i3)+px3(i,i4))
       p0zx4(i,it)=fourth*(px4(i,i1)+px4(i,i2)+px4(i,i3)+px4(i,i4))
       p0zx5(i,it)=fourth*(px5(i,i1)+px5(i,i2)+px5(i,i3)+px5(i,i4))
       p0zx6(i,it)=fourth*(px6(i,i1)+px6(i,i2)+px6(i,i3)+px6(i,i4))
       p0zx7(i,it)=fourth*(px7(i,i1)+px7(i,i2)+px7(i,i3)+px7(i,i4))
       p0zx8(i,it)=fourth*(px8(i,i1)+px8(i,i2)+px8(i,i3)+px8(i,i4))
       ENDDO
       DO i=1,nel
       p0xz1(i,it)=fourth*(pz1(i,i1)+pz1(i,i2)+pz1(i,i3)+pz1(i,i4))
       p0xz2(i,it)=fourth*(pz2(i,i1)+pz2(i,i2)+pz2(i,i3)+pz2(i,i4))
       p0xz3(i,it)=fourth*(pz3(i,i1)+pz3(i,i2)+pz3(i,i3)+pz3(i,i4))
       p0xz4(i,it)=fourth*(pz4(i,i1)+pz4(i,i2)+pz4(i,i3)+pz4(i,i4))
       p0xz5(i,it)=fourth*(pz5(i,i1)+pz5(i,i2)+pz5(i,i3)+pz5(i,i4))
       p0xz6(i,it)=fourth*(pz6(i,i1)+pz6(i,i2)+pz6(i,i3)+pz6(i,i4))
       p0xz7(i,it)=fourth*(pz7(i,i1)+pz7(i,i2)+pz7(i,i3)+pz7(i,i4))
       p0xz8(i,it)=fourth*(pz8(i,i1)+pz8(i,i2)+pz8(i,i3)+pz8(i,i4))
       ENDDO
      ENDDO
C     
      DO is=1,npts
        i1 = iyz(1,is)
        i2 = iyz(2,is)
        i3 = iyz(3,is)
        i4 = iyz(4,is)
       DO i=1,nel
       p0zy1(i,is)=fourth*(py1(i,i1)+py1(i,i2)+py1(i,i3)+py1(i,i4))
       p0zy2(i,is)=fourth*(py2(i,i1)+py2(i,i2)+py2(i,i3)+py2(i,i4))
       p0zy3(i,is)=fourth*(py3(i,i1)+py3(i,i2)+py3(i,i3)+py3(i,i4))
       p0zy4(i,is)=fourth*(py4(i,i1)+py4(i,i2)+py4(i,i3)+py4(i,i4))
       p0zy5(i,is)=fourth*(py5(i,i1)+py5(i,i2)+py5(i,i3)+py5(i,i4))
       p0zy6(i,is)=fourth*(py6(i,i1)+py6(i,i2)+py6(i,i3)+py6(i,i4))
       p0zy7(i,is)=fourth*(py7(i,i1)+py7(i,i2)+py7(i,i3)+py7(i,i4))
       p0zy8(i,is)=fourth*(py8(i,i1)+py8(i,i2)+py8(i,i3)+py8(i,i4))
       ENDDO
       DO i=1,nel
       p0yz1(i,is)=fourth*(pz1(i,i1)+pz1(i,i2)+pz1(i,i3)+pz1(i,i4))
       p0yz2(i,is)=fourth*(pz2(i,i1)+pz2(i,i2)+pz2(i,i3)+pz2(i,i4))
       p0yz3(i,is)=fourth*(pz3(i,i1)+pz3(i,i2)+pz3(i,i3)+pz3(i,i4))
       p0yz4(i,is)=fourth*(pz4(i,i1)+pz4(i,i2)+pz4(i,i3)+pz4(i,i4))
       p0yz5(i,is)=fourth*(pz5(i,i1)+pz5(i,i2)+pz5(i,i3)+pz5(i,i4))
       p0yz6(i,is)=fourth*(pz6(i,i1)+pz6(i,i2)+pz6(i,i3)+pz6(i,i4))
       p0yz7(i,is)=fourth*(pz7(i,i1)+pz7(i,i2)+pz7(i,i3)+pz7(i,i4))
       p0yz8(i,is)=fourth*(pz8(i,i1)+pz8(i,i2)+pz8(i,i3)+pz8(i,i4))
       ENDDO
      ENDDO
C
      RETURN
       END