disp_vel_saved_cload.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/.
!||====================================================================
!||    disp_vel_saved_cload   ../engine/source/assembly/disp_vel_saved_cload.F
!||--- called by ------------------------------------------------------
!||    resol                  ../engine/source/engine/resol.F
!||====================================================================
      SUBROUTINE disp_vel_saved_cload(V      ,D      ,VR    ,DR    ,IB     ,
     2                                DPL0CLD,VEL0CLD,NIBCLD,NCONLD,IRODDL ,
     3                                NUMNOD )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      !USE CONSTANT_MOD
      ! ,ONLY : THIRD,FOURTH
!      USE PINCHTYPE_MOD
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, INTENT(IN) :: NIBCLD,NCONLD,IRODDL,NUMNOD
      INTEGER, INTENT(IN) :: IB(NIBCLD,NCONLD)
C     REAL
      my_real, INTENT(IN) :: 
     .   d(3,numnod) ,v(3,numnod) ,dr(3,numnod)  ,vr(3,numnod)
      my_real, INTENT(INOUT) :: 
     .   dpl0cld(6,nconld),vel0cld(6,nconld)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  NL,N1,N2,N3,N4,ISK,IDIR,IFUN
C-----------------------------------------------
C     IBCL(1:NIBCLD,1:NUMCLD)              IPRES = IBCL(1:NIBCLD,NUMCCLD+1,NUMCLD+NUMPRES)         
C     1: Node Number                       1st node number of the segment
C     2: NS = 10*Noskew+Idir               2nd node number of the segment
!! Save displacement and velocity in the previous cycle
      !DO NL=1,NCLOAD
      DO nl=1,nconld
        n1 = ib(1,nl)
        n2 = ib(2,nl)
        n3 = ib(3,nl)
        n4 = ib(4,nl)
        ifun = ib(9,nl)
        isk = ib(2,nl)/10
        idir  = ib(2,nl)-10*isk
        IF (n4 == -1) THEN
          IF (ifun > 1) THEN
            IF (idir<=3) THEN
              dpl0cld(idir,nl) = d(idir,n1)
              vel0cld(idir,nl) = v(idir,n1)
            ELSEIF (idir<=6 .AND. iroddl/=0) THEN
              dpl0cld(idir,nl) = dr(idir-3,n1)
              vel0cld(idir,nl) = vr(idir-3,n1)
            ENDIF
          ENDIF
        ENDIF  
      ENDDO
C
      RETURN
      END