Tổng quan về luận án

Sự ra đời của vật liệu composite có cơ tính biến thiên (Functionally Graded Materials - FGM) từ sáng kiến của các nhà khoa học tại Viện Sendai (Nhật Bản) vào năm 1984 đã mở ra bước đột phá cho các kết cấu kỹ thuật hoạt động trong môi trường khắc nghiệt như hàng không vũ trụ, thiết bị năng lượng hạt nhân và quốc phòng. Bằng cách kết hợp giữa pha gốm (độ cứng cao, chịu nhiệt và chống oxy hóa vượt trội) và pha kim loại (độ bền cơ học cao, độ dẻo dai lớn), FGM loại bỏ được mặt phân cách đột ngột vốn là điểm yếu gây tập trung ứng suất và bong tách lớp trong composite truyền thống. Tuy nhiên, như luận án đã chỉ rõ: "Do có sự tồn tại của các vật liệu khác nhau... nên trong vật liệu FGM dễ tồn tại các điểm, các vùng khuyết tật và dễ sinh ra các lỗ rỗng và vết nứt trong vật liệu". Khi kết cấu tấm FGM chịu tải trọng cơ học lớn hoặc trường nhiệt độ phi tuyến, sự hiện diện của các vết nứt làm suy giảm nghiêm trọng độ cứng cục bộ, đẩy nhanh quá trình mất ổn định (buckling) và làm biến đổi sâu sắc phổ tần số dao động riêng (natural frequency).

Khoảng trống nghiên cứu (research gap) cốt lõi trong y văn hiện đại nằm ở chỗ: phần lớn các nghiên cứu về kết cấu FGM chỉ tập trung vào các kết cấu lý tưởng nguyên vẹn (như các công trình của Praveen & Reddy [11], Yang et al. [8], Thai & Kim [10]). Khi xem xét bài toán nứt, các phương pháp số truyền thống như Phương pháp Phần tử Hữu hạn mở rộng (XFEM - Moës et al. [19], Liu et al. [17], Natarajan et al. [64]) hay Kỹ thuật đóng vết nứt ảo (VCCT - Krueger [18]) đòi hỏi các thuật toán phức tạp để theo dõi hình học gián đoạn của vết nứt, làm giàu hàm dạng hoặc chia lại lưới (remeshing) liên tục. Đặc biệt, bài toán ổn định và dao động tự do của tấm FGM có chiều dày thay đổi (variable thickness), chứa vết nứt xuyên tâm với góc nghiêng bất kỳ, đặt trên nền đàn hồi Pasternak và chịu trường nhiệt độ phi tuyến vẫn chưa được giải quyết một cách tường minh và toàn diện.

Để lấp đầy khoảng trống này, luận án tập trung giải quyết 4 câu hỏi nghiên cứu và các giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1: Mô hình toán học nào có thể mô phỏng chính xác sự suy giảm độ cứng cục bộ do vết nứt mà không cần tái tạo lưới phần tử hữu hạn phức tạp? (Giả thuyết $H_1$: Lý thuyết trường pha Phase-Field kết hợp với hàm suy giảm năng lượng biến dạng liên tục có thể mô tả chính xác vết nứt tùy ý).
  • Câu hỏi nghiên cứu 2: Tác động tương hỗ giữa hình học vết nứt (chiều dài $c$, góc nghiêng $\alpha$) và sự thay đổi chiều dày (tuyến tính, hàm mũ, parabol) ảnh hưởng như thế nào đến tải trọng mất ổn định tới hạn $P_{cr}$ và hệ số ổn định $k_{cr}$? (Giả thuyết $H_2$: Vết nứt định hướng vuông góc với trục nén làm suy giảm độ cứng mạnh nhất, nhưng góc nghiêng lớn sẽ giúp tái phân bố ứng suất và tăng tải tới hạn).
  • Câu hỏi nghiên cứu 3: Nền đàn hồi hai tham số Pasternak đóng vai trò gì trong việc bù trừ tổn thất độ cứng do khuyết tật nứt đối với động học tấm FGM? (Giả thuyết $H_3$: Hệ số kháng cắt $k_s$ của nền Pasternak có ảnh hưởng vượt trội so với hệ số đàn hồi Winkler $k_w$ trong việc nâng cao tần số dao động tự nhiên).
  • Câu hỏi nghiên cứu 4: Trường nhiệt độ phi tuyến ảnh hưởng như thế nào đến phổ dao động riêng của tấm FGM có nứt khi xét đến cơ tính phụ thuộc nhiệt độ? (Giả thuyết $H_4$: Gradient nhiệt độ làm mềm pha gốm nhanh chóng, dẫn đến sụt giảm tần số riêng phi tuyến theo chỉ số tỷ phần thể tích $n$).

Khung lý thuyết của luận án tích hợp chặt chẽ giữa: (1) Lý thuyết biến dạng cắt bậc ba của Shi (Third-order Shear Deformation Theory - TSDT [86]) giúp thỏa mãn hoàn toàn điều kiện ứng suất tiếp bằng không tại hai bề mặt tự do mà không cần hệ số hiệu chỉnh cắt; (2) Lý thuyết trường pha (Phase-Field Fracture Theory - Francfort & Marigo [74], Borden et al. [26]) chuyển đổi bài toán gián đoạn hình học thành bài toán biến phân liên tục với biến trường $s \in [0, 1]$; (3) Mô hình nền đàn hồi Pasternak hai tham số ($k_w, k_s$); và (4) Quy luật phân bố cơ tính dạng hàm lũy thừa (P-FGM) kết hợp truyền nhiệt phi tuyến 1D qua chiều dày.

Đóng góp đột phá của công trình đã được định lượng hóa và kiểm chứng nghiêm ngặt: Xây dựng thành công bộ thuật toán và chương trình tính toán số trên nền tảng mã nguồn mở FreeFem++-cs; giải quyết chính xác bài toán ổn định và dao động với độ sai lệch so với thực nghiệm của Seifi et al. [80] chỉ từ 0,39% đến 8,97%; kết quả dao động nhiệt tương đồng tuyệt đối với nghiệm giải tích 3D của Huang et al. [54] và Fazzolari [56]. Toàn bộ kết quả đã được công bố trên 9 công trình khoa học uy tín, bao gồm 5 bài báo thuộc danh mục SCIE/SCI Q1 (trong đó có các tạp chí quốc tế hàng đầu như Thin-Walled Structures với Impact Factor từ 5 đến 11). Phạm vi nghiên cứu bao trùm kết cấu tấm vuông và chữ nhật ($L \times H = 0,24m \times 0,24m$; $0,5m \times 0,5m$), độ mảnh $L/h$ từ 20 đến 200, cấu tạo từ hợp kim nhôm - gốm ($Al/Al_2O_3$) và thép không gỉ - silicon nitride ($SUS304/Si_3N_4$), khảo sát tỷ lệ vết nứt $c/L$ từ 0,1 đến 0,8, góc nghiêng $\alpha$ từ $0^\circ$ đến $90^\circ$, dưới 7 điều kiện biên khác nhau (SSSS, CCCC, SCSC, CFCF, SSSF, SCSF, SCSS) và dải nhiệt độ $T_c = 300K - 600K$.

Literature Review và Positioning

Các công trình nghiên cứu về tấm FGM chịu nén và dao động đã phát triển mạnh mẽ qua ba trường phái lý thuyết chính:

  1. Trường phái phân tích ổn định và động lực học tấm FGM lý tưởng: Praveen & Reddy [11] đã đặt nền móng với phân tích nhiệt-đàn phi tuyến dựa trên lý thuyết biến dạng cắt bậc nhất (FSDT). Tiếp theo, Reddy [12] và Thai et al. [9, 10] phát triển các lời giải dạng khép kín (closed-form solutions) và lý thuyết cắt bậc cao (HSDT) cho tấm trên nền đàn hồi Pasternak. Nguyễn Đình Đức và cộng sự [33, 63] đã mở rộng phân tích dao động phi tuyến cho các tấm FGM có gân gia cường trong môi trường nhiệt độ phức tạp.
  2. Trường phái mô hình hóa vết nứt trong kết cấu tấm: Moës, Dolbow & Belytschko [19] phát triển XFEM cho phép mô phỏng vết nứt không cần chia lại lưới; Natarajan et al. [64] ứng dụng XFEM cho tấm FGM nứt; Loc Tran et al. [68] sử dụng phương pháp phân tích đẳng hình học mở rộng (XIGA); Huang et al. [65, 73] áp dụng phương pháp Ritz 3D; Joshi et al. [69] và Ankur Gupta et al. [70, 71] dùng lý thuyết giải tích kết hợp nhiễu loạn đa thang đo. Tuy nhiên, các kỹ thuật này đều gặp trở ngại lớn khi vết nứt có hình học phức tạp, phân nhánh hoặc khi chiều dày tấm biến thiên liên tục.
  3. Trường phái cơ học phá hủy trường pha (Phase-Field Fracture Mechanics): Khởi nguồn từ công trình biến phân của Francfort & Marigo [74] và định dạng số của Miehe, Hofacker & Ulmer [22, 75], Borden et al. [26, 76] đã chứng minh tính ưu việt của Phase-Field trong việc mô phỏng vết nứt động và phân nhánh. Đoàn Hồng Đức, Nguyễn Đình Đức và cộng sự [20, 21] đã tiên phong áp dụng Phase-Field cho bài toán nứt động trong epoxy chứa hạt thủy tinh biến thiên.

Trong y văn tồn tại hai luồng quan điểm tranh luận sâu sắc: Luồng quan điểm thứ nhất (ủng hộ phương pháp gián đoạn như XFEM/VCCT) cho rằng cần duy trì sự gián đoạn hình học sắc nét tại bề mặt vết nứt để thu được hệ số cường độ ứng suất ($K_I, K_{II}$) chính xác; luồng quan điểm thứ hai (ủng hộ lý thuyết trường pha) lập luận rằng việc xấp xỉ vết nứt bằng một trường suy giảm năng lượng biến thiên trơn (regularized diffuse crack topology) không chỉ loại bỏ kỳ dị toán học mà còn tối ưu hóa năng lượng toàn cục, đặc biệt thích hợp cho các bài toán phân tích trị riêng (eigenvalue problems) như mất ổn định và dao động tự do.

Vị trí học thuật của luận án được xác lập rõ ràng tại điểm giao thoa đột phá: Luận án là công trình tiên phong tích hợp thành công lý thuyết trường pha Phase-Field với lý thuyết tấm biến dạng cắt bậc ba của Shi [86] để khảo sát đồng thời bài toán ổn định và dao động tự do của tấm FGM có chiều dày thay đổi tùy ý, tương tác với nền đàn hồi và trường nhiệt độ phi tuyến.

So sánh đối chuẩn quốc tế:

  • So với nghiên cứu thực nghiệm của Seifi & Nafiseh [80] trên tấm đồng chất nứt chịu nén: Luận án dự báo lực tới hạn $P_{cr}$ với độ chính xác cao (mẫu nứt $c/H=0,3, \alpha=0^\circ$ sai lệch chỉ 0,39%; mẫu $c/H=0,5, \alpha=30^\circ$ sai lệch 1,21%), đồng thời giải thích được nguyên nhân sai số thực nghiệm 11,96% - 13,02% do dung sai chế tạo chiều dày $\pm 50\mu m$.
  • So với các nghiên cứu dao động nhiệt FGM của Huang et al. [54] và Fazzolari [56]: Tham số tần số không thứ nguyên $\bar{\omega}$ của luận án tại $T_c=400K, n=1$ đạt 6,2762, tiệm cận hoàn hảo với giá trị 6,2755 của Huang [54] và Fazzolari [56], khẳng định độ tin cậy tuyệt đối của thuật toán phần tử hữu hạn đề xuấ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 trực tiếp lý thuyết biến dạng cắt bậc ba của Shi [86] (TSDT) sang bài toán cơ học phá hủy trường pha và môi trường nhiệt - đàn hồi phức hợp. Cụ thể:

  • Mở rộng lý thuyết Shi [86] cho bài toán chiều dày biến thiên $h(x)$: Các phương trình động học chuyển vị $u(x,y,z,t), v(x,y,z,t), w(x,y,z,t)$ được dẫn xuất đầy đủ có xét đến các đạo hàm bậc nhất của chiều dày $\partial h/\partial x$, giải quyết triệt để hiện tượng nghẽn cắt (shear locking) mà không cần đưa vào các hệ số hiệu chỉnh cắt nhân tạo như trong lý thuyết Mindlin.
  • Thiết lập mô hình biến phân năng lượng liên hợp: Tích hợp hàm suy giảm trường pha $g(s) = s^2$ vào mật độ thế năng đàn hồi của tấm, phản ánh sự biến đổi liên tục từ trạng thái nguyên vẹn ($s=1$) qua trạng thái mềm hóa ($0 < s < 1$) đến trạng thái phá hủy hoàn toàn ($s=0$).

Mô hình lý thuyết được chuẩn hóa qua hệ mệnh đề:

  • Mệnh đề $P_1$: Trạng thái mất ổn định tĩnh của tấm FGM nứt được xác định qua bài toán biến phân cực tiểu hóa phiếm hàm năng lượng toàn phần $\Pi(d, s)$, dẫn đến phương trình trị riêng suy biến: $$(K^e + \lambda K_G^e)d = 0$$ trong đó $K^e$ là ma trận độ cứng phụ thuộc biến trường pha $s$, $K_G^e$ là ma trận độ cứng hình học và $\lambda$ là tải trọng tới hạn.
  • Mệnh đề $P_2$: Trường pha vết nứt $s(x,y)$ được chi phối bởi phương trình Euler-Lagrange dừng trên toàn miền $\Omega$: $$-\epsilon^2 \nabla^2 s + s = 1 - \frac{2\epsilon}{\pi G_c} \Psi(d)$$ với $\epsilon$ là tham số độ rộng vùng nứt, $G_c$ là tốc độ giải phóng năng lượng tới hạn Griffith, và $\Psi(d)$ là hàm mật độ năng lượng biến dạng nứt.
  • Mệnh đề $P_3$: Dao động riêng của tấm nứt trong môi trường nhiệt và nền đàn hồi là nghiệm của hệ phương trình động lực học Lagrange loại 2: $$(K^e - \omega^2 M^e)d = 0$$ trong đó ma trận khối lượng $M^e$ và ma trận độ cứng $K^e$ suy giảm trực tiếp theo phân bố của $s(x,y)$.

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

Khung phân tích của luận án là sự hợp nhất của 4 trụ cột lý thuyết:

  1. Cơ học vật liệu chức năng biến thiên (Power-law P-FGM): Tỷ phần thể tích gốm tuân theo $V_c(z) = (z/h + 1/2)^n$, với cơ tính hiệu dụng phụ thuộc nhiệt độ phi tuyến $P(z,T) = (P_c(T) - P_m(T))V_c(z) + P_m(T)$.
  2. Cơ học phá hủy trường pha liên tục: Thay thế đường gián đoạn hình học $\Gamma$ bằng tích phân thể tích trên miền $\Omega$ với độ dài điều hòa $\epsilon$. Tỷ số $H/\epsilon = 400$ được chứng minh bằng thực nghiệm số là giá trị tối ưu đảm bảo tính hội tụ và độc lập lưới.
  3. Lý thuyết nền đàn hồi Pasternak hai tham số: Năng lượng biến dạng tương tác nền gồm thành phần biến dạng đứng $k_w w^2$ và thành phần cắt $k_s [(\partial w/\partial x)^2 + (\partial w/\partial y)^2]$.
  4. Lý thuyết dẫn nhiệt phi tuyến 1D trạng thái dừng: Phân bố nhiệt độ $T(z)$ qua chiều dày được giải chính xác bằng chuỗi đa thức bậc cao thỏa mãn điều kiện biên nhiệt độ $T(h/2) = T_c$ và $T(-h/2) = T_m$.

Điều kiện biên và giới hạn áp dụng (Boundary conditions): Luận án thiết lập tường minh điều kiện hình học và tự nhiên cho 7 trường hợp liên kết biên (SSSS, CCCC, SCSC, CFCF, SSSF, SCSF, SCSS), áp dụng cho tấm mỏng đến tấm dày trung bình ($L/h = 20 - 200$), giả định biến dạng đàn hồi tuyến tính nhỏ trước khi mất ổn định, vết nứt xuyên suốt chiều dày tấm.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi quan điểm thực chứng (positivism) và cơ học tính toán định lượng (quantitative computational mechanics). Mô hình hóa đa cấp độ (multi-level design) được triển khai:

  • Cấp vi mô/vật liệu: Tính toán mô đun Young $E(z,T)$, hệ số Poisson $\nu(z,T)$, khối lượng riêng $\rho(z)$, hệ số giãn nở nhiệt $\alpha(z,T)$ và hệ số dẫn nhiệt $K(z,T)$ của từng lớp vi phân qua chiều dày $z \in [-h/2, h/2]$.
  • Cấp trung mô (meso): Xây dựng trường suy giảm độ cứng cục bộ thông qua biến trường pha $s(x,y)$.
  • Cấp vĩ mô (macro): Giải bài toán biên giá trị riêng để xác định các mode mất ổn định và mode dao động tự do toàn kết cấu.

Vật liệu khảo sát chính xác bao gồm:

  • Cặp $Al/Al_2O_3$: $E_m = 70 \text{ GPa}, \nu_m = 0,3, \rho_m = 2707 \text{ kg/m}^3$; $E_c = 380 \text{ GPa}, \nu_c = 0,3, \rho_c = 3800 \text{ kg/m}^3$.
  • Cặp $SUS304/Si_3N_4$: Các hệ số nhiệt độ phi tuyến $P_0, P_{-1}, P_1, P_2, P_3$ được xác định chính xác theo công thức hàm phi tuyến ở Bảng 1.2 của luận án (tại $T=300K$: $E_m = 207,79 \text{ GPa}, E_c = 322,27 \text{ GPa}$; $\alpha_m = 15,32 \times 10^{-6} \text{ K}^{-1}, \alpha_c = 5,87 \times 10^{-6} \text{ K}^{-1}$).

Quy trình nghiên cứu rigorous

Quy trình tính toán số được thiết lập tự động hóa hoàn toàn trên ngôn ngữ phần tử hữu hạn FreeFem++-cs:

  1. Rời rạc hóa miền tính toán: Sử dụng phần tử tam giác đẳng tham số 6 nút (isoparametric 6-node triangular element). Mỗi nút mang 5 bậc tự do: $${d_i} = {u_{0i}, v_{0i}, w_{0i}, \theta_{xi}, \theta_{yi}}^T$$ Tổng số bậc tự do của mỗi phần tử là 30.
  2. Chiến lược chia lưới thích nghi cục bộ (Adaptive Mesh Refinement): Toàn bộ tấm được chia lưới tam giác Delaunay chuẩn, riêng vùng lân cận vết nứt bán kính $r \le 2\epsilon$ được làm giàu mật độ phần tử với kích thước phần tử cực mịn ($h_e \le \epsilon/2$).
  3. Khảo sát độ hội tụ và độ tin cậy (Bảng 3.1b): Khảo sát lưới từ 761 phần tử đến 5942 phần tử cho tấm vuông $L=H=0,24m, h=0,012m, c/H=0,3, \alpha=0^\circ$:
    • 761 phần tử: $P_{cr} = 1537,99 \text{ N}$
    • 2973 phần tử: $P_{cr} = 1537,68 \text{ N}$
    • 4129 phần tử: $P_{cr} = 1537,45 \text{ N}$
    • 4944 phần tử: $P_{cr} = 1537,28 \text{ N}$
    • 5487 phần tử: $P_{cr} = 1537,12 \text{ N}$
    • 5942 phần tử: $P_{cr} = 1537,01 \text{ N}$ Lưới 4944 phần tử (với lưới biên $10 \times 10$ và làm giàu tại vết nứt) được lựa chọn tối ưu giữa độ chính xác và chi phí tính toán.

Data và phân tích

Thuật toán giải bài toán mất ổn định và dao động trải qua 4 bước đồng bộ:

  • Bước 1: Giải phương trình vi phân đạo hàm riêng elliptic để tìm phân bố trường pha $s(x,y)$ từ hình học vết nứt $c, \alpha$ thông qua hàm khoảng cách $d(x,y)$ và hằng số suy biến $B = 10^5$.
  • Bước 2: Tích phân số Gauss bậc cao trên từng phần tử để lập ma trận độ cứng suy giảm $K^e(s)$, ma trận hình học $K_G^e(s)$ và ma trận khối lượng $M^e(s)$.
  • Bước 3: Áp đặt các điều kiện biên động học chặt chẽ (áp đặt khử bậc tự do tại các nút biên).
  • Bước 4: Sử dụng thuật toán lặp không gian con (Subspace Iteration) kết hợp phương pháp Lanczos trong FreeFem++-cs để trích xuất chính xác 5 giá trị riêng nhỏ nhất ($\lambda_1 - \lambda_5$ hoặc $\omega_1 - \omega_5$) và các vector riêng tương ứng mô tả dạng mất ổn định/dạng dao động riêng.

Kiểm tra độ vững chắc (Robustness checks): Khảo sát tỷ số $H/\epsilon$ trong dải từ 100 đến 800 (Bảng 3.2). Kết quả cho thấy khi $H/\epsilon \ge 400$, hệ số ổn định $k_{cr}$ hoàn toàn tiệm cận và không đổi, loại bỏ nguy cơ phụ thuộc kích thước lưới phần tử.

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

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

  1. Quy luật suy giảm tải trọng mất ổn định theo chiều dài và góc nghiêng vết nứt (Bảng 3.1a, Hình 3.1):

    • Khi vết nứt nằm ngang ($\alpha = 0^\circ$, vuông góc phương nén), tăng chiều dài vết nứt $c/H$ từ 0,1 lên 0,5 làm lực mất ổn định tới hạn của tấm đồng chất giảm mạnh từ 1773 N xuống 1205 N (suy giảm 32,04%).
    • Khi tăng góc nghiêng vết nứt $\alpha$ từ $0^\circ$ lên $60^\circ$ đối với mẫu nứt lớn $c/H = 0,5$, lực tới hạn tăng vọt từ 1205 N lên 1723 N (tăng phục hồi 42,99%). Giải thích lý thuyết: Vết nứt nghiêng làm giảm diện tích mặt cắt trực giao chịu nén trực tiếp, chuyển hướng dòng truyền ứng suất chính dọc theo trục nén và hạn chế giải phóng năng lượng biến dạng nứt dạng Mode I, giúp kết cấu duy trì độ cứng uốn tổng thể.

  2. Ảnh hưởng phi tuyến của chỉ số thể tích FGM $n$ và kích thước nứt lên tần số dao động tự do (Hình 4.4, Hình 4.11):

    • Tần số dao động tự nhiên giảm phi tuyến rất nhanh trong dải $n \in [0, 2]$, sau đó giảm chậm dần khi $n > 2$. Khi $n$ tăng, hàm lượng kim loại mềm dẻo có mật độ khối lượng lớn ($\rho_m$) tăng lên trong khi hàm lượng gốm cứng nhẹ ($\rho_c, E_c$) giảm đi, làm giảm tỷ số độ cứng trên khối lượng ($K/M$).
    • Sự kết hợp đồng thời giữa tăng chỉ số $n$ và tăng chiều dài vết nứt $c/L$ tạo ra hiệu ứng suy giảm cộng hưởng kép đối với tần số riêng.

  3. Hiệu ứng tương tác đột phá với nền đàn hồi Pasternak (Bảng 4.8, Hình 4.8):

    • Đối với tấm FGM có chiều dày thay đổi đặt trên nền Pasternak ($k_w = 50, k_s = 10$), ngay cả khi vết nứt phát triển lớn ($c/L = 0,8$), tần số dao động tự do vẫn luôn cao hơn tần số của tấm nguyên vẹn không đặt trên nền.
    • Khi không có nền đàn hồi, tấm có chiều dày thay đổi có tần số nhỏ hơn tấm chiều dày không đổi ($h_{const}$) khoảng 4%; nhưng khi đặt trên nền đàn hồi, tần số của tấm chiều dày thay đổi lại lớn hơn tấm chiều dày không đổi tới gần 30%.
    • Tham số kháng cắt $k_s$ của nền Pasternak có vai trò áp đảo tham số Winkler $k_w$ trong việc khống chế chuyển vị ngang và nâng cao độ cứng uốn.

  4. Sự suy giảm động học trong môi trường nhiệt độ cao (Bảng 4.10, Hình 4.12):

    • Khi nhiệt độ bề mặt gốm $T_c$ tăng từ 300K lên 600K ($T_m = 300K$), ứng suất nén nhiệt dư sinh ra trong mặt phẳng kết hợp với sự suy giảm mô đun Young $E(T)$ làm tần số dao động riêng giảm sâu.
    • Đối với tấm có chiều dày biến thiên dạng parabol đối xứng $h(x) = h_0(4\beta\xi^2 - 4\beta\xi + 1)$, mặt cắt tại tâm ($\xi = 0,5$) có chiều dày nhỏ nhất $h_{mid}$ là vị trí xung yếu nhất. Khi vết nứt định hướng trùng với mặt cắt này, tần số dao động giảm mạnh nhất.

  5. Hiệu ứng sụp đổ độ cứng đa thông số (Bảng 4.11):

    • Khi đồng thời gia tăng tỷ số kích thước mặt bằng $L/H$, tỷ số suy giảm chiều dày tâm $h_0/h_{mid}$ và chiều dài vết nứt $c/L$, độ cứng tổng thể của tấm bị triệt tiêu nghiêm ngặt, khiến tần số dao động tự do sụt giảm trên 12 lần.

Implications đa chiều

  • Đóng góp cho lý thuyết: Xác lập tính đúng đắn của việc kết hợp lý thuyết biến dạng cắt bậc ba TSDT của Shi và lý thuyết trường pha Phase-Field trong phân tích trị riêng kết cấu composite dị tính phức tạp, cung cấp nền tảng toán học chuẩn xác cho cơ học phá hủy kết cấu tấm.
  • Đổi mới phương pháp luận: Thiết lập quy trình tính toán số hiệu năng cao, ổn định và tự động hóa trên FreeFem++-cs, mở ra phương pháp giải quyết các bài toán nứt đa vết nứt, nứt cong, nứt phân nhánh mà không vấp phải rào cản chia lại lưới.
  • Ứng dụng thực tiễn trong công nghiệp: Cung cấp công cụ chẩn đoán kỹ thuật không phá hủy (NDT) cho kỹ sư hàng không và năng lượng. Bằng cách đo đạc phổ tần số dao động riêng và tải trọng mất ổn định, kỹ sư có thể nhận diện và đánh giá mức độ nguy hiểm của vết nứt tiềm ẩn trong các tấm vách FGM chịu nhiệt.
  • Khuyến nghị thiết kế và chính sách tiêu chuẩn: Đề xuất bổ sung hệ số suy giảm độ cứng cục bộ do khuyết tật nứt vào các tiêu chuẩn thiết kế kết cấu FGM chịu nhiệt; khuyến cáo tăng cường các liên kết biên ngàm (CCCC) hoặc bố trí các lớp đệm đàn hồi kháng cắt để hạn chế tối đa nguy cơ mất ổn định sớm.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Mô hình trường pha đẳng hướng thuần túy chưa xét đến hiệu ứng đóng vết nứt phi tuyến (unilateral contact / crack closure effect) và ma sát giữa hai bờ vết nứt khi chịu tải trọng nén chu kỳ phức tạp.
  2. Giả thiết trường nhiệt độ phân bố 1D dừng qua chiều dày $z$, chưa xét đến gradient nhiệt độ 2D/3D cục bộ do hiện tượng cản trở truyền nhiệt tại chính khe nứt.
  3. Giả định vật liệu làm việc trong miền đàn hồi tuyến tính và phá hủy giòn, chưa mô phỏng vùng dẻo hóa cục bộ (plastic zone) tại đầu vết nứt trong các lớp giàu kim loại của tấm FGM.
  4. Nghiên cứu thực hiện trên kết cấu tấm phẳng, chưa mở rộng trực tiếp sang các kết cấu vỏ cong FGM (spherical shells, cylindrical panels).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Hướng 1: Phát triển mô hình Phase-Field dị tính bất đối xứng (anisotropic phase-field) kết hợp vùng dẻo hóa kim loại để mô phỏng chính xác sự lan truyền nứt động và phân nhánh trong tấm FGM chịu tải trọng va đập, xung kích nổ.
  • Hướng 2: Nghiên cứu ổn định và dao động của tấm FGM có gân gia cường nghiêng (oblique stiffeners) và tấm FGM rỗng xốp (porous FGM) chứa vết nứt.
  • Hướng 3: Mở rộng mô hình lý thuyết cho kết cấu vỏ thoải hai độ cong, vỏ trụ FGM có nứt đặt trên nền đàn hồi phi tuyến.
  • Hướng 4: Tích hợp vật liệu thông minh áp điện (piezoelectric) và vật liệu gia cường nano (graphene/CNTs) để điều khiển chủ động dao động và ngăn chặn sự phát triển của vết nứt.
  • Hướng 5: Thiết kế và thực hiện các thí nghiệm đo đạc dao động nhiệt độ cao và chụp ảnh cắt lớp quang học số (DIC/X-ray tomography) trên mẫu FGM chế tạo bằng công nghệ in 3D để kiểm chứng toàn diện mô hình số.

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

  • Tác động học thuật: Luận án đã công bố 5 bài báo trên các tạp chí SCIE/SCI Q1 hàng đầu thuộc nhà xuất bản Elsevier (Thin-Walled Structures) và Taylor & Francis. Đây là những đóng góp học thuật quan trọng, được cộng đồng cơ học tính toán quốc tế trích dẫn và phát triển rộng rãi trong các nghiên cứu về kết cấu tiên tiến.
  • Chuyển đổi công nghiệp: Cung cấp giải pháp tính toán mô phỏng trực tiếp cho ngành công nghiệp hàng không vũ trụ (vỏ tàu con thoi, cánh tên lửa), công nghiệp quốc phòng và các nhà máy điện hạt nhân, giúp tối ưu hóa khối lượng kết cấu mà vẫn đảm bảo hệ số an toàn phá hủy khắt khe.
  • Tác động tiêu chuẩn kỹ thuật: Tạo lập cơ sở dữ liệu số chuẩn mực phục vụ việc xây dựng quy chuẩn an toàn cho vật liệu composite biến thiên chịu nhiệt độ cao tại Việt Nam và quốc tế.
  • Lợi ích xã hội: Nâng cao độ bền và độ tin cậy của các công trình công nghệ cao, giảm thiểu nguy cơ tai nạn thảm khốc do phá hủy kết cấu đột ngột, thúc đẩy làm chủ công nghệ vật liệu tiên tiến tại Việt Nam.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận hệ thống phương trình vi phân biến phân hoàn chỉnh, mã nguồn FreeFem++-cs mẫu và phương pháp luận tích hợp TSDT với Phase-Field cho các bài toán cơ học vật rắn biến dạng.
  • Giảng viên và nhà nghiên cứu cao cấp: Sử dụng các phát hiện đột phá về tương tác nền đàn hồi - nhiệt độ - vết nứt làm tài liệu tham khảo chuyên sâu và phát triển các hướng nghiên cứu liên ngành.
  • Kỹ sư R&D trong công nghiệp: Ứng dụng các biểu đồ suy giảm lực tới hạn $P_{cr}$ và tần số $\bar{\omega}$ để thiết lập ngưỡng cảnh báo hư hại và quy trình bảo trì dự đoán cho kết cấu tấm FGM.
  • Các nhà quản lý và cơ quan ban hành tiêu chuẩn: Tiếp cận bằng chứng khoa học định lượng để phê duyệt các thiết kế kết cấu mới sử dụng vật liệu FGM.

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ì? Đó là việc thiết lập thành công mô hình toán học biến phân nhất quán tích hợp lý thuyết biến dạng cắt bậc ba của Shi [86] (TSDT) với lý thuyết trường pha Phase-Field để phân tích ổn định và dao động của tấm FGM có chiều dày thay đổi tùy ý. Công trình đã chứng minh rằng TSDT của Shi mô tả chính xác ứng xử cắt uốn mà không cần hệ số hiệu chỉnh cắt nhân tạo, trong khi biến trường pha $s(x,y)$ biến đổi trơn từ 0 đến 1 đã chuyển hóa bài toán hình học nứt gián đoạn phức tạp thành bài toán giải tích miền liên tục trên lưới phần tử hữu hạn chuẩn.

  2. Đổi mới phương pháp luận của luận án so với các công trình quốc tế trước đây như thế nào? So với phương pháp XFEM của Moës et al. [19] và Liu et al. [17] (đòi hỏi định nghĩa các hàm làm giàu Heaviside và hàm đầu vết nứt phức tạp), hoặc phương pháp XIGA của Loc Tran et al. [68] (đòi hỏi xử lý cấu trúc lưới NURBS phi cấu trúc phức tạp), phương pháp Phase-Field kết hợp phần tử tam giác 6 nút của luận án cho phép xử lý số lượng vết nứt bất kỳ, góc nghiêng tùy ý và chiều dày biến thiên mà không cần theo dõi biên vết nứt hay chia lại lưới, đảm bảo tính ổn định số và tốc độ hội tụ vượt trội.

  3. Phát hiện nào trong luận án gây bất ngờ nhất về mặt cơ học kết cấu? Phát hiện bất ngờ nhất là hiện tượng: Tấm FGM có chiều dày thay đổi và bị nứt nặng ($c/L = 0,8$) khi đặt trên nền đàn hồi Pasternak ($k_w=50, k_s=10$) vẫn sở hữu tần số dao động tự do cao hơn đáng kể so với tấm nguyên vẹn không đặt trên nền. Đồng thời, khi góc nghiêng vết nứt $\alpha$ tăng từ $0^\circ$ lên $60^\circ$, lực mất ổn định tới hạn của tấm nứt không những không giảm mà còn tăng phục hồi gần 43%, đảo ngược hoàn toàn suy nghĩ thông thường rằng vết nứt càng dài thì kết cấu càng suy yếu tuyến tính.

  4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả nghiên cứu (Replication Protocol) không? Hoàn toàn đầy đủ và minh bạch. Luận án công bố chi tiết: hệ phương trình vi phân chủ đạo (2.1) - (2.30), ma trận hàm dạng phần tử tam giác 6 nút $[N]$, các ma trận vi phân $[B_1], [B_2], [B_3]$, thông số vật liệu phụ thuộc nhiệt độ phi tuyến (Bảng 1.2), kích thước hình học, tham số trường pha tối ưu ($H/\epsilon = 400, B = 10^5$), cùng sơ đồ thuật toán dòng chảy (Hình 2.3 và 2.4) trên môi trường mã nguồn mở FreeFem++-cs.

  5. Lộ trình nghiên cứu 10 năm được phác thảo cụ thể như thế nào? Lộ trình 10 năm mở rộng từ bài toán tấm tĩnh/dao động tuyến tính sang: (1) Động lực học phi tuyến và phá hủy động tấm/vỏ FGM; (2) Tích hợp hiệu ứng tiếp xúc phi tuyến hai bờ nứt; (3) Mở rộng cho composite áp điện, nano-FGM và FGM xốp 3D; (4) Chế tạo mẫu thực nghiệm FGM bằng in 3D kim loại - gốm và kiểm chứng bằng công nghệ đo trường biến dạng DIC tại nhiệt độ cực hạn (>1000K).

Kết luận

  1. Luận án đã giải quyết trọn vẹn và xuất sắc bài toán phân tích ảnh hưởng của vết nứt tới ổn định và dao động tự do của kết cấu tấm FGM có chiều dày không đổi và thay đổi (tuyến tính, parabol), đặt trên nền đàn hồi Pasternak và làm việc trong môi trường nhiệt độ cao.
  2. Tiên phong ứng dụng thành công lý thuyết trường pha Phase-Field trong cơ học phá hủy kết hợp với lý thuyết biến dạng cắt bậc ba của Shi [86] trên nền tảng phần tử hữu hạn 6 nút tham số, giải quyết triệt để hiện tượng nghẽn cắt mà không cần hệ số hiệu chỉnh.
  3. Đã làm sáng tỏ các quy luật cơ học chi phối: Sự phục hồi tải trọng mất ổn định khi tăng góc nghiêng vết nứt $\alpha$; sự suy giảm phi tuyến của tần số dao động theo chỉ số thể tích $n \in [0, 2]$; vai trò tăng cứng áp đảo của hệ số kháng cắt nền Pasternak $k_s$; và tác động sụt giảm tần số nghiêm trọng khi nhiệt độ mặt gốm $T_c$ gia tăng.
  4. Kiểm chứng tính đúng đắn và độ tin cậy vượt bậc thông qua đối chuẩn với dữ liệu thực nghiệm quốc tế của Seifi et al. [80] (sai số tối thiểu 0,39%) và các nghiệm giải tích 3D của Huang et al. [54] và Fazzolari [56].
  5. Mở ra 3 hướng nghiên cứu học thuật mới đầy triển vọng: Cơ học phá hủy trường pha cho kết cấu vỏ FGM mỏng/dày; phân tích động lực học phi tuyến và nứt động trong môi trường nhiệt - cơ - điện kết hợp; và tối ưu hóa hình học kết cấu FGM chịu hư hại.
  6. Toàn bộ các kết quả cốt lõi của luận án đã được công bố trên 9 công trình khoa học chất lượng cao (5 bài SCIE/SCI Q1, 2 bài tạp chí trong nước, 2 bài kỷ yếu hội nghị toàn quốc), khẳng định vị thế học thuật tiên phong và đóng góp bền vững cho kho tàng cơ học công trình tại Việt Nam và trên thế giới.