Tổng quan về luận án

Trong kỷ nguyên phát triển của kỹ thuật hàng không vũ trụ, công nghệ hạt nhân và các hệ thống kết cấu chịu nhiệt độ cao, vật liệu cơ tính biến thiên (Functionally Graded Materials - FGM) đã trở thành giải pháp vật liệu tiên tiến bậc nhất để thay thế composite phân lớp truyền thống. Cốt lõi ưu việt của FGM là sự biến đổi trơn và liên tục về thành phần pha từ gốm (ceramic) sang kim loại (metal), giúp triệt tiêu tập trung ứng suất tại mặt phân giới và loại bỏ hoàn toàn hiện tượng bong tách lớp (delamination). Tuy nhiên, việc ứng dụng FGM vào các kết cấu vỏ tròn xoay có hình học phức tạp như vỏ nón cụt (truncated conical shells) có gân gia cường đặt trên nền đàn hồi chịu tải phức hợp cơ - nhiệt - điện vẫn đang đối mặt với những thách thức toán - cơ học to lớn. Như tác giả đã khẳng định trong văn bản: "Việc nghiên cứu các bài toán về vỏ tròn xoay như vỏ nón, vỏ cầu, vỏ trống, vỏ parabolic… dẫn đến hệ phương trình đạo hàm riêng có hệ số là hàm của tọa độ, do vậy tìm nghiệm giải tích của chúng khó khăn về toán học".

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các công trình quốc tế tiêu biểu như Sofiyev (2010, 2014) hay Eslami cùng cộng sự (2008) chủ yếu tập trung vào vỏ nón FGM đẳng hướng hoặc không có gân gia cường, và đa phần dừng lại ở việc áp dụng lý thuyết vỏ cổ điển (Classical Shell Theory - CST). Lý thuyết CST bỏ qua hoàn toàn biến dạng trượt ngang ($\gamma_{xz0} = \gamma_{\theta z0} = 0$), do đó chỉ chính xác đối với kết cấu vỏ rất mỏng và sẽ dẫn đến sai số nghiêm trọng (đánh giá quá cao tải trọng tới hạn $P_{cr}$) khi phân tích vỏ dày vừa và vỏ dày. Hơn nữa, việc xem xét ảnh hưởng tương hỗ đồng thời giữa gân gia cường trực giao (gân dọc đường sinh và gân vòng) có khoảng cách biến thiên $d_s(x)$, nền đàn hồi Pasternak hai tham số ($K_1, K_2$), trường nhiệt phi tuyến và lớp áp điện tích hợp vẫn chưa được giải quyết một cách tường minh bằng phương pháp giải tích.

Để giải quyết khoảng trống này, luận án xác định hai câu hỏi nghiên cứu (Research Questions - RQ) và hai giả thuyết tương ứng (Hypotheses - H):

  1. RQ1: Làm thế nào để thiết lập hệ phương trình ổn định và biểu thức hiển giải tích xác định tải trọng tới hạn tuyến tính, phi tuyến cho vỏ nón cụt FGM có gân gia cường khi tính đến biến dạng trượt ngang và sự thay đổi khoảng cách gân?
    • H1: Việc áp dụng lý thuyết biến dạng trượt bậc nhất (First-order Shear Deformation Theory - FSDT) kết hợp kỹ thuật san đều gân Lekhnitskii cải tiến và phương pháp Bubnov-Galerkin sẽ cho phép xây dựng nghiệm giải tích tường minh, giảm sai số cơ học từ 12-25% so với lý thuyết cổ điển CST đối với vỏ dày vừa ($R/h < 50$).
  2. RQ2: Quy luật tương tác đa trường (cơ - nhiệt - điện) và tác động của nền đàn hồi Pasternak ảnh hưởng như thế nào đến tần số dao động tự do và đáp ứng động lực phi tuyến của vỏ và panel nón cụt FGM?
    • H2: Sự tích hợp lớp áp điện với thế hiệu điều khiển $V_a$ cùng độ cứng trượt $K_2$ của nền Pasternak sẽ tạo ra hiệu ứng bù ứng suất màng, giúp triệt tiêu biên độ dao động phi tuyến và dịch chuyển tần số dao động riêng lên vùng an toàn.

Khung lý thuyết của công trình tích hợp lý thuyết biến dạng trượt Timoshenko-Mindlin, quan hệ biến dạng phi tuyến Von Kármán, kỹ thuật san đều gân Lekhnitskii, mô hình nền Pasternak và định luật Hooke nhiệt đàn hồi tổng quát. Phạm vi nghiên cứu bao gồm vỏ nón cụt FGM, panel nón cụt FGM, vỏ sandwich FGM 3 lớp với các quy luật phân bố thể tích: hàm lũy thừa (P-FGM), hàm sigmoid (S-FGM) và hàm mũ (E-FGM), khảo sát toàn diện từ góc bán đỉnh $\beta = 0^\circ \div 75^\circ$, tỷ số bán kính trên bề dày $R/h = 10 \div 200$, dưới các chế độ tải nén dọc trục $P$, áp lực ngoài $q$, trường nhiệt tăng đều $\Delta T$ và tăng tuyến tính $\Delta T(z)$.

Literature Review và Positioning

Nghiên cứu về cơ học kết cấu FGM đã trải qua các bước phát triển quan trọng với nhiều trường phái học thuật:

  1. Dòng nghiên cứu giải tích cho vỏ nón FGM không gia cường: Đứng đầu là Sofiyev cùng các cộng sự (2007, 2011, 2014) với các nghiên cứu mẫu mực về mất ổn định tĩnh, ổn định động và dao động phi tuyến của vỏ nón FGM chịu tải nén dọc trục và áp lực thủy tĩnh. Tuy nhiên, các mô hình của Sofiyev chủ yếu tiếp cận theo hàm ứng suất, áp dụng lý thuyết cổ điển CST và không xét đến hệ gân gia cường phức hợp.
  2. Dòng nghiên cứu số và bán giải tích: Nhóm tác giả Eslami, Naj và cộng sự (2008) đã sử dụng phương pháp vi phân cầu phương (Differential Quadrature Method - DQM) để khảo sát mất ổn định nhiệt của vỏ nón FGM dựa trên phương trình Donnell. Zhao và Liew (2011) ứng dụng phương pháp phần tử tự do kp-Ritz phân tích vỏ nón FGM; Tornabene (2012) áp dụng DQM cho bài toán dao động của vỏ tròn xoay. Điểm hạn chế của phương pháp số là không đưa ra được các công thức đóng (closed-form solutions) hiển thị rõ mối quan hệ định lượng giữa các thông số hình học và vật lý.
  3. Dòng nghiên cứu về kết cấu FGM có gân gia cường tại Việt Nam: Đào Văn Dũng và cộng sự (2013, 2014) đã tiên phong sử dụng kỹ thuật san đều gân Lekhnitskii phân tích vỏ nón cụt FGM chịu tải cơ và tải nhiệt; Nguyễn Đình Đức và cộng sự (2015, 2017) đã mở rộng nghiên cứu ổn định tĩnh và động lực học của vỏ và panel nón FGM trên nền đàn hồi. Dẫu vậy, hầu hết các công trình này đều bị đóng khung trong giả thuyết Kirchhoff-Love (CST).

Trong bức tranh tổng quan đó, tồn tại một tranh luận khoa học cốt lõi (academic debate): Một số quan điểm cho rằng với vỏ nón thoải hoặc vỏ mỏng ($R/h > 100$), lý thuyết CST là hoàn toàn đủ để mô tả ứng xử mất ổn định nhằm tối ưu hóa khối lượng tính toán. Ngược lại, trường phái biến dạng trượt bậc cao khẳng định rằng do tính chất không đồng nhất qua chiều dày của vật liệu FGM, sự tập trung ứng suất tiếp và độ trượt giữa các thớ vật liệu gốm - kim loại tạo ra độ mềm trượt lớn, khiến lý thuyết CST đánh giá sai lệch đáng kể khả năng chịu tải thực tế.

Luận án của Đỗ Quang Chấn định vị chính xác tại giao điểm này: phát triển hệ phương trình giải tích tường minh trên cơ sở FSDT cho vỏ nón cụt FGM có gân gia cường trực giao trên nền Pasternak. Luận án đối chiếu trực tiếp với hai nghiên cứu quốc tế điển hình:

  • So sánh với kết quả của Naj cùng cộng sự (2008) cho vỏ đẳng hướng và FGM không gia cường chịu tải nhiệt, kết quả của luận án cho độ lệch dưới 1.5%, đồng thời mở rộng thêm bài toán có gân gia cường và nền Pasternak mà Naj chưa giải quyết.
  • So sánh với mô hình của Baruch cùng cộng sự (1976) và Sofiyev (2011) về tải trọng tới hạn nén dọc trục $P_{cr}$, kết quả FSDT của luận án giải thích rõ hiện tượng suy giảm tải tới hạn khi góc bán đỉnh $\beta$ tăng và độ dày $h$ tăng, điều mà các công thức CST trước đây đã dự báo quá mức từ 8.4% đến 18.2%.

Đó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 căn bản lý thuyết biến dạng trượt bậc nhất (FSDT Timoshenko-Mindlin) và nguyên lý mở rộng Lekhnitskii vào không gian hình học phi phẳng của vỏ nón cụt FGM. Các đóng góp mang tính nền tảng bao gồm:

  • Mở rộng lý thuyết động học Mindlin cho vỏ tròn xoay FGM có gân: Trường chuyển vị được thiết lập chính xác theo $u_x = u + z\phi_x$, $u_\theta = v + z\phi_\theta$, $u_z = w$, trong đó $\phi_x, \phi_\theta$ là các góc xoay độc lập của pháp tuyến mặt giữa. Hệ số hiệu chỉnh cắt được xác định chuẩn xác $\kappa_x^2 = \kappa_\theta^2 = 5/6$.
  • Mô hình hóa toán học sự biến thiên mật độ gân dọc: Khác với vỏ trụ, khoảng cách giữa các gân dọc trên vỏ nón biến đổi liên tục theo tọa độ đường sinh: $d_s(x) = \frac{2\pi x \sin\beta}{n_s}$. Luận án đã tích hợp sự biến thiên này vào các độ cứng kéo $A_{ij}$, độ cứng uốn - nén liên hợp $B_{ij}$, và độ cứng uốn $D_{ij}$, chuyển hệ phương trình vi phân cấu trúc từ dạng cục bộ sang mô hình giải tích toàn cục tương đương.
  • Thiết lập mô hình giải tích trạng thái màng (membrane state): Trước khi xảy ra hiện tượng mất ổn định nhánh (bifurcation buckling), trường lực màng trước tới hạn được xác định chuẩn xác dưới dạng $N_{x0} = -\frac{P}{2\pi x \sin\beta \cos\beta} = -\frac{p x_0}{x \cos\beta}$, $N_{\theta 0} = 0$, $N_{x\theta 0} = 0$.

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

Khung phân tích của luận án là sự hợp nhất đồng bộ giữa ba lý thuyết cơ học lớn: (1) Lý thuyết phi tuyến Von Kármán mô tả chuyển vị uốn lớn; (2) Kỹ thuật đồng nhất hóa giải tích Lekhnitskii cho kết cấu gia cường; (3) Lý thuyết tương tác đa trường nhiệt - áp điện.

       ┌─────────────────────────────────────────────────────────────┐
       │                Vật liệu FGM (P-, S-, E-FGM)                 │
       │     E(z), α(z), K(z), ρ(z) - Tương quan Gốm/Kim loại        │
       └──────────────────────────────┬──────────────────────────────┘
                                      ▼
       ┌─────────────────────────────────────────────────────────────┐
       │            Động học FSDT & San đều gân Lekhnitskii          │
       │    ds(x) = 2πx sinβ / ns;  κ² = 5/6; Quan hệ Trượt Ngang    │
       └──────────────────────────────┬──────────────────────────────┘
                                      ▼
       ┌─────────────────────────────────────────────────────────────┐
       │             Hệ phương trình cân bằng phi tuyến              │
       │   5 PDE đạo hàm riêng có hệ số là hàm tọa độ: Rij(u,v,w,φ)  │
       └──────────────────────────────┬──────────────────────────────┘
                                      ▼
       ┌─────────────────────────────────────────────────────────────┐
       │         Biến đổi khử biến x & Phép chiếu Galerkin           │
       │   Nhân x^i để chuyển tích phân không sơ cấp -> đa thức      │
       └──────────────────────────────┬──────────────────────────────┘
                                      ▼
       ┌─────────────────────────────────────────────────────────────┐
       │ Nghiệm giải tích tải tới hạn (Pcr, Tcr) & ODE Động lực học  │
       │   Runge-Kutta bậc 4 -> Đáp ứng thời gian w(t) & Đồ thị Pha  │
       └─────────────────────────────────────────────────────────────┘

Điểm đột phá phương pháp luận là phép biến đổi toán học nhằm xử lý các tích phân phi sơ cấp dạng $\int \frac{1}{x^k} \sin\left(\frac{m\pi(x-x_0)}{L}\right) dx$. Tác giả đã nhân tương đương hệ phương trình vi phân với các thừa số tọa độ $x$ và $x^2$, chuyển đổi các biểu thức chứa mẫu số biến thiên thành các tích phân Galerkin chuẩn trên toàn miền $[x_0, x_0+L] \times [0, 2\pi]$, tạo tiền đề thiết lập định thức suy biến $5 \times 5$: $$\det |X_{ij}(P_{cr})| = 0 \quad (i, j = 1 \div 5)$$ cho phép suy ra biểu thức đóng hiển (closed-form expression) của tải trọng tới hạn $P_{cr}$ và độ biến thiên nhiệt tới hạn $T_{cr}$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng suy diễn chuẩn xác (deductive positivism), sử dụng phương pháp giải tích giải toán cơ học vật rắn. Cấu trúc mô hình hóa bao quát các cấp độ phân tích:

  • Cấp độ vi mô/thành phần: Quy tắc hỗn hợp (Voigt scheme) hiệu chỉnh theo hàm lũy thừa: $$P_{eff}(z) = (P_c - P_m)\left(\frac{2z+h}{2h}\right)^k + P_m$$ trong đó $k \ge 0$ là chỉ số tỷ phần thể tích gốm.
  • Cấp độ trung mô/cấu trúc lớp: Vỏ nón FGM đơn lớp có bề dày $h$, panel nón tích hợp lớp áp điện có bề dày $h_p$ chịu điện áp ngoài $V_a$, và vỏ sandwich 3 lớp với lõi FGM hoặc lõi gốm/kim loại đồng nhất.
  • Cấp độ vĩ mô/kết cấu - môi trường: Vỏ liên kết tựa đơn (simply supported edge conditions) đặt trên môi trường đàn hồi liên tục Pasternak với phản lực nền: $$q_f = K_1 w - K_2 \left(\frac{\partial^2 w}{\partial x^2} + \frac{1}{x}\frac{\partial w}{\partial x} + \frac{1}{x^2 \sin^2\beta}\frac{\partial^2 w}{\partial\theta^2}\right)$$

Quy trình nghiên cứu rigorous

Quy trình mô hình hóa và kiểm chứng tuân thủ các bước nghiêm ngặt:

  1. Thiết lập quan hệ động học và vật lý: Xây dựng phương trình trạng thái màng và phương trình biến dạng phi tuyến Von Kármán.
  2. Thiết lập hệ phương trình vi phân biến dạng: Hệ gồm 5 phương trình đạo hàm riêng liên kết phi tuyến giữa 5 hàm ẩn: chuyển vị dọc $u_1(x,\theta)$, chuyển vị vòng $v_1(x,\theta)$, độ võng hướng tâm $w_1(x,\theta)$, và hai góc xoay $\phi_{x1}(x,\theta), \phi_{\theta 1}(x,\theta)$.
  3. Lựa chọn hàm dạng nghiệm thỏa mãn điều kiện biên: $$u_1 = U \cos\left(\frac{m\pi(x-x_0)}{L}\right)\sin(n\theta), \quad v_1 = V \sin\left(\frac{m\pi(x-x_0)}{L}\right)\cos(n\theta)$$ $$w_1 = W \sin\left(\frac{m\pi(x-x_0)}{L}\right)\sin(n\theta), \quad \phi_{x1} = \Phi_1 \cos\left(\frac{m\pi(x-x_0)}{L}\right)\sin(n\theta)$$ $$\phi_{\theta 1} = \Phi_2 \frac{1}{x\sin\beta} \sin\left(\frac{m\pi(x-x_0)}{L}\right)\cos(n\theta)$$ trong đó $(m, n)$ lần lượt là số nửa bước sóng theo hướng đường sinh và số bước sóng theo hướng chu vi.
  4. Giải hệ phương trình: Áp dụng thuật toán Bubnov-Galerkin để đưa về phương trình đại số phi tuyến hoặc phương trình vi phân cấp hai theo thời gian. Giải đáp ứng động lực phi tuyến bằng thuật toán Runge-Kutta bậc 4 với bước tích phân thời gian $\Delta t = 10^{-4}\text{ s}$ đảm bảo tiêu chuẩn ổn định Courant-Friedrichs-Lewy.

Data và phân tích

Các bộ thông số vật liệu chuẩn được sử dụng trong tính toán số:

  • Hệ vật liệu Nhôm/Nhôm oxit ($Al/Al_2O_3$): Kim loại Nhôm ($E_m = 70\text{ GPa}$, $\rho_m = 2707\text{ kg/m}^3$, $\alpha_m = 23\times 10^{-6}\text{ K}^{-1}$, $K_m = 204\text{ W/mK}$); Gốm $Al_2O_3$ ($E_c = 380\text{ GPa}$, $\rho_c = 3800\text{ kg/m}^3$, $\alpha_c = 7.4\times 10^{-6}\text{ K}^{-1}$, $K_c = 10.4\text{ W/mK}$).
  • Hệ vật liệu Titan/Thép không gỉ và Gốm Silicon Nitride ($Si_3N_4$).
  • Tham số hình học: Tỷ số bán kính $R/h = 20 \div 200$, tỷ số chiều dài $L/R = 1 \div 5$, góc bán đỉnh $\beta = 10^\circ \div 75^\circ$.

Toàn bộ thuật toán giải tích và tính toán số được lập trình tự động trên nền tảng phần mềm toán học chuyên dụng Maple và Matlab. Độ hội tụ và tính ổn định thuật toán được kiểm tra lặp qua các dải mode sóng $(m, n)$ từ $(1, 1)$ đến $(10, 20)$ để trích xuất mode mất ổn định nguy hiểm nhất tương ứng với giá trị $P_{cr}$ hoặc $T_{cr}$ cực tiểu.

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 đột phá có giá trị thực nghiệm và định lượng cao:

  1. Hiệu ứng gia cường phi đối xứng của hệ gân trực giao: Bố trí gân ngoài (outside stiffeners) mang lại hiệu quả kháng uốn và tăng tải trọng tới hạn $P_{cr}$ vượt trội so với bố trí gân trong (inside stiffeners). Cụ thể, khi bố trí hệ gân trực giao $n_s = n_r = 15$ ở mặt ngoài, tải nén tới hạn tăng từ 1.85 đến 3.42 lần so với vỏ không gia cường. Luận án chỉ ra rằng gân đặt ở bề mặt giàu gốm (mặt chịu nén chính) triệt tiêu mômen uốn lệch tâm tốt hơn đáng kể.
  2. Sai số giới hạn của lý thuyết vỏ cổ điển (CST): Khi tỷ số bán kính trên bề dày $R/h \le 30$, lý thuyết CST đánh giá quá cao tải trọng tới hạn $P_{cr}$ từ 14.6% đến 22.8% so với lý thuyết biến dạng trượt FSDT. Hiện tượng này khẳng định tính bắt buộc của việc sử dụng FSDT trong thiết kế vỏ nón FGM chịu tải nặng.
  3. Quy luật chuyển dịch mode mất ổn định nhiệt: Dưới trường nhiệt tăng đều, vỏ mất ổn định ở số bước sóng vòng $n$ nhỏ hơn ($n = 3 \div 5$), trong khi dưới trường nhiệt tăng tuyến tính qua chiều dày ($\Delta T(z)$), ứng suất tập trung mạnh tại mặt gốm làm số bước sóng vòng tới hạn tăng vọt ($n = 8 \div 14$). Độ biến thiên nhiệt tới hạn $T_{cr}$ của vỏ có gân cao hơn vỏ không gân từ 45% đến 110%.
  4. Hiệu ứng triệt tiêu động lực của lớp áp điện chủ động: Đối với panel nón cụt FGM tích hợp vật liệu áp điện, việc cấp điện áp âm ($V_a < 0$) tạo ra trường ứng suất kéo trước, làm tăng tần số dao động tự do riêng $\omega$ lên tới 18.5%, đồng thời dập tắt biên độ dao động phi tuyến cưỡng bức theo thời gian nhanh hơn 40% so với trường hợp không có điều khiển chủ động.
  5. Vai trò kiềm chế của hệ số cắt nền Pasternak ($K_2$): Nền Pasternak tác động phi tuyến tích cực lên kết cấu; tham số cắt $K_2$ đóng góp tới 65% vào việc gia tăng tải tới hạn sau mất ổn định (postbuckling load capacity), vượt trội hoàn toàn so với mô hình nền Winkler đơn thuần ($K_1$).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung giải tích tổng quát FSDT cho lớp bài toán vỏ nón FGM có gân, đóng góp vào lý thuyết kết cấu phi tuyến đa trường.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình biến đổi tích phân Galerkin cho các hệ phương trình vi phân có hệ số phụ thuộc không gian, có thể chuyển giao áp dụng trực tiếp cho các hình học phức tạp khác như vỏ parabolic, hyperboloid hoặc cánh tuabin biến diện tích.
  • Về mặt kỹ thuật ứng dụng: Cung cấp các biểu đồ thiết kế (design charts) tải trọng tới hạn cho kỹ sư hàng không vũ trụ chế tạo vỏ tên lửa, đầu dẫn tên lửa vượt âm, thân tàu ngầm mini và bồn chứa áp lực nhiệt hạch.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu khách quan:

  • Giới hạn điều kiện biên: Nghiệm giải tích hiện tại tập trung cho điều kiện biên tựa đơn bốn cạnh (simply supported edges: $v_1 = w_1 = \phi_{\theta 1} = 0, N_{x1} = M_{x1} = 0$ tại $x = x_0, x_0 + L$). Các liên kết ngàm cứng hoàn toàn (clamped) hoặc biên tự do chưa thể giải hoàn toàn bằng biểu thức đóng do sự phức tạp của hàm dạng sóng.
  • Mô hình vật liệu: Chưa xét đến tính dị hướng tế vi của các hạt FGM định hướng trong quá trình in 3D (additive manufacturing) và ảnh hưởng của độ rỗng ngẫu nhiên (porosity) phát sinh do thiêu kết nhiệt.
  • Môi trường tải trọng: Tải nhiệt giả định ở trạng thái tĩnh hoặc giả tĩnh; chưa phân tích hiện tượng sốc nhiệt phi Fourier (hyperbolic heat conduction) sinh ra do xung nhiệt siêu ngắn.

Chương trình nghiên cứu tương lai cần phát triển 4 hướng:

  1. Mở rộng lời giải sang lý thuyết biến dạng trượt bậc cao (HSDT - Reddy) và lý thuyết biến dạng cắt vô sai tuyến tính (Sinusoidal Shear Deformation Theory).
  2. Phát triển mô hình phân tích vỏ nón FGM xốp (Porous FGM) có gân gia cường với cấu trúc rỗng chức năng đa chiều.
  3. Nghiên cứu thực nghiệm phá hủy và mất ổn định động trên mô hình in 3D FGM nón cụt để đối chuẩn trực tiếp với nghiệm giải tích.
  4. Mở rộng bài toán sang vỏ nón FGM gia cường bằng ống nano carbon (FG-CNTRC) hoặc graphene platelets (GPLs).

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

Công trình mang lại tác động học thuật và công nghệ sâu rộng:

  • Tác động học thuật: Dự báo tạo lập nguồn trích dẫn học thuật vững chắc trong phân nhóm tạp chí Cơ học và Kết cấu Quốc tế (Q1: Thin-Walled Structures, Composite Structures, Aerospace Science and Technology). Phương pháp tiếp cận giải tích của luận án là tài liệu quy chuẩn cho các nghiên cứu sinh chuyên ngành Cơ học tính toán và Cơ học vật rắn.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp vào Viện Công nghệ Hàng không Vũ trụ, Tập đoàn Công nghiệp - Viễn thông Quân đội (Viettel), và các tổ hợp đóng tàu ngầm/tàu chiến hiện đại nhằm tối ưu hóa tỷ số độ bền/khối lượng của kết cấu thân vỏ chịu nén ngoài và bức xạ nhiệt.
  • Lợi ích an ninh - xã hội: Tối ưu hóa việc sử dụng vật liệu hợp kim quý và gốm kỹ thuật cao, giảm tiêu hao nhiên liệu đẩy trong công nghệ phóng và nâng cao độ an toàn vận hành của lò phản ứng năng lượng nguyên tử quốc gia.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học: Khai thác các biểu thức hiển giải tích tường minh, hệ ma trận độ cứng $X_{ij}$, và quy trình giải tích số trên Maple/Matlab làm học liệu và cơ sở phát triển đề tài chuyên sâu.
  • Kỹ sư Thiết kế Hàng không & Quốc phòng: Sử dụng các đường cong tải - độ võng phi tuyến ($P - W/h$), biểu đồ tần số - biên độ để tính toán kiểm tra độ bền mỏi, độ ổn định vỏ nón trong buồng đốt tên lửa và chóp nón đầu đạn khí động.
  • Các nhà phát triển công nghệ vật liệu mới: Nhận diện được tỷ phần tối ưu $k$ và quy luật phân bố thể tích phù hợp để định hướng quy trình luyện kim bột và chế tạo đắp lớp (in 3D tự do).

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 hệ phương trình cân bằng và xác định tải trọng tới hạn dạng giải tích hiển (closed-form analytical solution) cho vỏ nón cụt FGM có gân gia cường dựa trên lý thuyết biến dạng trượt bậc nhất FSDT, đồng thời xử lý triệt để bài toán mật độ gân dọc biến thiên $d_s(x) = \frac{2\pi x\sin\beta}{n_s}$ mà chưa công trình giải tích nào trước đó giải quyết được.

2. Đột phá phương pháp luận so với các nghiên cứu trước đây? So với phương pháp CST của Đào Văn Dũng (2013) và phương pháp hàm ứng suất của Sofiyev (2011), luận án phát triển cách tiếp cận theo hàm chuyển vị đầy đủ 5 thành phần ($u_1, v_1, w_1, \phi_{x1}, \phi_{\theta 1}$). Đột phá toán học nằm ở việc thực hiện phép nhân khử biến $x^i$ biến đổi các tích phân suy biến không sơ cấp thành các tích phân Galerkin chuẩn, cho phép thu được nghiệm tường minh thay vì phải xấp xỉ số hoàn toàn.

3. Phát hiện bất ngờ nhất từ dữ liệu tính toán? Hiện tượng "đảo ngược tính kháng uốn" khi thay đổi trường nhiệt: Trong điều kiện nhiệt tăng đều, gân gia cường dọc đường sinh đóng vai trò chủ đạo trong việc nâng cao tải $P_{cr}$; nhưng khi chịu tải nhiệt tăng tuyến tính qua chiều dày với độ chênh lớn giữa hai mặt ($\Delta T > 300\text{ K}$), hệ gân vòng (rings) lại thể hiện vai trò kiềm chế biến dạng nhiệt mặt ngoài hiệu quả gấp 1.7 lần so với gân dọc.

4. Khả năng tái lập quy trình nghiên cứu (Replication Protocol)? Quy trình tính toán hoàn toàn minh bạch: Toàn bộ các hệ số giải tích $A_{ij}, B_{ij}, D_{ij}$, hệ số ma trận $R_{ij}$ và $X_{ij}$ đều được công bố chi tiết trong phụ lục toán học. Bất kỳ nhà nghiên cứu nào cũng có thể tái lập mã nguồn trên Matlab/Maple để cho ra kết quả trùng khớp $100%$ với các bảng số liệu đã công bố.

5. Lộ trình phát triển 10 năm tiếp theo của hướng nghiên cứu này? Lộ trình hướng tới: (1) Tích hợp học máy (Machine Learning/Physics-Informed Neural Networks - PINNs) để tối ưu hóa hình học gân topo phi tuyến; (2) Chế tạo thực nghiệm vật liệu FGM in 3D dạng composite thông minh 4D có khả năng tự thích ứng nhiệt - điện; (3) Phân tích động lực học va chạm vận tốc siêu cao (hypervelocity impact) cho các cấu trúc bay siêu vượt âm (Hypersonic Gliders).

Kết luận

  1. Luận án đã xây dựng thành công mô hình giải tích toàn diện phân tích ổn định tuyến tính, phi tuyến tĩnh và đáp ứng động lực của vỏ và panel nón cụt FGM có gân gia cường trên nền đàn hồi Pasternak chịu tải cơ - nhiệt - điện phức hợp.
  2. Thiết lập được hệ biểu thức hiển giải tích xác định chính xác tải nén tới hạn $P_{cr}$, tải nhiệt tới hạn $T_{cr}$, tần số dao động tự do $\omega$, và các đường cong đáp ứng thời gian uốn phi tuyến $W(t)$.
  3. Khẳng định tính ưu việt bắt buộc của lý thuyết FSDT so với lý thuyết cổ điển CST trong thiết kế vỏ FGM dày vừa, hạn chế rủi ro phá hủy kết cấu do đánh giá quá cao khả năng chịu lực.
  4. Chứng minh vai trò kiểm soát chủ động của lớp áp điện và tác động kháng uốn vượt trội của bố trí gân ngoài kết hợp nền đàn hồi Pasternak.
  5. Mở ra 3 nhánh nghiên cứu mới về vật liệu chức năng biến thiên: vỏ FGM có cấu trúc rỗng đa thang, tối ưu hóa topo gân gia cường thông minh, và cơ học tính toán đa trường cho các thiết bị bay siêu thanh thế hệ mới.
  6. Kết quả nghiên cứu tạo nền tảng lý thuyết vững chắc, làm phong phú kho tàng Cơ học vật rắn Việt Nam và đóng góp trực tiếp vào công nghệ thiết kế, chế tạo kết cấu hàng không vũ trụ hiện đại.