Tổng quan luận án

Công trình "Nghiên cứu đặc tính cơ học của đất yếu được xử lý bằng gia tải trước kết hợp thoát nước thẳng đứng" do nghiên cứu sinh Nguyễn Hồng Trường thực hiện năm 2017 tại Trường Đại học Thủy Lợi, dưới sự hướng dẫn khoa học của TS. Nguyễn Hữu Thái và TS. Nguyễn Tiếp Tân, thuộc chuyên ngành Địa kỹ thuật xây dựng (mã số 62-58-60-01).

Về tính cấp thiết, nền đất yếu phân bố phổ biến tại đồng bằng Bắc Bộ, đồng bằng Nam Bộ cũng như các dải ven sông, ven biển ở Việt Nam. Các lớp trầm tích này mang đặc tính cơ bản là bão hòa nước, hệ số rỗng lớn ($e > 1$), độ sệt cao ($I_L > 1$), khả năng chịu tải nhỏ, biến dạng lún lớn và diễn biến kéo dài theo thời gian. Giải pháp gia tải trước kết hợp bấc thấm (PVD) thường được áp dụng nhằm đẩy nhanh tốc độ thoát nước lỗ rỗng và tăng cường độ chống cắt. Tuy nhiên, hàng loạt sự cố công trình giao thông và dân dụng vẫn xảy ra trong thực tiễn (như sụt trượt, lún vượt dự báo, nứt gãy nền đường). Tác giả xác định ba nguyên nhân kỹ thuật cốt lõi tạo thành khoảng trống nghiên cứu:

  • Đánh giá độ cố kết (ĐCK) chưa chuẩn xác, dẫn đến đắp gia tải vượt quá tốc độ cố kết thực tế của nền; đồng thời thiếu giải pháp bố trí thiết bị quan trắc hợp lý khi số lượng đầu đo bị giới hạn.
  • Đánh giá sai sự gia tăng sức chống cắt không thoát nước ($S_u$) tại các thời điểm cố kết; việc áp dụng công thức tiêu chuẩn (như công thức V.6 trong 22TCN 262-2000) bộc lộ bất cập khi xét các giai đoạn cố kết trung gian; thiếu các hàm dự báo chỉ số nén ($C_c$) và chỉ số nở ($C_s$) phù hợp cho từng vùng trầm tích khi thiếu điều kiện thí nghiệm nén cố kết trực tiếp.
  • Việc áp dụng phương pháp phần tử hữu hạn (PTHH) còn gặp khó khăn do chưa xác định rõ mô hình đất nền phản ánh đúng ứng xử cơ học của đất sét yếu Việt Nam trong bài toán cố kết thấm.

Mục đích nghiên cứu của luận án được văn bản xác định cụ thể theo 4 nhiệm vụ:

  1. Đưa ra khuyến nghị về việc áp dụng phương pháp xác định ĐCK phù hợp cho các giai đoạn cố kết và đề nghị sơ đồ bố trí quan trắc hợp lý nhằm nâng cao hiệu quả đo đạc hiện trường.
  2. Lựa chọn các hàm dự báo chỉ tiêu cơ học ($C_c, C_s$) phù hợp với điều kiện địa chất các vùng đồng bằng Việt Nam khi không có số liệu thí nghiệm trực tiếp.
  3. Làm rõ phạm vi áp dụng của công thức (V.6) trong 22TCN 262-2000 về dự tính sức chống cắt không thoát nước ($S_u$) của đất được gia tải trước và đề xuất công thức tính $S_u$ áp dụng cho thời điểm cố kết bất kỳ.
  4. Lựa chọn mô hình đất phù hợp cho phân tích bài toán cố kết nền sét yếu trong tính toán bằng phương pháp phần tử hữu hạn (phần mềm Plaxis).

Đối tượng nghiên cứu bao gồm: lún cố kết, độ cố kết, các đặc tính cơ học của đất và các mô hình đất dùng trong tính toán địa kỹ thuật. Phạm vi nghiên cứu tập trung vào đất sét yếu điển hình tại đồng bằng Bắc Bộ và đồng bằng Nam Bộ; đất phục vụ thí nghiệm mô hình vật lý được lấy tại khu vực Yên Nghĩa (Hà Đông, Hà Nội); phạm vi biến dạng giới hạn ở quá trình nén lún cố kết thấm dưới tác dụng của tải trọng gia tải trước kết hợp bấc thấm.

Tổng quan tài liệu và vị trí của luận án

Trong Chương 1, tác giả đã hệ thống hóa các thành tựu lý thuyết và thực tiễn liên quan đến bài toán cố kết thoát nước thẳng đứng:

  • Lý thuyết cố kết kinh điển và phát triển: Khởi nguồn từ lý thuyết cố kết một chiều của Terzaghi và phương pháp xác định áp lực tiền cố kết ($p_c$) trên đường cong $e - \log p$ của Casagrande (1936). Lời giải kết hợp cố kết thấm phương đứng và phương ngang được Carrilo (1942) thiết lập qua phương trình $U_{vh} = 1 - (1 - U_h)(1 - U_v)$. Barron (1948) xây dựng lý thuyết cố kết đối xứng trục với hai giả thiết biến dạng đều và biến dạng tự do. Hansbo (1981, 1997) phát triển phương trình giải tích có xét đến đường kính vùng đất bị xáo trộn ($d_s$), hệ số thấm vùng xáo trộn ($k_s$) và sức cản dòng thấm của giếng thoát nước ($q_w, F_r$). Các công thức xác định vùng xáo trộn được điểm lại từ Jamiolkowski và cộng sự (1991), Holtz và Holm (1973), Akagi (1977).
  • Ứng dụng quốc tế và trong nước: Luận án điểm lại các công trình xử lý bấc thấm tại Singapore (Choa và cộng sự), Malaysia (Volders), Thái Lan (Bergado và cộng sự với đất sét mềm Bangkok), Philippines (Tominaga và cộng sự tại Vịnh Manila), Nhật Bản (Suzuki và Yamada), Hàn Quốc (cảng Busan, sân bay Incheon, khu công nghiệp Yeocheon, Kwangyang). Tại Việt Nam, giải pháp bấc thấm kết hợp gia tải nén trước được ghi nhận tại đường cao tốc Cầu Giẽ - Ninh Bình (đoạn nút giao Liêm Tuyền), cao tốc Hà Nội - Hải Phòng (gói thầu EX-9), Đại lộ Đông Tây TP. Hồ Chí Minh (đoạn Thủ Thiêm) cùng các tuyến Quốc lộ 1, Quốc lộ 5.
  • Nghiên cứu đặc trưng cơ học đất tại Việt Nam: Lê Bá Vinh và cộng sự (2007) chỉ ra sự thay đổi hệ số thấm ($k$) theo hệ số rỗng ($e$) làm kéo dài thời gian lún so với giả thiết $k$ không đổi; Nguyễn Thị Nụ nghiên cứu tỷ số hệ số cố kết ngang và đứng ($C_h/C_v$) của trầm tích $amQ_{IV}^{2-3}$ tại Kiên Giang, Tiền Giang, Sóc Trăng, Bến Tre, xác định $C_h/C_v = 6 - 7$ khi áp lực nén $< 0{,}5\text{ kPa}$ và $C_h/C_v = 3 - 3{,}7$ khi áp lực từ $0{,}5 - 4\text{ kPa}$; Dương Diệp Thúy và Phạm Quang Hưng phân tích tương quan 123 mẫu đất dính Hà Nội, xác lập quan hệ $C_c = 8{,}25 C_s$ với hệ số xác định $R^2 = 68{,}3%$; Nguyễn Viết Tình khảo sát 65 mẫu bùn sét phụ tầng Hải Hưng dưới tại Hà Nội, chỉ ra $C_h/C_v = 1{,}06 - 2{,}17$ và $k_h/k_v = 1{,}10 - 1{,}60$; Vũ Đức Sỹ đề xuất khái niệm chiều sâu vùng hoạt động cố kết theo thời gian ($Z_{at}$).
  • Tổng kết sự cố công trình: Tác giả dẫn chứng hàng loạt công trình gặp sự cố nền móng tại Việt Nam do gia tải đắp quá nhanh, dự báo lún sai lệch hoặc khảo sát địa chất không phát hiện túi bùn xen kẹp:
Tên công trình Mô tả thiết kế và địa chất Hiện tượng, sự cố Phân tích nguyên nhân
Km82+500 - Km83+500 Cao tốc Nội Bài - Lào Cai (2014) Nền đắp cao 7-9 m trên đất yếu; xử lý bằng bấc thấm Vết nứt, lún dạng vòng cung kéo dài gần 1 km Nền có túi bùn xen kẹp không được phát hiện; nền lún nhanh hơn dự kiến
Đường dẫn lên cầu Văn Thánh (2003) Nền đắp cao 3-4 m; xử lý bằng bấc thấm Sau 2 năm khai thác nền tiếp tục lún 2-3 cm/tháng, lún vượt 50 cm Rút ngắn thời gian đắp và thời gian chờ lún
Nam cầu Trìa Km732+100 QL 1A (2001) Nền đắp cao 5,4 m; sét xám đen dẻo chảy dày 9,6 m Lún sụt, trượt trồi về phía trái dài 50 m, vùng trượt rộng 26 m Đắp không theo dõi lún, không tính toán chiều cao đắp giới hạn
Km0+620 phía nam cầu vượt đường sắt Hoàng Long (1999) Nền đắp cao 7,5-9 m; bùn hữu cơ dày 4,2-7,4 m Đắp đến 6,8 m thì lún sụt 1,6-1,8 m dài 140 m, hai bên trồi đẩy Đắp tăng tải nhanh vượt quá tốc độ cố kết cần thiết
Đường dẫn cầu Đồng Niên QL 5 (1995) Nền đắp 9 m; sét pha cát dẻo chảy dày 19 m; bấc thấm sâu 17-19 m Độ lún thực tế 2,64-3,3 m sau 3 năm (dự báo 1,7 m); thời gian lún kéo dài trên 10 năm Tính toán, dự báo lún sai lệch so với thực tế
Cư xá Thanh Đa, Bình Thạnh, TP. HCM (2014) Móng băng gia cố cừ tràm dài 2,7 m, mật độ 16 cây/m² trên bùn yếu Lún không đều, góc nghiêng so với phương đứng xấp xỉ 10-20° Phương án móng không phù hợp điều kiện địa chất nền đất yếu
Kho cảng LPG Thị Vải, Bà Rịa - Vũng Tàu (2002-2005) San lấp dày 3,5 m; bùn sét yếu phân bố từ mặt đến sâu 8-23 m; bấc thấm gia tải Tốc độ lún ban đầu lớn, khoảng 1,1-1,2 m/tháng Khảo sát không phản ánh hết tầng đất yếu; thiết kế bỏ qua lún không đều của lớp san lấp

Vị trí của luận án: Kế thừa các lý thuyết cố kết đối xứng trục và mô hình đàn dẻo, luận án trực tiếp giải quyết sự thiếu đồng bộ giữa các phương pháp đánh giá ĐCK trên hiện trường, lấp đầy khoảng trống tính toán cường độ chống cắt $S_u$ ở trạng thái cố kết dở dang, và đưa ra khuyến nghị thực chứng về lựa chọn mô hình PTHH thích hợp cho đất sét yếu Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên khung lý thuyết cố kết thấm phi tuyến và cơ học đất tới hạn. Đất yếu được định nghĩa theo các tiêu chí: $I_L > 1$, $e > 1$, góc ma sát trong $\varphi < 10^\circ$, lực dính cắt nhanh không thoát nước $c < 15\text{ kPa}$, lực dính cắt cánh hiện trường $c_u < 35\text{ kPa}$, sức kháng mũi xuyên tĩnh $q_c < 0{,}1\text{ MPa}$, chỉ số SPT $N < 5$.

Về phương tiện thoát nước thẳng đứng, tác giả tổng hợp ưu nhược điểm kỹ thuật giữa các giải pháp:

Hạng mục Bấc thấm (PVD) Giếng cát Giếng cát có vỏ bọc
Nguyên lý chung Tạo đường thoát nước đứng, rút ngắn chiều dài đường thấm, tăng tốc độ cố kết Tạo đường thoát nước đứng, rút ngắn chiều dài đường thấm, tăng tốc độ cố kết Tạo đường thoát nước đứng, rút ngắn chiều dài đường thấm, tăng tốc độ cố kết
Ưu điểm Chất lượng đồng đều do chế tạo sẵn; ít xáo động nền khi cắm; thi công nhanh; hiệu quả kinh tế cao Tăng nhanh cường độ nền; ít nguy cơ bị cắt dòng thấm trong tầng đất không đồng nhất Tăng nhanh cường độ nền; ít xáo động hơn giếng cát thường; thi công nhanh; chống trượt taluy tốt
Nhược điểm Dòng chảy suy giảm nếu bấc bị xoắn gập do lún lớn; lưu lượng thoát nước giảm theo thời gian Tác động xáo trộn lớn khi thi công; giếng cát có thể bị đứt đoạn do bùn đất xâm nhập Khó thay đổi chiều sâu lắp đặt; kinh nghiệm thi công ở Việt Nam còn hạn chế

Đặc trưng cơ lý trung bình của đất yếu hai miền đồng bằng được tác giả sử dụng làm dữ liệu tham chiếu:

Khu vực địa tầng Loại đất Độ ẩm $W$ (%) Dung trọng tự nhiên $\gamma_w$ (g/cm³) Hệ số rỗng ban đầu $e_0$ Độ bão hòa $G, S$ (%) Giới hạn chảy $LL$ (%) Giới hạn dẻo $PL$ (%) Hệ số nén lún $a_{1-2}$ (cm²/kG) Mô đun biến dạng $E_0$ (kG/cm²) Góc ma sát $\varphi$ (độ) Lực dính $c$ (kG/cm²)
Đồng bằng Bắc Bộ Bùn sét lẫn than bùn 76 1,49 2,07 95,6 60 23 0,16 14 6 0,14
Đồng bằng Bắc Bộ Bùn sét pha cát 42 1,73 1,20 93,8 37 13 0,07 23 16 0,05
Đồng bằng Bắc Bộ Sét xám 42 1,87 1,00 100,0 46 21 0,05 32 10 0,22
Đồng bằng Bắc Bộ Bùn sét 62 1,58 1,71 94,8 55 24 0,06 29 11 0,16
Đồng bằng Nam Bộ Bùn sét ($ambQ_{IV}$) 65,4 1,59 1,81 88,1 59,5 37,0 0,13 12,5 5,8 0,10
Đồng bằng Nam Bộ Bùn á sét ($ambQ_{IV}$) 62,7 1,66 1,68 99,3 47,4 30,6 0,12 16,7 7,3 0,04

Phương pháp nghiên cứu thực nghiệm và mô phỏng của tác giả:

  • Phương pháp thực nghiệm mô hình vật lý (MHVL): Chế tạo 02 thùng mô hình vật lý chứa đất sét yếu khu vực Yên Nghĩa (Hà Nội) ở trạng thái bão hòa, lắp đặt bấc thấm đứng, vải địa kỹ thuật, đệm cát và gia tải theo 2 cấp tải trọng. Bố trí đầu đo áp lực nước lỗ rỗng dây rung (Vibrating Wire Piezometer) kết nối bộ đọc tự động qua phần mềm Logview tại các độ sâu 25 cm, 50 cm, 75 cm và các cự ly cách bấc thấm 15 cm, 50 cm; bố trí hệ thống mốc đo lún bề mặt và đo lún theo chiều sâu. Thí nghiệm vận hành đến khi nền đạt độ cố kết mục tiêu: MHVL 1 đạt ĐCK 50%, MHVL 2 đạt ĐCK 90%.
  • Phương pháp thí nghiệm trong phòng: Tiến hành thí nghiệm nén cố kết một trục (Oedometer) xác định chỉ số nén $C_c$, chỉ số nở $C_s$, áp lực tiền cố kết $p_c$, hệ số cố kết $C_v$; thí nghiệm thấm cột nước giảm dần xác định hệ số thấm đứng $k_v$ và hệ số thấm ngang $k_h$; thí nghiệm nén ba trục cố kết không thoát nước (CU); thí nghiệm cắt cánh hiện trường và cắt cánh trong phòng bằng thiết bị chuyên dụng (Model T0174/B1).
  • Phương pháp mô hình số: Sử dụng phương pháp phần tử hữu hạn thông qua phần mềm Plaxis để mô phỏng bài toán cố kết thấm phẳng và đối xứng trục, quy đổi hệ số thấm tương đương $k_{ve}$, so sánh các mô hình vật liệu: Đàn hồi tuyến tính (Linear Elastic), Mohr-Coulomb, Đàn dẻo tăng bền (Hardening Soil), Cam-Clay cải tiến (Modified Cam-Clay) và Soft Soil.
  • Nguồn dữ liệu: Mẫu đất Yên Nghĩa; dữ liệu khảo sát địa chất và số liệu quan trắc thực tế từ các dự án lớn: Nhà máy nhiệt điện Long Phú 1 (Sóc Trăng), Đường ô tô cao tốc Nội Bài - Lào Cai, Đường Tân Vũ - Lạch Huyện (Hải Phòng), Cao tốc Cầu Giẽ - Ninh Bình, Kho cảng LPG Thị Vải, Oriental Landmark, Tam Hòa (Nam Định), Cầu Kim Cương (Sài Gòn), Kè Cần Thơ, Thốt Nốt, Đồng Tháp.

Nội dung chính theo từng chương

Chương 1: Tổng quan về giải pháp thoát nước thẳng đứng xử lý nền đất yếu

Chương 1 trình bày tổng quan lý luận về nền đất dính mềm yếu và công nghệ xử lý bằng đường thoát nước thẳng đứng kết hợp gia tải trước. Tác giả phân tích địa tầng trầm tích Đệ Tứ của đồng bằng Bắc Bộ (hình thành từ phù sa sông Hồng, sông Thái Bình với hai tầng hạt thô bên dưới và hạt mịn, bùn, than bùn bên trên) và đồng bằng Nam Bộ (ba bậc trầm tích Holoxen với các phân khu chiều dày đất yếu từ dưới 15 m đến trên 30 m). Các phương pháp giải bài toán cố kết được phân tích gồm: đồ giải, giải tích và phương pháp số (sai phân hữu hạn, phần tử hữu hạn).

Tác giả đi sâu vào cơ sở lý thuyết bấc thấm với các công thức quy đổi đường kính tương đương $d_w = 2(a + b)/\pi$, đường kính ảnh hưởng $d_e = 1{,}05S$ (lưới tam giác) và $d_e = 1{,}13S$ (lưới hình vuông). Phương trình Hansbo (1981) được phân tích chi tiết với độ cố kết ngang $U_h = 1 - \exp(-8T_h / F)$, trong đó hệ số tổng hợp $F = F(n) + F_s + F_r$ phản ánh đồng thời khoảng cách bấc $F(n)$, vùng xáo trộn $F_s$ và sức cản giếng $F_r$. Chương này cũng tổng hợp các sự cố sụt trượt thực tế và kết luận rằng việc thiếu các công cụ dự tính chính xác độ cố kết, cường độ chống cắt gia tăng và lựa chọn mô hình tính toán số là nguyên nhân trực tiếp dẫn tới hư hỏng công trình.

Chương 2: Nghiên cứu thí nghiệm mô hình nền sét yếu được gia tải trước kết hợp thoát nước thẳng đứng

Chương 2 mô tả chi tiết công tác thiết kế, chế tạo và vận hành 02 mô hình vật lý (MHVL). Nền đất sét yếu Yên Nghĩa được nhào trộn, đưa vào thùng mô hình với các thông số ban đầu kiểm soát nghiêm ngặt. Hệ thống bấc thấm và thiết bị quan trắc (Vibrating Wire Piezometer, mốc đo lún) được lắp đặt tại các tọa độ không gian xác định. Quá trình gia tải được thực hiện qua hai cấp tải trọng. MHVL 1 được theo dõi đến độ cố kết 50%, MHVL 2 được theo dõi liên tục trong 165 ngày đến khi đạt độ cố kết trên 90%.

Trên cơ sở số liệu thu thập được, tác giả tiến hành tính toán và so sánh ĐCK theo hai hướng tiếp cận:

  • Phương pháp áp lực nước lỗ rỗng (ALNLR): Tính ĐCK cục bộ theo mức độ tiêu tán ALNLR thặng dư $U_u = 1 - u_t / u_{max}$.
  • Phương pháp dự báo lún: Tính ĐCK theo tỷ lệ độ lún tích lũy $U_s = S_t / S_\infty$, trong đó độ lún cuối cùng $S_\infty$ được xác định qua hai phương pháp ngoại suy phổ biến là phương pháp Asaoka và phương pháp Hyperbolic.

Tác giả mở rộng đối chiếu tính toán trên dữ liệu quan trắc thực tế của ba công trình lớn: Dự án Nhà máy nhiệt điện Long Phú 1 (Sóc Trăng), Dự án đường cao tốc Nội Bài - Lào Cai, và Dự án đường ô tô Tân Vũ - Lạch Huyện (Hải Phòng). Kết quả chỉ ra sự sai khác đáng kể giữa hai phương pháp: phương pháp quan trắc lún cho giá trị ĐCK cao hơn ở giai đoạn đầu nhưng phụ thuộc lớn vào thuật toán ngoại suy $S_\infty$; phương pháp ALNLR phản ánh cục bộ, nhạy cảm với bọt khí và quá trình tiêu tán nước không đồng đều quanh bấc thấm. Từ phân tích trường ứng suất và biến dạng, tác giả đề xuất giải pháp kỹ thuật: bố trí điểm quan trắc tại vị trí có ALNLR và độ lún đạt trị số lớn nhất trung bình trong ô lưới bấc thấm để giá trị đo được phản ánh chính xác nhất ĐCK trung bình của toàn nền.

Chương 3: Nghiên cứu các đặc tính cơ học của đất sét yếu

Chương 3 tập trung phân tích sự biến đổi các đặc trưng cơ lý của đất nền trước và sau khi xử lý gia tải từ hai MHVL và hàng loạt mẫu đất hiện trường. Các thí nghiệm ba trục nén cố kết không thoát nước (CU), thí nghiệm thấm ($k_v, k_h$) và nén Oedometer cho thấy sau khi cố kết, hệ số rỗng giảm, dung trọng tăng, lực dính hiệu quả $c'$ và góc ma sát trong hiệu quả $\varphi'$ đều tăng trưởng rõ rệt.

Nội dung trọng tâm thứ nhất của chương là đánh giá độ chính xác của các hàm dự báo chỉ số nén $C_c$ và chỉ số nở $C_s$. Tác giả thu thập và phân tích số liệu địa chất từ 13 dự án trên cả hai miền đất nước (Cầu Giẽ - Ninh Bình, Hà Nội - Hải Phòng, Oriental Landmark, Lạch Huyện, Tam Hòa, Nội Bài - Lào Cai, Cầu Kim Cương, Long Phú 1, LPG Thị Vải, Đồng Tháp, Cần Thơ, Thốt Nốt). Bằng việc so sánh sai số trung bình giữa kết quả tính theo các công thức kinh nghiệm (từ $LL, e_0, W$) với kết quả thí nghiệm nén cố kết trực tiếp, tác giả phân loại và kiến nghị các hàm dự báo có độ tin cậy cao nhất cho từng vùng địa tầng đồng bằng Bắc Bộ và Nam Bộ.

Nội dung trọng tâm thứ hai là nghiên cứu sức chống cắt không thoát nước ($S_u$). Tác giả làm rõ phạm vi áp dụng của công thức (V.6) trong Tiêu chuẩn 22TCN 262-2000: $$\Delta S_u = \Delta\sigma_v \cdot U \cdot \tan\varphi_{\text{cắt nhanh}}$$ Công thức quy phạm này bộc lộ sai số lớn khi dự tính $S_u$ ở các giai đoạn cố kết trung gian ($U < 80%$) do giả định quan hệ tuyến tính đơn giản. Dựa trên chuỗi số liệu thí nghiệm cắt cánh tại MHVL và hiện trường (dự án Tân Vũ, Long Phú 1), tác giả xây dựng công thức cải tiến xác định $S_u$ tại thời điểm cố kết bất kỳ: $$S_{u(t)} = S_{u0} + m \cdot (\Delta\sigma'_v)^n \cdot U^\alpha$$ Trong đó các tham số $m$ và $\alpha$ được chuẩn hóa qua phân tích hồi quy số liệu thực nghiệm cắt cánh. Kết quả kiểm nghiệm tại các mức cố kết $U = 50%$ và $U = 94{,}4%$ tại các dự án Tân Vũ và Long Phú 1 chứng minh công thức đề xuất tiệm cận rất sát với số liệu đo đạc thực tế, khắc phục được độ lệch lớn của công thức quy phạm.

Chương 4: Lựa chọn mô hình đất phù hợp cho phân tích cố kết thoát nước nền sét yếu

Chương 4 khảo sát cơ sở toán học và cơ học của các mô hình cấu tạo đất phổ biến trong phần mềm phần tử hữu hạn Plaxis: Mô hình đàn hồi tuyến tính (Linear Elastic), Mô hình Mohr-Coulomb đàn dẻo lý tưởng, Mô hình đàn dẻo tăng bền Hardening Soil, Mô hình Cam-Clay và Cam-Clay cải tiến (Modified Cam-Clay - MCC), và Mô hình Soft Soil. Tác giả phân tích đường cong quan hệ ứng suất - biến dạng thực nghiệm của các mẫu sét yếu điển hình Bắc Bộ và Nam Bộ, chỉ ra tính chất phi tuyến mạnh, sự suy giảm thể tích lỗ rỗng và hiện tượng tăng bền khi chịu nén cố kết.

Tác giả tiến hành mô phỏng số bài toán cố kết của MHVL trên phần mềm Plaxis bằng cách quy đổi lưới thấm bài toán phẳng qua hệ số thấm tương đương $k_{ve}$. Lần lượt cài đặt các mô hình đất Mohr-Coulomb, Soft Soil và Cam-Clay cải tiến để chạy bài toán cố kết theo thời gian. So sánh đường cong lún và biến thiên tiêu tán ALNLR giữa kết quả mô phỏng và số liệu quan trắc thực tế từ đầu đo cho thấy:

  • Mô hình Mohr-Coulomb không phản ánh đúng quá trình cố kết, dự báo sai lệch lớn về độ lún và tốc độ tiêu tán ALNLR do bỏ qua sự phụ thuộc của mô đun biến dạng vào mức ứng suất và sự thay đổi thể tích dẻo.
  • Mô hình Soft Soil cho kết quả tương đối khả quan về biến dạng nhưng gặp hạn chế trong việc mô tả chính xác áp lực tiền cố kết và đường ứng suất phức tạp.
  • Mô hình Cam-Clay cải tiến (MCC) cho kết quả tương thích nhất với số liệu quan trắc cả về diễn biến lún theo thời gian lẫn quy luật tiêu tan áp lực nước lỗ rỗng.

Để kiểm nghiệm tính đại diện của mô hình Cam-Clay cải tiến, tác giả áp dụng mô phỏng số cho 02 công trình quy mô lớn: Dự án đường ô tô Tân Vũ - Lạch Huyện (Hải Phòng) và Dự án cao tốc Nội Bài - Lào Cai. Kết quả phân tích khẳng định các thông số của mô hình MCC tái hiện sát với số liệu quan trắc lún hiện trường, chứng minh độ tin cậy khi đưa mô hình này vào thiết kế tính toán thực tế.

Kết quả và những đóng góp mới

Luận án đã đạt được các kết quả khoa học và đóng góp mới cụ thể sau:

  • Đóng góp mới về lý luận:
    1. Luận án làm rõ bản chất cơ học dẫn đến sự sai khác giữa hai phương pháp đánh giá độ cố kết (theo áp lực nước lỗ rỗng và theo dự báo lún quan trắc); xác lập cơ sở khoa học để xác định vị trí quan trắc đại diện có khả năng phản ánh chính xác nhất ĐCK trung bình của toàn bộ khối nền xử lý.
    2. Làm rõ phạm vi và giới hạn sai số của công thức quy phạm (V.6) trong 22TCN 262-2000; đề xuất công thức mới dự tính sức chống cắt không thoát nước ($S_u$) của nền đất sét yếu được gia tải trước tại bất kỳ thời điểm cố kết $U$ nào, đồng thời xác lập phương pháp thực nghiệm cắt cánh để xác định các hệ số $m$ và $\alpha$ trong công thức.
    3. Bằng sự kết hợp giữa mô hình vật lý và phương pháp phần tử hữu hạn, luận án chứng minh mô hình Cam-Clay cải tiến (Modified Cam-Clay) là mô hình cấu tạo phù hợp nhất trong số các mô hình đất thông dụng để phân tích bài toán cố kết thấm bằng bấc thấm trong điều kiện trầm tích sét yếu Việt Nam.
  • Đóng góp mới về thực tiễn và giải pháp:
    1. Đề xuất vị trí bố trí thiết bị quan trắc áp lực nước lỗ rỗng và đo lún hợp lý (tại điểm có trị số lớn nhất trung bình), giúp tối ưu hóa số lượng thiết bị đo đạc mà vẫn kiểm soát an toàn tốc độ đắp gia tải trên công trường.
    2. Cung cấp bộ hàm dự báo chỉ số nén $C_c$ và chỉ số nở $C_s$ đã được kiểm chứng sai số thực nghiệm, cho phép các kỹ sư áp dụng trực tiếp trong giai đoạn lập dự án và thiết kế sơ bộ khi chưa có điều kiện thí nghiệm nén Oedometer.
    3. Đưa ra quy trình phân tích số bằng Plaxis sử dụng mô hình Cam-Clay cải tiến với hệ thống thông số đầu vào được kiểm nghiệm qua các công trình thực tế cấp quốc gia (Tân Vũ - Lạch Huyện, Nội Bài - Lào Cai).

Hạn chế và hướng nghiên cứu tiếp

Về mặt hạn chế, văn bản luận án nêu rõ:

  • Kích thước hình học của thùng mô hình vật lý còn mang tính hữu hạn, chịu ảnh hưởng nhất định của hiệu ứng biên (dù tác giả đã thực hiện mô phỏng PTHH để đánh giá và giảm thiểu sai số biên).
  • Nghiên cứu tập trung giải quyết thành phần lún cố kết sơ cấp (lún thấm); chưa đi sâu định lượng biến dạng từ biến (lún cố kết thứ cấp $C_\alpha$) - một thành phần có thể chiếm tỷ trọng đáng kể đối với đất than bùn hoặc bùn sét chứa nhiều hàm lượng hữu cơ.
  • Số lượng công trình kiểm nghiệm thực tế mới tập trung chủ yếu ở một số dự án trọng điểm phía Bắc (Hải Phòng, Lào Cai) và Tây Nam Bộ (Sóc Trăng), chưa bao phủ toàn bộ các phân khu địa chất công trình khác.

Về hướng nghiên cứu tiếp theo:

  • Tiếp tục mở rộng các nghiên cứu thực nghiệm xét đến biến dạng từ biến của khung cốt hạt đất theo thời gian dài hạn sau khi kết thúc quá trình cố kết sơ cấp.
  • Tiếp tục thu thập số liệu cắt cánh tại nhiều vùng địa chất ven biển miền Trung và các tiểu vùng đồng bằng khác nhằm làm phong phú thêm ngân hàng dữ liệu cho các tham số $m, \alpha$ trong công thức tính $S_u$.
  • Nghiên cứu ứng dụng bài toán phần tử hữu hạn không gian 3 chiều (3D) mô phỏng chi tiết tương tác cục bộ giữa bấc thấm, lớp đệm cát và tải trọng đắp phức tạp.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo hữu ích cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học, giảng viên: Phục vụ công tác nghiên cứu và giảng dạy các chuyên đề chuyên sâu về Cơ học đất, Địa kỹ thuật công trình, Xử lý và cải tạo nền đất yếu tại các trường đại học khối kỹ thuật xây dựng, thủy lợi, giao thông.
  • Kỹ sư tư vấn thiết kế và giám sát công trình: Cung cấp phương pháp lựa chọn nhanh chỉ số nén lún $C_c, C_s$, công thức dự báo sự gia tăng sức chống cắt $S_u$ để lên biểu đồ tiến độ đắp gia tải nhiều giai đoạn; hướng dẫn thiết lập thông số mô hình Cam-Clay cải tiến trong phần mềm Plaxis.
  • Nhà thầu thi công và đơn vị quan trắc hiện trường: Cung cấp giải pháp bố trí đầu đo piezometer và mốc đo lún hợp lý nhằm xác định đúng thời điểm dỡ tải an toàn, tiết kiệm chi phí và phòng ngừa nguy cơ sụt trượt.

Các phần có giá trị tham khảo nổi bật bao gồm: Chương 2 với phương pháp luận phân tích ĐCK kết hợp giữa ALNLR và lún ngoại suy; Chương 3 với bộ dữ liệu so sánh sai số các hàm dự báo $C_c, C_s$ và công thức cải tiến tính $S_u$; Chương 4 với quy trình thiết lập bài toán phẳng tương đương mô phỏng bấc thấm trong Plaxis.

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

1. Theo số liệu trong luận án, đất sét yếu tại Việt Nam được định lượng qua những chỉ tiêu cơ lý đặc trưng nào?
Theo tổng kết của tác giả ở Chương 1, đất sét yếu tại Việt Nam có hệ số rỗng $e > 1$, độ sệt $I_L > 1$, góc ma sát trong $\varphi < 10^\circ$, lực dính cắt nhanh không thoát nước $c < 15\text{ kPa}$, lực dính cắt cánh hiện trường $c_u < 35\text{ kPa}$, sức kháng mũi xuyên tĩnh $q_c < 0{,}1\text{ MPa}$, và chỉ số SPT $N < 5$. Chiều dày tầng sét yếu thay đổi từ vài mét đến hơn 100 m tùy khu vực trầm tích.

2. Tại sao lại có sự sai khác giữa độ cố kết tính theo quan trắc lún và độ cố kết tính theo áp lực nước lỗ rỗng?
Luận án chỉ ra rằng phương pháp đo lún phụ thuộc vào thuật toán ngoại suy độ lún cuối cùng $S_\infty$ (qua phương pháp Asaoka hoặc Hyperbolic) và bị ảnh hưởng bởi độ lún tức thời cùng biến dạng từ biến chưa tách lọc. Ngược lại, phương pháp ALNLR phản ánh cục bộ quá trình tiêu tán áp lực nước tại vị trí đặt đầu đo, dễ bị sai số do bọt khí lọt vào thiết bị đo hoặc sự không đồng nhất của dòng thấm quanh bấc. Do đó, hai phương pháp thường không cho cùng một giá trị ở cùng một thời điểm.

3. Công thức (V.6) trong Tiêu chuẩn 22TCN 262-2000 có nhược điểm gì khi tính sức chống cắt không thoát nước gia tăng?
Công thức (V.6) trong tiêu chuẩn giả định sự gia tăng sức chống cắt tỷ lệ thuận tuyến tính đơn giản với độ cố kết $U$ và tải trọng gia tải, vốn chỉ phản ánh tương đối phù hợp khi độ cố kết đã đạt mức rất cao (kết thúc cố kết). Khi nền ở trạng thái dở dang ($U$ trung gian), công thức này cho kết quả chênh lệch lớn so với thí nghiệm cắt cánh thực tế. Luận án đã khắc phục bằng cách thiết lập công thức mới tích hợp số mũ $\alpha$ và hệ số $m$ rút ra từ thực nghiệm.

4. Vì sao mô hình Cam-Clay cải tiến (MCC) lại ưu việt hơn mô hình Mohr-Coulomb khi tính toán cố kết thấm bằng phần mềm Plaxis?
Mô hình Mohr-Coulomb giả thiết vật liệu đàn dẻo lý tưởng với mô đun đàn hồi cố định, không phản ánh được sự giảm thể tích lỗ rỗng, sự biến cứng của đất khi nén chặt và sự thay đổi của áp lực nước lỗ rỗng theo trạng thái ứng suất hiệu quả. Ngược lại, mô hình Cam-Clay cải tiến mô tả chính xác mối quan hệ phi tuyến thể tích - ứng suất, sự phụ thuộc của mặt dẻo vào áp lực tiền cố kết, từ đó cho kết quả lún và tiêu tán ALNLR trùng khớp với quan trắc thực tế.

5. Luận án ghi nhận quan hệ tương quan thực nghiệm nào giữa chỉ số nén $C_c$ và chỉ số nở $C_s$ đối với đất dính Hà Nội?
Dẫn chứng nghiên cứu của Dương Diệp Thúy và Phạm Quang Hưng trên 123 mẫu đất thí nghiệm tại khu vực Hà Nội, luận án ghi nhận tương quan tuyến tính giữa chỉ số nén đi và nén lại là $C_c = 8{,}25 C_s$ với hệ số tương quan xác định $R^2 = 68{,}3%$.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Hồng Trường đã hoàn thành toàn diện các mục tiêu nghiên cứu đặt ra, giải quyết các bất cập cốt lõi trong phân tích cơ học đất và công nghệ xử lý nền sét yếu bằng gia tải trước kết hợp thoát nước thẳng đứng. Bằng việc kết hợp thực nghiệm mô hình vật lý, thí nghiệm trong phòng, dữ liệu quan trắc hiện trường và mô phỏng số phần tử hữu hạn, tác giả đã làm sáng tỏ quy luật sai khác của các phương pháp đánh giá độ cố kết, đề xuất thành công công thức xác định sức chống cắt không thoát nước $S_u$ cho mọi giai đoạn cố kết, và chứng minh mô hình Cam-Clay cải tiến là công cụ số thích hợp nhất cho đất sét yếu Việt Nam. Các kết quả và đề xuất kỹ thuật trong công trình là đóng góp có ý nghĩa khoa học và thực tiễn cao đối với ngành Địa kỹ thuật xây dựng công trình tại Việt Nam.