Tổng quan về luận án

Sự phát triển vượt bậc của các ngành công nghiệp hàng không vũ trụ, kỹ thuật hàng hải, công nghiệp hạt nhân và quốc phòng đòi hỏi các kết cấu tấm và vỏ phải hoạt động bền bỉ trong môi trường tải trọng động phức tạp, nhiệt độ khắc nghiệt và áp suất biến đổi lớn. Vật liệu composite nhiều lớp truyền thống bộc lộ nhược điểm cố hữu là sự gián đoạn cơ tính tại mặt phân giới giữa các lớp, dễ dẫn đến hiện tượng tập trung ứng suất tiếp và gây phá hủy do tách lớp (delamination). Để khắc phục hạn chế này, vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) ra đời với đặc tính chuyển tiếp liên tục các pha thành phần gốm (ceramic) và kim loại (metal) theo chiều dày, giúp loại bỏ ứng suất tập trung bề mặt tiếp xúc và nâng cao khả năng chịu nhiệt. Luận án tiến sĩ cơ học của tác giả Lê Quang Vinh với tiêu đề "Nghiên cứu dao động tự do của kết cấu vỏ liên hợp bằng vật liệu có cơ tính biến thiên được bao quanh bởi nền đàn hồi" (Chuyên ngành: Cơ học, Mã số: 9440109, Đại học Bách khoa Hà Nội, 2022) do PGS.TS. Nguyễn Mạnh Cường và GS.TSKH. Nguyễn Đông Anh hướng dẫn đã giải quyết một mắt xích khoa học quan trọng trong việc mô hình hóa và tính toán động lực học của các hệ kết cấu phức hợp FGM trong môi trường tiếp xúc đàn hồi.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác định xuất phát từ thực tế: mặc dù dao động của vỏ tròn xoay đơn lẻ (trụ hoặc nón) làm bằng vật liệu FGM đã được một số nhóm tác giả nghiên cứu bằng phương pháp giải tích hoặc phương pháp phần tử hữu hạn (PTHH - FEM), nhưng việc phân tích động lực học dao động tự do của các hệ kết cấu vỏ FGM liên hợp phức tạp (như vỏ ghép nối nón - trụ, nón - vành - trụ, vỏ có gân gia cường dạng vành đàn hồi) đặt trong môi trường tương tác với nền đàn hồi hai hệ số Winkler-Pasternak bằng phương pháp Phần tử liên tục (Continuous Element Method - PTLT / Dynamic Stiffness Matrix Method - DSM) vẫn chưa được giải quyết một cách tường minh và hệ thống. Các phương pháp giải tích cổ điển gặp bế tắc khi xử lý hình học phức tạp, trong khi phương pháp PTHH truyền thống đòi hỏi chia lưới cực mịn, tốn kém tài nguyên tính toán và dễ phát sinh sai số tích lũy nghiêm trọng tại miền tần số trung bình và cao do hiện tượng "khóa số" (locking phenomena) và xấp xỉ chuyển vị tĩnh.

Luận án thiết lập 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học tương ứng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một mô hình phần tử liên tục chính xác cao dựa trên lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin cho các cấu trúc vỏ FGM đối xứng trục đơn lẻ (trụ, nón cụt, vành tròn) có xét đến tương tác với nền đàn hồi hai tham số?
    Giả thuyết 1 (H1): Véc-tơ trạng thái 10 thành phần kết hợp với ma trận vi phân trạng thái $10 \times 10$ cho phép thiết lập ma trận truyền và ma trận độ cứng động lực dạng giải tích - số chính xác tuyệt đối mà không cần chia nhỏ kết cấu thành nhiều phần tử hữu hạn.
  2. Câu hỏi nghiên cứu 2 (RQ2): Thuật toán ghép nối các ma trận độ cứng động lực theo cấu trúc nối tiếp (series) và cấu trúc song song phân nhánh chữ T (T-joint) có thể mô tả chính xác điều kiện tương thích chuyển vị và cân bằng nội lực tại các mặt phân giới tiếp giáp của vỏ FGM dạng bậc và vỏ có gân gia cường hay không?
    Giả thuyết 2 (H2): Ma trận độ cứng động lực tổng thể của hệ vỏ liên hợp và vỏ có gân gia cường dạng vành giải quyết triệt để bài toán biên phức tạp mà vẫn duy trì kích thước ma trận tối thiểu.
  3. Câu hỏi nghiên cứu 3 (RQ3): Phương pháp đường cong đáp ứng (Frequency Response Curve Method) có thể khắc phục được những hạn chế kinh điển của thuật toán Wittrick-Williams và việc tìm nghiệm định thức phi tuyến để trích xuất phổ tần số dao động tự do một cách ổn định hay không?
    Giả thuyết 3 (H3): Việc kích động một tải trọng xung điều hòa đơn vị $F = 1 \cdot e^{i\omega t}$ sẽ trực tiếp kích thích các đỉnh cộng hưởng trên đường cong chuyển vị - tần số, cho phép xác định chính xác các tần số riêng mà không bị bỏ sót nghiệm.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của: Lý thuyết vỏ biến dạng cắt bậc nhất Reissner-Mindlin (FSDT) với hệ số hiệu chỉnh cắt $f = 5/6$; Luật phân bố tỉ lệ thể tích 4 tham số tổng quát $V_c(z) = [1 - a(1/2 + z/h) + b(1/2 - z/h)^c]^p$; Mô hình nền đàn hồi hai tham số Winkler-Pasternak $p = k_w w - k_p \nabla^2 w$; và Phương pháp ma trận độ cứng động lực học (DSM/PTLT). Đóng góp mang tính đột phá của nghiên cứu là xây dựng thành công bộ công cụ tính toán số tin cậy trong môi trường MATLAB, chứng minh khả năng giảm thiểu số lượng phần tử xuống mức tối thiểu (1 phần tử cho mỗi phân đoạn hình học đồng nhất), tiết kiệm hơn 80% thời gian xử lý và lưu trữ so với PTHH 3D, đồng thời duy trì sai số dưới $0.5%$ so với các nghiệm giải tích đối chuẩn quốc tế. Phạm vi nghiên cứu bao quát các kết cấu vỏ tròn xoay từ vật liệu đẳng hướng kim loại đến 6 cặp tổ hợp FGM tiên tiến ($Si_3N_4/SUS304$, $ZrO_2/Ti6Al4V$, $ZrO_2/SUS304$, $Al_2O_3/Al$,...), khảo sát phổ biến thiên tham số nền $k_w \in [0, 5 \times 10^7]\text{ N/m}^3$, $k_p \in [0, 5 \times 10^6]\text{ N/m}$, góc nón $\alpha \in [0^\circ, 90^\circ]$, và số mũ thể tích $p \in [0, 10]$.

Literature Review và Positioning

Nghiên cứu động lực học kết cấu vỏ FGM trải qua quá trình tiến hóa từ các lý thuyết vỏ cổ điển đến các lý thuyết biến dạng cắt bậc cao và phương pháp số tiên tiến. Dòng nghiên cứu nền tảng khởi nguồn từ Sheng và Wang (2010, 2014), những người đã sử dụng nguyên lý Hamilton kết hợp lý thuyết phi tuyến von Kármán và lý thuyết FSDT để phân tích ứng xử dao động nhiệt - cơ của vỏ trụ FGM chứa dòng chảy chất lỏng. Song song đó, Xie và cộng sự (2014) sử dụng phương pháp truyền sóng Haar khảo sát dao động vỏ nón và vành tròn FGM. Tuy nhiên, các tiếp cận này chủ yếu giải quyết các kết cấu đơn lẻ và điều kiện biên lý tưởng.

Đối với bài toán kết cấu tương tác với môi trường nền đàn hồi, Shen và cộng sự (2010), Bagherizadeh và cộng sự (2011) đã phát triển lý thuyết biến dạng cắt bậc cao (HSDT) để khảo sát dao động và độ ổn định sau tới hạn của vỏ trụ FGM trên nền Pasternak. Trường phái nghiên cứu do Sofiyev và các cộng sự (Sofiyev, 2009; Deniz & Sofiyev, 2011; Sofiyev & Kuruoglu, 2015, 2016; Sofiyev & Schnack, 2012) dẫn đầu đã sử dụng lý thuyết vỏ mỏng Donnell và phương pháp Galerkin để giải hệ phương trình vi phân phi tuyến cho vỏ nón FGM trên nền đàn hồi hai hệ số. Một cuộc tranh luận học thuật lớn tồn tại giữa hai quan điểm:

  • Quan điểm 1: Nhóm nghiên cứu dựa trên lý thuyết vỏ cổ điển (CST) và lý thuyết Donnell (Sofiyev et al.) cho rằng phương pháp Galerkin với hàm chuyển vị xấp xỉ đơn giản là đủ để dự đoán tần số cơ bản của vỏ mỏng với chi phí tính toán thấp.
  • Quan điểm 2: Nhóm ủng hộ lý thuyết biến dạng cắt Reissner-Mindlin và HSDT (Tornabene, 2009, 2011; Sheng & Wang, 2014) phản biện rằng lý thuyết cổ điển bỏ qua biến dạng cắt ngang ($\gamma_{xz}, \gamma_{\theta z}$) và quán tính quay, dẫn đến việc đánh giá quá cao độ cứng kết cấu và gây sai số lớn (có thể vượt quá $15-20%$) khi chiều dày vỏ tăng ($h/R > 0.05$) hoặc ở các mode dao động bậc cao.

Trong phân tích kết cấu vỏ liên hợp phức tạp, các công trình quốc tế kinh điển bao gồm nghiên cứu của Caresta và Kessissoglou (2010) về vỏ ghép nối nón - trụ đẳng hướng bằng lý thuyết Flügge; Qu và cộng sự (2013, 2014) áp dụng phương pháp biến phân Fourier-Ritz cải tiến cho vỏ nón - trụ - cầu có gân gia cường; Patel và cộng sự (2006) sử dụng PTHH cho vỏ composite liên hợp; và Efraim & Eisenberger (2007) tiên phong ứng dụng ma trận độ cứng động lực cho vỏ nón - trụ đẳng hướng. Tại Việt Nam, các nhóm nghiên cứu mạnh của Đào Huy Bích, Đào Văn Dũng, Nguyễn Đình Đức, Trần Minh Tú đã công bố nhiều công trình xuất sắc về ổn định phi tuyến và động lực học của tấm/vỏ FGM có gân gia cường bằng phương pháp giải tích và phương pháp san gân (smeared stiffener technique); trong khi Trần Ích Thịnh và Nguyễn Mạnh Cường (2013-2019) phát triển phương pháp PTLT cho vỏ composite chứa chất lỏng.

TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VỎ VÀ PHƯƠNG PHÁP TÍNH TOÁN
===================================================================================
1. Lý thuyết vỏ cổ điển (CST / Donnell) [Sofiyev et al., Bích et al.]
2. Lý thuyết biến dạng cắt bậc nhất (FSDT / Reissner-Mindlin) [Tornabene, Sheng & Wang]
3. Phương pháp Phần tử hữu hạn truyền thống (FEM) [Qu et al., Patel et al.]
4. KHUNG PHÂN TÍCH LUẬN ÁN (PTLT / DSM + FSDT + Winkler-Pasternak + Ghép nối T-joint)
===================================================================================

Vị trí học thuật của luận án được xác lập rõ ràng tại giao điểm giữa lý thuyết biến dạng cắt bậc nhất FSDT và phương pháp PTLT chính xác. So với nghiên cứu của Sofiyev & Kuruoglu (2015) trên vỏ nón FGM nền Pasternak (vốn dựa vào phép xấp xỉ Galerkin trên phương trình Donnell), luận án đạt độ tổng quát hóa cao hơn khi mô tả đầy đủ hiệu ứng trượt cắt và quán tính quay. So với công trình của Qu và cộng sự (2013) dùng phương pháp biến phân Ritz (đòi hỏi lựa chọn chuỗi hàm số cơ sở phức tạp và tốn kém khi thay đổi điều kiện biên), phương pháp PTLT trong luận án cho phép xử lý mọi tổ hợp điều kiện biên (C-C, S-S, F-C, C-F) một cách tự nhiên thông qua việc áp đặt các thành phần trong véc-tơ trạng thái mà không cần định nghĩa lại hàm dạng toàn cục.

Đó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 và làm sâu sắc thêm các lý thuyết cơ học kết cấu FGM thông qua 4 đóng góp lý thuyết cốt lõi:

  1. Mở rộng lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin cho vỏ FGM đối xứng trục trên nền đàn hồi hai tham số: Thiết lập hệ phương trình vi phân chuyển động đầy đủ kết hợp tương tác màng - uốn - xoắn ($B_{ij}$) phát sinh do tính bất đối xứng của vật liệu FGM theo chiều dày. Mối liên hệ ứng suất - biến dạng tuân theo định luật Hooke tổng quát: $$\begin{pmatrix} \sigma_x \ \sigma_\theta \ \tau_{x\theta} \end{pmatrix} = \begin{bmatrix} Q_{11}(z) & Q_{12}(z) & 0 \ Q_{12}(z) & Q_{11}(z) & 0 \ 0 & 0 & Q_{66}(z) \end{bmatrix} \begin{pmatrix} \varepsilon_x \ \varepsilon_\theta \ \gamma_{x\theta} \end{pmatrix}$$ với $Q_{11}(z) = \frac{E(z)}{1-\nu^2(z)}$, $Q_{12}(z) = \frac{\nu(z)E(z)}{1-\nu^2(z)}$, $Q_{66}(z) = \frac{E(z)}{2[1+\nu(z)]}$.

  2. Mô hình hóa tương tác nền Pasternak vào động lực học vỏ: Luận án trích dẫn trực tiếp biểu thức áp lực nền:

    "Mối liên hệ giữa lực tác dụng và biến dạng của nền được biểu diễn như sau: $p = k_w w - k_p \nabla^2 w = k_w w - k_p \left(\frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2}\right)$" Việc đưa số hạng vi phân cấp hai của hệ số trượt $k_p$ vào hệ phương trình vi phân trạng thái trực tiếp làm biến đổi các phần tử của ma trận trạng thái $A_m(x,\omega)$, phản ánh chính xác hiệu ứng tăng độ cứng chống uốn cục bộ của kết cấu.

  3. Thiết lập mô hình phân bố thể tích 4 tham số tổng quát: Khảo sát cấu trúc vật liệu theo luật phân bố: $$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$$ Mô hình này bao quát cả hai dạng cấu hình $FGM_I$ (mặt đáy giàu gốm, mặt trên giàu kim loại) và $FGM_{II}$ (mặt đáy giàu kim loại, mặt trên giàu gốm), cho phép điều khiển linh hoạt gradient độ cứng và phân bố khối lượng theo phương pháp thiết kế vật liệu theo yêu cầu (tailored material design).

  4. Chuyển dịch mô hình tính toán từ phân tích xấp xỉ sang nghiệm số giải tích chính xác: Chứng minh rằng phương pháp ma trận độ cứng động lực học PTLT biểu diễn chính xác nghiệm giải tích của phương trình vi phân dao động dọc theo đường sinh của từng đoạn vỏ, loại bỏ hoàn toàn sai số xấp xỉ hàm dạng không gian vốn có của PTHH.

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết lớn: Lý thuyết môi trường liên tục đàn hồi FGM, Lý thuyết nền đàn hồi 2 tham số Pasternak, và Phương pháp ma trận độ cứng động lực dạng PTLT.

Quy trình phân tích được thiết kế nhất quán qua các điều kiện biên và mặt phân giới:

  • Điều kiện ghép nối nối tiếp (Series junction): Tại mặt tiếp giáp giữa phân đoạn 1 và phân đoạn 2 ($x = L_1$), điều kiện liên tục chuyển vị và cân bằng lực đòi hỏi: $$U_1(L_1) = U_2(0), \quad F_1(L_1) + F_2(0) = 0$$
  • Điều kiện ghép nối phân nhánh chữ T (T-junction): Áp dụng cho vỏ trụ có gân gia cường dạng vành, trong đó gân được mô hình hóa thành một phần tử vành tròn đàn hồi liên tục 2D. Tại nút ghép 3 nhánh (đoạn vỏ trước, đoạn vỏ sau và vành tròn), điều kiện tương thích và cân bằng được thiết lập: $$U_{\text{vỏ 1}} = U_{\text{vỏ 2}} = U_{\text{vành}}, \quad F_{\text{vỏ 1}} + F_{\text{vỏ 2}} + F_{\text{vành}} = 0$$
  • Điều kiện biên: Được xác định cụ thể trong miền không gian biến dạng tuyến tính bé, vật liệu đàn hồi trực hướng cục bộ, không xét đến phi tuyến hình học lớn hoặc hiện tượng từ biến, mỏi cơ học.

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ủ nghiêm ngặt quan điểm nhận thức luận của Chủ nghĩa thực chứng (Positivism) và Cơ học tính toán giải tích - số (Computational Mechanics). Bản chất của nghiên cứu là xây dựng mô hình toán học giải tích tất định từ các định luật cơ học cơ bản (nguyên lý d'Alembert, nguyên lý Hamilton) và kiểm chứng tính đúng đắn thông qua đối sánh chéo với các công trình chuẩn quốc tế.

Nghiên cứu áp dụng thiết kế đa cấp độ (Multi-level design):

  • Cấp độ vi mô: Mô hình hóa sự thay đổi liên tục của các đặc tính cơ lý $E(z), \nu(z), \rho(z)$ theo tọa độ chiều dày $z \in [-h/2, h/2]$.
  • Cấp độ trung mô: Tích phân giải tích qua chiều dày để thu nhận các ma trận độ cứng màng $[A]$, tương tác màng - uốn $[B]$, độ cứng uốn $[D]$, độ cứng cắt $[F]$ và các mômen quán tính khối lượng $I_0, I_1, I_2$: $$(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; \quad (I_0, I_1, I_2) = \int_{-h/2}^{h/2} \rho(z)(1, z, z^2)dz$$
  • Cấp độ vĩ mô: Thiết lập phương trình vi phân chủ đạo cho từng kết cấu vỏ tròn xoay phức tạp trong không gian 3 chiều đối xứng trục tương tác với môi trường nền đàn hồi.

Quy trình nghiên cứu rigorous

Luận án định nghĩa véc-tơ không gian trạng thái 10 thành phần bao gồm 5 thông số động học (chuyển vị và góc xoay) và 5 thông số lực học (lực màng, lực cắt và mô-men): $$y(x, \theta, z, t) = \left{ u_0, v_0, w_0, \varphi_x, \varphi_\theta, N_x, N_{x\theta}, Q_x, M_x, M_{x\theta} \right}^T$$

Khai triển nghiệm dao động điều hòa theo chuỗi Fourier theo phương chu vi $\theta$ với số sóng chu vi $m$: $$\begin{aligned} u_0(x,\theta,t) &= u_m(x)\cos(m\theta)e^{i\omega t}, & v_0(x,\theta,t) &= v_m(x)\sin(m\theta)e^{i\omega t}, & w_0(x,\theta,t) &= w_m(x)\cos(m\theta)e^{i\omega t} \ \varphi_x(x,\theta,t) &= \varphi_{xm}(x)\cos(m\theta)e^{i\omega t}, & \varphi_\theta(x,\theta,t) &= \varphi_{\theta m}(x)\sin(m\theta)e^{i\omega t}, & N_x(x,\theta,t) &= N_{xm}(x)\cos(m\theta)e^{i\omega t} \ N_{x\theta}(x,\theta,t) &= N_{x\theta m}(x)\sin(m\theta)e^{i\omega t}, & Q_x(x,\theta,t) &= Q_{xm}(x)\cos(m\theta)e^{i\omega t}, & M_x(x,\theta,t) &= M_{xm}(x)\cos(m\theta)e^{i\omega t} \ M_{x\theta}(x,\theta,t) &= M_{x\theta m}(x)\sin(m\theta)e^{i\omega t} \end{aligned}$$

Thay dạng nghiệm vào hệ phương trình vi phân cân bằng và quan hệ nội lực - chuyển vị, ta thu được hệ phương trình vi phân thường bậc nhất dạng ma trận: $$\frac{d y_m(x, \omega)}{dx} = A_m(x, \omega) y_m(x, \omega)$$ trong đó $A_m(x, \omega)$ là ma trận biến trạng thái kích thước $10 \times 10$ với các hệ số chứa đầy đủ các thông số hình học $R(x) = R_1 + x\sin\alpha$, thông số vật liệu và độ cứng nền đàn hồi ($k_w, k_p$).

Giải phương trình vi phân trên bằng phương pháp tích phân ma trận truyền trạng thái từ $x=0$ đến $x=L$, thu được ma trận truyền $T_m(L, \omega)$: $$y_m(L, \omega) = T_m(L, \omega) y_m(0, \omega)$$ Từ ma trận truyền $T_m$, thực hiện hoán vị và biến đổi đại số để phân tách véc-tơ trạng thái thành véc-tơ chuyển vị $U = {u_m, v_m, w_m, \varphi_{xm}, \varphi_{\theta m}}^T$ và véc-tơ lực $F = {N_{xm}, N_{x\theta m}, Q_{xm}, M_{xm}, M_{x\theta m}}^T$ tại hai đầu biên, từ đó thiết lập ma trận độ cứng động lực phần tử liên tục: $$\begin{pmatrix} F(0) \ F(L) \end{pmatrix} = [K(\omega)] \begin{pmatrix} U(0) \ U(L) \end{pmatrix}$$

Data và phân tích

Toàn bộ thuật toán được lập trình tự động hóa trên nền tảng phần mềm MATLAB. Thay vì sử dụng thuật toán Wittrick-Williams (vốn gặp khó khăn lớn về thời gian tính toán và điểm kỳ dị hàm siêu việt khi giải các hệ vỏ phức tạp) hoặc giải phương trình định thức phi tuyến $\det|K(\omega)| = 0$, luận án áp dụng xuất sắc Phương pháp đường cong đáp ứng tần số. Luận án trích dẫn nguyên lý then chốt:

"Một kích động xung bằng búa (mô phỏng tải trọng va chạm) qua phân tích Fourier theo lý thuyết sẽ kích thích tất cả các tần số dao động riêng của kết cấu... Tại các tần số riêng xảy ra cộng hưởng khiến cho chuyển vị và điện áp tăng vọt. Từ đường cong đáp ứng thu được giống như trên các máy đo ta xác định được các giá trị tần số dao động cần tìm là tọa độ các đỉnh của đường cong đáp ứng."

Thuật toán đặt một lực kích động điều hòa đơn vị tại một điểm biên $x=L$: $Q_x(L) = \frac{1}{R}\cos(m\theta)e^{i\omega t}$. Bằng cách quét dải tần số $\omega$ với bước tăng mịn ($\Delta\omega \le 0.1\text{ rad/s}$), hệ phương trình đại số tuyến tính $F(\omega) = [K(\omega)]U(\omega)$ được giải trực tiếp để tìm chuyển vị pháp tuyến $w(\omega)$. Tần số riêng chính là các hoành độ cực đại (cộng hưởng) trên đồ thị đáp ứng biên độ - tần số.

Tính chuẩn xác và độ tin cậy của phương pháp được kiểm chứng thông qua các bộ dữ liệu đối chuẩn:

  • Đặc tính vật liệu chuẩn:
    • Thép không gỉ (SUS304): $E_m = 207.7 \times 10^9\text{ N/m}^2$, $\nu_m = 0.3177$, $\rho_m = 8166\text{ kg/m}^3$.
    • Gốm Silicon Nitride ($Si_3N_4$): $E_c = 322.2 \times 10^9\text{ N/m}^2$, $\nu_c = 0.24$, $\rho_c = 2370\text{ kg/m}^3$.
    • Nhôm (Al): $E_m = 70.0 \times 10^9\text{ N/m}^2$, $\nu_m = 0.3$, $\rho_m = 2707\text{ kg/m}^3$.
    • Gốm Alumina ($Al_2O_3$): $E_c = 380 \times 10^9\text{ N/m}^2$, $\nu_c = 0.26$, $\rho_c = 3750\text{ kg/m}^3$.
    • Zirconia ($ZrO_2$): $E_c = 151.0 \times 10^9\text{ N/m}^2$, $\nu_c = 0.3$, $\rho_c = 3000\text{ kg/m}^3$.
    • Hợp kim Titanium ($Ti6Al4V$): $E_m = 105.7 \times 10^9\text{ N/m}^2$, $\nu_m = 0.298$, $\rho_m = 4429\text{ kg/m}^3$.
  • Độ tin cậy số liệu: So sánh trực tiếp kết quả tính toán tần số không thứ nguyên $\Omega = \omega R \sqrt{\rho_m(1-\nu_m^2)/E_m}$ của luận án với các công bố độc lập của Sofiyev & Schnack (2012) và Qu et al. (2013) cho thấy mức độ sai lệch tuyệt đối đều nhỏ hơn $0.35%$, khẳng định độ hội tụ và độ chính xác vượt trội của thuật toán PTLT.

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

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

Kết quả mô phỏng số và phân tích tham số toàn diện trong luận án đã mang lại 5 phát hiện mang tính bước ngoặt:

  1. Quy luật chi phối của nền đàn hồi Winkler-Pasternak lên phổ tần số riêng: Nền đàn hồi làm tăng đáng kể tần số dao động tự do của mọi loại kết cấu vỏ FGM do làm tăng độ cứng tương đương của toàn hệ thống. Đáng chú ý, hệ số cắt Pasternak $k_p$ có tác động làm tăng tần số mạnh hơn nhiều lần so với hệ số Winkler $k_w$ trên cùng một bậc độ lớn đơn vị chuẩn hóa. Khi $k_w$ tăng từ $0$ lên $5 \times 10^7\text{ N/m}^3$, tần số riêng cơ bản tăng từ $15%$ đến $35%$; tuy nhiên khi $k_p$ tăng từ $0$ lên $5 \times 10^6\text{ N/m}$, tần số dao động riêng có thể tăng vọt từ $50%$ đến hơn $120%$, đặc biệt ở các mode sóng bậc thấp.

  2. Ảnh hưởng phi tuyến của số mũ tỉ lệ thể tích $p$ và kiểu hàm phân bố ($FGM_I$ vs $FGM_{II}$): Khi số mũ thể tích $p$ tăng từ $0$ (100% gốm) đến vô cùng (100% kim loại), tần số dao động riêng giảm đơn điệu do mô đun đàn hồi của gốm lớn hơn kim loại ($E_c > E_m$). Tuy nhiên, tốc độ suy giảm tần số diễn ra mạnh nhất trong khoảng $0 \le p \le 2$ (chiếm hơn $75%$ tổng mức suy giảm độ cứng). Kiểu hàm $FGM_{II}$ (mặt ngoài giàu gốm) cho tần số dao động riêng cao hơn kiểu hàm $FGM_I$ khi kết cấu chịu tương tác nền bên ngoài, do lớp gốm cứng tiếp xúc trực tiếp với nền đàn hồi tạo ra hiệu ứng gia cường tiếp xúc vượt trội.

  3. Hiện tượng dịch chuyển mode dao động cơ bản (Mode switching): Ở các kết cấu vỏ nón và vỏ nón - trụ liên hợp, mode dao động ứng với tần số nhỏ nhất không cố định ở số sóng chu vi $m=1$ mà thay đổi phụ thuộc vào tỉ số hình học $L/R_1$, góc nón $\alpha$ và độ cứng nền. Khi góc nửa đỉnh nón $\alpha$ tăng từ $10^\circ$ lên $60^\circ$, tần số dao động tự do giảm mạnh do độ cứng chống uốn tổng thể của vỏ nón giảm theo phương dọc trục.

  4. Hiệu quả gia cường đột phá của gân dạng vành trên vỏ trụ FGM: Việc bổ sung các gân gia cường dạng vành (được mô hình hóa chính xác bằng phần tử vành đàn hồi 2D thông qua thuật toán ghép nối chữ T) làm tăng vọt tần số dao động riêng cơ bản (tăng từ $40%$ đến hơn $200%$ tùy thuộc vào số lượng và kích thước tiết diện gân), đồng thời chia nhỏ các dạng dao động cục bộ, ngăn ngừa hiệu quả hiện tượng mất ổn định động lực học dạng móp cục bộ.

  5. Sự vượt trội của phương pháp PTLT trên dải tần số cao: Khác với phương pháp PTHH truyền thống bị suy giảm độ chính xác nghiêm trọng khi tính toán các mode dao động bậc cao ($m \ge 8, n \ge 5$), phương pháp PTLT kết hợp đường cong đáp ứng duy trì độ chính xác ổn định tuyệt đối trên toàn bộ dải tần số từ 0 Hz đến trên 5000 Hz mà không cần tăng số lượng phần tử.

ẢNH HƯỞNG CỦA CÁC THAM SỐ ĐẾN TẦN SỐ DAO ĐỘNG RIÊNG (ω)
===================================================================================
Tham số tác động           Chiều hướng biến thiên      Mức độ tác động lên tần số riêng
-----------------------------------------------------------------------------------
Độ cứng nền Pasternak (k_p) Tăng (0 -> 5x10^6 N/m)      TĂNG RẤT MẠNH (+50% đến +120%)
Độ cứng nền Winkler (k_w)   Tăng (0 -> 5x10^7 N/m³)     TĂNG VỪA PHẢI (+15% đến +35%)
Số mũ thể tích FGM (p)      Tăng (0 -> 10)              GIẢM PHI TUYẾN (-20% đến -45%)
Góc nửa đỉnh nón (α)        Tăng (0° -> 60°)            GIẢM ĐƠN ĐIỆU (-30% đến -60%)
Gân gia cường vành (T-joint) Bổ sung gân đàn hồi        TĂNG ĐỘT BIẾN (+40% đến +200%)
Điều kiện biên (BCs)        Từ F-C -> S-S -> C-C        TĂNG BẬC THANG (C-C cao nhất)
===================================================================================

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp một khung toán học hoàn chỉnh, chính xác để phân tích động lực học kết cấu vỏ FGM tròn xoay, bổ sung vào kho tàng cơ học vật rắn biến dạng một phương pháp tiếp cận thanh lịch kết hợp ưu điểm của giải tích và số trị.
  • Về mặt phương pháp luận: Thuật toán ghép nối ma trận chữ T và phương pháp đường cong đáp ứng mở ra hướng tiếp cận tiêu chuẩn mới để xử lý các bài toán cơ học kết cấu dị hướng phức tạp, có thể mở rộng cho vật liệu áp điện FGM, vật liệu xốp có lỗ rỗng (porous FGM) hoặc kết cấu composite gia cường ống nano carbon (CNTRC).
  • Về mặt ứng dụng thực tiễn: Các biểu đồ và bảng dữ liệu số trong luận án cung cấp căn cứ kỹ thuật chính xác cho các kỹ sư thiết kế trong việc lựa chọn tối ưu hàm tỉ lệ thể tích $p$, góc côn $\alpha$, cấu hình bố trí gân và vật liệu nền nhằm chủ động điều khiển phổ tần số riêng, tránh hoàn toàn hiện tượng cộng hưởng phá hủy trong động cơ phản lực, thân vỏ tàu ngầm và ống dẫn áp lực cao.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu chính:

  1. Giới hạn giả thiết biến dạng bé tuyến tính: Nghiên cứu dừng lại ở bài toán dao động tự do đàn hồi tuyến tính biên độ bé, chưa xét đến hiệu ứng phi tuyến hình học von Kármán khi kết cấu chịu dao động biên độ lớn.
  2. Môi trường nhiệt độ đẳng nhiệt: Các đặc tính vật liệu FGM được coi là hằng số theo nhiệt độ trong quá trình dao động, chưa xét đến sự suy giảm mô đun đàn hồi khi làm việc trong trường nhiệt độ cao không đồng nhất ($E = E(z, T)$).
  3. Mô hình nền đàn hồi lý tưởng: Nền đàn hồi được giả thiết là đẳng hướng, đàn hồi tuyến tính thuần nhất dọc theo chiều dài, chưa mô hình hóa các trạng thái nền phi tuyến, nền viscoelastic (nhớt đàn hồi) hoặc nền rỗng cục bộ.
  4. Phương diện thực nghiệm: Các kết quả của luận án chủ yếu được kiểm chứng thông qua đối chuẩn số và giải tích với các công bố quốc tế uy tín, chưa tiến hành chế tạo mẫu thực nghiệm FGM và đo đạc dao động trực tiếp bằng máy bắn laser Doppler do điều kiện chế tạo FGM phức tạp.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng tiếp tục phát triển:

  • Mở rộng phương pháp PTLT để phân tích dao động phi tuyến hình học và phi tuyến vật liệu của vỏ liên hợp FGM.
  • Khảo sát bài toán tương tác đa trường Vật lý (Nhiệt - Cơ - Điện - Từ) cho vỏ FGM thông minh (Piezoelectric-FGM / Magneto-electro-elastic).
  • Nghiên cứu động lực học kết cấu vỏ FGM liên hợp tiếp xúc đồng thời với nền đàn hồi và chứa chất lỏng chuyển động bên trong (Fluid-Structure-Foundation Interaction).
  • Phát triển thuật toán tối ưu hóa đa mục tiêu (Multi-objective optimization) hình học và phân bố vật liệu để giảm thiểu khối lượng và triệt tiêu rung động.

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

Luận án tạo ra tác động học thuật và công nghệ sâu rộng:

  • Ảnh hưởng học thuật: Thiết lập chuẩn đối sánh số (benchmark numerical data) chất lượng cao cho cộng đồng cơ học tính toán tại Việt Nam và quốc tế; các công trình trích xuất từ luận án công bố trên các tạp chí chuyên ngành uy tín khẳng định năng lực nghiên cứu độc lập và trình độ tiệm cận quốc tế.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp vào quy trình tính toán thiết kế vỏ buồng đốt động cơ tuabin khí, thân vỏ tên lửa đẩy nhiên liệu hạt nhân, vỏ tàu ngầm chịu áp suất sâu và hệ thống đường ống công nghiệp hóa chất - dầu khí.
  • Lợi ích xã hội và an toàn công trình: Giúp nâng cao tuổi thọ công trình, ngăn ngừa thảm họa do mỏi dao động và cộng hưởng cơ học, tiết kiệm hàng triệu USD chi phí thử nghiệm thực địa nhờ độ tin cậy vượt trội của mô phỏng số.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận một tài liệu mẫu mực về phương pháp ma trận độ cứng động lực PTLT, lý thuyết biến dạng cắt bậc nhất FSDT và kỹ thuật lập trình số MATLAB chuyên sâu.
  • Các nhà khoa học và giảng viên cơ học kết cấu: Kế thừa thư viện phần tử liên tục mới và thuật toán ghép nối chữ T để giảng dạy và phát triển các bài toán dao động phức tạp hơn.
  • Kỹ sư R&D công nghiệp Hàng không, Hàng hải và Quốc phòng: Ứng dụng các công thức tính toán và dữ liệu tham số để tối ưu hóa thiết kế vỏ bọc FGM chịu tải trọng động khắc nghiệt.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất của luận án là việc xây dựng thành công mô hình 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 (FSDT) cho toàn bộ hệ thống vỏ tròn xoay liên hợp FGM (nón - trụ, nón - vành - trụ, vỏ có gân gia cường) tương tác với nền đàn hồi hai tham số Winkler-Pasternak. Nghiên cứu đã mở rộng lý thuyết vỏ FSDT từ bài toán vỏ đồng nhất/composite sang môi trường vật liệu có cơ tính biến thiên dị hướng liên tục theo luật 4 tham số, tích hợp đồng thời hiệu ứng tương tác màng - uốn - xoắn ($B_{ij}$) và áp lực trượt vi phân cấp hai của nền Pasternak vào ma trận vi phân trạng thái $10 \times 10$.

2. Sự đổi mới về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

So với nghiên cứu của Sofiyev & Kuruoglu (2015) (sử dụng lý thuyết vỏ mỏng cổ điển Donnell và phương pháp xấp xỉ Galerkin) và nghiên cứu của Qu và cộng sự (2013) (sử dụng phương pháp biến phân Fourier-Ritz cải tiến với lưới phần tử hữu hạn), phương pháp luận của luận án thể hiện 2 đột phá vượt bậc:

  • Thứ nhất, phương pháp PTLT của luận án không sử dụng hàm dạng xấp xỉ tĩnh mà giải trực tiếp nghiệm động lực học giải tích qua ma trận truyền $T(x,\omega)$, đảm bảo độ chính xác tuyệt đối mà chỉ cần duy nhất 1 phần tử liên tục cho mỗi phân đoạn hình học.
  • Thứ hai, thuật toán ghép nối song song dạng chữ T (T-joint) độc đáo cho phép mô hình hóa gân gia cường như một kết cấu vành tròn đàn hồi 2D thực thụ, khắc phục triệt để nhược điểm làm mịn gân giả định (smeared stiffener) của các phương pháp truyền thống.

3. Phát hiện bất ngờ nhất trong dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là vai trò chi phối áp đảo của hệ số trượt Pasternak $k_p$ so với hệ số Winkler $k_w$, cùng hiện tượng trượt mode dao động cơ bản ($m$) khi thay đổi độ cứng nền. Trong khi trực giác thông thường kỳ vọng hệ số đàn hồi đàn lò xo pháp tuyến $k_w$ sẽ đóng vai trò chính, số liệu tính toán chứng minh hệ số cắt $k_p$ làm tăng độ cứng chống uốn tổng thể mạnh gấp nhiều lần và làm thay đổi trật tự xuất hiện của các mode dao động riêng theo phương chu vi.

4. Luận án có cung cấp quy trình tái lập kết quả (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ và chi tiết. Luận án cung cấp tường minh:

  • Toàn bộ 100 phần tử giải tích của ma trận biến trạng thái $A_m(x,\omega)$ kích thước $10 \times 10$.
  • Các công thức tích phân xác định độ cứng $A_{ij}, B_{ij}, D_{ij}, F_{ij}$ và quán tính $I_0, I_1, I_2$.
  • Thuật toán chuyển đổi từ ma trận truyền $T_m$ sang ma trận độ cứng động lực $[K(\omega)]$.
  • Mã giả và quy trình quét tần số của phương pháp đường cong đáp ứng trong môi trường MATLAB, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập chính xác $100%$ các kết quả số đã công bố.

5. Kế hoạch nghiên cứu 10 năm được định hình như thế nào từ kết quả luận án?

Kế hoạch 10 năm tập trung vào 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1-3 năm): Tích hợp trường nhiệt độ phi tuyến và phân tích ổn định động lực học phi tuyến biên độ lớn cho vỏ FGM liên hợp.
  • Giai đoạn 2 (4-6 năm): Phát triển mô hình tương tác đa trường (Multiphysics) cho vật liệu FGM áp điện và composite nano xốp trên nền đàn hồi phi tuyến nhớt (viscoelastic foundation).
  • Giai đoạn 3 (7-10 năm): Chế tạo mẫu thực nghiệm FGM bằng công nghệ in 3D kim loại - gốm tiên tiến và thực hiện đo đạc dao động Laser Doppler Vibrometry để thương mại hóa gói phần mềm PTLT chuyên dụng cho công nghiệp quốc phòng và hàng không vũ trụ.

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, xuất sắc và có tính hệ thống cao, đóng góp những giá trị cốt lõi cho ngành Cơ học vật rắn biến dạng:

  1. Xây dựng thành công hệ thống lý thuyết và thư viện Phần tử liên tục (PTLT) mới chính xác cao cho các kết cấu vỏ tròn xoay đơn lẻ và liên hợp FGM tương tác với nền đàn hồi hai tham số Winkler-Pasternak.
  2. Phát triển 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, giải quyết triệt để bài toán biên phức tạp cho vỏ nón - trụ, nón - vành - trụ và vỏ trụ có gân gia cường dạng vành đàn hồi.
  3. Ứng dụng thành công phương pháp đường cong đáp ứng tần số, loại bỏ những khó khăn kinh điển của thuật toán Wittrick-Williams và phương pháp tìm nghiệm định thức phi tuyến.
  4. Khám phá và định lượng hóa toàn diện quy luật tác động của các tham số vật liệu FGM 4 hệ số, đặc tính nền Pasternak, góc nón và điều kiện biên lên phổ tần số tự do.
  5. Cung cấp bộ dữ liệu đối chuẩn tin cậy cao, mở ra 3 hướng nghiên cứu mới về dao động phi tuyến, tương tác đa trường và tối ưu hóa kết cấu FGM tiên tiến.

Công trình khẳng định bước tiến vượt bậc của cơ học tính toán Việt Nam, đóng góp thiết thực cho sự nghiệp công nghiệp hóa, hiện đại hóa và tự chủ công nghệ kỹ thuật cao.