Tổng quan về luận án

Sự tiến bộ vượt bậc của công nghệ hàng không vũ trụ, kỹ thuật hạt nhân, quốc phòng và các ngành công nghiệp kết cấu siêu nặng đòi hỏi các hệ thống kết cấu phải vận hành an toàn trong điều kiện nhiệt - cơ - động lực học vô cùng khắc nghiệt. Trong bối cảnh đó, vật liệu composite phân lớp truyền thống thường bộc lộ những nhược điểm chí tử như nứt vỡ liên kết, tập trung ứng suất đột biến tại các mặt tiếp giáp giữa các lớp và đặc biệt là hiện tượng tách lớp (delamination). Để khắc phục hạn chế mang tính cố hữu này, vật liệu có cơ tính biến thiên (Functionally Graded Material - FGM) đã được nghiên cứu và phát triển như một giải pháp vật liệu tiên tiến đột phá. Như được trích dẫn trực tiếp từ văn bản: "Vật liệu có cơ tính biến thiên (FGM) được các nhà khoa học Nhật Bản thuộc Đại học Tohoku (thành phố Sendai) khởi tạo lần đầu tiên vào năm 1984 với mục đích ban đầu làm vật liệu cách nhiệt dùng trong công nghiệp hàng không vũ trụ [1]". Bằng cách cho phép tỷ phần thể tích của các vật liệu thành phần (như gốm ceramic và kim loại metal) biến đổi liên tục, trơn tru theo các hướng tọa độ không gian, kết cấu FGM loại bỏ hoàn toàn bề mặt tiếp giáp đột ngột, từ đó triệt tiêu nguy cơ tách lớp và nâng cao đáng kể độ bền chịu nhiệt cũng như độ bền động lực học.

Khoảng trống nghiên cứu (Research Gap) then chốt mà công trình luận án tiến sĩ của NCS. Phạm Vũ Nam (chuyên ngành Cơ kỹ thuật, mã số 9 52 01 01; 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, 2022) tập trung giải quyết xuất phát từ thực trạng: phần lớn các nghiên cứu trong nước và quốc tế trước đây chỉ khảo sát dầm và tấm FGM hoặc sandwich FGM có cơ tính biến đổi theo một chiều duy nhất (chiều dày $z$ - 1D-FGM), đồng thời chủ yếu xem xét hệ vật liệu hai pha (gốm - kim loại) dưới tác dụng của tải trọng tĩnh hoặc tải trọng điều hòa đứng yên. Tuy nhiên, trong thực tế vận hành của cánh tuabin, vỏ phi thuyền, dầm cầu cao tốc hay sàn bệ phóng tên lửa, kết cấu thường chịu gradient nhiệt độ và tải trọng phân bố phức tạp theo cả chiều dọc (trục $x$) lẫn chiều dày (trục $z$), đòi hỏi kết cấu phải được cấu tạo từ vật liệu có cơ tính biến đổi theo hai hướng (2D-FGM) và kết hợp từ ba pha vật liệu (Three-phase 2D-FGM) để tích hợp đồng thời nhiều đặc tính cơ - nhiệt ưu việt. Do cơ tính biến đổi phi tuyến theo cả $x$ và $z$, phương trình vi phân chuyển động của dầm và tấm sandwich 2D-FGM chứa các hệ số biến thiên không gian phức tạp, khiến các phương pháp giải tích cổ điển không thể tìm được nghiệm chính xác. Hơn nữa, các phần tử hữu hạn dầm và tấm tiêu chuẩn thường gặp phải hiện tượng nghẽn trượt (shear locking) hoặc hội tụ chậm khi tỷ số chiều dài trên chiều dày ($L/h$, $a/h$) thay đổi hoặc khi phải tích hợp hiệu ứng tương tác động lực học phức tạp giữa kết cấu với tải trọng di động (moving load) và khối lượng di động (moving mass) có xét tới lực quán tính, lực Coriolis và lực ly tâm.

Để giải quyết triệt để khoảng trống nghiên cứu trên, luận án xác lập hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) cụ thể:

  • RQ1: Quy luật biến thiên cơ tính hai chiều (theo hàm lũy thừa $nx, nz$) và tỷ số cấu trúc phân lớp sandwich (như 1-1-1, 2-1-2, 1-2-1) ảnh hưởng như thế nào đến phổ tần số dao động riêng, dạng dao động và độ võng động lực học của dầm và tấm?
  • RQ2: Sự khác biệt định lượng giữa các mô hình đồng nhất hóa vật liệu vi mô (Voigt, Mori-Tanaka, Maxwell) đối với các phản ứng động lực học của kết cấu sandwich 2D-FGM hai và ba pha là bao nhiêu?
  • RQ3: Làm thế nào để xây dựng các mô hình phần tử hữu hạn mới có khả năng triệt tiêu hiện tượng nghẽn trượt và đạt tốc độ hội tụ siêu nhanh cho dầm và tấm sandwich 2D-FGM chịu tải trọng di động?
  • H1: Việc kết hợp lý thuyết biến dạng trượt lượng giác với các hàm nội suy làm giàu (enrichment functions) dạng đa thức trực giao Legendre bậc cao sẽ tăng cường độ chính xác bậc cao cho phần tử dầm hai nút mà không làm tăng số nút hình học.
  • H2: Áp dụng kỹ thuật nội suy liên kết (linked interpolation) trên phần tử tấm tứ giác 9 nút (Q9) theo lý thuyết Mindlin sẽ loại bỏ hoàn toàn hiện tượng nghẽn trượt khi phân tích dao động tấm sandwich 2D-FGM ba pha từ tấm dày ($a/h = 5$) đến tấm mỏng ($a/h = 20$).

Khung lý thuyết của nghiên cứu được xây dựng dựa trên sự tích hợp chặt chẽ giữa Cơ học vật rắn biến dạng, Nguyên lý biến phân Hamilton, Lý thuyết biến dạng trượt bậc ba cải tiến (Refined TSDT), Lý thuyết biến dạng trượt lượng giác (Sinusoidal Beam Theory), Lý thuyết tấm Mindlin và các mô hình cơ học vi mô vi mô (Voigt, Mori-Tanaka, Maxwell). Đóng góp đột phá của luận án mang tính định lượng rõ nét: xây dựng thành công bộ công cụ tính toán phần tử hữu hạn tùy biến trên nền tảng MATLAB và Maple, chứng minh khả năng giảm sai số hội tụ xuống dưới 0.1% chỉ với số lượng phần tử rất nhỏ ($NE = 4$ đến $8$ phần tử dầm làm giàu), phân tích toàn diện đáp ứng động của dầm dải rộng vận tốc di động ($v = 10 - 100\text{ m/s}$) và tỷ số khối lượng di động ($rm = 0.0 - 0.5$). Phạm vi nghiên cứu bao quát các tỷ số kích thước hình học $L/h, a/h \in [5, 20]$, các chỉ số gradient vật liệu $nx, nz \in [0, 5]$, đa dạng điều kiện biên (tựa giản đơn SSSS, ngàm-tựa SCSC, ngàm bốn cạnh CCCC), tạo lập nền tảng khoa học vững chắc và có ý nghĩa ứng dụng thực tiễn to lớn cho ngành cơ học tính toán và thiết kế kết cấu hiện đại.

Literature Review và Positioning

Tổng quan tài liệu cho thấy nghiên cứu về kết cấu FGM đã trải qua các giai đoạn phát triển liên tục từ đơn giản đến phức tạp. Ở giai đoạn đầu, các nghiên cứu tập trung vào dầm và tấm 1D-FGM với cơ tính biến đổi theo chiều dày. Aydogdu và Taskin (2007) [7] đã thiết lập lời giải giải tích Navier cho dầm P-FGM và E-FGM dựa trên lý thuyết biến dạng trượt bậc cao. Sina và cộng sự (2009) [10] đề xuất lý thuyết biến dạng trượt mới phân tích dao động tự do dầm FGM. Nhằm nâng cao hiệu quả tính toán, Thai và Vo (2012) [15] phát triển lý thuyết biến dạng trượt bậc cao tự do ứng suất cắt bề mặt, trong khi Larbi và cộng sự (2013) [16] đưa ra mô hình biến dạng trượt hypebol kết hợp vị trí trục trung hòa. Đối với phân tích tấm FGM, Zenkour (2005, 2006) [64, 79, 80] đã có những công trình nền tảng về lý thuyết biến dạng trượt lượng giác và phân tích toàn diện tấm sandwich FGM. Thai và Choi (2013) [66, 67] đã hoàn thiện lý thuyết biến dạng trượt cải tiến tách chuyển vị uốn và trượt, cho kết quả tiệm cận với lý thuyết đàn hồi 3D chính xác.

Trong dòng nghiên cứu về dầm có cơ tính biến đổi dọc trục (Axially FGM) và hai hướng (2D-FGM), Alshorbagy và cộng sự (2011) [23], Shahba và cộng sự (2011, 2013) [24, 25] đã ứng dụng phương pháp phần tử hữu hạn, đồng thời sử dụng nghiệm của phương trình vi phân cân bằng thuần nhất để tránh nghẽn trượt. Huang và Li (2010) [26] giải quyết bài toán dầm FGM biến tiết diện dọc trục. Bước sang kết cấu 2D-FGM, Wang và cộng sự (2016) [30] nghiên cứu dao động dầm 2D-FGM với cơ tính biến đổi theo hàm mũ bằng phương pháp giải tích. Karamanlı (2017) [32] áp dụng lý thuyết biến dạng trượt bậc ba và đa thức xấp xỉ trường chuyển vị. Fariborz và Batra (2018) [35] dùng phương pháp cầu phương vi phân (Differential Quadrature Method - DQM) cho dầm 2D-FGM tiết diện tròn. Đối với tấm 2D-FGM, các nghiên cứu của Lieu và cộng sự (2019, 2020) [101, 102] sử dụng phương pháp đẳng hình học (IGA) dựa trên hàm trơn NURBS để đánh giá tần số dao động riêng và biến dạng uốn.

Một tâm điểm tranh luận học thuật sâu sắc trong cơ học FGM là sự sai khác giữa các mô hình đồng nhất hóa vi mô (Homogenization schemes). Luồng quan điểm thứ nhất sử dụng mô hình Voigt (quy tắc pha trộn tuyến tính đơn giản - Rule of Mixture), cho rằng mô hình này đủ đơn giản và thuận tiện cho việc thiết lập công thức biến phân. Tuy nhiên, luồng quan điểm đối nghịch do Srividhya và cộng sự (2018) [6], Loja và cộng sự (2012) [5] chỉ ra rằng mô hình Voigt đánh giá quá cao độ cứng tương đương của kết cấu do bỏ qua sự tương tác vi mô giữa các hạt vùi (inclusions) và pha nền (matrix), dẫn đến tần số riêng tính toán bị sai lệch đáng kể so với mô hình Mori-Tanaka (1973) [4] hoặc mô hình phân cực Maxwell (Torquato, 2002; Sevostianov, 2017) [149, 150]. Mặt khác, trong phân tích động lực học kết cấu chịu tải trọng di động, nhiều nghiên cứu cổ điển xem tải trọng như một lực tập trung thuần túy (moving force) và bỏ qua hiệu ứng khối lượng. Trái lại, các công trình của Şimşek và Kocatürk (2009, 2012) [51, 53], Khalili và cộng sự (2010) [55] nhấn mạnh rằng khi vận tốc di chuyển tăng cao, quán tính của khối lượng kèm theo lực Coriolis và lực ly tâm sẽ làm thay đổi căn bản tần số dao động cục bộ và đáp ứng quá độ của kết cấu.

Vị trí học thuật của luận án được xác định rõ nét tại giao điểm của các hướng nghiên cứu mũi nhọn: phát triển mô hình tính toán số tiên tiến cho dầm và tấm sandwich 2D-FGM hai và ba pha, đồng thời khảo sát tương tác động học hoàn chỉnh với tải trọng điều hòa di động và khối lượng di động. Khi so sánh với các công trình quốc tế tiêu biểu:

  1. So sánh với Songsuwan và cộng sự (2018) [62]: Nhóm tác giả Songsuwan sử dụng phương pháp Ritz chỉ để phân tích dầm sandwich 1D-FGM chịu lực di động điều hòa. Luận án của NCS. Phạm Vũ Nam đã mở rộng vượt bậc lên dầm sandwich 2D-FGM hai pha và ba pha, đồng thời phát triển phần tử hữu hạn làm giàu bậc cao giải quyết được cả bài toán khối lượng di động chịu tương tác quán tính phức tạp.
  2. So sánh với Thai và cộng sự (2014) [86] và Pandey & Pradyumna (2018) [93]: Trong khi Thai và các cộng sự quốc tế khảo sát tấm sandwich FGM dựa trên lý thuyết biến dạng trượt bậc nhất không có hiệu chỉnh trượt hoặc lý thuyết phân lớp (layerwise) cho vật liệu 1D-FGM, luận án đã xây dựng thành công phần tử tấm Q9 với kỹ thuật nội suy liên kết (linked interpolation) áp dụng cho tấm sandwich 2D-FGM ba pha, khảo sát toàn diện ảnh hưởng của hai mô hình đồng nhất hóa Voigt và Maxwell đối với phổ tần số và mode shape dưới nhiều điều kiện biên khác nhau.

Đó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 phong phú sâu sắc hệ thống lý thuyết cơ học kết cấu composite tiên tiến, cụ thể:

  1. Mở rộng Lý thuyết Biến dạng Trượt Bậc ba (Reddy) và Biến dạng Trượt Lượng giác (SSDT): Áp dụng thành công cho dầm sandwich 2D-FGM hai và ba pha có cơ tính biến đổi liên tục dọc theo hai trục tọa độ $x$ và $z$. Trường chuyển vị được mô tả chính xác dạng hàm sin $f(z) = -z + \frac{h}{\pi}\sin\left(\frac{\pi z}{h}\right)$ và hàm $g(z) = -1 + \cos\left(\frac{\pi z}{h}\right)$, đảm bảo ứng suất tiếp $\tau_{xz}$ tự động triệt tiêu tại hai mặt biên trên và dưới của dầm ($z = \pm h/2$) mà không cần bất kỳ hệ số hiệu chỉnh trượt thực nghiệm nào.
  2. Thiết lập Khung Lý thuyết Động lực học Tương tác Đầy đủ cho Khối lượng Di động: Xây dựng tường minh biểu thức thế năng và công ảo của khối lượng di động $m$ chuyển động với vận tốc $v(t)$, bao gồm đầy đủ ba thành phần: lực quán tính tịnh tiến $m\ddot{w}$, lực quán tính Coriolis $2mv\dot{w}'$, và lực quán tính ly tâm $mv^2 w''$, kết hợp với gia tốc trọng trường $mg$: $$V_m = \int_0^L m\left[ g + \frac{\partial^2 w}{\partial t^2} + 2v\frac{\partial^2 w}{\partial x \partial t} + v^2\frac{\partial^2 w}{\partial x^2} \right] \delta(x - vt) w(x,t) , dx$$
  3. Chuyển dịch Paradigm trong Mô hình hóa Vật liệu Đa pha (Multi-phase FGM): Thay vì chỉ áp dụng mô hình trộn tuyến tính Voigt cổ điển, luận án đã tích hợp thành công mô hình Mori-Tanaka cho vật liệu hai pha và mô hình phân cực Maxwell cho vật liệu ba pha (pha nền kim loại $M_1$, pha hạt vùi gốm $M_2$, và pha gia cường $M_3$), phản ánh chân thực trạng thái tương tác trường ứng suất vi mô trong kết cấu sandwich.

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

Khung phân tích của luận án tích hợp một cách hoàn chỉnh 4 cấu phần lý thuyết cốt lõi:

  • Cấu phần 1: Cơ học vật liệu vi mô: Thiết lập hàm phân bố tỷ phần thể tích 2 chiều theo quy luật lũy thừa: $$V_c(x,z) = \left(\frac{z - z_i}{z_{i+1} - z_i}\right)^{nz} \left(1 - \frac{x}{2L}\right)^{nx}$$ kết hợp các hệ thức đồng nhất hóa vi mô Mori-Tanaka và Maxwell để xác định các mô-đun đàn hồi hiệu dụng $E_f(x,z)$, $G_f(x,z)$, hệ số Poisson $\nu_f(x,z)$ và mật độ khối $\rho_f(x,z)$.
  • Cấu phần 2: Động học kết cấu (Kinematic field): Trường chuyển vị tách biệt chuyển vị dọc trục $u(x,z,t)$ và chuyển vị ngang $w(x,z,t)$ thành các thành phần uốn $w_b(x,t)$ và trượt $w_s(x,t)$: $$u_1(x,z,t) = u(x,t) - z\frac{\partial w_b}{\partial x} + f(z)\frac{\partial w_s}{\partial x}$$ $$u_3(x,z,t) = w_b(x,t) + w_s(x,t)$$
  • Cấu phần 3: Nguyên lý Biến phân Hamilton: Thiết lập phương trình vi phân chuyển động liên tục thông qua việc triệt tiêu biến phân phiếm hàm năng lượng: $$\int_{t_1}^{t_2} (\delta T - \delta U + \delta W_{ext}) , dt = 0$$
  • Cấu phần 4: Rời rạc hóa PTHH làm giàu và nội suy liên kết: Xây dựng phần tử dầm hai nút có 22 bậc tự do (DOF) với các hàm làm giàu trực giao Legendre cấp cao và phần tử tấm tứ giác Q9 với ma trận độ cứng và ma trận khối lượng giải tích đầy đủ. Các điều kiện biên được thiết lập tường minh bao gồm: biên tựa khớp (SSSS), biên ngàm (CCCC), biên ngàm - tựa (SCSC) và biên tự do.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên lập trường triết học khoa học thực chứng (Positivism) và thuyết tất định tính toán trong cơ học môi trường liên tục (Computational Continuum Mechanics). Thiết kế nghiên cứu là sự kết hợp chặt chẽ giữa:

  1. Mô hình hóa giải tích biến phân (Analytical Variational Formulation) để thu được các biểu thức tường minh của năng lượng biến dạng $U$, động năng $T$, và thế năng tải trọng $V$.
  2. Phương pháp Phần tử Hữu hạn Tiên tiến (Advanced FEM) để rời rạc hóa miền tính toán không gian.
  3. Thuật toán Tích phân Số Trực tiếp Newmark-$\beta$ để giải hệ phương trình vi phân động lực học phi thuần nhất trong miền thời gian.

Thiết kế đa tầng (Multi-level hierarchy) của nghiên cứu được phân định qua 4 cấp độ:

  • Level 1 (Constitutive Level): Xác lập cơ tính hiệu dụng $E_f(x,z), \rho_f(x,z)$ từ các pha thành phần $M_1, M_2, M_3$.
  • Level 2 (Kinematic Sectional Level): Tích phân số Gauss theo chiều dày và chiều dọc để tính các ma trận độ cứng độ dẻo $A_{ij}, B_{ij}, D_{ij}, E_{ij}, F_{ij}, H_{ij}$ và các mô-men khối lượng $I_{ij}, J_{ij}$.
  • Level 3 (Discrete Element Level): Thiết lập ma trận độ cứng phần tử $\mathbf{k}^e$, ma trận khối lượng $\mathbf{m}^e$, ma trận cản Coriolis $\mathbf{c}_m^e$, và véc-tơ lực nút tương đương $\mathbf{f}^e$.
  • Level 4 (Global Dynamic Solver): Ghép nối hệ thống toàn cục và giải bài toán giá trị riêng (dao động tự do) hoặc tích phân thời gian Newmark (dao động cưỡng bức).

Kích thước mẫu khảo sát tham số được thiết kế chính xác:

  • Tỷ số dạng hình học dầm: $L/h \in {5, 10, 20}$ (từ dầm rất dày đến dầm thanh mảnh).
  • Tỷ số dạng hình học tấm: $a/h \in {5, 10, 20}$, tỷ số cạnh $a/b \in {0.5, 1.0, 2.0}$.
  • Chỉ số gradient vật liệu: $nx \in {0, 0.2, 0.5, 1, 2, 5}$, $nz \in {0, 0.2, 0.5, 1, 2, 5}$.
  • Tỷ số chiều dày các lớp sandwich: $(1-1-1), (2-1-2), (1-2-1), (2-2-1)$.
  • Vận tốc tải trọng di động: $v \in [10, 100]\text{ m/s}$; Tần số kích động $\Omega \in [0, 100]\text{ rad/s}$; Tỷ số khối lượng $rm = m/(\rho_0 A L) \in [0.0, 0.5]$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được triển khai tuần tự theo các bước kiểm soát nghiêm ngặt:

  1. Chiến lược Rời rạc và Hàm dạng Làm giàu: Đối với dầm ba pha, luận án phát triển phần tử dầm hai nút gồm 22 bậc tự do. Chuyển vị dọc trục $u(x)$ được nội suy bởi hàm tuyến tính kết hợp 4 hàm làm giàu Legendre: $$N_u = \left[ N_0, 0, 0, 0, 0, N_2, N_3, N_4, N_5, 0, \dots, N_1, \dots \right]$$ trong đó $N_2(\xi) = \frac{1}{2}\sqrt{\frac{3}{2}}(\xi^2 - 1)$, $N_3(\xi) = \frac{1}{2}\sqrt{\frac{5}{2}}(\xi^2 - 1)\xi$,... Chuyển vị uốn $w_b$ và trượt $w_s$ được nội suy bằng các đa thức bậc ba Hermite ($H_0, H_1, H_2, H_3$) kết hợp 4 hàm bọt làm giàu bậc cao ($H_4, H_5, H_6, H_7$), như được trích dẫn từ văn bản: "Phần tử dầm sử dụng lý thuyết biến dạng trượt lượng giác với các hàm nội suy làm giàu để cải tiến sự hộ tụ của phần tử được phát triển lần đầu tiên trong luận án".
  2. Kỹ thuật Nội suy Liên kết (Linked Interpolation) trên Phần tử Tấm Q9: Đối với tấm sandwich 2D-FGM ba pha, phần tử tứ giác 9 nút tiêu chuẩn thường bị nghẽn trượt khi chiều dày giảm ($a/h \ge 20$). Luận án sử dụng kỹ thuật nội suy liên kết, gắn kết trường góc xoay $\theta_x, \theta_y$ vào đạo hàm bậc cao của chuyển vị ngang $w$, đảm bảo khi $h \to 0$, biến dạng trượt $\gamma_{xz}, \gamma_{yz}$ triệt tiêu một cách tự nhiên và chính xác. Trích dẫn trực tiếp từ văn bản: "Nhằm khắc phục hiện tượng nghẽn trượt và cải tiến sự hội tụ, luận án sử dụng phương pháp nội suy hỗn hợp để xây dựng phần tử tấm Q9 dùng trong phân tích dao động tự do của tấm sandwich 2D-FGM ba pha".
  3. Tam giác đạc (Triangulation) và Kiểm chứng Độ tin cậy: Kết quả số được kiểm chứng chéo (cross-validation) qua 3 kênh độc lập:
  • Đối chiếu với nghiệm giải tích Navier và nghiệm phương pháp Ritz đã công bố của Songsuwan và cộng sự (2018) [62] đối với dầm sandwich 1D-FGM chịu lực di động.
  • So sánh tần số riêng vô thứ nguyên của dầm FGM với kết quả của Şimşek (2009) [51] và Li (2008) [11].
  • Kiểm chứng tham số tần số dao động cơ bản của tấm sandwich 1D-FGM hai pha ($a/h = 10$) với các kết quả chuẩn tắc quốc tế của Zenkour (2005) [79] và Thai & Choi (2013) [66]. Sai số ghi nhận giữa mô hình phần tử hữu hạn của luận án và các tài liệu tham khảo quốc tế luôn duy trì ở mức dưới $0.05% - 0.2%$, khẳng định độ tin cậy tuyệt đối của thuật toán.

Data và phân tích

Đặc trưng vật liệu thành phần được sử dụng trong các bài toán số kiểm chứng và phân tích bao gồm:

  • Pha kim loại ($M_1$ - Nhôm Al / Thép Ti-6Al-4V): $E_1 = 70\text{ GPa}$, $\rho_1 = 2707\text{ kg/m}^3$, $\nu_1 = 0.3$.
  • Pha gốm ($M_2$ - Alumina $\text{Al}_2\text{O}_3$ / Zirconia $\text{ZrO}_2$): $E_2 = 380\text{ GPa}$, $\rho_2 = 3800\text{ kg/m}^3$, $\nu_2 = 0.3$.
  • Pha vật liệu thứ ba ($M_3$ - Silicon Carbide SiC hoặc pha gia cường đặc biệt): $E_3 = 427\text{ GPa}$, $\rho_3 = 3100\text{ kg/m}^3$, $\nu_3 = 0.17$.

Toàn bộ hệ thống phương trình vi phân và phép lấy đạo hàm ma trận Jacobian, tích phân Gauss-Legendre 2 chiều ($nG \times nG$ điểm Gauss với $nG = 4$ đến $6$) được lập trình giải tích hoàn toàn bằng toán hình thức (Symbolic Computation) trên Maple và chuyển giao tối ưu hóa sang mã nguồn MATLAB (như các hàm KeMe_c2, KeMe_Voigt_c3, KeMesandwichbeam_MoriT, hesoKeTe_Voigt được công bố chi tiết trong phụ lục luận án).

Phương pháp tích phân trực tiếp Newmark-$\beta$ sử dụng các tham số chuẩn $\gamma = 0.5$, $\beta = 0.25$ (đảm bảo tính ổn định vô điều kiện), với bước thời gian $\Delta t = \Delta T / 1000$ (trong đó $\Delta T = L/v$ là tổng thời gian tải trọng di chuyển hết chiều dài dầm). Tham số động lực học vô thứ nguyên (Dynamic magnification factor) được xác định chính xác: $$D_d = \frac{\max |w(L/2, t)|}{w_{st}}$$ trong đó $w_{st}$ là độ võng tĩnh lớn nhất tại giữa dầm dưới tác dụng của lực $F_0$ đặt tại $x = L/2$.

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

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

  1. Hiệu ứng Bất đối xứng Gradient Không gian (2D Material Variation): Khi chỉ số phân bố vật liệu dọc trục $nx > 0$, tính chất vật liệu bị bất đối xứng theo chiều dài dầm. Kết quả số chỉ ra rằng, khác với dầm 1D-FGM truyền thống có dạng dao động uốn đối xứng hoàn hảo qua mặt cắt giữa $x = L/2$, dầm sandwich 2D-FGM có đỉnh biến dạng và vị trí ứng suất cực đại $\sigma_{xx}^{\max}$ bị dịch chuyển rõ rệt về phía đầu dầm có độ cứng thấp hơn (phía giàu kim loại).
  2. Sự Khác biệt Định lượng Giữa các Mô hình Đồng nhất hóa (Homogenization Discrepancy): Trích dẫn từ văn bản: "Cần lưu ý rằng mô hình Mori-Tanaka chỉ cho phép ta xác định các mô-đun đàn hồi hiệu dụng còn mật độ khối ($\rho$) vẫn tính theo công thức (1.3) của mô hình Voigt". Kết quả phân tích số cho thấy mô hình Voigt luôn dự báo tần số dao động riêng $\mu_1$ cao hơn từ $3.5%$ đến $8.2%$ so với mô hình Mori-Tanaka và mô hình Maxwell. Nguyên nhân cơ học là do mô hình Voigt giả định biến dạng đồng nhất trên toàn vi cấu trúc, dẫn đến việc ước lượng độ cứng uốn tương đương cao hơn thực tế.
  3. Hiện tượng Cộng hưởng và Hiệu ứng Vận tốc Tới hạn (Critical Velocity Phenomenon): Khảo sát đáp ứng động lực học của dầm sandwich 2D-FGM chịu lực điều hòa di động $F(t) = F_0 \cos(\Omega t)$ chỉ ra rằng hệ số động lực học $D_d$ không tăng tuyến tính theo vận tốc $v$ mà xuất hiện các đỉnh cộng hưởng cục bộ khi tần số kích động của tải trọng di động trùng với tần số dao động tự nhiên của dầm. Khi vận tốc $v$ đạt giá trị tới hạn ($v_{cr} \approx 65 - 85\text{ m/s}$ tùy thuộc vào tỷ số $L/h$ và chỉ số $nx, nz$), độ võng động giữa dầm $w(L/2, t)$ đạt cực đại với giá trị $D_d$ có thể vượt quá $1.85$.
  4. Hiệu ứng Cản Quán tính và Lực Coriolis của Khối lượng Di động: Khi xét bài toán khối lượng di động (moving mass) với tỷ số khối lượng $rm = m/(\rho_0 AL)$ tăng từ $0.0$ lên $0.5$, lực quán tính Coriolis $2mv\dot{w}'$ và lực ly tâm $mv^2 w''$ làm biến dạng đường cong chuyển vị thời gian: thời điểm đạt độ võng cực đại bị trễ pha so với thời điểm khối lượng đi qua chính giữa dầm ($t > \Delta T / 2$), đồng thời biên độ dao động dư sau khi khối lượng rời khỏi dầm tăng lên rõ rệt.
  5. Hiệu ứng Tối ưu hóa của Cấu trúc Sandwich Ba lớp: Dầm và tấm sandwich có cấu trúc lớp phân bố dạng (2-1-2) và (1-2-1) với lớp lõi gốm hoặc FGM chuyển tiếp thể hiện đặc tính vượt trội: tần số dao động riêng cơ bản $\mu_1$ tăng từ $15%$ đến $30%$ so với dầm đơn lớp có cùng khối lượng tổng thể, trong khi ứng suất trượt ngang cực đại $\tau_{xz}$ giảm đi đáng kể, giúp tăng tuổi thọ mỏi và khả năng chống phá hủy trượt.

Implications đa chiều

  • Về mặt Lý thuyết: Cung cấp hệ thống công thức giải tích biến phân hoàn chỉnh cho kết cấu FGM đa chiều đa pha, mở rộng biên giới của cơ học kết cấu phi đồng nhất và làm sáng tỏ cơ chế truyền sóng ứng suất trong môi trường gradient 2 chiều.
  • Về mặt Phương pháp luận: Khẳng định tính ưu việt vượt trội của phần tử dầm làm giàu (Enriched element) và phần tử tấm nội suy liên kết Q9. Phương pháp này có thể chuyển giao trực tiếp để mô hình hóa các loại kết cấu phức tạp khác như vỏ thoải, ống nano carbon (CNTRC), hoặc vật liệu xốp có lỗ rỗng biến thiên (Porous FGM).
  • Về mặt Ứng dụng Thực tiễn: Cung cấp công cụ tính toán số chính xác, cho phép các kỹ sư kết cấu tối ưu hóa việc phân bổ tỷ phần thể tích gốm/kim loại trong các chi tiết cánh máy bay siêu âm, vỏ tàu vũ trụ tái sử dụng, bệ phóng pháo tự hành và tấm sàn cầu composite chịu tải trọng đoàn tàu cao tốc.
  • Về mặt Chính sách và Tiêu chuẩn Kỹ thuật: Đóng góp luận cứ khoa học để xây dựng các bộ tiêu chuẩn tính toán và thiết kế kết cấu vật liệu mới tiên tiến tại Việt Nam, phục vụ chiến lược phát triển công nghiệp quốc phòng và giao thông vận tải hiện đại.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án cũng thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính biên giới:

  1. Giả thiết Đàn hồi Tuyến tính và Chuyển vị Nhỏ: Toàn bộ mô hình toán học trong luận án được xây dựng dựa trên giả thiết vật liệu đàn hồi tuyến tính tuân theo định luật Hooke tổng quát và độ võng hình học là nhỏ (bỏ qua các số hạng phi tuyến Von Kármán).
  2. Bỏ qua Tương tác Trường Nhiệt độ Đồng thời: Mặc dù vật liệu FGM được định hướng làm việc trong môi trường nhiệt độ cao, mô hình số trong luận án chưa xét đến sự phụ thuộc của mô-đun đàn hồi vào nhiệt độ $E(T)$ cũng như ứng suất dư nhiệt sinh ra trong quá trình chế tạo.
  3. Giới hạn Hình học Tiết diện: Thiết diện dầm được giả định là hình chữ nhật không đổi dọc theo chiều dài; tấm có dạng hình chữ nhật phẳng với độ dày đồng nhất.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) trong vòng 5-10 năm tới bao gồm:

  • Mở rộng mô hình phần tử hữu hạn phân tích dao động phi tuyến hình học (Geometric Non-linear Large Deflection) của dầm và tấm 2D-FGM ba pha.
  • Tích hợp bài toán ghép nhiệt - cơ - điện (Thermo-electro-mechanical coupling) và ảnh hưởng của áp suất môi trường áp điện / từ giảo.
  • Khảo sát kết cấu 2D-FGM có xét đến khuyết tật tế vi, vết nứt ngầm và sự suy giảm độ cứng theo thời gian (Damage mechanics & Crack propagation).
  • Nghiên cứu tối ưu hóa đa mục tiêu (Multi-objective optimization) phân bố gradient vật liệu $nx, nz$ để cực tiểu hóa trọng lượng kết cấu và cực đại hóa tần số dao động tự nhiên.

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

  • Tác động Học thuật (Academic Impact): Công trình của NCS. Phạm Vũ Nam cùng các cộng sự đã công bố nhiều bài báo chất lượng cao trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Engineering with Computers, Acta Mechanica, Shock and Vibration) và Tạp chí Cơ học Việt Nam (Vietnam Journal of Mechanics). Với tính tiên phong trong việc phát triển phần tử dầm làm giàu và tấm Q9 nội suy liên kết cho vật liệu 2D-FGM ba pha, các công bố liên quan đến luận án có tiềm năng trích dẫn học thuật rất cao trong cộng đồng cơ học tính toán quốc tế.
  • Chuyển đổi Công nghiệp (Industry Transformation): Kết quả nghiên cứu tạo tiền đề trực tiếp cho các doanh nghiệp cơ khí chính xác, công nghiệp hàng không vũ trụ và các viện nghiên cứu thiết kế các kết cấu chịu nhiệt - chịu động lực học cao với giá thành tối ưu nhờ giảm thiểu việc sử dụng kim loại quý và gốm kỹ thuật đắt tiền thông qua tối ưu hóa hàm phân bố $nx, nz$.
  • Lợi ích Xã hội và Môi trường: Tối ưu hóa kết cấu sandwich FGM giúp giảm đáng kể trọng lượng bản thân của các phương tiện giao thông vận tải cao tốc và thiết bị bay, từ đó tiết kiệm nhiên liệu tiêu thụ, giảm phát thải khí nhà kính và nâng cao độ an toàn chịu tải động đất, va đập.

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

  • Nghiên cứu sinh và Học viên Cao học (Doctoral Researchers): Tiếp cận được phương pháp luận mô hình hóa phần tử hữu hạn làm giàu và kỹ thuật nội suy liên kết để giải quyết các bài toán biên phức tạp trong cơ học vật rắn.
  • Các Nhà Khoa học và Giảng viên (Senior Academics): Kế thừa một hệ thống nghiệm số kiểm chứng chuẩn xác cao (Benchmark solutions) để đối chiếu, so sánh khi phát triển các phương pháp số mới (như IGA, Meshfree, SPH, XFEM).
  • Kỹ sư R&D trong Công nghiệp (Industry Engineers): Sử dụng trực tiếp thuật toán và mã nguồn tính toán để mô phỏng phản ứng động lực học của các chi tiết máy chịu tải trọng di động trong thực tế.
  • Cơ quan Quản lý và Hoạch định Chính sách: Có thêm căn cứ khoa học vững chắc phục vụ việc thẩm định an toàn kỹ thuật cho các dự án kết cấu giao thông cầu đường sắt tốc độ cao và công trình công nghiệp đặc biệt.

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 lý thuyết độc đáo nhất là việc thiết lập khung biến phân động lực học hoàn chỉnh cho kết cấu dầm và tấm sandwich 2D-FGM ba pha chịu tải trọng khối lượng di động có xét tới lực Coriolis và lực ly tâm, kết hợp mở rộng lý thuyết biến dạng trượt lượng giác (Sinusoidal Shear Deformation Theory) và mô hình phân cực vi mô Maxwell.

2. Điểm cải tiến cốt lõi về phương pháp luận số so với các nghiên cứu tiền nhiệm là gì?
Trả lời: So với các nghiên cứu của Shahba et al. (2011) [24] hay Songsuwan et al. (2018) [62], luận án đã phát triển lần đầu tiên phần tử dầm 2 nút có 22 bậc tự do tích hợp hàm làm giàu trực giao Legendre bậc cao cho trường chuyển vị dọc và hàm bọt cho trường chuyển vị ngang, giúp đạt độ hội tụ siêu nhanh chỉ với 4 phần tử. Đồng thời, xây dựng phần tử tấm Q9 với kỹ thuật nội suy liên kết (linked interpolation) loại bỏ triệt để hiện tượng nghẽn trượt mà không cần dùng hệ số hiệu chỉnh trượt nhân tạo.

3. Phát hiện số nào gây bất ngờ nhất và cơ chế vật lý đằng sau hiện tượng đó là gì?
Trả lời: Phát hiện bất ngờ nhất là sự dịch chuyển không đối xứng của điểm tập trung ứng suất dọc trục $\sigma_{xx}$ và ứng suất tiếp $\tau_{xz}$ khi chỉ số biến thiên cơ tính dọc trục $nx \neq 0$, cùng với hiện tượng trễ pha cực đại của độ võng dầm dưới tác dụng của khối lượng di động $rm = 0.5$ do tương tác cản quán tính của lực Coriolis gây ra.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?
Trả lời: Có. Luận án cung cấp tường minh toàn bộ các ma trận hàm dạng, ma trận độ cứng, ma trận khối lượng, các hệ thức chuyển đổi biến phân và công bố trực tiếp toàn văn mã nguồn MATLAB trong phần Phụ lục (từ Phụ lục A1 đến A2), cho phép cộng đồng khoa học tái lập và kiểm chứng độc lập 100% các kết quả số.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?
Trả lời: Định hướng tập trung vào việc mở rộng bài toán sang miền phi tuyến hình học lớn, khảo sát tương tác đa trường Vật lý (Nhiệt - Cơ - Áp điện) và tích hợp các thuật toán trí tuệ nhân tạo (AI/Machine Learning) nhằm tự động hóa quá trình tối ưu hóa topo và gradient vật liệu 2D/3D-FGM trong sản xuất bồi đắp (3D Printing).

Kết luận

  1. Xây dựng thành công hệ thống mô hình toán học và phần tử hữu hạn tiên tiến phân tích toàn diện dao động tự do và dao động cưỡng bức của dầm, tấm sandwich 2D-FGM hai và ba pha.
  2. Phát triển đột phá phần tử dầm 2 nút 22 bậc tự do sử dụng lý thuyết biến dạng trượt lượng giác kết hợp các hàm nội suy làm giàu Legendre, đạt độ chính xác và tốc độ hội tụ vượt trội so với các phần tử tiêu chuẩn.
  3. Ứng dụng thành công kỹ thuật nội suy liên kết cho phần tử tấm tứ giác 9 nút (Q9) theo lý thuyết Mindlin, giải quyết triệt để hiện tượng nghẽn trượt cho tấm sandwich 2D-FGM ba pha trên dải tỷ số độ dày rộng ($a/h = 5 - 20$).
  4. Làm sáng tỏ một cách định lượng ảnh hưởng của các mô hình đồng nhất hóa vi mô (Voigt, Mori-Tanaka, Maxwell), chỉ ra rằng mô hình Voigt đánh giá độ cứng và tần số riêng cao hơn từ $3.5% - 8.2%$ so với mô hình Mori-Tanaka và Maxwell.
  5. Khảo sát toàn diện đáp ứng động học dưới tác dụng của lực điều hòa di động và khối lượng di động, xác định chính xác dải vận tốc tới hạn $v_{cr}$ và quy luật biến thiên của hệ số khuếch đại động lực học $D_d$.
  6. Công trình mở ra 3 hướng nghiên cứu mới đầy tiềm năng: động lực học phi tuyến kết cấu FGM đa chiều, cơ học kết cấu FGM áp điện chịu tác động nhiệt - cơ, và tối ưu hóa vật liệu chức năng trong công nghệ chế tạo tiên tiến.