Télécharger la plaquette Cast3M


fr en

Télécharger norton_tra1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  norton_tra1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************************
  5. *                                                *
  6. *  TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT  *
  7. *  --------------------------------------------  *
  8. *              DE MATERIAU                       *
  9. *              -----------                       *
  10. *                                                *
  11. *  LOI DE COMPORTEMENT DE FLUAGE DE :            *
  12. *              NORTON                            *
  13. *                                                *
  14. *  TEST POUR DES ELEMENTS MASSIFS                *
  15. *  ------------------------------                *
  16. *                                                *
  17. *  MAILLAGE:                                     *
  18. *   UNE CUBE DE COTE L=1 M                       *
  19. *                                                *
  20. *  CHARGEMENT:                                   *
  21. *   ESSAI DE TRACTION                            *
  22. *   FORCES IMPOSEES                              * 
  23. *                                                * 
  24. ************************************************** 
  25. opti echo 0 dime 3 elem cub8 ;  
  26. *  
  27. *  Maillage  
  28. *  
  29. l1 = 1. ; 
  30. l2 = .1 ; 
  31. p1 = 0. 0. 0. ; 
  32. p2 = l2 0. 0. ; 
  33. p3 = l2 l2 0. ; 
  34. p4 = 0. l2 0. ; 
  35. p10 = 0. 0. l1 ;
  36. p20 = l2 0. l1 ;
  37. p30 = l2 l2 l1 ;
  38. p40 = 0. l2 l1 ; 
  39. li1 = d p1 1 p2 ; 
  40. li2 = d p2 1 p3 ; 
  41. li3 = d p3 1 p4 ;  
  42. li4 = d p4 1 p1 ; 
  43. li10 = d p10 1 p20 ;
  44. li20 = d p20 1 p30 ;
  45. li30 = d p30 1 p40 ; 
  46. li40 = d p40 1 p10 ;
  47. sur1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ; 
  48. sur10 = 'SURF' ( li10 'ET' li20 'ET' li30 'ET' li40 ) 'PLAN' ; 
  49. mail1 = sur1 'VOLU' 1 sur10 ; 
  50. *  
  51. *  Modele 
  52. *  
  53. AF1 = 1.E-45 ; 
  54. AF2 = 5. ; 
  55. AF3 = 1.2 ; 
  56. SMAX = 2.1E8 ; 
  57. youn0 = 2.1E11 ; 
  58. nu0 = .3 ; 
  59. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
  60.        'FLUAGE' 'NORTON' ; 
  61. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'AF1' AF1 
  62.        'AF2' AF2 'AF3' AF3 'SMAX' SMAX ; 
  63. *  
  64. *  Conditions aux limites 
  65. *  
  66. cl1 = 'BLOQ' sur1 'UZ' ; 
  67. cl2 = 'BLOQ' ( li1 'ET' li10 ) 'UY' ; 
  68. cl3 = 'BLOQ' ( li4 'ET' li40 ) 'UX' ; 
  69. cl0 = cl1 'ET' cl2 'ET' cl3 ; 
  70. *  
  71. *  Chargement 
  72. *  
  73. depmax = 5.E5 ;  
  74. dep2 = 'FORCE' sur10 'FX' 0. 'FY' 0. 'FZ' depmax ;  
  75. dep0 = dep2 ; 
  76. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-2 1000. )
  77. y ( 'PROG' 0. 1. 1. ) ; 
  78. cha0 = 'CHAR' 'MECA' dep0 ev0 ; 
  79. *    
  80. dt0 = .05 ;
  81. dt1 = 2. * dt0 ;
  82. dt2 = 2. * dt1 ;
  83. dt3 = 5. * dt2 ;
  84. tfin0 = 50. ; 
  85. *  
  86. *  Appel à PASAPAS 
  87. *  
  88. ta1 = 'TABLE' ; 
  89. ta1 .'MODELE' = mod0 ; 
  90. ta1 .'CARACTERISTIQUES' = mat0 ; 
  91. ta1 .'CHARGEMENT' = cha0 ; 
  92. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  93. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' dt0 dt1 ) 'ET'
  94. ( 'PROG' ( dt1 + dt0 ) 'PAS' dt1 dt2 ) 'ET'
  95. ( 'PROG' ( dt2 + dt1 ) 'PAS' dt2 dt3 ) 'ET'
  96. ( 'PROG' ( dt3 + dt2 ) 'PAS' dt3 tfin0 ) ; 
  97. * 
  98. TMASAU=table;
  99. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  100. TMASAU .'DEFTO'=VRAI;
  101. TMASAU .'DEFIN'=VRAI; 
  102. PASAPAS ta1 ;  
  103. * 
  104. *  Post traitement  
  105. *  
  106. dim0 = 'DIME' ta1 .'TEMPS' ; 
  107. eps10 = 'PROG' 0. ; 
  108. sigm10 = 'PROG' 0. ; 
  109. def10 = 'PROG' 0. ; 
  110. def_th2 = 'PROG' 0. ;
  111. sig_th2 = 'PROG' 0. ;
  112. err_e2 = 'PROG' 0. ;
  113. err_s2 = 'PROG' 0. ; 
  114. temp2 = 'PROG' 0. ; 
  115. rac0 = 1. ;  
  116. H0 = youn0 ; 
  117. * 
  118. prob0 = vrai ; 
  119. 'REPETER' bloc0 ( dim0 - 1 ) ; 
  120.   i1 = &bloc0 ; 
  121. *  
  122. *  Résultats calculés  
  123. * 
  124.   def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i1 ;  
  125.   depl0 = ta1 .'DEPLACEMENTS'.i1 ; 
  126.   sigm0 = ta1 .'CONTRAINTES'.i1 ;  
  127.   temp0 = ta1 .'TEMPS'.i1 ; 
  128.   var0 = ta1 .'VARIABLES_INTERNES'.i1 ; 
  129.   eps0 = 'EPSI' mod0 depl0 ; 
  130.   eps1 = 'EXTR' ( 'EXCO' eps0 'EPZZ' 'SCAL' ) 'SCAL' 1 1 1 ; 
  131.   sigm1 = 'EXTR' ( 'EXCO' sigm0 'SMZZ' 'SCAL' ) 'SCAL' 1 1 1 ; 
  132.   def1 = 'EXTR' ( 'EXCO' def0 'EIZZ' 'SCAL' ) 'SCAL' 1 1 1 ; 
  133.   eps10 = eps10 'ET' ( 'PROG' eps1 ) ; 
  134.   sigm10 = sigm10 'ET' ( 'PROG' sigm1 ) ; 
  135.   sigeq0 = 'VMIS' sigm0 mod0 ; 
  136.   sigeq1 = 'EXTR' sigeq0 'SCAL' 1 1 1 ; 
  137.   def10 = def10 'ET' ( 'PROG' def1 ) ; 
  138.   temp2 = temp2 'ET' ( 'PROG' temp0 ) ; 
  139.   var1 = 'EXTR' ( 'EXCO' var0 'EPSE' 'SCAL' ) 'SCAL' 1 1 1 ; 
  140.   sig0 = 'CHAN' 'CHPO' sigm0 mod0 ;   
  141.   sigx = 'EXCO' sig0 'SMXX' ; 
  142.   sigy = 'EXCO' sig0 'SMYY' ; 
  143.   sigz = 'EXCO' sig0 'SMZZ' ; 
  144. *titr ' Contraintes suivant X' ; 
  145. *trace sigx mail1 ; 
  146. *titr ' Contraintes suivant Y' ; 
  147. *trace sigy mail1 ; 
  148. *titr ' Contraintes suivant Z' ; 
  149. *trace sigz mail1 ; 
  150. * 
  151. *  Résultats analytiques  
  152. * 
  153. y2 = ( ( rac0 * depmax / ( l2 * l2 ) ) ** ( AF2 ) ) ; 
  154. y2 = y2 * ( temp0 ** AF3 ) ; 
  155. x2 = y2 * AF1 * rac0 ; 
  156. *  
  157. def_th1 = x2 ; 
  158. sig_th1 = H0 * ( eps1 - x2 ) ; 
  159. *  
  160. def_th2 = def_th2 'ET' ( 'PROG' def_th1 ) ; 
  161. sig_th2 = sig_th2 'ET' ( 'PROG' sig_th1 ) ;  
  162. *mess ' sig calculée      =' sigm1 ; 
  163. *mess ' sig theorique     =' sig_th1 ; 
  164. *mess ' def calculée      =' def1 ; 
  165. *mess ' def théorique     =' def_th1 ; 
  166. *  
  167. *  Erreurs 
  168. * 
  169. 'SI' ( def1 '>' ( 1.E-1 * eps1 ) ) ; 
  170.    err_e1 = ( 'ABS' ( def1 - def_th1 ) ) / def1 ; 
  171.    err_e1 = err_e1 * 100. ; 
  172. 'SINON' ; 
  173.    err_e1 = 0. ; 
  174. 'FINSI' ;  
  175. 'SI' ( sigm1 '>' 0. ) ; 
  176.    err_s1 = ( 'ABS' ( sigm1 - sig_th1 ) ) / sigm1 ;
  177.    err_s1 = err_s1 * 100. ; 
  178. 'SINON' ; 
  179.    err_s1 = 0. ; 
  180. 'FINSI' ; 
  181. * 
  182. err_e2 = err_e2 'ET' ( 'PROG' err_e1 ) ; 
  183. err_s2 = err_s2 'ET' ( 'PROG' err_s1 ) ; 
  184. *  
  185. 'SI' ( ( err_e1 '>' .12 ) 'OU' ( err_s1 '>' .12 ) ) ; 
  186.    prob0 = faux ; 
  187. 'FINSI' ; 
  188. *   
  189. 'MESS' ' Fin du calcul de restitution au temps ' temp0 ; 
  190. *mess ' ' ; 
  191. * 
  192. 'FIN' bloc0 ; 
  193. *  
  194. defo0 = 'DEFO' mail1 depl0 0. ;
  195. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  196. *titr ' deformee finale' ;
  197. *trace ( defo0 'ET' defo1 ) ;
  198. * 
  199. 'SI' ( prob0 'EGA' faux ) ; 
  200.    'ERRE' 5 ; 
  201. 'SINON' ; 
  202.    'ERRE' 0 ; 
  203. 'FINSI' ; 
  204. * 
  205. *  Tracés  
  206. * 
  207. evo1 = 'EVOL' 'MANU' 'epsp' def10 'sigm' sigm10 ; 
  208. evo2 = 'EVOL' 'MANU' 'epsp' def_th2 'sigm' sig_th2 ; 
  209. evo3 = 'EVOL' 'MANU' 'epsp' def10 'erreur_%' err_e2 ; 
  210. evo4 = 'EVOL' 'MANU' 'epsp' def10 'erreur_%' err_s2 ; 
  211. evo5 = 'EVOL' 'MANU' 'temps' temp2 'epsp' def10 ;  
  212. evo6 = 'EVOL' 'MANU' 'temps' temp2 'epsp' def_th2 ; 
  213. TAB = 'TABLE' ; 
  214. TAB. 1 = 'MARQ CROI ' ;
  215. TAB. 2 = 'MARQ PLUS ' ;
  216. TAB.'TITRE' = 'TABLE' ;
  217. TAB.'TITRE'. 1 = MOT 'CALCUL' ;
  218. TAB.'TITRE'. 2 = MOT 'THEORIE' ;
  219. *'DESSIN' ( evo1 'ET' evo2 )
  220. *'TITR' ' Comparaison calcul theorie'
  221. *'LEGE' TAB ;
  222. TAB = 'TABLE' ; 
  223. TAB. 1 = 'MARQ CROI ' ;
  224. *TAB. 2 = 'MARQ PLUS ' ;
  225. TAB.'TITRE' = 'TABLE' ;
  226. TAB.'TITRE'. 1 = MOT 'DEFORMATION' ;
  227. *TAB.'TITRE'. 2 = MOT 'CONTRAINTE' ;
  228. *'DESSIN' ( evo3 )
  229. *'TITR' ' Evolution des erreurs (en pourcentage) '
  230. *'LEGE' TAB ;
  231. TAB = 'TABLE' ;
  232. TAB. 1 = 'MARQ CROI ' ;
  233. TAB. 2 = 'MARQ PLUS ' ;
  234. TAB.'TITRE' = 'TABLE' ;
  235. TAB.'TITRE'. 1 = MOT 'CALCUL' ;
  236. TAB.'TITRE'. 2 = MOT 'THEORIE' ;
  237. *'DESSIN' ( evo5 'ET' evo6 )
  238. *'TITR' ' Comparaison calcul theorie, deformation-temps'
  239. *'LEGE' TAB ;
  240. *  
  241. 'FIN' ;    
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  

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