Tổng quan luận án

Công tác đánh giá trạng thái kỹ thuật và chẩn đoán hư hỏng kết cấu (Structural Damage Detection) trong khuôn khổ kiểm tra sức khỏe công trình (Structural Health Monitoring - SHM) là nhiệm vụ then chốt nhằm bảo đảm an toàn vận hành và nâng cao tuổi thọ của các công trình công nghiệp, dân dụng, hàng không và năng lượng. Trong các kỹ thuật dò tìm khuyết tật, phương pháp thử nghiệm động lực học (Dynamic Testing Method) dựa trên việc tạo kích động, đo đạc đáp ứng và xử lý số liệu được áp dụng rộng rãi. Sự tích hợp vật liệu áp điện (piezoelectric material) vào kết cấu cho phép cơ điện hóa và tự động hóa quá trình kích hoạt cũng như đo đạc phản ứng dao động.

Trong khi các bài toán chẩn đoán hư hỏng cho kết cấu vật liệu đồng nhất (homogeneous) đã được khảo sát tương đối toàn diện, bài toán đối vớ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 Material - FGM) có tích hợp phần tử áp điện vẫn còn nhiều khoảng trống nghiên cứu. Các phương pháp chẩn đoán vết nứt truyền thống dựa trên tần số riêng thường gặp hạn chế do độ nhạy thấp đối với khuyết tật cục bộ, không thể nhận diện vết nứt khi vị trí hư hỏng trùng với điểm nút dao động (nodal points), hoặc gặp hiện tượng đa nghiệm (không tồn tại nghiệm duy nhất) khi kết cấu có điều kiện biên đối xứng. Bên cạnh đó, việc sử dụng các cảm biến áp điện dạng miếng dán cục bộ (patches) ở khoảng cách xa vết nứt thường không thu được tín hiệu rõ ràng, trong khi các phương pháp dựa trên trở kháng cơ điện (EMI) dễ bị che lấp bởi tín hiệu điện trường.

Luận án này được thực hiện nhằm giải quyết các khoảng trống trên thông qua việc phát triển mô hình tính toán và đề xuất phương pháp chẩn đoán vết nứt trong dầm FGM sử dụng lớp áp điện phân bố liên tục.

Mục tiêu nghiên cứu của luận án gồm:

  1. Xây dựng mô hình độ cứng động cho phân tố và phần tử dầm FGM Timoshenko có lớp áp điện, phục vụ phân tích động lực học dầm nguyên vẹn có gắn các miếng áp điện theo mô hình dầm bậc.
  2. Xây dựng mô hình dầm FGM có vết nứt được phủ lớp áp điện, khảo sát sự phụ thuộc của các đặc trưng cơ học (tần số riêng, dạng dao động) và đặc trưng điện (điện tích cảm biến dao động) vào các tham số vết nứt và thông số vật liệu.
  3. Đề xuất và phát triển phương pháp chẩn đoán vết nứt trong dầm FGM bằng kỹ thuật đường đồng mức của đặc trưng điện tích cảm biến dao động.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Kết cấu dầm Timoshenko FGM có cơ lý tính biến thiên liên tục theo chiều dày theo quy luật hàm lũy thừa, có vết nứt ngang dạng nứt mở được mô hình hóa bằng cặp lò xo đàn hồi tương đương (lò xo dọc trục và lò xo xoắn uốn), được gắn các miếng áp điện cục bộ hoặc một lớp áp điện phân bố liên tục dọc theo chiều dài.
  • Phạm vi nghiên cứu: Kết cấu dạng dầm thanh phẳng trong xây dựng công trình; nghiên cứu lý thuyết giải tích mở rộng kết hợp mô phỏng số thông qua phương pháp độ cứng động (DSM) được lập trình và tính toán trong môi trường phần mềm MATLAB.

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu trong luận án được phân tích qua các nhánh nghiên cứu chính sau:

  1. Nghiên cứu chẩn đoán hư hỏng trên kết cấu đồng nhất và kết cấu composite: Các nghiên cứu nền tảng về phát hiện khuyết tật trên kết cấu đồng nhất thông qua tham số động lực học đã được thực hiện bởi Cawley và Adams [1], Fan và Qiao [2], Hou và Xia [3], Pandey và Samman [4], Salawu [5], Khiem và Lien [6]. Đối với kết cấu vật liệu composite và vật liệu có cơ tính biến thiên (FGM), nghiên cứu tổng quan của Montalvão và cộng sự [7] đã chỉ ra tính phức tạp trong ứng xử động lực học khi có sự không đồng nhất về vật liệu.

  2. Nghiên cứu dao động của kết cấu FGM: Các đặc trưng dao động của dầm và tấm FGM được nghiên cứu bài bản bởi Hu và Zhang [72], Li [73], Larbi và cộng sự [74], Sina và cộng sự [75], Su và Banerjee [76], Wang và cộng sự [77]. Bài toán dầm FGM có vết nứt được khảo sát bởi Yang và Chen [78, 79], Akbas [80], Aydin [81], Khiem và cộng sự [82, 83]. Đáng chú ý, Khiem và cộng sự [83] đã phát triển phương pháp độ cứng động cho dầm FGM có nhiều vết nứt, tạo tiền đề lý thuyết trực tiếp cho luận án này. Về chẩn đoán vết nứt trong dầm FGM, Yu và Chu [84] sử dụng phương pháp phần tử hữu hạn p-version kết hợp đường đồng mức tần số; Banerjee và cộng sự [85] sử dụng mô hình PTHH 3 chiều kết hợp thuật toán di truyền (GA); Nguyễn Tiến Khiêm và Nguyễn Ngọc Huyên [86] thiết lập phương trình tần số giải tích tường minh cho dầm Timoshenko FGM nhưng chỉ ra hạn chế khi vết nứt nằm tại điểm nút dao động.

  3. Nghiên cứu kết cấu tích hợp vật liệu áp điện và ứng dụng chẩn đoán: Sau phát hiện về hiệu ứng áp điện của anh em nhà Curie (1880), Crawley và cộng sự [13, 14], Lee và Moon [15] đã xây dựng lý thuyết về cảm biến và bộ kích hoạt áp điện phân bố cho dầm và tấm. Sunar và Rao [12] tổng kết các hướng ứng dụng gồm: dập tắt dao động (Bailey và Hubbard [17], Bendine và cộng sự [18]), hấp thụ năng lượng (Wang và cộng sự [19]), điều khiển hình dáng (Batra và Liang [20]) và chẩn đoán hư hỏng. Việc ứng dụng áp điện để sửa chữa kết cấu nứt được khảo sát bởi Ariaei và cộng sự [21], Kumar và cộng sự [22], Liu [23]. Phân tích tĩnh và động dầm/tấm gắn miếng áp điện được nghiên cứu bởi Wang và cộng sự [32], Yang và Lee [33].

Trong lĩnh vực chẩn đoán hư hỏng bằng vật liệu áp điện, Quek và cộng sự [60] chứng minh cảm biến áp điện kết hợp biến đổi wavelet nhạy hơn lá điện trở khi vết nứt có độ sâu dưới 20%. Zhao và cộng sự [61, 62, 63] cùng Jiang và cộng sự [64] phát triển phương pháp điều chỉnh quét (Scan Tuning Method). Misra [65] khảo sát dầm công xôn có vết nứt thở bằng các cặp cảm biến áp điện. Phương pháp trở kháng cơ điện (EMI) được nghiên cứu bởi Giurgiutiu và Rogers [66], Ritdumrongkul và Fujino [67], Lim và Soh [68], Wang và Zhu [69, 70], Na và Baek [58], Tseng và Naidu [71], Park và cộng sự [53]. Đối với dầm FGM có lớp áp điện, Khorramabadi và Nezamabadi [87], Li và cộng sự [88], Bendine và cộng sự [18] đã phân tích bài toán ổn định và điều khiển, trong khi Nguyễn Tiến Khiêm và cộng sự khảo sát ảnh hưởng của miếng vá áp điện đến tần số riêng.

Tại Việt Nam, các nghiên cứu về kết cấu áp điện chủ yếu tập trung vào phân tích tĩnh và động của tấm composite/FGM nhiều lớp chịu nhiệt, nền đàn hồi (nhóm GS. Trần Ích Thịnh, GS. Trần Minh Tú tại ĐHBKHN, ĐHXDHN [44, 45, 46, 47]; các nhóm tại ĐHQGHN [48, 49]; ĐHKT TP.HCM [50, 51]), chưa tập trung vào cơ chế chẩn đoán hư hỏng bằng tín hiệu điện tích phân bố.

Khoảng trống nghiên cứu được luận án lựa chọn giải quyết:

  • Thiếu mô hình độ cứng động lực học chính xác ở dải tần số cao cho dầm FGM Timoshenko có xét đến vị trí thực của trục trung hòa và tích hợp lớp/miếng áp điện.
  • Hạn chế của các phương pháp chẩn đoán hư hỏng hiện có: cảm biến cục bộ bị giảm độ nhạy khi ở xa khuyết tật; phương pháp tần số riêng truyền thống thất bại tại điểm nút dao động và đối xứng biên.
  • Chưa có công trình công bố về việc xây dựng cơ sở dữ liệu và thuật toán chẩn đoán vết nứt trong dầm FGM sử dụng đường đồng mức điện tích cảm biến dao động từ lớp áp điện phủ liên tục.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án kết hợp lý thuyết giải tích mở rộng và kỹ thuật mô phỏng tính toán số:

  1. Lý thuyết dầm Timoshenko FGM và hiệu ứng áp điện:
  • Quy luật biến thiên vật liệu: Các đặc tính vật liệu gồm mô đun đàn hồi $E(z)$, mô đun trượt $G(z)$ và khối lượng riêng $\rho(z)$ biến thiên theo chiều cao $z$ theo hàm lũy thừa: $$\Re(z) = \Re_b + (\Re_t - \Re_b)\left(\frac{z}{h} + 0.5\right)^n$$ trong đó chỉ số $t$ và $b$ chỉ lớp vật liệu mặt trên và mặt dưới dầm; $n$ là chỉ số gradient (chỉ số phân bố vật liệu).
  • Vị trí trục trung hòa thực tế: Do tính không đối xứng của vật liệu FGM, khoảng cách $h_0$ từ mặt phẳng giữa dầm đến trục trung hòa được xác định theo công thức: $$h_0 = \frac{n(r_e - 1)h}{2(n + 2)(n + r_e)}, \quad \text{với } r_e = \frac{E_t}{E_b}$$
  • Phương trình cấu tạo của lớp áp điện (polyme PVDF hoặc gốm PZT) xét theo lý thuyết áp điện tuyến tính 1D: $$\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$ là mô đun đàn hồi, $h_{13}$ là hằng số áp điện, $\beta_{33}^p$ là hằng số điện môi, $\in$ là cường độ điện trường và $D$ là mật độ dịch chuyển điện tích.
  • Điều kiện liên kết: Giả thiết tiếp xúc hoàn toàn (không trượt và không bong tách) giữa dầm chủ và lớp áp điện.
  1. Phương pháp độ cứng động (Dynamic Stiffness Method - DSM):
  • Áp dụng nguyên lý biến phân Hamilton $\int_{t_1}^{t_2} \delta (T - \Pi) dt = 0$ với tổng động năng $T = T_b + T_p$ và thế năng biến dạng $\Pi = \Pi_b + \Pi_p$ để thiết lập hệ phương trình vi phân dao động cơ - điện ghép cặp trong miền thời gian, sau đó chuyển sang miền tần số qua phép biến đổi Fourier: $$[\mathbf{A}]{\mathbf{Z}''(x,\omega)} + [\mathbf{B}]{\mathbf{Z}'(x,\omega)} + [\mathbf{C}]{\mathbf{Z}(x,\omega)} = 0$$
  • Giải phương trình vi phân bằng cách tìm nghiệm dạng $\mathbf{Z}_0 = \mathbf{d} e^{\lambda x}$, dẫn đến phương trình đặc trưng đại số bậc ba đối với $\eta = \lambda^2$.
  • Thiết lập quan hệ giữa vectơ lực nút ${\mathbf{P}_e(\omega)} = {N_1, M_1, Q_1, N_2, M_2, Q_2}^T$ và vectơ chuyển vị nút ${\mathbf{U}_e(\omega)} = {U_1, \Theta_1, W_1, U_2, \Theta_2, W_2}^T$ thông qua việc khử vectơ hằng số tích phân ${\mathbf{C}}$: $${\mathbf{U}_e(\omega)} = [\mathbb{Z}(\omega)]{\mathbf{C}}, \quad {\mathbf{P}_e(\omega)} = [\mathbb{Q}(\omega)]{\mathbf{C}} \implies [\mathbf{D}_e(\omega)] = [\mathbb{Q}(\omega)][\mathbb{Z}(\omega)]^{-1}$$
  • Ma trận độ cứng động toàn kết cấu được ghép nối theo quy tắc phần tử: $$[\mathbf{D}(\omega)] = \sum_{e=1}^{n_e} [\mathbf{T}_e]^T [\mathbf{D}_e(\omega)] [\mathbf{T}_e]$$ Phương trình tần số dao động riêng có dạng: $\det[\mathbf{D}(\omega)] = 0$.
  1. Mô hình vết nứt và đại lượng điện tích cảm biến dao động:
  • Vết nứt mở tại tọa độ $x = e$ với độ sâu $a$ được mô hình hóa thành một phân tố liên kết gồm hai lò xo đàn hồi tương đương: lò xo xoay uốn (mô tả độ giảm độ cứng chống uốn) và lò xo dọc trục (mô tả độ giảm độ cứng chịu kéo/nén).
  • Đại lượng điện tích cảm biến dao động (Modal Sensor Charge - ký hiệu $Q$ hoặc MPC) sinh ra trên bề mặt lớp áp điện khi dầm dao động ở dạng riêng thứ $j$ được xác định theo biểu thức: $$Q = \int D dA = b \int_0^L D dx = \frac{b h_{13}}{\beta_{33}^p} \int_0^L \left( u_0'(x) - \frac{h}{2} \theta'(x) \right) dx$$

Nội dung chính theo từng chương

Chương 1: Tổng quan và đặt vấn đề nghiên cứu

Chương 1 hệ thống hóa cơ sở lý thuyết về vật liệu áp điện (hiệu ứng thuận và ngược), phân loại vật liệu gốm PZT và polyme PVDF, cùng các ứng dụng trong kỹ thuật siêu âm, y tế, động cơ vi mô và thu hồi năng lượng. Tiếp đó, chương điểm luận chi tiết các công trình nghiên cứu trong và ngoài nước về dao động kết cấu FGM, kết cấu áp điện và các phương pháp chẩn đoán hư hỏng (phương pháp trở kháng cơ điện EMI, phương pháp điều chỉnh quét, phương pháp đường đồng mức tần số). Trên cơ sở phân tích các ưu điểm và hạn chế của phương pháp PTHH (bị giới hạn dải tần) và các phương pháp chẩn đoán dùng cảm biến rời rạc, chương xác lập tính cấp thiết, mục tiêu và 3 nội dung nghiên cứu chính của luận án.

Chương 2: Dao động của dầm FGM nguyên vẹn có miếng áp điện

Chương này tập trung thiết lập mô hình toán học và thuật toán số phân tích dao động tự do của dầm FGM nguyên vẹn có gắn các miếng áp điện:

  • Xây dựng hệ phương trình vi phân chuyển động dầm Timoshenko FGM có lớp áp điện, tính toán đầy đủ ảnh hưởng của độ lệch trục trung hòa $h_0$ và tính chất đàn hồi - điện ghép cặp ($E_p = C_{11}^p - h_{13}^2 / \beta_{33}^p$).
  • Xây dựng ma trận độ cứng động $[\mathbf{D}_e(\omega)]$ kích thước $6 \times 6$ cho phần tử dầm FGM áp điện và mô hình hóa dầm có miếng áp điện như một kết cấu dầm bậc đa phân tử.
  • Khảo sát số cho tham số tần số không thứ nguyên $\lambda = \omega (L^2/h)\sqrt{\rho_b/E_b}$ với các thông số vật liệu: $E_t = 390\text{ MPa}$, $\rho_t = 3960\text{ kg/m}^3$, $\mu_t = 0.3$.
  • So sánh kết quả tần số với công bố của Su và Banerjee [76] (khi $h_p = 0$) cho thấy sự phù hợp, xác nhận độ tin cậy của thuật toán DSM.
  • Đánh giá sự phụ thuộc của 5 mode tần số đầu vào tỷ số độ mảnh $L_b/h_b$, tỷ số chiều dày miếng dán $h_p/h_b$, tỷ số chiều dài $L_p/L_b$, vị trí gắn (đầu dầm, giữa dầm, suốt chiều dài) và chỉ số gradient vật liệu $n$ ($n = 0; 0.5; 1; 2; 5$).
Điều kiện biên dầm Vị trí gắn miếng áp điện Quy luật biến thiên của tần số chuẩn hóa theo độ dày $h_p/h_b$
Tựa đơn (SS) Tại các đầu dầm Biến thiên dạng phi tuyến (tương tự dạng parabol)
Tựa đơn (SS) Tại giữa dầm Giảm hoặc biến thiên không đơn điệu tùy thuộc mode dao động
Ngàm - Ngàm (CC) Tại đầu ngàm Tần số chuẩn hóa mode 1 tăng đơn điệu theo chiều dày $h_p$
Ngàm - Ngàm (CC) Tại giữa dầm Cả hai tần số chuẩn hóa đầu tiên giảm đơn điệu khi $h_p$ tăng
Ngàm - Tự do (CF) Tại đầu ngàm Tần số chuẩn hóa mode 1 tăng đơn điệu theo chiều dày $h_p$
Ngàm - Tự do (CF) Tại đầu tự do Cả hai tần số chuẩn hóa đầu tiên giảm đơn điệu khi $h_p$ tăng
Tất cả (SS, CC, CF) Toàn bộ chiều dài Tần số chuẩn hóa biến thiên dạng parabol đối với mọi điều kiện biên

Chương 3: Dao động của dầm FGM bị nứt có lớp áp điện

Chương 3 mở rộng bài toán sang dầm FGM có vết nứt và được phủ một lớp áp điện liên tục:

  • Mô hình hóa vết nứt mở ngang thông qua điều kiện gián đoạn chuyển vị và ma trận truyền qua phân tố lò xo tương đương.
  • Thiết lập quy trình tính toán nghiệm tổng quát và phát triển chương trình tính toán số trên MATLAB để xác định các đặc trưng động lực học cơ - điện.
  • Khảo sát sự thay đổi của các tần số riêng và các dạng dao động (bao gồm dạng dao động dọc trục A1 và các dạng dao động uốn B1-B4) khi thay đổi vị trí vết nứt $e/L$ và độ sâu vết nứt tương đối $a/h$.
  • Khảo sát chi tiết đại lượng điện tích cảm biến dao động $Q$ (MPC) cho 3 dạng riêng đầu tiên. Kết quả số chỉ ra:
    1. Tần số riêng giảm khi độ sâu vết nứt $a$ tăng; mức độ giảm phụ thuộc mạnh vào vị trí vết nứt $e$ và đạt cực đại tại các vị trí có mô men uốn lớn nhất.
    2. Điện tích cảm biến dao động $Q$ phản ánh độ biến thiên độ cong và góc xoay của dầm, thể hiện độ nhạy đối với sự xuất hiện của vết nứt, đồng thời mang dấu hiệu phân biệt rõ ràng giữa các vị trí khuyết tật khác nhau.
    3. Ảnh hưởng của chỉ số gradient vật liệu $n$ và chiều dày lớp áp điện $h_p$ đến phản ứng điện tích dao động được lượng hóa chi tiết, cung cấp dữ liệu đầu vào cho bài toán nhận dạng ngược.

Chương 4: Chẩn đoán vết nứt trong dầm FGM có lớp áp điện

Chương 4 trình bày thuật toán và kết quả giải bài toán ngược nhằm xác định vị trí $e$ và độ sâu $a$ của vết nứt:

  • Xây dựng cơ sở dữ liệu chẩn đoán vết nứt theo hai phương pháp:
    1. Cơ sở dữ liệu đường đồng mức tần số riêng: Dựa trên tập hợp các đường đẳng tần số trong không gian tham số $(e/L, a/h)$ tương ứng với từng mode dao động.
    2. Cơ sở dữ liệu đường đồng mức điện tích cảm biến dao động: Dựa trên tập hợp các đường đẳng điện tích cảm biến dao động $Q_j(e, a) = \text{const}$.
  • Nguyên lý chẩn đoán: Giao điểm duy nhất của các đường đồng mức từ ít nhất hai mode dao động khác nhau xác định chính xác tọa độ vị trí và độ sâu vết nứt.
  • Thử nghiệm số trên các dầm tựa đơn (SS), dầm ngàm hai đầu (CC) và dầm công xôn (CF) với hai trường hợp vết nứt:
    • Vết nứt tại giữa dầm ($e = L/2$) có độ sâu nứt $a/h = 30%$.
    • Vết nứt tại vị trí lệch tâm ($e = L/3$) có độ sâu nứt $a/h = 30%$.
  • So sánh hiệu quả giữa hai phương pháp: Phương pháp đường đồng mức tần số gặp hiện tượng đa giao điểm hoặc sai số lớn tại các điểm nút dao động và trường hợp dầm có biên đối xứng; trong khi đó, phương pháp đường đồng mức điện tích cảm biến dao động loại bỏ được tính đa nghiệm, xác định hội tụ chính xác vị trí và kích thước vết nứt.
Điều kiện biên dầm Trường hợp vết nứt khảo sát Kết quả chẩn đoán bằng đường đồng mức tần số Kết quả chẩn đoán bằng đường đồng mức điện tích áp điện
Tựa đơn (SS) $e = L/2, a/h = 30%$ Bị ảnh hưởng bởi tính đối xứng của dạng dao động Xác định chính xác tọa độ giao điểm $(e, a)$
Tựa đơn (SS) $e = L/3, a/h = 30%$ Nhận diện được vị trí nhưng có thể xuất hiện nghiệm ảo đối xứng Cho nghiệm hội tụ duy nhất tại vị trí $e = L/3$
Ngàm - Ngàm (CC) $e = L/2, a/h = 30%$ Trùng điểm nút của mode chẵn, gây khó khăn cho việc giao cắt Tín hiệu điện tích cho giao điểm rõ ràng
Ngàm - Ngàm (CC) $e = L/3, a/h = 30%$ Độ dốc đường đồng mức thấp tại một số dải Xác định chính xác vị trí và độ sâu 30%
Công xôn (CF) $e = L/2, a/h = 30%$ Cho giao điểm tương đối tốt Độ chính xác cao, loại bỏ nhiễu nghiệm
Công xôn (CF) $e = L/3, a/h = 30%$ Cho kết quả phù hợp Xác định chính xác tham số khuyết tật

Kết quả và những đóng góp mới

Các đóng góp mới về mặt khoa học và kỹ thuật được tác giả tự xác định trong luận án gồm:

  1. Đã xây dựng thành công mô hình độ cứng động (DSM) cho dầm FGM Timoshenko có lớp áp điện, có xét đến vị trí thực của trục trung hòa và hiệu ứng cơ điện liên hợp; áp dụng mô hình này để tính toán dao động tự do của dầm FGM có gắn các miếng áp điện cục bộ theo mô hình kết cấu dầm bậc.
  2. Đã xây dựng mô hình dầm FGM có vết nứt được phủ một lớp áp điện liên tục; thiết lập chương trình tính toán số phân tích toàn diện sự phụ thuộc của các đặc trưng động lực học (tần số riêng, dạng dao động uốn - dọc trục) và đặc trưng điện (điện tích cảm biến dao động) vào vị trí vết nứt, độ sâu vết nứt và chỉ số phân bố vật liệu.
  3. Đã đề xuất và phát triển phương pháp đường đồng mức áp điện (sử dụng đặc trưng điện tích cảm biến dao động sinh ra trên lớp áp điện phân bố) phục vụ chẩn đoán vết nứt trong dầm FGM; chứng minh phương pháp này giải quyết được vấn đề điểm nút dao động và tính đối xứng biên vốn là nhược điểm của phương pháp tần số riêng truyền thống.

Danh mục các công trình công bố của tác giả: Các kết quả nghiên cứu chính của luận án đã được công bố trong các công trình khoa học chuyên ngành, bao gồm:

  • 02 bài báo trên các tạp chí quốc tế thuộc danh mục ISI.
  • 01 bài báo trên Tạp chí Cơ học Việt Nam.
  • Các bài báo đăng toàn văn trong Tuyển tập Hội nghị Cơ học toàn quốc và Tuyển tập Hội nghị Cơ học vật rắn.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Về mô hình hình học và kết cấu: Luận án giới hạn phạm vi mô hình hóa ở kết cấu dầm thanh phẳng một chiều (1D) theo lý thuyết dầm Timoshenko; chưa mở rộng cho kết cấu tấm, vỏ FGM phức tạp hoặc hệ khung không gian 3D.
  • Về mô hình vật liệu: Quy luật biến thiên cơ tính FGM chỉ xét theo phương chiều dày dầm theo hàm số lũy thừa 1 hướng; chưa xét biến thiên đồng thời nhiều hướng hoặc ảnh hưởng của trường nhiệt độ và độ ẩm môi trường.
  • Về mô hình hư hỏng: Vết nứt được giả thiết là vết nứt mở ngang tuyến tính thông qua mô hình lò xo tương đương, chưa xét đến cơ chế phi tuyến của vết nứt đóng - mở (vết nứt thở) trong quá trình dao động.
  • Về khía cạnh thực nghiệm: Nghiên cứu dừng lại ở mức độ mô phỏng số giải tích trên MATLAB; việc chế tạo lớp cảm biến áp điện liên tục phủ kín chiều dài kết cấu trong thực tế đòi hỏi kỹ thuật gia công và xử lý tín hiệu thực tế phức tạp.

Hướng nghiên cứu phát triển tiếp:

  • Mở rộng phương pháp độ cứng động và kỹ thuật đường đồng mức điện tích áp điện cho kết cấu tấm FGM, vỏ FGM và dầm FGM có tiết diện thay đổi hoặc cơ tính biến thiên 2 hướng (2D-FGM).
  • Khảo sát mô hình vết nứt thở phi tuyến và bài toán chẩn đoán đa vết nứt đồng thời trên dầm áp điện.
  • Tiến hành các thử nghiệm thực nghiệm trên mô hình dầm thực tế để kiểm chứng và hiệu chỉnh độ chính xác của thuật toán chẩn đoán trước các tác động của nhiễu đo đạc.

Giá trị tham khảo

Luận án là tài liệu tham khảo hữu ích cho các đối tượng chuyên môn sau:

  • Nghiên cứu sinh, học viên cao học và nhà nghiên cứu thuộc các chuyên ngành Cơ kỹ thuật, Cơ học vật rắn biến dạng, Kỹ thuật xây dựng, Kỹ thuật cơ khí và Kỹ thuật công trình thủy.
  • Giảng viên đại học sử dụng làm tài liệu tham khảo chuyên đề cho các học phần Động lực học công trình nâng cao, Vật liệu và kết cấu thông minh, Phương pháp tính trong cơ học kết cấu.
  • Kỹ sư và chuyên gia trong lĩnh vực giám sát và chẩn đoán sức khỏe công trình (SHM), kiểm tra không phá hủy (NDT).

Các nội dung có giá trị tham khảo chuyên sâu gồm:

  • Quy trình toán học thiết lập ma trận độ cứng động $[\mathbf{D}_e(\omega)]$ cho phần tử dầm FGM Timoshenko có lớp áp điện và việc xác định chính xác trục trung hòa dầm FGM trong Chương 2.
  • Thuật toán giải bài toán dao động dầm FGM nứt tích hợp phân tố lò xo đàn hồi tương đương trong Chương 3.
  • Cơ sở dữ liệu và kỹ thuật chẩn đoán vết nứt bằng đường đồng mức điện tích cảm biến dao động trong Chương 4.

Câu hỏi thường gặp

  1. Trục trung hòa thực tế của dầm FGM được xác định như thế nào trong luận án? Trong dầm FGM, do mô đun đàn hồi biến thiên theo chiều cao dầm, trục trung hòa thực tế không trùng với trục hình học giữa dầm. Khoảng cách $h_0$ từ mặt phẳng giữa đến trục trung hòa được tính bằng công thức: $$h_0 = \frac{n(r_e - 1)h}{2(n + 2)(n + r_e)}$$ với $r_e = E_t / E_b$ (tỷ số mô đun đàn hồi giữa mặt trên và mặt dưới) và $n$ là chỉ số gradient phân bố vật liệu.

  2. Đại lượng điện tích cảm biến dao động $Q$ được định nghĩa và tính toán ra sao? Điện tích cảm biến dao động $Q$ (hay Modal Sensor Charge) là tổng điện tích sinh ra trên bề mặt lớp áp điện khi dầm dao động ở một dạng riêng xác định. Đại lượng này được tính bằng cách tích phân mật độ dịch chuyển điện tích $D$ trên toàn bộ diện tích bề mặt lớp áp điện: $$Q = \frac{b h_{13}}{\beta_{33}^p} \int_0^L \left( u_0'(x) - \frac{h}{2} \theta'(x) \right) dx$$ trong đó $b$ là chiều rộng dầm, $h_{13}$ là hằng số áp điện, $\beta_{33}^p$ là hằng số điện môi, $u_0(x)$ là chuyển vị dọc trục và $\theta(x)$ là góc xoay mặt cắt ngang.

  3. Phương pháp độ cứng động (DSM) có ưu điểm gì nổi bật so với phương pháp phần tử hữu hạn (PTHH) trong phân tích kết cấu áp điện? Phương pháp DSM sử dụng nghiệm giải tích chính xác của phương trình vi phân chuyển động làm hàm dạng, do đó ma trận độ cứng động $[\mathbf{D}(\omega)]$ chứa cả độ cứng và khối lượng phân bố liên tục. Điều này cho phép DSM tính toán chính xác dao động ở mọi dải tần số cao mà không cần chia nhỏ lưới phần tử như PTHH, rất phù hợp với đặc tính hoạt động tần số cao của vật liệu áp điện.

  4. Tại sao phương pháp đường đồng mức điện tích áp điện lại khắc phục được nhược điểm của phương pháp đường đồng mức tần số? Phương pháp đường đồng mức tần số truyền thống gặp khó khăn khi vết nứt nằm tại các điểm nút dao động (nơi vết nứt không làm thay đổi tần số tương ứng) hoặc trên kết cấu đối xứng (dẫn đến đa nghiệm vị trí). Ngược lại, điện tích cảm biến dao động $Q$ tích hợp toàn bộ biến dạng dọc trục và góc xoay trên suốt chiều dài dầm, do đó vẫn nhạy cảm tại điểm nút tần số và cung cấp các đường đồng mức độc lập, giúp tạo giao điểm duy nhất $(e, a)$.

  5. Những điều kiện biên và vị trí vết nứt nào đã được thử nghiệm số để kiểm chứng phương pháp chẩn đoán trong Chương 4? Luận án thực hiện thử nghiệm số trên ba điều kiện biên cơ bản gồm: dầm tựa đơn hai đầu (SS), dầm ngàm hai đầu (CC) và dầm công xôn ngàm - tự do (CF). Với mỗi điều kiện biên, hai vị trí vết nứt được khảo sát là vết nứt chính giữa dầm ($e = L/2$) và vết nứt lệch tâm ($e = L/3$), đều có độ sâu vết nứt tương đối là $a/h = 30%$.

Kết luận

Luận án tiến sĩ của tác giả Lưu Quỳnh Hường đã phát triển thành công phương pháp độ cứng động để mô hình hóa và phân tích ứng xử dao động cơ - điện của dầm FGM Timoshenko có vết nứt được tích hợp vật liệu áp điện. Công trình đã làm rõ quy luật biến thiên của tần số riêng, các dạng dao động uốn - dọc trục và đặc trưng điện tích cảm biến dao động theo các thông số vết nứt, kích thước hình học và chỉ số phân bố vật liệu. Đóng góp cốt lõi của luận án là đề xuất phương pháp đường đồng mức điện tích áp điện, mang lại giải pháp tin cậy để xác định chính xác vị trí và độ sâu vết nứt, khắc phục triệt để các hạn chế cố hữu của phương pháp chẩn đoán bằng tần số riêng truyền thống. Các kết quả đạt được tạo lập cơ sở lý thuyết và công cụ tính toán số hữu ích cho việc nghiên cứu, thiết kế và phát triển các hệ thống giám sát sức khỏe kết cấu thông minh làm từ vật liệu composite biến tính.