Télécharger la plaquette Cast3M


fr en

Télécharger timf1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  timf1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph='N';
  5. saut page;
  6. opti echo 1;
  7. ********************************************************************
  8. * Etude d'une poutre sur appuis simple chargee de facon repartie 
  9. * (on en etudie que la demi-longueur)
  10. *
  11. *
  12. *  symetrie                           appuis simple
  13. *
  14. *     |-------------------------------------*
  15. *                                           ^
  16. *     p1                                    p2
  17. ********************************************************************
  18. * Notations:
  19. *
  20. * t =epaisseur
  21. * L =longueur totale
  22. * I =moment d'inertie
  23. *
  24. * q =charge par unite de longueur
  25. * e =module d'Young
  26. * nu=Poisson
  27. ********************************************************************
  28. * Solution theorique:
  29. *
  30. * deflection/(qL**4/24eI) =  5/16 + (t/L)**2(1+nu)(1/2) 
  31. ********************************************************************
  32. * Objectif du test:
  33. *
  34. * on cherche a verifier la solution theorique pour la suite
  35. * d'epaisseur reduite suivante:
  36. *
  37. * t/L = .0001 .001 .01 .1  1. 10.
  38. ********************************************************************
  39. * 
  40. * 1) option de calcul
  41. opti dime 3 elem seg2;
  42. *
  43. * 2) quelques donees et calcul de quelques constantes
  44. *     pour la solution theorique
  45. *
  46. q=1.; e=1.e10; nu=.3; L=1.; 
  47. *
  48. qL4s24e=(q * (L ** 4))/(24. * e);
  49. cinqs16=5. / 16.;
  50. nup1s2 = (nu + 1.)/2.;
  51. *
  52. * 3) maillage
  53. *
  54. nel=50;
  55. p1 = 0 0 0; p2 = (L/2.) 0 0;
  56. d1= p1 d nel p2;
  57. *
  58. * 4) modele et materiaux
  59. *
  60. modti=MODE d1 mecanique elastique section timo ;
  61. *
  62. * 5) chargement reparti sur la poutre
  63. *
  64. j=0; repe loopel nel; j=j+1;
  65.   elj=d1 elem j;
  66.   pee=elj point init;
  67.   pde=elj point final;
  68.   vec=pde moin pee;
  69.   xl=norm vec;
  70. *
  71.   valf1 = (q * xl) * 0.5 ;
  72.   valf2 = valf1;
  73. *
  74.   chrtij=(manu 'CHPO' pee 1 fz valf1 nature discret)
  75.      et  (manu 'CHPO' pde 1 fz valf2 nature discret);
  76. *
  77.   si (j ega 1); chrti=chrtij;
  78.   sinon;        chrti=chrtij et chrti;
  79.   finsi;
  80. *
  81. fin loopel;
  82. *
  83. * 6) bloquage
  84. *
  85. bloqtot=bloq d1 'UY' 'RX' 'RZ';
  86. bloqp1 =bloq p1 'UX' 'RY';
  87. bloqp2 =bloq p2 'UZ';
  88. *
  89. * 7) definition des differentes epaisseurs
  90. *
  91. tt=prog .0001 .001 .01 .1  1. 10.;
  92. nt=dime tt;
  93. *
  94. * 8) boucle sur les epaisseurs
  95. *
  96. cerre=0;
  97. j=0; repe loopep nt; j=j+1;
  98.   t=extr tt j;
  99. *
  100. ** 8.1) etablissement des caracteristiques geometriques
  101. *
  102.   b=1.;
  103.   opti elem qua4;
  104.   ps1=((-1)*(b/2)) ((-10)*(t/2)) 0.;
  105.   ps2=((+1)*(b/2)) ((-10)*(t/2)) 0.;
  106.   ps3=     0            t       0.;
  107.   sect=(ps1 d 1 ps2) tran 1 ps3;
  108. *
  109.   mose=MODE sect mecanique elastique quas; 
  110.   mase=MATE mose 'YOUN' e  'NU'  nu;
  111.   case= 'CARA' mose 'ALPY' 1. 'ALPZ' 1.;
  112.   mcse=mase et case;
  113. *  lrig lmas=rmsf mose mcse;
  114. *
  115.   mater=MATE modti 'MODS' mose 'MATS' mcse ;
  116.   mc=mater;
  117. *
  118.   section=b*t; inery=(b/12) * (t ** 3);
  119. *
  120. ** 8.2) calcul de la rigidite
  121. *
  122.   rigti=rigi modti mc;
  123.   rigtti=rigti et bloqtot et bloqp1 et bloqp2;
  124. *
  125. ** 8.3) solution
  126. *
  127.   depti=reso rigtti chrti;
  128. *
  129. ** 8.4) deflections theoriques en p1
  130. *
  131.   qL4s24ei=qL4s24e / inery;
  132.   theorie=cinqs16 + (((t / L) ** 2) * nup1s2);
  133.   calcul=(extr depti 'UZ' p1) / qL4s24ei;
  134.   erreur=abs ((calcul - theorie) / theorie);
  135. *
  136. ** 8.5) comparaison et message
  137. *
  138.   mess 'deflection pour t/L=' (t/L);
  139.   mess 'calcul,theorie,erreur:  ' calcul theorie erreur;
  140.   mess '***********************';
  141. *
  142. ** 8.6) erreur
  143. *
  144.   si (erreur > 5.d-4); cerre=cerre+1; finsi;
  145. *
  146. fin loopep;
  147. si (cerre > 0); erre 5;
  148. sinon;          erre 0; finsi;
  149. fin;
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  

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