nastran计算过程中矩阵奇异的说明
(2014-08-02 12:22:40)
标签:
nastran矩阵奇异 |
分类: nastran |
前面几天贴了一些找到的资料,联系一下自己这几天的计算,整理一下自己的理解。
一般nastran出现矩阵奇异停止计算的原因是某些单元的柔度矩阵太大(意味着刚度太低),nastran认为产生这样大的柔度是不合理的,所以程序停止计算(并不是不能计算)。(当然你用param,bailout,-1这个参数的意义就是告诉nastran “哥们,不用考虑计算数据的合理性,算完它就可以!”)
之所以nastran要终止计算可能的原因是,程序认为边界条件有问题,缺少了约束(相当与Abaqus缺少约束产生无限位移的错误)。但实际情况是有些时候约束没有问题还是会产生这样的错误,主要的原因有2个。
1.是单元没有合并节点,导致有游离单元存在,这些单元没有被约束产生较大的柔度。这时候最好不用使用autospc还是仔细检查模型吧(仅指本人目前的工作的应用,当然如果是其他方面我的意见不一定正确),我个人感觉会对结果有影响。
2.是单元质量有问题,网上的说法是单元不协调。原因是在nastran中shell单元只有5个自由度,按照原理来说shell单元的第6个自由度应该是被和它相连的单元约束住的(也就是shell单元本身是不考虑法向转动的)。但是在计算整体刚度的时候,某些节点的第6个自由度没有了足够的刚度。(shell单元上没有定义,周围的单元又提供不了足够的约束,这种情况在理论上是不应该存在的,但有可能因为网格的原因产生这种问题)这时候可以考虑使用PARAM,K6ROT,100. 0参数,该参数的意义是增加shell单元的第6自由度,也就是在所有的shell单元上提供6个自由度。该参数貌似会增加计算的成本(具体增加多少计算时间不清楚),但可以解决不协调问题,因为所有的shell单元都有6个完整的自由度。