| 查看: 3598 | 回复: 12 | ||
[求助]
非线性问题，扭转振动方程matlab求解 已有1人参与
|
||
|
用matlab求解扭转振动方程，求的扭转角为什么不对呢？是程序的问题？还是非线性没处理好呢？ 方程，参数见图，还有程序，求指导！金币就3个，以后有了可以再补！谢谢  程序如下： %%%齿轮-转子系统扭转振动分析 单位 Kg m s function dx = niu_gou_yuan_Rg(t,x) TR = 30; TL = 6;%输入扭矩TR、输出扭矩TL I1 = 2.8; I2 = 2.6; I3 = 2.6; I4 = 2.6;%转动惯量 a= pi*20/180;%啮合角a m=3;%模数m Z2 = 14;Z3 = 14;% 齿数 Rb2 = (Z2*m*cos(a))*(1/2)/1000 ;%主动齿轮2的基圆半径,基圆db=d*cosa, Rb3 = (Z3*m*cos(a))*(1/2)/1000 ;%从动齿轮3的基圆半径 K1 = 4e5 ; C1 = 69;%扭转刚度及相应阻尼 K3 = 4e5 ; C3 = 63; %齿轮啮合误差激励ep ep=0.000315*sin(0.542*t-2*pi/9)+0.00002*sin(6.742*t-2*pi/9)-8.8003e-04 ; %齿轮啮合误差激励ep的导数 dep=1.7073e-04*cos(0.542*t-2*pi/9)+1.3484e-04*cos(6.742*t-2*pi/9) ; %齿轮间沿啮合线方向的弹性变形 xp = Rb2*x(3)-Rb3*x(5)-ep; dxp = Rb2*x(4)-Rb3*x(6)-dep; %啮合刚度K2n,啮合平均刚度Kp2,啮合波动刚度系数Kb2 %%齿轮副啮合频率we, 啮合相位角uk we=18; uk=0;Kb2 =1.2e8; Kp2 = 1.2e8; K2n= Kp2+Kb2*cos(we*t+uk); %啮合阻尼C2n C2n=137; b = 0.1e-3;%间隙 %啮合力 if xp > b Fp = K2n*(xp-b); elseif -b <= xp <= b; Fp = 0; else Fp = K2n*(xp+b); end %阻尼误差激励 Dp = C2n*dxp; F = Fp+ Dp; dx =[ x(2); (TR-K1*(x(1)-x(3))-C1*(x(2)-x(4)))/I1; x(4); (-Rb2*F-K1*(x(3)-x(1))-C1*(x(4)-x(2)))/I2; x(6); (Rb3*F-K3*(x(5)-x(7))-C3*(x(6)-x(8)))/I3; x(8); (-TL-K3*(x(7)-x(5))-C3*(x(8)-x(6)))/I4; ]; 求解程序： clc; x0=[0,0,0,0,0,0,0,0]; [t,y]=ode45(@niu_gou_yuan_Rg,[0,0.5],x0); figure plot(t,y(:,1)) xlabel(\'时间/s\'),ylabel(\'扭转角/rad\'),title(\'电机扭转角\') figure plot(t,y(:,3)) xlabel(\'时间/s\'),ylabel(\'扭转角/rad\'),title(\'主动齿轮扭转角\')  方程.png  模型.png  扭转角.png |
回复此楼» 猜你喜欢
【博士招生】南京林业大学车辆工程课题组招收26年9月入学“申请-考核制”博士研究生
已经有16人回复
-
求收留
已经有8人回复
-
金属材料论文润色/翻译怎么收费?
已经有249人回复
-
OLGA软件清管模拟 为什么老是清管方向相反 怎么修改
已经有0人回复
-
燕山大学纳米能源中心黄建宇课题组招生
已经有0人回复
-
中国地质大学（北京）深部探测与成像全国重点实验室诚邀全球英才申报 2026 年海外优青
已经有0人回复
-
青岛市直事业编，95硕，180 65，老家河南，爱好吉它唱歌
已经有0人回复
-
航空工程363分求调剂
已经有8人回复
-
求助调剂推荐
已经有0人回复
-
355 山科电气专硕 调剂推荐
已经有1人回复
-
11408 369分调剂求助
已经有2人回复
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
5楼2016-03-09 19:01:27
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
假装发呆8楼2016-03-09 19:37:53
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
 本帖仅楼主可见 |
2楼2016-03-09 12:28:51
3楼2016-03-09 14:28:19
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
4楼2016-03-09 18:58:35
6楼2016-03-09 19:13:53
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
7楼2016-03-09 19:23:41
9楼2016-03-09 19:41:29
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
|
本帖仅楼主可见 |
10楼2016-03-09 19:43:47
