Télécharger la plaquette Cast3M


fr en

Télécharger vpla3.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  vpla3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test vpla3.dgibi: jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7.  
  8. *******************************************************
  9. *           Test vpla3.dgibi: jeux de données         *
  10. *           ---------------------------------         *
  11. *******************************************************
  12.  
  13. *  POUR CALCUL COMPLET METTRE COMPLET À : VRAI;
  14.  
  15. COMPLET = faux;
  16.  
  17. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  18. * SINON SI GRAPH DIFFERENT DE N TOUS
  19. * LES GRAPHIQUES SONT AFFICHES
  20.  
  21. GRAPH = 'N' ;
  22.  
  23. SAUT PAGE;
  24. SI (NEG GRAPH 'N') ;
  25.   OPTI ECHO 1   ;
  26.   OPTI TRAC X   ;
  27. SINO ;
  28.   OPTI ECHO 1   ;
  29. FINSI ;
  30.  
  31. SAUT PAGE;
  32. *-----------------------------------------------------------------------
  33. *                           TEST VPLA3
  34. *
  35. * SORTIE DU DOMAINE ELASTIQUE ET PHASE VISCO-PLASTIQUE
  36. * COMPORTEMENT VISCO-PLASTIQUE - MODELE GUIONNET - CAS AXISYMETRIQUE
  37. * CYLINDRE CREUX SOUMIS A PRESSION INTERNE CONSTANTE ET
  38. *                         DEFORMATION AXIALE IMPOSEE CYCLIQUE VARIANT DE
  39. *                         0.35% A -0.35%
  40. * COMPARAISON DES BOUCLES SIGMA AXIALE FONCTION DE DEFORMATION AXIALE
  41. * AVEC RESULTATS EXPERIMENTAUX (RAPPORT DMT/93-434)
  42. * COMPARAISON DE LA DEFORMATION CIRCONFERENTIELLE PROGRESSIVE EN
  43. * FONCTION DU NOMBRE DE CYCLES AVEC RESULTATS EXPERIMENTAUX
  44. *-----------------------------------------------------------------------
  45. OPTI DIME 2 ELEM QUA8 MODE AXISYMETRIQUE ;
  46.  
  47. * -----DEFINITION DE LA GEOMETRIE
  48. RAY  = 8. ;
  49. EP   = 2. ;
  50. LONG = EP ;
  51.  
  52. P0 =  RAY      0.     ;
  53. P1 =  RAY      LONG   ;
  54. P2 = (RAY+EP)  LONG   ;
  55. P3 = (RAY+EP)  0.     ;
  56. *
  57. S01 = P0 DROIT 1 P1   ;
  58. S12 = P1 DROIT 1 P2   ;
  59. S23 = P2 DROIT 1 P3   ;
  60. S30 = P3 DROIT 1 P0   ;
  61. POUT= DALL S01 S12 S23 S30 'PLAN' ;
  62.  
  63.  
  64. SI (NEG GRAPH 'N') ;
  65.   TRAC 'QUAL' POUT ;
  66. FINSI ;
  67.  
  68. * -----DEFINITION DES CONDITIONS AUX LIMITES-----
  69. * -----     ET DES DEPLACEMENTS IMPOSES     -----
  70. CL1 = BLOQ UZ S12 ;
  71. CL2 = BLOQ UZ S30 ;
  72. CL  = CL1 ET CL2 ;
  73.  
  74. DELTAL=LONG * 0.35E-2;
  75. FL1 = DEPI CL1 DELTAL ;
  76.  
  77. MO = MODE POUT 'MECANIQUE' 'ELASTIQUE' 'VISCOPLASTIQUE' 'GUIONNET' ;
  78. MA = MATE MO   'YOUN' 192500.  'NU  ' 0.3     'N   ' 24     'K   ' 200.
  79.                'CK  ' 0.       'R0  ' 25.     'CD  ' 0.     'M   ' 0.8
  80.                'A   ' 24000.   'A1  ' 3.3E5   'C   ' 180.   'C1  ' 40.
  81.                'G   ' 0.       'R   ' 0.      'NN  ' 0.005  'C2  ' 40.
  82.                'Q   ' 1.466    'G1  ' 0.      'R1  ' 0.     'BETA' 1.1
  83.                'C0  ' 180.     'P1M0' 0.00333 ;
  84.  
  85. * -----CHARGEMENT EN TRACTION/COMPRESSION
  86.  
  87. *      MONTEE EN CHARGE ET 1ER CYCLE
  88. L1 = PROG 0. 1.  3. 5.   ;
  89. L2 = PROG 0. 1. -1. 1.   ;
  90.  
  91. *        AU TOTAL 10 CYCLES
  92. SI   COMPLET;
  93.   REPETER STRU1 9;
  94.     NL1 = DIME L1 ;
  95.     X   = EXTR L1 NL1 ;
  96.     L1  = L1 ET (PROG (X+2.) (X+4.)) ;
  97.     L2  = L2 ET (PROG -1 1)          ;
  98.   FIN STRU1        ;
  99. FINSI;
  100.  
  101. EV = EVOL MANU 'T' L1 'FT' L2 ;
  102. CHA1 = CHAR 'DIMP' FL1 EV ;
  103.  
  104. * -----CHARGEMENT EN PRESSION
  105. MOP = 'MODE' S01 'CHARGEMENT' 'PRESSION' ;
  106. MOTOT = MO 'ET' MOP ;
  107. *MATOT = MA ET MAP ;
  108. MATOT = MA ;
  109. *
  110. EVPRE = 'EVOL' 'MANU' 'TEMP' ('PROG' 0. 1. 100000.)
  111.                       'PRES' ('PROG' 0. 1. 1.) ;
  112. VALP = pres MOP 'PRES' -21.5 ;
  113. CHA2 = 'CHAR' 'PRES' VALP EVPRE ;
  114.  
  115. * -----MISE EN CHARGE
  116. LIS0 = PROG .2 .3 .4 PAS .05 .6 PAS .02 .8 PAS .01 1 ;
  117. SI COMPLET;
  118.  
  119. * -----PREMIER CYCLE
  120.   LIS1 = PROG 1.01 1.02 PAS 0.02 1.1 PAS 0.05 2.9 PAS 0.02 2.98 2.99 3.
  121.               3.01 3.02 PAS 0.02 3.1 PAS 0.05 4.9 PAS 0.02 4.98 4.99 5.;
  122.  
  123. * -----CYCLE 2 A 10
  124.   NLIS   = DIME LIS1    ;
  125.   LISCTE4=PROG NLIS * 4.;
  126.   LIS2 = LIS1 + LISCTE4 ;
  127.   LIS3 = LIS2 + LISCTE4 ;
  128.   LIS4 = LIS3 + LISCTE4 ;
  129.   LIS5 = LIS4 + LISCTE4 ;
  130.   LIS6 = LIS5 + LISCTE4 ;
  131.   LIS7 = LIS6 + LISCTE4 ;
  132.   LIS8 = LIS7 + LISCTE4 ;
  133.   LIS9 = LIS8 + LISCTE4 ;
  134.   LIS10= LIS9 + LISCTE4 ;
  135.  
  136.   LIS = LIS0 ET LIS1 ET LIS2 ET LIS3 ET LIS4 ET LIS5 ET LIS6 ET LIS7
  137.              ET LIS8 ET LIS9 ET LIS10;
  138. SINON;
  139.   LIS = LIS0;
  140. FINSI;
  141.  
  142. TAB = TABL ;
  143. TAB.'PRECISION'           = 1E-3;
  144. TAB.'BLOCAGES_MECANIQUES' = CL  ;
  145. TAB.'CARACTERISTIQUES'    = MATOT ;
  146. TAB.'CHARGEMENT'          = CHA1 'ET' CHA2 ;
  147. TAB.'TEMPS_CALCULES'      = LIS ;
  148. TAB.'MODELE'              = MOTOT  ;
  149. PASAPAS TAB ;
  150.  
  151. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  152. * COMPARAISON DES RESULTATS AVEC RESULTATS EXPERIMENTAUX
  153. *                RAPPORT DMT/93-434
  154. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  155. J = 0;
  156. PROGSIG = PROG;
  157. PROGEPS = PROG;
  158. SI COMPLET;
  159.   LECTI = LECT 37 85 133 181 229 277 325 469 517 949 997;
  160. SINON;
  161.   LECTI = LECT 37;
  162. FINSI;
  163.  
  164. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  165. * CONTRAINTE AXIALE EN FIN DE MONTEE EN CHARGE
  166. *                   EN MILIEU ET FIN DES CYCLES 1 2 ET 3
  167. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  168. SI COMPLET;
  169.   NFOIS = (DIME LECTI) - 4 ;
  170. SINON;
  171.   NFOIS = 1;
  172. FINSI;
  173. REPETER STRU1 NFOIS;
  174.   J = J + 1;
  175.   I = EXTR LECTI J;
  176.   SIG = TAB.CONTRAINTES.I;
  177.   SIG11 = EXTR SIG SMZZ 1 1 9;
  178.   PROGSIG = PROGSIG ET (PROG SIG11);
  179. FIN STRU1;
  180.  
  181. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  182. * DEFORMATION CIRCONFERENTIELLE MOYENNE (EN %) POUR CYCLES 5 ET 10
  183. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  184. SI COMPLET;
  185.   REPETER STRU2 2;
  186.     J = J + 1;
  187.     I1 = EXTR LECTI J    ;
  188.     I2 = EXTR LECTI (J+1);
  189.     EPS_1   = EPSI MOTOT (TAB.DEPLACEMENTS.I1);
  190.     EPS_2   = EPSI MOTOT (TAB.DEPLACEMENTS.I2);
  191.     EPS22_1 = ( (EXTR EPS_1 EPTT 1 1 9) + (EXTR EPS_1 EPTT 1 1 6) +
  192.                 (EXTR EPS_1 EPTT 1 1 3) ) / 3.;
  193.     EPS22_2 = ( (EXTR EPS_2 EPTT 1 1 9) + (EXTR EPS_2 EPTT 1 1 6) +
  194.                 (EXTR EPS_2 EPTT 1 1 3) ) / 3.;
  195.     EPS22 = (EPS22_1 + EPS22_2) / 2. ;
  196.     PROGEPS = PROGEPS ET (PROG EPS22);
  197.     J = J + 1;
  198.   FIN STRU2;
  199.  
  200.   PROGEPS = 100. * PROGEPS;
  201. FINSI;
  202.  
  203. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  204. *              RESULTATS EXPERIMENTAUX
  205. *                RAPPORT DMT/93-434
  206. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++
  207. SI COMPLET;
  208.   PROGSX = PROG 300 -300 330.3 -309.8 345.9 -312.8 351.9 ;
  209.   PROGEX = PROG 0.22 0.2675 ;
  210. SINON;
  211.   PROGSX = PROG 300 ;
  212. FINSI;
  213.  
  214. *++++++++++++++
  215. *   COMPARAISON
  216. *++++++++++++++
  217. ERR_SIG = ABS ( (PROGSIG - PROGSX) / PROGSX ) ;
  218. ERR_SIG = MAXI ERR_SIG                        ;
  219. MESS 'Erreur sur les Contraintes : ' ERR_SIG  ;
  220.  
  221. SI COMPLET ;
  222.   ERRADM = 0.18 ;
  223.  
  224.   ERR_EPS = ABS ( (PROGEPS - PROGEX) / PROGEX ) ;
  225.   ERR_EPS = MAXI ERR_EPS                        ;
  226.   MESS 'Erreur sur les Deformations : ' ERR_EPS ;
  227.   SI (ERR_EPS <EG 0.18);
  228.     ERRE 0;
  229.   SINON;
  230.     ERRE 5;
  231.   FINSI;
  232. SINON;
  233.  
  234.   ERRADM=0.03;
  235.   SI (ERR_SIG <EG ERRADM);
  236.     ERRE 0;
  237.   SINON;
  238.     ERRE 5;
  239.   FINSI;
  240. FINSI;
  241.  
  242. FIN;
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  

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