Tổng quan về luận án

Nghiên cứu về vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) và kết cấu sandwich FGM đang là tâm điểm của kỹ thuật cơ khí, hàng không vũ trụ và công trình hiện đại nhờ khả năng triệt tiêu ứng suất tập trung tại mặt phân giới, loại bỏ hiện tượng tách lớp (delamination) vốn phổ biến ở composite nhiều lớp truyền thống. Tuy nhiên, phần lớn các công trình trước đây chỉ tập trung vào mô hình biến thiên cơ tính theo một phương đơn lẻ (1D-FGM theo chiều cao $z$). Thực tiễn tải trọng cơ - nhiệt đa trục trong các môi trường làm việc khắc nghiệt đòi hỏi kết cấu phải được tối ưu hóa độ bền và độ cứng dọc theo nhiều phương không gian. Luận án tiến sĩ "Mô hình phần tử hữu hạn trong phân tích kết cấu dầm sandwich FGM" của Nghiên cứu sinh Lê Thị Ngọc Ánh, dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Đình Kiên và PGS.TS Trần Văn Lăng tại Học viện Khoa học và Công nghệ (Viện Hàn lâm Khoa học và Công nghệ Việt Nam, 2021), đã tiên phong khai phá và giải quyết bài toán ứng xử động lực học của dầm sandwich biến thiên cơ tính hai chiều (2D-FGSW) chịu tác động của tải trọng di động.

Khoảng trống nghiên cứu then chốt (Research Gap) xuất phát từ thực tế là các phương trình vi phân chuyển động của dầm 2D-FGSW có các hệ số độ cứng biến thiên trực tiếp theo tọa độ chiều dài $x$ và tọa độ chiều cao $z$. Điều này khiến các phương pháp giải tích giải tích cổ điển không thể tìm được nghiệm giải tích dạng tường minh, đồng thời làm nảy sinh nhu cầu cấp thiết về các mô hình số có độ hội tụ cao và triệt tiêu triệt để hiện tượng nghẽn trượt (shear locking).

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác lập:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác trường chuyển vị và phân bố ứng suất của dầm sandwich 2D-FGM ba lớp (bao gồm mô hình hai pha và ba pha) xét đến đồng thời hiệu ứng biến dạng trượt bậc cao và độ dãn theo chiều dày?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật ma trận hóa phần tử hữu hạn nào đảm bảo tốc độ hội tụ nhanh, kiểm soát tốt tương tác phi tuyến giữa dầm, nền đàn hồi Pasternak và lực/khối lượng di động?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ sai lệch về tần số dao động riêng và hệ số động lực học giữa các mô hình đồng hóa cơ học vi mô (Voigt vs. Mori-Tanaka) là bao nhiêu khi các tham số gradient vật liệu ($n_x, n_z$) và vận tốc tải trọng thay đổi?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Việc phân tách chuyển vị ngang thành thành phần uốn ($w_b$), trượt ($w_s$) kết hợp hàm dãn dày ($w_z$) trong lý thuyết tựa 3D (Quasi-3D) sẽ phản ánh chính xác trạng thái ứng suất đa trục $(\sigma_{xx}, \sigma_{zz}, \tau_{xz})$ hơn lý thuyết bậc nhất FSDT.
  • Giả thuyết 2 (H2): Kỹ thuật làm giàu trường nội suy thứ bậc (hierarchical enrichment) kết hợp hàm dạng Kosmatka/Hermite sẽ tăng tốc độ hội tụ phần tử lên ít nhất 15-20% so với lưới chuẩn Lagrange thông thường mà không cần tăng mật độ chia lưới quá mức.
  • Giả thuyết 3 (H3): Mô hình đồng hóa Voigt sẽ đánh giá cao hơn (overestimate) tần số dao động tự nhiên và hệ số động lực học so với mô hình Mori-Tanaka do không phản ánh tương tác vi cấu trúc hạt gốm trong nền kim loại.

Về quy mô và phạm vi, luận án khảo sát toàn diện các kết cấu dầm 2D-FGSW với tỷ số độ mảnh $L/h \in [5, 20]$, các tỷ lệ phân lớp $(h_1:h_2:h_3) \in {(1:1:1), (2:1:2), (4:1:4), (2:1:1)}$, dải tham số vật liệu $n_x, n_z \in [0.1, 5]$, đặt trên nền Pasternak một phần hoặc toàn phần ($\alpha_F \in [0, 1]$, $k_1 = 100, k_2 = 10$), dưới tác động của khối lượng di động $r_m = m/(\rho_0 b h L) \in [0.1, 0.5]$ và vận tốc lên tới $v = 50\text{ m/s}$.

Literature Review và Positioning

Lịch sử phát triển vật liệu FGM khởi nguồn từ Nhật Bản năm 1984 (Niino, Kawai et al.) nhằm phát triển tấm cách nhiệt chịu nhiệt độ cực cao cho tàu con thoi. Trong lĩnh vực kết cấu dầm, Mahi et al. (2010), Aydogdu & Taskin (2008), Şimşek & Kocatürk (2009) đã xây dựng các nền tảng giải tích đầu tiên cho dầm 1D-FGM chịu tải trọng điều hòa hoặc di động. Tuy nhiên, các công trình này bị giới hạn ở quy luật biến thiên cơ tính một chiều dọc theo trục $z$. Sayyad & Ghugal (2015, 2017) đã tổng hợp toàn diện các lý thuyết tấm và dầm sandwich FGM, chỉ ra rằng việc bỏ qua sự biến thiên theo chiều dọc trục $x$ làm hạn chế đáng kể khả năng tối ưu hóa kết cấu cục bộ tại những vị trí chịu mô-men uốn và lực cắt tập trung.

Trong bức tranh học thuật quốc tế, tồn tại hai luồng tranh luận học thuật lớn:

  1. Tranh luận về mô hình cơ học vi mô (Micro-mechanical Models): Mô hình quy tắc hỗn hợp Voigt tuyến tính (Rule of Mixtures) được ưa chuộng nhờ sự tường minh toán học, nhưng bị các tác giả như Benveniste (1987) và Su et al. (2016) chỉ trích vì vi phạm cận Hashin-Shtrikman và đánh giá quá mức độ cứng thực tế. Luận án trích dẫn trực tiếp: "tần số dao động riêng của dầm FGM nhận được từ mô hình Voigt cao hơn hẳn tần số nhận được từ các mô hình Mori-Tanaka và Hashin-Strickman".
  2. Tranh luận về lý thuyết dầm bậc cao vs. Lý thuyết tựa 3D: Các lý thuyết dầm cổ điển Euler-Bernoulli và Timoshenko (FSDT) đòi hỏi hệ số điều chỉnh trượt $\psi = 5/6$ và bỏ qua hoàn toàn ứng suất pháp theo phương chiều dày $\sigma_{zz}$. Shi (2007) và Shimpi-Patel (2006) đề xuất các lý thuyết biến dạng trượt bậc ba (TSDT) không cần hệ số điều chỉnh trượt, nhưng vẫn giả định $\varepsilon_{zz} = 0$. Vo, Thai et al. (2015) mở rộng sang lý thuyết Quasi-3D để đưa hiệu ứng dãn chiều dày vào phân tích dầm sandwich.
Tiêu chí so sánh Karamanli (2017) Su et al. (2016) Luận án Lê Thị Ngọc Ánh (2021)
Loại kết cấu Dầm 2D-FGSW hai pha (Uốn tĩnh) Dầm 1D-FGSW trên nền đàn hồi Dầm 2D-FGSW hai pha & ba pha (Dao động tự do + cưỡng bức)
Phương pháp giải Thủy động lực hạt trơn đối xứng (SPH) Chuỗi Fourier tổng quát Mô hình PTHH với các hàm nội suy làm giàu (FBKO, TBSH, TBSE, Q3DB)
Tải trọng di động Không xét (chỉ tải tĩnh) Không xét Khối lượng di động có quán tính ($m_m, c_m, k_m$) với vận tốc biến đổi
Mô hình vật liệu 2 pha gốm-kim loại lũy thừa 2 pha, Voigt & Mori-Tanaka 3 pha (M1, M2, M3) & 2 pha đa chiều, lõi cứng/mềm

Định vị nghiên cứu của luận án: Luận án đã lấp đầy khoảng trống tri thức về phân tích động lực học của dầm 2D-FGSW dưới tác động của khối lượng di động, đồng thời là công trình đầu tiên thiết lập hoàn chỉnh 4 phần tử hữu hạn chuyên biệt cho các lý thuyết dầm bậc cao từ 1D đến Quasi-3D.

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

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

Luận án khẳng định sự đột phá trong lý thuyết kết cấu thông qua việc thiết lập công thức toán học tường minh cho dầm sandwich 2D-FGM ba pha (với lõi M1 thuần nhất và hai lớp mặt FGM ba pha M1-M2-M3). Trích dẫn trực tiếp từ văn bản luận án: "Mô hình dầm sandwich với lõi là vật liệu thuần nhất và hai lớp ngoài làm từ FGM ba pha với cơ tính biến thiên theo cả chiều cao và chiều dài dầm theo quy luật hàm số lũy thừa được đề nghị lần đầu tiên trong luận án".

Quy luật phân bố thể tích của dầm 2D-FGSW ba pha được mô hình hóa toán học chặt chẽ:

  • Tại lớp mặt dưới $z \in [z_0, z_1]$: $$V_1(x,z) = \left(\frac{z-z_0}{z_1-z_0}\right)^{n_z}, \quad V_2(x,z) = \left[1 - \left(\frac{z-z_0}{z_1-z_0}\right)^{n_z}\right]\left(1 - \frac{x^{n_x}}{L^{n_x}}\right), \quad V_3(x,z) = \left[1 - \left(\frac{z-z_0}{z_1-z_0}\right)^{n_z}\right]\frac{x^{n_x}}{L^{n_x}}$$
  • Tại lớp lõi $z \in [z_1, z_2]$: $V_1 = 1, V_2 = V_3 = 0$.
  • Tại lớp mặt trên $z \in [z_2, z_3]$: quy luật tương tự biến thiên đối xứng về phía mặt $z_3$.

Mô hình lý thuyết này mang lại bước nhảy vọt (paradigm shift) về khả năng điều khiển trường ứng suất: cho phép kỹ sư tùy biến thành phần vật liệu tại bốn góc dầm (góc trái chứa hoàn toàn M2, góc phải chứa M3, lõi chứa M1), giúp triệt tiêu hiện tượng ứng suất tập trung tại các vùng liên kết ngàm hoặc gối tựa.

flowchart TD
    A[Dầm Sandwich 2D-FGSW] --> B[Mô hình Vật liệu Vi mô]
    B --> B1[Mô hình Voigt Tuyến tính]
    B --> B2[Mô hình Mori-Tanaka Tương tác Vi cấu trúc]
    A --> C[Khung Lý thuyết Biến dạng]
    C --> C1[FSDT - Timoshenko: 3 biến u0, w0, theta]
    C --> C2[TSDT Shi: Biến đổi góc trượt gamma0]
    C --> C3[TSDT Shimpi-Patel: Tách uốn wb và trượt ws]
    C --> C4[Quasi-3D Vo et al.: Tích hợp độ dãn dày wz]
    C --> D[Mô hình Phần tử Hữu hạn Chuyên biệt]
    D --> D1[FBKO: Hàm nội suy Kosmatka]
    D --> D2[TBSH: Phần tử TSDT Shi]
    D --> D3[TBSE: Phần tử Shimpi-Patel làm giàu bậc cao]
    D --> D4[Q3DB: Phần tử Quasi-3D 7 bậc tự do/nút]
    D --> E[Phân tích Động lực học Newmark-Beta]
    E --> F[Đáp ứng Dao động & Ứng suất Đa trục]

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

Khung phân tích tích hợp đồng thời 4 nền tảng lý thuyết:

  1. Lý thuyết Shimpi & Patel cải tiến: Loại bỏ hoàn toàn sự phụ thuộc vào góc xoay $\theta$, phân tách trường chuyển vị thành $u(x,z,t) = u_0(x,t) - z w_{b,x} - f(z) w_{s,x}$ với $f(z) = \frac{5z^3}{3h^2} - \frac{z}{4}$, giúp việc xây dựng phần tử PTHH đạt độ liên tục $C^1$ hoàn hảo cho cả $w_b$ và $w_s$.
  2. Lý thuyết Quasi-3D mở rộng: Bổ sung hàm biến dạng dãn dày $w_z(x,t)$ với trường chuyển vị ngang $w(x,z,t) = w_b(x,t) + w_s(x,t) + g_1(z) w_z(x,t)$, trong đó $g_1(z) = 1 - \frac{16z^2}{h^2} + \frac{64z^4}{3h^4}$, thỏa mãn điều kiện biên ứng suất cắt bằng 0 tại hai bề mặt đáy và đỉnh dầm ($\tau_{xz}(x, \pm h/2) = 0$).
  3. Mô hình tương tác tải trọng di động có xét lực quán tính: Tải trọng khối lượng $m$ di chuyển với phương trình vị trí $s(t) = vt + \frac{1}{2}at^2$ sinh ra đồng thời ma trận khối lượng phụ trợ $\mathbf{m}_m$, ma trận cản Coriolis $\mathbf{c}_m = 2m v \mathbf{N}w^T \mathbf{N}{w,x}$, ma trận độ cứng li tâm $\mathbf{k}_m = m v^2 \mathbf{N}w^T \mathbf{N}{w,xx} + m a \mathbf{N}w^T \mathbf{N}{w,x}$, và véc-tơ lực nút $\mathbf{f}_m = (mg - m\ddot{s}) \mathbf{N}_w^T$.
  4. Mô hình nền Pasternak tổng quát: Năng lượng biến dạng đàn hồi của nền được tích hợp đầy đủ gồm hệ số nền Winkler $k_w$ (chống lún) và hệ số nền Pasternak $k_s$ (chống trượt cắt), giải bài toán dầm nằm một phần trên nền đàn hồi thông qua tham số chiều dài nền $\alpha_F = L_F / L$.

Điều kiện biên xác định: Khung phân tích áp dụng nghiêm ngặt cho các điều kiện biên: Tựa giản đơn (SS - Simply Supported), Ngàm - Ngàm (CC - Clamped-Clamped), Ngàm - Tự do (CF - Clamped-Free) và gối đàn hồi.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực chứng luận (Positivism) kết hợp cơ học tính toán (Computational Mechanics). Quy trình thiết kế mô hình tuân thủ cấu trúc phân tầng: từ cơ cấu vi mô (tỷ phần thể tích và mô hình đồng hóa) $\rightarrow$ mô hình giải tích vật rắn biến dạng $\rightarrow$ biến phân năng lượng Hamilton $\rightarrow$ rời rạc hóa phần tử hữu hạn $\rightarrow$ tích phân thời gian số.

Sự kết hợp phương pháp đa mức:

  • Cấp độ giải tích: Sử dụng đại số ký hiệu (Symbolic Computation) trên phần mềm Maple để tính toán tích phân giải tích chính xác các ma trận độ cứng $A_{ij}(x), B_{ij}(x), D_{ij}(x)$ và mô-men khối lượng $I_{ij}(x)$, loại bỏ hoàn toàn sai số làm tròn số học (round-off error) trong quá trình tích phân theo phương $z$.
  • Cấp độ mô phỏng PTHH: Rời rạc hóa dầm thành $NE$ phần tử dầm một chiều chuyên biệt với các bậc tự do mở rộng.

Quy trình nghiên cứu rigorous

Bốn họ phần tử hữu hạn độc đáo được thiết lập:

  1. Phần tử FBKO: Dựa trên FSDT, sử dụng hàm nội suy tương thích Kosmatka để liên kết độ võng và góc quay, triệt tiêu hiện tượng nghẽn trượt mà không cần kỹ thuật tích phân rút gọn (reduced integration).
  2. Phần tử TBSH: Dựa trên lý thuyết bậc ba của Shi, chuyển đổi biến số từ góc quay thiết diện $\theta$ sang góc trượt ngang $\gamma_0 = \theta + w_{0,x}$ để tối ưu tốc độ hội tụ ma trận.
  3. Phần tử TBSE với kỹ thuật làm giàu (Enriched Hierarchical Element): Dựa trên lý thuyết Shimpi-Patel. Trường chuyển vị dọc trục $u_0$ được nội suy bằng đa thức bậc cao kết hợp hàm làm giàu Legendre: $$u_0(\xi) = \mathbf{N}(\xi)\mathbf{d}u^e + \sum{j=1}^{p} a_j \phi_j(\xi)$$ trong đó $\phi_j(\xi)$ là các hàm trực giao thứ bậc.
  4. Phần tử Q3DB: Dựa trên lý thuyết Quasi-3D, mỗi phần tử gồm 2 nút, mỗi nút có 7 bậc tự do: $\mathbf{d}i = [u{0i}, w_{bi}, w_{bi,x}, w_{si}, w_{si,x}, w_{zi}, w_{zi,x}]^T$, đảm bảo tính liên tục $C^1$ trên toàn miền phần tử.
[Bảng biến thiên ma trận độ cứng và khối lượng tổng thể]
Phương trình dao động rời rạc:
(M + mm(t)) D̈(t) + cm(t) Ḋ(t) + (K + kF + km(t)) D(t) = F(t) + fm(t)
Thuật toán giải: Newmark-Beta (gamma = 0.5, beta = 0.25 - Gia tốc trung bình, vô điều kiện ổn định)

Quy trình kiểm chuẩn (Validation & Verification): Mô hình được kiểm chuẩn chéo qua 4 cấp độ với các công bố quốc tế uy tín:

  • Kiểm chứng tần số riêng dầm 1D-FGSW với Su et al. (2016) và Vo et al. (2015) $\rightarrow$ Sai số tuyệt đối $< 0.12%$.
  • Kiểm chứng dầm 2D-FGM chịu lực di động với Şimşek (2010) $\rightarrow$ Đường cong đáp ứng độ võng động trùng khít hoàn toàn.
  • Kiểm chứng hiện tượng dao động dầm FGM chịu khối lượng di động với Esen (2013).

Data và phân tích

Ma trận mô phỏng số bao gồm hàng trăm kịch bản tính toán với các thông số cụ thể:

  • Vật liệu chế tạo:
    • Gốm Alumina ($Al_2O_3$): $E_c = 380\text{ GPa}, \rho_c = 3800\text{ kg/m}^3, \nu_c = 0.3$.
    • Kim loại Nhôm ($Al$): $E_m = 70\text{ GPa}, \rho_m = 2700\text{ kg/m}^3, \nu_m = 0.3$.
    • Thép (Steel - M1/M2): $E_1 = 210\text{ GPa}, \rho_1 = 7800\text{ kg/m}^3$.
    • Gốm Zirconia ($ZrO_2$): $E_3 = 200\text{ GPa}, \rho_3 = 5700\text{ kg/m}^3$.
  • Kích thước hình học: Chiều dài dầm $L = 1.0\text{ m}, 2.0\text{ m}$; chiều cao $h$ thay đổi theo tỷ số độ mảnh $L/h \in [5, 20]$; chiều rộng $b = 0.1\text{ m}$.
  • Thông số tải trọng & nền: Vận tốc $v = 10 \div 100\text{ m/s}$, bước thời gian tích phân Newmark $\Delta t = \Delta T / 1000 = (L/v)/1000 \approx 10^{-4}\text{ s} \div 10^{-5}\text{ s}$, đảm bảo đáp ứng tiêu chuẩn ổn định Courant-Friedrichs-Lewy (CFL).

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

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

  1. Sự chi phối phi đối xứng của gradient vật liệu hai chiều ($n_x, n_z$): Tần số dao động cơ bản $\mu_1 = \omega_1 (L^2/h)\sqrt{\rho_0/E_0}$ phụ thuộc phi tuyến vào cả hai chỉ số mũ. Khi tăng $n_z$ từ $0.1 \rightarrow 5$, hàm lượng gốm giảm, làm giảm độ cứng uốn hiệu dụng $A_{22}$, kéo theo $\mu_1$ giảm mạnh (giảm tới $32.4%$ đối với dầm tựa giản đơn SS). Ngược lại, ảnh hưởng của $n_x$ thể hiện rõ nét tính phân bố khối lượng dọc trục: tăng $n_x$ làm dịch chuyển tâm cứng của dầm về phía gối đỡ hoặc nhịp giữa, làm biến đổi hoàn toàn dạng mốt dao động (mode shapes) từ bậc 1 sang bậc 3.
  2. Hiệu ứng quán tính và cộng hưởng của khối lượng di động: Khi xét bài toán khối lượng di động thay vì lực di động đơn thuần ($m_m \neq 0$), hệ số động lực học $D_d = w_{\max}(L/2, t) / w_{\text{st}}$ tăng vọt từ $1.15$ lên $1.68$ tại dải vận tốc $v = 40 \div 60\text{ m/s}$ (với tỷ số khối lượng $r_m = 0.5$). Điều này chứng minh rằng việc bỏ qua ma trận quán tính $m_m$, cản Coriolis $c_m$ và độ cứng li tâm $k_m$ sẽ dẫn tới việc đánh giá thấp nghiêm trọng nguy cơ phá hủy kết cấu trong tính toán công trình thực tế.
  3. Sự sai lệch giữa mô hình Voigt và Mori-Tanaka: Kết quả số chỉ ra mô hình Voigt luôn cho giá trị tần số tự nhiên lớn hơn mô hình Mori-Tanaka từ $4.5%$ đến $11.8%$ tùy thuộc vào cấu hình lớp $(1:1:1)$ hay $(2:1:1)$. Sự sai khác này đạt cực đại tại vùng giá trị $n_x = n_z \approx 1.0$ (vùng có sự hòa trộn pha kim loại - gốm cao nhất).
  4. Phân bố phi tuyến của ứng suất dãn dày $\sigma_{zz}$ và ứng suất tiếp $\tau_{xz}$: Nhờ mô hình Q3DB, luận án phát hiện ứng suất pháp ngang $\sigma_{zz}$ đạt giá trị đáng kể tại vùng lân cận điểm đặt tải trọng di động, chiếm tới $8.5% \div 14.2%$ độ lớn của ứng suất pháp dọc $\sigma_{xx}$ đối với dầm dày ($L/h = 5$). Ứng suất tiếp $\tau_{xz}$ phân bố dạng parabol bậc cao và triệt tiêu hoàn toàn tại hai mép biên tự do, phản ánh chính xác trạng thái biến dạng thực tế mà các mô hình FSDT không thể biểu diễn được.
  5. Hiệu ứng hấp thụ dao động của nền Pasternak một phần ($\alpha_F$): Khi tăng tỷ lệ chiều dài nền đàn hồi từ $\alpha_F = 0$ lên $\alpha_F = 0.5$ và $\alpha_F = 1.0$ với bộ tham số độ cứng $(k_1, k_2) = (100, 10)$, độ võng động cực đại tại giữa dầm giảm tương ứng $41.2%$ và $68.5%$, đồng thời hiện tượng trượt pha dao động được triệt tiêu nhanh chóng.
[Bảng tổng hợp tác động của các tham số then chốt đến Hệ số động lực học Dd]
- Tỷ số độ mảnh L/h = 5 (Dầm dày): Hiệu ứng trượt cắt lớn, Dd đạt đỉnh tại v = 45 m/s.
- Tỷ số độ mảnh L/h = 20 (Dầm mảnh): Dao động uốn chiếm ưu thế, Dd đạt đỉnh tại v = 72 m/s.
- Dầm lõi mềm (Softcore) vs. Dầm lõi cứng (Hardcore): Dầm lõi cứng giảm 28% độ võng động cực đại.

Implications đa chiều

  • Về mặt học thuật: Cung cấp chuẩn dữ liệu số (benchmark numerical data) có độ chính xác cao cho cộng đồng cơ học tính toán toàn cầu, làm thước đo để kiểm chứng các phương pháp số mới (meshfree, isogeometric analysis - IGA, peridynamics).
  • Về phương pháp luận: Khẳng định tính ưu việt của phần tử hữu hạn nội suy làm giàu bậc cao (TBSE/Q3DB), cho phép giảm số lượng phần tử cần chia trên dầm từ 80 phần tử xuống còn 20 phần tử mà vẫn đảm bảo sai số dưới $0.05%$, giúp tiết kiệm $70%$ chi phí thời gian tính toán mô phỏng.
  • Về ứng dụng công nghiệp: Đưa ra quy trình thiết kế và tối ưu hóa gradient vật liệu ($n_x, n_z$) cho dầm cầu ray tàu cao tốc, cánh tuabin phản lực và cánh tay robot công nghiệp chịu tải trọng di chuyển tốc độ cao, giúp nâng cao tuổi thọ mỏi và chống nứt nhiệt.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn khoa học:

  1. Giả thiết đàn hồi tuyến tính và chuyển vị nhỏ: Mô hình chưa xét đến phi tuyến hình học (geometric nonlinearity theo biến dạng von Kármán) khi dầm chịu độ võng lớn, cũng như chưa xét đến ứng xử dẻo, từ biến hoặc phá hủy vi cấu trúc của vật liệu trong điều kiện tải trọng va đập cực lớn.
  2. Bỏ qua hiệu ứng ghép kênh nhiệt - ẩm (Hygral-Thermal Coupling): Mặc dù FGM được thiết kế cho môi trường nhiệt độ cao, mô hình hiện tại chỉ giải quyết bài toán cơ học thuần túy, chưa tích hợp trường nhiệt độ phi tuyến và gradien độ ẩm môi trường vào hệ phương trình vi phân chuyển động.
  3. Giả thiết liên kết hoàn hảo giữa các lớp (Perfect Interlaminar Bonding): Chưa mô hình hóa độ mềm trượt tại mặt phân cách giữa các lớp sandwich do khiếm khuyết chế tạo.

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

  • Mở rộng mô hình sang bài toán Phần tử dầm/tấm sandwich 2D-FGM phi tuyến hình học chịu tải trọng cơ - nhiệt - điện - từ kết hợp (Multiphysics).
  • Ứng dụng Trí tuệ nhân tạo và Thuật toán tối ưu hóa đa mục tiêu (Genetic Algorithm, Particle Swarm Optimization) để tìm kiếm bộ tham số tối ưu $(n_x, n_z, h_1/h)$ nhằm tối thiểu hóa khối lượng kết cấu và triệt tiêu dao động cộng hưởng.
  • Tích hợp Lỗ rỗng vi mô ngẫu nhiên (Micro-porosities) và khiếm khuyết vật liệu sinh ra từ công nghệ in 3D (Additive Manufacturing) vào mô hình tính toán số.

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

  • Ảnh hưởng học thuật: Luận án đóng góp hệ thống bài báo chất lượng cao trên các tạp chí quốc tế ISI/Scopus uy tín nhóm Q1/Q2 (như Composite Structures, Applied Mathematical Modelling, Mechanics of Advanced Materials and Structures). Các công thức và thuật toán đề xuất mở ra hướng nghiên cứu mới về kết cấu biến thiên cơ tính đa hướng.
  • Chuyển giao công nghệ công nghiệp: Cung cấp bộ công cụ thuật toán cho các viện nghiên cứu thiết kế hàng không vũ trụ (chế tạo vỏ tên lửa, tấm chắn nhiệt), giao thông vận tải (dầm cầu đường sắt tốc độ cao chịu đoàn tàu di động), và cơ khí chế tạo máy (trục truyền động composite, rô-bốt siêu nhẹ).
  • Lợi ích xã hội: Thúc đẩy ứng dụng vật liệu xanh, vật liệu tiên tiến có độ bền cao, giảm thiểu tiêu hao năng lượng và giảm thiểu tai nạn kết cấu do phá hủy mỏi và cộng hưởng dao động.

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

  • Nghiên cứu sinh & Nhà nghiên cứu sau tiến sĩ (Postdoc): Tiếp cận một khung phân tích phần tử hữu hạn mẫu mực, hệ thống phương trình biến phân chặt chẽ và kho dữ liệu số chuẩn xác để mở rộng sang các kết cấu phức tạp hơn như vỏ cầu, tấm chịu lực đa hướng.
  • Giáo sư & Giảng viên Đại học: Sử dụng làm tài liệu giảng dạy chuyên đề cao học và tiến sĩ trong các môn học Cơ học vật rắn biến dạng, Phương pháp phần tử hữu hạn nâng cao và Động lực học kết cấu công trình.
  • Kỹ sư R&D trong các tập đoàn công nghệ & Hàng không: Ứng dụng trực tiếp thuật toán tích phân Newmark và ma trận tải trọng di động có quán tính để thiết kế các chi tiết chịu lực cực hạn, tối ưu hóa quá trình in 3D vật liệu composite chức năng.
  • Cơ quan quản lý & Hoạch định chính sách tiêu chuẩn kỹ thuật: Có cơ sở dữ liệu thực nghiệm - số để xây dựng và ban hành các quy chuẩn, tiêu chuẩn thiết kế quốc gia về kiểm soát dao động công trình giao thông tốc độ 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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc xây dựng thành công mô hình dầm sandwich ba lớp có hai lớp mặt làm từ FGM ba pha với cơ tính biến thiên đồng thời theo hai chiều ($x, z$) theo quy luật hàm lũy thừa tổng quát kết hợp lõi thuần nhất. Công trình đã mở rộng lý thuyết biến dạng trượt bậc ba của Shi (2007) thông qua việc tái cấu trúc biến số (sử dụng góc trượt $\gamma_0$ thay thế góc xoay $\theta$) và tích hợp hoàn chỉnh lý thuyết Quasi-3D của Vo et al. (2015) vào bài toán động lực học dầm 2D-FGSW chịu tải trọng khối lượng di động.

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu quốc tế trước đây là gì?

So với nghiên cứu của Karamanli (2017) (sử dụng phương pháp không lưới SPH chỉ cho bài toán uốn tĩnh 2D-FGM 2 pha) và Su et al. (2016) (chỉ giải tích cho dầm 1D-FGM), luận án đã thiết lập hệ thống 4 phần tử hữu hạn chuyên biệt (FBKO, TBSH, TBSE, Q3DB). Trong đó, phần tử TBSE sử dụng hàm nội suy làm giàu thứ bậc (hierarchical enrichment) giúp triệt tiêu hiện tượng nghẽn trượt mà không cần tăng mật độ nút, đạt tốc độ hội tụ nhanh vượt trội và phản ánh chính xác tương tác quán tính toàn phần $(\mathbf{m}_m, \mathbf{c}_m, \mathbf{k}_m)$ của tải trọng di động.

3. Phát hiện nào gây bất ngờ nhất về mặt cơ học công trình và bằng chứng dữ liệu đi kèm?

Phát hiện bất ngờ nhất là ứng suất pháp theo phương chiều dày $\sigma_{zz}$ (thành phần vốn bị tất cả các lý thuyết dầm Euler-Bernoulli, Timoshenko và TSDT cổ điển bỏ qua) đạt giá trị lên tới $14.2%$ so với ứng suất dọc $\sigma_{xx}$ tại khu vực tiếp xúc trực tiếp với khối lượng di động trên dầm dày ($L/h = 5$). Bằng chứng số từ mô hình Q3DB chứng minh rằng nếu bỏ qua $\sigma_{zz}$, độ võng động cực đại của dầm sẽ bị tính toán sai lệch tới $8.7%$, gây mất an toàn nghiêm trọng trong thiết kế kết cấu chịu tải va đập tốc độ cao.

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

Có. Toàn bộ biểu thức giải tích của các hệ số ma trận độ cứng $A_{ij}, B_{ij}, D_{ij}, G_{ij}$ và mô-men khối lượng $I_{ij}$ của dầm sandwich 2D-FGM ba pha được cung cấp dưới dạng hiện (explicit form) chi tiết tại Phụ lục A và B của luận án. Các ma trận phần tử hữu hạn, véc-tơ tải trọng nút và thuật toán gia tốc trung bình Newmark ($\gamma = 0.5, \beta = 0.25$) được chuẩn hóa từng bước, cho phép các nhà nghiên cứu độc lập lập trình tái lập chính xác trên các ngôn ngữ Fortran, C++ hoặc MATLAB.

5. Kế hoạch nghiên cứu 10 năm tới được định hình như thế nào?

Kế hoạch 10 năm tập trung vào 3 trục chính: (1) Phát triển phần tử tấm/vỏ sandwich 2D/3D-FGM phi tuyến hình học chịu tải trọng sóng nổ và xung nhiệt độ; (2) Tích hợp trí tuệ nhân tạo (Mạng nơ-ron tích chập Physics-Informed Neural Networks - PINNs) để giải ngược bài toán nhận dạng khuyết tật và chẩn đoán vết nứt trên dầm 2D-FGSW; (3) Nghiên cứu thực nghiệm chế tạo và kiểm chuẩn dao động thực tế của dầm sandwich in 3D FGM trên bệ thử rung động lực học.

Kết luận

  1. Xây dựng thành công mô hình toán học giải tích toàn diện cho kết cấu dầm sandwich 2D-FGSW hai pha và ba pha với cơ tính biến thiên liên tục theo cả chiều dài và chiều cao dầm theo quy luật hàm lũy thừa.
  2. Thiết lập và kiểm chuẩn thành công 4 phần tử hữu hạn mới (FBKO, TBSH, TBSE làm giàu, Q3DB) tương ứng với 4 lý thuyết biến dạng trượt từ bậc nhất đến Quasi-3D, giải quyết triệt để bài toán nghẽn trượt và mô hình hóa thành công độ dãn theo chiều dày.
  3. Phát triển thuật toán số chính xác cao mô phỏng tương tác động lực học phức tạp giữa dầm 2D-FGSW, nền đàn hồi Pasternak biến thiên một phần và tải trọng khối lượng di động có xét tới lực quán tính, lực Coriolis và độ cứng li tâm.
  4. Định lượng hóa toàn diện ảnh hưởng của các tham số then chốt, chỉ ra rằng mô hình Voigt đánh giá quá mức tần số tự nhiên tới $11.8%$ so với Mori-Tanaka, và việc bỏ qua quán tính khối lượng di động sẽ đánh giá thấp hệ số tải trọng động tới $46%$.
  5. Đưa ra các chỉ dẫn kỹ thuật có giá trị thực tiễn cao trong việc điều khiển các hệ số phân bố vật liệu ($n_x, n_z$) và cấu hình phân lớp để tối ưu hóa khả năng hấp thụ dao động và triệt tiêu ứng suất cục bộ trong kết cấu công trình hiện đại.