Télécharger pas_epth.procedur

Retour à la liste

Numérotation des lignes :

  1. * PAS_EPTH PROCEDUR CB215821 17/05/22 21:15:05 9437
  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. WTAB.'ETREF' ='EPTH' MOREF WTAB.'MAT_REF' (TREF '-' TAREF) ;
  54. WTAB.'MAT_REF' ='MOT' 'SUPPRIME' ;
  55. 'FINS';
  56. 'SINO';
  57. 'SI' ('NEG' WTAB.'ETAT_REF' 'INCONNU' );
  58. MESS 'Dans PAS_EPTH, LOGIQUE fausse';
  59. ERRE 21;
  60. 'FINS';
  61. 'FINS';
  62.  
  63. ************************************************************************
  64. * Traitement courant 'ETHER2'
  65. ************************************************************************
  66. 'SI' ('EGA' ('TYPE' TCHP) 'CHPOINT');
  67. TEVAL = TCHP;
  68. TYPT = 'MOT' 'CHPOINT';
  69. 'SINO';
  70. 'SI' ('EGA' ('TYPE' TCHAM) 'MCHAML');
  71. TEVAL = TCHAM;
  72. TYPT = 'MOT' 'MCHAML';
  73. 'SINO';
  74. ERRE 37 'AVEC' 'CHPOINT';
  75. 'FINS';
  76. 'FINS';
  77. TAREF = WTAB.'TALPHA_REFERENCE' ;
  78. TYAREF = 'TYPE' TAREF ;
  79.  
  80. *On fait en sorte que la difference des temperatures soit realisable
  81. 'SI'(('EGA' TYPT 'CHPOINT') 'ET'
  82. ('EGA' TYAREF 'MCHAML' )) ;
  83. TEVAL = 'CHAN' 'CHAM' TEVAL MOREF 'TEMPERATURES' ;
  84. 'SINO';
  85. 'SI'(('EGA' TYPT 'MCHAML' ) 'ET'
  86. ('EGA' TYAREF 'CHPOINT')) ;
  87. TAREF = 'CHAN' 'CHAM' TAREF MOREF 'TEMPERATURES';
  88. 'FINS';
  89. 'FINS';
  90.  
  91. ETHER2 = ('EPTH' MODEVAL MATEVAL (TEVAL '-' TAREF)) '-'
  92. ('REDU' WTAB.'ETREF' MODEVAL ) ;
  93.  
  94.  
  95. 'SI' ('EGA' WTAB.'ETHER1' 'INCONNU');
  96. WTAB.'ETHER1'= ETHER2;
  97. 'FINS';
  98.  
  99.  
  100. ************************************************************************
  101. * Traitement de la formulation 'POREUX'
  102. ************************************************************************
  103. 'SI' WTAB.'POR1';
  104. MOPOR = 'EXTR' MODEVAL 'FORM' 'POREUX' ;
  105. MAPOR = 'REDU' MATEVAL MOPOR ;
  106. ZALM = ('EXCO' 'ALPM' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  107. ZBI = ('EXCO' 'COB' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  108. ZALPH = ('EXCO' 'ALPH' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  109.  
  110. ZFAC = 3.D0 '*' ( ZALM '-' ( ZBI '*' ZALPH ) );
  111. ZFAC = ZFAC 'CHAN' 'TYPE' 'SCALAIRE';
  112.  
  113. ZDT = 'REDU' (TEVAL '-' WTAB.'TEMPERATURE_REFERENCE') WTAB.'MAI_POR';
  114. ZDDT = 'CHAN' 'CHAM' ZDT MOPOR 'STRESSES' 'SCALAIRE';
  115. ZFAC = 'CHAN' 'STRESSES' MOPOR ZFAC ;
  116.  
  117. MSRTHER2 = ZFAC '*' ZDDT ;
  118.  
  119. 'SI' ('EGA' WTAB.'MSRTHER1' 'INCONNU');
  120. WTAB.'MSRTHER1'= MSRTHER2 ;
  121. 'FINS';
  122. 'SINO';
  123. MSRTHER2 = WTAB.'MSRTHER2';
  124. 'FINS';
  125.  
  126. 'FINP' ETHER2 MSRTHER2;
  127.  
  128.  
  129.  

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