Tổng quan về luận án

Nghiên cứu về sự tương tác động lực học giữa kết cấu móng cọc và môi trường đất nền dưới tác động của tải trọng ngang và sóng địa chấn là một trong những bài toán phức tạp bậc nhất của cơ học địa kỹ thuật hiện đại. Trong bối cảnh Việt Nam đẩy mạnh xây dựng các công trình hạ tầng quy mô lớn như nhà máy điện hạt nhân, công trình biển, đập thủy điện, cảng biển và các tổ hợp nhà cao tầng, nguy cơ rủi ro địa chấn từ các đới đứt gãy hoạt động mạnh (như đới đứt gãy Sông Mã với trận động đất Điện Biên năm 1935 đạt $M = 6.7$ Richter, đứt gãy Sơn La với trận động đất Tuần Giáo năm 1983 đạt $M = 6.8$ Richter, hay hiện tượng động đất kích thích tại khu vực thủy điện Sông Tranh 2) đặt ra yêu cầu cấp thiết phải hoàn thiện cơ sở lý thuyết và công cụ tính toán móng cọc chịu tải trọng phức tạp.

                  +-------------------------------------------------------+
                  |  TẢI TRỌNG TĨNH NGANG VÀ ĐỘNG ĐẤT (El Centro 1940)    |
                  +-------------------------------------------------------+
                                              |
                                              v
           +---------------------------------------------------------------------+
           |      NGUYÊN LÝ CỰC TRỊ GAUSS (GAUSS'S PRINCIPLE OF LEAST CONSTRAINT) |
           +---------------------------------------------------------------------+
                        /                                     \
                       /                                       \
                      v                                         v
+-------------------------------------------+ +-------------------------------------------+
| KHỐI ĐẤT NỀN ĐÀN HỒI LIÊN TỤC 3D          | | KẾT CẤU CỌC ĐƠN / DẦM TIMOSHENKO          |
| - Nửa không gian / Không gian vô hạn      | | - Xét biến dạng trượt ngang gamma_c       |
| - Phần tử hữu hạn khối 20 nút             | | - Góc xoay phi_c và độ võng u_c           |
| - Ma trận độ cứng và ma trận cản          | | - Nhân tử Lagrange liên kết biên lambda(z)|
+-------------------------------------------+ +-------------------------------------------+
                      \                                         /
                       \                                       /
                        v                                     v
           +---------------------------------------------------------------------+
           |      PHIẾM HÀM LƯỢNG CƯỠNG BỨC TỔNG THỂ: Z = Zd + Zc + ZAB = min     |
           +---------------------------------------------------------------------+
                                              |
                                              v
           +---------------------------------------------------------------------+
           | BỘ CHƯƠNG TRÌNH TÍNH TOÁN (Mstatic, Kstatic, KdynaS, KdynaP, KdynaPE)|
           | - Miền thời gian & miền tần số (Biến đổi DFT/FFT, Tích phân Duhamel)|
           | - Khảo sát sóng Love, phản ứng động học & tương tác quán tính       |
           +---------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong y văn xuất phát từ sự bất cập của hai nhóm phương pháp truyền thống: nhóm mô hình phản lực cục bộ (Winkler, đường cong $p-y$) bỏ qua tính liên tục của môi trường đất nền và hiệu ứng phân tán năng lượng sóng, trong khi nhóm mô hình phần tử hữu hạn đàn hồi liên tục 3D cổ điển đòi hỏi chi phí tính toán cực lớn và gặp khó khăn khi mô phỏng biên vô hạn trong bài toán động lực học. Luận án tiến sĩ kỹ thuật của tác giả Ngô Quốc Trinh với đề tài "Nghiên cứu sự làm việc của cọc chịu tải trọng ngang và tải trọng động đất" do PGS.TS. Vương Văn Thành và TS. Trần Hữu Hà hướng dẫn khoa học tại Trường Đại học Kiến trúc Hà Nội (2014) đã giải quyết triệt để ranh giới này.

Hệ thống câu hỏi và giả thuyết nghiên cứu được thiết lập cụ thể:

  1. Câu hỏi nghiên cứu 1 (Q1): Làm thế nào để mô hình hóa chính xác sự phân bố ứng suất và biến dạng của môi trường đất nền 3D bán vô hạn chịu tải trọng ngang mà không làm bùng nổ số bậc tự do tại các biên xa?
    • Giả thuyết 1 (H1): Ứng dụng Nguyên lý cực trị Gauss kết hợp hệ so sánh không gian đàn hồi vô hạn cho phép thiết lập phương trình biến phân chuyển vị thỏa mãn đồng thời phương trình vi phân cân bằng và các điều kiện biên vô cực.
  2. Câu hỏi nghiên cứu 2 (Q2): Biến dạng trượt ngang của cọc theo lý thuyết dầm Timoshenko ảnh hưởng như thế nào đến nội lực và chuyển vị của cọc khi làm việc đồng thời với môi trường đất nền liên tục?
    • Giả thuyết 2 (H2): Việc tích hợp lý thuyết dầm Timoshenko vào phiếm hàm năng lượng cưỡng bức Gauss thông qua nhân tử Lagrange sẽ làm tăng độ chính xác tính toán chuyển vị đầu cọc từ $8% - 14%$ so với lý thuyết dầm Bernoulli-Euler truyền thống.
  3. Câu hỏi nghiên cứu 3 (Q3): Quy luật lan truyền sóng cắt (sóng Love) và cơ chế phản ứng động lực học của cọc diễn ra như thế nào dưới tác động của phổ gia tốc động đất thực tế trong miền thời gian và miền tần số?
    • Giả thuyết 3 (H3): Phép biến đổi Fourier rời rạc (DFT/FFT) kết hợp tích phân Duhamel cho lời giải xung đơn vị cho phép bóc tách chính xác hiện tượng cộng hưởng tần số và vị trí tập trung mômen uốn cực đại trên thân cọc theo từng pha thời gian chấn động.

Khung lý thuyết của công trình tích hợp chặt chẽ giữa Nguyên lý cực trị Gauss (Gauss's Principle of Least Constraint), Lý thuyết dầm Timoshenko, Lý thuyết nửa không gian vô hạn đàn hồi (Mindlin, Boussinesq) và Lý thuyết truyền sóng trong môi trường đàn hồi (Love, Rayleigh). Phạm vi nghiên cứu bao quát các bài toán từ tải trọng tĩnh nằm ngang đặt tại đỉnh, thân và chân cọc đến tải trọng động điều hòa dải tần $0.5 - 60 \text{ Hz}$ và kích động động đất El Centro (1940) theo phương Bắc - Nam.

Literature Review và Positioning

Lịch sử nghiên cứu tương tác cọc - đất nền chịu tải trọng ngang chứng kiến sự cạnh tranh học thuật giữa hai trường phái tiếp cận cơ bản:

+----------------------------------------------------------------------------------------------------+
|                                    TIẾN TRÌNH TIẾP CẬN Y VĂN                                      |
+----------------------------------------------------------------------------------------------------+
| 1. Mô hình Phản lực Cục bộ (Winkler 1867, Zavriev 1962, Matlock 1970, Reese 1974)                  |
|    - Đất nền là các lò xo rời rạc độc lập (kh); đường cong p-y phi tuyến.                          |
|    - Ưu điểm: Đơn giản, giải nhanh.                                                                |
|    - Nhược điểm: Triệt tiêu hoàn toàn tương tác cắt và tính liên tục không gian 3D của đất nền.    |
+----------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
| 2. Mô hình Nửa không gian Đàn hồi Liên tục (Mindlin 1936, Poulos 1971, 1973, Banerjee 1978)       |
|    - Tích hợp phương trình Mindlin cho môi trường bán vô hạn đồng nhất; cọc dầm Bernoulli mỏng.   |
|    - Ưu điểm: Phản ánh tính liên tục của môi trường xung quanh.                                   |
|    - Nhược điểm: Bỏ qua biến dạng trượt của cọc; công thức tích phân biên phức tạp, khó mở rộng.  |
+----------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
| 3. Mô hình Động lực học & Trở kháng Động (Novak 1974, 1977, Gazetas 1991, Makris 1992)            |
|    - Tương tác động học (Kinematic) & Quán tính (Inertial); độ cứng phức K* = K_real + i*omega*C.  |
|    - Nhược điểm: Đa phần áp dụng cho cọc đàn hồi ngàm trong tầng đất đồng nhất lý tưởng.           |
+----------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
| 4. ĐỊNH VỊ CỦA LUẬN ÁN (Ngô Quốc Trinh, 2014)                                                      |
|    - Phương pháp Nguyên lý Cực trị Gauss (PPNLCT Gauss) hợp nhất Dầm Timoshenko + Nền Đàn hồi 3D.  |
|    - Khử biên vô hạn bằng hệ so sánh không gian đàn hồi; giải quyết tải trọng tĩnh, động, động đất.|
+----------------------------------------------------------------------------------------------------+

Trường phái thứ nhất dựa trên mô hình nền Winkler (1867), được Zavriev (1962) phát triển cho bài toán cọc trong công trình giao thông và sau đó được Matlock (1970) cùng Reese (1974) hoàn thiện dưới dạng phương pháp đường cong phản lực $p-y$. Mặc dù phương pháp $p-y$ phản ánh được tính phi tuyến cục bộ của đất nền gần mặt đất, điểm yếu cốt tử là xem các lớp đất hoạt động như những lò xo độc lập, triệt tiêu hoàn toàn ứng suất tiếp $\tau_{ij}$ và tính liên tục ba chiều của môi trường.

Trường phái thứ hai dựa trên mô hình môi trường đàn hồi liên tục, khởi xướng từ lời giải tích phân điểm đặt lực tập trung trong lòng bán không gian của Mindlin (1936), được Poulos (1971, 1973) ứng dụng để tính toán cọc chịu tải trọng ngang bằng phương pháp phân đoạn giải tích số. Tiếp nối hướng đi này, Novak (1974, 1977) và Gazetas (1991) đã mở rộng nghiên cứu sang bài toán động lực học thông qua các hàm độ cứng phức (impedance functions) bao gồm độ cứng kháng uốn $K_x$, độ cứng kháng xoay $K_\theta$, độ cứng liên hợp $K_{x-\theta}$ và các hệ số cản bức xạ $C$. Tuy nhiên, các nghiên cứu của Poulos và Novak đều giả định cọc tuân theo lý thuyết uốn thuần túy Euler-Bernoulli, hoàn toàn bỏ qua năng lượng biến dạng trượt ngang của tiết diện cọc.

Luận án của tác giả Ngô Quốc Trinh định vị chính xác tại điểm giao thoa giữa cơ học biến phân nâng cao và động lực học công trình ngầm:

  1. So sánh với nghiên cứu quốc tế của Poulos (1971): Trong khi Poulos sử dụng phương pháp giải tích tích phân biên giả định cọc mỏng đàn hồi trong nền đồng nhất, luận án giải quyết bài toán bằng Phương pháp Nguyên lý cực trị Gauss (PPNLCT Gauss) kết hợp phần tử dầm Timoshenko, cho phép khảo sát chính xác ứng xử của cả cọc ngắn, cọc cứng và cọc chịu uốn trong môi trường nhiều lớp đất có độ cứng biến thiên.
  2. So sánh với nghiên cứu thực nghiệm và số của Kim, O'Neill và Matlock (1987): Khi đối chiếu với kết quả thí nghiệm cọc chịu tải trọng ngang tăng dần từ $200 \text{ kN}$ đến $800 \text{ kN}$ của Kim, O'Neill và Matlock, phương pháp của luận án cho thấy sự tương đồng cao về dạng phân bố mômen uốn và chuyển vị ngang dọc thân cọc, đồng thời khắc phục được hiện tượng đánh giá quá cao độ võng đỉnh cọc của các mô hình lò xo Winkler cổ điển.

Đó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 việc ứng dụng Nguyên lý cực trị Gauss — vốn là một nguyên lý biến phân cơ bản của cơ học giải tích cổ điển do Carl Friedrich Gauss thiết lập năm 1829 — sang lĩnh vực cơ học môi trường liên tục và tương tác kết cấu - nền đất (Soil-Structure Interaction - SSI).

Bằng cách đưa khái niệm "Lượng cưỡng bức" (Constraint Quantity) vào trạng thái ứng suất - biến dạng của hệ cọc - đất, luận án đã mở rộng lý thuyết dầm Timoshenko (Timoshenko, 1921) để xây dựng phiếm hàm năng lượng có xét đồng thời biến dạng uốn $\chi_c$ và góc trượt ngang $\gamma_c$ của thân cọc:

$$Z_c = \frac{1}{2} \int_0^l \left[ E_c J \left( \frac{d\varphi_c}{dz} \right)^2 + \alpha G_c F \left( \frac{du_c}{dz} - \varphi_c \right)^2 \right] dz - \int_0^l q(z) u_c(z) dz$$

Trong đó:

  • $E_c, G_c$ lần lượt là môđun đàn hồi và môđun trượt của vật liệu cọc;
  • $J, F$ là mômen quán tính và diện tích tiết diện ngang của cọc;
  • $\varphi_c$ là góc quay của tiết diện do mômen uốn gây ra;
  • $\gamma_c = \frac{du_c}{dz} - \varphi_c$ là góc trượt ngang do lực cắt $Q$ sinh ra;
  • $\alpha$ là hệ số xét đến sự phân bố không đều của ứng suất tiếp trên tiết diện cọc.
+--------------------------------------------------------------------------------------------------------+
|                                     CẤU TRÚC PHIẾM HÀM CƯỠNG BỨC TỔNG THỂ                               |
+--------------------------------------------------------------------------------------------------------+
|                                              Z_total = min                                             |
|                                                    |                                                   |
|        +-------------------------------------------+------------------------------------------+        |
|        |                                           |                                          |        |
|        v                                           v                                          v        |
| +-------------------------------+         +--------------------------------+         +----------------+
| |        Zd (Khối đất)          |         |         Zc (Thân cọc)          |         |  ZAB (Biên AB) |
| | - Năng lượng biến dạng đàn hồi|         | - Biến dạng uốn (Ec, J, phi_c) |         | - Nhân tử      |
| | - Trường ứng suất liên tục    |         | - Biến dạng trượt (alpha, Gc, F|         |   Lagrange     |
| | - Phần tử khối 20 nút         |         | - Chuyển vị ngang cọc (u_c)    |         |   lambda(z)    |
+---------------------------------+         +--------------------------------+         +----------------+

Mô hình lý thuyết thiết lập mối quan hệ tương tác thông qua các mệnh đề (propositions):

  • Mệnh đề 1: Sự làm việc đồng thời giữa cọc và đất nền được bảo đảm tuyệt đối tại mặt tiếp xúc thông qua điều kiện liên tục chuyển vị $u_c(z) = u_d(z)$ gắn kèm thừa số Lagrange $\lambda(z)$.
  • Mệnh đề 2: Phiếm hàm lượng cưỡng bức tổng thể của hệ đạt giá trị cực tiểu thực sự khi và chỉ khi trường chuyển vị và trường ứng suất thỏa mãn phương trình vi phân cân bằng đàn hồi của Navier và phương trình truyền sóng Helmholtz.

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

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng: (1) Nguyên lý cực trị Gauss, (2) Lý thuyết dầm Timoshenko chịu uốn - trượt, và (3) Lý thuyết bán không gian đàn hồi vô hạn.

Điểm đột phá phương pháp luận thể hiện ở kỹ thuật "Hệ so sánh" (Comparison System). Thay vì phải xử lý trực tiếp điều kiện biên phức tạp của bán không gian vô hạn với vô số ẩn số biên, tác giả thiết lập bài toán tương tác thông qua một hệ so sánh là không gian vô hạn đàn hồi đã biết nghiệm giải tích cơ bản (lời giải Kelvin/Mindlin). Phiếm hàm cưỡng bức tổng thể của toàn hệ thống được biểu diễn tường minh:

$$Z = Z_d + Z_c + Z_{AB}$$

Trong đó $Z_d$ là lượng cưỡng bức của khối đất, $Z_c$ là lượng cưỡng bức của cọc, và $Z_{AB}$ là lượng cưỡng bức trên bề mặt tiếp xúc $AB$ có xét đến hàm phạt Lagrange:

$$Z_{AB} = \int_{\Omega_{AB}} \lambda(z) \left[ u_c(z) - u_d(x,y,z) \right] d\Omega$$

Điều kiện biên xác định (boundary conditions): Cọc ngàm hoặc tự do tại đỉnh ($z = 0$), chân cọc tự do hoặc tựa trên lớp đá cứng ($z = l$); môi trường đất nền mở rộng ra vô cực thỏa mãn điều kiện triệt tiêu ứng suất và chuyển vị tại vô tận ($\sigma_{ij} \to 0, u_i \to 0$ khi $r \to \infty$).

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng (Positivism) kết hợp triết lý cơ học tính toán diễn dịch (Deductive Computational Mechanics). Thiết kế nghiên cứu là sự kết hợp đa tầng giữa:

  1. Giải tích toán học thuần túy (biến phân nguyên lý Gauss);
  2. Rời rạc hóa số học bằng Phương pháp Phần tử Hữu hạn (FEM);
  3. Xử lý tín hiệu miền thời gian - tần số nâng cao (DFT/FFT và tích phân Duhamel).
+--------------------------------------------------------------------------------------------------------+
|                                   THIẾT KẾ PHƯƠNG PHÁP NGHIÊN CỨU ĐA TẦNG                              |
+--------------------------------------------------------------------------------------------------------+
|                                    TẦNG 1: TOÁN CƠ HỌC BIẾN PHÂN                                       |
|  - Phương trình vi phân Navier-Cauchy + Helmholtz.                                                     |
|  - Nguyên lý cực trị Gauss: Phiếm hàm lượng cưỡng bức Z.                                               |
+--------------------------------------------------------------------------------------------------------+
                                                    |
                                                    v
|                                    TẦNG 2: MÔ HÌNH HÓA RỜI RẠC SỐ                                      |
|  - Không gian đất 3D: Phần tử đẳng tham số Hexahedral 20 nút (Serendipity).                            |
|  - Kết cấu cọc: Phần tử dầm Timoshenko bậc cao (xét biến dạng trượt ngang gamma_c).                    |
|  - Liên kết tương tác: Kỹ thuật nhân tử Lagrange lambda(z).                                            |
+--------------------------------------------------------------------------------------------------------+
                                                    |
                                                    v
|                                    TẦNG 3: ĐỘNG LỰC HỌC CÔNG TRÌNH                                     |
|  - Gia tốc đồ El Centro 1940 (PGA = 0.348g, dt = 0.02s).                                               |
|  - Biến đổi DFT / FFT <---> IFFT; Tần số Nyquist f_N; Tích phân Duhamel.                               |
|  - Khảo sát sóng Love (G1/rho1 < G2/rho2) & Giảm chấn vật liệu zeta_h.                                 |
+--------------------------------------------------------------------------------------------------------+

Môi trường đất nền 3D được mô hình hóa bằng các phần tử khối chữ nhật 20 nút (20-node isoparametric serendipity brick elements) với 3 bậc tự do chuyển vị $(u, v, w)$ tại mỗi nút, cho phép mô phỏng chính xác gradien ứng suất phi tuyến cục bộ xung quanh thân cọc.

Quy trình nghiên cứu rigorous

Quy trình giải thuật số được kiểm chứng qua hệ thống 9 chương trình phần mềm chuyên dụng do tác giả tự lập trình bằng ngôn ngữ tính toán kỹ thuật:

+---------------------------------------------------------------------------------------------------------+
|                                 HỆ THỐNG 9 BỘ CHƯƠNG TRÌNH TÍNH TOÁN CỦA LUẬN ÁN                        |
+-----------------------------------+---------------------------------------------------------------------+
| Nhóm Tĩnh học Đất nền             | - Mstatic1: Tính khối đất chịu lực tĩnh ngang (Hệ SS bán vô hạn)    |
|                                   | - Kstatic1: Tính khối đất chịu lực tĩnh ngang (Hệ SS không gian vô hạn) |
+-----------------------------------+---------------------------------------------------------------------+
| Nhóm Tĩnh học Tương tác Cọc - Đất | - MstaticP1: Tính cọc đơn chịu tải tĩnh ngang (Hệ SS bán vô hạn)    |
|                                   | - KstaticP1: Tính cọc đơn chịu tải tĩnh ngang (Hệ SS vô hạn)        |
|                                   | - KstaticPLs: Tính cọc trong nền đàn hồi nhiều lớp (Multilayered)   |
+-----------------------------------+---------------------------------------------------------------------+
| Nhóm Động lực học & Địa chấn      | - KdynaS: Khảo sát dao động khối đất chịu tải trọng động điều hòa   |
|                                   | - KdynaL: Khảo sát sự hình thành và truyền sóng cắt (Sóng Love)     |
|                                   | - KdynaP: Tính cọc đơn chịu tải trọng động ngang                    |
|                                   | - KdynaPE: Tính cọc chịu tác động của sóng động đất (El Centro)    |
+-----------------------------------+---------------------------------------------------------------------+

Độ tin cậy (reliability) và tính hợp lệ xây dựng (construct validity) của thuật toán được bảo đảm thông qua quy trình kiểm chứng chéo ba bước:

  • Bước 1 (Cross-verification): So sánh kết quả chuyển vị giữa chương trình Mstatic1 (hệ so sánh nửa không gian vô hạn) và Kstatic1 (hệ so sánh không gian vô hạn), sai số tương đối tuyệt đối luôn nhỏ hơn $1.2%$.
  • Bước 2 (Benchmark so sánh): So sánh với lời giải giải tích của Poulos (1971) và Zavriev (1962). Kết quả tính chuyển vị lớn nhất đầu cọc theo PPNLCT Gauss có độ tương thích cao ($R^2 > 0.98$).
  • Bước 3 (Thực nghiệm kiểm chứng): Kiểm tra đối chứng với dữ liệu thí nghiệm hiện trường của Kim, O'Neill và Matlock (1987) dưới các cấp tải trọng ngang $200 \text{ kN}, 400 \text{ kN}, 600 \text{ kN}, 800 \text{ kN}$.

Data và phân tích

Phân tích động lực học công trình sử dụng chuỗi số liệu gia tốc đồ chuẩn của trận động đất El Centro (California, 1940) ghi nhận theo phương Bắc - Nam:

  • Độ dài bản ghi thời gian: $T = 30.0 \text{ s}$;
  • Bước rời rạc hóa thời gian: $\Delta t = 0.02 \text{ s}$;
  • Số lượng điểm dữ liệu: $N = 1500$ điểm;
  • Gia tốc cực đại nền (Peak Ground Acceleration - PGA): $a_{\max} = 0.348g \approx 3.41 \text{ m/s}^2$.

Quá trình chuyển đổi tín hiệu thực hiện qua phép biến đổi Fourier rời rạc (DFT) từ miền thời gian sang miền tần số, xác định dải tần số khảo sát $f \in [0.5, 60 \text{ Hz}]$ với tần số Nyquist $f_N = \frac{1}{2\Delta t} = 25 \text{ Hz}$.

Nghiên cứu áp dụng tích phân Duhamel trong miền tần số để giải phương trình dao động cưỡng bức có cản nhớt:

$$y(t) = \int_0^t p_{eff}(\tau) h(t - \tau) d\tau = \mathcal{F}^{-1} \left{ P_{eff}(\omega) \cdot H(\omega) \right}$$

Trong đó $H(\omega)$ là hàm truyền tần số của hệ cọc - đất có xét hệ số giảm chấn vật liệu $\zeta_h = 0.05$ ($5%$ cản quy ước của nền đất).

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

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

Nghiên cứu mang lại 5 phát hiện cốt lõi có giá trị khoa học cao:

+--------------------------------------------------------------------------------------------------------+
|                                    TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                          |
+--------------------------------------------------------------------------------------------------------+
| 1. HIỆU ỨNG TẬP TRUNG ỨNG SUẤT VÙNG NÔNG DƯỚI TẢI TĨNH                                                |
|    - 80% - 90% tải trọng ngang bị tiêu tán trong phạm vi độ sâu hữu hiệu L_e = (3 - 5)d.               |
|    - Mômen uốn cực đại M_max xuất hiện tại độ sâu z = (1.5 - 2.5)d.                                   |
+--------------------------------------------------------------------------------------------------------+
| 2. ĐÓNG GÓP CỦA BIẾN DẠNG TRƯỢT TIMOSHENKO                                                             |
|    - Tăng chuyển vị ngang đỉnh cọc từ 8.2% đến 14.1% so với Bernoulli-Euler.                           |
|    - Khắc phục nguy cơ đánh giá thiếu chuyển vị khi thiết kế cọc cứng hoặc cọc có L/d < 10.           |
+--------------------------------------------------------------------------------------------------------+
| 3. KHUẾCH ĐẠI ĐỘNG HỌC VÀ HIỆN TƯỢNG CỘNG HƯỞNG TẦN SỐ                                                |
|    - Đỉnh cộng hưởng biên độ chuyển vị xuất hiện rõ rệt tại tần số f = 4.9 Hz và f = 5.2 Hz.           |
|    - Biên độ chuyển vị ngang động lực học tăng gấp 2.2 - 3.1 lần so với trạng thái tĩnh tương đương.  |
+--------------------------------------------------------------------------------------------------------+
| 4. BẢY SÓNG CẮT VÀ ĐIỀU KIỆN HÌNH THÀNH SÓNG LOVE                                                      |
|    - Khi có sự tương phản vận tốc sóng cắt (G1/rho1 < G2/rho2), sóng Love hình thành mạnh ở lớp mặt.   |
|    - Chuyển vị ngang tại bề mặt tăng vọt từ 1.8 đến 2.6 lần so với nền đồng nhất không phân lớp.      |
+--------------------------------------------------------------------------------------------------------+
| 5. CƠ CHẾ PHẢN ỨNG NỘI LỰC THEO PHA THỜI GIAN ĐỘNG ĐẤT                                                 |
|    - Đỉnh mômen và lực cắt cực đại xuất hiện tức thời tại các mốc t = 0.56s, 3.12s, 8.24s, 18.48s.     |
|    - Điểm uốn nguy hiểm có xu hướng dịch chuyển sâu xuống thân cọc (z = 4d - 6d) dưới sóng địa chấn.  |
+--------------------------------------------------------------------------------------------------------+
  1. Quy luật tắt dần của chuyển vị và phân bố nội lực theo chiều sâu: Dưới tác dụng của tải trọng tĩnh ngang, $80% - 90%$ phản lực đất nền tập trung trong phạm vi độ sâu hữu hiệu $L_e \approx (3 \div 5)d$ (với $d$ là đường kính hoặc bề rộng cọc). Mômen uốn cực đại $M_{\max}$ luôn xuất hiện ở độ sâu $z \approx (1.5 \div 2.5)d$ đối với cọc ngàm tự do đầu cọc.
  2. Ảnh hưởng định lượng của biến dạng trượt Timoshenko: Biến dạng trượt $\gamma_c$ đóng góp làm tăng chuyển vị ngang đầu cọc từ $8.2%$ đến $14.1%$ tùy thuộc vào tỷ lệ độ mảnh $L/d$ và tỷ số độ cứng tương đối giữa cọc và đất $E_c/E_d$. Đối với cọc ngắn cứng ($L/d \le 8$), việc bỏ qua biến dạng trượt như các nghiên cứu trước đây dẫn đến sai số nghiêm trọng trong tính toán ổn định chuyển vị.
  3. Hiện tượng cộng hưởng tần số của hệ cọc - nền: Khảo sát dao động động lực học trong dải tần $0.5 - 60 \text{ Hz}$ chỉ ra rằng biên độ dao động ngang của cọc đạt giá trị đỉnh tại tần số dao động riêng thứ nhất $f_1 \approx 4.9 \text{ Hz}$ và tần số thứ hai $f_2 \approx 5.2 \text{ Hz}$. Tại các tần số này, chuyển vị động của cọc tăng vọt từ $2.2$ đến $3.1$ lần so với chuyển vị tĩnh tương đương.
  4. Cơ chế bẫy sóng và khuếch đại địa chấn của sóng Love: Khi tồn tại lớp đất mặt mềm phủ trên nền đá hoặc tầng đất cứng hơn ($G_1/\rho_1 < G_2/\rho_2$), hiện tượng giao thoa sóng cắt tạo nên sóng bề mặt Love. Chuyển vị ngang tại lớp đất bề mặt bị khuếch đại từ $1.8$ đến $2.6$ lần so với mô hình nền đồng nhất, làm gia tăng đột biến lực cắt quán tính tác dụng lên cổ cọc tại vị trí tiếp giáp đài móng.
  5. Đặc trưng đáp ứng động đất El Centro theo thời gian: Dưới kích động địa chấn El Centro, biểu đồ chuyển vị, lực cắt và mômen uốn dọc thân cọc ghi nhận các đỉnh phản ứng nguy hiểm nhất tại các thời điểm $t = 0.56 \text{ s}, 3.12 \text{ s}, 8.24 \text{ s}$ và $18.48 \text{ s}$. Đặc biệt, dưới tác động động đất, điểm xuất hiện mômen uốn cực đại có xu hướng dịch chuyển sâu hơn xuống thân cọc ($z \approx 4d \div 6d$) so với bài toán tải trọng tĩnh đặt tại đỉnh.

Implications đa chiều

  • Về mặt lý thuyết: Thiết lập thành công mô hình giải tích - số kết hợp biến phân Gauss cho bài toán tương tác đất - kết cấu (SSI), bổ sung một phương pháp luận chuẩn mực bên cạnh phương pháp sai phân hữu hạn (FDM) và phương pháp phần tử biên (BEM).
  • Về mặt phương pháp luận: Xây dựng quy trình xử lý bài toán vô hạn thông qua hệ so sánh không gian tương đương, giúp loại bỏ sự cần thiết của các phần tử vô hạn (infinite boundary elements) phức tạp trong mô phỏng cơ học địa kỹ thuật 3D.
  • Về ứng dụng thực tiễn: Cung cấp bộ công cụ tính toán chính xác để tối ưu hóa chiều dài cọc, đường kính cọc và bố trí cốt thép chịu uốn - cắt hợp lý dọc thân cọc. Kỹ sư có thể giảm bớt mật độ cốt thép ở phần chân cọc ($z > 8d$) và tập trung gia cường vùng đầu cọc ($z \le 5d$), giúp tiết kiệm $12% - 18%$ chi phí vật liệu cọc bê tông cốt thép trong các công trình kháng chấn.
  • Về chính sách và tiêu chuẩn thiết kế: Đóng góp cơ sở khoa học trực tiếp cho việc rà soát và cập nhật các tiêu chuẩn thiết kế công trình chịu động đất tại Việt Nam (như TCVN 9386:2012 / Eurocode 8), đặc biệt là quy định về hệ số tương tác động lực học giữa móng cọc và đất nền loại D, E (đất yếu, đất trầm tích ven sông, ven biển).

Limitations và Future Research

Nghiên cứu thể hiện tính khách quan học thuật thông qua việc thừa nhận 4 giới hạn cụ thể:

  1. Giả thiết vật liệu đàn hồi tuyến tính: Đất nền được giả thiết là môi trường đàn hồi tuyến tính vi phân hoặc bán đàn hồi liên tục. Mô hình chưa xét đến các quy luật ứng xử phi tuyến cấp cao khi biến dạng lớn như phá hủy dẻo Mohr-Coulomb, Drucker-Prager hoặc hiện tượng hóa lỏng cát (liquefaction) khi áp lực nước lỗ rỗng dâng cao trong động đất.
  2. Quy mô cọc đơn độc lập: Luận án tập trung chuyên sâu vào bài toán cọc đơn (single pile). Hiệu ứng nhóm cọc (pile group interaction), hiệu ứng che chắn sóng (shadowing effect) và sự tương tác phức tạp giữa các cọc thông qua đài cọc ngầm chưa được tích hợp trong cùng một phiếm hàm biến phân.
  3. Mặt tiếp xúc cọc - đất lý tưởng: Giả định tiếp xúc dính bám hoàn toàn tại biên giữa cọc và đất ($u_c = u_d$), chưa xét đến hiện tượng trượt phi tuyến (slipping) hoặc tách rời cục bộ (gapping/separation) thường xảy ra ở vùng nông sát mặt đất khi cọc chịu lực ngang lắc qua lại với biên độ lớn.
  4. Dữ liệu kích động địa chấn: Phân tích động đất mới thực hiện trên thành phần gia tốc nằm ngang theo một phương chính của trận động đất El Centro, chưa khảo sát tác động đồng thời của sóng động đất đa chiều không gian 3D (kết hợp đồng thời thành phần gia tốc đứng $a_z$ và hai phương ngang $a_x, a_y$).

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Hướng 1: Tích hợp mô hình đàn - dẻo phi tuyến (Elasto-plasticity) cho đất nền và phần tử tiếp xúc trượt/hở (Interface elements with Coulomb friction & gap opening) vào phiếm hàm cực trị Gauss.
  • Hướng 2: Mở rộng hệ chương trình tính toán cho đài cọc và nhóm cọc không gian 3D chịu tải trọng lệch tâm và xoắn.
  • Hướng 3: Nghiên cứu bài toán tương tác cọc - đất bão hòa nước nhiều pha theo lý thuyết môi trường xốp Biot dưới tác động của sóng động đất gây hóa lỏng.
  • Hướng 4: Tiến hành thí nghiệm mô hình bàn rung (shake table test) hoặc mô hình ly tâm (centrifuge test) quy mô lớn để hiệu chỉnh chính xác các tham số cản bức xạ và cản vật liệu trong điều kiện địa chất đặc thù tại Việt Nam.

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

+--------------------------------------------------------------------------------------------------------+
|                                    BẢN ĐỒ TÁC ĐỘNG VÀ ẢNH HƯỞNG TOÀN DIỆN                              |
+--------------------------------------------------------------------------------------------------------+
| HỌC THUẬT (ACADEMIC IMPACT)                                                                            |
| - Tiên phong ứng dụng Nguyên lý cực trị Gauss trong SSI tại Việt Nam.                                  |
| - Cung cấp bộ công cụ 9 chương trình tính toán mã nguồn mở nội bộ cho đào tạo Sau đại học.            |
+--------------------------------------------------------------------------------------------------------+
                                                    |
                                                    v
| CÔNG NGHIỆP XÂY DỰNG (INDUSTRY TRANSFORMATION)                                                         |
| - Tối ưu hóa thiết kế móng cọc công trình cầu cảng, nhà cao tầng, nhiệt điện/thủy điện.               |
| - Tiết kiệm 12% - 18% cốt thép thân cọc nhờ xác định chính xác biểu đồ bao mômen uốn và lực cắt.       |
+--------------------------------------------------------------------------------------------------------+
                                                    |
                                                    v
| AN TOÀN XÃ HỘI & CHÍNH SÁCH (POLICY & SOCIETAL BENEFIT)                                                |
| - Đảm bảo an toàn kháng chấn cho công trình trọng điểm quốc gia tại các vùng đứt gãy hoạt động.       |
| - Cơ sở khoa học bổ sung vào hệ thống Tiêu chuẩn Quốc gia TCVN và quy chuẩn xây dựng an toàn.          |
+--------------------------------------------------------------------------------------------------------+
  • Tác động học thuật: Luận án là một trong những công trình tiên phong tại Việt Nam ứng dụng thành công biến phân Nguyên lý cực trị Gauss vào cơ học kết cấu móng công trình ngầm. Các thuật toán và chương trình tính toán trong luận án tạo nền tảng vững chắc cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Xây dựng Dân dụng & Công nghiệp, Công trình ngầm và Cơ học Đất - Địa kỹ thuật công trình.
  • Chuyển đổi công nghiệp xây dựng: Cung cấp giải pháp tính toán kiểm tra độc lập có độ tin cậy cao cho các tổng công ty tư vấn thiết kế lớn khi thực hiện các dự án siêu cao tầng, cầu dây văng nhịp lớn, bến cảng nước sâu hoặc các công trình năng lượng trọng điểm tại những vùng có nguy cơ địa chấn cấp 7 - cấp 9.
  • Lợi ích kinh tế - xã hội: Việc mô phỏng chính xác ứng xử cọc - đất giúp ngăn ngừa sự cố phá hủy giòn kết cấu móng ngầm dưới động đất — loại sự cố cực kỳ tốn kém và gần như không thể khắc phục sau thảm họa. Đồng thời, việc tối ưu hóa vùng bố trí cốt thép dọc thân cọc mang lại hiệu quả kinh tế hàng chục tỷ đồng cho mỗi dự án quy mô lớn.

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

  • Nghiên cứu sinh và Giảng viên đại học: Tiếp cận một khung lý thuyết biến phân toán học chặt chẽ, các ma trận độ cứng phần tử dầm Timoshenko mở rộng và phương pháp khử biên vô hạn bằng hệ so sánh để phát triển các đề tài nghiên cứu chuyên sâu về cơ học tính toán.
  • Kỹ sư Kết cấu & Kỹ sư Địa kỹ thuật (R&D / Structural Engineers): Sở hữu quy trình tính toán rõ ràng, hiểu sâu sắc cơ chế phân phối tải trọng ngang và tác động khuếch đại sóng cắt địa chấn, từ đó thiết kế móng cọc an toàn, tránh việc áp dụng máy móc các công thức kinh nghiệm đơn giản hóa quá mức.
  • Chủ đầu tư và Ban Quản lý Dự án Trọng điểm: Nắm bắt cơ sở khoa học để thẩm tra, đánh giá mức độ an toàn và tối ưu hóa chi phí đầu tư xây dựng phần ngầm cho các công trình hạ tầng kỹ thuật quốc gia.
  • Cơ quan Quản lý Nhà nước (Bộ Xây dựng, Viện Vật lý Địa cầu): Sử dụng các kết quả nghiên cứu định lượng làm luận cứ khoa học để hoàn thiện bản đồ phân vùng rủi ro động đất gắn liền với điều kiện địa chất công trình và tiêu chuẩn thiết kế kháng chấn.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập phiếm hàm lượng cưỡng bức toàn phần $Z = Z_d + Z_c + Z_{AB}$ dựa trên Nguyên lý cực trị Gauss cho bài toán tương tác cọc - đất 3D liên tục có xét biến dạng trượt Timoshenko. Luận án đã mở rộng lý thuyết biến phân Gauss (vốn chủ yếu dùng trong động lực học hệ chất điểm và vật rắn) sang cơ học môi trường liên tục vô hạn, đồng thời mở rộng lý thuyết dầm Timoshenko cổ điển bằng cách tích hợp liên kết chuyển vị tiếp xúc thông qua phương pháp nhân tử Lagrange $\lambda(z)$.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu tiền nhiệm thể hiện ở điểm nào?

So với phương pháp tích phân nửa giải tích bán không gian của Poulos (1971) và phương pháp hàm trở kháng động lực học của Novak (1974), luận án đổi mới căn bản ở kỹ thuật chuyển đổi bài toán bán không gian về hệ so sánh không gian vô hạn đàn hồi. Giải pháp này cho phép tận dụng triệt để nghiệm cơ bản của không gian vô hạn để khử toàn bộ các điều kiện biên phức tạp ở vô cực mà không cần thiết lập lưới phần tử vô hạn (infinite elements) hay phần tử biên (BEM), giảm thiểu tối đa chi phí tính toán số.

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ và có ý nghĩa cảnh báo kỹ thuật cao nhất là: Vị trí xuất hiện mômen uốn cực đại ($M_{\max}$) của cọc dưới tác động của sóng động đất El Centro dịch chuyển sâu xuống lòng đất ($z \approx 4d \div 6d$) so với bài toán tải trọng tĩnh thông thường ($z \approx 1.5d \div 2.5d$). Điều này giải thích nguyên nhân vì sao nhiều công trình móng cọc bị gãy thân cọc ở độ sâu lớn trong các trận động đất lịch sử dù đã được thiết kế chịu lực cắt và mômen tĩnh rất lớn ở vùng đỉnh cọc sát đài móng.

+--------------------------------------------------------------------------------------------------------+
|                  SO SÁNH VỊ TRÍ MÔMEN UỐN CỰC ĐẠI: TẢI TRỌNG TĨNH VS TẢI TRỌNG ĐỘNG ĐẤT                 |
+--------------------------------------------------------------------------------------------------------+
| Độ sâu z                                                                                               |
|   0 +------+ <--- Đỉnh cọc (Mặt đất)                                                                  |
|     |      |                                                                                           |
|  2d | **** | <=== M_max dưới TẢI TĨNH NGANG (Tập trung vùng nông z = 1.5d - 2.5d)                      |
|     |      |                                                                                           |
|  4d |      |                                                                                           |
|     | #### | <=== M_max dưới ĐỘNG ĐẤT EL CENTRO (Dịch chuyển sâu xuống z = 4d - 6d do quán tính)       |
|  6d | #### |                                                                                           |
|     |      |                                                                                           |
|  8d |      |                                                                                           |
|     |      |                                                                                           |
|  L  +------+ <--- Mũi cọc / Chân cọc                                                                   |
+--------------------------------------------------------------------------------------------------------+

4. Luận án có cung cấp đầy đủ giao thức để tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp hệ thống giao thức tái lập hoàn chỉnh bao gồm: (1) Ma trận hàm nội suy và phương trình phần tử khối chữ nhật 20 nút; (2) Ma trận độ cứng phần tử dầm Timoshenko có xét biến dạng trượt; (3) Thuật toán biến đổi Fourier rời rạc DFT/FFT và sơ đồ khối tích phân Duhamel; (4) Mã nguồn chi tiết của 9 chương trình tính toán số từ Phụ lục 6 đến Phụ lục 14 (Mstatic1, Kstatic1, MstaticP1, KstaticP1, KstaticPLs, KdynaS, KdynaL, KdynaP, KdynaPE).

5. Định hướng nghiên cứu 10 năm tiếp theo từ nền tảng luận án được vạch ra như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào 3 mũi nhọn: (1) Tích hợp phi tuyến vật liệu đất nền đa pha (đất - nước - khí) và cơ chế hóa lỏng động đất; (2) Phát triển phần mềm mô phỏng không gian 3D tương tác toàn diện hệ Kết cấu bên trên - Đài móng - Nhóm cọc - Đất nền phân lớp dị hướng; (3) Số hóa dữ liệu gia tốc đồ các trận động đất đặc trưng tại Việt Nam để xây dựng phần mềm chuyên dụng hỗ trợ thiết kế kháng chấn tự động theo tiêu chuẩn quốc gia.

Kết luận

Luận án tiến sĩ của tác giả Ngô Quốc Trinh là một công trình khoa học công phu, mẫu mực và có tính đột phá cao trong lĩnh vực Kỹ thuật Xây dựng Dân dụng & Công nghiệp. Công trình ghi nhận 6 đóng góp mang tính nền tảng:

  1. Ứng dụng thành công Nguyên lý cực trị Gauss để giải quyết bài toán tương tác tĩnh và động giữa cọc và đất nền 3D liên tục;
  2. Tích hợp chuẩn xác biến dạng trượt ngang dầm Timoshenko vào phiếm hàm năng lượng cưỡng bức tổng thể thông qua kỹ thuật nhân tử Lagrange;
  3. Đề xuất phương pháp hệ so sánh không gian đàn hồi vô hạn, giải quyết triệt để bài toán biên vô cực trong cơ học địa kỹ thuật số;
  4. Làm sáng tỏ quy luật lan truyền sóng Love và hiện tượng khuếch đại chuyển vị trong nền đất phân lớp có độ cứng tương phản;
  5. Định lượng chính xác biên độ cộng hưởng tần số và cơ chế phân bố nội lực nguy hiểm dọc thân cọc dưới kích động động đất El Centro;
  6. Xây dựng và kiểm chứng thành công bộ 9 chương trình tính toán số chuyên dụng, phục vụ trực tiếp cho công tác nghiên cứu và thực hành thiết kế kỹ thuật công trình.

Công trình không chỉ mở ra hướng đi mới cho lý thuyết cơ học môi trường liên tục ứng dụng tại Việt Nam mà còn cung cấp cơ sở khoa học vững chắc phục vụ sự nghiệp công nghiệp hóa, hiện đại hóa và bảo đảm an toàn công trình quốc gia trước các hiểm họa thiên tai địa chấn.