Tài liệu Giáo trình mô hình hóa máy điện chương 1

1 CHƯƠNG 1: MÔ HÌNH HOÁ MÁY BIẾN ÁP §1. MÔ HÌNH MBA HAI DÂY QUẤN 1. Phương trình từ thông: Như ta đã biết trong phần máy điện, nếu bỏ qua dòng điện từ hoá ta có: i1W1 + i2W2 = 0 (1) i1 W = − 2 hay: i2 W1 Tỉ số biến đổi điện áp là: e 1 W1 ( dΦ m dt ) W1 = = (2) e 2 W2 ( dΦ m dt ) W2 Sau khi biến đổi ta có: e 1i 1 = − e 2 i 2 (3) Tổng trở của m.b.a sau khi quy đổi là: 2 W  Z 1 =  1  Z 2 (4) W  2 Từ thông trong mba bao gồm từ thông trong lõi thép Φ m , từ thông tản của cuộn sơ cấp Φ σ 1 và của cuộn thứ cấp Φ σ 2 . Như vậy từ thông của cuộn sơ cấp sẽ là: Φ 1 = Φ m + Φ σ 1 (5) và của cuộn dây thứ cấp: Φ 2 = Φ m + Φ σ 2 (6) Từ thông móc vòng với cuộn dây sơ cấp: λ 1 = W1Φ 1 = W1 ( Φ m + Φ σ 1 ) Từ thông tản Φ σ 1 tạo bởi s.t.đ của cuộn dây sơ cấp và từ thông hỗ cảm Φ (7) m tạo bởi s.t.đ của cả hai cuộn dây nên ta có thể viết lại biểu thức (7) dưới dạng:   λ 1 = W1  W1i1ρ σ 1 + ( W1i1 + W2 i 2 ) ρ m  = ( W12ρ σ 1 + W12ρ m ) i1 + W1W2ρ m i 2 (8) 14 2 43  14 2 43 1 4 44 2 4 4 43  1 4 44 2 4 4 43 L12 Φm L11  Φ σ1  Trong đó ρ σ 1 và ρ m là độ dẫn từ của mạch từ tản và của mạch từ chính (lõi thép). Tương tự, từ thông móc vòng với cuộn dây thứ cấp là: λ 2 = W2  ( W1i1 + W2 i 2 ) ρ m + W2 i 2ρ σ 2  = ( W22ρ σ 2 + W22ρ m ) i 2 + W1W2ρ m i1 14 2 43 (9) 1 4 44 2 4 4 43 L 22 Phương trình từ thông móc vòng với 2 cuộn dây được viết lại là: λ 1 = L11i1 + L12 i 2 λ 2 = L 21i1 + L 22 i 2 L 21 (10) (11) Trong đó L11 và L22 là hệ số tự cảm của các cuộn dây. L12 và L21 là hệ số hỗ cảm giữa chúng. Hệ số tự cảm của cuộn dây sơ cấp L 11 có thể được viết thành tổng của hệ số tự cảm ứng với từ trường tản và hệ số tự cảm ứng với thành phần từ hoá. Như vậy với i 2 = 0 ta có: 2 λ 1i L11 = i1 Trong đó Φ Tương tự: L 22 = 2= 0 m1 = W1 (Φ σ1 + Φ m1 ) i1 = W12ρ σ 1 + W12ρ m 123 123 Lσ 1 L m1 (12) = W1i1ρ m là phần từ thông trong lõi thép được tạo bởi dòng điện i1. λ 2i 1= 0 i2 = W2 (Φ σ2 i2 + Φ m2 ) = W22ρ σ 2 + W22ρ m 123 123 Lσ 2 L m2 (13) Trong đó Φ m 2 = W2 i 2ρ m là phần từ thông trong lõi thép được tạo bởi dòng điện i2. So sánh Lm1 và Lm2 ta có: 2 Lm 2 W  =  2  L m1  W1  (14) Tổng từ thông hõ cảm móc vòng với mỗi dây quấn có thể biểu diễn bằng biểu thức sau:    W2  W1Φ m = W1 ( Φ m1 + Φ m 2 ) = L m1  i1 + i2  (15) W 1   { i′2   2. Phương trình điện áp: S.đ.đ cảm ứng trong dây quấn sơ cấp là: dλ 1 di di e1 = − = −  L11 1 + L12 2  (16) dt dt dt   Thay L11 bằng (12) ta có: di W  d e 1 = − L σ 1 1 − L m 1  i 1 + 2 i 2  (17) dt dt  W1  W2 i 2 ta có: Sử dụng dòng điện quy đổi i′2 = W1 e1 = − L σ 1 d ( i1 + i′2 ) di1 − Lm1 dt dt Tương tự s.đ.đ cảm ứng trong cuộn dây thứ cấp là:  di dW e 2 = − L σ 2 2 − L m 2  1 i 1 + i 2  dt dt  W2  (18) (19) Quy đổi sang cuộn sơ cấp ta có: di′ d e ′2 = − L′σ 2 2 − L m 1 ( i 1 + i′2 ) (20) dt dt Phương trình cân bằng điện áp trên cuộn dây sơ cấp là; d ( i1 + i′2 ) di u1 = i1r1 − e1 = i1r1 + L σ 1 1 + L m1 (21) dt dt Và phương trình điện áp thứ cấp đã quy đổi có dạng: di′ d( i 1 + i′2 ) u′2 = i′2 r2′ + L′σ 2 2 + L m 1 (22) dt dt Từ các phương trình trên ta đưa ra được sơ đồ thay thế quen thuộc như đã biết trong giáo trình máy điện. Các thông số quy đổi của dây quấn thứ cấp là: 3 2 W  r2′ =  1  r2  W2  (23) W  L′σ 2 =  1  L σ 2  W2  (24) 2 i1 r1 u1 Lσ1 r2′ i1 + i′2 e1 Lm em L′σ 2 i′ 2 e′2 u′2 i2 u2 §2. MÔ PHỎNG MBA 1. Các phương trình: Ta sẽ sử dụng phương trình từ thông và điện áp để mô phỏng m.b.a. Tất nhiên sẽ có nhiều cách mô phỏng khác nhau. Ở đây ta sẽ dùng từ thông móc vòng với 2 dây quấn làm biến trạng thái. Lúc này phương trình điện áp có thể viết lại thành: 1 dΨ 1 u1 = i1r1 + (25) ω b dt 1 dΨ ′2 u′2 = i′2 r2′ + (26) ω b dt Trong đó Ψ 1 = ω b λ 1 , Ψ ′2 = ω b λ 2 và ω b là tần số cơ bản dùng để tính toán điện kháng. Do đó theo (21) và (22) ta có: Ψ 1 = ω b λ 1 = xσ 1i1 + Ψ m Ψ ′2 = ω b λ ′2 = x′σ 2 i′2 + Ψ m Ψ m = ω b L m 1 (i 1 + i′2 ) = x m 1 (i 1 + i′2 ) Dòng điện sơ cấp và thứ cấp có thể biểu diễn bằng biểu thức: Ψ − Ψm i1 = 1 xσ 1 Ψ′ − Ψm i′2 = 2 x′σ 2 Thay các dòng điện vào (29) ta có: Ψ m Ψ 1 − Ψ m Ψ ′2 − Ψ m = + x m1 xσ 1 x′σ 2 (27) (28) (29) (30) (31) (32) Hay:  1 1 + +   x m1 x σ 1 1 1 1 = + + x M x m1 x σ 1 Ψ Đặt: m Viết gọn lại ta có: 1  Ψ 1 Ψ ′2 + = x′σ 2  x σ 1 x′σ 2 1 x′σ 2 (33) (34) 4  Ψ Ψ′  = xM  1 + 2  (35)  x σ 1 x′σ 2  Cuối cùng ta nhận được:   Ψ − Ψ m   Ψ 1 =  ω b u1 − ω b r1  1 (36)   dt x σ 1       Ψ ′ − Ψ m   Ψ ′2 =  ω b u′2 − ω b r2′  2 (37)   dt x′σ 2     Tập hợp các phương trình (30), (31), (35), (36) và (37) tạo ra mô hình động học cơ bản của m.b.a hai dây quấn. Tính phi tuyến của mạch từ và tổn hao công suất có thể thêm vào khi cần thiết. Trong mô hình này, các từ thông móc vòng là biến bên trong. Điện áp trên các đầu là các biến vào và các dòng điện là các biến ra. Trong sơ đồ khối này, các biến vào là điện áp tức thời của các dây quấn. 2. Điều kiện tải: Sơ đồ trên sử dụng điện áp trên các cực làm đại lượng đầu vào để mô phỏng và tạo ra các đại lượng đầu ra là các dòng điện. Điện áp đưa vào cuộn sơ cấp u1 hoặc là giá trị cố định hoặc là nhận được từ việc mô phỏng các phần tử khác nối với dây quấn sơ cấp. Điều kiện ngắn mạch dây quấn thứ cấp dễ dàng được mô phỏng khi đặt điện áp thứ cấp u′2 = 0 . Mô phỏng điều kiện làm việc không tải khó hơn. Điều kiện không tải của m.b.a thể hiện bằng i′2 = 0 ở dây quấn thứ cấp. Khi thay 1 dΨ m vào (26) và (28) cho ta u′20 = . Để tránh việc thực hiện đạo hàm Ψm khi mô ω b dt dΨ 1 phỏng ta tính điện áp thứ cấp khi không tải từ giá trị trước khi tích phân nó để dt tạo ra Ψ1. Các quan hệ được sử dụng là các quan hệ giữa Ψm và Ψ1 trong (26), (28) và (29) với điều kiện i′2 = 0 : 1 dΨ m 1 x m1 dΨ 1 x m1 u′20 = = = ( u1 − i1r1 ) (38) ω b dt ω b xσ 1 + x m1 dt xσ 1 + x m1 Với một tải xác định trên cuộn thứ cấp, mô phỏng sẽ dễ dàng khi tải có thể biểu diễn bằng một tổng trở hay một tổng dẫn. Giả sử ta có một tải có dung lượng St và điện áp thứ cấp bằng điện áp định mức. Khi đó ta có tổng trở của phụ tải là: Ψ m ∫ ∫ 2 1  W1  U 22dm Z= =  = (G′ + jB′ )− 1  Y  W2  S t Tổng dẫn này có thể mô tả bằng mạch điện như hình a hay hình b. i′2 u′2 i′L i′2 i′R L’ R’ u′2 i′C i′R C’ R’ Phương trình của các tảia tương đương này có thể viết dưới dạng tích phân với điện b áp là đại lượng ra và dòng điện là đại lượng vào để bổ sung vào các phương trình ở đầu ra của dây quấn thứ cấp. 5 Khi tải có tính cảm ta mô tả bằng một mạch điện như hình a (các đại lượng đã quy đổi về sơ cấp). Điện áp thứ cấp là: u′2 = i′R R ′ = − ( i′2 + i′L ) R ′ (40) Trong đó i′2 là dòng điện đầu ra của dây quấn thứ cấp m.b.a nhận được khi mô phỏng và i′L là dòng điện nhận được khi tích phân điện áp trên L’. i′L = 1 u′2 dt = ω b B′ ∫ u′2 dt L′ ∫ (41) Khi tải có tính dung ta dùng sơ đồ mô tả tải như hình b. Điện áp thứ cấp sẽ có dạng: u′2 = 1 ω i′Cdt = b ∫ C′ B′ ∫  ′ u′2   − i 2 − ′  dt R   (42) i i′2 Khi cần liên kết hai hay nhiều phần của mạch ta nên thêm vào một điện iH Khối mô trở lớn hay một điện dung nhỏ để tạo Máy ′2 R u phỏng biến áp ra u′2 . Điện trở giả tưởng thêm vào RH H khác lớn hơn tổng trở của các phần tử thực của mạch cho phép ta tạo ra điện áp c u′2 mà không làm cho sai số mô phỏng tăng lên. Điện áp vào đòi hỏi để liên kết các modul bằng điện áp trên điện trở giả tưởng: u′2 = i H R H = − ( i′2 + i ) R H (43) Tương tự, điện dung nhỏ CL có thể được dùng thay cho RH và: 1 u′2 = − ( i′2 + i )dt (44) CL ∫ 3. Bão hoà của mạch từ: Các m.b.a thường làm việc trong trạng thái mạch từ bị bão hoà. Sự bão hoà mạch từ ảnh hưởng chủ yếu đến điện kháng hỗ cảm và ảnh hưởng ít hơn đến điện kháng tản. Điện kháng tản chỉ có thể xác định được khi biết cấu trúc của máy và điều này trong nhiều trường hợp không thể thực hiện được. Do đó khi mô phỏng, ta chỉ cần tính đến ảnh hưởng của sự bão hoà của mạch từ đến điện kháng hỗ cảm bằng cách: • Sử dụng giá trị điện kháng từ hoá bão hoà thích hợp tại mỗi thời điểm mô sat phỏng. Theo phương pháp này, ta cập nhật giá trị điện kháng từ hoá bão hoà x m 1 khi unsat mô phỏng bằng cách dùng tích của giá trị điện kháng từ hoá không bão hoà x m 1 và hệ số bão hoà ks. Cả hai giá trị này có thể được xác định xác định từ thí nghiệm không tải. Trong điều kiện không tải, điện áp rơi trên tổng trở của cuộn dây sơ cấp r1 + jx σ 1 nhỏ nên ta có thể bỏ qua. Do vậy U1 ≈ E1 = Em. Khi từ thông biến thiên hình sin, như ta thấy từ (25) và (29), Em = Ψm. Như vậy trục điện áp không tải trên đồ thị hình a cũng có thể coi như trục từ thông móc vòng Ψmm. Độ dốc của phần tuyến tính của đặc tính unsat không tải là giá trị không bão hoà của điện kháng không tải x m 1 . Điện kháng bão 6 sat hoà x m 1 ở một điện áp bất kì trên đặc tính không tải bằng độ dốc của đường thẳng nối điểm đó với gốc. Ψm Ψ msat Ψ munsat Ψ msat ∆Ψ a 45 0 i1 b Độ bão hoà có thể xác định bằng hệ số bão hoà: Ψ sat I unsat mhd k s = unsat = msat ks ≤ 1 Ψ mhd Im Nếu điện kháng từ hoá bão hoà hiệu dụng x ks = Ψ I sat mhd sat m Ψ I sat m unsat mhd x sat m1 = unsat x m1 sat m1 Ψ munsat (45) Ψ sat m được coi là tỉ số sat thì: Im (46) Đối với một số phương pháp mô phỏng, ví dụ mô phỏng tương tự, việc dùng các điện kháng hằng số trong (35) sẽ dễ dàng hơn dùng điện kháng biến đổi để tính đến sự bão hoà của mạch từ. Thông thường khi mô phỏng như vậy, giá trị hiện hành của unsat Ψ sat được tính m sẽ được xác định từ giá trị không bão hoà của từ thông hỗ cảm Ψ m unsat bằng cách dùng giá trị x m 1 . Ta sẽ xác định sai khác giữa giá trị từ thông bão hoà và không bão hoà: Ψ munsat = Ψ msat + ∆ Ψ m (47) Giá trị ∆Ψm dương trong góc phần tư thứ nhất nhưng âm trong góc phần tư thứ 3. unsat sat Quan hệ giữa ∆Ψm và Ψ m hay Ψ m có thể suy ra từ đường cong không tải của m.b.a. Như đã thấy từ hình a, với một dòng điện không tải cho i1 cho trước, ta có thể unsat sat xác định các giá trị tương ứng của Ψ m và Ψ m . Lặp lại các bước này cho các giá trị sat unsat i1 khác nhau ta có đường cong Ψ m = f( Ψ m ) và ∆ Ψ m (hình b). Xấp xỉ đường cong từ hoá bằng một số hàm giải tích. Muốn vậy ta phải xây dựng quan hệ hàm giữa giá trị biên độ của từ thông và giá trị biên độ của dòng điện. Vì thí nghiệm không tải thường được thực hiện bằng cách đưa điện áp hình sin vào dây quấn sơ cấp và bỏ qua điện áp rơi trên dây quấn nên từ thông trong lõi thép cũng được coi là biến thiên hình sin theo t và dòng điện từ hoá sẽ không hình sin. • Sử dụng quan hệ giữa giá trị từ thông móc vòng bão hoà và không bão hoà. Phương pháp này thích hợp khi ta chọn từ thông làm biến trạng thái. Để dễ hiểu ta thêm chỉ số phụ bên trên để phân biệt giữa giá trị từ thông hỗ cảm bão hoà và không • 7 bão hoà. Quan hệ giữa dòng điện với từ thông móc vòng bão hoà và không bão hoà thể hiện qua quan hệ với từ thông hỗ cảm. Ta viết lại (29): unsat ′ ′ Ψ munsat = ω b Lunsat (48) m1 ( i1 + i 2 ) = x m1 ( i1 + i 2 ) Tương tự, giá trị bão hoà của các dòng điện có thể tính theo từ thông móc vòng bão hoà: Ψ − Ψ msat i1 = 1 (49) xσ 1 Ψ 2′ − Ψ msat i′2 = (50) x′σ 2 Thay giá trị của các dòng điện vào (48) ta có: Ψ unsat Ψ 1 − Ψ sat Ψ ′2 − Ψ sat m m m = + (51) unsat x m1 xσ 1 x′σ 2 unsat Chú ý là các giá trị Ψ1 và Ψ 2′ trong (50) và (51) là các giá trị bão hoà. Thay Ψ m bằng Ψ msat + ∆ Ψ m và nhóm các số hạng Ψ msat ta có: ∆Ψ   Ψ Ψ′ Ψ msat = x M  1 + 2 − unsat  (52)  x σ 1 x′σ 2 x m 1  Trong đó giá trị x M cũng giống như trong phương trình (34) đối với trường hợp không bão hoà, nghĩa là: 1 1 1 1 = unsat + + (53) x M x m1 x σ 1 x′σ 2 Như vậy, muốn tính đến bão hoà, ta cần biết ∆Ψ ở vế phải của (52). Điều này được sat thực hiện nhờ quan hệ hàm giữa ∆Ψ và Ψ m . Sơ đồ mô phỏng việc tính toán này như hình sau. So sánh với sơ đồ đã có trước đây ta thấy sự thay đổi nằm ở số hạng cuối của sat (52) và một modul phụ cần để tính ∆Ψ từ Ψ m . Khi mô phỏng bằng máy tính số, giá trị hiện thời của ∆Ψ có thể được xác định bằng cách nội suy từ bảng số liệu hay đơn sat giản bằng quan hệ hàm gần đúng giữa ∆Ψ và Ψ m trong một phạm vi nào đó. sat Trong SIMULINK, bảng quan hệ ∆Ψ và Ψ m được thực hiện nhờ modul Look-up Table trong thư viện Nonlinear. Quan hệ vào-ra của modul Look-up Table được xác sat định bằng các mảng vào và ra có cùng độ dài. Quan hệ giữa ∆Ψ và Ψ m như hình sau có thế xấp xỉ bằng một hàm đơn giản. ∆Ψ ∆Ψ Đoạn hàm mũ Độ dốc A2 Ψ msat B1 B2 Độ dốc A1 B1 Độ dốc A2 B2 Ψ msat 8 sat Ta có hai ví dụ xấp xỉ ∆Ψ( Ψ m ) bằng 3 đoạn trong góc phần tư thứ nhất. Mô tả toán học của 3 đoạn là: sat Vùng tuyến tính ( Ψ m < B1): Trong phần không bão hoà: ∆Ψ = 0 (54) sat Vùng khuỷu cong (B1 < Ψ m < B2): Vùng này có tính phi tuyến cao. Nó có thể xấp xỉ bằng hàm: sat (55) ∆ Ψ = ae b ( Ψ m − B1 ) Trong đó hằng số b được xác định bằng cách cân bằng biểu thức với giá trị của ∆Ψ tại sat điểm Ψ m = B 2 , nghĩa là: (56) ∆ Ψ = ae b ( B2 − B1 ) sat sat Vùng bão hoà ( Ψ m > B2): Trong vùng này đường cong ∆Ψ( Ψ m ) được xấp xỉ bằng hàm tuyến tính: ∆ Ψ = A 2 ( Ψ msat − B 2 ) + ∆ Ψ ( B 2 ) (57) Mô tả toán học của cách xấp xỉ từng đoạn có thể biểu diễn bằng: ∆ Ψ = A 1 ( Ψ msat − B1 ) + A 2 ( Ψ msat − B 2 ) (58) sat Trong đó A1 sẽ bằng độ dốc 1 nếu Ψ m > B1 và bằng 0 trong các trường hợp khác; A 2 sat bằng (độ dốc 2 - độ dốc 1) nếu Ψ m > B2 và bằng 0 trong các trường hợp khác với độ dốc 1 và B1 là độ dốc và điểm gãy của đoạn thứ 2; độ dốc 2 và B2 là độ dốc và điểm gãy của đoạn thứ 3. Do Ψ m biến đổi, sự bão hoà khi Ψ m âm phải được tính bằng cách xấp xỉ hàm ∆ sat sat Ψ( Ψ m ) trong góc phần tư thứ 3. Với Ψ m < 0 độ dốc của phần tuyến tính không thay sat đổi do đó A không đổi nhưng dấu của điểm gãy B thay đổi theo Ψ m . Bây giờ ta xét đến đường cong bão hoà tính theo các giá trị tức thời. Đường cong từ hoá của m.b.a có được từ thí nghiệm không tải và thể hiện quan hệ U′2 = f(I1 ) . Do tất cả các biến dùng trong mô phỏng là các biến tức thời được quy đổi về sơ cấp nên ∆Ψ phải được biểu diễn bằng các biến tức thời quy đổi về dây quấn sơ cấp. Điện áp hiệu dụng thứ cấp đo được khi hở mạch dễ dàng quy đổi về sơ cấp bằng cách dùng tỉ số các vòng dây, nghĩa là: W  U10 =  1  U 20 (59)  W2  U Ψ Ψ Ψ Un Uk k Ψk Ψ1 Ψ2 Ψ1 Ik a In I sat m Ψ θ1 θ2 π/2 b θ 2 ik in c i 9 Trước hết ta vẽ đường cong không tải theo giá trị hiệu dụng (hình a). Các điểm được đánh số là 1, 1, 2 ,..., n. Điểm 0 nằm tại gốc, điểm1 ở cuối đoạn tuyến tính. Các điểm khác có thể phân bố gần đều trên đoạn bão hoà. Hình c cho thấy các điểm tương ứng trên đường cong giá trị sat tức thời của Ψ m theo i được xác định liên tiếp nhau, mỗi điểm một lần. Tương ứng với mỗi điểm trên đường cong không tải theo giá trị hiệu dụng ta khảo sát đường cong dòng điện không tải hiệu dụng của m.b.a. khi điện áp đưa vào hình sin có biên độ bằng 2 giá trị hiệu dụng của điện áp đặt vào như chỉ trên hình b đối với điểm thứ k. Với điện áp hình sin tần số ω, từ thông móc vòng tương ứng sẽ là hình sin và có giá trị biên độ là: Ψ k = 2U k k = 0, 1,..,n (60) Như vậy, các giá trị của Ψ 0 , Ψ 1 ..., Ψ n trong hình b có thể xác định từ quan hệ trên. Ngoại trừ điểm đầu tiên có i0 = 0, giá trị biên độ của các dòng điện từ hoá i1, i2,...,in khi này vẫn còn chưa biết. Chúng được xác định bằng cách cân bằng các biểu thức của giá trị hiệu dụng của dòng điện trong hình c với các giá trị hiệu dụng đo được tại các điểm tương ứng trong hình b. Khi số điểm được sử dụng đủ lớn và sự phân bố của chúng hợp lí, giá trị hiệu dụng của dòng điện khi điện áp kích thích hình sin có thể được xác định với độ chính xác chấp nhận được bằng cách dùng phương pháp tuyến tính hoá từng đoạn phần đường cong giữa hai điểm cạnh nhau như trên hình c. Gọi Kj là độ dốc của đoạn nối điểm thứ (j - 1) và điểm thứ j đo theo chiều đứng, nghĩa là: i j − i j− 1 Kj = j - 1, 2,..,n (61) Ψ j − Ψ j− 1 Giá trị của ik có thể biểu diễn bởi: ik = k ∑ j= 1 ( Kj Ψ j − Ψ j− 1 ) k = 1, 2,..,n (62) Bắt đầu với j = 1, giá trị biên độ của sóng từ thông móc vòng tương ứng với điểm 1 trên đặc tính không tải hình a là Ψ 1 = 2 U 1 . Với đoạn thẳng đầu tiên đi từ gốc biểu sat thị quan hệ Ψ m (i) , biểu thức giải tích của dòng điện tức thời là: i = K1Ψ1sinθ (63) Giả sử rằng điện áp là hình sin, từ thông móc vòng cũng sẽ hình sin. Khi bỏ qua từ trễ, dòng điện từ hoá sẽ có dạng sóng 1/4 hình sin. Như vậy ta chỉ cần khảo sát 1/4 sóng kích thích khi tính giá trị hiệu dụng. Ví dụ, đối với điểm thứ k, chúng ta chỉ cần khảo sát giá trị hiệu dụng của dòng điện nằm trong vùng gạch chéo như trên hình b. Với k = 1, ta có: I 12 = hay: K 1 = π 2 2 K 12 Ψ 12 2 ( K Ψ sin θ ) d θ = 1 1 π ∫0 2 (64) 2I 1 Ψ1 Tương tự, đối với điểm thứ 2 của hình a, ta sử dụng từ thông Ψ 2 = (65) 2 U 2 và có: 10 π  θ1  2 2  2 2 2 I 2 =  ∫ ( K 1 Ψ 2 sin θ ) dθ + ∫ [ K 1 Ψ 1 + K 2 ( Ψ 2 sin θ − Ψ 1 ) ] dθ  (66) π 0 θ1    −1 Trong đó θ 1 = sin ( Ψ 1 Ψ 2 ) . Do i1 = K1Ψ1 nên (66) có thể viết lại thành phương trình bậc 2 đối với K2: A 2 K 22 + B 2 K 2 + C 2 = 0 (67) Trong đó: A2 = π 2 ∫ ( Ψ 2 sin θ − θ1 B 2 = 2K 1 Ψ 1 Ψ 1 ) 2 dθ π 2 ∫ ( Ψ 2 sin θ − θ1 A2 > 0 Ψ 1 ) dθ B2 > 0 θ 1 π π C 2 = i 12  − θ 1  + ∫ ( K 1 Ψ 2 sin θ 1 ) 2 dθ − I 22 2 2  0 Và chỉ có duy nhất một giá trị dương của K2 là: B 22 − 4 A 2 C 2 K2 = 2A 2 Tương tự, với đoạn có độ dốc Kk ta có: A k K 2k + B k K k + C k = 0 − B2 + Trong đó: Ck = dk + ∑ ( K 2j A j + k− 1 j= 1 (68) ) K jB j + d j − (69) C2 > 0 (70) (71) (72) π 2 Ik 2 2 j− 1 j dj = i t t j = θ j − θ j− 1  Ψj  θ j = sin − 1  Ψk 1 sj = sin 2θ j − sin 2θ j− 1 2 g j = cos 2θ j − cos 2θ j− 1 ( Ψ k2 Aj = t j − s j + 2Ψ k Ψ 2 B j = − 2i j− 1 Ψ k g j + Ψ j− 1 t j ( ( ) ) ) j− 1 (73) gj + Ψ 2 j− 1 j t j = 1,...,k; 1 ≤ k ≤ n Bắt đầu với điểm ở gốc, nghĩa là k = 0, trong đó Ψ 0 = 0 , i 0 = 0 , θ 0 = 0 các giá trị của Kk với k = 1,..,n nhận được khi dùng liên tiếp (72) và (73) như đã thấy trước đây ở (71) với k = 2. Ta dùng file mginit.m, mgplt và smg.mdl dựa trên thuật toán trên để tìm giá sat trị từ thông Ψ m ( i ) từ đường cong không tải. 4. Các bài tập cần làm: a. Mô phỏng m.b.a một pha tuyến b. Mô phỏng m.b.a một pha phi tuyến tính: c. Mô phỏng m.b.a 3 pha nối Y/Y: