Télécharger la plaquette Cast3M


fr en

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