volpfv.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/.
!||====================================================================
!||    volpfv             ../engine/source/airbag/volpfv.F
!||--- called by ------------------------------------------------------
!||    monvol0            ../engine/source/airbag/monvol0.F
!||--- uses       -----------------------------------------------------
!||    finter_mixed_mod   ../engine/source/tools/finter_mixed.F90
!||    python_funct_mod   ../common_source/modules/python_mod.F90
!||====================================================================
      SUBROUTINE volpfv(IVOLU,RVOLU  ,VOL,FSAV,NPC,
     2                  TF   ,PMAIN, WFEXT,PYTHON,NFUNCT)
         USE python_funct_mod
         USE finter_mixed_mod
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      "com06_c.inc"
#include      "com08_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPC(*),IVOLU(*), PMAIN
      my_real tf(*),rvolu(*),fsav(*)
      DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
      TYPE(python_), intent(inout) :: PYTHON
      INTEGER, intent(in) :: NFUNCT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ITFUN,ITYPFUN
      my_real vol,v0,area,veps,vinc,pold,pext,vold,pres,q,qold,dv,scale,xfun
C-----------------------------------------------
C
C CALCUL DE LA PRESSION
C
      itfun  =ivolu(7)
      itypfun=ivolu(19)
      pext   =rvolu(3)
      v0     =rvolu(4)
      vinc   =rvolu(5)
      pold   =rvolu(12)
      vold   =rvolu(16)
      veps   =rvolu(17)
      area   =rvolu(18)
      qold   =rvolu(23)
      scale  =rvolu(15)
      vol = vol + veps
      v0  =  v0
      SELECT CASE (itypfun)
      CASE(0)
         xfun = (v0-vinc)/(vol-vinc)
      CASE(1, 3)
         xfun = tt
      CASE(2)
         xfun = (vol-vinc) / (v0-vinc)
      CASE DEFAULT
         
      END SELECT
      pres = scale * finter_mixed(python,nfunct,itfun,xfun,npc,tf)
      IF (itypfun == 3) THEN
         pres = pres * (v0-vinc) / (vol-vinc)
      ENDIF
      q = zero
      IF (ispmd+1==pmain) THEN
        dv=vol-vold
        wfext=wfext+(half*(q+qold+pres+pold))*dv
        fsav(1)=zero
        fsav(2)=vol
        fsav(3)=pres
        fsav(4)=area
        fsav(5)=zero
        fsav(6)=zero
        fsav(7)=zero
        fsav(8)=zero
        fsav(9)=zero
        fsav(10)=zero
        fsav(11)=zero
        fsav(12)=zero
      ENDIF
C
      rvolu(16)=vol
      rvolu(12)=pres
      rvolu(23)=q
C
      RETURN
      END