Télécharger la plaquette Cast3M


fr en

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 'LINE' ; 
  45. epstold = 'EPSI' mod0 depl1 'LINE' ; 
  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. *  
  188. *  
  189. *  Maillage 
  190. *  
  191. l1 = .5 ; 
  192. l2 = .01 ; 
  193. l3 = .01 ;  
  194. n1 = 50 ; 
  195. n2 = 1 ; 
  196. n3 = 1 ; 
  197. p1 = 0. 0. 0. ; 
  198. p2 = l1 0. 0. ; 
  199. p3 = l1 l2 0. ; 
  200. p4 = 0. l2 0. ; 
  201. p5 = 0. 0. l3 ; 
  202. p6 = l1 0. l3 ; 
  203. p7 = l1 l2 l3 ; 
  204. p8 = 0. l2 l3 ; 
  205. d1 = d p1 n3 p5 ; 
  206. d2 = d p5 n2 p8 ; 
  207. d3 = d p8 n3 p4 ; 
  208. d4 = d p4 n2 p1 ; 
  209. sur1 = 'DALL' d1 d2 d3 d4 ; 
  210. d5 = d p2 n3 p6 ; 
  211. d6 = d p6 n2 p7 ;  
  212. d7 = d p7 n3 p3 ; 
  213. d8 = d p3 n2 p2 ; 
  214. sur2 = 'DALL' d5 d6 d7 d8 ; 
  215. vol1 = sur1 'VOLU' n1 sur2 ;  
  216. titr ' Maillage du test de validation' ; 
  217. *trac vol1 ; 
  218. *  
  219. *  Deplacement maximal 
  220. * 
  221. depmax = .01 ;
  222. *   
  223. *  Courbe d'ecrouissage 
  224. *
  225. beta0 = 125. ; 
  226. sig0 = 410.E6 ; 
  227. xm = .1 ;       
  228. youn0 = 2.1E11 ; 
  229. epsp0 = 'PROG' 0. ; 
  230. deb0 = 'PROG' 0. ;  
  231. epsmax0 = depmax / l1 ;  
  232. pas0 = epsmax0 / 50 ; 
  233. pas1 = pas0 / 10. ; 
  234. epsmax1 = ( epsmax0 / 10. ) ; 
  235. epsp1 =  ( 'PROG' 0. 'PAS' pas1 epsmax1 ) 'ET'
  236. ( 'PROG' ( epsmax1 + pas0 )  'PAS' pas0 epsmax0 ) ; 
  237. *epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;  
  238. dim00 = 'DIME' epsp1 ;  
  239. un = 'PROG' dim00 '*' 1. ; 
  240. Y0 = sig0 ;  
  241. epsp2 = beta0 * epsp1 ; 
  242. epsp3 = epsp2 + un ; 
  243. epsp31 = log epsp3 ; 
  244. epsp32 = epsp31 * xm ;   
  245. epsp4 = exp epsp32 ; 
  246. Y1 = Y0 * epsp4 ; 
  247. *  
  248. epsp1 = epsp0 'ET' epsp1 ;  
  249. Y1 = deb0 'ET' Y1 ;   
  250. epsp2 = epsp1 + ( Y1 / youn0 ) ;   
  251. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ; 
  252. *  
  253. *    Modele de calcul 
  254. * 
  255. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ; 
  256. lsm1  = y1 enle 1 ;
  257. lep1  = (epsp2 enle 1) - (lsm1 / youn0) ;
  258. evec  = evol vert manu eps lep1 sig lsm1 ;
  259. TAB = 'TABLE' ;
  260. TAB. 1 = 'MARQ CROI ' ;
  261. TAB.'TITRE' = 'TABLE' ;
  262. TAB.'TITRE'. 1 = MOT '   TRACTION' ;
  263. TAB.'TITRE'. 2 = MOT 'ECROUISSAGE' ;
  264. si graph ;
  265.   'DESS' (trac0 et evec)
  266.   'TITR' ' Courbes de traction et d ecrouissage (vert)' 'LEGE' TAB ;
  267. fins ;
  268. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
  269.        'PLASTIQUE_ENDOM' 'ROUSSELIER' ;  
  270. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' 0.3      
  271.        'ECRO' evec 'F' 3.E-3 'FC' 1. 'D' 2. 'SIG1' 490.E6 ; 
  272. *  
  273. *    Conditions aux limites 
  274. *  
  275. cl1 = 'BLOQ' sur1 'UX' ; 
  276. cl2 = 'BLOQ' sur2 'UX' ; 
  277. cl3 = 'BLOQ' p1 'UX' 'UY' 'UZ' ;  
  278. cl4 = 'BLOQ' p5 'UX' 'UY' ; 
  279. cl5 = 'BLOQ' p4 'UX' 'UZ' ; 
  280. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 'ET' cl5 ;  
  281. *  
  282. *    Chargement 
  283. *  
  284. dep1 = 'DEPI' cl2 depmax ; 
  285. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ; 
  286. cha0 = 'CHAR' 'DIMP' dep1 ev0 ; 
  287. *  
  288. *  Temps du calcul  
  289. *  
  290. dt0 = 50. ; 
  291. tfin0 = 500. ;  
  292. *   
  293. *   
  294. *  Resolution par PASAPAS  
  295. *---------------------------------------
  296. *    
  297. ta1 = 'TABLE' ; 
  298. ta1 .'MODELE' = mod0 ; 
  299. ta1 .'CARACTERISTIQUES' = mat0 ; 
  300. ta1 .'BLOCAGES_MECANIQUES' = cl0 ; 
  301. ta1 .'CHARGEMENT' = cha0 ; 
  302. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' 1. 20. ) 'ET'  
  303. ( 'PROG' 25. 'PAS' 5. 45. ) 'ET' 
  304. ( 'PROG' 50. 'PAS' dt0 tfin0) ;  
  305. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' dt0 tfin0) ; 
  306. *
  307. TMASAU=table;
  308. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  309. TMASAU .'DEFTO'=VRAI;
  310. TMASAU .'DEFIN'=VRAI;  
  311. ta1 .'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  312. PASAPAS ta1 ; 
  313. *  
  314. *  
  315. *  Post-traitement 
  316. *-----------------------------------------
  317. * 
  318. *  Quelques donnees du calcul 
  319. *  
  320. mat0 = ta1 .'CARACTERISTIQUES' ; 
  321. youn1 = 'EXTR' ( 'EXCO' mat0 'YOUN' 'SCAL' ) 'SCAL' 1 1 1 ; 
  322. xnu0 = 'EXTR' ( 'EXCO' mat0 'NU' 'SCAL' ) 'SCAL' 1 1 1 ;
  323. G0 = youn1 / ( 2. * ( 1. + xnu0 ) ) ; 
  324. XD0 = 'EXTR' ( 'EXCO' mat0 'D' 'SCAL' ) 'SCAL' 1 1 1 ;  
  325. XF0 = 'EXTR' ( 'EXCO' mat0 'F' 'SCAL' ) 'SCAL' 1 1 1 ; 
  326. SIG = 'EXTR' ( 'EXCO' mat0 'SIG1' 'SCAL' ) 'SCAL' 1 1 1 ; 
  327. *  
  328. *  Post-traitement 
  329. *-----------------------------------------
  330. * 
  331. sigx1 = 'PROG' 0. ; 
  332. sigx2 = 'PROG' 0. ; 
  333. depx1 = 'PROG' 0. ; 
  334. depx2 = 'PROG' 0. ;  
  335. sig_1 = 'PROG' 0. ; 
  336. def_1 = 'PROG' 0. ; 
  337. bet_1 = 'PROG' 0. ;
  338. rho_1 = 'PROG' 1. ;
  339. def_th1 = 'PROG' 0. ; 
  340. sig_th1 = 'PROG' 0. ; 
  341. bet_th1 = 'PROG' 0. ; 
  342. rho_th1 = 'PROG' 1. ; 
  343. err_e1 = 'PROG' 0. ;  
  344. err_s1 = 'PROG' 0. ; 
  345. err_b1 = 'PROG' 0. ; 
  346. err_r1 = 'PROG' 0. ; 
  347. dim0 = 'DIME' ta1 .'TEMPS' - 1 ; 
  348. mod0 = ta1 .'MODELE' ; 
  349. sig1 = 0. ; 
  350. def1 = 0. ; 
  351. bet1 = 0. ; 
  352. rho1 = 1. ; 
  353. eps_el0 = sig0 / youn1 ; 
  354. mess ' eps_el0 = ' eps_el0 ;  
  355. 'REPETER' bloc1 dim0 ; 
  356.   i1 = &bloc1 ; 
  357.   temp0 = ta1 .'TEMPS'.i1 ; 
  358.   mat0 = ta1 .'CARACTERISTIQUES' ;
  359.   depl0 = ta1 .'DEPLACEMENTS'.i1 ; 
  360.   sigm0 = ta1 .'CONTRAINTES'.i1 ; 
  361.   depeq0 = ta1 .'VARIABLES_INTERNES'.i1 ;  
  362.   def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i1 ;  
  363.   sigx0 = 'MAXI' ( 'VMIS' mod0 sigm0 ) ;     
  364.   sigx00 = 'MINI' ( 'VMIS' mod0 sigm0 ) ; 
  365.   depx00 = 'MINI' ( 'EXCO' depeq0 'EPSE' 'SCAL' ) ;  
  366.   epsp0 = 'MAXI' ( 'EXCO' depeq0 'EPSE' 'SCAL' ) ; 
  367.   bet_0 = 'MAXI' ( 'EXCO' depeq0 'VHTA' 'SCAL' ) ; 
  368.   rho_0 = 'MAXI' ( 'EXCO' depeq0 'VHOT' 'SCAL' ) ;  
  369.   sig_00 = 'EXCO' sigm0 'SMXX' 'SCAL' ; 
  370.   sig_0 = 'EXTR' sig_00 'SCAL' 1 1 1 ;    
  371.   def_00 = 'EXCO' def0 'EIXX' 'SCAL' ; 
  372.   def_0 = 'EXTR' def_00 'SCAL' 1 1 1 ; 
  373.   sigx1 = sigx1 'ET' ( 'PROG' sigx0 ) ; 
  374.   sigx2 = sigx2 'ET' ( 'PROG' sigx00 ) ;  
  375.   depx1 = depx1 'ET' ( 'PROG' epsp0 ) ;  
  376.   depx2 = depx2 'ET' ( 'PROG' depx00 ) ; 
  377.   sig_1 = sig_1 'ET' ( 'PROG' sig_0 ) ; 
  378.   def_1 = def_1 'ET' ( 'PROG' def_0 ) ;  
  379.   bet_1 = bet_1 'ET' ( 'PROG' bet_0 ) ; 
  380.   rho_1 = rho_1 'ET' ( 'PROG' rho_0 ) ; 
  381. *  
  382. *  deformations totales 
  383. *  
  384.   eps0 = 'EPSI' mod0 depl0 'LINE' ; 
  385.   eps_x = 'MAXI' ( 'EXCO' eps0 'EPXX' ) ; 
  386. *  
  387. *  Calcul de la solution théorique  
  388. *  
  389. *B1 = SIG * XF0 * ( 'EXP' bet1 ) ;
  390. *B1 = B1 / ( 1. - XF0 + ( XF0 * ( 'EXP' bet1 ) ) ) ; 
  391. B1 = SIG * XF0 ; 
  392. *F2 = B1 * XD0 * ( 'EXP' ( sig1 / ( 3. * SIG ) ) ) ; 
  393. F2 = B1 * XD0 * ( 'EXP' ( sig0 / ( 3. * SIG ) ) ) ; 
  394. F1 = sig0 ; 
  395. YLIM = F1 - F2 ; 
  396. sigt = youn1 * ( eps_x - def1 ) ; 
  397. 'SI' ( sigt '>' YLIM ) ; 
  398.     sig1 = sig1 / rho1 ; 
  399.     sig_th0 def_th0 bet_th0 rho_th0 = PLROUSS ta1 beta0 
  400.     xm sig0 bet1 def1 sig1 i1 ;   
  401.     sig_th0 = sig_th0 * rho_th0 ; 
  402. 'SINON' ; 
  403.     sig_th0 = sigt ; 
  404.     def_th0 = def1 ; 
  405.     bet_th0 = bet1 ; 
  406.     rho_th0 = rho1 ; 
  407. 'FINSI' ; 
  408. *  
  409. sig1 = sig_th0 ; 
  410. def1 = def_th0 ; 
  411. bet1 = bet_th0 ; 
  412. rho1 = rho_th0 ; 
  413. *  
  414. def_th1 = def_th1 'ET' ( 'PROG' def_th0 ) ; 
  415. sig_th1 = sig_th1 'ET' ( 'PROG' sig_th0 ) ; 
  416. bet_th1 = bet_th1 'ET' ( 'PROG' bet_th0 ) ; 
  417. rho_th1 = rho_th1 'ET' ( 'PROG' rho_th0 ) ; 
  418. *  
  419. *  'MESS' ' deformation plastique theorique' def_th0 ; 
  420. *  'MESS' ' deformation plastique calculee ' def_0 ; 
  421. *  'MESS' ' contrainte théorique           ' sig_th0 ; 
  422. *  'MESS' ' contrainte calculee            ' sig_0 ;
  423. *  'MESS' ' endommagement théorique        ' bet_th0 ; 
  424. *  'MESS' ' endommagement calcule          ' bet_0 ; 
  425. *  'MESS' ' fraction de densité théorique  ' rho_th0 ;  
  426. *  'MESS' ' fraction de densité calculee   ' rho_0 ; 
  427. *  
  428. *  Calcul d'erreur  
  429. *  
  430. 'SI' ( def_0 '>' ( 1.E-1 * eps_el0 ) ) ; 
  431.      err_e0 = ( def_th0 - def_0 ) / def_0 ; 
  432.      err_e0 = 'ABS' err_e0 ; 
  433.      err_e0 = err_e0 * 100. ; 
  434. 'SINON' ; 
  435.      err_e0 = 0. ; 
  436. 'FINSI' ; 
  437. 'SI' ( sig_0 '>' 0. ) ; 
  438.      err_s0 = ( sig_th0 - sig_0 ) / sig_0 ; 
  439.      err_s0 = 'ABS' err_s0 ; 
  440.      err_s0 = err_s0 * 100. ; 
  441. 'SINON' ; 
  442.      err_s0 = 0. ; 
  443. 'FINSI' ; 
  444. 'SI' ( sig_th0 '>' ( ( ( 9. * sig0 ) + SIG1 ) / 10. ) ) ; 
  445.      err_b0 = ( bet_th0 - bet_0 ) / bet_0 ; 
  446.      err_b0 = 'ABS' err_b0 ; 
  447.      err_b0 = err_b0 * 100. ; 
  448. 'SINON' ; 
  449.      err_b0 = 0. ; 
  450. 'FINSI' ; 
  451. 'SI' ( err_e0 '>' 0. ) ; 
  452.      err_r0 = ( rho_th0 - rho_0 ) / rho_0 ; 
  453.      err_r0 = 'ABS' err_r0 ; 
  454.      err_r0 = err_r0 * 100. ; 
  455. 'SINON' ; 
  456.      err_r0 = 0. ; 
  457. 'FINSI' ; 
  458. *  
  459. err_e1 = err_e1 'ET' ( 'PROG' err_e0 ) ; 
  460. err_s1 = err_s1 'ET' ( 'PROG' err_s0 ) ; 
  461. err_b1 = err_b1 'ET' ( 'PROG' err_b0 ) ; 
  462. err_r1 = err_r1 'ET' ( 'PROG' err_r0 ) ; 
  463. * 
  464. 'SI' ( err_e0 '>' 2.5E-2 ) ; 
  465.    'ERRE' 5 ; 
  466. 'SINON' ; 
  467.    'ERRE' 0 ; 
  468. 'FINSI' ; 
  469. 'SI' ( err_s0 '>' 2.5E-2 ) ; 
  470.    'ERRE' 5 ; 
  471. 'SINON' ; 
  472.    'ERRE' 0 ;  
  473. 'FINSI' ; 
  474. 'SI' ( err_b0 '>' 6.5 ) ; 
  475.    'ERRE' 5 ; 
  476. 'SINON' ; 
  477.    'ERRE' 0 ;  
  478. 'FINSI' ; 
  479. 'SI' ( err_r0 '>' 2.5E-2 ) ; 
  480.    'ERRE' 5 ; 
  481. 'SINON' ; 
  482.    'ERRE' 0 ;  
  483. 'FINSI' ; 
  484. * 
  485.   'MESS' ' Fin du calcul de restitution, temps ' temp0 ;  
  486. *  'MESS' '  ' ; 
  487. *   
  488. 'FIN' bloc1 ; 
  489. *  
  490. *  Trace  
  491. *  
  492. trac1 = 'EVOL' 'ROUG' 'MANU' 'epsp' depx1 'sigm' sigx1 ; 
  493. TAB = 'TABLE' ;
  494. TAB.1 = 'MARQ CROI ' ;
  495. TAB.2 = 'MARQ PLUS ' ;
  496. TAB.'TITRE' = 'TABLE' ;
  497. TAB.'TITRE'. 1 = MOT 'CALC' ;
  498. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  499. *'DESSIN' ( trac1 'ET' trac00 )
  500. *'TITR' 'Contraintes - Deformation plastique'
  501. *'LEGE' TAB ;
  502. trac2 = 'EVOL' 'ROUG' 'MANU' 'epsp' depx2 'sigm' sigx2 ;  
  503. TAB = 'TABLE' ; 
  504. TAB.1 = 'MARQ CROI ' ;
  505. TAB.2 = 'MARQ PLUS ' ;
  506. TAB.'TITRE' = 'TABLE' ;
  507. TAB.'TITRE'. 1 = MOT 'MAXI' ;
  508. TAB.'TITRE'. 2 = MOT 'MINI' ;
  509. *'DESSIN' ( trac1 'ET' trac2 )
  510. *'TITR' 'Comparaison mini - maxi de la contrainte equivalente '
  511. *'LEGE' TAB ;
  512. trac3 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'sigm' sig_th1 ; 
  513. trac4 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_1 'sigm' sig_1 ;  
  514. TAB = 'TABLE' ;
  515. TAB.1 = 'MARQ CROI ' ;
  516. TAB.2 = 'MARQ PLUS ' ;
  517. TAB.'TITRE' = 'TABLE' ;
  518. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  519. TAB.'TITRE'. 2 = MOT 'CALCUL' ;
  520. *'DESSIN' ( trac3 'ET' trac4 )
  521. *'TITR' 'Comparaison calcul theorie '
  522. *'LEGE' TAB ;
  523. trac5 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_s1 ;
  524. trac6 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_e1 ;
  525. TAB = 'TABLE' ;
  526. TAB.1 = 'MARQ CROI ' ;
  527. TAB.2 = 'MARQ PLUS ' ;
  528. TAB.'TITRE' = 'TABLE' ;
  529. TAB.'TITRE'. 1 = MOT 'CONTRAINTE' ;
  530. TAB.'TITRE'. 2 = MOT 'DEFORMATION' ;
  531. *'DESSIN' ( trac5 'ET' trac6 )
  532. *'TITR' ' Evolution des erreurs '
  533. *'LEGE' TAB ;
  534. trac7 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_b1 ;
  535. trac8 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'erreur%' err_r1 ;
  536. TAB = 'TABLE' ;
  537. TAB.1 = 'MARQ CROI ' ;
  538. TAB.2 = 'MARQ PLUS ' ;
  539. TAB.'TITRE' = 'TABLE' ;
  540. TAB.'TITRE'. 1 = MOT 'ENDOMMAGAMENT' ;
  541. TAB.'TITRE'. 2 = MOT 'DENSITE' ;
  542. *'DESSIN' ( trac7 'ET' trac8 )
  543. *'TITR' ' Evolution des erreurs '
  544. *'LEGE' TAB ;
  545. TAB = 'TABLE' ;
  546. TAB.1 = 'MARQ CROI ' ;
  547. TAB.2 = 'MARQ PLUS ' ;
  548. TAB.'TITRE' = 'TABLE' ;
  549. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  550. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  551. *'DESSIN' ( trac3 'ET' trac00 )
  552. *'TITR' 'Comparaison theorie traction '
  553. *'LEGE' TAB ;
  554. trac10 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'beta'bet_1 ; 
  555. trac11 = 'EVOL' 'ROUG' 'MANU' 'epsp' def_th1 'beta'bet_th1 ; 
  556. TAB = 'TABLE' ;
  557. TAB.1 = 'MARQ CROI ' ;
  558. TAB.2 = 'MARQ PLUS ' ;
  559. TAB.'TITRE' = 'TABLE' ;
  560. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  561. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  562. *'DESSIN' ( trac11 'ET' trac10 )
  563. *'TITR' 'Comparaison endommagement theorique et calcule '
  564. *'LEGE' TAB ;
  565. TAB = 'TABLE' ; 
  566. TAB.1 = 'MARQ CROI ' ;
  567. TAB.2 = '  ' ;
  568. TAB.'TITRE' = 'TABLE' ;
  569. TAB.'TITRE'. 1 = MOT 'THEORIE' ;
  570. TAB.'TITRE'. 2 = MOT 'TRACTION' ;
  571. *'DESSIN' ( trac4 'ET' trac00 )
  572. *'TITR' 'Comparaison calcul traction '
  573. *'LEGE' TAB ;
  574. *  
  575. *    
  576. 'FIN' ;  
  577.  
  578.  
  579.  
  580.  
  581.  
  582.  
  583.  
  584.  
  585.  
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  
  595.  

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