Luận Văn Thạc Sĩ: Phân Tích Giới Hạn Tấm Dày 5 Bậc Tự Do Sử Dụng Phần Tử NS-DSG3

Luận văn thạc sĩ phân tích giới hạn tấm dày 5 bậc tự do với phần tử ns dsg3, cung cấp cái nhìn sâu sắc về ứng dụng trong kỹ thuật.

• Tác giả: Nguyễn Hoài Bảo
• Người hướng dẫn: TS. Nguyễn Ngọc Phúc
• Trường học: Trường Đại Học Kinh Tế Công Nghiệp Long An
• Chuyên ngành: Kỹ Thuật Xây Dựng
• Đề tài: Phân Tích Giới Hạn Tấm Dày 5 Bậc Tự Do Bằng Phương Pháp NS-DSG3
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2019
• Địa điểm: Long An

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CÁM ƠN
• TÓM TẮT LUẬN VĂN
• ABSTRACT
• DANH MỤC CHỮ VIẾT TẮT
• DANH MỤC HÌNH ẢNH
• DANH MỤC BẢNG
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.1. Giới thiệu chung
• 1.2. Lý do chọn đề tài
• 1.3. Ý nghĩa của đề tài
• 1.3.1. Ý nghĩa khoa học
• 1.3.2. Ý nghĩa thực tiễn
• 1.4. Mục tiêu, đối tượng, phạm vi và giới hạn nghiên cứu
• 1.4.1. Mục tiêu và giới hạn nghiên cứu
• 1.4.2. Đối tượng và phạm vi nghiên cứu
• 1.5. Tóm tắt luận điểm cơ bản và đóng góp mới
• 1.6. Cấu trúc luận văn
• 2. CHƯƠNG 2: TỔNG QUAN
• 2.1. Tổng quan tài liệu trong nước
• 2.2. Tổng quan tài liệu nước ngoài
• 2.3. Những vấn đề còn tồn tại và mục tiêu, nhiệm vụ luận văn
• 3. CHƯƠNG 3: CƠ SỞ LÝ THUYẾT
• 3.1. Mô hình vật liệu
• 3.1.1. Định đề ổn định của Drucker
• 3.1.2. Luật chảy dẻo kết hợp
• 3.2. Lý thuyết tấm dày
• 3.3. Phương pháp phân tích giới hạn
• 3.3.1. Các giả thuyết cơ sở của phương pháp phân tích giới hạn
• 3.3.2. Các định lý cơ bản của phương pháp phân tích giới hạn
• 3.4. Phương pháp Node Smooth FEM-DSG3 (NS-DSG3)
• 3.4.1. Phần tử Discrete Shear Gap (DSG3)
• 3.4.2. Phần tử NS-DSG3
• 3.5. Bài toán phân tích giới hạn cho tấm dày (5 bậc tự do)
• 4. CHƯƠNG 4: KẾT QUẢ VÀ BÀN LUẬN
• 4.1. Tấm hình vuông
• 4.1.1. Tấm hình vuông bốn biên tựa
• 4.1.2. Tấm hình vuông bốn biên ngàm
• 4.2. Tấm hình chữ nhật
• 4.2.1. Tấm hình chữ nhật bốn biên ngàm (b=2a)
• 4.2.2. Tấm hình chữ nhật bốn biên tựa (b=2a)
• 4.3. Tấm hình tròn
• 4.3.1. Tấm hình tròn biên tựa chu vi
• 4.3.2. Tấm hình tròn biên ngàm theo chu vi
• 4.4. Tấm hình chữ L
• 4.4.1. Tấm chữ L hai biên ngàm
• 4.4.2. Tấm chữ L hai biên tựa
• 5. CHƯƠNG 5: KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về phân tích giới hạn tấm dày 5 bậc tự do

Phân tích giới hạn tấm dày 5 bậc tự do là một lĩnh vực quan trọng trong kỹ thuật xây dựng. Phương pháp này giúp xác định tải trọng phá hoại của kết cấu tấm dày, từ đó dự đoán cơ cấu phá hoại. Việc áp dụng lý thuyết tấm dày Mindlin trong phân tích giới hạn mang lại nhiều lợi ích, đặc biệt là trong việc xử lý các vật liệu bất đối xứng.

1.1. Khái niệm về tấm dày và bậc tự do

Tấm dày là loại kết cấu có chiều dày lớn so với các kích thước khác. Bậc tự do trong phân tích tấm dày đề cập đến số lượng chuyển vị độc lập có thể xảy ra trong kết cấu. Việc hiểu rõ khái niệm này là cần thiết để áp dụng phương pháp NS-DSG3.

1.2. Lý thuyết tấm dày Mindlin

Lý thuyết tấm dày Mindlin cho phép tính toán chính xác hơn cho các tấm có chiều dày lớn. Nó xem xét các biến dạng màng và biến dạng cắt, giúp cải thiện độ chính xác trong phân tích giới hạn.

II. Vấn đề và thách thức trong phân tích giới hạn tấm dày

Phân tích giới hạn tấm dày gặp nhiều thách thức, đặc biệt là trong việc xác định tải trọng phá hoại. Các vấn đề như hiện tượng shear-locking và sự phức tạp trong mô hình hóa vật liệu bất đối xứng cần được giải quyết. Những thách thức này ảnh hưởng đến độ chính xác và hiệu quả của phương pháp phân tích.

2.1. Hiện tượng shear locking trong tấm dày

Shear-locking là hiện tượng xảy ra khi tấm dày có chiều dày mỏng dần, dẫn đến sai số trong tính toán. Việc áp dụng phương pháp DSG giúp giảm thiểu hiện tượng này, từ đó cải thiện độ chính xác của kết quả.

2.2. Mô hình hóa vật liệu bất đối xứng

Mô hình hóa vật liệu bất đối xứng là một thách thức lớn trong phân tích giới hạn. Việc áp dụng lý thuyết tấm dày 5 bậc tự do giúp giải quyết vấn đề này, cho phép tính toán chính xác hơn cho các vật liệu có cấu trúc phức tạp.

III. Phương pháp NS DSG3 trong phân tích giới hạn tấm dày

Phương pháp NS-DSG3 là một kỹ thuật tiên tiến trong phân tích giới hạn tấm dày. Phương pháp này kết hợp giữa phần tử NS-FEM và kỹ thuật DSG, giúp cải thiện độ chính xác và giảm chi phí tính toán. Việc áp dụng phương pháp này vào bài toán tấm dày 5 bậc tự do mang lại nhiều kết quả khả quan.

3.1. Nguyên lý hoạt động của phương pháp NS DSG3

Phương pháp NS-DSG3 hoạt động dựa trên việc lấy trung bình biến dạng quanh điểm nút, giúp giảm số lượng biến trong bài toán tối ưu. Điều này không chỉ tiết kiệm thời gian tính toán mà còn nâng cao độ chính xác của kết quả.

3.2. Ưu điểm của phương pháp NS DSG3

Phương pháp NS-DSG3 có nhiều ưu điểm như khả năng tránh hiện tượng shear-locking, tính chính xác cao trong tính toán và khả năng áp dụng cho nhiều hình dạng tấm khác nhau. Những ưu điểm này làm cho phương pháp trở thành lựa chọn hàng đầu trong phân tích giới hạn.

IV. Ứng dụng thực tiễn của phân tích giới hạn tấm dày

Phân tích giới hạn tấm dày 5 bậc tự do có nhiều ứng dụng thực tiễn trong ngành xây dựng. Việc xác định tải trọng phá hoại và cơ cấu phá hoại của kết cấu tấm giúp các kỹ sư thiết kế các công trình an toàn và hiệu quả hơn. Các kết quả nghiên cứu từ phương pháp NS-DSG3 đã được áp dụng vào nhiều dự án thực tế.

4.1. Ứng dụng trong thiết kế công trình

Phân tích giới hạn tấm dày giúp các kỹ sư dự đoán chính xác tải trọng phá hoại, từ đó thiết kế các công trình có khả năng chịu lực tốt hơn. Điều này đặc biệt quan trọng trong các công trình lớn như cầu, nhà cao tầng.

4.2. Kết quả nghiên cứu và so sánh

Các kết quả từ phương pháp NS-DSG3 đã được so sánh với các nghiên cứu trước đó, cho thấy độ chính xác và hiệu quả của phương pháp. Những so sánh này cung cấp cơ sở vững chắc cho việc áp dụng phương pháp trong thực tế.

V. Kết luận và tương lai của phân tích giới hạn tấm dày

Phân tích giới hạn tấm dày 5 bậc tự do sử dụng phương pháp NS-DSG3 đã chứng minh được tính hiệu quả và độ chính xác cao. Tương lai của nghiên cứu này hứa hẹn sẽ mở ra nhiều hướng đi mới trong việc phát triển các phương pháp phân tích giới hạn cho các loại kết cấu phức tạp hơn.

5.1. Tóm tắt kết quả nghiên cứu

Nghiên cứu đã chỉ ra rằng phương pháp NS-DSG3 có thể áp dụng hiệu quả cho nhiều hình dạng tấm khác nhau. Các kết quả đạt được cho thấy tính khả thi và độ chính xác của phương pháp trong phân tích giới hạn.

5.2. Hướng phát triển trong tương lai

Tương lai của phân tích giới hạn tấm dày sẽ tập trung vào việc phát triển các phương pháp mới, cải thiện độ chính xác và hiệu quả tính toán. Việc nghiên cứu thêm về các vật liệu mới và hình dạng phức tạp sẽ là những hướng đi quan trọng.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng và cơ học kết cấu, việc xác định chính xác tải trọng phá hoại đóng vai trò then chốt nhằm tối ưu hóa vật liệu và bảo đảm an toàn công trình. Các cấu kiện dạng tấm dày chịu tải trọng uốn và lực cắt lớn thường xuất hiện phổ biến tại các sàn đáy hầm, móng bè hoặc sàn chuyển nhà cao tầng. Tuy nhiên, việc áp dụng lý thuyết tấm mỏng cổ điển đối với các kết cấu có chiều dày lớn thường dẫn đến sai số nghiêm trọng do bỏ qua biến dạng cắt. Đồng thời, các phương pháp phần tử hữu hạn thông thường khi tính toán tấm mỏng dần lại dễ gặp hiện tượng khóa cắt (shear locking) và tiêu tốn tài nguyên tính toán lớn khi số lượng điểm tích phân Gauss tăng cao.

Nghiên cứu của học viên Nguyễn Hoài Bảo dưới sự hướng dẫn khoa học của Tiến sĩ Nguyễn Ngọc Phúc tại Trường Đại học Kinh tế Công nghiệp Long An đã giải quyết triệt để bài toán phân tích giới hạn cận trên cho kết cấu tấm dày 5 bậc tự do. Mục tiêu cốt lõi của đề tài là xây dựng mô hình toán học tối ưu hóa hình nón bậc hai kết hợp phần tử hữu hạn làm trơn trên điểm nút cải tiến kỹ thuật chênh lệch cắt rời rạc (NS-DSG3). Đề tài khảo sát toàn diện trên 4 dạng hình học tấm gồm tấm hình vuông kích thước 10 m x 10 m, tấm hình chữ nhật 5 m x 10 m, tấm hình tròn bán kính 6 m và tấm hình chữ L dưới các điều kiện biên ngàm và tựa đơn tự do, với dải độ mảnh biến thiên từ L/t = 1 đến L/t = 100 trên vật liệu thép có giới hạn chảy dẻo 250 MPa.

Ý nghĩa của nghiên cứu thể hiện qua việc giảm thiểu số lượng biến tối ưu về đúng bằng số nút lưới (từ 441 đến 1225 nút), qua đó tiết kiệm đáng kể thời gian xử lý dữ liệu và đạt độ chính xác cao với sai số chỉ từ 0,77% đến 1,86% so với các phương pháp không lưới tiên tiến.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng vững chắc của hai lý thuyết cơ học kinh điển: lý thuyết biến dạng cắt bậc nhất Mindlin-Reissner (FSDT) và lý thuyết phân tích giới hạn dẻo trực tiếp. Mô hình chuyển vị 5 bậc tự do bao gồm 3 thành phần chuyển vị thẳng theo các trục tọa độ và 2 thành phần góc xoay quanh trục, cho phép kể đến đầy đủ chuyển vị màng trong mặt phẳng cũng như biến dạng cắt ngoài mặt phẳng. Đây là ưu điểm vượt trội khi phân tích các kết cấu tấm có cấu tạo bất đối xứng theo chiều dày.

Mô hình vật liệu được giả định tuân theo ứng xử cứng - dẻo lý tưởng với giới hạn chảy dẻo danh định 250 MPa, thỏa mãn định đề ổn định Drucker và luật chảy dẻo kết hợp tương ứng với tiêu chuẩn chảy dẻo Von Mises trong trạng thái ứng suất phẳng. Ba khái niệm then chốt được tích hợp trong mô hình gồm: kỹ thuật Discrete Shear Gap (DSG3) nhằm triệt tiêu hoàn toàn hiện tượng khóa cắt trên phần tử tam giác 3 nút; phương pháp làm trơn biến dạng quanh điểm nút (NS-FEM) để lấy trung bình trường biến dạng trên các miền hỗ trợ nút; và quy hoạch nón bậc hai (Second-Order Cone Programming - SOCP) giúp giải quyết bài toán tối ưu năng lượng phi tuyến. Tích phân theo chiều dày tấm được thực hiện chính xác thông qua 3 điểm Gauss.

Phương pháp nghiên cứu

Nguồn dữ liệu của nghiên cứu được thiết lập từ mô hình số hóa không gian 2D trên lưới phần tử tam giác nền T3. Cỡ mẫu nghiên cứu được chia thành các mức lưới phân giải tăng dần gồm 441 nút (lưới 21x21), 729 nút (lưới 27x27) và 1225 nút (lưới 33x33). Phương pháp chọn mẫu chia lưới có tính chất đồng đều trên toàn bộ 4 miền hình học khảo sát (vuông, chữ nhật tỷ lệ b/a = 2, tròn và chữ L), tạo điều kiện kiểm chứng chặt chẽ tính hội tụ của nghiệm số khi độ dày tấm mỏng dần từ L/t = 1 đến L/t = 100.

Lý do lựa chọn phương pháp phân tích giới hạn cận trên kết hợp giải thuật SOCP là vì bài toán có thể thiết lập hàm mục tiêu cực tiểu hóa tổng năng lượng tiêu tán dẻo với điều kiện ràng buộc công ngoại lực được chuẩn hóa bằng 1. Việc chuyển đổi bài toán biến phân dẻo phi tuyến về dạng nón bậc hai cho phép tận dụng triệt để sức mạnh của bộ giải chuyên dụng Mosek trong môi trường Matlab, đảm bảo tốc độ hội tụ nhanh và không bị phụ thuộc vào bước lặp thời gian như các phân tích đàn dẻo từng bước thông thường. Toàn bộ quy trình tính toán số và đánh giá đối chuẩn được hoàn thành trong giai đoạn 2018-2019.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích số đã mang lại các phát hiện cốt lõi về hệ số tải trọng phá hoại giới hạn và cơ chế sụp đổ của kết cấu tấm:

• Đối với tấm hình vuông 4 biên ngàm chịu tải phân bố đều, hệ số tải trọng giới hạn đạt giá trị 40,19 đối với tấm mỏng (L/t = 100) và 9,04 đối với tấm dày (L/t = 1) trên hệ lưới 729 nút. Trường hợp 4 biên tựa đơn, hệ số tải trọng đạt 21,656 ở trạng thái tấm mỏng và 9,024 ở trạng thái tấm dày.
• Đối với tấm hình chữ nhật kích thước 5 m x 10 m (b = 2a) 4 biên ngàm, hệ số tải trọng giới hạn đạt mức 53,652 khi tấm mỏng (L/t = 100) và 13,940 khi tấm dày (L/t = 1). Sai số hội tụ giữa các bước tăng độ mảnh giảm mạnh từ 88,72% tại mức L/t = 2 xuống chỉ còn 5,40% tại mức L/t = 100.
• Đối với tấm hình chữ nhật 4 biên tựa đơn, hệ số tải trọng đạt 25,814 ở độ mảnh L/t = 100 và 13,686 ở độ mảnh L/t = 1, với mức độ chênh lệch hội tụ giữa các cấp độ mảnh giảm dần từ 73,58% xuống 0,18%.
• Đối với tấm sàn hình tròn bán kính 6 m, hệ số tải trọng giới hạn đạt 11,142 cho liên kết ngàm chu vi và 5,517 cho liên kết tựa chu vi khi tấm mỏng (2R/t = 100); khi tấm dày (2R/t = 1), giá trị tương ứng lần lượt là 2,429 và 2,426.

Thảo luận kết quả

Nguyên nhân dẫn đến sự khác biệt rõ rệt về hệ số tải trọng và cơ cấu phá hoại giữa tấm dày và tấm mỏng xuất phát từ cơ chế truyền lực nội tại. Khi tấm có chiều dày lớn (L/t = 1 đến L/t = 2), biến dạng cắt chiếm ưu thế tuyệt đối, khiến năng lượng tiêu tán dẻo tập trung cục bộ dọc theo các cạnh biên liên kết, dẫn đến hiện tượng phá hoại cắt trượt biên trước khi mô men uốn kịp phát triển toàn diện. Ngược lại, khi tấm mỏng dần (L/t đạt từ 40 đến 100), biến dạng uốn giữ vai trò chủ đạo, năng lượng dẻo phân tán sâu vào trong lòng bản, hình thành các đường chảy dẻo dạng rẽ quạt 45 độ tại góc đối với tấm vuông hoặc đường nứt xiên hướng tâm đối với tấm tròn.

Khi so sánh với các nghiên cứu công bố trước đây sử dụng phương pháp không lưới Galerkin (EFG) 5 bậc tự do của tác giả Hồ Thị Đoan Trang và phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-DSG3) của tác giả H. Le, kết quả của phương pháp NS-DSG3 thể hiện sự tương đồng xuất sắc. Độ sai lệch chỉ ở mức 0,77% đối với tấm chữ nhật ngàm, 1,80% đối với tấm chữ nhật tựa, và 1,86% đối với tấm vuông ngàm. Các tập dữ liệu tính toán này được trình bày trực quan thông qua biểu đồ đường cong hội tụ tải trọng theo độ mảnh và biểu đồ mật độ phân bố năng lượng tiêu tán dẻo 2D, chứng minh tính ổn định tuyệt đối và độ tin cậy của mô hình số.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu lý thuyết và thực nghiệm số, các đề xuất cụ thể được đưa ra nhằm thúc đẩy ứng dụng công nghệ tính toán tiên tiến vào thực tiễn kỹ thuật xây dựng:

• Tích hợp mô hình phần tử NS-DSG3 vào các phần mềm kết cấu chuyên dụng thương mại như SAP2000, ETABS và ABAQUS nhằm rút ngắn 30% đến 45% thời gian phân tích dẻo phi tuyến cho kết cấu sàn đáy và móng bè trong vòng 12 tháng tới, do các nhóm phát triển phần mềm và viện nghiên cứu công nghệ xây dựng chủ trì thực hiện.
• Chuẩn hóa quy trình thiết kế cấu kiện tấm dày bê tông cốt thép và tấm kim loại theo phương pháp phân tích giới hạn trực tiếp, hướng đến mục tiêu tiết kiệm 15% đến 20% khối lượng vật liệu chịu lực tại các dự án nhà cao tầng trong giai đoạn 2024-2026, do Bộ Xây dựng và các hiệp hội tư vấn thiết kế ban hành.
• Mở rộng ứng dụng kỹ thuật làm trơn nút NS-FEM cho bài toán kết cấu vật liệu composite lớp và vật liệu biến tính chức năng (FGM) có cấu tạo bất đối xứng, đặt mục tiêu nâng cao độ chính xác dự báo phá hủy đạt trên 98% trong thời gian nghiên cứu 18 tháng, do các phòng thí nghiệm trọng điểm cơ học công trình triển khai.
• Tổ chức các khóa đào tạo chuyên sâu về ứng dụng thuật toán tối ưu hóa hình nón bậc hai (SOCP) trên nền tảng ngôn ngữ tính toán hiện đại cho tối thiểu 500 kỹ sư kết cấu công trình hàng năm, do Hội Kết cấu và Công nghệ Xây dựng Việt Nam phối hợp cùng các trường đại học kỹ thuật tổ chức.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả khoa học của luận văn mang lại giá trị thực tiễn và học thuật cao cho nhiều nhóm độc giả chuyên ngành:

• Kỹ sư thiết kế kết cấu công trình: Ứng dụng trực tiếp thuật toán và bảng tra hệ số tải trọng để tính toán tối ưu khả năng chịu lực của các cấu kiện tấm dày từ 1,0 m đến 2,5 m như móng bè, nắp hầm ngầm, sàn chuyển tải và bể chứa công nghiệp mà không lo ngại sai số khóa cắt.
• Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Xây dựng và Cơ học kỹ thuật: Sử dụng luận văn làm tài liệu tham khảo mẫu mực về quy trình thiết lập bài toán tối ưu hình nón bậc hai (SOCP), phương pháp làm trơn biến dạng (Smoothed FEM) và kỹ thuật xử lý số với quy mô lưới từ 441 đến 1225 nút.
• Giảng viên và chuyên gia nghiên cứu tính toán cơ học: Tham khảo khung lý thuyết tấm Mindlin 5 bậc tự do để biên soạn bài giảng chuyên đề sau đại học, đồng thời phát triển các hướng nghiên cứu mở rộng cho bài toán dao động phi tuyến và cơ học phá hủy dẻo.
• Chuyên gia thẩm tra và kiểm định an toàn chất lượng công trình: Áp dụng quy luật phân bố năng lượng tiêu tán dẻo để dự đoán chính xác vị trí xung yếu, vết nứt cục bộ và cơ cấu sụp đổ tiềm ẩn trên các kết cấu tấm vượt nhịp trên 10 m chịu tải trọng nặng.

Câu hỏi thường gặp

Tại sao lý thuyết tấm dày 5 bậc tự do lại vượt trội hơn lý thuyết 3 bậc tự do trong phân tích giới hạn?
Lý thuyết 5 bậc tự do mô tả đồng thời cả chuyển vị màng trong mặt phẳng và chuyển vị uốn, cắt ngoài mặt phẳng. Điều này cho phép phân tích chính xác các tấm có cấu tạo bất đối xứng theo chiều dày, giúp loại bỏ sai số từ 15% đến 25% vốn thường xuất hiện trong lý thuyết 3 bậc tự do khi chiều dày tấm tăng cao với tỷ lệ L/t nhỏ hơn 10.

Hiện tượng khóa cắt (shear locking) trong tấm mỏng được triệt tiêu bằng cách nào?
Hiện tượng khóa cắt xuất hiện khi tấm mỏng dần đến mức L/t = 100 được xử lý triệt để nhờ tích hợp kỹ thuật Discrete Shear Gap (DSG3) trên phần tử tam giác 3 nút. Kỹ thuật này rời rạc hóa độc lập trường biến dạng cắt, giúp kết quả tính toán hội tụ mượt mà và duy trì sai số tải trọng dưới 1,85% trên toàn dải độ mảnh.

Ưu thế nổi bật của phương pháp tối ưu nón bậc hai (SOCP) so với quy hoạch tuyến tính (LP) là gì?
Thay vì phải xấp xỉ tuyến tính đa diện mặt dẻo phi tuyến gây sai số tích lũy từ 5% đến 10% như phương pháp LP cổ điển, thuật toán SOCP mô tả chính xác tiêu chuẩn Von Mises phi tuyến bậc hai. Khi kết hợp cùng bộ giải Mosek, thuật toán đạt tốc độ xử lý nhanh hơn 40% và giảm mạnh số biến tối ưu.

Cơ cấu phá hoại của tấm sàn dày khác biệt như thế nào so với tấm mỏng?
Ở tấm có độ dày lớn (L/t = 1), biến dạng cắt đóng vai trò quyết định, làm năng lượng tiêu tán dẻo tập trung hoàn toàn dọc theo chu vi biên liên kết gây phá hoại cắt cục bộ. Ngược lại, ở tấm mỏng (L/t = 100), biến dạng uốn chi phối tạo nên các đường chảy dẻo rẽ quạt 45 độ tại góc hoặc xiên chéo vào tâm bản.

Việc sử dụng 3 điểm Gauss dọc theo chiều dày tấm mang lại lợi ích gì về mặt tính toán?
Nghiên cứu áp dụng 3 điểm Gauss dọc theo chiều dày để lấy tích phân năng lượng tiêu tán dẻo. Cấu hình 3 điểm Gauss này đảm bảo độ chính xác tương đương các mô hình 5 hay 7 điểm Gauss với sai số dưới 1,0%, đồng thời tiết kiệm 35% dung lượng bộ nhớ máy tính trong quá trình tối ưu hóa.

Kết luận

• Luận văn đã xây dựng thành công mô hình phân tích giới hạn cận trên cho kết cấu tấm dày 5 bậc tự do theo tiêu chuẩn chảy dẻo Von Mises trên vật liệu có giới hạn chảy 250 MPa.
• Sự kết hợp giữa phần tử làm trơn nút NS-FEM và kỹ thuật chênh lệch cắt DSG3 đã loại bỏ hoàn toàn hiện tượng khóa cắt, bảo đảm độ chính xác cao trên toàn dải độ mảnh L/t từ 1 đến 100.
• Ứng dụng giải thuật quy hoạch nón bậc hai (SOCP) cùng bộ giải Mosek giúp tối ưu hóa số lượng biến số trên các hệ lưới từ 441 đến 1225 nút, đạt sai số ấn tượng dưới 1,86% so với các phương pháp không lưới.
• Nghiên cứu đã làm rõ quy luật biến chuyển cơ chế phá hoại dẻo từ phá hủy cắt trượt cục bộ dọc biên ở tấm dày sang hình thành các đường chảy dẻo rẽ quạt uốn ở tấm mỏng.
• Kết quả cung cấp hệ thống số liệu chuẩn mực về tải trọng giới hạn cho 4 dạng hình học tấm cơ bản gồm vuông 10 m, chữ nhật, tròn bán kính 6 m và tấm chữ L.

Đóng góp chính của luận văn là thiết lập một công cụ tính toán số tin cậy, hiệu quả cao cho bài toán phân tích giới hạn tấm dày 5 bậc tự do trong kỹ thuật xây dựng hiện đại. Trong giai đoạn 2024-2025, định hướng nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình cho kết cấu tấm composite phân lớp và tấm có lỗ khoét phức tạp. Kính mời các nhà khoa học, kỹ sư kết cấu và học viên chuyên ngành tham khảo toàn văn công trình để ứng dụng trực tiếp vào công tác nghiên cứu và thực tiễn thiết kế công trình.

Trích đoạn nội dung tài liệu

mở đầu đã trình bày vắn tắt lí do chọn đề tài cũng như các ý nghĩa khoa học và thực tiễn của đề tài. Qua đó, hướng nghiên cứu được định hình và bám sát với mục tiêu của luận văn. Hơn thế nữa, việc giới hạn phạm vi nghiên cứu là cần thiết phù hợp với thời gian làm luận văn. Tiếp theo sẽ là chương “Tổng quan” nhằm tìm hiểu các nghiên cứu của các nhà khoa học trong và ngoài nước gần với hướng nghiên cứu của đề tài.

Tổng quan tài liệu trong nước Nghiên cứu trong nước về phân tích giới hạn vẫn chưa được triển khai nhiều. Nguyen et al. (1998) đã thực hiện phân tích giới hạn và phân tích thích nghi dẻo và nứt cho kết cấu. Phân tích giới hạn cho khung được thực hiện bởi D.

Nguyen et al. Phân tích giới hạn sử dụng thuật toán đối ngẫu được tác giả C.V Le et al. (2006) cho tấm mỏng chịu uốn. Nguyen (2012) đã áp dụng phương pháp EFG, tiêu chuẩn dẻo Von Mises và chương trình hình nón bậc hai (SOCP) cho bài toán phân tích giới hạn tấm dày Mindlin-Reissner.

Cũng thời gian này, T. Doan (2012) đã sử dụng phương pháp FEM trơn dựa trên cạnh ES-DSG3 và chương trình hình nón bậc hai (SOCP) cũng cho bài toán phân tích giới hạn cận trên tấm dày Minlin-Reissner. Tổng quan tài liệu nước ngoài Trên thế giới, phương pháp phân tích giới hạn được phát triển từ đầu thế kỉ XX dựa trên mô hình vật liệu cứng – dẻo lý tưởng và trở thành công cụ rất hiệu quả để thiết kế và đánh giá độ an toàn của kết cấu. Những phát biểu hoàn chỉnh đầu tiên về lý thuyết phân tích giới hạn được giới thiệu bởi Drucker et al.

Sau đó, Prager (1972) và Martin (1975) đã có những đóng góp quan trọng để phát triển và áp dụng phương pháp phân tích giới hạn. Việc áp dụng lý thuyết phân tích giới hạn vào phân tích kết cấu được thực hiện bởi Hodge (1959, 1961, 1963). Một thời gian ngắn sau đó, Massonet & Save (1967) và Save & Massonet (1972) đã có những phát triển xa hơn. Trong vài thập niên qua, lý thuyết phân tích giới hạn thu hút sự quan tâm của nhiều nhà nghiên cứu như Charakbarty (1988), Chen & Han (1988), Lubliner (1990), Kamenjarzh (1996) và gần đây nhất là Jirasek & Bazant (2002) cho nhiều bài toán kỹ thuật.

Chương trình tuyến tính (LP) được sử dụng rộng rãi trong nhiều bài toán phân tích giới hạn vì thuật toán sử dụng đơn giản (Anderheggen & Knopfel, 1972; Cohn & Maier, 1979; Grierson, 1977; Nguyen Dang, 1984; Sloan, 1988). Một số tác giả đã sử dụng chương trình tuyến tính để tuyến tính một phần của mặt dẻo phi tuyến như Maier (1970), Laudiero (1972), Tin – Loi 7 (1970) và Christiansen (1996). Trở ngại lớn nhất khi thực hiện bài toán phân tích giới hạn đó là giải quyết bài toán tối ưu với ràng buộc phi tuyến tính. Với sự phát triển về các kỹ thuật toán tối ưu hiện nay đã có thể giải quyết được vấn đề trên nên các nghiên cứu gần đây có số lượng gia tăng đáng kể.

Những bài toán phân tích giới hạn sử dụng tiêu chuẩn dẻo có ràng buộc phi tuyến được giải quyết bằng các thuật toán tối ưu phi tuyến như thuật toán Sequential quadratic programming (SQP) với phương pháp Newton suy rộng cho tối ưu không ràng buộc (Gaudrat, 1991; Zouain et al., 1993; Liu et al. 1995; Andersen, 1996) hoặc thuật toán trực tiếp (Capsoni & Corradi, 1999). Thuật toán tối ưu chương trình hình nón bậc hai (SOCP) được sử dụng bởi Lobo et al. (1998), Andersen et al.

Le et al. Những vấn đề còn tồn tại và mục tiêu, nhiệm vụ luận văn Qua các nghiên cứu gần đây, sự phát triển của các lý thuyết tấm khác nhau và vật liệu ngày càng đa dạng. Điều này dẫn đến nhu cầu cần áp dụng các phương pháp khác nhau giúp mở rộng các trường hợp tính toán của bài toán phân tích giới hạn cho tấm dày. Với số lượng bậc tự do tăng lên đồng nghĩa với việc gia tăng các biến của bài toán.

Việc sử dụng các phương pháp khác nhau nhằm tạo điều kiện cho việc tận dụng các điểm mạnh của các phương pháp khác nhau để giải bài toán. Trong phạm vi luận văn này, chúng ta sẽ phân tích giới hạn tấm dày năm bậc tự do theo tiêu chuẩn chảy dẻo Von Mises theo cận trên với sự hỗ trợ của phương pháp NS-DSG3 bao gồm các giai đoạn:  Thiết lập công thức tính tiêu tán chảy dẻo tương ứng với tiêu chuẩn dẻo Von Mises. Sau đó thành lập bài toán tối ưu toán học cho việc phân tích giới hạn cận trên.  Chuyển đổi các công thức tìm được ở trên về dạng của bài toán tối ưu hình nón bậc hai.

 Thiết lập công thức xấp xỉ số dùng phương pháp NS-DSG3.  Lập trình và thực hiện tính toán số cho các bài toán tấm sàn đã được khảo sát trong các bài báo đã được xuất bản.  Phân tích và so sánh kết quả thu được với kết quả số khác. 8 Tổng kết: Chương này đã trình bày các nghiên cứu của các nhà khoa học trong và ngoài nước có liên quan đến hướng nghiên cứu của luận văn.

Khi đó tính mới và cấp thiết của luận văn được thể hiện một cách mạnh mẽ khi hướng nghiên cứu này được đăng ở các bài báo quốc tế. Qua đó các giai đoạn hình thành luận văn được thiết lập nhằm đảm bảo tiến độ thực hiện luận văn. Phần sau, các cơ sở lý thuyết được trình bày một cách rõ ràng. Hơn thế nữa, các lý thuyết của bài toán phân tích giới hạn cho tấm dày 5 bậc tự do được xây dựng.

Mô hình vật liệu Ứng xử của vật liệu được thể hiện qua đường cong quan hệ ứng suất – biến dạng.1, ứng xử của hai loại vật liệu: vật liệu mềm (cốt thép) và vật liệu dòn (bê tông) được trình bày. Mô hình vật liệu Việc phân tích kết cấu với ứng xử thật của vật liệu gặp phải nhiều khó khăn, vì vậy, trong lý thuyết dẻo, quy luật ứng xử của vật liệu được đưa về những mô hình lý tưởng. Sự lý tưởng hóa này làm cho việc phân tích kết cấu trở nên đơn giản, đặc biệt là khi sử dụng lý thuyết phân tích giới hạn, tải trọng phá hủy của kết cấu được xác định nhanh chóng mà bỏ qua quá trình hình thành cơ cấu phá hoại. Ứng xử của vật liệu có thể được đưa về nhiều mô hình khác nhau như: đàn – dẻo lý tưởng, cứng – dẻo lý tưởng, đàn – dẻo với tái bền tuyến tính, cứng – dẻo với tái bền tuyến tính… Trong luận văn này, giả sử rằng ứng xử của vật liệu tuân theo mô hình cứng – dẻo lý tưởng, nghĩa là bỏ qua giai đoạn đàn hồi, tái bền và mềm hóa, mối quan hệ ứng suất – biến dạng là đường thẳng nằm ngang.

Mô hình vật liệu cứng -dẻo lý tưởng 3. Định đề ổn định của Drucker “ Đối với một chu trình khép kín của tải trọng, công sinh ra không âm”.1) trong đó: σ - trường ứng suất hiện tại, σ0 - trường ứng suất ban đầu, ε - trường biến dạng. Ứng xử ổn định và không ổn định theo Drucker (a) vật liệu ổn định (b), (c) vật liệu không ổn định. Theo biểu thức (3.1) vật liệu có ứng xử ổn định gây ra công dương, vật liệu có ứng xử không ổn định gây ra công âm.3, vật liệu thỏa ứng xử (a) gọi là vật liệu ổn định (vật liệu chuẩn), vật liệu ứng xử như (b) và (c) gọi là vật liệu không ổn định (vật liệu không chuẩn).

Luật chảy dẻo kết hợp “Vectơ gia số biến dạng dẻo (vận tốc biến dạng dẻo) vuông góc với mặt tải và hướng ra ngoài mặt này”. Điều kiện: vectơ gia số biến dạng dẻo ε p phải nằm trên mặt chảy dẻo, tại điểm trơn hoặc nằm giữa những điểm kế cận tại một điểm góc. Tại những điểm lồi nghiêm ngặt (lồi và trơn), luật chảy dẻo kết hợp có thể biểu diễn dưới dạng:  i εp  μi (3.2) σ trong đó, μ là hệ số chảy dẻo, ψ là hàm chảy dẻo. Tại những điểm kì dị (không trơn), biểu thức (3.2) có thể biểu diễn như sau: n  i εp   μi (3.

Luật chảy dẻo kết hợp Tùy theo tiêu chuẩn dẻo áp dụng vào biểu thức (3.2) ta có luật chảy dẻo kết hợp tương ứng và năng lượng tiêu tán dẻo tương ứng. Lý thuyết tấm dày Lý thuyết cắt bậc nhất (Lý thuyết tấm dày Mindlin-Reissner): Đối với tấm dày, pháp tuyến của mặt trung hòa vẫn thẳng nhưng không nhất thiết phải vuông góc với mặt phẳng này trong suốt quá trình biến dạng. Lý thuyết của Mindlin-Reissner được sử dụng để tính toán biến dạng và ứng suất trong tấm dày trong khi lý thuyết Kirchhoff-Love được áp dụng cho các tấm mỏng. Lý thuyết của Mindlin giả định rằng có một sự biến đổi tuyến tính của 12 chuyển vị thông qua độ dày tấm nhưng độ dày tấm không thay đổi trong suốt quá trình biến dạng.

Điều này ngụ ý rằng, ứng suất pháp thông qua chiều dày tấm bị bỏ qua, trạng thái ứng suất phẳng xảy ra. Mặt khác, lý thuyết của Reissner giả định rằng ứng suất uốn là tuyến tính trong khi ứng suất cắt là bậc hai trên suốt chiều dày của tấm. Điều này dẫn đến một trạng thái khác mà ở đó chuyển vị trên suốt chiều dày tấm không nhất thiết phải tuyến tính và độ dày tấm có thể thay đổi trong quá trình biến dạng. Lý thuyết tấm dày Mindlin-Reissner là sự kết hợp lý thuyết của hai người và từ đây thường được gọi là lý thuyết biến dạng cắt bậc nhất.

Cùng xem xét miền thể tích trong R 2 với mặt phẳng giữa tấm (mặt trung bình) Hình 3. Kết cấu tấm và qui ước chiều ứng suất Trường chuyển vị theo lý thuyết FSDT [5] gồm 5 bậc tự do như sau: u  x , y, z   u0  x , y   z  x  x , y  v  x , y, z   v0  x , y   z  y  x , y  (3.4) w  x , y, z   w0  x , y  Biến dạng trong mặt phẳng được biểu hiện theo công thức T  xx  yy  xy   ε 0  zκ (3.5) Với biến dạng màng và biến dạng cong ε0   su0 (3.7) Biến dạng cắt 13 ε s  w  β (3. Lý thuyết phân tích giới hạn 3.