troumar 发表于 2012-5-12 15:43

共振情况下，正弦荷载作用下的位移失衡的数值解和理论解不吻合！

数值解我是用频响函数Hw与荷载时程的傅里叶变换后的频域值Pf的乘积得到的，即yf=Hw*Pf，再对该位移的频域值进行一次傅里叶逆变换得到位移时程yt，其中，频响函数Hw=(1/(1-(w/w1)^2+i*(2*c*(w/w1))))/k;其中，w和w1分别指荷载频率和结构自振频率，c阻尼比，k指结构刚度，i是虚数。理论解就是用的动力学中的稳态解，当荷载频率与结构自振频率不相等时，位移时程的数值解和理论解曲线都能吻合在一起，就只有当荷载频率和自振频率相等时，就出现了图中所示的情况，即便荷载频率与自振频率相差很小时都不会有下图的情形。希望各位高手帮帮忙，指出一下问题，不胜感激！http://fmn.rrimg.com/fmn056/20120511/1045/b_large_E0wt_668a00006cb61263.jpg当荷载频率和自振频率相差0.01的时候曲线图如下，该种情况，数值解与理论解吻合的非常好。http://fmn.rrimg.com/fmn064/20120511/1515/b_large_6AAp_3e8400009ee31261.jpg

dw04116 发表于 2012-5-14 09:20

回复 1 # troumar 的帖子

我看的不仔细，问个简单的，你这里的w1和k，m什么关系

troumar 发表于 2012-5-16 22:21

回复 2 # dw04116 的帖子

w1=(k/m)½

dw04116 发表于 2012-5-17 09:17

回复 3 # troumar 的帖子

这种“拍”现象貌似两自由度，固有频率很接近时才会发生。
是不是你的载荷就是两个频率很接近的简谐激励力呢？
（顺便请教个问题，你的傅里叶变换和逆变换如何在matlab中实现的）

dw04116 发表于 2012-5-17 09:27

回复 1 # troumar 的帖子

还有一种可能，激励频率接近固有频率，但是又不等于固有频率时，参考《机械振动与噪声学》P98，书中说的是无阻尼，但是我觉得有阻尼可能有类似的结论

troumar 发表于 2012-6-14 16:22

回复 4 # dw04116 的帖子

不好意思，很久没来了，呃，傅里叶的正逆变换我用的是fft和ifft命令

tiger0004 发表于 2012-7-5 15:43

好明显的拍现象。。。
页: [1]
查看完整版本: 共振情况下，正弦荷载作用下的位移失衡的数值解和理论解不吻合！