Télécharger thpl5.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : thpl5.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. MESS '*********************************************************';
  14. MESS '* COMMISSION V.P.C.S. : THERMO-PLASTICITE *';
  15. MESS '* *';
  16. MESS '*CYLINDRE SOUMIS A UNE REPARTITION LOGARITHMIQUE *';
  17. MESS '* DE TEMPERATURE. *';
  18. MESS '*CALCUL EN DEFORMATIONS PLANES *';
  19. MESS '* *';
  20. MESS '*FONCTIONS TESTEES : ELASTO-PLASTICITE AVEC ECROUISSAGE *';
  21. MESS '* SOUS CHARGEMENT THERMIQUE. *';
  22. MESS '* *';
  23. MESS '*********************************************************';
  24. *
  25. OPTI MODE PLAN DEFO;
  26. OPTI ELEM QUA8;
  27. OPTI TRAC BENS;
  28. GRAPH='N';
  29.  
  30. *++++++++++++++++++++++++++++++++*
  31. * GEOMETRIE DU CYLINDRE *
  32. *++++++++++++++++++++++++++++++++*
  33. *
  34. DENS 1. ; XA = 1. ; YA = 0. ; PA = XA YA;
  35. DENS 1.2; XB = 2. ; YB = 0. ; PB = XB YB;
  36. DENS 1. ; XC = 0. ; YC = 1. ; PC = XC YC;
  37. DENS 1.2; XD = 0. ; YD = 2. ; PD = XD YD;
  38. DENS 1. ; XO = 0. ; YO = 0. ; PO = XO YO ;
  39. DENS 1.1; pz = 0. 1.497; pm = 1.497 0. ;
  40.  
  41.  
  42. DAm = PA droi -7 pm;
  43. DmB = Pm DROI -3 PB;
  44. DDz = PD DROI -3 Pz;
  45. DzC = Pz DROI -7 pc;
  46. DAB = DAm et DmB;
  47. DDC = DDz et DzC;
  48. *
  49. si complet;
  50. N2 = 20;
  51. sinon;
  52. n2 = 10;
  53. finsi;
  54. SUR1 = DAB ROTA N2 90 PO;
  55.  
  56. *
  57. NOSUR1 = NBNO SUR1;
  58. ELSUR1 = NBEL SUR1;
  59. *
  60. ELIM (SUR1 ET DAB ET DDC) 0.0001;
  61. *
  62. SI(NEG GRAPH 'N');
  63. TITR ' SSNP05 : MAILLAGE PAR QUA8 (NOEUDS)';
  64. TRACE NOEUD SUR1;
  65. TITR ' SSNP05 : MAILLAGE PAR QUA8 (ELEMENTS)';
  66. TRACE ELEM SUR1;
  67. FINSI;
  68. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  69. * CONDITIONS AUX LIMITES *
  70. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  71. *
  72. CL1 = SYMT DEPL ROTA PA PB SUR1 0.0001;
  73. CL2 = SYMT DEPL ROTA PC PD SUR1 0.0001;
  74. CL = CL1 ET CL2 ;
  75. *
  76. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  77. * MODELES - MATERIAUX - CARACTERISTIQUES
  78. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  79. *
  80. MO = MODE sur1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;
  81.  
  82. MAT = MATE MO YOUN 200000 NU 0.3 ALPHA 12E-6 SIGY 300;
  83. *
  84. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  85. * CHAMPS RELATIFS A LA PROCEDURE *
  86. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  87. *
  88. * CHAMPS DE TEMPERATURE REPARTIE SUR SUR1 AUX DIFFERENTS
  89. * INSTANTS DE CHARGEMENT
  90. * 10 PAS DE CHARGEMENT ENTRE TE = 0C ET TE = 200C
  91. * (TE = TEMPERATURE EXTERIEURE).
  92. *
  93. *==================================================================*
  94. * TABLE DE TRAVAIL : CHPOINT DE TEMPERATURE A T0 = 0 ET T1 = 1. *
  95. *==================================================================*
  96. *
  97. T0 = 0.; T1 = 1. ; T2 = 2.; T3 = 3.; T4 = 4.; T5 = 5.;
  98. T6 = 6.; T7 = 7. ; T8 = 8.; T9 = 9.; T10 = 10.;
  99.  
  100. *
  101. * EXTRACTIONS DES CHPOINTS DES COORDONNEES X ET Y DES NOEUDS DE SUR1
  102. *
  103. CHPSURX = COOR 1 SUR1;
  104. CHPSURY = COOR 2 SUR1;
  105. CHPSUR = ((CHPSURX ** 2) + (CHPSURY ** 2)) ** 0.5;
  106. *
  107. * CREATION DU CHPOINT DE TEMPERATURE POUR TE = 1C
  108. *
  109. CONST1 = 1. / (LOG 2.);
  110. CHP2 = (LOG CHPSUR) * CONST1;
  111. *
  112. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2
  113. * UTILISATION DE L'OPERATEUR VARI
  114. *
  115. LIS10 = PROG 0. 200.;
  116. LIS20 = PROG 0. 200.;
  117. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;
  118. *
  119. CHP1 = VARI CHP2 EVOL2 T;
  120. *
  121. * CHARGEMENT DE LA TABLE DES CHPOINTS DE TEMPERATURE
  122. *
  123. TEMP = TABLE;WCH1 = TABLE;
  124. TEMP . 0 = T0;TEMP . 1 = T1;
  125. TEMP . 2 = T2;TEMP . 3 = T3;
  126. TEMP . 4 = T4;TEMP . 5 = T5;
  127. TEMP . 6 = T6;TEMP . 7 = T8;
  128. TEMP . 8 = T8;TEMP . 9 = T9;
  129. TEMP . 10 = T10;
  130.  
  131. WCH1.0 = MANU CHPO SUR1 1 T 0. ;
  132. WCH1.1 = CHP1 * 120.;
  133. WCH1.2 = CHP1 * 130.;
  134. WCH1.3 = CHP1 * 140.;
  135. WCH1.4 = CHP1 * 150.;
  136. WCH1.5 = CHP1 * 160.;
  137. WCH1.6 = CHP1 * 170.;
  138. WCH1.7 = CHP1 * 180.;
  139. WCH1.8 = CHP1 * 190.;
  140. WCH1.9 = CHP1 * 195.;
  141. WCH1.10 = CHP1 * 200.;
  142.  
  143. CHA2 = CHAR 'T' TEMP WCH1;
  144. *
  145. *
  146. *
  147. * LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1
  148. *
  149. si complet;
  150. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10;
  151. sinon;
  152. LIS1 = PROG T0 T1 T2;
  153. finsi;
  154.  
  155. *
  156. * CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE
  157. *
  158. si complet;
  159. LFORCE = PROG 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ;
  160. sinon;
  161. LFORCE = PROG 0. 0. 0.;
  162. finsi;
  163.  
  164. FOR1 = FORCE FX 1. PA;
  165. CHPO1 = RESU FOR1;
  166. EVOL1 = EVOL MANU TEMPS LIS1 FORCE LFORCE;
  167. CHA1 = CHAR 'MECA' CHPO1 EVOL1;
  168. *
  169. *
  170. TAB1 = TABLE;
  171. TAB1.'MODELE' = MO ;
  172. TAB1.'CARACTERISTIQUES' = MAT;
  173. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  174. TAB1.'TEMPS_CALCULES' = LIS1;
  175. TAB1.'TEMPS_SAUVES' = LIS1;
  176. TAB1.'BLOCAGES_MECANIQUES' = CL;
  177. TAB1.'PROCESSEURS' = 'AUTOMATIQUE';
  178. *
  179. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  180. * PROCEDURE PASAPAS *
  181. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  182. *
  183. PASAPAS TAB1;
  184. *
  185. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  186. * EXTRACTION et AFFICHAGE des RESULTATS *
  187. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  188. *
  189. MESS ' RESULTATS DU TEST SSNP05';
  190. MESS ' ------------------------';
  191. MESS ' Les resultats de reference sont la moyenne des calculs';
  192. MESS ' effectues par Framasoft et le Cetim';
  193.  
  194. *
  195. MESS ' Contraintes interpolees (MPa)';
  196. MESS ' -----------------------------';
  197. *
  198. * Interpolations des contraintes aux noeuds
  199. *
  200. si complet;
  201. NN10 = 10;
  202. sinon;
  203. nn10 = 2;
  204. finsi;
  205. CHEL1 = TAB1.CONTRAINTES.NN10;
  206. CHPO1 = CHAN CHPO MO CHEL1;
  207.  
  208. * Extraction des contraintes
  209. *
  210. si complet ;
  211. SMRR10 = 0.1731; SMRR30 = 0.3806; SMRR20 = 47.04;
  212. SMTT10 = 340.6; SMTT30 = -191.75 ; SMTT20 = -50.63;
  213. SMZZ10 = 114.1; SMZZ30 = -345.95; SMZZ20 = -284. ;
  214. sinon;
  215. SMRR10 = .36175; SMRR30 = -2.9683 ; SMRR20 = 35.943;
  216. SMTT10 =272.68; SMTT30 = -169.49; SMTT20 = -21.166;
  217. SMZZ10 = 81.913; SMZZ30 = -348.10; SMZZ20 = -177.17;
  218. finsi;
  219. SMRR1 = EXTR CHPO1 SMXX PA;
  220. SMRR3 = EXTR CHPO1 SMXX PB;
  221. SMRR2 = EXTR CHPO1 SMXX pm;
  222. *
  223. SMTT1 = EXTR CHPO1 SMYY PA;
  224. SMTT3 = EXTR CHPO1 SMYY PB;
  225. SMTT2 = EXTR CHPO1 SMYY pm;
  226. * ;
  227. SMZZ1 = EXTR CHPO1 SMZZ PA;
  228. SMZZ3 = EXTR CHPO1 SMZZ PB;
  229. SMZZ2 = EXTR CHPO1 SMZZ pm;
  230. *
  231. * Calcul des ecarts avec resultats de reference
  232. *
  233. ERGRR1 = 100 * (ABS ((SMRR1 - SMRR10) / SMRR10));
  234. ERGRR2 = 100 * (ABS ((SMRR2 - SMRR20) / SMRR20));
  235. ERGRR3 = 100 * (ABS ((SMRR3 - SMRR30) / SMRR30));
  236. *
  237. ERGTT1 = 100 * (ABS ((SMTT1 - SMTT10) / SMTT10));
  238. ERGTT2 = 100 * (ABS ((SMTT2 - SMTT20) / SMTT20));
  239. ERGTT3 = 100 * (ABS ((SMTT3 - SMTT30) / SMTT30));
  240. *
  241. ERGZZ1 = 100 * (ABS ((SMZZ1 - SMZZ10) / SMZZ10));
  242. ERGZZ2 = 100 * (ABS ((SMZZ2 - SMZZ20) / SMZZ20));
  243. ERGZZ3 = 100 * (ABS ((SMZZ3 - SMZZ30) / SMZZ30));
  244. *
  245. MESS ' Resultats pour r = 1. mm ';
  246. MESS ' ------------------------ ';
  247. *
  248. MESS ' Contrainte radiale de reference (smrr) :' SMRR10 'MPa';
  249. MESS ' Contrainte radiale calculee (smrr) :' SMRR1 'MPa';
  250. MESS ' Soit un ecart de :' ERGRR1 '%';
  251. MESS ' ECART NON-SIGNIFICATIF';
  252. *
  253. MESS ' Contrainte tangentielle de reference :' SMTT10 'MPa';
  254. MESS ' Contrainte tangentielle calculee (smtt) :' SMTT1 'MPa';
  255. MESS ' Soit un ecart de :' ERGTT1 '%';
  256. *
  257. MESS ' Contrainte axiale de reference :' SMZZ10 'MPa';
  258. MESS ' Contrainte axiale calculee (smzz) :' SMZZ1 'MPa';
  259. MESS ' Soit un ecart de :' ERGZZ1 '%';
  260. *
  261. MESS ' Resultats pour r = 1.497 mm ';
  262. MESS ' ---------------------------- ';
  263. *
  264. MESS ' Contrainte radiale de reference :' SMRR20 'MPa';
  265. MESS ' Contrainte radiale calculee (smrr) :' SMRR2 'MPa';
  266. MESS ' Soit un ecart de :' ERGRR2 '%';
  267. *
  268. MESS ' Contrainte tangentielle de reference :' SMTT20 'MPa';
  269. MESS ' Contrainte tangentielle calculee (smtt) :' SMTT2 'MPa';
  270. MESS ' Soit un ecart de :' ERGTT2 '%';
  271. *
  272. MESS ' Contrainte axiale de reference :' SMZZ20 'MPa';
  273. MESS ' Contrainte axiale calculee (smzz) :' SMZZ2 'MPa';
  274. MESS ' Soit un ecart de :' ERGZZ2 '%';
  275.  
  276. *
  277. MESS ' Resultats pour r = 2. mm ';
  278. MESS ' ------------------------ ';
  279. *
  280. MESS ' Contrainte radiale de reference :' SMRR30 'MPa';
  281. MESS ' Contrainte radiale calculee (smrr) :' SMRR3 'MPa';
  282. MESS ' Soit un ecart de :' ERGRR3 '%';
  283. MESS ' ECART NON-SIGNIFICATIF';
  284. *
  285. MESS ' Contrainte tangentielle de reference :' SMTT30 'MPa';
  286. MESS ' Contrainte tangentielle calculee (smtt) :' SMTT3 'MPa';
  287. MESS ' Soit un ecart de :' ERGTT3 '%';
  288. *
  289. MESS ' Contrainte axiale de reference :' SMZZ30 'MPa';
  290. MESS ' Contrainte axiale calculee (smzz) :' SMZZ3 'MPa';
  291. MESS ' Soit un ecart de :' ERGZZ3 '%';
  292. *
  293. MESS ' Deplacements radiaux';
  294. MESS ' --------------------';
  295. *
  296. * Extractions des deplacements radiaux
  297. *
  298. DEP1 = TAB1.DEPLACEMENTS.NN10;
  299. *
  300. UR1 = EXTR DEP1 UX PA;
  301. UR2 = EXTR DEP1 UX PB;
  302. UR3 = EXTR DEP1 UX pm;
  303. mess ur1 ur2 ur3;
  304. *
  305. si complet;
  306. UR10 = 0.1804E-2;
  307. UR20 = 0.386E-2;
  308. sinon;
  309. ur10 = 1.24001E-03 ;ur20 = 2.48285E-03;
  310. finsi;
  311. *UR30 = NON CALCULE
  312. *
  313. * Calcul des ecarts avec solution V.P.C.S.
  314. *
  315. ERGUR1 = 100 * ( ABS (( UR1 - UR10) / UR10 ));
  316. ERGUR2 = 100 * ( ABS (( UR2 - UR20) / UR20 ));
  317. *ERGUR3 = 100 * ( ABS (( UR3 - UR30) / UR30 ));
  318. *
  319. MESS ' Deplacement de reference radial pour r = 1. mm :' UR10 'mm';
  320. MESS ' Deplacement calcul radial pour r = 1. mm :' UR1 'mm';
  321. MESS ' Soit un ecart de :' ERGUR1 '%';
  322. SAUT 1 LIGN;;
  323. *
  324. *MESS ' Deplacement de reference radial pour r = 1.497 mm :' UR30 'mm';
  325. MESS ' Deplacement calcul radial pour r = 1.497 mm :' UR3 'mm';
  326. *MESS ' Soit un ecart de :' ERGUR3 '%';
  327. *
  328. MESS ' Deplacement de reference radial pour r = 2. mm :' UR20 'mm';
  329. MESS ' Deplacement calcul radial pour r = 2. mm :' UR2 'mm';
  330. MESS ' Soit un ecart de :' ERGUR2 '%';
  331. *
  332. MESS ' Nombre de noeuds : ' NOSUR1;
  333. MESS ' Nombre d elements : ' ELSUR1;
  334. *
  335. * Tracer de la deformee.
  336. *
  337. SI (NEG GRAPH 'N');
  338. DEF0 = DEFO SUR1 DEP1 0. ;
  339. DEF1 = DEFO SUR1 DEP1 10. ROUGE ;
  340. TITR ' SSNP05 : DEFORMEE (DEFO. PLANE)';
  341. TRAC (DEF0 ET DEF1);
  342. FINSI;
  343. *
  344. * tracer de la contrainte equivalente de Von-Mises
  345. *
  346. SI (NEG GRAPH 'N');
  347. CHEL2 = VMIS mo CHEL1;
  348. TITR ' TRACE DES CONTRAINTES DE VON-MISES';
  349. TRAC CHEL2 mo SUR1;
  350. LLIS2 = PROG 0. 290. 295. 300. 305. 310. 600.;
  351. TITR 'LIEU OU IL Y A PLASTIFICATION';
  352. TRAC LLIS2 CHEL2 mo SUR1;
  353. FINSI;
  354. *
  355. TEMPS;
  356. *
  357. *===================================================================*
  358. * Code fonctionnement *
  359. *===================================================================*
  360. *
  361. * Maximum des erreurs relatives
  362. *
  363. ERG1 = MAXI (PROG ergtt1 ergzz1 ergrr2 ergtt2 ergzz2 ergtt3 ergzz3);
  364. *
  365. * remarque : on ne tient pas compte des ecarts sur la contrainte
  366. * radiale pour r = 1.mm et r = 2.mm
  367. *
  368. *
  369. ERG2 = MAXI (PROG ERGur1 ergur2);
  370. ERGMAX = MAXI (PROG ERG1 ERG2);
  371. *
  372. * Il y a erreur si un des pourcentage erg*** est suprieur a 5%.
  373. *
  374. SI (ERGMAX <EG 5.);
  375. ERRE 0 ;
  376. SINON;
  377. ERRE 5 ;
  378. FINSI;
  379. *
  380. FIN;
  381. *
  382.  
  383.  
  384.  
  385.  
  386.  
  387.  
  388.  
  389.  
  390.  
  391.  
  392.  
  393.  
  394.  
  395.  

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