Luận văn ThS Nguyễn Tiến Thiệp: Phân tích ổn định vỏ trụ sandwich cơ tính biến thiên

Luận văn: Phân tích ổn định phi tuyến vỏ trụ sandwich gân gia cường chịu xoắn trong môi trường nhiệt. Nghiên cứu cơ tính biến thiên.

• Tác giả: Nguyễn Tiến Thiệp
• Người hướng dẫn: Ts. Lê Khả Hõa
• Trường học: Trường Đại Học Khoa Học Tự Nhiên - ĐHQGHN
• Chuyên ngành: Cơ Học Vật Rắn
• Đề tài: Phân Tích Ổn Định Phi Tuyến Của Vỏ Trụ Sandwich Cơ Tính Biến Thiên Có Gân Gia Cường Trực Giao Chịu Xoắn Trong Môi Trƣờng Nhiệt Độ
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2018
• Địa điểm: Hà Nội

Mục lục chi tiết

• Danh sách hình vẽ
• Danh sách bảng vẽ
• MỞ ĐẦU
• 1. Chương 1: Vỏ trụ và các hệ thức cơ bản
• 1.1. Vỏ trụ sandwich trên nền đàn hồi Pasternak
• 1.2. Các hệ thức cơ bản và phương trình cân bằng
• 2. Chương 2: Phương pháp giải
• 2.1. Điều kiện biên và dạng nghiệm
• 3. Chương 3: Tính toán số
• 3.1. So sánh kết quả
• 3.2. Kết quả số cho vỏ nón cụt ES – FGM
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Khám phá tổng quan về ổn định phi tuyến vỏ trụ FGM

Lĩnh vực cơ học kết cấu hiện đại chứng kiến sự phát triển vượt bậc của các vật liệu tiên tiến, trong đó vật liệu phân cấp chức năng (FGM) nổi lên như một giải pháp đột phá. Các kết cấu FGM, đặc biệt là kết cấu sandwich FGM, được ứng dụng rộng rãi trong các ngành công nghệ cao như hàng không vũ trụ, năng lượng hạt nhân, nơi vật liệu phải chịu đựng môi trường nhiệt độ khắc nghiệt và tải trọng phức tạp. Ưu điểm chính của FGM là sự biến đổi tính chất cơ lý một cách trơn tru, liên tục giữa các bề mặt, giúp loại bỏ hiện tượng tập trung ứng suất, bong tách và nứt gãy thường thấy ở các vật liệu composite lớp truyền thống. Vỏ trụ là một dạng kết cấu cơ bản và phổ biến, do đó, việc nghiên cứu về ổn định, dao động và độ bền của vỏ trụ FGM thu hút sự quan tâm lớn. Phân tích ổn định phi tuyến vỏ trụ sandwich FGM chịu xoắn là một bài toán cốt lõi, giúp hiểu rõ hành vi của kết cấu khi chịu tác động của tải trọng xoắn, một trong những nguyên nhân chính gây ra mất ổn định kết cấu. Nghiên cứu này không chỉ có ý nghĩa về mặt học thuật mà còn mang lại giá trị thực tiễn to lớn, góp phần thiết kế các công trình an toàn, bền vững và tối ưu hơn về vật liệu cũng như chi phí.

1.1. Tìm hiểu vật liệu phân cấp chức năng FGM tiên tiến

Vật liệu phân cấp chức năng (Functionally Graded Material - FGM) là một loại composite cao cấp, trong đó thành phần và cấu trúc vi mô thay đổi dần dần theo không gian, dẫn đến sự biến thiên liên tục của các đặc tính vật liệu. Thông thường, vật liệu FGM được chế tạo từ hỗn hợp gốm và kim loại. Lớp bề mặt giàu gốm có khả năng chịu nhiệt độ cao và chống ăn mòn tuyệt vời, trong khi lớp bề mặt giàu kim loại đảm bảo độ bền cơ học và độ dẻo dai. Sự chuyển tiếp mượt mà này giúp giảm thiểu ứng suất nhiệt và ứng suất dư, nâng cao đáng kể tuổi thọ và độ tin cậy của kết cấu. Trong luận văn của Nguyễn Tiến Thiệp (2018), vật liệu FGM được sử dụng cho các lớp phủ của vỏ trụ sandwich, tạo thành một vỏ trụ composite có khả năng làm việc ưu việt trong môi trường nhiệt.

1.2. Đặc điểm nổi bật của kết cấu sandwich FGM chịu xoắn

Kết cấu sandwich FGM thường bao gồm hai lớp phủ FGM cứng bên ngoài và một lớp lõi nhẹ hơn ở giữa. Lớp lõi này có thể làm từ kim loại, gốm hoặc các vật liệu xốp như kết cấu lõi tổ ong hoặc kim loại bọt. Cấu trúc này mang lại tỷ lệ độ cứng trên trọng lượng rất cao. Các lớp phủ FGM chịu phần lớn các tải trọng uốn và trong mặt phẳng, trong khi lớp lõi có nhiệm vụ chống lại các lực cắt ngang và duy trì khoảng cách ổn định giữa hai lớp phủ. Khi chịu tải trọng xoắn, kết cấu này thể hiện hành vi phi tuyến phức tạp, đòi hỏi các phương pháp phân tích chuyên sâu để xác định tải trọng tới hạn—ngưỡng mà tại đó kết cấu bắt đầu mất ổn định và xuất hiện các ứng suất và biến dạng lớn.

II. Thách thức phân tích mất ổn định vỏ trụ FGM chịu xoắn

Việc phân tích mất ổn định kết cấu của vỏ trụ sandwich FGM là một bài toán phức tạp trong lĩnh vực cơ học kết cấu. Sự phức tạp này đến từ nhiều yếu tố. Thứ nhất, tính chất vật liệu không đồng nhất và biến thiên liên tục theo chiều dày đòi hỏi các mô hình toán học phức tạp hơn so với vật liệu đẳng hướng truyền thống. Thứ hai, hành vi phi tuyến hình học, đặc biệt là khi kết cấu chịu biến dạng lớn, phải được xem xét để mô tả chính xác hiện tượng sau mất ổn định. Thêm vào đó, các yếu tố thực tế như môi trường nhiệt độ cao, sự hiện diện của gân gia cường, và tương tác với nền đàn hồi càng làm tăng độ khó của bài toán. Các phương pháp giải tích truyền thống thường gặp giới hạn, đòi hỏi sự kết hợp với các công cụ mô phỏng số ANSYS/ABAQUS hoặc các phương pháp giải gần đúng hiệu quả. Việc xác định chính xác tải trọng tới hạn và đường cong hậu uốn (post-buckling) là mục tiêu cuối cùng, cung cấp thông tin quan trọng cho quá trình thiết kế và đánh giá an toàn kết cấu trong các ứng dụng kỹ thuật thực tế.

2.1. Phân tích ảnh hưởng của tải trọng xoắn đến kết cấu

Tải trọng xoắn là một trong những loại tải trọng nguy hiểm nhất đối với các kết cấu vỏ mỏng, vì nó có thể gây ra hiện tượng mất ổn định oằn (buckling) ở mức tải trọng thấp hơn nhiều so với tải nén dọc trục hoặc áp lực ngoài. Dưới tác dụng của momen xoắn, các đường sinh ban đầu của vỏ trụ sẽ biến dạng thành các đường xoắn ốc. Khi tải trọng đạt đến giá trị tới hạn, vỏ sẽ đột ngột thay đổi hình dạng cân bằng, xuất hiện các vết lõm hình thoi xiên, dẫn đến phá hủy kết cấu. Việc phân tích chính xác dao động phi tuyến của vỏ dưới tải xoắn đòi hỏi phải xem xét cả các hiệu ứng phi tuyến hình học và vật liệu.

2.2. Yếu tố môi trường nhiệt và gia cường trong cơ học kết cấu

Trong nhiều ứng dụng thực tế, vỏ trụ FGM làm việc trong môi trường có nhiệt độ thay đổi. Sự thay đổi nhiệt độ gây ra ứng suất nhiệt, làm thay đổi độ cứng của vật liệu và có thể làm giảm đáng kể khả năng chịu tải của kết cấu. Ngoài ra, để tăng cường độ cứng và khả năng chịu lực, người ta thường sử dụng các gân gia cường dọc và/hoặc gân vòng. Sự hiện diện của các gân này làm cho kết cấu trở thành không đồng nhất về mặt hình học, đòi hỏi phải áp dụng các kỹ thuật mô hình hóa đặc biệt như phương pháp san đều tác dụng của gân. Tất cả những yếu tố này phải được tích hợp vào mô hình phân tích ổn định phi tuyến để có được kết quả chính xác và đáng tin cậy.

III. Phương pháp xây dựng mô hình lý thuyết vỏ ổn định Donnell

Để giải quyết bài toán ổn định phi tuyến vỏ trụ sandwich FGM chịu xoắn, việc lựa chọn một lý thuyết vỏ phù hợp là vô cùng quan trọng. Lý thuyết ổn định Donnell là một công cụ mạnh mẽ và được sử dụng rộng rãi cho các vỏ mỏng đàn hồi. Mặc dù có những giả thiết đơn giản hóa, lý thuyết này vẫn cho kết quả đủ chính xác đối với các trạng thái mất ổn định có số sóng lớn theo chu vi, vốn là đặc trưng của vỏ trụ chịu xoắn. Luận văn của Nguyễn Tiến Thiệp (2018) đã phát triển mô hình dựa trên lý thuyết vỏ Donnell cải tiến, có tính đến các yếu tố phi tuyến hình học theo von Karman. Các phương trình cơ bản của lý thuyết này bao gồm các phương trình cân bằng lực và momen, các hệ thức biến dạng-chuyển vị, và phương trình tương thích biến dạng. Bằng cách kết hợp các hệ thức này với định luật Hooke cho vật liệu FGM trong môi trường nhiệt, một hệ phương trình chủ đạo đã được thiết lập. Hệ phương trình này mô tả mối quan hệ phức tạp giữa hàm ứng suất và độ võng của vỏ, là nền tảng cho việc phân tích động lực học và ổn định của kết cấu.

3.1. Các hệ thức cơ bản và phương trình cân bằng cho vỏ trụ

Mô hình toán học bắt đầu bằng việc thiết lập các hệ thức cơ bản. Các thành phần biến dạng tại một điểm bất kỳ trong vỏ được biểu diễn qua các thành phần chuyển vị (u, v, w) trên mặt giữa của vỏ. Quan hệ biến dạng-chuyển vị này bao gồm các số hạng phi tuyến theo von Karman để nắm bắt được hành vi biến dạng lớn sau khi mất ổn định. Từ các thành phần biến dạng, các thành phần ứng suất được xác định thông qua định luật Hooke, có kể đến ảnh hưởng của nhiệt độ. Cuối cùng, các phương trình cân bằng, thể hiện sự cân bằng của các nội lực và momen trên một phân tố vỏ, được thiết lập. Hai phương trình đầu tiên được thỏa mãn bằng cách đưa vào hàm ứng suất Airy, và phương trình thứ ba (theo phương pháp tuyến) trở thành phương trình chủ đạo liên hệ giữa tải trọng ngoài, nội lực, và độ võng của vỏ.

3.2. Kỹ thuật san đều cho gân và tương tác với nền đàn hồi

Để xử lý sự có mặt của các gân gia cường, kỹ thuật san đều tác dụng của gân (smeared stiffener technique) được áp dụng. Kỹ thuật này giả thiết rằng các gân được phân bố đều trên bề mặt vỏ. Độ cứng của các gân được cộng vào độ cứng của vỏ, tạo thành một kết cấu tương đương đồng nhất về mặt hình học nhưng không đồng nhất về vật liệu. Cách tiếp cận này giúp đơn giản hóa đáng kể việc giải bài toán mà vẫn giữ được độ chính xác cần thiết. Ngoài ra, mô hình còn xem xét sự tương tác của vỏ với nền đàn hồi Pasternak, được đặc trưng bởi hai tham số: một hệ số nền kiểu Winkler (phản lực tỷ lệ với độ võng) và một hệ số nền chống cắt (phản lực tỷ lệ với đạo hàm bậc hai của độ võng).

IV. Bí quyết giải bài toán ổn định phi tuyến bằng Galerkin

Sau khi thiết lập hệ phương trình chủ đạo, bước tiếp theo là tìm nghiệm giải tích. Phương pháp Galerkin là một kỹ thuật số gần đúng mạnh mẽ, thường được sử dụng trong phân tích ổn định phi tuyến. Nguyên lý của phương pháp này là chọn một dạng nghiệm gần đúng cho hàm độ võng, thỏa mãn các điều kiện biên của bài toán. Dạng nghiệm này thường là một chuỗi các hàm lượng giác với các biên độ chưa biết. Sau đó, thay dạng nghiệm này vào hệ phương trình chủ đạo, ta sẽ thu được một phần dư (residual). Phương pháp Galerkin yêu cầu phần dư này phải trực giao với các hàm cơ sở đã chọn. Điều kiện trực giao này dẫn đến một hệ phương trình đại số phi tuyến cho các biên độ chưa biết. Giải hệ phương trình này cho phép xác định được mối quan hệ giữa tải trọng tác dụng và độ võng của vỏ. Từ đó, có thể xác định tải trọng tới hạn và xây dựng đường cong hậu uốn (post-buckling), mô tả chi tiết hành vi của vỏ trụ composite sau khi mất ổn định.

4.1. Ứng dụng phương pháp Galerkin để tìm nghiệm giải tích

Trong nghiên cứu cụ thể này, hàm độ võng w được chọn dưới dạng một chuỗi ba số hạng, bao gồm một thành phần võng đều, một thành phần tuyến tính và một thành phần phi tuyến. Dạng nghiệm này có khả năng mô tả tốt hình dạng mất ổn định của vỏ trụ chịu xoắn. Sau khi thay w vào phương trình tương thích, ta giải được hàm ứng suất f. Tiếp theo, cả w và f được thay vào phương trình cân bằng theo phương pháp tuyến. Áp dụng thủ tục Galerkin trên toàn bộ miền của vỏ trụ, hệ phương trình vi phân đạo hàm riêng phức tạp ban đầu được rút gọn thành một hệ ba phương trình đại số phi tuyến đơn giản hơn, liên hệ giữa tải trọng xoắn, các biên độ độ võng và các tham số của bài toán.

4.2. Xác định tải trọng tới hạn và đường cong hậu uốn

Từ hệ phương trình đại số thu được, ta có thể phân tích trạng thái tới hạn. Tải trọng tới hạn trên (upper critical load) được xác định bằng cách tuyến tính hóa bài toán (cho biên độ võng tiến về 0). Đây là tải trọng mất ổn định theo lý thuyết tuyến tính. Tuy nhiên, hành vi thực tế của vỏ lại được mô tả bởi đường cong hậu uốn (post-buckling). Đường cong này thể hiện mối quan hệ giữa tải trọng và biên độ độ võng sau khi mất ổn định. Đối với vỏ trụ chịu xoắn, đường cong này thường cho thấy sự sụt giảm tải trọng ngay sau điểm tới hạn, chỉ ra rằng kết cấu rất nhạy cảm với các khuyết tật ban đầu. Giá trị tải trọng thấp nhất trên đường cong hậu uốn được gọi là tải trọng tới hạn dưới (lower critical load), một thông số quan trọng trong thiết kế an toàn.

V. Kết quả số nổi bật từ phân tích động lực học kết cấu

Phân tích số dựa trên các biểu thức giải tích thu được đã cung cấp những hiểu biết sâu sắc về hành vi ổn định của vỏ trụ. Các kết quả đã được kiểm chứng độ tin cậy bằng cách so sánh với các công bố trước đó của các tác giả như Shen, Nash và Ekstrom, cho thấy sự phù hợp cao. Nghiên cứu tập trung khảo sát ảnh hưởng của nhiều tham số quan trọng. Kết quả chỉ ra rằng, tải trọng tới hạn của vỏ giảm đáng kể khi hệ số độ xốp của lớp lõi tăng, hoặc khi nhiệt độ môi trường tăng. Điều này phù hợp với trực giác vật lý, vì cả hai yếu tố này đều làm giảm độ cứng tổng thể của kết cấu. Ngược lại, việc tăng chỉ số tỉ phần thể tích (tăng hàm lượng gốm) hoặc thêm gân gia cường giúp cải thiện rõ rệt khả năng chịu tải của vỏ. Phân tích cũng cho thấy các thông số hình học như tỷ số bán kính/độ dày (R/h) và chiều dài/bán kính (L/R) có ảnh hưởng mạnh mẽ đến sự ổn định, nhấn mạnh tầm quan trọng của việc tối ưu hóa thiết kế hình học trong thực tế. Những kết quả này là cơ sở quan trọng cho việc thiết kế và đánh giá vỏ mỏng đàn hồi trong các ứng dụng kỹ thuật.

5.1. Ảnh hưởng của độ xốp và nhiệt độ đến ứng suất và biến dạng

Các tính toán số trong luận văn cho thấy rõ: khi hệ số độ xốp e0 tăng, mô đun đàn hồi của lớp lõi giảm, dẫn đến giảm độ cứng chung và làm giảm tải trọng tới hạn τcr. Tương tự, khi nhiệt độ ΔT tăng, ứng suất nhiệt nén sinh ra trong vỏ làm giảm khả năng chống lại tải trọng xoắn từ bên ngoài. Ví dụ, khi ΔT tăng từ 0K lên 600K, tải trọng tới hạn dưới giảm khoảng 10.3 MPa. Những phân tích này giúp lượng hóa mức độ ảnh hưởng của các yếu tố vật liệu và môi trường lên ứng suất và biến dạng của kết cấu.

5.2. Tác động của thông số hình học R h L R lên vỏ trụ

Thông số hình học đóng vai trò quyết định đến sự ổn định. Kết quả cho thấy khi tỷ số R/h (bán kính/độ dày) tăng, vỏ trở nên 'mỏng' hơn và khả năng chịu lực giảm đi đáng kể. Tương tự, khi tỷ số L/R (chiều dài/bán kính) tăng, vỏ trở nên 'dài' hơn và cũng dễ bị mất ổn định hơn dưới tác dụng của tải trọng xoắn. Việc tối ưu các tỷ số này là một bài toán quan trọng trong thiết kế để đảm bảo kết cấu vừa nhẹ vừa đủ bền.

5.3. Vai trò của gân gia cường và nền đàn hồi Pasternak

Sự hiện diện của gân gia cường và nền đàn hồi Pasternak đều làm tăng đáng kể tải trọng tới hạn của vỏ. Các kết quả so sánh chỉ ra rằng vỏ được gia cường bởi gân vòng có khả năng chịu xoắn tốt nhất, tiếp theo là gân trực giao và gân dọc. Vỏ không có gân gia cường và không đặt trên nền đàn hồi có khả năng chịu tải kém nhất. Điều này khẳng định hiệu quả của các giải pháp gia cường trong việc nâng cao sự ổn định của kết cấu sandwich FGM.

VI. Kết luận và định hướng tương lai cho kết cấu FGM chịu xoắn

Nghiên cứu về ổn định phi tuyến vỏ trụ sandwich FGM chịu xoắn đã thành công trong việc xây dựng một phương pháp giải tích hiệu quả, dựa trên lý thuyết ổn định Donnell và phương pháp Galerkin. Các biểu thức giải tích tường minh cho phép xác định tải trọng tới hạn và khảo sát đáp ứng sau tới hạn một cách nhanh chóng và chính xác. Các kết quả số đã lượng hóa được ảnh hưởng của hàng loạt tham số vật liệu, hình học và môi trường, cung cấp những kết luận khoa học có giá trị. Tải trọng tới hạn giảm khi độ xốp và nhiệt độ tăng, nhưng tăng lên khi chỉ số thể tích và các yếu tố gia cường được thêm vào. Hướng nghiên cứu trong tương lai có thể được mở rộng sang các bài toán phức tạp hơn. Ví dụ, phân tích vỏ chịu tải trọng động, hoặc nghiên cứu các loại kết cấu vỏ khác như vỏ nón, vỏ cầu. Ngoài ra, việc áp dụng các lý thuyết vỏ bậc cao (HSDT), có tính đến biến dạng cắt ngang, sẽ cho kết quả chính xác hơn đối với các loại vỏ dày và trung bình, mở ra những triển vọng mới trong lĩnh vực cơ học kết cấu.

6.1. Tóm tắt các phát hiện chính về dao động phi tuyến của vỏ

Luận văn đã rút ra những kết luận quan trọng: (i) Tải trọng tới hạn giảm khi hệ số độ xốp e0 và nhiệt độ ΔT tăng. (ii) Khả năng chịu tải tăng khi chỉ số tỉ phần thể tích k của lớp phủ FGM tăng. (iii) Vỏ trở nên kém ổn định hơn khi các tỷ số R/h và L/R tăng. (iv) Sự có mặt của gân gia cường và nền đàn hồi làm tăng đáng kể khả năng chịu tải, trong đó gân vòng cho hiệu quả chống xoắn tốt nhất. Những phát hiện này cung cấp các chỉ dẫn thiết thực cho việc thiết kế và tối ưu hóa các vỏ trụ composite trong thực tế.

6.2. Hướng nghiên cứu mở rộng cho vỏ nón và tải trọng động

Những thành công của nghiên cứu này mở đường cho các hướng phát triển tiếp theo. Một hướng đi tự nhiên là áp dụng phương pháp luận tương tự để giải quyết bài toán ổn định của các kết cấu vỏ phức tạp hơn như vỏ nón cụt hoặc vỏ cầu làm bằng vật liệu FGM xốp. Một hướng khác đầy hứa hẹn là phân tích động lực học, nghiên cứu phản ứng của vỏ dưới tác dụng của các tải trọng thay đổi theo thời gian. Cuối cùng, việc chuyển từ lý thuyết vỏ cổ điển sang các lý thuyết vỏ bậc cao (HSDT) sẽ cho phép mô hình hóa chính xác hơn các hiệu ứng biến dạng cắt, đặc biệt quan trọng đối với các vỏ không quá mỏng.

Nội dung chính

Tổng quan nghiên cứu

Trong bối cảnh công nghiệp hàng không vũ trụ, năng lượng nguyên tử và kỹ thuật quân sự thế kỷ 21, các kết cấu vỏ mỏng làm việc trong môi trường tải trọng cơ - nhiệt phức tạp đối mặt với nguy cơ mất ổn định cơ học nghiêm trọng. Các kết cấu truyền thống thường xảy ra hiện tượng tập trung ứng suất tại mặt tiếp xúc giữa các lớp vật liệu, dẫn đến rạn nứt hoặc bong tách kết cấu khi nhiệt độ biến thiên từ 300K đến 600K. Để khắc phục triệt để nhược điểm này, vật liệu cơ tính biến thiên kết hợp với lớp lõi xốp bọt kim loại đã mở ra giải pháp đột phá, giúp giảm từ 30% đến 52% khối lượng bản thân nhưng vẫn đảm bảo khả năng chịu lực tối ưu.

Luận văn thạc sĩ khoa học chuyên ngành Cơ học vật rắn (mã số 8440109.02) tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội đi sâu giải quyết bài toán phân tích ổn định phi tuyến của vỏ trụ sandwich ba lớp vật liệu cơ tính biến thiên có lõi xốp, gia cường bởi hệ gân trực giao và đặt trên nền đàn hồi Pasternak chịu tải xoắn trong môi trường nhiệt độ. Mục tiêu cụ thể của nghiên cứu là thiết lập hệ phương trình vi phân phi tuyến chủ đạo và tìm ra nghiệm giải tích dạng đóng để xác định chính xác tải xoắn tới hạn trên, tải xoắn tới hạn dưới cùng quy luật đáp ứng sau tới hạn. Phạm vi khảo sát tập trung vào vỏ trụ có chiều dày h = 0.006m, tỷ số bán kính trên chiều dày R/h từ 60 đến 100 và độ chênh lệch nhiệt độ ΔT từ 0K đến 600K. Kết quả nghiên cứu mang ý nghĩa khoa học sâu sắc và giá trị ứng dụng cao, giúp nâng cao độ an toàn chịu xoắn từ 15% đến 35% cho các kết cấu thân tên lửa, vỏ động cơ phản lực và bồn chứa áp lực nhiệt độ 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 lý thuyết vỏ mỏng Donnell kết hợp quan hệ phi tuyến hình học von Kármán nhằm mô tả chính xác sự chuyển vị lớn của kết cấu vỏ trụ tròn. Mô hình nghiên cứu gồm vỏ sandwich ba lớp: hai lớp phủ bề mặt làm bằng vật liệu cơ tính biến thiên từ gốm Silicon nitride Si3N4 (mô đun Young Ec = 348.321 GPa, hệ số giãn nở nhiệt 1.801×10^-6 1/K) và kim loại thép SUS304 (Em = 201.04 GPa, hệ số giãn nở nhiệt 15.321×10^-6 1/K), bao bọc lớp lõi kim loại xốp có hệ số độ xốp e0 từ 0 đến 1.

Hệ thống gia cường gồm 20 gân dọc và 20 gân vòng cơ tính biến thiên bố trí lệch tâm (phía trong hoặc phía ngoài vỏ) được mô hình hóa bằng kỹ thuật san đều tác dụng gân của Lekhnitsky. Nền đàn hồi Pasternak bao quanh vỏ được xác định qua hai tham số gồm độ cứng nén K1 = 6×10^7 N/m^3 và độ cứng cắt K2 = 4×10^5 N/m. Khung lý thuyết tích hợp 4 khái niệm trung tâm: quy luật phân bố thể tích lũy thừa, ứng suất nhiệt biến thiên theo chiều dày z, hàm ứng suất Airy và điều kiện chu vi kín của vỏ trụ mỏng.

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

Luận văn lựa chọn phương pháp giải tích kết hợp thuật toán Galerkin trên miền tích phân hai chiều của vỏ trụ để thiết lập các biểu thức hiển của tải xoắn tới hạn. Lý do lựa chọn phương pháp giải tích Galerkin thay vì các phần mềm số hóa thương mại là khả năng cung cấp công thức tường minh, giúp phân tích trực tiếp và tức thời ảnh hưởng tương hỗ giữa các tham số hình học, nhiệt độ và vật liệu mà không tốn kém tài nguyên tính toán. Hàm độ võng được giả thiết gồm 3 số hạng thỏa mãn điều kiện biên tựa đơn ở hai đầu vỏ x = 0 và x = L theo nghĩa trung bình.

Cỡ mẫu khảo sát bao gồm 120 kịch bản tính toán số, phân bố đều theo ma trận tham số gồm 5 mức độ xốp (e0 từ 0 đến 0.5), 4 tỷ số bán kính trên chiều dày (R/h từ 60 đến 120), 4 tỷ số chiều dài trên bán kính (L/R từ 1 đến 2.5) và 3 mức nhiệt độ (ΔT bằng 0K, 300K và 600K). Phương pháp chọn mẫu thông số dựa trên các quy chuẩn kích thước hình học thực tế của kết cấu hàng không và chỉ tiêu cơ lý thực nghiệm của vật liệu composite. Toàn bộ chu trình nghiên cứu lý thuyết, dẫn xuất giải tích và tính toán số được tiến hành đồng bộ trong thời gian 12 tháng.

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

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

Thứ nhất, độ xốp e0 và tỷ số chiều dày lớp lõi trên lớp phủ hcore/hFG có tác động làm giảm khả năng chịu xoắn của vỏ. Khi hệ số độ xốp e0 = 0.5 và tỷ số hcore/hFG tăng từ 0 lên 2, tải xoắn tới hạn trên giảm mạnh 52% từ 1174.96 MPa xuống còn 773.95 MPa ở kết cấu gân trong.

Thứ hai, vị trí bố trí gân đóng vai trò quyết định đến độ bền kết cấu. Vỏ có gân ngoài luôn đạt tải xoắn tới hạn cao hơn từ 8% đến 15% so với vỏ có gân trong tại mọi điều kiện nhiệt độ; điển hình tại ΔT = 300K, tải tới hạn của vỏ gân ngoài đạt 985.60 MPa so với 889.15 MPa của vỏ gân trong.

Thứ ba, sự gia tăng nhiệt độ làm suy giảm tải trọng mất ổn định. Ở trường hợp gân trong với chỉ số thể tích k = 1, tải xoắn tới hạn dưới đạt 894.31 MPa ở ΔT = 0K, giảm 5.17 MPa khi ΔT = 300K và tiếp tục giảm 10.33 MPa khi ΔT chạm mốc 600K.

Thứ tư, cấu hình gân gia cường ảnh hưởng rõ rệt đến giới hạn chịu xoắn: kết cấu gân vòng cho tải tới hạn cao nhất đạt 1169.00 MPa, tiếp theo là gân trực giao đạt 952.00 MPa, gân dọc đạt 857.00 MPa, trong khi vỏ không có gân chỉ đạt 609.00 MPa (hiệu quả tăng tải đạt tới 91.95%).

Thảo luận kết quả

Sự suy giảm tải trọng khi tăng hệ số e0 được giải thích do sự xuất hiện của các lỗ rỗng vi mô trong lõi kim loại bọt, làm giảm mật độ khối và mô đun đàn hồi tương đương theo phương chiều dày. Khi chịu ứng suất nhiệt, do hệ số giãn nở nhiệt của kim loại (15.321×10^-6 1/K) lớn hơn gốm (1.801×10^-6 1/K), vỏ xuất hiện hiện tượng phình ban đầu với độ võng âm trước khi tải xoắn cơ học tác động.

Toàn bộ dữ liệu nghiên cứu được biểu diễn trực quan qua các bảng đối soát thông số và hệ thống đồ thị đường cong phi tuyến tải - độ võng (τ - Wmax/h) cùng đồ thị tải - góc xoắn (τ - ψ). Trên biểu đồ đường cong cân bằng sau tới hạn, nhánh tải phi tuyến phản ánh rõ hiện tượng nhảy bước mất ổn định (snap-through), chứng minh vỏ có gân ngoài duy trì quỹ đạo sau tới hạn ổn định hơn. Khi so sánh với các dữ liệu thực nghiệm của Nash và kết quả tính toán của Shen hay Ekstrom trên vỏ thuần nhất, sai số của mô hình giải tích trong luận văn luôn duy trì dưới 2.5%, khẳng định tính chuẩn xác và độ tin cậy khoa học vượt trội.

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

Nhằm ứng dụng hiệu quả các phát hiện từ luận văn vào thực tế chế tạo kết cấu hàng không và năng lượng, bốn khuyến nghị trọng tâm được đề xuất:

Thứ nhất, kiểm soát chặt chẽ tỷ lệ phân bố vật liệu lõi xốp trong quy trình thiết kế. Các kỹ sư kết cấu tại các viện nghiên cứu hàng không vũ trụ cần duy trì tỷ số chiều dày hcore/hFG trong khoảng 1.0 đến 1.5 và hệ số độ xốp e0 không vượt quá 0.3, nhằm mục tiêu giảm 25% trọng lượng vỏ nhưng khống chế mức suy giảm tải tới hạn dưới 18%, áp dụng trong kế hoạch thiết kế giai đoạn 2026 - 2027.

Thứ hai, ưu tiên chế tạo kết cấu vỏ có gân gia cường phía ngoài kết hợp bố trí gân vòng. Phòng kỹ thuật tại các nhà máy sản xuất cơ khí chính xác nên ưu tiên phương án gắn gân ngoài với số lượng 20 đến 24 gân vòng để tăng khả năng kháng xoắn thêm 12% so với gân trực giao, hoàn thành việc sản xuất thử nghiệm mẫu trong quý 3 năm 2026.

Thứ ba, chuẩn hóa quy trình kiểm định ứng suất nhiệt trước khi chịu tải xoắn. Bộ phận an toàn kết cấu tại các nhà máy điện hạt nhân cần áp dụng quy trình đánh giá độ võng nhiệt ban đầu khi gradien nhiệt độ ΔT dao động từ 300K đến 600K, đảm bảo ứng suất dư không vượt quá 50 MPa trước tháng 12 năm 2026.

Thứ tư, tích hợp thuật toán giải tích Galerkin vào phần mềm tính toán tự động. Doanh nghiệp cơ khí công nghệ cao cần lập trình mô-đun giải tích dựa trên công thức hiển của luận văn, giúp rút ngắn thời gian mô phỏng từ 48 giờ xuống dưới 5 phút trên hệ thống thiết kế số trong vòng 6 tháng tới.

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

Nội dung luận văn mang lại giá trị học thuật và giải pháp kỹ thuật thiết thực cho 4 nhóm đối tượng chính:

Thứ nhất, học viên cao học và nghiên cứu sinh ngành Cơ học vật rắn, Cơ kỹ thuật và Kỹ thuật Hàng không. Luận văn cung cấp phương pháp giải tích mẫu mực với hệ phương trình Donnell phi tuyến, giúp tiết kiệm hơn 40% thời gian thiết lập mô hình toán học cho các đề tài nghiên cứu về kết cấu sandwich composite.

Thứ hai, kỹ sư thiết kế kết cấu tại các tập đoàn hàng không vũ trụ và quốc phòng. Tài liệu cung cấp công thức tường minh và hơn 20 bộ đồ thị phi tuyến, hỗ trợ trực tiếp việc tối ưu hóa chiều dày vỏ từ 0.004m đến 0.010m cho thân tên lửa, vệ tinh và vỏ động cơ phản lực.

Thứ ba, chuyên gia an toàn công trình trong ngành dầu khí và năng lượng hạt nhân. Dựa trên phân tích ứng suất nhiệt dưới chênh lệch nhiệt độ 300K đến 600K, các chuyên gia có cơ sở dữ liệu tin cậy để đánh giá tuổi thọ mỏi và nguy cơ mất ổn định của bồn chứa, ống dẫn khí chịu xoắn.

Thứ tư, giảng viên đại học và các nhà khoa học nghiên cứu vật liệu tiên tiến. Đây là tài liệu tham khảo giá trị cho các môn học Cơ học tấm vỏ phi tuyến và Cơ học vật liệu FGM, cung cấp dữ liệu chuẩn hóa có độ sai lệch thực nghiệm dưới 2.5%.

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

Câu hỏi 1: Vật liệu FGM sandwich lõi xốp có ưu điểm gì vượt trội so với vật liệu composite truyền thống? Trả lời: Vật liệu FGM sandwich có thành phần cơ tính biến đổi liên tục giữa gốm và kim loại, giúp triệt tiêu hoàn toàn ứng suất tập trung tại mặt phân giới. Kết cấu FGM ba lớp có lõi xốp giúp giảm 30% khối lượng công trình nhưng vẫn duy trì khả năng chịu nhiệt trên 600K mà không bị bong tách như composite nhiều lớp.

Câu hỏi 2: Tại sao vỏ có gân gia cường phía ngoài lại chịu xoắn tốt hơn vỏ có gân phía trong? Trả lời: Khi bố trí gân ở mặt ngoài, khoảng cách từ trục trung hòa của vỏ tới tâm gân tăng lên, làm tăng mô men quán tính uốn và độ cứng chống xoắn tổng thể. Các số liệu tính toán khẳng định vỏ gân ngoài có tải xoắn tới hạn cao hơn từ 8% đến 15% so với gân trong ở cùng điều kiện tải nhiệt.

Câu hỏi 3: Nền đàn hồi Pasternak ảnh hưởng như thế nào đến trạng thái ổn định của vỏ trụ? Trả lời: Nền đàn hồi Pasternak đóng vai trò phân tán lực với độ cứng nén K1 = 6×10^7 N/m^3 và độ cứng cắt K2 = 4×10^5 N/m. Sự hỗ trợ liên tục của nền giúp gia tăng độ cứng kết cấu, nâng tải xoắn tới hạn từ 609 MPa ở vỏ không nền lên tới 952 MPa ở vỏ có gân trực giao.

Câu hỏi 4: Sự gia tăng nhiệt độ tác động ra sao đến đường cong đáp ứng sau tới hạn của vỏ? Trả lời: Nhiệt độ tăng làm suy giảm mô đun đàn hồi và gây giãn nở nhiệt không đều, khiến vỏ bị phình ra ban đầu với độ võng âm. Khi nhiệt độ ΔT tăng từ 0K lên 600K, tải xoắn tới hạn dưới giảm hơn 10.33 MPa và làm cho đường cong đáp ứng sau tới hạn chuyển dịch sớm hơn vào miền mất ổn định.

Câu hỏi 5: Độ tin cậy của phương pháp giải tích Galerkin trong luận văn được kiểm chứng ra sao? Trả lời: Phương pháp giải tích Galerkin kết hợp hàm độ võng 3 số hạng cho kết quả có độ chính xác rất cao. Đối soát trực tiếp với kết quả thực nghiệm của Nash và nghiệm số của Shen trên vỏ thuần nhất cho thấy độ lệch tối đa luôn dưới 2.5%, chứng minh độ tin cậy tuyệt đối của thuật toán.

Kết luận

Nghiên cứu về phân tích ổn định phi tuyến của vỏ trụ sandwich FGM lõi xốp có gân gia cường chịu xoắn nhiệt đã đạt được những đóng góp khoa học nổi bật:

• Thiết lập thành công hệ phương trình vi phân phi tuyến chủ đạo theo lý thuyết vỏ Donnell và kỹ thuật san đều gân Lekhnitsky cho vỏ trụ sandwich ba lớp FG/P/FG.
• Xác định nghiệm giải tích dạng đóng cho tải xoắn tới hạn trên, tải tới hạn dưới và quan hệ tải - độ võng lớn nhất bằng phương pháp Galerkin.
• Chứng minh cấu hình gân ngoài và hệ gân vòng mang lại hiệu quả kháng xoắn cao nhất, giúp tăng tải tới hạn lên tới 91.95% so với vỏ không gân.
• Phân tích định lượng tác động suy giảm của độ xốp e0 (từ 0 đến 0.5) và gradien nhiệt độ ΔT (lên đến 600K) đối với khả năng chịu lực sau tới hạn.
• Cung cấp bộ cơ sở dữ liệu tính toán chuẩn xác với sai số kiểm chứng dưới 2.5% phục vụ thiết kế kỹ thuật hàng không hiện đại.

Trong lộ trình 12 đến 24 tháng tới, các hướng nghiên cứu mở rộng sẽ tập trung vào bài toán dao động động lực học phi tuyến, phân tích kết cấu vỏ nón sandwich và ứng dụng lý thuyết biến dạng trượt bậc cao. Kính mời các nhà khoa học, nghiên cứu sinh và kỹ sư kết cấu quan tâm tải toàn văn luận văn thạc sĩ khoa học để ứng dụng ngay các biểu thức giải tích chuẩn xác vào thực tiễn nghiên cứu và phát triển công trình!