声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2312|回复: 11

[其他相关] 怎样由频散关系曲线去判断是哪一种振动模态

[复制链接]
发表于 2007-11-29 16:31 | 显示全部楼层 |阅读模式

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

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

x
我通过解方程的数值解,得到几条频散关系曲线,但不知道每条曲线是代表哪种振动模态
回复
分享到:

使用道具 举报

发表于 2007-11-30 12:22 | 显示全部楼层
概念错误,什么叫频散关系?

模态的物理意义一般通过特征值或振型来判断.
发表于 2007-11-30 13:54 | 显示全部楼层

回复 #1 tangjikede 的帖子

看了几个帖子,确实还是不知道频散是什么意思
 楼主| 发表于 2007-11-30 15:55 | 显示全部楼层
频散关系是波数和波速之间的关系。我研究的是碳纳米管中的简谐波传导,有了频散关系后,不知道该怎样去判断波在纳米管中的传导是什么样的形式,为什么取圆波数m=1时,会有很多条解曲线。m=1意味着是弯曲振动吧。但是为什么会有很多条关系曲线
发表于 2007-11-30 16:12 | 显示全部楼层
I see。 但是振动问题中认为物体尺寸非常小,因而没有波动效应。 你这种不是振动问题,而是波动力学问题。是否存在振动模态这种概念,本人不知道。
发表于 2011-5-28 00:43 | 显示全部楼层
回复 1 # tangjikede 的帖子

你好,请问您的方程是怎么解出来的 我现在要画出这个频散曲线 可是不知道怎么做,请教
发表于 2011-7-3 22:11 | 显示全部楼层
回复 1 # tangjikede 的帖子

您好! 我直接求解行列式,应该是一个表达式,结果得到的是0,请问这是为什么?请高手指教!
发表于 2011-7-4 13:22 | 显示全部楼层
回复 7 # yuandianqi 的帖子

将位移波u=U*exp[i(kx-wt)]代入波动方程,会得到关于U的方程组,其中还有k(波数)、w(圆频率)是未知数。因为方程组是齐次方程组,若有非零解,系数行列式为零,这样就得到关于k、w的一个方程。解此方程,得到k、w关系,即频散关系。
此外还可得到相速度cp=w/k,群速度cg=dw/dk

评分

1

查看全部评分

发表于 2011-8-13 09:07 | 显示全部楼层
回复 8 # sunhu2003 的帖子

谢谢,我现在直接应用的资料上的结果,即把他的行列式直接写出,然后令他为0.得到的应该是一个关于w和cp相速度的高阶方程。请问下面应该这么画出w-cp的曲线呢?我得到的行列式为disperse。他的表达式为下面的形式:比较复杂
isperse=-w^6*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2)*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2)*(.669529e-9*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))-.1798160e-5*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))-.1798160e-5*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*cp^4*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4-.1798160e-5*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2+.975975e-3*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))-.1798160e-5*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2)))/cp^4
发表于 2011-8-13 15:45 | 显示全部楼层
回复 9 # yuandianqi 的帖子

哇塞,这个很复杂呀。解析解肯定是给不出了。
试试离散的吧。取一系列w值,看每个w值对应的这个行列式方程能不能解出cp(也许有什么命令可以解出哈)。

评分

1

查看全部评分

发表于 2011-8-13 22:41 | 显示全部楼层
回复 10 # sunhu2003 的帖子

好的,我再找找命令,试试。谢谢您
发表于 2011-8-19 21:32 | 显示全部楼层
本帖最后由 yuandianqi 于 2011-8-19 21:33 编辑

回复 8 # sunhu2003 的帖子
没出来啊,愁呀。还有一个问题,即使最后出来了,得到的也是很多的点,那可怎么确定到底怎么画线,怎么判断是哪个模态啊?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-19 20:49 , Processed in 0.128656 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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