Télécharger soret_12.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : soret_13.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. * CAS TEST soret_13.dgibi *
  7. ************************************************************************
  8. * Test effet Soret
  9. * - 3D TRID
  10. * - ELEMENTS FINIS TESTES : 'PRI6'
  11. * - Regime permanent
  12. * - C(Z0) = C0
  13. * - C(Z1) = C1
  14. * - GRAD(T) choisi lineaire (T=A*(Z**2) / 2.D0 ==> GRAD(T)= A*Z)
  15. * - CONDUCTION et SORET sont en concurrence
  16. ************************************************************************
  17.  
  18.  
  19. ************************************************************************
  20. * OPTIONS & REGLAGES
  21. ************************************************************************
  22. GRAPH = 'N' ;
  23. 'OPTI' 'DIME' 3 ;
  24. 'OPTI' 'ELEM' 'TRI3' ;
  25.  
  26. 'SI' ('NEG' GRAPH 'N');
  27. 'OPTI' 'TRAC' 'X';
  28. 'SINO';
  29. 'OPTI' 'TRAC' 'PSC';
  30. 'FINS';
  31.  
  32. NBE1 = 5 ;
  33. NBE2 = 5 ;
  34. NBE3 = 15 ;
  35.  
  36. K = -7.D0;
  37. A = 1.D0;
  38.  
  39. Z0 = -0.5 ;
  40. Z1 = 1. ;
  41.  
  42. C0 = 1. ;
  43. C1 = 2. ;
  44.  
  45. ************************************************************************
  46. * MAILLAGE
  47. ************************************************************************
  48. P1 = 0. 0. Z0;
  49. P2 = 0. 1. Z0;
  50. P3 = 1. 1. Z0;
  51. P4 = 1. 0. Z0;
  52.  
  53. L1 ='DROI' NBE1 P1 P2;
  54. L2 ='DROI' NBE1 P4 P3;
  55.  
  56. S1 ='REGL' NBE2 L1 L2;
  57. S2 = S1 'PLUS' (0. 0. (Z1 - Z0));
  58. L4 ='COTE' 4 S2;
  59. P5 ='POIN' L4 'FINAL';
  60.  
  61. 'OPTI' 'ELEM' 'PRI6';
  62. V1 ='VOLU' NBE3 S1 S2 ;
  63. ENV1='ENVE' V1;
  64. ARE1='ARET' ENV1;
  65. PAR1='CHAN' 'POI1' ARE1 ;
  66. M1 ='POIN' PAR1 'DROIT' P1 P5 1.D-5;
  67.  
  68. L6 = 'ELEM' ARE1 'APPU' 'STRI' M1;
  69. 'TRAC' 'QUAL' ENV1;
  70. 'TRAC' 'QUAL' ARE1;
  71.  
  72. ************************************************************************
  73. * MODELE & MATERIAU
  74. ************************************************************************
  75. MOD1 ='MODE' V1 'DIFFUSION' 'FICK' 'INCO' 'HI ' 'QHI ';
  76. MAT1 ='MATE' MOD1 'KD ' 1.D0 ;
  77. MAT2 ='MANU' 'CHML' MOD1 'SCAL' K 'RIGIDITE' ;
  78.  
  79. ************************************************************************
  80. * Champ de temperature pour que le GRAD(T) soit constant
  81. ************************************************************************
  82. COO3 ='COOR' 3 V1 ;
  83. EVO1 ='EVOL' 'BLEU' 'CHPO' COO3 'SCAL' L6 ;
  84. LRZ ='EXTR' EVO1 'ORDO' 1 ;
  85. TEMP1 ='CHAN' 'COMP' (A / 2.D0 * (COO3 ** 2)) 'T' ;
  86. 'TITR' 'Evolution de la TEMPERATURE avec (abscisse curviligne)';
  87. EVOT ='EVOL' 'BLEU' 'CHPO' TEMP1 'T' L6 ;
  88. 'DESS' EVOT ;
  89.  
  90. ************************************************************************
  91. * RIGIDITE
  92. ************************************************************************
  93. RIG1 ='COND' MOD1 MAT1;
  94. RIG2 ='SORE' MOD1 MAT1 MAT2 TEMP1;
  95.  
  96. ************************************************************************
  97. * CONDITIONS AUX LIMITES
  98. ************************************************************************
  99. CLIM1 ='BLOQ' 'HI' S1 ;
  100. CLIM2 ='BLOQ' 'HI' S2 ;
  101.  
  102.  
  103. CIMP1 ='DEPI' CLIM1 C0 ;
  104. CIMP2 ='DEPI' CLIM2 C1 ;
  105.  
  106. ************************************************************************
  107. * RESOLUTION
  108. ************************************************************************
  109. SOL1 ='RESO' (RIG1 ET RIG2 ET CLIM1 ET CLIM2) (CIMP1 ET CIMP2);
  110. 'TITR' 'Isovaleurs de concentration en HI';
  111. 'TRAC' SOL1 V1;
  112.  
  113. ************************************************************************
  114. * POST-TRAITEMENT
  115. ************************************************************************
  116. * Solution NUMERIQUE
  117. EVONUM ='EVOL' 'ROUG' 'CHPO' SOL1 'HI ' L6 ;
  118. LCNUM ='EXTR' EVONUM 'ORDO' 1 ;
  119. EVONUM ='EVOL' 'ROUG' 'MANU' 'X' LRZ 'HI' LCNUM;
  120.  
  121. * Solution ANALYTIQUE
  122. EVO2 ='EVOL' 'CHPO' TEMP1 'HI' L6 ;
  123. FLO1 = ((-2*K*A)**0.5)/2. ;
  124. XDENO= ( 'EXP'(K/2.*A*(LRZ**2)) ) *
  125. (('ERF'(FLO1*Z1)) - ('ERF'(FLO1*Z0)));
  126. LNUME= (('ERF'(FLO1*LRZ))*C1*('EXP'(K/2.*A*(Z1**2))))-
  127. (('ERF'(FLO1*LRZ))*C0*('EXP'(K/2.*A*(Z0**2))))+
  128. (('ERF'(FLO1*Z1 ))*C0*('EXP'(K/2.*A*(Z0**2))))-
  129. (('ERF'(FLO1*Z0 ))*C1*('EXP'(K/2.*A*(Z1**2))));
  130.  
  131. LCANA = LNUME / XDENO;
  132.  
  133. EVOANA ='EVOL' 'BLEU' 'MANU' 'X' LRZ 'HICA' LCANA;
  134. 'DESS' (EVONUM ET EVOANA);
  135.  
  136. * Comparaison
  137. LCDIFF = LCNUM - LCANA ;
  138. 'TITR' 'Difference entre la Solutions NUMERIQUE et ANALYTIQUE';
  139. EVODIF ='EVOL' 'VERT' 'MANU' 'X' LRZ 'HICA' LCDIFF;
  140. 'DESS' EVODIF;
  141.  
  142.  
  143. * Critere d'erreur
  144. 'OPTI' 'ECHO' 0;
  145. ERRABS ='MAXI' 'ABS' LCDIFF;
  146. ERRREL = ERRABS / (('MAXI' LCANA) - ('MINI' LCANA)) ;
  147. 'MESS' '*********************************************' ;
  148. 'MESS' '* Erreur MAXIMALE absolue :' ERRABS '*';
  149. 'MESS' '* Erreur MAXIMALE relative :' ERRREL '*';
  150. 'MESS' '*********************************************' ;
  151. 'SI' (ERRREL >EG 1.47D-2);
  152. 'ERRE' 5;
  153. 'FINS';
  154.  
  155. 'FIN';
  156.  
  157.  
  158.  

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