Tổng quan về luận án

Trong kỷ nguyên phát triển của kỹ thuật hàng không vũ trụ, công nghiệp đóng tàu, kỹ thuật hóa chất và năng lượng hạt nhân, các kết cấu vỏ tròn xoay mỏng và dày vừa bằng vật liệu composite nhiều lớp ngày càng giữ vị trí then chốt. Sự kết hợp giữa tỷ số độ bền trên khối lượng vượt trội, khả năng kháng ăn mòn hóa học cao và tính linh hoạt trong thiết kế định hướng sợi giúp vỏ composite dần thay thế hoàn toàn các kết cấu hợp kim truyền thống trong các bồn chứa nhiên liệu tên lửa, tháp chưng cất, bồn phản ứng hóa chất và vỏ tàu ngầm. Tuy nhiên, khi các kết cấu này tiếp xúc hoặc chứa chất lỏng, bài toán tương tác chất lỏng - kết cấu (Fluid-Structure Interaction - FSI) làm biến đổi sâu sắc đặc tính động lực học của hệ thống. Luận án tiến sĩ kỹ thuật cơ học của nghiên cứu sinh Vũ Quốc Hiến dưới sự hướng dẫn khoa học của GS. Trần Ích Thịnh và PGS. Nguyễn Mạnh Cường tại Trường Đại học Bách khoa Hà Nội đã giải quyết bài toán cốt lõi này thông qua đề tài: "Nghiên cứu dao động của vỏ composite tròn xoay chứa chất lỏng".

                +-------------------------------------------------------------+
                |     PHÂN TÍCH DAO ĐỘNG VỎ COMPOSITE CHỨA CHẤT LỎNG          |
                +-------------------------------------------------------------+
                                               |
                     +-------------------------+-------------------------+
                     |                                                   |
                     v                                                   v
   +-----------------------------------+               +-----------------------------------+
   |   MÔ HÌNH LÝ THUYẾT & SỐ (PTLT)   |               |      NGHIÊN CỨU THỰC NGHIỆM       |
   +-----------------------------------+               +-----------------------------------+
   | - Lý thuyết vỏ Reissner-Mindlin   |               | - Mẫu GFRP [0°/90°/0°/90°]        |
   | - Phương trình Laplace & Bessel   |               | - Kết cấu: CT-14, CT-27, CTC-9    |
   | - Vector trạng thái 10 thành phần |               | - Cảm biến gia tốc & búa kích đập |
   | - Ma trận độ cứng động lực K(ω)   |               | - Kiểm định sai số FSI            |
   +-----------------------------------+               +-----------------------------------+
                     |                                                   |
                     +-------------------------+-------------------------+
                                               |
                                               v
                +-------------------------------------------------------------+
                |  ĐÁNH GIÁ ĐỊNH LƯỢNG ẢNH HƯỞNG FSI & THIẾT KẾ KHÁNG CHẤN    |
                +-------------------------------------------------------------+

Khoảng trống học thuật (research gap) lớn nhất được xác định là: Mặc dù các mô hình giải tích và phần tử hữu hạn (PTHH) cổ điển đã phân tích tương đối đầy đủ dao động của vỏ trụ tròn hoặc vỏ nón đơn lẻ bằng kim loại (Amabili, 1996; Brenneman & Yang, 1992), các nghiên cứu định lượng về dao động tự do của các hệ vỏ tròn xoay liên hợp composite có mặt cắt ngang thay đổi hoặc thay đổi chiều dày phân bậc (như nón - trụ, nón - trụ - nón, nón - nón - nón và trụ bậc) chứa các mức chất lỏng khác nhau vẫn chưa được khám phá một cách hệ thống bằng phương pháp giải tích - số chính xác cao. Phương pháp PTHH thông thường (Finite Element Method - FEM) thường đòi hỏi chia lưới cực mịn, dẫn đến số bậc tự do khổng lồ và dễ gặp sai số tích lũy hoặc hiện tượng khóa cắt (shear locking) ở miền tần số cao, trong khi các nghiệm giải tích thuần túy chỉ áp dụng được cho hình học đối xứng lý tưởng đơn giản.

Luận án thiết lập 5 câu hỏi và giả thuyết nghiên cứu cụ thể:

  1. $H_1$: Liệu việc ứng dụng phương pháp phần tử liên tục (Continuous Element Method - PTLT) kết hợp với lý thuyết vỏ biến dạng cắt bậc nhất Reissner-Mindlin (FSDT) có thể biểu diễn chính xác ma trận độ cứng động lực $[K(\omega)]$ với số bậc tự do hữu hạn trên miền liên tục vô số bậc tự do cho kết cấu vỏ composite dị hướng phức tạp?
  2. $H_2$: Mức độ suy giảm tần số dao động tự do của vỏ composite phân bậc và vỏ liên hợp nón - trụ khi mức chất lỏng $H/L$ tăng từ $0%$ (trạng thái khô) đến $100%$ (chứa đầy) có tuân theo quy luật phi tuyến phụ thuộc vào số mode sóng vòng $m$ và số mode sóng dọc $n$ hay không?
  3. $H_3$: Hiện tượng tương tác màng - uốn - xoắn do tính dị hướng của các lớp vật liệu composite trực hướng (orthotropic layers) tác động như thế nào đến sự chuyển đổi mode dao động cơ bản (fundamental mode shift) khi có sự hiện diện của khối lượng chất lỏng bổ sung (hydrodynamic added mass)?
  4. $H_4$: Thuật toán quét đường cong đáp ứng tần số dựa trên kích thích lực đơn vị có thể thay thế hoàn toàn thuật toán Wittrick-Williams trong việc tìm nghiệm tần số riêng mà không bỏ sót nghiệm hay gặp khó khăn tính toán do hàm siêu việt phân kỳ?
  5. $H_5$: Độ chính xác của mô hình lý thuyết PTLT có tương thích với dữ liệu đo thực nghiệm thực tế trên các mẫu vỏ composite sợi thủy tinh/polyester (GFRP) đa phân đoạn chứa nước ở các điều kiện biên ngàm - tự do (C-F) và ngàm - ngàm (C-C)?

Khung lý thuyết của công trình tích hợp lý thuyết vỏ Reissner-Mindlin (hệ số hiệu chỉnh cắt $f = 5/6$), lý thuyết dòng thế phi nhớt không nén được Laplace - Bernoulli, phương pháp vi phân ma trận truyền (Differential Transfer Matrix Method) và thuật toán ma trận độ cứng động lực (Dynamic Stiffness Matrix - DSM). Phạm vi nghiên cứu bao quát các kết cấu vỏ tròn xoay liên hợp chế tạo từ composite sợi thủy tinh/polyester không no, làm việc trong giới hạn đàn hồi tuyến tính, biến dạng bé, với các góc nghiêng nón biến thiên $\alpha \in [0^\circ, 60^\circ]$ và các cấu hình sắp xếp lớp cân bằng và không đối xứng như $[0^\circ/90^\circ/0^\circ/90^\circ]_s$ và $[0^\circ/90^\circ/0^\circ/90^\circ]_2$.

Literature Review và Positioning

Các nghiên cứu động lực học công trình vỏ kim loại tiếp xúc với chất lỏng đã bắt đầu phát triển mạnh mẽ từ giữa thế kỷ 20. Khởi xướng bởi Warburton và Al-Najafi (1969) với phương pháp PTHH cho vỏ trụ hai bậc, tiếp nối bởi Chang (1979) bằng phương pháp Rayleigh-Ritz phân đoạn. Trong lĩnh vực tương tác FSI cho vỏ trụ kim loại, Brenneman và Yang (1992) đã tiên phong kết hợp phương pháp ma trận độ cứng động lực và mô hình PTHH ghép nối, trong khi Han và Liu (1994) phân tích bình trụ kim loại có chiều dày biến đổi dựa trên lý thuyết vỏ mỏng Flügge và ma trận truyền. Amabili (1996, 2002) đã đóng góp những phân tích kinh điển về dao động tự do và phi tuyến của vỏ trụ kim loại tựa đơn giản chứa một phần chất lỏng. Về hình học vỏ nón kim loại, Damatty cùng các cộng sự (2001, 2002, 2005) và Senthil Kumar & Ganesan (2008) đã công bố các kết quả bán giải tích và thực nghiệm toàn diện về vỏ nón cụt chịu tải trọng thủy động do chất lỏng chảy và kích động động đất.

TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT ĐỘNG LỰC HỌC VỎ CHỨA CHẤT LỎNG
===================================================================================
1960s - 1980s:  Warburton (1969), Chang (1979) 
                ---> Vỏ kim loại đơn giản, phương pháp Rayleigh-Ritz/FEM cổ điển.
-----------------------------------------------------------------------------------
1990s - 2000s:  Lakis & Laveau (1991, 1997), Amabili (1996, 2002), Damatty (2005)
                ---> Vỏ dị hướng, FSI giải tích và bán giải tích cho vỏ đơn.
-----------------------------------------------------------------------------------
2000s - 2010s:  Kulla (2003), Casimir (1996), Nguyễn Mạnh Cường & Casimir (2007)
                ---> Phương pháp phần tử liên tục (PTLT) cho kết cấu dầm, tấm và vỏ.
-----------------------------------------------------------------------------------
Luận án này:    Vũ Quốc Hiến, GS. Trần Ích Thịnh, PGS. Nguyễn Mạnh Cường
                ---> PTLT cho hệ vỏ liên hợp composite đa phân đoạn (nón-trụ-nón, 
                     trụ bậc) FSI + Thực nghiệm kiểm chứng đa mức chất lỏng.
===================================================================================

Đối với vỏ vật liệu composite dị hướng nhiều lớp, Lakis và Laveau (1991, 1997) cùng Toorani và Lakis (2001) đã đặt nền móng giải tích cho vỏ trụ composite hở chứa dòng chất lỏng chảy, có xét đến ảnh hưởng của biến dạng cắt ngang. Tiếp đó, Kochupillai, Ganesan và Padmanabhan (2002) cũng như Kadoli và Ganesan (2003) đã phát triển mô hình PTHH bán giải tích khảo sát ổn định nhiệt và dao động của vỏ composite tương tác với chất lỏng nhiệt độ cao. Tuy nhiên, các công trình này chủ yếu dừng lại ở các phân tử vỏ đơn lẻ hoặc dạng hình học hình trụ trơn. Đối với kết cấu ghép nối phức tạp, Qu et al. (2013) dùng phương pháp biến phân phân miền cho vỏ kim loại nón - trụ - vành tăng cứng; Caresta và Kessissoglou (2010) kết hợp chuỗi sóng và chuỗi bậc cao cho vỏ nón - trụ đẳng hướng; Kouchakzadeh và Shakouri (2014) khảo sát dao động vỏ nón - nón composite dựa trên lý thuyết Donnell và nguyên lý Hamilton nhưng chỉ dừng lại ở bài toán không chứa chất lỏng.

Trong bối cảnh đó, các trường phái tiếp cận tại Việt Nam do nhóm nghiên cứu của GS. Trần Ích Thịnh, PGS. Nguyễn Mạnh Cường (2013, 2016), Đỗ Văn Hiến và Nguyễn Xuân Hùng (2004), Nguyễn Tiến Khiêm và Trần Văn Liên (2004), cùng luận án của Tạ Thị Hiền (2014) đã từng bước thiết lập phương pháp ma trận độ cứng động lực cho dầm, tấm và các phần tử vỏ đơn. Định vị học thuật của luận án Vũ Quốc Hiến tạo ra bước nhảy vọt (paradigm leap) khi:

  1. Mở rộng thành công phương pháp PTLT giải quyết bài toán FSI cho các hệ vỏ ghép nối liên hợp đa diện phức hợp: trụ bậc, nón - trụ, nón - trụ - nón, nón - nón - nón.
  2. Vượt qua giới hạn của phương pháp giải tích cổ điển (vốn bất khả thi khi giải hệ phương trình vi phân biến thiên hệ số cho vỏ liên hợp) và khắc phục nhược điểm phụ thuộc mật độ lưới cũng như dung lượng bộ nhớ của PTHH truyền thống.
  3. Cung cấp bộ dữ liệu thực nghiệm FSI đa phân đoạn chưa từng được công bố trên các diễn đàn cơ học quốc tế cho kết cấu composite ghép nối nón - trụ và nón - trụ - nón chứa nước ở nhiều mức chất lỏng khác nhau.

Đó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 lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin và lý thuyết thủy động lực học cổ điển của Bessel - Bernoulli sang bài toán dao động tương tác FSI cho vỏ composite dị hướng ghép nối. Bằng cách thiết lập quan hệ chuyển vị 5 thành phần: $$u(x,\theta,z,t) = u_0(x,\theta,t) + z \phi_x(x,\theta,t)$$ $$v(x,\theta,z,t) = v_0(x,\theta,t) + z \phi_\theta(x,\theta,t)$$ $$w(x,\theta,z,t) = w_0(x,\theta,t)$$ nghiên cứu tính đến đầy đủ các góc xoay pháp tuyến $\phi_x, \phi_\theta$ và các hiệu ứng cắt ngang $Q_x, Q_\theta$, loại bỏ giả thiết Kirchhoff-Love vốn gây sai số lớn đối với vật liệu composite nhiều lớp có mô đun đàn hồi trượt ngang thấp ($G_{13}, G_{23} \ll E_1$).

                KHUNG QUAN HỆ CƠ HỌC VÀ MA TRẬN ĐỘ CỨNG N LỚP COMPOSITE
  +-----------------------------------------------------------------------------------+
  | {N} = [A]{ε} + [B]{κ}  --> Ma trận [A]: Độ cứng màng (Extensional stiffness)      |
  | {M} = [B]{ε} + [D]{κ}  --> Ma trận [B]: Tương tác màng-uốn-xoắn (Coupling)        |
  | {Q} = [F]{γ}           --> Ma trận [D]: Độ cứng uốn (Bending stiffness)           |
  |                            Ma trận [F]: Độ cứng cắt ngang (Transverse shear)      |
  +-----------------------------------------------------------------------------------+
                                            |
                                            v
  +-----------------------------------------------------------------------------------+
  |                     TƯƠNG TÁC THỦY ĐỘNG (FLUID-STRUCTURE INTERACTION)             |
  |  ∇²Φ = 0  -->  P(x,θ,t) = -ρ_f (∂Φ/∂t)|_r=R = -m_h (∂²w_0/∂t²)                    |
  |  Khối lượng quy dẫn chất lỏng: m_h = - (ρ_f / k_n) * [I_m(k_n R) / I_m'(k_n R)]   |
  +-----------------------------------------------------------------------------------+

Hệ phương trình chuyển động của vỏ trụ composite chứa chất lỏng được xây dựng với số hạng áp suất thủy động $P(x,\theta,t)$ tác dụng trực tiếp lên phương trình cân bằng lực cắt và chuyển vị pháp tuyến: $$\frac{\partial N_x}{\partial x} + \frac{1}{R}\frac{\partial N_{x\theta}}{\partial \theta} = I_0 \ddot{u}0 + I_1 \ddot{\phi}x$$ $$\frac{\partial N{x\theta}}{\partial x} + \frac{1}{R}\frac{\partial N\theta}{\partial \theta} + \frac{Q_\theta}{R} = I_0 \ddot{v}0 + I_1 \ddot{\phi}\theta$$ $$\frac{\partial Q_x}{\partial x} + \frac{1}{R}\frac{\partial Q_\theta}{\partial \theta} - \frac{N_\theta}{R} - P = I_0 \ddot{w}0$$ $$\frac{\partial M_x}{\partial x} + \frac{1}{R}\frac{\partial M{x\theta}}{\partial \theta} - Q_x = I_1 \ddot{u}0 + I_2 \ddot{\phi}x$$ $$\frac{\partial M{x\theta}}{\partial x} + \frac{1}{R}\frac{\partial M\theta}{\partial \theta} - Q_\theta = I_1 \ddot{v}0 + I_2 \ddot{\phi}\theta$$

Trong đó các hệ số quán tính $I_i = \sum_{k=1}^N \int_{z_{k-1}}^{z_k} \rho^{(k)} z^i dz$ ($i=0,1,2$). Khối lượng chất lỏng tác dụng lên thành vỏ được quy dẫn thông qua hàm thế vận tốc $\Phi(r,\theta,x,t) = W(r)f(x)\cos(m\theta)e^{j\omega t}$ thỏa mãn phương trình Laplace $\nabla^2 \Phi = 0$ và nghiệm hàm Bessel hiệu chỉnh loại 1 bậc $m$ ($I_m(k_n r)$), từ đó rút ra biểu thức khối lượng thủy động bổ sung $m_h$: $$m_h = -\frac{\rho_f}{k_n \frac{I_m'(k_n R)}{I_m(k_n R)}} = -\frac{\rho_f R}{m + k_n R \frac{I_{m+1}(k_n R)}{I_m(k_n R)}}$$

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

Khung phân tích của luận án kết hợp đồng thời ba trụ cột lý thuyết: Cơ học vật rắn biến dạng phi đẳng hướng (Anisotropic continuum mechanics), Thủy động lực học trường thế (Potential flow theory) và Lý thuyết ma trận độ cứng động lực dạng hàm mũ vi phân (Differential dynamic stiffness matrix formulation).

MA TRẬN TRUYỀN VÀ GHÉP NỐI ĐỘ CỨNG ĐỘNG LỰC HỆ VỎ LIÊN HỢP
===================================================================================
1. Phương trình trạng thái:
   (d/dx){y_m(x)} = [A_m(ω)] {y_m(x)}
   với {y_m} = {u_0, v_0, w_0, φ_x, φ_θ, N_x, N_xθ, Q_x, M_x, M_xθ}^T (10 thành phần)

2. Ma trận truyền qua hàm e-mũ (Casimir, 1996):
   {y_m(L)} = [T(L, ω)] {y_m(0)} = exp([A_m] * L) {y_m(0)}

3. Phân khối ma trận truyền và chuyển đổi sang ma trận độ cứng động lực:
   [T] = [ [T11, T12], [T21, T22] ]
   [K(ω)]_e = [ [-T12⁻¹ * T11,  T12⁻¹], 
                [T21 - T22 * T12⁻¹ * T11,  T22 * T12⁻¹] ]

4. Ghép nối tại các mặt cắt liên tục (Zhang & Xiang, 2007):
   {U_i^(j)} = {U_i^(j+1)},  {F_i^(j)} + {F_i^(j+1)} = {0}
   --> [K_global(ω)] * {U_global} = {F_kích_thích}
===================================================================================

Tính độc đáo của khung phân tích nằm ở việc vector hóa trạng thái cơ học thành vector 10 thành phần ${y_m(x)}$, biến hệ phương trình đạo hàm riêng liên kết phức tạp thành phương trình vi phân thường bậc nhất thuần nhất theo biến tọa độ đường sinh $x$: $$\frac{d{y_m(x)}}{dx} = [A_m(\omega)] {y_m(x)}$$

Ma trận vi phân $[A_m(\omega)]$ kích thước $10 \times 10$ được xác lập tường minh từ các hằng số đàn hồi thu gọn $Q_{ij}$, ma trận độ cứng màng $[A]$, mô men uốn $[D]$, tương tác $[B]$, độ cứng cắt $[F]$ và hệ số thủy động $m_h$. Điều kiện biên và điều kiện ghép nối liên tục tại các vị trí phân bậc hoặc thay đổi hình học (như tại giao diện nón - trụ) được thỏa mãn tuyệt đối thông qua việc cân bằng véc-tơ chuyển vị ${u_0, v_0, w_0, \phi_x, \phi_\theta}^T$ và véc-tơ nội lực ${N_x, N_{x\theta}, Q_x, M_x, M_{x\theta}}^T$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái duy thực phê phán và thực nghiệm chuẩn xác (positivist/critical realism quantitative paradigm), kết hợp chặt chẽ giữa mô hình hóa toán học giải tích - số giải thuật chính xác và kiểm chứng thực nghiệm vật lý trong phòng thí nghiệm đạt chuẩn quốc gia.

+-----------------------------------------------------------------------------------+
|                           SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU                            |
+-----------------------------------------------------------------------------------+
| [Mô hình hóa giải tích PTLT]              [Thực nghiệm dao động Modal Testing]   |
|   - Thiết lập vector trạng thái 10x10       - Mẫu GFRP: CT-14, CT-27, CTC-9       |
|   - Tích phân hàm mũ ma trận exp(Am*L)      - Kích thích búa gõ xung lực          |
|   - Lập trình thuật toán MATLAB             - Cảm biến gia tốc PCB Piezotronics   |
|   - Quét đường cong đáp ứng w(ω)            - Bộ thu thập Dactron & RT Pro        |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                  SO SÁNH VÀ HIỆU CHỈNH SAI SỐ LIÊN HỢP (PTLT vs EXP)              |
|          Độ lệch tần số: < 3.8% ở trạng thái khô; < 5.6% khi chứa đầy nước        |
+-----------------------------------------------------------------------------------+

Thiết kế đa cấp độ (multi-level design) bao gồm:

  1. Cấp độ phần tử đơn (Single Element Level): Phần tử vỏ trụ composite khô/chứa chất lỏng và vỏ nón composite khô/chứa chất lỏng.
  2. Cấp độ kết cấu liên hợp (Substructure Assembly Level): Ghép nối các ma trận độ cứng động lực cấp phần tử thành ma trận toàn cục $[K_{global}(\omega)]$ kích thước mở rộng ($15 \times 15, 20 \times 20, \dots$) tương ứng với số đoạn phân bậc và sự phân chia miền có/không chứa chất lỏng.
  3. Cấp độ thực nghiệm vật lý (Physical Validation Level): Đo đạc đáp ứng tần số trên các kết cấu chế tạo thực tế với các mức chất lỏng $H/L \in [0, 1]$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm 4 giai đoạn khép kín và kiểm soát sai số nghiêm ngặt:

TIẾN TRÌNH THỰC HIỆN THÍ NGHIỆM VÀ ĐO ĐẠC MODAL
===================================================================================
BƯỚC 1: Chế tạo mẫu composite (GFRP) bằng công nghệ cuốn/hút chân không tại 
        Viện Nghiên cứu & Chế tạo Tàu thủy - ĐH Thủy sản Nha Trang & ĐHBK Hà Nội.
        Vật liệu: Sợi thủy tinh E-glass / Nhựa polyester không no [0°/90°/0°/90°].
-----------------------------------------------------------------------------------
BƯỚC 2: Gá đặt mẫu trên khung giá cố định, thiết lập điều kiện biên ngàm (Clamped) 
        hoặc tự do (Free) tại Phòng thí nghiệm kiểm soát rung và ồn - Viện Cơ học VN.
-----------------------------------------------------------------------------------
BƯỚC 3: Bố trí các đầu đo gia tốc áp điện đa kênh dọc theo đường sinh và chu vi vỏ; 
        kết nối hệ thống khuếch đại tín hiệu và phân tích phổ FFT Dactron Focus.
-----------------------------------------------------------------------------------
BƯỚC 4: Kích thích dao động bằng búa tác động xung lực (Impulse Hammer) có gắn cảm 
        biến lực; ghi nhận hàm truyền tần số FRF (Frequency Response Function) 
        ở các mức chất lỏng H = 0, 1/4, 1/2, 3/4 và đầy nước (H = L).
===================================================================================

Tam giác đạc phương pháp (Methodological triangulation) được thiết lập thông qua việc so sánh đối chéo kết quả giữa ba nguồn độc lập:

  • Lời giải tính toán từ chương trình lập trình MATLAB của luận án dựa trên PTLT.
  • Lời giải mô phỏng số bằng phương pháp PTHH (ANSYS / ABAQUS) từ các nghiên cứu quốc tế đã công bố.
  • Dữ liệu đo đạc thực nghiệm trực tiếp tại Viện Cơ học Việt Nam.

Data và phân tích

Các thông số cơ lý tính của vật liệu composite sợi thủy tinh/polyester được đo đạc và xác định chuẩn xác:

  • Mô đun đàn hồi dọc sợi $E_1 = 18.5 \text{ GPa}$, mô đun ngang $E_2 = 5.6 \text{ GPa}$.
  • Mô đun trượt $G_{12} = G_{13} = 2.4 \text{ GPa}$, $G_{23} = 1.8 \text{ GPa}$.
  • Hệ số Poisson $\nu_{12} = 0.28, \nu_{21} = 0.085$.
  • Khối lượng riêng composite $\rho = 1850 \text{ kg/m}^3$; Khối lượng riêng của nước $\rho_{nước} = 1000 \text{ kg/m}^3$, của rượu $\rho_{rượu} = 789 \text{ kg/m}^3$.

Phương pháp tìm nghiệm tần số riêng áp dụng kỹ thuật đường cong đáp ứng động lực học (Frequency Response Curve): Đặt một xung lực kích thích điều hòa đơn vị ${Q}m = \sum{m=1}^\infty 1 \cdot \cos(m\theta)$ tại biên kết cấu ($x=L$). Giải hệ phương trình đại số tuyến tính $[K(\omega)]{U} = {F}$ trên miền quét tần số $\omega$ với bước quét mịn $\Delta \omega = 0.1 \text{ rad/s}$. Tại các điểm cực trị (resonance peaks) nơi chuyển vị hướng kính $w(\omega) \to \infty$, hoành độ chính xác tương ứng là tần số dao động tự do $\omega_{mn}$ của kết cấu.

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

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

+-----------------------------------------------------------------------------------+
|           TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ VỀ ĐẶC TRƯNG DAO ĐỘNG FSI                |
+-----------------------------------------------------------------------------------+
| 1. HIỆU ỨNG TẢI KHỐI LƯỢNG NƯỚC (HYDRODYNAMIC MASS EFFECT):                       |
|    - Khi chứa đầy nước (H/L = 1), tần số riêng cơ bản của vỏ composite giảm mạnh  |
|      từ 42% đến 65% so với trạng thái khô (H/L = 0).                              |
|    - Tốc độ suy giảm tần số diễn ra phi tuyến mạnh nhất ở giai đoạn H/L từ 0 đến 0.5|
|                                                                                   |
| 2. HIỆN TƯỢNG ĐẢO MODE DAO ĐỘNG RIÊNG THẤP NHẤT (MODE SWITCHING):                 |
|    - Đối với vỏ nón - trụ và nón - trụ - nón, tần số nhỏ nhất không xuất hiện ở   |
|      mode vòng m = 1 mà xuất hiện ở các mode sóng bậc cao m = 2, 3 hoặc 4         |
|      tùy thuộc vào tỷ lệ kích thước L/R và góc nghiêng nón α.                     |
|                                                                                   |
| 3. TÁC ĐỘNG CỦA HÌNH HỌC NÓN PHÂN KỲ VÀ HỘI TỤ:                                   |
|    - Kết cấu nón phân kỳ (nước chảy mở rộng) làm tăng độ cứng động lực cục bộ     |
|      và nâng cao tần số riêng so với kết cấu nón hội tụ cùng thể tích chứa.       |
|                                                                                   |
| 4. TÍNH DỊ HƯỚNG VÀ CẤU HÌNH XẾP LỚP:                                             |
|    - Cấu hình lớp đối xứng [0°/90°/0°/90°]_s triệt tiêu ma trận tương tác [B],    |
|      tạo tần số dao động cao hơn 12-18% so với cấu hình phi đối xứng cùng bề dày. |
+-----------------------------------------------------------------------------------+

  1. Hiệu ứng khối lượng nước gia tăng gây suy giảm tần số phi tuyến sâu sắc: Khi mức chất lỏng trong vỏ tăng dần từ $H = 0$ đến $H = L$, tần số dao động riêng của mọi dạng vỏ (trụ bậc, nón - trụ, nón - trụ - nón) đều suy giảm nghiêm trọng. Ví dụ với mẫu vỏ nón - trụ CT-27, tần số dao động tự do ở mode cơ bản giảm từ $184.5 \text{ Hz}$ (trạng thái khô) xuống còn $78.2 \text{ Hz}$ (chứa đầy nước), tương ứng mức sụt giảm $57.6%$. Tốc độ suy giảm đặc biệt dốc trong khoảng mức chất lỏng $H/L \in [0.25, 0.75]$.

  2. Quy luật phân bố mode vòng $m$ và mode dọc $n$ khác biệt hoàn toàn với kết cấu dầm/tấm: Trái ngược với kết cấu dầm hoặc tấm phẳng (nơi tần số thấp nhất luôn gắn với mode bậc nhất), ở vỏ composite tròn xoay mỏng, tần số dao động tự do nhỏ nhất ($\omega_{min}$) thường ứng với mode sóng vòng $m = 2, 3, 4$ hoặc $5$. Sự cạnh tranh giữa năng lượng biến dạng màng (chiếm ưu thế ở $m$ nhỏ) và năng lượng biến dạng uốn (chiếm ưu thế ở $m$ lớn) quyết định điểm cực tiểu của phổ tần số.

  3. Ảnh hưởng chi phối của góc nghiêng nón ($\alpha$) và hình học ghép nối: Khi góc nón $\alpha$ tăng từ $15^\circ$ đến $60^\circ$, độ cứng tổng thể của vỏ nón giảm theo phương dọc trục nhưng lại tăng diện tích tiếp xúc thủy động, dẫn đến tần số riêng giảm từ $15%$ đến $30%$. Tại các mặt cắt ghép nối nón - trụ hoặc trụ bậc, xuất hiện sự tập trung ứng suất động và bước nhảy về dạng dao động cục bộ (local mode localization).

  4. Sự tương thích vượt trội giữa lý thuyết PTLT và thực nghiệm: Các kết quả đo thực nghiệm trên mẫu CT-14, CT-27 và CTC-9 tại Viện Cơ học Việt Nam cho thấy độ trùng khớp tuyệt vời với đường cong đáp ứng thu được từ PTLT. Sai số tương đối giữa tính toán lý thuyết và thực nghiệm dao động trong biên độ rất nhỏ, từ $1.8%$ đến $5.6%$, khẳng định độ tin cậy tuyệt đối của phương pháp ma trận độ cứng động lực.

SO SÁNH TẦN SỐ DAO ĐỘNG TỰ DO (Hz) GIỮA PTLT VÀ THỰC NGHIỆM (MẪU CT-27)
+-------------------+--------------------+--------------------+--------------------+
| Mức nước trong vỏ |  Lý thuyết (PTLT)  |  Thực nghiệm (EXP) |    Sai số (%)      |
+-------------------+--------------------+--------------------+--------------------+
| H = 0 (Khô)       |      184.5 Hz      |      179.2 Hz      |       2.87%        |
| H = 1/4 L         |      152.1 Hz      |      147.5 Hz      |       3.02%        |
| H = 1/2 L         |      118.4 Hz      |      114.1 Hz      |       3.63%        |
| H = 3/4 L         |       92.7 Hz      |       88.9 Hz      |       4.10%        |
| H = L (Đầy nước)  |       78.2 Hz      |       74.5 Hz      |       4.73%        |
+-------------------+--------------------+--------------------+--------------------+

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án đã hoàn thiện cơ sở toán học cho phương pháp PTLT trong bài toán FSI vỏ composite đa phân đoạn, mở ra hướng tiếp cận chính xác cho các bài toán phân tích động lực học kết cấu đàn hồi tiếp xúc môi trường truyền sóng.
  • Về mặt phương pháp luận tính toán: Cung cấp giải pháp thay thế tối ưu cho các phần mềm thương mại đắt tiền. Thuật toán ma trận độ cứng động lực lập trình trên MATLAB có tốc độ giải nhanh hơn gấp 10-20 lần so với mô hình 3D Solid-Fluid FEM thông thường mà không đòi hỏi tài nguyên tính toán siêu máy tính.
  • Về mặt kỹ thuật thực tiễn và thiết kế kháng chấn: Cung cấp các biểu đồ thiết kế chuẩn (design charts) cho phép các kỹ sư hàng không, tàu thủy và hóa chất dự báo chính xác vùng tần số nguy hiểm, ngăn chặn triệt để thảm họa cộng hưởng dao động (resonance destruction) và hiện tượng dao động dập dềnh (sloshing-vibration coupling) trong các bồn chứa nhiên liệu lỏng khi vận hành hoặc khi chịu tải trọng động đất.

Limitations và Future Research

Mặc dù đạt được những đột phá khoa học xuất sắc, tác giả cũng thẳng thắn chỉ rõ các giới hạn biên (boundary limitations) của công trình:

  1. Mô hình chất lỏng lý tưởng: Chất lỏng được giả thiết là không nén được, phi nhớt và dòng chảy không xoáy ($\nabla \times \mathbf{v} = 0$). Các hiệu ứng nhớt, lớp biên thủy động sát thành vỏ và tổn hao năng lượng do ma sát chất lỏng chưa được tính đến.
  2. Dao động sóng mặt tự do (Free surface sloshing): Luận án bỏ qua dao động sóng mặt tự do và áp suất trọng trường tĩnh tác động phi tuyến lên chuyển vị ban đầu của vỏ.
  3. Giới hạn đàn hồi tuyến tính và góc nghiêng nón vừa phải: Chưa xét đến phi tuyến hình học (geometric nonlinearity) khi biên độ dao động lớn ($w \ge h$) và biến dạng trượt phi tuyến bậc cao qua chiều dày đối với các kết cấu composite siêu dày.

Chương trình nghiên cứu 5-10 năm tiếp theo (Future Research Agenda):

  • Mở rộng phương pháp PTLT cho vỏ composite chức năng biến đổi cơ tính (Functionally Graded Materials - FGM) và composite gia cường ống nano carbon (CNTRC) chứa chất lỏng chuyển động có áp lực dòng chảy.
  • Tích hợp mô hình chất lỏng nhớt Navier-Stokes và xét đến hiệu ứng tương tác sóng mặt tự do (sloshing coupling) chịu gia tốc động đất đa chiều.
  • Phát triển phần tử liên tục cho vỏ tròn xoay composite chịu tải trọng nhiệt độ cực hạn kết hợp áp suất âm động lực học.

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

+-----------------------------------------------------------------------------------+
|                        MA TRẬN TÁC ĐỘNG CỦA LUẬN ÁN                               |
+-----------------------------------------------------------------------------------+
| 1. GIÁO DỤC & HỌC THUẬT:                                                          |
|    - Cung cấp giáo trình chuyên đề Cơ học kết cấu vỏ & Composite tại ĐHBK Hà Nội. |
|    - Tiềm năng trích dẫn cao trong các tạp chí quốc tế hàng đầu (Q1/Q2 ISI/Scopus)|
|      như Composite Structures, Thin-Walled Structures, Journal of Fluids & Struct.|
|                                                                                   |
| 2. CÔNG NGHIỆP TÀU THỦY & DẦU KHÍ:                                                |
|    - Ứng dụng thiết kế bồn chứa hóa chất tại Nhà máy Hóa chất Việt Trì.           |
|    - Thiết kế tháp phản ứng Amon Nitrat Thái Bình và tháp composite lọc dầu.      |
|                                                                                   |
| 3. HÀNG KHÔNG VŨ TRỤ & QUỐC PHÒNG:                                                |
|    - Tính toán vỏ composite nhẹ cho thùng nhiên liệu tên lửa, vòm che radar       |
|      máy bay và thân tàu ngầm mini composite.                                     |
+-----------------------------------------------------------------------------------+

Công trình nghiên cứu tạo ra tác động sâu rộng đối với ngành cơ học tính toán tại Việt Nam. Các thuật toán và mô hình số của luận án đã được chuyển giao và ứng dụng thực tế trong việc thẩm tra an toàn động lực học cho các kết cấu bồn bể composite chứa hóa chất công nghiệp tại Nhà máy Hóa chất Việt Trì và tổ hợp sản xuất Amon Nitrat Thái Bình. Về mặt quốc tế, phương pháp PTLT giải FSI vỏ liên hợp composite định vị Việt Nam là một trong những nhóm nghiên cứu tiên phong trong khu vực Đông Nam Á làm chủ công nghệ mô phỏng giải tích - số kết cấu composite tương tác chất lỏng bậc cao.

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

+--------------------------+--------------------------------------------------------+
| Nhóm đối tượng           | Giá trị thụ hưởng cụ thể                               |
+--------------------------+--------------------------------------------------------+
| Nghiên cứu sinh &        | Khung phương pháp luận ma trận vi phân trạng thái 10x10|
| Học viên Cao học         | và thuật toán quét phổ tần số trên MATLAB.             |
+--------------------------+--------------------------------------------------------+
| Giáo sư & Các nhà        | Cơ sở lý thuyết vững chắc để phát triển các bài toán   |
| Cơ học tính toán         | tương tác đa trường (Multiphysics: Thermo-Fluid-Solid).|
+--------------------------+--------------------------------------------------------+
| Kỹ sư R&D Hàng không,    | Bộ công cụ tính toán nhanh để tối ưu hóa góc xếp lớp   |
| Đóng tàu & Dầu khí       | và tránh dải tần số cộng hưởng nguy hiểm.              |
+--------------------------+--------------------------------------------------------+
| Cơ quan Quản lý &        | Tiêu chuẩn định lượng an toàn dao động và kiểm định kết|
| Đăng kiểm Công trình     | cấu composite chứa chất lỏng dễ cháy nổ, độc hại.      |
+--------------------------+--------------------------------------------------------+

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công Ma trận độ cứng động lực chính xác $[K(\omega)]$ cho các kết cấu vỏ tròn xoay liên hợp composite đa phân đoạn chứa chất lỏng dựa trên việc mở rộng Lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin kết hợp hàm thế thủy động Bessel loại 1 bậc $m$ ($I_m(k_n r)$). Khác với các mô hình dầm hay tấm của Clough & Penzien (1975) hay Casimir (1996), luận án đã vector hóa 10 đại lượng lực và chuyển vị của vỏ mỏng/dày vừa dị hướng thành hệ phương trình trạng thái vi phân liên tục $\frac{d\mathbf{y}_m}{dx} = \mathbf{A}_m(\omega)\mathbf{y}_m$, cho phép tích phân trực tiếp thông qua hàm mũ ma trận $e^{\mathbf{A}_m x}$.

2. Phương pháp PTLT có bước đột phá nào khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?

  • So với nghiên cứu của Brenneman & Yang (1992) (dùng ma trận độ cứng động lực kết hợp PTHH cho vỏ kim loại đơn) và Qu et al. (2013) (dùng phương pháp biến phân phân đoạn): PTLT trong luận án này xử lý trực tiếp tính chất dị hướng phức tạp của composite nhiều lớp với ma trận tương tác màng - uốn $[B] \neq 0$ mà không cần rời rạc hóa bề mặt thành hàng chục ngàn phần tử như FEM, triệt tiêu hoàn toàn hiện tượng khóa cắt (shear locking).
  • So với Kouchakzadeh & Shakouri (2014) (chỉ giải vỏ nón - nón khô bằng lý thuyết Donnell cổ điển): Luận án giải quyết trọn vẹn bài toán FSI chứa chất lỏng đa mức cho 4 dạng vỏ liên hợp phức tạp (trụ bậc, nón - trụ, nón - trụ - nón, nón - nón - nón) bằng lý thuyết Mindlin chính xác hơn nhiều so với lý thuyết vỏ mỏng Donnell.

3. Phát hiện bất ngờ nhất (counter-intuitive finding) từ dữ liệu phân tích là gì?

Phát hiện bất ngờ nhất là: Tần số dao động tự do nhỏ nhất ($\omega_{min}$) của vỏ nón - trụ và nón - trụ - nón composite chứa nước không nằm ở mode sóng vòng thấp nhất ($m=1$) mà lại xuất hiện ở các mode sóng bậc cao ($m = 3$ hoặc $m = 4$). Điều này trái ngược với trực giác kỹ thuật thông thường ở kết cấu thanh dầm. Nguyên nhân cơ học là do ở $m=1$, năng lượng biến dạng màng (extensional strain energy) của vỏ tròn xoay chiếm tỷ trọng áp đảo làm kết cấu có độ cứng tương đương rất cao; khi $m$ tăng lên $3$ hoặc $4$, năng lượng màng giảm nhanh trong khi năng lượng uốn chưa tăng đủ lớn, tạo ra một "thung lũng tần số riêng" cực tiểu trước khi tần số tăng trở lại ở các mode $m \ge 5$ do năng lượng uốn chi phối.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) được cung cấp như thế nào?

Luận án cung cấp quy trình giải thuật minh bạch từng bước:

  1. Nhập ma trận hằng số đàn hồi vật liệu $E_1, E_2, G_{12}, \nu_{12}$ và tính ma trận độ cứng $A_{ij}, B_{ij}, D_{ij}, F_{ij}$.
  2. Thiết lập ma trận trạng thái $\mathbf{A}_m(\omega)$ kích thước $10 \times 10$ theo công thức tường minh (2.26).
  3. Tính ma trận truyền $\mathbf{T}(L_e, \omega) = \exp(\mathbf{A}_m L_e)$ bằng chuỗi hàm mũ ma trận Casimir.
  4. Phân khối ma trận $\mathbf{T}$ thành các khối $5 \times 5$ và chuyển đổi thành ma trận độ cứng phần tử $\mathbf{K}_e(\omega)$.
  5. Ghép nối theo điều kiện liên tục (2.36) thành $\mathbf{K}_{global}(\omega)$.
  6. Đặt lực kích thích đơn vị, giải tìm chuyển vị pháp tuyến $w(\omega)$ trên MATLAB để vẽ đường cong đáp ứng và trích xuất đỉnh cộng hưởng.

5. Kế hoạch nghiên cứu 10 năm được phác thảo ra sao?

Kế hoạch 10 năm tập trung vào 3 trọng tâm:

  1. Giai đoạn 1 (1-3 năm): Phát triển phần mềm tính toán đóng gói chuyên dụng PTLT-Shell-FSI có giao diện đồ họa (GUI) phục vụ các viện thiết kế tàu thủy và công trình biển.
  2. Giai đoạn 2 (4-6 năm): Mở rộng mô hình sang kết cấu vật liệu thông minh (Smart Structures) tích hợp các lớp áp điện (piezoelectric layers) để chủ động kiểm soát và giảm chấn dao động vỏ composite chứa chất lỏng.
  3. Giai đoạn 3 (7-10 năm): Tích hợp bài toán phá hủy mỏi, rão và tương tác nhiệt - thủy - cơ phi tuyến (Thermo-Hydro-Mechanical damage mechanics) cho vỏ composite hàng không vũ trụ tái sử dụng.

Kết luận

Luận án tiến sĩ của Vũ Quốc Hiến là một công trình khoa học mẫu mực, xuất sắc và toàn diện trong chuyên ngành Cơ học vật liệu và kết cấu. Sáu đóng góp cốt lõi định hình giá trị vượt thời gian của luận án bao gồm:

  1. Thiết lập hoàn chỉnh cơ sở lý thuyết và thuật toán Phương pháp Phần tử Liên tục (PTLT) dựa trên lý thuyết Reissner-Mindlin để phân tích dao động tự do của các hệ vỏ tròn xoay liên hợp composite lớp trực hướng có mặt cắt thay đổi.
  2. Giải quyết triệt để bài toán Tương tác Thủy động lực học (FSI) thông qua việc tích hợp chính xác hàm thế vận tốc Laplace - Bessel vào ma trận vi phân động lực học $10 \times 10$ của vỏ composite chứa chất lỏng.
  3. Phát triển thành công hệ thống chương trình tính toán số trên nền tảng MATLAB, cung cấp giải pháp xác định tần số riêng qua đường cong đáp ứng với tốc độ hội tụ nhanh, độ chính xác cao và tiết kiệm tài nguyên bộ nhớ vượt trội so với FEM.
  4. Khám phá và giải thích bản chất cơ học của các quy luật động lực học phi tuyến: Hiện tượng suy giảm tần số sâu sắc theo mức chất lỏng $H/L$, hiện tượng chuyển dịch mode dao động nguy hiểm ($m=3,4$) và ảnh hưởng của tính dị hướng màng - uốn $[B]$.
  5. Chế tạo thành công và thực hiện quy trình thực nghiệm đo dao động Modal Testing chính xác cao trên các mẫu vỏ composite nón - trụ (CT-14, CT-27) và nón - trụ - nón (CTC-9) chứa nước, tạo lập bộ dữ liệu thực nghiệm chuẩn mực kiểm chứng lý thuyết với sai số dưới $5.6%$.
  6. Định hình các khuyến nghị kỹ thuật then chốt cho ngành công nghiệp đóng tàu, bồn bể hóa chất, dầu khí và kỹ thuật hàng không vũ trụ tại Việt Nam, đặt nền tảng vững chắc cho các nghiên cứu cơ học kết cấu tiên tiến trong tương lai.