Télécharger uo2_cas2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : uo2_cas2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Test uo_cas2.dgibi: Jeux de donnees
  6. * ------------------------------------
  7. *
  8. 'OPTI' 'DIME' 3 'MODE' 'TRID' ;
  9. 'OPTI' 'ELEM' 'CUB8' ;
  10. 'OPTI' 'TRAC' 'PSC' ;
  11. 'TEMPS' 'ZERO' ;
  12. *
  13. **********************************************************
  14. * TEST DE VALIDATION *
  15. * ------------------ *
  16. * MODELE UO2 *
  17. * (COUPLAGE DES MODELES GATT_MONERIE ET OTTOSEN) *
  18. * UO2 STANDARD COMPRESSIBLE AVEC COUPLAGE STATIQUE *
  19. * *
  20. * MAILLAGE: *
  21. * CUBE *
  22. * *
  23. * CHARGEMENT: *
  24. * DEPLACEMENT VERTICAL IMPOSE EN TRACTION AVEC *
  25. * VITESSE CTE PAR MORCEAU (POSITIVE PUIS NEGATIVE *
  26. * PUIS POSITIVE) *
  27. * TEMPERATURE LINEAIRE CROISSANTE *
  28. * DENSITE DE FISSIONS LINEAIRE CROISSANTE *
  29. * OTTOSEN ACTIVE (TRACTION ELASTIQUE - FISSURATION - *
  30. * OUVERTURE - FERMETURE - COMPRESSION - REOUVERTURE -*
  31. * POURSUITE DE FISSURATION VERS RUPTURE) *
  32. * GONFLEMENT ACTIVE *
  33. * PAS DE VISCOPLASTICITE (GATT_MONERIE) *
  34. * PAS DE DENSIFICATION *
  35. **********************************************************
  36. *
  37. DESSI = FAUX ;
  38. *
  39. * -----------------------------*
  40. * MAILLAGE D UN CUBE DE COTE 1 *
  41. * -----------------------------*
  42. HS2 = 2E-5 ;
  43. *
  44. OO = 0. 0. 0. ;
  45. A1 = HS2 0. 0. ;
  46. A2 = HS2 0. HS2 ;
  47. A3 = 0. 0. HS2 ;
  48. *
  49. LB = 'DROIT' 1 OO A1 ;
  50. LD = 'DROIT' 1 A1 A2 ;
  51. LH = 'DROIT' 1 A2 A3 ;
  52. LG = 'DROIT' 1 A3 OO ;
  53. *
  54. SUA = 'DALL' LB LD LH LG 'PLAN' ;
  55. SU2 = SUA 'PLUS' (0. HS2 0.) ;
  56. VOL1 = 'VOLU' 1 SUA SU2 ;
  57. 'ELIM' (HS2/1000.) VOL1 ;
  58. *
  59. SUB = 'POIN' ('COOR' 3 VOL1) 'INFERIEUR' (HS2/1000.) ;
  60. SUH = 'POIN' ('COOR' 3 VOL1) 'SUPERIEUR' (HS2 - (HS2/1000.));
  61. LG = 'POIN' ('COOR' 1 VOL1) 'INFERIEUR' (HS2/1000.) ;
  62. * -----------------------------------------------------*
  63. * Conditions aux limites + deplacement impose sur CLZH *
  64. * -----------------------------------------------------*
  65. CLZB = 'BLOQ' UZ SUB ;
  66. CLZH = 'BLOQ' UZ SUH ;
  67. CLXG = 'BLOQ' UX LG ;
  68. CLYG = 'BLOQ' UY SUA ;
  69. CLT = CLZB 'ET' CLZH 'ET' CLXG 'ET' CLYG ;
  70. * -------*
  71. * Modele *
  72. * -------*
  73. MODL1= MODE VOL1 MECANIQUE ELASTIQUE VISCOPLASTIQUE UO2 ;
  74. *
  75. * Temperature moyenne
  76. * -------------------
  77. TE1 = 1420. + 273. ;
  78. *
  79. * Porosite initiale (cte au cours du calcul)
  80. * ------------------------------------------
  81. PO = 4.52E-2 ;
  82.  
  83. TA = @GATTPAR '/u2/castem/divers/fichier_gatt' ;
  84. *
  85. PP = 'MANU' 'CHML' MODL1 'T' TE1 'PORO' PO RIGIDITE ;
  86. *
  87. *
  88. * --------------------------------------*
  89. * 'YOUN' constant ou fonction de T et f *
  90. * --------------------------------------*
  91. MATREE = MATE MODL1 'YOUN' (TA.'YOUN') ;
  92. EE = 'VARI' 'NUAG' MODL1 MATREE PP ;
  93. * EE = (TA.'YOUN') ;
  94. MATREE = MATE MODL1 'YOUN' EE ;
  95. * ------------------------------------*
  96. * 'NU' constant ou fonction de T et f *
  97. * ------------------------------------*
  98. MATRUU = MATE MODL1 'NU' (TA.'NU') ;
  99. UU = 'VARI' 'NUAG' MODL1 MATRUU PP ;
  100. * UU = (TA.'NU') ;
  101. * --------------------------------------*
  102. * 'ALPH' constant *
  103. * --------------------------------------*
  104. XALFA = 1.319E-5 ;
  105. * -----------------------------------------------------------------*
  106. * Temperature de reference et temperature de reference pour 'ALPH' *
  107. * -----------------------------------------------------------------*
  108. TREF = TE1 ;
  109. TALPHA = 0. ;
  110. *TALPHA = TE1 ;
  111. *
  112. * -----------------------------*
  113. * Materiau partie GATT_MONERIE *
  114. * -----------------------------*
  115. * fluage primaire
  116. TA.'KPUO' = 0. ;
  117. * fluage secondaire (2 mecanismes)
  118. TA.'K1' = 0. ;
  119. TA.'K2' = 0. ;
  120. * fluage d"irradiation
  121. TA.'A' = 0. ;
  122. * densification/gonflement
  123. TA.'BUMI' = 5E3 * (TA.'BUMI') ;
  124. *
  125. MATR11 = 'MATE' MODL1 'YOUN' EE 'NU' UU 'RHO' (TA.'RHO') 'ALPH' XALFA ;
  126. *
  127. MATR12 = 'MATE' MODL1 'R' (TA.'R') 'DG0' (TA.'DG0') 'DG' (TA.'DG') 'K1' (TA.'K1') 'M1' (TA.'M1') 'Q1' (TA.'Q1') 'N1' (TA.'N1') 'K2' (TA.'K2') 'M2' (TA.'M2') 'Q2' (TA.'Q2') 'N2' (TA.'N2') 'OMEG' (TA.'OMEG') 'H' (TA.'H') 'Q' (TA.'Q') 'BETA' (TA.'BETA') 'K' (TA.'K') 'A' (TA.'A') 'Q3' (TA.'Q3') 'N3' (TA.'N3') 'CR' (TA.'CR') 'CR1' (TA.'CR1') 'CR2' (TA.'CR2') 'CR3' (TA.'CR3');
  128. *
  129. MATR13 = 'MATE' MODL1 'KP' (TA.'KPUO') 'AP' (TA.'AP') 'BP' (TA.'BP') 'QP' (TA.'QP') ;
  130. *
  131. MATR14 = 'MATE' MODL1 'ADEN' (TA.'ADEN') 'KGON' (TA.'KGON') 'POR0' (TA.'POR0') 'BUMI' (TA.'BUMI') 'EFIS' (TA.'EFIS') ;
  132. *
  133. * TYPE = 0. combustible UO2 sinon combustible AFA3GLAA
  134. * COMP = 0. combustible compressible sinon incompressible
  135. * DYN = 0. couplage statique sinon dynamique
  136. *
  137. * -----------------------------*
  138. * Materiau partie OTTOSEN *
  139. * -----------------------------*
  140. MATR16 = 'TAIL' MODL1 ;
  141. *
  142. * Resistance au cisaillement 'GS ' par defaut: 'YOUN'*1.8E-4
  143. * Ouverture a rupture 'WRUP' = 0 --> 'WRUP' = 2*'GFTR'/'LTR '
  144. * Relation bilineaire 'BILI' = 0 --> 'BILI' = 0 (ouverture)
  145. *
  146. XGFTR = 1.8 ; XGFTR1=XGFTR ; XGFTR2=XGFTR ; XGFTR3=XGFTR ;
  147. XLTR = 150E6 ; XLTR1 =XLTR ; XLTR2 =XLTR ; XLTR3 =XLTR ;
  148. XBTR = 0.2 ;
  149. *
  150. MATR17 = 'MATE' MODL1 'GFTR' XGFTR 'LTR ' XLTR 'BTR ' XBTR 'WRUP' 0. 'BILI' 0. 'GFT1' XGFTR1 'GFT2' XGFTR2 'GFT3' XGFTR3 'GS1 ' 0. 'GS2 ' 0. 'GS3 ' 0. 'LTR1' XLTR1 'LTR2' XLTR2 'LTR3' XLTR3 'WRU1' 0. 'WRU2' 0. 'WRU3' 0. 'BIL1' 0. 'BIL2' 0. 'BIL3' 0. ;
  151. *
  152. MATR1 = MATR11 'ET' MATR12 'ET' MATR13 'ET' MATR14 'ET' MATR16 'ET' MATR17 ;
  153. *
  154. *
  155. RHO0 = 1. - TA.'POR0' ;
  156. TAU0 = TA.'EFIS'*270./238./10950./RHO0 ;
  157. *
  158. *
  159. 'DEBPROC' TRINOME XA*'FLOTTANT' XB*'FLOTTANT' XC*'FLOTTANT' mau*'MOT';
  160. DISCRI = ( (XB**2.) - (4.*XA*XC) )**0.5 ;
  161. RACINE1 = ( (-1.*XB) + DISCRI ) / (2.*XA) ;
  162. RACINE2 = ( (-1.*XB) - DISCRI ) / (2.*XA) ;
  163. LRACINE = 'PROG' RACINE1 RACINE2 ;
  164. PRACINE = RACINE1 * RACINE2 ;
  165. SRACINE = RACINE1 + RACINE2 ;
  166. 'SI' (PRACINE '<' 0.) ;
  167. RACINE = 'MAXI' LRACINE ;
  168. 'SINON' ;
  169. RACINE = 'MINI' LRACINE ;
  170. 'FINSI' ;
  171. 'FINPROC' RACINE ;
  172. *
  173. *
  174. CONVSIG = 1E6 ;
  175. CONVTO = 24.*3600E3 ;
  176. *
  177. TVIT = 1.0E6 ;
  178. UZ_TVIT = HS2/200. ;
  179. VUZ1 = UZ_TVIT/TVIT ;
  180. VUZ2 = -0.5 * VUZ1 ;
  181. VUZ3 = 1.7 * VUZ1 ;
  182. *
  183. TBIF = 0.6E6 ;
  184. TFIN = 1.0E6 ;
  185. TFIN2 = 1.9E6 ;
  186. * ---------------------------------*
  187. * Chargement en deplacement impose *
  188. * ---------------------------------*
  189. DEPZ = 'DEPI' CLZH 1 ;
  190. UZ_TBIF = VUZ1 * TBIF ;
  191. UZ_TFIN = UZ_TBIF + (VUZ2*(TFIN-TBIF)) ;
  192. UZ_TFIN2 = UZ_TFIN + (VUZ3*(TFIN2-TFIN)) ;
  193. 'TITR' 'DEPLACEMENT IMPOSE' ;
  194. EVDEP = 'EVOL' 'MANU' 'Temps (s)' ('PROG' 0. TBIF TFIN TFIN2 ) 'Dep (m)' ('PROG' 0. UZ_TBIF UZ_TFIN UZ_TFIN2) ;
  195. LDEPX = 'EXTR' EVDEP 'ABSC' ;
  196. LDEPY = 'EXTR' EVDEP 'ORDO' ;
  197. si DESSI ; 'DESS' EVDEP ; 'FINSI' ;
  198. CHARMECA = 'CHAR' 'DIMP' DEPZ EVDEP ;
  199. * --------------------------*
  200. * Chargement en temperature *
  201. * --------------------------*
  202. TP1 = 1E-4 ;
  203. TE0 = 1320. + 273. ;
  204. TE2 = (TP1*TFIN2) + TE0 ;
  205. CHTEMP = 'MANU' 'CHPO' VOL1 1 'T' 1. ;
  206. EVT = 'EVOL' 'MANU' ('PROG' 0. TFIN2) ('PROG' TE0 TE2) ;
  207. CHARTEMP = 'CHAR' 'T' CHTEMP EVT ;
  208. * ---------------------------------*
  209. * Chargement en densite de fission *
  210. * ---------------------------------*
  211. PHIP1 = 3E14 ;
  212. PHI2 = PHIP1 * TFIN2 ;
  213. CHFISS = 'MANU' 'CHPO' VOL1 1 'DFIS' 1. ;
  214. EVF = 'EVOL' 'MANU' ('PROG' 0. TFIN2) ('PROG' 0. PHI2) ;
  215. CHARFISS = 'CHAR' 'DFIS' CHFISS EVF ;
  216. *
  217. * ---------------------*
  218. * -------------------- *
  219. * INSTANTS PRIVILEGIES *
  220. * -------------------- *
  221. * ---------------------*
  222. XA = (TA.'KGON' / 6.) * TAU0 * PHIP1 ;
  223. *
  224. LZ = 'MAXI' ('EXCO' MATR16 'LZZ') ;
  225. PENTZ = -1. * ( (XLTR1**2) / (2.*XGFTR1) ) ;
  226. JZ = 1./(LZ*PENTZ) ;
  227. LZLIM = ('MAXI' EE) / (('ABS' PENTZ) * (1.+('MAXI' UU))) ;
  228. WRUPT = (2.*XGFTR1) / XLTR1 ;
  229. ERUPT = WRUPT / LZ ;
  230. WFERM = XBTR * WRUPT ;
  231. EFERM = WFERM / LZ ;
  232. * ------------*
  233. * Fissuration *
  234. * ------------*
  235. ET1 = VUZ1/HS2 ;
  236. ET0 = 0. ;
  237. B0 = ET1 - (XALFA*TP1) ; XB = -1. * B0 ;
  238. XC = (XLTR1/('MAXI' EE)) - ET0 ;
  239. TFISS = TRINOME XA XB XC 'TFISS' ;
  240. * ----------*
  241. * Fermeture *
  242. * ----------*
  243. * Ouverture residuelle de fermeture
  244. ET1 = VUZ2 / HS2 ;
  245. ET0 = UZ_TBIF / HS2 ;
  246. ETH0 = XALFA*TP1*TBIF ;
  247. GF0 = XA*(TBIF**2) ;
  248. FZ0 = ET0-(ETH0+GF0) ;
  249. NUM = FZ0 + (XLTR1*JZ) ;
  250. DEN = (1./('MAXI' EE)) + JZ ;
  251. SMAX = NUM/DEN ;
  252. EMAX = JZ*(SMAX-XLTR1) ;
  253. WMAX = EMAX * LZ ;
  254. WFERM = XBTR * WMAX ;
  255. EFERM = WFERM / LZ ;
  256. JZ_PRIM = (EMAX - EFERM) / SMAX ;
  257. *
  258. B0 = ET1 - (XALFA*TP1) ; XB = -1. * B0 ;
  259. XC = (EFERM + (ET1*TBIF)) - ET0 ;
  260. TFERM = TRINOME XA XB XC 'TFERM' ;
  261. * ------------*
  262. * Reouverture *
  263. * ------------*
  264. ET1 = VUZ3 / HS2 ;
  265. ET0 = 'IPOL' TFERM LDEPX (LDEPY/HS2) ;
  266. ETH0 = XALFA*TP1*TFERM ;
  267. GF0 = XA*(TFERM**2) ;
  268. FZ0 = ET0-(ETH0+GF0) ;
  269. ET0 = UZ_TFIN / HS2 ;
  270. *
  271. B0 = ET1 - (XALFA*TP1) ; XB = -1. * B0 ;
  272. XC = (FZ0 + (ET1*TFIN)) - ET0 ;
  273. TREOUV = TRINOME XA XB XC 'TREOUV' ;
  274. * --------------------------- *
  275. * Poursuite de la fissuration *
  276. * --------------------------- *
  277. ET0 = 'IPOL' TREOUV LDEPX (LDEPY/HS2) ;
  278. ETH0 = XALFA*TP1*TREOUV ;
  279. GF0 = XA*(TREOUV**2) ;
  280. FZ0 = ET0-(ETH0+GF0) ;
  281. ET0 = UZ_TFIN / HS2 ;
  282. *
  283. B0 = ET1 - (XALFA*TP1) ; XB = -1. * B0 ;
  284. XC0 = (1./('MAXI' EE)) + JZ_PRIM ;
  285. XC0 = SMAX * XC0 ;
  286. XC = (FZ0 + XC0 + (ET1*TFIN)) - ET0 ;
  287. TBIF2 = TRINOME XA XB XC 'TBIF2' ;
  288. * -------------------*
  289. * Liste des instants *
  290. * -------------------*
  291. VPAS = TVIT/40. ;
  292. LTIME = ('PROG' 0. 'PAS' VPAS TBIF 'PAS' VPAS TFIN 'PAS' VPAS TFIN2) 'ET' ('PROG' TFERM TREOUV TBIF2) ;
  293. LTIME = 'ORDO' LTIME ;
  294. UN = 'PROG' ('DIME' LTIME)*1 ;
  295. * -------------------*
  296. * Deformation totale *
  297. * -------------------*
  298. EVEPS = EVDEP / HS2 ;
  299. LEPSX = 'EXTR' EVEPS 'ABSC' ;
  300. LEPSY = 'EXTR' EVEPS 'ORDO' ;
  301. LEPS = 'IPOL' LTIME LEPSX LEPSY ;
  302. 'TITR' 'DEFORMATION IMPOSEE' ;
  303. EVEPS = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Def' LEPS ;
  304. 'SI' DESSI ; 'DESS' EVEPS ; 'FINSI' ;
  305. * -------------------------------------*
  306. * Temperature et deformation thermique *
  307. * -------------------------------------*
  308. LT = (TP1*LTIME) + ((TE0-273.)*UN) ;
  309. 'TITR' 'TEMPERATURE' ;
  310. EVT = 'EVOL' 'MANU' 'Temps (s)' LTIME 'T (°C)' LT ;
  311. 'SI' DESSI ; 'DESS' EVT ; 'FINSI' ;
  312. *
  313. LTH = XALFA*TP1*LTIME ;
  314. 'TITR' 'DEFORMATION THERMIQUE' ;
  315. EVDTH = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Def Th' LTH ;
  316. 'SI' DESSI ; 'DESS' EVDTH ; 'FINSI' ;
  317. * -------------------*
  318. * Densite de fission *
  319. * -------------------*
  320. LPHI = PHIP1*LTIME ;
  321. 'TITR' 'DENSITE DE FISSION' ;
  322. EVF = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Phi (SI)' LPHI ;
  323. 'SI' DESSI ; 'DESS' EVF ; 'FINSI' ;
  324. * ------------------------------------------------*
  325. * Taux de combustion et deformation de gonflement *
  326. * ------------------------------------------------*
  327. LTO = (TAU0*PHIP1/2.)*(LTIME**2) ;
  328. 'TITR' 'TAUX DE COMBUSTION' ;
  329. EVTO_SI = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Taux(SI)' LTO ;
  330. 'SI' DESSI ; 'DESS' EVTO_SI ; 'FINSI' ;
  331. EVTO = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Taux(MWj/tU)' (LTO/CONVTO) ;
  332. 'SI' DESSI ; 'DESS' EVTO ; 'FINSI' ;
  333. *
  334. LGF = (TA.'KGON' / 3.) * LTO ;
  335. 'TITR' 'DEFORMATION DE GONFLEMENT' ;
  336. EVGF = 'EVOL' 'MANU' 'Temps (s)' LTIME 'Def Gonf' LGF ;
  337. 'SI' DESSI ; 'DESS' EVGF ; 'FINSI' ;
  338.  
  339. * ***************** *
  340. * AVANT FISSURATION *
  341. * ***************** *
  342. LTIME_AV = LTIME 'MASQ' 'EGINFE' TFISS ;
  343. LFZ = LEPS-(LTH+LGF) ;
  344. *
  345. LSIGZ = ('MAXI' EE) * LFZ ;
  346. SIGF = 'IPOL' TFISS LTIME (LSIGZ/CONVSIG) ;
  347. * ************************************** *
  348. * APRES FISSURATION ET AVANT BIFURCATION *
  349. * ************************************** *
  350. LTIME_AP = LTIME 'MASQ' 'SUPERIEUR' TFISS ;
  351. LJZ = JZ * LTIME_AP ;
  352. NUM = LFZ + (XLTR1*LJZ) ;
  353. LE1 = (1./('MAXI' EE))*UN ;
  354. DEN = LE1 + LJZ ;
  355. LSIGZ = NUM/DEN ;
  356. *
  357. LXLTR1 = XLTR1*UN ;
  358. LWZ = LTIME_AP * ( (LSIGZ-LXLTR1)/PENTZ ) ;
  359. *
  360. LEPSZF = LWZ / LZ ;
  361. *
  362. SIGZB = 'IPOL' TBIF LTIME (LSIGZ/CONVSIG) ;
  363. WB = 'IPOL' TBIF LTIME LWZ ;
  364. * ********************************* *
  365. * APRES BIFURCATION ET EN FERMETURE *
  366. * ********************************* *
  367. LTIME_AV = LTIME 'MASQ' 'EGINFE' TBIF ;
  368. LTIME_AP = LTIME 'MASQ' 'SUPERIEUR' TBIF ;
  369. LJZ_PRIM = JZ_PRIM * LTIME_AP ;
  370. LEFERM = EFERM * LTIME_AP ;
  371. NUM = LFZ - LEFERM ;
  372. DEN = LE1 + LJZ_PRIM ;
  373. LSIGZ = (LSIGZ * LTIME_AV) + ((NUM/DEN) * LTIME_AP) ;
  374. *
  375. LEMAX = EMAX * LTIME_AP ;
  376. LSMAX = SMAX * LTIME_AP ;
  377. LEPSZF = (LEPSZF * LTIME_AV) + LEMAX + (LJZ_PRIM * (LSIGZ - LSMAX)) ;
  378. *
  379. LWZ = LZ * LEPSZF ;
  380. *
  381. SIGZC = 'IPOL' TFERM LTIME (LSIGZ/CONVSIG) ;
  382. WFERMC = 'IPOL' TFERM LTIME LWZ ;
  383. * ****************************** *
  384. * APRES FERMETURE EN COMPRESSION *
  385. * ****************************** *
  386. LTIME_AV = LTIME 'MASQ' 'EGINFE' TFERM ;
  387. LTIME_AP = LTIME 'MASQ' 'SUPERIEUR' TFERM ;
  388. *
  389. ET0 = 'IPOL' TFERM LTIME LEPS ;
  390. ETH0 = XALFA*TP1*TFERM ;
  391. GF0 = XA*(TFERM**2) ;
  392. FZ0 = ET0-(ETH0+GF0) ;
  393. *
  394. LSIGZ_AP = ('MAXI' EE) * (LFZ - (FZ0 * LTIME_AP)) ;
  395. LSIGZ = (LSIGZ * LTIME_AV) + (LSIGZ_AP * LTIME_AP) ;
  396. *
  397. LEPSZF = (LEPSZF * LTIME_AV) + (EFERM * LTIME_AP) ;
  398. *
  399. LWZ = LZ * LEPSZF ;
  400. *
  401. SIGZC = 'IPOL' TREOUV LTIME (LSIGZ/CONVSIG) ;
  402. * ****************************************** *
  403. * REOUVERTURE AVANT POURSUITE DE FISSURATION *
  404. * ****************************************** *
  405. LTIME_AV = LTIME 'MASQ' 'EGINFE' TREOUV ;
  406. LTIME_AP = LTIME 'MASQ' 'SUPERIEUR' TREOUV ;
  407. *
  408. ET0 = 'IPOL' TREOUV LTIME LEPS ;
  409. ETH0 = XALFA*TP1*TREOUV ;
  410. GF0 = XA*(TREOUV**2) ;
  411. FZ0 = ET0-(ETH0+GF0) ;
  412. NUM = LFZ - (FZ0 * LTIME_AP) ;
  413. LJZ_PRIM = JZ_PRIM * LTIME_AP ;
  414. DEN = LE1 + LJZ_PRIM ;
  415. LSIGZ = (LSIGZ * LTIME_AV) + ((NUM/DEN) * LTIME_AP) ;
  416. *
  417. LEMAX = EMAX * LTIME_AP ;
  418. LSMAX = SMAX * LTIME_AP ;
  419. LEPSZF = (LEPSZF * LTIME_AV) + LEMAX + (LJZ_PRIM * (LSIGZ - LSMAX)) ;
  420. *
  421. LWZ = LZ * LEPSZF ;
  422. * ************************ *
  423. * POURSUITE DE FISSURATION *
  424. * ************************ *
  425. LTIME_AV = LTIME 'MASQ' 'EGINFE' TBIF2 ;
  426. LTIME_AP = LTIME 'MASQ' 'SUPERIEUR' TBIF2 ;
  427. *
  428. LJZ = JZ * LTIME_AP ;
  429. NUM = LFZ + (XLTR1*LJZ) ;
  430. DEN = LE1 + LJZ ;
  431. LSIGZ = (LSIGZ * LTIME_AV) + ((NUM/DEN) * LTIME_AP) ;
  432. *
  433. LWZ = (LWZ * LTIME_AV) + (LTIME_AP * ( (LSIGZ-LXLTR1)/PENTZ ));
  434. *
  435. LEPSZF = LWZ / LZ ;
  436. *
  437. SIGZC = 'IPOL' TBIF2 LTIME (LSIGZ/CONVSIG) ;
  438. WC = 'IPOL' TBIF2 LTIME LWZ ;
  439. *
  440. 'TITR' 'CONTRAINTE DANS LA DIRECTION Z' ;
  441. EVSIGZ = 'EVOL' 'ROUG' 'MANU' 'Temps (s)' LTIME 'SIGZ(MPa)' (LSIGZ/CONVSIG) ;
  442. 'SI' DESSI ; 'DESS' EVSIGZ; 'FINSI' ;
  443. 'TITR' 'OUVERTURE DANS LA DIRECTION Z' ;
  444. EVWZ = 'EVOL' 'ROUG' 'MANU' 'Temps (s)' LTIME 'WZ(m)' LWZ ;
  445. 'SI' DESSI ; 'DESS' EVWZ; 'FINSI' ;
  446. 'TITR' 'DEFORMATION DE FISSURATION DANS LA DIRECTION Z' ;
  447. EVEPSFZ = 'EVOL' 'ROUG' 'MANU' 'Temps (s)' LTIME 'EPSFZ' LEPSZF ;
  448. 'SI' DESSI ; 'DESS' EVEPSFZ; 'FINSI' ;
  449. 'TITR' 'RELATION CONTRAINTE/OUVERTURE DANS LA DIRECTION Z' ;
  450. EV_FISZ = 'EVOL' 'ROUG' 'MANU' 'WZ(m)' LWZ 'SIGZ(MPa)' (LSIGZ/CONVSIG);
  451. 'SI' DESSI ; 'DESS' EV_FISZ; 'FINSI' ;
  452. *-------------------*
  453. * Instants calcules *
  454. *-------------------*
  455. LIST1 = 'PROG' 0. 'PAS' 200. 3000. 'PAS' 1000. 10000. 'PAS' 10000. TFIN2 ;
  456. * ---------------------------------------*
  457. * Variables internes initiales 'PORO'=PO *
  458. * ---------------------------------------*
  459. VAR00 = 'ZERO' MODL1 'VARINTER' ;
  460. VAR01 = 'MANU' 'CHML' MODL1 'PORO' PO 'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  461. VAR0 = VAR00 + VAR01 ;
  462. *
  463. CHARTOT = CHARMECA 'ET' CHARTEMP 'ET' CHARFISS ;
  464. *
  465. TAB1 = TABLE ;
  466. * TAB1.'TEMPERATURES' = TABLE ;
  467. TAB1.'VARIABLES_INTERNES' = TABLE ;
  468. TAB1.'BLOCAGES_MECANIQUES' = CLT ;
  469. TAB1.'MODELE' = MODL1 ;
  470. TAB1.'CHARGEMENT' = CHARTOT ;
  471. TAB1.'VARIABLES_INTERNES' . 0 = VAR0 ;
  472. TAB1.'CARACTERISTIQUES' = MATR1 ;
  473. TAB1.'TEMPS_CALCULES' = LIST1 ;
  474. TAB1.'TEMPS_SAUVES' = LIST1 ;
  475. *TAB1.'TALPHA_REFERENCE' = TALPHA ;
  476. *
  477. TMASAU=table;
  478. tab1 . 'MES_SAUVEGARDES'=TMASAU;
  479. TMASAU .'DEFTO'=VRAI;
  480. TMASAU .'DEFIN'=VRAI;
  481.  
  482. PASAPAS TAB1 ;
  483. *
  484. *
  485. TT = TAB1 . 'TEMPS' ;
  486. SS = TAB1 . 'CONTRAINTES' ;
  487. VV = TAB1 . 'VARIABLES_INTERNES' ;
  488. IN = TAB1 . 'DEFORMATIONS_INELASTIQUES' ;
  489. NCONT = 'DIME' TAB1 . 'CONTRAINTES' ;
  490. *
  491. LZEIT = 'PROG' ;
  492. LSIGZ = 'PROG' ;
  493. LW1 = 'PROG' ;
  494. LW2 = 'PROG' ;
  495. LW3 = 'PROG' ;
  496. LEPSZF = 'PROG' ;
  497. *
  498. *
  499. 'REPE' BLOC NCONT ;
  500. *
  501. IND = &BLOC - 1 ;
  502. ZEIT = TT.IND ;
  503. *
  504. SIG = SS.IND ;
  505. SIGZ = 'EXCO' SIG 'SMZZ' ;
  506. SIGZ_MAX = 'MAXI' SIGZ ; SIGZ_MIN = 'MINI' SIGZ ;
  507. *
  508. VI = VV.IND ;
  509. EI = IN.IND ;
  510. *
  511. W1 = 'MAXI' ('EXCO' VI 'W1') ;
  512. W2 = 'MAXI' ('EXCO' VI 'W2') ;
  513. W3 = 'MAXI' ('EXCO' VI 'W3') ;
  514. EPSZF = 'MAXI' ('EXCO' EI 'EIZZ') ;
  515. *
  516. LZEIT = LZEIT 'ET' ('PROG' ZEIT) ;
  517. LSIGZ = LSIGZ 'ET' ('PROG' SIGZ_MAX) ;
  518. LW1 = LW1 'ET' ('PROG' W1) ;
  519. LW2 = LW2 'ET' ('PROG' W2) ;
  520. LW3 = LW3 'ET' ('PROG' W3) ;
  521. LEPSZF = LEPSZF 'ET' ('PROG' EPSZF) ;
  522. *
  523. 'FIN' BLOC ;
  524. *
  525. *
  526. * -----------------------*
  527. * Controle des resultats *
  528. * -----------------------*
  529. LSOUS = 'PROG' 1E3 'PAS' 2E3 9E3 1E5 8E5 1.5E6 1.6E6 ;
  530.  
  531. UN = 'PROG' ('DIME' LSOUS)*1 ;
  532. *
  533. SZ_CAL = 'IPOL' LSOUS LZEIT LSIGZ ;
  534. SZ_ANA = 'IPOL' LSOUS LTIME (('EXTR' EVSIGZ 'ORDO') * CONVSIG) ;
  535. EFZ_CAL = 'IPOL' LSOUS LZEIT LEPSZF ;
  536. EFZ_ANA = 'IPOL' LSOUS LTIME ('EXTR' EVEPSFZ 'ORDO') ;
  537. WZ_CAL = 'IPOL' LSOUS LZEIT LW1 ;
  538. WZ_ANA = 'IPOL' LSOUS LTIME ('EXTR' EVWZ 'ORDO') ;
  539. *
  540. PREX = 5E-2 ;
  541. T_SZ = (SZ_ANA - SZ_CAL) / (SZ_ANA + UN) ;
  542. T_SZ = 'MAXI' ('ABS' T_SZ) ;
  543. 'SI' ( T_SZ '&lt;EG' PREX ) ;
  544. 'ERRE' 0 ;
  545. 'SINO' ;
  546. 'MESS' 'ERREUR MAXIMALE CONTRAINTE SELON Z > ' PREX ;
  547. 'ERRE' 5 ;
  548. 'FINS' ;
  549. *
  550. WRUPT = (2.*XGFTR1) / XLTR1 ;
  551. ERUPT = WRUPT / LZ ;
  552. NONZERO = (ERUPT / 1E6) * UN ;
  553. T_FZ = (EFZ_ANA - EFZ_CAL) / (EFZ_ANA + NONZERO) ;
  554. T_FZ = 'MAXI' ('ABS' T_FZ) ;
  555. 'SI' ( T_FZ '&lt;EG' PREX ) ;
  556. 'ERRE' 0 ;
  557. 'SINO' ;
  558. 'MESS' 'ERREUR MAXIMALE DEFORMATION DE FISSURATION SELON Z > ' PREX ;
  559. 'ERRE' 5 ;
  560. 'FINS' ;
  561. *
  562. NONZERO = (WRUPT / 1E6) * UN ;
  563. T_WZ = (WZ_ANA - WZ_CAL) / (WZ_ANA + NONZERO) ;
  564. T_WZ = 'MAXI' ('ABS' T_WZ) ;
  565. 'SI' ( T_WZ '&lt;EG' PREX ) ;
  566. 'ERRE' 0 ;
  567. 'SINO' ;
  568. 'MESS' 'ERREUR MAXIMALE OUVERTURE SELON Z > ' PREX ;
  569. 'ERRE' 5 ;
  570. 'FINS' ;
  571. *
  572. *
  573. TAB = 'TABLE' ;
  574. TAB . 'TITRE' = 'TABLE' ;
  575. TAB . 'TITRE' . 1 = 'Solution Analytique' ;
  576. TAB . 'TITRE' . 2 = 'Solution Calculee' ;
  577. *
  578. *
  579. EVSIGZ0 = 'EVOL' 'MANU' 'Temps (s)' LZEIT 'SIGZ(MPa)' (LSIGZ/CONVSIG);
  580. 'SI' DESSI ;
  581. 'TITR' 'CONTRAINTE DANS LA DIRECTION Z' ;
  582. 'DESS' (EVSIGZ 'ET' EVSIGZ0) TAB 'LEGE' ;
  583. 'FINSI' ;
  584. *
  585. EVWZ0 = 'EVOL' 'MANU' 'Temps (s)' LZEIT 'WZ(m)' LW1 ;
  586. 'SI' DESSI ;
  587. 'TITR' 'OUVERTURE DANS LA DIRECTION Z' ;
  588. 'DESS' (EVWZ 'ET' EVWZ0) TAB 'LEGE' ;
  589. 'FINSI' ;
  590. *
  591. EVEPSFZ0 = 'EVOL' 'MANU' 'Temps (s)' LZEIT 'EPSFZ' LEPSZF ;
  592. 'SI' DESSI ;
  593. 'TITR' 'DEFORMATION DE FISSURATION DANS LA DIRECTION Z' ;
  594. 'DESS' (EVEPSFZ 'ET' EVEPSFZ0) TAB 'LEGE' ;
  595. 'FINSI' ;
  596. *
  597. EVFISZ = 'EVOL' 'MANU' 'WZ(m)' LW1 'SIGZ(MPa)' (LSIGZ/CONVSIG) ;
  598. 'SI' DESSI ;
  599. 'TITR' 'RELATION CONTRAINTE/OUVERTURE DANS LA DIRECTION Z' ;
  600. 'DESS' (EV_FISZ 'ET' EVFISZ) TAB 'LEGE' ;
  601. 'FINSI' ;
  602. *
  603. 'FIN';
  604.  
  605.  
  606.  
  607.  
  608.  
  609.  
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  

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