杜哈梅积分设计研究李伯维.docVIP

杜哈梅积分程序设计研究 李黎1，*李伯维1,吴洁2 （1.华中科技大学土木工程与力学学院 湖北 武汉430074；2.华中科技大学数学与统计学院 湖北 武汉430074） 摘 要：为了得到复杂荷载作用下结构的动力响应，本文基于杜哈梅积分方法，编写了FORTRAN数值积分程序。并应用解析解验证了程序的正确性。最后采用实测地震波，利用程序对于结构进行地震响应分析。 关键词：单自由度，杜哈梅积分，数值积分，程序设计 中图分类号: 文献标识码：A Study on Duhamel Integral programming LI Li1,LI Bowei1,Wu Jie2 (School of Civil Engineering and Mechanics, Huazhong University of Science Technology, Wuhan 430074, China) Abstract:In order to get the reaction of structure when under subjected to complex load. Based on Duhamel numerical integration method and programming by FORTRAN, verified the correctness of the program by analytical solutions that have been obtained. And analyzed a structure under real seismic action by this program. Key words: Single Degree-Of-Freedom;Duhamel numerical integration;Numerical integration;programming 1 引言 物体在持续外力作用下发生的振动称为受迫振动[1-3],结构在地震作用下的震动，也可认为是受迫振动的一种。对于受迫振动，我们可以通过求解结构平衡微分方程的方法得到结构的位移函数，也可以通过杜哈梅积分的方法进行求解[4]。然而微分方程的解法仅仅适用于以正弦荷载为代表的很小一部分荷载作用下的结构响应计算，对于一般的荷载作用，应用杜哈梅积分进行响应计算，变得更为有效。 对于杜哈梅积分，大多数荷载表达式较为复杂或者无法表达式，无法应用常规解析的方法，较为可行的方法是应用数值算法，但是要保证足够的精度。对于人工而言，数值算法工作量大，效率太低。但在计算机普遍应用的今天，我们可以通过计算机进行数值计算，可以高效地完成分析，甚至可以实现时程分析。 本文通过设计杜哈梅积分的程序，对结构力学教材中出现的杜哈梅积分进行数值化计算，可以输出结构的位移波形以及峰值。程序设置了地震分析模块，只需改变输入参数，并将荷载序列更换为加速度序列，即可完成地震作用下结构响应的分析。 2 程序设计与验证 2.1 杜哈梅积分程序设计 2.1.1杜哈梅积分 图1 杜哈梅积分示意图 设体系在t=0时处于静止状态。然后有瞬时冲量S作用。由于冲量S的作用，体系将生初速度0，但初位移为零。 式2-1 则有 式22 如图1所示，在时刻t=τ，冲量微元S=Pτ)dτ，将之代入式2-2并注意冲量P(τ)dτ引起的振动开始于时刻t=τ。 则有 式2-3 式2-4 m为体系集中质量，ω为体系自振频率； P0(τ)为荷载函数； ξ为体系阻尼比。 上式即为杜哈梅积分表达式。 由于地震记录的信息均为加速度序列，对于结构进行地震作用下分析，需对于上式进行适当变换，推导如下： 首先对于体系建立平衡微分方程如下式： 式2-5 式中：m为单自由度体系质量； 为地面加速度； 分别为体系质点位移、速度、加速度。 对式2-5进行移项，可以得到与强迫振动形式等效的运动方程： 式2-6 那么地震作用可以看作是在体系质量上集中作用一个的集中力，那么将2-4式中P0(τ)换为即可得到施加了加速度序列荷载的杜哈梅积分表达式，如式2-7所示： 式2-7 式中：所有符号意义同前，仅将荷载函数换为加速度函数； ξ为体系阻尼比。 通过式2.1.1即可处理输入荷载为加速度序列情况下结构的计算