有限元计算梁的模态

flyboyupc · 发表于 2008-7-22 16:53

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

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

x

小弟最近写了个有限元的程序，求解梁的模态
遇到这样的问题：
有限元离散出来的刚度矩阵肯定是奇异的
而且在利用子空间迭代求解时
要将刚度矩阵进行LTLD分解
这样必定会导致分解后的矩阵对角线元素为零
清华张雄老师《计算动力学》上提到，采用移轴法将刚度矩阵消除奇异性
但是在求解过程中仍然出现了矩阵对角线元素为零的问题
采用子空间迭代法求解固有频率

不知哪位大侠有过这方面的经验
请不吝赐教
小弟联系方式：QQ50944981
谢谢

欧阳中华 · 发表于 2008-7-22 18:11

.
一般结构有限元计算时，总刚度和总质量矩阵求广义特征值时，应该先对边界条件进行处理，这样就不存在半正定问题了.. ..

欧阳中华 · 发表于 2008-7-23 12:15

.
可以参考一下《结构计算程序设计》张汝清编著，重庆出版社，.. ..

lq12131010 · 发表于 2008-7-23 13:42

汗论坛里怎么这种让人很不爽的帖子这么多呢！
有没有互动讨论贴了？
以后改原始菜鸟求助中心吧·！

欧阳中华 · 发表于 2008-7-23 16:52

.
可以结合看看

  《有限元法程序设计》，杨菊生编，西安交通大学出版社

  《结构有限元分析》吴鸿庆编，中国铁道出版社

《Numerical Methods In Finite Analysis》 Bathe 和 Wilson （其中的子空间程序都很好用，包括分解和JACOBI都有..）

   这些书讲的结构有限元程序都是一个套路的，编程思想比较好，可以相互补充学习理解. ..

   学习程序也要从程序的思想理解，照搬是不可能的。

[ 本帖最后由欧阳中华于 2008-7-23 16:54 编辑 ]

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

[结构振动] 有限元计算梁的模态

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