Télécharger vibc.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ VIBC     NOTICE  BP208322  20/02/20    21:15:02     10528          
   2 :                                              DATE     20/02/20
   3 :                                              
   4 :   Operateur VIBC                           Voir aussi : VIBR PJBA
   5 :     --------------  
   6 :     --------------------
   7 :     | 1ere possibilite |
   8 :     --------------------
   9 :     
  10 :       BAS2 = VIBC  MASS1 RIG1 (AMOR1) (BAS1)  (ENT1);
  11 : 
  12 : 
  13 : 
  14 :     Objet :
  15 :     _______
  16 : 
  17 :     L'operateur VIBC recherche les valeurs propres complexes w (en Hz)
  18 :     et les vecteurs propres complexes X solutions de l'equation 
  19 :     fondamentale de la dynamique :
  20 :                         M q'' + C q' + K q = 0
  21 :                  avec   q(t) = X exp(i*2*pi*w*t)
  22 :                  
  23 :     Il resoud donc :
  24 :               [K + (i*2*pi*w)*C - (2*pi*w)**2 M] X = 0
  25 :     et fournit :
  26 :               X = X  + i X       et    w = w  + i w
  27 :                    R      I                 R      I
  28 :                    
  29 :     L'algorithme utilise est le QZ, ce qui implique que la taille du 
  30 :     probleme doit etre petite (on utilise generalement les matrices K, 
  31 :     C et M projetees sur la base de modes reels obtenue avec VIBR).
  32 : 
  33 :     Commentaire :
  34 :     _____________
  35 : 
  36 :     MASS1        : matrice de masse
  37 :                    (type RIGIDITE, sous-type MASSE)
  38 : 
  39 :     RIG1         : matrice de rigidite
  40 :                    (type RIGIDITE, sous-type RIGIDITE)
  41 : 
  42 :     AMOR1        : matrice d'amortissement
  43 :                    (type RIGIDITE, sous-type AMORTISSEMENT)
  44 : 
  45 :     Si le type des matrices correspond a MASSE, RIGIDITE et AMORTISSE-
  46 :     -MENT (dans le cas ou cette derniere est specifiee), elles sont 
  47 :     triees et leur ordre d'entree n'a pas d'importance. Sinon elles sont 
  48 :     traitees selon leur ordre d'entree.
  49 : 
  50 : 
  51 :     BAS1         : base de modes reels, sur laquelle les matrices ont 
  52 :                    ete eventuellement projetees (type TABLE, sous-type 
  53 :                    BASEMODA). Sa specification permet la restitution  
  54 :                    d'une base de modes complexes physique.
  55 :                    
  56 :     ENT1         : entier specifiant le nombre de couple de modes 
  57 :                    complexes de plus bas module a sortir. Par defaut,
  58 :                    tous sont fournis en sortie.
  59 : 
  60 :     BAS2         : objet resultat contenant les valeurs et les vecteurs
  61 :                    propres complexes (type TABLE, sous-type BASEMODA).
  62 : 
  63 : REMARQUE : structure de BAS2 (objet de type TABLE)
  64 : _________  
  65 : 
  66 :            
  67 :            BAS2.'SOUSTYPE'    = 'BASE_MODALE'
  68 :            BAS2.'MODES'       = TAB2  
  69 :            BAS2.'CONVERGENCE' = LOGIQ1  
  70 : 
  71 :            TAB2.'SOUSTYPE'    = 'BASE_DE_MODES'
  72 :            TAB2.'MAILLAGE'    = MAIL1
  73 :            TAB2.IMOD          = TAB3 (objet de type TABLE)
  74 : 
  75 :            TAB3.'SOUSTYPE'      = 'MODE_COMPLEXE'
  76 :            TAB3.'POINT_REPERE'  = PT1
  77 :            TAB3.'NUMERO_MODE'   = NUMOD
  78 :            TAB3.'FREQUENCE_REELLE'      = FLOT1
  79 :            TAB3.'FREQUENCE_IMAGINAIRE'  = FLOT2
  80 :            TAB3.'DEFORMEE_MODALE_REELLE'     = CH1
  81 :            TAB3.'DEFORMEE_MODALE_IMAGINAIRE' = CH2
  82 : 
  83 : 
  84 :       MAIL1 : support geometrique des modes (type MAILLAGE)
  85 :       
  86 :       IMOD  : nombre variant de 1 au nombre de modes calcules (ENTIER) 
  87 :       
  88 :       PT1   : point servant a reperer le mode (POINT)
  89 :       
  90 :       NUMOD : numero du mode (ENTIER)
  91 :  
  92 :       FLOT1 : partie reelle de la frequence propre (FLOTTANT)
  93 :       FLOT2 : partie imaginaire de la frequence propre (FLOTTANT)
  94 : 
  95 :       CH1   : partie reelle du vecteur propre (CHPOINT)
  96 :       CH2   : partie imaginaire du vecteur propre (CHPOINT)
  97 :       
  98 :       LOGIQ1 : logique indiquant si VIBC a converge
  99 : 
 100 : 
 101 :     --------------------
 102 :     | 2eme possibilite |
 103 :     --------------------
 104 : 
 105 :     
 106 :       BAS2 = VIBC  RIG1 RIG2 RIG3 RIG4 ;
 107 :       
 108 : 
 109 : 
 110 :     Objet :
 111 :     _______
 112 : 
 113 :     L'operateur VIBC recherche les valeurs propres complexes 
 114 :     de la matrice de monodromie A :
 115 :           [ A - \lambda I ] . X = 0
 116 :     A est fournie sous la forme de 4 RIGIDITEs elementaires :
 117 :           A = [ RIG1  RIG2 ] 
 118 :               [ RIG3  RIG4 ]
 119 :     L'algorithme utilise est le QZ, ce qui implique que la taille du 
 120 :     probleme est censee etre petite.
 121 : 
 122 :     Commentaire :
 123 :     _____________
 124 : 
 125 :     RIG1,2,3 et 4 : matrice de rigidite projetee
 126 :                    (type RIGIDITE, sous-type indifferent)
 127 : 
 128 :     BAS2  : objet TABLE resultat contenant les valeurs propres
 129 :          
 130 : 

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