Télécharger la plaquette Cast3M


fr en

Télécharger choctvf.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  choctvf.dgibi
  2. *
  3. *************************** choctvf.dgibi ***********************
  4. *
  5. *
  6. *----------------------------------------------------------------
  7. * Thermique transitoire : Choc thermique, test élémentaire pour
  8. * de l'opérateur LAPN en formulation Volumes Finis.
  9. *----------------------------------------------------------------
  10. * Référence : J.F.Saccadura, Initiation aux transferts thermiques
  11. * Paris, Technique et Documentation, 1982.
  12. *----------------------------------------------------------------
  13. * Auteur : F.DABBENE (TTMF)                                 12/97
  14. *----------------------------------------------------------------
  15. *
  16. *---------------------------------
  17. * Calcul de la solution analytique
  18. *---------------------------------
  19. *
  20. 'DEBP' TANAL GEO*MAILLAGE  TPS*FLOTTANT T0*FLOTTANT T1*FLOTTANT
  21.              LONG*FLOTTANT VK*FLOTTANT  EM*FLOTTANT ;
  22. *
  23. N      = 1 ;
  24. XX YY  = 'COOR' GEO ;
  25. TEMPFA = 'MANU' 'CHPO' GEO 1 'SCAL' 0.D0 ;
  26. 'REPE' BOUCINT 45 ;
  27.      TERMESIN = XX * N * 180. / LONG                           ;
  28.      TERMESIN = 'SIN' TERMESIN                                 ;
  29.      TERMEXP  = (N * PI / LONG)** 2. * (-1.)  * TPS * VK / EM  ;
  30.      TERMEXP  = 'EXP' TERMEXP                                  ;
  31.      TEMPFA   = TEMPFA + ( TERMESIN * TERMEXP / N )            ;
  32.      N = N + 2                                                 ;
  33. 'FIN' BOUCINT ;
  34. TEMPFA = ( TEMPFA * (T0 - T1) * 4. / PI ) + T1 ;
  35. 'FINP' TEMPFA ;
  36. *
  37. *-----------------
  38. * Post traitement
  39. *-----------------
  40. *
  41. 'DEBP' POPO AZ/'TABLE' ;
  42. TTRA = (RV . 'PASDETPS' . 'NUPASDT' - 1) * (RV . 'INCO' . 'DT') ;
  43. RV . 'PASDETPS' . 'NUPASDT' = RV . 'PASDETPS' . 'NUPASDT' + 1 ;
  44. *
  45. PANAC1 = TANAL (HYTOT . 'CENTRE') TTRA T0 T1 L VK C ;
  46. *PANAC1 = 'NOMC' 'SCAL' PANAC1 ;
  47. EPR1   = RV . 'INCO' . 'CN' ;
  48. ERRP1  = 'ABS' (EPR1 - PANAC1) ;
  49. MAXTP1 = 'MAXI' ERRP1 ;
  50. 'SI' ('EGA' GRAPH 'O') ;
  51.      AV3 = 'EVOL' 'ROUG' 'CHPO' PANAC1 'SCAL' DRMIC ;
  52.      LTITRE = 'TEXT' 'Temperature au temps' TTRA ;
  53.     'TITR' LTITRE ;
  54.      EV3 = 'EVOL' 'VERT' 'CHPO' EPR1 SCAL DRMIC ;
  55.     'DESS' (EV3 'ET' AV3) ;
  56.      LTITRE = 'TEXT' 'Erreur Absolue T temps' TTRA ;
  57.     'TITR' LTITRE ;
  58.      EV4 = 'EVOL' 'VERT' 'CHPO' ERRP1 'SCAL' DRMIC ;
  59.     'DESS' EV4 ;
  60. 'FINS' ;
  61. 'FINPROC' MAXTP1 ;
  62. *
  63. *------------------
  64. * Options générales
  65. *------------------
  66. *
  67. GRAPH = 'N'           ;
  68. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' 'ECHO' 0 'TRAC' 'PS' ;
  69. *
  70. *
  71. *=========
  72. * MAILLAGE
  73. *=========
  74. *
  75. *
  76. *- Création des points supports des DROITES
  77. *
  78. INUMX = 25                   ;
  79. INUMY = 1                    ;
  80. INUM1 = INUMX - 1            ;
  81. L     = 2.                   ;
  82. LS2   = L / 2.D0             ;
  83. DX1   = LS2 / INUMX / 2.     ;
  84. H     = 1.                   ;
  85. X0    = DX1 * -1.            ;
  86. X1    = X0 + LS2             ;
  87. Y0    = 0.                   ;
  88. Y1    = Y0 + H               ;
  89. Y01   = Y0 + Y1 * 0.5        ;
  90. XG    = X0 + DX1             ;
  91. XD    = X1 - DX1             ;
  92. *
  93. A1    = X0 Y0  ;
  94. A3    = X1 Y0  ;
  95. D1    = X0 Y1  ;
  96. D3    = X1 Y1  ;
  97. C1    = XG Y01 ;
  98. C3    = XD Y01 ;
  99. *
  100. *- Création des DROITES frontieres
  101. *
  102. DRBAS = A3 'DROI' INUMX A1 ;
  103. DRGAU = A1 'DROI' INUMY D1 ;
  104. DRHAU = D1 'DROI' INUMX D3 ;
  105. DRDRO = D3 'DROI' INUMY A3 ;
  106. PELIM = DX1 / (5. * INUMX) ;
  107. *
  108. *- Creation maillage GEOMETRIQUE
  109. *
  110. PTOT2 = 'DALL' DRBAS DRGAU DRHAU DRDRO ;
  111. *
  112. *- Creation maillage HYBRIDE y compris sous-objets (cond. limites)
  113. *
  114. HYTOT = 'DOMA' PTOT2 PELIM ;
  115. EXT1  = 'MANU' 'POI1' C1   ;
  116. DRMIC = C1 'DROI' INUM1 C3 ;
  117. 'ELIM' HYTOT.CENTRE DRMIC PELIM ;
  118. *
  119. *
  120. *================
  121. * INITIALISATIONS
  122. *================
  123. *
  124. *  T0   : Température au temps t=0 régnant uniformément sur la barre
  125. *  T1   : Température au temps t>0 régnant aux extrémités de la barre
  126. *  C    : Coefficient d'emmagasinement
  127. *  VK   : Valeur du coefficient de conductivité thermique
  128. *  TMAX : Temps final
  129. *  TSUP : Temps pour conditions aux limites
  130. *  DT   : Pas de temps
  131. *
  132. T0   = 0.D0    ;
  133. T1   = 100.D0  ;
  134. C    = 1.D0    ;
  135. VK   = 1.D0    ;
  136. TMIN = 0.D0    ;
  137. TMAX = 100.0D0 ;
  138. DT   = 0.005D0 ;
  139. *
  140. RV = 'EQEX' HYTOT 'ITMA' 20
  141.      'OPTI' 'VF' 'IMPL' 'CENTREE'
  142.      'ZONE' HYTOT 'OPER' 'POPO'
  143.      'ZONE' HYTOT 'OPER' 'DFDT' C 'CNM' 'DT' 'INCO' 'CN'
  144.      'ZONE' HYTOT 'OPER' 'LAPN' 'DIFF'       'INCO' 'CN'
  145.      'CLIM' 'CN' 'TIMP' EXT1 T1
  146. ;
  147.  
  148. VEXT1 = 'MANU' 'CHPO' EXT1 1 'SCAL' T1 ;
  149.  
  150. RV . 'INCO' = TABLE 'INCO' ;
  151. RV . 'INCO' . 'CN'   = 'KCHT' HYTOT 'SCAL' 'CENTRE' T0 VEXT1 ;
  152. RV . 'INCO' . 'CNM'  = 'KCHT' HYTOT 'SCAL' 'CENTRE' T0 VEXT1 ;
  153. RV . 'INCO' . 'DIFF' = 'KCHT' HYTOT 'SCAL' 'CENTRE' VK ;
  154. RV . 'INCO' . 'DT'   = DT ;
  155. *
  156. *=======================
  157. * Resolution transitoire
  158. *=======================
  159. *
  160. EXEC RV ;
  161. *
  162. *
  163. *- Gestion ERREUR
  164. *  Ecart inferieur à EPS0 degré à t=0.1 sec
  165. *
  166. VAL1 = POPO (RV.'1POPO') ;
  167. 'MESS' 'MAX ERREUR ABSOLUE :' VAL1 ;
  168. *
  169. EPS0 = 1. ;
  170. LOG1 = VAL1 > EPS0 ;
  171. 'SI' ( LOG1 ) ;
  172.    'ERRE' 5 ;
  173. 'SINO'     ;
  174.    'ERRE' 0 ;
  175. 'FINS'     ;
  176. 'FIN' ;
  177.  
  178.  
  179.  

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