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 'SOUCI';
  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.  
  159.  
  160.  

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