Tổng quan về luận án

Nghiên cứu cơ học kết cấu tiên tiến đóng vai trò then chốt trong việc tối ưu hóa độ an toàn và hiệu suất vận hành của các cấu kiện hàng không vũ trụ, lò phản ứng hạt nhân và kỹ thuật biển chịu tải trọng khắc nghiệt. Luận án tiến sĩ chuyên ngành Cơ Kỹ thuật (Mã số: 62 52 01 01) của tác giả Phạm Hồng Công, dưới sự hướng dẫn khoa học của GS. Nguyễn Đình Đức tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (2018), đã phát triển hệ thống phân tích phi tuyến tĩnh và động lực học của kết cấu tấm chữ nhật làm bằng vật liệu có cơ tính biến đổi (Functionally Graded Material - FGM) đặt trên nền đàn hồi. Điểm tiên phong của công trình là thiết lập lời giải giải tích tường minh cho các tấm FGM có gân gia cường lệch tâm (ES-FGM), tấm FGM đối xứng phân bố theo hàm Sigmoid (S-FGM), và tấm ES-FGM tích hợp lớp áp điện, trong điều kiện tính chất vật liệu của cả tấm và gân đều phụ thuộc phi tuyến vào nhiệt độ (Temperature-Dependent - T-D).

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định từ thực tiễn: phần lớn công bố quốc tế trước đó (như các nghiên cứu của Javaheri & Eslami, 2002; Sofiyev, 2010; Shen Hui-Shen, 2009) chủ yếu giới hạn ở tấm vỏ FGM trơn không gân, hoặc chỉ áp dụng lý thuyết tấm cổ điển (CPT) với giả định hệ số Poisson $\nu$ không đổi và tính chất vật liệu độc lập nhiệt độ (T-ID). Việc giải quyết bài toán phi tuyến cho tấm dày ES-FGM chịu đồng thời tải trọng nén cơ học, trường nhiệt độ phi tuyến và hiệu ứng điện trường thông qua phương pháp giải tích giải toán học chính xác kết hợp lý thuyết biến dạng trượt bậc nhất (FSDT) và bậc ba (TSDT) vẫn là bài toán mở đầy thách thức do tính cồng kềnh của hệ phương trình vi phân phi tuyến liên kết.

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

  1. Câu hỏi nghiên cứu 1 (RQ1): Sự tương tác phi tuyến giữa độ cứng gân lệch tâm, phản lực nền Pasternak hai tham số và tính chất biến đổi nhiệt độ ảnh hưởng như thế nào đến tải tới hạn và quỹ đạo sau tới hạn (postbuckling) của tấm mỏng và dày FGM?
  2. Câu hỏi nghiên cứu 2 (RQ2): Quy luật động lực học, tần số dao động tự do tuyến tính, quan hệ biên độ – tần số phi tuyến và đáp ứng thời gian – độ võng của tấm S-FGM và tấm ES-FGM áp điện diễn biến ra sao dưới tác dụng của tải trọng điều hòa ngoài?
  3. Giả thuyết nghiên cứu 1 (H1): Việc tích hợp gân trực giao lệch tâm kết hợp nền Pasternak sẽ gia tăng đáng kể khả năng chịu tải tới hạn nén và nhiệt, đồng thời bù đắp sự suy giảm mô đun đàn hồi do hiệu ứng suy thoái nhiệt gây ra.
  4. Giả thuyết nghiên cứu 2 (H2): Sự xuất hiện của điện áp điều khiển ngoài trên lớp áp điện kết hợp với biến dạng trượt bậc ba (TSDT) cho phép triệt tiêu hoặc kiểm soát tối ưu hiện tượng phách điều hòa (beat phenomenon) và độ võng động lực học.

Khung lý thuyết của luận án tích hợp: Lý thuyết tấm cổ điển Kirchhoff-Love (CPT), Lý thuyết biến dạng trượt bậc nhất Mindlin (FSDT), Lý thuyết biến dạng trượt bậc cao Reddy (TSDT), Quan hệ biến dạng phi tuyến hình học von Kármán, Tiêu chuẩn ổn định rẽ nhánh Euler-Poincaré, và Mô hình nền đàn hồi Pasternak hai hệ số ($k_1, k_2$). Đóng góp đột phá của luận án được lượng hóa rõ nét qua việc xây dựng thành công công thức gân tổng quát mở rộng từ kỹ thuật san đều Lekhnitsky có xét ứng suất nhiệt trong gân, cho phép nâng cao khả năng dự báo tải tới hạn chính xác hơn 12-18% so với mô hình bỏ qua biến dạng nhiệt của gân, đồng thời khảo sát định lượng toàn diện trên các kích thước tấm tiêu chuẩn ($a = b = 1.0\text{ m}$, tỷ lệ $a/h = 20 \div 100$).

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy sự tiến triển mạnh mẽ của cơ học vật liệu FGM kể từ khi được khởi xướng tại Viện Sendai, Nhật Bản vào năm 1984. Dòng nghiên cứu chính trên thế giới tập trung vào phân tích phi tuyến của tấm vỏ FGM không gân. Tiêu biểu là nhóm nghiên cứu của Shen Hui-Shen (2009, 2011, 2014) sử dụng phương pháp khai triển tiệm cận tham số bé kết hợp lý thuyết lớp biên để dự đoán ứng xử sau tới hạn của vỏ trụ và tấm nano composite FGM trong môi trường nhiệt. Cùng thời điểm, Sofiyev và Schnack (2012, 2014) triển khai FSDT và phương pháp Galerkin để đánh giá sự mất ổn định nhiệt đàn hồi của vỏ nón cụt FGM. Nhóm tác giả Eslami, Javaheri và Shariat (2002, 2005, 2006) áp dụng phương pháp hàm năng lượng giải tích cho tấm FGM chịu tải cơ - nhiệt dựa trên CPT và HSDT.

Tại Việt Nam, các công trình tiên phong của Nguyễn Đình Đức, Đào Huy Bích, Hoàng Văn Tùng, Đào Văn Dũng (2011–2016) đã đặt nền móng vững chắc cho việc giải quyết trạng thái sau tới hạn của kết cấu composite và FGM không hoàn hảo. Đào Huy Bích và cộng sự (2011) đề xuất mô hình gân trực giao đồng nhất để đảm bảo tính liên tục vật liệu. Đào Văn Dũng và Hoàng Thị Thiêm (2014, 2016) nghiên cứu tấm và vỏ FGM có gân chịu tải cơ bằng phương pháp hàm chuyển vị. Nguyễn Đình Đức và cộng sự (2013, 2015) lần đầu tiên đưa ra phương pháp xác định ứng suất nhiệt và sự biến đổi hình học trong gân dưới tác động của nhiệt độ.

Trong y văn tồn tại hai trường phái tranh luận chính: Trường phái thứ nhất (Shen Hui-Shen, Amabili) thiên về phương pháp xấp xỉ năng lượng và khai triển tiệm cận số cho kết cấu trơn, cho rằng việc gia cường bằng gân phức tạp hóa điều kiện biên và có thể giải quyết xấp xỉ qua phần tử hữu hạn. Trường phái thứ hai (Đào Huy Bích, Nguyễn Đình Đức, Đào Văn Dũng) khẳng định giải pháp giải tích tường minh sử dụng phương pháp hàm ứng suất và kỹ thuật san đều gân Lekhnitsky cải tiến là tối ưu để nắm bắt bản chất vật lý của sự mất ổn định rẽ nhánh mà không bị phụ thuộc vào lưới sai số số học. Luận án của NCS. Phạm Hồng Công định vị chính xác tại giao điểm này, mở rộng ranh giới nghiên cứu bằng cách kết hợp lý thuyết TSDT bậc cao với phương pháp hàm ứng suất Airy tổng quát, khắc phục triệt để nhược điểm của CPT đối với tấm dày, đồng thời giải quyết trọn vẹn bài toán ghép cặp cơ - nhiệt - điện mà các nghiên cứu quốc tế trước đây chưa công bố dạng nghiệm giải tí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 và thách thức trực tiếp các giới hạn của ba lý thuyết cơ học kinh điển: Lý thuyết tấm cổ điển Kirchhoff (CPT) được mở rộng cho bài toán biến thiên hệ số Poisson theo chiều dày $\nu = \nu(z)$; Lý thuyết biến dạng trượt bậc nhất Mindlin (FSDT) và Lý thuyết biến dạng trượt bậc ba Reddy (TSDT) được mở rộng thành công cho tấm dày chức năng có cấu trúc dị tính phức hợp (gân lệch tâm + lớp áp điện + nền đàn hồi hai thông số). Khung khái niệm (conceptual framework) của luận án thiết lập mối liên kết hữu cơ giữa phân bố tỷ lệ thể tích vật liệu (theo quy luật lũy thừa P-FGM, đối xứng Sigmoid S-FGM, và siêu việt E-FGM) với trường chuyển vị phi tuyến, tensor ứng suất nhiệt và lực tương tác tiếp xúc nền.

Mô hình lý thuyết dẫn xuất hệ 04 mệnh đề khoa học cơ bản:

  • Mệnh đề 1 (P1): Khi tỷ lệ chiều dày tấm $a/h < 50$, việc sử dụng CPT dẫn đến đánh giá quá cao độ cứng kết cấu và gây sai số tải tới hạn từ 8.5% đến 22.4% so với TSDT do bỏ qua hiệu ứng trượt ngang.
  • Mệnh đề 2 (P2): Ứng suất nhiệt trong hệ gân gia cường lệch tâm tỷ lệ thuận với gradient nhiệt độ và đóng vai trò như một trường lực nén/kéo ban đầu, làm dịch chuyển điểm rẽ nhánh ổn định của kết cấu.
  • Mệnh đề 3 (P3): Điện trường điều khiển $V_0$ tỷ lệ thuận với khả năng triệt tiêu dao động biên độ lớn thông qua mômen uốn áp điện kháng lại độ võng uốn cơ học.
  • Mệnh đề 4 (P4): Cấu trúc S-FGM đối xứng hai pha ceramic-kim loại-ceramic triệt tiêu hoàn toàn sự ghép cặp uốn-kéo ở mặt giữa, tạo ra độ ổn định động học vượt trội so với cấu trúc P-FGM bất đối xứng.

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

Sự đột phá phương pháp luận thể hiện qua việc tích hợp đồng thời ba cấu trúc lý thuyết: (1) Quy tắc hỗn hợp mở rộng Voigt phụ thuộc nhiệt độ cho vật liệu phân cấp; (2) Kỹ thuật san đều gân Lekhnitsky hiệu chỉnh tính đến biến dạng nhiệt hình học của tiết diện gân ($s_1^T, s_2^T, d_1^T, h_1^T$); và (3) Lý thuyết biến dạng trượt bậc cao TSDT không yêu cầu hệ số hiệu chỉnh trượt ($K_s = 5/6$).

Cách tiếp cận phân tích mới sử dụng phương pháp hàm ứng suất Airy $f(x,y)$ để chuyển đổi hệ phương trình cân bằng liên kết phi tuyến phức tạp thành hệ phương trình vi phân tương thích bậc bốn, cho phép áp dụng phép chiếu trực giao Galerkin để tìm nghiệm giải tích dạng tường minh:

$$Fx = b_{11}W + b_{12}\frac{W}{\overline{W}} + b_{13}\overline{W}^2 + b_{14}W\overline{W}^2$$

Điều kiện biên nghiên cứu được xác lập tường minh cho ba cấu hình kỹ thuật: Bốn cạnh tựa bản lề tự do dịch chuyển trong mặt phẳng (Freely Movable - FM), Bốn cạnh tựa bản lề cố định dịch chuyển trong mặt phẳng (Immovable - IM), và Cấu hình hỗn hợp hai cạnh FM kết hợp hai cạnh IM.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết lý khoa học thực chứng suy diễn (Deductive Positivism), kết hợp chặt chẽ giữa mô hình hóa giải tích toán lý thuyết (analytical mathematical modeling) và mô phỏng số học nghiệm thu (numerical computational verification). Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm ba cấp độ phân tích: Cấp độ vi mô (xác định tính chất vật liệu hiệu dụng thông qua quy luật hàm mũ và hàm phi tuyến nhiệt độ), Cấp độ trung mô (mô hình hóa tương tác cấu trúc giữa tấm FGM, gân trực giao và nền Pasternak), và Cấp độ vĩ mô (xác lập đáp ứng phi tuyến tĩnh và động lực học toàn cục của kết cấu).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu lý thuyết tuân thủ các bước chuẩn hóa nghiêm ngặt:

  1. Dữ liệu vật liệu đầu vào: Trích xuất chính xác bộ tham số nhiệt độ thực nghiệm của Silicon Nitride ($Si_3N_4$) và Thép không gỉ ($SUS304$) theo tiêu chuẩn Touloukian:
    • $Si_3N_4$: $P_0 = 348.43\times 10^9\text{ Pa}, 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}$.
    • $SUS304$: $P_0 = 201.04\times 10^9\text{ Pa}, 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$.
  2. Kiểm định chéo (Triangulation): Kết hợp đa dạng phương pháp lý thuyết (CPT, FSDT, TSDT) và so sánh kết quả suy biến (khi triệt tiêu gân và nền) với các công bố quốc tế uy tín.
  3. Độ tin cậy và chuẩn hóa: Độ hội tụ của thuật toán số Runge-Kutta bậc 4 đạt sai số tuyệt đối nhỏ hơn $\epsilon = 10^{-6}$, bước tích phân thời gian $\Delta t = 0.0001\text{ s}$ đảm bảo tính ổn định số học hoàn toàn.

Data và phân tích

Dữ liệu mô phỏng được lập trình tính toán tự động trên môi trường phần mềm tính toán ký hiệu và số trị khoa học. Kích thước hình học tiêu chuẩn: $a = b = 1.0\text{ m}$, chiều dày $h = 0.01 \div 0.05\text{ m}$; thông số gân: bề rộng $d_1 = d_2 = 0.005\text{ m}$, chiều cao $h_1 = h_2 = 0.015\text{ m}$, khoảng cách gân $s_1 = s_2 = 0.1\text{ m}$; hệ số nền đàn hồi không thứ nguyên $K_1 = k_1 a^4 / D_0 = 100 \div 300$, $K_2 = k_2 a^2 / D_0 = 10 \div 30$. Phân tích độ nhạy (robustness check) được thực hiện với các chỉ số phân bố thể tích $N \in [0, 5]$, tỷ số dạng hình học $a/b \in [0.5, 2.0]$, và tham số không hoàn hảo ban đầu $\mu = W_0/h \in [0, 0.5]$.

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

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

Công trình công bố 05 phát hiện khoa học mang tính đột phá với bằng chứng số liệu định lượng chính xác:

STT Phát hiện then chốt Dữ liệu định lượng & Bằng chứng thực nghiệm Ý nghĩa cơ học
01 Hiệu ứng gia cường của gân lệch tâm Tải nén tới hạn của tấm ES-FGM tăng từ 35.2% đến 64.8% so với tấm FGM trơn cùng độ dày danh định. Gân gia cường trực giao làm tăng đáng kể mômen quán tính kháng uốn mà không làm tăng nhiều khối lượng tổng thể kết cấu.
02 Tác động vượt trội của hệ số nền Pasternak Tăng hệ số cắt Pasternak $K_2$ từ 10 lên 20 làm tăng tần số dao động tự do cơ bản $\omega_m$ lên 28.4%, hiệu quả gấp 3.2 lần so với việc tăng hệ số Winkler $K_1$. Phản lực tiếp tuyến trượt của nền Pasternak có khả năng phân tán năng lượng biến dạng uốn hiệu quả hơn nhiều so với phản lực thẳng đứng độc lập.
03 Suy giảm độ cứng do nhiệt độ (T-D vs T-ID) Mô hình phụ thuộc nhiệt độ T-D dự báo tải tới hạn nhiệt thấp hơn từ 14.6% đến 21.3% so với mô hình T-ID tại mức tăng nhiệt độ $\Delta T = 300\text{ K}$. Bỏ qua sự suy giảm nhiệt của mô đun đàn hồi $E(T)$ sẽ đánh giá quá lạc quan khả năng chịu tải của tấm, gây nguy cơ mất an toàn kết cấu.
04 Hiện tượng phách điều hòa (Beat Phenomenon) Biên độ dao động bị điều chế mạnh mẽ khi tần số lực cưỡng bức $\Omega$ tiến gần tần số tự nhiên $\omega_{11}$ ($\Omega/\omega_{11} \approx 0.95 \div 1.05$). Hiện tượng tích tụ năng lượng cộng hưởng cục bộ đòi hỏi phải thiết lập cơ chế kiểm soát chủ động qua vật liệu áp điện.
05 Khả năng triệt tiêu dao động bằng điện áp áp điện Đặt điện áp điều khiển $V_0 = 200\text{ V}$ lên lớp áp điện làm giảm 32.5% biên độ độ võng cực đại trong phân tích động lực học thời gian thực. Lớp áp điện đóng vai trò cảm biến và cơ cấu chấp hành hiệu quả cao trong việc ổn định động học kết cấu thông minh.

Kết quả so sánh trạng thái sau tới hạn cho thấy đường cong độ võng – tải trọng của tấm không hoàn hảo ($\mu > 0$) bắt đầu biến dạng uốn ngay từ khi tải trọng $F_x > 0$, không tồn tại điểm rẽ nhánh đột ngột như tấm hoàn hảo lý tưởng ($\mu = 0$), phản ánh chính xác hành vi thực tế của kết cấu trong gia công chế tạo.

       Độ võng W/h
              Tới hạn ES-FGM         Tới hạn FGM trơn

Implications đa chiều

Về mặt lý thuyết, nghiên cứu hoàn thiện khung giải tích TSDT cho các cấu trúc FGM phức hợp, tạo nền tảng lý thuyết vững chắc cho cơ học vật liệu biến tính. Về mặt phương pháp luận, thuật toán kết hợp Galerkin và Runge-Kutta bậc 4 mở ra hướng tiếp cận chuẩn tắc cho các bài toán phi tuyến ghép cặp đa trường (cơ - nhiệt - điện). Về ứng dụng thực tiễn, kết quả cung cấp bộ thông số thiết kế tối ưu ($N = 1.0 \div 2.0$, bố trí khoảng cách gân $s/h = 10 \div 20$) cho các kỹ sư chế tạo vỏ động cơ phản lực không gian, tấm bảo vệ nhiệt trong lò phản ứng và kết cấu tàu ngầm thế hệ mới. Về mặt chính sách và tiêu chuẩn, luận án cung cấp căn cứ khoa học chuẩn xác để cập nhật các quy chuẩn kỹ thuật thiết kế kết cấu composite tiên tiến chịu nhiệt độ cao tại Việt Nam và khu vực.

Limitations và Future Research

Luận án thừa nhận một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn hình học và trường phân bố: Nghiên cứu tập trung vào tấm chữ nhật phẳng; chưa mở rộng cho các cấu trúc vỏ có độ cong kép phức tạp hoặc tấm có lỗ khoét hình học (cutouts).
  2. Giả định trường nhiệt độ: Trường nhiệt độ được giả thiết tăng đều hoặc phân bố một chiều dừng qua chiều dày; chưa xét đến trường nhiệt độ động không dừng (transient thermal shock) kết hợp truyền nhiệt bức xạ phi tuyến.
  3. Vật liệu gân: Gân gia cường được giả định là vật liệu đồng nhất đẳng hướng; chưa xét đến gân làm từ vật liệu composite cốt sợi nano (CNTRC) hoặc gân có cơ tính biến đổi đa hướng (2D-FGM).

Chương trình nghiên cứu 5-10 năm tới bao gồm các hướng ưu tiên:

  • Hướng 1: Phát triển mô hình giải tích cho tấm và vỏ FGM gia cường bằng ống nano carbon (CNT) và tấm nano graphene (GPL) phân bố theo hướng hàm số trên nền đàn hồi phi tuyến Winkler-Pasternak-Kerr.
  • Hướng 2: Mở rộng lý thuyết biến dạng trượt bậc cao tự do (quasi-3D shear deformation theories) để phân tích dao động ngẫu nhiên và ổn định động học dưới tác động của sóng xung kích nổ.
  • Hướng 3: Ứng dụng thuật toán tối ưu hóa đa mục tiêu (Machine Learning / Genetic Algorithms) để tối ưu hóa tỷ lệ thể tích $N$ và mật độ gân nhằm tối đa hóa tải tới hạn và giảm thiểu khối lượng cấu kiện.

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

Luận án tạo ra tác động học thuật sâu sắc với tiềm năng trích dẫn cao trong các tạp chí quốc tế hàng đầu danh mục ISI/Scopus (Q1) như Composite Structures, Thin-Walled Structures, Mechanics of Advanced Materials and Structures, Applied Mathematical Modelling. Trong công nghiệp, kết quả nghiên cứu thúc đẩy việc chuyển giao công nghệ thiết kế cánh đón gió tuabin khí công nghiệp và tấm chắn nhiệt tàu vũ trụ tại các viện nghiên cứu đầu ngành. Lợi ích kinh tế - xã hội được khẳng định thông qua việc tối ưu hóa lượng vật liệu ceramic quý hiếm đắt đỏ, giúp giảm từ 15% đến 25% chi phí sản xuất cấu kiện composite nhiệt độ cao mà vẫn bảo đảm hệ số an toàn chịu lực nghiêm ngặt.

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

Hệ thống kết quả nghiên cứu mang lại giá trị thiết thực cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp thiết lập phương trình vi phân phi tuyến giải tích tường minh, kỹ thuật xử lý hàm ứng suất Airy và giải thuật tích phân thời gian Runge-Kutta bậc 4.
  • Các nhà khoa học cơ học tính toán: Kế thừa khung lý thuyết TSDT mở rộng và mô hình gân Lekhnitsky cải tiến để giải quyết các bài toán kết cấu phức tạp hơn.
  • Kỹ sư R&D hàng không và quốc phòng: Sở hữu công thức đóng để tra cứu nhanh tải trọng tới hạn, tần số tự nhiên và đáp ứng biên độ động lực học khi thiết kế vỏ khí tài quân sự.
  • Cơ quan quản lý và hoạch định tiêu chuẩn: Tiếp cận dữ liệu thực chứng tin cậy để biên soạn tiêu chuẩn kiểm định độ bền nhiệt - cơ của vật liệu mới 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 mở rộng thành công Lý thuyết biến dạng trượt bậc ba của Reddy (TSDT) kết hợp hàm ứng suất Airy và kỹ thuật san đều gân Lekhnitsky có xét ứng suất nhiệt nội tại của gân và tấm phụ thuộc nhiệt độ (T-D). Đây là lời giải giải tích đóng đầu tiên giải quyết trọn vẹn sự tương tác phi tuyến cơ - nhiệt của tấm dày ES-FGM trên nền Pasternak.

  2. Cải tiến phương pháp luận so với các nghiên cứu trước đây như thế nào? So với nghiên cứu của Javaheri & Eslami (2002) chỉ dùng CPT cho tấm trơn và Đào Văn Dũng & Hoàng Thị Thiêm (2014) dùng phương pháp hàm chuyển vị cho tấm có gân chịu tải cơ, luận án thiết lập phương pháp hàm ứng suất toàn diện kết hợp giải thuật tích phân Runge-Kutta bậc 4, giải quyết đồng thời tải cơ, tải nhiệt và điều khiển điện áp áp điện với độ chính xác cao hơn hẳn.

  3. Phát hiện bất ngờ và thú vị nhất thu được từ dữ liệu là gì? Sự xuất hiện của hiện tượng "uốn trước tới hạn" đối với tấm không hoàn hảo chịu tải nhiệt: tấm bắt đầu võng ngay khi $\Delta T > 0$ dù chưa có tải nén dọc trục. Đồng thời, hệ số trượt cắt Pasternak ($K_2$) thể hiện vai trò ổn định kết cấu vượt trội gấp 3.2 lần so với độ cứng lò xo Winkler ($K_1$).

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Toàn bộ hệ thống phương trình chủ đạo, biểu thức giải tích tường minh cho các hệ số ($b_{ij}$, $L$, $I_{ij}$), bảng thông số vật liệu phụ thuộc nhiệt độ Touloukian của $Si_3N_4/SUS304$, và giải thuật lặp số đều được trình bày chi tiết trong văn bản và phụ lục, cho phép độc giả tái lập 100% kết quả số.

  5. Chương trình nghiên cứu 10 năm được phác thảo ra sao? Kế hoạch tập trung vào việc tích hợp vật liệu nano composite (graphene/CNT), mở rộng bài toán sang môi trường ăn mòn đa trường vật lý (hygro-thermo-electro-mechanical), và tích hợp trí tuệ nhân tạo (AI-driven structural optimization) để tự động hóa thiết kế kết cấu FGM thông minh.

Kết luận

Luận án tiến sĩ của NCS. Phạm Hồng Công đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 05 đóng góp cốt lõi:

  1. Xây dựng hoàn chỉnh hệ phương trình điều khiển giải tích cho tấm mỏng và dày ES-FGM trên nền Pasternak dựa trên CPT, FSDT và TSDT.
  2. Dẫn xuất công thức tường minh xác định lực nén tới hạn, nhiệt độ tới hạn và đường cong độ võng – tải trọng giai đoạn sau tới hạn cho cả kết cấu hoàn hảo và không hoàn hảo.
  3. Thiết lập hệ phương trình dao động phi tuyến và xác định chính xác tần số tự do tuyến tính, quan hệ tần số – biên độ và đáp ứng thời gian thực của tấm S-FGM và tấm ES-FGM áp điện.
  4. Chứng minh định lượng vai trò quyết định của hiệu ứng phụ thuộc nhiệt độ (T-D), mật độ gân gia cường và điện áp điều khiển áp điện trong việc kiểm soát độ võng và triệt tiêu dao động phách điều hòa.
  5. Cung cấp bộ công cụ tính toán số học tin cậy, được thẩm định chéo với các công bố quốc tế, phục vụ trực tiếp cho công tác nghiên cứu và thiết kế kỹ thuật công nghệ cao.

Công trình mở ra các hướng nghiên cứu mới về vật liệu nano chức năng và kết cấu thông minh, khẳng định vị thế học thuật của cơ học công trình Việt Nam trên trường quốc tế.