如何求解刚度矩阵?

计算科学 有限元 matlab
2021-12-13 14:35:37

我在这个领域很新。对于我的大学,我准备了一个刚度矩阵来解决一个项目组。这个矩阵由 450 个方程和 450 个未知数组成(这是一个 Matlab 脚本),我有相关的负载向量{L}450 行。因为我需要计算位移,你会建议我简单地做吗?{x}=[K]1{L}

1个回答

正如 Wolfgang Bangerth 在评论中提到的,对于这个问题,你想使用反斜杠 操作符,即

u = K \ L;

一般来说,你不想重新发明轮子,应该使用线性代数库来解决这类问题。

PS:请记住,如果您没有施加任何约束,您的刚度矩阵可能是奇异的。