Télécharger b52c.dgibi

Retour à la liste

Numérotation des lignes :

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

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