Tổng quan nghiên cứu

Ngành công nghiệp chế tạo hiện đại ghi nhận mức tăng trưởng ứng dụng vật liệu composite vượt hơn 15% mỗi năm nhờ các đặc tính vượt trội về tỷ trọng nhẹ, độ bền riêng cao và khả năng kháng ăn mòn. Tuy nhiên, trong các kết cấu kỹ thuật thực tế như vỏ máy bay, thân tàu thủy hay sàn chịu lực, vật liệu dạng tấm thường xuyên làm việc trong môi trường nhiệt độ biến đổi từ 30°C đến 100°C kết hợp tải trọng cơ học. Sự tương tác phức tạp giữa trường nhiệt và trường biến dạng khiến cấu trúc tấm mỏng và tấm trung bình có tỷ số bề dày trên chiều dài cạnh từ 1/100 đến 1/5 dễ bị võng lớn, suy giảm 20% đến 40% khả năng chịu tải và đối mặt với nguy cơ tách lớp vật liệu.

Mục tiêu trọng tâm của luận văn là nghiên cứu, mô hình hóa ứng xử cơ học của tấm composite nhiều lớp và xây dựng thuật toán tính toán chuyển vị, độ võng dưới tác dụng đồng thời của tải trọng và nhiệt độ. Nghiên cứu tập trung giải quyết bài toán biến dạng đàn hồi tuyến tính của vật liệu composite đẳng hướng ngang và trực hướng chịu uốn tĩnh trong môi trường nhiệt.

Được thực hiện trong khuôn khổ chương trình đào tạo thạc sĩ chuyên ngành Công nghệ chế tạo máy tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh giai đoạn 2013-2015 dưới sự hướng dẫn khoa học của Phó Giáo sư Nguyễn Hoài Sơn, luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc. Kết quả nghiên cứu cung cấp công cụ mô phỏng số tự lập trình giúp giảm khoảng 35% thời gian phân tích kết cấu sơ bộ, đạt độ tin cậy cao với sai số kiểm chứng dưới 3% so với các nghiệm giải tích kinh điển, góp phần nâng cao hiệu quả thiết kế và dự đoán tuổi thọ làm việc của chi tiết cơ khí.

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

Khung lý thuyết áp dụng

Luận văn xây dựng hệ thống tính toán dựa trên sự tích hợp của hai nền tảng lý thuyết cơ học cốt lõi và mô hình ứng suất phẳng hai chiều. Đầu tiên là lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin dành cho kết cấu tấm. Khác với lý thuyết tấm mỏng kinh điển của Kirchhoff vốn bỏ qua biến dạng cắt ngang, mô hình Mindlin giả thiết các đoạn thẳng vuông góc với mặt trung hòa vẫn thẳng sau biến dạng nhưng không nhất thiết phải duy trì phương vuông góc. Điều này cho phép tính toán chính xác các thành phần biến dạng trượt ngang gamma xz và gamma yz, mang ý nghĩa quyết định đối với vật liệu composite có mô đun trượt tương đối thấp.

Nền tảng thứ hai là lý thuyết tấm composite nhiều lớp liên tục. Trạng thái cơ nhiệt của từng lớp vật liệu lamina được mô tả qua quan hệ ứng suất và biến dạng mở rộng, tích hợp vector biến dạng nhiệt ban đầu tỷ lệ với hệ số giãn nở nhiệt alpha và độ biến thiên nhiệt độ Delta T. Các lớp vật liệu dị hướng được chuyển đổi từ hệ trục chính của sợi sang hệ trục tọa độ toàn cục của tấm thông qua ma trận quay biến dạng T epsilon. Cấu trúc tổng thể của tấm được đặc trưng bởi các đại lượng quan trọng: ma trận độ cứng màng A, ma trận liên kết tương tác màng - uốn - xoắn B, ma trận độ cứng uốn D và ma trận độ cứng cắt A'.

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

Nguồn dữ liệu nghiên cứu được thiết lập dựa trên quá trình mô phỏng số 4 bài toán kỹ thuật điển hình với các thông số vật liệu thực tế như sợi thủy tinh, sợi carbon kết hợp nền nhựa epoxy. Cỡ mẫu tính toán được rời rạc hóa trên mô hình lưới có cấu trúc gồm 100 phần tử và 121 điểm nút tính toán. Phương pháp chọn mẫu lưới phần tử hữu hạn đồng nhất theo ma trận 10x10 được lựa chọn nhằm đảm bảo sự cân bằng tối ưu giữa độ chính xác của trường chuyển vị 5 bậc tự do tại mỗi nút và chi phí giải thuật ma trận.

Phương pháp phân tích chủ đạo là phương pháp phần tử hữu hạn kết hợp tích phân số Gauss 2 điểm và hàm nội suy Lagrange bậc nhất trên phần tử tứ giác 4 nút Q4. Lý do lựa chọn phương pháp phần tử hữu hạn là khả năng xử lý xuất sắc các phương trình vi phân đạo hàm riêng của bài toán cơ học liên tục với điều kiện biên hình học phức tạp. Tác giả đã tự xây dựng toàn bộ mã nguồn thuật toán trên ngôn ngữ lập trình Matlab để thiết lập ma trận độ cứng tổng thể K, vector tải trọng cơ nhiệt P và giải hệ phương trình đại số tuyến tính tìm trường chuyển vị u, v, w. Toàn bộ quy trình nghiên cứu được triển khai chặt chẽ theo timeline 24 tháng, bao gồm các giai đoạn tổng hợp lý thuyết, lập trình giải thuật số, mô phỏng kiểm chứng và hiệu chuẩn sai số.

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 tính toán chi tiết và rút ra 4 phát hiện quan trọng về ứng xử cơ nhiệt của tấm composite:

Thứ nhất, nhiệt độ môi trường là tác nhân khuếch đại mạnh mẽ độ võng của tấm. Khi nhiệt độ làm việc tăng từ T = 0°C lên T = 100°C dưới cùng một mức tải trọng cơ học phân bố đều, độ võng cực đại w tại tâm tấm tăng từ 38% đến 65% tùy thuộc vào cấu hình vật liệu, phản ánh sự suy giảm rõ rệt của độ cứng chống uốn khi bị kích nhiệt.

Thứ hai, tỷ lệ thể tích pha cốt sợi Vf đóng vai trò then chốt trong việc kiểm soát biến dạng. Khi tăng hệ số thể tích sợi Vf từ mức 0.2 lên 0.6, độ võng của tấm composite giảm khoảng 45% ở điều kiện không nhiệt và giảm hơn 50% ở môi trường nhiệt độ T = 100°C, khẳng định sợi chịu lực là thành phần chủ lực kiềm chế sự giãn nở nhiệt của pha nền polyme.

Thứ ba, cấu trúc số lớp xếp chồng ảnh hưởng quyết định đến độ ổn định kết cấu. So sánh giữa tấm composite 2 lớp và tấm 8 lớp có cùng tổng chiều dày, kết cấu 8 lớp đối xứng giúp tăng mô men chống uốn tương đương lên 2.4 lần, làm giảm độ võng cực đại tới 52% nhờ triệt tiêu được các thành phần của ma trận liên kết uốn - màng B.

Thứ tư, quy luật phân bố trường chuyển vị không gian cho thấy độ võng lớn nhất luôn tập trung tại tâm tấm, trong khi độ võng tại vị trí biên giữa chỉ bằng khoảng 42% đến 48% giá trị cực đại, và tiệm cận về mức 0 mm tại các vị trí góc ngàm cố định.

Thảo luận kết quả

Nguyên nhân cốt lõi của sự gia tăng biến dạng trong môi trường nhiệt xuất phát từ sự không đồng nhất về hệ số giãn nở nhiệt giữa sợi và nhựa nền, tạo ra ứng suất nhiệt nội tại cục bộ. Khi các biên của tấm bị hạn chế chuyển vị bởi điều kiện ngàm, lực nén nhiệt sinh ra kết hợp với tải trọng ngang làm gia tăng mô men uốn tổng cộng. Hiện tượng này càng trở nên nghiêm trọng đối với các tấm composite có cấu trúc lớp bất đối xứng do xuất hiện hiệu ứng ghép uốn - màng.

So sánh với các nghiên cứu giải tích của Reddy và Panda trên cùng mô hình tấm Mindlin, kết quả tính toán từ chương trình Matlab của luận văn đạt sự tương đồng rất cao, với độ lệch tương đối chỉ dao động trong khoảng từ 1.8% đến 2.6%. Điều này chứng minh thuật toán số và ma trận phần tử tứ giác 4 nút được xây dựng có độ chính xác xuất sắc.

Toàn bộ tập dữ liệu nghiên cứu được biểu diễn chi tiết thông qua hệ thống 6 bảng số liệu từ Bảng 5.1 đến Bảng 5.6 và 6 đồ thị đặc tính từ Đồ thị 5.1 đến Đồ thị 5.6. Dữ liệu có thể được trình bày trực quan bằng biểu đồ tương quan đường cong giữa tỷ lệ sợi Vf và độ võng, hoặc bản đồ mật độ chuyển vị 3D, giúp các kỹ sư nhận diện trực quan vùng tập trung ứng suất và đưa ra phương án gia cường hợp lý.

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

Dựa trên kết quả phân tích ứng xử cơ nhiệt, luận văn đưa ra 4 khuyến nghị và giải pháp kỹ thuật có tính ứng dụng cao:

Thiết kế cấu hình lớp đối xứng: Áp dụng quy tắc xếp lớp đối xứng với số lớp tối thiểu từ 8 lớp trở lên cho các kết cấu composite làm việc trong dải nhiệt độ từ 60°C đến 100°C. Giải pháp này giúp triệt tiêu hoàn toàn ma trận tương tác uốn - màng B, đặt mục tiêu giảm thiểu trên 40% độ võng dư trong chu kỳ thiết kế thử nghiệm kéo dài 6 tháng. Chủ thể thực hiện là kỹ sư thiết kế kết cấu cơ khí và vật liệu.

Tối ưu hóa tỷ lệ thể tích sợi: Kiểm soát chặt chẽ tỷ lệ thể tích sợi Vf trong khoảng tối ưu từ 55% đến 65% trong quá trình gia công ép nguội hoặc đúc hút chân không. Mục tiêu nhằm nâng cao mô đun đàn hồi uốn của chi tiết thêm 25% mà không làm tăng khối lượng bản thân, thực hiện trong lộ trình cải tiến công nghệ sản xuất 9 tháng bởi phòng kỹ thuật của các nhà máy sản xuất composite.

Tích hợp công cụ tính toán số: Ứng dụng module lập trình phần tử hữu hạn trên Matlab của luận văn vào hệ thống hỗ trợ thiết kế sơ bộ nội bộ của doanh nghiệp. Mục tiêu tự động hóa 80% thời gian phân tích ứng xử nhiệt của tấm composite, hoàn thành việc chuyển giao và đóng gói phần mềm trong thời gian 12 tháng do bộ phận nghiên cứu và phát triển phối hợp với các chuyên gia tính toán số đảm nhiệm.

Mở rộng mô hình nghiên cứu phi tuyến: Triển khai các nghiên cứu mở rộng sang bài toán dao động riêng và phân tích ứng xử phi tuyến hình học khi tấm chịu tải nhiệt độ biến thiên theo thời gian từ âm 40°C đến 250°C. Kế hoạch nghiên cứu dài hạn 18 tháng do các nhóm nghiên cứu chuyên sâu tại viện nghiên cứu và trường đại học thực hiện.

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

Nội dung luận văn mang lại giá trị tham khảo thiết thực cho 4 nhóm đối tượng chính:

Học viên cao học và nghiên cứu sinh ngành Cơ kỹ thuật, Cơ học ứng dụng và Công nghệ chế tạo máy: Tài liệu cung cấp cơ sở toán học tường minh về lý thuyết tấm Mindlin, các bước thiết lập ma trận độ cứng phần tử hữu hạn và mã nguồn Matlab mẫu để phát triển các đề tài nghiên cứu chuyên sâu.

Kỹ sư thiết kế kết cấu trong ngành hàng không, đóng tàu và cơ khí chế tạo: Sử dụng các công thức chuyển đổi ứng suất, biến dạng và ma trận độ cứng để tính toán tối ưu hóa chiều dày, hướng đặt sợi và góc nghiêng của các tấm composite chịu tải nhiệt phức tạp.

Chuyên gia công nghệ vật liệu mới: Tham khảo các phân tích định lượng về sự thay đổi cơ tính theo tỷ lệ thể tích sợi và nhiệt độ môi trường, từ đó hoàn thiện công thức phối trộn vật liệu đáp ứng tiêu chuẩn làm việc khắc nghiệt.

Giảng viên các trường đại học kỹ thuật: Có thể sử dụng các ví dụ tính toán, sơ đồ phần tử và đồ thị kết quả của luận văn làm học liệu trực quan phục vụ giảng dạy các môn học Phương pháp phần tử hữu hạn, Cơ học vật liệu composite và Cơ học ứng suất phẳng.

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

Lý thuyết tấm biến dạng cắt bậc nhất Mindlin có ưu điểm gì vượt trội so với lý thuyết tấm mỏng Kirchhoff? Lý thuyết Mindlin tính đến sự thay đổi góc xoay của tiết diện ngang và biến dạng trượt cắt qua chiều dày, giúp phản ánh chính xác trạng thái chịu lực của tấm trung bình và tấm dày. Trong vật liệu composite có mô đun trượt G12 và G13 nhỏ, mô hình Mindlin cho kết quả chuyển vị tin cậy hơn với độ chính xác cao hơn 15% so với lý thuyết Kirchhoff.

Chương trình tính toán lập trình trên Matlab có ưu thế gì so với phần mềm thương mại đóng gói? Mã nguồn Matlab được xây dựng trực tiếp từ các phương trình vi phân và ma trận đàn hồi cơ sở, cho phép người dùng can thiệp linh hoạt vào từng biến số như hướng sợi, số lớp và điều kiện nhiệt. Điều này giúp tối ưu hóa thời gian tính toán cho lưới 100 phần tử xuống dưới vài giây và hỗ trợ tùy biến thuật toán dễ dàng.

Môi trường nhiệt độ ảnh hưởng như thế nào đến độ cứng của kết cấu composite? Nhiệt độ tăng cao làm giãn nở mạng lưới polyme trong pha nền, dẫn đến hiện tượng suy giảm mô đun đàn hồi E và mô đun trượt G của vật liệu. Khi tấm bị ngàm biên, sự giãn nở nhiệt bị ngăn cản sẽ sinh ra ứng suất nén nhiệt, làm giảm từ 20% đến 40% độ cứng chống uốn hiệu dụng của toàn bộ tấm.

Làm thế nào để hạn chế tối đa hiện tượng cong võng của tấm composite khi làm việc ở nhiệt độ cao? Kỹ sư nên lựa chọn cấu hình xếp lớp đối xứng để triệt tiêu ma trận liên kết uốn - màng B, đồng thời duy trì tỷ lệ thể tích sợi Vf trong khoảng 55% đến 60%. Việc tăng số lớp liên kết từ 4 lớp lên 8 lớp cũng giúp phân bố đều trường ứng suất và giảm trên 50% độ võng uốn cực đại.

Phương pháp phần tử hữu hạn trong luận văn có thể áp dụng cho các dạng điều kiện biên nào khác? Thuật toán được xây dựng tổng quát cho 5 bậc tự do tại mỗi nút, cho phép áp dụng linh hoạt cho mọi dạng điều kiện biên như tựa khớp đơn, ngàm 4 cạnh, ngàm đối xứng hoặc biên tự do thông qua việc khử hàng và cột tương ứng trong ma trận độ cứng tổng thể K.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán cơ học phức tạp về phân tích ứng xử của tấm composite nhiều lớp trong môi trường nhiệt độ, mang lại các đóng góp nổi bật:

  • Thiết lập thành công hệ thống phương trình quan hệ ứng suất, biến dạng và chuyển vị của tấm composite dựa trên lý thuyết tấm Mindlin và trạng thái ứng suất phẳng.
  • Xây dựng hoàn chỉnh chương trình tính toán số bằng ngôn ngữ Matlab sử dụng phần tử tứ giác 4 nút Q4 với độ chính xác cao và thời gian xử lý tối ưu.
  • Khảo sát định lượng tác động của nhiệt độ từ 0°C đến 100°C, chỉ ra mức độ gia tăng độ võng lên tới 65% khi chịu tải cơ nhiệt đồng thời.
  • Chứng minh hiệu quả của cấu trúc 8 lớp đối xứng và tỷ lệ thể tích sợi 60% trong việc triệt tiêu hiện tượng uốn - màng và giảm 52% độ võng cực đại.
  • Đóng góp nguồn dữ liệu tham khảo học thuật có giá trị khoa học và khả năng ứng dụng thực tiễn cao cho ngành cơ kỹ thuật tại Việt Nam.

Hướng nghiên cứu tiếp theo sẽ tập trung vào việc mở rộng bài toán sang miền phi tuyến và tải trọng động trong lộ trình 12 tháng tới. Các cá nhân, đơn vị quan tâm đến giải thuật và mô hình tính toán số có thể kết nối để cùng phát triển và ứng dụng công nghệ mô phỏng vật liệu composite tiên tiến vào thực tiễn sản xuất.