Newmark数值分析法matlab程序

clc
clear
M=[1,2,3,4];
m=diag(M); %质量对角线矩阵
k=[ 800 -800 0 0 ;
-800 2400 -1600 0 ;
0 -1600 4800 -3200 ;
0 0 -3200 8000 ]; %计算刚度矩阵
c=0.05m+0.02k;%c=2m0.05*sqrt(k/m) ;
f0= 100;
t1=5; %力的作用时间
nt=2000; %分2000步完成
dt=0.01; %时间步长
alfa=0.25; %γ=0.25
beta=0.5; %β=0.5

a0=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=dtbeta; %
d=zeros(4,nt); %初位移为0
v=zeros(4,nt); % 初速度为0
a=zeros(4,nt); % 初加速度为0
for i=2:nt
t=(i-1)dt;
if (t<t1),
f=[f0sin(4pit/t1);0;0;0]; %f=[300sin(6pit)-50cos(3pit);0;0;0]; %力作用时间内对结构进行加载
else
f=[0;0;0;0]; %力作用时间外结构不受力
end
ke=k+a0m+a1c; % 有效刚度矩阵
fe=f+m(a0d(:,i-1)+a2v(:,i-1)+a3a(:,i-1))+c(a1d(:,i-1)+a4v(:,i-1)+a5a(:,i-1));
% t+dt时刻的有效荷载
d(:,i)=ke\fe;%d(:,i)=inv(ke)fe; %求解t+dt时刻的位移
a(:,i)=a0(d(:,i)-d(:,i-1))-a2v(:,i-1)-a3a(:,i-1); %计算t+dt时刻的加速度
v(:,i)=v(:,i-1)+a6a(:,i-1)+a7*a(:,i); %计算t+dt时刻的速度
end

T=(0:dt:19.99); %离散系统dt为采样周期 19.99为终端时间
close all

figure %控制窗口数量
plot(T,d(1,:)) %绘制位移函数图像
title(‘各质点位移总图’) %添加标题为各质点位移总图
xlabel(‘s’) %对x轴进行标注为时间（s）
ylabel(‘m’) %对y轴进行标注为位移（m）
grid %显示画图中的个网线

figure
plot(T,a) %绘制加速度函数图像
title(‘加速度总图’) %添加标题为各质点加速度总图
legend(‘a1’,‘a2’,‘a3’,‘a4’)
xlabel(‘s’) %对x轴进行标注为时间（s）
ylabel(‘m/s^2’) %对y轴进行标注为加速度（m/s2）
grid

figure
plot(T,v) %绘制速度函数图像
title(‘速度总图’) %添加标题为各质点速度总图
legend(‘v1’,‘v2’,‘v3’,‘v4’)
xlabel(‘s’) %对x轴进行标注为时间（s）
ylabel(‘m/s)’) %对y轴进行标注为加速度（m/s）
grid

Newmark数值分析法matlab程序相关推荐

matlab灰关联度,[转载]灰色关联和加权关联分析法matlab程序
%by allen %灰色关联分析法和灰色加权关联分析法matlab程序,其区别主要在于求关联度是使用平均值法还是加权法 %平均值法为一般普通求法,加权值法可以根据侧重点不同进行分析 %matlab6 ...
matlab关联度排序模型,灰色关联和加权关联分析法matlab程序
%by allen %灰色关联分析法和灰色加权关联分析法matlab程序,其区别主要在于求关联度是使用平均值法还是加权法 %平均值法为一般普通求法,加权值法可以根据侧重点不同进行分析 %matlab6 ...
二次指数平滑法matlab程序_二次指数平滑法程序
二次指数平滑法程序线性指数平滑法 Matlab 程序,代码如下: 注: Data- 原始数据 s- 一次和二次平滑结果 at- 预测式中的 a 参数 bt- 预测式中的 b 参数 y1- 预测结果 ...
matlab用进退法写程序,进退法matlab程序
极小值点包含于区间 [ x(1) , x(3) ]或[ x(3) , x(1) ] (3)算法的 MATLAB 实现在 MATLAB 中编程实现的进退函数为: min JT 功能:用进退法求解.... ...
标号法matlab程序,matlab程序大全
上例 Matlab 程序如下: %使用表8.4 给出的小波基 ld=[0 0.0... matlab经典代码大全_计算机软件及应用_IT/计算机_专业资料.哈哈哈 MATLAB 显示正炫余炫图:plo ...
五步移相法matlab程序_用MATLAB软件对逆变电路移相调压进行仿真,30°,60°,90°的波形,并用分段分析法分析其原理_学小易找答案...
[简答题]Translate the following ad into (1) a TV comercial, and (2) a factual introduction. Marie Franc ...
《摄影测量学基础》光束法MATLAB程序
实验目的运用编程的方法进一步掌握光束法双向解析摄影测量. 实验要求用Matlab编写光束法双向解析摄影测量程序. 实验数据左片外方位元素近似值:Xs1 = 4999.0;Ys1 = 4999.0 ...
五步移相法matlab程序_基于结构光的相移法三维重建matlab
一.基本原理: 正弦条纹打在三维物体上,CCD记录到的条纹由于受到三维物体高度的调制而发生扭曲,扭曲的条纹(deformed fringe)实质上为原始条纹在物体具有高度存在的位置有了附加相位,各点的 ...
阈值法matlab程序,遗传算法优化BP神经网络权值和阈值的通用MATLAB源码
遗传算法优化神经网络有两种情况,一种是把训练好的神经网络作为黑箱函数,用遗传算法搜索该黑箱函数的最大值,另外一种情况,则是把遗传算法用于神经网络的训练,充分利用遗传算法全局搜索的特性,得到一个初始的权 ...

Newmark数值分析法matlab程序

Newmark数值分析法matlab程序相关推荐

最新文章

热门文章