Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình trên nền đất yếu, bài toán dự báo độ lún theo thời gian luôn là thách thức trọng tâm đối với các kỹ sư địa kỹ thuật. Tại các khu vực đồng bằng ven biển và đô thị mới, các lớp sét mềm bão hòa nước thường có chiều dày lớn, độ ẩm tự nhiên cao lên tới 75.7% và hệ số rỗng ban đầu đạt mức 2.104. Dưới tác dụng của tải trọng công trình hoặc tải trọng đắp, quá trình cố kết thoát nước diễn ra rất chậm chạp, có thể kéo dài từ vài chục năm đến hàng trăm năm, gây ra hiện tượng lún kéo dài trong suốt vòng đời khai thác công trình.

Vấn đề cốt lõi đặt ra là lý thuyết cố kết thấm một chiều kinh điển của Terzaghi vốn giả định hệ số cố kết, hệ số rỗng và hệ số thấm là các hằng số không đổi. Tuy nhiên, trong thực tế nén chặt, cấu trúc đất thay đổi liên tục khiến hệ số rỗng và hệ số thấm suy giảm phi tuyến theo sự gia tăng của ứng suất hữu hiệu. Sự khác biệt giữa giả thiết lý thuyết và bản chất cơ học thực tế dẫn đến độ lệch dự báo lún từ 15% đến 25%, gây khó khăn cho công tác kiểm soát tiến độ và an toàn công trình.

Luận văn thạc sĩ chuyên ngành Địa Kỹ thuật Xây dựng do học viên Bùi Anh Kiệt thực hiện dưới sự hướng dẫn của Phó Giáo sư Bùi Trường Sơ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 5 năm 2015) tập trung giải quyết bài toán cố kết thấm một chiều có xét đến sự biến thiên đồng thời của hệ số rỗng và hệ số thấm. Nghiên cứu sử dụng tập dữ liệu gồm 12 mẫu đất sét mềm nguyên dạng tại dự án Khu đô thị mới Nhà Bè, xây dựng mô hình toán học phi tuyến và lập trình mô phỏng trên phần mềm Mathcad. Kết quả nghiên cứu giúp giảm thiểu sai số tính toán xuống dưới 5%, đóng góp giải pháp nâng cao độ tin cậy trong thiết kế xử lý nền đất yếu.

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 cố kết thấm một chiều cổ điển của Karl Terzaghi (1925), kết hợp với các lý thuyết mở rộng về môi trường hai pha gồm cốt đất đàn dẻo và nước lỗ rỗng chứa bọt khí chịu nén tuyến tính theo quan điểm của Bishop (1950), Skempton (1954) và Schuurman. Trong đất loại sét bão hòa, bọt khí hòa tan hoặc bọc kín tạo nên tính nén ép của chất lỏng lỗ rỗng, đồng thời độ bền kiến trúc ban đầu của đất dao động trong khoảng 0.01 đến 0.02 kG/cm2 khiến ứng suất hữu hiệu hình thành ngay từ thời điểm đầu gia tải, sinh ra biến dạng lún tức thời.

Để phản ánh chính xác quy luật nén lún phi tuyến, nghiên cứu áp dụng mô hình hàm số mũ biểu diễn mối quan hệ giữa hệ số rỗng và ứng suất hữu hiệu, đồng thời thiết lập quy luật suy giảm của hệ số thấm theo hàm số mũ tương ứng với hệ số nén thể tích mv đạt 0.103 cm2/kG. Bằng cách sử dụng phép biến đổi biến số trung gian để tuyến tính hóa phương trình vi phân đạo hàm riêng phi tuyến bậc hai, tác giả đã tìm ra nghiệm giải tích dạng chuỗi Fourier cho bài toán phân bố áp lực nước lỗ rỗng thặng dư và độ lún theo thời gian, tích hợp hệ số cố kết thực nghiệm Cv đạt 0.276 x 10^-3 cm2/s.

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

Nguồn dữ liệu phục vụ nghiên cứu được trích xuất từ 12 mẫu đất sét mềm bão hòa nước nguyên dạng thu thập tại công trình thực tế thuộc Khu đô thị mới Nhà Bè, Thành phố Hồ Chí Minh. Phương pháp chọn mẫu có chủ đích được áp dụng nghiêm ngặt nhằm lấy các mẫu đất đại diện tại độ sâu tầng sét yếu chịu nén chính bằng ống mẫu thành mỏng nòng đôi, đảm bảo giữ nguyên kết cấu tự nhiên và độ ẩm ban đầu.

Toàn bộ thí nghiệm nén cố kết không nở hông (Oedometer) được tiến hành theo tiêu chuẩn quốc gia TCVN 4200-2012 trong dao vòng tiêu chuẩn có đường kính từ 50 đến 70 mm và chiều cao 20 mm. Các cấp áp lực nén được gia tải tuần tự từ 25, 50, 100, 200 đến 400 kPa (tương ứng từ 0.25 đến 4.0 kG/cm2). Biến dạng nén được ghi nhận tự động tại các mốc thời gian chuẩn: 15 giây, 30 giây, 1 phút, 2 phút, 4 phút, 8 phút, 15 phút, 30 phút, 1 giờ, 2 giờ, 3 giờ, 6 giờ, 12 giờ và 24 giờ cho mỗi cấp tải. Phương pháp Taylor (căn bậc hai thời gian) và Casagrande (logarit thời gian) được sử dụng đồng thời để xác định chính xác các đại lượng t90, t50, hệ số cố kết Cv và chỉ số nén Cc.

Lý do lựa chọn phương pháp phân tích giải tích kết hợp lập trình số trên phần mềm Mathcad là nhằm mô phỏng liên tục sự suy giảm của hệ số rỗng và hệ số thấm trong từng vi phân thời gian, khắc phục nhược điểm rời rạc của các phép tính thủ công truyền thống. Quá trình nghiên cứu được triển khai thực hiện liên tục từ tháng 8 năm 2014 đến tháng 5 năm 2015.

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

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

Thứ nhất, phân tích thống kê từ 12 mẫu thí nghiệm đất sét mềm Nhà Bè cho thấy hàm lượng hạt mịn (bột và sét) chiếm tỷ lệ áp đảo lên tới 97.5%, trong khi hạt cát chỉ chiếm 2.5%. Đất có độ ẩm tự nhiên 75.7%, giới hạn chảy LL đạt 83.0%, giới hạn dẻo PL là 30.0% và chỉ số dẻo PI lên đến 53.0%. Chỉ số nén lún Cc đạt giá trị cao 1.020 và chỉ số nở Cs là 0.167, phản ánh nền đất có tính nén lún rất mạnh và độ nhạy cảm cao với tải trọng ngoài.

Thứ hai, nghiên cứu đã thiết lập thành công phương trình hồi quy tương quan phi tuyến giữa hệ số rỗng e và ứng suất nén hữu hiệu theo dạng hàm số mũ: e = 2.0876 x exp(-0.192 x ứng suất hữu hiệu), với hệ số tương quan rất cao R^2 = 0.9975. Kết quả này chứng minh rằng hệ số rỗng giảm nhanh trong giai đoạn đầu gia tải từ 0.5 đến 1.0 kG/cm2 và giảm dần tốc độ khi ứng suất tăng lên mức 1.0 đến 2.0 kG/cm2.

Thứ ba, khi so sánh biến dạng lún mẫu đất 471-1 ở cấp áp lực 0.5 đến 1.0 kG/cm2 và 1.0 đến 2.0 kG/cm2, mô hình cố kết mở rộng có xét biến thiên hệ số rỗng và hệ số thấm cho đường cong lún bám sát hoàn toàn dữ liệu đo thực nghiệm tại thời điểm t100. Trong khi lý thuyết Terzaghi cổ điển tạo ra độ lệch lún vượt quá 20% do giả định hệ số thấm không đổi, mô hình cải tiến đã kéo giảm sai số dự báo xuống dưới mức 5%.

Thứ tư, nghiên cứu ghi nhận hệ số cố kết nén thứ cấp C-alpha đạt 0.0034 ở cấp áp lực 0.5 đến 1.0 kG/cm2. Điều này lý giải hiện tượng biến dạng từ biến vẫn tiếp diễn sau khi áp lực nước lỗ rỗng thặng dư đã tiêu tán hoàn toàn ở mốc 24 giờ.

Thảo luận kết quả

Cơ chế vật lý giải thích cho sự vượt trội của mô hình là trong giai đoạn đầu nén lún, sự thoát nước nhanh làm giảm nhanh thể tích lỗ rỗng, dẫn đến sự suy giảm đột ngột của hệ số thấm k. Khi hệ số thấm giảm, tốc độ thoát nước ở các giai đoạn sau bị hãm lại, làm quá trình tiêu tán áp lực nước lỗ rỗng diễn ra chậm hơn so với giả định tuyến tính của Terzaghi. Đồng thời, độ cứng khung hạt tăng dần theo hàm mũ giúp cốt đất chịu lực hữu hiệu sớm hơn.

Kết quả này hoàn toàn phù hợp với các nghiên cứu của Tsytovich và Abelev về ảnh hưởng của độ bền kiến trúc và tính nén ép của chất lỏng lỗ rỗng. Dữ liệu thực nghiệm được trình bày trực quan qua biểu đồ quan hệ e - log p, biểu đồ biến dạng theo căn bậc hai thời gian của Taylor và đường cong phân tán áp lực nước lỗ rỗng theo chiều sâu lớp đất. Việc thể hiện các đường cong tính toán lý thuyết song song với dữ liệu thực đo 24 giờ trên cùng một đồ thị Mathcad đã minh chứng rõ nét tính chính xác của mô hình đề xuất, cung cấp cơ sở tin cậy để tối ưu hóa thời gian gia tải trước cho các dự án san lấp mặt bằng.

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

Thứ nhất, chuẩn hóa quy trình tính toán phi tuyến: Các viện nghiên cứu, trường đại học và tổng công ty tư vấn thiết kế địa kỹ thuật cần đưa mô hình cố kết có xét đến sự thay đổi hệ số rỗng và hệ số thấm vào quy trình tính toán lún tiêu chuẩn trong vòng 6 đến 12 tháng tới. Giải pháp này nhằm nâng cao độ chính xác dự báo độ lún công trình thêm 15% đến 20%, hạn chế rủi ro nứt gãy công trình do lún lệch kéo dài.

Thứ hai, nâng cấp quy chuẩn thí nghiệm nén cố kết: Các phòng thí nghiệm chuyên ngành xây dựng (LAS-XD) cần thực hiện nghiêm túc tiêu chuẩn TCVN 4200-2012, áp dụng đầy đủ các cấp tải trọng từ 25 kPa đến 400 kPa và kéo dài thời gian theo dõi đến 24 giờ cho mỗi cấp nén. Cần bổ sung quy trình xác định hệ số nén thứ cấp C-alpha và chỉ số nén Cc trong vòng 3 tháng tới để cung cấp đầy đủ thông số đầu vào cho bài toán phi tuyến.

Thứ ba, phát triển công cụ phần mềm tính toán tự động: Đội ngũ kỹ sư địa kỹ thuật cần xây dựng các module tính toán tự động trên nền tảng Mathcad, Python hoặc tích hợp vào phần mềm phần tử hữu hạn thương mại trong thời hạn 6 tháng, rút ngắn thời gian xử lý dữ liệu xuống dưới 48 giờ cho mỗi hồ sơ thiết kế nền móng.

Thứ tư, thiết lập hệ thống quan trắc hiện trường đối chứng: Chủ đầu tư và các ban quản lý dự án tại khu vực Nhà Bè và các vùng đất yếu lân cận cần lắp đặt thiết bị đo áp lực nước lỗ rỗng (Piezometer) và mốc đo lún sâu liên tục trong 12 đến 24 tháng thi công. Dữ liệu quan trắc thực tế sẽ là căn cứ hiệu chỉnh các thông số hệ số nén nước lỗ rỗng mw và hệ số thấm ban đầu k0, đảm bảo an toàn tuyệt đối 100% cho công trình.

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

Thứ nhất, kỹ sư thiết kế địa kỹ thuật và công trình ngầm: Luận văn cung cấp công thức giải tích chuẩn xác để tính toán độ lún theo thời gian cho các dự án gia cố nền đất yếu bằng bấc thấm (PVD), cọc xi măng đất hoặc giếng cát, giúp tiết kiệm từ 10% đến 15% chi phí xử lý nền móng nhờ xác định chính xác thời gian dỡ tải gia tải trước.

Thứ hai, học viên cao học, nghiên cứu sinh và giảng viên ngành Xây dựng: Đây là tài liệu tham khảo học thuật giá trị cao về phương pháp giải phương trình vi phân phi tuyến trong cơ học đất, cung cấp phương pháp luận chặt chẽ kết hợp giữa lý thuyết môi trường nhiều pha và mô phỏng số trên Mathcad.

Thứ ba, chủ đầu tư và ban quản lý dự án hạ tầng đô thị: Tài liệu giúp các nhà quản lý đánh giá đúng mức độ rủi ro lún dài hạn tại các khu đô thị ven sông, khu vực san lấp lấn biển có hệ số rỗng ban đầu lớn hơn 2.0, từ đó hoạch định chính xác tiến độ thi công và kế hoạch bảo trì kết cấu hạ tầng.

Thứ tư, chuyên viên thí nghiệm và kiểm định chất lượng địa kỹ thuật: Cán bộ thí nghiệm có thể nắm vững quy trình xử lý số liệu nén cố kết Oedometer, phương pháp ngoại suy Casagrande - Taylor và cách xác định các chỉ tiêu cơ lý phức tạp phục vụ lập báo cáo khảo sát địa chất đạt độ tin cậy trên 95%.

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

Tại sao lý thuyết cố kết cổ điển của Terzaghi lại có sai số đáng kể khi áp dụng cho đất sét mềm Nhà Bè?

Lý thuyết Terzaghi giả thiết hệ số cố kết Cv, hệ số thấm k và hệ số rỗng e là hằng số. Trong thực tế, đất sét mềm Nhà Bè có hệ số rỗng ban đầu lên tới 2.104 và độ ẩm 75.7%. Khi chịu tải từ 0.5 đến 2.0 kG/cm2, đất bị nén chặt làm hệ số thấm suy giảm mạnh theo hàm mũ, khiến tốc độ thoát nước thực tế chậm hơn tính toán cổ điển từ 15% đến 25%.

Mối quan hệ giữa hệ số rỗng và ứng suất nén trong nghiên cứu được biểu diễn bằng công thức nào?

Dựa trên kết quả thí nghiệm nén không nở hông của 12 mẫu đất tại Nhà Bè, mối quan hệ này được biểu diễn bằng phương trình hàm mũ: e = 2.0876 x exp(-0.192 x ứng suất hữu hiệu), đạt hệ số tương quan R^2 = 0.9975. Phương trình này mô tả hoàn hảo quá trình giảm thể tích lỗ rỗng dưới tác dụng của tải trọng nén tăng dần.

Tính nén ép của nước lỗ rỗng chứa bọt khí ảnh hưởng như thế nào đến độ lún ban đầu?

Khi đất chưa bão hòa 100% hoặc chứa khí hòa tan, nước lỗ rỗng có tính nén ép với hệ số nén mw. Do đó, ngay thời điểm đặt tải t = 0, một phần tải trọng ngoài truyền trực tiếp lên khung hạt đất thay vì truyền toàn bộ lên nước lỗ rỗng, làm xuất hiện ứng suất hữu hiệu ban đầu và tạo ra độ lún tức thời chiếm khoảng 5% đến 10% tổng độ lún.

Thí nghiệm nén cố kết trong luận văn được thực hiện theo tiêu chuẩn kỹ thuật nào?

Thí nghiệm được thực hiện chuẩn xác theo tiêu chuẩn quốc gia TCVN 4200-2012 trong phòng thí nghiệm hợp chuẩn. Mẫu đất hình trụ kích thước chiều cao 20 mm được nén qua các cấp áp lực từ 25 đến 400 kPa, theo dõi độ dịch chuyển biến dạng liên tục từ 15 giây đến 24 giờ để xác định các chỉ tiêu Cc = 1.020 và Cv = 0.276 x 10^-3 cm2/s.

Mô hình đề xuất trong luận văn khác biệt như thế nào so với phương pháp Terzaghi tính từng bước?

Phương pháp Terzaghi tính từng bước chỉ cập nhật gián đoạn các thông số trung bình sau mỗi khoảng thời gian, gây ra tích lũy sai số. Ngược lại, mô hình giải tích của luận văn đưa trực tiếp quy luật biến thiên liên tục của hệ số rỗng và hệ số thấm vào phương trình vi phân nén chặt, giúp đường cong lún tiệm cận chính xác dữ liệu thực đo tại mốc t100 với sai số dưới 5%.

Kết luận

  • Luận văn đã phân tích toàn diện các hạn chế của lý thuyết cố kết một chiều Terzaghi khi áp dụng cho đất sét mềm có hệ số rỗng ban đầu cao e0 = 2.104 tại khu vực Nhà Bè.
  • Xây dựng thành công mô hình tương quan thực nghiệm hệ số rỗng theo hàm số mũ e = 2.0876 x exp(-0.192 x ứng suất hữu hiệu) với hệ số tương quan đạt R^2 = 0.9975 trên 12 mẫu đất nguyên dạng.
  • Thiết lập hệ phương trình vi phân cố kết mở rộng có xét đồng thời biến thiên hệ số rỗng, hệ số thấm và tính nén ép của nước lỗ rỗng chứa khí hòa tan.
  • Lập trình tính toán trên Mathcad chứng minh mô hình mới bám sát đường cong nén lún thực nghiệm 24 giờ, giảm thiểu sai số dự báo lún tại mốc t100 xuống dưới 5%.
  • Cung cấp cơ sở khoa học vững chắc phục vụ tính toán, thiết kế xử lý nền đất yếu cho các dự án giao thông, khu đô thị và công trình công nghiệp tại miền Nam.

Trong giai đoạn 12 đến 18 tháng tới, hướng nghiên cứu tiếp theo cần mở rộng bài toán sang không gian hai chiều và ba chiều, kết hợp dữ liệu quan trắc lún hiện trường dài hạn. Các đơn vị tư vấn và kỹ sư công trình hãy áp dụng ngay mô hình tính lún phi tuyến tiên tiến này vào các dự án thực tế để tối ưu hóa thiết kế nền móng và nâng cao chất lượng công trình!