m37init.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/.
!||====================================================================
!||    m37init        ../starter/source/materials/mat/mat037/m37init.F
!||--- called by ------------------------------------------------------
!||    matini         ../starter/source/materials/mat_share/matini.F
!||--- uses       -----------------------------------------------------
!||====================================================================
      SUBROUTINE m37init(
     .               IPM     ,PM        ,
     .               NGL     ,NUVAR     ,UVAR    ,UPARAM    ,X       ,
     .               MAT     ,IPARG     ,IFORM   ,IX        ,NIX     ,
     .               ILOADP    ,FACLOAD ,GBUF      ,NEL)
C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutine is initializing cell based on law37
C In case of /LOAD/HYDRO is attached to it then
C   Pressure is computed so that P(t=0) = P0 + DP 
C     where DP=rho0.g.h
C     and h is signed distance from water/air free surface
C     which is defined by B and n in input
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_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      "com04_c.inc"
#include      "param_c.inc"
#include      "vect01_c.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER :: IPM(NPROPMI,NUMMAT),NGL(MVSIZ)      ,MAT(MVSIZ), IPARG(NPARG)   ,IFORM,NIX,IX(NIX,*)
      my_real :: PM(NPROPM,NUMMAT),X(3,NUMNOD), UVAR(LLT,NUVAR)
      INTEGER,INTENT(IN) :: ILOADP(SIZLOADP,*),NEL
      my_real,INTENT(IN) :: FACLOAD(LFACLOAD,*)
      TYPE(g_bufel_), INTENT(INOUT),TARGET :: GBUF  
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  :: I,J,NUVAR, GG1, GG2, GG3, ISFLUID,NL
      INTEGER  :: NPH,IFLG,IE,IADBUF,NUPARAM, M_UID,M_IID,IGRAV
      my_real  :: MFRAC, A, MU1P1, MU1P2, RHO1, RHO2, RHO10, RHO20, C1, GAM, PMIN, P0, MU2P1, RHO,P,A1,PSH,DEPTH,Z(3)
      my_real  :: UPARAM(*),XBAS,YBAS,ZBAS,NX,NY,NZ,GRAV, BZ(3),DOTPROD, DELTA_P, RHO0_LIQ,RHO0_GAS, RHO_LIQ,RHO_GAS
      my_real  :: pold,height
      my_real  :: n(3,8)
C-----------------------------------------------
 
      !---------------------------------!
      ! CELL INIT.   (UVAR ARRAY)       !
      !---------------------------------!
      !UVAR(I,1) : massic percentage of liquid * global density  (rho1*V1/V : it needs to give liquid mass multiplying by element volume in aleconve.F)
      !UVAR(I,2) : density of gas
      !UVAR(I,3) : density of liquid
      !UVAR(I,4) : volumetric fraction of liquid
      !UVAR(I,5) : volumetric fraction of gas
      DO i=lft,llt
        ie        =  i+nft
        c1        =  uparam(4)
        gam       =  uparam(5)
        pmin      =  uparam(8)
        p0        =  uparam(9)
        rho10     =  uparam(11)
        rho20     =  uparam(12)
        a1        =  uparam(10)         
        p         = -p0  
        rho       =  rho10 * a1 + (one-a1)*rho20     
        IF(gam*c1>=em30)THEN !if Liquid and gas correctly defined
          mu1p1     = one       !(P/C1+ONE)
          mu2p1     = one       !( ONE+P/P0)**(ONE/GAM)
          rho1      = rho10*mu1p1
          rho2      = rho20*mu2p1
          a         = (rho-rho2)/(rho1-rho2)
          uvar(i,1) = a*rho1
          uvar(i,2) = rho2
          uvar(i,3) = rho1
          uvar(i,4) = a
          IF(uvar(i,4)<em20)uvar(i,4)=zero
          uvar(i,5) = one-uvar(i,4)
          !SOUNDSP(I)  = A*SQRT(C1/RHO1) + (ONE-A) * SQRT(GAM*P0/RHO2)  -> check sigeps37.F (starter)
        ELSE !boundary element
          uvar(i,3) = rho
        ENDIF
      ENDDO
 
C-----------------------------------------------
      RETURN
      END SUBROUTINE m37init