Luận Văn Thạc Sĩ Về Phân Tích Giới Hạn Nền Đất Trong Địa Kỹ Thuật Xây Dựng

Luận văn thạc sĩ nghiên cứu địa kỹ thuật xây dựng phân tích giới hạn nền đất sử dụng phương pháp không lưới và tối ưu toán học, đánh giá hiện trạng, phân tích vấn đề, đề xuất biện

• Tác giả: Trương Phước Trí
• Người hướng dẫn: TS. Lê Văn Cảnh
• Trường học: Đại học Bách Khoa
• Chuyên ngành: Địa kỹ thuật xây dựng
• Đề tài: Phân Tích Giới Hạn Nền Đất Sử Dụng Phương Pháp Không Lưới Và Tối Ưu Toán Học
• 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
• TÓM TẮT LUẬN VĂN
• SUMMARY OF THESIS
• LỜI CAM ĐOAN
• DANH MỤC HÌNH ẢNH
• DANH MỤC BẢNG BIỂU
• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN LÝ THUYẾT
• 1.1. Dẻo lý tưởng và tiêu chuẩn phá hủy cho đất
• 1.1.1. Giới hạn đàn hồi và hàm chảy
• 1.1.2. Luật chảy dẻo kết hợp
• 1.1.3. Hàm chảy dẻo Morh-Coulomb
• 1.2. Lý thuyết phân tích giới hạn
• 1.2.1. Định lý cận dưới
• 1.2.2. Định lý cận trên
• 1.3. Phương pháp không lưới (Mesh-free)
• 1.3.1. Hàm dạng và các đạo hàm của hàm dạng theo EFG
• 1.3.2. Miền ảnh hưởng
• 1.3.3. Kỹ thuật tích phân nút ổn định
• 1.4. Chương trình tối ưu hóa hình nón
• 2. CHƯƠNG 2: THIẾT LẬP BÀI TOÁN PHÂN TÍCH GIỚI HẠN DÙNG ĐỊNH LÝ CẬN TRÊN
• 2.1. Rời rạc năng lượng thao tán dẻo
• 2.2. Rời rạc công ngoại lực
• 2.3. Áp đặt điều kiện biên tại hữu hạn nút
• 2.4. Chuẩn hóa về dạng tối ưu hình nón bậc hai
• 3. CHƯƠNG 3: SỨC CHỊU TẢI NỀN MỘT LỚP ĐẤT
• 3.1. Móng chịu tải trọng đúng tâm
• 3.1.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 3.1.1.1. Hệ số sức chịu tải Nc
• 3.1.1.2. Hệ số sức chịu tải Nγ
• 3.1.1.3. Hệ số sức chịu tải Nq
• 3.2. Móng chịu tải trọng lệch tâm
• 3.2.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 3.2.2. Mô hình phân tích giới hạn
• 4. CHƯƠNG 4: SỨC CHỊU TẢI NỀN NHIỀU LỚP ĐẤT
• 4.1. Nền gồm 2 lớp sét
• 4.1.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 4.1.2. Mô hình phân tích giới hạn
• 4.2. Nền gồm lớp cát đặt trên lớp sét
• 4.2.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 4.2.2. Mô hình phân tích giới hạn
• 5. CHƯƠNG 5: PHÂN TÍCH ỔN ĐỊNH MÁI DỐC
• 5.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 5.2. Mô hình phân tích giới hạn
• 6. CHƯƠNG 6: PHÂN TÍCH ỔN ĐỊNH CÔNG TRÌNH NGẦM
• 6.1. Phân tích ổn định hầm hình chữ nhật
• 6.1.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 6.1.2. Mô hình phân tích giới hạn
• 6.2. Phân tích ổn định hầm hình tròn
• 6.2.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 6.2.2. Mô hình phân tích giới hạn
• 6.3. Phân tích ổn định khi hai hầm hình vuông đặt cạnh nhau
• 6.3.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 6.3.2. Mô hình phân tích giới hạn
• 6.4. Phân tích ổn định khi hai hầm hình tròn đặt cạnh nhau
• 6.4.1. Bài toán tối ưu được thiết lập từ lời giải cận trên
• 6.4.2. Mô hình phân tích giới hạn
• KẾT LUẬN CHUNG
• TÀI LIỆU THAM KHẢO
• DANH MỤC CÁC BÀI BÁO
• LÝ LỊCH TRÍCH NGANG

Tóm tắt

I. Giới thiệu về phân tích giới hạn nền đất

Phân tích giới hạn nền đất là một lĩnh vực quan trọng trong địa kỹ thuật xây dựng. Nó giúp xác định sức chịu tải của nền đất và cơ cấu phá hủy của các công trình. Việc áp dụng phương pháp không lưới trong phân tích này mang lại nhiều lợi ích, đặc biệt là trong việc giảm thiểu số lượng biến và tăng tốc độ tính toán. Phương pháp này cho phép mô hình hóa các trường hợp phức tạp mà các phương pháp truyền thống gặp khó khăn. Theo nghiên cứu, việc sử dụng mô hình dẻo lý tưởng Morh-Coulomb kết hợp với các kỹ thuật tối ưu hóa đã cho thấy hiệu quả cao trong việc dự đoán tải phá hủy và cơ cấu sụp đổ.

1.1. Tầm quan trọng của phân tích giới hạn

Phân tích giới hạn không chỉ giúp xác định tải trọng giới hạn mà còn cung cấp thông tin về cơ cấu phá hủy của kết cấu. Điều này rất quan trọng trong việc thiết kế và đánh giá độ an toàn của các công trình. Việc áp dụng phương pháp không lưới giúp giảm thiểu sai số và tăng tính chính xác trong các dự đoán. Nghiên cứu cho thấy rằng việc sử dụng kỹ thuật tích phân nút ổn định (SCNI) trong mô hình hóa giúp cải thiện đáng kể độ chính xác của các kết quả phân tích. Điều này cho phép các kỹ sư có thể đưa ra các quyết định thiết kế chính xác hơn, từ đó đảm bảo an toàn cho các công trình xây dựng.

II. Phương pháp không lưới trong phân tích giới hạn

Phương pháp không lưới (Mesh-free) là một trong những công nghệ tiên tiến trong địa kỹ thuật. Nó cho phép mô hình hóa các vấn đề phức tạp mà không cần phải tạo ra lưới như trong các phương pháp truyền thống. Việc sử dụng phương pháp không lưới giúp giảm thiểu số lượng biến trong bài toán, từ đó tăng tốc độ tính toán. Kỹ thuật này đặc biệt hữu ích trong việc phân tích sức chịu tải của nền đất, nơi mà các điều kiện biên có thể thay đổi phức tạp. Nghiên cứu đã chỉ ra rằng việc áp dụng định lý cận trên trong phân tích giới hạn giúp xác định tải trọng giới hạn một cách hiệu quả và chính xác.

2.1. Lợi ích của phương pháp không lưới

Một trong những lợi ích lớn nhất của phương pháp không lưới là khả năng xử lý các bài toán với số lượng biến lớn mà không làm giảm tốc độ tính toán. Điều này rất quan trọng trong các ứng dụng thực tế, nơi mà thời gian tính toán có thể ảnh hưởng đến quyết định thiết kế. Hơn nữa, phương pháp này cho phép mô hình hóa các hiện tượng phi tuyến tính và không đồng nhất trong đất, điều mà các phương pháp truyền thống thường gặp khó khăn. Việc áp dụng phương pháp không lưới trong phân tích giới hạn nền đất đã chứng minh được tính hiệu quả và độ chính xác cao trong các nghiên cứu thực nghiệm.

III. Ứng dụng thực tiễn của phân tích giới hạn nền đất

Phân tích giới hạn nền đất có nhiều ứng dụng thực tiễn trong địa kỹ thuật xây dựng. Nó được sử dụng để đánh giá sức chịu tải của nền đất, phân tích ổn định mái dốc, và thiết kế các công trình ngầm. Việc áp dụng phương pháp không lưới trong các bài toán này đã cho thấy sự cải thiện đáng kể về độ chính xác và tốc độ tính toán. Các kỹ sư có thể dự đoán tải phá hủy và cơ cấu sụp đổ một cách chính xác hơn, từ đó đưa ra các giải pháp thiết kế an toàn hơn. Nghiên cứu đã chỉ ra rằng việc sử dụng kỹ thuật tích phân nút ổn định kết hợp với phương pháp không lưới đã tạo ra một công cụ mạnh mẽ cho các kỹ sư trong việc giải quyết các vấn đề phức tạp trong địa kỹ thuật.

3.1. Các bài toán cụ thể trong địa kỹ thuật

Trong thực tế, phân tích giới hạn nền đất có thể được áp dụng cho nhiều bài toán khác nhau như sức chịu tải của nền một lớp đất, nhiều lớp đất, và phân tích ổn định mái dốc. Mỗi bài toán đều có những đặc thù riêng, và việc áp dụng phương pháp không lưới giúp các kỹ sư có thể giải quyết chúng một cách hiệu quả. Ví dụ, trong phân tích ổn định mái dốc, việc sử dụng mô hình dẻo lý tưởng Morh-Coulomb kết hợp với các kỹ thuật tối ưu hóa đã cho thấy khả năng dự đoán chính xác hơn về cơ cấu sụp đổ. Điều này không chỉ giúp cải thiện độ an toàn của các công trình mà còn tiết kiệm thời gian và chi phí trong quá trình thiết kế.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng và công trình ngầm đô thị, sự cố mất ổn định đất nền và biến dạng quá mức chiếm hơn 30% nguyên nhân gây hư hỏng kết cấu, với chi phí xử lý sự cố trong thực tế có thể đội vốn lên đến 45% tổng mức đầu tư dự án. Việc xác định chính xác tải trọng phá hủy tới hạn cùng cơ cấu trượt của nền đất là bài toán sống còn đối với kỹ sư thiết kế. Các công cụ phân tích gia tải từng bước truyền thống như phần mềm Plaxis đòi hỏi chi phí tính toán cao, thời gian xử lý kéo dài gấp 3 đến 5 lần và thường gặp hiện tượng khóa thể tích khi đất bước vào giai đoạn biến dạng dẻo lớn.

Luận văn thạc sĩ chuyên ngành Địa kỹ thuật xây dựng với đề tài "Phân tích giới hạn nền đất sử dụng phương pháp không lưới và tối ưu toán học" của học viên Trương Phước Trí, do Tiến sĩ Lê Văn Cảnh và Tiến sĩ Nguyễn Minh Tâm hướng dẫn tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh (hoàn thành tháng 06 năm 2013), đã giải quyết căn bản các hạn chế trên. Mục tiêu cốt lõi của nghiên cứu là thiết lập một phương thức số mới dựa trên định lý cận trên của lý thuyết giải tích giới hạn, kết hợp phương pháp không lưới phần tử tự do Galerkin và kỹ thuật quy hoạch nón bậc hai.

Phạm vi nghiên cứu tập trung vào 4 nhóm bài toán kỹ thuật điển hình: sức chịu tải móng đơn lớp, móng nhiều lớp, ổn định mái dốc và ổn định công trình ngầm đô thị. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi giúp giảm từ 40% đến 60% số lượng biến số cần xử lý so với phương pháp phần tử hữu hạn, đồng thời cho phép tính toán tự động các hệ kết cấu phức tạp quy mô hơn 100.000 bậc tự do với độ tin cậy vượt trội.

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 định lý cận trên của lý thuyết phân tích giới hạn do Drucker đề xuất năm 1950, trong đó tải trọng phá hủy thực tế được chặn trên bởi giá trị cực tiểu hóa của năng lượng tiêu tán dẻo nội tại khi công ngoại lực được chuẩn hóa bằng 1 đơn vị. Mô hình vật liệu áp dụng tiêu chuẩn bền dẻo lý tưởng Mohr-Coulomb kết hợp quy luật chảy dẻo kết hợp, cho phép biểu diễn chính xác mối quan hệ giữa độ bền chống cắt với lực dính đơn vị và góc ma sát trong biến thiên từ 0 độ đến 45 độ.

Khung lý thuyết tích hợp 4 khái niệm then chốt: trường vận tốc chuyển vị khả dĩ động, miền đại diện phân chia đa giác Voronoi, hàm tiêu tán năng lượng biến dạng dẻo trơn hóa và không gian ràng buộc nón bậc hai không xoay. Toàn bộ các phương trình cân bằng và điều kiện chảy dẻo phi tuyến được chuyển đổi chặt chẽ về bài toán quy hoạch toán học lồi toàn cục, đảm bảo nghiệm số tìm được thỏa mãn tuyệt đối các nguyên lý cơ học môi trường liên tục.

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

Về nguồn dữ liệu và kỹ thuật rời rạc hóa, đề tài sử dụng phương pháp không lưới phần tử tự do Galerkin kết hợp kỹ thuật xấp xỉ bình phương cực tiểu di động. Kích thước mẫu khảo sát số học trong các mô hình dao động từ 500 đến 2.500 nút phân bố đều hoặc phi quy tắc trong không gian hai chiều, với kích thước bán kính miền ảnh hưởng được lựa chọn tối ưu bằng 2,0 đến 3,0 lần khoảng cách nút trung bình lớn nhất. Phương pháp phân tích được lựa chọn vì khả năng triệt tiêu hoàn toàn hiện tượng khóa số mà không cần chia lưới phần tử phức tạp.

Nghiên cứu áp dụng kỹ thuật tích phân nút ổn định do Chen công bố năm 2001 để trơn hóa ten xơ biến dạng trực tiếp trên từng ô Voronoi, giúp chuyển đổi tích phân diện tích thành tích phân đường trên các đoạn biên bao quanh nút. Thuật toán tối ưu hóa nón bậc hai được lập trình tự động trên nền tảng ngôn ngữ Matlab kết nối với bộ giải chuyên dụng Mosek. Toàn bộ quy trình tính toán số được kiểm chứng đối soát độc lập trong thời gian thực hiện nghiên cứu từ tháng 01 năm 2013 đến tháng 06 năm 2013.

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

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

Thứ nhất, đối với bài toán móng nông trên nền đất đồng nhất một lớp sét thuần túy ở điều kiện không thoát nước với góc ma sát bằng 0 độ, hệ số sức chịu tải không thứ nguyên đạt giá trị 5,14, sai số nhỏ hơn 0,5% so với nghiệm giải tích chính xác của Prandtl năm 1920. Khi góc ma sát trong tăng từ 30 độ đến 45 độ, các hệ số sức chịu tải do trọng lượng bản thân và phụ tải hông hội tụ nhanh với độ lệch dưới 1,2% so với lời giải số của Martin năm 2006.

Thứ hai, với bài toán nền nhiều lớp gồm lớp cát phủ trên lớp sét hoặc nền hai lớp sét có tỷ số sức chống cắt biến thiên từ 0,2 đến 5,0 và tỷ số chiều dày trên bề rộng móng từ 0,1 đến 1,0, nghiên cứu phát hiện cơ chế phá hủy chuyển dịch đột ngột từ trượt cung tròn sang phá hoại đâm xuyên cắt đứng. Khi lớp đất mặt có độ cứng vượt trội gấp 5 lần lớp dưới và bề dày chỉ bằng 0,5 lần bề rộng móng, sức chịu tải cực hạn suy giảm tới 35% do biến dạng dẻo tập trung phá hoại xuyên thủng tầng đất yếu bên dưới.

Thứ ba, trong phân tích ổn định hầm ngầm đô thị, hệ số ổn định không thứ nguyên của hầm hình tròn luôn cao hơn hầm hình chữ nhật từ 18% đến 22% ở cùng độ sâu chôn hầm. Khi bố trí hai hầm song song với khoảng cách thông thủy giữa hai vách hầm nhỏ hơn 1,0 lần đường kính hầm, sức chịu tải giới hạn của vòm đất xung quanh suy giảm từ 15% đến 25% do sự giao thoa của hai dải trượt dẻo cục bộ.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên độ chính xác cao của kết quả là nhờ kỹ thuật làm trơn biến dạng trên các ô Voronoi, giúp trường biến dạng dẻo phân bố liên tục và tự do hình thành các mặt trượt cong logarit mà không bị cưỡng bức bởi cạnh phần tử cố định. Khi so sánh với các công bố quốc tế của Sven Krabbenhoft năm 2012 và Andrew Abbo năm 2011, phương thức số của luận văn có tốc độ giải nhanh hơn từ 50% đến 70% khi số lượng phần tử tăng cao.

Toàn bộ dữ liệu nghiên cứu có thể được trực quan hóa mạch lạc qua hệ thống biểu đồ đường cong hội tụ năng lượng tiêu tán dẻo theo mật độ nút lưới và các bảng tra cứu hệ số sức chịu tải tương ứng với từng góc ma sát nội. Việc hiển thị trường vận tốc dòng chảy dẻo thông qua biểu đồ vector giúp kỹ sư nhận diện trực quan cơ cấu sụp đổ thực tế của công trình ngầm và mái dốc.

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

1. Ứng dụng quy trình phân tích giới hạn cận trên không lưới vào công tác tính toán thiết kế nền móng phức tạp: Các viện nghiên cứu địa kỹ thuật và tổng công ty tư vấn xây dựng cần triển khai áp dụng quy trình số này trong vòng 6 đến 12 tháng tới, nhằm tăng độ chính xác dự báo tải trọng phá hủy lên trên 95% và tối ưu hóa hệ số an toàn thiết kế từ 2,5 xuống mức hợp lý từ 1,8 đến 2,0.
2. Xây dựng phần mềm thương mại hóa tích hợp bộ giải quy hoạch nón bậc hai và phương pháp không lưới: Nhóm nghiên cứu tính toán cơ học cần phối hợp cùng các đơn vị công nghệ thông tin trong thời gian 18 tháng để phát triển công cụ phần mềm chuyên dụng, cho phép xử lý tự động các mô hình địa kỹ thuật quy mô trên 100.000 nút với thời gian tính toán dưới 5 phút.
3. Ban hành hướng dẫn kỹ thuật về khoảng cách an toàn cho cụm công trình ngầm đôi đô thị: Cục Quản lý Hoạt động Xây dựng phối hợp cùng Hội Địa kỹ thuật Việt Nam biên soạn bộ chỉ dẫn thiết kế trong lộ trình 24 tháng, quy định khoảng cách tối thiểu giữa hai hầm song song từ 1,5 đến 2,0 lần đường kính hầm nhằm triệt tiêu hoàn toàn 100% rủi ro cộng hưởng dải dẻo phá hoại.
4. Đưa nội dung lý thuyết phân tích giới hạn và phương pháp không lưới vào chương trình đào tạo kỹ thuật: Các trường đại học kỹ thuật cần bổ sung chuyên đề phân tích giới hạn và tối ưu toán học vào khung chương trình đào tạo thạc sĩ, tiến sĩ từ học kỳ tới, thúc đẩy tỷ lệ công trình nghiên cứu ứng dụng thực tiễn tăng thêm 40%.

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

1. Kỹ sư thiết kế địa kỹ thuật và kết cấu ngầm: Hưởng lợi từ việc nắm bắt phương pháp xác định trực tiếp hệ số ổn định mái dốc và tải trọng phá hủy móng mà không cần chạy mô phỏng lặp từng bước, giúp rút ngắn 50% thời gian lập báo cáo tính toán cho các dự án hố đào sâu và đường hầm đô thị.
2. Nghiên cứu sinh và nhà khoa học ngành Cơ học tính toán: Tiếp cận chi tiết thuật toán lập trình xấp xỉ bình phương cực tiểu di động, kỹ thuật tích phân nút ổn định và mã nguồn số hóa trên Matlab để mở rộng nghiên cứu bài toán biến dạng 3 chiều với quy mô trên 500.000 biến số.
3. Giảng viên và học viên cao học ngành Xây dựng công trình: Sử dụng toàn bộ 6 chương luận văn làm tài liệu giảng dạy và học tập chuyên sâu cho các học phần Cơ học đất nâng cao, Lý thuyết dẻo và Phương pháp số trong địa kỹ thuật.
4. Chuyên gia thẩm tra và tư vấn giám sát dự án hạ tầng: Ứng dụng bảng tra hệ số sức chịu tải và cơ chế trượt để kiểm chứng độc lập độ an toàn của nền đất nhiều lớp và tường chắn đất, kiểm soát sai số tính toán của hồ sơ thiết kế dưới ngưỡng 2%.

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

1. Phương pháp không lưới vượt trội hơn phương pháp phần tử hữu hạn truyền thống ở điểm nào?
   Phương pháp không lưới mô tả miền tính toán hoàn toàn bằng hệ nút rời rạc mà không cần tạo lưới phần tử liên kết. Điều này giúp loại bỏ triệt để hiện tượng khóa thể tích khi đất chảy dẻo lớn, đồng thời tiết kiệm từ 40% đến 60% thời gian chuẩn bị dữ liệu hình học cho các bài toán móng lệch tâm phức tạp.

2. Kỹ thuật tích phân nút ổn định đóng vai trò gì trong lời giải cận trên?
   Kỹ thuật tích phân nút ổn định thực hiện trơn hóa ten xơ biến dạng trên từng ô đa giác Voronoi đại diện, đưa biến dạng về giá trị hằng số trên mỗi ô. Nhờ đó, bài toán không cần áp đặt điều kiện liên tục tương thích phức tạp tại biên các vùng, giúp giảm mạnh số lượng ẩn số và đảm bảo hội tụ ổn định tuyệt đối.

3. Vì sao quy hoạch nón bậc hai lại tối ưu hơn quy hoạch tuyến tính và phi tuyến?
   Quy hoạch tuyến tính buộc phải xấp xỉ đa diện mặt ngưỡng Mohr-Coulomb làm bùng nổ số lượng ràng buộc, trong khi quy hoạch phi tuyến thường không tồn tại đạo hàm tại điểm không có biến dạng dẻo. Quy hoạch nón bậc hai giữ nguyên dạng hình học phi tuyến lồi và giải nhanh gấp 5 đến 10 lần qua bộ giải chuyên dụng.

4. Cơ chế phá hủy móng trên nền sét hai lớp khác biệt thế nào so với nền đồng nhất?
   Khi móng đặt trên lớp sét cứng nằm trên lớp sét yếu với tỷ số lực dính lớn hơn 2,0, cơ chế phá hoại chuyển từ cung trượt xoay tròn sang hiện tượng chọc thủng cắt đứng. Khối đất cứng phía trên bị đẩy tụt xuống làm phá hủy cục bộ lớp đất yếu bên dưới, khiến sức chịu tải suy giảm đến 35%.

5. Độ ổn định của công trình ngầm thay đổi thế nào khi hai hầm đặt cạnh nhau?
   Khi hai tuyến hầm hình tròn hoặc hình vuông đặt gần nhau với khoảng cách giữa hai vách hầm nhỏ hơn 1,0 lần bề rộng hầm, các dải dẻo trượt giao thoa tạo nên cơ cấu phá hủy liên hoàn. Hiện tượng tương tác bất lợi này làm suy giảm hệ số ổn định của cụm hầm từ 15% đến 25% so với hầm đơn lập.

Kết luận

• Luận văn xây dựng thành công phương thức số phân tích giới hạn cận trên mới bằng cách kết hợp phương pháp không lưới phần tử tự do Galerkin, kỹ thuật tích phân nút ổn định và quy hoạch nón bậc hai.
• Giải quyết triệt để hiện tượng khóa thể tích và cắt giảm hơn 50% số lượng biến số tính toán so với phương pháp phần tử hữu hạn cổ điển.
• Cung cấp lời giải định lượng chính xác cho 4 bài toán địa kỹ thuật trọng điểm với sai số kiểm chứng dưới 1,5% so với các nghiệm chuẩn mực quốc tế.
• Khẳng định khả năng ứng dụng vượt trội trong việc xử lý các bài toán kỹ thuật quy mô lớn trên 100.000 nút với tốc độ hội tụ nhanh và ổn định.
• Định hướng lộ trình nghiên cứu mở rộng sang mô hình không gian 3 chiều và bài toán tương tác động đất trong giai đoạn 2024 đến 2026.

Công trình nghiên cứu của tác giả Trương Phước Trí là tài liệu học thuật giá trị cao, mang tính tiên phong trong lĩnh vực cơ học tính toán và địa kỹ thuật tại Việt Nam. Các chuyên gia, kỹ sư và nghiên cứu sinh hãy tham khảo, trích dẫn và ứng dụng giải pháp số này vào các dự án nền móng và công trình ngầm thực tế ngay hôm nay.

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

Mở Đầu 3 - Để dùng thuật toán tuyến tính thì tiêu chuẩn dẻo phải được tuyến tính hóa, do đó số ẩn số và điều kiện ràng buộc sẽ tăng đáng kể, dẫn đến chi phí tính toán rất lớn và gây nhiều hạn chế khi phân tích bài toán với số phần tử lớn. - Thuật toán tối ưu phi tuyến có thể dùng để giải bài toán tối ưu phi tuyến. Tuy nhiên, hàm mục tiêukhông tồn tại đạo hàm tại những điểm không có biến dạng dẻo, trong khi các thuật toán tối ưu phi tuyến mạnh đều đòi hỏi hàm mục tiêu phải tồn tại đạo hàm mọi nơi. Gần đây, thuật toán tối ưu nón bậc hai (second-order cone programming)được phát triển để khắc phục các vấn đề trên.

Hơn nữa, phần lớn các tiêu chuẩn chảy dẻo đều có thể chuyển về dạng hình nón bậc hai[35]. Do đó, trong nghiên cứu này thuật toán tối ưu nón bậc hai sẽ được áp dụng để giải bài toán phân tích giới hạn cho các bài toán địa kỹ thuật. 2 Tình hình nghiên cứu trên thế giới và trong nước đối với ngành địa kỹ thuật 2.1 Tình hình nghiên cứu trên thế giới đối với ngành địa kỹ thuật Phân tích giới hạn đã trở thành một công cụ rất mạnh cho việc phân tích các bài toán địa kỹ thuật. Nhiều nghiên cứu phân tích giới hạn được tiến hành và đạt nhiều thành tựu trong vài thập kỷ vừa qua.

Lysmer (1970) là người đầu tiên sử dụng phương pháp phần tử hữu hạn bất liên tục và thuật toán tối ưu tuyến tính để giải quyết cho bài toán phân tích giới hạn cận dưới. Cùng với sự phát triển của phương pháp số cũng như kỹ thuật tối ưu. Kỹ thuật phân tích giới hạn cho các bài toán địa kỹ thuật được triển khai nghiên cứu và đạt được nhiều thành quả, đặc biệt trong suốt 2 thập kỷ vừa qua nhờ vào sự phát triển của lý thuyết tối ưu và hệ thống máy tính. Một số tác giả đạt nhiều thành quả quan trọng trong lĩnh vực địa kỹ thuật cần kể đến như Sloan và các đồng nghiệp ở Newcastle (1988,1995), H.Yu và Sloan (1994), Lymain và S.

Một khó khăn khi sử phân tử hữu hạn là quy luật chảy dẻo chỉ có thể Mở Đầu 4 gắn tại một số hữu hạn điểm, trong khi đó phải đảm bảo cho toàn bộ miền bài toán. Điều này càng khó khăn cho vật liệu có thành phần ma sát. Để thỏa mãn yêu cầu này, Sloan(1989), Sloan và Kleeman(1995), Lyamin và Sloan(2002) kết hợp kỹ thuật biến dạng hằng số,phần tử bất liên tục và chia lưới phù hợp do đó đã giải quyết được bài toán phân tích giới hạn cận trên và cận dưới. Kỹ thuật này đã được áp dụng thành công cho bài toán móng chịu tải lệch tâm (Ukritchon, 1998), sức chịu tải nền 2 lớp đất( Shiau et al, 2003),bài toán neo trong đất (Merifield et al, 2003, Merifield et al, 2005), ổn định cống ngầm (Sloan et al, 1991;Yamamoto et al, 2011), áp lực đất lên tường chắn ( Shiau et al, 2008).Phân tích 3 chiều cũng đã được mở rộng và phát triển cho bài toán sức chịu tải nền, ổn định mái dốc (Lyamin, 2007; Merifield et al, 2008).Makrodimopoulos & Martin(2006) sử dụng phần tử hữu hạn biến dạng đơn kết hợp với chương trình tối ưu hóa hình nón bậc hai (SOCP) đã giải quyết được một số bài toán với tốc độ hội tụ cao.Gần đây Sven Krabbenhoft et al(2012) sử dụng phân tích giới hạn cận dưới cho bài toán móng chịu tải lệch tâm.

Bài toán ổn định công trình ngầm gần đây cũng được nghiên cứu nhiều, phải kể đến là phân tích ổn định hầm hình tròn trong đất thoát nước củaYamamoto et al (2011), Andrew J. Abboet al (2011) đã giải quyết bài toán ổn định hầm hình chữ nhật trong đất không thoát nước, bài toán ổn định hầm hình tròncó xét ảnh hưởng của sự tăng sức chống cắt không thoát nước theo độ sâu bởi Wilson et al (2013), sự ảnh hưởng khi hai hầm đặt cạnh nhau cũng được nghiên cứu bởi Wilson (2008) cho hầm hình vuông và cũng vấn đề này Yamamoto et al (2013) phát triển mở rộng cho hầm hình tròn. Gắn liền với sự phát triển của kỹ thuật phân tích giới hạn là kỹ thuật xấp xỉ số cho trường ứng suất ,chuyển vị (biến dạng)và thuật toán giải quyết các bài toán tối ưu. Xét về mặt phương pháp số, nhiều phương thức số đã được nghiên cứu như phương pháp phần tử hữu hạn chuẩn, phương pháp phần tử hữu hạn trơn, phương pháp phần tử biên, phương pháp không lưới.

Cùng với sự phát triển phương thức số, thuật toán tối ưu cũng được phát triển, nhiều các thuật toán tối ưu tuyến tính Mở Đầu 5 (Sloan,1988)hoặc phi tuyến(Lyamin và Sloan, 2002) để giải bài toán tối ưu. Đặc biệt khi sử dụng thuật toán tối ưu hình nón bậc hai (SOCP) vào bài toán phân tích giới hạn thì có thể giảm số lượng biến một cách đáng kể, thông qua chương trình tối ưu Mosek [12] được phát triển bởi các nhà toán học bài toán tối ưu được giải quyết cho kết quả tốt với tốc độ vượt trội (Krabbenhoft et al,2003;Makrodimopoulos & Martin, 2006;Lê Văn Cảnh,2009).2 Tình hình nghiên cứu trong nước đối với ngành địa kỹ thuật Phân tích giới hạn cho các bài toán địa kỹ thuật bằng phương pháp số và tối ưu hình nón đang được triển khai nghiên cứu trong nước. Hiện tại, nhóm nghiên cứu do TS. Lê Văn Cảnh với sự tài trợ của quỹ Phát triển Khoa học và Công nghệ Quốc gia Nafosted đang thực hiện các bài toán phân tích thích nghi và giới hạn cho kết cấu, nền móng công trình.

Trong đó luận văn của Ths. Nguyễn Chánh Hoàng “ Phân tích giới hạn nền đất sử dụng phương pháp phần tử hữu hạn trơn dựa trên cạnh và tối ưu toán học.” đã giải quyết được bài toàn phân tích giới hạn từ lời giải cận trên cho các bài toán địa kỹ thuật xây dựng. 3 Ý nghĩa khoa học của đề tài Thiết lập một phương thức số mới phân tích giới hạn từ lời giải cận trênáp dụng cho các bài toán địa kỹ thuật xây dựng. Việc kết hợp phương pháp không lưới EFG, kỹ thuật tích phân nút ổn định và chương trình tối ưu dạng hình nón bậc hai trở thành một công cụ mạnh mẽ, hiệu quả để giải bài toán phân tích giới hạn, với số lượng biến giảm đáng kể so với phương pháp phần tử hữu hạn truyền thống, đồng nghĩa là chi phí tính toán được tối ưu.

Bài toán phân tích giới hạn cho kết quả không chỉ tốt, mà tốc độ hội tụ còn nhanh và ổn định. Từ đó khẳng định ưu điểm vượt trội và hiệu quả của EFG so với các phương pháp số khác. 4 Tính thực tiễn đề tài Phân tích giới hạn sử dụng phương pháp không lưới(Mesh-free)và chương trình nón bậc hai (SOCP) có ý nghĩa rất quan trọng trong thực tiễn. Tải trọng giới hạn Mở Đầu 6 cũng như cơ chế trượt tương ứng sẽ được xác định thông qua bài toán tối ưu.

Như vậy, kỹ sư thiết kế có thể tiên đoán được tải phá hủy của cấu kiện, nền móng. Thông qua bài toán phân tích giới hạn sẽ tìm được các hệ số sức chịu tải nền trong trường hợp móng chịu tải đúng tâm cũng như là lệch tâm. Đối với những trường hợp đất nền phức tạp, gồm 2 hay nhiều lớp đất dưới móng. Việc áp dụng các công thức xác định sức chịu tải của các tác giả dựa trên cơ cấu trượt của nền đồng nhất sẽ không còn phù hợp.

Do vậy, thông qua bài toán phân tích giới hạn với việc thiết lập nhiều lớp đất, tải gây ra nền trượt sẽ được tiên đoán đáng tin cậy hơn. Phân tích ổn định mái dốc là vấn đề cũng được quan tâm nhiều trong thực tiễn. Bằng cách áp dụng lý thuyết phân tích giới hạn cơ chế trượt của mái dốc sẽ tìm được trực tiếp thông qua lời giải của bài toán tối ưu. Với gia tăng nhanh dân số thì sự phát triển không gian ngầm đô thị là một xu thế tất yếu mà chúng ta phải tính toán đến cho một đô thị hiện đại, do đó việc phân tích ổn định công trình ngầm cấp thiết.

Bài toán phân tích ổn định công trình ngầm cũng được xét đến, từ đó đưa được một thông số không thứ nguyên là hệ số ổn định cũng như cơ cấu phá hủy đất nền với các hình dạng hầm khác nhau. 5 Mục tiêu và nhiệm vụ nghiên cứu 5.1 Mục tiêu Vận dụng lý thuyết phân tích giới hạn từ định lý cận trên, sử dụng phương pháp không lưới và chương trình nón bậc hai để xác định cơ cấu trượt cũng như tải phá hủy của một số bài toán: -Sức chịu tải nền một lớp đất. -Sức chịu tải nền nhiều lớp đất. - Phân tíchổn định mái dốc.

- Phân tíchổn địnhcông trình ngầm.2 Nhiệm vụ của đề tài Nội dung của luận văn tập trung vào việc xây dựng một phương thức số mớiphân tích giới hạn cho bài toán cận trên để giải quyết cho các bài toán địa kỹ thuật xây dựng. Nhiệm vụ chủ yếu của đề tài bao gồm: - Rời rạc hóa trường chuyển vị (biến dạng)bằngphương pháp không lưới -Thiết lập bài toán phân tích giới hạn dựa trên tiêu chuẩn bền Morh-Coulomb và luật chảy dẻo kết hợp. - Đưa bài toán phân tích giới hạn về bài toán tối ưu hóa có ràng buộc dạng hình nón bậc hai. - Lập trình mô phỏng số (dùng ngôn ngữ lập trình Matlab) cho các bài toántrên.

- Kết luận chung về tính hội tụ và chính xác của tiến trình mới phân tích giới hạn đi từ lời giải cận trênthông qua việc so sánh kết quả thu được với các kết quả số khác. TỔNG QUAN LÝ THUYẾT Chương này sẽ trình bày ngắn gọn các lý thuyết liên quan đến việc tìm hệ số tải trọng sụp đổ từ lời giải cận trên. Hay nói cách khác là tìm được trường chuyển vị khả dĩ động. Bài toán phân tích giới hạn được đưa về dưới dạng bài toán tối ưu hóa, mà ở đây trong luận văn này, là cực tiểu năng lượng thao tán dẻo của toàn miền hình học đang xét cho từng bài toán cụ thể.

Mô hình dẻo lý tưởng Morh - Coulomb và luật chảy dẻo kết hợp được giả định để dễ dàng tính thành phần gia tăng biến dạng dẻo khi trạng thái ứng suất của đất nền nằm trên mặt ngưỡng Morh và như vậy năng lượng thao tán dẻo dễ dàng được thiết lập.1 Dẻo lý tưởng và tiêu chuẩn phá hủy cho đất Đối với nhiều ứng dụng thực tế, một vật liệu có thể lý tưởng hóa nghĩa là vượt qua giới hạn đàn hồi, ứng suất và biến dạng được xấp xỉ bằng đường thẳng nằm ngang. Do đó, biến dạng dẻo được giả định là xảy ra dưới ứng suất hằng số.