m17law.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com08_c.inc"
#include "param_c.inc"
#include "scr14_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	m17law (pm, off, sig, eint, rho, rk, re, vorti, wxx, wyy, wzz, voln, mat, vis, d1, d2, d3, d4, d5, d6, nel, alogey, ssp, rho0, tmu, amu, amu2, psh, pc, espe, c1, c2, c3, c4, c5, c6, df, dpdm, jpor)

Function/Subroutine Documentation

m17law()

subroutine m17law	(		pm,
			off,
			sig,
			eint,
			rho,
			rk,
			re,
			vorti,
			wxx,
			wyy,
			wzz,
			voln,
		integer, dimension(*)	mat,
			vis,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
		integer	nel,
			alogey,
			ssp,
			rho0,
			tmu,
			amu,
			amu2,
			psh,
			pc,
			espe,
			c1,
			c2,
			c3,
			c4,
			c5,
			c6,
			df,
			dpdm,
		integer, intent(in)	jpor )

Definition at line 28 of file m17law.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.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      "com01_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
#include      "scr14_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JPOR
      INTEGER MAT(*),NEL
      my_real
     .   pm(npropm,*), off(*), sig(nel,6), eint(*), rho(*), rk(*), re(*),
     .   vorti(*),wxx(*),wyy(*),wzz(*),voln(mvsiz),vis(*),
     .   d1(*), d2(*), d3(*), d4(*), d5(*), d6(*), alogey(*), ssp(*),
     .   rho0(*), tmu(*), amu(*), amu2(*), psh(*), pc(*), espe(*),
     .   c1(*), c2(*), c3(*), c4(*), c5(*), c6(*), df(*), dpdm(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, MX
      my_real
     .   vis2(mvsiz), 
     .   dav(mvsiz),
     .   yp0, cmu, ax, e, a, xmu,
     .   xm, xk, xe, yplus, rk2t, fac,
     .   rho0_1, vis_1, pc_1, c1_1, c2_1,
     .   c3_1, c4_1, c5_1, c6_1, psh_1     
C-----------------------------------------------
C
      ! SAUVEGARDE VORTICITE SI ANIM
      IF((anim_e(10)==1 .OR. anim_se(10)==1).AND. dt1/=0.)THEN
C      IF(ANIM_E(10)==1.AND. DT1/=0.)THEN
        fac=four/dt1
        IF(n2d==0)THEN
          DO 5 i=1,nel
 5        vorti(i)=fac*sqrt(wxx(i)**2+wyy(i)**2+wzz(i)**2)
        ELSE
          DO 6 i=1,nel
 6        vorti(i)=fac*wzz(i)
        ENDIF
      ENDIF
C
      mx     =mat(1)
C
      rho0_1=pm( 1,mx)
      vis_1 =pm(24,mx)
      pc_1  =pm(37,mx)
      c1_1  =pm(31,mx)-pm(88,mx)
      c2_1  =pm(32,mx)
      c3_1  =pm(33,mx)
      c4_1  =pm(34,mx)
      c5_1  =pm(35,mx)
      c6_1  =pm(36,mx)
      psh_1 =pm(88,mx)
C
      DO 10 i=1,nel
      rho0(i)=rho0_1
      vis(i) =vis_1
      pc(i)  =pc_1
      c1(i)  =c1_1
      c2(i)  =c2_1
      c3(i)  =c3_1
      c4(i)  =c4_1
      c5(i)  =c5_1
      c6(i)  =c6_1
      psh(i) =psh_1
   10 CONTINUE
C
      mx     =mat(1)
      IF(jpor/=2)THEN
      DO 20 i=1,nel
      yp0=pm(51,mx)
      cmu=pm(81,mx)
      ax =pm(47,mx)
      e  =pm(48,mx)
      a  =pm(49,mx)
      xmu      =rho(i)*vis(i)
      xm       =rho(i)*voln(i)
      xk       =rk(i)/xm
      xe       = max(em15,re(i)/xm)
      yplus    =cmu*xk**2/(ax*xe*vis(i))
      yplus    = max(yplus,yp0)
      alogey(i)= a * log(e*yplus)
      vis(i)   =xmu*ax*yplus/alogey(i)
      tmu(i)   =vis(i)-xmu
   20 CONTINUE
      ELSE
        DO 21 i=1,nel
         vis(i)=zero
         tmu(i)=zero 
   21   CONTINUE
      ENDIF
C
      DO 30 i=1,nel
      vis2(i)=two*vis(i)
   30 dav(i) =-(d1(i)+d2(i)+d3(i))*third
C------------------------------
C     CONTRAINTES DEVIATORIQUES
C------------------------------
      DO 40 i=1,nel
      sig(i,1)=vis2(i)*(d1(i)+dav(i))
      sig(i,2)=vis2(i)*(d2(i)+dav(i))
      sig(i,3)=vis2(i)*(d3(i)+dav(i))
      sig(i,4)=vis(i) *d4(i)
      sig(i,5)=vis(i) *d5(i)
   40 sig(i,6)=vis(i) *d6(i)
C
      DO 50 i=1,nel
      df(i)  =rho0(i)/rho(i)
      amu(i) =one/df(i)-one
      amu2(i)= max(zero,amu(i))**2
   50 espe(i)=df(i)*eint(i)/voln(i)
C
      DO 60 i=1,nel
      rk2t=two*rk(i)/(three*voln(i))
      dpdm(i) = dpdm(i)
     .     +(c5(i)+c6(i)*amu(i))*df(i)*df(i)*rk2t + rk2t*df(i)
   60 CONTINUE
C
      DO 70 i=1,nel
   70 ssp(i)=sqrt(abs(dpdm(i))/rho0(i))
C
      RETURN