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