zhijian19840707@163.com 这个关键是如何形成M,C,K,J的问题
简易参考钟一鄂的《转子动力学》第四、七、八章内容
分别是三种不同的方法,有详细的形成过程
然后用matlab的eig求解就行了
MATLAB程序
能不能给我一份MATLAB程序,最近正在研究!!还有没有利用matlab做振型和不平衡相应以及瞬态响应的??谢谢啦!!急用!!邮箱lovegirlmeijie@163.com 我最近在学机械动力学,也想要一份。谢谢!LOVERAPING@EYOU.COM 呵呵,找了好几年了。朋友能否给我一份?jinghuixie@hotmail.com回复 #3 21172485 的帖子
在作毕业设计,我也想要一份,谢谢了,gongchenglixuehf@163.com 自己动手吧,有半天时间就能搞定,而且还是考虑比较全面的程序 本帖最后由 VibInfo 于 2016-4-22 15:13 编辑原帖由 bysj_chao 于 2006-4-17 13:49 发表
不知道有没有哪位做转子临界转速分析的!希望提供点MATLAB方面的程序!谢谢?帮忙呀!
%转子系统前三阶临界转速法
clc
clear
%等截面轴参数
l1=0.12;
d=0.04;
A=pi*d*d/4;
%轮盘参数
D=0.5;
h=0.025;
%盘轴材料参数(忽略轴的质量)
a=1;
u=0.3;
rou=7800;
E=2.0e11;
G=E/(2*(1+u));
I=pi*(d^4)/64;
K1=2.0e7;
v1=6*E*I/(a*G*A*l1*l1);
mi=rou*pi*D^2/4*h;%轮盘的集质量
Jp=mi*D^2/8;
Jd=Jp/2;
Ji=Jp-Jd;
%参数的数组形式
L=;
M=;
K=;
v=;
J=;
k=0;
Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图';'第三阶频率的振型和弯矩图'];
for w=0:0.01:4000;
for i=1:15;
T(:,:,i)=
(L(i)^2)*(M(i)*w^2-K(i))/(2*E*I) 1+L(i)*J(i)*w^2/(E*I) L(i)/(E*I) L(i)^2/(2*E*I);
L(i)*(M(i)*w^2-K(i)) J(i)*w^2 1 L(i);
M(i)*w^2-K(i) 0 0 1];
end
H=T(:,:,1);
for i2=2:15;
H=T(:,:,i2)*H;
end
F=H(3,1)*H(4,2)-H(3,2)*H(4,1);
if F*(-1)^k < 0 %求解临界转速
k=k+1;
wi(k)=w;
w=wi(k)
ni(k)=wi(k)*30/pi;
end
end 谁有用传递矩阵求转子临界转速的MATLAB程序 ?
我想学习一下,谢谢!邮箱yangfeiyxl@tom.com 麻烦,也给我发一份,我真在为此头疼中。
mememuluan@163.com ,提前先谢谢了阿
我的邮箱是 gongchenglixuehf@163.com
我想要那个用vc编的程序,谢谢了 我也要一份matlab的,zhenruiwang@gmail.com谢谢了分享一下
分享一下matlab的传递矩阵程序. email:luxianxiao315@163.com 给我也发一份吧,谢谢! tianyarrrbbb@163.com程序
我也要一份,谢谢了286176658@qq.com