结构风致效应频域分析方法2数值算例.pdfVIP

第八届全国现代结构工程学术研讨会 结构风致效应频域分析方法(2)_数值算例 韩向科蒋露钱若军 (同济大学土木工程学院，上海200092) 握要：对于高层、大跨结构而言，风荷载往往是设计中的主要荷载。从摹奉的动力方程出发，把结构风致效应问题归结为动力方 程右端项(脉动风荷战)构造和方程求解这两方面的l’廿j题。运用新型高效的计算方法一虚拟激励法，进行了程序的编制， 并用于数值计算。 关键词；结构，风，风致效应。虚拟激励法 六、程序编制 究团队开发的AADS程序系统也进行了时程分析法的程序编制，与通用有限元软件的计算结果进行比较，取得 了良好的计算结果，并具有较高的效率。时程法的程序编制在诸多参考文献中都有阐述”捌，这里就不再赘述。 频域分析法中，根据前几章的理论公式和算法，采用Fortran语言，基于AADS程序系统进行了虚拟激励法 的程序编制。 依据上述虚拟激励法在抗风分析中的相关理论，结合研究团队共同开发的程序AADS进行了程序的编制。 与相关算例进行比较．取得了良好的计算结果。对包括网架网壳在内的各种结构进行数值模拟，结果表明虚拟激 励法是高效且有效的。 为了便于说明风振分析的程序编制，下面针对虚拟激励法的程序进行简单的介绍。 ．1、初始参数的计算。 ，． 载相关信息。 2、计算功率谱矩阵ls。。(动。 ： 3、乔列斯基(Cholesky)分解程序。。 将功率谱矩阵3阳(动按Choleskey分解方法转化成_s即(动=／：IALT；。 4、运用虚拟激励法计算响应功率谱矩阵。 疗。㈨=套日咖7瓯√石；得出u的功率谱矩阵ls明=∑玩‘以7。 k=l (5)峰值反应的估计。 r 计算位移响应最大值的均值或及反应的均方差盯，，-卅=f吒。 七、数值算例、 程序编制除了虚拟激励法的核心算法外，还涉及到大型方程组特征值求解、矩阵分解、数值积分等内容。 风致响应模块中，用户可以选择不同的风速自谱密度函数和相十函数。对于水平功率谱密度函数建议选取 系数按照规范或风洞试验确定后，可以通过数据读取直接进入程序。基本风压可按照规范或实测获得，程序按相 关规定自动生成基本风速。风压高度系数则只需给出支座高度，程序可按照规范自动进进行确定。相关计算理论 和参数取值，皆基于前面所述。 一 (一)程序正确性验证 文献【4】根据随机振动理论分析一个K6-8型单层球面网壳的风振响应，本文算例1采用相同算例进行计算， 634 工业建筑2008增刊 第八届全国现代结构工程学术研讨会 并与文献的计算结果进行对比，以此验证本文程序的正确性。 相干函数，基于准定常假定构造脉动风动力谱。计算时考虑其前150阶模态，积分上限为50rad／s。 图1网壳结构模型示意图 表1计算参数的取值 模型参数 参数取值 风荷载参数 参数取值 构件截面 m102×3 地貌类型 c类 屋面恒载 0．3I【N艋 基本风压 0．55kPa 支座约束条件 周边三向铰支 结构体形系数 按荷载规范 支座高度 lOm 风压高度系数 按荷载规范’ 跨度 401 矢高 4m 阻尼比 O．02 该网壳的前50阶自振频率见图2。由图可以看出该网壳自振频率非常密集，且成对出现。 选取Y方向对