| 查看: 1497 | 回复: 3 | ||||
[求助]
Newmark-β法计算转子的固有频率有什么优点？ 已有1人参与
|
| Newmark-β法计算转子的固有频率有什么优点？ @kerven |
回复此楼» 猜你喜欢
燕山大学先进结构材料团队招收金属材料相关方向 2026年硕士研究生
已经有0人回复
-
2026考研调剂上海工程技术大学材料学院
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有63人回复
-
2026上海工程技术大学李老师化学材料机械计算机方向招硕士7-10名
已经有0人回复
-
上海工程技术大学材料方向招收硕士7-10名
已经有10人回复
-
高温合金厂家
已经有0人回复
-
西安石油大学新能源学院接收材料类、能源动力类、机械类、计算机类等专业专硕调剂生！
已经有11人回复
-
西安石油大学新能源学院接收材料类、能源动力类、机械类、计算机类等专业专硕调剂生！
已经有11人回复
-
JMPT期刊有字数限制吗
已经有0人回复
-
使用高温共聚焦显微镜原位观察熔渣在钢表面的铺展情况
已经有0人回复
yangxiaoniao
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2914.7
- 散金: 1000
- 红花: 1
- 帖子: 452
- 在线: 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
