Tổng quan nghiên cứu

Trong các ngành kỹ thuật công nghệ cao như hàng không vũ trụ, công nghiệp năng lượng hạt nhân và chế tạo thiết bị chịu áp lực, kết cấu vỏ trụ mỏng phải hoạt động liên tục dưới tác động đồng thời của tải trọng cơ học phức tạp và môi trường nhiệt độ khắc nghiệt từ 300K đến 600K. Khi chịu ứng suất xoắn kết hợp gradient nhiệt độ cao, kết cấu rất dễ rơi vào trạng thái mất ổn định phi tuyến, gây hư hỏng đột ngột và tổn thất nghiêm trọng. Vật liệu composite nhiều lớp truyền thống thường bộc lộ nhược điểm nứt vỡ hoặc tách lớp do sự tập trung ứng suất tại bề mặt phân giới. Nhằm khắc phục hạn chế này, vật liệu cơ tính biến thiên (Functionally Graded Material - FGM) kết hợp cấu trúc lõi xốp kim loại (Porous Metal Foam) đã ra đời, giúp giảm khoảng 30% đến 50% trọng lượng kết cấu nhưng vẫn duy trì độ bền cơ lý và khả năng cách nhiệt tuyệt vời.

Mục tiêu trọng tâm của nghiên cứu là xây dựng lời giải giải tích tường minh để phân tích ứng xử ổn định phi tuyến và đáp ứng sau tới hạn của vỏ trụ sandwich ba lớp FGM có lõi xốp, được gia cường bởi hệ gân trực giao lệch tâm và đặt trên nền đàn hồi Pasternak hai tham số chịu tác dụng của tải xoắn trong môi trường nhiệt độ. Phạm vi nghiên cứu khảo sát các thông số hình học với bán kính $R$, chiều dày $h = 0.006$ m, tỉ số bán kính trên chiều dày $R/h$ từ 80 đến 120, tỉ số chiều dài trên bán kính $L/R$ từ 1 đến 2, được hoàn thiện tại Hà Nội vào năm 2018. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc cung cấp công thức giải tích chuẩn xác với sai số kiểm chứng dưới 2.02%, giúp các kỹ sư dự báo chính xác tải tới hạn trên ($\tau_{upper}$) và tải tới hạn dưới ($\tau_{lower}$), từ đó tối ưu hóa thiết kế vỏ mỏng chịu lực.

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 các lý thuyết cơ học kết cấu phi tuyến hiện đại kết hợp mô hình vật liệu tiên tiến:

  • Lý thuyết vỏ mỏng Donnell có xét đến tính phi tuyến hình học von Kármán: Mô tả chính xác mối quan hệ phi tuyến giữa các thành phần biến dạng và chuyển vị của vỏ trụ khi xảy ra biến dạng lớn.
  • Mô hình vật liệu FGM xốp ba lớp (FG/P/FG): Lớp phủ ngoài và lớp phủ trong làm bằng vật liệu FGM biến thiên liên tục giữa gốm Silicon nitride ($Si_3N_4$) với mô đun Young $E_c = 348.43$ GPa và thép kim loại SUS304 với $E_m = 201.04$ GPa; lớp lõi giữa là kim loại bọt có phân bố độ rỗng xốp đối xứng với hệ số độ xốp $e_0$ dao động từ 0 đến 0.5.
  • Kỹ thuật san đều tác dụng của gân Lekhnitsky: Quy đổi tác động cơ học cục bộ của hệ thống gân dọc và gân vòng lệch tâm thành độ cứng phân bố đều tương đương trên bề mặt vỏ.
  • Mô hình nền đàn hồi Pasternak hai tham số: Xét đến tương tác giữa vỏ trụ và môi trường đàn hồi bao quanh thông qua độ cứng chống lún $K_1 = 6 \times 10^7\text{ N/m}^3$ và độ cứng chống trượt $K_2 = 4 \times 10^5\text{ N/m}$.

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

Nguồn dữ liệu tính toán được thiết lập từ các thông số vật lý thực nghiệm chuẩn của hợp kim thép không gỉ SUS304 và gốm kỹ thuật $Si_3N_4$, với hệ số Poisson $\nu = 0.3$. Hệ gân gia cường gồm 20 gân dọc ($n_s = 20$) và 20 gân vòng ($n_r = 20$) có kích thước tiết diện $0.006\text{ m} \times 0.006\text{ m}$. Cỡ mẫu khảo sát bao gồm hàng trăm trường hợp biến thiên tham số hình học và cơ lý được phân tích có hệ thống nhằm đánh giá toàn diện mọi yếu tố ảnh hưởng.

Phương pháp giải tích Galerkin hai chiều được lựa chọn làm phương pháp phân tích chủ đạo vì cho phép thiết lập biểu thức hiển xác định tải tới hạn và đường cong quan hệ tải trọng - độ võng lớn nhất một cách tường minh, tiết kiệm thời gian tính toán hơn so với phương pháp phần tử hữu hạn số mà vẫn đạt độ tin cậy tuyệt đối. Hàm độ võng được giả thiết dạng sóng điều hòa thỏa mãn điều kiện biên tựa đơn ở hai đầu vỏ trụ và điều kiện chu vi kín. Tiến độ nghiên cứu được triển khai đồng bộ từ khâu thiết lập phương trình vi phân cân bằng, xây dựng thuật toán giải tích cho đến mô phỏng số kiểm chứng trong giai đoạn từ năm 2017 đến năm 2018 tại Khoa Toán - Cơ - Tin học, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội.

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

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

  1. Ảnh hưởng của độ xốp và tỉ số chiều dày lõi: Khi tăng hệ số độ xốp $e_0$ từ 0 lên 0.5 và tăng tỉ số chiều dày lõi trên lớp phủ $h_{core}/h_{FG}$ từ 0 lên 2, tải xoắn tới hạn trên $\tau_{upper}$ của vỏ có gân trong giảm mạnh khoảng 52%, từ 1174.96 MPa xuống còn khoảng 773.45 MPa.
  2. Ảnh hưởng của gradient nhiệt độ: Khi nhiệt độ môi trường $\Delta T$ tăng từ 0K lên 300K và 600K, tải xoắn tới hạn dưới $\tau_{lower}$ của vỏ giảm lần lượt khoảng 5.17 MPa và 10.34 MPa do ứng suất nhiệt trước nén làm mềm kết cấu.
  3. Hiệu quả của vị trí và cấu hình gân gia cường: Vỏ trụ bố trí gân ngoài luôn cho khả năng chịu tải xoắn cao hơn khoảng 15% đến 25% so với vỏ bố trí gân trong. Với cùng số lượng gân gia cường, vỏ có gân vòng chịu xoắn tốt nhất, tiếp theo là hệ gân trực giao, gân dọc và thấp nhất là vỏ không có gân.
  4. Tác động tích cực của nền đàn hồi Pasternak: Khi có sự xuất hiện của nền Pasternak với $K_1 = 6 \times 10^7\text{ N/m}^3$ và $K_2 = 4 \times 10^5\text{ N/m}$, tải xoắn tới hạn của vỏ tăng hơn 35% so với khi không đặt trên nền đàn hồi.

Thảo luận kết quả

Sự sụt giảm tải tới hạn khi tăng hệ số độ xốp $e_0$ và tỉ số $h_{core}/h_{FG}$ bắt nguồn từ việc mật độ vật liệu trong lớp lõi giảm, làm suy giảm trực tiếp mô đun đàn hồi Young tương đương và độ cứng chống xoắn của vỏ trụ. Khi nhiệt độ tăng, các thành phần ứng suất nhiệt xuất hiện gây biến dạng phồng ban đầu (độ võng âm) ngay cả khi chưa có lực xoắn cơ học tác dụng.

Khi kiểm chứng với các công bố khoa học quốc tế của Shen và Nash đối với vỏ trụ thuần nhất không gân, kết quả tính toán từ công thức giải tích của luận văn chỉ sai lệch từ 1.01% đến 2.02%, khẳng định độ tin cậy vượt trội của mô hình toán. Toàn bộ diễn biến làm việc sau tới hạn được minh họa rõ nét thông qua các đồ thị quan hệ tải xoắn - độ võng không thứ nguyên ($\tau - W_{max}/h$) và tải xoắn - góc xoắn ($\tau - \psi$). Các đường cong này chỉ ra rằng sau khi đạt điểm tới hạn trên $\tau_{upper}$, kết cấu xảy ra hiện tượng nhảy bước mất ổn định (snap-through) chuyển sang trạng thái cân bằng mới ở tải tới hạn dưới $\tau_{lower}$.

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

  1. Tối ưu hóa tỉ lệ thành phần gốm và độ xốp lõi: Khống chế hệ số độ xốp $e_0 \le 0.3$ và duy trì chỉ số tỉ phần thể tích $k \ge 1.0$ trong quá trình luyện kim bột để giữ hàm lượng gốm $Si_3N_4$ ở mức tối ưu, đảm bảo tải xoắn tới hạn trên luôn đạt trên 950 MPa. Chủ thể thực hiện: Các kỹ sư công nghệ vật liệu và chế tạo composite. Thời gian: Giai đoạn thiết kế thành phần vật liệu.
  2. Ưu tiên cấu hình gân gia cường mặt ngoài: Bố trí hệ gân vòng hoặc gân trực giao lệch tâm ở bề mặt ngoài của vỏ trụ nhằm nâng cao hiệu quả kháng mô-men xoắn lên ít nhất 20% so với cấu hình gân trong với cùng khối lượng vật liệu tiêu hao. Chủ thể thực hiện: Kỹ sư thiết kế kết cấu thân vỏ máy bay và tên lửa. Thời gian: 1 đến 3 tháng trong pha thiết kế kỹ thuật.
  3. Tích hợp đệm đàn hồi giảm chấn: Thiết kế bổ sung các lớp vật liệu chèn đàn hồi có độ cứng $K_1 \ge 6 \times 10^7\text{ N/m}^3$ xung quanh vỏ trụ để triệt tiêu biến dạng nhiệt và tăng độ cứng vững công trình lên trên 30%. Chủ thể thực hiện: Đơn vị thi công và lắp đặt công trình ngầm, ống dẫn áp lực. Thời gian: Pha thi công chế tạo.
  4. Phát triển phần mềm tính toán tự động: Xây dựng module tính toán tự động dựa trên các biểu thức giải tích tường minh của nghiên cứu để cắt giảm 70% thời gian phân tích tiền khả thi so với mô phỏng phần tử hữu hạn 3D thông thường. Chủ thể thực hiện: Nhóm nghiên cứu tính toán cơ học và doanh nghiệp phát triển phần mềm CAE. Thời gian: Triển khai trong 6 tháng.

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

  1. Học viên cao học và nghiên cứu sinh ngành Cơ học vật rắn: Nắm vững quy trình thiết lập phương trình phi tuyến Donnell, kỹ thuật san đều gân Lekhnitsky và thuật toán Galerkin trong bài toán tấm vỏ phức hợp.
  2. Kỹ sư thiết kế kết cấu hàng không vũ trụ và đóng tàu: Áp dụng các công thức giải tích tính toán nhanh sức bền và giới hạn ổn định xoắn của thân tên lửa, vỏ máy bay hoặc thân tàu ngầm composite.
  3. Các nhà khoa học phát triển vật liệu tiên tiến: Tham khảo mô hình phân bố cơ tính biến thiên và vi cấu trúc xốp kim loại để chế tạo các dòng vật liệu sandwich chịu nhiệt siêu nhẹ thế hệ mới.
  4. Giảng viên và chuyên gia cơ kỹ thuật: Sử dụng làm tài liệu giảng dạy chuyên đề sau đại học về cơ học vật liệu FGM, lý thuyết ổn định phi tuyến và phương pháp giải tích trong cơ học kết cấu.

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

Vật liệu FGM sandwich có lõi xốp mang lại ưu điểm vượt trội gì so với composite truyền thống? Vật liệu FGM có cơ tính biến thiên liên tục giúp triệt tiêu hoàn toàn sự tập trung ứng suất tiếp tuyến tại mặt phân giới, loại bỏ nguy cơ bong tách lớp. Lõi xốp kim loại giúp giảm khoảng 30% đến 50% khối lượng bản thân mà vẫn đảm bảo độ bền cơ học và cách nhiệt tốt.

Tại sao cần phân tích ổn định phi tuyến thay vì chỉ dùng lý thuyết tuyến tính? Dưới tác dụng của tải trọng xoắn lớn và nhiệt độ cao, vỏ trụ mỏng phát sinh các chuyển vị lớn vượt ra ngoài phạm vi tuyến tính. Lý thuyết phi tuyến von Kármán cho phép xác định chính xác hiện tượng biến dạng sau tới hạn và tải giới hạn dưới, tránh đánh giá quá cao độ an toàn của kết cấu.

Yếu tố nhiệt độ tác động như thế nào đến khả năng chịu xoắn của vỏ trụ? Sự gia tăng nhiệt độ $\Delta T$ gây ra ứng suất nén nhiệt trước và làm suy giảm mô đun đàn hồi của vật liệu. Khi nhiệt độ tăng 300K, tải tới hạn dưới của vỏ có thể giảm khoảng 5.17 MPa và làm vỏ bị biến dạng uốn phồng ban đầu trước cả khi có lực xoắn cơ học tác dụng.

Bố trí gân gia cường ở mặt ngoài hay mặt trong giúp tăng khả năng chịu xoắn tốt hơn? Kết quả nghiên cứu giải tích khẳng định gân bố trí ở mặt ngoài có khoảng cách xa trục trung hòa hơn, tạo ra mô-men kháng uốn và kháng xoắn lớn hơn, giúp nâng cao tải tới hạn của vỏ thêm khoảng 15% đến 25% so với khi bố trí gân ở mặt trong.

Phương pháp giải tích trong luận văn có độ tin cậy như thế nào khi so sánh với thực nghiệm? Công thức giải tích của luận văn khi kiểm chứng đối chiếu với các kết quả nghiên cứu độc lập của Shen và Nash cho thấy sai số chỉ dao động từ 1.01% đến 2.02%, đảm bảo độ chuẩn xác cao và có thể ứng dụng trực tiếp trong thực tế kỹ thuật.

Kết luận

  • Luận văn đã thiết lập thành công hệ phương trình giải tích tường minh xác định chính xác tải xoắn tới hạn trên, tải tới hạn dưới và đáp ứng sau tới hạn của vỏ trụ sandwich FGM lõi xốp có gân gia cường trên nền đàn hồi Pasternak.
  • Chứng minh rõ ràng quy luật suy giảm tải tới hạn khi tăng độ xốp $e_0$, tăng tỉ số chiều dày lõi $h_{core}/h_{FG}$, tăng nhiệt độ $\Delta T$, cũng như hiệu quả gia tăng độ cứng vượt trội của hệ gân ngoài và nền Pasternak.
  • Khẳng định tính đúng đắn và độ tin cậy của thuật toán giải tích với sai số so sánh chỉ dưới 2.02% so với các công bố quốc tế uy tín.
  • Kế hoạch phát triển tiếp theo trong 12 tháng tới bao gồm việc mở rộng bài toán sang tải trọng động lực học, xét đến hiệu ứng biến dạng trượt bậc cao và hình học vỏ nón cụt.
  • Quý bạn đọc, nhà nghiên cứu và kỹ sư quan tâm có thể khai thác trực tiếp các biểu thức giải tích trong công trình để ứng dụng vào thiết kế tối ưu các kết cấu vỏ mỏng chịu lực nhiệt phức tạp.