idege8.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/.
!||====================================================================
!||    idege8       ../engine/source/elements/solid/solide/idege8.F
!||--- called by ------------------------------------------------------
!||    sdlen_dege   ../engine/source/elements/solid/solide/sdlen_dege.F
!||    sdlen_sms    ../engine/source/elements/solid/solidez/sdlen_sms.F
!||====================================================================
      SUBROUTINE idege8(X1,X2,X3,X4,Y1,Y2,Y3,Y4,Z1,Z2,Z3,Z4,AMAX,FAC,IT)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   x1, x2, x3, x4, y1, y2, y3, y4,
     .   z1, z2, z3, z4,amax,fac
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IDE,IT
      my_real
     .   x12,y12,z12,x23,y23,z23,x34,y34,z34,x41,y41,z41,a,
     .   rx,ry,rz,sx,sy,sz,nx,ny,nz
C-----------------------------------------------
       rx = x2+x3-x1-x4
       ry = y2+y3-y1-y4
       rz = z2+z3-z1-z4
       sx = x3+x4-x1-x2
       sy = y3+y4-y1-y2
       sz = z3+z4-z1-z2
       nx = ry * sz - rz * sy 
       ny = rz * sx - rx * sz 
       nz = rx * sy - ry * sx 
C--------A: (2*AREA)^2-----       
       a = nx*nx+ny*ny+nz*nz
       ide = 0
       x12 = x2 - x1
       y12 = y2 - y1
       z12 = z2 - z1
       IF (x12==zero.AND.y12==zero.AND.z12==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x23 = x2 - x3
       y23 = y2 - y3
       z23 = z2 - z3
       IF (x23==zero.AND.y23==zero.AND.z23==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x34 = x4 - x3
       y34 = y4 - y3
       z34 = z4 - z3
       IF (x34==zero.AND.y34==zero.AND.z34==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x41 = x4 - x1
       y41 = y4 - y1
       z41 = z4 - z1
       IF (x41==zero.AND.y41==zero.AND.z41==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
   10   CONTINUE
C     
       IF (a > amax )THEN
              it=ide
          amax = a
             END IF
      RETURN
      END