Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Động đất là hiện tượng thiên tai nguy hiểm, gây ra những chuyển động bất lợi tại nền móng công trình và tiềm ẩn nguy cơ phá hoại nghiêm trọng kết cấu xây dựng. Trên lãnh thổ Việt Nam, mặc dù không nằm trong vành đai lửa địa chấn lớn của thế giới, nhưng thực tế tồn tại nhiều hệ thống đứt gãy hoạt động phức tạp như đới đứt gãy Lai Châu - Điện Biên, Sông Mã, Sơn La, Sông Hồng và Sông Cả. Dữ liệu ghi nhận từ đầu thế kỷ 20 đến nay cho thấy khu vực phía Bắc đã từng chịu 2 trận động đất cấp 8–9 theo thang MSK-64 (tương đương 6,7–6,8 độ Richter), tiêu biểu là trận động đất Điện Biên năm 1935 (6,7 độ Richter trên đới Sông Mã) và trận động đất Tuần Giáo năm 1983 (6,8 độ Richter trên đứt gãy Sơn La gây sạt lở diện rộng trong bán kính 35 km). Ngoài ra còn có hàng chục trận động đất cấp 7 (5,1–5,5 độ Richter) cùng hàng trăm chấn động yếu hơn, kết hợp với hiện tượng động đất kích thích ghi nhận tại khu vực thủy điện Sông Tranh 2 (Bắc Trà My, Quảng Nam).

Trong bối cảnh xây dựng nhiều công trình quy mô lớn như nhà máy điện hạt nhân, công trình ngoài khơi, đập thủy điện, cầu cảng lớn và nhà cao tầng, yêu cầu thiết kế kháng chấn cho móng cọc trở nên đặc biệt quan trọng. Hệ thống tiêu chuẩn thiết kế kháng chấn hiện hành tại Việt Nam thời điểm nghiên cứu (như tiêu chuẩn ngành 22 TCN 272-05 biên dịch từ AASHTO LRFD 1998, TCXDVN 375:2006 biên dịch từ Eurocode 8, 22 TCN 221-95 biên dịch từ XNiP II-7-81*, hay TCXD 205-1998) chủ yếu đưa ra các chỉ dẫn chung về giảm sức chịu tải hoặc xét lực quán tính tác dụng lên đầu cọc, nhưng còn thiếu các hướng dẫn cụ thể về phân tích tương tác động học giữa cọc và đất nền.

Khoảng trống nghiên cứu then chốt xuất phát từ việc tính toán trạng thái ứng suất - biến dạng của móng cọc chịu tải trọng ngang và động đất hiện nay chủ yếu sử dụng các mô hình đơn giản hóa (mô hình nền Winkler với lò xo tuyến tính/phi tuyến hoặc mô hình nửa không gian đàn hồi liên tục). Việc xác định các thông số thực nghiệm như hệ số phản lực nền $k_h$, đường cong $p-y$, hay hệ số cản nhớt bức xạ gặp nhiều khó khăn, thiếu độ chính xác và phụ thuộc chủ quan vào từng tác giả. Phần lớn phương pháp mới chỉ xét một chiều tác động của đất lên cọc mà chưa xem xét đồng thời ảnh hưởng ngược lại của cọc lên nền đất, đồng thời gặp trở ngại lớn khi xử lý điều kiện biên vô hạn trong bài toán truyền sóng địa chấn.

Mục tiêu nghiên cứu

Mục tiêu chính của luận án là xây dựng phương pháp lý thuyết nghiên cứu bài toán tương tác giữa cọc và nền đất, đồng thời lập hệ thống phần mềm tính toán nhằm xác định chính xác trạng thái ứng suất - biến dạng của cọc chịu tải trọng ngang và tải trọng động đất.

Nội dung nghiên cứu cụ thể được phân định thành 4 phần:

  1. Nghiên cứu trạng thái ứng suất - biến dạng của khối đất chịu tải trọng tĩnh nằm ngang.
  2. Nghiên cứu bài toán tương tác tĩnh học giữa cọc với nền đất khi chịu tải trọng tĩnh nằm ngang.
  3. Nghiên cứu bài toán tương tác động lực học giữa cọc với nền đất khi chịu tải trọng động nằm ngang và chịu tải trọng động đất trong miền tần số và miền thời gian.
  4. Xây dựng phần mềm tính toán mô phỏng cho các trường hợp nghiên cứu nêu trên.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Cọc đơn thẳng đứng nằm trong bán không gian vô hạn đàn hồi, chịu tác dụng của tải trọng tĩnh nằm ngang, tải trọng động nằm ngang bất kỳ và tải trọng sóng động đất.
  • Phạm vi nghiên cứu (giới hạn đề tài): Luận án tập trung trong môi trường đàn hồi, không tính toán trong các mô hình nền đất khác (đàn dẻo, đàn nhớt); không xét hiện tượng hóa lỏng của nền đất khi xảy ra động đất; không xét ảnh hưởng của áp lực nước lỗ rỗng trong nền đất bão hòa; và không nghiên cứu bài toán trạng thái giới hạn phá hoại của kết cấu cọc.

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

Các hướng nghiên cứu trước đây

Luận án hệ thống hóa các phương pháp nghiên cứu tương tác cọc - nền đất chịu tải trọng ngang theo hai nhóm lý thuyết chính:

1. Nhóm các phương pháp dựa trên mô hình nền Winkler

Nhóm phương pháp này xem cọc như dầm chịu uốn trên nền đàn hồi cục bộ, phản lực của đất $p$ tại một điểm tỷ lệ với chuyển vị ngang $y$ thông qua phương trình vi phân $E_c J \frac{d^4 y}{dz^4} + p = 0$:

  • Mô hình lò xo tuyến tính: Chang (1937) đưa ra phương pháp tính với hệ số nền hằng số theo độ sâu và đề xuất khái niệm chiều dài tới hạn của cọc dựa trên việc phân tích số liệu thực nghiệm của Feagin, xác định hệ số nền $k_h$ ứng với giá trị $E_s/3$. Các nghiên cứu tiếp theo của Reese và Matlock (1956), Zavriev (1962), Trần Bình (1968) tiếp tục phát triển hướng tính toán này.
  • Phương pháp đường cong $p-y$ (lò xo phi tuyến): Matlock (1970), Reese (1974), Bhushan (1979) mở rộng bài toán ra ngoài miền đàn hồi bằng cách thiết lập quan hệ phi tuyến $p = k_h f(y)$. Trong đất cát, phương trình đường cong chung có dạng $p = A p_u \tanh\left(\frac{k_h H y}{A p_u}\right)$. Trong đất sét, Matlock thiết lập quan hệ $p = 0{,}5 p_u (y / y_{50})^{1/3}$ với $y_{50} = 2{,}5 B \varepsilon_{50}$. Tiêu chuẩn của Viện Dầu khí Hoa Kỳ (API, 1993) đã ứng dụng đường cong $p-y$ cho móng cọc giàn khoan biển. Tuy nhiên, Murchison và O’Neill (1984) qua khảo sát thực nghiệm đã chỉ ra đường cong API có sai số lớn do bỏ sót cơ chế truyền lực cắt giữa các lò xo.
  • Mô hình động lực học và trở kháng cơ học:
    • Mô hình dầm trên lò xo phi tuyến và hộp nhớt trong các phần mềm SPASM (Matlock, 1978), NONSPS (Kagawa và Kraft, 1980), PAR (PMB Engineering, 1988), và PILE-PY của Thavaraj (2001) giải phương trình chuyển động $E_c J \frac{\partial^4 v}{\partial z^4} + \rho A (\frac{\partial^2 v}{\partial t^2} + \frac{\partial^2 v_g}{\partial t^2}) + c (\frac{\partial v}{\partial t} - \frac{\partial v_{ff}}{\partial t}) + k_h (v - v_{ff}) = 0$.
    • Quy trình 3 bước của Gazetas, Fan và Kaynia (1993): tách riêng phản ứng động học (kinematic interaction) của móng không khối lượng, xác định độ cứng phức/trở kháng cơ học ($K_x, K_{ry}, K_{x-ry}$ cùng hệ số cản $c_x$), sau đó phân tích tương tác quán tính (inertial interaction) của kết cấu công trình đặt trên hệ lò xo độ cứng phức.
    • Mô hình miền gần - miền xa của El Naggar và Novak (1995, 1996): sử dụng lò xo phi tuyến cho miền gần (Novak và Sheta) và độ cứng phức của bài toán biến dạng phẳng cho miền xa để mô phỏng sự bức xạ sóng.
+-----------------------------------------------------------------------------+
|               CÁC PHƯƠNG PHÁP TƯƠNG TÁC CỌC - ĐẤT NỀN                       |
+--------------------------------------+--------------------------------------+
|        Mô hình nền Winkler           |     Mô hình đàn hồi liên tục         |
|  (Lò xo tuyến tính & phi tuyến p-y)  | (Không gian / Bán không gian đàn hồi)|
+--------------------------------------+--------------------------------------+
| - Chang (1937), Zavriev (1962)       | - Lê Đức Thắng (1966, 1968)          |
| - Reese & Matlock (1956)             | - Poulos (1971a, 1971b)              |
| - Matlock (1970), Reese (1974)       | - Verruijt & Kooijman (1989)         |
| - Gazetas, Fan & Kaynia (1993)       | - Trần Hữu Hà (2006)                 |
| - El Naggar & Novak (1996)           | - Hướng phát triển của luận án       |
| - Thavaraj (2001 - PILE-PY)          |   (PPNLCT Gauss + Timoshenko 3D)     |
+--------------------------------------+--------------------------------------+

2. Nhóm các phương pháp dựa trên mô hình đàn hồi liên tục

  • Lê Đức Thắng (1966, 1968) là một trong những tác giả đầu tiên ứng dụng lời giải Mindlin cho lực tập trung nằm ngang tĩnh trong bán không gian vô hạn đàn hồi, thiết lập phương trình tương tác bằng cách chập chuyển vị của nền đất với chuyển vị của cọc.
  • Poulos (1971a,b) chia cọc thành $(n+1)$ phần tử, sử dụng lời giải Mindlin để tích phân chuyển vị đất và cho chập chuyển vị với cọc tại tâm các phân tố, bỏ qua ứng suất tiếp ngang trên mặt bên cọc.
  • Verruijt và Kooijman (1989) đề xuất mô hình cọc trong các lớp đất đàn hồi giải bằng thuật toán lặp 4 bước kết hợp giữa phương pháp phần tử hữu hạn cho các lớp đất và phương pháp sai phân hữu hạn cho thân cọc.
  • Trần Hữu Hà (2006) ứng dụng Phương pháp Nguyên lý cực trị Gauss (PPNLCT Gauss) sử dụng lời giải Mindlin cho bài toán tĩnh và lời giải không gian vô hạn đàn hồi cho bài toán động làm hệ so sánh. Tuy nhiên, công trình mới chỉ dừng lại ở cọc chịu lực dọc (cọc không chịu uốn) và chưa xét tải trọng ngang hay sóng động đất.

Khoảng trống luận án giải quyết

Luận án mở rộng việc ứng dụng PPNLCT Gauss kết hợp lý thuyết dầm Timoshenko (xét biến dạng trượt ngang của cọc) và phương pháp phần tử hữu hạn 3 chiều 20 nút cho môi trường đất. Cách tiếp cận này giúp khắc phục triệt để việc phải đặt các liên kết phụ giả định (lò xo, hộp nhớt) trên bề mặt tiếp xúc cọc - đất cũng như trên biên nhân tạo, tự động thỏa mãn điều kiện biên ở vô cùng trong cả bài toán tĩnh và bài toán động lực học truyền sóng địa chấn.

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

Cơ sở lý thuyết

  1. Phương pháp Nguyên lý cực trị Gauss (PPNLCT Gauss): Sử dụng phiếm hàm lượng cưỡng bức Gauss $Z$ đối với môi trường liên tục và hệ kết cấu. Bằng kỹ thuật sử dụng hệ so sánh (lời giải Mindlin cho bán không gian vô hạn đàn hồi trong bài toán tĩnh, hoặc lời giải Kelvin / lời giải xung đơn vị cho không gian vô hạn đàn hồi trong bài toán động), trạng thái ứng suất - biến dạng thực của hệ cần tính được thiết lập thông qua sai lệch lượng cưỡng bức cực tiểu.
  2. Lý thuyết dầm Timoshenko: Mô hình hóa cọc chịu uốn đồng thời có xét đến biến dạng trượt ngang, phản ánh chính xác trạng thái làm việc của cọc ngắn, cọc có độ cứng lớn hoặc chịu tải trọng động tần số cao.
  3. Lý thuyết động lực học và truyền sóng đàn hồi: Khảo sát sự lan truyền của sóng vật thể (sóng dọc P, sóng cắt S gồm SV và SH) và sóng bề mặt (sóng Rayleigh, sóng Love). Phân tích hiệu ứng khuếch đại dao động khi sóng truyền từ nền đá cứng qua các lớp trầm tích mềm lên mặt đất.
  4. Hệ số giảm chấn vật liệu (hysteretic damping / ma sát khô): Thay vì sử dụng hệ số cản nhớt phụ thuộc tần số thông thường, luận án sử dụng hệ số ma sát khô $\zeta_h$ (không phụ thuộc tần số), cho phép mô tả chính xác sự tiêu tán năng lượng trong môi trường đất và tạo tiền đề mở rộng xét biến dạng dẻo.
+-----------------------------------------------------------------------------+
|                 SƠ ĐỒ KHUNG PHÂN TÍCH LÝ THUYẾT VÀ SỐ                       |
+-----------------------------------------------------------------------------+
|  [Hệ so sánh Giải tích]             [Hệ cần tính Thực tế]                   |
|  - Tĩnh: Bán không gian Mindlin     - Khối đất 3D (Phần tử 20 nút)          |
|  - Động: Không gian vô hạn Kelvin   - Cọc Timoshenko (Uốn + Trượt ngang)    |
|             \                                /                              |
|              \                              /                               |
|               v                            v                                |
|        +------------------------------------------+                         |
|        |   PHIẾM HÀM LƯỢNG CƯỠNG BỨC GAUSS (Z)    |                         |
|        |   delta(Z) = 0 -> Hệ phương trình FEM    |                         |
|        +------------------------------------------+                         |
|                             |                                               |
|                             v                                               |
|        +------------------------------------------+                         |
|        |   PHÂN TÍCH ĐỘNG LỰC HỌC SỐ              |                         |
|        |   - Gia tốc đồ El Centro 1940            |                         |
|        |   - Tích phân Duhamel / DFT / FFT / IFFT |                         |
|        +------------------------------------------+                         |
|                             |                                               |
|                             v                                               |
|        +------------------------------------------+                         |
|        |   KẾT QUẢ: Ứng suất, Chuyển vị,          |                         |
|        |            Mô men uốn, Lực cắt cọc       |                         |
|        +------------------------------------------+                         |
+-----------------------------------------------------------------------------+

Phương pháp nghiên cứu và công cụ phân tích

  • Phương pháp số: Phương pháp phần tử hữu hạn (FEM). Môi trường đất được rời rạc hóa bằng phần tử khối 3 chiều 20 nút bậc hai; phần tử cọc sử dụng 2 nút cho trường chuyển vị ngang và 3 nút cho phân bố lực cắt.
  • Kỹ thuật xử lý tín hiệu động lực học:
    • Biến đổi Fourier rời rạc (DFT) và Biến đổi Fourier nhanh (FFT) để chuyển dữ liệu sóng động đất từ miền thời gian sang miền tần số.
    • Tích phân chập Duhamel để xác định phản ứng dao động trong miền tần số.
    • Biến đổi Fourier nhanh ngược (IFFT) để đưa kết quả phản ứng nội lực và chuyển vị về miền thời gian thực.
  • Nguồn dữ liệu thực nghiệm kiểm chứng: Gia tốc đồ thực tế theo hướng Bắc - Nam của trận động đất El Centro (1940); số liệu địa chất và kết quả thí nghiệm tải trọng cọc từ các nghiên cứu kinh điển của Kim, O'Neill, Matlock (1984), Poulos (1971), và Zavriev (1962).

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

Chương 1: Tổng quan các phương pháp nghiên cứu tương tác giữa cọc với nền đất chịu tải trọng ngang

Chương 1 tổng thuật chi tiết các khái niệm cơ bản về động chấn công trình: nguồn gốc động đất do đứt gãy kiến tạo (thuyết bật đàn hồi của H.F. Reid) và kiến tạo mảng (thuyết lục địa trôi của H. Hess); phân loại sóng động đất (sóng khối P, S và sóng mặt Rayleigh, Love); các thang đo cường độ địa chấn (thang định tính MM, MSK-64 và thang năng lượng Richter $M_s$ với quan hệ $\log E = 11{,}8 + 1{,}5 M_s$); các thông số chuyển động nền như gia tốc cực đại (PHA), vận tốc cực đại (PHV), chuyển vị cực đại (PHD).

Nội dung chương đi sâu phân tích ưu và nhược điểm của các phương pháp tính toán móng cọc chịu tải trọng ngang hiện nay (Winkler, đường cong $p-y$, mô hình đàn hồi liên tục Mindlin, phương pháp độ cứng phức của Gazetas, mô hình tương tác miền gần - miền xa của Novak). Trên cơ sở chỉ ra hạn chế của các phương pháp truyền thống trong việc xử lý điều kiện biên vô hạn và xét tương tác hai chiều cọc - đất, chương 1 xác lập cơ sở khoa học cho hướng tiếp cận của đề tài.

Chương 2: Nghiên cứu trạng thái ứng suất biến dạng của nền đất chịu tải trọng tĩnh nằm ngang

Chương 2 trình bày việc thiết lập các phương trình cơ bản của môi trường đàn hồi và phương trình truyền sóng thông qua PPNLCT Gauss. Tác giả xây dựng bài toán tương tác giữa khối đất hữu hạn cần tính với nửa không gian vô hạn đàn hồi hoặc không gian vô hạn đàn hồi (hệ so sánh).

Bằng phương pháp phần tử hữu hạn với phần tử khối 20 nút, chương trình giải số cho phép xác định chính xác trường ứng suất và chuyển vị của khối đất khi chịu lực ngang đặt tại mặt thoáng hoặc ở độ sâu bất kỳ dưới đáy khối đất. Kết quả kiểm tra số học chứng minh tính hội tụ và độ chính xác khi đối chiếu giữa hệ so sánh không gian vô hạn (lời giải Kelvin) và nửa không gian vô hạn (lời giải Mindlin), đồng thời đánh giá định lượng ảnh hưởng của trọng lượng bản thân đất nền đến chuyển vị ngang và chuyển vị đứng.

Chương 3: Nghiên cứu bài toán tương tác giữa cọc với nền đất chịu tải trọng tĩnh nằm ngang

Chương 3 xây dựng bài toán dầm Timoshenko chịu uốn có xét biến dạng trượt ngang theo PPNLCT Gauss và giải bằng FEM. Thuật toán tích hợp mô hình dầm cọc vào khối đất 3 chiều trong hệ so sánh nửa không gian đàn hồi Mindlin hoặc không gian đàn hồi Kelvin.

Chương tiến hành các khảo sát số nhằm kiểm chứng độ tin cậy của phương pháp:

  • Kiểm tra tính độc lập của kết quả khi thay đổi mô đun đàn hồi $E_0$ của hệ so sánh ($E_0 = 10\text{ MPa}$ và $E_0 = 20\text{ MPa}$).
  • So sánh kết quả chuyển vị và mô men uốn với các phương pháp cổ điển: phương pháp Zavriev (1962), phương pháp Poulos (1971), và phương pháp đường cong $p-y$ theo số liệu thí nghiệm của Kim, O'Neill, Matlock dưới các cấp tải trọng ngang $200\text{ kN}$, $400\text{ kN}$, $600\text{ kN}$, $800\text{ kN}$.
  • Khảo sát các tham số ảnh hưởng: chiều dài cọc (cọc ngắn $L = 4\text{ m}$ và cọc dài $L = 16\text{ m}$), điều kiện mũi cọc (cọc treo trong nền đồng nhất và cọc tựa trên lớp đá cứng), sự thay đổi mô đun đàn hồi của đất ($E_d = 5, 10, 15, 20\text{ MPa}$).
  • Khảo sát ảnh hưởng ngược lại của sự hiện diện thân cọc đến trường chuyển vị của đất nền xung quanh.
Tham số khảo sát trong Chương 3 Các trường hợp / Giá trị thực hiện tính toán
Mô đun đàn hồi hệ so sánh $E_0$ $10\text{ MPa}$; $20\text{ MPa}$
Chiều dài cọc đơn $L$ $4\text{ m}$ (cọc ngắn); $16\text{ m}$ (cọc dài)
Cấp tải trọng ngang tĩnh $P$ $200\text{ kN}$; $400\text{ kN}$; $600\text{ kN}$; $800\text{ kN}$
Điều kiện biên chân cọc Chân tự do trong nền đồng nhất; Chân ngàm tựa trên đá cứng
Mô đun đàn hồi của nền đất $E_d$ $5\text{ MPa}$; $10\text{ MPa}$; $15\text{ MPa}$; $20\text{ MPa}$
Chương trình phần mềm đối chuẩn MstaticP1, KstaticP1, KstaticPLs

Chương 4: Nghiên cứu bài toán tương tác giữa cọc với nền đất chịu tải trọng động nằm ngang và tải trọng động đất

Chương 4 mở rộng thuật toán sang bài toán động lực học. Tác giả sử dụng lời giải xung đơn vị của không gian vô hạn đàn hồi trong miền thời gian, miền biến đổi Laplace và miền tần số Fourier, kết hợp hệ số giảm chấn vật liệu $\zeta_h$.

Nội dung khảo sát động lực học bao gồm:

  • Phân tích dao động của khối đất chịu tải trọng điều hòa với các dải tần số $0{,}5 - 30\text{ Hz}$ (bước tần $0{,}5\text{ Hz}$) và $1{,}0 - 60\text{ Hz}$ (bước tần $1{,}0\text{ Hz}$), xác định biên độ chuyển vị ngang và góc pha tại lớp mặt và lớp đáy.
  • Mô phỏng quá trình truyền sóng cắt (sóng Love) từ tầng đá cứng bên dưới qua lớp đất mặt mềm hơn ($G_1/\rho_1 < G_2/\rho_2$), làm rõ hiện tượng khuếch đại dao động bề mặt theo phương vuông góc với phương truyền sóng khi thay đổi tỷ số mô đun ($E_2 = 4, 6, 8\text{ MPa}$ trên lớp $E_1 = 20\text{ MPa}$).
  • Khảo sát dao động của cọc đơn chịu tải trọng động nằm ngang, xác định chuyển vị ngang theo tần số tại đầu cọc, giữa cọc và chân cọc (tại các tần số cộng hưởng $4{,}9\text{ Hz}$ và $5{,}2\text{ Hz}$).
  • Khảo sát cọc chịu tác động trực tiếp của sóng động đất El Centro (1940). Phân tích diễn biến chuyển vị ngang, lực cắt và mô men uốn dọc theo chiều dài cọc tại các mốc thời gian tiêu biểu ($t = 0{,}56\text{ s}$; $3{,}12\text{ s}$; $8{,}24\text{ s}$; và $18{,}48\text{ s}$).
+-----------------------------------------------------------------------------+
|             HỆ THỐNG CÁC PHẦN MỀM TÍNH TOÁN ĐƯỢC TÁC GIẢ THIẾT LẬP          |
+---------------+-------------------------------------------------------------+
| Tên chương    | Chức năng và phạm vi tính toán                              |
| trình         |                                                             |
+---------------+-------------------------------------------------------------+
| Mstatic1      | Tính khối đất chịu tải trọng tĩnh nằm ngang (hệ so sánh     |
|               | Mindlin)                                                    |
| Kstatic1      | Tính khối đất chịu tải trọng tĩnh nằm ngang (hệ so sánh     |
|               | Kelvin)                                                     |
| MstaticP1     | Tính cọc đơn chịu tải tĩnh ngang (hệ so sánh Mindlin)       |
| KstaticP1     | Tính cọc đơn chịu tải tĩnh ngang (hệ so sánh Kelvin)        |
| KstaticPLs    | Tính cọc nằm trong nền đàn hồi nhiều lớp chịu tải tĩnh      |
|               | ngang                                                       |
| KdynaS        | Phân tích dao động khối đất chịu tải trọng động nằm ngang   |
| KdynaL        | Mô phỏng truyền sóng cắt Love trong nền đất phân tầng       |
| KdynaP        | Tính cọc đơn chịu tải trọng động nằm ngang điều hòa         |
| KdynaPE       | Phân tích phản ứng cọc chịu tải trọng sóng động đất         |
|               | El Centro                                                   |
+---------------+-------------------------------------------------------------+

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

Luận án tự đúc kết 6 đóng góp mới về mặt khoa học lý thuyết và công cụ ứng dụng:

  1. Quy đổi nghiệm giải tích cơ bản thông qua phương pháp số: Thông qua lời giải số bằng phương pháp phần tử hữu hạn, đã chuyển đổi thành công lời giải Kelvin về lời giải Mindlin. Điều này cho phép nhận được nghiệm của bài toán bán không gian vô hạn đàn hồi trực tiếp từ nghiệm của không gian vô hạn đàn hồi với tải trọng đặt tại vị trí bất kỳ trong lòng môi trường.
  2. Xây dựng hoàn chỉnh bài toán tương tác cọc - đất 3 chiều: Xây dựng thành công mô hình toán học giải quyết trọn vẹn bài toán tương tác tĩnh học và động lực học giữa cọc với nền đất chịu tải trọng tĩnh và tải trọng động đặt tại vị trí bất kỳ. Mô hình kết hợp phần tử hữu hạn khối 3 chiều 20 nút cho đất và phần tử dầm Timoshenko (2 nút cho chuyển vị, 3 nút cho lực cắt) cho cọc.
  3. Tự động thỏa mãn điều kiện biên không cần liên kết phụ: Phương pháp PPNLCT Gauss sử dụng hệ so sánh tự động thỏa mãn điều kiện biên bức xạ sóng ra vô cùng, điều kiện biên trên bề mặt khối đất chứa cọc và điều kiện tiếp xúc liên tục giữa cọc và đất nền mà không cần đặt thêm các phần tử liên kết nhân tạo (như hệ lò xo Winkler hay hộp cản nhớt).
  4. Ứng dụng hệ số giảm chấn vật liệu: Trong tính toán động lực học, thay vì dùng hệ số cản nhớt phụ thuộc tần số thông thường, luận án sử dụng hệ số giảm chấn vật liệu (hysteretic damping / ma sát khô). Hệ số này độc lập với tần số kích thích và tạo tiền đề lý thuyết thuận lợi để mở rộng xét hiện tượng biến dạng dẻo của đất.
  5. Mô phỏng thành công cơ chế truyền sóng bề mặt Love: Xây dựng được bài toán truyền sóng cắt từ tầng nền cứng lên các lớp đất mềm phía trên bằng cách xét đồng thời sóng cắt trong mặt phẳng ngang và mặt phẳng đứng. Lời giải số mô tả chính xác quy luật khuếch đại dao động bề mặt vuông góc với phương truyền sóng, hoàn toàn phù hợp với lý thuyết sóng Love.
  6. Xác định toàn diện nội lực và chuyển vị cọc dưới tải trọng động đất theo thời gian: Thiết lập quy trình phân tích động đất hoàn chỉnh: sử dụng tích phân Duhamel trong miền tần số và chuyển đổi sang miền thời gian bằng IFFT. Sử dụng gia tốc đồ thực tế của trận động đất El Centro (1940) làm đầu vào, xác định chính xác biểu đồ chuyển vị, lực cắt và mô men uốn tại mọi thời điểm và dọc theo toàn bộ chiều dài cọc. Đồng thời, tác giả đã lập trình hoàn chỉnh bộ 9 phần mềm tính toán trên nền tảng MATLAB.

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

Hạn chế của đề tài

Các giới hạn nghiên cứu được xác định rõ theo phạm vi do tác giả công bố:

  • Môi trường đất nền mới chỉ khảo sát trong giới hạn đàn hồi tuyến tính đẳng hướng, chưa xét đến các quy luật ứng xử phi tuyến phức tạp như đàn dẻo hay đàn nhớt.
  • Chưa xét đến hiện tượng hóa lỏng của đất cát bão hòa nước dưới tác động rung động động đất mạnh.
  • Chưa tính toán ảnh hưởng của áp lực nước lỗ rỗng và biến thiên mực nước ngầm trong đất bão hòa.
  • Chưa xét đến trạng thái giới hạn về khả năng chịu lực cực hạn và sự hình thành khớp dẻo trong kết cấu cọc.
  • Nghiên cứu tập trung vào bài toán cọc đơn thẳng đứng, chưa mở rộng giải bài toán đài cọc và nhóm cọc chịu tải trọng động đất phức tạp.

Hướng nghiên cứu tiếp

  • Mở rộng bài toán cho hệ nhóm cọc liên kết với đài móng và kết cấu công trình bên trên (xét tương tác toàn diện đất - móng cọc - kết cấu).
  • Đưa mô hình vật liệu phi tuyến đàn dẻo và tính chất tiêu tán năng lượng phi tuyến của đất vào thuật toán giải.
  • Bổ sung bài toán tương tác thủy cơ nhằm xét đến áp lực nước lỗ rỗng và đánh giá tiềm năng hóa lỏng của nền đất khi chịu chấn động địa chấn.
  • Nghiên cứu thực nghiệm mô hình bàn rung hoặc máy quay ly tâm để đối chứng trực tiếp với các kết quả phần mềm số đã xây dựng.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong các lĩnh vực Địa kỹ thuật công trình, Kỹ thuật xây dựng công trình ngầm, Động lực học công trình và Kỹ thuật kháng chấn:

  • Đối với nghiên cứu sinh và học viên cao học: Cung cấp phương pháp luận chặt chẽ trong việc ứng dụng Phương pháp Nguyên lý cực trị Gauss và phương pháp phần tử hữu hạn 3D để giải quyết các bài toán tương tác môi trường vô hạn mà không cần biên nhân tạo.
  • Đối với giảng viên và nhà khoa học: Là tài liệu tham khảo giá trị về lý thuyết sóng động đất (sóng Love, sóng Rayleigh), lý thuyết dầm Timoshenko, kỹ thuật chuyển đổi miền tần số - thời gian qua FFT/IFFT và tích phân chập Duhamel trong tính toán động lực học kết cấu.
  • Đối với kỹ sư thiết kế công trình ngầm, cầu cảng và công trình đặc biệt: Cung cấp cơ sở lý thuyết và các thuật toán tính toán chuyển vị, mô men uốn, lực cắt của móng cọc khi chịu tải trọng ngang và động đất, khắc phục các điểm thiếu hụt của phương pháp tra bảng hệ số nền thông thường. Hệ thống 9 mã nguồn phần mềm MATLAB trong phần phụ lục là công cụ tham khảo hữu ích cho việc phát triển phần mềm tính toán chuyên dụng.

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

1. Luận án giải quyết vấn đề điều kiện biên vô hạn của bài toán động lực học bằng cách nào?

Luận án sử dụng phương pháp hệ so sánh của PPNLCT Gauss. Đối với bài toán động, tác giả lấy lời giải động lực học xung đơn vị của không gian vô hạn đàn hồi làm hệ so sánh. Nhờ vậy, phương pháp tự động thỏa mãn điều kiện biên bức xạ sóng ra vô cực và điều kiện biên tiếp xúc cọc - đất mà không cần phải đặt các liên kết phụ nhân tạo như lò xo hay hộp nhớt cản sóng trên biên khối đất.

2. Sự khác biệt cốt lõi giữa mô hình cọc trong luận án và mô hình cọc của Poulos (1971) là gì?

Poulos (1971) chỉ xét mô hình 1 chiều, giả thiết đất không bị ảnh hưởng bởi sự có mặt của cọc, bỏ qua ứng suất tiếp ngang giữa đất và mặt bên cọc, đồng thời chưa chập chuyển vị tại đầu và mũi cọc. Ngược lại, luận án xây dựng bài toán 3 chiều đầy đủ, cọc xét theo lý thuyết Timoshenko (có biến dạng trượt ngang), xem xét đồng thời ảnh hưởng tương hỗ hai chiều giữa cọc và đất nền, tự động thỏa mãn liên tục chuyển vị và ứng suất tiếp xúc dọc thân cọc.

3. Tại sao luận án sử dụng hệ số giảm chấn vật liệu thay vì hệ số cản nhớt thông thường cho đất nền?

Hệ số cản nhớt thông thường phụ thuộc trực tiếp vào tần số dao động. Trong khi đó, hệ số giảm chấn vật liệu (hysteretic damping / ma sát khô) $\zeta_h$ phản ánh đúng bản chất ma sát nội và sự tiêu tán năng lượng cơ học trong khối đất, không phụ thuộc vào tần số kích thích, đồng thời tạo thuận lợi để mở rộng bài toán sang trạng thái biến dạng dẻo khi cần thiết.

4. Luận án đã sử dụng số liệu trận động đất nào để kiểm chứng phản ứng động lực học của cọc?

Luận án sử dụng dữ liệu gia tốc đồ thực tế theo phương Bắc - Nam của trận động đất El Centro (1940). Dữ liệu gia tốc theo thời gian được chuyển sang miền tần số bằng phép biến đổi Fourier rời rạc (DFT/FFT), phản ứng của hệ được giải bằng tích phân chập Duhamel, sau đó dùng IFFT để đưa về miền thời gian, cho phép vẽ biểu đồ chuyển vị ngang, mô men uốn và lực cắt của cọc tại các thời điểm $0{,}56\text{ s}$, $3{,}12\text{ s}$, $8{,}24\text{ s}$, và $18{,}48\text{ s}$.

5. Những phần mềm tính toán nào đã được xây dựng trong khuôn khổ luận án?

Tác giả đã lập trình 9 chương trình tính trên ngôn ngữ MATLAB (trình bày tại Phụ lục Quyển 2):

  • Mstatic1 và Kstatic1: Tính toán khối đất chịu tải trọng tĩnh nằm ngang sử dụng hệ so sánh Mindlin và Kelvin.
  • MstaticP1 và KstaticP1: Tính toán cọc đơn chịu tải trọng tĩnh nằm ngang sử dụng hệ so sánh Mindlin và Kelvin.
  • KstaticPLs: Tính cọc nằm trong nền đàn hồi nhiều lớp chịu tải trọng tĩnh nằm ngang.
  • KdynaS: Tính dao động của khối đất chịu tải trọng động nằm ngang.
  • KdynaL: Khảo sát hiện tượng truyền sóng cắt Love trong nền đất phân tầng.
  • KdynaP: Phân tích cọc đơn chịu tải trọng động nằm ngang điều hòa.
  • KdynaPE: Phân tích phản ứng của cọc chịu tác động của sóng động đất.

Kết luận

Luận án tiến sĩ của tác giả Ngô Quốc Trinh đã ứng dụng thành công Phương pháp Nguyên lý cực trị Gauss kết hợp lý thuyết dầm Timoshenko và phương pháp phần tử hữu hạn 3D để giải quyết toàn diện bài toán tương tác giữa cọc đơn và nền đất dưới tác dụng của tải trọng tĩnh ngang, tải trọng động ngang và sóng động đất. Công trình đóng góp lời giải số cho phép quy đổi từ nghiệm Kelvin sang nghiệm Mindlin, tự động thỏa mãn các điều kiện biên ở vô hạn mà không cần đặt biên hấp thụ hay lò xo nhân tạo. Kết quả nghiên cứu cùng hệ thống 9 chương trình phần mềm MATLAB đã thiết lập là cơ sở khoa học và công cụ tính toán có giá trị thực tiễn cao cho công tác thiết kế kháng chấn móng cọc công trình tại Việt Nam.