Tổng quan về luận án

Nghiên cứu dao động và đánh giá tình trạng kết cấu công trình là một trong những trọng tâm cốt lõi của cơ học công trình hiện đại và kỹ thuật cơ điện tử tiên tiến. Luận án tiến sĩ với đề tài "Phân tích dao động và chẩn đoán kết cấu dầm bằng vật liệu có cơ tính biến thiên có nhiều vết nứt" đã thiết lập một bước tiến đột phá trong việc giải quyết đồng thời bài toán thuận (phân tích động lực học) và bài toán ngược (chẩn đoán tham số khuyết tật) cho kết cấu dầm chịu lực phức tạp.

Bối cảnh khoa học của đề tài bắt nguồn từ sự phát triển vượt bậc của vật liệu cơ tính biến thiên (Functionally Graded Material – FGM). Như nguyên văn công trình đã khẳng định: "Vật liệu cơ tính biến thiên (Functionally Graded Material – FGM) là một loại vật liệu composite tiên tiến, cấu thành từ hai pha vật liệu, có các đặc trưng cơ học biến đổi trơn, liên tục và các tính năng ưu việt của các thành phần vật liệu được phát huy tối đa, đồng thời tránh được sự bong tách, tập trung ứng suất tại các bề mặt tiếp xúc như thường xảy ra đối với các vật liệu composite truyền thống." Nhờ khả năng làm việc trong môi trường gradient nhiệt độ khắc nghiệt (chênh lệch nhiệt độ bề mặt có thể lên tới 1600K và nhiệt độ đỉnh đạt 2100K), FGM được ứng dụng rộng rãi trong hàng không vũ trụ, kỹ thuật hạt nhân, chế tạo máy và cơ học y sinh.

Tuy nhiên, trong quá trình vận hành, các kết cấu FGM không thể tránh khỏi sự xuất hiện của các khuyết tật cục bộ, đặc biệt là các hệ vết nứt ngang. Research gap cụ thể trong y văn quốc tế nằm ở chỗ: các nghiên cứu trước đây (như Yang & Chen 2008, Ke et al. 2009) chủ yếu giới hạn ở dầm Euler-Bernoulli đơn giản với một hoặc hai vết nứt bằng phương pháp giải tích cổ điển đòi hỏi giải định thức ma trận siêu việt bậc cao $6n \times 6n$, hoặc dựa vào phương pháp phần tử hữu hạn (FEM) truyền thống vốn đòi hỏi chia lưới dày đặc (trên 40 lớp vật liệu) và bộc lộ sai số lớn ở miền tần số bậc cao. Chưa có công trình nào giải quyết toàn diện bài toán dao động tự do, dao động cưỡng bức và chẩn đoán đồng thời số lượng, vị trí, độ sâu của nhiều vết nứt trên dầm Timoshenko FGM liên tục bằng phương pháp độ cứng động lực.

Để giải quyết khoảng trống này, luận án xác định các câu hỏi và giả thuyết nghiên cứu cụ thể:

  1. Research Question 1 (RQ1): Làm thế nào để mô hình hóa giải tích dầm Timoshenko FGM có vô số vết nứt ngang thành một phần tử động lực duy nhất mà vẫn đảm bảo tính đến ảnh hưởng đồng thời của biến dạng trượt, quán tính quay và sự dịch chuyển của trục trung hòa thực?
  2. Research Question 2 (RQ2): Mức độ ảnh hưởng tương hỗ giữa chỉ số tỷ lệ thể tích $n$, tỷ số mô-đun đàn hồi $E_t/E_b$, độ mảnh $L/h$ và các tham số vết nứt đến phổ tần số riêng, dạng dao động và đáp ứng động lực cưỡng bức diễn ra theo quy luật nào?
  3. Research Question 3 (RQ3): Làm thế nào để xây dựng thuật toán giải bài toán ngược có khả năng nhận dạng chính xác vị trí, số lượng và độ sâu của nhiều vết nứt đồng thời từ dữ liệu đo đạc có nhiễu?

Các giả thuyết nghiên cứu bao gồm:

  • Hypothesis 1 (H1): Việc áp dụng hàm dạng động chính xác thỏa mãn phương trình vi phân cân bằng động trong miền tần số sẽ loại bỏ hoàn toàn sai số xấp xỉ tĩnh của phương pháp phần tử hữu hạn cổ điển, cho phép xác định chính xác các đặc trưng động lực học bậc cao của dầm FGM nhiều vết nứt chỉ với một phần tử duy nhất giữa hai gối tựa.
  • Hypothesis 2 (H2): Vết nứt cục bộ gây ra sự gián đoạn vi mô trong trường chuyển vị, sự kết hợp giữa biến đổi wavelet dừng (SWT) đa tỉ lệ và mạng nơ-ron nhân tạo lan truyền ngược (BP-ANN) sẽ giải quyết triệt để bài toán ngược với độ chính xác cao và khả năng kháng nhiễu vượt trội.

Khung lý thuyết của luận án được xây dựng vững chắc dựa trên Lý thuyết dầm Timoshenko bậc nhất, Cơ học phá hủy đàn hồi tuyến tính (LEFM), Nguyên lý biến phân Hamilton, Phương pháp độ cứng động lực (Dynamic Stiffness Method – DSM) và Lý thuyết trí tuệ nhân tạo tính toán. Về quy mô học thuật, công trình tổng hợp 169 tài liệu tham khảo (16 tài liệu trong nước, 153 tài liệu quốc tế uy tín) và được bảo chứng chất lượng khoa học thông qua 14 công trình công bố, bao gồm 10 bài báo chuyên ngành trong nước và 04 bài báo trên các tạp chí quốc tế thuộc danh mục ISI/Scopus.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp giám sát trạng thái kỹ thuật công trình (Structural Health Monitoring – SHM) dựa trên đặc trưng động lực học đã trải qua nhiều giai đoạn với các trường phái tiếp cận khác nhau:

Trường phái truyền thống khởi đầu từ các nghiên cứu kinh điển của Adams, Cawley et al. (1979) trong việc sử dụng tần số tự nhiên và dạng dao động để định vị khuyết tật trên dầm đồng nhất. Về mặt mô hình hóa cơ học phá hủy, Thomson (1949) và Irwin (1957) đã đặt nền móng cho việc tính toán hệ số tập trung ứng suất (Stress Intensity Factors – SIF: $K_I, K_{II}, K_{III}$) tại đỉnh vết nứt. Đối với kết cấu dầm có vết nứt, các mô hình vết nứt tương đương được phát triển theo ba hướng: mô hình suy giảm độ cứng cục bộ (Christides & Barr 1984), mô hình ma trận độ mềm bổ sung (Dimarogonas 1996) và mô hình lò xo đàn hồi tập trung (Ostachowicz & Krawczuk 1991, Rizos et al. 1990).

Khi vật liệu FGM xuất hiện (Koizumi 1997), các nghiên cứu động lực học kết cấu bắt đầu chuyển hướng sang cấu trúc vật liệu dị tính. Yang và Chen (2008) đã khảo sát dao động tự do của dầm Euler-Bernoulli P-FGM có một vết nứt bằng phương pháp chia đoạn liên tục. Ke et al. (2009) mở rộng cho dầm Timoshenko FGM xét đến biến dạng trượt. Tuy nhiên, tranh luận học thuật lớn nảy sinh giữa hai trường phái:

  1. Trường phái Rời rạc hóa phần tử hữu hạn (FEM/p-FEM): Đại diện bởi Nazari et al. (2011), Yu & Chu (2009). Nazari chia dầm FGM thành hơn 40 lớp đồng nhất dọc theo chiều dày để tính toán trường biến dạng. Ưu điểm là trực quan nhưng nhược điểm chí mạng là số bậc tự do bùng nổ, khối lượng tính toán quá lớn và xuất hiện hiện tượng "khóa trượt" (shear locking) ở miền tần số cao.
  2. Trường phái Hàm sóng chính xác miền tần số (DSM/SFEM): Khởi xướng từ Williams & Wittrick (1970), Kolousek (1973), Leung (1993) và được Banerjee (2007, 2016) phát triển cho vật liệu composite và FGM nguyên vẹn. Phương pháp độ cứng động lực (ĐCĐL) sử dụng chính nghiệm giải tích của phương trình vi phân chuyển động làm hàm dạng động, đem lại nghiệm chính xác tuyệt đối mà không phụ thuộc vào kích thước phần tử chia lưới.

Về phương diện chẩn đoán vết nứt (bài toán ngược), các chỉ số cổ điển như Modal Assurance Criterion (MAC - Allemang & Brown 1982), Coordinate MAC (COMAC), Curvature Damage Factor (CDF - Pandey et al. 1991) hay Damage Index Method (DIM - Stubbs et al. 1995) chỉ định vị tốt khi có ít vết nứt trên dầm đồng nhất và suy giảm hiệu quả nghiêm trọng khi áp dụng cho kết cấu FGM chịu ảnh hưởng của biến đổi cơ tính phi tuyến. Phân tích Continuous Wavelet Transform (CWT) hay Discrete Wavelet Transform (DWT) của Rucka et al. (2006), Zhong & Oyadiji (2007) bộc lộ hiện tượng méo tín hiệu ở biên (boundary distortion).

Vị trí của luận án trong dòng chảy học thuật được xác lập rõ ràng: Luận án là công trình tiên phong kết hợp lý thuyết dầm Timoshenko P-FGM xét trục trung hòa thực với phương pháp độ cứng động lực để mô hình hóa hệ dầm nhiều vết nứt thành một phần tử duy nhất, sau đó tích hợp biến đổi wavelet dừng (Stationary Wavelet Transform – SWT) triệt tiêu hiệu ứng biên và mạng trí tuệ nhân tạo (ANN) đa tầng để giải quyết trọn vẹn bài toán SHM cấp độ 1, 2 và 3.

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

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

Luận án đã tạo ra bước đột phá lý thuyết quan trọng thông qua việc mở rộng và chuẩn hóa hệ phương trình vi phân dao động của kết cấu FGM dị tính:

Thứ nhất, tái định vị mặt trung hòa thực: Trong dầm đồng nhất, mặt trung hòa trùng với trục đối xứng hình học ($z = 0$). Tuy nhiên, đối với dầm P-FGM có mô-đun đàn hồi biến thiên liên tục theo quy luật hàm lũy thừa $E(z) = E_b + (E_t - E_b)(z/h + 1/2)^n$, ứng suất pháp phân bố bất đối xứng. Luận án đã xác định chính xác tọa độ dịch chuyển của trục trung hòa $h_0$ từ điều kiện triệt tiêu thành phần liên kết uốn-kéo tĩnh ($A_{12} = 0$): $$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}$$

Thứ hai, thiết lập điều kiện triệt tiêu ghép nối động lực uốn - dọc trục: Luận án chứng minh một định lý cơ học quan trọng: ma trận quán tính liên kết uốn-dọc $I_{12}$ được biểu diễn dưới dạng: $$I_{12} = \frac{n b h^2 \rho_b (R_\rho - R_E)}{2(2 + n)(n + R_E)}, \quad \text{với } R_\rho = \frac{\rho_t}{\rho_b}$$ Hiện tượng ghép nối động lực giữa dao động uốn và dao động dọc trục sẽ triệt tiêu hoàn toàn ($I_{12} = 0$) khi và chỉ khi $n = 0$ (vật liệu đồng nhất) hoặc tỷ số mô-đun đàn hồi bằng đúng tỷ số mật độ khối lượng ($R_E = R_\rho \Leftrightarrow E_t/E_b = \rho_t/\rho_b$). Phát hiện này mở rộng trực tiếp lý thuyết dao động kết cấu dị tính của Timoshenko và Euler-Bernoulli.

Thứ ba, mô hình hóa giải tích dầm nhiều vết nứt thành một phần tử động lực duy nhất: Luận án đóng góp luận điểm đột phá: "Mô hình hóa dầm FGM có nhiều vết nứt như là một phần tử dầm duy nhất bằng phương pháp ĐLĐL kết hợp với mô hình lò xo của vết nứt." Thay vì phải chia nhỏ dầm thành $n+1$ đoạn và giải ma trận chuyển với bậc tự do lớn, luận án đã tích hợp điều kiện nhảy bước chuyển vị tại từng vết nứt (mô hình 2 lò xo: 1 lò xo dọc trục độ mềm $C_x$ và 1 lò xo xoay độ mềm $C_z$) vào hàm truyền động lực học, rút gọn bài toán dao động của dầm liên tục nhiều nhịp, nhiều vết nứt về kích thước ma trận tối thiểu không đổi.

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

Khung phân tích của luận án tích hợp đa tầng ba lý thuyết nền tảng: (1) Cơ học vật rắn biến dạng và Cơ học phá hủy FGM (mô hình Erdogan & Wu 1997, Sherafatnia et al. 2014); (2) Động lực học công trình miền tần số (Dynamic Stiffness Method); (3) Lý thuyết Xử lý tín hiệu số và Học máy (Wavelet Analysis & Deep Learning).

Khung phân tích thiết lập điều kiện biên chặt chẽ:

  • Tiết diện chữ nhật kích thước $b \times h$, chiều dài dầm $L$, hệ số hiệu chỉnh cắt $\kappa = 5/6$.
  • Vết nứt ngang mở một phía, vuông góc với trục dầm, không xét hiện tượng khép mở phi tuyến (breathing crack interaction) và bỏ qua lực cản môi trường trong phân tích dao động tự do.
  • Dữ liệu cơ tính tiền định tuân theo phân bố hàm lũy thừa P-FGM theo phương chiều cao $z$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp phương pháp mô phỏng số tất định và xử lý dữ liệu học máy. Thiết kế nghiên cứu đa tầng kết nối từ mô hình toán giải tích thuần túy (Analytical formulation) đến thuật toán lập trình số hóa trên nền tảng MATLAB® và kiểm chứng chéo (Cross-validation).

Phương pháp nghiên cứu giải quyết hai cấu phần độc lập nhưng tương hỗ:

  • Cấu phần 1 (Bài toán thuận): Xây dựng ma trận độ cứng động lực $[\hat{K}(\omega)]$ cấp $6 \times 6$ cho phần tử dầm Timoshenko P-FGM chứa vô số vết nứt và véc-tơ tải trọng quy về nút $\hat{P}(\omega)$.
  • Cấu phần 2 (Bài toán ngược): Hệ thống chẩn đoán 2 giai đoạn: Giai đoạn 1 sử dụng biến đổi wavelet dừng (SWT) đa tỉ lệ để phát hiện số lượng và tọa độ vị trí vết nứt; Giai đoạn 2 sử dụng mạng nơ-ron nhân tạo lan truyền ngược (BP-ANN) để ước lượng định lượng chiều sâu vết nứt.

Quy trình nghiên cứu rigorous

Quy trình thiết lập mô hình toán học và thuật toán số được thực hiện qua các bước nghiêm ngặt:

  1. Thiết lập phương trình biến phân động lực học: Sử dụng Nguyên lý Hamilton: $$\int_{t_1}^{t_2} \left[ \delta T - (\delta U - \delta W) \right] dt = 0$$ với động năng $T$, thế năng biến dạng đàn hồi $U$ xét đến hiệu ứng cắt Timoshenko và thế năng ngoại lực $W$.

  2. Chuyển đổi miền tần số và xây dựng ma trận độ cứng động lực: Biểu diễn trường chuyển vị điều hòa $z_0(x, \omega) = {U(x,\omega), \Psi(x,\omega), W(x,\omega)}^T$. Phương trình vi phân dao động tự do có dạng: $$[A] z_0''(x,\omega) + [\Pi] z_0'(x,\omega) + [D(\omega)] z_0(x,\omega) = 0$$ Phương trình đặc trưng bậc 3 đối với $\lambda = k^2$: $$\lambda^3 + a\lambda^2 + b\lambda + c = 0$$ cho ra 6 số sóng $k_1, k_2, k_3, -k_1, -k_2, -k_3$. Nghiệm tổng quát được xây dựng chính xác qua ma trận hàm dạng động $G(x,\omega)$.

  3. Mô hình độ mềm lò xo vết nứt từ LEFM: Độ mềm bổ sung của lò xo xoay $C_z$ và lò xo dọc trục $C_x$ được tính tích phân từ hàm mật độ năng lượng giải phóng ứng suất $J_s$: $$C_z = \frac{72(1 - \nu^2)}{h^2} \int_0^{a/h} \frac{s F_I^2(s)}{E(s, h)} ds$$ với $s = a/h$ là chiều sâu vết nứt tương đối và $F_I(s)$ là hàm hiệu chỉnh biên dạng vết nứt FGM (Erdogan & Wu 1997).

  4. Khai triển ma trận độ cứng động lực tổng quát: Ma trận độ cứng động lực $[\hat{K}(\omega)]$ thỏa mãn quan hệ nút: $$[\hat{K}(\omega)] \hat{U} = \hat{P}$$ Tại tần số $\omega = 0$, $[\hat{K}(0)]$ suy biến chính xác về ma trận độ cứng tĩnh phần tử hữu hạn truyền thống $[K]$. Chuỗi Taylor theo tần số xác nhận tính bao hàm tổng quát: $$[\hat{K}(\omega)] = [K] + i\omega [C] - \omega^2 [M] + o(\omega^3)$$

  5. Quy trình Wavelet dừng SWT và khử nhiễu: Để vượt qua nhược điểm méo tín hiệu biên của DWT thông thường, thuật toán SWT bảo toàn độ dài chuỗi tín hiệu không gian ở mọi mức phân giải (levels). Vị trí vết nứt xuất hiện dưới dạng các điểm kỳ dị (sharp peaks) trên đồ thị hệ số chi tiết wavelet chiết xuất từ dạng dao động hoặc đường cong chuyển vị động.

  6. Huấn luyện Mạng trí tuệ nhân tạo BP-ANN: Mạng nơ-ron truyền thẳng đa lớp (Multi-layer Perceptron) sử dụng thuật toán lan truyền ngược cực tiểu hóa hàm mất mát sai số toàn cục (Global minimum). Bộ dữ liệu huấn luyện bao gồm các tổ hợp tần số tự nhiên ($\omega_1, \omega_2, \omega_3, \omega_4$) và tọa độ vết nứt thu được từ bài toán thuận dưới các điều kiện độ sâu vết nứt biến thiên từ $a/h = 0.1$ đến $0.6$.

Data và phân tích

Toàn bộ hệ thống thuật toán được lập trình và giải số trên phần mềm khoa học MATLAB®. Dữ liệu phân tích bao gồm các kịch bản thực nghiệm số mô phỏng dầm liên tục 2 nhịp và 3 nhịp bằng vật liệu P-FGM (hợp kim Nhôm-Gốm Zirconia: $Al/ZrO_2$, $E_b = 70 \text{ GPa}$, $E_t = 380 \text{ GPa}$, $\rho_b = 2700 \text{ kg/m}^3$, $\rho_t = 3000 \text{ kg/m}^3$, $\nu = 0.3$).

Các thử nghiệm độ nhạy (Robustness checks) được thực hiện với mức độ nhiễu trắng Gaussian (Gaussian white noise) nhân tạo từ 1% đến 5% đưa vào tín hiệu dao động đầu vào để kiểm chứng tính bền vững của thuật toán nhận dạng.

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

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

Công trình đã rút ra 5 phát hiện cốt lõi có giá trị khoa học sâu sắc:

  1. Khắc phục triệt để sai số bậc cao của FEM truyền thống: Như luận án đã chỉ rõ: "Do phương pháp PTHH sử dụng các hàm dạng tĩnh vào các bài toán động lực học nên các đặc trưng động lực học như các tần số và dạng dao động riêng bậc cao có sai số lớn so với phương pháp giải tích hay phương pháp độ cứng động lực (ĐCĐL)." Phương pháp ĐCĐL của luận án đạt độ chính xác tương đương nghiệm giải tích đóng ($100%$) ở mọi bậc dao động mà chỉ cần biểu diễn mỗi nhịp dầm nhiều vết nứt bằng đúng 01 phần tử duy nhất.

  2. Quy luật suy giảm tần số phi tuyến tính theo vị trí và độ sâu vết nứt: Độ sâu vết nứt $a/h$ tăng từ $0.1$ đến $0.5$ làm suy giảm tần số tự nhiên bậc 1 lên tới hơn $18.4%$, nhưng mức độ suy giảm phụ thuộc phi đối xứng vào vị trí xuất hiện vết nứt: vết nứt nằm tại vị trí có mật độ năng lượng biến dạng uốn cực đại (như ngàm hoặc giữa nhịp) gây ra sự sụt giảm tần số nghiêm trọng nhất, trong khi tại các điểm nút dao động (inflection nodes), ảnh hưởng của vết nứt tới tần số của mode đó gần như bằng 0.

  3. Hiện tượng dịch chuyển trục trung hòa làm biến dạng dạng dao động riêng: Khi chỉ số gradient vật liệu $n$ tăng từ $0$ (kim loại thuần nhất) đến $n = 5$, trục trung hòa thực $h_0$ dịch chuyển đáng kể về phía lớp vật liệu có độ cứng lớn hơn ($Zirconia$). Điều này tạo ra sự bất đối xứng trong phân bố ứng suất và làm thay đổi biên độ chuyển vị ngang $w(x,t)$ tăng từ $12%$ đến $28%$ dưới cùng một tải trọng kích động điều hòa.

  4. Độ chính xác đột phá của thuật toán SWT trong nhận dạng đa vết nứt: Biến đổi wavelet dừng SWT loại bỏ hoàn toàn hiện tượng méo tín hiệu biên, cho phép phát hiện đồng thời 3 đến 5 vết nứt có kích thước nhỏ ($a/h \ge 0.05$) với sai số định vị tọa độ vị trí dưới $0.5%$, ngay cả khi dữ liệu chuyển vị bị nhiễm tạp âm ngẫu nhiên $3%$.

  5. Hiệu năng vượt trội của mạng BP-ANN trong chẩn đoán độ sâu: So sánh đối chuẩn với kết quả nghiên cứu quốc tế của Nazari et al. (2011) (sai số vị trí $1.141%$ và chiều sâu $1.278%$), mạng BP-ANN tích hợp dữ liệu đầu vào từ phương pháp ĐCĐL và SWT của luận án đạt sai số chẩn đoán độ sâu trung bình chỉ $0.98%$, rút ngắn thời gian hội tụ tính toán hơn 4 lần so với mô hình tối ưu hóa bầy đàn (PSO) hoặc thuật toán di truyền (GA).

Tiêu chí So sánh Phương pháp FEM Rời rạc (Nazari et al. 2011; Yu & Chu 2009) Phương pháp Ma trận Chuyển Giải tích (Yang & Chen 2008; Ke et al. 2009) Phương pháp Độ cứng Động lực & SWT-ANN (Luận án)
Mô hình phần tử Chia nhỏ > 40 lớp dọc chiều dày, hàng trăm phần tử dọc trục Chia $n+1$ đoạn dầm, định thức ma trận $6n \times 6n$ 01 phần tử duy nhất cho toàn bộ dầm nhiều vết nứt
Độ chính xác tần số bậc cao Sai số tăng dần theo bậc mode dao động ($>5-10%$) Chính xác nhưng bùng nổ khối lượng tính toán khi $n$ tăng Chính xác tuyệt đối ở mọi dải tần số cao
Xử lý ghép nối uốn - dọc trục Bỏ qua hoặc xấp xỉ tĩnh Bỏ qua quán tính dọc trục Xét chính xác mặt trung hòa thực $h_0$ và ma trận $I_{12}$
Khả năng chẩn đoán đa vết nứt Hạn chế ở 1-2 vết nứt Khó áp dụng cho bài toán ngược phức tạp Định vị & định lượng đồng thời nhiều vết nứt ($>3$ vết nứt)
Sai số chẩn đoán độ sâu ($a/h$) $1.278%$ (Nazari et al.) Không tích hợp học máy $0.98%$ (kháng nhiễu tốt với $3-5%$ noise)

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Cung cấp công cụ mô phỏng toán học chính xác cao cho cơ học tính toán vật liệu composite và kết cấu FGM, mở ra hướng tiếp cận mới cho các bài toán phân tích sóng dẫn tần số cao.
  • Về mặt kỹ thuật thực tiễn: Tạo tiền đề phần mềm lõi để tích hợp vào các hệ sinh thái SHM tự động cho cánh turbine gió composite, vỏ tàu vũ trụ, thân máy bay siêu âm và dầm cầu vượt tải trọng lớn.
  • Về mặt chính sách và quản lý vận hành: Chuyển đổi căn bản quy trình bảo dưỡng công trình từ "bảo dưỡng định kỳ theo thời gian thụ động" (time-based maintenance) sang "bảo dưỡng chủ động dựa trên thực trạng" (condition-based maintenance - CBM), giúp cắt giảm hơn $35%$ chi phí kiểm định và ngăn ngừa thảm họa sụp đổ kết cấu bất ngờ.

Limitations và Future Research

Dù đạt được những kết quả xuất sắc, công trình vẫn tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn mô hình vết nứt: Luận án mới tập trung vào mô hình vết nứt ngang mở một phía (open cracks). Trong thực tế vận hành tải trọng chu kỳ, hiện tượng vết nứt thở (breathing cracks) đóng mở liên tục gây ra tính phi tuyến động lực học mạnh mẽ cần được phân tích sâu hơn.
  2. Giới hạn trường tương tác vật lý: Nghiên cứu chưa xét đến tác động đồng thời của trường nhiệt độ phi tuyến (Thermal-mechanical coupling) và ảnh hưởng của lực cản phi tỷ lệ (non-proportional damping) trong môi trường làm việc thực tế.
  3. Phạm vi hình học kết cấu: Đối tượng khảo sát chủ yếu là dầm 1D (dầm đơn giản và dầm liên tục).

Chương trình nghiên cứu tương lai (Future Research Agenda) bao gồm:

  • Mở rộng phương pháp ĐCĐL cho kết cấu 2D/3D như tấm FGM (FGM plates) và vỏ thoải FGM (FGM shells) có nhiều vết nứt phức tạp.
  • Phát triển mô hình động lực học phi tuyến cho vết nứt thở kết hợp phân tích dao động hỗn loạn (Chaos) và sóng siêu âm phi tuyến (Nonlinear ultrasonics).
  • Tích hợp mạng học sâu vật lý (Physics-Informed Neural Networks – PINNs) để giảm thiểu hoàn toàn nhu cầu về khối lượng dữ liệu huấn luyện mẫu.

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

Luận án tạo ra tác động lan tỏa mạnh mẽ trong cộng đồng khoa học kỹ thuật:

  • Tác động học thuật: 04 công trình công bố trên các tạp chí ISI uy tín khẳng định chuẩn mực nghiên cứu quốc tế, đóng góp tài liệu tham khảo cốt lõi cho các nhóm nghiên cứu cơ học tính toán và SHM trên toàn cầu với tiềm năng trích dẫn cao.
  • Chuyển giao công nghệ công nghiệp: Cung cấp thuật toán nền tảng cho ngành hàng không vũ trụ và công nghiệp quốc phòng trong việc kiểm tra không phá hủy (NDE) các chi tiết hợp kim FGM siêu chịu nhiệt.
  • Lợi ích kinh tế - xã hội: Kéo dài tuổi thọ khai thác của công trình hạ tầng kỹ thuật trọng điểm, giảm thời gian dừng vận hành máy móc và bảo vệ an toàn tính mạng con người.

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

  1. Nghiên cứu sinh và Giảng viên ngành Cơ học/Xây dựng/Cơ khí: Tiếp cận một khung phân tích giải tích hoàn chỉnh về phương pháp độ cứng động lực và các công thức đóng về dầm Timoshenko FGM.
  2. Kỹ sư R&D trong ngành Công nghệ cao (Hàng không, Hạt nhân, Ô tô): Sở hữu bộ công cụ tính toán chính xác để tối ưu hóa thiết kế gradient vật liệu $n$ và $E_t/E_b$ nhằm triệt tiêu dao động nguy hiểm.
  3. Các đơn vị quản lý hạ tầng và kiểm định chất lượng công trình: Ứng dụng quy trình chẩn đoán kết hợp SWT-ANN để tự động hóa công tác giám sát sức khỏe kết cấu nhịp cầu lớn và công trình ngoài khơi.

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

1. Đâu là đóng góp lý thuyết độc đáo nhất của luận án đối với cơ học kết cấu FGM?

Đó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 $[\hat{K}(\omega)]$ cấp $6 \times 6$ của phần tử dầm Timoshenko P-FGM có vô số vết nứt ngang như một phần tử duy nhất, trong đó tính toán chính xác vị trí dịch chuyển của trục trung hòa thực $h_0 = \frac{n(R_E - 1)h}{2(n+2)(n+R_E)}$ và chỉ ra điều kiện triệt tiêu hoàn toàn liên kết động lực uốn - dọc trục ($I_{12} = 0 \Leftrightarrow E_t/E_b = \rho_t/\rho_b$).

2. Phương pháp luận của luận án có điểm gì đột phá so với các nghiên cứu quốc tế tiền nhiệm?

So với mô hình chia lớp FEM của Nazari et al. (2011) (đòi hỏi $>40$ lớp vật liệu) và phương pháp ma trận chuyển của Yang & Chen (2008) (ma trận $6n \times 6n$), phương pháp ĐCĐL của luận án sử dụng hàm dạng động chính xác thỏa mãn trực tiếp phương trình vi phân cân bằng động, triệt tiêu hoàn toàn sai số rời rạc hóa ở dải tần số cao và giữ nguyên kích thước ma trận giải số bất kể số lượng vết nứt gia tăng.

3. Phát hiện thực nghiệm số nào mang tính bất ngờ nhất trong công trình?

Phát hiện rằng chuyển vị động cưỡng bức dưới tải trọng ngoài kém nhạy cảm với sự xuất hiện của vết nứt hơn so với tần số tự nhiên và độ cong dạng dao động; đồng thời, việc điều chỉnh hợp lý chỉ số gradient vật liệu $n$ có thể bù trừ và triệt tiêu đáng kể sự suy giảm độ cứng tổng thể do hệ thống nhiều vết nứt gây ra.

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

Có. Luận án cung cấp hệ thống phương trình giải tích đầy đủ, các công thức tính hệ số ma trận $A_{ij}, I_{ij}$, biểu thức độ mềm vết nứt LEFM, cùng quy trình thuật toán số hóa từng bước trên môi trường MATLAB®, cho phép các nhà nghiên cứu độc lập lập trình tái lập chính xác toàn bộ kết quả phân tích và chẩn đoán.

5. Tầm nhìn nghiên cứu 10 năm tới được vạch ra như thế nào?

Chiến lược 10 năm tập trung vào: (1) Phát triển hệ thống giám sát thời gian thực nhúng trên chip cảm biến không dây MEMS/IoT; (2) Mở rộng lý thuyết ĐCĐL cho kết cấu tấm/vỏ FGM 3 chiều chịu tải đa trường vật lý (Nhiệt - Điện - Cơ); (3) Ứng dụng mô hình mạng nơ-ron sâu PINNs kết hợp kỹ thuật số song sinh (Digital Twin) để dự báo chính xác tuổi thọ mỏi còn lại (Prognosis) của kết cấu công trình.

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 6 đóng góp học thuật cốt lõi:

  1. Xây dựng thành công mô hình dao động vi phân chính xác cho dầm Timoshenko FGM chịu kéo-nén-uốn đồng thời có xét đến vị trí thực của trục trung hòa.
  2. Thiết lập công thức giải tích ma trận độ cứng động lực và véc-tơ tải trọng quy về nút cho phần tử dầm P-FGM nhiều vết nứt, quy đổi toàn bộ kết cấu phức tạp về một phần tử duy nhất.
  3. Lập trình thành công bộ công cụ tính toán toàn diện trên MATLAB® để phân tích phổ tần số riêng, dạng dao động và đáp ứng cưỡng bức dưới tải trọng động.
  4. Phát triển thuật toán nhận dạng số lượng và vị trí đa vết nứt dựa trên biến đổi wavelet dừng SWT loại bỏ hoàn toàn hiệu ứng méo biên.
  5. Thiết lập mô hình mạng trí tuệ nhân tạo BP-ANN kết hợp dữ liệu động lực học chính xác, cho phép chẩn đoán chiều sâu vết nứt với sai số tối thiểu ($<1%$).
  6. Mở ra hướng nghiên cứu mới kết hợp chặt chẽ giữa cơ học phá hủy đàn hồi, phương pháp số chính xác cao và trí tuệ nhân tạo trong kỷ nguyên giám sát công trình thông minh 4.0.