Télécharger la plaquette Cast3M


fr en

Télécharger thrm3d.eso

Retour à la liste

Numérotation des lignes : 

thrm3d

  1. C THRM3D    SOURCE    PV090527  23/01/27    21:16:05     11574          
  2.       subroutine thrm3d(teta1,nrjm,tetas,tetar,DT80,
  3.      # dth0,DTH,CTHP,CTHV,poro,vw,nrjw)
  4.  
  5. c     influence de la temperature sur les parametres de fluage 
  6. c     et endommagement thermique   
  7.  
  8. c     declarations externes
  9.       IMPLICIT INTEGER(I-N)
  10.       IMPLICIT REAL*8(A-H,O-Z)
  11.  
  12. c     variables externes
  13.       real*8 theta,nrjm,taum00,taum,tauk00,tauk,
  14.      #epsk00,epsk,tetas,tetar,Srw,poro,vw,nrjw
  15.  
  16. c     variables locales      
  17.       real*8 easurrm,easurrw,unsurtr,unsurts
  18.  
  19. c***********************************************************************
  20. c     reglage des activation thermiques pour le fluage
  21. c     +++ les teta sont en degres Celsius +++
  22. c***********************************************************************
  23.  
  24. c     calcul du terme d activation d Arrhenius pour le potentiel      
  25.       easurrm=nrjm/8.314D0
  26. c     cas de l'activation des viscostés
  27.       easurrw=nrjw/8.314D0
  28. c     la temperature de reference est tetar      
  29.       unsurtr=1.d0/(tetar+273.15d0)
  30. c     unsurts la temperature de seuil pour la modif du potentiel 
  31.       unsurts=1.d0/(tetas+273.15d0)
  32. c     calcul des coeff d activation thermique
  33.       unsurt=(1.D0/(teta1+273.15D0))   
  34. c     cas de l eau      
  35.       CTHV=exp(-easurrw*(unsurt-unsurtr)) 
  36. c     cas de l endommagement thermique 
  37.       if((dt80.gt.0.).and.(teta1.gt.tetas)) then      
  38.         xxx1=unsurts-unsurt      
  39.         xxx2=0.5D0*(xxx1+abs(xxx1))     
  40.         CTHP=exp(easurrm*xxx2)
  41.       else
  42.         CTHP=1.d0
  43.       end if
  44. c     saturation
  45.       Srw=max(min(vw/poro,1.d0),0.d0)    
  46.  
  47.  
  48. ************************************************************************
  49. *    endommagement thermique
  50. ************************************************************************
  51.       if((tetas.lt.80.d0).and.(DT80.gt.0.).and.(teta1.gt.tetas)) then
  52. c         prise en compte de l attenuation du potentiel en cas
  53. c         de sechage      
  54.           unsurt80=2.8316d-3
  55.           xxx180=unsurts-unsurt80
  56.           xxx280=0.5D0*(xxx180+abs(xxx180))
  57.           Cth80=exp(easurrm*xxx280)
  58.           Ath=1.d0/(Cth80-1.d0)*(srw*DT80/(1.d0-Srw*DT80))
  59.           Dth=1.d0-1.d0/(1.d0+Ath*(CTHP-1.d0))
  60.       else
  61.           Dth=0.d0
  62.       end if
  63.       DTH=max(dth0,dth)
  64. c      print*,'DTH ds therma t3d',dth
  65. c      print*,teta1,nrjm,tetas,tetar,DT80,
  66. c     # dth0,DTH,CTHP,CTHV
  67. ************************************************************************      
  68.  
  69.       return
  70.       end
  71.  
  72.  
  73.  
  74.  

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