等效黏性边界单元及黏性边界波动输入方法 Equivalent viscous boundary elements and the method of wave input for viscous boundary.pdfVIP

水 利 学 报 2008年lO月 SHUILIXUEBAO 第39卷第10期 文章编号：0559．9350(2008)1m1248—08 等效黏性边界单元及黏性边界波动输入方法 张顺福，王海波，杨会臣 (中国水利水电科学研究院工程抗震研究中心，北京100044) 摘要：为更好的利用大型通用软件进行结构一地基动力相互作用分析，本文在黏性边界的基础上提出黏性边界单 元的概念，通过构造相同阻尼矩阵的方式提出了能够实现黏性边界功效的特殊有限单元一等效黏性边界单元。等 效黏性边界单元能够较好的模拟地基的辐射阻尼效应，使黏性边界的实现更加方便。同时，为简化大型通用软件 中地震波的输入，提出了适合于黏性边界条件下水平成层地基中垂直入射地震波的输入方法。该方法不需要借助 其它辅助计算就可直接在各种大型通用有限元软件中用一个数值模型实现地震波动输入，具有很强的通用性。此 输入方法极大的减少了地震波动输入的数据准备工作，降低了操作的复杂程度。数值算例表明。等效黏性边界单 元及文中提出的地震波动输入是简单可靠的。 关键词：黏性边界；等效黏性边界单元；结构一地基动力相互作用；波动输入 中图分类号：TU973 文献标识码：A 1研究背景 结构一地基的动力相互作用是影响结构地震动响应的一个重要因素。在地震作用下，结构一地基 动力相互作用主要包含2个内容：地基辐射阻尼效应的模拟和地震荷载的输入。局部人工边界由于其 时空解耦性，在地基辐射阻尼效应模拟中获得广泛应用。1969年，rysmer等…基于一维波动理论提出黏 性边界。在人工边界上设置阻尼器吸收广义结构传向地基的波动能量。黏性边界的阻尼系数只与介质 的材料性质有关，而与波的频率无关，可在时域或频域内求解。当波以小角度外行时，黏性边界能够很 好的将能量透射出去。廖振鹏等心’31基于透射理论提出透射边界。Deeks等H1、刘晶波等b石1基于柱面波 动方程建立了能模拟半无限介质弹性恢复能力的二维黏弹性人工边界。刘晶波等1基于球面波动方程 将黏弹性人工边界推广到三维情形。基于矩阵等效原理，刘晶波等¨1提出采用普通有限单元模拟的二 维黏弹性边界单元，这种单元使黏弹性边界实施更为方便，是实施人工边界的一种新思路。与各种人工 边界条件相应的地震波动输入方法研究成果很多，但大多没有针对大型通用软件提出一种通用的、实施 简便的方法，这无疑增加了人工边界在结构地震动响应分析时的应用难度。 对于实际工程计算，黏性边界和黏弹性边界均有较好的吸能效果，两者所得计算结果差异很小旧1。 综合考虑计算精度和效率，本文的研究针对黏性边界进行。在黏性边界的基础上提出黏性边界单元，通 过构造相同阻尼矩阵的方式提出等效黏性边界单元。同时，提出一种黏性边界条件下水平成层地基中 地震波垂直入射时的输入方法，用一个数值模型实现了大型通用软件中地震波动的直接输入。从工程 实际的角度出发，文中讨论仅限于水平成层地基中地震波垂直入射的情况，并假定振动方向沿边界法向 的波为压缩波，沿切向的波为剪切波。 收稿日期：2007．10--10 作者简介：张顺福(1982一)，男。广西桂林人，硕士，主要从事结构抗震研究。E—mail：shuMu9227@sohu．corn —1248— 万方数据 2等效黏性边界单元及黏性边界波动输入方法 2．1黏性边界理论取介质1和介质2交界的法线方向为z轴，z轴正向从介质2指向介质1，介质交 加”(t，z)，则介质交界处的正应力满足¨刨 盯：l。：o=P】c，l[曲(t)一曲7(t)]=P2c，2t{J”(t) (1) 式中：n、c耐(i=1，2)分别为两种介质的密度以及介质中压缩波的传播速度。 由位移连续条件推得两种介质在交界面上的速度相等，从式(1)可以看出，交界面任一点上介质对 另一种介质的分布作用力与该点的运动速度成正比。利用阻尼器的力与速度成正比的性质，可在边界 法向施加一端与交界面连接，另一端固定的阻尼系数为P：Cp：的分布阻尼器代替介质2的作用，这就是 黏性边界的原理。黏性阻尼边界满足了边界处的应力和位移条件，实现了地震波的透射。对于垂直入 射的剪切波，则在边界切向施加阻尼系数为P：c。：的分布阻尼器，c．：为介质2中剪切