Télécharger la plaquette Cast3M


fr en

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 'TALP' 0. 'TREF' 0.;
  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 'TALP' 0. 'TREF' 0. '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 'TALP' 0. 'TREF' 0. '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 'TALP' 0. 'TREF' 0. '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 'TALP' 0. 'TREF' 0. '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.  
  313.  
  314.  

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