Télécharger la plaquette Cast3M


fr en

Télécharger tran11.dgibi

Retour à la liste

Numérotation des lignes : 

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

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