Accueil | Cast3M

   1 : $$$$ DYNAMIC  NOTICE  MB234859  17/10/02    21:15:09     9577           
   2 :                                              DATE     17/10/02
   3 :                                              
   4 :   Procedure DYNAMIC                       Voir aussi : DYNE PASAPAS
   5 :     -----------------  
   6 :     TAB2 = DYNAMIC TAB1 ;
   7 : 
   8 : 
   9 :     Objet :
  10 :     _______
  11 : 
  12 :     Cette procedure permet d'effectuer un calcul dynamique pas a pas.
  13 :     
  14 :     Elle peut utiliser les algorithmes :
  15 :     - Newmark centre (schema de l'acceleration moyenne) (par defaut)
  16 :     - HHT
  17 :     - alpha-generalise
  18 :     
  19 :     Les arguments d'entree et de sortie de la procedure sont des tables
  20 :     definies ci-apres.
  21 :     
  22 : 
  23 :     Commentaire :
  24 :     _____________
  25 : 
  26 :     
  27 :     +----------------------------------------------------------+
  28 :     | TAB1 : Objet de type TABLE contenant les donnes d'ENTREE |
  29 :     +----------------------------------------------------------+
  30 : 
  31 :       DONNEES DU PROBLEME
  32 :       ###################
  33 :     
  34 :       TAB1 . 'DEPL'   : deplacement initial     [CHPOINT] 
  35 :       TAB1 . 'VITE'   : vitesse initiale        [CHPOINT] 
  36 :       TAB1 . 'CHAR'   : chargement              [CHARGEME]
  37 :       TAB1 . 'RIGI'   : raideur                 [RIGIDITE]
  38 :       TAB1 . 'MASS'   : masse                   [RIGIDITE]
  39 :     ( TAB1 . 'AMOR' ) : amortissement           [RIGIDITE]
  40 : 
  41 :     
  42 :       TEMPS DE CALCUL
  43 :       ###############
  44 :     
  45 :       TAB1 . 'TEMPS_CALCULES' : [LISTREEL] strictement croissant 
  46 :                                 contenant tous les instants de 
  47 :                                 calcul (dont l'instant initial)
  48 :                         
  49 :                         
  50 :       PARAMETRES DE SAUVEGARDE
  51 :       ########################
  52 : 
  53 :       SORTIE GIBIANE
  54 :       --------------
  55 :       
  56 :     | TAB1 . 'PAS_SAUVES'       : [ENTIER] indiquant la cadence de la
  57 :     |                             sauvegarde ou [MOT] valant 'TOUS' ou
  58 :  OU |                             'FINAL' (par defaut, on conserve 1 
  59 :     |                             pas de temps sur 4)
  60 :     | TAB1 . 'TEMPS_SAUVES'     : [LISTREEL] des instants a sauver
  61 :     
  62 :     ( TAB1 . 'MAILLAGE_SAUVE' ) : [MAILLAGE] pouvant etre fourni pour
  63 :                                   limiter le support geometrique des
  64 :                                   donnees sauvegardees
  65 : 
  66 :     ( TAB1 . 'SAUV' ) : VRAI si l'on souhaite SAUVER a chaque pas
  67 :                         de maniere incrementale dans le fichier
  68 :                         a definir prealablement par la commande : 
  69 :                         "OPTI 'SAUV' nomfic ;"
  70 :                         
  71 :     ( TAB1 . 'ECON' ) : VRAI si l'on souhaite FANTomiser les resultats 
  72 :                         sauvegardes sur disque
  73 :                
  74 : 
  75 :       SORTIE VTK (PARAVIEW) => OPTIONNELLE
  76 :       ---------------------
  77 :       
  78 :     | TAB1 . 'PAS_SAUVES_VTK'   : idem ci-dessus
  79 :  OU | TAB1 . 'TEMPS_SAUVES_VTK'   (par defaut : aucune sortie VTK) 
  80 : 
  81 :       TAB1 . 'MAILLAGE_VTK'     : objet [MAILLAGE] ou [TABLE] passe a
  82 :                                   l'operateur SORT 'VTK' definissant la
  83 :                                   ou les geometries a sauvegarder
  84 :                                   (seulement si l'un des 2 indices
  85 :                                    ci-dessus est present)
  86 :                                   
  87 :     ( TAB1 . 'FICHIER_VTK' )    : [MOT] indiquant l'emplacement ou
  88 :                                   seront ecrits les fichiers VTK
  89 :                                   Beaucoup de fichiers peuvent etre
  90 :                                   crees => on recommande de les placer
  91 :                                   dans un sous-repertoire (par defaut,
  92 :                                   tout est dans le dossier courant)
  93 : 
  94 :                                   
  95 :       SORTIE CSV (TABLEUR) => OPTIONNELLE
  96 :       --------------------
  97 :                    
  98 :     | TAB1 . 'PAS_SAUVES_CSV'   : idem ci-dessus
  99 :  OU | TAB1 . 'TEMPS_SAUVES_CSV'   (par defaut : aucune sortie CSV) 
 100 :  
 101 :       TAB1 . 'MAILLAGE_CSV'     : objet [MAILLAGE] contenant les noeuds
 102 :                                   ou seront enregistrees les donnees
 103 :                                   (seulement si l'un des 2 indices
 104 :                                    ci-dessus est present)
 105 : 
 106 :     ( TAB1 . 'COMPOSANTES_CSV' ): [LISTMOTS] des composantes a conserver
 107 :                                   (par defaut, on conserve toutes les
 108 :                                    composantes disponibles)
 109 :                                   
 110 :     ( TAB1 . 'FICHIER_CSV' )    : [MOT] indiquant l'emplacement ou
 111 :                                   sera cree le fichier CSV (par defaut,
 112 :                                   on n'ecrit rien sur le disque)
 113 :  
 114 : 
 115 :       APPEL A DES PROCEDURES SPECIFIQUES
 116 :       ##################################
 117 :     
 118 :     ( TAB1 . 'CHARMECA' ) : VRAI si appelle la procedure CHARMECA
 119 :                             definie par ailleurs par l'utilisateur et
 120 :                             qui renvoie une table TCHAR avec un ou
 121 :                             plusieurs indices parmi :
 122 :                             
 123 :                               'ADDI_MATRICE' => ajout a l'operateur
 124 :                               'ADDI_SECOND'  => ajout au second membre
 125 :                                                 (composantes FLX)
 126 :                               'ADDI_KNL'     => ajout a la raideur
 127 :                               'ADDI_CNL'     => ajout a l'amortissement
 128 :                               'ADDI_FNL'     => ajout au second membre
 129 :                                                 (composantes hors FLX)
 130 :                                               
 131 :     ( TAB1 . 'VITEUNIL' ) : correction des vitesses lors d'un impact 
 132 :                             (VRAI par defaut, seulement possible avec
 133 :                              le schema de Newmark acceleration moyenne)
 134 :  
 135 : 
 136 :       SCHEMA D'INTEGRATION TEMPORELLE => OPTIONNEL
 137 :       ###############################
 138 :     
 139 :     | TAB1 . 'ALPHA_F'  : active le schema HHT et donne son parametre
 140 :     |                     \alpha
 141 :  OU |                        
 142 :     | TAB1 . 'RHO_INF'  : active le schema alpha-generalise et donne
 143 :     |                     son parametre \rho^{\intfy} (rayon spectral
 144 :     |                     a l'infini)
 145 : 
 146 :     
 147 :                        
 148 :     +-----------------------------------------------------------------+
 149 :     | TAB2 : Objet de type TABLE (aussi stocke dans TAB1.'RESULTATS') |
 150 :     |        contenant les objets GIBIANE sauvegardes a chaque pas de |
 151 :     |        temps demande                                            |
 152 :     +-----------------------------------------------------------------+
 153 : 
 154 :       TAB2 . I . 'TEMP' : temps        [FLOTTANT]
 155 :       TAB2 . I . 'DEPL' : deplacement  [CHPOINT]
 156 :       TAB2 . I . 'VITE' : vitesse      [CHPOINT]
 157 : 
 158 :     
 159 :                        
 160 :     +-----------------------------------------------------------------+
 161 :     | Un objet de type TABLE est stocke dans TAB1.'RESULTATS_CSV'     |
 162 :     | contenant autant de LISTREEL que de grandeurs demandees dans la |
 163 :     | sortie CSV (ainsi qu'un LISTREEL pour le temps)                 |
 164 :     +-----------------------------------------------------------------+
 165 :        
 166 :     
 167 :