24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

2267981137

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 107
帖子: 30
在线: 2.3小时
虫号: 4585241
注册: 2016-04-13
专业: 机械动力学

[求助] Newmark-β法计算转子的固有频率有什么优点？已有1人参与

Newmark-β法计算转子的固有频率有什么优点？ @kerven

回复此楼

» 猜你喜欢

招收硕，博士，博士后及联培：设备状态监测，可靠性与维护，人工智能，图像处理已经有0人回复
招收26年秋季入学博士生（北科大高精尖学院力学超材料/机器学习/增材制造相关方向）已经有18人回复
金属材料论文润色/翻译怎么收费? 已经有70人回复
国内权威又专业的知名论坛有哪些？已经有0人回复
有没有大佬推荐一下审稿快，稿费在1-2w之间的EI综合期刊，孩子急着毕业啊已经有10人回复
建筑、汽车、家具涂料如何合规？2025年四大权威第三方检测机构对比已经有0人回复

1楼2016-05-10 19:28:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

2267981137

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 107
帖子: 30
在线: 2.3小时
虫号: 4585241
注册: 2016-04-13
专业: 机械动力学

送红花一朵

引用回帖:

3楼: Originally posted by kingjinjing at 2016-05-11 10:07:53
特征值计算一般采用BLOCK lanczos法。。
而newmark是时间积分算法。。。不一样的。。。

E=210*10^9; %转子弹性模量
L=0.45; % 转轴长度
m1=0.8; % 转子质量
R=0.005; % 转轴半径
e=0.03; % 偏心块的偏心距
m2=0.002; %　偏心块质量
em1=m2*e/m1;
g1=9.8;
w=300; % 转速
dt=pi/500; %时域时间间隔
t0=0;
I=pi*R^4/4; %轴截面转动惯量
kz=48*E*I/L^3; %轴ｚ方向刚度
wn=(kz/m1)^0.5; %转子固有频率
q=m1*g1/kz; %转子的静位移

U=em1/q;
W=w/wn; %U和W是无量纲值

%NewmarkBeta算法
r=0.5;
a=0.25;
a1=1/(a*dt^2);
a2=-1/(a*dt);
a3=-(0.5/a-1);
b1=r/(a*dt);
b2=1-r/a;
b3=(1-0.5*r/a)*dt;

%无量纲化后各系数矩阵
m=eye(2); %无量纲化后，质量矩阵为单位矩阵
c=2*[0.01/W 0;0 0.01/W];%阻尼比取０．０１
k=[1/W^2 0;0 1/W^2];
g=[U*sin(t0);U*cos(t0)];

%初始位移、速度、加速度
u0=[0;0];
up0=[0.001;0.001];
%upp0=m^(-1)*(g-k*u0-c*up0);
upp0=g-k*u0-c*up0;

u=u0;
up=up0;
upp=upp0;
fid=fopen('zhuanzi.txt','w');
for t=(t0+dt):dt:300*pi
g=[U*sin(t);U*cos(t)];
K=a1*m+b1*c+k;
G=g+m*(a1*u-a2*up-a3*upp)+c*(b1*u-b2*up-b3*upp);
u1=K^(-1)*G;
up1=b1*(u1-u)+b2*up+b3*upp;
upp1=a1*(u1-u)+a2*up+a3*upp;
u=u1;
up=up1;
upp=upp1;
fprintf(fid,'%f %f %f\r\n',u1(1),u1(2),t); %'\r\n'数据换行数据分为3列
end
fclose(fid);

data=importdata('zhuanzi.txt');%%将zhuanzi.txt放到matlab的工作目录下，运行一下代码即可
data_x=data(:,1);
data_y=data(:,2);
data_t=data(:,3);
%mean_y=mean(data_y);%
%disp(['第三列，即Y坐标值的均值为' num2str(mean_y)])
x=q.*data(:,1) %x方向的实际振幅
y=q.*data(:,2) %y方向的实际振幅
figure(1)
subplot(2,1,1)
plot(data_t(100000:150000),data_x(100000:150000));%取时域的一部分,从100000个时间点开始到最后
grid
xlabel('t')
ylabel('x')
title('x方向的理论振幅');
subplot(2,1,2)
plot(data_t(100000:150000),x(100000:150000));%取时域的一部分,从100000个时间点开始到最后
grid
xlabel('t')
ylabel('x')
title('x方向的实际振幅');
figure(2)
plot(data_x(100000:150000),data_y(100000:150000));
xlabel('x');
ylabel('y');
title('轴心轨迹图');
能麻烦您看一下，我这个程序是计算什么的。这里面用到了Newmark。我只知道里面可以得到固有频率