Télécharger continu.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ CONTINU  NOTICE  BP208322  16/09/01    21:15:02     9010           
   2 :                                              DATE     16/09/01
   3 : 
   4 : 
   5 :     Procedure CONTINU                       Voir aussi : PASAPAS HBM
   6 :     ______________                                       CON_CALC AFT  
   7 :                                                                                 
   8 :     Objet :
   9 :     ______
  10 : 
  11 :     La procedure CONTINU propose de resoudre des problemes non-lineaires
  12 :     poses sous la forme d'equations algebriques (1) qui dependent
  13 :     d'un parametre (pseudo-temps noté t)
  14 :     par une methode de continuation par pseudo-longueur d'arc.
  15 :        
  16 :       R(U,t)  =  Fext(t) + Fnl(U) - Fint(U,\sigma)  = 0           (1.a)
  17 :       R(U,t)  =  Fext(t) + Fnl(U) + Z(t) U          = 0           (1.b)
  18 :     
  19 :     avec :
  20 :     
  21 :       R    : vecteur Residu
  22 :       U    : vecteur des inconnues
  23 :       Fext : vecteur des forces exterieures 
  24 :       Fint : vecteur des forces internes (=\int B^T \sigma)
  25 :       Fnl  : vecteur des forces non-lineaires
  26 :       Z    : matrice de raideur dynamique (voir la procedure HBM)
  27 :       t    : pseudo-temps    
  28 :     
  29 :     
  30 :     Le calcul est realise en 2 etapes :
  31 :     
  32 :     1. Pas predicteur :
  33 :        
  34 :        A partir d'une precedente position (U_n, t_n), le probleme 
  35 :        linearise est resolu (2) et une nouvelle position (U_p, t_p) est 
  36 :        trouvee en imposant la longueur de l'increment (3) egale a ds.
  37 :        
  38 :           dR/dU * dU_p = - dR/dt * dt_p                             (2)
  39 :           ds^2 = dt_p^2/dt_ref^2 + dU_p^T*dU_p / dU_ref^2           (3)
  40 :     
  41 :        avec : 
  42 :           dU_p = U_p - U_n
  43 :           dt_p = t_p - t_n
  44 :           ds   = 1 initialement, mais de valeur adaptative 
  45 :                  selon la difficulte de convergence
  46 : 
  47 :     2. Pas correcteurs :
  48 :        
  49 :        A partir de la position predite (U_p, t_p), des corrections 
  50 :        successives (4) sont realisees dans le plan orthognal a la 
  51 :        prediction jusqu'a  rendre le residu inferieur a une tolerance 
  52 :        donnee.
  53 :        
  54 :           [ dR/dU  dR/dt ] * (dU^(i))  =  (-R^(i-1))                (4)
  55 :           [ dU_p   dt_p  ]   (dt^(i))     ( 0      )
  56 :     
  57 :     
  58 :     
  59 :     Entree :    (on indique entre parentheses les entrees valables 
  60 :     _______      uniquement avec les problemens de type a ou b)
  61 :     
  62 :     TABHBM  =  TABLE
  63 :  
  64 :        . 'HBM'                        = VRAI pour indiquer que l'on 
  65 :                                         souhaite resoudre un probleme
  66 :                                         du type (b)
  67 :                                         (FAUX par defaut)
  68 :        . 'MODELE'                     = modele mecanique utilise (a)
  69 :        . 'CARACTERISTIQUES'           = materiau et caracteristiques (a)
  70 :        . | 'RIGIDITE_CONSTANTE' (a) | = raideur K (hors modele) 
  71 :          | 'RIGIDITE_HBM'       (b) |   
  72 :        . | 'AMORTISSEMENT_CONSTANT' (a)| = amortissement C
  73 :          | 'AMORTISSEMENT_HBM'      (b)|
  74 :        . | 'MASSE_CONSTANTE' (a) |    = masse M
  75 :          | 'MASSE_HBM'       (b) |
  76 :        . | 'BLOCAGES_MECANIQUES' (a) |= Kblocages (hors modele)
  77 :          | 'BLOCAGES_HBM'        (b) |
  78 :        . 'CHARGEMENT'                 = Fext(t)
  79 :        
  80 :        . 'MAXI_DEPLACEMENT'           = dU_ref
  81 :        . 'TEMPS_CALCULES'             = listreel de la discretisation 
  82 :                                         souhaitee pour les pseudo-pas de
  83 :                                         temps t 
  84 :                                         (les temps rellements converges
  85 :                                         seront differents et stockes a  
  86 :                                         l'indice TEMPS_PROG)
  87 :        . 'FREQUENCE' (b)              = evolution de la frequence 
  88 :                                         fondamentale du probleme w(t)
  89 :                                         (w(t)=t par defaut)
  90 :                                         voir remarque 1 pour les unites
  91 :                                         
  92 :        . 'GRANDS_DEPLACEMENTS'        = VRAI en grand deplacements (a)
  93 :        . 'PROCEDURE_CHAR_MECA'        = VRAI si forces non-lineaires
  94 :                                         a: pression suiveuse
  95 :                                         b: terme Fnl(U)
  96 :                                         voir remarque 2
  97 :        . 'PROCEDURE_FREQUENCE_TEMPS' (b) = 'AFT' pour l'utilisation de 
  98 :                                         la procedure AFT lors du calcul
  99 :                                         des efforts et Jacobienne non-
 100 :                                         lineaires
 101 :        . 'N_PT_TFR' (b)               = 2**N_PT_TFR points seront 
 102 :                                         utilises pour la discretisation
 103 :                                         temporelle lors de l'AFT
 104 :                                         
 105 :        . 'PAS_SAUVES'                 = entier N indiquant de sauver les 
 106 :                                         resultats (U,t,\sigma) dans une
 107 :                                         table tous les N pas
 108 :                                         
 109 :        . 'STABILITE'                  = listmots de mots-cles parmi :
 110 :                                       + DIAG pour sauver le nombre de 
 111 :                                         termes diagonaux negatifs 
 112 :                                         comptes lors de la factorisation
 113 :                                       + FLOQ pour le calcul des 
 114 :                                         exposants de Floquet (b)
 115 :        . |'RESULTATS'     (a) |       = table des resultats attendus
 116 :          |'RESULTATS_HBM' (b) |             
 117 :             . i . 'POINT_MESURE'        
 118 :                 . 'COMPOSANTES'     
 119 :                 . 'COULEUR'
 120 :                 . 'TITRE'
 121 :                                         
 122 :        . 'MAXITERATION'               = nombre maxi d'iterations par pas
 123 :                                         (24 par defaut)
 124 :        . 'NB_ITERATION'               = nombre d'iterations juge ideal
 125 :                                         (6 par defaut)
 126 :        . 'MAXIPAS'                    = nombre maxi de pas 
 127 :                                         (1000 par defaut)
 128 :        . 'PRECISION'                  = tolerance relative sur le residu
 129 :                                         (1.E-6 par defaut)
 130 :        . 'COMPOSANTES' 
 131 :             . | 'FORCE'     (a) |     = listmots des composantes a uti-
 132 :             . | 'FORCE_HBM' (b) |       -liser dans le produit scalaire
 133 :                                         definissant la norme du residu
 134 :  
 135 : 
 136 :     Sortie :
 137 :     _______
 138 :     
 139 :     TABHBM 
 140 :        . 'TEMPS'                      = table des pseudo-temps t
 141 :        . 'DEPLACEMENTS'               = table des inconnues    U
 142 :        . 'CONTRAINTES'                = table des contraintes \sigma (a)
 143 :        . 'RESULTATS_STABILITE'        
 144 :             . 'DIAG'                  = listreel des resultats de 
 145 :                                         DIAG(dR/dU)
 146 :             . 'FLOQ' 
 147 :                 . 'EXPOSANT_REEL'    |= table des parties relles et 
 148 :                 . 'EXPOSANT_IMAG'    |  imaginaires des exposants 
 149 :                                         de Floquet (b)
 150 :        . 'TEMPS_PROG'              = listreel des pseudo-temps calcules
 151 :        . | 'RESULTATS_HBM' (a) |
 152 :          | 'RESULTATS'     (b) | 
 153 :             . j . 'RESULTATS'      = listreel des resultats
 154 :                 . 'RESULTATS_EVOL' = evolution associee
 155 :                 
 156 : 
 157 :     Remarque 1 : probleme de type (b)
 158 :     ___________
 159 :      
 160 :     Les matrices ***_HBM peuvent etre generees par la procedure HBM.
 161 :     La frequence par defaut est en rad/s. Pour travailler en Hz, il faut
 162 :     multiplier AMORTISSEMENT_HBM par 2*pi et MASSE_HBM par (2*pi)**2.
 163 :     
 164 :       
 165 :     Remarque 2 :
 166 :     ___________
 167 :     
 168 :     Le calcul des differents termes (Residu, Raideur et leur derivee)
 169 :     est realise par la procedure CON_CALC definie dans le jeu de 
 170 :     donnees dgibi.
 171 :     Si PROCEDURE_CHAR_MECA = VRAI, CON_CALC appelle l'une des 
 172 :     procedures suivantes pour le calcul des termes non-lineaires :
 173 :     - CHARMECA (par defaut)
 174 :     - AFT      (si 'PROCEDURE_FREQUENCE_TEMPS' = mot 'AFT';) qui appelle
 175 :       a son tour CHARMECA qui realise le calcul des termes non-lineaires
 176 :       dans le domaine temporel.
 177 :                                         
 178 : 

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