变压器求电感

/title,TransformerAnalysis/nopr/comnominal/commagneticfieldinthecore:Hn=(N1I1+N2I2)/y=25/comfluxdensityinthecore:Bn=Hn/nui=0.5/comabsoluteenergyinthecore:nlene=HnBn/22xyz=0.0125/comabsolutecoenergyinthecore:nlcene=BnHn2xyz-nlene=0.0125/comcoil1currentincrementedby=0.2A/commagneticfieldinthecore:H11=Hn+N1dI1/y=45/comfluxdensityinthecore:B11=Bn+(H11-Hn)/nui=0.9/comenergyincrementinhalfcore:dW11=(B11-Bn)(H11-Hn)/2*xyz=0.004/comcoil2currentincrementedby=0.2A/commagneticfieldinthecore:H22=Hn+N2dI2/y=65/comfluxdensityinthecore:B22=Bn+(Hn-H22)/nui=1.3/comenergyincrementinhalfcore:dW22=(B22-Bn)(H22-Hn)/2*xyz=0.016/combothcurrentsincrementedby=0.2A/commagneticfieldinthecore:H12=Hn+(N1dI1+N2dI2)/y=85/comfluxdensityinthecore:B12=Bn+(Hn-H12)/nui=1.7/comenergyincrementinhalfcore:dW12=(B12-Bn)(H12-Hn)/2*xyz=0.036/cominductances/comselfcoil1:L11~2N1^2xz/(ynui)=0.4/comselfcoil2:L22~2N2^2xz/(ynui)=1.6/commutualbetweencoil1and2:L12~2N1N2xz/(ynui)=0.8/comfluxlinkages/comcoil1:psi1=2N1xzB0=0.1/comcoil2:psi2=2N2xzB0=0.2/com!geometrydatan=1!meshingparameterx1=0.1!width(xsize)ofcoil1x2=0.1!width(xsize)ofcoil2x=0.1!width(xsize)ofcorey=0.1!heightofcore,ysizeofwindowz=0.1!thicknessofironinzdirectionnui=50!absolutereluctivityofironn1=10!numberofturnsincoil1n2=20!numberofturnsincoil2!excitationdatausedbyLMATRIX.MACsymfac=2!symmetricfactorforinductancecomputationnc=2!numberofcoils*dim,cur,array,nc!nominalcurrentsofcoils*dim,coils,char,nc!namesofcoilcomponentscur(1)=0.2!nominalcurrentof1stcoilcoils(1)='coil1'!nameofcoil1componentcur(2)=0.025!nominalcurrentof2ndcoilcoils(2)='coil2'!nameofcoil2component!derivedauxiliaryparametersmu0=3.1415926*4.0e-7muri=1/nui/mu0!relativepermeabilityofironx3=x1+x2!xcoordinateoftherightofcoil2x4=x3+x!xcoordinateofmiddleofcore(symmetryplane)/prep7mp,murx,1,1!coil1et,1,117mp,murx,2,1!coil2mp,rsvx,2,1et,2,117,5 !SolenoidalFormulation　对导电区用AZ-VOLT自由度，对不导电区用AZ自由度．!5 -- SolidConductors(DCCurrent)!AZ,VOLTdegreesoffreedom:nonlinearsymmetricsolenoidal!formulationapplicabletostaticandtransientanalyses.!6 -- SolidConductors(DCCurrent)!AZ,VOLTdegreesoffreedom:linearunsymmetricsolenoidal!formulationapplicabletoharmonicanalyses.mp,murx,3,muri!ironet,3,117block,0,x1,0,y,0,z!coil1block,x1,x3,0,y,0,z!coil2block,x3,x4,0,y,0,z!corevglue,allvsel,s,loc,x,x1/2vatt,1,1,1!coil1volumeattributevsel,s,loc,x,x1+x2/2vatt,2,2,2!coil2volumeattributevsel,s,loc,x,x3+x/2vatt,3,3,3!ironvolumeattributevsel,allesize,,nvmesh,allnsel,s,loc,x,x4!fluxparallelDirichletatsymmplain,x=x4,z=0,z=znsel,a,loc,z,0nsel,a,loc,z,zd,all,az,0!homogeneousNeumannfluxnormatyoke,x=0,y=0,y=ynsel,all!currentdensityloadoncoil1esel,s,type,,1bfe,all,js,,,,n1*cur(1)/(x1*y)cm,coils(1),elemallsnsel,s,loc,z,0d,all,volt,0nsel,s,loc,z,zcp,1,volt,all!建立耦合allsnn1amps=node(x1,0,z)f,nn1amps,amps,-n2*cur(2) !除了加电流密度载荷外，还可以给一个块导体加总电流F,,amps !注：在加总电流之前需耦合节点VOLT自由度。nsel,s,node,,nn1ampscm,coils(2),nodensel,allfini/com/comobtainoperatingsolution/com/solucnvtol,amps,1,1e-6cnvtol,csg,1,1e-6solvefini/post1presol,blmatrix,symfac,'coil','cur','ind'!computeinductancematrixpsi1=ind(1,3)psi2=ind(2,3)L11=ind(1,1)L22=ind(2,2)L12=ind(1,2)/com,********Computedsolution*********/com*vwrite,psi1('psi1=',f10.4)*vwrite,psi2('psi2=',f10.4)*vwrite,L11('L11=',f10.4)*vwrite,L22('L22=',f10.4)*vwrite,L12('L12=',f10.4)Finish注：如果et,2,117,5改为et,2,117，则对绕组2电流的加载可以和绕组1一样，为esel,s,type,,2bfe,all,js,,,,-n2*cur(2)/(x2*y)cm,coils(2),elemalls其计算结果是一样的