| 查看: 3538 | 回复: 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 |
回复此楼» 猜你喜欢
燕山大学机械工程学院Marian Wiercigroch院士团队招收2026年博士5人
已经有0人回复
-
卢森堡大学拟招收2026博士生一名，CSC渠道， 协助完成研修计划
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有173人回复
-
华北电力大学王鹏课题组诚招博士/硕士
已经有3人回复
-
东北大学2026级申请考核制博士招生（安全工程、矿业工程、通风、能源动力等方向）
已经有19人回复
-
中国地质大学（北京）工程技术学院薛启龙课题组招收2026年博士生
已经有3人回复
-
招收硕，博士，博士后及联培：设备状态监测， 可靠性与维护，人工智能，图像处理
已经有0人回复
-
粉末热压生坯代加工
已经有0人回复
-
已编辑
已经有4人回复
-
湖南大学JQ团队段辉高教授课题组招光学、算法博士啦！2026级，2个名额，等你加入！
已经有11人回复
-
推荐一本新书，《半导体材料计算》vasp软件入门学习教程，配套参数文件
已经有0人回复
匿名
用户注销 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
假装发呆8楼2016-03-09 19:37:53
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 (幼儿园)
- 金币: 1601.4
- 红花: 5
- 帖子: 125
- 在线: 386.5小时
- 虫号: 0
- 注册: 2013-03-10
- 性别: MM
- 专业: 机械测试理论与技术
5楼2016-03-09 19:01:27