Tổng quan luận án
Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (mã số: 62 52 01 01) của nghiên cứu sinh Nguyễn Ngọc Huyên, được thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn khoa học của GS. Nguyễn Tiến Khiêm và TS. Nguyễn Đình Kiên (năm 2017), có tiêu đề: “Phân tích dao động và chẩn đoán vết nứt dầm FGM” (với kiến nghị chỉnh sửa tên thành “Phân tích phổ dao động của dầm FGM có vết nứt”).
1. Tính cấp thiết và khoảng trống nghiên cứu
Vật liệu cơ tính biến thiên (Functionally Graded Materials - FGM) là loại vật liệu composite cấu thành chủ yếu từ hai thành phần (như kim loại và gốm) với tỷ lệ thể tích biến đổi liên tục theo không gian, giúp loại bỏ hiện tượng tập trung ứng suất và tách lớp thường gặp trong composite phân lớp truyền thống. Do đặc tính chịu nhiệt độ cao, chống mài mòn và độ bền cao, kết cấu dầm FGM được ứng dụng trong các ngành công nghệ cao như hàng không vũ trụ, hạt nhân và cơ khí chế tạo.
Trong quá trình vận hành, kết cấu FGM có thể xuất hiện các khuyết tật hoặc vết nứt cục bộ. Sự xuất hiện của vết nứt làm thay đổi các đặc trưng động lực học của hệ. Mặc dù các bài toán dao động của kết cấu FGM nguyên vẹn đã được khảo sát nhiều bằng phương pháp giải tích gần đúng hoặc phương pháp phần tử hữu hạn (PTHH), các phương pháp này thường chỉ giải quyết tốt bài toán ở tần số thấp hoặc dao động cơ bản, gặp khó khăn khi phân tích ở miền tần số cao. Trong khi đó, tác động của vết nứt thể hiện rõ nét hơn ở dải tần số cao. Đồng thời, các công trình nghiên cứu dao động của dầm FGM chứa vết nứt và bài toán ngược nhận dạng vết nứt dựa trên các thông số dao động còn rất hạn chế. Do đó, việc xây dựng phương pháp tiếp cận trong miền tần số để phân tích dao động và nhận dạng vết nứt trên dầm FGM là một yêu cầu khoa học và kỹ thuật cần thiết.
2. Mục tiêu nghiên cứu
Mục tiêu chính của luận án là xây dựng cách tiếp cận phổ để phân tích dao động và nhận dạng dầm FGM có vết nứt. Cụ thể, các mục tiêu nghiên cứu được xác định theo văn bản gồm:
- Xây dựng các phương trình cơ bản mô tả dao động của dầm FGM trong miền tần số, bao gồm dao động tự do, dao động cưỡng bức và thiết lập ma trận độ cứng động.
- Nghiên cứu sự tương tác giữa dao động dọc trục và dao động uốn, xác định điều kiện không tương tác và khảo sát một số đặc trưng truyền sóng trong dầm FGM.
- Xây dựng mô hình dầm FGM chứa một vết nứt dừng sử dụng mô hình hai lò xo tương đương (lò xo dọc trục và lò xo xoắn); phân tích ảnh hưởng của vết nứt cùng các tham số vật liệu và hình học đến tần số riêng của dầm.
- Đề xuất thuật toán giải tích để chẩn đoán vị trí và chiều sâu vết nứt trong dầm FGM từ các tần số riêng và kiểm nghiệm bằng tính toán số.
3. Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Kết cấu dầm phẳng Timoshenko làm từ vật liệu FGM có cơ tính biến thiên liên tục theo chiều dày tuân theo quy luật hàm số lũy thừa (P-FGM), có xét đến vị trí thực của trục trung hòa và chứa một vết nứt cạnh mở không phát triển (vết nứt dừng).
- Phạm vi nghiên cứu: Mô hình hóa một chiều cho dầm đàn hồi tuyến tính; khảo sát các dạng liên kết biên cơ bản (tựa đơn hai đầu, ngàm hai đầu, công-xôn); phân tích dao động trong miền tần số và giải bài toán ngược chẩn đoán khuyết tật đơn bằng số liệu tần số riêng.
Tổng quan tài liệu và vị trí của luận án
1. Các hướng nghiên cứu trước được luận án điểm lại
Luận án đã hệ thống hóa các công trình nghiên cứu trong nước và quốc tế liên quan đến dao động dầm FGM và chẩn đoán vết nứt theo các nhóm phương pháp:
| Hướng nghiên cứu / Phương pháp |
Tác giả, công trình và năm |
Nội dung và phát hiện chính |
| Phương pháp giải tích |
Li (2008) [40] |
Bỏ qua chuyển vị dọc trục ($u=0$), phương trình dao động dầm FGM tương tự dầm Timoshenko đồng chất; tần số riêng tăng dần theo số lớp và đạt trạng thái FGM khi số lớp tiến tới vô cùng. |
|
Sina và cộng sự (2009) [49] |
Thiết lập hệ phương trình chuyển động đầy đủ cho cả ba chuyển vị ($u, \theta, w$), nhận được phương trình đặc trưng tổng quát nhưng chưa tính đến vị trí thực của mặt trung hòa. |
| Phương pháp Rayleigh-Ritz |
Pradhan và Chakraverty (2013) [44] |
Thiết lập phương trình $[K - \omega^2 M]q = 0$, khảo sát 5 tần số riêng đầu tiên phụ thuộc vào độ mảnh $L/h$ và số mũ $n$ với các điều kiện biên khác nhau. |
|
Simsek (2010) [47, 48] |
Khảo sát dầm FGM (nhôm/ô-xit nhôm) chịu tải di động theo 3 lý thuyết: Euler-Bernoulli (EBT), Timoshenko (TBT), biến dạng trượt bậc 3 (TSDT). Khi $L/h=20$, độ võng ba lý thuyết tương đương; khi $L/h=5$, TBT và TSDT tương đương nhưng khác EBT. |
| Phương pháp phần tử hữu hạn (PTHH) |
Kosmatka (1993) [39] |
Xây dựng hàm dạng đa thức từ nghiệm tĩnh cho dầm FGM. |
|
Chakraborty và cộng sự (2003) [9] |
Dùng hàm dạng bậc hai cho $u, \theta$ và bậc ba cho $w$ đối với dầm Timoshenko FGM chịu nhiệt; phát hiện tần số cắt của dầm FGM nằm giữa tần số cắt của hai vật liệu thành phần. |
|
Eltaher và cộng sự (2012) [16] |
Xây dựng mô hình PTHH dầm Euler-Bernoulli FGM có tính đến vị trí thực của mặt trung hòa. |
|
Tại Việt Nam: Thái Hữu Tài, Võ Phương Thức (2012) [52]; Nguyễn Trung Kiên và cs (2013) [42]; Võ Phương Thức và cs (2013, 2014) [53, 54]; Nguyễn Đình Kiên và cs (2008–2014) [21, 29–36]; Lê Thị Hà (2016) [1] |
Ứng dụng lý thuyết trượt bậc cao, trượt bậc nhất, phân tích dầm sandwich FGM, ứng xử phi tuyến dưới tải trọng di động và dầm FGM tiết diện thay đổi. |
| Phương pháp độ cứng động & Phần tử phổ |
Su và Banerjee (2013, 2014) [50, 51] |
Dùng phương pháp độ cứng động (DSM) cho dầm Euler-Bernoulli và Timoshenko FGM; chỉ ra mode dọc trục và uốn tương tác nhau; tần số giảm đơn điệu theo số mũ $n$. |
|
Chakraborty và Gopalakrishnan (2003) [8] |
Dùng phương pháp phần tử phổ (SEM); phát hiện 3 dạng sóng đàn hồi (dọc, ngang, cắt), trong đó sóng cắt chỉ xuất hiện ở tần số lớn hơn tần số cắt; tần số cắt tăng theo số mũ $n$. |
| Mô hình vết nứt & Dao động dầm nứt |
Erdogan và Wu (1997) [18]; Jin và Batra (1996) [23]; Yang và Chen (2008) [58] |
Mô hình hóa vết nứt bằng lò xo xoắn; tính hệ số tập trung ứng suất $K_I$ bằng chuỗi Chebyshev hoặc hàm xấp xỉ cho E-FGM (Ke và cs, 2009 [38]) và FGM (Karaagac và cs, 2011 [25]). |
|
Yang, Chen và cs (2008, 2009) [58, 59] |
Thiết lập phương trình tần số dầm Euler-Bernoulli FGM có vết nứt bằng định thức cấp $(6n \times 6n)$ do bỏ qua quán tính dọc trục; phân tích độ nhạy tần số riêng theo vị trí, độ sâu vết nứt. |
|
Wei và cộng sự (2012) [56] |
Dùng ma trận truyền rút gọn phương trình tần số dầm Timoshenko FGM nứt về định thức cấp $(3 \times 3)$; biến dạng trượt ảnh hưởng rõ khi $L/h < 10$. |
|
Aydin (2013) [6] |
Biểu thức giải tích dạng dao động dầm Euler-Bernoulli FGM nứt dẫn đến định thức cấp 3. |
| Nhận dạng vết nứt trong dầm FGM |
Yu và Chu (2009) [60] |
Dùng p-version FEM và phương pháp đường mức tần số; với sai số đo 2%, xác định vị trí vết nứt sai số tối đa 7%, độ sâu sai số 15%, chỉ áp dụng cho vết nứt $a/h > 20%$. |
|
Banerjee và cộng sự (2015) [7] |
Dùng thuật toán di truyền tối ưu hóa hàm mục tiêu $G(a, e) = \sum [f_i(a, e) - f_i^*]^2 \to \min$ tìm giao điểm các đường mức tần số. |
2. Khoảng trống luận án lựa chọn giải quyết
Từ tổng quan tài liệu, luận án xác định các khoảng trống sau:
- Các nghiên cứu trước về dầm FGM có vết nứt đa số bỏ qua quán tính dọc trục hoặc tách rời dao động uốn và dọc trục một cách nhân tạo, chưa xem xét đầy đủ sự kết cặp dao động khi xét đến vị trí thực của trục trung hòa.
- Chưa có mô hình giải tích kết hợp đồng thời cả hai thành phần độ mềm cục bộ (lò xo dọc trục và lò xo xoắn) để mô tả vết nứt trong dầm FGM Timoshenko.
- Phương pháp nhận dạng vết nứt hiện hành (như đường mức tần số kết hợp thuật toán di truyền) phụ thuộc vào tính toán số nặng, trong khi chưa có thuật toán giải tích tường minh xác định trực tiếp vị trí và độ sâu vết nứt từ một số lượng hữu hạn tần số riêng đo được.
Cơ sở lý thuyết và phương pháp nghiên cứu
1. Cơ sở lý thuyết và mô hình toán học
- Lý thuyết dầm Timoshenko biến dạng trượt bậc nhất: Trường chuyển vị của điểm trên tiết diện dầm được biểu diễn qua chuyển vị trục trung hòa $u(x, t), w(x, t)$ và góc xoay $\theta(x, t)$:
$$u(x, z, t) = u(x, t) - (z - h_0)\theta(x, t), \quad w(x, z, t) = w(x, t)$$
trong đó $h_0$ là khoảng cách từ trục đối xứng hình học (mặt giữa) đến trục trung hòa thực tế.
- Quy luật vật liệu P-FGM: Cơ tính biến thiên theo chiều cao $z \in [-h/2, h/2]$:
$$E(z) = E_b + (E_t - E_b)\left(\frac{z}{h} + \frac{1}{2}\right)^n, \quad \rho(z) = \rho_b + (\rho_t - \rho_b)\left(\frac{z}{h} + \frac{1}{2}\right)^n$$
$$G(z) = \frac{E(z)}{2(1 + \nu)}$$
với $E_t, E_b$ và $\rho_t, \rho_b$ lần lượt là mô đun đàn hồi và mật độ khối ở mặt trên và mặt dưới; $n$ là chỉ số phân bố vật liệu ($n \ge 0$).
- Vị trí trục trung hòa ($h_0$): Được xác định từ điều kiện triệt tiêu hợp lực pháp tuyến tĩnh:
$$h_0 = \frac{\int_{-h/2}^{h/2} E(z)z , dz}{\int_{-h/2}^{h/2} E(z) , dz} = \frac{n(R_e - 1)h}{2(n + 2)(n + R_e)}$$
với $R_e = E_t / E_b$. Khi $R_e = 1$ (vật liệu đồng chất), $h_0 = 0$.
- Mô hình vết nứt hai lò xo tương đương: Vết nứt dừng tại vị trí $x = e$ có độ sâu $a$ được mô hình hóa bằng lò xo dọc trục có độ cứng $T$ và lò xo xoắn có độ cứng $R$. Điều kiện nhảy chuyển vị và liên tục nội lực tại vết nứt:
$$u(e+0) - u(e-0) = \frac{N(e)}{T}, \quad \theta(e+0) - \theta(e-0) = \frac{M(e)}{R}, \quad w(e+0) = w(e-0)$$
$$N(e+0) = N(e-0) = N(e), \quad M(e+0) = M(e-0) = M(e), \quad Q(e+0) = Q(e-0)$$
2. Phương pháp nghiên cứu
Luận án kết hợp nghiên cứu lý thuyết giải tích và mô phỏng số:
- Phương pháp giải tích trong miền tần số: Sử dụng biến đổi Fourier để chuyển hệ phương trình vi phân chuyển động từ miền thời gian sang miền tần số, biểu diễn dưới dạng phương trình ma trận vi phân cấp 2:
$$A z'' + B z' + C(\omega) z = q(x, \omega)$$
với $z = {U, \Theta, W}^T$. Nghiệm tổng quát được thiết lập qua hệ hàm cơ bản và các nghiệm đặc trưng $\lambda_j = i k_j$.
- Phương pháp độ cứng động (DSM): Thiết lập ma trận liên hệ giữa biên độ nội lực nút và chuyển vị nút trong miền tần số cho phần tử dầm FGM.
- Phương pháp hàm truyền (FRF): Xây dựng biểu thức tường minh cho các hàm đáp ứng tần số đơn ($h_{aa}, h_{bb}$) và hàm đáp ứng tần số chéo ($h_{ab}, h_{ba}$).
- Thuật toán giải tích chẩn đoán vết nứt: Biến đổi định thức ma trận đặc trưng thành hệ phương trình đại số chứa các thông số độ mềm vết nứt $\gamma_1, \gamma_2$, thiết lập phương trình phi tuyến một biến xác định vị trí vết nứt $e$ và giải bài toán cực tiểu hóa để tìm độ sâu $a$.
- Tính toán và kiểm chứng số: Lập trình tính toán tần số riêng, dạng dao động và đáp ứng tần số; so sánh đối chiếu kết quả với các công bố của Banerjee và cộng sự [7], Yang và Chen [58] cùng các lời giải PTHH chuẩn.
Nội dung chính theo từng chương
Mở đầu
Phần mở đầu trình bày bối cảnh nghiên cứu về vật liệu FGM, sự cần thiết của việc phân tích dao động dầm FGM trong miền tần số cao, các mục tiêu, đối tượng, phạm vi nghiên cứu, phương pháp tiếp cận và tóm tắt 4 đóng góp mới của luận án cùng danh mục các công bố liên quan.
Chương 1: Tổng quan về dao động của dầm FGM
Chương 1 hệ thống hóa các kiến thức cơ sở và tài liệu tổng quan:
- Vật liệu FGM: Khái niệm, phân loại theo hàm phân bố thể tích: P-FGM (hàm lũy thừa), S-FGM (hàm Sigmoid 2 miền), E-FGM (hàm số mũ $E(z) = E_t e^{\beta(1/2 - z/h)}$ với $\beta = \ln(E_t/E_b)$); các ứng dụng chịu nhiệt, chống mài mòn trong hàng không vũ trụ và máy móc.
- Lý thuyết dầm: So sánh các giả thiết phân bố biến dạng trượt $\phi(z)$: Euler-Bernoulli ($\phi=0$), Timoshenko ($\phi=z$), biến dạng trượt bậc 2 parabol ($\phi = z(1 - 4z^2/(3h^2))$), dạng e-mũ ($\phi = z e^{-2(z/h)^2}$) và lý thuyết bậc cao (HSDBT).
- Tổng quan các phương pháp: Đánh giá ưu nhược điểm của các phương pháp giải tích, Rayleigh-Ritz, phần tử hữu hạn, độ cứng động và phần tử phổ.
- Vết nứt và chẩn đoán: Tổng quan cơ học phá hủy FGM, các công thức tính hệ số tập trung ứng suất $K_I$, các nghiên cứu dao động dầm nứt của Yang, Chen, Wei, Aydin và phương pháp đường mức tần số nhận dạng vết nứt của Yu, Chu [60] và Banerjee [7].
Chương 2: Dao động của dầm FGM Timoshenko có vết nứt
Chương 2 tập trung thiết lập khung lý thuyết giải tích tổng quát:
- Phương trình chuyển động trong miền tần số: Áp dụng nguyên lý Hamilton thiết lập hệ phương trình chuyển động trong miền thời gian có xét đến các hệ số quán tính $I_{11}, I_{12}, I_{22}$ và độ cứng $A_{11}, A_{12}, A_{22}, A_{33}$. Dùng biến đổi Fourier đưa về hệ vi phân miền tần số $A z'' + B z' + C z = q$. Phương trình đại số đặc trưng bậc 3 đối với $\lambda^2$:
$$(\lambda^2 A_{33} + \omega^2 I_{11})(\lambda^2 A_{12} + \omega^2 I_{12})^2 - (\lambda^2 A_{11} + \omega^2 I_{11})[(\lambda^2 A_{33} + \omega^2 I_{11})(\lambda^2 A_{22} + \omega^2 I_{22}) - \omega^2 I_{11} A_{33}] = 0$$
luôn có 3 cặp nghiệm $\pm k_1, \pm k_2, \pm k_3$.
- Hàm đáp ứng tần số (FRF) và Ma trận độ cứng động: Nghiệm tổng quát bao gồm nghiệm thuần nhất và nghiệm riêng dạng tích phân ma trận truyền $H_a, H_b$. Xác lập các hàm đáp ứng đơn $h_{aa}(\omega), h_{bb}(\omega)$ và đáp ứng chéo $h_{ab}(\omega), h_{ba}(\omega)$. Khử véc-tơ hằng số để thiết lập ma trận độ cứng động $K(\omega)$ kích thước $(6 \times 6)$ cho phần tử dầm FGM.
- Bài toán dầm có vết nứt: Biểu diễn nghiệm dạng $z(x) = z_0(x) + K(x-e) z_0(e)$, áp dụng điều kiện biên thu được phương trình tần số dạng định thức:
$$\Delta(\omega) = \det[B_{LL}(\omega)] = d_0(\omega) + \gamma_1 d_1(e, \omega) + \gamma_2 d_2(e, \omega) + \gamma_1 \gamma_2 d_{12}(e, \omega) = 0$$
- Thuật toán chẩn đoán vết nứt: Khi đo được 3 tần số riêng $\omega_1, \omega_2, \omega_3$, hệ phương trình được chuyển thành:
$$a_{k1}(e)\gamma_1 + a_{k2}(e)\gamma_2 + a_{k3}(e)\gamma_3 = b_{0k} \quad (k=1,2,3)$$
với $\gamma_3 = \gamma_1 \gamma_2$. Giải hệ đại số thu được $\gamma_1(e) = \Delta_1/\Delta, \gamma_2(e) = \Delta_2/\Delta, \gamma_3(e) = \Delta_3/\Delta$. Từ hệ thức $\gamma_3 = \gamma_1 \gamma_2$, thiết lập phương trình phi tuyến xác định vị trí vết nứt $e$:
$$\Delta(e) \equiv \Delta_1(e)\Delta_2(e) - \Delta \Delta_3(e) = 0$$
Nghiệm $\hat{e}$ cho phép tính trực tiếp độ lớn vết nứt $\hat{\gamma}_1, \hat{\gamma}_2$, từ đó tính chiều sâu vết nứt $a$ thông qua bài toán cực tiểu hóa $[F_1(a) - \hat{\gamma}_1]^2 + [F_2(a) - \hat{\gamma}_2]^2 \to \min$. Đối với mô hình một lò xo xoắn ($\gamma_1 = 0$), thuật toán chỉ cần 2 tần số riêng.
Chương 3: Sự tương tác giữa dao động dọc trục và dao động uốn trong dầm FGM
Chương 3 đi sâu khảo sát bản chất vật lý của sự kết cặp dao động:
- Điều kiện không tương tác: Xác định hệ số tương tác quán tính $I_{12}$ và độ cứng $A_{12}$. Luận án chứng minh rằng khi tỷ số mật độ khối bằng tỷ số mô đun đàn hồi:
$$R_\rho \equiv \frac{\rho_t}{\rho_b} = R_e \equiv \frac{E_t}{E_b} \quad \implies \quad r = \frac{R_e}{R_\rho} = 1$$
thì đồng thời $A_{12} = 0$ và $I_{12} = 0$. Khi đó, dao động dọc trục và dao động uốn hoàn toàn tách rời nhau. Luận án định nghĩa dầm FGM thỏa mãn điều kiện này là dầm FGM tỷ lệ (Proportional FGM Beam - PGF Beam). Dầm đồng chất là trường hợp riêng của dầm PGF khi $R_e = R_\rho = 1$.
- Đặc trưng sóng và tần số cắt: Phân tích 3 nhánh sóng truyền trong dầm: sóng dọc trục, sóng uốn và sóng cắt. Xác định tần số cắt (cut-off frequency) phụ thuộc vào tham số vật liệu và độ mảnh $L/h$. Tần số cắt phân tách miền sóng dừng và sóng truyền, có giá trị biến thiên phụ thuộc vào số mũ $n$ và tỷ số $R_e, R_\rho$.
Chương 4: Kết quả số và thảo luận
Chương 4 thực hiện tính toán số toàn diện để kiểm chứng và đánh giá:
- Tần số riêng chuẩn hóa và dạng dao động: Tính toán 5 tần số riêng đầu tiên $\bar{\omega}_j = \omega_j \frac{L^2}{h}\sqrt{\frac{\rho_b}{E_b}}$ cho các trường hợp độ mảnh $L/h \in [5, 100]$ và số mũ $n \in [0, 10]$. Khảo sát sự tráo đổi thứ tự giữa mode uốn và mode dọc trục khi độ mảnh thay đổi.
- Hàm đáp ứng tần số (FRF): Vẽ đồ thị đáp ứng đơn ($h_{aa}, h_{bb}$) và đáp ứng chéo ($h_{ab}, h_{ba}$) theo dải tần số. Đồ thị đáp ứng chéo triệt tiêu khi $r = 1$ (dầm PGF), khẳng định tính đúng đắn của điều kiện tách rời dao động.
- So sánh kiểm chứng: Đối chiếu tỷ số tần số riêng của dầm ngàm hai đầu có vết nứt sâu $a/h = 30%$ với kết quả PTHH của Banerjee và cộng sự [7] tại $E_b/E_t = 0.2$ và $E_b/E_t = 5.0$, cho thấy sự trùng khớp chặt chẽ giữa đường giải tích của luận án và đường PTHH. So sánh tần số cơ bản với kết quả của Yang và Chen [58] với vết nứt $a/h = 20%$.
- Khảo sát tham số ảnh hưởng của vết nứt:
- Tần số riêng giảm khi độ sâu vết nứt $a/h$ tăng từ $5%$ đến $30%$.
- Mức độ suy giảm tần số phụ thuộc vào vị trí vết nứt $e/L$; độ suy giảm lớn nhất khi vết nứt nằm tại các vị trí có mô men uốn hoặc lực dọc cực đại (đầu ngàm), và bằng 0 tại các nút dao động tương ứng.
- Khảo sát sự phụ thuộc vào tỷ số mô đun đàn hồi $R_e$ và độ mảnh $L/h \in [5, 50]$.
- Thử nghiệm chẩn đoán vết nứt: Áp dụng thuật toán giải tích cho dầm FGM ngàm hai đầu có các độ mảnh và thông số vết nứt khác nhau, xác định chính xác vị trí $e$ và chiều sâu $a$ từ các tần số riêng.
Kết quả và những đóng góp mới
1. Những đóng góp mới về mặt lý luận và khoa học
Theo công bố của tác giả, luận án đạt được 4 đóng góp mới chính:
- Xây dựng lời giải giải tích dao động riêng và ma trận độ cứng động: Thiết lập biểu thức giải tích tường minh cho tần số riêng, dạng dao động, hàm đáp ứng tần số (FRFs đơn và chéo) và ma trận độ cứng động kích thước $(6 \times 6)$ cho dầm Timoshenko FGM biến thiên theo hàm lũy thừa, có tính đến vị trí thực của trục trung hòa ($h_0 \neq 0$).
- Tìm ra điều kiện tách rời dao động dọc trục và dao động uốn: Chứng minh điều kiện cần và đủ để dao động dọc trục và uốn không tương tác là $R_e = R_\rho$ ($r = 1$); đề xuất khái niệm dầm FGM tỷ lệ (PGF Beam) và làm sáng tỏ quy luật truyền sóng đàn hồi trong dầm FGM phụ thuộc các tham số vật liệu.
- Mô hình hóa vết nứt bằng hệ hai lò xo: Xây dựng mô hình dầm Timoshenko FGM chứa vết nứt dừng mô tả đồng thời bởi lò xo dọc trục ($T$) và lò xo xoắn ($R$), dẫn xuất phương trình tần số dạng định thức hiển không đòi hỏi chia nhiều đoạn phần tử phức tạp.
- Đề xuất phương pháp giải tích chẩn đoán vết nứt: Xây dựng thuật toán giải tích trực tiếp xác định vị trí và chiều sâu của một vết nứt trong dầm FGM bằng 3 tần số riêng (hoặc 2 tần số riêng cho mô hình 1 lò xo xoắn), loại bỏ sự phụ thuộc vào các thuật toán tối ưu hóa tìm kiếm lặp phức tạp.
2. Kiến nghị và giải pháp kỹ thuật
- Áp dụng ma trận độ cứng động đã thiết lập để mở rộng tính toán dao động cho các kết cấu khung dầm FGM phức tạp trong miền tần số cao.
- Sử dụng mô hình hai lò xo tương đương và thuật toán giải tích 3 tần số làm công cụ kiểm tra không phá hủy (NDT), phục vụ giám sát sức khỏe kết cấu dầm FGM trong kỹ thuật công trình và cơ khí chế tạo.
Hạn chế và hướng nghiên cứu tiếp
Dựa trên phạm vi và nội dung văn bản luận án trình bày, các hạn chế và định hướng phát triển gồm:
- Hạn chế:
- Mô hình nghiên cứu giới hạn ở bài toán một chiều (kết cấu dầm phẳng), chưa mở rộng trực tiếp cho kết cấu hai chiều như tấm, vỏ FGM có vết nứt.
- Vết nứt được giả thiết là vết nứt cạnh mở đơn lẻ, dạng vết nứt dừng (không xét quá trình phát triển mỏi hay lan truyền vết nứt phi tuyến).
- Quy luật biến đổi cơ tính tập trung chủ yếu vào hàm lũy thừa (P-FGM) theo một phương chiều dày.
- Hướng nghiên cứu tiếp được đề xuất:
- Mở rộng phương pháp độ cứng động và hàm truyền cho kết cấu tấm và vỏ FGM có chứa khuyết tật/vết nứt.
- Khảo sát dầm FGM có nhiều vết nứt đồng thời và vết nứt có góc nghiêng hoặc vết nứt thở (mở/đóng phi tuyến).
- Nghiên cứu thực nghiệm kiểm chứng trên mẫu dầm FGM thực tế để đánh giá ảnh hưởng của nhiễu đo đạc môi trường đến độ chính xác của thuật toán chẩn đoán.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị chuyên môn cao trong lĩnh vực Cơ học vật rắn biến dạng và Cơ kỹ thuật:
- Đối với nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học kết cấu: Tham khảo phương pháp thiết lập phương trình chuyển động dầm FGM có tính đến vị trí trục trung hòa, kỹ thuật biến đổi Fourier trong miền tần số, phương pháp ma trận độ cứng động và lý thuyết phần tử phổ.
- Đối với giảng viên và nhà nghiên cứu: Cung cấp hệ thống công thức giải tích chuẩn xác về tần số riêng, hàm đáp ứng tần số đơn/chéo, điều kiện phân tách dao động (dầm PGF) và nghiệm giải tích dầm nứt làm cơ sở đối chuẩn (benchmark) cho các phần mềm tính toán PTHH số.
- Đối với kỹ sư phân tích và chẩn đoán kỹ thuật: Ứng dụng thuật toán nhận dạng vết nứt từ tần số riêng để xây dựng quy trình đánh giá khuyết tật kết cấu composite biến tính trong thực tế.
Câu hỏi thường gặp
1. Dầm FGM tỷ lệ (PGF Beam) là gì và điều kiện để dao động uốn tách rời khỏi dao động dọc trục trong dầm FGM?
Theo Chương 3 của luận án, dầm FGM tỷ lệ (PGF Beam) là dầm FGM mà tại đó tỷ số mô đun đàn hồi ở hai mặt ($R_e = E_t / E_b$) bằng tỷ số mật độ khối ($R_\rho = \rho_t / \rho_b$), tức hệ số tỷ lệ $r = R_e / R_\rho = 1$. Khi điều kiện này thỏa mãn, cả hệ số độ cứng tương tác $A_{12}$ và hệ số quán tính tương tác $I_{12}$ đồng thời triệt tiêu ($A_{12} = 0, I_{12} = 0$), làm cho dao động dọc trục và dao động uốn hoàn toàn tách rời nhau (không tương tác). Dầm đồng chất là một trường hợp riêng của dầm PGF khi $R_e = R_\rho = 1$.
2. Trục trung hòa thực tế của dầm P-FGM được xác định như thế nào và có vị trí ra sao so với trục giữa?
Trục trung hòa thực tế nằm cách mặt giữa một khoảng $h_0$, được xác định từ điều kiện triệt tiêu ứng suất pháp tĩnh tích phân trên tiết diện:
$$h_0 = \frac{n(R_e - 1)h}{2(n + 2)(n + R_e)}$$
Khoảng cách $h_0$ phụ thuộc vào chiều cao dầm $h$, số mũ phân bố vật liệu $n$ và tỷ số mô đun đàn hồi $R_e$. Trục trung hòa chỉ trùng với mặt giữa ($h_0 = 0$) khi dầm đồng chất ($R_e = 1$) hoặc $n = 0$.
3. Vết nứt trong dầm FGM Timoshenko được mô hình hóa bằng những yếu tố nào?
Luận án mô hình hóa một vết nứt cạnh mở bằng hệ hai lò xo tương đương tập trung tại mặt cắt bị nứt ($x = e$): một lò xo dọc trục có độ cứng $T$ (đặc trưng cho độ mềm giãn trục) và một lò xo xoắn có độ cứng $R$ (đặc trưng cho độ mềm uốn). Độ lớn không thứ nguyên của vết nứt được biểu diễn qua $\gamma_1 = EA/T$ và $\gamma_2 = EI/R$.
4. Thuật toán giải tích của luận án cần bao nhiêu tần số riêng để xác định vị trí và chiều sâu vết nứt?
- Đối với mô hình tổng quát sử dụng hai lò xo (dọc trục và xoắn), thuật toán yêu cầu đo 3 tần số riêng ($\omega_1, \omega_2, \omega_3$) để lập hệ 3 phương trình đại số, từ đó khử độ lớn vết nứt để thu được một phương trình phi tuyến duy nhất xác định vị trí vết nứt $e$, tiếp theo xác định chiều sâu $a$ qua nghiệm cực tiểu sai phương.
- Đối với mô hình đơn giản hóa sử dụng một lò xo xoắn ($\gamma_1 = 0$), thuật toán chỉ cần 2 tần số riêng ($\omega_1, \omega_2$).
5. Biến dạng trượt trong lý thuyết Timoshenko ảnh hưởng rõ rệt nhất đến dao động dầm FGM trong điều kiện nào?
Theo các kết quả phân tích số và tổng quan trong luận án, ảnh hưởng của biến dạng trượt thể hiện rõ rệt khi dầm có độ mảnh thấp (dầm ngắn/dày với tỷ số $L/h < 10$). Khi độ mảnh lớn ($L/h \ge 20$), kết quả tính theo lý thuyết Timoshenko và lý thuyết Euler-Bernoulli tiệm cận gần nhau đối với tần số cơ bản.
Kết luận
Luận án của tác giả Nguyễn Ngọc Huyên đã xây dựng một khung phân tích giải tích hoàn chỉnh cho bài toán dao động và truyền sóng của dầm Timoshenko làm bằng vật liệu FGM có xét đến vị trí thực của trục trung hòa. Nghiên cứu đã xác lập điều kiện phân tách dao động cho dầm FGM tỷ lệ và đề xuất mô hình hai lò xo tương đương để mô tả vết nứt dừng. Đặc biệt, phương pháp giải tích chẩn đoán vết nứt sử dụng 3 tần số riêng do luận án phát triển đã được kiểm chứng số chặt chẽ, mang lại lời giải trực tiếp, chính xác cho bài toán nhận dạng khuyết tật kết cấu FGM trong miền tần số.
