Télécharger thpl4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : thpl4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux;
  5. *
  6. *************************************************
  7. * pour calcul complet mettre complet à : vrai;
  8. ***************************************************
  9. * CAS TEST DU 92/03/20 PROVENANCE : TC1
  10. * CAS TEST DU 92/03/19 PROVENANCE : PHIL
  11. OPTI ECHO 0;
  12. *
  13. *********************************************************************
  14. * COMMISSION V.P.C.S. : THERMO-PLASTICITE *
  15. * THPL4.DGIBI *
  16. * *
  17. * SPHERE CREUSE EPAISSE EN THERMO-PLASTICITE *
  18. * CALCUL AXISYMETRIQUE *
  19. * *
  20. * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE SOUS *
  21. * CHARGEMENT THERMIQUE. *
  22. * *
  23. *********************************************************************
  24.  
  25. OPTION DIME 2 ;
  26. OPTION MODE AXIS;
  27. OPTION TRACE BENS;
  28. GRAPH='N';
  29. OPTION ELEM QUA8;
  30.  
  31. *
  32. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  33. * GEOMETRIE DU CYLINDRE *
  34. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  35. *
  36. DENS 1. ; RA = 1. ; ZA = 0. ; PA = RA ZA;
  37. dens 1.1; RM = 1.5; ZM = 0. ; PM1 = RM ZM;
  38. DENS 1.11; RB = 2. ; ZB = 0. ; PB = RB ZB;
  39. DENS 1. ; RC = 0. ; ZC = 1. ; PC = RC ZC;
  40. DENS 1.1; PM2 = 0. 1.5;
  41. DENS 1.11; RD = 0. ; ZD = 2. ; PD = RD ZD;
  42. DENS 1. ; RO = 0. ; ZO = 0. ; PO = RO ZO ;
  43. *
  44.  
  45. DAM1 = PA DROI -7 PM1; DM1B = PM1 D -4 PB; DAB = DAM1 et DM1B;
  46. DDM2 = PD DROI -4 PM2; DM2C = PM2 D -7 PC; DDC = DDM2 et DM2C;
  47. *
  48. N2 = 10;
  49. SUR1 = DAB ROTA N2 90 PO;
  50. *
  51. ELIM SUR1 0.0001;
  52. *
  53. SI(NEG GRAPH 'N');
  54. TITR ' SPHERE CREUSE EPAISSE : MAILLAGE';
  55. TRACE SUR1;
  56.  
  57. FINSI;
  58.  
  59. list ( NBNO SUR1);
  60. list ( NBEL SUR1);
  61.  
  62. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  63. * CONDITIONS AUX LIMITES *
  64. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  65. *
  66. CL1 = SYMT DEPL ROTA PA PB SUR1 0.0001;
  67. CL2 = SYMT DEPL ROTA PC PD SUR1 0.0001;
  68. CL = CL1 ET CL2 ;
  69. *
  70. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  71. * MODELES - MATERIAUX - CARACTERISTIQUES
  72. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  73. *
  74. MO = MODE SUR1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;
  75.  
  76. MAT = MATE MO YOUN 200e3 NU 0.3 ALPHA 1.E-5 SIGY 200 'TALP' 0. 'TREF' 0.;
  77. *
  78. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  79. * LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1
  80. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  81. *
  82. * LE NOMBRE DE PAS DE CALCUL (PAS DE CHARGEMENT), CORRESPOND AU
  83. * NOMBRE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS.
  84. * UN PAS ELASTIQUE (DEBUT PLASTIFICATION) JUSQU'A 100C.
  85. * UN PAS PLASTIQUE TOUT LES 10C APRES.
  86. *
  87. *
  88. T0 = 0.; T1 = 1.; T2 = 2.; T3 = 3.;
  89. T4 = 4.; T5 = 5.; T6 = 6.; T7 = 7.;
  90. T8 = 8.; T9 = 9.; T10 = 10.; T11 = 11.;
  91. T12 = 12.; T13 = 13.; T14 = 14.; T15 = 15.;
  92. T16 = 16.; T17 = 17.; T18 = 18.;
  93. *
  94. si complet;
  95. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12 T13
  96. T14 T15 T16 T17 T18;
  97. sinon;
  98. LIS1 = PROG T0 T1 ;
  99. finsi;
  100. *
  101. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  102. * CHAMPS RELATIFS A LA PROCEDURE NONLIN *
  103. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  104. *
  105. * CHAMPS DE TEMPERATURE REPARTIE SUR SUR1 AUX DIFFERENTS INSTANTS
  106. * DE LA LISTE OU L'ON VEUT DES RESULTATS (CF. LIST1)
  107. *
  108. *===================================================================*
  109. * *
  110. * CALCUL D'UN CHPOINT CHP1 DE TEMPERATURE POUR UN GRADIENT DE *
  111. * TEMPERATURE DE UN DEGRE CELSIUS, SOIT DT = 1C *
  112. * *
  113. * DONNEES : - DT EST CROISSANT DE 0 @ DTMAX, *
  114. * - POUR A=1. < R < B=2., LA REPARTITION DE TEMPERATURE *
  115. * EST : DT = T(A) - T(B) > 0. *
  116. * T(R) = DT/R * AB/(A-B) *
  117. * *
  118. * CALCUL POUR DT = 1C, DE LA REPARTITION DE TEMPERATURE *
  119. * SOIT T(R) = 2. / R *
  120. * *
  121. *===================================================================*
  122. *
  123. * EXTRACTION DES CHPOINTS DES COORDONNEES R ET Z DU MAILLAGE SUR1.
  124. *
  125. CHPSUR1 = COOR 1 SUR1;
  126. CHPSUR2 = COOR 2 SUR1;
  127. *
  128. * CREATION DU CHPOINT DE TEMPERATURE POUR DT = 1C
  129. *
  130. CHPSUR = ( (CHPSUR1**2 ) + (CHPSUR2**2) )**(-0.5);
  131. CHP2 = 2. * CHPSUR;
  132. *
  133. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2
  134. * UTILISATION DE L'OPERATEUR VARI
  135. *
  136. LIS10 = PROG 0. 1000.;
  137. LIS20 = PROG 0. 1000.;
  138. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;
  139. *
  140. CHP1 = VARI CHP2 EVOL2 T;
  141. *
  142. * CHARGEMENT DE LA TABLE DES CHPOINTS DE TEMPERATURE
  143. *
  144. *
  145. TEMP = TABLE;TEMPER = TABLE;
  146. TEMP . 0 = T0;TEMP . 1 = T1;
  147. TEMP . 2 = T2;TEMP . 3 = T3;
  148. TEMP . 4 = T4;TEMP . 5 = T5;
  149. TEMP . 6 = T6;TEMP . 7 = T8;
  150. TEMP . 8 = T8;TEMP . 9 = T9;
  151. TEMP . 10 = T10;TEMP . 11 = T11;
  152. TEMP . 12 = T12;TEMP . 13 = T13;
  153. TEMP . 14 = T14;TEMP . 15 = T15;
  154. TEMP . 16 = T16;TEMP . 17 = T17;
  155. TEMP . 18 = T18;
  156.  
  157. TEMPER . 0 = MANU CHPO SUR1 1 T 0.;DT = 98.;
  158. TEMPER . 1 = (CHP1 * DT);DT = 110.;
  159. TEMPER . 2 = (CHP1 * DT);DT = 120.;
  160. TEMPER . 3 = (CHP1 * DT);DT = 130.;
  161. TEMPER . 4 = (CHP1 * DT);DT = 140.;
  162. TEMPER . 5 = (CHP1 * DT);DT = 150.;
  163. TEMPER . 6 = (CHP1 * DT);DT = 160.;
  164. TEMPER . 7 = (CHP1 * DT);DT = 170.;
  165. TEMPER . 8 = (CHP1 * DT);DT = 180.;
  166. TEMPER . 9 = (CHP1 * DT);DT = 188.;
  167. TEMPER . 10 = (CHP1 * DT);DT = 196.;
  168. TEMPER . 11 = (CHP1 * DT);DT = 210.;
  169. TEMPER . 12 = (CHP1 * DT);DT = 220.;
  170. TEMPER . 13 = (CHP1 * DT);DT = 230.;
  171. TEMPER . 14 = (CHP1 * DT);DT = 240.;
  172. TEMPER . 15 = (CHP1 * DT);DT = 250.;
  173. TEMPER . 16 = (CHP1 * DT);DT = 255.;
  174. TEMPER . 17 = (CHP1 * DT);DT = 264.83;
  175. TEMPER . 18 = (CHP1 * DT);
  176.  
  177. CHA2 = CHAR 'T' TEMP TEMPER;
  178.  
  179. *
  180. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  181. * CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE *
  182. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  183. *
  184. LFORCE = PROG 0. 0.;
  185. LLIS1 = PROG 0. 18.;
  186. *
  187. FOR1 = FORCE FR 1. PA;
  188. CHPO1 = RESU FOR1;
  189. EVOL1 = EVOL MANU TEMPS LLIS1 FORCE LFORCE;
  190. CHA1 = CHAR 'MECA' CHPO1 EVOL1;
  191. *
  192. *
  193. TAB1 = TABLE;
  194. TAB1.'MODELE' = MO ;
  195. TAB1.'CARACTERISTIQUES' = MAT;
  196. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  197. TAB1.'TEMPS_CALCULES' = LIS1;
  198. TAB1.'TEMPS_SAUVES' = LIS1;
  199. TAB1.'BLOCAGES_MECANIQUES' = CL;
  200.  
  201. *
  202. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  203. * PROCEDURES PASAPAS *
  204. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  205. *
  206. PASAPAS TAB1;
  207.  
  208. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  209. * EXTRACTION DES RESULTATS *
  210. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  211. *
  212. MESS ' **************************************************************';
  213. MESS ' * COMMISSION V.P.C.S. : THERMO-PLASTICITE *';
  214. MESS ' * FICHIER : SSNA06PLAPN.DGIBI *';
  215. MESS ' * *';
  216. MESS ' * SPHERE CREUSE EPAISSE EN THERMO-PLASTICITE *';
  217. MESS ' * CALCUL AXISYMETRIQUE *';
  218. MESS ' * *';
  219. MESS ' * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE *';
  220. MESS ' * SOUS CHARGEMENT THERMIQUE. *';
  221. MESS ' * *';
  222. MESS ' **************************************************************';
  223. MESS ' * CAS PLASTICITE SANS ECROUISSAGE ISOTROPE LINEAIRE *';
  224. MESS ' **************************************************************';
  225. *
  226. MESS ' Les rsultats thoriques sont tirs de l article';
  227. MESS ' "Thermo-elastic behaviour of strain-hardening thich-walled';
  228. MESS ' sphere", D Durban, International journal of solids structures';
  229. MESS ' vol 19, nx7, pp 643-652, 1983';
  230. *===================================================================*
  231. * resultats a 98 degres : indice T1
  232. *===================================================================*
  233. *
  234. * on cherche pour r = 1.mm, 1.5mm, et 2.mm les valeurs suivantes :
  235. * smrr (contrainte radiale), smtt (contrainte tangentielle)
  236. * et ur (deplacement radial).
  237. * La temperature de 98 degres celsius correspond a la premiere
  238. * plastification.
  239. *
  240. NN1 = 1; NN2 = 2;
  241. NN3 = 3; NN4 = 4;
  242. NN5 = 5; NN6 = 6;
  243. NN7 = 7; NN8 = 8;
  244. NN9 = 9; NN10 = 10;
  245. NN11 = 11; NN12 = 12;
  246. NN13 = 13; NN14 = 14;
  247. NN15 = 15; NN16 = 16;
  248. NN17 = 17; NN18 = 18;
  249.  
  250. MESS ' *------------------------------------------------*';
  251. MESS ' * Resultats a 98 degres : premiere plastifiction *';
  252. MESS ' *------------------------------------------------*';
  253. *
  254. * Deplacements
  255. *
  256. DEP98 = TAB1.DEPLACEMENTS.NN1;
  257. *
  258. UR98 = TABLE; UR980 = TABLE;
  259. UR98.1 = EXTR DEP98 UR PA; UR980.1 = 1.26E-3;
  260. UR98.2 = EXTR DEP98 UR PM1; UR980.2 = 2.08E-3;
  261. UR98.3 = EXTR DEP98 UR PB; UR980.3 = 2.52E-3;
  262. *
  263. ERGUR98 = TABLE;
  264. I = 0;
  265. REPETER BLOC1 3;
  266. I = I + 1;
  267. SI (I &lt;EG 1) ; MESS ' r = 1.0 mm';FINSI;
  268. SI ((I &lt;EG 2) ET (I >EG 2)); MESS ' r = 1.5 mm';FINSI;
  269. SI (I >EG 3); MESS ' r = 2.0 mm'; FINSI;
  270. *
  271. ERGUR98.I = 100 * ( ABS (( UR98.I - UR980.I) / (UR980.I+1e-6)));
  272. MESS ' Deplacement radial theorique :' UR980.I 'mm';
  273. MESS ' Deplacement radial calcule :' UR98.I 'mm';
  274. MESS ' Soit un ecart de : ' ERGUR98.I '%';
  275. SAUT 1 LIGN;
  276. FIN BLOC1;
  277. *
  278. * Contraintes
  279. *
  280. CHEL1 = TAB1.CONTRAINTES.NN1;
  281. CONT98 = CHAN CHPO MO CHEL1;
  282. *
  283. SMRR98 = TABLE; SMRR980 = TABLE;
  284. SMRR98.1 = EXTR CONT98 SMRR PA; SMRR980.1 = 0.;
  285. SMRR98.2 = EXTR CONT98 SMRR PM1; SMRR980.2 = -38.52;
  286. SMRR98.3 = EXTR CONT98 SMRR PB; SMRR980.3 = 0.;
  287. *
  288. SMTT98 = TABLE; SMTT980 = TABLE;
  289. SMTT98.1 = EXTR CONT98 SMTT PA; SMTT980.1 = -200.;
  290. SMTT98.2 = EXTR CONT98 SMTT PM1; SMTT980.2 = 5.93;
  291. SMTT98.3 = EXTR CONT98 SMTT PB; SMTT980.3 = 80.;
  292. *
  293. ERGSRR98 = TABLE; ERGSTT98 = TABLE;
  294. I = 0;
  295. REPETER BLOC2 3;
  296. I = I + 1;
  297. SI (I &lt;EG 1); MESS ' r = 1.0 mm'; FINSI;
  298. SI ((I &lt;EG 2) ET (I >EG 2)); MESS ' r = 1.5 mm';FINSI;
  299. SI (I >EG 3); MESS ' r = 2.0 mm'; FINSI;
  300. *
  301. ERGSRR98.I = 100 * ( ABS (( SMRR98.I - SMRR980.I ) /
  302. ( SMRR980.I+1e-6 )));
  303.  
  304. MESS ' Contrainte radiale theorique :' SMRR980.I 'MPa';
  305. MESS ' Contrainte radiale calculee :' SMRR98.I 'MPa';
  306. MESS ' Soit un ecart de : ' ERGSRR98.I '%';
  307. si ((( abs smrr980.i ) &lt;eg 10.) et ((abs smrr98.i) &lt;eg 10.));
  308. mess 'CET ECART EST NON-SIGNIFICATIF';
  309. finsi;
  310.  
  311. SAUT 1 LIGN;
  312. *
  313. ERGSTT98.I = 100 * ( ABS (( SMTT98.I - SMTT980.I ) /
  314. ( SMTT980.I+1e-6)));
  315. MESS ' Contrainte tangentielle theorique :' SMTT980.I 'MPa';
  316. MESS ' Contrainte tangentielle calculee :' SMTT98.I 'MPa';
  317. MESS ' Soit un ecart de : ' ERGSTT98.I '%';
  318. si ((( abs smtt980.i ) &lt;eg 10.) et ((abs smtt98.i) &lt;eg 10.));
  319. mess 'CET ECART EST NON-SIGNIFICATIF';
  320. finsi;
  321.  
  322. SAUT 1 LIGN;
  323. FIN BLOC2;
  324. si complet;
  325. *
  326. *===================================================================*
  327. * resultats a 196 degres : indice T11
  328. *===================================================================*
  329. *
  330. * On cherche - pour r = 1. mm et r = 2. mm smtt et smrr,
  331. * - la valeur du rayon correspondant a la limite de
  332. * plastification, et pour ce rayon, on veut smtt et smrr. En
  333. * remarque on peut dire que smeq = 200.MPa = smtt + smrr.
  334. * - la valeur du rayon ou il y a un etat de contrainte
  335. * hydrostatique, et la valeur de cette contrainte.
  336. *
  337. MESS ' *------------------------*';
  338. MESS ' * Resultats a 196 degres *';
  339. MESS ' *------------------------*';
  340. *
  341. * smtt et smrr pour r = 1. mm et r = 2. mm
  342. *
  343. CHEL1 = TAB1.CONTRAINTES.NN11; CHV1 = VMIS MO CHEL1;
  344. CONT196 = CHAN CHPO MO chel1;
  345. *
  346. SMTT196 = TABLE; SMTT1960 = TABLE;
  347. ERTT196 = TABLE;
  348. SMTT196.1 = EXTR CONT196 SMTT PA; SMTT1960.1 = -200.;
  349. SMTT196.2 = EXTR CONT196 SMTT PB; SMTT1960.2 = 153.82;
  350. *
  351. SMRR196 = TABLE; SMRR1960 = TABLE;
  352. SMRR196.1 = EXTR CONT196 SMRR PA; SMRR1960.1 = 0.;
  353. SMRR196.2 = EXTR CONT196 SMRR PB;
  354. * SMRR1960.2 = 'NON-CONNUE';
  355. *
  356. MESS ' r = 1.0 mm';
  357. MESS ' Contrainte radiale theorique :' SMRR1960.1 'MPa';
  358. MESS ' Contrainte radiale calculee :' SMRR196.1 'MPa';
  359. si ((abs SMRR196.1) &lt;eg 5.);
  360. MESS 'ECART NON-SIGNIFICATIF';
  361. finsi;
  362. *
  363. MESS ' Contrainte tangentielle theorique :' SMTT1960.1 'MPa';
  364. MESS ' Contrainte tangentielle calculee :' SMTT196.1 'MPa';
  365. ERTT196.1 = 100 * ( ABS (( SMTT196.1 - SMTT1960.1 ) /
  366. ( SMTT1960.1+1e-6)));
  367. MESS ' Soit un ecart de : ' ERTT196.1 '%';
  368. *
  369. MESS ' r = 2.0 mm';
  370. * MESS ' Contrainte radiale theorique :' SMRR1960.1 'MPa';
  371. MESS ' Contrainte radiale calculee :' SMRR196.2 'MPa';
  372. *
  373. MESS ' Contrainte tangentielle theorique :' SMTT1960.2 'MPa';
  374. MESS ' Contrainte tangentielle calculee :' SMTT196.2 'MPa';
  375. ERTT196.2 = 100 * ( ABS (( SMTT196.2 - SMTT1960.2 ) /
  376. ( SMTT1960.2+1e-6)));
  377. MESS ' Soit un ecart de : ' ERTT196.2 '%';
  378. *
  379. * Recherche de la valeur du rayon de limite plastique
  380. *
  381. * On cherche cette valeur aux points de gauss pour les contraintes
  382. *
  383. * Creation d'un champs CHEL2 par elements des valeurs des rayons
  384. * (en spherique) aux points de gauss pour les contraintes. Alors le
  385. * champs CHEL2 a la meme structure que le champs de contrainte CHEL1.
  386. *
  387.  
  388. chsur1 = coor 1 chel1;
  389. chsur2 = coor 2 chel1;
  390. CHEL2 = ((CHSUR1 ** 2) + (CHSUR2 ** 2)) ** 0.5;
  391.  
  392.  
  393. * Boucle sur les elements numerotes de 1 a 11 de la sous-zone 1, la
  394. * valeur maximale etant inferieure a 1.5 mm. On affiche seulement le
  395. * rayon maximum (RAYON1MAX).
  396. *
  397. RAY1MAX = 1.;
  398. IELEM = 0;
  399. REPETER BLOC5 11;
  400. IELEM = IELEM + 1;
  401. *
  402. * Boucle sur les 9 points de gauss de l'element numerote IELEM
  403. *
  404. IPGAUS = 0;
  405. REPETER BLOC4 9;
  406. IPGAUS = IPGAUS + 1;
  407. SMTT1 = EXTR CHEL1 SMTT 1 IELEM IPGAUS;
  408. SMRR1 = EXTR CHEL1 SMRR 1 IELEM IPGAUS;
  409. SEQU1 = SMRR1 - SMTT1;
  410. *
  411. * test sur la contrainte equivalente : proche de 200 MPa ?
  412. *
  413. SI (( SEQU1 &lt;EG 205.) ET (SEQU1 >EG 195.));
  414. *
  415. * Dans ce cas, on cherche la valeur du rayon pour laquelle la contrainte
  416. * equivalente est proche de la contrainte limite plastique.
  417. *
  418. RAYON1 = EXTRA CHEL2 SCAL 1 IELEM IPGAUS;
  419. RAY1MAX = MAXI (PROG RAY1MAX RAYON1);
  420.  
  421. Si ( RAYON1 EGA RAY1MAX ) ;
  422. SMTTR1 = SMTT1;
  423. SMRRR1 = SMRR1;
  424. SMVMR1 = extr chv1 'SCAL' 1 IELEM IPGAUS;
  425. xel = extr CHSUR1 'SCAL' 1 IELEM IPGAUS;
  426. yel = extr CHSUR2 'SCAL' 1 IELEM IPGAUS;
  427.  
  428. finsi;
  429. FINSI;
  430. *
  431. * Fin boucle sur les points de gauss, et reinitialisation de IPGAUSS a 0.
  432. *
  433. FIN BLOC4;
  434. *
  435. * Fin de boucle sur les elements numerotes de 1 a 11;
  436. *
  437. FIN BLOC5;
  438. RAYTH1 = 1.1328;
  439. SMTTR10 = -249.88;
  440. SMRRR10 = -49.88;
  441. erra1 = 100*(abs ( ( RAYTH1 - RAY1MAX)/RAYTH1));
  442. erttra1 = 100*(abs ( (SMTTR1 - smttr10)/smttr10));
  443. errra1 = 100*(abs ( (SMRRR1 -smrrr10)/smrrr10));
  444.  
  445.  
  446. * Affichage de la limite plastique de la premiere zone.
  447. *
  448. MESS ' Rayon definisant la fin de la premiere zone plastique.';
  449. MESS ' ------------------------------------------------------';
  450. MESS ' Rayon1 calcul = ' RAY1MAX ' mm';
  451. MESS ' Rayon1 thorique = ' RAYTH1 ' mm';
  452. MESS ' ECART = ' erra1 '% ';
  453. MESS ' SMTT POUR CE RAYON =' SMTTR1 ' MPA';
  454. MESS ' SMTT thorique =' SMTTR10 ' MPA';
  455. MESS ' ECART =' erttra1 ' MPA';
  456. MESS ' SMRR POUR CE RAYON =' SMRRR1 ' MPA';
  457. MESS ' SMRR thorique =' SMRRR10 ' MPA';
  458. mESS ' ECART =' errra1 ' MPA';
  459. MESS ' CONTRAINTE DE VON MISES =' SMVMR1 'MPA';
  460. MESS 'LES COORDONNEES DU POINT DE GAUSS VALENT ';
  461. MESS ' X= ';list xel;
  462. MESS ' Y= ';list yel;
  463. SAUT 1 LIGN;
  464. *===================================================================*
  465. * resultats a 260 dgres : indice T18
  466. *===================================================================*
  467. *
  468. MESS ' *--------------------------------*';
  469. MESS ' * Resultats a 264 degres celsius *';
  470. MESS ' *--------------------------------*';
  471. MESS ' Apparition de la deuxieme zone plastique (r = 2.)';
  472. MESS ' On cherche la limite de la premiere zone plastique.';
  473. MESS ' La valeur theorique est proche de 1.1968mm';
  474. *
  475. * On cherche la limite de la premiere zone plastique.
  476. *
  477. CHEL1 = TAB1.CONTRAINTES.NN18; CP1 = chan CHPO mo chel1;
  478. cVM1 = VMIS mo chel1;
  479.  
  480.  
  481. RAY2MAX = 1.;
  482. IELEM = 0;
  483.  
  484.  
  485. REPETER BLOC5 10;
  486. IELEM = IELEM + 1;
  487. IPGAUS = 0;
  488. REPETER BLOC4 9;
  489. IPGAUS = IPGAUS + 1;
  490. SMTT1 = EXTR CHEL1 SMTT 1 IELEM IPGAUS;
  491. SMRR1 = EXTR CHEL1 SMRR 1 IELEM IPGAUS;
  492. SEQU1 = SMRR1 - SMTT1;
  493. *
  494. SI (( SEQU1 &lt;EG 201.) ET (SEQU1 >EG 199.));
  495. RAYON2 = EXTRA CHEL2 SCAL 1 IELEM IPGAUS;
  496. RAY2MAX = MAXI (PROG RAY2MAX RAYON2);
  497. * MESS ' Pour r =' RAYON2 'mm , sigma equivalente =' SEQU1 'MPa';
  498.  
  499. SI (RAYON2 EGA RAY2MAX);
  500. SMTTR2 = SMTT1;
  501. SMRRR2 = SMRR1;
  502. SMVMR2 = extr CVM1 'SCAL' 1 IELEM IPGAUS;
  503. xel = extr chsur1 'SCAL' 1 IELEM IPGAUS;
  504. yel = extr chsur2 'SCAL' 1 IELEM IPGAUS;
  505. FINSI;
  506.  
  507. FINSI;
  508. FIN BLOC4;
  509. FIN BLOC5;
  510. *
  511. * Affichage de la limite plastique de la premiere zone.
  512. *
  513. MESS ' Rayon definisant la fin de la premiere zone plastique.';
  514. MESS ' ------------------------------------------------------';
  515. MESS ' Rayon2 = ' RAY2MAX 'mm';
  516. RAY2TH = 1.1968;
  517. erra2 = 100*(abs (ray2th -ray2max)/ray2th);
  518. MESS ' RAYON2 thorique = ' RAY2TH ' mm';
  519. MESS 'ECART =' erra2;
  520.  
  521. *
  522. MESS 'AU POINT DE GAUSS,LES CONTRAINTES SONT';
  523. MESS 'SMTT =' SMTTR2 'MPa ;SMRR =' SMRRR2 'MPa';
  524. MESS 'SMVM =' SMVMR2 'MPa';
  525. MESS 'LES COORDONNEES DU POINT VALENT';
  526. MESS ' X= ' ;LIST xel;
  527. MESS ' Y= ' ;list yel;
  528.  
  529. TEMPS;
  530. *
  531.  
  532.  
  533. *
  534. * TRACE DES ISO-CONTRAINTES
  535. *
  536. CHELCONT = TABLE;
  537. CHELCONT.1 = TAB1.CONTRAINTES.NN1;
  538. CHELCONT.10 = TAB1.CONTRAINTES.NN10;
  539. CHELCONT.2 = TAB1.CONTRAINTES.NN2;
  540. CHELCONT.11 = TAB1.CONTRAINTES.NN11;
  541. CHELCONT.3 = TAB1.CONTRAINTES.NN3;
  542. CHELCONT.12 = TAB1.CONTRAINTES.NN12;
  543. CHELCONT.4 = TAB1.CONTRAINTES.NN4;
  544. CHELCONT.13 = TAB1.CONTRAINTES.NN13;
  545. CHELCONT.5 = TAB1.CONTRAINTES.NN5;
  546. CHELCONT.14 = TAB1.CONTRAINTES.NN14;
  547. CHELCONT.6 = TAB1.CONTRAINTES.NN6;
  548. CHELCONT.15 = TAB1.CONTRAINTES.NN15;
  549. CHELCONT.7 = TAB1.CONTRAINTES.NN7;
  550. CHELCONT.16 = TAB1.CONTRAINTES.NN16;
  551. CHELCONT.8 = TAB1.CONTRAINTES.NN8;
  552. CHELCONT.17 = TAB1.CONTRAINTES.NN17;
  553. CHELCONT.9 = TAB1.CONTRAINTES.NN9;
  554. CHELCONT.18 = TAB1.CONTRAINTES.NN18;
  555. *
  556. * RECHERCHE DES LIEUX OU IL Y A PLASTIFICATION : TRACEE
  557. SI (NEG GRAPH 'N');
  558.  
  559. CHEL2 = VMIS mo CHELCONT.1;
  560. TITR 'T= 98C : ISOCONTRAINTES DE VON-MISES';
  561. TRAC CHEL2 mo SUR1;
  562. *
  563. CHEL2 = VMIS mo CHELCONT.11;
  564. TITR 'T= 196C : ISOCONTRAINTES DE VON-MISES';
  565. TRAC CHEL2 mo SUR1;
  566. *
  567. CHEL2 = VMIS mo CHELCONT.18;
  568. TITR 'T= 265C : ISOCONTRAINTES DE VON-MISES';
  569. TRAC CHEL2 mo SUR1;
  570. *
  571. LLIS10 = PROG 195. 200. 205.;
  572. CHEL2 = VMIS mo CHELCONT.1;
  573. TITR 'T= 98C : ISOCONTRAINTES DE VON-MISES';
  574. TRAC LLIS10 CHEL2 mo SUR1;
  575. *
  576. CHEL2 = VMIS mo CHELCONT.11;
  577. TITR 'T= 196C : ISOCONTRAINTES DE VON-MISES';
  578. TRAC LLIS10 CHEL2 mo SUR1;
  579. *
  580. CHEL2 = VMIS mo CHELCONT.18;
  581. TITR 'T= 265C : ISOCONTRAINTES DE VON-MISES';
  582. TRAC LLIS10 CHEL2 mo SUR1;
  583. finsi;
  584. TEMPS ;
  585. finsi;
  586.  
  587. *****************************************************
  588. ****** CODE DE FONCTIONNEMENT *******
  589. *****************************************************
  590. * ON NE TIENT PAS COMPTE DES VALEURS NON-SIGNIFICATIVES *
  591. * IL Y A ERREUR SI UN DES POURCENTAGES EST SUPERIEUR A 5%*
  592. SUERG1 = PROG ERGUR98.1 ERGUR98.2 ERGUR98.3 ERGSRR98.2
  593. ERGSTT98.1 ERGSTT98.3;
  594. si complet;
  595. SUERG2 = PROG ERTT196.1 ERTT196.2 ERRA1 ERTTRA1 ERRRA1
  596. ERRA2;
  597. ERMA1 = MAXI SUERG1;
  598. ERMA2 = MAXI SUERG2;
  599. sinon;
  600. ERMA1 = MAXI SUERG1;
  601. ERMA2 = 0.;
  602. finsi;
  603. ERMA = MAXI ( PROG ERMA1 ERMA2);
  604. SI (ERMA &lt;EG 3.);
  605. ERRE 0;
  606. SINON;
  607. ERRE 5;
  608. FINSI;
  609. MESS 'L ECART MAXIMAL EST DE ' ERMA '%';
  610. FIN;
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  
  620.  
  621.  
  622.  

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