Thực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCMThực hành Phương pháp phần tử hữu hạn ( Fem) so đồ 0 Đại học mở TPHCM
BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Sơ đồ tính q= 12kN/m P= 18kN C B A Mã Đề: 094 với số liệu hình học tải trọng sau: - B= 32cm H= 18cm L= 5.5m P= 18kN q= 12kN/m E= 2,1.104 kN/m2 Trường hợp giải tay: Bước 1: Thực chia kết cấu thành phần tử hệ thống nút đánh số 1,2,3 Như kết cấu có phần tử Mỗi nút có bậc tự thành phần chuyển vị vng góc với trục dầm góc xoay mặt phẳng dầm Chú ý rằng, bậc tự biết ta đánh số Cách làm xem ta áp đặt điều kiên biên từ đầu 0 q1 2 q2 Ma trận số [b] xác định: 0 (1) b (2) - Bước 2,3,4: Xác định ma trận cứng phần tử vectơ tải phần tử, lắp ghép thành ma trận cứng tổng thể vectơ tải tổng thể 0 12 L 12 L 4L L L2 EI K1 12 6L L3 dx L2 K 2 12 L 12 4L 6L EI 12 L dx SVTT: Hồ Xuân Diệu 6L L2 6L L2 MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Sau thực hiên ghép nối phần tử, ta có: EI L2 L2 L3 dx K* L2 L2 EI L dx Vectơ tải phần tử: + Trường hợp tải trọng tập trung P đặt cách nút i khoảng cách a P1 P3 P P4 P2 j i 3a 2a L2 L3 P1 2a a a P2 T L L3 P1 N ( x a ) xP xP P3 3a 2a P4 L2 L3 a a3 L L3 Với a= L/4, số vào phương trình ta được: 3( L / 4) 2( L / 4)3 L2 L3 2( L / 4) ( L / 4)3 ( L / 4) L L3 P1 xP 3( L / 4) 2( L / 4)3 L2 L3 ( L / 4) ( L / 4)3 L L3 + Trường hợp tải phân bố q P1 27 P 32 PL 64 P 32 PL 64 0 q P3 P4 P2 SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng 3a 2a qL L L 2 P1 2a a qL2 a P2 L L3 12 P qdx P3 qL 3a 2a L 2 P4 L L qL2 a a 12 L L3 Lăp ghép vectơ tải tổng thể, ý phải kể thêm lực tập trung nút Đối với tốn này, khơng có lực tập trung nút nên P * - nut qL2 PL 64 12 P* 2 qL 12 Bước 5: Giải hệ thống phương trình K * q* P* EI q1 L q2 8q1 2q2 qL2 L PL 64 12 EI 2q1 4q2 8q1 2q2 qL2 PL 64 12 qL 12 qL2 L 12 EI 12(5.5) L x18 x5.5 64 12 EI 12(5.5) L 12 EI 1639 L 8q1 2q2 64 EI 121 L 2q1 4q2 EI 2607 L q1 448 EI 9383 L q2 896 EI Đây góc xoay B C 2q1 4q2 SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn - GVHD: Th.S Trần Trung Dũng Bước 6: Tìm moment phần tử M S q EI L3 6L 6L L2 6L 2 L2 2L 6L 4L 11.64 23.28 2607 L 448 EI M S q EI L3 6L 6L 4L 2L 6L 2L 6L 4L 2607 L 448 EI 2.33 30.25 9383 L 896 EI Nhận xét: Từ giá trị trên, ta dễ dàng vẽ biểu đồ moment uốn kết cấu Nhận thấy rằng, biểu đồ moment uốn nhận từ FEM chưa xác với nghiệm thật moment phần tử (1) phải đường cong bậc Đồng thời, kiểm tra điều kiên cân moment nút khơng thõa mãn Điều hồn tồn hiểu momnet uốn giải theo FEM xét đến chuyển vị nút gây rầm chưa kể đến tải trọng tác dụng phần tử gây Do đó, để kết nội lực xác, ta phải cộng thêm momnet gây bời lực phần tử xem nút bị gắn cứng Kết biểu đồ nội lực: SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng 23.28 2.33 11.64 30.25 BIỂ U ĐỒ MOMENT GIẢI FEM 30.25 30.25 13.92 4.64 6.96 45.375 BIỂ U ĐỒ MOMENT DO TẢI TRÊN PHẦN TỬ 27.92 2.28 9.87 45.375 BIỂ U ĐỒ MOMENT SAU KHI HIỆU CHỈNH SVTT: Hồ Xn Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng 38.07 + 8.84 + 9.16 - 27.92 BIỂU ĐỒ LỰC CẮT Q Trường hợp giải MATLAB: disp('CHUONG TRINH 2DTRUSS -PHAN TICH KET CAU DAN PHANG - TAC GIA THAY THS.LE VAN BINH') disp % nhap du lieu ve ket cau P=18; q=12; b=0.32; h=0.18; 10 S=R*s; 11 % 2.nhap du lieu toa nut 12 toadoxy(1,1)=0; toadoxy(1,2)=0; 13 toadoxy(2,1)=0; toadoxy(2,2)=0.1763; 14 toadoxy(3,1)=1; toadoxy(3,2)=0.1763; 15 toadoxy(4,1)=0; toadoxy(4,2)=2.924; 16 % dac trung vat lieu va hinh hoc 17 for i=1:P 18 F(i)=0.01; 19 end 20 E=2e8; 21 % ket noi phan tu theo cac diem nut 22 ketnoiphantu(1,1)=1; ketnoiphantu(1,2)=3; 23 ketnoiphantu(2,1)=3; ketnoiphantu(2,2)=4; 24 ketnoiphantu(3,1)=2; ketnoiphantu(3,2)=3; 25 %5.dieu kien bien 26 dieukienbien(1)=1; 27 giatridieukienbien(1)=0; 28 dieukienbien(2)=2; 29 giatridieukienbien(2)=0; 30 dieukienbien(3)=3; 31 giatridieukienbien(3)=0; 32 dieukienbien(4)=4; 33 giatridieukienbien(4)=0; SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 GVHD: Th.S Trần Trung Dũng dieukienbien(5)=7; giatridieukienbien(7)=0; dieukienbien(6)=8; giatridieukienbien(8)=0; %6.dat mang luu tru vecto va ma tran bat dau la zero P=zeros(S,1); K=zeros(S,S); b=zeros(r*s,1); qe=zeros(r*s,1); Ke=zeros(r*s,r*s); Se=zeros(1,r*s); lucdoc=zeros(N,1); % nhap du lieu tai nut P(3)=5; P(5)=14.848; P(6)=1.7365; % 8.lap he thong phuong trinh de giai fem for i=1:N nd(1)=ketnoiphantu(i,1); nd(2)=ketnoiphantu(i,2); x1=toadoxy(nd(1),1); y1=toadoxy(nd(1),2); x2=toadoxy(nd(2),1); y2=toadoxy(nd(2),2); L=sqrt((x2-x1)^2+(y2-y1)^2); l=(x2-x1)/L; m=(y2-y1)/L; area=F(i); disp(['Chi so cua phan tu thu',num2str(i),':']) b=matranchiso(nd) disp(['Ma tran cung cua phan tu thu',num2str(i),':']) [Ke]=matrancungphantu(E,L,area,l,m) K=lapghep(K,Ke,b); end % 9.khu dieu kien bien va giai phuong trinh dai so [K,P]=khudieukienbien(K,P,dieukienbien,giatridieukienbien); q=K\P; % 10.tinh luc doc cac dan for iel=1:N nd(1)=ketnoiphantu(iel,1); nd(2)=ketnoiphantu(iel,2); x1=toadoxy(nd(1),1); y1=toadoxy(nd(1),2); x2=toadoxy(nd(2),1); y2=toadoxy(nd(2),2); L=sqrt((x2-x1)^2+(y2-y1)^2); l=(x2-x1)/L; m=(y2-y1)/L; area=F(iel); b=matranchiso(nd); SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 GVHD: Th.S Trần Trung Dũng for i=1:(r*s) qe(i)=q(b(i)); end Se=matrantinhlucdoc(E,area,L,l,m); for j=1:(r*s) lucdoc(iel)=lucdoc(iel)+Se(j)*qe(j); end end disp(['VEC TO CHUYEN VI NUT TONG THE:']) q disp(['LUC DOC TRONG CAC PHAN TU:']) lucdoc % 11.ve thi chuyen vi nut figure(1) scale=10000; for i=1:N nd(1)=ketnoiphantu(i,1); nd(2)=ketnoiphantu(i,2); b=matranchiso(nd); x1=toadoxy(nd(1),1); x1c=x1+scale*q(b(1)); y1=toadoxy(nd(1),2); y1c=y1+scale*q(b(2)); x2=toadoxy(nd(2),1); x2c=x2+scale*q(b(3)); y2=toadoxy(nd(2),2); y2c=y2+scale*q(b(4)); hold on axis equal plot([x1 x2],[y1 y2],'b'); end pause for i=1:N nd(1)=ketnoiphantu(i,1); nd(2)=ketnoiphantu(i,2); b=matranchiso(nd); x1=toadoxy(nd(1),1); x1c=x1+scale*q(b(1)); y1=toadoxy(nd(1),2); y1c=y1+scale*q(b(2)); x2=toadoxy(nd(2),1); x2c=x2+scale*q(b(3)); y2=toadoxy(nd(2),2); y2c=y2+scale*q(b(4)); hold on axis equal plot([x1c x2c],[y1c y2c],'r '); end disp('QUA TRINH GIAI HOAN THANH-KET THUC CHUONG TRINH') SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng 133 function [K,P]=khudieukienbien(K,P,dieukienbien,giatridieukienbien) 134 n=length(dieukienbien); 135 sdof=size(K); 136 for i=1:n 137 c=dieukienbien(i); 138 for j=1:sdof 139 K(c,j)=0; 140 K(j,c)=0; 141 end 142 K(c,c)=1; 143 P(c)=giatridieukienbien(i); 144 end 145 function [K]=lapghep(K,Ke,b) 146 edof=length(b); 147 for i=1:edof 148 ii=b(i); 149 for j=1:edof 150 jj=b(j); 151 K(ii,jj)=K(ii,jj)+Ke(i,j); 152 end 153 end 154 function [b]=matranchiso(nd) 155 b(1)=2*nd(1)-1; 156 b(2)=2*nd(1); 157 b(3)=2*nd(2)-1; 158 b(4)=2*nd(2); 159 function [Ke]=matrancungphantu(E,L,A,l,m) 160 Ke=(E*A/L)*[l*l l*m -l*l -l*m; 161 l*m m*m -l*m -m*m; 162 -l*l -l*m l*l l*m; 163 -l*m -m*m l*m m*m]; 164 function [Se]=matrantinhlucdoc(E,area,L,l,m) 165 Se=(E*area/L)*[-l -m l m]; 166 167 Kết quả: 168 CHUONG TRINH 2DTRUSS -PHAN TICH KET CAU DAN PHANG - TAC GIA -THAY THS.LE VAN BINH 169 Chi so cua phan tu thu1: 170 171 b = 172 173 174 175 Ma tran cung cua phan tu thu1: 176 177 Ke = 178 179 1.0e+006 * SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang BTL Phương pháp phần tử hữu hạn 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 1.9103 0.3368 -1.9103 -0.3368 0.3368 0.0594 -0.3368 -0.0594 -1.9103 -0.3368 1.9103 0.3368 GVHD: Th.S Trần Trung Dũng -0.3368 -0.0594 0.3368 0.0594 Chi so cua phan tu thu2: b= Ma tran cung cua phan tu thu2: Ke = 1.0e+005 * 0.8000 -2.1982 -0.8000 2.1982 -2.1982 6.0399 2.1982 -6.0399 -0.8000 2.1982 0.8000 -2.1982 2.1982 -6.0399 -2.1982 6.0399 Chi so cua phan tu thu3: b= Ma tran cung cua phan tu thu3: Ke = 2000000 0 -2000000 0 -2000000 0 2000000 0 0 VEC TO CHUYEN VI NUT TONG THE: q= 1.0e-005 * SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 10 BTL Phương pháp phần tử hữu hạn 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 GVHD: Th.S Trần Trung Dũng 0 0.3663 0.1972 0 LUC DOC TRONG CAC PHAN TU: lucdoc = 7.7800 -0.4105 7.3265 Bài 2: clc clear format compact format short non=3; noe=2; nonpe=2; nodofpn=3; nodofos=non*nodofpn; coord(1,1)=0; coord(1,2)=0; coord(2,1)=0; coord(2,2)=1; coord(3,1)=-1.2; coord(3,2)=1; E=2e8; A=0.01; I=0.000015; elem(1,1)=1; elem(1,2)=2; elem(2,1)=2; elem(2,2)=3; nores=3; ixres(1)=7; vodof(1)=0; ixres(2)=8; vodof(2)=0; ixres(3)=9; vodof(3)=0; SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 11 BTL Phương pháp phần tử hữu hạn 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 GVHD: Th.S Trần Trung Dũng K=zeros(nonpe*nodofpn,nonpe*nodofpn); ix=zeros(nonpe*nodofpn,1); KOS=zeros(nodofos,nodofos); f=zeros(nodofos,1); q=zeros(nonpe*nodofpn,1); Se=zeros(1,noe*nodofpn); noiluc=zeros(noe,1); f(1)=26.25; f(2)=6; f(3)=-6; f(4)=-11.25; f(5)=13.5; f(6)=-2.75; f(7)=0; f(8)=7.5; f(9)=1.25; for ie=1:noe endoe(1)=elem(ie,1); endoe(2)=elem(ie,2); x1=coord(endoe(1),1); y1=coord(endoe(1),2); x2=coord(endoe(2),1); y2=coord(endoe(2),2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0; alpha=2*atan(1); else alpha=4*atan(1); end disp(['Ma tran cung cua phan tu thu',num2str(ie),':']) K= smoframelem(E,A,L,I,alpha) disp(['Chi so cua phan tu thu',num2str(ie),':']) ix=indexos(endoe,nonpe,nodofpn) disp(['Ma tran cung tong the',]) KOS=smos(KOS,K,ix) end [KOS,f]=proores(KOS,f,ixres,vodof); KOS(:,1)=KOS(:,1)*0.866+KOS(:,2)*0.5; KOS(:,2)=[]; KOS(1,9)=0.5; KOS(2,9)=-0.866; u=KOS\f; id=1:1:nodofos; disp(['vector chuyen vi tong the']) disp=[id' u] function [ix]=indexos(endoe,nonpe,nodofpn) k=0; for i=1:nonpe s=(endoe(i)-1)*nodofpn; SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 12 BTL Phương pháp phần tử hữu hạn 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 GVHD: Th.S Trần Trung Dũng for j=1:nodofpn k=k+1; ix(k)=s+j; end end function [KOS,f]=proores(KOS,f,ixres,vodof) m=length(ixres); n=length(KOS); for i=1:m c=ixres(i); for j=1:n KOS(c,j)=0; KOS(j,c)=0; end KOS(c,c)=1; f(c)=vodof(i); end function [K]=smoframelem(E,A,L,I,alpha) k=E/L*[A 0 -A 0; 12*I/L^2 6*I/L -12*I/L^2 6*I/L; 6*I/L 4*I -6*I/L 2*I; -A 0 A 0; -12*I/L^2 -6*I/L 12*I/L^2 -6*I/L; 6*I/L 2*I -6*I/L 4*I]; c=cos(alpha); s=sin(alpha); T=[c s 0 0; -s c 0 0; 0 0 0; 0 c s 0; 0 -s c 0; 0 0 1]; K=T'*k*T; function [KOS]=smos(KOS,K,ix) for i=1:length(ix) is=ix(i); for j=1:length(ix) js=ix(j); KOS(is,js)=KOS(is,js)+K(i,j); end end Kết quả: Ma tran cung cua phan tu thu1: K= 1.0e+006 * 0.0360 0.0000 -0.0180 -0.0360 SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 13 BTL Phương pháp phần tử hữu hạn 0.0000 2.0000 0.0000 -0.0000 -0.0180 0.0000 0.0120 0.0180 -0.0360 -0.0000 0.0180 0.0360 GVHD: Th.S Trần Trung Dũng Chi so cua phan tu thu1: ix = Ma tran cung tong the KOS = 1.0e+006 * 0.0360 0.0000 -0.0180 -0.0360 0 0.0000 2.0000 0.0000 -0.0000 0 -0.0180 0.0000 0.0120 0.0180 0 -0.0360 -0.0000 0.0180 0.0360 0 Ma tran cung cua phan tu thu2: K= 1.0e+006 * 1.6667 -0.0000 -0.0000 -1.6667 -0.0000 0.0208 -0.0125 0.0000 -0.0000 -0.0125 0.0100 0.0000 -1.6667 0.0000 1.6667 0.0000 0.0000 -0.0208 0.0125 -0.0000 -0.0000 -0.0125 0.0050 0.0000 Chi so cua phan tu thu2: ix = Ma tran cung tong the SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 14 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng KOS = 1.0e+006 * 0.0360 0.0000 -0.0180 -0.0360 -0.0000 -0.0180 0.0000 2.0000 0.0000 -0.0000 -2.0000 0.0000 -0.0180 0.0000 0.0120 0.0180 -0.0000 0.0060 -0.0360 -0.0000 0.0180 1.7027 -0.0000 0.0180 -0.0000 -2.0000 -0.0000 -0.0000 -0.0180 0.0000 0.0060 2.0208 -0.0125 0.0180 -0.0125 0.0220 vector chuyen vi tong the disp = 1.0000 0000 2.0000 0000 3.0000 0000 4.0000 -0.00102 5.0000 0000 6.0 0.00175 Trường hợp giải phần mềm SAP2000: SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 15 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Định nghĩa vật liệu Định nghĩa tiết diện SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 16 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Vẽ dầm D32x18 Gán tải trọng SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 17 BTL Phương pháp phần tử hữu hạn SVTT: Hồ Xuân Diệu GVHD: Th.S Trần Trung Dũng MSSV: 0851020047 Trang 18 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Hoàn tất việc gán tải trọng cho tiết diên Gán điều kiện biên SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 19 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Chọn bậc tự RY Sau giải ( bấm F5 ) kết quả: SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 20 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Biểu đồ moment Biểu đồ lực cắt SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 21 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng Chuyển vị điểm B Chuyển vị điểm C SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 22 BTL Phương pháp phần tử hữu hạn GVHD: Th.S Trần Trung Dũng So Sánh Kết Quả: Kết chuyển vị nút điểm có lực tập trung góc xoay gối Nút có lực tập trung gối Giải tay Matlab SAP2000 Sai số (%) -9.7x10-4 (rad) -0.00102 (rad) -0.00102 (rad) 1.7x10-3 (rad) 0.00175 (rad) +0.00194 (rad) 12 Bảng kết moment đầu phần tử Phần tử Giải tay Matlab SAP2000 Sai số (%) -2.28 (kNm) -2.59 (kNm) -2.64 (kNm) 13 -27.92 (kNm) -30.24 (kNm) -32.56 (kNm) 14 0 0 SVTT: Hồ Xuân Diệu MSSV: 0851020047 Trang 23 ... 0. 000 0 -0. 01 80 0 .00 00 0 .01 20 0 .01 80 -0. 000 0 0. 006 0 -0. 03 60 -0. 000 0 0. 01 80 1. 702 7 -0. 000 0 0. 01 80 -0. 000 0 -2 .00 00 -0. 000 0 -0. 000 0 -0. 01 80 0 .00 00 0 .00 60 2 .02 08 -0. 0125 0. 01 80 -0. 0125 0. 02 20 vector chuyen... -0. 01 80 -0. 03 60 0 0. 000 0 2 .00 00 0 .00 00 -0. 000 0 0 -0. 01 80 0 .00 00 0 .01 20 0 .01 80 0 -0. 03 60 -0. 000 0 0. 01 80 0 .03 60 0 Ma tran cung cua phan tu thu2: K= 1.0e +00 6 * 1.6667 -0. 000 0 -0. 000 0 -1.6667 -0. 000 0... -0. 000 0 -1.6667 -0. 000 0 0. 0 208 -0. 0125 0. 000 0 -0. 000 0 -0. 0125 0. 0 100 0. 000 0 -1.6667 0. 000 0 1.6667 0. 000 0 0. 000 0 -0. 0 208 0. 0125 -0. 000 0 -0. 000 0 -0. 0125 0. 005 0 0. 000 0 Chi so cua phan tu thu2: ix