Tổng quan về luận án

Công trình nghiên cứu tiến sĩ của tác giả Nguyễn Xuân Vũ, được thực hiện dưới sự hướng dẫn khoa học của PGS. Lương Văn Hải tại Trường Đại học Bách Khoa – Đại học Quốc Gia TP. Hồ Chí Minh (Chuyên ngành: Kỹ thuật xây dựng công trình dân dụng và công nghiệp, Mã số: 62580208), tập trung giải quyết bài toán cốt lõi trong cơ học công trình biển: phân tích ứng xử thủy đàn hồi (hydroelastic response) của kết cấu nổi siêu lớn (Very Large Floating Structures - VLFS) chịu tác động của tải trọng di động. Trong bối cảnh biến đổi khí hậu, nước biển dâng và áp lực mở rộng không gian sống, VLFS nổi lên như một giải pháp đột phá thay thế cho phương pháp lấn biển truyền thống nhờ tính linh hoạt tại vùng nước sâu, giảm thiểu tác động đến hệ sinh thái đáy biển và khả năng triệt tiêu chấn động động đất. Tuy nhiên, đặc thù hình học của VLFS – với chiều dài mặt bằng lên tới hàng trăm hoặc hàng nghìn mét nhưng chiều dày chỉ từ vài mét – khiến độ cứng chống uốn tổng thể rất nhỏ so với kích thước hình học, dẫn đến hiện tượng biến dạng uốn đàn hồi chiếm ưu thế hoàn toàn so với chuyển vị khối cứng.

Khoảng trống nghiên cứu (research gap) trọng yếu được xác định xuất phát từ những giới hạn nghiêm trọng của các phương pháp số hiện hữu. Phương pháp biến đổi Fourier (Fourier Transform Method - FTM) của Nevel (1970) và Nugroho et al. (1999) chỉ hiệu quả đối với các bài toán giải tích tuyến tính đơn giản và bế tắc khi xử lý các hệ động lực học di động nhiều bậc tự do hoặc xét tính phi tuyến tiếp xúc. Trong khi đó, phương pháp phần tử hữu hạn (Finite Element Method - FEM) và phương pháp kết hợp phần tử biên - phần tử hữu hạn (BEM-FEM) kinh điển (Endo & Yago, 1999; Kashiwagi, 2000; Cheng et al., 2001) gắn liền với hệ tọa độ cố định trong không gian. Khi tải trọng di chuyển trên tấm nổi rộng vô hạn, phương pháp BEM-FEM đòi hỏi phải chia lưới một miền tính toán khổng lồ dọc theo toàn bộ quãng đường di chuyển và liên tục cập nhật véc-tơ tải trọng tại từng bước thời gian, gây lãng phí tài nguyên máy tính nghiêm trọng và vấp phải sai số phản xạ sóng tại biên giả tạo.

Luận án thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học trọng tâm:

  1. Câu hỏi 1 (RQ1): Làm thế nào để xây dựng một phương pháp số lai hóa có khả năng mô phỏng tương tác thủy đàn hồi trong vùng nước sâu mà miền tính toán không bị phụ thuộc vào quãng đường di chuyển của tải trọng?
    • Giả thuyết 1 (H1): Việc tái lập phương pháp phần tử biên (BEM) cho miền chất lỏng kết hợp phần tử chuyển động (MEM) cho tấm kết cấu trong cùng một hệ tọa độ di động $(\xi, y, z)$ gắn với tải trọng sẽ loại bỏ hoàn toàn việc cập nhật vị trí tải và tối thiểu hóa lưới phần tử.
  2. Câu hỏi 2 (RQ2): Phương pháp phần tử chuyển động (MEM) có thể mở rộng đồng thời cho cả miền chất lỏng vùng nước nông và tương tác rắn - lỏng hay không?
    • Giả thuyết 2 (H2): Dựa trên lý thuyết nước nông tuyến tính, việc rời rạc hóa đồng thời cả miền tấm và miền nước bằng phần tử chuyển động MEM sẽ mang lại nghiệm động lực học chính xác cao với chi phí tính toán tối thiểu.
  3. Câu hỏi 3 (RQ3): Tính trực hướng (orthotropic properties) và góc định hướng của vật liệu cấu tạo nên VLFS ảnh hưởng như thế nào đến sóng uốn, hệ số khuếch đại động (DAF) và lực cản di chuyển?
    • Giả thuyết 3 (H3): Tỷ số mô-đun đàn hồi $E_1/E_2$ và góc nghiêng trục vật liệu $\theta$ làm thay đổi đáng kể phân bố sóng uốn, vùng lõm thủy đàn hồi và vận tốc tới hạn.
  4. Câu hỏi 4 (RQ4): Giải pháp kết cấu nào có thể làm suy giảm hiệu quả tác động của sóng nước thủy động lực học lên các phương tiện giao thông di chuyển trên bề mặt tấm nổi?
    • Giả thuyết 4 (H4): Cấu trúc tấm hai lớp kết nối bởi lớp đàn hồi Winkler (Winkler-type elastic core) sẽ đóng vai trò như một bộ tiêu tán năng lượng, giảm tải áp lực sóng truyền lên bề mặt tầng trên.

Khung lý thuyết của công trình được định hình vững chắc dựa trên sự giao thoa của lý thuyết tấm mỏng Kirchhoff-Love (đẳng hướng và trực hướng), lý thuyết dòng chảy thế phi xoáy thỏa mãn phương trình vi phân Laplace kết hợp phương trình Bernoulli tuyến tính hóa cho miền nước sâu, lý thuyết nước nông tuyến tính (shallow water/tidal theory), và nền tảng phương pháp phần tử chuyển động (MEM) do Koh (2003) khởi xướng. Nghiên cứu thực hiện khảo sát toàn diện trên mô hình sân bay nổi Mega-Float (MEF-300) với tải trọng hạ cánh của máy bay Boeing 747-400 và hệ tải tập trung di chuyển ở các dải vận tốc từ $20\text{ m/s}$ đến trên $51\text{ m/s}$, mang lại giá trị khoa học đột phá cho kỹ thuật công trình biển hiện đại.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử thủy đàn hồi của kết cấu nổi dưới tác động của tải trọng di động bắt nguồn từ các nghiên cứu cơ bản trên các lớp băng nổi ở vùng hàn đới. Nevel (1970) là người tiên phong sử dụng phương pháp FTM trên mô hình tấm Kirchhoff-Love để đánh giá nội lực của băng nổi. Hướng tiếp cận này tiếp tục được phát triển bởi Nugroho et al. (1999) cho bài toán đáp ứng trạng thái dừng của tấm nổi vô hạn, và bởi Kim & Webster (1998) khi ứng dụng mô hình dầm Bernoulli để xác định lực cản thủy đàn hồi tác dụng lên máy bay trong quá trình cất cánh trên sân bay nổi. Agarwala & Nair (2000) mở rộng FTM cho dầm Timoshenko-Mindlin dài vô hạn, trong khi Matiushina et al. (2016) kết hợp biến đổi Fourier và Laplace để phân tích độ võng và trạng thái ứng suất của tấm băng chịu tải trọng máy bay hạ cánh có gia tốc giảm dần. Điểm chung của trường phái giải tích này là tính chính xác cao đối với bài toán kinh điển nhưng bộc lộ giới hạn nghiêm trọng khi đối mặt với mô hình tải trọng động học phức tạp hoặc vật liệu có tính dị hướng.

Nhằm khắc phục rào cản trên, trường phái phương pháp số dựa trên rời rạc hóa không gian xuất hiện mạnh mẽ. Watanabe & Utsunomiya (1996) cùng Qiu (2001) đã ứng dụng FEM thuần túy để phân tích ứng xử tức thời của tấm nổi hữu hạn. Tuy nhiên, do thể tích chất lỏng của đại dương là vô hạn, việc chia lưới 3D toàn bộ trường lưu chất bằng FEM đòi hỏi số lượng phần tử khổng lồ, dẫn đến hiện tượng "bùng nổ" chi phí tính toán. Để giải quyết nút thắt này, Endo & Yago (1999) và Kashiwagi (2000) đã đề xuất phương pháp lai BEM-FEM, trong đó phương pháp phần tử biên (BEM) sử dụng hàm Green để tích phân phương trình dạng yếu chỉ trên mặt tiếp xúc rắn - lỏng, kết hợp với FEM mô hình hóa tấm kết cấu trong miền thời gian. Dù vậy, phương pháp BEM-FEM truyền thống vẫn xây dựng ma trận hệ thống trong hệ tọa độ cố định, đòi hỏi hàm Green phụ thuộc tần số phải được tính toán lại tại từng bước thời gian (Cheng et al., 2001; Ismail, 2006; An, 2008; Trường, 2011).

Song song đó, trong lĩnh vực công trình đường sắt cao tốc, Koh (2003) đã phát minh ra phương pháp phần tử chuyển động (MEM), tạo nên một cuộc cách mạng mô phỏng khi thiết lập hệ phần tử trong hệ quy chiếu chuyển động cùng vận tốc với tải trọng. Nghiên cứu của Ang et al. (2009), Cao et al. (2010), Xu et al. (2013), Lương et al. (2014) và Cao et al. (2017) đã chứng minh tính ưu việt tuyệt đối của MEM trong việc loại bỏ bài toán biên vô hạn của ray và tấm trên nền đàn hồi Pasternak/Kelvin.

                    ┌────────────────────────────────────────────────────────┐
                    │               KHẢO SÁT THỦY ĐÀN HỒI VLFS               │
                    └───────────────────────────┬────────────────────────────┘
                                                │
                 ┌──────────────────────────────┴──────────────────────────────┐
                 ▼                                                             ▼
┌─────────────────────────────────┐                           ┌─────────────────────────────────┐
│     TRƯỜNG PHÁI GIẢI TÍCH       │                           │       TRƯỜNG PHÁI SỐ HỌC        │
│  - FTM (Nevel 1970)             │                           │  - FEM 3D (Watanabe 1996)       │
│  - Dầm Bernoulli (Kim 1998)     │                           │  - BEM-FEM (Kashiwagi 2000)     │
│  - Fourier-Laplace              │                           │  - Nhược điểm: Lưới cố định,    │
│    (Matiushina 2016)            │                           │    chi phí tính toán bùng nổ    │
└────────────────┬────────────────┘                           └────────────────┬────────────────┘
                 │                                                             │
                 └──────────────────────────────┬──────────────────────────────┘
                                                ▼
                               ┌─────────────────────────────────┐
                               │     ĐỘT PHÁ CỦA LUẬN ÁN         │
                               │  - Vùng nước sâu: BEM-MEM       │
                               │  - Vùng nước nông: MEM-MEM      │
                               │  - Hệ tọa độ di động (ξ, y, z)  │
                               │  - Lớp cản ảo nhân tạo          │
                               │  - Cấu trúc tấm nổi hai lớp     │
                               └─────────────────────────────────┘

Công trình của Nguyễn Xuân Vũ định vị chuẩn xác khoảng trống học thuật: chưa từng có nghiên cứu nào kết hợp BEM và MEM trong hệ tọa độ di động để giải quyết bài toán tương tác thủy đàn hồi VLFS trong vùng nước sâu, cũng như chưa có công trình nào phát triển MEM cho cả miền lưu chất trong điều kiện nước nông. So với hai nghiên cứu quốc tế tiêu biểu là mô hình BEM-FEM miền thời gian của Kashiwagi (2000) và thực nghiệm sóng uốn trên băng của Takizawa (1985, 1988), phương pháp đề xuất trong Luận án đạt bước tiến vượt bậc: giảm thiểu tối đa kích thước lưới phần tử mà không phụ thuộc vào hành trình di chuyển của tải trọng, loại bỏ hoàn toàn nhu cầu cập nhật nút tiếp xúc và triệt tiêu sai số biên nhân tạo.

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

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

Luận án đã tạo ra bước chuyển dịch mô hình tính toán (paradigm shift) từ phương pháp lưới cố định sang phương pháp phần tử chuyển động toàn diện trong cơ học thủy đàn hồi.

Thứ nhất, công trình mở rộng lý thuyết tương tác rắn - lỏng (FSI) cổ điển sang hệ quy chiếu Galileo chuyển động đều với vận tốc $v$ dọc theo trục $x$. Trong hệ tọa độ mới $\xi = x - vt, \eta = y, \zeta = z$, toán tử đạo hàm theo thời gian được chuyển đổi sang đạo hàm toàn phần: $$\frac{\partial}{\partial t} \to \frac{\partial}{\partial t} - v \frac{\partial}{\partial \xi}$$ Quy tắc biến đổi vi phân này dẫn đến sự tái lập toàn bộ phương trình vi phân dao động của tấm Kirchhoff-Love và phương trình Laplace cho thế vận tốc lưu chất $\Phi(\xi, \eta, \zeta, t)$.

Thứ hai, nghiên cứu làm sâu sắc thêm lý thuyết sóng uốn thủy đàn hồi (flexural-gravity waves) của Takizawa (1988) bằng việc mô hình hóa chính xác quan hệ động lực học của vật liệu trực hướng với ma trận độ cứng biến đổi $Q_{ij}(\theta)$ theo góc xoay trục vật liệu $\theta$.

Thứ ba, tác giả đề xuất và chứng minh giải pháp lý thuyết cho mô hình tấm nổi hai lớp liên kết đàn hồi mềm loại Winkler: $$m_t \frac{\partial^2 w_t}{\partial t^2} + D_t \nabla^4 w_t + k_c (w_t - w_b) = p_t(\xi, \eta, t)$$ $$m_b \frac{\partial^2 w_b}{\partial t^2} + D_b \nabla^4 w_b + k_c (w_b - w_t) = p_{fluid}(\xi, \eta, t)$$ Mô hình lý thuyết này chứng minh sự tách biệt cơ học giữa bề mặt tiếp xúc giao thông bên trên ($w_t$) và tấm nổi đáy chịu áp lực sóng thủy động lực ($w_b$).

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

Khung phân tích của Luận án tích hợp 4 cấu phần lý thuyết cốt lõi: (1) Lý thuyết tấm mỏng đàn hồi tuyến tính Kirchhoff-Love (đẳng hướng và trực hướng), (2) Lý thuyết dòng chảy thế tuyến tính hóa Bernoulli-Laplace, (3) Lý thuyết nước nông dạng tuyến tính (Saint-Venant/thủy triều), và (4) Mô hình lớp cản hấp thụ sóng ảo nhân tạo (artificial damping layer).

┌─────────────────────────────────────────────────────────────────────────────────┐
│                    KHUNG PHÂN TÍCH THỦY ĐÀN HỒI BEM-MEM / MEM                   │
└────────────────────────────────────────┬────────────────────────────────────────┘
                                         │
     ┌───────────────────────────────────┼───────────────────────────────────┐
     ▼                                   ▼                                   ▼
┌─────────────────────────┐ ┌─────────────────────────┐ ┌─────────────────────────┐
│   KẾT CẤU TẤM (MEM)     │ │   LƯU CHẤT (BEM / MEM)  │ │   BIÊN HẤP THỤ SÓNG     │
│ - Tấm đẳng hướng        │ │ - Nước sâu: BEM thế năng│ │ - Lớp cản ảo nhân tạo   │
│ - Tấm trực hướng Q_ij   │ │   Laplace trong HQC (ξ) │ │   c(ξ) triệt tiêu năng  │
│ - Tấm 2 lớp Winkler k_c │ │ - Nước nông: MEM        │ │   lượng sóng phản xạ    │
└────────────┬────────────┘ └────────────┬────────────┘ └────────────┬────────────┘
             │                           │                           │
             └───────────────────────────┼───────────────────────────┘
                                         ▼
┌─────────────────────────────────────────────────────────────────────────────────┐
│                     HỆ PHƯƠNG TRÌNH TƯƠNG TÁC RẮN - LỎNG                        │
│               [M] d²U/dt² + [C] dU/dt + [K] U = F(t) - Q_int                    │
│               Giải bằng thuật toán tích phân từng bước Newmark                  │
└─────────────────────────────────────────────────────────────────────────────────┘

Điểm độc đáo nhất trong kỹ thuật số của khung phân tích là việc tích hợp lớp cản ảo $c_{v}(\xi)$ ở hai đầu miền tính toán hữu hạn: $$c_{v}(\xi) = c_0 \left(\frac{|\xi| - L_{core}}{L_{damp}}\right)^2 \quad \text{với } |\xi| > L_{core}$$ Kỹ thuật này cho phép sóng uốn phân kỳ lan truyền ra vô cực bị triệt tiêu hoàn toàn năng lượng trước khi chạm tới biên ngoài, ngăn chặn hiện tượng dội sóng phản xạ mà không làm biến dạng đáp ứng động lực học của tấm tại vị trí đặt tải trọng.

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

Thiết kế nghiên cứu

Luận án theo đuổi thế giới quan thực chứng (positivism) với phương pháp tiếp cận diễn dịch - mô phỏng số tất định (deterministic numerical simulation). Toàn bộ thuật toán được lập trình độc lập và tối ưu hóa trên nền tảng MATLAB, kết hợp kiểm chứng đa tầng với nghiệm giải tích FTM, phần mềm phần tử hữu hạn thương mại ANSYS/ABAQUS, và dữ liệu thực nghiệm hiện trường.

Thiết kế nghiên cứu được phân nhánh rõ ràng theo điều kiện thủy văn:

  • Nhánh 1 (Vùng nước sâu, $H/L \gg 0.1$): Phát triển phương pháp lai BEM-MEM. Tấm kết cấu được rời rạc bằng các phần tử dầm/tấm chuyển động MEM bậc cao, trong khi mặt thoáng chất lỏng bên dưới được rời rạc bằng các phần tử biên BEM trong cùng hệ trục di động $\xi$.
  • Nhánh 2 (Vùng nước nông, $H/L \le 0.05$): Phát triển phương pháp MEM thuần túy đa trường (fully coupled structural-fluid MEM). Cả tấm uốn và trường áp lực lưu chất 2D dạng nông được rời rạc hóa đồng thời bằng các phần tử chuyển động 4 nút.

Quy trình nghiên cứu rigorous

Quy trình thiết lập và giải bài toán tuân thủ tính chuẩn xác nghiêm ngặt theo các giao thức:

  1. Biến đổi hệ quy chiếu: Chuyển đổi toàn bộ phương trình cân bằng vi phân từ hệ Euler cố định $(x, y, z, t)$ sang hệ Lagrange-Euler di động $(\xi, y, z, t)$.
  2. Rời rạc hóa Galerkin: Áp dụng phương pháp phần tử hữu hạn chuyển động Galerkin để thành lập các ma trận phần tử: ma trận khối lượng $[M_e]$, ma trận cản $[C_e]$, ma trận độ cứng $[K_e]$, và ma trận tương tác rắn - lỏng $[Q_e]$.
  3. Khử kỳ dị tích phân phần tử biên: Trong BEM, sử dụng hàm Green cơ bản $G = 1/r + 1/r'$ và kỹ thuật tích phân giải tích kết hợp tích phân Gauss 8 điểm để triệt tiêu kỳ dị toán học khi điểm nguồn trùng điểm quan sát.
  4. Triệt tiêu sóng phản xạ: Tích hợp ma trận cản phụ thuộc tọa độ $[C_{damp}]$ của lớp cản nhân tạo tại hai biên $\xi = \pm L$.
  5. Giải phương trình động lực học: Áp dụng thuật toán tích phân từng bước Newmark-beta ($\gamma = 0.5, \beta = 0.25$) đảm bảo tính ổn định vô điều kiện (unconditionally stable) với bước thời gian kiểm soát $\Delta t \le 0.001\text{ s}$.

Data và phân tích

Mô hình nghiên cứu sử dụng các bộ tham số chuẩn quốc tế từ công trình thực tế của Viện Nghiên cứu Kỹ thuật Mega-Float Nhật Bản (TRAM):

  • Thông số sân bay nổi MEF-300: Chiều dài $L_x = 300\text{ m}$, chiều rộng $L_y = 60\text{ m}$, chiều dày $h_s = 2.0\text{ m}$, khối lượng riêng $\rho_s = 250\text{ kg/m}^2$, độ cứng chống uốn $D = 1.488 \times 10^7\text{ kN}\cdot\text{m}$.
  • Thông số lưu chất: Khối lượng riêng nước biển $\rho_w = 1025\text{ kg/m}^3$, gia tốc trọng trường $g = 9.81\text{ m/s}^2$, độ sâu đáy biển $H$ biến thiên từ $5\text{ m}$ (nước nông) đến $50\text{ m}$ và vô hạn.
  • Thông số tải trọng: Máy bay Boeing 747-400 có tổng tải trọng cất/hạ cánh $P_0 = 3.87 \times 10^5\text{ N}$ đến $3.9 \times 10^6\text{ N}$, độ cứng càng đáp $k_t = 1.0 \times 10^7\text{ N/m}$, hệ số cản càng $c_t = 1.0 \times 10^5\text{ N}\cdot\text{s/m}$, dải vận tốc chuyển động $v = 10 \div 100\text{ m/s}$.
  • Kiểm tra độ hội tụ lưới (Mesh convergence): Khảo sát chi tiết với kích thước phần tử $\Delta \xi$ giảm từ $5.0\text{ m}$ xuống $0.5\text{ m}$, chiều dài miền tính toán $L_{domain}$ từ $100\text{ m}$ đến $800\text{ m}$, xác lập miền tối ưu $L_{core} = 200\text{ m}, L_{damp} = 100\text{ m}$.

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

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

Luận án mang lại 5 phát hiện khoa học đột phá có giá trị cốt lõi:

Thứ nhất, sự thay đổi hình thái sóng uốn theo vận tốc và sự tồn tại của vận tốc tới hạn ($v_{cr}$). Khi vận tốc tải trọng $v < 20\text{ m/s}$ ($v \ll v_{cr}$), hình dạng độ võng của tấm nổi đối xứng quanh vị trí tải trọng, mang tính chất tựa tĩnh (quasi-static). Khi vận tốc tiệm cận vận tốc truyền sóng uốn nhỏ nhất $c_{min}$ ($v \approx 24 \div 37\text{ m/s}$), hai hệ sóng uốn xuất hiện đồng thời ở phía trước và phía sau tải trọng. Khi vận tốc siêu tới hạn ($v = 51\text{ m/s} > v_{cr}$), bước sóng phía trước bị nén ngắn lại trong khi bước sóng phía sau dãn dài ra, và hệ số khuếch đại động (DAF) đạt giá trị cực đại ($DAF > 1.85$).

Thứ hai, phát hiện về khoảng trễ tiếp xúc (lag distance) và vùng lõm thủy đàn hồi. Tâm của vùng lõm uốn lớn nhất không bao giờ trùng với điểm đặt tải trọng khi tải di chuyển với vận tốc cao. Điểm đặt tải luôn nằm lệch về phía trước của phễu võng một khoảng trễ $\delta \xi$, tạo nên một độ dốc liên tục ngay dưới bánh đáp, sinh ra lực cản sóng thủy đàn hồi bổ sung (hydroelastic wave resistance) tác động trực tiếp lên quá trình gia tốc của phương tiện.

Thứ ba, tương quan phi tuyến giữa độ sâu nước $H$ và vận tốc tới hạn $v_{cr}$. Trong điều kiện nước nông ($H \le 10\text{ m}$), vận tốc tới hạn $v_{cr}$ giảm mạnh tỉ lệ thuận với $\sqrt{gH}$. Khi độ sâu giảm, biên độ chuyển vị đứng $w_{max}$ tăng vọt do hiện tượng tắc nghẽn lưu chất cục bộ dưới đáy tấm nổi.

Thứ tư, ảnh hưởng của tính trực hướng vật liệu. Khi tỷ số mô-đun $E_1/E_2$ tăng từ $1.0$ lên $3.0$ dọc theo phương di chuyển tải trọng, độ cứng uốn hiệu dụng tăng làm giảm độ sâu vùng lõm $w_{max}$ tới $32%$, đồng thời đẩy vận tốc tới hạn $v_{cr}$ lên dải vận tốc cao hơn. Khi góc xoay sợi vật liệu $\theta$ thay đổi từ $0^\circ$ đến $90^\circ$, hình dạng vùng lõm uốn chuyển từ elip dọc sang elip ngang, làm thay đổi phân bố phản lực cản sóng lên tới $40%$.

Thứ năm, hiệu năng giảm chấn vượt trội của mô hình tấm hai lớp Winkler. Khi sử dụng kết cấu tấm kép với lớp đệm đàn hồi $k_c = 10^5 \div 10^7\text{ N/m}^3$, độ võng và gia tốc động lực học tại tấm trên (mặt đường) giảm từ $45%$ đến $70%$ so với tấm đáy trực tiếp tiếp xúc với nước, triệt tiêu hiệu quả hiện tượng cộng hưởng thủy đàn hồi lên phương tiện di chuyển.

Vận tốc khảo sát ($v$) Hình thái sóng uốn Hệ số DAF Hiện tượng vật lý nổi bật
$v = 20\text{ m/s}$ (Dưới tới hạn) Sóng cục bộ đối xứng $1.05 - 1.10$ Ứng xử tựa tĩnh, phễu võng nằm ngay dưới tải trọng
$v = 24\text{ m/s}$ (Tiệm cận tới hạn) Xuất hiện 2 hệ sóng trước/sau $1.35 - 1.45$ Bắt đầu xuất hiện khoảng trễ $\delta \xi$ giữa tải và tâm lõm
$v = 37\text{ m/s}$ (Gần tới hạn) Sóng uốn phát triển mạnh $1.60 - 1.75$ Sóng trước ngắn, sóng sau dài; lực cản uốn tăng nhanh
$v = 51\text{ m/s}$ (Vượt tới hạn) Sóng uốn phân kỳ bất đối xứng $1.85 - 2.15$ Đạt đỉnh cộng hưởng DAF; năng lượng sóng lan truyền xa

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung công cụ giải tích số hoàn chỉnh cho lĩnh vực thủy đàn hồi học phi tĩnh, đóng góp nền tảng cho việc chuẩn hóa các tiêu chuẩn thiết kế công trình nổi quy mô lớn.
  • Về mặt phương pháp luận: Thuật toán BEM-MEM và MEM trong hệ tọa độ di động cung cấp giải pháp mẫu mực có thể chuyển giao áp dụng cho nhiều bài toán cơ học tiếp xúc di động khác như đường sắt cao tốc trên nền đất yếu, đường hầm ngập nổi dưới biển (SFT - Submerged Floating Tunnels), hoặc dầm cầu chịu tải trọng xe nặng.
  • Về mặt thực tiễn kỹ thuật: Cung cấp biểu đồ tra cứu vận tốc tới hạn $v_{cr}$ và lực cản sóng cho kỹ sư thiết kế sân bay nổi Mega-Float, giúp tối ưu hóa chiều dày kết cấu, hướng cất/hạ cánh của đường băng so với trục sợi vật liệu, và lựa chọn độ cứng lớp đệm cách chấn.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật và điều kiện biên:

  1. Giới hạn vật lý của chất lỏng: Mô hình chất lỏng được giả định là lý tưởng (không nén được, không nhớt, dòng chảy thế phi xoáy). Do đó, các hiệu ứng tiêu tán năng lượng do nhớt, xoáy cục bộ tại mép tấm nổi và hiện tượng phi tuyến sóng bậc cao (sóng vỡ, bọt khí) chưa được mô hình hóa đầy đủ.
  2. Điều kiện tiếp xúc bánh xe - tấm: Tải trọng hạ cánh máy bay được mô hình hóa là hệ khối lượng - lò xo - giảm chấn chuyển động thẳng đều với giả định duy trì tiếp xúc liên tục, chưa xét đến hiện tượng va đập phi tuyến (wheel bounce/separation) khi bề mặt tấm có độ nhám hoặc gồ ghề ngẫu nhiên.
  3. Hình học biên và địa hình đáy: Nghiên cứu giới hạn trong điều kiện đáy biển bằng phẳng và tấm nổi phẳng tuyệt đối có biên rộng vô hạn hoặc chữ nhật lý tưởng, chưa bao quát địa hình đáy biển biến đổi phức tạp hoặc ảnh hưởng của hệ thống neo giữ (mooring systems) và tường chắn sóng (breakwaters).

Chương trình nghiên cứu 10 năm tiếp theo được vạch ra với 4 hướng mở rộng cụ thể:

  • Mở rộng phương pháp BEM-MEM cho bài toán chất lỏng phi tuyến hoàn toàn kết hợp phương pháp mặt nước tự do (Navier-Stokes VOF/SPH).
  • Nghiên cứu cơ chế tiếp xúc 3D có xét đến hiện tượng trượt, nảy bánh xe máy bay trên bề mặt đường băng nổi chịu sóng biển kết hợp tải trọng di động tăng/giảm tốc độ.
  • Tích hợp mô hình vật liệu composite FGM (Functionally Graded Materials) và kết cấu rỗng sandwich 3D tối ưu hóa topo cho các mô-đun VLFS thế hệ mới.
  • Phát triển phần mềm tính toán thủy đàn hồi thương mại hóa có giao diện người dùng (GUI) phục vụ chuyển giao công nghệ trực tiếp cho các tập đoàn xây dựng công trình biển.

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

Công trình luận án mang lại những tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Khẳng định vị thế tiên phong của cơ học tính toán Việt Nam trên trường quốc tế trong lĩnh vực mô phỏng công trình biển. Phương pháp BEM-MEM và MEM mở ra một nhánh nghiên cứu mới về cơ học môi trường liên tục trong hệ quy chiếu di động, với tiềm năng trích dẫn học thuật cao trên các tạp chí Q1 hàng đầu như Marine Structures, Ocean Engineering, và Applied Ocean Research.
  • Chuyển đổi công nghiệp: Tiết kiệm hàng ngàn giờ tính toán mô phỏng cho các viện thiết kế hàng hải và dầu khí. So với phương pháp BEM-FEM truyền thống, BEM-MEM giúp giảm dung lượng bộ nhớ RAM tiêu thụ hơn $80%$ và rút ngắn thời gian xử lý CPU từ hàng chục giờ xuống vài chục phút cho mỗi kịch bản tải trọng di động dài hạn.
  • Lợi ích chiến lược và xã hội: Cung cấp cơ sở khoa học vững chắc phục vụ chiến lược kinh tế biển quốc gia, thiết kế các đảo nhân tạo nổi, căn cứ hậu cần dịch vụ nghề cá, sân bay nổi bảo vệ chủ quyền biển đảo, và các trang trại điện mặt trời nổi ngoài khơi có tuổi thọ trên 50 năm.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận hệ thống công thức giải tích và mã nguồn thuật toán mẫu mực về BEM, MEM, và tương tác FSI trong hệ quy chiếu chuyển động.
  • Các nhà khoa học cơ học công trình biển: Sở hữu khung lý thuyết chuẩn xác để so sánh, phát triển các mô hình số bậc cao cho kết cấu nổi mềm.
  • Kỹ sư R&D và tư vấn thiết kế công trình biển: Ứng dụng trực tiếp công thức xác định vận tốc tới hạn $v_{cr}$, hệ số DAF và lực cản sóng uốn để tính toán kiểm tra độ bền mỏi và an toàn khai thác của sân bay nổi, cầu phao và đường băng nổi.
  • Các nhà hoạch định chính sách hạ tầng biển: Có cơ sở khoa học để phê duyệt các dự án đô thị nổi, cảng biển nổi thông minh thay thế giải pháp quai đê lấn biển gây tổn hại môi trường sinh thái.

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 độc đáo nhất là việc mở rộng phương pháp phần tử biên (BEM) và phương pháp phần tử chuyển động (MEM) từ hệ quy chiếu tĩnh sang hệ quy chiếu di động đồng bộ $(\xi, y, z)$ cho bài toán tương tác thủy đàn hồi 3D. Luận án đã mở rộng lý thuyết dòng chảy thế Laplace và lý thuyết tấm mỏng Kirchhoff-Love cổ điển, tích hợp thành công với lý thuyết cản ảo nhân tạo để giải quyết dứt điểm bài toán tấm vô hạn chịu tải di động mà không cần chia lưới toàn bộ hành trình tải.

2. Đột phá phương pháp luận của nghiên cứu khi so sánh với ít nhất hai công trình quốc tế kinh điển?

So với phương pháp biến đổi Fourier (FTM) của Nevel (1970) (vốn bất lực trước hệ động lực học phi tuyến nhiều bậc tự do) và phương pháp BEM-FEM miền thời gian của Kashiwagi (2000) (đòi hỏi lưới phần tử cố định khổng lồ và phải tính lại hàm Green theo thời gian), phương pháp BEM-MEM của Luận án chỉ cần tích phân hàm Green cơ bản một lần duy nhất, giải phóng hoàn toàn sự phụ thuộc giữa kích thước miền tính toán và quãng đường di chuyển của tải trọng, giúp giảm trên $80%$ tài nguyên tính toán máy tính.

3. Phát hiện thực nghiệm/mô phỏng gây bất ngờ nhất có số liệu minh chứng?

Phát hiện bất ngờ nhất là hiện tượng "khoảng trễ tiếp xúc" (lag distance) và đỉnh lực cản sóng uốn tại vận tốc siêu tới hạn $v = 51\text{ m/s}$. Tại dải vận tốc này, hệ số khuếch đại động $DAF$ tăng vọt lên mức $1.85 - 2.15$, và tâm vùng lõm bị tụt lại phía sau bánh đáp máy bay, tạo ra một dốc uốn nhân tạo liên tục buộc càng hạ cánh phải chịu phản lực cản uốn tương đương việc máy bay đang phải "leo dốc" trên mặt nước phẳng.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ. Luận án trình bày chi tiết toàn bộ các phương trình biến suy rộng, ma trận độ cứng phần tử $[K_e]$, ma trận cản $[C_e]$, ma trận khối lượng $[M_e]$, ma trận tương tác $[Q_e]$, công thức hàm Green, thông số lớp cản ảo $c_0, L_{core}, L_{damp}$, sơ đồ thuật toán Newmark từng bước, kèm theo đầy đủ bảng thông số hình học và cơ lý của kết cấu sân bay MEF-300 và máy bay Boeing 747-400.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào 3 giai đoạn: (1) Giai đoạn 1-3 năm: Mở rộng BEM-MEM cho mô hình chất lỏng phi tuyến Navier-Stokes và kết cấu tấm FGM; (2) Giai đoạn 4-6 năm: Phát triển bài toán tương tác va đập phi tuyến 3D giữa càng đáp máy bay bị trượt/nảy trên mặt đường băng nổi chịu kích động sóng biển ngẫu nhiên; (3) Giai đoạn 7-10 năm: Xây dựng bộ tiêu chuẩn thiết kế quốc gia và phần mềm chuyên dụng phục vụ công nghiệp hóa các đô thị nổi và sân bay nổi chiến lược tại Biển Đông.

Kết luận

  1. Xây dựng thành công phương pháp lai mới BEM-MEM trong hệ tọa độ di động, giải quyết triệt để bài toán phân tích thủy đàn hồi của kết cấu nổi siêu lớn (VLFS) trong vùng nước sâu chịu tải trọng di động.
  2. Phát triển toàn diện phương pháp phần tử chuyển động (MEM) cho cả miền tấm và miền lưu chất trong điều kiện nước nông tuyến tính, giúp tối ưu hóa tuyệt đối tài nguyên bộ nhớ và thời gian tính toán.
  3. Tích hợp thành công lớp cản ảo nhân tạo phụ thuộc không gian vào hệ phương trình động lực học, triệt tiêu hoàn toàn hiện tượng phản xạ sóng uốn giả tạo tại biên của miền tính toán hữu hạn.
  4. Làm sáng tỏ quy luật biến đổi hình thái sóng uốn, hệ số khuếch đại động DAF, khoảng trễ tiếp xúc và lực cản sóng uốn dưới ảnh hưởng của vận tốc tải trọng, độ sâu nước, tỷ số trực hướng $E_1/E_2$, và góc nghiêng sợi vật liệu $\theta$.
  5. Đề xuất và chứng minh tính ưu việt của cấu trúc tấm hai lớp liên kết đàn hồi Winkler trong việc triệt tiêu từ $45%$ đến $70%$ biên độ dao động và áp lực thủy động truyền lên bề mặt tầng trên.
  6. Đóng góp một di sản học thuật quan trọng với hệ thống công cụ số có độ tin cậy cao, mở đường cho việc ứng dụng kết cấu nổi siêu lớn trong chiến lược khai thác đại dương bền vững và bảo vệ chủ quyền biển đảo thế kỷ 21.