Phân tích động lực học tấm composite áp điện có gân gia cường chịu tải trọng khí động

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG HỌC VIỆN KỸ THUẬT QUÂN SỰ Nguyễn Ngọc Thủy PHÂN TÍCH ĐỘNG LỰC HỌC TẤM COMPOSITE ÁP ĐIỆN CÓ GÂN GIA CƯỜNG CHỊU TẢI TRỌNG KHÍ ĐỘNG LUẬN ÁN TIẾN SĨ KỸ THUẬT Hà Nội - 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG HỌC VIỆN KỸ THUẬT QUÂN SỰ Nguyễn Ngọc Thủy PHÂN TÍCH ĐỘNG LỰC HỌC TẤM COMPOSITE ÁP ĐIỆN CÓ GÂN GIA CƯỜNG CHỊU TẢI TRỌNG KHÍ ĐỘNG Chuyên ngành: Cơ kỹ thuật Mã ngành: 9.52.01.01 LUẬN ÁN TIẾN SĨ KỸ THUẬT Người hướng dẫn khoa học: GS.TS Nguyễn Thái Chung Hà Nội - 2020 LỜI CAM ĐOAN Tôi Nguyễn Ngọc Thủy xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết luận án trung thực chưa công bố công trỡnh no Tỏc gi Nguyn Ngc Thy LờI CảM ƠN Tác giả luận án xin bày tỏ lòng biết ơn chân thành GS.TS Nguyễn Thái Chung đà tận tình hớng dẫn, giúp đỡ cho nhiều dẫn khoa học có giá trị giúp cho tác giả hoàn thành luận án Tác giả trân trọng động viên, khuyến khích kiến thức khoa học nh chuyên môn mà Thầy hớng dẫn đà chia sẻ cho tác giả năm qua, giúp cho tác giả nâng cao lực chuyên môn phơng pháp nghiên cứu Tác giả trân trọng cảm ơn LÃnh đạo Học viện Kỹ thuật Quân sự, tập thể Bộ môn Cơ học vật rắn, Khoa Cơ khí, Phòng Sau đại học - Học viện Kỹ thuật quân sự, Phòng thí nghiệm Trọng điểm Quốc gia/ Đại học Bách khoa/Đại học Quốc gia TP Hồ Chí Minh đà tạo điều kiện thuận lợi, hợp tác trình nghiên cứu Tác giả xin trân trọng cảm ơn GS.TS.NGND Hoàng Xuân Lợng - Học viện Kỹ thuật Quân sự, GS.TSKH.NGND Đào Huy Bích - Đại học quốc gia Hà Nội, GS.TS Trần ích Thịnh - Đại học Bách khoa Hà Nội đà cung cấp cho tác giả nhiều tài liệu quý hiếm, kiến thức khoa học đại nhiều lời khuyên bổ ích, dẫn khoa học có giá trị để NCS hoàn thành luận án Cuối cùng, tác giả xin bày tỏ lòng cảm ơn Thủ trởng Tổng cục Công nghiệp Quốc phòng, Thủ trởng Cục Quản lý Công nghệ/Tổng cục CNQP đồng nghiệp, với ngời thân gia đình đà thông cảm, động viên chia sẻ khó khăn với tác giả suốt thời gian nghiên cứu để hoàn thành luận án Tác giả Nguyễn Ngọc Thủy iii MC LC Trang Lời cam đoan i Mục lục iii Danh mục ký hiệu, chữ viết tắt vii Danh mục bảng xiii Danh mục hình vẽ, đồ thị xiv MỞ ĐẦU CHƯƠNG TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Sơ lược tượng áp điện ứng dụng kỹ thuật 1.2 Tổng quan kết cấu composite áp điện 1.3 Các mơ hình lực khí động sử dụng tính tốn kết cấu 1.4 Tổng quan tình hình nghiên cứu kết cấu composite áp điện 11 1.5 Các kết đạt từ cơng trình công bố 20 1.6 Các vấn đề cần tiếp tục nghiên cứu 21 1.7 Kết luận rút từ tổng quan 22 CHƯƠNG PHÂN TÍCH PHI TUYẾN ĐỘNG LỰC HỌC TẤM COMPOSITE ÁP ĐIỆN CĨ GÂN GIA CƯỜNG CHỊU TẢI TRONG KHÍ ĐỘNG .23 2.1 Đặt vấn đề 23 2.2 Đặt toán giải thiết 23 2.3 Quan hệ ứng xử học composite lớp có gân gia cường .24 2.3.1 Quan hệ ứng xử học composite lớp 24 2.3.1.1 Quan hệ biến dạng chuyển vị 25 2.3.1.2 Quan hệ ứng suất biến dạng 30 2.3.1.3 Các thành phần nội lực 31 2.3.1.4 Các quan hệ ứng xử học composite lớp 32 2.3.2 Quan hệ ứng xử học gân gia cường 34 iv 2.3.2.1 Trường chuyển vị 34 2.3.2.2 Trường biến dạng 34 2.3.2.3 Trường ứng suất 36 2.4 Quan hệ ứng xử học composite áp điện có gân gia cường 36 2.4.1 Ứng xử lớp áp điện 37 2.4.2 Ứng xử n lớp composite m lớp áp điện 39 2.5 Thiết lập phương trình vi phân phi tuyến mô tả dao động phần tử composite áp điện có gân gia cường chịu tải trọng khí động 40 2.5.1 Phần tử composite có lớp áp điện 40 2.5.2 Phần tử composite áp điện có gân gia cường 55 2.5.3 Phần tử composite áp điện có gân gia cường chịu tải trọng khí động 60 2.5.4 Xây dựng ma trận tổng thể kết cấu từ ma trận phần tử .63 2.5.4.1 Ma trận tổng thể 63 2.5.4.2 Véc tơ tải trọng tổng thể 64 2.5.4.3 Phương trình mơ tả dao động hệ 65 2.6 Thuật tốn PTHH giải phương trình dao động composite áp điện có gân gia cường chịu tải trọng khí động 65 2.6.1 Bài toán dao động tự 65 2.6.2 Bài toán dao động cưỡng 66 2.7 Phân tích ổn định composite áp điện có gân gia cường chịu tác dụng lực khí động 70 2.7.1 Tiêu chuẩn ổn định động Budiansky-Roth 70 2.7.2 Phân tích ổn định composite áp điện có gân gia cường chịu tải trọng khí động 71 2.8 Giới thiệu chương trình kiểm tra mức độ tin cậy 71 2.8.1 Giới thiệu chương trình tính 71 2.8.2 Kiểm tra độ tin cậy chương trình 72 v 2.9 Kết luận chương 74 CHƯƠNG KHẢO SÁT SỐ VÀ THẢO LUẬN 75 3.1 Đặt vấn đề 75 3.2 Bài toán xuất phát 75 3.3 Ảnh hưởng số yếu tố đến dao động ổn định composite áp điện có gân gia cường chịu tải trọng khí động 81 3.3.1 Ảnh hưởng vận tốc dịng khí 81 3.3.2 Ảnh hưởng lớp áp điện 84 3.3.3 Ảnh hưởng tỷ số chiều dày bề rộng (h/W) .87 3.3.4 Ảnh hưởng điều kiện liên kết 90 3.3.5 Ảnh hưởng góc đặt cốt 94 3.3.6 Ảnh hưởng góc tới 97 3.3.7 Ảnh hưởng kích thước gân 103 3.3.8 Ảnh hưởng điện áp V áp đặt lên lớp áp điện 108 3.3.9 Ảnh hưởng kích thước áp điện 111 3.3.10 Ảnh hưởng tính chất cản 115 3.3.11 Miền ổn định điện áp góc đặt cốt thay đổi 118 3.4 Kết luận chương 120 CHƯƠNG NGHIÊN CỨU THỰC NGHIỆM 122 4.1 Đặt vấn đề 122 4.2 Mơ tả thí nghiệm 122 4.2.1 Thiết lập thí nghiệm 122 4.2.2 Thiết bị đo đáp ứng động 124 4.2.3 Bộ phát tín hiệu chuẩn 124 4.2.4 Bộ khuếch đại piezo tuyến tính EPA-104-230 125 4.2.5 Bộ thu nhận liệu kênh áp điện HnB75B 126 4.2.6 Máy sóng 126 vi 4.2.7 Thiết bị tạo gió (hầm gió) 127 4.3 Phương pháp xác định gia tốc, biến dạng kết cấu 128 4.3.1 Đo gia tốc 128 4.3.2 Đo biến dạng 129 4.4 Cơ sở phân tích xử lý kết thí nghiệm 129 4.5 Thí nghiệm kết thu 130 4.6 Kết luận chương 136 KẾT LUẬN VÀ KIẾN NGHỊ 138 DANH MỤC CÔNG TRÌNH CỦA TÁC GIẢ 140 TÀI LIỆU THAM KHẢO 141 PHỤ LỤC 154 Phụ lục Các biểu thức kết thí nghiệm chương 157 Phụ lục Mã nguồn chương trình SMART_STIFFENED_PLATE_2018 (SSP_2018) 170 vii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT Danh mục ký hiệu 1.1 Các ký hiệu chữ La tinh [A] Ma trận độ cứng màng composite có lớp áp điện [A]c Ma trận độ cứng màng composite Ae Ae Công ngoại lực phần tử Diện tích phần tử [B] Ma trận độ cứng màng - uốn - xoắn composite có lớp áp điện c Ma trận độ cứng màng - uốn - xoắn composite [B] L Ma trận quan hệ biến dạng tuyến tính – chuyển vị nút phần tử [B ] N Ma trận quan hệ biến dạng phi tuyến – chuyển vị nút phần tử [Bφ ] Ma trận tính điện trường [B ] ’ [C], [C ] Ma trận số độ cứng composite hệ tọa độ 1, 2, x, y, z Hệ số áp lực khí động Cp  C  Ma trận hệ số độ cứng vật liệu áp điện   C   A  [D] [D]  CeA  C   R   CeR   Ceair  c 160 Bảng Giá trị Max/Min gia tốc theo thời gian Giá trị gia tốc max/min vị trí gia tốc với f = 6.944 Hz, U = 10m/s 2.2 Trường hợp 2: - Khi đặt lệch góc 22,5°, vận tốc gió U = 5m/s 10 m/s Điện áp kích thích Vin = 9,50 V (đỉnh), tần số kích thích f = 6.981Hz, Kamp= 20, tần số lấy mẫu fsamp = 1000Hz Điện áp cấp 190.0 V Kết lần đo sau: 161 Hình Đáp ứng gia tốc theo thời gian lần đo liên tiếp với f = 6.981 Hz, U = 5m/s 162 Hình Đáp ứng gia tốc theo thời gian lần đo liên tiếp với f = 6.981 Hz, U= 10m/s 163 Bảng Giá trị Max/Min gia tốc theo thời gian Giá trị gia tốc max/min vị trí đáp ứng gia tốc với f = 6.981 Hz, U = 5m/s 164 Giá trị gia tốc max/min vị trí đáp ứng gia tốc với f = 6.981 Hz, U = 10 m/s 2.3 Trường hợp 3: Khi lệch góc 45 , vận tốc gió U = 5m/s 10m/s Điện áp kích Vin = 9,20V, tần số kích thích f = 6,993Hz, Kamp= 20, tần số lấy mẫu fsamp = 1000Hz Điện áp cấp 184,0V Kết lần đo sau: 165 Hình Đáp ứng gia tốc theo thời gian lần đo với f = 6,993Hz, U = 5m/s 166 Hình Đáp ứng gia tốc theo thời gian lần đo liên tiếp với f = 6.993 Hz, U= 10m/s 167 Bảng Giá trị Max/min theo thời gian đại lượng đo trên: Giá trị gia tốc max/min vị trí đáp ứng gia tốc với f = 6.993 Hz, U = 5m/s 168 Giá trị gia tốc max/min vị trí đáp ứng gia tốc với f = 6.993 Hz, U = 10 m/s 169 Phụ lục %CHUONG TRINH PHAN TICH DAO DONG CUA TAM COMPOSITE AP DIEN CO GAN GIA CƯƠNG % SMART_STIFFENED_PLATE_2018 (SSP-2018) - clear all; echo off; Thong so vao % % - For the plate and Graphite-Epoxy T300/976,PZT b=0.50; % Length of the plate[m] a=0.40; % Width of the plate [m] t_pzt=0.000075; % Thickness of the PZT layers) [m] t_cps=0.0008; % Thickness of the composite layers [m] n_cps=4; % Số lớp composite n_pzt=2; % Số lớp áp điện theta=45; h=n_cps*t_pzt+n_pzt*t_pzt; % Chiều dày composite bg=0.003; % Chiều rộng gân [m] hg=0.0048; % Chiều cao gân [m] V=50 %Voltage applied U=20; % Tốc độ gió Alfa=30; % Góc nghiêng composite deltat=0.001 %Time step t_tol=4.0; Tong thoi gian tinh deltat=0.005; Buoc thoi gian tich phan n_f=20; So luong tan so rieng xuat % -Graphite-Epoxy T300/976 E11=150; %Gpa E22=9; %Gpa E33=E22; %Gpa G12=7.1; G13=G12; %Gpa G23=2.5; - ν12=0.3; ν23=ν12; ν32=ν23; ρGE=1600; %kg/m3 PZT_5A %Gpa %Gpa % - E=63.0; G=24.2; ν=0.3; ρpzt=7600; d31=2.54*1e-10; d32=d31; p11=15*1e-9; p22=p11; p33=p22; %kg/m3 %m/V %m/V %F/m %F/m %F/m 170 % - Modeling - global nNode % So nut cua ket cau nDof % So bac tu cua ket cau nDof1 % So phan tu nElem Coords % Bang toa nut Dof % Bang danh so bac tu cua nut %(danh so lai bang cach toi uu bang nay) Edof % Bang danh so bac tu cua phan tu b2 % Gia tri ban dau nhist % Cac bac tu khao sat lay so lieu dau nMode % So dang dao dong rieng can phan tich % -% Cac bien trung gian % -function [Ne,Nex,Ney,Nexx,Neyy,Nexy]=plateshape(ex,ey,x,y) % -a=ex(3)-ex(1); b=ey(3)-ey(1); N=[1 x y x^2 x*y y^2 x^3 x^2*y x*y^2 y^3 x^3*y x*y^3]; Nx=[0 2*x y 3*x^2 2*x*y y^2 3*x^2*y y^3]; Nxx=[0 0 0 6*x 2*y 0 6*x*y 0]; Nxy=[0 0 0 2*x 2*y 3*x^2 3*y^2]; Ny=[0 x 2*y x^2 2*x*y 3*y^2 y 3*x*y^2]; Nyy=[0 0 0 0 2*x 6*y 6*x*y]; a2=a*a;a3=a2*a;b2=b*b;b3=b2*b;ab=a*b; % Ma tran C C=[ -a -b a2 ab b2 -a3 -a2*b -a*b2 -b3 a3*b a*b3; 0 -a -2*b a2 2*ab 3*b2 -a3 -3*a*b2; -1 2*a b -3*a2 -2*ab -b2 3*a2*b b3; a -b a2 -ab b2 a3 -a2*b a*b2 -b3 -a3*b -a*b3; 0 a -2*b a2 -2*ab 3*b2 a3 3*a*b2; -1 -2*a b -3*a2 2*ab -b2 3*a2*b b3; a b a2 ab b2 a3 a2*b a*b2 b3 a3*b a*b3; 0 a 2*b a2 2*ab 3*b2 a3 3*a*b2; -1 -2*a -b -3*a2 -2*ab -b2 -3*a2*b -b3; -a b a2 -ab b2 -a3 a2*b -a*b2 b3 -a3*b -a*b3; 0 -a 2*b a2 -2*ab 3*b2 -a3 -3*a*b2; -1 2*a -b -3*a2 2*ab -b2 -3*a2*b -b3]; E A rho nuy h end global % Cac ma tran cung, khoi luong %va can cua tam Piezo Surf % Index cua cac phan tu be mat chua ap dien Ex Ey % ep i 171 Elem % Index cua Node theo Phan tu mP cP kP % -function []=datain(filename) %Doc so lieu ket cau tu file vao cac bien tong the khai_bao_chung; % -%Doc file so lieu nut % -%fid = fopen(cat(2,filename,'.dat'),'r'); fid = fopen(cat(2,filename,'.txt'),'r'); if fid