frho3.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/.
!||====================================================================
!||    frho3     ../engine/source/elements/solid/solide/frho3.f
!||--- called by ------------------------------------------------------
!||    sforc3    ../engine/source/elements/solid/solide/sforc3.F
!||--- uses       -----------------------------------------------------
!||    ale_mod   ../common_source/modules/ale/ale_mod.F
!||====================================================================
      SUBROUTINE frho3(UVAR ,RHO0,DXX, DYY, DZZ,
     .                 FLUX,FLU1,VOLN,NGL,MAT,OFF,NEL)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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-----------------------------------------------
#include      "units_c.inc"
#include      "scr06_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NGL(*),MAT(*),NEL
      my_real uvar(nel,*),flux(mvsiz,6), flu1(*), voln(*), dxx(*),dyy(*),dzz(*),off(*),rho0(*)
      TARGET :: uvar
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real dvv, e0(nel),vavg,rv,rvp
      INTEGER I, J,COUNT,LIST(NEL),II
      my_real, DIMENSION(:), POINTER :: volo,rhon,eint,dvol
C-----------------------------------------------
      rhon => uvar(1:nel,1)   
      eint => uvar(1:nel,2)   
      volo => uvar(1:nel,3)   
      dvol => uvar(1:nel,4)   
C
C  incompressible flow avoir plus en cas de besoin
C      INCOMP = 1 not tested
      IF(ale%GLOBAL%INCOMP == 1)THEN
         e0 = zero ! 
         DO i=1,nel
          dvv=dt1*(dxx(i)+dyy(i)+dzz(i))
          rhon(i)=rhon(i)-dvv*rho0(i)
          vavg=half*(voln(i)+volo(i))
          dvol(i)=vavg*dvv
          eint(i)=eint(i)*volo(i)-e0(i)*dvv*vavg
         ENDDO
C compresssible flow
      ELSE
         DO i=1,nel
           rhon(i) = rhon(i)/voln(i)
          
           dvol(i) = voln(i) 
     .             - volo(i)+0.5*dt1*(flu1(i)+flux(i,1)+flux(i,2)
     .             +  flux(i,3)+flux(i,4)+flux(i,5)+flux(i,6) )
            volo(i)=voln(i)
           rhon(i) = max(rhon(i) , em30)
       
        ENDDO 
      ENDIF
C
        count=0
        DO  i=1,nel
          IF(rhon(i) <= zero .AND. off(i)/= zero)THEN
             count = count + 1
              list(count)=i
           ENDIF  
        ENDDO
        DO ii=1,count
         i=list(ii)
#include "lockon.inc"
        WRITE(iout,1000) ngl(i),rhon(i)
        WRITE(istdo,1000) ngl(i),rhon(i)
#include "lockoff.inc"
 1000   FORMAT('-STOP GAZ ',i8,' HAS NEGATIVE RHO',g12.5)
       ENDDO
C
      RETURN
      END SUBROUTINE frho3