Télécharger castest_lse2_litu.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : castest_lse2_litu.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************
  5. *** OPTIONS GENERALES ***
  6. **********************************************
  7. temps;
  8. *
  9. ** Graph = VRAI ;
  10. Graph = FAUX ;
  11. *
  12. option dime 3 elem seg2;
  13. OPTI EPSI LINEAIRE;
  14. *
  15. 'SI' Graph ;
  16. option trace 'PSC' ;
  17. titre 'MODE PROPRE D UN TUYAU REMPLI DE FLUIDE';
  18. 'FINSI' ;
  19. **********************************************
  20. *** DONNEES GENERALES DU MAILLAGE ***
  21. **********************************************
  22. ndiv = 15;
  23. rac2 = 1.4142136;
  24. x1 = 7*rac2/2.;
  25. x2 = 5*rac2;
  26. p1 = 0 0 0;
  27. p2 = x1 x1 0;
  28. p3 = x2 x2 0;
  29. oeil = 0 0 1000;
  30. *
  31. tuy1 = p1 D ndiv p2;
  32. tuy2 = p2 D ndiv p3;
  33. tuy = tuy1 et tuy2;
  34. *
  35. 'SI' Graph ; trace oeil tuy; 'FINSI' ;
  36. **********************************************
  37. *** DEFINITION DES MODELES ***
  38. **********************************************
  39. obcou1 = mode tuy1 liquide mecanique litu;
  40. obflu1 = mode tuy1 liquide lse2;
  41. obtuy1 = mode tuy1 mecanique tuya;
  42. obcou2 = mode tuy2 liquide mecanique litu;
  43. obflu2 = mode tuy2 liquide lse2;
  44. obtuy2 = mode tuy2 mecanique tuya;
  45. **********************************************
  46. *** CARACTERISTIQUES DU TUYAU ***
  47. **********************************************
  48. * vitesse du son modifie dans le cas E=2E+11
  49. *vit1 = 1355.2619;
  50. *vit2 = 1380.8619;
  51. *etuy= 2.E11;
  52. *--------------------------------------------*
  53. * vitesse du son modifie pour E=2E+08 *
  54. *--------------------------------------------*
  55. vit1 = 99.778516;
  56. vit2 = 111.49412;
  57. etuy= 2.E8;
  58. *--------------------------------------------*
  59. * premier troncon *
  60. *--------------------------------------------*
  61. macou1 = mate obcou1 rho 1000. rorf 1000. lcar 1.
  62. rayo .05;
  63. maflu1 = mate obflu1 rho 1000. rorf 1000. lcar 1.
  64. cson vit1 cref 1500. rayo .05;
  65. matuy1 = mate obtuy1 rho 7800. youn etuy nu .3
  66. rayo .055 epai .005;
  67. *--------------------------------------------*
  68. * deuxieme troncon *
  69. *--------------------------------------------*
  70. macou2 = mate obcou2 rho 1000. rorf 1000. lcar 1.
  71. rayo .08;
  72. maflu2 = mate obflu2 rho 1000. rorf 1000. lcar 1.
  73. cson vit2 cref 1500. rayo .08;
  74. matuy2 = mate obtuy2 rho 7800. youn etuy nu .3
  75. rayo .09 epai .01;
  76. **********************************************
  77. *** RIGIDITES ET MASSES ET COND. LIMITES ***
  78. **********************************************
  79. mas1 = mass obcou1 macou1;
  80. mas2 = mass obflu1 maflu1;
  81. mas3 = mass obtuy1 matuy1;
  82. mas4 = mass obcou2 macou2;
  83. mas5 = mass obflu2 maflu2;
  84. mas6 = mass obtuy2 matuy2;
  85. *
  86. mastot = mas1 et mas2 et mas3 et
  87. mas4 et mas5 et mas6;
  88. *
  89. rig2 = rigi obflu1 maflu1;
  90. rig3 = rigi obtuy1 matuy1;
  91. rig4 = rigi obflu2 maflu2;
  92. rig5 = rigi obtuy2 matuy2;
  93. *--------------------------------------------*
  94. * conditions limites *
  95. *--------------------------------------------*
  96. rig6 = bloq ux uy p1;
  97. rig7 = bloq uz rx ry rz tuy;
  98. rig8 = bloq 'P ' 'PI ' p1;
  99. *--------------------------------------------*
  100. rigtot = rig2 et rig3 et rig4 et rig5 et rig6 et
  101. rig7 et rig8;
  102. **********************************************
  103. *** RESOLUTION ET AFFICHAGE ***
  104. **********************************************
  105. res = vibr inte 1. 2. rigtot mastot tbas impr;
  106. *
  107. fre = res . modes . 1 . frequence;
  108. mess 'MODE & DE TRACTION COMPRESSION THEORIQUE 1.57 CALCULE 'fre;
  109. chsol = res . modes . 1 . deformee_modale;
  110. 'SI' Graph ;
  111. ev1 = evol chpo chsol 'P ' tuy;
  112. dess ev1;
  113. 'FINSI' ;
  114. **********************************************
  115. *** CODE DE BON FONCTIONNEMENT ***
  116. **********************************************
  117. sol = 1.57;
  118. err = 100*((abs(sol-fre)/sol));
  119. si (err < 5);
  120. erre 0;
  121. sinon;
  122. erre 5;
  123. finsi;
  124. *
  125. temps;
  126. fin;
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  

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