声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3262|回复: 6

[结构振动] 求助:有关结构振动控制Riccati方程的求解问题

[复制链接]
发表于 2005-10-23 17:20 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 VibInfo 于 2016-4-19 15:58 编辑

  我在做一个结构振动控制方面的课题,采用的是最优二次行控制,其实就是对结构动力学方程进行变换,变成那个dx=Ax+Bu的格式,A=[0,I; -K/M,o], 然后采用反馈控制的方式,这个问题在matlab里面似乎没有问题,我是在nastran里面作这个东西的,结果在使用Nastran里面编程后riccati方程的求解就是不收敛,不知道是身么问题。

  我做过很多中riccati方程的求解算法,都不能,所以问题应该出在第一步,初步怀疑这个A矩阵是不稳定的,不知道在此的同行们有没有对这个A有比较透彻的研究的,这个怀疑时候猜对了方向呢,如果是不稳定,那我又该怎么办呢?

  请各位在此的指点迷津,不胜感激

  我的qq号:23489744,欢迎探讨
回复
分享到:

使用道具 举报

发表于 2005-10-23 20:51 | 显示全部楼层

回复:(yun)求助:有关结构振动控制Riccati方程的求...

K和m分别是多少?
发表于 2005-10-23 20:54 | 显示全部楼层

回复:(yun)求助:有关结构振动控制Riccati方程的求...

另外,你在用MATLAB中使用什么求的?LQR函数?
 楼主| 发表于 2005-10-24 11:17 | 显示全部楼层
K,M分别是结构的刚镀矩阵和质量矩阵,在matlab中用的是lqr函数,谢谢
发表于 2005-10-24 12:31 | 显示全部楼层
你最好用钟万勰院士的精细积分做,可以得到计算机上的解析解
 楼主| 发表于 2005-10-25 21:26 | 显示全部楼层
用精细解发也不是办法啊,那个前提要求也是这个A矩阵特征值实部小于0,但是目前得到的这个矩阵的特征值实部很小,但是不全部小于0啊,也就是A并不稳定,该怎么作那?<BR>求求各位,看有没有什么其他的办法了
发表于 2005-11-6 10:17 | 显示全部楼层
按照你的说法,你这个系统应该是不稳定,能否把模型给得具体点
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-24 20:14 , Processed in 0.069888 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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