pasuka 发表于 2008-1-28 21:29:39

子空间迭代法和刚度矩阵的问题

建立的Mindlin板单元模型并推导刚度矩阵、质量矩阵,采用子空间迭代法计算自由-自由边界条件下前几阶固有频率的时候,出现错误,程序显示说:STIFFNESS MATRIX NOT POSITIVE DEFINITE 刚度矩阵非正定
参考《有限单元法》第三版p495页的说法,需要按照给定边界条件修改刚度矩阵K
请教大侠,应该怎么样进行修改呢?
谢谢

molen 发表于 2008-1-29 00:39:59

你的《有限单元法》书上肯定有施加边界条件的方法。看来你对有限元过程还未完全掌握,这样去写程序是比较危险的。
简单来说,是把约束的自由度相应的刚度阵中的行和列删除。

pasuka 发表于 2008-1-29 23:33:02

但是我老师告诉我,不需要矩阵降阶也是可以计算,只是会碰到特征值为0的刚体位移模态

dsxu1981_01 发表于 2008-1-30 18:03:38

pasuka老兄说的没错,某个方向不加约束会出现该方向的刚体位移模态
你可以ANSYS就能试出来。
一般我们编程都是处理掉约束的。
可以用乘大数或者划行划列法,都可以的。
这些方法一般有限元书上都会有的。
页: [1]
查看完整版本: 子空间迭代法和刚度矩阵的问题