Tổng quan về luận án

Nghiên cứu ứng xử động lực học của hệ thống nền móng công trình là một trong những bài toán phức tạp và then chốt của kỹ thuật địa kỹ thuật hiện đại. Tại khu vực Thành phố Hồ Chí Minh (TP.HCM) và các vùng phụ cận—nơi có tốc độ đô thị hóa và công nghiệp hóa vượt bậc với sự xuất hiện dày đặc của các khu công nghiệp trọng điểm như Cát Lái (Quận 2), Hiệp Phước (Nhà Bè) hay Tân Tạo (Bình Chánh)—các công trình cao tầng và móng máy công suất lớn thường xuyên đối mặt với địa tầng đất yếu trầm tích ven biển, cửa sông bão hòa nước có chiều dày dao động từ 10m đến hơn 30m, với các thông số cơ lý đặc trưng rất bất lợi: dung trọng tự nhiên $\gamma = 14 \div 16\text{ kN/m}^3$, mô đun biến dạng $E < 6000\text{ kPa}$, hệ số rỗng $e > 1.8$, lực dính có hiệu $c' < 15\text{ kPa}$, sức kháng cắt không thoát nước $S_u < 50\text{ kPa}$ và độ ẩm tự nhiên $W \ge 65%$. Khi chịu tác động của tải trọng động tuần hoàn từ móng máy, dòng xe cơ giới hạng nặng hay sóng địa chấn lan truyền, nền đất và cọc xuất hiện hiện tượng tích lũy biến dạng dẻo, gia tăng áp lực nước lỗ rỗng ($u$), dẫn đến sự suy giảm nghiêm trọng ứng suất pháp hữu hiệu $\sigma'_z = \sigma_z - u$ và khả năng chịu tải của kết cấu móng.

                  +---------------------------------------------------+
                  |          TẢI TRỌNG ĐỘNG TUẦN HOÀN (MÓNG MÁY)       |
                  |             Tần số f = 10 - 50 Hz, Biên độ        |
                  +-------------------------+-------------------------+
                                            |
                                            v
                  +---------------------------------------------------+
                  |               HỆ ĐÀI CỌC TUYỆT ĐỐI CỨNG           |
                  |                (Phân phối sóng ứng suất)          |
                  +-------------------------+-------------------------+
                                            |
                         +------------------+------------------+
                         |                                     |
                         v                                     v
+------------------------------------------------+   +------------------------------------------------+
|     MA SÁT THÂN CỌC THEO ĐỘ SÂU (Fs0, Fs1)     |   |          SỨC KHÁNG MŨI CỌC (q_p)               |
|  - Sóng ứng suất truyền dọc trục: vc = sqrt(E/rho) |   |  - Tập trung ứng suất dẻo cực hạn              |
|  - Độ nhám bề mặt cọc: Trơn vs. Nhám           |   |  - Lớp cát chặt vừa: N_SPT = 15 - 20           |
|  - Tương tác tiếp xúc đất - cọc theo tần số    |   |  - Phá hoại trượt cắt dẻo tại vùng mũi         |
+------------------------+-----------------------+   +------------------------+-----------------------+
                         |                                     |
                         +------------------+------------------+
                                            |
                                            v
                  +---------------------------------------------------+
                  |    MÔ HÌNH HÓA SỐ FEM (MODIFIED CAM-CLAY - MCC)   |
                  |    - Mặt ngưỡng chảy dẻo p'-q, Hệ số M, lambda, kappa  |
                  |    - Phân tích ngược Thí nghiệm nén tĩnh Load Test |
                  |    - Dự báo suy giảm sức chịu tải tại 250% P_TK    |
                  +---------------------------------------------------+

Khoảng trống nghiên cứu (Research Gap) lớn nhất hiện nay nằm ở sự thiếu hụt các công cụ phân tích tích hợp giữa mô hình toán học dẻo phi tuyến và thực nghiệm vật lý kiểm chứng đa tỉ lệ. Phần lớn các phương pháp thiết kế nền móng truyền thống chỉ áp dụng lý thuyết tĩnh quy chuẩn kết hợp với hệ số an toàn tĩnh ($FS = 2.0 \div 2.5$) hoặc nhân thêm các hệ số chiết giảm động lực học kinh nghiệm mang tính phỏng đoán, chưa phản ánh đúng quy luật phân bố lại lực dọc trục, sự huy động ma sát thành cọc ($R(z)$) và sức kháng mũi cọc ($q_p$) dưới tác động của dải tần số kích thích biến thiên. Nhằm giải quyết triệt để vấn đề này, luận án xác lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế phân bố ứng suất - biến dạng và sự suy giảm sức kháng ma sát đơn vị dọc thân cọc thay đổi như thế nào theo các dải tần số dao động ($f = 10 \div 50\text{ Hz}$) và tỷ lệ hình học $L/D$?
  • Câu hỏi nghiên cứu 2 (RQ2): Mô hình cấu tạo đất dẻo tái bền nào (Modified Cam-Clay, Soft Soil hay Hardening Soil) mô phỏng chính xác nhất lộ trình ứng suất hữu hiệu và biến dạng dẻo cực hạn xung quanh thân và mũi cọc so với thí nghiệm nén tĩnh hiện trường có gắn cảm biến biến dạng (Strain gauges)?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật tương quan định lượng giữa lực dọc trục, sức kháng mũi, ma sát bên đơn vị và độ lún phá hoại đối với các bề mặt cọc có độ nhám khác nhau dưới tải trọng động tuần hoàn được xác lập ra sao?
  • Giả thuyết khoa học 1 (H1): Dưới tác động của tải trọng động tuần hoàn, sức chịu tải cực hạn của cọc suy giảm phi tuyến tính theo tần số kích thích, đạt mức suy giảm cực đại khi tần số cưỡng bức tiệm cận tần số dao động tự nhiên của hệ cọc - đất.
  • Giả thuyết khoa học 2 (H2): Độ nhám bề mặt cọc chi phối trực tiếp đến gradient chuyển tiếp ứng suất cắt tiếp xúc $\tau$, trong đó cọc nhám duy trì khả năng huy động ma sát thành cao hơn và giảm thiểu tốc độ lún tích lũy so với cọc trơn.
  • Giả thuyết khoa học 3 (H3): Mô hình Modified Cam-Clay (MCC) với bộ thông số nén $\lambda$ và nở $\kappa$ hiệu chỉnh thông qua phân tích ngược từ thí nghiệm nén tĩnh đa chu kỳ cho phép dự báo tiệm cận trạng thái phá hoại cực hạn của cọc ở cấp tải $250% P_{TK}$ với độ lún giới hạn quy ước $U_y = 80\text{ mm}$.

Luận án thiết lập khung lý thuyết (Theoretical Framework) dựa trên sự phối hợp giữa Cơ học đất trạng thái tới hạn (Critical State Soil Mechanics - CSSM), lý thuyết truyền sóng ứng suất đàn hồi 1D (Wave Equation Theory) và lý thuyết tương tự động lực học (Cauchy-Froude Scaling Laws). Phạm vi thực nghiệm tập trung vào cọc đơn chịu nén đúng tâm ngàm vào lớp cát mịn đến trung trạng thái chặt vừa điển hình tại TP.HCM (chỉ số $N_{SPT} = 15 \div 20$, góc ma sát trong $\varphi = 26^\circ \div 27^\circ$), sử dụng mô hình vật lý thu nhỏ tỉ lệ 1:3 kết hợp đầu đo biến dạng đa tầng (SG0, SG1, SG2) và mô phỏng số phần tử hữu hạn (FEM) trên phần mềm Plaxis. Kết quả nghiên cứu mang lại đóng góp đột phá khi định lượng hóa chính xác sự phân bố lại nội lực và đề xuất hệ phương trình hồi quy xác định độ lún theo tần số phá hoại trong thực tế kỹ thuật.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử của cọc chịu tải trọng động ghi nhận sự phát triển qua nhiều giai đoạn với các trường phái tiếp cận khác nhau. Trường phái truyền sóng ứng suất cổ điển bắt đầu từ nghiên cứu nền tảng của Smith (1960), thiết lập phương trình vi phân truyền sóng trong thanh một chiều chịu va đập dọc trục, làm cơ sở cho phương pháp phân tích động biến dạng lớn PDA (Pile Driving Analyzer) và phần mềm CAPWAP. Carroll (1963) đã chứng minh bằng thực nghiệm rằng tốc độ biến dạng cắt có ảnh hưởng mang tính quyết định đến sức chống cắt không thoát nước của đất sét; cụ thể, khi tốc độ biến dạng tăng từ $50%$ đến $425%$, sức chống cắt động có thể gia tăng đáng kể so với điều kiện tĩnh: $$c_{u(\text{dynamic})} = (1.5 \div 2.0) , c_{u(\text{static})}$$ Ngược lại, đối với đất rời bão hòa, Vesic chỉ ra rằng góc ma sát hiệu quả chịu tác động suy giảm dưới tải trọng động, với góc ma sát động nhỏ nhất được xác định theo quan hệ: $$\varphi'{(\text{động})} = \varphi'{(\text{tĩnh})} - 2^\circ$$ Đồng thời, Casagrande và Shannon khi thí nghiệm nén nở hông với tải xung thời gian ngắn ($t_L = 0.02\text{s}$) đã ghi nhận sức chịu nén nở hông gia tăng đạt tỷ lệ: $$\frac{q_{u(\text{xung})}}{q_{u(\text{tĩnh})}} = 1.5 \div 2.0$$ trong khi mô đun biến dạng tức thời tăng gấp đôi so với tải tĩnh: $$E_{u(\text{tải tức thời})} = 2 , E_{u(\text{tải tĩnh})}$$

Trong phân tích động lực học móng máy, Das và Ramana (2011) đã hoàn thiện lý thuyết hệ lò xo - khối lượng - giảm chấn (Spring-Mass-Damper System) cho dao động cưỡng bức: $$m \ddot{z} + c \dot{z} + k z = Q_0 \sin(\omega t + \beta)$$ Từ đó xác định biên độ chuyển vị cực đại và lực lớn nhất tác dụng lên nền đất khi hệ thống đạt trạng thái cộng hưởng động.

+----------------------------------------------------------------------------------------------------+
|                                  MA TRẬN ĐỊNH VỊ HỌC THUẬT                                          |
+---------------------------+-------------------------------+----------------------------------------+
| Hướng nghiên cứu / Tác giả | Phương pháp / Mô hình         | Hạn chế tồn tại & Khoảng trống học thuật|
+---------------------------+-------------------------------+----------------------------------------+
| Smith (1960), CAPWAP      | Phân tích truyền sóng 1D      | Chỉ xét tải xung tức thời va đập búa;  |
|                           | (Wave equation)               | không mô phỏng tải tuần hoàn liên tục. |
+---------------------------+-------------------------------+----------------------------------------+
| Das & Ramana (2011),      | Hệ lò xo - khối lượng         | Giả định nền đất đàn hồi tuyến tính;   |
| Ahmed (2015)              | (Spring-Mass-Damper)          | bỏ qua biến dạng dẻo và suy giảm ma sát|
+---------------------------+-------------------------------+----------------------------------------+
| Wu (2020), Tolun (2020),  | Mô phỏng số ABAQUS / ANSYS    | Thiếu số liệu kiểm chứng từ mô hình    |
| Ravishankar (2014)        | 3D thuần túy                  | vật lý có gắn cảm biến biến dạng thân. |
+---------------------------+-------------------------------+----------------------------------------+
| Luận án này               | Tích hợp Mô hình vật lý 1:3   | Giải quyết toàn diện từ cơ chế vi mô,  |
| (Nguyễn Mạnh Tường)       | (Cauchy-Froude) + Phân tích   | dải tần số 10-50Hz, độ nhám cọc đến    |
|                           | ngược FEM (Modified Cam-Clay) | hệ phương trình tương quan thực tế.    |
+---------------------------+-------------------------------+----------------------------------------+

Tranh luận học thuật sâu sắc nhất trong lĩnh vực này tồn tại giữa hai quan điểm: (1) Quan điểm đàn hồi động lực học tuyến tính (Linear Dynamic Elasticity) do các nhà cơ học công trình đại diện, giả định đất nền là môi trường liên tục bán không gian đàn hồi đẳng hướng có cản nhớt; và (2) Quan điểm cơ học đất dẻo phi tuyến trạng thái tới hạn (Nonlinear Critical State Plasticity) được bảo vệ bởi các nhà địa kỹ thuật hiện đại, khẳng định rằng dưới tác động của chu kỳ tải rung, đất nền xảy ra hiện tượng trượt hạt, biến dạng thể tích dẻo ($\varepsilon_v^p$) và suy giảm độ cứng cắt động ($G/G_{max}$).

Luận án này định vị chính xác tại giao điểm của hai trường phái, vượt lên các nghiên cứu quốc tế gần đây như Ahmed (2015)—người chỉ khảo sát đáp ứng vận tốc - tần số của móng khối mà bỏ qua sự phân bố ma sát cục bộ dọc thân cọc—hay các mô phỏng của Tolun (2020) và Wu (2020) vốn thiếu các dữ liệu thực nghiệm vật lý có kiểm soát áp lực tiếp xúc và độ nhám cọc. Bằng cách tích hợp dữ liệu nén tĩnh hiện trường có gắn cảm biến biến dạng tại TP.HCM với mô hình vật lý tỉ lệ thu nhỏ, nghiên cứu đã giải mã được quy luật biến đổi của tỷ lệ kháng ma sát bên trên sức kháng mũi ($F_s / q_p$) dưới các dải tần số kích thích khác nhau.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng ứng dụng của Lý thuyết cơ học đất trạng thái tới hạn (Critical State Soil Mechanics) do Roscoe, Schofield và Wroth (1958) khởi xướng, cụ thể hóa qua mô hình Modified Cam-Clay (MCC) của Roscoe và Burland (1968), trong việc giải bài toán tương tác động lực học cọc - nền đất cát bão hòa nước. Thay vì chỉ áp dụng MCC cho đất sét cố kết bình thường trong điều kiện tĩnh, nghiên cứu đã hiệu chỉnh mặt ngưỡng chảy dẻo elip (Elliptical Yield Surface) trong không gian ứng suất $p'-q$: $$f(p', q, p_0') = q^2 - M^2 p' (p_0' - p') = 0$$ để mô tả chính xác quá trình suy giảm ứng suất hữu hiệu trung bình $p'$ do sự tích tụ áp lực nước lỗ rỗng và biến dạng thể tích dẻo quanh thân cọc khi chịu tải chu kỳ.

                  Ứng suất lệch q
                        ^
                        |             Đường trạng thái tới hạn (CSL): q = M.p'
                        |             /
                        |            /   Mặt ngưỡng chảy dẻo MCC
                        |           /  .---''''---.
                        |          /  /            \
                        |         /  /   Lộ trình   \
                        |        /  |   ứng suất    |
                        |       /   |   p' - q       |
                        |      /     \              /
                        |     /       \            /
                        |    /         '---....---'
                        +---+-----------------------> Ứng suất hữu hiệu trung bình p'
                            0          p'_f       p'_0

Hệ thống mệnh đề lý thuyết (Theoretical Propositions) được xác lập gồm:

  • Mệnh đề 1 ($P_1$): Tỷ số tham số nén và nở cải tiến $\lambda/\kappa$ không phải là hằng số cố định mà biến thiên phụ thuộc vào chu kỳ gia tải; qua 2 chu kỳ nén tĩnh và động lực học, tỷ số $\lambda/\kappa$ phản ánh mức độ định hướng lại của cấu trúc hạt cát và sự gia tăng độ cứng đàn hồi dư thừa.
  • Mệnh đề 2 ($P_2$): Dưới tác động của sóng ứng suất dọc trục truyền trong thân cọc tuân theo phương trình vi phân: $$\frac{\partial^2 u}{\partial t^2} = v_c^2 \frac{\partial^2 u}{\partial z^2} - \frac{R(z)}{\rho A_p}$$ sức kháng cắt tiếp xúc huy động tại mặt phân giới cọc - đất ($R(z)$) bị triệt tiêu từng phần do sự phá hủy liên kết ma sát tĩnh, chuyển hóa thành ma sát trượt động lực học phụ thuộc vào vận tốc biến dạng tiếp xúc.
  • Mệnh đề 3 ($P_3$): Trạng thái phá hoại của cọc chịu tải trọng động là sự kết hợp đồng thời giữa hiện tượng hóa lỏng cục bộ tại đới biến dạng dẻo quanh thân cọc và sự vượt quá sức chịu tải giới hạn dẻo quy ước tại mũi cọc ($U_y = 80\text{ mm}$ ở cấp tải $250% P_{TK}$).

Khung phân tích độc đáo

Khung phân tích của luận án được xây dựng dựa trên sự hợp nhất của ba nền tảng lý thuyết: (1) Lý thuyết trạng thái tới hạn với thông số độ dốc đường CSL ($M$), chỉ số nén ($\lambda$), chỉ số nở ($\kappa$), (2) Lý thuyết truyền sóng ứng suất dọc trục trong thanh đàn hồi ($v_c = \sqrt{E/\rho}$), và (3) Lý thuyết phân tích thứ nguyên và tương tự động lực học Cauchy-Froude.

Điều kiện biên (Boundary Conditions) được xác lập nghiêm ngặt: Nền đất được mô hình hóa là môi trường đồng nhất nhiều lớp với các thông số địa kỹ thuật từ hố khoan thực tế tại KCN Cát Lái; đài cọc được thiết kế tuyệt đối cứng nhằm triệt tiêu hoàn toàn biến dạng uốn cục bộ của đài; liên kết cọc - đất là liên kết tiếp xúc phi tuyến có xét đến độ nhám bề mặt ($\delta = \varphi_n / \varphi'$). Khung phân tích này cho phép cô lập và làm rõ ảnh hưởng riêng biệt của tần số kích thích ($f$), tỷ lệ chiều dài trên đường kính ($L/D$), và cấu trúc vi mô bề mặt cọc đối với sức chịu tải tổng thể.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý thực chứng thực nghiệm (Positivism) kết hợp với mô hình hóa tính toán kiểm chứng (Computational Realism). Thiết kế nghiên cứu đa cấp (Multi-level Research Design) được triển khai qua ba pha liên hoàn chặt chẽ:

  • Cấp độ 1 (Hiện trường): Thu thập số liệu địa kỹ thuật thực tế từ các hố khoan sâu, thí nghiệm xuyên tiêu chuẩn (SPT), xuyên tĩnh đo áp lực nước lỗ rỗng (CPTU) và thí nghiệm thử tải tĩnh cọc khoan nhồi/cọc ép có gắn cảm biến biến dạng (Strain gauges) tại TP.HCM.
  • Cấp độ 2 (Phần tử hữu hạn FEM): Mô phỏng số 2D/3D đối xứng trục trên phần mềm Plaxis, thực hiện phân tích ngược (Back-analysis) để tối ưu hóa bộ thông số đầu vào cho các mô hình đất nền MCC, Hardening Soil (HS) và Soft Soil (SS).
  • Cấp độ 3 (Mô hình vật lý thu nhỏ): Chế tạo hệ thống thí nghiệm mô hình vật lý tỉ lệ hình học $\alpha_L = 1:3$, tuân thủ các tiêu chuẩn tương tự động lực học Cauchy và Froude nhằm tái hiện chính xác trường ứng suất - biến dạng dưới tác động của tải trọng rung động tuần hoàn.
                    QUY TRÌNH NGHIÊN CƯÚ RIGOROUS ĐA CẤP
                    
+-------------------------------------------------------------------------------+
| PHA 1: THỰC NGHIỆM HIỆN TRƯỜNG TẠI TP.HCM                                      |
| - Khảo sát địa chất KCN Cát Lái, Hiệp Phước, Tân Tạo (SPT, CPTU, Hố khoan)    |
| - Thí nghiệm Load Test nén tĩnh 2 chu kỳ cọc thực tế                          |
| - Đo biến dạng thân cọc bằng Strain Gauges vi mô                               |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
| PHA 2: PHÂN TÍCH SỐ PHẦN TỬ HỮU HẠN (PLAXIS FEM)                              |
| - Xây dựng mô hình 2D đối xứng trục: MCC, Hardening Soil (HS), Soft Soil (SS) |
| - Phân tích ngược (Back-analysis) xác định bộ thông số tối ưu: M, lambda, kappa|
| - Mô phỏng nén phá hoại 250% P_TK và kiểm chứng chuyển vị giới hạn Uy = 80mm  |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
| PHA 3: MÔ HÌNH VẬT LÝ TỈ LỆ 1:3 (CAUCHY - FROUDE)                             |
| - Chế tạo thùng chứa đất, hệ khung phản lực, kích thủy lực & động cơ rung     |
| - Thí nghiệm chuỗi cọc L40, L50, L60 (bề mặt Trơn & Nhám)                     |
| - Đo biến dạng đa điểm (SG0, SG1, SG2), Gia tốc kế, Phân tích phổ FFT         |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
| ĐẦU RA: HỆ PHƯƠNG TRÌNH TƯƠNG QUAN LỰC - BIẾN DẠNG - ĐỘ LÚN - TẦN SỐ PHÁ HOẠI |
+-------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Hệ thống mô hình vật lý được chế tạo với độ chính xác cơ khí cao, bao gồm:

  • Thùng chứa đất thí nghiệm: Kích thước hình học được tính toán nhằm triệt tiêu hoàn toàn hiệu ứng phản xạ sóng từ thành thùng (Boundary reflection effects), chứa mẫu cát hạt mịn đến trung được đầm nén kiểm soát chặt chẽ dung trọng và độ ẩm.
  • Hệ thống gia tải động: Sử dụng động cơ tạo rung điều khiển vô cấp bằng biến tần kỹ thuật số, cho phép tạo ra tải trọng điều hòa dạng hình sin $Q(t) = Q_0 \sin(\omega t)$ với dải tần số quét chính xác từ $10\text{ Hz}$ đến $50\text{ Hz}$.
  • Hệ thống đo lường tích hợp: Bao gồm Loadcell đo tải trọng trục, biến trở tuyến tính (LVDT) đo chuyển vị lún đầu cọc, gia tốc kế áp điện gắn trên đài cọc và chuỗi cảm biến biến dạng vi mô (Strain gauges) gắn âm dọc theo chiều dài cọc tại các vị trí: đỉnh cọc (SG0), thân cọc (SG1) và mũi cọc (SG2). Tín hiệu được thu thập qua bộ xử lý đa kênh tốc độ cao kết nối máy tính.
  • Quy trình chuẩn hóa: Các cảm biến được hiệu chuẩn tuyến tính ($R^2 > 0.999$). Dữ liệu gia tốc và lực được phân tích miền tần số thông qua biến đổi Fourier nhanh (Fast Fourier Transform - FFT) để nhận diện chính xác tần số kích thích và đáp ứng phổ của hệ thống.
                    HỆ THỐNG THÍ NGHIỆM MÔ HÌNH VẬT LÝ TỈ LỆ
                    
                     +-----------------------------------+
                     |      ĐỘNG CƠ TẠO RUNG ĐIỀU TẦN    |
                     |         (Biến tần: 10 - 50 Hz)    |
                     +-----------------+-----------------+
                                       |
                                       v
                     +-----------------+-----------------+
                     |   LOADCELL + GIA TỐC KẾ ĐÀI CỌC   |
                     +-----------------+-----------------+
                                       |
                     +-----------------+-----------------+
    ====== Mặt đất   |       [SG0] Cảm biến đỉnh cọc     |
    |                +-----------------+-----------------+
    |                                  |
    | Cát mịn                          | Thân cọc L40/L50/L60
    | chặt vừa                         | (Đường kính D)
    | (gamma, e, Dr) +-----------------+-----------------+
    |                |       [SG1] Cảm biến thân cọc     |
    |                +-----------------+-----------------+
    |                                  |
    |                                  |
    |                +-----------------+-----------------+
    v Đáy thùng      |       [SG2] Cảm biến mũi cọc      |
                     +-----------------------------------+

Data và phân tích

Phân tích số phần tử hữu hạn được tiến hành đồng thời với việc khảo sát các chuỗi mẫu cọc có chiều dài thí nghiệm khác nhau: L40 ($L/D = 13.3$), L50 ($L/D = 16.7$), L60 ($L/D = 20.0$), kết hợp hai trạng thái bề mặt: cọc trơn và cọc nhám.

Trong mô phỏng FEM (Plaxis), bộ thông số mô hình Modified Cam-Clay (MCC) của nền đất được xác lập sau quá trình phân tích ngược từ thí nghiệm nén tĩnh 2 chu kỳ tải trọng:

  • Chỉ số nén $\lambda = 0.105 \div 0.135$
  • Chỉ số nở $\kappa = 0.018 \div 0.025$
  • Độ dốc đường trạng thái tới hạn $M = 1.05 \div 1.15$ (tương ứng góc ma sát $\varphi' \approx 26^\circ \div 28^\circ$)
  • Hệ số quá cố kết $OCR = 1.0 \div 1.5$

Quá trình mô phỏng nén phá hoại được đẩy lên mức tải trọng $250% P_{TK}$ với tiêu chí độ lún phá hoại quy ước $U_y = 80\text{ mm}$. Dữ liệu từ thực nghiệm mô hình và mô phỏng số được xử lý thống kê, lọc nhiễu tín hiệu số (Digital Filtering) và xây dựng các đường cong tương quan phi tuyến với độ tin cậy thống kê cao ($p\text{-value} < 0.01$).

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Sự dịch chuyển cơ chế truyền lực dọc trục và tái phân bố sức kháng theo tần số: Dưới tác động của tải trọng động tuần hoàn, sự phân bố lực dọc trong thân cọc biến đổi rõ rệt theo chiều sâu và tần số kích thích. Ở dải tần số thấp ($f < 20\text{ Hz}$), ma sát thân cọc đoạn trên (đo bởi SG0 và SG1) tiếp nhận phần lớn tải trọng truyền xuống ($65-75%$). Tuy nhiên, khi tần số tăng cao tiệm cận dải phá hoại ($f = 22 \div 28\text{ Hz}$ tùy loại cọc), lực ma sát thành cọc bị suy giảm đột ngột do hiện tượng trượt dẻo tiếp xúc, dẫn đến sóng ứng suất truyền thẳng xuống mũi cọc, làm gia tăng sức kháng mũi tức thời $q_p$ lên gấp $1.8-2.2$ lần trước khi toàn bộ hệ thống bị phá hoại do lún tụt.
  2. Ảnh hưởng quyết định của độ nhám bề mặt cọc: Kết quả thí nghiệm so sánh trực tiếp giữa chuỗi cọc nhám và cọc trơn (L40, L50, L60) chứng minh rằng cọc có bề mặt nhám duy trì sức kháng ma sát đơn vị ($F_s$) cao hơn từ $28%$ đến $38%$ so với cọc trơn ở cùng cấp tần số. Bề mặt nhám giúp mở rộng đới biến dạng dẻo cắt của các hạt cát xung quanh cọc, ngăn chặn hiện tượng trượt phẳng trơn trượt và nâng cao tần số phá hoại giới hạn thêm $3 - 5\text{ Hz}$.
  3. Ưu thế vượt trội của mô hình Modified Cam-Clay (MCC) trong mô phỏng nén tĩnh: So sánh kết quả mô phỏng số trên Plaxis với đường cong nén tĩnh hiện trường Load Test 2 chu kỳ cho thấy: Mô hình Soft Soil (SS) đánh giá quá cao độ lún dẻo ở các cấp tải ban đầu (sai số $> 25%$), mô hình Hardening Soil (HS) cho độ cứng đàn hồi quá mức ở chu kỳ dỡ tải (CK2_100%), trong khi mô hình Modified Cam-Clay (MCC) sau khi phân tích ngược cho đường cong quan hệ Tải trọng - Độ lún ($P-S$) trùng khít nhất với sai số biên dưới $6.5%$. Đặc biệt, tỷ số thông số $\lambda/\kappa$ được xác lập chính xác cho từng chu kỳ: Chu kỳ 1 ($\lambda/\kappa \approx 5.8$), Chu kỳ 2 ($\lambda/\kappa \approx 4.2$).
  4. Xác lập hệ phương trình tương quan thực nghiệm Độ lún - Tần số: Luận án đã thiết lập thành công hệ phương trình tương quan phi tuyến giữa độ lún tích lũy ($S$) và tần số rung động ($f$) cho cọc nguyên mẫu ngoài thực tế, phân tách rõ ràng theo độ nhám và tỷ lệ $L/D$:
  • Đối với cọc trơn ($L/D = 13.3 \div 20$): $$S = A_1 \cdot \exp\left(\frac{f}{B_1}\right) + S_0$$
  • Đối với cọc nhám ($L/D = 13.3 \div 20$): $$S = A_2 \cdot f^{B_2} + C_2$$ (trong đó $A_1, B_1, A_2, B_2, C_2$ là các hệ số thực nghiệm được chuẩn hóa theo tỷ lệ Cauchy-Froude).
   Độ lún tích lũy S (mm)
        ^                                      / Phá hoại lún tụt
        |                                     /  (Cọc trơn: f_cr = 20-22 Hz)
        |                                    /   (Cọc nhám: f_cr = 25-28 Hz)
        |                                   /
        |                       Cọc trơn   /     Cọc nhám
        |                             .---'         .---'
        |                           .'            .'
        |                         .'            .'
        |                     _.-'          _.-'
        |                 _.-'          _.-'
        +----------------+-------------+---------------------> Tần số f (Hz)
        0                10            20           30

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Nghiên cứu đã chứng minh tính khả thi của việc tích hợp lý thuyết trạng thái tới hạn (CSL/MCC) vào bài toán động lực học cọc, cung cấp bằng chứng thực nghiệm khẳng định sự thay đổi tỷ số ma sát bên/sức kháng mũi ($F_s / q_p$) là một hàm phụ thuộc phi tuyến vào tần số dao động cưỡng bức.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn mực kết hợp giữa phương pháp phân tích thứ nguyên Cauchy-Froude trong thí nghiệm mô hình vật lý và phân tích ngược số học FEM Plaxis, mở ra phương pháp tiếp cận kinh tế và tin cậy thay thế cho các thí nghiệm hiện trường quy mô lớn tốn kém.
  • Về mặt thực tiễn kỹ thuật: Cung cấp cho các kỹ sư thiết kế nền móng tại TP.HCM công cụ định lượng để tính toán sức chịu tải cọc đỡ móng máy, bệ máy turbin, trạm ép nén khí trong các khu công nghiệp, ngăn ngừa rủi ro phá hoại công trình do cộng hưởng dao động hoặc hóa lỏng cục bộ.
  • Về mặt chính sách và quy chuẩn: Đóng góp luận cứ khoa học thực nghiệm vững chắc để kiến nghị bổ sung các điều khoản kiểm soát biến dạng động và tần số giới hạn vào Tiêu chuẩn thiết kế móng cọc xây dựng tại vùng đất yếu đô thị.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu mang tính khách quan:

  1. Môi trường đất thí nghiệm: Nghiên cứu mô hình vật lý tập trung chủ yếu vào môi trường cát mịn đến trung trạng thái chặt vừa đồng nhất. Các hiện tượng tương tác phức tạp trong các tầng sét bùn xen kẹp hữu cơ dày (như lớp bùn sét bão hòa nước $e > 1.8$ đặc trưng tại Bình Chánh và Nhà Bè) chưa được mô phỏng toàn diện trong thùng thí nghiệm động.
  2. Cấu hình liên kết đài - cọc: Thí nghiệm mô hình sử dụng đài cọc tuyệt đối cứng và cọc đơn chịu tải dọc trục thẳng đứng, bỏ qua tương tác đàn dẻo của đài cọc mềm và thành phần tải trọng ngang/mô men uốn đồng thời tác dụng.
  3. Hiệu ứng nhóm cọc (Pile Group Effect): Nghiên cứu chưa khảo sát sự giao thoa sóng ứng suất và hiệu ứng bóng mờ (Shadowing effect) giữa các cọc trong nhóm cọc chịu tải trọng động đa nguồn.
  4. Ảnh hưởng lâu dài của hiện tượng từ biến (Creep): Tải trọng động thí nghiệm được thực hiện trong khoảng thời gian chu kỳ xác định, chưa xét đến biến dạng từ biến tích lũy qua nhiều năm khai thác công trình móng máy.

Chương trình nghiên cứu tương lai (Future Research Agenda) trong giai đoạn 5 - 10 năm tới cần tập trung vào:

  • Mở rộng mô hình vật lý và mô phỏng số cho hệ thống nhóm cọc ($2 \times 2, 3 \times 3$) trong môi trường đất nhiều lớp có xét đến tải trọng động lệch tâm và tải trọng ngang dao động.
  • Ứng dụng phương pháp phần tử rời rạc (Discrete Element Method - DEM) kết hợp FEM (DEM-FEM Coupling) để mô phỏng sự vỡ vụn hạt cát và tái cấu trúc vi mô tại mặt tiếp xúc cọc - đất dưới hàng triệu chu kỳ tải rung.
  • Nghiên cứu giải pháp gia cố cải tạo bề mặt cọc (như cọc có gờ, cọc mở rộng thân) nhằm triệt tiêu sóng ứng suất truyền xuống mũi và nâng cao tần số cộng hưởng an toàn cho công trình.

Tác động và ảnh hưởng

Luận án tạo ra những tác động đa diện và sâu rộng đối với cả giới học thuật lẫn thực tiễn sản xuất công nghiệp:

  • Tác động học thuật (Academic Impact): Đóng góp hệ thống dữ liệu thực nghiệm chất lượng cao về tương tác động lực học cọc - đất cát tại vùng hạ lưu lưu vực sông Đồng Nai - Sài Gòn. Các phương trình tương quan và dữ liệu phân tích phổ FFT dự kiến sẽ trở thành tài liệu tham khảo chuẩn mực, thúc đẩy các công bố quốc tế trên các tạp chí chuyên ngành địa kỹ thuật uy tín (như Computers and Geotechnics, Soil Dynamics and Earthquake Engineering).
  • Chuyển đổi công nghiệp (Industry Transformation): Ứng dụng trực tiếp cho các tổng thầu xây dựng, tư vấn thiết kế nền móng trong việc tối ưu hóa chiều dài cọc $L/D$, lựa chọn chủng loại cọc (cọc ly tâm dự ứng lực có gờ nhám thay vì cọc trơn), giúp giảm thiểu từ $15%$ đến $20%$ chi phí xử lý nền móng dư thừa do sử dụng hệ số an toàn tĩnh mù mờ.
  • Lợi ích kinh tế - xã hội: Bảo vệ an toàn tuyệt đối cho các tổ hợp nhà máy công nghiệp nặng, trạm biến áp, xưởng sản xuất công nghệ cao tại các KCN Cát Lái, Tân Tạo, Hiệp Phước; ngăn chặn nguy cơ lún nứt công trình lân cận do tải rung động lan truyền qua nền đất yếu đô thị.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật Địa kỹ thuật: Tiếp cận phương pháp luận mẫu mực về phân tích thứ nguyên Cauchy-Froude trong mô hình vật lý thu nhỏ, kết hợp phương pháp phân tích ngược tham số mô hình đất trạng thái tới hạn Modified Cam-Clay trên Plaxis.
  • Kỹ sư Thiết kế Kết cấu - Nền móng: Sở hữu bộ công thức tương quan thực nghiệm trực tiếp giữa Lực dọc - Độ lún - Tần số dao động ($f$) để tính toán chính xác sức chịu tải khả dụng của móng cọc chịu tải trọng máy rung động tuần hoàn.
  • Các Ban Quản lý Dự án & Doanh nghiệp Khai thác KCN: Có cơ sở khoa học để thẩm tra, phê duyệt biện pháp thi công và thiết kế móng thiết bị, kiểm soát rủi ro lún nứt công trình trong suốt vòng đời dự án.
  • Cơ quan Quản lý Nhà nước về Xây dựng: Nắm bắt luận cứ kỹ thuật để hoàn thiện quy chuẩn xây dựng, tiêu chuẩn kiểm định chất lượng móng cọc chịu tải trọng động tại các đô thị đặc thù trên nền đất yếu.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo và đột phá nhất của luận án là gì?
Đó là việc mở rộng thành công phạm vi ứng dụng của Mô hình cơ học đất trạng thái tới hạn Modified Cam-Clay (MCC) của Roscoe & Burland sang bài toán động lực học tương tác cọc - đất cát chịu tải chu kỳ. Luận án đã lượng hóa sự biến thiên của tỷ số tham số nén - nở $\lambda/\kappa$ qua các chu kỳ tải trọng nén tĩnh và xác lập được lộ trình ứng suất hữu hiệu $p'-q$ phản ánh chính xác trạng thái suy giảm ứng suất do gia tăng áp lực nước lỗ rỗng và biến dạng dẻo cực hạn ở cấp tải phá hoại $250% P_{TK}$.

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu quốc tế trước đây?
So với các nghiên cứu của Das & Ramana (2011) vốn chỉ dừng lại ở hệ lò xo đàn hồi tuyến tính đơn giản, hay các mô phỏng thuần túy số học của Tolun (2020) và Wu (2020), luận án đã thiết lập một quy trình thực nghiệm vật lý tỉ lệ 1:3 hoàn chỉnh tuân thủ nghiêm ngặt các định luật tương tự động lực học Cauchy-Froude. Hệ thống này được tích hợp cảm biến đo biến dạng đa tầng (SG0, SG1, SG2), gia tốc kế đài cọc và xử lý tín hiệu phổ FFT, cho phép đo đạc trực tiếp sự dịch chuyển của ma sát thành cọc và sức kháng mũi theo từng bước tần số $10 - 50\text{ Hz}$.

3. Phát hiện thực nghiệm nào mang tính bất ngờ và có giá trị khoa học cao nhất?
Phát hiện về sự biến đổi đột ngột của tỷ lệ ma sát bên trên sức kháng mũi ($F_s / q_p$) tại tần số cộng hưởng phá hoại: Ở trạng thái tĩnh và tần số thấp ($f < 20\text{ Hz}$), ma sát thân chịu phần lớn tải trọng, nhưng khi tần số kích thích đạt ngưỡng tới hạn ($22 - 28\text{ Hz}$), ma sát thân bị trượt dẻo và triệt tiêu nhanh chóng, dồn toàn bộ sóng ứng suất xung kích xuống mũi cọc, làm sức kháng mũi tức thời vọt tăng cục bộ trước khi đất nền quanh mũi bị phá hủy dẻo hoàn toàn.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không?
Có. Luận án trình bày chi tiết toàn bộ các bước thiết kế thứ nguyên, tỷ lệ chuyển đổi Cauchy-Froude cho bê tông cốt thép tỉ lệ 1:3, thông số cơ lý cát chuẩn, quy trình chế tạo thùng chứa đất, hệ khung phản lực, sơ đồ đấu nối mạch cầu strain gauge, quy trình hiệu chuẩn Loadcell và phương pháp thiết lập các phase phân tích trong phần mềm Plaxis FEM.

5. Định hướng nghiên cứu 10 năm tới được vạch ra như thế nào?
Định hướng tập trung vào việc mở rộng bài toán sang hệ nhóm cọc không gian 3D tương tác đa hướng, tích hợp mô hình hóa liên kết hạt DEM-FEM để giải mã cơ chế vi mô của sự suy giảm góc ma sát tiếp xúc $\varphi_n$, và phát triển các giải pháp kết cấu cọc thông minh giảm chấn tự động cho các siêu đô thị ven biển.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Mạnh Tường đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra, tạo nên dấu ấn học thuật sâu sắc với 6 đóng góp cốt lõi:

  1. Xác lập cơ sở lý luận toàn diện và kiểm chứng thực nghiệm thành công sự suy giảm sức chịu tải của cọc dưới tác động của tải trọng động tuần hoàn trong môi trường đất cát yếu đặc thù TP.HCM.
  2. Thiết kế, chế tạo hoàn chỉnh hệ thống mô hình vật lý tỉ lệ 1:3 theo chuẩn tương tự Cauchy-Froude, tích hợp hệ thống thu nhận dữ liệu biến dạng vi mô đa điểm dọc thân cọc.
  3. Chứng minh mô hình Modified Cam-Clay (MCC) sau khi phân tích ngược từ thí nghiệm nén tĩnh 2 chu kỳ là mô hình tối ưu nhất để mô phỏng trạng thái ứng suất - biến dạng dẻo và dự báo tải trọng phá hoại cực hạn của cọc so với các mô hình Hardening Soil và Soft Soil.
  4. Phát hiện quy luật phân bố lại nội lực dọc trục và sự chuyển dịch tỷ lệ huy động ma sát bên/sức kháng mũi theo các dải tần số kích thích từ $10\text{ Hz}$ đến $50\text{ Hz}$.
  5. Định lượng hóa ưu thế vượt trội của bề mặt cọc nhám trong việc nâng cao sức kháng cắt ma sát bên và kéo giãn giới hạn tần số phá hoại lún tụt so với cọc trơn.
  6. Đề xuất hệ phương trình hồi quy thực nghiệm xác định mối quan hệ Độ lún - Tần số phục vụ trực tiếp cho công tác tính toán, thiết kế và tối ưu hóa hệ móng cọc công trình công nghiệp trong thực tiễn xây dựng.