📄 源程序.txt
字号:
%输入纸17,tk.N2,W(i+2).N2的系数%
bc6=bc6+1;
end
% 下面为地基中的第N1列 %
A(column*(N1-3)+1,column*(N1-1)+1)=(1/2)*k1*(D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+2*EI2/h2)+(5/2)*EI2*k1/h2;
% t(N1-2,0)=(1/2)*k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)+(5/2)*EI2*k1/h2%
%输入纸18,t(N1-2).0,W(N1-1).0的系数%
for j=column*(N1-3)+2:column*(N1-2)-1
A(j,j+2*N2+2)=k1*D(N1-1,j-column*(N1-3));
% t(N1-2,j)=k1*D(N1-2,j)%
%输入论文54,ti.j,其中i=N1-2,j=1,2,...,N2-1%
end
A(column*(N1-2),column*N1)=(1/2)*k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)+(5/2)*EI4*k1/h2;
% t(N1-2,N2)=(1/2)*k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)+(5/2)*EI4*k1/h2%
%输入纸19,tN1-2.N2,W(N1-1).N2的系数%
% 下面为地基中的第N1+1列 %
A(column*(N1-2)+1,column*N1+1)=(1/2)*k1*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI2/h2)+4*EI2*k1/h2;
% t(N1-1,0)=(1/2)*k1*(D(N1-1,0)*(1-u(N1-1,0)*u(N1-1,0))+2*EI2/h2)+4*EI2*k1/h2%
%输入纸20,t(N1-1).0,WN1.0的系数%
for j=column*(N1-2)+2:column*(N1-1)-1
A(j,j+2*N2+2)=k1*D(N1,j-column*(N1-2));
% t(N1-1,j)=k1*D(N1-1,j)%
%输入论文54,ti.j,其中i=N1-1,j=1,2,...,N2-1%
end
A(column*(N1-1),column*row)=(1/2)*k1*(D(N1,column)*(1-u(N1,column)*u(N1,column))+2*EI4/h2)+4*EI4*k1/h2;
% t(N1-1,N2)=(1/2)*k1*(D(N1-1,N2)*(1-u(N1-1,N2)*u(N1-1,N2))+2*EI4/h2)+4*EI4*k1/h2%
%输入纸19,tN1-2.N2,W(N1-1).N2的系数%
%A矩阵上三角赋值结束%
A=A'+A;%A矩阵转置,并相加生成下三角%
%A对角线元素赋值%
% 输入A中的源数据表达式：对角线的元素
% 下面为地基中的第一列，第i列相对地基上的Wi.j所在列数，j=0，1，...，N2 %
A(1,1)=2*DL(1,1)*(1-ul(1,1))+1/2*k1*(D(2,1)*(1-u(2,1)*u(2,1))+2*EI2/h2)+1/2*k2*(D(1,2)*(1-u(1,2)*u(1,2))+2*EI1/h1)+2*EI2*k1/h2+2*EI1*k2/h1+1/4*C(1,1)*h12;
%a(0,0)=2*DL(1,1)*(1-ul(1,1))+1/2*k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)+1/2*k2*(D(0,1)*(1-u(0,1)*u(0,1))+2*EI1/h1)+2*EI2*k1/h2+2*EI1*k2/h1+1/4*C(0,0)*h1*h1*h2*h2 %
%输入纸1，a00,W00系数%
A(2,2)=k1*D(2,2)+2*DL(1,1)*(1-ul(1,1))+2*DL(1,2)*(1-ul(1,2))+1/2*k2*(D(1,3)*(1-u(1,3)*u(1,3))+2*EI1/h1)+2*k2*(D(1,2)*(1-u(1,2)*u(1,2))+2*EI1/h1)+25/2*EI1*k2/h1+1/2*C(1,2)*h12;
%a(0,1)=k1*D(1,1)+2*DL(1,1)*(1-ul(1,1))+2*DL(1,2)*(1-ul(1,2))+1/2*k2*(D(0,2)*(1-u(0,2)*u(0,2))+2*EI1/h1)+2*k2*(D(0,1)*(1-u(0,1)*u(0,1))+2*EI1/h1)+25/2*EI1*k2/h1+1/2*C(0,1)*h1*h1*h2*h2%
%输入纸2，a01,W01系数%
A(3,3)=k1*D(2,3)+2*DL(1,3)*(1-ul(1,3))+2*DL(1,2)*(1-ul(1,2))+1/2*k2*(D(1,4)*(1-u(1,4)*u(1,4))+2*EI1/h1)+1/2*k2*(D(1,2)*(1-u(1,2)*u(1,2))+2*EI1/h1)+2*k2*(D(1,3)*(1-u(1,3)*u(1,3))+2*EI1/h1)+8*EI1*k2/h1+1/2*C(1,3)*h12;
%a(0,2)=k1*D(1,2)+2*DL(1,3)*(1-ul(1,3))+2*DL(1,2)*(1-ul(1,2))+1/2*k2*(D(0,3)*(1-u(0,3)*u(0,3))+2*EI1/h1)+1/2*k2*(D(0,1)*(1-u(0,1)*u(0,1))+2*EI1/h1)+2*k2*(D(0,2)*(1-u(0,2)*u(0,2))+2*EI1/h1)+8*EI1*k2/h1+1/2*C(0,2)*h1*h1*h2*h2%
%输入纸3，a02,W02系数%
A(4,4)=k1*D(2,4)+2*(DL(1,4)*(1-ul(1,4))+DL(1,3)*(1-ul(1,3)))+k2*(1/2*D(1,5)*(1-u(1,5)*u(1,5))+1/2*D(1,3)*(1-u(1,3)*u(1,3))+2*D(1,4)*(1-u(1,4)*u(1,4))+6*EI1/h1)+1/2*(EI1*k2/h1+C(1,4)*h12);
%a(0,3)=k1*D(1,3)+2*DL(1,4)*(1-ul(1,4))+2*DL(1,3)*(1-ul(1,3))+1/2*k2*(D(0,4)*(1-u(0,4)*u(0,4))+2*EI1/h1)+1/2*k2*(D(0,2)*(1-u(0,2)*u(0,2))+2*EI1/h1)+2*k2*(D(0,3)*(1-u(0,3)*u(0,3))+2*EI1/h1)+1/2*EI1*k2/h1+1/2*C(0,3)*h1*h1*h2*h2%
%输入纸4，a03,W03系数%
%column为列节点数，通式%
for j=5:(N2-3)
A(j,j)=k1*D(2,j)+2*DL(1,j)*(1-ul(1,j))+2*DL(1,j-1)*(1-ul(1,j-1))+1/2*k2*(D(1,j+1)*(1-u(1,j+1)*u(1,j+1))+2*EI1/h1)+1/2*k2*(D(1,j-1)*(1-u(1,j-1)*u(1,j-1))+2*EI1/h1)+2*k2*(D(1,j)*(1-u(1,j)*u(1,j))+2*EI1/h1)+1/2*C(1,j)*h12;
%a(0,j)=k1*D(1,j)+2*DL(1,j+1)*(1-ul(1,j+1))+2*DL(1,j)*(1-ul(1,j))+1/2*k2*(D(0,j+1)*(1-u(0,j+1)*u(0,j+1))+2*EI1/h1)+1/2*k2*(D(0,j-1)*(1-u(0,j-1)*u(0,j-1))+2*EI1/h1)+2*k2*(D(0,j)*(1-u(0,j)*u(0,j))+2*EI1/h1)+1/2*C(0,j)*h1*h1*h2*h2 %
%输入纸5，a0j，W0j，j=4，...，N2-4的系数%
end
A((N2-2),(N2-2))=k1*D(2,N2-2)+2*DL(1,N2-2)*(1-ul(1,N2-2))+2*DL(1,N2-3)*(1-ul(1,N2-3))+1/2*k2*(D(1,N2-1)*(1-u(1,N2-1)*u(1,N2-1))+2*EI1/h1)+1/2*k2*(D(1,N2-3)*(1-u(1,N2-3)*u(1,N2-3))+2*EI1/h1)+2*k2*(D(1,N2-2)*(1-u(1,N2-2)*u(1,N2-2))+2*EI1/h1)+1/2*EI1*k2/h1+1/2*C(1,N2-2)*h12;
% a(0,N2-3)=k1*D(1,N2-3)+2*DL(1,N2-2)*(1-ul(1,N2-2))+2*DL(1,N2-3)*(1-ul(1,N2-3))+1/2*k2*(D(0,N2-2)*(1-u(0,N2-2)*u(0,N2-2))+2*EI1/h1)+1/2*k2*(D(0,N2-4)*(1-u(0,N2-4)*u(0,N2-4))+2*EI1/h1)+2*k2(D(0,N2-3)*(1-u(0,N2-3)*u(0,N2-3))+2*EI1/h1)+1/2*EI1*k2/h1+1/2*C(0,N2-3)*h1*h1*h2*h2 %
%输入纸6，a0(N2-3),W0(N2-3)系数%
%以下的没交换后两系数%
A((N2-1),(N2-1))=k1*D(2,N2-1)+2*DL(1,N2-1)*(1-ul(1,N2-1))+2*DL(1,N2-2)*(1-ul(1,N2-2))+1/2*k2*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI1/h1)+1/2*k2*(D(1,N2-2)*(1-u(1,N2-2)*u(1,N2-2))+2*EI1/h1)+2*k2*(D(1,N2-1)*(1-u(1,N2-1)*u(1,N2-1))+2*EI1/h1)+8*EI1*k2/h1+1/2*C(1,N2-1)*h12;
% a(0,N2-2)=k1*D(1,N2-2)+2*DL(1,N2-1)*(1-ul(1,N2-1))+2*DL(1,N2-2)*(1-ul(1,N2-2))+1/2*k2*(D(0,N2-1)*(1-u(0,N2-1)*u(0,N2-1))+2*EI1/h1)+1/2*k2*(D(0,N2-3)*(1-u(0,N2-3)*u(0,N2-3))+2*EI1/h1)+2*k2(D(0,N2-2)*(1-u(0,N2-2)*u(0,N2-2))+2*EI1/h1)+8*EI1*k2/h1+1/2*C(0,N2-2)*h1*h1*h2*h2 %
%输入纸7，a0(N2-2),W0(N2-2)系数%
A(N2,N2)=k1*D(2,N2)+2*DL(1,N2)*(1-ul(1,N2))+2*DL(1,N2-1)*(1-ul(1,N2-1))+1/2*k2*(D(1,N2-1)*(1-u(1,N2-1)*u(1,N2-1))+2*EI1/h1)+2*k2*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI1/h1)+25/2*EI1*k2/h1+1/2*C(1,N2)*h12;
%a(0,N2-1)=k1*D(1,N2-1)+2*DL(1,N2)*(1-ul(1,N2))+2*DL(1,N2-1)*(1-ul(1,N2-1))+1/2*k2*(D(0,N2-2)*(1-u(0,N2-2)*u(0,N2-2))+2*EI1/h1)+2*k2*(D(0,N2-1)*(1-u(0,N2-1)*u(0,N2-1))+2*EI1/h1)+25/2*EI1*k2/h1+1/2*C(0,N2-1)*h1*h1*h2*h2%
%输入纸8，a0(N2-1),W0(N2-1)系数%
%特殊点：%
A(column,column)=2*DL(1,N2)*(1-ul(1,N2))+1/2*k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)+1/2*k2*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI1/h1)+2*EI4*k1/h2+2*EI1*k2/h1+1/4*C(1,column)*h12;
% a(0,N2)=2*DL(1,N2)*(1-ul(1,N2))+1/2*k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)+1/2*k2*(D(0,N2-1)*(1-u(0,N2-1)*u(0,N2-1))+2*EI1/h1)+2*EI4*k1/h2+2*EI1*k2/h1+1/4*C(0,N2)*h1*h1*h2*h2%
%输入纸9，a0(N2),W0(N2)系数%
% 下面为地基中的第二列 %
A((N2+2),(N2+2))=k2*D(2,2)+2*DL(1,1)*(1-ul(1,1))+2*DL(2,1)*(1-ul(2,1))+1/2*k1*(D(3,1)*(1-u(3,1)*u(3,1))+2*EI2/h2)+2*k1*(D(2,1)*(1-u(2,1)*u(2,1))+2*EI2/h2)+25/2*EI2*k1/h2+1/2*C(2,1)*h12;
% a(1,0)=k2*D(1,1)+2*DL(1,1)*(1-ul(1,1))+2*DL(2,1)*(1-ul(2,1))+1/2*k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)+2*k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)+25/2*EI2*k1/h2+1/2*C(1,0)*h1*h1*h2*h2 %
%输入纸10，a10，W10的系数%
A((N2+3),(N2+3))=k1*D(3,2)+k2*D(2,3)+4*D(2,2)*(k1+k2+2*u(2,2))+2*DL(1,1)*(1-ul(1,1))+2*DL(1,2)*(1-ul(1,2))+2*DL(2,1)*(1-ul(2,1))+2*DL(2,2)*(1-ul(2,2))+C(2,2)*h12;
% a(1,1)=k1*D(2,1)+k2*D(1,2)+4*D(1,1)*(k1+k2+2*u(1,1))+2*DL(1,1)*(1-ul(1,1))+2*DL(1,2)*(1-ul(1,2))+2*DL(2,1)*(1-ul(2,1))+2*DL(2,2)*(1-ul(2,2))+C(1,1)*h1*h1*h2*h2%
%输入论文56，a11%
%N2为列节点末点的坐标，通式%
for j=(N2+4):(2*N2)
A(j,j)=k1*D(3,j-N2-1)+k2*(D(2,j-N2-2)+D(2,j-N2))+4*D(2,j-N2-1)*(k1+k2+2*u(2,j-N2-1))+2*DL(1,j-N2-2)*(1-ul(1,j-N2-2))+2*DL(1,j-N2-1)*(1-ul(1,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))+C(2,j-N2-1)*h12;
%a(1,j)=k1*D(2,j)+k2*(D(1,j-1)+D(1,j+1))+4*D(1,j)*(k1+k2+2*u(1,j))+2*DL(1,j)*(1-ul(1,j))+2*DL(1,j+1)*(1-ul(1,j+1))+2*DL(2,j)*(1-ul(2,j))+2*DL(2,j+1)*(1-ul(2,j+1))+C(1,j)*h1*h1*h2*h2%
%输入论文56，a1j，W0j,j=2,...,N2-2的系数%
%D(3,j)中j从3到N2-1,循环变量中的j与上面的j相差N2+1，同样对于DL，ul相差N2+2%
end
A((2*N2+1),(2*N2+1))=k1*D(3,N2)+k2*D(2,N2-1)+4*D(2,N2)*(k1+k2+2*u(2,N2))+2*DL(1,N2-1)*(1-ul(1,N2-1))+2*DL(1,N2)*(1-ul(1,N2))+2*DL(2,N2-1)*(1-ul(2,N2-1))+2*DL(2,N2)*(1-ul(2,N2))+C(2,N2)*h12;
% a(1,N2-1)=k1*D(2,N2-1)+k2*D(1,N2-2)+4*D(1,N2-1)*(k1+k2+2*u(1,N2-1))+2*DL(1,N2-1)*(1-ul(1,N2-1))+2*DL(1,N2)*(1-ul(1,N2))+2*DL(2,N2-1)*(1-ul(2,N2-1))+2*DL(2,N2)*(1-ul(2,N2))+C(1,N2-1)*h1*h1*h2*h2%
%输入论文56，a1N2-1%
A((2*N2+2),(2*N2+2))=k2*D(2,N2)+2*DL(1,N2)*(1-ul(1,N2))+2*DL(2,N2)*(1-ul(2,N2))+1/2*k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)+2*k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)+25/2*EI4*k1/h2+1/2*C(2,column)*h12;
% a(1,N2)=k2*D(1,N2-1)+2*DL(1,N2)*(1-ul(1,N2))+2*DL(2,N2)*(1-ul(2,N2))+1/2*k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)+2*k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)+25/2*EI4*k1/h2+1/2*C(1,N2)*h1*h1*h2*h2%
%输入纸11，a1N2，W1N2的系数%
% 下面为地基中的第三列 %
A((2*N2+3),(2*N2+3))=k2*D(3,2)+2*DL(2,1)*(1-ul(2,1))+2*DL(3,1)*(1-ul(3,1))+1/2*k1*(D(4,1)*(1-u(4,1)*u(4,1))+2*EI2/h2)+2*k1*(D(3,1)*(1-u(3,1)*u(3,1))+2*EI2/h2)+1/2*k1*(D(2,1)*(1-u(2,1)*u(2,1))+2*EI2/h2)+8*EI2*k1/h2+1/2*C(3,1)*h12;
% a(2,0)=k2*D(2,1)+2*DL(2,1)*(1-ul(2,1))+2*DL(3,1)*(1-ul(3,1))+1/2*k1*(D(3,0)*(1-u(3,0)*u(3,0))+2*EI2/h2)+2*k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)+1/2*k1*(D(1,0)*(1-u(1,0)*u(1,0))+2*EI2/h2)+8*EI2*k1/h2+1/2*C(2,0)*h1*h1*h2*h2 %
%输入纸12，a20，W20的系数%
A((2*N2+4),(2*N2+4))=k1*(D(2,2)+D(4,2))+k2*D(3,3)+4*D(3,2)*(k1+k2+2*u(3,2))+2*(DL(2,1)*(1-ul(2,1))+DL(2,2)*(1-ul(2,2))+DL(3,1)*(1-ul(3,1))+DL(3,2)*(1-ul(3,2)))+C(3,2)*h12;
% a(2,1)=k1*(D(1,1)+D(3,1))+k2*D(2,2)+4*D(2,1)*(k1+k2+2*u(2,1))+2*DL(2,1)*(1-ul(2,1))+2*DL(2,2)*(1-ul(2,2))+2*DL(3,1)*(1-ul(3,1))+2*DL(3,2)*(1-ul(3,2))+C(2,1)*h1*h1*h2*h2%
%输入论文57，ai1,i=2%
%N2为列节点末点的坐标，通式%
for j=(2*N2+5):(3*N2+1)
A(j,j)=k1*(D(2,j-2*N2-2)+D(4,j-2*N2-2))+k2*(D(3,j-2*N2-3)+D(3,j-2*N2-1))+4*D(3,j-2*N2-2)*(k1+k2+2*u(3,j-2*N2-2))+2*DL(2,j-2*N2-3)*(1-ul(2,j-2*N2-3))+2*DL(2,j-2*N2-2)*(1-ul(2,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))+C(3,j-2*N2-2)*h12;
%a(2,j)=k1*(D(1,j)+D(3,j)) +k2*(D(2,j-1) +D(2,j+1)) +4*D(2,j)*(k1+k2+2*u(2,j)) +2*DL(2,j)*(1-ul(2,j)) +2*DL(2,j+1)*(1-ul(2,j+1)) +2*DL(3,j)*(1-ul(3,j)) +2*DL(3,j+1)*(1-ul(3,j+1))+C(2,j)*h1*h1*h2*h2%
%输入论文57，aij，i=2,j=2,...,N2-2的系数%
end
A((3*N2+2),(3*N2+2))=k1*(D(2,N2)+D(4,N2))+k2*D(3,N2-1)+4*D(3,N2)*(k1+k2+2*u(3,N2))+2*DL(2,N2-1)*(1-ul(2,N2-1))+2*DL(2,N2)*(1-ul(2,N2))+2*DL(3,N2-1)*(1-ul(3,N2-1))+2*DL(3,N2)*(1-ul(3,N2))+C(3,N2)*h12;
% a(2,N2-1)=k1*(D(1,N2-1)+D(3,N2-1))+k2*D(2,N2-2)+4*D(2,N2-1)*(k1+k2+2*u(2,N2-1))+2*DL(2,N2-1)*(1-ul(2,N2-1))+2*DL(2,N2)*(1-ul(2,N2))+2*DL(3,N2-1)*(1-ul(3,N2-1))+2*DL(3,N2)*(1-ul(3,N2))+C(2,N2-1)*h1*h1*h2*h2 %
%输入论文57，aiN2-1,i=2%
A((3*N2+3),(3*N2+3))=k2*D(3,N2)+2*DL(2,N2)*(1-ul(2,N2))+2*DL(3,N2)*(1-ul(3,N2))+1/2*k1*(D(4,column)*(1-u(4,column)*u(4,column))+2*EI4/h2)+2*k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)+1/2*k1*(D(2,column)*(1-u(2,column)*u(2,column))+2*EI4/h2)+8*EI4*k1/h2+1/2*C(3,column)*h12;
% a(2,N2)=k2*D(2,N2-1)+2*DL(2,N2)*(1-ul(2,N2))+2*DL(3,N2)*(1-ul(3,N2))+1/2*k1*(D(3,N2)*(1-u(3,N2)*u(3,N2))+2*EI4/h2)+2*k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)+1/2*k1*(D(1,N2)*(1-u(1,N2)*u(1,N2))+2*EI4/h2)+8*EI4*k1/h2+1/2*C(2,N2)*h1*h1*h2*h2 %
%输入纸13，a2N2，W2N2的系数%
% 下面为地基中的第四列 %
A((3*N2+4),(3*N2+4))=k2*D(4,2)+2*DL(4,1)*(1-ul(4,1))+2*DL(3,1)*(1-ul(3,1))+1/2*k1*(D(5,1)*(1-u(5,1)*u(5,1))+2*EI2/h2)+1/2*k1*(D(3,1)*(1-u(3,1)*u(3,1))+2*EI2/h2)+2*k1*(D(4,1)*(1-u(4,1)*u(4,1))+2*EI2/h2)+1/2*EI2*k1/h2+1/2*C(4,1)*h12;
% a(3,0)=k2*D(3,1)+2*DL(4,1)*(1-ul(4,1))+2*DL(3,1)*(1-ul(3,1))+1/2*k1*(D(4,0)*(1-u(4,0)*u(4,0))+2*EI2/h2)+1/2*k1*(D(2,0)*(1-u(2,0)*u(2,0))+2*EI2/h2)+2*k1*(D(3,0)*(1-u(3,0)*u(3,0))+2*EI2/h2)+1/2*EI2*k1/h2+1/2*C(3,0)*h1*h1*h2*h2 %
%输入纸14，a30，W30的系数%
A((3*N2+5),(3*N2+5))=k1*(D(3,2)+D(5,2))+k2*D(4,3)+4*D(4,2)*(k1+k2+2*u(4,2))+2*DL(3,1)*(1-ul(3,1))+2*DL(3,2)*(1-ul(3,2))+2*DL(4,1)*(1-ul(4,1))+2*DL(4,2)*(1-ul(4,2))+C(4,2)*h12;
% a(3,1)=k1*(D(2,1)+D(4,1))+k2*D(3,2)+4*D(3,1)*(k1+k2+2*u(3,1))+2*DL(3,1)*(1-ul(3,1))+2*DL(3,2)*(1-ul(3,2))+2*DL(4,1)*(1-ul(4,1))+2*DL(4,2)*(1-ul(4,2))+C(3,1)*h1*h1*h2*h2%
%输入论文57，ai1,i=3%
%N2为列节点末点的坐标，通式%
for j=(3*N2+6):(4*N2+2)
A(j,j)=k1*(D(3,j-3*N2-3)+D(5,j-3*N2-3))+k2*(D(4,j-3*N2-4)+D(4,j-3*N2-2))+4*D(4,j-3*N2-3)*(k1+k2+2*u(4,j-3*N2-3))+2*DL(3,j-3*N2-4)*(1-ul(3,j-3*N2-4))+2*DL(3,j-3*N2-3)*(1-ul(3,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))+C(4,j-3*N2-3)*h12;
%a(3,j)=k1*(D(2,j)+D(4,j))+k2*(D(3,j-1)+D(3,j+1))+4*D(3,j)*(k1+k2+2*u(3,j))+2*DL(3,j)*(1-ul(3,j))+2*DL(3,j+1)*(1-ul(3,j+1))+2*DL(4,j)*(1-ul(4,j))+2*DL(4,j+1)*(1-ul(4,j+1))+C(3,j)*h1*h1*h2*h2 %
%输入论文57，aij，i=3,j=2,...,N2-2的系数%
end
A((4*N2+3),(4*N2+3))=k1*(D(3,N2)+D(5,N2))+k2*D(4,N2-1)+4*D(4,N2)*(k1+k2+2*u(4,N2))+2*(DL(3,N2-1)*(1-ul(3,N2-1))+DL(3,N2)*(1-ul(3,N2))+DL(4,N2-1)*(1-ul(4,N2-1))+DL(4,N2)*(1-ul(4,N2)))+C(4,N2)*h12;
% a(3,N2-1)=k1*(D(2,N2-1)+D(4,N2-1))+k2*D(3,N2-2)+4*D(3,N2-1)*(k1+k2+2*u(3,N2-1))+2*DL(3,N2-1)*(1-ul(3,N2-1))+2*DL(3,N2)*(1-ul(3,N2))+2*DL(4,N2-1)*(1-ul(4,N2-1))+2*DL(4,N2)*(1-ul(4,N2))+C(3,N2-1)*h1*h1*h2*h2 %
%输入论文57，aiN2-1,i=3%
A((4*N2+4),(4*N2+4))=k2*D(4,N2)+2*DL(4,N2)*(1-ul(4,N2))+2*DL(3,N2)*(1-ul(3,N2))+1/2*k1*(D(5,column)*(1-u(5,column)*u(5,column))+2*EI4/h2)+1/2*k1*(D(3,column)*(1-u(3,column)*u(3,column))+2*EI4/h2)+2*k1*(D(4,column)*(1-u(4,column)*u(4,column))+2*EI4/h2)+1/2*EI4*k1/h2+1/2*C(4,column)*h12;
% a(3,N2)=k2*D(3,N2-1)+2*DL(4,N2)*(1-ul(4,N2))+2*DL(3,N2)*(1-ul(3,N2))+1/2*k1*(D(4,N2)*(1-u(4,N2)*u(4,N2))+2*EI4/h2)+1/2*k1*(D(2,N2)*(1-u(2,N2)*u(2,N2))+2*EI4/h2)+2*k1*(D(3,N2)*(1-u(3,N2)*u(3,N2))+2*EI4/h2)+1/2*EI4*k1/h2+1/2*C(3,N2)*h1*h1*h2*h2 %
%输入纸15，a3N2，W3N2的系数%
% 下面为地基中的第五列到N1-3列 i=4，....N1-4 %
for i=(4*N2+5):column:((N1-4)*N2+N1-3)
A(i,i)=k2*D(i-bc0*N2,2)+2*DL(i-bc0*N2,1)*(1-ul(i-bc0*N2,1))+2*DL(i-bc0*N2-1,1)*(1-ul(i-bc0*N2-1,1))+1/2*k1*(D(i-bc0*N2+1,1)*(1-u(i-bc0*N2+1,1)*u(i-bc0*N2+1,1))+2*EI2/h2)+1/2*k1*(D(i-bc0*N2-1,1)*(1-u(i-bc0*N2-1,1)*u(i-bc0*N2-1,1))+2*EI2/h2)+2*k1*(D(i-bc0*N2,1)*(1-u(i-bc0*N2,1)*u(i-bc0*N2,1))+2*EI2/h2)+1/2*C(i-bc0*N2,1)*h12;
%a(i,0)=k2*D(i,1)+2*DL(i+1,1)*(1-ul(i+1,1))+2*DL(i,1)*(1-ul(i,1))+1/2*k1*(D(i+1,0)*(1-u(i+1,0)*u(i+1,0))+2*EI2/h2)+1/2*k1*(D(i-1,0)*(1-u(i-1,0)*u(i-1,0))+2*EI2/h2)+2*k1*(D(i,0)*(1-u(i,0)*u(i,0))+2*EI2/h2)+1/2*C(i,0)*h1*h1*h2*h2%
%输入纸16,ai0论文57%
%i从5到N1-3%
A(i+1,i+1)=k1*(D(i-bc0*N2-1,2)+D(i-bc0*N2+1,2))+k2*D(i-bc0*N2,3)+4*D(i-bc0*N2,2)*(k1+k2+2*u(i-bc0*N2,2))+2*DL(i-bc0*N2-1,1)*(1-ul(i-bc0*N2-1,1))+2*DL(i-bc0*N2-1,2)*(1-ul(i-bc0*N2-1,2))+2*DL(i-bc0*N2,1)*(1-ul(i-bc0*N2,1))+2*DL(i-bc0*N2,2)*(1-ul(i-bc0*N2,2))+C(i-bc0*N2,2)*h12;
% a(i,1)=k1*(D(i-1,1)+D(i+1,1))+k2*D(i,2)+4*D(i,1)*(k1+k2+2*u(i,1))+ 2*DL(i,1)*(1-ul(i,1))+2*DL(i,2)*(1-ul(i,2))+2*DL(i+1,1)*(1-ul(i+1,1))+2*DL(i+1,2)*(1-ul(i+1,2))+C(i,1)*h1*h1*h2*h2%
%输入论文57，ai1%
for j=i+2:i+N2-2
A(j,j)=k1*(D(i-bc0*N2-1,j-i+1)+D(i-bc0*N2+1,j-i+1))+k2*(D(i-bc0*N2,j-i)+D(i-bc0*N2,j-i+2))+4*D(i-bc0*N2,j-i+1)*(k1+k2+2*u(i-bc0*N2,j-i+1))+2*(DL(i-bc0*N2-1,j-i)*(1-ul(i-bc0*N2-1,j-i))+DL(i-bc0*N2-1,j-i+1)*(1-ul(i-bc0*N2-1,j-i+1))+DL(i-bc0*N2,j-i)*(1-ul(i-bc0*N2,j-i))+DL(i-bc0*N2,j-i+1)*(1-ul(i-bc0*N2,j-i+1)))+C(i-bc0*N2,j-i+1)*h12;
%a(i,j)=k1*(D(i-1,j)+D(i+1,j))+k2*(D(i,j-1)+D(i,j+1))+4*D(i,j)*(k1+k2+2*u(i,j))+ 2*DL(i,j)*(1-ul(i,j))+2*DL(i,j+1)*(1-ul(i,j+1))+2*DL(i+1,j)*(1-ul(i+1,j))+2*DL(i+1,j+1)*(1-ul(i+1,j+1))+C(i,j)*h1*h1*h2*h2%
%输入论文57,aij，j=2,...,N2-2%
end
A(i+N2-1,i+N2-1)=k1*(D(i-bc0*N2-1,N2)+D(i-bc0*N2+1,N2))+k2*D(i-bc0*N2,N2-1)+4*D(i-bc0*N2,N2)*(k1+k2+2*u(i-bc0*N2,N2))+2*DL(i-bc0*N2-1,N2-1)*(1-ul(i-bc0*N2-1,N2-1))+2*DL(i-bc0*N2-1,N2)*(1-ul(i-bc0*N2-1,N2))+2*DL(i-bc0*N2,N2-1)*(1-ul(i-bc0*N2,N2-1))+2*DL(i-bc0*N2,N2)*(1-ul(i-bc0*N2,N2))+C(i-bc0*N2,N2)*h12;
% a(i,N2-1)=k1*(D(i-1,N2-1)+D(i+1,N2-1))+k2*D(i,N2-2) +4*D(i,N2-1)*(k1+k2+2*u(i,N2-1)) +2*DL(i,N2-1)*(1-ul(i,N2-1))+2*DL(i,N2)*(1-ul(i,N2))+2*DL(i+1,N2-1)*(1-ul(i+1,N2-1))+2*DL(i+1,N2)*(1-ul(i+1,N2))+C(i,N2-1)*h1*h1*h2*h2%
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -