结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序.doc

《结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、中心差分法计算单自由度体系动力反映的报告前言基于叠加原理的时域积分法与频域积分法一样，都假设结构在在全部反应过程中都是线性的。而时域逐步积分法只是假设结构本构关系在一个微小的时间步距内是线性的，相当于分段直线来逼近实际的曲线。时域逐步积分法是结构动力问题中研究并应用广泛的课题。中心差分法是一种目前发展的一系列结构动力反应分析的时域逐步积分法的一种，时域逐步积分法还包括分段解析法、平均常加速度法、线性加速度法、和法等。中心差分法(centraldifferencemethod)原理中心差分法的基本思路将运

2、动方程中的速度向量和加速度向量用位移的某种组合来表示，将微分方程组的求解问题转化为代数方程组的求解问题，并在时间区间内求得每个微小时间区间的递推公式，进而求得整个时程的反应。中心差分法是一种显示的积分法，它基于用有限差分代替位移对时间的求导（即速度和加速度）。如果采用等时间步长，（为常数），则速度与加速度的中心差分近似为(1)(2)用表示位移，离散时间点的运动为：，，体系的运动方程为(3)将速度和加速度的差分近似公式(1)和(2)代入(3)中得出在时刻的运动方程，将方程整理得到由和表示的两步法的运动方程

3、(4)：(4)这样就可以根据及以前的时刻的运动求得时刻的运动。中心差分法属于两步法，用两步法计算时存在起步问题，必须要给出相邻的两个时刻的位移值，才能逐步计算。对于地震作用下结构的反应问题和一般的零初始条件下的动力问题，可以用(4)直接计算，因为总可以假设初始的两个时间点（一般取）的位移等于零。但是对应于非零初始条件或零时刻外荷载很大时，需要进行一定的分析，建立两个起步时刻（即）的位移值。假设给定的初始条件为(5)根据初始条件来确定。根据中心差分公式(6)消去得到的公式：(7)其中零时刻加速度值可以由时

4、的运动方程得到即(8)这样就可以根据初始条件得到，然后再将初始条件应用于公式(4)中，逐步求出不同时刻的运动。中心差分法分析时的具体计算步骤：(1)基本数据准备与初始条件计算已知：初始位移、和初始荷载值来计算和(2)计算等效刚度和中心差分法计算公式中的系数因此中心差分法计算公式可以表示为：(3)根据及以前的时刻的运动求得时刻的运动(4)下一步计算中用代替，对于线弹性体系重复第3步计算步骤，对于非线性弹性体系，重复第2和第3计算步骤。以上的中心差分法逐步计算公式具有2阶精度，即误差；并且是有条件稳定的，稳

5、定条件为：式中，为结构的自振周期，对于多自由度体系则为结构的最小自振周期。算例本算例根据结构动力学48页算例3.1数据编写，稳定条件为对于一个单层框架结构，假设楼板刚度无限大，且结构质量集中于楼层，其质量、刚度、阻尼系数，对结构施加动力荷载假设结构处于线弹性状态，用中心差分法计算结构的自由振动反应。采用语言编程，并以单自由度体系为例进行计算,设初位移和初速度，取不同的步长分别计算，以验证中心差分法的稳定条件。先计算，由稳定条件，而，则，所以本次计算取＝0.2，0.1，0.05分别进行计算。计算结果与分析

6、1）当时，可以得到位移u，速度v，加速度ac的时程曲线如下：2）当时，可以得到位移u，速度v，加速度ac的时程曲线如下：3）当时，可以得到如下提示：不满足稳定条件:dt<=Tn/pi,请重新输入符合稳定条件的时间步长dt。附录%m=质量；k=刚度；c=阻尼；u0=初始位移；v0=初始速度；all_time=所用时间；P0=荷载幅值；dt=时间步长；%u=位移；v=速度；ac=加速度；ek=等效刚度；p=荷载；ep=等效荷载；t=时间；clear%A0=input('请按格式和顺序输入初始矩阵，如A0=[

7、m,k,c,u0,v0,all_time,P0,dt]=[924064100.050307300*30.05]=');A0=[924064100.05020730000.05];m=A0(1,1);k=A0(1,2);c=A0(1,3);u0=A0(1,4);v0=A0(1,5);all_time=A0(1,6);P0=A0(1,7);dt=A0(1,8);ifdt>2*sqrt(m/k)%判断时间步长是否满足稳定条件disp('不满足稳定条件:dt<=Tn/pi,请重新输入符合稳定条件的时间步长dt'

8、)returnelseif0