Tổng quan luận án

Vật liệu có cơ tính biến thiên (Functionally Graded Material - FGM) là thế hệ vật liệu composite cải tiến, sở hữu các đặc tính cơ - lý thay đổi liên tục và trơn theo không gian. Nhờ loại bỏ các mặt phân chia pha đột ngột, 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 cục bộ và sự tách lớp. Khi kết hợp cấu trúc sandwich với vật liệu FGM, kết cấu đạt được tỷ số độ cứng trên khối lượng cao, khả năng chịu nhiệt và tối ưu hóa chi phí chế tạo.

Tuy nhiên, phần lớn các nghiên cứu trước đây chủ yếu tập trung vào kết cấu FGM biến đổi tính chất theo một chiều (chiều dày - 1D-FGM) và cấu tạo từ hai pha vật liệu gốm - kim loại. Nhằm đáp ứng yêu cầu chịu tải trọng cơ - nhiệt đa chiều phức tạp trong các môi trường kỹ thuật khắc nghiệt, các kết cấu có cơ tính biến đổi theo cả chiều dài và chiều dày (2D-FGM), cũng như được cấu tạo từ ba pha vật liệu thành phần (three-phase FGM) đã được phát triển. Do sự biến thiên cơ tính theo hai phương làm các hệ số trong phương trình vi phân chuyển động không còn là hằng số, các phương pháp giải tích truyền thống gặp nhiều trở ngại trong việc tìm nghiệm chính xác. Do đó, việc xây dựng các mô hình số dựa trên phương pháp phần tử hữu hạn (PTHH) là đòi hỏi cấp thiết nhằm phân tích toàn diện ứng xử động lực học của loại kết cấu này.

Mục tiêu nghiên cứu của luận án được xác định cụ thể như sau:

  1. Xây dựng các mô hình phần tử hữu hạn, thuật toán và chương trình tính toán số phục vụ phân tích dao động tự do và dao động cưỡng bức của dầm sandwich 2D-FGM hai pha và ba pha chịu tác động của tải trọng di động.
  2. Xây dựng mô hình phần tử hữu hạn tấm tứ giác 9 nút cải tiến nhằm phân tích bài toán dao động tự do của tấm sandwich 2D-FGM ba pha với các điều kiện biên khác nhau.
  3. Khảo sát, đánh giá định lượng ảnh hưởng của các mô hình đồng nhất hóa vật liệu (mô hình Voigt, mô hình Mori-Tanaka và mô hình Maxwell) cùng các tham số phân bố vật liệu, tỷ số hình học và tham số tải trọng đến tần số dao động riêng, đáp ứng động lực học và trường ứng suất của kết cấu dầm, tấm sandwich 2D-FGM.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Kết cấu dầm và tấm sandwich 2D-FGM hai pha và ba pha có tiết diện ngang đồng nhất, với tỷ phần thể tích các pha vật liệu thành phần biến thiên liên tục theo quy luật hàm lũy thừa (P-FGM) theo cả chiều dài ($x$) và chiều cao ($z$).
  • Phạm vi nghiên cứu:
    • Bài toán dao động tự do và dao động cưỡng bức của dầm sandwich 2D-FGM hai pha chịu lực điều hòa di động.
    • Bài toán dao động tự do và dao động cưỡng bức của dầm sandwich 2D-FGM ba pha chịu khối lượng di động (xét đến lực quán tính, lực Coriolis và lực ly tâm).
    • Bài toán dao động tự do của tấm sandwich 2D-FGM ba pha hình chữ nhật với các điều kiện biên tựa đơn (SSSS), ngàm - tựa xen kẽ (SCSC) và ngàm bốn cạnh (CCCC).
    • Giới hạn ứng xử cơ học trong miền đàn hồi tuyến tính, biến dạng nhỏ và chuyển vị nhỏ; không xét ảnh hưởng của nhiệt độ và độ ẩm môi trường đến tính chất vật liệu.
Ký hiệu thông số Ý nghĩa vật lý / hình học Ký hiệu thông số Ý nghĩa vật lý / hình học
$L, h, b$ Chiều dài, chiều cao, chiều rộng dầm $a, b, h$ Chiều dài, chiều rộng, chiều dày tấm
$n_x, n_z$ Chỉ số mũ phân bố vật liệu theo trục $x, z$ $E_f, G_f, \rho_f, \nu_f$ Mô-đun đàn hồi, mô-đun trượt, mật độ, hệ số Poisson hiệu dụng
$F_0, \Omega, v$ Biên độ, tần số, vận tốc lực di động $m, r_m$ Khối lượng hạt di động, tỷ số khối lượng ($m / \rho_0 A L$)
$\mu_i, \omega_i$ Tham số tần số thứ $i$, tần số dao động riêng $D_d$ Hệ số động lực học (Dynamic factor)

Tổng quan tài liệu và vị trí của luận án

Vật liệu 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, Nhật Bản [1]. Kể từ đó, việc mô hình hóa và phân tích ứng xử động lực học của kết cấu FGM đã thu hút sự chú ý của nhiều nhóm nghiên cứu trên thế giới và trong nước.

Các hướng nghiên cứu trước đây được luận án điểm luận bao gồm:

  • Dao động dầm FGM biến đổi theo chiều cao (1D-FGM): Aydogdu và Taskin (2007) [7] sử dụng nghiệm Navier phân tích dầm P-FGM và E-FGM; Sina và cộng sự (2009) [10] đề xuất lý thuyết biến dạng trượt mới; Li (2008) [11] giải bài toán giải tích cho dầm có cơ tính biến đổi tùy ý; Ying và cộng sự (2008) [12] dùng phương pháp không gian trạng thái; Giunta và cộng sự (2011) [13] áp dụng mô hình dầm thứ bậc; Su và Banerjee (2012) [14] dùng ma trận độ cứng động học; Thai và Vo (2012) [15], Larbi và cộng sự (2013) [16], Pradhan và Chakraverty (2014, 2015) [17, 18] phát triển các lý thuyết biến dạng trượt bậc cao không cần hệ số hiệu chỉnh trượt; Mahi và cộng sự (2010) [19], Wattanasakulpong và cộng sự (2011) [20] phân tích ảnh hưởng của nhiệt độ đến dao động dầm FGM.
  • Dầm FGM có cơ tính biến đổi dọc trục và dầm 2D-FGM: Alshorbagy và cộng sự (2011) [23], Shahba và cộng sự (2011) [24, 25] sử dụng PTHH với các hàm dạng rút ra từ nghiệm phương trình cân bằng; Huang và Li (2010) [26], Zhao và cộng sự (2012) [28] sử dụng đa thức Chebyshev cho dầm tiết diện thay đổi; Wang và cộng sự (2016) [30], Hao và Wei (2016) [31], Karamanlı (2017) [32] phân tích dầm 2D-FGM biến thiên theo hàm số mũ; Huynh và cộng sự (2018) [33] sử dụng phương pháp đẳng hình học (IGA); Rajasekaran và Khaniki (2018) [34] phân tích dầm 2D-FGM có vết nứt; Fariborz và Batra (2018) [35] dùng phương pháp cầu phương vi phân (CPVP).
  • Dầm sandwich FGM và dầm chịu tải trọng di động: Sayyad và Ghugal (2019) [37] tổng quan toàn diện về dầm sandwich FGM; Chakraborty và cộng sự (2003) [38], Bhangale và Ganesan (2006) [39], Pradhan và Murmu (2009) [40], Su và cộng sự (2014) [41], Osofero và cộng sự (2015) [46], Trinh và cộng sự (2016) [48], Kahya và Turan (2018) [49] nghiên cứu dao động tự do dầm sandwich FGM. Về bài toán tải trọng di động: Şimşek và Kocatürk (2009) [51], Şimşek (2010, 2012) [52, 53, 54], Khalili và cộng sự (2010) [55], Gan và cộng sự (2015) [57], Songsuwan và cộng sự (2018) [62], Wang và Wu (2018) [63] khảo sát đáp ứng dầm FGM chịu lực và khối lượng di động.
  • Dao động tấm FGM, tấm sandwich FGM và tấm 2D-FGM: Zenkour (2005, 2006) [64, 79], Thai và Choi (2013) [66], Matsunaga (2008) [70], Iurlaro và cộng sự (2014) [92], Pandey và Pradyumna (2018) [93]; Lieu và cộng sự (2019) [101, 102], Li và cộng sự (2020) [103] phân tích tấm 2D-FGM bằng phương pháp IGA; Esmaeilzadeh và Kadkhodayan (2019) [104] dùng phương pháp cân bằng động.
  • Các nghiên cứu trong nước: Luận án của Lê Thị Hà (2015) [105], Bùi Văn Tuyển (2018) [106], Trần Thị Thơm (2019) [107], Nguyễn Ngọc Huyên (2017) [108], Ngô Trọng Đức (2018) [109], Nguyễn Bá Duy (2019) [110], Lê Thị Ngọc Ánh (2022) [111]; các công bố của nhóm Nguyễn Trung Kiên và cộng sự (2013–2017) [112–116], Nguyễn Tiến Khiêm và cộng sự (2017, 2020) [117, 118], Trần Văn Liên và cộng sự (2017, 2019) [119, 120], Nguyễn Đình Kiên và cộng sự (2020, 2021) [148, 153].

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Mặc dù số lượng công trình về FGM rất lớn, các nghiên cứu về kết cấu sandwich FGM hầu như chỉ dừng lại ở bài toán biến thiên 1 chiều theo độ dày (1D-FGM) với hai pha vật liệu. Chưa có công trình nào thiết lập mô hình tính toán số hoàn chỉnh cho dầm và tấm sandwich 2D-FGM cấu tạo từ ba pha vật liệu thành phần, đồng thời xem xét đầy đủ các hiệu ứng động lực học của tải trọng di động (lực điều hòa di động, khối lượng di động có kể đến quán tính, lực Coriolis và lực ly tâm). Hơn nữa, việc phát triển các phần tử dầm làm giàu (enriched element) và phần tử tấm tứ giác Q9 sử dụng kỹ thuật nội suy liên kết (linked interpolation) nhằm triệt tiêu hiện tượng nghẽn trượt cho kết cấu sandwich 2D-FGM ba pha là hướng đi chưa được giải quyết trước đó.

Cơ sở lý thuyết và phương pháp nghiên cứu

Mô hình vật liệu và quy luật phân bố 2D-FGM

Luận án xem xét kết cấu sandwich P-FGM gồm ba lớp. Tỷ phần thể tích của các pha vật liệu ($V_i$) biến thiên đồng thời theo cả trục dọc $x$ và trục đứng $z$ theo hàm lũy thừa: $$V(x, z) = V_z(z) \cdot \left(1 - \frac{x}{2L}\right)^{n_x}$$ Trong đó $V_z(z)$ là quy luật biến thiên theo chiều cao tại từng lớp dầm/tấm với chỉ số mũ $n_z$.

Để xác định các tính chất cơ học hiệu dụng ($P_f$) của vật liệu từ các pha thành phần ($P_i$), luận án sử dụng ba mô hình đồng nhất hóa:

  1. Mô hình Voigt: $$P_f(x, z) = \sum_{i=1}^{k} P_i V_i(x, z)$$
  2. Mô hình Mori-Tanaka: Áp dụng cho FGM hai pha để xác định mô-đun đàn hồi hiệu dụng $E_f$ và hệ số Poisson $\nu_f$ thông qua mô-đun khối hiệu dụng $K_f$ và mô-đun trượt hiệu dụng $G_f$: $$E_f = \frac{9K_f G_f}{3K_f + G_f}, \quad \nu_f = \frac{3K_f - 2G_f}{6K_f + 2G_f}$$
  3. Mô hình Maxwell: Áp dụng để đánh giá các đặc trưng đàn hồi của hệ vật liệu FGM ba pha cấu tạo từ vật liệu nền và các hạt độn phân tán.

Lý thuyết biến dạng và thiết lập phương trình chuyển động

  • Dầm hai pha (Chương 2): Áp dụng lý thuyết biến dạng trượt bậc ba cải tiến (Reddy TSDT). Trường chuyển vị được biểu diễn: $$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)$$ Trong đó $w_b$ và $w_s$ lần lượt là thành phần độ võng do uốn và do trượt, hàm hình dáng trượt $f(z) = \frac{4z^3}{3h^2}$.
  • Dầm ba pha (Chương 3): Áp dụng lý thuyết biến dạng trượt lượng giác (Sinusoidal Shear Deformation Theory - SSDT) với các hàm phân bố trượt: $$f(z) = -z + \frac{h}{\pi}\sin\left(\frac{\pi z}{h}\right), \quad g(z) = -1 + \cos\left(\frac{\pi z}{h}\right)$$ Thế năng của khối lượng di động $m$ chuyển động với vận tốc $v$ được xây dựng có tính toán đầy đủ gia tốc toàn phần: $$a_z = \left(\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)_{x=vt}$$
  • Tấm ba pha (Chương 4): Áp dụng lý thuyết biến dạng trượt bậc nhất Mindlin (FSDT) với trường chuyển vị: $$u_1(x, y, z, t) = u(x, y, t) + z \theta_x(x, y, t)$$ $$u_2(x, y, z, t) = v(x, y, t) + z \theta_y(x, y, t)$$ $$u_3(x, y, z, t) = w(x, y, t)$$
  • Nguyên lý biến phân: Phương trình vi phân chuyển động tổng thể được thiết lập thông qua nguyên lý biến phân Hamilton: $$\int_{t_1}^{t_2} (\delta T - \delta U + \delta W) dt = 0$$

Kỹ thuật phần tử hữu hạn và thuật toán số

  • Phần tử dầm làm giàu (Enriched Beam Element): Phát triển phần tử dầm 2 nút với 22 bậc tự do, sử dụng các đa thức làm giàu bổ sung kết hợp đa thức Hermite để xấp xỉ chuyển vị dọc $u$, độ võng uốn $w_b$ và độ võng trượt $w_s$, giúp tăng tốc độ hội tụ số.
  • Phần tử tấm Q9 nội suy liên kết (Linked Interpolation): Xây dựng phần tử tấm tứ giác 9 nút (Q9) sử dụng trường nội suy liên kết giữa chuyển vị ngang $w$ và các góc xoay $\theta_x, \theta_y$ để loại bỏ hiện tượng khóa trượt (shear locking) khi tấm có tỷ số độ mảnh lớn ($a/h$ cao).
  • Thuật toán tích phân: Sử dụng thuật toán Newmark trực tiếp theo thời gian để giải phương trình động lực học dạng ma trận: $$\mathbf{M}\ddot{\mathbf{D}} + \mathbf{C}\dot{\mathbf{D}} + \mathbf{K}\mathbf{D} = \mathbf{F}(t)$$ Trong bài toán khối lượng di động, các ma trận khối lượng bổ sung $\mathbf{M}_m$, ma trận cản Coriolis $\mathbf{C}_m$ và ma trận độ cứng bổ sung $\mathbf{K}_m$ phụ thuộc vào vị trí tức thời của hạt tải trọng.
  • Công cụ thực thi: Toàn bộ thuật toán được lập trình tính toán giải tích symbolic và số trên nền tảng phần mềm Matlab và Maple.

Nội dung chính theo từng chương

Chương 1. Tổng quan

Chương 1 hệ thống hóa cơ sở khoa học về vật liệu FGM, các quy luật biến thiên hàm lũy thừa (P-FGM), hàm số mũ (E-FGM) và Sigmoid (S-FGM). Chương này phân tích chi tiết các mô hình đồng nhất hóa (Voigt, Mori-Tanaka, Maxwell) và tổng thuật các công trình nghiên cứu về dao động tự do, dao động cưỡng bức của dầm, tấm FGM 1D và 2D chịu các loại tải trọng tĩnh và động. Qua phân tích tổng quan các công trình trong và ngoài nước, chương 1 xác định rõ tính cấp thiết, khoảng trống học thuật và định hướng nghiên cứu trọng tâm của luận án.

Chương 2. Dầm sandwich 2D-FGM hai pha

Chương 2 nghiên cứu bài toán dao động tự do và dao động cưỡng bức của dầm sandwich 2D-FGM hai pha (gốm $\text{Al}_2\text{O}_3$ và kim loại $\text{Al}$) chịu lực điều hòa di động $F(t) = F_0 \cos(\Omega t)$. Mô hình toán học được xây dựng dựa trên lý thuyết biến dạng trượt bậc ba cải tiến.

Các kết quả và số liệu chính:

  • Kiểm chứng độ chính xác: So sánh độ võng tĩnh không thứ nguyên $w^*$ với lời giải giải tích và so sánh tần số cơ bản với công bố của Songsuwan và cộng sự (2018) [62], cho thấy sự trùng khít cao.
Thông số dầm ($L/h=10$) $n_x = 0, n_z = 0.5$ $n_x = 0.5, n_z = 0.5$ $n_x = 1.0, n_z = 0.5$ $n_x = 2.0, n_z = 0.5$
Tần số cơ bản $\mu_1$ (Mô hình Voigt) 2.8164 2.9542 3.0651 3.1985
Tần số cơ bản $\mu_1$ (Mô hình Mori-Tanaka) 2.7632 2.8981 3.0074 3.1382
Hệ số động lực học $D_d$ ($v=50\text{ m/s}, \Omega=0$) 1.1845 1.1523 1.1284 1.0965
  • Phát hiện quy luật cơ học:
    • Tần số dao động riêng tăng khi chỉ số vật liệu $n_x$ tăng (tăng hàm lượng pha gốm dọc trục) và giảm khi chỉ số $n_z$ tăng.
    • Tỷ số độ mảnh $L/h$ tăng từ 5 đến 20 làm tăng tần số dao động không thứ nguyên và làm tăng biên độ dao động động lực học tại giữa dầm.
    • Mô hình Voigt luôn cho giá trị tần số dao động riêng cao hơn mô hình Mori-Tanaka do mô hình Voigt đánh giá độ cứng hiệu dụng của kết cấu cao hơn.
    • Ứng suất pháp dọc trục $\sigma_{xx}$ và ứng suất trượt $\tau_{xz}$ phân bố phi tuyến mạnh theo chiều cao dầm và phụ thuộc rõ rệt vào tần số kích động $\Omega$.

Chương 3. Dầm sandwich 2D-FGM ba pha

Chương 3 mở rộng bài toán cho dầm sandwich 2D-FGM ba pha (gốm $\text{ZrO}_2$, gốm $\text{Al}_2\text{O}_3$ và kim loại $\text{Ti-6Al-4V}$) chịu khối lượng di động $m$. Mô hình dầm sử dụng lý thuyết biến dạng trượt lượng giác và phần tử dầm 2 nút làm giàu với 22 bậc tự do.

Các kết quả và số liệu chính:

  • Kiểm tra độ hội tụ của phần tử làm giàu: Đánh giá số phần tử rời rạc dầm ($NE$) cho thấy phần tử làm giàu hội tụ rất nhanh: chỉ cần $NE = 6 \div 8$ phần tử đã đạt độ chính xác tương đương phần tử chuẩn với $NE = 30 \div 40$.
Số phần tử ($NE$) Tham số tần số $\mu_1$ ($L/h=20$) Tham số động lực học $D_d$ ($r_m=0.5, v=50\text{ m/s}$)
2 2.9642 1.2845
4 2.9581 1.2612
6 2.9573 1.2584
8 2.9572 1.2580
16 2.9572 1.2580
  • Quy luật động lực học khối lượng di động:
    • Hiệu ứng khối lượng ($r_m$) làm giảm tần số dao động tức thời của hệ kết cấu và làm tăng mạnh độ võng động lực học cực đại $w(L/2, t)$ so với trường hợp lực di động thuần túy.
    • Khi vận tốc di động $v$ tăng, hệ số động lực học $D_d$ tăng đạt cực đại tại vùng vận tốc cộng hưởng, sau đó có xu hướng giảm.
    • So sánh hai mô hình đồng nhất hóa: Mô hình Voigt cho tần số riêng cao hơn mô hình Maxwell, đồng thời dẫn đến hệ số đáp ứng động lực học $D_d$ nhỏ hơn.

Chương 4. Dao động tự do của tấm sandwich 2D-FGM ba pha

Chương 4 thiết lập bài toán dao động tự do cho tấm chữ nhật sandwich 2D-FGM ba pha dựa trên lý thuyết biến dạng trượt bậc nhất Mindlin và phần tử tấm tứ giác 9 nút (Q9) nội suy liên kết.

Các kết quả và số liệu chính:

  • Kiểm chứng và loại bỏ nghẽn trượt: Phần tử Q9 nội suy liên kết cho kết quả tần số dao động riêng ổn định và chính xác cao ở cả miền tấm dày ($a/h = 5$) lẫn miền tấm rất mỏng ($a/h = 100$), loại bỏ hoàn toàn hiện tượng nghẽn trượt số.
  • Khảo sát tham số tần số riêng $\mu_1$ với các điều kiện biên:
Tỷ số $a/h$ Chỉ số $n_x = n_z$ Biên SSSS (Voigt / Maxwell) Biên SCSC (Voigt / Maxwell) Biên CCCC (Voigt / Maxwell)
5 0.5 0.0542 / 0.0518 0.0685 / 0.0654 0.0812 / 0.0776
5 1.0 0.0512 / 0.0489 0.0648 / 0.0619 0.0769 / 0.0735
20 0.5 0.0588 / 0.0562 0.0894 / 0.0855 0.1145 / 0.1095
20 1.0 0.0556 / 0.0531 0.0846 / 0.0809 0.1084 / 0.1037
  • Quy luật ứng xử của tấm:
    • Tần số dao động riêng tăng dần theo mức độ kiềm chế chuyển vị của điều kiện biên: $\mu_{\text{SSSS}} < \mu_{\text{SCSC}} < \mu_{\text{CCCC}}$.
    • Sự thay đổi tỷ số chiều dày giữa các lớp sandwich (ví dụ cấu trúc 1-1-1, 2-1-2, 1-2-1) ảnh hưởng trực tiếp đến độ cứng chống uốn tổng thể của tấm, trong đó cấu trúc lớp lõi dày hơn mang lại tần số riêng cao hơn khi chịu liên kết ngàm.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và phương pháp

  1. Xây dựng thành công các mô hình phần tử hữu hạn dầm 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 dùng cho phân tích dao động tự do và dao động cưỡng bức của dầm sandwich 2D-FGM hai pha và ba pha chịu tải trọng di động.
  2. Lần đầu tiên phát triển mô hình phần tử dầm 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 (22 bậc tự do), cho phép tăng tốc độ hội tụ số vượt bậc, giảm thiểu số lượng phần tử cần thiết trong mô phỏng động lực học dầm FGM chịu khối lượng di động có kể đến đầy đủ lực Coriolis và lực ly tâm.
  3. Xây dựng phần tử tấm tứ giác 9 nút (Q9) trên cơ sở lý thuyết Mindlin kết hợp kỹ thuật nội suy liên kết (linked interpolation) cho tấm sandwich 2D-FGM ba pha, khắc phục 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.
  4. Phân tích, so sánh một cách có hệ thống ảnh hưởng của các mô hình đồng nhất hóa vật liệu khác nhau (Voigt so với Mori-Tanaka đối với dầm hai pha; Voigt so với Maxwell đối với dầm và tấm ba pha) đến tần số dao động tự nhiên và đáp ứng động lực học của kết cấu sandwich 2D-FGM.

Đóng góp thực tiễn và kiến nghị kỹ thuật

  • Cung cấp bộ chương trình tính toán số tin cậy (trên nền tảng Matlab và Maple) phục vụ phân tích tĩnh và động lực học kết cấu composite FGM biến thiên đa chiều.
  • Thiết lập tập dữ liệu số chuẩn xác và hệ thống đồ thị quan hệ giữa các tham số phân bố vật liệu ($n_x, n_z$), tỷ số kích thước hình học ($L/h, a/h$), tham số tải trọng ($v, \Omega, r_m$) với các đặc trưng dao động, cung cấp căn cứ kỹ thuật trực tiếp cho việc thiết kế, tối ưu hóa độ bền và độ cứng của các kết cấu dạng dầm, tấm sandwich trong công trình cơ khí và công trình giao thông.

Hạn chế và hướng nghiên cứu tiếp

Dựa trên phạm vi và các giả thiết được xác định trong văn bản luận án, các hạn chế và định hướng phát triển tiếp theo bao gồm:

Hạn chế của nghiên cứu

  • Luận án giới hạn bài toán trong phạm vi lý thuyết đàn hồi tuyến tính, biến dạng nhỏ và chuyển vị nhỏ; chưa khảo sát các hiệu ứng phi tuyến hình học (chuyển vị lớn) hoặc phi tuyến vật liệu (ứng xử đàn - dẻo).
  • Bỏ qua ảnh hưởng của trường nhiệt độ và độ ẩm môi trường lên sự suy giảm cơ tính vật liệu cũng như trường ứng suất nhiệt dư trong kết cấu.
  • Chưa xét tới sự xuất hiện của các khuyết tật chế tạo như vết nứt tế vi hoặc độ xốp/lỗ rỗng vi mô (porosities) trong các pha vật liệu FGM.
  • Các kết quả nghiên cứu hoàn toàn dựa trên mô phỏng số và kiểm chứng chéo với các công bố giải tích/PTHH khác; chưa có điều kiện thực nghiệm kiểm chứng trên mẫu chế tạo thực tế.

Hướng nghiên cứu tiếp theo

  • Mở rộng mô hình tính toán cho bài toán dao động dầm, tấm sandwich 2D-FGM trong môi trường nhiệt - ẩm có xét đến sự phụ thuộc của tính chất vật liệu vào nhiệt độ.
  • Phát triển mô hình phần tử hữu hạn phi tuyến hình học để phân tích ứng xử động học phi tuyến và mất ổn định động của kết cấu.
  • Tích hợp ảnh hưởng của độ xốp ngẫu nhiên, vết nứt cục bộ và nghiên cứu bài toán chẩn đoán khuyết tật kết cấu dựa trên các dạng dao động.
  • Nghiên cứu ứng dụng các lớp vật liệu thông minh (như vật liệu áp điện - piezoelectric) gắn trên bề mặt kết cấu sandwich 2D-FGM để điều khiển chủ động dao động.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo trong lĩnh vực Cơ học vật rắn biến dạng, Cơ kỹ thuật và Kỹ thuật Cơ khí:

  • Đối với nghiên cứu sinh và học viên cao học: Luận án cung cấp quy trình toán học chặt chẽ từ việc áp dụng nguyên lý biến phân Hamilton, các lý thuyết biến dạng trượt bậc cao đến kỹ thuật lập trình phần tử hữu hạn nâng cao (phần tử làm giàu, phần tử nội suy liên kết Q9).
  • Đối với giảng viên và nhà nghiên cứu: Cung cấp phương pháp luận phân tích vật liệu composite FGM biến thiên hai chiều và mô hình hóa tương tác động lực học kết cấu - tải trọng di động (lực quán tính, Coriolis, ly tâm).
  • Đối với kỹ sư thiết kế công trình và cơ khí: Cung cấp số liệu định lượng về ảnh hưởng của quy luật phân bố vật liệu, giúp định hướng thiết kế tối ưu hóa tỷ trọng và khả năng chịu tải trọng động của các cấu kiện dạng tấm, dầm composite chịu tải di động tốc độ cao.
  • Phần đáng tham khảo nhất:
    • Kỹ thuật xây dựng ma trận độ cứng, ma trận khối lượng và ma trận quán tính di động trong Chương 2 và Chương 3.
    • Thuật toán phần tử Q9 nội suy liên kết trong Chương 4.
    • Mã nguồn lập trình các hàm tính toán phần tử trong Phụ lục A1 và A2.

Câu hỏi thường gặp

1. Tại sao luận án tập trung nghiên cứu vật liệu 2D-FGM và cấu trúc sandwich thay vì vật liệu 1D-FGM truyền thống?

Vật liệu 1D-FGM chỉ biến đổi tính chất cơ lý theo chiều dày, phù hợp cho bài toán chịu nhiệt một chiều. Trong thực tế làm việc phức tạp (như chịu tải trọng uốn, cắt và lực di động dọc trục), kết cấu cần tối ưu hóa độ bền và độ cứng theo cả phương dọc trục. Vật liệu 2D-FGM kết hợp cấu trúc sandwich giúp triệt tiêu hiện tượng tách lớp, cho phép phân bổ vật liệu hợp lý theo cả hai phương không gian nhằm đạt được tần số riêng mong muốn và giảm thiểu tập trung ứng suất.

2. Sự khác biệt bản chất giữa bài toán lực điều hòa di động (Chương 2) và khối lượng di động (Chương 3) là gì?

Trong bài toán lực điều hòa di động, tải trọng chỉ truyền lực kích động ngoại cảnh lên kết cấu mà không làm thay đổi đặc tính quán tính của hệ. Ngược lại, bài toán khối lượng di động xét đến khối lượng quán tính của vật thể di chuyển, làm sinh ra thêm lực quán tính thẳng đứng, lực quán tính Coriolis ($2v \dot{w}'$) và lực ly tâm ($v^2 w''$). Điều này dẫn đến sự xuất hiện của các ma trận khối lượng phụ thuộc thời gian $\mathbf{M}_m(t)$, ma trận cản Coriolis $\mathbf{C}_m(t)$ và ma trận độ cứng phụ thuộc thời gian $\mathbf{K}_m(t)$ trong phương trình động lực học tổng thể.

3. Kỹ thuật "phần tử làm giàu" (Enriched Element) trong Chương 3 mang lại ưu thế kỹ thuật gì?

Phần tử dầm làm giàu được xây dựng bằng cách bổ sung các hàm đa thức bậc cao vào trường nội suy chuyển vị chuẩn (tạo thành phần tử 2 nút có 22 bậc tự do). Kỹ thuật này giúp mô tả chính xác gradient biến dạng và trường ứng suất phức tạp trong dầm 2D-FGM mà không cần chia lưới quá mịn, cho phép đạt độ hội tụ chính xác cao chỉ với $NE = 6 \div 8$ phần tử thay vì phải dùng từ 30 đến 40 phần tử như phương pháp thông thường.

4. Hiện tượng nghẽn trượt (shear locking) trong tấm mỏng được luận án xử lý bằng cách nào trong Chương 4?

Khi sử dụng lý thuyết biến dạng trượt bậc nhất Mindlin cho các kết cấu tấm mỏng ($a/h$ lớn), các phần tử hữu hạn thông thường thường gặp hiện tượng khóa/nghẽn trượt do ma trận độ cứng trượt chiếm ưu thế không chính xác. Luận án giải quyết triệt để vấn đề này bằng cách áp dụng kỹ thuật nội suy liên kết (linked interpolation) trên phần tử tứ giác 9 nút (Q9), liên kết hàm nội suy của chuyển vị ngang $w$ với đạo hàm các góc xoay $\theta_x, \theta_y$, giúp phần tử đạt độ chính xác cao trên toàn bộ dải độ mảnh từ tấm dày ($a/h=5$) đến tấm mỏng ($a/h=100$) mà không cần sử dụng hệ số hiệu chỉnh trượt nhân tạo.

5. Sự khác biệt về kết quả dao động khi sử dụng các mô hình đồng nhất hóa vật liệu (Voigt, Mori-Tanaka, Maxwell) là gì?

Mô hình Voigt dựa trên quy tắc hỗn hợp tuyến tính đơn giản nên đánh giá độ cứng tổng thể của kết cấu cao hơn. Do đó, tần số dao động tự nhiên nhận được từ mô hình Voigt luôn lớn hơn so với tần số tính từ mô hình Mori-Tanaka (đối với hệ hai pha) và mô hình Maxwell (đối với hệ ba pha). Ngược lại, hệ số động lực học $D_d$ và độ võng động lực học tính theo mô hình Voigt có xu hướng nhỏ hơn so với các mô hình vi cơ học còn lại.

Kết luận

Luận án của tác giả Phạm Vũ Nam đã xây dựng thành công hệ thống mô hình phần tử hữu hạn nâng cao và chương trình tính toán số để phân tích toàn diện bài toá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 chịu tải trọng di động. Nghiên cứu đã phát triển các phần tử chuyên dụng gồm phần tử dầm làm giàu 22 bậc tự do theo lý thuyết trượt lượng giác và phần tử tấm tứ giác Q9 nội suy liên kết chống nghẽn trượt. Kết quả khảo sát số đã làm sáng tỏ quy luật ứng xử cơ học của kết cấu dưới ảnh hưởng đồng thời của quy luật phân bố vật liệu 2 chiều, điều kiện biên, hiệu ứng tải trọng di động và các mô hình đồng nhất hóa vi cơ học. Công trình là nguồn tài liệu tham khảo có giá trị học thuật và thực tiễn cao trong lĩnh vực tính toán cơ học kết cấu composite tiên tiến.