Télécharger b52c.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : b52c.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. 'OPTI' 'ECHO' 1 'DIME' 2 'ELEM' 'QUA8' 'MODE' 'AXIS' ;
  5. GRAPH = 'N' ;
  6.  
  7. **
  8. ** D3
  9. **
  10. ** D4 S1 D2
  11. **
  12. ** P1 D1 P2
  13. *
  14. * Tube
  15. X1 = 10. ;
  16. X2 = 11. ;
  17. P1 = X1 0. ;
  18. P2 = X2 0. ;
  19. D1 = 'DROI' 10 P1 P2 ;
  20. S1 = 'TRAN' D1 5 (0. 5.) ;
  21. D2 = 'COTE' 2 S1 ;
  22. D3 = 'COTE' 3 S1 ;
  23. D4 = 'COTE' 4 S1 ;
  24. PP = D2 'POIN' 'FINAL' ;
  25. ** THERMIQUE
  26. * Température imposée à l'extérieur et l'intérieur du cylindre
  27. T2 = 'TABLE' ;
  28. T4 = 'TABLE' ;
  29. T2.1 = 1400. ;
  30. T4.1 = 1400. ;
  31. T2.2 = 1400. ;
  32. T4.2 = 1405. ;
  33. T2.3 = 1400. ;
  34. T4.3 = 1410. ;
  35. T2.4 = 1400. ;
  36. T4.4 = 1415. ;
  37. T2.5 = 1400. ;
  38. T4.5 = 1420. ;
  39. T2.6 = 1400. ;
  40. T4.6 = 1425. ;
  41. T2.7 = 1400. ;
  42. T4.7 = 1430. ;
  43. * Facteur multiplicatif des temps
  44. MULT1 = 1.E2 ;
  45. * modèle
  46. MODS1 = 'MODE' S1 'THERMIQUE' 'ISOTROPE' ;
  47. **
  48. * matériau
  49. MATS1 = 'MATE' MODS1 'K' 5. ;
  50. * conductivité
  51. CONDS1 = 'COND' MODS1 MATS1 ;
  52. * Température imposée
  53. * sur la surface externe
  54. BLOQ2 = 'BLOQ' D2 'T' ;
  55. * sur la droite interne
  56. BLOQ4 = 'BLOQ' D4 'T' ;
  57. CHTER = 'TABLE' ;
  58. NT = 7 ;
  59. 'REPE' BOULOU NT ;
  60. LEA = (T4.&BOULOU - T2.&BOULOU) / (X1-X2) ;
  61. LEB = T2.&BOULOU - (LEA * X2) ;
  62. CHTER.&BOULOU = 'NOMC' 'T' (LEA * ('COOR' 1 S1) + LEB) ;
  63. 'FIN' BOULOU ;
  64. * MECANIQUE
  65. MODME = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'CERAMIQUE' ;
  66. * On suppose les coefficient de fluage constant pour l'instant
  67. QFLUGU = -5.3/8.61735E-5 ;
  68. NTFLUGU = 1. ;
  69. NFLUGU = 3.9 ;
  70. AFLUGU = ('EXP' 28.8)/(1.E6**3.9);
  71. * On met le terme dépendant de la température dans une EVOL
  72. YOYO = 1.E9 ;
  73. YOYO = 5.E11 ;
  74. YOYO = 1.E10 ;
  75. NUNU = 0.3 ;
  76. DILA = 9.E-26 ;
  77. RORO = 9000. ;
  78. SIMAX = 1.E6 ;
  79. QFLUGU = -5.3/8.61735E-5 ;
  80. NTFLUGU = 1. ;
  81. NFLUGU = 3.9 ;
  82. AFLUGU = ('EXP' 28.8)/(1.E6**3.9);
  83. ***********************************
  84. CHTER3 = 0. ;
  85. 'REPE' BOUOU 100 ;
  86. CHTER3 = CHTER3 + 20.;
  87. AVFLUK = AFLUGU * ('EXP' (QFLUGU * (CHTER3**-1.))) ;
  88. 'SI' ( CHTER3 '<' 1400.) ;
  89. AVFLUK = 1.E-40 ;
  90. 'FINS' ;
  91. 'SI' ('EGA' &BOUOU 1) ;
  92. PROT = 'PROG' CHTER3 ;
  93. PROAV = 'PROG' AVFLUK ;
  94. 'SINO' ;
  95. PROT = 'INSE' PROT &BOUOU CHTER3 ;
  96. PROAV = 'INSE' PROAV &BOUOU AVFLUK ;
  97. 'FINS' ;
  98. 'FIN' BOUOU ;
  99. AVFLUKE = 'EVOL' 'MANU' 'T' PROT 'AF1' PROAV;
  100. DILA = 9.E-06 ;
  101. MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU 'RHO' RORO 'ALPH' DILA 'TREF' 0. 'TALP' 0. 'SMAX' SIMAX 'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU 'TTRA' 0. 'ENDG' 100.01;
  102. MAT2 = 'TAIL' MODME ;
  103. * BLOCAGE DE D1
  104. BLOCO = 'BLOQ' 'UZ' D1 ;
  105. CHAMPI = 'DEPI' BLOCO 0. ;
  106. * BLOCAGE DES POINTS DE D3
  107. BLOCO3 = 'BLOQ' 'UZ' D3 ;
  108. * Table des temps
  109. TATEMPS = 'TABLE' ;
  110. TATEMPS.0 = 0. ;
  111. TATEMPS.1 = 1.*MULT1 ;
  112. TATEMPS.2 = 2.*MULT1 ;
  113. TATEMPS.3 = 3.*MULT1 ;
  114. TATEMPS.4 = 4.*MULT1 ;
  115. TATEMPS.5 = 5.*MULT1 ;
  116. TATEMPS.6 = 6.*MULT1 ;
  117. TATEMPS.7 = 7.*MULT1 ;
  118. TATEMPS.8 = 8.*MULT1 ;
  119. TATEMPS.9 = 9.*MULT1 ;
  120. TATEMPS.10 = 10.*MULT1 ;
  121. TATEMPS.11 = 11.*MULT1 ;
  122. TATEMPS.12 = 12.*MULT1 ;
  123. TATEMPS.13 = 13.*MULT1 ;
  124. TATEMPS.14 = 14.*MULT1 ;
  125. TATEMPS.15 = 15.*MULT1 ;
  126. TATEMPS.16 = 16.*MULT1 ;
  127. TATEMPS.17 = 17.*MULT1 ;
  128. TATEMPS.18 = 18.*MULT1 ;
  129. TATEMPS.19 = 19.*MULT1 ;
  130. * Table d'évolution de la température
  131. TABTER = 'TABLE' ;
  132. TABTER.0 = 'MANU' 'CHPO' S1 1 'T' 0. ;
  133. TABTER.1 = CHTER.1;
  134. TABTER.2 = CHTER.2;
  135. TABTER.3 = CHTER.3;
  136. TABTER.4 = CHTER.4;
  137. TABTER.5 = CHTER.5;
  138. TABTER.6 = CHTER.6;
  139. TABTER.7 = CHTER.7;
  140. TABTER.8 = CHTER.7;
  141. TABTER.9 = CHTER.7;
  142. TABTER.10 = CHTER.7;
  143. TABTER.11 = CHTER.7;
  144. TABTER.12 = CHTER.7;
  145. TABTER.13 = CHTER.7;
  146. TABTER.14 = CHTER.7;
  147. TABTER.15 = CHTER.7;
  148. TABTER.16 = CHTER.7;
  149. TABTER.17 = CHTER.7;
  150. TABTER.18 = CHTER.7;
  151. TABTER.19 = CHTER.7;
  152. * Table du chargement mécanique
  153. FORTOT = CHAMPI ;
  154. CHARMEC = 'TABLE' ;
  155. CHARMEC.0 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.0)));
  156. CHARMEC.1 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.1)));
  157. CHARMEC.2 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.2)));
  158. CHARMEC.3 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.3)));
  159. CHARMEC.4 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.4)));
  160. CHARMEC.5 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.5)));
  161. CHARMEC.6 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.6)));
  162. CHARMEC.7 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.7)));
  163. DILAT = (DILA * 5.*('MAXI' TABTER.7));
  164. CHARMEC.8 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+2.E-4) );
  165. CHARMEC.9 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+4.E-4));
  166. CHARMEC.10 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+6.E-4));
  167. CHARMEC.11 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+3.E-4));
  168. CHARMEC.12 = FORTOT 'ET' ('DEPI' BLOCO3 DILAT);
  169. CHARMEC.13 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-3.E-4));
  170. CHARMEC.14 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-6.E-4));
  171. CHARMEC.15 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-9.E-4));
  172. CHARMEC.16 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-12.E-4));
  173. CHARMEC.17 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-15.E-4));
  174. CHARMEC.18 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-18.E-4));
  175. CHARMEC.19 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-21.E-4));
  176. * Chargement temps force
  177. CHARFOR = 'CHAR' 'MECA' TATEMPS CHARMEC ;
  178. * Chargement temps température
  179. CHARGT = 'CHAR' 'T' TATEMPS TABTER ;
  180. TABA = 'TABLE' ;
  181. TABA.'CHARGEMENT' = CHARGT 'ET' CHARFOR ;
  182. LISTEMPS = 'PROG' 0. 1. 'PAS' 1. 19.;
  183. LISTEMPS = LISTEMPS*MULT1 ;
  184. TABA.'MODELE' = MODME ;
  185. TABA.'CARACTERISTIQUES' = MATME 'ET' MAT2;
  186. TABA.'TEMPS_CALCULES' = LISTEMPS ;
  187. TABA.'TEMPS_SAUVES' = LISTEMPS ;
  188. TABA.'BLOCAGES_MECANIQUES' = BLOCO 'ET' BLOCO3 ;
  189. TABA.'MAXITERATION' = 500 ;
  190. TABA.'PRECISION' = 4.E-4 ;
  191. TABA.'HYPOTHESE_DEFORMATION' = 'LINEAIRE';
  192. ***********************************
  193. PASAPAS TABA ;
  194. I = 'DIME' TABA.'TEMPS' ;
  195. * On essaie d'extraire les fissures
  196. * on extrait les variables internes
  197. FIFISS = TABA.'VARIABLES_INTERNES'.(I-1) ;
  198. VECI = 'VECT' FIFISS MODME .01 ;
  199. *
  200. 'REPE' BOULOU 3 ;
  201. 'SI' ('EGA' &BOULOU 1) ;
  202. 'SI' ('NEG' GRAPH 'N') ;
  203. 'TRAC' S1 VECI ;
  204. 'FINS' ;
  205. MA = 'ELEM' 30 S1 ;
  206. 'SINO' ;
  207. 'SI' ('EGA' &BOULOU 2) ;
  208. MA = 'ELEM' 21 S1 ;
  209. 'SINO' ;
  210. MA = 'ELEM' 23 S1 ;
  211. 'FINS' ;
  212. 'FINS' ;
  213. 'REPE' BOILA I ;
  214. * Contraintes
  215. SI1 = 'REDU' TABA.'CONTRAINTES'.(&BOILA-1) MA ;
  216. SILO2 = 'EXTR' SI1 'SMZZ' 1 1 2 ;
  217. 'MESS' 'SILO2 = ' SILO2 ;
  218. * Déformations
  219. DE1 = TABA.'DEPLACEMENTS'.(&BOILA-1) ;
  220. EEP = 'REDU' ('EPSI' MODME DE1 ) MA ;
  221. * On retranche la déformation thermique
  222. ETH = 'EPTH' MODME MATME TABTER.(&BOILA-1) ;
  223. EEPTH = 'REDU' ETH MA ;
  224. EEVOU = EEP - EEPTH ;
  225. EZDE2 = 'EXTR' EEVOU 'EPZZ' 1 1 2 ;
  226. TIMER = TABA.'TEMPS'.(&BOILA-1) ;
  227. 'MESS' 'EZDE2 ='EZDE2 ;
  228. 'SI' ('EGA' &BOILA 1) ;
  229. LOOX = 'PROG' EZDE2 ;
  230. LOOY = 'PROG' SILO2 ;
  231. LOOT = 'PROG' TIMER ;
  232. 'SINO' ;
  233. LOOX = 'INSE' LOOX &BOILA EZDE2;
  234. LOOY = 'INSE' LOOY &BOILA SILO2 ;
  235. LOOT = 'INSE' LOOT &BOILA TIMER ;
  236. 'FINS' ;
  237. 'FIN' BOILA ;
  238. 'SI' ('NEG' GRAPH 'N') ;
  239. TSIG = ('EVOL' 'MANU' 'TEMPS' LOOT 'SIGZZ' LOOY) ;
  240. 'DESS' TSIG ;
  241. TEPS = ('EVOL' 'MANU' 'TEMPS' LOOT 'EPSZZ' LOOX) ;
  242. 'DESS' TEPS ;
  243. TOLO = ('EVOL' 'MANU' 'EPSZZ' LOOX 'SIGZZ' LOOY) ;
  244. 'DESS' TOLO ;
  245. 'FINS' ;
  246. 'FIN' BOULOU ;
  247. 'MESS' '***************************************************' ;
  248. 'FIN' ;
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.  

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