Tổng quan về luận án

Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 62520101) của nghiên cứu sinh Vũ Thị Thùy Anh, 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), mang tiê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ông trình đặt trọng tâm vào việc giải quyết bài toán ứng xử cơ học phi tuyến phức tạp của các kết cấu vỏ thoải dạng cầu, vỏ cầu nhẫn (spherical toroidal shells) và mảnh cầu nhẫn chế tạo từ vật liệu có cơ tính biến đổi liên tục (Functionally Graded Materials - FGM) dưới tác động đồng thời của trường lực cơ học, tải trọng nhiệt độ và môi trường tương tác đàn hồi.

Trong bối cảnh kỹ thuật hiện đại đòi hỏi các kết cấu hàng không vũ trụ, hạt nhân, quốc phòng và công nghiệp đóng tàu phải hoạt động ở môi trường siêu khắc nghiệt, vật liệu composite phân lớp truyền thống thường bộc lộ nhược điểm trí mạng: hiện tượng tập trung ứng suất nhiệt tại bề mặt liên kết giữa các lớp vật liệu dị bản dẫn tới sự bóc tách (delamination) và phá hủy kết cấu. Sự ra đời của vật liệu composite FGM từ năm 1984 đã khắc phục hạn chế này nhờ sự biến đổi trơn tru và liên tục của các tỷ phần thể tích thành phần từ bề mặt kim loại dẻo dai sang bề mặt gốm chịu nhiệt cao. Tuy nhiên, tính chất bất đối xứng về cơ tính, sự phi tuyến hình học theo quan hệ Von Kármán – Donnell và tính nhạy cảm với các khuyết tật ban đầu đã đặt ra khoảng trống nghiên cứu (research gap) lớn: phần lớn các nghiên cứu quốc tế trước đó (như Shen H.S., Sofiyev A.H., Shahsiah & Eslami) tập trung chủ yếu vào vỏ trụ tròn, vỏ nón hoặc tấm phẳng, trong khi các công trình giải tích tường minh cho vỏ cầu thoải, vỏ cầu nhẫn và mảnh cầu nhẫn FGM/S-FGM có gân gia cường trên nền đàn hồi còn rất thưa thớt và chưa đồng bộ.

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

  • RQ1: Quy luật biến thiên hàm mũ (P-FGM) và hàm Sigmoid (S-FGM) ảnh hưởng như thế nào đến ứng xử ổn định phi tuyến tĩnh và khả năng chịu tải sau mất ổn định (post-buckling) của vỏ cầu thoải?
  • RQ2: Tương tác cơ – nhiệt kết hợp (nhiệt độ tăng đều và truyền nhiệt dừng phi tuyến 1D Fourier qua chiều dày) thay đổi đặc tính mất ổn định dạng cực trị (snap-through) như thế nào trong điều kiện biên tựa tự do (FM) và tựa cố định (IM)?
  • RQ3: Kỹ thuật san đều gân gia cường Lekhnitsky kết hợp với mô hình nền đàn hồi hai tham số Winkler – Pasternak có thể tối ưu hóa độ cứng chống vồng của vỏ cầu nhẫn và mảnh cầu nhẫn FGM ra sao?
  • H1: Hiện tượng hóp (snap-through) xảy ra khi và chỉ khi biệt thức hình học – vật liệu thỏa mãn điều kiện $B^2 - AC > 0$, trong đó biên độ bước nhảy tải trọng $\Delta q = q_u - q_l$ tỷ lệ thuận với độ sâu hình học $r_0/R$.
  • H2: Mô hình vật liệu S-FGM cấu trúc sandwich phân bố đối xứng gốm – kim loại – gốm loại bỏ hoàn toàn hiện tượng uốn kết hợp do mất đối xứng trục trung hòa ($E_2 = 0$), giúp tăng độ cứng tĩnh lên hơn 18% so với P-FGM cùng thể tích gốm.
  • H3: Thành phần độ cứng trượt Pasternak ($K_2$) giữ vai trò chi phối áp đảo so với độ cứng nén Winkler ($K_1$) trong việc triệt tiêu độ võng sau mất ổn định và nâng cao tải tới hạn nhiệt.

Khung lý thuyết của luận án được tích hợp từ Lý thuyết vỏ mỏng 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ệ phi tuyến hình học Von Kármán – Donnell, lý thuyết môi trường đàn hồi Winkler – Pasternak và nguyên lý san đều gân Lekhnitsky. Đóng góp đột phá của công trình là thiết lập hệ phương trình vi phân phi tuyến liên hệ trực tiếp giữa hàm ứng suất $F$ và hàm độ võng $w$, giải quyết hoàn toàn bằng phương pháp giải tích Bubnov – Galerkin để đưa ra các công thức đóng (closed-form expressions) tường minh cho 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à độ biến thiên nhiệt độ tới hạn ($\Delta T_{cr}$). Phạm vi nghiên cứu bao quát từ tỷ số bán kính cong trên chiều dày $R/h \in [70, 90]$, tỷ số hình học $r_0/R \in [0.1, 0.94]$, với hệ tham số vật liệu điển hình gồm Hợp kim nhôm/Thép không gỉ ($E_m = 70\text{ GPa}, \alpha_m = 23 \times 10^{-6}\text{ }^\circ\text{C}^{-1}, K_m = 204\text{ W/mK}$) và Silicon Nitride/Zirconia ($E_c = 380\text{ GPa}, \alpha_c = 7.4 \times 10^{-6}\text{ }^\circ\text{C}^{-1}, K_c = 7.4\text{ W/mK}$).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước phản ánh sự phát triển đa nhánh trong cơ học kết cấu FGM:

  1. Dòng nghiên cứu vỏ trụ và vỏ nón FGM: Huang H. và Han Q. (2008–2010) đi tiên phong trong việc sử dụng phương pháp hàm ứng suất khảo sát mất ổn định xoắn và nén dọc trục của vỏ trụ FGM. Sofiyev A.H. cùng các cộng sự (2007–2015) triển khai các nghiên cứu giải tích sâu rộng về vỏ nón cụt FGM chịu tải nén động lực học biến thiên theo thời gian, áp lực thủy tĩnh và tải xung trên nền đàn hồi Pasternak bằng cách áp dụng phương pháp Galerkin. Shen H.S. (2009–2014) sử dụng phương pháp tiệm cận tham số bé và đa thức Chebyshev kết hợp lý thuyết FSDT và HSDT bậc cao để đánh giá đáp ứng sau mất ổn định nhiệt của vỏ trụ FGM có xét đến khuyết tật hình học.
  2. Dòng nghiên cứu vỏ thoải và tấm FGM: Shahsiah R. và Eslami M.R. (2003–2005) nghiên cứu ổn định nhiệt tuyến tính của vỏ cầu thoải FGM dựa trên hệ phương trình Donnell-Mushtari-Vlasov với ba điều kiện biên nhiệt. Ganapathi M. (2007) khảo sát mất ổn định nhiệt rẽ nhánh cho vỏ cầu FGM tựa bản lề trong giả định độ võng nhỏ. Tại Việt Nam, nhóm nghiên cứu của GS. Đào Huy Bích, GS. TSKH. Nguyễn Đình Đức, Đào Văn Dũng và Hoàng Văn Tùng (2008–2016) đã dẫn đầu hướng giải tích cho vỏ thoải hai độ cong, vỏ trụ sandwich S-FGM áp điện và vỏ nón cụt FGM có gân gia cường.

Tranh luận học thuật lớn tồn tại giữa hai quan điểm:

  • Trường phái 1 (Classical Shell Theory vs. Higher-order Shear Theory): Nhiều nhà nghiên cứu cho rằng CST bỏ qua biến dạng trượt ngang nên sẽ ước tính quá cao tải tới hạn của kết cấu dày ($R/h < 20$). Tuy nhiên, đối với vỏ mỏng và thoải ($R/h \ge 70$), CST cho kết quả giải tích tường minh cực kỳ chính xác với chi phí tính toán thấp.
  • Trường phái 2 (Hệ số Poisson phụ thuộc tọa độ vs. Hệ số Poisson hằng số): Một số công bố yêu cầu tích phân phức tạp hệ số Poisson $\nu(z)$ biến thiên theo chiều dày. Nhóm Nguyễn Đình Đức và Đào Huy Bích (2013) đã chứng minh bằng số trị rằng sự biến thiên của $\nu(z)$ chỉ tạo ra sai số không đáng kể ($< 1.2%$) so với khi xem $\nu = \text{const} = 0.3$, giúp đơn giản hóa hệ phương trình cơ bản mà vẫn bảo toàn độ tin cậy vật lý.

Vị trí của luận án được xác định rõ rệt: Luận án không dừng lại ở các dạng hình học cơ bản mà tiên phong tiến sâu vào các hình dạng đặc thù chưa từng được giải quyết toàn diện bằng phương pháp giải tích thuần túy, bao gồm vỏ cầu nhẫn FGM, mảnh cầu nhẫn FGM có gân gia cường trực giao và vỏ cầu S-FGM có xét sự phụ thuộc nhiệt độ của cơ tính vật liệu $P(T) = P_0(1 + P_{-1}T^{-1} + P_1T + P_2T^2 + P_3T^3)$.

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

  • So với công trình của Prakash T. et al. (2006) sử dụng FSDT kết hợp phần tử hữu hạn (FEM) số trị để phân tích vỏ cầu FGM ngàm cứng, luận án đạt được các biểu thức giải tích tường minh hiển thị rõ mối liên hệ giữa các tham số hình học và hệ số nền $K_1, K_2$, cho phép đánh giá tức thời trạng thái sau vồng mà không gặp lỗi hội tụ số.
  • So với nghiên cứu của Chao C.C. và Lin C.C. (1998) về nắp cầu cạn composite trực hướng, luận án mở rộng mô hình sang composite biến đổi hàm cơ tính hai chiều, tích hợp mô hình nền tương tác đàn hồi và trường nhiệt độ phi tuyến dừng.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra những bước tiến quan trọng đối với lý thuyết phi tuyến của cơ học vật rắn biến dạng:

  1. Mở rộng lý thuyết mất ổn định vỏ mỏng Donnell sang lớp vật liệu có cơ tính biến thiên bất đối xứng (P-FGM) và đối xứng (S-FGM). Khác với vật liệu đẳng hướng đồng nhất nơi ma trận độ cứng liên kết uốn – kéo nén $B_{ij} = 0$, ở vật liệu P-FGM, sự xuất hiện của hệ số $E_2 = \int_{-h/2}^{h/2} z E(z) dz \neq 0$ dẫn đến hiện tượng ghép đôi uốn – kéo (bending-stretching coupling).
  2. Xây dựng mô hình toán học giải tích tường minh xác định bản chất mất ổn định loại 2 (mất ổn định kiểu cực trị / snap-through) thay vì chỉ dừng lại ở bài toán mất ổn định loại 1 (mất ổn định rẽ nhánh Euler).
  3. Thiết lập hệ mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề 1: Tồn tại ngưỡng hình học tới hạn $(r_0/R){th}$ mà tại đó nếu tỷ số độ thoải $r_0/R > (r_0/R){th}$, hệ phương trình cân bằng phi tuyến luôn tồn tại hai nghiệm cực trị thực $W_1$ và $W_2$, kích hoạt bước nhảy độ võng không liên tục.
    • Mệnh đề 2: Sự hiện diện của mô hình cấu trúc Sigmoid-FGM đối xứng ba lớp ($c-m-c$ hoặc $m-c-m$) triệt tiêu triệt để thành phần ghép đôi uốn-kéo ($E_2 = 0$), đưa phương trình vi phân chuyển vị về dạng tương đương đẳng hướng nhưng với độ cứng uốn hiệu dụng $D_{S-FGM} = \frac{E_1 E_3}{E_1 (1-\nu^2)}$ vượt trội.
    • Mệnh đề 3: Trường nhiệt độ tăng đều tác động như một tải trọng nén trước trong mặt phẳng, làm dịch chuyển đường cong cân bằng tải – độ võng xuống dưới và thu hẹp khoảng cách giữa $q_u$ và $q_l$.

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

Khung phân tích được xây dựng trên sự dung hợp của 4 trụ cột lý thuyết:

  1. Quy luật phân bố pha vật liệu: Hàm lũy thừa $V_c(z) = (\frac{2z+h}{2h})^k$ đối với P-FGM và quy luật hàm mũ hai phía (Sigmoid law) đối với S-FGM: $$V_m(z) = \frac{1}{2}\left(\frac{2z+h}{h}\right)^k \quad (-h/2 \le z \le 0); \quad V_m(z) = 1 - \frac{1}{2}\left(\frac{h-2z}{h}\right)^k \quad (0 \le z \le h/2)$$
  2. Hình học vi phân vỏ thoải: Biến đổi tọa độ kinh tuyến – vĩ tuyến sang hệ tọa độ bán kính đường tròn vĩ tuyến $r = R \sin\varphi \approx R\varphi$ (với giả định $\cos\varphi \approx 1, Rd\varphi \approx dr$).
  3. Mô hình tương tác nền Pasternak: Tải trọng phản lực nền $q_e = k_1 w - k_2 \nabla^2 w$, phản ánh đồng thời độ lún đàn hồi cục bộ và lực cắt tương hỗ giữa các phần tử nền.
  4. Kỹ thuật san đều gân Lekhnitsky: Biến đổi các gân gia cường rời rạc thành lớp vật liệu tương đương có độ cứng uốn và độ cứng chống xoắn bổ sung được tính bình quân trên chu kỳ bước gân: $$D_1 = \frac{E_1 E_3 - E_2^2}{E_1(1-\nu^2)} + \frac{E_{stiff} I_{stiff}}{d_{stiff}}$$

Điều kiện biên được xác định rõ ràng tại chu tuyến ngoài $r = r_0$:

  • Biên ngàm tự do dịch chuyển kinh tuyến (Freely Movable - FM): $w = 0, \frac{\partial w}{\partial r} = 0, N_r(r_0) = 0$.
  • Biên ngàm cố định hoàn toàn (Immovable - IM): $w = 0, \frac{\partial w}{\partial r} = 0, u(r_0) = 0 \Rightarrow \int_0^{r_0} \frac{\partial u}{\partial r} dr = 0$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan bản thể luận thực chứng (Positivism) và phương pháp luận diễn dịch cơ học thuần lý (Deductive Mechanics). Thiết kế nghiên cứu là sự kết hợp giữa mô hình hóa giải tích (Analytical Modeling) và phương pháp bán giải tích (Semi-analytical Approach) nhằm tạo ra các thuật toán đánh giá trạng thái cân bằng phi tuyến với độ tin cậy toán học tuyệt đối.

Sơ đồ quy trình nghiên cứu đa tầng:

  1. Thiết lập phương trình động học phi tuyến Von Kármán – Donnell cho trường biến dạng mặt giữa $\varepsilon_r^0, \varepsilon_\theta^0, \gamma_{r\theta}^0$ và độ biến thiên độ cong $\chi_r, \chi_\theta, \chi_{r\theta}$.
  2. Thiết lập định luật Hooke tổng quát cho vật liệu FGM có kể đến biến dạng nhiệt $\varepsilon_{th} = \alpha(z, T)\Delta T$.
  3. Biểu diễn nội lực $N_r, N_\theta, N_{r\theta}$ và mô men uốn $M_r, M_\theta, M_{r\theta}$ qua các hàm chuyển vị và hàm ứng suất Airy $F(r, \theta)$.
  4. Tích phân dọc theo chiều dày $z \in [-h/2, h/2]$ để xác lập các độ cứng kéo $E_1$, độ cứng liên kết $E_2$, độ cứng uốn $E_3$ và các số hạng nhiệt $P_m, P_b$.
  5. Áp dụng phép chiếu trực giao Bubnov – Galerkin trên miền lấy tích phân hình học của vỏ để chuyển hệ phương trình đạo hàm riêng phi tuyến thành phương trình đại số phi tuyến liên hệ giữa tải trọng và độ võng không thứ nguyên $W/h$.
[Mô hình FGM/S-FGM] + [Cơ tính P(T)]
          │
          ▼
[Động học Von Kármán - Donnell] + [Nền Pasternak] + [Gân Lekhnitsky]
          │
          ▼
[Hệ phương trình vi phân phi tuyến đạo hàm riêng (w, F)]
          │
          ▼
[Phương pháp Giải tích Bubnov - Galerkin]
          │
          ▼
[Phương trình Đại số Phi tuyến Tường minh q = f(W, K1, K2, T, k, μ)]
          │
          ▼
[Xác định Tải tới hạn qu, ql và Đường cong Snap-Through]

Quy trình nghiên cứu rigorous

Để bảo đảm tính chặt chẽ trong việc giải phương trình vi phân phi tuyến bậc bốn đối xứng trục:

  1. Hàm độ võng $w(r)$ và hàm dạng khuyết tật ban đầu $w^(r)$ được lựa chọn thỏa mãn chính xác các điều kiện biên ngàm: $$w(r) = W \left(1 - \frac{r^2}{r_0^2}\right)^2, \quad w^(r) = \mu W \left(1 - \frac{r^2}{r_0^2}\right)^2$$ trong đó $\mu \in [0, 1]$ là hệ số khuyết tật hình học ban đầu không thứ nguyên.
  2. Hàm ứng suất Airy $F(r)$ được xác định chính xác bằng cách thay $w(r)$ vào phương trình tương thích biến dạng và tích phân liên tiếp: $$\nabla^2_s \left(\nabla^2_s F\right) = -\frac{E_1}{R}\nabla^2_s w - \frac{E_1}{2} L(w, w + 2w^*)$$ Các hằng số tích phân được triệt tiêu dựa trên điều kiện biên hữu hạn về ứng suất tại đỉnh vỏ ($r = 0$) và điều kiện cân bằng biên tại chu tuyến ($r = r_0$).
  3. Kiểm tra tính hợp lệ (Validation & Triangulation): Kết quả giải tích của luận án được so sánh đối chiếu đa kênh:
    • Khi cho hệ số $k \to 0$ hoặc $k \to \infty$, kết quả suy biến chính xác về vỏ cầu đồng nhất đẳng hướng kim loại thuần túy hoặc gốm thuần túy.
    • Khi cho $K_1 = K_2 = 0$ và $\mu = 0$, biểu thức liên hệ tải trọng – độ võng trùng khớp hoàn toàn với công thức của Đào Huy Bích và Hoàng Văn Tùng (2010) [16].
    • Kết quả phân tích ổn định tấm tròn suy biến khớp với nghiệm của Ma L.S. & Wang T.J. (2003) [62].

Data và phân tích

Dữ liệu tính toán được thiết lập từ các thông số nhiệt động học thực nghiệm chuẩn của vật liệu FGM Si3N4/SUS304:

  • Thép không gỉ (SUS304): $E_m = 207.788\text{ GPa}$ tại $300\text{ K}$, hệ số dãn nở $\alpha_m = 15.321 \times 10^{-6}\text{ K}^{-1}$, hệ số dẫn nhiệt $K_m = 12.04\text{ W/mK}$.
  • Silicon Nitride ($Si_3N_4$): $E_c = 322.271\text{ GPa}$ tại $300\text{ K}$, hệ số dãn nở $\alpha_c = 5.872 \times 10^{-6}\text{ K}^{-1}$, hệ số dẫn nhiệt $K_c = 9.19\text{ W/mK}$.
  • Đối với trường hợp phân tích tham số chuẩn: $E_m = 70\text{ GPa}, E_c = 380\text{ GPa}, \nu = 0.3$.

Phương trình liên hệ phi tuyến giữa áp lực phân bố đều $q$ và độ võng cực đại tại tâm $W$ cho vỏ cầu thoải P-FGM trên nền đàn hồi có dạng giải tích tổng quát: $$q = \left[\frac{64 D}{R_0^4 h^4} + \frac{3 E_1}{7 R^2 h^2} + \frac{16 D K_1}{21 R_0^4 h^4} + \frac{40 D K_2}{7 R_0^2 R^2 h^4}\right] W - \left[\frac{1385 E_1}{693 R_0^2 R h^3}\right] W^2 + \left[\frac{848 E_1}{429 R_0^4 h^4}\right] W^3$$ với $R_0 = r_0/h, \bar{R} = R/h$.

Điều kiện cực trị xác định điểm tải vồng trên ($q_u$) và điểm tải vồng dưới ($q_l$) được thiết lập từ đạo hàm bậc nhất: $$\frac{dq}{dW} = A - 2BW + 3CW^2 = 0 \implies W_{1,2} = \frac{B \mp \sqrt{B^2 - 3AC}}{3C}$$ Biệt thức $B^2 - 3AC$ quyết định hình thái ứng xử cơ học của vỏ.

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

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

  1. Quy luật chi phối của tỷ số độ sâu hình học $r_0/R$ lên hiện tượng hóp (Snap-through): Khi tỷ số $r_0/R$ tăng từ $0.3$ lên $0.94$ (vỏ cầu càng sâu), tải trọng tới hạn trên $q_u$ tăng mạnh nhưng đồng thời biên độ bước nhảy sụt giảm tải trọng $\Delta q = q_u - q_l$ tăng theo hàm phi tuyến bậc ba. Điều này chỉ ra rằng các kết cấu vỏ cầu sâu có khả năng chịu lực nén ban đầu cao hơn nhưng khi mất ổn định sẽ xảy ra hiện tượng sụp đổ động (dynamic snap-through collapse) cực kỳ mãnh liệt và nguy hiểm.
  2. Hiệu năng vượt trội của cấu trúc S-FGM so với P-FGM: Trong cùng điều kiện biên và tỷ phần thể tích gốm, vỏ cầu S-FGM cấu trúc sandwich phân lớp chức năng đối xứng nâng cao tải trọng tới hạn tĩnh lên từ $15%$ đến $23.5%$ so với vỏ P-FGM tiêu chuẩn, đồng thời làm giảm độ võng dư sau mất ổn định tới $30%$.
  3. Tác động ổn định hóa vượt trội của tham số nền cắt Pasternak ($K_2$): Phân tích số trị tại bảng tham số $r_0/R = 0.94$ cho thấy khi tăng hệ số Pasternak $K_2$ từ $0$ lên $20$, tải trọng tới hạn $q_u$ tăng gấp hơn $1.8$ lần, trong khi việc tăng hệ số Winkler $K_1$ từ $0$ lên $100$ chỉ làm tăng tải tới hạn khoảng $12%$. Nền cắt Pasternak có tác dụng kiềm chế mạnh mẽ sự phát triển của các sóng uốn cục bộ.
  4. Ảnh hưởng tiêu cực của trường nhiệt độ và độ nhạy khuyết tật: Khi độ chênh lệch nhiệt độ $\Delta T$ tăng từ $0\text{ K}$ lên $300\text{ K}$, khả năng chịu tải cơ học tới hạn $q_u$ của vỏ cầu FGM giảm tuyến tính từ $28%$ đến $42%$. Độ không hoàn hảo hình học ban đầu $\mu = 0.1$ làm suy giảm tải trọng tới hạn tới $18.4%$ và chuyển dịch điểm mất ổn định cực trị sang trạng thái biến dạng đơn điệu bất lợi.
  5. Hiệu ứng gia cường tối ưu của hệ gân trực giao: Việc bố trí gân gia cường đồng nhất (gân gốm tại mặt giàu gốm, gân kim loại tại mặt giàu kim loại) theo phương pháp Lekhnitsky giúp tăng tải nén tới hạn của vỏ cầu nhẫn FGM lên từ $2.1$ đến $3.4$ lần mà chỉ làm gia tăng khối lượng kết cấu dưới $15%$.

Implications đa chiều

  • Về mặt lý thuyết: Luận án cung cấp hệ thống công thức hiển (explicit analytical formulas) hoàn chỉnh cho cơ học vỏ FGM, đóng góp cơ sở dữ liệu benchmark chuẩn xác để kiểm chứng các thuật toán số (FEM, XFEM, Isogeometric Analysis - IGA).
  • Về mặt phương pháp luận: Khẳng định tính khả thi và độ chính xác tuyệt đối của phương pháp Bubnov – Galerkin trong việc giải quyết các bài toán biên phi tuyến hình học bậc cao của vật liệu có cơ tính biến thiên liên tục.
  • Về mặt thực tiễn kỹ thuật: Đưa ra các biểu đồ thiết kế kỹ thuật (Design Nomograms) cho phép các kỹ sư hàng không và đóng tàu lựa chọn chỉ số thể tích tối ưu $k \in [0.5, 2.0]$, tỷ số độ dày $R/h$ và mật độ gân gia cường để triệt tiêu hoàn toàn hiện tượng hóp nguy hiểm.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp luận cứ khoa học để xây dựng các tiêu chuẩn tính toán thiết kế vỏ chịu áp lực cao bằng vật liệu composite tiên tiến trong các lò phản ứng nhiệt hạch thế hệ mới và tàu ngầm thám hiểm biển sâu.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính biên điều kiện:

  1. Luận án chủ yếu sử dụng Lý thuyết vỏ mỏng cổ điển (CST), do đó các kết quả tính toán có độ chính xác cao nhất đối với các kết cấu có $R/h \ge 50$. Đối với các vỏ rất dày ($R/h < 20$), hiệu ứng biến dạng trượt ngang cần được mô tả bằng lý thuyết bậc cao (HSDT) để tránh đánh giá quá cao tải tới hạn.
  2. Phương trình truyền nhiệt dừng một chiều Fourier qua chiều dày mới chỉ giải quyết được nghiệm giải tích dạng chuỗi tường minh cho trường hợp chỉ số tỷ lệ thể tích $k = 1$. Với các giá trị $k$ bất kỳ, bài toán truyền nhiệt vẫn phải dựa vào tích phân số hoặc xấp xỉ từng đoạn.
  3. Luận án tập trung vào miền biến dạng đàn hồi tuyến tính của vật liệu FGM, chưa xét đến hiện tượng dẻo hóa của pha kim loại ở nhiệt độ siêu cao ($T > 1000\text{ K}$) hoặc sự xuất hiện của các vết nứt tế vi phân tán trong pha gốm.

Chương trình nghiên cứu tiếp nối trong tương lai (Future Research Agenda):

  • Hướng 1: Mở rộng phân tích mất ổn định động lực học phi tuyến (Nonlinear Dynamic Buckling) của vỏ cầu FGM dưới tác động của tải sóng nổ xung kích (blast loads) và tải nén thủy tĩnh biến thiên tuần hoàn theo thời gian bằng phương pháp Runge – Kutta bậc 4.
  • Hướng 2: Tích hợp mô hình vật liệu FGM áp điện (Piezoelectric FGM) và vật liệu FGM có lỗ xốp (Porous FGM) để thiết kế các kết cấu vỏ cầu thông minh có khả năng tự cảm biến và điều khiển chủ động dao động/mất ổn định.
  • Hướng 3: Nghiên cứu bài toán tối ưu hóa đa mục tiêu hình học và cấu trúc vi mô của hệ gân gia cường dựa trên các thuật toán tiến hóa (Genetic Algorithm, Particle Swarm Optimization).
  • Hướng 4: Thực hiện các kiểm chứng thực nghiệm (Experimental Validation) trên các mẫu vỏ cầu FGM chế tạo bằng công nghệ in 3D laser kim loại/gốm trực tiếp (Direct Energy Deposition).

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

  • Tác động học thuật: Luận án đã công bố nhiều công trình trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Composite Structures, Thin-Walled Structures, Journal of Sandwich Structures & Materials), tạo lập nguồn trích dẫn học thuật quan trọng cho cộng đồng nghiên cứu cơ học tính toán vật liệu composite tiên tiến toàn cầu.
  • Chuyển đổi công nghiệp: Kết quả nghiên cứu được ứng dụng trực tiếp trong việc thiết kế vỏ nắp chịu áp lực của các thiết bị lặn biển sâu, bồn chứa khí hóa lỏng (LNG), thân máy bay siêu thanh và nắp bảo vệ đầu đạn tên lửa hành trình – nơi mà sự kết hợp giữa trọng lượng siêu nhẹ và khả năng kháng nhiệt là yếu tố sống còn.
  • Lợi ích kinh tế – xã hội: Việc tối ưu hóa chỉ số thể tích $k$ và bố trí gân gia cường hợp lý giúp tiết kiệm từ $20%$ đến $35%$ chi phí nguyên liệu gốm kỹ thuật đắt tiền mà vẫn duy trì hệ số an toàn kết cấu vượt trội, kéo dài tuổi thọ công trình và giảm thiểu nguy cơ tai nạn do sụp đổ kết cấu vỏ mỏng.

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

  • 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 chuyển vị – hàm ứng suất và quy trình giải tích trực giao Bubnov – Galerkin chuẩn mực cho kết cấu vỏ mỏng dị bản.
  • Các nhà khoa học cơ học kết cấu: Sử dụng các nghiệm đóng giải tích làm nghiệm chuẩn (benchmark analytical solutions) để phát triển và hiệu chỉnh các phần mềm phần tử hữu hạn thương mại như ABAQUS, ANSYS, COMSOL Multiphysics.
  • Kỹ sư R&D trong công nghiệp Hàng không & Quốc phòng: Nắm bắt quy luật tương tác cơ – nhiệt – nền đàn hồi để định cỡ kích thước hình học ($R, h, r_0$) và lựa chọn vật liệu FGM tối ưu trong các thiết kế chịu tải cực trị.
  • Cơ quan quản lý và Viện thiết kế công trình ngầm/biển: Có thêm cơ sở khoa học định lượng để xây dựng quy chuẩn an toàn cho kết cấu vỏ vòm composite trong môi trường xâm thực hóa – 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ì và đã mở rộng lý thuyết nào?
Trả lời: Đó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 liên hệ giữa tải trọng cơ – nhiệt – thông số nền đàn hồi và độ võng cho các hình dạng vỏ cầu đặc thù (vỏ cầu nhẫn và mảnh cầu nhẫn FGM). Luận án đã mở rộng Lý thuyết màng và uốn vỏ mỏng Donnell – Mushtari – Vlasov kết hợp với nguyên lý san đều gân Lekhnitsky cho môi trường vật liệu có tính chất biến đổi liên tục bất đối xứng, giải quyết tường minh sự mất đối xứng trục trung hòa ($E_2 \neq 0$) gây ra hiệu ứng uốn – kéo kết hợp.

2. Đổi mới phương pháp luận của luận án so với các công trình quốc tế cùng thời kỳ?
Trả lời: So với cách tiếp cận thuần túy số trị của Prakash et al. (2006) hay phương pháp tiệm cận nhiễu loạn của Shen H.S. (2009), luận án sử dụng phương pháp giải tích Bubnov – Galerkin dựa trên hàm dạng độ võng bậc bốn hai lớp thỏa mãn chính xác các điều kiện biên động học và tĩnh học ngàm cố định (IM) và ngàm trượt (FM). Cách tiếp cận này đưa ra được nghiệm hiển trực tiếp của tải trọng tới hạn $q_u, q_l$ mà không phụ thuộc vào bước lặp hội tụ phi tuyến Newton-Raphson.

3. Phát hiện số trị gây bất ngờ nhất có bằng chứng thực nghiệm/tính toán là gì?
Trả lời: Đó là tính chất hai mặt của độ sâu vỏ cầu $r_0/R$: khi tăng tỷ số $r_0/R$ từ $0.3$ lên $0.94$, khả năng chịu tải ban đầu $q_u$ tăng mạnh, nhưng bước nhảy sụp đổ $\Delta q = q_u - q_l$ tăng đột biến khiến kết cấu rơi vào trạng thái mất ổn định động sụp hóp (dynamic snap-through) cực kỳ nguy hiểm. Một phát hiện khác là tác động chế ngự mất ổn định của nền cắt Pasternak $K_2$ mạnh gấp nhiều lần so với nền đàn hồi Winkler $K_1$.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Trả lời: Toàn bộ quy trình giải tích được minh bạch hóa $100%$ trong văn bản: từ hệ phương trình vi phân tương thích (1.12), định luật Hooke (1.13), các tích phân độ cứng chiều dày $E_1, E_2, E_3$ (1.18), hàm độ võng giả định $w(r)$ (2.11), nghiệm giải tích tường minh của hàm Airy $F(r)$ (2.13), cho đến phương trình tải trọng – độ võng cuối cùng (2.16) và các phụ lục chi tiết. Bất kỳ nhà nghiên cứu nào cũng có thể lập trình lại trên phần mềm Maple, Mathematica hoặc MATLAB để tái tạo chính xác các đường cong cân bằng phi tuyến.

5. Chương trình nghiên cứu 10 năm được vạch ra từ luận án gồm những gì?
Trả lời: Lộ trình phát triển 10 năm bao gồm: (i) Chuyển đổi từ phân tích tĩnh sang phân tích động lực học phi tuyến và dao động hỗn loạn (Chaos) của vỏ FGM; (ii) Tích hợp trường đa vật lý: Nhiệt – Điện – Từ – Cơ (Electro-Magneto-Thermo-Mechanical); (iii) Áp dụng lý thuyết biến dạng trượt bậc cao HSDT và lý thuyết mỏi phi cổ điển (Nonlocal Strain Gradient Theory) cho các kết cấu vỏ micro/nano-FGM; (iv) Nghiên cứu thực nghiệm phá hủy trên vật liệu FGM chế tạo bằng công nghệ in 3D kim loại tiên tiến.

Kết luận

  1. Xây dựng thành công hệ phương trình cơ bản và thiết lập các công thức giải tích đóng tường minh cho bài toán ổn định phi tuyến tĩnh của vỏ cầu thoải, vỏ cầu nhẫn và mảnh cầu nhẫn làm bằng vật liệu composite P-FGM và S-FGM trên nền đàn hồi Winkler – Pasternak.
  2. Làm sáng tỏ cơ chế mất ổn định kiểu cực trị (snap-through) dưới tác động đồng thời của tải trọng cơ học và trường nhiệt độ phi tuyến (nhiệt độ tăng đều và truyền nhiệt dừng 1D Fourier), xác định chính xác các biểu thức tải tới hạn trên $q_u$, tải tới hạn dưới $q_l$ và bước nhảy sụt tải $\Delta q$.
  3. Khẳng định ưu thế vượt bậc của cấu trúc vật liệu chức năng phân bố Sigmoid (S-FGM) ba lớp trong việc loại trừ ứng suất tiếp xúc liên lớp và triệt tiêu thành phần biến dạng uốn – kéo ghép đôi ($E_2 = 0$), giúp nâng cao độ ổn định cơ nhiệt của kết cấu.
  4. Chứng minh vai trò gia cường chịu lực đặc biệt hiệu quả của hệ gân trực giao theo lý thuyết san đều Lekhnitsky và khả năng kiềm chế độ võng hữu hiệu của thành phần nền đàn hồi Pasternak ($K_2$).
  5. Mở ra 3 hướng nghiên cứu mới mang tính liên ngành: Ổn định động lực học phi tuyến vỏ FGM chịu tải sóng xung kích; Cơ học kết cấu thông minh Smart FGM tích hợp lớp áp điện; Tối ưu hóa vi cấu trúc FGM bằng trí tuệ nhân tạo và vật liệu in 3D.
  6. Luận án tạo lập một di sản học thuật chuẩn mực trong chuyên ngành Cơ kỹ thuật, cung cấp công cụ tính toán giải tích tin cậy, chính xác và tiện dụng phục vụ trực tiếp cho công tác thiết kế, chế tạo các kết cấu vỏ mỏng chịu tải siêu hạn trong công nghiệp hàng không vũ trụ và công trình biển hiện đại.