数值解法与幅频特性分析

siyanger · 发表于 2006-8-19 00:52

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

x

用数值法解出多自由度振动的响应，因为是离散的点，没有近似表达式，怎么能做出非线性的幅频曲线呢？

siyanger · 发表于 2006-8-19 00:56

如上图（其中，2.0应为1.0)。因为用数值法，解出来一个频率值对应一个幅值，怎么求得多值？还有跳跃。

siyanger · 发表于 2006-8-22 09:06

没有人做过类似的研究吗？谢谢啦！

无水1324 · 发表于 2006-8-29 12:38

同上问题!
高手指教!

siyanger · 发表于 2006-8-29 13:27

你做得怎么样了？多值就象上次你说的，不同的初值，解会不同的。但线性的话，不同的初值解是相同的。这个我试过了。跳跃是指解的值突然变大吧。

siyanger · 发表于 2006-8-29 13:31

还有，我做出来的解，想画个poincare图，按论坛上给的程序，怎么画都不对(因为线性非线性结果一样），你有没有做poincare图？

无水1324 · 发表于 2006-8-29 18:05

跳跃（见1楼图）应该是横坐标变化时，振幅由一个值，变成两个不同的值等的一个跳跃过程

poincare图我就是画的频闪图

siyanger · 发表于 2006-8-30 15:37

不好意思，频闪图怎么画啊？按周期取点？我按周期取点，纵横坐标值都一样大小，没有变化了。

无水1324 · 发表于 2006-8-30 17:12

间隔一个周期取点
当分岔参数r(i)时，位移(也可以是速度项)每过一个周期取一个点x(n1:T:n2,1)
然后plot(r, x)

无水1324 · 发表于 2006-8-30 17:13

还有你可以去图书馆借一个书《理论力学的计算机模拟》好像是这个
那里面有闪频的程序
我这里刚搞没有那个书

siyanger · 发表于 2006-8-31 20:29

原帖由 无水1324 于 2006-8-30 17:12 发表
间隔一个周期取点
当分岔参数r(i)时，位移(也可以是速度项)每过一个周期取一个点x(n1:T:n2,1)
然后plot(r, x)

我看一般poincare图都取位移和速度为坐标的。分岔参数一般取为激励频率，plot(r, x)也行吗？

无水1324 · 发表于 2006-9-1 14:20

不好意思说成分岔图了

闪频就是位移速度分别间隔一周期取值然后画图就可以了

siyanger · 发表于 2006-9-1 19:16

比如周期为T,这样画图:
t=t1:T:t2 ( t1 t2 为积分区间)
plot(x,dx)
我这样画出来好象不对．　

无水1324 · 发表于 2006-9-1 20:48

应该是
plot(x(n1:n:n2,1),x(n1:n:n2,2))
x(n1:n:n2,1)，为位移
x(n1:n:n2,2),速度
n为一个周期内所求的点
n1为画图的起点
n2 为终点

siyanger · 发表于 2006-9-1 23:02

原帖由 无水1324 于 2006-9-1 20:48 发表
应该是
plot(x(n1:n:n2,1),x(n1:n:n2,2))
x(n1:n:n2,1)，为位移
x(n1:n:n2,2),速度
n为一个周期内所求的点
n1为画图的起点
n2 为终点

还要请教一下,n为一个周期内所求的点怎么算啊?我试过了,n不同,画出来的图也不相同呢.多谢啦!

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

[动力学和稳定性] 数值解法与幅频特性分析