BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ CÔNG NGHIỆP LONG AN ------------------------------------------------------------------------ PHẠM THỊ ANH PHƯƠNG PHÂN TÍCH ỨNG XỬ BẤT ỔN ĐỊNH CỦA TẤM TĂNG CỨNG CHỊU TẢI NÉN VÀ ÁP LỰC NGANG LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây Dựng Long An, năm 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ CÔNG NGHIỆP LONG AN PHẠM THỊ ANH PHƯƠNG PHÂN TÍCH ỨNG XỬ BẤT ỔN ĐỊNH CỦA TẤM TĂNG CỨNG CHỊU TẢI NÉN VÀ ÁP LỰC NGANG LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây Dựng Mã ngành : 8.01 Người hướng dẫn khoa học: PGS. TRƯƠNG TÍCH THIỆN Long An, năm 2020 i LỜI CAM ĐOAN Tôi xin cam đoan luận văn này là công trình nghiên cứu của riêng tôi. Các số liệu, và kết quả trong luận văn là trung thực và chưa được công bố trong các tạp chí khoa học và công trình nào khác. Các thông tin số liệu trong luận văn này đều có nguồn gốc và được ghi chú rõ ràng. Tác giả (Ký và ghi rõ họ tên) ii LỜI CẢM ƠN Xin cám ơn Thầy PGS. Trương Tích Thiện, Thầy đã đưa ra gợi ý đầu tiên để hình thành nên ý tưởng của đề tài và đã làm tôi mạnh dạn tiếp cận với hướng nghiên cứu đồng thời, thầy là người đã tận tụy giúp tôi hệ thống hóa lại kiến thức quản lý và định lượng phân tích và hiểu biết thêm về nhiều điều mới trong quá trình nghiên cứu luận văn này. Để hoàn thành đề cương luận văn này, ngoài sự cố gắng và nỗ lực của bản thân, tôi đã nhận được nhiều sự giúp đỡ của tập thể và các cá nhân. Tôi xin tỏ lòng biết ơn đến tập thể và các cá nhân đã dành cho tôi sự giúp đỡ quý báu đó. Tôi xin chân thành cảm ơn Quý Thầy Cô Khoa Xây dựng trường Đại Học Kinh tế Công nghiệp Long An đã tạo điều kiện và giúp đỡ tôi trong quá trình học và nghiên cứu khoa học tại đây. Luận văn thạc sĩ đã hoàn thành trong thời gian quy định với sự nỗ lực của bản thân, tuy nhiên không thể không có những thiếu sót. Kính mong Quý Thầy Cô chỉ dẫn thêm để tôi bổ sung những kiến thức và hoàn thiện bản thân mình hơn. Xin trân trọng cảm ơn! Tác giả (Ký và ghi rõ họ tên) iii TÓM TẮT Tấm mỏng là dạng kết cấu dễ bị mất ổn định khi chịu tác động của lực nén và áp lực ngang. Trong luận văn này, phương pháp phần tử hữu hạn với nhiều ưu điểm trong phân tích kết cấu đã được chứng minh trong các nghiên cứu trước đó được sử dụng để phân tích hậu bất ổn định cho kết cấu tấm vật liệu đẳng hướng và tấm vật liệu composite. Các kết quả được tính toán thông qua chương trình ANSYS và được kiểm chứng so với kết quả giải tích. Luận văn đã thực hiện phân tích ứng xử của 2 mô hình tấm: tấm đẳng hướng và tấm Composite. Với mỗi mô hình tấm: - Luận văn phân tích tải tới hạn khi tấm chịu tải nén đơn trục, kết quả phân tích từ ANSYS được so sánh với kết quả giải tích. - Luận văn phân tích ảnh hưởng của tải trọng ngang đến khả năng ổn định của tấm. Kết quả cho thấy, khi có tải trọng ngang tác động thì tải tới hạn của tấm khi chịu nén đơn trục sẽ tăng lên. iv ASBTRACT Thin plates are susceptible to buckling under compression and lateral pressure. In this thesis, the finite element method (FEM) with many advantages in structural analysis has been demonstrated in previous studies to be used to analyze postbuckling behaviour of the isotropic plates and composite plates. The result is analyzed by ANSYS and compared to the theoretical results. Two plate models: isotropic and composite plate are analyzed in thesis. For each plate model: - The critical load of the isotropic plate which is uniaxial compression is analyzed, the results from ANSYS are compared with the analytical results. - The effect of horizontal load on the stability of the plate also being studied. The results show that, the critical load of the plate under uniaxial compression will increase. v MỤC LỤC LỜI CAM ĐOAN.iv MỤC LỤC . v DANH MỤC TỪ VIẾT TẮT TIẾNG VIỆT .vii DANH MỤC VIẾT TẮT TIẾNG ANH . viii DANH MỤC BẢNG BIỂU .ix DANH MỤC ĐỒ THỊ VÀ HÌNH VẼ .iv CHƯƠNG 1. TỔNG QUAN VỀ ĐỀ TÀI . Tính cần thiết của đề tài .1 Tình hình nghiên cứu trong nước . Tình hình nghiên cứu ở nước ngoài . Mục tiêu nghiên cứu . Đối tượng và phạm vi nghiên cứu . Câu hỏi nghiên cứu và giả thuyết nghiên cứu . Phương pháp nghiên cứu . Các khái niệm cơ bản . Trạng thái ổn định kết cấu . Trạng thái bất ổn định kết cấu . Lý thuyết tấm mỏng . Tấm đẳng hướng . Lý thuyết bất ổn định tấm mỏng chịu tải nén đơn trục . Tấm đẳng hướng . Lý thuyết bất ổn định tấm mỏng chịu tải nén đơn trục và lực theo phương ngang . Ứng xử biến dạng lớn của tấm . Phân tích bất ổn định của tấm chịu tải kết hợp. Phân tích sau bất ổn định của tấm chịu tải kết hợp .PHÂN TÍCH BẤT ỔN ĐỊNH VÀ HẬU BẤT ỔN ĐỊNH CỦA TẤM ĐẲNG HƯỚNG . Phân tích bất ổn định tấm chịu nén đơn trục .1 Thông số vật liệu. Điều kiện biên – tải trọng . Kết quả tính toán . Phân tích hậu bất ổn định tấm chịu tải kết hợp . Phân tích tuyến tính .2 Phân tích phi tuyến (biến dạng lớn) .BẤT ỔN ĐỊNH VÀ HẬU BẤT ỔN ĐỊNH CỦA TẤM VẬT LIỆU COMPOSITE . Phân tích bất ổn định tấm composite chịu nén đơn trục . Phân tích hậu bất ổn định tấm composite chịu tải kết hợp. Phân tích tuyến tính (bài toán biến dạng bé) . Phân tích phi tuyến (biến dạng lớn) . KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN . Những ưu điểm chính của luận văn . Những thiếu sót chính của luận văn. Hướng phát triển . 47 DANH MỤC TÀI LIỆU THAM KHẢO` . 48 vii DANH MỤC TỪ VIẾT TẮT TIẾNG VIỆT STT TỪ VIẾT TẮT VIẾT ĐẦY ĐỦ 1 PP PTHH Phương pháp phần tử hữu hạn 2 3 4 5 viii DANH MỤC VIẾT TẮT TIẾNG ANH STT TỪ VIẾT TẮT VIẾT ĐẦY ĐỦ 1 FEM Finite element method 2 ANSYS Analysis System ix DANH MỤC BẢNG BIỂU Bảng Trang Bảng 3. Thông số vật liệu tấm đẳng hướng . Kết quả tính toán tải bất ổn định . Giá trị lực tới hạn của tấm đẳng hướng dưới các tải phân bố . Giá trị tải phân bố thay đổi theo thời gian tác dụng lên tấm đẳng hướng .5: Giá trị tải tới hạn của tấm đẳng hướn trong trường hợp biến dạng lớn . Thông số vật liệu . Các thông số về lớp composite . Lực tới hạn của tấm composite chịu nén đơn trục . Kết quả tải bất ổn định tương ứng với các mức tải pz . Giá trị lực phân bố thay đổi theo thời gian tác dụng lên tấm đẳng hướng .6: Giá trị lực tới hạn của tấm composite trong trường hợp biến dạng lớn . 43 iv DANH MỤC ĐỒ THỊ VÀ HÌNH VẼ Hình Trang Hình 1. Mô hình tấm chịu nén và bị bất ổn định . Tấm hình chữ nhật chịu tải nén song trục và tải trọng ngang . Các thành phần nội lực trong tấm. Mô hình tấm hình chữ nhật chịu nén đơn trục . Mối quan hệ giữa hệ số bất ổn định và tỉ lệ kích thước tấm . Mô hình tấm composite . Điều kiện biên mô hình tấm khảo sát . Mô hình tấm khảo sát. Điều kiện biên bài toán . Mô hình bài toán tấm chịu nén đơn trục . Đồ thị đánh giá mức độ hội tụ lưới của tấm đẳng hướng . Kết quả 4 dạng biến dạng tương ứng với 4 tải bất ổn định đầu tiên . Dạng bất ổn định thứ nhất tương ứng với các mức tải . Đồ thị mô tả đặt tải phân bố vuông góc theo thời gian . Trường chuyển vị của tấm trước khi bất ổn định . Trường chuyển vị của tấm lúc xảy ra bất ổn định . Trường chuyển vị của tấm sau khi bất ổn định . Đồ thị tải nén- chuyển vị dọc trục với tải phân bố p0 của tấm đẳng hướng 32 Hình 3. Trường chuyển vị của tấm trước khi bất ổn định . Trường chuyển vị của tấm lúc xảy ra bất ổn định . Trường chuyển vị của tấm sau khi bất ổn định . Đồ thị tải nén- chuyển vị với tải phân bố p1 của tấm đẳng hướng . Trường chuyển vị của tấm trước khi bất ổn định . Trường chuyển vị của tấm lúc xảy ra bất ổn định . Trường chuyển vị của tấm sau khi bất ổn định .17: Đồ thị tải nén chuyển vị dọc trục với tải p2 của tấm đẳng hướng . Trường chuyển vị của tấm trước khi bất ổn định . Trường chuyển vị của tấm lúc xảy ra bất ổn định . Trường chuyển vị của tấm sau khi bất ổn định . Đồ thị tải nén - chuyển vị dọc trục với tải P3 của tấm đẳng hướng . Các mode bất ổn định . Mô hình bài toán . Trường chuyển vị của tấm tương ứng với các mức tải pz. Dạng biến dạng trước và sau bất ổn định tương ứng với mức tải p0 . Dạng biến dạng trước và sau bất ổn định tương ứng với mức tải p1 . Dạng biến dạng trước và sau bất ổn định tương ứng với mức tải p2 . Dạng biến dạng trước và sau bất ổn định tương ứng với mức tải p3 . TỔNG QUAN VỀ ĐỀ TÀI 1. Tính cần thiết của đề tài Trong xây dựng hiện đại ngày nay, các kết cấu dạng tấm hiện đang là một xu hướng phát triển tất yếu trong công trình kiến trúc. Việc sử dụng kết cấu dạng tấm đã tạo ra nhiều lựa chọn cho không gian kiến trúc hiện đại cho con người. Đặc biệt là kết cấu dạng tấm được tăng cứng (Stiffened Panels) với nhiều ưu điểm như khả năng chịu lực lớn (so với khối lượng tấm) nên được sử dụng rất phổ biến trong thực tiễn các ngành công nghiệp và dân dụng. Ở các nước phát triển như Mỹ, Nhật, Hàn Quốc… khi xây dựng nhà cao tầng, đặc biệt đối với nhà cao tầng với hệ khung bằng thép, người ta thường sử dụng hệ sàn liên hợp thép – bê tông có sử dụng tấm tôn hình dập nguội (sàn deck).
Tổng quan nghiên cứu
Trong lĩnh vực kỹ thuật xây dựng, kết cấu dạng tấm mỏng tăng cứng ngày càng được ứng dụng rộng rãi trong các công trình hiện đại nhờ khả năng chịu lực lớn so với khối lượng vật liệu. Theo báo cáo ngành, việc sử dụng tấm tăng cứng giúp giảm trọng lượng kết cấu, từ đó giảm tải trọng lên móng và đẩy nhanh tiến độ thi công. Tuy nhiên, do đặc tính mỏng và chịu tải nén cùng áp lực ngang, các tấm này dễ bị mất ổn định, gây nguy cơ sụp đổ đột ngột. Nghiên cứu này tập trung phân tích ứng xử bất ổn định và hậu bất ổn định của tấm tăng cứng chịu tải nén đơn trục và áp lực ngang, nhằm nâng cao độ chính xác trong thiết kế và đảm bảo an toàn kết cấu.
Mục tiêu cụ thể của luận văn là áp dụng phương pháp phần tử hữu hạn (FEM) để phân tích tải tới hạn và ảnh hưởng của tải trọng ngang đến khả năng ổn định của hai loại tấm: tấm đẳng hướng và tấm composite. Phạm vi nghiên cứu tập trung vào mô hình tấm mỏng trong điều kiện tải trọng nén đơn trục kết hợp với tải phân bố vuông góc, sử dụng phần mềm ANSYS để mô phỏng và so sánh với kết quả giải tích. Thời gian nghiên cứu được thực hiện trong năm 2020 tại Trường Đại học Kinh tế Công nghiệp Long An.
Nghiên cứu có ý nghĩa quan trọng trong việc cung cấp cơ sở khoa học cho việc thiết kế kết cấu tấm tăng cứng, góp phần nâng cao độ an toàn và hiệu quả kinh tế trong xây dựng công trình dân dụng và công nghiệp.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Luận văn dựa trên hai lý thuyết chính: lý thuyết tấm mỏng đẳng hướng và lý thuyết tấm composite. Lý thuyết tấm mỏng đẳng hướng sử dụng phương trình von Karman để mô tả trạng thái cân bằng và bất ổn định của tấm dưới tải nén đơn trục và tải phân bố vuông góc. Các phương trình cân bằng, tương thích và quan hệ ứng suất - biến dạng được xây dựng dựa trên giả thiết biến dạng nhỏ và biến dạng lớn, cho phép phân tích ứng xử phi tuyến của tấm.
Đối với tấm composite, luận văn áp dụng lý thuyết tấm composite cổ điển, trong đó vật liệu composite được mô tả là sự kết hợp của pha nền và pha cốt gia cường. Các ma trận độ cứng màng (A), uốn màng (B) và uốn (D) được sử dụng để xác định nội lực và moment uốn trong tấm. Phương trình cân bằng và tương thích của tấm composite tương tự tấm đẳng hướng nhưng có tính chất vật liệu phức tạp hơn do tính dị hướng.
Ba khái niệm chính được sử dụng trong nghiên cứu gồm: trạng thái ổn định và bất ổn định kết cấu, lý thuyết biến dạng lớn của tấm mỏng, và phân tích hậu bất ổn định (postbuckling) nhằm đánh giá ứng xử sau khi tấm bắt đầu mất ổn định.
Phương pháp nghiên cứu
Nguồn dữ liệu chính được thu thập từ các tài liệu khoa học trong và ngoài nước, các bài báo chuyên ngành về phân tích kết cấu tấm và phương pháp phần tử hữu hạn. Luận văn sử dụng phần mềm ANSYS với phần tử Shell181 để mô hình hóa và phân tích ứng xử của tấm đẳng hướng và tấm composite.
Cỡ mẫu mô hình bao gồm tấm hình chữ nhật với kích thước tiêu chuẩn, được chia lưới với mật độ từ 6×24 đến 12×48 phần tử để đánh giá độ hội tụ và sai số. Phương pháp chọn mẫu là chia lưới tinh dần nhằm đảm bảo độ chính xác và hiệu quả tính toán.
Phân tích được thực hiện theo hai giai đoạn: phân tích tuyến tính (biến dạng nhỏ) và phân tích phi tuyến (biến dạng lớn) sử dụng phương pháp Newton-Raphson trong ANSYS. Timeline nghiên cứu kéo dài trong năm 2020, bao gồm giai đoạn thu thập tài liệu, xây dựng mô hình, phân tích số liệu và đối chiếu kết quả với lý thuyết.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
-
Tải tới hạn của tấm đẳng hướng chịu nén đơn trục: Kết quả tính toán từ ANSYS với mật độ lưới 12×48 cho tải tới hạn là khoảng 17.15 MPa, sai số so với kết quả giải tích chỉ khoảng 0.4%, cho thấy độ chính xác cao của phương pháp phần tử hữu hạn. Sai số tăng lên đến 17.22% khi mật độ lưới giảm xuống 6×24.
-
Ảnh hưởng của tải trọng ngang đến tải tới hạn: Trong phân tích biến dạng nhỏ, tải trọng phân bố vuông góc không ảnh hưởng đến tải tới hạn của tấm đẳng hướng. Tuy nhiên, khi phân tích biến dạng lớn, tải trọng ngang làm tăng tải tới hạn lên khoảng 16.16 MPa, phù hợp với thực tế và cho thấy tầm quan trọng của phân tích phi tuyến.
-
Ứng xử hậu bất ổn định của tấm composite: Mô hình tấm composite cho thấy tải tới hạn chịu nén đơn trục thấp hơn so với tấm đẳng hướng do tính dị hướng của vật liệu. Các dạng biến dạng bất ổn định đầu tiên được mô phỏng cho thấy sự khác biệt rõ rệt giữa hai loại tấm, với tấm composite có xu hướng biến dạng phức tạp hơn.
-
Độ hội tụ của mô hình phần tử hữu hạn: Năng lượng biến dạng của tấm đẳng hướng hội tụ ổn định khi mật độ lưới đạt 12×48 phần tử, chứng tỏ không cần chia lưới quá mịn để đạt kết quả chính xác.
Thảo luận kết quả
Nguyên nhân tải trọng ngang làm tăng tải tới hạn trong phân tích biến dạng lớn là do sự tương tác giữa tải nén và tải phân bố vuông góc tạo ra hiệu ứng tăng cứng cho tấm. Kết quả này khác biệt với phân tích biến dạng nhỏ, nơi giả thiết tuyến tính làm giảm khả năng mô phỏng chính xác ứng xử thực tế.
So sánh với các nghiên cứu trước đây, kết quả phù hợp với báo cáo của ngành về việc sử dụng FEM trong phân tích kết cấu tấm composite và đẳng hướng. Việc sử dụng phần tử Shell181 trong ANSYS được chứng minh là phù hợp và hiệu quả cho bài toán này.
Dữ liệu có thể được trình bày qua biểu đồ tải nén - chuyển vị dọc trục, biểu đồ sai số theo mật độ lưới và hình ảnh mô phỏng các dạng biến dạng bất ổn định, giúp trực quan hóa sự khác biệt giữa các mô hình và điều kiện tải.
Đề xuất và khuyến nghị
-
Áp dụng phân tích phi tuyến trong thiết kế kết cấu tấm: Khuyến nghị các kỹ sư sử dụng phương pháp phân tích biến dạng lớn để đánh giá chính xác tải tới hạn và ứng xử hậu bất ổn định, nhằm tránh đánh giá sai khả năng chịu tải của kết cấu.
-
Tối ưu mật độ lưới trong mô hình FEM: Đề xuất sử dụng mật độ lưới khoảng 12×48 phần tử cho mô hình tấm mỏng để cân bằng giữa độ chính xác và thời gian tính toán, áp dụng trong các phần mềm phân tích kết cấu.
-
Phát triển mô hình tấm composite đa lớp: Khuyến khích nghiên cứu mở rộng phân tích cho các tấm composite đa lớp với các hướng sợi khác nhau để nâng cao độ tin cậy trong thiết kế vật liệu mới.
-
Đào tạo và nâng cao năng lực sử dụng phần mềm ANSYS: Đề xuất tổ chức các khóa đào tạo chuyên sâu về phân tích phần tử hữu hạn cho kỹ sư xây dựng nhằm nâng cao chất lượng thiết kế và kiểm tra kết cấu.
Các giải pháp trên nên được thực hiện trong vòng 1-2 năm tới bởi các đơn vị thiết kế, trường đại học và viện nghiên cứu chuyên ngành kỹ thuật xây dựng.
Đối tượng nên tham khảo luận văn
-
Kỹ sư thiết kế kết cấu: Nghiên cứu cung cấp cơ sở khoa học và phương pháp phân tích chính xác để thiết kế tấm tăng cứng chịu tải nén và áp lực ngang, giúp tối ưu hóa kết cấu và đảm bảo an toàn.
-
Giảng viên và sinh viên ngành kỹ thuật xây dựng: Luận văn là tài liệu tham khảo quý giá cho việc giảng dạy và học tập về lý thuyết tấm mỏng, phương pháp phần tử hữu hạn và phân tích hậu bất ổn định.
-
Nhà nghiên cứu vật liệu composite: Cung cấp kiến thức về ứng xử cơ học của tấm composite dưới tải trọng phức tạp, hỗ trợ phát triển vật liệu mới và ứng dụng trong công nghiệp.
-
Chuyên gia kiểm định và giám sát công trình: Giúp đánh giá chính xác khả năng chịu tải và ổn định của kết cấu tấm trong quá trình thi công và vận hành, từ đó đưa ra các biện pháp phòng ngừa sự cố.
Câu hỏi thường gặp
-
Phương pháp phần tử hữu hạn có phù hợp để phân tích ứng xử hậu bất ổn định của tấm không?
Có, phương pháp phần tử hữu hạn, đặc biệt sử dụng phần tử Shell181 trong ANSYS, cho kết quả chính xác và phù hợp với lý thuyết giải tích, giúp mô phỏng được cả biến dạng nhỏ và lớn. -
Tải trọng ngang ảnh hưởng như thế nào đến khả năng ổn định của tấm?
Tải trọng ngang làm tăng tải tới hạn của tấm khi phân tích biến dạng lớn, giúp tấm chịu được tải nén cao hơn so với trường hợp không có tải ngang. -
Mật độ lưới ảnh hưởng ra sao đến kết quả phân tích?
Mật độ lưới càng cao thì kết quả càng chính xác, tuy nhiên sau mức 12×48 phần tử, năng lượng biến dạng và tải tới hạn hội tụ ổn định, không cần chia lưới quá mịn. -
Tấm composite có ưu điểm gì so với tấm đẳng hướng trong ứng xử bất ổn định?
Tấm composite có khối lượng riêng nhỏ, độ bền cao và khả năng chống ăn mòn tốt, nhưng tải tới hạn thường thấp hơn do tính dị hướng và phức tạp trong ứng xử biến dạng. -
Làm thế nào để áp dụng kết quả nghiên cứu vào thiết kế thực tế?
Kết quả phân tích có thể được sử dụng để xác định tải tới hạn và các dạng biến dạng bất ổn định, từ đó thiết kế các tấm tăng cứng phù hợp với điều kiện tải trọng thực tế, đảm bảo an toàn và hiệu quả kinh tế.
Kết luận
- Luận văn đã thành công trong việc áp dụng phương pháp phần tử hữu hạn để phân tích ứng xử bất ổn định và hậu bất ổn định của tấm tăng cứng chịu tải nén và áp lực ngang.
- Kết quả mô phỏng trên ANSYS có độ chính xác cao, sai số dưới 1% so với lý thuyết giải tích khi sử dụng mật độ lưới phù hợp.
- Tải trọng ngang có ảnh hưởng tích cực đến tải tới hạn của tấm trong phân tích biến dạng lớn, điều này cần được tính đến trong thiết kế thực tế.
- Phân tích tấm composite cho thấy tính dị hướng ảnh hưởng đến ứng xử bất ổn định, đòi hỏi phương pháp phân tích chuyên sâu hơn.
- Đề xuất các giải pháp ứng dụng và phát triển nghiên cứu tiếp theo nhằm nâng cao độ tin cậy và hiệu quả trong thiết kế kết cấu tấm tăng cứng.
Tiếp theo, nghiên cứu sẽ mở rộng phân tích cho các loại tấm composite đa lớp và điều kiện tải trọng phức tạp hơn. Các kỹ sư và nhà nghiên cứu được khuyến khích áp dụng phương pháp và kết quả này để nâng cao chất lượng thiết kế kết cấu trong thực tế.