Télécharger la plaquette Cast3M


fr en

Télécharger rela_non_associee.dgibi

Retour à la liste

Numérotation des lignes : 

  1. *  fichier: rela_non_associee.dgibi
  2. *  test d'utilisation d'une relation non associee
  3. *
  4. *  on veut controler le deplacement d'un point intermediaire d'une barre en appliquant
  5. *  une force a l'extremite
  6. *  juste pour le fun on le traite en probleme unilateral et on teste la reaction
  7. *
  8. optio dime 2 elem qua4; dens 1;
  9.  
  10. *   maillage
  11. gau = (0 0)        ;
  12. pc = (1 0);
  13. droi = (3 0);
  14. e1 = gau d 1 pc;
  15. e2 = pc d 1 droi;
  16.  
  17. *   comportement
  18. mo = model (e1 et e2) mecanique elastique barr;
  19. ma = mater mo young 1 nu 1 sect 1;
  20.  
  21. *   raideur mecanique du systeme
  22. ri = rigi mo ma;
  23. bl0= (bloq ux    gau)                                  ;
  24. bl1= (bloq    uy gau)                                  ;
  25. bl2=                     (bloq uy pc)                  ;
  26. bl3=                                     (bloq uy droi);
  27. rit = ri et bl0 et bl1 et bl2 et bl3;
  28.  
  29. *   force de controle
  30. fo = force (-1 0) droi;
  31. fo = fo exco fx fx;
  32.  
  33. *   relation a imposer
  34. cp = manu chpo pc 1 ux 1;
  35. re = rela maxi cp dual fo ;
  36. list re;
  37.  
  38. *   second membre
  39. imp = (depi re 1) et (fo * -2) ;
  40.  
  41. *   resolution
  42. dep = resou imp (rit et re);
  43.  
  44. *  verif deplacement au noeud droi;
  45. *  on doit trouver 3
  46.  
  47. d3 = extr dep droi 'UX';
  48. si (((d3 - 3) abs) > 1d-10);
  49.   mess 'deplacement incorrect ' (d3 - 3);
  50.   erreur 5;
  51. finsi;
  52. pmult = extrai re maillage  mult ;
  53. l9 = (dep extrai valeur pmult 'LX') extrai 1;
  54. si (((l9 - 1) abs) > 1d-10);
  55.   mess 'reaction incorrecte ' (l9 - 1);
  56.   erreur 5;
  57. finsi;
  58. fin;
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  

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