Télécharger calp2.dgibi

Retour à la liste

Numérotation des lignes :

  1. *----------------------------------------------------------------------*
  2. * Cas-test du calcul de VMIS et des contraintes en peau des poutres *
  3. * avec CALP. *
  4. * *
  5. * On soumet l'extremite d'une poutre orientee suivant l'axe Ox a un *
  6. * deplacement suivant Oy+Oz, l'autre extremite etant encastree. *
  7. * *
  8. * On verifie le calcul de la contrainte de von Mises et de la compo- *
  9. * -sante SMXX des contraintes de peau calculee avec CALP. *
  10. * *
  11. * Le calcul est realise en DIME 3 ELEM SEG2. *
  12. *----------------------------------------------------------------------*
  13.  
  14. opti dime 3 elem seg2 ;
  15.  
  16. * Maillage :
  17. O1 = 0 0 0 ;
  18. X1 = 1 0 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. dx1 = 1. ;
  26. dy1 = 0.03 ;
  27. dz1 = 0.05 ;
  28. iry1 = 1. ;
  29. irz1 = 1. ;
  30. tor1 = 1. ;
  31. ma1 = mate mo1 youn 20.e9 nu 0.2 sect sec1 dx dx1 dy dy1 dz dz1 inry iry1 inrz irz1 tors tor1 ;
  32.  
  33. * CL / Chargement :
  34. cl1 = bloq depl rota O1 ;
  35. cl2 = bloq uy uz P2 ;
  36.  
  37. dcl2 = depi cl2 0.01 ;
  38.  
  39. * Resolution :
  40. K1 = rigi mo1 ma1 ;
  41. u1 = reso (K1 et cl1 et cl2) dcl2 ;
  42.  
  43. sig1 = sigm mo1 ma1 u1 ;
  44.  
  45. * VMIS :
  46. vmi1 = vmis mo1 ma1 sig1 ;
  47. *list vmi1 ;
  48.  
  49. * Calcul de VMIS selon la formule indiquee dans sa notice :
  50. * Seq = ((EFFX / SECT)**2 + (MOMX * DX / TORS)**2 + (MOMY * DZ / INRY)**2 + (MOMZ * DY / INRZ)**2)**0.5
  51. Seq1 = ((exco sig1 effx scal) / sec1 ** 2)
  52. + ((exco sig1 momx scal) * dx1 / tor1 ** 2)
  53. + ((exco sig1 momy scal) * dz1 / iry1 ** 2)
  54. + ((exco sig1 momz scal) * dy1 / irz1 ** 2) ;
  55. Seq1 = Seq1 ** 0.5 ;
  56. *List Seq1 ;
  57.  
  58. * Erreur calcul VMIS :
  59. ERR1 = maxi abs (vmi1 - seq1) ;
  60. list err1 ;
  61.  
  62. * SMXX en peau avec CALP :
  63. smp1 = calp sig1 mo1 ma1 ;
  64. smp1 = exco smp1 smxx smxx ;
  65. *list (exco smp1 smxx scal) ;
  66.  
  67. * Calcul de la contrainte de peau smxx ;
  68. * SMXX = ((EFFX / SECT) + (MOMY * DZ / INRY) - (MOMZ * DY / INRZ)
  69. smp2 = ((exco sig1 effx smxx) / sec1)
  70. + ((exco sig1 momy smxx) * dz1 / iry1)
  71. - ((exco sig1 momz smxx) * dy1 / irz1) ;
  72.  
  73. * Erreur calcul peau
  74. ERR2 = maxi abs (smp1 - smp2) ;
  75. list err2 ;
  76.  
  77. * Test de validation :
  78. ERR0 = maxi abs (prog err1 err2) ;
  79.  
  80. xpr1 = (vale prec) ** 0.5 ;
  81. si (err0 > xpr1) ;
  82. erre 5 ;
  83. fins ;
  84.  
  85. fin ;
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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