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

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