Télécharger pas_epth.procedur

Retour à la liste

Numérotation des lignes :

  1. * PAS_EPTH PROCEDUR CB215821 18/11/26 21:15:05 10010
  2. 'DEBP' PAS_EPTH PRECED*TABLE MODEVAL*MMODEL MATEVAL*MCHAML
  3. TCHP/CHPOINT TCHAM/MCHAML ;
  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. *
  17. * TCHP : Temperature CHPOINT a laquelle est evalué EPSILON_THERMIQUE
  18. * ou
  19. * TCHAM : Temperature MCHAML a laquelle est evalué EPSILON_THERMIQUE
  20. * Obligatoire pour les JOINTS
  21. * Necessaire pour certains éléments linéiques
  22. *
  23. * Les resultats sont :
  24. * ETHER2 : Deformation thermique dilatation
  25. * MSRTHER2 : Deformation thermique poreuse
  26. ************************************************************************
  27.  
  28. WTAB = PRECED. 'WTABLE';
  29.  
  30. ************************************************************************
  31. * Traitement de la reference : 'ETREF'
  32. ************************************************************************
  33. 'SI' ('EGA' WTAB.'ETAT_REF' 'SUPPRIME');
  34. * Deformation thermique de reference ETREF calculee si elle n'existe pas
  35. 'SI' ('EGA' WTAB.'ETREF' 'INCONNU') ;
  36. TREF = WTAB.'TEMPERATURE_REFERENCE';
  37. TAREF = WTAB.'TALPHA_REFERENCE' ;
  38. TYREF ='TYPE' TREF ;
  39. TYAREF ='TYPE' TAREF ;
  40.  
  41. * On fait en sorte que la difference des temperatures soit realisable
  42. MOREF = WTAB.'MOD_REF';
  43. 'SI'(('EGA' TYREF 'CHPOINT') 'ET'
  44. ('EGA' TYAREF 'MCHAML' )) ;
  45. TREF = 'CHAN' 'CHAM' TREF MOREF 'TEMPERATURES' ;
  46. 'SINO';
  47. 'SI'(('EGA' TYREF 'MCHAML' ) 'ET'
  48. ('EGA' TYAREF 'CHPOINT')) ;
  49. TAREF = 'CHAN' 'CHAM' TAREF MOREF 'TEMPERATURES';
  50. 'FINS';
  51. 'FINS';
  52.  
  53. DELTAT = TREF '-' TAREF ;
  54. TYDT ='TYPE' DELTAT;
  55. * Protection si Type FLOTTANT
  56. 'SI' ('EGA' TYDT 'FLOTTANT');
  57. DELTAT = (('ZERO' MOREF 'TEMPERATURES') ** 0) * DELTAT;
  58. 'FINS';
  59.  
  60. WTAB.'ETREF' ='EPTH' MOREF WTAB.'MAT_REF' DELTAT ;
  61. WTAB.'MAT_REF' ='MOT' 'SUPPRIME' ;
  62. 'FINS';
  63. 'SINO';
  64. 'SI' ('NEG' WTAB.'ETAT_REF' 'INCONNU' );
  65. MESS 'Dans PAS_EPTH, LOGIQUE fausse';
  66. ERRE 21;
  67. 'FINS';
  68. 'FINS';
  69.  
  70. ************************************************************************
  71. * Traitement courant 'ETHER2'
  72. ************************************************************************
  73. 'SI' ('EGA' ('TYPE' TCHP) 'CHPOINT');
  74. TEVAL = TCHP;
  75. TYPT = 'MOT' 'CHPOINT';
  76. 'SINO';
  77. 'SI' ('EGA' ('TYPE' TCHAM) 'MCHAML');
  78. TEVAL = TCHAM;
  79. TYPT = 'MOT' 'MCHAML';
  80. 'SINO';
  81. ERRE 37 'AVEC' 'CHPOINT';
  82. 'FINS';
  83. 'FINS';
  84. TAREF = WTAB.'TALPHA_REFERENCE' ;
  85. TYAREF = 'TYPE' TAREF ;
  86.  
  87. *On fait en sorte que la difference des temperatures soit realisable
  88. 'SI'(('EGA' TYPT 'CHPOINT') 'ET'
  89. ('EGA' TYAREF 'MCHAML' )) ;
  90. TEVAL = 'CHAN' 'CHAM' TEVAL MOREF 'TEMPERATURES' ;
  91. 'SINO';
  92. 'SI'(('EGA' TYPT 'MCHAML' ) 'ET'
  93. ('EGA' TYAREF 'CHPOINT')) ;
  94. TAREF = 'CHAN' 'CHAM' TAREF MOREF 'TEMPERATURES';
  95. 'FINS';
  96. 'FINS';
  97.  
  98. ETHER2 = ('EPTH' MODEVAL MATEVAL (TEVAL '-' TAREF)) '-'
  99. ('REDU' WTAB.'ETREF' MODEVAL ) ;
  100.  
  101.  
  102. 'SI' ('EGA' WTAB.'ETHER1' 'INCONNU');
  103. WTAB.'ETHER1'= ETHER2;
  104. 'FINS';
  105.  
  106.  
  107. ************************************************************************
  108. * Traitement de la formulation 'POREUX'
  109. ************************************************************************
  110. 'SI' WTAB.'POR1';
  111. MOPOR = 'EXTR' MODEVAL 'FORM' 'POREUX' ;
  112. MAPOR = 'REDU' MATEVAL MOPOR ;
  113. ZALM = ('EXCO' 'ALPM' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  114. ZBI = ('EXCO' 'COB' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  115. ZALPH = ('EXCO' 'ALPH' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  116.  
  117. ZFAC = 3.D0 '*' ( ZALM '-' ( ZBI '*' ZALPH ) );
  118. ZFAC = ZFAC 'CHAN' 'TYPE' 'SCALAIRE';
  119.  
  120. ZDT = 'REDU' (TEVAL '-' WTAB.'TEMPERATURE_REFERENCE') WTAB.'MAI_POR';
  121. ZDDT = 'CHAN' 'CHAM' ZDT MOPOR 'STRESSES' 'SCALAIRE';
  122. ZFAC = 'CHAN' 'STRESSES' MOPOR ZFAC ;
  123.  
  124. MSRTHER2 = ZFAC '*' ZDDT ;
  125.  
  126. 'SI' ('EGA' WTAB.'MSRTHER1' 'INCONNU');
  127. WTAB.'MSRTHER1'= MSRTHER2 ;
  128. 'FINS';
  129. 'SINO';
  130. MSRTHER2 = WTAB.'MSRTHER2';
  131. 'FINS';
  132.  
  133. 'FINP' ETHER2 MSRTHER2;
  134.  
  135.  
  136.  
  137.  

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