基于局部平稳法的粘滞阻尼被动控制结构 - 振动工程学报.doc

基于局部平稳法的粘滞阻尼被动控制结构抗震可靠度分析 孙香红1，孔凡 2, *，李书进2 （1. 长安大学 建筑工程学院，陕西 西安，7100612. 武汉理工大学 土木工程与建筑学院，湖北 武汉，430070） 摘要： 式中，为待分析信号的持时；为广义谐和小波中用于表示第阶小波的双指标因子； 为虚数符号；为了计算方便，认为每阶小波频宽相等；为时间平移因子且。将式进行Fourier变换后，可得到广义谐和小波的“盒状”频域表达 非平稳随机过程的局部平稳广义谐和小波模型 最初由Nason等人发展的、基于非抽样小波函数的非平稳向量随机过程模型，即局部平稳小波（Local Stationary Wavelet, LSW）模型可以写为[22] 式中，为非平稳向量随机过程的第个分量；为的维数；分别为尺度和时间平移指标；为相互独立的正交随机序列；为相应尺度和时间处的向量幅值系数；为尺度和时间平移后的非抽样小波基函数。 可以证明，当LSW模型中小波函数为广义谐和小波且各尺度小波频宽较窄时，式进一步可表示为 式中，为待分析信号在时间-频率子域和上的子向量信号；与为维相互独立的0均值随机向量，且有下式成立 其中，表示数学期望，为定义于局部时-频子域上向量随机过程的双边功率谱密度矩阵，其主对角线元素为的局部自谱，非对角元素为局部互谱，即 耗能减震结构的激励-响应功率谱关系 考虑在完全非平稳随机地震动作用下并可简化为剪切型集中质量模型的维多自由度系统 其中，分别为维质量、阻尼与刚度矩阵；为维质量向量；当为维相对位移向量时，各系数矩阵为 一般地，结构阻尼矩阵可用Rayleigh矩阵表示，即 其中，为Rayleigh阻尼系数，可由前二阶振型阻尼比得到。 粘滞液体阻尼器性能试验表明，虽然精确的阻尼器数学模型可由Maxwell模型[7]描述，但对于大多数在地震或强风作用下振动频率较低的工程结构而言，液体粘滞阻尼器提供的附加刚度几乎可以忽略。因此，粘滞阻尼能在工程结构感兴趣的频率范围内提供仅依赖于层间速度的阻尼力，即 其中，为设置于第层间的阻尼器提供的阻尼力；及为零频粘滞阻尼系数和流动指数，由具体器型号确定，文献[23]给出了确定粘滞阻尼系数的方法；当阻尼器水平安装且认为支撑刚度较大时，为第i层间位移。作为研究的初阶并为了突出重点，可不失一般性地认为阻尼器中的液体为流动指数的牛顿流体；当时，亦可通过能量等效原理将非线性阻尼等效为线性阻尼[24]。 因此，设置有线性粘滞阻尼器的被动控制结构，在随机地震动作用下简化多自由度的集中质量模型为 式中，为粘滞阻尼器附加的阻尼矩阵，当阻尼器满布于各层间时，其形式与结构自身刚度矩阵一致，即 当第i层无阻尼器设置时，取。 如将耗能减震系统随机动力位移响应和激励表示为式(14)的形式，同时注意到 式中，为0均值随机变量，且 其中，为完全非平稳随机地震动在局部时间-频率子域上的双边功率谱密度。进一步将式和代入耗能减震系统动力微分方程的分量形式，即 可得， 式中，为耗能减震结构总体阻尼矩阵，由原结构Rayleigh阻尼矩阵和阻尼器附加阻尼矩阵组成。将式两式分别乘以它们的共轭转置并求期望，同时引入式及式可得 将式化为简洁形式可得 式中， 为局部时间-频率子域上的传递函数。 式建立了在局部时-频子域上非平稳地震动功率谱与响应功率谱之间的关系。可见，式与经典线性平稳随机振动分析激励与响应功率谱关系的相似程度。这一点不足为奇：由于对非平稳随机过程在时-频区间上的划分，同时假定在子时间区间上的局部平稳性，直接导致了如式所示的功率谱关系式。基于如下判断，可认为局部平稳假定是合理的：对于大多数在随机地震动激励下的工程结构而言，由于激励相对于时间尺度的慢变性，从而使结构响应相对于此时间尺度也可认为是慢变的；更重要的是，尤其对于附加有一定阻尼装置的被动控制结构而言，阻尼器提供的等效临界阻尼比往往比原工程结构更大，从而使大阻尼比受控结构的瞬态响应可较快地耗散。 数值算例 为了验证局部平稳法在多自由度结构耗能减震中的应用效果，以一个算例进行分析。某六层框架结构，简化后可采用集中质量多自由度模型，其中各层集中质量向为、层间刚度系数以及相应模态频率如表1所示，粘滞阻尼器的线性粘滞阻尼系数取为N s/m。结构系统阻尼采用Rayleigh阻尼且前两阶模态阻尼比分别取为。在其中两层对称设置4个阻尼器：不失一般性，可考虑阻尼器水平放置且忽略支撑刚度影响。为了考察不同阻尼器数目及分布情况对受控结构的随机动力响应概率特征的影响，考虑以下四种阻尼器位置分布情况： (a) [1,1,1,1,1,1]； (b) [1,1,0,0,0,0]； (c) [0,0,0,0,1,1]；