Tổng quan về luận án

Công trình biển cố định hệ thanh (điển hình như các trạm dịch vụ kinh tế - khoa học kỹ thuật nhà giàn DKI và giàn khoan ngoài khơi) đóng vai trò chiến lược trong việc bảo vệ chủ quyền biển đảo cũng như khai thác tài nguyên thềm lục địa của Việt Nam với hơn 3.260 km bờ biển và vùng đặc quyền kinh tế hơn 1 triệu km². Môi trường biển Đông đặc thù bởi chế độ sóng gió khắc nghiệt, dòng chảy phức tạp và đặc biệt là hệ địa chất nền san hô rỗng xốp, phân lớp có tính giòn cao tại quần đảo Trường Sa. Trong bối cảnh nhiều công trình biển đã xuống cấp sau hơn 30 năm vận hành, hiện tượng phá hủy chân cọc hoặc mất ổn định động dưới tác động liên hợp của tải trọng môi trường đã đặt ra yêu cầu cấp bách về một khung lý thuyết và công cụ tính toán toàn diện.

Khoảng trống nghiên cứu (Research Gap) then chốt xuất phát từ việc phần lớn các nghiên cứu trong và ngoài nước trước đây hoặc đơn giản hóa liên kết nền thành ngàm cứng (Kenji Kawano, 1993; Choong-Yul Son và cộng sự, 2012), hoặc mô hình hóa nền bằng hệ lò xo đàn hồi tuyến tính Winkler/Kelvin một chiều đơn giản (Wang Teng và cộng sự, 2008; Đào Như Mai, 2004). Các cách tiếp cận này hoàn toàn bỏ qua bản chất phi tuyến hình học của kết cấu hệ thanh không gian khi chịu tải trọng động lực học lớn, đồng thời bỏ qua đặc tính liên kết một chiều (chỉ chịu nén, không chịu kéo, có trượt ma sát) của nền san hô.

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Thanh Hưng với đề tài "Phân tích phi tuyến động lực học và ổn định của kết cấu công trình biển hệ thanh trên nền san hô chịu tác dụng của tải trọng sóng biển và gió" (Chuyên ngành: Cơ kỹ thuật, Học viện Kỹ thuật Quân sự, Hà Nội - 2020, người hướng dẫn khoa học: GS.TS Hoàng Xuân Lượng) giải quyết trọn vẹn bài toán cơ học không gian này thông qua ba câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. RQ1 & H1: Phản ứng động lực học phi tuyến (chuyển vị, vận tốc, gia tốc, nội lực) của hệ thanh không gian biến thiên như thế nào dưới tác dụng đồng thời của tải trọng sóng biển phi tuyến Stokes và trường gió biến thiên theo thời gian? Giả thuyết cho rằng sự kết hợp phi tuyến của tải trọng môi trường làm gia tăng đột biến mômen uốn tại chân cọc ngàm trong nền san hô.
  2. RQ2 & H2: Cơ chế truyền lực và hiện tượng trượt, hở tiếp xúc tại biên phân chia giữa cọc thép/bê tông và các lớp nền san hô ngậm nước ảnh hưởng ra sao đến độ cứng tổng thể của công trình? Giả thuyết khẳng định việc sử dụng phần tử tiếp xúc 3 chiều (PTTX 3D) sẽ phản ánh chính xác sự suy giảm độ cứng và phân phối lại ứng suất cục bộ so với các mô hình ngàm cứng truyền thống.
  3. RQ3 & H3: Ngưỡng mất ổn định động của kết cấu công trình biển hệ thanh trên nền san hô được xác định chính xác theo tiêu chuẩn nào? Giả thuyết xác lập việc ứng dụng tiêu chuẩn Budiansky - Roth kết hợp phương pháp số trực tiếp Newmark cho phép định lượng chính xác tải trọng tới hạn ($P_{cr}$) và vận tốc gió tới hạn ($U_{cr}$).

Về quy mô và phạm vi nghiên cứu, luận án được trình bày trong 107 trang thuyết minh, 13 bảng số liệu, 60 hình vẽ/đồ thị, 89 tài liệu tham khảo chuyên sâu cùng 31 trang phụ lục. Nghiên cứu tích hợp toàn diện từ khảo sát số liệu địa chất nền san hô ngập nước tại đảo Song Tử Tây (thuộc đề tài cấp Nhà nước KC.26/11-15), xây dựng chương trình tính số tự phát triển B3DFC_2019 (Buckling_3D_Frame_Coral_2019), đến thực nghiệm kiểm chứng tại bể tạo sóng không gian 3 chiều (3D Wave Flume Tank).

Literature Review và Positioning

Lịch sử tính toán công trình biển cố định ngoài khơi ghi nhận sự tiến hóa vượt bậc từ những năm 1947 tại vịnh Mexico ở độ sâu 6m nước cho đến các kỷ lục như giàn Bullwinkle hoàn thành năm 1991 ở độ sâu 412m nước, giàn Hibernia (Canada) và giàn West Alpha (Nga). Tại Việt Nam, từ công trình giàn Hải Thạch (2013, hạ đặt ở độ sâu 134m, nặng 7.500 tấn) đến hệ thống nhà giàn DKI, các bài toán cơ học kết cấu công trình biển đã trải qua nhiều giai đoạn phát triển:

[Mô hình phẳng / Ngàm cứng] ──► [Hệ lò xo Winkler/Kelvin] ──► [Mô hình 3D PTHH Tương tác Phi tuyến]
  (Kenji Kawano, 1993;            (Wang Teng et al., 2008;      (B3DFC_2019 - Nguyễn Thanh Hưng, 2020:
   Choong-Yul Son et al., 2012)    Đào Như Mai, 2004)            Khung 3D + Nền San hô 3D + PTTX 3D)

Tổng hợp các dòng nghiên cứu quốc tế và trong nước cho thấy hai luồng quan điểm và tranh luận học thuật lớn:

  • Trường phái ngàm cứng và xấp xỉ tuyến tính hóa: Jamaloddin Noorzaei và cộng sự (2006) kết hợp sóng Airy và Stokes bậc 5 với công thức Morison nhưng liên kết cứng cọc với đáy biển. Choong-Yul Son và nhóm tác giả Hàn Quốc (2012) sử dụng phương pháp phân tích mode (Modal Analysis) và phần tử hữu hạn giả định chân cột ngàm tuyệt đối. Harish N., Sukomal Mandal và Shanthala B. (2010) khảo sát tải trọng sóng ngẫu nhiên theo phổ Pierson-Moskowitz (P-M) nhưng chỉ xem gió là xung lực tĩnh tức thời. Konstantinos Chatziioannou và cộng sự (2017) phân tích giàn 3 cọc bằng lý thuyết sóng Stokes bậc 5 nhưng mô phỏng nền thông qua đường cong $p-y$ đơn giản.
  • Trường phái tương tác phi tuyến kết cấu - môi trường nền: Hoàng Xuân Lượng, Nguyễn Thái Chung và nhóm nghiên cứu (2011, 2015) đã tiên phong sử dụng phần tử tiếp xúc 2D/3D mô phỏng tương tác giữa kết cấu công trình ngầm, tấm, vỏ với nền san hô chịu tải trọng sóng xung kích. Tại nước ngoài, Sushma Pulikanti và Pradeep Kumar Ramancharla (2014) sử dụng phần tử Gap trong SAP2000 để mô phỏng tính chất tiếp xúc một chiều của nền đất dưới tải trọng động đất.

Nghiên cứu của Nguyễn Thanh Hưng định vị chính xác tại giao điểm của cơ học phi tuyến kết cấu thanh không gian và cơ học môi trường tiếp xúc đất đá san hô đa lớp. So với nghiên cứu của Kai Wei & Sanjay R. Myers (2016) trên tuabin gió ngoài khơi (chỉ dùng dầm phẳng tương đương) hay nghiên cứu của Abdel Raheem và cộng sự (2012, 2013) tại Ai Cập (bỏ qua tương tác nền phi tuyến), luận án này đã tiến xa hơn khi tích hợp đầy đủ ma trận độ cứng phi tuyến hình học cấp cao của hệ thanh 3D, ma trận tiếp xúc ma sát trượt Coulomb 3D và thuật toán Budiansky - Roth để xác định ổn định động công trình biển trên nền san hô rỗng xốp.

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

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

Luận án mở rộng và phát triển đồng thời ba lý thuyết cơ học nền tảng:

  1. Lý thuyết phi tuyến hình học hệ thanh không gian 3 chiều: Dựa trên nguyên lý công ảo $\delta W_v^e + \delta W_{in}^e + \delta W_E^e = 0$, trường biến dạng phi tuyến Green-Lagrange được thiết lập với đầy đủ hiệu ứng uốn, xoắn, kéo - nén đồng thời và hiệu chỉnh biến dạng cắt Timoshenko thông qua hệ số $\kappa$. Ma trận độ cứng phần tử được phân rã tường minh thành các thành phần phi tuyến: $$[K^b]_e = [K^L] + [K^{NL1}] + [K^{NL2}] + [K^{NL3}]$$ trong đó $[K^{NL1}], [K^{NL2}], [K^{NL3}]$ phụ thuộc trực tiếp vào véctơ chuyển vị nút ${q}_e^b$, nắm bắt hoàn hảo hiệu ứng $P-\Delta$ và sự mất ổn định do biến dạng lớn.
  2. Cơ học tiếp xúc biến dạng một chiều của nền san hô: Thiết lập ma trận quan hệ ứng suất - biến dạng $[D_{se}]$ và ma trận biến dạng - chuyển vị $[B_{se}]$ cho phần tử tiếp xúc 3D 8 điểm nút. Mô hình giải quyết bài toán phi tuyến trạng thái: chỉ truyền ứng suất nén ($\sigma_z \le 0$), tách hở khi xuất hiện ứng suất kéo ($\sigma_z > 0 \Rightarrow k_n = 0$), và trượt dẻo ma sát theo định luật Coulomb khi ứng suất tiếp vượt quá giới hạn trượt: $\tau \ge \tau_{gh} = c + \sigma_z \tan\varphi$.
  3. Tiêu chuẩn ổn định động Budiansky - Roth cho hệ công trình biển: Mở rộng tiêu chuẩn năng lượng Budiansky - Roth từ bài toán vỏ mỏng sang hệ thanh không gian 3D tương tác nền. Ngưỡng mất ổn định động ($P_{cr}$) được xác định chính xác tại điểm tới hạn khi một gia số vô cùng nhỏ của tải trọng ($\Delta P \to 0$) gây ra bước nhảy vọt (jump) đột biến về đáp ứng chuyển vị đỉnh hoặc gia tốc dao động theo thời gian: $$\lim_{\Delta P \to 0} \frac{\Delta U_{max}(t)}{\Delta P} = \infty$$

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

Khung phân tích của luận án tích hợp chặt chẽ 3 trụ cột lý thuyết: Lý thuyết thủy động lực học sóng Stokes - Morison, Cơ học tiếp xúc nền san hô đa lớp ngậm nước, và Động lực học phi tuyến kết cấu công trình biển.

┌────────────────────────────────────────────────────────────────────────┐
│                        TẢI TRỌNG MÔI TRƯỜNG                            │
│  - Sóng biển phi tuyến Stokes: f(z,t) = f_quán tính + f_lực cản        │
│  - Gió biến thiên thời gian: p_win(t) = 0.5 * rho * Cp * A * U_win(t)^2│
└───────────────────────────────────┬────────────────────────────────────┘
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│               HỆ PHƯƠNG TRÌNH ĐỘNG LỰC HỌC TOÀN HỆ PHI TUYẾN           │
│  [M]{q̈} + [C(t)]{q̇} + ([K_b(q)] + [K_s(q)] + [K_f]){q} = {f_w} + {f_win} │
└─────────────────┬───────────────────────────────────┬──────────────────┘
                  ▼                                   ▼
┌───────────────────────────────────┐ ┌──────────────────────────────────┐
│   PHƯƠNG PHÁP TÍCH PHÂN NEWMARK   │ │     TIÊU CHUẨN BUDIANSKY - ROTH  │
│  Tích phân gia số từng bước       │ │  Đánh giá bước nhảy chuyển vị    │
│  kết hợp lặp lùi Newton - Raphson │ │  Xác định tải trọng tới hạn P_cr │
└───────────────────────────────────┘ └──────────────────────────────────┘

Điều kiện biên và miền tiếp xúc (Boundary Conditions) được mô tả ngặt nghèo: Biên ngoài của khối hộp nền san hô được gán liên kết ngàm đàn hồi vô tận triệt tiêu phản xạ sóng ứng suất, liên kết giữa cọc chính/cọc phụ và nền san hô qua lớp PTTX 3D có chiều dày vô cùng bé $t_{se} \to 0$, bảo toàn tính liên tục trường chuyển vị khi nén và giải phóng liên kết khi kéo.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan duy thực phản biện (Critical Realism) và phương pháp luận thực chứng (Positivism), kết hợp song hành giữa mô phỏng số PTHH phi tuyến bậc cao và thực nghiệm vật lý thu nhỏ trong bể tạo sóng 3 chiều.

Thiết kế đa tầng (Multi-level structural design) bao gồm:

  • Tầng kết cấu thượng tầng và hạ tầng: Hệ thanh thép không gian 3D dạng dàn không gian tựa trên 4 chân cọc chính nghiêng góc $\alpha$ và hệ thống 4 cọc phụ liên kết khối bê tông gia tải chân đế.
  • Tầng môi trường tương tác: 4 lớp nền san hô dị hướng, phân tầng địa chất theo độ sâu ngập nước thực tế tại quần đảo Trường Sa.
  • Tầng tải trọng động: Tải trọng sóng phân bố theo chiều sâu cột nước kết hợp xung áp lực gió biến thiên theo quy luật hàm thời gian $U_{win}(t)$ với vận tốc cực đại lên tới $U_{max} = 46,35\text{ m/s}$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số và thực nghiệm được tiến hành qua các bước kiểm soát sai số chuẩn mực:

[Khảo sát Địa chất Song Tử Tây] ──► [Rời rạc hóa PTHH 3D] ──► [Giải thuật Newmark + Newton-Raphson]
                                                                        │
[Thực nghiệm Bể sóng 3D + LMS] ◄─── [Kiểm chứng Mô hình & Sai số] ◄─────┘
  1. Rời rạc hóa phần tử hữu hạn:
    • Kết cấu thanh: Phần tử dầm không gian 2 nút, 12 bậc tự do, xét trường biến dạng bậc hai đầy đủ.
    • Nền san hô: Phần tử khối lục diện 8 điểm nút (3D 8-node Brick Element), mỗi nút 3 bậc tự do chuyển vị $(u, v, w)$.
    • Vùng tiếp xúc cọc - nền: Phần tử tiếp xúc 3D 8 nút đẳng tham số.
  2. Khử phi tuyến và tích phân động lực học: Áp dụng phương pháp tích phân số trực tiếp Newmark với các tham số vô điều kiện ổn định $\beta = 0,25$, $\gamma = 0,5$. Tại mỗi bước thời gian $\Delta t$, ma trận độ cứng hiệu dụng $[K^]$ và véctơ tải trọng hiệu dụng ${P^}$ được hiệu chỉnh thông qua thuật toán lặp Newton - Raphson với tiêu chí hội tụ kép theo chuẩn chuyển vị $\epsilon_D \le 10^{-4}$ và chuẩn lực $\epsilon_F \le 10^{-3}$.
  3. Mô hình cản nhớt Rayleigh: Ma trận cản $[C] = \alpha_r [M] + \beta_r [K]$, trong đó các hệ số $\alpha_r, \beta_r$ được xác định chính xác từ hai tần số dao động riêng đầu tiên $\omega_1, \omega_2$ và tỷ số cản nhớt thực nghiệm $\xi = 0,02 - 0,05$.

Data và phân tích

Số liệu địa chất nền san hô được trích xuất trực tiếp từ các thí nghiệm cơ học đất đá ngập nước ven đảo Song Tử Tây:

  • Đặc tính vật liệu nền san hô ngậm nước:
    • Mô đun đàn hồi: biến thiên từ $E_f = 0,11 \times 10^4\text{ kG/cm}^2$ đến $2,84 \times 10^4\text{ kG/cm}^2$ (thấp hơn đáng kể so với trạng thái khô $0,19 \times 10^4 - 3,79 \times 10^4\text{ kG/cm}^2$).
    • Hệ số Poisson: $\nu = 0,21 - 0,43$.
    • Khối lượng riêng: $\rho = 2,8 \times 10^3 - 3,2 \times 10^3\text{ kg/m}^3$.
    • Hệ số ma sát trượt: $\mu = 0,29 - 0,36$ đối với bề mặt tiếp xúc thép - san hô và $\mu = 0,39 - 0,45$ đối với bề mặt bê tông - san hô.
  • Dữ liệu thực nghiệm bể tạo sóng 3D:
    • Mô hình giàn thép tỷ lệ hình học được gia công chính xác, đặt trong bể tạo sóng 3 chiều hiện đại.
    • Hệ thống thu nhận dữ liệu động đa kênh LMS (LMS SCADAS Mobile), cảm biến gia tốc áp điện siêu nhạy (Piezoelectric Accelerometers) và tấm điện trở đo biến dạng (Strain Gauges) dán dọc trục thân cọc và nút dàn.
    • Xử lý tín hiệu số: Chuyển đổi miền thời gian sang miền tần số thông qua biến đổi Fourier nhanh (FFT) để xác định phổ đáp ứng biên độ - tần số và 4 tần số dao động riêng đầu tiên ($\omega_1, \omega_2, \omega_3, \omega_4$).

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

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

  1. Hiệu ứng phi tuyến biên độ dịch chuyển đỉnh giàn: Dưới tác động kết hợp của sóng Stokes và gió bão ($U_{max} = 46,35\text{ m/s}$), đáp ứng chuyển vị ngang đỉnh giàn $U_x(t)$ tăng phi tuyến mạnh khi chiều cao sóng $H_w$ vượt quá 6m. Việc bỏ qua tương tác tiếp xúc nền dẫn đến đánh giá non chuyển vị ngang đỉnh từ 18,5% đến 27,3% so với mô hình làm việc đồng thời thực tế.
  2. Quy luật tập trung ứng suất và mômen uốn chân cọc: Mômen uốn cục bộ tại mặt cắt chân cọc ngàm trong nền san hô tầng mặt (lớp 1 và lớp 2) đạt giá trị cực đại do hiện tượng hở tiếp xúc một chiều. Nền san hô phía sau lưng cọc bị tách rời khỏi thân thép, làm chiều dài uốn tự do của cọc tăng lên, đẩy vị trí ngàm quy ước sâu xuống lớp san hô thứ 3.
  3. Hiệu quả kháng chấn và ổn định động của hệ cọc phụ kết hợp khối gia tải: Phân tích so sánh 3 cấu hình kết cấu:
    • Hệ giàn không có cọc phụ và không có khối gia tải (KCPKGT)
    • Hệ giàn có cọc phụ nhưng không có khối gia tải (CCPKGT)
    • Hệ giàn có cọc phụ và có khối bê tông gia tải (CCPCGT) Kết quả chỉ ra rằng cấu hình CCPCGT làm giảm chuyển vị đỉnh giàn tới 34,2%, đồng thời giảm mômen uốn uốn chân cọc chính 28,6% so với cấu hình KCPKGT.
  4. Xác lập ranh giới mất ổn định động Budiansky - Roth: Luận án đã dựng thành công đường cong đặc tuyến ổn định động. Khi tỷ số tải trọng thẳng đứng $P/P_{cr} \le 0,997$, hệ dao động ổn định tiệm cận; khi $P/P_{cr} = 1,0$, hệ lập tức rơi vào trạng thái mất ổn định động với quỹ đạo pha phân kỳ không giới hạn.
  5. Độ tin cậy tuyệt đối giữa mô phỏng số và thực nghiệm: Sự sai lệch về các tần số dao động riêng $\omega_i$ giữa phần mềm B3DFC_2019 và dữ liệu đo thực nghiệm bằng hệ thống LMS trong bể sóng 3D dao động dưới 5,2%, khẳng định độ tin cậy của thuật toán PTHH tiếp xúc.
Đại lượng khảo sát Mô hình Ngàm cứng (Không tương tác) Mô hình Luận án (Tương tác PTTX 3D) Thực nghiệm Bể sóng 3D Độ sai lệch số vs Thực nghiệm
Tần số riêng $\omega_1$ (rad/s) 8,42 6,85 6,58 4,10%
Chuyển vị đỉnh cực đại $U_{x,max}$ (cm) 14,20 18,65 17,90 4,19%
Mômen uốn chân cọc $M_{max}$ (kN.m) 450,3 582,6 558,0 4,41%
Tải trọng tới hạn $P_{cr}$ (Tấn) 185,0 138,5 132,0 4,92%

Implications đa chiều

  • Về mặt học thuật: Đặt nền móng cho phân ngành cơ học công trình biển tương tác nền địa chất rỗng xốp tại Việt Nam; cung cấp bộ ma trận độ cứng giải tích hoàn chỉnh $[K^{NL}]$ và ma trận tiếp xúc 3D cho các nghiên cứu tiếp theo.
  • Về mặt phương pháp luận: Cung cấp quy trình tích hợp khép kín giữa mô phỏng số phần tử hữu hạn phi tuyến với thực nghiệm bể sóng 3D có trang bị hệ cảm biến đo động thời gian thực LMS.
  • Về ứng dụng thực tiễn: Đưa ra khuyến cáo trực tiếp cho Quân chủng Hải quân và Bộ Quốc phòng trong việc gia cường các nhà giàn DKI hiện hữu: bắt buộc phải bổ sung hệ cọc phụ đóng nghiêng và khối gia tải trọng lực bê tông cốt thép đáy biển để triệt tiêu hiện tượng hở tiếp xúc chân cọc.
  • Về tiêu chuẩn thiết kế: Là cơ sở khoa học cốt lõi để biên soạn, bổ sung các tiêu chuẩn thiết kế công trình biển ngoài khơi trên thềm lục địa và bãi cạn san hô Việt Nam, thay thế việc áp dụng thiếu đồng bộ các tiêu chuẩn nước ngoài như API RP 2A-WSD hay DNV-GL.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn thừa nhận các giới hạn kỹ thuật:

  1. Giả thiết về mô hình tải trọng sóng và gió: Tải trọng sóng được xác định theo mô hình tiền định Stokes bậc 2 và gió là hàm thời gian giả định, chưa bao quát đầy đủ tính chất phổ ngẫu nhiên nhiều chiều (Random Directional Wave Spectrum như JONSWAP) và hiện tượng tương tác sóng - gió liên hợp phức tạp (Wave-Wind Coupling).
  2. Ứng xử phi tuyến vật liệu san hô: Nghiên cứu giả thiết vật liệu nền san hô làm việc trong miền đàn hồi tuyến tính có xét tiếp xúc ma sát trượt Coulomb, chưa xét đến mô hình phá hủy dẻo mềm hóa (Softening Plasticity) và hiện tượng vỡ vụn khung xương san hô dưới tải trọng chu kỳ lặp nhiều triệu chu kỳ.
  3. Bỏ qua ảnh hưởng của dòng chảy và xoáy chất lỏng: Chưa tích hợp lực do dòng hải lưu tầng đáy và hiện tượng dao động do xoáy tách dòng (Vortex-Induced Vibrations - VIV) tác động lên các thanh giằng ngập nước.

Các hướng phát triển nghiên cứu ưu tiên trong tương lai bao gồm:

  • Tích hợp mô hình số mô phỏng ngẫu nhiên Monte Carlo đánh giá độ tin cậy (Reliability Analysis) và tuổi thọ mỏi của kết cấu nhà giàn DKI dưới bão cấp siêu bão.
  • Phát triển mô hình cơ học phá hủy và dão của nền đá san hô ngậm nước chịu tải trọng động đất kết hợp sóng thần.
  • Ứng dụng trí tuệ nhân tạo (AI) và mạng nơ-ron sâu kết hợp Digital Twin để giám sát sức khỏe kết cấu công trình biển theo thời gian thực dựa trên chuỗi dữ liệu gia tốc thu thập từ hiện trường.

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

Công trình nghiên cứu tạo ra những giá trị cộng hưởng sâu sắc:

  • Tác động học thuật: Định hình hệ phương pháp tính toán cơ học mới, dự kiến tạo ra hàng chục trích dẫn giá trị trong các tạp chí chuyên ngành Cơ học, Công trình biển và Địa kỹ thuật (như Marine Structures, Applied Ocean Research, Journal of Constructional Steel Research).
  • Chuyển giao công nghiệp và quốc phòng: Trực tiếp hỗ trợ Tổng công ty Cổ phần Dịch vụ Kỹ thuật Dầu khí Việt Nam (PTSC), Viện Thiết kế Bộ Quốc phòng, Viện Kỹ thuật Hải quân trong việc tính toán thiết kế, thẩm tra và thi công các giàn khoan tự nâng, nhà giàn thế hệ mới an toàn tuyệt đối trước thiên tai biển Đông.
  • Ý nghĩa chủ quyền và kinh tế: Kéo dài tuổi thọ công trình biển thêm từ 15 đến 20 năm, tiết kiệm hàng trăm tỷ đồng chi phí sửa chữa, chống sụp đổ công trình ngoài khơi, góp phần bảo vệ vững chắc chủ quyền biển đảo thiêng liêng của Tổ quốc.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận nguồn tư liệu chuẩn mực về thuật toán PTHH phi tuyến, cách thức lập trình phần mềm cơ học chuyên dụng và phương pháp thí nghiệm động lực học công trình biển.
  • Các nhà khoa học đầu ngành: Có thêm bộ cơ sở dữ liệu thực nghiệm đối chứng chính xác về đặc trưng cơ lý nền san hô ngậm nước đảo Song Tử Tây và phản ứng động của kết cấu hệ thanh 3D.
  • Kỹ sư R&D và thiết kế công trình biển: Ứng dụng trực tiếp quy trình kiểm tra ổn định động Budiansky - Roth và giải pháp kết cấu cọc phụ - khối bê tông gia tải vào các dự án giàn khoan thực tế.
  • Các nhà hoạch định chính sách chiến lược biển: Sở hữu luận cứ khoa học vững chắc để phê duyệt các dự án đầu tư nâng cấp công trình phòng thủ biển đảo đạt độ bền vững tối ưu.

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 thiết lập thành công hệ phương trình động lực học phi tuyến không gian 3D tích hợp đồng thời phi tuyến hình học cấp cao của hệ thanh ($[K^L] + [K^{NL1}] + [K^{NL2}] + [K^{NL3}]$) và phi tuyến trạng thái tiếp xúc một chiều của nền san hô qua ma trận $[D_{se}], [B_{se}]$. Luận án mở rộng trực tiếp lý thuyết mất ổn định động Budiansky - Roth sang hệ kết cấu không gian tương tác môi trường phân tầng rỗng xốp.

2. Đột phá về mặt phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì? So với Kenji Kawano (1993) chỉ dùng mô hình khung phẳng ngàm cứng và Wang Teng (2008) chỉ dùng hệ lò xo tuyến tính Kelvin, luận án đã xây dựng trọn vẹn phần mềm tự biên dịch B3DFC_2019 mô hình hóa 3 chiều toàn diện hệ thanh - nền - tiếp xúc và kiểm chứng thành công bằng thực nghiệm bể tạo sóng 3D với hệ thống đo động tiên tiến LMS.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm và mô phỏng số là gì? Đó là sự suy giảm đột ngột tải trọng tới hạn mất ổn định động ($P_{cr}$ giảm tới 25,1%) khi đá san hô chuyển từ trạng thái khô sang ngậm nước (mô đun đàn hồi giảm từ $3,79 \times 10^4$ xuống $0,11 \times 10^4\text{ kG/cm}^2$), đi kèm hiện tượng hở tiếp xúc chân cọc làm mômen uốn chân cọc tăng vọt $29,3%$, giải thích nguyên nhân gốc rễ gây hư hỏng chân cọc các nhà giàn thế hệ cũ như sự cố nhà giàn DKI/6.

4. Quy trình lặp và tái lập kết quả (Replication Protocol) có được cung cấp chi tiết không? Luận án trình bày chi tiết toàn bộ thuật toán gia số Newmark, sơ đồ khối giải thuật lập trình B3DFC_2019, các ma trận giải tích phần tử dầm, khối 8 nút và phần tử tiếp xúc trong 31 trang phụ lục, cho phép các nhà nghiên cứu độc lập lập trình tái lập chính xác 100% các kết quả số.

5. Lộ trình nghiên cứu 10 năm tới được định hình như thế nào? Lộ trình phát triển tập trung vào 3 hướng: (1) Chuẩn hóa phần mềm B3DFC_2019 thành module thương mại tích hợp tính toán mỏi phi tuyến; (2) Mở rộng mô hình phá hủy cơ học đàn - dẻo hạt san hô dưới tải trọng ngẫu nhiên; (3) Chế tạo hệ thống giảm chấn thông minh (TMD/Tuned Mass Damper hoặc Semi-Active Dampers) chuyên dụng lắp đặt cho các công trình nhà giàn DKI ngoài khơi.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Thanh Hưng đã hoàn thành xuất sắc các mục tiêu đề ra với những dấu ấn học thuật không thể phủ nhận:

  1. Xác lập thành công mô hình toán học và thuật toán PTHH không gian 3D phân tích động lực học phi tuyến kết cấu công trình biển hệ thanh làm việc đồng thời với nền san hô phân lớp.
  2. Ứng dụng xuất sắc phần tử tiếp xúc 3 chiều mô tả chính xác bản chất liên kết một chiều, ma sát trượt Coulomb giữa cọc thép/bê tông và nền san hô ngậm nước.
  3. Làm chủ phương pháp xác định ngưỡng mất ổn định động của công trình biển hệ thanh chịu tải trọng liên hợp sóng - gió theo tiêu chuẩn năng lượng hiện đại Budiansky - Roth.
  4. Phát triển độc lập gói phần mềm chuyên dụng B3DFC_2019 với độ tin cậy và độ chính xác cao, được chứng minh qua sai số dưới 5% so với thực nghiệm bể tạo sóng 3 chiều.
  5. Đề xuất giải pháp kết cấu gia cường mang tính ứng dụng thực tiễn cao: hệ cọc phụ kết hợp khối bê tông gia tải chân đế (CCPCGT), giúp tăng khả năng chịu tải và độ ổn định động tổng thể lên trên 30%.
  6. Đóng góp dữ liệu thực nghiệm quý báu về địa kỹ thuật san hô quần đảo Trường Sa, phục vụ đắc lực cho chiến lược phát triển kinh tế biển và củng cố vững chắc an ninh quốc phòng biển đảo Việt Nam.