Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình, việc xác định chính xác sức chịu tải của nền đất đóng vai trò quyết định đến độ an toàn kết cấu và chi phí đầu tư. Theo thống kê thực tế trong ngành địa kỹ thuật, các tính toán lý thuyết truyền thống thường có độ lệch an toàn từ 15% đến 30% so với điều kiện làm việc thực tế, trong khi các mô hình số nếu xử lý điểm kỳ dị không chuẩn xác có thể dẫn đến sai số biên độ lớn từ 10% đến 25%. Vấn đề cốt lõi đặt ra là mô hình vật liệu đàn dẻo Mohr-Coulomb tiêu chuẩn có bề mặt chảy dẻo dạng hình chóp lục giác trong không gian ứng suất chính, tại các vị trí cạnh và đỉnh xuất hiện các điểm dị thường không trơn khiến vector pháp tuyến dòng dẻo không xác định, gây hiện tượng phân kỳ toán học khi phân tích phần tử hữu hạn.

Luận văn thạc sĩ chuyên ngành Địa kỹ thuật xây dựng (mã số 8580211) thực hiện tại Trường Đại học Bách khoa – ĐHQG TP.HCM từ tháng 02/2022 đến tháng 06/2022 tập trung giải quyết bài toán phi tuyến vật liệu này. Mục tiêu nghiên cứu cụ thể bao gồm: hệ thống hóa lý thuyết sức chịu tải và cơ sở toán học của mô hình Mohr-Coulomb; đề xuất thuật toán hồi quy ứng suất kết hợp công thức ma trận mô đun tiếp tuyến cải tiến; xây dựng hoàn chỉnh chương trình máy tính trên ngôn ngữ MATLAB code; kiểm chứng độ chính xác qua bài toán thí nghiệm ba trục và bài toán móng băng chịu tải phân bố trong hai trường hợp móng cứng và móng mềm.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc làm chủ thuật toán số lõi, rút ngắn thời gian xử lý phương trình phi tuyến từ 20% đến 35% so với phương pháp cổ điển, đồng thời đảm bảo kết quả tính toán sức chịu tải cực hạn có độ tương đồng cao với phần mềm Plaxis Connect Edition V20 với sai số kiểm chứng dưới 3.8%.

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 cơ học vật rắn biến dạng và cơ học đất nâng cao, kết hợp hai khung lý thuyết trọng tâm:

  1. Lý thuyết đàn dẻo Mohr-Coulomb và định luật chảy dẻo Koiter (công bố năm 1953): Bề mặt chảy dẻo được định nghĩa thông qua các thông số cơ học cơ bản gồm lực dính c và góc ma sát trong phi. Khi ứng suất chạm đến trạng thái giới hạn dẻo, định luật Koiter cho phép thiết lập quan hệ gia số biến dạng dẻo trên đa bề mặt chảy.

  2. Lý thuyết sức chịu tải cổ điển Terzaghi (năm 1943): Cung cấp công thức giải tích xác định sức chịu tải cực hạn qua bộ ba hệ số Nc, Nq và Ngamma, đóng vai trò làm giá trị chuẩn đối chiếu.

Các khái niệm then chốt được phát triển trong mô hình gồm: bề mặt chảy dẻo đa diện không trơn (non-smooth yield surface), điểm dị thường (singularity point), thuật toán ánh xạ hồi quy lặp ẩn (implicit return mapping) và ma trận mô đun tiếp tuyến nhất quán cải tiến (consistent tangent modulus matrix).

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

Nghiên cứu sử dụng phương pháp mô phỏng số dựa trên phương pháp phần tử hữu hạn (FEM) phi tuyến kết hợp thuật toán lặp ẩn Newton-Raphson:

  • Nguồn dữ liệu và cỡ mẫu mô phỏng: Dữ liệu đầu vào thu thập từ các chỉ tiêu cơ lý tiêu chuẩn qua thí nghiệm cắt trực tiếp và nén ba trục. Cỡ mẫu số học được thiết lập trên miền tính toán bán không gian với kích thước chiều rộng 20 m đến 30 m và chiều sâu 10 m đến 15 m. Toàn bộ miền được rời rạc hóa thành 1247 nút và hơn 400 phần tử hữu hạn bậc cao.

  • Phương pháp chọn mẫu và chia lưới phần tử: Sử dụng phần tử tứ giác 8 nút đẳng tham số (Q8) kết hợp sơ đồ tích phân Gauss 2x2 và 3x3 điểm. Mật độ lưới được chia mịn cục bộ tại khu vực chân móng với bán kính ảnh hưởng từ 1.5B đến 3.0B (với B là bề rộng móng), nơi tập trung ứng suất cắt và gradient biến dạng dẻo lớn nhất.

  • Lý do lựa chọn phương pháp phân tích: Thuật toán return mapping lặp ẩn theo sơ đồ backward-Euler của de Souza Neto (2008) kết hợp phương pháp hình học của Clausen (2007) được lựa chọn vì tính ổn định số học vượt trội, bảo đảm hội tụ vô điều kiện, loại bỏ hiện tượng nhảy bước gradient và không bị suy biến khi góc Lode tiến sát giới hạn 29 độ đến 30 độ.

  • Tiến độ nghiên cứu: Toàn bộ quá trình từ tổng quan tài liệu, lập trình MATLAB, chạy thử nghiệm kiểm chứng và so sánh đối soát được hoàn thành trong chu kỳ 5 tháng nghiên cứu liên tục.

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

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

Quá trình chạy mô phỏng số và phân tích so sánh đã đưa ra 4 phát hiện quan trọng:

  1. Thuật toán return mapping lặp ẩn xử lý chính xác 4 trạng thái hồi quy ứng suất: hồi quy về một mặt phẳng, góc trái, góc phải và đỉnh chóp kim tự tháp Mohr-Coulomb, triệt tiêu hoàn toàn lỗi phân kỳ tại các điểm dị thường.

  2. Kiểm chứng bài toán mô phỏng nén ba trục theo số liệu chuẩn của Karaoulanis cho thấy biểu đồ tải trọng - chuyển vị của thuật toán đạt độ chính xác gần như tuyệt đối, sai số so với nghiệm giải tích dưới 1.2% ở cả trạng thái kéo và nén.

  3. Trong bài toán móng băng, hệ số tăng tải hội tụ ổn định đạt mức lambda = 234 đối với trường hợp móng mềm và lambda = 355 đối với trường hợp móng cứng, phản ánh rõ ràng sự phát triển của các vùng trượt dẻo cục bộ lan rộng thành mặt trượt tổng quát.

  4. So sánh sức chịu tải cực hạn tính toán bằng MATLAB code với phần mềm chuyên dụng Plaxis Connect Edition V20 và công thức giải tích Terzaghi cho thấy độ lệch kết quả chỉ dao động từ 2.1% đến 4.2%, chứng minh độ tin cậy cao của bộ mã nguồn tự phát triển.

Thảo luận kết quả

Sự tương đồng cao giữa thuật toán tự lập trình và phần mềm thương mại xuất phát từ việc áp dụng ma trận mô đun tiếp tuyến nhất quán cải tiến theo phương pháp hình học, giúp tốc độ hội tụ bậc hai của thuật toán Newton-Raphson được duy trì ổn định. Trung bình mỗi bước gia tải chỉ mất từ 4 đến 6 vòng lặp để đạt dung sai sai số lực dư nhỏ hơn 0.001.

Dữ liệu mô phỏng trong nghiên cứu được trực quan hóa rõ nét thông qua hệ thống biểu đồ và bảng biểu chi tiết:

  • Biểu đồ quan hệ ứng suất - chuyển vị tại các nút đặc trưng (nút 1010, nút 1224 và nút 1247) mô tả trực quan 3 giai đoạn ứng xử: giai đoạn đàn hồi tuyến tính ban đầu, giai đoạn chuyển tiếp dẻo cục bộ và giai đoạn dẻo chảy hoàn toàn khi đường cong chuyển vị tiệm cận phương ngang.

  • Các đồ thị đẳng trị ứng suất chính (sigma x, sigma y, sigma z, tau xy) và trường dịch chuyển tổng thể thể hiện rõ nét cơ chế phá hoại cắt trượt dạng cung trượt Prandtl dưới đáy móng.

  • Bảng tổng hợp so sánh thời gian xử lý chứng minh rằng việc áp dụng công thức ma trận cải tiến giúp giảm 18% đến 25% thời gian tính toán của CPU so với các thuật toán return mapping cổ điển. Kết quả này vượt trội hơn so với các phương pháp làm trơn gần đúng góc Lode vốn thường làm biến dạng mặt phá hoại thực tế.

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

Dựa trên kết quả nghiên cứu, 4 khuyến nghị hành động cụ thể được đề xuất nhằm nâng cao hiệu quả ứng dụng thực tiễn:

  1. Tích hợp module thuật toán return mapping cải tiến vào các công cụ tính toán địa kỹ thuật nội bộ của các đơn vị tư vấn thiết kế.
  • Mục tiêu: Giảm thiểu 25% thời gian phân tích kết cấu móng phức tạp.
  • Thời gian: Thực hiện trong lộ trình 6 đến 12 tháng.
  • Chủ thể: Đội ngũ kỹ sư R&D và kỹ sư địa kỹ thuật công trình.
  1. Mở rộng mã nguồn MATLAB từ mô hình Mohr-Coulomb sang các mô hình đất phi tuyến nâng cao như Hardening Soil Model và Cam-Clay cải tiến.
  • Mục tiêu: Tăng độ chính xác dự báo lún cố kết dài hạn đạt trên 92%.
  • Thời gian: Lộ trình nghiên cứu 12 đến 18 tháng.
  • Chủ thể: Các nhóm nghiên cứu chuyên sâu tại các trường đại học và viện chuyên ngành.
  1. Chuẩn hóa quy trình tự động phân tích độ ổn định công trình hố đào sâu và mái dốc có xét đến đa lớp đất.
  • Mục tiêu: Kiểm soát hệ số an toàn ổn định công trình trong ngưỡng an toàn FS từ 1.3 đến 1.5.
  • Thời gian: Áp dụng thử nghiệm trong vòng 6 tháng.
  • Chủ thể: Đơn vị tư vấn thẩm tra và giám sát địa kỹ thuật.
  1. Ứng dụng công nghệ tính toán song song trên nền tảng GPU để giải các bài toán không gian 3 chiều quy mô lớn.
  • Mục tiêu: Tối ưu hóa hiệu năng, xử lý các lưới phần tử có quy mô trên 10.000 phần tử trong thời gian dưới 15 phút.
  • Thời gian: Kế hoạch triển khai từ 9 đến 15 tháng.
  • Chủ thể: Các chuyên gia phát triển phần mềm mô phỏng cơ học tính toán.

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

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

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Địa kỹ thuật - Công trình xây dựng: Nắm bắt phương pháp luận nghiên cứu, kỹ thuật lập trình phần tử hữu hạn phi tuyến từ gốc và sử dụng mã nguồn MATLAB làm tài liệu đối chiếu cho các đề tài phát triển thuật toán số.

  2. Kỹ sư thiết kế nền móng và công trình ngầm tại các doanh nghiệp tư vấn: Hiểu rõ bản chất cơ học và thuật toán hồi quy ứng suất bên trong các phần mềm thương mại như Plaxis hay GeoStudio, từ đó thiết lập các thông số mô hình chính xác và tránh các lỗi phân kỳ số.

  3. Giảng viên và nhà khoa học trong lĩnh vực Cơ học đất - Nền móng: Sử dụng các ví dụ kiểm chứng số, thuật toán ma trận và sơ đồ khối trong luận văn làm học liệu giảng dạy cho các môn học Phương pháp phần tử hữu hạn nâng cao.

  4. Lập trình viên và kỹ sư phát triển phần mềm CAE/CAD: Tiếp cận cấu trúc mã nguồn giải phương trình phi tuyến đàn dẻo, phương pháp xây dựng ma trận tiếp tuyến nhất quán phục vụ việc phát triển các phần mềm kỹ thuật nội địa hóa.

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

Tại sao mô hình Mohr-Coulomb lại gây ra lỗi phân kỳ toán học tại các góc và đỉnh trong phân tích phần tử hữu hạn? Mặt năng suất Mohr-Coulomb trong không gian ứng suất chính có hình chóp lục giác với 6 cạnh và 1 đỉnh nhọn. Tại các vị trí biên này, gradient đạo hàm riêng không liên tục khiến vector pháp tuyến dòng dẻo không xác định được duy nhất, tạo ra các bước nhảy gradient gây lỗi chia cho 0 hoặc vòng lặp lặp lại vô hạn trong thuật toán giải số truyền thống.

Thuật toán return mapping lặp ẩn vượt trội hơn phương pháp lặp hiện như thế nào? Phương pháp lặp hiện (forward-Euler) chỉ ước lượng ứng suất theo hướng tiếp tuyến ban đầu nên dễ bị lệch xa khỏi mặt chảy dẻo khi bước tải lớn, đòi hỏi chia bước cực nhỏ từ 500 đến 1000 bước. Ngược lại, thuật toán lặp ẩn (backward-Euler) chiếu vuông góc ứng suất thử về mặt dẻo thực tế, đảm bảo tính hội tụ ổn định với chỉ 4 đến 6 bước lặp mỗi gia số tải.

Chương trình MATLAB code trong luận văn có thể thay thế phần mềm Plaxis trong thực tế không? Mã nguồn MATLAB được thiết kế hoàn chỉnh cho bài toán móng băng 2D phẳng nhằm mục đích nghiên cứu học thuật và làm chủ thuật toán lõi. Dù độ lệch kết quả so với Plaxis chỉ dưới 3.8%, để thay thế thương mại cần bổ sung thêm giao diện đồ họa, module tương tác đất - kết cấu và các bài toán 3D phức tạp.

Độ tin cậy của mô hình móng băng trong nghiên cứu này so với công thức Terzaghi cổ điển ra sao? Kết quả tính toán từ thuật toán đề xuất cho sức chịu tải cực hạn móng băng khớp với công thức Terzaghi với sai số biên độ từ 2.1% đến 4.5%. Điểm ưu việt của mô hình phần tử hữu hạn là mô phỏng được cả ứng xử chuyển vị phi tuyến và phân bố ứng suất cục bộ dưới đáy móng mà công thức Terzaghi không thể hiện được.

Thuật toán trong luận văn có thể mở rộng cho các bài toán hố đào sâu và ổn định mái dốc không? Hoàn toàn khả thi. Thuật toán return mapping và ma trận tiếp tuyến cải tiến là giải thuật cục bộ tại từng điểm tích phân Gauss, do đó có thể tích hợp trực tiếp vào bài toán phân tích ổn định hố đào sâu, tường chắn đất hay mái dốc bằng phương pháp giảm độ bền chống cắt (Strength Reduction Method).

Kết luận

  • Xây dựng thành công thuật toán return mapping lặp ẩn giải quyết triệt để vấn đề kỳ dị tại góc và đỉnh của bề mặt chảy dẻo Mohr-Coulomb.
  • Tối ưu hóa hiệu năng tính toán thông qua ma trận mô đun tiếp tuyến nhất quán cải tiến theo phương pháp hình học, rút ngắn từ 18% đến 25% thời gian xử lý.
  • Thiết lập hoàn chỉnh chương trình máy tính bằng ngôn ngữ MATLAB code tự động hóa quá trình phân tích phần tử hữu hạn phi tuyến cho bài toán nền móng.
  • Kiểm chứng thành công qua mô phỏng thí nghiệm ba trục và bài toán móng băng với độ chính xác đạt mức sai số chỉ từ 2.1% đến 4.2% so với Plaxis và Terzaghi.
  • Đóng góp nguồn tài liệu học thuật và mã nguồn chuẩn xác, tạo tiền đề làm chủ giải thuật phần mềm địa kỹ thuật tại Việt Nam.

Đóng góp lớn nhất của luận văn là đã chuyển hóa các lý thuyết toán cơ học phức tạp thành công cụ tính toán số cụ thể, trực quan và có độ chính xác cao. Kế hoạch tiếp theo trong giai đoạn 2023 - 2025 là tiếp tục phát triển mã nguồn mở rộng cho các bài toán tương tác đất - cọc không gian 3 chiều và bài toán đào hầm ngầm đô thị. Quý độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể khai thác các công thức và giải thuật trong luận văn để áp dụng vào các dự án phân tích địa kỹ thuật chuyên sâu.