Télécharger la plaquette Cast3M


fr en

Télécharger gatt_dpg2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  gatt_dpg2.dgibi
  2. *
  3. *     Test gatt_dpg2.dgibi: Jeux de données            
  4. *     ---------------------------------------         
  5. * 
  6. * REMARQUE : Cas-test identique a gatt_dpg.dgibi
  7. * ---------- mais avec elements BBAR.
  8. * 
  9. 'OPTI' 'DIME' 2 'MODE' 'PLAN' 'GENE' ;
  10. 'OPTI' 'ELEM' 'QUA8'                 ;
  11. 'OPTI' 'TRAC' 'PSC'                  ;
  12. 'TEMPS' 'ZERO'                       ;
  13. L = 'MOT' LIST ; F = 'MOT' FIN       ;
  14. *                                                                               
  15. ***********************************************************          
  16. *                TEST DE VALIDATION                       *
  17. *                ------------------                       *
  18. *               MODELE GATT_MONERIE                       *
  19. *  AFA3GLAA INCOMPRESSIBLE AVEC COUPLAGE DYNAMIQUE        *
  20. *                                                         *
  21. *  MAILLAGE:                                              *
  22. *     EPROUVETTE CARRE                                    *
  23. *                                                         *
  24. *  CHARGEMENT:                                            *
  25. *     CONTRAINTE DE COMPRESSION IMPOSEE (ESSAI DE FLUAGE) *
  26. *     TEMPERATURE CONSTANTE                               *
  27. *     DENSITE DE FISSIONS CONSTANTE                       *
  28. *     PRIMAIRE - MECANISMES 1 ET 2 - IRRADIATION          *
  29. *     DENSIFICATION - GONFLEMENT                          *
  30. ***********************************************************          
  31. *
  32. * repertoire des fichiers "divers"
  33. DIVERS = VENV 'CASTEM_DIVERS';
  34. *
  35. GRAPH    = FAUX ;
  36. LISTCOUR = VRAI ;
  37. *
  38. NE   = 1  ;
  39. H    = 1. ;
  40. *
  41. OO   = 0.              0.            ;
  42. A1   = H               0.            ;
  43. A2   = H               H             ;
  44. A3   = 0.              H             ;
  45. *
  46. LB   = 'DROIT' NE  OO A1             ;
  47. LD   = 'DROIT' NE  A1 A2             ;
  48. LH   = 'DROIT' NE  A2 A3             ;
  49. LG   = 'DROIT' NE  A3 OO             ;
  50. *
  51. SU1  = 'DALL' LB  LD  LH  LG  'PLAN' ;
  52. *
  53. MODL1= MODE SU1 MECANIQUE ELASTIQUE
  54.   VISCOPLASTIQUE GATT_MONERIE BBAR 'DPGE' OO      ;
  55. *
  56. * ----------------------*
  57. * Materiau GATT_MONERIE *
  58. * ----------------------*
  59. TE1  = 1300. + 273.  ;
  60. PO   = 4.9E-2        ;
  61. DF1  = 1.E18         ;
  62. *
  63. TA = @GATTPAR ('CHAINE' DIVERS '/fichier_gatt') ;
  64. *
  65. TA.'BP'   = 0.      ;
  66. TA.'POR0' = PO      ;
  67. TA.'CR'   = 0.16E-2 ;
  68. *
  69. MATREE = MATE MODL1 'YOUN' (TA.'YOUN') ; 
  70. PP = 'MANU' 'CHML' MODL1 'T' TE1 'PORO' PO RIGIDITE ; 
  71. EE = 'VARI' 'NUAG' MODL1 MATREE PP ;
  72. *
  73. MATR11 = 'MATE' MODL1  'YOUN' EE  'NU' 0.  'RHO' 10950.  'ALPH' 0. 'TALP' 0. 'TREF' 1. ;
  74. *
  75. MATR12 = 'MATE' MODL1 
  76.  'R'    (TA.'R')    'DG0' (TA.'DG0') 'DG'  (TA.'DGCR')
  77.  'K1'   (TA.'K1')   'M1'  (TA.'M1')  'Q1'  (TA.'Q1')  'N1'   (TA.'N1') 
  78.  'K2'   (TA.'K2')   'M2'  (TA.'M2')  'Q2'  (TA.'Q2')  'N2'   (TA.'N2') 
  79.  'OMEG' (TA.'OMEG') 'H'   (TA.'H')   'Q'   (TA.'Q')   'BETA' (TA.'BETA')
  80.  'K'    (TA.'K')    'A'   (TA.'A')   'Q3'  (TA.'Q3')  'N3'   (TA.'N3')
  81.  'CR'   (TA.'CR')   'CR1' (TA.'CR1') 'CR2' (TA.'CR2') 'CR3'  (TA.'CR3');
  82. *
  83. MATR13 = 'MATE' MODL1 
  84.  'KP'   (TA.'KPAF')   'AP' (TA.'AP')  'BP'  (TA.'BP')  'QP'  (TA.'QP') ;
  85. * 
  86. MATR14 = 'MATE' MODL1 
  87.  'ADEN' TA.'ADEN' 'KGON' TA.'KGON'. 
  88.  'POR0' (TA.'POR0') 'BUMI' (TA.'BUMI')  'EFIS' (TA.'EFIS') ;
  89. *
  90. * TYPE = 0. combustible UO2            sinon combustible AFA3GLAA
  91. * COMP = 0. combustible incompressible sinon compressible
  92. * DYN  = 0. couplage dynamique         sinon statique
  93. *
  94. MATR15 = 'MATE' MODL1 'TYPE' 1.   'COMP' 1.   'DYN' 1. ;
  95. *
  96. MATR16 = 'MATE' MODL1
  97.          'DYN1' (TA.'DYN1')   'DYN2' (TA.'DYN2')   'DYN3' (TA.'DYN3') ;
  98. *
  99. MATR1 = MATR11 'ET' MATR12 'ET' MATR13 'ET' MATR14 'ET'
  100.         MATR15 'ET' MATR16 ;
  101. *
  102. * -----------------------*
  103. * Conditions aux limites *
  104. * -----------------------*
  105. CLYB = 'BLOQ' 'UY' LB                  ;
  106. CLXG = 'BLOQ' 'UX' LG                  ;
  107. CLT  = CLYB 'ET' CLXG                  ;
  108. *
  109. TMIL = 1.E6 ;
  110. TFIN = 2.E6 ;
  111. *
  112. RHO0 = 1. - TA.'POR0' ;
  113. AKBU = TA.'EFIS'*270./238./10950./RHO0*DF1 ;
  114. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT 
  115. TBU = (TA.'BUMI' / AKBU) - 1. ;
  116. TTBU = ('ENTI' (TBU/10000.)) * 10000. ;
  117. *-------------------*
  118. * Instants calcules *
  119. *-------------------*
  120. 'SI' LISTCOUR ;
  121. *--------------
  122. LIST1 = 'PROG' 0  'PAS' 0.1 1 'PAS' 1 10 'PAS' 10 100 'PAS' 100 1000 
  123.                   'PAS' 500  15000 ;
  124. 'SINON' ;
  125. *--------
  126. LIST1 = 'PROG' 0  'PAS' 0.1 1 'PAS' 1 10 'PAS' 10 100 'PAS' 100 1000 
  127.                   'PAS' 500  30000 ;
  128. 'FINSI' ;
  129. *--------
  130. * ------------------------------*
  131. * Chargement cte en temperature *
  132. * ------------------------------*
  133. CHTEMP   = 'MANU' 'CHPO' SU1 1 'T' 1.           ;
  134. EVT      = 'EVOL' 'MANU' ('PROG' 0.  (2.*TFIN)) 
  135.                          ('PROG' TE1 TE1)       ;
  136. CHARTEMP = 'CHAR' 'T' CHTEMP EVT                ;
  137. * -------------------------------------*
  138. * Chargement cte en densite de fission *
  139. * -------------------------------------*
  140. CHFISS   = 'MANU' 'CHPO' SU1 1 'DFIS' 1.        ;
  141. EVF      = 'EVOL' 'MANU' ('PROG' 0.  (2.*TFIN)) 
  142.                          ('PROG' DF1 DF1)       ;
  143. CHARFISS = 'CHAR' 'DFIS' CHFISS EVF             ;
  144. * -----------------------*
  145. * Chargement en pression *
  146. * -----------------------*
  147. valpres  = 60E6                                    ;
  148. compr    = 'PRES' 'MASS' MODL1 valpres LH          ;
  149. EVP      = 'EVOL' 'MANU' ('PROG' 0.  1. (2.*TFIN)) 
  150.                          ('PROG' 0.  1.        1.) ;
  151. CHARMECA = 'CHAR' 'MECA' compr EVP                 ;
  152. *
  153. CHARTOT  = CHARMECA 'ET' CHARTEMP 'ET' CHARFISS ;
  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       ;
  167. TAB1.'MODELE' = MODL1                  ;
  168. TAB1.'CHARGEMENT' = CHARTOT            ;
  169. TAB1.'TEMPERATURES' . 0 = CHTEMP       ;
  170. TAB1.'VARIABLES_INTERNES' . 0 = VAR0   ; 
  171. TAB1.'CARACTERISTIQUES' = MATR1        ;
  172. TAB1.'TEMPS_CALCULES' = LIST1          ;
  173. TAB1.'TEMPS_SAUVES' = LIST1            ;
  174. *
  175. TMASAU                   = TABLE       ;
  176. tab1 . 'MES_SAUVEGARDES' = TMASAU      ;
  177. TMASAU .'DEFTO' = VRAI                 ;
  178. TMASAU .'DEFIN' = VRAI                 ;
  179. *
  180. PASAPAS TAB1                           ;
  181. *
  182. *
  183. *
  184. *****************************************
  185. *        CONTROLE DES RESULTATS         *
  186. *****************************************
  187. *
  188. *
  189. AP      = TA.'AP'        ;
  190. BP      = TA.'BP'        ;
  191. N1      = TA.'N1'        ;
  192. N2      = TA.'N2'        ;
  193. DG      = TA.'DGCR'      ;
  194. *
  195. * Facteurs multiplicatifs dus au dopage par le Chrome
  196. WC1 = 'TANH' ( (TA.'CR' - TA.'CR2') / TA.'CR3' ) ;
  197. WC1 = 1. + ( 0.5 * TA.'CR1' * (1. + WC1) )       ;
  198. CV  = 180. / PI                                  ;
  199. WC2 = 1. - ( 'COS' (CV * DG / TA.'DG0') )        ;
  200. WC2 = 2. * (TA.'DG0'**TA.'M2') * WC2             ;
  201. *
  202. * Calcul de TO pour la fonction de couplage dynamique
  203. TO = 'TANH' ( (TA.'DYN2' - TE1) / TA.'DYN3' ) ;
  204. TO = TA.'DYN1' * (1. + TO)                    ;
  205. TO = TO + 1.                                  ;
  206. *
  207. BUMI = TA.'BUMI' ;
  208. AKEVD = AKBU*TA.'KGON' ;
  209. ADEN = TA.'ADEN' ; 
  210. KGON = TA.'KGON' ;
  211. * -------------------------------------------------------------------*
  212. * Calcul de la cte AAAA intervenant ds le calcul de la def. de dens. *
  213. * -------------------------------------------------------------------*
  214. BUMAX0=60.D0*BUMI ;
  215. CRIT=1.D-10 ;
  216. 'REPE' BLOC 100 ;         
  217.   BUMAX = BUMI* ('EXP' (1. - (ADEN/(KGON*BUMAX0)))) ;
  218.   BUMAX = (0.2*BUMAX) + (0.8*BUMAX0) ;
  219.   TEST='ABS' ((BUMAX-BUMAX0)/BUMAX0) ;
  220.   'SI' ('<' TEST CRIT) ; 
  221.     'QUIT' BLOC ; 
  222.   'FINS' ;
  223.   BUMAX0=BUMAX ;
  224. 'FIN' BLOC ;
  225. AAAA = (-1.D0*RHO0*(ADEN-(KGON*BUMAX))) /
  226.                 ((1.D0+ADEN)*(LOG(BUMAX/BUMI)))  ;
  227. *
  228. * -----------------------*
  229. * Controle des resultats *
  230. * -----------------------*
  231. SS = TAB1 . 'CONTRAINTES'                   ;
  232. VV = TAB1 . 'VARIABLES_INTERNES'            ;
  233. IN = TAB1 . 'DEFORMATIONS_INELASTIQUES'     ;
  234. NCONT = ('DIME' (TAB1 . 'CONTRAINTES')) - 1 ; 
  235. *
  236. ERMAX1 = 0.        ;
  237. ERMAX2 = 0.        ;
  238. *
  239. Lzeit   = 'PROG' ;
  240. Lteta   = 'PROG' ;
  241. Lteta_c = 'PROG' ;
  242. Lepsv   = 'PROG' ;
  243. Lepsv_c = 'PROG' ;
  244. *
  245. ind0  = 9  ;
  246. zeit0 = 1. ;
  247. *
  248. epsv0 = 'MAXI' ('EXCO' IN.(ind0 + 1) 'EIYY') ; 
  249. *
  250. 'REPE' BLOC (NCONT - ind0) ;
  251. *
  252.   ind = ind0 + &BLOC ;
  253. *
  254. * ----------------------------------------------------
  255. * SMYY = valpres ET LES AUTRES COMPOSANTES SONT NULLES
  256. * ----------------------------------------------------
  257.   SMYY = -1. * valpres ;
  258.   SSM  = SMYY / 3.     ;
  259.   SSEQ = 'ABS' SMYY    ;
  260. * Deviateur des contraintes
  261.   SSPRIM2 = SMYY - SSM ;
  262. * Contrainte YY calculee
  263.   sy_c = 'MAXI' ('EXCO' SS.ind 'SMYY') ;
  264. * --------------------------------------------------------------
  265. * VITESSE DE DEFORMATION VISCO-PLASTIQUE EVP D"ORIGINE THERMIQUE
  266. * --------------------------------------------------------------
  267. *-- fluage primaire
  268. *
  269.   EVP02 = -1.5 * TA.'KPAF' * 
  270.    ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) * (SSEQ**AP) ;
  271. *
  272. *-- fluage secondaire (2 mecanismes)
  273. *
  274. * on teste la porosite :
  275.   PP = 'MAXI' ('EXCO' VV . &BLOC 'PORO')   ;
  276.   ECART1 = 'ABS' ((PP - PO) / PO)          ;
  277.   ERMAX1 = 'MAXI' ('PROG' ECART1 ERMAX1)   ;
  278. *
  279.   A1 = (N1* (PP**(-1./N1) - 1.)**(-2.*N1/(N1+1.)))      ;
  280.   B1 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N1/(N1+1.))) ;
  281.   B1 = B1 + (A1/4.)                                     ;
  282.   A1 = 0.                                               ;
  283. *
  284.   A2 = (N2* (PP**(-1./N2) - 1.)**(-2.*N2/(N2+1.)))      ;
  285.   B2 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N2/(N2+1.))) ;
  286.   B2 = B2 + (A2/4.)                                     ;
  287.   A2 = 0.                                               ;
  288. *
  289.   EVP12 = 0.5 * (TA.'K1' * (DG**TA.'M1') * WC1 *
  290.   ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) * 
  291.   (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  292.   ((A1*1.5*SSM) + (3.*B1*SSPRIM2))) ;
  293. *
  294.   EVP22 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) * 
  295.   (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  296.   ((A2*1.5*SSM) + (3.*B2*SSPRIM2))) ;
  297. * -----------------------------
  298. * FONCTION DE COUPLAGE STATIQUE
  299. * -----------------------------
  300.   TETA0 = 0.5 * TA.'BETA' *
  301.       (1. + 
  302.       ('TANH' ((TE1 - (TA.'OMEG' * (SSEQ**(-1.*TA.'Q')))) / TA.'H'))) ;
  303. * ------------------------------
  304. * FONCTION DE COUPLAGE DYNAMIQUE
  305. * ------------------------------
  306.   zeit  = TAB1.'TEMPS'.ind        ;
  307.   duree = zeit - zeit0            ; 
  308.   TETA = ( 1. + (zeit/TO) ) ** -1 ;
  309.   TETA = TETA0 * (1. - TETA)      ;
  310. * Fonction de couplage dynamique calculee
  311.   teta_c = 'MAXI' ('EXCO' VV.ind 'TETA') ;
  312. * ----------------------------------
  313. * VITESSE DE DEFORMATION IRRADIATION
  314. * ----------------------------------
  315.   EVIR2 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) * 
  316.        ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM2 ;
  317. * ------------------------------------
  318. * DEFORMATION VISCO-PLASTIQUE SELON YY
  319. * ------------------------------------
  320.  KFI      = 1.+(TA.'K'*DF1)                      ;
  321.  EPSVa    = EVIR2 + ( KFI * (EVP02+EVP12) )      ;
  322.  EPSVb    = KFI * (EVP22 - EVP12)                ;
  323.  INTGTETA = 'LOG' ( 1. + (duree/TO) )            ;
  324.  INTGTETA = TO * INTGTETA                        ;
  325.  INTGTETA = TETA0 * (duree - INTGTETA)           ;
  326.  epsv     = (EPSVa * duree) + (EPSVb * INTGTETA) ;
  327.  epsv     = epsv0 + epsv                         ;
  328. * Deformation viscoplastique calculee
  329.  epsv_c = 'MAXI' ('EXCO' IN.ind 'EIYY') ;
  330. * On teste la deformation viscoplastique :
  331.      ECART2 = 'ABS' ((epsv_c - epsv) / epsv)  ;
  332.      ERMAX2 = 'MAXI' ('PROG' ECART2 ERMAX2)   ;
  333. *
  334.   Lzeit   = Lzeit   'ET' ('PROG' zeit)   ;
  335.   Lteta   = Lteta   'ET' ('PROG' TETA)   ;
  336.   Lteta_c = Lteta_c 'ET' ('PROG' teta_c) ;
  337.   Lepsv   = Lepsv   'ET' ('PROG' epsv)   ;
  338.   Lepsv_c = Lepsv_c 'ET' ('PROG' epsv_c) ;
  339. *
  340. 'FIN' BLOC ;
  341. *
  342.   'SI' (ERMAX1 '&lt;EG' 1E-5)               ;
  343.      'ERRE' 0                            ;
  344.   'SINO'                                 ;
  345.      'MESS' 'POROSITE NON CONSTANTE'     ;
  346.      'MESS' 'ERREUR MAXIMALE POROSITE :' 
  347.              ERMAX1 '> 1E-5 '            ; 
  348.      'ERRE' 5                            ;
  349.   'FINS'                                 ;
  350. *
  351. 'SI' ( ERMAX2 '&lt;EG' 0.05)                                ;
  352.    'ERRE' 0                                              ;
  353. 'SINO'                                                   ;
  354.    'MESS' 'ERREUR MAXIMALE DEFORMATION VISCOPLASTIQUE :' 
  355.            ERMAX2 '> 0.05 '                              ; 
  356.    'ERRE' 5                                              ;
  357. 'FINS'                                                   ;
  358. *
  359. *
  360. 'SI' GRAPH ;
  361.  'TITR' 'Deformations planes general. : Fonction de couplage dynamique';
  362.  EVTETA    = 'EVOL' 'MANU' 'Temps' Lzeit 'TETA' Lteta   ;
  363.  EVTETA_C  = 'EVOL' 'MANU' 'Temps' Lzeit 'TETA' Lteta_c ;
  364.  TAD     = TABLE                                        ;
  365.  TAD . 1 = 'TIRC'                                       ;
  366.  TAD . 2 = 'TIRR'                                       ; 
  367.  TAD . 'TITRE' = TABLE                                  ;
  368.  TAD . 'TITRE' . 1 = 'TETA THEORIQUE'                   ;
  369.  TAD . 'TITRE' . 2 = 'TETA CALCULEE'                    ;
  370.  'DESS' (EVTETA 'ET' EVTETA_C) TAD 'LEGE'               ;
  371. *
  372.   'TITR' 'Deformations planes general. : Deformation viscoplastique' ;
  373.   EVEPSV    = 'EVOL' 'MANU' 'Temps' Lzeit 'EPSV' Lepsv               ;
  374.   EVEPSV_C  = 'EVOL' 'MANU' 'Temps' Lzeit 'EPSV' Lepsv_c             ;
  375.   TAD     = TABLE                                                    ;
  376.   TAD . 1 = 'TIRC'                                                   ;
  377.   TAD . 2 = 'TIRR'                                                   ; 
  378.   TAD . 'TITRE' = TABLE                                              ;
  379.   TAD . 'TITRE' . 1 = 'EPSV THEORIQUE'                               ;
  380.   TAD . 'TITRE' . 2 = 'EPSV CALCULEE'                                ;
  381.   'DESS' (EVEPSV 'ET' EVEPSV_C) TAD 'LEGE'                           ;
  382. 'FINS' ;
  383. *
  384. 'FIN';
  385.  
  386.  
  387.  
  388.  
  389.  
  390.  
  391.  
  392.  
  393.  
  394.  
  395.  
  396.  
  397.  

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