Accueil | Cast3M

   1 : $$$$ DECONV3D NOTICE  MB234859  17/10/02    21:15:07     9577           
   2 :                                              DATE     17/10/02
   3 : 
   4 :   procedure DECONV3D                          Voir aussi : DECONV
   5 :     ------------------  
   6 :    TABRESU =  DECONV3D COUCHE PC FOND_SOL MOD_SOL MAT_SOL DIR
   7 :                        GAMMAO F1 F2 (FC) (P_GAMMA)
   8 : 
   9 : 
  10 : 
  11 : Objet :
  12 : -------
  13 : 
  14 : Cette procedure permet d'effectuer des calculs sismiques d'interaction
  15 : sol-structure (ISS) en 3D par la methode des elements finis. Elle a
  16 : deux fonctions lorsque l'on n'évoque pas la methode de reduction de 
  17 : domaine de Bielak [1]:
  18 :   - Formuler la matrice d'amortissement correspondant a la frontiere
  19 :     absorbante visqueuse sur la bordure du maillage de sol,
  20 :   - Calculer le chargement sismique au cours du temps a appliquer sur
  21 :     la frontiere par la deconvolution du mouvement sismique
  22 :     (accelerogramme) donnee en surface libre du sol
  23 : 
  24 : Lorsque la methode de Bielak[1] est évoquée (nouvelle option, si 
  25 : COUCHE.'BLK' existe), le maillage du sol d'entrée (appelé zone interne)
  26 : peut être largement reduit. La procedure crée automatiquement une 
  27 : couche d'elements enveloppe (appelée zone intermediare) et calcule le
  28 : chargement sismique à appliquer sur cette zone. Dans ce cas, 
  29 : l'utilisateur doit créer lui-même une zone de sol externe a l'exterieure
  30 : de la zone intermediare et mettre sur sa bordure exterieure une 
  31 : frontiere absorbante à l'aide de l'operateur AMOR. Dans cette option, il 
  32 : est possible d'introduire des ondes SH inclinées comme chargement 
  33 : sismique.
  34 : 
  35 : La resolution du probleme peut s'effectuer ensuite dans le domaine
  36 : temporel sur l'ensemble du maillage sol-structure avec l'une des 
  37 : procedures d'integration suivantes :
  38 :   - procedure DYNAMIC en cas de comportement lineaire, possibilite
  39 :     d'inclure des liaisons unilaterales pour modeliser le decollement
  40 :     et le glissement du radier,
  41 :   - procedure PASAPAS en cas de comportement non lineaire du sol et de
  42 :     la structure.
  43 : 
  44 : Commentaire :
  45 : -------------
  46 : 
  47 : En entree :
  48 : -----------
  49 : 
  50 : COUCHE 'TABLE' : table a double indice
  51 :   COUCHE.'BLK' : si existe, deconvolution pour la methode de Bielak
  52 :   COUCHE.'BLK'.'EP_H' : epaisseur horizontale de la zone intermediare
  53 :                         verticale
  54 :   COUCHE.'BLK'.'EP_V' : epaisseur veriticle de la zone intermediare
  55 :                         horizontale
  56 :   COUCHE.I.'indice' : donnees pour la ieme couche du sol
  57 : 
  58 :     indice 'FRONTIERE'       'MAILLAGE' : frontiere verticale de la
  59 :                                           ieme couche
  60 :     indice 'MASSE_VOLUMIQUE' 'FLOTTANT' : masse volumique
  61 :     indice 'POISSON'         'FLOTTANT' : coefficient de Poisson
  62 :     indice 'YOUNG'           'FLOTTANT' : module d'Young
  63 :     indice 'AMORTISSEMENT    'FLOTTANT' : amortissement reduit
  64 : 
  65 : PC       'MAILLAGE' : point de reference situe au milieu du millage
  66 :                       de sol ou sur l'axe de symetrie si calcul sur
  67 :                       la moitie ou le quart du systeme
  68 : FOND_SOL 'MAILLAGE' : frontiere horizontale inferieure du sol
  69 : MOD_SOL  'MMODEL'   : modele du sol
  70 : MAT_SOL  'MCHAML'   : materiau du sol
  71 : DIR      'MOT'      : direction de l'acceleration GAMMAO
  72 :                       'UX' pour la direction X (onde SV)
  73 :                       'UY' pour la direction Y (onde SH)
  74 :                       'UZ' pour la direction Z (onde P)
  75 : GAMMAO   'EVOLUTIO' : acceleration du 'champ libre' en surface du sol
  76 :                       comportant une plage initiale a zero
  77 :                       acceleration sur au moins 100 pas de temps
  78 : F1, F2   'FLOTTANT' : frequences sur lesquelles l'amortissement
  79 :                     : reduit est ajuste suivant le modele de
  80 :                       RAYLEIGH
  81 : FC       'FLOTTANT' : frequence de coupure pour la deconvolution,
  82 :                       par defaut  FC = 50 Hz
  83 : P_GAMMA  'TABLE'    : points sur la frontiere verticale pour
  84 :                       lesquels on desire sortir les accelerogrammes
  85 : 
  86 :   indice 'ENTREE'   :'TABLE' : description de la nature de 
  87 :                      l'accelerogramme GAMMAO autre que sur la surface
  88 :                      libre (facultatif) ou description de l'onde 
  89 :                      incidente inclinée SH.
  90 :     sous-indice 'NATURE'   : nature du point de controle :
  91 :                             'MOT' INSIDE : dans le sol
  92 :                             'MOT' OUTCROP : outcrop du bedrock
  93 :                             'MOT' SH : Onde incidente SH inclinée dans 
  94 :                                   le cas de la methode Bielak, 
  95 :                                   DIR = UX : onde SH dans le plan YZ, 
  96 :                                   DIR = UY : onde SH dans le plan XZ.
  97 :     sous-indice 'ANGLE' : Angle d'incidence de l'onde SH si 
  98 :                          'NATURE' = SH
  99 : 
 100 :     sous-indice 'CONTROLE' :'MAILLAGE': 
 101 :                              si 'NATURE' = 'INSIDE', point de controle 
 102 :                              sur la frontiere verticale, 
 103 :                              si 'NATURE' = 'SH', point de controle à la
 104 :                              surface du sol pour lequel on impose 
 105 :                              l'accélérogramme du champ libre GAMMAO.
 106 : 
 107 :   indice 'I'        :'MAILLAGE': ieme (i = 1, 2, 3,...) points sur la
 108 :                       frontiere verticale pour lesquels on desire sortir
 109 :                       l'accelerogramme en champ libre (resultats de
 110 :                       deconvolution)
 111 : 
 112 : En sortie :
 113 : -----------
 114 : TABRESU 'TABLE' :  table qui contient les resultats du calcul
 115 : 
 116 :   indice 'CHAR'  'CHARGEMENT' : excitation sismique sur la frontiere
 117 :                                 du sol
 118 :   indice 'AMOR'  'RIGIDITE'   : frontiere absorbante (absent, si Bielak)
 119 :   indice 'DEFO'  'EVOLUTIO'   : deformation maximale du sol en
 120 :                                 fonction de la profondeur
 121 :   indice 'ACCE'  'TABLE'      : table qui contient les accelerations
 122 :                                 des points definis dans P_GAMMA
 123 :   indice 'PAS'   'FLOTTANT'
 124 :   indice 'FCDYN' 'FLOTTANT'   : pas de temps et frequence de coupure
 125 :                                 a utiliser pour le calcul dynamique
 126 :                                 de l'interaction sol-structure
 127 :                                 a l'aide de la procedure DYNAMIC ou
 128 :                                 PASAPAS.
 129 :   indice 'BLK_MAI', 'BLK_MOD', 'BLK_MAT', 'BLK_MAS', 'BLK_RIG', 
 130 :          'BLK_AMO' : maillage, modele, materiau, masse, rigidite et 
 131 :                      amortissement de la zone intermediaire cree par 
 132 :                      la procedure
 133 : 
 134 : Remarques :
 135 : -----------
 136 : 
 137 : 1) La procedure n'accepte pas les frontieres obliques.
 138 : 
 139 : 2) Le calcul ISS peut etre effectue sur la totalite, la moitie ou le
 140 : quart du systeme selon la symetrie du maillage. Il faut donner a la
 141 : procedure le maillage de frontiere correspondant.
 142 : 
 143 : 3) La methode de Bielak peut permettre de reduire le nombre d'elements 
 144 : du modele et le temps de calcul.
 145 : 
 146 : 4) Une onde SH incidente a l'angle azimutal quelconque peut être 
 147 : decomposée en deux onde SH respectivement dans les plans YZ et XZ.
 148 : 
 149 : Reference:
 150 : ----------
 151 : 
 152 : [1] Bielak et al "Domain reduction method for three dimensional 
 153 :     earthquake modeling in localized regions" Bulletin of the 
 154 :     Seismilogical Society of America, Vol. 93, No. 2, pp817, 
 155 :     April 2003
 156 :