Tổng quan về luận án

Nghiên cứu về ứng xử cơ học của các kết cấu tiên tiến hoạt động trong môi trường khắc nghiệt là một trong những trọng tâm cốt lõi của ngành Cơ học vật rắn đương đại. Kể từ khi được nhóm các nhà khoa học Nhật Bản đề xuất vào năm 1984, vật liệu có cơ tính biến đổi (Functionally Graded Materials - FGM) đã tạo ra một cuộc cách mạng trong thiết kế kết cấu chịu nhiệt độ cao nhờ loại bỏ hiện tượng tập trung ứng suất và hiện tượng tách lớp vốn tồn tại ở vật liệu composite phân lớp truyền thống. Trong số các dạng kết cấu không gian phức tạp, vỏ trống (toroidal shell segment) – một dạng vỏ tròn xoay hai độ cong sinh ra khi quay một cung tròn quanh trục nằm trong cùng mặt phẳng – đóng vai trò then chốt trong kỹ thuật hàng không vũ trụ (vỏ thân khí cụ bay, nắp mũi tên lửa), công nghiệp năng lượng (bình chịu áp, lò phản ứng hạt nhân) và kỹ thuật hàng hải (thân tàu ngầm, vỏ thiết bị lặn sâu).

Tuy nhiên, khoảng trống nghiên cứu (research gap) lớn tồn tại trong y văn quốc tế là phần lớn các công trình phân tích ổn định và động lực học của vỏ trống FGM trước đây chỉ dừng lại ở việc áp dụng Lý thuyết vỏ cổ điển (Classical Shell Theory - CST) hoặc Lý thuyết vỏ Donnell mỏng (Bích và cộng sự, 2013, 2016; Dũng và cộng sự, 2014, 2016; Thắng và Thời, 2016). Các lý thuyết này bỏ qua biến dạng trượt ngang ($\gamma_{xz} = \gamma_{yz} = 0$), dẫn đến sai số nghiêm trọng khi đánh giá các kết cấu vỏ trống có độ dày từ vừa đến lớn ($R/h \le 50$). Luận án tiến sĩ của Phạm Minh Vương (2023) với đề tài "Phân tích ổn định và động lực phi tuyến của vỏ trống có cơ tính biến thiên theo lý thuyết biến dạng trượt bậc cao", dưới sự hướng dẫn khoa học của GS. TSKH Nguyễn Đình Đức tại Trường Đại học Khoa học Tự nhiên - ĐHQGHN, đã giải quyết triệt để bài toán này bằng cách phát triển hệ phương trình giải tích tường minh hoàn chỉnh trên nền tảng Lý thuyết biến dạng trượt bậc ba (Third-order Shear Deformation Theory - TSDT) của Reddy và Liu (1987).

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

  • Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để thiết lập một khung giải tích chuẩn xác mô tả hành vi mất ổn định tĩnh phi tuyến (postbuckling) của vỏ trống FGM, vỏ sandwich FGM và vỏ trống có gân gia cường dưới tác dụng của các trường lực cơ - nhiệt phức hợp mà không cần hệ số hiệu chỉnh trượt?
  • Câu hỏi nghiên cứu 2 ($RQ_2$): Đặc trưng dao động phi tuyến và đáp ứng biên độ - tần số của vỏ trống FGM dày thay đổi ra sao khi chịu kích động áp lực ngoài biến thiên điều hòa theo thời gian?
  • Câu hỏi nghiên cứu 3 ($RQ_3$): Ngưỡng tải tới hạn động của vỏ trống FGM dưới tác dụng của tải nén dọc trục tăng tuyến tính theo thời gian được xác định chính xác theo tiêu chuẩn nào và chịu sự chi phối ra sao bởi tốc độ gia tải cùng các thông số hình học?
  • Giả thuyết nghiên cứu 1 ($H_1$): Việc áp dụng lý thuyết TSDT tích hợp phi tuyến hình học von Kármán sẽ khử bỏ sự đánh giá quá cao (overestimation) tải tới hạn của CST, mang lại lời giải hội tụ chính xác cho vỏ trống FGM dày.
  • Giả thuyết nghiên cứu 2 ($H_2$): Sự kết hợp giữa độ cong kinh tuyến (meridional curvature), nền đàn hồi Pasternak hai tham số và hệ thống gân gia cường trực giao sẽ tạo ra hiệu ứng tăng cứng phi tuyến vượt trội đối với cả vỏ trống lồi (convex) và vỏ trống lõm (concave).

Phạm vi nghiên cứu bao quát hệ vật liệu kim loại - gốm Si3N4/SUS304 với quy luật hàm lũy thừa (P-FGM), cấu trúc sandwich ba lớp và kết cấu gia cường gân trực giao, đặt trên mô hình nền đàn hồi Pasternak ($K_1, K_2$), thỏa mãn điều kiện biên tựa đơn bốn cạnh. Ý nghĩa khoa học và thực tiễn của công trình được định lượng thông qua việc công bố 10 công trình khoa học uy tín, trong đó có 6 bài báo trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/SCIE Q1 như Aerospace Science and Technology (IF = 5.5), Acta Mechanica, International Journal of Mechanics and Materials in Design (IF = 4.3), và Thin-Walled Structures (IF = 4.5).

Literature Review và Positioning

Lịch sử nghiên cứu cơ học vỏ trống bắt đầu từ các công trình kinh điển của Stein và McElman (1965, 1967) cùng Hutchinson (1967) tại Trung tâm Nghiên cứu NASA Langley và Đại học Harvard. Những tác giả này đã đặt nền móng cho giả thiết vỏ thoải (shallow shell assumption), coi vỏ trống là vỏ hai độ cong với bán kính cong kinh tuyến $a$ và bán kính xích đạo $R$. Tiếp sau đó, Stroud và Sykes (1969), Weingarten và cộng sự (1973), Oyesanya (2005) đã mở rộng nghiên cứu ổn định vỏ trống thuần nhất chịu nén và tải trọng ngẫu nhiên.

Sự xuất hiện của vật liệu composite FGM đã thúc đẩy làn sóng nghiên cứu mới từ đầu thập niên 2000. Shen (2004, 2012) tiên phong sử dụng phương pháp nhiễu để phân tích ổn định nhiệt và dao động biên độ lớn của vỏ trụ FGM. Praveen và Reddy (1998), Reddy (2000) thiết lập các mô hình phần tử hữu hạn dựa trên TSDT cho tấm FGM. Tại Việt Nam, nhóm nghiên cứu của Đào Văn Dũng, Nguyễn Đình Đức, Đào Huy Bích và Hoàng Văn Tùng đã có những đóng góp xuất sắc trong việc phân tích tấm và vỏ FGM. Tuy nhiên, khi chuyển sang kết cấu vỏ hai độ cong dạng vỏ trống, các nghiên cứu quốc tế và trong nước trước năm 2017 (như Bích và Ninh, 2016; Dũng và Vương, 2016) hầu như đều sử dụng CST.

Trong y văn tồn tại ba cuộc tranh luận học thuật (scholarly debates) lớn:

  1. CST đối kháng HSDT/TSDT: CST giả định đoạn pháp tuyến trước biến dạng vẫn thẳng và vuông góc với mặt giữa sau biến dạng ($\gamma_{xz} = \gamma_{yz} = 0$), làm triệt tiêu hiệu ứng trượt ngang. Đối với vỏ mỏng ($R/h \ge 200$), sai số giữa CST và TSDT là dưới 3%. Tuy nhiên, đối với vỏ dày trung bình ($R/h = 50$), CST đánh giá tải tới hạn cao hơn thực tế từ 12% đến 22%, tạo ra mối nguy hiểm tiềm tàng trong thiết kế chịu tải an toàn.
  2. Mô hình đồng hóa vi mô Voigt đối kháng Mori-Tanaka: Shen và Wang (2012) khi đánh giá dao động phi tuyến tấm FGM đã chỉ ra rằng mô hình Voigt và Mori-Tanaka cho kết quả tần số dao động và tải tới hạn chênh lệch nhau không đáng kể (dưới 2.5%), trong khi mô hình Voigt mang lại cấu trúc tích phân giải tích thuận lợi hơn rất nhiều để giải các hệ phương trình đạo hàm riêng phi tuyến.
  3. Dạng nghiệm Galerkin một số hạng đối kháng đa số hạng: Đối với bài toán xoắn và áp lực ngoài, trường biến dạng xuất hiện các mode sóng xiên và sóng kết hợp, đòi hỏi nghiệm hàm độ võng phải có dạng ba số hạng để nắm bắt chính xác hiện tượng mất ổn định rẽ nhánh.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Ali và Hasan (2019) (Proceedings of the Institution of Mechanical Engineers, Part C) khảo sát ổn định động của vỏ trống E-FGM dưới vận tốc nén không đổi bằng HSDT nhưng chỉ dừng lại ở phương pháp số và không có gân gia cường hay nền đàn hồi.
  • Nghiên cứu của Stein và McElman (1965) (AIAA Journal) phân tích ổn định vỏ trống thuần nhất bằng CST dưới tải cơ học thuần túy.

Luận án của Phạm Minh Vương đã định vị vượt trội khi xây dựng thành công lời giải giải tích giải trừ hoàn toàn hệ 5 phương trình đạo hàm riêng phi tuyến bậc cao về hệ 4 ẩn hàm, bao quát đầy đủ cả vỏ FGM đơn lớp, sandwich, có gân lệch tâm và tương tác nền Pasternak.

Đó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 biến dạng trượt bậc ba của Reddy và Liu (1987) cho bài toán vỏ trống hai độ cong FGM có xét đến phi tuyến hình học von Kármán. Trường chuyển vị được thiết lập chính xác tại mọi điểm $(x, y, z)$:

$$u(x, y, z) = u_0 - z \frac{\partial w_0}{\partial x} + z \left(1 - \frac{4z^2}{3h^2}\right) \phi_x$$

$$v(x, y, z) = v_0 \left(1 + \frac{z}{R}\right) - z \frac{\partial w_0}{\partial y} + z \left(1 - \frac{4z^2}{3h^2}\right) \phi_y$$

$$w(x, y, z) = w_0(x, y)$$

Mô hình này thỏa mãn tự nhiên điều kiện ứng suất tiếp bằng không tại hai bề mặt tự do $z = \pm h/2$, triệt tiêu hoàn toàn sự phụ thuộc vào hệ số hiệu chỉnh trượt $K_s$. Bằng cách tích hợp tenxơ nội lực bậc cao $P_{ij}$ và lực cắt bậc cao $R_i$, hệ thống phương trình cân bằng cơ học được chuyển đổi thành công sang hệ 4 phương trình đạo hàm riêng phi tuyến chứa hàm ứng suất Airy $F(x, y)$, độ võng $w(x, y)$ và hai góc xoay $\phi_x(x, y), \phi_y(x, y)$:

$$\Delta \Delta (\Phi_1 \phi_x + \Phi_2 \phi_y) + \mathcal{L}_1(w, F) = 0$$

$$\mathcal{L}_2(\phi_x, \phi_y, w, F) = q - K_1 w + K_2 \nabla^2 w$$

kết hợp với phương trình tương thích biến dạng:

$$\frac{1}{E_1} \Delta \Delta F = \left(\frac{\partial^2 w}{\partial x \partial y}\right)^2 - \frac{\partial^2 w}{\partial x^2}\frac{\partial^2 w}{\partial y^2} - \frac{1}{R}\frac{\partial^2 w}{\partial x^2} - \frac{1}{a}\frac{\partial^2 w}{\partial y^2}$$

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

Khung phân tích của luận án là sự hợp nhất của 5 lý thuyết và mô hình giải tích nền tảng:

  1. Lý thuyết biến dạng trượt bậc ba (Reddy & Liu, 1987): Mô hình hóa trạng thái ứng suất cắt phi tuyến qua chiều dày.
  2. Quan hệ biến dạng - chuyển vị phi tuyến von Kármán: Mô tả chính xác các biến dạng lớn và độ võng sau tới hạn.
  3. Mô hình nền đàn hồi hai tham số Pasternak (1954): Phản ánh tương tác cơ học đầy đủ thông qua độ cứng chống lún Winkler ($K_1$) và độ cứng lớp trượt cắt ($K_2$), mô tả qua lực phản chấn: $$q_e = K_1 w - K_2 \left(\frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2}\right)$$
  4. Kỹ thuật san đều gân Lekhnitskii (1968): Tính toán đóng góp độ cứng của hệ thống gân gia cường dọc (stringers) và gân vòng (rings) đặt lệch tâm vào các ma trận độ cứng tổng thể của kết cấu.
  5. Mô hình vật liệu P-FGM Voigt: Mô tả sự biến thiên cơ tính liên tục qua chiều dày $z \in [-h/2, h/2]$ với chỉ số thể tích $k$: $$E(z) = E_m + (E_c - E_m)\left(\frac{2z + h}{2h}\right)^k, \quad \rho(z) = \rho_m + (\rho_c - \rho_m)\left(\frac{2z + h}{2h}\right)^k$$

Điều kiện biên nghiên cứu được xác định chặt chẽ là liên kết tựa đơn (simply supported) bốn cạnh: $w = 0, M_x = 0, P_x = 0, \phi_y = 0$ tại các mép $x = 0, x = L$.

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 lập trường nhận thức luận Thực chứng (Positivism) kết hợp phương pháp suy diễn diễn dịch (Deductive-mechanistic paradigm) thuần lý thuyết giải tích và cơ học giải tích phi tuyến. Thiết kế nghiên cứu phân cấp đa tầng bao gồm:

  • Cấp độ vi mô/vật liệu: Mô hình hóa phân bố thể tích kim loại/gốm, tích phân xác định 13 hệ số cơ tính hiệu dụng $E_i$ ($i = 1 \dots 7$) và các tham số ứng suất nhiệt $\Phi_1, \Phi_2, \Phi_3, \Phi_4$.
  • Cấp độ trung mô/tiết diện: Tích phân tenxơ nội lực cơ bản ($N_{ij}, M_{ij}$) và nội lực bậc cao ($P_{ij}, Q_i, R_i$).
  • Cấp độ vĩ mô/kết cấu: Giải hệ phương trình vi phân đạo hàm riêng phi tuyến bằng phương pháp Galerkin và phương pháp cân bằng điều hòa.

Quy trình nghiên cứu rigorous

Quy trình giải tích được chuẩn hóa qua 4 bước:

  1. Thiết lập trường chuyển vị và hàm độ võng: Lựa chọn dạng nghiệm phù hợp cho từng bài toán. Với bài toán chịu áp lực ngoài và xoắn, hàm độ võng được khai triển 3 số hạng: $$w(x, y) = W_1 \sin(\alpha x)\sin(\beta y) + W_2 \sin^2(\alpha x) + W_3$$ Với bài toán nén dọc trục và bài toán động, sử dụng dạng nghiệm giải tích tường minh 1 số hạng: $w(x, y) = W \sin(m\pi x / L)\sin(n y / R)$.
  2. Xác định hàm ứng suất Airy $F(x, y)$ và góc xoay $\phi_x, \phi_y$: Giải chính xác phương trình tương thích biến dạng và hệ phương trình cân bằng bậc cao, dẫn xuất ra biểu thức giải tích của $F(x, y)$ chứa các số hạng phi tuyến bậc hai của độ võng $W^2$.
  3. Trực giao hóa Galerkin: Thực hiện phép tích phân biến phân trên toàn miền hình học của vỏ: $$\int_0^L \int_0^{2\pi R} \mathcal{R}(x, y) \sin\left(\frac{m\pi x}{L}\right)\sin\left(\frac{n y}{R}\right) dx dy = 0$$ dẫn đến hệ phương trình đại số phi tuyến bậc ba liên hệ giữa tải trọng và độ võng cực đại ($q - W_{max}/h$).
  4. Phân tích động lực học: Sử dụng phương pháp Runge-Kutta bậc 4 tích phân phương trình vi phân phi tuyến theo thời gian $t$. Áp dụng tiêu chuẩn Budiansky-Roth xác định điểm mất ổn định động tại điểm uốn (inflection point) của đường cong $W_{max}(t) - t$ nơi vận tốc biến dạng $d^2 W_{max} / dt^2$ đổi dấu đột ngột.

Data và phân tích

Thông số vật liệu chuẩn được sử dụng xuyên suốt là Gốm Silicon Nitride ($Si_3N_4$) và Kim loại Thép không gỉ (SUS304):

  • Gốm $Si_3N_4$: $E_c = 322.271 \text{ GPa}$, $\nu_c = 0.24$, $\alpha_c = 3.2 \times 10^{-6} \text{ K}^{-1}$, $\kappa_c = 2.09 \text{ W/mK}$, $\rho_c = 3000 \text{ kg/m}^3$.
  • Kim loại SUS304: $E_m = 207.788 \text{ GPa}$, $\nu_m = 0.318$, $\alpha_m = 12.46 \times 10^{-6} \text{ K}^{-1}$, $\kappa_m = 15.379 \text{ W/mK}$, $\rho_m = 8166 \text{ kg/m}^3$.
  • Hệ số Poisson trung bình $\nu = 0.3$.

Các chương trình tính toán số được lập trình trên nền tảng tính toán khoa học chuyên sâu, kiểm chứng tính hội tụ và độ tin cậy tuyệt đối thông qua việc so sánh đối chuẩn với các kết quả đã công bố của Shen (2004), Huang và Han (2010), Bích và Ninh (2016) trong các trường hợp suy biến về vỏ trụ ($a \to \infty$) và vỏ mỏng ($R/h \ge 200$).

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

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

  1. Hiệu ứng thành phần thể tích gốm ($k$): Tải tới hạn tĩnh và động tăng tỷ lệ thuận với hàm lượng gốm. Khi chỉ số thể tích giảm từ $k = 1$ xuống $k = 0.5$ (kết cấu giàu gốm hơn) tại tỷ số $R/h = 50$, tải tới hạn động nén dọc trục $P_{dcr}$ tăng 7.2% từ 1631.71 MPa lên 1748.65 MPa. Khi $k = 0.2$, tải tới hạn đạt tới 1879.21 MPa.
  2. Hiệu ứng độ dày hình học ($R/h$): Độ dày vỏ chi phối áp đảo khả năng kháng vồng. Với $k = 10$, khi vỏ dày lên tương ứng tỷ số $R/h$ giảm từ 200 xuống 100, tải tới hạn động tăng gần gấp đôi từ 472.96 MPa lên 943.44 MPa.
  3. Quy luật so sánh tải động - tải tĩnh ($P_{dcr} > P_{cr}$): Trong mọi trường hợp khảo sát, tải tới hạn động xác định theo tiêu chuẩn Budiansky-Roth luôn lớn hơn tải tới hạn tĩnh tương ứng. Ví dụ, tại $R/h = 50, k = 0.2$, tải động $P_{dcr} = 1879.21 \text{ MPa}$ so với tải tĩnh $P_{cr} = 1878.87 \text{ MPa}$; tại $R/h = 300, k = 10$, $P_{dcr} = 230.86 \text{ MPa}$ so với $P_{cr} = 230.16 \text{ MPa}$.
  4. Tính phi đối xứng nhạy cảm giữa vỏ lồi và vỏ lõm: Đối với vỏ trống lồi ($a > 0$), sự gia tăng bán kính kinh tuyến $a$ từ $5R$ đến $10R$ chỉ làm giảm nhẹ tải động từ 1633.78 MPa xuống 1631.69 MPa (chênh lệch dưới 0.13%). Ngược lại, đối với vỏ trống lõm ($a < 0$), kết cấu cực kỳ nhạy cảm với độ cong kinh tuyến: khi $a$ biến thiên từ $-10R$ sang $-20R$ (độ lõm giảm bớt), tải tới hạn động tăng vọt xấp xỉ 76% từ 487.23 MPa lên 855.43 MPa (tại $L = 100h$).
  5. Hiệu ứng gia cường của nền đàn hồi Pasternak và tốc độ đặt tải: Sự hiện diện của nền Pasternak làm tăng mạnh tải tới hạn. Tại $C = 100 \text{ MPa/s}$, từ trạng thái không nền $(K_1=0, K_2=0)$ lên nền Pasternak mức cao $(K_1 = 2 \times 10^8 \text{ N/m}^3, K_2 = 2 \times 10^5 \text{ N/m})$, tải tới hạn động tăng từ 1631.61 MPa lên 1836.27 MPa (tăng 12.54%). Tốc độ đặt tải $C$ tăng từ 10 MPa/s lên 300 MPa/s làm tải động tăng đơn điệu từ 1631.39 MPa lên 1632.36 MPa.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án hoàn thiện lý thuyết giải tích phi tuyến cho vỏ hai độ cong FGM dày, cung cấp hệ phương trình chuẩn mực có thể mở rộng trực tiếp cho các họ vật liệu dị hướng tiên tiến khác.
  • Về mặt phương pháp luận: Khẳng định tính ưu việt của phương pháp hàm ứng suất kết hợp Galerkin trong việc đưa bài toán TSDT 5 ẩn hàm phức tạp về dạng giải tích rút gọn có độ chính xác tương đương FEM nhưng tiết kiệm 95% chi phí tính toán.
  • Về mặt kỹ thuật thực tiễn: Cung cấp công thức hiển (explicit closed-form formulations) cho tải vồng và tần số dao động riêng, cho phép kỹ sư tính toán nhanh biên độ an toàn của cấu kiện vỏ trống trong tên lửa, tàu ngầm mà không cần chạy mô phỏng số cồng kềnh.
  • Về mặt tiêu chuẩn thiết kế: Khuyến cáo bắt buộc áp dụng lý thuyết TSDT cho các kết cấu vỏ FGM có $R/h < 100$ để tránh hiện tượng mất ổn định sớm do trượt cắ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 chính:

  1. Giả thiết vật liệu đàn hồi tuyến tính: Nghiên cứu chỉ xét biến dạng đàn hồi phi tuyến hình học, chưa tính đến biến dạng dẻo hoặc tính chất phi tuyến vật liệu của kim loại ở nhiệt độ siêu cao sát điểm nóng chảy.
  2. Điều kiện biên đơn giản hóa: Mới khảo sát điều kiện biên tựa đơn 4 cạnh, chưa mở rộng cho biên ngàm cứng (clamped), biên tự do hoặc liên kết đàn hồi bất đối xứng.
  3. Mô hình gân san đều: Phương pháp san đều Lekhnitskii giả định gân nhỏ và dày đặc, chưa mô phỏng được hiện tượng mất ổn định cục bộ của từng sườn gân đơn lẻ (local skin/stiffener buckling).
  4. Phân bố nhiệt độ đồng nhất: Giả thiết độ tăng nhiệt $\Delta T$ không đổi qua chiều dày, chưa giải phương trình truyền nhiệt Fourier phi tuyến 1D/2D qua bề dày thành vỏ.

Chương trình nghiên cứu 5 năm tiếp theo mở ra các hướng đột phá:

  • Mở rộng phân tích sang vật liệu nano composite gia cường ống nano carbon (FG-CNTRC) và tấm graphene (FG-GPLRC) hoặc cấu trúc xốp biến thiên (porous auxetic).
  • Nghiên cứu bài toán động lực học phá hủy vỏ trống có vết nứt dưới tác dụng của sóng xung kích nổ (blast loads) và trường điện từ coupled (electro-magneto-thermo-mechanics).
  • Áp dụng các phương pháp phần tử hữu hạn nâng cao đẳng hình học (Isogeometric Analysis - IGA) để mô hình hóa vỏ trống có lỗ khoét và liên kết biên thực tế.

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

  • Tác động học thuật: Hệ thống công thức giải tích của luận án đóng vai trò làm dữ liệu chuẩn (benchmark analytical data) phục vụ kiểm chứng các thuật toán số, phần mềm thương mại (Abaqus, Ansys) và các lý thuyết biến dạng trượt mới. Dự kiến các công bố từ luận án mang lại hàng trăm trích dẫn học thuật trên các tạp chí cơ học quốc tế uy tín.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp trong thiết kế khoang chịu áp của tàu lặn thám hiểm biển sâu, bình chứa nhiên liệu hydro lỏng trong hàng không vũ trụ và các vách ngăn bảo vệ lò phản ứng hạt nhân thế hệ mới, nơi vật liệu FGM Si3N4/SUS304 phát huy tối đa khả năng chịu nhiệt và áp suất nén.
  • Tác động xã hội và môi trường: Tối ưu hóa cấu trúc vỏ trống giúp giảm từ 15% đến 25% khối lượng vật liệu chế tạo khí cụ bay mà vẫn đảm bảo độ ổn định nhiệt - cơ, góp phần cắt giảm tiêu hao nhiên liệu và khí thải carbon trong công nghiệp hàng không.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận khung lý thuyết TSDT giải tích mẫu mực cho vỏ hai độ cong phức tạp, mở ra các hướng phát triển luận án tiến sĩ về vật liệu mới.
  • Kỹ sư R&D Hàng không vũ trụ và Quốc phòng: Sở hữu công cụ tính toán nhanh, chính xác tải vồng tĩnh, tải vồng động và vùng cộng hưởng phi tuyến của vỏ thiết bị bay.
  • Các tổ chức xây dựng tiêu chuẩn kỹ thuật: Bổ sung cơ sở khoa học để hiệu chỉnh các tiêu chuẩn an toàn kết cấu vỏ chịu áp lực ngoài và nhiệt độ cao.

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 tương thích và phương trình cân bằng giải tích phi tuyến hoàn chỉnh cho vỏ trống FGM hai độ cong dựa trên Lý thuyết biến dạng trượt bậc ba (TSDT) của Reddy, quy đổi bài toán 5 ẩn chuyển vị - góc xoay về hệ phương trình liên hệ giữa hàm ứng suất Airy $F(x, y)$, độ võng $w(x, y)$ và hai góc xoay $\phi_x, \phi_y$, từ đó tìm ra biểu thức giải tích tường minh cho tải vồng tĩnh và tần số dao động tự do.
  2. Tính đổi mới về phương pháp luận so với các nghiên cứu trước đây? So với Stein và McElman (1965) và Bích và Ninh (2016) vốn chỉ dùng lý thuyết vỏ cổ điển CST (bỏ qua trượt ngang), luận án giải quyết chính xác sự biến thiên phi tuyến của ứng suất trượt qua chiều dày mà không cần hệ số hiệu chỉnh trượt $K_s$. So với phương pháp số thuần túy (FEM) của Praveen và Reddy (1998), tiếp cận giải tích của luận án mang lại nghiệm giải tích hiển thị rõ bản chất vật lý và độ nhạy của từng tham số hình học.
  3. Phát hiện bất ngờ nhất từ dữ liệu số là gì? Sự nhạy cảm bất đối xứng cực lớn giữa vỏ trống lồi và vỏ trống lõm. Trong khi bán kính kinh tuyến $a$ hầu như không làm thay đổi tải tới hạn động của vỏ lồi ($< 0.13%$), thì đối với vỏ lõm, việc giảm độ cong từ $a = -10R$ sang $a = -20R$ làm tăng tải tới hạn động tới 76% (từ 487.23 MPa lên 855.43 MPa).
  4. Giao thức tái lập nghiên cứu (Replication protocol) có được cung cấp đầy đủ không? Hoàn toàn minh bạch và đầy đủ. Luận án cung cấp chi tiết toàn bộ các hệ số tích phân vật liệu hiệu dụng $E_1 \dots E_7$ (Phụ lục A1), các toán tử vi phân $A_1 \dots A_5$ (Phụ lục A2), biểu thức tường minh của hàm ứng suất và ma trận Galerkin từ Bảng A3 đến A13, cùng các thông số vật liệu $Si_3N_4$/SUS304 cụ thể, cho phép mọi nhà nghiên cứu tái lập chính xác 100% kết quả tính toán.
  5. Chương trình nghiên cứu 10 năm được vạch ra như thế nào? Định hướng mở rộng từ kết cấu FGM truyền thống sang cấu trúc siêu vật liệu (Auxetic metamaterials), vật liệu gia cường sợi graphene đa cấp độ, kết hợp các trường vật lý phức hợp (nhiệt - cơ - điện - từ), tính toán phá hủy phi tuyến và ứng dụng học máy (Machine Learning) để tối ưu hóa hình học vỏ trống đa mục tiêu.

Kết luận

Luận án tiến sĩ của Phạm Minh Vương đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính bước ngoặt:

  1. Xây dựng thành công khung lý thuyết giải tích giải quyết bài toán ổn định tĩnh, ổn định động và dao động phi tuyến cho kết cấu vỏ trống FGM dày dựa trên Lý thuyết biến dạng trượt bậc ba TSDT.
  2. Dẫn xuất trọn vẹn các biểu thức giải tích tường minh xác định tải tới hạn tĩnh, phương trình đường cong tải - độ võng sau tới hạn cho 4 bài toán ổn định tĩnh (chịu áp lực ngoài, tải xoắn, nén dọc trục kèm tải nhiệt, và vỏ có gân gia cường).
  3. Thiết lập biểu thức giải tích xác định tần số dao động tự do tuyến tính, quan hệ biên độ - tần số phi tuyến và phương trình vi phân dao động cưỡng bức dưới kích động điều hòa.
  4. Xác định chính xác tải tới hạn động của vỏ trống FGM chịu lực nén dọc trục tăng tuyến tính theo thời gian dựa trên tiêu chuẩn Budiansky-Roth, chứng minh quy luật $P_{dcr} > P_{cr}$.
  5. Định lượng chi tiết ảnh hưởng tương hỗ của chỉ số phân bố vật liệu $k$, tỷ số hình học $R/h, L/R, a/R$, hệ thống gân trực giao và nền đàn hồi Pasternak $(K_1, K_2)$ đến ứng xử cơ học của vỏ trống lồi và vỏ trống lõm.
  6. Mở ra 3 nhánh nghiên cứu mới trong cơ học kết cấu: phân tích động lực học phi tuyến vỏ hai độ cong vật liệu composite nano, tối ưu hóa kết cấu chịu nhiệt đa trường vật lý, và thiết lập chuẩn mực kiểm định cho các thuật toán phần tử hữu hạn phi tuyến tiên tiến trên phạm vi toàn cầu.