Télécharger la plaquette Cast3M


fr en

Télécharger ther3d.eso

Retour à la liste

Numérotation des lignes : 

ther3d

  1. C THER3D    SOURCE    FD218221  24/02/07    21:15:28     11834          
  2.        subroutine ther3d(xmat,nmat,vnmat,nstype,var0,varf,nvari,
  3.      # vnvari,ndtens,nbelas,teta1,teta2,tetaref,dtherm,kthermv,kthermp)
  4.  
  5. c      calcul des coeff thermiques pour mc3d
  6.  
  7. c      declarations externes
  8.        implicit integer(i-n)
  9.        implicit real*8(a-h,o-z)
  10.  
  11.        integer nmat,nstype,ndtens,nvari,nbelas
  12.        integer vnmat(nstype)
  13.        integer vnvari(nstype,2)
  14.        real*8 xmat(nmat)
  15.        real*8 var0(nvari),varf(nvari)
  16.        real*8 teta1,teta2,tetaref
  17.        real*8 dtherm,kthermv,kthermp
  18.  
  19. c      nrj activation pour l eau
  20.        real*8 nrjw 
  21.        parameter (nrjw=17000.d0)
  22.        real*8 teta,nrjf,tetas,DT80,dtherm0
  23.        integer num
  24.        real*8 vw,poro
  25.        real*8 CTHP,CTHV
  26.  
  27. c      on prend la temperature milieu de pas
  28.        teta=0.5d0*(teta1+teta2)
  29. c      recuperation des donnees
  30.  
  31. c      sous type 2=Mawxell
  32. c      sous num  3=nrjf
  33.        call xmat3d(nrjf,xmat,nmat,vnmat,nstype,nbelas,2,3)
  34.        call xmat3d(tetas,xmat,nmat,vnmat,nstype,nbelas,2,4)
  35.        call xmat3d(dt80,xmat,nmat,vnmat,nstype,nbelas,2,5)
  36.        call xmat3d(vw,xmat,nmat,vnmat,nstype,nbelas,1,3)
  37.        call xmat3d(poro,xmat,nmat,vnmat,nstype,nbelas,1,2)
  38. c      recup valeur et numero d une vari de type scalaire
  39.        call xvar3d(dtherm0,num,var0,nvari,vnvari,nstype,ndtens,1,6)
  40.  
  41.        call thrm3d(teta,nrjf,tetas,tetaref,DT80,
  42.      # dtherm0,dtherm,CTHp,CTHv,poro,vw,nrjw)
  43.  
  44. c      stockage de la vari calculee
  45.        varf(num)=dtherm
  46.        kthermv=CTHv
  47.        kthermp=CTHp
  48.  
  49.        return
  50.        end
  51.  
  52.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales