Tổng quan về luận án

Sự gia tăng dân số toàn cầu cùng tốc độ đô thị hóa nhanh chóng tại các quốc gia có đường bờ biển dài đã tạo ra áp lực khổng lồ lên quỹ đất đô thị và hạ tầng giao thông ven biển. Trong bối cảnh giải pháp lấn biển truyền thống chỉ khả thi về mặt kỹ thuật và kinh tế ở những vùng nước nông (độ sâu dưới 20 m), đồng thời gây ra những hệ lụy sinh thái nghiêm trọng như phá hủy hệ sinh thái ngập mặn và biến đổi đường bờ biển, kết cấu nổi siêu lớn (Very Large Floating Structures - VLFS) dạng tấm mỏng (Mega-Float) và dạng bán chìm (Mobile Offshore Base - MOB) đã nổi lên như một giải pháp công trình tiên phong mang tính chiến lược. VLFS sở hữu những ưu thế vượt trội: thi công nhanh theo dạng mô-đun hóa tiền chế, dễ dàng mở rộng hoặc tháo dỡ, triệt tiêu hoàn toàn tác động từ động đất truyền qua nền đất và thích ứng tự nhiên với hiện tượng mực nước biển dâng do biến đổi khí hậu toàn cầu.

Tuy nhiên, bài toán ứng xử thủy đàn hồi (hydroelasticity) của VLFS dưới tác động của tải trọng di động—điển hình là quá trình cất/hạ cánh của máy bay phản lực thương mại hạng nặng (như Boeing 747) hay các phương tiện giao thông vận tải—đang đặt ra thách thức mô phỏng số to lớn. Về mặt vật lý, do tỷ số giữa chiều dày ($h_s \approx 1\text{--}5\text{ m}$) và kích thước mặt bằng ($L, B$ lên đến hàng trăm hoặc hàng nghìn mét) rất nhỏ, kết cấu VLFS chịu biến dạng uốn đàn hồi chiếm ưu thế tuyệt đối so với chuyển vị khối cứng. Nền nước có độ cứng đẩy nổi nhỏ hơn rất nhiều so với nền đất đàn hồi, dẫn đến việc sóng uốn (flexural waves) hình thành dưới điểm đặt tải lan truyền với biên độ lớn và vươn xa trên một miền không gian rộng lớn trước khi suy giảm.

+-------------------------------------------------------------------------+
|                  HỆ THỐNG PHƯƠNG PHÁP NGHIÊN CỨU TRONG LUẬN ÁN          |
+-------------------------------------------------------------------------+
|                                                                         |
|      +-----------------------------------------------------------+      |
|      |               KẾT CẤU NỔI SIÊU LỚN (VLFS)                 |      |
|      |           Chịu Tải Trọng Di Động (Hệ Động Lực Học)        |      |
|      +-----------------------------+-----------------------------+      |
|                                    |                                    |
|                  +-----------------+-----------------+                  |
|                  |                                   |                  |
|                  v                                   v                  |
|         [ VÙNG NƯỚC SÂU / HỮU HẠN ]           [ VÙNG NƯỚC NÔNG ]        |
|                  |                                   |                  |
|                  v                                   v                  |
|         +-----------------+                 +-----------------+         |
|         | PHƯƠNG PHÁP     |                 | PHƯƠNG PHÁP     |         |
|         | BEM - MEM       |                 | THUẦN MEM       |         |
|         +-----------------+                 +-----------------+         |
|         | - BEM động chất |                 | - MEM cho tấm   |         |
|         |   lưu (Hệ HĐĐ)  |                 | - MEM cho nước  |         |
|         | - MEM cho tấm   |                 |   nông (Thủy    |         |
|         | - Ma trận tương |                 |   triều tuyến   |         |
|         |   tác di động   |                 |   tính)         |         |
|         +--------+--------+                 +--------+--------+         |
|                  |                                   |                  |
|                  +-----------------+-----------------+                  |
|                                    |                                    |
|                                    v                                    |
|      +-----------------------------------------------------------+      |
|      | - Triệt tiêu sóng phản xạ bằng Lớp Cản Ảo (Artificial     |      |
|      |   Damping Layer)                                          |      |
|      | - Tích phân Newmark trực tiếp theo thời gian              |      |
|      | - Giải pháp tấm 2 lớp lõi đàn hồi Winkler giảm chấn sóng  |      |
|      +-----------------------------------------------------------+      |
+-------------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong y văn quốc tế nằm ở sự hạn chế nghiêm trọng về mặt công cụ tính toán số:

  • Phương pháp biến đổi Fourier (Fourier Transform Method - FTM) chỉ giải quyết hiệu quả các mô hình tải trọng đơn giản chuyển động đều trên tấm vô hạn, hoàn toàn bất khả thi khi mở rộng cho các hệ động lực học phức tạp nhiều bậc tự do, xét đến tương tác phi tuyến lốp xe - mặt đường hoặc tải trọng có vận tốc biến thiên.
  • Phương pháp Phần tử Hữu hạn (FEM) và phương pháp ghép Phần tử Biên - Phần tử Hữu hạn truyền thống (BEM-FEM) trong hệ tọa độ Descartes cố định đòi hỏi miền chia lưới chất lưu và kết cấu cực lớn tương ứng với toàn bộ chiều dài đường chạy của tải trọng. Điều này dẫn đến sự bùng nổ tài nguyên tính toán (RAM, CPU time), đồng thời liên tục phải cập nhật vị trí điểm tiếp xúc tại từng bước thời gian và gặp phải sai số nghiêm trọng khi tải trọng di chuyển gần các vùng biên bị cắt cụt nhân tạo.

Luận án tiến sĩ của NCS. Nguyễn Xuân Vũ, dưới sự hướng dẫn khoa học của PGS. TS. 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 (2020), mang tên "Phát triển phương pháp phần tử biên và phần tử chuyển động trong phân tích ứng xử kết cấu nổi siêu lớn" (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) đã giải quyết triệt để các nút thắt trên thông qua hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để loại bỏ sự phụ thuộc kích thước lưới tính toán vào quãng đường di chuyển của tải trọng mà vẫn bảo đảm tính toán chính xác trường áp lực thủy động lực học trong miền nước sâu?
  • Câu hỏi nghiên cứu 2 (RQ2): Có thể thiết lập một mô hình phần tử chuyển động thuần túy (Pure MEM) đồng thời cho cả kết cấu tấm và miền chất lưu bên dưới trong điều kiện nước nông để đơn giản hóa quá trình tương tác rắn - lỏng hay không?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật ảnh hưởng của tính trực hướng vật liệu (orthotropy), góc xoay trục đàn hồi và độ sâu lớp nước đáy biển đến hệ số khuếch đại động (DAF), hình dạng vùng lõm biến dạng và lực cản cảm ứng lên phương tiện diễn ra như thế nào?
  • Câu hỏi nghiên cứu 4 (RQ4): Giải pháp kết cấu nào có khả năng triệt tiêu hoặc giảm thiểu hiệu quả tác động bất lợi của sóng nước đẩy ngược lên sàn xe/sân bay di chuyển bên trên?

Hệ thống giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Việc tái lập phương trình BEM cho chất lưu kết hợp MEM cho tấm trong hệ tọa độ chuyển động đồng bộ $(r, s, z)$ sẽ cố định hóa điểm đặt tải tại gốc tọa độ cục bộ, triệt tiêu nhu cầu cập nhật véc-tơ tải và giảm 60–80% khối lượng tính toán ma trận.
  • Giả thuyết 2 (H2): Dưới giả định nước nông tuyến tính (thuyết thủy triều), miền chất lưu có thể được rời rạc hóa trực tiếp bằng các phần tử chuyển động 2D kết nối trực tiếp với tấm thông qua ma trận tương tác rắn - lỏng di động.
  • Giả thuyết 3 (H3): Tồn tại một dải vận tốc tới hạn $V_{\text{cr}} \approx C_{\min}$ gây ra hiện tượng cộng hưởng thủy đàn hồi cực đại, và tính trực hướng vật liệu có thể định hướng lại sự phân bố năng lượng sóng uốn để làm giảm lực cản di động.
  • Giả thuyết 4 (H4): Mô hình tấm hai lớp liên kết bằng lõi đàn hồi Winkler đóng vai trò như một bộ lọc cơ học thông dải, hấp thụ áp lực động của nước và cách ly dao động cho bản mặt cầu phía trên.

Khung lý thuyết của công trình tích hợp đa tầng từ: 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, không nén được (phương trình Laplace và tích phân Bernoulli); Lý thuyết sóng nước nông dạng tuyến tính (thuyết sóng dài thủy triều); Lý thuyết nền đàn hồi Winkler; và Phương pháp Phần tử Chuyển động (Moving Element Method - MEM). Phạm vi nghiên cứu bao quát kết cấu tấm vô hạn nổi trên nước có độ sâu hữu hạn và nước nông phẳng đáy, chịu tác động của tải trọng tập trung hoặc hệ động lực học chuyển động thẳng đều với dải vận tốc từ dưới tới hạn đến siêu tới hạn ($V/C_{\min} = 0.2\text{--}1.5$).

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 siêu lớn chịu tải trọng di động đã trải qua bốn nhánh tiến hóa học thuật chính:

Nhánh thứ nhất khởi nguồn từ các phân tích giải tích sử dụng phép biến đổi không gian Fourier (FTM). Tiêu biểu như Nevel (1970) khi nghiên cứu nội lực tấm băng nổi vô hạn chịu tải di động; Kim và Webster (1998) phát triển mô hình dầm Bernoulli để dự báo lực cản sinh ra do độ võng của đường băng nổi tác dụng lên máy bay khi cất cánh; Agarwala và Nair (2006) áp dụng FTM cho dầm Timoshenko-Mindlin; Matiushina et al. (2016) kết hợp biến đổi Fourier và Laplace xét đến yếu tố giảm tốc khi máy bay hạ cánh; Hosking và Sneyd (1986) làm rõ mối quan hệ giữa vận tốc tải trọng và sự phân nhánh của các hệ sóng uốn. Điểm chung của dòng nghiên cứu FTM là giải quyết bài toán trong không gian phổ tuần hoàn, sau đó chuyển đổi ngược về miền không gian - thời gian thực. Tuy nhiên, FTM hoàn toàn bất lực trước các bài toán biên phức tạp, sự biến thiên độ cứng cục bộ hoặc hệ động lực học bánh xe phi tuyến.

Nhánh thứ hai tập trung vào việc số hóa miền tính toán thông qua phương pháp phần tử hữu hạn (FEM) và phương pháp ghép BEM-FEM. Watanabe và Utsunomiya (2003) phân tích tấm chữ nhật nổi chịu tải xung và tải di động bằng FEM; Qiu (2007) khảo sát tấm nổi trên nước sâu hữu hạn; Endo và Yago (1999) đề xuất mô hình BEM-FEM trong miền tần số rồi chuyển đổi về miền thời gian; Kashiwagi (2000) áp dụng phương pháp tổ hợp mode dao động trong miền thời gian để tính toán lực cản cất cánh của Boeing 747; Cheng et al. (2001) dùng hàm Green phụ thuộc thời gian; Li et al. (2003) và Ismail (2005) sử dụng hàm Green cơ bản độc lập với tần số; An (2010), Trương (2013), Thái (2014) lần lượt mở rộng BEM-FEM cho sân bay nổi có liên kết đàn hồi với đáy biển. Nhược điểm chí mạng của nhánh này là việc rời rạc hóa toàn bộ thể tích chất lưu (FEM) hoặc diện tích đáy kết cấu cố định (BEM) đòi hỏi số lượng phần tử khổng lồ, khiến thời gian tính toán tăng theo hàm mũ khi vận tốc tải trọng cao.

Nhánh thứ ba đánh dấu bước đột phá về mặt động học kết cấu khi Koh et al. (2003, 2006) phát minh phương pháp phần tử chuyển động (MEM) nhằm phân tích động lực học đường ray tàu cao tốc và tấm trên nền viscoelastic. Tiếp đó, Ang et al. (2011), Cao et al. (2010, 2018), Lei và Wang (2014), Xu et al. (2016), Lương et al. (2015) đã mở rộng MEM cho tấm Mindlin và tấm composite trực hướng trên nền Pasternak. Ưu điểm nổi trội của MEM là các phần tử được thiết lập trong hệ quy chiếu chuyển động cùng vận tốc $V$ với tải trọng. Tuy nhiên, toàn bộ các công trình trên chỉ giới hạn MEM cho kết cấu dầm/tấm thuần túy trên nền đàn hồi cơ học, chưa từng tích hợp MEM với phương trình chất lưu thủy động lực học có xét đến hiệu ứng khối lượng cộng thêm (added mass) và cản bức xạ (radiation damping).

Nhánh thứ tư nghiên cứu về hệ kết cấu nhiều lớp liên kết đàn hồi, khởi đầu từ Seelig và Hoppmann (1964), Oniszczuk (2000), Abu-Hilal (2006) và Mirzabeigy et al. (2016) trên các hệ dầm song song nối bởi tầng đàn hồi Winkler. Song, các mô hình này mới chỉ khảo sát dao động trong không khí, hoàn toàn bỏ qua tương tác động lực với môi trường chất lưu xung quanh.

+-------------------------------------------------------------------------+
|                  ĐỐI SÁNH CÁC PHƯƠNG PHÁP MÔ PHỎNG VLFS                 |
+-------------------+-----------------+-----------------+-----------------+
| Tiêu chí so sánh  | Phương pháp     | BEM - FEM       | BEM - MEM & MEM |
|                   | FTM Cổ điển     | Truyền thống    | (Luận án đề     |
|                   |                 |                 | xuất)           |
+-------------------+-----------------+-----------------+-----------------+
| 1. Hệ tọa độ      | Cố định         | Cố định         | Chuyển động     |
|    tính toán      | (Không gian phổ)| (Eulerian)      | đồng bộ (Moving)|
+-------------------+-----------------+-----------------+-----------------+
| 2. Kích thước     | Không áp dụng   | Phụ thuộc toàn  | Thu hẹp cố định |
|    lưới số        |                 | bộ đường chạy   | quanh tải trọng |
+-------------------+-----------------+-----------------+-----------------+
| 3. Cập nhật vị trí| Bất khả thi     | Bắt buộc tại    | Tải cố định tại |
|    tải trọng      | với tải phức tạp| từng bước $\Delta t$| gốc tọa độ      |
+-------------------+-----------------+-----------------+-----------------+
| 4. Xử lý biên cắt | Không xét       | Dễ nhiễu phản xạ| Lớp cản ảo nhân |
|    cụt (Truncation)                 | sóng biên       | tạo hấp thụ 100%|
+-------------------+-----------------+-----------------+-----------------+
| 5. Tải phi tuyến/ | Rất hạn chế     | Khả thi nhưng   | Khả thi, tối ưu |
|    Nhiều bậc TD   |                 | rất chậm        | tài nguyên      |
+-------------------+-----------------+-----------------+-----------------+
| 6. Tương tác thủy | Chất lưu lý     | BEM tần số hoặc | BEM di động hoặc|
|    động lực học   | tưởng vô hạn    | Green thời gian | MEM thủy triều  |
+-------------------+-----------------+-----------------+-----------------+

Luận án của Nguyễn Xuân Vũ đã định vị xuất sắc đóng góp khoa học tại điểm giao thoa giữa cơ học công trình biển và kỹ thuật mô phỏng tính toán tiên tiến. Công trình được đối chuẩn trực tiếp với hai nghiên cứu quốc tế kinh điển:

  1. Thí nghiệm hiện trường quy mô thực của Takizawa (1985, 1988) trên tấm băng nổi tại hồ Saroma (Nhật Bản) chịu tải trọng xe trượt tuyết di động để kiểm chứng quy luật lan truyền sóng uốn, vùng lõm biến dạng và hiện tượng trễ không gian;
  2. Dữ liệu tính toán sân bay nổi quy mô lớn Mega-Float MF-300 (chiều dài $L = 300\text{ m}$, rộng $B = 60\text{ m}$, chiều dày tương đương $h_s = 2\text{ m}$) của Hiệp hội Nghiên cứu Kỹ thuật Mega-Float Nhật Bản (TRAM) khi chịu máy bay hạ cánh.

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

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

Luận án mang lại những bước tiến lý thuyết đột phá thông qua việc thiết lập nền tảng toán học - cơ học mới cho bài toán thủy đàn hồi phụ thuộc thời gian:

Thứ nhất, thiết lập thành công hệ phương trình vi phân chuyển động thủy đàn hồi hoàn chỉnh trong hệ tọa độ di động Descartes $(r, s, z)$ tịnh tiến dọc trục $x$ với vận tốc không đổi $V$: $$r = x - Vt, \quad s = y, \quad z = z$$

Khi thực hiện phép biến đổi tọa độ Galois, các đạo hàm riêng theo thời gian $\frac{\partial}{\partial t}$ được thay thế bằng toán tử vi phân đối lưu toàn phần: $$\frac{\text{d}}{\text{d}t} = \frac{\partial}{\partial t} - V\frac{\partial}{\partial r}$$

Dưới hệ tọa độ mới, phương trình vi phân dao động của tấm mỏng trực hướng (Orthotropic Plate) nổi trên mặt nước chịu tác dụng của trường áp lực thủy động $p_z(r, s, t)$ và tải trọng di động $P_0$ có dạng toán học chuẩn tắc:

$$\begin{aligned} &\frac{1}{6}h_s^3 \left[ Q_{11}\frac{\partial^4 w}{\partial r^4} + 4Q_{16}\frac{\partial^4 w}{\partial r^3 \partial s} + 2(Q_{12} + 2Q_{66})\frac{\partial^4 w}{\partial r^2 \partial s^2} + 4Q_{26}\frac{\partial^4 w}{\partial r \partial s^3} + Q_{22}\frac{\partial^4 w}{\partial s^4} \right] \ &+ m\left( \frac{\partial^2 w}{\partial t^2} - 2V\frac{\partial^2 w}{\partial r \partial t} + V^2\frac{\partial^2 w}{\partial r^2} \right) + c_f\left( \frac{\partial w}{\partial t} - V\frac{\partial w}{\partial r} \right) + a\nabla^4 w = p_z(r, s, t) \end{aligned}$$

Trong đó, các hệ số độ cứng đàn hồi biến đổi $Q_{ij}$ phụ thuộc trực tiếp vào góc nghiêng trục sợi $\theta$ so với phương chuyển động thông qua định luật biến đổi tensor ứng suất bậc 4 của Reddy (2004).

Thứ hai, mở rộng phương trình thế vận tốc chất lưu $\Phi(r, s, z, t)$ thỏa mãn phương trình Laplace $\nabla^2\Phi = 0$ và phương trình tích phân áp lực Bernoulli tuyến tính hóa tại mặt tiếp xúc rắn - lỏng $z = 0$: $$p_z(r, s, t) = -\rho\left( \frac{\partial\Phi}{\partial t} - V\frac{\partial\Phi}{\partial r} + g w(r, s, t) \right)$$

Thứ ba, phát triển hệ phương trình vi phân cặp đôi cho mô hình kết cấu nổi hai lớp liên kết qua tầng đàn hồi Winkler có độ cứng $k_c$: $$\begin{aligned} m_t \frac{\text{d}^2 w_t}{\text{d}t^2} + a\nabla^4 w_t + D_t \nabla^4 w_t + k_c(w_t - w_b) &= P_0 \delta(r)\delta(s) \ m_b \frac{\text{d}^2 w_b}{\text{d}t^2} + a\nabla^4 w_b + D_b \nabla^4 w_b - k_c(w_t - w_b) &= -\rho\left( \frac{\text{d}\Phi}{\text{d}t} + gw_b \right) \end{aligned}$$

Mô hình này tạo ra sự chuyển dịch hệ hình (paradigm shift) khi chứng minh rằng lớp lõi đàn hồi đóng vai trò như một bộ cản chấn thủy động học (hydrodynamic damper), phân tán năng lượng sóng uốn và bảo vệ tấm mặt trên $w_t$ khỏi các xung kích áp lực nước $w_b$.

+-------------------------------------------------------------------------+
|                CƠ CHẾ GIẢM CHẤN CỦA MÔ HÌNH TẤM HAI LỚP                |
+-------------------------------------------------------------------------+
|                                                                         |
|        Tải Trọng Di Động P(t) [Máy bay / Xe cộ chuyển động]             |
|                          |                                              |
|                          v                                              |
|     +-------------------------------------------------------------+     |
|     |  TẤM TRÊN (Top Plate: h_t, E_t, D_t) -> Chuyển vị w_t GIẢM  |     |
|     +-------------------------------------------------------------+     |
|                          |                                              |
|                          v                                              |
|     +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~+     |
|     |  LỚP ĐỆM ĐÀN HỒI WINKLER (Elastic Core: Độ cứng k_c)        |     |
|     |  [Hấp thụ xung áp lực, triệt tiêu truyền chấn động lực]     |     |
|     +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~+     |
|                          |                                              |
|                          v                                              |
|     +-------------------------------------------------------------+     |
|     |  TẤM DƯỚI (Bottom Plate: h_b, E_b, D_b) -> Chịu uốn w_b     |     |
|     +-------------------------------------------------------------+     |
|                          |                                              |
|                          v  Áp lực đẩy nổi & quán tính chất lưu         |
|     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~     |
|                   VÙNG CHẤT LƯU BÊN DƯỚI (Nước biển rho)                |
+-------------------------------------------------------------------------+

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

Khung phân tích số của luận án kết hợp ba kỹ thuật tính toán độc nhất:

  1. Phương pháp BEM-MEM trong vùng nước sâu: Sử dụng phần tử chuyển động tấm kết hợp ma trận chuyển đổi biên chất lưu $[G_{dr}]$ và $[H]$. Toàn bộ hệ phương trình đại số tuyến tính được biểu diễn dưới dạng ma trận khối: $$\begin{bmatrix} [M] & [0] \ [0] & [0] \end{bmatrix} \begin{Bmatrix} {\ddot{w}_e} \ {\ddot{\Phi}_e} \end{Bmatrix} + \begin{bmatrix} [C] & [0] \ [0] & [0] \end{bmatrix} \begin{Bmatrix} {\dot{w}_e} \ {\dot{\Phi}_e} \end{Bmatrix} + \begin{bmatrix} [K] & [Q_1] \ [Q_2] & [H] \end{bmatrix} \begin{Bmatrix} {w_e} \ {\Phi_e} \end{Bmatrix} = \begin{Bmatrix} {P} \ {0} \end{Bmatrix}$$ Phương pháp này trích xuất trực tiếp chuyển vị thẳng đứng thông qua ma trận biến đổi $[L_1]$ và ma trận tải sóng $[L_2]$.

  2. Phương pháp thuần MEM cho vùng nước nông: Tận dụng lý thuyết sóng dài thủy triều: $$\frac{\partial \eta}{\partial t} + H \left( \frac{\partial u_f}{\partial x} + \frac{\partial v_f}{\partial y} \right) = 0$$ Rời rạc hóa đồng thời cả trường vận tốc chất lưu đáy và trường chuyển vị tấm bằng các hàm dạng phần tử chuyển động 4 nút bậc một, thiết lập ma trận chất lưu nước nông $[H_{sh}]$ và ma trận tương tác rắn - lỏng $[L_{dr2}]$.

  3. Lớp cản ảo nhân tạo (Artificial Damping Layer): Để giải quyết triệt để sự phản xạ sóng uốn tại hai đầu biên bị cắt cụt của miền tính toán hữu hạn $[-L_1, L_2]$, một lớp cản ảo phi tuyến có hệ số cản biến thiên $c_{vd}(r)$ được tích hợp vào phương trình tấm: $$c_{vd}(r) = \begin{cases} c_0 \left(\frac{|r| - L_{\text{core}}}{L_{\text{damp}}}\right)^2 & \text{khi } |r| > L_{\text{core}} \ 0 & \text{khi } |r| \le L_{\text{core}} \end{cases}$$ Kỹ thuật này triệt tiêu hoàn toàn năng lượng sóng phản xạ mà không làm sai lệch đáp ứng thủy đàn hồi cục bộ tại vùng lân cận điểm đặt tải.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt triết lý thực chứng (positivism) của cơ học tính toán (computational mechanics), vận hành qua quy trình kiểm chứng chéo ba trục (triangulation):

  • Trục 1: So sánh giải tích với nghiệm FTM kinh điển;
  • Trục 2: So sánh số với phần mềm thương mại và phương pháp BEM-FEM truyền thống;
  • Trục 3: Đối chiếu thực nghiệm vật lý với các bài đo đạc hiện trường của Takizawa.
+-------------------------------------------------------------------------+
|                  QUY TRÌNH KIỂM CHỨNG ĐỘ TIN CẬY (TRIANGULATION)        |
+-------------------------------------------------------------------------+
|                                                                         |
|                     [ KẾT QUẢ TỪ BEM-MEM / MEM ]                        |
|                                  |                                      |
|            +---------------------+---------------------+                |
|            |                                           |                |
|            v                                           v                |
|   [ TRỤC GIẢI TÍCH SỐ ]                     [ TRỤC THỰC NGHIỆM ]        |
|   - Nghiệm biến đổi Fourier (FTM)           - Thực nghiệm Takizawa      |
|   - Phương pháp BEM-FEM cố định               trên tấm băng tự nhiên    |
|   - Khảo sát hội tụ lưới & bước thời gian   - Thông số sân bay nổi      |
|     Delta t (Newmark-beta)                    Mega-Float MF-300         |
|            |                                           |                |
|            +---------------------+---------------------+                |
|                                  |                                      |
|                                  v                                      |
|               [ ĐÁNH GIÁ ĐỘ CHÍNH XÁC VÀ TIN CẬY CAO ]                  |
|                 - Sai số tương đối < 1.5%                               |
|                 - Tiết kiệm 60-80% thời gian tính toán                  |
+-------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Toàn bộ thuật toán được lập trình độc lập bằng ngôn ngữ tính toán ma trận khoa học MATLAB. Quy trình tính toán động lực học phi tuyến theo thời gian diễn ra qua các giai đoạn:

  1. Thiết lập lưới phần tử chuyển động cục bộ đồng chuyển động với tải trọng;
  2. Tích phân Gauss số học trên từng phần tử tấm và biên chất lưu để hình thành các ma trận độ cứng $[K_e]$, ma trận khối lượng $[M_e]$, ma trận cản $[C_e]$, ma trận tương tác $[Q_e]$ và véc-tơ tải ${P_e}$;
  3. Ghép nối thành các ma trận tổng thể có băng thông hẹp (banded sparse matrices);
  4. Áp dụng thuật toán tích phân Newmark trực tiếp theo thời gian với các tham số vô điều kiện ổn định ($\gamma = 0.5$, $\beta = 0.25$);
  5. Giải hệ phương trình đại số tại từng bước thời gian $\Delta t$ để thu được đồng thời chuyển vị thẳng đứng $w(r, s, t)$, góc xoay và áp suất đáy chất lưu.

Data và phân tích

Mô hình số sử dụng hệ thông số kỹ thuật tiêu chuẩn của sân bay nổi Mega-Float MF-300:

  • Kích thước hình học tấm: $L = 300\text{ m}$, $B = 60\text{ m}$, chiều dày $h_s = 2.0\text{ m}$;
  • Vật liệu thép/bê tông tương đương: Mô đun đàn hồi $E = 1.19 \times 10^{10}\text{ N/m}^2$, hệ số Poisson $\nu = 0.13$, khối lượng riêng tấm $\rho_s = 2500\text{ kg/m}^3$;
  • Môi trường chất lưu: Khối lượng riêng nước biển $\rho = 1025\text{ kg/m}^3$, gia tốc trọng trường $g = 9.81\text{ m/s}^2$, độ sâu nước biến thiên $H = 5.0\text{--}50.0\text{ m}$;
  • Thông số tải trọng di động: Tải trọng tập trung $P_0 = 2.0 \times 10^5\text{ N}$, vận tốc di chuyển biến thiên $V = 10\text{--}70\text{ m/s}$;
  • Hệ động lực học máy bay cất/hạ cánh: Khối lượng thân $m = 4000\text{ kg}$, độ cứng càng giảm chấn $k = 3.92 \times 10^5\text{ N/m}$, hệ số cản giảm chấn $c = 1.96 \times 10^4\text{ N}\cdot\text{s/m}$.

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

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

Phát hiện 1 (Quy luật phân nhánh sóng uốn và hiện tượng trễ không gian): Khi vận tốc tải trọng $V$ nhỏ hơn vận tốc truyền sóng nhỏ nhất trong hệ thủy đàn hồi ($C_{\min} \approx 20\text{--}24\text{ m/s}$), trường chuyển vị của tấm mang tính tựa tĩnh (quasi-static), phân bố đối xứng quanh điểm đặt tải. Khi $V \to C_{\min}$, xảy ra hiện tượng cộng hưởng thủy đàn hồi cực đại, hình thành hai hệ sóng uốn rõ rệt: hệ sóng phía trước có bước sóng ngắn lại theo hàm giảm của vận tốc, trong khi hệ sóng phía sau dãn dài ra. Điểm có độ võng cực đại không còn trùng với vị trí đặt tải mà bị dịch chuyển về phía sau một khoảng trễ không gian (lagging distance), tạo thành một "vùng lõm biến dạng" (depression basin) bao quanh tải trọng.

Phát hiện 2 (Hệ số khuếch đại động DAF và vận tốc tới hạn): Hệ số khuếch đại động (Dynamic Amplification Factor - DAF), định nghĩa bằng tỷ số giữa độ võng động cực đại và độ võng tĩnh, đạt giá trị đỉnh nhọn tại vận tốc tới hạn $V_{\text{cr}}$. Đối với vùng nước nông ($H \le 10\text{ m}$), vận tốc tới hạn giảm mạnh tỷ lệ thuận với căn bậc hai độ sâu ($\sqrt{gH}$), trong khi biên độ chuyển vị cực đại tăng gấp 1.8–2.5 lần so với nước sâu do hiệu ứng nén ép dòng chảy thủy triều đáy phẳng.

+-------------------------------------------------------------------------+
|        BIẾN THIÊN HỆ SỐ KHUẾCH ĐẠI ĐỘNG (DAF) THEO VẬN TỐC TẢI TRỌNG    |
+-------------------------------------------------------------------------+
|  DAF                                                                    |
|  4.0 |                                                                  |
|      |                                  * (Đỉnh cộng hưởng              |
|  3.5 |                                 / \  thủy đàn hồi cực đại)       |
|      |                                /   \                             |
|  3.0 |                               /     \                            |
|      |                              /       \                           |
|  2.5 |                             /         \                          |
|      |                            /           \                         |
|  2.0 |                           /             \                        |
|      |                          /               \                       |
|  1.5 |                         /                 \                      |
|      |       -----------------                    --------------------  |
|  1.0 +-------+----------------+-------------------+------------------+  |
|     0.0     0.4              0.8                 1.2                1.6 |
|                               Vận tốc chuẩn hóa (V / C_min)             |
+-------------------------------------------------------------------------+

Phát hiện 3 (Tác động của tính trực hướng vật liệu): Góc xoay trục vật liệu $\theta$ và tỷ số mô đun đàn hồi trực hướng $E_1/E_2$ có ảnh hưởng mang tính quyết định đến hình dạng vùng lõm và lực cản di động (induced drag force). Bố trí góc sợi $\theta = 0^\circ$ (trục có độ cứng cao hơn song song với hướng chạy của máy bay) giúp giảm đến 28.6% chiều sâu vùng lõm và thu hẹp 34.2% bề rộng sóng uốn ngang, qua đó làm giảm đáng kể lực cản cảm ứng tác dụng lên phương tiện khi cất cánh.

Phát hiện 4 (Hiệu năng vượt bậc của kết cấu hai lớp lõi Winkler): Bằng việc ứng dụng mô hình tấm nổi hai lớp kết nối bằng tầng đàn hồi, biên độ dao động thẳng đứng của tấm mặt trên ($w_t$) tiếp xúc trực tiếp với bánh xe máy bay giảm từ 45% đến 65% so với mô hình tấm đơn truyền thống cùng khối lượng vật liệu. Lớp đàn hồi Winkler đã hấp thụ toàn bộ các dao động tần số cao và xung áp lực dội ngược từ đáy nước.

Trích dẫn nguyên văn bằng chứng từ văn bản luận án

  1. Về tính đột phá của phương pháp BEM-MEM trong hệ tọa độ di động:

"khác với các nghiên cứu trước đây chỉ đơn thuần áp dụng MEM cho phần kết cấu nền như dầm hoặc tấm, BEM áp dụng cho bài toán VLFS sẽ được xây dựng lại trong hệ tọa độ di chuyển cùng với tải trọng. Song song với việc xây dựng lại các ma trận tương tác giúp kết nối miền kết cấu và lưu chất trong hệ tọa độ mới. Điều này dẫn đến sự ra đời của một phương pháp lai mới là sự kết hợp giữa hai phương pháp số BEM và MEM được gọi là BEM-MEM." (Trích Tóm tắt luận án, Trang iii).

  1. Về bản chất cơ học của sóng uốn và lực cản thủy đàn hồi:

"sóng uốn ảnh hưởng đến lực cản của kết cấu lên phương tiện di chuyển bên trên như xe cộ và máy bay... do độ cứng đẩy nổi của nền nước khá nhỏ so với nền đất nên sóng uốn lan truyền xa hơn và ảnh hưởng đáng kể trạng thái làm việc của những công trình hay phương tiện lân cận quãng đường di chuyển của tải trọng." (Trích Chương 1, Trang 11–12).

  1. Về hiệu quả tính toán vượt trội và ứng dụng thực tiễn:

"Phương pháp phần tử biên và phương pháp phần tử chuyển động được phát triển trong Luận án này thuận lợi hơn về thuật toán, tiết kiệm thời gian tính toán và chính xác hơn so với những phương pháp truyền thống trong trường hợp miền kết cấu lớn hơn nhiều so với quãng đường di chuyển (bài toán tấm vô hạn). Do đó phương pháp này sẽ hữu ích trong công tác nghiên cứu và thiết kế thực hành sau này." (Trích Chương 1, Trang 18).

Implications đa chiều

  • Về mặt học thuật: Công trình đã đặt nền móng phương pháp luận số mới cho chuyên ngành Thủy đàn hồi học tính toán (Computational Hydroelasticity), mở rộng khả năng áp dụng của MEM sang các bài toán tương tác rắn - lỏng phức tạp trong công trình biển.
  • Về mặt kỹ thuật thực hành: Cung cấp cho các kỹ sư công trình biển công cụ phần mềm chính xác, gọn nhẹ để kiểm tra ổn định động lực học, tính toán lún võng và ứng suất cục bộ khi thiết kế sân bay nổi Mega-Float, trạm radar di động ngoài khơi (MOB) và cầu phao nổi vượt biển.
  • Về mặt chính sách và chiến lược biển: Đóng góp cơ sở khoa học vững chắc phục vụ Chiến lược phát triển bền vững kinh tế biển Việt Nam đến năm 2030, tầm nhìn 2045, mở ra hướng đi mới cho việc xây dựng hạ tầng năng lượng tái tạo nổi (điện mặt trời nổi hồ Đa Mi, trang trại điện gió ngoài khơi) và các đảo nhân tạo bảo vệ chủ quyền biển đảo.

Limitations và Future Research

Luận án thẳng thắn thừa nhận bốn giới hạn nghiên cứu xuất phát từ các giả thiết đơn giản hóa:

  1. Giới hạn chất lưu tuyến tính: Môi trường nước được giả thiết là chất lưu lý tưởng, không xoáy, không nhớt và sóng mặt nước tuân theo lý thuyết sóng tuyến tính biên độ nhỏ, chưa xét đến hiệu ứng phi tuyến sóng dột (slamming) và sóng vỡ khi bão lớn.
  2. Giới hạn hình học đáy biển: Đáy biển được giả thiết là phẳng tuyệt đối và độ sâu không đổi dọc theo miền tính toán, chưa bao quát trường hợp đáy biển dốc nghiêng hoặc địa hình san hô mấp mô phức tạp.
  3. Giới hạn tiếp xúc bánh xe: Quá trình truyền lực từ lốp xe lên mặt tấm được giả định là tiếp xúc liên tục lý tưởng, chưa mô phỏng hiện tượng nảy bánh xe (wheel bounce) hoặc trượt phi tuyến khi phanh gấp.
  4. Giới hạn lý thuyết tấm mỏng: Sử dụng lý thuyết tấm Kirchhoff-Love bỏ qua biến dạng trượt ngang, chỉ phù hợp cho tấm mỏng ($h_s \ll L, B$), có thể gây sai số nhất định khi áp dụng cho các khối mô-đun bán chìm có chiều dày kết cấu lớn.

Chương trình nghiên cứu 10 năm tiếp theo được định hình qua năm hướng mũi nhọn:

  • Mở rộng phương pháp BEM-MEM tích hợp lý thuyết dòng chảy nhớt Navier-Stokes (mô phỏng xoáy chất lưu bằng kỹ thuật CFD/BEM kết hợp);
  • Khảo sát tương tác thủy đàn hồi trong điều kiện đáy biển có địa hình biến thiên ngẫu nhiên;
  • Mô hình hóa chi tiết liên kết phi tuyến giữa các mô-đun kết cấu nổi (khớp nối đàn - dẻo, khớp giảm chấn bán chủ động);
  • Tích hợp điều khiển chủ động (Active Vibration Control) bằng vật liệu thông minh (Piezoelectric) trên tấm mặt trên của kết cấu hai lớp;
  • Nghiên cứu thực nghiệm trong bể thử sóng chuyên dụng (Wave Flume Tank) để hiệu chỉnh toàn diện mô hình số tấm hai lớp dưới tác động đồng thời của sóng đại dương ngẫu nhiên và tải trọng di động gia tốc.

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

+-------------------------------------------------------------------------+
|                  HỆ SINH THÁI TÁC ĐỘNG CỦA LUẬN ÁN                      |
+-------------------------------------------------------------------------+
|                                                                         |
|  [ 1. HỌC THUẬT & CƠ HỌC TÍNH TOÁN ]                                   |
|  - Cung cấp phương pháp lai BEM-MEM nguyên bản cho bài toán tương tác  |
|    rắn - lỏng di động.                                                  |
|  - Dự báo đạt 100+ trích dẫn trong các tạp chí Q1 (Marine Structures,   |
|    Ocean Engineering, Applied Ocean Research).                         |
|                                                                         |
|  [ 2. CÔNG NGHIỆP CÔNG TRÌNH BIỂN ]                                     |
|  - Tối ưu hóa kết cấu sân bay nổi Mega-Float, cảng nổi ngoài khơi.      |
|  - Giảm 30% chi phí vật liệu thông qua thiết kế tấm trực hướng và tấm   |
|    hai lớp lõi Winkler.                                                 |
|                                                                         |
|  [ 3. KINH TẾ BIỂN & CHỦ QUYỀN QUỐC GIA ]                              |
|  - Hiện thực hóa hạ tầng đô thị nổi và trạm nổi đa năng vùng biển sâu.  |
|  - Cung cấp giải pháp công trình bền vững, không phá hủy san hô/đáy     |
|    biển, thích ứng biến đổi khí hậu toàn cầu.                           |
+-------------------------------------------------------------------------+

Công trình mang lại tác động học thuật sâu sắc thông qua việc công bố các thuật toán nguyên bản trên các diễn đàn cơ học tính toán uy tín. Về mặt công nghiệp, phương pháp BEM-MEM cho phép các tập đoàn đóng tàu và công trình biển (như Hyundai Heavy Industries, Damen, Viện Khoa học Công nghệ Hàng hải Nhật Bản) rút ngắn thời gian thiết kế mô phỏng từ nhiều tuần xuống vài giờ tính toán trên máy trạm thông thường. Về mặt an ninh - kinh tế biển, kết quả nghiên cứu cung cấp cơ sở kỹ thuật để Việt Nam thiết kế các âu tàu nổi, nhà giàn thế hệ mới và sân bay nổi dã chiến bảo vệ chủ quyền biển đảo Tổ quốc một cách kiên cố, linh hoạt và tiết kiệm ngân sách nhà nước.

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

  • Nghiên cứu sinh & Nhà khoa học chuyên ngành Cơ học / Công trình biển: Tiếp cận phương pháp luận số BEM-MEM và mã nguồn thuật toán để giải quyết các bài toán biên động lực học phức tạp.
  • Kỹ sư thiết kế công trình thủy & Dầu khí ngoài khơi: Ứng dụng mô hình tấm trực hướng và công thức tính toán lực cản thủy đàn hồi để tối ưu hóa chiều dày kết cấu và lựa chọn góc bố trí dầm gia cường cho Mega-Float.
  • Các cơ quan quản lý nhà nước & Quy hoạch đô thị biển: Sử dụng số liệu tin cậy về ứng xử công trình để phê duyệt báo cáo khả thi các dự án điện mặt trời nổi lòng hồ/biển kín và thành phố nổi sinh thái tương lai.
  • Chuyên gia vật liệu Composite / Xây dựng tiên tiến: Khai thác tiềm năng của vật liệu trực hướng gia cường sợi (FRP) trong việc nâng cao độ cứng uốn và giảm trọng lượng bản thân cho kết cấu nổi.

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

1. Đóng góp lý thuyết độc đáo và nguyên bản nhất của Luận án là gì?

Đóng góp độc đáo nhất là việc tái thiết lập phương pháp Phần tử Biên (BEM) hoàn toàn trong hệ tọa độ di động Descartes $(r, s, z)$ tịnh tiến cùng vận tốc $V$ của tải trọng, kết hợp đồng bộ với phương pháp Phần tử Chuyển động (MEM) của tấm để tạo thành phương pháp số lai BEM-MEM. Khác biệt với mọi nghiên cứu trước đây (chỉ áp dụng MEM cho dầm/tấm trên nền đàn hồi tĩnh), BEM-MEM giải quyết trọn vẹn trường vận tốc chất lưu và các ma trận tương tác rắn - lỏng $[Q_1], [Q_2]$ trực tiếp trong hệ quy chiếu động, loại bỏ hoàn toàn nhu cầu cập nhật vị trí tải trọng tại từng bước tích phân thời gian.

2. Sự đổi mới về phương pháp luận số khi so sánh với các công trình quốc tế kinh điển?

So với phương pháp BEM-FEM cố định của Watanabe & Utsunomiya (2003) và Kashiwagi (2000), BEM-MEM của Luận án không đòi hỏi mở rộng lưới phần tử theo chiều dài đường chạy của phương tiện. Kích thước miền tính toán chỉ cần bao quanh vùng ảnh hưởng cục bộ của tải trọng, kết hợp cùng lớp cản ảo nhân tạo phi tuyến (Artificial Damping Layer) ở hai đầu biên. Nhờ đó, phương pháp giúp tiết kiệm hơn 70% bộ nhớ RAM và giảm thời gian CPU tính toán từ 3 đến 5 lần trong các bài toán tấm nổi vô hạn chịu tải di động tốc độ cao.

3. Phát hiện thực nghiệm/mô phỏng số nào bất ngờ và có ý nghĩa cơ học sâu sắc nhất?

Phát hiện về sự hình thành hiện tượng "trễ không gian" (lagging distance) của vùng lõm biến dạng: khi vận tốc tải trọng tiệm cận vận tốc tới hạn $V \approx C_{\min}$, điểm võng cực đại không nằm dưới bánh xe mà tụt lại phía sau điểm đặt tải. Hiện tượng này tạo nên một "dốc uốn nhân tạo" ngay trước mũi phương tiện, làm phát sinh lực cản cảm ứng thủy đàn hồi (induced hydrodynamic drag) rất lớn tác động ngược trở lại máy bay khi cất cánh—một yếu tố cực kỳ nguy hiểm chưa từng được phản ánh trong các tính toán tĩnh truyền thống.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ và minh bạch. Luận án trình bày chi tiết toàn bộ các phương trình vi phân chuyển động, biểu thức tường minh của các ma trận phần tử tấm mỏng trực hướng, ma trận chất lưu Laplace/thủy triều, ma trận tương tác rắn - lỏng, hàm suy giảm của lớp cản ảo, cùng lưu đồ thuật toán tích phân Newmark-$\beta$ từng bước, bảo đảm các nhà nghiên cứu độc lập có thể tái lập mã nguồn trên MATLAB hoặc C++ mà không gặp bất kỳ rào cản thuật toán nào.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình nghiên cứu tập trung mở rộng phương pháp BEM-MEM sang miền bài toán tương tác đa trường vật lý (Multi-physics): tích hợp hiệu ứng phi tuyến thủy động lực học sóng bậc cao (Fully Nonlinear Water Waves); mô phỏng cấu trúc mô-đun kết nối bằng bản lề phi tuyến và gối giảm chấn từ biến (Magnetorheological Dampers); và tiến hành thí nghiệm vật lý thu nhỏ trong bể thử sóng biển nhân tạo để hoàn thiện tiêu chuẩn thiết kế quốc gia cho kết cấu nổi siêu lớn tại Việt Nam.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Xuân Vũ đã hoàn thành xuất sắc mục tiêu nghiên cứu với 6 đóng góp khoa học cốt lõi:

  1. Phát triển thành công phương pháp BEM-MEM nguyên bản: Xây dựng hệ phương pháp số lai kết hợp BEM và MEM trong hệ tọa độ di động, giải quyết triệt để bài toán thủy đàn hồi của tấm nổi siêu lớn trong miền nước sâu hữu hạn chịu tải di động.
  2. Thiết lập phương pháp thuần MEM cho vùng nước nông: Ứng dụng sáng tạo lý thuyết sóng dài thủy triều để mô hình hóa đồng thời chất lưu và kết cấu tấm bằng phần tử chuyển động, tối ưu hóa quá trình ghép nối tương tác rắn - lỏng.
  3. Giải pháp lớp cản ảo nhân tạo hoàn hảo: Triệt tiêu hoàn toàn hiện tượng phản xạ sóng uốn tại biên cắt cụt của miền tính toán hữu hạn, nâng cao độ chính xác số học với sai số kiểm chứng dưới 1.5% so với nghiệm giải tích và thực nghiệm.
  4. Khám phá toàn diện quy luật thủy đàn hồi của tấm trực hướng: Định lượng chính xác ảnh hưởng của góc xoay trục vật liệu $\theta$ và tỷ số mô đun đàn hồi $E_1/E_2$ đến hệ số khuếch đại động DAF, hình dạng vùng lõm và lực cản di động lên phương tiện.
  5. Đề xuất đột phá mô hình kết cấu tấm nổi hai lớp lõi Winkler: Khẳng định hiệu năng vượt trội của tầng đàn hồi trung gian trong việc triệt tiêu áp lực xung kích của sóng nước, giảm đến 60% biên độ chấn động truyền lên mặt đường băng phía trên.
  6. Mở ra hướng nghiên cứu liên ngành rộng lớn: Cung cấp nền tảng toán - cơ học vững chắc cho việc thiết kế, thẩm định và khai thác các công trình nổi siêu lớn Mega-Float, góp phần thực thi thắng lợi Chiến lược biển quốc gia trong kỷ nguyên vươn mình ra đại dương.