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

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