sfint43.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/.
!||====================================================================
!||    sfint43   ../engine/source/elements/solid/sconnect/sfint43.F
!||--- called by ------------------------------------------------------
!||    suser43   ../engine/source/elements/solid/sconnect/suser43.F
!||====================================================================
      SUBROUTINE sfint43(IPG   ,NPG   ,NEL   ,HH    ,AREAPG,SOFT  ,
     .     FX1   ,FX2   ,FX3   ,FX4   ,FX5   ,FX6   ,FX7   ,FX8   ,
     .     FY1   ,FY2   ,FY3   ,FY4   ,FY5   ,FY6   ,FY7   ,FY8   ,
     .     FZ1   ,FZ2   ,FZ3   ,FZ4   ,FZ5   ,FZ6   ,FZ7   ,FZ8   ,
     .     SIGNZZ,SIGNYZ,SIGNZX)
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-----------------------------------------------
      INTEGER IPG,NEL,NPG
      my_real, DIMENSION(NEL) :: AREAPG,SIGNZZ,SIGNYZ,SIGNZX,SOFT,
     .  FX1,FX2,FX3,FX4,FX5,FX6,FX7,FX8,FY1,FY2,FY3,FY4,FY5,FY6,FY7,FY8,         
     .  FZ1,FZ2,FZ3,FZ4,FZ5,FZ6,FZ7,FZ8
      my_real, DIMENSION(NPG,NPG) :: hh 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real FZZ,FYZ,FZX,RX,RY
      my_real, DIMENSION(NEL) :: DFX1,DFX2,DFX3,DFX4,DFX5,DFX6,DFX7,DFX8,
     .                           DFY1,DFY2,DFY3,DFY4,DFY5,DFY6,DFY7,DFY8,
     .                           dfz1,dfz2,dfz3,dfz4,dfz5,dfz6,dfz7,dfz8 
C=======================================================================
      DO i=1,nel
        fzz = signzz(i)*areapg(i)*soft(i)
        fyz = signyz(i)*areapg(i)*soft(i)
        fzx = signzx(i)*areapg(i)*soft(i)
C       inf
        dfx1(i) = fzx*hh(1,ipg)
        dfx2(i) = fzx*hh(2,ipg)
        dfx3(i) = fzx*hh(3,ipg)
        dfx4(i) = fzx*hh(4,ipg)
        dfy1(i) = fyz*hh(1,ipg)
        dfy2(i) = fyz*hh(2,ipg)
        dfy3(i) = fyz*hh(3,ipg)
        dfy4(i) = fyz*hh(4,ipg)
        dfz1(i) = fzz*hh(1,ipg)
        dfz2(i) = fzz*hh(2,ipg)
        dfz3(i) = fzz*hh(3,ipg)
        dfz4(i) = fzz*hh(4,ipg)
C       sup
        dfx5(i) =-fzx*hh(1,ipg)
        dfx6(i) =-fzx*hh(2,ipg)
        dfx7(i) =-fzx*hh(3,ipg)
        dfx8(i) =-fzx*hh(4,ipg)
        dfy5(i) =-fyz*hh(1,ipg)
        dfy6(i) =-fyz*hh(2,ipg)
        dfy7(i) =-fyz*hh(3,ipg)
        dfy8(i) =-fyz*hh(4,ipg)
        dfz5(i) =-fzz*hh(1,ipg)
        dfz6(i) =-fzz*hh(2,ipg)
        dfz7(i) =-fzz*hh(3,ipg)
        dfz8(i) =-fzz*hh(4,ipg)
      ENDDO
C-----------------------------
      DO i=1,nel
C       inf
        fx1(i) = fx1(i) + dfx1(i)
        fx2(i) = fx2(i) + dfx2(i)
        fx3(i) = fx3(i) + dfx3(i)
        fx4(i) = fx4(i) + dfx4(i)
        fy1(i) = fy1(i) + dfy1(i)
        fy2(i) = fy2(i) + dfy2(i)
        fy3(i) = fy3(i) + dfy3(i)
        fy4(i) = fy4(i) + dfy4(i)
        fz1(i) = fz1(i) + dfz1(i)
        fz2(i) = fz2(i) + dfz2(i)
        fz3(i) = fz3(i) + dfz3(i)
        fz4(i) = fz4(i) + dfz4(i)
C       sup             
        fx5(i) = fx5(i) + dfx5(i)
        fx6(i) = fx6(i) + dfx6(i)
        fx7(i) = fx7(i) + dfx7(i)
        fx8(i) = fx8(i) + dfx8(i)
        fy5(i) = fy5(i) + dfy5(i)
        fy6(i) = fy6(i) + dfy6(i)
        fy7(i) = fy7(i) + dfy7(i)
        fy8(i) = fy8(i) + dfy8(i)
        fz5(i) = fz5(i) + dfz5(i)
        fz6(i) = fz6(i) + dfz6(i)
        fz7(i) = fz7(i) + dfz7(i)
        fz8(i) = fz8(i) + dfz8(i)
      ENDDO
C-----------
      RETURN
      END