s8zdefot3_8F_source.html

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/.

!||====================================================================

!||    s8zdefot3   ../engine/source/elements/solid/solide8z/s8zdefot3.F

!||--- called by ------------------------------------------------------

!||    s8zforc3    ../engine/source/elements/solid/solide8z/s8zforc3.F

!||--- calls      -----------------------------------------------------

!||    jacob_j33   ../engine/source/elements/solid/solide8e/jacob_j33.F

!||====================================================================


      SUBROUTINE s8zdefot3(

     1   PX1,     PX2,     PX3,     PX4,

     2   PX5,     PX6,     PX7,     PX8,

     3   PY1,     PY2,     PY3,     PY4,

     4   PY5,     PY6,     PY7,     PY8,

     5   PZ1,     PZ2,     PZ3,     PZ4,

     6   PZ5,     PZ6,     PZ7,     PZ8,

     7   PXY1,    PXY2,    PXY3,    PXY4,

     8   PXY5,    PXY6,    PXY7,    PXY8,

     9   PYX1,    PYX2,    PYX3,    PYX4,

     A   PYX5,    PYX6,    PYX7,    PYX8,

     B   PXZ1,    PXZ2,    PXZ3,    PXZ4,

     C   PXZ5,    PXZ6,    PXZ7,    PXZ8,

     D   PZX1,    PZX2,    PZX3,    PZX4,

     E   PZX5,    PZX6,    PZX7,    PZX8,

     F   PYZ1,    PYZ2,    PYZ3,    PYZ4,

     G   PYZ5,    PYZ6,    PYZ7,    PYZ8,

     H   PZY1,    PZY2,    PZY3,    PZY4,

     I   PZY5,    PZY6,    PZY7,    PZY8,

     J   BXY1,    BXY2,    BXY3,    BXY4,

     K   BXY5,    BXY6,    BXY7,    BXY8,

     L   BYX1,    BYX2,    BYX3,    BYX4,

     M   BYX5,    BYX6,    BYX7,    BYX8,

     N   BXZ1,    BXZ2,    BXZ3,    BXZ4,

     O   BXZ5,    BXZ6,    BXZ7,    BXZ8,

     P   BZX1,    BZX2,    BZX3,    BZX4,

     Q   BZX5,    BZX6,    BZX7,    BZX8,

     R   BYZ1,    BYZ2,    BYZ3,    BYZ4,

     S   BYZ5,    BYZ6,    BYZ7,    BYZ8,

     T   BZY1,    BZY2,    BZY3,    BZY4,

     U   BZY5,    BZY6,    BZY7,    BZY8,

     V   VX1,     VX2,     VX3,     VX4,

     W   VX5,     VX6,     VX7,     VX8,

     X   VY1,     VY2,     VY3,     VY4,

     Y   VY5,     VY6,     VY7,     VY8,

     Z   VZ1,     VZ2,     VZ3,     VZ4,

     1   VZ5,     VZ6,     VZ7,     VZ8,

     2   DXX,     DXY,     DXZ,     DYX,

     3   DYY,     DYZ,     DZX,     DZY,

     4   DZZ,     ICP,     DET0,    JFAC,

     5   NU,      IDEG,    NEL)

C                                                                     12

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 ICP,IDEG(*)

C     REAL

      my_real

     .   VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),

     .   VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),

     .   VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),

     .   PX1(*), PX2(*), PX3(*), PX4(*),

     .   PX5(*), PX6(*), PX7(*), PX8(*),

     .   PY1(*), PY2(*), PY3(*), PY4(*),

     .   PY5(*), PY6(*), PY7(*), PY8(*),

     .   PZ1(*), PZ2(*), PZ3(*), PZ4(*),

     .   PZ5(*), PZ6(*), PZ7(*), PZ8(*),

     .   PXY1(*),PXY2(*),PXY3(*),PXY4(*),

     .   PXY5(*),PXY6(*),PXY7(*),PXY8(*),

     .   PYX1(*),PYX2(*),PYX3(*),PYX4(*),

     .   PYX5(*),PYX6(*),PYX7(*),PYX8(*),

     .   PXZ1(*),PXZ2(*),PXZ3(*),PXZ4(*),

     .   PXZ5(*),PXZ6(*),PXZ7(*),PXZ8(*),

     .   PZX1(*),PZX2(*),PZX3(*),PZX4(*),

     .   PZX5(*),PZX6(*),PZX7(*),PZX8(*),

     .   PYZ1(*),PYZ2(*),PYZ3(*),PYZ4(*),

     .   PYZ5(*),PYZ6(*),PYZ7(*),PYZ8(*),

     .   PZY1(*),PZY2(*),PZY3(*),PZY4(*),

     .   PZY5(*),PZY6(*),PZY7(*),PZY8(*),

     .   BXY1(*),BXY2(*),BXY3(*),BXY4(*),

     .   BXY5(*),BXY6(*),BXY7(*),BXY8(*),

     .   BYX1(*),BYX2(*),BYX3(*),BYX4(*),

     .   BYX5(*),BYX6(*),BYX7(*),BYX8(*),

     .   BXZ1(*),BXZ2(*),BXZ3(*),BXZ4(*),

     .   BXZ5(*),BXZ6(*),BXZ7(*),BXZ8(*),

     .   BZX1(*),BZX2(*),BZX3(*),BZX4(*),

     .   BZX5(*),BZX6(*),BZX7(*),BZX8(*),

     .   BYZ1(*),BYZ2(*),BYZ3(*),BYZ4(*),

     .   BYZ5(*),BYZ6(*),BYZ7(*),BYZ8(*),

     .   BZY1(*),BZY2(*),BZY3(*),BZY4(*),

     .   BZY5(*),BZY6(*),BZY7(*),BZY8(*),

     .   DXX(*), DXY(*), DXZ(*),

     .   DYX(*), DYY(*), DYZ(*),

     .   DZX(*), DZY(*), DZZ(*) ,DET0(*),JFAC(*),NU(*)

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER I

C     REAL

      my_real

     .   FAC,DET(MVSIZ),BASE

C-----------------------------------------------

       DXX(1:NEL) =px1(1:nel)*vx1(1:nel)+px2(1:nel)*vx2(1:nel)+px3(1:nel)*vx3(1:nel)+px4(1:nel)*vx4(1:nel)

     .        +px5(1:nel)*vx5(1:nel)+px6(1:nel)*vx6(1:nel)+px7(1:nel)*vx7(1:nel)+px8(1:nel)*vx8(1:nel)

       dyy(1:nel) =py1(1:nel)*vy1(1:nel)+py2(1:nel)*vy2(1:nel)+py3(1:nel)*vy3(1:nel)+py4(1:nel)*vy4(1:nel)

     .        +py5(1:nel)*vy5(1:nel)+py6(1:nel)*vy6(1:nel)+py7(1:nel)*vy7(1:nel)+py8(1:nel)*vy8(1:nel)

       dzz(1:nel) =pz1(1:nel)*vz1(1:nel)+pz2(1:nel)*vz2(1:nel)+pz3(1:nel)*vz3(1:nel)+pz4(1:nel)*vz4(1:nel)

     .        +pz5(1:nel)*vz5(1:nel)+pz6(1:nel)*vz6(1:nel)+pz7(1:nel)*vz7(1:nel)+pz8(1:nel)*vz8(1:nel)

       dxy(1:nel) =pxy1(1:nel)*vx1(1:nel)+pxy2(1:nel)*vx2(1:nel)

     +        +pxy3(1:nel)*vx3(1:nel)+pxy4(1:nel)*vx4(1:nel)

     +        +pxy5(1:nel)*vx5(1:nel)+pxy6(1:nel)*vx6(1:nel)

     +        +pxy7(1:nel)*vx7(1:nel)+pxy8(1:nel)*vx8(1:nel)

       dxz(1:nel) =pxz1(1:nel)*vx1(1:nel)+pxz2(1:nel)*vx2(1:nel)

     +        +pxz3(1:nel)*vx3(1:nel)+pxz4(1:nel)*vx4(1:nel)

     +        +pxz5(1:nel)*vx5(1:nel)+pxz6(1:nel)*vx6(1:nel)

     +        +pxz7(1:nel)*vx7(1:nel)+pxz8(1:nel)*vx8(1:nel)

       dyx(1:nel) =pyx1(1:nel)*vy1(1:nel)+pyx2(1:nel)*vy2(1:nel)

     +        +pyx3(1:nel)*vy3(1:nel)+pyx4(1:nel)*vy4(1:nel)

     +        +pyx5(1:nel)*vy5(1:nel)+pyx6(1:nel)*vy6(1:nel)

     +        +pyx7(1:nel)*vy7(1:nel)+pyx8(1:nel)*vy8(1:nel)

       dyz(1:nel) =pyz1(1:nel)*vy1(1:nel)+pyz2(1:nel)*vy2(1:nel)

     +        +pyz3(1:nel)*vy3(1:nel)+pyz4(1:nel)*vy4(1:nel)

     +        +pyz5(1:nel)*vy5(1:nel)+pyz6(1:nel)*vy6(1:nel)

     +        +pyz7(1:nel)*vy7(1:nel)+pyz8(1:nel)*vy8(1:nel)

       dzx(1:nel) =pzx1(1:nel)*vz1(1:nel)+pzx2(1:nel)*vz2(1:nel)

     +        +pzx3(1:nel)*vz3(1:nel)+pzx4(1:nel)*vz4(1:nel)

     +        +pzx5(1:nel)*vz5(1:nel)+pzx6(1:nel)*vz6(1:nel)

     +        +pzx7(1:nel)*vz7(1:nel)+pzx8(1:nel)*vz8(1:nel)

       dzy(1:nel) =pzy1(1:nel)*vz1(1:nel)+pzy2(1:nel)*vz2(1:nel)

     +        +pzy3(1:nel)*vz3(1:nel)+pzy4(1:nel)*vz4(1:nel)

     +        +pzy5(1:nel)*vz5(1:nel)+pzy6(1:nel)*vz6(1:nel)

     +        +pzy7(1:nel)*vz7(1:nel)+pzy8(1:nel)*vz8(1:nel)

C

      IF (icp==1) THEN

       CALL jacob_j33(

     1   det,     dxx,     dxy,     dxz,

     2   dyx,     dyy,     dyz,     dzx,

     3   dzy,     dzz,     nel)

C-----F'=FAC*F=MF'+[1]-->MF'=FAC*MF+(FAC-ONE)*[1]

        DO i=1,nel

         fac=two_third*nu(i)

         base = det0(i)/max(em20,det(i))

         jfac(i)=exp(fac*log(max(em20,base)))

         IF (ideg(i)>10) jfac(i) =one

        ENDDO

        dxx(1:nel) = jfac(1:nel)*dxx(1:nel)+jfac(1:nel)-one

        dyy(1:nel) = jfac(1:nel)*dyy(1:nel)+jfac(1:nel)-one

        dzz(1:nel) = jfac(1:nel)*dzz(1:nel)+jfac(1:nel)-one

        dxy(1:nel) = jfac(1:nel)*dxy(1:nel)

        dyx(1:nel) = jfac(1:nel)*dyx(1:nel)

        dzx(1:nel) = jfac(1:nel)*dzx(1:nel)

        dxz(1:nel) = jfac(1:nel)*dxz(1:nel)

        dyz(1:nel) = jfac(1:nel)*dyz(1:nel)

        dzy(1:nel) = jfac(1:nel)*dzy(1:nel)

      ELSE

        dyy(1:nel) =dyy(1:nel)+bxy1(1:nel)*vx1(1:nel)+bxy2(1:nel)*vx2(1:nel)

     +        +bxy3(1:nel)*vx3(1:nel)+bxy4(1:nel)*vx4(1:nel)

     +        +bxy5(1:nel)*vx5(1:nel)+bxy6(1:nel)*vx6(1:nel)

     +        +bxy7(1:nel)*vx7(1:nel)+bxy8(1:nel)*vx8(1:nel)

        dzz(1:nel) =dzz(1:nel)+bxz1(1:nel)*vx1(1:nel)+bxz2(1:nel)*vx2(1:nel)

     +        +bxz3(1:nel)*vx3(1:nel)+bxz4(1:nel)*vx4(1:nel)

     +        +bxz5(1:nel)*vx5(1:nel)+bxz6(1:nel)*vx6(1:nel)

     +        +bxz7(1:nel)*vx7(1:nel)+bxz8(1:nel)*vx8(1:nel)

        dxx(1:nel) =dxx(1:nel)+byx1(1:nel)*vy1(1:nel)+byx2(1:nel)*vy2(1:nel)

     +        +byx3(1:nel)*vy3(1:nel)+byx4(1:nel)*vy4(1:nel)

     +        +byx5(1:nel)*vy5(1:nel)+byx6(1:nel)*vy6(1:nel)

     +        +byx7(1:nel)*vy7(1:nel)+byx8(1:nel)*vy8(1:nel)

        dzz(1:nel) =dzz(1:nel)+byz1(1:nel)*vy1(1:nel)+byz2(1:nel)*vy2(1:nel)

     +        +byz3(1:nel)*vy3(1:nel)+byz4(1:nel)*vy4(1:nel)

     +        +byz5(1:nel)*vy5(1:nel)+byz6(1:nel)*vy6(1:nel)

     +        +byz7(1:nel)*vy7(1:nel)+byz8(1:nel)*vy8(1:nel)

        dxx(1:nel) =dxx(1:nel)+bzx1(1:nel)*vz1(1:nel)+bzx2(1:nel)*vz2(1:nel)

     +        +bzx3(1:nel)*vz3(1:nel)+bzx4(1:nel)*vz4(1:nel)

     +        +bzx5(1:nel)*vz5(1:nel)+bzx6(1:nel)*vz6(1:nel)

     +        +bzx7(1:nel)*vz7(1:nel)+bzx8(1:nel)*vz8(1:nel)

        dyy(1:nel) =dyy(1:nel)+bzy1(1:nel)*vz1(1:nel)+bzy2(1:nel)*vz2(1:nel)

     +        +bzy3(1:nel)*vz3(1:nel)+bzy4(1:nel)*vz4(1:nel)

     +        +bzy5(1:nel)*vz5(1:nel)+bzy6(1:nel)*vz6(1:nel)

     +        +bzy7(1:nel)*vz7(1:nel)+bzy8(1:nel)*vz8(1:nel)

      ENDIF !IF (ICP==1) THEN

C

      RETURN


      END

jacob_j33
subroutine jacob_j33(det, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, nel)
Definition jacob_j33.F:37

max
#define max(a, b)
Definition macros.h:21

s8zdefot3
subroutine s8zdefot3(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, icp, det0, jfac, nu, ideg, nel)
Definition s8zdefot3.F:71