源程序.txt

本软件研究了地基的计算方法

%输入纸3，//a0.3,W0.4的系数%
  A(4,6)=1/2*k2*(D(1,5)*(1-u(1,5)*u(1,5))+2*EI1/h1);
%a"(0,4)=1/2*k2*(D(0,4)*(1-u(0,4)*u(0,4))+2*EI1/h1) %
%输入纸4，//a0.4,W0.5的系数%
for j=5:N2-3
    A(j,j+2)=1/2*k2*(D(1,j+1)*(1-u(1,j+1)*u(1,j+1))+2*EI1/h1);
    %a"(0,j)=1/2*k2*(D(0,j)*(1-u(0,j)*u(0,j))+2*EI1/h1)%
    %输入纸5,//a0.k，W0.(j+2)的系数,其中K=j+1,即k=5,6,...,N2-3%
end
 A(N2-2,N2)=1/2*k2*(D(1,N2-1)*(1-u(1,N2-1)*u(1,N2-1))+2*EI1/h1)+5/2*EI1*k2/h1;
%a"(0,N2-2)=1/2*k2*(D(0,N2-2)*(1-u(0,N2-2)*u(0,N2-2))+2*EI1/h1)+5/2*EI1*k2/h1%
%输入纸6，//a0.(N2-2),W0.(N2-2)的系数%
A(N2-1,column)=1/2*k2*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI1/h1)+4*EI1*k2/h1;
% a"(0,N2-1)=1/2*k2*(D(0,N2-1)*(1-u(0,N2-1)*u(0,N2-1))+2*EI1/h1)+4*EI1*k2/h1%
%输入纸6，//a0.(N2-1),W0.(N2-1)的系数%

 %             下面为地基中的第二列到第N1列               %
for i=N2+2:column:((N1-1)*column+1)          %从N2+2处开始,步长为N2+1%
    for j=i:i+N2-2
        A(j,j+2)=k2*D(i-bc2*N2,j-i+2);
        %输入论文54,//aij，其中i=1,...,N1-1,j=1,...,N2-1%//aij=????????????
    end
    bc2=bc2+1;
end 

 %             下面为地基中的第N1+1列               %
A(N1*column+1,N1*column+3)=1/2*k2*(D(row,2)*(1-u(row,2)*u(row,2))+2*EI3/h1)+4*EI3*k2/h1;
%                 a"(N1,1)=1/2*k2*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI3/h1)+4*EI3*k2/h1%
%输入纸24，//aN1.1,WN1.2的系数%
A(N1*column+2,N1*column+4)=1/2*k2*(D(row,3)*(1-u(row,3)*u(row,3))+2*EI3/h1)+5/2*EI3*k2/h1;
%                 a"(N1,2)=1/2*k2*(D(N1,2)*(1-u(N1,2)*u(N1,2))+2*EI3/h1)+5/2*EI3*k2/h1 %
%输入纸25，//aN1.2,WN1.3的系数%
A(N1*column+3,N1*column+5)=1/2*k2*(D(row,4)*(1-u(row,4)*u(row,4))+2*EI3/h1);
%                 a"(N1,3)=1/2*k2*(D(N1,3)*(1-u(N1,3)*u(N1,3))+2*EI3/h1)%
%输入纸26，//aN1.3,WN1.4的系数%
A(N1*column+4,N1*column+6)=1/2*k2*(D(row,5)*(1-u(row,5)*u(row,5))+2*EI3/h1);
%                 a"(N1,4)=1/2*k2*(D(N1,4)*(1-u(N1,4)*u(N1,4))+2*EI3/h1) %
%输入纸27，//aN1.4,WN1.5的系数%
for j=N1*column+5:row*column-4
    A(j,j+2)=1/2*k2*(D(row,j-N1*column+1)*(1-u(row,j-N1*column+1)*u(row,j-N1*column+1))+2*EI3/h1);
   %a"(N1,j)=1/2*k2*(D(N1,j)*(1-u(N1,j)*u(N1,j))+2*EI3/h1)%
    %输入纸28,//aN1.k，WN1.(j+2)的系数,其中K=5,6,...,N2-3%
end
A(row*column-3,row*column-1)=1/2*k2*(D(row,N2-1)*(1-u(row,N2-1)*u(row,N2-1))+2*EI3/h1)+5/2*EI3*k2/h1;
%                      a"(N1,N2-2)=1/2*k2*(D(N1,N2-2)*(1-u(N1,N2-2)*u(N1,N2-2))+2*EI3/h1)+5/2*EI3*k2/h1 %
%输入纸29，//aN1.(N2-2),WN1.(N2-1)的系数%
A(row*column-2,row*column)=1/2*k2*(D(row,N2)*(1-u(row,N2)*u(row,N2))+2*EI3/h1)+4*EI3*k2/h1;
%                    a"(N1,N2-1)=1/2*k2*(D(N1,N2-1)*(1-u(N1,N2-1)*u(N1,N2-1))+2*EI3/h1)+4*EI3*k2/h1%
%输入纸29，//aN1.(N2-1),WN1.(N2-1)的系数%
%................................................................................................................................
%             下面为加梁EI1,EI2,EI3,EI4后增加的点               %
A(1,4)=(-EI1)*k2/h1;             %输入纸1，///a0.0,W0.3的系数%                                             
A(N2-2,column)=(-EI1)*k2/h1;       %输入纸1，///a0.0%

A(N1*column+1,N1*column+4)=(-EI3)*k2/h1; %输入纸24，///b0.0%
A(row*column-3,row*column)=(-EI3)*k2/h1; %输入纸24，///b0.0%

A(1,3*N2+4)=(-EI2)*k1/h2;       %输入纸1，k0.0,W3.0的系数%
A(column,4*N2+4)=(-EI4)*k1/h2;       %输入纸9，k0.1%

A((N1-3)*N2+N1-2,N1*column+1)=(-EI2)*k1/h2;       %输入纸1，k0.0%
A((N1-2)*N2+N1-2,row*column)=(-EI4)*k1/h2;       %输入纸9，k0.1%
%...........................................................................
%                         下面为S三对角阵跟T对角阵            %
 %                 输入S中的数据元素：对角线元素                  %
  %             下面为地基中的第2列               %
A(1,N2+2)=-2*DL(1,1)*(1-ul(1,1))-k1*(D(2,1)*(1-u(2,1)*u(2,1))+2*EI2/h2)-5*EI2*k1/h2;
%  S(1,0)=-2*DL(1,1)*(1-ul(1,1))-k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)-5*EI2*k1/h2%
%输入纸1，S1.0,W1.0的系数%
for j=2:N2
    A(j,j+column)=-2*D(2,j)*(k1+u(2,j))-2*DL(1,j-1)*(1-ul(1,j-1))-2*DL(1,j)*(1-ul(1,j));
    %    S(1,j)=-2*D(1,j)*(k1+u(1,j))-2*DL(1,j)*(1-ul(1,j))-2*DL(1,j+1)*(1-ul(1,j+1))%
    %输入论文55,S1.j，其中j=1,2,...,N2-1%
end
A(column,2*N2+2)=-2*DL(1,N2)*(1-ul(1,N2))-k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)-5*EI4*k1/h2;
%      S(1,N2)=-2*DL(1,N2)*(1-ul(1,N2))-k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)-5*EI4*k1/h2%
%输入纸9，S1.N2,W1.N2的系数%

%             下面为地基中的第3列               %
A(N2+2,2*N2+3)=-2*DL(2,1)*(1-ul(2,1))-k1*(D(2,1)*(1-u(2,1)*u(2,1))+D(3,1)*(1-u(3,1)*u(3,1))+4*EI2/h2)-10*EI2*k1/h2;
%       S(2,0)=-2*DL(2,1)*(1-ul(2,1))-k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)-k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)-10*EI2*k1/h2%
%输入纸10，S2.0,W2.0的系数%
for j=N2+3:2*column-1
    A(j,j+column)=-2*D(2,j-N2-1)*(k1+u(2,j-N2-1))-2*D(3,j-N2-1)*(k1+u(3,j-N2-1))-2*DL(2,j-N2-2)*(1-ul(2,j-N2-2))-2*DL(2,j-N2-1)*(1-ul(2,j-N2-1));
    %     S(2,j)=-2*D(1,j)*(k1+u(1,j))          -2*D(2,j)*(k1+u(2,j))          -2*DL(2,j)*(1-ul(2,j))          -2*DL(2,j+1)*(1-ul(2,j+1))%
    %输入论文55,Si.j，其中i=2,j=1,2,...,N2-1%
end
A(2*N2+2,3*N2+3)=-2*DL(2,N2)*(1-ul(2,N2))-k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)-k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)-10*EI4*k1/h2;
%        S(2,N2)=-2*DL(2,N2)*(1-ul(2,N2))-k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)      -k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)-10*EI4*k1/h2 %
%输入纸11，S2.N2,W2.N2的系数%

%             下面为地基中的第4列               %
A(2*N2+3,3*N2+4)=-2*DL(3,1)*(1-ul(3,1))-k1*(D(3,1)*(1-u(3,1)*u(3,1))+D(4,1)*(1-u(4,1)*u(4,1))+4*EI2/h2)-2*EI2*k1/h2;
%         S(3,0)=-2*DL(3,1)*(1-ul(3,1))-k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)-k1*(D(3,0)*(1-u(3,0)*u(3,0))+2*EI2/h2)-2*EI2*k1/h2%
%输入纸12，S3.0,W3.0的系数%
for j=2*N2+4:3*N2+2
    A(j,j+column)=-2*D(3,j-2*N2-2)*(k1+u(3,j-2*N2-2))-2*D(4,j-2*N2-2)*(k1+u(4,j-2*N2-2))-2*DL(3,j-2*N2-3)*(1-ul(3,j-2*N2-3))-2*DL(3,j-2*N2-2)*(1-ul(3,j-2*N2-2));
   %     S(3,j)=-2*D(2,j)*(k1+u(2,j))              -2*D(3,j)*(k1+u(3,j))              -2*DL(3,j)*(1-ul(3,j))              -2*DL(3,j+1)*(1-ul(3,j+1)) %
    %输入论文55,Si.j，其中i=3,j=1,2,...,N2-1%
end
A(3*N2+3,4*N2+4)=-2*DL(3,N2)*(1-ul(3,N2))-k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)-k1*(D(4,column)*(1-u(4,column)*u(4,column))+2*EI4/h2)-2*EI4*k1/h2  ;
%        S(3,N2)=-2*DL(3,N2)*(1-ul(3,N2))-k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)      -k1*(D(3,N2)*(1-u(3,N2)*u(3,N2))+2*EI4/h2)-2*EI4*k1/h2 %
%输入纸13，S3.N2,W3.N2的系数%

%             下面为地基中的第5列               %
A(3*N2+4,4*N2+5)=-2*DL(4,1)*(1-ul(4,1))-k1*(D(4,1)*(1-u(4,1)*u(4,1))+2*EI2/h2)-k1*(D(5,1)*(1-u(5,1)*u(5,1))+2*EI2/h2);
%         S(4,0)=-2*DL(4,1)*(1-ul(4,1))-k1*(D(3,0)*(1-u(3,0)*u(3,0))+2*EI2/h2)-k1*(D(4,0)*(1-u(4,0)*u(4,0))+2*EI2/h2) %
%输入纸14，S4.0,W4.0的系数%
for j=3*N2+5:4*N2+3
    A(j,j+column)=-2*D(4,j-3*N2-3)*(k1+u(4,j-3*N2-3))-2*D(5,j-3*N2-3)*(k1+u(5,j-3*N2-3))-2*DL(4,j-3*N2-4)*(1-ul(4,j-3*N2-4))-2*DL(4,j-3*N2-3)*(1-ul(4,j-3*N2-3));
    %    S(4,j)=-2*D(3,j)*(k1+u(3,j))              -2*D(4,j)*(k1+u(4,j))              -2*DL(4,j)*(1-ul(4,j))              -2*DL(4,j+1)*(1-ul(4,j+1)) %
    %输入论文55,Si.j，其中i=4,j=1,2,...,N2-1%
end
A(4*N2+4,5*N2+5)=-2*DL(4,N2)*(1-ul(4,N2))-k1*(D(4,column)*(1-u(4,column)*u(4,column))+2*EI4/h2)-k1*(D(5,column)*(1-u(5,column)*u(5,column))+2*EI4/h2);
%        S(4,N2)=-2*DL(4,N2)*(1-ul(4,N2))-k1*(D(3,N2)*(1-u(3,N2)*u(3,N2))+2*EI4/h2)      -k1*(D(4,N2)*(1-u(4,N2)*u(4,N2))+2*EI4/h2)%
%输入纸15，S4.N2,W4.N2的系数%

%             下面为地基中的第6列到第N1-2列   i=5,6...N1-3            %
for i=4*N2+5:column:N2*(N1-4)+N1-3
    A(i,i+column)=-2*DL(i-bc3*N2,1)*(1-ul(i-bc3*N2,1))-k1*(D(i-bc3*N2,1)*(1-u(i-bc3*N2,1)*u(i-bc3*N2,1))+D(i-bc3*N2+1,1)*(1-u(i-bc3*N2+1,1)*u(i-bc3*N2+1,1))+4*EI2/h2);
    %    S(i,0)=-2*DL(i,1)*(1-ul(i,1))              -k1*(D(i-1,0)*(1-u(i-1,0)*u(i-1,0))+2*EI2/h2)               -k1*(D(i,0)*(1-u(i,0)*u(i,0))+2*EI2/h2)%
    %输入纸16，Sk.0,W(i+1).0的系数,其中k=i+1，即k=5，...，N1-3%
    for j=i+1:i+N2-1
        A(j,j+column)=-2*D(i-bc3*N2,j-i+1)*(k1+u(i-bc3*N2,j-i+1))-2*D(i-bc3*N2+1,j-i+1)*(k1+u(i-bc3*N2+1,j-i+1))-2*DL(i-bc3*N2,j-i)*(1-ul(i-bc3*N2,j-i))-2*DL(i-bc3*N2,j-i+1)*(1-ul(i-bc3*N2,j-i+1));
        %    S(i,j)=-2*D(i-1,j)*(k1+u(i-1,j))                  -2*D(i,j)*(k1+u(i,j))                          -2*DL(i,j)*(1-ul(i,j))                  -2*DL(i,j+1)*(1-ul(i,j+1))%
        %输入论文55,Si.j，其中i=k,j=1,2,...,N2-1%
    end
    A(i+N2,i+2*N2+1)=-2*DL(i-bc3*N2,N2)*(1-ul(i-bc3*N2,N2))-k1*(D(i-bc3*N2,column)*(1-u(i-bc3*N2,column)*u(i-bc3*N2,column))+2*EI4/h2)-k1*(D(i-bc3*N2+1,column)*(1-u(i-bc3*N2+1,column)*u(i-bc3*N2+1,column))+2*EI4/h2);
    %        S(i,N2)=-2*DL(i,N2)*(1-ul(i,N2))              -k1*(D(i-1,N2)*(1-u(i-1,N2)*u(i-1,N2))+2*EI4/h2)                     -k1*(D(i,N2)*(1-u(i,N2)*u(i,N2))+2*EI4/h2)%
    %输入纸17，Sk.N2,W(i+1).N2的系数%
    bc3=bc3+1;
end

%             下面为地基中的第N1-1列               %
A((N1-3)*column+1,(N1-2)*column+1)=-2*DL(N1-2,1)*(1-ul(N1-2,1))-k1*(D(N1-2,1)*(1-u(N1-2,1)*u(N1-2,1))+D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+4*EI2/h2)-2*EI2*k1/h2;
%                        S(N1-2,0)=-2*DL(N1-2,1)*(1-ul(N1-2,1))-k1*(D(N1-3,0)*(1-u(N1-3,0)*u(N1-3,0))+2*EI2/h2)-k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)-2*EI2*k1/h2%
%输入纸18，S(N1-2).0,W(N1-2).0的系数%
for j=(N1-3)*column+2:(N1-2)*column-1
    A(j,j+column)=-2*D(N1-2,j-(N1-3)*column)*(k1+u(N1-2,j-(N1-3)*column))-2*D(N1-1,j-(N1-3)*column)*(k1+u(N1-1,j-(N1-3)*column))-2*DL(N1-2,j-(N1-3)*column-1)*(1-ul(N1-2,j-(N1-3)*column-1))-2*DL(N1-2,j-(N1-3)*column)*(1-ul(N1-2,j-(N1-3)*column));
    % S(N1-2,j)=-2*D(N1-3,j)*(k1+u(N1-3,j))                            -2*D(N1-2,j)*(k1+u(N1-2,j))                            -2*DL(N1-2,j)*(1-ul(N1-2,j))                                -2*DL(N1-2,j+1)*(1-ul(N1-2,j+1))%
    %输入论文55,Si.j，其中i=N1-2,j=1,2,...,N2-1%
end
A((N1-2)*column,(N1-1)*column)=-2*DL(N1-2,N2)*(1-ul(N1-2,N2))-k1*(D(N1-2,column)*(1-u(N1-2,column)*u(N1-2,column))+2*EI4/h2)-k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)-2*EI4*k1/h2;
%                   S(N1-2,N2)=-2*DL(N1-2,N2)*(1-ul(N1-2,N2))-k1*(D(N1-3,N2)*(1-u(N1-3,N2)*u(N1-3,N2))+2*EI4/h2)      -k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)-2*EI4*k1/h2%
%输入纸19,W(N1-2).N2的系数%

%             下面为地基中的第N1列               %
A((N1-2)*column+1,(N1-1)*column+1)=-2*DL(N1-1,1)*(1-ul(N1-1,1))-k1*(D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+2*EI2/h2)-k1*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI2/h2)-10*EI2*k1/h2;
%                        S(N1-1,0)=-2*DL(N1-1,1)*(1-ul(N1-1,1))-k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)-k1*(D(N1-1,0)*(1-u(N1-1,0)*u(N1-1,0))+2*EI2/h2)-10*EI2*k1/h2%
%输入纸20,S(N1-1).0,W(N1-1).0的系数%
for j=(N1-2)*column+2:(N1-1)*column-1
    A(j,j+column)=-2*D(N1-1,j-(N1-2)*column)*(k1+u(N1-1,j-(N1-2)*column))-2*D(N1,j-(N1-2)*column)*(k1+u(N1,j-(N1-2)*column))-2*DL(N1-1,j-(N1-2)*column-1)*(1-ul(N1-1,j-(N1-2)*column-1))-2*DL(N1-1,j-(N1-2)*column)*(1-ul(N1-1,j-(N1-2)*column));
    % S(N1-1,j)=-2*D(N1-2,j)*(k1+u(N1-2,j))                            -2*D(N1-1,j)*(k1+u(N1-1,j))                        -2*DL(N1-1,j)*(1-ul(N1-1,j))                                -2*DL(N1-1,j+1)*(1-ul(N1-1,j+1))%
    %输入论文55,Si.j，其中i=N1-1,j=1,2,...,N2-1%
end
A((N1-1)*column,N1*column)=-2*DL(N1-1,N2)*(1-ul(N1-1,N2))-k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)-k1*(D(N1,column)*(1-u(N1,column)*u(N1,column))+2*EI4/h2)-10*EI4*k1/h2 ;