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

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