Télécharger dilth.eso

Retour à la liste

Numérotation des lignes :

  1. C DILTH SOURCE STRU 08/12/22 21:15:09 6246
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE DILTH (TT,TTP, PDILT, VEPSTH)
  5. C --------------------------------------------------------------------------
  6. C Calcul des vitesses de dilatation thermique VEPSTH(IJ)
  7. C en fonction de :
  8. C - la temperature TT,
  9. C - sa derivee par rapport au temps TTP,
  10. C - et des parametres transmis dans PDILT,
  11. C pour un materiau orthotrope, dans la base principale d'orthotropie.
  12. C --------------------------------------------------------------------------
  13. IMPLICIT INTEGER (I-N)
  14. IMPLICIT REAL*8 (A-H, O-Z)
  15. DIMENSION PDILT(1:*)
  16. DIMENSION VEPSTH(1:*)
  17. ISOTH = NINT(PDILT(1))
  18. IP = 1
  19. IF (ISOTH.EQ.1) THEN
  20. C -----------------
  21. C Materiau isotrope
  22. C
  23. ALPHA = PDILT(2)+(PDILT(3)+PDILT(4)*TT)*TT
  24. DO IJ = 1,3
  25. VEPSTH(IJ) = ALPHA*TTP
  26. END DO
  27. ELSE IF (ISOTH.EQ.0) THEN
  28. C -------------------
  29. C Materiau orthotrope
  30. C
  31. DO IJ = 1,3
  32. VEPSTH(IJ) = (PDILT(IP+1)+(PDILT(IP+2)+PDILT(IP+3)*TT)*TT)*TTP
  33. IP = IP+3
  34. END DO
  35. ELSE IF (ISOTH.EQ.-1) THEN
  36. C ---------------------------
  37. C Pas de dilatation thermique
  38. C
  39. DO IJ = 1,3
  40. VEPSTH(IJ) = 0.
  41. END DO
  42. END IF
  43. C ------
  44. DO IJ = 4,6
  45. VEPSTH(IJ) = 0.
  46. END DO
  47. RETURN
  48. END
  49.  
  50.  

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