s8edefot12.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/.
!||====================================================================
!||    s8edefot12   ../engine/source/elements/solid/solide8e/s8edefot12.F
!||--- called by ------------------------------------------------------
!||    s8eforc3     ../engine/source/elements/solid/solide8e/s8eforc3.F
!||====================================================================
      SUBROUTINE s8edefot12(
     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   BXY1,    BXY2,    BXY3,    BXY4,
     8   BXY5,    BXY6,    BXY7,    BXY8,
     9   BYX1,    BYX2,    BYX3,    BYX4,
     A   BYX5,    BYX6,    BYX7,    BYX8,
     B   BXZ1,    BXZ2,    BXZ3,    BXZ4,
     C   BXZ5,    BXZ6,    BXZ7,    BXZ8,
     D   BZX1,    BZX2,    BZX3,    BZX4,
     E   BZX5,    BZX6,    BZX7,    BZX8,
     F   BYZ1,    BYZ2,    BYZ3,    BYZ4,
     G   BYZ5,    BYZ6,    BYZ7,    BYZ8,
     H   BZY1,    BZY2,    BZY3,    BZY4,
     I   BZY5,    BZY6,    BZY7,    BZY8,
     J   BXX1,    BXX2,    BXX3,    BXX4,
     K   BXX5,    BXX6,    BXX7,    BXX8,
     L   BYY1,    BYY2,    BYY3,    BYY4,
     M   BYY5,    BYY6,    BYY7,    BYY8,
     N   BZZ1,    BZZ2,    BZZ3,    BZZ4,
     O   BZZ5,    BZZ6,    BZZ7,    BZZ8,
     P   PXY1,    PXY2,    PXY3,    PXY4,
     Q   PXY5,    PXY6,    PXY7,    PXY8,
     R   PYX1,    PYX2,    PYX3,    PYX4,
     S   PYX5,    PYX6,    PYX7,    PYX8,
     T   PXZ1,    PXZ2,    PXZ3,    PXZ4,
     U   PXZ5,    PXZ6,    PXZ7,    PXZ8,
     V   PZX1,    PZX2,    PZX3,    PZX4,
     W   PZX5,    PZX6,    PZX7,    PZX8,
     X   PYZ1,    PYZ2,    PYZ3,    PYZ4,
     Y   PYZ5,    PYZ6,    PYZ7,    PYZ8,
     Z   PZY1,    PZY2,    PZY3,    PZY4,
     1   PZY5,    PZY6,    PZY7,    PZY8,
     2   VX1,     VX2,     VX3,     VX4,
     3   VX5,     VX6,     VX7,     VX8,
     4   VY1,     VY2,     VY3,     VY4,
     5   VY5,     VY6,     VY7,     VY8,
     6   VZ1,     VZ2,     VZ3,     VZ4,
     7   VZ5,     VZ6,     VZ7,     VZ8,
     8   R11,     R12,     R13,     R21,
     9   R22,     R23,     R31,     R32,
     A   R33,     DXX,     DXY,     DXZ,
     B   DYX,     DYY,     DYZ,     DZX,
     C   DZY,     DZZ,     ICP,     DET0,
     D   JFAC,    NU,      IDEGE,   OFFG,
     E   ISEL_V,  I_SH,    NEL,     ISMSTR,
     F   JCVT,    MTN ,    DIV0 )
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, INTENT(IN) :: ISMSTR
      INTEGER, INTENT(IN) :: JCVT,MTN
      INTEGER ISEL_V,I_SH,ICP,IDEGE(*)
C     REAL
      my_real
     .   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(MVSIZ),PXY2(MVSIZ),PXY3(MVSIZ),PXY4(MVSIZ),
     .   PXY5(MVSIZ),PXY6(MVSIZ),PXY7(MVSIZ),PXY8(MVSIZ),
     .   PYX1(MVSIZ),PYX2(MVSIZ),PYX3(MVSIZ),PYX4(MVSIZ),
     .   PYX5(MVSIZ),PYX6(MVSIZ),PYX7(MVSIZ),PYX8(MVSIZ),
     .   PXZ1(MVSIZ),PXZ2(MVSIZ),PXZ3(MVSIZ),PXZ4(MVSIZ),
     .   PXZ5(MVSIZ),PXZ6(MVSIZ),PXZ7(MVSIZ),PXZ8(MVSIZ),
     .   PZX1(MVSIZ),PZX2(MVSIZ),PZX3(MVSIZ),PZX4(MVSIZ),
     .   PZX5(MVSIZ),PZX6(MVSIZ),PZX7(MVSIZ),PZX8(MVSIZ),
     .   PYZ1(MVSIZ),PYZ2(MVSIZ),PYZ3(MVSIZ),PYZ4(MVSIZ),
     .   PYZ5(MVSIZ),PYZ6(MVSIZ),PYZ7(MVSIZ),PYZ8(MVSIZ),
     .   PZY1(MVSIZ),PZY2(MVSIZ),PZY3(MVSIZ),PZY4(MVSIZ),
     .   PZY5(MVSIZ),PZY6(MVSIZ),PZY7(MVSIZ),PZY8(MVSIZ),
     .   VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),
     .   VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),
     .   VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),
     .   R11(*), R12(*), R13(*),
     .   R21(*), R22(*), R23(*),
     .   R31(*), R32(*), R33(*),
     .   DXX(*), DXY(*), DXZ(*),
     .   DYX(*), DYY(*), DYZ(*),
     .   DZX(*), DZY(*), DZZ(*),
     .   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(*),
     .   BXX1(*),BXX2(*),BXX3(*),BXX4(*),
     .   BXX5(*),BXX6(*),BXX7(*),BXX8(*),
     .   BYY1(*),BYY2(*),BYY3(*),BYY4(*),
     .   BYY5(*),BYY6(*),BYY7(*),BYY8(*),
     .   BZZ1(*),BZZ2(*),BZZ3(*),BZZ4(*),
     .   bzz5(*),bzz6(*),bzz7(*),bzz8(*),offg(*),det0(*),jfac(*),nu(*)
      my_real, DIMENSION(MVSIZ), INTENT(IN) :: div0
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   VLX1(MVSIZ), VLX2(MVSIZ), VLX3(MVSIZ), VLX4(MVSIZ),
     .   VLX5(MVSIZ), VLX6(MVSIZ), VLX7(MVSIZ), VLX8(MVSIZ),
     .   VLY1(MVSIZ), VLY2(MVSIZ), VLY3(MVSIZ), VLY4(MVSIZ),
     .   VLY5(MVSIZ), VLY6(MVSIZ), VLY7(MVSIZ), VLY8(MVSIZ),
     .   VLZ1(MVSIZ), VLZ2(MVSIZ), VLZ3(MVSIZ), VLZ4(MVSIZ),
     .   VLZ5(MVSIZ), VLZ6(MVSIZ), VLZ7(MVSIZ), VLZ8(MVSIZ)
      my_real
     .   JAC_59_68, JAC_67_49, JAC_48_57,
     .   AJ11, AJ22, AJ33,DET(MVSIZ),FAC,BASE,DDIV
C-----------------------------------------------
      IF (JCVT > 0 ) THEN
        DO I=1,nel
         IF(offg(i)<=one) cycle
         vlx1(i)=r11(i)*vx1(i)+r21(i)*vy1(i)+r31(i)*vz1(i)
         vly1(i)=r12(i)*vx1(i)+r22(i)*vy1(i)+r32(i)*vz1(i)
         vlz1(i)=r13(i)*vx1(i)+r23(i)*vy1(i)+r33(i)*vz1(i)
         vlx2(i)=r11(i)*vx2(i)+r21(i)*vy2(i)+r31(i)*vz2(i)
         vly2(i)=r12(i)*vx2(i)+r22(i)*vy2(i)+r32(i)*vz2(i)
         vlz2(i)=r13(i)*vx2(i)+r23(i)*vy2(i)+r33(i)*vz2(i)
         vlx3(i)=r11(i)*vx3(i)+r21(i)*vy3(i)+r31(i)*vz3(i)
         vly3(i)=r12(i)*vx3(i)+r22(i)*vy3(i)+r32(i)*vz3(i)
         vlz3(i)=r13(i)*vx3(i)+r23(i)*vy3(i)+r33(i)*vz3(i)
         vlx4(i)=r11(i)*vx4(i)+r21(i)*vy4(i)+r31(i)*vz4(i)
         vly4(i)=r12(i)*vx4(i)+r22(i)*vy4(i)+r32(i)*vz4(i)
         vlz4(i)=r13(i)*vx4(i)+r23(i)*vy4(i)+r33(i)*vz4(i)
         vlx5(i)=r11(i)*vx5(i)+r21(i)*vy5(i)+r31(i)*vz5(i)
         vly5(i)=r12(i)*vx5(i)+r22(i)*vy5(i)+r32(i)*vz5(i)
         vlz5(i)=r13(i)*vx5(i)+r23(i)*vy5(i)+r33(i)*vz5(i)
         vlx6(i)=r11(i)*vx6(i)+r21(i)*vy6(i)+r31(i)*vz6(i)
         vly6(i)=r12(i)*vx6(i)+r22(i)*vy6(i)+r32(i)*vz6(i)
         vlz6(i)=r13(i)*vx6(i)+r23(i)*vy6(i)+r33(i)*vz6(i)
         vlx7(i)=r11(i)*vx7(i)+r21(i)*vy7(i)+r31(i)*vz7(i)
         vly7(i)=r12(i)*vx7(i)+r22(i)*vy7(i)+r32(i)*vz7(i)
         vlz7(i)=r13(i)*vx7(i)+r23(i)*vy7(i)+r33(i)*vz7(i)
         vlx8(i)=r11(i)*vx8(i)+r21(i)*vy8(i)+r31(i)*vz8(i)
         vly8(i)=r12(i)*vx8(i)+r22(i)*vy8(i)+r32(i)*vz8(i)
         vlz8(i)=r13(i)*vx8(i)+r23(i)*vy8(i)+r33(i)*vz8(i)
        ENDDO
      ELSE
         vlx1(1:nel)=vx1(1:nel)
         vly1(1:nel)=vy1(1:nel)
         vlz1(1:nel)=vz1(1:nel)
         vlx2(1:nel)=vx2(1:nel)
         vly2(1:nel)=vy2(1:nel)
         vlz2(1:nel)=vz2(1:nel)
         vlx3(1:nel)=vx3(1:nel)
         vly3(1:nel)=vy3(1:nel)
         vlz3(1:nel)=vz3(1:nel)
         vlx4(1:nel)=vx4(1:nel)
         vly4(1:nel)=vy4(1:nel)
         vlz4(1:nel)=vz4(1:nel)
         vlx5(1:nel)=vx5(1:nel)
         vly5(1:nel)=vy5(1:nel)
         vlz5(1:nel)=vz5(1:nel)
         vlx6(1:nel)=vx6(1:nel)
         vly6(1:nel)=vy6(1:nel)
         vlz6(1:nel)=vz6(1:nel)
         vlx7(1:nel)=vx7(1:nel)
         vly7(1:nel)=vy7(1:nel)
         vlz7(1:nel)=vz7(1:nel)
         vlx8(1:nel)=vx8(1:nel)
         vly8(1:nel)=vy8(1:nel)
         vlz8(1:nel)=vz8(1:nel)
      END IF !(JCVT > 0
C   
      IF (i_sh>0) THEN
       DO i=1,nel
        IF(offg(i)<=one) cycle
        dxy(i) =pxy1(i)*vlx1(i)+pxy2(i)*vlx2(i)
     +         +pxy3(i)*vlx3(i)+pxy4(i)*vlx4(i)
     +         +pxy5(i)*vlx5(i)+pxy6(i)*vlx6(i)
     +         +pxy7(i)*vlx7(i)+pxy8(i)*vlx8(i)
        dxz(i) =pxz1(i)*vlx1(i)+pxz2(i)*vlx2(i)
     +         +pxz3(i)*vlx3(i)+pxz4(i)*vlx4(i)
     +         +pxz5(i)*vlx5(i)+pxz6(i)*vlx6(i)
     +         +pxz7(i)*vlx7(i)+pxz8(i)*vlx8(i)
        dyx(i) =pyx1(i)*vly1(i)+pyx2(i)*vly2(i)
     +         +pyx3(i)*vly3(i)+pyx4(i)*vly4(i)
     +         +pyx5(i)*vly5(i)+pyx6(i)*vly6(i)
     +         +pyx7(i)*vly7(i)+pyx8(i)*vly8(i)
        dyz(i) =pyz1(i)*vly1(i)+pyz2(i)*vly2(i)
     +         +pyz3(i)*vly3(i)+pyz4(i)*vly4(i)
     +         +pyz5(i)*vly5(i)+pyz6(i)*vly6(i)
     +         +pyz7(i)*vly7(i)+pyz8(i)*vly8(i)
        dzx(i) =pzx1(i)*vlz1(i)+pzx2(i)*vlz2(i)
     +         +pzx3(i)*vlz3(i)+pzx4(i)*vlz4(i)
     +         +pzx5(i)*vlz5(i)+pzx6(i)*vlz6(i)
     +         +pzx7(i)*vlz7(i)+pzx8(i)*vlz8(i)
        dzy(i) =pzy1(i)*vlz1(i)+pzy2(i)*vlz2(i)
     +         +pzy3(i)*vlz3(i)+pzy4(i)*vlz4(i)
     +         +pzy5(i)*vlz5(i)+pzy6(i)*vlz6(i)
     +         +pzy7(i)*vlz7(i)+pzy8(i)*vlz8(i)
       ENDDO
      ELSE
       DO i=1,nel
        IF(offg(i)<=one) cycle
        dxy(i) =py1(i)*vlx1(i)+py2(i)*vlx2(i)+py3(i)*vlx3(i)+py4(i)*vlx4(i)
     +         +py5(i)*vlx5(i)+py6(i)*vlx6(i)+py7(i)*vlx7(i)+py8(i)*vlx8(i)
        dxz(i) =pz1(i)*vlx1(i)+pz2(i)*vlx2(i)+pz3(i)*vlx3(i)+pz4(i)*vlx4(i)
     +         +pz5(i)*vlx5(i)+pz6(i)*vlx6(i)+pz7(i)*vlx7(i)+pz8(i)*vlx8(i)
        dyx(i) =px1(i)*vly1(i)+px2(i)*vly2(i)+px3(i)*vly3(i)+px4(i)*vly4(i)
     +         +px5(i)*vly5(i)+px6(i)*vly6(i)+px7(i)*vly7(i)+px8(i)*vly8(i)
        dyz(i) =pz1(i)*vly1(i)+pz2(i)*vly2(i)+pz3(i)*vly3(i)+pz4(i)*vly4(i)
     +         +pz5(i)*vly5(i)+pz6(i)*vly6(i)+pz7(i)*vly7(i)+pz8(i)*vly8(i)
        dzx(i) =px1(i)*vlz1(i)+px2(i)*vlz2(i)+px3(i)*vlz3(i)+px4(i)*vlz4(i)
     +         +px5(i)*vlz5(i)+px6(i)*vlz6(i)+px7(i)*vlz7(i)+px8(i)*vlz8(i)
        dzy(i) =py1(i)*vlz1(i)+py2(i)*vlz2(i)+py3(i)*vlz3(i)+py4(i)*vlz4(i)
     +         +py5(i)*vlz5(i)+py6(i)*vlz6(i)+py7(i)*vlz7(i)+py8(i)*vlz8(i)
       ENDDO
      END IF !(ISMSTR==11) THEN
C     
       DO i=1,nel
        IF(offg(i)<=one) cycle
        dxx(i) =px1(i)*vlx1(i)+px2(i)*vlx2(i)+px3(i)*vlx3(i)+px4(i)*vlx4(i)
     .         +px5(i)*vlx5(i)+px6(i)*vlx6(i)+px7(i)*vlx7(i)+px8(i)*vlx8(i)
        dyy(i) =py1(i)*vly1(i)+py2(i)*vly2(i)+py3(i)*vly3(i)+py4(i)*vly4(i)
     .         +py5(i)*vly5(i)+py6(i)*vly6(i)+py7(i)*vly7(i)+py8(i)*vly8(i)
        dzz(i) =pz1(i)*vlz1(i)+pz2(i)*vlz2(i)+pz3(i)*vlz3(i)+pz4(i)*vlz4(i)
     .         +pz5(i)*vlz5(i)+pz6(i)*vlz6(i)+pz7(i)*vlz7(i)+pz8(i)*vlz8(i)
       ENDDO
       IF (isel_v>0) THEN
#include   "nofusion.inc"
         DO i=1,nel
          IF(offg(i)<=one) cycle
          dyy(i) =dyy(i)+bxy1(i)*vlx1(i)+bxy2(i)*vlx2(i)
     +                  +bxy3(i)*vlx3(i)+bxy4(i)*vlx4(i)
     +                  +bxy5(i)*vlx5(i)+bxy6(i)*vlx6(i)
     +                  +bxy7(i)*vlx7(i)+bxy8(i)*vlx8(i)
          dzz(i) =dzz(i)+bxz1(i)*vlx1(i)+bxz2(i)*vlx2(i)
     +                  +bxz3(i)*vlx3(i)+bxz4(i)*vlx4(i)
     +                  +bxz5(i)*vlx5(i)+bxz6(i)*vlx6(i)
     +                  +bxz7(i)*vlx7(i)+bxz8(i)*vlx8(i)
          dxx(i) =dxx(i)+byx1(i)*vly1(i)+byx2(i)*vly2(i)
     +                  +byx3(i)*vly3(i)+byx4(i)*vly4(i)
     +                  +byx5(i)*vly5(i)+byx6(i)*vly6(i)
     +                  +byx7(i)*vly7(i)+byx8(i)*vly8(i)
          dzz(i) =dzz(i)+byz1(i)*vly1(i)+byz2(i)*vly2(i)
     +                  +byz3(i)*vly3(i)+byz4(i)*vly4(i)
     +                  +byz5(i)*vly5(i)+byz6(i)*vly6(i)
     +                  +byz7(i)*vly7(i)+byz8(i)*vly8(i)
          dxx(i) =dxx(i)+bzx1(i)*vlz1(i)+bzx2(i)*vlz2(i)
     +                  +bzx3(i)*vlz3(i)+bzx4(i)*vlz4(i)
     +                  +bzx5(i)*vlz5(i)+bzx6(i)*vlz6(i)
     +                  +bzx7(i)*vlz7(i)+bzx8(i)*vlz8(i)
          dyy(i) =dyy(i)+bzy1(i)*vlz1(i)+bzy2(i)*vlz2(i)
     +                  +bzy3(i)*vlz3(i)+bzy4(i)*vlz4(i)
     +                  +bzy5(i)*vlz5(i)+bzy6(i)*vlz6(i)
     +                  +bzy7(i)*vlz7(i)+bzy8(i)*vlz8(i)
         ENDDO
      END IF !(ISEL_V>0) THEN
      IF (icp==1) THEN
       IF (mtn==1) THEN
        DO i=1,nel
         IF(idege(i)>10.OR.offg(i)<=one) cycle
         fac=two_third*nu(i)
         base = dxx(i)+dyy(i)+dzz(i)
         ddiv = fac*(div0(i)-base)
         dxx(i) = dxx(i)+ddiv
         dyy(i) = dyy(i)+ddiv
         dzz(i) = dzz(i)+ddiv
         jfac(i)= one   ! used in fint3 later
        ENDDO
       ELSE
        DO i=1,nel
         IF(offg(i)<=one) cycle
         aj11=dxx(i)+one
         aj22=dyy(i)+one
         aj33=dzz(i)+one
         jac_59_68=aj22*aj33-dyz(i)*dzy(i)
         jac_67_49=dyz(i)*dzx(i)-dyx(i)*aj33
         jac_48_57=dyx(i)*dzy(i)-aj22*dzx(i)
C       
         det(i)=aj11*jac_59_68+dxy(i)*jac_67_49+dxz(i)*jac_48_57
C-----F'=FAC*F=MF'+[1]-->MF'=FAC*MF+(FAC-ONE)*[1]
         fac=two_third*nu(i)
         base = det0(i)/max(em20,det(i))
         jfac(i)=exp(fac*log(max(em20,base)))
        ENDDO
        DO i=1,nel
         IF(idege(i)>10.OR.offg(i)<=one) cycle
         dxx(i) = jfac(i)*dxx(i)+jfac(i)-one
         dyy(i) = jfac(i)*dyy(i)+jfac(i)-one
         dzz(i) = jfac(i)*dzz(i)+jfac(i)-one
         dxy(i) = jfac(i)*dxy(i)
         dyx(i) = jfac(i)*dyx(i)
         dzx(i) = jfac(i)*dzx(i)
         dxz(i) = jfac(i)*dxz(i)
         dyz(i) = jfac(i)*dyz(i)
         dzy(i) = jfac(i)*dzy(i)
        ENDDO
       END IF !(MTN==1) THEN
      END IF !(ICP==1) THEN
C
      RETURN
      END