Tổng quan về luận án

Kiểm soát tiếng ồn và dao động trong các hệ thống giao thông vận tải hiện đại (hàng không, đường sắt cao tốc, tàu thủy composite và công nghiệp ô tô) là một trong những thách thức kỹ thuật cốt lõi của thế kỷ 21. Đặc biệt, trường âm bên trong khoang kín chịu sự chi phối mạnh mẽ bởi tương tác động lực học kết cấu - chất lưu (vibro-acoustic coupling) tại dải tần số thấp ($f < 400\text{ Hz}$). Trong bối cảnh vật liệu composite lớp cốt sợi polymer và kết cấu composite sandwich lõi xốp đàn hồi đang dần thay thế các kết cấu kim loại truyền thống nhờ tỷ dung cơ học vượt trội, việc thiết lập mô hình giải tích và dự đoán chính xác tổn thất truyền âm (Sound Transmission Loss - STL) trở thành nhiệm vụ cấp thiết. Luận án tiến sĩ cơ học của nghiên cứu sinh Phạm Ngọc Thành với đề tài "Nghiên cứu sự truyền âm qua tấm composite lớp cốt sợi" (chuyên ngành Cơ học, mã số 9440109, thực hiện tại Trường Đại học Bách Khoa Hà Nội dưới sự hướng dẫn khoa học của GS. TS Trần Ích Thịnh) là công trình tiên phong giải quyết toàn diện bài toán này.

Khoảng trống nghiên cứu (Research Gap) được xác định rõ rệt: trong suốt hơn năm thập kỷ qua, các nghiên cứu cơ học âm thanh kết cấu chủ yếu tập trung vào tấm kim loại đẳng hướng vô hạn hoặc hữu hạn (Maidanik, 1962; Roussos, 1984; Carneal & Fuller, 2004). Khi chuyển sang kết cấu composite cốt sợi trực hướng (orthotropic laminated composite panels) và tấm sandwich lõi xốp polymer (poroelastic polyurethane foam core), các mô hình hiện hữu trên thế giới đa phần giả định kích thước vô hạn để áp dụng phương pháp ma trận truyền cổ điển (Kuo et al., 2008; Bolton et al., 1996), bỏ qua hoàn toàn ảnh hưởng của kích thước hình học hữu hạn và các điều kiện biên thực tế (ngàm - Clamped và tựa bản lề - Simply Supported). Tại Việt Nam, các công trình trước đây chỉ dừng lại ở khảo sát thực nghiệm thuần túy hoặc tiếp cận gần đúng qua phân tích thống kê năng lượng SEA (Đinh Đức Tiến, 2017), chưa xây dựng được hệ phương trình giải tích tường minh kết hợp tương tác sóng ba chiều trong môi trường xốp Biot với lý thuyết tấm composite trực hướng hữu hạn.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H):

  1. RQ1: Ứng xử truyền âm qua tấm composite lớp cốt sợi trực hướng hữu hạn chịu các điều kiện biên ngàm và tựa bản lề thay đổi như thế nào dưới tác động của sóng âm điều hòa có góc tới ($\varphi$) và góc phương vị ($\theta$) tùy ý?
    • H1: Điều kiện biên ngàm (CCCC) làm dịch chuyển các tần số cộng hưởng uốn bậc thấp lên dải tần số cao hơn từ 15% đến 25% so với liên kết tựa bản lề (SSSS), làm thay đổi căn bản quy luật STL trong vùng kiểm soát độ cứng.
  2. RQ2: Cơ chế tương tác âm học trong hệ thống tấm kép composite ngăn cách bởi khoang không khí hữu hạn (panel - air cavity - panel) diễn ra như thế nào khi thay đổi bề dày khoang khí và cấu hình lớp sợi?
    • H2: Sự hình thành các mode dao động âm học trong khoang khí tạo ra hiện tượng cộng hưởng khối lượng - không khí - khối lượng (Mass-Air-Mass resonance), tạo ra vùng suy giảm STL cục bộ nhưng cải thiện mạnh mẽ STL dải tần số cao theo quy luật tăng $12\text{ dB/octave}$.
  3. RQ3: Làm thế nào để mô hình hóa chính xác sự tiêu tán năng lượng âm trong tấm composite sandwich có lớp da composite trực hướng và lớp lõi là vật liệu polymer xốp PU?
    • H3: Sự kết hợp giữa lý thuyết tấm mỏng cổ điển Kirchhoff và lý thuyết môi trường xốp Biot (xét đồng thời 2 sóng dọc nén $k_1, k_2$ và 1 sóng cắt trong hai pha rắn - lỏng tương tác) sẽ triệt tiêu hiện tượng sụt giảm STL tại tần số trùng (coincidence frequency) so với tấm đơn đồng nhất.
  4. RQ4: Độ tin cậy của thuật toán số giải tích tự lập trình trên nền tảng MATLAB được kiểm chứng và định lượng sai số như thế nào khi đối sánh trực tiếp với các phép đo thực nghiệm trong hệ thống phòng phát - phòng thu tiêu chuẩn?
    • H4: Sai số giữa mô hình giải tích chuỗi Fourier kép/Galerkin đề xuất và thực nghiệm đo trên các mẫu sandwich GFRP/PU foam tại Viện Nghiên cứu và Chế tạo Tàu thủy đạt mức dưới $3\text{ dB}$ trên toàn bộ dải tần số khảo sát $20\text{ Hz} - 10.000\text{ Hz}$.

Khung lý thuyết của luận án tích hợp chặt chẽ: Lý thuyết tấm mỏng Kirchhoff mở rộng cho vật liệu composite trực hướng, Lý thuyết truyền âm sóng phẳng trong trường tự do, Lý thuyết môi trường xốp đàn hồi Biot (1956) cho vật liệu xốp bão hòa chất lưu, và Phương pháp hàm Green kết hợp tích phân Rayleigh. Về mặt phạm vi, luận án nghiên cứu các kết cấu tấm phẳng chữ nhật hữu hạn kích thước $a \times b$, chiều dày tấm da $h_1, h_2$, chiều dày lõi $H$, chịu kích thích sóng âm phẳng trong không khí với dải tần số âm tiêu chuẩn từ $20\text{ Hz}$ đến $20.000\text{ Hz}$ (phân tích theo dải 1/3 Octave), mật độ không khí $\rho_0 = 1,21\text{ kg/m}^3$, vận tốc âm $c_0 = 340\text{ m/s}$, áp suất tham chiếu chuẩn $p_{\text{ref}} = 2 \times 10^{-5}\text{ N/m}^2$.

Literature Review và Positioning

Lịch sử nghiên cứu sự truyền âm qua kết cấu ngăn cách khởi nguồn từ các công trình kinh điển của Maidanik (1962) về ứng xử âm học thống kê của tấm kim loại có gân gia cường, chỉ ra mối liên hệ mật thiết giữa phổ gia tốc dao động kết cấu và phổ áp suất bức xạ âm. Craggs và Stead (1976) tiên phong ứng dụng phương pháp phần tử hữu hạn (FEM) để mô phỏng tương tác hệ thống phòng - tấm ngăn cách kim loại, song bị giới hạn bởi số bậc tự do khổng lồ và chỉ hội tụ tại tần số rất thấp. Roussos (1984) đã tạo ra bước đột phá lý thuyết khi thiết lập nghiệm giải tích cho tổn thất truyền âm qua tấm kim loại đẳng hướng và tấm composite trực hướng hữu hạn liên kết tựa bản lề 4 cạnh bằng tích phân hàm Green. Tuy nhiên, mô hình của Roussos bị giới hạn nghiêm ngặt ở điều kiện biên tựa đơn và chưa giải quyết được tấm kép hay tấm sandwich lõi xốp.

Trong phân hệ tấm kép có khoang khí, London (1950) đặt nền móng thực nghiệm đầu tiên về trở kháng âm học của tấm kép. Tiếp sau đó, Carneal và Fuller (2004) phát triển mô hình lý thuyết và thực nghiệm cho tấm nhôm kép hữu hạn tựa bản lề, sử dụng kỹ thuật điều khiển âm thanh kết cấu chủ động (ASAC) và phải nhân đôi độ cứng lý thuyết một cách gượng ép để xấp xỉ điều kiện biên ngàm thực nghiệm. Đối với kết cấu sandwich phức hợp, Bolton và cộng sự (1996) ứng dụng thành công lý thuyết môi trường xốp Biot (1956) để mô tả 3 dạng sóng lan truyền trong lớp xốp polyurethane (PU foam) kẹp giữa hai tấm nhôm đẳng hướng vô hạn. D'Alessandro và cộng sự (2013) cùng Isaac và cộng sự (2020) đã tổng kết toàn diện các trường phái tiếp cận STL, chỉ ra rằng phương pháp ma trận truyền (TMM) tuy nhanh nhưng chỉ áp dụng cho tấm vô hạn; phương pháp năng lượng thống kê (SEA) kém chính xác ở dải tần thấp nơi mật độ mode thưa thớt; trong khi phương pháp số FEM/BEM (Sgard et al., 2000; Panneton & Atalla, 1996) đòi hỏi tài nguyên tính toán cực lớn tại tần số trung và cao.

                  ┌────────────────────────────────────────────────────────┐
                  │          TIẾP CẬN LÝ THUYẾT TRUYỀN ÂM KẾT CẤU          │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
         ┌───────────────────────────────────┴───────────────────────────────────┐
         ▼                                                                       ▼
┌─────────────────────────────────┐                             ┌─────────────────────────────────┐
│     TRƯỜNG PHÁI TẤM VÔ HẠN      │                             │     TRƯỜNG PHÁI TẤM HỮU HẠN     │
│ (Infinite / Transfer Matrix)    │                             │   (Finite Boundary Mechanics)   │
├─────────────────────────────────┤                             ├─────────────────────────────────┤
│ • Kuo et al. (2008): TMM Laminated│                           │ • Roussos (1984): Chỉ SSSS đơn  │
│ • Bolton et al. (1996): Biot PU │                             │ • Carneal & Fuller (2004): ASAC │
│ • Lee & Kondo (2004): Equiv Fluid│                            │ • Sgard et al. (2000): FEM/BEM  │
├─────────────────────────────────┤                             ├─────────────────────────────────┤
│ KHUYẾT TẬT: Bỏ qua mode cục bộ  │                             │ KHUYẾT TẬT: FEM tốn chi phí số; │
│ và điều kiện biên thực tế.      │                             │ Roussos chưa làm sandwich/ngàm. │
└─────────────────────────────────┘                             └─────────────────────────────────┘
                                             │
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │         VỊ TRÍ ĐỊNH VỊ CỦA LUẬN ÁN (PHẠM NGỌC THÀNH)   │
                  ├────────────────────────────────────────────────────────┤
                  │ • Hợp nhất Kirchhoff CLPT + Biot Poroelastic + 3D Wave │
                  │ • Nghiệm tường minh cho cả SSSS và CCCC (Chuỗi kép)   │
                  │ • Khảo sát 3 hệ kết cấu: Tấm đơn, Tấm kép, Sandwich   │
                  │ • Kiểm chứng thực nghiệm phòng phát - thu tiêu chuẩn   │
                  └────────────────────────────────────────────────────────┘

Trong bức tranh toàn cảnh đó, luận án của NCS Phạm Ngọc Thành đã định vị chính xác vị trí học thuật tiên phong: bắc cầu nối hoàn chỉnh giữa cơ học môi trường xốp ba chiều Biot và cơ học kết cấu composite lớp trực hướng hữu hạn có điều kiện biên phức tạp.

Khi so sánh với hai nghiên cứu quốc tế tiêu biểu:

  1. So với Bolton, Shiau và Kang (1996) [J. Sound Vib.]: Công trình của Bolton chỉ khảo sát tấm nhôm đẳng hướng kích thước vô hạn, bỏ qua ảnh hưởng của kích thước hình học hữu hạn và tính dị hướng sợi. Luận án đã mở rộng bài toán sang tấm composite trực hướng hữu hạn với các góc định hướng sợi khác nhau, đồng thời xử lý chính xác điều kiện biên ngàm 4 cạnh thông qua việc khai triển hệ hàm dạng phi lượng giác.
  2. So với Carneal và Fuller (2004) [J. Acoust. Soc. Am.]: Nghiên cứu của Carneal chỉ xét tấm kép nhôm đẳng hướng và gặp bế tắc khi mô hình hóa điều kiện biên ngàm, buộc phải dùng hệ số hiệu chỉnh bán thực nghiệm. Luận án đã thiết lập hệ phương trình vi phân dao động ghép đầy đủ, đưa ra nghiệm giải tích trực tiếp cho cả biên ngàm (CCCC) và biên tựa (SSSS) đối với tấm kép composite trực hướng chứa khoang không khí.

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

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

Luận án đã mở rộng và thách thức các giới hạn của Lý thuyết tấm composite mỏng cổ điển (Classical Laminated Plate Theory - CLPT) và Lý thuyết môi trường xốp Biot bằng các đóng góp mang tính đột phá:

  • Mở rộng phương trình dao động vi phân ghép âm - kết cấu cho tấm trực hướng hữu hạn: Thiết lập hệ phương trình vi phân chuyển vị ngang $w(x,y,t)$ của tấm composite lớp trực hướng chịu kích thích áp suất âm điều hòa $p_i(x,y,z,t)$ với ma trận độ cứng uốn $D_{ij}$ ($i,j = 1, 2, 6$): $$D_{11}\frac{\partial^4 w}{\partial x^4} + 2(D_{12} + 2D_{66})\frac{\partial^4 w}{\partial x^2 \partial y^2} + D_{22}\frac{\partial^4 w}{\partial y^4} + m^* \frac{\partial^2 w}{\partial t^2} = p(x,y,0,t)$$ trong đó $m^* = \sum_{k=1}^N \rho_k h_k$ là mật độ khối lượng bề mặt, $p(x,y,0,t) = 2p_i(x,y,0,t) - p_r(x,y,0,t) - p_t(x,y,0,t)$ là áp suất âm tổng tác dụng lên tấm bao gồm áp suất sóng tới, sóng phản xạ và sóng truyền qua bức xạ.

  • Hợp nhất lý thuyết môi trường xốp Biot với điều kiện biên kết cấu hữu hạn: Trong lớp lõi xốp polymer (PU foam), sóng âm lan truyền dưới dạng 3 sóng riêng biệt: hai sóng giãn dọc truyền trong pha rắn và pha lỏng ($k_1, k_2$) và một sóng cắt ($k_3$). Luận án đã thiết lập hệ phương trình ứng suất - biến dạng và chuyển vị của khung rắn $\mathbf{u}^s$ và chất lưu $\mathbf{u}^f$: $$\sigma_{ij}^s = 2N e_{ij}^s + (A e^s + Q e^f)\delta_{ij}$$ $$\sigma_{ij}^f = (Q e^s + R e^f)\delta_{ij} = -\sigma p^f \delta_{ij}$$ với $A, N$ là các hệ số Lamé của khung rắn, $Q$ là hệ số ghép thể tích giữa hai pha, $R$ là mô đun bulk của chất lưu trong lỗ xốp, và $\sigma$ là độ rỗng (porosity). Luận án đã giải quyết thành công sự tương thích chuyển vị và liên tục ứng suất tại mặt phân cách giữa lớp da composite trực hướng và lớp lõi xốp Biot.

  • Thiết lập công thức tổng quát tính hệ số truyền công suất âm $\tau(\varphi, \theta)$ và tổn thất truyền âm $STL$: $$STL = -10 \log_{10} (\tau_d) = -10 \log_{10} \left( \frac{\Pi_t}{\Pi_i} \right)$$ trong đó công suất âm tới $\Pi_i = \frac{|p_0|^2 a b \cos\varphi}{2 \rho_0 c_0}$ và công suất âm truyền qua $\Pi_t = \frac{1}{2} \text{Re} \left[ \int_0^a \int_0^b p_t(x,y,0) \cdot v^*(x,y) , dx dy \right]$.

┌────────────────────────────────────────────────────────────────────────────────────────────────┐
│                             KHUNG PHÂN TÍCH LÝ THUYẾT ĐỒNG BỘ CỦA LUẬN ÁN                      │
└───────────────────────────────────────────────┬────────────────────────────────────────────────┘
                                                │
       ┌────────────────────────────────────────┼────────────────────────────────────────┐
       ▼                                        ▼                                        ▼
┌──────────────────────────────┐ ┌──────────────────────────────┐ ┌──────────────────────────────┐
│       PHÂN HỆ 1: TẤM ĐƠN     │ │       PHÂN HỆ 2: TẤM KÉP     │ │    PHÂN HỆ 3: TẤM SANDWICH   │
│ Composite Lớp Trực Hướng     │ │ Composite Khoang Khí Hữu Hạn │ │ Composite Lõi PU Xốp Biot    │
├──────────────────────────────┤ ├──────────────────────────────┤ ├──────────────────────────────┤
│ • Lý thuyết: Kirchhoff CLPT  │ │ • Ghép 2 tấm + Khí động học  │ │ • Tấm da: CLPT trực hướng    │
│ • Dao động uốn 2D dị hướng   │ │ • Mode áp suất khoang kín    │ │ • Lõi xốp: 3D Biot (2P + 1S) │
│ • Tích phân bức xạ Rayleigh  │ │ • Trở kháng âm khoang khí    │ │ • Ma trận ghép truyền biên   │
└──────────────────────────────┘ └──────────────────────────────┘ └──────────────────────────────┘
                                                │
                                                ▼
┌────────────────────────────────────────────────────────────────────────────────────────────────┐
│             PHƯƠNG PHÁP TOÁN HỌC ĐỒNG NHẤT: CHUỖI FOURIER KÉP & PHƯƠNG PHÁP GALERKIN           │
│  - Điều kiện biên tựa (SSSS): $\phi_{mn}(x,y) = \sin(m\pi x/a)\sin(n\pi y/b)$                 │
│  - Điều kiện biên ngàm (CCCC): $\phi_{mn}(x,y) = X_m(x)Y_n(y)$ (Hàm chùm tia rung)            │
│  - Thiết lập ma trận tổng quát $[Z_{mn,kl}] \{W_{kl}\} = \{F_{mn}\}$ $\rightarrow$ Giải trên MATLAB│
└────────────────────────────────────────────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp đồng thời ba trụ cột lý thuyết: (1) Cơ học vật liệu composite nhiều lớp trực hướng anisotropic; (2) Âm học sóng và lý thuyết bức xạ Helmholtz - Rayleigh; (3) Nhiệt động - âm học chất lưu nhớt trong môi trường xốp Biot.

Điểm độc đáo cốt lõi nằm ở việc xây dựng toán tử chuyển đổi ma trận trở kháng bức xạ kết cấu $[T_{ij}]$ cho phép thay thế linh hoạt điều kiện biên từ tựa đơn sang ngàm bốn cạnh mà không phá vỡ tính liên tục của trường áp suất âm phản xạ và truyền qua. Biên điều kiện (Boundary conditions) được kiểm soát nghiêm ngặt qua hệ số hàm dạng $Q_{i,mn}$ và hằng số điều kiện biên $f_{mn}$, đảm bảo mô hình hội tụ chính xác ngay cả khi tỷ số mô đun đàn hồi dọc/ngang $E_{11}/E_{22}$ biến thiên từ $1$ (đẳng hướng) đến $40$ (composite Graphite/Epoxy cao cấp).

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ủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp mô hình kiểm chứng thực nghiệm tất định. Quy trình nghiên cứu kết hợp đa tầng:

  1. Tầng mô hình hóa giải tích (Analytical Modeling): Thiết lập hệ phương trình vi phân đạo hàm riêng cấp 4 và cấp 6 mô tả dao động âm kết hợp.
  2. Tầng mô phỏng số (Computational Mechanics): Xây dựng các script giải thuật chuyên biệt trên môi trường MATLAB R2020a, xử lý tích phân 4 lớp phức tạp của hàm Green và giải hệ phương trình đại số tuyến tính phức kích thước lớn.
  3. Tầng thực nghiệm (Empirical Validation): Chế tạo mẫu thử thực tế, thiết kế và hiệu chuẩn hệ thống buồng đo vang âm tiêu chuẩn (reverberation room suite) tại Viện Nghiên cứu và Chế tạo Tàu thủy - Đại học Nha Trang.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các tiêu chuẩn đo lường âm học quốc tế khắt khe (tương thích ISO 10140, ISO 3741, ASTM C423 và ASTM E90):

  • Hiệu chuẩn buồng đo (Chamber Calibration): Thiết lập hệ thống gồm phòng phát (Source room) có thể tích lớn và phòng thu (Receiving room) có thể tích $V_2 = 16\text{ m}^3$, bề mặt được ốp gạch men phản xạ cao nhằm tối đa hóa thời gian vang $T$, bên ngoài bọc lớp xốp hấp thụ cách âm hoàn toàn với môi trường xung quanh. Khung gá mẫu thử có diện tích hữu hạn $A_p$ được cố định ngàm chặt bằng hệ gioăng cao su và khung thép nặng để triệt tiêu hiện tượng truyền âm sườn (flanking transmission).
  • Quy trình đo thực nghiệm: Tổn thất truyền âm $STL$ được xác định thông qua độ chênh lệch mức áp suất âm trung bình giữa phòng phát ($L_1$) và phòng thu ($L_2$) cùng hệ số hiệu chỉnh thời gian vang: $$STL = L_1 - L_2 + 10 \log_{10} \left( \frac{A_p \cdot T}{0,16 \cdot V_2} \right)$$
  • Triệt tiêu sai số và kiểm định độ tin cậy: Tiến hành hiệu chuẩn hệ thống đo bằng tấm thép tiêu chuẩn ($E = 210\text{ GPa}, \rho = 7800\text{ kg/m}^3$) và tấm nhôm tiêu chuẩn ($E = 70\text{ GPa}, \rho = 2700\text{ kg/m}^3$). Kết quả đo thực nghiệm trên tấm thép chuẩn hoàn toàn trùng khớp với đường cong lý thuyết quy luật khối lượng (Mass Law) với độ lệch chuẩn cực tiểu ($R^2 > 0,98$), khẳng định tính hợp thức của hệ thống thực nghiệm trước khi đo các mẫu composite.
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                           QUY TRÌNH THỰC NGHIỆM VÀ TÍNH TOÁN STL ĐỒNG BỘ                        │
└────────────────────────────────────────────────┬────────────────────────────────────────────────┘
                                                 │
      ┌──────────────────────────────────────────┴──────────────────────────────────────────┐
      ▼                                                                                     ▼
┌─────────────────────────────────────────┐                           ┌─────────────────────────────────────────┐
│     HỆ THỐNG PHÒNG ĐO THỰC NGHIỆM       │                           │     GIẢI THUẬT TÍNH TOÁN SỐ MATLAB      │
├─────────────────────────────────────────┤                           ├─────────────────────────────────────────┤
│ • Phòng phát: Nguồn âm đa hướng         │                           │ • Nhập cơ tính $E_1, E_2, G_{12}, \nu$  │
│ • Phòng thu: $V_2 = 16\text{ m}^3$      │                           │ • Nhập thông số Biot: $\sigma, \psi, N$ │
│ • Mẫu thử: Khung ngàm chuẩn $A_p$       │                           │ • Thiết lập ma trận trở kháng âm $[Z]$  │
│ • Đo $L_1, L_2$ & Thời gian vang $T$    │                           │ • Tích phân Gauss-Legendre 4 lớp        │
└─────────────────────────────────────────┘                           └─────────────────────────────────────────┘
                                                 │                                          │
                                                 └────────────────────┬─────────────────────┘
                                                                      │
                                                                      ▼
                                      ┌────────────────────────────────────────────────────────┐
                                      │   ĐỐI SOÁNH, ĐỊNH LƯỢNG SAI SỐ & TỐI ƯU HÓA KẾT CẤU    │
                                      │   (Độ lệch trung bình $< 2\text{ dB}$, $f \in [20-10k\text{ Hz}]$)│
                                      └────────────────────────────────────────────────────────┘

Data và phân tích

Luận án đã thu thập và xử lý ma trận số liệu cơ - âm phức tạp của các hệ vật liệu chế tạo thực tế:

  • Lớp da Composite (GFRP): Cốt sợi thủy tinh E-glass dệt (Woven Roving) kết hợp nhựa polyester không no, khối lượng riêng $\rho = 1750 - 1950\text{ kg/m}^3$, mô đun Young dọc sợi $E_{11} = 21,5 - 28,0\text{ GPa}$, $E_{22} = 7,2 - 8,5\text{ GPa}$, mô đun cắt $G_{12} = 3,5 - 4,2\text{ GPa}$, hệ số Poisson $\nu_{12} = 0,28$.
  • Lớp lõi xốp Polyurethane (PU Foam): Khảo sát 2 nhóm khối lượng riêng danh định $\rho_c = 100\text{ kg/m}^3$ ($E_s = 10\text{ MPa}, G_c = 4\text{ MPa}$) và $\rho_c = 190\text{ kg/m}^3$ ($E_s = 30\text{ MPa}, G_c = 12\text{ MPa}$), trở kháng dòng không khí $\psi = 2,5 \times 10^4\text{ Ns/m}^4$, độ rỗng $\sigma = 0,92$.
  • Ma trận mẫu thử thực nghiệm (Mẫu A, C, E, F, G, H, I, K):
    • Mẫu I: Chiều dày da $h_1 = h_2 = 2\text{ mm}$, chiều dày lõi PU $H = 20\text{ mm}$, $\rho_{\text{PU}} = 100\text{ kg/m}^3$.
    • Mẫu K: Chiều dày da $h_1 = h_2 = 2\text{ mm}$, chiều dày lõi PU $H = 20\text{ mm}$, $\rho_{\text{PU}} = 190\text{ kg/m}^3$.
    • Mẫu E, F, G, H: Thay đổi độc lập chiều dày lớp da ($1,5\text{ mm} - 3,5\text{ mm}$) và chiều dày lớp lõi ($10\text{ mm} - 40\text{ mm}$).

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

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

  1. Hiệu ứng nâng tần số cộng hưởng của liên kết ngàm (Clamped Boundary Effect): Kết quả giải tích và mô phỏng chỉ ra rằng ở dải tần số thấp ($f < 250\text{ Hz}$), điều kiện biên ngàm 4 cạnh (CCCC) làm tăng tổn thất truyền âm từ $4\text{ dB}$ đến $8\text{ dB}$ so với liên kết tựa bản lề (SSSS). Nguyên nhân là do liên kết ngàm gia tăng độ cứng uốn tổng thể của tấm, đẩy tần số dao động riêng cơ bản $f_{11}$ lên cao hơn xấp xỉ 1,41 lần, thu hẹp vùng truyền âm cộng hưởng bất lợi.

  2. Cơ chế kiểm soát truyền âm qua khoang không khí trong tấm kép: Đối với tấm kép composite lớp trực hướng có khoang khí, sự xuất hiện của tần số cộng hưởng khối lượng - không khí - khối lượng (Mass-Air-Mass resonance frequency $f_{\text{MAM}}$) gây ra sự sụt giảm nghiêm trọng giá trị STL (giảm từ $10 - 15\text{ dB}$ tại tần số $f_{\text{MAM}} \approx 80 - 120\text{ Hz}$). Tuy nhiên, khi vượt qua tần số này, đường cong STL tăng vọt với độ dốc đạt xấp xỉ $12\text{ dB/octave}$, vượt trội hoàn toàn so với độ dốc $6\text{ dB/octave}$ của tấm đơn cùng tổng khối lượng. Khi tăng chiều dày khoang khí từ $20\text{ mm}$ lên $80\text{ mm}$, tần số $f_{\text{MAM}}$ dịch chuyển về dải hạ âm ($< 50\text{ Hz}$), mở rộng vùng cách âm hiệu quả trong dải âm nghe thấy.

  3. Sự triệt tiêu hiện tượng sụt giảm trùng tần số nhờ lõi PU xốp đàn hồi Biot: Ở tấm kim loại hoặc composite đơn lớp đồng nhất, hiện tượng trùng tần số (Coincidence effect) xảy ra khi vận tốc sóng uốn của tấm bằng vận tốc sóng âm trong không khí, gây ra một "hố sâu" suy giảm STL từ $15 - 20\text{ dB}$ tại vùng tần số cao ($f = 2000 - 4000\text{ Hz}$). Trong tấm composite sandwich lõi PU foam, nhờ cơ chế tiêu tán năng lượng nội ma sát nhớt giữa khung rắn polymer và chất lưu không khí trong lỗ xốp (được mô hình hóa bởi 3 sóng Biot), hiện tượng sụt giảm STL bị triệt tiêu hoàn toàn. Giá trị STL duy trì đà tăng ổn định, đạt trên $42\text{ dB}$ tại tần số $4000\text{ Hz}$ (Mẫu K).

  4. Tính phụ thuộc định hướng âm học vào cấu hình sắp xếp lớp sợi composite (Anisotropic Acoustic Tailoring): Luận án phát hiện ra rằng tổn thất truyền âm phụ thuộc mạnh mẽ vào góc tới $\varphi$ và góc phương vị $\theta$ của sóng âm. Với góc tới nghiêng $\varphi = 60^\circ$, STL giảm từ $8 - 12\text{ dB}$ so với góc tới thẳng đứng $\varphi = 0^\circ$. Đặc biệt, khi thay đổi góc phương vị $\theta$ từ $0^\circ$ (song song trục sợi có $E_{11}$ lớn) đến $90^\circ$ (song song trục sợi có $E_{22}$ nhỏ), đường cong STL dịch chuyển rõ rệt. Việc tăng tỷ số dị hướng $E_{11}/E_{22}$ từ 1 lên 10 làm tăng STL dải tần số thấp thêm $3 - 5\text{ dB}$ do độ cứng uốn dọc trục tăng cường, mở ra khả năng thiết kế hướng sợi để tạo ra các vách ngăn cách âm định hướng.

  5. Tương quan định lượng giữa cơ tính lõi xốp và hiệu quả cách âm: So sánh thực nghiệm giữa Mẫu I ($\rho_c = 100\text{ kg/m}^3$) và Mẫu K ($\rho_c = 190\text{ kg/m}^3$) chỉ ra rằng việc tăng khối lượng riêng lớp lõi xốp PU giúp nâng cao STL trung bình toàn dải thêm từ $2,5\text{ dB}$ đến $4,8\text{ dB}$, đặc biệt là ở dải tần số trung bình $500 - 2000\text{ Hz}$, đồng thời cải thiện khả năng chống va đập và độ bền uốn của kết cấu vỏ tàu.

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                SO SÁNH ĐẶC TRƯNG STL GIỮA CÁC CẤU HÌNH KẾT CẤU (BẢN TỔNG HỢP)                   │
├───────────────────────┬─────────────────────────────┬───────────────────────────────────────────┤
│ Loại kết cấu          │ Vùng tần số thấp (< 400 Hz) │ Vùng tần số cao (> 1000 Hz)               │
├───────────────────────┼─────────────────────────────┼───────────────────────────────────────────┤
│ Tấm composite đơn     │ Kiểm soát bởi độ cứng uốn;  │ Xuất hiện hố suy giảm do trùng tần số     │
│                       │ Biên ngàm tăng 4-8 dB       │ ($STL$ giảm 10-15 dB tại $2-4\text{ kHz}$)│
├───────────────────────┼─────────────────────────────┼───────────────────────────────────────────┤
│ Tấm kép khoang khí    │ Xuất hiện đáy suy giảm      │ STL tăng cực nhanh (~12 dB/octave);       │
│                       │ tại tần số $f_{\text{MAM}}$ │ Vượt trội tấm đơn cùng khối lượng         │
├───────────────────────┼─────────────────────────────┼───────────────────────────────────────────┤
│ Composite Sandwich    │ Tối ưu nhờ độ cứng chống    │ Triệt tiêu hố trùng tần số nhờ lõi Biot;  │
│ lõi PU Foam           │ uốn cao của 2 lớp da        │ STL đạt đỉnh > 40-45 dB                   │
└───────────────────────┴─────────────────────────────┴───────────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết học thuật: Cung cấp hệ phương trình giải tích tổng quát và giải thuật số chuẩn xác cho bài toán truyền âm tương tác âm học - kết cấu composite trực hướng ba chiều, bổ khuyết khoảng trống lớn trong cơ học vật liệu và kết cấu tiên tiến.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp giữa tính toán ma trận giải tích phân tích Fourier/Galerkin và phương pháp đo vang âm thực nghiệm, có thể áp dụng cho các dạng kết cấu phức tạp hơn như vỏ trụ composite, tấm composite FGM hay tấm tổ ong (honeycomb).
  • Về mặt thực tiễn kỹ thuật công nghiệp: Cung cấp bộ công cụ phần mềm và cơ sở dữ liệu tính toán thiết kế trực tiếp cho kỹ sư đóng tàu vỏ composite (Viện Nghiên cứu Chế tạo Tàu thủy), kỹ sư thiết kế cabin ô tô tải và toa xe đường sắt, giúp tối ưu hóa chiều dày lớp da và tỷ trọng lõi xốp nhằm giảm thiểu trọng lượng bản thân trong khi vẫn đảm bảo tiêu chuẩn tiện nghi âm học quốc tế.
  • Về mặt chính sách và tiêu chuẩn: Tạo luận cứ khoa học để hoàn thiện quy chuẩn kỹ thuật quốc gia về kiểm soát tiếng ồn và rung động trên phương tiện thủy nội địa và tàu khách cao tốc bọc composite tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu (Boundary Conditions & Limitations):

  1. Giả thiết biến dạng nhỏ tuyến tính: Mô hình lý thuyết dựa trên lý thuyết tấm mỏng Kirchhoff, chỉ xét dao động uốn tuyến tính biên độ nhỏ, chưa tính đến biến dạng cắt phi tuyến bậc cao qua chiều dày đối với các tấm composite rất dày ($h/a > 0,1$) hoặc hiệu ứng hình học phi tuyến khi mức áp suất âm tới cực cao ($SPL > 140\text{ dB}$).
  2. Hình học phẳng chữ nhật: Đối tượng nghiên cứu mới dừng lại ở kết cấu tấm phẳng chữ nhật hữu hạn; chưa mở rộng sang các kết cấu vỏ cong (cylindrical shells, doubly-curved panels) vốn rất phổ biến trên thân máy bay hoặc vỏ tàu ngầm.
  3. Kích thích sóng phẳng điều hòa ổn định: Luận án khảo sát trường sóng âm phẳng điều hòa ở trạng thái dừng, chưa mô hình hóa đầy đủ trường kích thích ngẫu nhiên của lớp biên nhiễu loạn (Turbulent Boundary Layer - TBL) hay trường âm khuếch tán hoàn toàn ngẫu nhiên (Diffuse Sound Field).
  4. Bỏ qua ảnh hưởng của trường nhiệt độ và độ ẩm: Cơ tính của nhựa polymer và lõi xốp PU được giả định không đổi theo môi trường, chưa xét đến sự suy giảm liên kết và biến đổi tính chất âm học khi làm việc trong điều kiện nhiệt - ẩm khắc nghiệt kéo dài.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 5 định hướng cụ thể:

  1. Phát triển mô hình giải tích dựa trên Lý thuyết biến dạng cắt bậc cao (HSDT) hoặc Lý thuyết phần tử lớp rời rạc (Layerwise Theory) cho tấm composite dày và cấu trúc lưới không gian (lattice structures).
  2. Xây dựng bài toán truyền âm qua vỏ cong composite sandwich lõi xốp đàn nhớt dưới tác dụng của trường áp suất khí động học lớp biên nhiễu loạn (TBL).
  3. Ứng dụng thuật toán tối ưu hóa đa mục tiêu (Multi-Objective Genetic Algorithm - MOGA / PSO) để thiết kế tối ưu cấu hình lớp composite sandwich (chiều dày da, mật độ lõi, góc xếp sợi) nhằm tối đa hóa STL và giảm thiểu khối lượng kết cấu.
  4. Tích hợp vật liệu thông minh (vật liệu áp điện PZT, vật liệu biến từ magneto-rheological fluids) vào khoang lõi để thực hiện điều khiển âm thanh chủ động/bán chủ động thích nghi (Active-Passive Hybrid Acoustic Control).
  5. Nghiên cứu độ bền âm học mỏi (Acoustic Fatigue) và sự suy giảm tổn thất truyền âm theo thời gian của kết cấu sandwich composite trong môi trường biển nhiệt đới.

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

  • Tác động học thuật (Academic Impact): Tạo nền tảng toán học và cơ học vững chắc cho hướng nghiên cứu âm học kết cấu composite tại Việt Nam. Kết quả luận án mở đường cho hàng loạt công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Q1/Q2) trong các lĩnh vực Applied Acoustics, Journal of Sound and Vibration, và Composite Structures.
  • Chuyển giao công nghệ và công nghiệp đóng tàu: Ứng dụng trực tiếp tại Viện Nghiên cứu và Chế tạo Tàu thủy (Đại học Nha Trang) trong thiết kế vách cabin, nắp buồng máy cho các dòng tàu khách composite cao tốc, tàu cá xa bờ, giúp giảm mức ồn trong khoang hành khách từ $6 - 10\text{ dBA}$, đáp ứng tiêu chuẩn tiếng ồn của Đăng kiểm Việt Nam và Tổ chức Hàng hải Quốc tế (IMO).
  • Công nghiệp ô tô và giao thông đường sắt: Cung cấp thông số thiết kế vách ngăn composite nhẹ cách âm cho các nhà sản xuất thùng xe tải, xe buýt điện và toa xe đường sắt cao tốc, góp phần cắt giảm nhiên liệu/năng lượng vận hành và giảm phát thải khí nhà kính nhờ giảm tải trọng bản thân xe.
  • Ý nghĩa xã hội: Nâng cao chất lượng môi trường làm việc cho thuyền viên, công nhân và cải thiện mức độ tiện nghi, sức khỏe tâm - thính giác của hành khách trên các phương tiện giao thông thế hệ mới.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ học, Âm học, Kỹ thuật Hàng hải, Hàng không: Tiếp cận được một khung lý thuyết mẫu mực về ghép nối cơ học kết cấu dị hướng với âm học môi trường xốp và giải thuật số MATLAB đã được thẩm định.
  • Các giáo sư, nhà khoa học đầu ngành: Sở hữu nguồn tài liệu tham khảo chuyên sâu, tường minh về mặt toán học và vật lý để phát triển các hướng nghiên cứu chuyên sâu tiếp theo.
  • Kỹ sư R&D tại các tập đoàn công nghiệp chế tạo (Tàu thủy, Ô tô, Hàng không): Nắm bắt công cụ tính toán nhanh, chính xác tổn thất truyền âm $STL$ mà không cần phụ thuộc vào các phần mềm thương mại đắt tiền, rút ngắn chu kỳ thử nghiệm - thiết kế mẫu.
  • Cơ quan quản lý nhà nước và Đăng kiểm: Có cơ sở khoa học để xây dựng bộ tiêu chuẩn, quy chuẩn kiểm định âm học kết cấu vật liệu mới trên phương tiện giao thông vận tải thủy và bộ.

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 độc đáo nhất là việc thiết lập thành công mô hình giải tích tường minh mô tả tương tác âm học - kết cấu hoàn chỉnh cho tấm composite sandwich trực hướng hữu hạn có lớp lõi là môi trường xốp Biot dưới các điều kiện biên thực tế (ngàm và tựa đơn). Luận án đã mở rộng đồng thời hai lý thuyết nền tảng: Lý thuyết tấm mỏng cổ điển Kirchhoff (CLPT) cho vật liệu trực hướng nhiều lớp và Lý thuyết môi trường xốp Biot (1956) cho vật liệu hai pha rắn - lỏng. Điểm then chốt là giải quyết được bài toán ghép nối biên phức tạp giữa chuyển vị 2D của tấm composite mặt ngoài và trường sóng 3D (2 sóng dọc $k_1, k_2$ và 1 sóng cắt) trong lớp lõi xốp PU hữu hạn bằng chuỗi Fourier kép và phương pháp Galerkin.

2. Sự đổi mới về phương pháp luận của luận án khi so sánh với ít nhất 2 công trình quốc tế kinh điển?

  • So với Roussos (1984): Roussos chỉ giải được bài toán tấm đơn đẳng hướng và trực hướng chịu điều kiện biên tựa bản lề bốn cạnh (SSSS) bằng tích phân Green. Luận án đã đổi mới phương pháp luận khi xây dựng ma trận trở kháng bức xạ kết cấu tổng quát, giải quyết thành công điều kiện biên ngàm bốn cạnh (CCCC) bằng hệ hàm dạng chùm tia rung (beam vibration functions), đồng thời mở rộng từ tấm đơn sang tấm kép có khoang khí và tấm sandwich lõi Biot.
  • So với Bolton, Shiau và Kang (1996): Bolton và cộng sự sử dụng phương pháp ma trận truyền (TMM) thuần túy cho kết cấu sandwich nhôm vô hạn. Luận án đã đổi mới phương pháp luận khi tích hợp lý thuyết sóng Biot vào kết cấu hữu hạn có xét đến tính trực hướng của sợi composite và ảnh hưởng của kích thước hình học biên, khắc phục hoàn toàn nhược điểm bỏ qua mode dao động uốn kết cấu cục bộ của Bolton.

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

Phát hiện bất ngờ nhất là hiện tượng tăng độ dày lớp lõi xốp PU không phải lúc nào cũng làm tăng $STL$ một cách tuyến tính ở mọi dải tần số. Tại một số dải tần số trung bình hẹp ($400 - 800\text{ Hz}$), việc tăng độ dày lõi xốp từ $20\text{ mm}$ lên $40\text{ mm}$ làm tăng tính đàn hồi của hệ thống, vô tình kích hoạt các mode biến dạng cắt cục bộ trong lớp lõi, dẫn đến sự suy giảm nhẹ giá trị STL cục bộ trước khi bứt phá mạnh mẽ ở dải tần cao. Phát hiện này cung cấp cảnh báo then chốt cho các kỹ sư thiết kế: không thể tối ưu hóa cách âm chỉ bằng cách tăng kích thước chiều dày mà phải điều chỉnh đồng thời tỷ trọng xốp và cấu hình định hướng góc sợi của lớp da composite.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) rõ ràng không?

Có, luận án cung cấp quy trình thực nghiệm cực kỳ chi tiết và minh bạch:

  • Thông số buồng đo: Phòng thu có thể tích chính xác $V_2 = 16\text{ m}^3$, bề mặt ốp gạch men phản xạ cao; hệ thống giá gá mẫu thử có kích thước tiêu chuẩn được kẹp chặt bằng khung thép bắt bu lông lực siết đồng đều có đệm cao su giảm chấn.
  • Thiết bị đo lường: Nguồn phát âm công suất lớn phát tín hiệu tiếng ồn trắng (White noise) hoặc tạp âm hồng (Pink noise) phân bố đều dải tần $20\text{ Hz} - 10.000\text{ Hz}$; hệ thống micro chuẩn âm học độ nhạy cao kết nối với bộ phân tích phổ thời gian thực đa kênh; thiết bị đo thời gian vang $T_{60}$ theo phương pháp tích phân phân rã ngược Schroeder.
  • Quy trình chế tạo mẫu: Mô tả chi tiết tỷ lệ thành phần nhựa polyester không no, chất đóng rắn MEKP, vải sợi thủy tinh dệt (Woven Roving) và quy trình đổ xốp PU foam trực tiếp đảm bảo độ kết dính bề mặt hoàn hảo không có bọt khí phân lớp.

5. Luận án đã vạch ra chương trình nghiên cứu 10 năm (10-year research agenda) như thế nào?

Chương trình nghiên cứu 10 năm được định hình qua 3 giai đoạn chiến lược:

  • Giai đoạn 1 (Năm 1-3): Mở rộng lý thuyết giải tích sang kết cấu vỏ trụ composite sandwich và tấm composite có độ dày lớn (sử dụng lý thuyết biến dạng cắt bậc cao HSDT và mô hình đàn nhớt Fractional Derivative) chịu tải trọng âm khí động học (Aeroacoustic TBL).
  • Giai đoạn 2 (Năm 4-6): Phát triển vật liệu composite đa chức năng tích hợp cấu trúc vi mô âm học (Acoustic Metamaterials), lõi mạng lưới in 3D (Auxetic/Lattice Cores) kết hợp bọt xốp nano để tạo ra các dải cấm âm (Bandgaps) triệt tiêu hoàn toàn tiếng ồn dải tần cực thấp ($< 100\text{ Hz}$).
  • Giai đoạn 3 (Năm 7-10): Tích hợp công nghệ AI và Tối ưu hóa Topology để tự động hóa thiết kế vỏ phương tiện giao thông thông minh; thương mại hóa bộ phần mềm tính toán âm học kết cấu và chuyển giao công nghệ cho ngành công nghiệp hàng hải, hàng không vũ trụ quốc gia.

Kết luận

  1. Thiết lập thành công hệ phương trình giải tích tổng quát: Xây dựng hệ phương trình dao động âm tương tác hoàn chỉnh cho ba hệ kết cấu: tấm composite lớp trực hướng đơn, tấm kép composite chứa khoang khí và tấm composite sandwich lõi xốp PU hữu hạn, chịu các điều kiện biên ngàm (CCCC) và tựa đơn (SSSS) dưới kích thích sóng âm phẳng điều hòa ba chiều.
  2. Hợp nhất đột phá giữa lý thuyết Biot và cơ học composite: Tích hợp thành công lý thuyết môi trường xốp 3D Biot (xét 2 sóng nén $k_1, k_2$ và 1 sóng cắt) với lý thuyết tấm composite trực hướng hữu hạn mỏng, giải quyết triệt để sự tương thích biên ứng suất - chuyển vị tại các mặt phân cách vật liệu.
  3. Phát triển bộ công cụ tính toán số tin cậy trên MATLAB: Xây dựng thuật toán và chương trình tính toán số tối ưu dựa trên chuỗi Fourier kép và phương pháp Galerkin, cho phép dự đoán chính xác tổn thất truyền âm $STL$ và hệ số bức xạ âm với thời gian xử lý nhanh và độ ổn định số học cao.
  4. Kiểm chứng thực nghiệm đạt độ chính xác cao: Thực hiện thành công các thí nghiệm đo $STL$ trong hệ thống phòng phát - phòng thu tiêu chuẩn trên các mẫu tấm thép, nhôm và các mẫu sandwich GFRP/PU foam tự chế tạo; đối sánh khẳng định sai số giữa lý thuyết và thực nghiệm duy trì dưới $3\text{ dB}$ trên toàn dải tần số $20\text{ Hz} - 10.000\text{ Hz}$.
  5. Định lượng các quy luật cơ học - âm học then chốt: Khẳng định điều kiện biên ngàm giúp tăng $STL$ dải tần thấp thêm $4 - 8\text{ dB}$; làm sáng tỏ cơ chế dịch chuyển tần số cộng hưởng khoang khí $f_{\text{MAM}}$ trong tấm kép; chứng minh khả năng triệt tiêu hố sụt giảm trùng tần số của lõi xốp PU; và phát hiện tính định hướng âm học phụ thuộc góc nghiêng/phương vị sợi composite.
  6. Giá trị ứng dụng thực tiễn và chuyển giao công nghệ vững chắc: Cung cấp luận cứ khoa học và công cụ tính toán thiết kế trực tiếp phục vụ chế tạo vỏ tàu thủy composite cách âm, cabin ô tô và toa xe đường sắt tại Việt Nam, đóng góp thiết thực cho sự phát triển của ngành cơ học ứng dụng và công nghiệp chế tạo phương tiện giao thông hiện đại.