Tổng quan về luận án

Vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM), được khởi xướng lần đầu tiên vào năm 1984 bởi các nhà khoa học tại Đại học Tohoku (Sendai, Nhật Bản), đã đánh dấu bước chuyển dịch mang tính cách mạng trong khoa học vật liệu và cơ học kết cấu hiện đại. Nhờ cấu trúc vi mô với tỷ phần thể tích các pha thành phần biến đổi liên tục và trơn theo không gian, kết cấu FGM khắc phục triệt để các nhược điểm cố hữu của vật liệu composite nhiều lớp truyền thống như hiện tượng tập trung ứng suất tại mặt phân giới và nguy cơ tách lớp (delamination). Trong bối cảnh công nghệ chế tạo bồi đắp (additive manufacturing) và luyện kim bột phát triển vượt bậc, các cấu kiện dạng tấm và dầm sandwich FGM biến thiên cơ tính hai chiều (2D-FGM) với hai hoặc ba pha vật liệu thành phần (kim loại - gốm - chất nền) đã mở ra tiềm năng ứng dụng đột phá cho các môi trường vận hành cơ - nhiệt siêu khắc nghiệt trong ngành hàng không vũ trụ, năng lượng hạt nhân và cơ khí giao thông tốc độ cao.

Tuy nhiên, bức tranh học thuật quốc tế thời điểm trước công trình này bộc lộ một khoảng trống nghiên cứu sâu sắc (research gap): hầu hết các nghiên cứu về kết cấu sandwich FGM chỉ tập trung vào mô hình biến thiên một chiều theo phương chiều dày (1D-FGM dọc trục $z$). Các cấu kiện chịu tải trọng động lực học thực tế như cánh turbine, tấm chắn nhiệt hay dầm cầu đường sắt cao tốc phải đồng thời chịu gradient nhiệt - cơ phân bố phức tạp theo cả chiều cao ($z$) và chiều dài kết cấu ($x$). Hơn thế nữa, các nghiên cứu hiện có về 2D-FGM chủ yếu dừng lại ở cấu trúc hai pha dưới tải trọng tĩnh hoặc dao động tự do cơ bản; các bài toán dao động cưỡng bức của dầm sandwich 2D-FGM hai và ba pha chịu tải trọng di động động học (lực điều hòa di động, khối lượng di động có xét tới lực quán tính, lực Coriolis và lực ly tâm) kết hợp phân tích đa mô hình đồng nhất hóa vi cơ học (Voigt, Mori-Tanaka, Maxwell) gần như chưa được khai phá một cách toàn diện và tường minh.

Để giải quyết khoảng trống học thuật này, luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01) của nghiên cứu sinh Phạm Vũ Nam, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Đình Kiên và TS. Nguyễn Văn Chình 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, đã tập trung giải quyết 3 câu hỏi nghiên cứu cốt lõi:

  1. Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để xây dựng mô hình toán học và phần tử hữu hạn độ chính xác cao, khử triệt để hiện tượng nghẽn trượt (shear locking) cho bài toán phân tích dao động dầm và tấm sandwich 2D-FGM hai và ba pha?
  2. Câu hỏi nghiên cứu 2 ($RQ_2$): Gradient biến thiên vật liệu hai chiều ($n_x, n_z$) và tỷ số hình học lớp sandwich (cấu hình 1-1-1, 2-1-2, 1-2-1) chi phối như thế nào tới phổ tần số dao động tự nhiên và đáp ứng quá độ của kết cấu?
  3. Câu hỏi nghiên cứu 3 ($RQ_3$): Mức độ sai lệch định lượng giữa các mô hình đồng nhất hóa vi cơ học (Voigt, Mori-Tanaka, Maxwell) ảnh hưởng ra sao tới tần số dao động riêng và hệ số đáp ứng động lực học ($D_d$) khi chịu tải trọng di động có vận tốc cao?

Các giả thuyết nghiên cứu tương ứng được xác lập:

  • Giả thuyết $H_1$: Việc kết hợp lý thuyết biến dạng trượt bậc ba cải tiến (TSDT) và lý thuyết trượt lượng giác với các hàm nội suy làm giàu (enriched interpolation) trong phần tử dầm 2 nút 10 bậc tự do sẽ tối ưu hóa tốc độ hội tụ và mô tả chính xác ứng xử động học mà không cần sử dụng hệ số điều chỉnh trượt.
  • Giả thuyết $H_2$: Lực quán tính Coriolis ($2mv\dot{w}{,x}$) và lực ly tâm ($mv^2 w{,xx}$) từ khối lượng di động làm thay đổi đáng kể phân bố ứng suất và biên độ dao động cực đại so với giả định tải trọng lực tập trung đơn thuần.
  • Giả thuyết $H_3$: Mô hình Voigt đánh giá quá cao độ cứng kết cấu so với mô hình Mori-Tanaka và Maxwell, dẫn tới sự chênh lệch có hệ thống về tần số tự nhiên (từ 5% đến 15%) tùy thuộc vào chỉ số quy luật lũy thừa $n_x, n_z$.

Khung lý thuyết của công trình được định hình dựa trên sự tích hợp của Lý thuyết biến dạng trượt bậc ba cải tiến Bickford - Reddy, Lý thuyết biến dạng trượt lượng giác Touratier, Lý thuyết tấm Mindlin mở rộng với kỹ thuật nội suy liên kết (linked interpolation), Nguyên lý biến phân Hamilton và Phương pháp tích phân thời gian trực tiếp Newmark-$\beta$. Phạm vi nghiên cứu bao quát các kết cấu dầm sandwich 2D-FGM với tỷ số nhịp - chiều dày $L/h \in [5, 20]$, tấm sandwich 2D-FGM ba pha với tỷ số $a/h \in [5, 20]$, dải vận tốc tải trọng $v \le 80\text{ m/s}$, tham số khối lượng di động tương đối $r_m = m/(\rho_0 A L) \le 0.5$, và chỉ số biến thiên quy luật lũy thừa $n_x, n_z \in [0, 5]$.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử cơ học của kết cấu dầm và tấm FGM ghi nhận nhiều nhánh phát triển chính. Nhánh thứ nhất tập trung vào dao động tự do của dầm 1D-FGM có cơ tính biến đổi ngang, tiêu biểu với các công trình giải tích của Aydogdu và Taskin (2007), Sina và cộng sự (2009), Thai và Vo (2012) sử dụng lý thuyết biến dạng trượt bậc cao nhằm loại bỏ hệ số hiệu chỉnh trượt $k_s$. Nhánh thứ hai nghiên cứu dầm FGM có cơ tính biến đổi dọc (Axially FGM), nổi bật là các công bố của Shahba và cộng sự (2011), Huang và Li (2010), Alshorbagy và cộng sự (2011), trong đó các phương pháp phần tử hữu hạn và ma trận độ cứng động lực học được khai thác để giải quyết hệ phương trình vi phân có hệ số phụ thuộc tọa độ không gian. Nhánh thứ ba mở rộng sang kết cấu 2D-FGM, với các nghiên cứu giải tích hàm số mũ của Wang và cộng sự (2016), phương pháp Ritz đa thức của Karamanlı (2017), và phương pháp phần tử đẳng hình học (IGA) của Lieu và cộng sự (2018).

Trong dòng chảy học thuậ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 đồng nhất hóa vi cơ học: Quy tắc trộn Voigt (1889) giả định trường biến dạng đồng nhất bỏ qua tương tác vi mô giữa các pha, dẫn đến việc ước lượng quá cao độ cứng của vật liệu. Ngược lại, mô hình tự nhất quán Mori-Tanaka (1973) và mô hình Maxwell (1873) cho vật liệu đa pha xét tới tương tác vi mô dạng hạt vùi hình cầu trong môi trường nền liên tục, phản ánh sát thực hơn trạng thái ứng suất nội tại nhưng làm tăng tính phi tuyến và phức tạp trong tích phân ma trận độ cứng.
  2. Tranh luận về bậc của lý thuyết biến dạng trượt: Lý thuyết trượt bậc nhất Mindlin (FSDT) đòi hỏi hệ số hiệu chỉnh cắt phụ thuộc phức tạp vào hình học và quy luật biến thiên vật liệu; trong khi lý thuyết biến dạng trượt bậc cao (TSDT) và lý thuyết trượt lượng giác thỏa mãn tự nhiên điều kiện triệt tiêu ứng suất trượt trên hai bề mặt tự do $\tau_{xz}(x, \pm h/2) = 0$ mà không cần bất kỳ hệ số hiệu chỉnh nhân tạo nào.

Về mặt vị thế học thuật, công trình của NCS Phạm Vũ Nam đã tiến xa hơn các nghiên cứu quốc tế đương đại thông qua việc đối chuẩn trực tiếp với hai nghiên cứu cơ sở:

  • So sánh với nghiên cứu của Songsuwan và cộng sự (2015) [62] vốn chỉ khảo sát dầm sandwich 1D-FGM chịu lực di động bằng phương pháp Ritz; luận án đã mở rộng thành công lên cấu hình 2D-FGM ba pha chịu khối lượng di động toàn phần có xét tương tác động lực học phức tạp.
  • So sánh với công bố của Karamanlı (2018) [135] vốn chỉ phân tích tĩnh dầm sandwich 2D-FGM bằng phương pháp giải tích đa thức; luận án đã thiết lập một hệ thống phần tử hữu hạn 10-DOF đa năng có khả năng giải quyết chính xác cả bài toán dao động tự do lẫn đáp ứng quá độ phi dừng dưới tải trọng di động điều hòa và khối lượng di động.

Đó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 và làm sâu sắc thêm các lý thuyết cơ học kết cấu phi cổ điển thông qua các đóng góp cụ thể:

  • Mở rộng Lý thuyết biến dạng trượt bậc ba Bickford - Reddy cho dầm sandwich 2D-FGM hai pha bằng cách phân tách trường chuyển vị ngang thành thành phần uốn ($w_b$) và trượt ($w_s$): $$u_1(x, z, t) = u(x, t) - z w_{b,x} - z \left[ 1 - \frac{4}{3}\left(\frac{z}{h}\right)^2 \right] w_{s,x}$$ $$u_3(x, z, t) = w(x, t) = w_b(x, t) + w_s(x, t)$$ Trích dẫn nguyên văn từ luận án: "với việc lựa chọn trường chuyển vị hợp lý, ứng suất trượt trong lý thuyết này thỏa mãn điều kiện tự do ứng suất trên hai mặt của dầm, vì thế không cần đưa vào các hệ số điều chỉnh trượt".
  • Thiết lập hệ phương trình vi phân chuyển động hoàn chỉnh cho dầm sandwich 2D-FGM ba pha trên cơ sở Lý thuyết biến dạng trượt lượng giác Touratier, tích hợp chính xác các số hạng gia tốc đối lưu sinh ra do khối lượng di động $m$: gia tốc quán tính ($m\ddot{w}$), gia tốc Coriolis ($2mv\dot{w}{,x}$) và gia tốc ly tâm ($mv^2 w{,xx}$).
  • Xác lập 3 mệnh đề lý thuyết (Propositions) đóng vai trò nền tảng:
    • Mệnh đề 1: Gradient vật liệu hai chiều phá vỡ tính đối xứng của trường ứng suất dọc trục $\sigma_{xx}$, đẩy trục trung hòa dịch chuyển phi tuyến theo cả hai phương $x$ và $z$.
    • Mệnh đề 2: Lực Coriolis đóng vai trò như một lực cản động lực học không bảo toàn phụ thuộc vận tốc, làm dịch chuyển vị trí xuất hiện độ võng cực đại về phía nửa sau của nhịp dầm ($x > L/2$).
    • Mệnh đề 3: Tần số dao động tự nhiên của dầm/tấm sandwich 2D-FGM thỏa mãn tính đơn điệu giảm đối với các chỉ số mũ phân bố $n_x$ và $n_z$ khi tỷ phần thể tích kim loại tăng dần.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần lý thuyết: Cơ học môi trường liên tục biến phân Hamilton, Lý thuyết vi cơ học đồng nhất hóa vật liệu (Voigt, Mori-Tanaka, Maxwell), và Phương pháp phần tử hữu hạn phi cổ điển (Enriched/Linked Finite Element Methods). Mô hình tỷ phần thể tích hàm lũy thừa tổng quát cho hai lớp ngoài dầm sandwich 2D-FGM hai pha được định nghĩa: $$V_m(x, z) = \left( \frac{z - z_1}{z_0 - z_1} \right)^{n_z} \left( 1 - \frac{x}{2L} \right)^{n_x}, \quad z \in [z_0, z_1]$$ $$V_m(x, z) = 0, \quad z \in [z_1, z_2] \quad (\text{Lớp lõi thuần gốo})$$ $$V_m(x, z) = \left( \frac{z - z_2}{z_3 - z_2} \right)^{n_z} \left( 1 - \frac{x}{2L} \right)^{n_x}, \quad z \in [z_2, z_3]$$ với $V_c(x, z) = 1 - V_m(x, z)$. Điều kiện biên liên kết hoàn hảo (perfect bonding) được áp dụng tại các mặt phân giới $z_1, z_2$, bảo đảm tính liên tục tuyệt đối của trường chuyển vị và ứng suất cắt.

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ủ triết lý nghiên cứu thực chứng diễn dịch (Deductive Positivism) kết hợp chủ nghĩa hiện thực cơ chế (Mechanistic Realism). Thiết kế nghiên cứu tính toán được cấu trúc hóa theo mô hình đa tầng (multi-level computational design):

  • Cấp độ 1 (Micro-scale): Đồng nhất hóa tính chất đàn hồi hiệu dụng ($E_f, G_f, K_f, \nu_f, \rho_f$) dựa trên các mô hình vi cơ học giải tích.
  • Cấp độ 2 (Meso/Macro-scale): Xây dựng quan hệ biến dạng - ứng suất, thiết lập ma trận độ cứng độ mềm mặt cắt tích phân theo chiều dày và nhịp kết cấu.
  • Cấp độ 3 (Structural Dynamics): Thiết lập hệ phương trình đại số động lực học phụ thuộc thời gian và giải toán số trị.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các bước chuẩn hóa:

  1. Rời rạc hóa phần tử dầm hai pha: Phần tử dầm hai nút với 10 bậc tự do được phát triển, vec-tơ chuyển vị nút phần tử có dạng $\mathbf{d}e = { u_1, u_2, w{b1}, w_{b1,x}, w_{b2}, w_{b2,x}, w_{s1}, w_{s1,x}, w_{s2}, w_{s2,x} }^T$. Hàm nội suy Lagrange bậc nhất được dùng cho chuyển vị dọc trục $u$, đa thức Hermite bậc ba được dùng đồng thời cho $w_b$ và $w_s$.
  2. Rời rạc hóa phần tử dầm ba pha làm giàu: Sử dụng trường nội suy lượng giác bổ sung để làm giàu cơ sở hàm dạng Hermite, giúp mô tả chính xác gradient biến dạng cục bộ do khối lượng di động gây ra, loại bỏ hoàn toàn hiện tượng nghẽn trượt.
  3. Rời rạc hóa tấm sandwich 2D-FGM ba pha: Xây dựng phần tử tứ giác 9 nút (Q9) trên cơ sở lý thuyết tấm Mindlin, tích hợp phương pháp nội suy liên kết (linked interpolation) nhằm khử triệt để hiện tượng nghẽn trượt cắt khi tỷ số độ mảnh $a/h \to 20$.
  4. Thiết lập phương trình động lực học tổng thể: $$\mathbf{M}(t) \ddot{\mathbf{D}}(t) + \mathbf{C}(t) \dot{\mathbf{D}}(t) + \mathbf{K}(t) \mathbf{D}(t) = \mathbf{F}(t)$$ trong đó ma trận khối lượng $\mathbf{M}(t) = \mathbf{M}_{beam} + \mathbf{M}_m(t)$, ma trận cản Coriolis $\mathbf{C}(t) = \mathbf{C}m(t)$, và ma trận độ cứng $\mathbf{K}(t) = \mathbf{K}{beam} + \mathbf{K}_m(t)$ biến thiên liên tục theo vị trí tải trọng $x_m = vt$.
  5. Giải thuật số: Ứng dụng phương pháp tích phân trực tiếp Newmark-$\beta$ với tham số không tiêu tán năng lượng số ($\beta = 0.25, \gamma = 0.5$), bước thời gian chia nhỏ $\Delta t = \Delta T / 1000$ (với $\Delta T = L/v$).

Data và phân tích

Dữ liệu vật liệu thành phần được kiểm định chuẩn mực từ y văn quốc tế:

  • Pha gốm Zirconia ($ZrO_2$): $E_c = 151\text{ GPa}$, $\rho_c = 3000\text{ kg/m}^3$, $\nu_c = 0.3$.
  • Pha kim loại Nhôm (Al): $E_m = 70\text{ GPa}$, $\rho_m = 2707\text{ kg/m}^3$, $\nu_m = 0.3$.
  • Pha gốm Alumina ($Al_2O_3$): $E_1 = 380\text{ GPa}$, $\rho_1 = 3800\text{ kg/m}^3$, $\nu_1 = 0.3$.
  • Pha kim loại Thép (SUS304): $E_3 = 201.04\text{ GPa}$, $\rho_3 = 8166\text{ kg/m}^3$, $\nu_3 = 0.3$.
  • Pha Ti-6Al-4V: $E_2 = 105.7\text{ GPa}$, $\rho_2 = 4429\text{ kg/m}^3$, $\nu_2 = 0.2988$.

Các tích phân độ cứng và khối lượng dạng giải tích tường minh không thể biểu diễn trực tiếp do hệ số phụ thuộc tọa độ không gian $x$ được xử lý chính xác tuyệt đối bằng thuật toán tích phân số Cầu phương Gauss (Gauss-Legendre Quadrature) cấp cao. Toàn bộ mã nguồn phần tử hữu hạn và thuật toán Newmark được lập trình độc lập trên môi trường phần mềm MATLAB kết hợp hệ thống tính toán biểu thức hình thức Maple.

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

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

Kết quả mô phỏng số và phân tích thực nghiệm số trong luận án mang lại 5 phát hiện cốt lõi:

  1. Định lượng ảnh hưởng của gradient 2D ($n_x, n_z$): Khi tăng chỉ số mũ vật liệu $n_x$ và $n_z$ từ $0$ lên $5$, tỷ phần pha gốm giảm mạnh, dẫn tới tham số tần số dao động cơ bản không thứ nguyên $\mu_1 = \omega_1 (L^2/h)\sqrt{\rho_0/E_0}$ của dầm sandwich 2D-FGM hai pha ($L/h = 20$, cấu hình 2-1-2) giảm liên tục từ $5.1264$ xuống $3.4218$ (giảm xấp xỉ $33.25%$).
  2. Hiệu ứng vận tốc tới hạn và lực Coriolis: Dưới tác động của tải trọng di động, hệ số động lực học $D_d = \max |w(L/2, t)| / w_{static}$ đạt cực đại tại dải vận tốc cộng hưởng $v \approx 65 - 75\text{ m/s}$. Khối lượng di động ($r_m = 0.5$) làm tăng độ võng động cực đại thêm $24.8%$ so với mô hình tải trọng lực thuần túy do hiệu ứng quán tính kéo theo.
  3. Định lượng sai số giữa các mô hình đồng nhất hóa vi cơ học: Mô hình Voigt luôn dự báo tần số dao động riêng cao hơn mô hình Mori-Tanaka từ $4.2%$ đến $8.6%$ đối với dầm hai pha; và cao hơn mô hình Maxwell từ $6.1%$ đến $12.4%$ đối với dầm/tấm ba pha. Điều này khẳng định việc bỏ qua tương tác vi mô trong quy tắc trộn Voigt sẽ gây mất an toàn nghiêm trọng khi thiết kế dầm sandwich chịu tải động.
  4. Tối ưu hóa cấu hình lớp sandwich: Cấu hình phân lớp (2-1-2) với lõi thuần gốm cho khả năng kháng uốn động học và triệt tiêu ứng suất cắt tiếp tuyến $\tau_{xz}$ tại các mặt biên tốt hơn vượt bậc so với cấu hình (1-1-1) và (1-2-1), giảm ứng suất cắt đỉnh tới $18.5%$.
  5. Tính ổn định số và kháng nghẽn trượt của phần tử Q9: Phần tử tấm tứ giác 9 nút với nội suy liên kết duy trì độ chính xác hội tụ với sai số $< 0.12%$ trên toàn dải tỷ số hình học từ tấm dày ($a/h = 5$) đến tấm mỏng ($a/h = 20$).

Implications đa chiều

  • Về lý thuyết và phương pháp luận: Luận án cung cấp một hệ thống lý thuyết nhất quán kết hợp giữa cơ học môi trường liên tục bậc cao và vi cơ học hiện đại, mở đường cho việc phân tích các kết cấu composite biến thiên đa chiều phức tạp.
  • Về ứng dụng kỹ thuật và thực tiễn: Cung cấp công cụ tính toán định lượng chính xác cho các kỹ sư kết cấu trong việc tối ưu hóa tỷ phần thể tích gốm/kim loại nhằm thiết kế dầm cầu chịu đoàn tàu cao tốc, cánh nâng máy bay siêu âm và các tấm bảo vệ chịu xung nhiệt - cơ trong lò phản ứng hạt nhân.
  • Về chuẩn hóa thiết kế: Khuyến cáo các tiêu chuẩn kỹ thuật không nên sử dụng đơn thuần mô hình quy tắc trộn Voigt trong tính toán kết cấu sandwich FGM chịu tải trọng di động mà bắt buộc phải áp dụng các mô hình đồng nhất hóa vi cơ học bậc cao (Mori-Tanaka hoặc Maxwell) để tránh đánh giá thấp biên độ dao động thực tế.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu mang tính biên (boundary conditions):

  1. Giới hạn đàn hồi tuyến tính: Nghiên cứu chỉ giả định vật liệu làm việc trong miền đàn hồi tuyến tính Hooke với giả định biến dạng nhỏ, chưa xét đến tính phi tuyến hình học Von Kármán khi kết cấu chịu độ võng lớn.
  2. Chưa ngẫu hợp trường nhiệt độ động: Các tính chất cơ lý của vật liệu ($E, \nu, \rho$) được giả định không phụ thuộc vào sự thay đổi tức thời của nhiệt độ môi trường trong quá trình dao động.
  3. Giả định tiếp xúc lý tưởng giữa các lớp: Bỏ qua ứng suất bóc tách cục bộ và hiện tượng trượt vi mô tại mặt tiếp xúc giữa các lớp sandwich.
  4. Phạm vi hình học phẳng: Nghiên cứu giới hạn ở kết cấu dầm thẳng và tấm phẳng chữ nhật, chưa mở rộng cho kết cấu vỏ cong hoặc tấm có biên dạng phức tạp.

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

  • Hướng 1: Mở rộng mô hình phần tử hữu hạn phân tích động lực học phi tuyến hình học và phi tuyến vật liệu cho dầm/tấm 2D-FGM chịu tải trọng cơ - nhiệt ngẫu hợp.
  • Hướng 2: Tích hợp ảnh hưởng của độ rỗng vi mô (porosities) và khuyết tật hình học ngẫu nhiên trong quá trình chế tạo bồi đắp 3D-printing FGM.
  • Hướng 3: Phát triển các phần tử vỏ nón, vỏ trụ sandwich 2D-FGM ba pha phục vụ trực tiếp cho thân vỏ tên lửa và tàu không gian.
  • Hướng 4: Nghiên cứu bài toán điều khiển chủ động dao động kết cấu 2D-FGM tích hợp các lớp vật liệu áp điện thông minh (piezoelectric layers).

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

  • Ảnh hưởng học thuật (Academic Impact): Dự kiến tạo lập nguồn trích dẫn học thuật quan trọng cho các nhà nghiên cứu trong lĩnh vực cơ học tính toán kết cấu FGM, đặc biệt trong phân tích tải trọng di động và đồng nhất hóa vi cơ học.
  • Chuyển đổi công nghiệp (Industry Transformation): Hỗ trợ ngành hàng không vũ trụ và công nghiệp quốc phòng tối ưu hóa thiết kế vật liệu giảm trọng lượng kết cấu nhưng vẫn đảm bảo độ cứng vững động lực học và khả năng kháng mỏi cơ nhiệt.
  • Ảnh hưởng chính sách và quy chuẩn (Policy & Standards): Đóng góp cơ sở dữ liệu số chuẩn xác để xây dựng các hướng dẫn kỹ thuật thiết kế kết cấu FGM trong xây dựng công trình giao thông tốc độ cao tại Việt Nam.

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

  • Nghiên cứu sinh và học giả chuyên ngành Cơ học/Cơ kỹ thuật: Tiếp cận phương pháp lập trình phần tử hữu hạn làm giàu, khử nghẽn trượt và thuật toán Newmark tích hợp tải trọng động học phức tạp.
  • Kỹ sư R&D hàng không, vũ trụ và ô tô: Sử dụng trực tiếp quy trình tối ưu hóa tỷ lệ phân bố vật liệu ($n_x, n_z$) và cấu hình sandwich để nâng cao độ bền kết cấu.
  • Các nhà hoạch định tiêu chuẩn kỹ thuật kết cấu: Nắm bắt được biên độ an toàn và độ lệch định lượng giữa các mô hình vi cơ học khi xây dựng tiêu chuẩn tính toán kết cấu composite tiên tiến.

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?

Trả lời: Đóng góp độc đáo nhất là việc mở rộng thành công Lý thuyết biến dạng trượt lượng giác kết hợp hàm làm giàu (enriched interpolation) cho dầm sandwich 2D-FGM ba pha chịu khối lượng di động toàn phần. Mô hình đã tích hợp đồng thời và tường minh 3 thành phần động học: quán tính khối lượng di động ($m\ddot{w}$), lực cản Coriolis đối lưu ($2mv\dot{w}{,x}$), và lực ly tâm hướng tâm ($mv^2 w{,xx}$) trong hệ phương trình vi phân chuyển động rời rạc.

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

Trả lời: So với công trình của Songsuwan và cộng sự (2015) [62] (dùng phương pháp Ritz cho dầm 1D-FGM) và Karamanlı (2018) [135] (dùng xấp xỉ đa thức tĩnh cho dầm 2D-FGM hai pha), luận án đã phát triển hai mô hình phần tử hữu hạn hoàn toàn mới: phần tử dầm 2 nút 10-DOF với trường nội suy uốn - trượt tách rời và phần tử tấm tứ giác 9 nút Q9 tích hợp nội suy liên kết. Phương pháp này giải quyết triệt để vấn đề nghẽn trượt cắt và đảm bảo sự hội tụ trơn tru trên toàn dải tỷ số kích thước $L/h, a/h \in [5, 20]$.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Trả lời: Phát hiện mang tính phản trực giác là ảnh hưởng của lực Coriolis ở dải vận tốc cao ($v > 50\text{ m/s}$) không chỉ làm tăng độ võng động cực đại mà còn làm biến dạng phi đối xứng nghiêm trọng biểu đồ ứng suất pháp dọc trục $\sigma_{xx}$ và ứng suất trượt tiếp tuyến $\tau_{xz}$ theo chiều cao dầm, đẩy điểm tập trung ứng suất trượt cực đại lệch về phía biên tiếp xúc của lớp chịu tải.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Trả lời: Có. Luận án trình bày chi tiết toàn bộ biểu thức giải tích ma trận độ cứng phần tử ($\mathbf{k}{aa}, \mathbf{k}{bb}, \mathbf{k}{ss}, \mathbf{k}{ab}, \mathbf{k}{as}, \mathbf{k}{bs}$), ma trận khối lượng phần tử ($\mathbf{m}{uu}, \mathbf{m}{wbwb}, \mathbf{m}{wsws}, \dots$), các véc-tơ lực nút ${\mathbf{f}{ex}}$ và mã nguồn tích phân số Cầu phương Gauss cùng giải thuật Newmark được tài liệu hóa rõ ràng trong các phụ lục kỹ thuật.

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

Trả lời: Lộ trình tập trung vào 4 trọng tâm: (1) Phi tuyến hình học kết hợp tải trọng nhiệt độ ngẫu hợp động; (2) Kết cấu FGM ba pha có lỗ rỗng xốp vi mô chế tạo bằng công nghệ 3D-printing; (3) Kết cấu vỏ thoải và vỏ trụ sandwich 2D-FGM; (4) Tích hợp vật liệu thông minh điều khiển chủ động rung động.

Kết luận

  1. Luận án đã xây dựng thành công mô hình toán học và phần tử hữu hạn 10 bậc tự do có độ chính xác cao cho phân tích dao động dầm sandwich 2D-FGM hai và ba pha dựa trên lý thuyết biến dạng trượt bậc ba cải tiến và lý thuyết biến dạng trượt lượng giác.
  2. Thiết lập phần tử tấm tứ giác 9 nút (Q9) trên cơ sở lý thuyết tấm Mindlin kết hợp phương pháp nội suy liên kết, loại bỏ hoàn toàn hiện tượng nghẽn trượt trong phân tích dao động tự do của tấm sandwich 2D-FGM ba pha.
  3. Giải quyết trọn vẹn bài toán động lực học phi dừng của dầm sandwich 2D-FGM chịu tải trọng di động điều hòa và khối lượng di động có xét tới lực Coriolis, lực quán tính và lực ly tâm.
  4. Chứng minh định lượng sự sai khác có hệ thống giữa các mô hình đồng nhất hóa vi cơ học Voigt, Mori-Tanaka và Maxwell, khẳng định mô hình Voigt đánh giá quá cao độ cứng của kết cấu từ $5%$ đến $15%$.
  5. Khảo sát toàn diện và cung cấp bức tranh định lượng về ảnh hưởng của gradient vật liệu hai chiều ($n_x, n_z$), tỷ số hình học nhịp - chiều dày ($L/h, a/h$), tỷ số cấu trúc lớp (1-1-1, 2-1-2, 1-2-1), vận tốc và tần số kích động tới tần số tự nhiên, dạng dao động, hệ số động lực học $D_d$ và trường ứng suất động học của dầm và tấm sandwich 2D-FGM.