fvinjt8_1.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/.
!||====================================================================
!||    fvinjt8_1      ../engine/source/airbag/fvinjt8_1.F
!||--- called by ------------------------------------------------------
!||    fv_up_switch   ../engine/source/airbag/fv_up_switch.F
!||    fvbag1         ../engine/source/airbag/fvbag1.F
!||====================================================================
      SUBROUTINE fvinjt8_1(NJET   , IBAGJET , RBAGJET ,
     2                    IGEO,     GEO,     PM, 
     3                    IVOLU,   RVOLU, DMOUT, DHOUT)
C-----------------------------------------------
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NJET, IBAGJET(NIBJET,*), IGEO(NPROPGI,*), IVOLU(*)
C     REAL
      my_real
     .   dmout, dhout, rbagjet(nrbjet,*),
     .   geo(npropg,*), pm(npropm,*), rvolu(*)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IJ, IK, KK, I_INJ, I_TYPINJ, NGASES, I_GAS
      my_real
     .   r_igc1, rmwi, rmwg, rnm, rnmi, rnmg, mw, gmtot, gmi, dgmout, dgeout
C----------------------------------------------------------------------------
      r_igc1=pm(27,ivolu(66))
      rmwi  = rvolu(10)
C--------------------------
C     NOMBRE TOTAL DE MOLES
C--------------------------
      gmi = rvolu(11)
      rnm = gmi*rmwi
      DO ij=1,njet
        i_inj = iabs(ibagjet(13,ij))
        i_typinj = igeo(22,i_inj)
        ngases   = igeo(23,i_inj)
        DO ik=1,ngases
          IF (i_typinj==1) THEN
            i_gas = igeo(100+(ik-1)*3+1,i_inj)
          ELSE IF (i_typinj==2) THEN
            i_gas = igeo(100+(ik-1)*2+1,i_inj)
          END IF
          mw   = pm(20,i_gas)
          rmwg = r_igc1/mw
          kk=20+(ik-1)*4
          gmtot= rbagjet(kk+1,ij)
          rnmg = gmtot*rmwg
          rnm  = rnm+rnmg
        ENDDO
      ENDDO
      rvolu(14)=rnm
C-------------------------------
C     MASSE et ENTHALPY par GAZ
C-------------------------------
      rnmi=gmi*rmwi
      dgmout=dmout*rnmi/rnm
      dgeout=dhout*rnmi/rnm
C      DGEIN =DGMIN*RVOLU(63)
      rvolu(22)=dgeout
      rvolu(24)=dgmout
C      RVOLU(64)=DGMIN
C      RVOLU(65)=DGEIN
      rvolu(11)=rvolu(11)+dgmout*dt1
C
      DO ij=1,njet
         i_inj = iabs(ibagjet(13,ij))
         i_typinj = igeo(22,i_inj)
         ngases   = igeo(23,i_inj)
         rbagjet( 9,ij)=zero
         rbagjet(10,ij)=zero
         DO ik=1,ngases
          IF (i_typinj==1) THEN
            i_gas = igeo(100+(ik-1)*3+1,i_inj)
          ELSE IF (i_typinj==2) THEN
            i_gas = igeo(100+(ik-1)*2+1,i_inj)
          END IF
          mw   = pm(20,i_gas)
          rmwg = r_igc1/mw
          kk=20+(ik-1)*4
          gmtot= rbagjet(kk+1,ij)
          rnmg  =gmtot*rmwg
          dgmout=dmout*rnmg/rnm
          dgeout=dhout*rnmg/rnm
          rbagjet(kk+1,ij)=gmtot+dgmout*dt1
          rbagjet(kk+3,ij)=dgmout
          rbagjet(kk+4,ij)=dgeout
          rbagjet( 8,ij)=rbagjet( 8,ij)+dgmout*dt1
          rbagjet( 9,ij)=rbagjet( 9,ij)+dgmout
          rbagjet(10,ij)=rbagjet(10,ij)+dgeout
         ENDDO
      ENDDO
      RETURN
      END