Phân Tích Ứng Xử Của Dầm Composite Thành Mỏng Dùng Lý Thuyết Biến Dạng Cắt Bậc Nhất (T2020-77TĐ)

Tóm tắt nghiên cứu (Research Summary)

Nghiên cứu "Phân tích ứng xử của dầm composite thành mỏng dùng lý thuyết biến dạng cắt bậc nhất" (Mã số đề tài: T2020-77TĐ), do TS. Nguyễn Ngọc Dương chủ nhiệm cùng nhóm nghiên cứu tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh thực hiện, tập trung giải quyết bài toán phân tích ứng xử uốn tĩnh của dầm composite chức năng (Functionally Graded Material - FGM) có tiết diện thành mỏng hở (tiết diện chữ I) đặt trên nền đàn hồi hai thông số.

  • Câu hỏi nghiên cứu cốt lõi: Làm thế nào để mô hình hóa chính xác đồng thời ảnh hưởng của hiệu ứng biến dạng cắt (shear deformation) và tương tác nền đàn hồi lên ứng xử uốn của dầm FGM thành mỏng mà vẫn tối ưu hóa khối lượng tính toán?
  • Phương pháp tiếp cận: Đề xuất lý thuyết biến dạng cắt bậc nhất cải tiến (Improved First-Order Shear Deformation Theory) dựa trên kỹ thuật tách trường chuyển vị thành hai thành phần uốn ($b$) và cắt ($s$), kết hợp nguyên lý cực trị năng lượng Lagrange và phương pháp Ritz với các hàm dạng đa thức mũ thỏa mãn 4 điều kiện biên cơ bản (S-S, C-F, C-S, C-C).
  • Kết quả then chốt: Mô hình hội tụ nhanh ở số hàm dạng $m = 10$, phản ánh chính xác ảnh hưởng của tỷ số nhịp/chiều cao ($L/b_3$), chỉ số gradient vật liệu ($p$), phân bố bề dày gốm ($\eta$) và các tham số nền ($k_z, g_z$). Hiệu ứng biến dạng cắt thể hiện rõ nhất ở dầm ngàm hai đầu (C-C) và giảm dần khi độ mảnh tăng; thông số nền thứ hai ($g_z$) làm giảm độ võng của dầm từ 25% đến 38%.
  • Ý nghĩa thực tiễn: Cung cấp công cụ giải tích số tin cậy, tiết kiệm chi phí tính toán cho việc thiết kế các cấu kiện kết cấu composite siêu nhẹ chịu lực trong kỹ thuật hàng không, cơ khí chính xác và công trình biển.

Bối cảnh và tầm quan trọng (Context & Significance)

1. Hiện trạng tri thức và rào cản nghiên cứu

Vật liệu composite chức năng (FGM) ra đời nhằm khắc phục triệt để nhược điểm bong tách lớp (delamination) vốn rất phổ biến ở vật liệu composite phân lớp truyền thống (LCM). Nhờ sự biến đổi trơn, liên tục của thành phần cơ tính (từ kim loại sang gốm) theo các hướng xác định, FGM đem lại khả năng chịu nhiệt, chống mài mòn và phân bổ ứng suất tối ưu.

Trong tính toán kết cấu dầm thành mỏng, lý thuyết cổ điển Vlasov đóng vai trò nền tảng nhưng lại bỏ qua hoàn toàn biến dạng cắt, chỉ phù hợp với các dầm rất mảnh. Khi dầm làm việc ở dải kích thước trung bình hoặc chịu tải trọng lớn, biến dạng cắt cục bộ và hiện tượng vênh mặt cắt (warping) trở nên đáng kể, dẫn đến sai số lớn nếu áp dụng lý thuyết cổ điển.

2. Khoảng trống học thuật và tính cấp thiết

Mặc dù đã có nhiều công trình nghiên cứu về dầm FGM tiết diện chữ nhật hoặc dầm phân lớp LCM, các mô hình xem xét đồng thời hiệu ứng biến dạng cắt, tiết diện thành mỏng phức tạp (dạng chữ I) và tương tác nền đàn hồi hai thông số (Winkler-Pasternak) vẫn còn rất hạn chế. Hơn nữa, việc phát triển phương pháp giải tích - xấp xỉ Ritz cho dầm thành mỏng FGM chưa được phổ biến rộng rãi so với phương pháp phần tử hữu hạn (FEM) vốn đòi hỏi tài nguyên tính toán lớn.

3. Ý nghĩa và tiềm năng ứng dụng

Đề tài có tính thời sự cao, lấp đầy khoảng trống lý thuyết về tính toán dầm FGM thành mỏng tại Việt Nam. Kết quả nghiên cứu là nền tảng để tối ưu hóa thiết kế kết cấu dạng thanh mỏng trong các ngành công nghiệp đòi hỏi hiệu suất chịu lực/trọng lượng cao như hàng không vũ trụ, vỏ tàu thủy, giàn khoan ngoài khơi và cầu khẩu độ lớn.

Phương pháp luận và cách tiếp cận (Methodology & Approach)

Nghiên cứu xây dựng một quy trình giải tích số nhất quán, chặt chẽ từ mô hình động học đến thuật toán giải phương trình đại số:

1. Mô hình trường động học và chuyển vị tách biến

Hệ tọa độ gồm: hệ tọa độ Descartes tổng thể $(x, y, z)$, hệ tọa độ địa phương trên thành mỏng $(n, s, z)$ và hàm quạt $\omega(s)$ mô tả hiện tượng vênh do xoắn. Điểm đột phá là việc biểu diễn các chuyển vị $V(x)$, $W(x)$ và góc xoay $\Phi(x)$ dưới dạng tổng của hai thành phần uốn ($b$) và cắt ($s$):

$$\begin{aligned} V(x) &= V_b(x) + V_s(x) \ W(x) &= W_b(x) + W_s(x) \ \Phi(x) &= \Phi_b(x) + \Phi_s(x) \end{aligned}$$

Từ đó, trường chuyển vị dọc trục $u(s, x)$ được thiết lập thông qua đạo hàm uốn và hàm quạt, giảm thiểu độ phức tạp so với các lý thuyết bậc cao thông thường mà vẫn bảo toàn đầy đủ các hiệu ứng cắt $\gamma_{nx}$ và $\gamma_{sx}$.

2. Quan hệ ứng suất - biến dạng và phân bố vật liệu

Đề tài xét hai cấu hình phân bố vật liệu FGM trên tiết diện chữ I:

  • Tiết diện S1: Cánh trên và cánh dưới tuân theo hàm mũ phân lớp (Loại A), bản bụng dầm tuân theo hàm lũy thừa liên tục (Loại B).
  • Tiết diện S2: Toàn bộ cánh và bụng đều tuân theo phân bố Loại B, đảm bảo tính liên tục của mô đun đàn hồi tại các vị trí tiếp giáp giữa cánh và bụng dầm.

Mô đun đàn hồi hiệu dụng được chuẩn hóa:

$$E(n) = E_c V_c(n) + E_m [1 - V_c(n)]$$

3. Năng lượng hệ thống và thuật toán Ritz

Thế năng toàn phần của hệ bao gồm năng lượng biến dạng của dầm ($\Pi_E$), năng lượng biến dạng của nền đàn hồi hai thông số Pasternak ($\Pi_F$) và công của ngoại lực ($\Pi_W$):

$$\Pi = \Pi_E + \Pi_F - \Pi_W$$

Trong đó, tương tác nền đàn hồi bao gồm độ cứng lò xo Winkler ($k_y, k_z$), lớp cắt Pasternak ($g_y, g_z$) và đặc trưng kháng xoay ($k_t$). Áp dụng nguyên lý Lagrange $\frac{\partial \Pi}{\partial p_j} = 0$, hệ phương trình cân bằng tĩnh được đưa về dạng phương trình ma trận độ cứng:

$$[K] {p} = {F}$$

Với ${p} = {U_j, V_{bj}, W_{bj}, \Phi_{bj}, V_{sj}, W_{sj}, \Phi_{sj}}^T$, các hàm dạng Ritz được chọn là dạng đa thức kết hợp hàm mũ thỏa mãn hoàn toàn 4 điều kiện biên: Tựa đơn - tựa đơn (S-S), Ngàm - tự do (C-F), Ngàm - tựa đơn (C-S), Ngàm - ngàm (C-C).

Phát hiện chính của nghiên cứu (Key Findings)

Dựa trên các ví dụ số phân tích dầm chữ I chịu tải tập trung $P_z$ và tải phân bố đều $q_z$, nghiên cứu rút ra các kết quả quan trọng sau:

1. Tốc độ hội tụ và độ chính xác của mô hình

Nghiên cứu khảo sát sự hội tụ với các bậc hàm dạng Ritz từ $m = 2$ đến $m = 12$. Kết quả cho thấy lời giải hội tụ xuất sắc và ổn định tuyệt đối khi $m \ge 10$. So sánh với các công bố uy tín quốc tế của Kim & Lee (2015, 2017) cho dầm đồng nhất và dầm FGM:

  • Chuyển vị dầm đơn giản đồng nhất không xét biến dạng cắt (KXC) đạt 39.9465 mm (sai số $< 0.1%$ so với nghiệm tham chiếu).
  • Chuyển vị có xét biến dạng cắt (XC) phản ánh chính xác độ mềm tăng thêm của kết cấu.

2. Ảnh hưởng của tỷ số độ mảnh ($L/b_3$) và điều kiện biên

Hiệu ứng biến dạng cắt ($W_s / W \times 100%$) phụ thuộc chặt chẽ vào điều kiện liên kết biên và hình học dầm:

  • Ảnh hưởng cắt theo điều kiện biên: $\text{C-C} > \text{C-S} > \text{S-S} > \text{C-F}$. Dầm ngàm hai đầu (C-C) chịu ảnh hưởng của biến dạng cắt lớn nhất do độ cứng uốn bị kiềm chế tối đa, làm thành phần biến dạng cắt chiếm tỷ trọng cao trong tổng chuyển vị.
  • Quy luật suy giảm: Khi tỷ số $L/b_3$ tăng từ 10 lên 70 (dầm trở nên mảnh hơn), ảnh hưởng của biến dạng cắt giảm mạnh từ mức trên $15%$ xuống dưới $2%$.

3. Tác động của chỉ số gradient vật liệu ($p$) và phân bố bề dày gốm

  • Chỉ số thể tích $p$: Khi $p$ tăng từ 0 (thuần gốm) lên các giá trị cao hơn, hàm lượng kim loại tăng lên làm giảm độ cứng tổng thể của dầm, khiến chuyển vị giữa nhịp tăng đơn điệu ở tất cả các điều kiện biên.
  • Tỷ số độ cứng uốn/cắt ($E_{33}/E_{77}$): Biến dạng cắt tỷ lệ thuận với tỷ số $E_{33}/E_{77}$. Khi hàm lượng gốm ở bụng dầm tăng, hiệu ứng cắt giảm đáng kể; trong khi đó, tăng chiều dày gốm ở cánh dưới làm tăng ảnh hưởng của biến dạng cắt.

4. Tác động của các thông số nền đàn hồi ($k_z, g_z$)

  • Hệ số Winkler ($k_z$): Ít ảnh hưởng đến dầm ngắn ($L/b_3 \le 20$), nhưng trở nên rõ nét đối với dầm mảnh ($L/b_3 \ge 70$).
  • Hệ số Pasternak ($g_z$): Lớp cắt đàn hồi thứ hai có vai trò quyết định, làm giảm chuyển vị dầm C-C khoảng 25% và dầm C-S khoảng 38%. Dầm công xôn (C-F) là kết cấu nhạy cảm và hưởng lợi nhiều nhất từ độ cứng của nền đàn hồi.

Đóng góp khoa học và giá trị học thuật (Scientific Contributions)

Lĩnh vực Nội dung đóng góp cụ thể
Đóng góp lý thuyết Phát triển thành công lý thuyết dầm biến dạng cắt bậc nhất cải tiến cho dầm thành mỏng FGM; tích hợp linh hoạt cả hai trường hợp có xét (XC) và không xét (KXC) biến dạng cắt trong cùng một cấu trúc phương trình giải tích.
Cải tiến phương pháp Ứng dụng phương pháp Ritz kết hợp hàm dạng mũ cho bài toán dầm thành mỏng FGM nhiều ẩn; giúp giảm đáng kể thời gian lập trình và tài nguyên tính toán so với các phần tử 3D FEM phức tạp.
Giá trị thực tiễn Cung cấp dữ liệu chuẩn (benchmark data) và biểu đồ tham chiếu độ cứng cho các kỹ sư kết cấu khi lựa chọn quy cách dầm chữ I làm từ vật liệu composite biến tính.
Công bố quốc tế Đề tài đã hoàn thành sản phẩm phê duyệt với 01 bài báo khoa học đăng trên tạp chí quốc tế chuyên ngành thuộc danh mục Scimago Q1.

Đối tượng quan tâm và lợi ích ứng dụng (Target Audience)

  • Nhà nghiên cứu Cơ học tính toán & Kết cấu: Kế thừa khung lý thuyết tách biến để mở rộng sang các bài toán phân tích phi tuyến hình học, động lực học rung động tự do, ổn định mất ổn định nhiệt - ẩm (hygro-thermal buckling).
  • Kỹ sư Thiết kế Kết cấu Công trình & Cầu hầm: Sử dụng công thức tính nhanh độ võng dầm FGM đặt trên nền đất mềm để tối ưu hóa móng cọc, dầm giằng móng tiết diện mỏng.
  • Chuyên gia Ngành Hàng không, Tàu thủy & Xe hơi cao cấp: Ứng dụng mô hình dầm chữ I composite siêu bền - nhẹ nhằm cắt giảm tự trọng phương tiện, tăng khả năng tiết kiệm nhiên liệu và chịu tải va đập.

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

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Phát hiện then chốt là việc xác định định lượng ảnh hưởng của biến dạng cắt đối với dầm FGM thành mỏng (đạt cực đại ở dầm liên kết ngàm C-C) và chứng minh rằng thông số nền thứ hai ($g_z$ - lớp cắt đàn hồi Pasternak) đóng vai trò chi phối trong việc giảm độ võng của kết cấu (giảm tới 38%), trong khi thông số lò xo đơn thuần $k_z$ ít có tác dụng với dầm ngắn.

2. Điểm đặc biệt của phương pháp tách trường chuyển vị là gì?

Khác với FSDT truyền thống sử dụng các ẩn góc xoay độc lập dẫn đến hệ phương trình cồng kềnh, phương pháp này tách chuyển vị thành hai thành phần uốn ($b$) và cắt ($s$). Kỹ thuật này vừa giữ được hình thức toán học tường minh, vừa cho phép chuyển đổi tức thì giữa bài toán có xét và không xét biến dạng cắt mà không cần xây dựng lại toàn bộ ma trận độ cứng.

3. Kết quả của nghiên cứu có thể tổng quát hóa cho các dạng tiết diện khác không?

Có. Mặc dù bài báo minh họa cụ thể cho tiết diện chữ I đối xứng và không đối xứng, hệ thống công thức biến dạng và ma trận độ cứng thiết lập trong phụ lục hoàn toàn có thể mở rộng cho các tiết diện thành mỏng hở khác như chữ C, chữ Z, chữ U hoặc tiết diện hộp kín.

4. Đâu là các hướng phát triển tiếp theo của đề tài?

Các hướng nghiên cứu tiếp nối bao gồm:

  1. Mở rộng bài toán sang phân tích động lực học (dao động riêng, tải trọng di động).
  2. Phân tích ổn định uốn - xoắn phi tuyến dưới tác động đồng thời của nhiệt độ và độ ẩm.
  3. Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) để tối ưu hóa hàm phân bố vật liệu $p$ theo yêu cầu chịu lực đa mục tiêu.

5. Khả năng ứng dụng thực tế vào sản xuất công nghiệp như thế nào?

Kết quả giúp các kỹ sư lựa chọn chính xác tỷ lệ phân bố giữa pha gốm ($E_c$) và pha kim loại ($E_m$) trong quy trình luyện kim bột hoặc in 3D vật liệu FGM, đảm bảo cấu kiện dầm thành mỏng vừa có khả năng chịu uốn tốt, vừa không bị phá hủy giòn do biến dạng cắt cục bộ.

Kết luận (Conclusion)

Báo cáo nghiên cứu khoa học mã số T2020-77TĐ của TS. Nguyễn Ngọc Dương đã xây dựng thành công một lý thuyết biến dạng cắt bậc nhất cải tiến có tính ứng dụng cao cho dầm FGM thành mỏng tiết diện chữ I trên nền đàn hồi. Bằng việc kết hợp hài hòa giữa kỹ thuật tách biến chuyển vị uốn-cắt và thuật toán xấp xỉ Ritz, nghiên cứu mang lại lời giải chính xác, tin cậy và tiết kiệm chi phí tính toán.

Công trình không chỉ đóng góp một bài báo chất lượng cao thuộc nhóm Scimago Q1 cho nền khoa học nước nhà, mà còn mở ra hướng đi đầy triển vọng trong việc thiết kế, chế tạo các kết cấu composite chức năng tiên tiến cho các ngành công nghiệp mũi nhọn.