随机结构TMD系统风振响应概率密度演化分析.docVIP

随机结构TMD系统风振响应概率密度演化分析 　　摘要：本文采用概率密度演化方法，对比了随机结构在无控和TMD控制时结构动力响应和概率密度。结果表明，设置TMD后结构的安全性与舒适性大大改善，TMD对结构参数随机性具有一定的容错性。概率密度演化方法计算量小，精度高，与其他方法相比可获得更为精细的结构反应概率信息。 　　关键词：调谐阻尼器；概率密度演化法；风振响应；随机结构 　　0引言 　　工程结构在许多方面都包含着若干不确定因素，如材料参数、边界条件、阻尼与刚度等均存在随机性。这些随机因素会对结构动力响应分析结果带来不容忽略的影响，因此考虑结构参数的随机性非常重要。调谐质量阻尼器的概念是1909年由Ormondroyd和Den Hartog首先提出（Ormondroyd J，et al，1928），最早于1977年应用于美国244m高的Hancock大厦。调谐质量阻尼器一般由附加惯性质量，弹簧和阻尼元件组成，因其构造简单，易于安装，维护方便，经济适用，且无需外部能源，有着其他减震技术无可比拟的优点，广泛受到科研工作者和工程界的青睐（Rohman A，1984；王肇民，1994；李春祥，1999）。 　　目前对TMD的研究多为TMD的优化设计以及工程应用，主要关注确定性结构的动力响应对比；即便考虑随机性，也大都限于激励的随机性，对结构及减震装置的随机性分析不多。李杰提出的随机结构分析的概率密度演化方法可给出结构概率密度的时变过程，有助于深入了解结构的随机性（李杰等，2003；J.Li，et al，2004）。目前有关使用概率密度演化方法对结构动力响应研究的文献尚未见诸于随机结构TMD系统。本文采用线性滤波法中的自回归模型模拟作用在结构上的脉动风激励时程，将概率密度演化方法应用于随机结构―TMD系统，并考虑结构刚度与阻尼的随机性，以及TMD对结构概率密度的影响规律。 　　1 运动方程 　　?}动风荷载作用在结构上引起的振动，可视为是多点平稳随机激励下结构的随机振动问题。当结构在任意第j层安装一套TMD装置，假定TMD质量为md、刚度为kd、阻尼为Cd。结构在风荷载作用下的运动微分方程为： （1） 　　式中 为TMD对结构的作用力，表达式为： 　　（2） 　　TMD的运动方程为： （3） 　　式中xd和xj分别为TMD与结构第j层质量相对于地面的位移。 　　2脉动风模拟 　　线性滤波法中的自回归模型因计算量小、速度快，广泛用于随机振动和时间系列分析中（胡广书，1997；P M Mignolet，et al，1987；S H Jeong，et al，1997；J S Owen，et al，2001）。本节采用线性滤波法中的自回归模型模拟作用在结构上的脉动风激励时程。 　　模型阶次估计方法多基于估计预测误差功率。通常估计的预测误差功率随阶次的增加而减小，但单纯追求最小预测误差求得的阶次往往偏高，出现过拟合现象，这虽然提高了谱估计的分辨率，但会产生虚假谱峰，并且估计方差也会增大，同时阶数过高也会增加计算量，影响运算速度；反之，如果选取的阶次过低会导致谱估计平滑，不能达到精度要求。根据试算，本文AR模型阶数取为p=5。 　　脉动风荷载P的功率谱密度函数可用矩阵表示为： （4） 　　式中 为n×n阶常量矩阵，其元素为： （5） 　　式中AiAj为i，j点的迎风面积。 　　由于脉动风压均值为零，其协方差矩阵等于互相关矩阵，对R进行Cholesky分解可得： 　　（6） 　　式中矩阵C为一下三角矩阵。 　　顺风向脉动风压向量T可表达为： 　　（7） 　　式中，u（t）为n个互不相关的高斯随机过程组成的向量。 　　利用脉动风压乘以迎风面积，就可以得到式（1）中的脉动风激励时程 P（t）。根据我国规范，一般取Davenport脉动风速谱作为激励谱Sij。当已知结构各点高度处的平均风速Vj风荷载体型系数Us及迎风面积Aj后，就可以得出自荷载谱： 　　（8） 　　如果考虑各个随机激励之间的相关性，则引入相关系数，便可形成荷载功率谱矩阵。 　　3随机结构风振响应 　　一旦自荷载谱（8）已知，在构造的虚拟随机过程激励下，可很方便得到TMD相对于结构的位移为： 　　（9） 　　TMD相对于结构的位移均方差为： 　　（10） 　　式中， （11） 　　通过求解模态响应矢量z的协方差矩阵，可以得到结构每一质点的加速度响应方差为： 　　（12） 　　即可求出加速度均方差响应及TMD相对于第j质量层加速度均方差 。式中L为对角元素均为1的下三角矩阵，rj第i振型参与系数； 　　分别为第j振型向量； 为结构的阻尼比； 为TMD的阻尼比； μ为TMD与结构的质量比。q为所取振型数，当只考虑第一振型时，