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.5548077e7;
  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.  

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