求助：请大家帮我看一下

凌枫 · 发表于 2009-10-16 11:21

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

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

x

这是一个单自由度齿轮系统无冲击状态时的振动幅值问题无冲击状态时的量纲一化方程为：y'' + 2*e*u* y' + (1+2*e*cos(a*t)) * (y-1) = f0 + e* f1 * cos(a*t+b)
其中e=0.01为小参数，u=0.1为阻尼，f1=1，f0=1，b=pi/3，求系统振动幅值与频率比a的关系，当a=1，2时系统的幅值应较大，可是在用龙哥库塔法解的时候发现改变a 的值，方程的位移时间曲线没变化，请大家帮忙看一下是哪里出问题了？谢谢
function dy=wuchongji(t,y)
e=1;
u=0.1;
a=1;
b=0;
f0=0;
f1=1;
dy=[y(2);-2*e*u*y(2)-(1+2*e*cos(a*t))*(y(1)-1)+f0+e*f1*cos(a*t+b)];

[t,x]=ode45('wuchongji',[0:0.02:200],[0,1]);
plot(t,x(:,1));
就是改变a 的值，画出来的曲线却没有变化，不知咋回事？

无水1324 · 发表于 2009-10-18 21:18

a=1，2时系统的幅值应较大？
这个是怎么得到的？

凌枫 · 发表于 2009-10-20 20:43

因为a是频率比，在1、2处应该是主共振和次共振，所以幅值应较大，问题好像是积分时间短了

无水1324 · 发表于 2009-10-21 09:41

问题已经解决？

octopussheng · 发表于 2009-10-21 22:03

是否做过定量分析得出的结论？

凌枫 · 发表于 2009-10-21 22:06

没有做过定量分析，这是老板叫我审稿的一篇论文里面的，把积分时间取到5000后就可以了

凌枫 · 发表于 2009-10-21 22:12

顺便请教一个问题，ODE45里面能不能得到加速度时间响应？

无水1324 · 发表于 2009-10-27 08:37

可以
只是你根据方程再把，速度，位移代入可以得到的

凌枫 · 发表于 2009-10-27 18:32

嗯，谢谢无水指点

账号		自动登录	找回密码
密码			我要加入

[计算数学] 求助：请大家帮我看一下

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

回复板凳凌枫的帖子

回复 7楼凌枫的帖子

[计算数学] 求助：请大家帮我看一下

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

回复 板凳 凌枫 的帖子

回复 7楼 凌枫 的帖子

回复板凳凌枫的帖子

回复 7楼凌枫的帖子