flow_depla.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/.
!||====================================================================
!||    flow_depla   ../engine/source/assembly/flow_depla.F
!||--- called by ------------------------------------------------------
!||    resol        ../engine/source/engine/resol.F
!||====================================================================
      SUBROUTINE flow_depla(NALE,V    ,D    ,
     2                 NODFT,NODLT    )
C----6---------------------------------------------------------------7---------8
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
#include      "scr05_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODFT, NODLT, NALE(*)
      my_real d(3,*) ,v(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  :: N
      my_real :: vdt
      DOUBLE PRECISION VDT2
C-------------------------------------------------------------
C         CALCUL DES  DEPLACEMENTS TRANSLATIONS PUIS ROTATIONS
C-------------------------------------------------------------
 
      IF (iresp==1) THEN
#include      "vectorize.inc"
        DO n=nodft,nodlt
          IF(nale(n) == 3) THEN
             vdt2 = dt2*v(1,n)
             d(1,n)=d(1,n)+vdt2
cc             X(1,N)=X(1,N)+VDT2
cC             XDP(1,N)=XDP(1,N)+VDT2
C             X(1,N)=XDP(1,N)
             vdt2 = dt2*v(2,n)
             d(2,n)=d(2,n)+vdt2
cc             X(2,N)=X(2,N)+VDT2            
c             XDP(2,N)=XDP(2,N)+VDT2
c             X(2,N)=XDP(2,N)
             vdt2 = dt2*v(3,n)
             d(3,n)=d(3,n)+vdt2
cc             X(3,N)=X(3,N)+VDT2
C             XDP(3,N)=XDP(3,N)+VDT2
C             X(3,N)=XDP(3,N)
          ENDIF
        ENDDO
      ELSE
#include      "vectorize.inc"
        DO n=nodft,nodlt
          IF(nale(n) == 3) THEN
              vdt = dt2*v(1,n)
              d(1,n)=d(1,n)+vdt
cc              X(1,N)=X(1,N)+VDT
              vdt = dt2*v(2,n)
              d(2,n)=d(2,n)+vdt
cc              X(2,N)=X(2,N)+VDT
              vdt = dt2*v(3,n)
              d(3,n)=d(3,n)+vdt
cc              X(3,N)=X(3,N)+VDT
          ENDIF
        ENDDO
      ENDIF 
C
      RETURN
      END