Télécharger la plaquette Cast3M


fr en

Télécharger ohno2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ohno2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Ohno2.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. *
  8. **************************************************
  9. *                                                *  
  10. *  TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT  *  
  11. *  --------------------------------------------  * 
  12. *              DE MATERIAU                       * 
  13. *              -----------                       *   
  14. *                                                *  
  15. *  LOI DE COMPORTEMENT VISCOPLASTIQUE DE:        * 
  16. *   OHNO (CHABOCHE MODIFIE OHNO)                 * 
  17. *                                                * 
  18. *  TEST POUR DES ELEMENTS COQUE MINCE            *  
  19. *  ----------------------------------            *  
  20. *                                                *   
  21. *  MAILLAGES:                                    * 
  22. *   PREMIER MAILLAGE: ELEMENTS MASSIFS           *  
  23. *   DEUXIEME MAILLAGE: ELEMENTS COQUE MINCE      *  
  24. *   UN CYLINDRE CREUX AXISYMETRIQUE              * 
  25. *   LONGUEUR      L =.1 M                        *  
  26. *   RAYON INTERNE R =.05 M                       * 
  27. *   EPAISSEUR     E =.0005 M                     * 
  28. *                                                *  
  29. *  CHARGEMENT:                                   * 
  30. *   CYLINDRE ENCASTRE EN SON EXTREMITE BASSE     * 
  31. *   FORCES IMPOSEES SUR L'EXTREMITE DU HAUT:     *   
  32. *      FORCES RADIALES DIRIGEES VERS L'EXTERIEUR * 
  33. *                                                * 
  34. **************************************************  
  35. *
  36. complet = faux ;
  37. opti echo 0 dime 2 elem qua8 mode axis ; 
  38. *  
  39. *  
  40. *  Maillage 
  41. *  
  42. rayon0 = .05 ; 
  43. epais0 = .0005 ; 
  44. l1 = epais0 + rayon0 ; 
  45. l2 = .1 ; 
  46. n1 = 1 ;
  47. n2 = 200 ;
  48. p1 = rayon0 0. ; 
  49. p2 = l1 0. ; 
  50. p3 = l1 l2 ; 
  51. p4 = rayon0 l2 ; 
  52. d1 = d p1 n1 p2 ; 
  53. d2 = d p2 n2 p3 ; 
  54. d3 = d p3 n1 p4 ; 
  55. d4 = d p4 n2 p1 ; 
  56. vol1 = 'DALL' d1 d2 d3 d4 ; 
  57. *  
  58. *    Modele de calcul 
  59. * 
  60. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  61.        'VISCOPLASTIQUE' 'OHNO' ; 
  62. mat0 = 'MATE' mod0 'YOUN' 7.34E10 'NU' 0.33   
  63.        'RHO' 7.8E3 'ALPHA' 1.E-5 'N' 24  
  64.        'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35   
  65.        'ALF' 2.E6 'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6 
  66.        'R1' 4 'A2' 80.E6 'C2' 45 'BET2' 58480.E6   
  67.        'R2' 4 'PHI' 1. 'B' 12 'GAMA' 2.E-7 'M' 2 
  68.        'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.06     
  69.        'EXP1' 2. 'EXP2' 2. ; 
  70. *  
  71. *    Conditions aux limites 
  72. *  
  73. cl1 = 'BLOQ' d1 'UR' 'UZ' 'RT' 'RZ' ; 
  74. cl0 = cl1 ;  
  75. *  
  76. *    Chargement 
  77. *  
  78. pres0 = 1.E4 ; 
  79. p01 = ( rayon0 + ( epais0 / 2. ) ) l2 ;  
  80. p00 = vol1 'POIN' 'PROC' p01 ;   
  81. dep1 = 'FORCE' 'FR' pres0 p00 ;  
  82. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 10000. ) y 
  83.       ( 'PROG' 0. 1. ) ; 
  84. cha0 = 'CHAR' 'MECA' dep1 ev0 ; 
  85. vec1 = 'VECT' dep1 1.E-5 'FR' 'FZ' 'ROUG' ; 
  86. *  
  87. *  Temps du calcul  
  88. *  
  89. dt0 = 20. ; 
  90. si complet ;
  91. tfin0 = 700. ;  
  92. sinon ;
  93. tfin0 = 200. ;
  94. finsi ;
  95.  
  96. *   
  97. *   
  98. *  Resolution par PASAPAS  
  99. *---------------------------------------
  100. *    
  101. ta1 = 'TABLE' ; 
  102. ta1 .'MODELE' = mod0 ; 
  103. ta1 .'CARACTERISTIQUES' = mat0 ; 
  104. ta1 .'BLOCAGES_MECANIQUES' = cl0 ; 
  105. ta1 .'CHARGEMENT' = cha0 ; 
  106. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ; 
  107. *  
  108. PASAPAS ta1 ; 
  109. *  
  110. *  
  111. *  Post-traitement 
  112. *-----------------------------------------
  113. * 
  114. dim0 = 'DIME' ( ta1 .'TEMPS' ) ; 
  115. char2 = 'PROG' 0. ; 
  116. depl2 = 'PROG' 0. ; 
  117. defi2 = 'PROG' 0. ; 
  118. cont1 = 'PROG' 0. ; 
  119. *  
  120. reac0 = ta1 .'REACTIONS'.( dim0 - 1 ) ;  
  121. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ;  
  122. vec2 = 'VECT' reac0 1.E-2 'FR' 'FZ' 'ROUG' ; 
  123. defo0 = 'DEFO' vol1 depl0 0. 'ROUG' ; 
  124. defo1 = 'DEFO' vol1 depl0 1. ; 
  125. reac1 = 'REDU' reac0 p1 ; 
  126. reac2 = 'REDU' reac0 p2 ; 
  127. vec01 = 'VECT' reac1 1.E-2 'FR' 'FZ' 'ROUG' ; 
  128. vec02 = 'VECT' reac2 1.E-2 'FR' 'FZ' 'VERT' ; 
  129. *  
  130. 'REPETER' bloc0 ( dim0 - 1 ) ; 
  131.   i0 = &bloc0 ; 
  132.   reac0 = ta1 .'REACTIONS'.i0 ; 
  133.   depl0 = ta1 .'DEPLACEMENTS'.i0 ; 
  134.   sigm0 = ta1 .'CONTRAINTES'.i0 ; 
  135.   defi0 = ta1 .'VARIABLES_INTERNES'.i0 ; 
  136.   temp0 = ta1 .'TEMPS'.i0 ;   
  137.   depl1 = 'EXTR' depl0 'UR' p3 ; 
  138.   depl2 = depl2 'ET' ( 'PROG' depl1 ) ; 
  139.   defi1 = 'MAXI' ( 'EXCO' defi0 'EPSE' ) ; 
  140.   defi2 = defi2 'ET' ( 'PROG' defi1 ) ; 
  141. *  cont0 =  'MAXI' ( 'EXCO' 'SMXX' sigm0 ) ;   
  142. *  cont1 = cont1 'ET' ( 'PROG' cont0 ) ;
  143.   char1 = pres0 * temp0 / tfin0 ; 
  144.   char2 = char2 'ET' ( 'PROG' char1 ) ; 
  145. 'FIN' bloc0 ; 
  146. * 
  147. opti echo 0 dime 2 elem seg2 mode axis ; 
  148. *  
  149. *  
  150. *  Maillage 
  151. *  
  152. epais0 = .0005 ; 
  153. rayon0 = .05 ; 
  154. l10 = .1 ; 
  155. n10 = 400 ; 
  156. p10 = rayon0 0. ; 
  157. p20 = rayon0 l10 ; 
  158. mail1 = d p10 n10 p20 ; 
  159. *  
  160. *    Modele de calcul 
  161. * 
  162. mod00 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  163.        'VISCOPLASTIQUE' 'OHNO' 'COQ2' ; 
  164. mat00 = 'MATE' mod00 'YOUN' 7.34E10 'NU' 0.33   
  165.        'RHO' 7.8E3 'N' 24  
  166.        'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35   
  167.        'ALF' 2.E6 'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6  
  168.        'R1' 4 'A2' 80.E6 'C2' 45 'BET2' 58480.E6   
  169.        'R2' 4 'PHI' 1. 'B' 12 'GAMA' 2.E-7 'M' 2 
  170.        'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.06 
  171.        'EXP1' 2. 'EXP2' 2. ; 
  172. car00 = 'CARA' mod00 'EPAI' epais0 'CALF' ( 2. / 3. ) ;  
  173. *  
  174. *    Conditions aux limites 
  175. *  
  176. cl10 = 'BLOQ' p10 'UR' 'UZ' 'RT' 'RZ' ; 
  177. cl00 = cl10 ;  
  178. *  
  179. *    Chargement 
  180. *  
  181. pres0 = 1.E4 ; 
  182. dep10 = 'FORCE' 'FR' pres0 p20 ;  
  183. ev00 = 'EVOL' 'MANU' temps ( 'PROG' 0. 10000. ) y 
  184.        ( 'PROG' 0. 1. ) ; 
  185. cha00 = 'CHAR' 'MECA' dep10 ev00 ; 
  186. vec1 = 'VECT' dep10 1.E-5 'FR' 'FZ' ; 
  187. *  
  188. *  Temps du calcul  
  189. *  
  190. dt00 = 20. ; 
  191. si complet ;
  192. tfin00 = 700. ;
  193. sinon ;
  194. tfin00 = 200. ;
  195. finsi ;  
  196. *   
  197. *   
  198. *  Resolution par PASAPAS  
  199. *---------------------------------------
  200. *    
  201. ta10 = 'TABLE' ; 
  202. ta10 .'MODELE' = mod00 ; 
  203. ta10 .'CARACTERISTIQUES' = mat00 'ET' car00 ; 
  204. ta10 .'BLOCAGES_MECANIQUES' = cl00 ; 
  205. ta10 .'CHARGEMENT' = cha00 ; 
  206. ta10 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ; 
  207. *  
  208. PASAPAS ta10 ; 
  209. *  
  210. *  
  211. *  Post-traitement 
  212. *-----------------------------------------
  213. * 
  214. dim10 = 'DIME' ta10 .'TEMPS' ; 
  215. depl12 = 'PROG' 0. ; 
  216. *  
  217. reac0 = ta10 .'REACTIONS'.( dim10 - 1 ) ;  
  218. depl0 = ta10 .'DEPLACEMENTS'.( dim10 - 1 ) ; 
  219. vec2 = 'VECT' reac0 1.E-2 'FR' 'FZ' 'ROUG' ; 
  220. defo0 = 'DEFO' mail1 depl0 0. ;
  221. defo1 = 'DEFO' mail1 depl0 100. ;
  222. *   
  223. 'REPETER' bloc0 ( dim10 - 1 ) ; 
  224.   i10 = &bloc0 ; 
  225.   depl10 = ta10 .'DEPLACEMENTS'.i10 ; 
  226.   temp10 = ta10 .'TEMPS'.i10 ; 
  227.   depl11 = 'EXTR' depl10 'UR' p20 ; 
  228.   depl12 = depl12 'ET' ( 'PROG' depl11 ) ; 
  229. 'FIN' bloc0 ; 
  230. *  
  231. *  Erreur   
  232. * 
  233. prob1 = faux ; 
  234. 'REPETER' bloc0 ( dim10 - 1 ) ;
  235.   i10 = &bloc0 ;
  236.   depl1 = 'EXTR' depl2 i10 ; 
  237.   depl11 = 'EXTR' depl12 i10 ; 
  238. * 
  239. 'SI' (depl1 'NEG' 0. ) ; 
  240.   err_e0 = 'ABS' ( ( depl1 - depl11 ) / depl1 ) ; 
  241.   err_e0 = err_e0 * 100. ; 
  242. 'SINON' ; 
  243.   err_e0 = 0. ; 
  244. 'FINSI' ; 
  245. *  
  246. *
  247.   'SI' ( err_e0 '>' 3.E-1 ) ;
  248.      prob0 = vrai ;
  249.   'SINON' ;
  250.      prob0 = faux ;
  251.   'FINSI' ;
  252.   'SI' (prob0 'EGA' vrai ) ; 
  253.      prob1 = vrai ; 
  254.   'FINSI' ; 
  255. *
  256. 'FIN' bloc0 ;
  257. *
  258. 'SI' ( prob1 'EGA' vrai ) ;
  259.    'ERRE' 5 ;
  260. 'SINON' ;
  261.    'ERRE' 0 ;
  262. 'FINSI' ;
  263. *  
  264. *  
  265. 'FIN' ;   
  266.  
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  

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