源程序.txt

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

%               S(N1-1,N2)=-2*DL(N1-1,N2)*(1-ul(N1-1,N2))-k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)      -k1*(D(N1-1,N2)*(1-u(N1-1,N2)*u(N1-1,N2))+2*EI4/h2)-10*EI4*k1/h2%
%输入纸21,S(N1-1).N2,W(N1-1).N2的系数%

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

%             输入S中的数据元素：对角线上一行的元素                  %
  %             下面为地基中的第2列               %
A(1,N2+3)=2*DL(1,1)*(1-ul(1,1));
% S'(1,1)=2*DL(1,1)*(1-ul(1,1))%
%输入论文55,/S1.1%
for j=2:N2
    A(j,j+N2+2)=D(2,j)*u(2,j)+2*DL(1,j)*(1-ul(1,j));
    %   S'(1,j)=D(1,j-1)*u(1,j-1)+2*DL(1,j)*(1-ul(1,j))%
    %输入论文55,/S1.j,其中j=2，3，...，N2%
end
 
%             下面为地基中的第3列到第N1列   i=2,3,...,N1-1            %
for i=N2+2:column:column*(N1-2)+1
    A(i,i+N2+2)=D(i-bc4*N2,2)*u(i-bc4*N2,2)+2*DL(i-bc4*N2,1)*(1-ul(i-bc4*N2,1));
    %   S'(i,1)=D(i-1,1)*u(i-1,1)          +2*DL(i,1)*(1-ul(i,1))%
    %输入论文55,/Si.1%
    for j=i+1:i+N2-2
        A(j,j+N2+2)=D(i-bc4*N2,j-i+2)*u(i-bc4*N2,j-i+2)+D(i-bc4*N2+1,j-i+1)*u(i-bc4*N2+1,j-i+1)+2*DL(i-bc4*N2,j-i+1)*(1-ul(i-bc4*N2,j-i+1));
        %   S'(i,j)=D(i-1,j)*u(i-1,j)                  +D(i,j-1)*u(i,j-1)                      +2*DL(i,j)*(1-ul(i,j)) %
        %输入论文55,/Si.j，其中j=2，3，...，N2-1，i=2，3，...，N1-1%
    end
    A(i+N2-1,i+2*N2+1)=D(i-bc4*N2+1,N2)*u(i-bc4*N2+1,N2)+2*DL(i-bc4*N2,N2)*(1-ul(i-bc4*N2,N2));
    %         S'(i,N2)=D(i,N2-1)*u(i,N2-1)              +2*DL(i,N2)*(1-ul(i,N2))%
    %输入论文55,/[e-mail]
    bc4=bc4+1;
end
 
%             下面为地基中的第N1+1列               %
for j=column*(N1-1)+1:column*N1-2
    A(j,j+N2+2)=D(N1,j-column*(N1-1)+1)*u(N1,j-column*(N1-1)+1)+2*DL(N1,j-column*(N1-1))*(1-ul(N1,j-column*(N1-1))) ;
    %  S'(N1,j)=D(N1-1,j)*u(N1-1,j)                            +2*DL(N1,j)*(1-ul(N1,j))%
    %输入论文55,/SN1.j,其中j=1，2，...，N2-1%
end
A(column*N1-1,column*row)=2*DL(N1,N2)*(1-ul(N1,N2));
%                  S'(N1,N2)=2*DL(N1,N2)*(1-ul(N1,N2)) %
%输入论文55,/SN1.N2%

%             输入S中的数据元素：对角线下一行的元素                  %
  %             下面为地基中的第2列               %
for j=2:N2
    A(j,j+N2)=D(2,j)*u(2,j)+2*DL(1,j-1)*(1-ul(1,j-1)) ;
    %S''(1,j)=D(1,j)*u(1,j)+2*DL(1,j)*(1-ul(1,j))%
    %输入论文56,//S1.j,其中j=1，2，...，N2-1%
end
A(column,2*N2+1)=2*DL(1,N2)*(1-ul(1,N2));
%    S''(1,N2)=2*DL(1,N2)*(1-ul(1,N2))%
%输入论文56,//S1.N2%
 
%             下面为地基中的第3列到第N1列 i=2,3,...,N1-1              %
for i=N2+3:column:column*(N1-2)+2
    A(i,i+N2)=D(i-bc5*N2,2)*u(i-bc5*N2,2)+2*DL(i-bc5*N2-1,1)*(1-ul(i-bc5*N2-1,1));
    %S''(i,1)=D(i,1)*u(i,1)              +2*DL(i,1)*(1-ul(i,1))%
    %输入论文56,//Si.1%
    for j=i+1:i+N2-2
        A(j,j+N2)=D(i-bc5*N2-1,j-i+1)*u(i-bc5*N2-1,j-i+1)+D(i-bc5*N2,j-i+2)*u(i-bc5*N2,j-i+2)+2*DL(i-bc5*N2-1,j-i+1)*(1-ul(i-bc5*N2-1,j-i+1));
        %S''(i,j)=D(i-1,j-1)*u(i-1,j-1)                  +D(i,j)*u(i,j)                      +2*DL(i,j)*(1-ul(i,j))%
        %输入论文56,//Si.j，其中j=2，3，...，N2-1，i=2，3，...，N1-1%
    end
    A(i+N2-1,i+2*N2-1)=D(i-bc5*N2-1,N2)*u(i-bc5*N2-1,N2)+2*DL(i-bc5*N2-1,N2)*(1-ul(i-bc5*N2-1,N2));
    %        S''(i,N2)=D(i-1,N2-1)*u(i-1,N2-1)          +2*DL(i,N2)*(1-ul(i,N2))%
    %输入论文56,//Si.N2%
    bc5=bc5+1;
end
 
%             下面为地基中的第N1+1列               %
A(column*(N1-1)+2,column*N1+1)=2*DL(N1,1)*(1-ul(N1,1));
%                    S''(N1,1)=2*DL(N1,1)*(1-ul(N1,1))%
%输入论文56,//SN1.1%
for j=column*(N1-1)+3:column*N1
    A(j,j+N2)=D(N1,j-column*(N1-1)-1)*u(N1,j-column*(N1-1)-1)+2*DL(N1,j-column*(N1-1)-1)*(1-ul(N1,j-column*(N1-1)-1));
   %S''(N1,j)=D(N1-1,j-1)*u(N1-1,j-1)                        +2*DL(N1,j)*(1-ul(N1,j))%
    %输入论文56,//SN1.j,其中j=2，3，...，N2%
end


 %                 输入T中的数据元素：对角线元素                  %
 %             下面为地基中的第3列               %
A(1,2*N2+3)=(1/2)*k1*(D(2,1)*(1-u(2,1)*u(2,1))+2*EI2/h2)+4*EI2*k1/h2;
%    t(1,0)=(1/2)*k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)+4*EI2*k1/h2%
%输入纸1,t1.0,W2.0的系数%
for j=2:N2
    A(j,j+2*N2+2)=k1*D(2,j) ;
    %      t(1,j)=k1*D(1,j)%
    %输入论文54,ti.j,其中i=1,j=1,2,...,N2-1%
end 
A(column,3*N2+3)=(1/2)*k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)+4*EI4*k1/h2;
%      t(1,N2)=(1/2)*k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)+4*EI4*k1/h2%
%输入纸9,t1.N2,W2.N2的系数%

 %             下面为地基中的第4列               %
A(N2+2,3*N2+4)=(1/2)*k1*(D(3,1)*(1-u(3,1)*u(3,1))+2*EI2/h2)+(5/2)*EI2*k1/h2;
%       t(2,0)=(1/2)*k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)+(5/2)*EI2*k1/h2%
%输入纸10,t2.0,W3.0的系数%
for j=N2+3:2*column-1
    A(j,j+2*N2+2)=k1*D(3,j-N2-1);
    %      t(2,j)=k1*D(2,j)%
    %输入论文54,ti.j,其中i=2,j=1,2,...,N2-1%
end
A(2*N2+2,4*N2+4)=(1/2)*k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)+(5/2)*EI4*k1/h2;
%        t(2,N2)=(1/2)*k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)+(5/2)*EI4*k1/h2%
%输入纸11,t2.N2,W3.N2的系数%

%             下面为地基中的第5列               %
A(2*N2+3,4*N2+5)=(1/2)*k1*(D(4,1)*(1-u(4,1)*u(4,1))+2*EI2/h2);
%         t(3,0)=(1/2)*k1*(D(3,0)*(1-u(3,0)*u(3,0))+2*EI2/h2)%
%输入纸12,t3.0,W4.0的系数%
for j=2*N2+4:3*N2+2
    A(j,j+2*N2+2)=k1*D(4,j-2*N2-2);
    %      t(3,j)=k1*D(3,j)%
    %输入论文54,ti.j,其中i=3,j=1,2,...,N2-1%
end
A(3*N2+3,5*N2+5)=(1/2)*k1*(D(4,column)*(1-u(4,column)*u(4,column))+2*EI4/h2);
%        t(3,N2)=(1/2)*k1*(D(3,N2)*(1-u(3,N2)*u(3,N2))+2*EI4/h2)%
%输入纸13,t3.N2,W4.N2的系数%

%             下面为地基中的第6列               %
A(3*N2+4,5*N2+6)=(1/2)*k1*(D(5,1)*(1-u(5,1)*u(5,1))+2*EI2/h2);
%         t(4,0)=(1/2)*k1*(D(4,0)*(1-u(4,0)*u(4,0))+2*EI2/h2)%
%输入纸14,t4.0,W5.0的系数%
for j=3*N2+5:4*N2+3
    A(j,j+2*N2+2)=k1*D(5,j-3*N2-3);
    %      t(4,j)=k1*D(4,j) %
    %输入论文54,ti.j,其中i=4,j=1,2,...,N2-1%
end
A(4*N2+4,6*N2+6)=(1/2)*k1*(D(5,column)*(1-u(5,column)*u(5,column))+2*EI4/h2);
%        t(4,N2)=(1/2)*k1*(D(4,N2)*(1-u(4,N2)*u(4,N2))+2*EI4/h2)%
%输入纸15,t4.N2,W5.N2的系数%

%             下面为地基中的第7列到第N1-1列  i=5,6,...,N1-3             %
for i=4*N2+5:column:column*(N1-4)+1
    A(i,i+2*N2+2)=(1/2)*k1*(D(i-bc6*N2+1,1)*(1-u(i-bc6*N2+1,1)*u(i-bc6*N2+1,1))+2*EI2/h2);
    %      t(i,0)=(1/2)*k1*(D(i,0)*(1-u(i,0)*u(i,0))+2*EI2/h2)%
    %输入纸16,tk.0,W(i+2).0的系数,k=i+1,即k=5,6,...,N1-3%
    for j=i+1:i+N2-1
        A(j,j+2*N2+2)=k1*D(i-bc6*N2+1,j-i+1);
        %      t(i,j)=k1*D(i,j)%
        %输入论文54,ti.j,其中i=k,j=1,2,...,N2-1%
    end
    A(i+N2,i+3*N2+2)=(1/2)*k1*(D(i-bc6*N2+1,column)*(1-u(i-bc6*N2+1,column)*u(i-bc6*N2+1,column))+2*EI4/h2);
    %        t(i,N2)=(1/2)*k1*(D(i,N2)*(1-u(i,N2)*u(i,N2))+2*EI4/h2)%