此为用MATLAB编写程序,进行杆系结构有限元分析的一个实例。代码综合考虑了两端铰接链杆单元,两端刚接梁单元、一端刚接一端铰接梁单元。编程理论基础为结构力学中的矩阵位移法。经过多个算例的分析计算,程序基本可靠合理,但仍有一些需要改进的地方。
1、当结构中含有自由端时,节点与单元编号该如何处理;
2、后期可考虑如何将计算结果进行可视化;
3、是否能将前期结构信息的输入以GUI交互方式进行。
待计算结构简图
节点编号与单元编码
源程序MATLAB代码:
有限元与变分法大作业1——MATLAB程序代码备注:结构的单元类型可包含两端铰接的链杆、两端刚接的梁单元、一端刚接一端铰接的梁单元。1、EXCEL表格输入结构信息计算之前在同目录文件下的名为“imform_list.xls”的EXCEL表格中写入结构相关信息,工作表中有四个工作簿,分别为node、element、force、constaint。2、定义参数并读入EXCEL表格信息%定义参数E=2.06e5;%弹性模量A=10;%截面面积I=100;%截面惯性矩[informlist_node]=xlsread('imform_list',1,'A2:I100') %从Excel表中读入节点信息。[informlist_element]=xlsread('imform_list',2) %从Excel表中读入单元信息,包括单元号,单元连接的节点,节点坐标。[informlist_force]=xlsread('imform_list',3) %从Excel表中读入荷载信息。[informlist_constraint]=xlsread('imform_list',4,'A2:A100') %从Excel表中读入有约束的位移列阵信息。[informlist_unconstraint]=xlsread('imform_list',4,'B2:B100') %从Excel表中读入没有约束的位移列阵信息。N_node=informlist_node(end,1)%节点总数N_element=informlist_element(end,1)%单元总数N_dis=max(max(informlist_node( : ,4:6)))%未知位移分量的总个数(包括边界约束处位移)%绘制结构,以检查数据的正确性。for i=1:N_element axis([min(min(informlist_node(1:end,2)))-0.2,... max(max(informlist_node(1:end,2)))+0.2,... min(min(informlist_node(1:end,3)))-0.2,... max(max(informlist_node(1:end,3)))+0.2]); line(informlist_element(i,[5,7]),informlist_element(i,[6,8]),'Color','k','LineStyle','-')end3、计算各单元转换矩阵%计算各单元转换矩阵T=cell(N_element,1); %定义一个元胞数组,用来存放转换矩阵。for i=1:N_elementT(i)={Transtionmatrix(i,informlist_element)};end4、计算各单元的刚度矩阵%计算各单元刚度矩阵(在整体坐标系下)k=cell(N_element,1); %首先定义一个元胞数组,用来装单元刚度矩阵。for i=1:N_elementk(i)={element_Stiffness(E,A,I,informlist_element,i)};end5、建立整体刚度方程%KK清零,然后调用函数Node _Assembly进行刚度矩阵的集成组装。KK=zeros(N_dis)for i=1:N_elementKK=element_Assembly(KK,k{i},informlist_element,i)end6、生成综合节点荷载列阵Pe=EqLode(informlist_node,informlist_element) %整体坐标系下的非结点荷载Pp=informlist_force %整体坐标系结点荷载P1=Pe+Pp%综合结点荷载7、边界条件的处理及刚度方程求解采用高斯消去法进行求解,注意:MATLAB中的反斜线符号“\”就是采用高斯消去法。K=KK(informlist_unconstraint,informlist_unconstraint)P=P1(informlist_unconstraint,1)U=K\P %位移UU=zeros(N_dis,1);UU(informlist_unconstraint)=Uq=zeros(N_node,3) %所有的节点位移for i=1:N_node if isnan(informlist_node(i,6))~=1 q(i,1:3)=UU(informlist_node(i,4:6)) else q(i,1:2)=UU(informlist_node(i,4:5)) endend8、位移变形图绘制%变形之后的节点坐标nodepos_xy=zeros(N_element,6)for i=1:N_element nodepos_xy(i,1:3)=q(informlist_element(i,3),1:3) nodepos_xy(i,4:6)=q(informlist_element(i,4),1:3)endnodepos_xy1=nodepos_xy(:,[1,2,4,5])*(0.3/max(max(nodepos_xy)))+... informlist_element(1:end,5:8) %变形之后的节点坐标列阵,其中变形位移放大了适当的倍数。%绘制变形位移图axis([min(min(informlist_node(1:end,2)))-1,max(max(informlist_node(1:end,2)))+... 1,min(min(informlist_node(1:end,3)))-1,max(max(informlist_node(1:end,3)))+1]);for i=1:N_elementline(nodepos_xy1(i,[1,3]),nodepos_xy1(i,[2,4]),'Color','r','LineStyle','--')hold online(informlist_element(i,[5,7]),informlist_element(i,[6,8]),'Color','k','LineStyle','-')end9、计算各杆的内力绘制注:各杆的杆端内力为局部坐标系下的。FF=zeros(N_element,6) %定义一个矩阵,用来存储各单元内力for i=1:N_element if informlist_element(i,2)==1 FF(i,1:6)=T{i}*k{i}*nodepos_xy(i,1:6)' else FF(i,[1,2,4,5])=T{i}*k{i}*nodepos_xy(i,[1,2,4,5])' end end函数定义部分function T=Transtionmatrix(elementnumber,informlist_element)%该函数计算单元的转换矩阵%输入单元信息列阵informent_element、单元编号elementnumber%输出单元转换矩阵if informlist_element(elementnumber,2)==1 L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2); C=(informlist_element(elementnumber,7)-informlist_element(elementnumber,5))/L; S=(informlist_element(elementnumber,8)-informlist_element(elementnumber,6))/L; T=[C S 0 0 0 0; -S C 0 0 0 0;0 0 1 0 0 0;0 0 0 C S 0;0 0 0 -S C 0;0 0 0 0 0 1];else L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2); C=(informlist_element(elementnumber,7)-informlist_element(elementnumber,5))/L; S=(informlist_element(elementnumber,8)-informlist_element(elementnumber,6))/L; T=[C S 0 0; -S C 0 0;0 0 C S;0 0 -S C];endendfunction k=element_Stiffness(E,A,I,informlist_element,elementnumber)%该函数计算单元在整体坐标系下的刚度矩阵,根据单元类型分别计算。%输入弹性模量E,横截面积A,单元信息列表informlist_element%输出单元刚度矩阵,如果为链杆单元,单元刚度矩阵为k(4X4),如果为梁单元,单元刚度矩阵为k(6X6)L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2);if informlist_element(elementnumber,2)==1 k=Transtionmatrix(elementnumber,informlist_element)'*... [E*A/L 0 0 -E*A/L 0 0;0 12*E*I/L^3 6*E*I/L^2 0 -12*E*I/L^3 6*E*I/L^2;... 0 6*E*I/L^2 4*E*I/L 0 -6*E*I/L^2 2*E*I/L;-E*A/L 0 0 E*A/L 0 0;... 0 -12*E*I/L^3 -6*E*I/L^2 0 12*E*I/L^3 -6*E*I/L^2;... 0 6*E*I/L^2 2*E*I/L 0 -6*E*I/L^2 4*E*I/L]*... Transtionmatrix(elementnumber,informlist_element);else k=E*A/L*Transtionmatrix(elementnumber,informlist_element)'*... [1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0]*... Transtionmatrix(elementnumber,informlist_element);endendfunction z=element_Assembly(KK,k,informlist_element,m)%该函数进行单元刚度矩阵的组装%输入单元刚度矩阵k%输出整体刚度矩阵KKif informlist_element(m,2)==1 for n1=1:6 for n2=1:6 KK(informlist_element(m,n1+8),informlist_element(m,n2+8))=... KK(informlist_element(m,n1+8),informlist_element(m,n2+8))+k(n1,n2); end endelsefor n1=1:4 for n2=1:4 dof(1)=9; dof(2)=10; dof(3)=12; dof(4)=13; KK(informlist_element(m,dof(n1)),informlist_element(m,dof(n2)))=... KK(informlist_element(m,dof(n1)),informlist_element(m,dof(n2)))+k(n1,n2); endendendz=KK;endfunction Pe=EqLode(informlist_node,informlist_element)%该函数进行结构非节点荷载的等效节点荷载计算。%输入单元信息列表,返回结构在整体坐标系下的等效节点荷载列阵。Pe=zeros(max(max(informlist_node( : ,4:6))),1);for i=1:informlist_element(end,1) P1=informlist_element(i,16); L1=sqrt((informlist_element(i,8)-informlist_element(i,6))^2+... (informlist_element(i,7)-informlist_element(i,5))^2); if informlist_element(i,2)==1 if informlist_element(i,15)==0 Pe(informlist_element(i,9:14))=0; elseif informlist_element(i,15)==1 Pe(informlist_element(i,9:14))=Pe(informlist_element(i,9:14))+... Transtionmatrix(i,informlist_element)'*... [0 P1/2 -P1*L1/8 0 P1/2 P1*L1/8]'; else Pe(informlist_element(i,9:14))=Pe(informlist_element(i,9:14))+... Transtionmatrix(i,informlist_element)'*... [0 P1*L1/2 -P1*L1^2/12 0 P1*L1/2 P1*L1^2/12]'; end else if informlist_element(i,15)==0 Pe(informlist_element(i,[9 10 12 13]))=0; elseif informlist_element(i,15)==1 Pe(informlist_element(i,[9 10 12 13]))=... Pe(informlist_element(i,[9 10 12 13]))+... Transtionmatrix(i,informlist_element)'*[0 P1/2 0 P1/2]'; else Pe(informlist_element(i,[9 10 12 13]))=... Pe(informlist_element(i,[9 10 12 13]))+... Transtionmatrix(i,informlist_element)'*[0 P1*L1/2 0 P1*L1/2]'; end endendend
详细的计算迭代过程如下所示:
有限元与变分法大作业1——MATLAB程序代码备注:结构的单元类型可包含两端铰接的链杆、两端刚接的梁单元、一端刚接一端铰接的梁单元。 1、EXCEL表格输入结构信息计算之前在同目录文件下的名为“imform_list.xls”的EXCEL表格中写入结构相关信息,工作表中有四个工作簿,分别为node、element、force、constaint。2、定义参数并读入EXCEL表格信息%定义参数E=2.06e5;%弹性模量A=10;%截面面积I=100;%截面惯性矩[informlist_node]=xlsread('imform_list',1,'A2:I100') informlist_node = 38×6 1 0 10 1 2 3 2 10 10 4 5 6 3 10 10 4 5 7 4 10 0 8 9 10 5 20 10 11 12 13 6 25 10 14 15 16 7 30 10 17 18 19 8 35 10 20 21 22 9 45 10 23 24 25 10 45 0 26 27 28%从Excel表中读入节点信息。[informlist_element]=xlsread('imform_list',2) informlist_element = 27×16 1 1 1 2 0 10 10 10 1 2 3 4 5 6 2 -10 2 1 4 3 10 0 10 10 8 9 10 4 5 7 0 0 3 1 2 5 10 10 20 10 4 5 6 11 12 13 2 -10 4 1 5 6 20 10 25 10 11 12 13 14 15 16 2 -10 5 1 6 7 25 10 30 10 14 15 16 17 18 19 2 -10 6 1 7 8 30 10 35 10 17 18 19 20 21 22 2 -10 7 1 8 9 35 10 45 10 20 21 22 23 24 25 2 -10 8 1 10 9 45 0 45 10 26 27 28 23 24 25 0 0 9 1 9 11 45 10 55 10 23 24 25 29 30 31 0 0 10 1 11 12 55 10 60 10 29 30 31 32 33 34 0 0%从Excel表中读入单元信息,包括单元号,单元连接的节点,节点坐标。[informlist_force]=xlsread('imform_list',3) informlist_force = 62×1 0 0 0 0 0 0 0 0 0 0%从Excel表中读入荷载信息。[informlist_constraint]=xlsread('imform_list',4,'A2:A100') informlist_constraint = 13×1 1 2 8 9 10 26 27 28 45 46%从Excel表中读入有约束的位移列阵信息。[informlist_unconstraint]=xlsread('imform_list',4,'B2:B100') informlist_unconstraint = 49×1 3 4 5 6 7 11 12 13 14 15%从Excel表中读入没有约束的位移列阵信息。N_node=informlist_node(end,1)%节点总数N_node = 38N_element=informlist_element(end,1)%单元总数N_element = 27N_dis=max(max(informlist_node( : ,4:6)))%未知位移分量的总个数(包括边界约束处位移)N_dis = 62%绘制结构,以检查数据的正确性。for i=1:N_element axis([min(min(informlist_node(1:end,2)))-0.2,... max(max(informlist_node(1:end,2)))+0.2,... min(min(informlist_node(1:end,3)))-0.2,... max(max(informlist_node(1:end,3)))+0.2]); line(informlist_element(i,[5,7]),informlist_element(i,[6,8]),'Color','k','LineStyle','-')end3、计算各单元转换矩阵%计算各单元转换矩阵T=cell(N_element,1) %定义一个元胞数组,用来存放转换矩阵。T = 27×1 cell 数组 {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}for i=1:N_elementT(i)={Transtionmatrix(i,informlist_element)};end4、计算各单元的刚度矩阵%计算各单元刚度矩阵(在整体坐标系下)k=cell(N_element,1) %首先定义一个元胞数组,用来装单元刚度矩阵。k = 27×1 cell 数组 {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}for i=1:N_elementk(i)={element_Stiffness(E,A,I,informlist_element,i)}endL = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 247200 0 -1236000 -247200 0 -1236000 0 206000 0 0 -206000 0 -1236000 0 8240000 1236000 0 4120000 -247200 0 1236000 247200 0 1236000 0 -206000 0 0 206000 0 -1236000 0 4120000 1236000 0 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 247200 0 -1236000 -247200 0 -1236000 0 206000 0 0 -206000 0 -1236000 0 8240000 1236000 0 4120000 -247200 0 1236000 247200 0 1236000 0 -206000 0 0 206000 0 -1236000 0 4120000 1236000 0 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 5k = 6×6 412000 0 0 -412000 0 0 0 1977600 4944000 0 -1977600 4944000 0 4944000 16480000 0 -4944000 8240000 -412000 0 0 412000 0 0 0 -1977600 -4944000 0 1977600 -4944000 0 4944000 8240000 0 -4944000 16480000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 247200 0 -1236000 -247200 0 -1236000 0 206000 0 0 -206000 0 -1236000 0 8240000 1236000 0 4120000 -247200 0 1236000 247200 0 1236000 0 -206000 0 0 206000 0 -1236000 0 4120000 1236000 0 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 10k = 6×6 206000 0 0 -206000 0 0 0 247200 1236000 0 -247200 1236000 0 1236000 8240000 0 -1236000 4120000 -206000 0 0 206000 0 0 0 -247200 -1236000 0 247200 -1236000 0 1236000 4120000 0 -1236000 8240000k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 29.1548k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 23.3238k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 17.4929k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 11.6619k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 11.6619k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 17.4929k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 23.3238k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 29.1548k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {0×0 double} {0×0 double} {0×0 double} {0×0 double}L = 6k = 6×6107 × 0.1144 0 -0.3433 -0.1144 0 -0.3433 0 0.0343 0 0 -0.0343 0 -0.3433 0 1.3733 0.3433 0 0.6867 -0.1144 0 0.3433 0.1144 0 0.3433 0 -0.0343 0 0 0.0343 0 -0.3433 0 0.6867 0.3433 0 1.3733k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {6×6 double} {0×0 double} {0×0 double} {0×0 double}L = 3k = 6×6107 × 0.9156 0 -1.3733 -0.9156 0 -1.3733 0 0.0687 0 0 -0.0687 0 -1.3733 0 2.7467 1.3733 0 1.3733 -0.9156 0 1.3733 0.9156 0 1.3733 0 -0.0687 0 0 0.0687 0 -1.3733 0 1.3733 1.3733 0 2.7467k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {6×6 double} {6×6 double} {0×0 double} {0×0 double}L = 3k = 6×6107 × 0.9156 0 -1.3733 -0.9156 0 -1.3733 0 0.0687 0 0 -0.0687 0 -1.3733 0 2.7467 1.3733 0 1.3733 -0.9156 0 1.3733 0.9156 0 1.3733 0 -0.0687 0 0 0.0687 0 -1.3733 0 1.3733 1.3733 0 2.7467k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {6×6 double} {6×6 double} {6×6 double} {0×0 double}L = 3k = 6×6107 × 0.9156 0 -1.3733 -0.9156 0 -1.3733 0 0.0687 0 0 -0.0687 0 -1.3733 0 2.7467 1.3733 0 1.3733 -0.9156 0 1.3733 0.9156 0 1.3733 0 -0.0687 0 0 0.0687 0 -1.3733 0 1.3733 1.3733 0 2.7467k = 27×1 cell 数组 {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {4×4 double} {6×6 double} {6×6 double} {6×6 double} {6×6 double}5、建立整体刚度方程%KK清零,然后调用函数Node _Assembly进行刚度矩阵的集成组装。KK=zeros(N_dis)KK = 62×62 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0for i=1:N_elementKK=element_Assembly(KK,k{i},informlist_element,i)endKK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 247200 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 -1236000 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 453200 0 0 1236000 -247200 0 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 453200 -1236000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 -1236000 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 62×62107 × 0.0206 0 0 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0247 0.1236 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0.0659 0 0 0.1236 -0.0247 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0700 0 0 0 -0.0206 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 1.6480 0 0 0 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 -0.1236 0 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 0 0 -0.1236 0.0247 0 -0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.4120 -0.1236 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0KK = 206000 0 0 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 247200 1236000 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 8240000 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 659200 0 0 1236000 -247200 0 1236000 -206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 -1236000 0 700400 0 0 0 -206000 0 0 -247200 1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 4120000 0 0 16480000 0 0 0 0 0 -1236000 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 8240000 -1236000 0 4120000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -247200 0 0 -1236000 247200 0 -1236000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -206000 0 0 0 206000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1236000 0 0 4120000 -1236000 0 8240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 06、生成综合节点荷载列阵Pe=EqLode(informlist_node,informlist_element) %整体坐标系下的非结点荷载Pe = 0 -50.0000 83.3333 0 -50.0000 0 0 0 0 0Pp=informlist_force %结点荷载Pp = 0 0 0 0 0 0 0 0 0 0P1=Pe+Pp%综合结点荷载P1 = 0 -50.0000 83.3333 0 -50.0000 0 0 0 0 07、边界条件的处理及刚度方程求解采用高斯消去法进行求解,注意:MATLAB中的反斜线符号“\”就是采用高斯消去法。K=KK(informlist_unconstraint,informlist_unconstraint)K = 0.8240 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0659 0 0 0.1236 -0.0206 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.1236 0 0.0700 0 0 0 -0.0247 0.1236 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.4120 0 0 1.6480 0 0 -0.1236 0.4120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0 0 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0206 0 0 0 0.0670 0.0031 0 -0.0412 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0247 -0.1236 0 0.0031 0.2244 0.3708 0 -0.1978 0.4944 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1236 0.4120 0 0 0.3708 2.4720 0 -0.4944 0.8240 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0412 0 0 0.0889 0.0039 0 -0.0412 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.1978 -0.4944 0.0039 0.3979 0 0 -0.1978 0.4944 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0P=P1(informlist_unconstraint,1)P = 83.3333 0 -50.0000 0 0 0 0 -62.5000 0 0U=K\P %位移U = -0.0053 0.0226 -0.1483 -0.0137 -0.0034 0.0520 -0.2181 -0.0053 0.0565 -0.2192UU=zeros(N_dis,1);UU(informlist_unconstraint)=UUU = 0 0 -0.0053 0.0226 -0.1483 -0.0137 -0.0034 0 0 0q=zeros(N_node,3) %所有的节点位移q = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0for i=1:N_node if isnan(informlist_node(i,6))~=1 q(i,1:3)=UU(informlist_node(i,4:6)) else q(i,1:2)=UU(informlist_node(i,4:5)) endendq = 0 0 -0.5274 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0 0 0 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0 0 0 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0 0 0 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0q = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0.0226 -0.1483 -0.0034 0 0 0 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 08、位移变形图绘制%变形之后的节点坐标nodepos_xy=zeros(N_element,6)nodepos_xy = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0for i=1:N_element nodepos_xy(i,1:3)=q(informlist_element(i,3),1:3) nodepos_xy(i,4:6)=q(informlist_element(i,4),1:3)endnodepos_xy = 0 0 -0.5274 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0 0 0 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0 0 0 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0 0 0 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0 0 0nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy = 0 0 -0.0053 0.0226 -0.1483 -0.0137 0 0 0 0.0226 -0.1483 -0.0034 0.0226 -0.1483 -0.0137 0.0520 -0.2181 -0.0053 0.0520 -0.2181 -0.0053 0.0565 -0.2192 0.0042 0.0565 -0.2192 0.0042 0.0514 -0.1817 0.0106 0.0514 -0.1817 0.0106 0.0390 -0.1146 0.0164 0.0390 -0.1146 0.0164 0.0103 -0.0346 0.0004 0 0 0 0.0103 -0.0346 0.0004 0.0103 -0.0346 0.0004 0.0084 -0.0588 -0.0043 0.0084 -0.0588 -0.0043 0.0078 -0.0813 -0.0044nodepos_xy1=nodepos_xy(:,[1,2,4,5])*(0.3/max(max(nodepos_xy)))+... informlist_element(1:end,5:8) %变形之后的节点坐标列阵,其中变形位移放大了适当的倍数。nodepos_xy1 = 0 10.0000 10.1201 9.2122 10.0000 0 10.1201 9.2122 10.1201 9.2122 20.2762 8.8411 20.2762 8.8411 25.3000 8.8357 25.3000 8.8357 30.2730 9.0348 30.2730 9.0348 35.2074 9.3911 35.2074 9.3911 45.0547 9.8163 45.0000 0 45.0547 9.8163 45.0547 9.8163 55.0444 9.6875 55.0444 9.6875 60.0416 9.5682%绘制变形位移图axis([min(min(informlist_node(1:end,2)))-1,max(max(informlist_node(1:end,2)))+... 1,min(min(informlist_node(1:end,3)))-1,max(max(informlist_node(1:end,3)))+1]);for i=1:N_elementline(nodepos_xy1(i,[1,3]),nodepos_xy1(i,[2,4]),'Color','r','LineStyle','--')hold online(informlist_element(i,[5,7]),informlist_element(i,[6,8]),'Color','k','LineStyle','-')end8、计算各杆的内力绘制注:各杆的杆端内力为局部坐标系下的。FF=zeros(N_element,6) %定义一个矩阵,用来存储各单元内力FF = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0for i=1:N_element if informlist_element(i,2)==1 FF(i,1:6)=T{i}*k{i}*nodepos_xy(i,1:6)' else FF(i,[1,2,4,5])=T{i}*k{i}*nodepos_xy(i,[1,2,4,5])' end endFF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0 0 0 0 0 0 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0 0 0 0 0 0FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336FF = -4.6561 13.1912 83.3333 4.6561 -13.1912 48.5786 30.5462 1.3968 13.9684 -30.5462 -1.3968 0.0000 -6.0530 -6.2626 -48.5786 6.0530 6.2626 -14.0470 -1.8474 -3.7392 -48.4530 1.8474 3.7392 29.7568 2.0951 -1.3737 -29.7568 -2.0951 1.3737 22.8881 5.0825 0.4187 -22.8881 -5.0825 -0.4187 24.9818 5.9247 0.9240 37.5182 -5.9247 -0.9240 -28.2782 7.1236 3.0270 14.3288 -7.1236 -3.0270 15.9417 0.3974 1.1746 15.5196 -0.3974 -1.1746 -3.7737 0.2193 1.2814 3.7737 -0.2193 -1.2814 2.6336函数定义部分function T=Transtionmatrix(elementnumber,informlist_element)%该函数计算单元的转换矩阵%输入单元信息列阵informent_element、单元编号elementnumber%输出单元转换矩阵if informlist_element(elementnumber,2)==1 L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2); C=(informlist_element(elementnumber,7)-informlist_element(elementnumber,5))/L; S=(informlist_element(elementnumber,8)-informlist_element(elementnumber,6))/L; T=[C S 0 0 0 0; -S C 0 0 0 0;0 0 1 0 0 0;0 0 0 C S 0;0 0 0 -S C 0;0 0 0 0 0 1];else L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2); C=(informlist_element(elementnumber,7)-informlist_element(elementnumber,5))/L; S=(informlist_element(elementnumber,8)-informlist_element(elementnumber,6))/L; T=[C S 0 0; -S C 0 0;0 0 C S;0 0 -S C];endendfunction k=element_Stiffness(E,A,I,informlist_element,elementnumber)%该函数计算单元在整体坐标系下的刚度矩阵,根据单元类型分别计算。%输入弹性模量E,横截面积A,单元信息列表informlist_element%输出单元刚度矩阵,如果为链杆单元,单元刚度矩阵为k(4X4),如果为梁单元,单元刚度矩阵为k(6X6)L=sqrt((informlist_element(elementnumber,8)-informlist_element(elementnumber,6))^2+... (informlist_element(elementnumber,7)-informlist_element(elementnumber,5))^2)if informlist_element(elementnumber,2)==1 k=Transtionmatrix(elementnumber,informlist_element)'*... [E*A/L 0 0 -E*A/L 0 0;0 12*E*I/L^3 6*E*I/L^2 0 -12*E*I/L^3 6*E*I/L^2;... 0 6*E*I/L^2 4*E*I/L 0 -6*E*I/L^2 2*E*I/L;-E*A/L 0 0 E*A/L 0 0;... 0 -12*E*I/L^3 -6*E*I/L^2 0 12*E*I/L^3 -6*E*I/L^2;... 0 6*E*I/L^2 2*E*I/L 0 -6*E*I/L^2 4*E*I/L]*... Transtionmatrix(elementnumber,informlist_element)else k=E*A/L*Transtionmatrix(elementnumber,informlist_element)'*... [1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0]*... Transtionmatrix(elementnumber,informlist_element);endendfunction z=element_Assembly(KK,k,informlist_element,m)%该函数进行单元刚度矩阵的组装%输入单元刚度矩阵k%输出整体刚度矩阵KKif informlist_element(m,2)==1 for n1=1:6 for n2=1:6 KK(informlist_element(m,n1+8),informlist_element(m,n2+8))=... KK(informlist_element(m,n1+8),informlist_element(m,n2+8))+k(n1,n2); end endelsefor n1=1:4 for n2=1:4 dof(1)=9; dof(2)=10; dof(3)=12; dof(4)=13; KK(informlist_element(m,dof(n1)),informlist_element(m,dof(n2)))=... KK(informlist_element(m,dof(n1)),informlist_element(m,dof(n2)))+k(n1,n2); endendendz=KK;endfunction Pe=EqLode(informlist_node,informlist_element)%该函数进行结构非节点荷载的等效节点荷载计算。%输入单元信息列表,返回结构在整体坐标系下的等效节点荷载列阵。Pe=zeros(max(max(informlist_node( : ,4:6))),1);for i=1:informlist_element(end,1) P1=informlist_element(i,16); L1=sqrt((informlist_element(i,8)-informlist_element(i,6))^2+... (informlist_element(i,7)-informlist_element(i,5))^2); if informlist_element(i,2)==1 if informlist_element(i,15)==0 Pe(informlist_element(i,9:14))=0; elseif informlist_element(i,15)==1 Pe(informlist_element(i,9:14))=Pe(informlist_element(i,9:14))+... Transtionmatrix(i,informlist_element)'*... [0 P1/2 -P1*L1/8 0 P1/2 P1*L1/8]'; else Pe(informlist_element(i,9:14))=Pe(informlist_element(i,9:14))+... Transtionmatrix(i,informlist_element)'*... [0 P1*L1/2 -P1*L1^2/12 0 P1*L1/2 P1*L1^2/12]'; end else if informlist_element(i,15)==0 Pe(informlist_element(i,[9 10 12 13]))=0; elseif informlist_element(i,15)==1 Pe(informlist_element(i,[9 10 12 13]))=... Pe(informlist_element(i,[9 10 12 13]))+... Transtionmatrix(i,informlist_element)'*[0 P1/2 0 P1/2]'; else Pe(informlist_element(i,[9 10 12 13]))=... Pe(informlist_element(i,[9 10 12 13]))+... Transtionmatrix(i,informlist_element)'*[0 P1*L1/2 0 P1*L1/2]'; end endendend
#在最美的年华里记录那些敲代码的日子#