Télécharger pas_epth.procedur

Retour à la liste

Numérotation des lignes :

  1. * PAS_EPTH PROCEDUR CB215821 16/11/17 21:17:54 9179
  2. 'DEBP' PAS_EPTH PRECED*TABLE MODEVAL*MMODEL MATEVAL*MCHAML
  3. TEVAL*CHPOINT ;
  4. ************************************************************************
  5. * Cette procedure calcule la deformation thermique de reference si
  6. * necessaire et calcule de la deformation thermique du pas courant
  7. *
  8. * PRECED : La TABLE donnee en entree de PASAPAS. On a besoin des indices
  9. * - WTAB.'TEMPERATURE_REFERENCE'
  10. * - WTAB.'TALPHA_REFERENCE'
  11. * - WTAB.'MOD_REF'
  12. * - WTAB.'MAT_REF'
  13. * - WTAB.'ETREF' (On en a besoin et il peut etre calcule ici)
  14. * MODEVAL: Modele pour lequel on veut evaluer EPSILON_THERMIQUE
  15. * MATEVAL: Materiaux pour lequel on veut evaluer EPSILON_THERMIQUE
  16. * TEVAL : Temperature a laquelle on veut evaluer EPSILON_THERMIQUE
  17. *
  18. * Les resultats sont :
  19. * ETHER2 : Deformation thermique dilatation
  20. * MSRTHER2 : Deformation thermique poreuse
  21. ************************************************************************
  22.  
  23. WTAB = PRECED. 'WTABLE';
  24.  
  25. ************************************************************************
  26. * Traitement de la reference : 'ETREF'
  27. ************************************************************************
  28. 'SI' ('EGA' WTAB.'ETAT_REF' ('MOT' 'SUPPRIME'));
  29. * Deformation thermique de reference ETREF calculee si elle n'existe pas
  30. 'SI' ('EGA' WTAB.'ETREF' ('MOT' 'INCONNU')) ;
  31. TREF = WTAB.'TEMPERATURE_REFERENCE' ;
  32. TAREF = WTAB.'TALPHA_REFERENCE' ;
  33. TYREF = 'TYPE' TREF ;
  34. TYAREF = 'TYPE' TAREF ;
  35.  
  36. * On fait en sorte que la difference des temperatures soit realisable
  37. MOREF = WTAB.'MOD_REF';
  38. 'SI'(('EGA' TYREF ('MOT' 'CHPOINT')) 'ET'
  39. ('EGA' TYAREF ('MOT' 'MCHAML' )) ) ;
  40. TREF = 'CHAN' 'CHAM' TREF MOREF 'TEMPERATURES' ;
  41. 'SINO';
  42. 'SI'(('EGA' TYREF ('MOT' 'MCHAML' )) 'ET'
  43. ('EGA' TYAREF ('MOT' 'CHPOINT')) ) ;
  44. TAREF = 'CHAN' 'CHAM' TAREF MOREF 'TEMPERATURES';
  45. 'FINS';
  46. 'FINS';
  47.  
  48. WTAB.'ETREF' ='EPTH' MOREF WTAB.'MAT_REF' (TREF '-' TAREF) ;
  49. WTAB.'MAT_REF' ='MOT' 'SUPPRIME' ;
  50. 'FINS';
  51. 'SINO';
  52. 'SI' ('NEG' WTAB.'ETAT_REF' ('MOT' 'INCONNU'));
  53. MESS 'Dans PAS_EPTH, LOGIQUE fausse';
  54. ERRE 21;
  55. 'FINS';
  56. 'FINS';
  57.  
  58. ************************************************************************
  59. * Traitement courant 'ETHER2'
  60. ************************************************************************
  61. TAREF = WTAB.'TALPHA_REFERENCE' ;
  62. TYAREF = 'TYPE' TAREF ;
  63.  
  64. * On fait en sorte que la difference des temperatures soit realisable
  65. 'SI'('EGA' TYAREF ('MOT' 'MCHAML' )) ;
  66. TEVAL = 'CHAN' 'CHAM' TEVAL MODEVAL 'TEMPERATURES';
  67. 'FINS';
  68.  
  69. ETHER2 = ('EPTH' MODEVAL MATEVAL (TEVAL '-' TAREF)) '-'
  70. ('REDU' WTAB.'ETREF' MODEVAL ) ;
  71.  
  72.  
  73. 'SI' ('EGA' WTAB.'ETHER1' ('MOT' 'INCONNU'));
  74. WTAB.'ETHER1'= ETHER2;
  75. 'FINS';
  76.  
  77.  
  78. ************************************************************************
  79. * Traitement de la formulation 'POREUX'
  80. ************************************************************************
  81. 'SI' WTAB.'POR1';
  82. MOPOR = 'EXTR' MODEVAL 'FORM' 'POREUX' ;
  83. MAPOR = 'REDU' MATEVAL MOPOR ;
  84. ZALM = ('EXCO' 'ALPM' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  85. ZBI = ('EXCO' 'COB' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  86. ZALPH = ('EXCO' 'ALPH' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  87.  
  88. ZFAC = 3.D0 '*' ( ZALM '-' ( ZBI '*' ZALPH ) );
  89. ZFAC = ZFAC 'CHAN' 'TYPE' 'SCALAIRE';
  90.  
  91. ZDT = 'REDU' (TEVAL '-' WTAB.'TEMPERATURE_REFERENCE') WTAB.'MAI_POR';
  92. ZDDT = 'CHAN' 'CHAM' ZDT MOPOR 'STRESSES' 'SCALAIRE';
  93. ZFAC = 'CHAN' 'STRESSES' MOPOR ZFAC ;
  94.  
  95. MSRTHER2 = ZFAC '*' ZDDT ;
  96.  
  97. 'SI' ('EGA' WTAB.'MSRTHER1' ('MOT' 'INCONNU'));
  98. WTAB.'MSRTHER1'= MSRTHER2 ;
  99. 'FINS';
  100. 'SINO';
  101. MSRTHER2 = WTAB.'MSRTHER2';
  102. 'FINS';
  103.  
  104. 'FINP' ETHER2 MSRTHER2;
  105.  
  106.  
  107.  

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