Télécharger flam2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : flam2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5.  
  6. ************************************************************************
  7. * COMPARAISON DES CALCULS DE FLAMBAGE D'UN TUBE
  8. * SOUS PRESSION EXTERNE (D'APRèS R.J. GIBERT)
  9. ************************************************************************
  10.  
  11. *------------------- OPTIONS -----------------------
  12.  
  13. * OPTI DIME 3 ELEM QUA4 ECHO 0;
  14. OPTI DIME 3 ELEM QUA4 ;
  15. MESS ' COMPARAISON DES CALCULS DE FLAMBAGE';
  16. MESS 'D UN TUBE SOUS PRESSION EXTERNE MODéLISé';
  17. MESS ' AVEC COQ3 ET COQ4';
  18.  
  19. *------------- DEFINITION DU MAILLAGE -----------------
  20.  
  21. P0 = 0. 0. 0.;
  22. P1 = 0. 1. 0.;
  23. P2 = 0. 0. 1.;
  24. P3 = 0.08 0. 0.;
  25. P4 = 0.08 1. 0.;
  26. P5 = 0.08 0. 1.;
  27. C1 = CERC 20 P1 P0 P2;
  28. C2 = CERC 20 P4 P3 P5;
  29. L1 = D 1 P1 P4;
  30. L2 = D 1 P2 P5;
  31. OBJ4 = C1 TRAN 1 P3;
  32. ELIM (OBJ4 ET L1 ET L2 ET C1 ET C2) 1.E-3;
  33. OBJ3 = CHAN 'TRI3' OBJ4;
  34.  
  35. *------------- CONDITIONS AUX LIMITES -------------
  36.  
  37. BC1 = BLOQ C1 'UX' 'RY' 'RZ';
  38. BC2 = BLOQ C2 'UX' 'RY' 'RZ';
  39. * BL1 = BLOQ L1 'UZ' 'RX' 'RY';
  40. * BL2 = BLOQ L2 'UY' 'RX' 'RZ';
  41. BL1 = BLOQ L1 'UZ' 'RX';
  42. BL2 = BLOQ L2 'UY' 'RX';
  43. CL = BC1 ET BC2 ET BL1 ET BL2;
  44.  
  45. *------------- MODELE, MATERIAU, MATRICES -------------
  46.  
  47. * MODELE
  48. MOD3 = MODE OBJ3 MECANIQUE ELASTIQUE ISOTROPE COQ3;
  49. * MOD3 = MODE OBJ3 MECANIQUE ELASTIQUE ISOTROPE DKT;
  50. MOD4 = MODE OBJ4 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  51. *
  52. MOP1 = 'MODE' OBJ3 'CHARGEMENT' 'PRESSION' COQ3 ;
  53. MOP2 = 'MODE' OBJ4 'CHARGEMENT' 'PRESSION' COQ4 ;
  54. *
  55. * MATERIAU
  56. MAT3 = MATE MOD3 YOUN 2.E11 NU 0.3 EPAI 0.05;
  57. MAT4 = MATE MOD4 YOUN 2.E11 NU 0.3 EPAI 0.05;
  58. *
  59. MAP1 = 'PRES' MOP1 'PRES' -1000000. ;
  60. MAP2 = 'PRES' MOP2 'PRES' -1000000. ;
  61. *
  62. MOD3 = MOD3 ET MOP1 ; MAT3 = MAT3 ;
  63. MOD4 = MOD4 ET MOP2 ; MAT4 = MAT4 ;
  64. * RIGIDITé
  65. RIG3 = RIGI MOD3 MAT3;
  66. RIG4 = RIGI MOD4 MAT4;
  67.  
  68. *------------- CALCUL ELASTIQUE DES PRECONTRAINTES -------------
  69.  
  70. * CHARGEMENT
  71. PR3 = 'BSIG' MOP1 MAP1 ;
  72. PR4 = 'BSIG' MOP2 MAP2 ;
  73. *
  74. * CALCUL ELASTIQUE
  75. DEP3 = RESO (RIG3 ET CL) PR3;
  76. DEP4 = RESO (RIG4 ET CL) PR4;
  77.  
  78. * CALCUL DES CONTRAINTES
  79. SIG3 = SIGM MAT3 MOD3 DEP3;
  80. SIG4 = SIGM MAT4 MOD4 DEP4;
  81.  
  82. * MATRICE KS ET KP
  83. CHP3 = MANU 'CHPO' OBJ3 1 'SCAL' -1000000.;
  84. CHP4 = MANU 'CHPO' OBJ4 1 'SCAL' -1000000.;
  85. KS3 = KSIG SIG3 MOD3 mat3 FLAM;
  86. KS4 = KSIG SIG4 MOD4 mat4 FLAM;
  87. KP3 = KP MOP1 CHP3 FLAM;
  88. KP4 = KP MOP2 CHP4 FLAM;
  89. RAI3 = KP3 ET ((-1.)*KS3);
  90. RAI4 = KP4 ET ((-1.)*KS4);
  91. RAS3 = (-1.)*KS3;
  92. RAS4 = (-1.)*KS4;
  93.  
  94. *------------- RECHERCHE DE LA CHARGE CRITIQUE -------------
  95. wshift = 0.4;
  96. *REM: difficulte a converger pour VIBR IRAM pour wshift < 0.2 ...;
  97.  
  98. * calcul via VIBR PROCHE :
  99. * en prenant en compte Ksigma et KP :
  100. * wshift = 0.417;
  101. XMOD3 = VIBR PROCHE (PROG wshift) (RIG3 ET CL) RAI3 ;
  102. XMOD4 = VIBR PROCHE (PROG wshift) (RIG4 ET CL) RAI4 ;
  103. * en prenant en compte Ksigma seulement :
  104. * wshift = 0.482;
  105. XMOD3S = VIBR PROCHE (PROG wshift) (RIG3 ET CL) RAS3 ;
  106. XMOD4S = VIBR PROCHE (PROG wshift) (RIG4 ET CL) RAS4 ;
  107.  
  108. * calcul via VIBR IRAM :
  109. * en prenant en compte Ksigma et KP :
  110. * opti impi 5;
  111. TMOD3 = VIBR 'IRAM' wshift 1 (RIG3 ET CL) RAI3 ;
  112. TMOD4 = VIBR 'IRAM' wshift 1 (RIG4 ET CL) RAI4 ;
  113. * en prenant en compte Ksigma seulement :
  114. TMOD3S = VIBR 'IRAM' wshift 1 (RIG3 ET CL) RAS3 ;
  115. TMOD4S = VIBR 'IRAM' wshift 1 (RIG4 ET CL) RAS4 ;
  116. * opti impi 0;
  117.  
  118. * POST-TRAITEMENT :
  119. * solution theorique :
  120. XTH = ((2.E11/(4*(1 - 0.09)))*((0.05/1.)**3))/1000000.;
  121. * solutions calculees :
  122. FF3 = XMOD3.'MODES'. 1 .'FREQUENCE';
  123. FF4 = XMOD4.'MODES'. 1 .'FREQUENCE';
  124. FF3S = XMOD3S.'MODES'. 1 .'FREQUENCE';
  125. FF4S = XMOD4S.'MODES'. 1 .'FREQUENCE';
  126. IFF3 = TMOD3.'MODES'. 1 .'FREQUENCE';
  127. IFF4 = TMOD4.'MODES'. 1 .'FREQUENCE';
  128. IFF3S = TMOD3S.'MODES'. 1 .'FREQUENCE';
  129. IFF4S = TMOD4S.'MODES'. 1 .'FREQUENCE';
  130. X2PI = 2.D0 * PI;
  131. * solutions coq3 :
  132. XCOQ3 = (X2PI * (PROG FF3 IFF3 FF3S IFF3S)) ** 2;
  133. * solutions coq4 :
  134. XCOQ4 = (X2PI * (PROG FF4 IFF4 FF4S IFF4S)) ** 2;
  135. * erreur relative
  136. ERR3 = ABS (100*((extr XCOQ3 (lect 1 2)) - XTH)/XTH);
  137. ERR4 = ABS (100*((extr XCOQ4 (lect 1 2)) - XTH)/XTH);
  138.  
  139. MESS '--------------------------------------------------------------';
  140. MESS '--- VALEUR THEORIQUE : ' XTH;
  141. MESS '--------------------------------------------------------------';
  142. MESS '--- COQ3 : PROCH IRAM';
  143. MESS ' VALEUR CALCULéE AVEC KP : ' (extr XCOQ3 1) (extr XCOQ3 2);
  144. MESS ' ERREUR (%) : ' (extr ERR3 1) (extr ERR3 2) ;
  145. MESS ' VALEUR CALCULéE SANS KP : ' (extr XCOQ3 3) (extr XCOQ3 4);
  146. MESS ' VARIATION (THéORIQUE : 4/3) : '
  147. ((extr XCOQ3 3)/(extr XCOQ3 1)) ((extr XCOQ3 4)/(extr XCOQ3 2)) ;
  148. MESS '--------------------------------------------------------------';
  149. MESS '--- COQ3 : PROCH IRAM';
  150. MESS ' VALEUR CALCULéE AVEC KP : ' (extr XCOQ4 1) (extr XCOQ4 2);
  151. MESS ' ERREUR (%) : ' (extr ERR3 1) (extr ERR3 2) ;
  152. MESS ' VALEUR CALCULéE SANS KP : ' (extr XCOQ4 3) (extr XCOQ4 4);
  153. MESS ' VARIATION (THéORIQUE : 4/3) : '
  154. ((extr XCOQ4 3)/(extr XCOQ4 1)) ((extr XCOQ4 4)/(extr XCOQ4 2)) ;
  155. MESS '--------------------------------------------------------------';
  156.  
  157. *------------------ CODE DE BON FONCTIONNEMENT -------------------------
  158. SI (((MAXI ERR3 'ABS') < 5) ET ((MAXI ERR4 'ABS') < 5));
  159. ERRE 0;
  160. SINON;
  161. ERRE 5;
  162. FINSI;
  163.  
  164. FIN;
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  

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