Télécharger la plaquette Cast3M


fr en

Télécharger test_repr_modl_ther.dgibi

Retour à la liste

Numérotation des lignes : 

  1. ************************************************************************
  2. *  Cas test :  test_repr_modl_ther.dgibi
  3. *  --------
  4. *
  5. *  Categorie : Verification en cas de changement de modele
  6. *  ---------
  7. *
  8. *  Description :
  9. *  -----------
  10. *  plaque 2D sur laquelle on applique 2 ou 3 modeles selon etapes
  11. *    1- Diffusion de la chaleur
  12. *    2- Source de chaleur
  13. *    3a- Pas de convection a etapes 1 et 3
  14. *    3b- Convection avec exterieur a etape 2
  15. *
  16. ************************************************************************
  17.  
  18. GRAPH = FAUX ;
  19.  
  20. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' 'MODE' 'PLAN' ;
  21. 'OPTION' 'ECHO' 0 ;
  22.  
  23. *************** Geometrie *************************************************
  24. P1  =  0.     0.    ;
  25. P2  =  5.e-3  0.    ;
  26. P3  =  5.e-3 10.e-3 ;
  27. P4  =  0.    10.e-3 ;
  28.  
  29. n   = 20 ;
  30. L23 = P2 'DROI' n  P3 ;
  31. L14 = P1 'DROI' n  P4 ;
  32.  
  33. n   = 10 ;
  34. su1 = L14 'REGL' n L23 ;
  35.  
  36. L_BAS = ('COTE' 4 su1) 'COUL' 'BLEU' ;
  37. mpoi1 ='POIN' su1 'DROI' P1 P2 1.D-10 ;
  38. S_BAS = ('ELEM' su1 'APPU' 'LARG' MPOI1) 'COUL' 'VERT' ;
  39.  
  40. L_HAU = ('COTE' 2 su1) 'COUL' 'ORANGE' ;
  41. mpoi1 = 'POIN' su1 'DROI' P3 P4 1.D-10 ;
  42. S_HAU = ('ELEM' su1 'APPU' 'LARG' MPOI1) 'COUL' 'ROUG' ;
  43.  
  44. **************** Modeles et caracteristiques ***********************************
  45. modthe ='MODE' su1   'THERMIQUE' 'CONDUCTION' 'CONS' 'CON1' ;
  46. moconv ='MODE' L_HAU 'THERMIQUE' 'CONVECTION' 'CONS' 'COV1' ;
  47.  
  48. rho1   = 4500.    ; cp1    = 900.     ;
  49.  
  50. matthe = 'MATE' modthe 'C'  cp1   'K' 10. 'RHO' rho1 ;
  51. maconv = 'MATE' moconv 'H   ' 150. ;
  52.  
  53. *************** Pas de temps ***************************************************
  54. listt1 = prog 0. 'PAS' 10. 1000. ;
  55. listt2 = prog (1000.+20.) 'PAS' 20. 3000. ;
  56. listt3 = prog (3000.+10.) 'PAS' 10. 4000. ;
  57.  
  58. T_ini = 'MANU' 'CHPO' su1 1 'T'  0. 'NATURE' 'DIFFUS' ;
  59. 'SI' GRAPH ;
  60.   'TRAC' su1 T_ini 'TITR' 'Temperature initiale' ;
  61. 'FINS';
  62.  
  63. T_air = 25. ;
  64. EvTECO = 'EVOL' 'MANU' 'TEMPS' ('PROG' 0. 1.E+9)
  65.                        'T'     ('PROG' 1. 1.)  ;
  66. ChpTECO  = 'MANU' 'CHPO' L_HAU 1 'T' T_air      ;
  67. CharTECO = 'CHARGEMENT' 'TECO' ChpTECO EvTECO  ;
  68.  
  69. * CHARGEMENT
  70. srcth1  = 'SOURCE' ('REDU' modthe s_bas) s_bas 1.D7 ;
  71. EvS    = EVOL MANU 'TEMP' (PROG 0. 1.E+9) 'Q   ' (PROG 1. 1.) ;
  72. charths ='CHAR' 'Q'  srcth1 EvS ;
  73.  
  74. ltps = 'PROG'0. ; lT1 = 'PROG' 0. ; lT3 = 'PROG' 0. ;
  75.  
  76. * ETAPE 1 : On chauffe la piece fermee (adiabatiquement)
  77. xtab                       = TABL    ;
  78. xtab.'MODELE'              = modthe  ;
  79. xtab.'CARACTERISTIQUES'    = matthe  ;
  80. xtab.'CHARGEMENT'          = charths  ;
  81. xtab.'TEMPS_CALCULES'      = listt1  ;
  82.  
  83. xtab.'TEMPERATURES'        = table   ;
  84. xtab.'TEMPERATURES' . 0    = T_ini   ;
  85.  
  86. xtab.'PROCESSEURS'         = 'MOT' 'MONO_PROCESSEUR' ;
  87.  
  88. PASAPAS xtab ;
  89.  
  90. * Post-traitement fin etape 1
  91. II = ('DIMENSION' xtab.'TEMPS') - 1 ;
  92. ltps = ltps 'ET' xtab.'TEMPS'. II ;
  93. chp_T = 'EXCO' 'T   ' xtab.'TEMPERATURES'. II 'T   ' ;
  94. lT1 = lT1 'ET' ('EXTR' chp_T 'T' P1) ;
  95. lT3 = lT3 'ET' ('EXTR' chp_T 'T' P3) ;
  96.  
  97. * ETAPE 2 : On chauffe la piece qui a une fenetre ouverte (convection avec exterieure)
  98.  
  99. xtab.'MODELE'              = modthe 'ET' moconv  ;
  100. xtab.'CARACTERISTIQUES'    = matthe 'ET' maconv ;
  101. xtab.'CHARGEMENT'          = charths 'ET' charteco ;
  102. xtab.'TEMPS_CALCULES'      = listt2 ;
  103.  
  104. PASAPAS xtab ;
  105.  
  106. * Post-traitement fin etape 2
  107. II = ('DIMENSION' xtab.'TEMPS') - 1 ;
  108. ltps = ltps 'ET' xtab.'TEMPS'. II ;
  109. chp_T = 'EXCO' 'T   ' xtab.'TEMPERATURES'. II 'T   ' ;
  110. lT1 = lT1 'ET' ('EXTR' chp_T 'T' P1) ;
  111. lT3 = lT3 'ET' ('EXTR' chp_T 'T' P3) ;
  112.  
  113. * ETAPE 3 : On chauffe a nouveau la piece fermee (adiabatiquement)
  114.  
  115. xtab.'MODELE'              = modthe ;
  116. xtab.'CARACTERISTIQUES'    = matthe ;
  117. xtab.'CHARGEMENT'          = charths ;
  118. xtab.'TEMPS_CALCULES'      = listt3 ;
  119.  
  120. PASAPAS xtab ;
  121.  
  122. * Post-traitement fin etape 3
  123. II = ('DIMENSION' xtab.'TEMPS') - 1 ;
  124. ltps = ltps 'ET' xtab.'TEMPS'. II ;
  125. chp_T = 'EXCO' 'T   ' xtab.'TEMPERATURES'. II 'T   ' ;
  126. lT1 = lT1 'ET' ('EXTR' chp_T 'T' P1) ;
  127. lT3 = lT3 'ET' ('EXTR' chp_T 'T' P3) ;
  128.  
  129. * POST-TRAITEMENT FINAL
  130. *'OPTI' 'EPTR' 10 ;
  131.  
  132. 'SI' GRAPH ;
  133.   ev_p1 = EVOL 'ROUG' 'TEMP' xtab 'TEMPERATURES' 'T' P1 ;
  134.   ev_p3 = EVOL 'VERT' 'TEMP' xtab 'TEMPERATURES' 'T' P3 ;
  135.   dess (ev_p1 et ev_p3) ;
  136. 'FINSI' ;
  137.  
  138. lister ltps ; lister lT1 ; lister lT3 ;
  139.  
  140. err = 'MAXIMUM' (lT3 '-' ('PROG' 0. 122.626 58.399 183.524)) 'ABS' ;
  141. list err ;
  142. 'SI' ('>EG' err 1.E-3) ;
  143.   'MESS' ;
  144.   'MESS' 'ERREUR : Temperature(s) au point 3 erronee(s)' ;
  145.   'ERREUR' 5 ;
  146. 'FINSI' ;
  147.  
  148. 'FIN';
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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