sdlenmax.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/.
!||====================================================================
!||    sdlenmax   ../engine/source/elements/solid/solide/sdlenmax.F
!||--- called by ------------------------------------------------------
!||    s20forc3   ../engine/source/elements/solid/solide20/s20forc3.F
!||    s8eforc3   ../engine/source/elements/solid/solide8e/s8eforc3.F
!||    s8zforc3   ../engine/source/elements/solid/solide8z/s8zforc3.F
!||    sforc3     ../engine/source/elements/solid/solide/sforc3.F
!||    szforc3    ../engine/source/elements/solid/solidez/szforc3.F
!||====================================================================
      SUBROUTINE sdlenmax(LM_MAX, 
     1   X1,      X2,      X3,      X4,      
     2   X5,      X6,      X7,      X8,      
     3   Y1,      Y2,      Y3,      Y4,      
     4   Y5,      Y6,      Y7,      Y8,      
     5   Z1,      Z2,      Z3,      Z4,      
     6   Z5,      Z6,      Z7,      Z8,      
     7   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
      my_real, DIMENSION(MVSIZ) , INTENT(OUT)  :: LM_MAX 
      my_real, DIMENSION(MVSIZ) , INTENT(IN)  ::  
     .   X1, X2, X3, X4, X5, X6, X7, X8,
     .   Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,  
     .   Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, N
      my_real
     .   RX, RY, RZ, SX, SY, SZ, TX, TY, TZ,NORM
C=======================================================================
      DO I=1,nel
        rx = x1(i)+x2(i)+x5(i)+x6(i)-x3(i)-x4(i)-x7(i)-x8(i)
        ry = y1(i)+y2(i)+y5(i)+y6(i)-y3(i)-y4(i)-y7(i)-y8(i)
        rz = z1(i)+z2(i)+z5(i)+z6(i)-z3(i)-z4(i)-z7(i)-z8(i)
        sx = x5(i)+x6(i)+x7(i)+x8(i)-x1(i)-x2(i)-x3(i)-x4(i)
        sy = y5(i)+y6(i)+y7(i)+y8(i)-y1(i)-y2(i)-y3(i)-y4(i)
        sz = z5(i)+z6(i)+z7(i)+z8(i)-z1(i)-z2(i)-z3(i)-z4(i)
        tx = x3(i)+x2(i)+x7(i)+x6(i)-x1(i)-x4(i)-x5(i)-x8(i)
        ty = y3(i)+y2(i)+y7(i)+y6(i)-y1(i)-y4(i)-y5(i)-y8(i)
        tz = z3(i)+z2(i)+z7(i)+z6(i)-z1(i)-z4(i)-z5(i)-z8(i)
        norm = rx*rx + ry*ry + rz*rz
        lm_max(i) = norm
        norm = sx*sx + sy*sy + sz*sz
        lm_max(i) = max(lm_max(i),norm)
        norm = tx*tx + ty*ty + tz*tz
        lm_max(i) = max(lm_max(i),norm)
      ENDDO
      DO i=1,nel
        lm_max(i) = fourth*sqrt(lm_max(i))
      ENDDO 
C
      RETURN
      END