询问求解振动方程的计算方法

toes · 发表于 2006-3-6 15:08

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

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

x

本帖最后由 VibInfo 于 2016-4-21 15:15 编辑

　　本人现在要解一个4自由度的二阶振动方程组，想得到系统的响应，不知道各位有什么好的数值求解方法?如果有计算机程序请不要吝啬，嘿嘿^_^
谢谢

欧阳中华 · 发表于 2006-3-6 21:01

本帖最后由 VibInfo 于 2016-4-21 15:15 编辑

如果忽略阻尼，可以表示为系统质量矩阵和刚度矩阵，然后找个广义JACOBI子程序就出来了...[em05][em05][em05]

linqus · 发表于 2006-3-7 11:06

本帖最后由 VibInfo 于 2016-4-21 15:15 编辑

　　(1)写出状态方程，利用odes的数值分析理论求解;

　　(2)利用newmark假设，直接对二阶微分方程进行离散，求解。

toes · 发表于 2006-3-7 13:08

本帖最后由 VibInfo 于 2016-4-21 15:15 编辑

　　是有阻尼的。

　　我用matlab的ode45和自己编写的newmark法程序算了一个简单的2自由度2阶的常微分方程组，但从频域上看得到的位移响应中有刚体位移或者有一个定频率的响应，请问这是什么原因?

　　[em06]

qinyueying · 发表于 2006-3-7 21:22

有刚体位移说明你的刚度矩阵没有加上边界约束，你说有一个定频率的响应，不明白你是什么意思，可否说得详细点

henray · 发表于 2006-3-9 15:51

本帖最后由 VibInfo 于 2016-4-21 15:15 编辑

计算时发觉该“定频率”的响应与阻尼有关，可能是由于计算开始时的不稳定造成的，加大阻尼或增加采样时间应该是有用的。
不知各位有没有碰到过这样的情况？

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

[结构振动] 询问求解振动方程的计算方法

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