Bài viết trình bày ứng dụng phương pháp số để nghiên cứu một số kết cấu rôto, dạng sóng và xác định giá trị sức điện động không tải E0. Kết quả của bài báo sẽ giúp các nhà thiết kế tính toán và điều chỉnh thông số để đạt được giá trị E0 theo yêu cầu.
KHOA HỌC CÔNG NGHỆ P-ISSN 1859-3585 E-ISSN 2615-9619 NGHIÊN CỨU SỨC ĐIỆN ĐỘNG KHÔNG TẢI ĐỘNG CƠ ĐỒNG BỘ NAM CHÂM VĨNH CỬU KHỞI ĐỘNG TRỰC TIẾP CÔNG SUẤT NHỎ BẰNG PHƯƠNG PHÁP SỐ STUDYING PHASE BACK-EMF OF LOW POWER THREE PHASE LINE START PERMANENT MAGNET SYNCHRONOUS MOTORS BY A NUMERICAL TECHNIQUE Lê Anh Tuấn1,*, Bùi Đức Hùng2 KÝ HIỆU TÓM TẮT Động đồng nam châm vĩnh cửu khởi động trực tiếp (LSPMSM) có hiệu suất, hệ số công suất gần 1, đáp ứng tiêu chuẩn yêu cầu ngày cao động điện Vì vậy, LSPMSM chắn giải pháp thay phần động không đồng rơto lồng sóc tương lai Tuy có nhiều ưu điểm vận hành, LSPMSM bị ảnh hưởng từ thông số động kết cấu răng, rãnh, kích thước vật liệu nam châm vĩnh cửu, sức điện động khơng tải E0 Trong q trình thiết kế, việc lựa chọn cấu hình động khơng đảm bảo dẫn đến giá trị sức điện động không tải khơng phù hợp Điều dẫn đến hiệu suất động suy giảm, chí thấp cịn thấp động không đồng công suất Trong nội dung, báo ứng dụng phương pháp số để nghiên cứu số kết cấu rôto, dạng sóng xác định giá trị sức điện động khơng tải E0 Kết báo giúp nhà thiết kế tính tốn điều chỉnh thơng số để đạt giá trị E0 theo yêu cầu Ký hiệu Đơn vị Ý nghĩa E0 V Sức điện động không tải Br T Mật độ từ dư nam châm vĩnh cửu M1 Wb Biên độ thành phần từ thông Id, Iq A Dòng điện dọc trục ngang trục Pem W Công suất điện từ CHỮ VIẾT TẮT LSPMSM Động đồng nam châm vĩnh cửu khởi động trực tiếp NCVC Nam châm vĩnh cửu FDM, FEM Phương pháp sai phân phương pháp phần tử hữu hạn Từ khoá: Động đồng nam châm vĩnh cửu khởi động trực tiếp; nam châm vĩnh cửu; động đồng bộ; sức điện động ABSTRACT Line start permanent magnet synchronous motors with high efficiency and high power factor in steady state are meet the requirements of standards for high efficiency motors They will be a feasible solution for replacing partially squirrel cage induction motors (SCIM) which widely used at the present time Howerver, their performance are mainly depended on the configuration of their rotor configuration, sizes of permanent magnets, permanent magnet materials, phase back-EMF value (E0), In designing process, if the EMF value is not suitable for motor, the efficiency of the motor will be down Even, the efficiency of LSPMSM is lower than SCIM with the same power Therefore, this paper studies the effect of rotor’s configuration on LSPMSM’s phase back-EMF waveform and determines phase back-EMF value This study will help the motor designers computing E0 and modifying configuration of LSPMSM to reach appropriate E0 value Keywords: Line-Start Permanent Magnet Synchronous Motors; Permanent Magnet; Synchronous Motors; Electromotive Force Khoa Điện, Trường Đại học Công nghiệp Hà Nội Viện Điện, Trường Đại học Bách khoa Hà Nội * Email: leanhtuan0985@gmail.com Ngày nhận bài: 01/11/2021 Ngày nhận sửa sau phản biện: 25/12/2020 Ngày chấp nhận đăng: 26/02/2021 GIỚI THIỆU Động đồng nam châm vĩnh cửu chứng minh có ưu điểm vượt trội trình vận hành xác lập, ví dụ có hiệu suất biến đổi điện lớn, hệ số công suất gần 1, mật độ mômen lớn [1, 2] Do phát triển nhanh chóng cơng nghệ nam châm vĩnh cửu đất với mật độ từ trường cao, tích số lượng lớn, chủng loại nam châm Neodimium… nên hai thập kỷ gần thúc đẩy việc ứng dụng chủng loại nam châm lĩnh vực máy điện động có hiệu suất lượng cao Trong loại động sử dụng nam châm vĩnh cửu, LSPMSM lên giải pháp thay phần cho SCIM lĩnh vực công, nông nghiệp LSPMSM chất cải tiến từ động đồng nam châm vĩnh cửu động khơng đồng Nói cách khác lai tạp động khơng đồng động đồng cách đặt nam châm vĩnh cửu vào rôto SCIM Do vậy, LSPMSM sở hữu ưu điểm động đồng nam châm vĩnh cửu, cịn có ưu điểm SCIM khả tự khởi động trực tiếp Đối với LSPMSM khẳng định động có hiệu suất, hệ số cơng suất cao, thiết kế khơng 18 Tạp chí KHOA HỌC VÀ CÔNG NGHỆ ● Tập 57 - Số (02/2021) Website: https://tapchikhcn.haui.edu.vn SCIENCE - TECHNOLOGY P-ISSN 1859-3585 E-ISSN 2615-9619 tốt thơng số vận hành khơng đảm bảo, chí cịn khơng đạt tiêu chuẩn SCIM [3, 4] Một số yếu tố liên quan lớn đến hiệu suất hệ số công suất vận hành LSPMSM đại lượng sức điện động không tải E0 Sức điện động không tải sinh biến thiên từ thông nam châm vĩnh cửu gây dây quấn stato rơto quay Hiện nay, có số phương pháp để xác định sức điện động không tải E0 phương pháp giải tích, phương pháp số ứng dụng phân tích phần tử hữu hạn, phương pháp đo lường thực tế Phương pháp giải tích cho kết tính toán nhanh, cách xác định E0 cho LSPMSM thường áp dụng cơng thức tính tốn IPM với giả thiết rơto khơng bị bão hịa, khơng xét ảnh hưởng rơto tính tốn [7] Vì vậy, phương pháp giải tích thường cho kết với sai số cao Bên cạnh đó, phương pháp giải tích khơng cho phép người sử dụng quan sát dạng sóng sức điện động khơng tải Phương pháp đo lường thực tế cho kết xác nhất, nhiên phải chế tạo mẫu kèm biện pháp đo phức tạp Phương pháp số khẳng định có độ xác cao phân tích tốn điện từ trường máy điện [7, 10] Vì vậy, để nghiên cứu dạng sóng sức điện động khơng tải giá trị nó, báo tập trung ứng dụng phương pháp số khảo sát, tính tốn LSPMSM Bên cạnh đó, báo so sánh kết đạt với kết phương pháp số với kết thu từ phương pháp giải tích phương pháp thực nghiệm thực tế Trong nội dung nghiên cứu, báo khảo sát ba LSPMSM pha ứng với ba công suất khác nhau: 2,2kW, 5,5kW, 11kW Các LSPMSM hiệu chỉnh từ động không đồng Công ty Cổ phần Chế tạo động Hà Nội Trong hiệu chỉnh, kết cấu răng, rãnh, mạch từ, dây quấn, lồng sóc stato rơto giữ ngun, nhiên lõi thép rôto đặt thêm NCVC NdFe45 Từ kết đạt được, báo đánh giá rút số kết luận thiết kế vận hành để đảm bảo hiệu suất, hệ số công suất động XÁC ĐỊNH GIÁ TRỊ SỨC ĐIỆN ĐỘNG KHÔNG TẢI Ở CHẾ ĐỘ XÁC LẬP 2.1 Xác định sức điện động không tải theo công thức giải tích Ở điều kiện xác lập, sức điện động không tải E0 ảnh hưởng lớn đến hiệu suất hệ số công suất động Trong đó, sức điện động khơng tải xác định giá trị sức điện động quay động với tốc độ đồng trạng thái hở mạch không bao gồm thành phần sức điện động sinh dòng điện stato [5] Về lý thuyết E0 xác định từ cấu hình thiết kế động thơng qua phương trình giải tích thông số đầu vào: Mật độ từ dư NCVC, số vòng dây stato, chiều dài lõi từ,… [5, 6] Sức điện động không tải E0 LSPMSM xác định công thức (1), (2) Với giả thiết máy điện hoạt động chế độ xác lập không xảy bão hòa, người ta thường xác định sức điện động không tải LSPMSM tương tự với động đồng nam châm gắn chìm (IPM) [8, 9] E0 xác định sau: E0 2 k w1.Nph f M1 (1) Trong đó: kw1 hệ số dây quấn, Nph số vòng dây quấn pha, f tần số quay, M1 biên độ thành phần từ thông Biên độ thành phần từ thông xác định: M1 C 4.d.l Br sin( ) p. 1 (1 2 ) (2) Trong đó: d chiều rộng khối barier từ, l chiều dài lõi sắt rôto, Br mật độ từ dư nam châm vĩnh cửu hệ số cung cực từ C, , hệ số tính đến từ thơng rị, độ tập trung từ thông [9] 2.2 Xác định sức điện động không tải phương pháp số Có nhiều phương pháp số ứng dụng tính tốn tốn liên quan đến trường điện từ máy điện Các phương pháp phổ biến phương pháp sai phân hữu hạn (FDM - Finite Difference Method), phương pháp thành phần biên (BEM - Boundary Element Method), phương pháp phần tử hữu hạn (FEM - Finite Difference Method), phương pháp rời rạc (Discrete Element Method) Các phương pháp cho phép giải gần tốn phương trình vi phân riêng, ví dụ bốn phương trình Maxwell mơ tả định luật điện từ trường [6] Định luật Faraday cảm ứng điện từ: xE B t Định luật Ampere dòng điện chạy dây dẫn: xH J D t Định luật Gauss cho từ trường: B Hình Từ trường nam châm vĩnh cửu theo tiết diện ngang LSPMSM [8] Website: https://tapchikhcn.haui.edu.vn Định luận Gauss cho điện trường: D Vol 57 - No (Feb 2021) ● Journal of SCIENCE & TECHNOLOGY 19 KHOA HỌC CÔNG NGHỆ P-ISSN 1859-3585 E-ISSN 2615-9619 Trong phương pháp số, phương pháp phần tử hữu hạn khẳng định có độ xác cao ứng dụng giải gần phương trình vi tích phân toán liên quan đến điện từ trường [7] Vì vậy, nghiên cứu báo sử dụng phần mềm ứng dụng FEM để khảo sát sức điện động khơng tải E0 Bên cạnh đó, báo ứng dụng phần mềm Matlab để phân tích Fourier (FFT) tìm giá trị biên độ sóng thành phần sóng hài kèm sức điện động khơng tải Các kết sau so sánh với phương pháp giải tích đề cập mục 2.1 với phương pháp đo lường kiểm nghiệm thực tế Trong nội dung, báo trình bày phân tích ba LSPMSM có cơng suất tương ứng: 2,2kW; 7,5kW; 11kW, chi tiết mục bên Dải công suất LPMSM thực nghiệm, nhóm tác giả lựa chọn thuộc dải công suất nhỏ mà SCIM sử dụng phổ biến [8] Bốn phân vùng thể cực từ tương ứng từ trường NCVC sinh Hình Đường từ trường mơ LSPMSM 2,2kW phương pháp số không tải - Phân bố mật độ từ trường khe hở không khí LSPMSM 2,2kW hình XY Plot Curve Inf o Brad Setup1 : Transient Time='0ns' 0.50 0.25 0.00 -0.25 -0.50 -0.75 0.00 Nhóm tác giả ứng dụng phương pháp số để xác định sức điện động không tải loại động có cơng suất: 2,2 kW, 7,5kW 11kW pha, điện áp /Y-220/380V Các thơng số động cho bảng Bảng Các thơng số LSPMSM thử nghiệm Đường kính stato Đường kính ngồi stato Số rãnh stato Số vòng dây/pha Lớp dây quấn Khe hở khơng khí Số rãnh rơto Vật liệu NCVC Mật độ từ dư NCVC Bề rộng NCVC Độ dày NCVC Br wm lm T mm mm 2,2kW 104 170 36 312 0,5 28 NdFeN45 1,38 34 7,5kW 50.00 75.00 100.00 125.00 Kết mô LSPMSM với phần mềm ứng dụng phương pháp số mô động sau: - Đường từ trường mô động cơng suất 2,2kW biểu diễn hình - Từ kết quả, thấy đường từ trường mơ hình thành nên bốn phân vùng từ trường đối xứng 175.00 Hình cho thấy đường phân bố mật độ từ trường khe hở không khí lúc khơng phải dạng xung vng lý tưởng mà bị chia cắt theo kết cấu rãnh stato, rơto - Dạng sóng sức điện động khơng tải E0 chế độ xác lập LSPMSM 2,2kW mơ hình XY Plot LSPMSM_2.2kW_EMF 250.00 Curve Info InducedVoltage(PhaseA) Setup1 : Transient 11kW 252 180 140 272 36 48 156 144 1 0,35 0,5 36 40 NdFeN45 NdFeN45 1,38 1,38 44 74 5 150.00 Hình Phân bố mật độ từ trường LSPMSM 2,2kW không tải InducedVoltage(PhaseB) Setup1 : Transient InducedVoltage(PhaseC) Setup1 : Transient 125.00 Y1 [V] Công suất Ký Đơn vị hiệu D mm Dn mm Z1 N1 vòng Lớp g mm Z2 25.00 Dis tance [mm] DẠNG SÓNG VÀ GIÁ TRỊ SỨC ĐIỆN ĐỘNG KHÔNG TẢI E0 3.1 Áp dụng phương pháp số xác định E0 Maxwell2DDesign1 0.75 Brad Trong phương pháp số, phương pháp FDM FEM sử dụng phổ biến phân tích kết cấu điện từ máy điện Các phương pháp thực toán liên quan đến tượng vật lý cơ, nhiệt, từ điện Nói chung, hai phương pháp dựa rời rạc hóa khơng gian kết cấu thành bề mặt vuông tam giác nhỏ hay thành thể tích khối tứ diện, lục diện Trong đó, đỉnh hình khối gọi nút biểu diễn giá trị biến hệ Độ xác phương pháp liên quan trực tiếp đến độ mịn mắt lưới 0.00 -125.00 -250.00 0.00 10.00 20.00 30.00 Time [ms] 40.00 50.00 60.00 Hình Sức điện động khơng tải E0 LSPMSM 2,2kW Dạng sóng sức điện động khơng tải E0 pha khơng có dạng hình sin mà bị vạt đầu, dạng sóng sức điện động khơng tải E0 dao động điều hòa với tần số 50Hz Tuy nhiên, thực tế LSPMSM nói riêng động sử dụng nguồn điện xoay chiều nói chung, thiết kế để vận hành với sóng sin (điện áp, dịng điện) Vì dạng sóng điều hịa phân tích FFT để tìm giá trị biên độ sức điện động thành phần đồng thời phân tích đánh giá thành phần điều hịa bậc cao kèm 20 Tạp chí KHOA HỌC VÀ CÔNG NGHỆ ● Tập 57 - Số (02/2021) Website: https://tapchikhcn.haui.edu.vn SCIENCE - TECHNOLOGY P-ISSN 1859-3585 E-ISSN 2615-9619 3.2 Phân tích Fuorier dạng sóng sức điện động khơng tải E0 Từ kết mô giá trị chuỗi dạng sóng E0 thu từ mơ với Ansys/Maxwell, nhóm tác giả ứng dụng phần mềm Matlab để phân tích FFT tìm biên độ sóng biên độ thành phần sóng điều hịa kèm Thuật tốn tính ứng dụng phân tích FFT viết Maltab sau: >> L= length(E0_11kW) >> t=0:(1/L)*20*10^-3:20*10^-3*(L-1)/L; >> Fs=50*1000 ; rơto LSPMSM, việc áp dụng cơng thức tính tốn có sai số định ảnh hưởng răng, rãnh gây Sai số tính tốn sức điện động khơng tải FEM phương trình giải tích lớn ứng với động 11kW (tương ứng 7,1%), kết với sai số chấp nhận Bảng Sức điện động không tải E0 chế độ xác lập với FEM giải tích Động FEM Matlab Giải tích Sai số 2,2kW 173,9 185,1 6,1% 7,5kW 166,9 179,6 7,1% 11kW 198,7 214,8 7,5% 3.4 Kiểm nghiệm thực tế >> nfft=L; Để đánh giá thử nghiệm LSPMSM thực tế, nhóm tác giả thực nghiệm mơ hình LSPMSM 2,2kW cải tạo từ động không đồng pha 3K112-S4 Công ty Cổ phần chế tạo điện Hà Nội Sơ đồ thử nghiệm LSPMSM hình >> f=(0:nfft/2-1)*Fs/nfft; >> E0_2.2kW=fft(E0_2.2kW,nfft)/(L/2); >> E0_7.5kW=fft(E0_7.5kW,nfft)/(L/2); >> E0_11kW=fft(E0_11kW,nfft)/(L/2); >> mxE0_2.2kW=abs(E0_2.2kW); >> mxE0_7.5kW=abs(E0_7.5kW); >> mxE0_11kW=abs(E0_11kW); Kết phân tích FFT sóng điều hịa sức điện động khơng tải ứng với ba LSPMSM khảo sát hình 300 E0 2,2 kW E0 7,5 kW E0 11 kW Søc ®iƯn ®éng (V) 250 Hình Sơ đồ thử nghiệm LSPMSM 2,2kW Trong sơ đồ thử nghiệm, nghiên cứu sử dụng động đồng pha tốc độ 1.500 vg/ph động sơ cấp để kéo LSPMSM 2,2kW thử nghiệm quay với tốc độ đồng LSPMSM lúc trở thành máy phát Không nối với tải, tiến hành đo điện áp pha LSPMSM, giá trị đo sức điện động không tải E0 Mô hình triển khai đo lường thực tế hình 200 150 100 50 0 50 100 150 200 250 TÇn sè (Hz) 300 350 400 450 Hình Phân tích FFT dạng sóng E0 cho động khác Thành phần sức điện động LPSMSM thử nghiệm có sai khác cấu hình răng, rãnh kích thước nam châm vĩnh cửu lựa chọn cho động khác Khi phân tích FFT thấy thành phần sóng hài bậc lớn, trường hợp lớn LSPMSM 11kW lúc biên độ thành phần sóng bậc xấp xỉ 20% sóng 3.3 So sánh kết tính tốn E0 từ FEM giải tích Trong nghiên cứu, nhóm tác giả so sánh kết tính tốn sức điện động khơng tải E0 thu từ phương pháp giải tích đề cập mục 3.1 kết thu từ phân tích dạng sóng E0 với FEM Kết so sánh bảng Sai số phương pháp giải tích tính tốn với FEM giải thích ngun nhân sức điện động khơng tải tính theo phương pháp giải tích suy từ cách tính cho động đồng nam châm gắn chìm Tuy nhiên rôto IPM răng, rãnh Website: https://tapchikhcn.haui.edu.vn Hình Cấu hình rơto mơ hình thử nghiệm LSPMSM 2,2kW Đồng hồ vạn sử dụng đo điện áp VOM Digital Multimeter Victor VC890C+ cho phép đo giá trị hiệu dụng thực (true RMS) Kết thực tế đo so sánh với kết thu từ tính tốn giải tích phần mềm ứng dụng phương pháp FEM, kết bảng Bảng So sánh sai số phương pháp xác định E0 FEM Giải tích Đo lường thực tế Sức điện động (V) 173,9 185,1 165 Sai số so với đo lường thực tế 5,1% 10,9% Vol 57 - No (Feb 2021) ● Journal of SCIENCE & TECHNOLOGY 21 KHOA HỌC CÔNG NGHỆ Từ kết bảng thấy phương pháp xác định sức điện động không tải E0 FEM cho kết xác so với phương pháp giải tích Cụ thể trường hợp LSPMSM 2,2kW, so sánh kết tính tốn từ FEM thử nghiệm thực tế cho sai số 5,1% giá trị thu từ phương pháp giải tích sai số 10,9% KẾT LUẬN Bài báo nghiên cứu sức điện động không tải E0 LSPMSM phần mềm ứng dụng FEM, từ kết mơ báo phân tích FFT dạng sóng E0 với Matlab để xác định biên độ thành phần sóng thành phần sóng hài kèm Trong nội dung nghiên cứu, nhóm tác giả khảo sát sức điện động không tải E0 với ba LSPMSM khác có cơng suất tương ứng: 2,2kW; 7,5kW; 11kW Kết thu tính tốn sức điện động khơng tải E0 với FEM Matlab so sánh với kết tính tốn từ phương pháp giải tích Kết quả, so sánh cho thấy sai số hai phương pháp chấp nhận (sai số tính tốn hai phương pháp lớn ứng với LSPMSM 11kW, sai số 7,5%) Bên cạnh đó, nhóm tác giả tiến hành đo lường LSPMSM 2,2kW thử nghiệm để đánh giá kết với thực tế Kết cho thấy, sai số xác định E0 so với thực nghiệm phương pháp FEM Matlab 5,1% sai số tính tốn với phương pháp giải tích 10,9% Do kết luận, ứng dụng phương pháp FEM Matlab để xác định E0 LSPMSM xác so với phương pháp giải tích Như vậy, ứng dụng FEM Matlab công cụ hữu hiệu giúp người thiết kế xác định E0 với độ xác cao, từ điều chỉnh thơng số kết phù hợp Các biện pháp để điều chỉnh E0 thay đổi số vịng dây/pha, thay đổi vật liệu NCVC, kích thước NCVC thay đổi kết cấu răng, rãnh rôto P-ISSN 1859-3585 E-ISSN 2615-9619 [5] T.J.E Miller, 1989 Brushless Permanent-Magnet and Reluctance Motor Drives Claredon Press-Oxford [6] J.F.Gieras, M Wing, 2002 Permanent Magnet Motor Technology Design and Applications Marcel Dekker, Inc [7] V Elistatova, 2016 Optimal Design of Line Start Permanent Magnet Synchronous Motors of High Efficiency Doctoral Thesis, Université Lille-de-France [8] R T Ugale, Gaurav Singh, Sriniva Baka, B N Chaudhari, 2009 Effective Energy Conservation for the Agricultural Sector using Line Start Permanent Magnet Synchronous Motors IEEE Region 10 Conference TENCON, pp 1-5 [9] C.C Hwang, S.M Chang, C.T Pan, T.Y Chang, 2002 Estimation of Parameters of Interior Permanent Magnet Synchronous Motors Journal of Magnetism and Magnetic Materials 239, pp 600 -603 [10] CompleteMaxwell2D_V15, 2010 AUTHORS INFORMATION Le Anh Tuan1, Bui Duc Hung2 Faculty of Electrical Engineering, Hanoi University of Industry School of Electrical Engineering, Hanoi University of Science and Technology Bên cạnh đó, khảo sát sức điện động khơng tải E0 phương pháp FEM Matlab cịn cho phép người sử dụng mơ dạng sóng tính tốn thành phần sóng hài kèm Cụ thể, ba động thử nghiệm thành phần sóng hài bậc ba lớn đạt giá trị cao 7,5% biên độ sóng ứng với LSPMSM 11kW TÀI LIỆU THAM KHẢO [1] H Isfahani, S V Zadeh, M A Rahman, 2011 Evaluation of Synchronization Capability in Line Start Permanent Magnet Synchronous Motors IEEE International Electric Machines&Drives Conference, pp 1346-1350 [2] R Sadoughi, M Zare, M Aziz, 2015 Comparison Between Line Start PM Synchronous Motor and Induction Motor With Same Nominal Power and Same Pole Pairs When Fed by VF Control Drive International Journal Electronics Communication and Computer Engineering, Vol 6, Iss [3] A Vagati, G Franceschini, I Marongiu, G P Troglia, 1992 Design Criteria of High Performance Synchronous Reluctance Motors IEEE Industry Applications Society Annual Meeting [4] Abdolamir Nekoubin, 2011 Design a Line Start Synchronous Motor and Analysis Effect of the Rotor Structure on the Efficiency World Academy of Science, Engineering and Technology 57 22 Tạp chí KHOA HỌC VÀ CƠNG NGHỆ ● Tập 57 - Số (02/2021) Website: https://tapchikhcn.haui.edu.vn ... lập, sức điện động không tải E0 ảnh hưởng lớn đến hiệu suất hệ số cơng suất động Trong đó, sức điện động không tải xác định giá trị sức điện động quay động với tốc độ đồng trạng thái hở mạch không. .. động không tải sinh biến thiên từ thông nam châm vĩnh cửu gây dây quấn stato rơto quay Hiện nay, có số phương pháp để xác định sức điện động không tải E0 phương pháp giải tích, phương pháp số ứng... độ từ dư nam châm vĩnh cửu hệ số cung cực từ C, , hệ số tính đến từ thơng rị, độ tập trung từ thơng [9] 2.2 Xác định sức điện động không tải phương pháp số Có nhiều phương pháp số ứng dụng