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 = VRAI;
  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 0 ;
  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. MAP = 'MATE' MOP 'PRES' -21.5 ;
  107. MOTOT = MO 'ET' MOP ;
  108. MATOT = MA ET MAP ;
  109. *
  110. EVPRE = 'EVOL' 'MANU' 'TEMP' ('PROG' 0. 1. 100000.)
  111. 'PRES' ('PROG' 0. 1. 1.) ;
  112. VALP = 'CHAN' 'TYPE' ('EXCO' MAP 'PRES') 'CONTRAINTES' ;
  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.  

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