17 开关磁阻电机设计及控制.doc

毕业设计报告书课题名称：开关磁阻电机设计及控制专业班级：电机制造及运行维护报 告 人：指导老师：报告日期：年月日毕业设计(论文)任务书课题名称(包括主要技术参数)及要求:题 目: 开关磁阻电机设计及控制功 率：电源电压：相电压额定转速：额定效率：调速范围： 运行方式：连续运行绝缘等级: 课题内容及工作量:1. 了解开关磁阻电机原理、结构及控制2. 确定电机主要参数、主要尺寸3. 电磁方案计算4. 控制研究5. 编写设计说明书目 录前言第一章 开关型磁阻电机驱动系统的组成1.1 SRD基本结构1.2 SRD的工作原理1.3 SRD的特点1.4 SRD的研究动态第二章 SR电机工作原理和特点2.1 基本结构2.2 工作原理2.3 SR电机的特点第三章 开关磁阻电动机电磁设计3.1 设计特点3.2 确定SR电动机额定数据3.3 计算SR电动机主要尺寸3.4 计算SR电动机其他尺寸3.5 选择导通角和磁负荷，计算绕组匝数3.6 核算SR电动机的额定转矩、额定电流3.7 计算槽满率、导线线径及每相绕组的电阻3.8 参数计算和校核额定效率3.9 计算机辅助设计程序第四章 SRD的基本控制原理4.1 基本控制策略4.2 正反转控制4.3 SR电动机的起动运行4.4 低速斩波控制4.5 高速单脉冲控制第五章 SRD系统控制原理设计5.1 控制系统的总体设计5.2 控制器的设计5.3 电流检测器的设计5.4 位置检测器的设计5.5 功率变换器的设计第六章 SR电动机运行时的转矩脉动与噪声6.1 SR电动机的转矩脉动分析6.2 SR电动机的噪声分析第七章 开关磁阻电动机应用概况7.1 SRD系统在电动车方面的应用7.2 SRD系统在家用电器方面的应用第八章 结束语第九章 参考资料附录：参观实习报告摘 要：在此次设计中，主要是系统的介绍开关磁阻电动机驱动系统的基本结构和基本工作原理，并依此为基础简要的进行开关磁阻电动机驱动系统的设计；重点分析和详细的说明了开关磁阻电动机驱动系统的执行元件开关磁阻电动机的结构、控制原理和控制策略，并根据给定的电动机的额定数据进行电机的主要参数、主要尺寸、电磁方案的计算、性能的核算和控制方案研究关键词：开关磁阻电动机驱动系统 开关磁阻电动机 电磁设计 控制方案研究 转矩的脉动与噪声前言磁阻式电动机诞生于160年前，但在此后漫长时期内，它一直被认为是一种性能（效率、功率因数、利用系数等不高的电动机，故仅应用于少数小功率场所。70年代初，美国福特电动机（Ford Motor）公司研制出最早的开关磁阻电动机调速系统（Switched Reluctance Drive）。70年代中期，英国里兹（Leeds）大学和诺丁汉（Nottingham）大学共同研制以电动车辆为目标的开关磁阻电机调速系统。80年代以后，开关磁阻电机调速系统的性能指标高于其他传统的传动系统。因以结构简单、运作可靠及效率高等突出特点，成为交流电机调速系统、直流电机调速系统和无刷直流电机调速系统有力的竞争者。目前，开关磁阻电机已广泛应用于纺织、冶金、机械、石油、化工、航空、交通运输等行业。美国、英国、韩国已在家用电器中进行了应用，并取得了良好的效果。开关磁阻电动机是磁阻电机与电子开关驱动控制器组成一体的控制装置，又称开关磁阻电动机调速系统。它融新型电动机结构与现代电力电子技术、控制技术一体，兼有交流变频调速系统的电动机结构简单、坚固耐用、无刷、无整流子和直流调速系统的调速性能好、控制电路简单、输出特性变化灵活、能频繁正反转和起、制动、在负载和转速变化时均有较高效率以及能适应恶劣环境等优点，是一种性能价格比较高的无级调速系统。然而，作为SRD主体的磁阻电机，尽管电磁原理和结构相当简单，但磁路周期性变化并存在严重的局部饱和，其设计和性能分析与传统的电机差异较大，不可能完全按照传统的设计方法和理论进行该电机的电磁性能设计，因此在设计中，应给予足够的重视。在电机的控制方面，若对SR电机运行特性和控制方式缺乏了解，亦难以制定合理的控制策略，而控制策略的合理与不合理，关系到电机是否能达不到期望的运行性能，因此，必须在正确的理论指导下，设计出合理的控制策略。 由于这是第一次搞设计，设计经验不足，知识结构不够完备，在最大可能占有相关资料的前提下，完成了此次设计任务，因此存在许多不足之处。可能本设计结果和工厂的实际生产设计有所不同，但本次设计为以后到工厂搞设计积累了宝贵的经验。在这一个多月的设计中，由电气系老师谢卫才担任指导，在设计过程中谢老师给予了悉心的指导和热情的帮助，在此对谢老师谨致诚挚的谢意。本次设计难免有不少错误的地方，敬请各位老师提出加以指正。由于开关磁阻电动机是一个系统，而作为它的执行元件SR电动机其本身无法单独工作，因此，本次设计是从开关型磁阻电动机驱动系统出发。第一章 开关型磁阻电机驱动系统的组成1.1 基本结构开关型磁阻电机驱动系统（英文全称为Switched Reluctance Drive缩写为SRD）主要由开关型磁阻电机（英文全称为Swithd Reluctance Motors缩写为SR电动机）、功率变换器、控制器和位置检测器四部分组成，其基本框图如图1所示。其中SR电动机和位置检测器组成一体，功率变换器、控制器和电源是跟SR电动机相对独立的电子电气装置。 SR电动机是整个系统的执行部件，它采用双凸极结构，磁路的磁阻变化较大。定、转子极数不同，转子仅由叠片叠压而成，既无绕组也无永磁体；定子各极上绕有集中绕组，径向相对极的绕组串联，构成一相。SR电动机遵循“磁阻最小原理”，产生的磁拉力形成磁阻性质的电磁转矩而使转子转动。功率变换器向SR电动机提供运转所需要的能量，其本身由蓄电池或交流电源整流后得到的直流电供电。它主要由功率开关和续流等电路组成，实现电能的转换和能量回馈。它的主电路的结构形式与供电电压、电动机相数及主开关器件种类等有关。位置检测器是转子位置及速度等信号的提供者，是用来确定转子磁极的相对位置。控制器必须借助位置检测器获得的转子位置信息，以保证在恰当时刻接通或断开相应的相绕组，以产生不同的电磁转矩，获得不同的转速、转向及运行状态。控制器是系统的中枢，它综合处理速度指令、速度返馈信号及电流传感器、位置传感器的反馈信号，从而控制功率变换器中主开关器件的工作状态，实现对SR电动机运行状态的最终控制。1.2 SRD的工作原理由于不同转子位置处，其相电感L（、i）及L/是不同的，控制相绕组通电时刻即可改变电流的大小及波形，因此产生不同的转速、转向及进入不同的运行状态。SR电机电动机的控制器，根据位置检测器检测到定、转子间的相对位置，按照给定的运行指令导通相应的相绕组的度开关元件。相绕组中有电流流过的产生磁通，转子因受到使其向磁阻最小位置运动的趋而产生的磁拉力的作用而转动，直转子转到被其吸引的定子磁极相重合时为止。同时，控制器根据新的位置信息亲断当前相，而导通下一相，因而转子继续向下一个平衡位置转动。这样，控制器根据相应的位置信息按一定的控制逻辑，连续不断地导通和关断相应的主开关，就可以产生连续的同方向转矩，使转子在一定的转速下连续运行。再根据一定的控制策略，控制各相绕组的通断时刻以及绕组电流的大小，使系统处在最佳状态下运行。因此，位置闭环是SRD有别于步进驱动系统的重要标志之一。 1.3 SRD的特点开关磁阻电动机驱动系统，综合了交流电动机驱动系统和直流电动机驱动系统的优点，作为一种新型的调速系统，它的电动机本体结构、变换器型式以及控制方式都与众不同。其主要特点有：1 执行元件SR电动机结构简单，制造工序少、成本低、工作可靠、维修量小等优点。2 功率变换器主开关元件少，且每一个开关元件都与电动机绕组串联，从结构上排除了短路故障的可能，不会出现直通现象。3 系统具有良好的动态响应。4 系统的调速范围广、运行效率高、控制参数多、控制方式灵活、易实现正、反转控制及四象限运行。5 系统的噪声大于一般的电动机，且电动机与控制器间的联线较多，这是SRD系统的最度的缺点。1.4 SRD的研究动态性能优良的电气传动都需要调速电动机，先进的机电一体化调速电动机，正是顺应了这种必然发展趋势的需要。开关磁阻电动机调速系统作为一种方兴未艾的新技术，涉及到电机、微电子、电力电子、计算机、自动控制和机械工程等众多学科领域。从目前发展水平看，无论是在理论上还是在应用上，仍有许多问题有待进一步研究与解决。1. 进一步完善SR电动机设计理论，SR电动机的非线性使其性能的精确分析和计算较为困难。目前普遍采用二维非线性有限元法分析SR电动机内部的饱和磁场.其局限性主要表现在两方面：一是对以路为基础的设计方法研究不够。二是现有场的方法的精度亦有待提高，应计及端部效应，开展SR电动机三维场的研究。2加强铁心损耗的理论研究。SR电动机磁场特性的非线性导致相绕组供电电压和电流波形较为复杂，一般为单向脉动的非正弦波；定转子各部分铁心中的磁通密度变化规律也不相同，因此对定、转子铁心损耗的计算与测量颇为困难，面临的问题主要是如何建立准确、实用的铁心损耗计算模型和分析、测试手段。3 加强对转矩脉动及噪声的理论研究，SR电动机的转矩脉动及其引起的噪声是SRD一个颇为突出的缺点，这限制了其在诸如伺服驱动这类低速要求平稳且有一定静态转矩保持能力场合下的应用。因此，研究抑制SR电动机的振动和噪声也是改善SRD性能的重要课题之一，减少SR电动机振动和噪声的关键是在于如何减定子上的径向力的大小。从SR电动机自身的结构设计上，主要是合理设计形状、定子磁轭强度和电动机刚度，合理选择气隙、极弧参数；从控制角度看，主要是优选导通角和关断角及调节方案，尽可能调节好各相工作参数的对称性。4 进一步提高SRD中电动机、功率变换器及控制器三者之间的协调设计水平，SRD是典型的机电一体化装置，它的各部分间的联系非常密切，为使系统整体最优，就不应将各部分的设计院割裂开来，必须借助计算机辅助设计将电动机、功率变换器及控制器三者作为一个整体进行优化设计。目前尚停留在仅对特定类型的SRD分析核算程序的水平，只能完成局部的综合设计和个别参数优化。5 如何有效地回收绕组释放的磁场能量。因为SR电机的断开相时，铁心内部有剩磁，而SR电动机的转子在不停的旋转，切割了电磁线，由法拉第电磁感应定律产生了能量，此部分磁场的能量回馈给电源，由于此能量不是连续的，给能量的回收利用造成了很大的困难，主要如何在有效建立回馈方式。第二章 SR电机工作原理和特点SR电动机是SRD驱动系统的执行元件。它与传统的磁阻电动机相比，具有本质的区别，这是磁阻电动机和电力电子开关电路相结合产生的一种节能型变速电动机。它的结构和工作原理与传统的交直流电动机有着根本的区别。其主要表现在这种电动机不像传统电动机那样，依靠定、转子绕组电流产生的磁场的相互作用形成电磁转矩，而是依靠“磁阻最小原理”产生磁阻性质的转矩。因此，这种电机的原理和控制方法上与其他的电机肯定有不同之处。作为一种新型的电动机，在理论分析与实际应用方面还需要做更多、更细致的研究工作。2.1 基本结构图1所示是开关磁阻电机的典型结构原理图。从原理图上可以看出该电动机在结构上与步进电动机相似，其运行原理遵循“磁阻最小原理”磁通总是要沿磁阻最小的路径闭合，因磁场被扭曲而产生切向磁拉力的磁阻性质的电磁转矩。因此，它的结构原则是：转子旋转时磁路的磁阻要尽可能大的变化。所以SR电动机均采用双凸极结构。 图1 SR电动机结构原理图2.2 工作原理其运行原理的具体过程如下：当A相绕组电流控制开关S1、S2闭合时，A相励磁，而其他相不通电流，这是，电动机内建立以AA为轴线的磁场，磁场经过定子轭、定子磁极、气隙、转子磁极、转子轭等处闭合。由于通电时定、转子磁极轴线不重合，因此，此时通过气隙的磁力线是弯曲的。按照法拉第电磁力的观点，每根磁力线都被磁力线连着的定、转子沿磁力线方向相互吸引，因而产生的切向磁场力图使转子旋转到转子极轴线aa与定子极轴线AA的重合位置，产生磁阻性质的电磁转矩，使转子转动起来。当转子aa与AA轴线重合时，转子已达到稳定平衡位置，即A相定、转子磁极正对时，切向磁拉力消失，转子不再转动，此时在顺序给BCD相绕组通电（B、C、D各相绕组在图中未画出），则转子便按逆时针方向连续转动起来；反之，依次给BACD相绕组通电，则转子会沿顺时针方向转动.另外，当主开关管S1、S2导通时,A相绕组从直流电源V吸收电能;而当S1、S2关断时,绕组电流通过续流二极管D1、D2，将剩余能量回馈给电源V.因此，开关磁阻电机具有再生的能力，系统效率高。2.3 SR电机的特点对开关磁阻电机的理论研究和实践证明，该系统具有许多显著的特点： （1）电机结构简单、坚固，制造工艺简单，成本低，可工作于极高转速；定子线圈嵌放容易，端部短而牢固，工作可靠，能适用于各种恶劣、高温甚至强振动环境。 （2）损耗主要产生在定子，电机易于冷却；转子无永磁体，可允许有较高的温升.能四象限运行，具有较强的再生制动能力。 （3）转矩方向与电流方向无关，只要控制主开关器件的导通关断角，即可改变电动机的工作状态,无须辅助电力电子开关器件,降低系统成本。 （4）功率变换器不会出现直通故障，可靠性高。 （5）起动转矩大，低速性能好，无感应电动机在起动时所出现的冲击电流现象。在宽广的转速和功率范围内都具有高效率 （6）调速范围宽，控制灵活，易于实现各种特殊要求的转矩速度特性。 与当前广泛应用的变频调速感应电动机相比，开关磁阻电机在成本、效率、调速性能、单位体积功率、可靠性、散热性等都具有明显的优势或竞争力，图1-1是对几种常用变速传动系统各项主要经济指标所作的比较。 系统类型比较项目电 磁调速系统直 流调速系统PWM变频调速系统开关磁阻调速系统成 本 0.8 1.01.5 1.0 额定转速时效率（%）757677831/2额定转速效率（%）38656580电动机容量/体积 0.81.0 0.9 1.0控制能力0.31.0 0.5 0.9控制电路复杂性0.21.01.81.2 可靠性1.3 1.00.91.1 图1-1几种常用变速传动系统各项主经济指标所作的比较开关磁阻电机在宽广的速度和功率范围内都能保持较高的效率，这是变频调速感应电动机难以比拟的。感应电动机要取得与直流电机相近的调速特性需采用复杂的矢量控制系统，而开关磁阻电机通过调整开通角、关断角、电压和电流，可以得到不同负载要求的机械特性，控制简单、灵活，能容易地实现软启动和四象限运行，而且由于这是一种纯逻辑的控制方式，很容易智能化，通过修改软件调整电机工作特性满足不同应用要求。由于开关磁阻电机固有的转矩波动，可能导致较大的噪声和振动，事实上这种情况的发生往往与电机设计和控制的不合理相关，通过优化电机设计和控制策略，转矩波动和噪声完全可以得到有效的抑制. 与变频调速的另一类电机永磁无刷直流电动机相比，开关磁阻电机在成本、调速性能、可靠性和散热性等依然保持优势，但和变频感应电动机一样都存在励磁惩罚（既励磁电损耗）的问题，效率和单位体积功率稍显不足。但在一些应用场合，开关磁阻电机在效率和输出功率依然具有优势。第三 章 开关磁阻电动机电磁设计3.1设计特点本次设计在SR电动机的输出方程为基础上完成SR电动机结构尺寸、绕组匝数确定的，由于输出方程引入了电流系数，使之应用更可靠实用。与传统的电机设计相比，SR电机存在以下显著特点：（1）即使在结构尺寸、绕组匝数、端电压、线规确定之后，其电磁转矩还与控制参数相关。（2）电磁转矩由-轨迹包围的磁共能面积决定，磁化曲线族由不同转子位置的磁化曲线组成，磁路不均匀，计算工作量大、困难。（3）绕组匝数对最小电感值影响较大，影响电流峰值和电磁转矩的计算。（4）存在恒转矩区、恒功率区的调速范围，（5）效率计算困难。因此，在设计中采用了计算机辅助设计方法，完成SR电动机在额定点的电磁设计、效率计算和参数计算。3.2 确定SR电动机额定数据1功率：2电源电压：相电压3额定转速：4额定效率：5调速范围：恒转矩，恒功率6运行方式：连续运行7绝缘等级: 3.3 计算SR电动机的主要尺寸1相数 2定子极数3转子极数4绕组端电压5电磁功率 6 细长比7 电磁负荷初选8 转子外径 9 铁心叠长 10定子外径 为了利用Y180L6机座，取3.4 计算SR电动机其他尺寸1气隙2定、转子极弧、注：以下只列出各步骤的计算公式，而计算结果及说明在计算机辅助设计程序“主要尺寸计算”3定、转子极宽、4第二气隙5定、转子轭高、6轴径7定子槽深8有效铁心长度9转子极距根据以上结果，绘制出电机的结构尺寸说明图如下3.5 选择导通角和磁负荷，计算绕组匝数1控制参数 额定转速预取开通角关断角导通角2每相绕组串联匝数，取，重新校核注：以下各步骤的计算结果在计算机程序“电机校核程序”3.6 核算SR电动机的额定转矩、额定电流、额定效率1电流及转矩的计算2额定电磁转矩，是磁化曲线的面积3绕组电流有效值 4实际电流峰值注：以各步骤的结果在计算机程序“参数计算”。3.7 选择槽满率，计算导线线径及每相绕组的电阻1定子极间窗口面积2导线规格、导线面积、槽满率、预取电流密度由GB1193导线规格表查得3每槽导体净面积槽满率4电流密度5绕组平均匝长6每相绕组导线总长 7每相绕组电阻3.8 参数计算1铜重2定子铁心体积3转子铁心体积4铁心重5电负荷6铜耗 7。电机利用系数K 8电机的效率3.9 计算机辅助设计程序SR电机的计算机辅助设计，主要包括五个部分：一是主要尺寸的计算；二是a(最大电感)磁化曲线计算；三是u(最小电感)磁化曲线计算；四是电机校核程序；五是电机参数计算。/*主要尺寸的计算*/#include #include #define PI 3.1415main() double bps,bpr,gi,hcs,hcr,beta_sf,beta_rf,Hai; double Di,ds,Lfe,cta_c,tr,B_delta,Bps,Nph,Sw; double Nr=4,Ns=6,Da=132,Ds=260,La=155,N=1500; double g=0.4,beta_s=30,beta_r=32,U=280;printf(Bps=?); scanf(%lf,&Bps); beta_sf=beta_s*PI/180.0; beta_rf=beta_r*PI/180.0; bps=(Da+2*g)*sin(beta_sf/2.0); bpr=Da*sin(beta_rf/2.0); gi=bps/2.0; hcs=1.3*bps/2.0; hcr=1.4*bpr/2.0; Di=Da-2*(gi-g)-2*hcr; ds=(Ds-Da-2*g-2*hcs)/2.0; Lfe=0.95*La; tr=PI*Da/Nr; B_delta=0.805*bps*Bps/tr; cta_c=(2*180.0/Nr-beta_r)/2.0*(PI/180.0); Nph=(3.04*Nr*U*cta_c)/(N*B_delta*Da*0.001*1.05*La*0.001); Hai=(Ds/2.0-hcs)*(Ds/2.0-hcs)-(Da/2.0+g)*(Da/2.0+g); Sw=PI*Hai/2.0/Ns-bps*ds/2.0; printf(bps=%lftt bpr=%lftt gi=%lfttn,bps,bpr,gi); printf(hcs=%lftt hcr=%lftt Di=%lfttn,hcs,hcr,Di); printf(Lfe=%lftt ds=%lftt tr=%lfttn,Lfe,ds,tr); printf(B_delta=%lft Nph=%lftt Sw=%lftn,B_delta,Nph,Sw); Bps=?1.6bps=34.370179 bpr=36.383086 gi=17.185089hcs=22.340616 hcr=25.468160 Di=47.493501Lfe=147.250000 ds=41.259384 tr=103.669500B_delta=0.427018 Nph=125.233733 Sw=1171.032364根据以上结果，取定电机尺寸值(单位：mm)参 数尺 寸参 数尺 寸定子外径Ds260定子极弧s30转子外径Da132转子极弧r32铁心叠长La155定子轭高hcs22.34气隙g0.4转子轭高hcr25.46第二气隙gi17.4轴径Di47匝数130匝转子极距tr103.7定子槽深ds41.25校核B0.427T窗口面积Sw1171有效长度Lef147.2/*a(最大电感)磁化曲线计算*/#include #include #define PI 3.1415#define mu0 1.25e-6/*磁场强度H开始一元*/double LAG(double xx) double A,B,C; int i,j; double x22=0.4,0.42,0.45,0.5,0.54,0.57,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.85,1.9,2.1; double y22=1.37,1.40,1.46,1.56,1.64,1.70,1.77,2.03,2.39,2.91,3.63,4.60,6.10,8.14,11.4,17.2,30.7,61.5,110,141,180,260; C=0; for(j=0;j22;j+) for(i=0,A=1;i10.6|Fai3)printf(3Fai10.6);else beta_sf=PI*beta_s/180.0;beta_rf=PI*beta_r/180.0;lfe=la*0.95;bps=2.0*(Da/2.0+g)*sin(beta_sf/2.0);Sts=bps*lfe;bpr=2.0*Da/2.0*sin(beta_rf/2.0);Str=bpr*lfe;Scs=hcs*lfe;Scr=hcr*lfe;Sg=(beta_sf+beta_rf)*la*(Da+g)/2.0/2.0/2.0; /*气隙的计算*/lts=2.0*(Ds/2.0-hcs-Da/2.0-g);ltr=2.0*(gi-g);lg=2.0*g;lcs=(Ds-hcs)/2.0*PI;lcr=(Da/2.0-(gi-g)-hcr/2.0)*PI;fai_2=Fai/1000.0;Bts=fai_2/(Sts*0.000001);Btr=fai_2/(Str*0.000001)/1.2;Bg=fai_2/(Sg*0.000001)/2.25;Bcs=fai_2/(0.000001*Scs)/1.7;Bcr=fai_2/(0.000001*Scr)/2.0;a1=LAG(Bts)*lts/10.0;a2=LAG(Btr)*ltr/10.0;a3=LAG(Bcs)*lcs/10.0;a4=LAG(Bcr)*lcr/10.0;sgma_F=a1+a2+Bg*(lg/1000.0)/mu0+a3+a4;is_a=sgma_F/Nph;printf(is_a=%lft,is_a);运行结果如下：Fai=?6.0is_a=6.626409 Fai=?7.0is_a=8.622435 Fai=?8.0is_a=12.254323 Fai=?8.2is_a=13.569168 Fai=?8.5is_a=16.218487 Fai=?8.7is_a=18.593447 Fai=?9.0is_a=23.540281 Fai=?9.2is_a=27.993731 Fai=?9.4is_a=33.243781 Fai=?9.5is_a=36.146648 Fai=?9.6is_a=39.468180 Fai=?9.7is_a=43.855576 Fai=?9.8is_a=50.709707 Fai=?10结果说明：is_a是磁化电流。单位AFai是磁化电流所对应的磁通量。单位mWb注：以上数据是提供给面积计算。/*u(最小电感)磁化曲线计算*/#include #include #define PI 3.1415#define mu0 1.25e-6/*u齐位置电感计算*/double cta_u_Lu(double Nph) double tr,ts,P,W,V,Y,H,h,U,M,N,n,gF,Lu; double p1,p2,p3,p4,p5,pj,beta_sf,beta_rf,beta_s=30,beta_r=32; double D1=260,D2=132,La=155,Lfe=147.25,g=0.4,gi=17.4,hj1=22.34; double L_F,D,G,Gj,r,C,L_2D,xi,h1,h2,Nr=4,Ns=6; tr=2.0*PI/Nr; ts=2.0*PI/Ns; beta_sf=PI*beta_s/180.0; beta_rf=PI*beta_r/180.0; G=g/D1; Gj=gi/D1; D=D2/D1; C=hj1/D1; P=(D/2.0+G)*sin(beta_sf/2.0); W=(D/2.0+G)*tan(ts/2.0)-P; h1=sqrt(1.0/2.0-C)*(1.0/2.0-C)-P*P); h2=sqrt(1.0/2.0-C)*(1.0/2.0-C)-(W+P)*(W+P); V=h1-(D/2.0+G); Y=D/2.0+G-D/2.0*cos(tr-beta_rf)/2.0); H=D/2.0*sin(tr-beta_rf)/2.0)-P; h=H*D1; U=h1-h2; r=atan(W/U); M=U+V/2.0; N=Y+V/2.0; n=N*D1; p1=r*M*M*M*M/(4.0*W*(2.0*V+U)*W*(2.0*V+U); p2=(P+H)/Gj-2.0/PI*(tr-beta_rf); p3=2.0/PI*(log(N/H)+2.0*(N-H)*Y/W/V-(N*N-H*H)/(4.0*W*W*V*V)*(PI*W*V-2.0*Y*Y)-(N*N*N-H*H*H)*Y*PI/(6.0*W*W*V*V)+(N*N*N*N-H*H*H*H)*PI*PI/(64.0*W*W*V*V); p4=2.0/(tr-beta_rf)*log(2.0*tan(tr-beta_rf)+PI-(tr-beta_rf)/(2.0*tan(tr-beta_rf)+PI-2.0*(tr-beta_rf); p5=2.0/(PI-(tr-beta_rf)*log(2.0*Gj/H/(PI-(tr-beta_rf);pj=p1+p2+p3+p4+p5; gF=(PI/2.0*n+PI/2.0*h+(PI/2.0-(tr-beta_rf)/2.0)*h+2.0*gi)/5.0;xi=2.0/PI*(atan(2.0*n/gF)-gF/4.0/n*log(1.0+(2.0*n/gF)*(2.0*n/gF); L_F=La+2.0*n*(1-xi); L_2D=mu0*pj*Lfe*Nph*Nph*0.001; Lu=L_2D*(2.0*L_F/La-1.0); return Lu;main()double s,Fai,is_a,Nph0; printf(is_a=?); scanf(%lf,&is_a); s=is_a*cta_u_Lu(Nph); Fai=s/Nph; printf(Fai=%lft,Fai);运行结果如下：is_a=?6.626409Fai=0.000438 is_a=?8.622435Fai=0.000570 is_a=?12.254323Fai=0.000811 is_a=?13.569168Fai=0.000898 is_a=?16.218487Fai=0.001073 is_a=?18.593557Fai=0.001230 is_a=?23.540281Fai=0.001557 is_a=?27.993731Fai=0.001852 is_a=?33.243781Fai=0.002199 is_a=?36.146648Fai=0.002391 is_a=?39.468180Fai=0.002611 is_a=?43.855576Fai=0.002901 is_a=?50.709707Fai=0.003354结果说明：Fai是磁通量。单位mWb.is_a是磁通量所对应的磁化电流。单位A.注：以上数据是提供给面积计算。 /*电机校核程序*/#include #include #define PI 3.1415# define mu0 1.25e-6# define RO 0.0217/*磁化曲线面积的计算*/double W1(double b)double s;int c=4,i;double y15=0,0,0.006,0.007,0.008,0.0082,0.0087,0.009,0.0092,0.0093,0.0094,0.0095,0.0096,0.0097,0.0098;s=(y1-y14)/2.0; for(i=1;i8;i+) s=s+2.0*y2*i-1+y2*i; s=s*b/(3.0*c); return s;/*磁化曲线面积计算*/double W2(double b)double s;int c=7,i;double y15=0,0,0.000438,0.00057,0.000811,0.000898,0.001073,0.00123,0.001557,0.001852,0.002199,0.002391,0.002611,0.002901,0.003554;s=(y1-y14)/2.0; for(i=1;i=gi) gh=(Da/2.0+g)*cos(beta_sf/6.0)-(Da/2.0+g)*cos(beta_sf/2.0)+g;elsegh=(Da/2.0+g)*cos(beta_sf/2.0-beta_sf/7.0)-(Da/2.0+g)*cos(beta_sf/2.0)+g;w=(Da/2.0+g)*tan(PI/Ns)-bps/2.0;alph=PI/2.0-(asin(w+bps/2.0)/(Ds/2.0-hcs)+asin(bps/2.0)/(Ds/2.0-hcs)/2.0;u=w/tan(alph);v=sqrt(Ds/2.0-hcs)*(Ds/2.0-hcs)-(Da/2.0+g)*(Da/2.0+g)*tan(PI/Ns)* tan(PI/Ns)-(Da/2.0+g);hms=Ds/2.0-hcs-Da/2.0-g;t1=(gi-gh)/beta1;t2=h/sin(theta)-gi;t3=(gi-f)/beta2;m1=u+v/2.0;m2=u+2.0/3.0*v;n1=hms-m1;delta=PI/2.0+beta_sf/2.0-acos(n1/Da);p1=alph*m1*m1*m1*m1/(4.0*w*(2.0*v+u)*w*(2.0*v+u);p2_1=1.0/(PI/2.0+delta)*log(gh+(PI/2.0+delta)*n1)/gh);p2_2_1=n1*n1/(2.0*(PI/2.0+delta)-gh*n1/(PI/2.0+delta)*(PI/2.0+delta);p2_2_2=gh*gh/(PI/2.0+delta)*(PI/2.0+delta)*(PI/2.0+delta)*log(PI/2.0+delta)*n1+gh)/gh);p2_2=PI/(2.0*w*v)*(p2_2_1+p2_2_2);p2_3_1=n1*n1*n1/(3.0*(PI/2.0+delta)-gh*n1*n1/(2.0*(PI/2.0+delta)*(PI/2.0+delta);p2_3_2=gh*gh*n1/(PI/2.0+delta)*(PI/2.0+delta)*(PI/2.0+delta);p2_3_3=gh*gh/(PI/2.0+delta)*(PI/2.0+delta)*(PI/2.0+delta)*(PI/2.0+delta);p2_3_4=log(PI/2.0+delta)*n1+gh)/gh);p2_3_5=PI*PI/(16.0*w*v*v)*(n1*n1*n1*n1/(4.0*(PI/2.0+delta);p2_3=p2_3_5-PI*PI/(16.0*w*v*v)*gh/(PI/2.0+delta)*(p2_3_1+p2_3_2-p2_3_3*p2_3_4);p2=p2_1-p2_2+p2_3;p3=1.0/beta1*log(beta1*t1+gh)/gh);if(f=gi)b1=sqrt(gi+t2)*(gi+t2)-gi*gi); S5=(w*v-t2*b1/4.0)/(w*v); b2=sqrt(gi+hms-m2)*(gi+hms-m2)-gi*gi); S6=(w*v-v*b2/12.0-t2*b1/4.0)/(w*v); p4=(bps-t1)/gi; p5=S5*S5/(theta-ga_ma)*log(theta*(gi+t2)/(theta*gi+ga_ma*t2); p6=S6*S6/theta*log(theta*(gi+hms-m2)/(theta*(gi+t2); gF_1=(PI/2.0*n1+gh+delta*n1)+gh+(gh+beta1*t1)+2.0*gi+(gi*theta+ga_ma*t2); gF=1.0/8.0*(gF_1+(gi+t2)*theta+(gi+hms-m2)*theta);else S6p=1.0-v/(36.0*w*tan(PI/Ns); p4=(bps-t1-t3)/gi; p5=1.0/beta2*log(f+beta2*t3)/f); p6=2.0/PI*S6p*S6p*log(f+PI/2.0*(hms-m2)/f);gF_1=(PI/2.0*n1+gh+delta*n1)+gh+(gh+beta1*t1)+2.0*gi+(f+beta2*t3); gF=1.0/8.0*(gF_1+f+(f+PI/2.0*(hms-m2); p7=alph*m2*m2*m2*m2/(4.0*w*(2.0*v+u)*w*(2.0*v+u); sgma_p=p1+p2+p3+p4+p5+p6+p7; n=n1+gh;xi=2.0/PI*(atan(2.0*n/gF)-gF/4.0/n*log(1.0+(2.0*n/gF)*(2.0*n/gF); L_F=La+2*n*(1.0-xi); L_2D=2.0*mu0*Lfe*(Nph/2.0)*(Nph/2.0)*sgma_p*0.001; L1=L_2D*(2.0*L_F/La-1.0); return L1;/*u齐位置电感计算*/double cta_u_Lu(double Nph) double tr,ts,P,W,V,Y,H,h,U,M,N,n,gF,Lu; double p1,p2,p3,p4,p5,pj,beta_sf,beta_rf,beta_s=30,beta_r=32; double D1=260,D2=132,La=155,Lfe=147.25,g=0.4,gi=17.4,Nr=4,Ns=6; double L_F,D,G,Gj,r,C,L_2D,xi,h1,h2,hj1=22.34; tr=2.0*PI/Nr; ts=2.0*PI/Ns; beta_sf=PI*be