矩阵求逆时出现INF,怎么处理一下很简单的一段代码
clc;
E=0.217e12;
b=0.028;
t=0.003;
L=0.328;
A=b*t;
I=b*t^3/12;
n=9; %单元个数
%单元刚度矩阵
k1=E*A/L;
Ke=[k1 0 0 -k1 0 0;
0 12*E*I/L^3 6*E*I/L^2 0 -12*E*I/L^3 -6*E*I/L^2;
0 6*E*I/L^2 4*E*I/L 0 6*E*I/L^2 2*E*I/L;
-k1 0 0 k1 0 0;
0 -12*E*I/L^3 -6*E*I/L^2 0 12*E*I/L^3 -6*E*I/L^2;
0 6*E*I/L^2 2*E*I/L 0 -6*E*I/L^2 4*E*I/L];
%组集
K=zeros((n+1)*3);
for i=0:3:length(K)-length(Ke)
K(i+1:i+6,i+1:i+6)=K(i+1:i+6,i+1:i+6)+Ke;