Télécharger thpl3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : thpl3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 92/03/20 PROVENANCE : TC1
  5. * CAS TEST DU 92/03/19 PROVENANCE : PHIL
  6. opti echo 0;
  7. MESS ' ***************************************************************';
  8. MESS ' * COMMISSION V.P.C.S. : THERMO-PLASTICITE *';
  9. MESS ' * *';
  10. MESS ' * CYLINDRE SOUMIS A UNE REPARTITION LOGARITHMIQUE DE *';
  11. MESS ' * TEMPERATURE, CALCUL EN AXISYMETRIE. *';
  12. MESS ' * *';
  13. MESS ' * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE SOUS *';
  14. MESS ' * CHARGEMENT THERMIQUE. *';
  15. MESS ' * *';
  16. MESS ' ***************************************************************';
  17. MESS ' LES RESULTATS DE REFERENCE SONT LA MOYENNE DES CALCULS';
  18. MESS ' EFFECTUES PAR FRAMASOFT ET LE CETIM';
  19.  
  20. OPTION TRACE BENS;
  21. OPTION MODE AXIS;
  22. OPTION ELEM QUA8;
  23.  
  24. GRAPH='N';
  25. *
  26. *++++++++++++++++++++++++++++++++*
  27. * GEOMETRIE DU CYLINDRE *
  28. *++++++++++++++++++++++++++++++++*
  29. *
  30. DENS 1. ; RA = 1. ; ZA = 0. ; PA = RA ZA;
  31. DENS 1.2; RB = 2. ; ZB = 0. ; PB = RB ZB;
  32. dens 1.1; pm = 1.497 0;
  33.  
  34.  
  35. DENS 1. ; RC = 1. ; ZC = 0.1 ; PC = RC ZC;
  36. DENS 1.2; RD = 2. ; ZD = 0.1 ; PD = RD ZD;
  37. DENS 1.1; pz = 1.497 0.1;
  38. *
  39. DAm = PA DROIT -7 pm;
  40. DmB = pm droit -3 PB;
  41. DAB = DAm et DmB;
  42.  
  43. DBD = PB DROIT -1 PD;
  44. DDz = PD DROIT -3 Pz;
  45. DCA = PC DROIT -1 PA;
  46. DzC = Pz DROIT -7 PC;
  47. DDC = DDz et DzC;
  48.  
  49. *
  50. SUR1 = DALL DAB DBD DDC DCA PLAN ;
  51.  
  52.  
  53. ELIM SUR1 0.0001;
  54. *
  55. NOSUR1 = NBNO SUR1;
  56. ELSUR1 = NBEL SUR1;
  57. *
  58. SI (NEG GRAPH 'N');
  59. TITR ' SSNA04 : MAILLAGE PAR QUA8 (NOEUDS)';
  60. TRACE NOEUD SUR1;
  61. TITR ' SSNA04 : MAILLAGE PAR QUA8 (ELEMENTS)';
  62. TRACE ELEM SUR1;
  63. FINSI;
  64.  
  65. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  66. * CONDITIONS AUX LIMITES *
  67. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  68. *
  69. CL1 = BLOQUE DAB UZ ;
  70. CL2 = BLOQUE DDC UZ ;
  71. CL = CL1 ET CL2 ;
  72. *
  73. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  74. * MODELES - MATERIAUX - CARACTERISTIQUES
  75. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  76. *
  77. MO = MODE SUR1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;
  78.  
  79. MAT = MATE MO YOUN 200000 NU 0.3 ALPHA 12E-6 SIGY 300;
  80. *
  81. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  82. * CHAMPS RELATIFS A LA PROCEDURE PASAPAS *
  83. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  84. *
  85. *==============================================================*
  86. * CHAMPS DE TEMPERATURE REPARTIE SUR SUR1 *
  87. *==============================================================*
  88. *
  89. * 10 PAS DE CHARGEMENT ENTRE 0C ET 200C
  90. *
  91. T0 = 0.; T1 = 1. ; T2 = 2.; T3 = 3. ;
  92. T4 = 4.; T5 = 5. ; T6 = 6.; T7 = 7. ;
  93. T8 = 8.; T9 = 9. ; T10 = 10.;
  94. *
  95. * CREATION DU CHPOINT DE TEMPERATURE POUR TE = 1C (TEMP. EXTE.)
  96. *
  97. CHPSUR = COOR 1 SUR1;
  98. *
  99. CONST1 = 1. / (LOG 2.);
  100. CHP2 = (LOG CHPSUR) * CONST1;
  101.  
  102. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2
  103. * UTILISATION DE L'OPERATEUR VARI
  104. *
  105. LIS10 = PROG 0. 200.;
  106. LIS20 = PROG 0. 200.;
  107. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;
  108. *
  109. CHP1 = VARI CHP2 EVOL2 T;
  110.  
  111. *
  112. TEMP = TABLE;TEMPER = TABLE;
  113. TEMP . 0 = T0;TEMP . 1 = T1;
  114. TEMP . 2 = T2;TEMP . 3 = T3;
  115. TEMP . 4 = T4;TEMP . 5 = T5;
  116. TEMP . 6 = T6;TEMP . 7 = T8;
  117. TEMP . 8 = T8;TEMP . 9 = T9;
  118. TEMP . 10 = T10;
  119.  
  120. TEMPER . 0 = MANU CHPO SUR1 1 T 0.;
  121. TEMPER . 1 = (CHP1 * 20.);
  122. TEMPER . 2 = (CHP1 * 40.);
  123. TEMPER . 3 = (CHP1 * 60.);
  124. TEMPER . 4 = (CHP1 * 80.);
  125. TEMPER . 5 = (CHP1 * 100.);
  126. TEMPER . 6 = (CHP1 * 120.);
  127. TEMPER . 7 = (CHP1 * 140.);
  128. TEMPER . 8 = (CHP1 * 160.);
  129. TEMPER . 9 = (CHP1 * 180.);
  130. TEMPER . 10 = (CHP1 * 200.);
  131.  
  132. CHA2 = CHAR 'T' TEMP TEMPER;
  133.  
  134. *
  135. *================================================================
  136. * LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1
  137. *================================================================
  138. *
  139. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10;
  140. *
  141. *================================================================
  142. * CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE *
  143. *================================================================
  144. *
  145. LIS2 = PROG T0 T10;
  146. LFORCE = PROG 0. 0. ;
  147. FOR1 = FORCE FZ 1. PA;
  148. CHPO1 = RESU FOR1;
  149. EVOL1 = EVOL MANU TEMPS LIS2 FORCE LFORCE;
  150. CHA1 = CHAR 'MECA' CHPO1 EVOL1;
  151.  
  152. *================================================================
  153. * REMPLISSAGE DE LA TABLE POUR PASAPAS *
  154. *================================================================
  155. *
  156. TAB1 = TABLE;
  157. TAB1.'MODELE' = MO ;
  158. TAB1.'CARACTERISTIQUES' = MAT;
  159. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  160. TAB1.'TEMPS_CALCULES' = LIS1;
  161. TAB1.'TEMPS_SAUVES' = LIS1;
  162. TAB1.'BLOCAGES_MECANIQUES' = CL;
  163.  
  164. *
  165. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  166. * PROCEDURE PASAPAS *
  167. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  168. *
  169. PASAPAS TAB1;
  170. *
  171. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  172. * EXTRACTION DES RESULTATS *
  173. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  174.  
  175. NN = 10;
  176. CHEL1 = TAB1.CONTRAINTES.NN;
  177.  
  178.  
  179.  
  180. CHPO1 = CHAN CHPO MO CHEL1;
  181.  
  182. CHEL2 = VMIS MO CHEL1;
  183. *
  184. SI (NEG GRAPH 'N');
  185. TITR ' TRACE DES CONTRAINTES DE VON-MISES';
  186. TRAC CHEL2 MO SUR1;
  187. LLIS2 = PROG 0. 290. 295. 300. 305. 310. 600.;
  188. TITR 'LIEU OU IL Y A PLASTIFICATION';
  189. TRAC LLIS2 CHEL2 MO SUR1;
  190. FINSI;
  191. *
  192. SMRRA = EXTR CHPO1 SMRR PA; SMRRA1 = 0.1894;
  193. ERGRRA = ABS (100 * (SMRRA - SMRRA1) / SMRRA1) ;
  194. SMRRB = EXTR CHPO1 SMRR PB; SMRRB1 = 0.152 ;
  195. ERGRRB = ABS (100 * (SMRRB - SMRRB1) / SMRRB1) ;
  196. SMRRM = extr CHPO1 SMRR Pm; SMRRM1 = 47.09;
  197. ERGRRM = ABS (100 * (SMRRM - SMRRM1) / SMRRM1) ;
  198.  
  199.  
  200. SMZZA = EXTR CHPO1 SMZZ PA; SMZZA1 = 114.2 ;
  201. ERGZZA = ABS (100 * (SMZZA - SMZZA1) / SMZZA1) ;
  202. SMZZB = EXTR CHPO1 SMZZ PB; SMZZB1 = -345.9;
  203. ERGZZB = ABS (100 * (SMZZB - SMZZB1) / SMZZB1) ;
  204. SMZZm = EXTR CHPO1 SMZZ Pm; SMZZm1 = -285.1 ;
  205. ERGZZm = ABS (100 * (SMZZm - SMZZm1) / SMZZm1);
  206.  
  207. SMTTA = EXTR CHPO1 SMTT PA; SMTTA1 = 340.8 ;
  208. ERGTTA = ABS (100 * (SMTTA - SMTTA1) / SMTTA1) ;
  209. SMTTB = EXTR CHPO1 SMTT PB; SMTTB1 = -191. ;
  210. ERGTTB = ABS (100 * (SMTTB - SMTTB1) / SMTTB1) ;
  211. SMTTm = EXTR CHPO1 SMTT Pm; SMTTm1 = -51.52;
  212. ERGTTm = ABS (100 * (SMTTm - SMTTm1) / SMTTm1) ;
  213. opti echo 0;
  214.  
  215. MESS ' RESULTATS : TEST SSNA04/90 ';
  216. MESS ' -------------------------- ';
  217. MESS ' rA = 1; rM = 1.497; rB =2;';
  218.  
  219. MESS ' Resultats sur les contraintes interpolees (MPa)';
  220. MESS ' -----------------------------------------------';
  221. *
  222. MESS ' CONT. POINT REFERENCE CALCULEE ECART%';
  223. MESS ' SMRR A ' SMRRA1 SMRRA 'NON SIGNIFICATIF';
  224. MESS ' SMRR B ' SMRRB1 SMRRB 'NON SIGNIFICATIF';
  225. MESS ' SMRR M ' SMRRM1 SMRRM ERGRRM ;
  226.  
  227. MESS ' SMTT A ' SMTTA1 SMTTA ERGTTA ;
  228. MESS ' SMTT B ' SMTTB1 SMTTB ERGTTB ;
  229. MESS ' SMTT M ' SMTTm1 SMTTm ERGTTm ;
  230. *
  231. MESS ' SMZZ A ' SMZZA1 SMZZA ERGZZA ;
  232. MESS ' SMZZ B ' SMZZB1 SMZZB ERGZZB ;
  233. MESS ' SMZZ M ' SMZZm1 SMZZm ERGZZm ;
  234. SAUT 1 LIGN;
  235.  
  236. *
  237. DEP1 = TAB1.DEPLACEMENTS.NN;
  238. *
  239. SI (NEG GRAPH 'N');
  240. TITR ' DEFORMEE : CALCUL AXISYMETRIQUE (SSNA04)';
  241. DEF0 = DEFO SUR1 DEP1 0. ;
  242. DEF1 = DEFO SUR1 DEP1 20. ROUGE ;
  243. TRAC (DEF0 ET DEF1) ;
  244. FINSI;
  245. *
  246. URA = EXTR DEP1 UR PA;
  247. URB = EXTR DEP1 UR PB;
  248. URM = EXTR DEP1 UR Pm;
  249. *
  250. URA1 = 1.808E-03;
  251. URB1 = 3.859E-03;
  252. *
  253. ERGURA = ABS (100 * (URA - URA1) / URA1) ;
  254. ERGURB = ABS (100 * (URB - URB1) / URB1) ;
  255. *
  256. MESS ' Resultats sur les deplacements radiaux ';
  257. MESS ' -------------------------------------- ';
  258. *
  259. MESS ' Deplacement radial theorique pour r = 1 mm :' URA1 'mm';
  260. MESS ' Deplacement radial calcul pour r = 1 mm :' URA 'mm';
  261. MESS ' Soit une erreur relative de :' ERGURA '%';
  262.  
  263. MESS ' Deplacement radial calcul pour r = 1.497 mm : ' URM 'mm';
  264. *
  265. MESS ' Deplacement radial thorique pour r = 2 mm :' URB1 'mm';
  266. MESS ' Deplacement radial calcul pour r = 2 mm :' URB 'mm';
  267. MESS ' Soit une erreur relative de :' ERGURB '%';
  268. *
  269. *
  270. MESS ' Nombre de noeuds : ' NOSUR1;
  271. MESS ' Nombre d elements : ' ELSUR1;
  272. SAUT 1 LIGN ;
  273. *
  274. TEMPS;
  275. *
  276. SAUT 1 LIGN ;
  277. *
  278. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  279. * CODE FONCTIONNEMENT *
  280. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  281. *
  282. * Maximum des erreurs relatives
  283. *
  284. ERGMAX = MAXI (PROG ERGURA ERGURB ERGTTA ERGTTB ERGZZA ERGZZB
  285. ERGRRM ERGZZM ERGTTM);
  286. * Il y a erreur si un des pourcentage erg*** est superieur a 5%
  287. *
  288. SI (ERGMAX <EG 5.);
  289. ERRE 0 ;
  290. SINON;
  291. ERRE 5 ;
  292. FINSI;
  293. *
  294. FIN;
  295. *
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  

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