Télécharger gatt_cp.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gatt_cp.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Test gatt_cp.dgibi: Jeux de données
  6. * -----------------------------------
  7. *
  8. 'OPTI' 'DIME' 2 'MODE' 'PLAN' 'CONT' ;
  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 STANDARD COMPRESSIBLE AVEC COUPLAGE STATIQUE *
  19. * *
  20. * MAILLAGE: *
  21. * EPROUVETTE RECTANGULAIRE *
  22. * *
  23. * CHARGEMENT: *
  24. * DEPLACEMENT VERTICAL IMPOSE *
  25. * TRACTION - COMPRESSION *
  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. NA = 1 ;
  39. NR = 1 ;
  40. H = 6.2e-3 ;
  41. D = 5.e-3 ;
  42. R = D/2. ;
  43. HS2 = H/2. ;
  44. *
  45. OO = 0. 0. ;
  46. A1 = R 0. ;
  47. A2 = R HS2 ;
  48. A3 = 0. HS2 ;
  49. *
  50. LB = 'DROIT' NR OO A1 ;
  51. LD = 'DROIT' NA A1 A2 ;
  52. LH = 'DROIT' NR A2 A3 ;
  53. LG = 'DROIT' NA A3 OO ;
  54. *
  55. SU1 = 'DALL' LB LD LH LG 'PLAN' ;
  56. *
  57. MODL1= MODE SU1 MECANIQUE ELASTIQUE
  58. 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' = 1E3 * TA.'KPUO' ;
  70. *
  71. MATREE = MATE MODL1 'YOUN' (TA.'YOUN') ;
  72. PP = 'MANU' 'CHML' MODL1 'T' TE1 'PORO' PO RIGIDITE ;
  73. EE = 'VARI' 'NUAG' MODL1 MATREE PP ;
  74. *
  75. MATR11 = 'MATE' MODL1 'YOUN' EE 'NU' 0.3 'RHO' 10950. 'ALPH' 0. 'TALP' 0. 'TREF' 1.;
  76. MATR12 = 'MATE' MODL1
  77. 'R' (TA.'R') 'DG0' (TA.'DG0') 'DG' (TA.'DG')
  78. 'K1' (TA.'K1') 'M1' (TA.'M1') 'Q1' (TA.'Q1') 'N1' (TA.'N1')
  79. 'K2' (TA.'K2') 'M2' (TA.'M2') 'Q2' (TA.'Q2') 'N2' (TA.'N2')
  80. 'OMEG' (TA.'OMEG') 'H' (TA.'H') 'Q' (TA.'Q') 'BETA' (TA.'BETA')
  81. 'K' (TA.'K') 'A' (TA.'A') 'Q3' (TA.'Q3') 'N3' (TA.'N3')
  82. 'CR' (TA.'CR') 'CR1' (TA.'CR1') 'CR2' (TA.'CR2') 'CR3' (TA.'CR3');
  83. *
  84. MATR13 = 'MATE' MODL1 'KP' (TA.'KPUO') 'AP' (TA.'AP')
  85. 'BP' (TA.'BP') 'QP' (TA.'QP') ;
  86. *
  87. MATR14 = 'MATE' MODL1
  88. 'ADEN' (TA.'ADEN') 'KGON' (TA.'KGON')
  89. 'POR0' (TA.'POR0') 'BUMI' (TA.'BUMI') 'EFIS' (TA.'EFIS') ;
  90. *
  91. * TYPE = 0. combustible UO2 sinon combustible AFA3GLAA
  92. * COMP = 0. combustible incompressible sinon compressible
  93. * DYN = 0. couplage dynamique sinon statique
  94. *
  95. MATR1 = MATR11 'ET' MATR12 'ET' MATR13 'ET' MATR14 ;
  96. *
  97. * -----------------------------------------------------*
  98. * Conditions aux limites + deplacement impose sur CLYH *
  99. * -----------------------------------------------------*
  100. CLYB = 'BLOQ' 'UY' LB ;
  101. CLYH = 'BLOQ' 'UY' LH ;
  102. CLXG = 'BLOQ' 'UX' LG ;
  103. CLT = CLYB 'ET' CLYH 'ET' CLXG ;
  104. *
  105. TMIL = 1.E6 ;
  106. TFIN = 2.E6 ;
  107. *
  108. RHO0 = 1. - TA.'POR0' ;
  109. AKBU = TA.'EFIS'*270./238./10950./RHO0*DF1 ;
  110. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT
  111. TBU = (TA.'BUMI' / AKBU) - 1. ;
  112. TTBU = ('ENTI' (TBU/10000.)) * 10000. ;
  113. *-------------------*
  114. * Instants calcules *
  115. *-------------------*
  116. 'SI' LISTCOUR ;
  117. *--------------
  118. LIST1 = 'PROG' 0 'PAS' 50 300 'PAS' 100 1000 'PAS' 1000 10000
  119. 'PAS' 10000. TTBU
  120. TBU (TTBU+10000.) (TTBU+20000.)
  121. 'PAS' 10000. TMIL ;
  122. 'SINON' ;
  123. *--------
  124. LIST1 = 'PROG' 0 'PAS' 50 300 'PAS' 100 1000 'PAS' 1000 10000
  125. 'PAS' 10000. TTBU
  126. TBU (TTBU+10000.) (TTBU+20000.)
  127. 'PAS' 10000. TMIL
  128. 'PAS' 100. (TMIL + ((TFIN - TMIL) / 200.))
  129. 'PAS' 1000. (TMIL + ((TFIN - TMIL) / 20.))
  130. 'PAS' 2000. (TMIL + ((TFIN - TMIL) / 4.))
  131. 'PAS' 10000. TFIN ;
  132. 'FINSI' ;
  133. *--------
  134. * --------------------------*
  135. * Chargement en temperature *
  136. * --------------------------*
  137. CHTEMP = 'MANU' 'CHPO' SU1 1 'T' 1. ;
  138. EVT = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  139. ('PROG' TE1 TE1) ;
  140. CHARTEMP = 'CHAR' 'T' CHTEMP EVT ;
  141. * ---------------------------------*
  142. * Chargement en densite de fission *
  143. * ---------------------------------*
  144. CHFISS = 'MANU' 'CHPO' SU1 1 'DFIS' 1. ;
  145. EVF = 'EVOL' 'MANU' ('PROG' 0. (2.*TFIN))
  146. ('PROG' DF1 DF1) ;
  147. CHARFISS = 'CHAR' 'DFIS' CHFISS EVF ;
  148. * ---------------------------------*
  149. * Chargement en deplacement impose *
  150. * ---------------------------------*
  151. DEPZ = 'DEPI' CLYH 1 ;
  152. DELTAZ = HS2/200. ;
  153. EVEPS = 'EVOL' 'MANU'
  154. ('PROG' 0. TMIL TFIN (2.*TFIN))
  155. ('PROG' 0. DELTAZ 0. 0.) ;
  156. CHARMECA = 'CHAR' 'DIMP' DEPZ EVEPS ;
  157. *
  158. CHARTOT = CHARMECA 'ET' CHARTEMP 'ET' CHARFISS ;
  159. * ---------------------------------------*
  160. * Variables internes initiales 'PORO'=PO *
  161. * ---------------------------------------*
  162. VAR00 = 'ZERO' MODL1 'VARINTER' ;
  163. VAR01 = 'MANU' 'CHML' MODL1 'PORO' PO
  164. 'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  165. VAR0 = VAR00 + VAR01 ;
  166. *
  167. *
  168. TAB1 = TABLE ;
  169. TAB1.'TEMPERATURES' = TABLE ;
  170. TAB1.'VARIABLES_INTERNES'= TABLE ;
  171. TAB1.'BLOCAGES_MECANIQUES' = CLT ;
  172. TAB1.'MODELE' = MODL1 ;
  173. TAB1.'CHARGEMENT' = CHARTOT ;
  174. TAB1.'TEMPERATURES' . 0 = CHTEMP ;
  175. TAB1.'VARIABLES_INTERNES' . 0 = VAR0 ;
  176. TAB1.'CARACTERISTIQUES' = MATR1 ;
  177. TAB1.'TEMPS_CALCULES' = LIST1 ;
  178. TAB1.'TEMPS_SAUVES' = LIST1 ;
  179. *
  180. TMASAU = TABLE ;
  181. tab1 . 'MES_SAUVEGARDES' = TMASAU ;
  182. TMASAU .'DEFTO' = VRAI ;
  183. TMASAU .'DEFIN' = VRAI ;
  184. *
  185. PASAPAS TAB1 ;
  186. *
  187. *
  188. AP = TA.'AP' ;
  189. BP = TA.'BP' ;
  190. N1 = TA.'N1' ;
  191. N2 = TA.'N2' ;
  192.  
  193. * Facteurs multiplicatifs dus au dopage par le Chrome
  194. WC1 = 'TANH' ( (TA.'CR' - TA.'CR2') / TA.'CR3' ) ;
  195. WC1 = 1. + ( 0.5 * TA.'CR1' * (1. + WC1) ) ;
  196. CV = 180. / PI ;
  197. WC2 = 1. - ( 'COS' (CV * TA.'DG' / TA.'DG0') ) ;
  198. WC2 = 2. * (TA.'DG0'**TA.'M2') * WC2 ;
  199.  
  200. BUMI = TA.'BUMI' ;
  201. AKEVD = AKBU*TA.'KGON' ;
  202. ADEN = TA.'ADEN' ;
  203. KGON = TA.'KGON' ;
  204. * -------------------------------------------------------------------*
  205. * Calcul de la cte AAAA intervenant ds le calcul de la def. de dens. *
  206. * -------------------------------------------------------------------*
  207. BUMAX0=60.D0*BUMI ;
  208. CRIT=1.D-10 ;
  209. 'REPE' BLOC 100 ;
  210. BUMAX = BUMI* ('EXP' (1. - (ADEN/(KGON*BUMAX0)))) ;
  211. BUMAX = (0.2*BUMAX) + (0.8*BUMAX0) ;
  212. TEST='ABS' ((BUMAX-BUMAX0)/BUMAX0) ;
  213. 'SI' ('<' TEST CRIT) ;
  214. 'QUIT' BLOC ;
  215. 'FINS' ;
  216. BUMAX0=BUMAX ;
  217. 'FIN' BLOC ;
  218. AAAA = (-1.D0*RHO0*(ADEN-(KGON*BUMAX))) /
  219. ((1.D0+ADEN)*(LOG(BUMAX/BUMI))) ;
  220. *
  221. * -----------------------*
  222. * Controle des resultats *
  223. * -----------------------*
  224. SS = TAB1 . 'CONTRAINTES' ;
  225. VV = TAB1 . 'VARIABLES_INTERNES' ;
  226. IN = TAB1 . 'DEFORMATIONS_INELASTIQUES' ;
  227. NCONT = ('DIME' (TAB1 . 'CONTRAINTES')) - 1 ;
  228. *
  229. T1 = 0. ;
  230. SMYY1 = 0. ;
  231. FF1 = TA.'POR0' ;
  232. ERMAX1 = 0. ;
  233. ERMAX2 = 0. ;
  234. *
  235. LTT = 'PROG' ;
  236. LDFF = 'PROG' ;
  237. LEVFF = 'PROG' ;
  238. LEVP2 = 'PROG' ;
  239. LEVTH = 'PROG' ;
  240. LEVCALC = 'PROG' ;
  241. LSMYY = 'PROG' ;
  242. *
  243. 'REPE' BLOC NCONT ;
  244. *
  245. * -------------------------------------
  246. * VITESSE DE DEFORMATION THEORIQUE EVTH
  247. * -------------------------------------
  248. T2 = 'EXTR' (TAB1.'TEMPS_SAUVES') (&BLOC + 1);
  249. EVTH = DELTAZ / HS2 / TMIL ;
  250. 'SI' (> T2 (TMIL+0.1)) ;
  251. EVTH = -1. * DELTAZ / HS2 / (TFIN - TMIL) ;
  252. 'FINS' ;
  253. * ----------------------------------------------------
  254. * SMYY EST PREPONDERANTE DEVANT LES AUTRES CONTRAINTES
  255. * ----------------------------------------------------
  256. SMYY = 'MAXI' ('EXCO' (SS . &BLOC) 'SMYY') ;
  257. SSM = SMYY / 3. ;
  258. SSEQ = 'ABS' SMYY ;
  259. * Deviateur des contraintes
  260. SSPRIM1 = -1. * SMYY / 3. ;
  261. SSPRIM2 = 2. * SMYY / 3. ;
  262. SSPRIM3 = -1. * SMYY / 3. ;
  263. * ---------------------------------------------
  264. * DEFORMATION VISCO-PLASTIQUE EQUIVALENTE EVPEQ
  265. * ---------------------------------------------
  266. EVXX = 'MAXI' ('EXCO' IN.&BLOC 'EIXX') ;
  267. EVYY = 'MAXI' ('EXCO' IN.&BLOC 'EIYY') ;
  268. EVZZ = 'MAXI' ('EXCO' IN.&BLOC 'EIZZ') ;
  269. EVH = (EVXX+EVYY+EVZZ) / 3. ;
  270. EVPEQ = ((EVXX - EVH)**2) + ((EVYY - EVH)**2) + ((EVZZ - EVH)**2) ;
  271. EVPEQ = ((2./3.) * EVPEQ)**0.5 ;
  272. * ------------------------------------------
  273. * VITESSE DE DEFORMATION VISCO-PLASTIQUE EVP
  274. * ------------------------------------------
  275. *-- fluage primaire
  276. *
  277. 'SI' ('<' EVPEQ 1E-10) ; EVPEQ = 1E-10 ; 'FINSI' ;
  278. *
  279. EVP01 = 1.5 * TA.'KPUO' *
  280. ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) *
  281. ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM1 ;
  282. *
  283. EVP02 = 1.5 * TA.'KPUO' *
  284. ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) *
  285. ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM2 ;
  286. *
  287. EVP03 = 1.5 * TA.'KPUO' *
  288. ('EXP' (-1.*TA.'QP'/(TA.'R'*TE1))) *
  289. ( (SSEQ**(AP-1)) / (EVPEQ**BP) ) * SSPRIM3 ;
  290. *
  291. *-- fluage secondaire (2 mecanismes)
  292. *
  293. PP = 'MAXI' ('EXCO' VV . &BLOC 'PORO') ;
  294. *
  295. A1 = (N1* (PP**(-1./N1) - 1.)**(-2.*N1/(N1+1.))) ;
  296. B1 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N1/(N1+1.))) ;
  297. AB1 = (9. * A1) / ( (4.*B1) + A1 ) ;
  298. BA1 = B1 / ( B1 + (A1/4.) ) ;
  299. *
  300. GSSEQ = ((BA1 * (SSEQ**2)) + (AB1 * (SSM**2))) ** 0.5 ;
  301. *
  302. A2 = (N2* (PP**(-1./N2) - 1.)**(-2.*N2/(N2+1.))) ;
  303. B2 = (1. + (2.*PP/3.))/ ((1. - PP) **(2.*N2/(N2+1.))) ;
  304. *
  305. EVP11 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  306. ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) *
  307. (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  308. ((A1*1.5*SSM) + (3.*B1*SSPRIM1))) ;
  309. *
  310. EVP21 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  311. (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  312. ((A2*1.5*SSM) + (3.*B2*SSPRIM1))) ;
  313. *
  314. EVP12 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  315. ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) *
  316. (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  317. ((A1*1.5*SSM) + (3.*B1*SSPRIM2))) ;
  318. *
  319. EVP22 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  320. (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  321. ((A2*1.5*SSM) + (3.*B2*SSPRIM2))) ;
  322. *
  323. EVP13 = 0.5 * (TA.'K1' * (TA.'DG'**TA.'M1') * WC1 *
  324. ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) *
  325. (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1-1.)/2.)) *
  326. ((A1*1.5*SSM) + (3.*B1*SSPRIM3))) ;
  327. *
  328. EVP23 = 0.5 * (TA.'K2' * WC2 * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  329. (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2-1.)/2.)) *
  330. ((A2*1.5*SSM) + (3.*B2*SSPRIM3))) ;
  331. *
  332. PSI1 = TA.'K1' * WC1 * (TA.'DG'**TA.'M1')/(N1+1.) *
  333. ('EXP' (-1.*TA.'Q1'/(TA.'R'*TE1))) *
  334. (((A1*((1.5*SSM)**2)) + (B1*(SSEQ ** 2)))**((N1+1.)/2.)) ;
  335. *
  336. PSI2 = TA.'K2' * WC2/(N2+1.) * ('EXP' (-1.*TA.'Q2'/(TA.'R'*TE1))) *
  337. (((A2*((1.5*SSM)**2)) + (B2*(SSEQ ** 2)))**((N2+1.)/2.)) ;
  338. *
  339. EVPC1 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  340. (GSSEQ**(-1.*TA.'Q' - 2.)) *
  341. (1. -
  342. (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q')))) / TA.'H'))**2))
  343. * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) + (1.5*BA1*SSPRIM1) ) ;
  344. *
  345. EVPC2 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  346. (GSSEQ**(-1.*TA.'Q' - 2.)) *
  347. (1. -
  348. (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q')))) / TA.'H'))**2))
  349. * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) + (1.5*BA1*SSPRIM2) ) ;
  350. *
  351. EVPC3 = 0.5 * TA.'BETA' / TA.'H' * TA.'Q' * TA.'OMEG' *
  352. (GSSEQ**(-1.*TA.'Q' - 2.)) *
  353. (1. -
  354. (('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q')))) / TA.'H'))**2))
  355. * (PSI2 - PSI1) * ( ((AB1*SSM)/3.) + (1.5*BA1*SSPRIM3) ) ;
  356. *
  357. TETA = 0.5 * TA.'BETA' *
  358. (1. +
  359. ('TANH' ((TE1 - (TA.'OMEG' * (GSSEQ**(-1.*TA.'Q')))) / TA.'H'))) ;
  360. *
  361. EVP1 = EVP01 + ((1. - TETA) * EVP11) + (TETA * EVP21) + EVPC1 ;
  362. EVP2 = EVP02 + ((1. - TETA) * EVP12) + (TETA * EVP22) + EVPC2 ;
  363. EVP3 = EVP03 + ((1. - TETA) * EVP13) + (TETA * EVP23) + EVPC3 ;
  364. * -------------------------------------
  365. * VITESSE DE DEFORMATION ELASTIQUE EVEL
  366. * -------------------------------------
  367. DSMYY = (SMYY - SMYY1) / (T2 - T1) ;
  368. EVEL = DSMYY / ('MAXI' EE) ;
  369. * ----------------------------------
  370. * VITESSE DE DEFORMATION IRRADIATION
  371. * ----------------------------------
  372. EVIR1 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) *
  373. ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM1 ;
  374. EVIR2 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) *
  375. ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM2 ;
  376. EVIR3 = 1.5 * TA.'A' * DF1 * (SSEQ **(TA.'N3' - 1.)) *
  377. ('EXP' (-1.*TA.'Q3'/(TA.'R'*TE1))) * SSPRIM3 ;
  378. * ------------------------------------
  379. * VITESSE DE DEFORMATION DENSIFICATION
  380. * ------------------------------------
  381. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT
  382. BU = AKBU * T2 ;
  383. 'SI' (BU '>' BUMI) ;
  384. NUMER = (RHO0 * RHO0 * AKEVD) - (RHO0 * AAAA * AKEVD) +
  385. (RHO0 * AAAA * AKEVD * ('LOG' (BU / BUMI))) -
  386. (RHO0 * AAAA / T2) ;
  387. DENOM = (RHO0 + (AAAA * ('LOG' (BU / BUMI)))) *
  388. (RHO0 + (AAAA * ('LOG' (BU / BUMI)))) ;
  389. EVD = -1. * (AKEVD - (NUMER/DENOM)) / 3. ;
  390. 'SINO' ;
  391. EVD = 0. ;
  392. 'FINS' ;
  393. * ---------------------------------
  394. * VITESSE DE DEFORMATION GONFLEMENT
  395. * ---------------------------------
  396. * CE QUI SUIT EST VRAI CAR LE FLUX DE FISSIONS EST CONSTANT
  397. EVG = AKEVD/3. ;
  398. * -------------------------------
  399. * VITESSE DE DEFORMATION SELON YY
  400. * -------------------------------
  401. EVCALC = ((1.+(TA.'K'*DF1))*EVP2) + EVEL + EVIR2 +
  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. LEVP2 = LEVP2 'ET' ('PROG' EVP2) ;
  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. LSMYY = LSMYY 'ET' ('PROG' SMYY) ;
  423. *
  424. FF1 = FF ;
  425. SMYY1 = SMYY ;
  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' 'Contraintes planes : Vitesses de Deformation' ;
  448. EVEVP = 'EVOL' 'MANU' LTT 'Temps' LEVP2 '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' 'Contraintes planes : 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' 'Contraintes planes : Contraintes' ;
  471. EVSMYY = 'EVOL' 'MANU' LTT 'Temps' LSMYY 'Contrainte' ;
  472. 'DESS' EVSMYY ;
  473. 'FINS' ;
  474. *
  475. 'FIN';
  476.  
  477.  
  478.  
  479.  
  480.  
  481.  
  482.  
  483.  
  484.  
  485.  
  486.  
  487.  
  488.  
  489.  
  490.  

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