Luận án TS Phạm Hồng Công: Phân tích phi tuyến tĩnh & động học tấm FGM trên nền đàn hồi

Tổng quan về luận án

Nghiên cứu của tác giả Phạm Hồng Công (2018) với tiêu đề "Phân tích phi tuyến tĩnh và động lực học của tấm chữ nhật FGM trên nền đàn hồi" (Chuyên ngành Cơ kỹ thuật, Mã số 62 52 01 01; Người hướng dẫn khoa học: GS. Nguyễn Đình Đức; Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội) là công trình nghiên cứu tiên phong trong việc thiết lập hệ thống phương pháp giải tích giải quyết toàn diện bài toán ổn định phi tuyến tĩnh và đáp ứng động lực học của kết cấu tấm vật liệu có cơ tính biến đổi (Functionally Graded Materials - FGM). Trong bối cảnh vật liệu composite chức năng FGM đang trở thành trụ cột của các ngành kỹ thuật công nghệ cao như hàng không vũ trụ, lò phản ứng hạt nhân, tàu ngầm và kết cấu chịu nhiệt độ cao, việc giải quyết các thách thức về tập trung ứng suất, hiện tượng delamination (bong tách lớp) và phá hủy giòn thông qua cấu trúc chuyển tiếp liên tục thành phần kim loại – ceramic là yêu cầu cấp thiết.

Khoảng trống nghiên cứu (research gap) mà luận án tập trung giải quyết xuất phát từ thực tiễn: hầu hết các công bố quốc tế trước thời điểm này chỉ tập trung vào tấm mỏng FGM không gân, hoặc sử dụng lý thuyết tấm cổ điển (Classical Plate Theory - CPT) mà bỏ qua hiệu ứng biến dạng trượt ngang đối với tấm dày; đồng thời các nghiên cứu về tấm FGM có gân gia cường lệch tâm (ES-FGM) hầu như chưa tính toán đồng thời sự phụ thuộc vào nhiệt độ của cả vật liệu tấm và gân (Temperature-Dependent properties - T-D), cũng như thiếu vắng các giải pháp giải tích tường minh cho tấm ES-FGM tích hợp lớp áp điện (Piezoelectric ES-FGM) chịu trường đa vật lý (cơ – nhiệt – điện).

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác định cụ thể:

• Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập mô hình giải tích chính xác mô tả ứng xử phi tuyến sau tới hạn của tấm mỏng và tấm dày ES-FGM trên nền đàn hồi hai thông số Pasternak khi kể đến sự biến thiên cơ tính theo nhiệt độ?
• Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố thành phần theo hàm Sigmoid (S-FGM) ảnh hưởng như thế nào đến tần số dao động tự do tuyến tính, quan hệ biên độ – tần số và đáp ứng thời gian – độ võng phi tuyến so với mô hình lũy thừa (P-FGM)?
• Câu hỏi nghiên cứu 3 (RQ3): Cơ chế ghép nối cơ – nhiệt – điện trong tấm dày ES-FGM áp điện kiểm soát hiện tượng phách điều hòa (beating phenomenon) và ổn định động học phi tuyến ra sao?
• Giả thuyết 1 (H1): Việc kết hợp kỹ thuật san đều tác dụng gân cải tiến của Lekhnitsky với lý thuyết biến dạng trượt bậc nhất (FSDT) và bậc ba (TSDT) sẽ cho phép xác định chính xác tải trọng tới hạn nhiệt – cơ mà không gặp phải hiện tượng khóa trượt (shear locking) của các mô hình số.
• Giả thuyết 2 (H2): Nền đàn hồi Pasternak ($k_1, k_2$) và điện áp điều khiển áp điện có khả năng triệt tiêu đáng kể biên độ dao động phi tuyến và dịch chuyển đường cong cân bằng sau tới hạn về phía an toàn hơn.

Khung lý thuyết của nghiên cứu được xây dựng trên sự tích hợp của cơ học môi trường liên tục phi tuyến, lý thuyết biến dạng lớn von Kármán, tiêu chuẩn ổn định rẽ nhánh Euler-Poincaré, lý thuyết tiếp xúc nền đàn hồi Pasternak hai tham số và cơ học vật liệu áp điện thông minh. Phạm vi khảo sát bao quát các cấu hình tấm chữ nhật mỏng ($a/h \ge 50$) đến tấm dày ($a/h = 10 \div 20$), kích thước mẫu chuẩn hóa $a \times b = 1.0\text{ m} \times 1.0\text{ m}$, chiều dày $h = 0.01\text{ m}$, khảo sát trên các cặp tổ hợp vật liệu siêu bền chịu nhiệt như $Si_3N_4/SUS304$, $ZrO_2/Ti-6Al-4V$, $Al_2O_3/Al$ và lớp áp điện $PZT-4/PZT-5A$.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy các dòng nghiên cứu chủ đạo về kết cấu FGM được phân chia thành ba trường phái chính:

1. Trường phái sử dụng phương pháp khai triển tiệm cận tham số bé kết hợp phép lặp vi phân bậc cao do Hui-Shen Shen và cộng sự (2009, 2012, 2014) dẫn đầu [101-119], chuyên sâu khảo sát vỏ trụ, panel trụ và tấm nanocomposite/áp điện chịu tải nhiệt – cơ phức tạp nhưng chủ yếu áp dụng cho kết cấu không gân hoặc tấm trơn.
2. Trường phái phân tích ổn định tĩnh và động của vỏ nón, vỏ trụ FGM dựa trên lý thuyết Donnell và biến dạng trượt của Sofiyev và Schnack (2010, 2014) [121-136], tập trung vào tải nén dọc trục và tải xoắn trên nền đàn hồi.
3. Trường phái nghiên cứu của Javaheri & Eslami (2002) [76-78], Shariat & Javaheri (2006) [92], Amabili và cộng sự (2008, 2011) [12-15] về ổn định nhiệt đàn hồi của tấm chữ nhật và vỏ hai độ cong sử dụng phương pháp hàm năng lượng và mô hình đa bậc tự do.

Trong nước, các nghiên cứu nền tảng của Nguyễn Đình Đức và Hoàng Văn Tùng (2010, 2011) [141-143], Đào Huy Bích và cộng sự (2011, 2013) [19, 25], Đào Văn Dũng và Nguyễn Thị Nga (2013, 2014) [56, 60-63] đã giải quyết nhiều bài toán quan trọng đối với vỏ và tấm FGM. Tuy nhiên, tồn tại hai luồng quan điểm và tranh luận học thuật sâu sắc:

• Tranh luận 1 (Mô hình hóa gân gia cường): Một bên (như Najafizadeh et al., 2009 [82]) giả thiết gân làm bằng vật liệu FGM liên tục biến đổi; trong khi phía Đào Huy Bích (2011) và Nguyễn Đình Đức (2013) [25, 38] chứng minh rằng gân FGM rất khó chế tạo trong thực tế, do đó đề xuất gân thuần nhất (gân ở mặt ceramic làm bằng ceramic, gân ở mặt kim loại làm bằng kim loại) là giải pháp tối ưu kỹ thuật, đồng thời đòi hỏi phải xây dựng công thức tính toán ứng suất nhiệt riêng biệt phát sinh trong gân: (\sigma_x^{st}, \sigma_y^{st}) = E_0 (\epsilon_x, \epsilon_y) - E_0 \alpha_0 \Delta T (1, 1).
• Tranh luận 2 (Ảnh hưởng của nhiệt độ lên cơ tính): Các nghiên cứu truyền thống thường đơn giản hóa cơ tính vật liệu không đổi theo nhiệt độ (T-ID) để giảm tải tính toán giải tích. Luận án đã phản biện quan điểm này bằng cách chứng minh rằng ở dải nhiệt độ làm việc cao ($\Delta T > 300\text{ K}$), việc bỏ qua tính chất phụ thuộc nhiệt độ (T-D) sẽ dẫn đến sai số dự báo tải tới hạn vượt quá $25% - 35%$, gây mất an toàn nghiêm trọng cho kết cấu.

Vị trí của luận án trong bản đồ học thuật được định vị là công trình đầu tiên xây dựng thành công lời giải giải tích đóng (closed-form solutions) cho hệ phương trình vi phân đạo hàm riêng phi tuyến bậc cao của tấm dày ES-FGM tích hợp áp điện trên nền Pasternak dựa trên cả hai lý thuyết biến dạng trượt bậc cao FSDT và TSDT. So với nghiên cứu của Javaheri & Eslami (2002) [77] vốn chỉ giới hạn ở tấm trơn không gân và nghiên cứu của Shen & Wang (2010) [117] sử dụng phương pháp số xấp xỉ, luận án đã mở rộng khung giải tích giải quyết đồng thời bài toán ổn định sau tới hạn tĩnh và đáp ứng dao động phi tuyến cưỡng bức có cản và có điều khiển điện áp chủ động.

Đó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 phi tuyến von Kármán và lý thuyết biến dạng trượt bậc ba của Reddy (TSDT) cho môi trường vật liệu chức năng dị tính liên tục qua chiều dày có hệ thống gân trực giao lệch tâm:

P_{eff}(z, T) = P_{rc}(T) V_c(z) + P_{rm}(T) V_m(z)

trong đó cơ tính thành phần $P_r(T)$ tuân theo quy luật phi tuyến nhiệt độ của Touloukian:

P_r(T) = P_0 (P_{-1}T^{-1} + 1 + P_1 T + P_2 T^2 + P_3 T^3)

Mô hình lý thuyết đã chuyển đổi bài toán cơ học phức tạp thành hệ phương trình chủ đạo đối với hai hàm ẩn cơ bản: hàm độ võng $w(x,y,t)$ và hàm ứng suất Airy $f(x,y,t)$. Đóng góp lý thuyết mang tính đột phá bao gồm:

1. Thiết lập biểu thức giải tích tường minh xác định lực nén tới hạn $F_{xb}$, độ chênh nhiệt độ tới hạn $\Delta T_{xb}$, tần số dao động tự do riêng $\omega_{mn}$ và tần số dao động phi tuyến $\omega_{NL}$.
2. Xây dựng công thức tổng quát mô tả ứng xử sau tới hạn kiểu rẽ nhánh (bifurcation buckling) có kể đến độ không hoàn hảo hình học ban đầu $w^*(x,y) = \mu h \sin(\alpha_m x)\sin(\beta_n y)$ với $\mu \in [0, 1]$.
3. Mô hình hóa chính xác hiện tượng ghép nối cơ – nhiệt – điện trong tấm ES-FGM áp điện, giải thích cơ chế điện áp ngoài $V_0$ tác động ngược lại mômen uốn và lực màng để triệt tiêu dao động.

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

Khung phân tích của luận án là sự tích hợp đa tầng giữa:

• Lý thuyết tấm nhiều cấp độ: Lý thuyết tấm cổ điển CPT (cho tấm mỏng), Lý thuyết biến dạng trượt bậc nhất FSDT (Mindlin/Reissner) và Lý thuyết biến dạng trượt bậc ba TSDT (Reddy) với giả thiết trường chuyển vị thỏa mãn điều kiện ứng suất cắt bằng không tại hai bề mặt tự do $z = \pm h/2$ mà không cần đến hệ số hiệu chỉnh trượt $K_s$.
• Kỹ thuật san đều gân Lekhnitsky mở rộng: Chuyển đổi hệ thống gân rời rạc thành lớp vật liệu liên tục tương đương có xét đến độ cứng uốn, kéo/nén và khoảng cách độ lệch tâm $z_1^T, z_2^T$ bị biến dạng nhiệt:

d_1^T = d_1(1 + \alpha_m \Delta T), \quad h_1^T = h_1(1 + \alpha_m \Delta T)

• Mô hình nền đàn hồi hai tham số Pasternak: Thiết lập quan hệ phản lực nền $q_e(x, y) = k_1 w - k_2 \nabla^2 w$, trong đó $k_1$ đặc trưng cho độ cứng chống lún Winkler và $k_2$ phản ánh tương tác cắt giữa các phần tử nền.
• Điều kiện biên xác định tường minh: Khảo sát đầy đủ 3 trạng thái biên trong mặt phẳng: Bốn cạnh tựa bản lề tự do di chuyển (Freely Movable - FM), Bốn cạnh tựa bản lề cố định di chuyển (Immovable - IM) thỏa mãn theo nghĩa trung bình tích phân:

\int_0^b \int_0^a \frac{\partial u}{\partial x} dxdy = 0, \quad \int_0^a \int_0^b \frac{\partial v}{\partial y} dydx = 0

và trường hợp biên hỗn hợp hai cạnh tự do, hai cạnh cố định.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan khoa học thực chứng (Positivism) kết hợp phương pháp suy luận diễn dịch toán lý nghiêm ngặt (deductive mathematical physics). Phương pháp nghiên cứu chủ đạo là phương pháp bán giải tích (semi-analytical approach), kết hợp tối ưu giữa nghiệm giải tích tường minh (analytical closed-form formulation) và phương pháp số giải hệ phương trình vi phân phi tuyến (numerical integration).

Thiết kế nghiên cứu đa tầng (Multi-level design) bao gồm:

• Cấp độ 1 (Mô hình hóa vi mô - vi cơ học): Đồng nhất hóa tính chất vật liệu FGM theo quy luật phân bố hàm mũ (P-FGM), hàm Sigmoid (S-FGM mô hình I: Ceramic-Kim loại-Ceramic và mô hình II: Kim loại-Ceramic-Kim loại), hàm siêu việt (E-FGM) kết hợp quy luật phụ thuộc nhiệt độ Touloukian.
• Cấp độ 2 (Mô hình hóa vĩ mô kết cấu): Thiết lập phương trình cân bằng tĩnh, phương trình chuyển động động lực học và phương trình tương thích biến dạng cho tấm mỏng và tấm dày có gân lệch tâm trên nền Pasternak.
• Cấp độ 3 (Khử bậc và rời rạc hóa không gian): Sử dụng phương pháp hàm ứng suất Airy kết hợp phép trực giao hóa Galerkin để đưa hệ phương trình đạo hàm riêng phi tuyến về hệ phương trình đại số phi tuyến (đối với bài toán tĩnh) và hệ phương trình vi phân thường phi tuyến bậc hai theo thời gian dạng Duffing (đối với bài toán động lực học).
• Cấp độ 4 (Tích phân số thời gian): Áp dụng thuật toán Runge-Kutta bậc 4 với bước thời gian kiểm soát chặt chẽ ($\Delta t \le 10^{-4}\text{ s}$) để truy vết quỹ đạo thời gian – độ võng và khảo sát không gian pha (phase plane trajectories).

Quy trình nghiên cứu rigorous

Quy trình giải bài toán tĩnh và động lực học được chuẩn hóa theo 6 bước nghiêm ngặt:

1. Xây dựng quan hệ hình học và trường biến dạng: Áp dụng trường chuyển vị phi tuyến von Kármán kể đến biến dạng trượt ngang bậc nhất và bậc ba:

u(x, y, z, t) = u_0(x, y, t) + z \phi_x(x, y, t) - c_1 z^3 \left(\phi_x + \frac{\partial w}{\partial x}\right)

2. Thiết lập định luật thực nghiệm Hooke mở rộng: Tích hợp mô đun đàn hồi biến thiên $E(z, T)$, hệ số Poisson $\nu(z)$ hoặc $\nu = \text{const}$, hiệu ứng nhiệt $\alpha(z, T)\Delta T$ và tensor áp điện $e_{31}, e_{32}$.
3. Tính toán các ma trận độ cứng tích phân: Xác định chính xác các hệ số độ cứng kéo/nén $I_{1j}, I_{2j}, I_{3j}$, độ cứng uốn $D_{ij}$, độ cứng ghép nối kéo-uốn $B_{ij}$ và các tích phân nhiệt $\Phi_1, \Phi_2, \Phi_{1x}, \Phi_{1y}$.
4. Áp dụng kỹ thuật trực giao Galerkin: Nhân phương trình cân bằng với hàm mode trực giao $\sin(\alpha_m x)\sin(\beta_n y)$ và tích phân trên toàn miền $(0 \le x \le a; 0 \le y \le b)$ để triệt tiêu phần dư.
5. Thuật toán lặp nhiệt: Đối với bài toán phụ thuộc nhiệt độ T-D, sử dụng thuật toán lặp cập nhật ma trận độ cứng tại mỗi bước tăng nhiệt $\Delta T^{(k+1)} = \Delta T^{(k)} + \delta T$ với điều kiện hội tụ sai số tương đối $\epsilon = |\Delta T^{(k+1)} - \Delta T^{(k)}| / \Delta T^{(k)} < 10^{-5}$.
6. Kiểm chứng độ tin cậy (Triangulation & Verification): Đối chiếu độc lập với các công trình chuẩn mực quốc tế của Shen (2009), Javaheri & Eslami (2002), Đào Văn Dũng (2013, 2014) và Hoàng Văn Tùng (2011) trong các trường hợp suy biến (tấm đẳng hướng, tấm FGM không gân, tấm không phụ thuộc nhiệt độ).

Data và phân tích

Toàn bộ chương trình tính toán số, thuật toán lặp và tích phân Runge-Kutta bậc 4 được lập trình tự động hóa trên môi trường tính toán biểu thức đại số và ma trận chuyên dụng (MATLAB/Maple).

Các thông số vật liệu chuẩn được nạp vào hệ thống:

• Silicon Nitride ($Si_3N_4$ - Ceramic): $P_0(E) = 348.43\text{ GPa}$, $P_{-1} = 0$, $P_1 = -3.070 \times 10^{-4}\text{ K}^{-1}$, $P_2 = 2.160 \times 10^{-7}\text{ K}^{-2}$, $P_3 = -8.946 \times 10^{-11}\text{ K}^{-3}$; hệ số giãn nở nhiệt $\alpha_c = 5.8723 \times 10^{-6}\text{ K}^{-1}$; khối lượng riêng $\rho_c = 2370\text{ kg/m}^3$.
• Thép không gỉ (SUS304 - Metal): $P_0(E) = 201.04\text{ GPa}$, $P_{-1} = 0$, $P_1 = 3.079 \times 10^{-4}\text{ K}^{-1}$, $P_2 = -6.534 \times 10^{-7}\text{ K}^{-2}$, $P_3 = 0$; hệ số giãn nở nhiệt $\alpha_m = 15.321 \times 10^{-6}\text{ K}^{-1}$; khối lượng riêng $\rho_m = 8166\text{ kg/m}^3$.
• Lớp áp điện (PZT-4): $E_p = 81.3\text{ GPa}$, $\nu_p = 0.31$, $d_{31} = d_{32} = -1.22 \times 10^{-10}\text{ m/V}$, $\rho_p = 7500\text{ kg/m}^3$.
• Hệ số nền không thứ nguyên: $K_1 = k_1 a^4 / D_{110} \in [0, 300]$, $K_2 = k_2 a^2 / D_{110} \in [0, 30]$.

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

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

1. Hiệu ứng gia cường vượt trội của hệ thống gân lệch tâm (ES-FGM): Sự xuất hiện của gân gia cường trực giao làm tăng độ cứng chống uốn tổng thể của tấm lên từ $140% - 220%$ so với tấm FGM trơn cùng kích thước danh định, giúp nâng cao đáng kể tải nén tới hạn $F_{xb}$ và nhiệt độ tới hạn $\Delta T_{xb}$, đồng thời làm giảm biên độ võng cực đại trong giai đoạn sau tới hạn tới $45% - 60%$.
2. Ảnh hưởng suy giảm độ cứng do nhiệt độ (T-D Effect): Nghiên cứu phát hiện rằng khi nhiệt độ tăng từ $300\text{ K}$ lên $600\text{ K}$, mô đun đàn hồi hiệu dụng của tấm giảm từ $12% - 18%$, khiến đường cong tải trọng – độ võng sau tới hạn bị sụt giảm nghiêm trọng; tải nén tới hạn của tấm trong điều kiện T-D thấp hơn từ $15.4% - 28.6%$ so với mô hình T-ID truyền thống.
3. Cơ chế ổn định hóa vượt bậc của tham số nền cắt Pasternak ($K_2$): Kết quả phân tích độ nhạy chứng minh hệ số nền Pasternak $K_2$ (tương tác tiếp tuyến) có vai trò chi phối và triệt tiêu độ võng phi tuyến mạnh hơn gấp $3.5 - 5.0$ lần so với hệ số Winkler $K_1$ (chống lún đơn thuần). Khi tăng $K_2$ từ $0$ lên $20$, tần số dao động cơ bản của tấm tăng trên $85%$.
4. Kiểm soát chủ động hiện tượng phách điều hòa (Beating Phenomenon) bằng áp điện: Lần đầu tiên phát hiện và mô tả định lượng quy luật: việc đặt điện áp âm $V_0 < 0$ (ngược chiều phân cực) lên lớp áp điện sẽ sinh ra lực kéo phụ trợ trong mặt phẳng, làm tăng tần số dao động riêng, khống chế hiện tượng phách điều hòa và dập tắt dao động cộng hưởng cưỡng bức trong thời gian ngắn hơn $40%$.
5. Tính ưu việt của mô hình S-FGM đối xứng (Ceramic - Kim loại - Ceramic): Kết cấu tấm phân bố theo hàm Sigmoid mô hình I triệt tiêu hoàn toàn sự bất đối xứng uốn – kéo tại mặt giữa, giúp tấm không bị uốn sơ khởi khi chịu nhiệt độ tăng đều, đem lại độ ổn định nhiệt cao hơn $32%$ so với tấm phân bố theo hàm lũy thừa P-FGM thông thường.

Implications đa chiều

• Về mặt học thuật: Cung cấp hệ phương trình chuẩn tắc và bộ dữ liệu benchmark giải tích chính xác cao làm thước đo chuẩn để kiểm thử các phần mềm phần tử hữu hạn thương mại (Abaqus, ANSYS, COMSOL) trong việc tính toán kết cấu composite thông minh đa trường vật lý.
• Về mặt phương pháp luận: Mở ra phương pháp tiếp cận giải tích kết hợp hàm ứng suất Airy cải tiến cho các kết cấu tấm/vỏ phức tạp có gân gia cường lệch tâm mà trước đây bắt buộc phải dùng các mô hình phần tử hữu hạn tốn kém tài nguyên tính toán.
• Về mặt kỹ thuật ứng dụng: Cung cấp quy trình thiết kế tối ưu hóa tham số hình học gân (chiều cao $h_1, h_2$, bề rộng $d_1, d_2$, khoảng cách bước gân $s_1, s_2$) giúp tiết kiệm $20% - 35%$ khối lượng vật liệu cấu trúc mà vẫn đảm bảo khả năng chịu tải giới hạn trong vỏ tên lửa và buồng đốt động cơ phản lực.
• Về mặt công nghệ điều khiển: Đề xuất giải pháp tích hợp cảm biến/thiết bị chấp hành áp điện tự động bù sai lệch hình học và kiểm soát rung động thời gian thực cho các kết cấu hàng không làm việc trong môi trường nhiệt độ biến đổi khắc nghiệt.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu (boundary conditions & limitations):

1. Giả thiết gân trực giao và kích thước gân đồng nhất: Nghiên cứu giới hạn ở mô hình gân trực giao hình chữ nhật đặt đều đặn; chưa xét đến hệ gân xiên (anisogrid/isogrid) hoặc mạng gân tối ưu hóa topology phi cấu trúc.
2. Mô hình biến dạng nhiệt tĩnh trạng thái dừng: Trường nhiệt độ khảo sát là trường nhiệt tăng đều hoặc truyền nhiệt dừng 1D qua chiều dày theo phương trình Fourier; chưa xét đến quá trình sốc nhiệt không dừng (transient thermal shock) và hiệu ứng nhiệt đàn hồi liên kết bậc cao (coupled thermoelasticity).
3. Mô hình nền đàn hồi tuyến tính: Nền Pasternak được giả thiết phản ứng hoàn toàn đàn hồi tuyến tính hai thông số; chưa tính đến tính phi tuyến hình học của nền, tính từ biến (viscoelastic foundation) hay hiện tượng bong tách cục bộ giữa tấm và nền khi dao động biên độ lớn.
4. Điều kiện biên lý thuyết: Các biên tựa khớp bản lề (FM và IM) được xử lý hoàn hảo theo điều kiện giải tích; chưa mô phỏng các điều kiện biên đàn hồi bán ngàm (elastic restraint) hoặc biên tự do phức tạp trong kết cấu thực tế.

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

• Mở rộng phân tích sang vật liệu FGM gia cường ống nano carbon (Carbon Nanotube Reinforced Composite - CNTRC) và vật liệu xốp có lỗ rỗng biến thiên (Porous FGMs).
• Phát triển mô hình giải tích cho bài toán va chạm xung kích (impact dynamics), tải nổ áp lực cao (blast loads) và hiện tượng mất ổn định đàn dẻo phi tuyến (elasto-plastic postbuckling).
• Nghiên cứu cơ chế lan truyền vết nứt vi mô và hiện tượng mỏi nhiệt động học đa trục trong kết cấu ES-FGM áp điện.

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

• Tác động học thuật quốc tế: Các kết quả của luận án đã được công bố trên các tạp chí quốc tế ISI danh tiếng (như Composite Structures, Thin-Walled Structures, Journal of Sandwich Structures & Materials, Mechanics of Advanced Materials and Structures), đóng góp vào việc nâng cao chỉ số trích dẫn và vị thế của trường phái Cơ học tính toán Việt Nam trên bản đồ khoa học thế giới.
• Chuyển giao công nghệ công nghiệp: Cung cấp cơ sở khoa học trực tiếp cho các viện nghiên cứu hàng không vũ trụ, công nghiệp đóng tàu quân sự và năng lượng nguyên tử trong việc lựa chọn vật liệu và thiết kế tấm chắn nhiệt thế hệ mới.
• Hiệu quả xã hội và môi trường: Việc tối ưu hóa vật liệu FGM gia cường gân giúp giảm đáng kể trọng lượng kết cấu phương tiện bay, từ đó cắt giảm lượng tiêu thụ nhiên liệu hóa thạch và giảm phát thải khí nhà kính trong ngành giao thông vận tải hàng không.

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

• Nghiên cứu sinh và Giảng viên đại học: Tiếp cận khung lý thuyết giải tích giải quyết hệ phương trình vi phân phi tuyến bậc cao; làm tài liệu tham khảo mẫu mực cho phương pháp nghiên cứu cơ học vật rắn biến dạng.
• Kỹ sư R&D Hàng không – Không gian: Sử dụng các công thức giải tích đóng để tính toán sơ bộ nhanh chóng (rapid conceptual design) độ an toàn của cánh máy bay siêu thanh, vỏ khoang đổ bộ và tấm pin mặt trời vệ tinh.
• Các nhà hoạch định chính sách khoa học công nghệ: Có cơ sở dữ liệu thực nghiệm – tính toán để xây dựng tiêu chuẩn kỹ thuật quốc gia về thử nghiệm và ứng dụng vật liệu composite tiên tiến trong an ninh quốc phòng.

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 mở rộng thành công lý thuyết tấm biến dạng trượt bậc ba (TSDT) của Reddy kết hợp với kỹ thuật san đều gân Lekhnitsky cải tiến có xét đến sự biến dạng hình học và ứng suất nhiệt riêng biệt của gân để thiết lập hệ phương trình vi phân phi tuyến đóng giải bài toán ổn định và động lực học tấm dày ES-FGM áp điện trong môi trường nhiệt độ biến thiên (T-D).
2. Tính 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 phương pháp khai triển tiệm cận của Shen (2009) [101] và phương pháp năng lượng thuần túy của Javaheri & Eslami (2002) [76], luận án sử dụng phương pháp hàm ứng suất Airy kết hợp phép trực giao hóa Galerkin và thuật toán Runge-Kutta bậc 4, giải quyết tường minh được cả 3 chế độ biên trong mặt phẳng (FM, IM, hỗn hợp) và thu được nghiệm giải tích dạng đóng cho cả tấm mỏng (CPT) và tấm dày (FSDT, TSDT).
3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu là gì? Khi tăng nhiệt độ môi trường $\Delta T > 0$, các tấm có độ không hoàn hảo hình học ban đầu ($w^* \neq 0$) bị xuất hiện độ võng tĩnh ngay lập tức tại thời điểm tải nén cơ học bằng $0$ ($F_x = 0$), đồng thời làm thay đổi căn bản quy luật phân nhánh ổn định từ mất ổn định rẽ nhánh cổ điển sang mất ổn định theo dạng uốn - nén kết hợp liên tục không có điểm uốn gập sắc nét.
4. Quy trình tái lập kết quả (Replication Protocol) có được đảm bảo không? Hoàn toàn minh bạch và khả thi: Toàn bộ hệ thống ma trận hệ số độ cứng ($A_{ij}, B_{ij}, D_{ij}$), công thức tích phân nhiệt ($\Phi_i$), các bảng thông số vật liệu phụ thuộc nhiệt độ Touloukian ($P_0, P_1, P_2, P_3$) và thuật toán giải số Runge-Kutta đều được công bố chi tiết trong các công thức tường minh và phụ lục luận án.
5. Chương trình nghiên cứu 10 năm tới mở ra hướng đi nào? Phát triển lý thuyết cơ học đa tỷ lệ (multiscale mechanics) cho kết cấu FGM cấu trúc nano (nanoplate/nanoshell) có xét đến hiệu ứng kích thước phi cục bộ (nonlocal elasticity theory) kết hợp trí tuệ nhân tạo (AI/Machine Learning) để tối ưu hóa trường vật liệu gradient thời gian thực.

Kết luận

1. Xây dựng hoàn chỉnh hệ phương trình cơ bản và phương pháp giải tích phân tích phi tuyến tĩnh và động lực học của tấm chữ nhật mỏng và dày làm từ vật liệu FGM có và không có gân gia cường lệch tâm (ES-FGM), tấm FGM áp điện và tấm FGM đối xứng (S-FGM) trên nền đàn hồi Pasternak.
2. Chứng minh định lượng tính tất yếu của việc xét đến sự phụ thuộc nhiệt độ của cơ tính vật liệu (T-D) trong tính toán ổn định kết cấu chịu nhiệt độ cao.
3. Khám phá vai trò chi phối của hệ số nền cắt Pasternak $K_2$ và khả năng kiểm soát chủ động hiện tượng phách dao động bằng điện áp áp điện $V_0$.
4. Cung cấp hệ thống nghiệm giải tích đóng chính xác cao, đóng vai trò là dữ liệu chuẩn mực (benchmark data) phục vụ kiểm chứng các mô hình số và phần mềm tính toán kết cấu tiên tiến.
5. Đặt nền tảng lý thuyết vững chắc mở ra các hướng nghiên cứu hiện đại về vật liệu composite chức năng thông minh đa trường vật lý phục vụ công nghiệp công nghệ cao và quốc phòng an ninh.