Télécharger la plaquette Cast3M


fr en

Télécharger relacori.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  relacori.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph='N';
  5. ****************************************
  6. *
  7. *
  8. *  poutre maillée en coques et en massifs
  9. *
  10. * test de la relation de corps rigide
  11. ***********************************************
  13. opti elem cu20;
  14. p1 = 0 0 0 ; p2 = 0 1. 0 ;
  15. l1 = p1 d 2 p2 ;
  16. surf1 = l1 tran  5 (5 0 0 ) ;
  17.  
  18. l2 = cote 3 surf1 ;
  19. surf2 = l2 tran 1 (0 0 .5) ;
  20. surf2 = surf2 et (l2 tran 1 (0 0 -.5));
  21. vol1 = volu surf2 tran 8 (8. 0 0) ;
  22. vol2 = coul vol1 jaun ;
  23. mail1 = surf1 et vol2 ;
  24. pp = poin mail1 proche (13. .5 0) ;
  25.  
  26. * on tourne la poutre pour ne pas privilegier un repere particulier
  27. ax1 = -1 0.5 0 ; ax2 = -1 0.5 1 ; ax3 = -1 1 0 ;
  28. ax4 = 1 .5 0 ;
  29. depl mail1 tour 55. ax1 ax4 ;
  30. depl mail1 tour 30. ax1 ax2 ;
  31.  
  32. modl1 = mode surf1 mecanique coq8 ;
  33. mat1 = mate modl1 youn 2e11 nu .3 rho 8000. epai 1.;
  34. modl2 = mode vol2 mecanique cu20 ;
  35. mat2 = mate modl2 youn 2e11 nu .3 rho 8000.;
  36.  
  37. modlt = modl1 et modl2 ;
  38. matt = mat1 et mat2 ;
  39. rig1 = rigi modlt  matt ;
  40. mass1 = mass modlt  matt ;
  41. ch1 = manu chpo mail1 1 uZ 1. ;
  42. f1 = mass1 * ch1 ;
  43.  
  44. blo1 = bloq l1 depl rota ;
  45. blo2 = rela cori depl rota l2 surf2 ;
  46. rigt = blo1 et blo2 et rig1 ;
  47.  
  48. dep1 = reso f1 rigt ;
  49. def1 = (defo dep1 surf1 1000.) et
  50.                    (defo dep1 vol2 1000. jaun) ;
  51. si (ega 'O' graph);
  52. trac def1 cach;
  53. finsi;
  54.  
  55. x1 = extr dep1 uz pp ;
  56. ** solution analytique
  57. i1 = 1./12 ; q = 1.*8000;
  58. x2 = q*(13.**4)/2e11/8./i1 ;
  59. xe = (abs ((x1 - x2)/x2))*100 ;
  60.  
  61. mess 'deplacement suivant z';
  62. mess 'analytique'x2 'C2000' x1 'erreur' xe '%' ;
  63.  si (xe > .3) ; erre 5 ; sinon ; erre 0 ;finsi;
  64. fin ;
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  

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