Télécharger uo2s_cas2.dgibi

Retour à la liste

Numérotation des lignes :

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

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