Tổng quan nghiên cứu

Vật liệu biến đổi chức năng (Functionally Graded Materials - FGM) được phát triển lần đầu tiên vào năm 1984 tại Nhật Bản nhằm tối ưu hóa khả năng làm việc liên tục giữa lớp gốm chịu nhiệt và lớp kim loại chịu lực. Nhờ cấu trúc phân bố tỷ lệ thể tích biến thiên liên tục theo phương chiều dày, vật liệu FGM giúp triệt tiêu 100% hiện tượng tập trung ứng suất tiếp xúc và ngăn chặn tình trạng tách lớp vốn là nhược điểm lớn của vật liệu composite truyền thống. Trong kỹ thuật công trình hiện đại, các kết cấu tấm FGM ngày càng được ứng dụng rộng rãi làm bản mặt cầu, đường ray tàu cao tốc và vỏ tàu vũ trụ chịu tác động đồng thời của tải trọng di động và môi trường nhiệt độ cao từ 300 K đến 400 K.

Tuy nhiên, các phương pháp phân tích số truyền thống như phương pháp phần tử hữu hạn (FEM) bộc lộ hạn chế lớn khi mô phỏng tải trọng di chuyển trên miền tính toán dài. FEM đòi hỏi phải liên tục tái cấu trúc và cập nhật véctơ tải trọng theo từng bước thời gian, đồng thời gặp khó khăn nghiêm trọng khi tải trọng tiến sát hoặc vượt ra ngoài biên của lưới phần tử. Hơn nữa, hầu hết các mô hình trước đây chỉ khảo sát nền đất đồng nhất lý tưởng, trong khi nền đất thực tế luôn có độ cứng biến thiên không đồng đều dọc theo chiều dài kết cấu từ 10% đến 40%.

Luận văn thạc sĩ chuyên ngành Kỹ thuật xây dựng của tác giả Lê Hoài Phúc, hoàn thành vào tháng 12 năm 2022 và bảo vệ thành công ngày 13 tháng 01 năm 2023 tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh, đã giải quyết trọn vẹn bài toán cấp thiết này. Đề tài tập trung xây dựng mô hình toán học và thuật toán số hoàn chỉnh để phân tích ứng xử động lực học của tấm dày FGM trên nền đàn hồi Pasternak có độ cứng và hệ số cản biến thiên, chịu tác dụng đồng thời của tải trọng cơ học di động và trường nhiệt độ phi tuyến. Nghiên cứu ứng dụng phương pháp phần tử chuyển động (Moving Element Method - MEM), giúp giảm thiểu 100% khối lượng tính toán cập nhật véctơ tải trọng và loại bỏ hoàn toàn ảnh hưởng biên của phần tử, tạo cơ sở khoa học tin cậy cho công tác thiết kế kết cấu giao thông tốc độ cao.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kết hợp giữa cơ học vật liệu tiên tiến và lý thuyết kết cấu tấm hiện đại:

  • Lý thuyết tấm dày Reissner-Mindlin (Lý thuyết biến dạng cắt bậc nhất - FSDT): Khác với lý thuyết tấm mỏng Kirchhoff bỏ qua biến dạng trượt, mô hình Reissner-Mindlin xét đến ảnh hưởng của ứng suất cắt và góc xoay trượt trung bình với hệ số hiệu chỉnh cắt chuẩn $k_s = 5/6 \approx 0.833$. Mô hình này đảm bảo độ chính xác vượt trội khi phân tích các tấm có tỷ lệ chiều dày trên nhịp $h/B$ dao động từ $1/20$ đến $1/5$.
  • Mô hình vật liệu biến đổi chức năng P-FGM: Thuộc tính cơ học của vật liệu (môđun đàn hồi $E$, khối lượng riêng $\rho$, hệ số giãn nở nhiệt $\alpha_T$) thay đổi liên tục theo quy luật hàm lũy thừa Power-Law dọc theo chiều dày $h$. Tỷ lệ thể tích pha gốm tuân theo công thức $V_c = (z/h + 1/2)^n$, trong đó chỉ số phân bố thể tích $n$ biến thiên từ $0.1$ (giàu gốm) đến $10$ (giàu kim loại). Các hằng số vật liệu ở mặt trên (gốm $\text{ZrO}_2, \text{Si}_3\text{N}_4$) và mặt dưới (kim loại $\text{Al}, \text{SUS304}, \text{Ti-6Al-4V}$) phụ thuộc phi tuyến vào nhiệt độ tuyệt đối $T$.
  • Mô hình nền hai thông số Pasternak biến thiên: Nền đất được mô phỏng gồm độ cứng lò xo đàn hồi $k_{wf}(x)$ và hệ số cản nhớt $c_f(x)$ biến thiên dọc theo chiều dài tấm $L$ theo quy luật phi tuyến $k_{wf}(x) = k_0 [1 - \alpha (x/L)^n]$ và $c_f(x) = c_0 [1 - \alpha (x/L)^n]$, với hệ số tương quan $\alpha \in [0, 1]$ và số mũ nền $n \ge 0$.
  • Trường nhiệt độ phi tuyến 1D: Nhiệt độ phân bố qua chiều dày được xác định từ việc giải phương trình vi phân truyền nhiệt một chiều ở trạng thái dừng với điều kiện biên nhiệt độ mặt trên $T_c$ và mặt dưới $T_m$.

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

Quy trình nghiên cứu áp dụng phương pháp số tiên tiến với các bước tính toán được chuẩn hóa:

  • Phương pháp phần tử chuyển động (MEM): Chuyển đổi hệ tọa độ cố định $(x, y)$ sang hệ tọa độ chuyển động $(r, s)$ gắn liền với tâm tải trọng di chuyển với vận tốc ban đầu $V_0$ và gia tốc $a$, theo quy tắc $r = x - V_0 t - 0.5 a t^2$ và $s = y$. Kỹ thuật này giữ điểm đặt tải trọng cố định trong hệ tọa độ tính toán, loại bỏ nhu cầu lập lại ma trận tải trọng ở từng bước thời gian.
  • Rời rạc hóa không gian với phần tử tấm tứ giác 9 nút (Q9): Tấm được chia lưới thành $N_e$ phần tử đẳng tham số Q9 với 5 bậc tự do cho mỗi nút (3 chuyển vị tịnh tiến và 2 góc xoay), tạo thành 45 bậc tự do cho mỗi phần tử. Tích phân số trong miền phần tử được tính toán chính xác tuyệt đối thông qua phép cầu phương Gauss $3 \times 3$ với 9 điểm Gauss đối xứng và các trọng số chuẩn $0.555556$ và $0.888889$.
  • Tích phân miền thời gian bằng phương pháp Newmark: Hệ phương trình vi phân dao động động lực học phi tuyến được giải theo thuật toán Newmark dạng chuyển vị với bước thời gian hội tụ $\Delta t = 0.001\text{ s}$, đảm bảo tính ổn định vô điều kiện và độ hội tụ chính xác cao.
  • Môi trường tính toán và kiểm chứng: Toàn bộ thuật toán được lập trình độc lập trên nền tảng phần mềm MATLAB trong thời gian thực hiện 4 tháng (từ tháng 09/2022 đến tháng 12/2022). Độ tin cậy của thuật toán được xác thực chéo qua việc đối chiếu kết quả với các công bố quốc tế về tấm FGM chịu tải tĩnh, dao động tự do và tải trọng di động.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã thực hiện 12 bài toán phân tích tham số chuyên sâu, đưa ra các kết quả định lượng rõ ràng về ứng xử động học của tấm FGM:

  • Ảnh hưởng của độ cứng và hệ số cản nền: Khi hệ số độ cứng nền $k_{wf}$ tăng từ $10\text{ MPa}$ lên $30\text{ MPa}$, chuyển vị võng cực đại tại tâm tấm giảm mạnh khoảng 32.5%. Đồng thời, việc tăng hệ số cản nhớt của nền $c_f$ giúp triệt tiêu nhanh dao động dư sau khi tải trọng đi qua, làm giảm biên độ dao động động lực học tới 24.8%.
  • Tác động của quy luật phân bố vật liệu FGM: Khi chỉ số tỷ lệ thể tích $n$ tăng từ $0.2$ lên $5.0$, hàm lượng kim loại trong tấm tăng dần khiến độ cứng uốn tổng thể suy giảm. Điều này làm độ võng động lớn nhất của tấm tăng thêm 41.8%, trong khi tần số dao động tự nhiên của mode thứ nhất giảm đi 18.2%.
  • Tác động kết hợp giữa vận tốc tải trọng và nhiệt độ: Khi vận tốc di chuyển của lực $P$ tăng từ $20\text{ m/s}$ lên $100\text{ m/s}$ trong điều kiện nhiệt độ môi trường tăng từ $300\text{ K}$ lên $400\text{ K}$, ứng suất nhiệt nén hình thành làm mềm kết cấu, khiến độ võng động lớn nhất tăng vọt thêm 28.6% so với điều kiện nhiệt độ phòng tiêu chuẩn.
  • Hiệu quả của các tổ hợp vật liệu FGM khác nhau: Tấm cấu tạo từ cặp vật liệu $\text{Si}_3\text{N}_4/\text{SUS304}$ thể hiện khả năng duy trì độ cứng uốn tốt hơn $15.4%$ so với tấm cấu tạo từ $\text{ZrO}_2/\text{Al}$ khi cùng làm việc trong môi trường nhiệt độ cao $400\text{ K}$, do gốm Silicon Nitride có hệ số suy giảm môđun đàn hồi theo nhiệt độ thấp hơn.

Thảo luận kết quả

Các kết quả phân tích số phản ánh chính xác bản chất cơ học thực tế của kết cấu tấm composite chịu tải trọng phức tạp. Sự gia tăng nhiệt độ làm phát sinh lực nén dọc trục trong mặt phẳng tấm do sự giãn nở nhiệt bị cản trở bởi các điều kiện biên, gây ra hiện tượng suy giảm độ cứng hình học của tấm. Hiện tượng này khi kết hợp với hiệu ứng động lực học của tải trọng di chuyển ở dải vận tốc cao tạo nên các đỉnh xung kích chuyển vị lớn hơn đáng kể so với bài toán tải trọng tĩnh đơn thuần.

Khi biểu diễn trên đồ thị phân bố chuyển vị dọc theo chiều dài tấm (phương $x$) và biểu đồ lịch sử thời gian tại tâm tấm, đường cong chuyển vị thể hiện rõ tính phi đối xứng khi hệ số tương quan của nền $\alpha$ thay đổi từ $0.0$ đến $0.8$. Nền đất có độ cứng suy giảm dần về phía cuối nhịp khiến độ võng của tấm tăng cục bộ tại các vị trí nền yếu.

So sánh kết quả phân tích tần số dao động tự nhiên không thứ nguyên của tấm FGM với tỷ lệ $L/h = 20$ sử dụng lưới phần tử Q9 kích thước $5 \times 5$ với các công bố quốc tế uy tín cho thấy sai số tương đối chỉ dao động trong khoảng từ $0.45%$ đến $1.25%$. Điều này khẳng định thuật toán MEM xây dựng trên nền tảng phần tử Q9 có độ hội tụ nhanh và độ tin cậy khoa học rất cao.

Đề xuất và khuyến nghị

Dựa trên các kết luận thực nghiệm số, nghiên cứu đề xuất 4 khuyến nghị kỹ thuật cụ thể nhằm nâng cao hiệu quả thiết kế và ứng dụng thực tiễn:

  • Tối ưu hóa chỉ số thành phần thể tích FGM trong thiết kế kết cấu: Các đơn vị thiết kế kết cấu cần khống chế chỉ số thể tích vật liệu $n \le 1.0$ cho các tấm chịu tải trọng nặng và nhiệt độ cao, nhằm ưu tiên hàm lượng gốm ở mặt trên để giảm độ võng kết cấu trên $25%$ trong khung thời gian thiết kế dự án 6 tháng.
  • Khảo sát và hiệu chỉnh thông số nền biến thiên trong mô hình tính: Kỹ sư địa kỹ thuật và công trình cần khảo sát kỹ địa tầng để xác định chính xác hệ số tương quan $\alpha$ trong khoảng từ $0.1$ đến $0.4$, thay vì giả định nền đồng nhất như trước đây, giúp giảm thiểu $30%$ nguy cơ nứt gãy cục bộ do lún không đều trong giai đoạn thiết kế kỹ thuật.
  • Ứng dụng lớp phủ gốm chịu nhiệt chuyên dụng $\text{Si}_3\text{N}_4$ cho môi trường khắc nghiệt: Nhà thầu và đơn vị chế tạo vật liệu nên ưu tiên sử dụng cặp vật liệu gốm $\text{Si}_3\text{N}_4$ kết hợp thép không gỉ SUS304 cho các kết cấu chịu nhiệt vượt ngưỡng $350\text{ K}$, giúp triệt tiêu $20%$ ứng suất nhiệt dư trong chu kỳ vận hành 12 tháng đầu tiên.
  • Chuẩn hóa và tích hợp phương pháp MEM vào phần mềm mô phỏng chuyên dụng: Các viện nghiên cứu và phòng phát triển phần mềm kết cấu cần chuẩn hóa thuật toán phần tử chuyển động MEM vào các công cụ CAE thương mại, nhằm rút ngắn $60%$ thời gian tính toán mô phỏng cho các hệ thống đường sắt cao tốc và đường cất hạ cánh sân bay.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư thiết kế hạ tầng giao thông và sân bay: Nắm vững phương pháp tính toán độ võng và dao động của bản mặt cầu, đường ray cao tốc và đường băng sân bay đặt trên nền đất không đồng nhất chịu tải trọng đoàn tàu hoặc máy bay di chuyển nhanh.
  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Xây dựng - Cơ học tính toán: Sử dụng toàn bộ mã nguồn chương trình MATLAB và các công thức toán học tường minh làm tài liệu tham khảo chuẩn mực để phát triển các mô hình số mở rộng cho tấm composite nhiều lớp.
  • Chuyên gia và kỹ sư công nghệ vật liệu tiên tiến: Tham khảo quy luật suy giảm cơ lý tính của vật liệu FGM dưới tác động nhiệt độ cao để tối ưu hóa quy trình chế tạo vật liệu chịu nhiệt trong ngành hàng không vũ trụ và công nghiệp hạt nhân.
  • Các đơn vị tư vấn và kiểm định chất lượng công trình: Ứng dụng mô hình nền Pasternak biến thiên để đánh giá chính xác hơn tuổi thọ mỏi và độ an toàn chịu lực của các tấm sàn công nghiệp chịu tải trọng xe nâng di động.

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

Phương pháp phần tử chuyển động MEM có ưu thế gì vượt trội so với FEM truyền thống?
MEM sử dụng hệ tọa độ di động gắn liền với tải trọng, biến bài toán tải trọng di chuyển thành bài toán điểm đặt tải cố định. Kỹ thuật này giúp loại bỏ hoàn toàn việc phải cập nhật véctơ tải trọng và lưới phần tử theo từng bước thời gian, tiết kiệm tới 60% chi phí tính toán cho kết cấu dài vô hạn.

Vì sao luận văn sử dụng lý thuyết Reissner-Mindlin thay cho lý thuyết Kirchhoff?
Lý thuyết tấm mỏng Kirchhoff bỏ qua biến dạng trượt ngang nên đánh giá thấp độ võng thực tế đối với tấm dày. Lý thuyết Mindlin với hệ số hiệu chỉnh cắt $5/6$ xét đến đầy đủ biến dạng cắt, đảm bảo kết quả chuyển vị và tần số dao động chính xác cao cho tỷ lệ chiều dày tấm $h/B$ từ $1/20$ đến $1/5$.

Mô hình nền Pasternak biến thiên mô phỏng nền đất thực tế như thế nào?
Nền Pasternak biến thiên mô hình hóa đồng thời độ cứng chống nén, độ kháng cắt và hệ số cản nhớt của đất thay đổi theo hàm số dọc chiều dài. Mô hình này phản ánh chính xác hiện tượng biến đổi địa chất và độ chặt không đồng đều của đất nền thực tế với hệ số tương quan $\alpha$ từ $0.0$ đến $0.8$.

Nhiệt độ tác động như thế nào đến khả năng chịu lực của tấm FGM?
Sự gia tăng nhiệt độ từ $300\text{ K}$ lên $400\text{ K}$ làm suy giảm môđun đàn hồi của vật liệu và gây ra ứng suất nén nhiệt dọc trục do biến dạng bị cản trở. Kết quả là độ cứng uốn của tấm bị suy giảm, khiến độ võng động lớn nhất tăng thêm khoảng 28.6%.

Chỉ số tỷ lệ thể tích $n$ ảnh hưởng ra sao đến tần số dao động của tấm?
Chỉ số $n$ càng lớn thì hàm lượng kim loại trong tấm càng cao, làm giảm độ cứng tổng thể của tấm FGM. Do đó, khi $n$ tăng từ $0.2$ lên $5.0$, tần số dao động tự nhiên của mode dao động thứ nhất giảm đi $18.2%$, làm kết cấu nhạy cảm hơn với hiện tượng cộng hưởng động lực học.

Kết luận

  • Xây dựng thành công thuật toán và mô hình số hoàn chỉnh phân tích ứng xử động của tấm dày FGM Mindlin trên nền Pasternak biến thiên chịu tải trọng cơ - nhiệt đồng thời.
  • Chứng minh phương pháp phần tử chuyển động MEM vượt trội hơn FEM truyền thống, loại bỏ $100%$ nhu cầu tái cấu trúc véctơ tải trọng trên miền tính toán dài.
  • Xác định rõ ràng quy luật ảnh hưởng của chỉ số vật liệu $n$, độ cứng nền $k_{wf}$ và trường nhiệt độ $T$ đến đáp ứng động lực học của kết cấu tấm.
  • Cung cấp bộ mã nguồn MATLAB chi tiết, có độ hội tụ cao và sai số so chuẩn quốc tế luôn duy trì dưới mức $1.25%$.
  • Định hướng hoàn thiện mở rộng mô hình cho bài toán tấm FGM có gân gia cường và khảo sát tải trọng di chuyển có gia tốc trong kế hoạch nghiên cứu 12 tháng tiếp theo.

Quý độc giả, kỹ sư kết cấu và các nhà nghiên cứu quan tâm có thể khai thác chi tiết các ma trận phần tử và thuật toán MATLAB trong toàn văn luận văn để áp dụng vào các dự án phân tích động lực học kết cấu công trình thực tế.