我写信是弯箱梁桥的有限元代码和错误. .警告:矩阵接近奇异或严重了。结果可能是不准确的。RCOND = 9.359430 e-32。在Dec14 >(第579行)

4视图(30天)

答案(3)


沃尔特·罗伯森
沃尔特·罗伯森 2018年12月16日
与rcond小你可以确保你的矩阵是奇异的,不存在逆。如果你检查排名()它可能会告诉你没有满秩。
你需要再核对你的网状结构。可能的情况是,你的解决方案技术并不适用于您的情况。
你可以沉默的错误使用pinv相反如果发票,但是答案是不可能有意义的情况。

约翰D 'Errico
约翰D 'Errico 2018年12月16日
有限元的方法是识别原因有几个潜在的奇异矩阵。最明显的一个是你不可能建造它。
也许是可能的一个原因是你没有适当的约束问题。这是一个容易犯的错误。例如,考虑一个非常简单的桁架,但是你并没有固定的至少一个节点在空间桁架specfific位置。结果将是您可以翻译整个桁架在任何地方,通过添加一个简单的抵消整个桁架,不引入任何能量损失的形式桁架的变形。在这种情况下,将奇异方程组,奇异矩阵。这是一个问题吗?从某种意义上说,是的,因为系统的解不是唯一的。
实际上,pinv能解决问题。至少,它将有效地解决问题。更好的不过是解决问题的正确的约束集防止桁架处罚自由翻译,或者从能源自由旋转。

类别

找到更多的在结构分析帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!