Télécharger calp1.dgibi

Retour à la liste

Numérotation des lignes :

  1. *----------------------------------------------------------------------*
  2. * Cas-test du calcul de VMIS dans le cas des poutres. *
  3. * *
  4. * On soumet l'extremite d'une poutre orientee suivant l'axe Ox a un *
  5. * deplacement suivant Oy, l'autre extremite etant encastree. *
  6. * *
  7. * On verifie le calcul de la contrainte de von Mises. *
  8. * *
  9. * Le calcul est realise en DIME 2 ELEM SEG2. *
  10. *----------------------------------------------------------------------*
  11.  
  12. opti dime 2 elem seg2 ;
  13.  
  14. * Maillage :
  15. O1 = 0 0 ;
  16. X1 = 1 0 ;
  17. L1 = O1 droi 10 (10 * X1) ;
  18. P2 = L1 poin 2 ;
  19.  
  20. * Modele / caracteristiques :
  21. mo1 = mode L1 mecanique pout ;
  22. sec1 = 0.0025 ;
  23. dy1 = 0.03 ;
  24. irz1 = 1. ;
  25. tor1 = 1. ;
  26. ma1 = mate mo1 youn 20.e9 nu 0.2 sect sec1 dy dy1 inrz irz1 ;
  27.  
  28. * CL / Chargement :
  29. cl1 = bloq depl rota O1 ;
  30. cl2 = bloq uy P2 ;
  31.  
  32. dcl2 = depi cl2 0.01 ;
  33.  
  34. * Resolution :
  35. K1 = rigi mo1 ma1 ;
  36. u1 = reso (K1 et cl1 et cl2) dcl2 ;
  37.  
  38. sig1 = sigm mo1 ma1 u1 ;
  39.  
  40. * VMIS :
  41. vmi1 = vmis mo1 ma1 sig1 ;
  42. *list vmi1 ;
  43.  
  44. * Calcul de VMIS selon la formule indiquee dans sa notice :
  45. * Seq = ((EFFX / SECT)**2 + (MOMZ * DY / INRZ)**2)**0.5
  46. Seq1 = ((exco sig1 effx scal) / sec1 ** 2)
  47. + ((exco sig1 momz scal) * dy1 / irz1 ** 2) ;
  48. Seq1 = Seq1 ** 0.5 ;
  49. *List Seq1 ;
  50.  
  51. * Erreur calcul VMIS :
  52. ERR1 = maxi abs (vmi1 - seq1) ;
  53. list err1 ;
  54.  
  55. xpr1 = (vale prec) ** 0.5 ;
  56. si (err1 > xpr1) ;
  57. erre 5 ;
  58. fins ;
  59.  
  60. fin ;
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  

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