Notices | Cast3M

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