单层球面网壳结构风致响应计算的块里兹向量法.pdf

工 No．2 第20卷第2期 振 动 程 学 报 V01．20 2007年4月 ofVibration Apr．2007 Journal Engineering 单层球面网壳结构风致响应计算的块里兹向量法 李方慧1，倪振华2，沈世钊3 (1．黑龙江大学建筑工程学院，黑龙江哈尔滨150086；2．汕头大学土木工程系，广东汕头515063I 3．哈尔滨工业大学土木工程学院，黑龙江哈尔滨150090) 摘要；单层球面网壳结构自振频率密集，且高阶模态可能对结构的风致响应有较大贡献。利用传统模态叠加法计算 结构风致响应时需截取较多的模态，导致计算效率较低。本文提出的块里兹向量法所生成的向量都是被荷载所激 发的，因此仅使用少数向量即可获得准确的响应结果。块里兹向量迭代所需的多荷载模式可以利用本征正交分解 (POD)技术获得，这种技术将脉动风压场分解为与时间无关的协方差模态以及与空间位置无关的主坐标。根据风 洞实验同步测量的风压数据，利用本文计算了0．1和0．2矢跨比的单层球面网壳的风致响应。通过与传统模态叠加 法计算结果的对比分析说明了本文方法的有效性。 关键词：风致响应；块里兹向量法；本征正交分解法；风洞试验；单层网壳结构 中图分类号：TU393文献标识码：A 文章编号：1004—4523(2007)02．0128—05 引 言 1块里兹向量法 单层球面网壳结构的自振特性和风致响应特征 通常结构脉动风致响应分析的模态叠加法选用 与其它类型屋盖有较大不同。Nakayama首先指出的振型仅仅与结构的自振特性有关，而没有考虑风 该类屋盖可能存在高阶模态对结构风致响应贡献很 荷载的特性。在计算过程往往很难确定选用多少阶 大的现象[1]，一些学者进行了深入研究。Uematsu对模态可以达到要求的计算精度，如果截取模态数太 两种不同跨度、三个矢跨比的单层球面网壳结构在 少可能忽略高阶模态的贡献。此外，对于选用的模态 时域上计算风致响应，总结风致响应特征并采用阵 缺乏挑选标准，不管其对响应贡献大小均被包含在 风因子方法估计风荷载[2]。何艳丽采用模态补偿方 内，对响应贡献很小的模态也可能被考虑在内影响 法考虑了高阶模态对响应的贡献，以跨度50ITI的 了计算效率。为了克服传统模态叠加法的一些不足， K6—8型单层球面网壳为例分析[3]。黄明开采用里兹需要探寻一种高效计算结构风致响应的方法。1982 向量叠加法计算球壳屋盖风致响应，在里兹向量生 年，E．L．Wilson、袁明武等人发展了与荷载相关的 成过程中初始迭代向量是根据平均风荷载产生的静 里兹向量直接叠加法口]，该方法是计算结构动力响 力变形获得“]。上述针对高阶模态可能对响应贡献 应分析的一种快速算法，这些与质量矩阵正交的里 较大的情况给出的一些处理途径。 兹向量生成过程较为简单。生成的里兹向量不仅与 本文对屋盖风致响应计算的块里兹向量法进行 结构自振特性有关而且与作用结构上的动力荷载的 研究，在分析过程中考虑了脉动风荷载特性。利用本 空间分布特性相关，获得的少量模态均是被荷载激 征正交分解法(POD)提取脉动风荷载的多空间分发出都是对响应有贡献的[6~盯。该方法被广泛应用 布模式，初始的块里兹向量根据荷载空间分布模式 于地震响应分析，地震荷载可以明确的表示为一个 的静变形获得，POD技术为初始块里兹向量的确定 空间向量(与质量分布相关)和时间函数(与加速度 提供了有利工具。本文以0．1和0．2矢跨比K8型单 相关)的乘积，而对于脉动风荷载这种单荷载模式的 层球面网壳结构为例，利