Tổng quan nghiên cứu

Vật liệu composite đóng vai trò then chốt trong kỹ thuật hiện đại với lịch sử ứng dụng hơn 6 thập kỷ trong ngành hàng không vũ trụ và hơn 3 thập kỷ trong kỹ thuật xây dựng công trình nhờ tỷ lệ cường độ trên trọng lượng vượt trội cùng khả năng chống ăn mòn cao. Tuy nhiên, kết cấu composite nhiều lớp rất nhạy cảm với các dạng khuyết tật nội tại, đặc biệt là hiện tượng tách lớp phát sinh từ khâu chế tạo hoặc do tải trọng mỏi trong quá trình vận hành, làm suy giảm từ 20% đến 40% độ cứng tổng thể và đe dọa trực tiếp đến độ an toàn của toàn bộ công trình. Đề tài luận văn thạc sĩ chuyên ngành Xây dựng công trình dân dụng và công nghiệp của học viên Bùi Hồng Châu, dưới sự hướng dẫn khoa học của PGS. Nguyễn Thị Hiền Lương và TS. Bùi Quốc Tính tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh, tập trung nghiên cứu mô hình hóa và giải quyết bài toán truyền sóng nhằm phát hiện khuyết tật trong dầm composite.

Nghiên cứu được triển khai thực hiện từ ngày 18/08/2014 đến ngày 22/12/2014 và bảo vệ thành công vào tháng 01/2015. Mục tiêu trọng tâm của luận văn là phát triển mô hình phần tử hữu hạn phổ dựa trên biến đổi Wavelet để phân tích đặc tính truyền sóng trong dầm composite có tách lớp, đồng thời sử dụng các đáp ứng mô phỏng làm dữ liệu giải bài toán ngược định vị hư hại. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc rút ngắn hơn 50% thời gian tính toán so với phương pháp phần tử hữu hạn cổ điển, nâng cao độ chính xác nhận diện các vết tách lớp có kích thước từ 20mm đến 30mm và kiểm soát sai số thời gian tới của sóng phản xạ dưới mức 3%.

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 xây dựng trên nền tảng kết hợp chặt chẽ giữa cơ học kết cấu composite và toán học ứng dụng hiện đại:

  • Lý thuyết dầm Timoshenko bậc cao: Kế thừa mô hình dầm có xét đến ảnh hưởng của biến dạng trượt và quán tính quay, cho phép mô tả chính xác ứng xử động học của các lớp vật liệu ở dải tần số cao.
  • Lý thuyết phân tích đa phân giải và biến đổi Wavelet Daubechies: Sử dụng các hàm tỉ lệ Daubechies trực giao bậc $N = 4, 6, 12, 22$ để xấp xỉ trường chuyển vị theo thời gian, hỗ trợ giải quyết bài toán biên hữu hạn mà không bị ảnh hưởng bởi tính chu kỳ giả tạo.
  • Mô hình ống dẫn sóng Base-laminate và Sub-laminate: Chia dầm composite tại vùng tách lớp thành các phân đoạn dầm độc lập phía trên và phía dưới, thiết lập hệ cân bằng lực và liên tục chuyển vị tại mặt phân cách.
  • Các khái niệm cốt lõi: Phương pháp phần tử hữu hạn phổ Wavelet, bài toán đa thức trị riêng, hiện tượng tán sắc sóng, số sóng và vận tốc nhóm.

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

Nghiên cứu sử dụng phương pháp số giải tích kết hợp mô phỏng động lực học trên máy tính:

  • Nguồn dữ liệu: Dữ liệu khảo sát được xây dựng từ các thuộc tính cơ lý thực nghiệm chuẩn của vật liệu composite Graphite - Epoxy 4 lớp với mô đun đàn hồi dọc trục $E_{11}$, mô đun đàn hồi ngang $E_{22}$, mô đun trượt $G_{12}$ và mật độ khối xác định.
  • Cỡ mẫu và chọn mẫu nghiên cứu: Luận văn thiết lập 2 kịch bản kích thước tách lớp đại diện gồm $L_d = 20\text{ mm}$ và $L_d = 30\text{ mm}$ dọc theo trục dầm, kết hợp kiểm tra độ hội tụ qua 2 thang lấy mẫu thời gian rời rạc là $T_w = 512$ điểm và $T_w = 1024$ điểm. Phương pháp chọn mẫu có chủ đích giúp bao quát cả trường hợp kết cấu nguyên vẹn lẫn kết cấu có khuyết tật cục bộ ở các vị trí khác nhau.
  • Lý do lựa chọn phương pháp phân tích: Luận văn lựa chọn phương pháp phần tử hữu hạn phổ Wavelet thay vì phương pháp phần tử hữu hạn phổ Fourier truyền thống vì phương pháp phổ Fourier gặp giới hạn nghiêm trọng khi xử lý ống dẫn sóng ngắn và kết cấu hữu hạn không cản, gây ra hiện tượng phản xạ sóng giả định kỳ. Việc áp dụng hàm Daubechies giúp rút gọn chính xác phương trình đạo hàm riêng thành hệ phương trình đạo hàm thường và giải hiệu quả qua bài toán đa thức trị riêng.
  • Timeline nghiên cứu: Quá trình xây dựng mô hình toán, lập trình ma trận độ cứng và xử lý số liệu hoàn thành trong giai đoạn 4 tháng liên tục từ tháng 8 đến tháng 12 năm 2014.

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

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

Nghiên cứu đã chứng minh tính ưu việt vượt bậc của mô hình số thông qua các kết quả cụ thể:

  • Nhận diện chính xác sóng phản xạ từ vết tách lớp: Khi sóng kích thích dạng xung tác dụng dọc trục hoặc ngang dầm đi qua vùng khuyết tật, biên tách lớp tạo ra sóng phản xạ thứ cấp rõ nét; biên độ sóng phản xạ biến thiên tỷ lệ thuận với chiều dài vết tách lớp, làm suy giảm từ 15% đến 25% năng lượng sóng truyền qua khi chiều dài khuyết tật tăng từ $L_d = 20\text{ mm}$ lên $L_d = 30\text{ mm}$.
  • Triệt tiêu hoàn toàn hiện tượng tràn biên chu kỳ: Mô phỏng dầm công xôn với bậc hàm Daubechies $N = 22$ tại cửa sổ thời gian $T_w = 1024$ điểm cho thấy đáp ứng vận tốc dọc trục và chuyển vị ngang đạt độ ổn định tuyệt đối, triệt tiêu 100% hiện tượng phản xạ biên giả mà phương pháp phổ Fourier thường mắc phải.
  • Xác định rõ đặc tính vận tốc nhóm và tần số cắt: Vận tốc nhóm của mode sóng dọc trục trong dầm Graphite - Epoxy duy trì ổn định trong khoảng 4500 m/s đến 5200 m/s, trong khi mode sóng uốn thể hiện tính tán sắc rõ rệt ở dải tần số cao trên 100 kHz.
  • Tối ưu hóa hiệu năng tính toán: Nhờ kỹ thuật ngưng tụ các bậc tự do bên trong vùng tách lớp về 2 nút phần tử phổ chính, mô hình giúp giảm trên 60% số lượng bậc tự do tổng thể so với mô hình phần tử hữu hạn 2D mà vẫn đảm bảo độ chính xác trên 97%.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự thay đổi dạng sóng là do sự thay đổi đột ngột về độ cứng cục bộ và sự lệch trục trung hòa tại các phân đoạn Sub-laminate trong vùng tách lớp. Hiện tượng này làm thay đổi trở kháng cơ học của ống dẫn sóng, biến đổi một phần sóng dọc trục thành sóng uốn tán sắc.

So sánh với các nghiên cứu trước đây của Mujumdar, Suryanarayana và Nag, giải pháp sử dụng hàm Daubechies trong luận văn này đã khắc phục hoàn toàn điểm nghẽn về chiều dài kết cấu hữu hạn. Trong thực tế phân tích, toàn bộ dữ liệu phản xạ sóng có thể được trực quan hóa thông qua biểu đồ lịch sử thời gian của vận tốc hạt và đường cong tán sắc tần số - số sóng. Các biểu đồ này thể hiện rõ độ trễ thời gian tới của từng gói sóng, cung cấp cơ sở dữ liệu vững chắc cho bài toán ngược nhằm định vị tọa độ khuyết tật với độ lệch vị trí dưới 2%.

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

Dựa trên kết quả đạt được, luận văn đưa ra các khuyến nghị ứng dụng thực tiễn mang tính hệ thống:

  • Tích hợp mô-đun thuật toán WSFEM vào hệ thống giám sát sức khỏe kết cấu trực tuyến: Triển khai cài đặt giải thuật vào bộ vi xử lý của các hệ thống đo đạc để đạt mục tiêu phát hiện khuyết tật dầm sàn composite với độ chính xác trên 95%, dự kiến hoàn thành thử nghiệm trong vòng 12 tháng do các kỹ sư kết cấu và chuyên gia kiểm định công trình chủ trì.
  • Mở rộng kiểm chứng thực nghiệm trên kết cấu composite nhiều lớp: Tiến hành đo đạc thực tế trên các tấm dầm composite từ 8 đến 16 lớp có gắn cảm biến áp điện PZT nhằm rút ngắn thời gian xử lý tín hiệu phản xạ xuống dưới 5 giây, thực hiện trong khung thời gian 6 đến 9 tháng bởi các nhóm nghiên cứu phòng thí nghiệm vật liệu.
  • Nâng cấp giải thuật tự động tối ưu hóa bậc Daubechies: Điều chỉnh linh hoạt bậc hàm từ $N = 6$ đến $N = 22$ theo dải tần số kích thích nhằm tiết kiệm thêm 30% dung lượng bộ nhớ tính toán trong vòng 3 tháng, do các lập trình viên phần mềm mô phỏng kỹ thuật thực hiện.
  • Xây dựng bộ hướng dẫn tiêu chuẩn chẩn đoán hư hại cho công trình công nghiệp: Ban hành quy trình chuẩn đoán tách lớp dựa trên chỉ số thời gian tới và vận tốc sóng nhằm giảm tỷ lệ cảnh báo sai lệch xuống dưới 5% trong vòng 18 tháng, do các viện nghiên cứu chuyên ngành phối hợp với cơ quan quản lý tiêu chuẩn kỹ thuật ban hành.

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

Nội dung và phương pháp trong công trình này mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư kết cấu công trình dân dụng và công nghiệp: Tiếp cận phương pháp kiểm tra không phá hủy hiện đại để đánh giá độ an toàn của dầm, cột và sàn composite trong các công trình nhịp lớn hoặc nhà xưởng chịu tải trọng động.
  • Kỹ sư thiết kế và bảo dưỡng ngành hàng không, cơ khí: Vận dụng mô hình ống dẫn sóng để theo dõi vết nứt và tách lớp trên cánh máy bay, cánh quạt tuabin gió hoặc vỏ tàu thuyền composite trong suốt chu kỳ vận hành.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Cơ kỹ thuật: Sử dụng luận văn như tài liệu tham khảo chuyên sâu về kỹ thuật giải phương trình vi phân bằng biến đổi Wavelet Daubechies và phương pháp số tiên tiến.
  • Các đơn vị kiểm định chất lượng và quan trắc công trình: Khai thác nền tảng thuật toán để phát triển giải pháp cảm biến thông minh, giúp tiết kiệm từ 20% đến 30% ngân sách bảo trì định kỳ cho chủ đầu tư.

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

1. Phương pháp WSFEM có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn truyền thống?
Phương pháp WSFEM sử dụng lời giải giải tích của phương trình sóng làm hàm nội suy, cho phép mô phỏng chính xác các dao động tần số cao trên 100 kHz với kích thước phần tử lớn, giảm hơn 70% số lượng nút lưới so với phương pháp phần tử hữu hạn truyền thống.

2. Tại sao luận văn sử dụng hàm tỉ lệ Daubechies thay vì biến đổi Fourier nhanh?
Biến đổi Fourier giả định tín hiệu có tính tuần hoàn vô hạn nên gây ra hiện tượng nhiễu sóng giả định kỳ ở kết cấu ngắn. Hàm Daubechies có giá mang hữu hạn giúp giải quyết triệt để bài toán biên trong miền thời gian thực với độ chính xác cao.

3. Làm thế nào để xác định chính xác tọa độ vị trí vết tách lớp trong dầm?
Vị trí tách lớp được tính toán dựa trên tích số giữa vận tốc nhóm của sóng truyền và khoảng thời gian từ lúc phát xung đến khi cảm biến ghi nhận sóng phản xạ từ mép khuyết tật, sai số đo lường trong mô hình đạt dưới 2%.

4. Mô hình Base-laminate và Sub-laminate xử lý vùng tách lớp như thế nào?
Mô hình tách dầm tại vị trí khuyết tật thành hai phân đoạn riêng biệt trên và dưới mặt nứt, thiết lập 3 bậc tự do tại mỗi nút và sử dụng ma trận ngưng tụ để liên kết điều kiện biên với phần dầm nguyên vẹn.

5. Thuật toán trong luận văn có thể mở rộng cho các dạng kết cấu khác không?
Hoàn toàn có thể mở rộng. Giải thuật giải bài toán sóng bằng hàm Daubechies có thể áp dụng hiệu quả cho tấm composite 2 chiều, vỏ mỏng nhiều lớp hoặc kết cấu ống dẫn chịu áp lực cao trong thực tế.

Kết luận

  • Luận văn đã xây dựng thành công mô hình phần tử hữu hạn phổ Wavelet hoàn chỉnh để phân tích sóng truyền trong dầm composite có tách lớp.
  • Khắc phục triệt để các hạn chế về tràn biên và tính chu kỳ giả của phương pháp phần tử hữu hạn phổ Fourier trên kết cấu hữu hạn không cản.
  • Thiết lập quy trình chẩn đoán khuyết tật dựa trên thời gian tới của sóng phản xạ với sai số định vị khuyết tật thực tế dưới 2%.
  • Hoàn thành toàn diện các mục tiêu học thuật và thực nghiệm số trong khung thời gian 4 tháng với độ tin cậy và tính khoa học cao.
  • Đặt nền móng vững chắc cho hướng nghiên cứu ứng dụng phương pháp phần tử hữu hạn phổ Wavelet trong giám sát sức khỏe kết cấu tại Việt Nam.

Công trình đóng góp một công cụ tính toán mạnh mẽ cho lĩnh vực chẩn đoán công trình; các nhà nghiên cứu và kỹ sư quan tâm có thể tiếp tục mở rộng mô hình cho bài toán tấm vỏ không gian trong giai đoạn 12 tháng tiếp theo. Hãy tham khảo toàn văn nghiên cứu để khai thác các ma trận độ cứng phổ và thuật toán xử lý tín hiệu chuyên sâu.