Luận án sẽ tập trung giải quyết các nội dung chủ yếu sau: Nghiên cứu tổng quan về hiện tượng áp điện, ứng dụng, kết cấu composite áp điện, các mo hình lực khí động sử dụng trong tính toán kết cấu, tổng quan về tính toán kết cấu composite áp điện làm cơ sở giải quyết bài toán phân tích động lực học tấm composite áp điện có gân gia cườ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 TĨM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT Hà Nội – 2020 CƠNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI HỌC VIỆN KỸ THUẬT QUÂN SỰ - BỘ QUỐC PHÒNG Người hướng dẫn khoa học: GS.TS Nguyễn Thái Chung Phản biện 1: GS.TSKH Nguyễn Tiến Khiêm Phản biện 2: GS.TS Trần Ích Thịnh Phản biện 3: GS.TS Trần Minh Tú Luận án bảo vệ Hội đồng đánh giá luận án cấp Học viện theo định số 2880/QĐ-HV, ngày 03 tháng năm 2020 Giám đốc Học viện Kỹ thuật Quân sự, họp Học viện Kỹ thuật Quân vào hồi ngày tháng năm 2020 Có thể tìm hiểu luận án tại: - Thư viện Học viện Kỹ thuật Quân - Thư viện Quốc gia CÁC CƠNG TRÌNH ĐÃ CÔNG BỐ Nguyễn Thái Chung, Nguyễn Ngọc Thủy (2016), Nghiên cứu ổn định tuyến tính composite áp điện có gân gia cường chịu tác dụng tải trọng điều hịa, Tuyển tập cơng trình Hội nghị Khoa học toàn quốc “Vật liệu kết cấu Composite: Cơ học, Công nghệ Ứng dụng”, Đại học Nha Trang, trang 89-96 (ISBN:978-604-82-2026-6) Nguyễn Thái Chung, Nguyễn Ngọc Thủy (2017), Nghiên cứu điều khiển dao động tự composite áp điện có gân gia cường sử dụng lý thuyết bậc cao, Tuyển tập cơng trình Hội nghị Cơ học toàn quốc lần thứ X, Học viện KTQS, Tập 3, Quyển 1, trang 163-170 (ISBN: 978-604-913-721-1) Nguyễn Ngọc Thủy, Nguyễn Thái Chung (2017), Phân tích động lực học composite lớp áp điện có gân gia cường sử dụng lý thuyết bậc cao, Tuyển tập công trình Hội nghị Cơ học tồn quốc lần thứ X, Học viện KTQS, Tập 3, Quyển 2, trang 1197-1204 (ISBN: 978-604-913-722-8) Nguyễn Thái Chung, Nguyễn Ngọc Thủy (2018), Flutter Analysis for Smart Stiffened Plates using Distributed Piezoelectric Sensors and Actuators, Tạp chí Xây dựng Việt Nam, số 608, trang 208-211 (ISSN 0866-8762) Nguyen Thai Chung, Nguyen Ngoc Thuy, Duong Thi Ngoc Thu and Le Hai Chau (2019), Numerical and Experimental Analysis of the Dynamic Behavior of Piezoelectric Stiffened Composite Plates Subjected to Airflow, Mathematical Problems in Engineering (Article ID: 2697242), 10 page (SCIE) Nguyễn Ngọc Thủy, Nguyễn Thái Chung (2019), Experimental Analysis of the Dynamic Behaviors of Piezoelectric Stiffened Composite Plates, Tạp chí Nghiên cứu KH&CN Quân sự, số 59, trang 196-204 (ISSN 1859-1043) MỞ ĐẦU Tính cấp thiết đề tài Việc nghiên cứu tính tốn, ứng dụng kết cấu dạng tấm, dầm, vỏ composite áp điện nhiều nhà khoa học giới nước quan tâm nghiên cứu, tập trung chủ yếu vào hướng: tính tốn bền, dao động ổn định Tuy nhiên, việc ứng dụng kết cấu dạng tấm, dầm, vỏ composite có lớp áp điện, kết cấu có biện pháp gia cường vào lĩnh vực kỹ thuật đến cịn nghiên cứu tiềm sử dụng chúng lớn Do vậy, vấn đề “Phân tích động lực học composite áp điện có gân gia cường chịu tải trọng khí động” mà tác giả luận án tập trung giải có ý nghĩa khoa học thực tiễn Đối tượng phạm vi nghiên cứu luận án - Đối tượng nghiên cứu luận án: Tấm composite có gân gia cường gắn lớp miếng áp điện (trong lớp composite vật liệu đồng phương) chịu tải trọng khí động dịng khí chuyển động gây - Phạm vi nghiên cứu luận án: Động lực học phi tuyến ổn định kết cấu composite áp điện có gân gia cường Phương pháp nghiên cứu Kết hợp phương pháp nghiên cứu lý thuyết thực nghiệm Về lý thuyết, sử dụng lý thuyết mỏng Mindlin phương pháp PTHH Về thực nghiệm, xây dựng mơ hình tiến hành thí nghiệm đo đáp ứng động composite áp điện có gân gia cường chịu tải trọng khí động buồng thổi gió Cấu trúc luận án Luận án cấu trúc gồm: Mở đầu, chương, phần kết luận kiến nghị, tài liệu tham khảo phụ lục Mở đầu: Trình bày tính cấp thiết vấn đề nghiên cứu Chương Tổng quan vấn đề nghiên cứu Chương Phân tích phi tuyến động lực học composite áp điện có gân gia cường chịu tác dụng lực khí động Chương Khảo sát số thảo luận Chương Nghiên cứu thực nghiệm Kết luận kiến nghị: Trình bày kết luận án kiến nghị tác giả rút từ nội dung nghiên cứu CHƯƠNG TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU Trình bày tổng quan tượng áp điện, ứng dụng, kết cấu composite áp điện, mơ hình lực khí động, kết tính tốn kết cấu composite áp điện nước nước ngồi Từ cơng trình công bố, sở vấn đề cần tiếp tục nghiên cứu phát triển, tác giả luận án tập trung nghiên cứu: “Phân tích động lực học composite áp điện có gân gia cường chịu tải trọng khí động” Theo đó, luận án tập trung giải nội dung chủ yếu sau: Nghiên cứu tổng quan tượng áp điện, ứng dụng, kết cấu composite áp điện, mo hình lực khí động sử dụng tính tốn kết cấu, tổng quan tính tốn kết cấu composite áp điện làm sở giải tốn phân tích động lực học composite áp điện có gân gia cường Xây dựng thuật tốn PTHH chương trình máy tính phân tích dao động composite áp điện có gân gia cường chịu tác dụng tải trọng khí động theo thời gian, xem xét đến tổng cản hệ: cản áp điện, cản kết cấu Khảo sát số lớp toán khác với thông số thay đổi nhằm lựa chọn thông số hợp lý làm giảm dao động tăng khả ổn định cho composite áp điện có gân gia cường chịu tác dụng tải trọng khí động theo thời gian Nghiên cứu thực nghiệm, xác định phản ứng động composite áp điện có gân gia cường chịu tải trọng khí động buồng thổi CHƯƠNG PHÂN TÍCH PHI TUYẾN ĐỘNG LỰC HỌC CỦA TẤM COMPOSITE ÁP ĐIỆN CÓ GÂN GIA CƯỜNG CHỊU TẢI TRỌNG KHÍ ĐỘNG 2.1 Đặt vấn đề Trong chương này, tác giả luận án tập trung xây dựng mô hình, quan hệ ứng xử - điện, thuật tốn PTHH chương trình tính để giải toán đặt Tác giả luận án sử dụng lý thuyết biến dạng trượt bậc Mindlin quan hệ phi tuyến biến dạng - chuyển vị để thiết lập thuật toán điều khiển dao động composite áp điện sở áp dụng tính chất thuận nghịch vật liệu áp điện tác dụng tải trọng khí động 2.2 Đặt toán giả thiết Xét composite áp điện có gân gia cường chịu tải trọng khí động dịng khí di chuyển với vận tốc U tạo với mặt phẳng xoy góc α gây (Hình 2.1) U z y B I A Tấm composite Kích thích Mặt trung bình x I L Gân gia cường hp hs ` hg Cảm biến bg Hình 2.1 Mơ hình tốn Các giả thiết: Tấm gân liên kết bám dính tuyệt đối, thỏa mãn điều kiện lý thuyết biến dạng trượt bậc Mindlin, vật liệu gân làm việc giới hạn đàn hồi Mỗi lớp gân gia cường vật liệu composite đồng phương Vật liệu áp điện biến dạng bé, đàn hồi tuyến tính, hệ số đàn hồi, hệ số tương tác áp điện hệ số điện môi không thay đổi trình biến dạng Lớp áp điện bám dính tuyệt lớp composite Bỏ qua hiệu ứng rối dịng khí chấp nhận phương pháp độc lập, công tác dụng tải trọng 2.3 Quan hệ ứng xử học composite lớp có gân gia cường 2.3.1 Quan hệ ứng xử học composite lớp: 2.3.1.1 Quan hệ biến dạng chuyển vị: x y {L } N b xy b s s yz xz (2.11) 2.3.1.2 Quan hệ ứng suất biến dạng: c c k x Q Q12 11 y Q12 Q 22 xy Q16 Q 26 yz 0 xz k Q16 Q 26 Q66 0 C44 C45 x y xy C45 yz C55 k xz (2.12) 2.3.1.3 Các thành phần nội lực: a Lực màng: hk n c T c T N (x,y,t) N N N N dz, c c x y xy x y xy k (2.15) k1 h k1 b.Mô men uốn xoắn: c T Mc (x,y,t) Mc Mx My Mxy n hk k1 h c T k z x y xy dz, (2.16) k1 c Lực cắt: c Qx Qc (x,y,t) Qc Qy n hk k1 h k1 c c xz A44 A45 xz dz A A , yz 55 yz 45 k (2.17) 2.3.1.4 Các quan hệ ứng xử học composite lớp: c N A M B Q B D 0 0 F c o N s , (2.25) 2.3.2 Quan hệ ứng xử học gân gia cường 2.3.2.1 Trường chuyển vị: u u x, y, z, t u x, t z y x, t , v v x, y, z, t v x, t z x x, t , w w x, y, z, t w x, t y x x, t , (2.26) 2.3.2.2 Trường biến dạng: L NL Lx xNL NL Lxz xz , L NL xy xy (2.30) 2.3.2.3 Trường ứng suất: L NL x Q11 Q12 Q16 x Q11 Q12 Q16 x L NL Q Q Q Q Q Q xz 12 22 26 xz 12 22 26 xz (2.34) Q L NL xy 16 Q26 Q66 xy Q16 Q 26 Q66 xy 2.4 Quan hệ ứng xử học composite có lớp áp điện 2.4.1 Ứng xử học lớp áp điện: Quan hệ ứng suất - biến dạng: p k C k k ek Ek , (2.35) Khi áp điện tĩnh lên bề mặt lớp áp điện lớp áp điện sinh điện trường E, lúc ứng suất lớp áp điện thứ k có thành phần: ứng suất học ứng suất điện trường gây ra, ta có: (2.39) pk Mk Ek , 2.4.2 Ứng xử n lớp composite m lớp áp điện Phương trình ứng xử học gồm n lớp composite m lớp áp điện sau: N A M B Q 0 Dp e B D 0 e 0 0 F e e o N e , e s p Ek (2.43) 2.5 Thiết lập phương trình vi phân mô tả dao động phi tuyến phần tử composite áp điện có gân gia cường chịu tải trọng khí động 2.5.1 Phần tử composite có lớp áp điện: Phương trình vi phân mơ tả dao động phần tử composite có lớp áp điện sau: 1 e Meuu qeM CeR q eM Keuu Keu K Keu qeM K FeM e u e 1 (2.91) e c K Q , 2.5.2 Phần tử composite áp điện có gân gia cường Với giả thiết liên kết gân lý tưởng, chuyển vị mặt đồng với chuyển vị mặt gân Quan hệ chuyển vị vị trí nút chung thứ i phần tử gân phần tử biểu diễn sau: u 0gxi 1 v 0gxi 0 w 0gxi 0 gxi 0 gyi 0 hay: 0,5 h h gx 0 u 0i 0 0 v0i 0 w 0i 0 xi 0 0 yi q =T q gx i x i (2.107) (2.108) Quan hệ chuyển vị nút phần tử gân song song với trục ox với phần tử có chung nút có dạng: q = T q (2.110) T gx = Tx .diag 9,9 (2.111) gx e với: gx e Quan hệ chuyển vị nút phần tử gân song song với trục oy với phần tử xác định tương tự với trục ox Phương trình vi phân mơ tả dao động phần tử composite áp điện có gân gia cường sau: M M M epg M CeA Cepg q eM K epg uu q e R ME q e Fe (2.124) 2.5.3 Phần tử composite áp điện có gân gia cường chịu tải trọng khí động Trường hợp tổng quát, lực khí động tác dụng lên phần tử gồm hai thành phần: áp lực khí động phân bố tính theo mơ hình Scanlan áp lực khí động phân bố tính theo áp lực gió phân bố vng góc với diện tích Phương trình mơ tả lực khí động tác dụng lên phần tử sau: w * KH1 (K) Ucos B x C Usin * l U cos B KH (K) w a p a U cos * K H (K) x * w B x * KA (K) KA (K) 2 m Ucos B2 U cos Ucos a K 2A*3 (K)x (2.129) Phương trình mơ tả dao động phi tuyến phần tử composite áp điện có gân gia cường chịu tác dụng lực khí động lực học sau: e Kaeir qeM feMan Mepg qeM CeA Cepg Cair q eM Kepg ME uu R (2.138) 2.5.4 Xây dựng ma trận tổng thể kết cấu từ ma trận phần tử Sau tập hợp ma trận, véc tơ tải trọng tổng thể trình bày trên, ta có phương trình mơ tả dao động phi tuyến hình học composite áp điện có gân gia cường chịu tải trọng khí động lực học sau: M q CRAa q KMEa q R (2.145) 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 Luận án tập trung giải dạng tốn: Bài tốn dao động riêng tuyến tính, tần số riêng i, xác định giải phương trình: K LU 2 M U 0, (2.147) Tương ứng với tần số riêng i, véc tơ riêng {qi} lúc xác định phương trình: K LU i2 M U q i 0 (2.148) Với tốn dao động cưỡng phi tuyến, phương trình tác giả giải sở sử dụng phương pháp tích phân trực tiếp Newmark kết hợp với phương pháp lặp Newton – Raphson 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 2.7.1 Tiêu chuẩn ổn định động Budiansky – Roth Nội dung tiêu chuẩn phát biểu: Việc giải phương trình dao động ứng với điều kiện khác tải trọng cho ta đáp ứng động khác nhau, điều kiện đó, biên độ chuyển vị lớn hệ tăng theo thời gian có bước nhảy đột ngột, hệ thực dao động quanh vị trí cân khác với vị trí cân ban đầu (Hình 2.5) kết cấu xem ổn định, tải trọng tương ứng để xẩy tượng gọi lực tới hạn Thời điểm lân cận quanh vị trí xuất bước nhảy đột ngột chuyển vị lớn biểu đồ đáp ứng động theo thời gian gọi thời điểm kết cấu ổn định Hình 2.5 Biểu đồ dấu hiệu ổn định động theo tiêu chuẩn Budiansky - Roth 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 Sử dụng tiêu chuẩn ổn định động Budiansky – Roth [25], để xem xét ổn định động composite áp điện có gân gia cường chịu tải 10 Bài toán dao động riêng: Giải toán dao động riêng, tác giả nhận tần số riêng dạng dao động riêng, tần số riêng toán xuất phát với trường hợp có điện áp (50 V) khơng có điện áp (0V) Bảng 3.1 Bảng 3.1 Tám tần số riêng TT Tần số riêng f [Hz] f1 f2 f3 f4 f5 f6 f7 f8 50v (có áp điện) 18,81 52,96 117,07 185,80 277,02 329,33 398,70 408,79 0v (khơng có áp điện) 17,32 49,66 115,27 178,92 261,33 316,27 384,96 388,48 Bài toán động lực học: Kết đáp ứng độ võng w, vận tốc w , gia điểm A đáp ứng ứng suất điểm B hình 3.3, 3.6 tốc w -3 x 10 1.5 Ung suat XicmaB [N/m2] y Chuyen vi phuong dung tai A [m] -2 -4 0.5 -0.5 -1 -1.5 -6 -8 x 10 0.5 1.5 2.5 Thoi gian t[s] 3.5 -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.3 Đáp ứng độ võng wA Hình 3.6 Đáp ứng ứng suất By Bảng 3.2 Giá trị lớn đại lượng tính Chuyển vị Tham số w A 10-3 [m] Giá trị LN 6,923 Vận tốc A 10-2 w [m/s] 31,970 Gia tốc A w [m/s2] 36,832 Ứng suất y 10 [N/m2] 16,991 Biến dạng y 8,317.10-6 Nhận xét: Kết khảo sát cho thấy giá trị đáp ứng chuyển vị, vận tốc, gia tốc ứng suất có xu hướng giảm dần theo thời gian Điều cho thấy vận tốc dịng khí sinh lực khí động phạm vi làm việc ổn định 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 3.3.1 Ảnh hưởng vận tốc dịng khí Phân tích tốn với vận tốc dịng khí biến thiên từ 15m/s 11 (54km/h), đến 28m/s (100,8km/h) Kết đáp ứng độ võng ứng suất điểm tính hình 3.8 3.11 1.5 U = 15 m/s U = 20 m/s U = 25 m/s U = 28 m/s 0.1 x 10 Ung suat Xicmay [N/m2] 0.05 0.5 B Chuyen vi phuong dung tai A [m] 0.15 -0.5 -1 -1.5 U = 15 m/s U = 20 m/s U = 25 m/s U = 28 m/s -0.05 -2 -0.1 0.5 1.5 2.5 Thoi gian t[s] 3.5 -2.5 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.8 Đáp ứng độ võng w Hình 3.11 Đáp ứng ứng suất y Bảng 3.3 Giá trị lớn đại lượng tính thay đổi U Chuyển vị U [m/s] Vận tốc Gia tốc Ứng suất y 10 10 w [m/s ] w [m/s] [N/m2] 25,458 33,116 13,336 31,970 36,832 16,991 96,806 90,117 40,121 Đáp ứng tăng dần theo thời gian w 10-3[m] 15 20 25 28 5,321 6,923 16,904 Biến dạng y 10 -6 6,644 8,317 17,481 Nhận xét: Ảnh hưởng vận tốc U dịng khí tới dao động lớn, cụ thể: Khi U = 15m/s, dao động có xu hướng tắt dần, chuyển vị lớn wmax = 5,321.10-3 m Khi U = 25m/s, dao động bình ổn (biên độ khơng đổi) chuyển vị w = 16,904.10-3 m Tiếp tục tăng vận tốc U với U = 1m/s, nhận thấy vận Uth = 28m/s, độ võng điểm tính tăng theo thời gian, xem bị ổn định 3.3.2 Ảnh hưởng lớp áp điện Khảo sát composite có gân gia cường chịu tác dụng lực khí động có lớp áp điện khơng có lớp áp điện Kết đáp ứng độ võng ứng suất điểm tính hình 3.13 3.16 1.5 0.15 0.1 Ung suat XicmaBy [N/m2] Chuyen vi phuong dung tai A [m] U = 28 m/s (CAD) U = 24 m/s (KAD) x 10 0.05 0.5 -0.5 -1 -1.5 -0.05 -2 -0.1 0.2 0.4 0.6 0.8 Thoi gian t[s] 1.2 1.4 1.6 1.8 Hình 3.13 Đáp ứng độ võng w -2.5 U = 28 m/s (CAD) U = 24 m/s (KAD) 0.2 0.4 0.6 0.8 Thoi gian t[s] 1.2 1.4 1.6 1.8 Hình 3.16 Đáp ứng ứng suất y 12 Nhận xét: Từ kết khảo sát cho thấy, với trường hợp có lớp áp điện (CAD) composite có gân gia cường ổn định đến vận tốc dịng khí 28 m/s, composite có gân gia cường khơng có lớp áp điện (KAD) ổn định đến vận tốc dịng khí 24 m/s sau ổn định Qua cho thấy, lớp áp điện có tác dụng làm giảm dao động composite có gân gia cường tác dụng vận tốc dịng khí 3.3.3 Ảnh hưởng tỷ số chiều dày bề rộng (h/W) Giải toán với tỷ số h/W thay đổi từ 0,003 đến 0,010 (thay đổi chiều dày tấm) Kết đáp ứng chuyển vị ứng suất với tỷ số h/W khác hình 3.18, 3.21 thể bảng 3.4 0.02 0.01 Ung suat Xicmay [N/m2] B Chuyen vi phuong dung tai A [m] 0.03 -0.01 -0.02 -0.03 h/W h/W h/W h/W h/W -0.04 -0.05 -0.06 0.5 1.5 2.5 Thoi gian t[s] = 0,010 = 0,007 = 0,005 = 0,004 = 0,003 3.5 x 10 h/W h/W h/W h/W h/W = 0,010 = 0,007 = 0,005 = 0,004 = 0,003 -1 -2 -3 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.18 Đáp ứng độ võng w Hình 3.21 Đáp ứng ứng suất y Bảng 3.4 Giá trị lớn tham số khảo sát theo tỉ số h/W Tỉ số h/W w max 10-3 [m] max 10-2 [m/s] w max [m/s2] w y 10 [N/m ] -6 y 10 0,010 0,007 0,005 0,004 0,003 2,463 18,598 27,913 15,030 5,647 4,332 21,992 36,126 15,731 7,466 6,923 31,970 36,832 16,991 8,317 9,025 37,847 44,939 16,386 9,780 Tăng theo thời gian Nhận xét: Tỷ số chiều dày bề rộng (h/W) nhạy cảm với đáp ứng động tấm: Khi tỷ số h/W thay đổi từ 0,007 đến 0,004 giá trị chuyển vị lớn tăng lên 2,08 lần, giá trị ứng suất tăng không đáng kể (4,2%) Với giá trị h/W = 0,003, chuyển vị tăng theo thời gian có xu hướng khơng dừng, điều dẫn đến khả bị ổn định Điều làm sở cho tính tốn điều khiển dao động ổn định kết cấu chịu tác dụng lực khí động 3.3.4 Ảnh hưởng điều kiện liên kết Để xem xét quy luật học đáp ứng động lực học để từ cho phép đánh giá phù hợp thuật tốn độ tin cậy chương trình tính, phần tác giả giải tốn với trường hợp liên kết: 13 Ngàm cạnh ngắn, ngàm cạnh ngắn ngàm cạnh (1 cạnh ngắn tự do) Điểm xuất kết so sánh điểm C mặt điểm B cạnh bị ngàm trường hợp ngàm cạnh ngắn -3 x 10 x 10 Ngam canh Ngam canh Ngam canh B -1 -2 -3 Ngam canh Ngam canh Ngam canh 1.5 Ung suat Xicmay [N/m2] Chuyen vi phuong dung tai C [m] 0.5 -0.5 -1 -1.5 -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.23 Đáp ứng độ võng w Hình 3.26 Đáp ứng ứng suất y Bảng 3.5 Giá trị lớn tham số khảo sát theo điều kiện biên Điều kiện biên w Cmax 10-3 [m] Cmax 10-2[m/s] w Ngàm cạnh 2,386 Ngàm cạnh 0,106 Ngàm cạnh 0,089 14,693 5,915 3,413 Cmax [m/s2] w 19,739 11,700 6,796 By 104 [N/m2] 16,992 2,615 2,830 By 10-6 8,318 2,725 2,665 Nhận xét: Khi điều kiện liên kết thay đổi, đáp ứng động lực học thay đổi: Với trường hợp liên kết khảo sát, trường hợp ngàm cạnh đối xứng cho thấy ứng suất cạnh bị ngàm nhỏ nhất, chuyển vị trường hợp ngàm cạnh (1 cạnh tự do) nhỏ Điều hoàn toàn phù hợp với quy luật học 3.3.5 Ảnh hưởng góc đặt cốt Khảo sát tốn với góc đặt cốt [p/-0/0/0/-0/p], góc cốt biến thiên từ 00 đến 900 Kết biến thiên giá trị lớn chuyển vị ứng suất điểm tính thể bảng 3.6 hình 3.28, 3.31 Bảng 3.6 Biến thiên giá trị lớn đại lượng đến góc Chuyển vị Vận tốc Gia tốc Ứng suất Biến dạng o -6 [ ] -3 y 10 10-2[m/s] w [m/s2] y 10 [N/m ] w 10 [m] w 30 6,942 32,419 36,194 7,506 18,948 45 6,923 31,970 36,832 16,991 8,317 60 6,916 31,759 37,065 19,659 4,596 75 6,914 31,674 37,130 20,456 2,538 90 6,913 31,646 37,150 21,118 0,308 14 -3 x 10 6.95 2.2 Ung suat XicmaYm ax [N/m2] 6.945 6.94 6.935 1.8 1.6 B Chuyen vi thang dung lon nhat [m] x 10 6.93 6.925 6.92 1.4 1.2 0.8 6.915 6.91 30 40 50 60 70 Goc dat cot teta [o] 80 0.6 30 90 40 50 60 70 Goc dat cot teta [o] 80 90 Hình 3.28 Đáp ứng độ võng max Hình 3.31 Đáp ứng ứng suất max Nhận xét: Góc đặt cốt có ảnh hưởng rõ rệt đến độ cứng tấm, điều thông qua thay đổi giá trị chuyển vị, biến dạng theo đại lượng này.Qua khảo sát trên, cho thấy với góc đặt cốt [450 900] hợp lý 3.3.6 Ảnh hưởng góc tới dịng khí Giải tốn với trường hợp góc tới dịng khí thay đổi từ 00 đến 600, kết đáp ứng chuyển vị ứng suất thể đồ thị hình 3.33, 3.36 bảng 3.7 -3 x 10 x 10 Alpha = 0o Alpha = 10o Alpha = 30o Alpha = 60o Alpha = 20o -2 -0.5 -1 3.5 Alpha = 600 0.5 -1.5 1.5 2.5 Thoi gian t[s] Alpha = 300 -6 Alpha = 200 1.5 -4 0.5 Alpha = 100 2 -8 Alpha = 00 2.5 Ung suat XicmaB [N/m2] y Chuyen vi phuong dung tai A [m] 10 -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.33 Đáp ứng độ võng Hình 3.39 Đáp ứng ứng suất By Bảng 3.7 Giá trị lớn tham số khảo sát theo góc tới α Chuyển vị o α[ ] w 10-3[m] 10 20 30 40 50 60 1,342 4,329 5,712 6,923 4,711 2,685 2,930 Vận tốc 10w [m/s] 16,786 23,207 17,617 31,970 15,273 19,689 20,155 Gia tốc [m/s2] w 25,087 31,553 25,127 36,832 22,904 28,214 28,706 Ứng suất y 10 [N/m2] 3,941 11,054 13,872 16,991 11,756 7,263 7,784 Biến dạng y 10 -6 2,506 5,660 6,616 8,317 5,554 3,950 4,214 15 Nhận xét: Khi góc tới thay đổi từ 00 đến 600, đáp ứng động thay đổi phức tạp, thay đổi không theo quy luật rõ nét Trong khoảng khảo sát, góc tới thay đổi từ 100 đến 400 đáp ứng chuyển vị, ứng suất biến dạng đạt giá trị lớn 3.3.7 Ảnh hưởng kích thước gân Thay đổi tỷ số chiều cao/bề rộng (hg/bg) mặt cắt ngang gân cách giữ nguyên bg, thay đổi hg, kết đáp ứng chuyển vị ứng suất thể bảng 3.8 đồ thị hình 3.43 3.47 Bảng 3.8 Giá trị cực trị chuyển vị điểm tính thay đổi hg hg [m] 0,003 0,004 0,005 0,006 0,007 Chuyển vị Vận tốc -3 10-2[m/s] w 10 [m] w 7,419 34,231 7,190 31,609 6,923 31,970 6,700 30,985 6,188 27,540 Gia tốc Ứng suất Biến dạng -6 y 10 [m/s ] y 10 [N/m ] w 46,274 16,253 8,467 41,995 17,205 8,326 36,832 16,991 8,317 31,889 15,976 8,299 29,467 16,293 8,221 -3 x 10 x 10 1.5 -2 0.5 -0.5 -4 -1 -6 -1.5 -8 hg = 0,003 hg = 0,004 hg = 0,005 hg = 0,006 hg = 0,007 B Chuyen vi phuong dung tai A [m] 2.5 hg = 0,003 hg = 0,004 hg = 0,005 hg = 0,006 hg = 0,007 Ung suat Xicmay [N/m2] 10 0.5 1.5 2.5 3.5 Th oi g ian t[s] Hình 3.43 Đáp ứng chuyển vị wA -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.49 Đáp ứng ứng suất By Nhận xét: Khi tỷ số hg/bg tăng giá trị lớn chuyển vị, vận tốc gia tốc điểm tính giảm, ứng suất biến dạng điểm tính thay đổi không rõ quy luật Với thông số khảo sát tốn, nhận thấy chiều cao gân từ mm đến mm phù hợp 3.3.8 Ảnh hưởng điện áp V áp đặt lên lớp áp điện Tác giả tiến hành khảo sát toán với điện áp đặt lên áp điện V thay đổi từ 50V đến 200V Kết biến thiên chuyển dịch lớn điểm A biến thiên ứng suất pháp lớn điểm B thể bảng 3.9 đồ thị hình 3.53, 3.56 Bảng 3.9 Giá trị cực trị chuyển vị điểm tính thay đổi V 16 Chuyển vị Vận tốc V -2 (Vol) w.10 [m] w 10-2 [m/s] 50 6,923 31,970 100 5,776 19,936 150 4,928 17,666 200 4,103 15,897 Gia tốc [m/s2] w 36,832 15,003 14,673 14,053 Ứng suất y 10 [N/m ] 16,991 15,299 14,988 14,498 -3 10 x 10 2.5 V = 50V V = 100V V = 150V V = 200V B -2 0.5 -0.5 -4 -1 -6 -1.5 1.5 2.5 Thoi gian t[s] 3.5 V = 50V V = 100V V = 150V V = 200V 1.5 0.5 x 10 Ung suat Xicmay [N/m2] Chuyen vi phuong dung tai A [m] -8 Biến dạng -6 y 10 8,317 6,758 5,891 5,125 Hình 3.53 Đáp ứng chuyển vị wA -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.56 Đáp ứng ứng suất By Nhận xét: Điện áp áp đặt lên lớp áp điện có tác dụng làm giảm dao động tấm, cụ thể điện áp tăng biên độ dao động giảm Điều cho thấy khả điều khiển, tăng độ cứng composite áp điện có gân gia cường vật liệu áp điện 3.3.9 Ảnh hưởng kích thước áp điện Xem xét toán với việc thay áp điện phía phía composite xét miếng áp điện giống nhau, bố trí mặt mặt đối xứng qua mặt phẳng trung bình tấm, chiều dày miếng lớp áp điện composite áp điện xét, miếng áp điện có chiều rộng chiều rộng composite, chiều dài 2/3 chiều dài composite Kết quả, biến thiên chuyển vị lớn điểm A ứng suất pháp lớn điểm B thể bảng 3.10 đồ thị hình 3.58 3.61 Bảng 3.10 Giá trị lớn đại lượng thay đổi kích thước áp điện Trường hợp Tấm AD Miếng AD Chuyển vị Vận tốc Gia tốc w.10-3 [m] 6,923 7,379 10-2 [m/s] w [m/s2] w 31,970 36,169 36,832 48,777 Ứng suất Biến dạng -6 y 10 [N/m ] y 10 16,991 8,317 18,862 9,373 17 -3 x 10 Tam AD Mieng AD -2 -4 -6 -8 Tam AD Mieng AD 1.5 Ung suat XicmaB [N/m2] y Chuyen vi phuong dung tai A [m] x 10 0.5 -0.5 -1 -1.5 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.58 Đáp ứng chuyển vị wA -2 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.61 Đáp ứng ứng suất By Nhận xét: Kích thước lớp áp điện có tác dụng tăng “cứng” cho tấm, điều thể qua đáp ứng chuyển vị, vận tốc, gia tốc ứng suất điểm tính, giá trị giảm diện tích lớp áp điện tăng 3.3.10 Ảnh hưởng tính chất cản Xét trường hợp: - Trường hợp (TH1): Chỉ kể đến cản áp điện, không kể đến cản kết cấu ( = 0), gọi cản áp điện; - Trường hợp (TH2): Chỉ kể đến cản kết cấu, không kể đến cản áp điện (Gv = 0), gọi cản kết cấu; - Trường hợp (TH3): Kể đến tổng cản (cản kết cấu, cản áp điện) – toán xuất phát mục 3.2 trình bày trên, gọi tổng cản; Kết biến thiên chuyển vị, vận tốc, gia tốc điểm tính (điểm A B) thuộc thể bảng 3.11 hình 3.63 3.66 Bảng 3.11 Giá trị cực trị chuyển vị điểm tính thay đổi thay đổi thành phần lực cản Tính chất cản TH1 TH2 TH3 Chuyển vị w 10 [m] 8,068 6,923 Ứng suất y 10 10w w [m/s] [m/s2] [N/m2] Biên độ tăng dần theo thời gian 44,514 61,654 20,605 31,970 36,832 16,991 Vận tốc Gia tốc Biến dạng y 10 -6 10,666 8,317 18 0.1 TH1 TH2 TH3 0.5 0.06 Ung suat XicmaB [N/m2] y Chuyen vi phuong dung tai A [m] 0.08 x 10 0.04 0.02 -0.02 -0.5 -1 TH1 TH2 TH3 -0.04 -0.06 0.5 1.5 2.5 Thoi gian t[s] 3.5 -1.5 0.5 1.5 2.5 Thoi gian t[s] 3.5 Hình 3.63 Đáp ứng chuyển vị wA Hình 3.66 Đáp ứng ứng suất By Nhận xét: Với toán khảo sát, xét đến cản kết cấu đặc biệt phối hợp cản áp điện cản kết cấu (tổng cản) dao động bình ổn, quan niệm kể đến cản áp điện dao động với biên độ tăng dần theo thời gian, có tượng ổn định 3.3.11 Miền ổn định tham số thay đổi Miền ổn định trường hợp tham số thay đổi tập hợp tất điểm mặt phẳng toạ độ mà trục thể biến thiên tham số lựa chọn Giải toán thay đổi giá trị tham số này, cịn tham số khơng thay đổi, nhận điểm mặt phẳng toạ độ tô màu đậm (nếu hệ ổn định), để trống (nếu hệ ổn định) Khảo sát toán xuất phát với trường hợp: TH1: Góc đặt cốt điện áp áp đặt V thay đổi, biến thiên từ đến 1800 V biến thiên từ 0V đến 150V; TH2: Góc tới dịng khí vận tốc dịng khí U thay đổi, biến thiên từ 00 đến 900 U biến thiên từ 0m/s đến 50m/s Kết miền ổn định cho trường hợp thể hình 3.68 3.69 Hình 3.68 Miền ổn định ổn Hình 3.69 Miền ổn định ổn định theo V định theo U Nhận xét: Từ miền ổn định cho thấy, cho phép xác định giá trị hợp lý góc đặt cốt cho lớp composite điện áp đặt lên lớp áp 19 điện lựa chọn cặp giá trị góc vận tốc gió hợp lý để điều khiển khả ổn định cách chủ động 3.4 Kết luận chương - Khảo sát số nhiều lớp toán composite áp điện có gân gia cường chịu tác dụng lực khí động, thơng số vật liệu, kết cấu, liên kết thay đổi Kết khảo sát số cho thấy ảnh hưởng đáng kể tốc độ dịng khí, tính chất áp điện, tỷ lệ bề rộng chiều dày tấm, đến khả dao động ổn định Qua xác định tốc độ gió tới hạn để có giải pháp kết cấu cho phép điều khiển dao động ổn định chịu tác dụng lực khí động Ngồi ra, xây dựng miền ổn định cho với hai cặp thông số: -V; U-α để lựa chọn thông số tối ưu điều khiển dao động cách chủ động tính tốn kết cấu cho phù hợp với lực khí động điều kiện cụ thể - Các số liệu, nhận xét có tác dụng định hướng nghiên cứu đáp ứng động lực học, ổn định kết cấu composite áp điện có gân gia cường chịu tác dụng lực khí động áp dụng lĩnh vực kỹ thuật CHƯƠNG NGHIÊN CỨU THỰC NGHIỆM 4.1 Đặt vấn đề Trong chương tác giả luận án tiến hành nghiên cứu thực nghiệm nhằm xác định phản ứng động composite có miếng áp điện chịu tác dụng lực khí động Với kết thực nghiệm thu được, tác giả tiến hành so sánh với kết tính tốn lý thuyết chương trình tính lập chương nhằm góp phần khẳng định tính phù hợp thuật toán độ tin cậy chương trình lập 4.2 Mơ tả thí nghiệm 4.2.1 Thiết lập thí nghiệm Tấm composite có kích thước: 500x300x3,5 mm (axbxh), gia cường gân dọc (500x7mm) gân ngang (300x7mm) trực giao song song với cạnh Tấm gân có loại vật liệu composite Graphite/Epoxy T300/976 Trên có dán miếng áp điện PZT phía mặt khơng có gân Tấm liên kết ngàm cạnh ngắn khung giá đỡ chuyên dụng xoay chỉnh theo hướng gió tương ứng với góc Giá đỡ xem cứng tuyệt đối - Miếng áp điện loại PZT có kích thước 200900,1 mm 20 Hình 4.1 Thiết lập mơ hình thí nghiệm Thơng số vật liệu làm mơ hình thí nghiệm sau: - Tấm gân composite Graphite/Epoxy T300/976: E1 =178,2 GPa; E2= 9,78 Gpa; G12 = G13 = 7,25 GPa; G23 = 2,51 GPa; = 0,34; = 1600 kg/m3 - Tấm áp điện PZT G-1195N: E = 63,0 GPa; G = 24,2 GPa; = 0.3, = 7600 kg/m3; d31 = d32 = 2.54×10-10 m/V; p11 = p22 = p33 = 15×10-9 F/m 4.2.2 Thiết bị đo đáp ứng động Để xác định gia tốc biến dạng tấm, tác giả sử dụng cảm biến gia tốc piezo PZT liên kết trực tiếp lên tấm, vị trí cần đo (Hình 4.1) 4.2.3 Bộ phát tín hiệu chuẩn OG-422A Bộ phát tín hiệu chuẩn thiết bị phát điện áp (V) dạng hình sine cho phép thay đổi điện áp (V) với tần số điều chỉnh tùy chọn 4.2.4 Bộ huếch đại piezo tuyến tính EPA-104-230 Bộ khuếch đại dùng điều chỉnh điện áp vào, áp điện 4.2.5 Bộ thu nhận liệu kênh áp điện HnB75B Bộ thu nhận liệu kênh áp điện HnB75B cho phép thu nhận tín hiệu đo tương tự chuyển đổi sang tín hiệu số theo thời gian thực 4.2.6 Máy sóng (oscilloscope) Máy sóng ký hiệu GWinstek GDS-2104 có chức ghi lại đáp ứng gia tốc điểm đo lưu trữ liệu vào máy tính 4.2.7 Thiết bị tạo gió: Thiết bị tạo gió (hầm gió hở) sử dụng động 22 kw để tạo gió buồng thử có kích thước 1000x1000 mm Tốc độ thổi gió lớn đạt 15m/s 4.3 Phương pháp xác định gia tốc, biến dạng kết cấu 4.3.1 Đo gia tốc: Để đo gia tốc dao động vị trí thuộc tấm, tiến hành gắn đầu đo gia tốc cố định điểm cần đo, tùy theo phương gia tốc cần đo, đầu đo gia tốc gắn vào kết cấu cho phương trục đầu đo trùng với phương cần đo gia tốc Theo đó, tác dụng lực khí động, 21 dao động đồng thời đầu đo gia tốc dao động theo, tín hiệu đáp ứng gia tốc theo thời gian truyền khối thu, phân tích liệu đo Ứng với lần đo, tương ứng với đầu đo ta có số liệu thể đáp ứng gia tốc theo thời gian vị trí gắn đầu đo gia tốc 4.3.2 Đo biến dạng: Để có biến dạng điểm tấm, cần dán điện trở đo biến dạng vị trí đo Dưới tác dụng lực khí động, bị biến dạng làm cho PZT biến dạng theo, tín hiệu học chuyển thành tín hiệu số truyền khối thu, xử lý tín hiệu máy đo, kết có số liệu đáp ứng biến dạng theo thời gian điểm đo 4.4 Cơ sở phân tích xử lý kết thí nghiệm Việc thu thập, tích hợp số liệu (đáp ứng gia tốc, biến dạng theo thời gian) lần thí nghiệm mơđun lưu trữ, xử lý có thiết bị đo, số liệu thống kê cho 01 lần đo Để có kết đo đại diện, cần phải xử lý thống kê từ n lần đo, điều tác giả luận án thực sở lý thuyết xử lý số liệu xác suất thống kê 4.5 Thí nghiệm kết thí nghiệm Tiến hành thí nghiệm 05 lần đo (n = 5) để xác định số liệu đáp ứng gia tốc theo thời gian Sau xử lý số liệu theo phương pháp thống kê phần mềm DAS-100A chuyên dụng, ta có đáp ứng theo thời gian gia tốc điện áp vị trí đo Hình 4.11 Bố trí thiết bị đo thí nghiệm thực tế Kết đo: TH1 Tấm góc 0, vận tốc gió U = 5m/s 10m/s, điện áp kích Vin = 9.30V (đỉnh), tần số kích thích f = 6,944Hz, Kamp= 20, tần số lấy mẫu fsamp = 1000 Hz Điện áp cấp 186.0 V TH2 Tấm 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 22 TH3 Tấm lệch góc 450, 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 đo đáp ứng gia tốc theo thời gian lần đo với vận tốc gió U = 5m/s trường hợp thứ (TH1) thể hình 4.14 giá trị Max/Min đại lượng theo thời gian thể bảng 4.1 Trong luận án tác giả tiến hành đo lần cho trường hợp khác Hình 4.14 Đáp ứng gia tốc theo thời gian lần đo Bảng 4.1 Giá trị Max/Min theo thời gian đại lượng lần đo Kết so sánh đáp ứng gia tốc điểm đo với U = 5m/s, = 450 tính tốn lý thuyết tác giả (chương trình SSP-2018) thực nghiệm thể hình 4.15 Hình 4.15 Đáp ứng gia tốc điểm đo (U = 5m/s, = 450) Bảng 4.2 Giá trị lớn gia tốc (U = 5m/s góc gió khác nhau) Phương pháp Thí nghiệm SSP_2018 Thí nghiệm SSP_2018 Thí nghiệm SSP_2018 amax[m/s2] 0,050 0,0559 0,061 0,0675 0,072 0,0793 Sai khác[%] 450 11,8 22,50 10,7 00 10,1 23 Nhận xét: Kết nghiên cứu lý thuyết chương trình tính SSP_2018 tác giả lập kết thí nghiệm đáp ứng gia tốc điểm đo thuộc với góc hướng gió khác đồng dạng, với sai số lớn trường hợp 11,8% Điều cho thấy thuật tốn chương trình tính SSP_2018 tác giả lập chương có sở để tin cậy 4.6 Kết luận chương Qua việc thí nghiệm, tác giả rút số kết luận sau: - Đáp ứng gia tốc thực nghiệm có độ trễ so với tính tốn lý thuyết phù hợp với thực tế Tốc độ gió ảnh hưởng đến dao động tấm, góc dịng khí lớn dao động lại nhỏ so với góc nhỏ bị biến dạng trở nên “cứng” Ngồi ra, mơ hình thí nghiệm chưa giống hồn tồn so với mơ hình tính tốn lý thuyết nên có sai số đáp ứng gia tốc lý thuyết thực nghiệm (sai số lớn 11,8%), nhiên, sai số phạm vi chấp nhận - Kết nghiên cứu thực nghiệm tính tốn lý thuyết chương trình tính SSP_2018 mơ hình tương tự đồng dạng quy luật cho thấy độ tin cậy chương trình tính tác giả lập chương KẾT LUẬN VÀ KIẾN NGHỊ I Những đóng góp luận án: 1, Sử dụng phương pháp PTHH xây dựng thuật toán chương trình tính SSP_2018 phân tích phi tuyến động lực học composite áp điện có gân gia cường chịu tác dụng tải trọng khí động, xét đến tính chất cản tổng thể hệ Ngồi việc xác định đáp ứng phi tuyến động lực học, thuật tốn chương trình tính cịn cho phép xem khả ổn định tác dụng lực khí động Chương trình tính tác giả lập kiểm chứng cho thấy đảm bảo tin cậy Nội dung phần cơng bố cơng trình [1], [3], [4] tác giả luận án cộng 2, Khảo sát số nhiều lớp tốn với thơng số kết cấu, tải trọng, tính chất cản, liên kết, điện áp, đưa nhận xét, đánh giá định lượng có ý nghĩa khoa học thực tiễn lĩnh vực kỹ thuật Xây dựng miền ổn định động tấm, cho phép lựa chọn cặp thông số hợp lý để tăng khả ổn định giảm dao động cho Nội dung công bố cơng trình [2], [5] tác giả cộng 24 3, Bộ số liệu thực nghiệm nghiên cứu đáp ứng động lực học kết cấu composite áp điện có gân gia cường chịu tác dụng lực khí động góp phần khẳng định phù hợp thuật toán độ tin cậy chương trình tính tác giả lập Nội dung cơng bố cơng trình [6] tác giả cộng II Nhận xét kiến nghị: 1, Tính tốn phi tuyến động lực học cho phép xem xét toán ổn định kết cấu composite gia cường có tham gia lớp miếng áp điện chịu tải trọng động, có kể đến tính chất cản áp điện, cản kết cấu vấn đề khó phức tạp, qua kết luận án khẳng định tính tin cậy phương pháp chương trình tính Đây điều kiện tốt cho nghiên cứu kết cấu phức tạp như: composite áp điện lượn sóng, kết cấu vỏ composite áp điện có gân gia cường phục vụ cho lĩnh vực kỹ thuật 2, Ảnh hưởng tính chất cản, tính chất áp điện, vị trí miếng áp điện, kích thước gân, góc đặt cốt, đến làm việc composite áp điện có gân gia cường lớn, hầu hết đại lượng chuyển vị, vận tốc, gia tốc không thay đổi cách tuyến tính Cản áp điện kết hợp với cản kết cấu làm giảm dao động cách đáng kể nên tính tốn kết cấu dạng cần phải kể đến loại cản Miền ổn định composite áp điện có gân gia cường chịu tải trọng động có ý nghĩa thực tiễn, qua đánh giá khả ổn định hay không ổn định kết cấu 4, Nội dung nghiên cứu luận án phát triển theo hướng: - Phân tích dao động ổn định lượn sóng, vỏ composite áp điện có gân gia cường chịu tác dụng đồng thời tải trọng nhiệt độ - Phân tích dao động ổn định composite áp điện lượn sóng, vỏ composite áp điện có gân gia cường chịu tác dụng lực khí động hay lực khí động nhiệt độ - Nghiên cứu thực nghiệm kết cấu tấm, vỏ composite áp điện có biện pháp gia cường như: gân, lượn sóng, gấp nếp,v.v ... cho composite áp điện có gân gia cường chịu tác dụng tải trọng khí động theo thời gian Nghiên cứu thực nghiệm, xác định phản ứng động composite áp điện có gân gia cường chịu tải trọng khí động. .. cấu composite áp điện làm sở giải tốn phân tích động lực học composite áp điện có gân gia cường Xây dựng thuật tốn PTHH chương trình máy tính phân tích dao động composite áp điện có gân gia cường. .. vỏ composite áp điện có gân gia cường chịu tác dụng đồng thời tải trọng nhiệt độ - Phân tích dao động ổn định composite áp điện lượn sóng, vỏ composite áp điện có gân gia cường chịu tác dụng lực