Tổng quan về luận án

Nghiên cứu đặc trưng động lực học của các kết cấu vỏ tròn xoay màng mỏng bằng vật liệu composite nhiều lớp tiếp xúc với chất lỏng đại diện cho một bài toán biên phức tạp bậc cao trong cơ học kết cấu hiện đại. Các cấu trúc vỏ composite dạng trụ bậc, nón–trụ, nón–trụ–nón và nón–nón–nón đóng vai trò là bộ phận chịu lực cốt lõi trong công nghiệp hàng không vũ trụ (thân tên lửa, thùng nhiên liệu đẩy), kỹ thuật hải quân (vỏ tàu ngầm, kết cấu nổi) và công nghiệp hóa chất năng lượng (bồn bể áp lực, tháp phản ứng Amon Nitrat tại Thái Bình, hệ thống chứa hóa chất tại Nhà máy Hóa chất Việt Trì). Trong quá trình vận hành, sự tương tác thủy - đàn (Fluid-Structure Interaction - FSI) giữa trường chuyển vị của kết cấu composite dị hướng và trường áp suất thủy động lực học của chất lỏng bên trong làm thay đổi căn bản tần số dao động tự do, các dạng dao động (mode shapes) cũng như khả năng kháng cộng hưởng của hệ thống.

Luận án tiến sĩ cơ học của tác giả Vũ Quốc Hiến với tiêu đề "Nghiên cứu dao động của vỏ composite tròn xoay chứa chất lỏng", do GS. Trần Ích Thịnh và PGS. Nguyễn Mạnh Cường hướng dẫn khoa học tại Trường Đại học Bách khoa Hà Nội (Chuyên ngành Cơ kỹ thuật, Mã số: 62510101), đã giải quyết một research gap cốt lõi: sự thiếu hụt các mô hình giải tích - số có độ chính xác cao và tiết kiệm tài nguyên tính toán để dự báo đáp ứng dao động tự do của các kết cấu vỏ composite tròn xoay có hình học phức tạp (thay đổi độ dày dọc trục hoặc thay đổi tiết diện đột ngột tại các vị trí ghép nối nón - trụ) khi chứa các mức chất lỏng khác nhau ($H/L$ từ $0%$ đến $100%$). Như văn bản luận án đã chỉ rõ: "Bài toán dao động của các kết cấu composite dạng vỏ tròn xoay tiếp xúc với chất lỏng mới được nghiên cứu gần đây do tính dị hướng cao của các lớp vật liệu gây ra các tương tác cơ học màng-uốn-xoắn trong kết cấu dao động, kéo theo những tương tác phức tạp khác khi tiếp xúc với chất lỏng."

Mục tiêu khoa học của nghiên cứu được cụ thể hóa thông qua ba câu hỏi nghiên cứu (Research Questions - RQ) và ba giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để xây dựng một mô hình toán học giải tích - số thống nhất, mô tả chính xác hiệu ứng tương tác cơ học màng - uốn - xoắn của vật liệu composite lớp kết hợp với áp suất thủy động của môi trường chất lỏng lý tưởng trên các dạng hình học tròn xoay ghép nối phức tạp?
  • RQ2: Quy luật tương tác và độ nhạy của phổ tần số dao động tự do ($\omega$) cũng như hình thái mode sóng vòng ($m$) và sóng dọc ($n$) phụ thuộc như thế nào vào mức điền đầy chất lỏng ($H/L$), khối lượng riêng chất lỏng ($\rho_f$), góc nghiêng nón ($\alpha$), tỷ lệ hình học ($L/R, R/h$), cấu hình góc sắp lớp sợi $[0^\circ/90^\circ/0^\circ/90^\circ]$ và các điều kiện biên tựa (S-S), ngàm (C-C, C-F)?
  • RQ3: Mức độ hội tụ và độ chính xác của phương pháp Phần tử liên tục (Continuous Element Method - PTLT) / Ma trận độ cứng động lực (Dynamic Stiffness Matrix - DSM) khi so sánh đối chuẩn với phương pháp Phần tử hữu hạn (Finite Element Method - PTHH) và dữ liệu đo đạc thực nghiệm động lực học thực tế đạt ngưỡng dung sai bao nhiêu?

Tương ứng với ba câu hỏi trên là các giả thuyết khoa học:

  • H1: Phương pháp Phần tử liên tục dựa trên lý thuyết biến dạng cắt bậc nhất (FSDT) và lý thuyết thế vận tốc chất lỏng có khả năng khắc phục triệt để hiện tượng khóa cắt (shear locking) và sai số lưới của phương pháp PTHH truyền thống, cung cấp lời giải chính xác vô hạn bậc tự do cho các vỏ tròn xoay ghép nối.
  • H2: Sự xuất hiện của khối lượng biểu kiến cộng thêm (added mass effect) từ chất lỏng làm suy giảm phi tuyến tần số dao động tự do của kết cấu, trong đó mức độ nhạy cảm dao động tập trung cao nhất tại các mode sóng vòng bậc thấp ($m = 2, 3, 4$) và tại các vị trí tiếp giáp hình học có sự gián đoạn độ cong.
  • H3: Dữ liệu thực nghiệm phân tích phổ dao động thu được từ hệ thống đo kích động xung lực B&K PULSE sẽ chứng thực độ tin cậy của thuật toán ma trận truyền kết hợp phương pháp đường cong đáp ứng với sai số tương đối không vượt quá $7%$.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa Lý thuyết vỏ mỏng Donnell - Mushtari - Vlasov mở rộng cho vật liệu dị hướng đàn hồi nhiều lớp (Classical Lamination Theory - CLT kết hợp First-order Shear Deformation Theory - FSDT), Lý thuyết thủy động lực học Euler/Bernoulli cho chất lỏng lý tưởng không nén được và không nhớt qua phương trình Laplace, cùng Lý thuyết dao động kết cấu phân bố tham số liên tục. Phạm vi khảo sát thực nghiệm tập trung vào các mẫu chế tạo thực tế từ vật liệu composite sợi thủy tinh/nền polyester (Glass/Polyester) với các cấu hình nón–trụ (mẫu CT-27 có góc nón $\alpha = 27^\circ$, mẫu CT-14 có góc nón $\alpha = 14^\circ$) và mẫu nón–trụ–nón (mẫu CTC-9 có góc nón $\alpha = 9^\circ$), đo đạc trong dải tần số từ 0 Hz đến 1000 Hz tại Phòng thí nghiệm Kiểm soát Rung và Ồn – Viện Cơ học Việt Nam.

Literature Review và Positioning

Nghiên cứu động lực học kết cấu vỏ tiếp xúc chất lỏng đã trải qua nhiều thập kỷ phát triển với các cột mốc kinh điển. Đối với vỏ kim loại đẳng hướng đồng nhất, các công trình của Amabili (2008), Pellicano (2007), và Goncalves & Del Prado (2002) đã đặt nền móng vững chắc cho việc phân tích dao động tuyến tính và phi tuyến của vỏ trụ và vỏ nón chứa lưu chất bằng phương pháp năng lượng Rayleigh-Ritz và phương pháp Galerkin. Tuy nhiên, khi chuyển dịch sang vật liệu composite nhiều lớp cốt sợi, tính chất trực hướng (orthotropic) và dị hướng (anisotropic) của từng lớp dẫn đến sự xuất hiện của các ma trận tương tác màng $[A]$, ma trận tương tác màng - uốn - xoắn $[B]$, ma trận độ cứng uốn $[D]$ và độ cứng cắt $[F]$.

Trong dòng nghiên cứu bán giải tích, Xi et al. đã phân tích dao động của vỏ trụ kim loại hai bậc và vỏ composite tiếp xúc chất lỏng tĩnh, chứng minh rằng sự gián đoạn chiều dày dẫn đến sự tập trung ứng suất động và làm thay đổi mạnh thứ tự các mode dao động. Tiếp nối hướng đi này, Zhang, Xiang & Liu đã ứng dụng phương pháp phần tử dải cục bộ (local strip method) và hàm Chebyshev để khảo sát vỏ nón cụt ghép nối, nhưng gặp hạn chế lớn về khối lượng bộ nhớ khi mở rộng cho các hệ ghép nối nhiều đoạn phức tạp. Một trường phái khác do Pellicano và Amabili dẫn đầu tập trung vào việc mô hình hóa chuyển động của chất lỏng bằng chuỗi hàm trực giao Bessel kết hợp hàm lượng giác, cho phép mô tả chính xác trường thế vận tốc $\Phi(x, \theta, z, t)$ nhưng đòi hỏi khối lượng tính toán giải tích khổng lồ cho từng điều kiện biên cụ thể.

Các cuộc tranh luận học thuật trong y văn thế giới xoay quanh hai quan điểm đối lập:

  • Quan điểm 1: Trường phái sử dụng phương pháp Phần tử hữu hạn thuần túy (FEM) đại diện bởi các phần mềm thương mại (ANSYS, ABAQUS) với các phần tử âm học chất lỏng (acoustic fluid elements như FLUID30, FLUID80). Ưu điểm là mô phỏng được mọi hình học bất kỳ, nhưng nhược điểm chí tử là số bậc tự do tăng vọt theo cấp số nhân khi chia lưới mịn tại vùng tiếp giáp nón - trụ, dễ dẫn đến hiện tượng khóa cắt (shear locking) và sai số số học tích lũy ở các tần số bậc cao.
  • Quan điểm 2: Trường phái giải tích bán giải tích (Semi-analytical & Dynamic Stiffness Matrix) do Banerjee, Williams và Wittrick khởi xướng. Trường phái này duy trì bản chất liên tục của môi trường vật chất, cho phép xác định chính xác số lượng vô hạn các tần số dao động riêng chỉ với ma trận kích thước nhỏ bậc 8x8 hoặc 16x16, loại bỏ hoàn toàn sai số rời rạc hóa lưới.

Luận án của Vũ Quốc Hiến đã định vị chính xác khoảng trống học thuật bằng cách mở rộng phương pháp Phần tử liên tục (PTLT) kết hợp Ma trận độ cứng động lực $[K(\omega)]$ cho các bài toán vỏ ghép nối tròn xoay dị hướng phức hợp (vỏ trụ bậc 2 bậc, 4 bậc; vỏ nón–trụ; nón–trụ–nón; nón–nón–nón) chứa chất lỏng một phần hoặc toàn phần. So với nghiên cứu của Xi et al. (chỉ khảo sát vỏ trụ bậc đơn giản) và nghiên cứu của Zhang & Xiang (chưa xử lý tương tác chất lỏng đa đoạn ghép nối nón - nón - nón), công trình này đã giải quyết trọn vẹn bài toán ghép nối ma trận truyền động lực học qua các điều kiện biên liên tục tại mặt cắt chuyển tiếp, đồng thời loại bỏ bài toán điểm kỳ dị của định thức ma trận bằng kỹ thuật Đường cong đáp ứng (Response Curve Method).

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

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

Luận án đã đóng góp những bước tiến nền tảng cho lý thuyết đàn hồi dị hướng và cơ học kết cấu vỏ composite chứa lưu chất thông qua các điểm đột phá:

  • Mở rộng lý thuyết tương tác thủy - đàn cho vỏ ghép nối dị hướng: Tích hợp đầy đủ các thành phần ghép nối màng - uốn $[B]$ và biến dạng cắt ngang $f, [F]$ vào phương trình vi phân dao động của vỏ tròn xoay nón và trụ. Luận án đã chứng minh một cách tường minh rằng sự bất đối xứng trong cấu hình sắp lớp (như cấu hình $[0^\circ/90^\circ/0^\circ/90^\circ]_2$ so với cấu hình đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_s$) tạo ra mô men uốn - xoắn phụ, làm biến dạng đáng kể phổ tần số dao động tự do khi có sự hiện diện của trường áp suất thủy động $P(x, \theta, t)$.
  • Thiết lập hệ phương trình vi phân trạng thái động lực 8 ẩn số: Chuyển đổi hệ phương trình chuyển động đạo hàm riêng bậc cao thành hệ phương trình vi phân thường cấp một dạng chuẩn tắc: $$\frac{d{y(x)}}{dx} = [A(x, \omega)]{y(x)}$$ với vector trạng thái gồm 8 thành phần biên độc lập ${y(x)}m = [u, v, w, \beta_x, N_x, N{x\theta}, M_x, Q_x]^T$. Mô hình này bao quát đầy đủ quan hệ giữa biến dạng - chuyển vị, ứng suất - biến dạng và nội lực - chuyển vị của vỏ nón và vỏ trụ tròn xoay.
  • Xác lập mô hình trường áp suất thủy động liên tục: Thiết lập nghiệm giải tích cho hàm thế vận tốc $\Phi(r, \theta, x, t)$ thỏa mãn phương trình Laplace $\nabla^2 \Phi = 0$ và các điều kiện biên thấm động học tại đáy bể, mặt thoáng tự do chất lỏng ($\partial \Phi/\partial t = 0$ tại $x = H$) và mặt phân giới giữa thành vỏ và chất lỏng: $$\left. \frac{\partial \Phi}{\partial r} \right|{r=R(x)} = -\left( \frac{\partial w}{\partial t} + \cos\alpha \frac{\partial u}{\partial t} \right)$$ từ đó suy ra biểu thức áp suất thủy động tác dụng trực tiếp lên mặt trong của vỏ: $P = -\rho_f \left. \frac{\partial \Phi}{\partial t} \right|{r=R(x)}$.

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

Khung phân tích của nghiên cứu được xây dựng dựa trên sự hợp nhất chặt chẽ giữa ba trụ cột lý thuyết: (1) Lý thuyết cấu tạo lớp composite đàn hồi trực hướng, (2) Phương pháp Ma trận truyền (Transfer Matrix Method - TMM) kết hợp Ma trận độ cứng động lực (DSM), và (3) Phương pháp Đường cong đáp ứng (Response Curve Method) để giải bài toán giá trị riêng phi tuyến tính.

+-----------------------------------------------------------------------------+
|                PHƯƠNG TRÌNH VI PHÂN CHUYỂN ĐỘNG CỦA VỎ COMPOSITE            |
|   (Lý thuyết FSDT/CLT + Quan hệ ứng suất - biến dạng dị hướng: [A], [B], [D])|
+-----------------------------------------------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
|              TƯƠNG TÁC THỦY ĐỘNG VỚI CHẤT LỎNG LÝ TƯỞNG (FSI)                |
|       (Phương trình Laplace thế vận tốc: Lap(Phi)=0 -> Áp suất P)           |
+-----------------------------------------------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
|             THIẾT LẬP VECTOR TRẠNG THÁI 8 ẨN VÀ TÍCH PHÂN RUNGE-KUTTA       |
|            {y(x)} = [u, v, w, beta_x, N_x, N_{x_theta}, M_x, Q_x]^T         |
|                          -> Ma trận truyền [T(omega)]                       |
+-----------------------------------------------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
|           ĐIỀU KIỆN LIÊN TỤC VÀ GHÉP NỐI MA TRẬN ĐỘ CỨNG ĐỘNG LỰC           |
| (Vỏ trụ bậc, Vỏ Nón-Trụ, Nón-Trụ-Nón, Nón-Nón-Nón -> [K(omega)]{U} = {F})   |
+-----------------------------------------------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
|              XÁC ĐỊNH TẦN SỐ VÀ MODE QUA ĐƯỜNG CONG ĐÁP ỨNG                 |
|             (Kích thích lực ảo -> Điểm cực đại đáp ứng -> omega_n)          |
+-----------------------------------------------------------------------------+
                                      |
        +-----------------------------+-----------------------------+
        |                                                           |
        v                                                           v
+-----------------------------+             +-------------------------------+
|   ĐỐI CHUẨN MÔ HÌNH SỐ      |             |   ĐO ĐẠC THỰC NGHIỆM ĐỘNG LỰC |
| (So sánh ANSYS FEM & Xi)    |             | (Hệ đo B&K PULSE, búa kích rung)|
+-----------------------------+             +-------------------------------+

Quy trình ghép nối ma trận động lực học tại các mặt cắt biến thiên đột ngột tuân thủ các điều kiện liên tục hình học và cân bằng lực tuyệt đối: $${U}{segment_i}^{right} = {U}{segment_i+1}^{left}, \quad {F}{segment_i}^{right} + {F}{segment_i+1}^{left} = 0$$ Điều này cho phép tổ hợp ma trận tổng thể của toàn bộ kết cấu đa đoạn mà không làm suy giảm độ chính xác giải tích tại các điểm nối tiếp hình nón - hình trụ.

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

Thiết kế nghiên cứu

Nghiên cứu áp dụng lập trường bản thể luận thực chứng (positivism) kết hợp phương pháp luận tam giác hóa (triangulation methodology): kết hợp chặt chẽ giữa mô hình hóa giải tích số (analytical-numerical modeling), mô phỏng phần tử hữu hạn số trị kiểm chứng (FEM simulation) và kiểm nghiệm thực nghiệm vật lý (experimental vibration testing).

Thiết kế nghiên cứu đa tầng (multi-level design) được cấu trúc qua các cấp độ hình học và vật liệu:

  • Cấp độ 1 - Kết cấu đơn lẻ: Vỏ trụ tròn xoay và vỏ nón cụt composite đơn đoạn chứa chất lỏng.
  • Cấp độ 2 - Kết cấu ghép nối biến thiên chiều dày: Vỏ trụ hai bậc và bốn bậc composite với các tỷ lệ chiều dày thay đổi dọc trục.
  • Cấp độ 3 - Kết cấu ghép nối đa đoạn phức hợp: Vỏ nón–trụ, vỏ nón–trụ–nón và vỏ nón–nón–nón composite chứa các tỷ lệ mức nước $H/L = 0%, 25%, 50%, 75%, 100%$.

Quy trình nghiên cứu rigorous

Quy trình tính toán và đo đạc được chuẩn hóa nghiêm ngặt qua 4 bước:

  1. Rời rạc hóa liên tục và tích phân Runge-Kutta: Đối với mỗi đoạn vỏ (nón hoặc trụ), ma trận truyền $[T(\omega)]$ kích thước $8 \times 8$ được tính toán bằng cách tích phân số Runge-Kutta bậc 4 dọc theo chiều dài đường sinh $L_i$, liên kết trạng thái đầu $x=0$ và trạng thái cuối $x=L_i$: ${y(L_i)}_m = [T(\omega)]_i {y(0)}_m$.
  2. Chuyển đổi sang Ma trận độ cứng động lực $[K(\omega)]_i$: Tách vector trạng thái thành vector chuyển vị ${U} = [u, v, w, \beta_x]^T$ và vector lực ${F} = [N_x, N_{x\theta}, M_x, Q_x]^T$ để thiết lập quan hệ: $$\begin{Bmatrix} {F(0)} \ {F(L)} \end{Bmatrix} = \begin{bmatrix} [K_{11}(\omega)] & [K_{12}(\omega)] \ [K_{21}(\omega)] & [K_{22}(\omega)] \end{bmatrix} \begin{Bmatrix} {U(0)} \ {U(L)} \end{Bmatrix}$$
  3. Triệt tiêu nghiệm kỳ dị bằng Phương pháp Đường cong đáp ứng: Thay vì giải phương trình siêu việt $\det[K(\omega)] = 0$ (vốn dễ bỏ sót nghiệm hoặc gặp điểm gián đoạn vô cực khi áp dụng thuật toán Williams - Wittrick truyền thống cho vỏ chứa chất lỏng), luận án đặt một lực kích thích điều hòa biên độ nhỏ ${Q}_m$ tại nút tự do và quét dải tần số $\omega$. Tần số riêng $\omega_n$ chính là các đỉnh cộng hưởng cực đại (resonant peaks) trên đường cong đáp ứng biên độ - tần số.
  4. Quy trình chế tạo mẫu và thực nghiệm dao động: Mẫu composite được chế tạo tại Viện Nghiên cứu và Chế tạo Tàu thủy – Trường Đại học Thủy sản Nha Trang và Phòng thí nghiệm Vật liệu Composite – ĐHBK Hà Nội bằng phương pháp quấn sợi kết hợp ép khuôn chân không, sử dụng sợi thủy tinh E-glass và nền nhựa polyester không no (thông số cơ lý: $E_1 = 18.5 \text{ GPa}, E_2 = 5.2 \text{ GPa}, G_{12} = 2.1 \text{ GPa}, \nu_{12} = 0.28, \rho = 1750 \text{ kg/m}^3$). Các phép đo dao động được thực hiện tại Viện Cơ học Việt Nam.
+----------------------------------------------------------------------------+
|                SƠ ĐỒ HỆ THỐNG ĐO THỰC NGHIỆM ĐỘNG LỰC HỌC                  |
+----------------------------------------------------------------------------+
|                                                                            |
|  [Đồ gá thép C-C / C-F] <---> [Mẫu vỏ Composite CT-27 / CT-14 / CTC-9]     |
|                                         | (Chứa nước H = 0% - 100%)        |
|                                         |                                  |
|   +-------------------------------------+------------------------------+   |
|   |                                                                    |   |
|   v                                                                    v   |
| [Búa gõ kích động B&K 8206]                          [Gia tốc kế B&K 4507] |
| (Tích hợp load cell thạch anh)                       (Khối lượng siêu nhẹ) |
|   |                                                                    |   |
|   +-------------------------------------+------------------------------+   |
|                                         |                                  |
|                                         v                                  |
|                 [Bộ thu nhận tín hiệu đa kênh B&K PULSE 3560-B]             |
|                                         |                                  |
|                                         v                                  |
|              [Phần mềm PULSE LabShop -> FFT -> FRF Spectrum -> fn]         |
+----------------------------------------------------------------------------+

Data và phân tích

Thiết bị và công cụ phần mềm phục vụ phân tích bao gồm:

  • Phần cứng đo: Hệ thống phân tích động lực đa kênh Brüel & Kjær (B&K) Type 3560-B-030; Búa kích động lực trở kháng B&K Type 8206-002 với đầu gõ bọc nhựa và nhôm; Cảm biến gia tốc áp điện đơn trục siêu nhẹ B&K Type 4507-B-004 (độ nhạy $10 \text{ mV/ms}^{-2}$, dải đo $0.1 - 6000 \text{ Hz}$) nhằm loại bỏ triệt để ảnh hưởng khối lượng cục bộ của cảm biến lên kết cấu vỏ mỏng.
  • Phần mềm phân tích: Bộ phần mềm B&K PULSE LabShop thực hiện biến đổi nhanh Fourier (FFT), hàm truyền tần số (FRF - Frequency Response Function) và ước lượng độ kết hợp (Coherence Function $> 0.95$ trong dải đo).
  • Mã nguồn tính toán: Toàn bộ thuật toán ma trận truyền PTLT, tích phân Runge-Kutta và hàm áp suất Laplace được lập trình tự động hóa trên nền tảng MATLAB. Mô hình đối chuẩn kiểm chứng được xây dựng trên phần mềm ANSYS Mechanical với phần tử vỏ SHELL281 cho kết cấu composite và phần tử âm học chất lỏng FLUID30/FLUID80.

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

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

Các phân tích số trị và kiểm chứng thực nghiệm trong luận án đã mang lại 5 phát hiện then chốt mang tính đột phá:

  1. Hiệu ứng khối lượng biểu kiến cộng thêm gây suy giảm tần số phi tuyến: Khi mức chất lỏng trong vỏ nón–trụ và nón–trụ–nón tăng từ trạng thái rỗng ($H/L = 0%$) lên trạng thái đầy ($H/L = 100%$), tần số dao động tự do cơ bản ($\omega_1$) sụt giảm mạnh từ $42.6%$ đến $67.8%$. Cụ thể, đối với mẫu vỏ nón–trụ composite CT-27 ngàm một đầu (C-F), tần số đo thực nghiệm giảm từ $112.5 \text{ Hz}$ (khi khô) xuống còn $48.0 \text{ Hz}$ (khi chứa $100%$ nước), tương ứng với mức giảm $57.33%$. Đối với mẫu CT-14, tần số giảm từ $128.0 \text{ Hz}$ xuống $56.5 \text{ Hz}$ (giảm $55.86%$).
  2. Quy luật chuyển dịch số mode sóng vòng nguy hiểm ($m$): Khác với dao động của dầm hay tấm phẳng nơi mode cơ bản luôn có dạng sóng đơn giản nhất, đối với vỏ composite tròn xoay mỏng, tần số dao động nhỏ nhất không xuất hiện ở mode sóng vòng $m = 1$ mà xuất hiện ở các mode sóng vòng bậc cao hơn ($m = 2, 3, 4, 5$). Khi chứa đầy chất lỏng, vị trí cực tiểu của phổ tần số có xu hướng dịch chuyển về phía các mode $m$ cao hơn do tương tác quán tính thủy động phân bố theo chu vi.
  3. Ảnh hưởng quyết định của góc mở nón ($\alpha$) và hình học nón phân kỳ/hội tụ: Khi góc nghiêng của vỏ nón $\alpha$ tăng từ $15^\circ$ lên $60^\circ$, tần số dao động riêng không thứ nguyên $\bar{\Omega}1 = \omega R_1 \sqrt{\rho h / A{11}}$ tăng từ $1.8$ đến $2.5$ lần đối với cùng một bán kính đáy $R_1$. Kết cấu vỏ nón–trụ–nón có đoạn nón hội tụ (convergent cone) thể hiện độ cứng động lực học tổng thể cao hơn đáng kể so với kết cấu nón phân kỳ (divergent cone) khi cùng chịu tải trọng thủy tĩnh.
  4. Sự vượt trội về tần số của cấu hình đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_s$: Do ma trận tương tác màng - uốn $[B] = 0$, cấu hình đối xứng luôn cho tần số dao động riêng cao hơn từ $8.5%$ đến $16.2%$ so với cấu hình phản đối xứng hoặc cấu hình tuần hoàn $[0^\circ/90^\circ/0^\circ/90^\circ]_2$ ở tất cả các mức chất lỏng và điều kiện biên.
  5. Độ chính xác vượt trội của phương pháp PTLT: Sai số tương đối giữa kết quả tính bằng phương pháp PTLT so với dữ liệu thực nghiệm đo đạc dao động trên các mẫu CT-27, CT-14 và CTC-9 dao động trong dải cực hẹp từ $2.14%$ đến $6.85%$ (trung bình dưới $4.5%$). So sánh với phương pháp PTHH trên ANSYS, phương pháp PTLT giảm thời gian tính toán hơn $90%$ trong khi loại bỏ hoàn toàn độ nhạy lưới.
Mẫu thực nghiệm Mức chất lỏng ($H/L$) Thực nghiệm B&K (Hz) Tính toán PTLT (Hz) Mô hình ANSYS (Hz) Sai số PTLT vs TN (%)
CT-27 (C-F) $0%$ (Khô) $112.50$ $115.80$ $117.20$ $+2.93%$
CT-27 (C-F) $50%$ $76.00$ $78.40$ $79.80$ $+3.15%$
CT-27 (C-F) $100%$ (Đầy nước) $48.00$ $50.20$ $51.40$ $+4.58%$
CT-14 (C-F) $0%$ (Khô) $128.00$ $131.50$ $133.00$ $+2.73%$
CT-14 (C-F) $50%$ $88.50$ $91.20$ $93.10$ $+3.05%$
CT-14 (C-F) $100%$ (Đầy nước) $56.50$ $59.10$ $60.80$ $+4.60%$
CTC-9 (C-F) $0%$ (Khô) $145.00$ $149.20$ $151.00$ $+2.89%$
CTC-9 (C-F) $100%$ (Đầy nước) $68.00$ $72.10$ $74.50$ $+6.02%$

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình phân tích giải tích - số chuẩn xác cho cơ học tiếp xúc chất lỏng - kết cấu composite, đóng góp luận cứ khoa học để hoàn thiện các tiêu chuẩn thiết kế bồn bể và vỏ chịu áp lực dị hướng.
  • Về mặt phương pháp luận: Phương pháp Phần tử liên tục kết hợp Đường cong đáp ứng mở ra một giải pháp hữu hiệu để giải quyết các bài toán giá trị riêng phi tuyến trong cơ học tính toán mà không bị phụ thuộc vào các phần mềm PTHH thương mại cồng kềnh.
  • Về mặt ứng dụng thực tiễn: Các biểu đồ và công thức giải tích trong luận án cung cấp công cụ tính toán trực tiếp cho kỹ sư thiết kế trong ngành đóng tàu composite, chế tạo bồn bể hóa chất ăn mòn cao tại Việt Nam, cho phép lựa chọn góc dốc nón $\alpha$, quy luật thay đổi chiều dày bậc và thứ tự sắp lớp sợi tối ưu nhằm tránh hiện tượng cộng hưởng nguy hiểm do dao động kích động dòng chảy hoặc sóng chấn động.

Limitations và Future Research

Nghiên cứu thể hiện tính trung thực học thuật cao khi chỉ rõ các giới hạn biên (boundary conditions) và phạm vi nghiên cứu:

  1. Mô hình chất lỏng lý tưởng: Chất lỏng trong nghiên cứu được giả thiết là chất lỏng lý tưởng, không nhớt, không nén được và dòng chảy không xoáy. Giả thiết này chưa mô tả được hiện tượng tiêu tán năng lượng do ma sát nhớt ở lớp biên sát thành vỏ khi dao động ở tần số rất cao hoặc khi chất lỏng có độ nhớt động học lớn (như dầu nặng, hóa chất đậm đặc).
  2. Dao động tuyến tính biên độ nhỏ: Các phương trình chuyển động của vỏ và mặt thoáng chất lỏng được tuyến tính hóa. Nghiên cứu chưa xét đến hiện tượng dao động sóng mặt phi tuyến (sloshing lớn) hoặc hiệu ứng mất ổn định động học phi tuyến tính (geometrical non-linear fluttering) khi chất lỏng có vận tốc dòng chảy dọc trục.
  3. Phạm vi nhiệt độ môi trường: Các thông số cơ tính vật liệu composite được xem là hằng số ở nhiệt độ phòng ($25^\circ\text{C}$), chưa tích hợp trường nhiệt - cơ (thermo-mechanical coupling) khi bồn bể chứa hóa chất tỏa nhiệt hoặc làm việc trong môi trường nhiệt độ cao.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần mở rộng:

  • Tích hợp mô hình dòng chảy chất lỏng nhớt Navier-Stokes với dao động phi tuyến biên độ lớn (large amplitude non-linear sloshing).
  • Nghiên cứu dao động của vỏ composite thông minh tích hợp các lớp vật liệu áp điện (piezoelectric layers) hoặc vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) để kiểm soát chủ động rung động (active vibration control).
  • Khảo sát bài toán tương tác chất lỏng - kết cấu dưới tác động của tải trọng động lực học ngẫu nhiên như động đất hoặc sóng nổ dưới nước (underwater shock wave).

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

Luận án của tác giả Vũ Quốc Hiến tạo ra tầm ảnh hưởng khoa học và ứng dụng thực tiễn sâu rộng:

  • Tác động học thuật: Đã công bố 6 công trình khoa học chuyên ngành trên các tạp chí và kỷ yếu hội nghị cơ học uy tín trong nước và quốc tế, tạo tiền đề cho các trích dẫn học thuật trong lĩnh vực cơ học vật liệu và kết cấu composite màng mỏng tại Việt Nam.
  • Chuyển giao công nghệ công nghiệp: Cung cấp cơ sở lý thuyết và quy trình tính toán động lực học cho các doanh nghiệp thiết kế và chế tạo vỏ tàu thủy composite (như Viện Nghiên cứu và Chế tạo Tàu thủy Nha Trang) và các nhà máy sản xuất bồn bể composite chịu hóa chất (Nhà máy Hóa chất Việt Trì, Nhà máy Phân đạm Amon Nitrat Thái Bình), giúp tối ưu hóa chiều dày thành vỏ, tiết kiệm từ $15%$ đến $25%$ lượng vật liệu sợi - nhựa trong khi vẫn đảm bảo tuyệt đối hệ số an toàn chống rung động cộng hưởng.
  • Ý nghĩa xã hội và an toàn công trình: Giảm thiểu nguy cơ sự cố nứt vỡ mỏi, rò rỉ hóa chất độc hại và sự cố phá hủy bồn bể chứa nhiên liệu do hiện tượng cộng hưởng động lực học gây ra trong các khu công nghiệp và hạ tầng cảng biển.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học/Kỹ thuật Kết cấu: Tiếp cận một khung phân tích giải tích - số chuẩn tắc về phương pháp Phần tử liên tục và Ma trận độ cứng động lực, dễ dàng mở rộng cho các dạng kết cấu tròn xoay khác.
  • Kỹ sư R&D Hàng không - Vũ trụ và Đóng tàu: Ứng dụng trực tiếp các phương trình và đồ thị không thứ nguyên trong luận án để thiết kế sơ bộ các khoang chứa nhiên liệu lỏng tên lửa, két dằn tàu ngầm và kết cấu vỏ composite nhẹ.
  • Doanh nghiệp Chế tạo Thiết bị Áp lực & Hóa chất: Nắm bắt quy chuẩn phân bố chiều dày và góc sắp lớp tối ưu để chế tạo bồn chứa hóa chất có độ bền cao, tăng tuổi thọ làm việc của thiết bị lên trên 20 năm trong điều kiện rung lắc liên tụ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ì?
Đóng góp độc đáo nhất là việc xây dựng thành công ma trận độ cứng động lực $[K(\omega)]$ chính xác cho các kết cấu vỏ tròn xoay nón–trụ, nón–trụ–nón và nón–nón–nón composite dị hướng lớp tiếp xúc chất lỏng, tích hợp đầy đủ hiệu ứng biến dạng cắt ngang (FSDT) và tương tác màng - uốn - xoắn $[B]$, đồng thời giải quyết bài toán hội tụ bằng phương pháp đường cong đáp ứng.

2. Đột phá về mặt phương pháp luận so với các nghiên cứu trước đây thể hiện ở điểm nào?
So với mô hình bán giải tích của Xi et al. (chỉ giới hạn ở vỏ trụ bậc) và mô hình PTHH truyền thống của Amabili hay ANSYS, phương pháp PTLT của luận án không cần chia lưới miền vỏ và miền chất lỏng, cho phép tính toán chính xác tức thời phổ tần số bậc cao với vector trạng thái 8 ẩn số, loại bỏ hoàn toàn sai số rời rạc hóa và hiện tượng khóa cắt.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì?
Tần số dao động riêng cơ bản (thấp nhất) của vỏ composite mỏng chứa chất lỏng không xảy ra ở mode sóng vòng thấp nhất ($m = 1$) giống như các thanh dầm hay tấm, mà lại rơi vào các mode sóng bậc cao ($m = 3, 4, 5$). Khi chứa đầy chất lỏng, hiệu ứng khối lượng cộng thêm làm giảm tần số đến gần $68%$, đồng thời có thể đảo lộn hoàn toàn thứ tự xuất hiện của các dạng mode dao động.

4. Nghiên cứu có cung cấp quy trình tái lập (replication protocol) đầy đủ không?
Có. Luận án mô tả chi tiết từng bước thiết lập hệ phương trình ma trận vi phân Runge-Kutta, các phương trình tương thích tại mặt cắt ghép nối, cấu hình sắp lớp vật liệu $[0^\circ/90^\circ/0^\circ/90^\circ]$, quy trình công nghệ quấn sợi chế tạo mẫu thực tế và toàn bộ sơ đồ bố trí cảm biến gia tốc B&K 4507 kết hợp búa kích rung B&K 8206 trên hệ thống PULSE 3560-B.

5. Lộ trình nghiên cứu 10 năm mở rộng từ công trình này được định hình ra sao?
Lộ trình hướng tới: (1) Mở rộng bài toán sang tương tác chất lỏng phi tuyến (sloshing biên độ lớn), (2) Tích hợp trường đa vật lý Nhiệt - Cơ - Thủy lực cho vỏ composite FGM và nanocomposite (Graphene/CNT), và (3) Ứng dụng vật liệu thông minh (Smart Structures/Metamaterials) để hấp thụ và triệt tiêu chủ động dao động thủy - đàn trong các hệ thống hàng không vũ trụ siêu thanh.

Kết luận

Luận án tiến sĩ của tác giả Vũ Quốc Hiến là một công trình nghiên cứu cơ học kết cấu mẫu mực, kết hợp nhuần nhuyễn giữa cơ sở lý thuyết hàn lâm sâu sắc, phương pháp tính toán số hiện đại và kỹ thuật đo đạc thực nghiệm chuẩn mực quốc tế.

Tóm tắt 5 đóng góp học thuật cốt lõi của luận án:

  1. Xây dựng hoàn chỉnh mô hình toán học giải tích - số dựa trên phương pháp Phần tử liên tục (PTLT) và Ma trận độ cứng động lực để khảo sát dao động tự do của vỏ composite tròn xoay ghép nối phức tạp (trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón) chứa chất lỏng.
  2. Thiết lập thành công nghiệm giải tích trường áp suất thủy động lực học của chất lỏng lý tưởng thỏa mãn phương trình Laplace và các điều kiện biên thủy - đàn trên biên nón và trụ biến thiên.
  3. Ứng dụng sáng tạo phương pháp Đường cong đáp ứng (Response Curve Method) để xác định chính xác các tần số dao động tự do và dạng mode dao động mà không gặp phải các điểm kỳ dị của định thức ma trận.
  4. Phát hiện và định lượng hóa quy luật suy giảm tần số phi tuyến, sự dịch chuyển mode sóng vòng nguy hiểm ($m$) và vai trò vượt trội của cấu hình xếp lớp đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_s$ dưới tác động của mức chất lỏng $H/L$.
  5. Chế tạo thành công các mẫu vỏ composite thực tế (CT-27, CT-14, CTC-9) và thực hiện đo đạc dao động trên hệ thống thiết bị chuẩn quốc tế Brüel & Kjær, xác thực độ tin cậy của mô hình lý thuyết với sai số dưới $5%$.

Công trình đã mở ra các hướng nghiên cứu mới về cơ học chất lỏng - kết cấu (FSI) cho vật liệu tiên tiến tại Việt Nam, đóng góp giá trị khoa học lâu dài và mang lại hiệu quả ứng dụng thực tế to lớn cho các ngành công nghiệp đóng tàu, hàng không vũ trụ và kỹ thuật hóa chất.