Tổng quan về luận án

Nghiên cứu động lực học kết cấu vỏ tròn xoay bằng vật liệu composite phân lớp tiếp xúc với chất lỏng là một trong những chủ đề then chốt của cơ học tính toán hiện đại, đóng vai trò sống còn trong kỹ thuật hàng không vũ trụ, công nghiệp đóng tàu, kỹ thuật dầu khí và công nghiệp hóa chất hạt nhân. Khi vận hành trong điều kiện thực tế, như tác giả luận án đã nhấn mạnh: "Các đặc tính dao động của các kết cấu composite nói trên khi chứa chất lỏng bị thay đổi nhiều so với điều kiện làm việc trong không khí. Vì vậy, vấn đề nghiên cứu định lượng sự thay đổi của tần số và dạng dao động của các kết cấu vỏ composite lớp tròn xoay chứa chất lỏng sẽ có ý nghĩa khoa học và có vai trò quan trọng trong kỹ thuật, cụ thể là trong tính toán, thiết kế tối ưu các kết cấu."

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: trong khi các bài toán dao động của vỏ kim loại đẳng hướng tiếp xúc với chất lỏng đã được nghiên cứu sâu rộng từ 50–60 năm qua (như công trình của Amabili, 1996; Cho và cộng sự, 2002), thì các nghiên cứu định lượng về vỏ composite tròn xoay ghép nối đa đoạn phức tạp (như vỏ trụ bậc có chiều dày thay đổi, vỏ nón–trụ, nón–trụ–nón, nón–nón–nón) chứa chất lỏng vẫn còn rất hạn chế. Tính dị hướng cao của vật liệu composite phân lớp gây ra hiệu ứng tương tác cơ học màng – uốn – xoắn phức tạp (thể hiện qua ma trận $[B]$), kết hợp với tương tác thủy động lực học phi tuyến tại bề mặt phân giới kết cấu – lưu chất, dẫn đến hiện tượng suy giảm nghiêm trọng tần số dao động tự do và thay đổi hoàn toàn quy luật phân bố mode dao động.

Luận án tập trung giải quyết 5 câu hỏi và giả thuyết nghiên cứu chính:

  1. Câu hỏi 1 (Q1) & Giả thuyết 1 (H1): Làm thế nào để mô hình hóa chính xác dao động tự do của vỏ trụ bậc composite có độ dày thay đổi chứa chất lỏng mà không gặp phải sai số tích lũy do chia lưới rời rạc hóa của phương pháp phần tử hữu hạn (PTHH)? Giả thuyết H1: Phương pháp phần tử liên tục (PTLT) kết hợp hàm mũ ma trận giải tích sẽ loại bỏ hoàn toàn lỗi hội tụ lưới và cho nghiệm tần số chính xác cao.
  2. Câu hỏi 2 (Q2) & Giả thuyết 2 (H2): Hiệu ứng ghép nối hình học tại vị trí tiếp giáp của vỏ nón–trụ, nón–trụ–nón và nón–nón–nón composite ảnh hưởng như thế nào đến ma trận độ cứng động lực tổng thể khi chứa các mức chất lỏng $H$ khác nhau? Giả thuyết H2: Điều kiện liên tục về lực cắt, lực màng, mô men uốn và góc xoay tại mặt cắt tiếp giáp cho phép thiết lập ma trận độ cứng động lực ghép nối trực tiếp mà không cần dùng phần tử chuyển tiếp nhân tạo.
  3. Câu hỏi 3 (Q3) & Giả thuyết 3 (H3): Khối lượng riêng của chất lỏng $\rho_f$ và vật liệu vỏ $\rho_k$, cùng góc mở nón $\alpha$ tương tác ra sao đối với tần số riêng và số sóng mode vòng $m$? Giả thuyết H3: Hiệu ứng khối lượng ảo thủy động làm giảm phi tuyến tần số riêng, đồng thời dịch chuyển mode dao động nguy hiểm nhất sang các bậc sóng vòng $m$ cao hơn.
  4. Câu hỏi 4 (Q4) & Giả thuyết 4 (H4): Cấu hình phân lớp đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_s$ và không đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_2$ có mức độ triệt tiêu tương tác màng–uốn khác nhau như thế nào dưới tải trọng thủy động? Giả thuyết H4: Cấu hình đối xứng triệt tiêu hoàn toàn ma trận $[B]$, tạo độ cứng uốn hiệu dụng cao hơn và kháng dao động tốt hơn trong môi trường chất lỏng.
  5. Câu hỏi 5 (Q5) & Giả thuyết 5 (H5): Dữ liệu đo đạc thực nghiệm từ búa kích động và cảm biến gia tốc trên mẫu vỏ thực tế có kiểm chứng được mô hình giải tích PTLT trong miền sai số kỹ thuật cho phép ($< 8%$) hay không? Giả thuyết H5: Phương pháp đường cong đáp ứng thu nhận từ giải thuật PTLT phản ánh chính xác các đỉnh cộng hưởng đo được trên phổ FFT thực nghiệm.

Khung lý thuyết nền tảng tích hợp Lý thuyết biến dạng cắt bậc nhất Reissner–Mindlin (FSDT), Cơ học vật liệu composite dị hướng nhiều lớp, và Lý thuyết thế lưu chất Laplace–Bernoulli. Phạm vi nghiên cứu bao quát cả phân tích số lý thuyết trên các cấu hình vỏ tròn xoay 2 bậc, 4 bậc, vỏ nón–trụ, nón–trụ–nón, nón–nón–nón với các góc mở nón $\alpha_1, \alpha_2, \alpha_3$ biến thiên từ $15^\circ$ đến $90^\circ$, kết hợp khảo nghiệm thực nghiệm trên 3 bộ mẫu composite sợi thủy tinh/polyester không no (mẫu CT-27, CT-14, CTC-9).

Literature Review và Positioning

Nghiên cứu động lực học vỏ tròn xoay chứa chất lỏng trải qua ba dòng chảy học thuật chính. Dòng chảy thứ nhất tập trung vào vỏ kim loại đẳng hướng: Brenneman và Yang (1992) kết hợp ma trận độ cứng động lực với PTHH; Han và Liu (1994) dùng lý thuyết vỏ mỏng Flügge và phương pháp ma trận truyền; Amabili (1996, 2002) phát triển phương pháp Ritz giải tích phi tuyến biên độ lớn cho vỏ trụ chứa đầy nước; Cho và cộng sự (2002) sử dụng PTHH khảo sát vỏ có vách ngăn ngang chịu tải trọng động đất; Senthil Kumar và Ganesan (2008) phân tích mất ổn định thủy đàn hồi vỏ nón cụt; Qu và cộng sự (2013), Ma và cộng sự (2014) dùng phương pháp biến phân Fourier–Ritz giải vỏ nón–trụ kim loại với điều kiện biên tùy ý.

Dòng chảy thứ hai mở rộng sang kết cấu composite phân lớp: Lakis và Laveau (1991, 1997) phân tích phi tuyến vỏ trụ dị hướng chứa chất lỏng chảy; Toorani và Lakis (2001) khảo sát động lực học vỏ composite lớp hở chứa dòng chảy có xét đến biến dạng cắt; Kadoli và Ganesan (2003) nghiên cứu ổn định vỏ composite trụ tròn tương tác với chất lỏng nóng bằng PTHH bán giải tích; Patel và cộng sự (2000, 2014) phân tích vỏ composite nón–trụ bằng PTHH và mô hình Bert; Kouchakzadeh và Shakouri (2014) dùng lý thuyết vỏ thoải Donnell và nguyên lý Hamilton cho vỏ composite nón–nón. Tại Việt Nam, Trần Ích Thịnh, Nguyễn Mạnh Cường (2013, 2016) và Tạ Thị Hiền (2014) đã tiên phong áp dụng phương pháp PTLT cho tấm và vỏ trụ, nón composite đơn lẻ tiếp xúc lưu chất.

                      TIẾN TRÌNH HỌC THUẬT VỀ DAO ĐỘNG VỎ TRÒN XOAY
┌───────────────────────────────────────────────────────────────────────────────────┐
│ DÒNG CHẢY 1: VỎ KIM LOẠI ĐẲNG HƯỚNG                                               │
│ Flügge, Amabili (1996, 2002), Cho (2002), Senthil Kumar (2008), Qu (2013)         │
└────────────────────────────────────────┬──────────────────────────────────────────┘
                                         ▼
┌───────────────────────────────────────────────────────────────────────────────────┐
│ DÒNG CHẢY 2: VỎ COMPOSITE ĐƠN ĐOẠN / PTHH                                         │
│ Lakis & Laveau (1997), Toorani (2001), Kadoli (2003), Patel (2000), Tạ Thị Hiền    │
└────────────────────────────────────────┬──────────────────────────────────────────┘
                                         ▼
┌───────────────────────────────────────────────────────────────────────────────────┐
│ POSITIONING CỦA LUẬN ÁN: VỎ COMPOSITE GHÉP NỐI ĐA ĐOẠN + PTLT + THỰC NGHIỆM      │
│ Vỏ trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón chứa chất lỏng (FSDT + DSM + EMA)  │
└───────────────────────────────────────────────────────────────────────────────────┘

Tranh luận học thuật lớn tồn tại giữa hai trường phái tiếp cận: trường phái số rời rạc (PTHH truyền thống qua ANSYS, NASTRAN, ABAQUS) và trường phái giải tích/bán giải tích (Dynamic Stiffness Matrix / PTLT của Gorman, 1982; Casimir, 1996; Kulla, 2003). Trường phái PTHH vấp phải hạn chế nghiêm trọng về sai số khi chia lưới ở vùng tần số cao và khối lượng tính toán khổng lồ khi xử lý bài toán tương tác lưu chất – kết cấu đa bậc. Ngược lại, trường phái ma trận độ cứng động lực truyền thống gặp khó khăn khi tìm nghiệm định thức $|K(\omega)| = 0$ bằng thuật toán Wittrick–Williams do xuất hiện các điểm kỳ dị của hàm siêu việt.

Luận án định vị chính xác tại điểm giao thoa đột phá: phát triển phương pháp Phần tử liên tục (PTLT) dựa trên việc thiết lập véc tơ trạng thái 10 thành phần kết hợp nghiệm hàm mũ ma trận $T = \exp\left(\int A_m dx\right)$, giải quyết bài toán trị riêng qua Phương pháp đường cong đáp ứng (Response Curve Method). Luận án so sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu: nghiên cứu của Zhang và Xiang (2007) về vỏ trụ bậc kim loại bằng phương pháp không gian trạng thái và công trình của Efraim và Eisenberger (2006) về vỏ nón–trụ đẳng hướng. Kết quả của luận án mở rộng vượt bậc các mô hình trên khi giải quyết thành công bài toán vỏ composite đa lớp dị hướng cao có ghép nối đa đoạn phức tạp chứa lưu chấ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 sâu sắc Lý thuyết biến dạng cắt bậc nhất Reissner–Mindlin (FSDT) và Lý thuyết tấm vỏ composite cổ điển (CLPT) bằng việc thiết lập hệ phương trình vi phân liên kết 10 ẩn số động lực học, bao quát toàn diện các hiệu ứng tương tác màng ($[A]$), tương tác màng–uốn–xoắn ($[B]$), độ cứng uốn ($[D]$) và độ cứng cắt ngang ($[F]$) với hệ số hiệu chỉnh cắt $f = 5/6$:

$$\begin{cases} {N} = [A]{\varepsilon^0} + [B]{\kappa} \ {M} = [B]{\varepsilon^0} + [D]{\kappa} \ {Q} = f[F]{\gamma} \end{cases}$$

Sự chuyển dịch mô hình (paradigm shift) nằm ở việc thay thế phép xấp xỉ chuyển vị tĩnh từng phần tử trong PTHH bằng trường chuyển vị động chính xác trong miền tần số. Bằng cách biểu diễn véc tơ chuyển vị và nội lực dưới dạng chuỗi Fourier hội tụ theo tọa độ vòng $\theta$:

$$u_0(x,\theta,t) = u_m(x)\cos(m\theta)e^{i\omega t}, \quad v_0(x,\theta,t) = v_m(x)\sin(m\theta)e^{i\omega t}, \quad w_0(x,\theta,t) = w_m(x)\cos(m\theta)e^{i\omega t}$$

Phương trình chuyển động liên kết lưu chất – kết cấu được đưa về hệ phương trình trạng thái vi phân thường bậc nhất thuần nhất:

$$\frac{d{y}_m}{dx} = [A_m]{y}_m$$

Trong đó véc tơ trạng thái bao gồm 10 thành phần: ${y}m = {u_m, v_m, w_m, \varphi{xm}, \varphi_{\theta m}, N_{xm}, N_{x\theta m}, Q_{xm}, M_{xm}, M_{x\theta m}}^T$.

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

Khung phân tích của luận án là sự hợp nhất chặt chẽ giữa 3 lý thuyết nền tảng:

  1. Cơ học kết cấu vỏ Reissner–Mindlin: Xét đến ảnh hưởng của biến dạng trượt ngang và quán tính quay qua các đại lượng $I_0, I_1, I_2 = \sum_{k=1}^N \int_{z_{k-1}}^{z_k} \rho^{(k)} (1, z, z^2) dz$.
  2. Thủy động lực học thế lưu chất lý tưởng: Mô tả chuyển động của chất lỏng không nén được, không nhớt, không xoáy qua phương trình vi phân Laplace trong hệ tọa độ trụ $\nabla^2 \Phi(r,\theta,x,t) = 0$ và phương trình tích phân áp suất Bernoulli $P = -\rho_f \frac{\partial \Phi}{\partial t}\Big|_{r=R}$.
  3. Lý thuyết ma trận độ cứng động lực (DSM): Xây dựng ma trận truyền giải tích $T(L) = \exp([A_m]L)$, sau đó phân khối thành 4 ma trận con liên hệ giữa véc tơ chuyển vị ${U}$ và lực ${F}$ tại hai đầu biên:

$$\begin{Bmatrix} {U(L)} \ {F(L)} \end{Bmatrix} = \begin{bmatrix} [T_{11}] & [T_{12}] \ [T_{21}] & [T_{22}] \end{bmatrix} \begin{Bmatrix} {U(0)} \ {F(0)} \end{Bmatrix} \implies [K(\omega)] = \begin{bmatrix} -[T_{12}]^{-1}[T_{11}] & [T_{12}]^{-1} \ [T_{21}] - [T_{22}][T_{12}]^{-1}[T_{11}] & [T_{22}][T_{12}]^{-1} \end{bmatrix}$$

                KHUNG PHÂN TÍCH TƯƠNG TÁC THỦY ĐỘNG LỰC HỌC - KẾT CẤU
┌────────────────────────────────────────┐    ┌────────────────────────────────────────┐
│   CƠ HỌC VỎ COMPOSITE PHÂN LỚP FSDT    │    │      THỦY ĐỘNG LỰC HỌC LƯU CHẤT        │
│  - Quan hệ ứng suất - biến dạng        │    │  - Phương trình Laplace: ∇²Φ = 0       │
│  - Ma trận [A], [B], [D], [F]          │    │  - Áp suất Bernoulli: P = -ρf ∂Φ/∂t    │
│  - Quán tính quay I₀, I₁, I₂           │    │  - Nghiệm hàm Bessel hiệu chỉnh Iₘ, Kₘ │
└───────────────────┬────────────────────┘    └───────────────────┬────────────────────┘
                    │                                             │
                    └──────────────────────┬──────────────────────┘
                                           ▼
                    ┌─────────────────────────────────────────────┐
                    │    HỆ PHƯƠNG TRÌNH KHÔNG GIAN TRẠNG THÁI    │
                    │         d{y}ₘ/dx = [Aₘ(ω, m, P)]{y}ₘ        │
                    │   Khối lượng chất lỏng ảo: m*(ω, m, kn)     │
                    └──────────────────────┬──────────────────────┘
                                           ▼
                    ┌─────────────────────────────────────────────┐
                    │       MA TRẬN ĐỘ CỨNG ĐỘNG LỰC [K(ω)]       │
                    │     Ma trận truyền: T = exp(∫ Aₘ dx)        │
                    │  Ghép nối liên tục tại mặt cắt tiếp giáp    │
                    └──────────────────────┬──────────────────────┘
                                           ▼
                    ┌─────────────────────────────────────────────┐
                    │       PHƯƠNG PHÁP ĐƯỜNG CONG ĐÁP ỨNG        │
                    │   Kích thích lực đơn vị Qₓ(L) = cos(mθ)/πR  │
                    │    Quét miền tần số -> Xác định tần số ωₙ   │
                    └─────────────────────────────────────────────┘

Điểm biên giới lý thuyết (boundary conditions) được xác định rõ ràng: vật liệu làm việc hoàn toàn trong miền đàn hồi tuyến tính biến dạng bé; áp suất chất lỏng tác dụng vuông góc tuyệt đối với mặt trung bình của vỏ; số hạng áp suất tương tác thủy động được đưa trực tiếp vào phần tử đường chéo $A_{83}$ của ma trận ma trận trạng thái qua hệ số khối lượng ảo $m^*$:

$$A_{83} = \frac{c_6}{R^2} - I_0 \omega^2 - m^* \omega^2, \quad \text{với } m^* = \rho_f \frac{I_m(k_n R)}{k_n I_m'(k_n R)}$$

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 thực chứng toán học kết hợp kiểm chứng thực nghiệm (Positivist Empiricism). Luận án phát triển một kiến trúc tính toán đa tầng từ cấp độ phần tử cơ bản (vỏ trụ, vỏ nón), cấp độ đa đoạn ghép nối (trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón), đến cấp độ toàn thể kết cấu trong tương tác với môi trường chất lỏng đa pha.

                                  KIẾN TRÚC MÔ HÌNH HÓA ĐA CẤP
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ CẤP ĐỘ 1: PHẦN TỬ LIÊN TỤC CƠ BẢN (Ma trận 10x10)                                       │
│ [K_vỏ_trụ_khô]   [K_vỏ_trụ_chất_lỏng]   [K_vỏ_nón_khô]   [K_vỏ_nón_chất_lỏng]           │
└────────────────────────────────────────────┬────────────────────────────────────────────┘
                                             ▼
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ CẤP ĐỘ 2: GHÉP NỐI ĐIỀU KIỆN LIÊN TỤC MẶT CẮT (Ma trận 15x15 -> 30x30)                 │
│ N_xi = N_xj,  Q_xi = Q_xj,  M_xi = M_xj,  u_i = u_j,  w_i = w_j,  φ_xi = φ_xj            │
└────────────────────────────────────────────┬────────────────────────────────────────────┘
                                             ▼
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ CẤP ĐỘ 3: PHÂN TÍCH TẦN SỐ VÀ DẠNG DAO ĐỘNG QUA MATLAB                                  │
│ Quét dải tần số ω -> Giải {F} = [K(ω)]{U} -> Thu nhận đỉnh chuyển vị w_max              │
└────────────────────────────────────────────┬────────────────────────────────────────────┘
                                             ▼
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ CẤP ĐỘ 4: KHẢO NGHIỆM THỰC NGHIỆM KIỂM CHỨNG (EMA)                                     │
│ Mẫu CT-27, CT-14, CTC-9 -> Búa xung lực -> Cảm biến PCB -> Bộ thu FFT -> Đối chiếu PTLT│
└─────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu trải qua các bước chuẩn hóa nghiêm ngặt:

  1. Thiết lập thuật toán giải ma trận truyền: Ứng dụng giải thuật tính hàm mũ ma trận qua chuỗi lũy thừa mở rộng $e^{[X]} = \sum_{k=0}^\infty \frac{[X]^k}{k!}$, chia đoạn vi phân đường sinh để đảm bảo ma trận vi phân $[A_m]$ hội tụ tuyệt đối với độ chính xác số thực kép.

  2. Kỹ thuật ghép nối ma trận độ cứng động lực: Tại các mặt cắt tiếp giáp giữa các bậc trụ hoặc giữa đoạn nón và đoạn trụ, áp dụng điều kiện cân bằng nội lực và tương thích chuyển vị theo Zhang và Xiang (2007):

    $$N_{xi} = N_{xj}, \quad Q_{xi} = Q_{xj}, \quad N_{x\theta i} = N_{x\theta j}, \quad M_{xi} = M_{xj}, \quad M_{x\theta i} = M_{x\theta j}$$ $$u_i = u_j, \quad v_i = v_j, \quad w_i = w_j, \quad \varphi_{xi} = \varphi_{xj}, \quad \varphi_{\theta i} = \varphi_{\theta j}$$

    Ma trận tổng thể được mở rộng từ cỡ $10 \times 10$ cho 1 phần tử lên $15 \times 15$ cho 2 phần tử ghép nối và đạt $30 \times 30$ cho kết cấu vỏ trụ 4 bậc chứa 5 miền phần tử liên tục (gồm các đoạn khô và đoạn chứa chất lỏng).

  3. Quy trình thực nghiệm modal: Chế tạo mẫu thí nghiệm bằng phương pháp bọc khuôn lăn tay từ vải sợi thủy tinh dệt kết hợp nhựa polyester không no có bổ sung chất đóng rắn MEKP. Sử dụng đồ gá chuyên dụng kẹp chặt gốc vỏ để tạo điều kiện biên ngàm (Clamped - C), đầu trên để tự do (Free - F). Hệ thống thu thập tín hiệu gồm búa gõ kích động lực tích hợp cảm biến lực piezo, đầu đo gia tốc áp điện một trục có độ nhạy cao gắn tại các vị trí tối ưu dọc đường sinh, nối vào card thu thập dữ liệu đa kênh và xử lý phổ tần số qua phần mềm phân tích Fourier số.

                              SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM MODAL
┌─────────────────┐       ┌─────────────────┐       ┌─────────────────┐
│   BÚA KÍCH ĐỘNG │──────>│   MẪU VỎ NÓN-TRỤ│<──────│ CẢM BIẾN GIA TỐC│
│   (Force Sensor)│       │   (CT-27 / CT-14)│       │  (Accelerometer)│
└────────┬────────┘       └────────┬────────┘       └────────┬────────┘
         │                         │ Mực nước H              │
         ▼                         ▼                         ▼
┌─────────────────────────────────────────────────────────────────────┐
│                   BỘ THU THẬP TÍN HIỆU ĐA KÊNH                      │
└──────────────────────────────────┬──────────────────────────────────┘
                                   ▼
┌─────────────────────────────────────────────────────────────────────┐
│                 MÁY TÍNH XỬ LÝ PHỔ TẦN SỐ (FFT)                     │
│               So sánh đường cong FRF và kết quả PTLT                │
└─────────────────────────────────────────────────────────────────────┘

Data và phân tích

Toàn bộ thuật toán được lập trình tự động hóa hoàn toàn trong môi trường MATLAB.

  • Thông số vật liệu tính toán: Vỏ composite sợi thủy tinh/polyester có mô đun đàn hồi dọc sợi $E_1 = 38.5 \text{ GPa}$, mô đun ngang $E_2 = 8.4 \text{ GPa}$, mô đun trượt $G_{12} = G_{13} = 4.2 \text{ GPa}$, $G_{23} = 3.1 \text{ GPa}$, hệ số Poisson $\nu_{12} = 0.28$, khối lượng riêng $\rho = 1800 \text{ kg/m}^3$. Chất lỏng khảo sát gồm nước sạch ($\rho_f = 1000 \text{ kg/m}^3$) và cồn/rượu công nghiệp ($\rho_f = 789 \text{ kg/m}^3$).
  • Kiểm định độ tin cậy mô hình số: Đối chiếu tần số không thứ nguyên $\bar{\Omega} = \omega R \sqrt{\rho(1-\nu^2)/E}$ của vỏ trụ 2 bậc kim loại với kết quả công bố của Zhang và Xiang (2007) cho thấy sai số nhỏ hơn $0.05%$. So sánh tần số dao động nhỏ nhất $\bar{\Omega}1 = \omega_1 R_1 \sqrt{\rho h / A{11}}$ của vỏ nón–trụ composite với nghiên cứu của Patel và cộng sự (2000) tại các góc nón $\alpha = 0^\circ, 30^\circ, 60^\circ$ cho sự trùng khớp tuyệt đối về cả tần số và số sóng vòng $m$.

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

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

  1. Quy luật suy giảm tần số phi tuyến theo mực chất lỏng: Khi mức chất lỏng trong vỏ trụ bậc và vỏ nón–trụ tăng từ trạng thái rỗng ($H/L = 0$) đến đầy hoàn toàn ($H/L = 1.0$), tần số dao động cơ bản của kết cấu suy giảm mạnh từ $42%$ đến $68%$ do ảnh hưởng thống trị của khối lượng ảo tham gia dao động $m^*$.
  2. Sự dịch chuyển số sóng mode vòng nguy hiểm ($m$): Đối với vỏ nón–trụ và nón–trụ–nón composite, mode dao động ứng với tần số riêng nhỏ nhất không xuất hiện ở $m=1$ (dạng uốn tổng thể kiểu dầm) mà xuất hiện ở các mode sóng vòng bậc cao ($m = 3, 4, 5, 6$) tùy thuộc vào tỷ lệ hình học $L/R$ và góc mở nón $\alpha$. Khi chứa đầy chất lỏng, đỉnh đáp ứng tần số dịch chuyển sang các mode có $m$ lớn hơn.
  3. Ảnh hưởng vượt trội của cấu hình đối xứng: Trong vỏ trụ 4 bậc composite, cấu hình phân lớp đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_s$ mang lại tần số dao động tự do cao hơn từ $12.5%$ đến $18.3%$ so với cấu hình phản đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]_2$ dưới cùng điều kiện biên ngàm–tự do (C-F) và ngàm–ngàm (C-C), chứng minh tính hiệu quả của việc triệt tiêu độ cứng liên kết màng–uốn $[B] = 0$.
  4. Tác động của góc côn phân kỳ và hội tụ: Vỏ nón–trụ–nón có góc nón phân kỳ làm giảm độ cứng tổng thể nhanh hơn so với cấu hình nón hội tụ. Khi góc mở nón tăng từ $\alpha = 15^\circ$ lên $\alpha = 60^\circ$, tần số dao động riêng nhỏ nhất giảm tới $35.4%$.
  5. Kiểm chứng thực nghiệm xuất sắc: Dữ liệu đo đạc thực nghiệm trên các mẫu CT-27 (nón–trụ góc dốc $27^\circ$), CT-14 (góc dốc $14^\circ$) và CTC-9 (nón–trụ–nón) chứa các mực nước $0%, 25%, 50%, 75%, 100%$ cho thấy sự tương thích rất cao với đường cong đáp ứng lý thuyết PTLT. Sai số giữa thực nghiệm và mô hình tính dao động trong khoảng hẹp từ $2.1%$ đến $7.6%$, khẳng định độ tin cậy tuyệt đối của phương pháp đề xuất.
                              SO SÁNH TẦN SỐ DAO ĐỘNG RIÊNG (Hz)
                    KẾT QUẢ TÍNH PTLT VS THỰC NGHIỆM TRÊN MẪU CT-14
  Tần số (Hz)
   500 ┼─────────────────────────────────────────────────────────────────┐
       │     ▲ PTLT (Lý thuyết)                                          │
   400 ┼─────┼───────────────────▲───────────────────────────────────────┤
       │     ● Thực nghiệm (EMA) │                                       │
   300 ┼─────────────────────────┼───────────────────▲───────────────────┤
       │                         │                   │           ▲       │
   200 ┼─────────────────────────●───────────────────┼───────────┼───────┤
       │                                             ●           ●       │
   100 ┼─────────────────────────────────────────────────────────────────┤
       │                                                                 │
     0 └─────┬───────────────────┬───────────────────┬───────────┬───────┘
            0% (Rỗng)          25% Mực nước         50%         100% (Đầy)

Implications đa chiều

  • Đóng góp lý thuyết: Bổ sung vào kho tàng cơ học vật liệu composite một khung tính toán bán giải tích hoàn chỉnh, chứng minh tính chính xác của lý thuyết FSDT khi tích hợp trường chuyển vị động trong giải quyết bài toán FSI (Fluid-Structure Interaction).
  • Đổi mới phương pháp luận: Khẳng định phương pháp Phần tử liên tục (PTLT) kết hợp nghiệm hàm mũ ma trận và kỹ thuật đường cong đáp ứng vượt trội hoàn toàn so với thuật toán Wittrick–Williams cổ điển về tốc độ tính toán, tính ổn định và khả năng xử lý bài toán ghép nối phức tạp.
  • Ứng dụng thực tiễn: Cung cấp bộ công cụ tính toán chính xác cho các kỹ sư thiết kế bồn bể composite chứa hóa chất (như các công trình tại Nhà máy Hóa chất Việt Trì, Nhà máy Amon Nitrat Thái Bình), tháp làm mát, vòm che radar, thân tàu ngầm mini và bồn chứa nhiên liệu tên lửa, giúp tránh hiện tượng cộng hưởng phá hủy kết cấu khi vận hành.
  • Khuyến nghị chính sách kỹ thuật: Đề xuất ban hành tiêu chuẩn tính toán dao động thủy đàn hồi bắt buộc trong thiết kế chế tạo bồn bể composite chịu áp lực và chứa lưu chất tại Việt Nam.

Limitations và Future Research

Nghiên cứu thể hiện tính trung thực học thuật cao khi chỉ rõ 4 giới hạn nội tại:

  1. Giới hạn trường lưu chất: Chất lỏng mới được mô hình hóa ở trạng thái lý tưởng (không nén được, không nhớt, không xoáy) và dao động ở biên độ bé; chưa xét đến hiệu ứng sóng bề mặt phi tuyến mạnh (sloshing), hiện tượng tạo bọt khí (cavitation) hoặc chất lỏng có độ nhớt cao.
  2. Giới hạn trạng thái biến dạng: Mô hình FSDT giả thiết góc xoay pháp tuyến biến thiên tuyến tính theo chiều dày, chưa bao quát hiệu ứng biến dạng cắt bậc cao (HSDT) khi vỏ có tỷ lệ bề dày trên bán kính lớn ($h/R > 0.1$).
  3. Giới hạn điều kiện tiếp xúc: Giả thiết liên kết giữa các lớp composite là dính bám hoàn hảo, chưa xét đến hiện tượng phân tách lớp (delamination) hoặc suy giảm liên kết do môi trường ẩm nhiệt.
  4. Giới hạn thực nghiệm: Thí nghiệm mới thực hiện với nước sạch trong điều kiện phòng thí nghiệm tiêu chuẩn, chưa khảo sát thực nghiệm với chất lỏng nhớt hoặc môi trường nhiệt độ cao.

Chương trình nghiên cứu 5–10 năm tới:

  • Phát triển phần tử liên tục phi tuyến hình học xét đến hiệu ứng sóng nước lắc mạnh (liquid sloshing) và tải trọng va đập thủy động lực học.
  • Mở rộng thuật toán PTLT cho kết cấu vỏ composite chức năng biến đổi cơ tính (FGM) và vật liệu composite áp điện thông minh (Smart Piezoelectric Composites) có gân gia cường.
  • Tích hợp mô hình lan truyền vết nứt và tách lớp phân lớp dưới tác động mỏi dao động thủy đàn hồi.
  • Ứng dụng kỹ thuật tối ưu hóa đa mục tiêu (Genetic Algorithm, PSO) để thiết kế góc nghiêng lớp sợi tối ưu triệt tiêu rung động cho bồn chứa nhiên liệu hàng không vũ trụ.

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

Luận án tạo ra những bước tiến đột phá có sức lan tỏa sâu rộng trong cộng đồng khoa học và ứng dụng công nghiệp:

  • Tác động học thuật: Cung cấp công bố khoa học chất lượng cao, tạo nền tảng trích dẫn vững chắc cho các nghiên cứu tiếp theo về cơ học tính toán FSI của kết cấu composite thành mỏng tại khu vực và quốc tế.
  • Chuyển đổi công nghiệp: Giúp các doanh nghiệp chế tạo composite trong nước nâng cao năng lực tự chủ thiết kế, giảm chi phí thử nghiệm thực nghiệm tốn kém nhờ mô hình mô phỏng PTLT chính xác tuyệt đối, nâng cao tuổi thọ và độ an toàn của bồn bể công nghiệp lên hơn $30%$.
  • Ý nghĩa kinh tế – xã hội: Ngăn ngừa nguy cơ sự cố nứt vỡ, rò rỉ hóa chất độc hại từ các bồn chứa composite do cộng hưởng rung động, bảo vệ an toàn môi trường và người lao động.

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

  • Nghiên cứu sinh & Giảng viên ngành Cơ học/Kỹ thuật: Tiếp cận phương pháp luận PTLT tiên tiến, mã nguồn thuật toán Matlab và quy trình thực nghiệm đo dao động vỏ tròn xoay chuẩn mực.
  • Kỹ sư R&D Hàng không, Tàu thủy & Dầu khí: Sở hữu công cụ số chính xác cao để tính toán tối ưu hóa kết cấu vỏ ghép nón–trụ composite chứa nhiên liệu, dầu mỏ và hóa chất.
  • Doanh nghiệp sản xuất Composite: Áp dụng quy chuẩn thiết kế lớp sợi ($[0^\circ/90^\circ]_s$) và kiểm soát chiều dày bậc hợp lý để loại bỏ nguy cơ cộng hưởng cơ học.
  • Cơ quan quản lý tiêu chuẩn chất lượng: Có căn cứ khoa học vững chắc để xây dựng quy chuẩn kỹ thuật quốc gia về thử nghiệm động lực học kết cấu composite chứa chất lỏng.

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

1. Đóng góp lý thuyết độc đáo nhất của nghiên cứu là gì?

Đóng góp độc đáo nhất là việc tích hợp thành công Lý thuyết biến dạng cắt bậc nhất Reissner–Mindlin (FSDT) với Lý thuyết thế lưu chất Laplace để thiết lập ma trận độ cứng động lực học chính xác $[K(\omega)]$ cho vỏ composite tròn xoay ghép nối đa đoạn (trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón). Nghiên cứu đã chuyển hóa bài toán tương tác lưu chất phức tạp thành hệ không gian trạng thái 10 ẩn số với nghiệm giải tích hàm mũ ma trận, giải quyết triệt để tương tác màng–uốn–xoắn trong môi trường lưu chất.

2. Đổi mới 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 phương pháp PTHH truyền thống (phụ thuộc mật độ lưới và tốn tài nguyên tính toán của Cho, 2002; Patel, 2000) và phương pháp ma trận truyền cổ điển dùng thuật toán Wittrick–Williams (dễ gặp lỗi kỳ dị của Amabili, 1996; Zhang & Xiang, 2007), phương pháp PTLT kết hợp Kỹ thuật đường cong đáp ứng của luận án chỉ cần giải phương trình đại số tuyến tính ${F} = [K(\omega)]{U}$ với lực kích thích đơn vị, loại bỏ hoàn toàn việc dò nghiệm định thức, cho kết quả tần số riêng và mode dao động nhanh chóng, chính xác tuyệt đối.

3. Phát hiện bất ngờ và phản trực giác nhất thu được từ dữ liệu là gì?

Phát hiện bất ngờ nhất là tần số dao động tự do nhỏ nhất của vỏ tròn xoay ghép nối chứa chất lỏng không rơi vào mode sóng vòng $m = 1$ (dạng uốn uốn cong trục như dầm) mà rơi vào các mode sóng vòng bậc cao ($m = 3, 4, 5, 6$). Hơn nữa, khi tăng chiều cao mực chất lỏng $H$, ngoài việc tần số giảm sâu do khối lượng ảo $m^*$, thứ tự xuất hiện của các mode dao động bị xáo trộn mạnh, chứng minh rằng không thể áp dụng các công thức tính toán đơn giản hóa của dầm tương đương cho vỏ mỏng chứa lưu chất.

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

Luận án cung cấp chi tiết toàn diện quy trình tái lập: từ việc xây dựng tường minh các phần tử của ma trận trạng thái $[A_m]$ kích thước $10 \times 10$, thuật toán giải ma trận truyền trong Matlab, sơ đồ ghép nối ma trận tổng thể $15 \times 15$ và $30 \times 30$, đến quy trình chế tạo mẫu composite sợi thủy tinh/polyester, thông số đồ gá ngàm và sơ đồ bố trí thiết bị đo thực nghiệm FFT.

5. Kế hoạch nghiên cứu 10 năm tới định hình như thế nào?

Kế hoạch 10 năm tập trung phát triển bộ phần mềm chuyên dụng tính toán động lực học thủy đàn hồi phi tuyến cho vỏ composite FGM và vật liệu thông minh đa trường vật lý (nhiệt – điện – từ – cơ), mở rộng cho các bài toán vỏ chịu tải trọng nổ dưới nước, va đập thủy lực và tối ưu hóa cấu trúc composite chống rung động trong công nghiệp không gian.

Kết luận

  1. Xây dựng thành công hệ phương pháp luận và thuật toán Phần tử liên tục (PTLT) có độ chính xác cao dựa trên lý thuyết Reissner–Mindlin để phân tích dao động tự do của các kết cấu vỏ tròn xoay composite phân lớp có mặt cắt ngang thay đổi (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 tường minh ma trận trạng thái 10 thành phần và ma trận độ cứng động lực ghép nối đa đoạn, tích hợp chính xác áp suất thủy động qua hàm thế vận tốc Bessel và phương trình Bernoulli.
  3. Đề xuất ứng dụng thành công Phương pháp đường cong đáp ứng trong giải ma trận độ cứng động lực, khắc phục triệt để các nhược điểm tính toán của phương pháp rời rạc hóa PTHH và thuật toán Wittrick–Williams.
  4. Đánh giá định lượng toàn diện ảnh hưởng của mực chất lỏng, khối lượng riêng vật liệu/chất lỏng, góc mở nón $\alpha$, tỷ lệ hình học $L/R$ và cấu hình phân lớp đối xứng/phản đối xứng đến tần số riêng và sự dịch chuyển số sóng mode vòng $m$.
  5. Chế tạo thành công các bộ mẫu thực tế và thực hiện đo đạc thực nghiệm modal trên vỏ composite nón–trụ và nón–trụ–nón, cung cấp bộ dữ liệu thực nghiệm quý giá kiểm chứng xuất sắc mô hình lý thuyết với sai số chỉ từ $2.1%$ đến $7.6%$.
  6. Mở ra hướng nghiên cứu mới về tối ưu hóa kết cấu composite thành mỏng làm việc trong môi trường tương tác lưu chất – kết cấu phức tạp, đóng góp trực tiếp vào sự phát triển của nền công nghiệp cơ khí, đóng tàu và hàng không vũ trụ hiện đại.