Télécharger plas11.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas11.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Plas11.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTION ECHO 1;
  8. * ajout de option epsilon lineaire pour la precision des test!
  9. OPTION epsilon lineaire;
  10. GRAPH = 'N';
  11. SAUT PAGE;
  12. *
  13. *******************************************************
  14. * *
  15. * TEST : PLAS11 *
  16. * DEFORMATIONS PLANES GENERALISEES *
  17. * *
  18. * une poutre de section rectangulaire est *
  19. * soumise à une rotation imposée RX : *
  20. * *
  21. * Hauteur:4 Largeur:2 *
  22. * *
  23. * ''''''''' *
  24. * ' * ' *
  25. * ' | ' cette poutre est modélisée *
  26. * --------- avec une partie massive (1x1) *
  27. * | | une autre avec des éléments *
  28. * | -------->@RX COQ2 (long :.5 epais 2.) et *
  29. * | | la dernière avec des éléments *
  30. * --------- point POI1 : 4 point, 2 en *
  31. * ' | ' haut, 2 en bas (sect .5 ) *
  32. * ' * ' *
  33. * ''''''''' *
  34. * *
  35. * On étudie le comportement du moment résultant *
  36. * de cette rotation imposée jusqu'au domaine *
  37. * plastique *
  38. * *
  39. *******************************************************
  40. TITRE 'PLAS11 - SECTION DE POUTRE SOUMISE A UNE ROTATION IMPOSEE';
  41. OPTION ECHO 1;
  42. TEMPS;
  43. *
  44. *
  45. OPTI ELEM QUA8;
  46. *
  47. * on définit le point support des déformations
  48. * planes généralisées
  49. *
  50. PG = 0. 0.;PG1=PG PLUS (0. 0.);PG2 = PG PLUS (0. 0.);
  51. *
  52. OPTI MODE PLAN GENE ;
  53. *
  54. * quelques initialisations :
  55. *
  56. I = 1;
  57. EVL1 = PROG 0. ;
  58. EVL2 = PROG 0. ;
  59. EVL3 = PROG 0. ;
  60. *
  61. * et quelques constantes pour le calcul théorique :
  62. *
  63. * demi hauteur de la section
  64. H = 2.;
  65. * demi largeur de la section
  66. B = 1. ;
  67. * le module de young
  68. E1 = 2E11;
  69. * la limite élastique du matériau
  70. SILIM = 4.83E8;
  71. * le moment d'inertie par rapport a l'axe des X
  72. IX = (2.*B)*((2.*H)**3)/12.;
  73. * la valeur limite de la déformation élastique
  74. EPSEL = SILIM/E1;
  75. * la valeur limite du moment élastique
  76. MOMEL = E1*IX*EPSEL/H;
  77. * le moment plastique ultime
  78. MOLIM = 3*MOMEL/2.;
  79. *
  80. *======================================================
  81. * géometrie : maillage
  82. *======================================================
  83. *
  84. * dimension en mètres
  85. *
  86. *----------maillage partie massive-----------*
  87. *
  88. * points
  89. *
  90. P1 = 1. -1.;
  91. P2 = 1. 1.;
  92. P3 = -1. 1.;
  93. P4 = -1. -1.;
  94. *
  95. * droite et lignes
  96. *
  97. * coefficients de finesse du maillage
  98. *
  99. N = 4;
  100. M = 1;
  101. *
  102. D1 = P1 D N P2 ;
  103. D2 = P2 D M P3 ;
  104. D3 = P3 D N P4 ;
  105. D4 = P4 D M P1 ;
  106. *
  107. SURF1 = DALL D1 D2 D3 D4;
  108. *
  109. *----------maillage éléments coq2------------*
  110. *
  111. OPTI ELEM SEG2 ;
  112. *
  113. P5 = 0. 1. ;
  114. P6 = 0. 1.5 ;
  115. P7 = 0. -1. ;
  116. P8 = 0. -1.5 ;
  117. *
  118. D5 = P5 D 2 P6;
  119. D6 = P7 D 2 P8;
  120. COQU1 = D5 ET D6 ;
  121. *
  122. P9 = 0. 1.625;
  123. P10 = 0. 1.875;
  124. P11 = 0. -1.625;
  125. P12 = 0. -1.875;
  126. *
  127. *----------maillage éléments poi1------------*
  128. *
  129. PO11 = MANU POI1 P9;
  130. PO12 = MANU POI1 P10;
  131. PO13 = MANU POI1 P11;
  132. PO14 = MANU POI1 P12;
  133. PO1 = PO11 ET PO12 ET PO13 ET PO14;
  134. *
  135. SURTOT = SURF1 ET COQU1 ET PO1 ;
  136. *
  137. SI (NEG GRAPH 'N');
  138. TITR 'PLAS12 : MAILLAGE';
  139. TRAC QUAL SURTOT;
  140. FINSI;
  141. *
  142. ELTOT = NBEL SURTOT;
  143. NOTOT = NBNO SURTOT;
  144. *
  145. *
  146. *======================================================
  147. * modèle - matériau
  148. * rigidité
  149. *======================================================
  150. *
  151. *----------------partie massive-----------------*
  152. *
  153. MO1 = MODE SURF1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT QUA8 DPGE PG;
  154. MA1 = MATE MO1 YOUN E1 NU 0. SIGY SILIM;
  155. ma2= MATE MO1 YOUN 'TOTO' NU 0. SIGY SILIM;
  156. list ma2;
  157. *
  158. *----------------éléments coq2------------------*
  159. *
  160. MO2 = MODE COQU1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT COQ2 DPGE PG1;
  161. MAT2 = MATE MO2 YOUN E1 NU 0. SIGY SILIM;
  162. CAR2 = CARA MO2 EPAI 2. ;
  163. MA2 = MAT2 ET CAR2 ;
  164. *
  165. *-----------------éléments poi1-----------------*
  166. *
  167. MO3 = MODE PO1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT POI1 DPGE PG1;
  168. MAT3 = MATE MO3 YOUN E1 NU 0. SIGY SILIM;
  169. CAR3 = CARA MO3 SECT .5 ;
  170. MA3 = MAT3 ET CAR3;
  171. *
  172. MODE1 = MO1 ET MO2 ET MO3;
  173. MATE1 = MA1 ET MA2 ET MA3;
  174. *
  175. RI1 = RIGI MODE1 MATE1;
  176.  
  177. *
  178. *======================================================
  179. * conditions aux limites
  180. *======================================================
  181. *
  182. P15 = SURF1 POINT PROCHE (0. 0.);
  183. P16 = SURF1 POINT PROCHE (0. 1.);
  184. B1 = BLOQ DEPL P15;
  185. B2 = BLOQ UX P16;
  186. B3 = BLOQ DEPL PO1 ;
  187. B4 = BLOQ UX COQU1 ;
  188. B5 = BLOQ UY (P5 ET P7) ;
  189. B6 = BLOQ UZ RY (PG et pg1 );
  190. B7 = BLOQ RX ( PG et pg1 );
  191. BTOT = B1 ET B2 ET B3 ET B4 ET B5 ET B6 ET B7;
  192. RI = RI1 ET BTOT;
  193. *
  194. *======================================================
  195. * chargement
  196. *======================================================
  197. *
  198. TEMP1 = (PROG 4. PAS 2. 30.) ET (PROG 50.);
  199. LAM1 = (PROG 4. PAS 2. 30.) ET (PROG 50.);
  200. FAT1 = .0002;
  201. TITRE 'PLAS11 : CHARGEMENT';
  202. EV1 = EVOL MANU 'TEMPS' TEMP1 'RX(T)' LAM1 ;
  203. DE1 = DEPI B7 FAT1 ;
  204. CHRIX = CHAR 'DIMP' DE1 EV1 ;
  205. *
  206. SI (NEG GRAPH 'N');
  207. DESS EV1 ;
  208. FINSI;
  209. *
  210. *======================================================
  211. * calcul non linéaire
  212. *======================================================
  213. *
  214. opti sauv 'plas11.sauv';
  215. TAB1 = TABLE;
  216. TAB1.'BLOCAGES_MECANIQUES' = BTOT;
  217. TAB1.'MODELE' = MODE1;
  218. TAB1.'CARACTERISTIQUES' = MATE1;
  219. TAB1.'CHARGEMENT' = CHRIX;
  220. TAB1.'TEMPS_CALCULES' = TEMP1;
  221. tab1.'TEMPS_SAUVEGARDES' = TEMP1;
  222. tab1.'PRECISION' = 1.e-6;
  223.  
  224. PASAPAS TAB1 ;
  225.  
  226.  
  227.  
  228.  
  229. *
  230. * recalcul des trois dernier pas (pour tester l'option reprise de
  231. * pasapas
  232. *
  233. tab1.REPRISE=26.;
  234. tab1.'TEMPS_CALCULES'=prog 28. 30. 50.;
  235. tab1.'TEMPS_SAUVEGARDES'=prog 28. 30. 50.;
  236. pasapas tab1;
  237.  
  238.  
  239. *
  240. *======================================================
  241. * post-traitement : tracé d'une courbe et comparaison
  242. * avec la solution analytique
  243. *======================================================
  244. *
  245. DIM1 = DIMENSION TEMP1;
  246. *
  247. REPETER BLOC1 DIM1;
  248. *
  249. TPS = EXTR TEMP1 I;
  250. CH1 = PECHE TAB1 DEPLACEMENTS TPS;
  251. MOIMP = REAC RI CH1;
  252. *
  253. *EXTRACTION DU MOMENT MX CORRESPONDANT À LA ROTATION RX
  254. *
  255. MMX = EXTR MOIMP MX PG;MMX1= EXTR MOIMP MX PG1;MMX=MMX+MMX1;
  256. RRX = EXTR CH1 RX PG;RRX1= EXTR CH1 RX PG1;RRX=RRX+RRX1;
  257. EVL1 = EVL1 ET (PROG MMX );
  258. EVL2 = EVL2 ET (PROG RRX );
  259. *
  260. *---------------SOLUTION THEORIQUE--------------------*
  261. *
  262. EPSN = 2.*RRX ;
  263. *
  264. * quand on est dans le domaine plastique :
  265. *
  266. SI (EPSN > EPSEL) ;
  267. MOMN = (1 - ((( EPSEL / EPSN) ** 2) / 3)) * MOLIM ;
  268. EVL3 = EVL3 ET (PROG MOMN);
  269. *
  270. * quand on est dans le domaine élastique :
  271. *
  272. SINON;
  273. MOMN = E1 * IX * RRX ;
  274. EVL3 = EVL3 ET (PROG MOMN) ;
  275. FINSI;
  276. *
  277. I = I+1;
  278. *
  279. FIN BLOC1;
  280. *
  281. SI (NEG GRAPH 'N');
  282. *
  283. COURB = TABLE;
  284. COURB.1 = 'MARQ ETOI NOLI' ;
  285. *
  286. titre 'PLAS11 : COURBE MX EN FONCTION DE RX ';
  287. EV1 = EVOL MANU 'RX' EVL2 'MX' EVL1 ;
  288. EV2 = EVOL MANU 'RX' EVL2 'MX' EVL3 ;
  289. * courbe théorique : traits continus
  290. * point calculés : marques étoiles
  291. DESS (EV1 ET EV2) COURB LOGO;
  292. FINSI;
  293. *
  294. * calcul de l'écart avec la solution de référence
  295. *
  296. ERG1 = 100*(ABS((MOMN-MMX)/MOMN));
  297. *
  298. *======================================================
  299. * affichage des résultats
  300. *======================================================
  301. *
  302. MESS ' RESULTATS : TEST PLAS11 ';
  303. MESS '---------------------------';
  304. *
  305. MESS 'on cherche la valeur du moment calculé dans le ';
  306. MESS 'domaine plastique';
  307. MESS 'elle est comparée a une valeur obtenue analytiquement';
  308. MESS ' moment théorique :' MOMN 'MPa';
  309. MESS ' moment calculé :' MMX 'MPa';
  310. MESS ' Soit un écart de :' ERG1 '%';
  311. *
  312. MESS ' Nombre d elements :' ELTOT;
  313. MESS ' Nombre de noeuds :' NOTOT;
  314. MESS ' soit ' (NOTOT * 3) 'd.d.l.';
  315. *
  316. *
  317. *======================================================
  318. * code de fonctionnement
  319. *======================================================
  320. *
  321. * L'écart maximum entre valeur théorique et
  322. * calculée doit être inférieur a 0.5%.
  323. *
  324. SI (ERG1 <EG 0.5);
  325. ERRE 0;
  326. SINON;
  327. ERRE 5;
  328. FINSI;
  329. *
  330. *======================================================
  331. * Temps de calcul et fin
  332. *======================================================
  333. *
  334. TEMPS;
  335. *
  336. FIN;
  337. *
  338.  
  339.  
  340.  
  341.  
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  

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