Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 81 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
81
Dung lượng
1,69 MB
Nội dung
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC DÂN LẬP HẢI PHỊNG NGUYỄN BÁ THÀNH TÍNH TỐN NỘILỰCVÀCHUYỂNVỊCỦAHỆKHUNGCÓXÉTĐẾNBIẾNDẠNGTRƯỢTNGANGBẰNGPHƯƠNGPHÁPPHẦNTỬHỮUHẠN LUẬN VĂN THẠC SĨ KỸ THUẬT CHUYÊN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP MÃ SỐ: 60.58.02.08 NGƯỜI HƯỚNG DẪN KHOA HỌC TS ĐỖ TRỌNG QUANG HẢI PHỊNG - 2018 LỜI CAM ĐOAN Tơi xin cam đoan đề tài “Tính tốn nộilựcchuyểnvịhệkhungcóxétđếnbiếndạngtrượtngangphươngphápphầntửhữu hạn” đề tài nghiên cứu Những số liệu sử dụng phân tích luận văn hồn tồn trung thực có nguồn gốc rõ ràng Các kết nghiên cứu luận văn tơi tự tìm hiểu, phân tích cách trung thực, khách quan với thực tiễn Hải Phòng, ngày tháng năm 2018 Tác giả Nguyễn Bá Thành i LỜI CẢM ƠN Qua trình học tập nghiên cứu, giúp đỡ, thầy cô giáo trường Đại học Dân lập Hải Phòng, Khoa Đào tạo sau Đại học, tơi hồn thành chương trình học tập nghiên cứu luận văn với đề tài “Tính tốn nộilựcchuyểnvịhệkhungcóxétđếnbiếndạngtrượtngangphươngphápphầntửhữu hạn” Tôi xin trân trọng cảm ơn TS Đỗ Trọng Quang tạo điều kiện tận tình giúp đỡ tơi q trình nghiên cứu đề tài Bên cạnh đó, tơi xin chân thành gửi lời cảm ơn đến bạn bè, đồng nghiệp giúp đỡ tơi tìm hiểu nghiên cứu tài liệu liên quan q trình hồn thành luận văn tốt nghiệp Mặc dù nỗ lựccố gắng để hồn thành luận văn thời gian kiến thức hạn chế nên luận văn khơng tránh khỏi mặt tồn định Tôi mong nhận đóng góp ý kiến từ quý Thầy Cơ để hồn thiện tốt luận văn Xin trân trọng cảm ơn! Hải Phòng, ngày tháng Tác giả Nguyễn Bá Thành ii năm 2018 MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỞ ĐẦU CHƯƠNG BÀI TOÁNCƠ HỌC KẾT CẤU VÀ CÁC PHƯƠNGPHÁP GIẢI 1.1 Bài toán học kết cấu 1.2 Các phươngpháp giải 1.2.1 Phươngpháplực 1.2.2 Phươngphápchuyểnvị 1.2.3 Phươngpháp hỗn hợp phươngpháp liên hợp 1.2.4 Phươngpháp sai phânhữuhạn 1.2.5 Phươngpháp hỗn hợp sai phân – biếnphân CHƯƠNG 2: PHƯƠNGPHÁPPHẦNTỬHỮUHẠN 2.1 Phươngphápphầntửhữuhạn 2.1.1 Nội dung phươngphápphầntử hữa hạn theo mơ hình chuyểnvị 2.1.2 Cách xây dựng ma trận độ cứng phầntử chịu uốn 26 2.1.3 Cách xây dựng ma trận độ cứng tổng thể kết cấu 29 CHƯƠNG TÍNH TỐN KHUNG PHẲNG CÓXÉTĐẾNBIẾNDẠNGTRƯỢTNGANGBẰNGPHƯƠNGPHÁPPHẦNTỬHỮUHẠN 34 3.1 Lý thuyết dầm cóxétbiếndạngtrượtngang 34 3.2 Phươngphápphầntửhữuhạn tốn khungcóxétđếnbiếndạngtrượtngang 40 3.2.1 Hàm nội suy phầntử 41 3.2.1.1 Hàm nội suy chuyểnvị góc xoay hai nút đầu phầntử 41 3.2.1.2 Hàm nội suy lực cắt hai nút đầu phầntử 43 3.2.2 Ma trận độ cứng phầntử 44 iii 3.2.2.1 Trường hợp không xétbiếndạngtrượtngang 44 3.2.2.2 Trường hợp cóxétđếnbiếndạngtrượtngang 46 3.2.3 Ma trận độ cứng tổng thể 47 3.3 Các ví dụ tính tốn khungcóxétđếnbiếndạngtrượtngang 48 KẾT LUẬN VÀ KIẾN NGHỊ 72 KẾT LUẬN 72 KIẾN NGHỊ 73 Danh mục tài liệu tham khảo 74 iv MỞ ĐẦU Phươngphápphầntửhữuhạnphươngpháp xây dựng dựa ý tưởng rời rạc hóa cơng trình thành phầntử nhỏ (số phầntửhữu hạn) Các phầntử nhỏ nối lại với thông qua phương trình cân phương trình liên tục Để giải tốn học kết cấu, tiếp cận phươngpháp theo ba mơ hình gồm: Mơ hình chuyển vị, xem chuyểnvị đại lượng cần tìm hàm nội suy biểu diễn gần dạngphân bố chuyểnvịphần tử; Mô hình cân bằng, hàm nội suy biểu diễn gần dạngphân bố ứng suất hay nộilựcphầntử mơ hình hỗn hợp, coi đại lượng chuyểnvị ứng suất hai yếu tố độc lập riêng biệt Các hàm nội suy biểu diễn gần dạngphân bố chuyểnvị lẫn ứng suất phầntử Bài toán học kết cấu nói chung xây dựng theo bốn đường lối là: Xây dựng phương trình viphân cân phân tố; Phươngpháp lượng; Phươngpháp nguyên lý công ảo Phươngpháp sử dụng trực tiếp Phương trình Lagrange Các phươngpháp giải gồm có: Phươngpháp coi xác như, phươngpháp lực, phươngphápchuyển vị, phươngpháp hỗn hợp, phươngpháp liên hợp phươngpháp gần như: Phươngphápphầntửhữu hạn, phươngpháp sai phânhữu hạn, phươngpháp hỗn hợp sai phân - biếnphân Trong phươngphápphầntửhữuhạn dùng rộng rãi phổ biến tốn kỹ thuật nói chung tốn học kết cấu nói riêng Đối tượng, phươngpháp phạm vi nghiên cứu đề tài Trong luận văn này, tác giả sử dụng phươngphầntửhữuhạn để xây dựng giải toánkhung phẳng chịu uốn cóxétđếnbiếndạngtrượtngang chịu tác dụng tải trọng tĩnh Mục đích nghiên cứu đề tài “Tính tốn khung phẳng chịu uốn cóxétbiếndạngtrượtngang chịu tác dụng tải trọng tĩnhphươngphápphầntửhữu hạn” Nhiệm vụ nghiên cứu đề tài Tìm hiểu giới thiệu phươngpháp giải toán học kết cấu Trình bày lý thuyết dầm Euler - Bernoulli lý thuyết dầm cóxétđếnbiếndạngtrượtngang Trình bày phươngphápphầntửhữuhạn áp dụng để giải toánkhung phẳng, chịu tác dụng tải trọng tĩnh tập trung Lập chương trình máy tính điện tử cho tốn nêu CHƯƠNG BÀI TỐN CƠ HỌC KẾT CẤU VÀ CÁC PHƯƠNGPHÁP GIẢI Trong chương giới thiệu toán học kết cấu (bài toán tĩnh) phươngpháp giải thường dùng 1.1 Bài toán học kết cấu Bài toán học kết cấu nhằm xác định nộilựcchuyểnvịhệ thanh, tấm, vỏ tác dụng loại tải trọng, nhiệt độ, chuyểnvị cưỡng bức,…và chia làm hai loại: - Bài toántĩnh định: tốn có cấu tạo hình học bất biến hình đủ liên kết tựa với đất, liên kết xếp hợp lý, chịu loại tải trọng Để xác định nộilựcchuyểnvị cần dùng phương trình cân tĩnh học đủ; - Bài toán siêu tĩnh: toáncó cấu tạo hình học bất biến hình thừa liên kết (nội ngoại) chịu loại tải trọng, nhiệt độ, chuyểnvị cưỡng bức,…Để xác định nộilựcchuyểnvị ngồi phương trình cân ta phải bổ sung phương trình biếndạng Nếu tínhđến tận ứng suất, nói tốn học vật rắn biếndạngnói chung tốn học kết cấu nói riêng toán siêu tĩnh 1.2 Các phươngpháp giải Đã có nhiều phươngpháp để giải toán siêu tĩnh Hai phươngpháp truyền thống phươngpháplựcphươngphápchuyểnvị Khi sử dụng chúng thường phải giải hệphương trình đại số tuyến tính Số lượng phương trình tùy thuộc vào phươngphápphân tích Từphươngphápchuyểnvị ta có hai cách tính gần hay sử dụng H Cross G Kani Từ xuất máy tính điện tử, người ta bổ sung thêm phươngpháp số khác như: Phươngphápphầntửhữu hạn; Phươngpháp sai phânhữu hạn… 1.2.1 Phươngpháplực Trong hệ siêu tĩnh ta thay liên kết thừa lực chưa biết, giá trị chuyểnvịhệ tương ứng với vị trí phươnglực ẩn số thân lực nguyên nhân bên ngồi gây khơng Từ điều kiện ta lập hệphương trình đại số tuyến tính, giải hệ ta tìm ẩn số từ suy đại lượng cần tìm 1.2.2 Phươngphápchuyểnvị Khác với phươngpháp lực, phươngphápchuyểnvị lấy chuyểnvị nút làm ẩn Những chuyểnvị phải có giá trị cho phảnlực liên kết đặt thêm vào hệ thân chúng ngun nhân bên ngồi gây khơng Lập hệphương trình đại số tuyến tính thỏa mãn điều kiện giải hệ ta tìm ẩn, từ xác định đại lượng lại Hệphươngphápchuyểnvị giới hạn giải toán phụ thuộc vào số phầntử mẫu có sẵn 1.2.3 Phươngpháp hỗn hợp phươngpháp liên hợp Phươngpháp hỗn hợp, phươngpháp liên hợp kết hợp song song phươngpháplựcphươngphápchuyểnvị Trong phươngpháp ta chọn hệ theo phươngpháplực không loại bỏ hết liên kết thừa mà loại bỏ liên kết thuộc phận thích hợp với phươngpháp lực; chọn hệ theo phươngphápchuyểnvị không đặt đầy đủ liên kết phụ nhằm ngăn cản toànchuyểnvị nút mà đặt liên kết phụ nút thuộc phận thích hợp với phươngphápchuyểnvị Trường hợp đầu hệ siêu tĩnh, trường hợp sau hệ siêu động Trong hai cách nói trên, tốn ban đầu đưa hai toán độc lập: Một theo phươngpháplực theo phươngphápchuyểnvị 1.2.4 Phươngpháp sai phânhữuhạnPhươngpháp sai phânhữuhạn thay hệ liên tục mơ hình rời rạc, song hàm cần tìm (hàm mang đến cho phiếm hàm giá trị dừng), nhận giá trị gần số hữuhạn điểm miền tích phân, giá trị điểm trung gian xác định nhờ phươngpháp tích phânPhươngpháp cho lời giải số phương trình viphânchuyểnvịnộilực điểm nút Thông thường ta phải thay đạo hàm sai phân hàm nút Phương trình viphânchuyểnvịnộilực viết dạng sai phân nút, biểu thị quan hệchuyểnvị nút nút lân cận tác dụng ngoại lực 1.2.5 Phươngpháp hỗn hợp sai phân – biếnphân Kết hợp phươngpháp sai phân với phươngphápbiếnphân ta cóphươngpháp linh động hơn: Hoặc sai phân đạo hàm phương trình biếnphân sai phân theo phươngbiếnphân theo phương khác (đối với toán hai chiều) 96 768 - 768 - 768 768 - 96 - 96 16 96 K e - 96 96 0 - 0.0001 0.0001 96 0 - 96 0 - 0.0001 0.0001 16 0.0001 - 0.0001 0.0001 0 - 0.0001 0 - Ma trận độ cứng toàn dầm [K]: Ghép nối ma trận độ cứng phầntử [Ke] vào hệ tọa độ chung, ta ma trận độ cứng tổng thể toàn kết cấu [K(49x49)], khơng trình bày kích thước ma trận lớn - Véc tơ lực nút{F}: Trong ví dụ véc tơ cột 49 dòng, sau: Sau có [K] {F}: ta thành lập hệphương trình phươngphápphầntửhữu hạn, sau: [𝐾] = {F} Giải phương trình (e) ta nhận được: K 1F Theo ngơn ngữ lập trình Matlab ta viết: K \ F Kết chuyển vị, góc xoay nút Mơmen uốn dầm: Khi không xétđếnbiếndạngtrượt ngang, h/l=1/100 (npt=16) ta có x 10 -4 X: 15 Y: 0.0001453 x 10 -3 -1 -2 -5 -3 -10 -4 -15 -5 X: Y: -0.005864 X: Y: -0.001948 -20 10 -6 12 14 16 Hình 3.9a Đường độ võng cột 62 10 12 14 Hình 3.9b Đường độ võng dầm 16 0.12 X: 16 Y: 0.104 X: 16 Y: 0.1355 0.15 X: Y: 0.104 0.1 0.1 X: Y: 0.07251 0.08 0.05 0.06 0.04 -0.05 0.02 -0.1 -0.04 X: Y: -0.1302 -0.15 -0.02 X: Y: -0.03674 -0.2 10 12 14 16 10 12 14 16 Hình 3.10a Biểu đồ mơmen cột Hình 3.10b Biểu đồ mơmen dầm Khi xétđến ảnh hưởng biếndạngtrượt ngang, trường hợp (h/l=1/3) ta nhận kết sau: x 10 -3 -0.5 -0.002 -1 -1.5 -0.004 -2 -0.006 -2.5 -3 -0.008 -3.5 -4 -5 -0.01 X: Y: -0.004603 -4.5 X: Y: -0.01154 -0.012 10 12 14 16 10 12 14 16 Hình 3.11b Đường độ võng dầm Hình 3.11a Đường độ võng cột 0.12 0.15 X: Y: 0.104 0.1 0.08 X: 16 Y: 0.1335 0.1 X: 16 Y: 0.104 X: Y: 0.07448 0.05 0.06 0.04 -0.05 0.02 -0.1 -0.02 -0.04 X: Y: -0.1312 -0.15 X: Y: -0.03576 10 12 14 Hình 3.12a Biểu đồ mômen cột Nhận xét kết trên: 16 -0.2 10 12 14 16 Hình 3.12b Biểu đồ mơmen dầm Khi chia cột dầm thành 16 phần ta nhận kết trên, so sánh với kết theo lời giải giải tích có ta nhận sai số theo bảng sau: 63 BẢNG SO SÁNH MÔMEN UỐN TẠI CÁC TIẾT DIỆN CỘT VÀ DẦM Các tiết diện Lời giải số theo Lời giải số theo Chênh lệch phươngphápphươngpháp khơng cóxét cột dầm PTHH PTHH (h/l=1/3) BDTN (%) (h/l=1/100) Chân cột -0,0725 -0,0745 2,75 Giữa cột 0,0367 0,0367 Đầu cột -1,0400 -1,0400 Giữa dầm 0,1302 0,1312 0,76 Đầu phải dầm -0,1355 -0,1336 1,40 Ta thấy chênh lệch hai trường hợp khơng cóxétđến ảnh hưởng biếndạngtrượtngang không đáng kể Ví dụ 3.3: Khung siêu tĩnh bậc 2, hình 3.13 Xác định nộilựcchuyểnvịkhung chịu lực hình 2, độ cứng uốn EJ=const P Hình 3.13 Khung siêu tĩnh bậc SO DO DAM NGANG SO DO AN LUC CAT nqx2 SO DO AN LUC CAT 51 53 54 55 56 57 58 59 33 34 35 36 37 38 39 52 32 SO DO AN GOC XOAY nwx2 45 46 47 48 49 50 31 SO DO AN GOC XOAY SO DO COT PHAI SO DO AN CHUYENVI 44 25 26 27 28 29 30 SO DO AN CHUYENVI 24 43 nút3 nw3 nwx3 nqx3 SO DO NUT COT PHAI nw2 40 40 41 41 42 42 21 21 22 22 23 23 CHIEU DAI PHANTU 0 nút2 SO DO NUT DAM 20 SO DO AN LUC CAT 14 15 16 17 18 19 13 10 11 12 SO DO AN GOC XOAY SO DO AN CHUYENVI 3 2 SO DO NUT COT TRAI 1 SO DO COT TRAI nút1 nw1 nwx1 nqx1 4 3 CHIEU DAI PHANTU Hình 3.14 Sơ đồ rời rạc kết cấu Rời rạc hóa kết cấu dầm thành npt phần tử, hình 3.14 Các nút phầntử phải trùng với vị trí đặt lực tập trung, hay vị trí thay đổi tiết diện, chiều dài phầntử khác 64 Mỗi phầntửcó ẩn 𝑤1 , 1 , 𝑤2 , 2 , 𝑄1, 𝑄2 npt phầntử rời rạc tổng cộng có 6xnpt ẩn Khi chia khung thành phầntử thì: ta có ẩn chuyểnvị nwi, góc xoay nmxi, lực cắt nqxi sau: nw1 =[ 1 2 3 nwx1 =[5 10 11 nqx1 =[13 14 15 nw2 =[0 21 21 16 22 17 22 4] 18 23 12] 19 20] 23 0] nwx2 =[24 25 26 27 28 29 30 31] nqx2 =[32 33 34 35 36 37 38 39] nw3 =[0 40 40 41 41 42 42 43] nwx3 =[44 45 46 47 48 49 50 51] nqx3 =[52 53 54 55 56 57 58 59] Sau biết ẩn số thực dầm cột ta xây dựng ma trận độ cứng tổng thể khungcó kích thước K(nxn), ví dụ K(59X59) Bây xét điều kiện liên tục góc xoay phầntử Điều kiện liên tục góc xoay phầntử viết sau: dy Q dy Q i 5( npt 1) dx GF dx GF nut phantu truoc nut phantu sau i (a) Trong toán ta chia thành phầntử nên i=3(4-1)=9 điều kiện liên tục góc xoay phầntử Điều kiện biên viết sau: - Tại đầu ngàm chân cột trái 1, cột phải có góc xoay khơng: dy1 Q1 0 dx GF nut1 phantu 10 dy3 Q3 0 dx GF nut phantu 11 65 (b) Điều kiện mômen nút khớp bên trái góc xoay nút cứng giao cột dầm viết sau: Mômen nút cuối phầntử đầu cột momen nút đầu phầntử dầm khơng Góc xoay nút cuối phầntử đầu cột phải góc xoay nút cuối phầntử cuối dầm d y1 12 dx d y2 dQ2 0 GF dx nut dx2 GF dx nut1 dQ1 dy3 Q3 dy Q dx GF nut2 dx GF nut2 13 (c) Điều kiện chuyểnvịngang đầu cột trái 1, cột cột phải nhau: 14 y1 nut2 y3 nut2 Trong k(k=114) ẩn số tốn (có k ẩn số ), tổng số ẩn số tốn lúc (n+k), ma trận độ cứng phầntửlúc phải thêm k dòng k cột kích thước ma trận độ cứng tổng thể K n k,n k Chẳng hạnví dụ này, ta có n=98, k=24 tổng số ẩn toán n+k=59+14=73 ẩn Trong trường hợp ta xác định kích thước ma trận độ cứng tổng thể là: K[73x73] Như cuối ta thiết lập phương trình: K F 66 (d) F1 so hang n Fn đó: F ; 0 so hang k 1 1 n 1 2 k ẩn số toán Trong ví dụ 3.1 chia thành phần tử, h/l=1/100 (không xétbiếndạngtrượt ngang) ta có: - Ma trận độ cứng phầntử [Ke(6x6)], sau: 96 768 - 768 - 768 768 - 96 - 96 16 96 K e - 96 96 0 - 0.0001 0.0001 96 0 - 96 0 - 0.0001 0.0001 16 0.0001 - 0.0001 0.0001 0 - 0.0001 0 - Ma trận độ cứng toàn dầm [K]: Ghép nối ma trận độ cứng phầntử [Ke] vào hệ tọa độ chung, ta ma trận độ cứng tổng thể toàn kết cấu [K(122x122)], khơng trình bày kích thước ma trận lớn - Véc tơ lực nút{F}: Trong ví dụ véc tơ cột 122 dòng, sau: Sau có [K] {F}: ta thành lập hệphương trình phươngphápphầntửhữu hạn, sau: [𝐾] = {F} Giải phương trình (e) ta nhận được: K 1F Theo ngơn ngữ lập trình Matlab ta viết: K \ F 67 Kết chuyển vị, góc xoay nút Mơmen uốn dầm: - Khi không xétđếnbiếndạngtrượt ngang, h/l=1/100 Dưới lần đường độ võng biểu đồ mômen uốn cột dầm x 10 -0.005 -3 -1 -2 -0.01 -3 X: Y: -0.01564 -0.015 -4 -0.02 -5 X: 16 Y: -0.02088 -0.025 10 12 14 X: Y: -0.005178 -6 16 10 12 14 16 Hình 3.15b Đường độ võng dầm Hình 3.15a Đường độ võng cột trái 0.15 0.3 X: X: 16 Y: 0.1248 Y: 0.2501 0.25 0.1 0.2 0.15 0.05 0.1 0.05 0 -0.05 -0.05 X: Y: -0.1249 -0.1 -0.15 X: Y: -0.06258 -0.1 10 12 14 16 10 12 14 16 Hình 3.16b Biểu đồ mơmen dầm Hình 3.16a Biểu đồ mômen cột trái 0.15 X: 0.1 Y: 0.125 -0.005 0.05 X: Y: -0.01042 -0.01 -0.05 -0.015 -0.1 -0.02 X: 16 Y: -0.02088 -0.025 10 12 14 X: 16 Y: -0.1248 -0.15 -0.2 16 Hình 3.15c Đường độ võng cột phải 10 12 14 Hình 3.16c Biểu đồ mơmen cột phải 68 16 Nhận xét kết trên: - Khi không xétđến ảnh hưởng biếndạngtrượtngang (trường hợp h/l=1/100) Khi chia cột dầm thành 16 phần ta nhận kết trên, so sánh với kết xác theo lời giải giải tích ta nhận sai số theo bảng sau: BẢNG SO SÁNH MÔMEN UỐN TẠI CÁC TIẾT DIỆN CỘT VÀ DẦM Các tiết diện Lời giải số theo Lời giải cột 1,3 dầm phươngpháp xác Sai số % PTHH (h/l=1/100) Chân cột -0,2501 -0,2500 0,04 Giữa cột 0,1249 0,1250 0,08 Giữa dầm 0,0625 0,0625 Đầu phải dầm -0,1248 -0,1250 0,16 Đầu cột -0,1248 -0,1250 0,16 Chân cột 0,1250 0,1250 Ta thấy kết nhận gần trùng khớp với lời giải xác - Khi xétđến ảnh hưởng biếndạngtrượtngang (trường hợp h/l=1/3), ta nhận kết sau: 0 x 10 -3 -1 -0.005 -2 -0.01 -3 -4 -0.015 -5 -0.02 -6 X: Y: -0.0242 -7 X: 16 Y: -0.02712 -0.025 -0.03 -9 10 12 14 X: Y: -0.008128 -8 16 Hình 3.17a Đường độ võng cột trái 69 10 12 14 16 Hình 3.17b Đường độ võng dầm 0.15 0.3 X: Y: 0.2521 0.25 X: 16 Y: 0.1221 0.1 0.2 0.15 0.05 0.1 0.05 0 -0.05 -0.1 -0.15 X: Y: -0.06396 -0.05 X: Y: -0.124 10 12 14 16 -0.1 Hình 3.18a Biểu đồ mơmen cột trái 10 12 14 16 Hình 3.18b Biểu đồ mômen dầm 0.15 X: 0.1 Y: 0.1259 -0.005 0.05 -0.01 X: Y: -0.01332 -0.015 -0.05 -0.02 -0.1 X: 16 Y: -0.02712 -0.025 -0.03 10 12 14 X: 16 Y: -0.1221 -0.15 16 -0.2 10 12 14 16 Hình 3.18c Biểu đồ mơmen Hình 3.17c Đường độ võng cột phải cột phải BẢNG SO SÁNH MÔMEN UỐN TẠI CÁC TIẾT DIỆN CỘT VÀ DẦM Các tiết diện Lời giải số theo Lời giải số theo Chênh lệch phươngphápphươngpháp khơng cóxét cột 1,3 dầm PTHH PTHH BDTN (%) (h/l=1/100) (h/l=1/100) Chân cột -0,2501 -0,2521 0,80 Giữa cột 0,1249 0,1240 0,72 Giữa dầm 0,0625 0,0639 2,24 Đầu phải dầm -0,1248 -0,1221 2,16 Đầu cột -0,1248 -0,1221 2,16 Chân cột 0,1250 0,1259 0,72 70 Mômen thay đổi không đáng kể hai trường hợp khơng xétcóxétđến ảnh hưởng biếndạngtrượt ngang, bảng 3.3 BẢNG SO SÁNH ĐỘ VÕNG TẠI CÁC TIẾT DIỆN CỘT VÀ DẦM Các tiết diện Lời giải số theo Lời giải số theo Chênh lệch phươngphápphươngpháp khơng cóxét cột 1,3 dầm PTHH PTHH (h/l=1/3) BDTN (%) (h/l=1/100) Đầu cột 0,0208 0,0271 30,28 Giữa cột 0,0156 0,0242 55,12 Giữa dầm 0,0052 0,0081 55,76 Đầu cột 0,0208 0,0271 30,28 Giữa cột 0,0104 0,0133 27,88 Khi xétđến ảnh hưởng biếndạngtrượt độ võng tăng lên lớn, lớn dầm (55,76%), nhỏ cột (27,88%) 71 KẾT LUẬN VÀ KIẾN NGHỊ KẾT LUẬN Qua kết nghiên cứu từ chương, chương đến chương toánkhung phẳng cóxétđếnbiếndạngtrượtngang chịu tác dụng tải trọng tĩnh Tác giả rút kết luận sau: Bằngphươngphápphầntửhữu hạn, tác giả xác định nộilựcchuyểnvịkhung siêu tĩnh chịu tác dụng tải trọng tĩnhcó điều kiện biên khác Kết nộilựcchuyểnvị trùng khớp với kết nhận giải phươngphápcó tăng số lượng phầntử lên lớn 16 phầntử Ảnh hưởng việc xétbiếndạngtrượtngangđến kết nhận được, phụ thuộc vào toán, cách đặt tải, loại tải trọng điều kiện liên kết hay đối xứng , chẳng hạnví dụ 3.1 ta thấy kết nhận hai trường hợp khơng xétcóxétđến ảnh hưởng biếndạngtrượtngang (h/l=1/3) mômen uốn thay đổi tương đối lớn, hầu hết tăng đa số tiết diện, tăng lớn đầu phải dầm 2, tăng (47,08%) Nhưng ví dụ 3.2 3.3 ta thấy mômen thay đổi không đáng kể hai trường hợp, có khơng xétđến ảnh hưởng biếndạngtrượtngang Khi không xétđếnbiếndạngtrượtngang (h/l=1/100) hay (h0) kết hội tụ theo lý thuyết dầm truyền thống (lý thuyết dầm Euler Bernoulli) Khi rời rạc hóa kết cấu với số phầntử nhiều kết tiệm cận tới kết xác nhận từphươngpháp giải tích Đối với tốn khung chịu tải trọng tập trung để đạt nộilực xác cần chia dầm thành phầntử Trong ví dụ 3.3 cho thấy nộilực M khung hai trường hợp có không xétđến ảnh hưởng biếndạng trượt, thay đổi không đáng kể, chuyểnvị tiết diện khung thay đổi tương đối lớn, lớn tiết diện dầm (55,76 %) Điều cho thấy kết tốn khơng ảnh hưởng đáng kể đếntoán bền, lại ảnh hưởng lớn đếntoán cứng 72 KIẾN NGHỊ Sử dụng phươngphápphầntửhữuhạn để giải toán khác như: Dầm, khung, dàn, tấm, vỏ Nội dung nghiên cứu làm tài liệu tham khảo, nghiên cứu, học tập ứng dụng thực tế thiết kế kết cấu cơng trình 73 DANH MỤC TÀI LIỆU THAM KHẢO I Tiếng Việt [1] Lều Thọ Trình, Đỗ Văn Bình (2010), Cơ học cơng trình, Nhà xuất xây dựng [2] Lều Thọ Trình (2000), Cơ học kết cấu, Nhà xuất khoa học kỹ thuật [3] Nguyễn Văn Liên, Nguyễn Phương Thành, Đinh Trọng Bằng (2003), Giáo trình Sức bền vật liệu, Nhà xuất xây dựng, tái lần thứ 3, 330 trang [4] Nguyễn Văn Đạo (2001), Cơ học giải tích, Nhà xuất Đại học Quốc gia Hà nội, 337 trang [5] Nguyễn Trâm (2007), Phươngphápphầntửhữu hạn, Nhà xuất xây dựng [6] Đinh Văn Phong (2006), Phươngpháp số học, Nhà xuất khoa học kỹ thuật [7] Đồn Văn Duẩn (2015), Bài tốn học kết cấu dạng tổng quát, Tạp chí Xây dựng số 02 (Tr59-Tr61) [8] Đoàn Văn Duẩn (2015), Phươngpháp so sánh nghiên cứu nộilựcchuyểnvịhệ dầm, Tạp chí Xây dựng số 11 (Tr56-Tr58) [9] Đồn Văn Duẩn (2015), Tính tốn kết cấu khung chịu uốn phươngpháp so sánh, Tạp chí Xây dựng số 12 (Tr62-Tr64) [10] Timoshenko C.P, Voinópki- Krige X, (1971), Tấm Vỏ Người dịch, Phạm Hồng Giang, Vũ Thành Hải, Đoàn Hữu Quang, Nxb Khoa học kỹ thuật, Hà Nội IIi Tiếng Anh [11] Stephen P.Timoshenko-Jame M.Gere (1961), Theory of elastic stability, McGraw-Hill Book Company, Inc, New york – Toronto – London, 541 Tr 74 [12] William T.Thomson (1998), Theory of Vibration with Applications (Tái lần thứ 5) Stanley Thornes (Publishers) Ltd, 546 trang [13] Klaus – Jurgen Bathe (1996), Finite Element procedures Part one, Prentice – Hall International, Inc, 484 trang [14] Klaus – Jurgen Bathe (1996), Finite Element procedures Part two, Prentice – Hall International, Inc, 553 trang [15] Ray W.Clough, Joseph Penzien(1993), Dynamics of Structures (Tái lần thứ 2), McGraw-Hill Book Company, Inc, 738 trang [16] O.C Zienkiewicz-R.L Taylor (1991), The finite element method (four edition) Volume 2, McGraw-Hill Book Company, Inc, 807 trang [17] G.Korn-T.Korn (1961), Mathematical Handbook for sientists and Engineers, McGraw-Hill, New york (Bản dịch tiếng Nga, I.Bramovich chủ biên, Nhà xuất Nauka-Moscow, 1964) [18] Stephen P.Timoshenko-J Goodier (1970), Theory of elasticity, McGrawHill, New york (Bản dịch tiếng Nga, G Shapiro chủ biên, Nhà xuất Nauka-Moscow, 1979), 560 trang [19] D.R.J Owen, E.Hinton (1986), Finite Elements in Plasticity: Theory and Practice, Pineridge Press Lt [20] Lars Olovsson, Kjell Simonsson, Mattias Unosson (2006), Shear locking reduction in eight-node tri-linear solid finite elements, J ‘Computers @ Structures’,84, trg 476-484 [21] C.A.Brebbia, Techniques Theory J.C.F.Telles, L.C.Wrobel(1984), Boundary Element and Applications in Engineering Nxb Springer – Verlag.(Bản dịch tiếng Nga, 1987) [22] Chopra Anil K (1995) Dynamics of structures Prentice Hall, Englewood Cliffs, New – Jersey 07632 [23] Wilson Edward L Professor Emeritus of structural Engineering University of California at Berkeley (2002) Three – Dimensional Static and 75 Dynamic Analysis of structures, Inc Berkeley, California, USA Third edition, Reprint January [24] Wilson, E L., R L Taylor, W P Doherty and J Ghaboussi (1971) “Incompatible Displacement Models”, Proceedings, ORN Symposium on “Numerical and Computer Method in Structural Mechanics” University of Illinois, Urbana September Academic Press [25] Strang, G (1972) “Variational Crimes in the Finite Element Method” in “The Mathematical Foundations of the Finite Element Method” P.689 -710 (ed A.K Aziz) Academic Press [26] Irons, B M and O C Zienkiewicz (1968) “The isoparametric Finite Element System – A New Concept in Finite Element Analysis”, Proc Conf “Recent Advances in Stress Analysis” Royal Aeronautical Society London [27] Kolousek Vladimir, DSC Professor, Technical University, Pargue (1973) Dynamics in engineering structutes Butter worths London [28] Felippa Carlos A (2004) Introduction of finite element methods Department of Aerospace Engineering Sciences and Center for Aerospace Structures University of Colorado Boulder, Colorado 80309-0429, USA, Last updated Fall 76 ... TÍNH TỐN KHUNG PHẲNG CĨ XÉT ĐẾN BIẾN DẠNG TRƯỢT NGANG BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN 34 3.1 Lý thuyết dầm có xét biến dạng trượt ngang 34 3.2 Phương pháp phần tử hữu hạn toán khung có. .. phương pháp lực, phương pháp chuyển vị, phương pháp hỗn hợp, phương pháp liên hợp phương pháp gần như: Phương pháp phần tử hữu hạn, phương pháp sai phân hữu hạn, phương pháp hỗn hợp sai phân - biến. .. chuyển vị Sau luận văn trình nội dung phương pháp phần tử hữu hạn theo mơ hình chuyển vị 2.1.1 Nội dung phương pháp phần tử hữa hạn theo mơ hình chuyển vị Trong phương pháp phần tử hữu hạn - mơ