| 查看: 1612 | 回复: 3 | |||
| 【悬赏金币】回答本帖问题，作者霍笔小新将赠送您 5 个金币 | |||
[求助]
用求特征值的方法求无阻尼系统的固有频率时为复数，要怎么办
|
|||
|
用求特征值的方法求无阻尼系统的固有频率时一阶固有频率为复数，怎么取呢，取虚部的值吗。还是说我的刚度矩阵错了。程序如下 clear;clc m=3; z=22; theta0=pi/(2*z);%半齿厚角 rb=47.631/2000;%基圆半径 km=1.5972e6;%啮合刚度（单位N.M） kp1=0.5e7;%0.5e7;%支承刚度,单位N / m kp2=0.5e7;%0.5e7;%支承刚度,单位N / m kt1=7.855527131099762e+04;%0.254e8; kt2=7.750888522813271e+04;%0.254e8;%扭转刚度及阻尼 J1=0.00343577787; %转动惯量 J2=0.00373206031; Jl=75.263e-4; m1=4.51062; m2=4.55723; x=1:7; for k=1:m thetai(k)=theta0+2*pi*(k-1)/z-pi/6;%+wn*t; l1(k)=(sin(thetai(k))).^2; l2(k)=sin(thetai(k)).*cos(thetai(k)); l3(k)=(cos(thetai(k))).^2; l4(k)=sin(thetai(k)); l5(k)=cos(thetai(k)); end ll1=sum(l1); ll2=sum(l2); ll3=sum(l3); ll4=sum(l4); ll5=sum(l5); M = diag([m1 m1 m2 m2 J1 J2 Jl]); K = [kp1+km*ll1, -km*ll2, -km*ll1, km*ll2, -rb*km*ll4, rb*km*ll4, 0; -km*ll2, kp1+km*ll3, km*ll2, -km*ll3, rb*km*ll5, -rb*km*ll5, 0; ; -km*ll1, km*ll2, kp2+km*ll1, -km*ll2, rb*km*ll4, -rb*km*ll4, 0; km*ll2, -km*ll3, -km*ll2, kp2+km*ll3, -rb*km*ll5, rb*km*ll5, 0; -rb*km*ll4, rb*km*ll5, rb*km*ll4, -rb*km*ll5, kt1+rb^2*km, -rb^2*km, 0; rb*km*ll4, -rb*km*ll5, -rb*km*ll4, rb*km*ll5, -rb^2*km, kt2+rb^2*km, -kt2; 0, 0, 0, 0, 0, -kt2, kt2];%刚度矩阵 [V,D] = eig(K,M); % 求解特性值和特征向量 freqs = sqrt(diag(D))% 固有圆频率 modes = V; % 模态向量（振型） 计算结果如下 freqs = 1.0e+03 * 0.0000 + 0.2695i 1.0500 + 0.0000i 1.0501 + 0.0000i 1.1023 + 0.0000i 1.7730 + 0.0000i 4.8165 + 0.0000i 5.5921 + 0.0000i |
回复此楼» 猜你喜欢
pulverisette 7 planetary micro mill 高能球磨金属颗粒团聚和粘粉问题
已经有13人回复
-
欢迎参加第二届“X 射线晶体学（桂林）基础讲习班 及应用研讨会”（5.20-5.24）
已经有6人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有77人回复
-
欢迎参加第二届“X 射线晶体学（桂林）基础讲习班 及应用研讨会”（5.20-5.24）
已经有21人回复
-
第二届“X 射线晶体学（桂林）基础讲习班 及应用研讨会” (第二轮通知)
已经有10人回复
-
2025申博——配位聚合物（MOF）荧光性能研究方向
已经有16人回复
-
vasp计算能带报错，出现BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
已经有1人回复
-
求助：想看一下电子束蒸发系统（OHMIKER 50b）的样品腔图片
已经有0人回复
-
喷涂铜锰尖晶石求助
已经有1人回复
-
求助查看晶向和晶面的软件
已经有0人回复
-
设计了一套越障小车，不会设计电路控制系统。哪位有推荐的电路系统，求分享
已经有0人回复
曾经落伍了
铁杆木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 16133.8
- 散金: 2541
- 红花: 6
- 帖子: 3286
- 在线: 570小时
- 虫号: 4078713
- 注册: 2015-09-17
- 性别: GG
- 专业: 机械动力学
2楼2020-02-13 00:28:11
曾经落伍了
铁杆木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 16133.8
- 散金: 2541
- 红花: 6
- 帖子: 3286
- 在线: 570小时
- 虫号: 4078713
- 注册: 2015-09-17
- 性别: GG
- 专业: 机械动力学
3楼2020-02-13 00:28:50
4楼2020-02-13 08:44:52
30