Tính lựa chọn cấp bảo vệ dùng phần mềm matlab

Chương 4 TÍNH LỰA CHỌN CẤP BẢO VỆ DÙNG PHẦN MỀM MATLAB I . GIỚI THIỆU PHẦN MỀM MATLAB : Định nghĩa : MATLAB là môi trường tính toán kỹ thuật bằng máy tính. MATLAB tập hợp các phép phân tích số học, tính toán ma trận, xử lý tín hiệu và đồ hoạ trong môi trường dễ sử dụng trong đó các vấn đề và giải thuật được diễn tả bằng các biểu thức toán học mà không cần phải lập trình phức tạp như các ngôn ngữ khác . MATLAB cũng có những công cụ ứng dụng trong một số lĩnh vực gọi là các hộp công cụ (toolbox) rất quan trọng đối với hầu hết người sử dụng MATLAB. Các hộp công cụ bao gồm xử lý tín hiệu, thiết kế hệ thống điều khiển, tái tạo hệ thống động lực, nhận dạng hệ thống và các lĩnh vực khác. Cài đặt MATLAB vào máy tính: a. Yêu cầu của phần mềm : Phần mềm MATLAB được chứa trong đĩa CD hay bao gồm nhiều đĩa 1,4 MB: đĩa hệ thống và các đĩa chứa các hộp thoại công cụ. Do đó để cài đặt MATLAB, cần một máy tính với các yêu cầu sau : Máy tính IBM, Compaq hay hệ thống máy tính tương thích 100 % với một CPU Intel 386, 486, hay CPU hiệu Pentium. Một đĩa cứng với ít nhất 30 MB còn trống. Bộ nhớ trong phải tối thiểu là 8 MB. Một màn hình màu. Phần mềm Microsoft Window 3.X trở lên. Hệ điều hành MS-DOS version 3.1 hoặc mới hơn. Một mouse sử dụng được trong Window. b. Các bước cài đặt : Để máy tính sử dụng được MATLAB thì cấu hình tối thiểu bộ nhớ RAM phải có dung lượng 16 MB hay 32 MB, ổ cứng 1,2 GB trở lên, ổ đĩa CD để Install chương trình từ CD – ROM, nếu không có ổ CD thì nạp chương trình này vào máy từ ổ đĩa cứng phụ đã nạp sẵn phần mềm MATLAB. Khi để đĩa CD vào máy, chương trình Auto Run tự động chạy. Nhấp chuột vào INSTALL, hộp thoại Setup xuất hiện. Một vệt xanh xuất hiện tăng dần từ trái sang phải báo cho biết việc cài đặt đang tiến hành. Khi kết thúc quá trình thì hộp thoại Setup Type xuất hiện. Chương trình MATLAB cho 3 chế độ cài đặt: Typical – cài đặt những phần chung nhất thường được nhiều người sử dụng. Compact – cài đặt những phần hệ tối thiểu. Custom – cho phép người dùng chọn theo ý riêng. Vì MATLAB 5.3 chiếm dung lượng rất lớn, ở đây chọn Custom để chọn lựa các phần cần thiết cho phù hợp theo dung lượng ổ đĩa cứng đang còn trống. Nhấp chuột vào phần Custom – Destination Directory – mặc nhiên MATLAB chọn đường dẫn thư mục cài đặt ở C:\ program Files\ Odyssey\ -nếu đồng ý thì nhấp vào Next để tiếp tục. Hộp thoại Welcome xuất hiện. Nếu không đồng ý thì nhấp vào nút Browse để chỉ đường dẫn khác. Hộp thoại Choose Directory xuất hiện. Nhập đường dẫn thư mục cài đặt trong hộp Path ,nhấn OK tiếp tục cài đặt. Nhấn Next để tiếp tục cho đến khi hoàn thành quá trình cài đặt. c. Khởi động chương trình MATLAB: Trước tiên khởi động Windows, sau đó kích Mouse vào biểu tượng MATLAB để kích hoạt chương trình và cửa sổ vùng làm việc (workspace)của MATLAB. Từ đây có thể bắt đầu làm việc với chương trình MATLAB. II.TÍNH TOÁN LỰA CHỌN CẤP BẢO VỆ : Những tác hại do sét đánh trực tiếp gây ra, gây thiệt hại lớn về ngừơi và tàisản, do đó phải yêu cầu lựa chọn cấp bảo vệ. Trong hầu hết các trường hợp, phải cân nhắc các yếu tố sau: Công trình có nhiều người làm việc. Dịch vụ trong công trình hoạt động liên tục. Mật độ sét tại vùng xây dựng công trình. Chiều cao công trình. Công trình có hay không có chứa vật liệu cháy nổ, các di sản văn hoá không thể thay thế được. Tùy theo đặc điểm của công trình mà có thể lựa chọn cấp bảo vệ thích hợp. Tuỳ theo tỷ số giữa tần số sét chấp nhận được Nc và tần số sét tính toán Nd mà ta có thể quyết định lựa chọn cấp bảo vệ phù hợp. Mục đích: Tính toán lựa chọn cấp bảo vệ chống sét. Môi trường tính toán: Môi trường tính toán chống sét là phần mềm MATLAB 5.3. Các số liệu ban đầu: Mật độ phóng điện sét Ng là mật độ phóng điện sét hằng năm trên một km2 . Ng = Na / 2,2 Ngmax = 0,04 . Nk1,25 = Nk / 10 = 2Ng Trong đó : Ngmax : mật độ phóng điện sét lớn nhất Nk : số ngày dông trung bình trong năm ở một số địa phương của Việt Nam do tổng cục khí tượng thuỷ văn thống kê. Tần số sét đánh trực tiếp Nd : Trong một năm, tần số này được tính như sau: Nd = Ngmax . Ac . C1 . 10-6 ( / năm ) Trong đó : Ac : vùng tập trung tương đương của công trình, là vùng có xác suất bị sét đánh giống như công trình ở cùng một điều kiện sét. Đối với công trình có cấu trúc hình hộp chữ nhật với chiều dài là L, chiều rộng là W, chiều cao là H, thì vùng tập trung tương đương được tính: Ac = LW + 6H( L + W ) + 9pH2 C1 : hệ số phụ thuộc vào vị trí công trình Vị trí công trình C1 Công trình nằm trong vùng có các cấu trúc khác, cây cao ngang hoặc cao hơn công trình. 0,25 Công trình bao quanh bởi các công trình khác thấp hơn. 0,5 Công trình độc lập, không có các công trình khác trong phạm vi 3H. 1 Công trình nằm trên núi hay đồi cao. 1 Khi vùng tập trung tương đương của công trình bao cho cả những công trình lân cận thì các công trình lân cận khỏi cần phải tính toán kỹ. Tần số sét chấp nhận được Nc : Nc = 5,5 / C Trong đó : C = C2 . C3 . C4 . C5 Với : C2 : hệ số cấu trúc công trình C2 Kết cấu Kim loại Thường Dễ cháy Kim loại 0,5 1 2 Thường 1 1 2,5 Dễ cháy 2 2,5 3 C3 : hệ số phụ thuộc vào vật liệu chứa trong công trình C3 Không giá trị và không dễ cháy. 0,5 Giá trị bình thường và có khả năng bắt lửa. 1 Giá trị cao và đặc biệt bắt lửa. 2 Giá trị cao, bắt lửa cao, dễ gây nổ. 3 C4 : hệ số phụ thuộc tính năng của công trình C4 Không có người làm việc thường xuyên 0,5 Có người làm việc thường xuyên 1 Tập trung đông người, khó sơ tán 3 C5 : hệ số phụ thuộc vào tác hại của sét C5 Công trình không yêu cầu hoạt động liên tục, không có tác hại đối với môi trường 1 Công trình yêu cầu hoạt động liên tục ,có tác hại đối với môi trường 5 Có tác hại đối với môi trường 10 Nếu Nd Nc : công trình có thể không cần hệ thống chống sét. Nếu Nd > Nc : công trình cần thiết phải có hệ thống chống sét với hệ số E = 1 – Nc / Nd cho phép lựa chọn cấp bảo vệ cho công trình. Hệ số E Cấp bảo vệ lựa chọn Dòng xung đỉnh ( KA) Khoảng cách tiên đạo ( m ) E > 0,98 Cấp 1 + biện pháp bảo vệ bổ sung - - 0,95 [ E 0,98 Cấp 1 2,8 20 0,8 [ E < 0,95 Cấp 2 9,5 45 0 < E < 0,8 Cấp 3 14,7 60 Khi một hệ thống chống sét có hệ số E’ nhỏ hơn hệ số tính toán E thì phải sử dụng một số biện pháp bổ sung để tăng cường hiệu quả của các hệ thống, các biện pháp tiêu biểu là : Giảm điện áp bước và điện áp tiếp xúc. Giảm ảnh hưởng quá điện áp xung lên các thiết bị điện tử. Bắt đầu Lưu đồ: Nhập các thông số L, W, H, C1, C2, C3, C4, C5, Ngm Nc:=5.5/C2*C3*C4*C5 Ac:=L*W+6*H*(L+W)+9*p*H² Nd:=Ac*C1*Ngm*10-² Nd>Nc S Đ E:=1-Nc/Nd E:=1-Nc/Nd E:=1-Nc/Nd E>0.98 S 0.95£E £0.98 Đ S Cấp 1 + bổ sung Không cần bảo vệ 0.8£E<0.95 Đ Cấp 1 Đ S Cấp 2 Cấp 3 Kết thúc Hướng dẫn sử dụng: Khởi động MATLAB 5.3. Chuyển thư mục hiện hành sang ổ đĩa A bằng cách đánh c:\>cd A:\ ¿ Sau đó gọi file spkt bằng cách ngập lệnh spkt tại dấu nhắc của chương trình : A:\>spkt¿ Giao diện của chương trình xuất hiện như hình ve.õ Nhập mật khẩu " anh " vào ô mật khẩu và nhấn phím Enter (¿) Cửa sổ làm việc của chương trình xuất hiện như hình vẽ : Nhập các thông số cần thiết vào: chiều dài L, chiều rộng W, chiều cao H, mật độ phóng điện sét trong một năm trên một km Ngm, các hệ số của công trình: C1, C2, C3, C4, C5. Đưa con trỏ tới mục tính toán, ấn phím Enter (¿ ). Chương trình tự động tính toán và đưa ra kết quả. Các kết quả: Nhập vào các thông số: dài L=68m, rộng W=40m, cao H=11m, C1=1, C2=2, C3=3, C4=3, C5=10. Kết quả tính toán được hiển thị như sau: Thay đổi các thông số cần tính toán, dài L=44m, rộng W=25m, cao H=11m, c1=1, c2=2, C3=3, C4=3, C5=10 ta được kết quả như sau : Mã nguồn: Demo function demo() t=(1:2:15)'*pi/8; x=sin(t); y=cos(t); fill(x,y,'r'); axis equal; text(0,0,'CHỐNG SÉT','color',[1 1 1],... 'fontname','Vni-times','fontsize',40,... 'HorizontalAlignment','center'); Giao diện: *Tao cac giao dien pw=get(h1_Epw,'string'); if strcmp(pw,'anh') set(h1_Title1,'visible','off'); set(h1_CT,'visible','off'); set(h1_HCG,'visible','off'); set(h1_TPW,'visible','off'); set(h1_Epw,'visible','off'); set(h1_THE,'visible','on'); set(h1_Frame1,'visible','on'); set(h1_Frame2,'visible','on'); set(h1_Frame3,'visible','on'); set(h1_Frame4,'visible','on'); set(h1_Tlwh,'visible','on'); set(h1_TL,'visible','on'); set(h1_TW,'visible','on'); set(h1_TH,'visible','on'); set(h1_Ngmax,'visible','on'); set(h1_EL,'visible','on'); set(h1_EW,'visible','on'); set(h1_EH,'visible','on'); set(h1_ENgmax,'visible','on'); set(h1_TC,'visible','on'); set(h1_TC1,'visible','on'); set(h1_TC2,'visible','on'); set(h1_TC3,'visible','on'); set(h1_TC4,'visible','on'); set(h1_TC5,'visible','on'); set(h1_PC1,'visible','on'); set(h1_PC2,'visible','on'); set(h1_PC3,'visible','on'); set(h1_PC4,'visible','on'); set(h1_PC5,'visible','on'); set(h1_TKQ,'visible','on'); set(h1_TNd,'visible','on'); set(h1_TNc,'visible','on'); set(h1_TE,'visible','on'); set(h1_ENd,'visible','on'); set(h1_ENc,'visible','on'); set(h1_EE,'visible','on'); set(h1_BAction,'visible','on'); demo; end spkt load SPKT h0 = figure('Units','points', ... 'Color',[0 0 0], ... 'Colormap',mat0, ... 'FileName','SPKT', ... 'MenuBar','none', ... 'Name','HE CHUYEN GIA CHONG SET', ... 'NumberTitle','off', ... 'PaperPosition',[18 180 576 432], ... 'PaperUnits','points', ... 'Position',[0 21 480 324.75], ... 'Tag','Fig1', ... 'ToolBar','none',... 'DefaultaxesColor','none', ... 'DefaultaxesXColor',[1 0 1], ... 'DefaultaxesYColor',[1 0 1], ... 'DefaultaxesZColor',[1 0 1]); h1_axes = axes('Parent',h0, ... 'Units','pixels', ... 'CameraUpVector',[0 1 0], ... 'Color',[0 0 0], ... 'ColorOrder',mat1, ... 'Position',[232 45 386 314], ... 'Tag','Axes1', ... 'visible','off',... 'XColor',[1 0 1], ... 'YColor',[1 0 1], ... 'ZColor',[1 0 1]); h1_Title1 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0 0 0], ... 'FontName','Vni-times', ... 'FontSize',20, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[37.5 242.25 405.75 59.25], ... 'String','TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT KHOA ĐIỆN ĐIỆN TỬ', ... 'Style','text', ... 'Tag','StaticText1'); h1_CT = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0 0 0], ... 'FontName','Vni-times', ... 'FontSize',15, ... 'ForegroundColor',[0 1 0], ... 'ListboxTop',0, ... 'Position',[43.5 198 96.75 19.5], ... 'String','Chương trình', ... 'Style','text', ... 'Tag','StaticText2'); h1_HCG = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0 0 0], ... 'FontName','Vni-times', ... 'FontSize',30, ... 'ForegroundColor',[0 0 1], ... 'ListboxTop',0, ... 'Position',[57.75 129 395.25 37.5], ... 'String','XÁC ĐỊNH MỨC BẢO VỆ CHỐNG SÉT', ... 'Style','text', ... 'Tag','StaticText3'); h1_TPW = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0 0 0], ... 'FontName','Vni-times', ... 'FontSize',15, ... 'ForegroundColor',[0 1 0], ... 'ListboxTop',0, ... 'Position',[124.5 54.75 102 21.75], ... 'String','Nhập mật khẩu', ... 'Style','text', ... 'Tag','StaticText4'); h1_Epw = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[233.25 56.25 84.75 18], ... 'Style','edit', ... 'callback','giaodien',... 'Tag','EditText1'); h1_THE = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0 1 0], ... 'FontName','Vni-times', ... 'FontSize',18, ... 'ForegroundColor',[0.501960784313725 0 0.501960784313725], ... 'ListboxTop',0, ... 'Position',[162.75 287.25 306.75 25.5], ... 'String','XÁC ĐỊNH MỨC BẢO VỆ', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText5'); h1_Frame1 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.501960784313725], ... 'ListboxTop',0, ... 'Position',[3 7.5 150 315.75], ... 'Style','frame', ... 'visible','off',... 'Tag','Frame1'); h1_Frame2 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'ListboxTop',0, ... 'Position',[9 222.75 138.75 79.5], ... 'Style','frame', ... 'visible','off',... 'Tag','Frame2'); h1_Tlwh = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.501960784313725], ... 'FontName','Vni-times', ... 'FontSize',11, ... 'ForegroundColor',[0 1 0.250980392156863], ... 'ListboxTop',0, ... 'Position',[12.75 304.5 129 16.5], ... 'String','Nhập các thông số', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText6'); h1_TL = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 288 64.5 12], ... 'String','Thông số L', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText7'); h1_TW = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 269.25 64.5 12], ... 'String','Thông số W', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText7'); h1_TH = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 249 64.5 12], ... 'String','Thông số H', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText7'); h1_Ngmax = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[11.25 226.5 67.5 15.75], ... 'String','Thông số Ngm', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText7'); h1_EL = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[81.75 285.75 63.75 15], ... 'Style','edit', ... 'string','0',... 'visible','off',... 'Tag','EditText2'); h1_EW = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[81 267 63.75 15], ... 'Style','edit', ... 'string','0',... 'visible','off',... 'Tag','EditText2'); h1_EH = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[81 247.5 63.75 15], ... 'Style','edit', ... 'string','0',... 'visible','off',... 'Tag','EditText2'); h1_ENgmax = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[81 227.25 63.75 15], ... 'Style','edit', ... 'visible','off',... 'string','0',... 'Tag','EditText2'); h1_TC = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.501960784313725], ... 'FontName','Vni-times', ... 'FontSize',11, ... 'ForegroundColor',[0 1 0], ... 'ListboxTop',0, ... 'Position',[21 202.5 112.5 16.5], ... 'String','Chọn các thông số C', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText8'); h1_Frame3 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'ListboxTop',0, ... 'Position',[8.25 110.25 138.75 90.75], ... 'Style','frame', ... 'visible','off',... 'Tag','Frame3'); h1_TC1 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 183.75 60.75 13.5], ... 'String','Chọn C1', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText9'); h1_TC2 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[11.25 165.75 60.75 14.25], ... 'String','Chọn C2', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText9'); h1_TC3 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[11.25 149.25 60.75 12.75], ... 'String','Chọn C3', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText9'); h1_TC4 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 132 60.75 13.5], ... 'String','Chọn C4', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText9'); h1_TC5 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.250980392156863], ... 'FontName','Vni-times', ... 'FontSize',10, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[10.5 114.75 60.75 14.25], ... 'String','Chọn C5', ... 'Style','text', ... 'visible','off',... 'Tag','StaticText9'); h1_PC1 = uicontrol('Parent',h0, ... 'Units','points', ... 'ListboxTop',0, ... 'Position',[81 184.5 63.75 15.75], ... 'String','0.25|0.5|1.0|1.0', ... 'Style','popupmenu', ... 'Tag','PopupMenu1', ... 'visible','off',... 'Value',1); h1_PC2 = uicontrol('Parent',h0, ... 'Units','points', ... 'ListboxTop',0, ... 'Position',[81 166.5 63.75 15.75], ... 'String','0.5|1.0|2.0', ... 'Style','popupmenu', ... 'Tag','PopupMenu1', ... 'visible','off',... 'Value',1); h1_PC3 = uicontrol('Parent',h0, ... 'Units','points', ... 'ListboxTop',0, ... 'Position',[80.25 148.5 63.75 15.75], ... 'String','0.5|1.0|2.0|3.0', ... 'Style','popupmenu', ... 'Tag','PopupMenu1', ... 'visible','off',... 'Value',1); h1_PC4 = uicontrol('Parent',h0, ... 'Units','points', ... 'ListboxTop',0, ... 'Position',[80.25 130.5 63.75 15.75], ... 'String','0.5|1.0|3.0', ... 'Style','popupmenu', ... 'Tag','PopupMenu1', ... 'visible','off',... 'Value',1); h1_PC5 = uicontrol('Parent',h0, ... 'Units','points', ... 'ListboxTop',0, ... 'Position',[80.25 112.5 63.75 15.75], ... 'String','1.0|5.0|10.0', ... 'Style','popupmenu', ... 'Tag','PopupMenu1', ... 'visible','off',... 'Value',1); h1_TKQ = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.501960784313725 0.501960784313725 0.501960784313725], ... 'FontName','Vni-times', ... 'FontSize',11, ... 'ForegroundColor',[0 1 0], ... 'ListboxTop',0, ... 'Position',[37.5 96 73.5 13.5], ... 'String','Các kết quả', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText10'); h1_Frame4 = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 0], ... 'ListboxTop',0, ... 'Position',[8.25 11.25 138.75 81.75], ... 'Style','frame', ... 'visible','off',... 'Tag','Frame4'); h1_TNd = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 0], ... 'ForegroundColor',[1 0 0.501960784313725], ... 'ListboxTop',0, ... 'Position',[14.25 75 53.25 12.75], ... 'String','Nd :', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText11'); h1_TNc = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 0], ... 'ForegroundColor',[1 0 0.501960784313725], ... 'ListboxTop',0, ... 'Position',[13.5 56.25 53.25 12.75], ... 'String','Nc :', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText11'); h1_ENd = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[80.25 72.75 64.5 15], ... 'Style','edit', ... 'visible','off',... 'string','0',... 'Tag','EditText3'); h1_ENc = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[81 54.75 64.5 15], ... 'Style','edit', ... 'visible','off',... 'string','0',... 'Tag','EditText3'); h1_TE = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 0], ... 'ForegroundColor',[1 0 0.501960784313725], ... 'ListboxTop',0, ... 'Position',[13.5 36 53.25 12.75], ... 'String','E :', ... 'visible','off',... 'Style','text', ... 'Tag','StaticText11'); h1_EE = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[1 1 1], ... 'ListboxTop',0, ... 'Position',[80.25 36.75 64.5 15], ... 'Style','edit', ... 'visible','off',... 'string','0',... 'Tag','EditText3'); h1_BAction = uicontrol('Parent',h0, ... 'Units','points', ... 'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ... 'FontName','Vni-times', ... 'FontSize',15, ... 'ForegroundColor',[1 0 0], ... 'ListboxTop',0, ... 'Position',[9 13.5 137.25 21], ... 'String','Tính toán', ... 'callback','TinhToan,set(h1_axes,''visible'',''on'')',... 'visible','off',... 'Tag','Pushbutton1'); tính toan %Nhap cac thong so tinh toan %Tinh Nd L=get(h1_EL,'string'); L=str2num(L); W=get(h1_EW,'string'); W=str2num(W); H=get(h1_EH,'string'); H=str2num(H); Ngmax=get(h1_ENgmax,'string'); Ngmax=str2num(Ngmax); Ac=L*W+6*H*(L+W)+9*pi*H^2; %Chãn C1 choiceC1=get(h1_PC1,'value'); if choiceC1==1,C1=0.25 ; elseif choiceC1==2,C1=0.5 ; elseif choiceC1==3,C1=1.0; elseif choiceC1==4,C1=3.0 ; end Nd=Ngmax*Ac*C1*10^(-6); set(h1_ENd,'string',Nd); %TÛnh Nc %Chãn C2 choiceC2=get(h1_PC2,'value'); if choiceC2==1,C2=0.5 ; elseif choiceC2==2,C2=1.0; elseif choiceC2==3,C2=2.0 ; end %Chãn C3 choiceC3=get(h1_PC3,'value'); if choiceC3==1,C3=0.5 ; elseif choiceC3==2,C3=1.0; elseif choiceC3==3,C3=2.0 ; elseif choiceC3==4,C3=3.0 ; end %Chãn C4 choiceC4=get(h1_PC4,'value'); if choiceC4==1,C4=0.5 ; elseif choiceC4==2,C4=1.0; elseif choiceC4==3,C4=3.0 ; end %Chãn C5 choiceC5=get(h1_PC5,'value'); if choiceC5==1,C5=1.0 ; elseif choiceC5==2,C5=5.0; elseif choiceC5==3,C5=10.0 ; end C=C2*C3*C4*C5; Nc=5.5/C; set(h1_ENc,'string',Nc); if Nd<=Nc,set(h1_THE,'FontName','Vni-times', ... 'FontSize',15,'string','KHÔNG CẦN BẢO VỆ'); else E=1-Nc/Nd; set(h1_EE,'string',E); if (E>0.98),set(h1_THE,'FontName','Vni-times', ... 'FontSize',15,'string','BẢO VỆ CẤP 1 + BỔ SUNG'); elseif (E>=0.95)&(E<=0.98),set(h1_THE,'FontName','Vni-times', ... 'FontSize',15,'string','BẢO VỆ CẤP 1'); elseif (E>=0.8)&(E<0.95),set(h1_THE,'FontName','Vni-times', ... 'FontSize',15,'string','BẢO VỆ CẤP 2'); elseif (E>0)&(E<0.8),set(h1_THE,'FontName','Vni-times', ... 'FontSize',15,'string','BẢO VỆ CẤP 3'); end end ve(L,W,H); ve function ve(L,W,H) hh=H-2; %ChiỊu cao H subplot(6,1.52,3.1); xo=[0 0];yo=[H+2 2*H+2];plot(xo,yo,'--'); text(0,2*H,' \leftarrowH','fontsize',8,'color','r'); hold on %ChiỊu dªi L x1=[0 L];y1=[2*H+2 2*H+2];plot(x1,y1,'--'); text(L/2,2*H-1,' \downarrow L','fontsize',8,'color','r'); x2=[L L];y2=[2*H+2 H+2];plot(x2,y2,'--'); x3=[L 0];y3=[H+2 H+2];plot(x3,y3,'--'); %w Lª chiỊu rèng x4=[0 W];y4=[2*H+2 2.25*H+2];plot(x4,y4,'--'); x5=[W L+W];y5=[2.25*H+2 2.25*H+2];plot(x5,y5,'--'); x6=[L+W L];y6=[2.25*H+2 2*H+2];plot(x6,y6,'--'); text(L+1,2*H+1,'\uparrow W','fontsize',8,'color','r'); x7=[L+W L+W];y7=[2.25*H+2 1.25*H+2];plot(x7,y7,'--'); x8=[L+W L];y8=[1.25*H+2 H+2];plot(x8,y8,'--'); title('Kích thước toà nhà','Fontname','Vni-times',... 'fontsize',15,'color','y'); %axis off hold off global h subplot(2,1.65,3.3); h=rectangle('Position',[0,0,(6*H)+W,(6*H)+L],... 'Curvature',[0.8,0.4],... 'LineWidth',2,'FaceColor','r','EdgeColor','r'); title('Diện tích tương đương AC','Fontname','Vni-times',... 'fontsize',15,'color','y');