Tổng quan về luận án

Nghiên cứu về cơ học kết cấu vật liệu composite dạng thanh nhiều lớp (laminated composite beams) đóng vai trò nền tảng trong các ngành kỹ thuật hàng không, vũ trụ, cơ khí chế tạo và xây dựng công trình hiện đại. Với đặc tính ưu việt như "tỷ số độ cứng trên trọng lượng và độ bền trên trọng lượng cao, hệ số giãn nở nhiệt thấp, tuổi thọ mỏi được nâng cao và khả năng chống ăn mòn vượt trội" (trích xuất từ cơ sở lý thuyết tổng quan của luận án), dầm composite nhiều lớp đang dần thay thế các vật liệu kim loại truyền thống. Tuy nhiên, việc mô hình hóa chính xác ứng xử uốn (bending), ổn định/mất ổn định uốn dọc (buckling) và dao động tự do (free vibration) của kết cấu này đòi hỏi các mô hình giải tích và số học có độ tiệm cận cao với trạng thái ứng suất - biến dạng thực tế.

Khoảng trống nghiên cứu then chốt (specific research gap) được xác định qua công trình tổng quan của Ghugal và Shimpi [2]: phần lớn các lý thuyết dầm đơn lớp tương đương (Equivalent Single Layer Theories - ESLT) hiện hành như lý thuyết dầm cổ điển Euler-Bernoulli (Classical Beam Theory - CBT) hay lý thuyết biến dạng trượt bậc nhất Timoshenko (First-Order Beam Theory - FOBT) thường bỏ qua biến dạng trượt ngang hoặc đòi hỏi hệ số hiệu chỉnh trượt ($k_s$) phụ thuộc phức tạp vào điều kiện biên và hình học mặt cắt. Hơn nữa, ảnh hưởng của hiệu ứng Poisson và biến dạng pháp ngang ($\varepsilon_z$) hiếm khi được tích hợp đồng thời trong các phương trình cấu quan 1D. Ở quy mô vi mô (microbeams), các lý thuyết cơ học môi trường liên tục cổ điển hoàn toàn bất lực trong việc giải thích hiệu ứng phụ thuộc kích thước (size-dependent effects). Đối với phương pháp tính, phương pháp biến phân Ritz tuy có tính tổng quát cao đối với các điều kiện biên tùy ý nhưng việc lựa chọn hàm xấp xỉ (approximation functions) trước đây—chủ yếu là đa thức hoặc đa thức trực giao—thường dẫn đến tốc độ hội tụ chậm, chi phí tính toán cao hoặc mất ổn định số học (numerical instability).

Để giải quyết triệt để các hạn chế trên, luận án tập trung vào 4 câu hỏi và giả thuyết nghiên cứu chính:

  1. RQ1 & H1: Việc đề xuất chuỗi hàm lượng giác (trigonometric functions) và hàm lai (hybrid shape functions) mới trong phương pháp Ritz có đảm bảo thỏa mãn chính xác các điều kiện biên động học (kinematic boundary conditions) và đạt tốc độ hội tụ nhanh hơn các chuỗi đa thức truyền thống không?
  2. RQ2 & H2: Mô hình lý thuyết biến dạng trượt bậc cao (HOBT) và lý thuyết chuẩn 3D (Quasi-3D beam theory) khi tích hợp trực tiếp hiệu ứng Poisson vào phương trình đàn hồi có nâng cao độ chính xác dự báo ứng suất pháp $\bar{\sigma}{xx}$ và ứng suất trượt ngang $\bar{\sigma}{xz}$ so với nghiệm Navier cổ điển không?
  3. RQ3 & H3: Khi tích hợp lý thuyết ứng suất cặp hiệu chỉnh (Modified Couple Stress Theory - MCST) cùng tham số độ dài vật liệu (Material Length Scale Parameter - MLSP), ứng xử tĩnh và động lực học của dầm vi mô composite thay đổi ra sao theo tỷ số $L/h$ và góc xếp lớp sợi?
  4. RQ4 & H4: Ứng xử cơ nhiệt (thermo-mechanical) và đặc tính chịu tải của dầm thành mỏng tiết diện hở (I-beam, Channel-beam) chế tạo từ composite phân lớp hoặc vật liệu cơ tính biến thiên (FG sandwich) có thể được mô phỏng giải tích chính xác thông qua hệ phương trình Lagrange chuyển động tổng quát hay không?

Phạm vi nghiên cứu bao gồm việc xây dựng khung giải tích cho dầm tiết diện chữ nhật đặc, dầm tiết diện thành mỏng hở (I và Channel), dầm vi mô tổng quát và dầm FG sandwich dưới các điều kiện biên: Tựa đơn (Simply-supported: S-S), Ngàm - Tự do (Clamped-Free: C-F), Ngàm - Ngàm (Clamped-Clamped: C-C) và Ngàm - Tựa đơn (Clamped-Simply supported: C-S), khảo sát dải tham số tỷ số chiều dài trên chiều cao ($L/h$ từ 3 đến 50) cùng tỷ số dị hướng vật liệu ($E_1/E_2$ từ 10 đến 40).

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy sự phát triển của các lý thuyết kết cấu dầm composite phân nhánh thành ba trường phái chính:

  • Trường phái lý thuyết biến dạng dầm: Khởi đầu từ lý thuyết dầm Euler-Bernoulli (CBT) bỏ qua hoàn toàn trượt ngang, mở rộng sang lý thuyết Timoshenko (FOBT) [79-81]. Để loại bỏ hệ số hiệu chỉnh trượt nhân tạo, hàng loạt hàm phân bố trượt $f(z)$ bậc cao ra đời: hàm đa thức bậc ba của Reddy [84], Murthy [95], Levinson [94]; hàm trượt hyperbolic của Soldatos [98]; hàm lượng giác của Arya et al. [97], Touratier, Karama et al. [99] và Thai et al. [100]. Các mô hình Quasi-3D nâng cao được phát triển bởi Zenkour [36], Vo và Thai [22, 86, 88] nhằm tính thêm biến dạng pháp ngang $\varepsilon_z$.
  • Trường phái cơ học vi mô phi cổ điển: Khắc phục nhược điểm của môi trường liên tục quy ước khi phân tích kết cấu kích thước micro/nano, Eringen [3] đề xuất lý thuyết đàn hồi phi cục bộ (nonlocal elasticity theory), Mindlin và Aifantis đề xuất lý thuyết gradient biến dạng (strain gradient theory) [4-6], và Yang et al. [7] hoàn thiện lý thuyết ứng suất cặp hiệu chỉnh (MCST) chỉ chứa duy nhất một tham số độ dài vật liệu (MLSP).
  • Trường phái cấu trúc thành mỏng: Dựa trên nền tảng lý thuyết dầm thành mỏng đàn hồi của Vlasov [14], Bauld và Tzeng [15] đã mở rộng cho thanh composite thành mỏng tiết diện hở; Song và Librescu [16], Lee và Kim [17, 18] tiếp tục hoàn thiện mô hình động lực học và mất ổn định cho dầm chữ I và chữ U.

Sự mâu thuẫn và tranh luận học thuật sâu sắc tồn tại giữa hai quan điểm: (1) Việc bỏ qua biến dạng pháp ngang ($\varepsilon_z = 0$) và hiệu ứng Poisson (giả thiết trạng thái biến dạng phẳng $y-yz-xy = 0$) trong các nghiên cứu của Khdeir và Reddy [72, 104] có thực sự mang lại sai số chấp nhận được trong phân tích kết cấu dầm dày ($L/h \le 5$) hay không; (2) Phương pháp giải số như Phần tử hữu hạn (FEM) [20-28] và Lưới phi cấu trúc/Isogeometric [33, 34] tuy xử lý tốt biên phức tạp nhưng đòi hỏi tài nguyên tính toán lớn, trong khi nghiệm giải tích Navier [35, 36] bị giới hạn nghiêm ngặt chỉ giải được biên tựa đơn (S-S).

                      ┌─────────────────────────────────────────────────────────┐
                      │              LÝ THUYẾT DẦM COMPOSITE                    │
                      └──────────────────────────┬──────────────────────────────┘
                                                 │
         ┌───────────────────────────────────────┼────────────────────────────────────────┐
         ▼                                       ▼                                        ▼
┌──────────────────┐                   ┌──────────────────┐                     ┌──────────────────┐
│   CBT & FOBT     │                   │   HOBT Redefined │                     │ Quasi-3D & MCST  │
│(Euler-Bernoulli, │                   │(Reddy, Soldatos, │                     │(Zenkour, Vo-Thai,│
│   Timoshenko)    │                   │Karama, Arya, Thai│                     │Yang et al. MLSP) │
└────────┬─────────┘                   └────────┬─────────┘                     └────────┬─────────┘
         │                                      │                                        │
         └──────────────────────────────────────┼────────────────────────────────────────┘
                                                │
                                                ▼
                      ┌─────────────────────────────────────────────────────────┐
                      │  KHOẢNG TRỐNG: Hàm xấp xỉ Ritz đa thức kém ổn định số;  │
                      │  bỏ qua hiệu ứng Poisson, εz và hiệu ứng vi mô MCST.   │
                      └──────────────────────────┬──────────────────────────────┘
                                                 │
                                                 ▼
                      ┌─────────────────────────────────────────────────────────┐
                      │   ĐÓNG GÓP LUẬN ÁN: Hàm lượng giác/hybrid Ritz mới cho  │
                      │  HOBT/Quasi-3D/MCST & Dầm thành mỏng I/Channel FGM.     │
                      └─────────────────────────────────────────────────────────┘

So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  • So với Mantari và Canales [48], những người sử dụng hàm xấp xỉ lai đa thức - lượng giác (hybrid polynomial-trigonometric functions) gặp hiện tượng hội tụ chậm và kích thước ma trận mở rộng do phải dùng nhân tử Lagrange: giải pháp chuỗi lượng giác của luận án đạt nghiệm hội tụ đối với bài toán mất ổn định chỉ với số số hạng $m = 2$, và dao động tự do với $m = 12$.
  • So với nghiệm giải tích chính xác bằng không gian trạng thái (State-space approach) của Khdeir và Reddy [113] vốn chỉ giới hạn cho dầm xếp lớp trực giao (cross-ply), mô hình của luận án giải quyết trọn vẹn dầm góc tùy ý (angle-ply, arbitrary-ply) và dầm micro dưới mọi điều kiện biên 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 các định luật cơ bản của Cơ học vật liệu composite và Cơ học môi trường liên tục nâng cao thông qua ba trục đóng góp chính:

  1. Mở rộng lý thuyết biến dạng trượt bậc cao (HOBT): Tích hợp hàm phân bố ứng suất trượt dạng lượng giác và hàm mũ thỏa mãn điều kiện biên không ứng suất tiếp trên bề mặt tự do trên và dưới ($z = \pm h/2$): $$\tau_{xz}\left(x, \pm \frac{h}{2}\right) = 0$$ Loại bỏ hoàn toàn sự phụ thuộc vào hệ số trượt $k_s$, chứng minh được sự phân bố parabol chính xác của ứng suất tiếp qua chiều dày kết cấu.
  2. Thiết lập mô hình dầm vi mô phân lớp tổng quát (MGLCB) trên cơ sở MCST: Kết hợp tensor độ cong đối xứng $\chi_{ij}$ và tensor mômen ứng suất bậc cao $m_{ij}$ với tham số kích thước vi mô $l$ (MLSP). Luận án phát hiện cơ chế "tăng cứng kích thước" (size-dependent hardening effect), trong đó tần số dao động riêng và tải trọng tới hạn tăng phi tuyến tính rõ rệt khi kích thước mặt cắt dầm tiến về cấp độ micromet ($h \approx l$).
  3. Mô hình hóa ghép kênh dao động (Coupled Vibration Modes): Luận án chứng minh trên phương diện giải tích rằng "dầm đối xứng thể hiện dao động ghép đôi ($w_0, u_{1a}$) trong khi dầm phản đối xứng xuất hiện dao động ghép ba ($u_0, w_0, u_{1a}$)", làm sáng tỏ sự tương tác phức tạp giữa kéo - uốn - xoắn sinh ra do ma trận độ cứng liên hợp $B_{ij}$ và $B^s_{ij}$.

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

Khung phân tích của luận án là sự hợp nhất có hệ thống giữa ba trụ cột lý thuyết: (1) Lý thuyết tấm/dầm đơn lớp tương đương bậc cao; (2) Lý thuyết đàn hồi vi mô phi cổ điển (MCST); và (3) Nguyên lý biến phân năng lượng cực tiểu thông qua phương trình Lagrange:

$$\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\partial \Pi}{\partial \dot{p}_j}\right) - \frac{\partial \Pi}{\partial p_j} = 0$$

Trong đó, tổng thế năng toàn phần $\Pi = U_E - W + K$ được thiết lập từ biểu thức giải tích tường minh của nội năng biến dạng $U_E$, công của ngoại lực $W$ (tải nén dọc trục $N_0$ và tải phân bố đều $q$), cùng động năng $K$.

Đặc biệt, trong các phương trình cấu quan 1D, định luật Hooke dạng phẳng rút gọn được cải tiến để phản ánh chuẩn xác điều kiện ứng suất phẳng ($\sigma_y = \tau_{xy} = \tau_{yz} = 0$) qua các hệ số độ cứng hiệu chỉnh $C^_{11}$ và $C^_{55}$:

$$C^{11} = C{11} - \frac{C_{12}^2}{C_{22}}, \quad C^{55} = C{55} - \frac{C_{45}^2}{C_{44}}$$

Điều này tạo ra một ranh giới áp dụng chuẩn xác (boundary conditions for validity) giúp khắc phục triệt để sai số đánh giá quá cao độ cứng uốn khi sử dụng hệ số $Q_{11}$ thông thường mà không xét hiệu ứng Poisson.

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): Tiếp cận theo chủ nghĩa thực chứng diễn dịch (Deductive Positivism) và Cơ học tính toán duy lý (Rational Computational Mechanics), trong đó các quy luật cơ học được phát triển từ các tiên đề biến phân năng lượng liên tục, sau đó giải tích hóa bằng phương pháp số và thẩm định chéo với dữ liệu chuẩn (benchmark datasets).
  • Mô hình đa cấp độ (Multi-level Design):
    • Cấp độ vật liệu vi mô/vĩ mô: Từ đơn lớp chỉnh hướng (orthotropic lamina) $\rightarrow$ tấm nhiều lớp phân bố góc sợi $[\theta_k]$ $\rightarrow$ kết cấu vi mô chịu ảnh hưởng của trường ten-xơ vi ứng suất.
    • Cấp độ hình học kết cấu: Dầm tiết diện chữ nhật nguyên khối $\rightarrow$ Dầm thành mỏng tiết diện hở chữ I và chữ C (Channel) với sự tham gia của vật liệu FGM biến thiên cơ tính qua chiều dày theo hàm phân bố thể tích $V_c(z) = (z/h + 1/2)^p$.

Quy trình nghiên cứu rigorous

Quy trình xây dựng trường chuyển vị giải tích và thuật toán Ritz được thiết kế với độ chính xác cao:

  1. Thiết lập trường chuyển vị HOBT: $$u(x, z, t) = u_0(x, t) - z w_{0,x}(x, t) + f(z) u_{1a}(x, t)$$ $$w(x, z, t) = w_0(x, t)$$ Với hàm trượt tối ưu của Kaczkowski/Panc/Reissner: $$f(z) = \frac{5}{4}z\left(1 - \frac{4z^2}{3h^2}\right)$$
  2. Thiết lập chuỗi hàm lượng giác trực tiếp thỏa mãn các điều kiện biên động học:
    • Biên Tựa đơn (S-S): $$\phi_j(x) = \cos\left(\frac{j\pi x}{L}\right), \quad \psi_j(x) = \sin\left(\frac{j\pi x}{L}\right), \quad \varphi_j(x) = \cos\left(\frac{j\pi x}{L}\right)$$
    • Biên Ngàm - Tự do (C-F): $$\phi_j(x) = \sin\left(\frac{(2j-1)\pi x}{2L}\right), \quad \psi_j(x) = 1 - \cos\left(\frac{(2j-1)\pi x}{2L}\right), \quad \varphi_j(x) = \sin\left(\frac{(2j-1)\pi x}{2L}\right)$$
    • Biên Ngàm - Ngàm (C-C): $$\phi_j(x) = \sin\left(\frac{2j\pi x}{L}\right), \quad \psi_j(x) = \sin^2\left(\frac{j\pi x}{L}\right), \quad \varphi_j(x) = \sin\left(\frac{2j\pi x}{L}\right)$$
  3. Triển khai trường chuyển vị vào chuỗi Ritz: $$u_0(x, t) = \sum_{j=1}^m \phi_j(x) u_{0j} \mathrm{e}^{\mathrm{i}\omega t}, \quad w_0(x, t) = \sum_{j=1}^m \psi_j(x) w_{0j} \mathrm{e}^{\mathrm{i}\omega t}, \quad u_{1a}(x, t) = \sum_{j=1}^m \varphi_j(x) u_{1j} \mathrm{e}^{\mathrm{i}\omega t}$$
                ┌────────────────────────────────────────────────────────┐
                │     THIẾT LẬP NĂNG LƯỢNG TOÀN PHẦN (Π = UE - W + K)   │
                └──────────────────────────┬─────────────────────────────┘
                                           │
                                           ▼
                ┌────────────────────────────────────────────────────────┐
                │   CHỌN HÀM XẤP XỈ RITZ LƯỢNG GIÁC THỎA MÃN BIÊN (BCs)  │
                │        (S-S: cos/sin/cos; C-F: sin/1-cos/sin; ...)     │
                └──────────────────────────┬─────────────────────────────┘
                                           │
                                           ▼
                ┌────────────────────────────────────────────────────────┐
                │      ÁP DỤNG PHƯƠNG TRÌNH BIẾN PHÂN LAGRANGE           │
                │         d/dt(∂Π/∂(dpj/dt)) - ∂Π/∂pj = 0                │
                └──────────────────────────┬─────────────────────────────┘
                                           │
                                           ▼
                ┌────────────────────────────────────────────────────────┐
                │      HỆ PHƯƠNG TRÌNH ĐẠI SỐ TUYẾN TÍNH / RIÊNG         │
                │           ([K] - ω²[M]){Δ} = 0  hoặc  [K]{Δ} = {F}     │
                └──────────────────────────┬─────────────────────────────┘
                                           │
                                           ▼
                ┌────────────────────────────────────────────────────────┐
                │  KIỂM TRA HỘI TỤ & ĐỐI SÁNH VỚI DỮ LIỆU CHUẨN (BENCH) │
                │ (Hội tụ tại m=2 cho buckling, m=12 cho vibration)     │
                └────────────────────────────────────────────────────────┘

Data và phân tích

Hệ phương trình cân bằng tĩnh và bài toán riêng được giải trực tiếp trong môi trường tính toán số (MATLAB):

$$\begin{bmatrix} \mathbf{K}^{11} & \mathbf{K}^{12} & \mathbf{K}^{13} \ (\mathbf{K}^{12})^{\mathrm{T}} & \mathbf{K}^{22} & \mathbf{K}^{23} \ (\mathbf{K}^{13})^{\mathrm{T}} & (\mathbf{K}^{23})^{\mathrm{T}} & \mathbf{K}^{33} \end{bmatrix} \begin{Bmatrix} \mathbf{u}_0 \ \mathbf{w}0 \ \mathbf{u}{1a} \end{Bmatrix}

  • \omega^2 \begin{bmatrix} \mathbf{M}^{11} & \mathbf{M}^{12} & \mathbf{M}^{13} \ (\mathbf{M}^{12})^{\mathrm{T}} & \mathbf{M}^{22} & \mathbf{M}^{23} \ (\mathbf{M}^{13})^{\mathrm{T}} & (\mathbf{M}^{23})^{\mathrm{T}} & \mathbf{M}^{33} \end{bmatrix} \begin{Bmatrix} \mathbf{u}_0 \ \mathbf{w}0 \ \mathbf{u}{1a} \end{Bmatrix} = \begin{Bmatrix} \mathbf{0} \ \mathbf{F} \ \mathbf{0} \end{Bmatrix}$$

Dữ liệu đặc tính cơ lý của 4 hệ vật liệu chuẩn (Material Systems) được phân tích:

  • Vật liệu I (MAT I): $E_1/E_2 = 10, 20, 30, 40$; $G_{12}/E_2 = G_{13}/E_2 = 0.6$; $G_{23}/E_2 = 0.5$; $\nu_{12} = \nu_{13} = \nu_{23} = 0.25$.
  • Vật liệu II (MAT II): $E_1 = 144.8\text{ GPa}, E_2 = 9.65\text{ GPa}, G_{12} = G_{13} = 4.14\text{ GPa}, G_{23} = 3.45\text{ GPa}, \nu_{12} = 0.3, \rho = 1389.23\text{ kg/m}^3$.
  • Vật liệu III (MAT III): $E_1 = 181\text{ GPa}, E_2 = 10.3\text{ GPa}, G_{12} = G_{13} = 7.17\text{ GPa}, G_{23} = 3.78\text{ GPa}, \nu_{12} = 0.28, \rho = 1578\text{ kg/m}^3$.

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

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

Các kết quả phân tích số học cung cấp những bằng chứng định lượng chính xác:

Cấu hình dầm & Điều kiện biên Tỷ số $L/h$ Nghiệm đề xuất (Present) Mantari & Canales [20, 48] Khdeir & Reddy [72, 113] Vo & Thai (HOBT) [21, 22]
Độ võng $\bar{w}$ ($0^\circ/90^\circ/0^\circ$) S-S 5 2.665 - 2.665 2.666
Độ võng $\bar{w}$ ($0^\circ/90^\circ/0^\circ$) S-S 10 1.250 1.250 1.250 1.251
Độ võng $\bar{w}$ ($0^\circ/90^\circ/0^\circ$) S-S 20 0.897 0.897 0.897 0.897
Độ võng $\bar{w}$ ($0^\circ/90^\circ/0^\circ$) S-S 50 0.798 0.798 0.798 0.798
Tải uốn dọc tới hạn $\bar{N}_{cr}$ ($0^\circ/90^\circ/0^\circ$) S-S 5 8.906 8.796 8.906 8.666
Tải uốn dọc tới hạn $\bar{N}_{cr}$ ($0^\circ/90^\circ/0^\circ$) C-C 5 11.372 11.365 - -
Tải uốn dọc tới hạn $\bar{N}_{cr}$ ($0^\circ/90^\circ$) S-S 5 4.235 3.886 4.235 4.241
Ứng suất pháp $\bar{\sigma}_{xx}$ ($0^\circ/90^\circ/0^\circ$) 5 1.4311 - - 1.4311
Ứng suất trượt $\bar{\sigma}_{xz}$ ($0^\circ/90^\circ/0^\circ$) 5 0.9484 - - 0.9484
  1. Độ chính xác và tính tương thích tuyệt đối của trường ứng suất: Nghiên cứu chỉ ra rằng ở tỷ số kích thước nhỏ ($L/h = 5$), ứng suất trượt không thứ nguyên $\bar{\sigma}_{xz}$ đạt giá trị đỉnh tại trục trung hòa $z = 0$ và triệt tiêu hoàn toàn tại hai đáy dầm ($z = \pm h/2$), thỏa mãn tuyệt đối điều kiện biên tự do tự nhiên mà không cần bất kỳ hệ số hiệu chỉnh trượt nào.
  2. Hiện tượng suy giảm độ cứng uốn dọc do bất đối xứng: Đối với cấu hình lớp bất đối xứng $(0^\circ/90^\circ)$, hiệu ứng liên hợp kéo - uốn làm giảm tải trọng tới hạn từ $8.906$ xuống $4.235$ (ở $L/h = 5$) và tần số dao động cơ bản giảm gần 50% so với dầm đối xứng $(0^\circ/90^\circ/0^\circ)$.
  3. Hiện tượng khóa trượt và hiệu ứng vi mô kích thước (Size-dependent Stiffening): Khi tham số chiều dài vật liệu vi mô $h/l$ giảm (tức kích thước dầm thu nhỏ về thang micromet), độ võng không thứ nguyên giảm mạnh trong khi tần số tự nhiên tăng hơn 120% đối với dầm tựa đơn.
  4. Hiệu ứng góc xoay sợi $\theta$ trên dầm thành mỏng chữ I và Channel: Khi góc sợi thay đổi từ $0^\circ$ đến $90^\circ$, tần số dao động cơ bản của dầm hộp và dầm chữ I suy giảm phi tuyến theo dạng đường cong chữ S ngược, trong đó sự mất mát độ cứng uốn dọc diễn ra mạnh nhất trong khoảng góc $\theta \in [30^\circ, 60^\circ]$.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp chuẩn so sánh giải tích mới (new analytical benchmark solutions) cho cộng đồng cơ học vật rắn biến dạng, làm thước đo tin cậy để kiểm định các phần tử hữu hạn thương mại (Abaqus, Ansys) hoặc các phần tử HOBT/Isogeometric mới phát triển.
  • Về phương pháp luận: Chuỗi hàm lượng giác và hàm mũ lai trong phương pháp Ritz chứng minh tính ổn định tuyệt đối về mặt thuật toán, khắc phục hiện tượng ma trận độ cứng bị suy biến (ill-conditioned stiffness matrix) vốn thường gặp khi bậc chuỗi đa thức tăng cao.
  • Về mặt thực tiễn kỹ thuật: Cung cấp công cụ tính toán tức thời (real-time design tools) cho các kỹ sư thiết kế cánh máy bay composite, cánh tuabin gió, tay máy robot công nghiệp và hệ thống cảm biến/chấp hành vi cơ điện tử (MEMS/NEMS), giảm thiểu thời gian mô phỏng số 3D phức tạp từ hàng giờ xuống vài mili-giây.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật quan trọng, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Giả thiết tuyến tính hình học và vật liệu: Khung giải tích hiện tại giới hạn ở bài toán biến dạng bé đàn hồi tuyến tính (small strain linear elasticity), chưa xét đến ứng xử phi tuyến hình học von Kármán hay hiện tượng từ biến (creep) và phá hủy phân lớp (delamination) của vật liệu dẻo/nhiệt dẻo.
  2. Giới hạn trường nhiệt độ ổn định: Mô hình cơ nhiệt trong luận án khảo sát độ tăng nhiệt đồng đều hoặc phân bố tuyến tính qua chiều dày dầm, chưa tính đến truyền nhiệt phi tuyến phụ thuộc thời gian và môi trường ẩm phức tạp (hygrothermal analysis).
  3. Mô hình tiếp xúc giữa các lớp: Luận án áp dụng giả thiết dính bám hoàn hảo (perfect bonding) giữa các lớp đơn, chưa mô hình hóa sự trượt giữa các lớp (interlaminar slip) bằng lý thuyết lớp dính đàn hồi (cohesive zone modeling).
  4. Kiểm chứng thực nghiệm: Toàn bộ kết quả dựa trên mô hình giải tích toán học và thẩm định chéo với các công bố ISI uy tín, chưa tiến hành thí nghiệm đo biến dạng thực tế bằng tương quan ảnh kỹ thuật số (DIC) hoặc bẫy vi mô (micro-cantilever test).

Các hướng nghiên cứu tiếp nối trong 5-10 năm tới bao gồm:

  • Mở rộng phương pháp Ritz hàm lượng giác sang bài toán uốn/dao động phi tuyến phi điều hòa (nonlinear post-buckling & large-amplitude vibration).
  • Phân tích động lực học kết cấu dầm thông minh tích hợp lớp áp điện (piezoelectric) và vật liệu nhớ hình (SMA).
  • Nghiên cứu độ tin cậy và tối ưu hóa đa mục tiêu cấu trúc sắp xếp sợi (tailored fiber placement) cho dầm thành mỏng dưới tải trọng va đập động ngẫu nhiên.

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

Luận án của tác giả Ngoc-Duong Nguyen đã tạo ra dấu ấn học thuật mạnh mẽ thông qua 5 công trình công bố trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/Scopus (Q1) như Composite Structures, Composites Part B: Engineering, International Journal of Structural Stability and Dynamics, cùng hàng loạt bài báo trên các tạp chí chuyên ngành uy tín trong nước và hội nghị cơ học toàn quốc.

Ước tính các giải pháp giải tích đề xuất có khả năng thu hút hàng trăm trích dẫn học thuật trong các lĩnh vực Cơ học tính toán và Kỹ thuật vật liệu mới. Về tác động công nghiệp, phương pháp giúp tối ưu hóa trọng lượng kết cấu dầm chịu lực cho các mẫu máy bay hiện đại (như Boeing 787 với hơn 50% vật liệu composite), phương tiện không người lái (UAV), dầm cầu composite FRP trong giao thông vận tải và các thanh rung cộng hưởng siêu nhỏ trong thiết bị cảm biến y sinh học.

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

  • Nghiên cứu sinh & Nhà nghiên cứu cơ học: Tiếp cận nguồn dữ liệu chuẩn (benchmark data) chi tiết, minh bạch về các ma trận khối lượng, ma trận độ cứng và các hệ số tích phân giải tích qua chiều dày để kiểm chứng mô hình giải thuật số.
  • Giảng viên và chuyên gia học thuật: Sử dụng các dẫn xuất biến phân năng lượng và phương trình chuyển động Lagrange làm tài liệu tham khảo giảng dạy chuyên sâu cho các học phần Cơ học kết cấu nâng cao, Vật liệu composite và Phương pháp biến phân.
  • Kỹ sư R&D trong công nghiệp Hàng không & Ô tô: Áp dụng trực tiếp thuật toán và mã nguồn tính toán nhanh để xác định tải trọng phá hủy uốn dọc và dải tần số cộng hưởng, từ đó định cỡ tối ưu hóa góc dệt sợi mà không cần chạy các mô hình 3D FEM cồng kềnh.
  • Các nhà hoạch định chính sách xây dựng & hạ tầng: Có cơ sở dữ liệu định lượng về độ bền mỏi và khả năng chịu tải cơ nhiệt để xây dựng các tiêu chuẩn kỹ thuật quốc gia về ứng dụng dầm composite polymer cốt sợi (FRP) thay thế thép trong các công trình ven biển chịu ăn mòn cao.

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 lý thuyết độc đáo nhất là việc thiết lập thành công mô hình giải tích dầm vi mô composite tổng quát phân lớp tùy ý (MGLCB) trên nền tảng kết hợp Lý thuyết biến dạng trượt bậc cao (HOBT) với Lý thuyết ứng suất cặp hiệu chỉnh (MCST), đồng thời tích hợp trực tiếp hiệu ứng Poisson vào phương trình cấu quan đàn hồi phẳng rút gọn $C^*_{ij}$. Điều này giúp phản ánh chính xác cả hiệu ứng tập trung biến dạng trượt cục bộ lẫn hiệu ứng tăng cứng kích thước ở quy mô hiển vi.

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

So với nghiên cứu của Mantari và Canales [48] (dùng hàm lai đa thức - lượng giác kèm nhân tử Lagrange) và Khdeir - Reddy [104] (dùng phương pháp không gian trạng thái hạn chế ở biên tựa đơn), luận án đã phát triển bộ hàm xấp xỉ lượng giác thuần nhất thỏa mãn tự nhiên toàn bộ các điều kiện biên động học (S-S, C-F, C-C). Đột phá này cho phép ma trận độ cứng giữ nguyên kích thước tối thiểu, loại bỏ hiện tượng khóa bậc tự do và đạt tốc độ hội tụ siêu nhanh ($m = 2$ cho ổn định uốn dọc, $m = 12$ cho dao động tự do).

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

Phát hiện bất ngờ nhất là sự biến đổi cơ chế dao động riêng: trong khi dầm đối xứng $(0^\circ/90^\circ/0^\circ)$ chỉ xuất hiện dao động ghép đôi thuần túy giữa độ võng đứng $w_0$ và góc xoay $u_{1a}$, thì dầm bất đối xứng $(0^\circ/90^\circ)$ lập tức kích hoạt dao động ghép ba đồng thời giữa chuyển vị dọc trục $u_0$, độ võng ngang $w_0$ và góc xoay $u_{1a}$. Bằng chứng số học chỉ ra rằng hiện tượng ghép ba này làm suy giảm tải trọng mất ổn định tới hạn lên đến hơn 52% tại tỷ số $L/h = 5$.

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

Có. Luận án cung cấp tường minh 100% các thành phần ma trận độ cứng $\mathbf{K}{ij}$, ma trận khối lượng $\mathbf{M}{ij}$, các hệ số quán tính $I_0, I_1, I_2, J_1, J_2, K_2$, các hệ số tích phân độ cứng $A, B, D, B^s, D^s, H^s, A^s$, định nghĩa các hàm trượt $f(z)$ và chuỗi hàm xấp xỉ $\phi_j, \psi_j, \varphi_j$ cho từng loại điều kiện biên, tạo điều kiện thuận lợi cho việc tái lập hoàn toàn thuật toán trên mọi ngôn ngữ lập trình (MATLAB, Python, Fortran).

5. Lộ trình nghiên cứu 10 năm được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào 3 giai đoạn: (1) Mở rộng khung giải tích sang bài toán động lực học phi tuyến và phá hủy phân lớp (2026-2028); (2) Tích hợp trường đa vật lý: Cơ - Nhiệt - Điện - Từ trên các cấu trúc dầm/vỏ rỗng nano composite FG (2029-2032); (3) Sản xuất vi cấu trúc bằng công nghệ in 3D và tiến hành thử nghiệm cơ tính thực nghiệm quy mô nano/micro để hoàn thiện bộ tiêu chuẩn thiết kế MEMS (2033-2036).

Kết luận

  1. Đề xuất thành công họ hàm xấp xỉ lượng giác và hàm lai mới trong phương pháp Ritz, đạt độ chính xác tương đương nghiệm chính xác nhưng với chi phí tính toán và bậc tự do ma trận thấp hơn đáng kể.
  2. Xây dựng hoàn chỉnh mô hình giải tích cho dầm composite nhiều lớp dựa trên lý thuyết HOBT và Quasi-3D có tích hợp hiệu ứng Poisson, loại bỏ hoàn toàn hệ số hiệu chỉnh trượt nhân tạo.
  3. Tiên phong giải quyết bài toán cơ học dầm vi mô phân lớp tổng quát (MGLCB) bằng lý thuyết ứng suất cặp hiệu chỉnh (MCST), định lượng chính xác hiệu ứng tăng cứng phụ thuộc kích thước vật liệu.
  4. Mở rộng khung phân tích cho kết cấu dầm thành mỏng tiết diện hở (chữ I và chữ U) chế tạo từ vật liệu FGM và composite nhiều lớp dưới các điều kiện biên phức tạp.
  5. Cung cấp bộ dữ liệu chuẩn (benchmark data) chi tiết về ứng suất, độ võng, tải trọng uốn dọc tới hạn và tần số dao động tự do, phục vụ kiểm chuẩn cho các phương pháp số hiện đại.
  6. Mở ra ba nhánh nghiên cứu đột phá về cấu trúc vi mô, phân tích đa trường vật lý và tối ưu hóa hình học kết cấu thành mỏng trong kỷ nguyên vật liệu tiên tiến toàn cầu.