Accueil | Cast3M

   1 : $$$$ IMPE     NOTICE  MB234859  17/10/02    21:15:14     9577           
   2 :                                              DATE     17/10/02
   3 : 
   4 :      Operateur IMPE :     Voir aussi: MASS, RIGI, AMOR , CORIOLIS, GYRO
   5 :      ---------------                                  
   6 : 
   7 :  
   8 :       RIG1 = IMPE RIG2 (RIG3) (REEL1) (| 'MASSE'            | )  (FLAM); 
   9 :                                        | 'RAIDEUR'          |  
  10 :                                        | 'AMORTISSEMENT'    |  
  11 :                                        | 'QUELCONQUE'  LISREEL1  |  
  12 :   
  13 : 
  14 : Objet:
  15 : ------ 
  16 :   
  17 : 
  18 :  
  19 :   L'operateur IMPEDANCE calcule les matrices d impedance 
  20 :   de sous-type MASSE, RAIDEUR, ou AMORTISSEMENT  
  21 :                                                             
  22 : 
  23 :   INPUT :
  24 :   -------
  25 :  
  26 :   RIG2: Matrice de rigidite initiale dont les inconnues primales et 
  27 :         duales sont dans le domaine temporel (UX, UY, ..., FX, FY, ...)
  28 :       
  29 :   REEL1: REEL donnant la pulsation (en rad/s) pour laquelle est calculee
  30 :          l impedance (valeur par defaut: 1)
  31 :       
  32 :   FLAM: LOGIQUE (FAUX par defaut) specifiant si la matrice doit etre de
  33 :         sous type MASSE (necessaire si la matrice est utilise pour 
  34 :         un calcul de flambage) 
  35 : 
  36 :   'MASSE', 'RAIDEUR', 'AMORTISSEMENT', 'QUELCONQUE': Mot cle specifiant
  37 :   le type de matrice d impedance a creer
  38 : 
  39 :   L'option 'QUELCONQUE' necessite une LISREEL1 comprenant 4 reels:
  40 :   LISREEL1 = PROG REELA REELB REELC REELD; 
  41 :        
  42 :   Il est également possible de fournir 2 matrices en entree afin de 
  43 :   coupler les ddls de meme nom entre ces 2 matrices. 
  44 :   C'est par exemple le cas lorsqu'on cherche à coupler les déplacements
  45 :   symetriques et antisymétriques lors d'un calcul en mode de Fourier
  46 :   (cas des machines tournantes par ex.) avec :
  47 :   - RIG2 qui est la matrice calculee avec les ddls symetriques 
  48 :     (c'est a dire apres l'instruction OPTION MODE FOUR nhar)
  49 :   - RIG3 qui est la matrice calculee avec les ddls antisymetriques 
  50 :     (c'est a dire apres l'instruction OPTION MODE FOUR (-1*nhar))
  51 :  
  52 :   Seule l'option RAIDEUR est permise lorsqu'une matrice comprenant 
  53 :   des relations (i.e. faisant intervenir des multiplicateurs de Lagrange
  54 :   LX ) est fournie en entree.
  55 :        
  56 : 
  57 :   OUTPUT :
  58 :   --------
  59 : 
  60 :   RIG1: Objet RIGIDITE représentant la matrice d impedance 
  61 :         (de meme sous type que la RIG2 ou de sous-type MASSE si le 
  62 :         mot-clé 'FLAM' est précisé)
  63 :         
  64 :   La matrice RIG1 est creee par duplication de la matrice RIG1 sur
  65 :   les degres de liberte imaginaires (IUX, IUY, ..., IFX, IFY, ...)
  66 :   selon les règles suivantes :
  67 : 
  68 :   Cas 'RAIDEUR':     [  RIG2    0   ]
  69 :                      [   0     RIG2 ]
  70 : 
  71 :   Cas 'MASSE':       [ -REEL1^2*RIG2         0     ]
  72 :                      [       0       -REEL1^2*RIG2 ]
  73 : 
  74 :   Cas 'AMORTISSEMENT':   [      0     -REEL1*RIG2 ]
  75 :                          [ REEL1*RIG2      0      ]
  76 : 
  77 :   Cas 'QUELCONQUE':      [ REELA*RIG2     REELB*RIG2 ]
  78 :                          [ REELC*RIG2     REELD*RIG2 ]
  79 : 
  80 : 
  81 :   Dans le cas ou deux matrices sont fournies en entrée, RIG1 vaut :
  82 :   
  83 :   Cas 'RAIDEUR':     [  RIG2    0   ]
  84 :                      [   0     RIG3 ]
  85 : 
  86 :   Cas 'MASSE':       [ -REEL1^2*RIG2         0     ]
  87 :                      [       0       -REEL1^2*RIG3 ]
  88 : 
  89 :   Cas 'AMORTISSEMENT':   [      0     -REEL1*RIG2bar ]
  90 :                          [ REEL1*RIG3bar      0      ]
  91 :                          
  92 :   avec RIG2bar_{ij} = | - RIG2_{ij}  si j refere au ddl UT
  93 :                       |   RIG2_{ij}  pour tous les autres ddls.
  94 : 
  95 :                       
  96 :