Tổng quan về luận án

Công trình nghiên cứu tiến sĩ kỹ thuật mang tên "Phân tích ứng xử tĩnh và dao động kết cấu dầm sandwich có lớp lõi bằng vật liệu rỗng" (Chuyên ngành Kỹ thuật Xây dựng, Mã số: 9580201) do nghiên cứu sinh Đỗ Minh Đức thực hiện dưới sự hướng dẫn khoa học của PGS.TS Trần Quang Hưng và GS.TS Trần Minh Tú tại Trường Đại học Bách khoa – Đại học Đà Nẵng (2025) đại diện cho bước phát triển tiên phong trong cơ học tính toán công trình và kết cấu thông minh tại Việt Nam.

Bối cảnh khoa học của luận án xuất phát từ yêu cầu tối ưu hóa tỷ số độ cứng trên trọng lượng (stiffness-to-weight ratio), tăng cường khả năng tiêu tán năng lượng, cách âm, cách nhiệt và chịu nhiệt độ cao trong các ngành kỹ thuật mũi nhọn như hàng không vũ trụ, giao thông vận tải đường sắt tốc độ cao, 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 với lõi đồng nhất thường bộc lộ hạn chế nghiêm trọng về sự tập trung ứng suất cắt, hiện tượng phân lớp (delamination) và thiếu tính thích ứng chủ động. Luận án đã giải quyết trực diện khoảng trống nghiên cứu (research gap) khi kết hợp đồng thời lớp lõi vật liệu rỗng có cơ tính biến thiên được gia cường tiểu bản graphene (FGP-GPLRC) cùng hai lớp bề mặt bằng vật liệu có cơ tính biến thiên (FGM) hoặc lớp vật liệu thông minh áp điện (PEM) đóng vai trò cảm biến và cơ cấu chấp hành (sensor/actuator).

Các câu hỏi nghiên cứu cốt lõi được định hình gồm:

  1. RQ1: Quy luật phân bố lỗ rỗng và graphene (GPL) theo chiều dày ảnh hưởng định lượng như thế nào đến trường ứng suất, biến dạng và độ võng tĩnh của dầm sandwich đa lớp?
  2. RQ2: Đặc trưng dao động riêng và đáp ứng động lực học của dầm thay đổi ra sao dưới kích thích tải cơ - điện kết hợp và các dạng xung tải trọng khác nhau?
  3. RQ3: Hiệu quả triệt tiêu dao động chủ động của lớp phủ PEM thông qua thuật toán hồi tiếp vận tốc đạt được ở mức độ nào và nguyên nhân gây mất ổn định điều khiển (instability) là gì?
  4. RQ4: Độ chính xác và hiệu quả hội tụ của các phương pháp không lưới (MFM) 1D so với nghiệm giải tích Navier và kết quả thực nghiệm dầm thực tế như thế nào?

Các giả thuyết nghiên cứu chính được đặt ra và kiểm chứng:

  • H1: Việc tích hợp GPL phân bố tập trung ở bề mặt ngoài của lớp lõi (kiểu GPL-X/GPL-S) bù đắp vượt trội sự suy giảm độ cứng do các lỗ rỗng gây ra, nâng cao tần số dao động tự nhiên của dầm.
  • H2: Lý thuyết biến dạng cắt bậc cao của Reddy (HSDT) kết hợp nguyên lý Hamilton và phương pháp không lưới xây dựng hàm dạng $C^1$-Hermite mô tả chính xác ứng xử cơ - điện mà không cần sử dụng hệ số hiệu chỉnh cắt.
  • H3: Vòng điều khiển hồi tiếp vận tốc chủ động có khả năng triệt tiêu hoàn toàn hiện tượng cộng hưởng và phách điều hòa khi thiết kế bố trí các cặp điện cực áp điện phù hợp.

Khung lý thuyết của công trình tích hợp cơ học môi trường liên tục, mô hình vi cơ học hỗn hợp cải tiến Halpin–Tsai cho vật liệu nano GPL, lý thuyết áp điện tuyến tính và lý thuyết dầm bậc ba Reddy (HSDT 1D). Phạm vi nghiên cứu bao quát phân tích giải tích Navier, mô phỏng số không lưới (PIM, RPIM, MLS) trên nhiều điều kiện biên (SS, CC, CS, CF) và xác thực thực nghiệm trên mẫu dầm sandwich lõi bê tông hạt xốp EPS kết hợp tấm bề mặt SCG Cemboard từ nhà sản xuất NUCEWALL.

Literature Review và Positioning

Nghiên cứu về vật liệu có cơ tính biến thiên (FGM) khởi nguồn từ năm 1984 tại Nhật Bản trong các dự án hàng không vũ trụ nhằm đáp ứng khả năng chịu gradient nhiệt cực hạn "lên đến 1000 K trên mặt cắt có chiều dày bé hơn 10 mm" (Saleh et al. [5]). Từ các quy luật biến thiên hàm lũy thừa (power-law), hàm mũ (exponential) đến hàm sigmoid (Chi và Chung [9]), FGM đã giải quyết triệt để sự gián đoạn cơ tính tại mặt tiếp xúc lớp. Tiếp nối FGM, vật liệu rỗng biến thiên (FGP) với cấu trúc tế bào mở (open-cell) và tế bào kín (closed-cell) được tổng kết toàn diện bởi Ashby et al. [10] và Chen et al. [13], chứng minh khả năng giảm thiểu tải trọng bản thân và hấp thụ va đập vượt bậc.

Tuy nhiên, sự xuất hiện của các lỗ rỗng làm suy giảm đáng kể mô-đun đàn hồi tổng thể của kết cấu. Để khắc phục điểm yếu này, Kitipornchai et al. [19] đã mở ra bước ngoặt đột phá khi chứng minh việc "bổ sung 1% tỉ lệ khối lượng GPL, tần số dao động riêng cơ bản tăng lên đến hơn 41%, lực tới hạn tăng hơn 85%". Các nghiên cứu sau đó của Yang et al., Chen et al. [37] đã khai thác GPL trong các cấu trúc dầm và tấm đơn lớp.

Trong lĩnh vực kết cấu sandwich lõi FGP, các nghiên cứu quốc tế trước đây chủ yếu tập trung vào các lớp bề mặt đẳng hướng hoặc FGM thông thường:

  • Magnucki và Magnucka-Blandzi [30] phát triển lời giải giải tích uốn tĩnh cho dầm sandwich lõi FGP với hai mặt đẳng hướng.
  • Bamdad et al. [34] và Biglari et al. [36] áp dụng lý thuyết dầm quasi-3D khảo sát dao động dầm sandwich lõi FGP chịu tải trọng di động và trường từ tính.
  • Wang et al. [39, 43] sử dụng phương pháp Chebyshev–Ritz giải bài toán uốn phi tuyến dầm sandwich FGM/FGP.
  • Derikvand et al. [42] dùng phương pháp biến đổi vi phân (DTM) khảo sát mất ổn định của dầm cao sandwich FGM-FGP.

Tại Việt Nam, các nhóm nghiên cứu mạnh như Đặng Xuân Hùng et al. [68, 75], Chu Thanh Bình et al. [69, 70], Nguyễn Văn Long et al. [71], Vũ Thị An Ninh [72], Nguyễn Ngọc Dương et al. [73, 74], và Thái Hoàng Chiến et al. [81, 82] đã có những công bố quốc tế uy tín về dầm và tấm FGM/FGP đặt trên nền đàn hồi Pasternak/Winkler. Về phương pháp không lưới (MFM), các đóng góp tiêu biểu của Bùi Quốc Tính et al. [77], Tôn Thất Hoàng Lân [78], Vũ Văn Tân et al. [79, 80], và Trương Tích Thiện et al. [83, 84] đã khẳng định tính ưu việt của MFM trong cơ học tấm rỗng và tấm composite nứt.

Điểm định vị và bước tiến vượt bậc của luận án: Đa số các công trình trước đây chỉ khảo sát dầm FGP đơn lớp, hoặc dầm sandwich với lõi FGP thuần túy thiếu vắng hạt gia cường GPL, hoặc chưa tích hợp đồng thời trường vật liệu áp điện PEM để thiết lập hệ cơ điện thông minh kiểm soát dao động. Luận án của Đỗ Minh Đức đã tiên phong định vị nghiên cứu tại giao điểm của ba trụ cột: cấu trúc lõi siêu nhẹ FGP-GPLRC đa pha, lớp bề mặt thông minh PEM/FGM chịu nhiệt và điều khiển chủ động, cùng hệ thống phương pháp số không lưới đa dạng (PIM/RPIM/MLS) kết hợp kiểm chứng thực nghiệm trực tiếp.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp trực tiếp vào sự mở rộng lý thuyết dầm bậc cao của J.N. Reddy (Higher-Order Shear Deformation Theory - HSDT) bằng cách thiết lập mô hình tương tác cơ - điện đầy đủ cho kết cấu nhiều lớp không đối xứng. Bằng việc phân phối trường chuyển vị dọc trục $u(x,z,t)$ và độ võng $w(x,z,t)$ theo hàm phi tuyến bậc ba của tọa độ chiều dày $z$: $$u(x,z,t) = u_0(x,t) + z \theta(x,t) - \frac{4z^3}{3h^2}\left(\theta(x,t) + \frac{\partial w_0(x,t)}{\partial x}\right)$$ $$w(x,z,t) = w_0(x,t)$$ Lý thuyết này thỏa mãn điều kiện ứng suất tiếp bằng 0 tại các bề mặt tự do trên và dưới ($z = \pm h/2$), loại bỏ sự phụ thuộc vào hệ số hiệu chỉnh cắt vốn gây nhiều sai số trong lý thuyết dầm Timoshenko (FSDT).

+-------------------------------------------------------------------------+
|                  Lớp bề mặt trên (Top Facesheet): PEM / FGM             |
|                         (Chiều dày: h_f, z_2 đến z_3)                   |
+-------------------------------------------------------------------------+
|                                                                         |
|                 LỚP LÕI SIÊU NHẸ (Core): FGP-GPLRC                      |
|             (Phân bố lỗ rỗng: Poro-U, Poro-S, Poro-A)                  |
|             (Phân bố hạt nano: GPL-U, GPL-S, GPL-A)                     |
|                         (Chiều dày: h_c, z_1 đến z_2)                   |
|                                                                         |
+-------------------------------------------------------------------------+
|                 Lớp bề mặt dưới (Bottom Facesheet): PEM / FGM           |
|                         (Chiều dày: h_f, z_0 đến z_1)                   |
+-------------------------------------------------------------------------+

Mô hình hóa cơ tính hiệu dụng của lớp lõi FGP-GPLRC được liên kết thông qua quy tắc pha trộn mở rộng và mô hình vi cơ học Halpin–Tsai: $$E_{mf}(z) = E_0 [1 - e_0 \Phi(z)]$$ $$\rho(z) = \rho_0 [1 - e_p \Phi(z)]$$ trong đó $e_0$ và $e_p$ lần lượt là hệ số độ rỗng và hệ số độ rỗng khối lượng, $\Phi(z)$ là hàm phân bố độ rỗng (Poro-U: đều; Poro-S: đối xứng; Poro-A: bất đối xứng). Mô-đun đàn hồi hiệu dụng $E_c(z)$ sau khi gia cường GPL được xác định chính xác theo công thức Halpin–Tsai có tính đến kích thước hình học trung bình ($l_{GPL}, w_{GPL}, t_{GPL}$) và tỷ lệ thể tích $V_{GPL}(z)$.

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

Khung phân tích của luận án tích hợp thành công ba nền tảng lý thuyết:

  1. Nguyên lý biến phân Hamilton: $\int_{t_1}^{t_2} (\delta K - \delta U + \delta V) dt = 0$, thiết lập hệ phương trình vi phân đạo hàm riêng liên kết chuyển vị và điện thế tổng quát.
  2. Lý thuyết áp điện liên kết (Linear Piezoelectricity): Mô tả trường dịch chuyển điện tích $D_z = e_{31} \varepsilon_x + p_{33} E_z$ kết hợp hiệu ứng áp điện nghịch nhằm tạo ra mô-men uốn chủ động kháng lại dao động ngoại vi.
  3. Thuật toán điều khiển phản hồi vận tốc (Velocity Feedback Control): Tạo ra ma trận cản chủ động bổ sung $C_a = G_v \mathbf{K}{u\phi} \mathbf{K}{\phi\phi}^{-1} \mathbf{K}_{\phi u}$, trực tiếp can thiệp vào các mode dao động riêng.

Điều kiện biên và giới hạn áp dụng (boundary conditions) được xác lập chặt chẽ cho cả 4 trạng thái liên kết cơ học: Tựa đơn giản hai đầu (SS), Ngàm hai đầu (CC), Đầu ngàm đầu khớp (CS), và Ngàm - Tự do (CF), kết hợp với điều kiện hở mạch hoặc kín mạch điện trường mặt.

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

Thiết kế nghiên cứu

Luận án được thiết kế theo trường phái thực chứng thực nghiệm kết hợp duy lý diễn dịch (Positivist / Post-positivist Computational Paradigm). Quy trình kết hợp đa phương pháp (triangulation) bao gồm:

  • Lời giải giải tích chính xác (Analytical Navier Method): Phát triển chuỗi Fourier vô hạn cho các dầm tựa đơn giản SS, đóng vai trò là chuẩn đối sánh (benchmark).
  • Phương pháp số không lưới (Mesh-Free Methods - MFM): Xây dựng mô hình 1D giải quyết bài toán với điều kiện biên tùy ý và hình học phức tạp.
  • Nghiên cứu thực nghiệm vật lý (Physical Experiments): Đo đạc ứng xử tĩnh và động trên dầm sandwich thực tế để thẩm định mô hình toán.

Quy trình nghiên cứu rigorous

Quy trình tính toán không lưới MFM được phát triển với ba kỹ thuật xấp xỉ trường chuyển vị và điện thế:

  1. Phương pháp nội suy điểm dùng hàm đa thức cơ sở (Point Interpolation Method - PIM): Sử dụng cơ sở đa thức bậc thấp giải quyết bài toán qua ma trận nghịch đảo cục bộ.
  2. Phương pháp nội suy điểm dùng hàm bán kính cơ sở (Radial Point Interpolation Method - RPIM): Kết hợp hàm đa thức với hàm bán kính đa bậc hai (Multi-quadrics - MQ) có dạng $R_i(x) = (r_i^2 + (C \cdot d_c)^2)^q$, với tham số hình dạng tối ưu $C = 1.0 \div 3.0$ và $k_{sc} = 3.0 \div 4.0$, loại bỏ hoàn toàn hiện tượng suy biến ma trận hình thức (singularity).
  3. Phương pháp bình phương cực tiểu động (Moving Least Squares - MLS): Sử dụng hàm trọng số spline bậc 4 và kỹ thuật nhân tử Lagrange hoặc hàm phạt để thỏa mãn điều kiện biên thiết yếu.

Đặc biệt, luận án đề xuất kỹ thuật xây dựng hàm dạng $C^1$-Hermite duy nhất cho tất cả các biến chuyển vị cơ học ($u_0, w_0, \theta$) và hàm thế điện ($\phi_0$), đảm bảo tính liên tục của độ dốc biến dạng mà không làm tăng bậc của bài toán.

Quy trình tích hợp miền nền (background cells) sử dụng tích phân số Gauss 4 điểm và 6 điểm cho mỗi ô tính toán. Đối với bài toán động học theo thời gian, phương pháp tích phân từng bước Newmark với giả thiết gia tốc không đổi ($\gamma = 0.5, \beta = 0.25$) được lập trình giải hệ phương trình động lực học tổng thể: $$\mathbf{M} \ddot{\mathbf{u}} + (\mathbf{C} + \mathbf{C}_a) \dot{\mathbf{u}} + \mathbf{K} \mathbf{u} = \mathbf{F}_m(t) + \mathbf{F}_e(t)$$

                        QUY TRÌNH NGHIÊN CỨU TOÀN DIỆN
                        
    +--------------------------------------------------------------------+
    |                     MÔ HÌNH VẬT LÝ DẦM SANDWICH                    |
    |      Lõi: FGP-GPLRC (Poro-U/S/A, GPL-U/S/A) | Vỏ: FGM / PEM        |
    +---------------------------------+----------------------------------+
                                      |
            +-------------------------+-------------------------+
            |                                                   |
            v                                                   v
+-----------------------+                           +-----------------------+
|  NGHIÊN CỨU LÝ THUYẾT |                           |  NGHIÊN CỨU THỰC NGHIỆM|
|  - HSDT Reddy (1D)    |                           |  - Mẫu dầm NUCEWALL   |
|  - Nguyên lý Hamilton |                           |  - Lõi bê tông EPS    |
+-----------+-----------+                           |  - Mặt SCG Cemboard   |
            |                                       +-----------+-----------+
      +-----+-----+                                             |
      |           |                                             |
      v           v                                             |
+-----------+ +-----------+                                     |
| Giải tích | | Không lưới|                                     |
|  Navier   | | MFM (1D)  |                                     |
|  (SS BCs) | |PIM/RPIM/MLS|                                    |
+-----+-----+ +-----+-----+                                     |
      |             |                                           |
      +------+------+                                           |
             |                                                  |
             v                                                  v
    +--------------------------------------------------------------------+
    |      SO SÁNH - ĐỐI SOÁT - THẨM ĐỊNH (TRIANGULATION & VALIDATION)   |
    |      - Sai số võng tĩnh < 4.2% | Sai số tần số dao động < 3.8%     |
    +--------------------------------------------------------------------+

Data và phân tích

Nghiên cứu thực nghiệm được tiến hành nghiêm ngặt tại phòng thí nghiệm:

  • Mẫu dầm thực nghiệm: Dầm sandwich kích thước thực tế được cắt từ tấm tường nhẹ tiền chế NUCEWALL; chiều dài tính toán nhịp $L = 1100\text{ mm}$, chiều rộng $b = 100\text{ mm}$, chiều cao $h = 75\text{ mm}$ (lớp lõi bê tông hạt xốp EPS dày $hc = 57\text{ mm}$, hai lớp bề mặt Cemboard SCG Thái Lan dày $hf = 9\text{ mm}$).
  • Đặc tính vật liệu thực nghiệm: Bê tông xốp EPS có khối lượng riêng siêu nhẹ $\rho \approx 600 - 800\text{ kg/m}^3$, mô-đun đàn hồi $E_{core} \approx 650\text{ MPa}$; tấm SCG Cemboard có $E_{face} \approx 4500\text{ MPa}$.
  • Thiết bị đo lường & Phần mềm: Hệ thống thu nhận dữ liệu đa kênh tĩnh/động kết hợp cảm biến đo chuyển vị điện tử LVDT độ nhạy cao và đầu đo gia tốc áp điện kết nối máy tính qua phần mềm DCS-100A; phân tích phổ tần số dao động riêng bằng phép biến đổi Fourier nhanh (FFT) trên phần mềm chuyên dụng DAS-200A.

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

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

  1. Hiệu ứng cấu trúc rỗng đối xứng (Poro-S) kết hợp phân bố nano đối xứng (GPL-S): Kết quả giải tích và số chứng minh rằng kiểu phân bố lỗ rỗng đối xứng Poro-S (độ rỗng cực đại tại trục trung hòa $z=0$ và đặc chắc dần về hai phía bề mặt) kết hợp với kiểu phân bố GPL-S (hàm lượng graphene tập trung cao nhất ở lớp ngoài của lõi) mang lại hiệu quả chịu lực tối ưu nhất. Cụ thể, khi hệ số độ rỗng tăng từ $e_0 = 0.0$ lên $e_0 = 0.6$, độ võng tĩnh của dầm tăng nhưng nhờ có 1.0% khối lượng GPL-S, độ võng cực đại giảm trở lại đến 58.4% so với trường hợp lõi rỗng không gia cường, khôi phục hoàn toàn độ cứng uốn tổng thể.

  2. Sự triệt tiêu ứng suất cắt tiếp xúc mặt tiếp giáp: Khi sử dụng hai lớp bề mặt FGM biến thiên theo hàm lũy thừa ($k = 0.2 \div 2.0$), ứng suất tiếp $\tau_{xz}$ phân bố liên tục và trơn tuyệt đối qua các mặt phân chia lớp ($z = \pm h_c/2$). Điều này giải quyết tận gốc nhược điểm bong tách lớp kinh điển của dầm composite dán laminate truyền thống.

  3. Khả năng kiểm soát và dập tắt dao động chủ động qua lớp PEM: Dưới tác dụng của tải trọng xung hình chữ nhật và xung hình sin chu kỳ ngắn, khi kích hoạt hệ số hồi tiếp vận tốc $G_v = 1.0 \times 10^{-4} \div 5.0 \times 10^{-4}\text{ (V}\cdot\text{s/m)}$, biên độ dao động cưỡng bức của dầm giảm hơn 80% chỉ sau 3 chu kỳ dao động đầu tiên, đồng thời triệt tiêu hoàn toàn hiện tượng phách điều hòa (beating phenomenon) khi tần số kích thích tiệm cận tần số dao động tự nhiên của kết cấu ($\omega \approx \omega_1$).

  4. Hiện tượng mất ổn định thuật toán điều khiển (Control Instability) trên dầm bất đối xứng: Một phát hiện quan trọng có tính cảnh báo cao của luận án: Khi dầm có cấu hình phân bố vật liệu bất đối xứng (Poro-A hoặc GPL-A) với điều kiện biên ngàm hai đầu (CC), việc đặt các cặp điện cực sensor/actuator đối xứng ở hai mặt trên/dưới với hệ số khuếch đại $G_v > 0$ có thể kích hoạt hiện tượng "spillover" làm hệ thống mất ổn định động học (dao động phân kỳ). Luận án đã đưa ra giải pháp đột phá: tái thiết kế bố trí cặp sensor/actuator tập trung về một phía của dầm chủ, khôi phục tính ổn định tiệm cận hoàn toàn cho hệ thống điều khiển.

  5. Sự tương thích cao giữa thực nghiệm và mô hình MFM: Tần số dao động riêng thứ nhất $f_1$ đo đạc thực nghiệm qua phân tích FFT trên phần mềm DAS-200A đạt giá trị $42.50\text{ Hz}$, so sánh với tần số dự đoán từ mô hình số không lưới MFM là $44.12\text{ Hz}$, sai số tương đối chỉ ở mức $3.81%$. Độ võng uốn tĩnh đàn hồi giữa lý thuyết và đo đạc qua LVDT có độ tương đồng với sai số dưới $4.2%$, khẳng định độ tin cậy tuyệt đối của mô hình dầm HSDT 1D.

Tham số khảo sát Trường hợp dầm cơ sở (Poro-U, GPL-U) Trường hợp tối ưu (Poro-S, GPL-S) Mức độ cải thiện cơ học
Độ võng tĩnh giữa nhịp ($\bar{w}$) $1.452 \times 10^{-2}$ $0.604 \times 10^{-2}$ Giảm 58.4% biến dạng uốn
Tần số riêng cơ bản ($\bar{\omega}_1$) $3.892$ $5.518$ Tăng 41.8% độ cứng động lực học
Hệ số cản dập tắt dao động ($G_v$) $0.0\text{ (Không điều khiển)}$ $5.0 \times 10^{-4}\text{ V}\cdot\text{s/m}$ Giảm 82.5% biên độ sau 3 chu kỳ
Sai số thực nghiệm vs MFM ($f_1$) $44.12\text{ Hz (MFM)}$ $42.50\text{ Hz (Thực nghiệm)}$ Độ chính xác cao, sai số 3.81%

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định tính hiệu quả của việc kết hợp mô hình vi cơ học Halpin–Tsai với lý thuyết dịch chuyển cắt bậc cao Reddy trong việc mô tả vật liệu nano đa pha.
  • Về mặt phương pháp luận: Cung cấp bộ công cụ mã nguồn tính toán MFM (PIM/RPIM/MLS) dạng 1D hoàn chỉnh trên nền tảng MATLAB, cho phép giải quyết bài toán dầm đa lớp với chi phí tính toán chỉ bằng 15-20% so với mô hình phần tử hữu hạn 3D FEM truyền thống mà vẫn đạt độ chính xác tương đương.
  • Về ứng dụng thực tiễn: Cung cấp cơ sở khoa học trực tiếp cho các doanh nghiệp sản xuất cấu kiện tiền chế (như tấm tường sandwich NUCEWALL) tối ưu hóa tỷ lệ hạt xốp EPS, tăng cường độ chịu uốn và cách âm cho công trình xây dựng dân dụng và công nghiệp.

Limitations và Future Research

Mặc dù đạt được nhiều 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 nhất định:

  1. Nghiên cứu thực nghiệm mới chỉ dừng lại ở dầm sandwich vật liệu đẳng hướng (lõi bê tông EPS, mặt Cemboard), chưa thể chế tạo trực tiếp mẫu thực nghiệm dầm nano FGP-GPLRC tích hợp lớp áp điện PEM do hạn chế về công nghệ chế tạo vật liệu trong nước.
  2. Mô hình toán học mới khảo sát ứng xử tuyến tính hình học và đàn hồi tuyến tính của vật liệu, chưa xét đến biến dạng phi tuyến lớn von Kármán hay ứng xử từ biến/dẻo của vật liệu polymer xốp.
  3. Chưa tích hợp ảnh hưởng đồng thời của môi trường nhiệt độ cao (gradient nhiệt độ) và độ ẩm phức tạp (hygro-thermal environments) vào mô hình dao động động lực học.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Phát triển mô hình phân tích phi tuyến hình học và phi tuyến vật liệu cho tấm/vỏ sandwich FGP-GPLRC 2D chịu tải trọng nổ và va chạm vận tốc cao.
  • Nghiên cứu ứng xử phá hủy mỏi, cơ học rạn nứt tế bào rỗng và hiện tượng tách lớp cục bộ dưới tải trọng lặp chu kỳ lớn.
  • Áp dụng các thuật toán điều khiển thông minh phi tuyến hiện đại như Linear Quadratic Gaussian (LQG), $H_\infty$ và điều khiển thích nghi mờ (Fuzzy-Neural Control) để nâng cao độ bền vững của cơ cấu chấp hành PEM.
  • Ứng dụng công nghệ in 3D laser kim loại đa vật liệu để chế tạo mẫu dầm FGP-GPLRC thực thụ và kiểm chứng thực nghiệm đa trường vật lý.

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

Luận án tạo ra tác động sâu rộng trên nhiều bình diện:

  • Ảnh hưởng học thuật (Academic Impact): Đóng góp trực tiếp vào kho tàng cơ học tính toán vật liệu mới; các thuật toán MFM và nghiệm Navier đề xuất có tiềm năng trích dẫn cao trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/Scopus (Q1) như Composite Structures, Thin-Walled Structures, Engineering Structures, Applied Mathematical Modelling.
  • Chuyển đổi công nghiệp (Industry Transformation): Cung cấp giải pháp công nghệ trực tiếp cho ngành công nghiệp đóng tàu vỏ nhẹ, vỏ toa xe lửa cao tốc Bắc - Nam tương lai, cánh tuabin gió và các tấm sàn panel siêu nhẹ cách âm cách nhiệt trong xây dựng module tiền chế (modular construction).
  • Ý nghĩa xã hội & Môi trường: Thúc đẩy sử dụng vật liệu xây dựng xanh, giảm tiêu hao cốt liệu tự nhiên thông qua việc thay thế bằng hạt xốp EPS tái chế và bê tông rỗng, góp phần hiện thực hóa mục tiêu Net Zero Carbon trong ngành xây dựng.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận tài liệu tham khảo mẫu mực về quy trình toán học biến phân, kỹ thuật lập trình phương pháp không lưới (RPIM, MLS) và cách thức thiết lập đo đạc thực nghiệm động lực học.
  • Các nhà nghiên cứu cơ học & kỹ sư kết cấu: Sở hữu công thức giải tích dạng kín và các biểu đồ tham số chuẩn để tra cứu, kiểm tra nhanh độ võng, ứng suất và tần số dao động khi thiết kế kết cấu dạng thanh dầm sandwich.
  • Doanh nghiệp sản xuất vật liệu mới (R&D Departments): Nắm bắt nguyên lý phân bố tối ưu của độ rỗng và chất gia cường (Poro-S, GPL-S) để định hình công nghệ ép khuôn, in 3D hay đúc tấm panel composite đạt độ bền cơ học cao nhất với chi phí vật liệu nhỏ nhấ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ì?

Luận án đã mở rộng thành công lý thuyết dầm biến dạng cắt bậc cao HSDT của Reddy cho hệ dầm sandwich lai ba lớp (lõi FGP-GPLRC đa pha và hai mặt PEM/FGM) dưới dạng giải tích Navier và số không lưới MFM, thiết lập trọn vẹn mô hình tương tác đa trường cơ - điện - gradient vật liệu mà không cần sử dụng hệ số hiệu chỉnh cắt.

2. Điểm cải tiến cốt lõi của phương pháp số MFM trong nghiên cứu so với các công trình trước?

So với các nghiên cứu MFM truyền thống sử dụng FSDT (như Mesmoudi et al. [58, 59]), luận án đã xây dựng thành công kỹ thuật hàm dạng $C^1$-Hermite đồng nhất cho cả 3 thành phần chuyển vị và điện thế trên nền tảng RPIM với hàm bán kính MQ tối ưu ($C=2.0, k_{sc}=3.0$), khắc phục triệt để hiện tượng khóa cắt (shear locking) và hiện tượng suy biến ma trận nội suy.

3. Phát hiện bất ngờ nhất từ dữ liệu mô phỏng và thực nghiệm là gì?

Việc chỉ ra hiện tượng mất ổn định (instability) của thuật toán điều khiển hồi tiếp vận tốc chủ động $G_v$ khi cấu hình dầm có sự bất đối xứng vật liệu (Poro-A/GPL-A) và đề xuất giải pháp kỹ thuật đảo chiều cấu hình cảm biến/chấp hành để đảm bảo hệ thống luôn ổn định động học tuyệt đối.

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

Luận án cung cấp chi tiết toàn bộ thuật toán Newmark gia tốc không đổi, mã nguồn mẫu MATLAB xây dựng hàm dạng PIM, RPIM, MLS, thông số vật liệu chi tiết của dầm NUCEWALL, tấm Cemboard SCG, thiết lập thực nghiệm và cấu hình phần mềm DCS-100A/DAS-200A tại phần phụ lục.

5. Nghiên cứu có vạch ra lộ trình khoa học dài hạn 10 năm không?

Công trình đã phác thảo chi tiết lộ trình 10 năm từ mô hình dầm 1D tuyến tính tiến lên bản/vỏ 2D phi tuyến hình học, phân tích đa trường nhiệt - ẩm - cơ - điện và thương mại hóa sản phẩm tấm panel sandwich nhẹ chịu lực cao thông qua công nghệ chế tạo in 3D và vật liệu nano.

Kết luận

  1. Thiết lập thành công hệ thống nghiệm giải tích Navier và phương pháp số không lưới MFM (PIM, RPIM, MLS) dựa trên lý thuyết dầm biến dạng cắt bậc cao HSDT cho kết cấu dầm sandwich lõi rỗng biến thiên gia cường graphene (FGP-GPLRC) kết hợp lớp mặt FGM/PEM.
  2. Chứng minh định lượng rằng quy luật phân bố độ rỗng đối xứng Poro-S kết hợp phân bố nano GPL-S mang lại độ cứng uốn và tần số dao động riêng cao nhất, giúp giảm độ võng tĩnh tới 58.4% và tăng tần số dao động riêng cơ bản lên hơn 41.8%.
  3. Phát triển thành công thuật toán điều khiển chủ động qua lớp phủ áp điện PEM, dập tắt hơn 80% biên độ dao động dưới các dạng tải trọng xung cơ - điện phức tạp, đồng thời làm rõ và xử lý triệt để hiện tượng mất ổn định thuật toán trên kết cấu bất đối xứng.
  4. Thực hiện thành công chương trình thực nghiệm đo độ võng tĩnh và tần số dao động động lực học trên mẫu dầm sandwich bê tông hạt xốp EPS thực tế từ nhà sản xuất NUCEWALL, ghi nhận độ chính xác xuất sắc với sai số dưới 4.2% so với mô hình số lý thuyết.
  5. Mở ra hướng tiếp cận hoàn chỉnh, tin cậy và hiệu quả cao cho lĩnh vực cơ học tính toán công trình, đặt nền móng khoa học vững chắc cho việc tính toán, thiết kế và ứng dụng các kết cấu sandwich thông minh siêu nhẹ trong kỹ thuật xây dựng và công nghiệp hiện đại.