| 查看: 1598 | 回复: 3 | ||||
[求助]
Newmark-β法计算转子的固有频率有什么优点？ 已有1人参与
|
| Newmark-β法计算转子的固有频率有什么优点？ @kerven |
回复此楼» 猜你喜欢
一志愿齐工大 085501机械工程 304分求调剂
已经有2人回复
-
085501机械工程304分求调剂，一志愿齐工大
已经有15人回复
-
金属材料论文润色/翻译怎么收费?
已经有106人回复
-
304分085501机械工程 一志愿齐工大
已经有12人回复
-
求调剂304分085501机械工程 一志愿齐工大
已经有15人回复
-
西北工业大学机电学院 复合材料加工制造方向招聘博士后
已经有12人回复
-
大连海事大学轮机工程学院尚有少量调剂名额，如有兴趣的同学可以联系。
已经有33人回复
-
湖北汽车工业大学，来啦！汽车材料学院2026年招收调剂研究生计划
已经有0人回复
-
【模拟经典书】《Handbook of Materials Modeling》2020年 第2版
已经有5人回复
-
国自然申报交流祈福
已经有153人回复
-
哈尔滨理工大学材料加工工程专业26秋入学博士尚有少量招生名额
已经有0人回复
yangxiaoniao
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2929.7
- 散金: 1000
- 红花: 1
- 帖子: 455
- 在线: 102小时
- 虫号: 2667604
- 注册: 2013-09-20
- 性别: GG
- 专业: 机械动力学
2楼2016-05-10 22:28:35
kingjinjing
专家顾问 (著名写手)
-
专家经验: +159 - 机械EPI: 3
- 应助: 308 (大学生)
- 贵宾: 0.011
- 金币: 3935.9
- 散金: 5984
- 红花: 27
- 帖子: 1648
- 在线: 875.4小时
- 虫号: 1251690
- 注册: 2011-03-31
- 性别: GG
- 专业: 机械动力学
- 管辖: 机械
【答案】应助回帖
感谢参与，应助指数 +1
|
特征值计算一般采用BLOCK lanczos法。。 而newmark是时间积分算法。。。不一样的。。。 |
» 本帖已获得的红花（最新10朵）
3楼2016-05-11 10:07:53
送红花一朵 |
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。我只知道里面可以得到固有频率 |
4楼2016-05-12 14:33:03
20