Télécharger ther71.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ther71.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. SAUT PAGE ;
  6. *-------------------------------------------------------------------*
  7. * TEST THER7 *
  8. * *
  9. * TEST DES CL DE TEMPERATURE POUR COQ8 ET COQ6 *
  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. *
  18. * *
  19. * UNE PLAQUE RECTANGULAIRE CONSTITUEE D'{L{ments COQ8 et coq6 est *
  20. * soumise @ une temp{rature impos{e @ une de ses extr{mit{s et suc- *
  21. * essivement @ une condition de convection forc{e, @ une condition *
  22. * de flux impos{ et @ une condition de source volumique impos{e. *
  23. * *
  24. * Les r{sultats sont compares a ceux obtenus *
  25. * par un calcul {quivalent utilisant les elements massifs cu20 *
  26. * On compare les temperatures des points situ{s sur la face externe *
  27. * du cylindre , le long de la ligne AS2 DS2. *
  28. * Test elabore par Mr Dowlatyari 03/91. *
  29. * Test modifie le 20/05/92. *
  30. *-------------------------------------------------------------------*
  31. opti echo 1;
  32. ********************************************************************
  33. **************** 1/ calcul en massif *******************************
  34. ********************************************************************
  35. OPTION DIME 3 ELEM QUA8 ;
  36. *tri6 ;
  37. opti elem cu20;
  38.  
  39. opti trac bens;
  40. graph = 'N';
  41.  
  42. *
  43. *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS --
  44. *
  45. As1= 0.25 0. 0.; AS2 =0.25 1.5 0. ; AS3 = 0.25 3. 0.;
  46.  
  47. BS1= 0.25 0. 2.; BS2 =0.25 1.5 2.; BS3 = 0.25 3. 2.;
  48.  
  49. CS1= 0.25 0. 4.; CS2 =0.25 1.5 4; CS3 = 0.25 3. 4.;
  50.  
  51. DS1= 0.25 0. 6.; DS2 =0.25 1.5 6.; DS3 = 0.25 3. 6.;
  52.  
  53.  
  54. *------------- CREATION DES SURFACES ET VOLUMES ------------------ -
  55. *
  56. DAD = AS3 d 3 DS3;
  57. DDD = DS3 D 2 DS1;
  58. DDA = DS1 D 3 AS1;
  59. DAA = AS1 D 2 AS3;
  60. QTOTS = DALL DAD DDD DDA DAA;
  61. V = -0.5 0. 0.;
  62. QTOT = QTOTS VOLU 1 TRAN V;
  63.  
  64. OEIL = 5. 1.5 3.;
  65.  
  66. ge0 = face 1 qtot;
  67. scon = elem ge0 appu larg ddd;
  68.  
  69. ge1 = face 3 qtot;
  70. dfa = qtot poin proc ( 0. 1.5 0. );
  71. sinf = elem ge1 appu larg dfa;
  72.  
  73.  
  74. ssou = elem qtot appu larg ddd;
  75.  
  76.  
  77. si ( NEG GRAPH 'N' );
  78. trac oeil qtot;
  79. trac oeil scon;
  80. trac oeil sinf;
  81. finsi;
  82.  
  83. *--------------- DONNEES DU PROBLEME DE THERMIQUE -----------------------
  84. *
  85. *-------------- MODELISATION -------------------------------------------
  86. *
  87. MOD1 = MODE QTOT THERMIQUE ISOTROPE ;
  88. MOD2 = MODE SCON thermique CONVECTION ;
  89. *
  90. *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU ---------------
  91. *
  92. MAT1 = MATE MOD1 'K' 100. ;
  93. MAT2 = MATE MOD2 'H' 15. ;
  94. *
  95. *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------
  96. *
  97. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  98. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  99. *
  100. *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI -------------------
  101. *
  102. BB1 = BLOQUE SINF 'TINF' ;
  103. EE1 = DEPI BB1 10. ;
  104. BB2 = BLOQUE SINF 'TSUP' ;
  105. EE2 = DEPI BB2 10. ;
  106. BB3 = BLOQUE SINF 'T ' ;
  107. EE3 = DEPI BB3 10. ;
  108. BB = BB1 ET BB2 ET BB3 ;
  109. EE = EE1 ET EE2 ET EE3 ;
  110. *
  111. *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ----------
  112. *
  113. CNV1 = CONVECTION MOD2 MAT2 'T' 200. ;
  114. *
  115. *--------------- FLUX IMPOSE --------------------------------------------
  116. *
  117. FLU1 = FLUX MOD1 200. SCON;
  118. *
  119. *--------------- SOURCE DE CHALEUR --------------------------------------
  120. *
  121. S1 = SOURCE MOD1 MAT1 1000. SSOU ;
  122. *
  123. *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE -----------------
  124. *
  125. CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ;
  126. *
  127. *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS -------------
  128. *
  129. FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ;
  130. *
  131. *---------------- 1ERE RESOLUTION : CONVECTION --------------------------
  132. *
  133. CHTER1 = RESOUDRE CCC1 FFF1 ;
  134. *
  135. *---------------- 2EME RESOLUTION : FLUX --------------------------------
  136. *
  137. CHTER2 = RESOUDRE CCC2 FFF2 ;
  138. *
  139. *---------------- 3EME RESOLUTION: SOURCE -------------------------------
  140. *
  141. CHTER3 = RESOUDRE CCC2 FFF3 ;
  142. *
  143. *---------------- POST-TRAITEMENT: -
  144. *---------------- EXTRACTION DES VALEURS CARACTERISTIQUES -----------
  145. *
  146. *LISTE CHTER1 ;
  147. *LISTE CHTER2 ;
  148. *LISTE CHTER3 ;
  149. *
  150. AS2 = QTOT POIN PROC ( 0.25 1.5 0.);
  151. BS2 = QTOT POIN PROC ( 0.25 1.5 2.);
  152. CS2 = QTOT POIN PROC ( 0.25 1.5 4.);
  153. DS2 = QTOT POIN PROC ( 0.25 1.5 6.);
  154.  
  155.  
  156. T1Aq = EXTR CHTER1 T AS2 ;T1Bq =EXTR CHTER1 T BS2 ;
  157. T1Cq = EXTR CHTER1 T CS2 ;T1Dq =EXTR CHTER1 T DS2 ;
  158. T2Aq = EXTR CHTER2 T AS2 ;T2Bq =EXTR CHTER2 T BS2 ;
  159. T2Cq = EXTR CHTER2 T CS2 ;T2Dq =EXTR CHTER2 T DS2 ;
  160. T3Aq = EXTR CHTER3 T AS2 ;T3Bq =EXTR CHTER3 T BS2 ;
  161. T3Cq = EXTR CHTER3 T CS2 ;T3Dq =EXTR CHTER3 T DS2;
  162.  
  163.  
  164. LCON = PROG T1Aq T1Bq T1Cq T1Dq;
  165. LFLU = PROG T2Aq T2Bq T2Cq T2Dq;
  166. LSOU = PROG T3Aq T3Bq T3Cq T3Dq;
  167. *LIST LCON; LIST LFLU; LIST LSOU;
  168.  
  169. *******************************************************************
  170. ********************* 2/ CALCUL EN COQ6 ET COQ8 *******************
  171. *******************************************************************
  172.  
  173. SAUT PAGE ;
  174.  
  175. TITRE 'C.L. DE TEMPERATURE - ELEMENTS COQ8 , COQ6 ' ;
  176. OPTION DIME 3 ELEM QUA8 ;
  177.  
  178. *
  179. *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS --
  180. *
  181. A1 = 0. 0. 0. ; B1 = 0. 3.0 0. ;AB1 =0. 1.5 0. ;
  182. A2= 0. 0. 1. ; B2 = 0. 3. 1. ; AB2 = 0. 1.5 2. ;
  183. A3 = 0. 0. 2. ; B3 = 0. 3. 2. ; AB4 = 0. 1.5 4. ;
  184. A4 = 0. 0. 3. ; B4 = 0. 3. 3. ; AB5 = 0. 1.5 6. ;
  185. A5 = 0. 0. 4. ;B5 = 0. 3. 4. ; AB3= 0. 1.5 3. ;
  186. A6 = 0. 0. 5. ;B6 = 0. 3. 5. ;
  187. A7 = 0. 0. 6. ;B7 = 0. 3. 6. ;
  188. *
  189. *------------- CREATION DES QUADRANGLES ET TRIANGLES -------------------
  190. *
  191. QH01 = MANU QUA8 A1 AB1 B1 B2 B3 AB2 A3 A2 ;
  192. QH02 = MANU TRI6 A3 AB2 B3 B4 B5 AB3 ;
  193. QH03 = MANU TRI6 A3 AB3 B5 AB4 A5 A4 ;
  194. QH04 = MANU QUA8 A5 AB4 B5 B6 B7 AB5 A7 A6 ;
  195. *
  196. *
  197. QTOT = QH01 ET QH02 ET QH03 ET QH04 ;
  198. OEIL = 100. 100. 100. ;
  199. SI ( NEG GRAPH 'N' );
  200. TRACE OEIL QTOT;
  201. FINSI;
  202. QSOUR = QH04 ;
  203. *
  204. *
  205. L_INF = MANU SEG3 A1 AB1 B1 ;
  206. *
  207. S_COT = QH04;
  208. *
  209. * ATTENTION : il faut r{orienter les {l{ments (OPERATEUR ORIENT)
  210. * si on n'est pas sur que tous les {l{ments sont orient{s d'une
  211. * facon coherante
  212. *
  213. *--------------- DONNEES DU PROBLEME DE THERMIQUE -----------------------
  214. *
  215. *-------------- MODELISATION -------------------------------------------
  216. *
  217. MOD1 = MODE QTOT THERMIQUE ISOTROPE COQ8 COQ6 ;
  218. MOD2 = MODE S_COT thermique CONVECTION SUPERIEURE COQ8 COQ6 ;
  219. *
  220. *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU ---------------
  221. *
  222. MAT1 = MATE MOD1 'K' 100. 'EPAI' 0.5 ;
  223. MAT2 = MATE MOD2 'H' 15. ;
  224. *
  225. *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------
  226. *
  227. *OPTION DONN 5 ;
  228. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  229.  
  230. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  231. *
  232. *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI -------------------
  233. *
  234. BB1 = BLOQUE L_INF 'TINF' ;
  235. EE1 = DEPI BB1 10. ;
  236. BB2 = BLOQUE L_INF 'TSUP' ;
  237. EE2 = DEPI BB2 10. ;
  238. BB3 = BLOQUE L_INF 'T ' ;
  239. EE3 = DEPI BB3 10. ;
  240. BB = BB1 ET BB2 ;
  241. EE = EE1 ET EE2 ;
  242. *
  243. *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ----------
  244. *
  245. CNV1 = CONVECTION MOD2 MAT2 'TSUP' 200. ;
  246. *
  247. *--------------- FLUX IMPOSE --------------------------------------------
  248. *
  249. FLU1 = FLUX MOD1 200. S_COT SUPE ;
  250. *
  251. *--------------- SOURCE DE CHALEUR --------------------------------------
  252. *
  253. S1 = SOURCE MOD1 MAT1 1000. QSOUR ;
  254. *
  255. *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE -----------------
  256. *
  257. CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ;
  258. *
  259. *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS -------------
  260. *
  261. FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ;
  262. *FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;
  263. *
  264. *---------------- 1ERE RESOLUTION : CONVECTION --------------------------
  265. *
  266. *OPTION DONN 5 ;
  267. CHTER1 = RESOUDRE CCC1 FFF1 ;
  268. *
  269. *---------------- 2EME RESOLUTION : FLUX --------------------------------
  270. *
  271. CHTER2 = RESOUDRE CCC2 FFF2 ;
  272. *
  273. *---------------- 3EME RESOLUTION: SOURCE -------------------------------
  274. *
  275. CHTER3 = RESOUDRE CCC2 FFF3 ;
  276. *
  277. ************************************************************************-
  278. ****************** 3/COMPARAISON DES 2 CALCULS *************************
  279. ************************************************************************
  280. *LISTE CHTER1 ;
  281. *LISTE CHTER2 ;
  282. *LISTE CHTER3 ;
  283. *
  284. OEIL= 1000. 0. 0. ;
  285. *
  286. SI ( NEG GRAPH 'N');
  287.  
  288. TITRE 'ELEMENTS COQ6,COQ8 : CONVECTION' ;
  289. CHTER1I = 'EXCO' CHTER1 'TINF' 'NOID' 'TINF' ;
  290. TRACER OEIL QTOT CHTER1I ;
  291. TITRE 'ELEMENTS COQ6,COQ8 : FLUX IMPOSE' ;
  292. CHTER2I = 'EXCO' CHTER2 'TINF' 'NOID' 'TINF' ;
  293. TRACER OEIL QTOT CHTER2I ;
  294. TITRE 'ELEMENTS COQ6,COQ8 : SOURCE' ;
  295. CHTER3I = 'EXCO' CHTER3 'TINF' 'NOID' 'TINF' ;
  296. TRACER OEIL QTOT CHTER3I ;
  297. FINSI;
  298.  
  299. T1AB1 = EXTR CHTER1 TSUP AB1 ;T1AB2=EXTR CHTER1 TSUP AB2 ;
  300. T1AB3 = EXTR CHTER1 TSUP AB3 ;T1AB4=EXTR CHTER1 TSUP AB4 ;
  301. T1AB5 = EXTR CHTER1 TSUP AB5 ;
  302. T2AB1= EXTR CHTER2 TSUP AB1 ;T2AB2=EXTR CHTER2 TSUP AB2 ;
  303. T2AB3 = EXTR CHTER2 TSUP AB3 ;T2AB4=EXTR CHTER2 TSUP AB4 ;
  304. T2AB5 = EXTR CHTER2 TSUP AB5 ;
  305. T3AB1 = EXTR CHTER3 TSUP AB1 ;T3AB2=EXTR CHTER3 TSUP AB2 ;
  306. T3AB3 = EXTR CHTER3 TSUP AB3 ;T3AB4=EXTR CHTER3 TSUP AB4 ;
  307. T3AB5 = EXTR CHTER3 TSUP AB5 ;
  308. *
  309. si ( NEG GRAPH 'N');
  310.  
  311. DIST = PROG 0. 2. 4. 6. ;
  312. DIST1 = PROG 0. 2. 3. 4. 6. ;
  313. TEMPR1=PROG 10. 71.89 134.62 163.20 ;
  314. TEMPC1= PROG T1AB1 T1AB2 T1AB3 T1AB4 T1AB5 ;
  315. EV1R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' lcon ;
  316. EV1C = EVOL MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC1 ;
  317. TAB = 'TABLE' ;
  318. TAB.1 = 'MARQ PLUS' ;
  319. TAB.2 = 'MARQ ETOI' ;
  320. TITRE 'COQ6 ET COQ8 ;TEMPERATURE + CONVECTION ;+ : MASSIF,* : COQUE' ;
  321. DESSIN TAB (EV1R ET EV1C ) ;
  322. TEMPR2=PROG 10. 26.16 42.49 50.79 ;
  323. TEMPC2= PROG T2AB1 T2AB2 T2AB3 T2AB4 T2AB5 ;
  324. EV2R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' lflu;
  325. EV2C = EVOL MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC2 ;
  326. TITRE 'COQ6 ET COQ8;TEMPERATURE + FLUX IMPOSE ;+ : MASSIF,* : COQUE';
  327. DESSIN TAB (EV2R ET EV2C ) ;
  328. TEMPR3=PROG 10. 50. 90. 110. ;
  329. TEMPC3= PROG T3AB1 T3AB2 T3AB3 T3AB4 T3AB5 ;
  330. EV3R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' lsou;
  331. EV3C = EVOL MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC3 ;
  332. TITRE 'COQ6 ET COQ8 ;TEMPERATURE + SOURCE ;+ MASSIF,* : COQUE ';
  333. DESSIN TAB (EV3R ET EV3C ) ;
  334. FINSI;
  335.  
  336. *************************************************************
  337. *************** 4/ TEST DE FONCTIONNEMENT *******************
  338. *************************************************************
  339.  
  340. ERCA= ABS (( T1AB1-T1Aq)/T1AB1);
  341. ERCB= ABS (( T1AB2-T1Bq)/T1AB2);
  342. ERCC= ABS (( T1AB4-T1Cq)/T1AB4);
  343. ERCD= ABS (( T1AB5-T1Dq)/T1AB5);
  344. ERFA= ABS (( T2AB1-T2Aq)/T2AB1);
  345. ERFB= ABS (( T2AB2-T2Bq)/T2AB2);
  346. ERFC= ABS (( T2AB4-T2Cq)/T2AB4);
  347. ERFD= ABS (( T2AB5-T2Dq)/T2AB5);
  348. ERSA= ABS (( T3AB1-T3Aq)/T3AB1);
  349. ERSB= ABS (( T3AB2-T3Bq)/T3AB2);
  350. ERSC= ABS (( T3AB4-T3Cq)/T3AB4);
  351. ERSD= ABS (( T3AB5-T3Dq)/T3AB5);
  352.  
  353. ERTOT = PROG ERCA ERCB ERCC ERCD ERFA ERFB ERFC ERFD ERSA ERSB
  354. ERSC ERSD;
  355. ermax = maxi ertot;
  356. temps;
  357. SI ( ERMAX <EG 0.05 );
  358. ERRE 0;
  359. SINON;
  360. ERRE 5;
  361. FINSI;
  362. fin;
  363.  
  364.  
  365.  
  366. FIN ;
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  

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