Tổng quan về luận án

Nghiên cứu động lực học kết cấu thành mỏng làm bằng vật liệu composite lớp đóng vai trò then chốt trong kỹ thuật hiện đại. Sự kết hợp giữa tỷ gia độ bền/trọng lượng cao, khả năng chống ăn mòn hóa chất và tính linh hoạt trong gia công định hình khiến các kết cấu vỏ tròn xoay liên hợp—bao gồm vỏ trụ bậc, vỏ nón–trụ, vỏ nón–trụ–nón và vỏ nón–nón–nón—trở thành giải pháp kết cấu chiến lược trong công nghiệp hàng không vũ trụ, bồn bể hóa chất (điển hình tại Nhà máy Hóa chất Việt Trì, Nhà máy sản xuất Amon Nitrat Thái Bình), công trình biển và kỹ thuật năng lượng hạt nhân. Tuy nhiên, khi các kết cấu này chứa chất lỏng, hiện tượng tương tác thủy–đàn hồi (hydroelastic interaction) làm biến đổi sâu sắc các đặc trưng động lực học, gây suy giảm mạnh tần số dao động tự do và thay đổi trật tự các dạng mode dao động.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ: Các công trình quốc tế trước đây chủ yếu giải quyết bài toán dao động của vỏ kim loại đẳng hướng tiếp xúc chất lỏng (Amabili, 1996; Damatty et al., 2001, 2005) hoặc dừng lại ở các cấu trúc vỏ composite đơn lẻ có tiết diện không đổi (Lakis et al., 1991, 1997; Toorani & Lakis, 2001). Đối với các hệ vỏ composite tròn xoay liên hợp có mặt cắt ngang thay đổi đột ngột (trụ bậc) hoặc biến thiên liên tục dọc trục (hệ nón–trụ, nón–trụ–nón, nón–nón–nón) chứa các mực chất lỏng khác nhau, việc mô hình hóa gặp rào cản toán học lớn do tính dị hướng cơ học màng–uốn–xoắn–cắt kết hợp với điều kiện ghép nối động lực học phức tạp. Phương pháp phần tử hữu hạn (FEM) truyền thống bộc lộ nhược điểm tiêu tốn tài nguyên và dễ phát sinh sai số rời rạc hóa ở miền tần số cao, trong khi phương pháp giải tích thuần túy không thể áp dụng cho các hình học không đồng nhất.

Nhằm giải quyết triệt để bài toán trên, luận án tiến sĩ cơ học của tác giả Vũ Quốc Hiến (2017) với đề tài "Nghiên cứu dao động của vỏ composite tròn xoay chứa chất lỏng" (chuyên ngành Cơ kỹ thuật, mã số 62510101, Trường Đại học Bách khoa Hà Nội; người hướng dẫn: GS. Trần Ích Thịnh và PGS. Nguyễn Mạnh Cường) đã tập trung vào các câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một mô hình toán học giải tích–số chính xác cao dựa trên phương pháp phần tử liên tục (PTLT) nhằm xác định ma trận độ cứng động lực của các phân đoạn vỏ composite tròn xoay tiếp xúc chất lỏng lý tưởng?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ suy giảm tần số dao động riêng và sự dịch chuyển mode dao động vòng ($m$), mode dao động dọc ($n$) của hệ vỏ liên hợp biến thiên ra sao dưới tác động đồng thời của mức chất lỏng ($H/L$), góc côn ($\alpha$), tỷ lệ kích thước ($L/R$), điều kiện biên và cấu hình sắp xếp lớp composite?
  • Giả thuyết 1 (H1): Phương pháp ma trận độ cứng động lực (DSM/PTLT) kết hợp với kỹ thuật quét đường cong đáp ứng tần số cho phép triệt tiêu hoàn toàn sai số rời rạc hóa không gian của FEM, mang lại nghiệm dao động có độ hội tụ tuyệt đối với số bậc tự do tối thiểu.
  • Giả thuyết 2 (H2): Sự tương tác thủy động giữa áp suất chất lỏng (xác định qua hàm thế vận tốc Bessel) và thành vỏ composite gây ra hiệu ứng khối lượng bổ sung (added mass effect) phi tuyến tính, làm thay đổi dạng mode cơ bản từ $m=1$ sang các mode sóng vòng bậc cao hơn ($m \ge 2$) tùy thuộc vào mức độ đổ đầy chất lỏng.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa Lý thuyết biến dạng cắt bậc nhất của Reissner–Mindlin (FSDT) với Lý thuyết thủy động lực học chất lỏng lý tưởng (phương trình Laplace và phương trình Bernoulli không dừng). Phạm vi nghiên cứu bao quát các hệ vỏ trụ bậc (lên đến 4 bậc), vỏ ghép nối nón–trụ, nón–trụ–nón và nón–nón–nón bằng composite cốt sợi thủy tinh/nền nhựa polyester không no có cấu hình lớp $[0^\circ/90^\circ/0^\circ/90^\circ]$ hoặc đối xứng $[0^\circ/90^\circ]_s$, làm việc trong miền đàn hồi tuyến tính và chịu các điều kiện biên ngàm–tự do (C–F), ngàm–ngàm (C–C), ngàm–tựa (C–S).

Literature Review và Positioning

Lịch sử nghiên cứu dao động vỏ tròn xoay tương tác chất lỏng trải qua ba trường phái phương pháp luận chính:

Trường phái giải tích cổ điển đặt nền móng với các công trình của Amabili (1996, 2002), Lakis & Laveau (1991), Toorani & Lakis (2001). Các tác giả này ứng dụng phương pháp Ritz, Galerkin kết hợp với hàm sóng giải tích để giải phương trình vi phân chuyển động. Ưu điểm nổi bật là tính chính xác tuyệt đối về mặt toán học. Tuy nhiên, ranh giới hạn chế nghiêm trọng của trường phái này là chỉ giải quyết được các bài toán vỏ đơn lẻ có hình học đồng nhất và điều kiện biên đơn giản (như tựa đơn giản S–S hoặc ngàm hai đầu C–C). Khi hình học chuyển tiếp phức tạp (nón nối trụ, trụ đổi chiều dày), các hệ phương trình vi phân liên kết không còn khả năng tìm nghiệm giải tích tường minh.

Trường phái số hóa bằng phần tử hữu hạn (FEM) và phần tử biên (BEM) phát triển mạnh mẽ qua các nghiên cứu của Cho et al. (2002), Kim et al. (2002), Shekari et al. (2009, 2010), Patel et al. (2000, 2014) và Damatty et al. (2001, 2005). Phương pháp này chia lưới kết cấu thành các phần tử hữu hạn và sử dụng các hàm dạng chuyển vị tĩnh. Mặc dù rất linh hoạt trong việc mô hình hóa hình học bất kỳ, FEM gặp phải mâu thuẫn lớn: Khi phân tích miền tần số cao hoặc cấu trúc mỏng chịu tải thủy động, kích thước lưới phải cực mịn để tránh hiện tượng khóa cắt (shear locking) và sai số cắt cụt, dẫn đến kích thước ma trận bùng nổ, thời gian tính toán kéo dài và nguy cơ tích tụ sai số làm tròn số học.

Trường phái ma trận độ cứng động lực (DSM) hay Phương pháp phần tử liên tục (PTLT) khởi xướng bởi Clough (1975), Hallauer (1985), Kulla (1985, 2003), Gorman (1982, 1996) và Banerjee (1989, 1996, 2005) cho hệ dầm và tấm. Nguyen Manh Cuong & Casimir (2007), Tran Ich Thinh & Nguyen Manh Cuong (2013, 2016), và Ta Thi Hien (2014) đã tiên phong mở rộng PTLT sang phân tích vỏ trụ và vỏ nón kim loại/composite đơn lẻ. PTLT xem mỗi phân đoạn kết cấu là một "phần tử liên tục" có vô số bậc tự do, sử dụng chính nghiệm động lực học chính xác của phương trình vi phân chuyển động làm hàm dạng trong miền tần số.

                           KHÔNG GIAN LÝ THUYẾT & ĐỊNH VỊ NGHIÊN CỨU

Luận án của Vũ Quốc Hiến định vị đột phá tại giao điểm của cơ học vật liệu composite dị hướng và thủy động lực học tương tác: Lần đầu tiên thiết lập hoàn chỉnh thuật toán PTLT ghép nối cho chuỗi vỏ tròn xoay liên hợp đa phân đoạn (trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón) chứa chất lỏng, vượt lên trên các nghiên cứu quốc tế cùng thời (như Qu et al., 2013 chỉ phân tích vỏ trụ bậc kim loại bằng biến phân; Patel et al., 2014 chỉ tính nón–trụ composite khô bằng FEM).

Đó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 vỏ Reissner–Mindlin (FSDT) cho bài toán tương tác thủy–đàn hồi của kết cấu vỏ composite lớp dị hướng. Bằng cách thiết lập quan hệ ứng suất–biến dạng trên trạng thái biến dạng trượt bậc nhất, trường chuyển vị tại điểm có tọa độ $(x, \theta, z)$ được mô tả chính xác:

$$u(x, \theta, z, t) = u_0(x, \theta, t) + z \varphi_x(x, \theta, t)$$ $$v(x, \theta, z, t) = v_0(x, \theta, t) + z \varphi_\theta(x, \theta, t)$$ $$w(x, \theta, z, t) = w_0(x, \theta, t)$$

Trong đó $u_0, v_0, w_0$ là các chuyển vị tại mặt trung bình; $\varphi_x, \varphi_\theta$ là các góc xoay của pháp tuyến quanh trục $\theta$ và trục $x$. Quan hệ nội lực – chuyển vị tích hợp ma trận độ cứng vật liệu $N$ lớp:

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

Các hệ số độ cứng màng $[A]$, tương tác màng–uốn–xoắn $[B]$, độ cứng uốn $[D]$ và độ cứng cắt $[F]$ (với hệ số hiệu chỉnh cắt $f = 5/6$) được tích phân trực tiếp qua chiều dày vỏ từ ma trận độ cứng thu gọn $[Q_{ij}]$ của từng lớp trực hướng:

$$A_{ij} = \sum_{k=1}^N Q_{ij}^{(k)} (z_k - z_{k-1}), \quad B_{ij} = \frac{1}{2} \sum_{k=1}^N Q_{ij}^{(k)} (z_k^2 - z_{k-1}^2), \quad D_{ij} = \frac{1}{3} \sum_{k=1}^N Q_{ij}^{(k)} (z_k^3 - z_{k-1}^3)$$

Đóng góp lý thuyết then chốt là việc giải tích hóa áp suất thủy động $P(r, \theta, x, t)$ tác động lên thành vỏ từ nghiệm của phương trình sóng Laplace cho chất lỏng lý tưởng (không nén được, không nhớt, chuyển động không xoáy):

$$\nabla^2 \Phi = \frac{\partial^2 \Phi}{\partial r^2} + \frac{1}{r}\frac{\partial \Phi}{\partial r} + \frac{1}{r^2}\frac{\partial^2 \Phi}{\partial \theta^2} + \frac{\partial^2 \Phi}{\partial x^2} = 0$$

Áp dụng phương trình Bernoulli không dừng và điều kiện biên tiếp xúc động lực học $\left. \frac{\partial \Phi}{\partial r} \right|_{r=R} = \frac{\partial w_0}{\partial t}$, áp suất chất lỏng tác động lên bề mặt trong của vỏ được biểu diễn dưới dạng tường minh:

$$P(R, \theta, x, t) = m^* \omega^2 w_0 \cos(m\theta) e^{i\omega t}$$

Với hệ số khối lượng gia tăng thủy động $m^*$ được dẫn xuất thông qua hàm Bessel hiệu chỉnh loại 1 bậc $m$ ($I_m$):

$$m^* = \frac{\rho_f}{k_n} \frac{I_m(k_n R)}{I'm(k_n R)} = \frac{\rho_f}{k_n} \frac{I_m(k_n R)}{I{m-1}(k_n R) - \frac{m}{k_n R} I_m(k_n R)}$$

Số sóng dọc $k_n$ phụ thuộc nghiêm ngặt vào điều kiện biên của miền chất lỏng ($k_n = n\pi/L$ đối với biên ngàm–ngàm; $k_n = (2n-1)\pi/(2L)$ đối với biên ngàm–tự do).

                             KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN

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

Khung phân tích của luận án tích hợp ba trụ cột:

  1. Biểu diễn chuỗi Fourier đồng bộ: Triển khai các thành phần chuyển vị và nội lực theo hàm lượng giác hội tụ đều:

$${u_0, w_0, \varphi_x, N_x, M_x, Q_x} = \sum_{m=1}^\infty {u_m(x), w_m(x), \varphi_{xm}(x), N_{xm}(x), M_{xm}(x), Q_{xm}(x)} \cos(m\theta) e^{i\omega t}$$ $${v_0, \varphi_\theta, N_{x\theta}, M_{x\theta}, Q_\theta} = \sum_{m=1}^\infty {v_m(x), \varphi_{\theta m}(x), N_{x\theta m}(x), M_{x\theta m}(x), Q_{\theta m}(x)} \sin(m\theta) e^{i\omega t}$$

  1. Hệ phương trình vi phân vector trạng thái 10 thành phần: Chuyển đổi toàn bộ hệ 5 phương trình chuyển động vi phân đạo hàm riêng cấp hai thành hệ 10 phương trình vi phân thường cấp một:

$$\frac{d{y}_m(x)}{dx} = [A]_m {y}_m(x)$$

với ${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$ và $[A]_m$ là ma trận hệ số kích thước $10 \times 10$ chứa đầy đủ các đặc trưng hình học, cơ tính vật liệu $[A, B, D, F]$ và số hạng áp suất chất lỏng $m^*\omega^2$. 3. Điều kiện ghép nối liên tục nghiệm ngặt: Tại mặt cắt tiếp giáp $x = x_i$ giữa hai phân đoạn (giữa bậc trụ $i$ và $j$, hoặc giữa vỏ nón và vỏ trụ), 10 điều kiện tương thích lực và chuyển vị được thiết lập tuyệt đối:

$$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}$$ $$N_{xi} = N_{xj}, \quad N_{x\theta i} = N_{x\theta j}, \quad Q_{xi} = Q_{xj}, \quad M_{xi} = M_{xj}, \quad M_{x\theta i} = M_{x\theta j}$$

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism) với lập trường nhận thức luận tất định cơ học (Deterministic Mechanical Epistemology). Thiết kế nghiên cứu kết hợp đa tầng bậc:

  • Cấp độ vi mô: Mô hình hóa lớp vật liệu composite trực hướng đàn hồi với góc xếp sợi tùy biến.
  • Cấp độ cấu trúc: Thiết lập ma trận độ cứng động lực phân đoạn cho từng dạng hình học tròn xoay (vỏ trụ, vỏ nón).
  • Cấp độ ghép nối liên hợp: Hợp nhất các phân đoạn thành ma trận độ cứng tổng thể của hệ kết cấu vỏ phức hợp.
  • Cấp độ thực nghiệm: Kiểm chứng thực tế các mẫu chế tạo bằng kỹ thuật kích động xung lực và phân tích phổ tần số FFT.

Quy trình nghiên cứu rigorous

Quy trình tính toán số và thực nghiệm được chuẩn hóa qua các pha chặt chẽ:

                          QUY TRÌNH THỰC HIỆN NGHIÊN CỨU TOÀN DIỆN

Phương pháp đường cong đáp ứng giải quyết triệt để vấn đề tìm nghiệm của định thức siêu việt $\det|[K(\omega)]| = 0$, vốn gây tắc nghẽn hoặc bỏ sót nghiệm trong thuật toán Wittrick–Williams truyền thống khi ma trận $[K(\omega)]$ chứa các hàm lượng giác và hyperbolic phức tạp.

Data và phân tích

  • Chế tạo mẫu: Mẫu thí nghiệm chế tạo tại Phòng thí nghiệm Vật liệu Composite (Trường ĐHBK Hà Nội) và Viện Nghiên cứu & Chế tạo Tàu thủy (ĐH Thủy sản Nha Trang). Vật liệu sử dụng là composite sợi thủy tinh dệt kết hợp nhựa polyester không no, xếp lớp cân bằng $[0^\circ/90^\circ/0^\circ/90^\circ]$.
  • Thông số cơ lý vật liệu đo đạc:
    • Mô đun đàn hồi dọc sợi: $E_1 = 12.8 \times 10^9 \text{ N/m}^2$
    • Mô đun đàn hồi ngang sợi: $E_2 = 12.8 \times 10^9 \text{ N/m}^2$
    • Mô đun đàn hồi trượt: $G_{12} = 4.28 \times 10^9 \text{ N/m}^2$, $G_{13} = G_{23} = 2.14 \times 10^9 \text{ N/m}^2$
    • Hệ số Poisson: $\nu_{12} = 0.14$
    • Khối lượng riêng vỏ: $\rho_s = 1450 \text{ kg/m}^3$
    • Khối lượng riêng chất lỏng thử nghiệm: $\rho_{\text{nước}} = 1000 \text{ kg/m}^3$, $\rho_{\text{rượu}} = 790 \text{ kg/m}^3$
  • Danh mục mẫu thực nghiệm:
    • Mẫu CT-27: Vỏ nón–trụ, góc nón $\alpha = 27^\circ$, chiều dài nón $L_1 = 200 \text{ mm}$, chiều dài trụ $L_2 = 300 \text{ mm}$, bán kính $R_1 = 50 \text{ mm}, R_2 = 150 \text{ mm}$, chiều dày $h = 2.5 \text{ mm}$.
    • Mẫu CT-14: Vỏ nón–trụ, góc nón $\alpha = 14^\circ$, chiều dài nón $L_1 = 200 \text{ mm}$, chiều dài trụ $L_2 = 300 \text{ mm}$, bán kính $R_1 = 100 \text{ mm}, R_2 = 150 \text{ mm}$, chiều dày $h = 2.5 \text{ mm}$.
    • Mẫu CTC-9: Vỏ nón–trụ–nón, góc nón $\alpha_1 = 9^\circ, \alpha_2 = 9^\circ$, kích thước phân đoạn $L_1 = 150 \text{ mm}, L_2 = 200 \text{ mm}, L_3 = 150 \text{ mm}$, bán kính chuyển tiếp $R_1 = 125 \text{ mm}, R_2 = 150 \text{ mm}, R_3 = 125 \text{ mm}$, chiều dày $h = 2.5 \text{ mm}$.
  • Hệ thống thiết bị đo: Búa gõ xung lực chuẩn tích hợp cảm biến lực PCB Piezotronics, cảm biến gia tốc siêu nhẹ (khối lượng $< 1.5 \text{ g}$ nhằm triệt tiêu ảnh hưởng gia tải phụ), card thu thập dữ liệu đa kênh kết nối phần mềm phân tích phổ FFT chuyên dụng tại Phòng thí nghiệm Kiểm soát Rung và Ồn (Viện Cơ học Việt Nam).

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 $H/L$ tăng từ 0 (vỏ khô) lên 1.0 (chứa đầy), tần số dao động tự do của tất cả các dạng kết cấu vỏ composite giảm mạnh từ 35% đến 68%. Bằng chứng định lượng: Đối với mẫu nón–trụ CT-27 (biên C–F, mode $n=1, m=2$), tần số đo thực nghiệm giảm từ $185.0 \text{ Hz}$ (khi $H/L = 0$) xuống $142.5 \text{ Hz}$ ($H/L = 0.5$) và chỉ còn $107.5 \text{ Hz}$ ($H/L = 1.0$). Kết quả tính toán bằng PTLT tương ứng đạt $189.2 \text{ Hz}$ ($H/L = 0$), $145.8 \text{ Hz}$ ($H/L = 0.5$) và $110.1 \text{ Hz}$ ($H/L = 1.0$). Sai số tương đối giữa PTLT và thực nghiệm luôn duy trì ở mức cực thấp, dao động từ $2.27%$ đến $3.89%$.
  2. Hiện tượng chuyển dịch mode dao động cơ bản (Fundamental Mode Switching): Ở trạng thái vỏ khô ($H/L = 0$), tần số dao động nhỏ nhất thường xuất hiện ở mode sóng vòng thấp ($m = 1$ hoặc $m = 2$). Tuy nhiên, khi mức chất lỏng vượt quá $50%$ thể tích ($H/L \ge 0.5$), mode cơ bản có xu hướng chuyển dịch sang số sóng vòng cao hơn ($m = 3, 4$). Nguyên nhân vật lý là do áp suất thủy động phân bố theo hàm Bessel tác động mạnh nhất vào các biến dạng uốn cục bộ thành mỏng, làm suy giảm độ cứng uốn vòng nhanh hơn độ cứng uốn dọc.
  3. Tác động vượt trội của cấu hình góc nón và tính chất hình học: Đối với hệ vỏ nón–trụ–nón, kết cấu nón phân kỳ tạo ra độ cứng động lực cao hơn đáng kể so với kết cấu nón hội tụ. Tăng góc nón $\alpha$ dẫn đến việc nâng cao tần số riêng ở các mode uốn bậc thấp nhưng làm tăng mật độ các mode dao động cục bộ ở dải tần số trung bình.
  4. Ảnh hưởng của cấu hình sắp xếp lớp composite: Cấu hình lớp đối xứng $[0^\circ/90^\circ]_s$ triệt tiêu ma trận tương tác màng–uốn $[B] = 0$, giúp tần số dao động tự do cao hơn từ $8.5%$ đến $14.2%$ so với cấu hình lớp không đối xứng $[0^\circ/90^\circ]_2$ cùng chiều dày tổng thể do loại bỏ được hiện tượng ghép nối uốn–xoắn gây mất ổn định cục bộ.
  5. Tương quan giữa khối lượng riêng chất lỏng và tần số dao động: Khi thay đổi môi trường chất lỏng từ rượu ($\rho_{\text{rượu}} = 790 \text{ kg/m}^3$) sang nước ($\rho_{\text{nước}} = 1000 \text{ kg/m}^3$), tần số dao động tự do giảm trung bình $7.8% - 11.5%$, tuân theo tỷ lệ nghịch căn bậc hai của khối lượng riêng chất lỏng hiệu dụng: $\omega \propto 1/\sqrt{\rho_s h + m^*(\rho_f)}$.
Mẫu thử nghiệm & Cấu hình Mực chất lỏng ($H/L$) Thực nghiệm (Hz) Lý thuyết PTLT (Hz) Sai số tương đối (%)
CT-27 (Nón–Trụ, $\alpha=27^\circ$) $H/L = 0.00$ (Khô) 185.0 189.2 +2.27%
$H/L = 0.50$ 142.5 145.8 +2.31%
$H/L = 1.00$ (Đầy nước) 107.5 110.1 +2.41%
CT-14 (Nón–Trụ, $\alpha=14^\circ$) $H/L = 0.00$ (Khô) 212.5 218.4 +2.77%
$H/L = 0.50$ 165.0 170.2 +3.15%
$H/L = 1.00$ (Đầy nước) 125.0 129.6 +3.68%
CTC-9 (Nón–Trụ–Nón) $H/L = 0.00$ (Khô) 245.0 252.1 +2.89%
$H/L = 0.50$ 192.5 199.3 +3.53%
$H/L = 1.00$ (Đầy nước) 147.5 153.2 +3.86% 
    Tần số (Hz)
                 H/L = 0            H/L = 0.5          H/L = 1.0 (Mức chất lỏng)

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận chuẩn xác cao để tính toán dao động riêng của các cấu trúc vỏ thành mỏng phức tạp mà không cần chia lưới phần tử hữu hạn, mở ra hướng ứng dụng cho các bài toán tối ưu hóa động lực học kết cấu.
  • Về mặt kỹ thuật công nghiệp: Cung cấp bộ công cụ phần mềm tính toán trong môi trường Matlab có khả năng dự báo chính xác các dải tần số nguy hiểm, giúp kỹ sư thiết kế chủ động tránh được hiện tượng cộng hưởng phá hủy trong các bồn chứa hóa chất composite, tháp phản ứng và thùng nhiên liệu tên lửa/máy bay.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án có các giới hạn biên (boundary conditions) cần được mở rộng:

  1. Giả thiết chất lỏng lý tưởng: Mô hình chưa xét đến độ nhớt của chất lỏng và hiện tượng sóng mặt chuyển động phi tuyến (sloshing effect) khi chất lỏng dao động ở biên độ lớn.
  2. Miền biến dạng tuyến tính: Phân tích dừng lại ở lý thuyết biến dạng bé đàn hồi tuyến tính, chưa xét đến hiệu ứng phi tuyến hình học von Kármán khi kết cấu vỏ siêu mỏng làm việc dưới áp suất động lực cực đại.
  3. Hình học đối xứng trục: Thuật toán PTLT hiện tại chỉ áp dụng cho vỏ tròn xoay hoàn hảo, chưa tích hợp các khuyết tật hình học ban đầu (initial geometric imperfections) phát sinh trong quá trình gia công quấn sợi composite.

Hướng nghiên cứu tương lai:

  • Mở rộng thuật toán PTLT cho chất lỏng có độ nhớt (Navier–Stokes) kết hợp hiệu ứng dập dềnh mặt tự do (free surface sloshing).
  • Nghiên cứu dao động phi tuyến và phân tích ổn định động lực học (dynamic buckling) của vỏ composite liên hợp chịu tải trọng xung nổ hoặc động đất.
  • Ứng dụng vật liệu composite chức năng FGM (Functionally Graded Materials) và vật liệu thông minh (tích hợp các lớp áp điện Piezoelectric/SMA) để điều khiển chủ động rung động của vỏ chứa chất lỏng.

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

  • Tác động học thuật: Luận án tạo ra bộ dữ liệu thực nghiệm chuẩn (benchmark dataset) quý giá về dao động của vỏ nón–trụ và nón–trụ–nón composite chứa chất lỏng, phục vụ đối chứng cho các mô hình số và giải tích trong và ngoài nước.
  • Chuyển giao công nghệ công nghiệp: Các thuật toán và chương trình Matlab phát triển từ luận án có thể ứng dụng trực tiếp trong quy trình thiết kế bồn bể composite tại các đơn vị công nghiệp hóa chất, đóng tàu thủy và quốc phòng tại Việt Nam.
  • Hiệu quả kinh tế – xã hội: Việc làm chủ phương pháp tính toán giúp giảm thiểu hệ số an toàn dư thừa, tối ưu hóa chiều dày kết cấu, tiết kiệm chi phí vật liệu composite nhập khẩu đắt tiền, đồng thời triệt tiêu nguy cơ sự cố nứt vỡ bồn chứa hóa chất nguy hại, bảo vệ môi trường sống.

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

  • Nghiên cứu sinh & Giảng viên ngành Cơ kỹ thuật / Cơ học kết cấu: Tiếp cận một khung phân tích hoàn chỉnh từ giải tích vi phân vector trạng thái, ma trận truyền giải tích đến thuật toán quét đường cong đáp ứng tần số PTLT.
  • Kỹ sư R&D công nghiệp composite, đóng tàu, hàng không: Sở hữu công thức và mã nguồn tính toán tần số dao động tự do nhanh chóng, độ chính xác cao mà không cần đầu tư các gói phần mềm thương mại đắt tiền (ANSYS/Abaqus).
  • Cơ quan quản lý tiêu chuẩn & kiểm định chất lượng: Có cơ sở khoa học định lượng để xây dựng quy chuẩn an toàn động lực học cho bồn bể composite chịu tải chất lỏng tại Việt Nam.

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 thiết lập thành công ma trận độ cứng động lực $10 \times 10$ giải tích cho phân đoạn vỏ nón và vỏ trụ composite lớp dị hướng tiếp xúc chất lỏng lý tưởng. Khác biệt với các nghiên cứu trước chỉ xét vỏ kim loại đẳng hướng hoặc vỏ composite khô, luận án đã tích hợp trực tiếp số hạng áp suất thủy động giải tích $m^\omega^2$ (với $m^$ dẫn xuất từ hàm Bessel) vào ma trận vi phân trạng thái $[A]_m$, cho phép tìm nghiệm động lực học chính xác của hệ vỏ liên hợp mà không cần xấp xỉ trường chuyển vị tĩnh.

2. Phương pháp PTLT trong luận án đổi mới thế nào so với các nghiên cứu quốc tế của Damatty (2001) và Qu (2013)?

  • So với Damatty et al. (2001, 2005): Damatty sử dụng FEM với các phần tử vỏ rời rạc, đòi hỏi chia hàng nghìn phần tử để mô tả vùng tiếp giáp nón–trụ và phải sử dụng phương pháp lặp để giải bài toán tương tác chất lỏng. PTLT trong luận án chỉ cần duy nhất 1 phần tử liên tục cho mỗi phân đoạn đồng nhất, giải trực tiếp trong miền tần số, triệt tiêu sai số rời rạc và rút ngắn thời gian tính toán hơn 90%.
  • So với Qu et al. (2013): Qu sử dụng phương pháp biến phân phân đoạn kết hợp chuỗi Chebyshev cho vỏ trụ bậc kim loại khô. Luận án của Vũ Quốc Hiến đã giải quyết bài toán phức tạp hơn gấp nhiều lần: vỏ composite lớp dị hướng màng–uốn–xoắn, có hình học thay đổi liên hợp (nón nối trụ, trụ bậc) và tương tác trực tiếp với các mực chất lỏng khác nhau.

3. Phát hiện thực nghiệm nào mang tính bất ngờ nhất?

Phát hiện bất ngờ nhất là hiện tượng đảo lộn mode cơ bản khi mức chất lỏng $H/L$ vượt ngưỡng $0.5$. Ở trạng thái khô, kết cấu vỏ nón–trụ CT-27 dao động ở mode cơ bản $m=1$. Tuy nhiên, khi đổ đầy chất lỏng, do áp suất chất lỏng gia tải quán tính không đều theo chu vi, mode cơ bản ứng với tần số nhỏ nhất lại nhảy vọt sang mode sóng vòng $m=3$. Nếu kỹ sư thiết kế chỉ kiểm tra mode $m=1$ theo thói quen tính toán kết cấu khô, công trình sẽ đối mặt với nguy cơ cộng hưởng phá hủy nghiêm trọng ở tần số thấp hơn nhiều so với dự tính.

4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?

Luận án cung cấp đầy đủ và minh bạch toàn bộ giao thức tái lập:

  • Hệ phương trình vi phân trạng thái 10 thành phần với tất cả các hệ số ma trận $A_{ij}$ được viết tường minh.
  • Công thức ma trận truyền giải tích $T(L) = \exp([A]_m L)$ và thuật toán chuyển đổi sang ma trận độ cứng động lực $[K(\omega)]$.
  • Thông số cơ tính chi tiết của vật liệu composite sợi thủy tinh/polyester ($E_1, E_2, G_{12}, \nu_{12}, \rho$).
  • Kích thước hình học chi tiết của các mẫu thí nghiệm CT-27, CT-14, CTC-9 và quy trình đo dao động bằng búa xung lực và xử lý phổ FFT.

5. Lộ trình phát triển 10 năm cho hướng nghiên cứu này là gì?

  • Giai đoạn 1 (1–3 năm): Mở rộng mã nguồn PTLT cho bài toán vỏ composite chịu tải trọng động đất và sóng va đập thủy lực (water hammer), tích hợp thêm hiệu ứng dập dềnh mặt nước (sloshing).
  • Giai đoạn 2 (4–6 năm): Phát triển phần tử liên tục cho vật liệu FGM biến thiên cơ tính theo 3 chiều không gian và vật liệu composite gia cường hạt nano carbon/graphene.
  • Giai đoạn 3 (7–10 năm): Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) kết hợp với PTLT làm mô hình thế mạng (surrogate model) thời gian thực để chẩn đoán hư hỏng, vết nứt ngầm và giám sát sức khỏe kết cấu (Structural Health Monitoring - SHM) cho các bồn chứa nhiên liệu ngầm và vỏ tàu ngầm composite.

Kết luận

Luận án tiến sĩ của tác giả Vũ Quốc Hiến (2017) đại diện cho một bước tiến quan trọng trong lĩnh vực Cơ học kết cấu và Cơ học tính toán tại Việt Nam, mang lại các đóng góp then chốt:

  1. Xây dựng thành công thuật toán và chương trình tính toán số tin cậy cao trên môi trường Matlab bằng Phương pháp phần tử liên tục (PTLT) dựa trên Lý thuyết biến dạng cắt bậc nhất Reissner–Mindlin cho các hệ vỏ tròn xoay liên hợp composite lớp (vỏ trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón) chứa chất lỏng.
  2. Dẫn xuất hoàn chỉnh biểu thức giải tích cho hệ số khối lượng gia tăng thủy động $m^*$ thông qua hàm thế vận tốc Bessel, mô tả chính xác tương tác thủy–đàn hồi giữa chất lỏng lý tưởng và thành vỏ composite dị hướng.
  3. Ứng dụng thành công kỹ thuật quét đường cong đáp ứng tần số để xác định trực tiếp các tần số dao động tự do, khắc phục triệt để nhược điểm mất nghiệm và tốn thời gian của các thuật toán tìm nghiệm định thức truyền thống.
  4. Cung cấp bức tranh toàn diện và định lượng về ảnh hưởng của mức chất lỏng, góc nón, tỷ lệ kích thước $L/R$, điều kiện biên và cấu hình sắp xếp lớp đến đặc trưng dao động tự do của hệ vỏ composite tròn xoay.
  5. Thiết kế, chế tạo và thực hiện thành công chuỗi thí nghiệm đo tần số dao động riêng trên các mẫu vỏ composite liên hợp thực tế (CT-27, CT-14, CTC-9), chứng minh sự tương hợp xuất sắc với lý thuyết PTLT (sai số dưới $4%$).
  6. Mở ra hướng nghiên cứu mới về thiết kế tối ưu hóa động lực học và kiểm soát rung động cho các kết cấu composite mỏng chứa chất lỏng trong các ngành công nghiệp mũi nhọn.