结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序(优选.)VIP免费

word.(即速度和加速度)。如果(1)(2)(3)将方程整理得到'由'中心差分法计算单自由度体系动力反映的报告■*%刖言基于叠加原理的时域积分法与频域积分法一样，都假设结构在在全部反应过程中都是线性的。而时域逐步积分法只是假设结构本构关系在一个微小的时间步距内是线性的，相当于分段直线来逼近实际的曲线。时域逐步积分法是结构动力问题中研究并应用广泛的课题。中心差分法是一种目前发展的一系列结构动力反应分析的时域逐步积分法的一种，时域逐步积分法还包括分段解析法、平均常加速度法、线性加速度法、i-'乜JT...I一!<■.和h小：ii法等。中心差分法(centraldifferencemethod)原理中心差分法的基本思路将运动方程中的速度向量和加速度向量用位移的某种组合来表示，将微分方程组的求解问题转化为代数方程组的求解问题，并在时间区间内求得每个微小时间区间的递推公式，进而求得整个时程的反应。中心差分法是一种显示的积分法，它基于用有限差分代替位移对时间的求导采用等时间步长，血2二加(加为常数)，则速度与加速度的中心差分近似为2At..气-厂叫气二7At用•■表示位移，离散时间点的运动为：u严讥用，u严讥》),&产诃》)〔心卽体系的运动方程为rnii(t)++ku(6=P(t)将速度和加速度的差分近似公式⑴和⑵代入⑶中得出在：时刻的运动方程,和'I表示的两步法的运动方程(4)：这样就可以根据：及以前的时刻的运动求得：I时刻的运动。中心差分法属于两步法，用两步法计算时存在起步问题，必须要给出相邻的两个时刻的位移值，才能逐步计算。对于地震作用下结构的反应问题和一般的零初始条件下的动力问题，可以用(4)直接计算，因为总可以假设初始的两个时间点(一般取；--】)的位移等于零。但是对应于非零初始条件或零时刻外荷载很大时，需要进行一定的分析，建立两个起步时刻(即i—L-」)的位移值。假设给定的初始条件为word.it消去"得到「|的公式：it2=u0-Atu0+—其中零时刻加速度值一可以由「:时的运动方程得到即这样就可以根据初始条件得到'''，然后再将初始条件应用于公式(4)中,中心差分法分析时的具体计算步骤：(1)基本数据准备与初始条件计算已知：初始位移'和初始荷载值儿来计算九和「I(7)(8)逐步求出不同时刻的运(2)计算等效刚度和中心差分法计算公式中的系数,mc2ma=fc-—7因此中心差分法计算公式可以表示P.-aui-but(3)根据：及以前的时刻的运动求得：I时刻的运动Au0=«(o)'u0—u(O)根据初始条件来确定'I。根据中心差分公式word.（4）下一步计算中用…I代替，对于线弹性体系重复第3步计算步骤，对于非线性弹性体系，重复第2和第3计算步骤。以上的中心差分法逐步计算公式具有2阶精度，即误差■W；并且是有条件稳定的，稳定条件为:Af<—?r式中，匚为结构的自振周期，对于多自由度体系则为结构的最小自振周期。算例本算例根据结构动力学48页算例3.1数据编写，稳定条件为dt<=0.16s对于一个单层框架结构，假设楼板刚度无限大，且结构质量集中于楼层，其质量T沁；：注、刚度二八皿比阻尼系数—-，对结构施加动力荷载一-二；”一」心昭汛假设结构处于线弹性状态，用中心差分法计算结构的自由振动反应。采用；1汀■语言编程，并以单自由度体系为例进行计算，设初位移!''■N：I和初速度取不同的步长分别计算，以验证中心差分法的稳定条件At鼻。人T2先计算At，由稳定条件A<^n=，而3兀①n，所以本次计算取At=0.2,0.1，0.05分别进行计算。计算结果与分析1）当•I111'时，可以得到位移u,速度v,加速度ac的时程曲线如下:■■111i1iii•i11?斗Eieio12u111111f1ii11■1ieaio12i2）当」"■时，可以得到位移u,速度v,加速度ac的时程曲线如下:时间«2word.3）当时，可以得到如下提示：不满足稳定条件:dtv=Tn/pi,请重新输入符合稳定条件的时间步长dt。酒戾.的时程曲161620.C運lf时■间啊加题虫航的时程由«8亡〕It间（=〕加速虔检释产1姑理2:6S1D12IE13201FEBSS1D12U16时间闾担速扈n=的L•程曲孰fori=3:nn%将时间分步，采用等时间步长；%计算t的向量长度，得出步数；%设定存储u的矩阵；%设定存储v的矩阵；%设定存储ac的矩阵；%赋值向量第2项为u0；%赋值向量第2项为v0；%求出初始加速度ac0；%计算初始条件u-...