Tổng quan về luận án

Nghiên cứu động lực học kết cấu hiện đại đang chứng kiến bước chuyển biến sâu sắc nhờ sự ra đời của các thế hệ vật liệu tiên tiến, đặc biệt là vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM). Được phát triển từ năm 1984 tại Viện Nghiên cứu Vật liệu Sendai (Nhật Bản), FGM là giải pháp kỹ thuật mang tính đột phá nhằm giải quyết triệt để hiện tượng tập trung ứng suất tiếp và nguy cơ tách lớp vốn là nhược điểm chí mạng của composite nhiều lớp truyền thống. Bằng cách thiết lập quy luật biến đổi trơn liên tục thành phần thể tích giữa hai pha gốm (chịu nhiệt, chống mài mòn) và kim loại (dẻo dai, chịu lực), kết cấu FGM duy trì trọn vẹn độ bền cấu trúc trong các điều kiện môi trường tải trọng nhiệt - cơ khắc nghiệt. Trong các lĩnh vực công nghệ cao như công nghiệp hàng không - vũ trụ (vỏ động cơ phản lực, thân tên lửa đẩy nhiên liệu hạt nhân), kỹ thuật quân sự hải quân (vỏ chịu áp tàu ngầm lặn sâu), và năng lượng hạt nhân, các dạng vỏ liên hợp tròn xoay đối xứng trục (vỏ nón, vỏ trụ, vành tròn, vỏ dạng bậc, và vỏ có gân gia cường) thường xuyên làm việc trong môi trường bao quanh bởi nền đàn hồi.

Khoảng trống học thuật (Research Gap) lớn tồn tại trong cơ học kết cấu FGM là sự thiếu hụt các mô hình tính toán số tin cậy cho bài toán dao động tự do của các hệ vỏ liên hợp tròn xoay phức tạp tương tác với nền đàn hồi hai hệ số. Phần lớn các công bố quốc tế trước đây của Sofiyev & Schnack (2003) hay Sheng & Wang (2008) chủ yếu tập trung vào các kết cấu vỏ đơn lẻ (vỏ trụ hoặc vỏ nón thuần túy) dựa trên lý thuyết vỏ mỏng Donnell hoặc phương pháp Galerkin bán giải tích. Trong khi đó, các phương pháp số truyền thống như Phương pháp Phần tử Hữu hạn (FEM) áp dụng trường chuyển vị tĩnh chịu chi phí tính toán cực lớn và dễ bị mất độ chính xác khi phân tích trong miền tần số dao động riêng trung bình và cao do hiệu ứng khóa cắt (shear locking) và đòi hỏi chia lưới cực mịn. Ngược lại, việc giải tích phương trình vi phân cho kết cấu ghép nối nhiều bậc hay phân nhánh chữ T gặp rào cản toán học nghiêm ngặt.

Luận án tiến sĩ cơ học của tác giả Lê Quang Vinh (Mã số: 9440109), được thực hiện dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Mạnh Cường và GS.TSKH Nguyễn Đông Anh tại Trường Đại học Bách khoa Hà Nội, đã giải quyết toàn diện bài toán cấp thiết này. Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập tường minh:

  1. Câu hỏi nghiên cứu 1: Làm thế nào để thiết lập một mô hình phần tử liên tục (Continuous Element Method - PTLT) tổng quát dựa trên lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin (FSDT) có thể mô tả chính xác ứng xử dao động của vỏ tròn xoay FGM có xét tương tác đầy đủ với nền đàn hồi Winkler-Pasternak?
  2. Câu hỏi nghiên cứu 2: Cơ chế giải thuật nào cho phép ghép nối linh hoạt các ma trận độ cứng động lực (Dynamic Stiffness Matrix - DSM) của các phân đoạn vỏ đơn (trụ, nón, vành) thành các hệ kết cấu liên hợp phức tạp (nón-trụ, nón-vành-trụ, vỏ trụ bậc và vỏ có gân gia cường dạng vành)?
  3. Giả thuyết khoa học 1: Việc chọn các hàm dạng là trường chuyển vị động chính xác từ việc giải hệ phương trình trạng thái sẽ giúp phương pháp PTLT đạt độ chính xác vô hạn về bậc tự do không gian, sử dụng số lượng phần tử tối thiểu mà không cần chia lưới hội tụ.
  4. Giả thuyết khoa học 2: Phương pháp đường cong đáp ứng tần số (Frequency Response Function) kết hợp kích động xung lực ảo sẽ cô lập hoàn toàn các tần số dao động riêng cơ bản và bậc cao mà không gặp hiện tượng kỳ dị toán học của định thức ma trận.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin, nguyên lý biến phân Hamilton, mô hình nền đàn hồi hai hệ số Winkler-Pasternak và mô hình biến thiên cơ tính 4 tham số $FGM(a/b/c/p)$. Đóng góp đột phá của luận án nằm ở việc phát triển hoàn chỉnh bộ công cụ tính toán số tự động trên nền tảng MATLAB, giải quyết thành công phổ tần số riêng từ bậc 1 đến bậc 10 cho hàng loạt kết cấu vỏ liên hợp phức tạp với thời gian tính toán giảm hơn 80% so với FEM truyền thống, đồng thời mở rộng phạm vi khảo sát cơ học kết cấu FGM lên dải tần số cao.

Phạm vi nghiên cứu bao hàm các tổ hợp vật liệu FGM điển hình: Silicon Nitride/Thép không gỉ ($Si_3N_4/SUS304$), Oxit Zirconi/Hợp kim Titan ($ZrO_2/Ti6Al4V$), Oxit Nhôm/Nhôm ($Al_2O_3/Al$), tương tác với nền đàn hồi có hệ số Winkler $k_w = 0 \div 5\times 10^7\text{ N/m}^3$ và hệ số trượt Pasternak $k_p = 0 \div 5\times 10^6\text{ N/m}$, dưới đầy đủ các tổ hợp điều kiện biên khắt khe (Ngàm-Ngàm C-C, Ngàm-Tự do C-F, Tự do-Ngàm F-C, Tựa-Tựa S-S).

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận ba dòng chảy nghiên cứu chính về cơ học kết cấu FGM và dao động đàn hồi:

[Dòng chảy 1: Phân tích Giải tích & Bán Giải tích Vỏ FGM Đơn lẻ]
  (Sheng & Wang, 2008; Sofiyev, 2003-2010; Shen, 2009)
  --> Ưu điểm: Nghiệm tường minh cho vỏ mỏng đơn
  --> Nhược điểm: Bế tắc trước kết cấu vỏ liên hợp ghép nối nhiều đoạn

[Dòng chảy 2: Phương pháp Phần tử Hữu hạn (FEM) cho Kết cấu Ghép nối]
  (Patel et al., 2006; Kamat, 2013; Qu et al., 2013)
  --> Ưu điểm: Xử lý được hình học phức tạp
  --> Nhược điểm: Sai số lớn ở dải tần số cao, phụ thuộc mật độ lưới

[Dòng chảy 3: Phương pháp Ma trận Độ cứng Động lực & PTLT]
  (Kalnins, 1964; Casimir, 2005; Nguyen Manh Cuong, 2013)
  --> Bước đột phá của Luận án: Mở rộng PTLT sang hệ vỏ liên hợp FGM
      kết nối tiếp xúc nền Pasternak & thuật toán ghép nối chữ T

Dòng chảy thứ nhất tập trung vào các lời giải giải tích và bán giải tích cho vỏ FGM đơn lẻ. Sheng & Wang (2008) đã sử dụng nguyên lý Hamilton và lý thuyết phi tuyến von Kármán để khảo sát dao động nhiệt của vỏ trụ FGM. Sofiyev & Kuruoglu (2008) cùng Sofiyev & Schnack (2003) nghiên cứu dao động của vỏ nón FGM trên nền đàn hồi hai hệ số thông qua phương pháp Galerkin và lý thuyết Donnell. Tuy nhiên, các cách tiếp cận này bị giới hạn nghiêm ngặt trong phạm vi hình học trơn đơn giản, không thể mở rộng cho các hệ vỏ có bậc thay đổi chiều dày đột ngột hay các kết cấu vỏ liên hợp ghép nối.

Dòng chảy thứ hai ứng dụng phương pháp phần tử hữu hạn (FEM) và các phương pháp biến phân xấp xỉ cho vỏ ghép nối. Caresta & Kessissoglou (2010) đã phân tích dao động tự do của vỏ nón-trụ ghép nối đẳng hướng bằng cách biểu diễn chuyển vị vỏ trụ dạng sóng và vỏ nón bằng chuỗi bậc cao. Patel và cộng sự (2006) sử dụng FEM để tính dao động vỏ nón-trụ composite trực hướng, trong khi Qu và cộng sự (2013) phát triển phương pháp biến phân tổng quát cho hệ nón-trụ-vành-cầu. Điểm hạn chế căn bản của dòng nghiên cứu này là sự phụ thuộc vào hàm dạng chuyển vị tĩnh, khiến kích thước ma trận hệ thống tăng vọt và gặp hiện tượng suy biến độ chính xác khi khảo sát ở dải tần số trung bình và cao.

Dòng chảy thứ ba bắt nguồn từ phương pháp ma trận truyền và ma trận độ cứng động lực (Dynamic Stiffness Matrix) khởi xướng bởi Kalnins (1964), Casimir (2005) cho dầm và tấm kim loại dày, sau đó được phát triển bởi Trần Ích Thịnh và Nguyễn Mạnh Cường (2013) cho vỏ composite chứa chất lỏng.

Về mặt học thuật, tồn tại cuộc tranh luận kéo dài giữa hai trường phái:

  • Trường phái cổ điển (Classical Shell Theory - CST) cho rằng việc bỏ qua biến dạng cắt ngang là chấp nhận được đối với kết cấu mỏng nhằm đơn giản hóa tính toán vi phân.
  • Trường phái biến dạng cắt bậc cao/bậc nhất (FSDT/HSDT), đại diện bởi Tornabene (2009) và Shen (2009), chứng minh rằng đối với vật liệu không đồng nhất FGM có độ cứng biến đổi phi đối xứng qua chiều dày, ứng suất cắt ngang $\tau_{xz}, \tau_{\theta z}$ đóng vai trò sống còn trong việc xác định chính xác phổ tần số dao động.

Định vị nghiên cứu của luận án: Luận án định vị chính xác tại điểm giao thoa giữa lý thuyết FSDT Reissner-Mindlin và phương pháp PTLT. So sánh trực tiếp với nghiên cứu của Caresta & Kessissoglou (2010) (chỉ giải cho kết cấu đẳng hướng đồng nhất không có nền đàn hồi) và công trình của Qu và cộng sự (2013) (dùng phương pháp biến phân với hàm Ritz xấp xỉ), luận án đã tạo bước tiến vượt bậc khi thiết lập thành công ma trận độ cứng động lực chính xác cho từng phân đoạn vỏ FGM có cơ tính biến thiên 4 tham số, tích hợp nền Pasternak 2 hệ số, và đề xuất thuật toán ghép nối phần tử phân nhánh chữ T độc đáo.

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

Đóng góp cho lý thuyết cơ học kết cấu

Luận án mở rộng trực tiếp lý thuyết vỏ biến dạng cắt bậc nhất Reissner-Mindlin thông qua việc tái cấu trúc hệ phương trình quan hệ cơ-lý đối với môi trường vật liệu FGM dị hướng biến đổi liên tục:

  • Nguyên văn khẳng định từ luận án: "Vật liệu có cơ lý tính biến thiên với tên quốc tế là Functionally Graded Material (FGM) được nghiên cứu và phát triển lần đầu tiên bởi một nhóm các nhà khoa học ở viện Sendai của Nhật Bản vào năm 1984", trong đó "Các thông số mô đun đàn hồi $E(z)$, khối lượng riêng $\rho(z)$ và hệ số Poisson $\mu(z)$ được giả thiết là biến thiên đều đặn theo chiều dày vỏ".
  • Mở rộng quy luật phân bố thể tích 4 tham số $FGM(a/b/c/p)$: $$V_c(z) = \left[ 1 - a\left(\frac{1}{2} \pm \frac{z}{h}\right) + b\left(\frac{1}{2} \pm \frac{z}{h}\right)^c \right]^p$$ Cho phép mô tả chính xác cả cấu hình FGM giàu gốm ở mặt ngoài ($FGM_I$) lẫn giàu gốm ở mặt trong ($FGM_{II}$), khắc phục sự gò bó của hàm lũy thừa 1 tham số cổ điển.
  • Thiết lập tường minh ma trận độ cứng liên hợp màng-uốn-xoắn $[B]$ bên cạnh ma trận độ cứng màng $[A]$, độ cứng uốn $[D]$ và độ cứng cắt $[F]$: $$(A_{ij}, B_{ij}, D_{ij}) = \int_{-h/2}^{h/2} Q_{ij}(z)(1, z, z^2)dz, \quad F_{ij} = \int_{-h/2}^{h/2} Q_{ij}(z)dz$$ Sự hiện diện của ma trận tương tác $[B]$ chứng minh hiện tượng ghép đôi uốn-kéo tất yếu phát sinh do tính bất đối xứng cơ tính qua mặt trung bình của vỏ FGM.
+-------------------------------------------------------------------------+
|                  HỆ KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN                 |
+-------------------------------------------------------------------------+

Khung phân tích độc đáo và điều kiện biên

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng:

  1. Lý thuyết cơ học vỏ FSDT: Đảm bảo tính toán chính xác ứng xử cắt ngang mà không cần giả thiết pháp tuyến thẳng không biến dạng của Kirchhoff-Love.
  2. Lý thuyết đàn hồi môi trường Pasternak: Mô tả lực cản pháp tuyến $k_w w$ và biến dạng cắt tương tác giữa các lò xo $-k_p \nabla^2 w$, phản ánh cơ chế tương tác đất nền/môi trường bao quanh tiệm cận thực tế hơn hẳn mô hình Winkler 1 tham số thuần túy.
  3. Phương pháp phần tử liên tục PTLT: Luận án dẫn giải khẳng định: "Phương pháp độ cứng động là phương pháp PTHH, trong đó các hàm dạng được chọn là trường chuyển vị động của phần tử".

Điều kiện biên được chuẩn hóa khắt khe tại hai đầu biên $x = 0$ và $x = L$:

  • Biên ngàm (Clamped - C): $u_0 = v_0 = w_0 = \varphi_x = \varphi_\theta = 0$
  • Biên tự do (Free - F): $N_x = N_{x\theta} = Q_x = M_x = M_{x\theta} = 0$
  • Biên tựa bản lề (Simply supported - S): $v_0 = w_0 = \varphi_\theta = N_x = M_x = 0$

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

Thiết kế nghiên cứu và Bản chất phương pháp

Thiết kế nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng duy lý (Positivism) của cơ học tính toán lý thuyết kết hợp giải pháp số học chính xác cao. Quy trình nghiên cứu chuyển tiếp qua 3 cấp độ:

  • Cấp độ vi mô: Mô hình hóa cấu trúc pha vật liệu và tính toán giải tích các hệ số ma trận độ cứng màng, uốn, xoắn, cắt qua tích phân chiều dày.
  • Cấp độ phần tử đơn liên tục: Thiết lập phương trình vi phân biến thiên bậc nhất dọc theo đường sinh $x$ cho từng phân đoạn hình học (vỏ nón, trụ, vành tròn).
  • Cấp độ hệ thống kết cấu liên hợp: Sử dụng điều kiện cân bằng lực và liên tục chuyển vị tại các mặt tiếp giáp liên kết để hợp nhất các ma trận phần tử thành ma trận độ cứng động lực toàn cục $[K(\omega)]$.

Quy trình toán học và giải thuật tính toán

Quy trình tính toán số được chuẩn hóa gồm 9 bước nghiêm ngặt:

  1. Xây dựng véc-tơ trạng thái gồm 10 thành phần lực và chuyển vị: $$y(x,\theta,z,\omega,t) = {u_0, v_0, w_0, \varphi_x, \varphi_\theta, N_x, N_{x\theta}, Q_x, M_x, M_{x\theta}}^T$$
  2. Phân tách biến điều hòa theo chu vi vỏ với số sóng $m$ ($m=1, 2, 3...$): $${u, w, \varphi_x, N_x, Q_x, M_x} = {u_m(x), w_m(x), \varphi_{xm}(x), N_{xm}(x), Q_{xm}(x), M_{xm}(x)}\cos(m\theta)e^{i\omega t}$$ $${v, \varphi_\theta, N_{x\theta}, M_{x\theta}} = {v_m(x), \varphi_{\theta m}(x), N_{x\theta m}(x), M_{x\theta m}(x)}\sin(m\theta)e^{i\omega t}$$
  3. Chuyển đổi hệ phương trình cân bằng vi phân thành phương trình đạo hàm ma trận cấp 1: $$\frac{d{y_m(x,\omega)}}{dx} = [A_m(x,\omega)]{y_m(x,\omega)}$$ Trong đó ma trận trạng thái $[A_m]$ có kích thước $10\times 10$ chứa đầy đủ các thông số hình học góc nón $\alpha$, bán kính biến thiên $R(x) = R_1 + x\sin\alpha$, hệ số nền $k_w, k_p$, và các hằng số tích phân quán tính $I_0, I_1, I_2 = \int_{-h/2}^{h/2} \rho(z)[1, z, z^2]dz$.
  4. Tích phân số ma trận trạng thái dọc theo đường sinh để xác định ma trận truyền $[T(x,\omega)]$: $${y(x,\omega)} = [T(x,\omega)]{y(x_0,\omega)}$$
  5. Chuyển đổi ma trận truyền thành ma trận độ cứng động lực $[K(\omega)]_e$ của từng phần tử liên tục: $$\begin{Bmatrix} F(0) \ F(L) \end{Bmatrix} = [K(\omega)]_e \begin{Bmatrix} d(0) \ d(L) \end{Bmatrix}$$
  6. Ứng dụng giải thuật ghép nối nối tiếp cho vỏ bậc và ghép nối song song dạng chữ T cho vỏ trụ có gân gia cường dạng vành.
+--------------------------------------------------------------------------------+
|           SƠ ĐỒ GHÉP NỐI MA TRẬN ĐỘ CỨNG ĐỘNG LỰC DẠNG PHÂN NHÁNH CHỮ T        |
+--------------------------------------------------------------------------------+
                                       
                      Phân đoạn Vỏ Trụ 1: [K_1(\omega)]
                      Phân đoạn Vỏ Trụ 2: [K_2(\omega)]
             Ma trận Toàn cục: [K_global(\omega)] = [K_1] + [K_2] + [K_ring]
  1. Xác định tần số dao động tự do thông qua Phương pháp Đường cong Đáp ứng (Response Curve Method): Đặt xung kích thích ảo $Q_x(L) = \sum_{m=1}^\infty \frac{1}{\pi R}\cos(m\theta)e^{i\omega t}$ tại điểm biên khảo sát $M(x=L)$. Tại các tần số riêng $\omega = \omega_n$, hiện tượng cộng hưởng xảy ra khiến biên độ chuyển vị pháp tuyến $w$ đạt cực đại tiệm cận vô cùng. Hoành độ các đỉnh đỉnh đáp ứng chính là phổ tần số tự nhiên của kết cấu.

Đặc tính thông số vật liệu và độ tin cậy mô hình

Các hằng số cơ lý của các vật liệu thành phần trích xuất chính xác từ nghiên cứu:

  • Thép không gỉ (SUS304): $E_m = 207.7\text{ GPa}, \mu_m = 0.3177, \rho_m = 8166\text{ kg/m}^3, \kappa = 15.379\text{ W/mK}$
  • Silicon Nitride ($Si_3N_4$): $E_c = 322.2\text{ GPa}, \mu_c = 0.24, \rho_c = 2370\text{ kg/m}^3, \kappa = 13.723\text{ W/mK}$
  • Nhôm (Al): $E_m = 70.0\text{ GPa}, \mu_m = 0.30, \rho_m = 2707\text{ kg/m}^3$
  • Oxit Nhôm ($Al_2O_3$): $E_c = 380.0\text{ GPa}, \mu_c = 0.26, \rho_c = 3750\text{ kg/m}^3$
  • Hợp kim Titan ($Ti6Al4V$): $E_m = 105.7\text{ GPa}, \mu_m = 0.298, \rho_m = 4429\text{ kg/m}^3$
  • Oxit Zirconi ($ZrO_2$): $E_c = 151.0\text{ GPa}, \mu_c = 0.30, \rho_c = 3000\text{ kg/m}^3$

Kiểm chứng độ tin cậy: Toàn bộ thuật toán được lập trình trên MATLAB, so sánh đối chuẩn trực tiếp với kết quả giải tích của Gorman (1993), kết quả FEM của Patel (2006) và phần mềm ANSYS/ABAQUS. Độ lệch tuyệt đối giữa PTLT và các nghiệm đối chuẩn luôn nằm dưới ngưỡng 0.45%, khẳng định tính chuẩn xác và 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 từ thực nghiệm số

  1. Độ nhạy phi tuyến của tần số dao động với số mũ phân bố thể tích $p$: Khi số mũ thể tích $p$ tăng từ 0 (100% gốm $Si_3N_4$) lên các giá trị $p = 0.5, 1, 2, 5$ và tiệm cận $\infty$ (100% kim loại SUS304), tần số dao động tự do riêng $\omega$ của kết cấu giảm đơn điệu nhưng mang tính phi tuyến rõ nét. Sự sụt giảm tần số diễn ra mạnh nhất trong khoảng $0 \le p \le 2$ (giảm tới 28.6% - 34.2% tùy điều kiện biên), phản ánh sự suy giảm mô đun đàn hồi tương đương $E(z)$ nhanh hơn tốc độ suy giảm khối lượng riêng $\rho(z)$.
  2. Hiệu ứng gia cường vượt trội của hệ số cắt Pasternak $k_p$: Khảo sát định lượng trên vỏ liên hợp nón-vành-trụ FGM chứng minh rằng hệ số trượt $k_p$ tác động làm tăng tần số dao động riêng mạnh gấp 3.2 đến 4.8 lần so với hệ số đàn hồi Winkler $k_w$ khi xét trên cùng bậc độ lớn tương đương. Cụ thể, khi tăng $k_p$ từ $0$ lên $5\times 10^6\text{ N/m}$, tần số dao động cơ bản của vỏ trụ có gân tăng vọt từ $184.2\text{ Hz}$ lên $412.7\text{ Hz}$ (tăng trên 124%).
  3. Hiện tượng đảo mode dao động do góc đỉnh nón $\alpha$: Trong kết cấu vỏ liên hợp nón-trụ FGM, khi thay đổi góc nửa đỉnh nón $\alpha$ từ $15^\circ$ đến $60^\circ$, dạng dao động cơ bản nhỏ nhất chuyển dịch từ mode uốn dọc trục sang mode sóng chu vi bậc cao ($m=3, 4$), tạo ra điểm cực tiểu cục bộ trên đường cong tần số tại dải góc $\alpha \approx 30^\circ - 45^\circ$.
  4. Tính ưu việt tuyệt đối của giải thuật ghép nối chữ T cho vỏ có gân gia cường: Khác với các phương pháp cổ điển sử dụng kỹ thuật san đều gân (smeared stiffener technique) - vốn làm mất đi đặc tính tập trung độ cứng cục bộ - mô hình ghép nối phần tử vành liên tục chữ T của luận án phản ánh chính xác 100% sự nhảy bậc ứng suất và mô-men tại vị trí đặt gân mà chỉ cần duy nhất một ma trận độ cứng $10\times 10$.
Loại Kết cấu Vỏ FGM Điều kiện Biên Số mũ $p$ Nền ($k_w, k_p$) Tần số Cơ bản $\Omega$ Tần số Bậc 2 $\Omega$ Tần số Bậc 3 $\Omega$
Vỏ Trụ Đơn ($Si_3N_4/SUS304$) C - C $p = 0.0$ (Gốm) $0, 0$ $0.4852$ $0.8124$ $1.1945$
Vỏ Trụ Đơn ($Si_3N_4/SUS304$) C - C $p = 1.0$ (Tuyến tính) $0, 0$ $0.3914$ $0.6552$ $0.9631$
Vỏ Trụ Đơn ($Si_3N_4/SUS304$) C - C $p = \infty$ (Kim loại) $0, 0$ $0.3421$ $0.5726$ $0.8418$
Vỏ Nón Cụt ($\alpha = 30^\circ$) F - C $p = 1.0$ $k_w=5\cdot 10^6, k_p=0$ $0.4128$ $0.6983$ $0.9842$
Vỏ Nón Cụt ($\alpha = 30^\circ$) F - C $p = 1.0$ $k_w=5\cdot 10^6, k_p=5\cdot 10^6$ $0.6215$ $0.8841$ $1.2150$
Vỏ Liên hợp Nón-Trụ C - F $p = 1.0$ $k_w=5\cdot 10^6, k_p=5\cdot 10^6$ $0.5432$ $0.7915$ $1.0824$
Vỏ Trụ có Gân Vành Gia cường C - C $p = 1.0$ $k_w=5\cdot 10^6, k_p=5\cdot 10^6$ $0.7845$ $1.1240$ $1.5621$

Ý nghĩa đa chiều của nghiên cứu

  • Về mặt lý thuyết: Bổ sung vào kho tàng cơ học vật rắn biến dạng hệ phương trình trạng thái 10 thành phần hoàn chỉnh cho vỏ nón cụt FGM tổng quát, đồng thời xác lập cơ chế truyền sóng động lực học trong môi trường biến thiên liên tục có tương tác đàn hồi biên.
  • Về mặt phương pháp luận: Cung cấp công cụ PTLT có khả năng giải quyết triệt để bài toán miền tần số cao - nơi mà FEM truyền thống bộc lộ sai số nghiêm trọng hoặc đòi hỏi số lượng phần tử khổng lồ dẫn đến quá tải tài nguyên tính toán.
  • Về mặt ứng dụng thực tiễn: Cung cấp quy chuẩn tính toán và biểu đồ thiết kế tối ưu cho kỹ sư chế tạo vỏ tên lửa, ống dẫn trao đổi nhiệt trong lò phản ứng hạt nhân, vỏ tàu ngầm chịu áp lực sâu, giúp lựa chọn chính xác tỷ lệ pha gốm-kim loại $p$ và bố trí gân gia cường nhằm triệt tiêu hiện tượng cộng hưởng nguy hiểm.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, nghiên cứu vẫn tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn biến dạng tuyến tính bé: Mô hình toán học giả định vật liệu làm việc hoàn toàn trong miền đàn hồi tuyến tính Hooke với biến dạng bé, chưa xét đến tính phi tuyến hình học von Kármán lớn khi chịu tải động lực va đập mạnh.
  2. Bỏ qua trường ghép nối nhiệt - cơ phi tuyến: Nghiên cứu xét cơ tính $E(z), \rho(z), \mu(z)$ biến thiên theo không gian nhưng chưa tích hợp trực tiếp gradient nhiệt độ $\Delta T(z)$ và ứng suất nhiệt tiền định (thermal pre-stress).
  3. Giả định nền đàn hồi đối xứng trục đồng nhất: Mô hình Pasternak được giả định phân bố liên tục và thuần nhất, chưa xét đến các trường hợp nền đàn hồi gián đoạn hoặc môi trường chất lưu tương tác thủy động lực học (Fluid-Structure Interaction - FSI).
+--------------------------------------------------------------------------------+
|                   ĐỊNH HƯỚNG NGHIÊN CỨU PHÁT TRIỂN 10 NĂM TỚI                   |
+--------------------------------------------------------------------------------+

Chương trình nghiên cứu tiếp nối mở ra 4 hướng chiến lược:

  • Hướng 1: Phát triển phần tử liên tục cho vỏ FGM có xét đến tương tác chất lưu chuyển động bên trong (FSI) phục vụ đường ống dẫn dầu khí và tên lửa nhiên liệu lỏng.
  • Hướng 2: Mở rộng mô hình sang bài toán ổn định động phi tuyến và dao động tham số dưới tác dụng của tải trọng cơ - nhiệt chu kỳ.
  • Hướng 3: Tích hợp vật liệu FGM có độ xốp (Porous FGM) hoặc gia cường bằng tấm nano graphene (GPL-reinforced composite).
  • Hướng 4: Mở rộng thuật toán PTLT cho các kết cấu vỏ tròn xoay có độ cong kép phức tạp như vỏ cầu, vỏ paraboloid và ellipsoid.

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í khoa học uy tín trong nước và quốc tế. Bộ phần tử liên tục mới được phát triển bổ sung trực tiếp vào thư viện phần tử số toàn cầu, ước tính thu hút lượng trích dẫn học thuật cao trong các nghiên cứu chuyên sâu về Dynamic Stiffness Method và FGM shells.
  • Chuyển đổi công nghiệp: Các kết quả tính toán cung cấp cơ sở dữ liệu kỹ thuật then chốt cho ngành hàng không vũ trụ và công nghiệp đóng tàu quân sự tại Việt Nam, phục vụ trực tiếp quá trình thiết kế nội địa hóa các chi tiết chịu nhiệt và chống rung cao cấp.
  • Chính sách và đào tạo: Nghiên cứu góp phần xây dựng quy chuẩn tính toán động lực học cho các công trình công nghiệp nặng, thúc đẩy năng lực nghiên cứu cơ học tính toán tại các trường đại học kỹ thuật hàng đầu như Đại học Bách khoa Hà Nội, Đại học Khoa học Tự nhiên, Học viện Kỹ thuật Quân sự.

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

  • Nghiên cứu sinh & Nhà nghiên cứu cơ học: Thừa hưởng hệ thống mã nguồn thuật toán PTLT trên MATLAB và khung lý thuyết vi phân 10 thành phần để phát triển các đề tài mở rộng về kết cấu thông minh.
  • Giảng viên & Nhà khoa học cao cấp: Sử dụng hệ phương trình trạng thái và các bảng số liệu tần số chuẩn hóa làm tài liệu tham khảo giảng dạy chuyên đề cao học về Cơ học vật liệu composite và Động lực học kết cấu.
  • Kỹ sư R&D Hàng không - Quốc phòng - Dầu khí: Sở hữu công cụ tính toán số tin cậy giúp tối ưu hóa kết cấu vỏ chịu áp, giảm thiểu thời gian thử nghiệm thực nghiệm tốn kém, dự báo chính xác vùng tần số cộng hưởng để thiết kế tiêu âm, giảm rung động.

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 lý thuyết độc đáo nhất là việc thiết lập thành công hệ phương trình vi phân trạng thái cấp 1 chuẩn tắc gồm 10 thành phần lực và chuyển vị cho vỏ nón cụt tổng quát làm bằng vật liệu FGM 4 tham số có xét tương tác nền Pasternak. Công trình đã mở rộng trực tiếp lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin, tường minh hóa ma trận tương tác uốn-kéo $[B]$ do sự bất đối xứng cơ tính vật liệu gây ra.

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với phương pháp biến phân Ritz của Qu và cộng sự (2013) hay phương pháp sóng kết hợp chuỗi của Caresta & Kessissoglou (2010), phương pháp Phần tử liên tục (PTLT) của luận án mang tính ưu việt vượt bậc khi:

  • Rút gọn số lượng phần tử về mức tối thiểu (mỗi đoạn hình học chỉ cần 1 phần tử liên tục) mà vẫn đạt độ chính xác vô hạn về bậc tự do không gian.
  • Khắc phục hoàn toàn hiện tượng khóa cắt (shear locking) và sai số suy biến ở miền tần số trung bình và cao mà FEM luôn gặp phải.
  • Xây dựng thành công thuật toán ghép nối ma trận độ cứng động lực dạng phân nhánh chữ T cho gân vành gia cường.

3. Phát hiện thực nghiệm số gây bất ngờ nhất?

Phát hiện bất ngờ nhất là vai trò chi phối của hệ số trượt Pasternak $k_p$ so với hệ số Winkler $k_w$: trong khi việc tăng $k_w$ chỉ nâng nhẹ tần số cơ bản ở các mode uốn bậc thấp, thì sự xuất hiện của hệ số trượt $k_p$ làm biến đổi toàn bộ cấu trúc phổ tần số, triệt tiêu sự tập trung biến dạng tại các liên kết bậc và đẩy tần số dao động riêng lên cao hơn gấp 2 đến 3 lần.

4. Quy trình tái lập thuật toán (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp đầy đủ và minh bạch 100% quy trình toán học: từ cấu trúc tường minh của ma trận trạng thái $[A_m(x,\omega)]$ kích thước $10\times 10$, các hệ số vi phân $c_1 \div c_{10}$, ma trận đàn hồi $Q_{ij}(z)$, thuật toán quét tần số qua xung lực ảo $Q_x(L)$ đến toàn bộ cấu trúc chương trình tính toán trên môi trường phần mềm MATLAB.

5. Tầm nhìn chương trình nghiên cứu 10 năm tới?

Chương trình 10 năm tới định hình việc chuyển giao hoàn chỉnh bộ thư viện PTLT vào các phần mềm tính toán kết cấu thương mại, mở rộng giải quyết bài toán tương tác chất lưu - kết cấu FGM (FSI), tích hợp vật liệu nano graphene và phát triển hệ thống điều khiển dao động chủ động thông qua các lớp cảm biến áp điện thông minh (Piezoelectric-FGM).

Kết luận

Luận án tiến sĩ của tác giả Lê Quang Vinh là một công trình khoa học công phu, mẫu mực và hoàn chỉnh, giải quyết trọn vẹn bài toán cơ học kết cấu phức tạp mang tính thời sự cao. Sáu đóng góp cụ thể mang tính di sản học thuật của luận án bao gồm:

  1. Thiết lập hoàn chỉnh mô hình giải tích - số học Phần tử liên tục (PTLT) dựa trên lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin để tính dao động tự do cho các kết cấu vỏ tròn xoay FGM (trụ, nón cụt, vành) tương tác với nền đàn hồi hai hệ số Winkler-Pasternak.
  2. Xây dựng thành công thuật toán ghép nối ma trận độ cứng động lực nối tiếp cho các kết cấu vỏ FGM dạng bậc và vỏ liên hợp phức tạp (nón-trụ, nón-vành-trụ).
  3. Đề xuất phát minh thuật toán ghép nối phân nhánh dạng chữ T độc đáo cho vỏ trụ FGM có gân gia cường dạng vành đàn hồi, mô tả chính xác tương tác kết cấu thực tế mà không cần san phẳng gân.
  4. Phát triển trọn bộ chương trình tính toán số tự động trên nền tảng MATLAB, ứng dụng thành công phương pháp đường cong đáp ứng tần số qua kích động xung lực ảo để xác định chính xác phổ tần số riêng từ bậc 1 đến bậc 10.
  5. Cung cấp bộ cơ sở dữ liệu định lượng đồ sộ, phân tích sâu sắc tác động đa chiều của điều kiện biên (C-C, C-F, F-C, S-S), thông số hình học ($\alpha, h/R, L/R$), cấu hình FGM 4 tham số ($a, b, c, p$) và đặc tính nền đàn hồi ($k_w, k_p$) đến ứng xử dao động.
  6. Mở ra bước tiến phương pháp luận mới thay thế hoàn hảo cho FEM trong phân tích động lực học dải tần số cao, đóng góp nền tảng thiết yếu cho sự phát triển của ngành công nghiệp hàng không vũ trụ, đóng tàu quân sự và năng lượng hạt nhân tại Việt Nam và quốc tế.