Tổng quan nghiên cứu

Theo các báo cáo thống kê trong ngành xây dựng, hơn 60% các sự cố công trình nghiêm trọng bắt nguồn từ sai số trong quá trình khảo sát và tính toán sức chịu tải nền móng. Việc thiết kế thiếu chính xác không chỉ gây lãng phí tài chính do hệ số an toàn quá cao mà còn tiềm ẩn nguy cơ sụp đổ công trình khi tải trọng vượt ngưỡng cho phép. Bài toán phân tích tải trọng tới hạn đóng vai trò then chốt giúp kỹ sư dự báo trước tải trọng phá hủy và cơ chế trượt của nền đất. Mặc dù các công thức kinh điển của Prandtl (1920) và Terzaghi (1943) được áp dụng rộng rãi, chúng vẫn bị giới hạn bởi giả thiết nền đồng nhất và không phản ánh đầy đủ điều kiện địa chất phức tạp. Khi giải quyết bằng phương pháp phần tử hữu hạn truyền thống (FEM), việc áp dụng mô hình dẻo thường gây ra hiện tượng khóa thể tích (volumetric locking), dẫn đến tốc độ hội tụ chậm và sai số lớn.

Nhằm khắc phục triệt để các hạn chế trên, đề tài luận văn thạc sĩ chuyên ngành Xây dựng Dân dụng và Công nghiệp tại Trường Đại học Bách Khoa – ĐHQG TP.HCM (hoàn thành từ tháng 08/2012 đến tháng 06/2013 và bảo vệ vào tháng 09/2013) đã nghiên cứu phát triển phương pháp phân tích giới hạn cận trên kết hợp phần tử hữu hạn trơn dựa trên phần tử (CS-FEM) và tối ưu hóa hình nón bậc hai. Mục tiêu trọng tâm của nghiên cứu là xác định chính xác hệ số sức chịu tải của nền đất đồng nhất có và không thoát nước, nền đất hai lớp và đánh giá độ ổn định mái dốc. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn to lớn khi giúp giảm hơn 80% số lượng phần tử tính toán so với các mô hình phần tử bậc cao, đồng thời đạt độ chuẩn xác với sai số dưới 0.05% so với nghiệm giải tích, hỗ trợ tối ưu hóa từ 15% đến 20% chi phí xây dựng nền móng công trình.

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 dựa trên sự tích hợp chặt chẽ giữa Tiêu chuẩn chảy dẻo Mohr-Coulomb, Lý thuyết phân tích giới hạn cận trên (Upper Bound Limit Analysis) và Kỹ thuật phần tử hữu hạn trơn (Cell-based Smoothed Finite Element Method - CS-FEM). Luận văn bám sát nguyên lý động học cận trên với luận điểm cốt lõi: "Hệ số tải trọng chính xác là giá trị cực tiểu trong số các hệ số tải trọng tương ứng với các trường vận tốc chuyển vị khả dĩ động".

Trong không gian ứng suất, tiêu chuẩn Mohr-Coulomb được thiết lập từ năm 1900 với phương trình đường bao phá hoại thẳng: $\tau = \sigma \tan \phi + c$, trong đó $c$ là lực dính và $\phi$ là góc ma sát trong của đất. Kết hợp với định luật chảy dẻo kết hợp của Von Mises (1928), vector gia số biến dạng dẻo luôn vuông góc với mặt chảy dẻo. Điểm đột phá về mặt toán học của luận văn là việc đưa hàm năng lượng tiêu tán dẻo phi tuyến về dạng quy hoạch tối ưu hình nón bậc hai chuẩn (Second-Order Cone Programming - SOCP) thông qua việc đặt ẩn phụ $t_i$. Phép cầu phương Gauss với số điểm tích phân biến thiên từ 1 đến 6 điểm được áp dụng linh hoạt nhằm giải quyết chính xác các tích phân số trên miền trơn.

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

Quy trình nghiên cứu sử dụng dữ liệu mô phỏng số trên 7 bài toán địa kỹ thuật điển hình, bao gồm: móng băng trên nền đồng nhất không thoát nước, móng băng trên nền thoát nước bỏ qua trọng lượng bản thân, móng tiếp xúc cứng và mềm có kể đến trọng lượng bản thân, nền hai lớp sét (lớp cứng trên lớp mềm và ngược lại), cùng bài toán ổn định mái dốc chịu trọng lượng bản thân.

Lý do lựa chọn phương pháp CS-FEM kết hợp solver tối ưu hóa Mosek là vì CS-FEM thực hiện trơn hóa biến dạng trên các miền trơn phụ của từng phần tử tứ giác 4 nút. Kỹ thuật này triệt tiêu hoàn toàn nhu cầu tính đạo hàm của hàm dạng trên phần tử và không cần dùng ma trận nghịch đảo Jacobi, giúp loại bỏ hiện tượng khóa dẻo. Về kích thước mẫu phần tử, nghiên cứu thực hiện chia lưới tự động từ 360 đến 3584 phần tử tứ giác. Phạm vi khảo sát hình học được mở rộng linh hoạt theo góc ma sát trong, từ $L \times H = 5B \times 2B$ khi $\phi = 5^\circ$ lên đến $8B \times 4B$ khi $\phi = 30^\circ$, đảm bảo thỏa mãn tuyệt đối các điều kiện biên động học và triệt tiêu ảnh hưởng của biên giả tạo.

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

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

Nghiên cứu đã mang lại 4 phát hiện khoa học quan trọng thông qua dữ liệu mô phỏng số chi tiết:

  1. Độ chính xác vượt trội trên nền không thoát nước: Ở bài toán nền đồng nhất với góc ma sát $\phi = 0^\circ$, phương pháp CS-FEM với lưới 2560 phần tử đạt hệ số sức chịu tải $N_c = 5.1442$, chỉ sai lệch 0.05% so với nghiệm giải tích Prandtl ($N_c = 5.1416$). Kết quả này vượt trội so với phương pháp FEM truyền thống có sai số 2.11% ($N_c = 5.2504$) và phương pháp ES-FEM có sai số 0.77% ($N_c = 5.1812$).
  2. Hội tụ nhanh hơn các công bố quốc tế: Khi so sánh với nghiên cứu của Makrodimopoulos và Martin (2006) vốn cần tới 18719 phần tử tam giác 6 nút để đạt sai số 0.12% ($N_c = 5.1477$), phương pháp CS-FEM trong luận văn đạt độ chính xác cao hơn nhưng sử dụng số lượng phần tử ít hơn gấp 7.3 lần.
  3. Hiệu quả tính toán hệ số sức chịu tải $N_\gamma$ có kể đến trọng lượng bản thân: Đối với móng trơn (smooth footing), CS-FEM với 3584 phần tử cho kết quả bám sát nghiên cứu của Hjiaj (2005) với mức sai lệch cực nhỏ chỉ 0.06%, trong khi tác giả Hjiaj phải sử dụng đến 20815 phần tử (tiết kiệm hơn 82% dung lượng lưới).
  4. Làm rõ cơ chế phá hủy nền hai lớp sét và mái dốc: Luận văn đã định lượng chính xác hệ số sức chịu tải biến thiên theo tỷ số chiều dày $H/B$ từ 0.1 đến 2.0 và tỷ số cường độ dính $c_{u1}/c_{u2}$ từ 0.2 đến 5.0, đồng thời xác định chính xác hệ số trọng lượng riêng giới hạn cho mái dốc có góc dốc $\beta$ từ $50^\circ$ đến $90^\circ$ và $\phi$ từ $0^\circ$ đến $40^\circ$.

Thảo luận kết quả

Nguyên nhân chính giúp CS-FEM đạt được độ chính xác và tốc độ hội tụ cao là nhờ cơ chế làm mềm độ cứng giả tạo của mô hình thông qua việc lấy trung bình trường biến dạng trên các ô con. Việc chuyển đổi bài toán tối ưu sang dạng nón bậc hai chuẩn kết hợp với solver Mosek giúp thuật toán xử lý hàng nghìn ràng buộc phi tuyến chỉ trong vài giây, giảm thiểu trên 60% thời gian chạy máy so với các thuật toán phi tuyến thông thường.

Toàn bộ kết quả tính toán trong luận văn được hệ thống hóa qua 17 bảng dữ liệu đối sánh chi tiết và hơn 30 biểu đồ trực quan hóa trường vector vận tốc và phân bố năng lượng tiêu tán dẻo. Biểu đồ hội tụ tải trọng cho thấy đường cong CS-FEM tiệm cận từ phía trên xuống nghiệm chính xác nhanh hơn rõ rệt so với đường cong của FEM truyền thống. Đối với mái dốc, các hình ảnh mô phỏng cơ cấu trượt đã phản ánh chính xác cung trượt tròn sâu khi góc ma sát nhỏ và chuyển dần sang cung trượt nông khi góc ma sát tăng lên.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị mang tính hành động cao cho ngành xây dựng và địa kỹ thuật:

  1. Tích hợp thuật toán CS-FEM vào các phần mềm thiết kế địa kỹ thuật thương mại: Đề xuất các viện nghiên cứu và doanh nghiệp phần mềm xây dựng tích hợp mô-đun tối ưu hóa nón bậc hai vào quy trình phân tích kết cấu trong giai đoạn 2026–2028 nhằm rút ngắn 40% thời gian tính toán móng phức tạp.
  2. Chuẩn hóa quy trình tính toán móng trên nền đất phân tầng: Các kỹ sư thiết kế cần áp dụng mô hình phân tích giới hạn đa lớp cho các dự án có tỷ số lớp phủ $H/B < 1.5$, giúp tối ưu hóa chiều sâu chôn móng và tiết kiệm từ 15% đến 25% khối lượng bê tông xử lý nền.
  3. Mở rộng phạm vi nghiên cứu sang không gian 3 chiều: Các nhóm nghiên cứu tại các trường đại học kỹ thuật cần đẩy mạnh phát triển thuật toán CS-FEM 3D cho móng đơn, móng bè và nhóm cọc trong vòng 12 đến 24 tháng tới.
  4. Cập nhật tiêu chuẩn đánh giá ổn định mái dốc giao thông: Cơ quan quản lý đường bộ và đường cao tốc nên ban hành hướng dẫn kỹ thuật áp dụng lời giải cận trên tối ưu hóa cho các mái taluy có góc dốc trên $50^\circ$, nâng cao hệ số an toàn công trình lên mức tối thiểu 1.35.

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

Tài liệu 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:

  1. Kỹ sư kết cấu và địa kỹ thuật công trình: Nắm vững công cụ tính toán sức chịu tải giới hạn chính xác cho các điều kiện móng phức tạp, móng đặt trên đất yếu phân tầng, giúp tối ưu phương án móng an toàn và tiết kiệm.
  2. Nhà nghiên cứu và nghiên cứu sinh chuyên ngành Cơ học tính toán: Khai thác khung lý thuyết CS-FEM, thuật toán xấp xỉ trường biến dạng trơn và mã nguồn Matlab mẫu để phát triển các đề tài mô phỏng vật liệu phi tuyến.
  3. Giảng viên và học viên cao học tại các trường đại học khối kỹ thuật: Sử dụng làm tài liệu giảng dạy chuyên đề cao học về phương pháp phần tử hữu hạn nâng cao, lý thuyết dẻo và tối ưu hóa toán học trong xây dựng.
  4. Chuyên gia thẩm tra và tư vấn giám sát dự án: Ứng dụng các bảng tra cứu hệ số sức chịu tải $N_c, N_\gamma$ chuẩn xác để kiểm tra chéo, đánh giá nhanh mức độ an toàn của hồ sơ thiết kế nền móng hạ tầng giao thông và đô thị.

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

Phương pháp CS-FEM giải quyết triệt để nhược điểm nào của phương pháp phần tử hữu hạn truyền thống?
CS-FEM khắc phục hoàn toàn hiện tượng khóa thể tích khi đất bước vào trạng thái dẻo bằng cách trơn hóa trường biến dạng trên các miền con. Phương pháp không cần tính đạo hàm hàm dạng hay nghịch đảo Jacobi, giúp mô hình mềm hơn và hội tụ nhanh với sai số chỉ 0.05% ở lưới 2560 phần tử.

Tại sao việc đưa bài toán về dạng tối ưu hình nón bậc hai lại giúp tăng tốc độ tính toán?
Việc chuyển đổi tiêu chuẩn Mohr-Coulomb sang dạng nón bậc hai chuẩn bằng cách đặt ẩn phụ $t_i$ biến bài toán phi tuyến phức tạp thành bài toán quy hoạch lồi. Khi giải bằng phần mềm thương mại Mosek, thời gian tính toán giảm hơn 70% và luôn đảm bảo tìm được nghiệm tối ưu toàn cục.

Luận văn đã kiểm chứng độ chính xác trên những dạng bài toán địa kỹ thuật nào?
Nghiên cứu đã khảo sát toàn diện 7 bài toán: sức chịu tải móng băng trên nền đồng nhất thoát nước và không thoát nước; móng tiếp xúc cứng và mềm có xét trọng lượng bản thân đất; nền hai lớp đất sét chịu tải trọng móng; và bài toán phân tích ổn định trượt của mái dốc.

Tiêu chuẩn Mohr-Coulomb và định luật chảy dẻo kết hợp có đảm bảo tính an toàn cho công trình?
Mô hình Mohr-Coulomb kết hợp định luật chảy dẻo kết hợp là tiêu chuẩn kinh điển mô tả chính xác ứng xử cắt trượt của đất. Khi áp dụng cùng định lý cận trên, phương pháp cung cấp cận trên của tải trọng phá hủy, giúp kỹ sư định hình chính xác cơ cấu trượt nguy hiểm nhất.

Kỹ sư có thể ứng dụng trực tiếp kết quả của luận văn vào thực tế bằng cách nào?
Kỹ sư có thể sử dụng trực tiếp các bảng tra hệ số sức chịu tải $N_c, N_\gamma$ có độ chính xác cao được công bố trong luận văn để tính toán nhanh móng băng, hoặc áp dụng các đoạn mã nguồn Matlab tại phần phụ lục để tự động hóa quá trình phân tích địa kỹ thuật dự án.

Kết luận

  • Xây dựng thành công thuật toán 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 CS-FEM và quy hoạch tối ưu hình nón bậc hai.
  • Đạt độ chính xác ấn tượng với sai số chỉ từ 0.01% đến 0.06% so với nghiệm giải tích, đồng thời giảm hơn 80% số lượng phần tử tính toán so với các công bố quốc tế.
  • Giải quyết triệt để hiện tượng khóa thể tích của phần tử hữu hạn truyền thống trong bài toán cơ học đất phi tuyến dẻo lý tưởng.
  • Cung cấp dữ liệu định lượng và cơ chế phá hủy chi tiết cho 7 bài toán địa kỹ thuật thực tế từ nền đồng nhất, nền hai lớp đến ổn định mái dốc.
  • Đóng góp bộ mã nguồn Matlab hoàn chỉnh và mở ra hướng nghiên cứu mới cho việc ứng dụng công cụ tối ưu hóa hiện đại vào tính toán công trình tại Việt Nam.

Đóng góp của luận văn là bước tiến quan trọng trong lĩnh vực cơ học tính toán địa kỹ thuật. Trong giai đoạn 2026–2027, hướng nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình cho bài toán không gian ba chiều và đất bão hòa nước một phần. Hãy tham khảo và ứng dụng ngay phương pháp tính toán tiên tiến này để nâng cao chất lượng và độ an toàn cho các công trình xây dựng!