Luận án tiến sĩ về phân tích phi tuyến dầm sandwich FGP - Hương Quý Trường

Tổng quan về luận án

Công nghệ vật liệu composite đa lớp và kết cấu sandwich tiên tiến đang đóng vai trò then chốt trong ngành công nghiệp hàng không vũ trụ, kỹ thuật hàng hải và xây dựng công trình hiện đại. Kết cấu sandwich truyền thống gồm ba lớp với lõi nhẹ nằm giữa hai lớp bề mặt mỏng có độ cứng lớn thường bộc lộ hạn chế nghiêm trọng về hiện tượng tách lớp (delamination) do sự mất liên tục về cơ tính tại mặt tiếp xúc giữa các pha vật liệu. Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 9520101) của nghiên cứu sinh Hương Quý Trường tại Trường Đại học Xây dựng Hà Nội (2023), dưới sự hướng dẫn khoa học của TS. Đặng Xuân Hùng và GS. TS. Trần Minh Tú, đã tạo bước đột phá khi đề xuất mô hình dầm sandwich ba lớp thế hệ mới: hai lớp bề mặt làm bằng vật liệu có cơ tính biến thiên (FGM) và lớp lõi làm bằng vật liệu FGM rỗng (FGP) được gia cường bởi các phiến nano graphene (GPL), đặt trên nền đàn hồi Pasternak chịu liên kết biên đa dạng.

Khoảng trống nghiên cứu (research gap) mà luận án giải quyết xuất phát từ việc các nghiên cứu trước đây (như Frostig & Baruch [40, 41], Vo & Nguyen [106, 151], Chen et al. [23], Wang et al. [156]) chủ yếu tập trung vào phân tích tuyến tính hoặc chỉ khảo sát riêng rẽ từng trạng thái tĩnh hoặc động của dầm sandwich FGM lõi đẳng hướng, hoặc dầm lõi xốp chưa xét đến sự kết hợp đồng thời của pha nano GPL và lớp bề mặt FGM trong một khung lý thuyết biến dạng cắt tổng quát. Cụ thể, chưa có công trình nào thiết lập một lời giải bán giải tích hoàn chỉnh khảo sát đồng thời quan hệ phi tuyến tải - độ võng, trường ứng suất thực tế qua bề dày, tần số dao động tự do phi tuyến và đáp ứng chuyển vị cưỡng bức phi tuyến theo thời gian trong cùng một hệ thống tính toán nhất quán.

Các câu hỏi và giả thuyết nghiên cứu chính được xác lập rõ ràng:

1. RQ1: Làm thế nào để xây dựng một mô hình toán học tổng quát duy nhất có khả năng tích hợp linh hoạt các lý thuyết dầm từ bậc nhất (FSDBT) đến các lý thuyết biến dạng cắt bậc cao (PSDBT, TSDBT, HSDBT, ESDBT) cho dầm sandwich phức hợp FGM-FGP-GPL?
2. RQ2: Ảnh hưởng tương hỗ giữa quy luật phân bố lỗ rỗng (PD1, PD2, UPD), hệ số rỗng ($e_0 = 0,1 \div 0,9$), và quy luật phân tán GPL (Dạng A, Dạng B, Dạng C) với hàm lượng $W_{GPL} = 0% \div 2%$ tác động như thế nào đến độ cứng uốn và ứng xử cơ học phi tuyến của dầm?
3. RQ3: Sự tham gia của hiệu ứng phi tuyến hình học von Kármán và tương tác nền đàn hồi hai tham số Pasternak làm thay đổi các đặc trưng dao động tự do phi tuyến ($\Omega_{NL}/\Omega_L$) và quỹ đạo pha không gian trạng thái ra sao?
4. H1 (Giả thuyết 1): Việc gia cường một lượng rất nhỏ phiến nano GPL (khoảng $1,0% \div 2,0%$ khối lượng) vào lớp lõi FGP phân bố theo dạng đối xứng (Dạng A) sẽ triệt tiêu hoàn toàn sự suy giảm độ cứng do các lỗ rỗng gây ra, đồng thời giảm độ võng phi tuyến từ $30%$ đến $50%$.
5. H2 (Giả thuyết 2): Mô hình bán giải tích Pb-Ritz kết hợp thuật toán lặp Newton-Raphson và Runge-Kutta bậc 4 bảo đảm sự hội tụ đơn điệu nhanh và độ chính xác vượt trội so với các nghiệm giải tích Navier hay phần tử hữu hạn truyền thống trên mọi điều kiện biên từ ngàm (CC), khớp (SS), ngàm-khớp (CS) đến ngàm-tự do (CF).

Khung lý thuyết của luận án tích hợp chặt chẽ giữa mô hình vi cơ học Halpin-Tsai mở rộng, quy tắc hỗn hợp Voigt, mô hình vật liệu bọt xốp cấu trúc tế bào kín theo quy luật trường ngẫu nhiên Gauss (closed-cell Gaussian Random Field - GRF), nguyên lý biến phân Hamilton và hình học phi tuyến von Kármán. Phạm vi khảo sát trải rộng trên dải tham số kích thước dầm từ dày đến mỏng ($L/h = 5 \div 20$), các tỷ lệ phân chia chiều dày lớp $(1\text{-}1\text{-}1, 1\text{-}2\text{-}1, 1\text{-}4\text{-}1, 1\text{-}8\text{-}1)$, chỉ số thể tích FGM bề mặt ($p = 0 \div 5$), hệ số nền đàn hồi Winkler ($K_0 = 0 \div 100$) và hệ số cắt Pasternak ($J_0 = 0 \div 10$).

Literature Review và Positioning

Các nghiên cứu về kết cấu sandwich trong cơ học vật rắn biến dạng đã trải qua ba giai đoạn phát triển chính. Giai đoạn đầu tập trung vào kết cấu sandwich cổ điển với lớp bề mặt kim loại hoặc composite sợi và lớp lõi đẳng hướng, lõi mềm hoặc lõi bọt polymer (Frostig & Baruch [40, 41], Damanpack & Khalili [30], Salami et al. [128]). Nhược điểm cốt tử của dòng kết cấu này là sự bất liên tục về cơ tính tại mặt phân cách giữa các lớp, dẫn đến ứng suất tập trung cục bộ và nguy cơ tách lớp sớm khi làm việc dưới tải trọng chu kỳ hoặc gradient nhiệt độ cao.

Giai đoạn thứ hai đánh dấu sự xuất hiện của vật liệu có cơ tính biến thiên FGM (Zenkour [167], Vo & Nguyen [106, 151], Sayyad & Avhad [130], Simsek & Al-Shujairi [136]). Các nghiên cứu này đã giải quyết triệt để sự mất liên tục về vật liệu bằng cách cho tỷ phần pha gốm ($ZrO_2$, $Al_2O_3$) và kim loại ($Ti\text{-}6Al\text{-}4V$, $Al$) thay đổi liên tục theo quy luật hàm lũy thừa qua chiều cao tiết diện. Tuy nhiên, việc sử dụng lớp lõi FGM đặc thuần túy vẫn khiến trọng lượng bản thân kết cấu còn lớn, làm giảm hiệu suất chịu tải riêng trong các phương tiện bay.

Giai đoạn thứ ba hướng đến việc ứng dụng vật liệu FGM rỗng dạng bọt (FGP) kết hợp gia cường các cấu trúc nano (Kumar & Harsha [75], Wang et al. [155, 156], Chen et al. [23], Hung et al. [116]). Tuy nhiên, các tranh luận học thuật sâu sắc vẫn tồn tại ở hai khía cạnh chính:

1. Tranh luận về việc lựa chọn lý thuyết dầm: Lý thuyết dầm cổ điển Euler-Bernoulli (CBT) bỏ qua biến dạng trượt ngang nên cho kết quả thiếu chính xác đối với dầm dày hoặc dầm sandwich có lõi mềm. Lý thuyết dầm Timoshenko (FSDBT) kể đến trượt bậc nhất nhưng đòi hỏi hệ số hiệu chỉnh cắt $K_s$ (thường chọn $5/6$) vốn phụ thuộc phức tạp vào hình học và điều kiện biên. Các lý thuyết biến dạng cắt bậc cao (HSDT) như bậc ba của Reddy (PSDBT), hàm sin (TSDBT), hàm hyperbol (HSDBT) và hàm e-mũ của Karama (ESDBT) thỏa mãn điều kiện ứng suất tiếp triệt tiêu tại hai bề mặt tự do mà không cần hệ số hiệu chỉnh, nhưng mỗi tác giả trước đây lại xây dựng phương trình riêng rẽ, gây khó khăn cho việc đánh giá so sánh chéo.
2. Tranh luận về mô hình hóa cơ tính vi mô của lớp lõi rỗng: Nhiều tác giả sử dụng các giả định xốp đơn giản bỏ qua sự tương quan phi tuyến giữa khối lượng riêng và mô đun đàn hồi, dẫn đến sai lệch lớn khi hệ số rỗng $e_0 > 0,5$.

Luận án của NCS. Hương Quý Trường định vị một cách chính xác vào trung tâm của khoảng trống học thuật này. Bằng việc tích hợp mô hình closed-cell GRF của Roberts & Garboczi [123, 124] với hệ thức vi cơ học Halpin-Tsai và thiết lập một trường chuyển vị tổng quát đại diện cho 6 lý thuyết dầm khác nhau, nghiên cứu đã khắc phục toàn diện những hạn chế kể trên. Kết quả kiểm chứng được so sánh trực tiếp với hai nghiên cứu quốc tế điển hình: nghiên cứu của Zenkour [167] về dầm sandwich FGM trên nền đàn hồi và nghiên cứu của Chen et al. [23] về dao động phi tuyến dầm xốp FGM, cho thấy sự tương thích tuyệt đối ở các trường hợp suy biến đồng thời mở rộng năng lực phân tích sang dầm sandwich ba lớp FGM-FGP-GPL đầy đủ.

Đó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 phong phú thêm hệ thống lý thuyết tính toán kết cấu composite tiên tiến thông qua các đóng góp căn bản:

• Tổng quát hóa trường chuyển vị dầm: Thay vì xây dựng các hệ phương trình riêng biệt cho từng lý thuyết, luận án thiết lập trường chuyển vị dạng giải tích tổng quát: $$u(x,z,t) = u_0(x,t) - z\frac{\partial w_0(x,t)}{\partial x} + f(z)\theta(x,t)$$ $$w(x,z,t) = w_0(x,t)$$ trong đó hàm phân bố trượt $f(z)$ được tham số hóa linh hoạt cho phép chuyển đổi tức thời giữa FSDBT ($f(z)=z$), PSDBT ($f(z)=z - \frac{4}{3h^2}z^3$), TSDBT ($f(z)=\frac{h}{\pi}\sin(\frac{\pi z}{h})$), HSDBT ($f(z)=h\sinh(\frac{z}{h}) - z\cosh(\frac{1}{2})$) và ESDBT ($f(z)=z e^{-2(z/h)^2}$).
• Thiết lập mô hình cơ học phi tuyến von Kármán: Kể đến số hạng biến dạng phi tuyến $\frac{1}{2}(\frac{\partial w}{\partial x})^2$ trong trường biến dạng dọc trục: $$\varepsilon_x = \frac{\partial u_0}{\partial x} - z\frac{\partial^2 w_0}{\partial x^2} + f(z)\frac{\partial \theta}{\partial x} + \frac{1}{2}\left(\frac{\partial w_0}{\partial x}\right)^2$$ từ đó dẫn xuất hệ phương trình vi phân chuyển động phi tuyến thông qua việc giải phiếm hàm biến phân Hamilton: $$\delta \int_{t_1}^{t_2} \left[ T - (U_B + U_F - W) \right] dt = 0$$
• Cơ chế chuyển dịch mô thức (Paradigm Shift): Minh chứng bằng chứng lượng hóa rằng việc điều chỉnh quy luật biến thiên cơ tính không gian (spatial tailoring) thông qua phân bố GPL và bọt xốp có thể thay đổi hoàn toàn tính chất làm cứng (hardening behavior) phi tuyến của hệ thống mà không làm tăng thể tích kết cấu.

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

Khung phân tích của luận án là sự hợp nhất đa tầng của ba lý thuyết thành phần:

1. Lý thuyết pha FGM bề mặt: Tuân theo mô hình biến thiên hàm lũy thừa Voigt qua chiều cao $h_b$ và $h_t$: $$E_b(z) = E_m + (E_c - E_m)\left(\frac{z - z_1}{z_2 - z_1}\right)^p; \quad E_t(z) = E_m + (E_c - E_m)\left(\frac{z - z_4}{z_3 - z_4}\right)^p$$
2. Lý thuyết pha lõi FGP - Closed-cell GRF: Biểu diễn mối liên hệ giữa mô đun đàn hồi $E_{core}(z)$ và khối lượng riêng $\rho_{core}(z)$ thông qua số mũ phi tuyến $2,3$: $$\frac{E_{core}(z)}{E_1(z)} = \left(\frac{\rho_{core}(z)}{\rho_1(z)}\right)^{2,3} \left(\frac{1 + 0,121 \rho_{core}/\rho_1}{1,121}\right)$$ với ba quy luật phân bố lỗ rỗng: PD1 (đối xứng: $\cos(\frac{\pi z}{h_c})$), PD2 (bất đối xứng: $\cos(\frac{\pi z}{2h_c} + \frac{\pi}{4})$), và UPD (đều). Mối quan hệ bảo toàn khối lượng tổng thể dẫn xuất các hệ số phụ thuộc thực nghiệm chính xác, ví dụ khi $e_{01} = 0,5$, giá trị tương ứng của hệ số rỗng đều là $\alpha = 0,4638$ và $\alpha' = 0,6817$.
3. Lý thuyết nano composite GPL - Halpin-Tsai cải tiến: Xác định mô đun đàn hồi hiệu dụng $E_1(z)$ của pha nền kim loại gia cường GPL: $$E_1(z) = \frac{3}{8}\left[\frac{1 + \xi_L \eta_L V_{GPL}(z)}{1 - \eta_L V_{GPL}(z)}\right] E_m + \frac{5}{8}\left[\frac{1 + \xi_W \eta_W V_{GPL}(z)}{1 - \eta_W V_{GPL}(z)}\right] E_m$$ với các hệ số hình học nano $\xi_L = 2 L_{GPL}/t_{GPL}$ và $\xi_W = 2 w_{GPL}/t_{GPL}$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế trên lập trường triết học thực chứng (positivism) với phương pháp tiếp cận bán giải tích (semi-analytical approach) kết hợp mô phỏng số học chính xác cao. Cấu trúc thiết kế được phân tách thành hai cấp độ:

• Cấp độ vi mô (Micro-scale): Mô hình hóa cấu trúc tế bào xốp vi mô và sự định hướng phân tán của các phiến nano GPL ngẫu nhiên trong nền kim loại để xác định các đặc trưng cơ tính tương đương.
• Cấp độ vĩ mô (Macro-scale): Tích phân năng lượng biến dạng toàn phần của dầm sandwich ba lớp có kể đến năng lượng tương tác nền đàn hồi Pasternak ($U_F = \frac{1}{2}\int_0^L (K_W w^2 + K_P (\frac{\partial w}{\partial x})^2) dx$).

Quy trình nghiên cứu rigorous

Luận án triển khai phương pháp Pb-Ritz (Polynomial-boundary Ritz) bằng cách xấp xỉ các thành phần chuyển vị thông qua tích của chuỗi đa thức đại số bậc $m$ với hàm biên hình học $X(x)$: $$u_0(x,t) = X_u(x) \sum_{i=1}^m p_i(x) u_i(t); \quad w_0(x,t) = X_w(x) \sum_{i=1}^m p_i(x) w_i(t); \quad \theta(x,t) = X_\theta(x) \sum_{i=1}^m p_i(x) \theta_i(t)$$ trong đó hàm biên được định nghĩa chính xác theo điều kiện biên dầm: $X(x) = (x/L)^{p_0} (1 - x/L)^{p_1}$. Các số mũ $(p_0, p_1)$ nhận giá trị cụ thể: $(1, 1)$ cho dầm ngàm hai đầu (CC); $(0, 0)$ cho dầm tự do hai đầu; $(1, 0)$ cho dầm ngàm-tự do (CF); và $X_w(x) = (x/L)(1 - x/L)$ cho dầm khớp hai đầu (SS).

Quy trình khảo sát sự hội tụ được thực hiện nghiêm ngặt khi tăng dần số số hạng chuỗi đa thức từ $m = 4$ đến $m = 14$. Kết quả chứng minh rằng nghiệm số đạt trạng thái ổn định tuyệt đối với sai số tương đối $< 0,01%$ ngay từ bậc $m = 10 \div 12$. Độ tin cậy được kiểm chứng chéo thông qua việc tái lập chính xác các kết quả chuẩn quốc tế về dầm FGM của Zenkour [167] và dầm FGP của Chen et al. [23].

                +-------------------------------------------------------------+
                |             Mo hinh vat lieu dac thu da tang                |
                |  - Be mat: FGM (Voigt power-law index p)                    |
                |  - Loi: FGP (Closed-cell GRF: e0) + GPL (Halpin-Tsai: WGPL) |
                +-------------------------------------------------------------+
                                              |
                                              v
                +-------------------------------------------------------------+
                |        Truong chuyen vi va bien dang tong quat              |
                |  - Ly thuyet cat tong quat: CBT, FSDBT, PSDBT, TSDBT, ESDBT |
                |  - Bien dang phi tuyen hinh hoc von Karman                  |
                +-------------------------------------------------------------+
                                              |
                                              v
                +-------------------------------------------------------------+
                |          Thiet lap phuong trinh qua Bien phan Hamilton      |
                |               delta \int [T - (U_B + U_F - W)] dt = 0       |
                |          Phuong phap Pb-Ritz voi ham da thuc truc giao      |
                +-------------------------------------------------------------+
                                              |
                     +------------------------+------------------------+
                     |                                                 |
                     v                                                 v
   +------------------------------------+            +------------------------------------+
   |     Phan tich phi tuyen Tinh       |            |     Phan tich phi tuyen Dao dong   |
   | [K_L + K_NL(q)] * q = F            |            | M*q'' + [K_L + K_NL(q)]*q = F(t)   |
   | Giai bang: Newton-Raphson          |            | Dao dong tu do: Lap tri rieng      |
   | Output: Duong cong Tai - Do vong,  |            | Dao dong cuong buc: Runge-Kutta 4  |
   | Bien thien Ung suat phap, tiep     |            | Output: Tan so NL, Quy dao pha     |
   +------------------------------------+            +------------------------------------+

Data và phân tích

Toàn bộ thuật toán được tác giả tự lập trình hoàn chỉnh trên môi trường MATLAB (thể hiện chi tiết tại các Phụ lục B, C, D, E của luận án):

• Phân tích phi tuyến tĩnh: Giải hệ phương trình đại số phi tuyến $[K_L + K_{NL}(q)]q = F$ bằng thuật toán lặp Newton-Raphson với tiêu chí hội tụ sai số gia số chuyển vị $|\Delta q^{(k)}| / |q^{(k)}| \le 10^{-6}$.
• Phân tích dao động tự do phi tuyến: Sử dụng thuật toán lặp gần đúng dần (iterative eigenvalue solver) để trích xuất tỷ số tần số phi tuyến trên tuyến tính $\Omega_{NL}/\Omega_L$ tại các biên độ dao động đỉnh $W_{max}/h$ biến thiên từ $0,2$ đến $2,0$.
• Phân tích đáp ứng động cưỡng bức: Tích phân từng bước theo thời gian bằng phương pháp Runge-Kutta bậc 4 (RK4) với bước thời gian cực nhỏ $\Delta t = 10^{-4}\text{ s}$ để vẽ chính xác đáp ứng chuyển vị thời gian và quỹ đạo pha $(w, \dot{w})$.

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

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

1. Quy luật làm cứng phi tuyến vượt trội: Khi biên độ dao động không thứ nguyên $W_{max}/h$ tăng từ $0,2$ lên $2,0$, tỷ số tần số dao động $\Omega_{NL}/\Omega_L$ tăng mạnh từ $1,02$ lên $1,78$ (đối với biên CC, dầm sandwich $1\text{-}8\text{-}1$, $e_0=0,5$, $W_{GPL}=1%$). Điều này chứng minh kết cấu thể hiện đặc tính lò xo làm cứng (hardening spring behavior) rất rõ nét do ảnh hưởng kéo dài màng của biến dạng phi tuyến von Kármán.
2. Hiệu ứng tối ưu của cặp cấu hình PD1 - Dạng A: Dạng phân bố lỗ rỗng đối xứng PD1 (tập trung rỗng tại trục trung hòa) kết hợp với quy luật gia cường GPL Dạng A (tập trung GPL ở hai biên ngoài của lớp lõi tiếp giáp bề mặt FGM) mang lại độ cứng uốn cực đại. Độ võng tĩnh không thứ nguyên $\bar{w}$ của cấu hình này giảm tới $42,6%$ so với phân bố đều (UPD - Dạng C) dưới cùng một tải trọng phân bố $P_0 = 10$.
3. Cơ chế bù đắp độ cứng của hạt nano GPL: Việc bổ sung chỉ $W_{GPL} = 1,0%$ khối lượng GPL làm tăng tần số dao động riêng cơ bản của dầm lên hơn $118%$ so với trường hợp không gia cường ($W_{GPL} = 0%$), qua đó vô hiệu hóa hoàn toàn tác động làm mềm kết cấu của hệ số rỗng $e_0 = 0,5$.
4. Ảnh hưởng quyết định của tỷ lệ bề dày và điều kiện biên: Cấu hình dầm $1\text{-}8\text{-}1$ (lớp lõi dày gấp 8 lần tổng hai lớp bề mặt) thể hiện khả năng hấp thụ năng lượng động và tỷ số tần số phi tuyến cao nhất. Điều kiện biên ngàm hai đầu (CC) có độ cứng lớn nhất, làm suy giảm độ võng cực đại tới $78,5%$ so với điều kiện biên khớp (SS).
5. Hiện tượng dịch pha và phi tuyến hóa đáp ứng cưỡng bức: Khi tần số kích thích $\Omega$ tiệm cận tần số dao động riêng tuyến tính $\omega_L$, đáp ứng chuyển vị tuyến tính xảy ra hiện tượng cộng hưởng không giới hạn, trong khi đáp ứng phi tuyến duy trì biên độ hữu hạn và đường cong pha tạo thành chu trình giới hạn (limit cycle) ổn định trên không gian trạng thái.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp một khung toán học nhất thể hóa chính xác cho bài toán cơ học phi tuyến vật liệu composite nano dị hướng, đóng vai trò làm chuẩn tham chiếu cho việc hiệu chỉnh các mô hình số phức tạp.
• Về phương pháp luận: Chứng minh tính hiệu quả và độ ổn định cao của phương pháp bán giải tích Pb-Ritz trong việc giải quyết các bài toán phi tuyến hình học đa tham số, tránh được hiện tượng khóa cắt (shear locking) vốn thường gặp trong phương pháp phần tử hữu hạn bậc thấp.
• Về ứng dụng thực tiễn: Cung cấp hướng dẫn định lượng chính xác cho các kỹ sư kết cấu trong việc thiết kế tối ưu các cấu kiện chịu tải nhiệt-cơ nặng trong cánh máy bay siêu thanh, vỏ tàu ngầm và kết cấu sàn chịu rung động trong kỹ thuật xây dựng dân dụng và công nghiệp.

Limitations và Future Research

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

1. Giới hạn về liên kết tiếp xúc (Interfacial bonding): Luận án giả thiết các lớp liên kết hoàn hảo tuyệt đối (perfect bonding). Trong thực tế, dưới tải trọng mỏi chu kỳ biên độ lớn, hiện tượng trượt vi mô và suy giảm độ dính bám giữa lớp bề mặt FGM và lớp lõi FGP có thể xuất hiện.
2. Giới hạn về trường nhiệt độ và độ ẩm (Hygro-thermal effects): Mô hình chưa xét đến sự biến đổi cơ tính của vật liệu nền và GPL theo hàm nhiệt độ phi tuyến và độ ẩm môi trường hoạt động.
3. Giới hạn về hình học kết cấu: Nghiên cứu giới hạn ở phần tử dầm một chiều (1D).

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng chính:

• Phát triển mô hình cho kết cấu tấm/vỏ hai chiều (2D plate/shell) FGM-FGP-GPL chịu tác dụng của trường đa vật lý (nhiệt - điện - từ - cơ).
• Tích hợp lý thuyết đàn hồi không cục bộ (Nonlocal Elasticity Theory) và lý thuyết gradient biến dạng (Strain Gradient Theory) để phân tích các cấu kiện dầm/tấm nano FGP-GPL trong hệ thống vi cơ điện tử (MEMS/NEMS).
• Khảo sát thực nghiệm phá hủy uốn tĩnh và rung động thực tế để xây dựng hệ số hiệu chỉnh thực nghiệm cho mô hình Halpin-Tsai tế bào kín.
• Nghiên cứu cơ chế lan truyền vết nứt động và độ bền mỏi đa trục của dầm sandwich FGP-GPL.

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

• Tác động học thuật: Kết quả của luận án đóng góp các bộ số liệu chuẩn (benchmark data) có giá trị cao cho cộng đồng cơ học tính toán quốc tế. Dự báo công trình sẽ thu hút trích dẫn khoa học mạnh mẽ trên các tạp chí quốc tế chuyên ngành uy tín (ISI/Scopus Q1 như Composite Structures, Engineering Analysis with Boundary Elements, Thin-Walled Structures).
• Chuyển đổi công nghiệp: Định hình giải pháp vật liệu nhẹ mới cho ngành hàng không vũ trụ và giao thông vận tải, cho phép giảm từ $20%$ đến $35%$ trọng lượng kết cấu chịu lực mà vẫn đảm bảo đầy đủ các chỉ tiêu độ bền uốn và ổn định động học.
• Tác động kinh tế - xã hội: Giảm tiêu hao nhiên liệu và khí thải carbon của các phương tiện vận tải thông qua việc tối ưu hóa khối lượng cấu kiện; nâng cao tuổi thọ và độ an toàn của công trình xây dựng trước các tác động động đất và tải trọng gió bão.

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

• Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật/Cơ học kết cấu: Thụ hưởng một khung lý thuyết dầm tổng quát mẫu mực và mã nguồn tính toán MATLAB đã được thẩm định hoàn chỉnh.
• Các nhà khoa học và chuyên gia cơ học tính toán: Sử dụng các công thức giải tích và số liệu kiểm chứng để phát triển các phần tử hữu hạn mới hoặc các phương pháp không lưới tiên tiến.
• Kỹ sư thiết kế R&D tại các viện nghiên cứu hàng không, hàng hải và quốc phòng: Tiếp cận trực tiếp quy trình tối ưu hóa phân bố vi mô lỗ rỗng và hạt nano GPL để tạo ra các sản phẩm sandwich có khả năng chống va đập và tiêu tán năng lượng rung động vượt trội.
• Các nhà hoạch định chính sách và tiêu chuẩn hóa xây dựng: Sử dụng căn cứ khoa học định lượng từ luận án để xây dựng các tiêu chuẩn kỹ thuật thiết kế kết cấu composite thế hệ mới tại Việt Nam.

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

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

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích tổng quát dựa trên nguyên lý biến phân Hamilton kết hợp phi tuyến hình học von Kármán, bao quát đồng thời 6 lý thuyết dầm biến dạng cắt khác nhau (CBT, FSDBT, PSDBT, TSDBT, HSDBT, ESDBT) cho cấu trúc sandwich phức hợp 3 lớp FGM-FGP-GPL trên nền đàn hồi Pasternak. Luận án đã mở rộng mô hình vi cơ học Halpin-Tsai kết hợp quy tắc closed-cell GRF để giải quyết chính xác sự biến thiên cơ tính đa hướng qua chiều cao tiết diện.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu trước đây thể hiện ở đâu?

So với nghiên cứu của Frostig & Baruch [40] (chỉ dùng giải tích cho lõi mềm), Zenkour [167] (chỉ dùng nghiệm Navier cho dầm tựa đơn giản SS) và Chen et al. [23] (phương pháp Ritz cho dầm FGM xốp đơn giản), luận án đã phát triển phương pháp Pb-Ritz đa thức tổng quát hóa xử lý hoàn hảo mọi điều kiện biên phức tạp (CC, SS, CS, CF) mà không gặp hiện tượng khóa cắt. Luận án kết hợp đồng thời ba thuật toán số cao cấp: lặp Newton-Raphson cho tĩnh phi tuyến, lặp trị riêng phi tuyến cho dao động tự do và Runge-Kutta bậc 4 cho đáp ứng động lực học.

3. Phát hiện bất ngờ nhất từ các kết quả số là gì?

Phát hiện bất ngờ nhất là tính chất "làm cứng động lực học" phi tuyến rất mạnh của dầm ngàm hai đầu (CC) có lớp lõi phân bố theo dạng PD1 gia cường GPL Dạng A. Khi biên độ dao động đạt $W_{max}/h = 2,0$, tỷ số tần số $\Omega_{NL}/\Omega_L$ đạt tới $1,78$, trong khi ở dầm khớp hai đầu (SS) tỷ số này chỉ đạt $1,34$. Điều này chứng minh rằng việc kiềm chế chuyển vị dọc trục tại biên ngàm đã kích hoạt hiệu ứng lực kéo màng (membrane tension effect) cực lớn, làm tăng độ cứng toàn hệ thống một cách vượt bậc.

4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ các ma trận độ cứng tuyến tính $[K_L]$, ma trận độ cứng phi tuyến $[K_{NL}]$, ma trận khối lượng $[M]$ và véc tơ tải trọng $[F]$. Đặc biệt, tác giả công khai các chương trình nguồn viết trên nền MATLAB từ Phụ lục A đến Phụ lục E, bao gồm chương trình tính cơ tính hiệu dụng vi mô, chương trình phân tích tĩnh phi tuyến, chương trình dao động tự do và chương trình phân tích đáp ứng động lực học.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?

Chương trình nghiên cứu 10 năm tập trung vào lộ trình: (1) Mở rộng lý thuyết sang tấm và vỏ composite FGP-GPL chịu tải nhiệt-cơ-điện kết hợp; (2) Ứng dụng lý thuyết đàn hồi vi mô (Micropolar/Nonlocal Theory) để nghiên cứu cấu kiện nano; (3) Tối ưu hóa topo vi cấu trúc lỗ rỗng bằng trí tuệ nhân tạo (AI-driven structural optimization); (4) Tiến hành thử nghiệm chế tạo mẫu in 3D bọt kim loại gia cường graphene và kiểm chứng phá hủy thực tế.

Kết luận

1. Luận án đã giải quyết thành công và trọn vẹn bài toán phân tích phi tuyến tĩnh và dao động của dầm sandwich ba lớp FGM-FGP-GPL đặt trên nền đàn hồi Pasternak bằng phương pháp bán giải tích Pb-Ritz chính xác.
2. Xây dựng một trường chuyển vị tổng quát đại diện duy nhất, cho phép chuyển đổi linh hoạt và so sánh định lượng trực tiếp giữa 6 lý thuyết dầm biến dạng cắt khác nhau từ cổ điển đến bậc cao.
3. Thiết lập hệ thống thuật toán và mã nguồn MATLAB độc lập có độ tin cậy và tốc độ hội tụ cao, giải quyết đồng thời bài toán tĩnh phi tuyến (Newton-Raphson), dao động tự do phi tuyến (lặp trị riêng) và dao động cưỡng bức (Runge-Kutta 4).
4. Khám phá và lượng hóa quy luật tương tác tối ưu: phân bố lỗ rỗng đối xứng PD1 kết hợp với phân bố nano GPL Dạng A và hàm lượng $W_{GPL} = 1,0% \div 2,0%$ mang lại độ cứng uốn và tần số dao động cao nhất cho kết cấu sandwich.
5. Chứng minh rõ nét hiệu ứng lò xo làm cứng phi tuyến hình học von Kármán và vai trò kiểm soát chuyển vị của tương tác nền đàn hồi Pasternak dưới các điều kiện biên khác nhau.
6. Mở ra ba hướng nghiên cứu mới đầy triển vọng về kết cấu tấm/vỏ nano composite đa trường vật lý, đóng góp nền tảng lý thuyết và thực tiễn vững chắc cho sự phát triển của ngành Cơ kỹ thuật và Công nghệ vật liệu tiên tiến tại Việt Nam và trên thế giới.