Tổng quan về luận án

Nghiên cứu dao động của các kết cấu vỏ tròn xoay phức hợp hoạt động trong môi trường đàn hồi là một trong những bài toán kinh điển nhưng luôn mang tính thời sự cấp thiết của cơ học vật rắn biến dạng và kỹ thuật công trình hiện đại. Trong các kết cấu tiên tiến thuộc ngành hàng không vũ trụ (vỏ tên lửa, thân máy bay siêu âm), công trình biển (thân tàu ngầm, ống dẫn ngầm đáy biển), và kỹ thuật năng lượng hạt nhân (vỏ lò phản ứng, bình áp lực), vật liệu composite truyền thống dạng lớp thường bộc lộ hạn chế nghiêm trọng về hiện tượng tập trung ứng suất tiếp và nguy cơ bong tách lớp (delamination) do sự gián đoạn cơ tính tại các mặt tiếp giáp. Sự ra đời của vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM), khởi xướng từ năm 1984 tại Nhật Bản, đã mang lại bước đột phá nhờ khả năng phân bố trơn và liên tục các pha vật liệu (như gốm $Si_3N_4, ZrO_2, Al_2O_3$ kết hợp với kim loại $SUS304, Ti6Al4V, Al, Ni$) theo chiều dày thành vỏ, loại bỏ hoàn toàn bề mặt phân cách đột ngột, tối ưu hóa khả năng chịu nhiệt và phân tán ứng suất động.

Tuy nhiên, khoảng trống nghiên cứu (research gap) lớn trong cơ học kết cấu FGM hiện nay là: phần lớn các công trình công bố quốc tế tập trung vào các kết cấu vỏ đơn lẻ, đồng nhất về hình học (như nghiên cứu của Sofiyev [2009-2012] cho vỏ nón, Sheng & Wang [2008-2010] cho vỏ trụ) hoặc sử dụng phương pháp phần tử hữu hạn (PTHH) kinh điển vốn đòi hỏi chia lưới cực mịn, phát sinh chi phí tính toán khổng lồ và dễ gặp hiện tượng khóa cắt (shear locking) ở miền tần số trung bình và cao. Ngược lại, các kết cấu vỏ liên hợp phức tạp trong thực tế như vỏ nón - trụ (conical-cylindrical joined shells), nón - vành - trụ (conical-annular-cylindrical joined shells), vỏ có chiều dày dạng bậc và đặc biệt là vỏ trụ FGM gia cường bằng các gân vành đàn hồi tương tác với nền đàn hồi hai thông số Winkler - Pasternak hầu như chưa có một công cụ giải tích - số chính xác cao nào dựa trên phương pháp Phần tử liên tục (Continuous Element Method - CEM / Ma trận độ cứng động lực Dynamic Stiffness Method - DSM).

Luận án tiến sĩ cơ học của nghiên cứu sinh Lê Quang Vinh, 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 trọn vẹn khoảng trống học thuật này. Các câu hỏi và giả thuyết nghiên cứu được thiết lập cụ thể:

  • Câu hỏi nghiên cứu 1 (Q1): Làm thế nào để xây dựng ma trận độ cứng động lực chính xác cho từng phần tử vỏ tròn xoay FGM đơn (trụ, nón cụt, vành tròn) dựa trên lý thuyết biến dạng cắt bậc nhất của Reissner - Mindlin có tương tác với nền đàn hồi hai thông số?
  • Câu hỏi nghiên cứu 2 (Q2): Thuật toán ghép nối ma trận độ cứng động lực dạng nối tiếp và dạng phân nhánh chữ T (T-junction assembly) cần được thiết lập như thế nào để mô hình hóa chính xác các kết cấu liên hợp dạng bậc, ghép nối đa phân đoạn và vỏ có gân gia cường đàn hồi mà không phải dùng các giả thiết san đều gân kinh điển?
  • Câu hỏi nghiên cứu 3 (Q3): Quy luật biến thiên 4 tham số $FGM(a/b/c/p)$, góc nửa đỉnh nón $\alpha$, và các hệ số nền Winkler ($k_w$), Pasternak ($k_p$) tương tác phi tuyến như thế nào đến phổ tần số dao động tự do?
  • Giả thuyết nghiên cứu 1 (H1): Phương pháp Phần tử liên tục với số lượng phần tử tối thiểu (bằng đúng số đoạn kết cấu thực tế) có thể triệt tiêu sai số rời rạc hóa hình học của PTHH và cho độ chính xác tiệm cận tuyệt đối ở cả dải tần số thấp, trung bình và cao.
  • Giả thuyết nghiên cứu 2 (H2): Mô hình hóa gân gia cường dưới dạng một phần tử vành tròn đàn hồi liên tục ghép nối song song dạng chữ T sẽ phản ánh chính xác tương tác lực - mô men cục bộ hơn hẳn phương pháp san gân (smeared stiffener technique).

Khung lý thuyết của công trình tích hợp lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin (FSDT), mô hình nền đàn hồi hai thông số Winkler - Pasternak ($k_w, k_p$), luật phân bố tỉ phần thể tích 4 tham số $FGM(a/b/c/p)$, và phương pháp ma trận độ cứng động lực kết hợp phương pháp đường cong đáp ứng (frequency response curve method). Phạm vi khảo sát trải rộng trên nhiều hệ vật liệu công nghiệp ($Si_3N_4/SUS304, ZrO_2/Ti6Al4V, Al_2O_3/Al, Si_3N_4/Ni$), các điều kiện biên đa dạng (Ngàm - Ngàm C-C, Tựa - Tựa S-S, Tự do - Ngàm F-C, Ngàm - Tự do C-F), với thông số nền đàn hồi $k_w \in [0, 5\times 10^7]\text{ N/m}^3$ và $k_p \in [0, 5\times 10^8]\text{ N/m}$.

Literature Review và Positioning

Lịch sử nghiên cứu động lực học kết cấu vỏ tròn xoay ghi nhận sự phân nhánh rõ rệt giữa các trường phái tiếp cận lý thuyết và số trị:

Luồng nghiên cứu thứ nhất tập trung vào phương pháp giải tích và bán giải tích trên kết cấu FGM đơn giản. Shen [2009] và Bagherizadeh cùng cộng sự [2011] đã sử dụng lý thuyết biến dạng cắt bậc cao (HSDT) và phương pháp tham số bé để khảo sát mất ổn định và dao động vỏ trụ FGM. Sofiyev [2009-2012] cùng Kuruoglu [2011] ứng dụng lý thuyết màng Donnell và phương pháp Galerkin để giải dao động phi tuyến của vỏ nón cụt FGM trên nền đàn hồi. Điểm hạn chế cốt tử của luồng nghiên cứu này là chỉ áp dụng được cho các điều kiện biên đơn giản (thường là tựa đơn S-S) và bất khả thi khi mở rộng cho các hình học vỏ liên hợp phức tạp hoặc có bước nhảy tiết diện.

Luồng nghiên cứu thứ hai ứng dụng phương pháp phần tử hữu hạn (PTHH) và các phương pháp biến phân cải biên. Qu cùng cộng sự [2013, 2014] đã công bố phương pháp biến phân tổng quát cho hệ vỏ liên hợp nón - trụ - vành và nón - trụ - cầu FGM có gân gia cường. Ma cùng cộng sự [2014] áp dụng phương pháp Fourier - Ritz mở rộng cho vỏ nón - trụ ghép nối. Tuy nhiên, PTHH truyền thống (như triển khai trên ANSYS, ABAQUS) luôn đối mặt với hiện tượng số bậc tự do bùng nổ khi nâng cao dải tần số khảo sát, dẫn đến sai số tích lũy do hàm dạng tĩnh không mô tả được sóng dao động cục bộ có bước sóng ngắn.

Luồng nghiên cứu thứ ba là sự phát triển của Ma trận độ cứng động lực (DSM) và Phần tử liên tục (CEM). Khởi nguồn từ các nghiên cứu kinh điển về dầm của Williams & Kennedy [1987], Banerjee [1989-1997], phương pháp được mở rộng cho tấm dày kim loại bởi Gorman [1993, 2000] và Casimir [2003] (nền tảng của phần mềm ETAPE thuộc Hải quân Pháp), sau đó là tấm composite bởi Boscolo [2012]. Đối với vỏ tròn xoay kim loại và composite, Nguyen Manh Cuong & Casimir [2009], Tran Ich Thinh & Nguyen Manh Cuong [2013, 2017], Vu Quoc Hien [2016], và Le Thi Bich Nam [2018] đã thiết lập thành công ma trận độ cứng động lực cho vỏ trụ composite chứa chất lỏng và vỏ liên hợp composite trên nền đàn hồi. Tại Việt Nam, trường phái nghiên cứu của Đào Huy Bích [2009-2016] và Nguyễn Đình Đức [2011-2017] đạt nhiều thành tựu vang dội trong việc phân tích phi tuyến tĩnh và động của vỏ FGM có gân gia cường nhưng chủ yếu dựa trên kỹ thuật san đều gân và lý thuyết Donnell giải tích.

Về mặt định vị học thuật, luận án của Lê Quang Vinh là công trình tiên phong trên thế giới tích hợp thành công phương pháp Phần tử liên tục cho các kết cấu vỏ liên hợp FGM đa phân đoạn (nón - trụ, nón - vành - trụ, trụ bậc, nón bậc) và vỏ trụ FGM có gân gia cường đàn hồi tương tác với nền đàn hồi Pasternak. Khi so sánh với nghiên cứu quốc tế của Kouchakzadeh & Shakouri [2014] (nghiên cứu vỏ composite nón - nón ghép nối dùng lý thuyết Donnell) và nghiên cứu của Caresta & Kessissoglou [2010] (phân tích vỏ nón - trụ đẳng hướng dùng chuỗi sóng bậc cao), mô hình CEM của luận án thể hiện tính vượt trội về khả năng mô phỏng trực tiếp trường chuyển vị động học mà không cần xấp xỉ chuỗi hình học phức tạp, đồng thời giữ nguyên vẹn liên kết liên tục lực - mô men tại các đường biên ghép nối.

                  ┌────────────────────────────────────────────────────────┐
                  │                 TIẾP CẬN KINH ĐIỂN                     │
                  │  - Giải tích (Galerkin, Ritz, Donnell, Flügge)         │
                  │  - Phần tử hữu hạn (ANSYS, ABAQUS - Hàm dạng tĩnh)     │
                  │  - Kỹ thuật san đều gân (Smeared stiffener)            │
                  └──────────────────────────┬─────────────────────────────┘
                                             │  Hạn chế: Bùng nổ bậc tự do,
                                             │  khóa cắt ở tần số cao,
                                             │  mất liên tục ứng suất cục bộ
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │         ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN (CEM/DSM)         │
                  │  - Hàm dạng là nghiệm động lực học chính xác           │
                  │  - Số lượng phần tử tối thiểu (1 phân đoạn = 1 PTLT)   │
                  │  - Ghép nối chữ T mô hình hóa gân vành đàn hồi rời rạc │
                  │  - Nền Pasternak 2 thông số tích hợp trực tiếp         │
                  └────────────────────────────────────────────────────────┘

Đó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 biên giới của Lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin khi áp dụng cho môi trường vật chất dị hướng biến thiên liên tục theo không gian 3 chiều của vỏ tròn xoay. Công trình đã chính xác hóa hệ phương trình vi phân dao động của phần tử vỏ FGM khi chịu đồng thời lực phục hồi đàn hồi Winkler ($k_w$) và lực cắt tương tác mặt nền Pasternak ($k_p$).

Mô hình lý thuyết thiết lập luật phân bố tỉ phần thể tích gốm $V_c(z)$ theo quy luật lũy thừa tổng quát 4 tham số: $$V_c(z) = \left[ 1 - a\left( \frac{1}{2} + \frac{z}{h} \right) + b\left( \frac{1}{2} + \frac{z}{h} \right)^c \right]^p \quad (\text{kiểu } FGMI)$$ $$V_c(z) = \left[ 1 - a\left( \frac{1}{2} - \frac{z}{h} \right) + b\left( \frac{1}{2} - \frac{z}{h} \right)^c \right]^p \quad (\text{kiểu } FGMII)$$

Từ đó, các đặc trưng cơ học hiệu dụng gồm mô đun đàn hồi $E(z)$, khối lượng riêng $\rho(z)$, và hệ số Poisson $\mu(z)$ được xác định liên tục theo chiều dày: $$E(z) = (E_c - E_m)V_c(z) + E_m, \quad \rho(z) = (\rho_c - \rho_m)V_c(z) + \rho_m, \quad \mu(z) = (\mu_c - \mu_m)V_c(z) + \mu_m$$

Hệ thống quan hệ ứng suất - biến dạng tích hợp đầy đủ ma trận độ cứng màng $[A]$, ma trận tương tác uốn - màng $[B]$, ma trận độ cứng uốn $[D]$, và ma trận độ cứng cắt trượt $[F]$: $$\begin{pmatrix} {N} \ {M} \end{pmatrix} = \begin{bmatrix} [A] & [B] \ [B] & [D] \end{bmatrix} \begin{pmatrix} {\varepsilon^0} \ {\kappa} \end{pmatrix}, \quad {Q} = [F]{\gamma}$$

Bốn mệnh đề lý thuyết cốt lõi được chứng minh:

  • Mệnh đề 1 (P1): Sự xuất hiện của ma trận ghép $[B] \neq 0$ trong vật liệu FGM bất đối xứng tạo ra tương tác uốn - kéo nén nội tại, làm biến đổi hoàn toàn phổ mode dao động so với vỏ kim loại đồng chất.
  • Mệnh đề 2 (P2): Ma trận độ cứng động lực $[K(\omega)]$ của phần tử liên tục là ma trận đối xứng, thực, và chứa nghiệm hàm sóng chính xác, cho phép biểu diễn miền liên tục vô hạn bậc tự do bằng một quan hệ phần tử biên hữu hạn.
  • Mệnh đề 3 (P3): Thuật toán ghép nối song song chữ T bảo toàn tuyệt đối tính liên tục về chuyển vị ($u, v, w, \beta_x, \beta_\theta$) và cân bằng lực tại nút ghép 3 nhánh giữa vỏ trụ và gân vành gia cường.
  • Mệnh đề 4 (P4): Ảnh hưởng làm cứng kết cấu của nền Pasternak thông qua hệ số trượt $k_p$ có tính chất phi tuyến bậc cao theo góc dốc đường sinh nón $\alpha$.
   ┌────────────────────────────────────────────────────────────────────────┐
   │               KHUNG PHÂN TÍCH MA TRẬN ĐỘ CỨNG ĐỘNG LỰC                │
   └───────────────────────────────────┬────────────────────────────────────┘
                                       │
        ┌──────────────────────────────┴──────────────────────────────┐
        ▼                                                             ▼
┌──────────────────────────────┐                            ┌──────────────────────────────┐
│  MÔ HÌNH HỌC & VẬT LIỆU      │                            │  CƠ HỌC MÔI TRƯỜNG LIÊN TỤC  │
│  - Vỏ: Nón, Trụ, Vành        │                            │  - Reissner-Mindlin FSDT     │
│  - 4 tham số FGM(a/b/c/p)    │                            │  - Nền Winkler-Pasternak     │
│  - Ma trận [A], [B], [D], [F]│                            │  - Cân bằng động lực học     │
└──────────────┬───────────────┘                            └──────────────┬───────────────┘
               │                                                           │
               └──────────────────────────────┬────────────────────────────┘
                                              ▼
               ┌───────────────────────────────────────────────────────────┐
               │    THIẾT LẬP VÉC TƠ TRẠNG THÁI 10 THÀNH PHẦN {y(x, ω)}    │
               │         dy/dx = [A(x, ω)]y ──► {y(x)} = [T(x, ω)]{y(0)}   │
               └──────────────────────────────┬────────────────────────────┘
                                              ▼
               ┌───────────────────────────────────────────────────────────┐
               │         MA TRẬN ĐỘ CỨNG ĐỘNG LỰC PHẦN TỬ [K(ω)]           │
               │             {F(ω)} = [K(ω)] {q(ω)}                        │
               └──────────────────────────────┬────────────────────────────┘
                                              ▼
               ┌───────────────────────────────────────────────────────────┐
               │      GHÉP NỐI TOÀN CỤC & ĐƯỜNG CONG ĐÁP ỨNG TẦN SỐ        │
               │   - Ghép nối tiếp (Stepped/Joined) & Chữ T (Ring-Stiffened)│
               │   - Tìm điểm cực trị đáp ứng ──► Xác định chính xác f(Hz) │
               └───────────────────────────────────────────────────────────┘

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

Điểm sáng tạo đặc biệt của khung phân tích nằm ở việc chuyển đổi toàn bộ bài toán biên giá trị động lực học phức tạp thành hệ phương trình vi phân thường bậc nhất thuần nhất trong không gian trạng thái. Véc tơ trạng thái 10 thành phần được lựa chọn chặt chẽ: $${y(x, \theta, z, \omega, t)} = {u, v, w, \beta_x, \beta_\theta, N_x, N_{x\theta}, Q_x, M_x, M_{x\theta}}^T$$

Phương trình vi phân trạng thái có dạng: $$\frac{d{y(x, \omega)}}{dx} = [A(x, \omega)]{y(x, \omega)}$$

Giải hệ phương trình này bằng tích phân số phân đoạn ma trận vi phân truyền sóng, luận án thiết lập ma trận truyền $[T(x, \omega)]$ thỏa mãn: $${y(x, \omega)} = [T(x, \omega)]{y(x_0, \omega)}$$ Từ ma trận truyền $[T(x, \omega)]$, ma trận độ cứng động lực phần tử $[K(\omega)]$ được suy dẫn chính xác, liên kết trực tiếp véc tơ lực nút ${F(\omega)}$ và véc tơ chuyển vị nút ${q(\omega)}$.

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

Thiết kế nghiên cứu

Luận án vận dụng thế giới quan duy thực tất định (positivism/mechanistic determinism) của cơ học môi trường liên tục giải tích kết hợp với thuật toán số trị chính xác cao trên máy tính. Bản chất của thiết kế nghiên cứu là thiết kế số cấp bậc đa tầng (hierarchical multi-level computational design):

  • Cấp độ 1: Phần tử liên tục cơ bản vỏ đơn (nón cụt, trụ, vành tròn phẳng).
  • Cấp độ 2: Kết cấu vỏ liên hợp dạng bậc (nón bậc 4 phân đoạn, trụ bậc 2-4 phân đoạn, vành bậc 2-3 phân đoạn).
  • Cấp độ 3: Kết cấu vỏ liên hợp ghép nối không đồng dạng hình học (nón - trụ, nón - vành - trụ).
  • Cấp độ 4: Kết cấu vỏ trụ FGM gia cường bằng các gân vành đàn hồi rời rạc đặt trong nền đàn hồi Pasternak.

Quy trình nghiên cứu rigorous

Quy trình tính toán số được lập trình toàn diện trong môi trường MATLAB thông qua các bộ mã nguồn chuyên dụng do tác giả tự phát triển: ThesisFGMJoinedshellNT, ThesisFGMJoinedshellNTNF, và ThesisFGMJoinedshellTG.

Quy trình giải bài toán tần số riêng tự do bao gồm các bước nghiêm ngặt:

  1. Xác định hàm phân bố vật liệu $V_c(z)$ và tính tích phân theo chiều dày để thu được các ma trận độ cứng vật liệu $[A], [B], [D], [F]$.
  2. Xây dựng ma trận động học trạng thái $[A(x, \omega)]$ dọc theo biến tọa độ đường sinh $x$.
  3. Tính toán ma trận truyền $[T(\omega)]$ qua phương pháp Peano - Baker hoặc tích phân Runge - Kutta bậc 4.
  4. Chuyển đổi ma trận truyền thành ma trận độ cứng động lực $[K(\omega)]^{(e)}$ cho từng phần tử.
  5. Ghép nối ma trận độ cứng động lực toàn cục $[K(\omega)]_{\text{global}}$ theo điều kiện liên tục chuyển vị và cân bằng nội lực tại các mặt cắt tiếp giáp.
  6. Xác định tần số dao động tự do thông qua Phương pháp đường cong đáp ứng (Frequency Response Curve Method): Đặt một lực kích thích điều hòa đơn vị tại vị trí biên tự do hoặc điểm tiếp giáp, quét dải tần số $\omega$ để xác định các đỉnh cộng hưởng đáp ứng chuyển vị biên độ cực đại. Phương pháp này vượt qua triệt để các khó khăn về hàm siêu việt phân kỳ mà thuật toán Williams - Wittrick kinh điển hay gặp phải khi giải phương trình định thức $\det |K(\omega)| = 0$.

Độ tin cậy và tính hợp thức (validity & reliability) của phương pháp được kiểm chứng chéo (triangulation) thông qua việc so sánh kết quả với:

  • Nghiên cứu giải tích của Kouchakzadeh & Shakouri [2014] cho vỏ ghép nối nón - trụ kim loại.
  • Mô hình PTHH 3D chạy trên phần mềm thương mại ANSYS sử dụng phần tử vỏ SHELL281 với lưới chia cực mịn.
  • Các kết quả công bố của Sofiyev, Gorman, Qu, và Trần Ích Thịnh. Sai số tương đối luôn duy trì ở mức dưới $0.5% - 1.2%$, khẳng định tính chuẩn xác tuyệt đối của phương pháp đề xuất.
       ┌───────────────────────────────────────────────────────────────┐
       │   ĐIỀU KIỆN BIÊN LIÊN KẾT NGOÀI:                              │
       │   C-C (Ngàm-Ngàm) | F-C (Tự do-Ngàm) | S-S (Tựa-Tựa)          │
       └───────────────────────────────┬───────────────────────────────┘
                                       │
                                       ▼
       ┌───────────────────────────────────────────────────────────────┐
       │           ĐIỀU KIỆN GHÉP NỐI TẠI MẶT TIẾP GIÁP                │
       │                                                               │
       │  1. Ghép nối tiếp (Vỏ nón - Vỏ trụ - Vành bậc):               │
       │     {q_phải} = [T_chuyển trục] {q_trái}                       │
       │     {F_phải} + [T_chuyển trục] {F_trái} = {0}                 │
       │                                                               │
       │  2. Ghép phân nhánh chữ T (Vỏ trụ - Gân vành đàn hồi):        │
       │     {q_vỏ 1} = {q_vỏ 2} = {q_gân}                             │
       │     {F_vỏ 1} + {F_vỏ 2} + {F_gân} = {0}                       │
       └───────────────────────────────┬───────────────────────────────┘
                                       │
                                       ▼
       ┌───────────────────────────────────────────────────────────────┐
       │             HỆ PHƯƠNG TRÌNH ĐỘ CỨNG ĐỘNG LỰC                  │
       │          [K(ω)_toàn cục] {q_toàn cục} = {F_kích thích}        │
       └───────────────────────────────────────────────────────────────┘

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

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

Dữ liệu bằng chứng trích xuất trực tiếp từ luận án đối với kết cấu vỏ trụ FGM có gân gia cường ($R_1 = 1\text{ m}; L_1 = L_2 = 2\text{ m}; h = 0.1\text{ m}; R_r = 0.1\text{ m}; h_r = 0.1\text{ m}$; vật liệu $Si_3N_4/Ni$ với $E_c = 322.27\text{ GPa}, \rho_c = 2370\text{ kg/m}^3, \mu_c = 0.24$ và $E_m = 205.98\text{ GPa}, \rho_m = 8900\text{ kg/m}^3, \mu_m = 0.31$; điều kiện biên F-C; hàm FGMI $a=1, b=0.5, c=2, p$):

| Bậc dao động (Mode) | $p = 0.5$ (Hz) | $p = 2$ (Hz) | $p = 5$ (Hz) | $p = 20$ (Hz) | $p = 100$ (Hz) | Mức suy giảm tần số khi tăng $p$ | | :--- | :--- | :--- | :--- | :--- | :--- | | Mode 1 | 158.2 | 116.7 | 97.3 | 86.1 | 83.3 | -47.35% | | Mode 2 | 172.0 | 126.0 | 104.0 | 93.0 | 91.0 | -47.09% | | Mode 3 | 345.0 | 254.0 | 213.0 | 190.0 | 182.0 | -47.25% | | Mode 4 | 358.0 | 262.0 | 218.0 | 195.0 | 189.0 | -47.21% | | Mode 5 | 458.0 | 337.0 | 281.0 | 251.0 | 242.0 | -47.16% | | Mode 10| 687.0 | 502.0 | 440.0 | 393.0 | 377.0 | -45.12% | | Mode 12| 831.0 | 611.0 | 548.0 | 491.0 | 475.0 | -42.84% |

Bốn phát hiện khoa học đột phá:

  1. Quy luật chi phối phi tuyến của chỉ số thể tích $p$: Khi số mũ $p$ tăng từ 0.5 lên 100, tần số tự nhiên của kết cấu giảm đơn điệu và phi tuyến gần $48%$. Nguyên nhân vật lý xuất phát từ việc gia tăng $p$ khiến cấu trúc chuyển dịch từ trạng thái giàu gốm (độ cứng riêng $E/\rho$ cực cao của $Si_3N_4$) sang trạng thái giàu kim loại (mô đun đàn hồi thấp hơn nhưng khối lượng riêng của Nickel tăng gấp 3.75 lần, từ $2370\text{ kg/m}^3$ lên $8900\text{ kg/m}^3$).
  2. Hiện tượng ngưỡng nhạy cảm bậc thang của nền Pasternak: Luận án phát hiện quy luật phản ứng hai pha đối với các hệ số nền:
    • Khi $0 < k_w < 1\times 10^8\text{ N/m}^3$, sự thay đổi độ cứng Winkler ảnh hưởng không đáng kể đến tần số riêng; chỉ khi $k_w \ge 1\times 10^8\text{ N/m}^3$ (và đặc biệt rõ ở $k_w \ge 2.5\times 10^7\text{ N/m}^3$ đối với tương tác kép), tần số dao động mới tăng vọt.
    • Đối với hệ số trượt Pasternak $k_p$, khi $0 < k_p < 2.5\times 10^7\text{ N/m}$, tác động của nền là thứ yếu. Tuy nhiên, khi $k_p \ge 2.5\times 10^7\text{ N/m}$ (khảo sát ở các mức $5\times 10^6, 2.5\times 10^7, 1\times 10^8\text{ N/m}$), các đường cong đáp ứng tách biệt rõ rệt và tần số dao động tăng trưởng với độ dốc rất lớn.
  3. Phản ứng hình học phân kỳ giữa các hệ vỏ liên hợp: Khi tăng góc nửa đỉnh nón $\alpha$, tần số tự nhiên của kết cấu vỏ liên hợp trụ - nón có xu hướng tăng dần do hiệu ứng vòm tăng cứng hình học tại phần đáy nón, nhưng đối với vỏ liên hợp trụ - vành - nón, tần số lại có xu hướng giảm dần do sự suy giảm độ cứng chống uốn cục bộ tại vị trí liên kết với vành chuyển tiếp.
  4. Tính ưu việt của mô hình gân vành đàn hồi chữ T: Khác với phương pháp san gân truyền thống vốn làm mịn hóa độ cứng và bỏ qua tương tác ứng suất cắt cục bộ, mô hình ghép nối chữ T cho phép gân vành tự do biến dạng uốn - xoắn độc lập, phản ánh chính xác các mode dao động cục bộ giữa các gân (inter-ring localized modes) ở dải tần số cao.

Implications đa chiều

  • Ý nghĩa học thuật: Bổ sung vào kho tàng cơ học tính toán thế giới thư viện các Phần tử liên tục FGM hoàn chỉnh cho vỏ nón, trụ, vành, vỏ dạng bậc và vỏ liên hợp phức tạp.
  • Đổi mới phương pháp luận: Thiết lập quy trình tính toán động lực học hiệu năng cao, giảm số lượng bậc tự do từ hàng trăm nghìn (trong PTHH) xuống một số lượng tối thiểu các phần tử liên tục mà không làm suy giảm độ chính xác.
  • Ứng dụng công nghiệp & Quốc phòng: Cung cấp công cụ tính toán chính xác để thiết kế kết cấu vỏ tàu ngầm lặn sâu chịu áp lực ngoài bọc lớp đệm đàn hồi, vỏ đầu đạn tên lửa siêu vượt âm chống cộng hưởng rung động nhiệt - cơ, và vỏ bình phản ứng hạt nhân thế hệ mới.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  • Giả thiết biến dạng đàn hồi tuyến tính và chuyển vị bé, chưa xét đến các hiệu ứng phi tuyến hình học lớn (geometric non-linearity theo nghĩa von Kármán) trong quá trình dao động biên độ lớn.
  • Bỏ qua biến dạng nở ngang theo chiều dày ($\varepsilon_z = 0$) do áp dụng lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin, điều này có thể gây ra sai số nhỏ đối với các kết cấu vỏ cực dày.
  • Nền đàn hồi Winkler - Pasternak được giả định là môi trường đàn hồi tuyến tính liên tục, chưa xét đến yếu tố khối lượng tham gia dao động của nền hoặc lực cản nhớt (damping).
  • Chưa xét đến sự suy giảm cơ tính vật liệu phụ thuộc trực tiếp vào trường nhiệt độ môi trường (temperature-dependent properties).

Chương trình nghiên cứu tiếp nối 5 hướng cụ thể:

  1. Phát triển thuật toán Phần tử liên tục phi tuyến cho kết cấu vỏ FGM chịu kích động tham số lớn và tải trọng cơ - nhiệt kết hợp.
  2. Xây dựng mô hình PTLT tương tác thủy động học (Fluid - Structure Interaction - FSI) cho vỏ liên hợp FGM chứa chất lỏng hoặc chuyển động trong chất lỏng có gân gia cường.
  3. Mở rộng lý thuyết phần tử liên tục bậc cao (HSDT / Zig-zag theory) để mô phỏng chính xác trường ứng suất xuyên lớp cho vỏ FGM composite nhiều lớp lai (hybrid FGM).
  4. Tích hợp mô hình nền đàn hồi ba thông số phi tuyến (Kerr foundation) có xét đến quán tính và độ cản của đất/nền.
  5. Nghiên cứu tối ưu hóa đa mục tiêu quy luật phân bố 4 tham số $FGM(a/b/c/p)$ để chủ động điều khiển và dập tắt dao động cưỡng bức.

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

Công trình tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Tác động học thuật: Mở ra hướng tiếp cận mới trong cơ học tính toán kết cấu FGM tại Việt Nam, đóng góp 14 công trình khoa học chất lượng cao công bố tại các tạp chí chuyên ngành uy tín (Vietnam Journal of Mechanics, Journal of Science and Technology) và các kỷ yếu hội nghị cơ học toàn quốc/quốc tế. Ước tính các kết quả số của luận án sẽ trở thành dữ liệu chuẩn (benchmark data) thu hút hàng trăm trích dẫn từ cộng đồng nghiên cứu cơ học kết cấu FGM toàn cầu.
  • Chuyển đổi công nghệ công nghiệp: Cung cấp thuật toán tối ưu hóa thiết kế kết cấu thành mỏng trong ngành đóng tàu ngầm composite/FGM, hàng không vũ trụ và chế tạo bồn bể chịu áp lực cao.
  • Nâng cao năng lực kỹ thuật quốc gia: Đặt nền móng tự chủ về công nghệ tính toán, mô phỏng động lực học kết cấu tiên tiến mà không phụ thuộc hoàn toàn vào các phần mềm thương mại đắt đỏ của nước ngoài.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận mã nguồn thuật toán CEM/DSM, phương pháp ghép nối ma trận độ cứng động lực dạng chữ T và phương pháp đường cong đáp ứng để giải quyết các bài toán biên phức tạp.
  • Các nhà khoa học đầu ngành: Kế thừa bộ dữ liệu benchmark chính xác tuyệt đối cho các hệ vỏ nón - trụ - vành FGM đặt trên nền đàn hồi để kiểm chứng các phương pháp số mới (IGA, XFEM, Meshless).
  • Kỹ sư R&D hàng không, vũ trụ, hải quân: Sở hữu công cụ tính toán tần số dao động riêng nhanh chóng, chính xác để thiết kế kết cấu chống cộng hưởng, giảm chấn động và tối ưu hóa trọng lượng.
  • Các nhà hoạch định chính sách khoa học: Có cơ sở khoa học để định hướng đầu tư phát triển vật liệu mới FGM trong chiến lược phát triển công nghiệp quốc phòng và công nghiệp nặng quốc gia.

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 độc đáo nhất là việc mở rộng thành công Lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin kết hợp với Phương pháp Phần tử liên tục (CEM/DSM) cho hệ kết cấu vỏ liên hợp phức hợp FGM 4 tham số có gân gia cường đàn hồi tương tác với nền Pasternak hai thông số. Công trình đã chuyển hóa toàn diện phương trình vi phân dao động của môi trường vật chất biến thiên liên tục thành ma trận độ cứng động lực đối xứng chính xác tuyệt đối.

2. Đổi mới phương pháp luận nổi bật khi so sánh với ít nhất hai nghiên cứu quốc tế tiền nhiệm? So với nghiên cứu của Qu cùng cộng sự [2013] (sử dụng phương pháp biến phân với hàm dạng đa thức) và Kouchakzadeh & Shakouri [2014] (sử dụng lý thuyết Donnell giải tích cho vỏ nón - nón ghép nối), phương pháp CEM của luận án:

  • Không bị giới hạn bởi điều kiện biên, giải quyết hoàn hảo mọi tổ hợp biên (C-C, F-C, S-S, C-F).
  • Không yêu cầu chia nhỏ lưới hình học, triệt tiêu sai số làm mềm/làm cứng nhân tạo của hàm dạng tĩnh trong PTHH.
  • Mô hình hóa gân gia cường dạng vành đàn hồi bằng thuật toán ghép nối chữ T chuẩn xác thay vì xấp xỉ san đều gân.

3. Phát hiện thực nghiệm/số trị nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ? Phát hiện bất ngờ nhất là hiện tượng đảo chiều đáp ứng tần số khi thay đổi góc nửa đỉnh nón $\alpha$: Trong khi góc $\alpha$ tăng làm tăng tần số riêng của vỏ liên hợp nón - trụ (do tăng cứng vòm tại đáy nón), thì lại làm giảm tần số riêng của vỏ liên hợp trụ - vành - nón. Đồng thời, sự tồn tại của ngưỡng $k_p \ge 2.5\times 10^7\text{ N/m}$ của nền Pasternak tạo ra bước nhảy vọt về tần số dao động riêng của toàn kết cấu.

4. Quy trình lặp lại nghiên cứu (Replication protocol) có được cung cấp đầy đủ không? Quy trình được cung cấp trọn vẹn với các bước thuật toán chi tiết từ thiết lập véc tơ trạng thái 10 thành phần, ma trận liên hệ vi phân $[A(x,\omega)]$, tích phân ma trận truyền $[T(\omega)]$, ghép nối ma trận độ cứng động lực toàn cục và xác định tần số qua quét phổ đường cong đáp ứng trên ba bộ chương trình MATLAB độc lập (ThesisFGMJoinedshellNT, ThesisFGMJoinedshellNTNF, ThesisFGMJoinedshellTG).

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào? Lộ trình hướng tới: (1) Mô hình hóa FGM phi tuyến hình học biên độ lớn; (2) Tương tác đa trường vật lý FSI (nhiệt - thủy - khí động học); (3) Tối ưu hóa vi cấu trúc FGM in 3D; (4) Điều khiển chủ động dao động vỏ FGM tích hợp các lớp áp điện/từ giảo thông minh.

Kết luận

  1. Xây dựng thành công phương pháp và thuật toán Phần tử liên tục (CEM/DSM) mô phỏng chính xác dao động tự do của các kết cấu vỏ tròn xoay FGM đơn và liên hợp phức tạp trên nền đàn hồi Winkler - Pasternak.
  2. Đề xuất thuật toán ghép nối ma trận độ cứng động lực nối tiếp và song song dạng chữ T độc đáo, cho phép mô hình hóa gân gia cường dạng vành đàn hồi một cách chân thực nhất mà không cần san gân.
  3. Làm sáng tỏ định lượng toàn diện ảnh hưởng của quy luật phân bố 4 tham số $FGM(a/b/c/p)$, điều kiện biên, thông số hình học và hệ số nền ($k_w, k_p$) đến phổ tần số dao động tự do của hệ vỏ.
  4. Chứng minh tính ưu việt vượt trội của CEM so với PTHH truyền thống về tốc độ tính toán, dung lượng bộ nhớ và độ chính xác ở mọi dải tần số khảo sát.
  5. Mở ra 3 nhánh nghiên cứu chuyên sâu mới: Động lực học vỏ liên hợp FGM phi tuyến, tương tác thủy động vỏ FGM chứa chất lỏng, và tối ưu hóa hình học kết cấu FGM đa vùng.
  6. Thiết lập hệ thống dữ liệu chuẩn mực có giá trị khoa học và thực tiễn cao, đóng góp nền tảng vững chắc cho việc thiết kế, chế tạo các kết cấu tiên tiến trong công nghiệp quốc phòng, hàng không vũ trụ và năng lượng tại Việt Nam cũng như trên trường quốc tế.