Télécharger nano3d.eso

Retour à la liste

Numérotation des lignes :

nano3d
  1. C NANO3D SOURCE FD218221 24/02/07 21:15:21 11834
  2. subroutine nano3d(xmat,nmat,vnmat,nstype,var0,varf,nvari,vnvari,
  3. # ndtens,nbelas,teta1,teta2,tetaref,ppas,dt,kthermv,knanov)
  4.  
  5. c calcul des coeff hydrique 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,mfr
  17. logical ppas
  18. real*8 kthermv,knanov
  19.  
  20. c nrj activation pour l eau
  21. real*8 nrjw
  22. parameter (nrjw=17000.d0)
  23. real*8 wcsh0,pwcsh0,mdtt,tdtt,wdtt,pdtt,wcshf,pwcshf,pwcap
  24. integer nump,numw,numc
  25.  
  26.  
  27. call xmat3d(mdtt,xmat,nmat,vnmat,nstype,nbelas,3,1)
  28. call xmat3d(tdtt,xmat,nmat,vnmat,nstype,nbelas,3,2)
  29. call xmat3d(wdtt,xmat,nmat,vnmat,nstype,nbelas,3,3)
  30. call xmat3d(pdtt,xmat,nmat,vnmat,nstype,nbelas,3,4)
  31.  
  32. call xvar3d(pwcsh0,nump,var0,nvari,vnvari,nstype,ndtens,3,1)
  33. call xvar3d(wcsh0,numw,var0,nvari,vnvari,nstype,ndtens,3,2)
  34. call xvar3d(pwcap,numc,varf,nvari,vnvari,nstype,ndtens,1,3)
  35.  
  36. if(ppas) then
  37. c initialisation de la quantite d'eau dans les CSH a la valeur
  38. c de reference
  39. wcsh0=wdtt
  40. pwcsh0=0.d0
  41. end if
  42.  
  43. c prise en compte de la pression dans les csh pour la DTT
  44. c cf. F. Manzoni, T. Vidal, A. Sellier, X. Bourbon, and G. Camps,
  45. c “On the origins of transient thermal deformation of concrete,”
  46. c Cem. Concr. Compos., vol. 107, 2020,
  47. c https://doi.org/10.1016/j.cemconcomp.2019.103508
  48.  
  49. call wdtt3d(wcsh0,wcshf,tdtt,mdtt,wdtt,pdtt,kthermv,
  50. # teta1,teta2,pwcap,dt,knanov,pwcsh0,pwcshf,tetaref,ierr1)
  51. if(ierr1.eq.1) then
  52. print*,'ERREUR calcul de la DTT dans nano3d'
  53. ierr1=1
  54. stop
  55. end if
  56.  
  57. varf(nump)=pwcshf
  58. varf(numw)=wcshf
  59.  
  60. return
  61. end
  62.  
  63.  
  64.  

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