flyash 发表于 2007-6-20 20:55:07

关于三维棱边单元法静态分析的转矩误差问题

我用的SOLID117，作 3－D静态分析，在如图中所得到的转矩理论上应该是正值，为什么反而得到的转矩为负值呢？而且，不同位置的转矩应该是正弦分布的，至少是近似正弦分布，可为什么得到的值总是不对呢？

我用的EMFT求转矩，求的是作用在线圈上的转矩～线圈上加载电流，带铁心
到底问题会出现在哪儿呢？程序查了好多遍，网格也细化过，最多的时候划了10万多个点，可结果还是有问题
下面的是永磁体，上面的是线圈，中间有气隙
请指教，谢谢

程序如下：
FINISH
/CLEAR
/PREP7
/PNUM,VOLUME,1
EMUNIT,MKS


INNERRADIUS=0.06
PMTHICKNESS=0.04
THETA=22.5
T=-8
FAI=30
AIRGAP=0.001      
A=0.0054               
B=0.0104               
HEIGHTCOIL=0.041275      

HC=930000               
N=1500                  
I=1                     

ACOIL=(B-A)*HEIGHTCOIL   
JDENS=N*I/ACOIL      




*AFUN,DEG
LOCAL,11,0,,,,,-90
WPCSYS,11
CYL4,0,0,INNERRADIUS,-FAI/2,INNERRADIUS+PMTHICKNESS,FAI/2
CSYS,0
WPCSYS,0
K,1000,0,0,-1
K,1001,0,0,1
VROTAT,1,,,,,,1000,1001,90,1
KDELE,1000,1001,1
CSYS,11
WPCSYS,11
WPROTA,0,0,-THETA
RECTNG,-0.15,0.15,-0.15,0.15
WPROTA,0,0,-THETA
RECTNG,-0.15,0.15,-0.15,0.15
WPROTA,0,0,-THETA
RECTNG,-0.15,0.15,-0.15,0.15
ASEL,S,AREA,,7,9,1,1
VSBA,ALL,ALL,,DELETE,DELETE
NUMCMP,VOLU
WPROTA,0,0,3*THETA
CSYS,0
WPCSYS,0



WPROTA,2*THETA+T
WPOFFS,INNERRADIUS+PMTHICKNESS+AIRGAP
WPROTA,,,90
CYL4,0,0,0,0,A,360,HEIGHTCOIL
CYL4,0,0,A,0,B,360,HEIGHTCOIL
VGLUE,5,6
WPAVE,0,0,0
WPROTA,0,0,-90
WPROTA,-2*THETA-T




CSYS,11
WPCSYS,11
CYL4,0,0,0,-FAI/2-15,0.3,FAI/2+15
CSYS,0
WPCSYS,0
K,1000,0,0,-1
K,1001,0,0,1
VROTAT,7,,,,,,1000,1001,90,1
KDELE,1000,1001,1
VSEL,ALL
VOVLAP,ALL
NUMCMP,VOLU


ET,1,SOLID117

MP,MURX,1,1
MP,MURX,2,1
MP,MGXX,2,HC
MP,MURX,3,1
MP,MGXX,3,-HC
MP,MURX,4,1
MP,MURX,5,1000

LOCAL,12,0,,,,THETA/2
VSEL,S,VOLU,,1
VATT,2,0,1,12

LOCAL,13,0,,,,3/2*THETA
VSEL,S,VOLU,,3
VATT,3,0,1,13

LOCAL,14,0,,,,5/2*THETA
VSEL,S,VOLU,,4
VATT,2,0,1,14

LOCAL,15,0,,,,7/2*THETA
VSEL,S,VOLU,,2
VATT,3,0,1,15

CSYS,2
VSEL,S,VOLU,,7
VATT,1,0,1

VSEL,S,VOLU,,6
VATT,4,0,1

VSEL,S,VOLU,,5
VATT,5,0,1


CSYS,2
/PNUM,MAT,1
VPLOT

ALLSEL,ALL
smrtsize,3
MSHAPE,1,3D
VMESH,ALL
EPLOT


!施加电流密度
LOCAL,16,1,(INNERRADIUS+PMTHICKNESS+AIRGAP)*COS(2*THETA+T),(INNERRADIUS+PMTHICKNESS+AIRGAP)*SIN(2*THETA+T),0,2*THETA+T,0,90
VSEL,S,VOLU,,6
ESLV
EMOD,ALL,ESYS,16      
BFE,ALL,JS,,,JDENS


CSYS,2
ALLSEL
FINISH


/SOLU
ANTYPE,STATIC
!边界条件
CSYS,2
ASEL,S,LOC,X,0.3
NSLA               
D,ALL,AZ,0
ASEL,S,AREA,,30
ASEL,A,AREA,,31
NSLA
D,ALL,AZ,0


ALLSEL,ALL
MAGSOLV
SAVE
FINISH

/POST1
ESEL,S,MAT,,4
NSLE
ESEL,ALL
RSYS,2
EMFT

ljzhit 发表于 2007-6-21 09:18:14

我做永磁电机分析求转矩时,求的也不对,也找不到原因.不知为什么?有哪位大虾能说说求转矩的方法吗？

grandjiang 发表于 2007-6-28 09:26:22

我也是做电机分析的,同样计算出的结果是在不同位置的受力有正有负,百思不得其解,
我的QQ是279255056,望有见解或在做3D电机分析的朋友们给点帮助

hbbdyzf 发表于 2015-9-28 10:49:26

我也出现同样问题
页: [1]
查看完整版本: 关于三维棱边单元法静态分析的转矩误差问题