一种Duhamel积分的递推算法及在结构风振响应时域计算中的应用.pdfVIP

第十四届全国结构风工程学术会议论文集 一种Duhamel积分的递推算法及在结构风振响应时域计算中的应用 王国砚 钱浩 同济大学航空航天与力学学院，上海200092 摘要本文建立一种线性结构动力响应时域计算的基于Duhamel积分的递推算法．首先，针对单自由度结构模 型建立该递推算法的一般算式，并验证了该算式的计算精度和计算效率。然后，针对多自由度结构模型建立基于振 型分解法的递推算式，并通过算例验证了该算式的正确性．最后，通过一个数值算例演示了本文算法在结构风振响 应时域计算中的应用． 关键词Duhamel积分递推算法风振响应时域计算 1 引 言 由于时域方法的一些独特优点，近年来结构风振响应的时域计算日益受到关注。例如，Chen、 细长结构颤振响应计算的时域计算方法；美国国家标准技术研究所(NIST)51：发出了用于计算刚、柔 性结构风振时域响应的数据库辅助设计软件(DAD)扪。 用较多的数值方法4)，5、。该方法的优点是适用于线性和非线性结构；缺点是，由于线性加速度假定， 使得这些方法的解在本质上属于近似解，当积分步数较多时，计算结果会出现周期变长和振幅变小 等现象41。为此，许多学者致力于时域数值算法的研究。例如，林家浩6’和钟万勰7’先后提出一种时 程积分的高精度数值方法。 种递推方法，并认为对于分段线性激励而言，该方法计算结果为精确解。随后，Dempsey等9)在 and 的算式一致。Sun 本文针对分段线性激励的情形建立一种Duhamel积分的递推算法，具有精度高和效率高等特点。 实际工程结构的动力学问题应该说仍以线性为主；在风洞试验中，由于风荷载是在每个时间步上采 样的，这样的荷载可假定为是分段线性的。在这种情况下，用本文方法可快速求出高精度的结果。 2基于Duhamel积分的递推算法 图1所示为一线性单自由度结构，动力学方程可以用下式表达： 夕+2分矽+r02y=FO) (1) 在零初始条件下，激励F(t)产生的位移响应可由下面的Duhamel积分表达： 1 ． y(t)=÷【F(力e一矧卜f)sinaf(t一力df (2) 通过式(2)可以算出任意时刻t的位移响应，以此作为递推算法在该时间步的初值，再做 算下去。这就是Duhamel积分递推算法的主要思想。 1041 第十四届全国结构风工程学术会议论文集 L——一啪 ■∞ 图1单自由度结构示意图 本文根据上述递推算法的思想，将，州时刻的位移响应表不成F式： 少(‘+。)=X(tl+1)一艺(‘+。) (3) 式中 X(‘+。)=4(‘+。)·厶(‘+。) (3a) 4(“)2古e讽sin∥“ (3a-I) 厶(“)=t+lF(力e缈cosco订r (3a-2) Y2(tf+1)=A2(‘+1)·12(ff+1) (3b) A2(t,+，)2古e嘲+Icosd“ (3b-1) a)订r