Télécharger la plaquette Cast3M


fr en

Télécharger flurevi.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : flurevi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTION ECHO 0;                                                                  
  5. GRAPH='N';
  6. SAUT PAGE ;                                                                     
  7.  
  8. *                                                                               
  9. *******************************************************       
  10. *                                                                               
  11. *  TEST DU MODELE DE FLUAGE DE N.REVIRON
  12. *                                                                               
  13. *******************************************************       
  14. *
  15. a1 = 0.48;
  16. b1 = 25.e6;
  17. slmax = 1.;
  18.  
  19. * parametre de fluage propre
  20. * fluage propre
  21. k1 = 776.7e9;
  22. k2 = 472.e9;
  23. k3 = 98.1e9;
  24. eta1 = 2.44e12*1.e9;
  25.  
  26. opti 'DIME' 2 'MODE' axis 'ELEM' qua8;
  27.  
  28. *
  29. * GENERATION DU MAILLAGE
  30. *
  31.  
  32. h = 5.e-3;
  33. ri = 0.;
  34. re = ri + 8.e-2;
  35.  
  36. p1= ri 0.;
  37. p2= re 0.;
  38. p3= re h;
  39. p4= ri h;
  40.  
  41. dini1 = 1.003;
  42. dfin1 = 0.997;
  43.  
  44. nelem1 = 1;
  45. nelem2 = -30;
  46.  
  47. l1 = droi nelem2 p1 p2 'DINI' dini1 'DFIN' dfin1;
  48. l2 = droi nelem1 p2 P3 ;
  49. l3 = droi nelem2 p3 p4 'DINI' dfin1 'DFIN' dini1;
  50. l4 = droi nelem1 p4 p1 ;
  51.  
  52. surf1 = DALLER l1 l2 l3 l4;
  53.  
  54. SI(NEG GRAPH 'N');
  55.  trac surf1 ;
  56. FINSI ;
  57.  
  58. *
  59. * DISCRETISATION TEOMPORELLE
  60. *
  61. * 1h 1j 30j 700j
  62. * nombre de pas de temps entre les temps indiqués
  63. *
  64.  
  65. nb1 =  2 ;
  66. nb2 =  20 ;
  67. nb3 =  20 ;
  68. nb4 =  30 ;
  69.  
  70. nbt1 = nb1 + nb2 + nb3 + nb4;
  71.  
  72. temps1 = 0.1*3600.;
  73. temps2 = 2.99*24.*3600.;
  74. temps3 = 3.*24.*3600.;
  75. temps4 = 4.*24.*3600.;
  76.  
  77.  
  78. ltemps1 = PROG 0. 
  79.  pas (temps1/nb1) temps1
  80.  pas ((temps2-temps1)/nb2) temps2
  81.  pas ((temps3-temps2)/nb3) temps3
  82.  pas ((temps4-temps3)/nb4) temps4;
  83.  
  84. * ---------------------------
  85. * CALCUL MECANIQUE
  86. * ---------------------------
  87.  
  88. * Modele 
  89. mod1 = MODE surf1 'MECANIQUE' 'ELASTIQUE' 'VISCOPLASTIQUE' 'FLUTRA';
  90.  
  91.  
  92. * Module d Young
  93. eyou1 = 1.e20;
  94. *33.7e9;
  95. * Coefficient de Poisson
  96. nu1 = 0.248;
  97. * Seuil de fissuration en traction
  98. kt01 = 3.5e6/eyou1 ;
  99.  
  100. atram1 = -0.2;
  101. btram1 = 5.e3;
  102. acomm1 = 14.4;
  103. bcomm1 = 920.;
  104.  
  105. jour1 = 24. * 3600.;
  106.  
  107. mat1 = MATE mod1 'YOUN' eyou1 'NU  ' nu1 'KTR0' kt01
  108.  'ALPHA' 1.22e-5
  109.  'ATRA' atram1 'BTRA' btram1 'ACOM' acomm1 'BCOM' bcomm1 
  110.  'EKV1' k1 'TKV1' (0.1*jour1)  'EKV2' k2 'TKV2' (1.*jour1)
  111.  'EKV3' k3 'TKV3' (10.*jour1)  'EKV4' 1.e30  'TKV4' (100.*jour1)
  112.  'EKV5' 1.e30  'TKV5' (1000.*jour1)  'EKV6' 1.e30 'TKV6' (10.*jour1)
  113.  'ETA1' eta1 'NUPO' nu1 'BETA' 0. 'XFDE' 3.e-10
  114.  'VGEA' a1 'VGEB' b1 'SLMA' slmax;
  115.  
  116. * conditions aux limites mecaniques
  117.  
  118. rig1 = BLOQ 'UR' l4;
  119. rig2 = BLOQ 'UZ' l1;
  120.  
  121. * Evolution du chargement de traction
  122.  
  123. vcont1 = -1.e6;
  124. cont1 = PRES 'MASS' mod1 l3 vcont1;
  125. pre1 = PROG (nb1+1)*1. (nb2+nb3)*1. nb4*1.e-6;
  126. ev1 = EVOL 'MANU' 'Temps' ltemps1 'Pression' pre1;
  127.  
  128. SI(NEG GRAPH 'N');
  129.  DESS EV1 ;
  130. FINSI ;
  131.  
  132. cham1 = CHAR 'MECA' cont1 ev1;
  133.  
  134. * ------------------------
  135. * calcul pasapas
  136. * ------------------------
  137.  
  138. tab1=table;
  139.  
  140. tab1.caracteristiques=  mat1;
  141. tab1.chargement=cham1;
  142. tab1.modele= mod1;
  143. tab1.temps_calcules=ltemps1;
  144.  
  145. tab1.blocages_mecaniques=rig1 ET rig2;
  146. tab1.mova='D';
  147. *  pas de controle de positivite
  148. tab1.'STABILITE' = faux;
  149. tab1.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  150. pasapas tab1;
  151.  
  152.  
  153. ind1='INDEX' (TAB1.'CONTRAINTES') ;
  154. nfin1 = ('DIME' ind1) - 1;
  155.  
  156. p5= p3;
  157.  
  158. n3 = 0;
  159. dep1 = PROG;
  160. tmps = PROG;
  161.  
  162. REPETER BOUCLE (nfin1 + 1);
  163.  
  164. tmps1 = tab1.temps.n3;
  165. tmps1 = tmps1 / (24.*3600).;
  166. tmps = tmps ET (PROG tmps1);
  167.  
  168. val1 = EXTR tab1.deplacements.n3 'UZ' p5;
  169. val1 = -1. * val1 * 1.e12 / (vcont1 * (h));
  170. dep1 = dep1 ET (PROG val1);
  171.  
  172.  
  173. n3 = n3 + 1;
  174.  
  175. FIN BOUCLE;
  176.  
  177. evdep1 = EVOL 'ROUG' 'MANU' 'temps jours' tmps 'Deplact 1e6' dep1;
  178.  
  179. SI(NEG GRAPH 'N');
  180.     DESS evdep1;
  181. FINSI ;
  182. *
  183. TEMPS ;
  184. *             CODE FONCTIONNEMENT
  185.  
  186. DEPREF = PROG   0. 1.65860E-02  4.93848E-02   1.4608       2.0792
  187.        2.4920       2.8345       3.1404       3.4200       3.6778
  188.        3.9168       4.1396       4.3480       4.5439       4.7287
  189.        4.9038       5.0702       5.2290       5.3809       5.5267
  190.        5.6670       5.8024       5.9333       5.9337       5.9341 ;
  191. DEPREF = DEPREF ET ( PROG
  192.        5.9346       5.9350       5.9354       5.9359       5.9363
  193.        5.9367       5.9372       5.9376       5.9380       5.9384
  194.        5.9389       5.9393       5.9397       5.9402       5.9406 );
  195. DEPREF = DEPREF ET ( PROG
  196.        5.9410       5.9415       5.9419       5.7260       5.3419
  197.        5.0478       4.8190       4.6373       4.4899       4.3676
  198.        4.2638       4.1736       4.0937       4.0216       3.9554 );
  199. DEPREF = DEPREF ET ( PROG
  200.        3.8940       3.8363       3.7816       3.7296       3.6797
  201.        3.6317       3.5855       3.5407       3.4974       3.4554
  202.        3.4147       3.3751       3.3366       3.2991       3.2627
  203.        3.2273       3.1928       3.1592 ) ;
  204.  
  205. VMOY = 0. ;
  206. NVAL = DIME DEPREF ;
  207. REPE LAB NVAL  ;
  208.   VMOY = VMOY + ( EXTR DEPREF &LAB) ;
  209. FIN LAB;
  210. VMOY = VMOY / NVAL ;
  211. PRECIS = 2.E-5 ;
  212.  
  213. RESI1=(MAXI ABS (DEP1 - DEPREF))/ VMOY ;
  214. MESS '    ';
  215. MESS '    ';
  216. MESS '    ';
  217. mess ' RESIDU ' RESI1  '  precision ' PRECIS ;
  218. MESS '    ';
  219.  
  220. SI(RESI1 >  PRECIS);
  221. EVAL = EVOL 'ROUG' 'MANU' 'REFERENCE' DEPREF 'CALCUL' DEP1;
  222. LIST EVAL ;
  223.     ERRE  5;
  224. FINSI;
  225. FIN;
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  

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