Ứng dụng mô hình đường dây truyền tải và phương pháp phân sai phần hữu hạn vào việc nghiên cứu quá độ trong hệ thống nối đất

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA HỒ ĐỨC HUẤN ỨNG DỤNG MƠ HÌNH ĐƯỜNG DÂY TRUYỀN TẢI VÀ PHƯƠNG PHÁP SAI PHÂN HỮU HẠN VÀO VIỆC NGHIÊN CỨU QUÁ ĐỘ TRONG HỆ THỐNG NỐI ĐẤT CHUYÊN NGÀNH : THIẾT BỊ, MẠNG VÀ NHÀ MÁY ĐIỆN Mà SỐ NGÀNH : 60 52 50 LUẬN VĂN THẠC SỸ TP HỒ CHÍ MINH, THÁNG 11 NĂM 2008 CƠNG TRÌNH HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: TS HOÀNG VIỆT Cán chấm nhận xét 1: Cán chấm nhận xét 2: Luận văn thạc sĩ bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày tháng năm 2008 Có thể tham khảo luận văn THƯ VIỆN TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA Xà HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH ĐỘC LẬP – TỰ DO – HẠNH PHÚC Tp, HCM, ngày tháng năm 2008 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : HỒ ĐỨC HUẤN Phái: Nam Ngày, tháng, năm sinh: 30/04/1976 Nơi sinh: Gia Lai Chuyên ngành: Thiết Bị, Mạng Nhà Máy Điện MSHV:01806482 I - TÊN ĐỀ TÀI: Ứng dụng mơ hình đường dây truyền tải phương pháp sai phân hữu hạn vào việc nghiên cứu độ hệ thống nối đất II - NHIỆM VỤ VÀ NỘI DUNG: Tổng quan : a) Về sét tình hình dơng sét Việt Nam b) Về nhân tố ảnh hưởng đến hệ thống nối đất trạm biến áp đường dây truyền tải tản dịng sét c) Về phương pháp tính toán tổng trở tản xung lưới nối đất Chọn phương pháp tính tốn tổng trở tản xung lưới nối đất có ý đến ảnh hưởng tượng ion hóa đất Lập trình tính tốn tổng trở tản xung lưới nối đất có số lưới khác nhau, dịng sét vào lưới vị trí khác Đánh giá kết quả, so sánh kiểm chứng với kết tính tốn tác giả khác, với phương pháp tính tốn khác Kết luận ưu nhược điểm phương pháp hướng phát triển đề tài III - NGÀY GIAO NHIỆM VỤ IV - NGÀY HOÀN THÀNH NHIỆM VỤ: 30/11/2008 V - CÁN BỘ HƯỚNG DẪN: TS HOÀNG VIỆT CÁN BỘ HƯỚNG DẪN CN BỘ MÔN QL CHUYÊN NGÀNH TS Hoàng Việt Nội dung đề cương luận văn thạc sĩ hội đồng chuyên ngành thông qua Ngày tháng năm 2008 TRƯỞNG PHÒNG ĐT – SĐH TRƯỞNG KHOA QL NGÀNH Để hồn thành chương trình cao học luận văn tốt nghiệp, nhận quan tâm, tạo điều kiện thuận lợi giúp đỡ tận tình Ban giám hiệu, phịng đào tạo sau đại học, khoa Điện – Điện tử thầy cô giáo trường Đại học Bách khoa thành phố Hồ Chí Minh Tơi xin chân thành cảm ơn quý thầy cô giảng dạy chuyên ngành truyền đạt cho nhiều kiến thức sâu rộng, mẻ lĩnh vực Thiết bị - Mạng Nhà máy điện Tơi xin bày tỏ lịng biết ơn sâu sắc đến thầy TS Hoàng Việt - Cán hướng dẫn khoa học dành nhiều thời gian, tận tâm trực tiếp dẫn chu đáo, truyền đạt nhiều kinh nghiệm quý báu để thực đề tài nghiên cứu khoa học hồn thành luận văn tốt nghiệp Tơi xin chân thành cảm ơn Ban giám đốc, phòng chức Điện lực Gia Lai nói chung, lãnh đạo đồng nghiệp phân xưởng phát điện Biển Hồ nói riêng tạo điều kiện, giúp đỡ, chia sẻ, động viên tơi suốt q trình học tập, thực đề tài, hồn thành luận văn tốt nghiệp Tơi xin gởi lời cảm ơn đến bạn lớp cao học khóa 2006 chuyên ngành Thiết bị - Mạng Nhà máy điện dành nhiều động viên, chia sẻ, giúp đỡ tơi suốt q trình học tập thực đề tài Xin chân thành ghi tâm tình cảm, cơng lao, cổ vũ động viên thân tình gia đình bè bạn, nguồn động lực lớn để tiếp sức cho phấn đấu vươn lên bước đường cơng tác Xin trân trọng cảm ơn! Hồ Đức Huấn ================================================================ TĨM TẮT LUẬN VĂN Nước ta chịu khí hậu nhiệt đới gió mùa, bị ảnh hưởng khí hậu lục địa lẫn khí hậu biển nên có lượng mưa bão kèm giông sét hàng năm lớn Số ngày giơng trung bình khoảng 100 ngày/năm số giơng trung bình 250 giờ/năm, cường độ sét lớn đo Việt nam 98kA (8/20ms) Hàng năm dông sét thường tập trung từ tháng 05 - tháng 11 Do đặc điểm khí hậu nêu nên ngành điện Việt Nam gặp nhiều khó khăn việc giảm thiểu thiệt hại sét gây Các thiết bị thu sét không phát huy hiệu bị phản tác dụng khơng có hệ thống nối đất tốt để tản dòng sét vào đất giữ cho điện phần tử nối đất không cao để hạn chế phóng điện ngược từ phần tử đến phận mang điện thiết bị điện khác Khi tản dịng điện sét, tượng ion hóa đất (hay tượng phóng điện đất) ảnh hưởng điện cảm hệ thống nối đất tượng vật lý chủ yếu quan tâm nghiên cứu Các cơng trình nghiên cứu giới lĩnh vực tượng ion hóa đất dừng lại khảo sát mơ hình vật lý Điểm luận văn sử dụng mơ hình đường dây truyền tải phương pháp sai phân hữu hạn, cách rời rạc hóa thông số (R,L,C,G) phân đoạn dẫn để khảo sát tượng ion hóa hệ thống nối đất sở sử dụng mối quan hệ cường độ điện trường ngưỡng gây nên tượng ion hóa đất thay đổi bán kính biểu kiến phân đoạn dẫn xảy tượng ion hóa lưới nối đất Kết tính tốn chương trình tính tốn cho thấy cần thiết phải tính đến tượng ion hóa đất thiết kế hệ thống nối đất, đặt biệt đất có điện trở suất cao kích thước lưới nối đất khoảng 100x100m2 trở xuống Điều với điều kiện nêu trên, tượng ion hóa xảy mãnh liệt làm tăng kích thước biểu kiến dẫn tạo nên lưới đặc biệt có điểm nối với ================================================================ ================================================================ điểm sét vào lưới, gây nên giảm điện trở tản xung đáng kể hệ thống nối đất Việc nghiên cứu tượng ion hóa hệ thống nối đất góp phần làm rõ chất, chế vật lý tượng khả giảm chi phí xây dựng hệ thống nối đất tính tốn ảnh hưởng tượng lên điện trở tản xung hệ thống Hướng phát triển đề tài : Trong q trình nghiên cứu đề tài, khơng có điều kiện nghiên cứu thực nghiệm nên đề tài sử dụng kết nghiên cứu dẫn nối đất để áp dụng vào tính tốn q độ lưới nối đất Điều làm cho kết tính tốn chương trình bị hạn chế, chưa phản ảnh xác giá trị vật lý thực lưới nối đất Để hoàn chỉnh đề tài cần có thêm nghiên cứu thực nghiệm bổ sung yếu tố ảnh hưởng trình độ vật lý xảy lưới vào thông số (R,L,C,G) phân đoạn dẫn theo kích thước, hình dạng hình học lưới ================================================================ ================================================================ MỤC LỤC Chương : Tổng quan sét, mơ hình hệ thống nối đất phương pháp số 1.1 Sơ lược hình thành sét , tham số chủ yếu cường độ hoạt động sét 1.1.1 Sơ lược hình thành sét 1.1.2 Các tham số chủ yếu sét 1.1.3 Mật độ sét đánh trung bình: 1.2.Tình hình dơng sét Việt nam 1.2.1 Đặc điểm tình hình dơng sét Việt Nam 1.2.2 Phân hóa mùa dơng lãnh thổ Việt Nam 1.2.3 Khả dự báo sét Việt Nam 10 1.3 Sự cần thiết việc nghiên cứu, tính tốn tổng trở tản xung hệ thống nối đất 12 1.4 Các yếu tố ảnh hưởng đến hệ thống nối đất tản dòng sét 13 1.4.1 Điện trở tản xung dạng nối đất tập trung 13 1.4.2 Điện trở tản xung nối đất kéo dài 13 1.4.2.1 Khi bỏ qua tượng ion hóa đất 13 1.4.2.2 Khi có q trình phóng điện đất 14 1.5 Tổng quan mơ hình hệ thống nối đất phương pháp số 14 1.5.1 Sơ lược hướng tiếp cận mơ hình tính tốn hệ thống nối đất thu sét 14 1.5.2 Cách tiếp cận mạch điện (Circuit approach) 16 1.5.3 Cách tiếp cận trường điện từ :( Electromagnetic field approach) 19 1.5.4 Cách tiếp cận hỗn hợp : (Hybrid approach) 20 1.5 Cách tiếp cận theo mơ hình đường dây truyền tải : (Transmission line approach) 20 1.6 Giới thiệu phương pháp số phổ biến 21 ================================================================ 1.6.1 Phương pháp phần tử hữu hạn ( FEM-Finite Element Method) 21 1.6.2 Phương pháp sai phân hữu hạn : (FDM- Finite Differential Method) 22 1.7 Hướng nghiên cứu đề tài 22 Chương : Hiện tượng ion hóa đất mơ hình hệ thống nối đất 24 2.1.Giới thiệu 24 2.2 Vấn đề điện trở suất dư (residual resistivity) vùng đất bị ion hóa 24 2.2.1 Các kết thí nghiệm liên quan đến điện trở suất dư (residual resistivity) vùng đất bị ion hóa 24 2.2.2 Kết luận điện trở suất dư vùng đất bị ion hóa 33 2.3 Vấn đề tượng ion hóa đất điện áp xung sét 34 2.4 Quan điểm hình dạng vùng ion hóa nối đất 37 2.5 Nhận xét chung tượng ion hóa đất sét 40 Chương Ứng dụng phương pháp sai phân hữu hạn vào mơ hình đường dây truyền tải nghiên cứu độ hệ thống nối đất 41 3.1 Phương pháp sai phân hữu hạn 41 3.2 Mơ hình đường dây truyền tải hệ thống nối đất 42 3.3 Áp dụng phương pháp sai phân vào mơ hình đường dây truyền tải 44 3.4 Cách tính điện áp, dịng điện dẫn chôn đất xét đến tượng ion hóa đất 46 3.5 Cách xác định dòng điện tản 47 3.6 Cách tính tốn điện trở tản xung hệ thống nối đất 47 3.7 Quan hệ điểm nút lưới 48 3.8 Thuật giải tổng quát cho trường lưới nối đất có xét đến q trình ion hóa 49 ================================================================ 3.9 Thông số đầu vào đầu chương trình tính tốn q độ 49 Chương : Kết tính tốn q độ lưới nối đất 51 4.1 Kết thu nối đất dạng 51 4.1.1 Trường hợp : So sánh với kết [29] 51 4.1.2 Trường hợp 53 4.1.3 Trường hợp 55 4.2 Tính tốn q độ chưa xét đến tượng ion hóa đất 57 4.2.1 Kết thu nối đất lưới 1x1 58 4.2.2 Kết thu nối đất lưới 2x2 60 4.2.3 Kết thu nối đất lưới 4x4 65 4.2.4 Kết thu nối đất lưới 6x6 70 4.3.Kết tính tốn q độ xét thay đổi điện trở suất đất 76 4.3.1 Lưới nối đất 2x2 76 4.3.2 Lưới nối đất 4x4 77 4.3.3 Lưới nối đất 6x6 78 4.4 Kết tính tốn q độ có xét đến tượng ion hóa đất 79 4.4.1 Trường hợp 79 4.4.2 Trường hợp lưới 4x4 82 4.4.2 Trường hợp lưới 6x6 84 4.5.Kết luận& đánh giá chung 86 Phụ lục 87 Chương : Tổng quan sét, mơ hình hệ thống nối đất phương pháp số ================================================================ Chương : TỔNG QUAN VỀ SÉT, CÁC MƠ HÌNH HỆ THỐNG NỐI ĐẤT VÀ PHƯƠNG PHÁP SỐ 1.1 Sơ lược hình thành sét , tham số chủ yếu cường độ hoạt động sét [1] 1.1.1 Sơ lược hình thành sét : Sét thực chất dạng phóng điện tia lửa khơng khí với khoảng cách phóng điện lớn Chiều dài trung bình khe sét khoảng ÷ km, phần lớn chiều dài phát triển đám mây dơng Q trình phóng điện sét tương tự q trình phóng điện tia lửa điện trường không đồng với khoảng cách phóng điện lớn Chính tương tự cho phép mơ sét phịng thí nghiệm để nghiên cứu qui luật nghiên cứu biện pháp bảo vệ chống sét Q trình phóng điện sét gồm giai đoạn chủ yếu : Hình 1.1:Các giai đoạn phóng điện biến thiên dịng điện sét theo thời gian (a) Giai đoạn phóng điện tiên đạo (b) Tia tiên đạo đến gần mặt đất hình thành khu vực ion hóa mãnh liệt (c) Giai đoạn phóng điện ngược hay phóng điện chủ yếu (d) Phóng điện chủ yếu kết thúc, dịng sét đạt giá trị cực đại Giai đoạn phóng điện tiên đạo đợt : Trang 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m of 17 I(41,1,n+1)=A(41,1,n)*I(41,1,n)+B(41,1,n)*(V(17,n)-U(41,2,n)); I(41,Sm,n+1)=A(41,Sm,n)*I(41,Sm,n)+B(41,Sm,n)*(U(41,Sm,n)-V(26,n)); I(42,1,n+1)=A(42,1,n)*I(42,1,n)+B(42,1,n)*(V(26,n)-U(42,2,n)); I(42,Sm,n+1)=A(42,Sm,n)*I(42,Sm,n)+B(42,Sm,n)*(U(42,Sm,n)-V(27,n)); I(43,1,n+1)=A(43,1,n)*I(43,1,n)+B(43,1,n)*(V(18,n)-U(43,2,n)); I(43,Sm,n+1)=A(43,Sm,n)*I(43,Sm,n)+B(43,Sm,n)*(U(43,Sm,n)-V(27,n)); I(44,1,n+1)=A(44,1,n)*I(44,1,n)+B(44,1,n)*(V(27,n)-U(44,2,n)); I(44,Sm,n+1)=A(44,Sm,n)*I(44,Sm,n)+B(44,Sm,n)*(U(44,Sm,n)-V(28,n)); I(45,1,n+1)=A(45,1,n)*I(45,1,n)+B(45,1,n)*(V(19,n)-U(45,2,n)); I(45,Sm,n+1)=A(45,1,n)*I(45,Sm,n)+B(45,1,n)*(U(45,Sm,n)-V(28,n)); I(46,1,n+1)=A(46,1,n)*I(46,1,n)+B(46,1,n)*(V(28,n)-U(46,2,n)); I(46,Sm,n+1)=A(46,Sm,n)*I(46,Sm,n)+B(46,Sm,n)*(U(46,Sm,n)-V(29,n)); I(47,1,n+1)=A(47,1,n)*I(47,1,n)+B(47,1,n)*(V(20,n)-U(47,2,n)); I(47,Sm,n+1)=A(47,Sm,n)*I(47,Sm,n)+B(47,Sm,n)*(U(47,Sm,n)-V(29,n)); I(48,1,n+1)=A(48,1,n)*I(48,1,n)+B(48,1,n)*(V(29,n)-U(48,2,n)); I(48,Sm,n+1)=A(48,Sm,n)*I(48,Sm,n)+B(48,Sm,n)*(U(48,Sm,n)-V(30,n)); I(49,1,n+1)=A(49,1,n)*I(49,1,n)+B(49,1,n)*(V(21,n)-U(49,2,n)); I(49,Sm,n+1)=A(49,Sm,n)*I(49,Sm,n)+B(49,Sm,n)*(U(49,Sm,n)-V(30,n)); I(50,1,n+1)=A(50,1,n)*I(50,1,n)+B(50,1,n)*(V(30,n)-U(50,2,n)); I(50,Sm,n+1)=A(50,Sm,n)*I(50,Sm,n)+B(50,Sm,n)*(U(50,Sm,n)-V(31,n)); I(51,1,n+1)=A(51,1,n)*I(51,1,n)+B(51,1,n)*(V(32,n)-U(51,2,n)); I(51,Sm,n+1)=A(51,Sm,n)*I(51,Sm,n)+B(51,Sm,n)*(U(51,Sm,n)-V(31,n)); I(52,1,n+1)=A(52,1,n)*I(52,1,n)+B(52,1,n)*(V(21,n)-U(52,2,n)); I(52,Sm,n+1)=A(52,Sm,n)*I(52,Sm,n)+B(52,Sm,n)*(U(52,Sm,n)-V(32,n)); I(53,1,n+1)=A(53,1,n)*I(53,1,n)+B(53,1,n)*(V(33,n)-U(53,2,n)); I(53,Sm,n+1)=A(53,Sm,n)*I(53,Sm,n)+B(53,Sm,n)*(U(53,Sm,n)-V(32,n)); I(54,1,n+1)=A(54,1,n)*I(54,1,n)+B(54,1,n)*(V(22,n)-U(54,2,n)); I(54,Sm,n+1)=A(54,Sm,n)*I(54,Sm,n)+B(54,Sm,n)*(U(54,Sm,n)-V(33,n)); I(55,1,n+1)=A(55,1,n)*I(55,1,n)+B(55,1,n)*(V(34,n)-U(55,2,n)); I(55,Sm,n+1)=A(55,Sm,n)*I(55,Sm,n)+B(55,Sm,n)*(U(55,Sm,n)-V(33,n)); I(56,1,n+1)=A(56,1,n)*I(56,1,n)+B(56,1,n)*(V(23,n)-U(56,2,n)); I(56,Sm,n+1)=A(56,Sm,n)*I(56,Sm,n)+B(56,Sm,n)*(U(56,Sm,n)-V(34,n)); I(57,1,n+1)=A(57,1,n)*I(57,1,n)+B(57,1,n)*(V(35,n)-U(57,2,n)); I(57,Sm,n+1)=A(57,Sm,n)*I(57,Sm,n)+B(57,Sm,n)*(U(57,Sm,n)-V(34,n)); I(58,1,n+1)=A(58,1,n)*I(58,1,n)+B(58,1,n)*(V(24,n)-U(58,2,n)); I(58,Sm,n+1)=A(58,Sm,n)*I(58,Sm,n)+B(58,Sm,n)*(U(58,Sm,n)-V(35,n)); 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m of 17 I(59,1,n+1)=A(59,1,n)*I(59,1,n)+B(59,1,n)*(V(36,n)-U(59,2,n)); I(59,Sm,n+1)=A(59,Sm,n)*I(59,Sm,n)+B(59,Sm,n)*(U(59,Sm,n)-V(35,n)); I(60,1,n+1)=A(60,1,n)*I(60,1,n)+B(60,1,n)*(V(25,n)-U(60,2,n)); I(60,Sm,n+1)=A(60,Sm,n)*I(60,Sm,n)+B(60,Sm,n)*(U(60,Sm,n)-V(36,n)); I(61,1,n+1)=A(61,1,n)*I(61,1,n)+B(61,1,n)*(V(26,n)-U(61,2,n)); I(61,Sm,n+1)=A(61,Sm,n)*I(61,Sm,n)+B(61,Sm,n)*(U(61,Sm,n)-V(37,n)); I(62,1,n+1)=A(62,1,n)*I(62,1,n)+B(62,1,n)*(V(37,n)-U(62,2,n)); I(62,Sm,n+1)=A(62,Sm,n)*I(62,Sm,n)+B(62,Sm,n)*(U(62,Sm,n)-V(38,n)); I(63,1,n+1)=A(63,1,n)*I(63,1,n)+B(63,1,n)*(V(27,n)-U(63,2,n)); I(63,Sm,n+1)=A(63,Sm,n)*I(63,Sm,n)+B(63,Sm,n)*(U(63,Sm,n)-V(38,n)); I(64,1,n+1)=A(64,1,n)*I(64,1,n)+B(64,1,n)*(V(38,n)-U(64,2,n)); I(64,Sm,n+1)=A(64,Sm,n)*I(64,Sm,n)+B(64,Sm,n)*(U(64,Sm,n)-V(39,n)); I(65,1,n+1)=A(65,1,n)*I(65,1,n)+B(65,1,n)*(V(28,n)-U(65,2,n)); I(65,Sm,n+1)=A(65,Sm,n)*I(65,Sm,n)+B(65,Sm,n)*(U(65,Sm,n)-V(39,n)); I(66,1,n+1)=A(66,1,n)*I(66,1,n)+B(66,1,n)*(V(39,n)-U(66,2,n)); I(66,Sm,n+1)=A(66,Sm,n)*I(66,Sm,n)+B(66,Sm,n)*(U(66,Sm,n)-V(40,n)); I(67,1,n+1)=A(67,1,n)*I(67,1,n)+B(67,1,n)*(V(29,n)-U(67,2,n)); I(67,Sm,n+1)=A(67,Sm,n)*I(67,Sm,n)+B(67,Sm,n)*(U(67,Sm,n)-V(40,n)); I(68,1,n+1)=A(68,1,n)*I(68,1,n)+B(68,1,n)*(V(40,n)-U(68,2,n)); I(68,Sm,n+1)=A(68,Sm,n)*I(68,Sm,n)+B(68,Sm,n)*(U(68,Sm,n)-V(41,n)); I(69,1,n+1)=A(69,1,n)*I(69,1,n)+B(69,1,n)*(V(30,n)-U(69,2,n)); I(69,Sm,n+1)=A(69,Sm,n)*I(69,Sm,n)+B(69,Sm,n)*(U(69,Sm,n)-V(41,n)); I(70,1,n+1)=A(70,1,n)*I(70,1,n)+B(70,1,n)*(V(41,n)-U(70,2,n)); I(70,Sm,n+1)=A(70,Sm,n)*I(70,Sm,n)+B(70,Sm,n)*(U(70,Sm,n)-V(42,n)); I(71,1,n+1)=A(71,1,n)*I(71,1,n)+B(71,1,n)*(V(31,n)-U(71,2,n)); I(71,Sm,n+1)=A(71,Sm,n)*I(71,Sm,n)+B(71,Sm,n)*(U(71,Sm,n)-V(42,n)); I(72,1,n+1)=A(72,1,n)*I(72,1,n)+B(72,1,n)*(V(42,n)-U(72,2,n)); I(72,Sm,n+1)=A(72,Sm,n)*I(72,Sm,n)+B(72,Sm,n)*(U(72,Sm,n)-V(43,n)); I(73,1,n+1)=A(73,1,n)*I(73,1,n)+B(73,1,n)*(V(44,n)-U(73,2,n)); I(73,Sm,n+1)=A(73,Sm,n)*I(73,Sm,n)+B(73,Sm,n)*(U(73,Sm,n)-V(43,n)); I(74,1,n+1)=A(74,1,n)*I(74,1,n)+B(74,1,n)*(V(31,n)-U(74,2,n)); I(74,Sm,n+1)=A(74,Sm,n)*I(74,Sm,n)+B(74,Sm,n)*(U(74,Sm,n)-V(44,n)); I(75,1,n+1)=A(75,1,n)*I(75,1,n)+B(75,1,n)*(V(45,n)-U(75,2,n)); I(75,Sm,n+1)=A(75,Sm,n)*I(75,Sm,n)+B(75,Sm,n)*(U(75,Sm,n)-V(44,n)); I(76,1,n+1)=A(76,1,n)*I(76,1,n)+B(76,1,n)*(V(32,n)-U(76,2,n)); I(76,Sm,n+1)=A(76,Sm,n)*I(76,Sm,n)+B(76,Sm,n)*(U(76,Sm,n)-V(45,n)); I(77,1,n+1)=A(77,1,n)*I(77,1,n)+B(77,1,n)*(V(46,n)-U(77,2,n)); 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m of 17 I(77,Sm,n+1)=A(77,Sm,n)*I(77,Sm,n)+B(77,Sm,n)*(U(77,Sm,n)-V(45,n)); I(78,1,n+1)=A(78,1,n)*I(78,1,n)+B(78,1,n)*(V(33,n)-U(78,2,n)); I(78,Sm,n+1)=A(78,Sm,n)*I(78,Sm,n)+B(78,Sm,n)*(U(78,Sm,n)-V(46,n)); I(79,1,n+1)=A(79,1,n)*I(79,1,n)+B(79,1,n)*(V(47,n)-U(79,2,n)); I(79,Sm,n+1)=A(79,Sm,n)*I(79,Sm,n)+B(79,Sm,n)*(U(79,Sm,n)-V(46,n)); I(80,1,n+1)=A(80,1,n)*I(80,1,n)+B(80,1,n)*(V(34,n)-U(80,2,n)); I(80,Sm,n+1)=A(80,Sm,n)*I(80,Sm,n)+B(80,Sm,n)*(U(80,Sm,n)-V(47,n)); I(81,1,n+1)=A(81,1,n)*I(81,1,n)+B(81,1,n)*(V(48,n)-U(81,2,n)); I(81,Sm,n+1)=A(81,Sm,n)*I(81,Sm,n)+B(81,Sm,n)*(U(81,Sm,n)-V(47,n)); I(82,1,n+1)=A(82,1,n)*I(82,1,n)+B(82,1,n)*(V(35,n)-U(82,2,n)); I(82,Sm,n+1)=A(82,Sm,n)*I(82,Sm,n)+B(82,Sm,n)*(U(82,Sm,n)-V(48,n)); I(83,1,n+1)=A(83,1,n)*I(83,1,n)+B(83,1,n)*(V(49,n)-U(83,2,n)); I(83,Sm,n+1)=A(83,Sm,n)*I(83,Sm,n)+B(83,Sm,n)*(U(83,Sm,n)-V(48,n)); I(84,1,n+1)=A(84,1,n)*I(84,1,n)+B(84,1,n)*(V(36,n)-U(84,2,n)); I(84,Sm,n+1)=A(84,Sm,n)*I(84,Sm,n)+B(84,Sm,n)*(U(84,Sm,n)-V(49,n)); %************************************************************* for nh=1:Snh for m=2:Sm U(nh,m,n+1)=CC(nh,m,n)*U(nh,m,n)+D(nh,m,n)*(I(nh,m-1,n+1)-I(nh,m, n+1)); end end %******************************* V(1,n+1)=CCn(1,n)*V(1,n)+ 0.5*Dn(1,n)*(Is(1,n)-I(1,1,n+1)-I(4,1,n+1)-I (1,1,n)-I(4,1,n)); V(2,n+1)=CCn(2,n)*V(2,n)+ 0.5*Dn(2,n)*(Is(2,n)+I(1,Sm,n+1)-I(5,1,n+1)-I (2,1,n+1)+I(1,Sm,n)-I(5,1,n)-I(2,1,n)); V(3,n+1)=CCn(3,n)*V(3,n)+ 0.5*Dn(3,n)*(Is(3,n)+I(2,Sm,n+1)+I(3,Sm,n+1)I(7,1,n+1)-I(10,1,n+1)+I(2,Sm,n)+I(3,Sm,n)-I(7,1,n)-I(10,1,n)); V(4,n+1)=CCn(4,n)*V(4,n)+ 0.5*Dn(4,n)*(Is(4,n)+I(4,Sm,n+1)-I(3,1,n+1)-I (12,1,n+1)+I(4,Sm,n)-I(3,1,n)-I(12,1,n)); V(5,n+1)=CCn(5,n)*V(5,n)+ 0.5*Dn(5,n)*(Is(5,n)+I(5,Sm,n+1)-I(6,1,n+1)I(13,1,n+1)+I(5,Sm,n)-I(6,1,n)- I(13,1,n)); V(6,n+1)=CCn(6,n)*V(6,n)+ 0.5*Dn(6,n)*(Is(6,n)+I(6,Sm,n+1)+I(7,Sm,n+1)I(8,1,n+1) -I(15,1,n+1)+I(6,Sm,n)+I(7,Sm,n)-I(8,1,n) -I(15,1,n)); V(7,n+1)=CCn(7,n)*V(7,n)+ 0.5*Dn(7,n)*(Is(7,n)+I(8,Sm,n+1)+I(9,Sm,n+1)I(17,1,n+1) -I(20,1,n+1)+I(8,Sm,n)+I(9,Sm,n)- I(17,1,n) -I(20,1,n)); V(8,n+1)=CCn(8,n)*V(8,n)+ 0.5*Dn(8,n)*(Is(8,n)+I(10,Sm,n+1)+I(11,Sm, n+1)-I(9,1,n+1)- I(22,1,n+1)+I(10,Sm,n)+I(11,Sm,n)-I(9,1,n)- I(22,1,n)); V(9,n+1)=CCn(9,n)*V(9,n)+ 0.5*Dn(9,n)*(Is(9,n)+I(12,Sm,n+1)-I(11,1,n+1) - I(24,1,n+1)+I(12,Sm,n)-I(11,1,n)- I(24,1,n)); V(10,n+1)=CCn(10,n)*V(10,n)+ 0.5*Dn(10,n)*(Is(10,n)+ I(13,Sm,n+1)-I (25,1,n+1)-I(14,1,n+1)+ I(13,Sm,n)-I(25,1,n)-I(14,1,n)); V(11,n+1)=CCn(11,n)*V(11,n)+ 0.5*Dn(11,n)*(Is(11,n)+I(14,Sm,n+1) +I(15, Sm,n+1)-I(27,1,n+1)-I(16,1,n+1)+I(14,Sm,n) +I(15,Sm,n)-I(27,1,n)-I(16,1,n)); 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m 10 of 17 V(12,n+1)=CCn(12,n)*V(12,n)+ 0.5*Dn(12,n)*(Is(12,n)+I(16,Sm,n+1)+I(17, Sm,n+1)-I(29,1,n+1)-I(18,1,n+1)+I(16,Sm,n)+I(17,Sm,n)-I(29,1,n)-I(18,1,n)); V(13,n+1)=CCn(13,n)*V(13,n)+ 0.5*Dn(13,n)*(Is(13,n)+I(18,Sm,n+1)+I(19, Sm,n+1)-I(31,1,n+1)-I(34,1,n+1)+I(18,Sm,n)+I(19,Sm,n)-I(31,1,n)-I(34,1,n)); V(14,n+1)=CCn(14,n)*V(14,n)+ 0.5*Dn(14,n)*(Is(14,n)+I(20,Sm,n+1) +I(21, Sm,n+1)-I(19,1,n+1)-I(36,1,n+1)+I(20,Sm,n) +I(21,Sm,n)-I(19,1,n)-I(36,1,n)); V(15,n+1)=CCn(15,n)*V(15,n)+ 0.5*Dn(15,n)*(Is(15,n)+I(22,Sm,n+1)+I(23, Sm,n+1)-I(21,1,n+1)+I(38,1,n+1)+I(22,Sm,n)+I(23,Sm,n)-I(21,1,n)+I(38,1,n)); V(16,n+1)=CCn(16,n)*V(16,n)+ 0.5*Dn(16,n)*(Is(16,n)+I(24,Sm,n+1)-I (23,1,n+1)-I(40,1,n+1)+I(24,Sm,n)-I(23,1,n)-I(40,1,n)); V(17,n+1)=CCn(17,n)*V(17,n)+ 0.5*Dn(17,n)*(Is(17,n)+I(25,Sm,n+1)-I (41,1,n+1)-I(26,1,n+1)+I(25,Sm,n)-I(41,1,n)-I(26,1,n)); V(18,n+1)=CCn(18,n)*V(18,n)+ 0.5*Dn(18,n)*(Is(18,n)+I(26,Sm,n+1)+I(27, Sm,n+1)-I(43,1,n+1)- I(28,1,n+1)+I(26,Sm,n)+I(27,Sm,n)-I(43,1,n)- I(28,1,n)); V(19,n+1)=CCn(19,n)*V(19,n)+ 0.5*Dn(19,n)*(Is(19,n)+ I(28,Sm,n+1)+I(29, Sm,n+1)-I(30,1,n+1)- I(45,1,n+1)+ I(28,Sm,n)+I(29,Sm,n)-I(30,1,n)- I(45,1,n)); V(20,n+1)=CCn(20,n)*V(20,n)+ 0.5*Dn(20,n)*(Is(20,n)+ I(30,Sm,n+1)+I(31, Sm,n+1)-I(47,1,n+1)- I(32,1,n+1)+ I(30,Sm,n)+I(31,Sm,n)-I(47,1,n)- I(32,1,n)); V(21,n+1)=CCn(21,n)*V(21,n)+ 0.5*Dn(21,n)*(Is(21,n)+ I(32,Sm,n+1)+I(33, Sm,n+1)-I(49,1,n+1)- I(52,1,n+1)+ I(32,Sm,n)+I(33,Sm,n)-I(49,1,n)- I(52,1,n)); V(22,n+1)=CCn(22,n)*V(22,n)+ 0.5*Dn(22,n)*(Is(22,n) + I(34,Sm,n+1)+I (35,Sm,n+1)-I(33,1,n+1)- I(54,1,n+1)+ I(34,Sm,n)+I(35,Sm,n)-I(33,1,n)- I(54,1,n)); V(23,n+1)=CCn(23,n)*V(23,n)+ 0.5*Dn(23,n)*(Is(23,n)+I(36,Sm,n+1)+I(37, Sm,n+1)-I(35,1,n+1)-I(56,1,n+1)+I(36,Sm,n)+I(37,Sm,n)-I(35,1,n)-I(56,1,n)); V(24,n+1)=CCn(24,n)*V(24,n)+ 0.5*Dn(24,n)*(Is(24,n) + I(38,Sm,n+1)+I (39,Sm,n+1)-I(37,1,n+1)- I(58,1,n+1)+ I(38,Sm,n)+I(39,Sm,n)-I(37,1,n)- I(58,1,n)); V(25,n+1)=CCn(25,n)*V(25,n)+ 0.5*Dn(25,n)*(Is(25,n)+I(40,Sm,n+1)-I (39,1,n+1) -I(60,1,n+1)+I(40,Sm,n)-I(39,1,n) -I(60,1,n)); V(26,n+1)=CCn(26,n)*V(26,n)+ 0.5*Dn(26,n)*(Is(26,n)+I(41,Sm,n+1)-I(61, Sm,n+1)-I(42,1,n+1)+I(41,Sm,n)-I(61,Sm,n)-I(42,1,n)); V(27,n+1)=CCn(27,n)*V(27,n)+ 0.5*Dn(27,n)*(Is(27,n) + I(42,Sm,n+1)+I (43,Sm,n+1)-I(63,1,n+1)- I(44,1,n+1)+ I(42,Sm,n)+I(43,Sm,n)-I(63,1,n)- I(44,1,n)); V(28,n+1)=CCn(28,n)*V(28,n)+ 0.5*Dn(28,n)*(Is(28,n) + I(44,Sm,n+1)+I (45,Sm,n+1)-I(65,1,n+1)- I(46,1,n+1)+ I(44,Sm,n)+I(45,Sm,n)-I(65,1,n)- I(46,1,n)); V(29,n+1)=CCn(29,n)*V(29,n)+ 0.5*Dn(29,n)*(Is(29,n) + I(46,Sm,n+1)+I (47,Sm,n+1)-I(67,1,n+1)- I(48,1,n+1)+ I(46,Sm,n)+I(47,Sm,n)-I(67,1,n)- I(48,1,n)); V(30,n+1)=CCn(30,n)*V(30,n)+ 0.5*Dn(30,n)*(Is(30,n) + I(48,Sm,n+1)+I (49,Sm,n+1)-I(69,1,n+1)- I(50,1,n+1)+ I(48,Sm,n)+I(49,Sm,n)-I(69,1,n)- I(50,1,n)); V(31,n+1)=CCn(31,n)*V(31,n)+ 0.5*Dn(31,n)*(Is(31,n) + I(50,Sm,n+1)+I (51,Sm,n+1)-I(71,1,n+1)- I(74,1,n+1)+ I(50,Sm,n)+I(51,Sm,n)-I(71,1,n)- I(74,1,n)); V(32,n+1)=CCn(32,n)*V(32,n)+ 0.5*Dn(32,n)*(Is(32,n) + I(52,Sm,n+1)+I (53,Sm,n+1)-I(51,1,n+1)- I(76,1,n+1)+ I(52,Sm,n)+I(53,Sm,n)-I(51,1,n)- I(76,1,n)); V(33,n+1)=CCn(33,n)*V(33,n)+ 0.5*Dn(33,n)*(Is(33,n) + I(54,Sm,n+1)+I (55,Sm,n+1)-I(53,1,n+1)-I(78,1,n+1)+ I(54,Sm,n)+I(55,Sm,n)-I(53,1,n)-I(78,1,n)); V(34,n+1)=CCn(34,n)*V(34,n)+ 0.5*Dn(34,n)*(Is(34,n) + I(56,Sm,n+1)+I (57,Sm,n+1)-I(55,1,n+1)- I(80,1,n+1)+ I(56,Sm,n)+I(57,Sm,n)-I(55,1,n)- I(80,1,n)); V(35,n+1)=CCn(35,n)*V(35,n)+ 0.5*Dn(35,n)*(Is(35,n) + I(58,Sm,n+1)+I (59,Sm,n+1)-I(57,1,n+1)- I(82,1,n+1)+ I(58,Sm,n)+I(59,Sm,n)-I(57,1,n)- I(82,1,n)); V(36,n+1)=CCn(36,n)*V(36,n)+ 0.5*Dn(36,n)*(Is(36,n)+I(60,Sm,n+1)-I (59,1,n+1)-I(84,1,n+1)+I(60,Sm,n)-I(59,1,n)-I(84,1,n)); V(37,n+1)=CCn(37,n)*V(37,n)+ 0.5*Dn(37,n)*(Is(37,n)+I(61,Sm,n+1) -I (62,1,n+1)+I(61,Sm,n) -I(62,1,n)); V(38,n+1)=CCn(38,n)*V(38,n)+ 0.5*Dn(38,n)*(Is(38,n)+I(62,Sm,n+1)+I(63, Sm,n+1)-I(64,1,n+1)+I(62,Sm,n)+I(63,Sm,n)-I(64,1,n)); V(39,n+1)=CCn(39,n)*V(39,n)+ 0.5*Dn(39,n)*(Is(39,n) +I(64,Sm,n+1)+I(65, 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m 11 of 17 Sm,n+1)-I(66,1,n+1)+I(64,Sm,n)+I(65,Sm,n)-I(66,1,n)); V(40,n+1)=CCn(40,n)*V(40,n)+ 0.5*Dn(40,n)*(Is(40,n) +I(66,Sm,n+1)+I(67, Sm,n+1)-I(68,1,n+1)+I(66,Sm,n)+I(67,Sm,n)-I(68,1,n)); V(41,n+1)=CCn(41,n)*V(41,n)+ 0.5*Dn(41,n)*(Is(41,n) +I(68,Sm,n+1)+I(69, Sm,n+1)-I(70,1,n+1)+I(68,Sm,n)+I(69,Sm,n)-I(70,1,n)); V(42,n+1)=CCn(42,n)*V(42,n)+ 0.5*Dn(42,n)*(Is(42,n) +I(70,Sm,n+1)+I(71, Sm,n+1)-I(72,1,n+1)+I(70,Sm,n)+I(71,Sm,n)-I(72,1,n)); V(43,n+1)=CCn(43,n)*V(43,n)+ 0.5*Dn(43,n)*(Is(43,n)+I(72,Sm,n+1)+I(73, Sm,n+1)+I(72,Sm,n)+I(73,Sm,n)); V(44,n+1)=CCn(44,n)*V(44,n)+ 0.5*Dn(44,n)*(Is(44,n)+I(74,Sm,n+1)+I(75, Sm,n+1)-I(73,1,n+1)+I(74,Sm,n)+I(75,Sm,n)-I(73,1,n)); V(45,n+1)=CCn(45,n)*V(45,n)+ 0.5*Dn(45,n)*(Is(45,n)+I(76,Sm,n+1)+I(77, Sm,n+1)-I(75,1,n+1)+I(76,Sm,n)+I(77,Sm,n)-I(75,1,n)); V(46,n+1)=CCn(46,n)*V(46,n)+ 0.5*Dn(46,n)*(Is(46,n)+I(78,Sm,n+1)+I(79, Sm,n+1)-I(77,1,n+1)+I(78,Sm,n)+I(79,Sm,n)-I(77,1,n)); V(47,n+1)=CCn(47,n)*V(47,n)+ 0.5*Dn(47,n)*(Is(47,n)+ I(80,Sm,n+1)+I(81, Sm,n+1)-I(79,1,n+1)+ I(80,Sm,n)+I(81,Sm,n)-I(79,1,n)); V(48,n+1)=CCn(48,n)*V(48,n)+ 0.5*Dn(48,n)*(Is(48,n)+I(82,Sm,n+1)+I(83, Sm,n+1)-I(81,1,n+1)+I(82,Sm,n)+I(83,Sm,n)-I(81,1,n)); V(49,n+1)=CCn(49,n)*V(49,n)+ 0.5*Dn(49,n)*(Is(49,n)+I(84,Sm,n+1)-I (83,1,n+1)+I(84,Sm,n)-I(83,1,n)); %******************************* % Gan lai gia tri dau va cuoi moi de tinh dien ap trung % binh tren moi U(1,1,n+1)=V(1,n+1); U(4,1,n+1)=V(1,n+1); U(1,Sm+1,n+1)=V(2,n+1); U(5,1,n+1) =V(2,n+1); U(2,1,n+1) =V(2,n+1); U(2,Sm+1,n+1)=V(3,n+1); U(3,Sm+1,n+1)=V(3,n+1); U(7,1,n+1) =V(3,n+1); U(10,1,n+1)=V(3,n+1); U(4,Sm+1,n+1)=V(4,n+1); U(3,1,n+1) =V(4,n+1); U(12,1,n+1)=V(4,n+1); U(5,Sm+1,n+1)=V(5,n+1); U(6,1,n+1) =V(5,n+1); U(13,1,n+1)=V(5,n+1); U(6,Sm+1,n+1)=V(6,n+1); U(7,Sm+1,n+1)=V(6,n+1); U(8,1,n+1) =V(6,n+1); U(15,1,n+1)=V(6,n+1); U(8,Sm+1,n+1)=V(7,n+1); U(9,Sm+1,n+1)=V(7,n+1); U(17,1,n+1)=V(7,n+1); U(20,1,n+1)=V(7,n+1); 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m U(10,Sm+1,n+1)=V(8,n+1); U(11,Sm+1,n+1)=V(8,n+1); U(9,1,n+1) =V(8,n+1); U(22,1,n+1) =V(8,n+1); U(12,Sm+1,n+1)=V(9,n+1); U(11,1,n+1) =V(9,n+1); U(24,1,n+1) =V(9,n+1); U(13,Sm+1,n+1)=V(10,n+1); U(25,1,n+1) =V(10,n+1); U(14,1,n+1) =V(10,n+1); U(14,Sm+1,n+1)=V(11,n+1); U(15,Sm+1,n+1)=V(11,n+1); U(27,1,n+1) =V(11,n+1); U(16,1,n+1) =V(11,n+1); U(16,Sm+1,n+1)=V(12,n+1); U(17,Sm+1,n+1)=V(12,n+1); U(29,1,n+1) =V(12,n+1); U(18,1,n+1) =V(12,n+1); U(18,Sm+1,n+1)=V(13,n+1); U(19,Sm+1,n+1)=V(13,n+1); U(31,1,n+1) =V(13,n+1); U(34,1,n+1) =V(13,n+1); U(20,Sm+1,n+1)=V(14,n+1); U(21,Sm+1,n+1)=V(14,n+1); U(19,1,n+1) =V(14,n+1); U(36,1,n+1) =V(14,n+1); U(22,Sm+1,n+1)=V(15,n+1); U(23,Sm+1,n+1)=V(15,n+1); U(21,1,n+1) =V(15,n+1); U(38,1,n+1) =V(15,n+1); U(24,Sm+1,n+1)=V(16,n+1); U(23,1,n+1) =V(16,n+1); U(40,1,n+1) =V(16,n+1); U(25,Sm+1,n+1)=V(17,n+1); U(41,1,n+1) =V(17,n+1); U(26,1,n+1) =V(17,n+1); U(26,Sm+1,n+1)=V(18,n+1); U(27,Sm+1,n+1)=V(18,n+1); U(43,1,n+1) =V(18,n+1); U(28,1,n+1) =V(18,n+1); U(28,Sm+1,n+1)=V(19,n+1); U(29,Sm+1,n+1)=V(19,n+1); U(30,1,n+1) =V(19,n+1); U(45,1,n+1) =V(19,n+1); 12 of 17 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m U(30,Sm+1,n+1)=V(20,n+1); U(31,Sm+1,n+1)=V(20,n+1); U(47,1,n+1) =V(20,n+1); U(32,1,n+1) =V(20,n+1); U(32,Sm+1,n+1)=V(21,n+1); U(33,Sm+1,n+1)=V(21,n+1); U(49,1,n+1) =V(21,n+1); U(52,1,n+1) =V(21,n+1); U(34,Sm+1,n+1)=V(22,n+1); U(35,Sm+1,n+1)=V(22,n+1); U(33,1,n+1) =V(22,n+1); U(54,1,n+1) =V(22,n+1); U(36,Sm+1,n+1)=V(23,n+1); U(37,Sm+1,n+1)=V(23,n+1); U(35,1,n+1) =V(23,n+1); U(56,1,n+1) =V(23,n+1); U(38,Sm+1,n+1)=V(24,n+1); U(39,Sm+1,n+1)=V(24,n+1); U(37,1,n+1) =V(24,n+1); U(58,1,n+1) =V(24,n+1); U(40,Sm+1,n+1)=V(25,n+1); U(39,1,n+1) =V(25,n+1); U(60,1,n+1) =V(25,n+1); U(41,Sm+1,n+1)=V(26,n+1); U(61,Sm+1,n+1)=V(26,n+1); U(42,1,n+1) =V(26,n+1); U(42,Sm+1,n+1)=V(27,n+1); U(43,Sm+1,n+1)=V(27,n+1); U(63,1,n+1) =V(27,n+1); U(44,1,n+1) =V(27,n+1); U(44,Sm+1,n+1)=V(28,n+1); U(45,Sm+1,n+1)=V(28,n+1); U(65,1,n+1) =V(28,n+1); U(46,1,n+1) =V(28,n+1); U(46,Sm+1,n+1)=V(29,n+1); U(47,Sm+1,n+1)=V(29,n+1); U(67,1,n+1) =V(29,n+1); U(48,1,n+1) =V(29,n+1); U(48,Sm+1,n+1)=V(30,n+1); U(49,Sm+1,n+1)=V(30,n+1); U(69,1,n+1) =V(30,n+1); U(50,1,n+1) =V(30,n+1); U(50,Sm+1,n+1)=V(31,n+1); 13 of 17 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m U(51,Sm+1,n+1)=V(31,n+1); U(71,1,n+1) =V(31,n+1); U(74,1,n+1) =V(31,n+1); U(52,Sm+1,n+1)=V(32,n+1); U(53,Sm+1,n+1)=V(32,n+1); U(51,1,n+1) =V(32,n+1); U(76,1,n+1) =V(32,n+1); U(54,Sm+1,n+1)=V(33,n+1); U(55,Sm+1,n+1)=V(33,n+1); U(53,1,n+1) =V(33,n+1); U(78,1,n+1) =V(33,n+1); U(56,Sm+1,n+1)=V(34,n+1); U(57,Sm+1,n+1)=V(34,n+1); U(55,1,n+1) =V(34,n+1); U(80,1,n+1) =V(34,n+1); U(58,Sm+1,n+1)=V(35,n+1); U(59,Sm+1,n+1)=V(35,n+1); U(57,1,n+1) =V(35,n+1); U(82,1,n+1) =V(35,n+1); U(60,Sm+1,n+1)=V(36,n+1); U(59,1,n+1) =V(36,n+1); U(84,1,n+1) =V(36,n+1); U(61,Sm+1,n+1)=V(37,n+1); U(62,1,n+1) =V(37,n+1); U(62,Sm+1,n+1)=V(38,n+1); U(63,Sm+1,n+1)=V(38,n+1); U(64,1,n+1) =V(38,n+1); U(64,Sm+1,n+1)=V(39,n+1); U(65,Sm+1,n+1)=V(39,n+1); U(66,1,n+1) =V(39,n+1); U(66,Sm+1,n+1)=V(40,n+1); U(67,Sm+1,n+1)=V(40,n+1); U(68,1,n+1) =V(40,n+1); U(68,Sm+1,n+1)=V(41,n+1); U(69,Sm+1,n+1)=V(41,n+1); U(70,1,n+1) =V(41,n+1); U(70,Sm+1,n+1)=V(42,n+1); U(71,Sm+1,n+1)=V(42,n+1); U(72,1,n+1) =V(42,n+1); U(72,Sm+1,n+1)=V(43,n+1); U(73,Sm+1,n+1)=V(43,n+1); U(74,Sm+1,n+1)=V(44,n+1); 14 of 17 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m 15 of 17 U(75,Sm+1,n+1)=V(44,n+1); U(73,1,n+1) =V(44,n+1); U(76,Sm+1,n+1)=V(45,n+1); U(77,Sm+1,n+1)=V(45,n+1); U(75,1,n+1) =V(45,n+1); U(78,Sm+1,n+1)=V(46,n+1); U(79,Sm+1,n+1)=V(46,n+1); U(77,1,n+1) =V(46,n+1); U(80,Sm+1,n+1)=V(47,n+1); U(81,Sm+1,n+1)=V(47,n+1); U(79,1,n+1) =V(47,n+1); U(82,Sm+1,n+1)=V(48,n+1); U(83,Sm+1,n+1)=V(48,n+1); U(81,1,n+1) =V(48,n+1); U(84,Sm+1,n+1)=V(49,n+1); U(83,1,n+1) =V(49,n+1); for nh=1:Snh for m=1: Sm % Xac dinh ap trung binh tren cac segment Vi Vi(nh,m,n+1)=0.5*(U(nh,m,n+1)+U(nh,m+1,n+1)); % Xac dinh dong dien tan tren cac segment Itan(nh,m,n+1)=G(nh,m,n)* Vi(nh,m,n+1); end end for nh=1:Snh for m=1: Sm if Itan(nh,m,n+1)> IO for k=n+1:N+1 if r(nh,m,k)1 % Tai cac vi tri giao giua hai segment G(nh,m,n+1)=0.5*(Gi(nh,m-1,n+1)+Gi(nh,m,n+1)); C(nh,m,n+1)=0.5*(Ci(nh,m-1,n+1)+Ci(nh,m,n+1)); CC(nh,m,n+1)=(1-G(nh,m,n+1)*dt/(2*C(nh,m,n+1)))/(1+G(nh,m,n+1)*dt/(2*C(nh,m, n+1))); D(nh,m,n+1)=1.0/((G(nh,m,n+1)/2+C(nh,m,n+1)/dt)*dx); end end end end for n=1:N+1 Ztan(1,n)=V(chon,n)/Iset(1,n); end tg=[0:N]*dt; t=1:N; figure(1) hold on; grid on; plot(tg(t),V(13,t)','r'); xlabel('Thoi gian'); ylabel('Dien Ap (kV)'); title('DIEN AP V13 LUOI 6x6 KHI DONG SET VAO TAI 13'); %legend('V1','V10','V11','V12','V13',-1); figure(2) hold on; grid on; for n=1:N+1 Ztan(1,n)=V(chon,n)/Iset(1,n); end plot(tg(t),Ztan(1,t)','r'); xlabel('Thoi gian (s)'); ylabel('Dien tro (ohm)'); title('DIEN TRO TAN XUNG LUOI 6X6 KHI DONG SET VAO VI TRI 13'); figure(3) hold on; grid on; for m=1: Sm Ztan(m,t)=I(31,m,t); end plot(tg(t),Ztan(1:Sm,t)','r'); plot(tg(t),IO','b'); xlabel('Thoi gian (s)'); ylabel('Dong dien (kA)'); title('PHAN BO DONG DIEN TAN TREN THANH SO 31 LUOI 6X6 KHI SET VAO NUT 13'); figure(4) hold on; grid on; for m=1: Sm 11/29/08 7:14 PM E:\6x6_ion\Luoi_6x6_thv13_ionhoa.m Ztan(m,t)=r(31,m,t); end plot(tg(t),Ztan(1:Sm,t)','r'); %plot(tg(t),IO','b'); xlabel('Thoi gian (s)'); ylabel('ban kinh phan doan (m)'); title('BAN KINH CAC PHAN DOAN THANH SO 31 LUOI 6X6 KHI SET VAO NUT 13'); 17 of 17 =================================================================== TÀI LIỆU THAM KHẢO [1] Hoàng Việt, Kỹ thuật điện cao áp – Tập : Quá điện áp hệ thống điện Nhà xuất Đại học Quốc Gia TpHCM, 2004 [2] Tiêu chuẩn xây dựng Việt Nam 46 : 2007 : chống sét cho cơng trình xây dựng Hướng dẫn thiết kế, kiểm tra bảo trì hệ Thống [3] Các tài liệu tình hình giơng sét Việt Nam Viện Vật Lý Địa Cầu [4] Yaquing Liu, Transient response of Grounding system caused by Lightning : Modelling and Experiments, Acta Universitatis Upsaliensis Uppsala 2004 [5] L V Bewley, “ Theory and tests of the counterpoise”, Elec Engr., Vol 53, pp 1163-1172, Aug 1934 [6] A P Meliopoulos and M G Moharam, “Transient analysis of grounding systems”, IEEE Trans Power Apparatus and Systems, Vol.PAS-102, No.2, pp 389399, 1983 [7]M Ramamoorty, M M Babu Narayanan, S Parameswaran et al., ”Transient performance of grounding grids”, IEEE Trans Power Delivery, Vol 4, pp 2053-2059, October 1989 [8] A Geri, ”Behaviour of grounding systems excited by high impulse currents: the model and its validation”, IEEE Trans Power Delivery, Vol 14, No 3, pp 1008 1017, July 1999 [9] A F Otero, J Cidras and J L del Alamo, “Frequency-dependent grounding system calculation by means of a conventional nodal analysis technique”, IEEE Tansactions on Power Delivery, Vol 14, No 3, pp 873-878, July, 1999 [10] J Cidras, A F Otero and C Garrido, “Nodal frequency analysis of grounding systems considering the soil ionization effect”, IEEE Trans Power Delivery, Vol 15, No.1, pp 103-107, Jan 2000 [11] L Grcev and F Dawalibi, ”An electromagnetic model for transients in grounding system”, IEEE Trans Power Delivery, Vol 5, pp 1773-1781, November 1990 [12] L Grcev, ” Computation of transient voltages near complex grounding systems caused by lightning currents”, Proceedings of IEEE 1992 International Symposium on EMC, 92CH3169-0, p 393-399 [13] L D Grcev and F E Menter, ”Transient electromagnetic fields near large earthing systems”, IEEE Trans Magnetics, Vol 32, No 3, May 1996, p 1525-1528 =================================================================== =================================================================== [14] L D Grcev, ”Computer analysis of transient voltages in large grounding systems”, IEEE Trans Power Delivery, Vol 11, pp 815-823, April 1996 [15] B Nekhoul, C Cuerin, P Labie, G Meunier and R Feuillet, “A finiteelement method for calculating the electromagnetic fields generated by substation grounding systems”, IEEE Transactions on Magnetics, Vol 31, No 3, pp 2150-2153, May, 1995 [16] B Nekhoul, P Labie, F X Zgainski and G Meunier, “Calculating the impedance of a grounding system”, IEEE Transactions on Magnetics,Vol 32, No 3, pp 15091512, May, 1996 [17] F Dawalibi, ” Electromagnetic fields generated by overhead and buried short conductors, part I- single conductor”, IEEE Trans On Power Delivery, Vol PWRD-1, No 4, 1986, p 105-111 [18] F Dawalibi, ” Electromagnetic fields generated by overhead and buried short conductors, part II- ground networks”, IEEE Trans On Power Delivery, Vol PWRD1, No 4, 1986, p 112-119 [19] R Andolfato, L Bernardi and L Fellin, “ Aerial and grounding system analysis by the shifting complex images method”, IEEE Trans Power Delivery, Vol 15, No 3, pp 10011009, July 2000 [20] Y Liu, N Theethayi, R Gonzalez and R Thottappillil, “The residual resistivity in soil ionization region around grounding system for different experimental results”, IEEE International Symposium on EMC, Boston, paper no TH-PM-2-4, Aug 18th22nd, 2003, pp 794-799 [21] P L Bellaschi, R E Arminggton and A E Snowden, ” Impulse and 60-cycle characteristics of driven grounds II ”, AIEE Transactions, Vol 61, p 349-363, 1942 [22] A C Liew and M Darveniza, ” Dynamic model of impulse characteristics of concentrated earths ”, Proc IEE, Vol 121, No 2, p 123-135, Feb 1974 [23] Y Liu, N Theethayi, R Gonzalez and R Thottappillil, “The residual resistivity in soil ionization region around grounding system for different experimental results”, IEEE International Symposium on EMC, Boston, paper no TH-PM-2-4, Aug 18th22nd, 2003, pp 794-799 [24] M Loboda, Z Pochanke, “ Experimental study of electric properties of soil with impuse current injections”, 18th ICLP Munich, p 191-198, 1985 [25] I.F Gonos and I.A Stathopulos, Soil ionisation under lightning impulse voltages, IEEE Proc.-Sci Meas Technol., Vol 151, No 5, September 2004 =================================================================== =================================================================== [26]Gao Yanquing, ZENG Rong, HE Jinliang, LIANG Xidong, “ Loss Transmission line Model of Grounding Electrodes Considering Soil Ionization of Lightning Impulse, 0-7803-7277, IEEE 2002 [27] Yaqing Liu, Nelson Theethayi, and Rajeev Thottappillli, “ An Engineering Model for Trasient Analysis of Grounding System Under Lightning Strikes : Nonuniform Transmission – Line Approach IEEE Trans, Power Delivery, Vol.20, No 2, pp.722730 2005 [28] Vu Phan Tu, Nguyen Nhat Nam, Huynh Quoc Viet, “Finite Difference Method and ATP/EMTP for selected EMC problem in power systems:modeling of grounding systems transients”, International Symposium on Electrical & Electronics Engineering 2007 - Oct 24, 25 2007 - HCM City, Vietnam .[29] M I Lorentzou, N D Hatziargyriou, Senior Member, IEEE, and B C Papadias, Life Fellow, IEEE “Time Domain Analysis of Grounding Electrodes Impulse Response”, IEEE transactions on power delivery, vol 18, no 2, april 2003 =================================================================== LÝ LỊCH TRÍCH NGANG Họ tên : Hồ Đức Huấn Ngày, tháng, năm sinh : 30/04/1976 Nơi sinh : Gia Lai Địa liên lạc : 200 - Cách mạng tháng - Thành phố Pleiku - Tỉnh Gia Lai Điện thoại di động : 0963535222 Điện thoại nhà : 059.3827173 ™ Quá trình học tập cơng tác: Năm 1994-1999: sinh viên khóa 1994 – Khoa Điện – Điện tử - Trường Đại học Bách khoa thành phố Hồ Chí Minh Tháng 05/2000 đến nay: cơng tác Phân xưởng phát điện Biển Hồ - Điện lực Gia Lai Năm 2006-2008: học viên cao học khóa 2006 – Ngành Thiết bị, Mạng Nhà máy điện – Trường Đại học Bách khoa thành phố Hồ Chí Minh ... đất sét 40 Chương Ứng dụng phương pháp sai phân hữu hạn vào mơ hình đường dây truyền tải nghiên cứu độ hệ thống nối đất 41 3.1 Phương pháp sai phân hữu hạn 41 3.2 Mơ hình. .. nghiên cứu độ hệ thống nối đất ================================================================ Chương ỨNG DỤNG PHƯƠNG PHÁP SAI PHÂN HỮU HẠN VÀO MƠ HÌNH ĐƯỜNG DÂY TRUYỀN TẢI NGHIÊN CỨU QUÁ ĐỘ CỦA HỆ... I - TÊN ĐỀ TÀI: Ứng dụng mơ hình đường dây truyền tải phương pháp sai phân hữu hạn vào việc nghiên cứu độ hệ thống nối đất II - NHIỆM VỤ VÀ NỘI DUNG: Tổng quan : a) Về sét tình hình dơng sét Việt