没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB用中心差分法实现多自由度的振动问题
5星 · 超过95%的资源 需积分: 50 121 下载量 155 浏览量
2012-12-06
21:16:52
上传
评论 7
收藏 3KB TXT 举报
温馨提示
试读
2页
MATLAB用中心差分法实现多自由度的振动问题。非常好用。
资源推荐
资源详情
资源评论
clear all
M=input('输入质量矩阵M=');
C=input('输入阻尼矩阵C=');
K=input('输入刚度矩阵K=');
P=input('输入外部载荷矩阵P=');
u0=input('输入初始位移矩阵u0=');
v0=input('输入初始速度矩阵v0=');
time=input('输入模拟时间time=');
dt=input('输入时间步长dt=');
[m,h]=size(M); %m表示质量矩阵的行向量个数(自由度的个数)h表示其列向量个数
n=time/dt; %计算时间步长
u=zeros(m,floor(n)+3); %定义中间位移矩阵
v=zeros(m,floor(n)+2); %定义中间速度矩阵
ac=zeros(m,floor(n)+2); %定义中间加速度矩阵
U=zeros(m,floor(n)+1); %定义位移输出矩阵
V=zeros(m,floor(n)+1); %定义速度输出矩阵
AC=zeros(m,floor(n)+1); %定义加速度输出矩阵
Ke=zeros(m,1); %定义系数列向量Ke,且Ke=M/dt^2+C/2*dt
a=zeros(m,1); %定义系数列向量a,且a=K-2M/dt^2
b=zeros(m,1); %定义系数列向量b,且b=M/dt^2-C/2*dt
Pe=zeros(m,1); %定义系数列向量Pe,且Pe=-a*u(i)-b*u(i-1)
for j=1:1:m; %控制矩阵的行数
u(j,2)=u0(1,j); %将所有初始位移赋予位移矩阵第二列
v(j,2)=v0(1,j); %将所有初始速度赋予速度矩阵第二列
Ke(j,1)=M(j,j)/(dt^2)+C(j,j)/(2*dt); %计算列向量Ke
a(j,1)=K(j,j)-2*M(j,j)/dt^2; %计算列向量a
b(j,1)=M(j,j)/(dt^2)-C(j,j)/(2*dt); %计算列向量b
end
for g=1:1:m; %控制矩阵的行数,自由度的个数
for i=3:1:floor(n)+3; %控制矩阵的列数
M=input('输入质量矩阵M=');
C=input('输入阻尼矩阵C=');
K=input('输入刚度矩阵K=');
P=input('输入外部载荷矩阵P=');
u0=input('输入初始位移矩阵u0=');
v0=input('输入初始速度矩阵v0=');
time=input('输入模拟时间time=');
dt=input('输入时间步长dt=');
[m,h]=size(M); %m表示质量矩阵的行向量个数(自由度的个数)h表示其列向量个数
n=time/dt; %计算时间步长
u=zeros(m,floor(n)+3); %定义中间位移矩阵
v=zeros(m,floor(n)+2); %定义中间速度矩阵
ac=zeros(m,floor(n)+2); %定义中间加速度矩阵
U=zeros(m,floor(n)+1); %定义位移输出矩阵
V=zeros(m,floor(n)+1); %定义速度输出矩阵
AC=zeros(m,floor(n)+1); %定义加速度输出矩阵
Ke=zeros(m,1); %定义系数列向量Ke,且Ke=M/dt^2+C/2*dt
a=zeros(m,1); %定义系数列向量a,且a=K-2M/dt^2
b=zeros(m,1); %定义系数列向量b,且b=M/dt^2-C/2*dt
Pe=zeros(m,1); %定义系数列向量Pe,且Pe=-a*u(i)-b*u(i-1)
for j=1:1:m; %控制矩阵的行数
u(j,2)=u0(1,j); %将所有初始位移赋予位移矩阵第二列
v(j,2)=v0(1,j); %将所有初始速度赋予速度矩阵第二列
Ke(j,1)=M(j,j)/(dt^2)+C(j,j)/(2*dt); %计算列向量Ke
a(j,1)=K(j,j)-2*M(j,j)/dt^2; %计算列向量a
b(j,1)=M(j,j)/(dt^2)-C(j,j)/(2*dt); %计算列向量b
end
for g=1:1:m; %控制矩阵的行数,自由度的个数
for i=3:1:floor(n)+3; %控制矩阵的列数
资源评论
- u0111788902013-06-23似乎没有外部激励只是一个有初值的自由振动?
- u0105703282013-05-05资源还不错,挺有用的。
- haohaoxuexiliangjian2013-01-12资源很好,对我很有帮助
- dongyang39352016-12-12注释很详细,但是需要先看一下原理
- bej4062014-06-29资源很好,对我很有帮助
zeng3336
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功