声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 6182|回复: 34

[非线性振动] 齿轮参数振动问题(单自由度)

  [复制链接]
发表于 2011-4-1 10:31 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
1.gif
其中的无穷大,编程时我取了1000
function dx=eq9(t,x)
    dx=zeros(2,1);
    dx(1)=x(2);   
    a=0;
    for n=1:1:1000
        b=(sin(1.88*n)/n)*cos(0.0175*n*t);
        a=a+b;
    end
    c=1+0.254*a;
    dx(2)=0.74-3.3*10^(-4)*sqrt(c)*x(2)-c*x(1);

tspan=linspace(700,1800,10000);
[t,x]=ode45('eq9',tspan,[0 0]);
plot(t,x(:,1),'-');

得到的响应怎么会是这样的呢?求大家帮看看程序是不是有问题啊?谢谢。其中方程做过无量纲处理。啮合刚度是方波。
untitled.jpg
下面这张图是我把阻尼取得比较大的时候,得到的响应。

untitled1.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2011-4-1 10:48 | 显示全部楼层
在线跪求解答。
 楼主| 发表于 2011-4-1 10:50 | 显示全部楼层
跪求解答
发表于 2011-4-1 11:05 | 显示全部楼层
本帖最后由 meiyongyuandeze 于 2011-4-1 11:07 编辑

我运行了下你的程序,感觉没什么直接的错误。我没做过齿轮的研究,但仅从你的图上可以看出系统的平衡点的位置在变化,看看你能不能从物理上解释的通。
 楼主| 发表于 2011-4-1 11:07 | 显示全部楼层
回复 4 # meiyongyuandeze 的帖子

请您帮我看看程序有没有错,好吧?谢谢了
发表于 2011-4-1 11:08 | 显示全部楼层
我还在运行,在看,请问你的方程系数的量级差别大吗?
 楼主| 发表于 2011-4-1 11:13 | 显示全部楼层
方程在上面,应该不大吧。阻尼项可能比较小,但也不会影响吧
 楼主| 发表于 2011-4-1 15:41 | 显示全部楼层
回复 6 # meiyongyuandeze 的帖子

方程在上面,应该不大吧。阻尼项可能比较小,但也不会影响吧
 楼主| 发表于 2011-4-1 16:40 | 显示全部楼层
大家帮我看看程序有没有错啊,谢谢啦
 楼主| 发表于 2011-4-1 19:30 | 显示全部楼层
 楼主| 发表于 2011-4-1 19:38 | 显示全部楼层
都过节去了啊
 楼主| 发表于 2011-4-2 09:47 | 显示全部楼层
求帮助
发表于 2011-4-2 10:02 | 显示全部楼层
程序感觉没有问题,如果是做的齿轮系统的模型,那么你的参数与实际的齿轮参数是怎么样的近似过程,这个问题搞清楚了你就会感觉没有什么错误的。
发表于 2011-4-2 12:00 | 显示全部楼层
今天有重新帮你运行了下,感觉程序应该是没错误,出现这种结果可能是方程系数的问题,可能是你本身的这个动力学方程可能有些许问题吧!
 楼主| 发表于 2011-4-3 16:46 | 显示全部楼层
回复 13 # 无水1324 的帖子

谢谢! 请问我这个方程怎么求加速度对时间的曲线啊 即:加速度-时间历程曲线
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-24 07:45 , Processed in 0.099536 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表