Télécharger la plaquette Cast3M


fr en

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