newmark法程序法计算多自由度体系动力响应

《newmark法程序法计算多自由度体系动力响应》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.用matlab编程实现法计算多自由度体系的动力响应用matlab编程实现法计算多自由度体系的动力响应一、法的基本原理Newmark-b法是一种逐步积分的方法，避免了任何叠加的应用，能很好的适应非线性的反应分析。Newmark-b法假定：(1-1)(1-2)式中，b和g是按积分的精度和稳定性要求进行调整的参数。当b=0.5，g=0.25时，为常平均加速度法，即假定从t到t+Dt时刻的速度不变，取为常数...。研究表明，当b≥0.5，g≥0.25(0.5+b)2时，Newmark-b法是一种无条件稳定的格式。由式(2-141)和式(2

2、-142)可得到用及，，表示的，表达式，即有(1-3)(1-4)考虑t+Dt时刻的振动微分方程为：(1-5)将式(2-143)、式(2-144)代入(2-145)，得到关于ut+Dt的方程(1-6)式中求解式(2-146)可得，然后由式(2-143)和式(2-144)可解出和。由此，Newmark-b法的计算步骤如下：1.初始计算：（1）形成刚度矩阵[K]、质量矩阵[M]和阻尼矩阵[C]；（2）给定初始值，和；（3）选择积分步长Dt、参数b、g，并计算积分常数，，，，，，，；（4）形成有效刚度矩阵；...2.对每个时间步的计算：（1

3、）计算t+Dt时刻的有效荷载：（2）求解t+Dt时刻的位移：（3）计算t+Dt时刻的速度和加速度：Newmark-b方法是一种无条件稳定的隐式积分格式，时间步长Dt的大小不影响解的稳定性，Dt的选择主要根据解的精度确定。二、本文用法计算的基本问题四层框架结构在顶部受一个简谐荷载的作用，力的作用时间=5s，计算响应的时间为100s，分2000步完成。阻尼矩阵由Rayleigh阻尼构造。具体数据如下图：...图一：结构基本计算简图二、计算法的源程序clcclearm=[1,2,3,4];m=diag(m);%计算质量矩阵k=[800-8

4、0000;-8002400-16000;0-16004800-3200;00-32008000];%计算刚度矩阵c=2*m*0.05*sqrt(k/m);t1=5;%力的作用时间nt=2000;%分2000步完成dt=0.01;%时间步长alfa=0.25;%g=0.25beta=0.5;%b=0.5a0=1/alfa/dt/dt;%a1=beta/alfa/dt;%...a2=1/alfa/dt;%a3=1/2/alfa-1;%a4=beta/alfa-1;%a5=dt/2*(beta/alfa-2);%a6=dt*(1-beta

5、);%a7=dt*beta;%d=zeros(4,nt);%初位移为0v=zeros(4,nt);%初速度为0a=zeros(4,nt);%初加速度为0fori=2:ntt=(i-1)*dt;if(t

6、-1)+a4*v(:,i-1)+a5*a(:,i-1));%t+Dt时刻的有效荷载d(:,i)=inv(ke)*fe;%求解t+Dt时刻的位移a(:,i)=a0*(d(:,i)-d(:,i-1))-a2*v(:,i-1)-a3*a(:,i-1);%计算t+Dt时刻的加速度v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(:,i);%计算t+Dt时刻的速度endT=[0:dt:19.99];%离散系统dt为采样周期19.99为终端时间closeallfigure%控制窗口数量plot(T,[d])%绘制位移函数图像ti

7、tle('¸÷ÖʵãλÒÆ×Üͼ')%添加标题为各质点位移总图legend('ÖʵãÒ»','Öʵã¶þ','ÖʵãÈý','ÖʵãËÄ')%添加图例的标注xlabel('ʱ¼ä(s)')%对x轴进行标注为时间（s）ylabel('λÒÆ(m)')%对y轴进行标注为位移（m）grid%显示画图中的个网线figureplot(T,[a])%绘制加速度函数图像...title('¸÷Öʵã¼ÓËÙ¶È×Üͼ')%添加标题为各质点加速度总图legend('ÖʵãÒ»','Öʵã¶þ','ÖʵãÈý','Öʵã

8、ËÄ')xlabel('ʱ¼ä(s)')%对x轴进行标注为时间（s）ylabel('¼ÓËÙ¶È(m/s^2)')%对y轴进行标注为加速度（m/s2）grid%figureplot(T,[v])%绘制速度函数图像title('¸÷Ö