Télécharger gdep2ma.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gdep2ma.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * portes tournant autour d'un axe maillé en massif
  6. *
  7. opti dime 2 elem TRI3 ; dens 4;
  8. pa = 0 0 ; pb = 4 0 ; liab = pa d pb;
  9. suto = liab trans 4 ( 0 10);
  10. pq = cote 3 suto ;
  11. opti dime 3; optio elem cub8;
  12. voto = suto volu tran 1 (0 0 1);
  13.  
  14.  
  15. mo= voto model mecanique elastique pri6;
  16. EP = 1.;
  17. ma= mater mo young 10000 nu 0.3 ;
  18. bli = bloq pa UX ;
  19. blf = bloq liab UY UZ ;
  20. blro = bloq suto RY RZ ; bl = bli et blf et blro;
  21. bly = bloq pq uy;
  22. blz = bloq pq uz;
  23. ang = 11.25;
  24. gtz =10. *( prog 0. ( sin ang ) ( sin ( 2. * ang) )
  25. ( sin ( 3. * ang) )
  26. ( sin ( 4. * ang)) ( sin ( 5. * ang)) ( sin ( 6. * ang))
  27. ( sin ( 7. * ang)) ( sin ( 8. * ang)) ( sin ( 9. * ang))
  28. ( sin ( 10. * ang)) ) ;
  29.  
  30. gty =10. *( prog 0. ( cos ang - 1.) ( cos ( 2. * ang) - 1.)
  31. ( cos ( 3. * ang) - 1.)
  32. (cos (4.*ang) - 1.) (cos (5.*ang) - 1.) (cos (6.*ang) - 1.)
  33. (cos (7.*ang) - 1.) (cos (8.*ang) - 1.) (cos (9.*ang) - 1.)
  34. (cos (10.*ang) - 1.) ) ;
  35. fz = depimp blz 1. ;
  36. fy = depimp bly 1.;
  37. ft = prog 0. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10.;
  38. evz = evol manu abs ft ord gtz;
  39. evy = evol manu abs ft ord gty;
  40. chz = charge dimp fz evz;chy = charge dimp fy evy;
  41. ch = chz et chy;
  42. tb = table;
  43. tb . grands_deplacements = vrai;
  44. tb . blocages_mecaniques = blz et bly et bl;
  45. *tb . automatique = vrai; tb . autocrit = 100000;
  46. lit = prog 1. 2. 3. 4. ;
  47. tb . chargement = ch;
  48. tb . modele = mo;
  49. tb . caracteristiques = ma;
  50. tb . temps_calcules = lit;
  51. tb . precision = 1e-8;
  52. pasapas tb;
  53. co1 = tb . contraintes . 4;
  54. cma = maxi co1;
  55. mess ' contrainte maximum ' cma;
  56. si ( cma > 1.e-2);
  57. erre 5;
  58. finsi;
  59. mess '************* test correct **************';
  60. fin;
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  

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