Tính toán điện trở nối đất tự nhiên của hệ thống dây chống sét – cột điện

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC LẠC HỒNG  TÍNH TỐN ĐIỆN TRỞ NỐI ĐẤT TỰ NHIÊN CỦA HỆ THỐNG DÂY CHỐNG SÉT – CỘT ĐIỆN LUẬN VĂN THẠC SĨ KỸ THUẬT ĐIỆN ĐỒNG NAI, NĂM 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC LẠC HỒNG  TÍNH TỐN ĐIỆN TRỞ NỐI ĐẤT TỰ NHIÊN CỦA HỆ THỐNG DÂY CHỐNG SÉT – CỘT ĐIỆN CHUYÊN NGÀNH: KỸ THUẬT ĐIỆN MÃ SỐ: 8520201 LUẬN VĂN THẠC SĨ KỸ THUẬT ĐIỆN ĐỒNG NAI, NĂM 2020 LỜI CẢM ƠN Trong q trình hồn thành đề tài nghiên cứu cách thuận lợi, tác giả nhận giúp đỡ, hướng dẫn nhiệt tình Thầy/Cơ giáo trường Đại Học Lạc Hồng Với tình cảm chân thành, tác giả bày tỏ lòng biết ơn trực tiếp đến Thầy TS Nguyễn Thanh Sơn phòng Sau Đại Học, Khoa Kỹ Thuật Điện – Trường Đại Học Thầy/Cô giáo tham gia quản lý, giảng dạy giúp đỡ cho tác giả suốt trình học tập nghiên cứu Mặc dù có nhiều cố gắng suốt trình thực đề tài, song cịn có mặt hạn chế, thiếu sót Tác giả mong nhận ý kiến đóng góp dẫn thầy cô giáo Cuối cùng, tác giả xin gửi lời biết ơn sâu sắc đến cha, mẹ ln quan tâm chăm sóc động viên học viên Chân thành cảm ơn bạn học viên khích lệ tinh thần để học viên theo đuổi hồn thành tốt khóa học Xin chân thành cảm ơn! Đồng Nai, ngày 31 tháng 12 năm 2020 LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu riêng tác giả Các số liệu, kết nêu trong luận văn trung thực chưa công bố nơi khác Các thơng tin trích dẫn nguồn gốc rõ ràng Kết tính tốn dựa tiêu chuẩn xây dựng hành Nếu không điều nêu trên, tác giả xin hoàn toàn chịu trách nhiệm đề tài Đồng Nai, ngày 31 tháng 12 năm 2020 TĨM TẮT LUẬN VĂN Tính toán điện trở nối đất tự nhiên hệ thống dây chống sét cột điện luôn nhà thiết kế đường dây truyền tải, trạm biến áp cao quan tâm Hệ thống nối đất nầy góp phần vào việc làm giảm nhỏ giá trị nối đất an toàn trạm biến áp Luận văn xây dựng đề xuất số cơng thức tính điện trở nối đất tự nhiên hệ thống dây chống sét cột điện xác so với cơng trình cơng bố [1, 2, 3, ] Đồng thời có khảo sát đến ảnh hưởng loại dây chống sét mà cơng trình nêu khơng đề cập đến Sau đó, so sánh kết báo với kết công bố MỤC LỤC TRANG PHỤ BÌA LỜI CẢM ƠN LỜI CAM ĐOAN TĨM TẮT LUẬN VĂN MỤC LỤC DANH MỤC HÌNH VẼ DANH MỤC BẢNG BIỂU CHƯƠNG MỞ ĐẦU 1 Đặt vấn đề: Nhiệm vụ luận án: Các bước tiến hành: Điểm luận án ý nghĩa khoa học đề tài Giá trị thực tiễn đề tài: Nội Dung Luận Án CHƯƠNG : CÁC VẤN ĐỀ VỀ AN TOÀN TRONG HỆ THỐNG ĐIỆN I Tổng Quan Về Hệ Thống Nối Đất : Hệ Thống Nối Đất : 1.1 Nối Đất Làm Việc : 1.2 Nối Đất An Toàn Hay Nối Đất Bảo Vệ : 1.2.1 Đặc Tính Của Dòng Điện Chạy Trong Đất : 1.2.2 Tác Dụng Bảo Vệ Của Hệ Thống Nối Đất : 1.3 Nối Đất Chống Sét : Các Vấn Đề Về Thiết Kế Hệ Thống Nối Đất: .8 2.1 Nguyên Tắc Chung : 2.2 Hệ Thống Nối Đất Của Trạm Biến Áp : 10 2.3 Thiết Kế Hệ Thống Nối Đất Theo Tiêu Chuẩn IEEE: 12 2.4 Dòng Điện Sự Cố Chạy Qua Hệ Thống Nối Đất Của Trạm: 13 2.4.1 Ảnh Hưởng Của Dây Chống Sét Và Dây Trung Tính: 14 2.4.2 Ảnh Hưởng Của Cáp Và Hệ Thống ống Dẫn: 14 2.5 Vị Trí Sự Cố Gây Ra Giá Trị Cự Đại Của Dòng Sự Cố Chạy Qua Hệ Thống Nối Đất Của Trạm: .14 II Bảo Vệ Chống Sét Cho Các Phần Tử Hệ Thống Điện : .15 Bảo Vệ Chống Sét Đường Dây Tải Điện : 15 1.1 Quá Điện Áp Cảm Ứng Đường Dây : 16 1.2 Sét Đánh Trên Đường Dây Không Treo Dây Chống Sét : 16 1.3 Đối Với Đường Dây Có Treo Dây Chống Sét : .16 1.4 Khi Có Hai Dây Chống Sét : 17 Bảo Vệ Chống Sét Trạm Biến Áp : 17 Bảo Vệ Chống Sét Nhà Máy Điện : .20 CHƯƠNG : TÍNH TỐN ĐIỆN TRỞ NỐI ĐẤT TỰ NHIÊN CỦA HỆ THỐNG DÂY CHỐNG SÉT-CỘT ĐIỆN .23 I Sơ lý thuyết đường dây dài : .23 Lý thuyết đường dây dài : .23 1.1 Khái Niệm Về Mạch Thông Số Rải: 23 1.2 Các Thông Số Đơn Vị Của Đường Dây Dài : .26 1.3 Phương Trình Đường Dây Dài : 27 1.3.1 Nghiệm Của Phương Trình Đường Dây Dài Với Tác Động Sin : 29 1.3.2 Các thơng số sóng đường dây dài : 32 II Mơ Hình Tốn Hệ Truyền Tải Điện Đơn Dây Chống Sét : .34 Tính Tổng Trở Thevenin Nhìn Từ Đầu Đường Dây Với Mơ Hình N-Mạch Pi Có Các Thông Số ZS Và ZP Giống Nhau : 35 Tính Tổng Trở Thevenin Nhìn Từ Đầu Đường Dây Với Mơ Hình N-Mạch Pi Có Các Thơng Số ZS, ZP Giống Nhau Và Cuối Đường Dây Có Mắc Thêm Một Phần Tử Z1 41 Đề Xuất : .44 Tính Toán Gần Đúng: 45 CHƯƠNG : PHƯƠNG PHÁP TÍNH NHANH ĐIỆN TRỞ NỐI ĐẤT TỰ NHIÊN HỆ DÂY CHỐNG SÉT-CỘT ĐIỆN 48 Phương Pháp Tính : .48 1.1 Công Thức Tổng Quát : 48 1.2 Xác Định Khoảng Vượt Tới Hạn ( n ): 49 th Đề Xuất Biểu Thức Giải Tích Đơn Giản Xác Định Điện Trở Dây Chống Sét Cột Điện 56 2.1 Công Thức Đơn Giản Xác Định Điện Trở Dây Chống Sét-Cột Điện Khi n  n : 57 th 2.2 Công Thức Đơn Giản Xác Định Điện Trở Dây Chống Sét-Cột Điện Khi n  n : 63 th CHƯƠNG : KẾT LUẬN .71 Cơng Thức Tổng Qt Tính Tổng Trở Tương Đương Hệ Dây Chống Sét – Cột Điện 71 1.1 Tính Tổng Trở Thevenin Nhìn Từ Đầu Đường Dây Với Mơ Hình N-Mạch Pi Có Các Thông Số ZS Và ZP Giống Nhau 72 1.2 Tính Tổng Trở Thevenin Nhìn Từ Đầu Đường Dây Với Mơ Hình N-Mạch Pi Có Các Thơng Số ZS Và ZP Giống Nhau Và Cuối Đường Dây Có Mắc Thêm Một Phần Tử Z1 73 Công Thức Tính Tốn Nhanh Điện Trở Nối Đất Tự Nhiên Của Hệ Thống Dây Chống Sét – Cột Điện: .75 TÀI LIỆU THAM KHẢO PHỤ LỤC I,II,III, IV DANH MỤC HÌNH VẼ Hình 1.1 Đặc tính dịng điện chạy đất Hình 1.2 Giải thích vai trị bảo vệ nối đất Hình 1.3 Sơ đồ bảo vệ trạm biến áp 35 – 110kV 18 Hình 1.4 Sơ đồ bảo vệ trạm biến áp cụt 19 Hình 1.5 Sơ đồ bảo vệ chống sét nhà máy điện 20 Hình 1.6 Sơ đồ giải thích vai trị chống sét tụ điện C 20 Hình 1.7 Sơ đồ chống sét nhà máy điện có cáp xuất tuyến 22 Hình 2.1 Sơ đồ đường dây dài 26 Hình 2.2 Sơ đồ tương đương đoạn dây x 27 Hình 2.3 Mạch điện hình thang thông số tập trung thay cho dây chống sét điểm nối đất cột 34 Hình 2.4 Mơ hình mạch n- pi đường dây chống sét 35 Hình 2.5 Mơ hình mạch n- pi sau phân tích thành n - T 35 Hình 2.6 Mơ hình mạng cửa tương cho n mạng hình T 37 Hình 2.7 Mạng cửa tương đương 38 Hình 2.8 Mạch hình T tương đương n mạch hình T 38 Hình 2.9 Mạch điện tương đương đường dây 39 Hình 2.10 Mơ hình mạch hình pi đường dây chống sét có Z1 cuối đường dây 41 Hình 2.11 Mơ hình mạch hình pi sau phân tích thành mạch hình T 42 Hình 2.12 Mạch điện tương đương đường dây 42 Hình 2.13 Mạch tương đương sau biến đổi 42 Hình 3.1 Mơ hình đơn giản hệ thống dây chống sét – cột điện 48 Hình 3.2 Quan hệ tổng trở theo khoảng vượt Z  f (n) 49 Hình 3.3 Quan hệ nth theo điện trở cột điện Rc dây  C 50 52 Hình 3.4 Quan hệ nth theo điện trở cột điện Rc dây  C 70 53 Hình 3.5 Quan hệ nth theo điện trở cột điện Rc dây  C 90 54 Hình 3.6 Quan hệ nth theo điện trở cột điện Rc dây  C 90 56 Hình 4.1 Mạch điện hình thang thơng số tập trung thay cho dây chống sét điểm nối đất cột 71 Hình 4.2 Mơ hình mạch n- pi đường dây chống sét 72 Hình 4.3 Mơ hình mạch n- pi sau phân tích thành n - T 72 Hình 4.4 Mơ hình mạng cửa tương cho n mạng hình T 72 Hình 4.5 Mơ hình mạch hình pi đường dây chống sét có Z1 cuối đường dây 73 Hình 4.6 Mơ hình mạch hình pi sau phân tích thành mạch hình T 74 Hình 4.7 Mạch điện tương đương đường dây 74 Bảng 3.14.a Số liệu nhận từ (3.9) (3.1)Khi n  nth Khoảng vượt 250m, dây  C 90, Rcs = 0.4275 Rc  10, nth  27 Rc  5, nth  20 Z  (%) Zn Z Rc  15, nth  33  (%) Zn Z  (%) n Zn 5.4275 5.8934 8.5844 10.4275 11.1906 7.3180 15.4275 16.4204 6.4358 3.0300 3.1850 5.1144 5.5321 5.8278 5.3451 8.0329 8.4398 5.0658 2.3142 2.3769 2.7111 3.9892 4.1380 3.7281 5.6589 5.8779 3.8708 2.0095 2.0333 1.1871 3.2791 3.3604 2.4793 4.5363 4.6666 2.8736 1.8609 1.8660 0.2752 2.8969 2.9421 1.5589 3.9105 3.9915 2.0720 1.7837 1.7792 -0.2513 2.6737 2.6979 0.9039 3.5293 3.5804 1.4470 1.7422 1.7326 -0.5488 2.5371 2.5485 0.4492 3.2846 3.3165 0.9710 1.7195 1.7072 -0.7147 2.4512 2.4546 0.1392 3.1221 3.1413 0.6152 1.7070 1.6932 -0.8066 2.3961 2.3945 -0.0695 3.0117 3.0223 0.3528 10 1.7000 1.6855 -0.8573 2.3605 2.3555 -0.2090 2.9356 2.9403 0.1613 11 1.6962 1.6812 -0.8852 2.3372 2.3301 -0.3015 2.8826 2.8833 0.0227 12 1.6940 1.6788 -0.9004 2.3219 2.3135 -0.3628 2.8455 2.8433 -0.0772 13 1.6928 1.6774 -0.9088 2.3119 2.3026 -0.4032 2.8192 2.8151 -0.1489 14 1.6922 1.6767 -0.9134 2.3053 2.2954 -0.4297 2.8007 2.7951 -0.2001 15 1.6918 1.6763 -0.9159 2.3009 2.2906 -0.4472 2.7876 2.7810 -0.2367 16 1.6916 1.6761 -0.9173 2.2980 2.2875 -0.4587 2.7782 2.7709 -0.2628 17 1.6915 1.6759 -0.9181 2.2961 2.2854 -0.4663 2.7716 2.7638 -0.2813 18 1.6914 1.6759 -0.9185 2.2948 2.2840 -0.4712 2.7668 2.7587 -0.2945 19 1.6914 1.6758 -0.9187 2.2940 2.2831 -0.4745 2.7634 2.7550 -0.3039 20 1.6913 1.6758 -0.9188 2.2935 2.2825 -0.4766 2.7610 2.7525 -0.3105 21 2.2931 2.2821 -0.4780 2.7593 2.7506 -0.3153 22 2.2928 2.2819 -0.4789 2.7581 2.7493 -0.3186 23 2.2927 2.2817 -0.4796 2.7572 2.7484 -0.3210 24 2.2926 2.2816 -0.4799 2.7566 2.7477 -0.3227 25 2.2925 2.2815 -0.4802 2.7562 2.7472 -0.3239 26 2.2925 2.2815 -0.4804 2.7558 2.7469 -0.3248 2.2924 2.7556 2.7466 -0.3254 28 2.7554 2.7465 -0.3258 29 2.7553 2.7463 -0.3261 30 2.7553 2.7463 -0.3263 31 2.7552 2.7462 -0.3265 32 2.7552 2.7462 -0.3266 33 2.7551 2.7461 -0.3267 27 Rc  20, nth  34 n  (%) 2.2814 -0.4805 Rc  25, nth  37  (%) Z 20.4275 21.6147 5.8119 25.4275 26.7862 5.3436 30.4275 31.9415 4.9756 10.5332 11.0351 4.7646 13.0335 13.6196 4.4968 15.5336 16.1962 4.2653 7.3270 7.6070 3.8208 8.9946 9.3290 3.7182 10.6619 11.0461 3.6032 5.7900 5.9637 2.9998 7.0422 7.2551 3.0230 8.2937 8.5427 3.0015 4.9176 5.0311 2.3074 5.9220 6.0652 2.4185 6.9250 7.0958 2.4669 4.3746 4.4506 1.7385 5.2153 5.3147 1.9047 6.0537 6.1749 2.0016 4.0170 4.0684 1.2811 4.7426 4.8127 1.4764 5.4648 5.5524 1.6036 3.7726 3.8073 0.9198 4.4139 4.4636 1.1253 5.0502 5.1143 1.2685 3.6014 3.6244 0.6383 4.1790 4.2142 0.8413 4.7500 4.7971 0.9899 10 3.4794 3.4940 0.4215 4.0080 4.0326 0.6142 4.5283 4.5627 0.7607 11 3.3913 3.3999 0.2558 3.8817 3.8986 0.4340 4.3619 4.3869 0.5740 12 3.3271 3.3314 0.1300 3.7875 3.7986 0.2922 4.2357 4.2536 0.4228 13 3.2801 3.2812 0.0350 3.7167 3.7234 0.1811 4.1391 4.1516 0.3012 14 3.2454 3.2442 -0.0365 3.6632 3.6666 0.0945 4.0648 4.0731 0.2038 15 3.2198 3.2169 -0.0901 3.6225 3.6235 0.0273 4.0073 4.0123 0.1262 16 3.2008 3.1967 -0.1303 3.5915 3.5906 -0.0248 3.9626 3.9651 0.0644 17 3.1867 3.1816 -0.1603 3.5679 3.5656 -0.0651 3.9277 3.9283 0.0155 18 3.1763 3.1705 -0.1827 3.5498 3.5464 -0.0962 3.9005 3.8996 -0.0233 19 3.1685 3.1621 -0.1994 3.5359 3.5317 -0.1202 3.8792 3.8772 -0.0540 20 3.1626 3.1559 -0.2118 3.5253 3.5204 -0.1387 3.8626 3.8595 -0.0781 21 3.1583 3.1513 -0.2211 3.5171 3.5117 -0.1529 3.8495 3.8457 -0.0972 Zn Z  (%) Zn Zn Z Rc  30, nth  43 22 3.1551 3.1479 -0.2280 3.5108 3.5051 -0.1638 3.8392 3.8349 -0.1122 23 3.1527 3.1453 -0.2331 3.5060 3.5000 -0.1722 3.8311 3.8264 -0.1241 24 3.1509 3.1434 -0.2369 3.5023 3.4960 -0.1787 3.8248 3.8197 -0.1334 25 3.1495 3.1420 -0.2398 3.4994 3.4930 -0.1836 3.8198 3.8144 -0.1407 26 3.1485 3.1409 -0.2419 3.4972 3.4907 -0.1874 3.8159 3.8103 -0.1465 27 3.1478 3.1401 -0.2435 3.4956 3.4889 -0.1904 3.8128 3.8070 -0.1510 28 3.1472 3.1395 -0.2446 3.4943 3.4875 -0.1926 3.8103 3.8044 -0.1546 29 3.1468 3.1391 -0.2455 3.4933 3.4865 -0.1943 3.8084 3.8024 -0.1574 30 3.1465 3.1388 -0.2462 3.4925 3.4857 -0.1956 3.8069 3.8008 -0.1596 31 3.1463 3.1385 -0.2466 3.4919 3.4850 -0.1967 3.8057 3.7996 -0.1613 32 3.1461 3.1383 -0.2470 3.4914 3.4845 -0.1974 3.8048 3.7986 -0.1627 33 3.1460 3.1382 -0.2473 3.4911 3.4842 -0.1980 3.8041 3.7978 -0.1637 34 3.1459 3.1381 -0.2475 3.4908 3.4839 -0.1985 3.8035 3.7972 -0.1646 35 3.1458 3.1380 -0.2476 3.4906 3.4837 -0.1989 3.8030 3.7967 -0.1652 36 3.1457 3.1380 -0.2477 3.4904 3.4835 -0.1991 3.8027 3.7964 -0.1658 37 3.1457 3.1379 -0.2478 3.4903 3.4834 -0.1993 3.8024 3.7961 -0.1662 38 3.1457 3.1379 -0.2479 3.4902 3.4833 -0.1995 3.8022 3.7958 -0.1665 39 3.1457 3.1379 -0.2479 3.4901 3.4832 -0.1996 3.8020 3.7956 -0.1667 40 3.1456 3.1378 -0.2479 3.4901 3.4831 -0.1997 3.8018 3.7955 -0.1669 41 3.4901 3.4831 -0.1998 3.8017 3.7954 -0.1671 42 3.4900 3.4830 -0.1998 3.8016 3.7953 -0.1672 43 3.4900 3.4830 -0.1999 3.8016 3.7952 -0.1673 44 3.4900 3.4830 -0.1999 3.8015 3.7952 -0.1674 45 3.4900 3.4830 -0.1999 3.8015 3.7951 -0.1674 46 3.4899 3.4830 -0.1999 3.8014 3.7951 -0.1675 47 3.8014 3.7951 -0.1675 48 3.8014 3.7950 -0.1676 49 3.8014 3.7950 -0.1676 50 3.8014 3.7950 -0.1676 PHỤ LỤC III CHƯƠNG TRÌNH MƠ PHỎNG KẾT QUẢ TÍNH TỐN KHOẢNG VƯỢT TỚI HẠN clc; clear all; disp('CHUONG TRINH MO PHONG KET QUA TINH TOAN KHOANG VUOT TOI HAN') % Nhap thong so dau vao n = 0;% So khoang vuot Z1 = 0; Z2 = 0; Z3 = 0; Z4 = 0; Z5 = 0; Z6 = 0; Rc1 = 5;% Dien tro cot ung voi Rc = Rc2 = 10; Rc3 = 15; Rc4 = 20; Rc5 = 25; Rc6 = 30; Rcs = input('Nhap dien tro day chong set Rcs='); fprintf(' n Rc=5 Rc=10 Rc=15 for I = 1:50; n = n+1; % Tinh gia tri b cong thuc 3.1 b1 = Rcs+2*Rc1; b2 = Rcs+2*Rc2; b3 = Rcs+2*Rc3; b4 = Rcs+2*Rc4; b5 = Rcs+2*Rc5; b6 = Rcs+2*Rc6; Rc=20 Rc=25 Rc=30 \n') % Tinh gia tri alpha1 cong thuc 3.1 a11= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc1))/2; a12= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc2))/2; a13= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc3))/2; a14= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc4))/2; a15= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc5))/2; a16= (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc6))/2; % Tinh gia tri alpha2 cong thuc 3.1 a21= -(Rcs+a11); a22= -(Rcs+a12); a23= -(Rcs+a13); a24= -(Rcs+a14); a25= -(Rcs+a15); a26= -(Rcs+a16); % Tinh tong tro Z tuong ung voi tung gia tri Rc Z1 = (a11*(b1-sqrt(Rcs*Rcs+4*Rcs*Rc1))^na21*(b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n)/((b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n-(b1sqrt(Rcs*Rcs+4*Rcs*Rc1))^n); Z2 = (a12*(b2-sqrt(Rcs*Rcs+4*Rcs*Rc2))^na22*(b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n)/((b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n-(b2sqrt(Rcs*Rcs+4*Rcs*Rc2))^n); Z3 = (a13*(b3-sqrt(Rcs*Rcs+4*Rcs*Rc3))^na23*(b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n)/((b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n-(b3sqrt(Rcs*Rcs+4*Rcs*Rc3))^n); Z4 = (a14*(b4-sqrt(Rcs*Rcs+4*Rcs*Rc4))^na24*(b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n)/((b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n-(b4sqrt(Rcs*Rcs+4*Rcs*Rc4))^n); Z5 = (a15*(b5-sqrt(Rcs*Rcs+4*Rcs*Rc5))^na25*(b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n)/((b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n-(b5sqrt(Rcs*Rcs+4*Rcs*Rc5))^n); Z6 = (a16*(b6-sqrt(Rcs*Rcs+4*Rcs*Rc6))^na26*(b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n)/((b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n-(b6sqrt(Rcs*Rcs+4*Rcs*Rc6))^n); disp([n, Z1, Z2, Z3, Z4, Z5, Z6]) end PHỤ LỤC IV CHƯƠNG TRÌNH MƠ PHỎNG KẾT QUẢ SO SÁNH SAI SỐ GIỮA CÔNG THỨC (3.6) VÀ (3.1) clc; clear all; disp('CHUONG TRINH SO SANH SAI SO GIUA (3.6) VA (3.1)') n = 0; Rc1 = 5;% Dien tro cot ung voi Rc = Rc2 = 10; Rc3 = 15; Rc4 = 20; Rc5 = 25; Rc6 = 30; Rcs = input(' Nhap dien tro day chong set Rcs='); fprintf(' Rc=5 Rc=10 Rc=15 fprintf(' n (3.1) (3.6) saiso (3.1) (3.6) saiso (3.1) (3.6) saiso (3.1) (3.6) saiso \n') for I = 1:50; n = n+1; % Tinh gia tri b cong thuc 3.1 b1 = Rcs+2*Rc1; b2 = Rcs+2*Rc2; b3 = Rcs+2*Rc3; b4 = Rcs+2*Rc4; b5 = Rcs+2*Rc5; b6 = Rcs+2*Rc6; % Tinh gia tri alpha1 cong thuc 3.1 a11 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc1))/2; a12 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc2))/2; a13 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc3))/2; a14 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc4))/2; a15 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc5))/2; a16 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc6))/2; Rc=20 (3.1) (3.6) saiso Rc=25 Rc=30 \n') (3.1) (3.6) saiso % Tinh gia tri alpha2 cong thuc 3.1 a21 = -(Rcs+a11); a22 = -(Rcs+a12); a23 = -(Rcs+a13); a24 = -(Rcs+a14); a25 = -(Rcs+a15); a26 = -(Rcs+a16); % Tinh tong tro Z theo cong thuc (3.1) tuong ung voi tung gia tri Rc Z1 = (a11*(b1-sqrt(Rcs*Rcs+4*Rcs*Rc1))^na21*(b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n)/((b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n-(b1sqrt(Rcs*Rcs+4*Rcs*Rc1))^n); Z2 = (a12*(b2-sqrt(Rcs*Rcs+4*Rcs*Rc2))^na22*(b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n)/((b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n-(b2sqrt(Rcs*Rcs+4*Rcs*Rc2))^n); Z3 = (a13*(b3-sqrt(Rcs*Rcs+4*Rcs*Rc3))^na23*(b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n)/((b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n-(b3sqrt(Rcs*Rcs+4*Rcs*Rc3))^n); Z4 = (a14*(b4-sqrt(Rcs*Rcs+4*Rcs*Rc4))^na24*(b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n)/((b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n-(b4sqrt(Rcs*Rcs+4*Rcs*Rc4))^n); Z5 = (a15*(b5-sqrt(Rcs*Rcs+4*Rcs*Rc5))^na25*(b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n)/((b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n-(b5sqrt(Rcs*Rcs+4*Rcs*Rc5))^n); Z6 = (a16*(b6-sqrt(Rcs*Rcs+4*Rcs*Rc6))^na26*(b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n)/((b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n-(b6sqrt(Rcs*Rcs+4*Rcs*Rc6))^n); % Tinh tong tro Z theo cong thuc (3.6)tuong ung voi tung gia tri Rc Z11 = Rcs*(1/2+sqrt(1/4+Rc1/Rcs)); Z12 = Rcs*(1/2+sqrt(1/4+Rc2/Rcs)); Z13 = Rcs*(1/2+sqrt(1/4+Rc3/Rcs)); Z14 = Rcs*(1/2+sqrt(1/4+Rc4/Rcs)); Z15 = Rcs*(1/2+sqrt(1/4+Rc5/Rcs)); Z16 = Rcs*(1/2+sqrt(1/4+Rc6/Rcs)); % Tinh sai so saiso1 = ((Z11-Z1)/Z1)*100; saiso2 = ((Z12-Z2)/Z2)*100; saiso3 = ((Z13-Z3)/Z3)*100; saiso4 = ((Z14-Z4)/Z4)*100; saiso5 = ((Z15-Z5)/Z5)*100; saiso6 = ((Z16-Z6)/Z6)*100; disp([n,Z1,Z11,saiso1,Z2,Z12,saiso2,Z3,Z13,saiso3,Z4,Z14,saiso4,Z5,Z15,saiso5,Z6,Z16 ,saiso6]) end PHỤ LỤC V CHƯƠNG TRÌNH MƠ PHỎNG KẾT QUẢ SO SÁNH SAI SỐ GIỮA CÔNG THỨC (3.8) VÀ (3.1) clc; clear all; disp('CHUONG TRINH SO SANH SAI SO GIUA (3.8) VA (3.1)') n = 0; Rc1 = 5;% Dien tro cot ung voi Rc = Rc2 = 10; Rc3 = 15; Rc4 = 20; Rc5 = 25; Rc6 = 30; Rcs = input(' Nhap dien tro day chong set Rcs='); fprintf(' Rc=5 Rc=10 Rc=15 fprintf(' n (3.1) (3.8) saiso (3.1) (3.8) saiso (3.1) (3.8) saiso (3.1) (3.8) saiso \n') for I = 1:50; n = n+1; % Tinh gia tri b cong thuc 3.1 b1 = Rcs+2*Rc1; b2 = Rcs+2*Rc2; b3 = Rcs+2*Rc3; b4 = Rcs+2*Rc4; b5 = Rcs+2*Rc5; b6 = Rcs+2*Rc6; % Tinh gia tri alpha1 cong thuc 3.1 a11 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc1))/2; a12 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc2))/2; a13 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc3))/2; a14 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc4))/2; Rc=20 (3.1) (3.8) saiso Rc=25 Rc=30 \n') (3.1) (3.8) saiso a15 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc5))/2; a16 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc6))/2; % Tinh gia tri alpha2 cong thuc 3.1 a21 = -(Rcs+a11); a22 = -(Rcs+a12); a23 = -(Rcs+a13); a24 = -(Rcs+a14); a25 = -(Rcs+a15); a26 = -(Rcs+a16); % Tinh tong tro Z theo cong thuc (3.1) tuong ung voi tung gia tri Rc Z1 = (a11*(b1-sqrt(Rcs*Rcs+4*Rcs*Rc1))^na21*(b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n)/((b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n-(b1sqrt(Rcs*Rcs+4*Rcs*Rc1))^n); Z2 = (a12*(b2-sqrt(Rcs*Rcs+4*Rcs*Rc2))^na22*(b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n)/((b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n-(b2sqrt(Rcs*Rcs+4*Rcs*Rc2))^n); Z3 = (a13*(b3-sqrt(Rcs*Rcs+4*Rcs*Rc3))^na23*(b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n)/((b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n-(b3sqrt(Rcs*Rcs+4*Rcs*Rc3))^n); Z4 = (a14*(b4-sqrt(Rcs*Rcs+4*Rcs*Rc4))^na24*(b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n)/((b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n-(b4sqrt(Rcs*Rcs+4*Rcs*Rc4))^n); Z5 = (a15*(b5-sqrt(Rcs*Rcs+4*Rcs*Rc5))^na25*(b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n)/((b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n-(b5sqrt(Rcs*Rcs+4*Rcs*Rc5))^n); Z6 = (a16*(b6-sqrt(Rcs*Rcs+4*Rcs*Rc6))^na26*(b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n)/((b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n-(b6sqrt(Rcs*Rcs+4*Rcs*Rc6))^n); % Tinh tong tro Z theo cong thuc (3.7)tuong ung voi tung gia tri Rc Z11 = sqrt(Rcs*Rc1)*coth((n*sqrt(Rcs/Rc1))); Z12 = sqrt(Rcs*Rc2)*coth((n*sqrt(Rcs/Rc2))); Z13 = sqrt(Rcs*Rc3)*coth((n*sqrt(Rcs/Rc3))); Z14 = sqrt(Rcs*Rc4)*coth((n*sqrt(Rcs/Rc4))); Z15 = sqrt(Rcs*Rc5)*coth((n*sqrt(Rcs/Rc5))); Z16 = sqrt(Rcs*Rc6)*coth((n*sqrt(Rcs/Rc6))); % Tinh sai so saiso1 = ((Z11-Z1)/Z1)*100; saiso2 = ((Z12-Z2)/Z2)*100; saiso3 = ((Z13-Z3)/Z3)*100; saiso4 = ((Z14-Z4)/Z4)*100; saiso5 = ((Z15-Z5)/Z5)*100; saiso6 = ((Z16-Z6)/Z6)*100; disp([n,Z1,Z11,saiso1,Z2,Z12,saiso2,Z3,Z13,saiso3,Z4,Z14,saiso4,Z5,Z15,saiso5,Z6,Z16 ,saiso6]) end PHỤ LỤC VI CHƯƠNG TRÌNH MƠ PHỎNG KẾT QUẢ SO SÁNH SAI SỐ GIỮA CÔNG THỨC (3.9) VÀ (3.1) clc; clear all; disp('CHUONG TRINH SO SANH SAI SO GIUA (3.9) VA (3.1)') n = 0; Rc1 = 5;% Dien tro cot ung voi Rc = Rc2 = 10; Rc3 = 15; Rc4 = 20; Rc5 = 25; Rc6 = 30; Rcs = input(' Nhap dien tro day chong set Rcs='); fprintf(' Rc=5 Rc=10 Rc=15 fprintf(' n (3.1) (3.9) saiso (3.1) (3.9) saiso (3.1) (3.9) saiso (3.1) (3.9) saiso \n') for I = 1:50; n = n+1; % Tinh gia tri b cong thuc 3.1 b1 = Rcs+2*Rc1; b2 = Rcs+2*Rc2; b3 = Rcs+2*Rc3; b4 = Rcs+2*Rc4; b5 = Rcs+2*Rc5; b6 = Rcs+2*Rc6; % Tinh gia tri alpha1 cong thuc 3.1 a11 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc1))/2; a12 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc2))/2; a13 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc3))/2; a14 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc4))/2; a15 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc5))/2; Rc=20 (3.1) (3.9) saiso Rc=25 Rc=30 \n') (3.1) (3.9) saiso a16 = (-Rcs+sqrt(Rcs*Rcs+4*Rcs*Rc6))/2; % Tinh gia tri alpha2 cong thuc 3.1 a21 = -(Rcs+a11); a22 = -(Rcs+a12); a23 = -(Rcs+a13); a24 = -(Rcs+a14); a25 = -(Rcs+a15); a26 = -(Rcs+a16); % Tinh tong tro Z theo cong thuc (3.1) tuong ung voi tung gia tri Rc Z1 = (a11*(b1-sqrt(Rcs*Rcs+4*Rcs*Rc1))^na21*(b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n)/((b1+sqrt(Rcs*Rcs+4*Rcs*Rc1))^n-(b1sqrt(Rcs*Rcs+4*Rcs*Rc1))^n); Z2 = (a12*(b2-sqrt(Rcs*Rcs+4*Rcs*Rc2))^na22*(b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n)/((b2+sqrt(Rcs*Rcs+4*Rcs*Rc2))^n-(b2sqrt(Rcs*Rcs+4*Rcs*Rc2))^n); Z3 = (a13*(b3-sqrt(Rcs*Rcs+4*Rcs*Rc3))^na23*(b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n)/((b3+sqrt(Rcs*Rcs+4*Rcs*Rc3))^n-(b3sqrt(Rcs*Rcs+4*Rcs*Rc3))^n); Z4 = (a14*(b4-sqrt(Rcs*Rcs+4*Rcs*Rc4))^na24*(b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n)/((b4+sqrt(Rcs*Rcs+4*Rcs*Rc4))^n-(b4sqrt(Rcs*Rcs+4*Rcs*Rc4))^n); Z5 = (a15*(b5-sqrt(Rcs*Rcs+4*Rcs*Rc5))^na25*(b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n)/((b5+sqrt(Rcs*Rcs+4*Rcs*Rc5))^n-(b5sqrt(Rcs*Rcs+4*Rcs*Rc5))^n); Z6 = (a16*(b6-sqrt(Rcs*Rcs+4*Rcs*Rc6))^na26*(b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n)/((b6+sqrt(Rcs*Rcs+4*Rcs*Rc6))^n-(b6sqrt(Rcs*Rcs+4*Rcs*Rc6))^n); % Tinh tong tro Z theo cong thuc (3.9)tuong ung voi tung gia tri Rc Z11 = (sqrt(Rcs/(4*Rc1))+1)*sqrt(Rcs*Rc1)*coth((n*sqrt(Rcs/Rc1))); Z12 = (sqrt(Rcs/(4*Rc2))+1)*sqrt(Rcs*Rc2)*coth((n*sqrt(Rcs/Rc2))); Z13 = (sqrt(Rcs/(4*Rc3))+1)*sqrt(Rcs*Rc3)*coth((n*sqrt(Rcs/Rc3))); Z14 = (sqrt(Rcs/(4*Rc4))+1)*sqrt(Rcs*Rc4)*coth((n*sqrt(Rcs/Rc4))); Z15 = (sqrt(Rcs/(4*Rc5))+1)*sqrt(Rcs*Rc5)*coth((n*sqrt(Rcs/Rc5))); Z16 = (sqrt(Rcs/(4*Rc6))+1)*sqrt(Rcs*Rc6)*coth((n*sqrt(Rcs/Rc6))); % Tinh sai so saiso1 = ((Z11-Z1)/Z1)*100; saiso2 = ((Z12-Z2)/Z2)*100; saiso3 = ((Z13-Z3)/Z3)*100; saiso4 = ((Z14-Z4)/Z4)*100; saiso5 = ((Z15-Z5)/Z5)*100; saiso6 = ((Z16-Z6)/Z6)*100; disp([n,Z1,Z11,saiso1,Z2,Z12,saiso2,Z3,Z13,saiso3,Z4,Z14,saiso4,Z5,Z15,saiso5,Z6,Z16 ,saiso6]) end