Tổng quan về luận án

Nghiên cứu về kết cấu vật liệu có cơ tính biến đổi (Functionally Graded Materials - FGM) đại diện cho bước chuyển dịch quan trọng trong kỹ thuật kết cấu tiên tiến, nhằm giải quyết triệt để các hạn chế cố hữu của vật liệu composite phân lớp truyền thống như hiện tượng tách lớp (delamination) và tập trung ứng suất tại mặt phân giới. Luận án tiến sĩ "Mô hình phần tử hữu hạn trong phân tích dao động của dầm có cơ tính biến đổi theo hai chiều" do nghiên cứu sinh Trần Thị Thơm thực hiện (chuyên ngành Cơ học vật rắn, mã số 9440107, dưới sự hướng dẫn của PGS. Nguyễn Đình Kiên và PGS. Nguyễn Xuân Thà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, 2019) là công trình tiên phong giải quyết bài toán cơ học phức tạp đối với kết cấu dầm có cơ tính biến đổi đồng thời theo hai chiều không gian (2D-FGM).

Trong kỹ thuật hàng không vũ trụ, năng lượng hạt nhân và thiết bị phản ứng nhiệt, kết cấu thường xuyên chịu các gradient nhiệt độ và tải trọng cơ học biến thiên phức tạp theo cả chiều cao lẫn chiều dài dầm. Tuy nhiên, phần lớn các công trình công bố trước đây chỉ khảo sát dầm FGM một chiều (1D-FGM) với cơ tính thay đổi dọc trục (axially graded) hoặc theo chiều cao (transverse graded) dựa trên hai pha vật liệu cơ bản. Luận án đã xác định khoảng trống nghiên cứu (research gap) trọng yếu: Thiếu hụt các mô hình toán học giải tích tường minh và mô hình số tin cậy để phân tích dao động tự do, dao động cưỡng bức của dầm FGM 4 pha vật liệu (2 gốm, 2 kim loại) có cơ tính biến thiên hai chiều $(x, z)$ theo quy luật hàm lũy thừa kép kết hợp với sự thay đổi tiết diện (dầm thon) và ảnh hưởng phi tuyến của nhiệt độ môi trường.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:

  • RQ1: Làm thế nào để xây dựng mô hình phần tử hữu hạn (PTHH) chính xác cao cho dầm 2D-FGM 4 pha chịu trường nhiệt phụ thuộc nhiệt độ mà không gặp hiện tượng khóa trượt (shear locking)?

  • RQ2: Quy luật phân bố vật liệu hai chiều $(n_x, n_z)$ tác động như thế nào đến các mode dao động riêng, ma trận độ cứng tương hỗ kéo-uốn $(A_{12})$, xoắn-uốn $(A_{44})$ và tần số dao động tự do?

  • RQ3: Dạng thon của thiết diện ngang $(A(x), I(x))$ theo ba quy luật hình học (Dạng A, B, C) tương tác ra sao với gradient vật liệu 2D trong việc cải thiện đáp ứng động lực học?

  • RQ4: Vận tốc của lực di động $(v)$ và gradient nhiệt $(\Delta T)$ ảnh hưởng như thế nào đến tham số động lực học cực đại $\max(D_d)$ và sự phân bố ứng suất pháp $\sigma_{xx}$, ứng suất tiếp $\tau_{xz}$ theo chiều cao dầm?

  • H1: Việc biểu diễn trường chuyển vị theo lý thuyết biến dạng trượt bậc ba cải tiến (ITSDT) sử dụng góc trượt ngang $\gamma_0$ làm biến độc lập sẽ loại bỏ hoàn toàn hệ số điều chỉnh trượt $\psi$ và tăng tốc độ hội tụ số so với cách tiếp cận theo góc quay $\theta$.

  • H2: Gradient vật liệu dọc trục $(n_x)$ kết hợp với gradient theo chiều cao $(n_z)$ tạo ra sự dịch chuyển đáng kể của trục trung hòa vật lý, làm xuất hiện các ghép nối động lực học phi đối xứng giữa dao động dọc trục và dao động uốn.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của mô hình đồng hóa vật liệu Voigt đa pha, lý thuyết biến dạng trượt bậc nhất (FSDT), lý thuyết biến dạng trượt bậc ba cải tiến (ITSDT) của Shi (2007), nguyên lý biến phân Hamilton và thuật toán tích phân thời gian Newmark gia tốc trung bình. Nghiên cứu thực hiện trên phạm vi dầm 2D-FGM tuyến tính đàn hồi, tỷ số độ mảnh $L/h \in [5, 100]$, dầm thon với tham số tiết diện $c \in [0, 0.8]$, chịu sự gia tăng nhiệt độ $\Delta T$ từ $0\text{K}$ đến $80\text{K}$ (nhiệt độ tham chiếu $T_0 = 300\text{K}$) dưới tác dụng của tải trọng di động với dải vận tốc $v$ từ $10\text{m/s}$ đến $150\text{m/s}$.

Literature Review và Positioning

Phân tích lịch sử nghiên cứu cơ học kết cấu FGM cho thấy sự phát triển qua ba giai đoạn chính:

  1. Giai đoạn nền tảng 1D-FGM giải tích: Aydogdu và Taskin (2007) áp dụng phương pháp chuỗi Navier giải bài toán dao động dầm FGM đơn giản; Sankar (2001) đưa ra nghiệm đàn hồi chính xác cho dầm FGM chịu tải hình sin; Şimşek và Kocatürk (2009) sử dụng đa thức Lagrange giải bài toán dầm Euler-Bernoulli FGM chịu tải di động; Wattanasakulpong và cộng sự (2012) áp dụng lý thuyết trượt bậc ba Shi phân tích mất ổn định nhiệt. Hạn chế cốt lõi của các mô hình này là chỉ giải quyết được hệ số phương trình vi phân dạng hằng số trong không gian 1D.
  2. Giai đoạn phương pháp số cho dầm FGM phức tạp: Shahba và cộng sự (2011), Gan và Nguyễn Đình Kiên (2013-2015), Lê Thị Hà (2015), Bùi Văn Tuyển (2016) phát triển các phần tử dầm Timoshenko và dầm thứ bậc phân tích dầm thon 1D-FGM, khung FGM và dầm có lỗ rỗng xốp vi mô. Phương pháp cầu phương vi phân (CPVP) cũng được Malekzadeh (2008), Ebrahimi và cộng sự (2015) triển khai cho bài toán nhiệt.
  3. Giai đoạn tiếp cận 2D-FGM hiện đại: Şimşek (2015, 2016) phân tích dầm 2D-FGM 2 pha chịu tải di động bằng hàm Ritz; Hao và Wei (2016) dùng phương pháp ma trận độ cứng động lực học; Lezgy-Nazargah (2016) sử dụng hàm NURBS đẳng hình học tính toán ứng suất nhiệt; Shafiei và cộng sự (2017) nghiên cứu dầm micro/nano 2D-FGM bằng CPVP.
                  TIẾN TRÌNH NGHIÊN CỨU DẦM FGM VÀ VỊ TRÍ CỦA LUẬN ÁN

Trong bức tranh tổng quan đó, hai cuộc tranh luận học thuật lớn đang diễn ra:

  • Tranh luận về hệ số điều chỉnh trượt $\psi$: Phái ủng hộ FSDT (Timoshenko) duy trì hệ số $\psi = 5/6$, trong khi các nhà nghiên cứu cơ học vật liệu biến tính chỉ ra rằng $\psi$ trong FGM là một hàm phụ thuộc vị trí và gradient vật liệu, không thể là hằng số. Phái lý thuyết bậc cao (TSDT, ITSDT) khẳng định cần triệt tiêu hoàn toàn sự phụ thuộc vào $\psi$ bằng trường biến dạng trượt phân bố bậc hai theo chiều cao dầm.
  • Tranh luận về biến số độc lập trong công thức ITSDT: Việc sử dụng góc quay pháp tuyến $\theta$ đòi hỏi đạo hàm bậc cao $w_{0,xx}$, dẫn đến yêu cầu tính liên tục $C^1$ cho chuyển vị uốn và làm phức tạp ma trận phần tử. Việc thay thế bằng góc trượt ngang $\gamma_0 = \theta + w_{0,x}$ mở ra khả năng đơn giản hóa bậc tự do và kiểm soát trực tiếp năng lượng trượt.

So sánh định vị học thuật với hai nghiên cứu quốc tế tiêu biểu:

  1. So với nghiên cứu của Şimşek (2015, 2016): Şimşek sử dụng phương pháp giải tích Ritz với hàm cơ sở đa thức cho dầm 2D-FGM 2 pha thuần nhất về tiết diện. Luận án vượt trội khi phát triển thành công mô hình phần tử hữu hạn tổng quát cho dầm 4 pha, giải quyết được đồng thời dầm có tiết diện biến thiên (dầm thon) và các điều kiện biên phức tạp mà phương pháp Ritz gặp khó khăn về hội tụ.
  2. So với công trình của Lezgy-Nazargah (2016): Lezgy-Nazargah sử dụng phương pháp đẳng hình học NURBS cho vật liệu biến thiên theo hàm số mũ. Mô hình PTHH của luận án (đặc biệt là phần tử $TBS\gamma$ và $FBHi$) đạt hiệu quả tính toán cao hơn về chi phí ma trận, áp dụng trực tiếp cho hàm lũy thừa kép thực tế hơn trong chế tạo luyện kim bột và tính toán chính xác ứng suất nhiệt động dưới tải trọ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à hoàn thiện lý thuyết dầm đàn hồi phi đồng nhất trên các phương diện:

  1. Mô hình hóa quy luật vật liệu 2D-FGM 4 pha: Thiết lập hàm phân bố tỷ phần thể tích tổng quát cho 4 pha vật liệu gồm 2 gốm ($C_1: \text{Al}2\text{O}3, C_2: \text{ZrO}2$) và 2 kim loại ($M_1: \text{SUS304}, M_2: \text{Al}$ hoặc $\text{Ti-6Al-4V}$): $$V{C1} = \left(\frac{z}{h} + \frac{1}{2}\right)^{n_z} \left(1 - \frac{x}{L}\right)^{n_x}, \quad V{C2} = \left(\frac{z}{h} + \frac{1}{2}\right)^{n_z} \left(\frac{x}{L}\right)^{n_x}$$ $$V{M1} = \left[1 - \left(\frac{z}{h} + \frac{1}{2}\right)^{n_z}\right] \left(1 - \frac{x}{L}\right)^{n_x}, \quad V_{M2} = \left[1 - \left(\frac{z}{h} + \frac{1}{2}\right)^{n_z}\right] \left(\frac{x}{L}\right)^{n_x}$$ Thỏa mãn điều kiện chuẩn hóa $\sum V_i = 1$. Cơ tính hiệu dụng $P(x, z, T)$ được nội suy chính xác qua mô hình Voigt phụ thuộc nhiệt độ phi tuyến: $$P(x, z, T) = [P_{C1}(T) - P_{M1}(T)]\left(\frac{z}{h}+\frac{1}{2}\right)^{n_z}\left(1-\frac{x}{L}\right)^{n_x} + P_{M1}(T)\left(1-\frac{x}{L}\right)^{n_x} + [P_{C2}(T) - P_{M2}(T)]\left(\frac{z}{h}+\frac{1}{2}\right)^{n_z}\left(\frac{x}{L}\right)^{n_x} + P_{M2}(T)\left(\frac{x}{L}\right)^{n_x}$$ với $P(T) = P_0(P_{-1}T^{-1} + 1 + P_1T + P_2T^2 + P_3T^3)$.

  2. Mở rộng lý thuyết biến dạng trượt bậc ba cải tiến (ITSDT): Xây dựng hệ phương trình biến dạng và ứng suất không gian hai chiều, thiết lập ma trận độ cứng độ đàn hồi tổng quát với các thành phần ghép nối bậc cao $(A_{34}, A_{44}, A_{66})$ và ma trận độ cứng trượt bậc cao $(B_{11}, B_{22}, B_{44})$: $$U_B = \frac{1}{2}\int_0^L \left[ A_{11}\varepsilon_m^2 + 2A_{12}\varepsilon_m\varepsilon_b + A_{22}\varepsilon_b^2 - 2A_{34}\varepsilon_m\varepsilon_{hs} - 2A_{44}\varepsilon_b\varepsilon_{hs} + A_{66}\varepsilon_{hs}^2 + \left(\frac{25}{16}B_{11} - \frac{5}{2h^2}B_{22} + \frac{1}{h^4}B_{44}\right)\gamma_0^2 \right] dx$$

                       KHUNG PHÂN TÍCH LÝ THUYẾT VÀ TÍCH HỢP HỆ THỐNG

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

Sự kết hợp giữa cơ học giải tích biến phân và đại số ký hiệu hình thức (Symbolic Computation với phần mềm Maple) cho phép tác giả suy dẫn chính xác các ma trận độ cứng phần tử $[k^e]$, ma trận khối lượng $[m^e]$, và ma trận độ cứng hình học do ứng suất nhiệt $[k_T^e]$ dưới dạng tích phân giải tích từng phần, loại bỏ triệt để sai số tích phân số Gauss khi cơ tính thay đổi liên tục dọc theo chiều dài phần tử.

Các điều kiện biên cơ học tường minh:

  • Ngàm cứng (Clamped - C): $u_0 = w_0 = \theta = 0$ (hoặc $\gamma_0 = w_{0,x}$) tại $x=0, L$.
  • Tựa giản đơn (Simply Supported - S): $u_0 = w_0 = M_{xx} = 0$.
  • Đầu tự do (Free - F): Lực dọc trục $N_{xx} = 0$, mô-men uốn $M_{xx} = 0$, lực cắt $V_x = 0$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ lập trường nhận thức luận thực chứng (positivist paradigm) và phương pháp suy diễn định lượng trong cơ học tính toán (Computational Mechanics). Thiết kế nghiên cứu đa cấp độ (Multi-level modeling framework) được triển khai qua 4 bước:

                      QUY TRÌNH THIẾT KẾ NGHIÊN CỨU VÀ TÍNH TOÁN

Quy trình nghiên cứu và chi tiết 4 mô hình phần tử

Để xử lý triệt để hiện tượng khóa trượt và đảm bảo tính liên tục của trường chuyển vị, luận án đã phát triển 4 mô hình phần tử chuyên biệt:

  1. Mô hình FBKo (FSDT - Hàm dạng Kosmatka): Phần tử dầm 2 nút, mỗi nút 3 bậc tự do $(u_0, w_0, \theta)$. Trường chuyển vị ngang $w_0$ và góc quay $\theta$ được ghép nối thông qua nghiệm giải tích của phương trình tĩnh học dầm Timoshenko chịu uốn thuần túy: $$N_w = \begin{bmatrix} 1 - \frac{x}{l} - \frac{x(l-x)(l-2x)}{(1+\Phi)l^3} & \dots \end{bmatrix}, \quad \Phi = \frac{12 E I}{\psi G A l^2}$$ Cách tiếp cận này loại bỏ hoàn toàn hiện tượng khóa trượt khi dầm trở nên siêu mỏng ($L/h \to 100$).
  2. Mô hình FBHi (FSDT - Hàm dạng thứ bậc): Sử dụng các đa thức trực giao Legendre bậc cao làm hàm dạng bổ sung bên trong phần tử, cho phép tăng bậc nội suy (p-refinement) mà không cần chia nhỏ lưới (h-refinement), tối ưu hóa ma trận độ cứng cho dầm có tiết diện biến thiên $A(x), I(x)$.
  3. Mô hình TBS$\theta$ (ITSDT - Biểu diễn theo $\theta$): Phần tử gồm các bậc tự do $(u_0, w_0, w_{0,x}, \theta)$ tại mỗi nút, sử dụng hàm Hermite bậc ba $C^1$ cho $w_0$ và hàm tuyến tính $C^0$ cho $u_0, \theta$.
  4. Mô hình TBS$\gamma$ (ITSDT - Biểu diễn theo $\gamma_0$): Phần tử 2 nút với các bậc tự do $(u_0, w_0, w_{0,x}, \gamma_0)$. Đây là điểm nhấn đột phá của luận án: việc đưa biến dạng trượt tại mặt giữa $\gamma_0$ làm ẩn số nút độc lập cho phép kiểm soát trực tiếp trường biến dạng trượt thực tế, nâng cao độ chính xác khi phân tích dầm dày ($L/h = 5$).

Dữ liệu vật liệu và công cụ tính toán

Dữ liệu đặc tính cơ - nhiệt của 4 pha vật liệu được trích xuất chính xác từ thực nghiệm quốc tế:

  • Alumina ($\text{Al}_2\text{O}_3$): $E_0 = 349.55\text{ GPa}, \rho = 3800\text{ kg/m}^3, \nu = 0.26, \alpha_0 = 10.00 \times 10^{-6}\text{ K}^{-1}$
  • Zirconia ($\text{ZrO}_2$): $E_0 = 244.27\text{ GPa}, \rho = 3000\text{ kg/m}^3, \nu = 0.288, \alpha_0 = 12.766 \times 10^{-6}\text{ K}^{-1}$
  • Thép không gỉ ($\text{SUS304}$): $E_0 = 201.04\text{ GPa}, \rho = 7800\text{ kg/m}^3, \nu = 0.3, \alpha_0 = 12.33 \times 10^{-6}\text{ K}^{-1}$
  • Hợp kim Titan ($\text{Ti-6Al-4V}$): $E_0 = 122.56\text{ GPa}, \rho = 4429\text{ kg/m}^3, \nu = 0.2888, \alpha_0 = 7.5788 \times 10^{-6}\text{ K}^{-1}$
  • Nhôm ($\text{Al}$): $E_0 = 70\text{ GPa}, \rho = 2707\text{ kg/m}^3, \nu = 0.3$

Hệ thống tính toán tích hợp phần mềm đại số ký hiệu Maple để thiết lập dạng tường minh của ma trận độ cứng và ma trận khối lượng, sau đó chuyển giao sang mã nguồn số học tối ưu hóa trên nền tảng Fortran/MATLAB để giải hệ phương trình động lực học Newmark với bước thời gian $\Delta t = 10^{-4}\text{ s}$.

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

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

  1. Sự bất đối xứng trong tương tác gradient vật liệu 2D: Tần số tự nhiên không thứ nguyên $\mu_1 = \omega_1 L^2 \sqrt{\rho_0 A_0 / E_0 I_0}$ chịu sự chi phối bất đối xứng mạnh mẽ từ cặp tham số $(n_x, n_z)$. Khi tăng số mũ phân bố theo chiều cao $n_z$ từ $0$ lên $5$, tỷ phần kim loại tăng lên làm giảm độ cứng chống uốn hiệu dụng $A_{22}$, dẫn đến $\mu_1$ giảm tới $38.4%$. Ngược lại, biến thiên của tham số dọc trục $n_x$ tạo ra sự tái phân bố khối lượng quán tính $I_{11}(x)$ và độ cứng dọc $A_{11}(x)$, làm thay đổi vị trí các điểm nút và bụng sóng của các mode dao động riêng bậc cao.

  2. Hiện tượng suy giảm tần số do ứng suất nhiệt và nhiệt độ: Sự gia tăng nhiệt độ môi trường $\Delta T$ gây ra tác động kép: làm suy giảm mô-đun đàn hồi $E(x,z,T)$ của vật liệu thành phần và sinh ra lực nén dọc trục nhiệt ban đầu $\sigma_{xx}^T$. Đối với dầm ngàm hai đầu (C-C), khi $\Delta T$ tăng từ $0\text{K}$ lên $80\text{K}$, tham số tần số cơ bản $\mu_1$ sụt giảm nghiêm trọng từ $22.37$ xuống còn $16.12$ (giảm $27.9%$). Nếu tiếp tục tăng nhiệt độ tới ngưỡng tới hạn $T_{cr}$, $\mu_1$ tiến dần về $0$, báo hiệu hiện tượng mất ổn định nhiệt (thermal buckling).

  3. Ưu thế vượt trội của phần tử $TBS\gamma$ trong phân tích dầm dày: Khảo sát đối chuẩn trên dầm tựa đơn (S-S) với tỷ số độ mảnh $L/h = 5$ (dầm có hiệu ứng trượt rất lớn) chỉ ra rằng mô hình cổ điển CBT đánh giá quá cao tần số riêng tới $18.2%$ do bỏ qua trượt. Mô hình $TBS\gamma$ đạt độ chính xác gần như tuyệt đối so với lời giải giải tích đàn hồi 2D chính xác, với sai số dưới $0.15%$ chỉ với $n_E = 20$ phần tử, chứng minh khả năng mô phỏng biến dạng trượt bậc ba vượt trội so với FSDT truyền thống.

  4. Tương tác động lực học phức tạp dưới lực di động: Tham số động lực học $D_d = w_{0,\max}(L/2, t) / w_{st}$ phụ thuộc phi tuyến vào vận tốc di động $v$. Luận án phát hiện hiện tượng cộng hưởng động lực học xuất hiện tại dải vận tốc tới hạn $v_{cr} \approx 80 - 110\text{ m/s}$ tùy thuộc vào dạng thon. Dầm thon Dạng C ($A(x) = A_0(1-cx/L)^2$) cho thấy khả năng giảm độ võng động cực đại tới $24.6%$ so với dầm đều có cùng khối lượng, khẳng định tiềm năng tối ưu hóa hình học kết hợp gradient vật liệu.

            SO SÁNH THAM SỐ TẦN SỐ CƠ BẢN μ1 CỦA DẦM S-S 2D-FGM (L/h = 20, ΔT = 50K)

           ẢNH HƯỞNG CỦA VẬN TỐC LỰC DI ĐỘNG v ĐẾN THAM SỐ ĐỘNG LỰC HỌC Dd (L/h = 20)

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp hệ phương trình chuyển động giải tích hoàn chỉnh cho dầm 2D-FGM 4 pha, làm cơ sở mở rộng cho bài toán tấm và vỏ 2D-FGM.
  • Về mặt phương pháp luận: Khẳng định tính ưu việt của việc sử dụng biến độc lập $\gamma_0$ trong ITSDT, cung cấp một phương pháp luận số học hiệu quả để xây dựng các phần tử hữu hạn bậc cao không bị khóa trượt.
  • Về mặt ứng dụng thực tiễn: Cho phép các kỹ sư hàng không vũ trụ thiết kế các cánh dầm tuabin, vỏ tàu con thoi chịu nhiệt độ cực hạn với khả năng chủ động điều chỉnh gradient $(n_x, n_z)$ để triệt tiêu dao động cộng hưởng và tối ưu hóa khối lượng kết cấu.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu cần tiếp tục hoàn thiện:

  1. Giả thiết tuyến tính hình học: Nghiên cứu giới hạn ở bài toán biến dạng bé, bỏ qua số hạng phi tuyến von Kármán $(\frac{1}{2}w_{0,x}^2)$, do đó chưa mô tả được đáp ứng sau mất ổn định (post-buckling) và dao động biên độ lớn.
  2. Mô hình trường nhiệt độ đồng nhất: Giả thiết nhiệt độ tăng đều $\Delta T$ trên toàn bộ thể tích dầm chưa phản ánh đầy đủ trường nhiệt độ 2D phi tuyến thực tế sinh ra do sự dẫn nhiệt dừng theo phương trình Fourier hai chiều.
  3. Mô hình tải trọng di động giản lược: Tải trọng được mô hình hóa dưới dạng một lực tập trung $P$ di động với vận tốc không đổi, bỏ qua quán tính của vật thể chuyển động, lực ma sát tiếp xúc và hiệu ứng gia tốc/giảm tốc (Coriolis và gia tốc ly tâm).
  4. Bỏ qua khuyết tật vi mô: Chưa tích hợp ảnh hưởng của lỗ xốp vi mô (porosity) hình thành trong quá trình thiêu kết luyện kim bột FGM.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển mô hình phần tử dầm phi tuyến hình học hoàn toàn (Co-rotational hoặc Total Lagrangian) cho dầm 2D-FGM.
  • Giải bài toán truyền nhiệt 2D ổn định và không ổn định để thiết lập trường nhiệt độ thực $T(x,z)$ trước khi phân tích dao động.
  • Mở rộng bài toán tương tác động lực học xe - cầu FGM dưới đoàn tải trọng ngẫu nhiên.
  • Nghiên cứu tối ưu hóa cấu trúc topo (Topology Optimization) phân bố gradient 2D nhằm cực đại hóa tần số riêng và giảm thiểu ứng suất nhiệt cục bộ.

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

  • Tác động học thuật: Các công bố từ luận án trên các tạp chí quốc tế uy tín (ISI/Scopus như Composite Structures, Latin American Journal of Solids and Structures) và tạp chí chuyên ngành trong nước (Tạp chí Cơ học) đã tạo ra một nhánh nghiên cứu mới về kết cấu 2D-FGM tại Việt Nam, thu hút nhiều trích dẫn trong cộng đồng cơ học tính toán.
  • Chuyển giao công nghiệp và R&D: Cung cấp công cụ phần mềm tính toán độc lập, tin cậy cho các viện nghiên cứu thiết kế hàng không vũ trụ (Viện Công nghệ Không gian), ngành công nghiệp quốc phòng và giao thông vận tải đường sắt tốc độ cao.
  • Ảnh hưởng chính sách và đào tạo: Đóng góp tài liệu học thuật chuyên sâu phục vụ công tác đào tạo sau đại học ngành Cơ học vật rắn, Cơ kỹ thuật và Khoa học vật liệu tại Học viện Khoa học và Công nghệ và các trường đại học kỹ thuật hàng đầu.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận xây dựng phần tử hữu hạn nâng cao, kỹ thuật tính toán symbolic kết hợp giải thuật số Newmark.
  • Các nhà khoa học đầu ngành: Kế thừa khung lý thuyết dầm 2D-FGM để mở rộng cho các bài toán phức tạp hơn như tấm 2D-FGM, vỏ thoải, dầm nano có hiệu ứng kích thước (size-dependent effects).
  • Kỹ sư thiết kế kết cấu tiên tiến: Sử dụng các bảng số liệu chuẩn hóa và biểu đồ tham số trong luận án làm tài liệu tra cứu trực tiếp khi tối ưu hóa kết cấu chịu nhiệt độ cao và tải trọng động.

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 độc đáo nhất là việc mở rộng thành công Lý thuyết biến dạng trượt bậc ba cải tiến (ITSDT) của Shi (2007) cho bài toán dầm 2D-FGM 4 pha vật liệu có tiết diện biến thiên. Luận án đã tái cấu trúc lại toàn bộ trường biến dạng trượt bằng cách sử dụng góc trượt ngang $\gamma_0$ làm biến số độc lập, loại bỏ hoàn toàn hệ số điều chỉnh trượt $\psi$, đồng thời suy dẫn chính xác hệ phương trình vi phân chuyển động và biểu thức năng lượng biến dạng đàn hồi $U_B$ chứa đầy đủ các số hạng độ cứng tương hỗ bậc cao $(A_{34}, A_{44}, A_{66}, B_{11}, B_{22}, B_{44})$.

2. Sự đổi mới trong phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đó?

So với nghiên cứu của Şimşek (2015, 2016) (chỉ áp dụng phương pháp Ritz cho dầm 2 pha đều) và Lezgy-Nazargah (2016) (dùng phương pháp NURBS cho dầm hàm mũ), luận án đã:

  • Xây dựng 4 mô hình phần tử hữu hạn đa dạng (FBKo, FBHi, TBS$\theta$, TBS$\gamma$), cho phép rời rạc hóa linh hoạt mọi miền hình học dầm thon (3 dạng thon A, B, C) với điều kiện biên tùy ý.
  • Sử dụng hàm nội suy thứ bậc (FBHi) và hàm Kosmatka (FBKo) khắc phục triệt để hiện tượng khóa trượt số, đạt tốc độ hội tụ nhanh hơn $30%$ về số lượng phần tử cần thiết để đạt cùng độ chính xác.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu số là gì?

Phát hiện bất ngờ nhất là hiện tượng đảo chiều phân bố ứng suất pháp $\sigma_{xx}$ theo chiều cao dầm khi thay đổi tỷ lệ gradient $(n_x, n_z)$ dưới tải trọng di động vận tốc cao ($v = 100\text{ m/s}$). Khác với dầm 1D-FGM hoặc dầm thuần nhất (nơi ứng suất kéo/nén phân bố đối xứng qua trục hình học), dầm 2D-FGM xuất hiện sự dịch chuyển phi tuyến của vị trí ứng suất triệt tiêu dọc theo chiều dài dầm. Tại một số tiết diện gần ngàm, ứng suất pháp cực đại không nằm ở thớ đáy ngoài cùng mà dịch chuyển vào bên trong do tương tác giữa mô-đun đàn hồi cục bộ $E(x,z,T)$ và gradient biến dạng uốn-trượt bậc cao.

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

Có. Luận án cung cấp giao thức tính toán minh bạch và chi tiết:

  • Toàn bộ công thức ma trận phần tử giải tích được dẫn xuất tường minh.
  • Bảng hệ số nhiệt động lực học $P_0, P_{-1}, P_1, P_2, P_3$ cho từng vật liệu ($\text{Al}_2\text{O}_3, \text{ZrO}_2, \text{SUS304}, \text{Ti-6Al-4V}, \text{Al}$).
  • Sơ đồ khối thuật toán phân tích dao động tự do (Hình 3.2) và quy trình tích phân Newmark tính đáp ứng động lực học dưới tải di động (Hình 3.3).

5. Chương trình nghiên cứu 10 năm được định hình như thế nào?

  1. Năm 1-3: Mở rộng mô hình PTHH 2D-FGM sang bài toán phi tuyến hình học lớn (Large Displacement/Rotation) và phân tích sau mất ổn định nhiệt.
  2. Năm 4-6: Tích hợp mô hình dẫn nhiệt 2D Fourier/Non-Fourier và hiệu ứng kích thước micro/nano (Nonlocal Strain Gradient Elasticity) cho các bộ vi cơ điện tử tử (MEMS/NEMS).
  3. Năm 7-10: Phát triển thuật toán học máy (Machine Learning/Physics-Informed Neural Networks - PINNs) kết hợp tối ưu hóa đa mục tiêu topology dầm/tấm 2D-FGM trong các kết cấu hàng không thế hệ mới.

Kết luận

Luận án tiến sĩ của NCS. Trần Thị Thơm là một công trình khoa học mẫu mực, hoàn chỉnh và có tính đột phá cao trong lĩnh vực Cơ học vật rắn và Cơ học tính toán kết cấu FGM:

  1. Đề xuất và thiết lập thành công mô hình dầm 2D-FGM 4 pha vật liệu đầu tiên trong y văn, mô tả chính xác quy luật phân bố hàm lũy thừa kép theo hai hướng không gian kết hợp ảnh hưởng nhiệt độ phi tuyến.
  2. Phát triển 4 mô hình phần tử hữu hạn tiên tiến (FBKo, FBHi, TBS$\theta$, TBS$\gamma$), trong đó mô hình $TBS\gamma$ dựa trên biến độc lập $\gamma_0$ của lý thuyết biến dạng trượt bậc ba cải tiến Shi đại diện cho bước tiến đột phá về độ chính xác và tính ổn định số.
  3. Khảo sát toàn diện ảnh hưởng của hình học dầm thon (3 dạng thon A, B, C), xác định quy luật tương tác phức tạp giữa tham số thiết diện $c$, gradient vật liệu $(n_x, n_z)$ và độ mảnh $L/h$ tới phổ tần số dao động riêng.
  4. Làm sáng tỏ bản chất ứng xử động lực học dưới tải trọ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 sóng ứng suất động trong môi trường nhiệt độ cao.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Cơ học kết cấu 2D-FGM phi tuyến, tối ưu hóa vật liệu chức năng đa hướng và tính toán đa tỷ lệ cho kết cấu nhiệt hàng không vũ trụ.
  6. Giá trị kế thừa vững chắc: Cung cấp bộ công cụ thuật toán số và dữ liệu chuẩn đối chuẩn (benchmark data) có độ tin cậy tuyệt đối, đóng góp quan trọng vào sự phát triển của ngành Cơ học tính toán tại Việt Nam và hội nhập quốc tế.