| 查看: 1416 | 回复: 3 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
Newmark-β法计算转子的固有频率有什么优点？已有1人参与
|
||||
| Newmark-β法计算转子的固有频率有什么优点？ @kerven |
回复此楼» 猜你喜欢
各位未来的青A大佬，帮忙看看不中评审意见，是啥原因5B吗
已经有7人回复
-
靶材制备加工
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有269人回复
-
哈工大深圳杰青课题组诚招2026土木/材料/力学/海洋等专业博士生（2025.9.28日截止）
已经有12人回复
-
求aspen adsorption软件安装包破解版
已经有0人回复
-
重庆工商大学-废油资源化技术与装备教育部工程研究中心2026年硕士研究生招生
已经有41人回复
-
二苯甲酮合成
已经有1人回复
-
求教邻氯酰氯在氯化锌催化剂作用下与对氯苯胺反应能得到什么产物
已经有1人回复
-
求流动传热博导
已经有0人回复
-
硅溶胶和硅酸钠水玻璃干燥后的结构为什么出现明显差异？
已经有6人回复
-
犹他大学机械工程 Tianli Feng课题组博士生招生
已经有0人回复
kingjinjing
专家顾问 (著名写手)
-
专家经验: +159 - 机械EPI: 3
- 应助: 308 (大学生)
- 贵宾: 0.011
- 金币: 3935.9
- 散金: 5984
- 红花: 27
- 帖子: 1648
- 在线: 875.4小时
- 虫号: 1251690
- 注册: 2011-03-31
- 性别: GG
- 专业: 机械动力学
- 管辖: 机械
【答案】应助回帖
感谢参与，应助指数 +1
|
特征值计算一般采用BLOCK lanczos法。。 而newmark是时间积分算法。。。不一样的。。。 |
3楼2016-05-11 10:07:53
yangxiaoniao
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2677.3
- 散金: 1000
- 红花: 1
- 帖子: 404
- 在线: 101.8小时
- 虫号: 2667604
- 注册: 2013-09-20
- 性别: GG
- 专业: 机械动力学
2楼2016-05-10 22:28:35
送红花一朵|
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