有限元八节点六面体单元，用MATLAB编出总刚的程序，代入动力学方程求解，显示奇异，如何求解呢？

总刚度矩阵本来就具有对称、奇异、阶数等于总自由度数等特点。
注意平衡方程[K]{u}={F}只是用到了力的平衡，而有限元法又是”位移解法“，当不施加位移约束条件时，必然不能求得唯一解，所以刚度矩阵奇异。
求解的常规做法有：(1)将已知位移约束的方程划去，从而求得未知位移，就可以求出未知节点力。
(2)冲大数法，求出节点位移，之后容易得到节点力