Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 74 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
74
Dung lượng
1,7 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 VĂN TRƢỜNG PHƢƠNG PHÁPPHẦNTỬHỮUHẠNĐỐIVỚIBÀITOÁNDẦMLIÊNTỤCCHỊUTẢITRỌNGTĨNHTẬPTRUNG Chuyên ngành: Kỹ thuật Xây dựng Công trình Dân dụng & Công nghiệp Mã số: 60.58.02.08 LUẬN VĂN THẠC SỸ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC: GS.TSKH HÀ HUY CƢƠNG Hải Phòng, 2017 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng Các số liệu, kết luận văn trung thực chưa công bố công trình khác Tác giả luận văn Nguyễn Văn Trƣờng LỜI CẢM ƠN Tác giả luận văn xin trân trọng bày tỏ lòng biết ơn sâu sắc GS.TSKH Hà Huy Cương ý tưởng khoa học độc đáo, bảo sâu sắc phươngpháp nguyên lý cực trị Gauss chia sẻ kiến thức học, toán học uyên bác Giáo sư Giáo sư tận tình giúp đỡ cho nhiều dẫn khoa học có giá trị thường xuyên động viên, tạo điều kiện thuận lợi, giúp đỡ tác giả suốt trình học tập, nghiên cứu hoàn thành luận văn Tác giả xin chân thành cảm ơn nhà khoa học, chuyên gia trường Đại học Dân lập Hải phòng tạo điều kiện giúp đỡ, quan tâm góp ý cho luận văn hoàn thiện Tác giả xin trân trọng cảm ơn cán bộ, giáo viên Khoa xây dựng, Phòng đào tạo Đại học Sau đại học - trường Đại học Dân lập Hải phòng, đồng nghiệp tạo điều kiện thuận lợi, giúp đỡ tác giả trình nghiên cứu hoàn thành luận văn Tác giả luận văn Nguyễn Văn Trƣờng MỤC LỤC MỞ ĐẦU CHƢƠNG 1.BÀI TOÁN CƠ HỌC KẾT CẤU VÀ CÁC PHƢƠNG PHÁP GIẢI 1.1 Bàitoán học kết cấu 1.2 Các phươngpháp giải 1.2.1 Phươngpháp lực 1.2.2 Phươngpháp chuyển vị 1.2.3 Phươngpháp hỗn hợp phươngphápliê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ến phân CHƢƠNG 2: PHƢƠNG PHÁ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ển vị 2.1.1.1 Rời rạc hoá miền khảo sát 2.1.1.2 Chọn hàm xấp xỉ 2.1.1.3 Xây dựng phương trình cân phần tử, thiết lập ma trận độ cứng K e vectơ tảitrọng nút Fe phầntử thứ e 2.1.1.4 Ghép nối phầntử xây dựng phương trình cân toàn hệ .12 2.1.1.5: Sử lý điều kiện biên toán 21 2.1.1.6 Giải hệ phương trình cân 27 2.1.1.7 Xác định nội lực 27 2.1.2 Cách xây dựng ma trận độ cứng phầntửchịu uốn 28 2.1.3 Cách xây dựng ma trận độ cứng tổng thể kết cấu 30 CHƢƠNG PHƢƠNG PHÁPPHẦNTỬHỮUHẠNĐỐIVỚIDẦMCHỊU UỐN 35 3.1 Lý thuyết dầm Euler – Bernoulli [ ] 35 3.1.1 Dầmchịu uốn túy phẳng 35 3.1.2 Dầmchịu uốn ngang phẳng 38 3.2.Giải toándầmliêntụcphươngphápphầntửhữuhạn 44 3.2.1.Tính toándầmliêntục 44 KẾT LUẬN VÀ KIẾN NGHỊ 65 KẾT LUẬN 65 Danh môc tµi liÖu tham kh¶o 66 MỞ ĐẦU Bàitoá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 vi phâ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áp chuyển vị, phươngpháp hỗn hợp, phươngphápliê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ến phânPhươ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êntục Để giải toá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ển vị đại lượng cần tìm hàm nội suy biểu diễn gần dạng phân bố chuyển vị phần tử; Mô hình cân bằng, hàm nội suy biểu diễn gần dạng phân bố ứng suất hay nội lực phầntử mô hình hỗn hợp, coi đại lượng chuyển vị ứ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ạng phân bố chuyển vị lẫn ứng suất phầntửĐối tƣợng, phƣơng pháp phạm vi nghiên cứu đề tàiTrong luận văn này, tác giả sử dụng phươngphầntửhữuhạn theo mô hình chuyển vị để xây dựng giải toándầmliêntụcchịu tác dụng tảitrọngtĩnhtậptrung Mục đích nghiên cứu đề tài “Xác định nội lực chuyển vị dầmliêntụcchịutảitrọngtĩnhtậptrungphươ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 Trình bày phươngphápphầntửhữuhạn áp dụng để giải toándầmliên tục, chịu tác dụng tảitrọngtĩnhtậptrung Lập chương trình máy tính điện tử cho toán nêu CHƢƠNG BÀITOÁN CƠ HỌC KẾT CẤU VÀ CÁC PHƢƠNG PHÁ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àitoán học kết cấu Bàitoán học kết cấu nhằm xác định nội lực chuyển vị hệ thanh, tấm, vỏ tác dụng loại tải trọng, nhiệt độ, chuyển vị cưỡng bức,…và chia làm hai loại: - Bàitoántĩnh định: toá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ảitrọng Để xác định nội lực chuyển vị cần dùng phương trình cân tĩnh học đủ; - Bàitoán siêu tĩnh: toán có 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ển vị cưỡng bức,…Để xác định nội lực chuyển vị phương trình cân ta phải bổ sung phương trình biến dạng Nếu tính đến tận ứng suất, nói toán học vật rắn biến dạng nói chung toán học kết cấu nói riêng toán siêu tĩnh 1.2 Các phƣơng phá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áp lực phươngpháp chuyển vị 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áp chuyển vị 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ƣơng pháp lực Trong hệ siêu tĩnh ta thay liên kết thừa lực chưa biết, giá trị chuyển vị hệ tương ứng với vị trí phương lực ẩn số thân lực nguyên nhân bên 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ƣơng pháp chuyển vị Khác vớiphươngpháp lực, phươngpháp chuyển vị lấy chuyển vị nút làm ẩn Những chuyển vị phải có giá trị cho phản lực liên kết đặt thêm vào hệ thân chúng nguyên nhân bên 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áp chuyển vị 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ƣơng pháp hỗn hợp phƣơng phápliên hợp Phươngpháp hỗn hợp, phươngphápliên hợp kết hợp song song phươngpháp lực phươngpháp chuyển vị Trongphươngpháp ta chọn hệ theo phươngpháp lự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ớiphươngpháp lực; chọn hệ theo phươngpháp chuyển vị không đặt đầy đủ liên kết phụ nhằm ngăn cản toàn chuyển vị nút mà đặt liên kết phụ nút thuộc phận thích hợp vớiphươngpháp chuyển vị 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, toán ban đầu đưa hai toán độc lập: Một theo phươngpháp lực theo phươngpháp chuyển vị 1.2.4 Phƣơng pháp sai phânhữuhạnPhươngpháp sai phânhữuhạn thay hệ liêntụ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 vi phân chuyển vị nội lự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 vi phân chuyển vị nội lực viết dạng sai phân nút, biểu thị quan hệ chuyển vị nút nút lân cận tác dụng ngoại lực 1.2.5 Phƣơng pháp hỗn hợp sai phân – biến phân Kết hợp phươngpháp sai phânvớiphươngpháp biến phân ta có phươngpháp linh động hơn: Hoặc sai phân đạo hàm phương trình biến phân sai phân theo phương biến phân theo phương khác (đối vớitoán hai chiều) Ví dụ 3.3: Dầm hai nhịp (hình 3.5) Xác định nội lực chuyển vị dầmliêntục hai nhịp có tổng chiều dài nhịp l , độ cứng uốn EJ, chịutải SO DO DAMtrọngtậptrung P nhịp 2, hình 3.5a Rời rạc hóa kết cấu dầm thành n pt 2 n 10 ngx nút SO DO NUT DAMphầntử Các nút phầntử phải trùngvớ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 Mỗi phầntử có ẩn 1 0 W SO DO AN CHUYEN VI SO DO AN GOC XOAY n pt phầntử rời rạc tổng cộng có n pt ẩn Hình 3.5 Dầm hai nhịp Nhưng cần đảm bảo liêntục chuyển vị chuyển vị nút cuối phầntử thứ e chuyển vị nút đầu phầntử thứ e 1 nên số bậc tự nhỏ n pt Khi giải ta cần đảm bảo điều kiện liêntục chuyển vị điều kiện liêntục góc xoay xét cách cách đưa vào điều kiện ràng buộc Ví dụ dầm (ví dụ 3.5a) ta chia thành phầntử (hình 3.5b) Khi chia dầm thành phầntử số nút dầm 5, thứ tựtừ trái sang phải [1, 2, 3, 4, 5] (hình 3.5b), số ẩn chuyển vị nw=, thứ tựtừ trái sang phải [1, 2] (hình 3.5c), ẩn chuyển vị hai đầu vị trí gối trung gian dầm không, ẩn góc xoay ngx=8, thứ tựtừ trái sang phải [3, 4, 5, 6, 7, 8, 9, 10] (hình 3.5d) Như vậy, tổng cộng số ẩn 10 ẩn < 4x4=16 ẩn Gọi ma trận n w ma trận chuyển vị có kích thước n w n pt ,2 ma trận có n pt hàng cột chứa ẩn số chuyển vị nút phầntử (hình 3.1) nw (1, :) 0 1; nw (2, :) 1 0 ; nw (3, :) 0 2; nw (4, :) 2 0 nw 0 1 0 2 0 55 Gọi ma trận ngx ma trận chuyển vị có kích thước ngx(npt,2) ma trận có n pt hàng cột chứa ẩn số góc xoay nút phầntử (hình 3.5) ngx (1, :) 3 4; ngx (2, :) 5 6 ; ngx (3, :) 7 8; ngx (4, :) 9 10 ngx 3 10 Sau biết ẩn số thực dầm ta xây dựng độ cứng tổng thể dầm (có nhiều cách ghép nối phầntử khác nhau, tùy vào trình độ lập trình người nên tác giả không trình bày chi tiết cách ghép nối phầntử lại để ma trận độ cứng toàndầm xem code mô đun chương trình tác giả) Nếu toán có nw ẩn số chuyển vị n gx ẩn số góc xoay ma trận độ cứng dầm K có kích thước (nxn), K n,n với n=(nw+ngx) Như ví dụ 3.3, n=10 Bây xét điều kiện liêntục góc xoay phầntử Điều kiện liêntục góc xoay phầntử viết sau: dyi dx nut dyi 1 0 dx nut1 (a) nut1 dy2 dy3 hay: 2 0 dx dx nut nut1 dy dy 3 dx dx nut nut1 dy1 dy dx nut dx 1 (b) Trong i ẩn số toán (có k ẩn số), tổng số ẩn số toá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 K n k,n k Gọi k góc xoay nút phầntử trước, k góc xoay nút phầntử sau ta có hệ số ma trận độ cứng K: k n i,k1 2 (i k) ; k n i,k x x (c) 56 2 (i k) ; k k ,n i (d) x x Nếu có hai phầntử có điều kiện góc xoay, có n pt phầntử k k1 ,n i có 2n pt điều kiện liêntục góc xoay phầntử Như cuối ta thiết lập phương trình: K F F1 so hang n Fn đó: F ; 0 so hang k (e) 1 1 n ẩn số toán 1 2 k Trong ví dụ 3.3 chia thành phần tử, ta có: - Ma trận độ cứng phầntử [Ke], nhƣ sau: [ ] [ ] - Ma trận độ cứng toàndầ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 sau: 57 1536 - 96 - 96 96 96 K 1056 1536 0 0 - 96 - 96 96 96 0 0 - 96 16 0 0 0 0 - 128 - 96 16 0 0 0 0 - 136 96 0 16 0 0 -1 0 0 96 0 16 0 0 0 0 - 96 - 96 0 0 0 0 16 8 16 0 0 0 -1 0 0 0 96 0 0 0 16 0 -1 0 96 0 0 0 16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 - Véc tơ lực nút{ }: F 0 0 0 0 0 0 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: 58 1 0.0000 - 0.0022 W2 - 0.0003 W x Pl ; 0.0063 3 x Pl W4 0.0017 0.0012 4 0.0110 5 Mômen uốn dầm: M1 - 0.0130 M - 0.0297 M M - 0.0557 x Pl M 0.0971 4 M Ta thấy kết trên: - Về mômen gần trùng khớp với kết giải xác theo phươngpháp giải tích: + Tại gối trung gian: x 10 -4 -5 -10 -15 + Tại dầm: -20 10 12 14 Hình 3.6a Đường độ võng dầm - Về chuyển vị nhịp trùng khớp với kết giải xác theo phươngpháp giải tích: Hình 3.6a Biểu đồ M Nhận xét: Nếu ta rời rạc hóa dầm thành 16 phầntử kết trùng khớp với kết xác nhận phươngpháp giải tích Khi biểu đồ mômen uốn đường độ võng dầm hình 3.6: 59 16 Ví dụ 3.4: Dầm hai nhịp (hình 3.7) Xác định nội lực chuyển vị dầm hai nhịp chiều dài nhịp l , độ cứng uốn EJ, chịu mômen tậptrung M nhịp lực tậptrung P nhịp 2, SO DO DAM hình 3.7a 2 n 10 ngx nút SO DO NUT DAM Rời rạc hóa kết cấu dầm thành n pt 1 0 W SO DO AN CHUYEN VI phầntử Các nút phầntử phải trùngvớ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 Mỗi phầntử có ẩn SO DO AN GOC XOAY Hình 3.7 Dầm hai nhịp n pt phầntử rời rạc tổng cộng có n pt ẩn Nhưng cần đảm bảo liêntục chuyển vị chuyển vị nút cuối phầntử thứ e chuyển vị nút đầu phầntử thứ e 1 nên số bậc tự nhỏ n pt Khi giải ta cần đảm bảo điều kiện liêntục chuyển vị điều kiện liêntục góc xoay xét cách cách đưa vào điều kiện ràng buộc Ví dụ dầm (ví dụ 3.1a) ta chia thành phầntử (hình 3.1b) Như vậy, tổng cộng số ẩn 11 ẩn < 4x4=16 ẩn Gọi ma trận n w ma trận chuyển vị có kích thước n w n pt ,2 ma trận có n pt hàng cột chứa ẩn số chuyển vị nút phầntử (hình 3.1) nw (1, :) 0 1; ngx (2, :) 1 2 ; ngx (3, :) 2 3; ngx (4, :) 3 4 nw 0 1 2 3 4 Gọi ma trận n ma trận chuyển vị có kích thước n n pt ,2 ma trận có n pt hàng cột chứa ẩn số góc xoay nút phầntử (hình 3.5) 60 ngx (1, :) 5 6; ngx (2, :) 7 8 ; ngx (3, :) 9 10; ngx (4, :) 11 12 ngx 5 10 11 12 Sau biết ẩn số thực ta xây dựng độ cứng tổng thể (có nhiều cách ghép nối phầntử khác nhau, tùy vào trình độ lập trình người nên tác giả không trình bày chi tiết cách ghép nối phầntử lại để ma trận độ cứng toàndầm xem code mô đun chương trình tác giả) Nếu toán có n cv ẩn số chuyển vị n gx ẩn số góc xoay ma trận độ cứng dầm K có kích thước (nxn), K n,n với n n cv n gx Như ví dụ 3.3, n 11 Bây xét điều kiện liêntục góc xoay phầntử Điều kiện liêntục góc xoay phầntử viết sau: dyi dx nut dyi 1 0 dx nut1 nut1 dy2 dy3 2 0 dx nut dx nut1 dy3 dy4 3 0 dx dx nut nut1 (a) dy1 dy dx nut dx 1 hay: (b) Trong i ẩn số toán (có k ẩn số), tổng số ẩn số toá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 K n k,n k Gọi k góc xoay nút phầntử trước, k góc xoay nút phầntử sau ta có hệ số ma trận độ cứng K: 2 (i k) ; k n i,k x x 2 k k1 ,n i (i k) ; k k ,n i x x k n i,k1 (c) (d) 61 Nếu có hai phầntử có điều kiện góc xoay, có n pt phầntử có 2n pt điều kiện liêntục góc xoay phầntử Như cuối ta thiết lập phương trình: K F F1 so hang n Fn đó: F ; 0 so hang k (e) 1 1 n ẩn số toán 1 2 k Trong ví dụ 3.1 chia thành phần tử, ta có: - Ma trận độ cứng phầntử [Ke], nhƣ sau: [ ] [ ] - Ma trận độ cứng toàndầ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 sau: 1536 0 1536 - 96 96 96 96 - 96 K - 96 96 96 0 0 0 0 1056 Vé - 96 - 96 0 16 8 16 0 0 0 0 0 0 0 0 - 128 - 136 96 0 16 0 0 -1 0 0 96 0 16 0 0 0 0 - 96 0 0 16 0 -1 0 0 - 96 0 0 16 0 0 0 96 0 0 0 16 0 -1 0 96 0 0 0 16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 62 c tơ lực nút{ }: F 1 0 0 0 0 0 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ị mô men uốn chia dầm thành 160 phầntử nhƣ sau: W24 0.0007733 W - 0.001984 W 53 x ql ; W94 0.0004531 W135 - 0.0002352 M1 - 0.2144 Trai M 40 0.4606 Phai - 0.5052 M M 40 x ql M 80 0.1607 M - 0.04464 120 M160 0.0000 Ta thấy kết trên: - Về mômen gối trung gian nhịp thứ trùng khớp với kết giải xác theo phươngpháp 63 giải tích: - Momen ngàm nhịp thứ gần trùng khớp với kết xác - Về chuyển vị kết trùng khớp với kết giải xác theo phươngpháp giải tích: Biểu đồ mômen uốn đƣờng độ võng dầm nhƣ hình 3.8: x 10 -3 0.5 -0.5 -1 -1.5 -2 20 40 60 80 100 120 140 160 Hình 3.8a Đường độ võng dầm Hình 3.8a Biểu đồ M Nhận xét: Qua ví dụ thấy rằng, dầmchịutảitrọng Momen tậptrung để nhận kết xác ta phải chia dầm thành nhiều phầntử 64 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ándầmliêntụcchịu tác dụng tảitrọngtĩnhtậptrung Tác giả rút kết luận sau: Trình bày phươngpháp giải toán học kết cấu Trình bày phươngphápphầntửhữuhạntoán học kết cấu Đã trình bày toándầmchịu uốn theo lý thuyết dầm Euler Bernoulli Bằng phươngphápphầntửhữu hạn, tác giả xác định nội lực chuyển vị dầmliêntụcchịutảitrọngtĩnhtậptrung có điều kiện biên khác Kết nội lực chuyển vị trùng khớp với kết nhận giải phươngpháp có 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 Đốivớitoándầmliêntụcchịutảitrọngtĩnhtậptrung (lực tập trung) để đạt chuyển vị xác cần chia dầm thành từ đến phần tử, để tìm nội lực xác cần chia dầm thành từ đến 16 phầntửĐốivớidầmchịu mômen tậptrung cần chia dầm thành 100 đến 160 phầntử nhận kết xác giải phươngpháp giải tích 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ỏ 65 Danh mục tài liệu tham khảo I TIẾNG VIỆT [1] Hà Huy Cương (2005), Phươngpháp nguyên lý cực trị Gauss, Tạp chí Khoa học kỹ thuật, IV/ Tr 112 118 [2] 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 [3] Phạm Văn Trung (2006), PhươngphápTínhtoán hệ dây mái treo, Luận án Tiến sỹ kỹ thuật [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 Văn Đạo, Trần Kim Chi, Nguyễn Dũng (2005), Nhập môn Động lực học phi tuyến chuyển động hỗn độn Nhà xuất Đại học Quốc gia Hà nội [6] Đoàn Văn Duẩn (2007), Phươngpháp nguyên lý Cực trị Gauss toán ổn định công trình, Luận văn thạc sỹ kỹ thuật [7] Đoàn Văn Duẩn (2010), Phươngphápphầntửhữuhạn nghiên cứu ổn định uốn dọc thanh, Tạp chí kết cấu Công nghệ xây dựng, số 05, Qúy IV(Tr30-Tr36) [8] Đoàn Văn Duẩn (2011), Nghiên cứu ổn định đàn hồi hệ thanh, Luận án Tiến sỹ kỹ thuật [9] Đoàn Văn Duẩn (2012), Phươngpháptínhtoán dây mềm, Tạp chí kết cấu công nghệ Xây dựng số 09, Qúy II (Tr56-Tr61) [10] Đoàn Văn Duẩn (2014), Phươngpháp chuyển vị cưỡng giải toán trị riêng véc tơ riêng, Tạp chí Xây dựng số 11 (Tr82-Tr84) 66 [11] Đoàn Văn Duẩn (2015), Bàitoán học kết cấu dạng tổng quát, Tạp chí Xây dựng số 02 (Tr59-Tr61) [12] Đoàn Văn Duẩn (2015), Phươngpháp so sánh nghiên cứu nội lực chuyển vị hệ dầm, Tạp chí Xây dựng số 11 (Tr56-Tr58) [13] Đoàn Văn Duẩn (2015), Tínhtoá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) [14] Trần Thị Kim Huế (2005), Phươngpháp nguyên lý Cực trị Gauss toán học kết cấu, Luận văn thạc sỹ kỹ thuật [15] Nguyễn Thị Liên (2006), Phươngpháp nguyên lý Cực trị Gauss toán động lực học công trình, Luận văn thạc sỹ kỹ thuật [16] 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 II TIẾNG PHÁP [17] Robert L‟Hermite (1974), Flambage et Stabilité – Le flambage élastique des pièces droites, édition Eyrolles, Paris IIi TIẾNG ANH [18] Stephen P.Timoshenko-Jame M.Gere (1961), Theory of elastic stability, McGraw-Hill Book Company, Inc, New york – Toronto – London, 541 Tr [19] William T.Thomson (1998), Theory of Vibration with Applications (Tái lần thứ 5) Stanley Thornes (Publishers) Ltd, 546 trang [20] Klaus – Jurgen Bathe (1996), Finite Element procedures Part one, Prentice – Hall International, Inc, 484 trang [21] Klaus – Jurgen Bathe (1996), Finite Element procedures Part two, Prentice – Hall International, Inc, 553 trang [22] Ray W.Clough, Joseph Penzien(1993), Dynamics of Structures (Tái lần thứ 2), McGraw-Hill Book Company, Inc, 738 trang 67 [23] O.C Zienkiewicz-R.L Taylor (1991), The finite element method (four edition) Volume 2, McGraw-Hill Book Company, Inc, 807 trang [24] 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) [25] 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 [26] D.R.J Owen, E.Hinton (1986), Finite Elements in Plasticity: Theory and Practice, Pineridge Press Lt [27] 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 [28] 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) [29] Chopra Anil K (1995) Dynamics of structures Prentice Hall, Englewood Cliffs, New – Jersey 07632 [30] Wilson Edward L Professor Emeritus of structural Engineering University of California at Berkeley (2002) Three – Dimensional Static and Dynamic Analysis of structures, Inc Berkeley, California, USA Third edition, Reprint January [31] 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 68 [32] 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 [33] 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 [34] Kolousek Vladimir, DSC Professor, Technical University, Pargue (1973) Dynamics in engineering structutes Butter worths London [35] 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 69 ... Các phương pháp giải gồm có: Phương pháp coi xác như, 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. .. PHƢƠNG PHÁP PHẦN TỬ HỮU HẠN Trong chương trình bày số khái niệm phương pháp phần tử hữu hạn, để phục vụ cho việc xây dựng toán xác định nội lực chuyển vị cho dầm liên tục chịu tải trọng tĩnh tập trung. .. phương pháp phần tử hữu hạn chương 2.1 Phƣơng pháp phần tử hữu hạn Phương pháp phần tử hữu hạn phương pháp số đặc biệt có hiệu để tìm dạng gần hàm chưa biết miền xác định V Tuy nhiên phương pháp