Télécharger plas7.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas7.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Plas7.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. * CAS TEST DU 91/06/13 PROVENANCE : TEST
  8.  
  9. *******************************************************
  10. * Test plas7.dgibi: Jeux de données *
  11. * --------------------------------- *
  12. *******************************************************
  13.  
  14. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  15. * SINON SI GRAPH DIFFERENT DE N TOUS
  16. * LES GRAPHIQUES SONT AFFICHES
  17.  
  18. GRAPH = 'N' ;
  19.  
  20. SAUT PAGE;
  21. SI (NEG GRAPH 'N') ;
  22. OPTI ECHO 1 ;
  23. OPTI TRAC X ;
  24. SINO ;
  25. OPTI ECHO 0 ;
  26. FINSI ;
  27.  
  28. SAUT PAGE;
  29. *
  30. *******************************************************
  31. * *
  32. * TEST PLAS7 *
  33. * *
  34. * RECTANGLE EN FLEXION PURE (DEFORMATIONS PLANES, *
  35. * PLASTICITE PARFAITE) *
  36. * *
  37. * Test V.P.C.S. numero SSNP11/89 *
  38. * *
  39. * References : R. Hill, Mathematical theory of *
  40. * plasticity, Oxford, Oxford Clavendon Press, 1950. *
  41. * *
  42. * Description : *
  43. * ----------- *
  44. * *
  45. * axe y *
  46. * D | A *
  47. * |----- ---- UA epaisseur a = 12.5 mm *
  48. * | | / hauteur b = 100.0 mm *
  49. * | | / *
  50. * | | / *
  51. * | |/ *
  52. * O |-----/-----> axe x *
  53. * | /| *
  54. * | / | *
  55. * | / | *
  56. * | / | *
  57. * ----- *
  58. * C UB B *
  59. * *
  60. * Chargement : Deplacement impose variant *
  61. * lineairement de UA a -UA entre A et B, *
  62. * UA croissant jusqu'a : UAmax = 0.5 mm . *
  63. * *
  64. * Conditions aux limites : *
  65. * sur CD : u = 0. *
  66. * au point O : v = 0. *
  67. * *
  68. * Proprietes des materiaux : *
  69. * Materiau elastique, parfaitement plastique *
  70. * E = 210 000 MPa , NU = 0.3 , SIGY = 483 MPa *
  71. * *
  72. *******************************************************
  73.  
  74. *
  75. OPTI ECHO 1;
  76. *
  77. OPTI MODE PLAN DEFO;
  78. OPTI ELEM QUA4;
  79. *
  80. TEMPS;
  81. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  82. * geometrie : maillage *
  83. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  84. *
  85. PA = 12.5E-3 50.E-3;
  86. PB = 12.5E-3 -50.E-3;
  87. PC = 0. -50.E-3;
  88. PD = 0. 50.E-3;
  89. *
  90. NB1 = 32; NB2 = 4;
  91. D1 = PA DROI NB1 PB;
  92. D2 = PB DROI NB2 PC;
  93. D3 = PC DROI NB1 PD;
  94. D4 = PD DROI NB2 PA;
  95. *
  96. SUR1 = DALL D1 D2 D3 D4 PLAN;
  97. *
  98. SI (NEG GRAPH 'N');
  99. TITR 'PLAS7 : MAILLAGE';
  100. TRAC 'QUAL' SUR1;
  101. FINSI;
  102. *
  103. NBNO1 = NBNO SUR1;
  104. NBEL1 = NBEL SUR1;
  105. *
  106. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  107. * definition du modele, des caracteristiques *
  108. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  109. *
  110. MOD1 = MODE SUR1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT QUA4;
  111. MAT1 = MATE MOD1 YOUN 2.1E11 NU 0.3 SIGY 483.E6;
  112. *
  113. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  114. * conditions aux limites *
  115. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  116. *
  117. CL1 = BLOQ UX D3;
  118. CL2 = BLOQ UY (SUR1 POIN PROC (0. 0.));
  119. CL3 = BLOQ UX D1;
  120. CL = CL1 ET CL2 ET CL3;
  121. *
  122. * Creation du chpoint de deplacement impose variable :
  123. * deplacement impose variant lineairement de u a -u
  124. * entre A et B, u croissant jusqu'a 0.5 mm .
  125. *
  126. I = 0.;
  127. LISD1 = PROG 1.;
  128. REPETER BLOC10 NB1;
  129. I = I + 1.;
  130. YY = 1. - (2. * (I / NB1));
  131. LISD1 = LISD1 ET (PROG YY);
  132. FIN BLOC10;
  133. *
  134. LISD1 = LISD1 * 1.E-3;
  135. CHPO1 = MANU CHPO (CHAN 'POI1' D1) 1 UX LISD1;
  136. *
  137. CHA1 = DEPI CL3 CHPO1;
  138. *
  139. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  140. * objets relatifs a la procedure PASAPAS *
  141. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  142. *
  143. TAB1 = TABLE;
  144. *
  145. * On prend comme liste des pseudo-temps
  146. * le deplacement u.
  147. *
  148. LREEL1 = PROG 0.029435 0.1 0.2 0.3 0.4 0.5 1.0;
  149. *
  150. * Chargement
  151. *
  152. EVOL1 = EVOL MANU UA LREEL1 COEFF LREEL1;
  153. CHAR1 = CHAR DIMP CHA1 EVOL1;
  154. *
  155. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  156. * resolution *
  157. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  158. *
  159. TAB1.'CARACTERISTIQUES' = MAT1;
  160. TAB1.'CHARGEMENT' = CHAR1;
  161. TAB1.'MODELE' = MOD1;
  162. TAB1.'TEMPS_CALCULES' = LREEL1;
  163. TAB1.'BLOCAGES_MECANIQUES' = CL;
  164. PASAPAS TAB1;
  165. *
  166. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  167. * extraction des resultats *
  168. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  169. *
  170. * Initialisation
  171. *
  172. MOM1 = TABL;
  173. LISD2 = LISD1 * 50.;
  174. MOM2 = TABL;
  175. ERG = TABL;
  176. *
  177. MOM2.1 = 90.57;
  178. MOM2.2 = 135.2;
  179. MOM2.3 = 138.4;
  180. MOM2.4 = 139.0;
  181. MOM2.5 = 139.2;
  182. MOM2.6 = 139.3;
  183. MOM2.7 = 139.4;
  184. *
  185. MOM2.8 = 543.4;
  186. MOM2.9 = 557.7;
  187.  
  188. *Boucle pour chaque pas de chargement : fonction de UA.
  189. I = 0;
  190. REPETER BLOC1 7;
  191. I = I + 1;
  192. CONT1 = TAB1.CONTRAINTES.I;
  193. DEPL1 = TAB1.DEPLACEMENTS.I;
  194. CHPO10 = CHAN CHPO MOD1 CONT1;
  195. * Trace facultatif de la deformee
  196. * (verification du chargement)
  197. SI (I EGA 7);
  198. SI (NEG GRAPH 'N');
  199. TITR 'PLAS7 : DEFORMEE Increment :' I;
  200. DEF0 = DEFO SUR1 DEPL1 0.;
  201. DEF1 = DEFO SUR1 DEPL1 2. VERT;
  202. TRAC SUR1 (DEF0 ET DEF1);
  203. FINSI;
  204. FINSI;
  205.  
  206. * Calcul du moment de flexion du aux contraintes
  207. * smxx s'exercant le long de AB.
  208. LIST2 = PROG;
  209. J = 0;
  210. NB = NB1 + 1;
  211. REPETER BLOC4 NB;
  212. J = J + 1;
  213. D = EXTR LISD2 J;
  214. POI1 = SUR1 POIN PROC (12.5E-3 D);
  215. SM = EXTRA CHPO10 SMXX POI1;
  216. LIST2 = LIST2 ET (PROG SM);
  217. FIN BLOC4;
  218.  
  219. LIST3 = LIST2 * LISD2;
  220. TITR 'CONTRAINTE SMXX Increment :' I;
  221. ESMXX = EVOL MANU Y LISD2 SMXX LIST2;
  222. TITR 'PRODUIT Y * SMXX Increment :' I;
  223. EYSMXX = EVOL MANU Y LISD2 YSMXX LIST3;
  224.  
  225. SI (I EGA 7);
  226. SI (NEG GRAPH 'N');
  227. DESS ESMXX;
  228. DESS EYSMXX;
  229. FINSI;
  230. FINSI;
  231.  
  232. MOM1.I = INTG EYSMXX;
  233. MOM1.I = EXTR MOM1.I 1;
  234. MOM1.I = -1.E-4 * MOM1.I;
  235.  
  236. * Calcul de l'ecart entre theorie et calcul.
  237. ERG.I = 100. * (ABS ((MOM1.I - MOM2.I) / MOM2.I));
  238.  
  239. * Extraction de la contrainte maximale SMXX
  240. * pour les pas 1 et 7.
  241. SI (I EGA 1);
  242. MOM1.8 = EXTR CHPO10 SMXX PA;
  243. MOM1.8 = MOM1.8 / 1.E6;
  244. ERG.8 = 100. * (ABS ((MOM1.8 - MOM2.8) / MOM2.8));
  245. FINSI;
  246.  
  247. SI (I EGA 7);
  248. MOM1.9 = EXTR CHPO10 SMXX (D1 POIN 2);
  249. MOM1.9 = MOM1.9 / 1.E6;
  250. ERG.9 = 100. * (ABS ((MOM1.9 - MOM2.9) / MOM2.9));
  251. FINSI;
  252.  
  253. SI (NEG GRAPH 'N');
  254. TITR 'SMXX Increment :' I;
  255. TRAC MOD1 (EXCO CONT1 SMXX);
  256. FINSI;
  257.  
  258. FIN BLOC1;
  259.  
  260. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  261. * affichage des resultats *
  262. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  263. MESS 'RESULTATS : PLAS7';
  264. MESS '-------------------';
  265.  
  266. MESS 'Moment de flexion du aux contraintes SMXX le long de AB.';
  267. MESS '--------------------------------------------------------';
  268.  
  269. MESS ' UAmax = 0.029435 mm :';
  270. MESS ' ------------------- ';
  271. MESS ' Moment theorique :' MOM2.1 'kN.m';
  272. MESS ' Moment calcule :' MOM1.1 'kN.m';
  273. MESS ' Soit un ecart de :' ERG.1 '%';
  274.  
  275. MESS ' UAmax = 0.1 mm :';
  276. MESS ' -------------- ';
  277. MESS ' Moment theorique :' MOM2.2 'kN.m';
  278. MESS ' Moment calcule :' MOM1.2 'kN.m';
  279. MESS ' Soit un ecart de :' ERG.2 '%';
  280.  
  281. MESS ' UAmax = 0.2 mm :';
  282. MESS ' --------------';
  283. MESS ' Moment theorique :' MOM2.3 'kN.m';
  284. MESS ' Moment calcule :' MOM1.3 'kN.m';
  285. MESS ' Soit un ecart de :' ERG.3 '%';
  286.  
  287. MESS ' UAmax = 0.3 mm :';
  288. MESS ' --------------';
  289. MESS ' Moment theorique :' MOM2.4 'kN.m';
  290. MESS ' Moment calcule :' MOM1.4 'kN.m';
  291. MESS ' Soit un ecart de :' ERG.4 '%';
  292.  
  293. MESS ' UAmax = 0.4 mm :';
  294. MESS ' --------------';
  295. MESS ' Moment theorique :' MOM2.5 'kN.m';
  296. MESS ' Moment calcule :' MOM1.5 'kN.m';
  297. MESS ' Soit un ecart de :' ERG.5 '%';
  298.  
  299. MESS ' UAmax = 0.5 mm :';
  300. MESS ' --------------';
  301. MESS ' Moment theorique :' MOM2.6 'kN.m';
  302. MESS ' Moment calcule :' MOM1.6 'kN.m';
  303. MESS ' Soit un ecart de :' ERG.6 '%';
  304.  
  305. MESS ' UAmax = 1.0 mm :';
  306. MESS ' ---------------';
  307. MESS ' Moment theorique :' MOM2.7 'kN.m';
  308. MESS ' Moment calcule :' MOM1.7 'kN.m';
  309. MESS ' Soit un ecart de :' ERG.7 '%';
  310.  
  311. MESS 'Contrainte maximale smxx :';
  312. MESS '------------------------ :';
  313. MESS ' UAmax = 0.029435 mm';
  314. MESS ' -------------------';
  315. MESS ' Valeur theorique :' MOM2.8 'MPa';
  316. MESS ' Valeur calculee :' MOM1.8 'MPa';
  317. MESS ' Soit un ecart de :' ERG.8 '%';
  318.  
  319. MESS ' UAmax = 1.0 mm';
  320. MESS ' --------------';
  321. MESS ' Valeur theorique :' MOM2.9 'MPa';
  322. MESS ' Valeur calculee :' MOM1.9 'MPa';
  323. MESS ' Soit un ecart de :' ERG.9 '%';
  324.  
  325. MESS 'Donnees sur le maillage';
  326. MESS '-----------------------';
  327. MESS ' Nombre d elements :' NBEL1 'QUA4';
  328. MESS ' Nombre de noeuds :' NBNO1 ;
  329.  
  330. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  331. * code fonctionnement *
  332. *+++++++++++++++++++++++++++++++++++++++++++++++++++++*
  333. ERGMAX = MAXI (PROG ERG.1 ERG.2 ERG.3 ERG.4 ERG.5 ERG.6 ERG.7 ERG.8 ERG.9);
  334.  
  335. SI (ERGMAX <EG 5.0);
  336. ERRE 0;
  337. SINON;
  338. ERRE 5;
  339. FINSI;
  340.  
  341. TEMPS;
  342.  
  343. FIN;
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  

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