Télécharger plas_incomp.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas_incomp.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. SAUT PAGE;
  5. ************************************************************************
  6. * CUBE EN TRACTION UNIAXIALE
  7. * Plastique parfait
  8. * verification de l'incompressibilite de l'ecoulement plastique
  9. *
  10. *
  11. optio elem cub8;
  12.  
  13. SAUT PAGE;
  14. TITRE 'CUBE EN TRACTION UNIAXIALE';
  15. OPTION DIME 3 ;
  16.  
  17. graph = 'N';
  18.  
  19. TEMPS ;
  20. DENS 1;
  21. *-----------------------------------------------------------------------
  22. ***** DEFINITION DE LA GEOMETRIE
  23. *-----------------------------------------------------------------------
  24.  
  25. p1 = 0 0 0; p2 = 0 1 0; l1 = p1 d 1 p2; s1 = l1 tran 1 (0 0 1);
  26. s2 = s1 plus (1 0 0); v1 = volu 1 s1 s2;
  27. en = enve v1;
  28. e1 = en elem 1; e2 = en elem 2; e3 = en elem 3; e4 = en elem 4; e5 = en elem 5; e6 = en elem 6;
  29.  
  30. bl1 = bloq ux s1; bl2 = bloq ux s2;
  31. bl3 = bloq uy uz p1; bl4 = bloq uz p2;
  32. bl5 = bloq uz (e2 et e5);
  33. blt = bl1 et bl2 et bl3 et bl4 et bl5;
  34.  
  35.  
  36. *-----------------------------------------------------------------------
  37. * VALEURS IMPOSEES DES DEPLACEMENTS
  38. *-----------------------------------------------------------------------
  39. di = depi bl1 -1 ;
  40. *-----------------------------------------------------------------------
  41. ***** CALCUL MODELE
  42. *-----------------------------------------------------------------------
  43. MO = MODE V1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' CUB8;
  44. MAT1 = MATE MO YOUN 2.5E11 SIGY 5.0E6 NU .499 ;
  45.  
  46.  
  47. *-----------------------------------------------------------------------
  48. ***** CHAMPS RELATIFS A L'OPERATEUR PASAPAS
  49. *-----------------------------------------------------------------------
  50. LP1 = PROG 0. 20.;
  51. LP2 = PROG 0. 20.;
  52. EV = EVOL MANU T LP1 F(T) LP2;
  53. CHA1 = CHAR 'DIMP' Di EV;
  54. LIS1 = PROG 0. PAS 1d-1 2 ;
  55. *-----------------------------------------------------------------------
  56. ***** CALCUL NON LINEAIRE
  57. *-----------------------------------------------------------------------
  58. TAB1 = TABLE;
  59. TAB1.'BLOCAGES_MECANIQUES' = blt;
  60. TAB1.'CARACTERISTIQUES' = MAT1;
  61. TAB1.'TEMPS_CALCULES' = LIS1;
  62. TAB1.'MODELE' = MO;
  63. TAB1.'CHARGEMENT' = CHA1;
  64. TAB1.'GRANDS_DEPLACEMENTS' = VRAI ;
  65. TAB1.'SOUS_INCREMENT' = 1 ;
  66. *TAB1.'LAGRANGIEN' = 'MOT' 'FIN_PAS';
  67. TAB1.'LAGRANGIEN' = 'MOT' 'REACTUALISE';
  68. TAB1.'PRECISION' = 1d-9 ;
  69. PASAPAS TAB1;
  70.  
  71. form tab1.wtable.for0;
  72. repete bou (dime tab1.deplacements - 2) ;
  73. dd = tab1.deplacements.&bou;
  74. form dd;
  75. mes = mesu volu v1;
  76. list mes;
  77. form tab1.wtable.for0;
  78. fin bou;
  79.  
  80. list (mes - 1);
  81. si (abs (mes - 1) > 1d-7) ; erreur 5; finsi;
  82. fin;
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  

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