自己建的一个2D场路耦合的例子,模型建立没有什么问题,求解时总出现错误
Error:Element 1 has negative or zero cross section specified.check real constant table11
请高手指教,是不是我的求解程序有问题?
下面是我的程序
/PREP7 !开始前处理程序
/UNITS, SI !选取国际单位
ET,1,53,0 ! AIR
ET,2,53,3 ! STATOR COIL(CIRCUIT-COUPLED STRANDED COIL)
ET,3,53,0 ! 铁心和衔铁
EMUNIT, MKS !选取(米、千克、秒)电磁单位制
MP,MURX,1,1 !AIR
MP,MURX,2,1 !STATOR WINDING
mp,rsvx,2,0.0434e-6 !导体电阻率
mp,murx,3,2000
LOCAL,11,0 !选取直角坐标系统
csys,11
K, 1,-0.05, 0.000
K, 2,-0.05, 0.100
K, 3, 0.05, 0.100
K, 4, 0.05, 0.000
K, 5,-0.02, 0.010
K, 6,-0.02, 0.030
K, 7, 0.02, 0.030
K, 8, 0.02, 0.010
K, 9,-0.02, 0.050
K, 10,-0.02, 0.070
K, 11, 0.02, 0.070
K, 12, 0.02, 0.050
K, 13,-0.04, 0.030
K, 14,-0.04, 0. 075
K, 15,-0.02, 0.075
K, 16, 0.02, 0.075
K, 17, 0.04, 0.075
K, 18, 0.04, 0.030
K, 19,-0.04, 0.080
K, 20,-0.04, 0.090
K, 21, 0.04, 0.090
K, 22, 0.04, 0.080
LSTR, 1, 2 ! L1 第一条线(Straight line)
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 1
LSTR, 5, 6 ! L5
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 5
LSTR, 9, 10
LSTR, 10, 11 ! L10
LSTR, 11, 12
LSTR, 12, 9
LSTR, 6, 13
LSTR, 13, 14
LSTR, 14, 15 ! L15
LSTR, 15, 10
LSTR, 11, 16
LSTR, 16, 17
LSTR, 17, 18
LSTR, 18, 7 ! L20
LSTR, 19, 20
LSTR, 20, 21
LSTR, 21, 22
LSTR, 22, 19
LSTR, 14, 19 ! L25
LSTR, 15, 16
LSTR, 17, 22 ! L27
LESIZE,1,,,20 !求解区边界,分得粗一些
LESIZE,2,,,20
LESIZE,3,,,20
LESIZE,4,,,20
LESIZE,5,,,5 !绕组下部,分得稍密一些
LESIZE,6,,,10
LESIZE,7,,,5
LESIZE,8,,,10
LESIZE,9,,,5 !绕组上部,分得稍密一些
LESIZE,10,,,20 !气隙部分要分得最密
LESIZE,11,,,5
LESIZE,12,,,10
LESIZE,13,,,5 !铁心
LESIZE,14,,,10
LESIZE,15,,,10 !气隙部分要分得最密
LESIZE,16,,,3
LESIZE,17,,,3 ! 与L16相同
LESIZE,18,,,10 ! 与L15相同
LESIZE,19,,,10 ! 与L14相同
LESIZE,20,,,5 ! 与L13相同
LESIZE,21,,,4 ! 衔铁部分要分得密
LESIZE,22,,,30
LESIZE,23,,,4
LESIZE,24,,,40
LESIZE,25,,,3
LESIZE,26,,,20
LESIZE,27,,,3
AL, 5, 6, 7, 8 ! A1 绕组下部
AL, 9, 10, 11, 12 ! A2 绕组上部
LSEL,S,,,13,16,1,1 ! 选一组线 13,14,15,16
LSEL,A,,,17,20,1 ! 加选17,18,19,20
LSEL,A,,,6,12,3 ! 加选6,9,12
LSEL,A,,,11 ! 加选11
AL,ALL ! 由所选的所有线段组成面A3
! 产生衔铁和气隙部分面积
ALLSEL ! 选取所有点线面
AL,21,22,23,24 ! A4 衔铁
AL,15,25,24,27,18,26 ! A5 主气隙
AL,10,16,26,17 ! A6 线圈与衔铁间气隙
AL,1,2,3,4 ! A7 4条外边界线构成
AOVLAP,ALL ! 所有面重叠后产生一个新面A8
! A7 面不再存在
allsel
NUMCMP,ALL
asel,s,,,5,7,1,1
AATT,1,,1
asel,s,,,1,,,1
AATT,2,11,2
asum
*get,aa,area,,area
r,11,aa,150,,1,0.8 !下线圈单元的实常数
asel,s,,,2,,,1
AATT,2,11,2
asum
*get,aa,area,,area
r,13,aa,150,,-1,0.8 !上线圈单元的实常数
asel,s,,,3,4,1,1
AATT,3,,3
ALLSEL
AMESH,ALL
finish
/prep7
local,11,0
csys,11
et,10,124,4,1 !独立电压源(用于建立单相电源)
r,10,0,115,60 !电压115V,相位为0,f=60Hz
et,11,124,5 !电路耦合绞线圈,用于线圈
r,11
r,13
et,12,124,0 !电阻,用于端环阻抗
r,12,1e-10
!(线圈自由度耦合)
asel,s,,,1,,,1
nsla,,1
cp,1,curr,all
cp,2,emf,all
asel,s,,,2,,,1
nsla,,1
cp,3,curr,all
cp,4,emf,all
allsel
numcmp,node
*get,nmax,node,,num,max
nmax1=nmax+1
nmax2=nmax+2
nmax3=nmax+3
nmax4=nmax+4
nmax11=nmax+5
nmax22=nmax+6
nmax33=nmax+7
nmax44=nmax+8
nmax5=nmax+9
asel,s,,,1,,,1
nsla,,1
*get,nn1,node,,num,min
xx=nx(nn1)
yy=ny(nn1)
n,nmax1,xx,yy,0.05
n,nmax11,xx,yy,0.10
n,nmax2,xx,yy,-0.05
n,nmax22,xx,yy,-0.10
asel,s,,,2,,,1
nsla,,1
*get,nn2,node,,num,min
xx=nx(nn2)
yy=ny(nn2)
n,nmax3,xx,yy,0.05
n,nmax33,xx,yy,0.10
n,nmax4,xx,yy,-0.05
n,nmax44,xx,yy,-0.10
xx=nx(nmax22)
yy1=ny(nmax22)
yy2=ny(nmax44)
yy=0.5*(yy1+yy2)
n,nmax5,xx,yy,0.10
allsel
type,11
real,11
e,nmax2,nmax1,nn1
real,13
e,nmax3,nmax4,nn2
type,12
real,12
e,nmax1,nmax11
e,nmax2,nmax22
e,nmax3,nmax33
e,nmax4,nmax44
e,nmax22,nmax44
type,10
real,10
e,nmax33,nmax11,nmax5
allsel
save
finish
/solu
antype,trans
trnopt,full
allsel
tm_start=1e-5
tm_end=30
tm_incr=1.5
*do,tm,tm_start,tm_end,tm_incr
time,tm
csys,11
LSEL,S,,,1,4,1,1 !选择L1,L2,L3,L4
NSLL,S !选择线上的所有节点
D,ALL,AZ,0 !对边界节点赋值AZ=0
allsel
d,nmax11,volt,0 !接地
allsel
neqit,50
cnvtol,a,,1e-5
solve
finish |