Nghiên Cứu Mô Phỏng Đánh Giá Cơ Chế Phá Huỷ Của Cấu Trúc Composite Máy Bay Khi Va Chạm

Tài liệu nghiên cứu Xây dựng mô phỏng đánh giá ơ chế phá hủy của cấu trúc composite của máy bay khi va chạm, tổng hợp lý thuyết và thực hành, cung cấp kiến thức chuyên sâu về kỹ

• Tác giả: Lã Tiến Thành
• Người hướng dẫn: TS. Lê Thị Tuyết Nhung
• Trường học: Trường Đại Học Bách Khoa Hà Nội
• Chuyên ngành: Kỹ Thuật Cơ Khí Động Lực
• Đề tài: Nghiên Cứu Mô Phỏng Va Chạm Của Cấu Trúc Composite Máy Bay
• 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

• LỜI MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Tai nạn trong hàng không
• 1.2. Chứng nhận bird strikes
• 1.3. Va chạm của chim
• 1.3.1. Va chạm ban đầu
• 1.3.2. Giải phóng áp suất
• 1.3.3. Dòng ổn định
• 1.3.4. Kết thúc dòng
• 1.4. Chuyển hóa động lượng
• 1.5. Thời gian va chạm
• 1.6. Lực va chạm trung bình
• 2. CHƯƠNG II: LÝ THUYẾT
• 2.1. Phương pháp mô phỏng
• 2.1.1. Phương pháp Lagrangian
• 2.1.2. Phương pháp Eulerian
• 2.1.3. Phương pháp Coupled Eulerian Lagrangian (CEL)
• 2.1.4. Phương pháp Smooth Particle Hydrodynamics (SPH)
• 2.2. Phương trình điều khiển
• 2.2.1. Bảo toàn khối lượng
• 2.2.2. Bảo toàn động lượng tuyến tính
• 2.2.3. Bảo toàn momem động lượng
• 2.2.4. Bảo toàn năng lượng
• 2.2.5. Nguyên tắc Vitual Work
• 2.2.6. Phương pháp CEL
• 2.2.7. Phương pháp Smoothed Particle Hydrodynamics (SPH)
• 3. CHƯƠNG III: MÔ PHỎNG
• 3.1. Mô hình mép vào cánh
• 3.1.1. Mô hình vật liệu composite và tiêu chuẩn phá hủy
• 3.1.2. Mô hình vật liệu lõi
• 3.1.3. Mô hình chim
• 3.1.4. Xác nhận mô hình
• 3.1.4.1. Cấu trúc sandwich
• 3.1.4.2. Phương pháp mô phỏng
• 3.1.4.2.1. Cánh sử dụng lõi ACG-1/4
• 3.1.4.2.2. Cánh sử dụng lõi ACG-1
• 3.1.4.3. Cân bằng năng lượng
• 4. CHƯƠNG IV: KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Nghiên Cứu Mô Phỏng Va Chạm Cấu Trúc Composite

Nghiên cứu mô phỏng va chạm cấu trúc composite máy bay ngày càng trở nên quan trọng do sự gia tăng sử dụng vật liệu composite trong ngành hàng không. Các phần mềm như Abaqus/Explicit đóng vai trò then chốt trong việc đánh giá thiết kế, giảm chi phí và thời gian thử nghiệm. Yêu cầu về khả năng chịu chim đâm (birdstrikes) là một trong những yếu tố được quan tâm hàng đầu. Nghiên cứu này tập trung vào việc xây dựng mô hình số có khả năng đánh giá cơ chế phá hủy của vật liệu composite cấu trúc sandwich khi chịu va chạm ở vận tốc cao, cụ thể là mép vào cánh. Các tai nạn hàng không có thể xảy ra do va chạm với các vật thể bên ngoài, đặc biệt là chim. Việc va chạm với chim không chỉ gây nguy hiểm cho an toàn bay mà còn gây ra chi phí bảo trì đáng kể.

1.1. Tổng Quan Về Tai Nạn Hàng Không Do Va Chạm

Trong suốt vòng đời của máy bay, luôn tồn tại nguy cơ va chạm với các đối tượng bên ngoài. Theo tiêu chuẩn hàng không, thuật ngữ "birdstrike" dùng để chỉ va chạm giữa chim và các bộ phận của máy bay như kính chắn gió, thân, cánh và động cơ. Nguy cơ tai nạn cao nhất ở khu vực sân bay trong giai đoạn cất cánh và hạ cánh, đặc biệt vào sáng sớm và chiều muộn. Mức độ nghiêm trọng của va chạm chim đã tăng lên do sự gia tăng vận tải hàng không và tốc độ máy bay. Sự cố birdstrike không chỉ ảnh hưởng đến an toàn bay mà còn phát sinh chi phí bảo trì đáng kể. Các tổ chức quốc tế như the Birdstrike Committee USA đã được thành lập để nghiên cứu và giảm thiểu rủi ro này. Theo thống kê, hơn 60,000 vụ va chạm giữa động vật hoang dã và máy bay dân dụng đã được ghi nhận tại Mỹ từ năm 1990 đến 2005. Thiệt hại hàng năm ước tính vượt quá 530,000 giờ bay và 614 triệu đô la.

1.2. Các Tiêu Chuẩn Chứng Nhận Khả Năng Chống Va Chạm

Việc kiểm tra khả năng chịu va chạm với chim được thực hiện theo các quy định hàng không Liên bang (FAR) Phần 23, 25 và 33. Trước đây, thử nghiệm sử dụng gà sống bắn vào cấu trúc. Tuy nhiên, hiện nay, mô hình chim tổng hợp với kích thước và trọng lượng phù hợp được sử dụng. Một chương trình kiểm tra kính chắn gió bao gồm nhiều lần bắn chim vào các vị trí khác nhau trên kính và khung. Để được chứng nhận, kính chắn gió và khung phải chứng minh rằng phi công không bị thương, các mảnh vỡ không gây nguy hiểm và cấu trúc vẫn giữ được áp suất cabin sau va chạm. Các camera tốc độ cao được sử dụng để ghi lại quá trình phá hủy. Đối với cánh máy bay, các kiểm tra được thực hiện ở mép vào, đảm bảo chim không xuyên qua hoặc gây hư hỏng nghiêm trọng cho xà dọc trước của cánh. Tương tự, các cánh đuôi cũng được kiểm tra để đảm bảo an toàn. Đối với động cơ, động cơ phải tiếp tục hoạt động với 75% lực đẩy trong 5 phút sau va chạm và không bị bắt lửa hoặc vỡ sau va chạm.

II. Thách Thức Nghiên Cứu Mô Phỏng Va Chạm Cấu Trúc Máy Bay

Nghiên cứu mô phỏng va chạm cấu trúc composite máy bay đối mặt với nhiều thách thức. Mô hình hóa chính xác vật liệu composite và cơ chế phá hủy của chúng là một vấn đề phức tạp. Các phương pháp mô phỏng khác nhau, như Lagrangian, Eulerian, CEL và SPH, có những ưu điểm và nhược điểm riêng, đòi hỏi sự lựa chọn phù hợp cho từng bài toán cụ thể. Việc xác định các thông số vật liệu chính xác và điều kiện biên phù hợp cũng rất quan trọng để đảm bảo tính chính xác của kết quả mô phỏng. Ngoài ra, mô phỏng va chạm ở vận tốc cao đòi hỏi khả năng tính toán lớn và các thuật toán hiệu quả.

2.1. Khó Khăn Trong Mô Hình Hóa Vật Liệu Composite

Vật liệu composite có cấu trúc phức tạp và tính chất không đồng nhất, gây khó khăn trong việc mô hình hóa. Các cơ chế phá hủy của composite, như phân lớp, nứt ma trận và đứt sợi, cần được mô tả chính xác trong mô hình. Việc xác định các thông số vật liệu cần thiết cho mô hình, như độ bền, mô đun đàn hồi và các tham số phá hủy, cũng là một thách thức. Hơn nữa, các mô hình phá hủy khác nhau có thể cho kết quả khác nhau, đòi hỏi sự lựa chọn và hiệu chỉnh cẩn thận.

2.2. Lựa Chọn Phương Pháp Mô Phỏng Va Chạm Phù Hợp

Việc lựa chọn phương pháp mô phỏng va chạm phù hợp là rất quan trọng. Phương pháp Lagrangian phù hợp cho các bài toán với biến dạng nhỏ và không có phá hủy lớn. Phương pháp Eulerian phù hợp cho các bài toán với biến dạng lớn và dòng chảy vật liệu. Phương pháp CEL kết hợp ưu điểm của cả hai phương pháp Lagrangian và Eulerian, phù hợp cho các bài toán tương tác giữa cấu trúc và chất lỏng. Phương pháp SPH là một phương pháp không lưới, phù hợp cho các bài toán với biến dạng cực lớn và phá hủy mạnh. Việc lựa chọn phương pháp phù hợp phụ thuộc vào đặc điểm của bài toán và khả năng tính toán.

III. Phương Pháp Mô Phỏng Va Chạm Bằng CEL và SPH Trong Abaqus

Luận văn so sánh hai phương pháp mô phỏng phổ biến là CEL (Coupled Eulerian Lagrangian) và SPH (Smooth Particle Hydrodynamics) trong phần mềm Abaqus/Explicit. CEL kết hợp khả năng mô tả cấu trúc chính xác của Lagrangian với khả năng xử lý biến dạng lớn của Eulerian. SPH là phương pháp không lưới, thích hợp cho các bài toán va chạm với biến dạng cực lớn và phá hủy phức tạp. Việc so sánh hai phương pháp này giúp xác định phương pháp phù hợp nhất cho việc mô phỏng va chạm cấu trúc composite sandwich.

3.1. Ưu Điểm và Nhược Điểm của Phương Pháp CEL

Phương pháp CEL cho phép mô tả chính xác sự tương tác giữa vật thể Lagrangian (cấu trúc composite) và vật thể Eulerian (chim). Nó có khả năng xử lý biến dạng lớn của chim mà vẫn duy trì độ chính xác của cấu trúc composite. Tuy nhiên, CEL có thể tốn nhiều thời gian tính toán và yêu cầu kỹ năng thiết lập mô hình tốt. Độ chính xác của CEL phụ thuộc vào kích thước phần tử và sự phù hợp của lưới Eulerian.

3.2. Ưu Điểm và Nhược Điểm của Phương Pháp SPH

Phương pháp SPH có khả năng xử lý biến dạng cực lớn và phá hủy phức tạp một cách tự nhiên. Nó không yêu cầu lưới và dễ dàng thiết lập mô hình. Tuy nhiên, SPH có thể kém chính xác hơn CEL trong việc mô tả cấu trúc và sự tương tác giữa các vật thể. Kết quả của SPH có thể nhạy cảm với các tham số mô hình và kích thước hạt.

3.3. Xác Thực Mô Hình Bằng Thí Nghiệm và Mô Phỏng Chuẩn

Để đảm bảo tính chính xác của mô hình, cần phải xác thực mô hình bằng cách so sánh kết quả mô phỏng với kết quả thí nghiệm hoặc các mô phỏng chuẩn khác. Trong luận văn, mô hình chim được xác thực bằng cách mô phỏng bắn chim vào tấm nhôm phẳng và so sánh kết quả với các nghiên cứu trước đó. Mô hình cấu trúc sandwich được xác thực bằng cách mô phỏng bắn quả cầu thép vào tấm sandwich và so sánh kết quả với thí nghiệm.

IV. Ứng Dụng Mô Phỏng Va Chạm Chim Vào Mép Cánh Composite

Sau khi xác thực mô hình, luận văn ứng dụng phương pháp mô phỏng để nghiên cứu va chạm chim vào mép cánh composite. Các mô phỏng được thực hiện với hai loại lõi khác nhau để so sánh ảnh hưởng của lõi đến khả năng chịu va chạm. Kết quả mô phỏng cho thấy cơ chế phá hủy của cấu trúc sandwich và khả năng hấp thụ năng lượng của nó khi va chạm với chim. Mô phỏng này là quan trọng để đánh giá thiết kế của cánh máy bay và đảm bảo an toàn.

4.1. Mô Hình Hóa Mép Cánh Composite và Lõi Sandwich

Mép cánh composite được mô hình hóa bằng vật liệu composite lớp lớp với các tham số vật liệu phù hợp. Lõi sandwich được mô hình hóa bằng vật liệu tổ ong với các tính chất đàn hồi và dẻo. Mô hình phải bao gồm các chi tiết quan trọng như kích thước, hình dạng và cách bố trí các lớp composite và lõi.

4.2. So Sánh Kết Quả Mô Phỏng Với Các Loại Lõi Khác Nhau

Kết quả mô phỏng được so sánh với hai loại lõi khác nhau để đánh giá ảnh hưởng của lõi đến khả năng chịu va chạm. Các thông số so sánh bao gồm độ võng, ứng suất, biến dạng và năng lượng hấp thụ. So sánh này giúp xác định loại lõi nào có khả năng bảo vệ mép cánh tốt hơn.

V. Phân Tích Kết Quả Mô Phỏng và Cơ Chế Phá Hủy Composite

Kết quả mô phỏng va chạm cung cấp thông tin chi tiết về cơ chế phá hủy của cấu trúc composite sandwich. Các cơ chế phá hủy chính bao gồm nứt ma trận, phân lớp, đứt sợi và phá hủy lõi. Phân tích này giúp hiểu rõ cách cấu trúc hấp thụ năng lượng va chạm và cách nó bị phá hủy. Thông tin này có thể được sử dụng để cải thiện thiết kế và tăng cường khả năng chịu va chạm.

5.1. Xác Định Các Cơ Chế Phá Hủy Chính Trong Quá Trình Va Chạm

Quá trình va chạm dẫn đến các cơ chế phá hủy khác nhau trong vật liệu composite. Nứt ma trận xảy ra khi ứng suất vượt quá độ bền của ma trận. Phân lớp xảy ra khi các lớp composite tách rời nhau. Đứt sợi xảy ra khi ứng suất vượt quá độ bền của sợi. Phá hủy lõi xảy ra khi lõi sandwich bị nén hoặc cắt. Xác định các cơ chế phá hủy chính giúp hiểu rõ cách cấu trúc bị phá hủy.

5.2. Đánh Giá Khả Năng Hấp Thụ Năng Lượng Của Cấu Trúc

Khả năng hấp thụ năng lượng là một thước đo quan trọng về khả năng chịu va chạm của cấu trúc. Khả năng hấp thụ năng lượng của cấu trúc phụ thuộc vào vật liệu, hình dạng và cơ chế phá hủy. Đánh giá khả năng hấp thụ năng lượng giúp xác định cách cải thiện thiết kế để tăng cường khả năng chịu va chạm.

VI. Kết Luận và Hướng Nghiên Cứu Tương Lai Về Mô Phỏng Va Chạm

Nghiên cứu mô phỏng va chạm là một công cụ quan trọng để đánh giá thiết kế và tăng cường khả năng chịu va chạm của cấu trúc composite máy bay. Luận văn đã trình bày một phương pháp mô phỏng hiệu quả bằng cách sử dụng CEL và SPH trong Abaqus/Explicit. Kết quả mô phỏng cung cấp thông tin chi tiết về cơ chế phá hủy và khả năng hấp thụ năng lượng của cấu trúc sandwich. Trong tương lai, nghiên cứu có thể được mở rộng để xem xét các yếu tố khác như ảnh hưởng của nhiệt độ, độ ẩm và các điều kiện môi trường khác.

6.1. Tóm Tắt Kết Quả và Đóng Góp Của Nghiên Cứu

Nghiên cứu đã trình bày một phương pháp mô phỏng va chạm hiệu quả bằng cách sử dụng CEL và SPH trong Abaqus/Explicit. Nghiên cứu đã so sánh hai phương pháp này và xác định phương pháp phù hợp nhất cho việc mô phỏng va chạm cấu trúc composite sandwich. Nghiên cứu đã ứng dụng phương pháp mô phỏng để nghiên cứu va chạm chim vào mép cánh composite và phân tích cơ chế phá hủy. Nghiên cứu cung cấp thông tin hữu ích cho việc thiết kế và cải thiện khả năng chịu va chạm của cấu trúc composite máy bay.

6.2. Hướng Nghiên Cứu Phát Triển Trong Tương Lai

Trong tương lai, nghiên cứu có thể được mở rộng để xem xét các yếu tố khác như ảnh hưởng của nhiệt độ, độ ẩm và các điều kiện môi trường khác. Nghiên cứu cũng có thể được mở rộng để xem xét các loại vật liệu composite và cấu trúc khác. Ngoài ra, có thể phát triển các mô hình phá hủy tiên tiến hơn để mô tả chính xác hơn cơ chế phá hủy của vật liệu composite.

Nội dung chính

Tổng quan nghiên cứu

Trong ngành hàng không, va chạm giữa chim và máy bay (bird strike) là một trong những nguyên nhân chính gây ra tai nạn và thiệt hại nghiêm trọng cho cấu trúc tàu bay. Theo báo cáo của ngành, từ năm 1990 đến 2005, hơn 60.000 vụ va chạm với động vật hoang dã đã được ghi nhận tại Mỹ, gây thiệt hại ước tính lên đến 614 triệu đô la Mỹ và làm lãng phí hơn 530.000 giờ bay. Vấn đề này càng trở nên nghiêm trọng khi tốc độ và mật độ vận tải hàng không tăng cao, đặc biệt trong giai đoạn cất cánh và hạ cánh tại sân bay. Mục tiêu nghiên cứu là xây dựng mô phỏng đánh giá cơ chế phá hủy của cấu trúc composite sandwich ở mép vào cánh máy bay khi chịu va chạm chim với vận tốc cao, nhằm hỗ trợ thiết kế và chứng nhận an toàn cho các bộ phận tàu bay. Nghiên cứu tập trung vào mô phỏng va chạm chim với vận tốc khoảng 150 m/s, sử dụng phần mềm Abaqus/Explicit, so sánh hai phương pháp mô phỏng Coupled Eulerian Lagrangian (CEL) và Smooth Particle Hydrodynamics (SPH), đồng thời phân tích ảnh hưởng của hai loại lõi tổ ong khác nhau trong cấu trúc sandwich. Phạm vi nghiên cứu được giới hạn trong mô hình mép vào cánh với vật liệu composite carbon và lõi tổ ong nhôm, áp dụng cho điều kiện va chạm thực tế tại các sân bay. Kết quả nghiên cứu có ý nghĩa quan trọng trong việc giảm chi phí thử nghiệm thực tế, nâng cao độ tin cậy và an toàn của tàu bay trước các sự cố va chạm chim.

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 dựa trên lý thuyết va chạm cơ thể mềm (soft body impact) và thủy động lực học, trong đó chim được xem như một vật thể có tính chất chất lỏng do vận tốc va chạm cao vượt quá giới hạn đàn hồi của vật liệu. Quá trình va chạm được phân thành bốn pha: va chạm ban đầu, giải phóng áp suất, dòng ổn định và kết thúc dòng. Áp suất va chạm được mô tả bằng phương trình sốc Hugoniot và phương trình trạng thái Mie-Gruneisen, cho phép tính toán áp suất và nội năng trong quá trình va chạm. Về vật liệu, cấu trúc composite sandwich được mô hình hóa dựa trên tiêu chuẩn phá hủy Hashin, đánh giá các chế độ phá hủy sợi và nền trong vật liệu CFRP. Lõi tổ ong nhôm được mô phỏng với mô hình CRUSHABLE FOAM, phản ánh ứng xử nén và nghiền của lõi. Hai phương pháp mô phỏng chính được áp dụng là Coupled Eulerian Lagrangian (CEL) và Smooth Particle Hydrodynamics (SPH). Phương pháp CEL kết hợp ưu điểm của lưới Lagrangian (cho vật rắn) và Eulerian (cho chất lỏng), giúp mô phỏng tương tác phức tạp giữa chim và cấu trúc. Phương pháp SPH là phương pháp không lưới, mô phỏng vật liệu dưới dạng các hạt tương tác, phù hợp với các bài toán biến dạng lớn và va chạm phức tạp.

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

Nguồn dữ liệu nghiên cứu bao gồm số liệu thực nghiệm về va chạm chim với tấm nhôm và tấm composite sandwich, cùng các thông số cơ tính vật liệu CFRP AS4/8552, lõi tổ ong nhôm 3003, và nhôm Al 6061-T6. Cỡ mẫu mô phỏng gồm hàng trăm nghìn phần tử: 426.120 phần tử Eulerian cho chim, 10.750 phần tử rắn cho lõi và xà dọc, 35.780 phần tử vỏ cho lớp composite. Phương pháp chọn mẫu là mô phỏng số học phần tử hữu hạn với phần mềm Abaqus/Explicit, sử dụng hai kỹ thuật mô phỏng CEL và SPH để so sánh hiệu quả và độ chính xác. Timeline nghiên cứu kéo dài trong năm 2018, bao gồm các bước: xây dựng mô hình chim và cấu trúc, xác nhận mô hình bằng mô phỏng va chạm quả cầu thép và tấm nhôm, thực hiện mô phỏng va chạm chim với mép vào cánh, phân tích kết quả và so sánh hai loại lõi tổ ong. Phân tích dữ liệu tập trung vào các biến chuyển vị, ứng suất, thiệt hại vật liệu và cân bằng năng lượng trong quá trình va chạm.

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

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

1. Xác nhận mô hình chim và vật liệu: Mô phỏng va chạm chim với tấm nhôm ở vận tốc 150 m/s cho thấy phương pháp CEL cho kết quả chuyển vị lớn nhất tại tâm tấm là 39,4 mm, gần với kết quả thực nghiệm 41 mm, trong khi phương pháp SPH cho kết quả 46,1 mm, thể hiện độ chính xác cao của phương pháp CEL.
2. So sánh phương pháp mô phỏng: Phương pháp CEL có ưu điểm về độ chính xác và mô phỏng dòng chảy thực tế hơn, tuy nhiên thời gian tính toán lâu hơn so với SPH. Phương pháp SPH có thời gian tính toán nhanh hơn nhưng đòi hỏi thời gian tạo lưới lâu và chia lưới phức tạp.
3. Ảnh hưởng của lõi tổ ong: Hai loại lõi tổ ong ACG-1/4 và ACG-1 được so sánh trong mô phỏng va chạm chim với mép vào cánh. Kết quả cho thấy lõi ACG-1 có khả năng hấp thụ năng lượng va chạm cao hơn, thể hiện qua các đường cong ứng suất-biến dạng và cân bằng năng lượng, giúp giảm thiểu thiệt hại cấu trúc.
4. Cơ chế phá hủy vật liệu composite: Tiêu chuẩn Hashin được áp dụng để mô phỏng các chế độ phá hủy sợi và nền, cho thấy sự phát triển thiệt hại theo chuyển vị tương đương, với các phần tử bị xóa khi đạt ngưỡng phá hủy, phản ánh chính xác quá trình hư hỏng vật liệu trong va chạm.

Thảo luận kết quả

Nguyên nhân kết quả mô phỏng CEL gần với thực nghiệm hơn là do phương pháp này kết hợp ưu điểm của lưới Eulerian và Lagrangian, giúp mô phỏng chính xác tương tác giữa chim (vật thể mềm) và cấu trúc cứng. Kết quả này phù hợp với các nghiên cứu trước đây trong lĩnh vực mô phỏng va chạm cơ thể mềm. Việc lõi tổ ong ACG-1 hấp thụ năng lượng tốt hơn được giải thích bởi kích thước ô tổ ong lớn hơn, tạo ra diện tích liên kết lớn và khả năng chịu nén tốt hơn, đồng thời giảm thiểu biến dạng cục bộ. Kết quả này có ý nghĩa quan trọng trong thiết kế cấu trúc composite chịu va chạm, giúp lựa chọn vật liệu lõi phù hợp để nâng cao độ bền và an toàn. Dữ liệu có thể được trình bày qua biểu đồ ứng suất-biến dạng, đồ thị cân bằng năng lượng và hình ảnh mô phỏng chuyển vị để minh họa rõ ràng cơ chế phá hủy và hiệu quả hấp thụ năng lượng của các cấu trúc khác nhau.

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

1. Áp dụng phương pháp mô phỏng CEL trong thiết kế: Khuyến nghị các nhà thiết kế và kỹ sư hàng không sử dụng phương pháp CEL để mô phỏng va chạm chim với cấu trúc composite nhằm nâng cao độ chính xác và tin cậy của kết quả, đặc biệt trong giai đoạn thiết kế và chứng nhận sản phẩm.
2. Lựa chọn lõi tổ ong phù hợp: Đề xuất ưu tiên sử dụng lõi tổ ong có kích thước ô lớn hơn như ACG-1 để tăng khả năng hấp thụ năng lượng va chạm, giảm thiểu thiệt hại cấu trúc, đồng thời cân nhắc chi phí và trọng lượng trong thiết kế.
3. Tăng cường thử nghiệm mô phỏng kết hợp thực nghiệm: Khuyến khích thực hiện các thử nghiệm thực tế kết hợp với mô phỏng số để xác nhận và hiệu chỉnh mô hình, đảm bảo tính ứng dụng cao trong thực tế sản xuất và vận hành.
4. Phát triển phần mềm mô phỏng chuyên biệt: Đề xuất nghiên cứu và phát triển các công cụ mô phỏng chuyên biệt cho va chạm chim trong hàng không, tích hợp các mô hình vật liệu tiên tiến và thuật toán tối ưu để giảm thời gian tính toán và nâng cao hiệu quả phân tích.

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

1. Kỹ sư thiết kế hàng không: Giúp hiểu rõ cơ chế phá hủy của cấu trúc composite khi va chạm chim, từ đó tối ưu thiết kế mép vào cánh và các bộ phận chịu lực khác.
2. Nhà nghiên cứu vật liệu composite: Cung cấp dữ liệu và mô hình vật liệu composite sandwich và lõi tổ ong, hỗ trợ phát triển vật liệu mới có khả năng chịu va chạm tốt hơn.
3. Chuyên gia mô phỏng và phân tích kết cấu: Tham khảo phương pháp mô phỏng CEL và SPH, áp dụng trong các bài toán tương tác chất lỏng - cấu trúc phức tạp.
4. Cơ quan chứng nhận hàng không: Sử dụng kết quả nghiên cứu để xây dựng tiêu chuẩn và quy trình kiểm tra va chạm chim, đảm bảo an toàn bay và giảm thiểu chi phí bảo trì.

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

1. Phương pháp CEL và SPH khác nhau như thế nào trong mô phỏng va chạm chim?
   Phương pháp CEL kết hợp lưới Eulerian và Lagrangian, mô phỏng chính xác tương tác giữa vật rắn và chất lỏng, cho kết quả gần thực nghiệm hơn nhưng tốn thời gian tính toán. SPH là phương pháp không lưới, mô phỏng vật liệu dưới dạng các hạt, nhanh hơn nhưng phức tạp trong tạo lưới và có thể kém chính xác hơn.

2. Tại sao lại chọn mô hình vật liệu composite dựa trên tiêu chuẩn Hashin?
   Tiêu chuẩn Hashin cho phép đánh giá chính xác các chế độ phá hủy sợi và nền trong composite, phù hợp với đặc tính vật liệu CFRP và giúp mô phỏng quá trình hư hỏng vật liệu dưới tải va chạm.

3. Lõi tổ ong ảnh hưởng thế nào đến khả năng chịu va chạm của cấu trúc?
   Lõi tổ ong có kích thước ô lớn hơn tạo ra diện tích liên kết lớn, khả năng chịu nén và hấp thụ năng lượng tốt hơn, giúp giảm thiểu thiệt hại khi va chạm chim.

4. Mô hình chim được xây dựng như thế nào trong nghiên cứu?
   Chim được mô phỏng đơn giản bằng hình trụ có đầu bán cầu, vật liệu thay thế bằng khối lượng nước tương ứng, sử dụng phương trình trạng thái Mie-Gruneisen để mô phỏng phản ứng thủy động lực.

5. Kết quả mô phỏng có thể ứng dụng thực tế như thế nào?
   Kết quả giúp các nhà thiết kế và cơ quan chứng nhận đánh giá khả năng chịu va chạm của cấu trúc composite, từ đó tối ưu thiết kế, giảm chi phí thử nghiệm thực tế và nâng cao an toàn bay.

Kết luận

• Nghiên cứu đã xây dựng thành công mô hình mô phỏng va chạm chim với cấu trúc composite sandwich mép vào cánh máy bay sử dụng phần mềm Abaqus/Explicit.
• Phương pháp CEL cho kết quả mô phỏng chính xác và gần với thực nghiệm hơn so với SPH, phù hợp cho các bài toán tương tác chất lỏng - cấu trúc phức tạp.
• Lõi tổ ong ACG-1 thể hiện khả năng hấp thụ năng lượng va chạm tốt hơn so với lõi ACG-1/4, góp phần giảm thiểu thiệt hại cấu trúc.
• Tiêu chuẩn phá hủy Hashin được áp dụng hiệu quả trong mô phỏng quá trình hư hỏng vật liệu composite dưới tác động va chạm.
• Đề xuất tiếp tục phát triển mô hình, mở rộng phạm vi nghiên cứu và ứng dụng trong thiết kế, chứng nhận an toàn hàng không.

Next steps: Triển khai thử nghiệm thực tế kết hợp mô phỏng, phát triển công cụ mô phỏng chuyên biệt, và áp dụng kết quả vào thiết kế sản phẩm mới.

Các nhà nghiên cứu và kỹ sư trong lĩnh vực hàng không nên áp dụng phương pháp mô phỏng CEL và tiêu chuẩn Hashin để nâng cao hiệu quả thiết kế và đảm bảo an toàn bay trong tương lai.

Trích đoạn nội dung tài liệu

Chương I: Tổng quan ^ -(( = ^ ` 'B abZ [ //////////////////////////////@ABcC ' là bán kính mặt cắt vật phóng và [/là góc va chạm Hình 8 Va chạm với góc xiên 1.6 Lực va chạm trung bình Hiệu ứng chính của hiện tượng chim đâm là truyền năng lượng cho cấu trúc tàu bay bị tác động. Nó có thể được tính toán xấp xỉ. Sau va chạm thay đổi động năng của chim xác định bởi phương trình : A 4d# = e = ;' = VT &///////////////////////////////////@ABfC K Với e là công ; là lực tác động ' là khoảng cánh lực thực hiện V là khối lượng chim T là vận tốc tàu bay Lực tác dụng lên tàu bay : 4d# VT & ;= = ////////////////////////////////////@ABgC ^-(( K^-(( Trong trường hợp vật đích là tấm phẳng có thể biến dạng, lực truyền từ chim phụ thuộc nhiều vào hình dạng và độ lớn biến dạng. 11 Chương II: Lý thuyết CHƯƠNG II: LÝ THUYẾT 2.1 Phương pháp mô phỏng Mặc dù đã có nhiều nghiên cứu về tác động của cơ thể mềm nhưng không có phương pháp chuẩn nào có sẵn để phân tích các vấn đề va chạm tương tác cấu trúc chất lỏng.

Hiện tại có bốn phương pháp mô hình hoá, như Lagrangian, Eulerian, Couuler Lagrangian Eulerian (CEL), và Smooth Particle Hydrodynamics, được sử dụng cho các phân tích thiệt hại do va chạm. Mỗi phương pháp đều có điểm mạnh và điểm yếu riêng.1 Phương pháp Lagrangian Phương pháp Lagrangian chia thể tích thành một số lượng lớn các hình khối nhỏ được gọi là phần tử, và nhìn chung nó phù hợp để mô tả vật liệu rắn (solid) chịu tác động, khi đó lưới số học di chuyển và bị méo như hình 9. Hình 9: Mô tả biến đổi của lưới Lagrangian Khác biệt chính với các phương pháp khác là lựa chọn hệ tọa độ tham chiếu để mô tả sự thay đổi. Phương pháp Lagrangian sử dụng hệ tọa độ vật liệu để tham chiếu.

Các nút nằm trong lưới liên kết với các hạt vật liệu dưới một sự giám sát, sau đó mỗi nút trong lưới thay đổi theo mỗi hạt riêng biệt. Tính chất vật liệu theo thời gian cũng được mô tả tốt với cánh tiếp cận này vì khả năng dễ dàng ghi lại vật liệu theo thời gian. Hơn nữa mặt phẳng tự do và bề mặt vật liệu được xác định vị trí tại vùng biên và được giữ vững tốt trong suốt quá trình tính toán, điều này giúp đặt điều kiện biên dễ dàng. Lưới Lagrangian đặc biệt thích hợp cho mô tả ứng xử của vật rắn, đó cũng chính là nhược điểm chính, do tính chất của các dạng biến dạng lớn, lưới số học có 12 Chương II: Lý thuyết thể bị quá méo khiến bước thời gian nhỏ và có thể mất độ chính xác.

Trong trường hợp đó bộ giải số học chỉ có thể được tiến hành cho một điểm chắc chắn trước khi lưới Lagragian méo gây ra phân tích bị dừng lại do bước thời gian rất nhỏ. Đối với phương pháp này vật liệu di chuyển theo lưới, nếu vật liệu chịu biến dạng lớn, lưới cũng sẽ chịu biến dạng ngang bằng và nó khiến kết quả không chính xác và mất ổn định số học ( lưới rối và/hoặc thể tích lưới là âm). Tuy nhiên, chương trình có thể bổ sung thêm các đặc tính, như “rezoning” và “erosion” có thể áp dụng cho lưới Lagrangian để mở rộng phân tích. Với phần tử erosion, phần tử méo bị xóa khỏi mesh bị xóa khỏi lưới bằng cách đặt một giới hạn độ bền dẻo.

Xóa phần tử vượt quá giá trị ngưỡng biến dạng dẻo giải quyết cả vấn đề thể tích âm lẫn bước thời gian bị giảm. Mặt khác, việc xóa phần tử cũng xóa khối lượng và năng lượng biến dạng của cấu trúc, vì vậy không đảm bảo bảo toàn khối lượng và năng lượng. Quy trình rezoning đòi hỏi từng bước chia lưới lại của phần bị méo và nó có thể khiến tăng thời gian bộ giải và lỗi số học liên quan đến các xấp xỉ. Vì lý do đó kỹ thuật này hoạt động tốt nhất thì cần có thêm các bước chia lưới lại.2 Phương pháp Eulerian Trong công thức Eulerian, lưới tham chiếu được coi như thể tích điều khiển, tức là lưới vẫn cố định, và một vật liệu nghiên cứu chảy qua lưới.

Lưới không di chuyển, không có khả năng biến dạng lưới, đó là một lợi thế lớn. Công thức này hầu hết được áp dụng cho các ứng dụng chất lỏng. Tuy nhiên, công thức Eulerian yêu cầu tính toán nhiều hơn Lagrangian, dẫn đến thời gian mô phỏng dài hơn. Ngoài ra, rất khó để theo dõi các giao diện vật liệu, và lịch sử của các biến vật liệu.

Trong Hình 10, mô phỏng tác động cơ thể mềm sử dụng công thức Eulerian được hiển thị ở các thời điểm khác nhau. Hình 10 Phương pháp Eulerian 13 Chương II: Lý thuyết 2.3 Phương pháp Coupled Eulerian Lagrangian (CEL) CEL là sự kết hợp của Lagrangian, và Eulerian. Mục đích chính của phương pháp CEL là sử dụng các lợi thế của lưới Lagrangian và Eulerian. Đối với vấn đề kết cấu - chất lỏng tổng quát, lưới Lagrangian được sử dụng để phân biệt kết cấu, trong khi lưới Eulerian được sử dụng để phân biệt một chất lỏng.

Giao diện giữa kết cấu và chất lỏng có thể được biểu diễn bằng cách sử dụng ranh giới của miền Lagrangian. Mặt khác, lưới Eulerian, đại diện cho chất lỏng có thể bị biến dạng lớn, không có vấn đề liên quan đến lưới, và các biến dạng phần tử. Hạn chế duy nhất của phương pháp CEL là thời gian tính toán dài hơn. Trong hình 11, một mô phỏng tác động cơ thể mềm sử dụng công thức CEL ở những thời điểm khác nhau của thời gian.

Hình 11 Mô hình trong CEL 2.4 Phương pháp Smooth Particle Hydrodynamics (SPH) Phương pháp không lưới SPH là một phương pháp tính toán kết cấu được gây dựng độc lập vào năm 1977 bởi Monaghan Gingold và Lucy, ban đầu để giải quyết các bài toán trong lĩnh vực thiên văn học. Tuy nhiên, ngày nay nó đã được sử dụng rộng rãi trong các lĩnh vực tính toán khác như cơ học chất lỏng, va chạm kết cấu. Phương pháp không lưới SPH chia kết cấu tính toán thành một nhóm các nút phần tử (hạt). Mỗi nút phần tử đại diện bởi vị trí, khối lượng, vận tốc, ứng suất.

Mỗi hạt 14 Chương II: Lý thuyết có một vùng ảnh hưởng của nó lên các hạt khác. Ảnh hưởng này của hạt lên các hạt khác được diễn tả bằng một hàm kernel. Mức độ ảnh hưởng phụ thuộc vào khoảng cách từ hạt đó đến hạt đang xét. Ảnh hưởng là lớn nhất tại trung tâm và bằng không tại biên của vùng ảnh hưởng.

Hình 12 Mô hình trong phương pháp SPH 2.2 Phương trình điều khiển Mô phỏng va chạm của chim với vận tốc cao là một vấn đề phi tuyến với cả động học phi tuyến và vật liệu phi tuyến do tỉ lệ biến dạng cao, thay đổi hình dạng lớn và biến dạng không đàn hồi. Một trong những công trình nghiên cứu tiên phong trong nghiên cứu về chim đâm được thực hiện bởi Wilbeck[4], ông đề xuất rằng thuyết thủy động lực có thể áp dụng cho vật phóng trong va chạm khi ứng suất của vật phóng lớn hơn rất nhiều so với độ bền của vật liệu làm vật phóng.Ông đã thành công trong việc chứng minh vật phóng có ứng xử như chất lỏng vì độ bền không đáng kể, còn gọi là cơ thể mềm (soft body). Bước đầu trong việc áp dụng thuyết thủy động lực, Wilbeck tách tải trọng va chạm khỏi phản hồi của vật đích bằng cánh giả sử vật đích là tuyệt đối cứng, sau đó bao gồm các hiệu ứng biến dạng của vật đích. Phương pháp này ước tính tải trọng va chạm sử dụng hệ thức sốc 1-D trong chế độ sốc, và giải bằng phương trình Bernoulli nén 1-D cho chế độ dòng chảy ổn định.

Tuy nhiên , với sự phát triển của công nghệ máy tính ngày nay, có thể dễ dàng giải các phương trình điều khiển phi tuyến 3-D. Phần này sẽ giới thiệu về phương 15 Chương II: Lý thuyết trình điều khiển 3-D phi tuyến đầy đủ như bảo toàn khối lượng, bảo toàn động lượng tuyến tính, bảo toàn mô men động lượng và bảo toàn năng lượng.1 Bảo toàn khối lượng Hình 13 Vật thể trước và sau biến dạng Phương trình bảo toàn khối lượng : >h = >0 ///////////@KBAC Với >0 và > là khối lượng riêng ban đầu và sau biến dạng. J là định thức Jacobian, cân bằng với định thức của gradient biến dạng F: mnk h = ija/@;kl C = / /////////////////////@KBKC mol Trong phương trình trên, nk là vị trí hiện tại của điểm, ol /là đại diện cho tọa độ của điểm trong khung gốc hoặc khung tham chiếu. Cả i và j nằm trong khoảng 1 đến 3 trong 3-D.

16 Chương II: Lý thuyết 2.2 Bảo toàn động lượng tuyến tính Với các điều kiện biên phù hợp , bảo toàn động lượng tuyến tính biểu diễn bằng: pkl%l ` />qk = >Tr////////////////////@KBQC Với pkl là ứng suất Cauchy q k lực tác dụng lên mỗi đơn vị khối lượng T k là vận tốc của vật hiện tại Điều kiện biên phù hợp có ba dạng : a. Điều kiện biên kéo : pkl sl = tk ////////////////////@KBUC Trong đó tk là thành phần vector kéo sl là thành phần tiếp tuyến mặt phẳng. Điều kiện biên chuyển vị : nk ok % t = uk t /////////////////////@KB\C Với uk là các thành phần của vector chuyển vị c. Điều kiện biên tiếp xúc Tại mặt phẳng tiếp xúc, thành phần chuyển vị của hai mặt phẳng tiếp xúc phải đáp ứng các rằng buộc : x nkv w nk svk y z////////////////////@KB_C α và β là hai mặt phẳng tiếp xúc Trong phương trình trên, khi các ràng buộc cân bằng, thành phần tiếp tuyến của kéo trên các bề mặt tiếp xúc phải bằng nhau.

Điều kiện này nghĩa là : x x x pklv skv slv w pkl sk sl = z////////////////////@KBcC 17 Chương II: Lý thuyết 2.3 Bảo toàn momem động lượng Trong trường hợp không có cặp vật thể, bảo toàn momen động lượng là tensor ứng suất Cauchy đối xứng, pkl = pkl ///////////////////@KBfC 2.4 Bảo toàn năng lượng Phương trình bảo toàn năng lượng dung để tính toán nội năng e sử dụng trong phương tình trạng thái để thu được mối quan hệ áp suất- khối lượng riêng của vật liệu. Bảo toàn năng lượng có thể viết dưới dạng: >N = pkl 1r{ ` >qk Tk ////////////////////@KBgC Với 1r{ là tensor tỉ số biến dạng Tỉ số biến dạng (strain) giống với tỉ số thay đổi hình dạng (deformation) ukl. Thu được từ gradient vận tốc mTk R kl = /////////////////////////@KBAzC mnl A ukl = @R ` Rlk /C////////////@KBAAC K kl Thành phần đối xứng lệch : A ekl = R w Rlk / ////////@KBAKC/ K kl ekl còn được gọi là tensor xoay.