Télécharger la plaquette Cast3M


fr en

Télécharger prot1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  prot1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. *                                                                  *
  6. *    Test de la procedure PROT : Projection de Temperature d'un    *
  7. *    massif sur une coque                                          *
  8. *                                                                  *
  9. ********************************************************************
  10. *
  11. OPTI ECHO 0;
  12. GRAPH = 'N';
  13. *
  15. *
  16. *==================================================================*
  17. *             geometrie : maillage                                 *
  18. *==================================================================*
  19. *
  20. * Dimension en metres
  21. *
  22. * Points
  23. *
  24. P1 = -3. -1. 4.; P1E = -3.01 -1.01 4.;
  25. P3 =  3. -1. 4.; P3E =  3.01 -1.01 4.;
  26. P4 = -3.  1. 4.; P4E = -3.01  1.01 4.; 
  27. P6 =  3.  1. 4.; P6E =  3.01  1.01 4.;
  28. P01 = 0.  0. 5.215 ;
  29. P02 = 0.  0. 4.995  ;
  30. C0 =  0.  0. 0.;
  31.  
  32. OPTI ELEM SEG3;
  33. L1 = P1 CERC 10 C0 P3 CERC 4 C0 P6 CERC 10 C0 P4 CERC 4 C0 ;
  34.  
  35. OPTI ELEM SEG2;
  36.  
  37. L1L = P1E CERC 12 C0 P3E CERC 5 C0 P6E CERC 12 C0 P4E CERC 5 C0 ;
  38. L2 = P3E CERC 5 C0 P6E  ;
  39. L3 = P1E CERC 5 C0 P4E  ;
  40.  
  41. L1B = L1L PROJ 'CONI' C0 'SPHE' C0 P01;             
  42. L1C = L1L PROJ 'CONI' C0 'SPHE' C0 P02;             
  43. L2B = L2 PROJ 'CONI' C0 'SPHE' C0 P01;             
  44. L3C = L3 PROJ 'CONI' C0 'SPHE' C0 P02;             
  45.  
  46. *  Creation des objets
  47.  
  48. OPTI ELEM TRI6 ;
  49. PLAN1 = SURF  L1 SPHERIQUE C0 ;
  50.  
  51. OPTI ELEM TRI3 ;
  52.  
  53. PLAN1C = SURF L1C SPHERIQUE C0 ;
  54. PLAN1B = PLAN1C PROJ 'CONI' C0 'SPHE' C0 P01;
  55.  
  56. OPTI ELEM PRI6 ;
  57. GEOV = PLAN1C VOLU PLAN1B  ;
  58. ELIM 1.E-3 (L1B ET L2B ET GEOV ET L3C);
  59.  
  60. SI (NEG GRAPH 'N') ;
  61.   trace (PLAN1 ET GEOV); 
  62. FINSI ;
  63.  
  64. * Definition du modele coque
  65. *----------------------------
  66.  
  67. MO1 = MODELI PLAN1 THERMIQUE COQ6 ;
  68. MO2 = MODELI PLAN1 THERMIQUE COQ6 ;
  69.  
  70. MO = MO1 ET MO2 ;
  71.  
  72. * Definition du materiau et de l'epaisseur
  73.  
  74. MA1 = MATE MO1 EPAI 0.1 EXCE .5E-1 ;
  75. MA2 = MATE MO2 EPAI 0.1 EXCE -.5E-1 ;
  76. MA  = MA1 ET MA2 ;
  77.  
  78. * Definition du modele massif
  79. *-----------------------------
  80.  
  81. MOV = MODELI GEOV THERMIQUE  ;
  82.  
  83. * Definition du materiau et de l'epaisseur
  84.  
  85. MAV = MATE MOV K 15. 'C' 460. RHO 8000. ;
  86.  
  87. * Definition des conditions aux limites
  88. *---------------------------------------
  89.  
  90. CL1 = BLOQ T L2B ;
  91. CL2 = BLOQ T L3C ;
  92.  
  93. RIT =  CL1 ET CL2;
  94.  
  95. * Temperature imposee
  96. *---------------------
  97.  
  98. TEM1 = DEPI CL1 500.;
  99. TEM2 = DEPI CL2 0.;
  100.  
  101. * Matrice conductivite
  102. *----------------------
  103.  
  104. CON1 = CONDUCTIVITE MOV MAV ;
  105. CONT = CON1 ET RIT ;
  106.  
  107. * Second membre
  108. *---------------
  109.  
  110. TEMT = TEM1 ET TEM2 ;
  111.  
  112. * Resolution du systeme
  113. *-----------------------
  114.  
  115. CHT1 = RESO CONT TEMT ;
  116.  
  117. CHT2 = CHAN CHAM CHT1 GEOV ;
  118.  
  119. * Projection sur la coque
  120. *-------------------------
  121.  
  122. CHT31 = CALP MO1 MA1 CHT2 ;
  123. CHT32 = CALP  MO2 MA2 CHT2 ;
  124. CHP31 = CHANGE 'CHPO'  MO1 CHT31;                        
  125. CHP32 = CHANGE 'CHPO'  MO2 CHT32;                        
  126.  
  127. CHT31T = EXCO CHT31 T ; 
  128. CHT31TI = EXCO CHT31 TINF ; 
  129. CHT31TS = EXCO CHT31 TSUP ; 
  130.  
  131. CHT32T = EXCO CHT32 T ; 
  132. CHT32TI = EXCO CHT32 TINF ; 
  133. CHT32TS = EXCO CHT32 TSUP ; 
  134.  
  135. SI (NEG GRAPH 'N') ;
  136.   trace CHT32TI MO2 ;
  137.   trace CHT32T MO2 ;
  138.   trace CHT32TS MO2 ;
  139.   trace CHT31TI MO1 ;
  140.   trace CHT31T MO1 ;
  141.   trace CHT31TS MO1 ;
  142. FINSI ;
  143.  
  144.  
  145. T31S3 = EXTR CHP31 TSUP P3;
  146. T32I1 = EXTR CHP32 TINF P1;
  147. *
  148. * Les temperatures TINF sur la coque associee au modele MO1 doit
  149. * etre nulle au point P3 ( Condition  au limite), et la temperature
  150. * TSUP sur la coque associee au modele MO2 doit etre egale a 500 en P1
  151. *
  152.  
  153. DIF = (ABS(T32I1)) + (ABS(T31S3 - 500.)) ;
  154. SI (DIF <EG 5.) ;
  155.   ERRE 0 ;
  156. SINO ;
  157.   ERRE 5 ;
  158. FINSI ;
  159.  
  160. FIN ;
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  

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