Télécharger uo2_cas4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : uo2_cas4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Test uo2_cas4.dgibi: Jeux de données
  6. * ---------------------------------------
  7. *
  8. 'OPTI' 'DIME' 2 'MODE' 'PLAN' 'GENE' ;
  9. 'OPTI' 'ELEM' 'QUA4' ;
  10. **'OPTI' 'ECHO' 0 ;
  11. 'TEMPS' 'ZERO' ;
  12. L = 'MOT' LIST ; F = 'MOT' FIN ;
  13. *
  14. ***********************************************************
  15. * TEST DE VALIDATION *
  16. * ------------------ *
  17. * MODELE UO2 *
  18. * *
  19. * MAILLAGE: *
  20. * EPROUVETTE CARRE *
  21. * *
  22. * CHARGEMENT: *
  23. * DEPLACEMENT IMPOSE *
  24. * TEMPERATURE *
  25. * POROSITE CONSTANTE *
  26. * DENSITE DE FISSIONS CONSTANTE *
  27. * PRIMAIRE - MECANISMES 1 ET 2 - IRRADIATION *
  28. * DENSIFICATION - GONFLEMENT *
  29. ***********************************************************
  30. *
  31. H = 1. ;
  32. *
  33.  
  34. OO = 0. 0. ;
  35. A1 = H 0. ;
  36. A2 = H H ;
  37. A3 = 0. H ;
  38. *
  39. LB = 'DROIT' 2 OO A1 ;
  40. LD = 'DROIT' 5 A1 A2 ;
  41. LH = 'DROIT' 3 A2 A3 ;
  42. LG = 'DROIT' 10 A3 OO ;
  43.  
  44. *
  45. SU1 = surf ( LB et LD et LH et LG ) 'PLAN' ;
  46. **trac su1 ;
  47.  
  48. *
  49. * -------*
  50. * Modele *
  51. * -------*
  52. MODL1= MODE SU1 MECANIQUE ELASTIQUE VISCOPLASTIQUE UO2 'DPGE' OO ;
  53. *
  54. * ----------------------*
  55. * Materiau GATT_MONERIE *
  56. * ----------------------*
  57. TE1 = 1300. + 273. ;
  58. PO = 4.9E-2 ;
  59. DF1 = 1.E18 ;
  60. *
  61. TA = @GATTPAR '/u2/castem/divers/fichier_gatt' ;
  62. *
  63. TA.'POR0' = PO ;
  64.  
  65. MATR11 = 'MATE' MODL1 'YOUN' TA.YOUN 'NU' (TA.NU) 'RHO' (TA.RHO)
  66. 'ALPH' (TA.ALPH) ;
  67. *
  68. MATR12 = 'MATE' MODL1
  69. 'R' (TA.'R') 'DG0' (TA.'DG0') 'DG' (TA.'DGCR')
  70. 'K1' (TA.'K1') 'M1' (TA.'M1') 'Q1' (TA.'Q1') 'N1' (TA.'N1')
  71. 'K2' (TA.'K2') 'M2' (TA.'M2') 'Q2' (TA.'Q2') 'N2' (TA.'N2')
  72. 'OMEG' (TA.'OMEG') 'H' (TA.'H') 'Q' (TA.'Q') 'BETA' (TA.'BETA')
  73. 'K' (TA.'K') 'A' (TA.'A') 'Q3' (TA.'Q3') 'N3' (TA.'N3')
  74. 'CR' (TA.'CR') 'CR1' (TA.'CR1') 'CR2' (TA.'CR2') 'CR3' (TA.'CR3');
  75. *
  76. MATR13 = 'MATE' MODL1
  77. 'KP' (TA.'KPUO') 'AP' (TA.'AP') 'BP' (TA.'BP') 'QP' (TA.'QP') ;
  78. *
  79. bumi = 1.e10*(TA.'BUMI') ;
  80. MATR14 = 'MATE' MODL1
  81. 'ADEN' TA.'ADEN' 'KGON' TA.'KGON'.
  82. 'POR0' (TA.'POR0') 'BUMI' bumi 'EFIS' (TA.'EFIS') ;
  83.  
  84. *
  85. * -----------------------------*
  86. * Materiau partie OTTOSEN *
  87. * -----------------------------*
  88. MATR17 = 'TAIL' MODL1 ;
  89. *
  90. * Resistance au cisaillement 'GS ' par defaut: 'YOUN'*1.8E-4
  91. * Ouverture a rupture 'WRUP' = 0 --> 'WRUP' = 2*'GFTR'/'LTR '
  92. * Relation bilineaire 'BILI' ou 'EPSB' = 0 --> 'BILI' = 0 (ouverture)
  93. *
  94. XGFTR = 100.8 ; XGFTR1=XGFTR ; XGFTR2=XGFTR ;
  95. XLTR = 5.E6 ; XLTR1 =XLTR ; XLTR2 =XLTR ;
  96. * Deformation a rupture 'EPSR' par defaut: 3*'LTR '/'YOUN'
  97. * XYOUN = 'MAXI' ('VARI' 'NUAG' MODL1 MATREE PP) ;
  98. * NEPSR = 3. ; XEPSR = NEPSR * (XLTR / XYOUN) ;
  99. *
  100. MATR18 = 'MATE' MODL1 'GFTR' XGFTR 'LTR ' XLTR 'BTR ' 0.2
  101. 'WRUP' 0. 'BILI' 0.
  102. 'GFT1' XGFTR1 'GFT2' XGFTR2
  103. 'GS1 ' 0. 'GS2 ' 0.
  104. 'LTR1' XLTR1 'LTR2' XLTR2
  105. 'WRU1' 0. 'WRU2' 0.
  106. 'BIL1' 0. 'BIL2' 0. ;
  107. *
  108. MATR1 = MATR11 'ET' MATR12 'ET' MATR13 'ET' MATR14 'ET'
  109. MATR17 'ET' MATR18 ;
  110. *
  111. * -----------------------*
  112. * Conditions aux limites *
  113. * -----------------------*
  114. CLYB = 'BLOQ' 'UY' LB ;
  115. CLXG = 'BLOQ' 'UX' LG ;
  116. CLT = CLYB 'ET' CLXG ;
  117. *
  118. TMIL = 1.E6 ;
  119. TFIN = 2.E6 ;
  120. *
  121.  
  122.  
  123. * ------------------------------*
  124. * Chargement en temperature *
  125. * ------------------------------*
  126. CHTEMP = 'MANU' 'CHPO' SU1 1 'T' 1. ;
  127. EVT = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  128. ('PROG' 400. 500.) ;
  129. CHARTEMP = 'CHAR' 'T' CHTEMP EVT ;
  130.  
  131.  
  132. chpor0 = 'MANU' 'CHPO' SU1 1 'PORO' 1.;
  133. EVT = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  134. ('PROG' 0.049 0.049) ;
  135. CHARPOR = 'CHAR' 'PORO' chpor0 EVT ;
  136. * -------------------------------------*
  137. * Chargement cte en densite de fission *
  138. * -------------------------------------*
  139. CHFISS = 'MANU' 'CHPO' SU1 1 'DFIS' 1. ;
  140. EVF = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  141. ('PROG' DF1 DF1) ;
  142. CHARFISS = 'CHAR' 'DFIS' CHFISS EVF ;
  143. * ---------------------------------*
  144. * Chargement en deplacement impose *
  145. * ---------------------------------*
  146.  
  147. EVP = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  148. ('PROG' 0. 1.) ;
  149. CL3 = 'BLOQ' LH 'UY' ;
  150. DEP3 = 'DEPI' cl3 100. ;
  151. CHARMECA = 'CHAR' 'DIMP' dep3 evp ;
  152.  
  153. CHARTOT = CHARMECA 'ET' CHARTEMP 'ET' CHARFISS 'ET' CHARPOR ;
  154. * ---------------------------------------*
  155. * Variables internes initiales 'PORO'=PO *
  156. * ---------------------------------------*
  157. VAR00 = 'ZERO' MODL1 'VARINTER' ;
  158. VAR01 = 'MANU' 'CHML' MODL1 'PORO' PO
  159. 'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  160. VAR0 = VAR00 + VAR01 ;
  161. *
  162. *
  163. TAB1 = TABLE ;
  164. TAB1.'TEMPERATURES' = TABLE ;
  165. TAB1.'VARIABLES_INTERNES'= TABLE ;
  166. TAB1.'BLOCAGES_MECANIQUES' = CLT et cl3;
  167. TAB1.'MODELE' = MODL1 ;
  168. TAB1.'CHARGEMENT' = CHARTOT ;
  169. TAB1.'TEMPERATURES' . 0 = 'MANU' 'CHPO' SU1 1 'T' 400. ;
  170. TAB1.'VARIABLES_INTERNES' . 0 = VAR0 ;
  171. TAB1.'CARACTERISTIQUES' = MATR1 ;
  172. TAB1.'TEMPS_CALCULES' = 'PROG' 0 PAS 0.1 1.2 ;
  173.  
  174. PASAPAS TAB1 ;
  175.  
  176. **-- Verification de la valeur du taux de combustion (variable interne BU)
  177.  
  178. * Valeur theorique du taux de combustion
  179. * bu_th = 3.2e-11*0.270*df1*t/(0.238*10950.*(1.-0.049))
  180.  
  181. AA = 3.2e-11*0.270*DF1/(0.238*10950.*(1.-0.049)) ;
  182.  
  183. **-- Taux de combustion a l instant 1.1 s (pas de non convergence)
  184.  
  185. BU_CALC1 = 'MINI' ('EXCO' TAB1.'VARIABLES_INTERNES' . 11 'BU') ;
  186. BU_CALC2 = 'MAXI' ('EXCO' TAB1.'VARIABLES_INTERNES' . 11 'BU') ;
  187.  
  188. ECART1 = 'ABS' (BU_CALC1-BU_CALC2) / BU_CALC1 ;
  189.  
  190. BU_TH = AA*1.1 ;
  191. ECART2 = 'ABS' (BU_CALC1-BU_TH) / BU_TH ;
  192.  
  193. **-- Taux de combustion a l instant 1.2 s (apres non convergence)
  194.  
  195. BU_CALC = 'MAXI' ('EXCO' TAB1.'VARIABLES_INTERNES' . 12 'BU') ;
  196. BU_TH = AA*1.2 ;
  197. ECART3 = 'ABS' (BU_CALC-BU_TH) / BU_TH ;
  198.  
  199. 'SI' (ECART1 '<EG' 1.E-8) ;
  200. 'ERRE' 0 ;
  201. 'SINON' ;
  202. 'ERRE' 5 ;
  203. 'FINSI' ;
  204.  
  205. 'SI' (ECART2 '<EG' 1.E-8) ;
  206. 'ERRE' 0 ;
  207. 'SINON' ;
  208. 'ERRE' 5 ;
  209. 'FINSI' ;
  210.  
  211. 'SI' (ECART3 '<EG' 1.E-8) ;
  212. 'ERRE' 0 ;
  213. 'SINON' ;
  214. 'ERRE' 5 ;
  215. 'FINSI' ;
  216. 'FIN' ;
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  

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