Tổng quan về luận án

Nghiên cứu về cơ học kết cấu tiên tiến và vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) đóng vai trò then chốt trong sự phát triển của các ngành công nghiệp mũi nhọn như hàng không vũ trụ, kỹ thuật quân sự, công nghệ hạt nhân và công trình biển. Luận án tiến sĩ 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" do nghiên cứu sinh Phạm Minh Vương 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 Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (chuyên ngành Cơ học vật rắn, mã số: 9440109.02) là một công trình khoa học tiên phong, giải quyết toàn diện bài toán cơ học phi tuyến của kết cấu vỏ hai độ cong có dạng hình học phức tạp.

Vỏ trống (barrel shell) là kết cấu vỏ tròn xoay hai độ cong được tạo thành khi một cung tròn bán kính $a$ quay quanh trục đối xứng nằm trong cùng mặt phẳng, đặc trưng bởi bán kính xích đạo $R$ và chiều dài $L$. Khi bán kính lớn $a \to \infty$, vỏ trống suy biến thành vỏ trụ tròn tiêu chuẩn. Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các công trình trước đây của Stein và McElman (1965), Hutchinson (1968), Bích và cộng sự (2012, 2014) chủ yếu khảo sát vỏ trống thuần nhất hoặc vỏ FGM dựa trên lý thuyết vỏ cổ điển (Classical Shell Theory - CST). Lý thuyết CST bỏ qua biến dạng trượt ngang ($\gamma_{xz} = \gamma_{yz} = 0$), do đó dẫn đến sai số đáng kể khi đánh giá ứng xử của các kết cấu vỏ có độ dày trung bình hoặc vỏ composite phân lớp chịu tải trọng phức tạp. Sự thiếu vắng một khung phân tích giải tích hoàn chỉnh dựa trên lý thuyết biến dạng trượt bậc ba (Third-order Shear Deformation Theory - TSDT) cho vỏ trống FGM chính là điểm nghẽn lý thuyết mà luận án tập trung khắc phục.

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ẽ:

  1. $RQ_1$: Bằng cách nào xây dựng được hệ phương trình vi phân chủ đạo khép kín mô tả ứng xử phi tuyến của vỏ trống FGM dày có xét đến biến dạng trượt bậc cao và tương tác nền đàn hồi Pasternak?
  2. $RQ_2$: Tải trọng tới hạn tĩnh (áp lực ngoài $q_{cr}$, mô men xoắn $\tau_{cr}$, lực nén dọc trục $P_{cr}$, tải nhiệt $T_{cr}$) và đường cong cân bằng sau vồng (postbuckling) chịu chi phối như thế nào bởi tham số độ cong $a/R$, tỷ phần thể tích $k$, và hệ thống gân gia cường?
  3. $RQ_3$: Tải trọng tới hạn động ($P_{dcr}$) của vỏ trống FGM dưới tác dụng của lực nén phụ thuộc thời gian biến đổi ra sao theo tiêu chuẩn mất ổn định động Budiansky-Roth?
  4. $RQ_4$: Đặc trưng dao động phi tuyến và hiện tượng cộng hưởng của vỏ trống FGM chịu kích động điều hòa chịu ảnh hưởng như thế nào từ các tham số cơ - nhiệt - hình học?

Tương ứng với đó là các giả thuyết:

  • $H_1$: Việc áp dụng TSDT của Reddy và Liu sẽ loại bỏ nhu cầu sử dụng hệ số hiệu chỉnh trượt, phản ánh chính xác sự phân bố phi tuyến của ứng suất trượt ngang qua chiều dày vỏ.
  • $H_2$: Vỏ trống lồi ($a > 0$, độ cong Gauss dương) có độ cứng và khả năng chịu tải tới hạn vượt trội so với vỏ trụ ($a \to \infty$) và vỏ trống lõm ($a < 0$, độ cong Gauss âm).
  • $H_3$: Tải trọng động tới hạn $P_{dcr}$ phụ thuộc phi tuyến vào tốc độ gia tải $C$, trong đó quán tính kết cấu tạo ra hiệu ứng trễ làm tăng khả năng chịu tải tức thời.

Khung lý thuyết của luận án tích hợp lý thuyết biến dạng trượt bậc ba Reddy - Liu (1995), quan hệ động học phi tuyến hình học von Kármán, mô hình đồng nhất hóa vật liệu Voigt, mô hình nền đàn hồi hai tham số Pasternak và kỹ thuật san đều gân gia cường Lekhnitskii. Luận án mang lại đóng góp đột phá khi thiết lập thành công các công thức giải tích tường minh xác định tải vồng và tần số dao động tự do, tạo công cụ tính toán chính xác, tiết kiệm chi phí thử nghiệm thực nghiệm cho kỹ thuật thiết kế kết cấu FGM chịu môi trường nhiệt - cơ khắc nghiệt.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử cơ học của kết cấu vỏ trống bắt đầu từ các công trình đặt nền móng của Manuel Stein và John A. McElman (1965) tại NASA Langley Research Center về dao động và ổn định của vỏ trống thuần nhất chịu tải cơ học. Tiếp nối hướng đi này, John W. Hutchinson (1968) tại Đại học Harvard đã mở rộng phân tích cho vỏ trống thoải có độ không hoàn hảo hình học ban đầu. Các nghiên cứu kinh điển tiếp theo của Weingarten và cộng sự (1968), Stroud và Sykes (1969), Oyesanya (2008) chủ yếu tập trung tối ưu hóa trọng lượng và phân tích độ mất ổn định ngẫu nhiên của vỏ trống mỏng thuần nhất sử dụng CST.

Về dòng nghiên cứu vật liệu FGM, từ đề xuất ban đầu của các nhà khoa học Nhật Bản năm 1984, lý thuyết đồng nhất hóa và phân tích kết cấu FGM phát triển mạnh mẽ qua các công trình của Praveen và Reddy (1998), Shen (2009, 2011), Shariyat (2009). Đối với phân tích tấm vỏ FGM, các nghiên cứu trong nước của Đào Văn Dũng và Nguyễn Đình Đức (2012, 2014, 2016) đã tiên phong áp dụng CST và HSDT cho vỏ trụ, vỏ nón có gân gia cường trong môi trường nhiệt độ cao. Tuy nhiên, khi khảo sát riêng kết cấu vỏ trống FGM, các công trình của Bích và Ninh (2012, 2014), Thắng và Thời (2014) hoàn toàn dựa trên nền tảng lý thuyết vỏ cổ điển CST hoặc lý thuyết Donnell cải tiến, vốn bỏ qua hoàn toàn biến dạng trượt ngang.

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

  • Trường phái CST (Classical Shell Theory): Giữ nguyên giả thuyết Kirchhoff-Love, giả định đường pháp tuyến với mặt phẳng trung hòa luôn thẳng và vuông góc với mặt trung hòa sau khi biến dạng. Cách tiếp cận này giúp đơn giản hóa hệ phương trình vi phân chỉ còn 2 phương trình với 2 ẩn hàm ($w, F$), nhưng gây ra sự đánh giá quá cao (overestimation) độ cứng kết cấu, dẫn đến sai số tải trọng tới hạn lên tới 15% - 25% đối với vỏ có tỷ số $R/h < 30$.
  • Trường phái HSDT/TSDT (Higher-order Shear Deformation Theory): Đại diện bởi J.N. Reddy (1995, 2004), giả thiết trường chuyển vị có dạng hàm bậc ba của tọa độ chiều dày $z$. Lý thuyết này thỏa mãn điều kiện ứng suất tiếp bằng không tại hai bề mặt tự do của vỏ mà không cần dùng hệ số hiệu chỉnh trượt ($K_s = 5/6$ như trong FSDT của Mindlin/Reissner), nhưng dẫn đến hệ phương trình đạo hàm riêng 5 ẩn vô cùng phức tạp ($u, v, w, \phi_x, \phi_y$).

Luận án của Phạm Minh Vương định vị chính xác tại điểm giao thoa giữa độ chính xác cơ học cao cấp và tính khả thi của giải pháp giải tích. Bằng việc kết hợp hàm ứng suất Airy $F(x,y)$ với hai hàm góc xoay ($\phi_x, \phi_y$) và độ võng ($w$), tác giả đã rút gọn hệ 5 phương trình chuyển vị phức tạp về hệ 4 phương trình vi phân đạo hàm riêng chủ đạo. Khi so sánh với các nghiên cứu quốc tế như công trình của Huang và Han (2010) cho vỏ trụ FGM hoặc nghiên cứu của Ali và Hasan (2018) cho vỏ trống E-FGM, luận án thể hiện bước tiến vượt bậc khi khảo sát đồng thời cả 4 dạng tải trọng tĩnh (áp lực ngoài, mô men xoắn, lực nén dọc trục, tải nhiệt), kết cấu FGM sandwich, hệ thống gân gia cường trực giao và đáp ứng động lực học phi tuyến theo thời gian.

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

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

Luận án mang lại bốn đóng góp học thuật mang tính đột phá cho lý thuyết cơ học vật rắn:

  1. Mở rộng lý thuyết TSDT cho vỏ hai độ cong FGM: Luận án đã mở rộng thành công lý thuyết biến dạng trượt bậc ba của Reddy và Liu sang hệ tọa độ cong cho vỏ trống thoải ($R_0 = R - a(1 - \sin\phi) \approx R$ khi $\sin\phi \approx 1$). Trường chuyển vị được thiết lập chính xác: $$u^(x,y,z) = u\left(1 - \frac{z}{a}\right) + z\phi_x - \frac{4z^3}{3h^2}(w_{,x} + \phi_x)$$ $$v^(x,y,z) = v\left(1 - \frac{z}{R}\right) + z\phi_y - \frac{4z^3}{3h^2}(w_{,y} + \phi_y)$$ $$w^*(x,y,z) = w(x,y)$$ Biểu thức này loại trừ hoàn toàn các thành phần ứng suất trượt ngang giả tạo ở biên trên và biên dưới ($z = \pm h/2$).

  2. Thiết lập hệ phương trình vi phân đạo hàm riêng phi tuyến 4 ẩn: Tác giả đã dẫn xuất thành công hệ 4 phương trình vi phân phi tuyến bậc cao gồm 3 phương trình cân bằng và 1 phương trình tương thích biến dạng: $$\Delta^2 F = E_1 \left(w_{,xy}^2 - w_{,xx}w_{,yy} - \frac{w_{,yy}}{a} - \frac{w_{,xx}}{R}\right)$$ kết hợp cùng các phương trình liên hệ giữa lực cắt bậc cao ($R_x, R_y$), mô men bậc cao ($P_x, P_y, P_{xy}$) và các góc xoay pháp tuyến ($\phi_x, \phi_y$). Đây là một kết quả mới, chưa từng được công bố trong các tài liệu cơ học vỏ trước đây.

  3. Mô hình hóa giải tích vỏ trống FGM có gân gia cường: Áp dụng quan điểm san đều Lekhnitskii, luận án đã lượng hóa đóng góp độ cứng của $n_s$ gân dọc (chiều rộng $b_s$, chiều cao $h_s$) và $n_r$ gân vòng (chiều rộng $b_r$, chiều cao $h_r$) vào ma trận độ cứng tổng thể của kết cấu FGM, cho phép đánh giá chính xác tác động tương hỗ giữa cơ tính biến thiên liên tục và tính trực giao cục bộ của gân.

  4. Khái quát hóa điều kiện ổn định động theo tiêu chuẩn Budiansky-Roth: Luận án chứng minh một cách giải tích rằng điểm uốn (inflection point) trên đường cong độ võng cực đại theo thời gian $w_{max}(t)$ là tiêu chuẩn định lượng khách quan nhất để xác định tải trọng tới hạn động $P_{dcr}$, giải quyết triệt để tính mơ hồ của các tiêu chuẩn chuyển vị tới hạn truyền thống.

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

Khung phân tích của luận án tích hợp đồng bộ ba lý thuyết cơ học trụ cột:

  • Lý thuyết vật liệu P-FGM: Tỷ phần thể tích gốm $V_c(z) = \left(\frac{2z+h}{2h}\right)^k$ và kim loại $V_m(z) = 1 - V_c(z)$ phân bố theo quy luật lũy thừa qua chiều dày ($-h/2 \le z \le h/2$), với $k \in [0, \infty)$ là chỉ số phân bố. Mô hình Voigt tuyến tính hóa các đặc trưng cơ - nhiệt hữu hiệu: $$P_{eff}(z) = P_m + (P_c - P_m)\left(\frac{2z+h}{2h}\right)^k$$ với $P \in {E, \alpha, \rho, \kappa}$.
  • Mô hình nền tương tác Pasternak hai tham số: Áp lực phản lực của nền đàn hồi được mô hình hóa qua toán tử vi phân bậc hai: $$q_{sf} = K_1 w - K_2 \nabla^2 w = K_1 w - K_2 \left(\frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2}\right)$$ trong đó $K_1$ (N/m³) là độ cứng Winkler đặc trưng cho biến dạng nén thể tích và $K_2$ (N/m) là độ cứng lớp cắt Pasternak đặc trưng cho tương tác trượt ngang giữa các phần tử nền.
  • Tiếp cận nghiệm xấp xỉ dạng ba số hạng (Three-term displacement field): $$w(x,y) = W_0 + W_1 \sin(Mx)\sin(Ny) + W_2 \sin^2(Mx)$$ với $M = m\pi/L$, $N = n/R$, ($m, n \in \mathbb{N}^*$). Dạng nghiệm này cho phép mô tả chính xác dạng vồng phi đối xứng và hiệu ứng uốn cục bộ của vỏ dưới tác dụng của tải trọng áp lực ngoài và tải xoắn thuần túy, vượt qua hạn chế của dạng nghiệm một số hạng truyền thống.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan thực chứng (positivism) và lập trường nhận thức luận duy lý cơ học (rationalist mechanics). Thiết kế nghiên cứu là sự kết hợp chuẩn xác giữa phương pháp giải tích thuần túy (analytical approach) cho các bài toán tĩnh và phương pháp nửa giải tích (semi-analytical approach) cho các bài toán động lực học.

Quy mô khảo sát hình học và cấu trúc vật liệu được xác định cụ thể:

  • Mô hình vỏ trống lồi ($a > 0$) và vỏ trống lõm ($a < 0$) với bán kính cong $a/R \in [2, 10]$ và $a \to \infty$ (trường hợp suy biến vỏ trụ).
  • Tỷ số kích thước hình học: Chiều dài trên bán kính $L/R \in [1, 5]$; độ mỏng $R/h \in [20, 100]$ (từ vỏ dày vừa đến vỏ mỏng).
  • Cấu trúc vật liệu: Vỏ P-FGM đơn lớp, vỏ FGM sandwich ba lớp (lõi FGM mặt gốm/kim loại hoặc lõi thuần nhất mặt FGM với tỷ số bề dày $h_f/h \in [0.1, 0.4]$), và vỏ FGM gia cường gân trực giao.

Quy trình nghiên cứu rigorous

Quy trình giải tích được thực hiện theo 5 giai đoạn nghiêm ngặt:

  1. Thiết lập quan hệ biến dạng - chuyển vị và ứng suất - biến dạng: Áp dụng hệ thức von Kármán cho vỏ thoải và định luật Hooke tổng quát cho môi trường dị tính cơ - nhiệt.
  2. Khép kín hệ phương trình vi phân: Tích phân các thành phần ứng suất qua chiều dày để thu được các lực màng ($N_x, N_y, N_{xy}$), mô men uốn ($M_x, M_y, M_{xy}$), mô men bậc cao ($P_x, P_y, P_{xy}$) và lực cắt bậc cao ($Q_x, Q_y, R_x, R_y$). Đưa vào hàm ứng suất Airy $F(x,y)$ để thỏa mãn đồng nhất các phương trình cân bằng lực màng trong mặt phẳng.
  3. Xác định dạng nghiệm giải tích tương thích: Dựa trên dạng nghiệm giả định của độ võng $w(x,y)$, tiến hành thế vào phương trình tương thích biến dạng và các phương trình liên hệ góc xoay để giải ra dạng tường minh của $F(x,y)$, $\phi_x(x,y)$ và $\phi_y(x,y)$.
  4. Trực giao hóa Galerkin: Các điều kiện biên tựa đơn được thỏa mãn trực tiếp hoặc theo nghĩa tích phân trung bình: $$\int_0^{2\pi R} \int_0^L M_{x,x} , dx dy = 0, \quad \int_0^{2\pi R} \int_0^L N_{x,x} , dx dy = 0$$ Áp dụng phép chiếu trực giao Galerkin trên toàn miền xác định: $$\int_0^{2\pi R} \int_0^L \Phi(x,y) \sin(Mx)\sin(Ny) , dx dy = 0$$
  5. Thiết lập biểu thức tải tới hạn và phân tích phi tuyến: Biểu thức tải trọng tới hạn tại điểm rẽ nhánh được suy dẫn chính xác bằng cách lấy giới hạn khi biên độ võng không đều tiến về không: $$P_{cr} = \lim_{W_1 \to 0} P(W_1)$$

Data và phân tích

Các tham số vật liệu thành phần được trích xuất từ dữ liệu thực nghiệm chuẩn quốc tế:

  • Gốm Silicon Nitride ($Si_3N_4$): Mô đun đàn hồi $E_c = 322.4 \times 10^9 \text{ N/m}^2$, hệ số Poisson $\nu = 0.24$, hệ số giãn nở nhiệt $\alpha_c = 1.2 \times 10^{-6} \text{ K}^{-1}$, hệ số dẫn nhiệt $\kappa_c = 2.09 \text{ W/mK}$, khối lượng riêng $\rho_c = 3000 \text{ kg/m}^3$.
  • Kim loại Thép không gỉ (SUS304): Mô đun đàn hồi $E_m = 207.788 \times 10^9 \text{ N/m}^2$, hệ số Poisson $\nu = 0.24$, hệ số giãn nở nhiệt $\alpha_m = 15.3 \times 10^{-6} \text{ K}^{-1}$, hệ số dẫn nhiệt $\kappa_m = 12.08 \text{ W/mK}$, khối lượng riêng $\rho_m = 8166 \text{ kg/m}^3$.

Phương pháp số hóa và thuật toán: Toàn bộ hệ phương trình phi tuyến đại số và phương trình vi phân phi tuyến theo thời gian được lập trình giải trên môi trường tính toán biểu thức tượng trưng MATLAB/Maple. Đối với bài toán động lực học phi tuyến, phương pháp Runge-Kutta bậc 4-5 thích ứng thời gian được áp dụng để vẽ quỹ đạo nghiệm phase-plane và đường cong đáp ứng $w_{max}(t)$.

Độ tin cậy (validity) của mô hình được kiểm chứng chéo (cross-validation) bằng cách so sánh kết quả khi suy biến ($a \to \infty$) với các công bố uy tín quốc tế của Shen (2009) cho vỏ trụ FGM, Huang và Han (2010), cũng như các nghiệm giải tích của Đào Văn Dũng (2012). Sai số giữa mô hình TSDT của luận án và các tài liệu tham khảo trong trường hợp suy biến đều nhỏ hơn 1.5%, chứng minh độ tin cậy tuyệt đối của thuật toán.

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

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

  1. Hiệu ứng gia cường hình học của độ cong Gauss dương: Nghiên cứu chỉ ra rằng tỷ số bán kính $a/R$ có ảnh hưởng mang tính quyết định đến khả năng chịu lực. Vỏ trống lồi ($a > 0$) có tải trọng tới hạn tĩnh ($q_{cr}, \tau_{cr}, P_{cr}$) và tải tới hạn động ($P_{dcr}$) cao hơn đáng kể so với vỏ trụ tròn ($a \to \infty$) và vượt trội hoàn toàn so với vỏ trống lõm ($a < 0$). Khi tỷ số $a/R$ giảm từ 10 xuống 2 (độ lồi tăng lên), tải trọng nén tới hạn $P_{cr}$ của vỏ trống lồi tăng tới 34.8%, trong khi vỏ trống lõm bị suy giảm khả năng chịu tải tương ứng do xuất hiện vùng tập trung ứng suất uốn cục bộ.

  2. Quy luật suy giảm độ cứng theo chỉ số tỷ phần thể tích $k$: Khi chỉ số $k$ tăng từ 0 (thuần gốm $Si_3N_4$) lên $k = 0.5, 1, 5$ và tiến tới $\infty$ (thuần kim loại SUS304), tải trọng tới hạn tĩnh và tần số dao động tự do tuyến tính $\omega$ suy giảm đơn điệu. Cụ thể, tải trọng tới hạn áp lực ngoài $q_{cr}$ của vỏ thuần gốm cao hơn khoảng 55.2% so với vỏ thuần kim loại. Nguyên nhân vật lý là do gốm có mô đun đàn hồi $E_c = 322.4 \text{ GPa}$ lớn hơn 1.55 lần so với thép không gỉ ($E_m = 207.79 \text{ GPa}$).

  3. Ảnh hưởng vượt trội của độ cứng lớp cắt Pasternak ($K_2$) so với Winkler ($K_1$): Khảo sát số chứng minh rằng việc tăng tham số nền $K_2$ làm tăng tải trọng tới hạn tĩnh và động hiệu quả hơn gấp nhiều lần so với việc tăng tham số $K_1$. Lớp cắt $K_2$ tạo ra liên kết chống trượt ngang giữa các phần tử vỏ, làm biến đổi mode vồng từ dạng sóng dài sang dạng sóng ngắn cục bộ ($m, n$ tăng lên), từ đó triệt tiêu sự phát triển mất ổn định vồng lớn sau tới hạn.

  4. Hiệu ứng trễ quán tính trong mất ổn định động: 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 $P(t) = C \cdot t$, tải trọng tới hạn động $P_{dcr}$ luôn lớn hơn tải trọng tới hạn tĩnh $P_{cr}$ ($P_{dcr} > P_{cr}$). Tốc độ đặt tải $C$ (Pa/s) càng lớn thì giá trị $P_{dcr}$ đạt được càng cao. Hiện tượng này được giải thích bằng cơ chế trễ động lực học (dynamic delay): vận tốc biến dạng cao khiến kết cấu không kịp tích lũy biến dạng uốn ở mode cơ bản trước khi vượt qua ngưỡng cân bằng tĩnh.

  5. Tối ưu hóa cấu trúc sandwich và hệ thống gân gia cường: Đối với vỏ trống FGM sandwich, cấu trúc có lớp lõi FGM và hai bề mặt thuần gốm cho khả năng kháng vồng nhiệt $\Delta T_{cr}$ tối ưu nhất. Đối với vỏ có gân, việc bố trí gân vòng kết hợp gân dọc giúp tăng tải trọng tới hạn $q_{cr}$ lên tới 2.5 - 3.8 lần so với vỏ trơn cùng khối lượng vật liệu, đồng thời chuyển đường cong cân bằng sau vồng từ dạng rẽ nhánh không ổn định (unstable snap-through) sang dạng ổn định tăng tải.

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định giới hạn ứng dụng bắt buộc của lý thuyết trượt bậc cao TSDT đối với kết cấu vỏ hai độ cong FGM có chiều dày trung bình ($R/h \le 50$), đặt nền tảng giải tích chuẩn mực để hiệu chỉnh các phần mềm phần tử hữu hạn thương mại (như ABAQUS, ANSYS).
  • Về phương pháp luận: Đóng góp quy trình giải tích nửa tường minh kết hợp giữa hàm ứng suất Airy và phép chiếu Galerkin đa số hạng, mở rộng năng lực giải quyết các bài toán biên phi tuyến bậc cao trong cơ học giải tích.
  • Về ứng dụng thực tiễn: Cung cấp bộ công thức tường minh tính toán tải trọng giới hạn cho các kỹ sư thiết kế thân tàu ngầm mini (chịu áp lực ngoài), vỏ buồng đốt động cơ tên lửa và tuabin khí (chịu nhiệt - nén kết hợp), ống dẫn dầu áp lực cao và vỏ máy bay tàng hình.

Limitations và Future Research

Mặc dù đạt được những kết quả học thuật xuất sắc, luận án vẫn tồn tại một số giới hạn nghiên cứu (boundary conditions) cần được nhìn nhận khách quan:

  1. Giả thiết vỏ thoải: Luận án sử dụng giả thiết độ cong thoải của Stein và McElman ($\sin\phi \approx 1, \cos\phi \approx 0$), do đó mô hình chỉ phản ánh chính xác nhất cho các vỏ có góc mở đáy không quá lớn.
  2. Quy luật vật liệu 1D-FGM: Cơ tính chỉ biến thiên theo một phương chiều dày $z$ (1D P-FGM), chưa xét đến các quy luật biến thiên không gian hai hoặc ba chiều (2D/3D-FGM) hay vật liệu FGM xốp có lỗ rỗng (porosities).
  3. Điều kiện biên: Phân tích chủ yếu giới hạn ở điều kiện biên tựa đơn bốn cạnh (simply supported), chưa mở rộng giải tích tường minh cho điều kiện biên ngàm chặt (clamped) hoặc biên tự do phức tạp.
  4. Phụ thuộc nhiệt độ: Một số bài toán tĩnh chưa tích hợp đầy đủ sự suy giảm mô đun đàn hồi theo nhiệt độ phi tuyến phi tuyến tính $E(T)$.

Chương trình nghiên cứu mở rộng trong giai đoạn tiếp theo (5-10 năm tới) bao gồm:

  • Mở rộng hệ phương trình TSDT cho vỏ trống sâu (deep barrel shells) không sử dụng giả thiết thoải.
  • Nghiên cứu vỏ trống làm từ vật liệu composite tiên tiến mới: Polyme gia cường ống nano carbon (FG-CNTRC) hoặc tấm nano graphene (FG-GRC).
  • Phân tích động lực học kết cấu vỏ trống FGM tương tác thủy động lực (chứa chất lỏng hoặc ngâm trong chất lỏng vô hạn).
  • Ứng dụng phương pháp cầu phương vi phân (Differential Quadrature Method - DQM) để giải bài toán với các điều kiện biên hỗn hợp tùy ý.

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

Công trình của Phạm Minh Vương tạo ra tác động học thuật và công nghệ sâu rộng:

  • Tác động học thuật: Các bài báo khoa học xuất bản từ luận án trên các tạp chí quốc tế thuộc danh mục ISI/Scopus uy tín (như Thin-Walled Structures, Composite Structures, Journal of Sandwich Structures & Materials) đóng góp nguồn trích dẫn học thuật chất lượng cao, khẳng định vị thế của trường phái Cơ học tính toán Việt Nam (đứng đầu bởi nhóm nghiên cứu của GS. TSKH Nguyễn Đình Đức) trên trường quốc tế.
  • Chuyển giao công nghiệp & Quốc phòng: Bộ công thức giải tích trong luận án cho phép xây dựng các module phần mềm chuyên dụng tính toán nhanh độ bền kết cấu, giúp các kỹ sư hàng không vũ trụ và công nghiệp quốc phòng tối ưu hóa chiều dày và tỷ lệ gốm/kim loại mà không cần chạy mô phỏng FEM tốn kém hàng trăm giờ tính toán.
  • Lợi ích kinh tế - xã hội: Giảm thiểu trọng lượng kết cấu FGM từ 20% - 30% so với kết cấu hợp kim truyền thống, giúp tiết kiệm nhiên liệu trong vận tải hàng không và nâng cao độ an toàn chịu nhiệt trong các lò phản ứng năng lượng hạt nhân thế hệ mới.

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

  • Nghiên cứu sinh và Giảng viên Cơ học: Tiếp cận hệ thống phương trình giải tích chuẩn mực, làm tài liệu tham khảo cao cấp để phát triển các đề tài luận án về kết cấu composite thông minh áp điện, vật liệu biến tính đa pha.
  • Kỹ sư R&D Công nghiệp Hàng không & Vũ trụ: Ứng dụng ngay các công thức hiển để định kích thước nắp mũi tàu vũ trụ, cánh turbine khí, tối ưu hóa tỷ phần thể tích $k$ đạt độ bền nhiệt cơ lý tưởng.
  • Ngành Công nghiệp Chế tạo Quốc phòng: Ứng dụng mô hình vỏ FGM có gân chịu áp lực ngoài để chế tạo các khoang áp lực của thiết bị lặn không người lái (UUV) và thân xe bọc thép chịu nhiệ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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độ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 ba (TSDT) của Reddy và Liu (1995) cho kết cấu vỏ tròn xoay hai độ cong có dạng vỏ trống FGM thoải. Bằng việc kết hợp hàm ứng suất Airy $F(x,y)$, luận án đã rút gọn hệ 5 phương trình đạo hàm riêng phi tuyến phức tạp về hệ 4 phương trình khép kín với 4 ẩn hàm ($w, F, \phi_x, \phi_y$), loại bỏ hoàn toàn hệ số hiệu chỉnh trượt $K_s$ và mô tả chính xác ứng suất tiếp phân bố phi tuyến parabol qua chiều dày.

2. Đột phá về phương pháp luận so với ít nhất 2 nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Stein và McElman (1965) và Bích và Ninh (2014) vốn sử dụng lý thuyết vỏ cổ điển CST (bỏ qua trượt ngang) với nghiệm 1 số hạng, luận án đã đột phá khi:

  • Sử dụng lý thuyết trượt bậc cao TSDT kết hợp nghiệm trường chuyển vị 3 số hạng $w(x,y) = W_0 + W_1\sin(Mx)\sin(Ny) + W_2\sin^2(Mx)$, cho phép nắm bắt chính xác sự biến dạng bất đối xứng và tương tác xoắn - uốn phức tạp.
  • So với phương pháp năng lượng Ritz số trị thuần túy của Huang và Han (2010) trên vỏ trụ FGM, luận án giải quyết thành công bài toán vỏ trống hai độ cong bằng giải pháp giải tích tường minh giải được bằng phương pháp trực giao Galerkin.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì và cơ sở vật lý giải thích?

Phát hiện bất ngờ nhất là ứng xử cơ học dị biệt của vỏ trống lõm ($a < 0$): Mặc dù có cùng diện tích bề mặt và thể tích vật liệu với vỏ trống lồi ($a > 0$), vỏ trống lõm lại có tải trọng tới hạn tĩnh và động sụt giảm nghiêm trọng, đồng thời đường cong biên độ dao động xuất hiện các vùng mất ổn định hỗn loạn sớm hơn nhiều. Cơ sở vật lý là do độ cong Gauss âm ($K = -1/(|a|R) < 0$) tạo ra trạng thái ứng suất màng dạng hyperbol, gây ra sự tập trung ứng suất cắt cục bộ mạnh mẽ, triệt tiêu khả năng chịu nén tổng thể của kết cấu.

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

Có. Luận án cung cấp hệ thống công thức giải tích hoàn chỉnh, chi tiết từ các hệ số tích phân độ cứng $E_i$ ($i = 1-7$), các tham số tương tác nhiệt $\Phi_1, \Phi_2, \Phi_4$, ma trận hệ số Galerkin $A_i, B_i, C_i, D_i$, cùng các bước thiết lập nghiệm giải tích và thông số vật liệu đầu vào ($Si_3N_4$, SUS304). Bất kỳ nhà nghiên cứu nào cũng có thể tái lập 100% các kết quả khảo sát số thông qua các phần mềm giải tích ký hiệu như Maple hoặc MATLAB.

5. Khung chương trình nghiên cứu 10 năm tới được vạch ra như thế nào?

Luận án định hình chương trình nghiên cứu dài hạn tập trung vào:

  1. Phát triển lý thuyết vỏ trượt bậc cao phi cục bộ (Nonlocal HSDT) cho vỏ trống FGM kích thước micro/nano.
  2. Tích hợp cơ chế lan truyền vết nứt nhiệt - cơ phi tuyến (thermo-mechanical fracture mechanics) trong vỏ FGM.
  3. Nghiên cứu ứng xử biến dạng lớn của vỏ trống FGM áp điện - từ - nhiệt thông minh thích ứng (Smart Magneto-Electro-Elastic FGM shells).

Kết luận

Luận án tiến sĩ của Phạm Minh Vương là một công trình khoa học xuất sắc, mẫu mực về phương pháp nghiên cứu giải tích cơ học vật rắn.

Các đóng góp cốt lõi được đúc kết cụ thể:

  1. Thiết lập thành công hệ 4 phương trình vi phân đạo hàm riêng chủ đạo mô tả ứng xử phi tuyến của vỏ trống FGM hai độ cong dựa trên lý thuyết biến dạng trượt bậc ba Reddy - Liu.
  2. Dẫn xuất bộ công thức giải tích tường minh xác định tải trọng tới hạn tĩnh ($q_{cr}, \tau_{cr}, P_{cr}, T_{cr}$) và quan hệ tải - độ võng sau vồng cho vỏ trống FGM trơn, FGM sandwich và FGM có gân gia cường trên nền Pasternak.
  3. Xác định chính xác tải trọng tới hạn động $P_{dcr}$ dưới tác dụng của tải nén biến thiên thời gian theo tiêu chuẩn Budiansky-Roth, chứng minh bản chất của hiệu ứng trễ quán tính.
  4. Phân tích toàn diện động lực học phi tuyến và dao động cưỡng bức điều hòa, làm sáng tỏ ảnh hưởng của các tham số hình học ($a/R, L/R, R/h$) và quy luật vật liệu $k$.
  5. Định lượng hóa vai trò vượt trội của độ cong lồi ($a > 0$) và độ cứng lớp cắt Pasternak ($K_2$) trong việc nâng cao độ bền ổn định của kết cấu.

Công trình không chỉ hoàn thiện bức tranh lý thuyết về cơ học vỏ hai độ cong FGM mà còn mở ra các nhánh nghiên cứu mới về vật liệu nano composite đa chức năng, đóng góp trực tiếp vào sự nghiệp phát triển khoa học công nghệ vật liệu và kỹ thuật công trình hiện đại.