电机设计matlab程序
matlab2010b电机仿真
matlab2010b电机仿真在MATLAB 2010b中进行电机仿真可以使用Simulink和Simscape Power Systems工具箱。
下面是一个简单的示例，演示了如何使用Simulink建立并运行一个电机仿真模型。
1. 打开MATLAB 2010b软件。
2. 在Simulink库浏览器中找到电机建模组件。
可以在“Simulink”标签下的“Electrical”部分找到一些相关组件，如“Induction Motor”和“DC Motor”等。
3. 双击相应的组件，将其拖动到模型编辑器中。
可以使用这些组件来构建一个电机模型。
4. 连接电机建模组件的输入和输出端口。
例如，可以将一个输入信号传递到电机的控制端口，将输出信号连接到电机的运动端口。
5. 配置电机的参数。
双击电机组件，可以打开参数对话框，并设置电机的参数，如转矩、速度、电压等。
6. 添加其他必要的组件和连接，以完成电机模型的搭建。
7. 单击模型编辑器中的“运行”按钮，开始仿真电机模型。
可以通过观察仿真结果和信号波形来分析电机的行为和性能。
注意：MATLAB 2010b版本可能需要安装额外的工具箱才能进行电机仿真。
可以在软件安装目录下的“toolbox”文件夹中查找相关的工具箱，并确保其已经安装和激活。
此外，Simscape Power Systems工具箱也提供了一系列电机模型和仿真组件，可以用于建立更精细和复杂的电机系统模型。
可以按照类似的步骤，使用Simscape Power Systems工具箱中的电机组件进行仿真。
希望以上信息对您有帮助！。
MATLAB用于电机电磁计算编程
采 取 的步骤 是 ： ｉ
（） １ 建立 新 文件 “ｗ 编 制 函数 ： ｎｔｎＹ＝ａ ； ａ ”， ｆ ｃｉ ｗ ｕ ｏ （ｉｘ） 其 中 ｘ 和 ｘ 对应 于计 算 中 的 和 Ｆ １ ｘ， ， ｊ ｉ ｊ Ｅ。 ：
文件 中 编 人 了 三 种 材 料 ： １一Ｄ ４ ； １ ２ Ｗ５ ０ ３ ５一 ！
使 用 ＭＡ Ｌ Ｂ编写 电机设 计 程 序 ， 常 编写 一 ＴＡ 通 个主程 序文件 和 一个 辅 助 文件 夹 （ 内放 置需 要 的各 通 过 ｐｔ（ 线 ） 置 使 其 链 接 能 通 过 主 程 序 方 便 ａ 路 ｈ 设 加 以 引用 。编 写 主程序 的过程 可 以将一 个例 题边 输 入边 验 算 ， 时 检 验 程 序 。 ＭＡ Ｌ Ｂ 没 有 Ｆ Ｒ 及 ＴＡ Ｏ． Ｔ Ａ Ｃ等 语 言等 要 求 的严格 的语 法 规 则 ， Ｒ Ｎ、 因此 没 有 多少计 算机语 言基 础 的人都 能很 快掌 握 和使用 。
ＭＡ Ｌ Ｂ是 一 种 功能 强 大 的解 释 性语 言 ， ＴＡ 它提
供 丰富 的运 算符 及 数 据结 构 、 图形 功 能 、 功能 箱 ， 远 超过 Ｆ Ｒ Ｒ Ｎ、 Ｏ Ｔ Ａ Ｃ等语 言 ； 它 的灵 活 简便 又超 过 而 Ｂ ＳＣ等语 言 。ＭＡ Ｌ Ｂ语 法 限制 不 严 格 ， 别 适 ＡＩ ＴＡ 特 用于 电机 电磁 计 算 的计 算 机 程 序 ， 序 的 自由度 和 程 可移植 性 好 ， 可 通 过 Ｎ ｔｂｏ 程 序 移 到 Ｗｏｄ 还 ｏ ｏｋ将 ｅ ｒ 环境 中享用 。 以下就程 序 编制 的一些 心得 与读 者进
摘 等等 。
matlab-电机控制综合仿真任务书
任务一：利用Simpowersystem里面自带的DC电机模块，完成他励直流电机单闭环调速仿真，速度调节用PI控制方法，要求封装PI模块，给定速度100rad/s，负载由空载到1s时跳变到20N。
调节不同的PI参数，观察仿真结果总结速度波形、转矩波形的变化规律。
另外要求将scope图中的4条曲线参数导出到工作空间，并用subplot和plot函数画在同一个窗口中，每个子图加上对应的标题。
写报告时必须仿真模型图和绘制的曲线图电机相关参数如下：任务二：利用Simulink 基本模块搭建他励直流电机闭环调速系统直流电机子模块，根据以下电机数学模型搭建e M a T C I =Φe L m d T T B J dtωω--= 电机模块要求封装，参数20.05kg m J =⋅，0.02N m s m B =⋅⋅，165m C =，0.01Wb f Φ=点击封装模块时输入。
PI 控制模块要求封装，调节不同的PI 参数，得到不同的速度波形。
写报告时，要给出仿真模型和不同PI 参数的速度波形。
任务三：修改任务二的模型，使其PI 控制参数来自外部变量设定，将速度偏差的仿真结果导出到工作空间（workspace ），用m 文件编写一段循环程序调用任务二的仿真模型（调用函数为sim ），在m 文件中实现对PI 参数的调整，处理数据，得到满意的控制参数，绘制仿真结果曲线图。
M 文件中包含：设定PI 参数――调用仿真模型――速度偏差的最后100个数据点平均值是否小于某个设定值eps ――不满足，则修改PI 参数（用循环的方式修改）……最后用plot 画出速度和转矩波形，用subplot 将两个图画在一个窗口。
求数组长度函数length()，求平均值函数mean()，Subplot ，plot任务四：异步电机调速方法仿真从Help－Demos－SimPowerSystems－General Demos中调出AC Motor Drive - Vector Control (discrete) （power_acdrive.mdl）和Asynchronous Machine （power_pwm.mdl），运行看看结果，并找到相应的模块在哪，并试着自己建立同样的模型。
基于MATLAB的直流电机双闭环调速系统设计
摘要直流电机双闭环调速系统是一个复杂的自动控制系统，是目前直流调速系统中的主流设备，具有调速范围宽、平稳性好、稳速精度高等优点，在理论和实践方面都是比较成熟的系统，在拖动领域中发挥着极其重要的作用。
由于直流电机双闭环调速是各种电机调速系统的基础，本文从直流电机的基本特性进行双闭环直流电机设计，最后用实际系统进行工程设计，并采用MATLAB/SIMULINK进行仿真。
对于直流电机双闭环调速系统，在设计和调试过程中有大量的参数需要计算和调整，运用传统的设计方法工作量大，系统调试困难。
本文对直流电机双闭环调速系统进行辅助设计，选择调节器结构，进行参数计算和近似校验，根据给出和计算出的相应参数，建立起制动、抗电网电压扰动和抗负载扰动的MATLAB/SIMULINK仿真模型，分析转速和电流的仿真波形，并进行调试，使直流电机双闭环调速系统趋于完善、合理。
仿真结果证明了该方法的可行性和合理性。
关键词：直流调速，双闭环系统，电流调节器，转速调节器，计算机仿真AbstractThe double closed loop direct current velocity modulation system is a complex automatic control system, is in the present direct current velocity modulation system mainstream equipment, has the velocity modulation scope width, the stability is good, the steady fast precision higher merit, in the theory and the practice aspect all is the quite mature system, in drives in the domain to play the extremely vital role.Because the direct current machine double closed loop velocity modulation is each kind of electrical machinery velocity modulation system foundation, from the direct current machine basic characteristic conducts the double closed loop direct current machine design method research, finally uses the actual system to carry on the engineering design, and uses MATLAB/Simulink to carry on the simulation.Regarding the double closed loop direct current velocity modulation system, has the massive parameters in the design and the debugging process to need to calculate and to adjust, the utilization tradition design method work load is big, system debugging difficulty.This article carries on the assistance design to the double closed loop direct current velocity modulation system, chooses the regulator structure, carries on the parameter computation and the approximate verification, according to produces the corresponding parameter which and calculates, establishes applies the brake, the anti- electrical network voltage perturbation and the anti-load perturbation MATLAB/Simulink simulation model, the analysis rotational speed and the electric current simulation profile, and carries on the debugging, enable the double closed loop direct current velocity modulation system to tend to the consummation, is reasonable. The simulation result has proven this method feasibility and the rationality.Keywords:Speed control of DC-drivers，Double-closed-loop，Current regulator，Speed regulator，Computer simulation目录摘要 (I)Abstract........................................................... I I 目录........................................................... I II 1．绪论. (1)1.1 课题背景 (1)1.2 课题研究的目的和意义 (1)1.3 论文的主要内容 (2)2．直流电动机调速系统 (4)2.1 直流电动机简介 (4)2.1.1 直流电动机的工作原理 (4)2.1.2 直流电动机的运行特性 (5)2.1.3 直流电动机的起动与调速 (6)2.2 转速控制的要求和调速指标 (7)3．方案选择及系统工作原理 (9)3.1 电动机参数及设计要求 (9)3.2 方案选择及系统框图 (9)3.2.1 方案一：直流电机单闭环调速系统 (9)3.2.2 方案二：直流电机双闭环调速系统 (10)3.2.3 方案三：双闭环脉宽调速系统 (11)3.3 系统工作原理简介 (11)3.3.1 双闭环调速系统静态特性 (11)3.3.2 双闭环系统启动过程分析 (14)3.3.3 双闭环调速系统的动态抗扰动性能 (17)3.3.4 双闭环调速系统中两个调节器的作用 (18)4．双闭环调速系统的设计 (19)4.1 双闭环直流调速系统总体设计方案 (19)4.2 主电路设计与参数计算 (20)4.2.1 主电路原理图 (20)4.2.2 整流变压器的设计 (21)4.2.3 晶闸管元件选择 (23)4.2.4 电抗器参数的计算 (24)4.2.5 励磁电路 (26)4.2.6 三相桥式全控整流电路 (26)4.2.7 晶闸管触发电路 (28)4.3 直流调速系统的保护 (30)4.3.1 过电压保护 (30)4.3.2 电流保护 (33)4.4 控制电路设计 (34)4.4.1 电流调节器的设计 (35)4.4.2 转速调节器的设计 (37)5．调速系统的仿真 (40)5.1.1 MATLAB简介 (40)5.1.2 MATLAB的安装 (41)5.1.3 MATLAB的启动运行 (41)5.1.4 MATLAB的帮助文件 (41)5.1.5 MATLAB所定义的特殊变量及其意义 (41)5.1.6 MATLAB工具箱及SIMULINK简介 (42)5.2 调速系统仿真模型的建立 (43)5.3 仿真结果 (43)5.4 仿真结果分析 (46)结论 (47)参考文献 (48)致谢 (49)附录 A (50)1．绪论1.1 课题背景直流调速是现代电力拖动自动控制系统中发展较早的技术。
matlab_simulink_永磁同步电机_概述及解释说明
matlab simulink 永磁同步电机概述及解释说明1. 引言1.1 概述在电力传动领域中，永磁同步电机已成为一种重要的电机类型。
相比于传统的感应电机和直流电机，永磁同步电机具有高效率、高功率密度和较低的维护成本等优势。
随着现代工业对能源效率和环境保护的日益重视，永磁同步电机在工业应用中得到了广泛的推广和应用。
本文将介绍永磁同步电机及其与Matlab Simulink的结合。
首先，我们将简要介绍Matlab Simulink软件以及其在工程领域中的应用。
接下来，我们将详细介绍永磁同步电机的基本原理、结构特点以及在工业中的实际应用情况。
然后，我们将重点讲解如何使用Matlab Simulink建模永磁同步电机，并通过仿真设计过程详解该方法的具体操作步骤。
最后，我们将分析仿真结果，评估永磁同步电机性能以及控制策略调整优化方法论述与解释。
1.2 文章结构本文共分为五个部分：引言、Matlab Simulink简介、永磁同步电机简介、Matlab Simulink建模永磁同步电机原理及方法解析以及结论与展望。
在引言部分，我们将概述本文的主要内容和结构安排，为读者提供一个整体的框架。
接下来的各个部分将逐一介绍Matlab Simulink软件、永磁同步电机以及它们之间的关联，并详细解释如何使用Matlab Simulink建模永磁同步电机以及评估其性能和优化控制策略。
最后，我们将总结全文观点并对未来永磁同步电机建模与控制策略设计进行展望。
1.3 目的本文的目的是介绍Matlab Simulink和永磁同步电机，并阐述它们之间的关系。
通过对Matlab Simulink建模永磁同步电机过程的详细解释，读者可以了解到使用该软件进行系统建模和仿真的好处，并且理解永磁同步电机在工业中的应用情况以及其优势和局限性。
此外，我们还将分享一些调整优化方法，帮助读者评估永磁同步电机性能并设计出更高效的控制策略。
通过本文的阅读，读者将对Matlab Simulink和永磁同步电机有更深入的了解，并对未来的相关研究和应用有所展望。
MATLAB电机控制综合仿真实验
MATLAB电机控制综合仿真实验一、他励直流电机单闭环调速仿真实验要求：利用Simpowersystem里面自带的DC电机模块，完成他励直流电机单闭环调速仿真，速度调节用PI控制方法，要求封装PI模块，给定速度100rad/s，负载由空载到1s时跳变到20N。
调节不同的PI参数，观察仿真结果总结速度波形、转矩波形的变化规律（PI参数和超调量、稳定时间、稳态误差、振荡次数）。
另外要求将scope图中的4条曲线参数导出到工作空间，并用subplot和plot 函数画在同一个窗口中，每个子图加上对应的标题。
电机相关参数的设置图：仿真原理图：在仿真试验中需要按照实验要求对PI控制器子系统进行封装，然后更改Kp、Ki参数值的大小。
封装PI模块图如下：Plot绘图程序：>>subplot(411)>> plot(t,W,'r'),title('转速')>> subplot(412)>> plot(t,Ia,'b'),title('电枢电流')>> subplot(413)>> plot(t,Te,'g'),title('转矩')>> subplot(414)>> plot(t,If,'y'),title('励磁电流')速度调节用PI控制方法，给定速度100rad/s，负载由空载到1s 时跳变到20N，调节不同的PI参数，从PI模块封装中调节，修改不同的参数Ki 、Kp观察仿真结果。
Ki=100, Kp=5;050100w (r a d /s )00.51 1.52 2.53 3.54 4.55-2000200I a (A )-202I f (A )-1000100T e (N .m )Ki=2, Kp=1;w (r a d /s)I a (A)00.51 1.52 2.53 3.54 4.55I f (A)00.51 1.52 2.53 3.54 4.55T e (N .m )二、 他励直流电机闭环调速系统仿真实验要求：利用Simulink 基本模块搭建他励直流电机闭环调速系统直流电机子模块，根据以下电机数学模型搭建：电磁转矩公式：e M a T C I =Φ 动力学平衡方程：e L m d T T B J dtωω--=电机模块要求封装，参数20.05kg m J =⋅，0.02N m s m B =⋅⋅，165m C =，0.01Wb f Φ=，恒定负载T L =20N 点击封装模块时输入。
Matlab在电机设计与控制中的应用方法总结
Matlab在电机设计与控制中的应用方法总结一、引言电机设计与控制是电力系统领域的重要研究方向。
随着计算机技术的发展，Matlab作为一种强大的数学软件工具，成为电机设计与控制中不可或缺的工具。
本文将总结Matlab在电机设计与控制中的应用方法，并探讨其在该领域的优势和潜力。
二、电机设计1. 参数计算与仿真电机设计的第一步是参数计算和仿真。
Matlab提供了丰富的数学计算和仿真工具，可用于电机参数的计算和模拟。
例如，可以通过Matlab的符号计算功能，快速且准确地计算电机各种参数，如电机的电感、电阻、转速等。
同时，Matlab的仿真工具Simulink可以方便地建立电机的仿真模型，并进行系统级的仿真分析，有助于优化电机设计。
2. 磁场建模与优化Matlab在磁场建模与优化方面具有强大的功能。
通过使用Matlab，可以建立电机的二维和三维磁场模型，并进行电磁场分析。
同时，Matlab还提供了各种优化算法，如遗传算法和粒子群算法，可以在电机设计中应用，以寻找最优的设计方案。
三、电机控制1. 控制算法设计与仿真Matlab提供了丰富的控制算法工具箱，可以用于电机控制系统的算法设计和仿真。
例如，PID控制器是电机控制中常用的算法之一，通过使用Matlab中的控制系统工具箱，可以方便地设计和调试PID控制器，并进行闭环仿真分析。
此外，Matlab还支持其他高级控制算法，如模糊控制和自适应控制，可以根据实际需求选择适用的控制算法。
2. 控制系统性能评估Matlab可以用于电机控制系统性能的评估。
通过分析系统的响应曲线和频率特性，可以评估系统的稳定性、响应时间、抗干扰能力等性能指标。
Matlab提供了强大的信号处理和频谱分析工具，可以对电机控制系统的信号进行分析，以评估系统性能。
四、应用案例以一台直流电机的设计与控制为例，介绍Matlab的应用方法。
1. 参数计算与仿真首先，使用Matlab的符号计算功能，计算直流电机的电感、电阻等参数。
交流电动机模型与MATLAB实现
直接转矩控制设计
利用Matlab的Simulink工具搭建 直接转矩控制系统，包括转矩和 磁通观测器、滞环比较器、逆变 器模型等，进行仿真分析。
滑模控制设计
在Matlab中编写滑模控制算法， 构建滑模面和控制律，通过仿真 验证滑模控制在交流电动机中的 有效性。
电动机参数调整
根据实验或实际数据，调整电动机模型的参数， 如电阻、电感、转动惯量等，以提高模型的准确 性。
控制策略优化
针对特定应用需求，优化控制策略，如速度控制、 转矩控制等，以改善电动机的性能。
04
交流电动机控制策略及 Matlab实现
控制策略概述
矢量控制（Vector Control）
通过变换电流来控制电动机的转矩和磁通，实现高性能的调速控制。
智能化控制
引入智能控制算法，如神经网络、模糊控制等，实现电动机性能的智 能化优化和控制。
06
实验验证与结果分析
实验平台搭建及参数配置
实验平台搭建
1
2
选择适当的交流电动机，如感应电动机或永磁同 步电动机。
3
搭建控制电路，包括功率变换器（如逆变器）和 控制器。
实验平台搭建及参数配置
实验平台搭建及参数配置
性能优化措施探讨
优化设计
从电动机结构、材料、制造工艺等方面进行优化设计，提高电动机本 身的性能。
控制策略优化
改进控制算法，提高控制精度和响应速度，使电动机在更宽的运行范 围内保持良好的性能。
参数辨识与自适应控制
通过在线参数辨识和自适应控制方法，实时调整控制器参数，以适应 电动机参数变化和外部扰动，提高系统鲁棒性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%电机设计程序clear allformat short em1=3;p=2;f=50%1.额定功率PN=5.5*10^3 ;%2.额定电压（单位V，三角形接法）UN=380;UN0=380;%3.功电流(单位A)IKW=PN/(m1*UN0)%4.效率eta按照技术条件的规定eta=0.875eta=0.855 ;%5.功率因数cos(phi) =0.84,按照技术条件的规定cos(phi)=0.84phi=acos(0.84);cos(phi);%6.极对数p=2p=2;%7.定转子槽数：每极每相槽数取整数。
参考类似规格电机取q1=3,则Z1=2m1pq1，再查表10-8选Z2=32，并采用转子斜槽。
q1=3;Z1=2*m1*p*q1Z2=32 ;%8.定转子每极槽数Zp1=Z1/(2*p)Zp2=Z2/(2*p)%9.确定电机的主要尺寸；一般可参考类似电机的主要尺寸来确定Di1和lef.现按10-2中的KE1=0.0108*log(PN/1000)-0.013*p+0.931P1=KE1*PN/(eta*cos(phi))alphap1=0.68;KNm1=1.10;Kdp1=0.96;A1=25000;Bdelta1=0.69;n1=1450;V=(6.1/(alphap1*KNm1*Kdp1))*(1/(A1*Bdelta1 ))*(P1/n1)D1=0.21;%铁心的有效长度Di1=0.136;lef =V/((Di1)^2)%气隙的确定%参考类似产品或由经验公式（10-10a），得lt=0.115;delta =0.0004lef=lt + 2*deltaD2=Di1-2*delta%转子内径先按转轴直径决定（以后再校验转子轭部磁密）Di2=0.048 ;%11.极距tautau =pi*Di1/(2*p)%12.定子齿距t1t1=(pi*Di1/Z1)%转子齿距t2t2=(pi*D2/Z2)bsk=0.01187;%15.设计定子绕组Nphi11=eta*cos(phi)*pi*Di1*A1/(m1*IKW)%取并联支路a1=1,由式（10-15），可得每槽导体数a1=1;Ns1=47%16.每相串联导体数Nphi1Nphi1=Ns1*Z1/(m1*a1)%每相串联匝数N1N1=Nphi1/2%17.绕组线规设计%初选定子电密J11=5.0A/mm^2,由式（10-16），计算导线并绕根数和每根导线面积的乘积。
J11=5.0;%其中定子电流初步估计值I11=IKW/(eta*cos(phi))Nt1Ac11=I11/(a1*J11)%在附录二中选用截面积相近的铜线，高强度漆包线Nt1=2;d1=0.00090;d=0.00095;Ac11=0.0000006362;Nt1Ac11=0.0000006362*2;%18. 设定子槽形%因定子绕组为原导线散嵌，故采用梨形槽，齿部平行。
初步取Bt11=1.4T,按式（10-18）估计定子齿宽Bt11=1.493641;KFe=0.95 ;bt1=t1*Bdelta1/(KFe*Bt11)%初步取Bj11=1.25T,按式（10-19），估计定子轭部计算高度Bj11=1.25;hj11=tau*alphap1*Bdelta1/(2*KFe*Bj11)%按齿宽和定子轭部计算高度的估算值作出定子槽形如图10-24，槽口尺寸参考类似产品决定，齿宽计算如下：b01=0.0035;b11=0.0067;h01=0.0008;h11=(b11-b01)/2*tan(pi/6);h21=0.0145-h11r21=0.0044;bi11=pi*(Di1+2*h01+2*h21+2*h11)/Z1-2*r21;bi12=pi*(Di1+2*h01+2*h11)/Z1-b11;bt1=0.5*(bi12+bi11)%齿部基本平行，齿宽%19.槽满率%曹面积hs1=h11+h21;h=0.002 ;As=((2*r21+b11)/2)*(hs1-h)+pi*(r21)^2/2%按附录三，槽绝缘采用DMD复合绝缘，,槽楔为h=2mm层压板，槽绝缘占面积Aideltai=0.00023;Ai=deltai*(2*hs1+pi*r21)%槽的有效面积Aef=As-Ai%槽满率sf=Nt1*Ns1*d^2/Aef%20.绕组系数alpha1=(p*360/Z1)Kp1=1;alpha=alpha1*pi/180,q=q1;Kd1=sin(q*alpha/2)/(q*sin(alpha/2))Kdp1=Kd1*Kp1%每相有效串联导体数NefNef=Nphi1*Kdp1%21.设计转子槽形与转子绕组%按式（10-39），预计转子导条电流：Ki=0.89 ;I21=Ki*I11*(3*Nphi1*Kdp1)/Z2%初步取转子导条电密JB1=3.5A/mm^2,于是导条截面积JB1=3.5;AB1=I21/JB1%初步取Bt21=1.3T,估算转子齿宽Bt21=1.3;bt21=(t2/10*Bdelta1)/(KFe*Bt21)%初步取Bj21=1.25T,估计转子轭部计算高度Bj21=1.25;hj21=tau*alphap1*Bdelta1/(2*KFe*Bj21)%为获得比较好的起动性能，采用平行槽，作槽形图如图10-25所示，取槽口尺寸b02=1mm,h02=0.5mmb02=0.001;h02=0.0005 ;%齿壁不平行的槽形的宽度（按3-3）计算如下：b12=0.0055;h12=0.0011;h22=0.026-0.0011;h12=(b12-b02)/2*tan(pi/6);h22=0.023-h12;b22=0.003;a=0.026;b=0.010;bt213=pi*(D2-2*2/3*(h02+h12+h22))/Z2-(b22+(b12-b22)/3)%导条截面积（转子槽面积）AB=(b02+b12)/2*h12+(b12+b22)/2*h22%按式（11-41）估计端环电流IR1=I21*Z2/(2*p*pi)%端环所需面积JR1=0.6*JB1;AR11=IR1/JB1AR11=(a+b)/2*(h02+h12+h22)*10^6AR1=input('输入端环面积AR1=')%其中端环电密JR1=0.6,JB1=2.1A/mm^2.按照工艺要求有所需面积确定端环内径及厚度如图10-25b,得端环面积AE=40*10^(-6)%(二)磁路计算%22KE1=0.918;epsilonl1=1-0.918;E1=(1-epsilonl1)*UN%23Ks(1)=1.2;KNm(1)=1.095;precision=1;while precision >0.01for n=1:3KNm(n)=input('输入波形系数KNm(n)=')alphap(n)=input('输入极弧系数alphap(n)=')phi(n)=E1/(4*KNm(n)*Kdp1*f*N1)%24At1=KFe*lt*bt1*Zp1At2=KFe*lt*bt213*Zp2%25hs1=h21+h01+h11+r21;hs2=h02+h12+h22;hj11=(D1-Di1)/2-hs1+r21/3hj21=(D2-Di2)/2-hs2Aj1=KFe*lt*hj11Aj2=KFe*lt*hj21%26Adelta=tau*lef%27Fs(n)=input('输入Fs(n)=')%28Bdelta(n)=Fs(n)*phi(n)/Adelta%29Bt1=Fs(n)*phi(n)/At1%30Bt2=Fs(n)*phi(n)/At2%31Ht1(n)=input('输入Ht1(n)=');Ht2(n)=input('输入Ht2(n)=');%32Kdelta1=t1*100*(4.4*delta*100+0.75*b01*100)/(t1*100*(4.4*delta*100+0.75*b01* 100)-(b01*100)^2)Kdelta2=t2*100*(4.4*delta*100+0.75*b02*100)/(t2*100*(4.4*delta*100+0.75*b02* 100)-(b02*100)^2)Kdelta=Kdelta1*Kdelta2deltaef=Kdelta*delta%33Lt1=(h11+h21)+1/3*r21Lt2=(h12+h22)%34Lj11=pi*(D1-hj11)/(4*p)Lj21=pi*(Di2+hj21)/(4*p)%35mu0=4*pi*10^(-7);Fdelta(n)=Kdelta*Bdelta(n)*delta/mu0%36Ft1(n)=Ht1(n)*Lt1*100Ft2(n)=Ht2(n)*Lt2*100%37Ks(2*n)=(Fdelta(n)+Ft1(n)+Ft2(n))/Fdelta(n)Ks(2*n+1)=Ks(2*n)-(Ks(2*n)-Ks(2*n-1))/3precision=abs((Ks(2*n-1)-Ks(2*n))/Ks(2*n))endendKs=Ks(6);alphap1=alphap(3);phi=phi(3);Bdelta=Bdelta(3);Fdelta=Fdelta(3);Ft1=Ft1( 3);Ft2=Ft2(3);KNm=KNm(3);Fs=Fs(3);Ht1=Ht1(3);Ht2=Ht2(3);%迭代的出来的结论%38Bj1=0.5*phi/Aj1%39Bj2=0.5*phi/Aj2%40Hj1=input('输入Hj1=')Hj2=input('输入Hj2=')%41hj11/tauCj1=input('根据Bj1和hj11/tau输入Cj1=')Fj1=Cj1*Hj1*Lj11*10^2hj21/tauCj2=input('根据Bj2和hj21/tau输入Cj2=')Fj2=Cj2*Hj2*Lj21*10^2%42F0=Fdelta+Ft1+Ft2+Fj1+Fj2%43Im=2*p*F0/(0.9*m1*N1*Kdp1)%44Ims=Im/IKW%45Xms=(4*f* mu0*m1/pi)*((N1*Kdp1)^2/(Ks*p))*lef*(tau/deltaef)Xmss=Xms*IKW/UN%（三)参数计算%46beta=(2*8+1*7)/(3*9)tauy=pi*(Di1+2*(h01+h11)+h21+r21)*beta/(2*p)d1=0.015;lB=lt+2*d1K0=1.2;l0=lB+K0*tauy%47lE=2*d1+K0*tauy%48Cx=4*pi*f*mu0*(N1*Kdp1)^2*lef*PN/(m1*p*UN^2)%49KU1=1.0;KL1=1.0;lambdaU1=h01/b01+2*h11/(b01+b11)h212r21=h21/(2*r21),b112r21=b11/(2*r21)lambdaL1=input('输入lambdaL1=')lambdas1=KU1*lambdaU1+KL1*lambdaL1%50Xs1s=2*m1*p/(Z1*Kdp1^2)*(lt/lef)*lambdas1*Cx%51sigmas=0.0129;Xdelta1s=(m1/(pi^2)*(tau/deltaef)*(sigmas/((Kdp1^2)*Ks)))*Cx %52XE1s=0.47*(lE-0.64*tauy)*Cx/(lef*Kdp1^2)%53Xsigma1s=Xs1s+Xdelta1s+XE1s%54%55lambdaU2=h02/b02hb2=h22/b22,b1b2=b12/b22lambdaL=input('输入lambdaL=');lambdaL2=2*h12/(b02+b12)+lambdaLlambdas2=lambdaU2+lambdaL2%56Xs2s=2*m1*p/Z2*(lt/lef)*lambdas2%57Z22p=Z2/(2*p)sigmaR=input('根据Z22p输入sigmaR=')Xdelta2s=m1*tau*sigmaR/(pi^2*deltaef*Ks)%58DR=0.107;XE2s=(0.2523*Z2*DR/(2*p*lef*2*p))*(2*m1*p/Z2)%59Xsks=0.5*(bsk/t2)^2*Xdelta2s%60Xsigma2s=(Xs2s+Xdelta2s+Xsks+XE2s)*Cx%61Xsigmas=Xsigma1s+Xsigma2s%62rhow=0.0217*10^(-6);R1=rhow*(2*N1*l0)/(Nt1*Ac11*a1)%63R1s=R1*IKW/UN%64C=1.05;rhoCu=8.9*10^3;GCu=C*l0*Ns1*Z1*Ac11*Nt1*rhoCurhoFe=7.8*10^3;GFe=KFe*lt*(D1+delta)^2*rhoFe%65 AR1代替ARKB=1.04;rhow1=0.0434*10^(-6);R21=rhow1*(KB*lt/AB+Z2*DR/(2*pi*p^2*AR1*10^(-6)))*(4*m1*(N1*Kdp1)^2/Z 2)RB1=rhow1*(KB*lt/AB)*(4*m1*(N1*Kdp1)^2/Z2)RBs=RB1*IKW/UNRR1=rhow1*Z2*DR/(2*pi*p^2*AR1*10^(-6))*(4*m1*(N1*Kdp1)^2/Z2)RRs=RR1*IKW/UNR2s=R21*IKW/UN%66eta(1)=0.855;I1Ps(1)=1/eta(1);precisioneta=1;while precisioneta>0.005for m=1:2I1Ps(m)=1/eta(2*m-1)%67sigma1=1+Xsigma1s/XmssIXs(m)=sigma1*Xsigmas*(I1Ps(m))^2*(1+(sigma1*Xsigmas*I1Ps(m))^2)%68I1Qs(m)=Ims+IXs(m)%69epsilonL(m)=I1Ps(m)*R1s+I1Qs(m)*Xsigma1sKE=1-(I1Ps(m)*R1s+I1Qs(m)*Xsigma1s)%70epsilon0=Ims*Xsigma1s1-epsilon0%71Bt10(m)=(1-epsilon0)/(1-epsilonL(m))*Bt1Ht10(m)=input('输入Ht10(m)=')Bt20(m)=(1-epsilon0)/(1-epsilonL(m))*Bt2Ht20(m)=input('输入Ht20(m)=')%73Bj10(m)=(1-epsilon0)/(1-epsilonL(m))*Bj1Hj10(m)=input('输入Hj10(m)=')%74Bj20(m)=(1-epsilon0)/(1-epsilonL(m))*Bj2Hj20(m)=input('输入Hj20(m)=')%75Bdelta0(m)=(1-epsilon0)/(1-epsilonL(m))*Bdelta%76Ft10(m)=Ht10(m)*Lt1*100%77Ft20(m)=Ht20(m)*Lt2*100%78hj11/tau,Bj10Cj1(m)=input('输入Cj1(m)=')Fj10(m)=Cj1(m)*Hj10(m)*Lj11*100%79hj21/tau,Bj20Cj2(m)=input('输入Cj2(m)=')Fj20(m)=Cj2(m)*Hj20(m)*Lj21*100%80Fdelta0(m)=Kdelta*delta*Bdelta0(m)/mu0%81F00(m)=Fdelta0(m)+Fj20(m)+Fj10(m)+Ft20(m)+Ft10(m) %82Im0=2*p*F00(m)/(0.9*m1*N1*Kdp1)%(四)工作性能计算%83I1s(m)=(I1Ps(m)^2+I1Qs(m)^2)^(0.5)I1(m)=I1s(m)*IKW%84J1(m)=(I1(m)/(a1*Nt1*Ac11))/10^6%85A1(m)=m1*Nphi1*I1(m)/(pi*Di1)%86I2s(m)=(I1Ps(m)^2+IXs(m)^2)^(0.5)I2(m)=I2s(m)*IKW*m1*Nphi1*Kdp1/Z2IR(m)=I2(m)*Z2/(2*pi*p)%87JB(m)=I2(m)/(AB*10^6)JR(m)=IR(m)/(AR1)%88pCu1s(m)=I1s(m)^2*R1spCu1(m)=pCu1s(m)*PN%89pA12s(m)=I2s(m)^2*R2spA12(m)=pA12s(m)*PN%90pss=0.02;ps=pss*PN%91pfw=(3/p)^2*(D1)^4*10^4pfws=pfw/PN%92Bj10phej(m)=input('输入phej(m)=')k2=2;Gj=4*p*Aj1*Lj11*rhoFepFej(m)=k2*phej(m)*GjBt10phet(m)=input('输入phet(m)=')Gt=2*p*At1*Lt1*rhoFek1=2;pFet(m)=k1*phet(m)*GtpFe(m)=pFej(m)+pFet(m)pFes(m)=pFe(m)/PN%93sigmaps(m)=pCu1s(m)+pA12s(m)+pss+pfws+pFes(m) %94PN1s(m)=1+sigmaps(m)%95eta(2*m)=1-sigmaps(m)/PN1s(m)eta(2*m+1)=eta(2*m)+(eta(2*m)-eta(2*m-1))/5precisioneta=abs((eta(2*m)-eta(2*m-1))/eta(2*m))endend%96I1s=I1s(2);I1Ps=I1Ps(2);phi=acos(I1Ps/I1s);cos(phi)%97pFers=((1-1/2)*pFej+(1-1/2.5)*pFet)*10^(-3)/PN*1000sN=pA12s/(1+pA12s+pFers+pss+pfws)%98nN=60*f/p*(1-sN)%99Tms=(1-sN)/(2*(R1s+(R1s^2+Xsigmas^2)^(0.5)))%(五)启动性能%100K(1)=3;%101precisionist=1;K(1)=3;Ist(1)=K(1)*Tms*IKWwhile precisionist>0.03for k=1:3Fst(k)=Ist(2*k-1)*Ns1/a1*0.707*(KU1+Kd1^2*Kp1*Z1/Z2)*(1-epsilon0)^(0.5) beta0=0.64+2.5*(delta/(t1+t2))^(0.5)BL=mu0*Fst(k)/(2*delta*beta0)%102KZ(k)=input('输入启动漏抗饱和KZ（k）=')%103cs1(k)=(t1-b01)*(1-KZ(k))cs2(k)=(t2-b02)*(1-KZ(k))%105deltalambdaU1(k)=(h01+0.58*h11)/b01*(cs1(k)/(cs1(k)+1.5*b01)) lambdas1st(k)=KU1*(lambdaU1-deltalambdaU1(k))+KL1*lambdaL1 %106Xs1sts(k)=lambdas1st(k)/lambdas1*Xs1s%107Xdelta1sts(k)=KZ(k)*Xdelta1s%108Xsigma1sts(k)=Xs1sts(k)+Xdelta1sts(k)+XE1s%109hB=h22;xi=1.987*10^(-3)*hB*(f/rhow1)^(0.5)%110b1b2=b22/b12;KF(k)=input('输入电阻增加系数KF(k)=')Kx(k)=input('输入电抗增加系数Kx(k)=')%111deltalambdaU2(k)=h02/b02*(cs2(k)/(cs2(k)+b02))lambdas2st(k)=(lambdaU2-deltalambdaU2(k))+Kx(k)*lambdaL2%112Xs2sts(k)=lambdas2st(k)*Xs2s/lambdas2%113Xdelta2sts(k)=KZ(k)*Xdelta2s%114Xksts(k)=KZ(k)*Xsks%115Xsigma2sts(k)=(Xs2sts(k)+Xdelta2sts(k)+Xksts(k)+XE2s)*Cx%116Xsigmasts(k)=Xsigma1sts(k)+Xsigma2sts(k)%117R2sts(k)=KF(k)*RBs+RRs%118Rsts(k)=R1s+R2sts(k)Zsts(k)=(Rsts(k)^2+(Xsigmasts(k))^2)^(0.5)%120Ist(2*k)=IKW/Zsts(k)precisionist=abs((Ist(2*k)-Ist(2*k-1))/Ist(2*k))Ist(2*k+1)=input('输入Ist(2*k+1)=')endend%121ist=Ist(6)/I1(2)%122Tsts=R2sts(3)/Zsts(3)^2*(1-sN)。