Télécharger ther7or.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ther7or.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. * Test Ther7or.dgibi: Jeux de données *
  7. * ----------------------------------- *
  8. * *
  9. SAUT PAGE ;
  10. *-----------------------------------------------------*
  11. * TEST THER7OR *
  12. * *
  13. * TEST DES CL DE TEMPERATURE POUR COQ8 ET COQ6 *
  14. * ORTHOTROPE *
  15. * -------------------------------------------- *
  16. * *
  17. * TEMPERATURE IMPOSEE + CONVECTION *
  18. * + FLUX *
  19. * + SOURCE *
  20. * *
  21. * Ce test permet de vérifier le bon fonctionnement *
  22. * des divers opérateurs thermiques de CASTEM2000 en *
  23. * ORTHOTROPIE *
  24. * *
  25. * UNE PLAQUE RECTANGULAIRE CONSTITUEE D'éléments COQ8 *
  26. * et coq6 est soumise à une température imposée à une *
  27. * de ses extrémités et sucessivement à une condition *
  28. * de convection forcée, à une condition de flux impo- *
  29. * -sé et à une condition de source volumique imposée. *
  30. * *
  31. * *
  32. * Il s'agit du meme probleme que pour le test ther71. *
  33. * Mais on utilise un materiau orthotrope au lieu d'un *
  34. * materiau isotrope avec la meme conductivite suivant *
  35. * l'axe A1A7. Or le gradient de temperature etant *
  36. * parallele a cet axe, on doit obtenir un resultat *
  37. * identique. *
  38. * On compare donc les resultats obtenus a ceux du *
  39. * test ther71 pour les elements coques. *
  40. * *
  41. * Test elabore par Mr Dowlatyari 03/91. *
  42. * Test modifie pour evolution 21/05/92. *
  43. * *
  44. *-----------------------------------------------------*
  45. opti echo 1;
  46. opti trac bens;
  47. graph = 'N';
  48. SAUT PAGE ;
  49. *
  50. TITRE 'C.L. DE TEMPERATURE - ELEMENTS COQ8 , COQ6 ' ;
  51. OPTION DIME 3 ELEM QUA8 ;
  52. TEMPS ;
  53. *
  54. *CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS
  55. *
  56. A1 = 0. 0. 0. ; B1 = 0. 3.0 0. ;AB1 =0. 1.5 0. ;
  57. A2= 0. 0. 1. ; B2 = 0. 3. 1. ; AB2 = 0. 1.5 2. ;
  58. A3 = 0. 0. 2. ; B3 = 0. 3. 2. ; AB4 = 0. 1.5 4. ;
  59. A4 = 0. 0. 3. ; B4 = 0. 3. 3. ; AB5 = 0. 1.5 6. ;
  60. A5 = 0. 0. 4. ;B5 = 0. 3. 4. ; AB3= 0. 1.5 3. ;
  61. A6 = 0. 0. 5. ;B6 = 0. 3. 5. ;
  62. A7 = 0. 0. 6. ;B7 = 0. 3. 6. ;
  63. M=0. 0. 0. ;N=2. 0. 0. ;
  64. *
  65. ANG=20.0 ;
  66. A1=A1 TOUR ANG M N ;A2=A2 TOUR ANG M N ;
  67. A3=A3 TOUR ANG M N ;A4=A4 TOUR ANG M N ;
  68. A5=A5 TOUR ANG M N ;A6=A6 TOUR ANG M N ;
  69. A7=A7 TOUR ANG M N ;B1=B1 TOUR ANG M N ;
  70. B2=B2 TOUR ANG M N ;B3=B3 TOUR ANG M N ;
  71. B4=B4 TOUR ANG M N ;B5=B5 TOUR ANG M N ;
  72. B6=B6 TOUR ANG M N ;B7=B7 TOUR ANG M N ;
  73. AB1=AB1 TOUR ANG M N ;AB2=AB2 TOUR ANG M N ;
  74. AB3=AB3 TOUR ANG M N ;AB4=AB4 TOUR ANG M N ;
  75. AB5=AB5 TOUR ANG M N ;
  76. *------- CREATION DES QUADRANGLES ET TRIANGLES -------
  77. *
  78. QH01 = MANU QUA8 A1 AB1 B1 B2 B3 AB2 A3 A2 ;
  79. QH02 = MANU TRI6 A3 AB2 B3 B4 B5 AB3 ;
  80. QH03 = MANU TRI6 A3 AB3 B5 AB4 A5 A4 ;
  81. QH04 = MANU QUA8 A5 AB4 B5 B6 B7 AB5 A7 A6 ;
  82. *
  83. *
  84. QTOT = QH01 ET QH02 ET QH03 ET QH04 ;
  85.  
  86. si ( NEG GRAPH 'N' );
  87. OEIL = 100. 100. 100. ;
  88. TRACE OEIL QTOT;
  89. finsi;
  90. QSOUR = QH04 ;
  91. *
  92. *
  93. L_INF = MANU SEG3 A1 AB1 B1 ;
  94. *
  95. S_COT = QH04;
  96. *
  97. * ATTENTION : il faut réorienter les éléments
  98. * (OPERATEUR ORIENT) si on n'est pas sur que tous les
  99. * éléments sont orientés d'une façon cohérante
  100. *
  101. *----------- DONNEES DU PROBLEME DE THERMIQUE --------
  102. *
  103. *-------------- MODELISATION ------------------------
  104. *
  105. MOD1 = MODE QTOT THERMIQUE ORTHOTROPE COQ8 COQ6 ;
  106. MOD2 = MODE S_COT thermique CONVECTION SUPERIEURE COQ8 COQ6 ;
  107. *
  108. *------ DONNEES DES CARACTERISTIQUES DU MATERIAU ------
  109. *
  110. MAT1 = MATE MOD1 'DIRECTION' A2 'PARALLELE' 'K1' 100. 'K2' 1000. 'K3' 100. EPAI 0.5 ;
  111. MAT2 = MATE MOD2 'H' 15. ;
  112. *
  113. *------- CREATION DES MATRICES DE CONDUCTIVITE --------
  114. *
  115. *OPTION DONN 5 ;
  116.  
  117. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  118. CND2 = CONDUCTIVITE MOD2 MAT2 ;
  119. *
  120. *-------- TEMPERATURES IMPOSEES: BLOQUE + DEPI --------
  121. *
  122. BB1 = BLOQUE L_INF 'TINF' ;
  123. EE1 = DEPI BB1 10. ;
  124. BB2 = BLOQUE L_INF 'TSUP' ;
  125. EE2 = DEPI BB2 10. ;
  126. BB3 = BLOQUE L_INF 'T ' ;
  127. EE3 = DEPI BB3 10. ;
  128. BB = BB1 ET BB2 ET BB3 ;
  129.  
  130. EE = EE1 ET EE2 ET EE3 ;
  131.  
  132. *
  133. *--- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ----
  134. *
  135. CNV1 = CONVECTION MOD2 MAT2 'TSUP' 200.;
  136. *
  137. *--------------- FLUX IMPOSE --------------------------<-div>=/di>
  138. *
  139. FLU1 = FLUX MOD1 200. S_COT SUPE ;
  140. *
  141. *--------------- SOURCE DE CHALEUR --------------------
  142. *
  143. list qsour;
  144. S1 = SOURCE MOD1 MAT1 1000. QSOUR ;
  145. list s1;
  146. *
  147. *-------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE -------
  148. *
  149. CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ;
  150. *
  151. *------ ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS -----
  152. *
  153. FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ;
  154. *FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;
  155. *
  156. *----------- 1ERE RESOLUTION : CONVECTION -------------
  157. *
  158. *OPTION DONN 5 ;
  159. CHTER1 = RESOUDRE CCC1 FFF1 ;
  160. *
  161. *---------------- 2EME RESOLUTION : FLUX --------------
  162. *
  163. CHTER2 = RESOUDRE CCC2 FFF2 ;
  164. *
  165. *---------------- 3EME RESOLUTION: SOURCE -------------
  166. *
  167. CHTER3 = RESOUDRE CCC2 FFF3 ;
  168. *
  169. *----- POST-TRAITEMENT: TRACE DES CHAMPS RESULTATS ----
  170. *------ ET EXTRACTION DES VALEURS CARACTERISTIQUES ----
  171. *
  172. *LISTE CHTER1 ;
  173. *LISTE CHTER2 ;
  174. *LISTE CHTER3 ;
  175. *
  176. OEIL= 1000. 0. 0. ;
  177. *
  178. SI (NEG GRAPH 'N');
  179. TITRE 'ELEMENTS COQ6,COQ8 : CONVECTION' ;
  180. CHTER1I = 'EXCO' CHTER1 'TINF' 'NOID' 'TINF' ;
  181. TRACER OEIL QTOT CHTER1I ;
  182. TITRE 'ELEMENTS COQ6,COQ8 : FLUX IMPOSE' ;
  183. CHTER2I = 'EXCO' CHTER2 'TINF' 'NOID' 'TINF' ;
  184. TRACER OEIL QTOT CHTER2I ;
  185. TITRE 'ELEMENTS COQ6,COQ8 : SOURCE' ;
  186. CHTER3I = 'EXCO' CHTER3 'TINF' 'NOID' 'TINF' ;
  187. TRACER OEIL QTOT CHTER3I ;
  188. FINSI;
  189.  
  190. T1AB1 = EXTR CHTER1 TSUP AB1 ;
  191. T1AB2=EXTR CHTER1 TSUP AB2 ;
  192. T1AB3 = EXTR CHTER1 TSUP AB3 ;
  193. T1AB4=EXTR CHTER1 TSUP AB4 ;
  194. T1AB5 = EXTR CHTER1 TSUP AB5 ;
  195. T2AB1= EXTR CHTER2 TSUP AB1 ;
  196. T2AB2=EXTR CHTER2 TSUP AB2 ;
  197. T2AB3 = EXTR CHTER2 TSUP AB3 ;
  198. T2AB4=EXTR CHTER2 TSUP AB4 ;
  199. T2AB5 = EXTR CHTER2 TSUP AB5 ;
  200. T3AB1 = EXTR CHTER3 TSUP AB1 ;
  201. T3AB2=EXTR CHTER3 TSUP AB2 ;
  202. T3AB3 = EXTR CHTER3 TSUP AB3 ;
  203. T3AB4=EXTR CHTER3 TSUP AB4 ;
  204. T3AB5 = EXTR CHTER3 TSUP AB5 ;
  205. *
  206.  
  207. SI ( NEG GRAPH 'N') ;
  208.  
  209. DIST = PROG 0. 2. 4. 6. ;
  210. DIST1 = PROG 0. 2. 3. 4. 6. ;
  211. TEMPR1=PROG 10. 71.44 132.03 158.86 ;
  212. TEMPC1= PROG T1AB1 T1AB2 T1AB3 T1AB4 T1AB5 ;
  213. EV1R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPR1 ;
  214. EV1C = EVOL =span style="color: #0033FF;">MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC1;
  215. TAB = 'TABLE' ;
  216. TAB.1 = 'MARQ PLUS' ;
  217. TAB.2 = 'MARQ ETOI' ;
  218. TITRE 'COQ6 ET COQ8 ;TEMPERATURE + CONVECTION ;+ : MASSIF,* : COQUE' ;
  219. DESSIN TAB (EV1R ET EV1C ) ;
  220. TEMPR2=PROG 10. 26.58 42.878 50.984 ;
  221. TEMPC2= PROG T2AB1 T2AB2 T2AB3 T2AB4 T2AB5 ;
  222. EV2R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPR2 ;
  223. EV2C = EVOL MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC2;
  224. TITRE 'COQ6 ET COQ8;TEMPERATURE + FLUX IMPOSE ;+ : MASSIF,* : COQUE';
  225. DESSIN TAB (EV2R ET EV2C ) ;
  226. TEMPR3=PROG 10. 51.43 91.75 111.75 ;
  227. TEMPC3= PROG T3AB1 T3AB2 T3AB3 T3AB4 T3AB5 ;
  228. EV3R = EVOL MANU 'DISTANCE' DIST 'TEMPERATURE' TEMPR3 ;
  229. EV3C = EVOL MANU 'DISTANCE' DIST1 'TEMPERATURE' TEMPC3;
  230. TITRE 'COQ6 ET COQ8 ;TEMPERATURE + SOURCE ;+ MASSIF,* : COQUE ';
  231. DESSIN TAB (EV3R ET EV3C ) ;
  232. TEMPS ;
  233. FINSI;
  234. T1Aq = 10.; T1Bq = 71.44; t1Cq = 132.03; T1DQ = 158.86;
  235. T2Aq = 10.; T2Bq = 26.58; t2Cq = 42.878; T2DQ = 50.984;
  236. T3AQ = 10.; T3BQ = 58.38; T3CQ = 105.34 ; T3DQ = 123.33;
  237.  
  238.  
  239.  
  240. *******************************************************
  241. ************** TEST DE FONCTIONNEMENT ****************
  242. *******************************************************
  243. mess t1ab1 T1Aq T1AB2 T1Bq;
  244. mess T1AB4 T1Cq T1AB5 T1Dq;
  245. mess T2AB1 T2Aq T2AB2 T2Bq;
  246. mess T2AB4 T2Cq T2AB5 T2Dq;
  247. mess T3AB1 T3Aq T3AB2 T3Bq;
  248. mess T3AB4 T3Cq T3AB5 T3Dq;
  249. ERCA= ABS (( T1AB1-T1Aq)/T1AB1);
  250. ERCB= ABS (( T1AB2-T1Bq)/T1AB2);
  251. ERCC= ABS (( T1AB4-T1Cq)/T1AB4);
  252. ERCD= ABS (( T1AB5-T1Dq)/T1AB5);
  253. ERFA= ABS (( T2AB1-T2Aq)/T2AB1);
  254. ERFB= ABS (( T2AB2-T2Bq)/T2AB2);
  255. ERFC= ABS (( T2AB4-T2Cq)/T2AB4);
  256. ERFD= ABS (( T2AB5-T2Dq)/T2AB5);
  257. ERSA= ABS (( T3AB1-T3Aq)/T3AB1);
  258. ERSB= ABS (( T3AB2-T3Bq)/T3AB2);
  259. ERSC= ABS (( T3AB4-T3Cq)/T3AB4);
  260. ERSD= ABS (( T3AB5-T3Dq)/T3AB5);
  261.  
  262. ERTOT = PROG ERCA ERCB ERCC ERCD ERFA ERFB ERFC ERFD ERSA ERSB ERSC ERSD;
  263. ermax = maxi ertot;
  264. temps;
  265. SI ( ERMAX &lt;EG 0.05 );
  266. ERRE 0;
  267. SINON;
  268. ERRE 5;
  269. FINSI;
  270. fin;
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  

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