Télécharger la plaquette Cast3M


fr en

Télécharger grot1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  grot1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. opti dime 2 elem qua4 mode plan defo;
  5. opti echo 0;
  6.  
  7. ********************************************
  8. *   grandes rotations sur un element 2D-DP
  9. *    calcul elastique lineaire
  10. ********************************************
  11.  
  12.  
  13. l = 1000.;
  14. p1 = 0 l;
  15. p3 = 0 0;
  16.  
  17. d1 = droit 1 p3 p1;
  18. plaq = d1 tran 1 (l 0);
  19. p2 = (enve plaq) poin proc (l l) ;
  20. p4 = (enve plaq) poin proc (l 0);
  21.  
  22. *-------------------------------------------------
  23. *   creation du modele
  24. *-------------------------------------------------
  25.  
  26. objaf = MODE plaq mecanique elastique plastique parfait;
  27.  
  28. *-------------------------------------------------
  29. *   caracteristiques du materiau
  30. *-------------------------------------------------
  31.  
  32. mat = MATE objaf YOUNG 200000. NU 0.3 SIGY 10E6 ;
  33.  
  34. *--------------------------------------------------
  35. *   calcul des rigidites elementaires
  36. *             et 
  37. *   definition des blocages
  38. *--------------------------------------------------
  39.  
  40.  
  41. *    deplacement impose
  42. cdl1 = bloq UX p1;
  43. cdep1 = depi cdl1 (-1000.);
  44. cdl11 = bloq UY p1;
  45. cdep11 = depi cdl11 (-1000.);
  46.  
  47. cdl2 = bloq UX p2;
  48. cdep2 = depi cdl2 (-2000.);
  49. cdl21 = bloq UY p2;
  50. cdep21 = depi cdl21 0;
  51.  
  52. cdl3 = bloq UX p3;
  53. cdep3 = depi cdl3 0;
  54. cdl31 = bloq UY p3;
  55. cdep31 = depi cdl31 0;
  56.  
  57. cdl4 = bloq UX p4;
  58. cdl41 = bloq UY p4;
  59. cdep4 = depi cdl4 (-1000.);
  60. cdep41 = depi cdl41 (1000.);
  61.  
  62. cdlt = cdl1 et cdl11 et cdl2 et cdl21 et cdl3 et cdl31 
  63.        et cdl4 et cdl41;
  64.  
  65. fo1 = cdep1 et cdep11 et cdep2 et cdep21 et cdep3 et cdep31
  66.       et cdep4 et cdep41 ;
  67.  
  68. cdep41 = depi cdl41 1100.;
  69. cdep21 = depi cdl21 100.;
  70.  
  71. fo2 = cdep1 et cdep11 et cdep21 et cdep3 et cdep31
  72.       et cdep4 et cdep41 et cdep2;
  73.  
  74.  
  75. *   creation du chargement
  76. li1 = prog 0 19 20;
  77.  
  78.  
  79. li3 = prog 0 1 0;
  80. li4 = prog 0 0 1;
  81. li5 = prog 0 pas 19 19 pas 1 20;
  82. evt2 = evol manu t li1 f li3;
  83. evt3 = evol manu t li1 f li4;
  84.  
  85. cha1 = char 'DIMP' fo1 evt2;
  86. cha2 = char 'DIMP' fo2 evt3;
  87.  
  88. *----------------------------------------------------
  89. *             resolution avec pasapas
  90. *----------------------------------------------------
  91.  
  92. tab1 = table;
  93. tab1.'GRANDS_DEPLACEMENTS' = vrai;
  94. tab1.'BLOCAGES_MECANIQUES' = cdlt;
  95. tab1.'MODELE' = objaf;
  96. tab1.'CARACTERISTIQUES' = mat;
  97. tab1.'CHARGEMENT' = cha1 et cha2;
  98. tab1.'TEMPS_CALCULES' = li5;
  99. pasapas tab1;
  100.  
  101. pp = react cdl21 (tab1.'DEPLACEMENTS'.(2));
  102. ppa = 1.34310E+07;
  103. ppy = maxi pp;
  104. err1 = abs ((ppy - ppa)/ppa);
  105. SI (ERR1 <  1.E-2);
  106.     ERRE 0;
  107. SINO;
  108.     ERRE 5;
  109. FINSI;
  110.  
  111. fin;
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  

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