Tổng quan về luận án

Nghiên cứu cơ học kết cấu hiện đại đang chứng kiến bước chuyển mình mạnh mẽ nhằm đáp ứng các yêu cầu khắt khe về tối ưu hóa trọng lượng, khả năng chịu lực vượt trội và độ bền nhiệt trong các ngành công nghiệp mũi nhọn như hàng không - vũ trụ, đóng tàu biển, quốc phòng và kỹ thuật xây dựng công trình biển. Kết cấu sandwich truyền thống gồm ba lớp với cấu hình hai lớp bề mặt mỏng chịu lực kéo/nén và lớp lõi dày có trọng lượng nhẹ chủ yếu chịu lực cắt đã chứng minh được hiệu quả chịu uốn và hấp thụ năng lượng cao. Tuy nhiên, các rào cản kỹ thuật cố hữu như sự tập trung ứng suất tại bề mặt tiếp giáp gây bong tách lớp (delamination) và giới hạn chịu nhiệt của vật liệu composite polymer thông thường đã đặt ra bài toán cấp thiết về việc phát triển các thế hệ vật liệu tiên tiến mới.

Luận án tiến sĩ kỹ thuật mang tên "Phân tích phi tuyến tĩnh và dao động của dầm sandwich FGP gia cường GPL" do Nghiên cứu sinh Hương Quý Trường thực hiện tại Trường Đại học Xây dựng Hà Nội (chuyên ngành Cơ kỹ thuật, mã số: 9520101; 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ú) đã tiên phong giải quyết trọn vẹn khoảng trống khoa học này. Luận án đặt trọng tâm vào việc thiết lập mô hình kết cấu dầm sandwich tiên tiến tích hợp hai lớp bề mặt bằng vật liệu có cơ tính biến thiên (Functionally Graded Material - FGM) hệ gốm - kim loại và lớp lõi xốp biến thiên (Functionally Graded Porous - FGP) dạng bọt kim loại được gia cường bằng các phiến nano graphene (Graphene Platelets - GPL).

+-------------------------------------------------------------------------+
| LỚP BỀ MẶT TRÊN: Vật liệu FGM (Gốm/Kim loại) - Biến thiên hàm mũ p      |
+-------------------------------------------------------------------------+
| LỚP LÕI TRUNG TÂM: Vật liệu rỗng FGP dạng bọt kim loại gia cường GPL   |
| (Mô hình tế bào kín GRF; Phân bố độ rỗng PD1/PD2/UPD; Gia cường A/B/C)  |
+-------------------------------------------------------------------------+
| LỚP BỀ MẶT DƯỚI: Vật liệu FGM (Kim loại/Gốm) - Biến thiên hàm mũ p      |
+-------------------------------------------------------------------------+
         |                                                 |
   [Biên: CC / CS / SS / CF]                     [Nền đàn hồi Pasternak: K0, J0]

Research Gaps cụ thể từ y văn

Qua tổng quan các công trình nghiên cứu trong và ngoài nước, tác giả xác định ba khoảng trống học thuật (research gaps) mang tính bản lề:

  1. Thiếu vắng mô hình kết cấu sandwich đa cấu tử tích hợp toàn diện: Mặc dù các nghiên cứu trước đây (như Vo và Nguyen, 2013, 2015; Chen và cs., 2016; Wang và cs., 2018) đã khảo sát dầm FGM thuần nhất hoặc dầm FGP đơn lẻ, chưa có công trình nào thiết lập một mô hình hoàn chỉnh kết hợp đồng thời hai lớp bề mặt FGM và lớp lõi xốp FGP gia cường GPL phân tán theo không gian.
  2. Khoảng trống về phân tích ứng xử phi tuyến hình học tổng thể: Phần lớn các nghiên cứu hiện hành chỉ dừng lại ở phân tích tuyến tính với giả thiết chuyển vị bé hoặc khảo sát đơn lẻ một chỉ tiêu cơ học riêng biệt mà chưa đánh giá đồng thời ứng xử phi tuyến tĩnh (quan hệ phi tuyến tải - độ võng, tải - mô men) và đáp ứng dao động phi tuyến (tần số dao động tự do phi tuyến phụ thuộc biên độ và đáp ứng chuyển vị cưỡng bức phi tuyến theo miền thời gian).
  3. Sự phân mảnh trong các mô hình lý thuyết dầm: Việc tính toán dầm sandwich nhiều lớp thường đòi hỏi các mô hình biến dạng cắt khác nhau, tạo nên sự phân tán trong các công thức giải tích và gây khó khăn khi so sánh định lượng độ chính xác giữa các lý thuyết biến dạng cắt bậc cao.

Câu hỏi nghiên cứu và Giả thuyết khoa học

Luận án được định hình xoay quanh 4 câu hỏi nghiên cứu ($RQ$) và 4 giả thuyết khoa học ($H$) tương ứng:

  • $RQ_1$: Quy luật phân bố lỗ rỗng và cơ chế phân tán phiến nano GPL trong lớp lõi FGP ảnh hưởng như thế nào đến các hằng số đàn hồi hiệu dụng và khả năng chịu tải của kết cấu dầm sandwich?
    • $H_1$: Sự kết hợp giữa quy luật phân bố lỗ rỗng đối xứng (PD1) và phân tán GPL tập trung tại vùng biên lớp lõi (Dạng A) sẽ tối đa hóa mô men quán tính uốn, bù đắp hoàn toàn sự suy giảm độ cứng do độ rỗng gây ra mà không làm gia tăng trọng lượng dầm.
  • $RQ_2$: Làm thế nào để xây dựng một khung phân tích bán giải tích hợp nhất có khả năng bao quát đồng thời các lý thuyết dầm từ bậc nhất đến bậc cao khi xét đến hiệu ứng phi tuyến hình học von Kármán?
    • $H_2$: Việc tích hợp trường chuyển vị tổng quát với phương pháp đa thức Pb-Ritz dựa trên nguyên lý biến phân Hamilton sẽ thiết lập được hệ phương trình đại số phi tuyến và phương trình vi phân chuyển động có độ hội tụ cao cho mọi điều kiện biên cơ học.
  • $RQ_3$: Mức độ ảnh hưởng của tương tác kết cấu - nền đàn hồi hai tham số Pasternak và tỷ số độ mảnh hình học $L/h$ đến ứng xử phi tuyến uốn và các thành phần ứng suất nội tại của dầm sandwich là gì?
    • $H_3$: Độ cứng chống trượt của nền Pasternak ($J_0$) sẽ chế ngự đáng kể biên độ võng phi tuyến và làm thay đổi mạnh mẽ sự phân bố ứng suất pháp $\sigma_x$, ứng suất tiếp $\tau_{xz}$ theo chiều cao tiết diện dầm.
  • $RQ_4$: Tính phi tuyến hình học chi phối như thế nào đến đặc tính dao động tự do và hiện tượng cộng hưởng động lực học của dầm sandwich khi chịu tải trọng điều hòa di động hoặc phân bố?
    • $H_4$: Hiệu ứng cứng hóa hình học (hardening spring behavior) sẽ làm tăng tỷ số tần số phi tuyến $\omega_{NL}/\omega_L$ theo biên độ dao động, đồng thời làm lệch pha và giới hạn biên độ đáp ứng chuyển vị cưỡng bức khi tần số kích thích tiến gần tần số dao động riêng.

Khung lý thuyết và Phạm vi nghiên cứu

Luận án xây dựng khung lý thuyết dựa trên sự hợp nhất giữa:

  • Mô hình đồng nhất hóa vi cơ học Halpin-Tsai cải tiến xác định mô đun đàn hồi $E_1$ của pha nanocomposite gia cường GPL;
  • Mô hình cấu trúc tế bào kín theo trường ngẫu nhiên chuẩn Gaussian Random Field (GRF) của Roberts & Garboczi (2000, 2002) mô tả các tính chất cơ học hiệu dụng của lớp bọt kim loại xốp FGP;
  • Mô hình hàm lũy thừa Voigt mô tả sự biến thiên liên tục của tỷ phần thể tích gốm/kim loại trong hai lớp bề mặt FGM;
  • Lý thuyết dầm biến dạng cắt tổng quát tích hợp biến dạng phi tuyến von Kármán và nguyên lý biến phân Hamilton.

Phạm vi và quy mô khảo sát: Luận án tiến hành phân tích toàn diện trên dầm sandwich với các tỷ lệ chiều dày lớp (1-1-1, 1-2-1, 1-4-1, 1-8-1), tỷ số độ mảnh $L/h \in [5, 20]$, hệ số lỗ rỗng $e_0 \in [0.1, 0.9]$, tỷ trọng khối lượng GPL $W_{GPL} \in [0%, 2%]$, chỉ số tỷ phần thể tích FGM $p \in [0, 5]$, chịu 4 dạng điều kiện biên kinh điển: Ngàm - Ngàm (CC), Khớp - Khớp (SS), Ngàm - Khớp (CS), Ngàm - Tự do (CF) và đặt trên nền đàn hồi Pasternak ($K_0 \in [0, 100], J_0 \in [0, 10]$).

Literature Review và Positioning

Tổng quan các dòng nghiên cứu chính

Lịch sử phát triển của cơ học kết cấu sandwich ghi nhận sự tiến hóa từ các vật liệu đẳng hướng cổ điển sang các hệ composite phân lớp chức năng cao cấp:

[1940s-1980s]               [1990s-2010s]                 [2010s-2020s]                  [Luận án (2023)]
Sandwich Cổ Điển  ------->  Sandwich FGM       ------->   Sandwich FGP/GPL   ------->   Sandwich FGM-FGP/GPL
(Gỗ, Balsa, Nhôm,           (Gốm/Kim loại,                (Bọt xốp kim loại,            (Mô hình 3 lớp hợp nhất,
 lõi xốp polymer)            chống bong tách,              gia cường nano CNT/GPL,       phi tuyến von Kármán,
 Frostig & Baruch (1990)     Vo & Nguyen (2013)            Chen (2016), Wang (2018)      Pb-Ritz & HSDT tổng quát)
  1. Dòng nghiên cứu kết cấu sandwich truyền thống và lõi mềm (Flexible/Soft Core): Các công trình kinh điển của Frostig & Baruch (1990, 1993) đã đặt nền móng cho việc phân tích dầm sandwich có lớp lõi mềm chịu uốn và dao động riêng dựa trên mô hình dầm Timoshenko giải tích. Salami và cs. (2014) cùng Shen và cs. (2013) tiếp tục mở rộng phân tích tĩnh và thực nghiệm uốn dầm sandwich lõi xốp theo lý thuyết biến dạng cắt bậc cao.
  2. Dòng nghiên cứu kết cấu sandwich FGM: Nhằm khắc phục hiện tượng tách lớp do sự bất liên tục về tính chất nhiệt - cơ học tại mặt phân giới, vật liệu FGM được đưa vào kết cấu. Vo và Nguyen (2013, 2015) đã ứng dụng lý thuyết tựa đàn hồi ba chiều (quasi-3D) và chuyển vị bậc cao dạng hyperbolic để khảo sát tĩnh, dao động và mất ổn định của dầm sandwich FGM. Zenkour và cs. (2014) khảo sát dầm sandwich bề mặt FGM nhớt đàn hồi trên nền đàn hồi Pasternak bằng hàm chuyển vị dạng sin. Sayyad và Avhad (2017) phát triển lời giải Navier cho dầm sandwich FGM theo lý thuyết trượt hyperbolic.
  3. Dòng nghiên cứu vật liệu rỗng FGP và gia cường hạt Nano (GPL/CNT): Sự kết hợp giữa cấu trúc bọt rỗng để giảm trọng lượng và gia cường hạt nano để bù độ cứng tạo ra bước đột phá mới. Kitipornchai và cs. (2017), Chen và cs. (2016) đã tiên phong công bố các nghiên cứu về dao động phi tuyến của dầm FGP đơn lớp. Wang và cs. (2018, 2020) nghiên cứu uốn phi tuyến dầm sandwich lõi bọt kim loại gia cường GPL bề mặt bằng phương pháp Chebyshev-Ritz.

Các tranh luận học thuật chính (Academic Debates)

Trong y văn cơ học kết cấu composite, hai luồng quan điểm đối lập lớn tồn tại:

  • Tranh luận 1: Lý thuyết đơn lớp tương đương (Equivalent Single Layer - ESL) đối đầu Lý thuyết đa lớp (Layerwise Theory - LW): Luồng quan điểm ủng hộ LW cho rằng cần thiết lập trường chuyển vị riêng biệt cho từng lớp để mô tả chính xác sự nhảy bậc của đạo hàm biến dạng tại mặt phân giới. Ngược lại, trường phái ESL chỉ ra rằng đối với vật liệu cơ tính biến thiên liên tục (FGM/FGP), mô hình ESL kết hợp các hàm biến dạng cắt bậc cao hoàn toàn thỏa mãn các điều kiện biên về ứng suất cắt triệt tiêu tại hai bề mặt tự do $\sigma_{xz}(\pm h/2) = 0$ mà không làm bùng nổ số lượng bậc tự do tính toán, giảm thiểu đáng kể chi phí bộ nhớ máy tính.
  • Tranh luận 2: Lý thuyết Timoshenko (FSDBT) với hệ số hiệu chỉnh cắt đối đầu Lý thuyết biến dạng cắt bậc cao (HSDT): FSDBT đòi hỏi hệ số hiệu chỉnh cắt $K_s$ (thường lấy $5/6$), tuy nhiên giá trị này phụ thuộc phức tạp vào tỷ lệ hình học, điều kiện biên và quy luật biến thiên cơ tính của vật liệu phân lớp. Việc sử dụng HSDT (Reddy, hàm lượng giác, hàm hyperbolic, hàm e-mũ) loại bỏ hoàn toàn sự phụ thuộc vào $K_s$, mang lại nghiệm ứng suất tiếp phân bố dạng parabol chính xác qua chiều dày tiết diện.

Định vị học thuật và Bước tiến so với các nghiên cứu quốc tế

So với các nghiên cứu quốc tế tiêu biểu như công trình của Frostig & Baruch (1990) (chỉ dừng lại ở lõi đẳng hướng tuyến tính), Chen và cs. (2016) (chỉ khảo sát dầm FGP một lớp theo Timoshenko) và Wang và cs. (2018) (khảo sát dầm sandwich lõi FGP bề mặt GPLRC bằng phương pháp Chebyshev-Ritz nhưng chưa tích hợp lớp vỏ FGM), luận án của Hương Quý Trường (2023) đã thiết lập một bước tiến vượt bậc:

  • Xây dựng mô hình dầm sandwich 3 lớp bất đối xứng hoàn chỉnh: Hai lớp biên FGM hàm lũy thừa và lớp lõi FGP gia cường GPL đa quy luật phân bố.
  • Khung giải thuật Pb-Ritz tổng quát hóa, cho phép chuyển đổi linh hoạt giữa 6 lý thuyết dầm khác nhau chỉ bằng việc thay đổi hàm định dạng $f(z)$, đồng thời giải quyết triệt để bài toán phi tuyến hình học von Kármán.

Đó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 lý thuyết đàn hồi phi tuyến và cơ học vật liệu composite phân lớp thông qua các đóng góp lý thuyết nền tảng:

                  NGUYÊN LÝ BIẾN PHÂN HAMILTON
                  \int_{t1}^{t2} (\delta T - \delta U + \delta W) dt = 0
                                |
             +------------------+------------------+
             |                                     |
    TRƯỜNG CHUYỂN VỊ TỔNG QUÁT             BIẾN DẠNG PHI TUYẾN VON KÁRMÁN
    u(x,z,t) = u0 - z*(dw0/dx) + f(z)*phi_x     epsilon_x = du/dx + 0.5*(dw/dx)^2
    w(x,z,t) = w0(x,t)                          gamma_xz = dw/dx + du/dz
             |                                     |
             +------------------+------------------+
                                |
                     KHUNG NĂNG LƯỢNG PB-RITZ
       Xấp xỉ chuyển vị: u0(x) = P(x)*q_u; w0(x) = P(x)*q_w; phi_x(x) = P(x)*q_phi
                                |
             +------------------+------------------+
             |                                     |
    TĨNH PHI TUYẾN (Newton-Raphson)        DAO ĐỘNG PHI TUYẾN (RK4 / Lặp trị riêng)
    ([K_L] + [K_NL(q)]) {q} = {F}          [M]{q''} + ([K_L] + [K_NL(q)]){q} = {F(t)}

  1. Thiết lập Khung động học dầm tổng quát (Unified Kinematic Framework): Trường chuyển vị của dầm được biểu diễn dưới dạng tổng quát: $$u(x, z, t) = u_0(x, t) - z \frac{\partial w_0(x, t)}{\partial x} + f(z) \varphi_x(x, t)$$ $$w(x, z, t) = w_0(x, t)$$ Trong đó, hàm phân bố biến dạng cắt $f(z)$ được mô đun hóa để đại diện chính xác cho:

    • Lý thuyết dầm Timoshenko (FSDBT): $f(z) = z$
    • Lý thuyết dầm bậc ba Reddy (PSDBT): $f(z) = z \left(1 - \frac{4z^2}{3h^2}\right)$
    • Lý thuyết dầm dạng hàm Sin (TSDBT): $f(z) = \frac{h}{\pi} \sin\left(\frac{\pi z}{h}\right)$
    • Lý thuyết dầm dạng hàm Sinh (HSDBT): $f(z) = h \sinh\left(\frac{z}{h}\right) - z \cosh\left(\frac{1}{2}\right)$
    • Lý thuyết dầm dạng hàm e-mũ (ESDBT): $f(z) = z \exp\left(-2\left(\frac{z}{h}\right)^2\right)$

  2. Mô hình hóa vi cơ học đa tỷ lệ (Multi-scale Micromechanical Constitutive Model):

    • Đối với hai lớp bề mặt FGM, các hằng số vật liệu hiệu dụng (Mô đun Young $E$, hệ số Poisson $\nu$, khối lượng riêng $\rho$, mô đun trượt $G$) tuân theo quy tắc Voigt:

      "Sử dụng mô hình Voigt [132], hai lớp bề mặt bằng vật liệu FGM với các hằng số vật liệu $H$ được giả thiết biến thiên liên tục theo tọa độ chiều cao dầm như sau: $H(z) = V_m H_m + V_c H_c$." (Trích dẫn Luận án, Chương 2, tr. 25). Quy luật biến thiên ở lớp bề mặt dưới: $$E_b(z) = E_m + (E_c - E_m)\left(\frac{z - z_1}{z_2 - z_1}\right)^p; \quad \rho_b(z) = \rho_m + (\rho_c - \rho_m)\left(\frac{z - z_1}{z_2 - z_1}\right)^p$$

    • Đối với lớp lõi FGP gia cường GPL, cấu trúc lỗ rỗng tế bào kín được mô hình hóa theo sơ đồ GRF:

      "Sử dụng mô hình vật liệu có cấu trúc tế bào kín phân bố theo quy luật trường ngẫu nhiên chuẩn (closed-cell cellular under Gaussian Random Field - GRF scheme) trong tài liệu [124], ta có: $\frac{E_{core}(z)}{E_1} = \left(\frac{\rho_{core}(z)}{\rho_1}\right)^{2.3} + 0.121\left(\frac{\rho_{core}(z)}{\rho_1}\right) / 1.121$." (Trích dẫn Luận án, Chương 2, tr. 28).

    • Mô đun đàn hồi $E_1$ của nền kim loại chứa GPL khi chưa có lỗ rỗng được xác định chính xác thông qua mô hình Halpin-Tsai: $$E_1 = \left[ \frac{3}{8}\left(\frac{1 + \xi_L \eta_L V_{GPL}}{1 - \eta_L V_{GPL}}\right) + \frac{5}{8}\left(\frac{1 + \xi_w \eta_w V_{GPL}}{1 - \eta_w V_{GPL}}\right) \right] E_m$$ Với các thông số hình học hạt nano GPL: $\xi_L = 2 L_{GPL}/t_{GPL}$, $\xi_w = 2 w_{GPL}/t_{GPL}$.

  3. Thiết lập Mối quan hệ Biến dạng - Ứng suất Phi tuyến: Tích hợp thành phần phi tuyến hình học von Kármán: $$\varepsilon_x = \frac{\partial u}{\partial x} + \frac{1}{2}\left(\frac{\partial w}{\partial x}\right)^2 = \frac{\partial u_0}{\partial x} - z \frac{\partial^2 w_0}{\partial x^2} + f(z)\frac{\partial \varphi_x}{\partial x} + \frac{1}{2}\left(\frac{\partial w_0}{\partial x}\right)^2$$ $$\gamma_{xz} = \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x} = f'(z) \varphi_x$$

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

Phương pháp năng lượng biến phân dựa trên nguyên lý Hamilton: $$\int_{t_1}^{t_2} \left( \delta T - \delta U_B - \delta U_F + \delta W \right) dt = 0$$ Trong đó: $T$ là động năng kết cấu kể đến quán tính quay, $U_B$ là thế năng biến dạng đàn hồi nội tại, $U_F$ là thế năng phản lực nền đàn hồi Pasternak ($U_F = \frac{1}{2} \int_0^L (K_W w_0^2 + K_P (\frac{\partial w_0}{\partial x})^2) dx$), và $W$ là công ngoại lực của tải trọng phân bố $P_0$.

Khung phân tích thiết lập ma trận phần tử giải tích đồng nhất, trong đó các phương trình vi phân chuyển động phi tuyến bậc hai liên kết được dẫn xuất chính xác mà không gặp phải hiện tượng khóa cắt (shear locking) nhờ việc xấp xỉ trường chuyển vị độc lập qua chuỗi đa thức Pb-Ritz.

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

Thiết kế nghiên cứu

  • Triết lý nghiên cứu (Research Philosophy): Nghiên cứu được xây dựng trên nền tảng chủ nghĩa thực chứng thực nghiệm số (Computational Positivism) kết hợp cơ học môi trường liên tục tất định.
  • Tiếp cận Bán giải tích (Semi-analytical Approach): Sử dụng phương pháp đa thức Ritz với hàm biên phụ thuộc (Pb-Ritz method - Polynomial boundary Ritz). Chuyển vị chưa biết được xấp xỉ bằng chuỗi đa thức nhân với hàm biên hình học thỏa mãn tiên nghiệm các điều kiện biên động học: $$u_0(x) = \phi_u(x) \sum_{i=1}^m c_{ui} P_i(x); \quad w_0(x) = \phi_w(x) \sum_{j=1}^m c_{wj} P_j(x); \quad \varphi_x(x) = \phi_\varphi(x) \sum_{k=1}^m c_{\varphi k} P_k(x)$$ Trong đó, các hàm biên $\phi(x) = (x/L)^\alpha (1 - x/L)^\beta$ với các chỉ số $\alpha, \beta \in {0, 1, 2}$ được gán tự động tương ứng với từng điều kiện liên kết (SS, CC, CS, CF).
+----------------------------------------------------------------------------------------------------+
|                                  QUY TRÌNH THUẬT TOÁN TÍNH TOÁN                                    |
+----------------------------------------------------------------------------------------------------+
| 1. KHỞI TẠO CƠ TÍNH: Voigt (FGM) + GRF & Halpin-Tsai (FGP + GPL) -> Ma trận độ cứng A, B, D, E, F, H|
+----------------------------------------------------------------------------------------------------+
                                                  |
+----------------------------------------------------------------------------------------------------+
| 2. RỜI RẠC HÓA PB-RITZ: Thiết lập các ma trận [K_L], [K_NL(q)], [M] và véc-tơ ngoại lực {F}       |
+----------------------------------------------------------------------------------------------------+
                                                  |
         +----------------------------------------+----------------------------------------+
         |                                        |                                        |
         v                                        v                                        v
+-----------------------+        +--------------------------------+        +-----------------------+
|  TĨNH PHI TUYẾN       |        |  DAO ĐỘNG TỰ DO PHI TUYẾN      |        |  ĐÁP ỨNG ĐỘNG CƯỠNG BỨC|
|  Thuật toán:          |        |  Thuật toán:                   |        |  Thuật toán:          |
|  Newton-Raphson       |        |  Lặp gần đúng dần trị riêng    |        |  Runge-Kutta bậc 4    |
|  Sai số: eps <= 1e-6  |        |  Xác định: omega_NL / omega_L  |        |  Bước thời gian dt    |
+-----------------------+        +--------------------------------+        +-----------------------+

Quy trình nghiên cứu Rigorous và Kiểm soát độ tin cậy

  1. Khảo sát tính hội tụ của nghiệm (Convergence Analysis): Số lượng số hạng đa thức $m$ trong phép khai triển chuỗi Pb-Ritz được khảo sát chi tiết từ $m = 4$ đến $m = 18$. Kết quả trong Luận án chứng minh rằng nghiệm độ võng không thứ nguyên $\bar{w}$, mô men uốn $\bar{M}$ và tần số dao động tự do $\bar{\omega}$ đạt mức độ hội tụ tuyệt đối với sai số biến thiên dưới $0.01%$ khi $m \ge 12$. Do đó, cấu hình $m = 14$ được chọn làm chuẩn cho toàn bộ các phép tính số.
  2. Quy trình giải bài toán Tĩnh phi tuyến: Sử dụng thuật toán lặp Newton-Raphson đa chiều với điều kiện dừng nghiêm ngặt: $$\frac{|\mathbf{q}^{(r+1)} - \mathbf{q}^{(r)}|}{|\mathbf{q}^{(r+1)}|} \le 10^{-6}$$
  3. Quy trình giải bài toán Dao động tự do phi tuyến: Áp dụng phương pháp lặp xấp xỉ liên tiếp xác định tần số dao động tự do phi tuyến $\omega_{NL}$ tại các biên độ dao động đỉnh $w_{max}/h \in [0, 2.0]$.
  4. Quy trình giải bài toán Đáp ứng động phi tuyến: Sử dụng thuật toán tích phân thời gian tường minh Runge-Kutta bậc 4 (RK4) với bước thời gian $\Delta t = 10^{-4}\text{ s}$ nhằm bắt trọn quỹ đạo pha (phase portraits) và đáp ứng miền thời gian ổn định.

Kiểm chứng độ tin cậy độc lập (Validation)

Độ chính xác của thuật toán và mã nguồn MATLAB tự xây dựng được kiểm chứng chéo thông qua 5 bài toán đối chuẩn quốc tế:

  • Kiểm chứng độ võng không thứ nguyên của dầm đẳng hướng chịu tải phân bố đều trên nền Pasternak so sánh với kết quả giải tích của Timoshenko (1921) và Reddy (1984) (sai số $< 0.05%$).
  • Kiểm chứng độ võng và ứng suất pháp $\bar{\sigma}x$, ứng suất tiếp $\bar{\tau}{xz}$ của dầm sandwich FGM so sánh với nghiệm giải tích 3D của Vo và cs. (2015) và Zenkour (2006).
  • Kiểm chứng tần số dao động tự do phi tuyến $\omega_{NL}/\omega_L$ với công bố của Chen và cs. (2016) cho dầm FGP (kết quả trùng khớp đến 4 chữ số thập phân).

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

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

+------------------------------------------------------------------------------------------------------+
| TÓM TẮT PHÁT HIỆN THỰC NGHIỆM SỐ ĐỘT PHÁ CỦA LUẬN ÁN                                                 |
+------------------------------------------------------------------------------------------------------+
| 1. HIỆU QUẢ GIA CƯỜNG GPL: 1.0% wt GPL (Dạng A) giúp GIẢM 58.4% ĐỘ VÕNG & TĂNG 84.2% TẦN SỐ TỰ DO.   |
| 2. TỐI ƯU HÓA LỖ RỖNG: Phân bố PD1 (đối xứng) vượt trội hoàn toàn PD2 (bất đối xứng) và UPD (đều).  |
| 3. CỨNG HÓA HÌNH HỌC PHI TUYẾN: Tỷ số omega_NL/omega_L > 1.45 khi w_max/h = 1.0 (Biên CC, L/h=20).   |
| 4. CHẾ NGỰ NỀN PASTERNAK: Nền cắt J0 = 10 triệt tiêu 72.3% dao động cộng hưởng biên độ lớn.         |
| 5. GIỚI HẠN LÝ THUYẾT: FSDBT sai lệch 8-12% ở dầm dày (L/h=5); HSDT mô tả hoàn hảo ứng suất tiếp.    |
+------------------------------------------------------------------------------------------------------+

  1. Hiệu ứng gia cường vượt trội của phiến nano GPL (GPL Reinforcement Synergy): Kết quả số chỉ ra rằng việc bổ sung chỉ $1.0%$ khối lượng GPL ($W_{GPL} = 1%$) vào lớp lõi xốp bọt kim loại với quy luật phân tán Dạng A (tập trung đối xứng tại hai mặt biên của lớp lõi) giúp làm giảm độ võng phi tuyến tĩnh lên tới $58.4%$ và tăng tần số dao động riêng cơ bản lên $84.2%$ so với kết cấu chưa gia cường ($W_{GPL} = 0%$). Điều này chứng minh cơ chế truyền ứng suất cắt cực kỳ hiệu quả của cấu trúc nano graphene dạng phiến 2D so với ống nano carbon 1D (CNT).

  2. Ảnh hưởng quyết định của Quy luật phân bố độ rỗng (Porosity Distribution Mechanics): Trong ba dạng phân bố độ rỗng của lớp lõi:

    • Dạng PD1 (Phân bố không đều đối xứng - mật độ lỗ rỗng lớn nhất ở trục trung hòa và đặc dần về hai phía biên) mang lại độ cứng uốn tổng thể cao nhất.
    • Dạng PD2 (Phân bố bất đối xứng) gây ra hiện tượng lệch trục trung hòa hình học, làm phát sinh mô men uốn phụ và làm gia tăng độ võng cục bộ thêm $24.7%$.
    • Dạng UPD (Phân bố đều) cho độ cứng trung gian nhưng không tối ưu hóa được tỷ lệ độ bền trên trọng lượng.

  3. Hiện tượng Cứng hóa Hình học Phi tuyến (Nonlinear Hardening Spring Behavior): Trong phân tích dao động tự do phi tuyến, tỷ số tần số $\omega_{NL}/\omega_L$ luôn lớn hơn $1.0$ và tăng phi tuyến tính theo tỷ số biên độ dao động đỉnh $w_{max}/h$. Đối với dầm liên kết hai đầu ngàm (CC), hiệu ứng cứng hóa thể hiện mạnh nhất: khi $w_{max}/h = 1.0$, tần số dao động phi tuyến cao hơn tần số tuyến tính đến $45.8%$ do lực kéo màng (in-plane membrane tension) xuất hiện khi chuyển vị lớn.

  4. Vai trò dập tắt dao động và ổn định hóa của nền đàn hồi Pasternak: Hệ số độ cứng cắt $J_0$ của nền Pasternak đóng vai trò quyết định trong việc hạn chế hiện tượng mất ổn định động lực học. Khi tăng tham số $J_0$ từ $0$ lên $10$, độ võng phi tuyến tĩnh của dầm giảm $63.5%$, đồng thời các đỉnh cộng hưởng trong biểu đồ đáp ứng tần số - biên độ bị triệt tiêu mạnh mẽ, ngăn chặn nguy cơ phá hủy mỏi kết cấu dưới tác dụng của tải kích thích điều hòa.

  5. Đánh giá định lượng sự sai khác giữa các Lý thuyết dầm: Đối với dầm có tỷ số độ mảnh lớn ($L/h \ge 20$), kết quả tính toán giữa lý thuyết dầm cổ điển, Timoshenko (FSDBT) và các lý thuyết biến dạng cắt bậc cao (PSDBT, TSDBT, HSDBT, ESDBT) có sự sai khác không đáng kể ($< 2.5%$). Tuy nhiên, đối với dầm dày và dầm ngắn ($L/h = 5$), lý thuyết FSDBT đánh giá quá cao độ cứng của dầm (sai số độ võng lên tới $11.8%$) do không xét đến quy luật biến thiên phi tuyến của ứng suất trượt qua chiều cao tiết diện. Các lý thuyết bậc cao PSDBT, TSDBT và ESDBT cho kết quả trùng khít nhau và mô tả chính xác điều kiện biên ứng suất tiếp triệt tiêu tại các thớ ngoài cùng.

Implications đa chiều

  • Ý nghĩa Lý thuyết và Học thuật: Luận án mở rộng biên giới của cơ học tính toán vật liệu composite bằng việc thiết lập một khung toán học nhất thể hóa cho kết cấu đa tầng dị hướng biến thiên liên tục có xét đến biến dạng lớn von Kármán.
  • Ý nghĩa Phương pháp luận: Cung cấp bộ công cụ thuật toán Pb-Ritz bán giải tích bằng mã nguồn MATLAB hoàn chỉnh, có thể tái lập và mở rộng trực tiếp cho các bài toán phân tích tấm sandwich, vỏ hình trụ và kết cấu không gian phức tạp.
  • Ứng dụng Thực tiễn và Kỹ thuật Công trình: Cung cấp các biểu đồ thiết kế chuẩn hóa (design charts) cho phép các kỹ sư kết cấu lựa chọn cấu hình tối ưu về tỷ lệ chiều dày (cấu hình 1-4-1 và 1-8-1 cho thấy hiệu quả giảm trọng lượng tối ưu nhất mà vẫn đảm bảo độ võng cho phép), hệ số rỗng $e_0 \le 0.5$ và hàm lượng gia cường nano GPL $W_{GPL} = 1.0%$.

Limitations và Future Research

Mặc dù luận án đã đạt được những kết quả học thuật xuất sắc và giải quyết trọn vẹn các mục tiêu nghiên cứu đề ra, tác giả vẫn thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions & limitations) mang tính khách quan:

+-------------------------------------------------------------------------+
|                  GIỚI HẠN VÀ ĐỊNH HƯỚNG PHÁT TRIỂN                      |
+-------------------------------------------------------------------------+
| GIỚI HẠN HIỆN TẠI:                                                      |
| 1. Mô hình hình học 1D dầm phẳng (chưa xét ứng xử 3D/tấm/vỏ phức tạp).  |
| 2. Giả thiết đẳng nhiệt (chưa xét ghép cặp Nhiệt - Điện - Áp điện).     |
| 3. Vật liệu giả định phân tán hoàn hảo (bỏ qua kết tụ nano GPL).        |
|                                                                         |
| ĐỊNH HƯỚNG MỞ RỘNG (10-YEAR ROADMAP):                                   |
| -> Mở rộng phần tử tấm/vỏ sandwich FGM-FGP 2D & 3D.                    |
| -> Mô hình hóa ghép cặp đa trường Đa vật lý (Thermo-Electro-Mechanical).|
| -> Phân tích phá hủy động học, lan truyền vết nứt & Stochastic FE.     |
+-------------------------------------------------------------------------+
  1. Giới hạn về Mô hình Động học Một chiều (1D Beam Kinematics): Nghiên cứu tập trung vào mô hình kết cấu dạng dầm một chiều, giả thiết chiều rộng $b$ không đổi và biến dạng theo phương ngang $y$ được bỏ qua. Mặc dù lý thuyết dầm tổng quát mô tả rất tốt ứng xử uốn phẳng, nó chưa thể hiện đầy đủ các hiệu ứng tập trung ứng suất không gian 3D tại các góc cạnh liên kết phức tạp.
  2. Giả thiết về Môi trường Đẳng nhiệt (Isothermal Conditions): Toàn bộ các phân tích phi tuyến tĩnh và dao động được thực hiện ở nhiệt độ phòng chuẩn ($T = 300\text{ K}$). Ảnh hưởng của gradient nhiệt độ môi trường khắc nghiệt và tải trọng nhiệt kết hợp (thermo-mechanical coupling) chưa được tích hợp trong mô hình toán học này.
  3. Giả thiết Phân tán Lý tưởng của Hạt Nano GPL: Mô hình Halpin-Tsai giả định các phiến nano GPL được phân tán hoàn toàn đồng đều trong chất nền kim loại ở quy mô vi mô. Trong thực tế chế tạo luyện kim bột hoặc in 3D, hiện tượng vón cục (nanoparticle agglomeration) của GPL có thể xảy ra, làm suy giảm một phần cơ tính lý thuyết của vật liệu.

Định hướng nghiên cứu tương lai (Future Research Agenda)

  • Mở rộng khung phân tích Pb-Ritz sang bài toán phân tích phi tuyến tĩnh, mất ổn định và dao động của tấm và vỏ thoải sandwich FGM-FGP gia cường GPL trong không gian hai chiều (2D plate/shell formulation).
  • Nghiên cứu bài toán ghép cặp Đa vật lý (Multiphysics Coupling): Nhiệt - Cơ - Điện - Từ đối với dầm sandwich có tích hợp các lớp bề mặt áp điện (piezoelectric layers) đóng vai trò cảm biến và cơ cấu chấp hành điều khiển dao động chủ động.
  • Phát triển mô hình Cơ học phá hủy và Động học nứt (Fracture & Damage Mechanics): Khảo sát sự hình thành vết nứt vi mô trong ma trận bọt xốp kim loại dưới tải trọng mỏi chu kỳ cao.
  • Áp dụng các phương pháp mô phỏng ngẫu nhiên (Stochastic Finite Element Method - SFEM) để đánh giá độ tin cậy kết cấu khi tính chất vật liệu của GPL và lỗ rỗng có sự biến thiên ngẫu nhiên do sai số chế tạo.

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

  • Tác động Học thuật (Academic Impact): Luận án cung cấp một hệ thống dữ liệu số chuẩn mực cao (benchmark numerical data) với hàng trăm bảng biểu chi tiết về tần số dao động riêng, tỷ số tần số phi tuyến, độ võng và ứng suất không thứ nguyên. Các kết quả này được dự báo sẽ trở thành tài liệu tham khảo cốt lõi cho các nhóm nghiên cứu cơ học tính toán trên thế giới khi 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 (meshfree methods).
  • Chuyển đổi Công nghiệp Vật liệu Mới (Industrial Innovation): Kết quả nghiên cứu chứng minh tính khả thi tuyệt đối của việc ứng dụng vật liệu siêu nhẹ bọt kim loại gia cường graphene trong thiết kế cánh máy bay không người lái (UAV), tấm vỏ tàu cao tốc, thân xe bọc thép chống va đập và kết cấu giàn khoan ngoài khơi. Việc tiết kiệm 30-40% trọng lượng kết cấu trong khi vẫn đảm bảo độ cứng uốn giúp giảm tiêu hao nhiên liệu và khí thải carbon trong giao thông vận tải.
  • Đóng góp cho Tiêu chuẩn Thiết kế Xây dựng & Cơ khí: Khung lý thuyết và các công thức giải tích trong luận án là cơ sở khoa học vững chắc để cập nhật và hoàn thiện các tiêu chuẩn thiết kế kết cấu composite nhiều lớp tiên tiến tại Việt Nam và khu vực.

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

+-------------------------------------------------------------------------------------------------------+
| MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI TỪ CÔNG TRÌNH NGHIÊN CỨU                                                 |
+-------------------------------------------------------------------------------------------------------+
| 1. NGHIÊN CỨU SINH & HỌC GIẢ TRẺ: Khung lý thuyết mẫu mực, giải thuật Pb-Ritz & code mẫu MATLAB.      |
| 2. CÁC NHÀ KHOA HỌC ĐẦU NGÀNH: Khung động học dầm tổng quát, chuẩn đối sánh kiểm định mô hình mới.    |
| 3. KỸ SƯ THIẾT KẾ R&D HÀNG KHÔNG/TÀU THỦY: Bộ biểu đồ thiết kế tối ưu hóa cấu trúc sandwich FGP/GPL.  |
| 4. DOANH NGHIỆP VẬT LIỆU NANO COMPOSITE: Bằng chứng định lượng về tỷ lệ pha GPL tối ưu (1% wt).      |
+-------------------------------------------------------------------------------------------------------+
  • Nghiên cứu sinh và Giảng viên ngành Cơ kỹ thuật, Cơ học Xây dựng: Tiếp cận một cấu trúc luận án mẫu mực, phương pháp biến phân chặt chẽ và hệ thống thuật toán lặp số hoàn chỉnh được lập trình chi tiết trên MATLAB (các Phụ lục B, C, D, E).
  • Các Nhà nghiên cứu Vật liệu Tiên tiến: Nắm bắt được cơ chế tương tác vi mô giữa hạt độn nano 2D (GPL) và cấu trúc ma trận xốp tế bào kín GRF, làm cơ sở định hướng công nghệ chế tạo vật liệu trong phòng thí nghiệm.
  • Kỹ sư Thiết kế Kết cấu Hàng không, Quốc phòng và Công trình biển: Khai thác trực tiếp các bộ thông số hình học tối ưu hóa ($L/h$, $h_c/h$, $p$, $e_0$, $W_{GPL}$) để thiết kế các cấu kiện dầm/sườn chịu uốn và hấp thụ xung lực va đập đạt hiệu năng cơ lý tối đa.

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 lý thuyết độc đáo nhất của luận án là việc xây dựng thành công một mô hình giải tích tổng quát hợp nhất (Unified Generalized Analytical Formulation) kết hợp đồng thời lý thuyết dầm biến dạng cắt bậc cao (HSDT) và quan hệ biến dạng phi tuyến von Kármán cho kết cấu sandwich 3 lớp dị hướng phức hợp (FGM-FGP/GPL). Luận án đã mở rộng trực tiếp lý thuyết dầm biến dạng cắt của Timoshenko (1921) và Reddy (1984) lên hệ kết cấu có cơ tính biến thiên phi tuyến liên tục theo hai phương thức đồng thời: biến thiên pha vật liệu gốm - kim loại (lớp vỏ) và biến thiên vi cấu trúc xốp bọt kim loại gia cường hạt nano 2D GPL (lớp lõi).

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

So với các nghiên cứu trước đây (như nghiên cứu của Sokolinsky và cs., 2004 dùng thí nghiệm và HSDT cho lõi PVC; Chen và cs., 2016 dùng phương pháp Ritz cho dầm FGP thuần nhất; Vo và cs., 2015 dùng nghiệm Navier giải tích bị bó hẹp trong điều kiện biên khớp SS), phương pháp nghiên cứu của luận án có 2 điểm đột phá:

  • Tính tổng quát hóa điều kiện biên: Áp dụng phương pháp Pb-Ritz với chuỗi đa thức độc lập nhân với hàm biên hình học, cho phép giải quyết trơn tru mọi tổ hợp điều kiện biên cơ học phức tạp (CC, SS, CS, CF) mà phương pháp giải tích Navier không thể thực hiện được.
  • Khung thuật toán số đa tầng: Kết hợp đồng bộ ba thuật toán số cao cấp: Thuật toán Newton-Raphson cho bài toán phi tuyến tĩnh, thuật toán lặp xấp xỉ liên tiếp cho dao động tự do phi tuyến và phương pháp Runge-Kutta bậc 4 (RK4) cho bài toán đáp ứng động lực học miền thời gian.

3. Phát hiện thực nghiệm số nào mang tính bất ngờ nhất trong công trình?

Phát hiện bất ngờ nhất nằm ở hiệu ứng đảo ngược ảnh hưởng của độ rỗng khi được định hướng vi cấu trúc: Mặc dù độ rỗng tăng cao ($e_0 = 0.5 \div 0.8$) làm suy giảm mạnh mô đun đàn hồi thể tích của vật liệu gốc, nhưng nếu độ rỗng này được phân bố theo quy luật đối xứng PD1 (tập trung tại trục trung hòa) kết hợp với hàm lượng nhỏ GPL $1%$ phân tán dạng A (tập trung tại các thớ xa trục trung hòa nhất của lớp lõi), độ cứng uốn kháng phi tuyến của dầm sandwich tăng lên vượt trội so với dầm hoàn toàn đặc không có lỗ rỗng. Cấu trúc xốp đã hoạt động như một bộ khung đẩy các lớp vật liệu gia cường GPL ra xa trục trung hòa, tối ưu hóa triệt để mô men quán tính mặt cắt ngang với hiệu quả trọng lượng nhẹ nhất.

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

Có, hoàn toàn đầy đủ và minh bạch. Luận án cung cấp toàn bộ hệ thống biểu thức giải tích tường minh từ Chương 2 đến Chương 4, đồng thời đính kèm toàn bộ mã nguồn chương trình tính toán viết trên ngôn ngữ MATLAB trong phần Phụ lục:

  • Phụ lục A: Thuật toán chi tiết phương pháp giải lặp Newton-Raphson.
  • Phụ lục B: Chương trình MATLAB xác định tính chất hiệu dụng vật liệu vi mô FGM và FGP-GPL.
  • Phụ lục C: Chương trình MATLAB phân tích uốn phi tuyến tĩnh và nội lực.
  • Phụ lục D: Chương trình MATLAB tính toán tần số dao động tự do tuyến tính và phi tuyến.
  • Phụ lục E: Chương trình MATLAB tích phân thời gian Runge-Kutta bậc 4 cho bài toán đáp ứng động lực học.

5. Lộ trình nghiên cứu 10 năm tới từ nền tảng của luận án được vạch ra như thế nào?

Lộ trình phát triển 10 năm (2023–2033) được định hình qua 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1-3 năm): Mở rộng khung Pb-Ritz sang phân tích tấm chữ nhật, tấm nghiêng và vỏ thoải composite sandwich FGM-FGP gia cường GPL chịu tải trọng nhiệt cơ kết hợp.
  • Giai đoạn 2 (4-6 năm): Xây dựng phần mềm mô phỏng phần tử hữu hạn phi tuyến chuyên dụng (User Material Subroutine - UMAT/VUMAT trong ABAQUS) tích hợp mô hình vi cơ học đa tỷ lệ cho kết cấu sandwich hàng không.
  • Giai đoạn 3 (7-10 năm): Chế tạo thực nghiệm mẫu vật liệu bằng công nghệ in 3D kim loại và kiểm chứng phá hủy va đập vận tốc cao (ballistic impact) trên các kết cấu sandwich thực tế.

Kết luận

Luận án tiến sĩ của tác giả Hương Quý Trường là một công trình nghiên cứu khoa học xuất sắc, chuẩn mực, có tính hệ thống cao và đóng góp những giá trị học thuật cốt lõi cho chuyên ngành Cơ kỹ thuật:

  1. Thiết lập thành công mô hình kết cấu dầm sandwich FGP gia cường GPL đa tầng tiên tiến, giải quyết trọn vẹn mâu thuẫn giữa yêu cầu giảm trọng lượng kết cấu và duy trì độ cứng chịu lực cao, đồng thời triệt tiêu nguy cơ tách lớp phân giới nhờ các lớp bề mặt FGM.
  2. Xây dựng khung lý thuyết dầm tổng quát nhất thể hóa, cho phép bao quát và so sánh định lượng độ chính xác của 6 lý thuyết dầm biến dạng cắt khác nhau (CBT, FSDBT, PSDBT, TSDBT, HSDBT, ESDBT) kết hợp với biến dạng phi tuyến von Kármán.
  3. Phát triển bộ công cụ tính toán bán giải tích Pb-Ritz chính xác cao, kết hợp các thuật toán số tiên tiến (Newton-Raphson, lặp xấp xỉ liên tiếp, Runge-Kutta bậc 4), khắc phục triệt để các hạn chế về điều kiện biên của các phương pháp giải tích cổ điển.
  4. Khám phá và chứng minh định lượng cơ chế tối ưu hóa vi cấu trúc: Quy luật độ rỗng đối xứng PD1 kết hợp phân tán GPL Dạng A với hàm lượng $1.0%$ mang lại hiệu năng cơ học uốn và dao động phi tuyến tối ưu nhất.
  5. Cung cấp hệ thống dữ liệu số kiểm chứng toàn diện và bộ mã nguồn MATLAB mở, đóng góp thiết thực cho cộng đồng nghiên cứu cơ học tính toán và mở ra các hướng ứng dụng đột phá trong công nghiệp hàng không vũ trụ và xây dựng hiện đại.