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

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