求：二自由度振动方程的解法

wanyeqing2003

原帖由 咕噜噜 于 2007-6-24 20:06 发表
其是我觉得有时候不必要过分强调字母的含义，抱歉，虽然理论上大家都认为k是刚度，c是阻尼，不过这里我是把它倒过来看了，这是一般方程的理解问题，不要太拘泥于字面的东西了

是的。我也倾向这样的解释。只是前面有网友提出此方程不好解耦，使我偶发联想，随口答来。

至于最近没有到访小郝的专栏，实在抱歉。最近由点忙，出差去了一趟济南，再加上工作比较多。以后有时间，我会到你的栏目拜访，因为我的分叉问题只研究了一半，我会继续做完它的。到时候少不了请你帮忙。
谢谢你的邀请。

joeddong

谢谢
wanyeqing2003
无水1324
VibrationMaster
咕噜噜
drbethune （排名不分先后）
等
诸位帮忙，小弟看了2天书，传几个版本的解法上来，帮忙再看一下还有问题不？
多谢了！:victory: :victory: :victory:

[ 本帖最后由 joeddong 于 2007-6-25 09:12 编辑 ]

joeddong

用这种方法求解的时候出现了一个问题
就是用复数 代替激振力
的时候，用的是欧拉方程中的虚部，
再最后方程求解后的反带过程是不是也一样
取复数的虚部就可以了呢？
这样解可以不？


[ 本帖最后由 joeddong 于 2007-6-25 10:02 编辑 ]

joeddong

这是化为单自由度求解过程，
出现了一点问题，即
在解方程（10）的时候，
得到的结果太复杂，不实用，操作性差。

[ 本帖最后由 joeddong 于 2007-6-25 10:02 编辑 ]

咕噜噜

大致看了一下，你这个方法是有局限性的，你这么做出来得可是通解，不是特解，不过也可以看出长期效果，但是我不知是不是一定可行，但是从你的算法上来说问题不是很大，有时间我在看一下，顺便帮你找个其他思路

wanyeqing2003

是的。用矩阵方法做比较麻烦。不过用matlab是可以实现的。

如果用振型分解法就会简单一些。

我想有阻尼受迫振动，特解应该是可以的。

joeddong

原帖由 咕噜噜 于 2007-6-25 11:59 发表
大致看了一下，你这个方法是有局限性的，你这么做出来得可是通解，不是特解，不过也可以看出长期效果，但是我不知是不是一定可行，但是从你的算法上来说问题不是很大，有时间我在看一下，顺便帮你找个其他思路

多谢，就是觉得太麻烦了!
不具有操作性！

joeddong

原帖由 wanyeqing2003 于 2007-6-25 12:07 发表
是的。用矩阵方法做比较麻烦。不过用matlab是可以实现的。

如果用振型分解法就会简单一些。

我想有阻尼受迫振动，特解应该是可以的。

对振型分解法不是很熟悉，正在看书！
多谢了哦，不懂的再向你请教哦！

wanyeqing2003

原帖由 joeddong 于 2007-6-25 12:29 发表


对振型分解法不是很熟悉，正在看书！
多谢了哦，不懂的再向你请教哦！

振型分解法是先求特征值，再由特征值计算特征向量，然后根据特征向量作正交变换，也就是解耦。
在这样的正交坐标下，计算振动响应。最后再通过坐标变换（其实是反变换）计算出实际坐标下的物理量，也就是所要计算的结果。

咕噜噜

我已经告诉你那本参考书了啊，找找看，那本书上就有你这个问题，很简单
因为说实话我现在事情也多，不过你可能自己看的话我还是鼓励你自己解决，虽然你有点着急，不过这个的确不需要多少时间就能看完，另你要是有系统动力学的基础那几乎就不用半天就可以，matlab编程求特征值什么的有问题可以到时候再上来讨论一下

long773910

dddd