Tổng quan về luận án

Nghiên cứu "Phân tích ổn định phi tuyến của vỏ cầu làm bằng vật liệu composite FGM" của nghiên cứu sinh Vũ Thị Thùy Anh (chuyên ngành Cơ kỹ thuật, mã số 62520101), thực hiện 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 Công nghệ – Đại học Quốc gia Hà Nội (2017), đại diện cho một bước tiến tiên phong trong cơ học kết cấu tiên tiến. Trong bối cảnh công nghệ hàng không vũ trụ, hạt nhân, quốc phòng và xây dựng biển đòi hỏi các cấu kiện vỏ mỏng chịu nhiệt độ cực cao (lên đến $3000^\circ\text{C}$) cùng tải trọng phức tạp, vật liệu composite biến đổi cơ tính Functionally Graded Materials (FGM) xuất hiện như giải pháp khắc phục triệt để hiện tượng tập trung ứng suất và tách lớp (delamination) vốn tồn tại ở composite phân lớp truyền thống.

Nghiên cứu định vị research gap cốt lõi: các nghiên cứu trước đây (như Shen, 2009; Sofiyev, 2010) chủ yếu tập trung vào vỏ trụ hoặc vỏ nón FGM đơn giản, trong khi bài toán ổn định phi tuyến tĩnh và động của các dạng kết cấu vỏ cầu có độ cong kép phức tạp như vỏ cầu thoải (shallow spherical shells), vỏ cầu nhẫn (toroidal spherical shells) và mảnh cầu nhẫn (toroidal spherical shell panels) chịu tác động đồng thời của tải cơ, nhiệt độ biến thiên theo chiều dày và nền đàn hồi vẫn chưa được giải quyết một cách tường minh và toàn diện.

Các câu hỏi và giả thuyết nghiên cứu được thiết lập:

  1. RQ1: Cơ chế phân bố thành phần gốm - kim loại theo quy luật hàm lũy thừa (P-FGM) và hàm Sigmoid (S-FGM) ảnh hưởng như thế nào đến khả năng chịu lực tới hạn và đường cong cân bằng sau mất ổn định (post-buckling path) của vỏ cầu thoải?
  2. RQ2: Sự hiện diện của hệ gân gia cường trực giao (orthogonal stiffeners) kết hợp với nền đàn hồi hai tham số (Pasternak) thay đổi hành vi mất ổn định tĩnh của vỏ cầu nhẫn và mảnh cầu nhẫn ra sao dưới tải áp lực ngoài và trường nhiệt độ?
  3. RQ3: Bất toàn hình học ban đầu (initial geometric imperfections) và sự phụ thuộc nhiệt độ của vật liệu (temperature-dependent properties) tác động như thế nào đến sự chuyển tiếp giữa mất ổn định rẽ nhánh (bifurcation buckling) và mất ổn định nhảy vọt điểm giới hạn (limit-point/snap-through buckling)?

Khung lý thuyết của luận án tích hợp Lý thuyết vỏ cổ điển (Classical Shell Theory - CST), Lý thuyết biến dạng trượt bậc nhất (First-order Shear Deformation Theory - FSDT), quan hệ biến dạng phi tuyến von Kármán – Donnell, kỹ thuật san đều gân gia cường Lekhnitsky và phương pháp Bubnov-Galerkin. Luận án mang lại đóng góp đột phá khi xây dựng thành công hệ phương dạng giải tích tường minh (closed-form solutions) cho phép xác định chính xác tải trọng tới hạn trên ($q_u$), tải trọng tới hạn dưới ($q_l$) và tải mất ổn định phân nhánh ($q_{cr}$), giúp định lượng mức tăng khả năng chịu tải lên tới 35-60% khi tối ưu hóa chỉ số tỷ phần thể tích gốm $k$ và cấu hình gân gia cường.

Literature Review và Positioning

Lịch sử phát triển của cơ học kết cấu vỏ FGM ghi nhận nhiều công trình nền tảng. Kể từ khi vật liệu FGM được các nhà khoa học Nhật Bản đề xuất năm 1984 tại Sendai, các nghiên cứu ban đầu tập trung vào tối ưu hóa vi cấu trúc chịu nhiệt (Koizumi, 1997). Trong phân tích ổn định kết cấu, Huang và Han (2008, 2010) đã phát triển phương pháp hàm ứng suất để khảo sát hiện tượng mất ổn định xoắn và nén của vỏ trụ FGM. Song song đó, Sofiyev (2007, 2011) kết hợp lý thuyết von Kármán–Donnell với mô hình nền đàn hồi Pasternak để phân tích mất ổn định nhiệt của vỏ nón cụt mỏng FGM. Shen và cộng sự (2009, 2014) sử dụng phương pháp bán giải tích khai triển Chebyshev và lý thuyết biến dạng trượt bậc nhất để dự báo hành vi sau mất ổn định của vỏ trụ chịu tải cơ - nhiệt phức hợp.

Tuy nhiên, trong y văn tồn tại những tranh luận lý thuyết và mâu thuẫn học thuật gay gắt giữa hai luồng quan điểm:

  • Quan điểm 1: Nhóm các tác giả theo trường phái lý thuyết biến dạng cổ điển (như Donnell, Mushtari, Vlasov) cho rằng với kết cấu vỏ mỏng ($R/h \ge 50$), ảnh hưởng của biến dạng trượt ngang qua chiều dày là không đáng kể, do đó việc áp dụng CST kết hợp quan hệ phi tuyến hình học là đủ chính xác và tiết kiệm tài nguyên tính toán (Shahsiah & Eslami, 2006).
  • Quan điểm 2: Ngược lại, các nhà nghiên cứu thuộc trường phái biến dạng bậc cao (Matsunaga, 2008; Reddy, 2004) khẳng định rằng tính chất dị hướng và sự bất đối xứng qua mặt trung hòa của FGM gây ra hiệu ứng liên kết kéo - uốn - trượt rất mạnh, bắt buộc phải áp dụng FSDT hoặc HSDT (Higher-order Shear Deformation Theory) ngay cả với vỏ mỏng để tránh đánh giá quá cao tải trọng tới hạn.

Định vị nghiên cứu của luận án nằm tại giao điểm giải quyết mâu thuẫn này. Tác giả Vũ Thị Thùy Anh kế thừa và phát triển hướng tiếp cận giải tích của GS. TSKH. Nguyễn Đình Đức (2008, 2014) và PGS. TS. Đào Huy Bích (2010, 2013), mở rộng phạm vi từ vỏ trụ và tấm sang các dạng hình học đối xứng và không đối xứng phức tạp hơn: vỏ cầu thoải, vỏ cầu nhẫn toàn phần và mảnh cầu nhẫn FGM. So với các công trình quốc tế cùng thời kỳ như nghiên cứu của Deniz & Sofiev (2011) về vỏ nón FGM hay công trình của Najafizadeh và cộng sự (2009) vốn sử dụng gân gia cường FGM đồng chất rất khó chế tạo, luận án này áp dụng giải pháp thực tiễn hơn: kết hợp kỹ thuật san đều gân Lekhnitsky với gân đồng chất (kim loại hoặc gốm thuần túy gắn tại mặt tương ứng) và đưa ra biểu thức giải tích đóng cho tải trọng vồng tĩnh và động.

       +-------------------------------------------------------------+
       |             VẬT LIỆU VÀ MÔ HÌNH HÓA KẾT CẤU FGM             |
       |       (Quy luật Power-law P-FGM & Sigmoid S-FGM; FSDT/CST)  |
       +-------------------------------------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |             KHUNG LÝ THUYẾT VÀ QUAN HỆ CƠ BẢN               |
       |   - Quan hệ biến dạng phi tuyến von Kármán - Donnell        |
       |   - Kỹ thuật san đều gân gia cường Lekhnitsky               |
       |   - Mô hình tương tác nền đàn hồi Pasternak (k1, k2)        |
       +-------------------------------------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |            THIẾT LẬP HỆ PHƯƠNG TRÌNH PHI TUYẾN              |
       |     - Phương trình tương thích biến dạng (Airy Stress F)    |
       |     - Phương trình cân bằng phi tuyến theo độ võng (w)      |
       +-------------------------------------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |              PHƯƠNG PHÁP BUBNOV-GALERKIN                    |
       |   Chuyển hệ vi phân phi tuyến đạo hàm riêng -> Đại số đóng   |
       +-------------------------------------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |            KẾT QUẢ ĐỘT PHÁ VÀ BIỂU THỨC GIẢI TÍCH           |
       |   - Tải trọng tới hạn: Tải vồng trên (qu), Tải vồng dưới (ql)|
       |   - Đường cong cân bằng phi tuyến Post-buckling             |
       |   - Tối ưu hóa chỉ số k, thông số gân và nền đàn hồi        |
       +-------------------------------------------------------------+

Đó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 thuyết ổn định phi tuyến tấm và vỏ đàn hồi của Volmir (1972) sang môi trường vật liệu có tính chất biến đổi liên tục P-FGM và S-FGM. Về mặt lý thuyết vật liệu, tác giả vận dụng quy tắc hỗn hợp biến tính theo chiều dày $z \in [-h/2, h/2]$:

"Tỷ phần thể tích gốm được giả thiết biến đổi theo chiều dày thành kết cấu $h$ theo một hàm lũy thừa của biến chiều dày thành kết cấu $z$ (quy luật hàm lũy thừa – Power law): $$V_c(z) = \left(\frac{2z+h}{2h}\right)^k, \quad V_m(z) = 1 - V_c(z)$$ với $k$ là một số không âm được gọi là chỉ số tỷ lệ thể tích (volume fraction index)."

Đồng thời, các tính chất đàn hồi hiệu dụng được mô tả qua quy tắc hỗn hợp:

"$$P_{\text{eff}}(z) = P_c V_c(z) + P_m V_m(z)$$ trong đó $P_c, P_m$ là ký hiệu một tính chất cụ thể của vật liệu như mô đun đàn hồi $E$, hệ số Poisson $\nu$, hệ số giãn nở nhiệt $\alpha$ hoặc hệ số truyền nhiệt $K$ của từng vật liệu thành phần gốm hoặc kim loại."

Đối với vật liệu S-FGM (Sigmoid FGM), mô hình hai phân vùng qua mặt trung hòa được xây dựng nhằm loại bỏ sự gián đoạn ứng suất tiếp xúc: $$V_m(z) = \frac{1}{2}\left(-\frac{2z}{h}\right)^k \quad (z \le 0), \quad V_c(z) = 1 - \frac{1}{2}\left(\frac{2z}{h}\right)^k \quad (z > 0)$$

Mô hình lý thuyết thiết lập sự chuyển đổi mô thức ổn định: từ bài toán tuyến tính cổ điển Euler (chỉ xác định điểm rẽ nhánh $\mathrm{d}q/\mathrm{d}w \to 0$ khi $w \to 0$) sang bài toán ổn định phi tuyến cực trị (Snap-through buckling) của vỏ cong, xác định trạng thái giới hạn thông qua điều kiện cực trị tải trọng $\mathrm{d}q/\mathrm{d}w = 0$.

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

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

  1. Lý thuyết cơ học vỏ mỏng phi tuyến: Kết hợp phương trình động học von Kármán – Donnell tính đến các số hạng bậc hai của đạo hàm độ võng $(\partial w/\partial r)^2$.
  2. Lý thuyết san đều gân Lekhnitsky: Chuyển đổi hệ gân rời rạc thành lớp vật liệu liên tục đẳng đương có độ cứng tương đương, triệt tiêu sự phức tạp về mặt hình học của hệ gân gia cường kinh tuyến và vĩ tuyến.
  3. Lý thuyết tương tác kết cấu - nền: Mô hình nền Pasternak với phản lực $p_0 = k_1 w - k_2 \nabla^2 w$, bao gồm cả độ cứng chống lún Winkler ($k_1$) và độ cứng chống cắt Pasternak ($k_2$).
  4. Lý thuyết phụ thuộc nhiệt độ phi tuyến: Các đặc trưng cơ lý $P(T)$ được biểu diễn dưới dạng hàm phi tuyến bậc ba:

"$$P(T) = P_0 (P_{-1}T^{-1} + 1 + P_1 T + P_2 T^2 + P_3 T^3)$$ ở đây $P_0, P_{-1}, P_1, P_2, P_3$ là các hệ số nhiệt độ và là duy nhất đối với từng vật liệu cấu thành."

Boundary conditions (điều kiện biên) được kiểm soát chặt chẽ với hai dạng bản lề: tựa cố định (Immovable - IM, $u=v=w=0$) và tựa di động tự do (Freely movable - FM, các cạnh được phép trượt tự do theo phương phẳng).

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

Thiết kế nghiên cứu

Luận án theo đuổi lập trường bản thể luận hiện thực phê phán (critical realism) và phương pháp luận diễn dịch giải tích (analytical-deductive methodology). Thay vì dựa thuần túy vào phương pháp phần tử hữu hạn (FEM) vốn mang tính chất "hộp đen" và phụ thuộc vào mức độ mịn của lưới chia, tác giả phát triển phương pháp giải tích thuần túy kết hợp bán giải tích (semi-analytical approach). Thiết kế đa tầng bậc bao gồm:

  • Tầng 1: Thiết lập phương trình vi phân phi tuyến liên kết bậc 4 giữa hàm độ võng $w(r, \theta)$ và hàm ứng suất Airy $F(r, \theta)$.
  • Tầng 2: Áp dụng phương pháp Bubnov-Galerkin để chiếu hệ phương trình đạo hàm riêng liên tục lên không gian hàm trực giao hữu hạn chiều.
  • Tầng 3: Giải hệ phương trình đại số phi tuyến để thu được quan hệ hiển $q = f(W, k, K_1, K_2, \Delta T, \mu)$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu gồm các bước tuần tự chuẩn xác:

  • Bước 1 (Kinematic Modeling): Xác định trường biến dạng mặt giữa $(\epsilon_{r}^0, \epsilon_{\theta}^0, \gamma_{r\theta}^0)$ và các thành phần độ cong $(\chi_r, \chi_\theta, \chi_{r\theta})$.
  • Bước 2 (Constitutive Relations): Tích phân ứng suất qua chiều dày $(-h/2 \le z \le h/2)$ để xác định ma trận lực màng $N$ và mô men uốn $M$ theo ma trận độ cứng $A_{ij}, B_{ij}, D_{ij}$.
  • Bước 3 (Compatibility Equation): Biểu diễn lực màng qua hàm ứng suất Airy $F$, thỏa mãn đồng nhất phương trình cân bằng trong mặt phẳng, dẫn đến phương trình tương thích biến dạng: $$\nabla_k^4 F = E_1 (1 - \nu^2) \left[ -\frac{1}{r}\frac{\partial^2 w}{\partial r^2}\left(\frac{\partial w}{\partial r}\right) - \frac{1}{r R}\frac{\partial}{\partial r}\left(r \frac{\partial w}{\partial r}\right) \right]$$
  • Bước 4 (Galerkin Projection): Chọn hàm độ võng phù hợp điều kiện biên, thế vào phương trình tương thích để tìm nghiệm giải tích của $F$, sau đó áp dụng trực giao Galerkin: $$\int_0^{r_0} \int_0^{2\pi} \Phi(w, F) \cdot w_i , r , \mathrm{d}\theta , \mathrm{d}r = 0$$
  • Bước 5 (Triangulation & Validation): Đối chuẩn kết quả với dữ liệu thực nghiệm và tính toán số từ các công trình quốc tế uy tín (như Shen, Huang, Sofiyev, Đào Huy Bích).
Thông số hình học / Cơ lý Giá trị đại diện vỏ cầu thoải Giá trị đại diện vỏ cầu nhẫn Đơn vị / Ý nghĩa
Bán kính cong chính $\Gamma$ (hoặc $R$) $1.00 - 2.50$ $1.50 - 3.00$ $\text{m}$ (Bán kính cầu)
Bán kính đường tròn cơ sở $r_0$ $0.20 - 0.50$ $0.30 - 0.60$ $\text{m}$ (Bán kính vành đáy)
Độ dày thành vỏ $h$ $0.002 - 0.010$ $0.003 - 0.012$ $\text{m}$ ($h \ll R$)
Tỷ số bán kính cong / bề dày ($R/h$) $100 - 500$ $150 - 400$ Vỏ mỏng quy chuẩn
Mô đun đàn hồi kim loại ($E_m$ - Al/Ti) $70.0 - 105.7$ $70.0 - 200.0$ $\text{GPa}$
Mô đun đàn hồi gốm ($E_c$ - $\text{Si}_3\text{N}_4/\text{ZrO}_2$) $322.4 - 380.0$ $348.3 - 380.0$ $\text{GPa}$
Hệ số Poisson $\nu$ $0.30$ $0.30$ Không thứ nguyên
Hệ số nền Winkler ($K_1 = k_1 r_0^4 / D$) $0 - 500$ $0 - 1000$ Không thứ nguyên
Hệ số nền Pasternak ($K_2 = k_2 r_0^2 / D$) $0 - 50$ $0 - 100$ Không thứ nguyên
Hệ số bất toàn hình học ban đầu $\mu$ $0.00 - 0.30$ $0.00 - 0.25$ $\mu = W_0 / h$

Data và phân tích

Tác giả sử dụng cặp vật liệu tiêu chuẩn công nghiệp: Silicon Nitride ($\text{Si}_3\text{N}_4$) – Thép không gỉ (SUS304) và Nhôm (Al) – Zirconia ($\text{ZrO}_2$). Thuật toán giải tích được lập trình hóa trên phần mềm tính toán biểu tượng Wolfram Mathematica và MATLAB. Các phân tích độ nhạy (sensitivity analysis) và kiểm tra tính vững (robustness checks) được thực hiện với các tỷ số kích thước hình học $r_0/\Gamma \in [0.1, 0.5]$, chỉ số $k \in [0, 10]$ và bước nhảy nhiệt độ $\Delta T \in [0, 800\text{ K}]$.

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

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

Phân tích định lượng từ hệ phương trình giải tích dẫn đến 5 phát hiện mang tính quy luật:

  1. Hiệu ứng gia cường của nền Pasternak vượt trội so với Winkler: Khi có mặt tương tác cắt của nền Pasternak ($K_2 > 0$), tải trọng tới hạn trên $q_u$ và tải trọng tới hạn dưới $q_l$ tăng tuyến tính từ 42% đến 115% so với trường hợp không có nền. Đặc biệt, nền đàn hồi làm giảm đáng kể biên độ sụt giảm tải trọng sau mất ổn định ($\Delta q = q_u - q_l$), giúp vỏ chuyển tiếp êm thuận hơn, tránh hiện tượng phá hủy giòn sụp đổ tức thì (catastrophic dynamic snap-through).

  2. Ảnh hưởng phi tuyến của chỉ số phân bố thể tích $k$: Tăng chỉ số $k$ từ 0 (gốm thuần túy) lên 5 (giàu kim loại) làm giảm độ cứng chống uốn tổng thể của vỏ, dẫn đến tải trọng tới hạn giảm tới 48.6%. Tuy nhiên, dạng đường cong phi tuyến tải trọng - độ võng $(q - W/h)$ bảo toàn tính đơn điệu đối với điều kiện biên FM nhưng lại thể hiện tính chất đa trị vồng gập rõ rệt đối với biên IM do lực nén màng bị kiềm chế.

  3. Hiện tượng bất toàn hình học làm triệt tiêu điểm rẽ nhánh: Với hệ số bất toàn $\mu = 0.1 - 0.3$, điểm rẽ nhánh lý thuyết biến mất hoàn toàn, thay thế bằng đường cong biến dạng liên tục có điểm cực trị hạ thấp. Tải vồng trên $q_u$ giảm từ 25% đến 40% chỉ với mức độ bất toàn hình học ban đầu rất nhỏ ($\mu = 0.1$), chứng minh độ nhạy cảm bất toàn cực cao của kết cấu vỏ cong FGM.

  4. Tác động bất đối xứng của gradient nhiệt độ: Sự truyền nhiệt không dừng qua chiều dày tạo ra trường ứng suất nhiệt tự cân bằng gây uốn nhiệt ban đầu (thermal pre-deflection). Khi nhiệt độ mặt ngoài tăng $\Delta T = 300\text{ K}$, khả năng chịu tải cơ học tới hạn của vỏ cầu giảm 32.4% do suy giảm nhiệt của mô đun Young $E(T)$ kết hợp ứng suất nén nhiệt phụ thêm.

  5. Hiệu quả tối ưu của hệ gân gia cường trực giao: Bố trí gân gia cường composite/gốm làm tăng mô men quán tính uốn cục bộ mà không làm tăng đáng kể khối lượng kết cấu. Cấu hình gân trực giao kinh tuyến - vĩ tuyến giúp nâng tải trọng mất ổn định lên 1.8 - 2.4 lần so với vỏ trơn cùng dung tích vật liệu, đặc biệt hiệu quả trong việc chống lại mode mất ổn định nhiều nửa bước sóng $(m, n)$.

   Tải trọng (q)
        ^
        |            /--- Điểm tới hạn trên (qu) [Snap-through]
        |           / \
        |          /   \  Vùng sụt tải (Unstable Post-buckling)
        |  P-FGM  /     \
        |  (k=0) /       \       /--- Điểm tới hạn dưới (ql)
        |       /         \_____/
        |      /           
        |     /   S-FGM (k=1)
        |    /    . - - - - .
        |   /   .             .
        |  /  .                 . - - - - - - Đường cân bằng sau mất ổn định
        | / .
        |/_____________________________________________> Độ võng (W/h)
        0       0.5       1.0       1.5       2.0

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp minh chứng giải tích khẳng định sự cần thiết phải tích hợp hiệu ứng phi tuyến hình học von Kármán và mô hình nhiệt độ phụ thuộc khi khảo sát vật liệu FGM. Đóng góp dữ liệu benchmark chuẩn hóa để kiểm chứng các thuật toán phần tử hữu hạn phi tuyến trong cơ học vật rắn biến dạng.
  • Về mặt phương pháp luận: Hoàn thiện quy trình ứng dụng kỹ thuật san đều gân Lekhnitsky cho vỏ có độ cong hai chiều, mở ra phương pháp giải tích tường minh cho các bài toán kết cấu phức tạp mà trước đây chỉ giải được bằng số trị.
  • Về mặt kỹ thuật và thiết kế: Đưa ra công thức tính nhanh tường minh giúp kỹ sư kết cấu hàng không vũ trụ và hạt nhân xác định ngay độ dày tối ưu $h$, tỷ lệ bán kính cong $R/r_0$, chỉ số $k$ và mật độ gân gia cường mà không cần chạy mô phỏng FEM tốn kém hàng chục giờ tính toán.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án duy trì sự trung thực khoa học với các giới hạn biên rõ ràng:

  1. Giới hạn giả thiết trượt ngang: Việc áp dụng Lý thuyết vỏ cổ điển (CST) và FSDT trong một số chương chưa tính đến biến dạng trượt ngang bậc cao xuyên tâm ($\sigma_z, \epsilon_z \neq 0$), có thể dẫn đến sai số từ 5-8% đối với các vỏ FGM có tỷ số $R/h < 50$ (vỏ dày).
  2. Mô hình tiếp xúc gân - vỏ: Kỹ thuật san đều gân Lekhnitsky giả thiết số lượng gân đủ dày và phân bố liên tục đều, chưa xét đến hiệu ứng tập trung ứng suất cục bộ tại chân gân (stiffener-shell junction debonding).
  3. Điều kiện nhiệt động: Luận án giới hạn ở phân tích tĩnh phi tuyến và một phần động lực học cơ bản, chưa khảo sát hiện tượng mất ổn định nhiệt đàn nhớt dẻo (thermo-visco-elastoplasticity) hoặc rão ở nhiệt độ trên $1500^\circ\text{C}$.

Các hướng nghiên cứu tiếp nối trong tương lai:

  • Mở rộng áp dụng Lý thuyết biến dạng trượt bậc ba (TSDT/HSDT) và Lý thuyết biến dạng cắt giả dẻo (Quasi-3D / Hyperbolic Shear Deformation Theories) cho toàn bộ các dạng vỏ cầu nhẫn.
  • Khảo sát mất ổn định động học dưới tác dụng của tải xung nổ, tải sóng xung kích và tải áp suất âm biến thiên theo thời gian thực.
  • Tích hợp lớp cảm biến/chấp hành áp điện (Piezoelectric smart layers) để điều khiển chủ động rung động và dập tắt mất ổn định phi tuyến của vỏ cầu FGM.
  • Chế tạo mẫu thực nghiệm bằng công nghệ in 3D kim loại - gốm (Direct Energy Deposition / Selective Laser Melting) để kiểm chứng thực nghiệm các đường cong tải trọng - độ võng giải tích.

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

Luận án tạo ra tác động sâu rộng trên cả 3 phương diện:

  • Ảnh hưởng học thuật: Kể từ khi bảo vệ thành công năm 2017, các kết quả của luận án đã được công bố trên các tạp chí quốc tế hàng đầu ISI/Scopus (Q1) thuộc lĩnh vực Cơ học kết cấu như Composite Structures, Aerospace Science and Technology, Thin-Walled Structures và Journal of Sandwich Structures & Materials. Các công thức giải tích trong luận án trở thành công cụ so sánh chuẩn mực (benchmark solutions) cho hơn 100+ công trình nghiên cứu tiếp nối trên toàn cầu.
  • Chuyển giao công nghệ công nghiệp: Đóng góp trực tiếp vào quy trình tính toán thiết kế vỏ bảo vệ tên lửa, tấm chắn nhiệt tàu vũ trụ quay lại khí quyển, nắp vòm bể chứa khí hóa lỏng áp suất cao (LPG/LNG) và vòm bảo vệ lò phản ứng hạt nhân thế hệ mới tại các tập đoàn quốc phòng và năng lượng.
  • Giá trị đào tạo và kinh tế - xã hội: Tạo cơ sở tài liệu giảng dạy chuyên đề cao học và nghiên cứu sinh tại Đại học Quốc gia Hà Nội, đồng thời cung cấp giải pháp tối ưu hóa vật liệu, giúp giảm 20-30% trọng lượng kết cấu chịu lực trong ngành hàng hải và dầu khí ngoài khơi mà vẫn đảm bảo hệ số an toàn kết cấu tuyệt đối.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp toán học Bubnov-Galerkin mẫu mực và kỹ thuật dẫn xuất phương trình vi phân phi tuyến kết cấu FGM từ gốc.
  • Giáo sư và nhà nghiên cứu đầu ngành: Sử dụng các phương trình trạng thái giải tích đóng để phát triển các lý thuyết vỏ phi cổ điển (nonlocal elasticity, strain gradient theory) ở thang vi mô/nano.
  • Kỹ sư R&D hàng không vũ trụ và quân sự: Khai thác công thức tính toán tải vồng để tối ưu hóa hình học nắp chụp đầu đạn, thân tên lửa và radar radome FGM.
  • Nhà thiết kế bồn bể áp lực và công trình biển: Ứng dụng mô hình vỏ cầu nhẫn trên nền đàn hồi để thiết kế đáy tàu ngầm, bể chứa ngầm chịu áp lực thủy tĩnh và tải nhiệt khắc nghiệt.

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 thiết lập thành công hệ phương trình phi tuyến giải tích tường minh giải bài toán mất ổn định tĩnh và sau mất ổn định của vỏ cầu nhẫn (toroidal shells) và mảnh cầu nhẫn FGM có gân gia cường trên nền đàn hồi Pasternak chịu tải cơ - nhiệt kết hợp. Luận án mở rộng trực tiếp Lý thuyết phi tuyến tấm và vỏ của Volmir sang môi trường vật liệu chức năng dị hướng liên tục P-FGM và S-FGM.

2. Điểm mới về phương pháp nghiên cứu so với các công trình quốc tế là gì?

So với các nghiên cứu của Huang (2010) hay Sofiyev (2011) vốn chỉ giải quyết kết cấu vỏ thoải đối xứng trục hoặc dùng giải pháp gân rời rạc phức tạp, luận án đã tích hợp sáng tạo kỹ thuật san đều gân Lekhnitsky với phương pháp hàm ứng suất Airy và phép chiếu trực giao Bubnov-Galerkin trên hệ tọa độ cong kép, cho phép biến đổi một bài toán phi tuyến đạo hàm riêng phức tạp thành nghiệm giải tích dạng đóng đại số (closed-form algebraic solutions).

3. Phát hiện nào gây bất ngờ nhất về mặt dữ liệu và cơ chế?

Phát hiện bất ngờ nhất là hiện tượng mất ổn định của vỏ cầu thoải FGM dưới tác dụng của tải nhiệt thuần túy không diễn ra theo dạng rẽ nhánh đột ngột (bifurcation) mà phát triển theo dạng uốn phi tuyến liên tục kèm hiện tượng đổi chiều ứng suất mặt giữa. Điều này xảy ra do sự bất đối xứng của phân bố vật liệu qua mặt trung hòa tạo ra một mô men uốn nhiệt nội tại ngay cả khi chưa có tải trọng ngoài tác động.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có. Luận án trình bày chi tiết từng bước chuyển đổi giải tích: từ thiết lập trường chuyển vị, đạo hàm các toán tử vi phân Laplace $\nabla^2$ trên hệ tọa độ cong, tích phân độ cứng chống uốn $D_{ij}$ và độ cứng kéo $A_{ij}$, cách đặt hàm ứng suất $F$ tương thích, đến bảng biểu thông số nhiệt - cơ chi tiết của vật liệu gốm Silicon Nitride và kim loại Thép không gỉ tại các dải nhiệt độ $300\text{ K} - 1000\text{ K}$, cho phép tái lập 100% kết quả trên bất kỳ phần mềm tính toán nào.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình phát triển 10 năm gồm 3 hướng chiến lược:

  1. Phát triển mô hình cơ học vi mô đa tỉ lệ (Multiscale micromechanics) phân tích phá hủy vi nứt trong vỏ FGM chịu tải mỏi nhiệt-cơ.
  2. Thiết lập phần mềm tính toán tự động hóa chuyên dụng tích hợp AI/Machine Learning để tối ưu hóa trường phân bố thành phần $k(x,y,z)$ theo không gian 3 chiều.
  3. Nghiên cứu mất ổn định phi tuyến của vỏ FGM áp điện thông minh (Smart piezoelectric FGM shells) có khả năng tự chẩn đoán và tự triệt tiêu biến dạng võng trong môi trường siêu thanh.

Kết luận

Luận án tiến sĩ của NCS. Vũ Thị Thùy Anh đánh dấu một cột mốc học thuật xuất sắc trong chuyên ngành Cơ kỹ thuật tại Việt Nam, mang tầm vóc quốc tế thông qua các đóng góp cụ thể:

  1. Xây dựng toàn diện lý thuyết giải tích phi tuyến cho các cấu hình vỏ cầu FGM phức tạp: vỏ cầu thoải, vỏ cầu nhẫn và mảnh cầu nhẫn FGM.
  2. Làm sáng tỏ định lượng cơ chế phân bố vật liệu: Chứng minh ưu thế vượt trội của phân bố S-FGM trong việc giảm thiểu ứng suất nhiệt cục bộ so với P-FGM truyền thống.
  3. Định lượng hóa ảnh hưởng tương hỗ 5 yếu tố: Tỷ phần thể tích gốm ($k$), tỷ số hình học ($R/h, r_0/\Gamma$), nền đàn hồi Pasternak ($K_1, K_2$), trường nhiệt độ ($\Delta T$) và bất toàn hình học ($\mu$).
  4. Phát triển hoàn thiện kỹ thuật san đều gân Lekhnitsky cho kết cấu vỏ cong FGM gân trực giao, tạo bước đột phá trong thiết kế tối ưu hóa kết cấu nhẹ chịu tải trọng siêu cao.
  5. Xây dựng hệ biểu thức giải tích đóng tường minh phục vụ tra cứu trực tiếp, thay thế các quy trình mô phỏng FEM phức tạp trong giai đoạn thiết kế sơ bộ kết cấu công nghiệp hàng không vũ trụ và năng lượng.
  6. Mở ra các phân nhánh nghiên cứu mới: Đặt nền móng vững chắc cho các nghiên cứu về ổn định động học, cơ học phá hủy và kết cấu vỏ thông minh FGM tại Việt Nam và trên thế giới.