Télécharger la plaquette Cast3M


fr en

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

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