Tổng quan về luận án

Công trình biển cố định hệ thanh ngoài khơi, tiêu biểu là hệ thống nhà giàn DKI và các giàn khoan khai thác dầu khí tại thềm lục địa và quần đảo Trường Sa của Việt Nam, giữ vai trò chiến lược sống còn đối với chủ quyền an ninh biển đảo và phát triển kinh tế biển. Đặt trong môi trường hải dương nhiệt đới khắc nghiệt, các kết cấu này chịu tác động đồng thời của tổ hợp tải trọng thủy động lực sóng biển, bão gió cực hạn và tải trọng sàn công tác lớn. Đặc biệt, nền móng của công trình được đặt trực tiếp trên tầng địa chất đá san hô ngập nước—một loại vật liệu địa kỹ thuật sinh hóa có tính dị hướng, độ rỗng cao, ứng xử giòn và khả năng chịu kéo gần như bằng không. Khi các công trình tiến vào giai đoạn khai thác lâu năm, hiện tượng suy giảm độ cứng và nguy cơ mất ổn định động lực học phi tuyến dưới tác động sóng gió cực hạn trở thành một hiểm họa nghiêm trọng dẫn đến sụp đổ công trình.

Khoảng trống nghiên cứu then chốt (research gap) mà luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (mã số 9.01) của tác giả Nguyễn Thanh Hưng, dưới sự hướng dẫn khoa học của GS.TS Hoàng Xuân Lượng tại Học viện Kỹ thuật Quân sự (2020), định vị giải quyết là: Sự thiếu vắng một mô hình tính toán giải tích - số không gian 3 chiều (3D) đồng bộ và tường minh có khả năng mô tả chính xác tương tác động lực học phi tuyến giữa kết cấu hệ thanh và nền san hô phân lớp ngậm nước, kết hợp tiêu chuẩn đánh giá mất ổn định động lực học công trình biển dưới tải trọng tổ hợp sóng phi tuyến và gió biến thiên theo thời gian.

Để giải quyết triệt để khoảng trống này, luận án tập trung vào ba câu hỏi nghiên cứu cốt lõi:

  1. Câu hỏi 1: Mô hình toán học nào cho phép mô tả chính xác ứng xử phi tuyến hình học của thanh không gian 3D kết hợp với ứng xử tiếp xúc một chiều (unilateral contact) và ma sát trượt Coulomb tại mặt phân giới cọc - nền san hô?
  2. Câu hỏi 2: Tiêu chuẩn ổn định động nào phản ánh chân thực trạng thái tới hạn sụp đổ của kết cấu công trình biển hệ thanh khi chịu đồng thời tải trọng sóng Stokes bậc 2 và gió động $U_{max} = 46,35\text{ m/s}$?
  3. Câu hỏi 3: Mức độ ảnh hưởng định lượng của cọc phụ (cọc váy), khối bê tông gia tải đỉnh giàn và sự suy giảm môđun đàn hồi nền san hô ngậm nước ($E = 0,11 \cdot 10^4\text{ kG/cm}^2 - 2,84 \cdot 10^4\text{ kG/cm}^2$) đến tần số dao động riêng và biên độ đáp ứng chuyển vị đỉnh giàn là bao nhiêu?

Hai giả thuyết khoa học chính được đặt ra:

  • Giả thuyết $H_1$: Ứng xử phi tuyến tiếp xúc một chiều tại giao diện cọc - nền san hô làm mềm hóa hệ kết cấu, dẫn đến việc giảm tần số dao động riêng từ 8% - 15% và hạ thấp đáng kể ngưỡng tải trọng mất ổn định động lực học so với các mô hình ngàm cứng hoặc mô hình lò xo Winkler tuyến tính truyền thống.
  • Giả thuyết $H_2$: Bố trí hệ cọc phụ xiên kết hợp tối ưu hóa khối lượng gia tải sàn đỉnh sẽ làm tăng mômen kháng uốn tại mặt cắt chân cọc chính từ 20% - 35%, đồng thời triệt tiêu hiện tượng tách hở tiếp xúc và ngăn chặn phân kỳ động lực học theo tiêu chuẩn Budiansky - Roth.

Khung lý thuyết của luận án tích hợp phương pháp phần tử hữu hạn phi tuyến (Nonlinear FEM), lý thuyết sóng phi tuyến Stokes, phương trình Morison hiệu chỉnh, lý thuyết tiếp xúc biến dạng trượt Coulomb và tiêu chuẩn ổn định động Budiansky - Roth. Luận án nghiên cứu mô hình không gian đại diện cho nhà giàn DKI tại bãi cạn san hô ngập nước đảo Song Tử Tây (thuộc đề tài cấp Nhà nước KC.26/11-15), kết hợp thực nghiệm kiểm chứng quy mô trên mô hình vật lý trong bể tạo sóng 3 chiều hiện đại.

Literature Review và Positioning

Nghiên cứu ứng xử động lực học và độ ổn định của công trình biển cố định hệ thanh đã trải qua nhiều giai đoạn phát triển với những trường phái tiếp cận khác nhau trong y văn quốc tế và trong nước:

Trường phái thứ nhất tập trung vào việc mô hình hóa kết cấu hệ thanh ngàm cứng đáy biển hoặc thay thế nền bằng liên kết lò xo tuyến tính đơn giản. Kenji Kawano [61] cùng Po-Yen Chang và các cộng sự [68] đã phát triển các mô hình khung phẳng để phân tích đáp ứng ngẫu nhiên của giàn khoan dưới tác dụng của phổ sóng biển, nhưng bỏ qua hoàn toàn tính tương tác không gian 3 chiều và biến dạng phi tuyến của nền móng. Jamaloddin Noorzaei và cộng sự [54] đã ứng dụng lý thuyết sóng Stokes bậc 5 và phương trình Morison để tính toán lực sóng tác dụng lên hệ thanh ống thép trụ tròn nhưng vẫn giả định liên kết ngàm tuyệt đối với đáy biển. Wang Teng, Huajun Li và Kuithua Wang [86] đã cải tiến bài toán bằng cách thay thế nền san hô bằng hệ lò xo Kelvin kết hợp tuyến tính hóa thống kê, tuy nhiên mô hình này không thể hiện được tính bất đối xứng trong khả năng chịu lực của nền san hô (chịu nén nhưng không chịu kéo).

Trường phái thứ hai mở rộng phân tích tổ hợp tải trọng sóng - gió và động đất. Choong-Yul Son và cộng sự [44] đã kết hợp phương pháp giải tích và PTHH 3D phân tích mode dao động dưới tác dụng đồng thời của sóng và gió. Harish N., Sukomal Mandal và Shanthala B. [51] đã phân tích động lực học công trình biển hệ thanh sử dụng phổ sóng Pierson-Moskowitz (P-M) kết hợp tải trọng xung gió tức thời. Abdel Raheem, Elsayed M. Abdel Aal, Aly G. Abdel Shafy và Fayez K. Abdel Seed [75] đã sử dụng lý thuyết sóng Stokes bậc 5 và FEM 3D khảo sát nội lực cọc chính nhưng chưa xét đến tính phi tuyến hình học của kết cấu thượng tầng và sự trượt tróc tại biên cọc - nền. Konstantinos Chatziioannou, Vanessa Katsardi, Apostolos Koukouselis và Euripidis Mistakidis [62] đã phân tích giàn jacket 3 chân với đường cong tương tác $p-y$, nhận thấy chuyển vị chân cột phụ thuộc mạnh vào chiều sâu ngàm và độ cứng phi tuyến của đất nền.

Tại Việt Nam, các nghiên cứu tiên phong của Nguyễn Đông Anh [2], [31], Hoàng Xuân Lượng và Nguyễn Thái Chung [5], [14], [18], Nguyễn Tiến Khiêm [11], Phạm Khắc Hùng [10], Đào Như Mai [22], Lê Anh Tuấn [26] và Lê Hoàng Anh [1] đã từng bước tiếp cận bài toán tương tác công trình biển - nền san hô. Tuy nhiên, các công bố trước đây chủ yếu dừng lại ở bài toán tĩnh, bài toán phẳng 2D, hoặc chỉ khảo sát đáp ứng dao động cưỡng bức tuyến tính dưới tải trọng xung nổ hoặc tải trọng điều hòa đơn giản.

Luận án của NCS. Nguyễn Thanh Hưng tạo ra bước đột phá định vị khoa học bằng cách tích hợp toàn diện 4 yếu tố then chốt chưa từng được giải quyết đồng thời trong bất kỳ công trình nào:

  1. Mô hình hóa kết cấu thanh không gian 3D với 6 bậc tự do/nút, xét đến ma trận độ cứng phi tuyến hình học đầy đủ $[K_{NL1}]$, $[K_{NL2}]$, $[K_{NL3}]$;
  2. Mô hình hóa môi trường nền san hô phân lớp ngậm nước 3 chiều bằng phần tử lục diện 8 nút đẳng tham số;
  3. Thiết lập phần tử tiếp xúc 3D (PTTX 3D) mô tả chính xác liên kết một chiều (chỉ chịu nén, không chịu kéo) và cơ chế trượt giới hạn Coulomb giữa vỏ cọc thép và nền san hô;
  4. Ứng dụng tiêu chuẩn Budiansky - Roth để xác định chính xác ngưỡng mất ổn định động lực học phi tuyến của công trình dưới tác dụng đồng thời của tải trọng sóng Stokes bậc 2 và trường gió động biến thiên theo thời gian $U_{win}(t)$.

So với hai nghiên cứu quốc tế điển hình là Abdel Raheem et al. [75] (chỉ phân tích đáp ứng động lực học tuyến tính của khung jacket liên kết cứng) và Nam-Il Kim & Dong-Ho Choi [66] (chỉ khảo sát ổn định tĩnh của khung thép phẳng), luận án đã vượt trội khi thiết lập một hệ thống giải thuật PTHH không gian hoàn chỉnh, giải quyết trọn vẹn cả bài toán dao động phi tuyến lẫn bài toán ổn định động có xét đến cơ chế trượt tróc nền móng thực tế.

Đó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à làm phong phú sâu sắc lý thuyết động lực học công trình biển và cơ học tiếp xúc thông qua 4 đóng góp lý thuyết cốt lõi:

  • Mở rộng lý thuyết ổn định động lực học phi tuyến cho hệ thanh không gian tương tác nền: Luận án đã phát triển thành công tiêu chuẩn Budiansky - Roth vốn ban đầu dùng cho kết cấu vỏ mỏng và vòm đàn hồi sang bài toán kết cấu không gian phức tạp tương tác môi trường phân lớp. Bằng cách khảo sát bước nhảy đột biến của đáp ứng chuyển vị ngang cực đại đỉnh giàn $U_x$ theo thông số gia tải không thứ nguyên $P/P_{cr}$, luận án đã định lượng hóa ranh giới phân kỳ ổn định động lực học khi tải trọng đạt ngưỡng $P/P_{cr} = 1,0$.
  • Bổ sung toán học cho ma trận độ cứng phần tử thanh không gian 3D phi tuyến: Dựa trên trường chuyển vị đầy đủ xét đến góc xoay lớn và biến dạng cắt hiệu chỉnh qua hệ số $\alpha_s$, luận án đã thiết lập tường minh các ma trận độ cứng phi tuyến cấp 1, cấp 2 và cấp 3 ($[K_{NL1}], [K_{NL2}], [K_{NL3}]$) xuất phát từ nguyên lý công ảo biến phân: $$\delta W_v^e = -\delta {q_e^b}^T \left( [K_b^L] + [K_{NL1}] + [K_{NL2}] + [K_{NL3}] \right) {q_e^b}$$
  • Thiết lập lý thuyết phần tử tiếp xúc 3D phi tuyến cho nền san hô: Xây dựng mô hình ứng xử phi tuyến cho lớp tiếp xúc giữa cọc và nền san hô dựa trên ma trận độ cứng tiếp xúc $[K_s^e]$ và ma trận đàn dẻo $[D_{se}]$. Mô hình phản ánh đúng bản chất vật lý: khi biến dạng pháp tuyến $\varepsilon_z > 0$ (kéo), độ cứng pháp tuyến $k_{nz} \to 0$ (tách hở tiếp xúc); khi $\varepsilon_z \le 0$ (nén), độ cứng pháp tuyến đạt giá trị chống nén $k_n$; khi ứng suất tiếp vượt qua ứng suất trượt giới hạn $\tau_{gh}$, hiện tượng trượt Coulomb xảy ra với độ cứng tiếp tuyến suy giảm $k_{res}$.
  • Xác lập hệ phương trình vi phân dao động phi tuyến toàn phần của hệ tương tác: $$[M]{\ddot{q}} + [C]{\dot{q}} + \left( [K_L] + [K_{NL}({q})] + [K_{se}({q})] \right){q} = {f_{wave}(t)} + {f_{wind}(t)}$$ trong đó $[C]$ là ma trận cản nhớt Rayleigh phụ thuộc thời gian $[C] = \alpha_r [M] + \beta_r [K]$, phản ánh tiêu tán năng lượng của cả kết cấu thép lẫn nền san hô ngậm nước.

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

Khung phân tích của luận án tích hợp đa tầng ba lý thuyết cơ bản: lý thuyết thủy động lực học sóng phi tuyến (Stokes bậc 2 kết hợp công thức Morison), lý thuyết khí động lực học công trình và lý thuyết cơ học môi trường liên tục tiếp xúc.

Sự độc đáo của khung phân tích nằm ở giải thuật lặp phi tuyến kép:

  1. Vòng lặp trong (Không gian tiếp xúc): Tại mỗi bước tích phân thời gian $\Delta t$, sử dụng phương pháp lặp biến dạng gia tăng để kiểm tra trạng thái tiếp xúc (dính kết - trượt dẻo - tách hở) của từng phần tử tiếp xúc 3D xung quanh cọc chính và cọc phụ;
  2. Vòng lặp ngoài (Phi tuyến hình học và hội tụ động học): Ứng dụng phương pháp tích phân từng bước Newmark kết hợp thuật toán lặp sai phân tự do Newton - Raphson cải biên với hai điều kiện kiểm soát dung sai nghiêm ngặt: sai số chuyển vị $\epsilon_D \le 10^{-4}$ và sai số lực nút mất cân bằng $\epsilon_F \le 10^{-3}$.

Điều kiện biên không gian được xử lý chính xác thông qua mô hình miền nghiên cứu hình hộp của nền san hô, kích thước biên ngoài được tối ưu hóa qua giải thuật lặp sao cho ảnh hưởng sóng phản xạ tại biên tự do triệt tiêu hoàn toàn và trường ứng suất biên tiến về trạng thái tự nhiên.

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ủ triết lý nghiên cứu thực chứng nghiêm ngặt (positivism), kết hợp chặt chẽ giữa mô hình hóa toán học giải tích, tính toán mô phỏng số PTHH nâng cao và kiểm chứng thực nghiệm vật lý tỷ lệ trên mô hình bể tạo sóng không gian 3 chiều.

Thiết kế nghiên cứu được chia thành 3 phân hệ liên hoàn:

  • Phân hệ 1 (Mô hình hóa số học): Rời rạc hóa kết cấu nhà giàn không gian bằng phần tử thanh dầm 3D (2 nút, 12 bậc tự do); rời rạc hóa nền san hô 3 lớp bằng phần tử khối lục diện 8 nút (mỗi nút 3 bậc tự do chuyển vị); mô hình hóa bề mặt cọc - đất bằng phần tử tiếp xúc 3D không độ dày.
  • Phân hệ 2 (Phát triển phần mềm chuyên dụng): Lập trình hoàn chỉnh bộ chương trình tính toán số mang tên B3DFC_2019 (Buckling_3D_Frame_Coral_2019) bằng ngôn ngữ lập trình khoa học hiệu năng cao, tích hợp thuật toán giải phi tuyến Newmark - Newton Raphson.
  • Phân hệ 3 (Thực nghiệm kiểm chứng): Chế tạo mô hình giàn thép tỷ lệ hình học đặt trong bể tạo sóng 3D tại phòng thí nghiệm chuyên ngành, gắn hệ thống cảm biến gia tốc và biến dạng động điện trở kết nối với bộ thu thập và phân tích dữ liệu đa kênh LMS.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm các giao thức đo đạc và kiểm định có độ tin cậy khoa học cao:

  • Xác lập thông số đầu vào địa kỹ thuật: Sử dụng số liệu thí nghiệm cơ lý nền san hô ngập nước tại bãi cạn đảo Song Tử Tây (thuộc đề tài KC.26/11-15):
    • Lớp 1 (San hô bở rời vụn vụn gắn kết yếu): Chiều dày $h_1 = 3,0\text{ m}$, $E_1 = 0,11 \cdot 10^4\text{ kG/cm}^2$, hệ số Poisson $\nu_1 = 0,38$, khối lượng riêng $\rho_1 = 2,8 \cdot 10^3\text{ kg/m}^3$.
    • Lớp 2 (Đá san hô gắn kết trung bình): Chiều dày $h_2 = 5,0\text{ m}$, $E_2 = 1,45 \cdot 10^4\text{ kG/cm}^2$, $\nu_2 = 0,32$, $\rho_2 = 2,95 \cdot 10^3\text{ kg/m}^3$.
    • Lớp 3 (Đá san hô phong hóa nhẹ đến nguyên khối): Chiều dày $h_3 = 12,0\text{ m}$, $E_3 = 2,84 \cdot 10^4\text{ kG/cm}^2$, $\nu_3 = 0,25$, $\rho_3 = 3,1 \cdot 10^3\text{ kg/m}^3$.
    • Hệ số ma sát tiếp xúc: Giữa cọc thép và nền san hô $\mu = 0,32$; ứng suất trượt tới hạn $\tau_{gh} = 1,85\text{ kG/cm}^2$.
  • Xác lập tải trọng thủy khí động học:
    • Tải trọng sóng: Mô hình sóng phi tuyến Stokes bậc 2 với dải chiều cao sóng khảo sát $H_w = 3\text{ m} \div 11\text{ m}$, chu kỳ sóng $T_w = 7,5\text{ s} \div 12,0\text{ s}$, tính toán qua phương trình Morison: $$f(z,t) = \frac{1}{2} C_D \rho_{water} D |u_x| u_x + C_M \rho_{water} \frac{\pi D^2}{4} a_x$$ với hệ số cản $C_D = 1,2$, hệ số quán tính $C_M = 2,0$.
    • Tải trọng gió: Vận tốc gió cực hạn dạng xung biến thiên thời gian với giá trị đỉnh $U_{max} = 46,35\text{ m/s}$ (tương đương bão cấp 14-15), xác định theo công thức: $$p_{win}(t) = \frac{1}{2} \rho_{air} C_p A U_{win}^2(t) \cos^2\beta$$ với $\rho_{air} = 1,225\text{ kg/m}^3$, hệ số áp lực gió $C_p = 1,3$.
  • Quy trình thực nghiệm bể sóng 3D:
    • Cảm biến gia tốc 3 chiều áp điện (Piezoelectric Accelerometers) với độ nhạy $100\text{ mV/g}$, dải đo tần số $0,5 - 3000\text{ Hz}$.
    • Cảm biến đo biến dạng vi phân (Strain Gauges dạng điện trở lá $120,\Omega$) dán tại các mặt cắt chân cọc chính và cọc phụ để đo mômen uốn động.
    • Hệ thống máy đo động đa kênh LMS Scadas Mobile kết hợp phần mềm LMS Test.Lab xử lý tín hiệu FFT phân tích phổ biên độ - tần số và mode shape.

Data và phân tích

Luận án đã thực hiện kiểm tra độ tin cậy và chuẩn hóa chương trình tính B3DFC_2019 thông qua hai bước so sánh độc lập:

  1. Kiểm tra tính đúng đắn với lời giải giải tích và bài toán chuẩn quốc tế: So sánh kết quả 4 tần số dao động riêng đầu tiên với mô hình kiểm chuẩn của Kenji Kawano [61] và SAP2000 cho thấy sai số tương đối cực đại $< 1,85%$.
  2. So sánh đối chứng với kết quả thực nghiệm bể tạo sóng 3D:
    • Tần số dao động riêng thứ nhất theo tính toán số của B3DFC_2019 là $f_1 = 1,42\text{ Hz}$, kết quả đo thực nghiệm trên hệ thống LMS là $f_{1,exp} = 1,38\text{ Hz}$ (sai số chỉ đạt $2,89%$).
    • Biên độ gia tốc cực đại tại đỉnh giàn khi có khối gia tải: kết quả mô phỏng đạt $\ddot{U}{x,max} = 0,382\text{ m/s}^2$, kết quả đo thực nghiệm đạt $\ddot{U}{x,exp} = 0,395\text{ m/s}^2$ (độ trùng khớp đạt $96,7%$).

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

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

Nghiên cứu mang lại 5 phát hiện khoa học đột phá có giá trị thực tiễn và lý thuyết sâu sắc:

  1. Quy luật suy giảm tần số riêng do tính chất tiếp xúc phi tuyến của nền san hô: Khi xét đến tương tác tiếp xúc một chiều (PTTX 3D) thay vì ngàm cứng vô hạn, tần số dao động riêng bậc 1 của công trình nhà giàn giảm từ $1,78\text{ Hz}$ xuống $1,42\text{ Hz}$ (giảm $20,22%$). Điều này chứng minh rằng việc giả định ngàm cứng truyền thống đã đánh giá quá cao độ cứng của hệ, dẫn đến nguy cơ công trình rơi vào miền cộng hưởng nguy hiểm với tần số sóng biển ($0,1\text{ Hz} - 0,25\text{ Hz}$) trong các kịch bản biển động mạnh.
  2. Phát hiện về cơ chế truyền lực và phân phối mômen uốn tại chân cọc: Đối với hệ kết cấu có cọc phụ (CCPCGT và CCPKGT), cọc phụ xiên đóng vai trò như một thanh chống mềm gánh tới $34,6%$ tổng mômen uốn chân công trình, làm giảm mômen uốn chân cọc chính từ $245,6\text{ kNm}$ xuống còn $160,8\text{ kNm}$ ($34,53%$). Điều này giải thích nguyên nhân hư hỏng gãy nứt chân cọc tại các nhà giàn đời cũ không có cọc phụ (như sự cố chân cọc nhà giàn DKI/6).
  3. Hiệu ứng phi tuyến của khối bê tông gia tải đỉnh giàn: Việc bổ sung khối gia tải sàn công tác làm giảm chuyển vị ngang cực đại đỉnh giàn $U_x$ từ $18,4\text{ cm}$ xuống còn $11,2\text{ cm}$ (giảm $39,13%$) do tăng quán tính chống lật, tuy nhiên lại làm tăng mômen uốn uốn chân cọc lên $14,8%$ và làm giảm tần số dao động cơ bản của hệ.
  4. Xác định ngưỡng mất ổn định động lực học Budiansky - Roth: Dưới tác dụng của tổ hợp tải trọng sóng Stokes bậc 2 ($H_w = 9,5\text{ m}$) và gió động cực hạn ($U_{max} = 46,35\text{ m/s}$), kết cấu duy trì trạng thái ổn định đàn hồi ở mức tải $P/P_{cr} = 0,997$. Khi tải trọng vượt qua $P/P_{cr} = 1,0$, đáp ứng chuyển vị ngang và góc xoay đỉnh giàn tăng vọt phi tuyến dạng phân kỳ bùng nổ (dynamic jump), minh chứng cho sự mất ổn định động lực học toàn thể do phá hủy dẻo tiếp xúc đồng thời tại lớp san hô trên cùng.
  5. Độ nhạy của ổn định động với môđun đàn hồi nền ngậm nước: Khi môđun đàn hồi $E_3$ của lớp san hô chịu lực chính giảm $50%$ (do hiện tượng ngậm nước bão hòa hoặc xâm thực hóa học), tải trọng tới hạn mất ổn định $P_{cr}$ của công trình suy giảm tới $28,4%$, và chiều cao sóng tới hạn gây phá hủy giảm từ $H_{w,cr} = 10,8\text{ m}$ xuống còn $8,1\text{ m}$.

Implications đa chiều

  • Về mặt lý thuyết: Bác bỏ giả thiết ngàm cứng đáy biển trong tính toán công trình biển trên nền địa chất san hô; khẳng định tính bắt buộc của việc xét ma trận độ cứng phi tuyến tiếp xúc $[K_{se}]$ và phi tuyến hình học $[K_{NL}]$ trong phân tích dao động và ổn định công trình biển.
  • Về phương pháp luận: Cung cấp giải thuật số tích hợp toàn diện giữa phương pháp Newmark cải tiến và tiêu chuẩn Budiansky - Roth, có thể mở rộng áp dụng cho các dạng kết cấu ngoài khơi khác như chân đế điện gió ngoài khơi (Offshore Wind Turbines), giàn TLP, giàn nửa chìm nửa nổi.
  • Về ứng dụng thực tiễn kỹ thuật: Đưa ra các khuyến cáo định lượng cụ thể phục vụ trực tiếp cho Bộ Quốc phòng và ngành dầu khí:
    1. Tuyệt đối không thiết kế chân cọc cắm đơn lập không có cọc phụ trên nền san hô ngậm nước;
    2. Chiều sâu ngàm cọc chính tối thiểu phải xuyên sâu qua lớp san hô mềm số 1 và số 2, cắm sâu vào lớp số 3 ít nhất $6,0\text{ m}$;
    3. Bắt buộc phải bố trí các vòng jacket gia cường liên kết giữa cọc chính và cọc phụ tại vùng nước dao động triều để triệt tiêu tập trung ứng suất.
  • Về chính sách và tiêu chuẩn xây dựng: Đóng góp cơ sở khoa học trực tiếp để hoàn thiện quy chuẩn thiết kế công trình quốc phòng ven biển và hải đảo Việt Nam trong bối cảnh biến đổi khí hậu gia tăng tần suất siêu bão.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính khách quan:

  1. Giới hạn cơ học vật liệu nền: Mô hình ứng xử của đá san hô vẫn giả định là đàn hồi tuyến tính chịu nén đơn thuần kết hợp mặt tiếp xúc Coulomb, chưa mô phỏng được quá trình dập vỡ vi cấu trúc (grain crushing) và tích lũy biến dạng dẻo trượt thể tích dưới tải trọng trùng phục hàng triệu chu kỳ sóng.
  2. Giới hạn tương tác thủy động lực: Luận án sử dụng phương trình Morison với các hệ số $C_D, C_M$ không đổi dọc theo chiều dài cọc, chưa xét đến hiện tượng xoáy tách dòng cục bộ (Vortex-Induced Vibrations - VIV) và sóng phủ tương tác trực tiếp lên sàn đáy (Wave Slamming/Inundation).
  3. Quy mô mô hình thực nghiệm: Do điều kiện kích thước bể sóng 3D, mô hình thí nghiệm phải thu nhỏ theo tỷ lệ hình học Froude, dẫn đến hiệu ứng tỷ lệ Reynolds đối với lực nhớt của chất lỏng chưa được thỏa mãn tuyệt đối.

Chương trình nghiên cứu 10 năm tiếp theo được định hình qua 4 hướng mở rộng:

  • Hướng 1: Phát triển mô hình cấu tạo đàn - dẻo phá hủy (elasto-plastic damage model) đặc thù cho vật liệu san hô chịu tải trọng động trùng phục đa trục.
  • Hướng 2: Ứng dụng phương pháp tương tác chất lỏng - kết cấu đầy đủ (Two-way FSI - Fluid-Structure Interaction) kết hợp mô phỏng số thủy khí động học (CFD) để đánh giá sóng đổ và tải trọng va đập sóng cực hạn.
  • Hướng 3: Tích hợp thuật toán tối ưu hóa đa mục tiêu và hệ thống điều khiển dao động chủ động/bán chủ động (TMD, AMD) vào kết cấu nhà giàn DKI thế hệ mới.
  • Hướng 4: Nghiên cứu suy giảm mỏi và đánh giá độ tin cậy rủi ro theo thời gian thực (real-time structural health monitoring) dựa trên chuỗi số liệu cảm biến IoT gắn trên các công trình biển thực tế.

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

Nghiên cứu của luận án tạo ra chuỗi giá trị tác động sâu rộng:

  • Ảnh hưởng học thuật: Thiết lập một mốc chuẩn khoa học mới trong lĩnh vực Cơ học công trình biển tại Việt Nam; cung cấp công cụ phần mềm nội địa hóa B3DFC_2019 có độ chính xác cao tương đương các phần mềm thương mại quốc tế (SACS, SESAM, ANSYS) nhưng chuyên biệt hóa tối ưu cho nền địa chất san hô Biển Đông.
  • Chuyển đổi ngành công nghiệp và xây dựng công trình biển: Cung cấp phương pháp luận tính toán kiểm tra an toàn kết cấu cho hàng trăm công trình dân sự, dàn khoan dầu khí của Vietsovpetro, PVN và các công trình kinh tế biển đảo.
  • Ý nghĩa an ninh quốc phòng: Trực tiếp phục vụ công tác thiết kế, gia cường, nâng cấp hệ thống nhà giàn DKI trên thềm lục địa phía Nam và quần đảo Trường Sa, đảm bảo công trình đứng vững trước các siêu bão cấp 15-16, bảo vệ an toàn tính mạng cho cán bộ chiến sĩ làm nhiệm vụ bảo vệ chủ quyền biển đảo Tổ quốc.

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

  • Nghiên cứu sinh và học giả chuyên ngành Cơ kỹ thuật/Công trình biển: Tiếp cận hệ thống phương trình vi phân và ma trận độ cứng phi tuyến không gian đầy đủ, làm tài liệu tham khảo mẫu mực cho phương pháp số PTHH phi tuyến.
  • Các viện thiết kế và kỹ sư công trình biển (Viện Kỹ thuật Hải quân, Viện Thiết kế Bộ Quốc phòng, Viện Dầu khí Việt Nam): Ứng dụng trực tiếp thuật toán và phần mềm B3DFC_2019 để tính toán kiểm tra ổn định công trình mới và thiết kế giải pháp gia cố chân cọc các công trình hiện hữu.
  • Cơ quan quản lý nhà nước và hoạch định chính sách (Bộ Quốc phòng, Bộ Giao thông Vận tải, Bộ Tài nguyên và Môi trường): Có cơ sở khoa học định lượng để ban hành các tiêu chuẩn kỹ thuật quốc gia về khảo sát địa chất và thiết kế móng công trình trên nền san hô ngập nước.

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? Trả lời: Đóng góp độc đáo nhất là việc tích hợp thành công tiêu chuẩn ổn định động Budiansky - Roth vào bài toán dao động phi tuyến không gian 3D của hệ kết cấu thanh tương tác nền san hô có xét đến liên kết tiếp xúc một chiều (PTTX 3D). Luận án đã mở rộng lý thuyết ổn định công trình biển truyền thống (vốn chỉ xét ổn định tĩnh Euler hoặc khung phẳng ngàm cứng) thành lý thuyết ổn định động lực học phi tuyến toàn diện, xác định được ngưỡng phân kỳ bùng nổ chuyển vị khi tỷ số tải trọng đạt $P/P_{cr} = 1,0$.

2. Đột phá phương pháp luận của luận án so với các công trình quốc tế đi trước là gì? Trả lời: So với nghiên cứu của Abdel Raheem et al. [75] (chỉ sử dụng mô hình tuyến tính, không tương tác nền) và Kai Wei & Myers [58] (dùng mô hình dầm phẳng đơn giản), luận án đã thiết lập một hệ phần tử hữu hạn không gian 3 chiều đồng thời gồm: phần tử thanh phi tuyến hình học ($[K_{NL1}], [K_{NL2}], [K_{NL3}]$), phần tử khối nền san hô 8 nút, và phần tử tiếp xúc 3D mô tả ma sát Coulomb và hiện tượng tróc hở tiếp xúc ($k_{nz} \to 0$ khi chịu kéo). Toàn bộ được giải tự động qua chương trình tự lập trình B3DFC_2019 kết hợp thuật toán lặp Newmark - Newton Raphson chính xác.

3. Phát hiện bất ngờ nhất (counter-intuitive) rút ra từ dữ liệu mô phỏng và thực nghiệm là gì? Trả lời: Việc gia tải khối bê tông lớn lên đỉnh giàn (sàn công tác) tưởng chừng làm kết cấu nặng hơn và dễ sụp đổ hơn, nhưng thực tế lại giúp giảm chuyển vị ngang cực đại đỉnh giàn tới $39,13%$ dưới tải trọng bão gió. Khối lượng tập trung lớn đã tạo ra hiệu ứng quán tính chống lật mạnh mẽ, đóng vai trò như một bộ tiêu tán năng lượng thụ động làm ghìm đỉnh giàn, dù cái giá phải trả là mômen uốn chân cọc chính tăng cục bộ $14,8%$.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Luận án cung cấp đầy đủ, chi tiết và minh bạch toàn bộ ma trận giải tích phần tử trong các phụ lục toán học, quy luật quan hệ ứng suất - biến dạng của phần tử tiếp xúc, thông số cơ lý chi tiết 3 lớp nền san hô ngập nước tại Song Tử Tây, giản đồ xung gió $U_{win}(t)$ với $U_{max} = 46,35\text{ m/s}$, và sơ đồ khối thuật toán của chương trình B3DFC_2019, cho phép các nhà nghiên cứu độc lập tái lập 100% kết quả tính toán số.

5. Chương trình nghị sự nghiên cứu 10 năm tới giải quyết vấn đề cốt lõi nào? Trả lời: Giải quyết bài toán phá hủy mỏi tích lũy đa trục của liên kết nút ống thép dưới tác động ngẫu nhiên của phổ sóng biển Biển Đông, kết hợp mô hình hóa quá trình suy thoái cơ lý hóa cục bộ của nền đá san hô ngậm nước dưới tác động của biến đổi khí hậu và nước biển dâng.

Kết luận

  1. Thiết lập thành công mô hình toán học không gian 3D hoàn chỉnh mô tả đồng thời ứng xử phi tuyến hình học của kết cấu công trình biển hệ thanh và ứng xử tiếp xúc một chiều trượt ma sát Coulomb với nền san hô phân lớp ngậm nước.
  2. Phát triển bộ chương trình tính toán số chuyên dụng B3DFC_2019, tích hợp giải thuật Newmark phi tuyến và tiêu chuẩn ổn định động Budiansky - Roth, đã được kiểm chuẩn chính xác tuyệt đối qua nghiệm giải tích quốc tế và đo đạc thực nghiệm bể sóng 3D (sai số $< 2,89%$).
  3. Định lượng hóa vai trò chịu lực xung yếu của hệ cọc phụ xiên, chứng minh khả năng chia sẻ tới $34,6%$ mômen uốn chân công trình, làm cơ sở khoa học bắt buộc cho quy chuẩn thiết kế kháng bão nhà giàn DKI thế hệ mới.
  4. Khẳng định sự suy giảm nghiêm trọng của độ ổn định công trình khi nền san hô bão hòa nước, làm giảm tải trọng tới hạn $P_{cr}$ tới $28,4%$ và hạ thấp chiều cao sóng giới hạn an toàn.
  5. Mở ra 3 hướng nghiên cứu học thuật tiên phong: (i) Động lực học phá hủy nền móng san hô dưới tải trọng chu kỳ cao; (ii) Tương tác thủy khí động FSI phi tuyến bậc cao; (iii) Hệ thống giám sát sức khỏe kết cấu công trình biển thời gian thực phục vụ bảo vệ chủ quyền biển đảo bền vững.