电机绕组反电势为什么会这么大 附求反电势代码

cqlyj2008

高手帮我看看我做的同步发电机磁场az图以及后处理定子感应相电压，


为什么峰值会有16000那么大（实际8980左右）。我的铁芯场0.98m，根据az图大致可以算出每极磁通0.46左右，那么定子电压有效值就e=-4.44fφnk，比求出来的小多了。这是为什么了，期待高手为我解答啊，百般感谢。以下附上求反电势源代码：
*dim,az1,array,z1
*dim,az2,array,z1
*dim,elem,array,z1
*dim,elemarea,array,z1
*dim,Fi0,array,z1
*dim,Fi1,array,z1
*dim,Fi2,array,z1
*dim,E,table,z1+1,2
*dim,E1,array,z1
*dim,E2,array,z1
E(0,0)=1e-5
E(0,1)=1
E(0,2)=2
/POST1
rsys,1
*do,ii,1,3*q1,1
    asel,s,area,,aa_sc1(ii)                                         !线圈上层边
    esla,s
    *get,elemnum,elem,0,count
    *get,elemmin,elem,0,num,min
    *do,k,1,elemnum,1
        elem(k)=elemmin
     esel,s,,,elem(k)
    *get,elema,elem,elem(k),area
     elemarea(k)=elema
     nsle,s,all
    *get,nodenum,node,0,count
    *get,nodemax,node,0,num,max
    *get,nodeazm,node,nodemax,A,z
    sumaz=nodeazm
    *do,jj,1,nodenum-1,1
        nodenext=nodemax        
        *get,nodenext,node,nodemax,nxtl
        *get,nodeaz,node,nodenext,A,z
        sumaz=sumaz+nodeaz   
    *enddo
    az1(k)=sumaz/nodenum                                      !求线圈边单元平均磁势
    Fi0(k)=az1(k)*elemarea(k)
    elemmin=elemmin+1
    *enddo
   Fi1(ii)=0
   *do,k,1,elemnum,1
      Fi1(ii)=Fi1(ii)+Fi0(k)                                   !求线圈边平均磁势
   *enddo
   
   star=0
   *if,ii+y1,gt,z1/4,then                                      !因为只建了1/4之一模型
       star=1
   *endif
    asel,s,area,,aa_sc2(ii+y1-star*z1/4)                       !线圈下层边
    esla,s
    *get,elemnum,elem,0,count
    *get,elemmin,elem,0,num,min
    *do,k,1,elemnum,1
        elem(k)=elemmin
     esel,s,,,elem(k)
    *get,elema,elem,elem(k),area
     elemarea(k)=elema
     nsle,s,all
    *get,nodenum,node,0,count
    *get,nodemax,node,0,num,max
    *get,nodeazm,node,nodemax,A,z
    sumaz=nodeazm
    *do,jj,1,nodenum-1,1
        nodenext=nodemax        
        *get,nodenext,node,nodemax,nxtl
        *get,nodeaz,node,nodenext,A,z
        sumaz=sumaz+nodeaz   
    *enddo
    az2(k)=sumaz/nodenum                                      !求线圈边一个单元平均磁势
    Fi0(k)=az2(k)*elemarea(k)
    elemmin=elemmin+1
    *enddo
   Fi2(ii)=0                        
   *do,k,1,elemnum,1
      Fi2(ii)=Fi2(ii)+Fi0(k)                                   !求线圈边平均磁势           
   *enddo
*enddo
!====求取A相线圈上下层边的磁链和====
*do,ii,1,2*p,1
*if,mod(ii,2),eq,0,then
  *do,i,1,3*q1,1
   *if,i+q1,gt,z1/4,then
       flag1=-1
       flag=1
   *else
       flag1=1
       flag=0
   *endif
       E1(i+(ii-1)*3*q1)=-lt*sn1/2/carea1*(flag1*Fi1(i+q1-flag*z1/4)-Fi1(i))     !上层边磁链变化量
   *if,i+q1,gt,z1/4+3,then
      flag1=-1
   *else
      flag1=1
   *endif
   *if,i,gt,3,then
       E2(i+(ii-1)*3*q1)=-lt*sn1/2/carea1*(-flag1*Fi2(i+q1-flag*z1/4)+Fi2(i))    !下层边磁链变化量
     *else
       E2(i+(ii-1)*3*q1)=-lt*sn1/2/carea1*(-flag1*Fi2(i+q1-flag*z1/4)-Fi2(i))   
     *endif
  *enddo
*else
  *do,i,1,3*q1,1
   *if,i+q1,gt,z1/4,then
       flag1=-1
       flag=1
   *else
       flag1=1
       flag=0
   *endif
       E1(i+(ii-1)*3*q1)=lt*sn1/2/carea1*(flag1*Fi1(i+q1-flag*z1/4)-Fi1(i))     !上层边磁链变化量
   *if,i+q1,gt,z1/4+3,then
      flag1=-1
   *else
      flag1=1
   *endif
   *if,i,gt,3,then
       E2(i+(ii-1)*3*q1)=lt*sn1/2/carea1*(-flag1*Fi2(i+q1-flag*z1/4)+Fi2(i))    !下层边磁链变化量
     *else
       E2(i+(ii-1)*3*q1)=lt*sn1/2/carea1*(-flag1*Fi2(i+q1-flag*z1/4)-Fi2(i))   
     *endif
  *enddo
*endif
*enddo
!===用于求旋转整个槽距角时的绕组电压值====
*do,i,1,z1+1,1
   E(i,1)=p*(i-1)*360/z1
   *if,i,gt,z1,then
     flag=1
   *else
     flag=0
   *endif
   E(i,2)=-2*z1*f*(E1(i-flag*z1)-E2(i-flag*z1))!因为每次旋转一个槽距角，所以时间1/z1/f,再乘以并联支路数a=2                    
*enddo
!====绘制电压波形图====
/grid,1
/axlab,x,E-agnle (o)
/axlab,y,phase A voltage (V)
*vscf,xmin,min,E(1,1)
*vscf,xmax,max,E(1,1)
/xrange,xmin,xmax
*vscf,ymin,min,E(1,2)
*vscf,ymax,max,E(1,2)
/yrange,(nint(ymin/1000)-5)*1000,(nint(ymax/1000)+5)*1000
/color,curve,orange,1
*vplot,E(1,1),E(1,2)

fossiler

你确定并联支路在计算反电势的时候是要乘以并联支路数。

cqlyj2008

2# fossiler
因为我每次只算了一个极相组线圈的磁链，所以得乘以2P，因为并联支路数是2，所以最后该乘以2，是这样吧？望您指正啊

fossiler

磁链=磁通*匝数
1匝导线不管是由多少并联支路组成，都应该按照1匝来计算吧。就像两个电压源并联以后你不能说他的外特性就是二倍的电动势吧。
那么你就要确定你的两个支路是并联还是串联，另外你在加载支路电流的时候是否加载正确，是以支路电流加的还是以相电流加的。

cqlyj2008

4# fossiler
我算的是4极同步发电机，在转子线圈上加的电流密度，现在求的是定子相电压。定子绕组是两路并联，我上面算的磁链只是一个极相组的，也就是6个线圈的磁链，实际上应该算24个线圈的磁链再除以并联支路数，那么也就是12个线圈的磁链了。所以上面求出的6个线圈的磁链乘以2，应该是这样吧。fossiler 兄，望您批评指正。

pengchangjun

e=-4.44fφnk
φ为基波磁通，气隙磁密傅立叶分解得到，
这公式计算完全正确，误差小于1%。

cqlyj2008

6# pengchangjun
这个公式不是有效值的计算公式吗？我要得到的是每时刻的电压波形 啊