sesame_tools.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/.
!||====================================================================
!||    tb2si1               ../starter/source/materials/eos/sesame_tools.F
!||--- called by ------------------------------------------------------
!||    hm_read_eos_sesame   ../starter/source/materials/eos/hm_read_eos_sesame.F
!||====================================================================
      SUBROUTINE tb2si1(DENS,PRES,ENER,NR,NT)
C-----------------------------------------------
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 NR,NT
      my_real
     .   dens(*),pres(nr,*),ener(nr,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real
     .   dtab2si,ptab2si,etab2si
C
C---------------------------------------------------------------
C Passage des unites table Sesame type 301-306 aux unites SI
C    1 Mg/m^3 --> 1000 Kg/m^3
C    1 GPa    --> 1E9 Kg/(m.s^2)
C    1 MJ/Kg  --> 1E6 m^2/s^2
C
C DTAB2SI : coef. de passage des masses volumiques tables aux unites SI
C PTAB2SI : coef. de passage des pressions tables aux unites SI
C ETAB2SI : coef. de passage des energies libres tables aux unites SI
C
      dtab2si = 1000.
      ptab2si = 1.e09
      etab2si = 1.e06
C
      DO i=1,nr
          dens(i)=dens(i)*dtab2si
      ENDDO
C
      DO  i=1,nr
        DO  j=1,nt
          pres(i,j)=pres(i,j)*ptab2si
          ener(i,j)=ener(i,j)*etab2si
        ENDDO
      ENDDO
C
      RETURN
      END
 
 
 
!||====================================================================
!||    tbusr1               ../starter/source/materials/eos/sesame_tools.F
!||--- called by ------------------------------------------------------
!||    hm_read_eos_sesame   ../starter/source/materials/eos/hm_read_eos_sesame.F
!||--- uses       -----------------------------------------------------
!||====================================================================
      SUBROUTINE tbusr1(R,P,E,NR,NT,UNITAB)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
      USE unitab_mod
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NR,NT
      TYPE(unit_type_), INTENT(IN) :: UNITAB
      my_real
     .   r(*),p(nr,*),e(nr,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real alscal,amscal,atscal,dscal,pscal,escal
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      amscal=unitab%FAC_M_WORK
      alscal=unitab%FAC_L_WORK
      atscal=unitab%FAC_T_WORK
C----------------------------------------------------------------
C Passage des unites table Sesame type 301-306 aux unites
C utilisateur passer obligatoirement avant l'appel de cette
C routine par la routine TB2SI1
C coefficients de passage des unites SI aux unites utilisateur
C M --> masse
C L --> longueur
C T --> temps
C masse volumique : M.L^-3
C pression        : M.l^-1.T^-2
C energie         : L^2.T^-2
C L'energy array est en fait equivalent a la chaleur specifique 
C de transition de phase. Cf Aide Memoire de Physique, B. YAVORSKI
C ET A.DETLAF, p 927.
C
      dscal = (alscal**3)/amscal
      pscal = alscal*(atscal**2)/amscal
      escal = (atscal**2)/(alscal**2)
C
      DO  i=1,nr
          r(i)=r(i)*dscal
      ENDDO
C
      DO  i=1,nr
        DO  j=1,nt
          p(i,j)=p(i,j)*pscal
          e(i,j)=e(i,j)*escal
        ENDDO
      ENDDO
C
      RETURN
      END