Télécharger flua1t.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : flua1t.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux;
  5. *
  6. *************************************************
  7. * pour calcul complet mettre complet à : vrai;
  8. *************************************************
  9. *------------------------------------------------------------------;
  10. * ;
  11. TITRE 'CREEP DATA FOR SS 304 LN (RCC-MR) USING 3-D FORMULATION' ;
  12. ****************************************************************
  13. * *
  14. * TEST DE FLUAGE SUIVANT LE MODELE RCC-MR POUR ACIERS 304-SS *
  15. * POUR UN CUBE SOUMIS A UNE PRESSION SUIVANT UNE FACE. * *
  16. * *
  17. * COMPARAISON AVEC LA FORMULE ANALYTIQUE DE FLUAGE *
  18. * *
  19. * TEST ELABORE PAR Mr MILLARD. *
  20. * TEST MODIFIE POUR EVOLUTION LE 7/07/92 *
  21. * *
  22. *------------------------------------------------------------------;
  23. OPTION ECHO 0 ;
  24. GRAPH = 'N';
  25. OPTION DIME 3 ELEM CUB8 ;
  26. DENS 1. ;
  27. OPTI TRAC BENS;
  28. *__________________________________________________________________;
  29. * ;
  30. * SKETCH OF GEOMETRY ;
  31. * ;
  32. *__________________________________________________________________;
  33. * ;
  34. * ;
  35. * .____.____. ;
  36. * /. /| ;
  37. * / . PZ / | ;
  38. * PY/__._.____/... ;
  39. * | . | / ;
  40. * |. |/ ;
  41. * P0.____.____/ ;
  42. * LOX PX ;
  43. * ;
  44. * ;
  45. * ;
  46. *__________________________________________________________________;
  47. * GEOMETRY DETAIL ;
  48. *__________________________________________________________________;
  49. P0 = 0 0 0 ;
  50. PX = 1 0 0 ;
  51. PY = 0 1 0 ;
  52. PZ = 0 0 1 ;
  53. L0X = P0 DROI 1 PX ;
  54. FAC1 = L0X TRAN 1 PY ;
  55. VOLM = FAC1 VOLU 1 TRAN PZ ;
  56. OEIL = 20 -30 10 ;
  57. FAC4 = FACE VOLM 2 ;
  58.  
  59. SI (NEG GRAPH 'N');
  60. TRACE OEIL VOLM ;
  61. FINSI;
  62. *__________________________________________________________________;
  63. * BOUNDARY CONDITIONS ;
  64. *__________________________________________________________________;
  65. FAC2 = VOLM POIN PLAN P0 PX PZ (0.1) ;
  66. FAC3 = VOLM POIN PLAN P0 PY PZ (0.1) ;
  67. CL1 = BLOQ UZ FAC1 ;
  68. CL2 = BLOQ UY FAC2 ;
  69. CL3 = BLOQ UX FAC3 ;
  70. CLY = CL1 ET CL2 ET CL3 ;
  71. *__________________________________________________________________;
  72. * MATERIAL DATA ;
  73. *__________________________________________________________________;
  74. MO = MODE VOLM 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'RCCMR_304' ;
  75.  
  76. MAT = MATE MO YOUN 2.00E5 NU 0.3 RHO 7.8E3 ALPH 20.0E-06 SMAX 50.0
  77. BF1 0.222E-6 BF2 0.036797 BF3 6.0
  78. RF1 0.4051E-4 RF2 0.031829 RF3 3.5
  79. AF1 7.544324 SF1 0.2266E-2 AF2 -0.257143E-2
  80. AF3 0.62202E-4 SF2 41.34 ;
  81.  
  82. *__________________________________________________________________;
  83. * EXTERNAL LOADING ;
  84. * ;
  85. *------------------------------------------------------------------;
  86. LOAD = PRES MASS MO -200.0 FAC4 ;
  87. * ;
  88. *__________________________________________________________________;
  89. * STIFFNESS COMPUTATION ;
  90. *__________________________________________________________________;
  91. RI1 = RIGI MO MAT ;
  92. RI2 = RI1 ET CLY ;
  93. DD = RESOL RI2 LOAD ;
  94. SS = SIGMA MAT MO DD ;
  95. SXX = EXTR SS SMZZ 1 1 1 ;
  96. *__________________________________________________________________;
  97. * INITIALISATION DE LA TABLE ;
  98. *__________________________________________________________________;
  99. LI1 = PROG 0.0 100.0 ;
  100. LI2 = PROG 1.0 1.0 ;
  101. EV = EVOL MANU T LI1 F(T) LI2 ;
  102. CHA = CHAR 'MECA' LOAD EV ;
  103. si complet;
  104. LIS = PROG 0.01 PAS 0.01 0.1 0.2 PAS 0.1 1.0 2.0 PAS 1.0 10.0 ;
  105. sinon;
  106. LIS = PROG 0.01 PAS 0.01 0.1 ;
  107. finsi;
  108. TAB = TABLE ;
  109. TAB . 'CONTRAINTES' = TABLE ;
  110. TAB . 'DEPLACEMENTS' = TABLE ;
  111. TAB . 'CONTRAINTES' . 0 = SS ;
  112. TMASAU=table;
  113. tab . 'MES_SAUVEGARDES'=TMASAU;
  114. TMASAU .'DEFTO'=VRAI;
  115. TMASAU .'DEFIN'=VRAI;
  116. TAB . 'DEPLACEMENTS' . 0 = DD ;
  117. TAB.'BLOCAGES_MECANIQUES' = CLY ;
  118. TAB.'MODELE' = MO ;
  119. TAB.'CARACTERISTIQUES' = MAT ;
  120. TAB.'TEMPS_CALCULES' = LIS ;
  121. TAB.'CHARGEMENT' = CHA ;
  122.  
  123. PASAPAS TAB ;
  124. *__________________________________________________________________;
  125. * BOUCLE DE CALCUL ANALYTIQUE ET D'EXTRACTION
  126. * DES RESULTATS
  127. *__________________________________________________________________;
  128. TAB2 = TAB.CONTRAINTES ;
  129. TAB3 = TAB.DEFORMATIONS_INELASTIQUES ;
  130.  
  131. NPAS = 'DIME' LIS ;
  132. I = 0 ; ;
  133. LTT = PROG 0.000 ;
  134. LSXX = PROG SXX ;
  135. LEXX = PROG 000.0 ;
  136. ldefeq = prog 0.;
  137. S = 200.0;
  138. LEF= prog 0.00;
  139.  
  140. 'REPETER' AFFICHE NPAS ;
  141.  
  142.  
  143. I = I + 1 ;
  144. T = 'EXTRAIRE' LIS I ;
  145. SIGM = TAB.CONTRAINTES.I ; ;
  146. *IST SIGM ;
  147. EPSM = TAB.DEFORMATIONS_INELASTIQUES.I ;
  148. LTT = LTT ET (PROG T ) ;
  149. SXX = EXTRAIRE SIGM SMZZ 1 1 1 ;
  150. S = SXX;
  151.  
  152. EXX = EXTRAIRE EPSM EIZZ 1 1 1 ;
  153.  
  154.  
  155. LSXX = LSXX ET (PROG SXX ) ;
  156. LEXX = LEXX ET (PROG EXX ) ;
  157. BF1 = .222E-6; BF2 = .036797; BF3 = 6.0;
  158. RF1 = .4051E-4; RF2 = .031829; RF3 = 3.5;
  159. AF1 = 7.544324; AF2 = -.257143E-2; AF3 = .62202E-4;
  160. SF1 = .2266E-2; SF2 = 41.34;
  161.  
  162. ******* CALCUL ANALYTIQUE ************
  163. HB = BF2*S/BF3;
  164. HR = RF2*S/RF3;
  165.  
  166.  
  167. B = BF1*((((EXP(HB))-(EXP(HB*(-1.))))/2)**BF3);
  168.  
  169. R = RF1*((((EXP(HR))-(EXP(HR*(-1.))))/2)**RF3);
  170.  
  171. A1 = AF1*B/R;
  172.  
  173. GS=SF1/RF1*R;
  174.  
  175. SI ( S > SF2 );
  176. A2 = AF2+(AF3*S;
  177. SINON ;
  178. A2=0.;
  179. FINSI;
  180. EF = (A1*(1-(EXP(-1.*R*T)))) + (A2*(1-(EXP(-1.*gS*T)))) +(B*t);
  181.  
  182. DEFEQ = EXTRAIRE TAB.VARIABLES_INTERNES.I EPSE 1 1 1;
  183.  
  184. ldefeq = ldefeq et ( prog defeq);
  185.  
  186. LEF = LEF ET (PROG EF);
  187.  
  188. 'FIN' AFFICHE ;
  189.  
  190. I = 1;
  191. LER = PROG 0.;
  192.  
  193. REPETER BOUCLE ( dime LIS);
  194.  
  195. i = i + 1;
  196. exxi = extraire lexx i;
  197. efi = extraire lef i;
  198. eri = abs((exxi-efi)/efi*100);
  199. ler = ler et (prog eri);
  200.  
  201. FIN BOUCLE;
  202. LIST LER;
  203.  
  204. *__________________________________________________________________;
  205. EV1 = EVOL MANU LTT LSXX ;
  206. EV2 = EVOL MANU LTT LEXX ;
  207.  
  208. SI ( NEG GRAPH 'N');
  209.  
  210. DESS EV1 ;
  211. DESS EV2 ;
  212. FINSI;
  213. *__________________________________________________________________;
  214.  
  215. ermax = maxi ler;
  216. LIST ERMAX;
  217. SI (ermax <EG 5.);
  218. ERRE 0;
  219. SINON;
  220. ERRE 5;
  221. FINSI
  222.  
  223.  
  224. FIN; ;
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  

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