钢结构设计软件USSCAD地震反应谱计算研究.pdfVIP

钢结构 2012年增刊第27卷总第167期 杜新喜邓文刘茂青 (武汉大学土木建筑工程学院，武汉430072) 提要：在钢结构结构设计中，地震作用的计算方法一般以Duhamel积分运算和直接积分法两种方法较为常见。本文通 过MATLAB编程计算，针对同一个实际的地震波文件，对两种方法的精度和效率做出了定量分析。研究结果表明，从 精度和效率上比较，直接积分法具有快速高精度的特点。同时，这一计算方法也在基于AutoCAD二次开发的通用钢结 构设计软件(USSCAD)中得到了实现。 关键词：地震反应谱；直接积分法；Duhamel积分；MATLAB；USSCAD 一引言 由于AutoCAD的强大功能、易用性和开放的体系结构，其已成为世界上最广泛的CAD设计软件 包之一…。结构设计行业因其独特的行业标准，在使用AutoCAD的过程中根据自身特点进行了定制和 的。但这些软件在计算结构地震作用时，都无法对针对某一具体地震波而进行地震反应计算。考虑到当 前USSCAD等软件在这一方面的缺陷，为此有必要对各种计算地震反应谱的分析方法做出定量比较， 以供结构设计人员在算法设计上做出正确抉择。 积分法和直接积分法做出了分析比较。同时将两种结果与太平洋地震工程研究中心(PEER)的BISPEC 软件计算结果进行了对比。研究结果表明，当前USSCAD软件中的地震反应谱计算精度高、速度快。 二Duhamel积分法 1Duhamel积分法的原理 对于单自由度地震作用的求解，根据质点力的平衡方程，对于单自由度体系，有如下的二阶线性非 齐次微分运动方程‘3]： 37+2考coJ+c02x=一《 (2—1) 舯∞=后“：去 方程(2．1)的数学上的解析解即为杜哈梅(Duhamel)积分式，当阻尼较小时其形式如下： 帮助文档【5]。 2Duhamel积分法的MATLAB实现 198 钢结构 2012年增刊第27卷总第167期 的部分程序代码及解释： for ％从0．1～5s取100个点进行反应谱计算 k=l：100； T(k)=0．05术k； ％周期 w=2木pi／T(k)； ％圆频率 ％加速度文件的时间数列 t=0：dt：(1ength(at)-1、木dt； ％形成卷积数组hw hw=exp(一1木kesi术w木t)．半sin(w水t)； y=conv(a_t，hw)．水(dt半w)；％卷积计算 Sa(k)=max(abs(y))； ％绝对值取最大得到最大地震反应 end； ％结束对周期的循环求解 三直接积分法 直接积分法有多种算法，本文采用线性加速度法，即假设在At时间间隔内，加速度线性变化。在 实际工程中，地震时地面运动加速度艺(f)是关于时刻f=胛×At的一组离散数据∞1。因此，采用线性加 地面运动在时刻f+△f的速度量g@+△f)和位移Xg(H△f)可以按照泰勒展开式展开为吐 {耋：二：：三耋善；二{吴：：主：竺二妻{三：；i耄：二三竺+吉。△，，，x。+△，， c3—1， 硝Ⅲf)-硝f)+半引卅半硝Ⅲ力 (3-2) 蹦Ⅲ归叫卅(％∽+竽铋卅竿叫m幻 同理，对于相对加速度戈， 仍旧假设为线性变化‘61，则有以下关系式： hM刊+半艳，+半艄心，