Télécharger flam3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : flam3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * COMPARAISON DES CALCULS DE FLAMBAGE D'UN TUBE
  5. * SOUS PRESSION INTERNE (TUBE DE SACLAY)
  6. OPTI DIME 3 ELEM QUA4 ECHO 1;
  7. MESS ' COMPARAISON DES CALCULS DE FLAMBAGE';
  8. MESS 'D UN TUBE SOUS PRESSION INTERNE MODéLISé';
  9. MESS ' AVEC COQ3 ET COQ4 (TUBE DE SACLAY)';
  10. L = 10. ;
  11. T = .01 ;
  12. R = .5 ;
  13. E = 2E11 ;
  14. NV = 30 ;
  15. NH = 8 ;
  16. ************************
  17. P1 = R 0 0 ;
  18. P2 = R 0 L ;
  19. L1 = P1 D NV P2 ;
  20. AX1 = 0 0 0 ;
  21. AX2 = 0 0 L;
  22. SURF4 = ROTA L1 AX1 AX2 NH 90. ;
  23. SURF3 = CHAN 'TRI3' SURF4;
  24. SURF8 = CHAN 'QUADRATIQUE' SURF4 ;
  25. * MODéLES, MATéRIAUX ET RIGIDITé
  26. MODL3 = MODE SURF3 MECANIQUE COQ3;
  27. MODL4 = MODE SURF4 MECANIQUE COQ4;
  28. MODL8 = MODE SURF8 MECANIQUE COQ8;
  29. MAT3 = MATE MODL3 YOUN E NU .30 RHO 8000 EPAI T ;
  30. MAT4 = MATE MODL4 YOUN E NU .30 RHO 8000 EPAI T ;
  31. MAT8 = MATE MODL8 YOUN E NU .30 RHO 8000 EPAI T ;
  32. *
  33. MOP1 = 'MODE' SURF3 'CHARGEMENT' 'PRESSION' COQ3 ;
  34. MOP2 = 'MODE' SURF4 'CHARGEMENT' 'PRESSION' COQ4 ;
  35. MOP3 = 'MODE' SURF8 'CHARGEMENT' 'PRESSION' COQ8 ;
  36. MAP1 = 'PRES' MOP1 'PRES' -100. ;
  37. MAP2 = 'PRES' MOP2 'PRES' -100. ;
  38. MAP3 = 'PRES' MOP3 'PRES' -100. ;
  39. MAP3B = 'CARA' MOP3 'EPAI' T ;
  40. MODL3 = MODL3 ET MOP1 ; MAT3 = MAT3 ;
  41. MODL4 = MODL4 ET MOP2 ; MAT4 = MAT4 ;
  42. MODL8 = MODL8 ET MOP3 ; MAT8 = MAT8 ;
  43. *
  44. RIG3 = RIGI MODL3 MAT3 ;
  45. RIG4 = RIGI MODL4 MAT4 ;
  46. RIG8 = RIGI MODL8 MAT8 ;
  47. I1 = PI*(R**3)*T ;
  48. S1 = 2*PI*R*T ;
  49. * CONDITION AUX LIMITES
  50. LL41 = COTE 2 SURF4 ;
  51. BL41 = SYMT DEPL ROTA (0 0 0) (1 0 0) (0 1 0) SURF4 .01 ;
  52. BL42 = SYMT DEPL ROTA AX1 AX2 P1 SURF4 .01 ;
  53. BL43 = SYMT DEPL ROTA AX1 AX2 (0 1 0) SURF4 .01;
  54. F4 = 'BSIG' MOP2 MAP2 ;
  55. LL31 = COTE 2 SURF4 ;
  56. BL31 = SYMT DEPL ROTA (0 0 0) (1 0 0) (0 1 0) SURF3 .01 ;
  57. BL32 = SYMT DEPL ROTA AX1 AX2 P1 SURF3 .01 ;
  58. BL33 = SYMT DEPL ROTA AX1 AX2 (0 1 0) SURF3 .01;
  59. F3 = 'BSIG' MOP1 MAP1 ;
  60. BL81 = SYMT DEPL ROTA (0 0 0) (1 0 0) (0 1 0) SURF8 .01 ;
  61. BL82 = SYMT DEPL ROTA AX1 AX2 P1 SURF8 .01 ;
  62. BL83 = SYMT DEPL ROTA AX1 AX2 (0 1 0) SURF8 .01;
  63. F8 = 'BSIG' MOP3 MAP3 MAP3B ;
  64. * CALCULS éLASTIQUES
  65. DEP3 = RESO (RIG3 ET BL31 ET BL32 ET BL33) F3 ;
  66. SIG3 = SIGM MODL3 MAT3 DEP3 ;
  67. KSIG3 = KSIG MODL3 MAT3 SIG3 FLAM ;
  68. DEP4 = RESO (RIG4 ET BL41 ET BL42 ET BL43) F4 ;
  69. SIG4 = SIGM MODL4 MAT4 DEP4 ;
  70. KSIG4 = KSIG MODL4 MAT4 SIG4 FLAM ;
  71. DEP8 = RESO (RIG8 ET BL81 ET BL82 ET BL83) F8 ;
  72. SIG8 = SIGM MODL8 MAT8 DEP8 ;
  73. KSIG8 = KSIG MODL8 MAT8 SIG8 FLAM ;
  74. * CONDITIONS AUX LIMITES POUR LE CALCUL DE FLAMBAGE
  75. BL33 = ANTI DEPL ROTA AX1 AX2 (0 1 0) SURF3 .01;
  76. PP0 = POIN SURF3 PROCHE (0 R L) ;
  77. BL34 = BLOQ UX PP0 ;
  78. RIGT3 = RIG3 ET BL31 ET BL32 ET BL33 ET BL34 ;
  79. BL43 = ANTI DEPL ROTA AX1 AX2 (0 1 0) SURF4 .01;
  80. PP0 = POIN SURF4 PROCHE (0 R L) ;
  81. BL44 = BLOQ UX PP0 ;
  82. RIGT4 = RIG4 ET BL41 ET BL42 ET BL43 ET BL44 ;
  83. BL83 = ANTI DEPL ROTA AX1 AX2 (0 1 0) SURF8 .01;
  84. PP0 = POIN SURF8 PROCHE (0 R L) ;
  85. BL84 = BLOQ UX PP0 ;
  86. RIGT8 = RIG8 ET BL81 ET BL82 ET BL83 ET BL84 ;
  87. * CALCUL DE FLAMBAGE
  88. KP3 = KP MOP1 (MANU CHPO SURF3 1 P -100.) FLAM ;
  89. TAB3 = VIBR PROCHE (PROG 79.) RIGT3
  90. ((1*KP3) ET (-1*KSIG3) ) ;
  91. KP4 = KP MOP2 (MANU CHPO SURF4 1 P -100.) FLAM ;
  92. TAB4 = VIBR PROCHE (PROG 79.) RIGT4
  93. ((1*KP4) ET (-1*KSIG4) ) ;
  94. KP8 = KP MOP3 (MANU CHPO SURF8 1 P -100.) FLAM ;
  95. TAB8 = VIBR PROCHE (PROG 79.) RIGT8
  96. ((1*KP8) ET (-1*KSIG8) ) ;
  97. DEP3 = TAB3.MODES. 1 .DEFORMEE_MODALE ;
  98. DEP4 = TAB4.MODES. 1 .DEFORMEE_MODALE ;
  99. DEP8 = TAB8.MODES. 1 .DEFORMEE_MODALE ;
  100. DEF3 = DEFO DEP3 SURF3 1. ROUGE ;
  101. DEF4 = DEFO DEP4 SURF4 1. VERT ;
  102. DEF8 = DEFO DEP8 SURF8 1. BLEU ;
  103. DEF0 = DEFO DEP3 SURF3 0. ;
  104. *TRAC (DEF3 ET DEF4 ET DEF0) ;
  105. F3 = TAB3.MODES. 1 .FREQUENCE ;
  106. LL3 = (2*PI*F3)**2 ;
  107. F4 = TAB4.MODES. 1 .FREQUENCE ;
  108. LL4 = (2*PI*F4)**2 ;
  109. F8 = TAB8.MODES. 1 .FREQUENCE ;
  110. LL8 = (2*PI*F8)**2 ;
  111. PCOQ3 = LL3*100 ;
  112. PCOQ4 = LL4*100 ;
  113. PCOQ8 = LL8*100 ;
  114. PTH = PI**2*E*R*T/4/L/L ;
  115. ERR3 = ABS (100.*(PTH - PCOQ3)/PCOQ3);
  116. ERR4 = ABS (100.*(PTH - PCOQ4)/PCOQ4);
  117. ERR8 = ABS (100.*(PTH - PCOQ8)/PCOQ8);
  118. MESS 'CALCUL THéORIQUE : ' PTH;
  119. MESS 'CALCUL NUMéRIQUE (COQ3) : ' PCOQ3;
  120. MESS 'ERREUR : ' ERR3 '%';
  121. MESS 'CALCUL NUMéRIQUE (COQ4) : ' PCOQ4;
  122. MESS 'ERREUR : ' ERR4 '%';
  123. MESS 'CALCUL NUMéRIQUE (COQ8) : ' PCOQ8;
  124. MESS 'ERREUR : ' ERR8 '%';
  125. *------------------ CODE DE BON FONCTIONNEMENT -------------------------
  126. SI ((ERR3 < 5) ET (ERR4 < 5) ET (ERR8 < 8));
  127. ERRE 0;
  128. SINON;
  129. ERRE 5;
  130. FINSI;
  131. FIN;
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  

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