Télécharger la plaquette Cast3M


fr en

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