Tổng quan về luận án

Nghiên cứu về cơ học kết cấu tiên tiến đóng vai trò then chốt trong sự phát triển của các ngành công nghiệp hàng không vũ trụ, năng lượng nguyên tử, kỹ thuật hàng hải và hệ thống vi cơ điện tử (MEMS/NEMS). Trong môi trường vận hành khắc nghiệt, các cấu kiện composite phân lớp truyền thống thường bộc lộ nhược điểm cố hữu là sự gián đoạn cơ tính tại bề mặt tiếp xúc giữa các lớp, dẫn đến hiện tượng tập trung ứng suất cục bộ cao và gây ra phá hủy tách lớp (delamination). Để khắc phục hạn chế này, vật liệu cơ tính biến thiên (Functionally Graded Materials - FGMs) được nghiên cứu ứng dụng nhằm tạo ra sự biến đổi liên tục, trơn tru của các đặc tính nhiệt cơ học theo chiều không gian. Luận án tiến sĩ kỹ thuật "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" của NCS. Nguyễn Bá Duy (2019), dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Trung Kiên và TS. Võ Phương Thức tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh, là công trình tiên phong giải quyết bài toán cơ học đa trường phức tạp này.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở việc thiếu vắng các mô hình lý thuyết bậc cao thống nhất, vừa phản ánh chính xác phân bố phi tuyến của ứng suất tiếp và hiệu ứng nở hông (thickness stretching), vừa tích hợp đồng thời trường biến dạng đa vật lý (cơ học, nhiệt độ phi tuyến Fourier, độ ẩm không đồng nhất) và hiệu ứng quy mô kích thước (size-dependent effects) trên dầm vi mô/nano với các điều kiện biên tùy ý. Các câu hỏi và giả thuyết nghiên cứu chính bao gồm:

  1. RQ1: Làm thế nào để xây dựng một lý thuyết biến dạng trượt bậc cao tổng quát (Unified HSBT) và lý thuyết chuẩn 3D (Quasi-3D) tối ưu hóa số lượng biến số trường chuyển vị mà không cần hệ số hiệu chỉnh trượt (shear correction factor)?
  2. RQ2: Mức độ suy giảm độ cứng và dịch chuyển tần số dao động riêng, tải trọng mất ổn định (buckling) của dầm sandwich FGM dưới tác động đồng thời của các quy luật phân bố nhiệt phi tuyến (NLTR) và độ ẩm (NLMR) diễn ra như thế nào?
  3. RQ3: Hiệu ứng kích thước vi mô (thông qua thuyết ứng suất cặp hiệu chỉnh MCST) và kích thước nano (thông qua thuyết đàn hồi phi cục bộ Eringen) chi phối phản ứng nhiệt - cơ của dầm FGM ra sao dưới các điều kiện biên khác nhau?
  4. H1: Hàm dạng trượt hyperbol ngược mới $f(z) = \sinh^{-1}(rz/h) - \frac{8rz^3}{3h^3(r^2+4)}$ tích hợp tham số hiệu chỉnh $r$ sẽ mô tả chính xác hơn trường ứng suất tiếp không ứng suất mặt tự do so với các mô hình kinh điển của Reddy hay Touratier.
  5. H2: Sự ghép nối giữa phương pháp Ritz sử dụng hàm dạng lai mới (hybrid shape functions) và nguyên lý biến phân Hamilton cho phép xác định nghiệm giải tích - số chính xác cao, loại bỏ triệt để hiện tượng khóa trượt (shear locking) trong phương pháp phần tử hữu hạn (FEM).

Khung lý thuyết của công trình được định vị trên nền tảng cơ học môi trường liên tục, cơ học vật liệu FGM (gốm kim loại $Si_3N_4/SUS304$, $Al/Al_2O_3$), lý thuyết đàn hồi phi cục bộ Eringen và lý thuyết ứng suất cặp hiệu chỉnh (MCST). Phạm vi nghiên cứu bao quát từ dầm đơn lớp (Type A), dầm sandwich lõi đồng tính - mặt FGM (Type B), dầm sandwich mặt đồng tính - lõi FGM (Type C) đến dầm biến thiên hai chiều (Type D) với tỷ lệ nhịp - chiều cao ($L/h$) từ 5 đến 100, mang lại giá trị định lượng vượt bậc cho phân tích độ bền và ổn định công trình.

Literature Review và Positioning

Tổng quan y văn cơ học kết cấu composite FGM ghi nhận sự phát triển qua nhiều giai đoạn lý thuyết. Khởi đầu từ lý thuyết dầm cổ điển Euler-Bernoulli (Classical Beam Theory - CBT), giả thiết mặt cắt phẳng vuông góc với trục trung hòa được duy trì sau biến dạng đã triệt tiêu hoàn toàn biến dạng trượt ngang ($\gamma_{xz} = 0$), khiến CBT chỉ áp dụng chính xác cho dầm mảnh ($L/h \ge 20$) như đã được chỉ ra bởi Ghugal và Shimpi (2001). Nhằm khắc phục hạn chế này, lý thuyết biến dạng trượt bậc nhất (First-order Shear Deformation Theory - FSDT) của Timoshenko đã tính đến trượt ngang nhưng lại giả định ứng suất tiếp phân bố đều qua chiều dày, bắt buộc phải sử dụng hệ số hiệu chỉnh trượt phụ thuộc vào hình học tiết diện.

Các lý thuyết biến dạng trượt bậc cao (Higher-order Shear Deformation Beam Theories - HSDTs) sau đó được phát triển mạnh mẽ nhằm thỏa mãn điều kiện ứng suất tiếp bằng 0 tại hai bề mặt tự do trên và dưới của dầm ($z = \pm h/2$). Điển hình là mô hình đa thức bậc 3 của Reddy (1984), mô hình lượng giác của Touratier (1991), mô hình hàm mũ của Karama và cộng sự (2003), cùng các phát triển gần đây của Sayyad và Ghugal (2015), Mantari và cộng sự (2012). Tuy nhiên, các cuộc tranh luận học thuật nảy sinh xung quanh hai trường phái:

  • Trường phái hàm dạng đại số/lượng giác tiêu chuẩn: Ưu tiên tính đơn giản giải tích nhưng khó phản ánh chính xác sự thay đổi đột ngột gradient cơ tính tại vùng chuyển tiếp pha vật liệu.
  • Trường phái chuẩn 3D (Quasi-3D): Bổ sung biến dạng pháp ngang ($\varepsilon_{zz} \neq 0$) để nắm bắt ứng xử 3 chiều, song thường làm tăng đáng kể số lượng ẩn số chuyển vị, gây khó khăn cho việc thiết lập ma trận độ cứng phần tử hữu hạn.

Trong bối cảnh quốc tế, các nghiên cứu của Simsek (2010, 2016) đã ứng dụng phương pháp Ritz và DQM phân tích dao động tự do và mất ổn định của dầm FGM 2D, hay nghiên cứu của Pradhan và Murmu (2009) về dao động nhiệt cơ dầm FGM. Tuy nhiên, các công trình này chủ yếu dừng lại ở điều kiện biên đơn giản (Simply-Supported - S-S), bỏ qua trường nhiệt phi tuyến nghiệm từ phương trình truyền nhiệt Fourier 1D, hoặc chưa khảo sát toàn diện tương tác đồng thời giữa độ ẩm và nhiệt độ kết hợp hiệu ứng kích thước vi mô/nano.

Luận án của NCS. Nguyễn Bá Duy định vị sự đóng góp vượt trội bằng cách lấp đầy khoảng trống này: thiết lập mô hình HSDT và Quasi-3D 3 biến số tối ưu, xây dựng hàm trượt hyperbol ngược mới, phát triển hàm dạng lai cho phương pháp Ritz đa dạng điều kiện biên (C-C, C-F, C-S, C-H, S-S), và tích hợp tham số chiều dài vật liệu $l$ (MCST) cùng tham số phi cục bộ $e_0a$ (Eringen) trong môi trường cơ - thủy - nhiệt phức hợ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 và làm sâu sắc thêm các lý thuyết cơ học môi trường liên tục và cơ học kết cấu tiên tiến:

  1. Mở rộng lý thuyết biến dạng trượt bậc cao (HSDTs): Thiết lập công thức biến phân tổng quát cho trường chuyển vị $u(x,z,t) = u_0(x,t) - zw_{0,x}(x,t) + f(z)\phi(x,t)$ và mô hình phân rã chuyển vị uốn - trượt $u(x,z,t) = u_0(x,t) - zw_{b,x}(x,t) - [f(z) - z]w_{s,x}(x,t)$, giúp giảm bậc tự do nhưng vẫn đảm bảo tính tương thích biến dạng cao cấp.
  2. Đột phá với hàm dạng trượt hyperbol ngược mới: Đề xuất dạng giải tích hàm trượt $f(z) = \sinh^{-1}(rz/h) - \frac{8rz^3}{3h^3(r^2+4)}$ với đạo hàm $g(z) = \frac{r}{\sqrt{r^2z^2+h^2}} - \frac{8rz^2}{h^3(r^2+4)}$. Hàm số này kết hợp hài hòa giữa đa thức bậc ba mô tả tính đồng nhất tổng thể và thành phần hyperbol ngược phản ánh sự thay đổi gradient cơ tính của FGM thông qua tham số hiệu chỉnh $r$.
  3. Mở rộng lý thuyết chuẩn 3D (Quasi-3D): Bổ sung hàm $g(z)$ vào chuyển vị đứng $w(x,z,t) = w_0(x,t) + g(z)w_z(x,t)$, cho phép tính toán trực tiếp ứng suất pháp ngang $\sigma_{zz}$, giải quyết trọn vẹn hiện tượng tập trung ứng suất 3 chiều trong dầm sandwich dày ($L/h = 5$).
  4. Phát triển lý thuyết kích thước phụ thuộc đa quy mô: Tích hợp thành công tensor độ cong đối xứng $\chi_{ij} = \frac{1}{2}(\theta_{i,j} + \theta_{j,i})$ trong MCST của Yang và cộng sự (2002) vào biểu thức thế năng biến dạng $U = \frac{1}{2}\int_V (\sigma_{ij}\varepsilon_{ij} + m_{ij}\chi_{ij})dV$, chứng minh rõ ràng vai trò của tham số tỉ lệ chiều dài vật liệu $l$ trong việc tăng cường độ cứng kết cấu vi mô.

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

Khung phân tích được xây dựng dựa trên sự liên kết đa tầng giữa ba lý thuyết nền tảng:

  • Lý thuyết đồng hóa vi mô FGM: Sử dụng luật phân bố lũy thừa (Power-law Type A, B, C), hàm Sigmoid và hàm mũ (Exponential) để nội suy modul đàn hồi $E(z)$, hệ số Poisson $\nu(z)$, khối lượng riêng $\rho(z)$, hệ số giãn nở nhiệt $\alpha(z)$, hệ số giãn nở ẩm $\beta(z)$ và hệ số dẫn nhiệt $k(z)$.
  • Lý thuyết cơ - thủy - nhiệt phi tuyến: Thiết lập phương trình phân bố nhiệt độ $T(z)$ giải từ phương trình dẫn nhiệt Fourier 1 chiều trạng thái dừng: $$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 trường độ ẩm biến thiên phi tuyến hình sin dạng $C(z) = C_b + \frac{C_t - C_b}{2}\left(1 - \cos\frac{\pi(2z+h)}{2h}\right)$.
  • Phương pháp biến phân năng lượng và chuyển vị giải tích: Áp dụng nguyên lý d'Alembert và nguyên lý Hamilton $\int_{t_1}^{t_2} (\delta U + \delta V - \delta K) dt = 0$ để thiết lập hệ phương trình vi phân chuyển động chính xác.
graph TD
    A["Mô hình cấu trúc FGM<br>(Gốm - Kim loại: Si3N4/SUS304, Al/Al2O3)"] --> B["Trường tải trọng Cơ - Thủy - Nhiệt<br>(UTR, LTR, Fourier NLTR, NLMR)"]
    A --> C["Hiệu ứng kích thước vi mô/nano<br>(MCST: Tham số l | Eringen: Tham số e0a)"]
    B --> D["Lý thuyết động học dầm bậc cao<br>(Unified HSBT & Quasi-3D | Hàm trượt sinh^-1)"]
    C --> D
    D --> E["Nguyên lý Hamilton & Biến phân năng lượng<br>delta(U + V - K) = 0"]
    E --> F["Phương pháp giải số/giải tích<br>(Ritz với hàm lai mới | FEM phần tử 2 nút)"]
    F --> G["Ứng xử kết cấu chuẩn xác<br>(Uốn tĩnh, Dao động tự do omega, Mất ổn định Ncr, Tcr)"]

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý thực chứng (positivism) với phương pháp tiếp cận diễn dịch toán lý nghiêm ngặt. Hệ thống mô phỏng số và giải tích được thiết kế đa cấp độ (Multi-level design), từ cấp độ vật liệu (xác định đặc tính phụ thuộc nhiệt độ $P(T) = P_0(P_{-1}T^{-1} + 1 + P_1T + P_2T^2 + P_3T^3)$), cấp độ vi cấu trúc (tính toán hiệu ứng kích thước hạt/chiều dài nội tại), đến cấp độ kết cấu vĩ mô (dầm đơn lớp, dầm sandwich 3 lớp với các tỷ lệ chiều dày 1-0-1, 1-1-1, 1-2-1, 2-1-2, 2-2-1).

Quy trình nghiên cứu rigorous

Quy trình tính toán thực hiện qua các bước chuẩn hóa:

  1. Thiết lập quan hệ ứng suất - biến dạng đa trường: $$\sigma_{xx} = Q_{11}(z) [\varepsilon_{xx} - \alpha(z,T)\Delta T(z) - \beta(z,C)\Delta C(z)]$$ $$\sigma_{xz} = Q_{55}(z) \gamma_{xz}, \quad \sigma_{zz} = Q_{22}(z) [\varepsilon_{zz} - \alpha(z,T)\Delta T(z) - \beta(z,C)\Delta C(z)]$$
  2. Xây dựng hệ hàm dạng lai mới (New Hybrid Shape Functions) cho phương pháp Ritz: Để khắc phục nhược điểm hội tụ chậm của đa thức thuần túy hoặc hàm lượng giác đơn lẻ khi giải quyết các điều kiện biên phức tạp (C-C, C-F, C-S, C-H), luận án đã phát triển bộ hàm trực giao kết hợp: $$X_m(x) = \sum_{j=1}^{N_r} c_j \Phi_j(x)$$ thỏa mãn tuyệt đối các điều kiện biên động học tại $x = 0$ và $x = L$.
  3. Kiểm tra độ hội tụ và độ tin cậy: Số lượng số hạng trong chuỗi Ritz ($N$) được khảo sát hội tụ từ $N=4$ đến $N=14$. Kết quả chứng minh từ $N=10$ đến $N=12$, nghiệm số đạt độ ổn định tuyệt đối với sai số nhỏ hơn $0.01%$.
  4. Mô hình hóa phần tử hữu hạn (FEM): Xây dựng phần tử dầm 2 nút bậc cao, mỗi nút chứa các bậc tự do chuyển vị $(u_0, w_b, w_{b,x}, w_s, w_{s,x})$, sử dụng hàm nội suy Hermite bậc ba cho chuyển vị uốn và trượt, triệt tiêu hoàn toàn hiện tượng khóa trượt số.

Data và phân tích

Dữ liệu vật liệu thực nghiệm chuẩn quốc tế được áp dụng trong toàn bộ các phân tích:

  • Cặp vật liệu gốm - kim loại $Si_3N_4/SUS304$: $E_c = 348.43\text{ GPa}$, $\rho_c = 2370\text{ kg/m}^3$, $\alpha_c = 5.8723 \times 10^{-6}\text{ K}^{-1}$, $k_c = 9.19\text{ W/mK}$; $E_m = 201.04\text{ GPa}$, $\rho_m = 8166\text{ kg/m}^3$, $\alpha_m = 12.330 \times 10^{-6}\text{ K}^{-1}$, $k_m = 12.04\text{ W/mK}$.
  • Cặp vật liệu $Al/Al_2O_3$: $E_m = 70\text{ GPa}$, $\rho_m = 2702\text{ kg/m}^3$; $E_c = 380\text{ GPa}$, $\rho_c = 3800\text{ kg/m}^3$.
  • Kỹ thuật tính toán số được lập trình toàn diện trên môi trường MATLAB, xử lý hệ phương trình đại số phi tuyến và bài toán giá trị riêng suy rộng $[K]{\Delta} - \omega^2[M]{\Delta} = 0$ hoặc $([K] - N_{cr}[K_g]){\Delta} = 0$.

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

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

  1. Độ chính xác vượt trội của hàm trượt hyperbol ngược mới: Kết quả đối sánh tại Bảng 3.4 và 3.5 cho thấy tần số tự nhiên không thứ nguyên $\bar{\omega} = \omega \frac{L^2}{h}\sqrt{\frac{\rho_0}{E_0}}$ và tải trọng tới hạn $N_{cr}$ của mô hình HSBT đề xuất tiệm cận nghiệm giải tích 2D Elasticity chính xác hơn so với mô hình Reddy và Touratier, với sai số tương đối giảm dưới $0.15%$ trên toàn dải chỉ số thể tích $p \in [0, 10]$.
  2. Ảnh hưởng suy giảm độ cứng nhiệt - ẩm phi tuyến: Dưới tác động tăng nhiệt độ phi tuyến Fourier (NLTR) và tăng độ ẩm hình sin (NLMR), độ cứng tổng thể của dầm dầm giảm nghiêm trọng. Ví dụ với dầm $Si_3N_4/SUS304$ ($p=1, L/h=20$), khi độ tăng nhiệt độ $\Delta T$ tăng từ $0\text{ K}$ lên $100\text{ K}$ kết hợp $\Delta C = 1%$, tần số dao động cơ bản giảm hơn $28.4%$, và nhiệt độ mất ổn định tới hạn $T_{cr}$ của phân bố Fourier thấp hơn đáng kể so với giả định tăng nhiệt tuyến tính (LTR) truyền thống.
  3. Quy luật chi phối của cấu trúc phân lớp Sandwich: Đối với dầm sandwich Type B (mặt FGM, lõi cứng isotropic), khả năng chịu lực mất ổn định $N_{cr}$ cao hơn từ $1.35$ đến $1.82$ lần so với dầm sandwich lõi mềm (soft core) cùng tỷ lệ thể tích. Phân bố ứng suất pháp $\bar{\sigma}{xx}$ và ứng suất tiếp $\bar{\sigma}{xz}$ qua chiều dày thể hiện tính liên tục hoàn hảo tại các mặt tiếp giáp, triệt tiêu hoàn toàn bước nhảy ứng suất (stress jump) vốn là nguyên nhân gây phá hủy liên kết.
  4. Hiện tượng "làm cứng" do hiệu ứng kích thước vi mô (MCST): Khi xem xét tham số tỉ lệ chiều dài vật liệu $l/h$ tăng từ $0$ lên $0.6$ (Hình 5.5), tần số dao động riêng không thứ nguyên tăng vọt lên tới $45 - 60%$ tùy thuộc vào điều kiện biên (đặc biệt nhạy cảm với biên ngàm C-C). Điều này chứng minh lý thuyết cổ điển đánh giá thấp đáng kể độ cứng của các cấu kiện dầm vi cơ điện tử.
  5. Hiệu ứng "làm mềm" của tham số phi cục bộ Eringen ($e_0a$): Ngược lại với MCST, khi tham số phi cục bộ $\mu = (e_0a)^2$ tăng từ $0$ đến $4\text{ nm}^2$ trên dầm nano (Bảng 5.4 - 5.6), tần số dao động tự do lại giảm rõ rệt, do các tương tác tầm xa giữa các nguyên tử làm suy giảm độ cứng liên kết vĩ mô.

Implications đa chiều

  • Về mặt học thuật & lý thuyết: Bổ sung vào hệ thống lý thuyết đàn hồi phi cổ điển một khung toán học thống nhất, chứng minh sự tương tác đối nghịch nhưng bổ trợ giữa hiệu ứng ứng suất cặp (làm cứng) và hiệu ứng phi cục bộ (làm mềm).
  • Về mặt phương pháp luận: Khẳng định tính ưu việt của hàm dạng lai trong phương pháp Ritz và quy trình phần tử hữu hạn không khóa trượt, có thể chuyển giao trực tiếp để phân tích tấm, vỏ composite FGM phức tạp.
  • Về mặt ứng dụng kỹ thuật: Cung cấp biểu đồ và bảng tra cứu thông số kỹ thuật chuẩn xác cho các kỹ sư hàng không vũ trụ trong thiết kế tấm chắn nhiệt cho tàu con thoi khi hồi quyển (chịu gradient nhiệt độ từ $300\text{ K}$ đến trên $1000\text{ K}$), cũng như thiết kế cánh dầm cảm biến vi cơ điện tử (MEMS resonators) đạt độ ổn định tần số cộng hưởng cao.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, công trình vẫn tồn tại một số giới hạn học thuật:

  • Giả thiết phân bố trường nhiệt độ và độ ẩm mới chỉ dừng lại ở bài toán 1 chiều (theo chiều dày $z$), chưa giải quyết trọn vẹn trường nhiệt 3 chiều có đối lưu và bức xạ phức tạp.
  • Chưa tích hợp biến dạng lớn phi tuyến hình học von Kármán vào bài toán dầm vi mô/nano chịu tải trọng động lực học va đập.
  • Kiểm chứng thực nghiệm (experimental validation) của vật liệu nano FGM còn hạn chế do điều kiện chế tạo mẫu nano đa pha có gradient liên tục đòi hỏi công nghệ luyện kim bột nano tối tân.

Chương trình nghiên cứu tương lai cần tập trung vào 4 hướng chiến lược:

  1. Phát triển mô hình tấm và vỏ sandwich FGM 2D/3D tích hợp áp điện (piezoelectric) và từ giảo (magnetostrictive) thông minh.
  2. Ứng dụng lý thuyết biến dạng gradient bậc cao (Strain Gradient Theory) kết hợp cơ học phá hủy liên tục để mô phỏng sự lan truyền vết nứt mỏi nhiệt.
  3. Thiết lập thuật toán tối ưu hóa cấu trúc FGM đa mục tiêu (vừa tối đa hóa $N_{cr}$, vừa tối thiểu hóa khối lượng và độ võng) bằng trí tuệ nhân tạo (AI/Machine Learning).
  4. Nghiên cứu thực nghiệm kiểm chứng trên vật liệu FGM in 3D (Additive Manufacturing) đa kim loại - ceramic.

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

Luận án của TS. Nguyễn Bá Duy đã tạo ra những ảnh hưởng học thuật và thực tiễn sâu sắc. Các bài báo công bố trên các tạp chí quốc tế ISI/Scopus uy tín (như Composite Structures, Aerospace Science and Technology, European Journal of Mechanics - A/Solids) đã thu hút nhiều trích dẫn khoa học, khẳng định vị thế của nhóm nghiên cứu GACES (Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh) trên trường quốc tế.

Về công nghiệp và xã hội, các mô hình tính toán giúp giảm thiểu chi phí thử nghiệm phá hủy vật liệu đắt tiền trong ngành hàng không và năng lượng. Việc kiểm soát chính xác ứng xử mất ổn định nhiệt cơ góp phần nâng cao hệ số an toàn cho các kết cấu tuabin khí thế hệ mới, lò phản ứng hạt nhân và các thiết bị MEMS y sinh, thúc đẩy sự phát triển bền vững của ngành công nghệ vật liệu tiên tiến tại Việt Nam và trên thế giới.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận mã nguồn thuật toán giải tích, phương pháp thiết lập hàm Ritz lai và quy trình xử lý bài toán biến phân cơ học đa trường.
  • Giáo sư & Chuyên gia nghiên cứu cao cấp: Khai thác các đóng góp lý thuyết về hàm trượt hyperbol ngược mới và mô hình Quasi-3D để mở rộng cho các bài toán kết cấu phức tạp hơn.
  • Kỹ sư R&D công nghiệp: Ứng dụng các công thức giải tích đóng và dữ liệu thiết kế để chế tạo các cấu kiện dầm/cánh composite FGM chịu nhiệt độ và độ ẩm cao.
  • Nhà hoạch định chính sách công nghệ: Có thêm cơ sở khoa học định lượng để đầu tư chiến lược vào các phòng thí nghiệm chế tạo vật liệu tiên tiến và công nghệ vi cơ điện tử 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 đề xuất hàm dạng trượt hyperbol ngược mới $f(z) = \sinh^{-1}(rz/h) - \frac{8rz^3}{3h^3(r^2+4)}$ tích hợp trong mô hình dầm Quasi-3D 3 biến số. Công trình này mở rộng trực tiếp lý thuyết trượt bậc cao kinh điển của Reddy (1984) và Touratier (1991), giải quyết hoàn hảo bài toán biến dạng pháp ngang $\varepsilon_{zz} \neq 0$ mà không cần đến hệ số hiệu chỉnh trượt, đồng thời cho phép tinh chỉnh nghiệm chính xác thông qua tham số $r$.

2. Sự đổi mới về 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 Navier chỉ giải được biên gối tựa (S-S) của Zenkour (2006) hay phương pháp phần tử hữu hạn truyền thống dễ bị khóa trượt của Chakraborty và cộng sự (2003), luận án đã phát triển bộ hàm dạng lai trực giao mới áp dụng cho phương pháp Ritz. Sự cải tiến này cho phép giải quyết chính xác tuyệt đối mọi tổ hợp điều kiện biên (C-C, C-F, C-S, C-H, S-S) với tốc độ hội tụ siêu việt (chỉ cần $N=10$ số hạng) và chi phí tính toán chỉ bằng một phần nhỏ so với FEM thương mại (như ANSYS/ABAQUS).

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Phát hiện bất ngờ nhất là sự đối kháng cơ học giữa hai trường tải trọng: trong khi gradient nhiệt độ phi tuyến (NLTR) làm mềm kết cấu và giảm mạnh tần số dao động tự do, thì hiệu ứng kích thước vi mô (MCST) với tham số $l/h = 0.6$ lại có khả năng "bù trừ độ cứng", làm tăng tần số dao động lên tới hơn $50%$. Điều này chứng minh rằng ở quy mô vi cơ điện tử, hiệu ứng kích thước có thể áp đảo sự suy thoái nhiệt cơ thông thường.

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

Có. Luận án cung cấp chi tiết toàn bộ các phương trình vi phân chủ đạo, các biểu thức giải tích ma trận khối lượng $[M]$, ma trận độ cứng tĩnh $[K]$, ma trận độ cứng hình học $[K_g]$, bảng hệ số vật liệu phụ thuộc nhiệt độ cho $Si_3N_4/SUS304$ và $Al/Al_2O_3$, cùng cấu trúc thuật toán giải trên MATLAB, cho phép tái lập kết quả với độ lệch số học $\le 0.01%$.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm xác định rõ: Giai đoạn 1-3 năm mở rộng sang tấm/vỏ FGM 2 chiều chịu tải va đập động lực học phi tuyến; Giai đoạn 4-6 năm tích hợp vật liệu nano lai (Graphene/CNT reinforced FGM) và cấu trúc xốp (porous FGMs); Giai đoạn 7-10 năm hoàn thiện công nghệ chế tạo thực nghiệm in 3D FGM và thử nghiệm độ bền mỏi trong buồng mô phỏng không gian nhiệt - chân không cực hạn.

Kết luận

  1. Xây dựng thành công hệ thống lý thuyết dầm biến dạng trượt bậc cao tổng quát (Unified HSBT) và dầm chuẩn 3D (Quasi-3D) với hàm dạng trượt hyperbol ngược mới có tham số điều chỉnh $r$, triệt tiêu hoàn toàn sự phụ thuộc vào hệ số hiệu chỉnh trượt.
  2. Giải quyết trọn vẹn bài toán tương tác cơ - thủy - nhiệt phức hợp, chỉ ra sự suy giảm độ cứng phi tuyến rõ rệt của dầm FGM dưới tác động đồng thời của trường nhiệt độ Fourier (NLTR) và độ ẩm hình sin (NLMR).
  3. Làm sáng tỏ quy luật quy mô kích thước đa tầng: chứng minh định lượng hiệu ứng tăng cứng của lý thuyết ứng suất cặp hiệu chỉnh (MCST) trên dầm vi mô và hiệu ứng làm mềm của lý thuyết phi cục bộ Eringen trên dầm nano dưới đa dạng điều kiện biên.
  4. Phát triển thành công phương pháp số hiệu năng cao kết hợp giữa hàm dạng lai trong phương pháp Ritz và mô hình phần tử hữu hạn 2 nút không khóa trượt.
  5. Cung cấp bộ cơ sở dữ liệu số chuẩn mực quốc tế về ứng suất, độ võng, tải trọng mất ổn định $N_{cr}$ và tần số tự nhiên $\bar{\omega}$ cho 4 cấu hình dầm FGM (Type A, B, C, D), phục vụ trực tiếp cho công tác tính toán, thiết kế và tối ưu hóa kết cấu tiên tiến trong kỹ thuật hiện đại.