Télécharger ther4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ther4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************
  6. * Test ther4.dgibi: jeux de données *
  7. * --------------------------------- *
  8. *******************************************
  9.  
  10. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  11. * SINON SI GRAPH DIFFERENT DE N TOUS
  12. * LES GRAPHIQUES SONT AFFICHES
  13.  
  14. GRAPH = 'N' ;
  15.  
  16. SAUT PAGE;
  17. SI (NEG GRAPH 'N') ;
  18. OPTI ECHO 1 ;
  19. OPTI TRAC PSC ;
  20. SINO ;
  21. OPTI ECHO 0 ;
  22. FINSI ;
  23.  
  24. SAUT PAGE;
  25. *-----------------------------------------*
  26. * THER4 *
  27. * *
  28. * TEST DES CL DE TEMPERATURE POUR *
  29. * PY13 CU20 TE10 PR15 TRIDIM: *
  30. * --------------------------------------*
  31. * *
  32. * TEMPERATURE IMPOSEE + CONVECTION *
  33. * + FLUX *
  34. * + SOURCE *
  35. * *
  36. * Cet exemple permet de tester les *
  37. * conditions aux limites de température *
  38. * pour des élements massifs 3D quadratiques
  39. * *
  40. * Un parallélépipède rectangle constitué *
  41. * d'éléments PY13, CU20, TE10 et PR15 est *
  42. * soumis à une température imposée à sa *
  43. * face supérieure et successivement à une *
  44. * convection forcée, à un flux imposé et à*
  45. * une source volumique de chaleur. *
  46. * *
  47. * Les résultats sont présentés sous forme *
  48. * d'isovaleurs qui doivent etre parallèles*
  49. * aux faces supérieure et inférieure. *
  50. * *
  51. * 21 DECEMBRE 1988 *
  52. *-----------------------------------------*
  53.  
  54. *--- DISPOSITION DES NOEUDS DES ELEMENTS:
  55. *
  56. * A4__AF4___F4___DF4___D4
  57. * / / /|
  58. * AB4 EF4 CD4 D34
  59. * / / / |
  60. * B4___BE4___E4___CE4___C4 D3
  61. * | 6PY13 | 1CU20 | /
  62. * B34 E34 C34 CD3
  63. * | | |/
  64. * B3___BE3___E3___CE3___C3
  65. *
  66. * IN:CENTRE DE
  67. * (AN BN EN FN AN+1 BN+1 EN+1 FN+1)
  68. *
  69. * JN:CENTRE DE
  70. * (DN CN EN FN DN+1 CN+1 EN+1 FN+1)
  71. *
  72. * A3__AF3___F3___DF3___D3
  73. * / / /|
  74. * AB3 EF3 DE3 CD3 D23
  75. * / / / |
  76. * B3___BE3___E3___CE3___C3 D2
  77. * | 1CU20 | 4PY13 | /
  78. * B23 E23 + C23 CD2
  79. * | | 4TE10 |/
  80. * B2___BE2___E2___CE2___C2
  81. *
  82. * A2__AF2___F2___DF2___D2
  83. * / / /|
  84. * AB2 EF2 DE2 CD2 D12
  85. * / / / |
  86. * B2___BE2___E2___CE2___C2 D1
  87. * | 6PY13 | 2PR15 | /
  88. * B12 E12 C12 CD1
  89. * | | |/
  90. * B1___BE1___E1___CE1___C1
  91. *
  92. *------------------------------------------
  93. *--- OPTIONS GENERALES DE CALCUL:
  94. *
  95. TITR 'THER4 - C.L. TEMPERATURE -
  96. ELEMENTS QUADRATIQUES' ;
  97. TEMPS ;
  98.  
  99. *--- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS
  100. A1=0. 0. 0.; B1=1. 0. 0.;
  101. C1=1. 4. 0.; D1=0. 4. 0.;
  102. A2=0. 0. 2.; B2=1. 0. 2.;
  103. C2=1. 4. 2.; D2=0. 4. 2.;
  104. A3=0. 0. 4.; B3=1. 0. 4.;
  105. C3=1. 4. 4.; D3=0. 4. 4.;
  106. A4=0. 0. 6.; B4=1. 0. 6.;
  107. C4=1. 4. 6.; D4=0. 4. 6.;
  108. E1=1. 2. 0.; F1=0. 2. 0.;
  109. E2=1. 2. 2.; F2=0. 2. 2.;
  110. E3=1. 2. 4.; F3=0. 2. 4.;
  111. E4=1. 2. 6.; F4=0. 2. 6.;
  112. I1= .5 1. 1.; J2=.5 3. 3.;
  113. I3= .5 1. 5.;
  114.  
  115. A12=0. 0. 1.; B12=1. 0. 1.;
  116. C12=1. 4. 1.; D12=0. 4. 1.;
  117. A23=0. 0. 3.; B23=1. 0. 3.;
  118. C23=1. 4. 3.; D23=0 4. 3.;
  119. A34=0. 0. 5.; B34=1. 0. 5.;
  120. C34=1. 4. 5.; D34=0. 4. 5.;
  121. E12=1. 2. 1.; F12=0. 2. 1.;
  122. E23=1. 2. 3.; F23=0. 2. 3.;
  123. E34=1. 2. 5.; F34=0. 2. 5.;
  124.  
  125. AB1=.5 0. 0.; AB2=.5 0. 2.;
  126. AB3=.5 0. 4.; AB4=.5 0. 6.;
  127. BE1=1. 1. 0.; BE2=1. 1. 2.;
  128. BE3=1. 1. 4.; BE4=1. 1. 6.;
  129. CE1=1. 3. 0.; CE2=1. 3. 2.;
  130. CE3=1. 3. 4.; CE4=1. 3. 6.;
  131. CD1=.5 4. 0.; CD2=.5 4. 2.;
  132. CD3=.5 4. 4.; CD4=.5 4. 6.;
  133. DF1=0. 3. 0.; DF2=0. 3. 2.;
  134. DF3=0. 3. 4.; DF4=0. 3. 6.;
  135. AF1=0. 1. 0.; AF2=0. 1. 2.;
  136. AF3=0. 1. 4.; AF4=0. 1. 6.;
  137. EF1=.5 2. 0.; EF2=.5 2. 2.;
  138. EF3=.5 2. 4.; EF4=.5 2. 6.;
  139. MEC1=1. 3. 1.; MFD1=0. 3. 1.;
  140.  
  141. *--- MILIEUX DES SEGMENTS (XN IN) ET (XN JN)
  142. AI1= .25 .5 .5; BI1= .75 .5 .5;
  143. EI1= .75 1.5 .5; FI1= .25 1.5 .5;
  144. EJ2= .75 2.5 2.5; FJ2= .25 2.5 2.5;
  145. DJ2= .25 3.5 2.5; CJ2= .75 3.5 2.5;
  146. AI3= .25 .5 4.5; BI3= .75 .5 4.5;
  147. EI3= .75 1.5 4.5; FI3= .25 1.5 4.5;
  148.  
  149. *--- MILIEUX DES SEGMENTS (XN+1 IN) ET (XN+1 JN)
  150. A2I1=.25 .5 1.5; B2I1=.75 .5 1.5;
  151. E2I1=.75 1.5 1.5; F2I1=.25 1.5 1.5;
  152. E3J2=.75 2.5 3.5; F3J2=.25 2.5 3.5;
  153. D3J2=.25 3.5 3.5; C3J2=.75 3.5 3.5;
  154. A4I3=.25 .5 5.5; B4I3=.75 .5 5.5;
  155. E4I3=.75 1.5 5.5; F4I3=.25 1.5 5.5;
  156. DE2 =.5 3. 2. ; DE3 =.5 3. 4. ;
  157.  
  158. *--- CREATION DES PYRAMIDES A 13 NOEUDS:
  159. PP01 = MANU PY13
  160. A1 AB1 B1 BE1 E1 EF1 F1 AF1
  161. AI1 BI1 EI1 FI1 I1 ;
  162. PP02 = MANU PY13
  163. E1 EF1 F1 F12 F2 EF2 E2 E12
  164. EI1 FI1 F2I1 E2I1 I1 ;
  165. PP03 = MANU PY13
  166. E2 EF2 F2 AF2 A2 AB2 B2 BE2
  167. E2I1 F2I1 A2I1 B2I1 I1 ;
  168. PP04 = MANU PY13
  169. A2 AB2 B2 B12 B1 AB1 A1 A12
  170. A2I1 B2I1 BI1 AI1 I1 ;
  171. PP05 = MANU PY13
  172. B1 BE1 E1 E12 E2 BE2 B2 B12
  173. BI1 EI1 E2I1 B2I1 I1 ;
  174. PP06 = MANU PY13
  175. A1 AF1 F1 F12 F2 AF2 A2 A12
  176. AI1 FI1 F2I1 A2I1 I1 ;
  177. PP07 = MANU PY13
  178. E2 EF2 F2 F23 F3 EF3 E3 E23
  179. EJ2 FJ2 F3J2 E3J2 J2 ;
  180. PP08 = MANU PY13
  181. D3 CD3 C3 C23 C2 CD2 D2 D23
  182. D3J2 C3J2 CJ2 DJ2 J2 ;
  183. PP09 = MANU PY13
  184. C2 CE2 E2 E23 E3 CE3 C3 C23
  185. CJ2 EJ2 E3J2 C3J2 J2 ;
  186. PP10 = MANU PY13
  187. D2 DF2 F2 F23 F3 DF3 D3 D23
  188. DJ2 FJ2 F3J2 D3J2 J2 ;
  189. PP11 = MANU PY13
  190. A3 AB3 B3 BE3 E3 EF3 F3 AF3
  191. AI3 BI3 EI3 FI3 I3 ;
  192. PP12 = MANU PY13
  193. E3 EF3 F3 F34 F4 EF4 E4 E34
  194. EI3 FI3 F4I3 E4I3 I3 ;
  195. PP13 = MANU PY13
  196. E4 EF4 F4 AF4 A4 AB4 B4 BE4
  197. E4I3 F4I3 A4I3 B4I3 I3 ;
  198. PP14 = MANU PY13
  199. A4 AB4 B4 B34 B3 AB3 A3 A34
  200. A4I3 B4I3 BI3 AI3 I3 ;
  201. PP15 = MANU PY13
  202. B3 BE3 E3 E34 E4 BE4 B4 B34
  203. BI3 EI3 E4I3 B4I3 I3 ;
  204. PP16 = MANU PY13
  205. A3 AF3 F3 F34 F4 AF4 A4 A34
  206. AI3 FI3 F4I3 A4I3 I3 ;
  207.  
  208. *--- VOLUME COMPOSE DES 16 PYRAMIDES
  209. * A 13 NOEUDS:
  210. VOL1 = PP01 ET PP02 ET PP03 ET PP04 ET
  211. PP05 ET PP06 ET PP07 ET PP08 ET
  212. PP09 ET PP10 ET PP11 ET PP12 ET
  213. PP13 ET PP14 ET PP15 ET PP16 ;
  214.  
  215. *--- CREATION DES TETRAEDRES
  216. * A 10 NOEUDS:
  217. TT01 = MANU TE10
  218. D2 CD2 C2 CE2 E2 DE2 DJ2 CJ2
  219. EJ2 J2 ;
  220. TT02 = MANU TE10
  221. D2 DF2 F2 EF2 E2 DE2 DJ2 FJ2
  222. EJ2 J2 ;
  223. TT03 = MANU TE10
  224. D3 CD3 C3 CE3 E3 DE3 D3J2 C3J2
  225. E3J2 J2 ;
  226. TT04 = MANU TE10
  227. D3 DF3 F3 EF3 E3 DE3 D3J2 F3J2
  228. E3J2 J2 ;
  229.  
  230. *--- VOLUME COMPOSE DES 4 TETRAEDRES
  231. * A 10 NOEUDS:
  232. VOL2 = TT01 ET TT02 ET TT03 ET TT04 ;
  233.  
  234. *--- CREATION DES PRISMES A 15 NOEUDS:
  235. PM01 = MANU PR15 E2 E12 E1 MEC1 C2 CE2 EF2
  236. EF1 CD2
  237. F2 F12 F1 MFD1 D2 DF2 ;
  238. PM02 = MANU PR15 E1 CE1 C1 C12 C2 MEC1 EF1
  239. CD1 CD2
  240. F1 DF1 D1 D12 D2 MFD1 ;
  241.  
  242. *--- VOLUME COMPOSE DES 2 PRISMES
  243. * A 15 NOEUDS:
  244. VOL4 = PM01 ET PM02 ;
  245.  
  246. *--- CREATION DES CUBES A 20 NOEUDS:
  247. CB02 = MANU CU20
  248. E2 BE2 B2 AB2 A2 AF2 F2 EF2 E23 B23 A23 F23
  249. E3 BE3 B3 AB3 A3 AF3 F3 EF3 ;
  250. CB03 = MANU CU20
  251. E3 CE3 C3 CD3 D3 DF3 F3 EF3 E34 C34 D34 F34
  252. E4 CE4 C4 CD4 D4 DF4 F4 EF4 ;
  253.  
  254. *--- VOLUME COMPOSE DES 2 CUBES
  255. * A 20 NOEUDS:
  256. VOL3 = CB02 ET CB03 ;
  257.  
  258. *--- REUNION DES VOLUMES ELEMENTAIRES:
  259. VOLTOT = VOL1 ET VOL2 ET VOL3 ET VOL4 ;
  260. VOLSOU = PP11 ET PP12 ET PP13 ET PP14 ET
  261. PP15 ET PP16 ET CB03 ;
  262.  
  263. *--- SURFACE INFERIEURE:
  264. S_INF1 = MANU QUA8 A1 AB1 B1 BE1 E1 EF1 F1
  265. AF1 ;
  266. S_INF2 = MANU QUA8 C1 CD1 D1 DF1 F1 EF1 E1
  267. CE1 ;
  268. S_INF = S_INF1 ET S_INF2 ;
  269.  
  270. *--- SURFACE SUPERIEURE:
  271. S_SUP1 = MANU QUA8 A4 AB4 B4 BE4 E4 EF4 F4
  272. AF4 ;
  273. S_SUP2 = MANU QUA8 C4 CD4 D4 DF4 F4 EF4 E4
  274. CE4 ;
  275. S_SUP = S_SUP1 ET S_SUP2 ;
  276.  
  277. OEIL1 = 15 12 10 ;
  278. SI(NEG GRAPH 'N');
  279. TRAC 'QUAL' OEIL1 VOLTOT ;
  280. FINSI;
  281.  
  282. *--- DONNEES DU PROBLEME DE THERMIQUE:
  283. *--- MODELISATION:
  284. MOD1 = MODE VOLTOT THERMIQUE ISOTROPE ;
  285. MOD2 = MODE S_SUP THERMIQUE CONVECTION ;
  286.  
  287. *--DONNEE DES CARACTERISTIQUES DU MATERIAU:
  288. KC = 100. ;
  289. MAT1 = MATE MOD1 'K' KC ;
  290. HC = 15. ;
  291. MAT2 = MATE MOD2 'H' HC ;
  292.  
  293. *--- CREATION DES MATRICES DE CONDUCTIVITE:
  294. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  295. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  296.  
  297. *--- TEMPERATURES IMPOSEES: BLOQUE + DEPI
  298. BB1 = BLOQ S_INF 'T' ;
  299. TI = 10. ;
  300. EE1 = DEPI BB1 TI ;
  301.  
  302. *--- FLUX CONVECTIF A LA SURFACE:
  303. TEX = 200. ;
  304. CNV1 = CONVECTION MOD2 MAT2 'T' TEX ;
  305.  
  306. *--- FLUX IMPOSE:
  307. PHI = 20. ;
  308. FLU1 = FLUX MOD1 PHI S_SUP ;
  309.  
  310. *--- SOURCE VOLUMIQUE:
  311. QVOL = 2. ;
  312. S1 = SOUR MOD1 QVOL VOLSOU ;
  313.  
  314. *-- ASSEMBLAGE DES TERMES DE CONDUCTIVITE:
  315. CCC1 = CND1 ET CND2 ET BB1 ;
  316. CCC2 = CND1 ET BB1 ;
  317.  
  318. *ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS:
  319. FFF1 = EE1 ET CNV1 ;
  320. FFF2 = EE1 ET FLU1 ;
  321. FFF3 = EE1 ET S1 ;
  322.  
  323. *--- 1ERE RESOLUTION:
  324. * TEMPERATURE IMPOSEE + CONVECTION
  325. CHTER1 = RESO CCC1 FFF1 ;
  326.  
  327. *--- POST-TRAITEMENT
  328. Z1 = COOR 3 F1;
  329. Z2 = COOR 3 F4;
  330.  
  331. *Valeur Théorique de la température en Z2
  332. TREF1=((KC * TI) + (HC * Z2 * TEX)) /
  333. ((KC ) + (HC * Z2 )) ;
  334.  
  335. TETA1 = EXTR CHTER1 T F4 ;
  336. RESI1=100. * (ABS((TETA1-TREF1)/TREF1));
  337.  
  338. MESS 'Temperature theorique :' TREF1' C';
  339. MESS 'Temperature calculee :' TETA1' C';
  340. MESS ' Soit un ecart de : ' RESI1 '%';
  341. RESITOT = PROG RESI1 ;
  342.  
  343. OEIL = 100. 0. 0. ;
  344. SI(NEG GRAPH 'N');
  345. TITR 'ELEMENTS PYR5 CUB8 TET4 PRI6 :
  346. CONVECTION' ;
  347. TRAC OEIL VOLTOT CHTER1 ;
  348. FINSI;
  349.  
  350. *--- 2EME RESOLUTION:
  351. * TEMPERATURE IMPOSEE + FLUX
  352. CHTER2 = RESO CCC2 FFF2 ;
  353.  
  354. TREF2= TI + (PHI * Z2 / KC );
  355.  
  356. TETA2 = EXTR CHTER2 T F4 ;
  357. RESI2=100. * (ABS((TETA2-TREF2)/TREF2));
  358. MESS 'Temperature theorique :' TREF2' C';
  359. MESS 'Temperature calculee :' TETA2' C';
  360. MESS ' Soit un ecart de : ' RESI2 '%';
  361. RESITOT = RESITOT ET (PROG RESI2) ;
  362.  
  363. SI(NEG GRAPH 'N');
  364. TITR 'ELEMENTS PYR5 CUB8 TET4 PRI6 :
  365. FLUX';
  366. TRAC OEIL VOLTOT CHTER2 ;
  367. FINSI;
  368.  
  369. *--- 3EME RESOLUTION:
  370. * TEMPERATURE IMPOSEE + SOURCE
  371. CHTER3 = RESO CCC2 FFF3 ;
  372.  
  373. TREF3= TI - (QVOL * (Z2**2) / (2 * KC)) +
  374. (QVOL * (Z2**2) / ( KC));
  375.  
  376. TETA3 = EXTR CHTER3 T F4 ;
  377. RESI3=100. * (ABS((TETA3-TREF3)/TREF3));
  378. MESS 'Temperature theorique :' TREF3' C';
  379. MESS 'Temperature calculee :' TETA3' C';
  380. MESS ' Soit un ecart de : ' RESI3 '%';
  381. RESITOT = RESITOT ET (PROG RESI3) ;
  382.  
  383. SI(NEG GRAPH 'N');
  384. TITR 'ELEMENTS PYR5 CUB8 TET4 PRI6 :
  385. SOURCE';
  386. TRAC OEIL VOLTOT CHTER3 ;
  387. FINSI;
  388.  
  389. * CODE DE FONCTIONNEMENT
  390. LIST RESITOT;
  391. SI((MAXI RESITOT) <EG 1.5);
  392. ERRE 0;
  393. SINO;
  394. ERRE 5;
  395. FINSI;
  396. FIN;
  397.  
  398.  
  399.  
  400.  
  401.  

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