双目标控制选择Rayleigh阻尼系数方法.docVIP

双目标控制选择Rayleigh阻尼系数方法 　　摘要： 合理选择确定Rayleigh阻尼矩阵比例阻尼系数的振型频率对于准确计算场地地震响应有重要影响。提出以土层表面的位移反应峰值误差与加速度反应峰值误差的加权平方和为控制目标，基于完全二次项平方根振型组合方法建立控制目标与Rayleigh阻尼系数的函数关系，以控制目标最小为原则形成求解Rayleigh阻尼系数的方法。随机选择28条具有代表性的地震波，通过对某典型深厚土层模型进行地震反应分析，并与其他研究成果进行对比，验证所建议方法的精度及适用性。关键词： 地震响应分析； Rayleigh阻尼； 阻尼系数； 双目标控制； 深覆盖土层 　　中图分类号： TU311.3； O322文献标志码： A文章编号：2017 　　DOI：10.16385/ki.issn2017.04.015 　　1概述 　　结构体系阻尼的形成涉及较多的影响因素，其形成机理极为复杂，目前尚难以定量。基于黏滞阻尼模型的Rayleigh阻尼矩阵由于数学处理的方便而得到广泛应用。Rayleigh阻尼矩阵是由结构体系的质量矩阵和刚度矩阵的线性组合而成C=αM+βK（1）式中M，C和K分别是结构体系的质量矩阵、阻尼矩阵和刚度矩阵，α为质量比例系数，β为刚度比例系数。根据振型的正交性，由Rayleigh阻尼矩阵计算所得的动力系统各阶振型阻尼比可表示为ζi=α2ωi+βωi2（2）式中ωi是第i阶振型的圆频率，ζi为第i阶振型阻尼比的近似计算值。通过选取结构的两个振型频率（ωj，ωk）及对应振型阻尼比（ζj，ζk），可确定质量比例系数和刚度比例系数。这样，除了这两阶振型之外，由Rayleigh阻尼矩阵确定的其他振型阻尼比都与真实的振型阻尼比不一致，由此，会高估或低估这些振型对结构动力响应的贡献，使得?Y构动力响应的计算结果出现偏差。 　　工程计算和研究工作中广泛使用Rayleigh阻尼矩阵，但是关于Rayleigh阻尼矩阵对结构动力响应计算结果带来影响的认识并不深入。其原因之一是：以往大量工程结构的基频高于或接近于外部动力作用主要分量的激励频率，而结构的低阶振型对结构总动力响应起主要控制作用，因此往往选取结构的两个低阶振型来确立Rayleigh阻尼矩阵，由于低阶振型的计算阻尼比与真实阻尼比相一致，Rayleigh阻尼矩阵带来的计算误差就很小。楼梦麟等[12]和马俊玲等[3]的研究表明：当结构基频远低于外部动力作用主要分量的激励频率时，仅选择结构的两个低阶振型对应的频率和阻尼比来计算Rayleigh阻尼矩阵的比例系数，将使结构的动力反应被低估。为改进Rayleigh阻尼矩阵的数值建模，有研究工作[47]选择结构或者地震波的某些特定频率（及对应的阻尼比）来确定Rayleigh阻尼矩阵的比例系数，这些特定频率包括结构基频、其他对结构反应有显著贡献的高阶振型频率、激励时程的反应谱（或傅氏谱）的峰值或谱曲线重心对应的频率、场地特征频率等。 　　为避免选取参考频率的经验性和任意性，刘红石[8]提出采用加权最小二乘法从振型阻尼比误差最小的角度计算Rayleigh阻尼矩阵的比例系数，建议以权重系数反映各阶振型对动力反应贡献的差异，但没有提出相应的权重系数。Yang等[9]提出以峰值应变能系数γ2iMiS2vi（其中γi为第i阶振型的参与系数，Mi为第i阶振型的广义质量，Svi为第i阶振型的速度反应谱值）作为权重系数，而潘旦光等[10]建议以γ2iφ2kiS′2ui（其中φki为第i阶振型向量的第k个元素，S′ui为第i阶振型的位移反应谱对振型阻尼比的导数）作为权重系数。这些方法都是双参数优化方法，董云等[11]提出基于结构基频的单参数优化方法，不仅考虑结构基本振型对结构地震反应的重要性，也使得计算过程简化，便于工程应用。除加权最小二乘法外，Spears等[12]提出以控制∑ni=1MeiAi（其中Mei为第i阶振型的有效质量，Ai为频率与第i阶振型频率一致的单自由度体系的最大加速度响应）接近于零为目标，通过迭代方法求解Rayleigh阻尼矩阵的比例系数。 　　第4期王淮峰，等： 双目标控制选择Rayleigh阻尼系数的方法振 动 工 程 学 报第30卷重大工程场地的地震安全性评价工作中通常需要对土层的地震反应进行分析。而深覆盖土层场地基频往往较低（如上海约300 m厚的土层，其基频约0.5 Hz），远低于常见基岩地震波主要分量的激励频率。因此，要较准确地计算此类场地的地震响应，首先需要合理选择确定Rayleigh阻尼矩阵比例系数的参考振型频率。本文作者[13]曾基于完全二次项平方根振型组合方法（CQC），以土层表面加速度反应峰值的误差为控制目标，同时将土层基频作为参