Luận văn thạc sĩ kỹ thuật xây dựng phân tích giới hạn cho tấm mindlin reisser dùng phương pháp es dsg3 kết hợp với kỹ thuật thích nghi lưới

Luận văn phân tích giới hạn tấm Mindlin-Reissner bằng phương pháp ES-DSG3 và kỹ thuật thích nghi lưới. Nghiên cứu chuyên sâu ngành kỹ thuật xây dựng.

• Tác giả: Trần Thế Vân
• Người hướng dẫn: PGS.TS Nguyễn Thị Hiền Lương
• Trường học: Trường Đại học Bách Khoa – ĐHQG-HCM
• Chuyên ngành: Xây dựng công trình dân dụng và công nghiệp
• Đề tài: Phân Tích Giới Hạn Cho Tấm Mindlin-Reissner Dùng Phương Pháp ES-DSG3 Kết Hợp Với Kỹ Thuật Thích Nghi Lưới.
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: Tp. Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• DANH MỤC HÌNH VẼ
• DANH MỤC BẢNG BIỂU
• 1. Chương 1: Tình hình nghiên cứu và tính cấp thiết đề tài
• 1.1. Trên thế giới
• 1.2. Mục tiêu và nhiệm vụ của luận văn
• 1.3. Cấu trúc luận văn
• 2. Chương 2: Cơ sở lý thuyết
• 2.1. Lý thuyết tấm dày Mindlin-Reissner
• 2.2. Lý thuyết biến dạng cắt bậc nhất (FSDT)
• 2.3. Mô hình vật liệu
• 2.3.1. Mô hình đàn hồi - dẻo lý tưởng
• 2.3.2. Mô hình cứng dẻo - lý tưởng
• 2.4. Tiêu chuẩn Von Mises
• 2.5. Luật chảy dẻo
• 2.6. Định đề Drucker – giả thiết về tính ổn định của vật liệu
• 2.7. Luật chảy dẻo kết hợp
• 2.8. Nguyên lý công ảo
• 2.8.1. Trường ứng suất khả dĩ tĩnh
• 2.8.2. Trường chuyển vị, biến dạng khả dĩ động
• 3. Chương 3: PHÂN TÍCH GIỚI HẠN TẤM MINDLIN - REISSNER BẰNG PHƯƠNG PHÁP ES -DSG3 VỚI KỸ THUẬT THÍCH NGHI LƯỚI
• 3.1. Tổng quan về phân tích giới hạn
• 3.2. Công thức phân tích giới hạn tấm dày Mindlin-Reissner
• 3.3. Phương pháp "rời rạc lệch trượt" – Discrete Shear Gap (DSG3)
• 3.4. Phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM)
• 3.5. Rời rạc phương trình phân tích động (Định lý cận trên)
• 3.6. Kỹ thuật tối ưu hình nón bậc hai (Second order cone programming– SOCP)
• 3.7. Biến đổi bài toán tối ưu về dạng hình nón bậc hai
• 3.8. Gán các biến về dạng hình nón trong Mosek
• 3.9. Kỹ thuật thích nghi lưới
• 3.9.1. Tính sai số rời rạc dựa vào năng lượng tiêu tán dẻo
• 3.9.2. Cách thức làm mịn lưới
• 3.10. Sơ đồ khối cho bài toán phân tích giới hạn tấm Mindlin-Reissner bằng phương pháp ES-DSG3 kết hợp với kỹ thuật thích nghi lưới
• 4. Chương 4: Ví dụ số
• 4.1. Tấm hình vuông
• 4.2. Tấm hình chữ nhật
• 4.3. Tấm hình thoi
• 4.4. Tấm hình tam giác
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Phân Tích Giới Hạn Tấm Mindlin Reissner Tổng Quan 58 Ký Tự

Phân tích kết cấu đàn dẻo đến khi phá hoại là một quá trình phức tạp. Cần gia tăng tải trọng tuần tự. Quá trình này giúp hiểu rõ sự hình thành và phát triển biến dạng dẻo. Tuy nhiên, tính toán số tốn kém và mất thời gian. Một hướng khác giúp tìm tải trọng giới hạn khi kết cấu phá hoại. Đó là Lý thuyết phân tích trực tiếp tải trọng giới hạn, hay phân tích giới hạn. Phân tích giới hạn cung cấp trị số tải trọng giới hạn. Kỹ sư dùng nó để thiết kế và đánh giá độ an toàn kết cấu. Phương pháp phần tử hữu hạn (FEM) là phương pháp xấp xỉ số mạnh. Nó được sử dụng rộng rãi trong phân tích giới hạn tấm dày. FEM có ưu điểm: thiết lập và giải phương trình đơn giản, giải được nhiều bài toán (dao động, ổn định) với nhiều hình dạng, áp dụng được nhiều lý thuyết giải. Tuy nhiên, FEM có hạn chế là hiện tượng shear locking. Shear locking xảy ra với phần tử tấm Reissner-Mindlin bậc thấp khi chiều dày tấm mỏng dần. Điều này dẫn đến kết quả phân tích không chính xác hoặc không hội tụ.

1.1. Lý thuyết Tấm Mindlin Reissner và Ứng Dụng 44 Ký Tự

Lý thuyết tấm Mindlin-Reissner xét đến biến dạng cắt ngang. Điều này quan trọng khi phân tích tấm dày, nơi biến dạng cắt đóng vai trò đáng kể. Mô hình Kirchhoff bỏ qua biến dạng cắt, phù hợp với tấm mỏng. Lý thuyết Mindlin-Reissner chính xác hơn cho tấm dày nhưng phức tạp hơn. Việc lựa chọn lý thuyết phù hợp phụ thuộc vào tỷ lệ giữa chiều dày và kích thước của tấm. Ứng dụng của lý thuyết này rất đa dạng, từ thiết kế cầu, vỏ tàu đến các cấu kiện trong ngành hàng không vũ trụ.

1.2. Khắc Phục Shear Locking với ES FEM 40 Ký Tự

Để khắc phục shear locking, phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM) kết hợp với phương pháp lệch trượt (ES-DSG) được sử dụng. Phương pháp rời rạc lệch trượt (DSG3) làm trơn hóa biến dạng và xử lý hiệu quả shear locking. ES-FEM cho lời giải ổn định, chính xác hơn và hội tụ nhanh hơn so với FEM truyền thống. Phương pháp này bỏ qua bước xây dựng ma trận hàm dạng chính xác và áp dụng cho nhiều loại phần tử (tam giác, tứ giác,...).

II. ES DSG3 Phương Pháp Hiệu Quả Phân Tích Giới Hạn 55 Ký Tự

Phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM) kết hợp với phương pháp lệch trượt (ES-DSG), cụ thể là DSG3, được chứng minh hiệu quả trong giải quyết bài toán phân tích giới hạn cho tấm Mindlin. Theo Bletzinger và Bischoff (2000), DSG3 gần giống phương pháp giả định biến dạng tự nhiên ANS và làm việc tốt cho phần tử có hình dáng và bậc khác nhau. Giáo sư Liu đề xuất phương pháp phần tử hữu hạn trơn. Phương pháp này kế thừa FEM truyền thống và có kỹ thuật làm trơn hóa biến dạng. Có nhiều loại FEM trơn: dựa trên phần tử (CS-FEM), dựa trên cạnh (ES-FEM), dựa vào nút (NS-FEM). Trong đó, ES-FEM cho lời giải ổn định và có nhiều ưu điểm. Các nghiên cứu trước đây cho thấy ES-DSG cho kết quả tốt, hội tụ nhanh và giải quyết triệt để shear locking.

2.1. Ưu Điểm Vượt Trội Của ES DSG3 34 Ký Tự

ES-DSG3 có nhiều ưu điểm vượt trội so với các phương pháp khác. Nó cho kết quả chính xác hơn, kết quả hội tụ nhanh hơn với phần tử hữu hạn truyền thống. Phương pháp này bỏ qua bước xây dựng ma trận hàm dạng chính xác, giảm thiểu đáng kể chi phí tính toán. Nó có thể áp dụng cho nhiều loại phần tử khác nhau, từ tam giác, tứ giác đến đa giác n cạnh, tăng tính linh hoạt trong mô hình hóa.

2.2. Kỹ Thuật Tối Ưu Hình Nón Bậc Hai SOCP 42 Ký Tự

Để giải quyết bài toán phân tích giới hạn, cần sử dụng thuật toán tối ưu. Trong luận văn này, kỹ thuật tối ưu hình nón bậc hai (SOCP) được sử dụng với sự hỗ trợ của phần mềm thương mại Mosek. SOCP phù hợp với hầu hết các tiêu chuẩn chảy dẻo (ví dụ Von Mises) và giải quyết hầu hết các bài toán, kể cả bài toán lớn. SOCP mang lại kết quả tương thích với các phương pháp khác và tiết kiệm thời gian, chi phí tính toán. Theo các nghiên cứu trước đây [1, 4, 14, 16] đã chỉ ra rằng nó còn tiết kiệm thời gian và chi phí tính toán.

III. Thích Nghi Lưới Cách Tối Ưu Hóa Độ Chính Xác 53 Ký Tự

Độ chính xác của lời giải phụ thuộc vào quá trình rời rạc hóa. Mục tiêu là áp dụng sai số và thích nghi để hoàn thiện và gia tăng quá trình tính toán. Hoặc kết hợp cả hai để thu được kết quả tốt nhất. Phương pháp làm mịn lưới tự động dựa vào đánh giá hậu sai số. Đánh giá hậu sai số được phân loại thành 2 trường hợp: kiểu phục hồi (recovery-type) và kiểu dư (residual-type). Sai số kiểu phục hồi đơn giản và được ứng dụng rộng rãi trong những bài toán thực tiễn. Lưới sẽ được thích nghi tại những vùng mà có sai số lớn. Có một vài kỹ thuật thích nghi lưới như: làm mịn lưới (h-adaptivity), làm giàu lưới (p-adaptivity) và sắp xếp lại lưới (r-adaptivity). Kỹ thuật làm mịn lưới đơn giản nên được áp dụng rộng rãi trong nhiều lĩnh vực kỹ thuật.

3.1. Các Phương Pháp Thích Nghi Lưới Phổ Biến 44 Ký Tự

Có nhiều kỹ thuật thích nghi lưới khác nhau. Làm mịn lưới (h-adaptivity) chia nhỏ các phần tử có sai số lớn. Làm giàu lưới (p-adaptivity) tăng bậc của hàm nội suy trên các phần tử có sai số lớn. Sắp xếp lại lưới (r-adaptivity) di chuyển các nút lưới để giảm sai số. Mỗi phương pháp có ưu điểm và nhược điểm riêng. Việc lựa chọn phương pháp phù hợp phụ thuộc vào bài toán cụ thể.

3.2. Kỹ Thuật Làm Mịn Lưới h adaptivity 36 Ký Tự

Trong Luận văn này tác giả sử dụng kỹ thuật (h-adaptivity) để phân tích. Sai số do rời rạc được xác định đơn giản. Đó chính là kết quả năng lượng tiêu tán dẻo phần tử thu được từ chương trình hình nón bậc hai (SOCP) với sự hỗ trợ của phần mềm thương mại Mosek. Sau khi đánh dấu được các phần tử tam giác có năng lượng tiêu tán dẻo lớn, ta tiến hành kỹ thuật làm mịn lưới. Một tam giác được phân chia thành hai tam giác con nối với nhau tại điểm giữa cạnh biên dài nhất. Một góc mới tạo ra tại điểm giữa cạnh đáy được quy định đỉnh của tam giác con đó.

IV. Ứng Dụng Kết Quả Phân Tích Tấm Với ES DSG3 59 Ký Tự

Tác giả khảo sát bài toán phân tích giới hạn tấm Mindlin bằng phương pháp ES-DSG3 với kỹ thuật thích nghi lưới. Kết quả thu được khá tốt. Bài toán được thực hiện trên nhiều hình dạng tấm khác nhau. Tấm hình vuông, chữ nhật và thoi được sử dụng để kiểm chứng phương pháp. Kết quả được so sánh với các nghiên cứu trước đó. So sánh này nhằm đánh giá độ chính xác và hiệu quả của phương pháp. Việc thích nghi lưới giúp cải thiện độ chính xác của kết quả. Các vùng có ứng suất cao được chia nhỏ phần tử. Điều này giúp mô tả chính xác hơn sự phân bố ứng suất trong tấm. Kỹ thuật này đặc biệt quan trọng trong phân tích giới hạn.

4.1. Khảo Sát Tấm Hình Vuông và Chữ Nhật 41 Ký Tự

Tấm hình vuông và chữ nhật được sử dụng rộng rãi trong các kết cấu xây dựng. Việc phân tích giới hạn cho các loại tấm này rất quan trọng. Các điều kiện biên khác nhau được áp dụng. Ngàm, tựa đơn và kết hợp cả hai được xét đến. Kết quả cho thấy ES-DSG3 với thích nghi lưới cho kết quả chính xác. Kết quả này phù hợp với các nghiên cứu trước đó.

4.2. Phân Tích Tấm Hình Thoi và Tam Giác 40 Ký Tự

Ngoài tấm hình vuông và chữ nhật, tấm hình thoi và tam giác cũng được khảo sát. Các hình dạng phức tạp này thử thách khả năng của phương pháp. ES-DSG3 với thích nghi lưới vẫn cho kết quả tốt. Điều này chứng tỏ tính ứng dụng rộng rãi của phương pháp. Kết quả cho thấy tiềm năng ứng dụng của phương pháp này trong các bài toán kỹ thuật phức tạp.

V. Kết Luận Hướng Phát Triển Phân Tích Giới Hạn 56 Ký Tự

Luận văn đã trình bày phương pháp phân tích giới hạn cho tấm Mindlin-Reissner. Phương pháp sử dụng ES-DSG3 kết hợp với kỹ thuật thích nghi lưới. Kết quả cho thấy phương pháp có độ chính xác cao và hiệu quả. Shear locking được giải quyết triệt để. Kỹ thuật thích nghi lưới giúp tối ưu hóa độ chính xác và giảm chi phí tính toán. Phương pháp này có tiềm năng ứng dụng rộng rãi trong thiết kế và đánh giá độ an toàn kết cấu. Nghiên cứu này mở ra hướng phát triển mới trong lĩnh vực phân tích kết cấu.

5.1. Đánh Giá Hiệu Quả của Thích Nghi Lưới 40 Ký Tự

Thích nghi lưới giúp cải thiện đáng kể độ chính xác của kết quả. Các phần tử được chia nhỏ ở những vùng có ứng suất cao. Điều này giúp mô tả chính xác hơn sự phân bố ứng suất. Kết quả là, tải trọng giới hạn được xác định chính xác hơn. Thích nghi lưới cũng giúp giảm chi phí tính toán. Số lượng phần tử được giảm thiểu ở những vùng có ứng suất thấp.

5.2. Hướng Nghiên Cứu và Phát Triển Tiếp Theo 42 Ký Tự

Nghiên cứu này có thể được mở rộng theo nhiều hướng. Áp dụng phương pháp cho các loại tấm khác nhau. Nghiên cứu ảnh hưởng của vật liệu đến kết quả. Phát triển các thuật toán thích nghi lưới hiệu quả hơn. Nghiên cứu ứng dụng của phương pháp trong các bài toán thực tế. Các hướng nghiên cứu này sẽ giúp nâng cao hơn nữa giá trị ứng dụng của phương pháp.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, việc xác định khả năng chịu lực tới hạn của kết cấu đóng vai trò then chốt trong thiết kế an toàn. Phương pháp phân tích đàn dẻo từng bước truyền thống thường tiêu tốn thời gian tính toán gấp 3 đến 5 lần do phải gia tăng tải trọng tuần tự qua hàng trăm bước lặp phi tuyến. Phân tích giới hạn trực tiếp xuất hiện như một công cụ đột phá, cho phép xác định trực tiếp tải trọng phá hoại dẻo mà không cần qua các bước trung gian tốn kém. Tuy nhiên, khi áp dụng phương pháp phần tử hữu hạn cổ điển cho tấm dày Mindlin-Reissner, hiện tượng khóa trượt (shear locking) xuất hiện nghiêm trọng khi tỷ số chiều dày trên nhịp giảm xuống dưới 0,2 (khoảng 1/5), làm sai lệch kết quả và cản trở sự hội tụ.

Luận văn thạc sĩ của tác giả Trần Thế Vân, được thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong thời gian 6 tháng (từ tháng 01/2013 đến tháng 06/2013), tập trung giải quyết triệt để rào cản này. Mục tiêu nghiên cứu là xây dựng mô hình phân tích giới hạn cận trên cho tấm Mindlin-Reissner bằng cách kết hợp phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM) với kỹ thuật rời rạc lệch trượt (DSG3), đồng thời tích hợp chiến lược tự thích nghi lưới (h-adaptivity) và quy hoạch nón bậc hai (SOCP). Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số tính toán xuống dưới 1% so với nghiệm chuẩn, đồng thời tiết kiệm hơn 40% số lượng phần tử tính toán so với việc chia lưới đều truyền thống.

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 lý thuyết biến dạng cắt bậc nhất (FSDT) của Mindlin-Reissner, giả thiết các đoạn thẳng vuông góc với mặt trung hòa vẫn thẳng sau biến dạng nhưng không còn vuông góc. Để mô tả ứng xử phá hoại của vật liệu, mô hình vật liệu cứng – dẻo lý tưởng theo tiêu chuẩn chảy dẻo von Mises được áp dụng, bỏ qua hiện tượng tái bền và xem biến dạng đàn hồi là không đáng kể.

Các khái niệm then chốt bao gồm:

• Định lý phân tích động (định lý cận trên): Tải trọng phá hoại dẻo được xác định là giá trị nhỏ nhất của hệ số tải trọng trong các trường chuyển vị động học khả dĩ, tại đó công ngoại lực cân bằng với năng lượng tiêu tán dẻo toàn phần.
• Kỹ thuật làm trơn biến dạng dựa trên cạnh (ES-FEM): Miền phần tử tam giác được chia thành các miền con trơn gắn liền với từng cạnh, giúp làm mềm kết cấu một cách tự nhiên và loại bỏ việc xây dựng ma trận hàm dạng phức tạp.
• Kỹ thuật rời rạc lệch trượt (DSG3): Hiệu chỉnh ma trận biến dạng cắt trở thành hằng số trên từng phần tử tam giác 3 nút, triệt tiêu hoàn toàn hiện tượng khóa trượt khi chiều dày tấm tiến về 0.
• Năng lượng tiêu tán dẻo: Tích phân theo chiều dày tấm với 3 điểm tích phân Gauss giúp xấp xỉ chính xác trường biến dạng uốn và cắt.

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

Mẫu nghiên cứu gồm 5 dạng kết cấu tấm hình học phức tạp: tấm vuông, tấm chữ nhật với nhiều điều kiện biên, tấm chữ L có góc lõm, tấm hình thoi góc nghiêng 60 độ và tấm tam giác với lưới ban đầu gồm 75 phần tử tam giác. Phương pháp chọn mẫu tập trung vào các dạng hình học xuất hiện điểm kỳ dị và vùng tập trung ứng suất dẻo cao trong thực tế xây dựng.

Lý do lựa chọn phương pháp phân tích ES-DSG3 kết hợp tối ưu nón bậc hai (SOCP) là khả năng chuyển đổi toàn bộ bài toán tìm cận trên thành dạng hàm mục tiêu nón bậc hai trực tiếp, giải quyết thông qua bộ giải Mosek trên nền tảng Matlab. Thay vì phải giải các hệ phương trình phi tuyến phức tạp, thuật toán điểm trong của SOCP giải quyết bài toán cực tiểu hóa năng lượng tiêu tán dẻo với tốc độ vượt trội.

Chiến lược thích nghi lưới h-adaptivity sử dụng năng lượng tiêu tán dẻo trên từng phần tử làm chỉ số đánh giá hậu sai số rời rạc. Thuật toán phân ngưỡng Dörfler với hệ số tối ưu 0,5 được sử dụng để lọc ra các phần tử có năng lượng tiêu tán lớn nhất, sau đó tiến hành làm mịn tự động thông qua kỹ thuật phân chia đôi qua đỉnh gần nhất (newest vertex bisection). Toàn bộ quy trình tính toán được hoàn thiện và kiểm chứng trong 6 tháng thực hiện đề tài.

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

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

Thứ nhất, mô hình ES-DSG3 loại bỏ hoàn toàn hiện tượng khóa trượt trên toàn bộ dải độ mảnh khảo sát. Khi kiểm tra tấm vuông và chữ nhật với tỷ số nhịp trên chiều dày tăng từ 10 lên 100 (tấm mỏng dần), hệ số tải trọng giới hạn cận trên vẫn hội tụ ổn định với sai số tương đối dưới 0,8% so với nghiệm giải tích, trong khi phần tử hữu hạn tiêu chuẩn bị khóa cứng hoàn toàn.

Thứ hai, kỹ thuật thích nghi lưới h-adaptivity chứng minh hiệu quả phân bố tài nguyên tính toán vượt bậc. Thay vì làm mịn lưới đồng loạt, các phần tử tự động tập hợp dày đặc tại các dải khớp dẻo, khu vực ngàm biên và các điểm kỳ dị góc nhọn (như góc trong của tấm chữ L và góc tù tấm hình thoi 60 độ). Kết quả cho thấy phương pháp thích nghi lưới giúp giảm 45% đến 55% tổng số bậc tự do cần thiết để đạt cùng độ chính xác so với phương pháp chia lưới đều.

Thứ ba, nghiên cứu làm sáng tỏ ảnh hưởng của các điều kiện biên lên khả năng chịu tải dẻo của tấm chữ nhật. Khảo sát trên 4 tổ hợp biên gồm ngàm toàn phần (C-C-C-C), tựa đơn toàn phần (S-S-S-S), hai ngàm - hai tự do (C-F-C-F) và ba ngàm - một tự do (C-C-C-F) cho thấy khi thay đổi tỷ số hai cạnh từ 1,0 đến 2,0, hệ số tải trọng phá hoại suy giảm phi tuyến từ 35% đến 50%.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp đạt độ chính xác cao là nhờ kỹ thuật ES-FEM tạo ra độ mềm phù hợp cho mô hình, bù đắp xu hướng quá cứng của phương pháp phần tử hữu hạn truyền thống, trong khi DSG3 cân bằng trường biến dạng cắt. Khi so sánh với các công bố quốc tế của Vicente da Silva và các nghiên cứu của Bletzinger, độ sai lệch của hệ số tải trọng giới hạn trong luận văn chỉ dao động trong khoảng 0,5% đến 1,2%, khẳng định độ tin cậy tuyệt đối của thuật toán.

Dữ liệu nghiên cứu được biểu diễn trực quan và sinh động thông qua các bảng so sánh sai số tương đối qua 5 vòng lặp làm mịn lưới, biểu đồ đường biểu diễn mối quan hệ giữa hệ số tải trọng giới hạn với tỷ số độ mảnh 2L/h, cùng các biểu đồ phân bố đường mức chuyển vị và đám mây mật độ năng lượng tiêu tán dẻo. Những hình ảnh trực quan này giúp người kỹ sư nhận diện rõ ràng cơ chế phá hoại cục bộ và quỹ đạo phát triển khớp dẻo trong kết cấu tấm.

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

Thứ nhất, tích hợp module tính toán ES-DSG3 vào các phần mềm thương mại: Các doanh nghiệp công nghệ xây dựng cần hợp tác cùng nhóm nghiên cứu để đóng gói thuật toán thành plugin phân tích giới hạn trực tiếp, đặt mục tiêu rút ngắn 30% thời gian phân tích kết cấu tấm vỏ phức tạp trong lộ trình 12 tháng tới.

Thứ hai, chuẩn hóa quy trình phân tích giới hạn thích nghi lưới trong tiêu chuẩn thiết kế: Các cơ quan quản lý và viện nghiên cứu chuyên ngành xây dựng cần đưa quy trình phân tích giới hạn cận trên vào cẩm nang hướng dẫn tính toán cấu kiện sàn đáy hầm, móng bè và sàn nấm, hướng tới mục tiêu tối ưu hóa và tiết kiệm 15% khối lượng vật liệu bê tông cốt thép trong giai đoạn 2024 - 2026.

Thứ ba, mở rộng mô hình cho kết cấu vật liệu nhiều lớp và vật liệu biến tính (FGM): Các nhà khoa học cơ học tính toán cần tiếp tục phát triển công thức ES-DSG3 cho các tấm composite chịu uốn và cắt lớn, hướng tới nâng cao độ chính xác dự báo trạng thái phá hoại lên trên 98% trong vòng 18 tháng.

Thứ tư, đổi mới chương trình đào tạo kỹ sư và thạc sĩ kỹ thuật công trình: Các trường đại học kỹ thuật cần bổ sung chuyên đề tối ưu hóa nón bậc hai (SOCP) và phương pháp phần tử hữu hạn trơn vào học phần Cơ học kết cấu nâng cao, phấn đấu đạt 80% học viên cao học làm chủ kỹ năng lập trình mô phỏng số sau 1 học kỳ.

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

• Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Ứng dụng phương pháp để đánh giá chính xác sức chịu tải giới hạn của bản sàn chịu tải trọng nặng, sàn chuyển, nắp hầm và kết cấu thép tấm, giúp tối ưu hóa tiết diện và đảm bảo an toàn tuyệt đối cho công trình.
• Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng: Khai thác cơ sở lý thuyết hoàn chỉnh về phần tử hữu hạn trơn ES-FEM, kỹ thuật lệch trượt DSG3 và bài toán tối ưu hóa nón bậc hai để phát triển các đề tài nghiên cứu chuyên sâu về cơ học phi tuyến.
• Giảng viên và nhà nghiên cứu tại các trường đại học khối kỹ thuật: Sử dụng luận văn làm tài liệu tham khảo chất lượng cao cho các môn học Phương pháp số trong cơ học công trình, Lý thuyết dẻo và Tối ưu hóa kết cấu.
• Chuyên viên phát triển phần mềm mô phỏng kỹ thuật (CAE/FEA): Tham khảo thuật toán thích nghi lưới tự động h-adaptivity và mã nguồn mô phỏng trên nền tảng Matlab - Mosek để tích hợp vào các phần mềm phân tích kết cấu thế hệ mới.

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

Phương pháp ES-DSG3 giải quyết hiện tượng khóa trượt (shear locking) trong tấm mỏng như thế nào? Trong phần tử tam giác bậc thấp tiêu chuẩn, biến dạng cắt xuất hiện giả tạo khiến độ cứng tăng vọt khi chiều dày tấm giảm. ES-DSG3 sử dụng phép biến đổi rời rạc lệch trượt để chuyển ma trận biến dạng cắt thành hằng số, kết hợp làm trơn biến dạng trên từng cạnh, giúp triệt tiêu hoàn toàn khóa trượt ngay cả khi tỷ số độ mảnh đạt 100.

Tại sao bài toán phân tích giới hạn lại được biến đổi về dạng quy hoạch nón bậc hai (SOCP)? Tiêu chuẩn chảy dẻo von Mises có dạng phương trình bậc hai phi tuyến trong không gian ứng suất. Việc chuyển đổi năng lượng tiêu tán dẻo về ràng buộc nón bậc hai cho phép tận dụng thuật toán điểm trong cực mạnh của phần mềm Mosek, giúp tìm ra nghiệm tối ưu toàn cục chỉ sau 10 đến 15 vòng lặp mà không gặp sự cố phân kỳ.

Chỉ số nào được sử dụng để đánh giá sai số trong kỹ thuật thích nghi lưới của luận văn? Luận văn sử dụng trực tiếp giá trị năng lượng tiêu tán dẻo thu được từ bài toán SOCP trên từng phần tử làm sai số do rời rạc. Cách tiếp cận này loại bỏ các bước tính toán đạo hàm hậu sai số phức tạp, cho phép xác định chính xác các phần tử cần phân chia mịn với hệ số phân ngưỡng 0,5.

Phương pháp có ứng dụng hiệu quả cho các tấm có hình dạng bất kỳ hay không? Phương pháp hoạt động tuyệt vời trên các hình học phức tạp. Luận văn đã kiểm chứng thành công trên 5 cấu hình hình học từ tấm vuông, chữ nhật đến tấm chữ L có góc lõm và tấm tam giác 75 phần tử ban đầu, chứng minh khả năng tự động tập trung lưới chính xác tại các điểm tập trung ứng suất cao.

Lợi ích thực tế lớn nhất của phân tích giới hạn trực tiếp so với phân tích gia tải từng bước là gì? Phân tích gia tải từng bước đòi hỏi giải hàng trăm hệ phương trình đàn dẻo phi tuyến theo thời gian để theo dõi biến dạng, rất tốn kém tài nguyên. Phân tích giới hạn trực tiếp xác định ngay lập tức khả năng chịu tải cực hạn trong 1 lần tối ưu duy nhất, giúp rút ngắn hơn 70% thời gian tính toán.

Kết luận

• Xây dựng thành công quy trình phân tích giới hạn cận trên cho tấm dày Mindlin-Reissner bằng sự kết hợp giữa ES-FEM và DSG3.
• Triệt tiêu hoàn toàn hiện tượng khóa trượt đối với các kết cấu tấm có tỷ số độ mảnh cao đến 100.
• Phát triển kỹ thuật thích nghi lưới h-adaptivity dựa trên năng lượng tiêu tán dẻo, giúp giảm hơn 45% số phần tử tính toán mà vẫn đảm bảo độ hội tụ chính xác.
• Biến đổi bài toán tối ưu dẻo về dạng quy hoạch nón bậc hai (SOCP) và giải quyết hiệu quả bằng công cụ Mosek.
• Kiểm chứng độ chính xác vượt trội qua 5 dạng hình học và 4 điều kiện biên kết cấu khác nhau.

Đóng góp chính của luận văn là đã thiết lập một phương pháp số hoàn chỉnh, ổn định và hiệu quả cao để xác định tải trọng giới hạn phá hoại của kết cấu tấm công trình. Kế hoạch tiếp theo là mở rộng thuật toán sang bài toán tấm vỏ composite nhiều lớp chịu tải trọng động trong vòng 12 đến 24 tháng tới. Hãy đón đầu và ứng dụng ngay phương pháp mô phỏng số ES-DSG3 thích nghi lưới để nâng tầm độ chính xác và tối ưu hóa chi phí cho các dự án xây dựng hiện đại!

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

Chương 1. Tình hình nghiên cứu và tính cấp thiết đề tài. Trên thế giới. Mục tiêu và nhiệm vụ của luận văn:.

Cấu trúc luận văn:.1 Lý thuyết tấm dày Mindlin-Reissner .2 Lý thuyết biến dạng cắt bậc nhất (FSDT) .3 Mô hình vật liệu. Mô hình đàn hồi - dẻo lý tưởng. Mô hình cứng dẻo - lý tưởng .4 Tiêu chuẩn Von Mises .5 Luật chảy dẻo .6 Định đề Drucker – giả thiết về tính ổn định của vật liệu .7 Luật chảy dẻo kết hợp .8 Nguyên lý công ảo. Trường ứng suất khả dĩ tĩnh.

Trường chuyển vị, biến dạng khả dĩ động. PHÂN TÍCH GIỚI HẠN TẤM MINDLIN - REISSNER BẰNG PHƯƠNG PHÁP ES -DSG3 VỚI KỸ THUẬT THÍCH NGHI LƯỚI. Tổng quan về phân tích giới hạn. Công thức phân tích giới hạn tấm dày Mindlin-Reissner.

Phương pháp "rời rạc lệch trượt" – Discrete Shear Gap (DSG3):. Phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM):. Rời rạc phương trình phân tích động (Định lý cận trên). Kỹ thuật tối ưu hình nón bậc hai (Second order cone programming– SOCP):.

Biến đổi bài toán tối ưu về dạng hình nón bậc hai:. Gán các biến về dạng hình nón trong Mosek:. Kỹ thuật thích nghi lưới .1 Tính sai số rời rạc dựa vào năng lượng tiêu tán dẻo.2 Cách thức làm mịn lưới. Sơ đồ khối cho bài toán phân tích giới hạn tấm Mindlin-Reissner bằng phương pháp ES-DSG3 kết hợp với kỹ thuật thích nghi lưới.

Tấm hình vuông. Tấm hình chữ nhật. Tấm hình thoi. Tấm hình tam giác.

KẾT LUẬN VÀ KIẾN NGHỊ. 72 TÀI LIỆU THAM KHẢO. 77 v DANH MỤC HÌNH VẼ Hình 1. Các kiểu thích nghi lưới.

Mô hình Kirchhoff và mô hình Mindlin – Reissner. Lý thuyết biến dạng. Mối liên hệ    trong mô hình đàn - dẻo lý tưởng. Mối liên hệ    trong mô hình cứng - dẻo lý tưởng.

Ứng xử ổn định và không ổn định theo Drucker. Hình học của luật chảy dẻo kết hợp. Vật thể cân bằng. Phần tử tam giác 3 nút.

Phần tử tam giác 3 nút và tọa độ địa phương. Phân chia miền trơn  k. Phân chia phần tử tam giác. Điều kiện biên bài toán tấm hình vuông.

Phân chia miền bài toán tấm hình vuông. Mô phỏng kết quả làm mịn lưới tấm hình vuông biên ngàm. Mô phỏng năng lượng tiêu tán dẻo tấm hình vuông biên ngàm. Mô phỏng trường chuyển vị tấm hình vuông biên ngàm.

Đánh giá hiệu quả sử dụng kỹ thuật thích nghi lưới cho bài toán tấm hình vuông biên ngàm. Hệ số tải trọng giới hạn   tấm hình vuông biên ngàm. Mô phỏng kết quả làm mịn lưới tấm hình vuông biên tựa đơn. Mô phỏng năng lượng tiêu tán dẻo tấm hình vuông biên tựa đơn.

Mô phỏng trường chuyển vị tấm hình vuông biên tựa đơn. Đánh giá hiệu quả sử dụng kỹ thuật thích nghi lưới cho bài toán tấm hình vuông biên tựa đơn. Hệ số tải trọng giới hạn   tấm hình vuông biên tựa đơn. Điều kiện biên cho bài toán tấm hình chữ nhật.

Phân chia miền bài toán tấm hình chữ nhật. Mô phỏng kết quả làm mịn lưới tấm hình chữ nhật biên ngàm. Mô phỏng năng lượng tiêu tán dẻo tấm hình chữ nhật biên ngàm. Mô phỏng trường chuyển vị tấm hình chữ nhật biên ngàm.

Đánh giá hiệu quả sử dụng kỹ thuật thích nghi lưới cho bài toán tấm hình chữ nhật biên ngàm. Hệ số tải trọng giới hạn   tấm hình chữ nhật biên ngàm. Mô phỏng kết quả làm mịn lưới tấm hìnhchữ nhật biên tựa đơn. Mô phỏng năng lượng tiêu tán dẻo tấm hình chữ nhật biên tựa đơn.

Mô phỏng trường chuyển vị tấm hình chữ nhật biên tựa đơn. Đánh giá hiệu quả sử dụng kỹ thuật thích nghi lưới cho bài toán tấm hình chữ nhật biên tựa đơn. Hệ số tải trọng giới hạn   tấm hình chữ nhật biên tựa đơn. Mô phỏng kết quả làm mịn lưới tấm hình chữ nhật biên 2 ngàm – 2 tự do.

Mô phỏng năng lượng tiêu tán dẻo tấm hình chữ nhật 2 ngàm – 2 tự do. Mô phỏng trường chuyển vị tấm hình chữ nhật biên 2 ngàm – 2 tự do. Mô phỏng kết quả làm mịn lưới tấm hình chữ nhật biên 3 ngàm – 1 tự do. Mô phỏng năng lượng tiêu tán dẻo tấm hình chữ nhật 3 ngàm – 1 tự do.

Mô phỏng trường chuyển vị tấm hình chữ nhật biên 3 ngàm – 1 tự do. Quan hệ giữa hệ số   và số phần tử ứng với các điều kiện biên. Quan hệ giữa hệ số   và tỷ số 2H/h ứng với các điều kiện biên. Quan hệ giữa   và tỷ số 2L/2H.

Bài toán tấm hình L. Mô phỏng kết quả làm mịn lưới tấm hình L biên tựa đơn. Mô phỏng năng lượng tiêu tán dẻo tấm hình L biên tựa đơn. Mô phỏng trường chuyển vị tấm hình L biên tựa đơn.

Quan hệ giữa hệ số   và số phần tử tấm hình L biên tựa đơn. Mô phỏng kết quả làm mịn lưới tấm hình L biên ngàm. Mô phỏng năng lượng tiêu tán dẻo tấm hình L biên ngàm. Mô phỏng trường chuyển vị tấm hình L biên ngàm.

Tấm hình thoi   600. Phân chia miền bài toán tấm hình thoi. Mô phỏng kết quả làm mịn lưới tấm hình thoi biên tựa đơn. Mô phỏng năng lượng tiêu tán dẻo tấm hình thoi biên tựa đơn.

Mô phỏng trường chuyển vị tấm hình thoi biên tựa đơn. Mô phỏng kết quả làm mịn lưới tấm hình thoi biên ngàm. Mô phỏng năng lượng tiêu tán dẻo tấm hình thoi biên ngàm. Mô phỏng trường chuyển vị tấm hình thoi biên ngàm.

Quan hệ giữa hệ số   và góc  tấm hình thoi biên tựa đơn. Quan hệ giữa   và góc  tấm hình thoi biên ngàm. Phân chia miền cho tấm hình tam giác với 75 phần tử tam giác ban đầu. Mô phỏng kết quả làm mịn lưới tấm hình tam giác biên tựa đơn.

Mô phỏng năng lượng tiêu tán dẻo tấm hình tam giác biên tựa đơn. Mô phỏng trường chuyển vị tấm hình tam giác biên tựa đơn. Hệ số tải trọng giới hạn   tấm hìnhtam giác biên tựa đơn. 71 ix DANH MỤC BẢNG BIỂU Bảng 4.

Đánh giá tính hiệu quả của phương pháp pháp (2L/h=100) bài toán tấm vuông biên ngàm. Bảng so sánh hệ số tải trọng giới hạn   trong bài toán tấm hình vuông biên ngàm. Đánh giá hiệu quả của phương pháp (2L/h=100) bài toán tấm vuông biên tựa đơn. Bảng so sánh hệ số tải trọng giới hạn   trong bài toán tấm hình vuông biên tựa đơn.

Hệ số tải trọng giới hạn   tấm hình vuông chiều dày tấm giảm dần. Đánh giá hiệu quả của phương pháp (2H/h=100) bài toán tấm hình chữ nhật biên ngàm. Bảng so sánh hệ số tải trọng giới hạn   trong bài toán tấm chữ nhật biên ngàm. Đánh giá hiệu quả của phương pháp (2H/h=100) bài toán tấm hình chữ nhật biên tựa đơn.

Bảng so sánh hệ số tải trọng giới hạn   trong bài toán hình chữ nhật biên tựa đơn. Bảng so sánh hệ số tải trọng giới hạn   trong bài toán tấm hình chữ nhật với các điều kiện biên khác nhau. Hệ số tải trọng giới hạn   cho tấm hình chữ nhật với các điều kiện biên khác nhau. Hệ số tải trọng giới hạn   trong bài toán hình chữ nhật tỷ số hai cạnh thay đổi.

Bảng so sánh hệ số tải trọng giới hạn   tấm hình L biên tựa đơnchịu tải phân bố đều. Hệ số tải trọng giới hạn   tấm hình L biên ngàm chịu tải phân bố đều. Bảng so sánh hệ số tải trọng giới hạn   trong bài toán tấm hình thoi chịu tải phân bố đều. Hệ số tải trọng giới hạn   trong bài toán tấm hình tam giác biên tựa đơnchịu tải phân bố đều.

Tổng quan: Việc phân tích kết cấu đàn dẻo cho đến khi phá hoại là một quá trình phức tạp do phải tiến hành từng bước với việc gia tăng tuần tự tải trọng. Quá trình phân tích trên cho phép ta hiểu biết được toàn bộ quá trình hình thành và phát triển biến dạng dẻo nhưng không có lợi về mặt tính toán số và tiêu tốn không ít thời gian xử lý. Do vậy một hướng khác giúp ta tìm được tải trọng giới hạn khi kết cấu đạt đến trạng thái phá hoại mà không cần phải qua các bước trung gian đó chính là “Lý thuyết phân tích trực tiếp tải trọng giới hạn” hay còn gọi là “Phân tích giới hạn”. Phân tích giới hạn đóng vai trò quan trọng vì đã cung cấp cho người kỹ sư trị số của tải trọng giới hạn làm cơ sở cho việc thiết kế và đánh giá độ an toàn cho kết cấu.

Phương pháp phần tử hữu hạn (Finite Element Method – FEM) là một trong những phương pháp xấp xỉ số mạnh và được sử dụng rộng rãi trong tất cả bài toán nói chung và phân tích giới hạn cho tấm dày nói riêng vì có những ưu điểm sau:  Đơn giản trong các bước thiết lập và giải phương trình.  Giải được hầu hết các bài toán: dao động, ổn định,… với nhiều dạng hình học khác nhau như: tròn, tam giác, tứ giác,….  Áp dụng được nhiều lý thuyết giải khác nhau. Tuy nhiên hạn chế của FEM là hiện tượng “shear locking” sẽ xảy ra đối với phần tử tấm Reissner-Mindlin bậc thấp (low oder) khi chiều dày của tấm mỏng dần về zero do đó dẫn đến kết quả phân tích không chính xác hay không hội tụ.

Một trong những phương pháp số được đề xuất có thể đáp ứng các yêu cầu trên là phương pháp phần tử hữu hạn trơn dựa trên cạnh (Edge-based smoothed finite element method – ES FEM) kết hợp với phương pháp lệch trượt (gọi tắt là ES-DSG) được sử dụng để giải quyết bài toán phân tích giới hạn cho tấm Mindlin. Phương pháp rời rạc lệch trượt (Discrete shear gap – DSG3) được đề xuất bởi tác giả Bletzinger và Bischoff (2000) [26]. Phương pháp “rời rạc lệch trượt” gần giống phương pháp giả định biến dạng tự nhiên ANS (Assumed Natured Strain) và 2 làm việc tốt cho phần tử có hình dáng và bậc khác nhau. Hạn chế của phương pháp này là tính ổn định chưa tốt.

Phương pháp phần tử hữu hạn trơn được đề xuất bởi Giáo sư Liu, phương pháp này mang tính kế thừa của phương pháp phần tử hữu hạn truyền thống, điểm mới của nó là kỹ thuật làm trơn hóa các biến dạng (Strain smoothing technique) mục đích làm cho kết quả chính xác hơn.