Télécharger rousselier.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rousselier.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *****************************************************
  5. * *
  6. * TEST DE VALIDATION DE LA LOI *
  7. * D'ENDOMMAGEMENT DUCTILE DE *
  8. * ROUSSELIER *
  9. * *
  10. * TEST: *
  11. * UN BARREAU EST CHARGE EN TRACTION *
  12. * LE CHARGEMENT EST DES DEPLACEMENTS IMPOSES *
  13. * CALCUL EN MASSIF 3D *
  14. * *
  15. * MAILLAGE: *
  16. * UN BARREAU DE SECTION CARREE DE .5 M SUR .01 M *
  17. * 50 ELEMENTS *
  18. * *
  19. *****************************************************
  20. opti echo 1 ;
  21. graph = faux ;
  22. *
  23. 'DEBPROC' PLROUSS ta1*'TABLE' beta0*'FLOTTANT'
  24. xm0*'FLOTTANT' sig0*'FLOTTANT'
  25. bet1*'FLOTTANT' def1*'FLOTTANT'
  26. sig1*'FLOTTANT' i1*'ENTIER' ;
  27. *
  28. * nombre d'itérations
  29. *
  30. niter = 1 ;
  31. *
  32. t0 = ta1 .'TEMPS'.i1 ;
  33. t1 = ta1 .'TEMPS'.( i1 - 1 ) ;
  34. dt0 = ( t0 - t1 ) / niter ;
  35. mat0 = ta1 .'CARACTERISTIQUES' ;
  36. mod0 = ta1. 'MODELE' ;
  37. depl0 = ta1 .'DEPLACEMENTS'.i1 ;
  38. depl1 = ta1 .'DEPLACEMENTS'.( i1 - 1 ) ;
  39. youn0 = 'MAXI' ( 'EXCO' mat0 'YOUN' ) ;
  40. xnu = 'MAXI' ( 'EXCO' mat0 'NU' ) ;
  41. F0 = 'MAXI' ( 'EXCO' mat0 'F' ) ;
  42. XD0 = 'MAXI' ( 'EXCO' mat0 'D' ) ;
  43. SIG = 'MAXI' ( 'EXCO' mat0 'SIG1' ) ;
  44. epst = 'EPSI' mod0 depl0 ;
  45. epstold = 'EPSI' mod0 depl1 ;
  46. def_x = 'MAXI' ( 'EXCO' epst 'EPXX' ) ;
  47. def_y = 'MAXI' ( 'EXCO' epst 'EPYY' ) ;
  48. def_z = 'MAXI' ( 'EXCO' epst 'EPZZ' ) ;
  49. def_xold = 'MAXI' ( 'EXCO' epstold 'EPXX' ) ;
  50. def_yold = 'MAXI' ( 'EXCO' epstold 'EPYY' ) ;
  51. def_zold = 'MAXI' ( 'EXCO' epstold 'EPZZ' ) ;
  52. def00 = ( def_x - def_xold ) / niter ;
  53. XK0 = youn0 / ( 3. * ( 1. - ( 2. * xnu ) ) ) ;
  54. *
  55. * trace de la déformation totale
  56. *
  57. treps1 = def_x + def_y + def_z ;
  58. trepsold = def_xold + def_yold + def_zold ;
  59. treps00 = ( treps1 - trepsold ) / niter ;
  60. *
  61. def_int = def1 ;
  62. bet_int = bet1 ;
  63. sig_int = sig1 ;
  64. *
  65. 'REPETER' bloc1 niter ;
  66. i10 = &bloc1 ;
  67. *
  68. treps0 = trepsold + ( i10 * treps00 ) ;
  69. def0 = def_xold + ( i10 * def00 ) ;
  70. *
  71. x0 = def_int ;
  72. x2 = x0 ;
  73. sigt = ( youn0 * ( def0 - x2 ) ) - sig_int ;
  74. cont0 = ( youn0 * ( def0 - x2 ) ) ;
  75. cont0 = cont0 / 3. ;
  76. dpeq1 = XD0 * ( 'EXP' ( cont0 / SIG ) ) ;
  77. *bet2 = bet_int +
  78. *( ( x2 - def_int ) * dpeq1 ) +
  79. *( ( ( x2 - def_int ) * ( x2 - def_int ) / 2. ) *
  80. *dpeq1 * sigt / ( 3. * SIG ) ) ;
  81. tr2 = ( 1. - ( 2. * xnu ) ) * ( def0 - x2 ) ;
  82. tr1 = treps0 - tr2 ;
  83. 'SI' ( tr1 '>' 1.E-5 ) ;
  84. bet2 = 'LOG' ( ( ( 'EXP' tr1 ) - 1. + F0 ) / F0 ) ;
  85. 'SINON' ;
  86. bet2 = tr1 / F0 ;
  87. 'FINSI' ;
  88. B1 = SIG * F0 * ( 'EXP' bet2 ) ;
  89. B1 = B1 / ( 1. - F0 + ( F0 * ( 'EXP' bet2 ) ) ) ;
  90. tr1 = tr1 / 3. ;
  91. y2 = sig0 * ( ( 1. + ( beta0 * ( x2 - tr1 ) ) ) ** xm0 ) ;
  92. y2 = y2 -
  93. ( B1 * XD0 * ( 'EXP' ( youn0 * ( def0 - x2 )
  94. / ( SIG * 3. ) ) ) ) ;
  95. y2 = y2 - ( youn0 * ( def0 - x2 ) ) ;
  96. y0 = y2 ;
  97. *
  98. * Cas élastique
  99. *
  100. 'SI' ( y0 '>' 0. ) ;
  101. def_th0 = def_int ;
  102. bet_th0 = bet_int ;
  103. sig_th0 = youn0 * ( def0 - def_th0 ) ;
  104. 'ITERER' bloc1 ;
  105. 'FINSI' ;
  106. *
  107. x1 = def0 ;
  108. x2 = x1 ;
  109. sigt = ( youn0 * ( def0 - x2 ) ) - sig_int ;
  110. cont0 = ( youn0 * ( def0 - x2 ) ) ;
  111. cont0 = cont0 / 3. ;
  112. dpeq1 = XD0 * ( 'EXP' ( cont0 / SIG ) ) ;
  113. tr2 = ( 1. - ( 2. * xnu ) ) * ( def0 - x2 ) ;
  114. tr1 = treps0 - tr2 ;
  115. 'SI' ( tr1 '>' 1.E-5 ) ;
  116. bet2 = 'LOG' ( ( ( 'EXP' tr1 ) - 1. + F0 ) / F0 ) ;
  117. 'SINON' ;
  118. bet2 = tr1 / F0 ;
  119. 'FINSI' ;
  120. B1 = SIG * F0 * ( 'EXP' bet2 ) ;
  121. B1 = B1 / ( 1. - F0 + ( F0 * ( 'EXP' bet2 ) ) ) ;
  122. tr1 = tr1 / 3. ;
  123. y2 = sig0 * ( ( 1. + ( beta0 * ( x2 - tr1 ) ) ) ** xm0 ) ;
  124. y2 = y2 -
  125. ( B1 * XD0 * ( 'EXP' ( youn0 * ( def0 - x2 )
  126. / ( SIG * 3. ) ) ) ) ;
  127. y2 = y2 - ( youn0 * ( def0 - x2 ) ) ;
  128. y1 = y2 ;
  129. err1 = 'ABS' ( y0 - y1 ) ;
  130. *
  131. 'REPETER' bloc0 ;
  132. x2 = ( x0 + x1 ) / 2. ;
  133. sigt = ( youn0 * ( def0 - x2 ) ) - sig_int ;
  134. cont0 = ( youn0 * ( def0 - x2 ) ) ;
  135. cont0 = cont0 / 3. ;
  136. dpeq1 = XD0 * ( 'EXP' ( cont0 / SIG ) ) ;
  137. tr2 = ( 1. - ( 2. * xnu ) ) * ( def0 - x2 ) ;
  138. tr1 = treps0 - tr2 ;
  139. 'SI' ( tr1 '>' 1.E-5 ) ;
  140. bet2 = 'LOG' ( ( ( 'EXP' tr1 ) - 1. + F0 ) / F0 ) ;
  141. 'SINON' ;
  142. bet2 = tr1 / F0 ;
  143. 'FINSI' ;
  144. B1 = SIG * F0 * ( 'EXP' bet2 ) ;
  145. B1 = B1 / ( 1. - F0 + ( F0 * ( 'EXP' bet2 ) ) ) ;
  146. tr1 = tr1 / 3. ;
  147. y2 = sig0 * ( ( 1. + ( beta0 * ( x2 - tr1 ) ) ) ** xm0 ) ;
  148. y2 = y2 -
  149. ( B1 * XD0 * ( 'EXP' ( youn0 * ( def0 - x2 )
  150. / ( SIG * 3. ) ) ) ) ;
  151. y2 = y2 - ( youn0 * ( def0 - x2 ) ) ;
  152. *
  153. 'SI' ( y2 '>' ( 1.E-4 ) ) ;
  154. x1 = x2 ;
  155. 'SINON' ;
  156. 'SI' ( y2 '<' ( -1.E-4 ) ) ;
  157. x0 = x2 ;
  158. 'SINON' ;
  159. def_th0 = x2 ;
  160. bet_th0 = bet2 ;
  161. sig_th0 = youn0 * ( def0 - def_th0 ) ;
  162. 'QUITTER' bloc0 ;
  163. 'FINSI' ;
  164. 'FINSI' ;
  165. 'FIN' bloc0 ;
  166. *
  167. * Calcul du nouvel endommagement
  168. *
  169. *
  170. def_int = def_th0 ;
  171. bet_int = bet_th0 ;
  172. sig_int = sig_th0 ;
  173. *
  174. *'MESS' ' Fin iteration interne no' i10 ;
  175. *
  176. 'FIN' bloc1 ;
  177. *
  178. *bet_th0 = bet2 ;
  179. *
  180. rho_th0 = 1. / ( 1. - F0 + ( F0 * ( 'EXP' bet_th0 ) ) ) ;
  181. *
  182. 'FINPROC' sig_th0 def_th0 bet_th0 rho_th0 ;
  183. *
  184. *
  185. *================================================================
  186. opti echo 0 dime 3 elem cub8 ;
  187. * ajout de option epsilon lineaire pour la precision des test!
  188. OPTION epsilon lineaire;
  189. *
  190. *
  191. * Maillage
  192. *
  193. l1 = .5 ;
  194. l2 = .01 ;
  195. l3 = .01 ;
  196. n1 = 50 ;
  197. n2 = 1 ;
  198. n3 = 1 ;
  199. p1 = 0. 0. 0. ;
  200. p2 = l1 0. 0. ;
  201. p3 = l1 l2 0. ;
  202. p4 = 0. l2 0. ;
  203. p5 = 0. 0. l3 ;
  204. p6 = l1 0. l3 ;
  205. p7 = l1 l2 l3 ;
  206. p8 = 0. l2 l3 ;
  207. d1 = d p1 n3 p5 ;
  208. d2 = d p5 n2 p8 ;
  209. d3 = d p8 n3 p4 ;
  210. d4 = d p4 n2 p1 ;
  211. sur1 = 'DALL' d1 d2 d3 d4 ;
  212. d5 = d p2 n3 p6 ;
  213. d6 = d p6 n2 p7 ;
  214. d7 = d p7 n3 p3 ;
  215. d8 = d p3 n2 p2 ;
  216. sur2 = 'DALL' d5 d6 d7 d8 ;
  217. vol1 = sur1 'VOLU' n1 sur2 ;
  218. titr ' Maillage du test de validation' ;
  219. *trac vol1 ;
  220. *
  221. * Deplacement maximal
  222. *
  223. depmax = .01 ;
  224. *
  225. * Courbe d'ecrouissage
  226. *
  227. beta0 = 125. ;
  228. sig0 = 410.E6 ;
  229. xm = .1 ;
  230. youn0 = 2.1E11 ;
  231. epsp0 = 'PROG' 0. ;
  232. deb0 = 'PROG' 0. ;
  233. epsmax0 = depmax / l1 ;
  234. pas0 = epsmax0 / 50 ;
  235. pas1 = pas0 / 10. ;
  236. epsmax1 = ( epsmax0 / 10. ) ;
  237. epsp1 = ( 'PROG' 0. 'PAS' pas1 epsmax1 ) 'ET'
  238. ( 'PROG' ( epsmax1 + pas0 ) 'PAS' pas0 epsmax0 ) ;
  239. *epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  240. dim00 = 'DIME' epsp1 ;
  241. un = 'PROG' dim00 '*' 1. ;
  242. Y0 = sig0 ;
  243. epsp2 = beta0 * epsp1 ;
  244. epsp3 = epsp2 + un ;
  245. epsp31 = log epsp3 ;
  246. epsp32 = epsp31 * xm ;
  247. epsp4 = exp epsp32 ;
  248. Y1 = Y0 * epsp4 ;
  249. *
  250. epsp1 = epsp0 'ET' epsp1 ;
  251. Y1 = deb0 'ET' Y1 ;
  252. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  253. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ;
  254. *
  255. * Modele de calcul
  256. *
  257. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ;
  258. lsm1 = y1 enle 1 ;
  259. lep1 = (epsp2 enle 1) - (lsm1 / youn0) ;
  260. evec = evol vert manu eps lep1 sig lsm1 ;
  261. TAB = 'TABLE' ;
  262. TAB. 1 = 'MARQ CROI ' ;
  263. TAB.'TITRE' = 'TABLE' ;
  264. TAB.'TITRE'. 1 = MOT ' TRACTION' ;
  265. TAB.'TITRE'. 2 = MOT 'ECROUISSAGE' ;
  266. si graph ;
  267. 'DESS' (trac0 et evec)
  268. 'TITR' ' Courbes de traction et d ecrouissage (vert)' 'LEGE' TAB ;
  269. fins ;
  270. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  271. 'PLASTIQUE_ENDOM' 'ROUSSELIER' ;
  272. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' 0.3
  273. 'ECRO' evec 'F' 3.E-3 'FC' 1. 'D' 2. 'SIG1' 490.E6 ;
  274. *
  275. * Conditions aux limites
  276. *
  277. cl1 = 'BLOQ' sur1 'UX' ;
  278. cl2 = 'BLOQ' sur2 'UX' ;
  279. cl3 = 'BLOQ' p1 'UX' 'UY' 'UZ' ;
  280. cl4 = 'BLOQ' p5 'UX' 'UY' ;
  281. cl5 = 'BLOQ' p4 'UX' 'UZ' ;
  282. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 'ET' cl5 ;
  283. *
  284. * Chargement
  285. *
  286. dep1 = 'DEPI' cl2 depmax ;
  287. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ;
  288. cha0 = 'CHAR' 'DIMP' dep1 ev0 ;
  289. *
  290. * Temps du calcul
  291. *
  292. dt0 = 50. ;
  293. tfin0 = 500. ;
  294. *
  295. *
  296. * Resolution par PASAPAS
  297. *---------------------------------------
  298. *
  299. ta1 = 'TABLE' ;
  300. ta1 .'MODELE' = mod0 ;
  301. ta1 .'CARACTERISTIQUES' = mat0 ;
  302. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  303. ta1 .'CHARGEMENT' = cha0 ;
  304. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' 1. 20. ) 'ET'
  305. ( 'PROG' 25. 'PAS' 5. 45. ) 'ET'
  306. ( 'PROG' 50. 'PAS' dt0 tfin0) ;
  307. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' dt0 tfin0) ;
  308. *
  309. TMASAU=table;
  310. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  311. TMASAU .'DEFTO'=VRAI;
  312. TMASAU .'DEFIN'=VRAI;
  313. PASAPAS ta1 ;
  314. *
  315. *
  316. * Post-traitement
  317. *-----------------------------------------
  318. *
  319. * Quelques donnees du calcul
  320. *
  321. mat0 = ta1 .'CARACTERISTIQUES' ;
  322. youn1 = 'EXTR' ( 'EXCO' mat0 'YOUN' 'SCAL' ) 'SCAL' 1 1 1 ;
  323. xnu0 = 'EXTR' ( 'EXCO' mat0 'NU' 'SCAL' ) 'SCAL' 1 1 1 ;
  324. G0 = youn1 / ( 2. * ( 1. + xnu0 ) ) ;
  325. XD0 = 'EXTR' ( 'EXCO' mat0 'D' 'SCAL' ) 'SCAL' 1 1 1 ;
  326. XF0 = 'EXTR' ( 'EXCO' mat0 'F' 'SCAL' ) 'SCAL' 1 1 1 ;
  327. SIG = 'EXTR' ( 'EXCO' mat0 'SIG1' 'SCAL' ) 'SCAL' 1 1 1 ;
  328. *
  329. * Post-traitement
  330. *-----------------------------------------
  331. *
  332. sigx1 = 'PROG' 0. ;
  333. sigx2 = 'PROG' 0. ;
  334. depx1 = 'PROG' 0. ;
  335. depx2 = 'PROG' 0. ;
  336. sig_1 = 'PROG' 0. ;
  337. def_1 = 'PROG' 0. ;
  338. bet_1 = 'PROG' 0. ;
  339. rho_1 = 'PROG' 1. ;
  340. def_th1 = 'PROG' 0. ;
  341. sig_th1 = 'PROG' 0. ;
  342. bet_th1 = 'PROG' 0. ;
  343. rho_th1 = 'PROG' 1. ;
  344. err_e1 = 'PROG' 0. ;
  345. err_s1 = 'PROG' 0. ;
  346. err_b1 = 'PROG' 0. ;
  347. err_r1 = 'PROG' 0. ;
  348. dim0 = 'DIME' ta1 .'TEMPS' - 1 ;
  349. mod0 = ta1 .'MODELE' ;
  350. sig1 = 0. ;
  351. def1 = 0. ;
  352. bet1 = 0. ;
  353. rho1 = 1. ;
  354. eps_el0 = sig0 / youn1 ;
  355. mess ' eps_el0 = ' eps_el0 ;
  356. 'REPETER' bloc1 dim0 ;
  357. i1 = &bloc1 ;
  358. temp0 = ta1 .'TEMPS'.i1 ;
  359. mat0 = ta1 .'CARACTERISTIQUES' ;
  360. depl0 = ta1 .'DEPLACEMENTS'.i1 ;
  361. sigm0 = ta1 .'CONTRAINTES'.i1 ;
  362. depeq0 = ta1 .'VARIABLES_INTERNES'.i1 ;
  363. def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i1 ;
  364. sigx0 = 'MAXI' ( 'VMIS' mod0 sigm0 ) ;
  365. sigx00 = 'MINI' ( 'VMIS' mod0 sigm0 ) ;
  366. depx00 = 'MINI' ( 'EXCO' depeq0 'EPSE' 'SCAL' ) ;
  367. epsp0 = 'MAXI' ( 'EXCO' depeq0 'EPSE' 'SCAL' ) ;
  368. bet_0 = 'MAXI' ( 'EXCO' depeq0 'VHTA' 'SCAL' ) ;
  369. rho_0 = 'MAXI' ( 'EXCO' depeq0 'VHOT' 'SCAL' ) ;
  370. sig_00 = 'EXCO' sigm0 'SMXX' 'SCAL' ;
  371. sig_0 = 'EXTR' sig_00 'SCAL' 1 1 1 ;
  372. def_00 = 'EXCO' def0 'EIXX' 'SCAL' ;
  373. def_0 = 'EXTR' def_00 'SCAL' 1 1 1 ;
  374. sigx1 = sigx1 'ET' ( 'PROG' sigx0 ) ;
  375. sigx2 = sigx2 'ET' ( 'PROG' sigx00 ) ;
  376. depx1 = depx1 'ET' ( 'PROG' epsp0 ) ;
  377. depx2 = depx2 'ET' ( 'PROG' depx00 ) ;
  378. sig_1 = sig_1 'ET' ( 'PROG' sig_0 ) ;
  379. def_1 = def_1 'ET' ( 'PROG' def_0 ) ;
  380. bet_1 = bet_1 'ET' ( 'PROG' bet_0 ) ;
  381. rho_1 = rho_1 'ET' ( 'PROG' rho_0 ) ;
  382. *
  383. * deformations totales
  384. *
  385. eps0 = 'EPSI' mod0 depl0 ;
  386. eps_x = 'MAXI' ( 'EXCO' eps0 'EPXX' ) ;
  387. *
  388. * Calcul de la solution théorique
  389. *
  390. *B1 = SIG * XF0 * ( 'EXP' bet1 ) ;
  391. *B1 = B1 / ( 1. - XF0 + ( XF0 * ( 'EXP' bet1 ) ) ) ;
  392. B1 = SIG * XF0 ;
  393. *F2 = B1 * XD0 * ( 'EXP' ( sig1 / ( 3. * SIG ) ) ) ;
  394. F2 = B1 * XD0 * ( 'EXP' ( sig0 / ( 3. * SIG ) ) ) ;
  395. F1 = sig0 ;
  396. YLIM = F1 - F2 ;
  397. sigt = youn1 * ( eps_x - def1 ) ;
  398. 'SI' ( sigt '>' YLIM ) ;
  399. sig1 = sig1 / rho1 ;
  400. sig_th0 def_th0 bet_th0 rho_th0 = PLROUSS ta1 beta0
  401. xm sig0 bet1 def1 sig1 i1 ;
  402. sig_th0 = sig_th0 * rho_th0 ;
  403. 'SINON' ;
  404. sig_th0 = sigt ;
  405. def_th0 = def1 ;
  406. bet_th0 = bet1 ;
  407. rho_th0 = rho1 ;
  408. 'FINSI' ;
  409. *
  410. sig1 = sig_th0 ;
  411. def1 = def_th0 ;
  412. bet1 = bet_th0 ;
  413. rho1 = rho_th0 ;
  414. *
  415. def_th1 = def_th1 'ET' ( 'PROG' def_th0 ) ;
  416. sig_th1 = sig_th1 'ET' ( 'PROG' sig_th0 ) ;
  417. bet_th1 = bet_th1 'ET' ( 'PROG' bet_th0 ) ;
  418. rho_th1 = rho_th1 'ET' ( 'PROG' rho_th0 ) ;
  419. *
  420. * 'MESS' ' deformation plastique theorique' def_th0 ;
  421. * 'MESS' ' deformation plastique calculee ' def_0 ;
  422. * 'MESS' ' contrainte théorique ' sig_th0 ;
  423. * 'MESS' ' contrainte calculee ' sig_0 ;
  424. * 'MESS' ' endommagement théorique ' bet_th0 ;
  425. * 'MESS' ' endommagement calcule ' bet_0 ;
  426. * 'MESS' ' fraction de densité théorique ' rho_th0 ;
  427. * 'MESS' ' fraction de densité calculee ' rho_0 ;
  428. *
  429. * Calcul d'erreur
  430. *
  431. 'SI' ( def_0 '>' ( 1.E-1 * eps_el0 ) ) ;
  432. err_e0 = ( def_th0 - def_0 ) / def_0 ;
  433. err_e0 = 'ABS' err_e0 ;
  434. err_e0 = err_e0 * 100. ;
  435. 'SINON' ;
  436. err_e0 = 0. ;
  437. 'FINSI' ;
  438. 'SI' ( sig_0 '>' 0. ) ;
  439. err_s0 = ( sig_th0 - sig_0 ) / sig_0 ;
  440. err_s0 = 'ABS' err_s0 ;
  441. err_s0 = err_s0 * 100. ;
  442. 'SINON' ;
  443. err_s0 = 0. ;
  444. 'FINSI' ;
  445. 'SI' ( sig_th0 '>' ( ( ( 9. * sig0 ) + SIG1 ) / 10. ) ) ;
  446. err_b0 = ( bet_th0 - bet_0 ) / bet_0 ;
  447. err_b0 = 'ABS' err_b0 ;
  448. err_b0 = err_b0 * 100. ;
  449. 'SINON' ;
  450. err_b0 = 0. ;
  451. 'FINSI' ;
  452. 'SI' ( err_e0 '>' 0. ) ;
  453. err_r0 = ( rho_th0 - rho_0 ) / rho_0 ;
  454. err_r0 = 'ABS' err_r0 ;
  455. err_r0 = err_r0 * 100. ;
  456. 'SINON' ;
  457. err_r0 = 0. ;
  458. 'FINSI' ;
  459. *
  460. err_e1 = err_e1 'ET' ( 'PROG' err_e0 ) ;
  461. err_s1 = err_s1 'ET' ( 'PROG' err_s0 ) ;
  462. err_b1 = err_b1 'ET' ( 'PROG' err_b0 ) ;
  463. err_r1 = err_r1 'ET' ( 'PROG' err_r0 ) ;
  464. *
  465. 'SI' ( err_e0 '>' 2.5E-2 ) ;
  466. 'ERRE' 5 ;
  467. 'SINON' ;
  468. 'ERRE' 0 ;
  469. 'FINSI' ;
  470. 'SI' ( err_s0 '>' 2.5E-2 ) ;
  471. 'ERRE' 5 ;
  472. 'SINON' ;
  473. 'ERRE' 0 ;
  474. 'FINSI' ;
  475. 'SI' ( err_b0 '>' 6.5 ) ;
  476. 'ERRE' 5 ;
  477. 'SINON' ;
  478. 'ERRE' 0 ;
  479. 'FINSI' ;
  480. 'SI' ( err_r0 '>' 2.5E-2 ) ;
  481. 'ERRE' 5 ;
  482. 'SINON' ;
  483. 'ERRE' 0 ;
  484. 'FINSI' ;
  485. *
  486. 'MESS' ' Fin du calcul de restitution, temps ' temp0 ;
  487. * 'MESS' ' ' ;
  488. *
  489. 'FIN' bloc1 ;
  490. *
  491. * Trace
  492. *
  493. trac1 = 'EVOL' 'ROUG' 'MANU' 'epsp' depx1 'sigm' sigx1 ;
  494. TAB = 'TABLE' ;
  495. TAB.1 = 'MARQ CROI ' ;
  496. TAB.2 = 'MARQ PLUS ' ;
  497. TAB.'TITRE' = 'TABLE' ;
  498. TAB.'TITRE'. 1 = MOT 'CALC' ;
  499. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  500. *'DESSIN' ( trac1 'ET' trac00 )
  501. *'TITR' 'Contraintes - Deformation plastique'
  502. *'LEGE' TAB ;
  503. trac2 = 'EVOL' 'ROUG' 'MANU' 'epsp' depx2 'sigm' sigx2 ;
  504. TAB = 'TABLE' ;
  505. TAB.1 = 'MARQ CROI ' ;
  506. TAB.2 = 'MARQ PLUS ' ;
  507. TAB.'TITRE' = 'TABLE' ;
  508. TAB.'TITRE'. 1 = MOT 'MAXI' ;
  509. TAB.'TITRE'. 2 = MOT 'MINI' ;
  510. *'DESSIN' ( trac1 'ET' trac2 )
  511. *'TITR' 'Comparaison mini - maxi de la contrainte equivalente '
  512. *'LEGE' TAB ;
  513. trac3 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'sigm' sig_th1 ;
  514. trac4 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_1 'sigm' sig_1 ;
  515. TAB = 'TABLE' ;
  516. TAB.1 = 'MARQ CROI ' ;
  517. TAB.2 = 'MARQ PLUS ' ;
  518. TAB.'TITRE' = 'TABLE' ;
  519. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  520. TAB.'TITRE'. 2 = MOT 'CALCUL' ;
  521. *'DESSIN' ( trac3 'ET' trac4 )
  522. *'TITR' 'Comparaison calcul theorie '
  523. *'LEGE' TAB ;
  524. trac5 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_s1 ;
  525. trac6 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_e1 ;
  526. TAB = 'TABLE' ;
  527. TAB.1 = 'MARQ CROI ' ;
  528. TAB.2 = 'MARQ PLUS ' ;
  529. TAB.'TITRE' = 'TABLE' ;
  530. TAB.'TITRE'. 1 = MOT 'CONTRAINTE' ;
  531. TAB.'TITRE'. 2 = MOT 'DEFORMATION' ;
  532. *'DESSIN' ( trac5 'ET' trac6 )
  533. *'TITR' ' Evolution des erreurs '
  534. *'LEGE' TAB ;
  535. trac7 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_b1 ;
  536. trac8 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_r1 ;
  537. TAB = 'TABLE' ;
  538. TAB.1 = 'MARQ CROI ' ;
  539. TAB.2 = 'MARQ PLUS ' ;
  540. TAB.'TITRE' = 'TABLE' ;
  541. TAB.'TITRE'. 1 = MOT 'ENDOMMAGAMENT' ;
  542. TAB.'TITRE'. 2 = MOT 'DENSITE' ;
  543. *'DESSIN' ( trac7 'ET' trac8 )
  544. *'TITR' ' Evolution des erreurs '
  545. *'LEGE' TAB ;
  546. TAB = 'TABLE' ;
  547. TAB.1 = 'MARQ CROI ' ;
  548. TAB.2 = 'MARQ PLUS ' ;
  549. TAB.'TITRE' = 'TABLE' ;
  550. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  551. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  552. *'DESSIN' ( trac3 'ET' trac00 )
  553. *'TITR' 'Comparaison theorie traction '
  554. *'LEGE' TAB ;
  555. trac10 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'beta'bet_1 ;
  556. trac11 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'beta'bet_th1 ;
  557. TAB = 'TABLE' ;
  558. TAB.1 = 'MARQ CROI ' ;
  559. TAB.2 = 'MARQ PLUS ' ;
  560. TAB.'TITRE' = 'TABLE' ;
  561. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  562. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  563. *'DESSIN' ( trac11 'ET' trac10 )
  564. *'TITR' 'Comparaison endommagement theorique et calcule '
  565. *'LEGE' TAB ;
  566. TAB = 'TABLE' ;
  567. TAB.1 = 'MARQ CROI ' ;
  568. TAB.2 = ' ' ;
  569. TAB.'TITRE' = 'TABLE' ;
  570. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  571. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  572. *'DESSIN' ( trac4 'ET' trac00 )
  573. *'TITR' 'Comparaison calcul traction '
  574. *'LEGE' TAB ;
  575. *
  576. *
  577. 'FIN' ;
  578.  
  579.  
  580.  
  581.  
  582.  
  583.  
  584.  
  585.  
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  

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