Télécharger la plaquette Cast3M


fr en

Télécharger fuite_fissure.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  fuite_fissure.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *   Test fuite_fissure.dgibi: Jeux de données         *
  5. *    ----------------------------------------         *
  6. *                                                     *
  7. *******************************************************
  8. *****   CALCUL DU DEBIT DE FUITE D'UN             *****
  9. *****   MELANGE D'AIR SEC                         *****
  10. *****   A TRAVERS UNE FISSURE TRAVERSANTE         *****
  11. *****   POUR UNE DIFFERENCE DE PRESSION IMPOSEE   *****
  12. *****   écoulement monodimensionnel, diphasique   *****
  13. *****   homogène, en régime permanent             *****
  14. *******************************************************
  15. *****   AUTEUR   : HELENE SIMON                   *****
  16. *****   DATE     : 26/10/2005                     *****
  17. *******************************************************
  18. opti echo 1 dime 2 elem qua4 ;
  19. *une seule ligne : la fissure est parallele au cote l12
  20. l = 1./(2.**0.5) ;
  21. p1=0. 0. ; p2= l   l ;
  22. l12 = d p1 p2 dini 0.01 dfin 0.01  ;
  23. *
  24. ouv = 'PROG' (nbno l12)*200e-6 ;
  25. ep = manu 'CHPO' (chan 'POI1' l12) 1 'OUV' ouv natu 'DIFFUS';
  26. tp = manu chpo l12  1 'T' 120. natu 'DIFFUS' ;
  27. *modele
  28. MO1 = 'MODELISER' l12 FISSURE MASS PARF POISEU_COLEBROOK;
  29. *champ de prop materielles
  30. MAT1 = 'MATERIAU' MO1 'RUGO' 40e-6  ;
  31. *conditions aux limites
  32. TAB_CL = table ;
  33. TAB_CL.'PRESSION_TOTALE_AMONT' = 5e5 ;
  34. TAB_CL.'PRESSION_VAPEUR_AMONT' = 0. ;
  35. TAB_CL.'TEMPERATURE_AMONT' = 170. ;
  36. TAB_CL.'PRESSION_TOTALE_AVAL' = 1e5 ;
  37. TAB_CL.'TEMPERATURE_PAROI' = tp ;
  38. TAB_CL.'OUVERTURE' = ep ;
  39. *
  40. res = fiss MO1 MAT1 TAB_CL ;
  41. *
  42. Qn = 'EXTRAIRE' res 'Q' P2 ;
  43. Qref = 1.74459E-02 ;
  44. 'MESSAGE' 'debit calcule = ' Qn 'debit reference = ' Qref;
  45. ere1=(Qn-Qref)/Qref ;
  46. ere='ABS' ere1 ;
  47. SI (ere  <EG 1.d-5 ) ;
  48.    ERRE 0;
  49. SINON;
  50.    ERRE 5;
  51. FINSI;
  52. 'FIN' ;
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  

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