Đang tải... (xem toàn văn)
Bài viết trình bày các kết quả nghiên cứu mới về việc xác định vết nứt trong các kết cấu hệ dầm bằng vật liệu có cơ tính biến thiên (FGM) dựa trên phân tích wavelet dừng (SWT) các dạng dao động riêng có kể đến ảnh hưởng của nhiễu trắng Gausian.
Tạp chí Khoa học Cơng nghệ Xây dựng NUCE 2018 12 (7): 20–33 XÁC ĐỊNH VỊ TRÍ VẾT NỨT TRONG DẦM FGM BẰNG PHÂN TÍCH WAVELET DỪNG CÁC DẠNG DAO ĐỘNG RIÊNG Ngô Trọng Đứca , Trần Văn Liêna,∗, Nguyễn Thị Hườnga a Khoa Xây dựng Dân dụng Công nghiệp, Trường Đại học Xây dựng, 55 đường Giải Phóng, quận Hai Bà Trưng, Hà Nội, Việt Nam Nhận ngày 23/10/2018, Sửa xong 26/11/2018, Chấp nhận đăng 29/11/2018 Tóm tắt Bài báo trình bày kết nghiên cứu việc xác định vết nứt kết cấu hệ dầm vật liệu có tính biến thiên (FGM) dựa phân tích wavelet dừng (SWT) dạng dao động riêng có kể đến ảnh hưởng nhiễu trắng Gausian Các dạng dao động riêng xác định từ mơ hình phần tử FGM chịu kéo, nén uốn có nhiều vết nứt theo mơ hình lò xo phương pháp độ cứng động lực (DSM) Kết nghiên cứu cho thấy phương pháp đề xuất phương pháp khả thi hiệu Từ khoá: vết nứt; FGM; DSM; SWT; dạng dao động riêng DETERMINATION OF CRACK LOCATION IN BEAM USING USING STATIONARY WAVELET TRANSFORM OF MODE SHAPES Abstract This paper proposed new crack identification method on multiple cracked beams made of functionally graded material (FGM) by using stationary wavelet transform (SWT) of measured mode shapes This study also investigated the influence of Gaussian noise to SWT Cracks were modelled as equivalent springs and mode shapes are obtained from multiple cracked FGM beam element model The investigated results show that crack identification method by using SWT of mode shapes is efficient and realizable Keywords: crack; FGM; DSM; SWT; mode shapes https://doi.org/10.31814/stce.nuce2018-12(7)-03 c 2018 Trường Đại học Xây dựng (NUCE) Tổng quan Vật liệu có tính biến thiên (FGM) làm từ hỗn hợp gốm kim loại, có khả chịu nhiệt, chống ăn mòn gốm độ bền học làm giảm ứng suất nhiệt kim loại Đây loại vật liệu composite hệ mới, có tiềm ứng dụng rộng rãi nhiều lĩnh vực kỹ thuật môi trường nhiệt độ cao hàng khơng vũ trụ, lò phản ứng hạt nhân, máy phát điện hay công nghiệp ô tô Với ưu điểm đó, vật liệu FGM sử dụng cho kết cấu, chi tiết quan trọng thường xuyên làm việc môi trường khắc nghiệt Vết nứt hay hư hỏng xuất kết cấu FGM ảnh hưởng nghiêm trọng đến khả làm việc cấu kiện Do đó, vấn đề xác định vết nứt kết cấu FGM thực cần thiết thu hút ý nhiều nhà nghiên cứu nước Hầu hết nghiên cứu xác định vị trí vết nứt kết cấu sử dụng phương pháp thí nghiệm khơng phá hủy dựa đặc trưng động lực tần số, dạng dao động riêng, hàm phản ứng tần số, Các đặc trưng động lực học ∗ Tác giả Địa e-mail: lientv@nuce.edu.vn (Liên, T V.) 20 Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng dầm FGM có vết nứt xác định phương pháp giải tích [1–6], phương pháp bán giải tích Galerkin [7], phương pháp phần tử hữu hạn (FEM) [8–10], phương pháp độ cứng động lực (DSM) [11–14] Trong số phương pháp khác sử dụng để xác định vị trí số lượng vết nứt kết cấu [15], phương pháp phân tích wavelet có hiệu lớn, đặc biệt nhận dạng hư hỏng, vết nứt nhỏ kết cấu Liew Wang [16] chứng minh vị trí vết nứt dầm đơn giản xác định từ phân tích wavelet ứng xử dao động riêng đo dọc theo dầm thời điểm định Nghiên cứu tiếp tục phát triển Wang Deng [17] cho trường hợp phản ứng xung dầm với điều kiện biên khác Chang Chen [18] đề xuất phương pháp để xác định vị trí độ sâu dầm có nhiều vết nứt dựa phân tích khơng gian wavelet dạng dao động riêng Việc dự đoán độ sâu vết nứt đơn giản hóa cách sử dụng vị trí vết nứt dự đoán tần số dao động riêng Zhong Oyadiji [19] chứng minh phân tích wavelet dừng (SWT) công cụ hiệu cho việc chẩn đoán vết nứt dựa dạng riêng kết cấu dầm Hầu hết nghiên cứu công bố sử dụng phân tích wavelet để chẩn đốn vết nứt kết cấu dạng dầm, (SWT) cơng chẩnchỉđốn dựaqua dạng riêng kết cấu dầm nghiên cứu vềcụ hệhiệu kết cấu hệcho thanhviệc thường đượcvết thựcnứt hiệnchỉ thông FEM Hầu hếtĐối cácvới nghiên cứu cơng bố sử dụng phân tích wavelet để vết nứt kết vật liệu FGM, Banerjee cộng [10] đưa hai phương chẩn pháp đoán khác để nhận cấu dạng dạng vết dầm, nghiên cứu hệ kết cấu hệ chỉcác thựcvẽhiện thôngđường qua FEM nứtnhững dầm Timoshenko FGM Ở phương phápthường thứ nhất, tác giả lên Đối đồng với vật liệu FGM, Banerjee cộng [10] đưa hai phương pháp khác để nhận mức tần số dựa vào vị trí kích thước vết nứt, giao điểm đường đồng mức dạng vết nứt dầmlà Timoshenko FGM phương pháp nhất,Phương tác giảthứvẽhailên đường dạng dao động sở để dự đốn vị trí Ở chiều sâuthứ vết nứt pháp dựanhững mơ đồnghình mức số ứng dựa(RSM) vào vịhồi tríquy kích tối thước vếtbằng nứt,thuật giao điểm các(GA) đường đồngNam, mức mặttần phản ưu hóa tốn gen di truyền Tại Việt từngtác dạng dao động sở để dự đoán vị trí chiều sâu vết nứt Phương pháp thứ hai dựa giả Khiêm Huyên [20] đề xuất phương pháp chẩn đoán vết nứt dầm Timoshenko FGM mơ hình mặt phản ứng (RSM) hồi quy tối ưu hóa thuật tốn gen di truyền (GA) dựa biểu thức giải tích tần số dao động riêng Các phương pháp có điểm Tại Việtchung Nam,làtác Khiêm, Nguyễn Huyên [20] phương dựagiả trênNguyễn tần sốTiến dao động riêng đầu tiên,Ngọc nhiên thựcđãtế đề việcxuất xác định xác pháp chẩnnhiều đoántần vết nứt dầm Timoshenko FGM dựa biểu thức giải tích tần số dao động số dao động riêng khó khăn, vết nứt số vị trí dầm gần không riêng Các phương pháp có điểm chung dựa tần số dao động riêng đầu tiên, thay đổi tần số dao động riêng định, dẫn đến khó khăn việc sử dụng tần số để chẩn nhiên thực tế việc xác định xác nhiều tần số dao động riêng khó khăn, vết đốn vết nứt nứt số vị trí dầm gần không thay đổi tần số dao động riêng định, dẫn đến Nội dung báo sử dụng phân tích SWT để chẩn đốn vết nứt kết cấu khó khăn việc sử dụng tần số để chẩn đoán vết nứt dầm FGM dầm đơn giản, dầm liên tục nhiều nhịp, Đầu tiên, DSM sử dụng để xây dựng Nội mơ dung tử bàidầm báocónày vết sử nứt, dụng chẩn vết kết cấu hình phần nhiều từ phân tínhtích tốnSWT để xácđốn tần số nứt dạng daocác động dầmriêng FGMcủa đơn FGM giản,có dầm liênvếttục nhiều Đầudạng tiên,dao DSM kếtdầm cấu dầm nhiều nứt Tiếp nhịp, theo, từ độngđược riêngsử tìmdụng được,để tácxây giả dựng mơ hình phần tử dầm vết vị nứt, vàSWT dạng dao đề xuất phương phápcó để nhiều chẩn đốn trí từ vếtđó nứttính tốn kết cấu dầmchính FGM xác bằngcác cáchtần sử số dụng độngtrên riêng kết cấudạng dầmdao FGM nhiềuNgoài vết nứt Tiếp theo, từ dạng tìm được, tín hiệu đo độngcóriêng nghiên cứu xét đếndao ảnhđộng hưởngriêng nhiễu tác giả đề đến xuấthệmột phương đốn vị trí vết nứt kết cấu dầm FGM cách sử đo đạc số chi tiết củapháp phân để tíchchẩn wavelet dụng SWT tín hiệu đo dạng dao động riêng Ngồi nghiên cứu xét đến ảnh hưởng nhiễu đo đạc đến hệ số chi tiết phân tích wavelet Xác định dạng dao động riêng dầm có nhiều vết nứt DSM Xác định dạng dao động riêng dầm có nhiều vết nứt DSM 2.1 Các hệ thức 2.1 Các hệ thức z Et Gt rt µt Trục trung hòa h x Eb Gb rb µb b Hình 1: Dầm FGM có nhiều vết nứt Hình Dầm FGM có nhiều vết nứt Xét dầm có chiều dài L, tiết diện chữ nhật với kích thước A = b ´ h chế tạo từ vật liệu FGM (Hình 1) với hàm đặc trưng vật liệu có dạng lũy thừa 21 [21] n z h h {R( z )} = {Rb } + {Rt - Rb } ổỗ + ư÷ ; - £ z £ 2 èh 2ø (1) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng Xét dầm có chiều dài L, tiết diện chữ nhật với kích thước A = b × h chế tạo từ vật liệu FGM (Hình 1) với hàm đặc trưng vật liệu có dạng lũy thừa [21] {R(z)} = {Rb } + {Rt − Rb } z + h n h h ;− ≤ z ≤ 2 (1) R mô-đun đàn hồi E, mô-đun trượt G mật độ khối lượng ρ; số t b ký hiệu vật liệu lớp lớp dưới; z tọa độ tính từ mặt dầm Giả thiết biến dạng bé, chuyển vị dầm Timoshenko điểm tiết diện có dạng u(x, z, t) = u0 (x, t) − (z − h0 ) θ(x, t); w(x, z, t) = w0 (x, t) (2) dó u0 (x, t), w0 (x, t) chuyển vị điểm trục trung hòa; h0 khoảng cách từ trục trung hòa đến trục x; θ góc xoay tiết diện quanh trục y Cách xác định h0 xét ảnh hưởng vị trí thực đường trục trung hòa đến dao động dầm FGM trình bày [21] Ký hiệu (A11 , A12 , A22 ) = E(z) 1, z − h0 , (z − h0 )2 dA; A33 = η A (I11 , I12 , I22 ) = G(z)dA A (3) ρ(z) 1, z − h0 , (z − h0 ) dA A ma trận, véc-tơ A11 −A12 [A] = −A12 A22 0 A33 0 0 A [Π] ; = 33 −A33 −ω2 I12 ω I11 2 [D] ; = −ω I12 ω I22 − A33 0 ω2 I11 ∞ {z} = {U, Θ, W} = {u0 (x, t), θ(x, t), w0 (x, t)}T e−iωt dt T −∞ (4) Trong miền tần số, phương trình dao động dầm có dạng [21] [A] z + [Π] z + [D] {z} = {0} (5) z đạo hàm theo không gian hàm z Nghiệm tổng qt phương trình (5) khơng có vế phải viết dạng (6) {z0 (x, ω)} = [G(x, ω)]{C} với {C} = (C1 , , C6 )T số độc lập α1 ek1 x α2 ek2 x α3 ek3 x α1 e−k1 x α2 e−k2 x α3 e−k3 x k1 x k2 x k3 x −k1 x −k2 x [G(x, ω)] = e e e e e e−k3 x β1 ek1 x β2 ek2 x β3 ek3 x −β1 e−k1 x −β2 e−k2 x −β3 e−k3 x (7) Đối với dầm có vết nứt tiết diện có tọa độ e, vết nứt mơ hình hóa hai lò xo: lò xo dọc có độ cứng T lò xo xoắn có độ cứng R (Hình 2) Điều kiện liên tục vị trí vết nứt [21] U(e + 0) = U(e − 0) + γ1 U x (e); Θ(e + 0) = Θ(e − 0) + γ2 Θ x (e); W(e + 0) = W(e − 0) U x (e + 0) = U x (e − 0); Θ x (e + 0) = Θ x (e − 0); W x (e + 0) = W x (e − 0) + γ2 Θ x (e) 22 (8) {z ( x, w)} = [G( x, w )]{C} (6) với {C} = (C1 , ,C ) số độc lập T éa 1e k1x a e k2 x a e k3 x a 1e - k1x a e - k2 x a e - k3 x ù ê k1x ú [G(x, w )] = ê e e k2 x e k3 x e -k1x e -k2 x e - k3 x ú Liên,êT V.k1và học Công k3 x x cs / Tạp b e k2 x chí b eKhoa b1e -k1x nghệ - b 2Xây e -k2 xdựng - b e -k3 x úû ë b1e (7) R a h a) b) T Hình 2: Dầm FGM với vết nứt mở mơ hình hai lò xo tương đương Dầm FGM mởvếtvànứt mơđược hìnhmơ haihình lò xo tương Đối với dầm Hình có vết2.nứt tiết diệnvới cóvết tọanứt độ e, hóa bằngđương hai lò xo: lò xo dọc có độ cứng T lò xo xoắn có độ cứng R (hình 2) Điều kiện liên tục vị trí vết nứt [21] U (e + 0) = U (e - 0) + g U ¢ (e) ; Q(e + 0) = Q(e - 0) + g Q¢ (e) ; W (e + 0) = W (e - 0) x x Các tham số γ1 , γ2 (8) là1 hàm tham số vật liệu mô-đun đàn hồi, chiều(8) cao dầm, U x¢ (e + 0) = U x¢ (e - 0) ; Q¢x (e + 0) = Q¢x (e - 0) ;Wx¢ (e + 0) = Wx¢ (e - 0) + g Q¢x (e) hệ số nở ngang Poisson, bao gồm trường hợp dầm đồng Et = Eb = E0 hay RE = Các tham số g1, g2 (8) hàm tham số vật liệu mô-đun đàn hồi, chiều cao dầm, hệ số nở ngang dầm γđồng E =E =E hay R 2=1 γ1 = Poisson, A11 /T = bao 2π gồm − v2trường hσ1 fhợp (s); = A22 /Rt =b6π − vE hσ2 f2 (s) g = A11 / T = 2p (1 -n )hs f1 (s) ; g = A22 / R = 6p (1 -n )hs f (s) (9) (9) đótrong σ1 (RE , n) = ( R + n) 3RE ++nn R2R +n R + n E2 ö+ n Et 2424 æ3R s21 ((R RE ,En)+=n) ;E σ (R ; s, n) - − E E a++n αE + R a ÷ ; RE α=2 ; R = Et ( R= E , n) = ỗ E E R + 1 + n R + 3(3 + n ) + n + n Eb E ( )( ) è ø n E E (RE + 1) (1 + n) RE + 3(3 + n) 2+n 1+ Eb f1 ( s) = s (0.6272 - 0.17248s + 5.92134 s - 10.7054 2 s + 31.5685s4 - 67.47 s +5 f1 (s) = s (0,6272 − 0,17248s + 5,92134s − 10,7054s + 31,5685s − 67,47s + + 139.123s - 146.682s + 92.3552 s8 ) + 139,123s6 − 146,682s + 92,3552s8 ) (10) f ( s ) = s (0.6272 - 1.04533s + 4.5948s - 9.9736s + 20.2948s - 33.0351s + f2 (s) = s2 (0,6272 − 1,04533s + 4,5948s2+−47.1063 9,9736s + 20,2948s4 − 33,0351s5 + s - 40.7556s + 19.6s8 ) + 47,1063s6 − 40,7556s7 + 19,6s8 ) Với a = + h0 h , s = a h a chiều sâu vết nứt Ký hiệu {z c ( x, w )} nghiệm riêng (5) thỏa mãn điều kiện ban đầu (10) với α = 1/2 + h0 /h, s = a/h a chiều sâu vết Tnứt Ký hiệu {zc (x, ω)}T nghiệm riêng (5) thỏa (11) {z c (0)} = (g 1U x¢ (e),g Q¢x (e), 0) ; {z¢c (0)} = (0, 0, g Q¢x (e)) mãn điều kiện ban đầu {zc (0)} = γ1 U x (e), γ2 Θ x (e), T ; zc (0) = 0, 0, γ2 Θ x (e) T (11) {z0 (e)} = U0 (e) Θ0 (e) W0 (e) T (12) Ta nhận [21] {zc (x)} = [Φ(x)] [Σ] z0 (e) = [Gc (x)] z0 (e) (13) [Gc (x)] ma trận 3× α1 cosh k1 x α2 cosh k21 x α3 cosh k3 x cosh k21 x cosh k3 x [Gc (x)] = cosh k1 x β1 sinh k1 x β2 sinh k2 x β3 sinh k3 x δ11 δ12 δ13 δ21 δ22 δ23 δ31 δ32 δ33 γ1 0 γ2 γ2 δ11 = (k3 β3 − k2 β2 )/∆; δ12 = (α3 k2 β2 − α2 k3 β3 )/∆; δ13 = (α2 − α3 )/∆ (14) δ21 = (k1 β1 − k3 β3 )/∆; δ22 = (α1 k3 β3 − α3 k1 β1 )/∆; δ23 = (α3 − α1 )/∆ δ31 = (k2 β2 − k1 β1 )/∆; δ32 = (α2 k1 β1 − α1 k2 β2 )/∆; δ33 = (α1 − α2 )/∆ ∆ = k1 β1 (α2 − α3 ) + k2 β2 (α3 − α1 ) + k3 β3 (α1 − α2 ) Ta đưa vào ma trận hàm vết nứt G(x) = [Gc (x)] : x > [0] :x≤0 23 (15) [Gc(x)] ma trận 3×3 éa cosh k1 x a cosh k 21 x a cosh k x ù éd 11 d 12 [G c ( x)] = êê cosh k1 x cosh k 21 x cosh k x úú × êêd 21 d 22 êë b1 sinh k1 x b sinh k x b sinh k x úû êëd 31 d 32 Đối với dầm có nhiều vết d Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây (k b - k b ) D ; d 12 = (a k b - a k b ) D ; d 13 = (a d 11 =dựng d 21 = (k1 b1 - k b ) D ; d 22 = (a k b - a k1 b1 ) D ; d 23 = (a nứt, ta nhận nghiệm tổng quát phương trình (5) dạng d 31 = (k b - k1 b1 ) D ; d 32 = (a k1 b1 - a k b ) D ; d 33 = (a n D = k1 b1 (a - a ) + k b (a - a ) + k b (a - a {zc (x, ω)} = {z0 (x, ω)} + j=1 G x − ej µj Ta đưa vào ma trận hàm vết nứt (16) [G ( x)] = ìí[G [(0x]) ] ::xx£>00 c với µ j véc-tơ × ỵ Đối với dầm có nhiều vết nứt, ta nhận nghiệm tổng quát phươn j−1 µ j = z0 e j + G e j − ek k=1 n µk ; j = 1, 2, 3, , n (17) {zc ( x, w)} = {z ( x, w)} + å éëG( x - e j ) ùû.{μ j} j =1 Thay biểu thức nghiệm tổng quát (6) vào (16), ta nhận với {μ j } véc-tơ 3×1 j -1 n { μ j } = {z¢0 (e j )}+ å [G ¢(e j - ek )].{μ k } ; G x − e j χ˜ j {C} = [Ψ(x, ω)] {C} (18) {zc (x, ω)} = [G(x, ω)] + k =1 Thay biểu thức nghiệm tổng quát (6) vào (16), ta nhận j=1 G x − ej ỉ n ö è j =1 ø {z c ( x, w )} = ỗ [G( x, w)] + ộởG( x - e j ) ùû éëχ! j ùû ÷ {C} = éëΨ ( x, j−1 n [Ψ(x, ω)] = [G(x, ω)] + j =1 χ˜ j = G e j + G e j − ek χ˜ k ; j = 1, 2, 3, , n j -1 n k=1 ëéΨ ( x, w )ûù = [G( x, w )] + å ëéG( x - e j ) ûù ëéχ! j ûù ; ëéχ! j ûù = ëéG¢(e j ) ûù + å ëéG¢(e j - e j=1 (19) k =1 j =1 Ký hiệu tọa độ nút, lực nút phần tử chịu uốn kéo, nén đồng thời Hình Ký hiệu tọa độ nút, lực nút phần tử chịu uốn kéo, nén đ ˆ e = {U1 , Θ1 , W1 , U2 , Θ2 , W2 }T ; U T M1 ,2W, 1Q {Pe } = {N1 , M1{,UˆQe }1=, N {U21,, Q ,U22}, Q ,W2 }T ; {Pe }(20) = {N1 , M , Q1 , N , M Ta nhận [21] ˆ e (ω) · U ˆe {Pe (ω)} = K [−BF (Y) x=0 ] [BF (Ψ) x=L ] với [BF ] toán tử ma trận [Ψ(0, ω)] [Ψ(L, ω)] L N1 i M1 [ ] x N2 ˆ m K e phần tử dầm j W1 W2 −1 U1 Q1 (22) {Pe (w Q2 Q1 (21) ˆ e ma trận độ cứng động lực phần K tử dầm FGM có nhiều vết nứt ˆe = K ta nhận [21 y Q2 M2 ˆ ] = é[- B F [K e ê [B ë F U2 với [BF] tốn Hình 3: Thanh chịu uốn kéo, nén đồng thời Hình Thanh chịu uốn kéo, nén đồng thời A11 ∂ x −A12 ∂ x [BF ] = −A12 ∂ x A22 ∂ x −A33 A33 ∂ x (23) Sau đó, việc lắp ghép ma trận độ cứng động lực véc-tơ tải trọng quy nút phần tử vào ma trận độ cứng động lực véc-tơ tải trọng quy nút kết cấu thực khơng khác phương pháp phần tử hữu hạn Bài toán dao động riêng giải phương trình sau ˆ ˆ = {0} [K(ω)]{ U} (24) Với tần số dao động riêng {ω} = {ω1 ω2 ωn } có từ phương trình ˆ det [K(ω)] =0 (25) Mỗi nghiệm ω j phương trình tần số dao động riêng kết cấu ứng với dạng dao động riêng φ j (x) = C 0j Ψ x, ω j [Ψ(0, ω)] [Ψ(L, ω)] −1 ˆj U ˆ j nghiệm chuẩn hóa (24) ứng với ω j C 0j số bất kỳ, U 24 (26) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng Phân tích wavelet SWT Phân tích wavelet bắt đầu việc chọn hàm wavelet ψ(x) (gọi wavelet mẹ) Phân tích wavelet liên tục (CWT) định nghĩa [21] ∞ C(a, b) = √ a ∞ x−b f (x)ψ dx = a −∞ f (x)ψa,b (x)dx (27) −∞ a > tỷ lệ b tham số dịch mức; ψa,b (x) hàm số x−b ψa,b (x) = √ ψ a a (28) Kết CWT hệ số wavelet C(a, b) thể tương quan hàm số wavelet tín hiệu phân tích f (x) Vì thế, thay đổi đột ngột f (x) tạo hệ số wavelet có biên độ lớn, đặc điểm để xây dựng phương pháp nhận dạng vết nứt dựa phân tích wavelet tín hiệu Tín hiệu ban đầu f (x) tái tạo từ hệ số wavelet C(a, b) ∞ ∞ f (x) = Kψ C(a, b)ψa,b (x) dbda a2 (29) −∞ −∞ số Kψ phụ thuộc vào loại wavelet Giả thiết hệ số wavelet C(a, b) có giá trị tỷ lệ a < a0 , dành riêng cho thành phần tần số cao tín hiệu, với a > a0 , xem nhiễu Trong trường hợp này, tín hiệu tái tạo cần phần bù tín hiệu tương ứng với a > a0 Để thực việc này, người ta đưa vào hàm φ(x) khác gọi “hàm tỷ lệ”, thu hệ số wavelet ∞ D (a0 , b) = √ a0 −∞ ∞ x−b f (x)φ dx = a0 f (x)φa0 ,b (x)dx (30) −∞ Hàm tỷ lệ cần thiết cho tính tốn số Thay (29), tín hiệu ban đầu f (x) tái tạo từ a0 ∞ ∞ dbda C(a, b)ψa,b (x) + Kψ a0 a f (x) = Kψ a=0 b=−∞ D (a0 , b) φa0 ,b (x)db (31) b=−∞ Một nhược điểm CWT trình phân tích cho số lượng lớn hệ số wavelet C(a, b) Để giảm khối lượng tính tốn, phân tích wavelet rời rạc (DWT) sử dụng tỷ lệ rời rạc tham số dịch mức dạng cặp số: a = j ; b = k2 j j k số nguyên, j mức dyadic Phân tích wavelet rời rạc DWT sau ∞ C j,k = ∞ f (x)ψ x − k dx = − j/2 −j −∞ f (x)ψ j,k (x)dx (32) −∞ ψi,k (x) hàm wavelet rời rạc ψ j,k (x) = 2− j/2 ψ 2− j x − k 25 (33) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng Thay cho (29), tín hiệu DWT tái tạo từ hệ số wavelet C j,k ∞ ∞ f (x) = C j,k 2− j/2 ψ 2− j x − k (34) j=−∞ k=−∞ Tín hiệu qua nhiều lọc (gồm có lọc cao lọc thấp) để tách lấy thành phần tần số cao tần số thấp tương ứng Thay cho (31), tín hiệu DWT biểu diễn hàm xấp xỉ chi tiết sau J ∞ ∞ f (x) = cD j (k)ψ j,k (x) + cA J (k)φ j,k (x) = D j (x) + A j (x) (35) j=−∞ k=∞ j≤J k=−∞ với A j (x) hàm xấp xỉ mức J; Di (x) chi tiết mức j ≤ J ∞ D j (x) = ∞ cD j (k)ψ j,k (x); A j (x) = k=∞ cA J (k)φ j,k (x) (36) f (x)φ J,k (x)dx (37) k=∞ cD j cA j theo thứ tự hệ số chi tiết hệ số xấp xỉ ∞ cD J (k) = ∞ f (x)ψ J,k (x)dx; cA J (k) = −∞ −∞ Dưới đây, tác giả quan tâm nhiều tới tín hiệu chi tiết Nếu f (x) tín hiệu ứng xử kết cấu đường độ võng, tín hiệu D j (x) chứa thơng tin cần thiết để phát vết nứt Tuy nhiên DWT cổ điển có nhược điểm khơng phải biến đổi bất biến theo thời gian Điều có nghĩa là, với tín hiệu tuần hồn, DWT phiên dịch tín hiệu gốc f (x) dịch DWT tín hiệu gốc f (x) Để khắc phục vấn đề này, ta sử dụng đến phần dư phân tích tín hiệu [19] ∞ ∞ D j,k = − j/2 x−k dx; f (x)ψ 2j A j,k = − j/2 f (x)φ x−k dx 2j (38) −∞ −∞ Hệ số xấp xỉ chi tiết xác định theo (37) gọi phân tích wavelet dừng (SWT) Cần phải ý SWT liệu gốc khơng bị tiêu hao nhiều Có nghĩa là, kích thước dãy số liệu biến đối sau SWT không bị cắt phần Ngược lại, DWT, kích thước dãy liệu sau biến đổi nửa so với kích thước tín hiệu gốc Do đó, DWT q trình phân giải tín hiệu mà kết cho nghèo nàn tín hiệu gốc Trong đó, SWT q trình tách tín hiệu kết nhiều tín hiệu gốc Do đó, hệ số chi tiết từ phân tích DWT có thơng tin so với phân tích SWT Bởi vậy, SWT có tiềm lớn với thuận lợi tách nhận dạng điểm bật tín hiệu xác định hư hỏng kết cấu 3.1 Nhiễu đo đạc khử nhiễu Thực tế, liệu dạng dao động riêng kết cấu có vết nứt gồm phần [19] y = yintact + ynoise + ycrack 26 (39) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng y dạng dao động riêng đo lường được; yintact dạng dao động riêng kết cấu khơng có vết nứt khơng có nhiễu; ynoise thành phần nhiễu có mặt trường hợp kết cấu có nứt khơng có nứt mà ta khơng biết giá trị hay không đo lường được; ycrack thành phần tín hiệu xuất thêm vào có vết nứt tồn kết cấu có vết nứt Về mặt lý thuyết, liệu dạng dao động phân tích SWT thành hai thành phần: Thành phần thứ hệ số xấp xỉ có chứa yintact , đường cong trơn; Thành phần thứ hai hệ số chi tiết gồm ynoise ycrack lưu ý thành phần nhiễu bao gồm hai thành phần kết cấu khơng nứt kết cấu có nứt phân tách vào hệ số chi tiết Hệ số chi tiết cung cấp thơng tin nhận biết tồn vết nứt Để thuận lợi cho nhận dạng vết nứt việc lựa chọn phương pháp giảm thiểu nhiễu hay khử nhiễu cho tín hiệu kết cấu khơng có nứt có nứt cần thiết Thơng thường, quy trình khử nhiễu phân tích wavelet thực thơng qua đặt ngưỡng (thresholding) định sẵn Có hai ngưỡng hay sử dụng hard-thresholding (ngưỡng cứng) soft-thresholding (ngưỡng mềm) Hàm hard-thresholding định nghĩa ηth = x |x| > th |x| < th th ngưỡng (40) Hàm ngưỡng mềm soft-thresholding định nghĩa ηth = sgn(x) max(|x| − th, 0) (41) Hàm soft-thresholding sử dụng nhiều việc khử nhiễu Đối với hàm ngưỡng mức ngưỡng th tính tốn theo cơng thức th = σ log N (42) N chiều dài tín hiệu, σ độ lệch chuẩn nhiễu Trong hai trường hợp kết cấu có vết nứt kết cấu khơng có vết nứt, nhiễu có tần số số cao tần số thấp Đối với kết cấu có vết nứt, nhiễu tác động thêm vào tín hiệu hệ số chi tiết phát sinh có vết nứt Như nhận dạng vết nứt kết cấu, tác động mức nhiễu vào tần số thấp bỏ qua coi có hệ số xấp xỉ kết phân tích wavelet, không ảnh hưởng đến việc nhận dạng vết nứt theo hệ số chi tiết Để mô số liệu đo thực nghiệm, ta sử dụng nhiễu trắng Gaussian cộng thêm vào dạng dao động kết cấu [21] − (x−µ) f (x) = √ e 2σ2 (43) σ 2π Thông thường đường cong Gaussian phụ thuộc vào giá trị kỳ vọng µ phương sai σ2 , đánh giá theo trị số SNR (Signal to Noise Ratio) SNR = 20 log10 norm(Signal) norm(Noise) (dB) (44) norm chuẩn đo tín hiệu f (x) có độ dài N s N 1/2 s norm( f ) = | f (xi )|2 i=1 27 (45) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng Kết số, nhận xét đánh giá với N s số điểm rời rạc mẫu tín hiệu f (x) Chỉ số SNR cao tương ứng với diện nhiễu 3.1 Dầm đơn giản nhỏ ngược lại Như vậy, véc-tơ tín hiệu nhiễu có dạng: Xét dầm FGM Timoshenko hai đầu gối tựa cố định với tham số: Et=70GPa; rt=2780kg/m3; ynr vànorm intact ) L=1,0m; b=0,1m; h=0,1m µt=0,33; Et/Eb=2; rb=7850kg/m3y; noise µb=0,33; n=0,5 kích(ythước: = (46) (0.05×SNR) norm (ynr ) 10 Hình biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết nứt ytại x1=0,2m tính từ nút bên trái với độ sâu vết nứt 10%, 20% 30% ứng với số điểm đo với nr = Rand size yintact,std véc-tơ giá trị giả lập ngẫu nhiên phân bố đoạn (0, 1), 50 điểm (hình 4a-c), 200 điểm (hình 4g-i) Trong tính tốn đây, ta có chiều dài tương ứng100 với điểm chiều (hình dài tín4d-f), hiệu gốc chọn số điểm đo 100 điểm Tín hiệu đầu vào dạng dao động riêng dầm FGM có vết nứt cộng thêm nhiễu với mức nhiễu SNR khác Hình thể hệ số chi tiết nhận xétcho đánh giá động riêng thứ (5a-c), thứ hai (5d-f) thứ ba (5g-i) đốiKết với phânsố,tích SWT dạng dao dầm đơn giản FGM với mức nhiễu 75, 80 90dB Ta nhận thấy dạng dao động riêng cao 4.1 Dầm đơn giản ảnh hưởng nhiễu giảm dần Để nhận dạng vị trí vết nứt, mức nhiễu dạngXét daodầm động riêng phải dạng hai,ρ giá trị FGM Timoshenko hailớn đầuhơn gối80dB, tựa cốtrong định với cácvới tham số:dao Et động = 70thứ GPa; t = 2780 là3 75dB Dưới đây, để đơn giản ta giả thiết3 khơng có nhiễu dạng dao động riêng Hình kg/m ; µt = 0,33; Et /Eb = 2; ρb = 7850 kg/m ; µb = 0,33; n = 0,5 kích thước: L = 1,0 m; b = biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết nứt 0,1 m; h =tính 0,1từ m.nút bên trái với độ sâu 30% số mũ n thay đổi n=0,1; 1; 10 x1=0,2m a) b) c) d) e) f) g) h) i) Hình 4: Hệ số chi tiết SWT ba dạng dao động riêng dầm FGM có vết nứt vị Hình Hệ số chi tiết SWT ba dạng dao động riêng dầm FGM có vết nứt vị trí 0,2 m; trí 0,2m; độ sâu vết nứt 10%; 20%; 30% với số điểm đo 50(a-c); 100(d-e) 200 điểm(g-i) độ sâu vết nứt 10%; 20%; 30% với số điểm đo 50(a-c); 100(d-e) 200 điểm(g-i) Hình biểu đồ hệ số chi tiết SWT loại db4 28 cho dạng dao động dầm có vết nứt x1=0,2m tính từ nút bên trái với độ sâu 30% tỷ số Et/Eb thay đổi 0,5; 1; Hình biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết nứt cách khoảngLiên, 0,2m với nứt thay nghệ đổi Xây 10%, T V độ cs /sâu Tạp vết chí Khoa học Cơng dựng 20% 30% Số điểm đo hình 5-7 100 điểm Hình biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết Ta nhận thấy: nứt x1 = 0,2m tính từ nút bên trái với độ sâu vết nứt 10%, 20% 30% ứng với số điểm đo 50 -điểm Biểu đồ hệ số chi tiết các(Hình dạng dao động khác nhau(Hình có4(g)–(i)) điểm gián đoạncác tính vị trítốn vết nứt (Hình 4(a)–(c)), 100của điểm 4(d)–(f)), 200 điểm Trong ta chọncủa số điểm đo làđồ100 điểm dạng -đây,Giá trị đỉnh biểu tăng lên Tín hiệu chiềuđầu sâuvào vếtlànứt tăng dao lên, động nghĩariêng vếtđầu nứttiên lớncủa dầm biểu FGM nứt tiết cộng đồcó hệvết số chi rõthêm nét nhiễu với mức nhiễu SNR khác Hình thể hệ số chi tiết phân tích SWT cho dạng dao động riêng thứ (Hình 5(a)–(c)), thứ hai (Hình 5(d)–(f)) - Đối với vết nứt có độ sâu vị trí khác biên độ đỉnh thứ (Hìnhbiên 5(g)–(i)) củalớn dầm đơnvới giản FGM 75, 90 Ta nhận khácbanhau, độ đỉnh ứng vị trí vết với nứtmức làm nhiễu thay đổi lớn80 dạng daodB động Như thấy vậy, dạng dao động riêng cao ảnh hưởng nhiễu giảm dần Để nhận dạng biên độ đỉnh phụ thuộc vào vị trí vết nứt mà phụ thuộc độ sâu vết nứt vị trí vết nứt, mức nhiễu dạng dao động riêng phải lớn 80 dB, với dạng - Khi tăng số lượng điểm lấy mẫu ảnh hưởng chiều sâu vết nứt đến dạng gián đoạn dao động thứ hai, giá trị 75 dB Dưới đây, để đơn giản ta giả thiết khơng có nhiễu biểu đồ hệ số chi tiết ổn định, đồng thời giá trị tuyệt đối biên độ tăng lên đáng kể dạng dao xảy độngrariêng Hình biểuthu đồhẹp hệ số vùng gián đoạn lại.chi tiết SWT loại db4 cho dạng dao động dầm có vết nứt x1 = 0,2 m tính từ nút bên trái với độ sâu 30% số mũ n thay đổi -là nKhi số mũ n nhỏ tỷ số Et/Eb lớn dầm nhạy cảm với vết nứt, biểu đồ hệ = 0,1; 1; 10 số chi tiết có bước nhảy lớn a) b) c) d) e) f) g) h) i) Hình 5: Hệ số chi tiết SWT dạng dao động riêng thứ (a-c), thứ hai (d-f) thứ ba (g-i) Hình 5.dHệ chi tiết đốiởvới dạng daotừđộng riêng (a-c),vàthứ hainhiễu (d-f) (g-i) dầm ầmsố FGM có 1SWT vết nứt vị trí 0,2m nút trái vớithứ độnhất sâu 30% mức 75,thứ 80 ba 90dB FGM có vết nứt vị trí 0, m từ nút trái với độ sâu 30% mức nhiễu 75, 80 90 dB Hình biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết 29 a) b) c) d) f) e) Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng nứt x1 = 0,2 m tính từ nút bên trái với độ sâu 30% tỷ số Et /Eb thay đổi 0,5; 1; Hình biểu đồ hệ số chi tiết SWT loại db4 cho dạng dao động dầm có vết nứt cách khoảng 0,2 m với độ sâu vết nứt thay đổi 10%, 20% 30% Số điểm đo Hình 5–7 100 điểm g) h) i) a) b) c) Hình 6.cóHệ với độ ba dạng daonứt động riêng dầm có vết nứt vị số vếtchi nứttiếtở SWT vị trí đối 0,2m, sâu vết 30% vàđầu số tiên mũ thay đổiđơn lần giản lượt FGM n=0,1; 1; 10 b) n = 0,1; 1; 10 c) trí 0,2 m, độ sâu a) vết nứt 30% số mũ thay đổi Hình 6: Hệ số chi tiết SWT ba dạng dao động riêng dầm đơn giản FGM có vết nứt vị trí 0,2m, độ a) sâu vết nứt 30% số mũ thay b) đổi n=0,1; 1; 10 c) Hình 6: Hệ số chi tiết SWT ba dạng dao động riêng dầm đơn giản FGM có vết nứt vị trí 0,2m, độ sâu vết nứt 30% số mũ thay đổi n=0,1; 1; 10 a) b) c) Hình 7.cóHệ1 số tiếtở SWT với ba dao nứt động30% riêngvàđầu dầm đơn giản FGM có vết nứt vị vếtchi nứt vị tríđối 0,2m, độdạng sâu vết tỷ tiên số Ecủa t/Eb thay đổi 0,5; 1; c) trí 0,2 m, độ sâua)vết nứt 30% tỷ số Et /Eb thay đổib)lần lượt 0,5; 1; Hình 7: Hệ số chi tiết SWT ba dạng dao động riêng dầm đơn giản FGM có vết nứt vị trí 0,2m, độ sâu vết nứt 30% tỷ số Et/Eb thay đổi 0,5; 1; a) b) c) Hình 8: Hệ số chi tiết SWT ba dạng dao động riêng dầm đơn giản FGM có vết nứt phân bố đềua)trên dầm độ sâu vết nứt đềub) thay đổi 10%, 20%, 30% c) M 3.2 liênchi tục nhiều nhịpvới ba dạng dao động riêng dầm đơn giản FGM có vết nứt phân HìnhDầm Hệ có số vết tiết nứtSWT phânđối bố dầm độ sâu vết nứt thay đổi 10%, 20%, 30% bố dầm độ sâu vết nứt thay đổi 10%, 20%, 30% h 3.2 Dầm liên tục nhiều nhịp L2=1,2m L3=0,6m b L1=0,7m h Hình 9: Dầm liên tục có vết nứt 30 L2=1,2m L3=0,6m b L1=0,7m Xét dầm liên tục FGM với tham số vật liệu sau: Et=70GPa; rt=2780kg/m3; µt=0,33; Hình Dầm liênvật tụcliệu có vết nứt Kích thước tiết diện: b=0,1m, Et/Eb=2; rb=7850kg/m3; µb=0.33 số mũ9:đặc trưng n=0,5 h=0,1m Xét dầm(hình liên 9) tục FGM với tham số vật liệu sau: Et=70GPa; rt=2780kg/m3; µt=0,33; Liên, T V cs / Tạp chí Khoa học Công nghệ Xây dựng Ta nhận thấy: - Biểu đồ hệ số chi tiết dạng dao động khác có điểm gián đoạn vị trí vết nứt - Giá trị đỉnh biểu đồ tăng lên chiều sâu vết nứt tăng lên, nghĩa vết nứt lớn biểu đồ hệ số chi tiết rõ nét - Đối với vết nứt có độ sâu vị trí khác biên độ đỉnh khác nhau, biên độ đỉnh lớn ứng với vị trí vết nứt làm thay đổi lớn dạng dao động Như vậy, biên độ đỉnh phụ thuộc khơng vào vị trí vết nứt mà phụ thuộc độ sâu vết nứt - Khi tăng số lượng điểm lấy mẫu ảnh hưởng chiều sâu vết nứt đến dạng gián đoạn biểu đồ hệ số chi tiết ổn định, đồng thời giá trị tuyệt đối biên độ tăng lên đáng kể vùng xảy gián đoạn thu hẹp lại - Khi số mũ n nhỏ tỷ số Et /Eb lớn dầm nhạy cảm với vết nứt, biểu đồ hệ số chi tiết có bước nhảy lớn 4.2 Dầm liên tục nhiều nhịp Xét dầm liên tục FGM với tham số vật liệu sau: Et = 70 GPa; ρt = 2780 kg/m3 ; µt = 0,33; Et /Eb = 2; ρb = 7850 kg/m3 ; µb = 0,33 số mũ đặc trưng vật liệu n = 0,5 Kích thước tiết diện: b = 0,1m, h =liên 0,1m (Hình 9).nhịp 3.2 Dầm tụcnhiều nhiềunhịp 3.2 Dầm liên tục L1L =0,7m 1=0,7m L2=1,2m L2=1,2m h h L3=0,6mb b L3=0,6m Hình 9: Dầm liênliên tụccó cóvết vết nứt nứt Hình Dầm tục cónứt vết Hình 9: Dầm liên tục Xét số số vậtvật liệuliệu nhưnhư sau:sau: Et=70GPa; rt=2780kg/m ; µt3=0,33; Xét dầm dầmliên liêntục tụcFGM FGMvới vớicáccáctham tham Et=70GPa; rt=2780kg/m ; µt=0,33; Hình 10 biểu đồ 3 hệ số chi tiết wavelet ba dạng dao động dầm có vết nứt, vị trí EEt/E =2; r =7850kg/m ; µ =0.33 số mũ đặc trưng vật liệu n=0,5 Kích thước tiết diện: b=0,1m, b b b t/Eb=2; rb=7850kg/m ; µb=0.33 số mũ đặc trưng vật liệu n=0,5 Kích thước tiết diện: b=0,1m, vết nứt(hình x1 9) =9) 0,2 m, tính từ nút bên trái (trên nhịp dầm đầu tiên), chiều sâu vết nứt thay đổi h=0,1m h=0,1m (hình 10%, 20%, 30% Hình 11 biểu đồ hệ số chi tiết wavelet ba dạng dao động Hình 10 làlàbiểu đồđồhệhệsốsốchichi tiếttiết wavelet ba ba dạng daodao động đầuđầu tiêntiên của dầmdầm có 1có vết1 nứt, vị trívịvết Hình 10 biểu wavelet dạng động vết nứt, trí vết dầm có vết nứt, vị trí vết nứt x = 0,2 m, tính từ nút thứ hai bên trái (trên nhịp dầm thứ chiều nứt x =0,2m, tính từ nút bên trái (trên nhịp dầm đầu tiên), chiều sâu vết nứt thay đổi2), lần nứt x1 1=0,2m, tính từ nút bên trái (trên nhịp dầm đầu tiên), chiều sâu vết nứt thay đổi lần sâu vết nứt thay đổi lần Hình lượt 10%, 20%, 30% 12 biểu đồ số chi wavelet ba dạng lượt là biểu đồđồ hệ hệ sốHình chichi tiết wavelet bahệ dạng daotiết động đầu dao lượtlàlà10%, 10%,20%, 20%,30% 30% Hình1111 biểu số tiết wavelet ba dạng dao động tiên động dầm có vết nứt, vị trí vết nứt x = 0,2 m, tính từ nút thứ ba (trên nhịp dầm thứ dầm có vết nứt, vị trí vết nứt x =0,2m, tính từ nút thứ hai bên trái (trên nhịp dầm thứ 2), chiều 1 dầm có vết nứt, vị trí vết nứt x1=0,2m, tính từ nút thứ hai bên trái (trên nhịp dầm thứ 2), chiều sâu vết nứt thay đổi 10%, 20%, 30% Hình 12 biểu đồ hệ số chi tiết wavelet ba dạng 3), chiều sâu vết nứt thay đổi 10%, 20%, 30% s vết nứt thay đổi 10%, 20%, 30% Hình 12 biểu đồ hệ số chi tiết wavelet ba dạng dao động dầm có vết nứt, vị trí vết nứt x1=0,2m, tính từ nút thứ ba (trên nhịp dầm thứ 3), chiều sâu vết nứt thay đổi 10%, 20%, 30% â u 10 b) c) a) c) a) với ba dạng dao động riêng đầub)tiên dầm liên tục FGM Hình 10: Hệ số chi tiết SWT đối có vết nứt vị trí 0,2m nhịp dầm đầu tiên, độ sâu vết nứt thay đổi a/h=10%, 20%, 30% có vết nứt vị trí 0,2m nhịp dầm đầu tiên, độ sâu vết nứt thay đổi a/h=10%, 20%, Hình 10 Hệ số chi tiết SWT ba dạng dao động riêng dầm liên tục FGM có vết nứt vị trí 0,2 m nhịp dầm đầu tiên, độ sâu vết nứt thay đổi a/h =10%, 20%, 30% Ta nhận thấy: - Tương tự với trường hợp dầm đơn giản, biểu đồ hệ số chi tiết có điểm gián đoạn vị trí vết nứt 31 b) a) c) Hình 10: Hệ số chi tiết SWTa)đối với ba dạng dao động riêng b) dầm liên tục FGM c) Liên, T V cs / Tạp chítiên, Khoađộ họcsâu Cơng nghệ Xây dựng cóHình vết10: nứt vị trí nhịp dầm vết nứtđầu thay đổi a/h=10%, 30% Hệở số chi0,2m tiết SWT bađầu dạng dao động riêng tiên dầm liên20%, tục FGM có vết nứt vị trí 0,2m nhịp dầm đầu tiên, độ sâu vết nứt thay đổi a/h=10%, 20%, 30% a) b) c) Hình 11: Hệ số chi tiết SWT ba dạng dao động riêng dầm liên tục FGM có a) b) c) vết nứt vị trí 0,2m nhịp dầm thứ hai, độ sâu vết nứt thay đổi a/h=10%, 20%, 30% Hình 11: Hệ số chi tiết SWT ba dạng dao động riêng dầm liên tục FGM có Hình 11 Hệ số chi tiết SWT ba dạng dao động riêng dầm liên tục FGM có vết nứt vị vết nứt ởtrí vị0,2 trí m 0,2m dầmhai, thứđộhai, a/h=10%, 20%, 30% trêntrên nhịpnhịp dầm thứ sâuđộ vếtsâu nứt vết thaynứt đổi thay a/h =đổi 10%, 20%, 30% b) c) a) Hình 12: Hệ số chi tiết SWTa)đối với ba dạng dao động riêng b) dầm liên tục FGM c) có vết trí tiết 0,2m đối nhịpvới dầm độđộng sâu vết nứtđầu thaytiên đổicủa a/h=10%, 30% Hình 12:nứt Hệởsốvịchi SWT ba thứ dạngba, dao riêng dầm liên20%, tục FGM vếtsốnứt vị trí 0,2m nhịp dầm ba,riêng độ sâu Hìnhcó 12.1 Hệ chiởtiết SWT đốitrên với ba dạng daothứ động đầuvết tiênnứt củathay dầm đổi liêna/h=10%, tục FGM có20%, vết30% nứt vị m trênnhịp nhịp dầm thứ ba, độ sâu vết nứt thay đổi a/h = 10%, 20%, 30% 3.2 Dầm liên trí tục0,2 nhiều 3.2.- Giá Dầmtrịliên tục nhiều nhịpđồ tăng lên chiều sâu vết nứt tăng lên, nghĩa vết nứt lớn biểu đỉnh biểu h đồ hệ số chi tiết rõ nét L2=1,2m L3=0,6m b h L1=0,7m - Đối với vết nứt có độ sâu vị trí nhau, nhịp dầm khác L2=1,2m L3khác =0,6m L1=0,7m b 9: Dầm tục có biên độ đỉnh khác nhau, Hình biên độ đỉnh liên lớn ứng vớivết vị nứt trí vết nứt làm thay đổi lớn dạng Xétđộng dầm Như liên vậy, tục FGM cácphụ tham số khơng vậtDầm liệu nhưtục sau: Etrít=70GPa; rt=2780kg/m ; µđộ Hình 9: liên có nứtnứt mà t=0,33; dao biên độvới đỉnh thuộc vào vị vết vết phụ thuộc sâu 3 E /E =2; r =7850kg/m ; µ =0.33 số mũ đặc trưng vật liệu n=0,5 Kích thước tiết diện: b=0,1m, t nứt bdầm liên b b tham số vật liệu sau: Et=70GPa; rt=2780kg/m ; µt=0,33; vết Xét tục FGM với h=0,1m 9) Et/Eb=2;(hình rb=7850kg/m ; µb=0.33 số mũ đặc trưng vật liệu n=0,5 Kích thước tiết diện: b=0,1m, h=0,1m (hình biểu9).đồ hệ số chi tiết wavelet ba dạng dao động dầm có vết nứt, vị trí vết 5.Hình Kết10 luận nứt 10 x1=0,2m, nútchiđầu bên trái dầm đầu chiều nứtnứt, thayvịđổi Hình biểu tính đồ hệtừsố tiếttiên wavelet ba (trên dạng nhịp dao động đầu tiên), tiên dầmsâu cóvết vết trí lần vết lượt 10%, 20%, 30% Hình 11 biểu đồ hệ số chi tiết wavelet ba dạng dao động Trong báo này, tác giả trình bày kết nghiên cứu việc xác định vết nứt nứt x1=0,2m, tính từ nút bên trái (trên nhịp dầm đầu tiên), chiều sâu vết nứt thay đổi lần dầm vết nứt, vị30% tríđơn vếtHình nứt 11 xlà nút bênvết trái (trên dầm 2), chiều 1=0,2m, dầm FGM (dầm giản, dầm liên nhiều cóhai nhiều dựanhịp trênđộng phânthứ tích SWT lượt có 10%, 20%, biểutục đồtính hệ từ số nhịp) chi thứ tiết wavelet banứt dạng dao sâu vết nứt thay đổi là 10%, 20%, 30% Hình 12 làhai biểu đồtrái hệ số chi tiếtdầm wavelet badạng dạng dầm cócác vết nứt, vị động trí vết nứt tính từ nút thứnhiễu bên (trên nhịp thứcác 2), chiều dạng dao riêng cóxkể đến ảnh hưởng trắng Gausian Tần số 1=0,2m, dao động dầm có vết nứt, vị trí vết nứt x =0,2m, tính từ nút thứ ba (trên nhịp dầm sâuđộng vết nứt thay đổixác lầnđịnh lượttừlàmơ 10%, 20%, Hìnhđàn 12hồi biểu số chiuốn tiếtcówavelet ba nứt dạng dao riêng hình phần30% tử chịu đồ kéo,hệnén nhiều vết thứ 3), chiều sâu vết nứt thay đổi 10%, 20%, 30% dao mơ động đầulòtiên dầm có 1pháp vết nứt, vị trí động vết nứt x1=0,2m, tính từ thứthấy ba (trên nhịp dầm theo hình xo phương độ cứng lực.làKết nghiên cứunút cho phương pháp thứ 3), chiều sâu vết nứt thay đổi 10%, 20%, 30% có tính khả thi, hiệu áp dụng vào thực tế Lời cảm ơn Tác giả chân thành cảm ơn hỗ trợ tài Quỹ phát triển khoa học công nghệ quốc gia (NAFOSTED) cho đề tài mã số 107.02-2017.301 10 10 32 Liên, T V cs / Tạp chí Khoa học Cơng nghệ Xây dựng Tài liệu tham khảo [1] Ke, L L., Yang, J., Kitipornchai, S., Xiang, Y (2009) Flexural vibration and elastic buckling of a cracked Timoshenko beam made of functionally graded materials Mechanics of Advanced Materials and Structures, 16(6):488–502 [2] Yang, J., Chen, Y., Xiang, Y., Jia, X L (2008) Free and forced vibration of cracked inhomogeneous beams under an axial force and a moving load Journal of Sound and Vibration, 312(1-2):166–181 [3] Aydin, K (2013) Free vibration of functionally graded beams with arbitrary number of surface cracks European Journal of Mechanics-A/Solids, 42:112–124 [4] Yang, J., Chen, Y (2008) Free vibration and buckling analyses of functionally graded beams with edge cracks Composite Structures, 83(1):48–60 [5] Wei, D., Liu, Y., Xiang, Z (2012) An analytical method for free vibration analysis of functionally graded beams with edge cracks Journal of Sound and Vibration, 331(7):1686–1700 [6] Sherafatnia, K., Farrahi, G., Faghidian, S A (2013) Analytic approach to free vibration and buckling analysis of functionally graded beams with edge cracks using four engineering beam theories International Journal of Engineering-Transactions C: Aspects, 27(6):979–990 [7] Kitipornchai, S., Ke, L L., Yang, J., Xiang, Y (2009) Nonlinear vibration of edge cracked functionally graded Timoshenko beams Journal of Sound and Vibration, 324(3-5):962–982 [8] Yu, Z., Chu, F (2009) Identification of crack in functionally graded material beams using the p-version of finite element method Journal of Sound and Vibration, 325(1-2):69–84 [9] Akbas¸, S¸ D (2013) Free vibration characteristics of edge cracked functionally graded beams by using finite element method International Journal of Engineering Trends and Technology, 4(10):4590–4597 [10] Banerjee, A., Panigrahi, B., Pohit, G (2016) Crack modelling and detection in Timoshenko FGM beam under transverse vibration using frequency contour and response surface model with GA Nondestructive Testing and Evaluation, 31(2):142–164 [11] Su, H., Banerjee, J (2015) Development of dynamic stiffness method for free vibration of functionally graded Timoshenko beams Computers & Structures, 147:107–116 [12] Lien, T V., Duc, N T., Khiem, N T (2017) Mode shape analysis of multiple cracked functionally graded beam-like structures by using dynamic stiffness method Vietnam Journal of Mechanics, 39(3):215–228 [13] Lien, T V., Duc, N T., Khiem, N T (2017) Free vibration analysis of multiple cracked functionally graded Timoshenko beams Latin American Journal of Solids and Structures, 14(9):1752–1766 [14] Khiem, N T., Lien, T V (2002) The dynamic stiffness matrix method in forced vibration analysis of multiple-cracked beam Journal of Sound and Vibration, 254(3):541–555 [15] Sohn, H., Farrar, C R., Hemez, F., Czarnecki, J (2003) A review of structural health monitoring literature: 1996–2001 Los Alamos National Laboratory, USA [16] Liew, K M., Wang, Q (1998) Application of wavelet theory for crack identification in structures Journal of Engineering Mechanics, 124(2):152–157 [17] Wang, Q., Deng, X (1999) Damage detection with spatial wavelets International Journal of Solids and Structures, 36(23):3443–3468 [18] Chang, C C., Chen, L W (2005) Detection of the location and size of cracks in the multiple cracked beam by spatial wavelet based approach Mechanical Systems and Signal Processing, 19(1):139–155 [19] Zhong, S., Oyadiji, S O (2007) Crack detection in simply supported beams without baseline modal parameters by stationary wavelet transform Mechanical Systems and Signal Processing, 21(4):1853– 1884 [20] Khiem, N T., Huyen, N N (2017) A method for crack identification in functionally graded Timoshenko beam Nondestructive Testing and Evaluation, 32(3):319–341 [21] Liên, T V., Khiêm, N T (2017) Phương pháp độ cứng động lực phân tích chẩn đoán kết cấu Nhà xuất Xây dựng 33 ... đồđồhệhệsốsốchichi tiếttiết wavelet ba ba dạng daodao động đầuđầu tiêntiên của dầmdầm có 1có vết1 nứt, vị trív vết Hình 10 biểu wavelet dạng động vết nứt, trí vết dầm có vết nứt, vị trí vết nứt x = 0,2 m,... để xác định vị trí độ sâu dầm có nhiều vết nứt dựa phân tích khơng gian wavelet dạng dao động riêng Việc dự đoán độ sâu vết nứt đơn giản hóa cách sử dụng vị trí vết nứt dự đốn tần số dao động riêng. .. 4: Hệ số chi tiết SWT ba dạng dao động riêng dầm FGM có vết nứt vị Hình Hệ số chi tiết SWT ba dạng dao động riêng dầm FGM có vết nứt vị trí 0,2 m; trí 0,2m; độ sâu vết nứt 10%; 20%; 30% với số