Télécharger la plaquette Cast3M


fr en

Télécharger tran9.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  tran9.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet =faux;
  5. *
  6.  
  7. *************************************************
  8. *  pour calcul complet mettre complet à : vrai;
  9. *************************************************
  10. GRAPH = 'N';
  11. SAUT PAGE ;
  12. OPTION  ECHO 1;
  13. *--------------------------------------------------------------------*
  14. *                                                                    *
  15. *  TEST TRAN9                                                        *
  16. *                                                                    *
  17. *--- Probl}me :                                                      *
  18. *                                                                    *
  19. * Cylindre soumis a un choc froid sur la peau interne                *
  20. *                                                                    *
  21. *                                                                    *
  22. *--- Description de la g{om{trie :                                   *
  23. *                                                                    *
  24. *                                                                    *
  25. *                           LI2                                      *
  26. *       ----------------------------------------------------         *
  27. *      | 1|                  2                              |        *
  28. *  LI1 |  |                                                 | LI3    *
  29. *       ----------------------------------------------------         *
  30. *                          LI4                                       *
  31. *                                                                    *
  32. *     Partie 1 : revetement                                          *
  33. *     Partie 2 : meteriau de base                                    *
  34. *                                                                    *
  35. * ---- Echange de surface :                                          *
  36. *                                                                    *
  37. *     flux nul sur LI2 ,LI3 et LI4                                   *
  38. *     convection forcee sur LI1                                      *
  39. *                                                                    *
  40. * ----- Materiau  et temperature exterieure :                        *
  41. *                                                                    *
  42. *     conductivite K et capacite C*RHO  dependent de la temperature  *
  43. *     coeffiecient d'echange H et la temperature exterieure Text     *
  44. *     dependent du temps                                             *
  45. *                                                                    *
  46. *------Resultats et comparaison  :                                   *
  47. *     on compare les profils de temperature a travers l'epaisseur    *
  48. *     pour des differents instants a ceux obtenus par le code TRIO   *
  49. *                                                                    *
  50. **********************************************************************
  51. *-----TITRE 'CHOC FROID SUR LA PEAU INTERNE D'UN CYLINDRE ' ;
  52. *OPTION  ECHO 0 ;
  53. OPTI DIME 2 ELEM QUA4 MODE AXIS;
  54. *
  55. *--------------------- Cr{ation de la g{om{trie -------------------------
  56. *
  57. DENS 1.;
  58. R=1994.;
  59. P1=R 0.;P2=R 5.;
  60. LI1=P1 D 1 P2;
  61. OBJ1=LI1 TRAN 1 (1.5 0.);LI2=COTE 2 OBJ1;LI4=INVE (COTE 4 OBJ1);
  62. OBJ2=(INVE (COTE 3 OBJ1)) TRAN 6 (6. 0.);LI2=LI2 ET (COTE 2 OBJ2);
  63. LI4=LI4 ET (INVE (COTE 4 OBJ2));REVET=OBJ1 ET OBJ2;
  64. OBJ3=(INVE (COTE 3 OBJ2)) TRAN 12 (12. 0.);LI2=LI2 ET (COTE 2 OBJ3);
  65. LI4=LI4 ET (INVE (COTE 4 OBJ3));
  66. OBJ4=(INVE (COTE 3 OBJ3)) TRAN 1 (2. 0.);LI2=LI2 ET (COTE 2 OBJ4);
  67. LI4=LI4 ET (INVE (COTE 4 OBJ4));
  68. OBJ5=(INVE (COTE 3 OBJ4)) TRAN 2 (6. 0.);LI2=LI2 ET (COTE 2 OBJ5);
  69. LI4=LI4 ET (INVE (COTE 4 OBJ5));
  70. OBJ6=(INVE (COTE 3 OBJ5)) TRAN 2 (10. 0.);LI2=LI2 ET (COTE 2 OBJ6);
  71. LI4=LI4 ET (INVE (COTE 4 OBJ6));
  72. OBJ7=(INVE (COTE 3 OBJ6)) TRAN 17 (170. 0.);LI2=LI2 ET (COTE 2 OBJ7);
  73. LI4=LI4 ET (INVE (COTE 4 OBJ7));LI3=COTE 3 OBJ7;
  74. MBASE=OBJ3 ET OBJ4 ET OBJ5 ET OBJ6 ET OBJ7;
  75. REVET=REVET COUL BLEU;MBASE=MBASE COUL VERT;
  76. CUVE=REVET ET MBASE;
  77. *
  78. *SI(NEG GRAPH 'N');
  79. *   TRAC CUVE ;
  80. *FINSI;
  81. *
  82. *--------------- CREATION DE LA TABLE THERMIQUE  -----------------------
  83. *
  84. ETAB=TABLE;
  85. CHTER0=MANU CHPO CUVE 1 T 287.;
  86. ETAB.'TEMPERATURES' = TABLE;
  87. ETAB.'TEMPERATURES' . 0 = CHTER0; 
  88. *
  89. *---------VARIATIONS DES PROPRIETES DES MATERIAUX AVEC LA TEMPERATURE
  90. *
  91. LTEMP=PROG 7. 21.1 37.8 93.3 148.9 204.4 260. 287.8 315.6 600.;
  92. LCOND1=PROG 1.438 1.445 1.454 1.54 1.618 1.696 1.770 1.809 1.852 2.18;
  93. LCOND2=PROG 3.729 3.773 3.825 3.963 4.05 4.050 4.015 3.963 3.929 3.30;
  94. LCOND1=LCOND1 * 1.E-2;LCOND2=LCOND2 * 1.E-2;
  95. LCAPA1=PROG 3.713 3.738 3.768 3.856 3.940 4.032 4.136 4.162 4.204 4.85;
  96. LCAPA2=PROG 3.420 3.480 3.550 3.760 3.930 4.100 4.270 4.360 4.430 5.41;
  97. LCAPA1=LCAPA1 * 1.E-3;LCAPA2=LCAPA2 * 1.E-3;
  98. EVOCOND1=EVOL MANU 'T' LTEMP 'K' LCOND1;
  99. EVOCOND2=EVOL MANU 'T' LTEMP 'K' LCOND2;
  100. EVOCAPA1=EVOL MANU 'T' LTEMP 'C' LCAPA1;
  101. EVOCAPA2=EVOL MANU 'T' LTEMP 'C' LCAPA2;
  102. LIST EVOCAPA2;
  103. MOD1= MODE REVET THERMIQUE ISOTROPE ;
  104. MOD2= MODE MBASE THERMIQUE ISOTROPE ;
  105.  
  106. MAT2 = 'MATE' MOD1 'K' EVOCOND1 'RHO' 1.;
  107. MAT3 = 'MATE' MOD2 'K' EVOCOND2 'RHO' 1.;
  108. MAT4 = 'MATE' MOD1 'C' EVOCAPA1 'RHO' 1.;
  109. MAT5 = 'MATE' MOD2 'C' EVOCAPA2 'RHO' 1.;
  110. *
  111. *----- COEFFICIENTS D'ECHANGE ET TEMPERATURES EXTERIEURES
  112. *
  113. LINST=PROG -1. 0. 50. 100. 300. 525. 600. 800. 1300. 1800. 2800.
  114. 3800. 4800. 5800. 5801.;
  115. LCECH=PROG 146.645 146.645 75.805 36.25 20.58 4.31 4.12 3.64 2.45 1.6
  116. 1.2 1. 1. 1. 1.;LCECH=LCECH * 1.E-3;
  117. LTEXT=PROG 287. 287. 292. 287. 274. 251. 229.5 176.3 66.6 27.9 9.6 7.3 
  118. 7. 7. 7.; 
  119. MODCONV= MODE LI1 THERMIQUE CONVECTION           ;
  120. EVOCV  = EVOL MANU 'TEMP' LINST  'H' LCECH    ;
  121. EVOTEXT= EVOL MANU 'TEMPS' LINST 'T' LTEXT    ;
  122.  
  123. CH_TE= 'MANU' 'CHPO' LI1 1 'T' 1. ;
  124. CHA1 = 'CHAR' 'TECO' CH_TE EVOTEXT;
  125.  
  126. MAT1 = 'MATE' MODCONV  'H' EVOCV;
  127. MATT = MAT1 ET MAT2 ET MAT3 ET MAT4 ET MAT5;
  128. ETAB.'CARACTERISTIQUES' = MATT;
  129. *
  130. *---- OBJET MODELE DE L'ENSEMBLE
  131. *
  132. OBJMOD = MOD1 ET MOD2 ET MODCONV  ;
  133. ETAB.'MODELE'=OBJMOD;
  134. *
  135. * -------  DEFINITION DU PAS DE TEMPS
  136. *
  137. si complet;
  138.  LIS = PROG 25. 'PAS' 25. 2300.;
  139. sinon;
  140.  LIS = PROG 50. PAS 50. 700.;
  141. finsi;
  142.  
  143. ETAB.'TEMPS_SAUVES' = LIS;
  144. ETAB.'TEMPS_CALCULES' = LIS;
  145.  
  146. ETAB.'CHARGEMENT' = CHA1;
  147. ETAB.'CRITERE' = 1E-5;
  148. ETAB.'PROCEDURE_THERMIQUE' = DUPONT;
  149. *
  150. * ---- APPEL A LA PROCEDURE PASAPAS
  151. *
  152. PASAPAS ETAB;
  153. *
  154. *------- POST TRAITEMENT
  155. *  RESULTATS SOUS FORME DE LISTING ET DESSINS
  156. *
  157.  
  158. I=0;
  159. nfoi = dime etab.'TEMPERATURES'   -  1; 
  160. REPETER BBBB nfoi;
  161. I=I + 1;
  162. TI=ETAB.'TEMPERATURES'.I;
  163.  
  164. TITRE '25 PROFILS DE TEMPERATURE A TRAVERS L EPAISSEUR t1=100,dt=200 ' ;
  165. EVII=EVOL 'CHPO' TI 'T' LI2; EVI = EXTR EVII 'COUR' 1 ;
  166. TITRE '25 PROFILS DE TEMPERATURE A TRAVERS L EPAISSEUR t1=100,dt=200 ' ;
  167. SI (I EGA 1);EVT=EVI;SINON;EVT=EVT ET EVI;FINSI;
  168. aa = 300. ;
  169. si (ABS ( ( etab .'TEMPS'.I - aa ) / aa) < 0.001 ) ;
  170. *mess ' saisie de ev300 ';
  171.  ev300=evi;fins;
  172. aa = 700;
  173. si (ABS ( ( etab .'TEMPS'.I - aa ) / aa) < 0.001 ) ;
  174. *mess ' saisie de ev700 ';
  175.  ev700=evi;fins;
  176. aa = 1100.;
  177. si (ABS ( ( etab .'TEMPS'.I - aa ) / aa) < 0.001 ) ;
  178. *mess ' saisie de ev110 ';
  179.  ev1100=evi;fins;
  180. aa = 2300.;
  181. si (ABS ( ( etab .'TEMPS'.I - aa ) / aa) < 0.001 ) ;
  182. *mess ' saisie de ev230 ';
  183.  ev2300=evi;fins;
  184. FIN BBBB;
  185. TITRE '25 PROFILS DE TEMPERATURE A TRAVERS L EPAISSEUR t1=100,dt=200 ' ;
  186. SI(NEG GRAPH 'N');
  187. *DESS EVT XBOR -50. 246.6 YBOR 0. 506.;
  188. DESS EVT XBOR -50. 250. YBOR 0. 350.;
  189. FINSI ;
  190. *
  191. *     COMPARAISON AVEC TRIO
  192. *
  193. *     LISTE DE TEMPERATURES  POUR T=300
  194. CAST1=EXTR EV300 'ORDO' 1    ;
  195. *cast1 = ev300;
  196. TRI1A=PROG
  197.  .27454906738653E+03 .27547912984181E+03 .27608124977709E+03
  198.  .27666950892879E+03 .27724431491926E+03 .27780607379700E+03
  199.  .27835518984928E+03 .27889206541745E+03 .27912904995927E+03
  200.  .27936067393693E+03 .27958701640942E+03 .27980815632783E+03
  201.  .28002417250759E+03 .28023514360122E+03 .28044114807151E+03
  202.  .28064226416523E+03 .28083856988731E+03 .28103014297556E+03
  203.  .28121706087590E+03 .28139940071813E+03 .28175065142581E+03 ;
  204. TRI1B=PROG
  205.  .28224509463441E+03 .28270242350016E+03 .28338752739051E+03
  206.  .28398397971857E+03 .28494443707211E+03 .28564692423383E+03
  207.  .28614701633725E+03 .28649224032164E+03 .28672211596256E+03
  208.  .28686858019860E+03 .28695669260044E+03 .28700551809609E+03
  209.  .28702909087935E+03 .28703737835863E+03 .28703718314855E+03
  210.  .28703294154444E+03 .28702739610385E+03 .28702213597865E+03
  211.  .28701801035175E+03 .28701542733999E+03 .28701455326992E+03 ;
  212. TRI1 = TRI1A ET TRI1B ;
  213. *     LISTE DE TEMPERATURES  POUR T=700
  214. CAST2=EXTR EV700 'ORDO' 1    ;
  215. * cast2 = ev700;
  216. TRI2A=PROG
  217.  .21554690597296E+03 .21977264321412E+03 .22251537853267E+03
  218.  .22520081984070E+03 .22783086154177E+03 .23040736716155E+03
  219.  .23293216971470E+03 .23540707203812E+03 .23645556043979E+03
  220.  .23748515345522E+03 .23849610932029E+03 .23948868564869E+03
  221.  .24046313932818E+03 .24141972641952E+03 .24235870205816E+03
  222.  .24328032035872E+03 .24418483432249E+03 .24507249574764E+03
  223.  .24594355514258E+03 .24679826164223E+03 .24845959556347E+03 ;
  224. TRI2B=PROG
  225.  .25083506287881E+03 .25307648621451E+03 .25653137481945E+03
  226.  .25965980006019E+03 .26504734663421E+03 .26943519883087E+03
  227.  .27299220748966E+03 .27586447176620E+03 .27817605388500E+03
  228.  .28003069368933E+03 .28151406519849E+03 .28269623067575E+03
  229.  .28363404722261E+03 .28437336661707E+03 .28495094072340E+03
  230.  .28539600090402E+03 .28573151934705E+03 .28597518368239E+03
  231.  .28614012621910E+03 .28623544927535E+03 .28626658215282E+03 ;
  232. TRI2 = TRI2A ET TRI2B ;
  233. *     LISTE DE TEMPERATURES  POUR T=1100
  234. si complet;
  235. CAST3=EXTR EV1100 'ORDO' 1    ;
  236. finsi;
  237. *cast3 = ev1100;
  238. TRI3A=PROG
  239.  .13489895413085E+03 .14152078087484E+03 .14585629635509E+03
  240.  .15013023226939E+03 .15434395121700E+03 .15849880467037E+03
  241.  .16259611374034E+03 .16663717033181E+03 .16825908659624E+03
  242.  .16986423236885E+03 .17145269670093E+03 .17302456936035E+03
  243.  .17457994081398E+03 .17611890221003E+03 .17764154536035E+03
  244.  .17914796272278E+03 .18063824738354E+03 .18211249303955E+03
  245.  .18357079398090E+03 .18501324507326E+03 .18785096822197E+03 ;
  246. TRI3B=PROG
  247.  .19199105477766E+03 .19599372132888E+03 .20236658084107E+03
  248.  .20837981739814E+03 .21938402518124E+03 .22910613889455E+03
  249.  .23764345899695E+03 .24509403542563E+03 .25155488355154E+03
  250.  .25712036580151E+03 .26188088918392E+03 .26592566797397E+03
  251.  .26933023583926E+03 .27216406343684E+03 .27448971486688E+03
  252.  .27636223011978E+03 .27782870437921E+03 .27892802922215E+03
  253.  .27969075855912E+03 .28013906367270E+03 .28028674645156E+03  ;
  254. TRI3 = TRI3A ET TRI3B ;
  255. *     LISTE DE TEMPERATURES  POUR T=2300
  256. si complet;
  257. CAST4=EXTR EV2300 'ORDO' 1   ;
  258. finsi;
  259. *cast4 = ev2300;
  260. TRI4A=PROG
  261.  .53565241215241E+02 .58487407171063E+02 .61745534731690E+02
  262.  .64984895392304E+02 .68205395123323E+02 .71406942107464E+02
  263.  .74589446709785E+02 .77752821447898E+02 .78969186323051E+02
  264.  .80181197280593E+02 .81388823375051E+02 .82592033864805E+02
  265.  .83790798213789E+02 .84985086093064E+02 .86174867382398E+02
  266.  .87360112171791E+02 .88540790763000E+02 .89716873670953E+02
  267.  .90888331625200E+02 .92055135571369E+02 .94374990472378E+02 ;
  268. TRI4B=PROG
  269.  .97821455001972E+02 .10122716950969E+03 .10681045302737E+03
  270.  .11227431915653E+03 .12282882732222E+03 .13286320239439E+03
  271.  .14235190823359E+03 .15127177921024E+03 .15962631384529E+03
  272.  .16739359114038E+03 .17455059943268E+03 .18107708900021E+03
  273.  .18695550479151E+03 .19217090461829E+03 .19671085955413E+03
  274.  .20056533924575E+03 .20372657052909E+03 .20618958261782E+03
  275.  .20795028877166E+03 .20900642127171E+03 .20935793566180E+03  ;
  276. TRI4 = TRI4A ET TRI4B ;
  277. si ( non complet);
  278. cast3 = tri3; cast4 = tri4;
  279. finsi;
  280. *
  281. *  CODE DE FONCTIONNEMENT
  282. *
  283. NB = DIME CAST1 ; I=0 ; RESI1=0.0 ;
  284. RESI2=0.0 ;RESI3=0.0 ;RESI4=0.0 ;
  285. xermax= 0.;
  286. REPETER BOUC1 NB ;
  287. I = I+1  ;
  288. THC1 = EXTR CAST1 I ;
  289. THT1 = EXTR TRI1  I ;
  290. err= ABS((THC1-THT1)/THT1); si ( err > xermax ) ;xermax=err; fins;
  291. RESI1= err + RESI1 ;
  292. THC2 = EXTR CAST2 I ;
  293. THT2 = EXTR TRI2  I ;
  294. err= ABS((THC2-THT2)/THT2); si ( err > xermax ) ;xermax=err; fins;
  295. RESI2= err + RESI2 ;
  296. THC3 = EXTR CAST3 I ;
  297. THT3 = EXTR TRI3  I ;
  298. err=ABS((THC3-THT3)/THT3);si ( err > xermax ) ;xermax=err; fins;
  299. RESI3= err + RESI3 ;
  300. THC4 = EXTR CAST4 I ;
  301. THT4 = EXTR TRI4  I ;
  302. err=ABS((THC4-THT4)/THT4);si ( err > xermax ) ;xermax=err; fins;
  303. RESI4= err + RESI4 ;
  304. FIN BOUC1 ;
  305. resi1 = resi1 / NB;
  306. resi2 = resi2 / NB;
  307. resi3 = resi3 / NB;
  308. resi4 = resi4 / NB;
  309. mess '   ';
  310. mess '   ';
  311. mess '   ';
  312. mess ' erreur max en tout point pour les instants' ;
  313. mess '   ';
  314. mess ' 300  700 1100 ET 2300  secondes : '  xermax;
  315. mess '   ';
  316. mess '   ';
  317. mess '   ';
  318. mess ' erreur moyenne pour les 4 courbes  : ';
  319. mess '   ';
  320. mess ' ' resi1  resi2 resi3 resi4 ;
  321. mess '   ';
  322. mess '   ';
  323. si complet ;
  324. erma = 1.E-2;
  325. sinon;
  326. erma = 1.5E-2;
  327. finsi;
  328. SI (RESI1 >EG  erma ) ;
  329.     ERRE  5 ;
  330. FINSI ;
  331.  
  332. SI (RESI2 >EG  erma ) ;
  333.     ERRE  5 ;
  334. FINSI ;
  335. SI (RESI3 >EG  erma ) ;
  336.     ERRE  5 ;
  337. FINSI ;
  338. SI (RESI4 >EG  erma ) ;
  339.     ERRE  5 ;
  340. FINSI ;
  341. TEMPS ;
  342. FIN;
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  
  362.  

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