Luận án tính toán vỏ thoải composite lớp chịu tác dụng đồng thời tải trọng sóng xung kích và nhiệt độ

Vật liệu tiên tiến đóng vai trò then chốt trong kỹ thuật kết cấu hiện đại. Các nghiên cứu gần đây tập trung mạnh vào kết cấu vỏ thoải FGM và tấm composite nhiều lớp. Vật liệu FGM có cơ tính biến thiên liên tục giúp giảm thiểu tập trung ứng suất tại mặt phân giới. Kết cấu composite lớp lại cho phép gia cường hướng chịu lực cục bộ theo yêu cầu. Nghiên cứu thực nghiệm kết hợp mô phỏng số đang là xu thế chủ đạo toàn cầu. Nhiều công trình tập trung giải quyết bài toán cơ nhiệt phi tuyến phức tạp. Việc đánh giá độ bền khi chịu tải động lực học nổ mang ý nghĩa then chốt. Cơ sở lý thuyết vững chắc giúp nâng cao tuổi thọ vận hành thực tế.

Tải trọng nổ sinh ra sóng áp suất biến đổi đột ngột theo thời gian. Sóng xung kích truyền qua môi trường tạo ra xung lực cực lớn lên bề mặt tấm vỏ. Biên độ áp suất đỉnh tăng vọt trong micro giây rồi suy giảm nhanh chóng. Hiện tượng này kích thích các bậc dao động bậc cao trong kết cấu mỏng. Vỏ thoải dễ xuất hiện biến dạng lớn cục bộ và hư hại vật liệu. Việc mô hình hóa chính xác hàm áp suất theo thời gian là điều kiện tiên quyết. Các mô hình sóng Friedlander cải tiến thường được áp dụng để mô tả quá trình nổ. Phân tích đúng tải trọng nổ giúp bảo vệ an toàn cho kết cấu bảo vệ quân sự.

Môi trường nhiệt độ cao làm thay đổi đáng kể đặc tính vật liệu composite. Ứng xử cơ nhiệt trở nên phi tuyến phức tạp dưới tác động nhiệt cục bộ. Gradien nhiệt gây ra ứng suất nhiệt nội tại và làm suy giảm mô đun đàn hồi. Khi kết hợp cùng tải trọng va đập động, nguy cơ phá hủy tăng gấp nhiều lần. Sự bất đối xứng trong trường nhiệt độ làm tăng độ cong vênh bề mặt. Bài toán ghép trường cơ và nhiệt đòi hỏi các thuật toán tính toán chính xác cao. Việc phân lập ảnh hưởng của nhiệt độ và sóng áp lực giúp tối ưu hóa lớp cách nhiệt. Đây là tiền đề cho giải pháp kết cấu composite sandwich bền vững.

Mô hình chuyển vị được xây dựng dựa trên lý thuyết biến dạng cắt bậc nhất FSDT. Lý thuyết này tính đến ảnh hưởng của góc quay pháp tuyến và biến dạng trượt ngang. Biến dạng cắt ngang đóng vai trò đặc biệt quan trọng đối với vật liệu composite nhiều lớp. Trường chuyển vị được biểu diễn qua năm thành phần chuyển vị độc lập tại mặt trung hòa. Sự đơn giản hóa này duy trì độ chính xác cao nhưng tiết kiệm đáng kể tài nguyên tính toán. Quan hệ biến dạng và chuyển vị phi tuyến Von-Kármán được tích hợp chặt chẽ. Hệ thức chuyển vị đáp ứng hoàn hảo cho các bài toán vỏ thoải có độ cong ban đầu nhỏ.

Quan hệ ứng suất và biến dạng tuân theo định luật Hooke tổng quát cho vật liệu trực giao. Ảnh hưởng của nhiệt độ được biểu diễn qua thành phần giãn nở nhiệt theo từng hướng sợi. Ứng suất nhiệt xuất hiện khi chuyển vị tự do bị cản trở bởi liên kết biên. Các thành phần nội lực bao gồm lực kéo nén, mô men uốn và lực cắt trượt ngang. Tích phân ứng suất dọc theo chiều dày vỏ xác định ma trận độ cứng liên kết A, B, D. Hiện tượng phi tuyến hình học được xem xét đầy đủ để phản ánh đúng biến dạng lớn. Quan hệ ứng xử cơ nhiệt chính xác giúp dự đoán chính xác trạng thái phá hủy của vật liệu.

Phương trình năng lượng được thiết lập dựa trên nguyên lý công ảo Lagrange. Trường nhiệt độ phân bố qua chiều dày vỏ được giả thiết dạng tuyến tính hoặc phi tuyến. Công của ngoại lực bao gồm áp lực sóng xung kích và tải nhiệt độ bề mặt. Thế năng biến dạng đàn hồi tích lũy toàn bộ năng lượng biến dạng cơ và nhiệt. Động năng của kết cấu tính đến vận tốc chuyển vị tịnh tiến và quán tính quay của các phần tử. Sự cân bằng năng lượng toàn phần dẫn đến hệ phương trình vi phân dao động hoàn chỉnh. Phương pháp này đảm bảo tính ổn định và hội tụ của các thuật toán số tiếp theo.

Liên kết giữa các lớp sandwich được mô phỏng thông qua mô hình lớp đàn hồi tương đương. Hệ số độ cứng cắt của chốt liên kết đóng vai trò truyền tải lực cắt giữa các lớp. Mô hình xem xét chuyển vị trượt tương đối tại các bề mặt tiếp xúc. Chuyển vị trượt gây ra tổn hao năng lượng và làm thay đổi tần số dao động riêng. Thế năng biến dạng của hệ thống chốt được bổ sung trực tiếp vào thế năng tổng thể. Thuật toán cho phép điều chỉnh mật độ chốt và góc bố trí chốt linh hoạt. Giải pháp liên kết chốt gia tăng độ cứng chống uốn cho panel sandwich đáng kể.

Tác động đồng thời của nhiệt độ cao và xung lực nổ dễ gây mất ổn định động lực học. Tải trọng nhiệt tạo ra lực nén trước trong mặt phẳng trung hòa của vỏ. Khi biên độ sóng áp suất vượt ngưỡng tới hạn, kết cấu rơi vào trạng thái mất ổn định cục bộ. Biến dạng chuyển vị tăng vọt ngoài mặt phẳng trong thời gian cực ngắn. Hiện tượng bật ngược động lực học gây nguy cơ sụp đổ hoàn toàn tấm vỏ. Tiêu chuẩn Budiansky-Roth được áp dụng để xác định tải trọng tới hạn mất ổn định. Việc kiểm soát độ cứng cắt của chốt giúp nâng cao ngưỡng ổn định của tấm panel sandwich.

Hệ phương trình vi phân dao động cưỡng bức của panel cầu sandwich có dạng ma trận cấp hai. Các thành phần quán tính, cản đàn hồi và độ cứng phi tuyến được ghép nối chặt chẽ. Lực kích động bao gồm tải trọng sóng xung kích biến thiên và ứng suất nhiệt tĩnh. Độ cong hai chiều của panel cầu tạo ra sự tương tác phức tạp giữa kéo nén và uốn. Việc đưa mô hình chốt liên kết vào phương trình giúp phản ánh sát thực tế làm việc. Phương trình vi phân được chuẩn hóa để thuận tiện cho việc giải số. Nghiên cứu cung cấp công cụ dự báo dao động chính xác cho kết cấu vỏ thoải dạng cầu.

Mô hình phần tử vỏ thoải tám nút hoặc chín nút đẳng tham số được áp dụng phổ biến. Các hàm dạng Lagrange bậc cao bảo đảm tính liên tục của trường chuyển vị qua biên phần tử. Véc tơ biến dạng phần tử được tính toán từ đạo hàm hàm dạng theo tọa độ tự nhiên. Phép tích phân số Gauss toàn phần hoặc rút gọn được sử dụng để tính ma trận phần tử. Mô hình phần tử vỏ thoải composite phản ánh chính xác hiệu ứng cắt ngang bậc nhất FSDT. Các liên kết biên cơ bản như ngàm tựa đơn hoặc tự do được thiết lập linh hoạt. Lưới phần tử mịn tại các vùng tập trung ứng suất giúp tăng độ tin cậy.

Quá trình ghép nối ma trận phần tử tạo nên hệ phương trình dao động toàn cục của kết cấu. Ma trận độ cứng tiếp tuyến tổng hợp gồm độ cứng tuyến tính, phi tuyến và ma trận hình học nhiệt. Ma trận khối lượng kết cấu được tính toán theo mô hình khối lượng phân bố nhất quán. Véc tơ tải trọng nút bao gồm xung lực nổ và tải trọng tương đương do ứng suất nhiệt. Phương trình vi phân mô tả trạng thái dao động phi tuyến được giải theo từng bước thời gian nhỏ. Ma trận suy biến do nhiệt được cập nhật liên tục theo trường nhiệt độ phân bố. Kết quả thu được cung cấp lịch sử chuyển vị và vận tốc dao động chi tiết.

Độ tin cậy của thuật toán được kiểm chứng thông qua so sánh với các công bố quốc tế uy tín. Các bài toán chuẩn bao gồm panel trụ và panel cầu composite chịu tải tĩnh và tải động. Sai số giữa kết quả mô phỏng bằng chương trình tự phát triển và tài liệu tham khảo dưới 3%. Tần số dao động riêng và biên độ đáp ứng chuyển vị có sự trùng khớp rất cao. Chương trình tính toán thể hiện sự ổn định tuyệt đối trong cả miền đàn hồi tuyến tính và phi tuyến. Khả năng mô phỏng chính xác tạo cơ sở vững chắc cho các nghiên cứu khảo sát tham số mở rộng. Độ tin cậy thuật toán được khẳng định hoàn toàn.

Tỷ số giữa cạnh và chiều dày (a/h) ảnh hưởng trực tiếp đến độ cứng chống uốn của vỏ thoải. Khi tỷ số a/h tăng, kết cấu trở nên mỏng hơn và tần số dao động riêng giảm mạnh. Tỷ số cạnh a/b quyết định dạng dao động riêng và vị trí tập trung ứng suất uốn cực đại. Tỷ số chiều dày lớp lõi và lớp mặt (hc/ht) chi phối năng lượng hấp thụ biến dạng trượt. Việc tăng tỷ số hc/ht giúp giảm khối lượng tổng thể nhưng vẫn duy trì độ bền uốn mong muốn. Sự kết hợp kích thước hình học hợp lý giúp tránh hiện tượng cộng hưởng nguy hiểm với tần số sóng xung kích.

Độ cứng cắt ks của chốt liên kết đóng vai trò then chốt trong việc kiểm soát chuyển vị trượt. Khi tăng hệ số ks, sự liên kết giữa các lớp trở nên hoàn hảo, làm tăng tần số dao động riêng. Ngược lại, chốt có độ cứng cắt thấp làm tăng biến dạng trượt và giảm độ cứng tổng thể. Tỷ số mô đun đàn hồi giữa lớp ngoài và lớp giữa (Ec/Et, Eb/Et) ảnh hưởng sâu sắc đến phân bố ứng suất. Lớp mặt có mô đun đàn hồi cao giúp chịu tải nén và uốn chủ đạo hiệu quả. Tối ưu hóa thông số chốt và vật liệu mặt giúp nâng cao khả năng chịu lực toàn diện.

Tải nổ kết hợp môi trường nhiệt độ cao gây ra đáp ứng động lực học phức tạp nhất. Biên độ chuyển vị tại tâm vỏ đạt giá trị đỉnh chỉ sau vài mili giây kích nổ. Nhiệt độ môi trường tăng làm kéo dài chu kỳ dao động và gia tăng chuyển vị dư của vỏ thoải. Hiện tượng dao động phi tuyến thể hiện rõ nét qua sự dịch chuyển của đường cân bằng động. Sự hiện diện của chốt liên kết giúp giảm chấn động và hạn chế tối đa nguy cơ phá hủy tách lớp. Kết quả khẳng định tầm quan trọng của việc gia cường chốt trong môi trường nhiệt và tải trọng va đập nổ.