Tổng quan nghiên cứu

Trong thực tế thi công công trình, hơn 70% các sự cố hư hỏng nghiêm trọng bắt nguồn từ những sai sót trong quá trình khảo sát và tính toán sức chịu tải của nền móng. Việc dự báo thiếu chính xác tải trọng giới hạn không chỉ dẫn đến nguy cơ sụp đổ kết cấu mà còn gây lãng phí từ 20% đến 35% chi phí đầu tư do các giải pháp thiết kế quá thiên về an toàn. Các công thức giải tích kinh điển của Prandtl đề xuất năm 1920 hay Terzaghi năm 1943 chỉ áp dụng hiệu quả cho nền đất đồng nhất lý tưởng và bộc lộ nhiều hạn chế khi xử lý các bài toán địa kỹ thuật phức tạp như đất nhiều lớp hoặc mái dốc chịu tải.

Nghiên cứu này tập trung giải quyết bài toán xác định tải trọng giới hạn và cơ cấu phá hoại của nền đất thông qua lý thuyết phân tích giới hạn cận trên kết hợp phương pháp phần tử hữu hạn trơn dựa trên ô phần tử (CS-FEM) và công cụ tối ưu hóa hình nón bậc hai (SOCP). Mục tiêu cụ thể là xây dựng quy trình tính toán tự động, chính xác và tối ưu hóa thời gian mô phỏng cho các dạng bài toán: sức chịu tải của nền đồng nhất có hoặc không thoát nước, nền sét hai lớp không đồng nhất và độ ổn định của mái dốc đất.

Được thực hiện trong giai đoạn từ tháng 8 năm 2012 đến tháng 6 năm 2013 tại Thành phố Hồ Chí Minh, công trình mang ý nghĩa khoa học và thực tiễn sâu sắc. Giải pháp mới giúp kiểm soát sai số tính toán tải trọng tới hạn dưới mức 0,05% so với nghiệm giải tích, đồng thời cắt giảm hơn 80% số lượng phần tử cần thiết so với các phương pháp phần tử hữu hạn thông thường, mang lại công cụ đắc lực cho công tác thiết kế địa kỹ thuật hiện đạ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 kết hợp 3 lý thuyết cơ học và toán học mũi nhọn. Đầu tiên là tiêu chuẩn chảy dẻo Mohr-Coulomb với luật chảy dẻo kết hợp, đóng vai trò xác định điều kiện chuyển trạng thái từ đàn hồi sang biến dạng dẻo lý tưởng của vật liệu đất dựa trên hai thông số cơ bản: lực dính kết c và góc ma sát trong phi. Tiếp theo, định lý phân tích giới hạn cận trên được áp dụng để xác định hệ số tải trọng sụp đổ thông qua việc cực tiểu hóa công suất tiêu tán năng lượng dẻo bên trong miền thể tích đất dựa trên trường vận tốc khả dĩ động.

Khái niệm cốt lõi thứ ba là kỹ thuật làm trơn biến dạng trong phương pháp CS-FEM. Khác với phương pháp phần tử hữu hạn truyền thống (FEM) vốn tính biến dạng qua đạo hàm của các hàm dạng, CS-FEM chia mỗi phần tử tứ giác thành các miền trơn phụ và tính biến dạng trung bình thông qua tích phân đường trên biên. Cuối cùng, toàn bộ bài toán tối ưu phi tuyến được chuyển đổi chuẩn xác về dạng quy hoạch nón bậc hai (SOCP) bằng việc đưa vào biến phụ t, cho phép tận dụng các thuật toán điểm trong cực kỳ mạnh mẽ để tìm nghiệm tối ưu toàn cục.

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

Nguồn dữ liệu của nghiên cứu được thu thập và tổng hợp từ hệ thống mô phỏng số trên các bài toán chuẩn quốc tế đã được công bố trong cơ học đất. Cỡ mẫu nghiên cứu bao gồm 7 chuỗi mô hình phân tích với mật độ lưới phần tử tăng dần từ 360, 640, 1000, 1440, 1960, 2560 đến 3584 phần tử tứ giác bậc thấp, tương ứng với các kích thước miền tính toán biến thiên từ 5B x 2B đến 5B x 4B nhằm loại trừ hoàn toàn ảnh hưởng của hiệu ứng biên.

Phương pháp chọn mẫu số liệu được thiết lập dựa trên nguyên tắc đối xứng hình học và điều kiện biên động học chuẩn xác của móng băng và mái dốc. Lý do lựa chọn phương pháp phân tích CS-FEM kết hợp thuật toán tối ưu hóa hình nón bậc hai chạy trên nền tảng solver thương mại Mosek là vì phương pháp này triệt tiêu hoàn toàn hiện tượng khóa thể tích (volumetric locking) thường gặp trong phần tử hữu hạn bậc thấp khi đất chảy dẻo, đồng thời nâng cao tốc độ giải ma trận từ 3 đến 5 lần so với các phương pháp tối ưu hóa phi tuyến thông thường. Toàn bộ chu trình thuật toán được lập trình và xử lý số liệu hoàn chỉnh trên môi trường Matlab trong suốt 10 tháng nghiên cứu.

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

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

Quá trình mô phỏng số đã mang lại 4 phát hiện quan trọng về mặt định lượng:

Thứ nhất, đối với bài toán nền đồng nhất không thoát nước (phi = 0 độ), lời giải CS-FEM với lưới 2560 phần tử đạt giá trị hệ số sức chịu tải Nc = 5,1440. Kết quả này chỉ sai lệch 0,05% so với nghiệm giải tích chính xác của Prandtl là 5,1416. Trong khi đó, phương pháp FEM truyền thống cho sai số lên tới 2,10% (Nc = 5,2478) và phương pháp phần tử hữu hạn trơn dựa trên cạnh ES-FEM sai số 0,77% (Nc = 5,1812).

Thứ hai, trong trường hợp nền đất thoát nước với góc ma sát trong biến thiên từ 5 độ đến 30 độ, CS-FEM cho kết quả trùng khít với đường cong giải tích. Cụ thể tại phi = 30 độ, CS-FEM đạt Nc = 30,1430 (sai số chỉ 0,01%), vượt trội rõ rệt so với kết quả của Sloan và Kleeman đạt Nc = 31,47 (sai số 4,40%) và ES-FEM đạt Nc = 31,02.

Thứ ba, khi phân tích hệ số Ngamma có xét trọng lượng bản thân của đất với góc ma sát từ 5 độ đến 40 độ, mô hình CS-FEM sử dụng 3584 phần tử đạt độ hội tụ tương đương với kết quả của Hjiaj vốn phải sử dụng tới 20815 phần tử tam giác gián đoạn, với mức chênh lệch kết quả cực kỳ nhỏ chỉ 0,06%.

Thứ tư, đối với bài toán nền sét hai lớp và ổn định mái dốc có góc dốc từ 50 độ đến 90 độ, phương pháp đã xác định chính xác sự thay đổi của cơ cấu trượt từ phá hoại cục bộ sang phá hoại tổng thể khi thay đổi tỷ số chiều dày lớp đất H/B và tỷ số sức kháng cắt không thoát nước cu1/cu2.

Thảo luận kết quả

Nguyên nhân căn bản giúp phương pháp CS-FEM đạt được độ chính xác vượt trội nằm ở kỹ thuật trơn hóa trường biến dạng. Bằng cách lấy tích phân trên biên của các miền trơn, độ cứng của mô hình số được làm mềm một cách tự nhiên, giúp triệt tiêu hiện tượng nghẽn dẻo mà không cần phải sử dụng các phần tử bậc cao phức tạp vốn tiêu tốn tài nguyên tính toán.

Trong các báo cáo kỹ thuật, dữ liệu hội tụ của phương pháp nên được trình bày trực quan thông qua đồ thị so sánh sai số tương đối giữa các phương pháp số với nghiệm giải tích theo số lượng phần tử. Bên cạnh đó, cơ chế phá hoại của đất nền được mô tả sinh động qua bản đồ trường vector vận tốc chuyển vị và biểu đồ phân bố mật độ năng lượng tiêu tán dẻo. Các biểu đồ này phản ánh rõ nét vùng nêm đất trượt cứng hình tam giác dưới đáy móng và cung trượt logarit lan tỏa ra hai bên sườn móng, hoàn toàn phù hợp với lý thuyết dẻo kinh điển nhưng có độ phân giải biên sắc nét hơn.

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

Nhằm chuyển hóa các kết quả nghiên cứu thành giá trị thực tiễn trong ngành kỹ thuật xây dựng, 4 giải pháp trọng tâm được khuyến nghị thực hiện:

Thứ nhất, tích hợp thuật toán phân tích giới hạn CS-FEM kết hợp tối ưu hóa hình nón bậc hai vào các phần mềm tính toán địa kỹ thuật chuyên dụng. Mục tiêu là giúp các kỹ sư rút ngắn từ 40% đến 60% thời gian phân tích sức chịu tải của các hệ móng phức tạp so với phương pháp phần tử hữu hạn đàn dẻo phi tuyến từng bước truyền thống. Thời gian triển khai module này ước tính trong vòng 6 đến 12 tháng do các đơn vị phát triển phần mềm xây dựng chủ trì.

Thứ hai, thiết lập quy chuẩn kiểm tra kích thước miền khảo sát biến dạng tự động trong phân tích số. Kỹ sư thiết kế cần mở rộng kích thước hình học mô hình từ 5B x 2B lên mức tối thiểu 5B x 4B khi góc ma sát trong của đất tăng từ 15 độ lên 30 độ trở lên, nhằm đảm bảo sai số biên luôn được kiểm soát dưới ngưỡng 1,0%.

Thứ ba, chuẩn hóa quy trình thu thập dữ liệu cơ lý đất nền đầu vào thông qua các thí nghiệm nén ba trục tiêu chuẩn trong 3 tháng đầu của giai đoạn khảo sát địa chất, đặc biệt chú trọng xác định chuẩn xác tương quan lực dính c và góc ma sát trong phi trong điều kiện có và không thoát nước.

Thứ tư, các viện nghiên cứu và trường đại học chuyên ngành xây dựng cần tiếp tục tài trợ kinh phí để mở rộng hướng nghiên cứu này cho bài toán nền móng không gian 3 chiều và bài toán tương tác động lực học đất - kết cấu trong vòng 12 đến 24 tháng tới.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

Nhóm 1 là các kỹ sư kết cấu và địa kỹ thuật tại các công ty tư vấn thiết kế xây dựng. Luận văn cung cấp phương pháp luận chuẩn xác để đánh giá sức chịu tải giới hạn của móng băng, móng bè đặt trên nền đất yếu nhiều lớp hoặc gần mái dốc, giúp tối ưu hóa kích thước hình học móng và tiết kiệm từ 15% đến 25% khối lượng bê tông cốt thép.

Nhóm 2 là các nhà phát triển phần mềm mô phỏng kỹ thuật công trình (CAE). Tài liệu cung cấp chi tiết thuật toán rời rạc hóa CS-FEM, ma trận biến dạng trơn và cách thiết lập bài toán tối ưu nón bậc hai để lập trình tích hợp vào các phần mềm tính toán địa kỹ thuật thế hệ mới.

Nhóm 3 là giảng viên và nghiên cứu sinh chuyên ngành Cơ học đất - Nền móng và Địa kỹ thuật công trình. Luận văn là tài liệu tham khảo giá trị về ứng dụng lý thuyết dẻo hiện đại, định lý cận trên và phương pháp không lưới cải tiến trong cơ học môi trường liên tục.

Nhóm 4 là sinh viên và học viên cao học ngành Kỹ thuật Xây dựng dân dụng và công nghiệp. Đây là tài liệu mẫu mực về phương pháp nghiên cứu khoa học số, cách thức tổ chức dữ liệu mô phỏng và kỹ năng lập trình giải bài toán cơ học bằng công cụ Matlab và solver tối ưu hóa thương mại.

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

Phương pháp CS-FEM có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn truyền thống trong cơ học đất? Phương pháp CS-FEM sử dụng kỹ thuật trơn hóa biến dạng trên từng ô phần tử, biến trường biến dạng thành hằng số trên mỗi miền trơn. Điều này giúp tính toán trực tiếp mà không cần lấy đạo hàm hàm dạng phức tạp, triệt tiêu hoàn toàn hiện tượng khóa thể tích và giảm sai số tính toán xuống mức 0,05% so với 2,10% của FEM truyền thống ở cùng mật độ lưới 2560 phần tử.

Tại sao nghiên cứu lại lựa chọn phương pháp phân tích giới hạn cận trên thay vì phương pháp cân bằng giới hạn cổ điển? Phương pháp cân bằng giới hạn cổ điển của Terzaghi hay Fellenius đòi hỏi phải giả định trước hình dạng mặt trượt hình tròn hoặc cung phẳng, dễ dẫn đến sai lệch lớn trong đất không đồng nhất. Ngược lại, phân tích giới hạn cận trên tự động tìm kiếm cơ chế phá hoại nguy hiểm nhất thông qua trường vận tốc khả dĩ, đảm bảo tính an toàn cao và độ tin cậy tuyệt đối cho công trình.

Việc chuyển đổi bài toán tối ưu sang dạng hình nón bậc hai (SOCP) mang lại lợi ích gì về mặt tính toán? Hàm mục tiêu tiêu tán năng lượng dẻo theo tiêu chuẩn Mohr-Coulomb có tính chất phi tuyến phức tạp. Việc đưa vào biến phụ t để chuyển đổi về bài toán tối ưu hình nón bậc hai chuẩn cho phép sử dụng các thuật toán điểm trong chuyên dụng của Mosek, giúp tốc độ giải toán nhanh hơn từ 3 đến 5 lần và loại bỏ nguy cơ không hội tụ.

Thuật toán trong luận văn có thể áp dụng cho bài toán nền đất gồm nhiều lớp có tính chất cơ lý khác nhau không? Hoàn toàn có thể. Luận văn đã chứng minh tính hiệu quả vượt trội khi khảo sát nền hai lớp đất sét với tỷ số sức kháng cắt cu1/cu2 biến thiên từ 0,2 đến 5,0. Thuật toán mô tả chính xác hiện tượng mặt trượt bị bẻ gãy khúc hoặc giam hãm cục bộ tại mặt tiếp xúc giữa hai lớp đất mà không gặp lỗi số học.

Độ chính xác của phương pháp khi tính toán hệ số sức chịu tải Ngamma có kể đến trọng lượng bản thân đất đạt mức nào? Phương pháp đạt độ chính xác tương đương các nghiên cứu quốc tế hàng đầu. Khi so sánh với kết quả của Hjiaj công bố năm 2005 sử dụng 20815 phần tử tam giác, mô hình CS-FEM chỉ cần 3584 phần tử tứ giác nhưng vẫn đạt kết quả với độ lệch chỉ 0,06%, giúp tiết kiệm hơn 82% bộ nhớ và thời gian tính toán.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn nổi bật:

  • Xây dựng thành công quy trình kết hợp phương pháp phần tử hữu hạn trơn CS-FEM và lý thuyết phân tích giới hạn cận trên cho bài toán địa kỹ thuật.
  • Thiết lập phương pháp chuyển đổi hàm tiêu tán năng lượng dẻo Mohr-Coulomb về dạng tối ưu hóa hình nón bậc hai chuẩn xác và hiệu quả.
  • Đạt độ chính xác vượt trội với sai số tải trọng tới hạn chỉ 0,05% đối với nền không thoát nước và 0,01% đối với nền thoát nước.
  • Tiết kiệm hơn 80% số lượng phần tử tính toán so với các phương pháp phần tử hữu hạn truyền thống trong bài toán xét trọng lượng bản thân đất.
  • Mô phỏng trực quan và chính xác cơ chế phá hoại trượt dẻo của nền đất hai lớp và mái dốc có góc nghiêng từ 50 độ đến 90 độ.

Đóng góp lớn nhất của công trình là cung cấp một công cụ tính toán số mạnh mẽ, giải quyết triệt để sự xung đột giữa độ chính xác cơ học và khối lượng tính toán trong thiết kế nền móng. Kế hoạch phát triển trong 12 tháng tới tập trung vào việc tự động hóa lưới thích nghi và mở rộng bài toán sang không gian 3D. Các đơn vị tư vấn và kỹ sư thiết kế được khuyến khích ứng dụng ngay thuật toán này để nâng cao độ an toàn và tối ưu hóa chi phí cho các công trình xây dựng.