ansys中如何算损耗
众所周知,ansys有静态,谐波,瞬态三磁场种场求解方式。Bertotti提出过常系数三相式铁耗计算模型: Pfe=Kh*f*B^2+Ke*f^2*B^2+Ka*f^1。 5*B^1。5
式中,第1项为磁滞损耗,第二项为涡流损耗,第三项为巴克豪森越变产生的损耗,即附加损耗。其中Kh
,
Ke,
Ka均为相应的损耗系数;Ke有想应的计算公式,Kh
,Ka可通过硅钢片实测损耗数据拟合求得。 ansoft就是用这个算的。
问题1:硅钢片电阻率的设置。定子采用叠片铁心时,片间垫有绝缘纸,那么应该设置x,y,z三个方向的电阻率才对。...全部
众所周知,ansys有静态,谐波,瞬态三磁场种场求解方式。Bertotti提出过常系数三相式铁耗计算模型: Pfe=Kh*f*B^2+Ke*f^2*B^2+Ka*f^1。
5*B^1。5
式中,第1项为磁滞损耗,第二项为涡流损耗,第三项为巴克豪森越变产生的损耗,即附加损耗。其中Kh
,
Ke,
Ka均为相应的损耗系数;Ke有想应的计算公式,Kh
,Ka可通过硅钢片实测损耗数据拟合求得。
ansoft就是用这个算的。
问题1:硅钢片电阻率的设置。定子采用叠片铁心时,片间垫有绝缘纸,那么应该设置x,y,z三个方向的电阻率才对。z方向这值多大呢?
问题2:感应电机用谐波场计算,通过powerh宏可以算的单位时间内的功率损失,但是这个损失考虑了磁滞效应没(应该没有吧)。
问题3:永磁电机因含有永磁体,不能用谐波场。可以用面电流来模拟永磁体。这里就涉及到瞬态求解损耗的问题。用时步有限元计算一个周期内的磁场分布,在/post1中可以用etable取出eddy(涡流密度),也可以取出jheat(焦耳热),不知道这两者之间的关系是,未验证??
那么在每一步我们都可以通过公式p=jt*jt*电阻率*体积算出涡流损耗。
那么在一个周期内就可以得到n个时刻的涡流损耗p,那么电机的涡流损耗是不是所有的p加起来除以T(周期)呢?
问题4:我们是不是也可以像ansoft那也,算完一个周期后,提取单元磁密,傅里叶分解。
输入损耗曲线,然后分别计算三种损耗,再加起来。但这样实现起来好像比较复杂,计算量也大,做过的高手请指点指点。
欢迎大家讨论!!!
以下是本人用ansys计算铁耗的一点思路,大概是用时步有限元算完一个周期的磁场后,通过post26提取个单元磁密波形,再进行相关处理。
其中几个log还未完成,欢迎大家发表自己的意见!!期待和您讨论
/post26
t=0。5 !硅钢片厚度
ro=9。
8*1e3 !硅钢片密度
gama=2。33*1e7 !硅钢片电导率
ke=pi**2*gama*t**2/6/ro !涡流损耗系数
/input,plcurve,log !硅钢片的损耗曲线
/input,fitting,log !拟合曲线,求取损耗系数kh,ka
allsel
csys,1
esel,s,mat,,4 !选取定子硅钢片
*get,enum,elem,0,count
*get,emin,elem,0,num,min
*dim,Bx,array,360,1
*dim,Ploss,array,10
*dim,Plosum,array,enum
*do,i,1,enum
*get,volu,elem,emin,volu !提取单元体积
esol,2,emin,,b,x !获取径向磁密波形
vget,Bx(1,1),2,0,0 !将径向B值存入数组中
/input,definition,log !定义数组
/input,fourier,log !傅里叶分解磁密
Bj=fundamental_coeff !基波幅值
Pf=kh*f*Bj**2+ke*f*Bj**2+ka*f*Bj**1。
5 !基波损耗
Pf=Pf*volu*ro*t !单元基波损耗
plosum(i)=Pf
*do,ii,1,10
Bm(ii)=top10(ii) !提取10次最大幅值的谐波,暂定10次
F(ii)=top10_mode(ii)
Ploss(ii)=kh*F(ii)*f*Bm(ii)**2+ke*(F(ii)*f*Bm(ii))**2+ka*(F(ii)*f*Bm(ii))**1。
5 !常系数三项式损耗模型
Ploss(ii)=Ploss(ii)*volu*ro*t
plosum(i)=Ploss(ii)+plosum(i) !将各次谐波损耗叠加
*enddo
*get,emin,elem,emin,nxth
*enddo
*do,i,1,enum !总损耗
Plsum=plsum+plosum(i)
*enddo。
收起