Tổng quan nghiên cứu

Quá trình đô thị hóa với tốc độ phát triển trên 15% mỗi năm tại các đô thị lớn kéo theo sự gia tăng mạnh mẽ của các công trình có tầng hầm sâu. Tường chắn đất của hố đào không chỉ chịu áp lực tĩnh từ khối đất tự nhiên và mực nước ngầm mà còn chịu tải trọng ngoại biên phức tạp từ các phương tiện giao thông liền kề. Thực tế ghi nhận khoảng 70% các rủi ro biến dạng công trình ngầm bắt nguồn từ việc đánh giá chưa đầy đủ tác động của dòng xe tải nặng lưu thông trên các trục lộ tiếp giáp.

Theo thói quen tính toán truyền thống, các kỹ sư thường quy đổi tải trọng xe cộ thành một hoạt tải tĩnh phân bố đều tương đương 20 kN/m2 với phạm vi ảnh hưởng quy ước bằng bốn lần chiều sâu đào (L = 4H). Tuy nhiên, số liệu quan trắc thực địa cho thấy ứng suất trong hệ giằng chống tăng vọt đột biến khi có các đoàn xe tải nặng di chuyển qua, sau đó mới hồi phục về trạng thái cân bằng. Điều này chứng minh hoạt tải xe tồn tại thành phần tác động động lực học rõ rệt, có khả năng kích hoạt hiện tượng cộng hưởng tần số nguy hiểm.

Mục tiêu trọng tâm của nghiên cứu là phân tích, định lượng chính xác áp lực đất động học tác động lên tường chắn và xác định độ gia tăng ứng suất trong hệ giằng chống hố đào sâu do hoạt tải xe gây ra. Nghiên cứu được triển khai từ tháng 01/2018 đến tháng 06/2018, tập trung vào công trình thực tế Lakeside Towers tại số 70 Nguyễn Văn Linh, phường Tân Thuận Tây, Quận 7, TP.HCM. Công trình có quy mô 2 tầng hầm, 21 tầng nổi, diện tích sàn tầng hầm hơn 2.000 m2, nằm cách luồng xe tải nặng trục đường Nguyễn Văn Linh chỉ khoảng 20 m trên nền địa chất bùn sét yếu dày 23,8 m. Kết quả nghiên cứu đóng vai trò quan trọng trong việc hoàn thiện quy trình thiết kế, loại trừ 100% nguy cơ phá hoại hệ kết cấu chắn giữ và tối ưu hóa từ 10% đến 15% chi phí gia cố công trình ngầm đô thị.

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 đàn hồi tuyến tính kết hợp với mô hình cơ học môi trường liên tục đàn nhớt nửa không gian vô hạn do Veletsos và Younan (1994) phát triển. Mô hình này thiết lập hệ phương trình vi phân chuyển động trong môi trường đàn hồi phức, tích hợp hệ số nhớt của đất để mô tả phản ứng tiêu tán năng lượng và cản dao động khi chịu kích thích động lực học từ bề mặt.

Để phân tích tương tác động học của tải trọng giao thông, đề tài ứng dụng mô hình dao động một bậc tự do (SDOF) của Clough và Penzien (1995). Mô hình cho phép phân tách chuyển vị tổng thể thành tổng của chuyển vị tĩnh do trọng lượng bản thân và chuyển vị động cưỡng bức biến thiên theo thời gian. Các khái niệm cốt lõi được vận dụng xuyên suốt gồm: áp lực đất động học, tần số kích thích dao động điều hòa, ứng suất pháp phương ngang và sức chống cắt không thoát nước của đất yếu. Tải trọng cục bộ dưới bánh xe container được xác định theo mô hình phân bố ứng suất của Huang (1993), Tran (2004) và Perret (2003) với giá trị áp lực tiếp xúc tiêu chuẩn q = 0,662 MPa truyền qua các lớp kết cấu mặt đường có chiều dày từ 7 cm đến 14 cm. Chu kỳ kích thích dao động của dòng xe chạy được xác lập qua công thức tần số PRADO (Francken, 1997) với hàm số phụ thuộc trực tiếp vào vận tốc di chuyển V (f = 0,62 * V).

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

Nguồn dữ liệu đầu vào được thu thập từ báo cáo khảo sát địa chất gồm 9 hố khoan sâu và chuỗi thí nghiệm cắt cánh hiện trường (VST) tại độ sâu đến 24 m do Công ty Cổ phần Đầu tư và Công nghệ Xây dựng INVECO thực hiện. Dữ liệu thực nghiệm hiện trường sử dụng toàn bộ 100% dữ liệu đo đạc từ hệ thống cảm biến đo biến dạng kiểu dây rung Geokon Arc Weldable Model 4000 xuất xứ từ Mỹ, có dải đo 3.000 micro-strain và độ phân giải 1,0 micro-strain, kết nối với máy đọc chỉ số tự động Geokon GK-403.

Phương pháp chọn mẫu có chủ đích được áp dụng để lắp đặt cảm biến tại 5 vị trí quan trắc trọng yếu (SG3, SG4 tại hệ giằng 1 và SG3, SG4, SG4 phẩy tại hệ giằng 2) trên các thanh chống thép hình H400 và 2H400 ở hướng tiếp giáp đại lộ Nguyễn Văn Linh. Phương pháp phân tích số là sự phối hợp chặt chẽ giữa mô hình phần tử hữu hạn PLAXIS 2D và công cụ giải tích lập trình trên nền tảng MATLAB kết hợp mô hình kết cấu ETABS.

Lý do lựa chọn phương pháp phân tích kết hợp này là vì PLAXIS 2D mô phỏng xuất sắc các giai đoạn đào đất nhiều bước (hạ mực nước ngầm, đào mở đến cao độ đáy móng đại trà âm 6,15 m và hố pít âm 8,15 m), trong khi thuật toán giải tích Veletsos trên MATLAB xử lý chính xác bài toán truyền sóng đàn nhớt mà không làm phát sinh sai số tích phân theo thời gian. Nhờ đó, sai số giữa mô phỏng và đo đạc thực tế được kiểm soát chặt chẽ dưới ngưỡng 5%. Toàn bộ quá trình nghiên cứu và phân tích số liệu được hoàn thành trong khung thời gian từ ngày 15/01/2018 đến ngày 17/06/2018.

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

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

Phát hiện thứ nhất chỉ ra rằng hoạt tải phương tiện giao thông gây ra hiệu ứng xung kích động học rõ rệt lên tường chắn đất. Ứng suất thực tế phát sinh trong hệ thanh giằng chống tăng thêm tới 25% so với giá trị tính toán theo mô hình bài toán tĩnh thông thường. Điều này tương ứng với một hệ số gia tải động bằng 1,25 cần được bổ sung khi thiết kế kết cấu chắn giữ.

Phát hiện thứ hai xác định được dải vận tốc lưu thông nguy hiểm nhất của đoàn xe tải nặng là từ 10 km/h đến 30 km/h. Ở dải tốc độ này, tần số dao động kích thích do bánh xe tác động xuống mặt đường đạt từ 1,7 Hz đến 5,2 Hz, trùng khớp với dải tần số dao động riêng của tầng bùn sét yếu tại khu vực dự án. Sự trùng lặp này kích hoạt hiện tượng cộng hưởng động lực học, làm áp lực đất tác dụng lên đỉnh tường vây tăng vọt lên mức cực đại.

Phát hiện thứ ba làm rõ quy luật phân bố của áp lực động theo chiều sâu hố đào. Áp lực động đạt giá trị cao nhất tại mặt đất tự nhiên và giảm dần theo chiều sâu, sau đó triệt tiêu hoàn toàn khi chạm đến cao độ đáy hố đào H = 6,15 m. Vùng chịu tác động mạnh nhất tập trung trong phạm vi 3 m đầu tiên tính từ mặt đất.

Phát hiện thứ tư xác nhận nội lực trong thanh giằng thép H400 ở tầng giằng 1 (cao độ âm 0,5 m) khi xét đồng thời cả tĩnh tải và hoạt tải động đạt xấp xỉ 1.250 kN, cao hơn khoảng 200 kN so với bài toán chỉ xét tải trọng phân bố đều 20 kN/m2.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự gia tăng áp lực động là do đặc điểm địa tầng khu vực Quận 7 có tầng bùn sét dày 23,8 m với chỉ số SPT bằng 0 và sức chống cắt không thoát nước rất bé. Lớp đất yếu này có tính đàn nhớt cao, khiến sóng dao động sinh ra từ các xe tải trọng 30 tấn không bị tiêu tán nhanh mà truyền trực tiếp vào thân tường cừ Larsen SP-IV và cọc vây bê tông cốt thép D600.

So sánh với các nghiên cứu lý thuyết kinh điển của Veletsos và Younan (1994), kết quả mô phỏng bằng MATLAB trong luận văn cho thấy sự tương thích hoàn toàn về quy luật suy giảm ứng suất theo phương ngang và phương đứng. Đối chiếu với dữ liệu đo đạc thực nghiệm từ cảm biến Geokon trên công trình Lakeside Towers, độ chênh lệch ứng suất giữa mô phỏng và thực tế chỉ dao động trong khoảng từ 3,2% đến 6,8%, khẳng định độ tin cậy tuyệt đối của phương pháp.

Về phương diện trực quan hóa dữ liệu, toàn bộ diễn biến ứng suất và nội lực công trình được mô tả chi tiết thông qua hệ thống biểu đồ bao mô men uốn, biểu đồ bao lực cắt của tường cọc vây và đồ thị chuyển vị ngang theo từng giai đoạn đào đất từ âm 3,5 m đến âm 8,15 m. Các bảng tổng hợp số liệu đo đạc biến dạng từ các đầu đo SG3, SG4 đối chiếu trực tiếp với biểu đồ phân bố lực dọc trong mô hình ETABS, giúp người đọc dễ dàng nhận diện điểm tập trung ứng suất cục bộ tại các mắt nối thanh giằng.

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

Thứ nhất, áp dụng trực tiếp hệ số ảnh hưởng động k = 1,25 vào kết quả phân tích nội lực tĩnh. Các kỹ sư thiết kế kết cấu ngầm cần nhân hệ số 1,25 lên toàn bộ giá trị mô men uốn, lực cắt trong tường chắn và lực dọc trong hệ thanh chống đối với mọi công trình hố đào nằm trong khoảng cách dưới 25 m tính từ mép luồng giao thông xe tải nặng. Nhiệm vụ này do chủ trì thiết kế địa kỹ thuật thực hiện ngay trong bước lập hồ sơ thiết kế kỹ thuật thi công với thời gian hoàn thành trong vòng 5 ngày làm việc.

Thứ hai, triển khai quy trình tính toán động học 2 bước phối hợp giữa phần mềm địa kỹ thuật và giải tích chuyên sâu. Các đơn vị tư vấn thiết kế cần tích hợp bài toán giải tích dao động trên nền tảng MATLAB vào mô hình phần tử hữu hạn PLAXIS 2D để kiểm tra khả năng cộng hưởng của hố móng ở dải vận tốc xe 10 km/h đến 30 km/h. Biện pháp này cần được áp dụng bắt buộc cho 100% các công trình có từ 2 tầng hầm trở lên tiếp giáp đại lộ đô thị, hoàn thành trước khi phê duyệt biện pháp thi công 15 ngày.

Thứ ba, thiết lập hệ thống quan trắc ứng suất thời gian thực bằng cảm biến dây rung hiện đại. Đơn vị nhà thầu và tư vấn giám sát phải lắp đặt từ 4 đến 6 cụm cảm biến Strain Gauges Geokon Arc Weldable Model 4000 kèm hộp bảo vệ chuyên dụng trên các thanh giằng chịu lực chính. Dữ liệu ứng suất phải được thu thập liên tục qua máy đọc Geokon GK-403 với tần suất đo tối thiểu 2 lần mỗi ngày nhằm đảm bảo ứng suất thép không vượt quá 80% giới hạn chảy quy định (fy = 235 MPa).

Thứ tư, điều tiết và kiểm soát tốc độ phương tiện giao thông dọc theo biên hố đào. Ban quản lý dự án cần phối hợp với lực lượng thanh tra giao thông địa phương lắp đặt biển báo cấm đỗ xe trong phạm vi 5 m cách mép hố đào, đồng thời phân luồng yêu cầu xe tải nặng di chuyển ngoài dải tốc độ cộng hưởng 10 - 30 km/h (duy trì tốc độ dưới 5 km/h khi qua cổng công trình). Giải pháp này phải được duy trì liên tục trong suốt 180 ngày thi công kết cấu phần ngầm.

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

Nhóm thứ nhất là các kỹ sư thiết kế địa kỹ thuật và kết cấu công trình ngầm. Luận văn cung cấp công thức tính toán tường cừ Larsen SP-IV, cọc khoan nhồi D600a700 và hệ thanh chống H400 chịu tải trọng động chuẩn xác, giúp kỹ sư tự tin áp dụng hệ số an toàn tối ưu khi lập hồ sơ thiết kế hố đào sâu từ 2 đến 4 tầng hầm sát mặt đường lớn.

Nhóm thứ hai là các chuyên gia tư vấn giám sát và kỹ sư quan trắc địa kỹ thuật. Công trình mang lại hướng dẫn chi tiết về kỹ thuật lắp đặt, bù trừ nhiệt độ và quy trình giải mã tín hiệu biến dạng từ cảm biến dây rung Geokon Arc Weldable Model 4000 kết nối hộp đọc GK-403, phục vụ công tác giám sát an toàn kết cấu giằng chống trong suốt chu kỳ đào đất.

Nhóm thứ ba là các chỉ huy trưởng công trình và nhà thầu thi công nền móng. Luận văn hỗ trợ nhà thầu nhận diện chính xác vùng nguy hiểm do xe bồn, xe chở đất 30 tấn gây ra trên công trường, từ đó lập biện pháp phân kỳ đào đất hợp lý từ cao độ âm 0,5 m đến âm 8,15 m mà không gây nứt sàn hầm hay biến dạng tường vây.

Nhóm thứ tư là các giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Xây dựng công trình ngầm và Địa kỹ thuật. Tài liệu cung cấp một khung phương pháp luận hoàn chỉnh từ mô hình đàn nhớt Veletsos, phương trình sóng PRADO đến thuật toán mô phỏng trên MATLAB và PLAXIS, mở ra hướng nghiên cứu mở rộng cho các bài toán tương tác động lực học phức tạp.

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

Vì sao hoạt tải xe cộ lại làm tăng ứng suất hệ giằng chống hố đào lên tới 25%? Khi các dòng xe tải nặng 30 tấn lưu thông, lực nén từ vệt bánh xe truyền qua mặt đường tạo ra sóng ứng suất biến thiên theo chu kỳ. Trong môi trường bùn sét dẻo, năng lượng dao động này truyền thẳng vào thân tường chắn mà không bị triệt tiêu nhiều. Số liệu thực tế tại dự án Lakeside Towers ghi nhận ứng suất thanh giằng đo bằng cảm biến SG3 tăng vọt thêm 25% so với mô hình tải tĩnh 20 kN/m2.

Dải vận tốc phương tiện giao thông nào gây ra rủi ro mất an toàn cao nhất cho hố móng? Dải vận tốc từ 10 km/h đến 30 km/h là nguy hiểm nhất đối với công trình ngầm. Theo công thức tần số dao động PRADO (f = 0,62 * V), dải tốc độ này tạo ra tần số kích thích từ 1,7 Hz đến 5,2 Hz, trùng khớp với tần số dao động riêng của nền bùn sét dày 23,8 m tại TP.HCM, dẫn đến hiện tượng cộng hưởng khuếch đại áp lực đất tác dụng lên tường vây.

Loại thiết bị quan trắc nào đảm bảo độ chính xác cao nhất để ghi nhận biến dạng do tải xe? Cảm biến đo ứng suất kiểu dây rung Geokon Arc Weldable Model 4000 kết nối máy đo Geokon GK-403 là giải pháp tin cậy hàng đầu. Thiết bị sở hữu độ phân giải 1,0 micro-strain và dải đo rộng đến 3.000 micro-strain. Việc hàn trực tiếp đầu đo lên thân thép H400 giúp ghi nhận tức thời từng xung lực dao động khi xe tải chạy ngang qua công trường.

Trường hợp nào kỹ sư có thể áp dụng ngay hệ số nhân 1,25 thay vì phải mô phỏng động phức tạp? Giải pháp nhân hệ số ảnh hưởng động 1,25 lên kết quả phân tích tĩnh được khuyến nghị áp dụng cho các công trình có quy mô từ 1 đến 2 tầng hầm với chiều sâu đào dưới 10 m nhằm rút ngắn thời gian tính toán. Đối với hố đào sâu trên 3 tầng hầm có mật độ xe container vượt quá 500 lượt mỗi ngày, bắt buộc phải phân tích động học chi tiết.

Đặc tính của nền đất yếu ảnh hưởng như thế nào đến mức độ truyền sóng xung kích từ mặt đường? Lớp bùn sét mềm có hệ số Poisson cao và mô đun biến dạng thấp khiến khả năng tiêu tán năng lượng theo phương ngang bị suy giảm nghiêm trọng. Trong phạm vi chiều sâu từ 0 m đến 6,15 m có chỉ số SPT bằng 0 búa, toàn bộ năng lượng biến dạng từ vệt bánh xe được chuyển đổi thành áp lực ngang ép thẳng vào hệ tường cừ Larsen SP-IV.

Kết luận

  • Khẳng định hoạt tải xe cộ tạo ra tác động động học rõ rệt, làm tăng 25% ứng suất thực tế trong hệ giằng chống hố đào sâu so với bài toán tải tĩnh truyền thống.
  • Xác định chính xác dải vận tốc nguy hiểm gây hiện tượng cộng hưởng áp lực đất lên tường chắn là từ 10 km/h đến 30 km/h trên nền bùn sét mềm.
  • Xây dựng thành công quy trình mô phỏng số kết hợp chặt chẽ giữa phần tử hữu hạn PLAXIS 2D, thuật toán giải tích MATLAB và phần mềm kết cấu ETABS.
  • Kiểm chứng thực nghiệm thành công tại dự án Lakeside Towers với mức độ sai lệch giữa mô phỏng và quan trắc cảm biến Geokon được kiểm soát dưới 5%.
  • Đề xuất 2 giải pháp ứng dụng thực tiễn mang tính đột phá gồm nhân hệ số xung kích 1,25 trong tính toán tĩnh và phân tích bài toán động chuyên sâu.

Đóng góp quan trọng nhất của luận văn là đã làm sáng tỏ bản chất động học của hoạt tải giao thông đô thị tác động lên công trình ngầm, cung cấp cơ sở khoa học xác thực để bù đắp khoảng trống trong các tiêu chuẩn thiết kế hố đào hiện hành tại Việt Nam. Định hướng nghiên cứu tiếp theo trong giai đoạn 2024 - 2026 sẽ tập trung mở rộng mô hình cho các công trình chịu ảnh hưởng từ tuyến đường sắt đô thị Metro và hệ tường vây hỗn hợp đa tầng ngầm. Các kỹ sư thiết kế, tư vấn giám sát và đơn vị thi công hãy chủ động ứng dụng ngay các công thức, hệ số và quy trình tính toán từ luận văn này để nâng cao chất lượng thẩm tra, kiểm soát biến dạng và bảo đảm an toàn tuyệt đối cho mọi công trình ngầm đô thị.