Télécharger ther51.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ther51.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. SAUT PAGE ;
  6. *-------------------------------------------------------------------*
  7. * TEST THER51 *
  8. * *
  9. * TEST DES CL DE TEMPERATURE POUR L'ELEMENT COQ2 *
  10. * ------------------------------------------------------------ *
  11. * *
  12. * TEMPERATURE IMPOSEE + CONVECTION *
  13. * + FLUX *
  14. * + SOURCE *
  15. * *
  16. * Ce test permet de v{rifier le bon fonctionnement des divers *
  17. * op{rateurs thermiques de CASTEM2000 dans le cas des {l{ments *
  18. * coq2 . *
  19. * *
  20. * Un cylindre modelis{ par des {L{ments coq2 est soumis @ une *
  21. * temp{rature impos{e @ une de ses extr{mit{s et successivement *
  22. * @ une condition de convection forc{e, @ une condition de flux *
  23. * impos{ et @ une condition de source volumique impos{e. *
  24. * *
  25. * Les r{sultats sont compares a ceux obtenus *
  26. * par un calcul {quivalent utilisant les {l{ments massifs. En *
  27. * particulier on compare la temperature des points situ{s a la face *
  28. * exterieure du cylindre . *
  29. * Test elabore par Mr Dowlatyari 03/91, *
  30. * Test modifie le 13/05/92 *
  31. *-------------------------------------------------------------------*
  32. opti echo 1;
  33. graph = 'N';
  34. opti trac bens;
  35. SAUT PAGE ;
  36. ********************************************************************** -
  37. *************** 1/CALCUL EN ELEMENTS MASSIFS ************************
  38. **********************************************************************
  39.  
  40.  
  41. OPTION DIME 2 ELEM QUA4 MODE AXISYMETRIQUE ;
  42. TEMPS ;
  43. *
  44. *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS --
  45. *
  46. A1 = 19.75 0. ; A2 = 20.25 0.;
  47. B1 = 19.75 2. ; B2 = 20.25 2.;
  48. C1 = 19.75 4. ; C2 = 20.25 4.;
  49. D1 = 19.75 6. ; D2 = 20.25 6.;
  50. *
  51. *------------- CREATION DES SEGMENTS ET RECTANGLES --------------------
  52. *
  53. L1 = A2 D 1 B2;
  54. L2 = B2 D 1 C2;
  55. L3 = C2 D 1 D2;
  56. *
  57. LTOT = L1 ET L2 ET L3 ;
  58. LINF = A1 D 1 A2;
  59. RECSUP = MANU QUA4 C1 C2 D2 D1;
  60. RECMI = MANU QUA4 B1 B2 C2 C1;
  61. RECINF = MANU QUA4 A1 A2 B2 B1;
  62. *
  63. LTOT = RECSUP ET RECMI ET RECINF;
  64. si (NEG GRAPH 'N');
  65. TRAC LTOT;
  66. finsi;
  67. L_COT= L3 ;
  68. * ATTENTION : il faut r{orienter les {l{ments (OPERATEUR ORIENT)
  69. * si on n'est pas sºr que tous les {l{ments sont orient{s d'une
  70. * fa\on coh{rante
  71. *
  72. *--------------- DONNEES DU PROBLEME DE THERMIQUE -----------------------
  73. *
  74. *-------------- MODELISATION -------------------------------------------
  75. *
  76. MOD1 = MODE LTOT THERMIQUE ISOTROPE ;
  77. MOD2 = MODE L_COT thermique CONVECTION ;
  78. *
  79. *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU ---------------
  80. *
  81. MAT1 = MATE MOD1 'K' 100. ;
  82. MAT2 = MATE MOD2 'H' 15. ;
  83. *
  84. *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------
  85. *
  86. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  87. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  88. *
  89. *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI -------------------
  90. *
  91. BB1 = BLOQUE LINF 'TINF' ;
  92. EE1 = DEPI BB1 10. ;
  93. BB2 = BLOQUE LINF 'TSUP' ;
  94. EE2 = DEPI BB2 10. ;
  95. BB3 = BLOQUE LINF 'T ' ;
  96. EE3 = DEPI BB3 10. ;
  97. BB=BB1 ET BB2 ET BB3 ;
  98. EE=EE1 ET EE2 ET EE3 ;
  99. *
  100. *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ----------
  101. *
  102. CNV1 = CONVECTION MOD2 MAT2 'T' 200. ;
  103. *
  104. *--------------- FLUX IMPOSE --------------------------------------------
  105. *
  106. FLU1 = FLUX MOD1 200. L_COT ;
  107. *
  108. *--------------- SOURCE DE CHALEUR --------------------------------------
  109. *
  110. S1 = SOURCE MOD1 MAT1 1000. RECSUP;
  111. *
  112. *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE -----------------
  113. *
  114. CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ;
  115. *
  116. *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS -------------
  117. *
  118. FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ;
  119. *
  120. *---------------- 1ERE RESOLUTION : CONVECTION --------------------------
  121. *
  122. CHTER1 = RESOUDRE CCC1 FFF1 ;
  123. *
  124. *---------------- 2EME RESOLUTION : FLUX --------------------------------
  125. *
  126. CHTER2 = RESOUDRE CCC2 FFF2 ;
  127. *
  128. *---------------- 3EME RESOLUTION: SOURCE -------------------------------
  129. *
  130. CHTER3 = RESOUDRE CCC2 FFF3 ;
  131. *
  132. *---------------- POST-TRAITEMENT:
  133. *--------------- EXTRACTION DES VALEURS CARACTERISTIQUES ------------
  134. *-------------ET TRACE DES COURBES COMPARATIVES------------------------
  135. *
  136.  
  137.  
  138. *LISTE CHTER1 ;
  139. *LISTE CHTER2 ;
  140. *LISTE CHTER3 ;
  141. *T1A1 = EXTR CHTER1 T A1 ;T1B1 =EXTR CHTER1 T B1 ;
  142. *T1C1 = EXTR CHTER1 T C1 ;T1D1 =EXTR CHTER1 T D1 ;
  143. *T2A1 = EXTR CHTER2 T A1 ;T2B1 =EXTR CHTER2 T B1 ;
  144. *T2C1 = EXTR CHTER2 T C1 ;T2D1 =EXTR CHTER2 T D1 ;
  145. *T3A1 = EXTR CHTER3 T A1 ;T3B1 =EXTR CHTER3 T B1 ;
  146. *T3C1 = EXTR CHTER3 T C1 ;T3D1 =EXTR CHTER3 T D1 ;
  147.  
  148. T1A2 = EXTR CHTER1 T A2 ;T1B2 =EXTR CHTER1 T B2 ;
  149. T1C2 = EXTR CHTER1 T C2 ;T1D2 =EXTR CHTER1 T D2 ;
  150. T2A2 = EXTR CHTER2 T A2 ;T2B2 =EXTR CHTER2 T B2 ;
  151. T2C2 = EXTR CHTER2 T C2 ;T2D2 =EXTR CHTER2 T D2 ;
  152. T3A2 = EXTR CHTER3 T A2 ;T3B2 =EXTR CHTER3 T B2 ;
  153. T3C2 = EXTR CHTER3 T C2 ;T3D2 =EXTR CHTER3 T D2 ;
  154.  
  155. *T1AM = ( T1A2 + T1A1)/2;T1BM = ( T1B2 + T1B1)/2;
  156. *T1CM = ( T1C2 + T1C1)/2;T1DM = ( T1D2 + T1D1)/2;
  157. *T2AM = ( T2A2 + T2A1)/2;T2BM = ( T2B2 + T2B1)/2;
  158. *T2CM = ( T2C2 + T2C1)/2;T2DM = ( T2D2 + T2D1)/2;
  159. *T3AM = ( T3A2 + T3A1)/2;T3BM = ( T3B2 + T3B1)/2;
  160. *T3CM = ( T3C2 + T3C1)/2;T3DM = ( T3D2 + T3D1)/2;
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167. TCON = PROG T1A2 T1B2 T1C2 T1D2;
  168. *LIST TCON;
  169. TFLUX = PROG T2A2 T2B2 T2C2 T2D2;
  170. *LIST TFLUX;
  171. TSOU = PROG T3A2 T3B2 T3C2 T3D2;
  172. *LIST TSOU;
  173.  
  174. *********************************************************************
  175. ***************** 2/ CALCUL EN COQUES ******************************
  176. *********************************************************************
  177.  
  178.  
  179. *
  180. TITRE 'C.L. DE TEMPERATURE - ELEMENTS COQ2' ;
  181. OPTION DIME 2 ELEM SEG2 MODE AXISYMETRIQUE ;
  182. TEMPS ;
  183. *
  184. *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS --
  185. *
  186. A1 = 20. 0. ;
  187. B1 = 20. 2. ;
  188. C1 = 20. 4. ;
  189. D1 = 20. 6. ;
  190. *
  191. *------------- CREATION DES SEGMENTS ---------------------------------
  192. *
  193. L1= MANU SEG2 A1 B1 ;
  194. L2 = MANU SEG2 B1 C1;
  195. L3 = MANU SEG2 C1 D1 ;
  196. *
  197. LTOT = L1 ET L2 ET L3 ;
  198.  
  199. *
  200. *
  201. P_INF = A1 ;
  202. *
  203. L_COT= L3 ;
  204. * ATTENTION : il faut r{orienter les {l{ments (OPERATEUR ORIENT)
  205. * si on n'est pas sºr que tous les {l{ments sont orient{s d'une
  206. * fa\on coh{rante
  207. *
  208. *--------------- DONNEES DU PROBLEME DE THERMIQUE -----------------------
  209. *
  210. *-------------- MODELISATION -------------------------------------------
  211. *
  212. MOD1 = MODE LTOT THERMIQUE ISOTROPE COQ2 ;
  213. MOD2 = MODE L_COT thermique CONVECTION 'SUPERIEURE' COQ2 ;
  214. *
  215. *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU ---------------
  216. *
  217. MAT1 = MATE MOD1 'K' 100. 'EPAI' 0.5 ;
  218. MAT2 = MATE MOD2 'H' 15. ;
  219. *
  220. *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------
  221. *
  222. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  223. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  224. *
  225. *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI -------------------
  226. *
  227. BB1 = BLOQUE P_INF 'TINF' ;
  228. EE1 = DEPI BB1 10. ;
  229. BB2 = BLOQUE P_INF 'TSUP' ;
  230. EE2 = DEPI BB2 10. ;
  231. BB3 = BLOQUE P_INF 'T ' ;
  232. EE3 = DEPI BB3 10. ;
  233. BB=BB1 ET BB2 ET BB3 ;
  234. EE=EE1 ET EE2 ET EE3 ;
  235. *
  236. *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ----------
  237. *
  238. CNV1 = CONVECTION MOD2 MAT2 'TSUP' 200. ;
  239. *
  240. *--------------- FLUX IMPOSE --------------------------------------------
  241. *
  242. FLU1 = FLUX MOD1 200. L_COT 'SUPE' ;
  243. *
  244. *--------------- SOURCE DE CHALEUR --------------------------------------
  245. *
  246. S1 = SOURCE MOD1 MAT1 1000. L3 ;
  247. *
  248. *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE -----------------
  249. *
  250. CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ;
  251. *
  252. *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS -------------
  253. *
  254. FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ;
  255. *
  256. *---------------- 1ERE RESOLUTION : CONVECTION --------------------------
  257. *
  258. CHTER1 = RESOUDRE CCC1 FFF1 ;
  259. *
  260. *---------------- 2EME RESOLUTION : FLUX --------------------------------
  261. *
  262. CHTER2 = RESOUDRE CCC2 FFF2 ;
  263. *
  264. *---------------- 3EME RESOLUTION: SOURCE -------------------------------
  265. *
  266. CHTER3 = RESOUDRE CCC2 FFF3 ;
  267. *
  268. *---------------- POST-TRAITEMENT:
  269. *--------------- EXTRACTION DES VALEURS CARACTERISTIQUES ------------
  270. *-------------ET TRACE DES COURBES COMPARATIVES------------------------
  271. *
  272. *
  273. *ISTE CHTER1 ;
  274. *LISTE CHTER2 ;
  275. *LISTE CHTER3 ;
  276. T1A = EXTR CHTER1 T A1 ;T1B=EXTR CHTER1 T B1 ;
  277. T1C = EXTR CHTER1 T C1 ;T1D=EXTR CHTER1 T D1 ;
  278. T2A = EXTR CHTER2 T A1 ;T2B=EXTR CHTER2 T B1 ;
  279. T2C = EXTR CHTER2 T C1 ;T2D=EXTR CHTER2 T D1 ;
  280. T3A = EXTR CHTER3 T A1 ;T3B=EXTR CHTER3 T B1 ;
  281. T3C = EXTR CHTER3 T C1 ;T3D=EXTR CHTER3 T D1 ;
  282.  
  283.  
  284. ******************************************************************
  285. ************** 3/ COMPARAISON ENTRE LES DEUX CALCULS *************
  286. ******************************************************************
  287.  
  288.  
  289.  
  290. DIST = PROG 0. 2. 4. 6. ;
  291.  
  292. si (NEG GRAPH 'N');
  293. TEMPC1= PROG T1A T1B T1C T1D ;
  294. EV1R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TCON ;
  295. EV1C = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPC1 ;
  296. TAB = 'TABLE' ;
  297. TAB.1 = 'MARQ PLUS' ;
  298. TAB.2 = 'MARQ ETOI' ;
  299. TITRE 'ELEMENT COQ2 ;TEMPERATURE + CONVECTION ;+ : MASSIF,* : COQUE' ;
  300. DESSIN TAB (EV1R ET EV1C ) ;
  301.  
  302. TEMPC2= PROG T2A T2B T2C T2D ;
  303. EV2R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TFLUX;
  304. EV2C = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPC2 ;
  305. TITRE 'ELEMENT COQ2 ;TEMPERATURE + FLUX IMPOSE ;+ : MASSIF,* : COQUE';
  306. DESSIN TAB (EV2R ET EV2C ) ;
  307.  
  308. TEMPC3= PROG T3A T3B T3C T3D ;
  309. EV3R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TSOU;
  310. EV3C = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPC3 ;
  311. TITRE 'ELEMENT COQ2 ;TEMPERATURE + SOURCE ;+ MASSIF,* : COQUE ';
  312. DESSIN TAB (EV3R ET EV3C ) ;
  313. finsi;
  314. ******************************************************************
  315. ***************** 4/ TEST DE FONCTIONNEMENT **********************
  316. ******************************************************************
  317.  
  318. ERCA= ABS (( T1A-T1A2)/T1A);
  319. ERCB= ABS (( T1B-T1B2)/T1B);
  320. ERCC= ABS (( T1C-T1C2)/T1C);
  321. ERCD= ABS (( T1D-T1D2)/T1D);
  322. ERFA= ABS (( T2A-T2A2)/T2A);
  323. ERFB= ABS (( T2B-T2B2)/T2B);
  324. ERFC= ABS (( T2C-T2C2)/T2C);
  325. ERFD= ABS (( T2D-T2D2)/T2D);
  326. ERSA= ABS (( T3A-T3A2)/T3A);
  327. ERSB= ABS (( T3B-T3B2)/T3B);
  328. ERSC= ABS (( T3C-T3C2)/T3C);
  329. ERSD= ABS (( T3D-T3D2)/T3D);
  330.  
  331. ERTOT = PROG ERCA ERCB ERCC ERCD ERFA ERFB ERFC ERFD ERSA ERSB
  332. ERSC ERSD;
  333. ERMAX = MAXI ERTOT;
  334. temps;
  335. SI ( ERMAX <EG 0.05 );
  336. ERRE 0;
  337. SINON;
  338. ERRE 5;
  339. FINSI;
  340.  
  341. FIN;
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  
  356.  
  357.  

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