【求助】MATLAB结构有限元分析，总提示矩阵奇异

glenn

各位MATLAB达人，大家好，请教一个问题！欢迎指导！！{:3_54:}
参考徐斌老师的《MATLAB有限元结构动力学分析与工程应用》一书。
编写了一小段计算简支梁的变形和频率的程序，发现几个问题，归纳如下：
1.采用书中提供的FrameElement21和FrameElement22两种单元（分别对应欧拉梁和铁木辛柯梁）；欧拉梁可以计算出变形，铁木辛柯梁在用kk\mm计算节点位移的时候提示有奇异矩阵。 PS：因此，下面的问题都是用欧拉梁计算时碰到的问题。
2.尝试改过密度，改了之后发现静变形没有一点变化。。我去看了FrameElement21这个计算单元刚度阵和单元质量阵的子程序，发现密度应该会影响质量阵，，然而事实上计算结果（静变形）没有变化。PS：频率及总质量阵有所变化，这里我想到了，不考虑重力确实不会发生变化。
3.计算阶跃载荷动力响应的时候，加阻尼这一步很难找到合适的数字，假设阻尼阵为质量阵和刚度阵的线性组合，质量阵和刚度阵的系数ac，bc应该是有人提出过经验公式，但是我做的时候只有把这两个数加得极小才不会提示过阻尼（小到0.0000002这样），并且算出来的响应比静态变形大很多（按照理论的话应该是和静变形一样大）。

目前发现问题就是这些！
现将程序附上。

pasuka

那本书不值得阅读，拼凑出来的三流评职称书籍
建议网上找找A. J. M. Ferreira的《MATLAB Codes for Finite Element Analysis》，浅显易懂

glenn

pasuka 发表于 2012-3-28 21:48
那本书不值得阅读，拼凑出来的三流评职称书籍
建议网上找找A. J. M. Ferreira的《MATLAB Codes for Finite  ...

好的，谢谢你！{:3_57:}

希望再了解下，比例阻尼一般aM+bK中的a和b取多少呀   {:3_66:}