有限元计算梁的模态
小弟最近写了个有限元的程序，求解梁的模态遇到这样的问题：
有限元离散出来的刚度矩阵肯定是奇异的
而且在利用子空间迭代求解时
要将刚度矩阵进行LTLD分解
这样必定会导致分解后的矩阵对角线元素为零
清华张雄老师《计算动力学》上提到，采用移轴法将刚度矩阵消除奇异性
但是在求解过程中仍然出现了矩阵对角线元素为零的问题
采用子空间迭代法求解固有频率
不知哪位大侠有过这方面的经验
请不吝赐教
小弟联系方式：QQ50944981
谢谢 .
一般结构有限元计算时，总刚度和总质量矩阵求广义特征值时，应该先对边界条件进行处理，这样就不存在半正定问题了.. .. .
可以参考一下《结构计算程序设计》张汝清编著，重庆出版社，.. .. 汗论坛里怎么这种让人很不爽的帖子这么多呢！
有没有互动讨论贴了？
以后改原始菜鸟求助中心吧·！ .
可以结合看看
《有限元法程序设计》，杨菊生编，西安交通大学出版社
《结构有限元分析》吴鸿庆编，中国铁道出版社
《Numerical Methods In Finite Analysis》 Bathe 和 Wilson （其中的子空间程序都很好用，包括分解和JACOBI都有..）
这些书讲的结构有限元程序都是一个套路的，编程思想比较好，可以相互补充学习理解. ..
学习程序也要从程序的思想理解，照搬是不可能的。
[ 本帖最后由 欧阳中华 于 2008-7-23 16:54 编辑 ]
页:
[1]