Télécharger la plaquette Cast3M


fr en

Télécharger gatt_axi.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  gatt_axi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *     Test gatt_axi.dgibi: Jeux de données            
  6. *     ---------------------------------------         
  7. * 
  8. 'OPTI' 'DIME' 2 'MODE' 'AXIS'        ;
  9. 'OPTI' 'ELEM' 'QUA8'                 ;
  10. 'OPTI' 'TRAC' 'PSC'                  ;
  11. 'TEMPS' 'ZERO'                       ;
  12. L = 'MOT' LIST ; F = 'MOT' FIN       ;
  13. *
  14. *******************************************************          
  15. *                TEST DE VALIDATION                   *
  16. *                ------------------                   *
  17. *  MODELE GATT_MONERIE                                *
  18. *  UO2 DOPE CHROME   Cr=0.15% TAILLE DE GRAIN=7 µm    *
  19. *               COMPRESSIBLE AVEC COUPLAGE STATIQUE   *
  20. *                                                     *
  21. *  MAILLAGE:                                          *
  22. *     EPROUVETTE RECTANGULAIRE                        * 
  23. *                                                     *
  24. *  CHARGEMENT:                                        *
  25. *     DEPLACEMENT VERTICAL IMPOSE                     *
  26. *       TRACTION - COMPRESSION                        *
  27. *     TEMPERATURE CONSTANTE                           *
  28. *     DENSITE DE FISSIONS CONSTANTE                   *
  29. *     PRIMAIRE - MECANISMES 1 ET 2 - IRRADIATION      * 
  30. *     DENSIFICATION - GONFLEMENT                      *
  31. *******************************************************          
  32. *
  33. * repertoire des fichiers "divers"
  34. DIVERS = VENV 'CASTEM_DIVERS';
  35. *
  36. GRAPH    = FAUX ;
  37. LISTCOUR = VRAI ;
  38. *
  39. NA   = 1                             ;
  40. NR   = 1                             ;
  41. H    = 6.2e-3                        ;
  42. D    = 5.e-3                         ;
  43. R    = D/2.                          ;
  44. HS2  = H/2.                          ;
  45. *
  46. OO   = 0.              0.            ;
  47. A1   = R               0.            ;
  48. A2   = R               HS2           ;
  49. A3   = 0.              HS2           ;
  50. *
  51. LB   = 'DROIT' NR  OO A1             ;
  52. LD   = 'DROIT' NA  A1 A2             ;
  53. LH   = 'DROIT' NR  A2 A3             ;
  54. LG   = 'DROIT' NA  A3 OO             ;
  55. *
  56. SU1  = 'DALL' LB  LD  LH  LG  'PLAN' ;
  57. *
  58. MODL1= MODE SU1 MECANIQUE ELASTIQUE VISCOPLASTIQUE GATT_MONERIE  ;
  59. *
  60. * ----------------------*
  61. * Materiau GATT_MONERIE *
  62. * ----------------------*
  63. TE1  = 1400. + 273.  ;
  64. PO   = 4.52E-2       ;
  65. DF1  = 1.E18         ;
  66. TA = @GATTPAR ('CHAINE' DIVERS '/fichier_gatt') ;
  67. *
  68. * fluage primaire rendu significatif pour UO2 standard
  69.  TA.'KPUO'   = 1E12 * TA.'KPUO' ;
  70. *
  71. TA.'DG' = TA.'DG0' ;
  72. TA.'CR' = 0.15E-2  ;
  73. * 
  74. MATREE = MATE MODL1 'YOUN' (TA.'YOUN') ; 
  75. PP = 'MANU' 'CHML' MODL1 'T' TE1 'PORO' PO RIGIDITE ; 
  76. EE = 'VARI' 'NUAG' MODL1 MATREE PP ;
  77. *
  78. MATR11 = 'MATE' MODL1  'YOUN' EE   'NU' 0.3   'RHO' 10950.  'ALPH' 0. 'TALP' 0. 'TREF' 1. ;
  79. *
  80. MATR12 = 'MATE' MODL1 
  81.  'R'    (TA.'R')    'DG0' (TA.'DG0') 'DG'  (TA.'DG')
  82.  'K1'   (TA.'K1')   'M1'  (TA.'M1')  'Q1'  (TA.'Q1')  'N1'   (TA.'N1') 
  83.  'K2'   (TA.'K2')   'M2'  (TA.'M2')  'Q2'  (TA.'Q2')  'N2'   (TA.'N2') 
  84.  'OMEG' (TA.'OMEG') 'H'   (TA.'H')   'Q'   (TA.'Q')   'BETA' (TA.'BETA')
  85.  'K'    (TA.'K')    'A'   (TA.'A')   'Q3'  (TA.'Q3')  'N3'   (TA.'N3')
  86.  'CR'   (TA.'CR')   'CR1' (TA.'CR1') 'CR2' (TA.'CR2') 'CR3'  (TA.'CR3');
  87. *
  88. MATR13 = 'MATE' MODL1 
  89.  'KP'  (TA.'KPUO')   'AP' (TA.'AP')  'BP' (TA.'BP')  'QP'  (TA.'QP') ;
  90. *
  91. MATR14 = 'MATE' MODL1 
  92.  'ADEN' TA.'ADEN' 'KGON' TA.'KGON'. 
  93.  'POR0' (TA.'POR0') 'BUMI' (TA.'BUMI')  'EFIS' (TA.'EFIS') ;
  94. *
  95. * TYPE = 0. combustible UO2            sinon combustible AFA3GLAA
  96. * COMP = 0. combustible compressible   sinon incompressible
  97. * DYN  = 0. couplage statique          sinon dynamique
  98. *
  99. MATR1 = MATR11 'ET' MATR12 'ET' MATR13 'ET' MATR14 ;
  100. *
  101. * -----------------------------------------------------*
  102. * Conditions aux limites + deplacement impose sur CLZH *
  103. * -----------------------------------------------------*
  104. CLZB = 'BLOQ' 'UZ' LB                  ;
  105. CLZH = 'BLOQ' 'UZ' LH                  ;
  106. CLO  = 'BLOQ' 'UR' OO                  ;
  107. CLT  = CLZB 'ET' CLZH 'ET' CLO         ;
  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.01 1     'PAS' 10 100 
  123.                    'PAS' 100 1000   'PAS' 500  20000
  124.                    'PAS' 2000 TBU   'PAS' 1000  (TTBU+10000.) ;
  125. 'SINON' ;
  126. *--------
  127. LIST1 = 'PROG' 0   'PAS' 0.01 1     'PAS' 10 100 
  128.                    'PAS' 100 1000   'PAS' 500  20000
  129.                    'PAS' 2000 TBU   'PAS' 1000  (TTBU+10000.) ;
  130. LIST1 = LIST1 ET (PROG 
  131.                    'PAS' 1000. TMIL
  132.                    'PAS' 1000. (TMIL + ((TFIN - TMIL) / 200.))
  133.                    'PAS' 1000. (TMIL + ((TFIN - TMIL) / 20.))
  134.                    'PAS' 2000. (TMIL + ((TFIN - TMIL) / 4.))
  135.                    'PAS' 2000. TFIN) ;
  136. 'FINSI' ;
  137. *--------
  138. * --------------------------*
  139. * Chargement en temperature *
  140. * --------------------------*
  141. CHTEMP   = 'MANU' 'CHPO' SU1 1 'T' 1.          ;
  142. EVT      = 'EVOL' 'MANU' ('PROG' 0.  (2.*TFIN)) 
  143.                     ('PROG' TE1 TE1)           ;
  144. CHARTEMP = 'CHAR' 'T' CHTEMP EVT               ;
  145. * ---------------------------------*
  146. * Chargement en densite de fission *
  147. * ---------------------------------*
  148. CHFISS   = 'MANU' 'CHPO' SU1 1 'DFIS' 1.       ;
  149. EVF      = 'EVOL' 'MANU' ('PROG' 0.  (2.*TFIN)) 
  150.                     ('PROG' DF1 DF1)           ;
  151. CHARFISS = 'CHAR' 'DFIS' CHFISS EVF            ;
  152. * ---------------------------------*
  153. * Chargement en deplacement impose *
  154. * ---------------------------------*
  155. DEPZ   = 'DEPI' CLZH 1                         ; 
  156. DELTAZ = HS2/200.                              ;
  157. EVEPS  = 'EVOL' 'MANU' 
  158.           ('PROG' 0. TMIL   TFIN (2.*TFIN))
  159.           ('PROG' 0. DELTAZ 0.    0.)          ;
  160. CHARMECA = 'CHAR' 'DIMP' DEPZ EVEPS            ;
  161. *
  162. CHARTOT  = CHARMECA 'ET' CHARTEMP 'ET' CHARFISS ;
  163. * ---------------------------------------*
  164. * Variables internes initiales 'PORO'=PO *
  165. * ---------------------------------------*
  166. VAR00 = 'ZERO' MODL1 'VARINTER'                                    ;
  167. VAR01 = 'MANU' 'CHML' MODL1 'PORO' PO 
  168.                             'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  169. VAR0  = VAR00 + VAR01                                              ;
  170. *
  171. TAB1 = TABLE                           ;
  172. TAB1.'TEMPERATURES' = TABLE            ;
  173. TAB1.'VARIABLES_INTERNES'= TABLE       ;
  174. TAB1.'BLOCAGES_MECANIQUES' = CLT       ;
  175. TAB1.'MODELE' = MODL1                  ;
  176. TAB1.'CHARGEMENT' = CHARTOT            ;
  177. TAB1.'TEMPERATURES' . 0 = CHTEMP       ;
  178. TAB1.'VARIABLES_INTERNES' . 0 = VAR0   ; 
  179. TAB1.'CARACTERISTIQUES' = MATR1        ;
  180. TAB1.'TEMPS_CALCULES' = LIST1          ;
  181. TAB1.'TEMPS_SAUVES' = LIST1            ;
  182. *
  183. TMASAU                   = TABLE       ;
  184. tab1 . 'MES_SAUVEGARDES' = TMASAU      ;
  185. TMASAU .'DEFTO' = VRAI                 ;
  186. TMASAU .'DEFIN' = VRAI                 ;
  187. *
  188. PASAPAS TAB1 ;
  189. *
  190. AP    = TA.'AP'        ;
  191. BP    = TA.'BP'        ;
  192. N1    = TA.'N1'        ;
  193. N2    = TA.'N2'        ;
  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 * TA.'DG' / TA.'DG0') )   ;
  200. WC2 = 2. * (TA.'DG0'**TA.'M2') * WC2             ;
  201.  
  202. BUMI  = TA.'BUMI'      ;
  203. AKEVD = AKBU*TA.'KGON' ;
  204. ADEN  = TA.'ADEN'      ; 
  205. KGON  = TA.'KGON'      ;
  206. * -------------------------------------------------------------------*
  207. * Calcul de la cte AAAA intervenant ds le calcul de la def. de dens. *
  208. * -------------------------------------------------------------------*
  209. BUMAX0=60.D0*BUMI ;
  210. CRIT=1.D-10 ;
  211. 'REPE' BLOC 100 ;         
  212.   BUMAX = BUMI* ('EXP' (1. - (ADEN/(KGON*BUMAX0)))) ;
  213.   BUMAX = (0.2*BUMAX) + (0.8*BUMAX0) ;
  214.   TEST='ABS' ((BUMAX-BUMAX0)/BUMAX0) ;
  215.   'SI' ('<' TEST CRIT) ; 
  216.     'QUIT' BLOC ; 
  217.   'FINS' ;
  218.   BUMAX0=BUMAX ;
  219. 'FIN' BLOC ;
  220. AAAA = (-1.D0*RHO0*(ADEN-(KGON*BUMAX))) /
  221.                 ((1.D0+ADEN)*(LOG(BUMAX/BUMI)))  ;
  222. *
  223. * -----------------------*
  224. * Controle des resultats *
  225. * -----------------------*
  226. SS = TAB1 . 'CONTRAINTES'                   ;
  227. VV = TAB1 . 'VARIABLES_INTERNES'            ;
  228. IN = TAB1 . 'DEFORMATIONS_INELASTIQUES'     ;
  229. NCONT = ('DIME' (TAB1 . 'CONTRAINTES')) - 1 ; 
  230. *
  231. T1     = 0.        ;
  232. SMZZ1  = 0.        ;
  233. FF1    = TA.'POR0' ;
  234. ERMAX1 = 0.        ;
  235. ERMAX2 = 0.        ;
  236. *
  237. LTT     = 'PROG' ;
  238. LDFF    = 'PROG' ;
  239. LEVFF   = 'PROG' ;
  240. LEVP3   = 'PROG' ;
  241. LEVTH   = 'PROG' ;
  242. LEVCALC = 'PROG' ;
  243. LSMZZ   = 'PROG' ;
  244. * 
  245. 'REPE' BLOC NCONT ;
  246. *
  247. * -------------------------------------
  248. * VITESSE DE DEFORMATION THEORIQUE EVTH
  249. * -------------------------------------
  250.   T2   = 'EXTR' (TAB1.'TEMPS_SAUVES') (&BLOC + 1);
  251.   EVTH = DELTAZ / HS2 / TMIL                     ;
  252.   'SI' (> T2 (TMIL+0.1))                         ;
  253.     EVTH = -1. * DELTAZ / HS2 / (TFIN - TMIL)    ; 
  254.   'FINS'                                         ;
  255. * ----------------------------------------------------
  256. * SMZZ EST PREPONDERANTE DEVANT LES AUTRES CONTRAINTES
  257. * ----------------------------------------------------
  258.   SMZZ = 'MAXI' ('EXCO' (SS . &BLOC) 'SMZZ') ;
  259.   SSM  = SMZZ / 3.                           ;
  260.   SSEQ = 'ABS' SMZZ                          ;
  261. * Deviateur des contraintes
  262.   SSPRIM1 = -1. * SMZZ / 3                   ;
  263.   SSPRIM2 = -1. * SMZZ / 3                   ;
  264.   SSPRIM3 =  2. * SMZZ / 3                   ;
  265. * ---------------------------------------------
  266. * DEFORMATION VISCO-PLASTIQUE EQUIVALENTE EVPEQ
  267. * ---------------------------------------------
  268.   EVRR  = 'MAXI' ('EXCO' IN.&BLOC 'EIRR')                           ;
  269.   EVTT  = 'MAXI' ('EXCO' IN.&BLOC 'EITT')                           ;
  270.   EVZZ  = 'MAXI' ('EXCO' IN.&BLOC 'EIZZ')                           ;
  271.   EVH   = (EVRR+EVTT+EVZZ) / 3.                                     ;
  272.   EVPEQ = ((EVRR - EVH)**2) + ((EVTT - EVH)**2) + ((EVZZ - EVH)**2) ;
  273.   EVPEQ = ((2./3.) * EVPEQ)**0.5                                    ;
  274. * ------------------------------------------
  275. * VITESSE DE DEFORMATION VISCO-PLASTIQUE EVP
  276. * ------------------------------------------
  277. *-- fluage primaire
  278. *
  279. 'SI' ('<' EVPEQ 1E-10) ; EVPEQ = 1E-10 ; 'FINSI' ;
  280. *
  281.   EVP01 = 1.5 * TA.'KPUO' * 
  282.   ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) * 
  283.   ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM1 ;
  284. *
  285.   EVP02 = 1.5 * TA.'KPUO' * 
  286.   ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) * 
  287.   ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM2 ;
  288. *
  289.   EVP03 = 1.5 * TA.'KPUO' * 
  290.   ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) * 
  291.   ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM3 ;
  292. *
  293. *-- fluage secondaire (2 mecanismes)
  294. *
  295.   PP = 'MAXI' ('EXCO' VV . &BLOC 'PORO') ;
  296. *
  297.   A1 = (N1* (PP**(-1./N1) - 1.)**(-2.*N1/(N1+1.)))      ;
  298.   B1 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N1/(N1+1.))) ;
  299.   AB1 = (9. * A1) / ( (4.*B1) + A1 )                    ;
  300.   BA1 = B1 / ( B1 + (A1/4.) )                           ;
  301. *
  302.   GSSEQ = ((BA1 * (SSEQ**2)) + (AB1 * (SSM**2))) ** 0.5 ;
  303. *
  304.   A2 = (N2* (PP**(-1./N2) - 1.)**(-2.*N2/(N2+1.))) ;
  305.   B2 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N2/(N2+1.))) ;
  306. *
  307.   EVP11 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  308.   ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) * 
  309.   (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  310.   ((A1*1.5*SSM) + (3.*B1*SSPRIM1))) ; 
  311. *
  312.   EVP21 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  313.   (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  314.   ((A2*1.5*SSM) + (3.*B2*SSPRIM1))) ; 
  315. *
  316.   EVP12 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  317.   ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) * 
  318.   (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  319.   ((A1*1.5*SSM) + (3.*B1*SSPRIM2))) ; 
  320. *
  321.   EVP22 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) * 
  322.   (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  323.   ((A2*1.5*SSM) + (3.*B2*SSPRIM2))) ; 
  324. *
  325.   EVP13 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  326.   ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) * 
  327.   (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  328.   ((A1*1.5*SSM) + (3.*B1*SSPRIM3))) ; 
  329. *
  330.   EVP23 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) * 
  331.   (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  332.   ((A2*1.5*SSM) + (3.*B2*SSPRIM3))) ; 
  333. *
  334.   PSI1 = TA.'K1' * WC1 * (TA.'DG'**TA.'M1')/(N1+1.) *
  335.   ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) * 
  336.   (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1+1.)/2.)) ;
  337. *
  338.   PSI2 = TA.'K2' * WC2/(N2+1.) * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  339.   (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2+1.)/2.)) ;
  340. *
  341.   EVPC1 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  342.   (GSSEQ**(-1.*TA.'Q' - 2.)) *
  343.       (1. -
  344.   (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q'))))  / TA.'H'))**2))
  345.    * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) +  (1.5*BA1*SSPRIM1) ) ;
  346. *
  347.   EVPC2 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  348.   (GSSEQ**(-1.*TA.'Q' - 2.)) *
  349.       (1. -
  350.   (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q'))))  / TA.'H'))**2))
  351.    * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) +  (1.5*BA1*SSPRIM2) ) ;
  352. *
  353.   EVPC3 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  354.   (GSSEQ**(-1.*TA.'Q' - 2.)) *
  355.       (1. -
  356.   (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q'))))  / TA.'H'))**2))
  357.    * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) +  (1.5*BA1*SSPRIM3) ) ;
  358. *
  359.   TETA = 0.5 * TA.'BETA' *
  360.       (1. + 
  361.       ('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q')))) / TA.'H'))) ;
  362. *
  363.   EVP1 = EVP01 + ((1. - TETA) * EVP11) + (TETA * EVP21) + EVPC1 ; 
  364.   EVP2 = EVP02 + ((1. - TETA) * EVP12) + (TETA * EVP22) + EVPC2 ; 
  365.   EVP3 = EVP03 + ((1. - TETA) * EVP13) + (TETA * EVP23) + EVPC3 ;
  366. * -------------------------------------
  367. * VITESSE DE DEFORMATION ELASTIQUE EVEL
  368. * ------------------------------------- 
  369.   DSMZZ = (SMZZ - SMZZ1) / (T2 - T1) ;
  370.   EVEL  = DSMZZ / ('MAXI' EE)        ;
  371. * ----------------------------------
  372. * VITESSE DE DEFORMATION IRRADIATION
  373. * ----------------------------------
  374.   EVIR1 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) * 
  375.        ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM1 ;
  376.   EVIR2 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) * 
  377.        ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM2 ;
  378.   EVIR3 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) * 
  379.        ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM3 ;
  380. * ------------------------------------
  381. * VITESSE DE DEFORMATION DENSIFICATION
  382. * ------------------------------------
  383. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT
  384.   BU = AKBU * T2 ;
  385.   'SI' (BU '>' BUMI) ;
  386.     NUMER = (RHO0 * RHO0 * AKEVD) - (RHO0 * AAAA * AKEVD) +
  387.             (RHO0 * AAAA * AKEVD * ('LOG' (BU / BUMI))) -
  388.             (RHO0 * AAAA / T2) ; 
  389.     DENOM = (RHO0 + (AAAA * ('LOG' (BU / BUMI)))) * 
  390.             (RHO0 + (AAAA * ('LOG' (BU / BUMI)))) ;
  391.     EVD = -1. * (AKEVD - (NUMER/DENOM)) / 3. ;
  392.   'SINO' ;
  393.     EVD = 0. ;
  394.   'FINS' ;
  395. * --------------------------------- 
  396. * VITESSE DE DEFORMATION GONFLEMENT
  397. * ---------------------------------
  398. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT
  399.   EVG = AKEVD/3. ;
  400. *
  401.   EVCALC = ((1.+(TA.'K'*DF1))*EVP3) + EVEL + EVIR3 + 
  402.                        EVD + EVG ;
  403.   ECART1 = (EVCALC - EVTH) / EVTH ;
  404.   ERMAX1 = 'MAXI' ('PROG' ('ABS' ECART1) ERMAX1) ;
  405. * -------------------
  406. * VITESSE DE POROSITE
  407. * -------------------
  408.   FF = 'MAXI' ('EXCO' (VV . &BLOC) 'PORO') ;
  409.   DFF = (FF - FF1) / (T2 - T1) ; 
  410.   EVFF = (3.*EVD) + 
  411.          ((1. - FF) * (1.+(TA.'K'*DF1))*(EVP1+EVP2+EVP3)) +
  412.          ((1. - FF) * (EVIR1 + EVIR2 + EVIR3)) ;
  413.   ECART2 = ('ABS' (DFF - EVFF)) / FF ; 
  414.   ERMAX2 = 'MAXI' ('PROG' ('ABS' ECART2) ERMAX2) ;
  415. *
  416.   LTT     = LTT     'ET' ('PROG' T2)     ;
  417.   LEVP3   = LEVP3   'ET' ('PROG' EVP3)   ;
  418.   LEVCALC = LEVCALC 'ET' ('PROG' EVCALC) ;
  419.   LEVTH   = LEVTH   'ET' ('PROG' EVTH)   ;
  420.   LDFF    = LDFF    'ET' ('PROG' DFF)    ;
  421.   LEVFF   = LEVFF   'ET' ('PROG' EVFF)   ;
  422.   LSMZZ   = LSMZZ   'ET' ('PROG' SMZZ)   ;
  423. *
  424.   FF1   = FF ;      
  425.   SMZZ1 = SMZZ ;
  426.   T1 = T2 ; 
  427. *
  428. 'FIN' BLOC ;
  429. *
  430. 'SI' ( ERMAX1 '&lt;EG' 0.05 );
  431.    'ERRE' 0;
  432. 'SINO';
  433.    'MESS' 'ERREUR MAXIMALE VITESSE DEFORMATION :' 
  434.            ERMAX1 '> 0.05 ' ; 
  435.    'ERRE' 5;
  436. 'FINS';
  437. *
  438. 'SI' ( ERMAX2 '&lt;EG' 0.05 );
  439.    'ERRE' 0;
  440. 'SINO';
  441.    'MESS' 'ERREUR MAXIMALE POROSITE :' 
  442.            ERMAX2 '> 0.05 ' ; 
  443.    'ERRE' 5;
  444. 'FINS';
  445. *
  446. 'SI' GRAPH ;
  447.   'TITR' 'Axisymetrique : Vitesses de Deformation'      ;
  448.   EVEVP  = 'EVOL' 'MANU' LTT 'Temps'  LEVP3   'Vitesse' ;
  449.   EVEVC  = 'EVOL' 'MANU' LTT 'Temps'  LEVCALC 'Vitesse' ;
  450.   EVEVTH = 'EVOL' 'MANU' LTT 'Temps'  LEVTH   'Vitesse' ;
  451.   TAD     = TABLE                                       ;
  452.   TAD . 1 = 'TIRC'                                      ;
  453.   TAD . 3 = 'TIRR'                                      ; 
  454.   TAD . 'TITRE' = TABLE                                 ;
  455.   TAD . 'TITRE' . 1 = 'Viscoplastique Analytique'       ;
  456.   TAD . 'TITRE' . 2 = 'Totale Calculee'                 ;
  457.   TAD . 'TITRE' . 3 = 'Totale Analytique'               ;
  458.   'DESS' (EVEVP 'ET' EVEVC 'ET' EVEVTH) TAD 'LEGE'      ;
  459. *
  460.   'TITR' 'Axisymetrique : Vitesses de variation de la porosite' ;
  461.   EVDFF  = 'EVOL' 'MANU' LTT 'Temps' LDFF 'Vitesse'             ;
  462.   EVEVFF = 'EVOL' 'MANU' LTT 'Temps' LEVFF 'Vitesse'            ;
  463.   TAF     = TABLE                                               ;
  464.   TAF . 2 = 'TIRR'                                              ; 
  465.   TAF . 'TITRE' = TABLE                                         ;
  466.   TAF . 'TITRE' . 1 = 'Calculee'                                ;
  467.   TAF . 'TITRE' . 2 =  'Analytique'                             ;
  468.   'DESS' (EVDFF 'ET' EVEVFF) TAF 'LEGE'                         ;
  469. *
  470.   'TITR' 'Axisymetrique : Contrainte'                      ;
  471.   EVSMYY = 'EVOL' 'MANU' LTT 'Temps'   LSMZZ 'Contrainte'  ;
  472.   'DESS' EVSMYY                                            ;
  473. 'FINS' ;
  474. *
  475. 'FIN';
  476.  
  477.  
  478.  
  479.  
  480.  
  481.  
  482.  
  483.  
  484.  
  485.  
  486.  
  487.  
  488.  
  489.  
  490.  
  491.  
  492.  

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