Télécharger la plaquette Cast3M


fr en

Télécharger prot.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  prot.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. OPTI ELEM SEG2;
  16. *
  17. *==================================================================*
  18. *             geometrie : maillage                                 *
  19. *==================================================================*
  20. *
  21. * Dimension en metres
  22. *
  23. * Points
  24. *
  25. P1 = -3. -1. 5.; P1B = -3. -1. 5.5; P1C = -3. -1. 4.5;
  26. P2 =  0. -1. 5.; P2B =  0. -1. 5.5; P2C =  0. -1. 4.5;
  27. P3 =  3. -1. 5.; P3B =  3. -1. 5.5; P3C =  3. -1. 4.5; 
  28. P4 = -3.  1. 5.; P4B = -3.  1. 5.5; P4C = -3.  1. 4.5; 
  29. P5 =  0.  1. 5.; P5B =  0.  1. 5.5; P5C =  0.  1. 4.5; 
  30. P6 =  3.  1. 5.; P6B =  3.  1. 5.5; P6C =  3.  1. 4.5; 
  31.  
  32. L1 = P1 D 5 P2 D 3 P5 D 5 P4 D 3  ;
  33. L2 = P6 D 5 P5 D 3 P2 D 5 P3 ;
  34. L3 = P3 D 3 P6  ;
  35. ELIM 1.E-3 (L1 ET L2);
  36. L1B = P1B D 5 P3B D 3 P6B D 5 P4B D 3  ;
  37. L2B = P3B D 3 P6B  ;
  38. ELIM 1.E-3 (L1B ET L2B);
  39. *L1C = P1C D 5 P3C D 3 P6C D 5 P4C D 3 ;
  40. L2C = P1C D 3 P4C  ;
  41. *ELIM 1.E-3 (L1C ET L2C);
  42.  
  43. *  Creation des objets
  44.  
  45. OPTI ELEM TRI6 ;
  46. PLAN1 = SURF  L1 PLANE ;
  47. OPTI ELEM QUA8 ;
  48. PLAN2 = SURF (L2 ET L3) PLANE ;
  49.  
  50. OPTI ELEM TRI3 ;
  51.  
  52. PLAN1B = SURF L1B PLANE ;
  53. PLAN1C = PLAN1B 'PLUS' (0. 0. -1.) ;
  54. ELIM PLAN1C L2C 0.00001;
  55. *PLAN1C = SURF L1C PLANE ;
  56.  
  57. GEO  = (PLAN1 ET PLAN2) COUL ROUG ;
  58. ELIM GEO 1.E-7;
  59.  
  60. OPTI ELEM PRI6 ;
  61.  
  62. GEOV = PLAN1B VOLU PLAN1C ;
  63.  
  64. * Definition du modele coque
  65. *----------------------------
  66.  
  67. MO1 = MODELI PLAN1 THERMIQUE COQ6 ;
  68. MO2 = MODELI PLAN1 THERMIQUE COQ6 ;
  69. MO3 = MODELI PLAN2 THERMIQUE COQ8 ;
  70. MO4 = MODELI PLAN2 THERMIQUE COQ8 ;
  71.  
  72. MO13 =  MO1 ET MO3 ;
  73. MO24 =  MO2 ET MO4 ;
  74. MO = MO13 ET MO24 ;
  75.  
  76. * Definition du materiau et de l'epaisseur
  77.  
  78. MA1 = MATE MO13 EPAI 0.5 EXCE 2.5E-1 ;
  79. MA2 = MATE MO24 EPAI 0.5 EXCE -2.5E-1 ;
  80. MA  = MA1 ET MA2 ;
  81.  
  82. * Definition du modele massif
  83. *-----------------------------
  84.  
  85. MOV = MODELI GEOV THERMIQUE  ;
  86.  
  87. * Definition du materiau et de l'epaisseur
  88.  
  89. MAV = MATE MOV K 15. 'C' 460. RHO 8000. ;
  90.  
  91. * Definition des conditions aux limites
  92. *---------------------------------------
  93.  
  94. CL1 = BLOQ T L2B ;
  95. CL2 = BLOQ T L2C ;
  96.  
  97. RIT =  CL1 ET CL2;
  98.  
  99. * Temperature imposee
  100. *---------------------
  101.  
  102. TEM1 = DEPI CL1 500.;
  103. TEM2 = DEPI CL2 0.;
  104.  
  105. * Matrice conductivite
  106. *----------------------
  107.  
  108. CON1 = CONDUCTIVITE MOV MAV ;
  109. CONT = CON1 ET RIT ;
  110.  
  111. * Second membre
  112. *---------------
  113.  
  114. TEMT = TEM1 ET TEM2 ;
  115.  
  116. * Resolution du systeme
  117. *-----------------------
  118.  
  119. CHT1 = RESO CONT TEMT ;
  120.  
  121. CHT2 = CHAN CHAM CHT1 GEOV ;
  122.  
  123. *trace CHT1 (PLAN1B ET PLAN1C) ;
  124.  
  125. * Projection sur la coque
  126. *-------------------------
  127.  
  128. *CHT3 = PROI MO MA CHT2 ;
  129.  
  130. CHT31 = CALP  MO13 MA1 CHT2 ;
  131. CHT32 = CALP  MO24 MA2 CHT2 ;
  132. CHP31 = CHANGE 'CHPO'  MO13 CHT31;                        
  133. CHP32 = CHANGE 'CHPO'  MO24 CHT32;                        
  134.  
  135. CHT31T = EXCO CHT31 T ; 
  136. CHT31TI = EXCO CHT31 TINF ; 
  137. CHT31TS = EXCO CHT31 TSUP ; 
  138. CHT32T = EXCO CHT32 T ; 
  139. CHT32TI = EXCO CHT32 TINF ; 
  140. CHT32TS = EXCO CHT32 TSUP ; 
  141. *
  142. SI (NEG GRAPH 'N') ;
  143.   trace CHT32TI MO13 ;
  144.   trace CHT32T MO13 ;
  145.   trace CHT32TS MO13 ;
  146.   trace CHT31TI MO13 ;
  147.   trace CHT31T MO13 ;
  148.   trace CHT31TS MO13 ;
  149. FINSI ;
  150. *
  151. T31S3 = EXTR CHP31 TSUP P3;
  152. T32I1 = EXTR CHP32 TINF P1;
  153. *
  154. * Les temperatures TINF sur la coque associee au modele MO13 doit
  155. * etre nulle au point P3 ( Condition  au limite), et la temperature
  156. * TSUP sur la coque associee au modele MO24 doit etre egale a 500 en P1
  157. *
  158. DIF = (ABS(T32I1)) + (ABS(T31S3 - 500.)) ;
  159. SI (DIF <EG 1.E-2) ;
  160.   ERRE 0 ;
  161. SINO ;
  162.   ERRE 5 ;
  163. FINSI ;
  164.  
  165. FIN ;
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  

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