Luận văn thạc sĩ VNU UET tính toán tham số TMD theo phương pháp biến đổi laplace

ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VIỆN CƠ HỌC PHAN THỊ TRÀ MY TÍNH TỐN THAM SỐ TMD THEO PHƯƠNG PHÁP BIẾN ĐỔI LAPLACE LUẬN VĂN THẠC SĨ Hà Nội - 2011 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CƠNG NGHỆ VIỆN CƠ HỌC PHAN THỊ TRÀ MY TÍNH TOÁN THAM SỐ TMD THEO PHƯƠNG PHÁP BIẾN ĐỔI LAPLACE Ngành: Cơ học Chuyên ngành: Cơ học vật thể rắn Mã số: 60 44 21 LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC: GS TSKH Nguyễn Đông Anh Hà Nội - 2011 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VIỆN CƠ HỌC PHAN THỊ TRÀ MY TÍNH TỐN THAM SỐ TMD THEO PHƯƠNG PHÁP BIẾN ĐỔI LAPLACE LUẬN VĂN THẠC SĨ Hà Nội - 2011 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VIỆN CƠ HỌC PHAN THỊ TRÀ MY TÍNH TỐN THAM SỐ TMD THEO PHƯƠNG PHÁP BIẾN ĐỔI LAPLACE Ngành: Cơ học Chuyên ngành: Cơ học vật thể rắn Mã số: 60 44 21 LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC: GS TSKH Nguyễn Đông Anh Hà Nội - 2011 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỤC LỤC Mở đầu Chương 1: Giới thiệu thiết bị tiêu tán lượng khối lượng dạng TMD 1.1 Giảm dao động thiết bị tiêu tán lượng 1.2 Ứng dụng thiết bị TMD 1.2.1 TMD dạng chuyển động theo hướng trực giao 1.2.2 TMD dạng lắc đơn 1.2.3 TMD dạng lắc nhiều tầng 1.2.4 TMD dạng lắc cho kết cấu mảnh 1.2.5 TMD giảm dao động cầu 1.2.6 Một số nghiên cứu thiết bị TTNL TMD Việt Nam 1.3 Giới thiệu phương pháp Laplace 11 1.4 Kết luận chương 14 Chương 2: Tính toán tham số tối ưu TMD theo phương pháp biến đổi Laplace 15 2.1 Chuyển động hệ chưa có TMD 15 2.2 Tính tốn tham số tối ưu TMD phương pháp biến đổi Laplace 17 2.2.1 Thiết lập hệ phương trình vi phân chuyển động hệ 17 2.2.2 Tính tốn tham số tối ưu TMD 19 2.2.3 Giải tìm nghiệm hệ phương trình 27 2.3 Kết luận chương 31 Chương 3: Mơ dao động hệ có gắn TMD 33 3.1 Mô dao động hệ 33 3.1.1 Chuyển động hệ chưa lắp TMD 34 3.1.2 Chuyển động hệ lắp thêm TMD 35 3.2 So sánh với số phương pháp khác 37 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 3.3 Kết luận chương 41 Kết luận 43 Tài liệu tham khảo 45 Phụ lục 46 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CÁC KÝ HIỆU VIẾT TẮT BĐL: HMD: MATLAB: TBTTNL: Biến đổi Laplace Hybrid mass damper Matrix laboratory Thiết bị tiêu tán lượng TTNL: Tiêu tán lượng TMD: Tuned mass damper LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỞ ĐẦU Trong trình làm việc tác động ngoại lực, máy móc thiết bị, kết cấu phát sinh dao động không mong muốn Những dao động ảnh hưởng xấu đến điều kiện làm việc làm giảm độ bền chi tiết phận máy, kết cấu Thông thường để giảm dao động có hại, người ta tăng độ cứng vững kết cấu khử nguyên nhân gây dao động Tuy nhiên người ta sử dụng phương pháp vì, phí tốn để tăng độ cứng vững kết cấu, khử hết nguyên nhân gây dao động Trong vài thập kỷ trở lại đây, giới phát triển công nghệ sử dụng thiết bị tiêu tán lượng (TBTTNL) để giảm dao động TBTTNL thiết bị lắp vào hệ kỹ thuật với mục đích hấp thụ chuyển phần lượng dao động hệ thành dạng lượng khác, từ làm giảm dao động va chạm có hại hệ kỹ thuật Phương pháp chừng mực định bổ xung cho hai phương án Đối tượng phạm vi nghiên cứu Nghiên cứu ứng dụng TBTTNL vào cơng trình địi hỏi có nhiều giai đoạn Trong phân thành hai giai đoạn chính: - Nghiên cứu dao động cơng trình - Nghiên cứu dao động TBTTNL Từ cơng trình, ta áp dụng dạng TBTTNL cách phù hợp Vì kiến trúc kết cấu đa dạng nên hình dáng cấu trúc thiết bị khơng phần phong phú Trong luận văn ta xét đến dao động hệ bậc tự có gắn thiết bị TTNL dạng TMD Phương pháp nghiên cứu Sử dụng phép biến đổi Laplace xây dựng phương trình chuyển động, từ tìm tham số tối ưu tỷ số cản tỷ số tần số tối ưu LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Sử dụng phần mềm Matlab tính tốn số cho hệ So sánh với phương pháp điểm cố định phương pháp cực tiểu moment bậc để đánh giá độ tin cậy phương pháp luận văn Mục đích luận văn Mục đích luận văn xây dựng phương pháp tính toán mới, xác định tham số đặc trưng hệ Chứng minh tính đắn phương pháp giải tìm nghiệm Bố cục luận văn Luận văn bao gồm phần mở đầu, chương phần kết luận - Chương 1: Giới thiệu thiết bị tiêu tán lượng khối lượng TMD phương pháp biến đổi Laplace - Chương 2: Tính tốn tham số tối ưu TMD theo phương pháp biến đổi Laplace - Chương 3: Mơ dao động hệ có gắn TMD LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG GIỚI THIỆU THIẾT BỊ TIÊU TÁN NĂNG LƯỢNG DẠNG KHỐI LƯỢNG TMD VÀ PHƯƠNG PHÁP BIẾN ĐỔI LAPLACE 1.1 Giảm dao động thiết bị tiêu tán lượng Dao động tượng phổ biến tự nhiên kỹ thuật Các phương tiện giao thông vận tải, nhà cao tầng, cầu giao thông… hệ dao động kỹ thuật Dao động có ích có hại tuỳ thuộc vào trường hợp cụ thể Đối với hệ kỹ thuật, nhiều dạng dao động có hại, khơng mong muốn xảy ra: phương tiện giao thơng chịu kích động mặt đường, tàu thuỷ cơng trình ngồi khơi chịu tác động sóng gió, tháp vơ tuyến, cao ốc chịu tác động gió động đất, cầu giao thông nhịp lớn chịu tác động phương tiện vận tải, cầu treo chịu tải trọng gió bão, thiết bị tuốc bin hoạt động với tốc độ cao… Các loại dao động có hại ngày nguy hiểm có ảnh hưởng xấu đến điều kiện làm việc, làm giảm độ bền máy móc, thiết bị, cơng trình Đồng thời, với phát triển khoa học kỹ thuật, tăng lên quy mô kết cấu, tốc độ máy móc cường độ kích động ngồi; cấp thiết việc giảm giá thành cơng trình lớn; u cầu cao an tồn cho cơng trình quan trọng mà cơng nghệ giảm dao động quan tâm hàng đầu nhiều quan nghiên cứu ứng dụng Trước đây, phương pháp phổ biến giảm dao động tăng độ cứng vững cơng trình khử ngun nhân gây dao động Tuy nhiên, nhiều trường hợp, chi phí để tăng độ cứng vững kết cấu tốn kết cấu cơng trình thiết kế, xây dựng theo quan điểm cổ điển dựa vào vững thân, tải tĩnh đặt kết cấu, khơng có khả thích ứng với thay đổi ln ln lực tác dụng từ ngồi mơi trường Vì vấn đề đặt phải thiết kế phận để lắp vào kết cấu, làm giảm dao động có hại cách tốt làm tăng độ an toàn, hiệu tuổi thọ cơng trình Vì vài thập kỷ trở lại đây, giới phát triển công nghệ sử dụng thiết bị tiêu tán lượng (TBTTNL) để giảm dao động Các phương pháp giảm dao động va chạm thường dựa sở làm giảm cách ly lực kích động, thay đổi tần số riêng để tránh cộng hưởng, tăng cản cho LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 56 if get(findobj(findobj('tag','mainmanu'),'tag','ch_act'),'value') if get(findobj(findobj('tag','mainmanu'),'tag','ch_opt'),'value') k2=k1*m2*m1/((m1+m2)^2); d2=m2*sqrt(k2/(m1+m2)); else k2=xx(5); d2=xx(6); end Umin=activectrl(xx); hw = waitbar(0,'Please wait '); [t3,y3]=ADAM(time_evl,[xo,dxo,x20,dx20]); close(hw); else t3=0;y3=0; end %================================================================= ========= %\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ function tai=taingaunhien(leng_sig,noisuy) x=rand(1,leng_sig)-0.5; tai=interp(x,noisuy); %///////////////////////////////////////////////////////////// function [t,y]=ADAM(t_in,yo) global step_nn tai step_nn=1; y(1,:)=yo; t(1)=t_in(1); ht=t_in(2); N=t_in(3)/ht; bacj=length(yo); % tich phan bang Runge-Kutta -for i=1:5 % tinh k1 for j=1:bacj k1(j)=f(j,t(i),y(i,:)); end % tinh k2 for j=1:bacj k2(j)=f(j,t(i)+ht/2,y(i,:)+k1*ht/2); end % tinh k3 for j=1:bacj k3(j)=f(j,t(i)+ht/2,y(i,:)+k2*ht/2); end LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 57 % tinh k4 for j=1:bacj k4(j)=f(j,t(i)+ht,y(i,:)+k3*ht); end % tinh y y(i+1,:)=y(i,:)+ht/6*(k1+2*k2+2*k3+k4); t(i+1)=t(i)+ht; step_nn=i; end % Tich phan Adam for i=5:N-1 for j=1:bacj y(i+1,j)=y(i,j)+ht*(55*f(j,t(i),y(i,:))-59*f(j,t(i-1),y(i-1,:))+37*f(j,t(i-2),y(i-2,:))-9*f(j,t(i3),y(i-3,:)))/24; end t(i+1)=t(i)+ht; step_nn=i; if mod(i,3)==0 waitbar(i/N); end end %///////////////////////////////////////////////////////////// %===============ham F 1DOF========================== function dydt = f(chiso,t,y) global m1 m2 d2 k2 k1 d1 p f g tai step_nn Lo_sw Umin switch Lo_sw case P=0; case P=p*sin(2*pi*f*t); case P=p*cos(2*pi*f*t); case P=p*tai(step_nn); otherwise msgbox('error!!!'); end % P=0; if length(y)==2 dydt_temp(1)=y(2); dydt_temp(2)=(P-d1*y(2)-k1*y(1))/m1; dydt=dydt_temp(chiso); else LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 58 U=Umin(1)*y(1)+Umin(2)*y(2)+Umin(3)*y(3)+Umin(4)*y(4); dydt_temp(1)=y(2); dydt_temp(2)=-k1/m1*y(1)-d1/m1*y(2)+k2/m1*y(3)+d2/m1*y(4)+P/m1-2*U/m1; dydt_temp(3)=y(4); dydt_temp(4)=k1/m1*y(1)+d1/m1*y(2)+(-k2/m1-1/m2*k2)*y(3)+(-d2/m11/m2*d2)*y(4)+1/m2*U-1/m1*P+2/m1*U; dydt=dydt_temp(chiso); end %============================================== %///////////////////////////////////////////////////////////// function fTMDn if get(findobj(findobj('tag','mainmanu'),'tag','cd_TMDn'),'value') set(findobj(findobj('tag','mainmanu'),'tag','cd_TMDn'),'value',1); set(findobj(findobj('tag','mainmanu'),'tag','var16'),'enable','on','BackgroundColor',[1 1]); set(findobj(findobj('tag','mainmanu'),'tag','var17'),'enable','on','BackgroundColor',[1 1]); set(findobj(findobj('tag','mainmanu'),'tag','var18'),'enable','on','BackgroundColor',[1 1]); % set(findobj(findobj('tag','mainmanu'),'tag','var9'),'enable','off','BackgroundColor',[1 0.8 0.8]); else set(findobj(findobj('tag','mainmanu'),'tag','var16'),'enable','off','BackgroundColor',[1 0.8 0.8]); set(findobj(findobj('tag','mainmanu'),'tag','var17'),'enable','off','BackgroundColor',[1 0.8 0.8]); set(findobj(findobj('tag','mainmanu'),'tag','var18'),'enable','off','BackgroundColor',[1 0.8 0.8]); % set(findobj(findobj('tag','mainmanu'),'tag','var9'),'enable','on','BackgroundColor',[1 1]); End function funct1(action) if(strcmp(action,'exit')) close all; elseif(strcmp(action,'reset')) for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var1'),'String',250000); % m1 set(findobj(findobj('tag','mainmanu'),'Tag','var2'),'String',15000); % m2 % set(findobj(findobj('tag','mainmanu'),'Tag','var16'),'String',3000,'enable','off','BackgroundColor ',[1 8]); % md % set(findobj(findobj('tag','mainmanu'),'Tag','var17'),'String',4,'enable','off','BackgroundColor',[1 8]); % md % set(findobj(findobj('tag','mainmanu'),'Tag','var18'),'String',38,'enable','off','BackgroundColor',[ 8]); % md set(findobj(findobj('tag','mainmanu'),'Tag','var3'),'String',1e+07); % k1 set(findobj(findobj('tag','mainmanu'),'Tag','var4'),'String',55000); % d1 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 59 set(findobj(findobj('tag','mainmanu'),'Tag','var5'),'String',533997.864,'enable','off','Background Color',[1 8]); % k2 =7e+005 set(findobj(findobj('tag','mainmanu'),'Tag','var6'),'String',21293.0602,'enable','off','Background Color',[1 8]); % d2=60000 set(findobj(findobj('tag','mainmanu'),'Tag','var7'),'String',5e+5); % p set(findobj(findobj('tag','mainmanu'),'Tag','var8'),'String',0.1); % om set(findobj(findobj('tag','mainmanu'),'Tag','var9'),'String',0.01); % xo set(findobj(findobj('tag','mainmanu'),'Tag','var10'),'String',0); % dxo set(findobj(findobj('tag','mainmanu'),'Tag','var11'),'String',30); % t set(findobj(findobj('tag','mainmanu'),'Tag','T_binh'),'enable','off','value',0); set(findobj(findobj('tag','mainmanu'),'Tag','num_r'),'enable','off','String',0); set(findobj(findobj('tag','mainmanu'),'Tag','num_st'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ch_1dof'),'value',1,'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','ch_2dof'),'value',1,'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','ch_opt'),'value',1,'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','ch_act'),'value',1,'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','ch4'),'value',0,'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','ch6'),'value',0); elseif(strcmp(action,'ngaunhien')) if get(findobj(findobj('tag','mainmanu'),'tag','taitrong'),'value')==4 set(findobj(findobj('tag','mainmanu'),'Tag','T_binh'),'enable','on','value',0); set(findobj(findobj('tag','mainmanu'),'Tag','num_r'),'enable','on','String',20); set(findobj(findobj('tag','mainmanu'),'Tag','num_st'),'enable','on'); else set(findobj(findobj('tag','mainmanu'),'Tag','T_binh'),'enable','off','value',0); set(findobj(findobj('tag','mainmanu'),'Tag','num_r'),'enable','off','String',0); set(findobj(findobj('tag','mainmanu'),'Tag','num_st'),'enable','off'); end elseif(strcmp(action,'return')) delete(findobj('tag','result')); elseif(strcmp(action,'khaosat')) set(findobj(findobj('tag','mainmanu'),'Tag','bool_ks'),'value',1); set(findobj(findobj('tag','mainmanu'),'Tag','ks_k2'),'enable','on'); funct1('ks_k2') elseif(strcmp(action,'codinh')) set(findobj(findobj('tag','mainmanu'),'Tag','bool_ks'),'value',0); set(findobj(findobj('tag','mainmanu'),'Tag','ks_k2'),'enable','off'); % set(findobj(findobj('tag','mainmanu'),'Tag','ks_d2'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_start_st'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step_st'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end_st'),'enable','off'); set(findobj(findobj('tag','mainmanu'),'Tag','bool_cd'),'value',1); LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 60 set(findobj(findobj('tag','mainmanu'),'Tag','cd_1dof'),'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','cd_2dof'),'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','cd_tktu'),'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','cd_sosanh'),'enable','on'); set(findobj(findobj('tag','mainmanu'),'Tag','cd_TMDn'),'enable','on'); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end if get(findobj(findobj('tag','mainmanu'),'tag','cd_tktu'),'value') set(findobj(findobj('tag','mainmanu'),'Tag','var8'),'enable','off','BackgroundColor',[1 0.8 0.8]); set(findobj(findobj('tag','mainmanu'),'Tag','var9'),'enable','off','BackgroundColor',[1 0.8 0.8]); end elseif(strcmp(action,'ks_k2')) ksat=get(findobj(findobj('tag','mainmanu'),'tag','ks_k2'),'value'); switch ksat case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',200000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',20000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',300000); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var1'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',150000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',20000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',250000); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var15'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',10000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',2000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',20000); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var2'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',2e+010); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',2e+09); LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 61 set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',4e+010); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var3'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',35000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',5000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',100000); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var4'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',35); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',2); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',55); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var5'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',4); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',1); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',10); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var6'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',35); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',2); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',55); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var7'),'enable','off','BackgroundColor',[1 0.8 0.8]); case set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',400000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',40000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',800000); for i=1:15 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 62 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var8'),'enable','off','BackgroundColor',[1 0.8 0.8]); case 10 set(findobj(findobj('tag','mainmanu'),'Tag','ks_start'),'String',10000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_step'),'String',2000); set(findobj(findobj('tag','mainmanu'),'Tag','ks_end'),'String',30000); for i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); end set(findobj(findobj('tag','mainmanu'),'Tag','var9'),'enable','off','BackgroundColor',[1 0.8 0.8]); end End function funct2 clc for i=1:11 texttg=['var' num2str(i)]; hand_edit=findobj(findobj('tag','mainmanu'),'tag',texttg); xx(i)= str2num(get(hand_edit,'String')); end [t1,y1,t2,y2,t3,y3]=equ(xx); size(t1) size(y1) size(t2) size(y2) inp.t1=t1; inp.y1=y1; inp.t2=t2; inp.y2=y2; inp.t3=t3; inp.y3=y3; inp.xx=xx; % aTMDe=(out1.y2(2:length(out1.y2),4)-out1.y2(1:length(out1.y2)-1,4))/(out1.t2(2)-out1.t2(1)); I1=out1.y2(1:length(out1.y2)-1,2).*aTMDe; figure;plot(out1.t1(1:length(out1.t1)-1),I1); set(gcf,'color',[1 1]);set(gca,'xlim',[0,40]); function khungnhin max_x=get(findobj(gcf,'tag','slid'),'value'); if get(findobj(findobj('tag','mainmanu'),'Tag','bool_ks'),'value') result1; else set(gca,'xlim',[0 max_x]); end LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 63 function Y = myfun(x) Y(:,1) = 200*sin(x(:))./x(:); Y(:,2) = x(:).^2; function out=neurap(kin,din) global kk1 dd1 kk1=kin;dd1=din; k=0; out=zeros(2,1); out1=zeros(2,1); out1(1)=kin;out1(2)=din; nn=40000;epsilon=0.005;n=1;dkdung=1; h=waitbar(0,'Dang tinh c2, k2 '); while dkdung>epsilon && n1 lege={'Đáp ứng không lắp TMD','Đờng bao đáp ứng không lắp TMD','Đáp ứng không lắp TMD', 'Đờng bao đáp ứng không lắp TMD','Đáp ứng có lắp TMD','Đờng bao đáp ứng có lắp TMDe', 'Đáp ứng có lắp TMDe','Đờng bao đáp ứng có lắp TMDe', 'Dịch chuyển TMDe','Đờng bao dịch chuyển TMDe', 'Đáp ứng có lắp Active','Đờng bao đáp ứng có lắp Active','Đáp ứng có lắp Active', 'Đờng bao đáp ứng có lắp Active','Dịch chuyển Active','Đờng bao dịch chuyển Active', 'Đáp ứng có lắp TMDn','Đờng bao đáp ứng có lắp TMDn','Đáp ứng có lắp TMDn', 'Đờng bao đáp ứng có lắp TMDn','Dịch chuyển TMDn','Đờng bao dịch chuyển TMDn', 'Đáp ứng cố định TMD','Đờng bao đáp ứng cố định TMD','Đáp ứng cố định TMD', 'Đờng bao đáp ứng cố định TMD'}; lg=legend(lege(lgd)); set(lg,'fontname','vk sans serif','fontsize',8); end max_x=get(findobj(gcf,'tag','slid'),'value'); set(gca,'xlim',[0 max_x]); function result farb1=[1 1]; f_g2=[0.05 0.06 0.9 0.85]; menu2=figure; set(menu2,'Units','normalized', 'Pos',f_g2,'color',farb1, 'NumberTitle','off', 'menubar','figure', 'menubar','figure', 'tag','result', 'Userdata',0,'Name','hien thi ket qua'); set(menu2,'DefaultUicontrolFontunits','points', 'DefaultUicontrolFontsize',9, LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 68 'DefaultUicontrolUnits','normalized', 'DefaultUicontrolHorizontalAlignment','left', 'DefaultUicontrolBackgroundColor',farb1, 'DefaultUicontrolMin',0, 'DefaultUicontrolMax',1, 'DefaultUicontrolStyle','text'); % lef=0.05;dx=0.05;dy=0.017; step=0.036;dis=0.05; j=0.9; % h1 = axes('Position',[lef+0.17 72 72],'fontname','vk sans serif','fontsize',8); uicontrol('pos',[lef+0.12,0.94,dx,2.5*dy],'string','(m)','fontsize',12); uicontrol('pos',[lef+0.87,0.2,dx,2.5*dy],'string','t (s)','fontsize',12); sli_max=str2num(get(findobj(findobj('tag','mainmanu'),'tag','var11'),'String')); uicontrol('style','slider','pos',[lef+0.17,.07,.72 05],'BackgroundColor',[1 1], 'Tag','slid','Min',1,'Max',sli_max,'SliderStep',[0.01 02],'value',sli_max,'callback','khungnhin') % % uicontrol('pos',[0.35,0.95,9*dx,2*dy],'string','Hieu qua giam dao dong cua TMDe:'); uicontrol('pos',[0.6,0.95,0.8*dx,2*dy],'string','','tag','p1'); uicontrol('pos',[0.68,0.95,1.2*dx,2*dy],'string','TMDs:'); uicontrol('pos',[0.74,0.95,0.8*dx,2*dy],'string','','tag','p2'); uicontrol('pos',[0.81,0.95,1.5*dx,2*dy],'string','TMDn:'); uicontrol('pos',[0.86,0.95,0.8*dx,2*dy],'string','','tag','p3'); lef=lef-0.03; j=0.96; if get(findobj(findobj('tag','mainmanu'),'tag','ch_1dof'),'value'); uicontrol('style','frame','pos',[lef-0.015,j-7*dy,0.15,7*dy]); uicontrol('pos',[lef-0.01,j-1*dy,2.6*dx,2*dy],'string','HƯ kh«ng cã TMD');j=j-3*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x1', 'tag','x1','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_x1', 'tag','A_x1','callback','result1');j=j-1.8*dy; Lo_sw=get(findobj(findobj('tag','mainmanu'),'tag','taitrong'),'value'); uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','dx1', 'tag','dx1','callback','result1') uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_dx1', 'tag','A_dx1','callback','result1');j=j-1.8*dy; if Lo_sw==2 | Lo_sw==3 uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x1f', 'tag','x1_f','callback','result1') uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','Ax1f', 'tag','A_x1_f','callback','result1'); end j=j-2*dy; LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 69 end if get(findobj(findobj('tag','mainmanu'),'tag','ch_2dof'),'value') uicontrol('style','frame','pos',[lef-0.015,j-11*dy,0.15,11*dy]); uicontrol('pos',[lef-0.01,j-1*dy,2.4*dx,2*dy],'string','TMD Năng lợng');j=j-3*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x1', 'tag','x1dk','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_x1', 'tag','A_x1dk','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','dx1', 'tag','dx1dk','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_dx1', 'tag','A_dx1dk','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x2', 'tag','x2dk','callback','result1') uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_x2', 'tag','A_x2dk','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','dx2', 'tag','dx2dk','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_dx2', 'tag','A_dx2dk','callback','result1');j=j-1.8*dy; Lo_sw=get(findobj(findobj('tag','mainmanu'),'tag','taitrong'),'value'); if Lo_sw==2 | Lo_sw==3 uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','xf', 'tag','x1dk_f','callback','result1') uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_xf', 'tag','A_x1dk_f','callback','result1'); end j=j-2*dy; end if get(findobj(findobj('tag','mainmanu'),'tag','ch_act'),'value') uicontrol('style','frame','pos',[lef-0.015,j-11*dy,0.15,11*dy]); uicontrol('pos',[lef-0.01,j-1*dy,2.8*dx,2*dy],'string','§iỊu khiĨn tÝch cùc');j=j-3*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x1', 'tag','x1dk_tr','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_x1', 'tag','A_x1dk_tr','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','dx1', 'tag','dx1dk_tr','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_dx1', 'tag','A_dx1dk_tr','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','x2', 'tag','x2dk_tr','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_x2', 'tag','A_x2dk_tr','callback','result1');j=j-1.8*dy; uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','dx2', LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 70 'tag','dx2dk_tr','callback','result1'); uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_dx2', 'tag','A_dx2dk_tr','callback','result1');j=j-1.8*dy; Lo_sw=get(findobj(findobj('tag','mainmanu'),'tag','taitrong'),'value'); if Lo_sw==2 | Lo_sw==3 uicontrol('style','checkbox','pos',[lef,j,dx,2*dy],'string','xf', 'tag','x1dk_ftr','callback','result1') uicontrol('style','checkbox','pos',[lef+1.2*dx,j,1.2*dx,2*dy],'string','A_xf', 'tag','A_x1dk_ftr','callback','result1'); end j=j-2*dy; end j=j-2*dy; uicontrol('style','pushbutton','string','Return', 'pos',[lef-0.01,j,2.5*dx,2.5*dy],'callback','Funct1(''return'')'); function toiuu if get(findobj(findobj('tag','mainmanu'),'tag','ch_opt'),'value') set(findobj(findobj('tag','mainmanu'),'tag','ch_opt'),'value',1); set(findobj(findobj('tag','mainmanu'),'tag','var5'),'enable','off','BackgroundColor',[1 0.8 0.8]); set(findobj(findobj('tag','mainmanu'),'tag','var6'),'enable','off','BackgroundColor',[1 0.8 0.8]); else set(findobj(findobj('tag','mainmanu'),'tag','var5'),'enable','on','BackgroundColor',[1 1]); set(findobj(findobj('tag','mainmanu'),'tag','var6'),'enable','on','BackgroundColor',[1 1]); end clc; mar='o.svp'; ff1='1-m^2*eta*gamma^3-3*m*eta*gamma-m*eta^2'; for i=2:2:10 ff2=subs(ff1,'m',0.01*i); ezplot(ff2,[0,120,0,120]);hold on; df=findobj('type','line'); set(df(1),'Marker',mar(i/2),'markersize',5); end df=findobj('type','line') set(df,'LineStyle',':'); legend('m=2%','m=4%','m=6%','m=8%','m=10%',3); ff3='1-m*eta*gamma'; for i=2:2:10 ff4=subs(ff3,'m',0.01*i); ezplot(ff4,[0,120,0,120]);hold on; end set(gca,'xlim',[0,10],'ylim',[0,10]); set(gcf,'color',[1 1]); title('Miền ổn định') LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ... Chương 2: Tính tốn tham số tối ưu TMD theo phương pháp biến đổi Laplace 15 2.1 Chuyển động hệ chưa có TMD 15 2.2 Tính tốn tham số tối ưu TMD phương pháp biến đổi Laplace ... cậy phương pháp luận văn Mục đích luận văn Mục đích luận văn xây dựng phương pháp tính tốn mới, xác định tham số đặc trưng hệ Chứng minh tính đắn phương pháp giải tìm nghiệm Bố cục luận văn Luận. .. Nhận thấy hai tham số tìm cách sử dụng phương pháp phép biến đổi Laplace trùng với tham số tìm phương pháp cân cực Việc sử dụng phép biến đổi Laplace để tính tốn thiết kế thơng số cho TMD chưa có