Độ nhạy cảm của độ võng tấm có vết nứt bằng phương pháp biến đổi wavelet

ĐẠI HỌC QUỐC QIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA  PHẠM NGỌC TIẾN ĐỘ NHẠY CẢM CỦA ĐỘ VÕNG TẤM CÓ VẾT NỨT BẰNG PHƯƠNG PHÁP BIẾN ĐỔI WAVELET CHUYÊN NGÀNH: XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP MÃ SỐ NGÀNH: 60.58.20 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH THÁNG 07 NĂM 2007 LUẬN VĂN ĐƯỢC HÒAN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH Nhận xét cán hướng dẫn: PGS.TS NGUYỄN THỊ HIỀN LƯƠNG Nhận xét cán phản biện 1: ……………………………………… Nhận xét cán phản biện 2: ……………………………………… Luận văn thạc sĩ bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ Trường Đại học Bách Khoa, ngày… tháng… năm 2007 Có thể tìm hiểu luận văn thư viện trường Đại học Bách Khoa 608TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH Độc lập – Tự – Hạnh phúc - Tp Hồ Chí Minh, ngày …tháng …năm 2007 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên: PHẠM NGỌC TIẾN Phái: Nam Ngày, tháng, năm sinh: 21-08-1977 Nơi sinh: Bình Định Khóa: 2005 Mã học viên: 02105521 Chuyên ngành: XDDD&CN Mã số ngành: 60.58.20 I TÊN ĐỀ TÀI ĐỘ NHẠY CẢM CỦA ĐỘ VÕNG TẤM CÓ VẾT NỨT BẰNG PHƯƠNG PHÁP BIẾN ĐỔI WAVELET II NHIỆM VỤ VÀ NỘI DUNG Luận văn khảo sát tính nhạy cảm độ võng chịu uốn có vết nứt thông qua kỹ thuật biến đổi wavelet Việc giải tốn để tìm giá trị độ võng cho thực phương pháp phần tử hữu hạn kết hợp với mơ hình phần tử đẳng tham số tứ giác bậc hai nút mô hình phần tử đẳng tham số tam giác Barsoum (nhằm mơ tả tính kỳ dị đáy vết nứt) Các kết giải kiểm chứng phần mềm Ansys Sau đó, giá trị độ võng phân tích lại cơng cụ wavelet Từ rút số nhận xét kiến nghị III NGÀY GIAO NHIỆM VỤ: 05-02-2007 IV NGÀY HOÀN THÀNH: 05-07-2007 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS.TS NGUYỄN THỊ HIỀN LƯƠNG CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ NGÀNH PGS.TS Nguyễn Thị Hiền Lương Nội dung đề cương luận văn thạc sĩ Hội Đồng Chun Ngành thơng qua Ngày…….tháng…….năm 2007 PHỊNG ĐÀO TẠO SĐH KHOA QUẢN LÝ NGÀNH LỜI CẢM ƠN Trong năm theo học chương trình đào tạo Thạc sỹ Trường Đại học Bách Khoa Tp Hồ Chí Minh, cuối vượt qua chặng đường sau hồn tất luận văn tốt nghiệp Có kết này, trước tiên xin chân thành cảm ơn Cô PGS.TS NGUYỄN THỊ HIỀN LƯƠNG tận tình hướng dẫn, uốn nắn, động viên tơi thời gian thực đề tài luận văn Cô giúp tác giả hình thành nên phong cách làm việc khoa học, hướng dẫn tác giả bước đường nghiên cứu khoa học Xin gởi đến Cô lời biết ơn sâu sắc Tôi xin chân thành gởi lời cảm ơn đến q thầy phụ trách chương trình đào tạo Thạc sỹ truyền đạt kiến thức quí báu, cung cấp nguồn tài liệu kịp thời Đó hành trang tơi đường đời sau Tôi xin cảm ơn bạn bè, đồng ngiệp, người thân động viên giúp đỡ tơi suốt q trình thực luận văn Cuối cùng, tơi xin bày tỏ lịng biết ơn vơ hạn đến với cha mẹ anh em gia đình, người ln quan tâm, động viên, chia sẻ với tác giả mặt sống Xin chân thành cảm ơn! TÓM TẮT LUẬN VĂN TÓM TẮT Các khuyết tật, đặc biệt vết nứt, có ảnh hưởng lớn đến khả làm việc tuổi thọ cho cơng trình Trong số cơng trình nghiên cứu có liên quan, việc phân tích ảnh hưởng vết nứt dựa thay đổi đặc trưng động học thu kết hữu ích việc chẩn đốn tình trạng làm việc cho kết cấu Tuy nhiên, trình vận dụng đặc trưng tồn số hạn chế định Ở khía cạnh khác, luận văn ý đến việc khảo sát kết cấu chịu uốn có vết nứt liên quan đến độ võng, thông qua phép biến đổi wavelet Việc giải tốn để tìm giá trị độ võng cho tiến hành sở phương pháp Phần tử hữu hạn dựa mơ hình phần tử đẳng tham số tứ giác bậc hai nút mơ hình phần tử đẳng tham số tam giác Barsoum (nhằm mô tả tính kỳ dị đáy vết nứt) Các kết giải kiểm chứng phần mềm Ansys Sau đó, giá trị độ võng biến đổi lại cơng cụ wavelet để phân tích độ nhạy cảm độ võng chiều dài, vị trí vết nứt Để hỗ trợ cho cơng việc khảo sát nêu trên, tác giả nghiên cứu, giải vấn đề có liên quan trình bày chương luận văn sau: Chương 1: Tổng quan Trình bày sơ lược khuyết tật, đặc biệt vết nứt kết cấu Quá trình phát triển tốn Chẩn đốn kỹ thuật, lịch sử phát triển toán Tấm chịu uốn, lịch sử phát triển phương pháp Phần tử hữu hạn lịch sử phát triển Wavelet Trên sở nêu lên ý nghĩa khoa học thực tiễn đề tài Chương 2: Lý thuyết học phá huỷ, lý thuyết chịu uốn phương pháp PTHH Trình bày sơ lược kiến thức học phá huỷ đặc trưng đàn hồi, hình thành biến dạng dẻo đáy vết nứt, từ lựa chọn mơ hình phần tử mơ tả làm việc cho vùng lân cận đáy vết nứt Các sở lý thuyết toán chịu uốn theo lý thuyết cổ điển Kirchhoff lý thuyết Reissner-Mindlin, sở lý thuyết cho việc thiết lập ma trận độ cứng vectơ tải cho tốn chịu uốn theo mơ hình phần tử đẳng tham số bậc hai dựa lý thuyết Reissner-Mindlin Chương 3: Phương pháp biến đổi wavelet Trình bày sở lý thuyết wavelet, sở phân tích wavelet, q trình biến đổi tái tạo dùng phân tích tín hiệu chiều, khả áp dụng hộp công cụ wavelet phần mềm Matlab dùng cho phân tích tín hiệu Chương 4: Xác định độ võng phương pháp PTHH TĨM TẮT LUẬN VĂN Trình bày tốn cụ thể áp dụng phương pháp PTHH vừa thiết lập để tìm tập hợp giá trị độ võng cho kết cấu Kiểm chứng kết giải phần mềm Ansys Chương 5: Biến đổi wavelet cho độ võng Trên sở số liệu độ võng giải chương 4, dùng phép biến đổi wavelet để khảo sát độ nhạy cảm độ võng cho số trường hợp cụ thể Chương 6: Kết luận kiến nghị Từ kết thu được, tác giả đưa kết luận, vấn đề tồn chưa giải đề xuất hướng phát triển đề tài SUMMARY SUMMARY Defects, especially cracks affect strongly on working and life of structures Among recent research works which deal with this area, the affection analysis of the cracks based on the change of the dynamic properties has brought useful results in estimating the working condition of structures Nevertheless, the using of dynamic properties has still exits some restrictions On the other side, the thesis concentrates on studying the bending plate with a crack which involves in the transverse displacements by the wavelet transform technique The calculation to find the transverse displacements is proceeded by the finite element method based on both the 8-node quadrilateral isoparametric element and the isoparametric element of Barsoum These results are examined by the Ansys software Thereafter, the transverse displacements transformed by the wavelet are used to analyse their sensitive to the length, the location of the crack To support to the above studying, some related problems which are studied and solved are presented in chapters: Chapter 1: General introduction Presenting primarily about defects, particularly cracks in structures The developing of the estimating problem, the history of the bending plate, the history of the finite element method and the history of wavelet Thereupon, the scientific and practical import of the thesis are to set up Chapter 2: The theory of the fracture mechanics, the bending plate and the finite element method Presenting primarily about the basic knowledge of the fracture mechanics such as elastic properties, the origination of the plastic region at the edged crack then selecting the element model describes the working of the singular region The basic theories of the bending plate are followed the Kirchhoff’s classical theory and the ReissnerMindlin’s theory, the establishment of the stiffness matrix and the load vector for the bending plate by the 8-node quadrilateral isoparametric element model obtained from the Reissner-Mindlin’s theory Chapter 3: The wavelet transform method Presenting primarily about wavelet knowledge, the basis of wavelet analysis, the process of wavelet transforms and wavelet reconstruction used in 1-dimension signal analysis, the application of the wavelet toolbox with Matlab for analysing signals Chapter 4: The determination to the transverse displacements of a bending plate using the finite element method SUMMARY Presenting some concrete problems which are applied to the established finite element method deals with the transverse displacements of the bending plate The results are compared with the one that obtained from the Ansys software plate Chapter 5: Wavelet transforms for the transverse displacements of the bending Base on the numerical results of the transverse displacements obtained from chapter 4, using the wavelet transform to investigate sensitivity of the transverse displacements in some specific cases Chapter 6: Conclusions and perspectives From the obtained results, the advantages and limitations of the proposed approach as well as some problems that has not yet resolved are taken out and the further developments of the thesis are presented and discussed MỤC LỤC MỤC LỤC LỜI CẢM ƠN TÓM TẮT LUẬN VĂN Chương 1: TỔNG QUAN 11 1.1 Ý nghĩa khoa học thực tiễn luận văn - 11 1.2 Lịch sử phát triển toán chẩn đoán kỹ thuật 13 1.3 Lịch sử phát triển toán chịu uốn - 14 1.4 Phương pháp phần tử hữu hạn (PTHH) - 15 1.5 Lịch sử phát triển wavelet 16 1.6 Nhiệm vụ luận văn - 16 Chương 2: LÝ THUYẾT CƠ HỌC PHÁ HUỶ, LÝ THUYẾT TẤM CHỊU UỐN VÀ PHƯƠNG PHÁP PTHH 18 2.1 Lý thuyết học phá hủy 18 2.2 Lý thuyết chịu uốn 26 2.3 Phương pháp PTHH cho toán chịu uốn - 34 Chương 3: PHƯƠNG PHÁP BIẾN ĐỔI WAVELET 43 3.1 Cơ sở lý thuyết wavelet 43 3.2 Cơ sở phân tích wavelet - 46 3.3 Biến đổi wavelet cho tín hiệu chiều - 48 3.4 Tái tạo wavelet cho tín hiệu chiều - 53 3.5 Giới thiệu họ wavelet hộp công cụ phần mềm Matlab - 57 Chương 4: XÁC ĐỊNH ĐỘ VÕNG TẤM BẰNG PHƯƠNG PHÁP PTHH 61 4.1 Sơ đồ khối chương trình xác định độ võng 62 4.2 Khảo sát hội tụ cho tốn khơng nứt 63 4.3 So sánh kết tính tốn độ võng lớn cho không nứt - 66 4.4 So sánh kết tính tốn độ võng lớn cho có vết nứt - 68 4.5 Kết luận chương 75 Chương 5: BIẾN ĐỔI WAVELET CHO ĐỘ VÕNG TẤM 76 5.1 Sơ đồ khối chương trình biến đổi wavelet - 76 5.2 Khảo sát nhạy cảm độ võng 77 MỤC LỤC 5.3 Kết luận chương 102 Chương 6: KẾT LUẬN VÀ KIẾN NGHỊ 103 6.1 Kết luận 103 6.2 Kiến nghị hướng phát triển 104 TÀI LIỆU THAM KHẢO 105 PHỤ LỤC B axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; C=cwt(TRAN_N12,Scale,Wname); for j=1:length(C) Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on;C=cwt(TRAN_N16,Scale,Wname); for j=1:length(C) Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News BÀI TOÁN % -% BIEN DOI WAVELET LIEN TUC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletC_Lc_ynut10_2tt_ngam.m % (Chieu dai vet nut thay doi, vi tri vet nut ynut=10cm) % -clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm P=2.5; % Luc tap trung P=2.5(N) tac dung tai cac vi tri % (x=5cm,y=35cm) va (x=35cm,y=35cm) h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y ynut=10; % Tung vet nut y =10cm Trang 68 PHỤ LỤC B % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N4,TRAN_N8,TRAN_N12,TRAN_N16]=Lc_ynut10_2tt_ngam; Wname=input('Chon loai wavelet phan tich:'); Scale=input('Chon ty le phan tich:'); figure(1); hold on; Cu=cwt(TRAN_KN,Scale,Wname); for j=1:length(Cu) Yuw(1,j)=Cu(1,j)/max(max(Cu),-min(Cu)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 1]); axis([-0.1 1.1 -1.1 1.1]); grid on; title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; C=cwt(TRAN_N4,Scale,Wname); for j=1:length(C) Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(3); hold on; C=cwt(TRAN_N8,Scale,Wname); for j=1:length(C) Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=8cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News figure(4); hold on; C=cwt(TRAN_N12,Scale,Wname); for j=1:length(C) Trang 69 PHỤ LỤC B Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on;C=cwt(TRAN_N16,Scale,Wname); for j=1:length(C) Yw(1,j)=C(1,j)/max(max(C),-min(C)); end plot(AY,Yw','-','LineWidth',1,'Color',[0 0.5 1]); title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so wavelet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News Trang 70 PHỤ LỤC B B3 CHƯƠNG TRÌNH BIẾN ĐỔI WAVELET RỜI RẠC CHO ĐỘ VÕNG TẤM BÀI TOÁN % -% BIEN DOI WAVELET ROI RAC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletD_Lc4_ynut_2tt.m % (Chieu dai vet nut Lc=4cm, vi tri vet nut ynut thay doi) % -clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm Lc=4; % Chieu dai vet nut Lc=4cm P=2.5; % Luc tap trung P=2.5(N) tac dung tai cac vi tri % (x=5cm,y=35cm) va (x=35cm,y=35cm) h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N2,TRAN_N10,TRAN_N20,TRAN_N30,TRAN_N38]=Lc4_ynut_2tt; Wname=input('Chon loai wavelet phan tich:'); figure(1); hold on; [Cu,Lu]=wavedec(TRAN_KN,1,Wname); Du1=wrcoef('d',Cu,Lu,Wname,1); for j=1:length(Du1) Yuw(1,j)=Du1(1,j)/max(max(Du1),-min(Du1)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 0]); title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; [C,L]=wavedec(TRAN_N2,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=2cm','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; Trang 71 PHỤ LỤC B xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(3); hold on; [C,L]=wavedec(TRAN_N10,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; [C,L]=wavedec(TRAN_N20,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=20cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on; [C,L]=wavedec(TRAN_N30,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=30cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News News figure(6); hold on; [C,L]=wavedec(TRAN_N38,1,Wname); D1=wrcoef('d',C,L,Wname,1); Trang 72 PHỤ LỤC B for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=38cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News BÀI TOÁN % -% BIEN DOI WAVELET ROI RAC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletD_Lc4_ynut_2tt.m % (Chieu dai vet nut thay doi, vi tri vet nut ynut=10cm) % -clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm P=2.5; % Luc tap trung P=2.5(N) tac dung tai cac vi tri % (x=5cm,y=35cm) va (x=35cm,y=35cm) h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y ynut=10; % Tung vet nut y =10cm % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N4,TRAN_N8,TRAN_N12,TRAN_N16]=Lc_ynut10_2tt; Wname=input('Chon loai wavelet phan tich:'); figure(1); hold on; [Cu,Lu]=wavedec(TRAN_KN,1,Wname); Du1=wrcoef('d',Cu,Lu,Wname,1); for j=1:length(Du1) Yuw(1,j)=Du1(1,j)/max(max(Du1),-min(Du1)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 0]); title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; Trang 73 PHỤ LỤC B [C,L]=wavedec(TRAN_N4,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(3); hold on; [C,L]=wavedec(TRAN_N8,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=8cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; [C,L]=wavedec(TRAN_N12,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News News figure(5); hold on; [C,L]=wavedec(TRAN_N16,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); Trang 74 PHỤ LỤC B title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News BÀI TOÁN % -% BIEN DOI WAVELET ROI RAC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletD_Lc_ynut10_pb.m % (Chieu dai vet nut Lc=4cm, vi tri vet nut ynut thay doi) % -clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm q=0.001; % Luc phan bo deu voi cuong q=0.001N/cm2 h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y ynut=10; % Tung vet nut y =10cm % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N4,TRAN_N8,TRAN_N12,TRAN_N16]=Lc_ynut10_pb; Wname=input('Chon loai wavelet phan tich:'); figure(1); hold on; [Cu,Lu]=wavedec(TRAN_KN,1,Wname); Du1=wrcoef('d',Cu,Lu,Wname,1); for j=1:length(Du1) Yuw(1,j)=Du1(1,j)/max(max(Du1),-min(Du1)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 0]); title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; [C,L]=wavedec(TRAN_N4,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end Trang 75 PHỤ LỤC B plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(3); hold on; [C,L]=wavedec(TRAN_N8,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=8cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; [C,L]=wavedec(TRAN_N12,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on; [C,L]=wavedec(TRAN_N16,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News News News Trang 76 PHỤ LỤC B BÀI TOÁN % -% BIEN DOI WAVELET ROI RAC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletD_Lc4_ynut_1tt.m % (Chieu dai vet nut thay doi, vi tri vet nut ynut=10cm) % -clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm P=5; % Luc tap trung P=5N tac dung tai vi tri x=20cm,y=35cm h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y ynut=10; % Tung vet nut y =10cm % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N4,TRAN_N8,TRAN_N12,TRAN_N16]=Lc_ynut10_1tt; Wname=input('Chon loai wavelet phan tich:'); figure(1); hold on; [Cu,Lu]=wavedec(TRAN_KN,1,Wname); Du1=wrcoef('d',Cu,Lu,Wname,1); for j=1:length(Du1) Yuw(1,j)=Du1(1,j)/max(max(Du1),-min(Du1)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 0]); title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; [C,L]=wavedec(TRAN_N4,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News Trang 77 PHỤ LỤC B figure(3); hold on; [C,L]=wavedec(TRAN_N8,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=8cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; [C,L]=wavedec(TRAN_N12,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on; [C,L]=wavedec(TRAN_N16,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News News BÀI TOÁN % -% BIEN DOI WAVELET ROI RAC % DOI VOI TRUONG HOP BAI TOAN % File: WaveletD_Lc_ynut10_2tt_ngam.m % (Chieu dai vet nut thay doi, vi tri vet nut ynut=10cm) % Trang 78 PHỤ LỤC B clear all,clc; format short; % DU LIEU CUA BAI TOAN -L=40; % Chieu dai tam L=40cm W=40; % Chieu rong tam W=40cm P=2.5; % Luc tap trung P=2.5(N) tac dung tai cac vi tri % (x=5cm,y=35cm) va (x=35cm,y=35cm) h=0.05; % Chieu day tam h=0.05cm E=7.1e6; % Modun dan hoi Young E=7.1e6N/cm2 nuy=0.3; % He so Poisson=0.3 nx=20; % So phan tu theo phuong x ny=20; % So phan tu theo phuong y ynut=10; % Tung vet nut y =10cm % -BIEN DOI WAVELET LIEN TUC -AY=[0:(1/13)/W:1]; [TRAN_KN,TRAN_N4,TRAN_N8,TRAN_N12,TRAN_N16]=Lc_ynut10_2tt_ngam; Wname=input('Chon loai wavelet phan tich:'); figure(1); hold on; [Cu,Lu]=wavedec(TRAN_KN,1,Wname); Du1=wrcoef('d',Cu,Lu,Wname,1); for j=1:length(Du1) Yuw(1,j)=Du1(1,j)/max(max(Du1),-min(Du1)); end plot(AY,Yuw,'-','LineWidth',1,'Color',[0 0.5 0]); title('TAM KHONG NUT','Fontname','Times News Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(2); hold on; [C,L]=wavedec(TRAN_N4,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=4cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News figure(3); hold on; [C,L]=wavedec(TRAN_N8,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Trang 79 PHỤ LỤC B Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=8cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(4); hold on; [C,L]=wavedec(TRAN_N12,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=12cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) figure(5); hold on; [C,L]=wavedec(TRAN_N16,1,Wname); D1=wrcoef('d',C,L,Wname,1); for j=1:length(D1) Yw(1,j)=D1(1,j)/max(max(D1),-min(D1)); end plot(AY,Yw','-','LineWidth',1,'Color',[1 0]); title('TAM NUT CO Lc=16cm, ynut=10cm','Fontname','Times Roman','Color','black','Fontsize',12) axis([-0.1 1.1 -1.1 1.1]); grid on; xlabel('Vi tri chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) ylabel('He so chi tiet chuan hoa','Fontname','Times News Roman','Color','black','Fontsize',12) News News News Trang 80 PHỤ LỤC B B4 CÁC HÀM PHỤ % function [CVKN,CVN2,CVN10,CVN20,CVN30,CVN38]=Lc4_ynut_2tt % Mang chua vong cua tam can bien doi wavelet % function [CVKN,CVN4,CVN8,CVN12,CVN16]=Lc_ynut10_2tt % Mang chua vong cua tam can bien doi wavelet % function [CVKN,CVN4,CVN8,CVN12,CVN16]=Lc_ynut10_pb % Mang chua vong cua tam can bien doi wavelet % function [CVKN,CVN4,CVN8,CVN12,CVN16]=Lc_ynut10_1tt % Mang chua vong cua tam can bien doi wavelet % function [CVKN,CVN4,CVN8,CVN12,CVN16]=Lc_ynut10_2tt_ngam % Mang chua vong cua tam can bien doi wavelet % - Trang 81 LÝ LỊCH TRÍCH NGANG Họ tên: PHẠM NGỌC TIẾN Phái: Nam Ngày sinh: 21 – 08 – 1977 Nơi sinh: Bình Định Địa liên lạc: Khu tập thể giáo viên-Trường Cao đẳng Xây dựng số 3, 24-Nguyễn Du-Phường 7-Tp Tuy Hòa-Tỉnh Phú Yên Điện thoại liên lạc: 0906728521 QUÁ TRÌNH ĐÀO TẠO Từ 1996 đến 2001: Học đại học Trường Đại học Bách Khoa Tp Hồ Chí Minh Từ 2005 đến 2007: Học cao học Trường Đại học Bách Khoa Tp Hồ Chí Minh Q TRÌNH CƠNG TÁC Từ 2001 đến 2002: Làm việc cơng ty Xây dựng Nam Thành-Tp Hồ Chí Minh Từ 2003 đến nay: Công tác Trường Cao đẳng Xây dựng số 3, Tp Tuy Hòa, Phú Yên ... TÀI ĐỘ NHẠY CẢM CỦA ĐỘ VÕNG TẤM CÓ VẾT NỨT BẰNG PHƯƠNG PHÁP BIẾN ĐỔI WAVELET II NHIỆM VỤ VÀ NỘI DUNG Luận văn khảo sát tính nhạy cảm độ võng chịu uốn có vết nứt thông qua kỹ thuật biến đổi wavelet. .. tích nhạy cảm độ võng dầm có vết nứt kỹ thuật wavelet Kết độ võng dầm cần thiết cho việc biến đổi wavelet xác định phương pháp phần tử hữu hạn A Gentile A Messine [22], xác định mở rộng vết nứt. .. Y.Y Lee, H.F Lam, R.K.L Su [17], phân tích nhạy cảm độ võng có vết nứt phương pháp biến đổi wavelet Các số liệu độ võng cần thiết cho việc biến đổi wavelet xác định từ kết bán giải tích đo đạc