Phân tích dao động của tấm có cơ tính biến thiên chịu tác dụng của lực khí động và nhiệt độ

Chia sẻ: Thi Thi | Ngày: | Loại File: PDF | Số trang:6

Thêm vào BST

Báo xấu

72
lượt xem 2
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài báo trình bày kết quả phân tích dao động của tấm có cơ tính biến thiên chịu tác dụng của lực khí động và nhiệt độ. Theo đó, hệ phương trình vi phân mô tả dao động của tấm được các tác giả giải trên cơ sở tích phân trực tiếp Newmark. Chương trình tính cụ thể hóa thuật toán, phân tích bài toán được các tác giả viết trong môi trường Matlab. Các kết quả nghiên cứu góp phần phục vụ việc thiết kế và sửa chữa lớp vỏ của thân và cánh của các thiết bị bay.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Phân tích dao động của tấm có cơ tính biến thiên chịu tác dụng của lực khí động và nhiệt độ

Lê Thúc Định và Đtg Tạp chí KHOA HỌC & CÔNG NGHỆ 122(08): 15 - 20 PHÂN TÍCH DAO ĐỘNG CỦA TẤM CÓ CƠ TÍNH BIẾN THIÊN CHỊU TÁC DỤNG CỦA LỰC KHÍ ĐỘNG VÀ NHIỆT ĐỘ Lê Thúc Định1*, Vũ Quốc Trụ1, Trần Thị Hương2 1Học viện Kỹ thuật Quân sự, 2Trường Cao đẳng Kỹ thuật Lý Tự Trọng TÓM TẮT Bài báo trình bày kết quả phân tích dao động của tấm có cơ tính biến thiên chịu tác dụng của lực khí động và nhiệt độ. Theo đó, hệ phương trình vi phân mô tả dao động của tấm được các tác giả giải trên cơ sở tích phân trực tiếp Newmark. Chương trình tính cụ thể hóa thuật toán, phân tích bài toán được các tác giả viết trong môi trường Matlab. Các kết quả nghiên cứu góp phần phục vụ việc thiết kế và sửa chữa lớp vỏ của thân và cánh của các thiết bị bay. Từ khóa: vật liệu có cơ tính biến thiên, gốm, kim loại, dao động, lực khí động, nhiệt độ ĐẶT VẤN ĐỀ* Panel flutter là hiện tượng dao động tự kích gây mất ổn định cục bộ của tấm hoặc vỏ mỏng đàn hồi chịu tác dụng của lực khí động. Đối với các thiết bị bay có tốc độ lớn, việc nghiên cứu hiện tượng panel flutter là rất cần thiết cho việc thiết kế, chế tạo lớp vỏ của chúng. NỘI DUNG BÀI TOÁN Giới thiệu vật liệu có cơ tính biến thiên Vật liệu có cơ tính biến thiên (FGM) thường gặp là loại 2 thành phần, nó là hỗn hợp của gốm (ceramic) và kim loại (metal) với tỷ lệ thể tích của các thành phần biến đổi trơn, liên tục theo chiều dày thành kết cấu và là hàm lũy thừa của biến chiều dày z. k  z 1 Vc (z)     h 2 Vm (z)  1  Vc (z), (0  k  ) (1) trong đó: k- chỉ số mũ tỷ lệ thể tích; Vc (z) , Vm (z) - tỉ lệ thể tích của thành phần gốm và kim loại tương ứng; h - chiều dày thành kết cấu; z là trục tọa độ hướng theo chiều dày của tấm. Tính chất hiệu dụng của vật liệu được xác định theo biểu thức sau: k  z 1 Pe   Pc  Pm      Pm h 2 (2) với Pe, Pc, Pm là tính chất hiệu dụng của vật liệu FGM, gốm và kim loại tương ứng. Sự phụ thuộc vào nhiệt độ của tính chất các vật liệu thành phần (gốm, kim loại) [4] [5] [6]: * Tel: 0982 140560, Email: ledinhvhp@gmail.com P(T)  P0  P1T 1  1  P1T  P2T 2  P3T 3  (3) trong đó: P0, P-1, P1, P2, và P3 là các hệ số của nhiệt độ; T - nhiệt độ (K). Do đó, tính chất hiệu dụng của vật liệu FGM theo tọa độ và nhiệt độ được xác định như sau:  2z  h  Pe  T,z   Pc  T   Pm  T     Pm  T   2h  k (4) Phân bố nhiệt độ theo chiều dày của tấm: T  z   Tm   Tc  Tm    z  trong đó   z  là hàm truyền nhiệt được xác định [3] [6]: k 1  2z  h   K cm  2z  h        2h   k  1 K m  2h     2k 1 2 K 2z  h     cm      2k  1 K 2  2h   m   3k  1  1 K 3cm  2z  h    z       3  C   3k  1 K m  2h     4k 1 4 K cm  2z  h       4   4k  1 K  2h   m   5k  1   K 5cm  2z  h     5    5k  1 K m  2h   (5) với C 1  2 K cm K cm    k  1 K m  2k  1 K m2 4 K 3cm K cm K 5cm   3 4  3k  1 K m  4k  1 K m 5k  1 K 5m (6) K cm  K c  K m Mô hình bài toán và các giả thiết Xét tấm vật liệu có cơ tính biến thiên hình chữ nhật, kích thước như hình vẽ (hình 1), 15 Lê Thúc Định và Đtg Tạp chí KHOA HỌC & CÔNG NGHỆ mặt trên của tấm chịu tác dụng của dòng khí vượt âm với vận tốc V hướng theo trục x (trùng với một cạnh tấm), mặt dưới chịu áp suất không đổi bằng áp suất tĩnh của dòng khí tự do. Bề mặt giàu gốm có nhiệt độ lớn hơn bề mặt giàu kim loại. Kết cấu được xét là tấm mỏng đàn hồi. Bỏ qua biến dạng cắt ngang 122(08): 15 - 20 của phần tử: u, v, w T   N q e  (9) trong đó: u,v - chuyển vị màng; w - chuyển vị uốn (độ võng); [N] - ma trận hàm dạng; {q e} véc tơ chuyển vị nút của phần tử. y w4   xz   yz  0 z v4 x4 u4 y4 w1 y1 Hình 1. Mô hình bài toán w3 v3 y3 x3 u3 w2 v1 x1 u1 v2 x2 u2 y2 x Hình 2. Phần tử tấm phẳng hình chữ nhật 4 nút Lực khí động Quan hệ biến dạng - chuyển vị Với mô hình của bài toán, tác giả lựa chọn mô hình lực khí động theo lý thuyết piston bậc nhất 0: Véc tơ các thành phần biến dạng và độ cong của phần tử:  g D w D w  Pa    a 110   110  4  a  t a 3 x   0 (7) với 1  D 2  v2 q  a ;   M 2  1; 0   1104  ; 2  ha  Ca M  M 2  2   2  a a a v  M2  2  2qa 3 ;g  ; a 2 3 h0 D110  1 h 2 (8)  2w   u   2     x   x  2  v      w     m  ; m     ; k    2  (10)  k   y   y   u v   2w     2   xy   y x  trong đó:  m  , k tương ứng là véc tơ biến dạng màng, véc tơ độ cong. trong đó: Pa - lực khí động, v - vận tốc dòng khí, M - số Mach; q - áp suất khí động, a khối lượng riêng của dòng khí;  - áp suất khí động không thứ nguyên ; D110 - độ cứng trụ của tấm; Ca - hệ số cản khí động; g a - hệ số cản khí động không thứ nguyên; 0 - tần số quy ước;  - khối lượng riêng của tấm; h chiều dày của tấm; a - chiều dài của tấm; w chuyển vị của tấm theo phương z. Biểu diễn biến dạng theo chuyển vị nút:  N    A   M      B   B m   NT   D  k  MT  Các phương trình cơ bản trong đó: N  Nx , Ny , Nxy  - lực màng; {}  [B0 ]{qe} (11) với [B0 ] là ma trận tính biến dạng. Quan hệ ứng suất - biến dạng Quan hệ giữa ứng suất và biến dạng có xét đến ảnh hưởng của nhiệt độ: (12) T Quan hệ chuyển vị - chuyển vị nút M  Mx ,My ,Mxy - mô men uốn và xoắn; Chọn loại phần tử phẳng hình chữ nhật 4 nút, mỗi nút có 5 bậc tự do (hình 2). Chuyển vị tại một điểm bất kỳ trên mặt trung bình của phần tử được biểu diễn thông qua chuyển vị nút [A], [B], [D] là ma trận độ cứng màng, ma trận độ cứng tương tác màng-uốn-xoắn và ma trận độ cứng uốn tương ứng;  N T  , M T  lực màng và mô men uốn, xoắn do biến dạng nhiệt. 16 T Lê Thúc Định và Đtg h Tạp chí KHOA HỌC & CÔNG NGHỆ màng-uốn:  A, B, D    E 1, z, z 2 dz 2 h [K emu ]   [B0m ]T [B][B0u ]dS 2   z   N T  ,M T    1, z  E   z  T  z dz h  0  2   h S S - Ma trận cản khí động của phần tử: (13) 0  Aei   0  chuẩn  T0  300K  . Thiết lập phương trình chuyển động Áp dụng nguyên lý công khả dĩ đối với các lực tác dụng lên phần tử (lực đàn hồi, lực nhiệt, lực khí động và lực quán tính), ta có: (14) W=Wint  Wext  0 Công khả dĩ của nội lực:  Wint    m  N  k M dS  T T  q e   K e q e   q e  PTe  T T  w  hw+Pa   u  hu    Wext    dS    v  hv   S (16)  M q   A q    A q  ed e ei e Thay (15), (16) vào (14) và biến đổi ta được phương trình dao động của phần tử tấm dưới tác dụng của lực khí động và nhiệt độ như sau:  Me q e      A ei   ga q e   PTe   Aed q e     0   Ke   (17) (18) [K em ],[K eu ] - ma trận độ cứng của trạng thái màng, ma trận độ cứng của trạng thái uốn: [K em ]   [B0m ]T [A][B0m ]dS S [K eu ]   [B0u ]T [D][B0u ]dS  M   M e    0em  D11 T  N  N u  dS 3  u a x S ;  A eiu    (22) 0  M  eu  (23) với  Mem    h  Nm   Nm  dS; Meu    h  Nu   Nu  dS T T S S - Véc tơ tải trọng do thay đổi nhiệt độ: {P } {PTe }   Tem  (20x1) {PTeu }  (24) {PTem}   [B0m ]T [NT ]dS; {PTeu }   [B0u ]T [MT ]dS S (19) S [K emu ],[K eum ] - ma trận độ cứng tương tác (12x1) S - q e , q e : Véc tơ vận tốc, véc tơ gia tốc nút của phần tử. Thực hiện ghép nối các ma trận và véc tơ phần tử  M e  ,  K e  ,  A ed  ,  A ei  ,  PTe  ,  q e  ,  q e  ,  q e  ta được các ma trận và véc tơ tổng thể tương ứng của toàn bộ tấm  M  ,  K  ,  A d  ,  Ai  ,  PT  ,  q  ,  q  ,  q  . Khi đó, phương trình dao động của tấm FGM chịu tác dụng của tải trọng khí động và nhiệt độ có dạng: Mq  trong đó: - Ma trận độ cứng tuyến tính của phần tử:  [K ] [K ] Ke   [K em ] [Kemu]  eu   eum 0   Aeiu  - Ma trận khối lượng của phần tử [Me]: (8x1) Công khả dĩ của các ngoại lực lực (lực khí động và lực quán tính): e (21) với (15) S e 0  D110 T ; A    N   Nu  dS  Aedu  edu S a 4 u - Ma trận độ cứng do ảnh hưởng khí động của phần tử: k   q e  0  Aed   0   2z  h  E  T, z    E c  T   E m  T      Em T   2h  T  z   T  z   T0 , T0 là nhiệt độ T (20) [K eum ]   [B0u ]T [B][B0m ]dS  [K emu ]T 2   1  0  E  T, z    E 0   1 1  2  1   0 0  2    122(08): 15 - 20 ga Ad q    Ai   K q  PT  0 (25) Phương pháp giải bài toán Viết lại phương trình (25) ở dạng sau: 0  qm   ga 0 0  qm   Mu  qu   0 0 Adu  qu   (26)  0 0    K m   K mu   q m  PTm            0  Aiu     K um   K u    q u   PTu      M m     0 Khai triển (26) ta được hai phương trình sau: (27)  M m q m    K m q m    K mu q u   PTm  17 Lê Thúc Định và Đtg  M u q u   Tạp chí KHOA HỌC & CÔNG NGHỆ ga  Adu q u      Aiu   K u q u   0   K um q m   PTu  (28) Thay (29) vào (28) và biến đổi ta nhận được:    Aiu    K u    g  q u   Mu qu   a  Adu qu    1  0 (30)    K um  K m   K mu   1  PTu    K um  K m  PTm  Khi xét đến ảnh hưởng của cản kết cấu, phương trình (30) trở thành:  M u q u   Ckc   Ckd q u   1 (31)     A iu    K u    K um  K m   K mu q u   1 Các thông số của dòng khí: a = 1.225 kg/m3, a  = 340.3 m/s. Nhiệt độ mặt giàu gốm Tc Theo Vonmir, do quán tính màng rất nhỏ nên có thể bỏ qua. Do vậy, từ (27) ta rút ra được: 1 (29) qm   K m  PTm   K mu q u   PTu    K um  K m  122(08): 15 - 20 =320K, nhiệt độ mặt giàu kim loại Tm = 300K. Tấm được chia lưới đều 8×8, nhận được 64 phần tử và 81 nút. Tỉ số cản kết cấu   0,01 . Bước tích phân t = 1/100000s. Giải bài toán dao động riêng Giải bài toán dao động riêng ta thu được 2 dạng riêng đầu tiên (hình 3, hình 4) tương ứng với 2 tần số riêng đầu tiên: 1 = 1,23.103 Hz, 2 = 2,15.103 Hz. PTm  trong đó: Ckd  - ma trận cản khí động; Ckc  - ma trận cản kết cấu được xác định như sau: Ckc     M u     K u  với   (32) 2 ;   12 ; 1, 2 1  2 Hình 3. Dạng riêng thứ nhất ( 1 = 1,23.103 Hz) là hai tần số dao động riêng đầu tiên;  là tỉ số cản kết cấu. Để giải phương trình (31) ta sử dụng phương pháp tích phân trực tiếp Newmark và lập trình bằng ngôn ngữ Matlab. TÍNH TOÁN SỐ Xét tấm có cơ tính biến thiên 2 thành phần là gốm (Si3N4) và thép không gỉ (SUS304). Tấm hình chữ nhật kích thước (0.40×0.32×0.002)m, được ngàm cả 4 cạnh. Kết quả phân tích dao động của tấm Hình 4. Dạng riêng thứ hai ( 2 = 1,23.103 Hz) Hình 5. Dao động tắt dần của tấm với  = 820 (tại nút 41) 18 Lê Thúc Định và Đtg Tạp chí KHOA HỌC & CÔNG NGHỆ 122(08): 15 - 20 Hình 6. Dao động tuần hoàn của tấm với  = 879,5 (tại nút 41) Hình 7. Dao động với biên độ tăng dần của tấm với  = 882 (tại nút 41) * Nhận xét: - Khi  nhỏ (hình 5) lực khí động chưa đủ lớn nên dao động của tấm là dao động tắt dần. - Khi  tăng dần tới lân cận giá trị tới hạn (hình 6) tấm có xu hướng dao động tuần hoàn do xác lập trạng thái cân bằng động lực giữa lực khí động với đáp ứng tuyến tính của tấm. - Khi  vượt quá giá trị tới hạn (hình 7), dao động của tấm có biên độ tăng dần, tấm có khả năng mất ổn định. Khảo sát đáp ứng chuyển vị ngang lớn nhất (wmax/h) theo áp suất khí động không thứ nguyên và nhiệt độ * Nhận xét: Khi nhiệt độ tăng làm cho độ cứng uốn của tấm có cơ tính biên thiên giảm. Do đó, làm tăng giá trị chuyển vị ngang lớn nhất (wmax/h). Khi lực khí động càng lớn thì ảnh hưởng của nhiệt độ đến chuyển vị ngang tăng càng mạnh. KẾT LUẬN Bài báo đã đạt được các kết quả như sau: - Thiết lập và giải hệ phương trình vi phân mô tả dao động của kết cấu tấm FGM chịu tác dụng của lực khí động và nhiệt độ. - Xây dựng được chương trình tính trong môi trường Matlab. - Áp dụng tính toán số để xác định đáp ứng động của kết cấu tấm có cơ tính biến thiên cho 3 trường hợp khác nhau của lực khí động. Kết quả cho thấy tính chất dao động của tấm phụ thuộc vào lực khí động thông qua áp suất khí động không thứ nguyên, trên cơ sở đó có thể đánh giá khả năng ổn định của tấm. Hình 8. Đáp ứng chuyển vị ngang wmax/h theo  và Tc - Khảo sát ảnh hưởng của chuyển vị ngang lớn nhất theo lực khí động (thông qua áp suất 19