Télécharger fluage_fibre_norton_2.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * Comparaison du modele de fluage de Norton : *
  3. * - modele poutre a fibre VS modele massif *
  4. * - chargement uniaxial, traction puis compression *
  5. * - en force imposee *
  6. ************************************************************************
  7.  
  8.  
  9.  
  10. ** Options generales
  11. OPTI 'DIME' 3 'ELEM' 'CU20' 'ECHO' 0 ;
  12. itrac = FAUX ;
  13.  
  14.  
  15. ** Parametres geometrie (poutre a section rectangulaire)
  16. a = 0.05 ;
  17. b = 0.02 ;
  18. l = 1. ;
  19. se = a * b ;
  20.  
  21.  
  22. ** Nombres d'elements
  23. nea = 1 ;
  24. neb = 1 ;
  25. nel = 1 ;
  26.  
  27.  
  28. ** Parametres materiau
  29. yo = 1.E8 ;
  30. nu = 0.3 ;
  31. af1 = 3.E-13 ;
  32. af2 = 1.2 ;
  33. af3 = 1.4 ;
  34.  
  35.  
  36. ** Parametres chargement
  37. ftrac = 10000. ;
  38. tl = 200. ;
  39.  
  40.  
  41. ** Maillage et modele volumique
  42. p0 = 0. 0. 0. ;
  43. p1 = 0. a 0. ;
  44. l1 = DROI nea p0 p1 ;
  45. s1 = l1 TRAN neb (0. 0. b) ;
  46. p2 = s1 POIN 'PROC' (0. 0. b) ;
  47. v1 = s1 VOLU 'TRAN' nel (l 0. 0.) ;
  48. env1 = ENVE v1 ;
  49. are1 = ARET v1 ;
  50. p3 = v1 POIN 'PROC' (l 0. 0.) ;
  51. s2 = v1 FACE 2 ;
  52. mov = MODE v1 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'NORTON' ;
  53. mav = MATE mov 'YOUN' yo 'NU' nu 'AF1' af1 'AF2' af2 'AF3' af3 ;
  54.  
  55.  
  56. ** Maillage et modele de section
  57. OPTI 'ELEM' 'CUB8' ;
  58. lig1 = DROI nea ((-0.5 * a) (-0.5 * b) 0.) ((0.5 * a) (-0.5 * b) 0.) ;
  59. s3 = lig1 TRAN neb (0. b 0.) ;
  60. s3 = SURF (CONT s3) 'PLAN' ;
  61. mos = MODE s3 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'NORTON' 'QUAS' 'TRIS' ;
  62. mas = MATE mos 'YOUN' yo 'NU' nu 'AF1' af1 'AF2' af2 'AF3' af3 'SMAX' (yo / 1000.) 'ALPY' 0.66 'ALPZ' 0.66 ;
  63.  
  64.  
  65. ** Maillage et modele de poutre TIMO
  66. lf = DROI nel p0 p3 ;
  67. mop = MODE lf 'MECANIQUE' 'ELASTIQUE' 'SECTION' 'PLASTIQUE' 'SECTION' 'TIMO' ;
  68. map = MATE mop 'MODS' mos 'MATS' mas 'VECT' (0. 1. 0.) ;
  69.  
  70.  
  71. ** Conditions aux limites
  72. tl0 = 1.E-5 ;
  73. tl2 = 2. * tl ;
  74. tl3 = 3. * tl ;
  75. tl4 = 4. * tl ;
  76. evcha = EVOL 'MANU' (PROG 0. tl0 tl (tl + tl0) tl2 (tl2 + tl0) tl3 (tl3 + tl0) tl4)
  77. (PROG 0. 0.999 1. 0.001 0. -0.999 -1. -0.001 0.) ;
  78.  
  79. * pour le modele volumique
  80. * encastrement
  81. bl1 = (BLOQ 'UX' s1) ET (BLOQ 'UY' 'UZ' p0) ET (BLOQ 'UZ' p1) ;
  82. * force imposee
  83. ft = FSUR 'MASS' mov s2 ((ftrac / se) 0. 0.) ;
  84. chaft = CHAR 'MECA' ft evcha ;
  85.  
  86. * pour le modele TIMO
  87. * encastrement
  88. bl1p = BLOQ 'DEPL' 'ROTA' p0 ;
  89. * force imposee
  90. ftp = FORC (ftrac 0. 0.) p3 ;
  91. chaftp = CHAR 'MECA' ftp evcha ;
  92.  
  93.  
  94. ** Resolution
  95. xpas = tl / 20. ;
  96. ltc = PROG tl0 'PAS' xpas tl
  97. (tl + tl0) 'PAS' xpas tl2
  98. (tl2 + tl0) 'PAS' xpas tl3
  99. (tl3 + tl0) 'PAS' xpas tl4 ;
  100.  
  101. tv = TABL ;
  102. tv . 'MODELE' = mov ;
  103. tv . 'CARACTERISTIQUES' = mav ;
  104. tv . 'BLOCAGES_MECANIQUES' = bl1 ;
  105. tv . 'CHARGEMENT' = chaft ;
  106. tv . 'TEMPS_CALCULES' = ltc ;
  107. PASAPAS tv ;
  108. itp1 = TEMP 'ENTI' 'HORL' ;
  109.  
  110. tp = TABL ;
  111. tp . 'MODELE' = mop ;
  112. tp . 'CARACTERISTIQUES' = map ;
  113. chmsg0 = ZERO mos 'CONTRAIN' ;
  114. chmvi0 = ZERO mos 'VARINTER' ;
  115. tp . 'VARIABLES_INTERNES' = TABL ;
  116. tp . 'VARIABLES_INTERNES' . 0 = MANU 'CHML' mop 'VONS' chmsg0 'VAIS' chmvi0 'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  117. tp . 'BLOCAGES_MECANIQUES' = bl1p ;
  118. tp . 'CHARGEMENT' = chaftp ;
  119. tp . 'TEMPS_CALCULES' = ltc ;
  120. PASAPAS tp ;
  121. itp11 = TEMP 'ENTI' 'HORL' ;
  122.  
  123. MESS 'Temps horloge : ' ;
  124. MESS (CHAI '- modele massif :' ' ' itp1) ;
  125. MESS (CHAI '- modele poutre :' ' ' itp11) ;
  126.  
  127.  
  128. ** Post traitement
  129. n1 = DIME (tv . 'TEMPS') ;
  130. ltps = PROG ;
  131. * deplacements de l'extremite
  132. luv = EXTR (EVOL 'TEMP' tv 'DEPLACEMENTS' 'UX' p3) 'ORDO' ;
  133. lup = EXTR (EVOL 'TEMP' tp 'DEPLACEMENTS' 'UX' p3) 'ORDO' ;
  134. * force de reaction a l'encastrement
  135. lfv = PROG ;
  136. REPE b1 n1 ;
  137. i1 = &b1 - 1 ;
  138. ltps = ltps ET (tv . 'TEMPS' . i1) ;
  139. fv = 0. ;
  140. SI (NEG i1 0) ;
  141. fv = MAXI (EXCO 'FX' (RESU (REDU (tv . 'REACTIONS' . i1) s1))) ;
  142. FINSI ;
  143. lfv = lfv ET fv ;
  144. FIN b1 ;
  145. lfp = EXTR (EVOL 'TEMP' tp 'REACTIONS' 'FX' p0) 'ORDO' ;
  146. * defomration inelastique
  147. lev = PROG ;
  148. lep = PROG ;
  149. vol1 = MESU v1 ;
  150. REPE b1 n1 ;
  151. i1 = &b1 - 1 ;
  152. tps1 = tv . 'TEMPS' . i1 ;
  153. lev = lev ET ((INTG mov (tv . 'VARIABLES_INTERNES' . i1) 'EPSE') / vol1) ;
  154. uv = tv . 'DEPLACEMENTS' . i1 ;
  155. up = tp . 'DEPLACEMENTS' . i1 ;
  156. tabv = TABL ;
  157. tabv . 'DEPLACEMENTS' = TABL ;
  158. tabv . 'DEPLACEMENTS' . 1 = up ;
  159. tabv . 'VONS' = TABL ;
  160. tabv . 'VONS' . 1 = EXCO (tp . 'VARIABLES_INTERNES' . i1) 'VONS' ;
  161. tabv . 'VAIS' = TABL ;
  162. tabv . 'VAIS' . 1 = EXCO (tp . 'VARIABLES_INTERNES' . i1) 'VAIS' ;
  163. mail1 = POUT2MAS mop map 'GAUSS' tabv ;
  164. vips = tabv . 'VAIS_3D' . 1 ;
  165. mobid = MODE mail1 'MECANIQUE' ;
  166. cham1 = CHAN 'CHAM' mobid vips ;
  167. epsep = (INTG mobid cham1 'EPSE') * l / (NBEL lf) / se ;
  168. lep = lep ET epsep ;
  169. FIN b1 ;
  170.  
  171.  
  172. ** Analyse/trace des resultats
  173. tleg = TABL ;
  174. tleg . 1 = MOT 'MARQ LOSA NOLI' ;
  175. tleg . 2 = MOT 'MARQ ROND NOLI' ;
  176. tleg . 'TITRE' = TABL ;
  177. tleg . 'TITRE' . 1 = 'Modele 3D massif' ;
  178. tleg . 'TITRE' . 2 = 'Modele poutre fibre' ;
  179. MESS 'Ecart relatif max.' ;
  180. evuv = EVOL 'VERT' 'MANU' 'Temps' ltps 'Deplacement' luv ;
  181. evup = EVOL 'ROUG' 'MANU' 'Temps' ltps 'Deplacement' lup ;
  182. difu = (MAXI (ABS (lup - luv))) / (MAXI (ABS luv)) ;
  183. MESS '-- deplacement :' difu ;
  184. SI itrac ;
  185. DESS (evuv ET evup) 'TITR' 'Deplacement vs. Temps' 'LEGE' tleg ;
  186. FINSI ;
  187. evfv = EVOL 'VERT' 'MANU' 'Temps' ltps 'Force' lfv ;
  188. evfp = EVOL 'ROUG' 'MANU' 'Temps' ltps 'Force' lfp ;
  189. diff = (MAXI (ABS (lfp - lfv))) / (MAXI (ABS lfv)) ;
  190. MESS '-- effort :' diff ;
  191. SI itrac ;
  192. DESS (evfv ET evfp) 'TITR' 'Force vs. Temps' 'LEGE' tleg ;
  193. FINSI ;
  194. evev = EVOL 'VERT' 'MANU' 'TEMPS' ltps 'EPSE' lev ;
  195. evep = EVOL 'ROUG' 'MANU' 'TEMPS' ltps 'ESPE' lep ;
  196. dife = (MAXI (ABS (lep - lev))) / (MAXI (ABS lev)) ;
  197. MESS '-- def. fluage :' dife ;
  198. SI itrac ;
  199. DESS (evev ET evep) 'TITR' 'Deformation non lin. (EPSE) vs Temps' 'LEGE' tleg ;
  200. FINSI ;
  201.  
  202.  
  203. ** Erreur si l'ecart relatif est trop eleve
  204. lerr = PROG difu diff dife ;
  205. errmax = MAXI lerr ;
  206. SI (errmax > 6000.) ;
  207. MESS 'Echec du cas test !' ;
  208. ERRE 5 ;
  209. SINON ;
  210. MESS 'Succes du cas test !' ;
  211. FINSI ;
  212.  
  213.  
  214. FIN ;
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  

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