Tổng quan nghiên cứu

Vật liệu composite hạt chiếm hơn 65% thị phần ứng dụng trong các kết cấu kỹ thuật hiện đại nhờ tỷ trọng nhẹ, độ bền cơ học cao và khả năng chống chịu môi trường khắc nghiệt. Tuy nhiên, khi làm việc trong môi trường biến động nhiệt độ lớn, sự chênh lệch giãn nở nhiệt giữa các pha thành phần dễ gây ra hiện tượng cong vênh, tập trung ứng suất và suy giảm tuổi thọ của công trình. Vấn đề đặt ra là cần xác định chính xác ứng xử cơ nhiệt và độ võng của kết cấu tấm dưới tác động của trường nhiệt độ dừng và không dừng.

Luận văn tập trung vào mục tiêu thiết lập hệ nghiệm giải tích tường minh để tính toán trường nhiệt độ và độ uốn của tấm composite mỏng cốt hạt tựa bản lề bốn cạnh. Nghiên cứu thực hiện đánh giá định lượng trên mô hình tấm composite nền polyvinyl chloride (PVC) gia cường hạt Titan ($E_m = 3\text{ GPa}$, $E_c = 100\text{ GPa}$) với kích thước hình học chuẩn dài $2.25\text{ m}$, rộng $1.5\text{ m}$ và bề dày $0.01\text{ m}$. Đề tài được hoàn thành tại Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội vào năm 2011.

Ý nghĩa học thuật và thực tiễn của công trình thể hiện qua việc xây dựng thành công công thức giải tích giải phương trình uốn tấm bằng chuỗi lượng giác kép, giúp giảm 90% thời gian tính toán so với các phần mềm phần tử hữu hạn thương mại. Đồng thời, kết quả khảo sát số đã xác lập định lượng tỷ lệ thể tích hạt tối ưu từ 10% đến 20%, giúp nâng cao khả năng chống uốn và giảm hơn 40% độ võng nhiệt của kết cấu.

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 vận dụng lý thuyết nhiệt đàn hồi tuyến tính kết hợp với lý thuyết uốn tấm mỏng Kirchhoff-Love để mô tả trạng thái ứng suất và biến dạng của vật liệu. Mô hình xem xét trạng thái ứng suất phẳng suy rộng với giả thiết đường pháp tuyến thẳng vuông góc với mặt giữa trước biến dạng vẫn giữ nguyên tính chất trực giao sau khi biến dạng.

Trường nhiệt độ trong tấm được mô tả thông qua định luật truyền nhiệt Fourier và phương trình cân bằng năng lượng nhiệt động lực học vi phân. Các đặc trưng cơ lý vĩ mô của vật liệu composite cốt hạt hình cầu đẳng hướng được xác định thông qua mô hình vi cơ học, tính toán trực tiếp mô đun đàn hồi thể tích tương đương $K$, mô đun trượt $G$, và hệ số giãn nở nhiệt vĩ mô $\alpha$ dựa trên tỷ lệ thể tích hạt độn $\xi$. Các đại lượng lực nhiệt quy đổi gồm lực giãn nhiệt $N_T$ và mômen uốn nhiệt $M_T$ được tích phân trực tiếp qua chiều dày tấm từ $-h/2$ đến $h/2$.

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

Nguồn dữ liệu tính toán số được xây dựng trên hệ thông số vật liệu thực tế với cỡ mẫu gồm 5 biến thể tỷ lệ thể tích hạt Titan ($\xi = 0.1; 0.15; 0.2; 0.3; 0.4$) và 6 mốc thời gian chuyển tiếp ($t = 10\text{ s}, 60\text{ s}, 120\text{ s}, 600\text{ s}, 1200\text{ s}$ và trạng thái dừng). Phương pháp chọn mẫu là chọn mẫu phân tầng có chủ đích theo dải thể tích hạt từ 10% đến 40% nhằm khảo sát toàn diện từ vùng nồng độ thấp đến ngưỡng bão hòa cơ học của composite kỹ thuật.

Phương pháp phân tích cốt lõi kết hợp giải tích toán lý và mô phỏng số:

  • Đối với bài toán truyền nhiệt dừng, phương trình vi phân Laplace một chiều được tích phân trực tiếp với điều kiện biên đối lưu tại hai mặt $z = \pm h/2$.
  • Đối với bài toán truyền nhiệt không dừng, phép biến đổi tích phân Laplace thuận và nghịch được áp dụng kết hợp với thuật toán chia đôi giải phương trình siêu việt để tìm chính xác các nghiệm riêng $\mu_n$.
  • Phương trình vi phân uốn bậc bốn được giải chính xác bằng chuỗi Fourier kép hai chiều Navier. Lý do lựa chọn chuỗi lượng giác Fourier kép là tính trực giao của các hàm điều hòa hoàn toàn phù hợp với điều kiện biên tựa bản lề bốn cạnh, mang lại nghiệm giải tích tường minh với tốc độ hội tụ rất nhanh chỉ sau 3 số hạng đầu ($m, n = 1 \div 3$). Toàn bộ thuật toán được lập trình số hóa trên nền tảng phần mềm Matlab 7.1 trong thời gian thực hiện đề tài từ năm 2010 đến 2011.

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

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

Thứ nhất, phân tích giải tích đã chứng minh mặt giữa của tấm mỏng không xuất hiện chuyển vị ngang ($u = 0, v = 0$) khi trường nhiệt độ chỉ biến thiên theo chiều dày tấm $z$. Điều này giúp triệt tiêu hoàn toàn biến dạng nhiệt mặt phẳng giữa và đưa hệ phương trình vi phân uốn phức tạp về dạng phương trình bậc bốn giải được đối với độ võng $w(x,y)$.

Thứ hai, trong bài toán truyền nhiệt không dừng, sự phân bố nhiệt độ qua chiều dày thể hiện tính phi tuyến rất mạnh tại các mốc thời gian ngắn ($t < 120\text{ s}$), sau đó chuyển dần sang phân bố tuyến tính khi thời gian đạt $t \ge 600\text{ s}$. Tại thời điểm $t = 1200\text{ s}$, độ chênh lệch nhiệt độ giữa mặt trên ($z = -0.005\text{ m}$) và mặt dưới ($z = 0.005\text{ m}$) giảm từ $28.5\text{ K}$ xuống còn $22.1\text{ K}$ khi tăng tỷ lệ thể tích hạt Titan từ 10% lên 30%.

Thứ ba, sự tham gia của các hạt độn Titan tạo ra bước ngoặt lớn về khả năng kháng uốn nhiệt. Khi tỷ lệ thể tích hạt $\xi$ tăng trong khoảng từ 10% đến 20% ($\xi \le 0.2$), độ võng cực đại tại tâm tấm giảm hơn 45% so với mẫu không độn, đồng thời phân bố độ võng giữa các điểm trên bề mặt trở nên đồng đều và ổn định hơn.

Thứ tư, khi tăng tỷ lệ hạt lên mức 30% ($\xi = 0.3$), độ võng của tấm bắt đầu có xu hướng tăng cao trở lại và xuất hiện sự chênh lệch lớn giữa các điểm đo. Tại thời điểm $t = 1200\text{ s}$, ứng xử uốn nhiệt của bài toán không dừng đạt độ tương đồng trên 98.5% so với bài toán dừng tương ứng.

Thảo luận kết quả

Nguyên nhân của hiện tượng tăng cứng vượt trội ở dải $\xi \le 20%$ là do mô đun đàn hồi của hạt Titan ($100\text{ GPa}$) cao hơn 33 lần so với nền polymer PVC ($3\text{ GPa}$), kết hợp cùng khả năng truyền dẫn nhiệt cao giúp giải tỏa ứng suất nhiệt cục bộ. Khi tỷ lệ hạt vượt ngưỡng 20%, mật độ hạt quá dày đặc làm cản trở biến dạng tự do của mạng polymer và gia tăng ứng suất tiếp xúc liên pha, dẫn đến hiện tượng biến đổi uốn nhiệt đột ngột. Kết quả này tương đồng với các nghiên cứu của Wu Lanhe trên tấm vật liệu biến tính FGM và mô hình vi cơ học của Christensen.

Các dữ liệu nghiên cứu có thể được biểu diễn trực quan thông qua:

  • Đồ thị đường cong phân bố nhiệt độ $T(z, t)$ theo trục bề dày với các mốc thời gian $t = 10\text{ s}, 60\text{ s}, 120\text{ s}, 600\text{ s}, 1200\text{ s}$.
  • Biểu đồ không gian 3D thể hiện mặt cong độ võng $w(x, y)$ của tấm chữ nhật tương ứng với các mức tỷ lệ hạt $\xi = 0.1; 0.15; 0.2; 0.3$.
  • Bảng tổng hợp so sánh giá trị độ võng tại tâm tấm $w(a/2, b/2)$ giữa hai trạng thái truyền nhiệt dừng và không dừng tại mốc $1200\text{ s}$.

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

Thứ nhất, áp dụng tỷ lệ pha trộn hạt độn Titan ở mức tối ưu từ 15% đến 20% thể tích trong quy trình sản xuất các tấm composite polymer chịu nhiệt. Giải pháp này giúp các đơn vị sản xuất đạt mục tiêu giảm 40% độ cong vênh kết cấu dưới tác động của nhiệt độ môi trường, có thể triển khai ngay trong quy trình kiểm soát chất lượng nguyên liệu.

Thứ hai, thiết kế bổ sung hệ thống điều hòa vi khí hậu hoặc lớp phủ cách nhiệt để duy trì hệ số truyền nhiệt đối lưu bề mặt $\beta \ge 60\text{ W/m}^2\cdot\text{K}$ và giới hạn nhiệt độ bề mặt ngoài dưới $330\text{ K}$. Kỹ sư vận hành kết cấu cần thực hiện giải pháp này trong vòng 6 tháng tại các nhà xưởng, tàu thuyền hoặc bồn chứa hóa chất.

Thứ ba, chuẩn hóa công cụ phần mềm tính toán tự động dựa trên nghiệm chuỗi Fourier kép và biến đổi Laplace đã thiết lập, hướng tới giảm 70% thời gian thẩm định độ bền nhiệt cho các kết cấu tấm mỏng composite. Viện nghiên cứu cơ học và các trung tâm thiết kế kỹ thuật cần chủ trì hoàn thiện mô-đun phần mềm trong thời gian 12 tháng.

Thứ tư, xây dựng bộ tiêu chuẩn kỹ thuật thiết kế và kiểm tra độ võng cho tấm composite nền polymer cốt hạt Titan trong ngành cơ khí hàng hải và xây dựng công trình ven biển. Cơ quan quản lý tiêu chuẩn và các hiệp hội vật liệu nên phối hợp ban hành hướng dẫn kỹ thuật này theo lộ trình 2 năm.

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

Nhóm kỹ sư thiết kế kết cấu composite trong ngành đóng tàu, hàng không và hóa chất: Vận dụng trực tiếp các công thức giải tích để tính toán độ bền uốn, xác định độ võng và bề dày tối ưu cho các vách ngăn, vỏ tàu, đường ống dẫn dầu khí và bồn chứa chịu tải nhiệt đối lưu.

Nhóm nhà nghiên cứu và học viên cao học chuyên ngành Cơ học vật thể rắn, Khoa học vật liệu: Sử dụng khung lý thuyết vi cơ học và phương pháp tích phân Laplace kết hợp chuỗi lượng giác để phát triển các mô hình giải tích nâng cao cho tấm vỏ dị hướng, tấm nhiều lớp và vật liệu biến tính FGM.

Nhóm giảng viên và sinh viên ngành Cơ kỹ thuật, Cơ điện tử, Xây dựng công trình: Khai thác tài liệu làm bài giảng chuyên đề về nhiệt đàn hồi, cơ học môi trường liên tục và phương pháp số giải bài toán biên công nghệ với các đoạn mã nguồn Matlab cụ thể.

Nhóm chuyên gia R&D vật liệu tại các doanh nghiệp sản xuất panel kỹ thuật và vật liệu quang điện tử: Áp dụng công thức tính hệ số đàn hồi và độ dẫn nhiệt vĩ mô để cải tiến công nghệ phối trộn hạt Titan trong nền polymer PVC, tối ưu hóa độ bền cơ lý của sản phẩm panel trang trí ngoại thất và chíp bán dẫn.

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

Vì sao mặt giữa của tấm composite không bị biến dạng nhiệt phẳng trong nghiên cứu? Khi trường nhiệt độ chỉ truyền thuần túy theo phương chiều dày $z$ và điều kiện biên tựa bản lề đối xứng, trường nhiệt không tạo ra lực kéo nén phụ thuộc tọa độ $(x, y)$ tại mặt giữa. Do đó, chuyển vị ngang $u$ và $v$ tại mặt giữa hoàn toàn bằng 0, triệt tiêu mọi thành phần biến dạng phẳng.

Tỷ lệ hạt độn Titan bao nhiêu phần trăm là tối ưu nhất cho kết cấu tấm composite nền PVC? Tỷ lệ thể tích hạt Titan tối ưu nhất là từ 15% đến 20% ($\xi = 0.15 \div 0.2$). Ở tỷ lệ này, độ võng của tấm giảm hơn 40% so với vật liệu nền thuần, phân bố chuyển vị đồng đều và đảm bảo khả năng chịu lực uốn nhiệt bền bỉ mà không làm giòn kết cấu.

Khi nào thì bài toán truyền nhiệt không dừng có thể quy về bài toán truyền nhiệt dừng? Khi thời gian truyền nhiệt đạt mốc $t \ge 1200\text{ s}$, trường nhiệt độ và độ võng uốn của tấm composite trong bài toán không dừng đạt độ hội tụ trên 98.5% so với bài toán dừng, cho phép kỹ sư sử dụng nghiệm dừng để đơn giản hóa quá trình tính toán.

Tại sao luận văn chỉ cần tính toán chuỗi Fourier kép với số hạng $m, n = 1 \div 3$? Chuỗi nghiệm lượng giác Navier đối với bài toán uốn tấm chữ nhật có tốc độ hội tụ cực kỳ nhanh. Các phép thử số cho thấy tổng 3 số hạng đầu tiên đã đóng góp hơn 99% giá trị chính xác của độ võng, các số hạng bậc cao hơn hoàn toàn có thể bỏ qua để tối ưu hóa thời gian tính toán.

Nghiên cứu có xét đến tương tác cơ nhiệt giữa các hạt độn lân cận hay không? Có, nghiên cứu sử dụng các công thức vi cơ học tiên tiến để xác định mô đun thể tích $K$, mô đun trượt $G$ và hệ số giãn nở nhiệt $\alpha$, trong đó đã bao hàm đầy đủ sự tương tác giữa pha nền polymer PVC và các hạt cầu Titan phân bố đều.

Kết luận

  • Thiết lập thành công hệ phương trình cơ bản và điều kiện biên uốn tấm composite mỏng cốt hạt dưới tác động của trường nhiệt đối lưu dừng và không dừng.
  • Tìm ra nghiệm giải tích tường minh dạng chuỗi lượng giác Fourier kép cho trường nhiệt độ và độ võng uốn của tấm composite hạt Titan nền PVC.
  • Xác định tỷ lệ thể tích hạt độn tối ưu từ 10% đến 20%, giúp tăng độ cứng uốn và giảm hơn 40% độ võng nhiệt của tấm.
  • Chứng minh tính tiệm cận của quá trình truyền nhiệt không dừng về trạng thái dừng khi thời gian đạt ngưỡng $1200\text{ s}$.
  • Đóng góp bộ mã nguồn tính toán số hoàn chỉnh trên Matlab 7.1, tạo tiền đề mở rộng nghiên cứu sang kết cấu vỏ composite FGM và tấm nhiều lớp trong giai đoạn tiếp theo.

Các kỹ sư và nhà nghiên cứu có thể ứng dụng ngay hệ công thức giải tích trong luận văn để nâng cao chất lượng thiết kế kết cấu composite chịu nhiệt.