Nghiên cứu phương pháp dãi hữu hạn và ứng dụng để khảo sát dao động của tấm có sườn

BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH ĐẠI HỌC BÁCH KHOA TP.HỒ CHÍ MINH CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập – Tự Do – Hạnh Phúc -oOo - Nhiệm vụ luận án cao học Họ tên học viên Ngày, tháng, năm sinh Chuyên ngành : Lê Hiền Anh : 13-08-1975 : Xây Dựng DD & CN Phái Nơi sinh Khoá : nam : TP Đà Nẵng : 10 Tên đề tài: Lập trình tính toán kết cấu phương pháp dải hữu hạn ứng dụng việc tính dao động có sườn Nhiệm vụ nội dung:  Nghiên cứu sở lý thuyết phương pháp dải hữu hạn  Xây dựng phần mềm tính toán  ng dụng khảo sát ảnh hưởng đại lượng hình học có sườn đến dao động  Kết luận kiến nghị phương pháp dải hữu hạn Ngày giao nhiệm vụ: Ngày hoàn thành nhiệm vụ: Họ tên cán hướng dẫn: PGS.TS Chu Quốc Thắng Họ tên cán phản biện 1: Họ tên cán phản biện 2: Cán hướng dẫn Cán phản biện Cán phản biện PGS.TS Chu Quốc Thắng Nội dung đề cương luận án cao học thông qua hội đồng chuyên ngành Phòng Quản Lý Khoa Học – Sau Đại Học Ngày tháng năm 2003 Chủ nhiệm ngành PGS.TS Chu Quốc Thắng CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH Nhận xét cán hướng dẫn :PGS.TS Chu Quốc Thắng Nhận xét cán phản biện 1: Nhận xét cán phản biện 2: Luận án cao học bảo vệ Hội Đồng Bảo Vệ Luận n Cao Học trường Đại Học Bách Khoa Thành Phố Hồ Chí Minh, ngày tháng năm 2003 Có thể tìm hiểu luận án thư viện trường đại học Bách Khoa Thành Phố Hồ Chí Minh LỜI CẢM ƠN Xin chân thành cảm ơn cha mẹ tôi, người thân thương tôi, sinh thành nuôi dưỡng đến ngày hôm Xin cảm ơn PGS TS Chu Quốc Thắng hướng dẫn tận tình cho hoàn thành luận văn tốt nghiệp Xin cảm ơn bạn bè; đồng nghiệp động viên chia với khó khăn lúc thực luận văn Lê Hiền Anh TÓM TẮT Phương pháp dải hữu hạn (FSM) phương pháp phát triển từ phương pháp phần tử hữu hạn (FEM) áp dụng riêng cho kết cấu có hình học đặc biệt : kết cấu có hình dạng thẳng, tiết diện không đổi chịu tải trọng tương đối đơn giản … FSM nghiên cứu phát triển từ cuối năm 1960 đến phát triển Tuy nhiên Việt Nam phương pháp chưa phổ biến Nghiên cứu đề tài học viên Trần Văn Bình thực đề cập đến vấn đề tính toán chuyển vị ứng suất toán ứng suất phẳng Luận văn tác giả nghiên cứu nhằm giải toán: Tìm chuyển vị; ứng suất đối toán kết cấu chịu uốn ứng suất phẳng đồng thời Tần số dao động riêng dạng dao động tương ứng Tìm tải trọng tới hạn dạng ổn định Luận văn gồm chương: Chương - Chương tổng quan : Chương – Cơ sở lý thuyết phương pháp Chương – Ví dụ minh hoạ Chương – Khảo sát ảnh hưởng sườn gia cường đến dao động chử nhật Chương – kết luận; kiến nghị Một số ký hiệu viết tắt  a : khoảng cách sườn gia cường  b : bề rộng phần tử dải  c : tham số hàm chuyển vị tổng quát  f : hàm chuyển vị tổng quát  h : chiều dày phần tử  i : biến đếm; số phức  j : biến đếm  k : hệ số động  l : chiều dài khoảng nút  m : số đường nút  n : số phần tử  p : tải phân bố đường  q : tải phân bố diện tích  r : chuyển vị góc quanh trục y (chuyển vị đại diện cho đường nút)  s : số thứ tự chuổi xét  t : thời gian  u : chuyển vị phương x  v : chuyển vị phương y  w : chuyển vị phương z  x : tọa độ địa phương  y : tọa độ địa phương  z : tọa độ địa phương  Ai : diện tích mặt cắt phần tử lăng trụ thứ i  A : diện tích lấy tích phân  B : ma trận tính biến dạng  C : ma trận hệ số cản nhớt  D : ma trận tính chất vật liệu  E : mô dun đàn hồi Young  F : hàm dạng tổng quát  G : ma trận hiệu chỉnh (khi giải toán dao động cưỡng bức)  K : ma trận độ cứng  KG : ma trận độ cứng hình học  L : chiều dài phần tử  M : ma trận khối lượng  N : hàm nội suy phương x  P : tải tập trung  Q : hàm tải theo thới gian  S : tổng số chuổi cần tính toán  T : ma trận chuyển trục  U : lượng biến dạng  V : thể tích lấy tích phân  X : chuổi hàm nội suy phương x (bài toán lớp hữu hạn)  Y : chuổi hàm nội suy phương y   góc hệ toạ độ địa phương hệ toạ độ tổng thể   tỷ số tần số lực cưỡng tần số riêng   độ biến thiên chuyển vị góc (thành phần thêm vào phần tử bậc cao HO2)   vectochuyển vị đại diện   vecto biến dạng   vecto riêng mô tả mode dao động  xy thành phần biến dạng góc = 2xy   : hệ số xác định tải trọng tới hạn   hệ số chu kỳ chuổi   hệ số nở hông   khối lượng riêng   vecto ứng suất   biến thời gian (biến tạm dấu tích phân)   : tham số đặc tính chuổi (trong công thức chuổi hàm phương y)   tần số riêng     thế biến dạng toàn phần   : tần số tương đương   chuyển vị tương đương độ biến thiên Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG CHƯƠNG 1: TỔNG QUAN 1 GIỚI THIỆU PHƯƠNG PHÁP DẢI HỮU HẠN (FINITE STRIP METHOD – FSM) QUÁ TRÌNH PHÁT TRIỂN CỦA FSM .2 CHƯƠNG 2: CƠ SỞ LÝ THUYẾT CỦA PHƯƠNG PHÁP .5 PHẦN TỬ DẢI TRONG HỆ TOẠ ĐỘ DECAC PHƯƠNG PHÁP LỚP HỮU HẠN 19 PHƯƠNG PHÁP LĂNG TRỤ HỮU HẠN 20 PHÂN TÍCH DAO ĐỘNG 22 PHÂN TÍCH ỔN ĐỊNH BẰNG FSM 23 BÀI TOÁN ĐỘNG LỰC HỌC GIẢI BẰNG FSM 27 FSM viết tronh hệ toạ độ trụ 32 NHẬN XÉT PHƯƠNG PHÁP .34 CHƯƠNG 3: VÍ DỤ MINH HOẠ 35 BÀI TOÁN TÌM ỨNG SUẤT, CHUYỂN VỊ (VÍ DỤ 1) .35 BÀI TOÁN TÌM TẦN SỐ DAO ĐỘNG RIÊNG CHO DẦM (VÍ DỤ 2) .38 BÀI TOÁN TÌM TẦN SỐ DAO ĐỘNG RIÊNG CHO TẤM CÓ SƯỜN (VÍ DỤ 3) 39 BÀI TOÁN TÌM TẢI TRỌNG TỚI HẠN (VÍ DỤ 4) 44 TÌM CHUYỂN VỊ TẤM CÓ SƯỜN (NHƯ VÍ DỤ 3) CHỊU TÁC DỤNG TẢI TRỌNG ĐIỀU H OÀ – SO SÁNH VỚI BÀI TOÁN TĨNH CÓ CÙNG CƯỜNG ĐỘ (VÍ DỤ 5) 45 TÌM CHUYỂN VỊ TẤM CÓ SƯỜN CHỊU TÁC DỤNG TẢI TRỌNG BẤT KỲ(VÍ DỤ 6) 49 CHƯƠNG 4: KHẢO SÁT ẢNH HƯỞNG CỦA SƯỜN GIA CØNG ĐẾN DAO ĐỘNG CỦA TẤM CHỬ NHẬT 51 ĐẶT VẤN ĐỀ .51 KHẢO SÁT TẤM KHÔNG SƯỜN .52 KHẢO SÁT ẢNH HƯỞNG CỦA CHIỀU CAO SƯỜN ĐẾN TẦN SỐ DAO ĐỘNG RIÊNG VÀ MODE SHAPE TƯƠNG ỨNG 53 KHẢO SÁT ẢNH HƯỞNG CỦA KHOẢNG CÁCH CÁC SƯỜN 58 KHẢO SÁT ẢNH HƯỞNG TỶ LỆ GIỮA CHIỀU DÀI & CHIỀU DÀY TẤM ĐẾN TẦN SỐ RIÊNG THỨ NHẤT 62 KHẢO SÁT ẢNH HƯỞNG CỦA CHIỀU DÀY SƯỜN ĐỐI VỚI TẦN SỐ RIÊNG THỨ NHẤT 64 CHƯƠNG 5: GIỚI THIỆU CHƯƠNG TRÌNH 68 THUẬT TOÁN 68 THỰC HIỆN : LÊ HIỀN ANH Trang Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG SỬ DỤNG CHƯƠNG TRÌNH .70 CHƯƠNG 6: KẾT LUẬN VÀ KIẾN NGHỊ .76 KẾT LUẬN 77 KIẾN NGHỊ 78 TÀI LIỆU THAM KHẢO 80 THỰC HIỆN : LÊ HIỀN ANH Trang CHƯƠNG 1: TỔNG QUAN GIỚI THIỆU PHƯƠNG PHÁP DẢI HỮU HẠN (FINITE STRIP METHOD – FSM) FSM trường hợp đặc biệt phương pháp phần tử hữu hạn (finite element method – FEM) sử dụng phần tử chiều, điều kiện tương thích Nếu FEM người ta xấp xỉ hàm chuyển vị đa thức nội suy theo tất phương, FSM hàm chuyển vị xấp xỉ việc kết hợp chuổi hàm đa thức nội suy Do xuất phát từ việc giải toán chử nhật nên phương pháp có tên FSM (strip có nghóa dải) Thực FSM bao gồm phương pháp: lớp hữu hạn lăng trụ hữu hạn Đặc điểm FSM kết cấu phải có đặc điểm hình học không đổi theo phương mô tả chuổi hàm Cho đến FSM thành công việc phân tích toán sau:  Bài toán tấm; vỏ; kết cấu chiều có mặt cắt không đổi theo phương – sử dụng phần tử tương tự phần tử shell FEM với điều kiện biên khác  Bài toán tổ hợp kết cấu (hệ dầm sàn …)  Bài toán đàn hồi; phi đàn hồi  Tương tự sử dụng phần tử chiều (lăng trụ hữu hạn)  Bài toán nhiều lớp (lớp hữu hạn)  Cũng toán giải với phần tử bậc cao FEM  p dụng hệ toạ độ decac toạ độ trụ chí hệ toạ độ không vuông góc để giải thẳng, cong, xiên …  Sử dụng phối hợp nhiều loại hệ trục toạ độ để giải toán dải không đều; lăng trụ không …  Tấm cong ghềnh có hình chiếu chử nhật; hình bình hành hay hình quạt  Cho phép giải toán động lực học toán ổn định  Cho phép giải toán phi tuyến vật liệu; phi tuyến hình học  Cho phép kết hợp FEM FSM Thông qua phạm vi ứng dụng nêu trên, ta lầm tưởng FSM thực tất toán mà FEM giải Thực sử dụng FSM ta phải nhớ kết cấu phải có hình dạng hình học đơn giản để việc mô không cần dùng nhiều số hạng chuổi Chính FSM thường hay dùng tính toán cầu Trang Trong FSM phần chuổi hàm phần đáng để quan tâm nghiên cứu nhất, phần khác biệt FEM FSM Trong loại chuổi hàm nghiên cứu, chuổi lượng giác chuổi Spline loại chuổi phổ biến FSM sử dụng chuổi lượng giác, gọi phương pháp bán giải tích (semi analysis), phương pháp nêu bật tất đặc điểm FSM Tuy nhiên phạm vi ứng dụng không rộng nên thường sử dụng để nghiên cứu FSM với chuổi Spline thực cách mạng FSM giúp nhà khoa học mỡ rộng phạm vi ứng dụng phương pháp QUÁ TRÌNH PHÁT TRIỂN CỦA FSM Phương pháp dải hữu hạn (FSM) Y.K Cheung khởi xướng Tác phẩm xuất năm 1968 đề cập đến lời giải cho toán tựa đơn Năm 1969 Rowel Orden đưa lời giải cho toán cầu chữ nhật cách độc lập với Y.K.Cheung Sau phương pháp nghiên cứu nhiều nước tập trung vào cuối thập niên 1960 đầu thập niên 1970 Riêng việc nghiên cứu FSM để giải toán thiết kế cầu có nhiều thành tựu: • Tấm chữ nhật với điều kiện biên tổng quát Y.K Cheung năm 1968 • Cầu dầm hộp Y.K.Cheung 1969 • Lớp hữu hạn để giải toán dày Y.K.Cheung Chakrabarti năm 1971 • Lăng trụ hữu hạn để giải cầu hộp dày Zienkiewicz Too năm 1972 • Tấm cong tròn cầu dầm hộp (Y.K Cheung năm 1969 ; Y.K.Cheung M.S.Cheung năm 1971) • Tấm gấp khúc cầu dầm hộp ( Brown Ghali, 1972 1975) • Cầu liên tục, tiếp cận mềm (tương tự phương pháp lực) M.S.Cheung, Y.K.Cheung Ghali năm 1970 • Cầu dầm hộp liên tục có sườn ngang, tiếp cận mềm Loo năm 1975 • Cầu cầu dầm hộp tiếp cận cứng Wu Cheung năm 1974; Decourt Cheung năm 1978) • Tính dao động tự (M.S.Cheung Cheung năm 1971) • Phản ứng động cầu tác động tải trọng di động Smith năm 1973 • Tính ổn định Premieniecki năm 1973 Wittrick Plank năm 1974 Trang Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG function P=Nngampt(px,py,pz,b,L,s,xc,xd,zc,zd); syms x y real a=sin(s)-sinh(s)/(cos(s)-cosh(s)); if xc==xd if zc>zd px1=pz; pz1=-px; else px1=-pz; pz1=px; end else beta=atan((zc-zd)/(xc-xd)); px1=px*cos(beta)+pz*sin(beta); pz1=pz*cos(beta)-px*sin(beta); end x1=x/b; P(1)=int(int((1-x1)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*px1,x,0,b),y,0,L); P(2)=int(int((1-x1)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*py,x,0,b),y,0,L/2) -int(int((1-x1)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*py,x,0,b),y,L/2,L); P(3)=int(int((1-3*x1^2+2*x1^3)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*pz1,x,0,b),y,0,L); P(4)=int(int(x*(1-2*x1+x1^2)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*pz1,x,0,b),y,0,L); P(5)=int(int(x1*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*px1,x,0,b),y,0,L); P(6)=int(int(x1*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*py,x,0,b),y,0,L/2) -int(int(x1*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*py,x,0,b),y,L/2,L); P(7)=int(int((3*x1^2-2*x1^3)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*pz1,x,0,b),y,0,L); P(8)=int(int(x*(x1^2-x1)*(sin(s*y/L)-sinh(s*y/L)-a*(cos(s*y/L)cosh(s*y/L)))*pz1,x,0,b),y,0,L); 26 file Nptcong.m function P=Nptcong(px,py,pz,b,L,s,xc,xd,zc,zd); syms x y real if xc==xd if zc>zd px1=pz; pz1=-px; else px1=-pz; pz1=px; end THỰC HIỆN : LÊ HIỀN ANH Trang 45 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG else beta=atan((zc-zd)/(xc-xd)); px1=px*cos(beta)+pz*sin(beta); pz1=pz*cos(beta)-px*sin(beta); end x1=x/b; P(1)=int(int((1-x1)*sin(s*pi*y/L)*xd*px1,x,0,b),y,0,L); P(2)=int(int((1-x1)*cos(s*pi*y/L)*xd*py,x,0,b),y,0,L/2) -int(int((1-x1)*cos(s*pi*y/L)*xd*py,x,0,b),y,L/2,L); P(3)=int(int((1-3*x1^2+2*x1^3)*sin(s*pi*y/L)*xd*pz1,x,0,b),y,0,L); P(4)=int(int(x*(1-2*x1+x1^2)*sin(s*pi*y/L)*xd*pz1,x,0,b),y,0,L); P(5)=int(int(x1*sin(s*pi*y/L)*xd*px1,x,0,b),y,0,L); P(6)=int(int(x1*cos(s*pi*y/L)*xd*py,x,0,b),y,0,L/2) -int(int(x1*cos(s*pi*y/L)*xd*py,x,0,b),y,L/2,L); P(7)=int(int((3*x1^2-2*x1^3)*sin(s*pi*y/L)*xd*pz1,x,0,b),y,0,L); P(8)=int(int(x*(x1^2-x1)*sin(s*pi*y/L)*xd*pz1,x,0,b),y,0,L); 27 file Mtong.m function M=Mtong(ro,t,L,m,n,X,Z,nutdau,nutcuoi,s); M=zeros(4*m,4*m); for i=1:n % n la so phan tu b=sqrt((X(nutcuoi(i))-X(nutdau(i)))^2+(Z(nutcuoi(i))-Z(nutdau(i)))^2); T=Doitruc(X(nutcuoi(i)),X(nutdau(i)),Z(nutcuoi(i)),Z(nutdau(i))); % Gan bien chuyen vi u1(i)=4*nutdau(i)-3; v1(i)=4*nutdau(i)-2; w1(i)=4*nutdau(i)-1; r1(i)=4*nutdau(i); u2(i)=4*nutcuoi(i)-3; v2(i)=4*nutcuoi(i)-2; w2(i)=4*nutcuoi(i)-1; r2(i)=4*nutcuoi(i); % Tinh K tong the cua phan tu M1=Mpt(t(i),b,L,ro); mx=T*M1*T'; M(u1(i),u1(i))=M(u1(i),u1(i))+mx(1,1); M(u1(i),v1(i))=M(u1(i),v1(i))+mx(1,2); M(u1(i),w1(i))=M(u1(i),w1(i))+mx(1,3); M(u1(i),r1(i))=M(u1(i),r1(i))+mx(1,4); M(u1(i),u2(i))=M(u1(i),u2(i))+mx(1,5); M(u1(i),v2(i))=M(u1(i),v2(i))+mx(1,6); M(u1(i),w2(i))=M(u1(i),w2(i))+mx(1,7); M(u1(i),r2(i))=M(u1(i),r2(i))+mx(1,8); M(v1(i),u1(i))=M(v1(i),u1(i))+mx(2,1); M(v1(i),v1(i))=M(v1(i),v1(i))+mx(2,2); M(v1(i),w1(i))=M(v1(i),w1(i))+mx(2,3); M(v1(i),r1(i))=M(v1(i),r1(i))+mx(2,4); THỰC HIỆN : LÊ HIỀN ANH Trang 46 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG M(v1(i),u2(i))=M(v1(i),u2(i))+mx(2,5); M(v1(i),v2(i))=M(v1(i),v2(i))+mx(2,6); M(v1(i),w2(i))=M(v1(i),w2(i))+mx(2,7); M(v1(i),r2(i))=M(v1(i),r2(i))+mx(2,8); M(w1(i),u1(i))=M(w1(i),u1(i))+mx(3,1); M(w1(i),v1(i))=M(w1(i),v1(i))+mx(3,2); M(w1(i),w1(i))=M(w1(i),w1(i))+mx(3,3); M(w1(i),r1(i))=M(w1(i),r1(i))+mx(3,4); M(w1(i),u2(i))=M(w1(i),u2(i))+mx(3,5); M(w1(i),v2(i))=M(w1(i),v2(i))+mx(3,6); M(w1(i),w2(i))=M(w1(i),w2(i))+mx(3,7); M(w1(i),r2(i))=M(w1(i),r2(i))+mx(3,8); M(r1(i),u1(i))=M(r1(i),u1(i))+mx(4,1); M(r1(i),v1(i))=M(r1(i),v1(i))+mx(4,2); M(r1(i),w1(i))=M(r1(i),w1(i))+mx(4,3); M(r1(i),r1(i))=M(r1(i),r1(i))+mx(4,4); M(r1(i),u2(i))=M(r1(i),u2(i))+mx(4,5); M(r1(i),v2(i))=M(r1(i),v2(i))+mx(4,6); M(r1(i),w2(i))=M(r1(i),w2(i))+mx(4,7); M(r1(i),r2(i))=M(r1(i),r2(i))+mx(4,8); M(u2(i),u1(i))=M(u2(i),u1(i))+mx(5,1); M(u2(i),v1(i))=M(u2(i),v1(i))+mx(5,2); M(u2(i),w1(i))=M(u2(i),w1(i))+mx(5,3); M(u2(i),r1(i))=M(u2(i),r1(i))+mx(5,4); M(u2(i),u2(i))=M(u2(i),u2(i))+mx(5,5); M(u2(i),v2(i))=M(u2(i),v2(i))+mx(5,6); M(u2(i),w2(i))=M(u2(i),w2(i))+mx(5,7); M(u2(i),r2(i))=M(u2(i),r2(i))+mx(5,8); M(v2(i),u1(i))=M(v2(i),u1(i))+mx(6,1); M(v2(i),v1(i))=M(v2(i),v1(i))+mx(6,2); M(v2(i),w1(i))=M(v2(i),w1(i))+mx(6,3); M(v2(i),r1(i))=M(v2(i),r1(i))+mx(6,4); M(v2(i),u2(i))=M(v2(i),u2(i))+mx(6,5); M(v2(i),v2(i))=M(v2(i),v2(i))+mx(6,6); M(v2(i),w2(i))=M(v2(i),w2(i))+mx(6,7); M(v2(i),r2(i))=M(v2(i),r2(i))+mx(6,8); M(w2(i),u1(i))=M(w2(i),u1(i))+mx(7,1); M(w2(i),v1(i))=M(w2(i),v1(i))+mx(7,2); M(w2(i),w1(i))=M(w2(i),w1(i))+mx(7,3); M(w2(i),r1(i))=M(w2(i),r1(i))+mx(7,4); M(w2(i),u2(i))=M(w2(i),u2(i))+mx(7,5); M(w2(i),v2(i))=M(w2(i),v2(i))+mx(7,6); M(w2(i),w2(i))=M(w2(i),w2(i))+mx(7,7); M(w2(i),r2(i))=M(w2(i),r2(i))+mx(7,8); M(r2(i),u1(i))=M(r2(i),u1(i))+mx(8,1); M(r2(i),v1(i))=M(r2(i),v1(i))+mx(8,2); M(r2(i),w1(i))=M(r2(i),w1(i))+mx(8,3); THỰC HIỆN : LÊ HIỀN ANH Trang 47 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THAÉNG M(r2(i),r1(i))=M(r2(i),r1(i))+mx(8,4); M(r2(i),u2(i))=M(r2(i),u2(i))+mx(8,5); M(r2(i),v2(i))=M(r2(i),v2(i))+mx(8,6); M(r2(i),w2(i))=M(r2(i),w2(i))+mx(8,7); M(r2(i),r2(i))=M(r2(i),r2(i))+mx(8,8); end %end for n 28 file Mtongcong.m function M=Mtongcong(ro,t,L,m,n,X,Z,nutdau,nutcuoi,s); M=zeros(4*m,4*m); for i=1:n % n la so phan tu b=sqrt((X(nutcuoi(i))-X(nutdau(i)))^2+(Z(nutcuoi(i))-Z(nutdau(i)))^2); T=Doitruccong(X(nutcuoi(i)),X(nutdau(i)),Z(nutcuoi(i)),Z(nutdau(i))); % Gan bien chuyen vi u1(i)=4*nutdau(i)-3; v1(i)=4*nutdau(i)-2; w1(i)=4*nutdau(i)-1; r1(i)=4*nutdau(i); u2(i)=4*nutcuoi(i)-3; v2(i)=4*nutcuoi(i)-2; w2(i)=4*nutcuoi(i)-1; r2(i)=4*nutcuoi(i); % Tinh K tong the cua phan tu M1=Mptcong(t(i),b,L,ro,X(nutdau(i))); mx=T*M1*T'; M(u1(i),u1(i))=M(u1(i),u1(i))+mx(1,1); M(u1(i),v1(i))=M(u1(i),v1(i))+mx(1,2); M(u1(i),w1(i))=M(u1(i),w1(i))+mx(1,3); M(u1(i),r1(i))=M(u1(i),r1(i))+mx(1,4); M(u1(i),u2(i))=M(u1(i),u2(i))+mx(1,5); M(u1(i),v2(i))=M(u1(i),v2(i))+mx(1,6); M(u1(i),w2(i))=M(u1(i),w2(i))+mx(1,7); M(u1(i),r2(i))=M(u1(i),r2(i))+mx(1,8); M(v1(i),u1(i))=M(v1(i),u1(i))+mx(2,1); M(v1(i),v1(i))=M(v1(i),v1(i))+mx(2,2); M(v1(i),w1(i))=M(v1(i),w1(i))+mx(2,3); M(v1(i),r1(i))=M(v1(i),r1(i))+mx(2,4); M(v1(i),u2(i))=M(v1(i),u2(i))+mx(2,5); M(v1(i),v2(i))=M(v1(i),v2(i))+mx(2,6); M(v1(i),w2(i))=M(v1(i),w2(i))+mx(2,7); M(v1(i),r2(i))=M(v1(i),r2(i))+mx(2,8); M(w1(i),u1(i))=M(w1(i),u1(i))+mx(3,1); M(w1(i),v1(i))=M(w1(i),v1(i))+mx(3,2); M(w1(i),w1(i))=M(w1(i),w1(i))+mx(3,3); M(w1(i),r1(i))=M(w1(i),r1(i))+mx(3,4); M(w1(i),u2(i))=M(w1(i),u2(i))+mx(3,5); M(w1(i),v2(i))=M(w1(i),v2(i))+mx(3,6); M(w1(i),w2(i))=M(w1(i),w2(i))+mx(3,7); THỰC HIỆN : LÊ HIỀN ANH Trang 48 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG M(w1(i),r2(i))=M(w1(i),r2(i))+mx(3,8); M(r1(i),u1(i))=M(r1(i),u1(i))+mx(4,1); M(r1(i),v1(i))=M(r1(i),v1(i))+mx(4,2); M(r1(i),w1(i))=M(r1(i),w1(i))+mx(4,3); M(r1(i),r1(i))=M(r1(i),r1(i))+mx(4,4); M(r1(i),u2(i))=M(r1(i),u2(i))+mx(4,5); M(r1(i),v2(i))=M(r1(i),v2(i))+mx(4,6); M(r1(i),w2(i))=M(r1(i),w2(i))+mx(4,7); M(r1(i),r2(i))=M(r1(i),r2(i))+mx(4,8); M(u2(i),u1(i))=M(u2(i),u1(i))+mx(5,1); M(u2(i),v1(i))=M(u2(i),v1(i))+mx(5,2); M(u2(i),w1(i))=M(u2(i),w1(i))+mx(5,3); M(u2(i),r1(i))=M(u2(i),r1(i))+mx(5,4); M(u2(i),u2(i))=M(u2(i),u2(i))+mx(5,5); M(u2(i),v2(i))=M(u2(i),v2(i))+mx(5,6); M(u2(i),w2(i))=M(u2(i),w2(i))+mx(5,7); M(u2(i),r2(i))=M(u2(i),r2(i))+mx(5,8); M(v2(i),u1(i))=M(v2(i),u1(i))+mx(6,1); M(v2(i),v1(i))=M(v2(i),v1(i))+mx(6,2); M(v2(i),w1(i))=M(v2(i),w1(i))+mx(6,3); M(v2(i),r1(i))=M(v2(i),r1(i))+mx(6,4); M(v2(i),u2(i))=M(v2(i),u2(i))+mx(6,5); M(v2(i),v2(i))=M(v2(i),v2(i))+mx(6,6); M(v2(i),w2(i))=M(v2(i),w2(i))+mx(6,7); M(v2(i),r2(i))=M(v2(i),r2(i))+mx(6,8); M(w2(i),u1(i))=M(w2(i),u1(i))+mx(7,1); M(w2(i),v1(i))=M(w2(i),v1(i))+mx(7,2); M(w2(i),w1(i))=M(w2(i),w1(i))+mx(7,3); M(w2(i),r1(i))=M(w2(i),r1(i))+mx(7,4); M(w2(i),u2(i))=M(w2(i),u2(i))+mx(7,5); M(w2(i),v2(i))=M(w2(i),v2(i))+mx(7,6); M(w2(i),w2(i))=M(w2(i),w2(i))+mx(7,7); M(w2(i),r2(i))=M(w2(i),r2(i))+mx(7,8); M(r2(i),u1(i))=M(r2(i),u1(i))+mx(8,1); M(r2(i),v1(i))=M(r2(i),v1(i))+mx(8,2); M(r2(i),w1(i))=M(r2(i),w1(i))+mx(8,3); M(r2(i),r1(i))=M(r2(i),r1(i))+mx(8,4); M(r2(i),u2(i))=M(r2(i),u2(i))+mx(8,5); M(r2(i),v2(i))=M(r2(i),v2(i))+mx(8,6); M(r2(i),w2(i))=M(r2(i),w2(i))+mx(8,7); M(r2(i),r2(i))=M(r2(i),r2(i))+mx(8,8); end %end for n 29 file Mpt.m function M=Mpt(t,b,L,ro);% khong phu thuoc s M=zeros(8,8); M(1,1)=1/6*ro*t*b*L; THỰC HIỆN : LÊ HIỀN ANH Trang 49 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THAÉNG M(1,5)=1/12*ro*t*b*L; M(2,2)=1/6*ro*t*b*L; M(2,6)=1/12*ro*t*b*L; M(5,1)=M(1,5); M(5,5)=1/6*ro*t*b*L; M(6,2)=M(2,6); M(6,6)=1/6*ro*t*b*L; M(3,3)=13/70*ro*t*b*L; M(3,4)=11/420*ro*t*b^2*L; M(3,7)=9/140*ro*t*b*L; M(3,8)=-13/840*ro*t*b^2*L; M(4,3)=M(3,4); M(4,4)=1/210*ro*t*b^3*L; M(4,7)=13/840*ro*t*b^2*L; M(4,8)=-1/280*ro*t*b^3*L; M(7,3)=M(3,7); M(7,4)=M(4,7); M(7,7)=13/70*ro*t*b*L; M(7,8)=-11/420*ro*t*b^2*L; M(8,3)=M(3,8); M(8,4)=M(4,8); M(8,7)=M(7,8); M(8,8)=1/210*ro*t*b^3*L; 30 file Mptcong.m function M=Mptcong(t,b,L,ro,r);% khong phu thuoc s p=pi; M=zeros(8,8); M(1,1)=1/6*ro*t*r*b*L; M(1,2)=0; M(1,3)=7/40*ro*t*r*b*L; M(1,4)=1/40*ro*t*r*b^2*L; M(1,5)=1/12*ro*t*r*b*L; M(1,6)=0; M(1,7)=3/40*ro*t*r*b*L; M(1,8)=0; M(2,1)=0; M(2,2)=1/6*ro*t*r*b*L; M(2,3)=0; M(2,4)=0; M(2,5)=0; M(2,6)=1/12*ro*t*r*b*L; M(2,7)=0; M(2,8)=3/40*ro*t*r*b*L; THỰC HIỆN : LÊ HIỀN ANH Trang 50 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG M(3,1)=7/40*ro*t*r*b*L; M(3,2)=0; M(3,3)=13/70*ro*t*r*b*L; M(3,4)=11/420*ro*t*r*b^2*L; M(3,5)=3/40*ro*t*r*b*L; M(3,6)=0; M(3,7)=9/140*ro*t*r*b*L; M(3,8)=0; M(4,1)=1/40*ro*t*r*b^2*L; M(4,2)=0; M(4,3)=11/420*ro*t*r*b^2*L; M(4,4)=1/210*ro*t*r*b^3*L; M(4,5)=1/60*ro*t*r*b^2*L; M(4,6)=0; M(4,7)=13/840*ro*t*r*b^2*L; M(4,8)=0; M(5,1)=1/12*ro*t*r*b*L; M(5,2)=0; M(5,3)=3/40*ro*t*r*b*L; M(5,4)=1/60*ro*t*r*b^2*L; M(5,5)=1/6*ro*t*r*b*L; M(5,6)=0; M(5,7)=7/40*ro*t*r*b*L; M(5,8)=0; M(6,1)=0; M(6,2)=1/12*ro*t*r*b*L; M(6,3)=0; M(6,4)=0; M(6,5)=0; M(6,6)=1/6*ro*t*r*b*L; M(6,7)=0; M(6,8)=7/40*ro*t*r*b*L; M(7,1)=3/40*ro*t*r*b*L; M(7,2)=0; M(7,3)=9/140*ro*t*r*b*L; M(7,4)=13/840*ro*t*r*b^2*L; M(7,5)=7/40*ro*t*r*b*L; M(7,6)=0; M(7,7)=13/70*ro*t*r*b*L; M(7,8)=0; M(8,1)=0; M(8,2)=3/40*ro*t*r*b*L; M(8,3)=0; M(8,4)=0; M(8,5)=0; M(8,6)=7/40*ro*t*r*b*L; M(8,7)=0; THỰC HIỆN : LÊ HIỀN ANH Trang 51 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG M(8,8)=13/70*ro*t*r*b*L; 31 file KGtong.m function KQ=KGtong(E,v,L,t,m,n,s,X,Z,nutdau,nutcuoi,cvi); KQ=zeros(4*m,4*m); for i=1:n % n la so phan tu b=sqrt((X(nutcuoi(i))-X(nutdau(i)))^2+(Z(nutcuoi(i))-Z(nutdau(i)))^2); %T=Doitruc(X(nutcuoi(i)),X(nutdau(i)),Z(nutcuoi(i)),Z(nutdau(i))); % Gan bien chuyen vi u1(i)=4*nutdau(i)-3; v1(i)=4*nutdau(i)-2; w1(i)=4*nutdau(i)-1; r1(i)=4*nutdau(i); u2(i)=4*nutcuoi(i)-3; v2(i)=4*nutcuoi(i)-2; w2(i)=4*nutcuoi(i)-1; r2(i)=4*nutcuoi(i); % Tinh K tong the cua phan tu K=KGpt1(E,v,L,b,t(i),s,cvi(u1(i)),cvi(v1(i)),cvi(w1(i)),cvi(r1(i)),cvi(u2(i )),cvi(v2(i)),cvi(w2(i)),cvi(r2(i))); %K=T'*K1*T; KQ(u1(i),u1(i))=KQ(u1(i),u1(i))+K(1,1); KQ(u1(i),v1(i))=KQ(u1(i),v1(i))+K(1,2); KQ(u1(i),w1(i))=KQ(u1(i),w1(i))+K(1,3); KQ(u1(i),r1(i))=KQ(u1(i),r1(i))+K(1,4); KQ(u1(i),u2(i))=KQ(u1(i),u2(i))+K(1,5); KQ(u1(i),v2(i))=KQ(u1(i),v2(i))+K(1,6); KQ(u1(i),w2(i))=KQ(u1(i),w2(i))+K(1,7); KQ(u1(i),r2(i))=KQ(u1(i),r2(i))+K(1,8); KQ(v1(i),u1(i))=KQ(v1(i),u1(i))+K(2,1); KQ(v1(i),v1(i))=KQ(v1(i),v1(i))+K(2,2); KQ(v1(i),w1(i))=KQ(v1(i),w1(i))+K(2,3); KQ(v1(i),r1(i))=KQ(v1(i),r1(i))+K(2,4); KQ(v1(i),u2(i))=KQ(v1(i),u2(i))+K(2,5); KQ(v1(i),v2(i))=KQ(v1(i),v2(i))+K(2,6); KQ(v1(i),w2(i))=KQ(v1(i),w2(i))+K(2,7); KQ(v1(i),r2(i))=KQ(v1(i),r2(i))+K(2,8); KQ(w1(i),u1(i))=KQ(w1(i),u1(i))+K(3,1); KQ(w1(i),v1(i))=KQ(w1(i),v1(i))+K(3,2); KQ(w1(i),w1(i))=KQ(w1(i),w1(i))+K(3,3); KQ(w1(i),r1(i))=KQ(w1(i),r1(i))+K(3,4); KQ(w1(i),u2(i))=KQ(w1(i),u2(i))+K(3,5); KQ(w1(i),v2(i))=KQ(w1(i),v2(i))+K(3,6); KQ(w1(i),w2(i))=KQ(w1(i),w2(i))+K(3,7); KQ(w1(i),r2(i))=KQ(w1(i),r2(i))+K(3,8); KQ(r1(i),u1(i))=KQ(r1(i),u1(i))+K(4,1); KQ(r1(i),v1(i))=KQ(r1(i),v1(i))+K(4,2); THỰC HIỆN : LÊ HIỀN ANH Trang 52 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG KQ(r1(i),w1(i))=KQ(r1(i),w1(i))+K(4,3); KQ(r1(i),r1(i))=KQ(r1(i),r1(i))+K(4,4); KQ(r1(i),u2(i))=KQ(r1(i),u2(i))+K(4,5); KQ(r1(i),v2(i))=KQ(r1(i),v2(i))+K(4,6); KQ(r1(i),w2(i))=KQ(r1(i),w2(i))+K(4,7); KQ(r1(i),r2(i))=KQ(r1(i),r2(i))+K(4,8); KQ(u2(i),u1(i))=KQ(u2(i),u1(i))+K(5,1); KQ(u2(i),v1(i))=KQ(u2(i),v1(i))+K(5,2); KQ(u2(i),w1(i))=KQ(u2(i),w1(i))+K(5,3); KQ(u2(i),r1(i))=KQ(u2(i),r1(i))+K(5,4); KQ(u2(i),u2(i))=KQ(u2(i),u2(i))+K(5,5); KQ(u2(i),v2(i))=KQ(u2(i),v2(i))+K(5,6); KQ(u2(i),w2(i))=KQ(u2(i),w2(i))+K(5,7); KQ(u2(i),r2(i))=KQ(u2(i),r2(i))+K(5,8); KQ(v2(i),u1(i))=KQ(v2(i),u1(i))+K(6,1); KQ(v2(i),v1(i))=KQ(v2(i),v1(i))+K(6,2); KQ(v2(i),w1(i))=KQ(v2(i),w1(i))+K(6,3); KQ(v2(i),r1(i))=KQ(v2(i),r1(i))+K(6,4); KQ(v2(i),u2(i))=KQ(v2(i),u2(i))+K(6,5); KQ(v2(i),v2(i))=KQ(v2(i),v2(i))+K(6,6); KQ(v2(i),w2(i))=KQ(v2(i),w2(i))+K(6,7); KQ(v2(i),r2(i))=KQ(v2(i),r2(i))+K(6,8); KQ(w2(i),u1(i))=KQ(w2(i),u1(i))+K(7,1); KQ(w2(i),v1(i))=KQ(w2(i),v1(i))+K(7,2); KQ(w2(i),w1(i))=KQ(w2(i),w1(i))+K(7,3); KQ(w2(i),r1(i))=KQ(w2(i),r1(i))+K(7,4); KQ(w2(i),u2(i))=KQ(w2(i),u2(i))+K(7,5); KQ(w2(i),v2(i))=KQ(w2(i),v2(i))+K(7,6); KQ(w2(i),w2(i))=KQ(w2(i),w2(i))+K(7,7); KQ(w2(i),r2(i))=KQ(w2(i),r2(i))+K(7,8); KQ(r2(i),u1(i))=KQ(r2(i),u1(i))+K(8,1); KQ(r2(i),v1(i))=KQ(r2(i),v1(i))+K(8,2); KQ(r2(i),w1(i))=KQ(r2(i),w1(i))+K(8,3); KQ(r2(i),r1(i))=KQ(r2(i),r1(i))+K(8,4); KQ(r2(i),u2(i))=KQ(r2(i),u2(i))+K(8,5); KQ(r2(i),v2(i))=KQ(r2(i),v2(i))+K(8,6); KQ(r2(i),w2(i))=KQ(r2(i),w2(i))+K(8,7); KQ(r2(i),r2(i))=KQ(r2(i),r2(i))+K(8,8); end %end for n 32 file KGpt.m (khai trieån KG theo M.S.Cheung) function K=KGpt(E,v,L,b,t,s,u1,v1,w1,r1,u2,v2,w2,r2); K=zeros(8,8); K(3,3)=-1/210*t*E*(-336*L^3*u1-12*pi^3*b^3*s^3*v240*pi^3*b^3*s^3*v1+336*L^3*u2-8*pi^2*b^2*s^2*v*u1*L-8*pi^2*b^2*s^2*u1*L26*pi^2*b^2*s^2*u2*L-98*v*v2*L^2*b*s*pi238*v*v1*b*s*pi*L^2+70*v1*b*s*pi*L^270*v2*b*s*pi*L^2+78*pi^2*b^2*s^2*v*u2*L)/b^2/(-1+v^2)/L^2/s/pi; THỰC HIỆN : LÊ HIỀN ANH Trang 53 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG K(3,4)=-1/1260*t*E*(-168*v*v2*L^2*b*s*pi+66*pi^2*b^2*s^2*v*u2*L66*pi^2*b^2*s^2*v*u1*L+22*pi^2*b^2*s^2*u1*L-168*L^3*u122*pi^2*b^2*s^2*u2*L+168*L^3*u2-30*pi^3*b^3*s^3*v1-14*pi^3*b^3*s^3*v2)/b/(1+v^2)/L^2/s/pi; K(3,7)=1/140*t*E*(-112*v*v2*L^2*b*s*pi+18*pi^2*b^2*s^2*v*u1*L18*pi^2*b^2*s^2*v*u2*L-112*v*v1*b*s*pi*L^2224*L^3*u1+224*L^3*u2+6*pi^2*b^2*s^2*u2*L+6*pi^3*b^3*s^3*v1+6*pi^3*b^3*s^3* v2-6*pi^2*b^2*s^2*u1*L)/b^2/(-1+v^2)/L^2/s/pi; K(3,8)=-1/2520*t*E*(-78*pi^2*b^2*s^2*v*u2*L+78*pi^2*b^2*s^2*v*u1*L336*v*v1*b*s*pi*L^2+26*pi^2*b^2*s^2*u2*L+336*L^3*u2-26*pi^2*b^2*s^2*u1*L336*L^3*u1+24*pi^3*b^3*s^3*v2+28*pi^3*b^3*s^3*v1)/b/(-1+v^2)/L^2/s/pi; K(4,3)=K(3,4); K(4,4)=-1/2520*t*E*(-24*pi^2*b^2*s^2*v*u1*L112*v*v2*L^2*b*s*pi+24*pi^2*b^2*s^2*v*u2*L-336*v*v1*b*s*pi*L^2448*L^3*u1+448*L^3*u2-10*pi^3*b^3*s^3*v16*pi^3*b^3*s^3*v2+8*pi^2*b^2*s^2*u1*L-8*pi^2*b^2*s^2*u2*L)/(1+v^2)/L^2/s/pi; K(4,7)=1/2520*t*E*(-336*v*v2*L^2*b*s*pi78*pi^2*b^2*s^2*v*u2*L+78*pi^2*b^2*s^2*v*u1*L+26*pi^2*b^2*s^2*u2*L336*L^3*u126*pi^2*b^2*s^2*u1*L+336*L^3*u2+28*pi^3*b^3*s^3*v2+24*pi^3*b^3*s^3*v1)/b/(1+v^2)/L^2/s/pi; K(4,8)=1/2520*t*E*(-18*pi^2*b^2*s^2*v*u1*L56*v*v1*b*s*pi*L^2+18*pi^2*b^2*s^2*v*u2*L-56*v*v2*L^2*b*s*pi6*pi^3*b^3*s^3*v1-112*L^3*u1+112*L^3*u2-6*pi^3*b^3*s^3*v26*pi^2*b^2*s^2*u2*L+6*pi^2*b^2*s^2*u1*L)/(-1+v^2)/L^2/s/pi; K(7,3)=K(3,7); K(7,4)=K(4,7); K(7,7)=-1/210*t*E*(-336*L^3*u1-40*pi^3*b^3*s^3*v212*pi^3*b^3*s^3*v1+336*L^3*u278*pi^2*b^2*s^2*v*u1*L+26*pi^2*b^2*s^2*u1*L+14*pi^3*b^2*s^2*u2*L238*v*v2*L^2*b*s*pi-98*v*v1*b*s*pi*L^270*v1*b*s*pi*L^2+70*v2*b*s*pi*L^2+8*pi^2*b^2*s^2*v*u2*L)/b^2/(1+v^2)/L^2/s/pi; K(7,8)=1/1260*t*E*(66*pi^2*b^2*s^2*v*u2*L-66*pi^2*b^2*s^2*v*u1*L168*v*v1*b*s*pi*L^2-22*pi^2*b^2*s^2*u2*L-168*L^3*u1+168*L^3*u214*pi^3*b^3*s^3*v1-30*pi^3*b^3*s^3*v2+22*pi^2*b^2*s^2*u1*L)/b/(1+v^2)/L^2/s/pi; K(8,3)=K(3,8); K(8,4)=K(4,8); K(8,7)=K(7,8); K(8,8)=-1/2520*t*E*(24*pi^2*b^2*s^2*v*u2*L-336*v*v2*L^2*b*s*pi24*pi^2*b^2*s^2*v*u1*L-112*v*v1*b*s*pi*L^2-448*L^3*u1+448*L^3*u28*pi^2*b^2*s^2*u2*L+8*pi^2*b^2*s^2*u1*L-10*pi^3*b^3*s^3*v26*pi^3*b^3*s^3*v1)/(-1+v^2)/L^2/s/pi; 33 file KGpt1.m (khai triển KG cách đạo hàm trực tiếp) function K=KGpt1(E,v,L,b,t,s,u1,v1,w1,r1,u2,v2,w2,r2); K=zeros(8,8); p=pi; THỰC HIỆN : LÊ HIỀN ANH Trang 54 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG K(1,1)=0; K(1,2)=0; K(1,3)=-1/2520*E*(-4032*(-1)^s*L^2*w2+4032*L^2*w2-4032*L^2*w1+270*b^2*v*((1)^s)^3*s^2*p^2*w2+110*b^3*v*((-1)^s)^3*s^2*p^2*r1+360*b^2*v*((1)^s)^3*s^2*p^2*w1-336*L^2*r1*b+4032*(-1)^s*L^2*w1-65*b^3*v*((1)^s)^3*s^2*p^2*r2+336*(-1)^s*L^2*r2*b+336*(-1)^s*L^2*r1*b264*b^2*s^2*p^2*w1+54*b^2*s^2*p^2*w2-13*b^3*s^2*p^2*r2110*b^3*v*s^2*p^2*r1-360*b^2*v*s^2*p^2*w1270*b^2*v*s^2*p^2*w2+22*b^3*s^2*p^2*r1+65*b^3*v*s^2*p^2*r2336*L^2*r2*b+264*b^2*((-1)^s)^3*s^2*p^2*w1+13*b^3*((-1)^s)^3*s^2*p^2*r222*b^3*((-1)^s)^3*s^2*p^2*r1-54*b^2*((-1)^s)^3*s^2*p^2*w2)/(1+v)/L/p/s/(v1)/b^2*t; K(1,4)=-1/2520*E*(-336*(-1)^s*L^2*w2+336*L^2*w2-336*L^2*w1+65*b^2*v*((1)^s)^3*s^2*p^2*w2+20*b^3*v*((-1)^s)^3*s^2*p^2*r1+110*b^2*v*((1)^s)^3*s^2*p^2*w1-448*L^2*r1*b+336*(-1)^s*L^2*w1-15*b^3*v*((1)^s)^3*s^2*p^2*r2-112*(-1)^s*L^2*r2*b+448*(1)^s*L^2*r1*b+22*b^2*s^2*p^2*w1+13*b^2*s^2*p^2*w2-3*b^3*s^2*p^2*r220*b^3*v*s^2*p^2*r1-110*b^2*v*s^2*p^2*w165*b^2*v*s^2*p^2*w2+4*b^3*s^2*p^2*r1+15*b^3*v*s^2*p^2*r2+112*L^2*r2*b22*b^2*((-1)^s)^3*s^2*p^2*w1+3*b^3*((-1)^s)^3*s^2*p^2*r2-4*b^3*((1)^s)^3*s^2*p^2*r1-13*b^2*((-1)^s)^3*s^2*p^2*w2)/(-1+v^2)/L/b/s/p*t; K(1,5)=0; K(1,6)=0; K(1,7)=-1/2520*E*(4032*(-1)^s*L^2*w2-4032*L^2*w2+4032*L^2*w1+780*b^2*v*((1)^s)^3*s^2*p^2*w2+65*b^3*v*((-1)^s)^3*s^2*p^2*r1+270*b^2*v*((1)^s)^3*s^2*p^2*w1+336*L^2*r1*b-4032*(-1)^s*L^2*w1-110*b^3*v*((1)^s)^3*s^2*p^2*r2-336*(-1)^s*L^2*r2*b-336*(1)^s*L^2*r1*b+54*b^2*s^2*p^2*w1+156*b^2*s^2*p^2*w2-22*b^3*s^2*p^2*r265*b^3*v*s^2*p^2*r1-270*b^2*v*s^2*p^2*w1780*b^2*v*s^2*p^2*w2+13*b^3*s^2*p^2*r1+110*b^3*v*s^2*p^2*r2+336*L^2*r2*b54*b^2*((-1)^s)^3*s^2*p^2*w1+22*b^3*((-1)^s)^3*s^2*p^2*r2-13*b^3*((1)^s)^3*s^2*p^2*r1-156*b^2*((-1)^s)^3*s^2*p^2*w2)/b^2/p/s/(-1+v^2)/L*t; K(1,8)=1/2520*E*(336*(-1)^s*L^2*w2-336*L^2*w2+336*L^2*w1+110*b^2*v*((1)^s)^3*s^2*p^2*w2+15*b^3*v*((-1)^s)^3*s^2*p^2*r1+65*b^2*v*((1)^s)^3*s^2*p^2*w1-112*L^2*r1*b-336*(-1)^s*L^2*w1-20*b^3*v*((1)^s)^3*s^2*p^2*r2-448*(-1)^s*L^2*r2*b+112*(1)^s*L^2*r1*b+13*b^2*s^2*p^2*w1+22*b^2*s^2*p^2*w2-4*b^3*s^2*p^2*r215*b^3*v*s^2*p^2*r1-65*b^2*v*s^2*p^2*w1110*b^2*v*s^2*p^2*w2+3*b^3*s^2*p^2*r1+20*b^3*v*s^2*p^2*r2+448*L^2*r2*b13*b^2*((-1)^s)^3*s^2*p^2*w1+4*b^3*((-1)^s)^3*s^2*p^2*r2-3*b^3*((1)^s)^3*s^2*p^2*r1-22*b^2*((-1)^s)^3*s^2*p^2*w2)/(-1+v^2)/L/b/s/p*t; K(2,1)=0; K(2,2)=0; K(2,3)=-1/630*E*(504*(-1)^s*v*w1*L^2+105*v*w1*((-1)^s)^3*L^2-504*(1)^s*v*w2*L^2+84*(-1)^s*v*r2*b*L^2-105*w1*((-1)^s)^3*L^2-7*b^3*((1)^s)^3*s^2*p^2*r2+15*b^3*((-1)^s)^3*s^2*p^2*r1+120*b^2*((1)^s)^3*s^2*p^2*w1+27*b^2*((-1)^s)^3*s^2*p^2*w2+504*v*w2*L^2-84*v*r2*b*L^2609*v*w1*L^2+105*L^2*w1-15*b^3*s^2*p^2*r1-120*b^2*s^2*p^2*w127*b^2*s^2*p^2*w2+7*b^3*s^2*p^2*r2)/b/L^2/(-1+v^2)*t; K(2,4)=-1/1260*E*b*(-28*(-1)^s*v*r2*L^2-3*b^2*((-1)^s)^3*s^2*p^2*r2+168*(1)^s*v*r1*L^2+12*b*s^2*p^2*((-1)^s)^3*w2+5*b^2*((1)^s)^3*s^2*p^2*r1+30*b*s^2*p^2*((-1)^s)^3*w1-168*v*r1*L^2+28*v*r2*L^2THỰC HIỆN : LÊ HIỀN ANH Trang 55 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG 5*b^2*s^2*p^2*r1-30*b*s^2*p^2*w112*b*s^2*p^2*w2+3*b^2*s^2*p^2*r2)/L^2/(1+v)/(v-1)*t; K(2,5)=0; K(2,6)=0; K(2,7)=-1/630*E*(504*(-1)^s*v*w2*L^2-504*(-1)^s*v*w1*L^2-105*v*((1)^s)^3*w2*L^2-84*(-1)^s*v*r2*b*L^2+6*b^3*((-1)^s)^3*s^2*p^2*r1+105*((1)^s)^3*w2*L^2+27*b^2*((-1)^s)^3*s^2*p^2*w1+36*b^2*((-1)^s)^3*s^2*p^2*w27*b^3*((-1)^s)^3*s^2*p^2*r2+504*v*w1*L^2+84*v*r2*b*L^2-399*v*w2*L^2105*L^2*w2-36*b^2*s^2*p^2*w2+7*b^3*s^2*p^2*r2-27*b^2*s^2*p^2*w16*b^3*s^2*p^2*r1)/b/L^2/(-1+v^2)*t; K(2,8)=1/1260*E*(-3*b^3*((-1)^s)^3*s^2*p^2*r2+14*b^2*((1)^s)^3*s^2*p^2*w2+14*b^2*((-1)^s)^3*s^2*p^2*w1+3*b^3*((1)^s)^3*s^2*p^2*r1-56*(-1)^s*v*r2*b*L^2+168*(-1)^s*v*w2*L^2+28*(1)^s*v*r1*L^2*b-168*(-1)^s*v*w1*L^2-28*v*r1*L^2*b168*v*w2*L^2+56*v*r2*b*L^2+168*v*w1*L^2-14*b^2*s^2*p^2*w2+3*b^3*s^2*p^2*r214*b^2*s^2*p^2*w1-3*b^3*s^2*p^2*r1)/L^2/(-1+v^2)*t; K(3,1)=K(1,3); K(3,2)=K(2,3); K(3,3)=1/1260*E*(-78*b^2*s^2*p^2*u2*((-1)^s)^3*L^2-72*b^3*s^3*p^3*v2*((1)^s)^3*L-132*b^2*((-1)^s)^3*s^2*p^2*u1*L^2-240*b^3*s^3*p^3*v1*((1)^s)^3*L+390*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2-1008*(-1)^s*v*v1*L^3*b*s*p1008*(-1)^s*v*v2*L^3*b*s*p-210*((-1)^s)^3*v2*L^3*b*s*p+210*((1)^s)^3*v1*L^3*b*s*p-180*b^2*s^2*p^2*v*u1*((-1)^s)^3*L^2+210*v*((1)^s)^3*v2*L^3*b*s*p-210*v*((-1)^s)^3*v1*L^3*b*s*p-2016*(1)^s*L^4*u1+1218*v*v1*L^3*b*s*p+2016*(-1)^s*L^4*u22016*L^4*u2+2016*L^4*u1+798*v*v2*L^3*b*s*p210*v1*L^3*b*s*p+210*v2*L^3*b*s*p+132*b^2*s^2*p^2*u1*L^2+78*b^2*s^2*p^2*u2* L^2+240*b^3*s^3*p^3*v1*L390*b^2*s^2*p^2*v*u2*L^2+72*b^3*s^3*p^3*v2*L+180*b^2*s^2*p^2*v*u1*L^2)/L^3/ b^2/p/s/(-1+v^2)*t; K(3,4)=1/1260*E*(55*b^2*s^2*p^2*v*u2*L^2+11*b^2*s^2*p^2*u2*L^2+30*b^3*s^3*p^3*v1*L+168*L^4*u 1+14*b^3*s^3*p^3*v2*L-168*L^4*u2+168*(-1)^s*L^4*u2+11*b^2*s^2*p^2*((1)^s)^3*u1*L^2-11*b^2*s^2*p^2*u2*((-1)^s)^3*L^2-55*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2+55*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2+168*v*v2*L^3*b*s*p30*b^3*s^3*p^3*v1*((-1)^s)^3*L-14*b^3*s^3*p^3*v2*((-1)^s)^3*L11*b^2*s^2*p^2*u1*L^2-168*(-1)^s*L^4*u1-168*(1)^s*v*v2*L^3*b*s*p+55*b^2*s^2*p^2*v*u1*L^2)/(-1+v^2)/L^3/b/s/p*t; K(3,5)=1/2520*E*(-4032*(-1)^s*L^2*w2+4032*L^2*w2-4032*L^2*w1+270*b^2*v*((1)^s)^3*s^2*p^2*w2+110*b^3*v*((-1)^s)^3*s^2*p^2*r1+780*b^2*v*((1)^s)^3*s^2*p^2*w1-336*L^2*r1*b+4032*(-1)^s*L^2*w1-65*b^3*v*((1)^s)^3*s^2*p^2*r2+336*(-1)^s*L^2*r2*b+336*(1)^s*L^2*r1*b+156*b^2*s^2*p^2*w1+54*b^2*s^2*p^2*w2-13*b^3*s^2*p^2*r2110*b^3*v*s^2*p^2*r1-780*b^2*v*s^2*p^2*w1270*b^2*v*s^2*p^2*w2+22*b^3*s^2*p^2*r1+65*b^3*v*s^2*p^2*r2-336*L^2*r2*b156*b^2*s^2*p^2*((-1)^s)^3*w1+13*b^3*((-1)^s)^3*s^2*p^2*r2-22*b^3*((1)^s)^3*s^2*p^2*r1-54*b^2*s^2*p^2*((-1)^s)^3*w2)/b^2/p/s/(-1+v^2)/L*t; K(3,6)=-1/630*E*(-504*(-1)^s*v*w2*L^2+504*(-1)^s*v*w1*L^2+84*(1)^s*v*r1*L^2*b-105*v*w1*((-1)^s)^3*L^2-6*b^3*((1)^s)^3*s^2*p^2*r2+7*b^3*((-1)^s)^3*s^2*p^2*r1+105*w1*((1)^s)^3*L^2+36*b^2*s^2*p^2*((-1)^s)^3*w1+27*b^2*s^2*p^2*((1)^s)^3*w2+504*v*w2*L^2-84*v*r1*L^2*b-399*v*w1*L^2-105*L^2*w1THỰC HIỆN : LÊ HIỀN ANH Trang 56 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG 27*b^2*s^2*p^2*w2-36*b^2*s^2*p^2*w1+6*b^3*s^2*p^2*r27*b^3*s^2*p^2*r1)/b/L^2/(-1+v^2)*t; K(3,7)=1/1260*E*(54*b^3*s^3*p^3*v1*L135*b^2*s^2*p^2*v*u2*L^2+54*b^3*s^3*p^3*v2*L+135*b^2*s^2*p^2*v*u1*L^22016*L^4*u1+1008*(-1)^s*v*v1*L^3*b*s*p-1008*v*v1*L^3*b*s*p+2016*L^4*u227*b^2*s^2*p^2*u2*((-1)^s)^3*L^2+27*b^2*s^2*p^2*((-1)^s)^3*u1*L^2+1008*(1)^s*v*v2*L^3*b*s*p+135*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^254*b^3*s^3*p^3*v1*((-1)^s)^3*L-54*b^3*s^3*p^3*v2*((-1)^s)^3*L-2016*(1)^s*L^4*u2+2016*(-1)^s*L^4*u1-135*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2+27*b^2*s^2*p^2*u2*L^2-1008*v*v2*L^3*b*s*p27*b^2*s^2*p^2*u1*L^2)/L^3/b^2/p/s/(1+v)/(v-1)*t; K(3,8)=-1/2520*E*(-65*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2+13*b^2*s^2*p^2*u2*L^2+28*b^3*s^3*p^3*v1*L+336*L^4*u265*b^2*s^2*p^2*v*u2*L^2+24*b^3*s^3*p^3*v2*L-336*v*v1*L^3*b*s*p336*L^4*u1+336*(-1)^s*v*v1*L^3*b*s*p+13*b^2*s^2*p^2*((-1)^s)^3*u1*L^213*b^2*s^2*p^2*u2*((-1)^s)^3*L^2-28*b^3*s^3*p^3*v1*((-1)^s)^3*L-336*(1)^s*L^4*u2+65*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2+65*b^2*s^2*p^2*v*u1*L^224*b^3*s^3*p^3*v2*((-1)^s)^3*L-13*b^2*s^2*p^2*u1*L^2+336*(1)^s*L^4*u1)/p/s/b/(1+v)/L^3/(v-1)*t; K(4,1)=K(1,4); K(4,2)=K(2,4); K(4,3)=K(3,4); K(4,4)=1/2520*E*(-4*b^2*s^2*p^2*u2*((-1)^s)^3*L^2+10*b^3*s^3*p^3*v1*L20*b^2*s^2*p^2*v*u2*L^2+112*v*v2*L^3*b*s*p+20*b^2*s^2*p^2*v*u1*L^2-448*(1)^s*L^4*u1-336*(-1)^s*v*v1*L^3*b*s*p+336*v*v1*L^3*b*s*p448*L^4*u2+448*L^4*u14*b^2*s^2*p^2*u1*L^2+4*b^2*s^2*p^2*u2*L^2+6*b^3*s^3*p^3*v2*L+4*b^2*s^2*p^2* ((-1)^s)^3*u1*L^2+448*(-1)^s*L^4*u2-112*(-1)^s*v*v2*L^3*b*s*p10*b^3*s^3*p^3*v1*((-1)^s)^3*L-6*b^3*s^3*p^3*v2*((1)^s)^3*L+20*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2-20*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2)/L^3/p/s/(-1+v^2)*t; K(4,5)=1/2520*E*(-336*(-1)^s*L^2*w2+336*L^2*w2-336*L^2*w1+65*b^2*v*((1)^s)^3*s^2*p^2*w2+20*b^3*v*((-1)^s)^3*s^2*p^2*r1+110*b^2*v*((1)^s)^3*s^2*p^2*w1-448*L^2*r1*b+336*(-1)^s*L^2*w1-15*b^3*v*((1)^s)^3*s^2*p^2*r2-112*(-1)^s*L^2*r2*b+448*(1)^s*L^2*r1*b+22*b^2*s^2*p^2*w1+13*b^2*s^2*p^2*w2-3*b^3*s^2*p^2*r220*b^3*v*s^2*p^2*r1-110*b^2*v*s^2*p^2*w165*b^2*v*s^2*p^2*w2+4*b^3*s^2*p^2*r1+15*b^3*v*s^2*p^2*r2+112*L^2*r2*b22*b^2*s^2*p^2*((-1)^s)^3*w1+3*b^3*((-1)^s)^3*s^2*p^2*r2-4*b^3*((1)^s)^3*s^2*p^2*r1-13*b^2*((-1)^s)^3*s^2*p^2*w2)/(-1+v^2)/L/b/s/p*t; K(4,6)=-1/1260*E*(-168*(-1)^s*v*w2*L^2-28*(-1)^s*v*r2*b*L^2+56*(1)^s*v*r1*L^2*b+168*(-1)^s*v*w1*L^2+3*b^3*((1)^s)^3*s^2*p^2*r1+14*b^2*s^2*p^2*((-1)^s)^3*w1+14*b^2*((1)^s)^3*s^2*p^2*w2-3*b^3*((-1)^s)^3*s^2*p^2*r256*v*r1*L^2*b+28*v*r2*b*L^2+168*v*w2*L^2-168*v*w1*L^2-14*b^2*s^2*p^2*w214*b^2*s^2*p^2*w1-3*b^3*s^2*p^2*r1+3*b^3*s^2*p^2*r2)/L^2/(1+v)/(v-1)*t; K(4,7)=1/2520*E*(-13*b^2*s^2*p^2*u2*((-1)^s)^3*L^2+336*L^4*u2-336*L^4*u1336*(-1)^s*L^4*u2+28*b^3*s^3*p^3*v2*L+65*b^2*s^2*p^2*v*u1*L^2+13*b^2*((1)^s)^3*s^2*p^2*u1*L^2+336*(-1)^s*L^4*u1-24*b^3*s^3*p^3*v1*((-1)^s)^3*L28*b^3*s^3*p^3*v2*((-1)^s)^3*L+65*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2336*v*v2*L^3*b*s*p-65*b^2*s^2*p^2*v*u2*L^2-65*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2THỰC HIỆN : LÊ HIỀN ANH Trang 57 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG 13*b^2*s^2*p^2*u1*L^2+24*b^3*s^3*p^3*v1*L+13*b^2*s^2*p^2*u2*L^2+336*(1)^s*v*v2*L^3*b*s*p)/(-1+v^2)/L^3/b/s/p*t; K(4,8)=-1/2520*E*(15*b^2*s^2*p^2*v*u2*L^2+6*b^3*s^3*p^3*v2*L+56*v*v2*L^3*b*s*p+112*(1)^s*L^4*u2-112*(-1)^s*L^4*u1-56*(-1)^s*v*v1*L^3*b*s*p+56*v*v1*L^3*b*s*p112*L^4*u2+112*L^4*u1+6*b^3*s^3*p^3*v1*L-3*b^2*s^2*p^2*u2*((-1)^s)^3*L^23*b^2*s^2*p^2*u1*L^2+3*b^2*((1)^s)^3*s^2*p^2*u1*L^2+15*b^2*s^2*p^2*v*u1*L^2-56*(-1)^s*v*v2*L^3*b*s*p6*b^3*s^3*p^3*v1*((-1)^s)^3*L-6*b^3*s^3*p^3*v2*((1)^s)^3*L+3*b^2*s^2*p^2*u2*L^2+15*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^215*b^2*s^2*p^2*v*u1*((-1)^s)^3*L^2)/L^3/p/s/(-1+v^2)*t; K(5,1)=K(1,5); K(5,2)=K(2,5); K(5,3)=K(3,5); K(5,4)=K(4,5); K(5,5)=0; K(5,6)=0; K(5,7)=1/2520*E*(4032*(-1)^s*L^2*w2-4032*L^2*w2+4032*L^2*w1+360*b^2*v*((1)^s)^3*s^2*p^2*w2+65*b^3*v*((-1)^s)^3*s^2*p^2*r1+270*b^2*v*((1)^s)^3*s^2*p^2*w1+336*L^2*r1*b-4032*(-1)^s*L^2*w1-110*b^3*v*((1)^s)^3*s^2*p^2*r2-336*(-1)^s*L^2*r2*b-336*(1)^s*L^2*r1*b+54*b^2*s^2*p^2*w1-264*b^2*s^2*p^2*w2-22*b^3*s^2*p^2*r265*b^3*v*s^2*p^2*r1-270*b^2*v*s^2*p^2*w1360*b^2*v*s^2*p^2*w2+13*b^3*s^2*p^2*r1+110*b^3*v*s^2*p^2*r2+336*L^2*r2*b54*b^2*((-1)^s)^3*s^2*p^2*w1+22*b^3*((-1)^s)^3*s^2*p^2*r2-13*b^3*((1)^s)^3*s^2*p^2*r1+264*b^2*((-1)^s)^3*s^2*p^2*w2)/(1+v)/L/p/s/(v-1)/b^2*t; K(5,8)=-1/2520*E*(336*(-1)^s*L^2*w2-336*L^2*w2+336*L^2*w1+110*b^2*v*((1)^s)^3*s^2*p^2*w2+15*b^3*v*((-1)^s)^3*s^2*p^2*r1+65*b^2*v*((1)^s)^3*s^2*p^2*w1-112*L^2*r1*b-336*(-1)^s*L^2*w1-20*b^3*v*((1)^s)^3*s^2*p^2*r2-448*(-1)^s*L^2*r2*b+112*(1)^s*L^2*r1*b+13*b^2*s^2*p^2*w1+22*b^2*s^2*p^2*w2-4*b^3*s^2*p^2*r215*b^3*v*s^2*p^2*r1-65*b^2*v*s^2*p^2*w1110*b^2*v*s^2*p^2*w2+3*b^3*s^2*p^2*r1+20*b^3*v*s^2*p^2*r2+448*L^2*r2*b13*b^2*((-1)^s)^3*s^2*p^2*w1+4*b^3*((-1)^s)^3*s^2*p^2*r2-3*b^3*((1)^s)^3*s^2*p^2*r1-22*b^2*((-1)^s)^3*s^2*p^2*w2)/(1+v)/L/b/p/s/(v-1)*t; K(6,1)=K(1,6); K(6,2)=K(2,6); K(6,3)=K(3,6); K(6,4)=K(4,6); K(6,5)=K(5,6); K(6,6)=0; K(6,7)=-1/630*E*(-84*(-1)^s*v*r1*L^2*b-504*(-1)^s*v*w1*L^2+105*v*((1)^s)^3*w2*L^2+504*(-1)^s*v*w2*L^2+27*b^2*((-1)^s)^3*s^2*p^2*w1+120*b^2*((1)^s)^3*s^2*p^2*w2-105*((-1)^s)^3*w2*L^2-15*b^3*((1)^s)^3*s^2*p^2*r2+7*b^3*((-1)^s)^3*s^2*p^2*r1+504*v*w1*L^2609*v*w2*L^2+84*v*r1*L^2*b+105*L^2*w2+15*b^3*s^2*p^2*r2-120*b^2*s^2*p^2*w227*b^2*s^2*p^2*w1-7*b^3*s^2*p^2*r1)/b/L^2/(-1+v^2)*t; K(6,8)=1/1260*E*b*(12*b*s^2*p^2*((-1)^s)^3*w1+30*b*s^2*p^2*((1)^s)^3*w2+3*b^2*((-1)^s)^3*s^2*p^2*r1-5*b^2*((-1)^s)^3*s^2*p^2*r2-168*(1)^s*v*r2*L^2+28*(-1)^s*v*r1*L^2-28*v*r1*L^2-3*b^2*s^2*p^2*r1- THỰC HIỆN : LÊ HIỀN ANH Trang 58 Luận văn tốt nghiệp cao hoc ngành XD DD&CN K10 GVHD:PGS.TS.CHU QUỐC THẮNG 30*b*s^2*p^2*w2-12*b*s^2*p^2*w1+168*v*r2*L^2+5*b^2*s^2*p^2*r2)/L^2/(1+v^2)*t; K(7,1)=K(1,7); K(7,2)=K(2,7); K(7,3)=K(3,7); K(7,4)=K(4,7); K(7,5)=K(5,7); K(7,6)=K(6,7); K(7,7)=1/1260*E*(-2016*(-1)^s*L^4*u1+2016*(-1)^s*L^4*u2+798*v*v1*L^3*b*s*p2016*L^4*u2+2016*L^4*u1+132*b^2*s^2*p^2*u2*((-1)^s)^3*L^2240*b^3*s^3*p^3*v2*((-1)^s)^3*L+1218*v*v2*L^3*b*s*p-72*b^3*s^3*p^3*v1*((1)^s)^3*L+78*b^2*((-1)^s)^3*s^2*p^2*u1*L^2+180*b^2*s^2*p^2*v*u2*((1)^s)^3*L^2+210*v1*L^3*b*s*p-210*v2*L^3*b*s*p-1008*(-1)^s*v*v1*L^3*b*s*p1008*(-1)^s*v*v2*L^3*b*s*p-210*((-1)^s)^3*v1*L^3*b*s*p390*b^2*s^2*p^2*v*u1*((-1)^s)^3*L^2+210*((-1)^s)^3*v2*L^3*b*s*p78*b^2*s^2*p^2*u1*L^2-132*b^2*s^2*p^2*u2*L^2+72*b^3*s^3*p^3*v1*L180*b^2*s^2*p^2*v*u2*L^2+240*b^3*s^3*p^3*v2*L+390*b^2*s^2*p^2*v*u1*L^2210*v*((-1)^s)^3*v2*L^3*b*s*p+210*v*((-1)^s)^3*v1*L^3*b*s*p)/L^3/b^2/p/s/(1+v^2)*t; K(7,8)=-1/1260*E*(-168*(-1)^s*v*v1*L^3*b*s*p+168*v*v1*L^3*b*s*p168*L^4*u2+168*L^4*u1-168*(-1)^s*L^4*u1-11*b^2*s^2*p^2*u2*((1)^s)^3*L^2+14*b^3*s^3*p^3*v1*L+11*b^2*((-1)^s)^3*s^2*p^2*u1*L^2+168*(1)^s*L^4*u2+55*b^2*s^2*p^2*v*u1*L^2-14*b^3*s^3*p^3*v1*((-1)^s)^3*L30*b^3*s^3*p^3*v2*((-1)^s)^3*L+11*b^2*s^2*p^2*u2*L^2+30*b^3*s^3*p^3*v2*L11*b^2*s^2*p^2*u1*L^2+55*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^255*b^2*s^2*p^2*v*u1*((-1)^s)^3*L^2-55*b^2*s^2*p^2*v*u2*L^2)/(1+v^2)/L^3/b/s/p*t; K(8,1)=K(1,8); K(8,2)=K(2,8); K(8,3)=K(3,8); K(8,4)=K(4,8); K(8,5)=K(5,8); K(8,6)=K(6,8); K(8,7)=K(7,8); K(8,8)=1/2520*E*(-4*b^2*s^2*p^2*u2*((-1)^s)^3*L^2+6*b^3*s^3*p^3*v1*L20*b^2*s^2*p^2*v*u2*L^2+336*v*v2*L^3*b*s*p+20*b^2*s^2*p^2*v*u1*L^2-448*(1)^s*L^4*u1-112*(-1)^s*v*v1*L^3*b*s*p+112*v*v1*L^3*b*s*p448*L^4*u2+448*L^4*u14*b^2*s^2*p^2*u1*L^2+4*b^2*s^2*p^2*u2*L^2+10*b^3*s^3*p^3*v2*L+4*b^2*s^2*p^2 *((-1)^s)^3*u1*L^2+448*(-1)^s*L^4*u2-336*(-1)^s*v*v2*L^3*b*s*p6*b^3*s^3*p^3*v1*((-1)^s)^3*L-10*b^3*s^3*p^3*v2*((1)^s)^3*L+20*b^2*s^2*p^2*v*u2*((-1)^s)^3*L^2-20*b^2*s^2*p^2*v*u1*((1)^s)^3*L^2)/(1+v)/p/s/L^3/(v-1)*t; K=K/2; THỰC HIỆN : LÊ HIỀN ANH Trang 59 ... .51 KHẢO SÁT TẤM KHÔNG SƯỜN .52 KHẢO SÁT ẢNH HƯỞNG CỦA CHIỀU CAO SƯỜN ĐẾN TẦN SỐ DAO ĐỘNG RIÊNG VÀ MODE SHAPE TƯƠNG ỨNG 53 KHẢO SÁT ẢNH HƯỞNG CỦA KHOẢNG CÁCH CÁC SƯỜN... TOÁN TĨNH CÓ CÙNG CƯỜNG ĐỘ (VÍ DỤ 5) 45 TÌM CHUYỂN VỊ TẤM CÓ SƯỜN CHỊU TÁC DỤNG TẢI TRỌNG BẤT KỲ(VÍ DỤ 6) 49 CHƯƠNG 4: KHẢO SÁT ẢNH HƯỞNG CỦA SƯỜN GIA CØNG ĐẾN DAO ĐỘNG CỦA TẤM CHỬ NHẬT... nghiệp động viên chia với khó khăn lúc thực luận văn Lê Hiền Anh TÓM TẮT Phương pháp dải hữu hạn (FSM) phương pháp phát triển từ phương pháp phần tử hữu hạn (FEM) áp dụng riêng cho kết cấu có hình