ANSYS低频电磁场的场路耦合分析应用实例.doc

【精品文档】如有侵权，请联系网站删除，仅供学习与交流ANSYS低频电磁场的场路耦合分析应用实例.精品文档.ANSYS低频电磁场的场路耦合分析功能及 感应电机场路耦合分析应用实例ANSYS程序的电磁场分析功能分为低频和高频这两个主要的分析模块，在低频电磁场分析部分，其功能涵盖静态磁场分析、谐波（交流）磁场分析、瞬态磁场分析、静电场分析、电流传导分析、电路分析、场路耦合分析等内容。本文重点介绍场路耦合分析的一些应用技巧，并以一个两极三相感应电机启动状态的场路耦合分析为例予以简要说明。对于其它分析类型，可参考ANSYS相应的技术说明手册。1、 择合适的单元对于场路耦合分析，场分析模型中有两种类型的终端条件实现与电路的联接：电路供电绞线圈和电路供电块导体。在二维状态下，场分析中与电路联接的部分必须用53号单元，而三维时必须用97号单元，这两种单元的第一关键选项设置为3时，表示电路供电绞线圈，其第一关键选项设置为4时，表示电路供电块导体。在路分析模型中的124号电路单元里，也有两种类型的单元选项与上面两种形式对应，以实现“联接”，当其第一关键选项设置为5时，为二维或三维绞线圈单元，设置为6时，为二维块导体单元，设置为7时，为三维块导体单元。当这些电路单元的类型设置与场单元的一致、并共享某些节点后，即实现了联接，如下所述。2、 实现联接对于二维场路耦合分析，场分析模型中的每一个电路耦合绞线圈（53号单元的Keyopt(1)=3）和电路耦合块导体（53号单元的Keyopt(1)=4）上，必须分别耦合所有节点的电流自由度（CURR）和电动势自由度（EMF）。在定义电路中的绞线圈单元（124号单元的Keyopt(1)=5）或块导体单元（124号单元的Keyopt(1)=6）时（这两种单元都只有I、J、K三个节点），其K节点必须为场模型中相应绞线圈或相应块导体上的一个节点。对于三维场路耦合分析，场分析模型中的每一个电路耦合绞线圈（97号单元的Keyopt(1)=3）上，必须耦合所有节点的电流自由度（CURR）和电动势自由度（EMF）。在定义电路中的绞线圈单元（124号单元的Keyopt(1)=5）时（此时该单元共有I、J、K三个节点），其K节点必须为场模型中相应绞线圈上的一个节点。三维块导体的定义和联接要稍微麻烦些，首先，场模型中的块导体上，其两个与路联接的端面需分别设置场路耦合界面标志（MCI分别等于+1和-1）；其次，在定义电路中的块导体单元（124号单元的Keyopt(1)=7）时（此时该单元共有I、J、K、L四个节点），其K节点必须为前面MCI=-1的那个面上的一个节点，其L节点必须为前面MCI=+1的那个面上的一个节点；最后，再进行一系列自由度耦合操作：一是耦合电路单元的I节点和K节点所在面上所有节点的电压自由度（VOLT），二是耦合L节点所在面上所有节点的电压自由度（VOLT）（注意此处不含电路单元的J节点），三是分别耦合K节点和L节点所在面上所有节点的电流自由度（CURR）。3、 其它部分的模型建立前面已讲了场路耦合分析中的最重要部分：如何实现场模型与路模型的联接。一旦联接完成，其它部分的模型就可很方便的建立了，比如电阻、电感、电源等其它电路元件，利用ANSYS的专用电路建模程序，直接建立即可。需要指出的是，对于场分析模型而言，坐标的概念很重要，这一点不言而喻；而对于路分析模型而言，其节点或单元坐标没有什么实际意义，重要的是路的联接关系。因此，在定义电路单元的节点时，其坐标值原则上可以任意设定，但为了在图形显示中观察起来更方便，定义电路单元的节点坐标可遵循一定的规则，比如，对于二维场路耦合分析，电路绞线圈单元的I和J节点的X和Y坐标可以设置为其K节点的X和Y坐标（K节点是场模型上的一个节点），而I、J节点的Z坐标可分别设为+z和-z（此处，z为一任意值），这样，场路模型的耦合关系在图形显示上就非常直观了，余者类推。这样的模型建立方式通过ANSYS的APDL语言可以很方便地完成。4、 两极三相感应电机场路耦合分析应用实例在电机设计中，有大量的经典设计程序可以利用，这些经典程序都是基于磁路的方法，但由于非线性（铁芯饱和）及结构复杂性等因素的存在，磁路计算无法得到准确的结果，包括定子槽漏抗、励磁感抗、转子槽漏抗等等，而基于有限元方法的场分析对于这些参数的计算就很有优势。从另外一个方面，电机的端部由于其联线方式的复杂性，用场分析的方式基本上没有办法实现，只有采用路来代替，而路的参数（比如端部阻抗和漏抗等）可用经典电机设计程序计算（由于端部不存在铁芯影响，经典计算很准确）。这种现代计算方法与经典计算方法的结合、场分析方法与路分析方法的结合，双方优势互补，可以获得高精度的计算结果。相应的分析过程、程序设计和使用说明如下! 两极三相感应电机启动状态场路耦合电磁场分析程序! (第一段程序) !本段程序的目的 ! 1 输入电机各参数 ! 2 建立电机几何模型 ! 3 划分场分析部分的网格 ! (场路耦合分析的电路连接部分由后面一段程序完成) ! (单纯的场分析用此段程序即可完成全面的建模任务) !finish/clear/uis,msgpop,3/prep7csys,0btol,1e-10*afun,degmultipro,'start',6*cset, 1, 3,d1, 'Outer Diameter of Stator(mm)',950 !定子外直径*cset, 4, 6,d2, 'Inner Diameter of Stator(mm)',520 !定子内直径*cset, 7, 9,gap, 'Gap between Stator and Rotor(mm)',5.5 !气隙高度*cset,10,12,nslots,'Slot Number of Stator',48 !定子槽数*cset,13,15,nslotr,'Slot Number of Rotor', 40 !转子槽数*cset,16,18,d4, 'Diameter of Rotor Axis(mm)', 280 !转子轴直径multipro,'end'r1=d1/2/1000 ! 定子外半径,转化为国际单位制r2=d2/2/1000 ! 定子内半径,转化为国际单位制gap=gap/1000 ! 定转子间气隙高度,转化为国际单位制r3=r2-gap ! 转子外半径,转化为国际单位制r4=d4/2/1000 ! 转子轴半径,转化为国际单位制delta=gap/1e5 !一个用于确定选择区域的极小量multipro,'start',3*cset, 1, 3,nhole, 'Hole Number of Rotor', 20 !转子轴向孔的数目*cset, 4, 6,dhole, 'Main Diameter of Holes', 395 !转子轴向孔圆心处的直径*cset, 7, 9,dh, 'Hole Diameter',20 !转子轴向孔本身的直径multipro,'end'rhole=dhole/2/1000 ! 转子轴向孔圆心处的半径,转化为国际单位制rh=dh/2/1000 ! 转子轴向孔半径,转化为国际单位制!定义定子槽尺寸参数!(如果采用与前面一样的方式,这些数据也可以由菜单界面输入)ss_h=72/1000 !槽高,转化为国际单位制(下同)ss_w=17/1000 !槽宽ss_ci=2/1000 !绕组层间绝缘厚度ss_turn=48/2 !定子每个绕组导线根数ss_cln=2 !并绕导线根数ss_clw=5.6/1000 !导线宽度ss_clh=2/1000 !导线高度ss_cli=0.5/1000 !导线间漆膜厚度ss_cll=600/1000 !绕组长度(定子铁芯长度)ss_ch=(ss_clh+ss_cli)*(ss_turn/ss_cln) !每个绕组高度ss_cw=(ss_clw+ss_cli)*ss_cln !每个绕组宽度!定义转子槽尺寸参数sr_h1=7/1000 !小缝高度sr_w1=2/1000 !小缝宽度sr_h2=14/1000 !中缝高度sr_w2=7/1000 !中缝宽度sr_h3=30/1000 !槽高sr_w31=17/1000 !槽顶宽sr_w32=12/1000 !槽底宽sr_cll=690/1000 !导条轴向长度(直线段长度)!sr_ci=4.8/2/1000 !槽绝缘厚度!计算定子线圈参数,用于单元实常数定义care1=ss_ch*ss_cw !绕组横截面积turn1=ss_turn/ss_cln !绕组匝数(注意,要除以并绕根数!)leng1=ss_cll !绕组轴向长度fill1=(ss_clw*ss_clh)*ss_turn/care1 !绕组填充系数!计算转子导条参数,用于单元实常数定义hh=sr_h3-2*sr_ci !梯形导条高度ww1=sr_w31-2*sr_ci !梯形导条上底宽ww2=sr_w32-2*sr_ci !梯形导条下底宽care2=(ww1+ww2)*hh/2 !梯形导条横截面积leng2=sr_cll !导条轴向长度!定义单元类型et,1,53,0 !AZ自由度,用于空气和绝缘区域et,2,53,3 !AZ,CURR,EMF自由度,电路耦合绞线圈,用于定子线圈导体et,3,53,4 !AZ,CURR,EMF自由度,电路耦合块导体,用于转子槽内导条et,4,53,0 !AZ自由度,用于铁芯区域r,2,care1,turn1,leng1,fill1 !定子线圈单元的实常数r,3,care2,leng2, !转子导条单元的实常数!定义材料参数mp,murx,1,1 !空气导磁率,用于空气和绝缘mp,murx,2,1 !导体导磁率,用于定子线圈mp,rsvx,2,0.0434e-6 !导体电阻率mp,murx,3,1 !导体导磁率,用于转子导条mp,rsvx,3,0.0434e-6 !导体电阻率!mp,murx,4,2000 !调试程序时,用线性材料,省时间些tb,bh,4,18 !铁芯,输入BH曲线tbpt,31.85,0.1tbpt,45.28,0.2tbpt,55.73,0.3tbpt,64.49,0.4tbpt,72.45,0.5tbpt,76.83,0.6tbpt,79.62,0.7tbpt,87.58,0.8tbpt,104.3,0.9tbpt,121.0,1.0tbpt,143.3,1.1tbpt,189.5,1.2tbpt,262.7,1.3tbpt,429.9,1.4tbpt,995.2,1.5tbpt,2547.8,1.6tbpt,5095.5,1.7tbpt,7643.3,1.78!建立几何模型!(建立几何模型时,充分考虑了如下因素:! 1 快捷(尽量采用简单的(虽然这样稍微增加了编程工作量)! 2 便于划分网格(尤其是定子与转子间的气隙处,无论结构形式如何,! 本程序都可映射生成沿径向均匀分布的层数任意的网格而不影响! 气隙附近的网格密度分布)cyl4,r4 !第一个圆形面(转子轴)cyl4,r4,r3 !第二个圆环面(转子铁芯)cyl4,r2,r1 !第三个圆环面(定子铁芯)! 说明: 由于考虑到电机外壳的漏磁很小,无需建立电机外壳外面的空气单元,! 也无需建立远场吸收单元来考虑电机外壳漏磁! 建立定子几何模型xo=sqrt(r2*2-(ss_w/2)*2)wpave,xoblc4,-ss_w/2,ss_h,ss_wblc4,ss_h-ss_ci,ss_w/2,-ss_ch,-ss_wblc4,ss_h-ss_ch-ss_ci*2,ss_w/2,-ss_ch,-ss_wblc4,ss_h,ss_cw/2,-(ss_ch*2+ss_ci*2),-ss_cwasel,s,loc,x,r2,r2+ss_hcsys,1lsel,s,loc,x,r1-delta,r1+deltaasll,alsel,allaovlap,allasel,s,loc,x,xo,r2adele,all,1asel,s,loc,x,r2,r2+ss_hlslalsel,r,loc,x,r2-delta,r2+deltalcomb,alllsel,s,loc,x,r1-delta,r1+deltaaslllslalsel,r,loc,x,r2-delta,r2+deltalcomb,allallselasel,s,loc,x,r2,r2+ss_h!(下面这种常用的方式运行时间较长,约几分钟)!agen,nslots,all,360/nslots!lsel,s,loc,x,r1-delta,r1+delta!asll,a!aovlap,all!(改用如下方式,虽然程序较长,但运行时间只有几秒!转子建模与此相同)cm,aa1,areaagen,2,allcmsel,u,aa1aadd,allagen,nslots,all,360/nslotscm,aa2,arealsel,s,loc,x,r1-delta,r1+deltaasllcm,aa3,areacmsel,a,aa2allsel,below,areaasba,aa3,aa2cm,aa3,areacmsel,s,aa1agen,nslots,all,360/nslotscmsel,a,aa3allsel,below,areanummrg,kp! 建立转子几何模型csys,0xo=sqrt(r3*2-(sr_w1/2)*2)wpave,xoblc4,sr_w1/2,-sr_h1,-sr_w1xo=xo-sr_h1wpave,xoblc4,sr_w2/2,-sr_h2,-sr_w2asel,s,loc,x,xo-sr_h2,xo+sr_h1asel,r,loc,y,-sr_w2/2,sr_w2/2allsel,below,areaaglue,allx1=xoy1=sr_w2/2z1=0nk=kp(x1,y1,z1)kmodif,nk,x1-(sr_w2/2-sr_w1/2),y1,z1y1=-sr_w2/2nk=kp(x1,y1,z1)kmodif,nk,x1-(sr_w2/2-sr_w1/2),y1,z1xo=xo-sr_h2wpave,xoblc4,sr_w31/2,-sr_h3,-sr_w31x1=xo-sr_h3y1=sr_w31/2z1=0nk=kp(x1,y1,z1)kmodif,nk,x1,y1-(sr_w31/2-sr_w32/2),z1y1=-sr_w31/2nk=kp(x1,y1,z1)kmodif,nk,x1,y1+(sr_w31/2-sr_w32/2),z1asel,a,area,2allsel,below,areaaovlap,allasel,s,loc,x,xo-sr_h3,xo+sr_h2asel,s,loc,x,r3-(sr_h1+sr_h2+sr_h3),r3-sr_h1asel,r,loc,y,-sr_w31/2,sr_w31/2allsel,below,areaaadd,allasel,s,loc,x,r3-(sr_h1+sr_h2+sr_h3),r3asel,r,loc,y,-sr_w31/2,sr_w31/2cm,aa1,arealslacsys,1lsel,r,loc,x,r3-delta,r3+delta*get,nl,line,count*if,nl,eq,2,thenlcomb,all*endifcm,ll1,lineallsel,below,arealsel,s,loc,x,r3-delta,r3+deltacmsel,u,ll1lcomb,allcmsel,s,aa1agen,2,allcmsel,u,aa1aadd,allagen,nslotr,all,360/nslotrcm,aa2,arealsel,s,loc,x,r3-delta,r3+deltaasllcmsel,u,aa1cmsel,u,aa2cm,aa3,areacmsel,a,aa2allsel,below,areaasba,aa3,aa2cm,aa3,areacmsel,s,aa1agen,nslotr,all,360/nslotrcmsel,a,aa3allsel,below,areanummrg,kpwpave,0,0,0!创建气隙圆环面allselnummrg,kpnumcmp,kpnumcmp,linenumcmp,area*get,na,area,countcsys,1lsel,s,loc,x,r3-delta,r3+deltaal,alllsel,s,loc,x,r2-delta,r2+deltaal,allasba,na+2,na+1allsellplot!切割气隙圆环,以便于保证气隙厚度上的网格数csys,1ksel,s,loc,x,r3-delta,r3+delta*get,nk3,kp,count*get,nk3min,kp,num,mincm,kk3,kpksel,s,loc,x,r2-delta,r2+delta*get,nk2,kp,count*get,nk2min,kp,num,mincm,kk2,kpcmsel,a,kk3cm,kk23,kp!(先连接圆环内近似径向的线,以免生成短线)da=gap/(2*3.1416*r3)*360/3 !与气隙尺寸相关的一个小角度lsel,s,loc,x,r2-delta,r2+deltalsel,a,loc,x,r3-delta,r3+delta*get,nl32,line,count*get,nl32min,line,num,mindsys,1*do,i,1,nl32da1=abs(ly(nl32min,0)-ly(nl32min,1)/3 !与最短线相关的一个小角度*if,da1,lt,da,thenda=da1*endifnl32min=lsnext(nl32min)*enddo*do,i,1,nk3cmsel,s,kk3ynk3=ky(nk3min)nk2min1=nk2min*do,j,1,nk2cmsel,s,kk23ynk2=ky(nk2min1)da32=abs(ynk3-ynk2)*if,da32,le,da,thenlstr,nk3min,nk2min1*endifcmsel,s,kk2nk2min1=kpnext(nk2min1)*enddocmsel,s,kk3nk3min=kpnext(nk3min)*enddo!(创建其它径向切割线)lsel,s,loc,x,r3+delta,r2-deltaksllcm,kk4,kpcmsel,s,kk2cmsel,u,kk4*get,nk2,kp,count*get,nk2min,kp,num,mincm,kk2,kpcmsel,s,kk3cmsel,u,kk4*get,nk3,kp,count*get,nk3min,kp,num,mincm,kk3,kpcsys,1dsys,1allselnumcmp,kpnumcmp,linenumcmp,area*get,nkmax,kp,num,maxcmsel,s,kk2*do,i,1,nk2ynk2=ky(nk2min)k,nkmax+i,r3-delta,ynk2,0lstr,nkmax+i,nk2mincmsel,s,kk2nk2min=kpnext(nk2min)*enddoallselnumcmp,kpnumcmp,linenumcmp,area*get,nkmax,kp,num,maxcmsel,s,kk3*do,i,1,nk3ynk3=ky(nk3min)k,nkmax+i,r2+delta,ynk3,0lstr,nkmax+i,nk3mincmsel,s,kk3nk3min=kpnext(nk3min)*enddo!(完成切割)lsel,s,loc,x,r2-delta,r2+deltalsel,a,loc,x,r3-delta,r3+deltaasll,s,1lsel,s,loc,x,r3+delta,r2-deltaasbl,all,all!创建转子轴向孔csys,1dsys,0xo=rholeyo=360/(2*nslotr)zo=0wpave,xo,yo,zocyl4,rh,0,90cyl4,rh,90,180cyl4,rh,180,270cyl4,rh,270,360ksel,s,loc,x,xoksel,r,loc,y,yolslkasllallsel,below,areanummrg,kp!(如上复杂操作的目的是为了后面有良好的网格划分)cm,aa1,areaagen,2,allcmsel,u,aa1aadd,allagen,nhole,all,360/nholecm,aa2,arealsel,s,loc,x,r4asllcm,aa3,areacmsel,a,aa2allsel,below,areaasba,aa3,aa2cm,aa3,areacmsel,s,aa1agen,nhole,all,360/nholecmsel,a,aa3allsel,below,areanummrg,kpsave!到此为止，几何建模完成!下面设置网格划分控制!先将面形成组件csys,1lsel,s,loc,x,r4asll,s,1csys,0wpcsys,-1,0wprota,90asbw,allwprota,90asbw,allwpcsys,-1,0cm,a_axis,area !转子轴csys,1lsel,s,loc,x,r4asllcmsel,u,a_axiscm,a_rotor,area !转子铁芯rmin=rhole-rh+deltarmax=rhole+rh-deltaasel,s,loc,x,rmin,rmaxcm,a_hole,area !转子轴向孔asel,s,loc,x,r3-sr_h1-sr_h2-sr_h3,r3-sr_h1cm,a_rc,area !转子槽中的导体asel,s,loc,x,r3-sr_h1,r3-deltacm,a_ri,area !转子槽中的气隙asel,s,loc,x,r3+delta,r2-deltacm,a_gap,area !气隙lsel,s,loc,x,r1asllcm,a_stator,area !定子铁芯lslaasllcmsel,u,a_gapcmsel,u,a_statorrmax=r2+ss_h-ss_ch-ss_cirmin=rmax-ss_ciasel,a,loc,x,rmin,rmaxcm,a_si,area !定子槽中的绝缘和空气rmax=r2+ss_h-ss_cirmin=rmax-ss_ci-2*ss_chasel,s,loc,x,rmin,rmaxcmsel,u,a_sicm,a_sc,area !定子槽中的导体!定义各面的单元属性:材料,实常数,单元类别cmsel,s,a_axiscmsel,a,a_holecmsel,a,a_ricmsel,a,a_gapcmsel,a,a_siaatt,1,1,1cmsel,s,a_scaatt,2,2,2cmsel,s,a_rcaatt,3,3,3cmsel,s,a_rotorcmsel,a,a_statoraatt,4,4,4!定义网格尺寸esize_gap=3 !气隙径向网格划分3层esize_r_gap=2 !转子槽气隙周向划分2层esize_min=ss_cw/1.8 !除气隙及转子导条外的最小网格尺寸esize_min1=esize_min/2 !转子导条网格尺寸（考虑到集肤效应，网格尺寸尽量小）esize_mid=5*esize_min !中等尺度的网格尺寸esize_max=15*esize_min !最大的网格尺寸mshkey,2 !尽量用映射网格划分!(划分气隙网格)cmsel,s,a_rilslalsel,r,loc,x,r3-delta,r3+deltacm,ll1,linelesize,all,esize_r_gap,1,1,nolsel,s,loc,x,r3-delta,r3+deltacmsel,u,ll1lesize,all,esize_min,1,1,nolsel,s,loc,x,r3+delta,r2-deltalesize,all,esize_gap,1,1,nocmsel,s,a_gaplslaamesh,all!(划分转子槽中的气隙的网格)rmin=r3-sr_h1lsel,s,loc,x,rmin-delta,rmin+deltacm,ll2,linecmsel,s,a_rilslacmsel,u,ll1cmsel,u,ll2lesize,all,esize_min,1,1,nocmsel,ll2lesize,all,esize_r_gap,1,1,nocmsel,s,a_rilslaamesh,all!(划分定子槽中导体的网格)cmsel,s,a_sclslalesize,all,esize_min,1,1,noamesh,all!(划分转子槽中导体的网格)cmsel,s,a_rclslalesize,all,esize_min1,1,1,noamesh,all!(划分定子槽中的绝缘和空气的网格)cmsel,s,a_silslalesize,all,esize_min,1,1,noamesh,all!(划分转子轴向孔处的网格)cmsel,s,a_holelslalesize,all,esize_min,1,1,noamesh,all!(划分转子轴的网格)cmsel,s,a_axislslalesize,all,esize_mid,1,1,noamesh,all!(划分转子铁芯的网格)cmsel,s,a_rotoramesh,all!(划分定子铁芯的网格)lsel,s,loc,x,r1lesize,all,esize_max,1,1,nocmsel,s,a_statoramesh,allallselsave!至此,有限元部分的模型建立完毕! (第二段程序) !本程序的目的: 定义电路模型及其与有限元模型的连接!finish/prep7local,11,1csys,11cscir,11,1et,11,124,6 !电路耦合块导体(用于转子导条与电路的耦合)r,11et,12,124,0 !电阻,用于模拟转子端环阻抗r,12,1e-10et,13,124,1 !电感,用于模拟转子端环漏抗r,13,1e-12et,21,124,5 !电路耦合绞线圈,用于外层线圈(半径较大者)r,21et,22,124,5 !电路耦合绞线圈,用于内层线圈(靠近气隙者)r,22et,31,124,1 !电感(用于定子端部连接,模拟每个线圈的端部电感值)r,31,1e-4,0 !(暂将电感设为极小值,初始电流定义为零)et,32,124,0 !电阻(用于定子端部连接,模拟每个线圈的端部电阻值)r,32,0.014 !电阻值设定与线圈内部直线段电阻值接近(计算可得每个线圈直线段电阻为0.01395)et,33,124,0 !电阻(用于模拟导线连接,其电阻值设为一极小值)r,33,1e-10,0et,41,124,4,0 !独立电压源(用于建立三相电源)r,41,6000,0 !A相,电压6000V,相位为0r,42,6000,120 !B相,电压6000V,相位为120r,43,6000,-120 !C相,电压6000V,相位为-120!首先定义定子线圈(电路耦合绞线圈单元)和转子导条(电路耦合2D块导体)的自由度耦合rmax=r2+ss_h-ss_cirmin=rmax-ss_chasel,s,loc,x,rmin,rmaxcmsel,r,a_sccm,a_sc1,areacmsel,s,a_sccmsel,u,a_sc1cm,a_sc2,area!(定子线圈自由度耦合)da_s=360/nslots*dim,aa_sc1,array,nslots !获取面的编号,以备后用*dim,aa_sc2,array,nslots*do,i,1,nslotsang1=(i-1)*da_scmsel,s,a_sc1asel,r,loc,y,ang1*get,aa,area,num,minaa_sc1(i)=aa !外层(半径较大者)线圈面nsla,1cp,next,curr,allcp,next,emf,allcmsel,s,a_sc2asel,r,loc,y,ang1*get,aa,area,num,minaa_sc2(i)=aa !内层(靠近气隙者)线圈面nsla,1cp,next,curr,all !耦合CURR自由度cp,next,emf,all !耦合EMF自由度*enddo!(转子导条自由度耦合)da_s=360/nslotr*dim,aa_rc,array,nslotr !获取面的编号,以备后用*do,i,1,nslotrang1=(i-1)*da_scmsel,s,a_rcasel,r,loc,y,ang1*get,aa,area,num,minaa_rc(i)=aa !转子导条面nsla,1cp,next,curr,allcp,next,emf,all*enddo!下面定义转子导条与电路的耦合 - 二维电路耦合块导体allselnumcmp,node*get,nmax,node,num,max*dim,nn_rc1,array,nslotr*dim,nn_rc2,array,nslotrtype,11real,11zz=r3/2*do,i,1,nslotrasel,s,area,aa_rc(i)nsla,1*get,nn3,node,num,minnmax=nmax+1nn_rc1(i)=nmaxn,nmax,nx(nn3),ny(nn3),zznmax=nmax+1nn_rc2(i)=nmaxn,nmax,nx(nn3),ny(nn3),-zze,nn_rc1(i),nn_rc2(i),nn3 !建立电路耦合块导体单元*enddo*get,nmax,node,num,maxxx=nx(nn_rc1(1)yy=ny(nn_rc1(1)dyy=360/nslotryy=yy-dyy/2nmax=nmax+1nn1=nmaxn,nn1,xx,yy,zznmax=nmax+1nn2=nmaxn,nn2,xx,yy,-zztype,12real,12e,nn_rc1(nslotr),nn1e,nn_rc2(nslotr),nn2type,13real,13e,nn1,nn_rc1(1)e,nn2,nn_rc2(1)*do,i,1,nslotr-1yy=yy+dyynmax=nmax+1nn1=nmaxn,nn1,xx,yy,zznmax=nmax+1nn2=nmaxn,nn2,xx,yy,-zztype,12real,12e,nn_rc1(i),nn1 !创建端环电阻e,nn_rc2(i),nn2type,13real,13e,nn1,nn_rc1(i+1) !创建端环漏抗e,nn2,nn_rc2(i+1)*enddo!下面定义定子绕组与电路的耦合 - 电路耦合绞线圈!(先定义绞线圈)allselnumcmp,node*get,nmax,node,num,max*dim,nn_sc11,array,nslots !外层线圈上节点*dim,nn_sc12,array,nslots !外层线圈下节点*dim,nn_sc21,array,nslots !内层线圈上节点*dim,nn_sc22,array,nslots !内层线圈下节点type,21real,21zz=r1/2*do,i,1,nslotsasel,s,area,aa_sc1(i)nsla,1*get,nn3,node,num,minnmax=nmax+1nn_sc11(i)=nmaxn,nmax,nx(nn3),ny(nn3),zznmax=nmax+1nn_sc12(i)=nmaxn,nmax,nx(nn3),ny(nn3),-zze,nn_sc11(i),nn_sc12(i),nn3 !建立外层电路耦合绞线圈单元*end