Đánh giá một số phương pháp chẩn đoán hư hỏng kết cấu dầm từ kết quả phân tích dao động

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐỖ NGỌC TRIỀU ĐÁNH GIÁ MỘT SỐ PHƯƠNG PHÁP CHẨN ĐOÁN HƯ HỎNG KẾT CẤU DẦM TỪ KẾT QUẢ PHÂN TÍCH DAO ĐỘNG Chuyên ngành : Xây dựng cơng trình dân dụng cơng nghiệp Mã số ngành : 605820 LUẬN VĂN THẠC SĨ Công trình hồn thành tại: Trường Đại học Bách Khoa – ĐHQG - Tp.HCM Cán hướng dẫn khoa học: TS HỒ ĐỨC DUY Cán chấm nhận xét 1: PGS TS NGUYỄN THỊ HIỀN LƯƠNG Cán chấm nhận xét 2: PGS TS BÙI CÔNG THÀNH Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa – Đại học Quốc gia Tp.HCM ngày 30 tháng năm 2014 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: PGS TS CHU QUỐC THẮNG PGS TS BÙI CÔNG THÀNH PGS TS NGUYỄN THỊ HIỀN LƯƠNG PGS TS NGÔ HỮU CƯỜNG TS HỒ ĐỨC DUY Xác nhận Chủ tịch Hội đồng đánh giá luận văn Trưởng Khoa quản lý chuyên ngành sau luận văn sữa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG ĐẠI HỌC QUỐC GIA TP HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc Lập – Tự Do – Hạnh Phúc -oo - NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: ĐỖ NGỌC TRIỀU MSHV: 11210252 Ngày, tháng, năm sinh: 26/06/1988 Nơi sinh: Cần Thơ Chun ngành: Xây dựng cơng trình DD&CN Mã số ngành: 60.58.20 I TÊN ĐỀ TÀI: ĐÁNH GIÁ MỘT SỐ PHƯƠNG PHÁP CHẨN ĐOÁN HƯ HỎNG KẾT CẤU DẦM TỪ KẾT QUẢ PHÂN TÍCH DAO ĐỘNG II NHIỆM VỤ LUẬN VĂN: Đánh giá so sánh phương pháp chẩn đoán hư hỏng khác kết phân tích dao động Kết so sánh với kết phân tích thực nghiệm từ nghiên cứu tác giả khác - Phương pháp 1: Phương pháp dựa thay đổi tần số - Phương pháp 2: Phương pháp thay đổi dạng dao động - Phương pháp 3: Phương pháp dựa độ cong dạng dao động - Phương pháp 4: Phương pháp dựa thay đổi độ cứng tần số - Phương pháp 5: Phương pháp dựa lượng biến dạng III NGÀY GIAO NHIỆM VỤ: IV NGÀY HOÀN THÀNH NHIỆM VỤ: V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN : TS HỒ ĐỨC DUY Tp Hồ Chí Minh, ngày … tháng … năm 2014 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO TS HỒ ĐỨC DUY TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG LỜI CẢM ƠN Đầu tiên xin chân thành cảm ơn Thầy hướng dẫn T.S Hồ Đức Duy, người tận tình hướng dẫn tơi từ ngày bắt đầu công việc nghiên cứu khoa học lúc hồn thành luận văn Thầy có lời khuyên dẫn tận tình vấn đề đề tài, Thầy động viên nhiều lúc tơi gặp khó khăn Những kiến thức mà Thầy truyền đạt lòng tận tình Thầy tơi ghi nhớ Hình ảnh Thầy gương người Thầy đáng kính nghiệp giáo dục hình mẫu để cố gắng học tập Tôi gởi lời cảm ơn chân thành đến Ban giám hiệu Trường đại học Bách Khoa Thành Phố Hồ Chí Minh, Thầy trực tiếp tham gia giảng dạy truyền đạt kiến thức phương pháp học tập, nghiên cứu Tôi chân thành cảm ơn giúp đỡ anh chị học viên khóa 2011 bạn bè Sau cùng, tơi muốn tỏ lịng biết ơn đến cha mẹ, người thân gia đình hỗ trợ sát cánh bên tơi, động viên tơi suốt q trình học tập làm luận văn trường TP Hồ Chí Minh, tháng năm 2014 Đỗ Ngọc Triều i TÓM TẮT LUẬN VĂN Luận văn nghiên cứu lĩnh vực chẩn đoán hư hỏng kết cấu Đánh giá phương pháp chẩn đoán tốt chẩn đoán hư hỏng dầm đơn giản, từ phân tích tính hiệu tốn thực tế Các phương pháp chẩn đoán dựa kết phân tích dao động Đầu tiên, phương pháp chẩn đoán áp dụng dầm mô phần tử không gian hai chiều, sau mở rộng chẩn đốn dầm mô phần tử khối không gian ba chiều Các phương pháp phân tích sử dụng luận văn sau:  Phương pháp dựa thay đổi tần số  Phương pháp thay đổi dạng dao động  Phương pháp dựa độ cong dạng dao động  Phương pháp dựa thay đổi độ cứng tần số  Phương pháp dựa lượng biến dạng ii LỜI CAM ĐOAN Tôi xin cam đoan luận văn tơi thực hướng dẫn khoa học Thầy TS Hồ Đức Duy Đỗ Ngọc Triều Tháng năm 2014 iii MỤC LỤC Lời cảm ơn i Tóm tắt luận văn ii Lời cam đoan iii Mục lục hình ảnh xi Mục lục bảng xxi Danh mục kí hiệu xxiv Chương GIỚI THIỆU 1.1 Đặt vấn đề 1.2 Mục tiêu luận văn 1.3 Cấu trúc luận văn Chương CƠ SỞ LÝ THUYẾT 2.1 Giới thiệu chung 2.1.1 Trên giới 2.1.2 Tại Việt Nam 2.1.3 Tính cần thiết nghiên cứu 10 2.2 Cơ sở lý thuyết 10 2.2.1 Phương pháp dựa thay đổi tần số (Frequency Change Based Damage Detection Method) 10 2.2.1.1 Giới thiệu phương pháp 10 2.2.1.2 Công thức đánh giá 11 2.2.1.3 Trình tự thực 11 2.2.2 Phương pháp thay đổi dạng dao động (Modal Assurance Criterion) 12 iv 2.2.2.1 Giới thiệu phương pháp 12 2.2.2.2 Công thức đánh giá 13 2.2.2.3 Trình tự thực 14 2.2.3 Phương pháp dựa độ cong dạng dao động (Mode Shape Curvature - Based Damage Detection Method) 15 2.2.3.1 Giới thiệu phương pháp 15 2.2.3.2 Công thức đánh giá 15 2.2.3.3 Trình tự thực 17 2.2.4 Phương pháp dựa thay đổi độ cứng tần số (Stiffness Frequency Change – Based Damage Detection Method) 17 2.2.4.1 Giới thiệu phương pháp 17 2.2.4.2 Công thức đánh giá 18 2.2.4.3 Trình tự thực 22 2.2.5 Phương pháp dựa lượng biến dạng (Modal Strain Energy - Based Damage Detection Method) 22 2.2.5.1 Giới thiệu phương pháp 22 2.2.5.2 Công thức đánh giá 23 2.2.5.3 Trình tự thực 27 Chương CÁC BÀI TOÁN ỨNG DỤNG 28 3.1 Giới thiệu chung 28 3.2 Bài toán - Chẩn đoán dầm mô phần tử 29 3.2.1 Thiết lập toán 29 3.2.1.1 Mô dầm không hư hỏng 29 3.2.1.2 Mô dầm với trường hợp hư hỏng 30 v 3.2.2 Kết toán – Giảm độ cứng 5% vị trí cách gối tựa 0.5 lần chiều dài dầm 32 3.2.2.1 Số liệu ban đầu 32 3.2.2.2 Phương pháp dựa thay đổi tần số (Frequency Change Based Damage Detection Method) 34 3.2.2.3 Phương pháp thay đổi dạng dao động (Modal Assurance Criterion) 35 3.2.2.4 Phương pháp dựa độ cong dạng dao động (Mode Shape Curvature - Based Damage Detection Method) 36 3.2.2.5 Phương pháp dựa thay đổi độ cứng tần số (Stiffness - Frequency Change – Based Damage Detection Method) 38 3.2.2.6 Phương pháp dựa lượng biến dạng (Modal Strain Energy - Based Damage Detection Method) 40 3.2.3 Kết toán – Giảm độ cứng 10% 50% vị trí cách gối tựa 0.5 lần chiều dài dầm 42 3.2.3.1 Số liệu ban đầu 42 3.2.3.2 Phương pháp dựa thay đổi tần số (Frequency Change Based Damage Detection Method) 45 3.2.3.3 Phương pháp thay đổi dạng dao động (Modal Assurance Criterion) 45 3.2.3.4 Phương pháp dựa độ cong dạng dao động (Mode Shape Curvature - Based Damage Detection Method) 47 3.2.3.5 Phương pháp dựa thay đổi độ cứng tần số (Stiffness - Frequency Change – Based Damage Detection Method) 49 3.2.3.6 Phương pháp dựa lượng biến dạng (Modal Strain Energy - Based Damage Detection Method) 52 vi 3.2.4 Kết toán – Giảm độ cứng 5% vị trí cách gối tựa 0.25 lần chiều dài dầm 55 3.2.4.1 Số liệu ban đầu 55 3.2.4.2 Phương pháp dựa thay đổi tần số (Frequency Change Based Damage Detection Method) 57 3.2.4.3 Phương pháp thay đổi dạng dao động (Modal Assurance Criterion) 58 3.2.4.4 Phương pháp dựa độ cong dạng dao động (Mode Shape Curvature - Based Damage Detection Method) 58 3.2.4.5 Phương pháp dựa thay đổi độ cứng tần số (Stiffness - Frequency Change – Based Damage Detection Method) 60 3.2.4.6 Phương pháp dựa lượng biến dạng (Modal Strain Energy - Based Damage Detection Method) 61 3.2.5 Kết toán - Giảm độ cứng 10% 50% vị trí cách gối tựa 0.25 lần chiều dài dầm 63 3.2.5.1 Số liệu ban đầu 63 3.2.5.2 Phương pháp dựa thay đổi tần số (Frequency Change Based Damage Detection Method) 66 3.2.5.3 Phương pháp thay đổi dạng dao động (Modal Assurance Criterion) 66 3.2.5.4 Phương pháp dựa độ cong dạng dao động (Mode Shape Curvature - Based Damage Detection Method) 68 3.2.5.5 Phương pháp dựa thay đổi độ cứng tần số (Stiffness - Frequency Change – Based Damage Detection Method) 70 3.2.5.6 Phương pháp dựa lượng biến dạng (Modal Strain Energy - Based Damage Detection Method) 73 vii Chương 5: Kết luận 5.2 KIẾN NGHỊ Dựa vào phân tích nhận xét trên, tác giả kiến nghị hướng phát triển đề tài:  Phát triển dự đốn mức độ hư hỏng dầm bê tơng cốt thép dầm thép dựa phương pháp chẩn đoán tốt  Dựa vào đánh giá phương pháp tốt cho việc chẩn đoán hư hỏng dầm, ta nên tiếp tục phát triển phương pháp cho dạng kết cấu phức tạp chẩn đoán cho kết cấu khung  Phát triển dự đốn vị trí hư hỏng hệ khung bê tông cốt thép hệ giàn thép  Phát triển dự đoán hư hỏng xảy hệ thép nối bu-long hàn ghép vị trí hàn nối 128 Tài liệu tham khảo TÀI LIỆU THAM KHẢO [1] R D Adams et al, “A Vibration Technique for Nondestructively Assessing the Integrity of Structures,” Journal of Mechanical Engineering Science, vol 20, pp 93–100, Apr.1978 [2] P Cawley, R D Adams “The Location of Defects in Structures from Measurements of Natural Frequencies,” Journal of Strain Analysis, vol.14 No2 1979 [3] H Sato “Free Vibration of Beams with Abrupt Changes of Cross Section,” Journal of Sound and Vibration, 89(1), 59-64 1983 [4] M M F Yuen “A Numerical Study of the Eigenparameters of a Damaged Cantilever,” Journal of Sound and Vibration, 103(3), 301-310 1985 [5] Pandey AK, Biswas M, Samman MM “Damage Detection from Changes in Curvature Mode Shape,” Journal of Sound and Vibration, 145(2):321-32 1991 [6] Y Liang, K Choy, Jialou Hu “Detection of Cracks in Beam Structures Using Measurements of Natural Frequencies,” Journal of the Franklin Institute, 00160032 1991 [7] Ostachowicz WM, Krawczuk M “Analysis of the Effect of Cracks on the Natural Frequencies of a Cantilever Beam,” Journal of Sound and Vibration, 150: 191-201 1991 [8] Jialou Hu, Robert Y Liang “An Integrated Appoach to Detection of Cracks Using Vibration Characteristics,” Department of Civil Engineering, University of Akron, Akron, OH 44325-3905 U.S.A 1993 [9] Norris Stubbs, Jeong-Tae-Kim “Field Verification of a Nondestructive Damage Localization and Severity Estimation Algorithm,” Proceedings of the 13th International Modal Analysis Conference, pp 210-218 1995 129 Tài liệu tham khảo [10] B P Nandwana, S K Maiti “Modelling of Vibration of Beam in Presence of Inclined Edge or Internal Crack for Its Possible Detection Based on Frequency Measurements,” Engineering Fracture Mechanics, vol 58, No 3, pp 193-205 1997 [11] O.S Salawu “Detection of Structural Damage Through Changes in Frequency : a Review,” Engineering Strucures, vol 19, No 9, pp 718-723 1997 [12] E.I Shifrin, R Routolo “Natural Frequencies of a Beam with an Arbitrary Number of Cracks,” Journal of Sound and Vibration, 222(3), 409-423 1999 [13] M N Cerri, F Vestroni “Detection of Damage in Beams Subjected to Diffused Cracking,” Journal of Sound and Vibration, 234(2), 259-276 2000 [14] P Cornwell et al “Application of the Strain Energy Damage Detection Method to Plate Like Structures,” Proceedings of SPIE, the International Society for Optical Engineering, vol 3089, no 2, pp 1312-1318, 1997 [15] Lars Rigner “Modal Assurance Criteria Value for Two Orthogonal Modal Vectors,” S-405 08 Gothenburg Sweden [16] D Fotsch, D J Ewins “Application of MAC in the Frequency Domain”, Imperial College of Science, Technology and Medicine, London SW7 2BX, United Kingdom [17] D.P Patil, S.K Maiti “Detection of Multiple Cracks Using Frequency Measurements,” Engineering Fracture Mechanics, 70 1553-1572 2003 [18] Randall J Allemang “The Modal Assurance Criterion – Twenty Years of Use and Abuse,” Sound and vibration, 2003 [19] Jeong-Tae Kim, Yeon-Sun Ryu, Hyun-Man Cho, Norris Stubbs “Damage Identification in Beam-Type Structures: Frequency-Based Method vs Mode-ShapeBased Method,” Engineering Structures, 25 57-67 2003 130 Tài liệu tham khảo [20] E Douka, G Bamnios, A Trochidis “A method for determining the location and depth of cracks in double – cracked beams,” Science Direct, 65 997-1008 2004 [21] D P Patil, S K Maiti “Experimental Verification of a Method of Detection of Multiple Cracks in Beams Based on Frequency Measurements,” Science Direct, 281 439-451 2005 [22] Y J Yan, L Cheng, Z Y Wu, L H Yam “Development in Vibration-based Structural Damage Detection Technique,” Science Direct, 2198-221 2007 [23] Jialai Wang, Pizhong Qiao “On Irregularity-based Damage Detection Method for Cracked Beams,” Science Direct, 688-704 2008 [24] Le Xuan Hang, Nguyen Thi Hien Luong “Phân Tích Chẩn Đốn Dầm Đàn Hồi Có Nhiều Vết Nứt,” Tạp Chí Phát Triển KH&CN, Tập 12, Số 18 2009 [25] Jeong-Tae Kim, Jae-Hyung Park, Dong-Soo Hong, Woo-Sun Park “Hybrid Health Monitoring of Prestressed Concrete Girder Bridges by Sequential Vibrationimpedance Approaches,” Engineering Structures, 32 115-128 2010 [26] Varinder S Kanwar, Naveen Kwatra, Pankaj Aggarwal, Ramesh P Singh “Use of Vibration Measurements in Health Monitoring of Reinforced Concrete Buildings,” International Journal of Structural Integrity, 1757-9864 2010 [27] A.Dixit, S.Hanagud “Single Beam Analysis of Damage Beams Verified Using a Strain Energy Based Damage Measure,” International Journal of Solids and Structures, 48 592-602 2011 [28] W.Hong, Z.S.Wu, C.Q Yang, C.F.Wan,G Wu, Y.F.Zhang “Condition assessment of reinforced concrete beams using dynamic data measured with distributed long – gage marco – strain sensors,” Journal of Sound and Vibration, vol.331 2764 – 2782 2012 [29] Miroslav Pastor, Michal Binda, Tomas Harcarik “Modal Assurance Criterion,” Sciverse ScienceDirect, 48 543-548 2012 131 Tài liệu tham khảo [30] S.M Seyedpoor “A Two Stage Method for Structural Damage Detection Using a Modal Strain Energy Based Index and Particle Swarm Optimization,” International Journal of Non – Linear Mechanics, 47 1-8 2012 [31] Liheng Wang, Xiyuan Zhou, He Liu, Weiming Yan “Damage detection of RC beams based on experiment and analysis of nonlinear dynamic characteristics,” Construction and Building Materials, 29 420-427 2012 [32] Ming-Hui Hu, Shan-Tung Tu, Fu-Zhen Xuan, Chun-Ming Xia, Hui-He Shao, “Strain Energy Numerical Technique for Structural Damage Detection,” Applied Mathematics and Computation, 219 2424-2431 2012 [33] Duc Duy Ho, Jeong Tae Kim “Prestress – Force Estimation in PSC Girder Using Modal Parameters and System Identification,” Advances in Structural Engineering, vol 15 No.6 2012 [34] Duc Duy Ho, Khac Duy Nguyen, Jeong Tae Kim “Damage Detection in Beam Type Structures Via PZT’s Dual Piezeolectric Responses,” Smart Structures and Systems, vol 11, No 217-240 2013 [35] Sy Dzung Nguyen, Kieu Nhi Ngo, Quang Thinh Tran, Seung Bok Choi “A New Method for Beam – Damage – Diagnosis Using Adaptive Fuzzy Neural Structure and Wavelet Analysis,” Mechanical Systems and Signal Processing, 2013 [36] Model Code 2010 – Final complete draft Volume [37] ABAQUS Analysis User’s Manual, Version 6.10, Dassault Systemes Simulia Corp, Providence, RI, USA 2010 [38] ABAQUS Theory Manual, Version 6.10, Dassault Systems Simulia Corp, Providence, RI, USA 2010 132 Mã nguồn phần mềm MATLAB MÃ NGUỒN PHẦN MỀM MATLAB %% Luan van %% Do Ngoc Trieu %% MSHV: 11210252 clear all; clc; %% 0.0 NHAP THONG SO DAU VAO chieudai= ; %Nhap chieu dai dam (m) chieudai= '; cao= ; %Nhap chieu cao dam (m)= '; rong= ; %'Nhap chieu rong dam (m)= '; Em= ; %Nhap modun dan hoi (N/m2)= '; pro= ; %Nhap khoi luong rieng cua vat lieu (kg/m^3) = ; mmqt=(cao^3*rong)/12; %Tinh mo ment quan tinh cua dam docung=mmqt*Em; %Do cung dientich=cao*rong; %Dien tich mat cat ngang dam ND=1; %So luong hu hong tren dam zo = 1.5; %He so tin cay: 1.5 = 93.3% = 97.7% dx_sp= ; % Do tang cua duong noi suy (m) % Xac dinh vi tri hu hong thuc te – Dung de ve duong chuan vi tri hu hong xd1= ; %Vi tri dau phan tu – hu hong thuc te xd2= ; %Vi tri cuoi phan tu - hu hong thuc te %% 0.0.1 Nhan du lieu tu thu muc excel Chuyen vi cua dam (m) input1=xlsread('filekhonghuhong.xlsx',1); % Load file excel cua dam khong hu hong input2=xlsread(‘filehuhong.xlsx',1); % Load file excel cua dam hu hong %% 0.0.2 Nhan du lieu tu thu muc excel Tan so tu nhien cua dam -(Hz) input3=xlsread('filetansokhonghuhong.xlsx',1); % Tan so dam khong hu hong input4=xlsread('filetansohuhong.xlsx',1); % Tan so dam hu hong %% 0.0.3 Loc du lieu dau vao mode =max(input1(:,5)); % So dang dao dong sodiem=max(input1(:,1)); % So diem chia tren dam li=chieudai/(sodiem-1); % Chieu dai tung phan tu tren dam freq_non=input3(:,3); % Tan so dam khong hu hong (Hz) freq_dam=input4(:,3); % Tan so dam hu hong (Hz) circfreq_non=input3(:,4); % Tan so goc dam khong hu hong (rad/s) circfreq_dam=input4(:,4); % Tan so goc dam hu hong (rad/s) %% CAC DU LIEU %% 1.1.1 Du lieu dau vao tan so tu nhien (Hz) freq=[freq_non,freq_dam] ; % Ma tran tan so [khonghuhong,huhong] [r c]=size(freq); % Loc kich thuoc ma tran tan so dentaq=[]; %% 1.1.2 Tinh lech tan so goc (rad/s) circfreq=[circfreq_non,circfreq_dam]; % Ma tran tan so goc [khonghuhong,huhong] [rw cw]=size(circfreq); % Loc kich thuoc ma tran tan so goc %% 1.2.1 Du lieu dau vao dam khong hu hong m1 = 0; for i1 = 1:mode k1 = 0; m1 = m1 + 1; for j = 1:length(input1(:,1)) if input1(j,5) == m1 k1 = k1 + 1; A(k1,i1) = input1(j,8); %A: matran chuyenvi dam khong hu hong end end end for h=1:mode if A(2,h)>0 A(:,h)=A(:,h)*(-1); else A(:,h)=A(:,h); %A: matran chuyenvi dam khong hu hong end end %% 1.2.2 Du lieu dau vao dam hu hong m2 = 0; for i2 = 1:mode k2 = 0; m2 = m2 + 1; for j2 = 1:length(input2(:,1)) 133 Mã nguồn phần mềm MATLAB if input2(j2,5) == m2 k2 = k2 + 1; D(k2,i2) = input2(j2,8); %D: matran chuyenvi dam hu hong end end end for k=1:mode if D(2,k)>0 D(:,k)=D(:,k)*(-1); else D(:,k)=D(:,k); %D: matran chuyenvi dam hu hong end end %% 1.2.3 Chuan hoa chuyen vi dam khong hu hong va dam hu hong for i8 = 1:mode for j8 = 1:sodiem norm_non(j8,i8) = A(j8,i8)/norm(A(:,i8)); %Chuan hoa KO HUHONG norm_dam(j8,i8) = D(j8,i8)/norm(D(:,i8)); %Chuan hoa HUHONG end end %% 1.3 Lap ma tran li – ma tran dai phan tu k3=0; for i3=0:sodiem-1 k3=k3+1; ll=li.*i3; L(k3,1)=ll; %L: matran dai end %% 1.4.1 Dao ham bac II - Bien doi duong cong dang dao dong – dam khong hu hong for i5=1:mode for j5=2:sodiem-1 c1=(A(j5+1,i5)-2.*A(j5,i5)+A(j5-1,i5))./(li.^2); curv_non(j5,i5)=c1; % curv_non: ma tran duong cong dang dao dong – dam khong hu hong end curv_non(1,i5)=0; % Gan dieu kien bien - dam don gian curv_non(sodiem,i5)=0; % Gan dieu kien bien - dam don gian end %% 1.4.2 Dao ham bac II - Bien doi duong cong dang dao dong – dam hu hong for i6=1:mode for j6=2:sodiem-1 c2=(D(j6+1,i6)-2.*D(j6,i6)+D(j6-1,i6))./(li.^2); curv_dam(j6,i6)=c2; % curv_dam: ma tran duong cong dang dao dong – dam hu hong end curv_dam(1,i6)=0; % Gan dieu kien bien curv_dam(sodiem,i6)=0; % Gan dieu kien bien end %% 1.5 Chuan hoa duong cong dang dao dong for i7 = 1:mode for j7 = 1:sodiem % Chuan hoa duong cong dang dao dong cho truong hop dam khong hu hong curvnon_norm(j7,i7) = curv_non(j7,i7)/norm(curv_non(:,i7)); % Chuan hoa duong cong dang dao dong cho truong hop dam hu hong curvdam_norm(j7,i7) = curv_dam(j7,i7)/norm( curv_dam(:,i7)); end end %% 1.6 Noi suy duong cong dang dao dong xx =0:dx_sp:chieudai; temp1 = []; temp2 = []; for is = 1:mode y1 = spline(L,curvnon_norm(:,is),xx); temp1(:,is) = y1; y2 = spline(L,curvdam_norm(:,is),xx); temp2(:,is) = y2; end curvnon_sp = temp1; %Noi suy duong cong dam khong hu hong curvdam_sp = temp2; %Noi suy duong cong dam hu hong [node_sp cot] = size(curvnon_sp(:,1)); %% 2.2 CAC PHUONG PHAP CHAN DOAN %% 2.2.1 PHUONG PHAP DUA TREN SU THAY DOI TAN SO (FREQUENCY CHANGE BASED DAMAGE DETECTION METHOD) 134 Mã nguồn phần mềm MATLAB % Tinh phan tram lech tan so for i=1:r dentaq(i,1)=(freq(i,1)-freq(i,2)); percen_q(i,1)=(freq(i,1)-freq(i,2))/freq(i,1)*100; %Phan tram lech tan so denta_q2(i,1)=(freq(i,1)^2-freq(i,2)^2)/(freq(i,1)^2); end % Tinh phan tram lech tan so goc dentaw=[]; for iw=1:rw dentaw(iw,1)=(circfreq(iw,1)-circfreq(iw,2)); %Phan tram lech tan so goc denta_w2(iw,1)=(circfreq(iw,1)^2-circfreq(iw,2)^2)/(circfreq(iw,1)^2); end %% 2.2.2 PHUONG PHAP DO CONG DANG DAO DONG (MODAL ASSURANCE CRITERION) %%2.2.2.1 Tinh gia tri so M.A.C cho tuong quan dam hu hong – dam khong hu hong tAX3=[]; mA3=[]; mX3=[]; for n3=1:mode tAX3(1,n3)=((norm_dam(:,n3))'*norm_non(:,n3)); mA3(1,n3)=(norm_dam(:,n3))'*norm_dam(:,n3); mX3(1,n3)=(norm_non(:,n3))'*norm_non(:,n3); end MAC_AX3=(sum(tAX3(1,:)))^2/(sum(mA3(1,:))*sum(mX3(1,:))); %Gia tri tuong quan hai dam %2.2.2.2 Tinh ma tran so M.A.C cho cung mot mo hinh dam – dam khong hu hong for ii1=1:5 for jj1=1:5 t3=(norm_non(:,ii1)'*norm_non(:,jj1))^2; ms5=norm_non(:,jj1)'*norm_non(:,jj1); ms6=norm_non(:,ii1)'*norm_non(:,ii1); mac3(ii1,jj1)=t3/(ms5*ms6); end end mac3; %Ma tran so M.A.C cho tuong quan cac dang dao dong – dam khong hu hong %2.2.2.3 Tinh ma tran so M.A.C cho tuong quan giua dam khong hu hong – dam hu hong for ib=1:5 for jb=1:5 t4=(norm_dam(:,ib)'*norm_non(:,jb))^2; ms7=norm_non(:,jb)'*norm_non(:,jb); ms8=norm_dam(:,ib)'*norm_dam(:,ib); mac4(ib,jb)=t4/(ms7*ms8); end end mac4; %Ma tran so M.A.C cho tuong quan dam khong hu hong – dam hu hong %% 2.2.2.3 PHUONG PHAP DO CONG DANG DAO DONG (MODE SHAPE CUVARTURE - BASED DAMAGE DETECTION METHOD) differ=abs(abs(curv_dam)-abs(curv_non)); %Ham su khac biet giua cac dang dao dong %% 2.2.2.4 PHUONG PHAP DUA TREN SU THAY DOI DO CUNG VA TAN SO (STIFFNESS – FREQUENCY CHANGE – BASED DAMAGE DETECTION METHOD) f1=percen_q(1,1); %tao giam tan so - dang dao dong f2=percen_q(2,1); %tao giam tan so - dang dao dong f3=percen_q(3,1); %tao giam tan so - dang dao dong f4=percen_q(4,1); %tao giam tan so - dang dao dong f5=percen_q(5,1); %tao giam tan so - dang dao dong x_qsin=[0.01:0.01:1]; Kr1 =(docung.*((sin(1.*pi.*x_qsin)).^2))./(chieudai.*f1); %Gia tri cung phan tu – dang dao dong Kr2 =(docung.*((sin(2.*pi.*x_qsin)).^2))./(chieudai.*f2); %Gia tri cung phan tu – dang dao dong Kr3 =(docung.*((sin(3.*pi.*x_qsin)).^2))./(chieudai.*f3); %Gia tri cung phan tu – dang dao dong Kr4 =(docung.*((sin(4.*pi.*x_qsin)).^2))./(chieudai.*f4); %Gia tri cung phan tu – dang dao dong Kr5 =(docung.*((sin(5.*pi.*x_qsin)).^2))./(chieudai.*f5); %Gia tri cung phan tu – dang dao dong %% 2.2.2.5 PHUONG PHAP NANG LUONG BIEN DANG (MODAL STRAIN ENERGY –BASED DAMAGE DETECTION METHOD) %% 2.2.2.5.1 Gia tri chi so vung hu hong mode1_sp=[curvnon_sp(:,1) curvdam_sp(:,1)]; mode2_sp=[curvnon_sp(:,2) curvdam_sp(:,2)]; mode3_sp=[curvnon_sp(:,3) curvdam_sp(:,3)]; mode4_sp=[curvnon_sp(:,4) curvdam_sp(:,4)]; mode5_sp=[curvnon_sp(:,5) curvdam_sp(:,5)]; [node no_data] = size(mode1_sp); for i = 1:no_data for j = 1:node_sp-1 135 Mã nguồn phần mềm MATLAB MSE1(j,i) = (mode1_sp(j,i)^2 + mode1_sp(j+1,i)^2)*dx_sp/2; MSE2(j,i) = (mode2_sp(j,i)^2 + mode2_sp(j+1,i)^2)*dx_sp/2; MSE3(j,i) = (mode3_sp(j,i)^2 + mode3_sp(j+1,i)^2)*dx_sp/2; MSE4(j,i) = (mode4_sp(j,i)^2 + mode4_sp(j+1,i)^2)*dx_sp/2; MSE5(j,i) = (mode5_sp(j,i)^2 + mode5_sp(j+1,i)^2)*dx_sp/2; xxx(j) = (j - 0.5)*dx_sp; end end for i = 1:no_data MSE1(:,i) = MSE1(:,i) / sum(MSE1(:,i)); MSE2(:,i) = MSE2(:,i) / sum(MSE2(:,i)); MSE3(:,i) = MSE3(:,i) / sum(MSE3(:,i)); MSE4(:,i) = MSE4(:,i) / sum(MSE4(:,i)); MSE5(:,i) = MSE5(:,i) / sum(MSE5(:,i)); end %% 2.2.2.5.1 Gia tri hu hai Zi for i = 2:no_data for j = 1:node_sp-1 beta1(j,i-1) = (MSE1(j,i)+1) / ((sum(MSE1(:,1)/ND*denta_q2(1)) + MSE1(j,1))+1); beta2(j,i-1) = (MSE2(j,i)+1) / ((sum(MSE2(:,1)/ND*denta_q2(2)) + MSE2(j,1))+1); beta3(j,i-1) = (MSE3(j,i)+1) / ((sum(MSE3(:,1)/ND*denta_q2(3)) + MSE3(j,1))+1); beta4(j,i-1) = (MSE4(j,i)+1) / ((sum(MSE4(:,1)/ND*denta_q2(4)) + MSE4(j,1))+1); beta5(j,i-1) = (MSE5(j,i)+1) / ((sum(MSE5(:,1)/ND*denta_q2(5)) + MSE5(j,1))+1); beta6(j,i-1) = ((MSE1(j,i)+MSE2(j,i))+1)/(((sum(MSE1(:,1)/ND*denta_q2(1)) + MSE1(j,1))+(sum(MSE2(:,1)/ND*denta_q2(2)) + MSE2(j,1)))+1); beta7(j,i-1) = ((MSE1(j,i)+MSE3(j,i))+1)/(((sum(MSE1(:,1)/ND*denta_q2(1)) + MSE1(j,1))+(sum(MSE3(:,1)/ND*denta_q2(3)) + MSE3(j,1)))+1); beta8(j,i-1) = ((MSE2(j,i)+MSE3(j,i))+1)/(((sum(MSE2(:,1)/ND*denta_q2(2)) + MSE2(j,1))+(sum(MSE3(:,1)/ND*denta_q2(3)) + MSE3(j,1)))+1); beta9(j,i-1) = ((MSE1(j,i)+MSE1(j,i)+MSE3(j,i))+1)/(((sum(MSE1(:,1)/ND*denta_q2(1)) + MSE1(j,1))+(sum(MSE2(:,1)/ND*denta_q2(2)) + MSE2(j,1))+(sum(MSE3(:,1)/ND*denta_q2(3)) + MSE3(j,1)))+1); end Z1(:,i-1) = (beta1(:,i-1) - mean(beta1(:,i-1))) / std(beta1(:,i-1)); Z2(:,i-1) = (beta2(:,i-1) - mean(beta2(:,i-1))) / std(beta2(:,i-1)); Z3(:,i-1) = (beta3(:,i-1) - mean(beta3(:,i-1))) / std(beta3(:,i-1)); Z4(:,i-1) = (beta4(:,i-1) - mean(beta4(:,i-1))) / std(beta4(:,i-1)); Z5(:,i-1) = (beta5(:,i-1) - mean(beta5(:,i-1))) / std(beta5(:,i-1)); Z6(:,i-1) = (beta6(:,i-1) - mean(beta6(:,i-1))) / std(beta6(:,i-1)); end % Note: Z1 = Mode1, Z2 = Mode2, Z3 = Mode3,Z4 = Mode4,Z5 = Mode5, Z6 = Mode1 + Mode2 for i = 1:length(xxx) Zo(i) = zo; %Vi tri hu hai gia dinh end %% 3.0 VE DO THI disp('NHAP SO DE CHON PHUONG PHAP’); disp('0-Ve cac dang dao dong'); disp('1-Phuong phap su thay doi tan so'); disp('2-Phuong phap thay doi dang dao dong'); disp('3-Phuong phap cong dang dao dong'); disp('4-Phuong phap su thay doi cung va tan so'); disp('5-Phuong phap nang luong bien dang'); so_tieuchuan=input('Nhap so = ','s'); switch so_tieuchuan %% 3.0.1 Ve cac dang dao dong % 3.0.1 Dang dao dong % Ve dang dao dong1 %%%%%%%%%%%%%% case '0' figure (1); plot(L,A(:,1),'b-'); hold on plot(L,D(:,1),'k-.','LineWidth',2); grid on title(' Mode1 ','fontsize',13); xlabel('Dimension-m','fontsize',13); ylabel('Displacement','fontsize',13); set(gcf, 'color', 'white'); legend('Intact','Damage'); grid on % 3.0.2 Dang dao dong % Ve dang dao dong %%%%%%%%%%%%%% figure (2); plot(L,A(:,2),'b-'); hold on plot(L,D(:,2),'k-.','LineWidth',2); grid on 136 Mã nguồn phần mềm MATLAB title(' Mode2 ','fontsize',13); xlabel('Dimension-m','fontsize',13); ylabel('Displacement','fontsize',13); set(gcf, 'color', 'white'); legend('Intact','Damage'); grid on % 3.0.3 Dang dao dong % Ve dang dao dong %%%%%%%%%%%%%% figure (3); plot(L,A(:,3),'b-'); hold on plot(L,D(:,3),'k-.','LineWidth',2); grid on title(' Mode3 ','fontsize',13); xlabel('Dimension-m','fontsize',13); ylabel('Displacement','fontsize',13); set(gcf, 'color', 'white'); legend('Intact','Damage'); grid on % 3.0.4 Dang dao dong % Ve dang dao dong %%%%%%%%%%%%%% figure (4); plot(L,A(:,4),'b-'); hold on plot(L,D(:,4),'k-.','LineWidth',2); grid on title(' Mode4 ','fontsize',13); xlabel('Dimension-m','fontsize',13); ylabel('Displacement','fontsize',13); set(gcf, 'color', 'white'); legend('Intact','Damage'); grid on % 3.0.5 Dang dao dong % Ve dang dao dong %%%%%%%%%%%%%% figure (5); plot(L,A(:,5),'b-'); hold on plot(L,D(:,5),'k-.','LineWidth',2); grid on title(' Mode5 ','fontsize',13); xlabel('Dimension-m','fontsize',13); ylabel('Displacement','fontsize',13); set(gcf, 'color', 'white'); legend('Intact','Damage'); grid on %% 3.1 Ve phan tram lech tan so – Phuong phap dua tren su thay doi tan so case '1' figure(1); freq(:,1); x_q=1:5; plot(x_q,freq(1:5,1),'b-',x_q,freq(1:5,2),'k-.'); legend('Intact','Damage'); title('Criterion Frequency','fontsize',13); xlabel('Mode','fontsize',13); ylabel('Frequency','fontsize',13); set(gcf, 'color', 'white'); h=gca; %Gan cac ki hieu vao vi tri mong muon set(h,'XTick',1:0.5:5); set(h,'XTickLabel',{'1','','2','','3','','4','','5'}); %Gan ten cho cac ki hieu grid on figure(2); dolech=[]; for dl=1:5 dolech(dl,1)=freq(dl,1)-freq(dl,2); end bar(dolech,0.2); hold on set(gcf, 'color', 'white'); title('Frequency Deflection','fontsize',13); xlabel('Mode','fontsize',13); ylabel('Frequency','fontsize',13); grid on figure(3); dolechts=[]; for dl=1:5 dolechts(dl,1)=((freq(dl,1)-freq(dl,2))/freq(dl,1))*100; end bar(dolechts,0.2); hold on 137 Mã nguồn phần mềm MATLAB set(gcf, 'color', 'white'); title('% Frequency Deflection','fontsize',13); xlabel('Mode','fontsize',13); ylabel('% Frequency','fontsize',13); grid on %% 3.2 Xuat gia tri so M.A.C va gia tri ma tran M.A.C - Phuong phap thay doi dang dao dong case'2' disp('Gia tri tuong dong giua mo hinh dam khong hu hong – dam hu hong); MAC_AX3 disp('Ma tran so M.A.C giua cac dang dao dong – dam khong hu hong'); mac3 disp('Ma tran so M.A.C giua mo hinh dam khong hu hong – dam hu hong'); mac4 %% 3.3 Ve thi phuong phap cong dang dao dong case '3' %% 3.3.1 Do thi ba chieu matran0=ones(sodiem,mode); figure(1); plot3(L,0.*matran0(:,1),differ(:,1),'k-','LineWidth',1.5); hold on%mode1 plot3(L,1.*matran0(:,1),differ(:,2),'b-','LineWidth',1.5); hold on%mode2 plot3(L,2.*matran0(:,1),differ(:,3),'r-','LineWidth',1.5); hold on%mode3 plot3(L,3.*matran0(:,1),differ(:,4),'k-','LineWidth',1.5); hold on%mode4 plot3(L,4.*matran0(:,1),differ(:,5),'b-','LineWidth',1.5); hold on%mode5 % Gan diem hu hong max1=max(differ(:,1)); max2=max(differ(:,2)); max3=max(differ(:,3)); max4=max(differ(:,4)); max5=max(differ(:,5)); plot3(4.5,0:0.1:4,0,'k-*','LineWidth',0.2); hold on plot3(5,0:0.1:4,0,'k-*','LineWidth',0.2); hold on %Thiet lap lai truc toa set(gcf, 'color', 'white'); h=gca; %Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'});%Gan ten cho cac ki hieu %Thiet lap truc Y set(h,'YTick',[0,1,2,3,4,5]); %Vi tri dat ki hieu set(h,'YTickLabel',{'Mode1','Mode2','Mode3','Mode4','Mode5',''}); %Gan ten cho cac ki hieu %Gan ten truc tren thi xlabel('Point number','fontsize',13), ylabel('Mode','fontsize',13),zlabel('Difference','fontsize',13); title('Criterion Curvature Mode Shape','fontsize',13); grid on %% 3.3.2 Do thi hai chieu cho cac dang dao dong %% figure(2); %Hinh dang dao dong plot(L,differ(:,1),'k-','linewidth',1.5); hold on plot([xd1 xd2],[0+max(differ(:,1))/50 0+max(differ(:,1))/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); h=gca; % Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'});%Gan ten cho cac ki hieu % Gan ten truc X-Y xlabel('Point number','fontsize',13); ylabel('Difference','fontsize',13); title('Criterion Curvature Mode Shape Mode 1','fontsize',13); grid on %% figure(3); %Hinh dang dao dong plot(L,differ(:,2),'b-','linewidth',1.5); hold on plot([xd1 xd2],[0+max(differ(:,2))/50 0+max(differ(:,2))/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); h=gca; % Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'}); % Gan ten truc X-Y xlabel('Point number','fontsize',13); ylabel('Difference','fontsize',13); 138 Mã nguồn phần mềm MATLAB title('Criterion Curvature Mode Shape Mode 2','fontsize',13); grid on %% figure(4); %Hinh dang dao dong plot(L,differ(:,3),'r-*','linewidth',1.5); hold on plot([xd1 xd2],[0+max(differ(:,3))/50 0+max(differ(:,3))/50],'-k','LineWidth',3); hold on h=gca; set(gcf, 'color', 'white'); % Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'});%Gan ten cho cac ki hieu % Gan ten truc X-Y xlabel('Point number','fontsize',13); ylabel('Difference','fontsize',13); title('Criterion Curvature Mode Shape Mode 3','fontsize',13); grid on %% figure(5); %Hinh dang dao dong plot(L,differ(:,4),'k-','linewidth',1.5); hold on plot([xd1 xd2],[0+max(differ(:,4))/50 0+max(differ(:,4))/50],'-k','LineWidth',3); hold on h=gca; set(gcf, 'color', 'white'); % Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'});%Gan ten cho cac ki hieu % Gan ten truc X-Y xlabel('Point number','fontsize',13); ylabel('Difference','fontsize',13); title('Criterion Curvature Mode Shape Mode 4','fontsize',13); grid on %% figure(6); %Hinh dang dao dong plot(L,differ(:,5),'b-','linewidth',1.5); hold on plot([xd1 xd2],[0+max(differ(:,5))/50 0+max(differ(:,5))/50],'-k','LineWidth',3); hold on h=gca; set(gcf, 'color', 'white'); % Thiet lap truc X set(h,'XTick',L); %Vi tri dat ki hieu set(h,'XTickLabel',{'1','','3','','5','','7','','9','','11','','13','','15','','17','','19','','21'});%Gan ten cho cac ki hieu % Gan ten truc X-Y xlabel('Point number','fontsize',13); ylabel('Difference','fontsize',13); title('Criterion Curvature Mode Shape Mode 5','fontsize',13); grid on %% 3.4 Ve thi phuong phap cong dang dao dong case '4' % Ve figure(1); plot(x_qsin,Kr1,'b') ;hold on plot(x_qsin,Kr2,'r') ;hold on plot(x_qsin,Kr3,'k') ;hold on plot(x_qsin,Kr4,'b') ;hold on plot(x_qsin,Kr5,'b-.','LineWidth',2) ;hold on % Ve vi tri hu hai gia dinh – Thiet lap truc X-Y plot([xd1*0.1 xd2*0.1],[0+max(Kr1)/50 0+max(Kr1)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); set(gca,'fontsize',13); title('Frequency Criterion'); xlabel('Li/L');ylabel('Kr'); set(gca,'XTick',0:0.1:1); %Vi tri dat ki hieu set(gca,'XTickLabel',{'0','','0.2','','0.4','','0.6','','0.8','','1',}); %Gan ten cho cac ki hieu legend('Mode1','Mode2','Mode3','Mode4','Mode5'); grid on %% 3.5 Ve thi phuong phap nang luong bien dang case '6' % Ve dang dao dong figure(1); % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 jz1=0; for i=1:length(Z1) if (Z1(i,1)-1.5)>10^-2 139 Mã nguồn phần mềm MATLAB jz1=jz1+1; xz1(jz1,1)=xxx(i); end end x1_min=min(xz1); x1_max=max(xz1); % Ve lim_y = [-2 5]; lim_y = [-2 5]; text1='Z1 (Mode 1)'; set(gca,'fontsize',13); plot(xxx,Z1,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z1)/50 min(lim_y)+max(Z1)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text1); hold on plot([x1_min x1_min], lim_y,'-black',[x1_max x1_max], lim_y,'-black'); hold off % Ve dang dao dong figure(2); % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 jz2=0; for i=1:length(Z2) if (Z2(i,1)-1.5)>10^-2 jz2=jz2+1; xz2(jz2,1)=xxx(i); end end x2_min=min(xz2); x2_max=max(xz2); % Ve lim_y = [-2 5]; text2=' Z2 (Mode 2)'; set(gca,'fontsize',13); plot(xxx,Z2,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z2)/50 min(lim_y)+max(Z2)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text2); hold on plot([x2_min x2_min], lim_y,'-black',[x2_max x2_max], lim_y,'-black'); hold off % Ve dang dao dong figure(3); % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 jz3=0; for i=1:length(Z3) if (Z3(i,1)-1.5)>10^-2 jz3=jz3+1; xz3(jz3,1)=xxx(i); end end x3_min=min(xz3); x3_max=max(xz3); % Ve lim_y = [-2 5]; text3=' Z3 (Mode 3)'; set(gca,'fontsize',13); plot(xxx,Z3,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z3)/50 min(lim_y)+max(Z3)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text3); hold on plot([x3_min x3_min], lim_y,'-black',[x3_max x3_max], lim_y,'-black'); hold off % Ve dang dao dong % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 figure(4); jz4=0; 140 Mã nguồn phần mềm MATLAB for i=1:length(Z4) if (Z4(i,1)-1.5)>10^-2 jz4=jz4+1; xz4(jz4,1)=xxx(i); end end x4_min=min(xz4); x4_max=max(xz4); % Ve lim_y = [-3 5]; text4=' Z4 (Mode 4)'; set(gca,'fontsize',13); plot(xxx,Z4,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z4)/50 min(lim_y)+max(Z4)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text4); hold on plot([x4_min x4_min], lim_y,'-black',[x4_max x4_max], lim_y,'-black'); hold off % Ve dang dao dong figure(5); % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 jz5=0; for i=1:length(Z5) if (Z5(i,1)-1.5)>10^(-2) jz5=jz5+1; xz5(jz5,1)=xxx(i); end end x5_min=min(xz5); x5_max=max(xz5); % Ve lim_y = [-2 5]; text5=' Z5 (Mode 5)'; set(gca,'fontsize',13); plot(xxx,Z5,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z5)/50 min(lim_y)+max(Z5)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text5); hold on plot([x5_min x5_min], lim_y,'-black',[x5_max x5_max], lim_y,'-black'); hold off % Ve dang dao dong figure(6); % Tinh vi tri giao cua duong cong va duong thang di qua truc tung gia tri 1.5 jz6=0; for i=1:length(Z6) if (Z6(i,1)-1.5)>10^-2 jz6=jz6+1; xz6(jz6,1)=xxx(i); end end x6_min=min(xz6); x6_max=max(xz6); % Ve text6=' Z6 (Mode1+Mode2)'; set(gca,'fontsize',13); plot(xxx,Z6,'-b',xxx,Zo,' r'); hold on plot([xd1 xd2],[min(lim_y)+max(Z6)/50 min(lim_y)+max(Z6)/50],'-k','LineWidth',3); hold on set(gcf, 'color', 'white'); xlabel('Length (m)');ylabel('Normalized Damage Index'); title(text6); hold on plot([x6_min x6_min], lim_y,'-black',[x6_max x6_max], lim_y,'-black'); hold off end %% ket thuc ma nguon MATLAB 141 LÍ LỊCH TRÍCH NGANG Họ tên: ĐỖ NGỌC TRIỀU Phái: Nam Ngày sinh: 26 tháng năm 1988 Nơi sinh: Cần Thơ Địa liên lạc: 98/46 Trần Hưng Đạo, Q Ninh Kiều, TP Cần Thơ Điện thoại: 0978 228 298 Email: ngoctrieudo@gmail.com QUÁ TRÌNH ĐÀO TẠO Văn Chuyên ngành đào tạo Thời gian Năm Nơi đào tạo đào tạo Kỹ sư Xây dựng dân dụng công nghiệp 4.5 năm 2006-2011 Trường Đại Học Cần Thơ Thạc sĩ Xây dựng dân dụng công nghiệp năm 2011-2014 Trường ĐH Bách Khoa TP.HCM Q TRÌNH CƠNG TÁC Thời gian Từ 09/2013 07/2014 Tên công ty Đến 07/2014 Cơng ty CP Cơ Khí Điện Máy Camemco Cần Thơ Sở Tài Chính – TP Cần Thơ Bộ Cơng việc phận giao Phịng Kỹ thuật Phịng Đầu Tư Người tham chiếu PGĐ Đỗ Thành Quận Kỹ sư (0913617353) Kỹ sư PTP Lê Thị Hoàng Yến (0919811464) ... Mã số ngành: 60.58.20 I TÊN ĐỀ TÀI: ĐÁNH GIÁ MỘT SỐ PHƯƠNG PHÁP CHẨN ĐOÁN HƯ HỎNG KẾT CẤU DẦM TỪ KẾT QUẢ PHÂN TÍCH DAO ĐỘNG II NHIỆM VỤ LUẬN VĂN: Đánh giá so sánh phương pháp chẩn đoán hư hỏng. .. lĩnh vực chẩn đoán hư hỏng kết cấu Đánh giá phương pháp chẩn đoán tốt chẩn đoán hư hỏng dầm đơn giản, từ phân tích tính hiệu toán thực tế Các phương pháp chẩn đoán dựa kết phân tích dao động Đầu... Thu giá trị tần số dạng dao động, từ kết phân tích dao động dầm hư hỏng dầm không hư hỏng phần mềm phần tử hữu hạn  Tính giá trị hiệu tần số ( f n ) dầm không hư hỏng dầm hư hỏng dạng dao động