MATLAB矩阵位移法

用MATLAB写的矩阵位移法代码，可以解决很多结构力学的问题，结构工程电算化的起点。

for i=1n
    wk=qk（EAIl（i））
    wl=ql（a（i））;
    kk2=wkwl;
    kk1（6i-56i）=kk2;
    kk=wl‘kk2
    k（6i-56i）=kk;
end
qqq=max（max（r））;
K=zeros（mm）;
KR=zeros（qqq-mm）;
for p=1n
    aa=r（p）;
    for i=16
        for j=16
            if aa（i）=m & aa（j）=m
                kkk=k（6p-56p）;
                K（aa（i）aa（j））=K（aa（i）aa（j））+kkk（ij）;
            end
            if aa（i）=m+1 & aa（j）=m
                kkk=k（6p-56p）;
                KR（aa（i）-maa（j））=KR（aa（i）-maa（j））+kkk（ij）;
            end
        end
    end
end
K
KR
Fell=zeros（1qqq）;
for p=1n                                                
    bb=r（p）;
    tt=ql（a（p））;
    ll=tt‘Fg（p）;
    lll=ll‘;
    for i=16                                                       
        Fell（1bb（i））=Fell（1bb（i））+lll（1i）;                                                                                                                                                                                                                                                                               
    end
end
Fel=-Fell‘;
Fej=Fed+Fel;
Fejj=Fej‘;
Fe111=Fejj（1m）;
Fe1=Fe111‘
Fe222=Fejj（m+1qqq）;
Fe2=Fe222‘
d=KFe1
d1=d‘;
d1（m+1qqq）=0;
for p=1n                                                
    cc=r（p）;
    for i=16                                                       
        dd（1i）=d1（1cc（i））;                                                                                                                                                                                                                                                                            
    end
    p
    F=kk1（6p-56p）dd‘+Fg（p）
end
FR=KRd-Fe2
定义函数qk
function y=qk（EAIl）
y=[EAl 0 0 -EAl 0 0;0 12EIlll -6EIll 0 -12EIlll -6EIll;0 -6EIll 4EIl 0 6EIll 2EIl;-EAl 0 0 EAl 0 0;0 -12EIlll 6EIll 0 12EIlll 6EIll;0 -6EIll 2EIl 0 6EIll 4EIl];
定义函数ql
function yy=ql（a）
yy=[cos（a） sin（a） 0 0 0 0;-sin（a） cos（a） 0 0 0 0;0 0 1 0 0 0;0 0 0 cos（a） sin（a） 0;0 0 0 -sin（a） cos（a） 0;0 0 0 0 0 1];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2093  2009-05-12 11:58  矩阵位移法.m

----------- ---------  ---------- -----  ----

                 2093                    1