$$$$ DECONV3D NOTICE MB234859 17/10/02 21:15:07 9577 DATE 17/10/02 procedure DECONV3D Voir aussi : DECONV ------------------ TABRESU = DECONV3D COUCHE PC FOND_SOL MOD_SOL MAT_SOL DIR GAMMAO F1 F2 (FC) (P_GAMMA) Objet : ------- Cette procedure permet d'effectuer des calculs sismiques d'interaction sol-structure (ISS) en 3D par la methode des elements finis. Elle a deux fonctions lorsque l'on n'évoque pas la methode de reduction de domaine de Bielak [1]: - Formuler la matrice d'amortissement correspondant a la frontiere absorbante visqueuse sur la bordure du maillage de sol, - Calculer le chargement sismique au cours du temps a appliquer sur la frontiere par la deconvolution du mouvement sismique (accelerogramme) donnee en surface libre du sol Lorsque la methode de Bielak[1] est évoquée (nouvelle option, si COUCHE.'BLK' existe), le maillage du sol d'entrée (appelé zone interne) peut être largement reduit. La procedure crée automatiquement une couche d'elements enveloppe (appelée zone intermediare) et calcule le chargement sismique à appliquer sur cette zone. Dans ce cas, l'utilisateur doit créer lui-même une zone de sol externe a l'exterieure de la zone intermediare et mettre sur sa bordure exterieure une frontiere absorbante à l'aide de l'operateur AMOR. Dans cette option, il est possible d'introduire des ondes SH inclinées comme chargement sismique. La resolution du probleme peut s'effectuer ensuite dans le domaine temporel sur l'ensemble du maillage sol-structure avec l'une des procedures d'integration suivantes : - procedure DYNAMIC en cas de comportement lineaire, possibilite d'inclure des liaisons unilaterales pour modeliser le decollement et le glissement du radier, - procedure PASAPAS en cas de comportement non lineaire du sol et de la structure. Commentaire : ------------- En entree : ----------- COUCHE 'TABLE' : table a double indice COUCHE.'BLK' : si existe, deconvolution pour la methode de Bielak COUCHE.'BLK'.'EP_H' : epaisseur horizontale de la zone intermediare verticale COUCHE.'BLK'.'EP_V' : epaisseur veriticle de la zone intermediare horizontale COUCHE.I.'indice' : donnees pour la ieme couche du sol indice 'FRONTIERE' 'MAILLAGE' : frontiere verticale de la ieme couche indice 'MASSE_VOLUMIQUE' 'FLOTTANT' : masse volumique indice 'POISSON' 'FLOTTANT' : coefficient de Poisson indice 'YOUNG' 'FLOTTANT' : module d'Young indice 'AMORTISSEMENT 'FLOTTANT' : amortissement reduit PC 'MAILLAGE' : point de reference situe au milieu du millage de sol ou sur l'axe de symetrie si calcul sur la moitie ou le quart du systeme FOND_SOL 'MAILLAGE' : frontiere horizontale inferieure du sol MOD_SOL 'MMODEL' : modele du sol MAT_SOL 'MCHAML' : materiau du sol DIR 'MOT' : direction de l'acceleration GAMMAO 'UX' pour la direction X (onde SV) 'UY' pour la direction Y (onde SH) 'UZ' pour la direction Z (onde P) GAMMAO 'EVOLUTIO' : acceleration du 'champ libre' en surface du sol comportant une plage initiale a zero acceleration sur au moins 100 pas de temps F1, F2 'FLOTTANT' : frequences sur lesquelles l'amortissement : reduit est ajuste suivant le modele de RAYLEIGH FC 'FLOTTANT' : frequence de coupure pour la deconvolution, par defaut FC = 50 Hz P_GAMMA 'TABLE' : points sur la frontiere verticale pour lesquels on desire sortir les accelerogrammes indice 'ENTREE' :'TABLE' : description de la nature de l'accelerogramme GAMMAO autre que sur la surface libre (facultatif) ou description de l'onde incidente inclinée SH. sous-indice 'NATURE' : nature du point de controle : 'MOT' INSIDE : dans le sol 'MOT' OUTCROP : outcrop du bedrock 'MOT' SH : Onde incidente SH inclinée dans le cas de la methode Bielak, DIR = UX : onde SH dans le plan YZ, DIR = UY : onde SH dans le plan XZ. sous-indice 'ANGLE' : Angle d'incidence de l'onde SH si 'NATURE' = SH sous-indice 'CONTROLE' :'MAILLAGE': si 'NATURE' = 'INSIDE', point de controle sur la frontiere verticale, si 'NATURE' = 'SH', point de controle à la surface du sol pour lequel on impose l'accélérogramme du champ libre GAMMAO. indice 'I' :'MAILLAGE': ieme (i = 1, 2, 3,...) points sur la frontiere verticale pour lesquels on desire sortir l'accelerogramme en champ libre (resultats de deconvolution) En sortie : ----------- TABRESU 'TABLE' : table qui contient les resultats du calcul indice 'CHAR' 'CHARGEMENT' : excitation sismique sur la frontiere du sol indice 'AMOR' 'RIGIDITE' : frontiere absorbante (absent, si Bielak) indice 'DEFO' 'EVOLUTIO' : deformation maximale du sol en fonction de la profondeur indice 'ACCE' 'TABLE' : table qui contient les accelerations des points definis dans P_GAMMA indice 'PAS' 'FLOTTANT' indice 'FCDYN' 'FLOTTANT' : pas de temps et frequence de coupure a utiliser pour le calcul dynamique de l'interaction sol-structure a l'aide de la procedure DYNAMIC ou PASAPAS. indice 'BLK_MAI', 'BLK_MOD', 'BLK_MAT', 'BLK_MAS', 'BLK_RIG', 'BLK_AMO' : maillage, modele, materiau, masse, rigidite et amortissement de la zone intermediaire cree par la procedure Remarques : ----------- 1) La procedure n'accepte pas les frontieres obliques. 2) Le calcul ISS peut etre effectue sur la totalite, la moitie ou le quart du systeme selon la symetrie du maillage. Il faut donner a la procedure le maillage de frontiere correspondant. 3) La methode de Bielak peut permettre de reduire le nombre d'elements du modele et le temps de calcul. 4) Une onde SH incidente a l'angle azimutal quelconque peut être decomposée en deux onde SH respectivement dans les plans YZ et XZ. Reference: ---------- [1] Bielak et al "Domain reduction method for three dimensional earthquake modeling in localized regions" Bulletin of the Seismilogical Society of America, Vol. 93, No. 2, pp817, April 2003
© Cast3M 2003 - Tous droits réservés.
Mentions légales