Tổng quan về luận án

Nghiên cứu về độ bền, ổn định và đáp ứng động lực học của kết cấu vỏ làm từ vật liệu có cơ tính biến đổi (Functionally Graded Materials - FGM) là một trong những hướng đi tiên phong của Cơ học vật rắn biến dạng hiện đại. Ra đời lần đầu tại Sendai, Nhật Bản vào năm 1984, FGM được tạo thành từ hỗn hợp gốm (ceramic) và kim loại (metal) với tỷ phần thể tích biến đổi liên tục qua chiều dày, giúp triệt tiêu tập trung ứng suất và khắc phục hiện tượng bong tách phân lớp vốn tồn tại cố hữu trong vật liệu composite truyền thống. Trong khi các kết cấu dạng tấm và vỏ trụ FGM đã được khảo sát tương đối phong phú, bài toán vỏ tròn xoay mà điển hình là vỏ nón cụt FGM có gân gia cường đặt trên nền đàn hồi vẫn tồn tại rào cản toán học phức tạp: hệ phương trình cân bằng và chuyển động là hệ phương trình đạo hàm riêng phi tuyến có hệ số biến thiên theo tọa độ $x$. Luận án tiến sĩ của nghiên cứu sinh Đỗ Quang Chấn với đề tài "Phân tích ổn định và đáp ứng động lực của vỏ nón cụt FGM" do GS.TSKH. Nguyễn Đình Đức và PGS.TS. Vũ Đỗ Long hướng dẫn tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2019) đã giải quyết trọn vẹn bài toán giải tích này.

Khoảng trống nghiên cứu (research gap) trọng yếu được xác định: Các nghiên cứu quốc tế đương thời của Sofiyev (2010–2014) hay Naj và Eslami (2008) chủ yếu tập trung vào vỏ nón cụt FGM không gân, bỏ qua biến dạng trượt ngang (chỉ dùng lý thuyết vỏ cổ điển - CST), hoặc mới dừng ở phương pháp số xấp xỉ mà thiếu vắng các nghiệm giải tích tường minh cho kết cấu vỏ nón cụt FGM gia cường gân trực giao lệch tâm chịu tác dụng đồng thời của trường cơ - nhiệt - điện và tương tác nền đàn hồi.

Luận án thiết lập 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học cốt lõi:

  • RQ1: Làm thế nào để thiết lập mô hình giải tích chính xác xác định tải trọng tới hạn tĩnh ($P_{cr}$, $T_{cr}$) của 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)?
    • H1: Mô hình FSDT kết hợp kỹ thuật san đều gân Lekhnitskii cải tiến sẽ dự báo tải tới hạn chính xác hơn CST từ 8.5% đến 22.4% đối với các vỏ có độ dày trung bình ($R/h \le 30$).
  • RQ2: Ứng xử sau mất ổn định (postbuckling) phi tuyến của vỏ nón cụt FGM gia cường dưới tải nén dọc trục và tải nhiệt biến thiên chịu ảnh hưởng như thế nào bởi khoảng cách gân thay đổi $d_s(x)$?
    • H2: Sự biến thiên khoảng cách gân dọc đường sinh theo tọa độ $x$ làm thay đổi đáng kể độ cứng cục bộ, và đường cong tải - độ võng phi tuyến ($P - W/h$) có dạng cực trị phụ thuộc vào chỉ số tỷ phần thể tích $k$.
  • RQ3: Tác động ghép kênh cơ - nhiệt - điện lên tần số dao động tự do và đáp ứng động lực phi tuyến của panel nón cụt FGM áp điện trên nền Pasternak diễn ra theo quy luật nào?
    • H3: Điện áp ngoài $V_a$ và độ cứng lớp cắt $K_2$ của nền đàn hồi có khả năng kiểm soát chủ động biên độ dao động phi tuyến và dải tần số cộng hưởng của panel.
  • RQ4: Ảnh hưởng của 5 cấu trúc vỏ sandwich nón cụt FGM (lõi kim loại, lõi gốm, lớp mặt FGM) đến khả năng chịu nhiệt tới hạn $T_{cr}$?
    • H4: Cấu trúc sandwich lõi gốm với lớp mặt FGM đem lại khả năng chịu nhiệt bức xạ cực hạn cao nhất so với các mô hình vỏ FGM đơn lớp thông thường.

Khung lý thuyết của luận án tích hợp: Lý thuyết biến dạng trượt bậc nhất Mindlin-Timoshenko (FSDT), quan hệ phi tuyến hình học Von Kármán, quan điểm san đều tác dụng gân của Lekhnitskii, mô hình nền hai tham số Pasternak ($K_1, K_2$) và tiêu chuẩn cân bằng lân cận Euler. Phạm vi khảo sát trải rộng trên nhiều cấu hình hình học (góc bán đỉnh $\beta = 10^\circ \div 75^\circ$, tỷ số $R/h = 20 \div 100$, $L/R = 1 \div 5$), các tổ hợp vật liệu FGM điển hình ($\text{Al/Al}_2\text{O}_3$, $\text{Ti-6Al-4V/ZrO}_2$, Thép không gỉ/Silicon nitride) với chỉ số $k = 0 \div \infty$.

Literature Review và Positioning

Nghiên cứu về cơ học kết cấu FGM chia thành 3 dòng tư tưởng chính:

  1. Dòng nghiên cứu giải tích trên tấm và vỏ mỏng cổ điển: Tiêu biểu là các công trình của Reddy (2000), Shen (2009), Bích và cộng sự (2012) [174], Dũng và cộng sự (2013, 2014) [175-177]. Dòng này chủ yếu sử dụng lý thuyết vỏ cổ điển (CST), bỏ qua góc xoay trượt ngang $\gamma_{xz}, \gamma_{\theta z} = 0$, dẫn đến việc đánh giá quá cao (overestimate) độ cứng kết cấu đối với vỏ dày và dày vừa.
  2. Dòng nghiên cứu bán giải tích và số xấp xỉ: Đại diện bởi Sofiyev và cộng sự (2007–2016) [137–155], Tornabene (2011) [159], Viola (2013) [161], Jam và Kiani (2015) [172]. Nhóm này áp dụng phương pháp vi phân cầu phương (DQM), phương pháp không lưới (kp-Ritz) hoặc phần tử hữu hạn (FEM). Mặc dù giải quyết được tính bất đối xứng, cách tiếp cận này phụ thuộc dung sai lưới và không cung cấp được biểu thức hiển tường minh giải tích để khảo sát trực tiếp liên hệ tham số.
  3. Dòng nghiên cứu ổn định nhiệt - cơ vỏ nón FGM của Eslami và nhóm tác giả Iran: Eslami, Naj, Baruch (2008) [156, 187] khảo sát vỏ nón cụt FGM dưới tải nhiệt tăng đều và tuyến tính nhưng hoàn toàn bỏ qua hệ gân gia cường và tác động đàn hồi của môi trường bao quanh.

Tranh luận học thuật cốt lõi diễn ra giữa quan điểm của Sofiyev [146] (cho rằng phương pháp đặt bài toán theo hàm ứng suất kết hợp Galerkin đơn mode là tối ưu cho vỏ nón) và trường phái của Nguyễn Đình Đức [181-184] (khẳng định phương pháp đặt theo hệ 5 chuyển vị $u, v, w, \phi_x, \phi_\theta$ theo FSDT mới mô tả chính xác ứng xử tiếp xúc gân - vỏ và trượt cắt ngang). Luận án của Đỗ Quang Chấn định vị chính xác tại điểm giao thoa này: kế thừa tính linh hoạt của mô hình dịch chuyển 5 thành phần, đồng thời giải quyết triệt để bài toán hệ phương trình vi phân biến thiên phức tạp bằng thuật toán tích phân trọng số nhân tử $x^i$ độc đáo.

So sánh đối chuẩn với 2 nghiên cứu quốc tế kinh điển:

  • So sánh 1 (với Naj et al. [156] và Baruch et al. [187]): Khi suy biến mô hình luận án về vỏ đẳng hướng thuần nhất không gân chịu nén dọc trục, sai số tương đối của tải tới hạn $P_{cr}$ không vượt quá 0.85%, khẳng định tính chuẩn xác tuyệt đối của phương trình trạng thái màng tiền tới hạn $N_x^0 = -P / (2\pi x \sin\beta \cos\beta)$.
  • So sánh 2 (với Sofiyev [146]): Khi so sánh quan hệ biên độ - tần số dao động phi tuyến của vỏ nón FGM, quỹ đạo pha và đường cong suy giảm tần số phi tuyến của luận án hoàn toàn đồng nhất với công bố quốc tế của Sofiyev, đồng thời mở rộng thêm ảnh hưởng của hệ gân gia cường $n_s, n_r$ và hệ số nền cắt $K_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 lý thuyết tấm vỏ phi tuyến Mindlin-Timoshenko và nguyên lý năng lượng cực tiểu cho lớp kết cấu vỏ nón cụt FGM có độ cong biến thiên:

  • Xây dựng hệ phương trình cân bằng phi tuyến 5 ẩn chuyển vị: Thiết lập hệ 5 phương trình đạo hàm riêng phi tuyến bậc cao mô tả đầy đủ các hiệu ứng trượt cắt ngang ($\phi_x, \phi_\theta$), hiệu ứng xoắn màng ($k_{x\theta}$), và tương tác nén - uốn - cắt kết hợp.
  • Mô hình hóa chính xác hàm phân bố vật liệu: Định thức hóa 3 quy luật biến thiên cơ tính: Power-law (P-FGM), Sigmoid (S-FGM) và Exponential (E-FGM), trong đó mô đun đàn hồi $E(z)$, hệ số giãn nở nhiệt $\alpha(z)$ và hệ số dẫn nhiệt $K(z)$ biến thiên trơn theo tọa độ $z$: $$E(z) = E_m + (E_c - E_m)\left(\frac{2z+h}{2h}\right)^k$$
  • Thiết lập tiêu chuẩn mất ổn định Euler - Galerkin: Chuyển đổi thành công hệ phương trình phi tuyến về hệ đại số thuần nhất $5 \times 5$ thông qua nghiệm dạng chuỗi điều hòa thỏa mãn biên tựa đơn, từ đó thiết lập biểu thức hiển xác định định thức đặc trưng triệt tiêu $\det[X_{ij}(P_{cr})] = 0$.

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

Khung phân tích của luận án kết hợp đồng thời 4 trụ cột cơ học:

  1. Lý thuyết biến dạng trượt bậc nhất (FSDT) với hệ số hiệu chỉnh cắt $\kappa_x^2 = \kappa_\theta^2 = 5/6$.
  2. Kỹ thuật san đều gân Lekhnitskii mở rộng: Coi hệ gân dọc ($n_s$) và gân vòng ($n_r$) phân bố liên tục trên bề mặt, với khoảng cách gân dọc biến thiên $d_s(x) = (2\pi x \sin\beta)/n_s$, phản ánh trung thực bản chất hình học nón cụt.
  3. Mô hình tương tác nền hai tham số Pasternak: Phản lực nền $q_f = K_1 w - K_2 \nabla^2 w$, mô phỏng cả sức kháng lún nén Winkler ($K_1$) và sức kháng trượt tiếp tuyến ($K_2$).
  4. Hiệu ứng áp điện nghịch đảo: Tích hợp các số hạng lực điện trường $E_z = V_a / h_p$ vào ma trận ứng suất - biến dạng tổng quát.

Biên điều kiện áp dụng (boundary conditions): Kết cấu tựa đơn bốn cạnh ($S-S-S-S$), biến dạng nhỏ và độ võng nhỏ vừa thỏa mãn giả thiết von Kármán.

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 định lượng (Positivism) kết hợp giải tích toán - cơ suy diễn chặt chẽ. Hệ phương trình vi sai được xây dựng từ phương trình cân bằng phần tử vỏ, định luật Hooke tổng quát cho vật liệu dị hướng cục bộ và các hệ thức hình học biến dạng - chuyển vị.

Quy trình nghiên cứu rigorous

Quy trình tính toán giải tích bao gồm 4 giai đoạn chuẩn hóa:

  1. Giai đoạn 1 - Xác định trạng thái màng tiền tới hạn: Giải phương trình cân bằng phi mômen để tìm nghiệm chính xác của trường ứng suất màng $N_x^0, N_\theta^0, N_{x\theta}^0$ dưới tải trọng cơ học và gradient nhiệt độ.
  2. Giai đoạn 2 - Thiết lập phương trình gia số ổn định: Dựa trên tiêu chuẩn cân bằng lân cận, đặt $u = u_0 + u_1, \dots, w = w_0 + w_1$ để phân tách trạng thái cân bằng dừng và trạng thái nhiễu loạn lân cận.
  3. Giai đoạn 3 - Xử lý điểm kỳ dị toán học bằng tích phân Galerkin nhân tử trọng số: Để khắc phục triệt để biến số $x$ ở mẫu số vốn dẫn tới các tích phân phi sơ cấp dạng $\int \frac{\sin(m\pi x/L)}{x^p} dx$, luận án nhân hai vế phương trình với $x^1$ và $x^2$ trước khi áp dụng trực giao hóa Bubnov-Galerkin trên toàn miền tọa độ $[x_0, x_0 + L] \times [0, 2\pi]$.
  4. Giai đoạn 4 - Giải động lực học phi tuyến bằng Runge-Kutta bậc 4: Đưa phương trình vi phân chuyển động phi tuyến theo thời gian về hệ phương trình vi phân thường bậc nhất và tích phân số theo bước thời gian cực nhỏ $\Delta t = 10^{-5}\text{ s}$ trên môi trường Matlab.
[Mô hình hình học & Vật liệu FGM] 
[Hệ thức biến dạng FSDT & San đều gân Lekhnitskii]
[Hệ 5 Phương trình đạo hàm riêng phi tuyến]
[Nhân tử trọng số x^i]     [Trạng thái màng tiền tới hạn]
[Phương pháp Bubnov-Galerkin]
[Ổn định tĩnh (P_cr, T_cr)] [Động lực phi tuyến Runge-Kutta]

Data và phân tích

Thông số cơ tính vật liệu chuẩn được sử dụng:

  • Nhôm ($\text{Al}$): $E_m = 70.0\text{ GPa}$, $\nu = 0.3$, $\rho_m = 2707\text{ kg/m}^3$, $\alpha_m = 23.0 \times 10^{-6}\text{ K}^{-1}$, $K_m = 204\text{ W/mK}$.
  • Gốm Alumina ($\text{Al}_2\text{O}_3$): $E_c = 380.0\text{ GPa}$, $\nu = 0.3$, $\rho_c = 3800\text{ kg/m}^3$, $\alpha_c = 7.4 \times 10^{-6}\text{ K}^{-1}$, $K_c = 10.4\text{ W/mK}$.
  • Titan hợp kim ($\text{Ti-6Al-4V}$): $E_m = 105.7\text{ GPa}$, $\alpha_m = 7.57 \times 10^{-6}\text{ K}^{-1}$.
  • Zirconia ($\text{ZrO}_2$): $E_c = 168.06\text{ GPa}$, $\alpha_c = 8.95 \times 10^{-6}\text{ K}^{-1}$.

Các công cụ phần mềm tính toán: Xây dựng bộ mã nguồn giải tích tự động trên Maple để rút gọn biểu thức ma trận $X_{ij}$ và lập trình thuật toán số trên Matlab để truy xuất đồ thị đáp ứng động lực $W/h - t$ và đường cong tải - độ võng $P - W/h$.

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

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

  1. Hiệu ứng gia cường phi tuyến của hệ gân trực giao: Sự xuất hiện của gân gia cường dọc ($n_s = 15$) và gân vòng ($n_r = 15$) làm tăng tải nén dọc trục tới hạn $P_{cr}$ lên từ 2.4 đến 3.8 lần so với vỏ không gân tương đương. Đặc biệt, gân đặt ở bề mặt ngoài (rich-ceramic) cho khả năng chịu tải cao hơn 14.6% so với gân đặt ở bề mặt trong do mômen tĩnh đối với mặt trung hòa thay đổi có lợi.
  2. Sự chênh lệch bản chất giữa lý thuyết FSDT và CST: Đối với vỏ nón cụt dày vừa ($R/h = 20$), lý thuyết vỏ cổ điển (CST) đánh giá sai lệch tải tới hạn $P_{cr}$ cao hơn tới 18.7% so với lý thuyết FSDT do bỏ qua góc trượt ngang. Khi góc bán đỉnh $\beta$ tăng từ $15^\circ$ lên $60^\circ$, độ chênh lệch giữa hai lý thuyết càng gia tăng rõ rệt, chứng minh tính bắt buộc phải áp dụng FSDT trong thiết kế thực tế.
  3. Cơ chế chịu nhiệt tới hạn $T_{cr}$ phi trực quan: Khi tải nhiệt tăng tuyến tính qua chiều dày, chênh lệch nhiệt độ tới hạn $\Delta T_{cr}$ của vỏ nón FGM có gân cao gấp 1.85 lần so với trường hợp gia nhiệt tăng đều. Điều này giải thích bởi gradient nhiệt độ tạo ra nội ứng suất uốn ngược chiều với biến dạng phồng tới hạn của vỏ.
  4. Vai trò triệt tiêu rung động của nền đàn hồi Pasternak: Hệ số nền cắt $K_2$ đóng vai trò vượt trội so với hệ số nén Winkler $K_1$. Việc tăng $K_2$ từ $0$ lên $5 \times 10^5\text{ N/m}$ làm tăng tần số dao động riêng cơ bản của panel nón FGM lên tới 42.3%, đồng thời rút ngắn 65% thời gian dập tắt dao động biên độ lớn.
  5. Khả năng điều khiển chủ động của lớp áp điện: Đặt điện áp kích hoạt $V_a > 0$ (ngược chiều phân cực) tạo ra ứng suất kéo trước trong màng, giúp dịch chuyển đường cong đáp ứng động lực $W/h - t$ về miền biên độ an toàn, triệt tiêu hiện tượng cộng hưởng phi tuyến.

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ nghiệm giải tích dạng chuẩn (benchmark closed-form solutions) làm thước đo kiểm chứng cho tất cả các phần mềm phần tử hữu hạn thương mại (Abaqus, Ansys) khi phân tích kết cấu FGM phức tạp.
  • Về mặt phương pháp luận: Kỹ thuật nhân trọng số giải tích để xử lý phương trình vi phân có hệ số hàm số mở ra hướng đi mới giải quyết các kết cấu hình học phức tạp khác như vỏ hyperboloid, vỏ paraboloid hay vỏ trống FGM.
  • Về mặt công nghệ ứng dụng:
    • Hàng không vũ trụ: Ứng dụng trực tiếp vào tối ưu hóa vỏ mũi nón tên lửa đẩy, thân tàu con thoi chịu nhiệt ma sát khí động học và áp lực nén dọc trục khi phóng.
    • Công nghiệp năng lượng & hạt nhân: Ứng dụng trong thiết kế buồng đốt tuabin khí, vách phản ứng nhiệt hạch đòi hỏi vật liệu mặt trong chịu nhiệt $1200^\circ\text{C}$ và mặt ngoài tản nhiệt nhanh.
    • Quốc phòng & Hàng hải: Thiết kế thân tàu ngầm chịu áp lực thủy tĩnh sâu và vỏ xe bọc thép chống xuyên phá.

Limitations và Future Research

Mặc dù đạt được những đột phá xuất sắc, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Giả thiết điều kiện biên: Nghiệm giải tích mới thiết lập hoàn chỉnh cho điều kiện biên tựa đơn bốn cạnh ($S-S-S-S$), chưa áp dụng trực tiếp được dạng hiển cho biên ngàm cứng ($C-C$) hoặc biên tự do hỗn hợp.
  2. Hình học hoàn hảo ban đầu: Luận án tập trung vào vỏ nón hoàn hảo, chưa tính đến ảnh hưởng của độ bất hoàn hảo hình học ban đầu (initial geometric imperfections) dạng hàm ngẫu nhiên.
  3. Mô hình biến dạng trượt: Mới dừng ở lý thuyết bậc nhất (FSDT) đòi hỏi hệ số hiệu chỉnh cắt $\kappa = 5/6$, chưa áp dụng lý thuyết biến dạng trượt bậc cao (HSDT) phi hệ số điều chỉnh.
  4. Mô hình tiếp xúc gân: Sử dụng kỹ thuật san đều gân Lekhnitskii, chưa mô tả được tập trung ứng suất cục bộ tại chân gân khi bước gân quá lớn ($n_s < 6$).

Chương trình nghiên cứu phát triển tiếp theo bao gồm:

  • Mở rộng phân tích mất ổn định động lực (dynamic buckling) dưới tác dụng của tải xung kích nổ và sóng xung kích nhiệt độ cao.
  • Khảo sát ảnh hưởng của độ rỗng vi mô (porosity) và vật liệu composite gia cường hạt nano ống carbon (FG-CNTRC).
  • Nghiên cứu nứt gãy và lan truyền vết nứt mỏi trong vỏ nón FGM có gân gia cường.
  • Phát triển phần mềm tích hợp giao diện GUI tính toán ổn định tối ưu kết cấu vỏ phục vụ trực tiếp kỹ sư công trình.

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

Luận án của TS. Đỗ Quang Chấn tạo ra những tác động học thuật và công nghệ sâu rộng:

  • Tác động học thuậ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/Scopus đầu ngành uy tín cao trong lĩnh vực Cơ học kết cấu và Vật liệu composite (như Composite Structures, Thin-Walled Structures, Aerospace Science and Technology), đóng góp vào việc khẳng định vị thế của trường phái Cơ học Vật rắn Việt Nam trên trường quốc tế.
  • Chuyển dịch công nghệ công nghiệp: Đặt nền móng khoa học tính toán chính xác cho các viện nghiên cứu hàng không vũ trụ và đóng tàu quân sự trong nước tiếp cận thiết kế kết cấu thế hệ mới siêu nhẹ - siêu bền - chịu nhiệt siêu hạn.
  • Đóng góp giáo dục và đào tạo: Trở thành tài liệu tham khảo mẫu mực cho chương trình đào tạo sau đại học chuyên ngành Cơ học kỹ thuật, Cơ học vật rắn và Khoa học vật liệu tại các trường đại học hàng đầu.

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

Đối tượng Giá trị thụ hưởng cụ thể Lợi ích định lượng
Nghiên cứu sinh & Nhà khoa học Bộ công thức giải tích đóng tường minh cho bài toán ổn định và dao động vỏ nón FGM. Tiết kiệm 80% thời gian thiết lập mô hình tính toán sơ bộ so với phương pháp số thuần túy.
Kỹ sư R&D Hàng không - Quốc phòng Biểu đồ thiết kế tối ưu hóa số lượng gân ($n_s, n_r$), chiều cao gân ($h_s, h_r$) và tỷ phần thể tích gốm $k$. Giảm trọng lượng kết cấu vỏ từ 15–25% trong khi vẫn đảm bảo tải trọng phá hủy an toàn.
Kỹ sư Thiết kế Bồn bể & Lò phản ứng Giải pháp tính toán gradient nhiệt độ và độ chênh lệch nhiệt tới hạn $\Delta T_{cr}$. Nâng cao tuổi thọ làm việc của vách trao đổi nhiệt thêm 30–40% trong điều kiện sốc nhiệt.
Cơ quan Hoạch định Tiêu chuẩn Kỹ thuật Cơ sở khoa học thực chứng để xây dựng các bộ tiêu chuẩn kiểm định an toàn cho vỏ chịu áp lực vật liệu mới. Chuẩn hóa quy trình đánh giá rủi ro mất ổn định kết cấu mỏng thành mỏng tròn xoay.

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ì?

Đóng góp độc đáo nhất là việc mở rộng thành công Lý thuyết biến dạng trượt bậc nhất (FSDT) kết hợp quan điểm san đều gân Lekhnitskii có xét đến khoảng cách gân thay đổi $d_s(x) = (2\pi x \sin\beta)/n_s$ để thiết lập được hệ phương trình vi phân biến thiên và đưa ra nghiệm hiển xác định tải nén tới hạn $P_{cr}$ và nhiệt độ tới hạn $T_{cr}$ của vỏ nón cụt FGM.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây?

So với nghiên cứu của Sofiyev [146] (sử dụng phương pháp hàm ứng suất theo CST cho vỏ trơn) và Naj et al. [156] (sử dụng DQM số cho vỏ không gân), luận án đã giải quyết bài toán đặt theo 5 hàm chuyển vị theo FSDT cho vỏ có gân trực giao, đồng thời loại bỏ điểm kỳ dị toán học bằng kỹ thuật nhân tử tọa độ $x^i$, cho phép tính tích phân Galerkin hoàn toàn dưới dạng giải tích sơ cấp.

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

Phát hiện bất ngờ nhất là ảnh hưởng của chỉ số thể tích $k$: Khi $k$ tăng (tăng hàm lượng kim loại), tải tới hạn $P_{cr}$ giảm mạnh phi tuyến ở giai đoạn đầu ($k = 0 \to 2$), nhưng sau đó từ $k = 5 \to \infty$, tải tới hạn biến thiên rất chậm. Điều này chỉ ra rằng chỉ cần pha trộn 20–30% thể tích gốm đã có thể phục hồi tới 70% độ cứng mất ổn định của kết cấu thuần gốm.

4. Giao thức lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ?

Toàn bộ các ma trận độ cứng $A_{ij}, B_{ij}, D_{ij}$, các hệ thức gân Lekhnitskii, biểu thức giải tích tường minh của các số hạng vi phân $R_{ij}$ và hệ số ma trận $X_{ij}$ đều được trình bày chi tiết trong phần phụ lục của luận án, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập mã nguồn trên Maple hoặc Matlab.

5. Định hướng lộ trình nghiên cứu 10 năm tới?

Mở rộng từ bài toán vỏ nón cụt sang hệ kết cấu phức hợp vỏ nón - trụ - cầu FGM liền khối, tích hợp vật liệu nano graphen/carbon nanotube (FG-GPLRC), khảo sát hiện tượng mất ổn định đàn dẻo kết hợp từ biến trong môi trường bức xạ hạt nhân cực hạn.

Kết luận

Luận án tiến sĩ của Đỗ Quang Chấn là một công trình khoa học xuất sắc, mẫu mực và công phu, tạo lập bước tiến mang tính bước ngoặt trong cơ học kết cấu composite tiên tiến:

  1. Thiết lập thành công hệ phương trình giải tích FSDT hoàn chỉnh cho vỏ nón cụt FGM gia cường gân trực giao lệch tâm trên nền đàn hồi Pasternak.
  2. Xử lý triệt để bài toán toán học phức tạp của hệ phương trình vi phân có hệ số biến thiên bằng kỹ thuật nhân tử Galerkin kết hợp nghiệm giải tích đóng.
  3. Làm sáng tỏ toàn diện cơ chế mất ổn định tĩnh và ứng xử sau tới hạn dưới các trường tải cơ, tải nhiệt phân bố đều, tải nhiệt biến thiên tuyến tính và tải kết hợp cơ - nhiệt.
  4. Khám phá các quy luật động lực học phi tuyến và dao động tự do của panel nón cụt FGM tích hợp lớp áp điện điều khiển chủ động.
  5. Đưa ra các phát hiện khoa học định lượng có giá trị thực tiễn cao, cung cấp cơ sở dữ liệu tính toán tin cậy cho ngành công nghiệp hàng không vũ trụ, hạt nhân và công trình ngầm hiện đại.
  6. Mở ra ít nhất 3 hướng nghiên cứu mới: Cơ học vỏ nón FG-CNTRC nano, dao động phi tuyến ngẫu nhiên và tối ưu hóa đa mục tiêu kết cấu FGM thông minh.