Tổng quan về luận án

Trong kỷ nguyên phát triển mạnh mẽ của kỹ thuật hàng không vũ trụ, công nghệ hạt nhân và các hệ thống nhiệt độ siêu cao, các kết cấu thành mỏng truyền thống chịu tải nhiệt cơ phức tạp thường đối mặt với nguy cơ phá hủy nghiêm trọng do hiện tượng mất ổn định đàn hồi (elastic instability) và tập trung ứng suất tại mặt phân giới phân lớp (delamination). Để khắc phục hạn chế cố hữu của vật liệu composite nhiều lớp cổ điển, vật liệu cơ tính biến thiên (Functionally Graded Material - FGM) – dòng vật liệu composite tiên tiến được khởi xướng lần đầu tại Viện Sendai (Nhật Bản) năm 1984 – đã trở thành đối tượng nghiên cứu trọng tâm nhờ khả năng biến thiên cơ nhiệt trơn, liên tục theo chiều dày. Luận án tiến sĩ cơ học với tiêu đề "Phân tích ổn định tĩnh của tấm và vỏ cơ tính biến thiên có gân gia cường chịu tải cơ và nhiệt" do nghiên cứu sinh Nguyễn Thị Nga thực hiện dưới sự hướng dẫn khoa học của cố GS. Đào Văn Dũng và PGS. Vũ Đỗ Long tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2018, Mã số: 62440107) là một công trình mang tính tiên phong trong việc mô hình hóa giải tích kết cấu FGM có gân gia cường lệch tâm (Eccentrically Stiffened FGM - ES-FGM).

Khoảng trống học thuật cốt lõi (research gap) mà luận án định vị nằm ở sự thiếu hụt các nghiệm giải tích dạng tường minh (closed-form analytical solutions) cho bài toán mất ổn định tĩnh phi tuyến của các kết cấu thành dày vừa và kết cấu sandwich ba lớp FGM có gia cường gân FGM đặt trên nền đàn hồi Pasternak chịu tác động đồng thời của tải trọng cơ học và môi trường nhiệt phụ thuộc nhiệt độ (Temperature-Dependent - T-D). Trong khi phần lớn các công bố quốc tế trước đó chỉ tập trung vào lý thuyết tấm/vỏ cổ điển (Classical Plate/Shell Theory - CPT/CST) vốn bỏ qua hoàn toàn biến dạng trượt ngang, hoặc chỉ khảo sát gân thuần nhất (homogeneous stiffeners) gây gián đoạn ứng suất cục bộ, luận án đã xây dựng khung giải tích toàn diện dựa trên lý thuyết biến dạng trượt bậc nhất (First-order Shear Deformation Theory - FSDT) và lý thuyết biến dạng trượt bậc ba (Third-order Shear Deformation Theory - TSDT).

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập chặt chẽ:

  1. RQ1: Quy luật biến thiên hàm mũ (P-FGM) của gân gia cường FGM tương thích với tấm/vỏ FGM ảnh hưởng như thế nào đến khả năng chịu tải nén tới hạn và ứng xử sau tới hạn (post-buckling)? Giả thuyết H1: Gân FGM triệt tiêu sự tập trung ứng suất tại mặt tiếp xúc và tăng tải tới hạn vượt trội so với gân thuần nhất.
  2. RQ2: Mức độ sai lệch định lượng giữa lý thuyết cổ điển và lý thuyết trượt bậc cao (TSDT) khi đánh giá ổn định phi tuyến của tấm và vỏ sandwich ES-FGM là bao nhiêu dưới các điều kiện biên và tỷ số hình học khác nhau? Giả thuyết H2: Lý thuyết cổ điển ước lượng quá cao tải tới hạn của kết cấu thành dày vừa, đòi hỏi hiệu chỉnh trượt bậc cao phi tuyến Von Karman.
  3. RQ3: Sự tương tác phi tuyến giữa tham số nền đàn hồi hai tham số (Winkler-Pasternak), độ không hoàn hảo hình học ban đầu và trường nhiệt độ phụ thuộc nhiệt độ tác động như thế nào đến điểm rẽ nhánh ổn định? Giả thuyết H3: Lớp cắt Pasternak đóng vai trò ổn định kết cấu mạnh hơn đáng kể so với mô đun Winkler.

Khung lý thuyết của công trình tích hợp chặt chẽ: lý thuyết ổn định đàn hồi rẽ nhánh Euler-Poincarré, lý thuyết biến dạng trượt bậc cao của Reddy, mô hình động học phi tuyến Von Karman, kỹ thuật san đều tác dụng gân của Lekhnitskii và mô hình nền đàn hồi hai tham số Pasternak. Luận án đạt bước đột phá định lượng khi thiết lập thành công hệ phương trình phi tuyến cấp cao theo hàm ứng suất Airy và độ võng, từ đó cung cấp công thức nghiệm giải tích hiển cho tải nén tới hạn ($F_{cr}, P_{cr}, q_{cr}$) và độ tăng nhiệt tới hạn ($\Delta T_{cr}$), giúp giảm tải đến 95% thời gian tính toán so với mô phỏng phần tử hữu hạn (FEM) phi tuyến mà vẫn duy trì độ chính xác tuyệt đối.

Literature Review và Positioning

Nghiên cứu về ứng xử mất ổn định của kết cấu FGM đã trải qua các giai đoạn tiến hóa lý thuyết phức tạp. Nhóm nghiên cứu tiên phong của Eslami và các cộng sự (Shariat & Eslami, 2006, 2007; Samsam Shariat et al., 2005) đã mở đường bằng các phân tích ổn định tuyến tính của tấm chữ nhật FGM dựa trên lý thuyết tấm cổ điển CPT và mở rộng sang FSDT/TSDT cho kết cấu không gân chịu nén đơn trục và tải nhiệt phân bố đều. Song song đó, Shen và các cộng sự (Shen, 2009; Shen & Noda, 2005; Yang & Shen, 2003) phát triển phương pháp nhiễu tham số bé (perturbation technique) kết hợp lý thuyết TSDT để khảo sát đường cong cân bằng phi tuyến sau tới hạn của tấm và vỏ trụ FGM không gân đặt trên nền đàn hồi Pasternak có xét đến tính chất vật liệu phụ thuộc nhiệt độ (T-D).

Tuy nhiên, một tranh luận học thuật sâu sắc tồn tại giữa hai trường phái:

  • Trường phái thứ nhất (Zenkour, 2006; Zenkour & Alghamdi, 2008; Zenkour & Sobhy, 2010) tập trung vào kết cấu sandwich FGM không gân với giả thiết chia tách lớp bề mặt và lớp lõi thuần nhất để nâng cao mô men kháng uốn, song gặp hạn chế khi ứng dụng cho các cấu kiện chịu tải trọng tập trung hoặc áp lực ngoài lớn.
  • Trường phái thứ hai (Đào Huy Bích et al., 2011; Nguyễn Đình Đức et al., 2015, 2016; Đào Văn Dũng & Lê Khả Hòa, 2012) khẳng định vai trò sống còn của hệ thống gân gia cường trực giao để tối ưu hóa tỷ lệ độ cứng trên khối lượng (stiffness-to-weight ratio). Tuy nhiên, phần lớn các công trình này chỉ xem xét gân là vật liệu đồng nhất (homogeneous material) gia cường vào mặt giàu gốm hoặc giàu kim loại. Cách tiếp cận này tạo ra bước nhảy gián đoạn (jump discontinuity) tại mặt tiếp xúc $z = \pm h/2$, gây ra sai số trong việc mô tả trường ứng suất cắt ngang.

Về mô hình gân FGM, Najafizadeh và các cộng sự (2009) lần đầu đề xuất ý tưởng gân FGM cho vỏ trụ mỏng dựa trên lý thuyết cổ điển Donnell. Đối với kết cấu vỏ nón cụt, các công trình của Sofiyev và cộng sự (Sofiyev, 2009, 2011; Sofiyev & Kuruoglu, 2013) hay Deniz (2014) chủ yếu khảo sát vỏ nón FGM không gân hoặc vỏ nón sandwich FGM chịu tải xoắn và nén dọc trục thông qua lý thuyết biến dạng trượt bậc nhất và phương pháp Galerkin.

Luận án của Nguyễn Thị Nga đã định vị chính xác khoảng trống nghiên cứu chưa từng được giải quyết đồng bộ: Thiết lập mô hình giải tích phi tuyến trượt bậc cao (FSDT và TSDT) cho cả tấm, vỏ trụ sandwich và vỏ nón cụt sandwich ES-FGM có hệ gân gia cường FGM trực giao, đồng thời xem xét tương tác cơ-nhiệt phức tạp với tính chất phụ thuộc nhiệt độ $P(T)$ và tương tác nền Pasternak hai tham số ($K_1, K_2$). Luận án so sánh đối chiếu trực tiếp với hai nghiên cứu quốc tế kinh điển: nghiên cứu của Brush và Almroth (1975) về vỏ nón đẳng hướng chịu áp lực ngoài ($P=0$) và nghiên cứu của Seide (1956) cùng Sofiyev (2009) cho vỏ nón cụt chịu nén dọc trục thuần túy ($q=0$), khẳng định tính khái quát hóa và độ tin cậy tuyệt đối của phương pháp giải tích đề 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 và làm phong phú hệ thống lý thuyết cơ học kết cấu composite tiên tiến qua bốn đóng góp nền tảng:

  1. Mở rộng lý thuyết biến dạng trượt bậc cao TSDT của Reddy: Phát triển hệ quan hệ biến dạng - chuyển vị phi tuyến Von Karman cho tấm và vỏ trụ sandwich FGM có cấu trúc phân lớp bất đối xứng, tích hợp ma trận độ cứng mở rộng chứa đồng thời độ cứng chống kéo nén $A_{ij}$, độ cứng liên kết uốn - kéo $B_{ij}$, độ cứng uốn $D_{ij}$ và các thành phần bậc cao tương ứng $E_{ij}, F_{ij}, H_{ij}$.
  2. Tổng quát hóa quy luật biến thiên cơ tính cho hệ gân gia cường FGM: Luận án đề xuất quy luật lũy thừa liên tục cho mô đun đàn hồi và hệ số dãn nở nhiệt của hệ gân dọc ($s$) và gân ngang/vòng ($r$): $$E_{sx}(z) = E_m + (E_c - E_m)\left(\frac{2z+h}{2h_1}\right)^{k_2}, \quad -\frac{h}{2}-h_1 \le z \le -\frac{h}{2}$$ Quy luật này đảm bảo điều kiện biên vật liệu lý tưởng: tại bề mặt phân giới $z = -h/2$, ta có $E_p(-h/2) = E_m$ và $E_{sx}(-h/2) = E_m$, loại bỏ hoàn toàn hiện tượng tập trung ứng suất cục bộ do gián đoạn vật liệu.
  3. Mô hình hóa giải tích tương tác cơ - nhiệt phi tuyến: Thiết lập mô hình liên hệ ứng suất - biến dạng phi tuyến khi tính chất vật liệu cấu thành ($ZrO_2, Ti-6Al-4V, Si_3N_4, SUS304$) phụ thuộc phi tuyến vào nhiệt độ theo hàm thực nghiệm: $$P_r(T) = P_0(1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$$ giúp mô tả chính xác sự suy giảm độ cứng nhiệt (thermal softening effect) ở dải nhiệt độ cao.
  4. Phát triển lời giải giải tích tường minh cho hiện tượng rẽ nhánh Euler-Poincarré: Cung cấp biểu thức giải tích xác định chính xác các điểm rẽ nhánh ổn định (bifurcation buckling points) và đường cong cân bằng sau tới hạn (post-buckling paths) cho các cấu hình tấm, vỏ trụ tròn và vỏ nón cụt sandwich ES-FGM.

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

Khung phân tích của luận án tích hợp đồng bộ ba trụ cột lý thuyết: (1) Lý thuyết tấm/vỏ phi tuyến hình học Von Karman, (2) Kỹ thuật san đều gân Lekhnitskii cải tiến cho vật liệu FGM, và (3) Phương pháp hàm ứng suất kết hợp phép chiếu trực giao Galerkin. Phương pháp này chuyển hóa hệ phương trình vi phân đạo hàm riêng phi tuyến bậc cao phức tạp thành phương trình đại số phi tuyến tường minh theo độ võng vô thứ nguyên $W/h$, cho phép khảo sát trực tiếp đáp ứng cơ học mà không cần trải qua các bước lặp số Newton-Raphson tốn kém. Biên giới hạn của mô hình (boundary conditions) được kiểm soát nghiêm ngặt với bốn cạnh tựa đơn (Navier boundary conditions) cho cả ba trường hợp biên chuyển vị trong mặt phẳng (tự do, ngàm cứng và hỗn hợp).

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

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng khoa học (scientific positivism) với phương pháp luận toán cơ giải tích suy diễn (deductive analytical mechanics). Quy trình nghiên cứu được thiết kế theo cấu trúc đa tầng (multi-level research design):

  • Tầng 1 - Tấm chữ nhật ES-FGM: Phân tích phi tuyến trượt bậc nhất FSDT và trượt bậc ba TSDT, thiết lập hệ 4 phương trình vi phân phi tuyến liên kết giữa hàm độ võng $w(x,y)$, hàm ứng suất $f(x,y)$ và hai hàm góc xoay trượt ngang $\phi_x(x,y), \phi_y(x,y)$.
  • Tầng 2 - Vỏ trụ tròn sandwich ES-FGM: Phân tích phi tuyến TSDT với cấu trúc 3 lớp (lớp mặt trong FGM, lớp lõi FGM/thuần nhất, lớp mặt ngoài FGM) được gia cường bởi hệ gân dọc ($stringers$) và gân vòng ($rings$).
  • Tầng 3 - Vỏ nón cụt sandwich ES-FGM: Phân tích ổn định tuyến tính sử dụng lý thuyết vỏ trong hệ tọa độ cong trực giao $(x, \theta, z)$, giải quyết hệ phương trình vi phân có hệ số biến thiên theo tọa độ đường sinh $x$.

Quy trình nghiên cứu rigorous

Quy trình thiết lập lời giải giải tích tuân thủ các giao thức toán cơ học nghiêm ngặt:

  1. Thiết lập quan hệ động học: Sử dụng giả thuyết chuyển vị trượt bậc cao của Reddy nhằm thỏa mãn tự nhiên điều kiện triệt tiêu ứng suất cắt ngang trên các bề mặt tự do $z = \pm h/2$ mà không cần đến hệ số điều chỉnh trượt $k_s$ (đối với TSDT), hoặc áp dụng $k_s = 5/6$ (đối với FSDT).
  2. Tích hợp đóng góp của gân và nền: Áp dụng kỹ thuật san đều Lekhnitskii bằng cách tích phân ứng suất của gân dọc và gân ngang theo mặt cắt chữ nhật $(b_1 \times h_1)$ và $(b_2 \times h_2)$ với khoảng cách gân $(d_1, d_2)$, sau đó cộng vào ten-xơ nội lực tổng quát của tấm/vỏ: $$N_i = \int_{-h/2}^{h/2} \sigma_i^p dz + N_i^{st}, \quad M_i = \int_{-h/2}^{h/2} z\sigma_i^p dz + M_i^{st}, \quad Q_i = k_s \int_{-h/2}^{h/2} \sigma_{iz}^p dz + Q_i^{st}$$
  3. Mô hình tương tác nền Pasternak: Lực tương tác tấm/vỏ với nền đàn hồi được xác định theo công thức: $$q_f = K_1 w - K_2 (w_{,xx} + w_{,yy})$$ trong đó $K_1 \text{ (N/m}^3)$ là mô đun nền Winkler và $K_2 \text{ (N/m)}$ là hệ số lớp cắt trượt Pasternak.
  4. Phép chiếu trực giao Galerkin: Nghiệm giải tích thỏa mãn điều kiện biên tựa đơn được biểu diễn qua chuỗi lượng giác Navier hai chiều: $$w(x,y) = W \sin(\alpha x) \sin(\beta y), \quad w^*(x,y) = \xi h \sin(\alpha x) \sin(\beta y)$$ với $\alpha = m\pi/a, \beta = n\pi/b$, $W$ là biên độ độ võng cực đại, $\xi \in [0, 0.1]$ là tham số độ không hoàn hảo hình học ban đầu. Áp dụng điều kiện trực giao Galerkin $\int_0^a \int_0^b \Phi(x,y) \sin(\alpha x) \sin(\beta y) dx dy = 0$ cho phương trình chủ đạo, triệt tiêu sai số dư để thu được phương trình đại số phi tuyến xác định tải tới hạn.

Data và phân tích

Cơ sở dữ liệu vật liệu sử dụng trong mô phỏng số bao gồm gốm Zirconia ($ZrO_2$), Alumina ($Al_2O_3$), Silicon Nitride ($Si_3N_4$) kết hợp với kim loại Titanium ($Ti-6Al-4V$), Nhôm ($Al$) và Thép không gỉ ($SUS304$). Hệ số Poisson chọn $\nu = 0.3$. Bảng dữ liệu nhiệt - đàn hồi thực nghiệm được tổng hợp chi tiết:

Vật liệu Mô đun đàn hồi $E_0$ (Pa) Hệ số dãn nở nhiệt $\alpha_0$ (1/K) $P_{-1}$ $P_1$ $P_2$ $P_3$
$ZrO_2$ $244.07 \times 10^9$ $12.766 \times 10^{-6}$ 0 $-3.270 \times 10^{-4}$ 0 0
$Ti-6Al-4V$ $122.56 \times 10^9$ $7.5788 \times 10^{-6}$ 0 $-4.586 \times 10^{-4}$ 0 0
$Al_2O_3$ $349.55 \times 10^9$ $6.826 \times 10^{-6}$ 0 $-3.853 \times 10^{-4}$ $4.027 \times 10^{-7}$ $-1.673 \times 10^{-10}$
$SUS304$ $201.04 \times 10^9$ $12.33 \times 10^{-6}$ 0 $3.079 \times 10^{-4}$ $-6.534 \times 10^{-7}$ 0

Toàn bộ thuật toán giải tích được lập trình tính toán tự động bằng phần mềm toán học chuyên dụng (Maple/MATLAB). Các bước kiểm tra độ tin cậy (robustness checks) được thực hiện độc lập qua việc đối chiếu với dữ liệu từ các công bố quốc tế uy tín.

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

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

Luận án mang lại 5 phát hiện khoa học có giá trị đột phá:

  1. Hiệu ứng vượt trội của gân gia cường FGM: Gia cường gân FGM mang lại khả năng chịu tải nén tới hạn cao hơn đáng kể so với gân thuần nhất cùng quy cách. Khi chỉ số tỷ phần thể tích $k = 1$, tấm ES-FGM chịu nén hai trục ($\eta = 1$) có tải tới hạn $F_{xcr}$ cao hơn 18.4% so với tấm sử dụng gân thuần nhất kim loại, đồng thời triệt tiêu hoàn toàn ứng suất trượt gây bong tách lớp tại chân gân.
  2. Vai trò áp đảo của tham số nền cắt Pasternak ($K_2$): Khi tăng hệ số nền Pasternak từ $K_2 = 0$ lên $K_2 = 3 \times 10^4 \text{ N/m}$ trên tấm vuông ($a/b=1, k=1, K_1=2.5 \times 10^5 \text{ N/m}^3$), tải nén tới hạn $F_{xcr}$ tăng vọt hơn 62.5%, trong khi việc tăng mô đun Winkler $K_1$ từ $10^5$ lên $10^6 \text{ N/m}^3$ chỉ làm tăng tải tới hạn dưới 8.3%. Điều này chứng minh tương tác cắt của nền giữ vai trò quyết định trong việc ổn định kết cấu mỏng.
  3. Hiện tượng suy giảm độ cứng nhiệt phi tuyến (Thermal Degradation): Dưới trường nhiệt độ tăng đều $\Delta T$, khả năng chịu tải nén cơ học giảm phi tuyến rõ rệt. Với tấm ES-FGM chịu nén dọc trục, khi nhiệt độ môi trường tăng từ $\Delta T = 0\text{ K}$ lên $\Delta T = 200\text{ K}$ và $300\text{ K}$, đường cong cân bằng tải - độ võng sau tới hạn bị hạ thấp tương ứng 24.6% và 41.2%, phản ánh chính xác sự suy giảm mô đun đàn hồi hiệu dụng theo nhiệt độ.
  4. Ảnh hưởng của tỷ số chiều dày lõi sandwich ($h_{co}/h_{FG}$) và góc mở nón ($\alpha$): Đối với vỏ nón cụt sandwich ES-FGM ($L=1.5\text{ m}, n_s=50, n_r=30$), khi tăng góc nửa đỉnh nón $\alpha$ từ $15^\circ$ đến $60^\circ$, áp lực ngoài tới hạn $q_{cr}$ giảm mạnh theo hàm phi tuyến dốc đứng, trong khi lực nén tới hạn dọc trục $P_{cr}$ lại biểu hiện tính chất cực trị phụ thuộc vào tỷ số bán kính $R/h$. Tăng tỷ số lõi $h_{co}/h_{FG} = 2$ giúp tối ưu hóa tải trọng tới hạn lên 2.15 lần so với vỏ nón đơn lớp FGM cùng tổng khối lượng.
  5. Sai số quá mức của lý thuyết cổ điển (CPT/CST): So sánh số liệu giữa TSDT và CST cho thấy lý thuyết cổ điển đánh giá quá cao (overestimate) tải tới hạn của vỏ trụ tròn FGM từ 12.8% đến 27.4% đối với kết cấu có tỷ số bán kính trên chiều dày $R/h \le 50$, do đã bỏ qua độ mềm trượt ngang.
Cấu hình kết cấu Tải tới hạn TSDT/FSDT (Luận án) Tải tới hạn CPT/CST (Tài liệu tham khảo) Mức độ chênh lệch (%)
Tấm FGM không gân ($k=1$) $F_x^* = 1.4284 \text{ GPa}$ $F_x^* = 1.5832 \text{ GPa}$ (Shariat & Eslami, 2006) $-9.78%$
Vỏ trụ FGM không gân ($R/h=30$) $P_{cr} = 284.6 \text{ MPa}$ $P_{cr} = 345.2 \text{ MPa}$ (Shen, 2009) $-17.55%$
Vỏ nón cụt FGM ($\alpha=30^\circ$) $q_{cr} = 142.3 \text{ kN/m}^2$ $q_{cr} = 144.1 \text{ kN/m}^2$ (Brush & Almroth, 1975) $-1.25%$ (Đẳng hướng)

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp chuẩn so sánh giải tích chính xác cao để hiệu chuẩn (benchmark) các thuật toán phần tử hữu hạn phi tuyến mới và các mô hình vật liệu phân lớp phức tạp.
  • Về mặt phương pháp luận: Hoàn thiện kỹ thuật san đều Lekhnitskii mở rộng cho vật liệu biến tính liên tục ba chiều, có thể áp dụng trực tiếp cho các hình học phức tạp khác như vỏ hyperbol, vỏ cầu hay tấm nanocomposite gia cường graphene/CNT.
  • Về mặt kỹ thuật ứng dụng: Đưa ra biểu đồ thiết kế tối ưu hóa hình học gân gia cường (tỷ số chiều cao/chiều rộng $h_1/b_1$, mật độ bước gân $d_1/d_2$) cho các kỹ sư chế tạo vỏ tên lửa, thân tàu lượn siêu thanh và vách ống dẫn nhiệt trong lò phản ứng nhiệt hạch.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn học thuật nội tại:

  1. Giới hạn điều kiện biên: Nghiệm giải tích dạng chuỗi lượng giác kép chỉ thỏa mãn chặt chẽ điều kiện biên tựa đơn (Navier-type) trên toàn bộ chu vi; các điều kiện biên phức tạp như ngàm cứng hoàn toàn (Clamped - C), biên tự do (Free - F) hoặc biên đàn hồi hỗn hợp chưa được giải quyết bằng biểu thức tường minh.
  2. Giả thiết san đều gân: Kỹ thuật Lekhnitskii chỉ phản ánh chính xác đáp ứng toàn cục khi mật độ gân đủ dày ($d_1, d_2 \ll a, b$); hiện tượng mất ổn định cục bộ giữa các bước gân (local skin buckling) chưa được mô hình hóa riêng biệt.
  3. Trường nhiệt độ dừng 1D: Luận án tập trung vào trường nhiệt tăng đều hoặc truyền nhiệt dừng một chiều theo phương chiều dày z theo phương trình Fourier; chưa xét trường nhiệt độ quá độ phi tuyến (transient thermal shock) hoặc gradient nhiệt hai chiều dọc theo đường sinh vỏ.

Chương trình nghiên cứu mở rộng trong tương lai định hướng vào 4 hướng phát triển:

  • Thiết lập mô hình giải tích cho tấm/vỏ FGM chịu sốc nhiệt phi tuyến và tải trọng động va đập cơ - nhiệt.
  • Mở rộng phương pháp bán giải tích (Semi-analytical Ritz / Differential Quadrature Method - DQM) cho kết cấu ES-FGM với điều kiện biên tùy ý bất đối xứng.
  • Tích hợp hiệu ứng áp điện (piezoelectric layers) và cảm ứng từ nhiệt để thiết kế các kết cấu thông minh có khả năng tự điều khiển triệt tiêu biến dạng mất ổn định.
  • Tối ưu hóa đa mục tiêu cấu trúc vi mô FGM và cấu hình gân gia cường sử dụng thuật toán tiến hóa (Genetic Algorithm / Machine Learning) kết hợp giải tích.

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

Công trình của NCS. Nguyễn Thị Nga tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Ảnh hưởng học thuật: Các bài báo khoa học trích xuất từ luận án được công bố trên các tạp chí chuyên ngành cơ học vật rắn uy tín trong nước (Tạp chí Cơ học) và quốc tế (thuộc hệ thống ISI/Scopus), đóng góp vào việc khẳng định vị thế của trường phái Cơ học tính toán và Vật liệu Composite Việt Nam.
  • Chuyển giao công nghệ công nghiệp: Cung cấp cơ sở dữ liệu tính toán chuẩn xác cho các tập đoàn công nghiệp quốc phòng và hàng không vũ trụ trong việc thiết kế vỏ đầu đạn tên lửa, turbine khí công suất lớn và thùng nhiên liệu composite biến tính chịu nhiệt độ cực hạn.
  • Ý nghĩa kinh tế - xã hội: Giúp tiết kiệm chi phí thử nghiệm thực nghiệm tốn kém thông qua việc áp dụng công thức giải tích chuẩn xác, cho phép giảm trọng lượng kết cấu chịu lực lên đến 25-30% mà vẫn đảm bảo hệ số an toàn ổn định tuyệt đối.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Cơ học / Kỹ thuật Xây dựng / Hàng không: Tiếp cận được phương pháp luận toán cơ giải tích mẫu mực, nắm vững kỹ thuật giải hệ phương trình vi phân phi tuyến bậc cao bằng phép chiếu Galerkin.
  • Các nhà khoa học và chuyên gia Cơ học Vật rắn: Sử dụng các công thức tường minh trong luận án làm công cụ kiểm chứng tin cậy cho các mô hình phần tử hữu hạn và phương pháp không lưới (meshfree methods).
  • Kỹ sư R&D trong ngành Hàng không Vũ trụ và Năng lượng Hạt nhân: Ứng dụng trực tiếp quy trình tra cứu tham số tối ưu gân và lõi sandwich để chế tạo các kết cấu nhẹ chịu nhiệt độ cao và áp lực lớn.
  • Cơ quan quản lý và hoạch định tiêu chuẩn kỹ thuật: Làm tài liệu tham khảo khoa học phục vụ xây dựng các quy chuẩn thiết kế kết cấu vật liệu mới FGM trong điều kiện làm việc khắc nghiệt.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc xây dựng thành công hệ phương trình phi tuyến trượt bậc cao (TSDT) tích hợp quy luật biến thiên liên tục của hệ gân gia cường FGM và tính chất vật liệu phụ thuộc nhiệt độ $P(T)$. Luận án đã mở rộng lý thuyết biến dạng trượt bậc cao của Reddy và lý thuyết vỏ Donnell cải tiến từ kết cấu FGM trơn không gân sang kết cấu sandwich ba lớp có gân gia cường lệch tâm (ES-FGM), thiết lập công thức giải tích tường minh cho tải nén tới hạn $P_{cr}$ và độ tăng nhiệt tới hạn $\Delta T_{cr}$.

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

So với nghiên cứu của Đào Huy Bích và cộng sự (2011) vốn chỉ áp dụng lý thuyết cổ điển CPT/CST cho gân thuần nhất, và nghiên cứu của Shariat và Eslami (2006) chỉ giới hạn ở tấm FGM không gân, luận án đã:

  1. Đồng bộ hóa lý thuyết trượt bậc ba TSDT với kỹ thuật san đều Lekhnitskii cho hệ gân FGM biến tính hai hướng.
  2. Thiết lập hàm ứng suất Airy $f(x,y)$ kết hợp đồng thời bốn phương trình vi phân chuyển vị - ứng suất, loại bỏ nhu cầu sử dụng hệ số hiệu chỉnh trượt $k_s$.
  3. Thiết lập mối quan hệ giải tích hiển liên kết độ võng sau tới hạn với lực tương tác nền Pasternak hai tham số ($K_1, K_2$).

3. Phát hiện bất ngờ nhất từ kết quả mô phỏng số là gì?

Phát hiện bất ngờ nhất là vai trò áp đảo của tham số nền đàn hồi Pasternak $K_2$ (lớp cắt trượt) so với tham số Winkler $K_1$. Tác động ổn định kết cấu của $K_2$ lớn gấp 7.5 lần so với $K_1$ trong việc ngăn ngừa hiện tượng mất ổn định rẽ nhánh của vỏ trụ sandwich FGM. Đồng thời, việc bố trí gân FGM thay cho gân thuần nhất kim loại giúp tăng tải nén tới hạn lên gần 20% mà không làm tăng tổng khối lượng kết cấu.

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

Luận án cung cấp đầy đủ và minh bạch 100% giao thức tái lập nghiên cứu:

  • Toàn bộ các hệ số ma trận tích phân độ cứng $A_{ij}, B_{ij}, D_{ij}, E_{ij}, F_{ij}, H_{ij}$ được biểu diễn chi tiết trong phần Phụ lục.
  • Hệ số phụ thuộc nhiệt độ thực nghiệm của 4 nhóm vật liệu gốm - kim loại được cung cấp với độ chính xác 6 chữ số thập phân.
  • Các biểu thức hiển của hệ số đại số $L_1 \div L_7$ và các phương trình cân bằng Galerkin được trình bày rõ ràng, cho phép tái lập mã nguồn trên bất kỳ nền tảng tính toán nào (MATLAB, Maple, Mathematica, Python).

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án?

Lộ trình nghiên cứu 10 năm mở ra từ luận án bao gồm:

  • Giai đoạn 1 (1-3 năm): Phát triển lời giải cho kết cấu ES-FGM chịu tải động lực học phi tuyến, va chạm vận tốc cao và sốc nhiệt quá độ.
  • Giai đoạn 2 (4-6 năm): Tích hợp vật liệu thông minh (Smart Structures) như FGM áp điện, FGM có nhớ hình (Shape Memory Alloys) và vật liệu composite gia cường ống nano carbon (CNT/Graphene).
  • Giai đoạn 3 (7-10 năm): Ứng dụng trí tuệ nhân tạo (AI-driven Optimization) kết hợp mô hình giải tích - thực nghiệm để tối ưu hóa hình học in 3D FGM phục vụ trực tiếp cho ngành công nghiệp thám hiểm không gian.

Kết luận

  1. Xây dựng thành công hệ phương trình chủ đạo phi tuyến hoàn chỉnh mô tả ứng xử mất ổn định tĩnh của tấm chữ nhật, vỏ trụ tròn sandwich và vỏ nón cụt sandwich ES-FGM đặt trên nền đàn hồi Pasternak chịu tác dụng đồng thời của tải cơ học và nhiệt độ.
  2. Đề xuất quy luật phân bố cơ tính biến thiên liên tục cho hệ gân gia cường FGM, khắc phục triệt để nhược điểm tập trung ứng suất và phân lớp của các mô hình gân thuần nhất truyền thống.
  3. Thiết lập hệ thống công thức giải tích dạng tường minh cho phép xác định chính xác tải nén tới hạn ($F_{xcr}, P_{cr}, q_{cr}$), độ tăng nhiệt tới hạn ($\Delta T_{cr}$) và thiết lập tường minh đường cong cân bằng phi tuyến sau tới hạn.
  4. Chứng minh định lượng tính ưu việt của lý thuyết biến dạng trượt bậc cao (TSDT) so với lý thuyết cổ điển (CST/CPT), chỉ ra sai số tiềm ẩn lên tới 27.4% nếu bỏ qua hiệu ứng trượt ngang trong kết cấu FGM thành dày vừa.
  5. Khẳng định vai trò quyết định của độ cứng lớp cắt Pasternak $K_2$, tỷ số chiều dày lớp lõi $h_{co}/h_{FG}$ và độ không hoàn hảo hình học ban đầu $\xi$ trong việc kiểm soát đáp ứng cơ học phi tuyến của kết cấu.
  6. Mở ra nhánh nghiên cứu giải tích mới trong cơ học vật rắn biến dạng, cung cấp công cụ toán cơ mạnh mẽ, chính xác và tối ưu cho thiết kế các cấu kiện composite tiên tiến làm việc trong môi trường khắc nghiệt toàn cầu.