Télécharger la plaquette Cast3M


fr en

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 2  ;
  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 0. pas 2   4.;
  47. tb . chargement = ch;
  48. tb . modele = mo;
  49. tb . caracteristiques = ma;
  50. tb . temps_calcules = lit;                                           
  51. tb . precision      = 1e-8 ;
  52. *  le predicteur hpp marche super bien vu qu'il n'y a pas d'effort
  53. tb . predicteur     = hpp  ;
  54. pasapas tb;       
  55. co1 = tb . contraintes . 2 ;
  56. cma = maxi co1;
  57. mess ' contrainte maximum ' cma;
  58. si ( cma  >  1.e-2); 
  59.        erre  5;
  60. finsi;
  61. mess '************* test correct  **************';
  62. fin;
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  

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