Télécharger la plaquette Cast3M


fr en

Télécharger calp1.dgibi

Retour à la liste

Numérotation des lignes : 

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

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