Télécharger vibr15.dgibi

Retour à la liste

Numérotation des lignes :

  1. **********************************************************************
  2. * test elementaire de VIBC
  3. * syntaxe : [ A - \lambda I ] . X = 0
  4. * avec A = rigidite elementaire réelle symetrique
  5. **********************************************************************
  6.  
  7. opti dime 2 elem seg2;
  8. p1 = 0. 0.;
  9. p2 = 1. 0.;
  10. mesh = droi 1 p1 p2;
  11. prA= PROG 1.0000 0.5000 0.3333 0.2500
  12. 0.5000 1.0000 0.6667 0.5000
  13. 0.3333 0.6667 1.0000 0.7500
  14. 0.2500 0.5000 0.7500 1.0000;
  15. moprim = MOTS 'UX' 'UY';
  16. modual = MOTS 'FX' 'FY';
  17. A = MANU 'RIGI' mesh moprim prA;
  18.  
  19. tmodA = VIBC A;
  20.  
  21. * resultats https://www.wolframalpha.com:
  22. * λ_1≈2.53616, v_1≈(0.742622, 1.04602, 1.10333, 1)
  23. * λ_2≈0.848229, v_2≈(-1.73861, -0.402678, 0.645626, 1)
  24. * λ_3≈0.407833, v_3≈(0.881891, -1.47979, -0.0969911, 1)
  25. * λ_4≈0.207775, v_4≈(-0.132821, 0.693238, -1.47418, 1)
  26.  
  27. * on pourrait comparer, mais on va se contenter de calculer les résidus
  28.  
  29. prMaxRes1 = PROG;
  30. prMaxYvp = PROG;
  31.  
  32. nmod1 = 4;
  33. REPE Bmod1 nmod1;
  34. xvp = tmodA . 'MODES' . &Bmod1 . 'VALEUR_PROPRE';
  35. Yvp = tmodA . 'MODES' . &Bmod1 . 'DEFORMEE_MODALE';
  36. * MESS ' ';
  37. MESS '\l_'&Bmod1 ' = ' xvp;
  38. * LIST Yvp;
  39. Res1 = (A * Yvp) - (xvp * (EXCO Yvp moprim modual));
  40. * LIST Res1;
  41. prmaxYvp = prmaxYvp ET (MAXI 'ABS' Yvp);
  42. prmaxRes1 = prmaxRes1 ET (MAXI 'ABS' Res1);
  43. FIN Bmod1;
  44.  
  45.  
  46. xRes1 = (MAXI prmaxRes1) / (MAXI prmaxYvp);
  47. * rem : on met une tolérance très serrée
  48. xTol1 = 10. * (VALE 'PREC');
  49. MESS xRes1 xTol1;
  50.  
  51. SI (xRes1 > xTol1); ERRE 5; FINSI;
  52.  
  53. FIN;
  54.  
  55.  
  56.  

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