Télécharger la plaquette Cast3M


fr en

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