此為用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
#在最美的年華裡記錄那些敲代碼的日子#