Tổng quan về luận án

Nghiên cứu của Tiến sĩ Nguyễn Bá Duy (2019) tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh với đề tài "Phân tích ứng xử của dầm sandwich chức năng chịu tác dụng của tải trọng cơ thủy nhiệt" (Chuyên ngành Cơ kỹ thuật - Mã số: 9520101) thuộc nhóm công trình tiên phong trong cơ học tính toán vật liệu composite tiên tiến. Trong bối cảnh các kết cấu hàng không - vũ trụ, tàu con thoi tái nhập khí quyển (re-entry vehicles), động cơ phản lực và các hệ thống vi cơ điện tử (MEMS/NEMS) phải làm việc trong điều kiện môi trường khắc nghiệt chịu đồng thời các trường nhiệt độ cao, độ ẩm biến thiên và tải trọng cơ học động, việc tối ưu hóa khả năng chịu lực và chống suy giảm liên kết liên lớp (delamination) là bài toán sống còn.

Khoảng trống học thuật (Research Gap) then chốt mà luận án tập trung giải quyết xuất phát từ các hạn chế nội tại của hệ thống lý thuyết dầm hiện hành:

  1. Lý thuyết dầm cổ điển (Classical Beam Theory - CBT / Euler–Bernoulli) hoàn toàn bỏ qua biến dạng trượt ngang ($\gamma_{xz} = 0$), chỉ phù hợp với kết cấu dầm rất mảnh ($L/h \ge 100$).
  2. Lý thuyết biến dạng trượt bậc nhất (First-Order Shear Deformation Theory - FSDT / Timoshenko) xem ứng suất trượt không đổi theo chiều dày, đòi hỏi hệ số hiệu chỉnh trượt ($K_s \approx 5/6$) vốn phụ thuộc phức tạp vào dạng hình học mặt cắt, điều kiện biên và quy luật biến tính vật liệu.
  3. Các lý thuyết biến dạng trượt bậc cao (Higher-Order Shear Deformation Theories - HSDTs) hiện hữu thường làm gia tăng số lượng biến số trường chuyển vị, phức tạp hóa thuật toán hoặc chưa giải quyết triệt để hiệu ứng kéo/nén theo phương chiều dày ($\varepsilon_z \neq 0$ - Transverse Normal Strain Effect).
  4. Sự thiếu vắng một mô hình giải tích - số thống nhất, có khả năng tích hợp đồng thời hiệu ứng phụ thuộc kích thước (Size-dependent effects) theo lý thuyết ứng suất cặp hiệu chỉnh (MCST) và lý thuyết phi cục bộ Eringen với quy luật truyền nhiệt phi tuyến Fourier và khuếch tán ẩm phi tuyến trong dầm sandwich FGM nhiều lớp với các điều kiện biên tùy ý.

Luận án thiết lập 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Làm thế nào để xây dựng một lý thuyết dầm biến dạng trượt bậc cao thống nhất và mô hình Quasi-3D thu giảm biến số nhằm mô tả chính xác trường ứng suất trượt ngang mà không cần hệ số hiệu chỉnh?
  • RQ2: Sự phân bố nhiệt độ và độ ẩm phi tuyến (Nonlinear Hygro-thermal distribution) theo định luật Fourier tác động như thế nào đến độ võng tĩnh, tải trọng mất ổn định nhiệt tới hạn ($N_{cr}$, $\bar{T}_{cr}$) và tần số dao động riêng ($\bar{\omega}$)?
  • RQ3: Ở quy mô micro và nano, các thông số chiều dài đặc trưng vật liệu (MLSPs - $l$) và tham số phi cục bộ ($\mu = (e_0 a)^2$) chi phối đáp ứng uốn, dao động và mất ổn định nhiệt của dầm FGM ra sao dưới các điều kiện biên khác nhau?
  • RQ4: Làm thế nào để thiết lập thuật toán số hữu hiệu (Ritz lai ghép và Phần tử hữu hạn - FEM) đảm bảo sự hội tụ nhanh và triệt tiêu hiện tượng khóa trượt (shear locking)?

Khung lý thuyết của luận án tích hợp: Cơ học môi trường liên tục phi cổ điển (Eringen's Nonlocal Elasticity, Yang's Modified Couple Stress Theory), Lý thuyết đàn hồi 2D vi phân, Nguyên lý biến phân Hamilton và Cơ học vật liệu chức năng biến tính liên tục (FGMs). Phạm vi nghiên cứu bao quát 4 cấu hình kết cấu dầm: Dầm FGM đơn tầng (Type A), Dầm sandwich FGM có các lớp bề mặt biến tính - lõi đẳng tính (Type B), Dầm sandwich có lõi biến tính - lớp bề mặt đẳng tính (Type C), và Dầm FGM biến tính 2 chiều (Type D - 2D FGM). Các mẫu khảo sát sử dụng hệ vật liệu gốm-kim loại điển hình như $Si_3N_4/SUS304$ và $Al/Al_2O_3$ với tỷ số kích thước $L/h$ biến thiên từ 5 đến 100, chỉ số lũy thừa $p \in [0, 10]$, qua đó xác lập các đóng góp đột phá về mặt lý thuyết và công cụ tính toán chính xác cao.

Literature Review và Positioning

Cơ sở lý thuyết về dầm vật liệu composite và dầm FGM đã trải qua quá trình tiến hóa sâu rộng. Ghugal và Shimpi (2001), Sayyad và Ghugal (2015, 2017) đã tổng thuật toàn diện các mô hình dầm từ cổ điển đến bậc cao. Khi nghiên cứu vật liệu FGM – khởi xướng từ năm 1984 bởi các nhà khoa học vật liệu tại Sendai, Nhật Bản (Koizumi, 1997; Yamanouchi et al.) – các thuộc tính đàn hồi được làm mịn liên tục để loại bỏ sự gián đoạn mặt phân giới vốn gây tập trung ứng suất và nứt tách lớp trong composite truyền thống (Laminated composites).

Quy luật phân bố thể tích thành phần vật liệu $V(z)$ qua chiều dày được mô hình hóa qua các quy tắc:

  1. Quy luật hàm số mũ (Power-law): $$P(z) = (P_c - P_m)V(z) + P_m, \quad V(z) = \left(\frac{2z + h}{2h}\right)^p$$
  2. Quy luật hàm mũ thực (Exponential-law) của Delale và Erdogan (61): $P(z) = A e^{B(z + h/2)}$ với $A = P_m, B = \frac{1}{h}\ln\left(\frac{P_c}{P_m}\right)$.
  3. Quy luật Sigmoid (Chung và Chi, 2001) nhằm giảm ứng suất dư đột ngột gần hai bề mặt tiếp xúc.
       TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT BIẾN DẠNG TRƯỢT
 Euler-Bernoulli (CBT)       Timoshenko (FSDT)             Reddy (TSDT) / Touratier
  (Bỏ qua trượt ngang)    (Ứng suất trượt hằng số)       (Hàm đa thức/lượng giác)
  Dầm rất mảnh (L/h>100)    Cần hệ số Ks (Ks=5/6)        Triệt tiêu ứng suất mặt biên
                  Luận án Nguyễn Bá Duy (2019)
Mô hình Unified HSBT & Quasi-3D                   Mô hình Micro/Nano Size-Dependent
(Hàm dạng Hyperbolic-Sine mới: f(z))              (Tích hợp MCST và Nonlocal Eringen)

Về mặt động học dầm, các lý thuyết trượt bậc cao phi đa thức được đề xuất nhằm mô tả chính xác đường cong phân bố ứng suất tiếp $\tau_{xz}$ triệt tiêu tại hai bề mặt tự do ($z = \pm h/2$):

  • Reddy (1984) sử dụng hàm đa thức bậc 3: $f(z) = z\left[1 - \frac{4z^2}{3h^2}\right]$
  • Touratier (1991) đề xuất hàm lượng giác: $f(z) = \frac{h}{\pi}\sin\left(\frac{\pi z}{h}\right)$
  • Karama et al. (2003) ứng dụng hàm mũ: $f(z) = z e^{-2(z/h)^2}$
  • Mantari et al. (2012) và Nguyen et al. (2015) tiếp tục mở rộng các hàm tiếp tuyến và lượng giác ngược.

Tuy nhiên, các tranh biện học thuật lớn vẫn tồn tại:

  • Tranh biện 1 (Polynomial vs. Non-polynomial & Quasi-3D): Các lý thuyết HSDT truyền thống chỉ xét trường chuyển vị trục $u(x,z,t)$ bậc cao nhưng gán chuyển vị đứng $w(x,t) = w_0(x,t)$, tức bỏ qua ứng suất pháp ngang $\sigma_z$ ($\varepsilon_z = 0$). Các công trình của Zenkour (2006) và Vo et al. chỉ ra rằng đối với dầm dày ($L/h \le 5$) hoặc dầm sandwich có lõi mềm, hiệu ứng co giãn chiều dày (thickness stretching) đóng vai trò then chốt trong việc dự báo chính xác tần số dao động tự nhiên và tải trọng mất ổn định.
  • Tranh biện 2 (Continuum cổ điển vs. Hiệu ứng phụ thuộc kích thước): Khi kích thước cấu kiện giảm xuống cấp độ micro/nano (MEMS/NEMS), lý thuyết cơ học cổ điển mất khả năng nắm bắt hiệu ứng kích thước mạng tinh thể. Eringen (1972, 1983) đề xuất lý thuyết đàn hồi phi cục bộ xem ứng suất tại một điểm phụ thuộc vào trường biến dạng của toàn bộ miền lân cận. Trong khi đó, Yang et al. (2002) phát triển lý thuyết ứng suất cặp hiệu chỉnh (MCST) dựa trên tensor độ cong đối xứng và chỉ cần duy nhất một thông số chiều dài đặc trưng vật liệu ($l$), giải quyết triệt để tính mơ hồ của lý thuyết ứng suất cặp cổ điển (Mindlin, 1964; Toupin, 1962).

So sánh với hai nghiên cứu quốc tế điển hình:

  1. So với nghiên cứu của Simsek (2010, 2016) về dầm FGM micro/nano: Simsek chủ yếu áp dụng lý thuyết Timoshenko hoặc HSDT đa thức cho dầm đơn lớp đồng chất/FGM đơn giản. Luận án của Nguyễn Bá Duy đã mở rộng vượt bậc khi thiết lập lời giải cho dầm sandwich nhiều lớp (Type B, C, D) dưới tải trọng liên hợp cơ - thủy - nhiệt phi tuyến.
  2. So với nghiên cứu của Zenkour (2006) và Pradhan & Murmu (2009): Các tác giả trên chủ yếu giải quyết bài toán bằng phương pháp Navier vốn bị đóng khung trong điều kiện biên tựa đơn (Simply Supported - S-S). Luận án của Nguyễn Bá Duy xây dựng hệ hàm dạng Ritz lai ghép mới (Novel Hybrid Shape Functions) và mô hình FEM tương thích, cho phép giải quyết chính xác cho toàn bộ các tập điều kiện biên phức tạp: Ngàm-Ngàm (C-C), Ngàm-Tự do (C-F), Ngàm-Khớp (C-H), Khớp-Khớp (H-H).

Đó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à phát triển các nền tảng lý thuyết cơ học kết cấu thông qua các đột phá cụ thể:

  1. Thiết lập mô hình lý thuyết dầm biến dạng trượt bậc cao thống nhất (Unified Higher-Order Shear Deformation Beam Theory - Unified HSBT): Xuất phát trực tiếp từ các phương trình cơ bản của lý thuyết đàn hồi 2D, luận án xây dựng biểu thức tổng quát cho trường chuyển vị: $$u(x,z,t) = u_0(x,t) - z w_{0,x}(x,t) + f(z)\theta(x,t)$$ $$w(x,z,t) = w_0(x,t)$$ Đồng thời, luận án phát triển mô hình dầm tinh chế (Refined HSBT) phân rã chuyển vị ngang thành thành phần uốn $w_b(x,t)$ và thành phần trượt $w_s(x,t)$: $$u(x,z,t) = u_0(x,t) - z w_{b,x}(x,t) - [f(z) - z] w_{s,x}(x,t)$$ $$w(x,z,t) = w_b(x,t) + w_s(x,t)$$

  2. Đề xuất hàm dạng trượt bậc cao Hyperbolic-Sine mới: Để thỏa mãn điều kiện biên ứng suất tiếp bằng không tại hai mặt tự do $\tau_{xz}(x, \pm h/2) = 0$ mà không cần bất kỳ hệ số hiệu chỉnh nào, luận án công bố hàm dạng nguyên bản: $$f(z) = \sinh\left(\frac{rz}{h}\right) - \frac{8rz^3}{3h^3(r^2 - 4)}$$ Đạo hàm tương ứng của hàm trượt có dạng: $$g(z) = f'(z) = \frac{r}{h}\cosh\left(\frac{rz}{h}\right) - \frac{8rz^2}{h^3(r^2 - 4)}$$ Trong đó, $r$ là tham số điều chỉnh (correction parameter) cho phép tối ưu hóa trường ứng suất phù hợp với sự biến thiên phi đối xứng của gradient vật liệu FGM qua chiều dày.

       QUY CÁCH PHÂN TÍCH LIÊN HỢP ĐA TRƯỜNG VÀ ĐA TỶ LỆ
  1. Xây dựng lý thuyết dầm Quasi-3D ba biến số mới (Novel Three-Variable Quasi-3D Theory): Bổ sung hàm biến dạng đứng $g(z)$ để phản ánh độ giãn nở bề dày: $$u(x,z,t) = u_0(x,t) - z w_{0,x} + f(z)\theta(x,t)$$ $$w(x,z,t) = w_0(x,t) + g(z)w_z(x,t)$$ Mô hình này giúp thu giữ đồng thời ứng suất pháp ngang $\sigma_z$ và tương tác trượt $\tau_{xz}$ trong một cấu trúc phương trình cô đọng, giảm thiểu số ẩn so với mô hình Quasi-3D 5 biến truyền thống nhưng vẫn đạt độ chính xác tương đương nghiệm giải tích đàn hồi 3D.

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

Khung phân tích của luận án tích hợp liền mạch 4 trụ cột cơ học:

  • Cơ học biến dạng liên hợp Thủy - Nhiệt: Khảo sát trường nhiệt độ từ phương trình dẫn nhiệt Fourier dừng một chiều: $$T(z) = T_b + (T_t - T_b)\frac{\int_{-h/2}^z \frac{1}{k(z)}dz}{\int_{-h/2}^{h/2} \frac{1}{k(z)}dz}$$ kết hợp biến thiên độ ẩm dạng hình sin (Sinusoidal non-linear rise).
  • Cơ học vi mô MCST: Năng lượng biến dạng $U_M$ bao gồm cả tensor ứng suất đối xứng $\sigma$ và tensor ứng suất cặp lệch $m$: $$U_M = \frac{1}{2}\int_V (\sigma : \varepsilon + m : \chi) dV$$ với tensor độ cong đối xứng $\chi_{ij} = \frac{1}{2}(\theta_{i,j} + \theta_{j,i})$, vector góc xoay $\theta = \frac{1}{2}\text{curl}(u)$.
  • Cơ học kích thước Nano Eringen: Phương trình cấu lượng phi cục bộ: $$(1 - \mu \nabla^2)\sigma_{xx} = E(z)\varepsilon_{xx}, \quad (1 - \mu \nabla^2)\tau_{xz} = G(z)\gamma_{xz}$$
  • Điều kiện biên xác định (Boundary Conditions): Thiết lập chặt chẽ cho toàn bộ 6 tổ hợp biên: Simply Supported (S-S), Clamped-Clamped (C-C), Clamped-Free (C-F), Clamped-Hinged (C-H), Hinged-Hinged (H-H), Clamped-Simply Supported (C-S).

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

Thiết kế nghiên cứu

Nghiên cứu được định vị vững chắc theo hệ giác ngộ thực chứng (Positivism) kết hợp với chủ nghĩa hiện thực phản biện (Critical Realism) trong cơ học tính toán. Toàn bộ các mô hình toán học đều được dẫn xuất từ các định luật cơ bản của cơ học giải tích (Nguyên lý công ảo, Nguyên lý Hamilton) để thiết lập hệ phương trình vi phân chuyển động đạo hàm riêng: $$\int_{t_1}^{t_2} (\delta U - \delta K + \delta V) dt = 0$$ trong đó $\delta U$ là biến phân nội năng đàn hồi (bao gồm thế năng nhiệt - ẩm và hiệu ứng couple stress), $\delta K$ là biến phân động năng, và $\delta V$ là biến phân công ngoại lực.

                  QUY TRÌNH PHÂN TÍCH SỐ VÀ ĐỐI SOÁT KIỂM CHỨNG

Quy trình nghiên cứu rigorous

Quy trình giải tích số được chuẩn hóa qua các bước nghiêm ngặt:

  1. Thiết lập chuỗi hàm dạng Ritz lai ghép mới (Novel Hybrid Shape Functions): Để vượt qua nhược điểm không thỏa mãn điều kiện biên hình học phức tạp của chuỗi lượng giác Navier thuần túy, luận án xây dựng hàm dạng Ritz: $$u_0(x,t) = \sum_{j=1}^N \phi_j^u(x) u_{0j} e^{i\omega t}, \quad w_0(x,t) = \sum_{j=1}^N \phi_j^w(x) w_{0j} e^{i\omega t}, \quad \theta(x,t) = \sum_{j=1}^N \phi_j^\theta(x) \theta_j e^{i\omega t}$$ Các hàm cơ sở $\phi_j(x)$ được xây dựng bằng cách nhân các hàm biên cơ bản thỏa mãn điều kiện biên hình học tại hai đầu dầm với hệ đa thức trực giao Gram-Schmidt hoặc đa thức Chebyshev cải biên.
  2. Khử hiện tượng khóa trượt trong mô hình Phần tử hữu hạn (FEM): Xây dựng phần tử dầm 2 nút với các bậc tự do chuyển vị được nội suy bằng hàm Hermite bậc 3 cho thành phần uốn $w_b$ (đảm bảo tính liên tục $C^1$) và hàm Lagrange tuyến tính cho thành phần trượt $w_s$ và chuyển vị dọc trục $u_0$ (đảm bảo tính liên tục $C^0$).
  3. Tam giác đạc phương pháp (Methodological Triangulation): Kết quả tính toán từ hàm dạng Ritz lai ghép được đối soát chéo độc lập với:
    • Nghiệm giải tích chính xác Navier (đối với biên S-S).
    • Mô hình Phần tử hữu hạn (FEM) tự lập trình.
    • Các công bố quốc tế chuẩn mực (Reddy, Simsek, Zenkour, Sayyad & Ghugal).

Data và phân tích

  • Cơ sở dữ liệu vật liệu: Khảo sát các cặp vật liệu gốm-kim loại với cơ tính phụ thuộc nhiệt độ (Temperature-Dependent - TD) được tính toán theo mô hình Touloukian (1967): $$P(T) = P_0 (P_{-1}T^{-1} + 1 + P_1 T + P_2 T^2 + P_3 T^3)$$ Bao gồm: Gốm Silicon Nitride ($Si_3N_4$), Alumina ($Al_2O_3$), Zirconia ($ZrO_2$) liên kết với Kim loại Thép không gỉ (SUS304), Nhôm ($Al$).
  • Phần mềm và thuật toán: Phát triển mã nguồn chuyên dụng trên nền tảng MATLAB (do nhóm GACES và các cộng sự hỗ trợ tối ưu hóa), giải bài toán trị riêng tổng quát: $$([K] - \omega^2 [M]){\Delta} = {0} \quad \text{(Dao động tự do)}$$ $$([K] - N_{cr} [K_G]){\Delta} = {0} \quad \text{(Mất ổn định uốn/nhiệt)}$$
  • Kiểm tra độ hội tụ (Convergence & Robustness Checks): Khảo sát sự hội tụ của tần số dao động riêng không thứ nguyên $\bar{\omega} = \omega \frac{L^2}{h}\sqrt{\frac{\rho_m}{E_m}}$ và tải trọng mất ổn định $N_{cr} = \frac{P_{cr} L^2}{E_m I}$ với số lượng số hạng hàm Ritz $N$ tăng dần từ 6 đến 16. Kết quả hội tụ tuyệt đối đạt được ở mức sai số $\le 0.05%$ khi $N \ge 10$.

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

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

  1. Hiệu ứng phi tuyến của Gradient vật liệu ($p$) lên độ cứng tổng thể: Khi chỉ số lũy thừa $p$ tăng từ 0 (gốm thuần túy) lên 10 (giàu kim loại), tần số dao động tự nhiên không thứ nguyên $\bar{\omega}$ và tải trọng mất ổn định tới hạn $N_{cr}$ suy giảm phi tuyến mạnh mẽ. Đối với dầm FGM đơn tầng ($Si_3N_4/SUS304$, $L/h = 20$), tần số cơ bản giảm hơn $42.6%$ khi $p$ tăng từ 0 lên 5 do tính mô-đun đàn hồi của gốm ($E_c \approx 348.43 \text{ GPa}$) vượt trội so với kim loại ($E_m \approx 201.04 \text{ GPa}$).
  2. Vai trò sống còn của biến dạng pháp ngang $\varepsilon_z$ trong dầm dày: Ở dầm dày ($L/h = 5$), mô hình Quasi-3D mới của luận án dự báo độ võng tĩnh không thứ nguyên $\bar{w}$ lớn hơn từ $12.4%$ đến $18.2%$ so với lý thuyết HSDT truyền thống và lớn hơn $25.7%$ so với lý thuyết cổ điển CBT. Điều này chứng minh việc bỏ qua hiệu ứng co giãn chiều dày ở dầm sandwich lõi mềm dẫn đến đánh giá non (thiếu an toàn) biến dạng công trình.
  3. Sự vượt trội của mô hình truyền nhiệt Fourier (NLTR): Nhiệt độ mất ổn định tới hạn $\bar{T}_{cr}$ tính toán theo quy luật truyền nhiệt phi tuyến Fourier thấp hơn đáng kể (từ $8.5%$ đến $15.3%$) so với giả định phân bố nhiệt tuyến tính (LTR) hoặc đều (UTR). Điều này khẳng định giả định gradient nhiệt tuyến tính kinh điển thường đánh giá quá cao khả năng chịu nhiệt của dầm FGM, tiềm ẩn nguy cơ phá hủy nhiệt cục bộ.
  4. Quy luật đối nghịch giữa hiệu ứng phi cục bộ nano và hiệu ứng couple stress micro:
    • Ở cấp độ nano, tham số phi cục bộ $\mu$ gia tăng làm mềm hóa kết cấu (Stiffness-softening effect), dẫn đến tần số dao động riêng $\bar{\omega}$ giảm rõ rệt.
    • Ở cấp độ micro, thông số chiều dài đặc trưng vật liệu $l/h$ (MCST) gia tăng lại làm cứng hóa kết cấu (Stiffness-hardening effect), làm tăng vọt tần số dao động và đẩy ngưỡng nhiệt độ mất ổn định $\bar{T}_{cr}$ lên cao gấp 1.5 - 2.2 lần so với lý thuyết đàn hồi cổ điển.
  5. Hiệu năng vượt trội của cấu hình Sandwich Type B (Lõi đẳng tính - Bề mặt FGM): Dầm sandwich Type B với tỷ lệ bề dày (1-2-1) thể hiện khả năng phân tán ứng suất cắt tiếp $\tau_{xz}$ trơn tru qua mặt phân giới, triệt tiêu hiện tượng tập trung ứng suất vốn xuất hiện gay gắt ở cấu hình Type C.
Cấu hình dầm / Tỷ số CBT (Cổ điển) FSDT (Timoshenko) HSDT Reddy Luận án (Unified HSBT / Quasi-3D)
Dầm dày ($L/h = 5, p = 1$) - Tần số $\bar{\omega}$ 5.2341 4.8762 4.8821 4.8895
Dầm dày ($L/h = 5$) - Tải trọng $N_{cr}$ 24.672 19.845 19.921 19.968
Dầm mảnh ($L/h = 100$) - Tần số $\bar{\omega}$ 9.8690 9.8652 9.8654 9.8655
Ứng suất trượt ngang $\tau_{xz}(\pm h/2)$ $\neq 0$ (Vi phạm) Hằng số (Vi phạm) $= 0$ (Thỏa mãn) $= 0$ (Thỏa mãn chính xác)
Hiệu ứng biến dạng ngang $\varepsilon_z$ Bỏ qua ($\varepsilon_z=0$) Bỏ qua ($\varepsilon_z=0$) Bỏ qua ($\varepsilon_z=0$) Thu giữ đầy đủ ($\varepsilon_z \neq 0$)

Implications đa chiều

  • Về mặt học thuật: Cung cấp hệ phương trình tường minh và thư viện hàm dạng trượt mới, đặt nền móng chuẩn hóa cho các nghiên cứu tiếp nối về tấm và vỏ composite thông minh tích hợp vật liệu áp điện hoặc hợp kim nhớ hình (SMA).
  • Về phương pháp luận: Thiết lập quy trình Ritz lai ghép và FEM 2 nút phi khóa trượt, ứng dụng hiệu quả cho phân tích kết cấu vi mô phức tạp mà không cần dùng đến các phần tử 3D Solid tốn kém tài nguyên tính toán.
  • Về ứng dụng thực tiễn: Cung cấp công cụ thiết kế chuẩn xác tấm chắn nhiệt cho phương tiện hàng không vũ trụ siêu thanh, cánh turbine phản lực khí, và các cảm biến/bộ kích hoạt vi cơ điện tử (MEMS Resonators), giúp ngăn ngừa hiện tượng mất ổn định nhiệt đàn hồi (thermal aeroelastic flutter).

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (Boundary Conditions & Limitations):

  1. Mô hình hóa hình học: Nghiên cứu tập trung vào dầm 1D (tiết diện chữ nhật không đổi), chưa mở rộng trực tiếp cho dầm có tiết diện biến đổi (tapered beams) hoặc kết cấu tấm/vỏ thoải 3D chịu trường xoắn phức tạp.
  2. Giả định liên kết mặt phân giới: Giả thiết các lớp vật liệu FGM và lõi sandwich liên kết dính bám hoàn hảo (perfect bonding), chưa xét đến suy giảm liên kết mềm (imperfect interface) hoặc sự phát triển của vết nứt tế vi theo thời gian.
  3. Mô hình vật liệu tất định: Các thông số cơ tính gốm và kim loại được giả định là biến thiên trơn tất định, chưa tích hợp yếu tố bất định ngẫu nhiên (stochastic/probabilistic analysis) do sai số chế tạo bột thiêu kết (powder metallurgy process).
  4. Phạm vi phi tuyến: Chủ yếu tập trung vào ứng xử tĩnh, mất ổn định và dao động tuyến tính tự do dưới gradient nhiệt-ẩm; chưa khảo sát ứng xử động lực phi tuyến hình học lớn (von Kármán large deflection) dưới tải trọng va đập cưỡng bức.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng lý thuyết: Tích hợp lý thuyết biến dạng góc trượt bậc cao vào bài toán phi tuyến hình học lớn và động lực học phi tuyến của tấm/vỏ FGM có gân tăng cứng.
  • Phát triển vật liệu mới: Khảo sát dầm FGM gia cường ống nano carbon (CNTs) hoặc hạt nano Graphene (GPLs) trong môi trường thủy nhiệt.
  • Nâng cấp phương pháp tính: Áp dụng phương pháp đẳng hình học (Isogeometric Analysis - IGA) dựa trên hàm B-splines/NURBS để mô phỏng chính xác hình học và đảm bảo tính liên tục trơn bậc cao $C^{k-1}$.
  • Nghiên cứu thực nghiệm: Chế tạo mẫu vi dầm FGM bằng công nghệ in 3D laser kim loại-gốm và đo kiểm dao động bằng hệ thống quét Laser Doppler Vibrometer để hiệu chuẩn chính xác tham số chiều dài vật liệu $l$ và hệ số $\mu$.

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

  • Tác động học thuật (Academic Impact): Luận án đã công bố nhiều công trình trên các tạp chí quốc tế ISI/Scopus uy tín chuyên ngành cơ học tính toán (như Composite Structures, European Journal of Mechanics - A/Solids...). Khung lý thuyết và mã tính MATLAB của luận án tạo tiền đề cho các nhóm nghiên cứu về cơ học vật liệu tiên tiến tại Việt Nam (GACES, HUTECH, TDTU, HCMUTE) trích dẫn và phát triển chuyên sâu.
  • Chuyển đổi công nghiệp (Industry Transformation): Định hình nguyên lý thiết kế tối ưu cấu trúc sandwich FGM cho ngành hàng không, công nghiệp đóng tàu quân sự và chế tạo thiết bị năng lượng nhiệt độ cao, giúp giảm khối lượng cấu kiện từ $15-30%$ trong khi vẫn duy trì khả năng chịu tải cơ - nhiệt vượt trội.
  • Chính sách và Tiêu chuẩn hóa (Policy & Standards): Đóng góp cơ sở dữ liệu tính toán tin cậy làm tài liệu tham chuẩn cho việc xây dựng các tiêu chuẩn kỹ thuật thiết kế kết cấu vật liệu tiên tiến chịu nhiệt tại Việt Nam trong kỷ nguyên công nghiệp 4.0.

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

  • Nghiên cứu sinh & Học viên Cao học (Doctoral/Master Researchers): Tiếp cận một khung lý thuyết mạch lạc từ vi phân đàn hồi 2D đến dầm Quasi-3D, nắm vững kỹ thuật lập trình Ritz lai ghép và phần tử hữu hạn tinh chế.
  • Các Giáo sư & Chuyên gia Cơ học tính toán: Sử dụng các kết quả kiểm chuẩn số (benchmark data) trong luận án để kiểm tra tính đúng đắn của các phương pháp số mới (S-FEM, Meshless, IGA).
  • Kỹ sư R&D Hàng không - Vũ trụ & Ô tô: Vận dụng trực tiếp các biểu đồ phân tích mất ổn định nhiệt để lựa chọn chỉ số gradient $p$ và tỷ lệ bề dày tối ưu cho các tấm chắn nhiệt và ống xả phản lực.
  • Nhà phát triển phần mềm CAE: Thu nhận thuật toán phần tử dầm bậc cao không khóa trượt để tích hợp vào các module phân tích vật liệu FGM trong các phần mềm tính toán kết cấu thương mại.

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 thiết lập mô hình lý thuyết dầm biến dạng trượt bậc cao thống nhất (Unified HSBT) và dầm Quasi-3D 3 biến số với hàm dạng hyperbolic-sine mới: $$f(z) = \sinh\left(\frac{rz}{h}\right) - \frac{8rz^3}{3h^3(r^2 - 4)}$$ Công trình này mở rộng trực tiếp lý thuyết biến dạng trượt bậc ba của Reddy (1984) và lý thuyết dầm Quasi-3D của Zenkour (2006). Điểm đột phá là loại bỏ hoàn toàn hệ số hiệu chỉnh trượt $K_s$, mô tả ứng suất trượt ngang $\tau_{xz}$ thỏa mãn chính xác điều kiện biên tự do trên hai bề mặt, đồng thời thu giữ được hiệu ứng biến dạng pháp ngang $\varepsilon_z \neq 0$ qua một số lượng biến số trường chuyển vị tối thiểu (3 biến so với 5 biến ở các mô hình Quasi-3D thông thường).

2. Đổi mới phương pháp luận của luận án thể hiện thế nào so với các nghiên cứu quốc tế trước đây?

So với các nghiên cứu của Simsek (2010) (dùng Navier cho dầm FGM đơn giản) và Pradhan & Murmu (2009) (dùng DQM nhưng gặp khó khăn về tính ổn định lưới), luận án đổi mới toàn diện ở hai khía cạnh:

  1. Xây dựng hệ hàm dạng Ritz lai ghép mới (Novel Hybrid Shape Functions) kết hợp hàm biên động học với chuỗi đa thức trực giao, cho phép giải quyết chính xác bài toán dao động và mất ổn định trên toàn bộ các tổ hợp điều kiện biên (S-S, C-C, C-F, C-H, H-H, C-S) mà Navier không giải được.
  2. Xây dựng mô hình FEM 2 nút tinh chế phân rã trường uốn - trượt, sử dụng phép nội suy hỗn hợp Hermite ($C^1$) - Lagrange ($C^0$), triệt tiêu hoàn toàn hiện tượng khóa trượt (shear locking) ở tỷ số dầm siêu mỏng ($L/h = 100$).

3. Phát hiện bất ngờ nhất về mặt số liệu trong nghiên cứu là gì?

Phát hiện bất ngờ và có ý nghĩa cảnh báo kỹ thuật sâu sắc nhất là: Mô hình truyền nhiệt phi tuyến Fourier (NLTR) chỉ ra rằng nhiệt độ mất ổn định tới hạn $\bar{T}_{cr}$ của dầm FGM thấp hơn tới $15.3%$ so với mô hình nhiệt tuyến tính (LTR). Trước đây, nhiều nghiên cứu quốc tế đơn giản hóa bài toán bằng cách giả định nhiệt độ phân bố tuyến tính qua chiều dày. Kết quả của luận án chứng minh rằng sự biến thiên phi tuyến của hệ số dẫn nhiệt $k(z)$ làm tích tụ mật độ nhiệt lớn hơn ở lớp kim loại dưới đáy, đẩy kết cấu vào trạng thái mất ổn định nhiệt sớm hơn nhiều so với dự báo của các mô hình kinh điển.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết:

  • Toàn bộ ma trận độ cứng $[K]$, ma trận khối lượng $[M]$, ma trận độ cứng hình học $[K_G]$ và vector tải nhiệt/ẩm dưới dạng tường minh tích phân qua chiều dày các hệ số $A_{ij}, B_{ij}, D_{ij}, E_{ij}, F_{ij}, H_{ij}$.
  • Bảng thông số vật liệu phụ thuộc nhiệt độ ($Si_3N_4, SUS304, Al, Al_2O_3$) với các hệ số đa thức $P_{-1}, P_0, P_1, P_2, P_3$ chuẩn hóa.
  • Thuật toán số giải bài toán trị riêng tổng quát trong môi trường MATLAB và bảng kiểm chứng số (Benchmark tables) với độ chính xác 4 chữ số thập phân cho từng mode dao động và tải trọng tới hạn.

5. Chương trình nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?

  1. Năm 1-3: Mở rộng lý thuyết Unified HSBT và Quasi-3D cho bài toán phân tích phi tuyến hình học lớn (Geometric Nonlinearity), phân tích sau mất ổn định (Post-buckling) và đáp ứng động lực va đập của dầm/tấm sandwich FGM.
  2. Năm 4-6: Tích hợp vật liệu thông minh (Graphene Platelets - GPLs, vật liệu áp điện Piezoelectric, Magnetostrictive) và mô hình hóa liên kết suy giảm liên lớp (Imperfect Interfacial Bonding).
  3. Năm 7-10: Phát triển công nghệ tính toán đẳng hình học (IGA) 3D, tích hợp thuật toán tối ưu hóa đa mục tiêu bằng trí tuệ nhân tạo (Machine Learning/AI surrogate modeling) để thiết kế tối ưu hóa vi cấu trúc FGM cho ngành công nghiệp vũ trụ và thiết bị y sinh.

Kết luận

Luận án Tiến sĩ của tác giả Nguyễn Bá Duy là một công trình nghiên cứu xuất sắc, hoàn chỉnh và có tính học thuật chuyên sâu cao trong lĩnh vực Cơ học vật liệu và kết cấu composite biến tính. Các đóng góp cốt lõi được đúc kết qua các điểm then chốt:

  1. Phát triển thành công hệ lý thuyết dầm biến dạng trượt bậc cao thống nhất (Unified HSBT) và lý thuyết Quasi-3D 3 biến số mới, sử dụng hàm phân bố trượt dạng hyperbolic-sine nguyên bản, triệt tiêu sự phụ thuộc vào hệ số hiệu chỉnh trượt và thu giữ chính xác hiệu ứng biến dạng pháp ngang $\varepsilon_z$.
  2. Mô hình hóa toàn diện tương tác liên hợp cơ - thủy - nhiệt phi tuyến, khẳng định tính tất yếu của việc áp dụng định luật truyền nhiệt Fourier và khuếch tán ẩm phi tuyến trong tính toán an toàn nhiệt đàn hồi cho dầm FGM sandwich.
  3. Tiên phong tích hợp đồng thời hai lý thuyết cơ học kích thước nhỏ (MCST cho cấp micro và Nonlocal Eringen cho cấp nano), luận giải sâu sắc cơ chế đối nghịch giữa hiệu ứng "cứng hóa" và "mềm hóa" kết cấu vi mô dưới các điều kiện biên đa dạng.
  4. Đột phá về phương pháp giải tích số với hệ hàm dạng Ritz lai ghép mới và mô hình FEM 2 nút phi khóa trượt, cung cấp công cụ tính toán hội tụ nhanh, ổn định và chính xác vượt trội cho các kỹ sư và nhà nghiên cứu.
  5. Mở ra các hướng nghiên cứu liên ngành mới về tối ưu hóa kết cấu composite thông minh, tạo nền tảng lý thuyết vững chắc cho các ứng dụng công nghệ cao trong ngành hàng không - vũ trụ và công nghệ vi hệ thống MEMS/NEMS toàn cầu.