Télécharger te35.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : te35.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' 'MODE' 'PLAN' 'CONT' ;
  5. **
  6. **
  7. **
  8. ** D4 S1 D2
  9. **
  10. ** P1 D1 P2
  11. *
  12. GRAPH = 'N' ;
  13. *
  14. * OPTIONS DE CALCUL
  15. * CHOIMEC = 1 OTTOSEN
  16. * CHOIMEC = 2 CERAMIQUE
  17. * CHOIMEC = 3 FLUAGE
  18. CHOIMEC = 2;
  19. * DANS LE CAS OU ON A LE MODELE CERAMIQUE ON PEUT INACTIVER OTTOSEN
  20. * INACTO = VRAI
  21. * OU INACTIVER LE FLUAGE
  22. * INACTF = VRAI
  23. 'SI' ('EGA' CHOIMEC 2) ;
  24. INACTO = FAUX ;
  25. INACTF = FAUX ;
  26. 'FINS' ;
  27. CHOICAS = 1;
  28. * CHOICAS = 1
  29. * ON CASSE EN DIMINUANT LA TEMPERATURE AU DELA DU TEMPS 5E-5
  30. * On passe de 1500 à 1400
  31. * CHOICAS = 2
  32. * ON CASSE EN AUGMENTATION LA VITESSE DE DEFORMATION /2
  33. * CHOICAS = 3 On reste à T = cte
  34. * AU DELA DU TEMPS 5E-5
  35. * facteur de multiplication des temps qui vont de 1 à 9 MULT1
  36. MULT1 = 1.E-5 ;
  37. *
  38. X1 = 10. ;
  39. X2 = 11. ;
  40. P1 = X1 0. ;
  41. P2 = X2 0. ;
  42. D1 = 'DROI' 10 P1 P2 ;
  43. S1 = 'TRAN' D1 10 (0. 5.) ;
  44. 'SI' ('NEG' GRAPH 'N') ;
  45. 'TRAC' S1 ;
  46. 'FINS' ;
  47. D2 = 'COTE' 2 S1 ;
  48. D3 = 'COTE' 3 S1 ;
  49. D4 = 'COTE' 4 S1 ;
  50. PP = D2 'POIN' 'FINAL' ;
  51. ** THERMIQUE
  52. * Température imposée à l'extérieur et l'intérieur de la
  53. * meule
  54. TCHOIX = 1500. ;
  55. CHTER = 'MANU' 'CHPO' S1 1 'T' TCHOIX ;
  56. 'SI' ('EGA' CHOICAS 1 );
  57. TCHOIX2 = 1400. ;
  58. CHTER2 = 'MANU' 'CHPO' S1 1 'T' TCHOIX2 ;
  59. 'SINO' ;
  60. CHTER2 = CHTER ;
  61. 'FINS' ;
  62. **********************************************
  63. ************** MODELE ET MATERIAUX MECANIQUE
  64. **********************************************
  65. *
  66. YOYO = 5.E11 ;
  67. NUNU = 0.3 ;
  68. DILA = 9.E-26 ;
  69. RORO = 9000. ;
  70. * OTTOSEN
  71. 'SI' ('EGA' CHOIMEC 1) ;
  72. MODME = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'OTTOSEN' ;
  73. MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU
  74. 'RHO' RORO 'ALPH' DILA;
  75. MAT2 = 'TAIL' MODME ;
  76. 'FINS' ;
  77. 'SI' ( 'OU' ('EGA' CHOIMEC 2) ('EGA' CHOIMEC 3)) ;
  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. CHTER3 = 0. ;
  84. 'REPE' BOUOU 100 ;
  85. CHTER3 = CHTER3 + 20.;
  86. AVFLUK = AFLUGU * ('EXP' (QFLUGU * (CHTER3**-1.))) ;
  87. 'SI' ('<' CHTER3 1400.) ;
  88. AVFLUK = 1.E-40 ;
  89. 'FINS' ;
  90. 'SI' ('EGA' &BOUOU 1) ;
  91. PROT = 'PROG' CHTER3 ;
  92. PROAV = 'PROG' AVFLUK ;
  93. 'SINO' ;
  94. PROT = 'INSE' PROT &BOUOU CHTER3 ;
  95. PROAV = 'INSE' PROAV &BOUOU AVFLUK ;
  96. 'FINS' ;
  97. 'FIN' BOUOU ;
  98. AVFLUKE = 'EVOL' 'MANU' 'T' PROT 'AF1' PROAV;
  99. * CERAMIQUE
  100. 'SI' ('EGA' CHOIMEC 2) ;
  101. MODME = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'CERAMIQUE' ;
  102. 'SI' ('OU' INACTO INACTF) ;
  103. 'SI' INACTO ;
  104. MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU
  105. 'RHO' RORO 'ALPH' DILA 'SMAX' SIMAX
  106. 'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU
  107. 'LTR' 1.E20 'GS' 1.E20 'TTRA' 0. 'ENDG' 100.01;
  108. 'FINS' ;
  109. 'SI' INACTF ;
  110. AVFLUK = 1.E-40 ;
  111. MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU
  112. 'RHO' RORO 'ALPH' DILA 'SMAX' SIMAX
  113. 'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU
  114. 'TTRA' 0. 'ENDG' 100.01;
  115. 'FINS' ;
  116. 'SINO' ;
  117. MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU
  118. 'RHO' RORO 'ALPH' DILA 'SMAX' SIMAX
  119. 'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU
  120. 'TTRA' 0. 'ENDG' 100.01;
  121. 'FINS' ;
  122. MAT2 = 'TAIL' MODME ;
  123. 'FINS' ;
  124. 'SI' ('EGA' CHOIMEC 3) ;
  125. MODME = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'NORTON' ;
  126. MATME = 'MATE' MODME 'YOUN' YOYO 'NU'NUNU
  127. 'RHO' RORO 'ALPH' DILA 'SMAX' SIMAX
  128. 'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU ;
  129. 'FINS' ;
  130. 'FINS' ;
  131. ***************************************************
  132. ******** BLOQUAGES MECA
  133. *********************************************
  134. * BLOCAGE DE D1
  135. BLOCO = 'BLOQ' 'UY' D1 ;
  136. * BLOCAGE DES POINTS DE D3
  137. BLOCO3 = 'BLOQ' 'UY' D3 ;
  138. BLOP = 'BLOQ' 'UX' 'UY' P1 ;
  139. *
  140. 'SI' ('EGA' CHOICAS 2) ;
  141. FACTI = 10. ;
  142. 'SINO' ;
  143. FACTI = 1.;
  144. 'FINS' ;
  145. DITER = DILA * 1500.*5. ;
  146. 'MESS' 'DITER = ' DITER ;
  147. CHAMPI = 'DEPI' BLOCO 0. ;
  148. CHAMPI31 = 'DEPI' BLOCO3 DITER ;
  149. CHAMPI32 = 'DEPI' BLOCO3 (DITER + 0.5E-3) ;
  150. CHAMPI33 = 'DEPI' BLOCO3 (DITER + 1.E-3) ;
  151. CHAMPI34 = 'DEPI' BLOCO3 (DITER + 1.5E-3) ;
  152. CHAMPI35 = 'DEPI' BLOCO3 (DITER + 2.E-3) ;
  153. *
  154. HILO = 2.E-3 + (0.5E-3*FACTI) ;
  155. CHAMPI36 = 'DEPI' BLOCO3 (DITER + HILO) ;
  156. HILO = HILO + (0.5E-3*FACTI) ;
  157. CHAMPI37 = 'DEPI' BLOCO3 (DITER + HILO) ;
  158. HILO = HILO + (0.5E-3*FACTI) ;
  159. CHAMPI38 = 'DEPI' BLOCO3 (DITER + HILO) ;
  160. HILO = HILO + (0.5E-3*FACTI) ;
  161. CHAMPI39 = 'DEPI' BLOCO3 (DITER + HILO) ;
  162. * Table des temps
  163. TATEMPS = 'TABLE' ;
  164. TATEMPS.0 = 0. ;
  165. TATEMPS.1 = 1.*MULT1 ;
  166. TATEMPS.2 = 2.*MULT1 ;
  167. TATEMPS.3 = 3.*MULT1 ;
  168. TATEMPS.4 = 4.*MULT1 ;
  169. TATEMPS.5 = 5.*MULT1 ;
  170. TATEMPS.6 = 6.*MULT1 ;
  171. TATEMPS.7 = 7.*MULT1 ;
  172. TATEMPS.8 = 8.*MULT1 ;
  173. TATEMPS.9 = 9.*MULT1 ;
  174. TATEMPS.10 = 10.*MULT1 ;
  175. TATEMPS.11 = 11.*MULT1 ;
  176. TATEMPS.12 = 12.*MULT1 ;
  177. * Table d'évolution de la température
  178. TABTER = 'TABLE' ;
  179. TABTER.0 = 'MANU' 'CHPO' S1 1 'T' 0. ;
  180. TABTER.1 = CHTER;
  181. TABTER.2 = CHTER;
  182. TABTER.3 = CHTER;
  183. TABTER.4 = CHTER;
  184. TABTER.5 = CHTER;
  185. TABTER.6 = CHTER2;
  186. TABTER.7 = CHTER2;
  187. TABTER.8 = CHTER2;
  188. TABTER.9 = CHTER2;
  189. TABTER.9 = CHTER2;
  190. TABTER.10 = CHTER2;
  191. TABTER.11 = CHTER2;
  192. TABTER.12 = CHTER2;
  193. * Table du chargement mécanique
  194. FORTOT = CHAMPI;
  195. CHARMEC = 'TABLE' ;
  196. CHARMEC.0 = FORTOT ;
  197. CHARMEC.1 = FORTOT 'ET' CHAMPI31;
  198. CHARMEC.2 = FORTOT 'ET' CHAMPI32;
  199. CHARMEC.3 = FORTOT 'ET' CHAMPI33;
  200. CHARMEC.4 = FORTOT 'ET' CHAMPI34;
  201. CHARMEC.5 = FORTOT 'ET' CHAMPI35;
  202. CHARMEC.6 = FORTOT 'ET' CHAMPI36;
  203. CHARMEC.7 = FORTOT 'ET' CHAMPI37;
  204. CHARMEC.8 = FORTOT 'ET' CHAMPI38;
  205. CHARMEC.9 = FORTOT 'ET' CHAMPI39;
  206. CHARMEC.10 = FORTOT 'ET' CHAMPI33;
  207. CHARMEC.11 = FORTOT 'ET' CHAMPI33;
  208. CHARMEC.12 = FORTOT 'ET' CHAMPI33;
  209. * Chargement temps force
  210. CHARFOR = 'CHAR' 'MECA' TATEMPS CHARMEC ;
  211. * Chargement temps température
  212. CHARGT = 'CHAR' 'T' TATEMPS TABTER ;
  213. TABA = 'TABLE' ;
  214. TABA.'CHARGEMENT' = CHARGT 'ET' CHARFOR ;
  215. LISTEMPS = 'PROG' 0. 1. 'PAS' 0.2 5. 'PAS' 0.01 5.2 'PAS' 0.2 9.;
  216. LISTEMPS = LISTEMPS * MULT1 ;
  217. TABA.'MODELE' = MODME ;
  218. 'SI' ('OU' ('EGA' CHOIMEC 1) ('EGA' CHOIMEC 2)) ;
  219. TABA.'CARACTERISTIQUES' = MATME 'ET' MAT2;
  220. 'SINO' ;
  221. TABA.'CARACTERISTIQUES' = MATME ;
  222. 'FINS' ;
  223. TABA.'TEMPS_CALCULES' = LISTEMPS ;
  224. TABA.'TEMPS_SAUVES' = LISTEMPS ;
  225. TABA.'PRECISON' = 1.E-4 ;
  226. TABA.'BLOCAGES_MECANIQUES' = BLOCO 'ET' BLOCO3 'ET' BLOP;
  227. TABA.'MAXITERATION' = 500 ;
  228. TABA.'PROCESSEURS' = 'AUTOMATIQUE';
  229. PASAPAS TABA ;
  230. I = 'DIME' TABA.'TEMPS_SAUVES';
  231. *
  232. 'SI' ('OU' ('EGA' CHOIMEC 1) ('EGA' CHOIMEC 2)) ;
  233. * On d'extrait les fissures
  234. * on extrait les variables internes
  235. FIFISS = TABA.'VARIABLES_INTERNES'.(I-1) ;
  236. VECI = 'VECT' FIFISS MODME .01 ;
  237. *
  238. 'SI' ('NEG' GRAPH 'N') ;
  239. 'TRAC' S1 VECI ;
  240. 'FINS' ;
  241. 'FINS' ;
  242.  
  243. MA = 'ELEM' 61 S1 ;
  244. 'REPE' BOILA I ;
  245. * Contraintes
  246. SI1 = 'REDU' TABA.'CONTRAINTES'.(&BOILA-1) MA ;
  247. SILO2 = 'EXTR' SI1 'SMXX' 1 1 2 ;
  248. SILO3 = 'EXTR' SI1 'SMYY' 1 1 2 ;
  249. SILO4 = 'EXTR' SI1 'SMZZ' 1 1 2 ;
  250. SILO5 = 'EXTR' SI1 'SMXY' 1 1 2 ;
  251. * Déformations
  252. DE1 = TABA.'DEPLACEMENTS'.(&BOILA-1) ;
  253. EEP = 'REDU' ('EPSI' MODME DE1 ) MA ;
  254. * On retranche la déformation thermique
  255. TATA = 'TIRE' CHARGT (TABA.'TEMPS'.(&BOILA-1)) ;
  256. ETH = 'EPTH' MODME MATME TATA;
  257. EEPTH = 'REDU' ETH MA ;
  258. EEVOU = EEP - EEPTH ;
  259. EZDE2 = 'EXTR' EEVOU 'EPXX' 1 1 2 ;
  260. EZDE3 = 'EXTR' EEVOU 'EPYY' 1 1 2 ;
  261. EZDE4 = 'EXTR' EEVOU 'EPZZ' 1 1 2 ;
  262. EZDE5 = 'EXTR' EEVOU 'GAXY' 1 1 2 ;
  263. TIMER = TABA.'TEMPS'.(&BOILA-1) ;
  264. 'MESS' 'EZDE2 ='EZDE2 ;
  265. 'SI' ('EGA' &BOILA 1) ;
  266. LOOX2 = 'PROG' EZDE2 ;
  267. LOOX3 = 'PROG' EZDE3 ;
  268. LOOX4 = 'PROG' EZDE4 ;
  269. LOOX5 = 'PROG' EZDE5 ;
  270. LOOY2 = 'PROG' SILO2 ;
  271. LOOY3 = 'PROG' SILO3 ;
  272. LOOY4 = 'PROG' SILO4 ;
  273. LOOY5 = 'PROG' SILO5 ;
  274. LOOT = 'PROG' TIMER ;
  275. 'SINO' ;
  276. LOOX2 = 'INSE' LOOX2 &BOILA EZDE2 ;
  277. LOOX3 = 'INSE' LOOX3 &BOILA EZDE3 ;
  278. LOOX4 = 'INSE' LOOX4 &BOILA EZDE4 ;
  279. LOOX5 = 'INSE' LOOX5 &BOILA EZDE5 ;
  280. LOOY2 = 'INSE' LOOY2 &BOILA SILO2 ;
  281. LOOY3 = 'INSE' LOOY3 &BOILA SILO3 ;
  282. LOOY4 = 'INSE' LOOY4 &BOILA SILO4 ;
  283. LOOY5 = 'INSE' LOOY5 &BOILA SILO5 ;
  284. LOOT = 'INSE' LOOT &BOILA TIMER ;
  285. 'FINS' ;
  286. 'FIN' BOILA ;
  287. 'SI' ('NEG' GRAPH 'N') ;
  288. TSIG = 'EVOL' 'MANU' 'TEMPS' LOOT 'SIGMAXX' LOOY2 ;
  289. 'DESS' TSIG ;
  290. TSIG = 'EVOL' 'MANU' 'TEMPS' LOOT 'SIGMAYY' LOOY3 ;
  291. 'DESS' TSIG ;
  292. TSIG = 'EVOL' 'MANU' 'TEMPS' LOOT 'SIGMAZZ' LOOY4 ;
  293. 'DESS' TSIG ;
  294. TSIGES = 'EVOL' 'MANU' 'EPSILYY' LOOX3 'SIGMAYY' LOOY3 ;
  295. 'DESS' TSIGES ;
  296. TEPS = 'EVOL' 'MANU' 'TEMPS' LOOT 'EPSILXX' LOOX2 ;
  297. 'DESS' TEPS ;
  298. TEPS = 'EVOL' 'MANU' 'TEMPS' LOOT 'EPSILYY' LOOX3 ;
  299. 'DESS' TEPS ;
  300. TEPS = 'EVOL' 'MANU' 'TEMPS' LOOT 'EPSILZZ' LOOX4 ;
  301. 'DESS' TEPS ;
  302. 'FINS' ;
  303. 'MESS' '***************************************************' ;
  304. 'FIN' ;
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  

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