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 ĐẠ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 ĐẠ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 ĐẠ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 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 .3 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 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 2.2.1 Thiết lập hệ phương trình vi phân chuyển động hệ 17 2.2.2 tối ưu TMD 2.2.3 hệ phương trình 2.3 17 Tính tốn tham số 19 Giải tìm nghiệm 27 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 hệ chưa lắp TMD Chuyển động 34 3.1.2 hệ lắp thêm TMD 3.2 Chuyển động 35 So sánh với số phương pháp khác 37 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 CÁC KÝ HIỆU VIẾT TẮT M U Trong trình làm vi c d i tác thi t b , k t c u s ng nh h ti t c ng nh b ng i ta có th ng Tuy nhiên khơng ph i ng pháp vì, ho c ph i chi phí t khơng th kh Trong vài th p k d ng thi t b tiêu tán n ng l thi t b l ng dao va ch m có h i c a h nh có th it c l p vào h b ng ph m vi nghiên c Nghiên c u nhi u giai o n Trong ó có - Nghiên c u dao - Nghiên c u dao T ó i v i t ng cơng trình, ta có th áp d ng t ng d ng TBTTNL m t cách phù h p Vì ki n trúc c a k t c u r t a d ng nên hình dáng c u trúc c a thi t b c ng không ph n phong phú Trong lu n v n ta s xét n s dao ng c a h m t 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 ó tìm tham s t i u t s c n t s ng phơ ng trình chuy n ng, t t n s t i u 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 c a phơ ng pháp lu n v n M c ích c a lu n v n M c ích c a 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 c a h Ch ng minh tính úng n c a phơ ng pháp gi i tìm nghi m B c c c a lu n v n Lu n v n bao g m ph n m - Chơ ng 1: Gi i thi u thi t b tiêu tán n ng l ng kh i l ng phơ ng pháp bi n i Laplace - Chơ ng 2: Tính tốn tham s Laplace - Chơ ng 3: Mô ph ng dao 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 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(i-3),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 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]); End % set(findobj(findobj('tag','mainmanu'),'tag','var9'),'enable','on','BackgroundColor',[1 1]); 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); set(findobj(findobj('tag','mainmanu'),'Tag','var2'),'String',15000); % 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); set(findobj(findobj('tag','mainmanu'),'Tag','var4'),'String',55000); % k1 % d1 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]); set(findobj(findobj('tag','mainmanu'),'Tag','var7'),'String',5e+5); set(findobj(findobj('tag','mainmanu'),'Tag','var8'),'String',0.1); set(findobj(findobj('tag','mainmanu'),'Tag','var9'),'String',0.01); set(findobj(findobj('tag','mainmanu'),'Tag','var10'),'String',0); set(findobj(findobj('tag','mainmanu'),'Tag','var11'),'String',30); 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); 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); i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); for 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); i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); for 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); 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); i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); for 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); i=1:15 for 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); i=1:15 set(findobj(findobj('tag','mainmanu'),'Tag',['var' num2str(i)]),'enable','on','BackgroundColor',[1 1]); for end 0.8]); set(findobj(findobj('tag','mainmanu'),'Tag','var9'),'enable','off','BackgroundColor',[1 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 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, 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; 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=j3*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',[lef0.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', 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 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') 1]); ... 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. .. ĐẠ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... ĐẠ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