Télécharger la plaquette Cast3M


fr en

Télécharger elas18.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  elas18.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. saut page;
  5. ********************************************************************
  6. * Etude d'une poutre encastree a une extremite et chargee a l'autre
  7. *    on néglige l'énergie associée à l'effort tranchant
  8. *    La flexion se fait dans le plan xoy
  9. *
  10. *  encastrement                       extremite chargee
  11. *
  12. *     |-------------------------------------*
  13. *                                           ^
  14. *     p1                                    p2
  15. ********************************************************************
  16. * Notations:
  17. *
  18. * h =hauteur
  19. * b = largeur
  20. * L =longueur totale
  21. * I = b * h**3 / 12
  22. * P =charge appliquee en P2
  23. * E =module d'Young
  24. * nu=Poisson
  25. ********************************************************************
  26. * Objectif du test:
  27. * on cherche à verifier la solution theorique en flexion pure
  28. * donnant la flèche en bout
  29. ********************************************************************
  30. * 
  31. * 1) option de calcul
  32. opti dime 3 elem seg2;
  33. *
  34. * 2) quelques donnees et calcul de quelques constantes
  35. *     pour la solution theorique
  36. *
  37. h=1.; E=1.e10; nu=.3; L=10.; b = 0.5 ;
  38. P = 50000. ;
  39. h3 = h ** 3 ;
  40. L3 = l ** 3 ;
  41. inerz = ( b * h3 ) /12. ;
  42. inery = ( (b ** 3 ) * h ) /12. ;
  43. *
  44. *
  45. * 3) maillage
  46. *
  47. nel=5;
  48. p1 = 0 0 0; p2 = L 0 0;
  49. d1= p1 d nel p2;
  50. *
  51. * 4) modele et materiaux
  52. *    et  caracteristiques geometriques
  53. *
  54. modti=MODE d1 mecanique elastique pout;
  55. mater=MATE modti 'YOUN' E 'NU' nu;
  56. carac= 'CARA' modti 'SECT' ( b * h )  'INRY' inery 'INRZ' inerz
  57.                     'TORS' (inery + inerz)   VECT ( 0 1 0 );
  58. mc=mater et carac;
  59. *
  60. * 5) chargement
  61. *
  62. chrti = force 'FY' P  p2 ;
  63. *
  64. * 6) bloquage
  65. *
  66. bloqtot=bloq depl rota p1;
  67. *
  68. *
  69. * 7) calcul de la rigidite
  70. *
  71.   rigti=rigi modti mc;
  72.   rigtti=rigti et bloqtot ;
  73. *
  74. * 8) solution
  75. *
  76.   depti=reso rigtti chrti;
  77. *
  78. * 9) calcul de la solution analytique
  79. *
  80.   theorie=  ( P * L3 ) / ( ( 3. * E ) * inerz ) ;
  81.   calcul= extr depti 'UY' p2 ;
  82.   erreur=abs ((calcul - theorie) / theorie);
  83. *
  84. * 10) comparaison et message
  85. *
  86.   mess 'calcul,theorie,erreur:  ' calcul theorie erreur;
  87.   mess '***********************';
  88.   si (erreur > 5.d-4);
  89.        erre 5;
  90.   sinon;
  91.        erre 0;
  92.   finsi;
  93. fin;
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  

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