我写信是弯箱梁桥的有限元代码和错误. .警告:矩阵接近奇异或严重了。结果可能是不准确的。RCOND = 9.359430 e-32。在Dec14 >(第579行)
4视图(30天)
显示旧的评论
X_matrix =发票(K_modi) * F_matrix;
答案(3)
约翰D 'Errico
2018年12月16日
有限元的方法是识别原因有几个潜在的奇异矩阵。最明显的一个是你不可能建造它。
也许是可能的一个原因是你没有适当的约束问题。这是一个容易犯的错误。例如,考虑一个非常简单的桁架,但是你并没有固定的至少一个节点在空间桁架specfific位置。结果将是您可以翻译整个桁架在任何地方,通过添加一个简单的抵消整个桁架,不引入任何能量损失的形式桁架的变形。在这种情况下,将奇异方程组,奇异矩阵。这是一个问题吗?从某种意义上说,是的,因为系统的解不是唯一的。
实际上,pinv能解决问题。至少,它将有效地解决问题。更好的不过是解决问题的正确的约束集防止桁架处罚自由翻译,或者从能源自由旋转。