Tổng quan về luận án

Kiểm tra sức khỏe công trình (Structural Health Monitoring - SHM) và chẩn đoán hư hỏng kết cấu đóng vai trò quyết định trong việc bảo đảm an toàn vận hành, phòng ngừa sự cố thảm khốc và tối ưu hóa vòng đời khai thác của các cơ sở hạ tầng trọng yếu như nhà máy điện hạt nhân, công trình biển, hàng không vũ trụ và cầu nhịp lớn. Trong bối cảnh vật liệu composite tiên tiến ngày càng được ứng dụng rộng rãi, kết cấu làm từ vật liệu có cơ lý tính biến thiên liên tục (Functionally Graded Materials - FGM) nổi lên như một giải pháp ưu việt giúp loại bỏ ứng suất tập trung tại mặt phân giới. Tuy nhiên, tính chất bất đồng nhất và bất đối xứng theo chiều dày của dầm FGM dẫn đến hiện tượng ghép cặp cơ học phức tạp giữa biến dạng uốn và biến dạng dọc trục, khiến việc phát hiện khuyết tật nứt trở nên vô cùng khó khăn đối với các kỹ thuật SHM truyền thống.

Khoảng trống nghiên cứu (research gap) then chốt nằm ở chỗ: các phương pháp chẩn đoán vết nứt dựa thuần túy trên sự suy giảm tần số dao động riêng đối với dầm FGM thường gặp bế tắc tại các "điểm nút tần số" (frequency nodal points) – nơi vết nứt xuất hiện nhưng không làm biến đổi tần số của mode dao động tương ứng – đồng thời không thể phân định vị trí hư hỏng khi kết cấu có điều kiện biên đối xứng. Bên cạnh đó, các phương pháp số như Phần tử hữu hạn (FEM) bộc lộ sự hạn chế nghiêm trọng về bậc tự do và độ chính xác khi mô phỏng ở dải tần số cao.

Luận án tiến sĩ "Ứng dụng vật liệu áp điện trong đánh giá trạng thái kỹ thuật công trình" của tác giả Lưu Quỳnh Hường (Chuyên ngành: Cơ kỹ thuật, Mã số: 9 52 01 01; Người hướng dẫn: GS. Nguyễn Tiến Khiêm, TS. Trần Thanh Hải; Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam, 2021) đã giải quyết triệt để bài toán này. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  • RQ1: Làm thế nào để thiết lập một mô hình giải tích mở rộng giải quyết chính xác bài toán ghép cặp cơ - điện trong dầm FGM Timoshenko có gắn lớp/miếng áp điện ở mọi dải tần số?
  • RQ2: Ảnh hưởng của vị trí, kích thước miếng áp điện và quy luật biến thiên vật liệu (chỉ số gradient $n$) đến các tham số động lực học là gì?
  • RQ3: Sự xuất hiện của vết nứt làm thay đổi trường cơ điện như thế nào và làm sao tận dụng điện tích cảm biến dao động (Modal Sensor Output Charge - MPC) để xây dựng thuật toán nhận dạng hư hỏng đơn định?
  • H1: Phương pháp độ cứng động lực (Dynamic Stiffness Method - DSM) có thể khắc phục triệt để sự suy giảm độ chính xác ở tần số cao của phương pháp PTHH đối với kết cấu FGM áp điện có bậc.
  • H2: Việc kết hợp đường đồng mức tần số dao động riêng với đường đồng mức điện tích cảm biến áp điện phân bố cho phép xác định chính xác, duy nhất tọa độ và độ sâu vết nứt, loại bỏ hoàn toàn điểm mù tại các nút dao động.

Nghiên cứu tích hợp lý thuyết dầm cắt bậc cao Timoshenko, nguyên lý biến phân Hamilton, cơ học phá hủy đàn hồi tuyến tính và lý thuyết áp điện phân bố. Luận án đạt đóng góp đột phá khi phát triển thành công phương pháp đường đồng mức điện tích cảm biến dao động, nâng cao độ chính xác chẩn đoán vết nứt với sai số vị trí dưới 1% ngay cả trong điều kiện biên đối xứng phức tạp.

Literature Review và Positioning

Nghiên cứu về nhận dạng hư hỏng kết cấu khởi đầu từ các công trình kinh điển của Cawley và Adams (1979), Salawu (1997), Pandey và Samman (1994) trên kết cấu đồng nhất (homogeneous). Đối với kết cấu composite và FGM, các mô hình dao động được phát triển bởi Hu và Zhang (2003), Li (2008), Larbi và cộng sự (2013), Su và Banerjee (2015). Nghiên cứu về dao động dầm FGM có vết nứt sau đó được mở rộng bởi Yang và Chen (2008), Aydin (2013), Khiem và Lien (2001).

Trong lĩnh vực vật liệu thông minh, hiệu ứng áp điện phát hiện bởi Pierre và Jacques Curie (1880) đã được Crawley và de Luis (1987) phát triển thành lý thuyết kích hoạt áp điện phân bố (distributed actuator). Lee và Moon (1990) xây dựng lý thuyết cảm biến và bộ kích hoạt dạng mode (modal sensor/actuator), chứng minh khả năng đo đạc và kích thích các mode dao động riêng biệt của kết cấu dầm tấm.

                      ┌────────────────────────────────────────┐
                      │    TIẾN TRÌNH LÝ THUYẾT & POSITIONING  │
                      └──────────────────┬─────────────────────┘
                                         │
       ┌─────────────────────────────────┴─────────────────────────────────┐
       ▼                                                                   ▼
┌───────────────────────────────┐                       ┌─────────────────────────────────────┐
│  Chẩn đoán hư hỏng dầm FGM    │                       │     Kết cấu & Cảm biến áp điện      │
│  - Yang & Chen (2008)         │                       │  - Crawley & de Luis (1987)         │
│  - Yu & Chu (2009): p-FEM     │                       │  - Lee & Moon (1990): Modal sensor  │
│  - Banerjee et al. (2014): GA │                       │  - Park et al. (2003): EMI dải cao  │
│  - Khiem & Huyen (2017)       │                       │  - Sunar & Rao (1999)               │
└──────────────┬────────────────┘                       └──────────────────┬──────────────────┘
               │                                                           │
               │   ┌───────────────────────────────────────────────────┐   │
               └──►│ BẾ TẮC: Điểm nút tần số & Đa nghiệm biên đối xứng │◄──┘
                   └─────────────────────────┬─────────────────────────┘
                                             │
                                             ▼
                   ┌───────────────────────────────────────────────────┐
                   │    POSITIONING ĐỘT PHÁ CỦA LUẬN ÁN (2021)         │
                   │    - Dynamic Stiffness Method (DSM) cho FGM       │
                   │    - True Neutral Axis Shift: h0 = f(n, re)       │
                   │    - Modal Sensor Charge Contour Method           │
                   └───────────────────────────────────────────────────┘

Tranh luận học thuật cốt lõi tồn tại giữa hai trường phái:

  1. Trường phái chẩn đoán dựa trên trở kháng cơ điện (EMI) (Giurgiutiu & Rogers, 2000; Park et al., 2003; Na & Baek, 2018): Khai thác đáp ứng cục bộ tần số cao (130–400 kHz) thông qua chỉ số RMSD và MAPD. Tuy nhiên, phương pháp này có bán kính cảm nhận hạn chế (0.4m đối với composite, 2m đối với kim loại), tín hiệu hư hỏng dễ bị che lấp bởi nhiễu điện từ trường và biến thiên nhiệt độ môi trường.
  2. Trường phái tham số dao động toàn cục (Modal parameters) (Zhao et al., 2015; Jiang et al., 2017): Tận dụng điện áp cảm biến để tăng độ nhạy của tần số riêng. Dẫu vậy, phương pháp vẫn vấp phải rào cản khi vết nứt nằm ở xa đầu dò hoặc nằm ngay tại các điểm nút tần số (nodal points).

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

  • Yu và Chu (2009): Sử dụng phương pháp phần tử hữu hạn mở rộng bậc cao ($p$-version FEM) kết hợp đường đồng mức tần số trên dầm FGM công-xôn. Hạn chế: khi có sai số đo đạc thực nghiệm 1–2%, sai số chẩn đoán vị trí vết nứt bị khuếch đại lên tới 15%.
  • Banerjee và cộng sự (2014): Sử dụng mô hình PTHH 3 chiều kết hợp thuật toán di truyền (Genetic Algorithm - GA) để dò tìm vết nứt. Dù độ chính xác được cải thiện (sai số cực đại 6.4%), phương pháp đòi hỏi chi phí tính toán cực lớn và không thiết lập được biểu thức giải tích tường minh cho trường ghép cặp điện - cơ.

Luận án của Lưu Quỳnh Hường định vị đột phá bằng cách tích hợp trực tiếp lớp cảm biến áp điện phân bố phủ kín chiều dài dầm FGM, xây dựng mô hình ma trận độ cứng động (DSM) chính xác tuyệt đối ở mọi dải tần số, kết hợp chỉ số điện tích cảm biến $Q$ nhằm loại bỏ hoàn toàn các điểm mù chẩn đoán mà Yu & Chu (2009) cũng như Khiem & Huyên (2017) chưa thể giải quyế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 và hoàn thiện hệ thống lý thuyết dao động kết cấu thông minh thông qua các đóng góp căn bản:

  • Mở rộng lý thuyết dầm cắt Timoshenko cho vật liệu FGM bất đối xứng: Xác định chính xác độ dịch chuyển của trục trung hòa thực $h_0$, khắc phục giả định giản đơn coi trục trung hòa trùng với trục hình học ở các nghiên cứu tiền đề.
  • Thiết lập hệ phương trình vi phân dao động ghép cặp điện - cơ: Áp dụng nguyên lý Hamilton dẫn xuất tường minh quan hệ giữa chuyển vị dọc trục $u_0(x,t)$, góc xoay tiết diện $\theta(x,t)$, độ võng $w_0(x,t)$ và chuyển vị điện $D(x,t)$ của lớp áp điện.
  • Xây dựng mô hình tán xạ vết nứt kép (Axial - Rotational Spring Model): Mô tả vết nứt hở trong dầm FGM bằng cặp lò xo đàn hồi cục bộ (lò xo dọc trục độ cứng $K_t$ và lò xo xoay độ cứng $K_r$), tính toán thông qua giải tích hàm giải phóng năng lượng biến dạng (Strain Energy Release Rate) và tích phân biến dạng đàn hồi tuyến tính.

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

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết: (1) Cơ học vật liệu biến thiên cơ lý tính FGM theo hàm lũy thừa; (2) Điện từ học và nhiệt đàn hồi áp điện phân bố PZT/PVDF; (3) Lý thuyết sóng động lực học và độ cứng động lực (DSM).

Mô hình vật liệu FGM biến thiên theo chiều dày được mô tả bởi quy luật: $$\Re(z) = \Re_b + (\Re_t - \Re_b)\left(\frac{z}{h} + 0.5\right)^n$$ trong đó $\Re \in {E, G, \rho}$, chỉ số $t$ và $b$ tương ứng với lớp bề mặt trên (top - gốm) và dưới (bottom - kim loại); $n$ là chỉ số gradient phân bố thể tích ($n \ge 0$).

Vị trí trục trung hòa thực $h_0$ được chứng minh có biểu thức giải tích chặt chẽ: $$h_0 = \frac{n(r_e - 1)h}{2(n + 2)(n + r_e)}; \quad r_e = \frac{E_t}{E_b}$$

Phương trình cấu hoán áp điện của lớp cảm biến PZT: $$\sigma_{px} = C_{11}^p \varepsilon_{px} - h_{13} D; \quad \tau_p = C_{55}^p \gamma_p; \quad \in = -h_{13} \varepsilon_{px} + \beta_{33}^p D$$ trong đó $C_{11}^p, h_{13}, \beta_{33}^p$ lần lượt là mô đun đàn hồi, hằng số áp điện và hằng số điện môi; $D$ là thông lượng điện dịch chuyển dọc trục phân cực.

Biểu thức xác định điện tích đầu ra $Q$ của lớp cảm biến phân bố trên toàn bộ chiều dài dầm: $$Q = \int_A D dA = b \int_0^L D dx = \frac{b h_{13}}{\beta_{33}^p} \int_0^L \left(u_0' - \frac{h}{2}\theta'\right) dx$$

Điều kiện biên và tính tương thích tiếp xúc hoàn toàn (không trượt, không tách lớp) giữa dầm FGM chủ và lớp áp điện: $$u\left(x, \frac{h_b}{2}, t\right) = u_p\left(x, -\frac{h_p}{2}, t\right); \quad w\left(x, \frac{h_b}{2}, t\right) = w_p\left(x, -\frac{h_p}{2}, t\right)$$

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

Thiết kế nghiên cứu

Luận án quán triệt thế giới quan thực chứng duy lý kết hợp chủ nghĩa hiện thực phản biện (Critical Realism) trong cơ học công trình: phát triển mô hình toán học giải tích mở rộng (Semi-analytical Dynamic Stiffness Formulation), sau đó kiểm chứng chéo thông qua mô phỏng số học trường dữ liệu động lực học đa chiều trên môi trường MATLAB.

Thiết kế nghiên cứu đa tầng bao gồm:

  1. Tầng 1 (Phần tử vi phân cơ điện): Dẫn xuất hệ phương trình vi phân đạo hàm riêng cấp 2 liên kết theo không gian và thời gian trong miền tần số.
  2. Tầng 2 (Phần tử dầm bậc FGM áp điện DSM): Thành lập ma trận độ cứng động chính xác $[D_e(\omega)]$, ghép nối kết cấu nhiều bậc theo thuật toán phần tử phổ.
  3. Tầng 3 (Bài toán ngược chẩn đoán hư hỏng): Xây dựng không gian nghiệm 2 tham số (vị trí vết nứt $e/L$ và độ sâu vết nứt $a/h$) bằng phương pháp giao thoa đường đồng mức đa biến (Multi-contour intersections).
 ┌─────────────────────────────────────────────────────────────────────────┐
 │               QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU TOÀN DIỆN                │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
 ┌────────────────────────────────────┴────────────────────────────────────┐
 │ BƯỚC 1: MÔ HÌNH HÓA ĐỘ CỨNG ĐỘNG LỰC (DSM)                              │
 │ - Lý thuyết dầm Timoshenko FGM + Neutral axis shift h0                  │
 │ - Ghép cặp liên tục dầm chủ - Lớp áp điện PZT                           │
 │ - Nguyên lý biến phân Hamilton -> Hệ PTPT vi phân miền tần số           │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
 ┌────────────────────────────────────┴────────────────────────────────────┐
 │ BƯỚC 2: THIẾT LẬP NGHIỆM GIẢI TÍCH & MA TRẬN PHẦN TỬ                    │
 │ - Phương trình đặc trưng bậc 3: det[λ²A + λB + C] = 0                   │
 │ - 6 hàm dạng chính xác: {Z0(x,ω)} = [G0(x,ω)]{C}                        │
 │ - Ma trận độ cứng động phần tử: [De(ω)] = [Q(ω)][Z(ω)]⁻¹                │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
 ┌────────────────────────────────────┴────────────────────────────────────┐
 │ BƯỚC 3: MÔ HÌNH HÓA VẾT NỨT & BÀI TOÁN THUẬN                            │
 │ - Cặp lò xo tương đương (Kt, Kr) tại vị trí vết nứt e                   │
 │ - Ghép ma trận toàn hệ: [D(ω)]{U} = {P} -> det[D(ω)] = 0                │
 │ - Xuất trường dao động: Tần số riêng ω_i & Điện tích cảm biến Q_i       │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
 ┌────────────────────────────────────┴────────────────────────────────────┐
 │ BƯỚC 4: THUẬT TOÁN BÀI TOÁN NGƯỢC CHẨN ĐOÁN HƯ HỎNG                     │
 │ - Xây dựng cơ sở dữ liệu đường đồng mức tần số f(e, a)                  │
 │ - Xây dựng cơ sở dữ liệu đường đồng mức điện tích cảm biến Q(e, a)      │
 │ - Giao điểm đường đồng mức đa mode -> Xác định DUY NHẤT (e, a)          │
 └─────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thiết lập ma trận độ cứng động lực được triển khai qua các bước:

  1. Áp dụng biến đổi Fourier đưa hệ phương trình chuyển động về dạng ma trận trong miền tần số: $$[A]{Z''(x,\omega)} + [B]{Z'(x,\omega)} + [C]{Z(x,\omega)} = 0$$ với ${Z(x,\omega)} = {u_0, \theta, w_0}^T$.
  2. Tìm nghiệm dạng ${Z_0} = {d}e^{\lambda x}$, dẫn đến phương trình đặc trưng bậc 3 đối với $\eta = \lambda^2$: $$\det[\lambda^2 A + \lambda B + C] = 0$$ cho ra 6 nghiệm riêng $\lambda_1, \dots, \lambda_6$ và thiết lập ma trận nghiệm tổng quát: $${Z_0(x,\omega)} = [G_0(x,\omega)]{C}$$
  3. Thiết lập quan hệ giữa vectơ chuyển vị nút ${U_e(\omega)}$ và lực nút ${P_e(\omega)}$: $${U_e(\omega)} = [\mathbb{Z}(\omega)]{C}; \quad {P_e(\omega)} = [\mathbb{Q}(\omega)]{C}$$ Từ đó triệt tiêu vectơ hằng số tích phân ${C}$, thu được ma trận độ cứng động lực phần tử: $$[D_e(\omega)] = [\mathbb{Q}(\omega)][\mathbb{Z}(\omega)]^{-1}$$
  4. Ma trận độ cứng động toàn kết cấu được ghép nối chính xác: $$[D(\omega)] = \sum_{e=1}^{N_e} [T_e]^T [D_e(\omega)][T_e]$$ Phương trình tần số dao động tự do không suy biến: $\det[D(\omega)] = 0$.

Data và phân tích

Thông số vật liệu kiểm chuẩn trong nghiên cứu:

  • Vật liệu FGM (Alumina/Thép): $E_t = 390\text{ GPa}$, $\rho_t = 3960\text{ kg/m}^3$, $\mu_t = 0.25$; $E_b = 210\text{ GPa}$, $\rho_b = 7800\text{ kg/m}^3$, $\mu_b = 0.31$.
  • Lớp áp điện PZT-4: $\rho_p = 7750\text{ kg/m}^3$, $C_{11}^p = 139\text{ GPa}$, $C_{55}^p = 25.6\text{ GPa}$, $h_{13} = -7.0 \times 10^8\text{ V/m}$, $\beta_{33}^p = 1.15 \times 10^8\text{ m/F}$.
  • Dải khảo sát hình học: Tỷ số độ mảnh $L_b/h_b \in {10, 20, 30, 50, 100}$; tỷ số chiều dày áp điện $h_p/h_b \in [0, 0.5]$; tỷ số chiều dài miếng vá áp điện $L_p/L_b \in [1/6, 1.0]$; chỉ số gradient vật liệu $n \in {0.2, 0.5, 1, 2, 5}$.
  • Điều kiện biên: Tựa đơn hai đầu (SS), Ngàm - Ngàm (CC), Ngàm - Tự do (CF).

Độ tin cậy của thuật toán DSM được kiểm chuẩn đối sánh trực tiếp với kết quả công bố của Su và Banerjee (2015) trên dầm FGM nguyên vẹn không áp điện ($h_p = 0$).

Điều kiện biên Mode Su & Banerjee (2015) Luận án (DSM, $h_p=0$) Độ lệch tương đối (%)
Dầm SS ($L/h=10, n=2$) Mode 1 2.5714 2.5711 0.012%
Mode 2 9.3871 9.3860 0.012%
Mode 3 18.946 18.945 0.005%
Dầm CC ($L/h=10, n=2$) Mode 1 5.2536 5.2530 0.011%
Mode 2 12.873 12.871 0.016%
Mode 3 22.842 22.840 0.009%
Dầm CF ($L/h=10, n=2$) Mode 1 0.9782 0.9780 0.020%
Mode 2 5.2341 5.2338 0.006%
Mode 3 12.981 12.980 0.008%

Sự khác biệt cực nhỏ ($\le 0.02%$) xuất phát từ việc Su & Banerjee bỏ qua độ dịch chuyển trục trung hòa thực $h_0$, trong khi luận án tính toán chính xác hình học biến dạng bất đối xứng.

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

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

Nghiên cứu mang lại 5 phát hiện cốt lõi có giá trị khoa học vượt bậc:

┌─────────────────────────────────────────────────────────────────────────┐
│                    5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                      │
├─────────────────────────────────────────────────────────────────────────┤
│ 1. Tần số chuẩn hóa biến thiên dạng Parabol theo chiều dày áp điện      │
│    -> Tồn tại điểm chuyển tiếp tối ưu giữa độ cứng uốn và khối lượng.   │
├─────────────────────────────────────────────────────────────────────────┤
│ 2. Tần số riêng dầm trơn giảm nhanh hơn dầm có áp điện khi n tăng       │
│    -> Lớp áp điện PZT đóng vai trò ổn định và bù trừ độ cứng suy giảm.  │
├─────────────────────────────────────────────────────────────────────────┤
│ 3. Vị trí miếng dán tạo hiệu ứng đảo chiều động lực học                 │
│    -> Dán ở ngàm: tần số tăng đơn điệu; dán ở giữa/tự do: tần số giảm.  │
├─────────────────────────────────────────────────────────────────────────┤
│ 4. Điện tích cảm biến Q nhạy cảm vượt trội với vết nứt so với tần số ω  │
│    -> Đạo hàm dạng riêng (góc xoay) tại vết nứt kích phát đột biến Q.   │
├─────────────────────────────────────────────────────────────────────────┤
│ 5. Giao thoa đường đồng mức Q loại bỏ 100% điểm mù tại các nút tần số   │
│    -> Phân giải đơn trị tuyệt đối bài toán ngược với biên đối xứng.    │
└─────────────────────────────────────────────────────────────────────────┘
  1. Quy luật biến thiên parabol của tần số chuẩn hóa theo bề dày áp điện: Tần số chuẩn hóa (tỷ số giữa tần số riêng của dầm có áp điện và dầm nguyên vẹn) biến thiên phi tuyến dạng parabol đối với hầu hết các điều kiện biên. Khi tăng $h_p/h_b$, sự gia tăng độ cứng uốn lấn át tải trọng bản thân ở giai đoạn đầu, nhưng hiệu ứng quán tính khối lượng áp điện ($\rho_p = 7750\text{ kg/m}^3$) sẽ chiếm ưu thế ở chiều dày lớn, khiến tần số đạt cực trị cục bộ.
  2. Hiệu ứng bù trừ suy giảm độ cứng theo chỉ số gradient: Khi chỉ số gradient $n$ tăng từ 0 lên 5 (hàm lượng kim loại tăng lên, mô đun đàn hồi giảm), tần số riêng của dầm giảm mạnh. Tuy nhiên, tần số riêng của dầm có dán lớp áp điện giảm chậm hơn đáng kể so với dầm trần, khẳng định lớp phủ áp điện đóng vai trò gia cường thụ động và ổn định hóa động lực học cho kết cấu FGM.
  3. Sự phụ thuộc phi đối xứng theo vị trí lắp đặt miếng áp điện:
    • Dán tại vị trí ngàm (vùng mô men uốn cực đại): Tần số chuẩn hóa bậc 1 của dầm CC và CF tăng đơn điệu theo độ dày miếng dán.
    • Dán tại giữa nhịp hoặc đầu tự do: Tần số chuẩn hóa bậc 1 và bậc 2 giảm đơn điệu do miếng dán chủ yếu đóng vai trò gia tải khối lượng tập trung thay vì đóng góp độ cứng uốn.
  4. Độ nhạy vượt trội của điện tích cảm biến dao động ($Q$) trước sự xuất hiện vết nứt: Khác với tần số riêng dao động $\omega$ (vốn là đại lượng tích phân toàn cục ít nhạy với khuyết tật nhỏ), điện tích $Q$ tích phân trực tiếp từ góc xoay $\theta(x)$ và biến dạng dọc trục $u_0'(x)$. Tại vị trí vết nứt hở, góc xoay mặt cắt ngang chịu bước nhảy đột ngột, làm thay đổi mạnh mẽ mật độ điện tích cảm ứng bề mặt.
  5. Khắc phục triệt để điểm mù tại các nút tần số và tính đa nghiệm biên đối xứng:
    • Khi vết nứt nằm tại nút dao động của mode $i$ (nơi $\Delta\omega_i = 0$), đường đồng mức tần số của mode đó suy biến thành một dải vô định. Ngược lại, điện tích cảm biến dao động $Q_i$ vẫn biến đổi rõ rệt do độ cong và góc xoay tại nút không triệt tiêu đồng thời.
    • Thử nghiệm chẩn đoán số với vết nứt tại $e = L/3$ và độ sâu $a/h = 30%$ trên dầm SS, CC và CF cho thấy: việc giao thoa giữa các đường đồng mức điện tích cảm biến 3 mode đầu xác định chính xác giao điểm duy nhất $(e=0.333L, a=0.30h)$ với sai số xác định vị trí $\le 0.85%$ và sai số độ sâu $\le 1.2%$.

Implications đa chiều

  • Lý thuyết: Đặt nền móng giải tích chính xác cho cơ học tiếp xúc điện - đàn hồi trên kết cấu bất đồng nhất FGM chịu hư hỏng cục bộ.
  • Phương pháp luận: Mở rộng năng lực của DSM sang lĩnh vực SHM, thay thế quy trình mô phỏng số nặng nề bằng thuật toán phần tử phổ tốc độ cao.
  • Thực tiễn & Công nghệ: Cung cấp nguyên lý thiết kế cảm biến phân bố thông minh trên cánh tuabin điện gió, cánh máy bay composite và ống dẫn hạt nhân, cho phép tự động hóa hoàn toàn quá trình giám sát trực tuyến (Real-time Online Monitoring).

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn nội tại cần được tiếp tục hoàn thiện:

  1. Mô hình vết nứt đàn hồi tuyến tính đơn giản hóa: Vết nứt được giả định là dạng nứt hở một phía (open crack) và mô hình hóa bằng cặp lò xo tập trung tuyến tính. Trong thực tế vận hành tải trọng chu kỳ, vết nứt thường là dạng "nứt thở" (breathing crack), gây ra hiện tượng phi tuyến động học tiếp xúc mở - đóng tuần hoàn.
  2. Giới hạn trường biến dạng 1D/2D Timoshenko: Mô hình chưa xét đến hiệu ứng ứng suất 3 chiều tại lân cận đỉnh vết nứt (stress singularity) và hiệu ứng phân tách lớp (delamination) giữa bề mặt FGM và lớp áp điện dưới tải trọng rung động cường độ lớn.
  3. Bỏ qua ảnh hưởng của trường nhiệt độ và độ ẩm môi trường: Tính chất áp điện ($h_{13}, \beta_{33}$) và cơ lý tính FGM phụ thuộc đáng kể vào nhiệt độ môi trường, có thể tạo ra các trôi dạt tín hiệu (thermal drift) nếu không có thuật toán bù nhiệt.
  4. Quy mô thực nghiệm: Kết quả chẩn đoán vết nứt trong luận án chủ yếu dựa trên mô phỏng số và dữ liệu giải tích mở rộng độ chính xác cao; cần mở rộng thực nghiệm vật lý trên bàn rung quy mô công nghiệp.

Chương trình nghiên cứu 5–10 năm tiếp theo:

  • Mở rộng thuật toán DSM cho bài toán vết nứt thở phi tuyến (Nonlinear breathing crack DSM).
  • Khảo sát tương tác Đa trường (Thermo-Electro-Mechanical Coupled Fields) trên kết cấu vỏ và tấm FGM đa hướng (2D/3D-FGM).
  • Tích hợp mạng nơ-ron học sâu (Physics-Informed Neural Networks - PINNs) với cơ sở dữ liệu đường đồng mức điện tích áp điện để chẩn đoán hư hỏng tự động đa điểm nứt theo thời gian thực.

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

  • Học thuật: Công trình đóng góp 02 bài báo quốc tế uy tín thuộc danh mục ISI/Scopus, 01 bài báo trên Tạp chí Cơ học Việt Nam và nhiều báo cáo hội nghị toàn quốc, thiết lập tài liệu tham khảo nền tảng cho cộng đồng cơ học tính toán trong nước và quốc tế.
  • Chuyển đổi công nghiệp: Định hình giải pháp SHM thế hệ mới cho ngành công nghiệp hàng không vũ trụ, năng lượng tái tạo (điện gió ngoài khơi) và giao thông đường sắt cao tốc thông qua cảm biến áp điện màng mỏng tích hợp.
  • Chính sách & An toàn công trình: Cung cấp luận cứ khoa học để xây dựng tiêu chuẩn kỹ thuật về kiểm định an toàn và quan trắc sức khỏe tự động cho các công trình hạ tầng trọng điểm quốc gia.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận phương pháp độ cứng động lực (DSM) và mã nguồn thuật toán giải tích tường minh cho kết cấu thông minh.
  • Chuyên gia Cơ kỹ thuật & Kết cấu: Sở hữu công cụ tính toán dao động tần số cao chính xác tuyệt đối mà không gặp bẫy suy biến ma trận của phương pháp PTHH.
  • Kỹ sư R&D công nghiệp: Ứng dụng quy trình chẩn đoán vết nứt bằng đường đồng mức điện tích áp điện vào thiết kế các hệ thống giám sát sức khỏe kết cấu tự động (Smart SHM Systems).
  • Cơ quan quản lý an toàn hạ tầng: Thụ hưởng công nghệ giảm thiểu chi phí kiểm định định kỳ, chuyển dịch từ bảo trì thụ động sang bảo trì dự đoán (Predictive Maintenance).

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 ghép cặp điện - cơ cho dầm FGM Timoshenko có vết nứt gắn lớp áp điện phân bố, tính toán chính xác độ dịch chuyển trục trung hòa thực $h_0 = n(r_e - 1)h / [2(n+2)(n+r_e)]$. Công trình đã mở rộng lý thuyết cảm biến dạng mode của Lee & Moon (1990) từ vật liệu đồng nhất sang môi trường composite FGM biến thiên cơ tính phức tạp, đồng thời phát triển phương pháp độ cứng động lực (DSM) của Banerjee cho hệ cơ điện đa pha.

2. Đột phá phương pháp luận so với các nghiên cứu tiền tiêu biểu?

So với nghiên cứu $p$-FEM của Yu & Chu (2009) (sai số 15% khi có 1–2% nhiễu) và mô hình 3D-FEM kết hợp GA của Banerjee et al. (2014) (tốn kém chi phí tính toán), phương pháp DSM kết hợp đường đồng mức điện tích cảm biến dao động của luận án giải quyết bài toán thuần nhất bằng nghiệm giải tích dạng đóng. Phương pháp này giảm thiểu thời gian tính toán hàng trăm lần, không phụ thuộc vào kích thước chia lưới phần tử và triệt tiêu sai số làm tròn ở miền tần số cao.

3. Phát hiện bất ngờ nhất có dữ liệu thực nghiệm/mô phỏng chứng minh?

Phát hiện bất ngờ nhất là: Điện tích cảm biến dao động $Q$ không bị triệt tiêu độ nhạy tại các điểm nút tần số (Frequency Nodal Points). Trong khi các phương pháp chẩn đoán truyền thống hoàn toàn "mù" khi vết nứt xuất hiện tại vị trí nút dao động (nơi $\Delta\omega = 0$), điện tích cảm biến $Q$ vẫn ghi nhận bước nhảy tín hiệu đột biến do sự gián đoạn của đạo hàm góc xoay $\theta'(x)$, bảo đảm khả năng phát hiện vết nứt đạt độ chính xác 100%.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ?

Luận án cung cấp đầy đủ: (1) Hệ phương trình vi phân và ma trận cấu tử $[A], [B], [C]$; (2) Thuật toán thành lập ma trận độ cứng động $[D_e(\omega)]$; (3) Toàn bộ hằng số vật liệu đàn hồi, áp điện (PZT-4, Alumina/Steel) và sơ đồ khối thuật toán tính toán trong môi trường MATLAB, bảo đảm khả năng tái lập độc lập 100%.

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

Chiến lược 10 năm tập trung vào 3 trụ cột: (1) Chuyển đổi mô hình vết nứt từ 1D sang mô hình phá hủy 3D đa vết nứt tương tác; (2) Tích hợp công nghệ cảm biến áp điện không dây IoT thu nhận năng lượng tự thân (Energy Harvesting); (3) Xây dựng hệ thống bản sao số (Digital Twin) kết hợp AI cho các công trình siêu trường, siêu trọng.

Kết luận

  1. Xây dựng thành công hệ phương trình dao động ghép cặp điện - cơ chính xác cho dầm Timoshenko FGM có xét đến vị trí thực của trục trung hòa và hiệu ứng áp điện phân bố.
  2. Phát triển hoàn chỉnh ma trận độ cứng động lực (DSM) cho phân tố dầm FGM áp điện, giải quyết triệt để bài toán dao động tần số cao vượt trội so với phương pháp phần tử hữu hạn truyền thống.
  3. Thiết lập mô hình giải tích dầm FGM có vết nứt mô tả bằng cặp lò xo tương đương dọc trục và xoay, cho phép khảo sát toàn diện trường cơ điện phụ thuộc vào vị trí và độ sâu nứt.
  4. Phát minh phương pháp đường đồng mức điện tích cảm biến dao động (Modal Sensor Output Charge Contour Method), loại bỏ hoàn toàn hiện tượng đa nghiệm và điểm mù tại nút dao động.
  5. Kiểm chứng độ chính xác chẩn đoán vượt bậc với sai số vị trí vết nứt $< 1%$ trên các điều kiện biên đối xứng và không đối xứng (SS, CC, CF).
  6. Mở ra nhánh nghiên cứu liên ngành mới kết hợp Cơ học tính toán, Vật liệu thông minh và Công nghệ chẩn đoán kết cấu tự động (Smart SHM), nâng cao độ an toàn và tuổi thọ cho các công trình hạ tầng hiện đại.