newmark法程序法計(jì)算多自由度體系動(dòng)力響應(yīng)

1、-. z用matlab編程實(shí)現(xiàn)法計(jì)算多自由度體系的動(dòng)力響應(yīng)用matlab編程實(shí)現(xiàn)法計(jì)算多自由度體系的動(dòng)力響應(yīng)一、法的根本原理Newmark-法是一種逐步積分的方法，防止了任何疊加的應(yīng)用，能很好的適應(yīng)非線性的反響分析。Newmark-法假定： (1-1) (1-2)式中，和是按積分的精度和穩(wěn)定性要求進(jìn)展調(diào)整的參數(shù)。當(dāng)=0.5，=0.25時(shí)，為常平均加速度法，即假定從t到t+t時(shí)刻的速度不變，取為常數(shù)。研究說(shuō)明，當(dāng)0.5，0.25(0.5+)2時(shí)，Newmark-法是一種無(wú)條件穩(wěn)定的格式。由式(2-141)和式(2-142)可得到用及，表示的，表達(dá)式，即有 (1-3) (1-4)考慮t+t時(shí)刻的振

2、動(dòng)微分方程為： (1-5)將式(2-143)、式(2-144) 代入(2-145)，得到關(guān)于ut+t的方程 (1-6)式中求解式(2-146)可得，然后由式(2-143)和式(2-144)可解出和。由此，Newmark-法的計(jì)算步驟如下：1.初始計(jì)算：1形成剛度矩陣K、質(zhì)量矩陣M和阻尼矩陣C；2給定初始值，和；3選擇積分步長(zhǎng)t、參數(shù)、，并計(jì)算積分常數(shù)，；4形成有效剛度矩陣；2.對(duì)每個(gè)時(shí)間步的計(jì)算：1計(jì)算t+t時(shí)刻的有效荷載：2求解t+t時(shí)刻的位移：3計(jì)算t+t時(shí)刻的速度和加速度：Newmark-方法是一種無(wú)條件穩(wěn)定的隱式積分格式，時(shí)間步長(zhǎng)t的大小不影響解的穩(wěn)定性，t的選擇主要根據(jù)解的精度確定。

3、本文用法計(jì)算的根本問(wèn)題四層框架構(gòu)造在頂部受一個(gè)簡(jiǎn)諧荷載的作用，力的作用時(shí)間=5s，計(jì)算響應(yīng)的時(shí)間為100s，分2000步完成。阻尼矩陣由Rayleigh阻尼構(gòu)造。具體數(shù)據(jù)如下列圖：圖一：構(gòu)造根本計(jì)算簡(jiǎn)圖計(jì)算法的源程序clcclearm=1,2,3,4;m=diag(m); %計(jì)算質(zhì)量矩陣k=800 -800 0 0; -800 2400 -1600 0; 0 -1600 4800 -3200; 0 0 -3200 8000;%計(jì)算剛度矩陣c=2*m*0.05*sqrt(k/m) ; t1=5; %力的作用時(shí)間nt=2000; %分2000步完成dt=0.01; %時(shí)間步長(zhǎng)alfa=0.25;

4、%=0.25beta=0.5; %=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); %a7=dt*beta; % d=zeros(4,nt); %初位移為0 v=zeros(4,nt); % 初速度為0 a=zeros(4,nt); % 初加速度為0 for i=2:nt t=(i-1)*dt; if (tt1), f=300*sin(6*pi*t)-50*cos(3*pi*t);0;0;

5、0; %力作用時(shí)間對(duì)構(gòu)造進(jìn)展加載 else f=0;0;0;0; %力作用時(shí)間外構(gòu)造不受力 end ke=k+a0*m+a1*c; % 有效剛度矩陣fe=f+m*(a0*d(:,i-1)+a2*v(:,i-1)+a3*a(:,i-1)+c*(a1*d(:,i-1)+a4*v(:,i-1)+a5*a(:,i-1); % t+t時(shí)刻的有效荷載 d(:,i)=inv(ke)*fe; %求解t+t時(shí)刻的位移 a(:,i)=a0*(d(:,i)-d(:,i-1)-a2*v(:,i-1)-a3*a(:,i-1); %計(jì)算t+t時(shí)刻的加速度 v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(

6、:,i); %計(jì)算t+t時(shí)刻的速度endT=0:dt:19.99; %離散系統(tǒng)dt為采樣周期 19.99為終端時(shí)間close allfigure %控制窗口數(shù)量plot(T,d) %繪制位移函數(shù)圖像title() %添加標(biāo)題為各質(zhì)點(diǎn)位移總圖 legend(,) %添加圖例的標(biāo)注*label(s) %對(duì)*軸進(jìn)展標(biāo)注為時(shí)間sylabel(m) %對(duì)y軸進(jìn)展標(biāo)注為位移mgrid %顯示畫圖中的個(gè)網(wǎng)線figure plot(T,a) %繪制加速度函數(shù)圖像title() %添加標(biāo)題為各質(zhì)點(diǎn)加速度總圖legend(,)*label(s) %對(duì)*軸進(jìn)展標(biāo)注為時(shí)間sylabel(m/s2) %對(duì)y軸進(jìn)展標(biāo)注為加速度m/s2grid%figure plot(T,v) %繪制速度函數(shù)圖像title() %添加標(biāo)題為各質(zhì)點(diǎn)速度總圖legend(,)*label(s) %對(duì)*軸進(jìn)展標(biāo)注為時(shí)間sylabel(m/s2) %對(duì)y軸進(jìn)展標(biāo)注為加速度m/s2grid計(jì)算結(jié)果截圖最后程序分別計(jì)算出四個(gè)質(zhì)點(diǎn)的位移、速度、加