Tổng quan về luận án

Trong kỷ nguyên phát triển mạnh mẽ của khoa học vật liệu và kỹ thuật công trình thế kỷ XXI, 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à tích hợp chức năng tự điều khiển thông minh đã trở thành bài toán trung tâm của cơ học kết cấu hiện đại. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng (Mã số: 9580201) của nghiên cứu sinh Đỗ Minh Đức, được thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Trần Quang Hưng (Trường Đại học Bách khoa – Đại học Đà Nẵng) và GS.TS. Trần Minh Tú (Trường Đại học Xây dựng Hà Nội) với đề tài: "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" (2025), đã giải quyết một cách hệ thống những thách thức cốt lõi nêu trên.

flowchart TD
    A["Kết cấu Dầm Sandwich Đa Năng Tiên Tiến"] --> B["Lớp Bề Mặt (Face Sheets)"]
    A --> C["Lớp Lõi Chịu Cắt (Core Layer)"]
    
    B --> B1["Vật liệu FGM: Chống chịu nhiệt độ cao, ngăn bong tách lớp"]
    B --> B2["Vật liệu PEM: Cảm biến và điều khiển dao động chủ động"]
    
    C --> C1["Vật liệu FGP: Giảm 30-50% trọng lượng, tiêu tán xung lực"]
    C --> C2["Gia cường GPL: Bổ sung 1% wt tăng >41% tần số dao động"]
    
    A --> D["Khung Phương Pháp Luận Tích Hợp"]
    D --> D1["Giải tích Navier (Closed-form)"]
    D --> D2["Phương pháp Không lưới MFM (PIM, RPIM, MLS)"]
    D --> D3["Thực nghiệm Thực tế (EPS Concrete + Cemboard)"]

Bối cảnh khoa học của luận án xuất phát từ sự chuyển dịch từ kết cấu nhiều lớp truyền thống (laminated composites) sang kết cấu lai thông minh đa chức năng. Kết cấu sandwich kinh điển với lớp lõi nhẹ kẹp giữa hai bản mặt mỏng cường độ cao đóng vai trò tương tự như tiết diện chữ I trong cơ học kết cấu: "lớp lõi dày có trọng lượng nhẹ đặt giữa hai lớp bề mặt làm bằng vật liệu có cường độ cao... hai lớp vật liệu bề mặt đưa ra xa đường trung hòa, đóng vai trò chịu lực chính... lớp lõi đóng vai trò là lớp cấu tạo do vậy thường làm từ vật liệu nhẹ, rỗng để giảm trọng lượng kết cấu và tăng tính cách âm, cách nhiệt, hấp thụ năng lượng". Tuy nhiên, sự xuất hiện của các gradient ứng suất lớn tại mặt tiếp xúc giữa các lớp đồng nhất dễ dẫn đến hiện tượng phá hủy phân lớp (delamination). Để khắc phục triệt để điểm nghẽn cơ học này, luận án tiên phong phát triển mô hình dầm sandwich lai tích hợp:

  1. Lớp lõi bằng vật liệu rỗng có cơ tính biến thiên gia cường phiến nano graphene (Functionally Graded Porous - Graphene Platelet Reinforced Composite: FGP-GPLRC).
  2. Hai lớp bề mặt bằng vật liệu có cơ tính biến thiên (FGM) kháng nhiệt hoặc vật liệu áp điện thông minh (PEM) có khả năng cảm biến và điều khiển chủ động dao động.

Khoảng trống nghiên cứu (research gap) được xác định chính xác từ việc tổng thuật hơn 120 công trình quốc tế uy tín: (1) Thiếu vắng các mô hình giải tích dạng kín (closed-form solutions) dự đoán chính xác phản ứng tĩnh và động lực học của dầm sandwich đa lớp chứa lõi FGP-GPLRC kết hợp lớp bề mặt FGM/PEM chịu đồng thời trường lực cơ học và điện trường; (2) Chưa có khung tính toán không lưới (Mesh-Free Method - MFM) 1D thống nhất áp dụng kỹ thuật nội suy $C^1$-Hermite cho các điều kiện biên phức tạp; (3) Hiện tượng mất ổn định thuật toán điều khiển dao động chủ động (active vibration control instability) trên kết cấu nano-composite phi đối xứng chưa từng được nhận diện và xử lý triệt để; (4) Thiếu sự kiểm chứng thực nghiệm độc lập giữa lý thuyết phân tích bậc cao và vật liệu bê tông rỗng nhẹ ứng dụng thực tế.

Để giải quyết trọn vẹn khoảng trống trên, 4 câu hỏi nghiên cứu ($RQ$) và 4 giả thuyết khoa học ($H$) được định hình chặt chẽ:

  • $RQ_1$: Quy luật phân bố độ rỗng và hướng định vị của phiến nano GPL ảnh hưởng như thế nào đến độ cứng tổng thể, trạng thái ứng suất và phổ tần số dao động riêng của dầm sandwich?
    $\rightarrow \mathbf{H_1:}$ Sự kết hợp giữa phân bố độ rỗng đối xứng (Poro-S) và phân bố GPL đối xứng (GPL-S) sẽ tối ưu hóa mô-men kháng uốn, bù đắp hoàn toàn sự suy giảm độ cứng do các lỗ rỗng gây ra mà vẫn bảo toàn tỷ trọng siêu nhẹ.
  • $RQ_2$: Lời giải giải tích Navier trên nền tảng lý thuyết biến dạng cắt bậc cao (HSDT) có thể biểu diễn trọn vẹn tương tác cơ - điện (electro-mechanical coupling) và hiệu ứng điều khiển chủ động hay không?
    $\rightarrow \mathbf{H_2:}$ Mô hình giải tích HSDT kết hợp nguyên lý biến phân Hamilton sẽ thiết lập được nghiệm giải tích dạng kín chính xác cao mà không cần đến hệ số hiệu chỉnh cắt (shear correction factor).
  • $RQ_3$: Phương pháp không lưới MFM với các kỹ thuật nội suy PIM, RPIM, MLS có ưu thế vượt trội như thế nào so với FEM truyền thống trong việc phân tích các điều kiện biên bất kỳ?
    $\rightarrow \mathbf{H_3:}$ Kỹ thuật nội suy $C^1$-Hermite trong không gian 1D sẽ triệt tiêu hoàn toàn hiện tượng khóa cắt (shear locking) và bảo đảm tính liên tục bậc cao của các hàm chuyển vị.
  • $RQ_4$: Cơ chế điều khiển phản hồi vận tốc (velocity feedback control) qua lớp phủ áp điện PEM có bị giới hạn bởi tính bất đối xứng vật liệu hay không, và độ tin cậy thực tế của mô hình lý thuyết khi đối chứng với vật liệu bê tông xốp EPS thực tế đạt mức độ nào?
    $\rightarrow \mathbf{H_4:}$ Tính bất đối xứng vật liệu (GPL-A) sẽ kích hoạt sự ghép nối uốn - nén gây mất ổn định điều khiển khi áp dụng cảm biến/kích hoạt hai mặt đối xứng; việc tái cấu trúc cặp sensor/actuator về một phía sẽ khôi phục hoàn toàn tính ổn định động học.

Phạm vi nghiên cứu bao quát từ mô hình hóa vi cơ học (micromechanics) với lý thuyết Halpin-Tsai và quy tắc trộn mở rộng, lý thuyết biến dạng cắt bậc ba của Reddy, phương pháp tích phân thời gian Newmark gia tốc không đổi ($\beta = 1/4, \gamma = 1/2$), cho đến thực nghiệm trên dầm sandwich bê tông hạt xốp EPS tấm Cemboard SCG do Nucewall sản xuất. Nghiên cứu mang tính bước ngoặt khi định lượng hóa tác động tăng cường: chỉ với $1%$ tỷ lệ khối lượng GPL gia cường đã giúp tần số dao động riêng tăng vọt trên $41%$ và khả năng chịu tải tới hạn tăng hơn $85%$.

Literature Review và Positioning

Phân tích tổng quan y văn quốc tế cho thấy sự tiến hóa của 4 dòng chảy học thuật chính yếu định hình nên lĩnh vực nghiên cứu kết cấu sandwich tiên tiến:

timeline
    title Tiến Trình Phát Triển Lý Thuyết & Định Vị Học Thuật Của Luận Án
    1984 : Khởi xướng vật liệu FGM tại Nhật Bản (Koizumi, Niino)
    2000 : Lý thuyết vi cơ học Halpin-Tsai & Quy tắc trộn FGP (Ashby et al.)
    2017 : Đột phá gia cường GPL trong dầm rỗng FGP (Kitipornchai et al.)
    2020 : Mô phỏng bán giải tích Chebyshev-Ritz & HSDT (Wang et al., Xin & Kiani)
    2025 : Luận Án Đỗ Minh Đức : Tích hợp Giải tích Navier + MFM C1-Hermite + Active PEM Stability + Thực nghiệm EPS
  1. Dòng nghiên cứu vật liệu FGM và vật liệu rỗng FGP: Kể từ phát minh đột phá năm 1984 của các nhà khoa học Nhật Bản trong dự án hàng không vũ trụ nhằm chế tạo vật liệu chịu gradient nhiệt $1000\text{ K}$ trên tiết diện mỏng dưới $10\text{ mm}$ (Saleh et al. [5]), FGM đã giải quyết triệt để sự tập trung ứng suất. Sayyad và Ghugal [8] tổng hợp toàn diện các quy luật biến thiên hàm lũy thừa (power-law), hàm mũ (exponential) và hàm sigmoid đối xứng/phản xứng. Đối với vật liệu rỗng, Ashby et al. [10] và Chen et al. [13] đã xác lập quan hệ phi tuyến giữa khối lượng riêng và mô-đun đàn hồi theo cấu trúc lỗ rỗng mở (open-cell: $E(z)/E_0 = [\rho(z)/\rho_0]^2$) và lỗ rỗng kín (closed-cell: $E(z)/E_0 = [\rho(z)/\rho_0]^3$).
  2. Dòng nghiên cứu gia cường hạt nano carbon và graphene: Kitipornchai et al. [19] mở ra bước đột phá khi chứng minh việc đưa phiến phẳng graphene (GPL) vào ma trận xốp tạo ra sự gia tăng độ cứng vượt bậc. Các nghiên cứu tiếp nối của Wang et al. [39, 43] và Zheng et al. [46] áp dụng mô hình Halpin-Tsai để đánh giá đáp ứng động phi tuyến của dầm composite gia cường GPL chịu tải trọng va chạm và khối lượng di động.
  3. Dòng nghiên cứu phương pháp giải tích, bán giải tích và phương pháp số: 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 chịu tải phân bố. Tiếp đó, Xin và Kiani [35] dùng lời giải Navier và lý thuyết quasi-3D khảo sát tần số dao động của dầm sandwich lõi bọt kim loại trên nền đàn hồi. Nhằm khắc phục hạn chế về điều kiện biên của phương pháp giải tích, các phương pháp bán giải tích như Chebyshev-Ritz (Wang et al. [39]), Gram-Schmidt-Ritz (Srikarun et al. [41], Songsuwan et al. [44]), Differential Transform Method - DTM (Derikvand et al. [42]) và phương pháp số không lưới Mesh-Free (Amirani et al. [48], Gu & Liu [50], Mesmoudi et al. [58, 59]) được khai phá mạnh mẽ.
  4. Dòng nghiên cứu kết cấu thông minh tích hợp vật liệu áp điện (PEM): Bắt nguồn từ khám phá của Pierre Curie (1880) và sự phát triển của Preumont [20], Bamdad et al. [34] đã tích hợp các lớp PEM để khảo sát dao động và ổn định từ - nhiệt - cơ.
Tiêu chí Khoa học Nghiên cứu của Xin & Kiani (2020) [35] Nghiên cứu của Wang et al. (2021) [39] Luận án Tiến sĩ Đỗ Minh Đức (2025)
Mô hình cấu trúc cốt lõi Lõi bọt kim loại FGP đơn thuần, mặt đẳng hướng Lõi FGP, 2 lớp bề mặt FGM chịu tải di động Lõi FGP-GPLRC đa pha, 2 bề mặt FGM/PEM đa chức năng
Lý thuyết dầm sử dụng Quasi-3D Beam Theory HSDT biến dạng cắt mới + Chebyshev-Ritz Reddy HSDT bậc 3 thống nhất + Cặp liên hệ áp điện $e_{31}, p_{33}, d_{31}$
Phương pháp tính toán Giải tích Navier (chỉ xét biên SS) Bán giải tích Chebyshev-Ritz Đa phương pháp: Giải tích Navier + MFM ($C^1$-Hermite PIM, RPIM, MLS)
Phân tích điều khiển chủ động Không xem xét Không xem xét Kiểm soát dao động chủ động, phát hiện và xử lý mất ổn định $G_v$
Kiểm chứng thực nghiệm Không có (chỉ mô phỏng số) Không có (chỉ so sánh bán giải tích) Thực nghiệm toàn diện: Đo võng LVDT, gia tốc kế, FFT trên mẫu dầm EPS

Trong nước, các nhóm nghiên cứu mạnh của Đặ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], Bùi Quốc Tính et al. [77], Vũ Văn Tân et al. [79, 80], và Thái Hoàng Chiến et al. [81, 82] đã đạt được nhiều thành tựu lớn trên kết cấu tấm/vỏ FGM/FGP. Tuy nhiên, luận án của NCS. Đỗ Minh Đức đã định vị chính xác vị thế tiên phong bằng việc thiết lập một cầu nối hoàn chỉnh khép kín: Lý thuyết vi cơ học nano $\rightarrow$ Lời giải giải tích dạng kín Navier $\rightarrow$ Thuật toán không lưới MFM $C^1$-Hermite cho biên bất kỳ $\rightarrow$ Thuật toán điều khiển chủ động ổn định hóa $\rightarrow$ Thực nghiệm kiểm chứng trên kết cấu xây dựng thực tế.

Đó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 biên giới lý thuyết của cơ học vật liệu composite và kết cấu thông minh qua 3 đóng góp mang tính nền tảng:

  1. Mở rộng lý thuyết dầm biến dạng cắt bậc cao Reddy (HSDT) cho môi trường liên tục đa pha FGP-GPLRC/PEM:
    Trường chuyển vị được thiết lập chặt chẽ không đòi hỏi hệ số hiệu chỉnh cắt: $$\begin{cases} u(x,z,t) = u_0(x,t) + z \theta(x,t) - c_1 z^3 \left( \theta(x,t) + \dfrac{\partial w_0(x,t)}{\partial x} \right) \ w(x,z,t) = w_0(x,t) \end{cases}$$ trong đó $c_1 = 4/(3h^2)$, $u_0, w_0$ là chuyển vị dọc trục và độ võng tại mặt trung bình ($z=0$), $\theta$ là góc xoay của tiết diện ngang.

  2. Thiết lập mô hình giải tích dự đoán mô-đun đàn hồi $E(z)$ và khối lượng riêng $\rho(z)$ biến thiên:
    Theo quy luật phân bố độ rỗng liên tục: $$E(z) = E_0 \left[1 - e_0 \alpha(z)\right], \quad \rho(z) = \rho_0 \left[1 - e_p \alpha(z)\right]$$ với $e_0 = 1 - E_{\min}/E_0$ là hệ số độ rỗng, $e_p = 1 - \sqrt{1 - e_0}$ (cho cấu trúc lỗ rỗng mở) hoặc $e_p = 1 - \sqrt[3]{1 - e_0}$ (cho lỗ rỗng kín). Mô-đun đàn hồi của vật liệu không chứa lỗ rỗng $E_0$ được tính toán chính xác thông qua mô hình vi cơ học Halpin-Tsai: $$E_0 = \dfrac{3}{8}\left(\dfrac{1 + \xi_L \eta_L V_{\text{GPL}}}{1 - \eta_L V_{\text{GPL}}}\right) E_{\text{mf}} + \dfrac{5}{8}\left(\dfrac{1 + \xi_W \eta_W V_{\text{GPL}}}{1 - \eta_W V_{\text{GPL}}}\right) E_{\text{mf}}$$ với các hệ số hình học nano $\xi_L = 2(l_{\text{GPL}}/t_{\text{GPL}})$, $\xi_W = 2(w_{\text{GPL}}/t_{\text{GPL}})$.

  3. Phát hiện và chứng minh hiện tượng mất ổn định điều khiển dao động do bất đối xứng vật liệu:
    Luận án chỉ ra rằng khi dầm có phân bố GPL bất đối xứng (GPL-A) chịu điều khiển hồi tiếp vận tốc $V_a(x,t) = -G_v \dfrac{\partial w_0}{\partial t}$, ma trận cản chủ động bổ sung $C_a$ làm mất tính đối xứng của hệ phương trình trạng thái khi đặt cặp sensor/actuator ở hai mặt đối xứng, dẫn đến các nghiệm có phần thực dương tại độ lợi $G_v = 5 \times 10^{-5}$. Luận án đề xuất giải pháp thiết kế dời toàn bộ cặp sensor/actuator về một phía mặt dầm chủ, chứng minh bằng toán học sự khôi phục hoàn toàn tính ổn định tiệm cận của hệ thống.

classDiagram
    class AnalyticalFramework {
        +HamiltonPrinciple()
        +ReddyHSDT()
        +NavierSolution()
        +NewmarkIntegration()
    }
    class Micromechanics {
        +HalpinTsaiModel()
        +RuleOfMixtures()
        +PorosityDistributions(Poro-U, Poro-S, Poro-A)
        +GPLDistributions(GPL-U, GPL-S, GPL-A)
    }
    class SmartControl {
        +PiezoelectricConstitutive()
        +ActiveVelocityFeedback(Gv)
        +AsymmetricInstabilityRemedy()
    }
    class MeshFreeMethod {
        +PointInterpolation(PIM)
        +RadialBasisFunction(RPIM)
        +MovingLeastSquares(MLS)
        +HermiteC1Formulation()
    }
    AnalyticalFramework --> Micromechanics : Integrates
    AnalyticalFramework --> SmartControl : Couples
    AnalyticalFramework <|-- MeshFreeMethod : Generalizes Boundary Conditions

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

Khung phân tích của luận án tích hợp đồng thời 4 trụ cột lý thuyết: (1) Lý thuyết biến dạng cắt bậc cao Reddy HSDT; (2) Lý thuyết vi cơ học Halpin-Tsai; (3) Lý thuyết điện động học áp điện với các phương trình trạng thái áp điện trực tiếp: $$\sigma_x = Q_{11} \varepsilon_x - e_{31} E_z, \quad D_z = e_{31} \varepsilon_x + p_{33} E_z$$ và (4) Nguyên lý biến phân Hamilton toàn diện: $$\int_{t_1}^{t_2} \left( \delta K - \delta U + \delta V \right) dt = 0$$ trong đó $\delta K$ là biến phân động năng, $\delta U$ là biến phân thế năng biến dạng và năng lượng điện trường, $\delta V$ là biến phân công của ngoại lực cơ học phân bố $q_0(x,t)$ và điện áp kích hoạt $V_0(x,t)$.

Hệ phương trình vi phân chuyển động liên kết cơ - điện tổng quát thu được có dạng: $$\mathbf{M} \ddot{\mathbf{d}}(t) + (\mathbf{C} + \mathbf{C}_a) \dot{\mathbf{d}}(t) + \mathbf{K} \mathbf{d}(t) = \mathbf{F}_m(t) + \mathbf{F}_e(t)$$ Khung phân tích xác định rõ các điều kiện biên: Tựa đơn giản hai đầu (SS), Ngàm hai đầu (CC), Đầu ngàm - đầu khớp (CS), Đầu ngàm - đầu tự do (CF), bảo đảm tính chuẩn mực tuyệt đối trong cơ học tính toán.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt quan điểm nhận thức luận thực chứng (Positivism) kết hợp với chủ nghĩa thực dụng tính toán (Computational Pragmatism). Thiết kế nghiên cứu là sự kết hợp đa phương pháp (tri-partite methodology) gồm 3 nhánh tương hỗ: Giải tích dạng kín $\rightarrow$ Phương pháp không lưới MFM $\rightarrow$ Thực nghiệm vật lý.

graph LR
    subgraph Theoretical Branch
        T1["Reddy HSDT & Hamilton"] --> T2["Lời Giải Navier (Closed-Form)"]
    end
    
    subgraph Numerical Branch
        N1["Rời rạc hóa nút (No Mesh)"] --> N2["Nội suy C1-Hermite: PIM, RPIM, MLS"]
        N2 --> N3["Tích phân Gauss Background Cells"]
        N3 --> N4["Newmark Solver"]
    end
    
    subgraph Experimental Branch
        E1["Mẫu dầm Nucewall (EPS + Cemboard)"] --> E2["Hệ thống LVDT & Gia Tốc Kế"]
        E2 --> E3["Kyowa DCS-100A & DAS-200A FFT"]
    end
    
    T2 <-->|Cross Validation| N4
    N4 <-->|Empirical Benchmark| E3

Quy trình nghiên cứu rigorous

Quy trình số hóa sử dụng phương pháp không lưới (Mesh-Free Method - MFM) 1D độc đáo với 3 kỹ thuật xây dựng hàm dạng (shape functions):

  1. Phương pháp nội suy điểm đa thức (PIM): Sử dụng đa thức cơ sở $\mathbf{p}^T(x) = [1, x, x^2, \dots]$.
  2. Phương pháp nội suy hàm bán kính cơ sở (RPIM): Kết hợp hàm bán kính đa bậc hai nhiều nhánh (Multiquadric - MQ): $R_i(x) = (r_i^2 + C^2)^{k_{sc}}$ với các tham số hình dạng tối ưu được khảo sát hội tụ ($C = 0.1\sim 0.5, k_{sc} = 0.5\sim 1.5$).
  3. Phương pháp bình phương cực tiểu động (MLS): Dùng hàm trọng số spline bậc 4 trên miền ảnh hưởng (influence domain) chứa các nút $N_s$.

Kỹ thuật nội suy $C^1$-Hermite duy nhất cho phép xấp xỉ đồng thời chuyển vị cơ bản và đạo hàm bậc nhất tại các nút, đáp ứng hoàn hảo tính khả vi liên tục của trường chuyển vị mà không cần tăng số bậc tự do cục bộ. Quá trình tích phân số trên các ô lưới nền (background cells) được tối ưu hóa với số điểm Gauss khảo sát từ 2 đến 8 điểm, đảm bảo triệt tiêu năng lượng ảo và sai số tích phân.

Quy trình thực nghiệm được thiết lập bài bản tại phòng thí nghiệm chuyên dụng:

  • Đối tượng: Mẫu dầm sandwich kích thước $L = 1200\text{ mm}$, $b = 100\text{ mm}$, $h = 90\text{ mm}$, lớp lõi bê tông hạt xốp EPS dày $h_c = 70\text{ mm}$ ($\rho = 850\text{ kg/m}^3, E = 1.25\text{ GPa}$), hai lớp bề mặt tấm Cemboard SCG Thái Lan dày $h_f = 10\text{ mm}$ ($\rho = 1350\text{ kg/m}^3, E = 4.5\text{ GPa}$).
  • Thiết bị thu nhận dữ liệu: Cảm biến chuyển vị điện tử vi sai LVDT, gia tốc kế phân giải cao, hệ thống đo động đa kênh Kyowa DCS-100A, phần mềm xử lý tín hiệu và biến đổi Fourier nhanh DAS-200A (Fast Fourier Transform - FFT).

Data và phân tích

Toàn bộ thuật toán giải tích Navier, phần mềm MFM nội suy PIM/RPIM/MLS và thuật giải tích phân thời gian Newmark được lập trình độc lập trên nền tảng MATLAB (mã nguồn chi tiết tại Phụ lục B, C, D của luận án). Khảo sát độ hội tụ được tiến hành nghiêm ngặt: số số hạng chuỗi Fourier $n$ trong nghiệm Navier tăng từ $1$ đến $19$ (đạt độ ổn định 6 chữ số thập phân tại $n \ge 9$); số lượng nút MFM phân bố từ $N = 11, 21, 41$ đến $81$ nút (sai số tương đối dưới $0.05%$).

Kiểm chứng chéo (robustness checks) giữa nghiệm Navier, phần tử hữu hạn FEM tài liệu tham khảo và phương pháp MFM cho thấy sự trùng khít tuyệt đối với sai số cực đại không vượt quá $0.42%$. Kết quả đo dao động thực nghiệm qua phân tích phổ FFT xác định tần số dao động riêng thứ nhất $f_1 = 34.21\text{ Hz}$, hoàn toàn tương thích với dự đoán lý thuyết ($f_1^{\text{th}} = 35.56\text{ Hz}$, sai số chỉ $3.8%$).

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

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

  1. Hiệu ứng tối ưu hóa cấu trúc của phân bố Poro-S và GPL-S:
    Kết quả số chứng minh cấu hình độ rỗng đối xứng (Poro-S: lỗ rỗng tập trung nhiều nhất ở trục trung hòa, đặc chắc ở biên) kết hợp với phiến graphene phân bố tập trung ở hai mép ngoài (GPL-S) tạo ra mô-men kháng uốn cực đại. Độ võng không thứ nguyên giữa nhịp $\bar{w}$ giảm tới $48.6%$ so với trường hợp phân bố đều (Poro-U, GPL-U).
graph TD
    subgraph Cơ Chế Tối Ưu Hóa Độ Cứng
        P1["Phân bố Poro-S: Lõi rỗng ở tâm, đặc ở biên"] --> OPT["TỐI ƯU HÓA MÔ-MEN KHÁNG UỐN"]
        P2["Phân bố GPL-S: Nano GPL tập trung ở mép ngoài"] --> OPT
        OPT --> RES1["Giảm độ võng tĩnh: 48.6%"]
        OPT --> RES2["Tăng tần số cơ bản: >41%"]
        OPT --> RES3["Tiết kiệm khối lượng: 30-50%"]
    end
  1. Khả năng bù trừ độ cứng siêu việt của hạt nano GPL:
    Trong điều kiện hệ số độ rỗng tăng cao đến $e_0 = 0.6$ (làm giảm mô-đun đàn hồi hiệu dụng tới $60%$), việc chỉ cần bổ sung $1.0\text{ wt}%$ nano GPL đã làm tăng tần số dao động cơ bản không thứ nguyên $\bar{\omega}_1$ lên hơn $41%$, bù đắp hoàn toàn sự suy giảm độ cứng cơ học do độ rỗng tạo ra.
  2. Hiệu ứng triệt tiêu biến dạng nhờ ghép nối cơ - điện áp điện:
    Khi đặt cặp điện áp mặt trái dấu $V_0 = \pm 100\text{ V}$ lên các lớp PEM, điện trường $E_z$ sinh ra biến dạng ngược chiều với biến dạng uốn do tải trọng ngoài $q_0 = -100\text{ N/m}^2$, cho phép chủ động điều chỉnh độ võng tĩnh của dầm về gần bằng $0$.
  3. Hiện tượng phách điều hòa và cộng hưởng bị triệt tiêu bằng điều khiển hồi tiếp vận tốc:
    Với hệ số khuếch đại $G_v = 1 \times 10^{-4}$, biên độ dao động quá độ dưới tải trọng xung chữ nhật và xung hình sin bị dập tắt chỉ sau chưa đầy $0.15\text{ giây}$. Hiện tượng phách (beating) và cộng hưởng động lực học bị loại trừ hoàn toàn.
  4. Phát hiện hiện tượng phân kỳ mất ổn định động học:
    Khi cấu hình dầm có lớp lõi phân bố bất đối xứng GPL-A, việc điều khiển đối xứng tạo ra mô-men kích động ngược pha dẫn đến dao động phân kỳ vô hạn ($G_v = 5 \times 10^{-5}$); luận án chứng minh giải pháp bố trí cụm cảm biến/kích hoạt một phía triệt tiêu hoàn toàn sự mất ổn định này.

Implications đa chiều

  • Về mặt học thuật: Luận án đặt nền móng lý thuyết cho việc phân tích các kết cấu thông minh đa pha có gradient cơ tính phức tạp, cung cấp bộ dữ liệu benchmark giải tích chuẩn xác để cộng đồng cơ học tính toán kiểm chứng các phần mềm thương mại.
  • Về mặt phương pháp: Kỹ thuật không lưới $C^1$-Hermite MFM 1D phát triển trong luận án là công cụ tính toán mạnh mẽ, loại bỏ hoàn toàn chi phí tạo lưới phức tạp của FEM, mở đường cho việc áp dụng trên các kết cấu tấm và vỏ mỏng 2D/3D.
  • Về mặt kỹ thuật công nghệ: Mở ra giải pháp chế tạo các cấu kiện sàn, vách ngăn siêu nhẹ, panel cách âm, vỏ tàu cao tốc, thân máy bay không người lái (UAV) có khả năng tự triệt tiêu rung chấn và chịu xung kích vượt trội.
  • Về mặt chính sách và quy chuẩn: Cung cấp luận cứ khoa học thực nghiệm vững chắc để các cơ quan quản lý (Bộ Xây dựng, Viện Khoa học Công nghệ Xây dựng IBST) xây dựng tiêu chuẩn kỹ thuật thiết kế kết cấu composite bê tông nhẹ và panel sandwich thế hệ mới tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn chỉ rõ 4 giới hạn nghiên cứu:

  1. Mô hình dầm 1D: Chưa khảo sát đầy đủ trường ứng suất 3D phức tạp tại các góc mép tự do và hiện tượng tập trung ứng suất cục bộ dưới tải trọng tập trung cao.
  2. Giả thiết liên kết lớp lý tưởng: Coi sự bám dính giữa lớp bề mặt và lớp lõi là hoàn hảo, chưa mô hình hóa phi tuyến trượt dính (interlayer slip) hoặc phá hủy tách lớp theo thời gian.
  3. Môi trường đẳng nhiệt: Chưa xét đến tương tác nhiệt - điện - cơ đồng thời trong điều kiện cháy hoặc gradient nhiệt cực hạn.
  4. Vật liệu đàn hồi tuyến tính: Chưa xét đến biến dạng dẻo hoặc tổn thương tích lũy của hạt xốp EPS khi chịu tải trọng va đập lặp lại chu kỳ cao.

Chương trình nghiên cứu 5-10 năm tới mở ra 4 hướng phát triển mũi nhọn:

  • Mở rộng mô hình lý thuyết và MFM sang kết cấu tấm, vỏ sandwich FGP-GPLRC đa chiều (2D/3D plates and shells).
  • Tích hợp mô hình cơ học phá hủy liên tục (Continuum Damage Mechanics) và phần tử dính kết (Cohesive Zone Model - CZM) để mô phỏng sự tách lớp.
  • Khảo sát ứng xử nhiệt - cơ phi tuyến hình học (geometrically non-linear thermo-mechanical behavior) dưới điều kiện biên ngẫu nhiên.
  • Nghiên cứu công nghệ in 3D vật liệu FGP sinh học (PLA lattice) kết hợp cảm biến sợi quang (FBG) để giám sát sức khỏe kết cấu theo thời gian thực (Structural Health Monitoring - SHM).

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

mindmap
  root((Tác Động Của Luận Án))
    Học Thuật Quốc Tế
      Công bố Q1/Scopus
      Chỉ số trích dẫn cao
      Bộ số liệu Benchmark giải tích
    Công Nghiệp & Xây Dựng
      Panel tường nhẹ EPS Nucewall
      Vỏ tàu cao tốc & Radome hàng không
      Tiết kiệm 30-50% chi phí kết cấu
    Môi Trường & Xã Hội
      Giảm phát thải Carbon xi măng
      Cách âm cách nhiệt giảm tải HVAC
      Công trình xanh bền vững
    Quy Chuẩn & Đào Tạo
      Giáo trình Cơ học tính toán cao học
      Tiêu chuẩn thiết kế TCVN vật liệu mới
  • Tác động học thuật quốc tế: Các công trình công bố từ luận án trên các tạp chí chuyên ngành quốc tế uy tín (ISI/Scopus Q1/Q2) cung cấp giải pháp giải tích và không lưới tiên phong, ước tính thu hút sự quan tâm và trích dẫn lớn từ cộng đồng cơ học tính toán toàn cầu.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp giải pháp tấm sandwich lõi bê tông nhẹ EPS gia cường tại Công ty Cổ phần Tường nhẹ NUCEWALL, giúp tối ưu hóa chiều dày vách ngăn, giảm tải trọng bản thân công trình từ $30 - 50%$, rút ngắn tiến độ thi công lắp ghép.
  • Hiệu quả xã hội và môi trường: Sử dụng hạt xốp tái sinh và vật liệu nhẹ giúp giảm thiểu tiêu thụ clinker xi măng, giảm phát thải khí nhà kính $\text{CO}_2$, nâng cao hiệu quả tiết kiệm năng lượng điều hòa không khí nhờ hệ số truyền nhiệt thấp của lõi FGP.

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

Nhóm Đối Tượng Giá Trị Cốt Lõi Nhận Được Chỉ Số Lợi Ích Định Lượng
Nghiên cứu sinh & Nhà nghiên cứu trẻ Tiếp cận hệ phương trình vi phân biến phân, mã nguồn MATLAB MFM chuẩn tắc Tiết kiệm 6-12 tháng phát triển thuật toán mô phỏng kết cấu
Giáo sư & Chuyên gia Cơ học tính toán Bộ nghiệm giải tích Navier dạng kín và kỹ thuật nội suy $C^1$-Hermite Dữ liệu đối chuẩn (benchmark) tuyệt đối cho các mô hình số 2D/3D
Kỹ sư R&D Doanh nghiệp Xây dựng/Vật liệu Công thức vi cơ học Halpin-Tsai xác định tỷ lệ phối trộn GPL/EPS tối ưu Tăng $41%$ độ cứng với chỉ $1%$ chi phí phụ gia nano
Cơ quan Quản lý Quy chuẩn Xây dựng Căn cứ khoa học thực nghiệm để biên soạn tiêu chuẩn thiết kế panel nhẹ Hoàn thiện khung quy chuẩn TCVN cho vật liệu composite xây dựng

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 hệ phương trình vi phân chủ đạo liên kết cơ - điện hoàn chỉnh cho dầm sandwich lõi FGP-GPLRC tích hợp bề mặt PEM trên cơ sở lý thuyết dầm biến dạng cắt bậc cao Reddy (HSDT) mà không cần sử dụng hệ số hiệu chỉnh cắt. Luận án đã mở rộng lý thuyết dầm đàn hồi kinh điển sang môi trường vật liệu xốp biến thiên liên tục đa pha gia cường nano, giải quyết trọn vẹn sự phụ thuộc không gian của mô-đun Young $E(z)$, hệ số Poisson $\nu(z)$, khối lượng riêng $\rho(z)$ và ten-sơ áp điện $e_{31}(z)$.

2. Đột phá phương pháp luận của nghiên cứu so với các công trình quốc tế tiền nhiệm?

So với các nghiên cứu của Xin & Kiani [35] (chỉ giải tích Navier cho dầm tựa đơn giản SS) và Garg et al. [45] (sử dụng FEM truyền thống dễ bị khóa cắt và phụ thuộc lưới), luận án đã phát triển phương pháp không lưới MFM 1D với kỹ thuật nội suy $C^1$-Hermite thống nhất cho cả ba phương pháp nội suy PIM, RPIM và MLS. Kỹ thuật này cho phép mô phỏng chính xác tuyệt đối dầm sandwich đa lớp dưới mọi điều kiện biên phức tạp (CC, CS, CF, SS) với độ hội tụ trơn tru, triệt tiêu hoàn toàn sự méo mó lưới và khóa cắt số học.

3. Phát hiện bất ngờ nhất trong quá trình tính toán và cơ sở giải thích cơ học?

Phát hiện bất ngờ nhất là sự mất ổn định động học (dynamic instability / flutter-like divergence) của thuật toán điều khiển hồi tiếp vận tốc khi dầm có phân bố GPL bất đối xứng (GPL-A) với cấu hình cảm biến/kích hoạt đối xứng ở hai mặt. Cơ sở cơ học là do phân bố GPL bất đối xứng làm lệch trục trung hòa hình học khỏi trục trung hòa vật lý, tạo ra sự ghép nối uốn - dọc trục (axial-bending coupling). Điện áp điều khiển sinh ra lực dọc trục kích động kết hợp mô-men làm triệt tiêu độ cản tương đương, khiến hệ thống mất ổn định tại $G_v = 5 \times 10^{-5}$. Luận án đã giải quyết hiện tượng này bằng giải pháp tái định vị cặp sensor/actuator về cùng một phía của dầm.

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

Hoàn toàn đầy đủ và minh bạch. Toàn bộ mã nguồn MATLAB xây dựng hàm dạng PIM (Phụ lục B), hàm dạng RPIM (Phụ lục C), hàm dạng MLS (Phụ lục D), cùng quy trình tích phân Newmark (Phụ lục A) và toàn bộ thông số hình học, cơ lý tính của vật liệu hạt xốp EPS, tấm Cemboard SCG, cảm biến LVDT, gia tốc kế được công khai chi tiết trong văn bản luận án.

5. Lộ trình phát triển học thuật 10 năm tiếp theo từ kết quả luận án?

Lộ trình 10 năm tập trung vào 3 giai đoạn chiến lược:

  • Năm 1 - 3: Nâng cấp thuật toán MFM $C^1$-Hermite sang phân tích tấm và vỏ mỏng thoải (shallow shells) FGP-GPLRC 2D chịu tải trọng nhiệt - cơ phi tuyến.
  • Năm 4 - 6: Phát triển mô hình học máy (Physics-Informed Neural Networks - PINNs) kết hợp thuật toán tối ưu hóa đa mục tiêu (Multi-objective Genetic Algorithm) để tự động thiết kế gradient độ rỗng và mật độ GPL theo yêu cầu tải trọng thực tế.
  • Năm 7 - 10: Chuyển giao công nghệ sản xuất cấu kiện thông minh tự giám sát và dập tắt rung chấn tích hợp cảm biến áp điện in 3D quy mô công nghiệp trong các dự án hạ tầng giao thông tốc độ cao và nhà siêu cao tầng.

Kết luận

  1. Xây dựng thành công hệ thống mô hình giải tích dạng kín Navier trên cơ sở lý thuyết biến dạng cắt bậc cao Reddy HSDT, mô phỏng chính xác tuyệt đối ứng xử tĩnh, dao động tự do, dao động cưỡng bức và đáp ứng cơ - điện của dầm sandwich lõi FGP-GPLRC.
  2. Phát triển hoàn chỉnh công cụ tính toán số bằng phương pháp không lưới MFM ($C^1$-Hermite PIM, RPIM, MLS), chứng minh tính ưu việt vượt trội trong việc phân tích dầm sandwich với các điều kiện biên động học bất kỳ (SS, CC, CS, CF) mà không phụ thuộc vào lưới phần tử.
  3. Lần đầu tiên phát hiện và đề xuất giải pháp triệt tiêu hiện tượng mất ổn định thuật toán điều khiển dao động chủ động trên kết cấu nanocomposite có phân bố vật liệu bất đối xứng, hoàn thiện lý thuyết điều khiển kết cấu thông minh.
  4. Xác lập quy luật tối ưu hóa cấu trúc: Khẳng định cấu hình kết hợp Poro-S (lõi rỗng ở giữa) và GPL-S (nano graphene tập trung ở mép ngoài) đem lại độ cứng kháng uốn cực đại, cho phép giảm $30 - 50%$ trọng lượng công trình trong khi nâng cao trên $41%$ tần số dao động riêng cơ bản.
  5. Thực hiện thành công thực nghiệm kiểm chứng độc lập trên mẫu dầm sandwich lõi bê tông hạt xốp EPS tấm Cemboard SCG thực tế, xác nhận độ tin cậy vượt trội của mô hình lý thuyết với sai số dưới $3.8%$.
  6. Mở ra 3 nhánh nghiên cứu liên ngành đột phá: Cơ học tính toán không lưới cho kết cấu thông minh nano, công nghệ chế tạo cấu kiện bê tông siêu nhẹ cách âm cách nhiệt đa năng, và hệ thống kiểm soát dao động chủ động thích ứng cho các công trình tương lai.