永磁电机设计程序1

#include "stdio.h" %头文件,用于调用标准输入输出函库数 #include "iostream.h" #include "math.h" %头文件,用于调用数学函数库中的函数 const double pi=3.1415926; %定义实数 int i,n; float x7[21], y7[21]; %定义实型数组 char in[15],out[15]; %定义字符型数组 float p2,p,m1,u,f,pfw,d1,di1,g,di2,l,x,sa,sx;%p2为输出功率，P为极数，m1为相数, U为相电压, F频率, PFW为风扇摩擦损耗, D1定子外径, DI1定子内径, G为气隙, DI2转子内径, L铁心长度, x为转子槽型（5种，1为平口槽，2为梯形槽，3为刀形槽，4为半凸型槽，5为凸型槽）, sa定子绕组类型（为1单层绕组，为2为双绕组）, sx为单层情况下绕组的三种类型（1为单层同心式，2为单层链式，3为单层交叉式）。 float q1,b01,hs0,bs1,bs2,h1,fys,hsp,alph;%Q1代表定子槽数, B01槽口宽, HS0槽口矩形高度, BS1为1/3处槽宽, BS2槽底宽, H1定子槽型尺寸, fys代表定子槽体（1为平底槽，2为圆底槽）, HSP槽楔厚度, ALPH定子槽间的角度 float z1,n1,d,n2,dd,a,lef,sef;%Z1定子线圈每圈匝数, N1每匝几个导体并联, D直径, N2每匝几个导体串联, DD导体直径, A为定子绕组并联支路数 float yy,q2,hr0,b02,fyr,bsk1,sr,lo,alph2;%q2转子槽数,HR0转子槽口矩形高度, b02转子槽口宽度, fyr是转子槽型（1平底，2为圆底）,BSK1代表转子有对斜槽（1表示斜一个定子槽的距离）, SR转子两个端环截面积, LO转子电阻系数（铜）, ALPH2转子槽楔角 float br1,br2,hr12,br3,br4,hr2,hr3,rr1,rr2;%BR1转子梯形槽上面宽度,BR2转子梯形槽下面宽度, br3,br4,HR12,hr3" 3.2, 3.2, 4.5,2.4,15, 20.6 float yd[10],akesai; float q,qq,rd1,rr,qp1,qp2,lfe,kfe=0.95,le,d2,bsk,ikw;% Bsk转子斜槽距 int iopen=1,mal=1; char al[]="zz-AL",cu[]="zz-Cu",bq[]="zz-BQ",d23[]="D23.dat"; float lob,lor,kb,alpha,alpha2,ld,hs,hs1,hs2,hr1,hrr,hrtotal,ty,dr,bet,bet1,ycg; float lbt1x,lbt1s,lbt1,lbt2sx,lbt2ss,v,lbt2xx,lbt2xs,lbt2x,lbt2s,lbt2,v1,r;% LBT1s= 定子齿上部的宽度 LBT1x定子齿下部的宽度，LBT2ss转子大齿上部宽度 LBT2sx转子大齿下部宽度 LBT2xs转子小齿上部宽度 LBT2xx转子小齿下部宽度 float sinrf,cs,ls,ks,lz,zfi,ci,kw1,kw2,se,se1,se2,se3,sf,sc,r1,r11,gcu; float gfe,xb,kd1,kp1,kdp1,v2,v3,xgms,k,sb,sb1,sb2,sb3,rbar,rring,r21,tp,t2; % Gfe铁的重量，kd1 绕组分布系数 kp1绕组短矩系数 kdp1绕组系数，Tp极距 ，t1定子齿距，t2转子齿距 ， float kc1,t1,kc2,ge,cx,xgmr,fy,oh,ob1,ob2,lu1,kl1,ku1,xe11,xs11,ll1,ol,ls1; float am,hc2p,sc1,hc1p,st1,sg,st2,st2x,st2s,kep,ike,ieta,ps1,pfw1,vc1,vt1,lc1,lc2,sc2; st1定子铁心一个极下齿部截面积，st2转子铁心一个极下齿部截面积，sc1定子轭部截面积，sc2转子轭部截面积，Ls1定子槽漏磁导= 1.5879， float ht1,hfe,att1,ft1,bfe,bt1,fg,atg,bg,fsp,fi,att2x,att2s,ht2s,ht2x,ft2,bt2,bt2x,bt2s,att2,ht2; ht1定子齿磁场强度，ht2转子齿磁场强度，hc1定子轭部磁场强度，hc2转子轭部磁场强度 ，Ft1定子齿磁位降，Ft2转子齿磁位降， float xsk2,xe21,xd21,xd11,im,im1,ak1,lj2p,lj1p,hc2,fc2,atc2,bc2,fc1,hc1,c1,c2,atc1,hcp,fs,ft,bc1; Fc1定子轭部磁位降，FC2转子轭部磁位降，c1定子轭部长度 c2转子轭部长度， float fu2,ls22,oll,i2p,lu2,oh3,ll2,ob3,ob4,oh2,oh1,ob0,pfe1,etap,pc1,bc10,pt1,bt10,km,kep0,x21,xx; float xs21,x11,att10,bg0,atg0,a1,do1,del1,i1,i11,eta,xgmap,xgmap1,p11,pcu21,pcu11,ir,i21,i2,ke,ir1; KE电势系数，I2起动计算 I21= 1.1175， float ix1,ip1,att2x0,att2s0,bt2x0,bt2s0,att20,bt20,im0,im10,hl0,atc20,bc20,atc10,hpr,obr,hpx,ohx; % Im0空载时励磁电流 float ohr,oh123,kst,tm,sn,cosfi,q0,delb,delr,a2,ll2st,hr4,rst1,krr,hr4p,bpr,sb3pr,sb31,abeta,hstr; % TM最大转矩，SN转差，COSFI功率因数，delb导条电流密度 ，delr端环电流密度， float bstr,sb2p,bprc,hprc,hpxc,kec,tst,ist,zst,xst,x2st,xskst,xd2st,xs2st,ls2st,i2t,lu2st,x1st,xd1st; % Ist起动电流倍数，TST起动转矩倍数，Rst启动总阻抗 xs2st转子槽漏抗 xd2st转子谐波漏抗 xskst转子斜槽漏抗 ZST启动总阻抗 X1ST定子侧启动电抗 X2ST转子侧启动电抗 XST启动总电抗 float xs1st,ls1st,cs1,cs2,kz,blp,bl,bl1,istp,br2p,hstx,bpxc,v11,b,kr1,kr11,v4,kr2,kec22,shkec,sikec,chkec,fikec; KZ启动漏磁路饱和系数，BL假想气隙磁密， float wo,wo1,vo,vo1,uo,uo1,ol0,ol1,ol11,ol12,ol13,ol21,ol22,ol23,olp,kkec,krf; float hl=0,xm1=0,bc2x=0,bc2x0=0,ol2,ol3,kr5p,kr4p,kr5,kr4,osb,s23,os3,os2,os1; FILE *fp1,*fp2,*fp3; double fnk1(double b2,double b1) { double aa; aa=0.009+0.324*b2/b1+0.053*(b2/b1-1)*(b2/b1-1); return (aa); } double fnk4(double b2,double b1) { double aa; aa=0.5+(b1/b2-1+log(b2/b1))/((b1 / b2 - 1)*(b1 / b2 - 1)); return (aa); } double fnk5(double b2,double b1) { double aa; aa=log(b2/b1)/(1-b1/b2); return (aa); } void sub3640(void) { for(i=0;i{ if(bfe<=x7[i+1]) { goto label3650; } } label3650: if(x==0) { goto label3680; } if(x>n-1) { i=n-2;goto label3680; } if(((x7[i]+x7[i+1])/2)>bfe) { i=i-1; } label3680: uo1=(bfe-x7[i+1])*(bfe-x7[i+2]); uo=uo1/(x7[i]-x7[i+1])/(x7[i]-x7[i+2]); vo1=(bfe-x7[i])*(bfe-x7[i+2]); vo=vo1/(x7[i+1]-x7[i])/(x7[i+1]-x7[i+2]); wo1=(bfe-x7[i])*(bfe-x7[i+1]); wo=wo1/(x7[i+2]-x7[i])/(x7[i+2]-x7[i+1]); hfe=vo*y7[i+1]+uo*y7[i]+wo*y7[i+2]; }

void sub3720(void) { if(p==2) { c1=0.86*hcp/tp-bfe/1.57+1.27; c2=0.86-0.4*bfe; } else { if(p==4) { c1=0.65*hcp/tp-bfe/1.57+1.25; c2=0.95-0.4*bfe; } else { if(p>=6) {