Luận án Xác định các tham số trong bài toán chẩn đoán kết cấu bằng phương pháp động để cải tiến công tác quản lý công trình cầu

BỘ GIÁO DỤC VÀ ĐÀO TẠO 
TRƢỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI 
---------------------- 
NGUYỄN TIẾN MINH 
XÁC ĐỊNH CÁC THAM SỐ TRONG BÀI TOÁN 
CHẨN ĐOÁN KẾT CẤU BẰNG PHƢƠNG PHÁP 
ĐỘNG ĐỂ CẢI TIẾN CÔNG TÁC QUẢN LÝ 
CÔNG TRÌNH CẦU 
LUẬN ÁN TIẾN SĨ KỸ THUẬT 
HÀ NỘI – 2017 
 BỘ GIÁO DỤC VÀ ĐÀO TẠO 
TRƢỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI 
---------------------- 
NGUYỄN TIẾN MINH 
XÁC ĐỊNH CÁC THAM SỐ TRONG BÀI TOÁN 
CHẨN ĐOÁN KẾT CẤU BẰNG PHƢƠNG PHÁP 
ĐỘNG ĐỂ CẢI TIẾN CÔNG TÁC QUẢN LÝ 
CÔNG TRÌNH CẦU 
Chuyên ngành :Kỹ thuật Xây dựng Cầu - Hầm 
Mã số : 62.58.25.05 
LUẬN ÁN TIẾN SĨ KỸ THUẬT 
NGƢỜI HƢỚNG DẪN KHOA HỌC: 
1. PGS.TS. Nguyễn Ngọc Long 
2. PGS.TS. Trần Đức Nhiệm 
HÀ NỘI – 2017 
 i 
LỜI CAM ĐOAN 
Tôi xin cam đoan đây là công trình nghiên cứu của riêng tôi. Các số liệu, kết 
quả nêu trong luận án là trung thực và chƣa đƣợc ai công bố trong bất kỳ công trình 
nào khác. 
 Hà Nội, ngày tháng năm 
Tác giả 
Nguyễn Tiến Minh 
 ii 
MỤC LỤC 
LỜI CAM ĐOAN ........................................................................................................... i 
MỤC LỤC ...................................................................................................................... ii 
DANH MỤC CÁC CHỮ VIẾT TẮT ............................................................................ v 
DANH MỤC CÁC BẢNG BIỂU ................................................................................ vii 
DANH MỤC CÁC HÌNH VẼ .................................................................................... viii 
MỞ ĐẦU ........................................................................................................................ 1 
CHƢƠNG 1TỔNG QUAN VỀ CHẨN ĐOÁN KẾT CẤU BẰNG PHƢƠNG PHÁP 
DAO ĐỘNG ................................................................................................................... 5 
1.1. Giới thiệu chung về bài toán chẩn đoán kỹ thuật công trình bằng phƣơng pháp 
dao động ..................................................................................................................... 5 
1.1.1. Khái niệm về chẩn đoán công trình và chẩn đoán kết cấu bằng phƣơng 
pháp dao động ......................................................................................................... 5 
1.1.2. Các phƣơng pháp giải bài toán chẩn đoán kỹ thuật công trình bằng phƣơng 
pháp dao động ......................................................................................................... 6 
1.2. Tình hình nghiên cứu trên thế giới về chẩn đoán kết cấu bằng phƣơng pháp dao 
động ............................................................................................................................ 7 
1.3. Tình hình nghiên cứu ở Việt Nam về chẩn đoán kết cấu bằng phƣơng pháp dao 
động .......................................................................................................................... 11 
1.4. Đo dao động trong điều kiện khai thác và tổng quan về lý thuyết nhận dạng dao 
động .......................................................................................................................... 14 
1.4.1. Đo dao động trong điều kiện khai thác ....................................................... 14 
1.4.2. Tổng quan về lý thuyết nhận dạng dao động kết cấu cầu ........................... 15 
1.4.3. Lƣới bố trí điểm đo trên KCN cầu .............................................................. 18 
1.4.4. Công nghệ cảm biến ................................................................................... 20 
Kết luận Chƣơng 1 ................................................................................................... 21 
CHƢƠNG 2 CƠ SỞ LÝ THUYẾT VỀ DAO ĐỘNG KẾT CẤU, PHƢƠNG PHÁP 
CHẨN ĐOÁN ĐỘNG VÀ CÁC THAM SỐ SỬ DỤNG TRONG BÀI TOÁN CHẨN 
ĐOÁN ĐỘNG .............................................................................................................. 23 
2.1. Cơ sở lý thuyết về dao động kết cấu. ................................................................ 23 
2.1.1. Phƣơng trình vi phân dao động của hệ một bậc tự do ................................ 23 
2.1.2. Dao động của hệ nhiều bậc tự do................................................................ 27 
2.1.3. Các tham số đặc trƣng dao động................................................................ 27 
2.2. Các phƣơng pháp chẩn đoán kết cấu dựa trên dao động ................................... 29 
2.2.1. Phƣơng pháp dựa trên sự thay đổi tần số .................................................... 29 
2.2.2. Phƣơng pháp dựa trên độ mềm biểu kiến ................................................... 30 
2.2.3. Phƣơng pháp dựa trên năng lƣợng biến dạng hình thức ............................. 31 
2.2.4. Phƣơng pháp dựa trên độ cong đàn hồi ...................................................... 31 
2.2.5. Phƣơng pháp dựa trên độ cong hình dạng mode ........................................ 32 
2.2.6. Phƣơng pháp dựa trên độ cong bề mặt do tải trọng rải đều ........................ 33 
 iii 
2.2.7. Phƣơng pháp dựa trên sự thay đổi độ cứng ................................................ 33 
2.2.8. Đánh giá và lựa chọn phƣơng pháp cho bài toán chẩn đoán KCN bằng dao 
động ...................................................................................................................... 35 
Kết luận chƣơng 2 .................................................................................................... 38 
CHƢƠNG 3 THỰC NGHIỆM ĐO DAO ĐỘNG MỘT SỐ KCN CẦU TRÊN ĐỊA 
BÀN TP HÀ NỘI VÀ XÂY DỰNG MÔ HÌNH CHẨN ĐOÁN ĐỘNG KẾT CẤU 
CẦU ............................................................................................................................. 39 
3.1. Lựa chọn một số KCN cầu điển hình trên địa bàn TP Hà Nội để đo dao động ....... 39 
3.1.1. Hiện trạng hệ thống cầu trên địa bàn TP Hà Nội........................................ 39 
3.1.2. Một số công trình cầu điển hình áp dụng phƣơng pháp đo dao động ........ 41 
3.2. Xác định các tham số đặc trƣng dao động của các KCN cầu đƣợc lựa chọn ở 
thời điểm ban đầu. .................................................................................................... 42 
3.2.1. Xây dựng mô hình PTHH KCN cầu ........................................................... 42 
3.2.2. Tính toán các tham số đặc trƣng dao động của KCN cầu .......................... 44 
3.3. Thực nghiệm đo dao động 6 KCN cầu trên địa bàn TP Hà Nội. ...................... 49 
3.3.1. Trình tự thực hiện đo dao động .................................................................. 49 
3.3.2. Thiết bị đo ................................................................................................... 50 
3.3.3. Bố trí điểm đo dao động ............................................................................. 51 
3.3.4. Kết quả đo dao động KCN của 6 cầu ......................................................... 56 
3.3.5. So sánh kết quả tính và kết quả đo dao động 6 cầu .................................... 72 
3.4. Phân tích ảnh hƣởng của các hƣ hỏng đến đặc trƣng dao động của KCN cầu BT
 .................................................................................................................................. 74 
3.4.1. Miêu tả mô hình .......................................................................................... 74 
3.4.2. Bê tông bị suy giảm mô đun đàn hồi (E) .................................................... 76 
3.4.3. KCN có vết nứt hoặc hƣ hỏng .................................................................... 79 
3.5. Phân tích các đặc điểm của kết cấu nhạy cảm với sự thay đổi đặc trƣng dao 
động .......................................................................................................................... 84 
3.5.1. Thuộc tính cản ............................................................................................ 84 
3.5.2. Thuộc tính độ cứng ..................................................................................... 85 
3.5.3. Thuộc tính khối lƣợng ................................................................................ 85 
3.5.4. Ảnh hƣởng do khối lƣợng của phần kết cấu phụ ........................................ 86 
3.5.5. Phân tích ảnh hƣởng của độ cứng gối cầu đến đặc trƣng dao động ........... 86 
3.6. Xác định vị trí hƣ hỏng trên KCN cầu dầm bằng phƣơng pháp chẩn đoán động
 .................................................................................................................................. 88 
3.7. So sánh chi phí thử tải theo phƣơng pháp tĩnh và phƣơng pháp động .............. 92 
Kết luận Chƣơng 3 ................................................................................................... 93 
CHƢƠNG 4 ĐỀ XUẤT ỨNG DỤNG PHƢƠNG PHÁP ĐO DAO ĐỘNG VÀ CHẨN 
ĐOÁN ĐỘNG VÀO CÔNG TÁC QUẢN LÝ CẦU CỦA TP HÀ NỘI .................... 95 
4.1. Đề xuất tích hợp bổ sung một số đặc trƣng dao động vào hệ thống các tham số 
cần theo dõi, đo đạc và kiểm tra trong công tác quản lý khai thác cầu TP Hà Nội. 95 
 iv 
4.2. Xây dựng chƣơng trình máy tính dùng để chẩn đoán KCN cầu bằng phƣơng 
pháp dao động dựa trên các tham số đã lựa chọn. .................................................... 95 
4.3. Xây dựng quy trình quản lý và chẩn đoán KCN bằng phƣơng pháp dao động ....... 99 
4.3.1. Thu thập hồ sơ ............................................................................................ 99 
4.3.2. Khảo sát hiện trạng kết cấu, đánh giá dựa trên quan sát bằng mắt thƣờng 99 
4.3.3. Lập mô hình kết cấu ban đầu ...................................................................... 99 
4.3.4. Bố trí điểm đo dao động trên KCN cầu .................................................... 100 
4.3.5. Đo dao động KCN cầu .............................................................................. 101 
4.3.6. Hiệu chỉnh mô hình KCN cầu .................................................................. 102 
4.3.7. Xác định và theo dõi các tham số đặc trƣng dao động của KCN. ............ 103 
4.3.8. Quản lý thiết bị máy móc, quản lý kỹ thuật máy ...................................... 103 
4.3.9. Quản lý số liệu đo và tổ chức thực hiện ................................................... 104 
4.3.10. Đề xuất xây dựng ngƣỡng cảnh báo ....................................................... 104 
Kết luận chƣơng 4 .................................................................................................. 104 
KẾT LUẬN VÀ KIẾN NGHỊ ................................................................................... 106 
TÀI LIỆU THAM KHẢO .......................................................................................... 108 
DANH MỤC CÁC CÔNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN ........................ 1 
PHỤ LỤC ....................................................................................................................... 2 
 v 
DANH MỤC CÁC CHỮ VIẾT TẮT 
AASHTO American Association of State Highway and Transportation Officials 
(Hiệp hội Giao thông và Vận tải đƣờng bộ Hoa Kỳ) 
ASCE American Society of Civil Engineers (Hiệp hội Kỹ sƣ xây dựng Hoa 
Kỳ) 
BT Bê tông 
BTCT Bê tông cốt thép 
BTDUL Bê tông dự ứng lực 
CMIF Các chức năng chỉ định mode phức hợp 
CSI Phƣơng pháp nhận dạng không gian con - ngẫu nhiên 
CSI/ref Phƣơng pháp nhận dạng không gian con dựa trên tham chiếu 
CWT Biến đổi wavelet liên tục 
DAQ Hệ thống thu thập dữ liệu 
DOF Bậc tự do 
DUL Dự ứng lực 
DWT Biến đổi wavelet rời rạc 
EMA (FVT) Đo dao động cƣỡng bức 
EMD Phƣơng pháp phân tích dạng dao động cơ bản 
FEA Finite Element Analysis (phân tích phần tử hữu hạn) 
FEM Finite Element Method (phƣơng pháp phần tử hữu hạn) 
GTVT Giao thông Vận tải 
KC Kết cấu 
KCN Kết cấu nhịp 
LRFD Load Resistance Factor Design (thiết kế theo hệ số tải trọng và sức 
kháng) 
MFC Sự thay đổi độ mềm biểu kiến kết cấu 
OMA (AVT) Đo dao động trong điều kiện khai thác 
OMAX Đo dao động tổng hợp 
PP Phƣơng pháp chọn đỉnh 
PSD Mật độ phổ công suất 
PTHH Phần tử hữu hạn 
 vi 
SHM Theo dõi tình trạng sức khỏe kết cấu 
SSI Phƣơng pháp nhận dạng không gian con ngẫu nhiên 
SWT Biến đổi wavelet dừng 
TP Thành phố 
WPT Biến đổi wavelet packet 
RFV Residual Force Vector Method - phƣơng pháp véc tơ lực dƣ 
 vii 
DANH MỤC CÁC BẢNG BIỂU 
TT Tên bảng Trang 
Bảng 3.1. Thống kê số liệu công trình cầu trên địa bàn TP Hà Nội ......................... 39 
Bảng 3.2. Mƣời mode dao động đầu tiên của mô hình PTHH.................................. 45 
Bảng 3.3. Sáu mode dao động tính toán đầu tiên của cầu Kiến Hƣng ...................... 47 
Bảng 3.4. Sáu mode dao động tính toán đầu tiên của cầu Cống Thần ...................... 47 
Bảng 3.5. Khai báo tọa độ các nút dầm công son ..................................................... 62 
Bảng 3.6. Các tần số tìm đƣợc từ dữ liệu đo dao động dầm công son ...................... 66 
Bảng 3.7. Kết quả xử lý số liệu dao động KCN cầu Kiến Hƣng .............................. 69 
Bảng 3.8. Kết quả xử lý số liệu dao động KCN cầu Phùng Xá ................................ 69 
Bảng 3.9. Kết quả xử lý số liệu dao động KCN cầu Cống Thần .............................. 70 
Bảng 3.10. Kết quả xử lý số liệu dao động KCN cầu Tế Tiêu ................................. 70 
Bảng 3.11. Kết quả xử lý số liệu dao động KCN cầu La Khê .................................. 71 
Bảng 3.12. Kết quả xử lý số liệu đo KCN cầu Giẽ ................................................... 71 
Bảng 3.13. So sánh tần số tính và đo của cầu Kiến Hƣng ........................................ 73 
Bảng 3.14. So sánh tần số tính và đo cầu Cống Thần ............................................... 73 
Bảng 3.15. So sánh tần số tính và đo của cầu La Khê .............................................. 74 
Bảng 3.16. So sánh tần số tính và đo của cầu Giẽ .................................................... 74 
Bảng 3.17. Sự thay đổi tần số khi giảm mô đun đàn hồi của BT bản mặt cầu. .............. 76 
Bảng 3.18. Sự thay đổi tần số khi giảm mô đun đàn hồi của BT dầm. ..................... 77 
Bảng 3.19. Sự thay đổi tần số khi giảm mô đun đàn hồi của BT bản mặt cầu 
và dầm ...................................................................................................... 77 
Bảng 3.20. Sự thay đổi tần số khi giảm mô đun đàn hồi của BT dầm. ..................... 82 
Bảng 3.21. Tổng hợp chi phí kiểm định trung bình cho 1 KCN cầu dầm ................ 92 
 viii 
DANH MỤC CÁC HÌNH VẼ 
TT Tên hình Trang 
Hình 1.1. Minh họa biểu đồ ổn định dải tần số đƣợc xử lý từ dữ liệu đo dao 
động .................................................................................................................. 16 
Hình 1.2. Minh họa vị trí lắp các cảm biến tham chiếu (khoanh tròn) trên 
cầu. .................................................................................................................... 19 
Hình 1.3. B ... C 3.2 để xử lý số liệu đo và thực hiện nhận dạng dao động 
của kết cấu. Kết quả đƣợc thể hiện ở các bảng và hình dƣới đây. 
 Bảng A12. Kết quả đo dao động KCN N1 cầu Giẽ 
Mode 
No. 
Tần số 
f (Hz) 
Chu 
kỳ 
T(s) 
Hệ số 
cản ξ 
(%) 
MPC 
(-) 
MP 
(o) 
MPD 
(o) 
Dạng dao động 
1 4.12 0.23 0.73 0.98 0.23 5.70 Dạng uốn thứ 1 
2 5.39 0.20 0.65 0.99 2.66 5.10 Dạng uốn thứ 2 
3 13.56 0.07 0.64 0.88 4.20 6.50 Dạng xoắn thứ 1 
4 14.92 0.06 0.72 0.92 3.21 5.40 Dạng uốn thứ 3 
5 15.26 0.06 0.45 0.97 0.90 3.50 Dạng uốn thứ 4 
6 16.33 0.06 0.58 0.88 4.67 6.40 Dạng xoắn thứ 2 
Bảng A13. Kết quả đo dao động KCN N2 cầu Giẽ 
Mode 
No. 
Tần số 
f (Hz) 
Chu 
kỳ 
T(s) 
Hệ số 
cản ξ 
(%) 
MPC 
(-) 
MP 
(o) 
MPD 
(o) 
Dạng dao động 
1 4.45 0.22 0.73 0.95 0.25 4.50 Dạng uốn thứ 1 
2 4.92 0.20 0.65 0.97 2.47 5.50 Dạng uốn thứ 2 
3 14.13 0.07 0.64 0.88 4.70 3.60 Dạng xoắn thứ 1 
4 15.67 0.06 0.72 0.91 3.50 4.70 Dạng uốn thứ 3 
5 16.93 0.06 0.45 0.89 0.97 2.10 Dạng uốn thứ 4 
6 17.20 0.06 0.58 0.88 3.27 4.80 Dạng xoắn thứ 2 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A35: Mode dao động thứ 1: f= 4.34 (Hz) 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A36: Mode dao động thứ 2: f= 4.93 (Hz) 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A37: Mode dao động thứ 3: f= 14.26 (Hz) 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A38: Mode dao động thứ 4: f= 15.69 (Hz) 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A39: Mode dao động thứ 5: f= 16.78 (Hz) 
Mode dao động thực đo Kết quả mô hình phần mềm 
Hình A40: Mode dao động thứ 6: f= 17.19 (Hz) 
8. Code của chƣơng trình “Chẩn đoán KCN cầu – MFC Version 1.0” 
function varargout = MFC_v1(varargin) 
% MFC_V1 MATLAB code for MFC_v1.fig 
% MFC_V1, by itself, creates a new MFC_V1 or raises the 
existing 
% singleton*. 
% 
% H = MFC_V1 returns the handle to a new MFC_V1 or the 
handle to 
% the existing singleton*. 
% 
% MFC_V1('CALLBACK',hObject,eventData,handles,...) calls 
the local 
% function named CALLBACK in MFC_V1.M with the given 
input arguments. 
% 
% MFC_V1('Property','Value',...) creates a new MFC_V1 or 
raises the 
% existing singleton*. Starting from the left, property 
value pairs are 
% applied to the GUI before MFC_v1_OpeningFcn gets 
called. An 
% unrecognized property name or invalid value makes 
property application 
% stop. All inputs are passed to MFC_v1_OpeningFcn via 
varargin. 
% 
% *See GUI Options on GUIDE's Tools menu. Choose "GUI 
allows only one 
% instance to run (singleton)". 
% 
 % See also: GUIDE, GUIDATA, GUIHANDLES 
% Edit the above text to modify the response to help MFC_v1 
% Last Modified by GUIDE v2.5 08-Mar-2017 12:58:56 
% Begin initialization code - DO NOT EDIT 
gui_Singleton = 1; 
gui_State = struct('gui_Name', mfilename, ... 
 'gui_Singleton', gui_Singleton, ... 
 'gui_OpeningFcn', @MFC_v1_OpeningFcn, ... 
 'gui_OutputFcn', @MFC_v1_OutputFcn, ... 
 'gui_LayoutFcn', [] , ... 
 'gui_Callback', []); 
if nargin && ischar(varargin{1}) 
 gui_State.gui_Callback = str2func(varargin{1}); 
end 
if nargout 
 [varargout{1:nargout}] = gui_mainfcn(gui_State, 
varargin{:}); 
else 
 gui_mainfcn(gui_State, varargin{:}); 
end 
% End initialization code - DO NOT EDIT 
% --- Executes just before MFC_v1 is made visible. 
function MFC_v1_OpeningFcn(hObject, eventdata, handles, 
varargin) 
% This function has no output args, see OutputFcn. 
% hObject handle to figure 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% varargin command line arguments to MFC_v1 (see VARARGIN) 
% Choose default command line output for MFC_v1 
handles.output = hObject; 
 % Update handles structure 
guidata(hObject, handles); 
% This sets up the initial plot - only do when we are 
invisible 
% so window can get raised using MFC_v1. 
if strcmp(get(hObject,'Visible'),'off') 
 %plot(rand(5)); 
 plot(membrane); 
%axes(handles.axes1); 
%view([-55 50]); 
%cla; 
end 
% UIWAIT makes MFC_v1 wait for user response (see UIRESUME) 
% uiwait(handles.figure1); 
% --- Outputs from this function are returned to the command 
line. 
function varargout = MFC_v1_OutputFcn(hObject, eventdata, 
handles) 
% varargout cell array for returning output args (see 
VARARGOUT); 
% hObject handle to figure 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% Get default command line output from handles structure 
varargout{1} = handles.output; 
% -----------------------------------------------------------
--------- 
function FileMenu_Callback(hObject, eventdata, handles) 
% hObject handle to FileMenu (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
 % handles structure with handles and user data (see 
GUIDATA) 
% -----------------------------------------------------------
--------- 
function OpenMenuItem_Callback(hObject, eventdata, handles) 
% hObject handle to OpenMenuItem (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
file = uigetfile('*.fig'); 
if ~isequal(file, 0) 
 open(file); 
end 
% -----------------------------------------------------------
--------- 
function PrintMenuItem_Callback(hObject, eventdata, handles) 
% hObject handle to PrintMenuItem (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
printdlg(handles.figure1) 
% -----------------------------------------------------------
--------- 
function CloseMenuItem_Callback(hObject, eventdata, handles) 
% hObject handle to CloseMenuItem (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
selection = questdlg(['Close ' get(handles.figure1,'Name') 
'?'],... 
 ['Close ' get(handles.figure1,'Name') 
'...'],... 
 'Yes','No','Yes'); 
 if strcmp(selection,'No') 
 return; 
end 
delete(handles.figure1) 
% --- Executes on selection change in popupmenu1. 
function popupmenu1_Callback(hObject, eventdata, handles) 
% hObject handle to popupmenu1 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% Hints: contents = get(hObject,'String') returns popupmenu1 
contents as cell array 
% contents{get(hObject,'Value')} returns selected item 
from popupmenu1 
% --- Executes during object creation, after setting all 
properties. 
function popupmenu1_CreateFcn(hObject, eventdata, handles) 
% hObject handle to popupmenu1 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles empty - handles not created until after all 
CreateFcns called 
% Hint: popupmenu controls usually have a white background on 
Windows. 
% See ISPC and COMPUTER. 
if ispc && isequal(get(hObject,'BackgroundColor'), 
get(0,'defaultUicontrolBackgroundColor')) 
 set(hObject,'BackgroundColor','white'); 
end 
 set(hObject, 'String', {'plot(rand(5))', 
'plot(sin(1:0.01:25))', 'bar(1:.5:10)', 'plot(membrane)', 
'surf(peaks)'}); 
% --- Executes on button press in pushbutton2. 
function pushbutton2_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton2 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
[filename1,filepath1]=uigetfile({'*.mat','MATLAB data'},'Chon 
file'); 
% --- Executes on button press in pushbutton3. 
function pushbutton3_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton3 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
[filename1,filepath1]=uigetfile({'*.*','All Files'},'Chon 
file'); 
global fh 
fh=load([filepath1 filename1]); 
% --- Executes on button press in pushbutton5. 
function pushbutton5_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton5 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
[filename1,filepath1]=uigetfile({'*.*','All Files'},'Chon 
file'); 
global phih 
phih=load([filepath1 filename1]); 
% --- Executes on button press in pushbutton6. 
 function pushbutton6_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton6 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
global fh 
global fd 
global phih 
global phid 
omega=2*pi*fh; 
lambda=1./omega./omega; 
global n_nodes 
global n_modes 
n_nodes=length(phih(:,1)); 
n_modes=length(phih(1,:)); 
L=[]; 
for k=1:n_modes 
 L(k,k)=lambda(k); 
end 
Fh=phih*L*phih' 
omegad=2*pi*fd; 
lambdad=1./omegad./omegad; 
Ld=[]; 
for k=1:n_modes 
 Ld(k,k)=lambdad(k); 
end 
Fd=phid*Ld*phid' 
dF=Fd-Fh 
X=ones(n_nodes,1) 
global MFC1 
global MFC 
global x y 
MFC1=dF*X 
MFC=[]; 
for k=1:n_nodes 
 MFC=[MFC dF(k,k)] 
end 
axes(handles.axes1); 
 cla; 
% Normalize MFC for La Khe 
Ld=0; 
ngirder=0; 
spacing=0; 
St = get(handles.edit1,'String'); 
Ld=str2num(St); 
ngirder = str2num(get(handles.edit2,'String')); 
spacing = str2num(get(handles.edit3,'String')); 
%spacing = 2.4; 
%ngirder = 6; 
ndiv_ = 24; 
x=0:Ld/ndiv_:Ld; 
y=0:1:ngirder-1; 
y=y*spacing; 
MFC=reshape(MFC,ndiv_+1,ngirder); 
MFC=MFC'; 
global nMFC 
nMFC=MFC; 
for k=1:ngirder 
 mx=max(abs(nMFC(k,:))); 
 nMFC(k,:)=nMFC(k,:)./mx; 
end 
global MFCI 
MFCI=[]; 
MFCI=nMFC.*abs(nMFC); 
surf(x,y,nMFC) 
xlabel('Doc cau [m]') 
ylabel('Ngang cau [m]') 
zlabel('Chi so MFC') 
view([-55 50]); 
rotate3d 
%plot(membrane) 
% --- Executes on button press in pushbutton7. 
function pushbutton7_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton7 (see GCBO) 
 % eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
axes(handles.axes1); 
cla; 
grid on 
hold on 
global MFC 
global nMFC x y 
for k=1:6 
 plot3(x,y(k)+zeros(1,25),nMFC(k,:),'LineWidth',2) 
end 
xlabel('Doc cau [m]') 
ylabel('Ngang cau [m]') 
zlabel('Chi so MFC') 
view([-55 50]); 
rotate3d 
hold off 
% --- Executes on button press in pushbutton8. 
function pushbutton8_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton8 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
[filename1,filepath1]=uigetfile({'*.*','All Files'},'Chon 
file'); 
global fd 
fd=load([filepath1 filename1]); 
% --- Executes on button press in pushbutton9. 
function pushbutton9_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton9 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
 [filename1,filepath1]=uigetfile({'*.*','All Files'},'Chon 
file'); 
global phid 
phid=load([filepath1 filename1]); 
function edit1_Callback(hObject, eventdata, handles) 
% hObject handle to edit1 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% Hints: get(hObject,'String') returns contents of edit1 as 
text 
% str2double(get(hObject,'String')) returns contents 
of edit1 as a double 
% --- Executes during object creation, after setting all 
properties. 
function edit1_CreateFcn(hObject, eventdata, handles) 
% hObject handle to edit1 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles empty - handles not created until after all 
CreateFcns called 
% Hint: edit controls usually have a white background on 
Windows. 
% See ISPC and COMPUTER. 
if ispc && isequal(get(hObject,'BackgroundColor'), 
get(0,'defaultUicontrolBackgroundColor')) 
 set(hObject,'BackgroundColor','white'); 
end 
function edit2_Callback(hObject, eventdata, handles) 
 % hObject handle to edit2 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% Hints: get(hObject,'String') returns contents of edit2 as 
text 
% str2double(get(hObject,'String')) returns contents 
of edit2 as a double 
% --- Executes during object creation, after setting all 
properties. 
function edit2_CreateFcn(hObject, eventdata, handles) 
% hObject handle to edit2 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles empty - handles not created until after all 
CreateFcns called 
% Hint: edit controls usually have a white background on 
Windows. 
% See ISPC and COMPUTER. 
if ispc && isequal(get(hObject,'BackgroundColor'), 
get(0,'defaultUicontrolBackgroundColor')) 
 set(hObject,'BackgroundColor','white'); 
end 
function edit3_Callback(hObject, eventdata, handles) 
% hObject handle to edit3 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
% Hints: get(hObject,'String') returns contents of edit3 as 
text 
 % str2double(get(hObject,'String')) returns contents 
of edit3 as a double 
% --- Executes during object creation, after setting all 
properties. 
function edit3_CreateFcn(hObject, eventdata, handles) 
% hObject handle to edit3 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles empty - handles not created until after all 
CreateFcns called 
% Hint: edit controls usually have a white background on 
Windows. 
% See ISPC and COMPUTER. 
if ispc && isequal(get(hObject,'BackgroundColor'), 
get(0,'defaultUicontrolBackgroundColor')) 
 set(hObject,'BackgroundColor','white'); 
end 
% --- Executes on button press in pushbutton10. 
function pushbutton10_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton10 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA) 
h=msgbox({'Chuong trinh chan doan KCN cau MFC - Version 1.0' 
'Tac gia:' ' ThS. Nguyen Tien Minh' ' TS. Do Anh Tu'},... 
 'Gioi thieu'); 
% --- Executes on button press in pushbutton11. 
function pushbutton11_Callback(hObject, eventdata, handles) 
% hObject handle to pushbutton11 (see GCBO) 
% eventdata reserved - to be defined in a future version of 
MATLAB 
% handles structure with handles and user data (see 
GUIDATA)