连续梁的matlab编程,矩阵法解内力matlab实现(3跨连续梁)

E=4176000; %E modous unit:

k/in^4

%% TO calcuate K1

L=25; I=20;

K=[12/L^2 6/L -12/L^2 6/L;6/L 4

-6/L 2;-12/L^2 -6/L 12/L^2 -6/L;6/L 2 -6/L 4];

K1=K*E*I/L;

%% TO calcuate K2

L=30; I=30;

K=[12/L^2 6/L -12/L^2 6/L;6/L 4

-6/L 2;-12/L^2 -6/L 12/L^2 -6/L;6/L 2 -6/L 4];

K2=K*E*I/L;

%% TO calcuate K3

L=30; I=30;

K=[12/L^2 6/L -12/L^2 6/L;6/L 4

-6/L 2;-12/L^2 -6/L 12/L^2 -6/L;6/L 2 -6/L 4];

K3=K*E*I/L;

%% for calculate

Kglobal

Kglobal1=zeros(8,8);

Kglobal2=zeros(8,8);

Kglobal3=zeros(8,8);

Kglobal1(1:4,1:4)=K1;

Kglobal2(3:6,3:6)=K2;

Kglobal3(5:8,5:8)=K3;

Kglobal=Kglobal1+Kglobal2+Kglobal3;

%% for f, Kdis

f4=-45.833; f6=75;